应用自然电场法研究地下水流场

讨论了不同水文地质条件下自然电场“8”字图特征,对评价输水管线等渗漏引起的环境问题具有指导意义。     

Salt transport on islands in the Okavango Delta: Numerical investigations

This study uses a numerical model to investigate the groundwater flow and salt transport mechanisms below islands in the Okavango Delta. Continuous evapotranspiration on the islands results in accumulation of solutes and the formation of a saline boundary layer, which may eventually become unstable. A novel Lagrangian method is employed in this study and compared to other numerical methods. The numerical results support the geophysical observations of density fingering on Thata Island. However, the process is slow and it takes some hundreds of years until density fingering is triggered. The results are sensitive to changes of the hydraulic gradient and the evapotranspiration rate. Small changes may lead to different plume developments. Results further demonstrate that density effects may be entirely overridden by lateral flow on islands embedded in a sufficiently high regional hydraulic gradient.     

A spatially and temporally adaptive solution of Richards’ equation

Efficient, robust simulation of groundwater flow in the unsaturated zone remains computationally expensive, especially for problems characterized by sharp fronts in both space and time. Standard approaches that employ uniform spatial and temporal discretizations for the numerical solution of these problems lead to inefficient and expensive simulations. In this work, we solve Richards’ equation using adaptive methods in both space and time. Spatial adaption is based upon a coarse grid solve and a gradient error indicator using a fixed-order approximation. Temporal adaption is accomplished using variable order, variable step size approximations based upon the backward difference formulas up to fifth order. Since the advantages of similar adaptive methods in time are now established, we evaluate our method by comparison with a uniform spatial discretization that is adaptive in time for four different one-dimensional test problems. The numerical results demonstrate that the proposed method provides a robust and efficient alternative to standard approaches for simulating variably saturated flow in one spatial dimension.     

Integrated surface–groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model

Interactions between surface and groundwater are a key component of the hydrologic budget on the watershed scale. Models that honor these interactions are commonly based on the conductance concept that presumes a distinct interface at the land surface, separating the surface from the subsurface domain. These types of models link the subsurface and surface domains via an exchange flux that depends upon the magnitude and direction of the hydraulic gradient across the interface and a proportionality constant (a measure of the hydraulic connectivity). Because experimental evidence of such a distinct interface is often lacking in field systems, there is a need for a more general coupled modeling approach.     

Influence of uncertain boundary conditions and model structure on flood inundation predictions

In this study, the GLUE methodology is applied to establish the sensitivity of flood inundation predictions to uncertainty of the upstream boundary condition and bridges within the modelled region. An understanding of such uncertainties is essential to improve flood forecasting and floodplain mapping. The model has been evaluated on a large data set. This paper shows uncertainty of the upstream boundary can have significant impact on the model results, exceeding the importance of model parameter uncertainty in some areas. However, this depends on the hydraulic conditions in the reach e.g. internal boundary conditions and, for example, the amount of backwater within the modelled region. The type of bridge implementation can have local effects, which is strongly influenced by the bridge geometry (in this case the area of the culvert). However, the type of bridge will not merely influence the model performance within the region of the structure, but also other evaluation criteria such as the travel time. This also highlights the difficulties in establishing which parameters have to be more closely examined in order to achieve better fits. In this study no parameter set or model implementation that fulfils all evaluation criteria could be established. We propose four different approaches to this problem: closer investigation of anomalies; introduction of local parameters; increasing the size of acceptable error bounds; and resorting to local model evaluation. Moreover, we show that it can be advantageous to decouple the classification into behavioural and non-behavioural model data/parameter sets from the calculation of uncertainty bounds.     

Pore scale modeling of immiscible and miscible fluid flows using smoothed particle hydrodynamics

A numerical model based on smoothed particle hydrodynamics (SPH) was developed and used to simulate immiscible and miscible fluid flows in porous media and to study effects of pore scale heterogeneity and anisotropy on such flows.     

A bed‐load transport model for rough turbulent open‐channel flows on plane beds

Data from flume studies are used to develop a model for predicting bed‐load transport rates in rough turbulent two‐dimensional open‐channel flows moving well sorted non‐cohesive sediments over plane mobile beds. The object is not to predict transport rates in natural channel flows but rather to provide a standard against which measured bed‐load transport rates influenced by factors such as bed forms, bed armouring, or limited sediment availability may be compared in order to assess the impact of these factors on bed‐load transport rates. The model is based on a revised version of Bagnold's basic energy equation i_bs_b = e_bω, where i_b is the immersed bed‐load transport rate, ω is flow power per unit area, e_b is the efficiency coefficient, and s_b is the stress coefficient defined as the ratio of the tangential bed shear stress caused by grain collisions and fluid drag to the immersed weight of the bed load. Expressions are developed for s_b and e_b in terms of G, a normalized measure of sediment transport stage, and these expressions are substituted into the revised energy equation to obtain the bed‐load transport equation i_b = ω G ^3·4. This equation applies regardless of the mode of bed‐load transport (i.e. saltation or sheet flow) and reduces to i_b = ω where G approaches 1 in the sheet‐flow regime. That i_b = ω does not mean that all the available power is dissipated in transporting the bed load. Rather, it reflects the fact that i_b is a transport rate that must be multiplied by s_b to become a work rate before it can be compared with ω. It follows that the proportion of ω that is dissipated in the transport of bed load is i_bs_b/ω, which is approximately 0·6 when i_b = ω. It is suggested that this remarkably high transport efficiency is achieved in sheet flow (1) because the ratio of grain‐to‐grain to grain‐to‐bed collisions increases with bed shear stress, and (2) because on average much more momentum is lost in a grain‐to‐bed collision than in a grain‐to‐grain one. Copyright © 2006 John Wiley & Sons, Ltd.     

裂隙岩质边坡地下水渗流数值模拟

本文从反倾向岩质高边坡的渗流特征和渗流普遍规律出发,根据渗流力学知识,建立数学模型,通过渗流数值模拟计算分析了渗透压力、水头变化规律以及对岩质边坡稳定性的影响.     

东海天然气水合物地热研究及其环境意义

依据地热资料研究天然气水合物稳定带厚度在东海海域的分布情况。东海在地质构造上位于新生代环太平洋构造带西部边缘岛弧的内侧,又是欧亚板块、太平洋板块和菲律宾海板块的相互作用带。依据国际热流委员会(IHFC)提供的东海地热数据,经过统计确定出该区域的热流分布,热流平均值为121·0mW/m2,最小值为73·0mW/m2,最大值为168·0mW/m2。同时利用天然气水合物温压模型计算了稳定带厚度,数据显示稳定带厚度平均值为92·2m,最小值为1·4m,最大值为190·6m,薄于其他已经发现的海洋天然气水合物稳定带厚度(约400m)。天然气水合物大部分分布在条件适宜的陆坡和岛坡上,冲绳海槽底部水合物稳定带厚度相对较薄。统计分析表明本区热流值与水合物稳定带厚度相关性很差,相关系数仅有0·12。这是由于天然气水合物所在海域水深较浅时,海底温度的变化迫使运算所应用的非线性方程影响因子迅速积累,从而导致相关系数降低。最后结合东海陆坡的地质条件,探讨了在天然气水合物存在的情况下,陆坡失稳的可能性及其造成的环境影响。     

基于GIS的滑坡、泥石流灾害危险性区划关键问题研究

随着GIS技术的引入,滑坡、泥石流灾害危险性区划的效率和准确性得以大大提高。依据工程地质类比原则,在灾害学理论指导下,结合专家打分、层次分析、人工神经网络、信息量、Logistic回归、统计量等模型方法,以MAPGIS软件为平台,利用C++语言开发了滑坡、泥石流灾害危险性区划评价分析系统;并重点探讨了GIS支持下的滑坡、泥石流灾害危险性区划过程中的因子分析、模型选取、模型复合、单元划分、系统集成、结果评价等关键问题,建立了一整套基于GIS的滑坡、泥石流灾害评价方法体系。应用该系统对长江三峡库区和辽宁省鞍山市分别开展了滑坡、泥石流灾害危险性区划研究,取得了较好的效果。