非均质含水层中渗流与溶质运移研究进展简

Natural aquifer heterogeneity controls the groundwater flow and solute transport, and how to accurately quantify the flow and solute transport in heterogeneous aquifers has received wide attention by many scholars, and has become a hot research topic in earth science. Theoretically, a systematic review is given by the following aspect: flow and solute transport model, moment analysis, multi scale analysis. The resolved and remained issues for scale conversion in hydrogeology research are pointed out. Secondly, recent advances of heterogeneous field test, uncertainty and velocity connectivity are analyzed. Finally, the geophysical inversion of aquifer heterogeneity, stochastic theory and development of stochastic simulation software, scale conversion and uncertainty of velocity connectivity, and the relationship between heterogeneity and hydrogeological condition on the major four aspects of the future research direction is pointed out.     

A Simple Approach to Account for Radial Flow and Boundary Conditions When Kriging Hydraulic Head Fields for Confined Aquifers

The estimation and mapping of realistic hydraulic head fields, hence of flow paths, is a major goal of many hydrogeological studies. The most widely used method to obtain reliable head fields is the inverse approach. This approach relies on the numerical approximation of the flow equation and requires specifying boundary conditions and the transmissivity of each grid element. Boundary conditions are often unknown or poorly known, yet they impose a strong signature on the head fields obtained by inverse analysis. A simpler alternative to the inverse approach is the direct kriging of the head field using the measurements obtained at observation wells. The kriging must be modified to incorporate the available information. Use of the dual kriging formalism enables simultaneously estimating the head field, the aquifer mean transmissivity, and the regional hydraulic gradient from head data in steady or transient state conditions. In transient state conditions, an estimate of the storage coefficient can be obtained. We test the approach on simple analytical cases, on synthetic cases with solutions obtained numerically using a finite element flow simulator, and on a real aquifer. For homogeneous aquifers, infinite or bounded, the kriging estimate retrieves the exact solution of the head field, the exact hydrogeological parameters and the flow net. With heterogeneous aquifers, kriging accurately estimates the head field with prediction errors of the same magnitude as typical head measurement errors. The transmissivities are also accurately estimated by kriging. Moreover, if inversion is required, the kriged head along boundaries can be used as realistic boundary conditions for flow simulation.     

含水层中考虑水—岩间阳离子交换的溶质运移方程

导出了能反映含水层中水—岩间阳离子交换作用的溶质运移方程,并提出了一个新的海水入侵数学模型.模型以交换阳离子Na+和Ca2+作为模拟因子,考虑了海水入侵过程中水—岩间的Na+—Ca2+交换,模型用于模拟山东省龙口市黄河营地区的海水入侵,取得了令人满意的结果,从而证明导出的溶质运移方程是正确的,建立的模型是可靠的.     

岩溶裂隙地下水流数学模型求解的有限体积法及应用

岩溶裂隙介质可视为双重孔隙度介质。根据水流连续性原理、质量守衡和达西定律,建立了裂隙-岩溶承压含水层非稳定流双重介质数学模型,采用有限体积法对其求解,并以深圳龙岗区为例对模型进行验证。结果表明,计算模型与实际水文地质条件比较接近,计算水位与实测水位相吻合。     

地下水溶质迁移模拟研究进展

在地下水的相关研究中,农药和石油等地下水污染、土地盐碱化、海水入侵等诸多实际问题主要的研究方法都涉及地下水溶质迁移模拟. 相比地下水水流模拟的相对完善,对溶质迁移的模拟比较薄弱且迁移过程本身复杂性较高,目前地下水溶质迁移的研究工作还处在全面发展的阶段. 文中阐述了反映地下水溶质迁移机理和过程的数学模型,综述了溶质迁移模拟在地下水污染物防治、土地盐碱化、海水入侵、石油和放射性废物扩散等问题的诸多应用,归类了目前溶质迁移模拟所使用的对流迁移、对流-弥散模拟等主要数值方法,并对这些方法的优缺点和应用实例做了总结. 最后,分析了目前溶质迁移模拟中存在的不足,展望了未来在参数确定、裂隙介质运移机理和多相介质条件下运移模拟可能取得的突破.     

水箱-管道系统溶质运移实验研究及其岩溶水文地质意义

为了研究岩溶管道中溶潭对溶质运移的影响,在实验室内构建水箱-管道系统,在不同管道结构和水流条件下进行定量示踪实验并得到相应的穿透曲线(BTCs);采用Qtracer2软件分析溶质运移参数,采用滞后系数R分析实验结果与一维经典对流弥散方程解析解之间的差别。实验结果显示:随着水箱数量的增加,示踪剂(NaCl)峰值质量浓度逐渐降低,弥散系数和弥散度逐渐增加,穿透曲线拖尾逐渐增长,表明水箱的瞬态存储使溶质运移滞后;与不对称水箱相比,对称水箱BTC拖尾较长;峰现时间随着不对称水箱数量的增多明显滞后;出口流量增加时,弥散度减小,BTC拖尾变短。一维经典对流弥散方程解析解仅对单管道最大流量条件下的BTC拟合较好,对流量较小的单管道和水箱-管道系统的BTC拟合较差,需研究适用的模型解释其拖尾现象。     

Optimal Groundwater Development in Coastal Aquifers Near Beihai, China

INTRODUCTIONThe city of Beihai,located on the south coast ofGuangxi,China,relies heavily on groundwater for its potablewater supply and agricultural irrigation.With rapid increasein population (for instance,from 134 0 0 0 in 1987to 47930 0in1995 ) and in developm ent program s,the demand for freshwater has been growing. Approxim ately 170 0 0 0 m3/ d ofgroundwater has been pumped from the productive coastalaquifers in recent years.Contamination of the fresh water inthe coastal aquifers b…     

云南蒙自断陷盆地岩溶地下水流向研究

岩溶地区地下水径流特征复杂。以云南蒙自断陷盆地为例,通过开展水文地质调查、示踪试验、综合研究等,揭示了不同水文地质单元岩溶地下水径流方向:蒙自盆地内东山—蒙自(Ⅰ-1)、大屯(Ⅰ-2),蒙自盆地东南黑龙潭(Ⅰ-3)单元内岩溶地下水整体由南东流向北西;新安所(Ⅰ-4)单元内由南向北;蒙自盆地南部8号暗河(Ⅱ-1)、9号暗河(Ⅱ-2)两个水文地质单元内岩溶地下水由北流向南,对研究西南典型岩溶区地下水径流特征、合理评价和开发地下水资源具有重要的理论意义和实际应用价值。     

水文地球化学信息对岩溶地下水流动系统特征的指示意义--以山西神头泉域为例

神头泉域岩溶水的水化学特征显示,岩溶水样的总溶解固体质量浓度、Sr2 质量浓度及宏量组分质量浓度均由泉域北、西、南部向排泄区增大,与岩溶水的流动方向一致.根据c(Sr)/c(Ca)与总溶解固体质量浓度可将神头泉群划分为3个泉组,对应于3个不同的岩溶水子系统,岩溶水样的聚类分析结果与这一结论相符合.研究表明水文地球化学信息对地下水系统的水动力特征分析具有指示意义.     

桂林东区峰林平原岩溶地下水示踪试验与分析

为查明桂林东区峰林平原地下水系统的结构特征,在岩溶水文地质调查的基础上,通过人工取样并采用高分辨率野外荧光仪,对桂林东区峰林平原区投放的示踪剂荧光素钠的接收情况进行检测。结果表明：峰林平原区地下水的出露比较多,主要以沟渠、湖塘、民井等为主,其中1号点初现时间为第22天,最大视流速为19.5 m/d,平均视流速4.9 m/d,说明本试验段内地下径流为典型的面状流态,岩溶含水介质比较均匀。峰林平原区岩溶含水层结构岩溶化程度很高,呈网络状,存在统一的地下水位,岩溶发育较均一,同时也存在较大的裂隙甚至较大规模的管道和地下河系统,容易形成岩溶地下水优先流。     

晋陕峡谷地区北段岩溶地下水流特征分析

晋陕峡谷地区北段岩溶地下水流系统比较复杂,尽管在该地区开展了不少地下水资源的勘查工作,但仍然存在不少有争议的问题.本次研究利用水文地球化学方法对研究区内的清水河地区、准格尔地区、偏关-龙口地区和天桥地区的岩溶地下水系统特征以及黄河水对岩溶地下水系统的影响进行了分析,对区内岩溶地下水流系统取得较为清晰的认识.     

宁夏南部“南北古脊梁”岩溶裂隙水流系统分析

运用构造控水分析、水化学同位素等方法,对宁南“南北古脊梁”岩溶裂隙水系统进行了深入的分析与讨论。该区储水空间以岩溶裂隙为主,岩溶水的空间分布明显受SN向大型断裂构造控制。岩溶地下水以大气降水起源为主,并表现为多元水混合而成。水质分布呈现南优北劣的分带特征,北部水一岩相互作用形式为溶滤一蒸发浓缩型,呈高矿化咸水,南部为溶滤一混合型,呈低矿化淡水。根据水动力场和水化学场特征,划分了3个相对独立的岩溶水子系统。     

Dealing with the Spatial Synthetic Heterogeneity of Aquifers in the North China Plain: A Case Study of Luancheng County in Hebei Province

The complexity of alluvial-pluvial fan depositional systems makes the detailed characterization of their heterogeneity difficult, yet such a detailed characterization is commonly needed for construction of reliable groundwater models. Traditional models mainly focus on using a single aquifer property to qualitatively or semi-quantitatively characterize the heterogeneity of aquifer, so that they are unable to quantitatively reflect the synthetic heterogeneity of all aquifer properties. In this paper, we propose the heterogeneity synthetic index (HSI) for quantitative characterization of synthetic heterogeneity of an aquifer. The proposed calculation process involves four steps: (1) estimation of the hydraulic conductivity of a sediment sample using the cloud-Markov model, (2) establishment of the sedimentary microfacies distribution model through the Markov chain, (3) characterization of the distribution model of hydrogeological parameters using the improved sequential simulation method according to the “facies-controlled modeling” technique, and (4) application of the entropy weight method to calculate the weight coefficient of the above aquifer properties. The HSI of an aquifer is calculated by superposition of these models according to the corresponding weight coefficient. This approach was applied to the Luancheng aquifer deposit in the southeast Hutuo River alluvial-pluvial fan in the North China Plain (NCP). The results have demonstrated that aquifer 3 which was formed in the middle Pleistocene has the strongest heterogeneity, with an HSI of 0.25–0.75. Aquifer 4 formed in the early Pleistocene shows an intermediate heterogeneity, with the HSI ranging 0.35–0.75. The weakest heterogeneity was found in aquifers 1 and 2 formed in the Holocene and late Pleistocene, with HSI values of 0.40–0.75 and 0.40–0.80, respectively. The heterogeneity of all the four aquifers is relatively strong in the radial direction of the Huai River alluvial-pluvial fan due to the abrupt change of microfacies. In contrast, in the radial direction of the Hutuo River alluvial-pluvial fan, the microfacies change mildly, and the continuity of hydrogeological parameters is better, which has resulted in weaker heterogeneity of the four aquifers in this direction. Findings suggest that the sedimentary environment has significant effects on the aquifer heterogeneity. Considering that there are many aquifer properties, HSI can quantitatively characterize the synthetic heterogeneity of the aquifer and describe the influence of each aquifer property on the synthetic heterogeneity of the aquifer according to its weight coefficient. Thus the HSI approach can be successfully used to deal with the spatial heterogeneity of aquifer and provide a foundation for studies on contaminant transport.     

岩溶含水系统降水入渗补给研究进展

岩溶含水系统中赋存着丰富的优质地下水,而大气降水是浅部可供开采的岩溶地下水的最主要补给来源。受岩溶含水系统各向异性、不均一性和直接观测难度大等因素的影响,降水入渗补给量的计算是一个非常复杂的过程。确定岩溶含水系统的汇水范围是降水入渗补给计算的首要问题,示踪法与经验公式法被证明是最有效的两种方法。降水入渗补给量的计算方法主要包括水文过程线法分析法、氧同位素法、氯质量平衡法、基于GIS的多变量综合分析法和模型法。本文对目前岩溶含水系统降水入渗补给计算方法的关键点和适用条件进行了总结和对比,同时指出大气降水物理化学性质的时空特征以及水-岩反应可以作为未来研究岩溶含水系统降水补给的研究方向之一。     

Sensitivity Analysis of a Combined Groundwater Flow and Solute Transport Model Using Local-Grid Refinement: A Case Study

Combining groundwater flow models with solute transport models represents a common challenge in groundwater resources assessments and contaminant transport modeling. Groundwater flow models are usually constructed at somewhat larger scales (involving a coarser discretization) to include natural boundary conditions. They are commonly calibrated using observed groundwater levels and flows (if available). The groundwater solute transport models may be constructed at a smaller scale with finer discretization than the flow models in order to accurately delineate the solute source and the modeled target, to capture any heterogeneity that may affect contaminant migration, and to minimize numerical dispersion while still maintaining a reasonable computing time. The solution that is explored here is based on defining a finer grid subdomain within a larger coarser domain. The local-grid refinement (LGR) implemented in the Modular 3D finite-difference ground-water flow model (MODFLOW) code has such a provision to simulate groundwater flow in two nested grids: a higher-resolution sub-grid within a coarse grid. Under the premise that the interface between both models was well defined, a comprehensive sensitivity and uncertainty analysis was performed whereby the effect of a parameter perturbation in a coarser-grid model on transport predictions using a higher-resolution grid was quantified. This approach was tested for a groundwater flow and solute transport analysis in support of a safety evaluation of the future Belgian near-surface radioactive waste disposal facility. Our reference coarse-grid groundwater flow model was coupled with a smaller fine sub-grid model in two different ways. While the reference flow model was calibrated using observed groundwater levels at a scale commensurate with that of the coarse-grid model, the fine sub-grid model was used to run a solute transport simulation quantifying concentrations in a hypothetical well nearby the disposal facility. When LGR coupling was compared to a one-way coupling, LGR was found to provide a smoother flow solution resulting in a more CPU-efficient transport solution. Parameter sensitivities performed with the groundwater flow model resulted in sensitivities at the head observation locations. These sensitivities identified the recharge as the most sensitive parameter, with the hydraulic conductivity of the upper aquifer as the second most sensitive parameter in regard to calculated groundwater heads. Based on one-percent sensitivity maps, the spatial distribution of the observations with the highest sensitivities is slightly different for the upper aquifer hydraulic conductivity than for recharge. Sensitivity analyses were further performed to assess the prediction scaled sensitivities for hypothetical contaminant concentrations using the combined groundwater flow and solute transport models. Including all pertinent parameters into the sensitivity analysis identified the hydraulic conductivity of the upper aquifer as the most sensitive parameter with regard to the prediction of contaminant concentrations.     

非均质含水层中地下水年龄数值模拟

地下水年龄数值模拟是近年来国际地下水领域的一个全新研究方向。在简单介绍地下水年龄数值模拟相关理论和方法基础上,通过一个简单算例模拟对比了稳定流条件下不同非均质含水层中的水流运动和年龄分布,探讨了夹层透水性的强弱对地下水年龄分布的影响。结果表明,在地下水系统保持补径排通量不变的前提下,弱透水夹层的存在对于特定区域地下水的深部循环和更新具有积极作用。同时,在模拟区域深部含水层的年龄分布时,不能忽略其弥散作用。地下水年龄模拟为研究区域尺度上的地下水循环和地下水资源可持续利用提供了一种新的途径。     

基于CFP的岩溶管道流数值模拟研究——以桂林寨底地下河子系统为例

为了描述岩溶多重含水介质存在非达西流问题,研究管道流CFP模型在岩溶区数值模拟中的应用。在分析管道流CFP基本原理的基础上,通过建立概念模型算例(5个管道,上下临界雷诺数分别为2 000和4 000),探讨不同管道水文地质参数对出口流量及水流状态的影响,并对参数敏感性进行分析,最后将CFPM1模型应用于广西寨底岩溶地下河系统,探讨在实际应用中的管道流模型适用性。概念模型结果表明,参数敏感度从大到小依次为:管道直径、水力梯度、管道渗透系数、弯曲度及粗糙度,流量与雷诺数随管道参数增大而增大,水流状态从层流变为紊流。寨底岩溶管道流模型中3个观测井及地下河出口实测流量表明CFPM1能较好地模拟出岩溶地下水水位变化趋势。得出结论:CFPM1管道流模型允许岩溶管道与基岩含水层进行水流交换,能够较好地刻画岩溶区管道介质非达西流特征,但管道参数获取较为困难,且不能刻画岩溶管道形态变化。该方法具有一定通用性,可为实际应用提供指导。     

硝基苯在含水层中迁移转化影响因素的定量研究

通过室内模拟实验,确定了影响硝基苯在含水层中迁移转化的主要因素,并采用解析解方法计算出各因素的影响水平。研究表明,硝基苯在淤泥介质中的吸附能力大于砂砾介质,其在淤泥介质中的阻滞系数为2.22,在砂砾介质中的阻滞系数为1.74。根据一级动力学方程计算得出,硝基苯在淤泥介质中的自然衰减总速率常数为0.063 5/h,在砂砾介质中的自然衰减总速率常数为0.0059/h。而且,不论是在淤泥介质还是在砂砾介质中,生物降解作用均是硝基苯迁移转化过程的主要影响因素,其影响水平占硝基苯自然衰减能力的50%左右。     

重庆岩溶地下水氢氧稳定同位素地球化学特征

重庆地区分布有380条岩溶地下河,是重庆市重要的水资源。为掌握岩溶地下河水稳定同位素地球化学特征及其环境意义,研究了重庆市不同地区51条地下河水体的稳定同位素地球化学特征。研究表明,重庆市岩溶地下河旱、雨季δ18O、δD值均沿大气降水线分布,表明地下河水均起源于大气降水。受雨季降水云团运动规律(环流效应)和区域地形的影响,地下河水δ18O、δD 值雨季表现出渝东北地区(渝西地区,渝东地区)<渝东南地区的明显区域分布规律(“<”表示偏负于),旱季由于地下河水在含水层中运动较慢,δ18O、δD值的区域性规律不明显,且由于具有较雨季长的滞留时间,导致其d-excess值明显小于雨季。利用岩溶地下水δ18O值和区域高程建立了二者之间的二元回归模型,揭示了重庆岩溶地下河水旱季δ18O值随高度的变化率为-0.34‰/100 m,雨季为-0.31‰/100 m,这对于区域水循环研究具有重要意义。     

大同煤田岩溶地下水的水文地质特征分析

根据大同煤田岩溶地下水化学特征资料、钻孔揭露岩溶发育的资料、以及矿井出水资料，结合岩溶地下水区域水文地质条件、岩溶地下水的水温与同位素特征，对大同煤田岩溶地下水的富水层位与找水方向作了进一步的探讨与分析。