纵向弥散作用与渗透介质非均质性定量关系的模拟研究

在地下水污染模拟预报中,弥散参数是很难确定的一个模型参数。因实验室小尺度弥散规律一般不能用于大尺度弥散过程,而野外示踪试验却耗资大、周期长,限制了其实用性。文中利用随机数值模拟手段、基于随机理论的蒙特卡罗方法及序贯高斯模拟技术来生成渗透系数随机场,并研究渗透系数对数场的方差、相关长度以及变异函数类型在不同尺度上对纵向弥散度的影响,进而建立纵向弥散度与随机分布渗透系数场的方差和相关长度的统计定量关系,并与Gelhar理论计算结果进行比较。数值模拟结果表明,经过一定迁移距离后纵向弥散度与随机分布渗透系数对数场的方差和相关长度具有良好的线性统计关系,与Gelhar理论公式表达的关系类型类似。但对于较大的方差,纵向弥散度模拟结果明显大于Gelhar理论计算值,而对于较大相关长度在迁移距离不很大时,纵向弥散度模拟结果明显小于Gelhar理论计算值。本研究可为野外大尺度地下水污染预报模型中水动力弥散参数的确定提供方法借鉴。     

裂隙网络管道模型弥散试验

为了探求不同裂隙几何参数对裂隙网络溶质运移的影响,基于离散裂隙网络思想和优势流、沟槽流理论,建立裂隙管道网络概念模型,搭建不同管径、不同连通方式的管道网络试验装置,进行渗流和溶质运移实验。运用应用广泛的模拟软件CHEMFLO-2000建立等效多孔介质模型,拟合不同几何参数下等效弥散度,定量刻画不同管道网络几何参数对溶质运移的影响,讨论了不同管径、连通方式等与等效弥散度之间的关系。通过进一步分析得出:在连通方式相同的情况下,不同管径的管道网络等效弥散度存在差异,但是差异不大;溶质在小管径中的穿透时间短于大管径,穿透曲线缓和程度相差不大;管道网络连通方式越复杂,其等效弥散度越大、对溶质运移的影响越大、穿透曲线越缓和、路径越长,等效弥散度越大;用这种等效弥散度的方法表征管道网络对溶质运移的影响,与多孔介质弥散度具有相似性;管道数目、管道面数目与等效弥散度成正相关关系,且等效弥散度随尺度的增大而增加。     

An experimental and theoretical study of the seepage migration of suspended particles with different sizes

This study experimentally investigates the effect of particle size, particle concentration and flow velocity on the migration of suspended particles of size 1.02–47 μm in porous media. The results show that at the same flow velocity, the peak values of the breakthrough curves decrease and corresponding pore volumes increase slightly with increasing particles size. The migration velocity of smaller suspended particles is even greater than water flow velocity, which is attributed to the size exclusion effect. With increase of the injected particle concentration, the deposition coefficients of small single particles increase at first and then tend to a steady state or even decrease slightly, explained by the maximum retention concentration. The dispersivity of small particles decreases with increasing velocity. However, at a high flow velocity, the hydrodynamic dispersivity becomes increasingly dominant with the increase of particle size. The deposition coefficients for large-sized particles are higher than those for small-sized particles, which is attributed to considerable mass removal due to straining. An analytical solution, considering the release effect of sorbed particles, is developed to account for the one-dimensional flow and dispersive effect using a source function method, and then three transport parameters—dispersivity, deposition coefficient and release coefficient—are fitted using the experimental results. Finally, suspended-particle migration is predicted by the proposed model for short-time constant-concentration injection and repeated three-pulse injection. Overall, particle size has a significant effect on the seepage migration parameters of suspended particles in porous media such as the particle velocity, dispersivity and deposition coefficient.     

Assessment of key transport parameters in a karst system under different dynamic conditions based on tracer experiments: the Jeita karst system,Lebanon

Artificial tracer experiments were conducted in the mature karst system of Jeita (Lebanon) under various flow conditions using surface and subsurface tracer injection points, to determine the variation of transport parameters (attenuation of peak concentration, velocity, transit times, dispersivity, and proportion of immobile and mobile regions) along fast and slow flow pathways. Tracer breakthrough curves (TBCs) observed at the karst spring were interpreted using a two-region nonequilibrium approach (2RNEM) to account for the skewness in the TBCs’ long tailings. The conduit test results revealed a discharge threshold in the system dynamics, beyond which the transport parameters vary significantly. The polynomial relationship between transport velocity and discharge can be related to the variation of the conduit’s cross-sectional area. Longitudinal dispersivity in the conduit system is not a constant value (α?=?7–10 m) and decreases linearly with increasing flow rate because of dilution effects. Additionally, the proportion of immobile regions (arising from conduit irregularities) increases with decreasing water level in the conduit system. From tracer tests with injection at the surface, longitudinal dispersivity values are found to be large (8–27 m). The tailing observed in some TBCs is generated in the unsaturated zone before the tracer actually arrives at the major subsurface conduit draining the system. This work allows the estimation and prediction of the key transport parameters in karst aquifers. It shows that these parameters vary with time and flow dynamics, and they reflect the geometry of the flow pathway and the origin of infiltrating (potentially contaminated) recharge.     

土柱变密度溶质运移试验推求水平及垂直方向弥散度研究

基于溶质运移对流弥散理论的变密度海水入侵模型广泛用于海水入侵研究,而水动力弥散系数是影响模型模拟效果的关键性参量之一。利用传统土柱溶质运移试验结合旁侧抽水,采用数值反演法成功获取了水平及垂直两方向弥散度。相对于传统方法,在不增加试验复杂度的前提下,同时推求了不同方向的弥散度,提高了试验效率,节约了试验成本,可广泛用于测定水动力弥散系数等参数。     

Spatial resolution of transport parameters in a subtropical karst conduit system during dry and wet seasons

Karst aquifers are characterized by a high degree of hydrologic variability and spatial heterogeneity of transport parameters. Tracer tests allow the quantification of these parameters, but conventional point-to-point experiments fail to capture spatiotemporal variations of flow and transport. The goal of this study was to elucidate the spatial distribution of transport parameters in a karst conduit system at different flow conditions. Therefore, six tracer tests were conducted in an active and accessible cave system in Vietnam during dry and wet seasons. Injections and monitoring were done at five sites along the flow system: a swallow hole, two sites inside the cave, and two springs draining the system. Breakthrough curves (BTCs) were modeled with CXTFIT software using the one-dimensional advection-dispersion model and the two-region nonequilibrium model. In order to obtain transport parameters in the individual sections of the system, a multi-pulse injection approach was used, which was realized by using the BTCs from one section as input functions for the next section. Major findings include: (1) In the entire system, mean flow velocities increase from 183 to 1,043 m/h with increasing discharge, while (2) the proportion of immobile fluid regions decrease; (3) the lowest dispersivity was found at intermediate discharge; (4) in the individual cave sections, flow velocities decrease along the flow direction, related to decreasing gradients, while (5) dispersivity is highest in the middle section of the cave. The obtained results provide a valuable basis for the development of an adapted water management strategy for a projected water-supply system.     

Mass transport modelling for assessement of groundwater contamination around Mathura oil refinery, Mathura, Uttar Pradesh, India

Mathura oil refinery was commissioned during 1977 and effluent storage ponds were constructed at the same time. These storage ponds receive wastewater from the refinery at a rate of 10000 m³/day. After treatment, waste water is discharged through a 3-km pipeline to a stream leading to the Yamuna river. The groundwater-monitoring as well as water-quality monitoring was carried at 24 observation wells in the refinery site during 1997. The water quality measurements indicated total dissolved chloride and sulphate concentration of native groundwater as 400 mg/l, whereas elevated levels up to 600 mg/l were found at the wells close to polishing ponds. Thus combined transport of chloride and sulphate was simulated in the mass transport model. A three-dimensional flow, pathlines and mass transport model of the aquifer system were constructed to analyze the impact of seepage from polishing ponds contaminating the groundwater regime. The permeability of aquifer varies between 1.5–2.5 m/day. The porosity of formation was assumed as 0.2. The constant head and constant concentration boundaries were assigned to the nodes representing effluent storage ponds. Longitudinal dispersivity of 100 m, horizontal transverse dispersivity of 10 m and vertical transverse dispersivity of 0.01 m were assumed. The mass transport model was calibrated for 20 years by comparing total dissolved chloride and sulphate concentrations from 1997. The model predictions indicate further migration of contaminants on the east of effluent ponds in future. Received: 4 January 1999 · Accepted: 12 July 1999     

Migration of lignosulfonates in a karstic-fractured-porous aquifer: History and prognosis for a Zn-Pb mine,Pomorzany, southern Poland

Lignosulfonate lyes, which had been deposited by a paper factory for about 50 years into Quaternary sediments and accumulated in a quasistagnant depression, became mobile when a depression cone caused by mine activity in a deeper aquifer reached the erosion window close to the disposal site. Their transport parameters in the karstic-fractured-porous carbonate aquifer are shown to be mainly governed by diffusion between the mobile water in karstic channels and fractures and the stagnant water in matrix, although the hydraulic conductivity is shown to depend mainly on the karstic system. The main transport velocity is, with a sufficient accuracy, related to the hydraulic conductivity by a simple model in which matrix porosity is the governing rock parameter. The dispersivity, although also governed by matrix diffusion, depends strongly on several parameters that are difficult to determine.     

单井注抽试验测算地下水流速的数值分析

为了分析注入阶段地下水流速及含水层弥散度对单井注抽试验解析模型测算地下水流速的影响机理,通过采用GMS(groundwater modeling system)软件建立了单井注抽试验数值模型,通过与Leap and Kaplan(1988)近似解析模型计算结果对比分析上述因素对解析模型计算结果的影响.研究结果表明:地下...     

西宁市贵德县河滨公园林场第四系含水层弥散试验研究

地下水污染问题日益严重,研究溶质运移的弥散理论开始应用于实际问题。建立地下水溶质运移模型,对地下水中污染物的运移及发展趋势进行准确预测,是对地下水进行保护、对地下水污染进行控制的基础。而弥散参数的确定则是地下水溶质运移模型建立的关键环节之一,直接影响着模型预测结果的精度和准确性。 对西宁市贵德县地下水污染的水质运移规律进行分析,在贵德县河滨公园林场采用径向收敛流水动力弥散理论方法进行了第四系含水层现场弥散试验,计算了试验场地潜水含水层的弥散度,获得纵向弥散度(a_L)为0.843~0.998 cm,横向弥散度(a_T)经验推断值为0.17~0.20 cm,为进一步建立该地区的地下水溶质运移模型、预测地下水污染的发展趋势和评价该地区地下水环境质量提供了数据参考。     

Use of mathematical modeling to investigate inter-aquifer contamination by organic-rich water through an unplugged well (central Wielkopolska, Poland)

Hydrogeological research was performed to establish the origin of organic-rich water in a Quaternary aquifer (Wielkopolska Buried Valley aquifer, Poland). The spatial distribution of such water was determined by means of multilevel sampling with a packer. The nature, spatial distribution and chemistry of the organic-rich water suggest hydraulic connection between the Quaternary aquifer and the underlying Neogene aquifer. This connection may be enabled by an old abandoned, improperly plugged well, and must therefore be considered as an artificial hydraulic connection between hydrogeologically separated aquifers. For the verification of this hypothesis, a groundwater contaminant transport model was constructed. The results of contaminant transport modeling allow the contaminant transport parameters (transverse and longitudinal dispersivity) to be identified in the field. The numerical groundwater model was used for the prognosis of organic-rich water remediation, and for specification of water-resource management principles in the region where organic-rich water occurs. The results of this research enable the verification of assumptions regarding complex groundwater flow conditions in the boundary zone of a buried valley.     

Effect of tidal fluctuations on transient dispersion of simulated contaminant concentrations in coastal aquifers

Variable-density groundwater models require extensive computational resources, particularly for simulations representing short-term hydrologic variability such as tidal fluctuations. Saltwater-intrusion models usually neglect tidal fluctuations and this may introduce errors in simulated concentrations. The effects of tides on simulated concentrations in a coastal aquifer were assessed. Three analyses are reported: in the first, simulations with and without tides were compared for three different dispersivity values. Tides do not significantly affect the transfer of a hypothetical contaminant into the ocean; however, the concentration difference between tidal and non-tidal simulations could be as much as 15%. In the second analysis, the dispersivity value for the model without tides was increased in a zone near the ocean boundary. By slightly increasing dispersivity in this zone, the maximum concentration difference between the simulations with and without tides was reduced to as low as 7%. In the last analysis, an apparent dispersivity value was calculated for each model cell using the simulated velocity variations from the model with tides. Use of apparent dispersivity values in models with a constant ocean boundary seems to provide a reasonable approach for approximating tidal effects in simulations where explicit representation of tidal fluctuations is not feasible.     

Inverse modeling of variable-density groundwater flow in a semi-arid area in Iran using a genetic algorithm

Flow and mass transport parameter estimation was done by creating an inverse model of a seawater intrusion system using a genetic algorithm (GA) method as the optimization procedure. Firstly, the SEAWAT code was used for the forward solution part and then a program was written in MATLAB for coupling the forward and inverse processes. The auto-calibration objective function was defined with the root mean square errors (RMSE) between the observed and the simulated values. A simple GA was used to minimize the RMSE criterion. The methodology was applied to a coastal aquifer with heterogeneous formations in a semi-arid area near salty Tashk Lake (electrical conductivity 61,420 µS/cm), Fars province, Iran. In the last two decades, the overexploitation of groundwater has caused a major water level drawdown and, consequently, salt-water intrusion. Firstly, flow and transport parameters (hydraulic conductivity, porosity, specific storage coefficient and longitudinal dispersivity) were estimated simultaneously in steady-state and, secondly, in the developed code, these results were used as initial values of the parameters in transient-state. Results show a good match for observed and simulated data. It can be concluded that GA is a helpful tool for automatic calibration of variable density fluid systems such as seawater intrusion cases.     

Laboratory investigation and simulation of breakthrough curves in karst conduits with pools

A series of laboratory experiments are performed under various hydrological conditions to analyze the effect of pools in pipes on breakthrough curves (BTCs). The BTCs are generated after instantaneous injections of NaCl tracer solution. In order to test the feasibility of reproducing the BTCs and obtain transport parameters, three modeling approaches have been applied: the equilibrium model, the linear graphical method and the two-region nonequilibrium model. The investigation results show that pools induce tailing of the BTCs, and the shapes of BTCs depend on pool geometries and hydrological conditions. The simulations reveal that the two-region nonequilibrium model yields the best fits to experimental BTCs because the model can describe the transient storage in pools by the partition coefficient and the mass transfer coefficient. The model parameters indicate that pools produce high dispersion. The increased tailing occurs mainly because the partition coefficient decreases, as the number of pools increases. When comparing the tracer BTCs obtained using the two types of pools with the same size, the more appreciable BTC tails that occur for symmetrical pools likely result mainly from the less intense exchange between the water in the pools and the water in the pipe, because the partition coefficients for the two types of pools are virtually identical. Dispersivity values decrease as flow rates increase; however, the trend in dispersion is not clear. The reduced tailing is attributed to a decrease in immobile water with increasing flow rate. It provides evidence for hydrodynamically controlled tailing effects.     

Probabilistic analysis of groundwater and radionuclide transport model from near surface disposal facilities

As per the regulatory requirements controlling the disposal of radioactive waste, the performance of waste disposal facilities needs to be assessed quantitatively using predictive models. This estimates the potential impact of disposal on the environment and public health. Near Surface Disposal Facilities (NSDFs), constructed to contain the low level radioactive waste are considered to model the radionuclide migration from the system to the geo-sphere. The radiation dose experienced by an individual through drinking water pathway is the endpoint of assessment of the model. A three dimensional groundwater contaminant transport model with a decaying source is modelled numerically to determine the radiation dose for short-lived and long-lived radionuclides. The consideration of uncertainties constitutes an intrinsic part of modelling. The uncertain input parameters include porosity, longitudinal dispersivity, transverse dispersivity, diffusion coefficient and distribution coefficient. The uncertainty propagation and quantification is carried out using collocation based stochastic response surface method (CSRSM). To run the simulations for the huge set of input, a code is developed using built-in python interface in the numerical model. The results are processed further to obtain the sensitive parameters affecting the output concentrations. Further, the probability of radiation dose exceeding permissible value is estimated by subset simulation.     

Magnesium in the marine sedimentary environment: Mg---NH4 ion exchange

The release of exchangeable Mg in marine sediments from displacement by ammonium ions was estimated by way of experimentally determining the parameters that govern this ion-exchange equilibrium on solid geochemical phases: smectite, humic acid, illite and opal.

We showed that: (a) both the conditional selectivity constant as well as the solid concentration are important parameters in determining the relative contribution of ammonium-exchangeable Mg from smectite, organic matter, illite and opal; and (b) that, except in the cases where opal or organic matter concentrations are very high, the clays are the dominant carrier phases for labile Mg which is exchangeable by ammonium.

A model, based on the sum of the contributions from the major geochemical phases present in the sediment reliably predicts the amount of Mg released by exchange with ammonium in marine sediments.     

强含水层的浆液低压运移及其井筒堵水应用

采用液体溶质水动力弥散方程及其相应分析方法描述治理破裂井筒注浆的浆液运移过程,通过坐标变换推导出径向水流条件下考虑吸附作用的浆液浓度解析公式。同时,结合一个煤矿风井堵水工程实例,介绍了如何根据现场观测资料,采用标准曲线法分析流动浆液吸附性参数和指导注浆进程的应用过程。      

岩溶热储井间连通性的示踪研究

岩溶热储储层的不均匀性强,采灌井之间的连通性不易确定。示踪技术可以将运移参数量化,有效刻画储层流体的特征,研究回灌井和开采井之间的水力联系,包括导水通道,流体流速等信息,对长期回灌可能引起的开采井的冷却进行预测,是岩溶热储井间连通性研究十分有效的技术。本文以华北牛驼镇地热田雄县地区为例,针对蓟县系雾迷山组岩溶热储,采用荧光素钠示踪剂,进行了1口井注入,10口井观测的群井示踪试验。采用裂隙介质溶质运移模型,对示踪试验数据进行了解释,得到优势通道的长度、渗透流速、纵向弥散度、回收率等储层性质,获得了试验区内采灌井之间的连通特征。对调整采灌井布局提出了建议。     

Simulation of nitrate transport and wastewater seepage in groundwater flow system

Groundwater quantity and quality modeling is one approach for optimal use of available water resources in arid and semiarid regions. This study was conducted to simulate flow treatment and nitrate transport on Shahrekord aquifer using three-dimensional solute transport model and geographical information system. Hydraulic conductivity, specific yield and recharge values in flow simulation process and effect molecular diffusion coefficient, longitudinal dispersivity and distribution coefficient in quality model were calibrated. 120 water samples during July 2007 to June 2008 were collected monthly from 10 wells and measurements of nitrate were carried out. The results show that the developed model is successfully used to simulate flow path and nitrate transport in saturated porous media. The highest values of nitrate occur along Bahram–Abad village and the surroundings. The groundwater quality in the area represents a complex system, which is affected by different factors of pollution, such as urban wastewater and leaching of agricultural lands.     

Analytical hydrologic models and the design of policy instruments for groundwater-quality management

This paper explores how analytical hydrologic models can inform the effective design and choice of policy instruments to manage groundwater quality by coupling a social-planner’s problem of optimal groundwater-quality management with analytical solutions from the hydrology literature. A theoretical analysis is performed in order to characterize the properties of an optimal emissions policy. The model is then applied in a numerical analysis of groundwater contamination by chloride from highway deicers, demonstrating the relevance of the theoretical results to practical management settings. This analytical approach can help determine which policy instruments are likely to be effective in controlling groundwater pollution, especially if costly numerical groundwater models are not available for the aquifer in question. Unlike previous economic studies of groundwater contamination, this approach defines optimal emissions policies as a function of geophysical parameters employed by hydrologists such as distance between source and sink, groundwater velocity, and aquifer dispersivity. The theoretical section of this paper also demonstrates that the relationship between geophysical parameters and optimal emissions levels may be ambiguous.