岩溶区多重介质水流模型研究进展

数值模拟技术广泛应用于地下水资源评价工作中,为水资源合理开发利用提供科学依据。本文针对岩溶区的数值模拟技术展开研究,综述了岩溶区地下水流建模技术的主要方法和研究进展,介绍了SWMM模型、UGRFLOW模型、CAVE模型、CFP模型等多重介质模型的建模方式和运行原理。岩溶区地下水赋存介质复杂,管道、裂隙和空隙多重介质并存,水流特征多样,等效连续介质模型不能很好地刻画岩溶区复杂的地下水流特征。随着对岩溶水流系统研究的深入,二重和三重介质模型逐渐应用于岩溶区地下水流模拟工作中,取得较好成果,但仍存在如忽略中宽裂隙的导水作用、岩溶管道设置单一、模型适用性不强等问题。今后应从解析岩溶含水介质的结构特征入手,关注中宽裂隙在地下水流系统中的导水作用和水流特征,改进管道流模块以更好地刻画复杂的地下河管道,改进不同介质间水流交换量算法,提高模型精确性和适用性。      

An empirical theory for gravitationally unstable flow in porous media

In this paper, we follow a similar procedure as proposed by Koval (SPE J 3(2):145–154, 1963) to analytically model CO₂ transfer between the overriding carbon dioxide layer and the brine layer below it. We show that a very thin diffusive layer on top separates the interface from a gravitationally unstable convective flow layer below it. Flow in the gravitationally unstable layer is described by the theory of Koval, a theory that is widely used and which describes miscible displacement as a pseudo two-phase flow problem. The pseudo two-phase flow problem provides the average concentration of CO₂ in the brine as a function of distance. We find that downstream of the diffusive layer, the solution of the convective part of the model, is a rarefaction solution that starts at the saturation corresponding to the highest value of the fractional-flow function. The model uses two free parameters, viz., a dilution factor and a gravity fingering index. A comparison of the Koval model with the horizontally averaged concentrations obtained from 2-D numerical simulations provides a correlation for the two parameters with the Rayleigh number. The obtained scaling relations can be used in numerical simulators to account for the density-driven natural convection, which cannot be currently captured because the grid cells are typically orders of magnitude larger than the wavelength of the initial fingers. The method can be applied both for storage of greenhouse gases in aquifers and for EOR processes using carbon dioxide or other solvents.     

自由跌流的水力特性

运用解析函数的Riemann-Hilbert混合边值问题理论,导出了计算重力作用下二维跌水流动的边界积分方程,并在物理平面上进行了数值求解.计算结果与Rouse试验资料进行了比较.     

Numerical modelling of the saline interface in coastal karstic aquifers within a conceptual model uncertainty framework

Numerical modelling is increasingly used as a tool for improving management strategies in aquifers and to support the design of comprehensive projects considering natural and anthropogenic processes. Overall, numerical simulation in karstic aquifers poses a major scientific challenge due to the non-Darcian groundwater flow dynamics. In specific cases, the equivalent porous medium approach has shown acceptable results, particularly in poorly karstified aquifers with regional/subregional scales such as this case. The Yucatan coastal karstic aquifer (Mexico) has been defined as a complex regional heterogeneous system, partially confined, thus allowing the discussion of multiple conceptual models. In this research, a two-dimensional numerical model of flow and transport was implemented using SEAWAT for the NW Yucatan aquifer. Four likely conceptual models were audited, calibrated and verified using hydrogeological field data, to select the best one, considering their fit and complexity. The numerical model accuracy was evaluated using the root-mean-square error, Nash Sutcliffe efficiency and the Pearson coefficient. The Akaike information criterion and Bayesian information criterion were included for evaluating the complexity of the numerical models. In addition, the signal of tide propagation into the aquifer was assessed as a proxy to improve the numerical calibration process. Results show that the most complex numerical model has a better calibration than the simpler models, but the model accuracy is worse when compared to less complex numerical models in the verification exercise. This research offers enhancement in the knowledge of numerical modelling in heterogeneous coastal aquifers within a conceptual-model uncertainty setting.

     

利用有限单元伽勒金法反求非线性流情况下砾石层的渗透系数

地下水流问题的数字计算机的应用,现虽处于高度发展阶段,但在采用这种数学物理模型时,首先要知道含水层的导水性能。过去利用水位观察资料或抽水试验资料计算水文地质参数的方法,只能适用于简单的含水体系,由于高度复杂的地质体并且间隙水流为非线性的这类参数如何确定还没有解决,如果利用近期发展起来的有限单元法来反求参数是一种理想的选择,特别是与古典的伽勒金法结合起来,既可简化数学处理过程,并有利于解非线性流问题。古典的伽勒金法原来只能适用于简易的分析域,而且要求介质为均一,但与空间离散的有限单元法结合起来,还便于解非均质而且边界不规则的问题。     

On the time scales of nonlinear instability in miscible displacement porous media flow

In this paper, we analyze the time scales associated with instable fingering induced by density contrasts in miscible displacement porous media flow. We perform numerical simulations of a two-dimensional domain with boundaries that are closed to flow and identify the three regimes of the dynamics, namely the development of a stable diffusive boundary layer, the onset and growth of instabilities, and the fully nonlinear dynamics. Special focus is given to the onset of the fully nonlinear regime. The results are generic in the sense that there are no parameters in the non-dimensional model problem. Large ensembles are studied and an error estimate is given based on the combined effect of numerical errors and sampling errors. The nonlinear time scales show a dependence on the size of initial perturbations. We estimate this size for three formations used for CO₂ storage and find that the onset of enhanced convective mixing is considerably delayed compared with the linear onset time.     

Groundwater origin and flow along selected transects in Ethiopian rift volcanic aquifers

The disruption of lithologies by cross-cutting faults and the variability in volcanic structures make the hydrogeology of the rifted volcanic terrain in Ethiopia very complex. Along two transects, selected due to their hydrogeologic characteristics, groundwater flow, depth of circulation and geochemical evolution have been conceptualized. The groundwater flow continuity between the high rainfall plateau bounding the rift and the rift valley aquifers depends principally on the nature of the bounding faults. Up to 50% of recharge to the rift aquifers comes from the plateau as groundwater inflow where the rift is cross cut by transverse fault zones. Recharge from the mountains is found to be insignificant where the rift is bounded by marginal grabens; channel loss and local precipitation are the principal sources of recharge to the rift aquifers in such cases. At a regional scale, there is a clear zonation in the geochemical compositions of groundwaters, the result of aquifer matrix composition differences. The environmental isotope results show that the majority of the aquifers contain modern groundwaters. In a few localities, particularly in thermal groundwaters representing deeper circulation, palaeo-groundwaters have been identified. Deeper groundwaters in the rift floor have a uniform ¹⁴C age ranging between 2,300 and 3,000 years.     

采用复合单元法模拟裂隙多孔介质变饱和流动

采用复合单元法建立了模拟裂隙多孔介质变饱和流动的数值模型。该模型具有以下特点：裂隙不需要离散成特定单元，而是根据几何位置插入到孔隙基质单元中形成复合单元；在复合单元中，分别建立裂隙流和孔隙基质流的计算方程，二者通过裂隙?基质界面产生联系并整合成复合单元方程；复合单元方程具有和常规有限单元方程相同的格式，因此，可以使用常规有限单元方程的求解技术。采用欠松弛迭代、集中质量矩阵以及自适应时步调节等技术，开发了裂隙多孔介质变饱和流动计算程序。通过模拟一维干土入渗和复杂裂隙含水层内的流动问题，验证了该模型的合理性和适用性。模拟结果为进一步认识非饱和裂隙含水层地下水流动特性提供了理论依据。     

含水岩组概化的累积导水系数法

地下水流系统理论和数值模拟技术分别是水文地质学的基本理论和技术方法,含水岩组的概化是地下水流系统分析和地下水数值模拟的重要基础,直接影响着数值模拟和水流系统分析的精度和可信度.为提高含水岩组概化的精度和可信度,提出一种含水岩组概化的新方法,即累积导水系数法.依据岩层厚度与渗透系数乘积累积值随深度的变化,以及水文地质剖面岩性分布的整体特征,概化含水介质结构.以玛纳斯河流域为例,应用该方法概化流域内的岩性剖面,结合GMS软件中TINS模块构建水文地质结构模型.结果表明,应用该方法概化后的含水层结构具有较好的合理性和仿真性,建立的三维模型很好地显示了研究区含水介质的空间展布特征,为建立地下水流模型奠定了良好的基础.      

山东淄博市地下水资源评价及其合理开发利用研究

本文采用有限单元法进行地下水资源评价以解决评价区供水面积大、精度要求高的矛盾。首先根据实际的开采量和长期观测孔的水位资料,校正数学模型,反求有关参数;然后给定开采方案进行未来水头的预测;最后研究地下水资源的合理开发利用问题。利用有限单元法进行水源地或单一含水层的地下水资源评价已被证明是先进的、有效的。本文是对水文地质条件非常复杂的大面积的区域,用有限单元法进行地下水资源评价的一种尝试,可供类似地区参考。     

基于GMS的三维TOUGH2模型及模拟

GMS和TOUGH2均是应用广泛的地下流动系统数值模拟软件,其中,GMS界面友好,功能强大;而TOUGH2虽具有强大的数值计算能力,但缺乏友好的可视化前后处理界面。本文借助GMS强大的前处理能力,基于概念模型建立了三维复杂模型,把GMS/MODFLOW三维数值模型(包括网格数据、岩性数据、初始和边界数据等)转化为TOUGH2数值模型进行数值计算。通过两个计算实例(含水层水平没有起伏和有起伏)对比分析了GMS/MODFLOW和TOUGH2计算结果的差异。结果显示,该方法可以快速地建立刻画复杂地质条件的TOUGH2模型,计算结果与GMS/MODFLOW差异很小,说明两个软件均有很高的可信度;同时,该方法发挥了两个软件各自的优势,为进行更为复杂的多相流动数值模拟提供了可行性。     

利用岸边水头动态确定含水层在临海方向上的边界

本文建立了潮汐效率与越流共同作用条件下承压水头动态的数学模型.利用复变函数的基本原理,考虑到地下水流动的连续性,推导出了相应的承压水头动态公式.以此动态公式为基础,提出了直接利用承压水头的潮汐动态来确定延伸到海底的含水层长度(边界条件)的计算方法.最后结合广西北海市的地下水头长观资料,对上述方法进行了应用和分析研究.     

Temporal variation of chemical and isotopic signals in major discharges of an alpine karst aquifer in Turkey: implications with respect to response of karst aquifers to recharge

Proper management of karst aquifers requires a better understanding of flow and transport mechanisms in these systems. Flow in karst aquifers is inherently very complex due to the non-linear and non-stationary relationship between recharge and discharge. Information on this relationship has been acquired for a large (1,000 km²), mountainous (>3,500 m asl) karst aquifer with a deep unsaturated zone (>2,000 m) in the Aladaglar mountain range of south-central Turkey. All major discharges from the aquifer, which drain almost all the recharge, have been observed periodically for specific electrical conductivity, tritium and oxygen-18 variations during a period of 12 months. Observations reveal that the system’s response to recharge depends strongly on the competition between the infiltration and drainage velocities. These velocities, which are controlled by variables such as the time of precipitation, time of infiltration, intensity, and continuity of recharge, determine the degree of dominance of different types of flow mechanisms in the aquifer. Bypass, well-mixed and piston flow mechanisms are used to explain the response of the aquifer to the spatio-temporal variations in recharge. It appears that the aquifer switches among these flow mechanisms depending on the prevailing recharge mode and the competition between infiltration and drainage velocities.     

Reconstruction and visualisation of melt topology in veined microdioritic enclaves

A number of xenocrystic microdioritic enclaves within the Ross of Mull Granite contain narrow veins of monzonitic material that have formed through the selective partial fusion of granitic magma mechanically incorporated into the microdiorite magma. Field evidence suggests that, as melting of the granitic material ensued, a monzonitic melt formed and exfiltrated from the surrounding microdioritic matrix resulting in Saffman–Taylor instabilities and the development of viscous fingering phenomena. The monzonitic melt network is now preserved as monzonitic veins. Mechanical serial sectioning through a representative veined microdioritic enclave shows that: (a) the topology of the veins is complex but crudely polyhedral, (b) there is porosity (monzonitic melt fraction) variation throughout the length of the veined enclave and (c) that the monzonitic veins are interconnected in three dimensions. Using image analysis software, the total porosity of individual veined enclave sections has been estimated, along with the variation in porosity with depth. Three dimensional rendering and animation software has been used to reconstruct and visualise the monzonitic melt network using Apple QuickTime™ software. The above method provides a cheap and powerful way of visualising remnant macroscopic transport properties (such as connectivity and melt topology) in partially molten igneous rocks.Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Using impacts of deep-level mining to research karst hydrology—a Darcy-based approach to predict the future of dried-up dolomitic springs in the Far West Rand goldfield (South Africa). Part 2: predicting inter-compartmental flow and final groundwater tables

Some of the world’s deepest goldmines operate below dolomitic karst aquifers in the Far West Rand (FWR) goldfield, South Africa. Associated impacts include the continuous dewatering of karst aquifers for over six decades and irreversible changes of the hydrogeological setting. Affecting an area of approximately 400 km² by drawing down the water table up to 700 m, these impacts, and the large amounts of data generated in the process, are used as unique research opportunities to better understand the complex karst hydrology. The focus of this study is on predicting final water table elevations in rewatered aquifers after mining ceases taking the fact that mines hydraulically linked previously disconnected aquifers into account. While part 1 of this series develops the conceptual model, this second part utilises large sets of pertinent data to calculate actual flow rates for predicting the fate of dried up springs after mine closure. Following a Darcy-based approach first applied by Swart et al. (Environ Geol 44:751–770, 2003a) it is not only predicted that the springs will flow again but also shown that linear relationships exist between flow rates through a combined system of karst-fractured aquifers overlying the mine void and the associated hydraulic head driving them. This suggests that—at this scale—porous media-based equations can be meaningfully used to predict flow in non-porous media.     

Evaluation of alternative conceptual models for groundwater modelling

This study evaluates the alternative conceptual models for groundwater modelling. A true model was created with a synthetic alluvial fan-plain hydrogeological framework. Various alternative conceptual models were evaluated for groundwater flow simulations. The first alternative model is a single aquifer layer model; the second alternative model is a 3-layer aquifer model; and the third model is a 5-layer model consisting of 3 aquifers separated by 2 aquitards. All models could fit very well to the observations with optimized values of hydraulic conductivities. However, the single aquifer layer model can only compute water balance components with good accuracy. The 3-layer aquifer model can be used for water balance computation and groundwater head simulation with small errors. The 5-layer model is capable of simulating water budget, groundwater head distribution and travel times with high accuracy. Multi-model analysis found only the 3rd alternative model superior.     

Karst terranes and environmental aspects

Karst is a complex geological phenomenon that relates to terranes composed of limestone, dolomites, gypsum, halite, or other soluble rocks. Protection of groundwater in karst against pollution is needed because of the high velocity of its flow (several hundreds of thousands of meters per day) and where polluted materials are carried without being filtered. Protection of karstic aquifers against pathogens and rapidly degradable chemicals is carried out for the catchment areas with an estimated delay time of 60 days and for a distance of more than 30 m. Tracer methods have recently become quite useful in karst regions; a tracer may respond as a multiimpulse because of different flow lines, flow velocities, and water aquifers. Therefore, if a total water balance has to be established, nearly all answer-back impulses should be considered.     

地下水污染脆弱性评价方法

地下水污染脆弱性是指污染物自顶部含水层以上某一位置到达地下水系统中某一特定位置的趋势和可能性,进一步分为固有脆弱性和特殊脆弱性。地下水污染脆弱性受地下水流系统和地球化学系统的影响和控制。其主要评价方法有主观分级评价法、统计或基于过程的评价法和综合评价法三大类。中国地下水污染脆弱性评价已有很好的工作基础,评价工作中应以地下水系统为单元,以饮用水井、集中供水水源地、区域含水层系统的补给区为重点保护目标,评价方法应综合区域地下水流系统的过程分析和指数评价方法,并利用已有的区域水质资料进行检验,增强评价结论的科学性和可靠性。     

Geoelectrical Schlumberger investigation for characterizing the hydrogeological conditions using GIS in Omalur Taluk,Salem District,Tamil Nadu,India

The increasing demand for freshwater has necessitated the exploration for new sources of groundwater, particularly in hard rock terrain, where groundwater is a vital source of freshwater. A fast, cost effective, and economical way of exploration is to study and analyze geophysical resistivity survey data. The present study area Omalur taluk, Salem District, Tamil Nadu, India, is overlain by Archaean crystalline metamorphic complex. The study area is a characteristic region of unconfined aquifer system. The potential for occurrence of groundwater in the study areas was classified as very good, good, moderate, and poor by interpreting the subsurface geophysical investigations, namely vertical electrical soundings, were carried out to delineate potential water-bearing zones. The studies reveal that the groundwater potential of shallow aquifers is due to weathered zone very low resistivity and very high thickness and the potential of deeper aquifers is determined by fracture zone very low resistivity and very high thickness area. By using conventional GIS method, the spatial distribution maps for different layer (top soil, weathered zone, first fracture zone, and second fracture zone) thicknesses were prepared. The geoelectrical approach was successfully applied in the study area and can be therefore easily adopted for similar environments.     

Hydrogeological model of the Baltic Artesian Basin

The Baltic Artesian Basin (BAB) is a complex multi-layered hydrogeological system in the south-eastern Baltic covering about 480,000 km². The aim of this study is to develop a closed hydrogeological mathematical model for the BAB. Heterogeneous geological data from different sources were used to build the geometry of the model, i.e. geological maps and stratigraphic information from around 20,000 boreholes. The finite element method was used for the calculation of the steady-state three-dimensional (3D) flow of unconfined groundwater. The 24-layer model was divided into about 1,000,000 finite elements. A simple recharge model was applied to describe the rate of infiltration, and the discharge was set at the water-supply wells. Variable hydraulic conductivities were used for the upper (Quaternary) deposits, while constant hydraulic conductivity values were assumed for the deeper layers. The model was calibrated on the statistically weighted borehole water-level measurements, applying L-BFGS-B (automatic parameter optimization method) for the hydraulic conductivities of each layer. The principal flows inside the BAB and the integral flow parameters were analyzed. The modeling results suggest that deeper aquifers are characterized by strong southeast–northwest groundwater flow, which is altered by the local topography in the upper, active water-exchange aquifers.