基于CFP的岩溶管道流溶质运移数值模拟研究

多重岩溶含水介质的复杂性导致岩溶地下水流动及溶质运移的数学模拟成为地下水研究难点之一。为了探讨岩溶多重含水介质中地下水流溶质运移特征,文章构建了管道流CFP水流模型和MT3DMS溶质运移三维耦合数值模型。在阐述管道流CFP和MT3DMS基本原理的基础上,通过建立水文地质概念模型算例（1个落水洞、4个直管道）,探讨岩溶管道水流及溶质运移规律,分析讨论不同水文地质参数对浓度穿透曲线的影响。研究结果表明:管道流CFP模型能够刻画岩溶管道与基岩裂隙水流交换特征,MT3DMS模型能够模拟穿透曲线的拖尾现象,符合实际岩溶区特征。随着水力梯度、管道直径及管道渗透系数增大,孔隙度减小,浓度曲线峰值越大,峰值到达时间越快,浓度穿透曲线越对称。得出结论:耦合CFP水流模型和MT3DMS溶质运移模型能够刻画岩溶管道流溶质运移规律,为研究岩溶复杂介质污染物运移特征提供一种思路和途径。     

弯曲岩溶管道溶质运移的尺度效应研究

岩溶管道溶质运移的尺度效应研究对穿透曲线的正确解译非常重要,但目前针对单一弯曲管道中溶质运移尺度效应的研究仍比较缺乏。文章将岩溶管道和溶潭分别概化为透明软管和水箱,基于前期建立的水箱-管道系统(简称“管道系统”),在水箱下游设置不同长度的弯曲管道,通过示踪试验研究管道运移尺度对穿透曲线的影响,并采用暂时存储模型模拟试验曲线。结果表明：(1)随着水箱下游管道长度的增加,峰值质量浓度逐渐缓慢降低,穿透曲线上升段斜率无明显变化,穿透曲线拖尾逐渐缩短,表明运移管道长度增加对溶质运移的影响大于下游管道弯曲;(2)穿透曲线偏度系数、后段溶质羽穿透时间和溶质羽穿透时间与管道系统长度呈良好的负线性相关关系(R²>0.96);(3)当对称和不对称水箱管道系统长度分别增加至154.5 m和164.3 m时,偏度系数接近0,穿透曲线分布接近对称;(4)弥散系数、存储区截面积和交换系数与管道系统长度呈良好的负线性相关关系,当对称和不对称水箱管道系统长度分别增加至159.9 m和178.1 m时,存储区截面积接近0,水箱导致的溶质运移滞后效应基本消失。研究结果对野外岩溶管道穿透曲线的...     

岩溶管道与裂隙介质间溶质暂态存储机制

溶质暂态存储是岩溶地下水溶质运移过程中的普遍现象。为揭示岩溶管道与裂隙介质间溶质暂态存储机制,本文构建室内管道-裂隙物理模型,开展集中补给条件下的定量示踪试验,运用双区对流弥散模型实现溶质运移过程模拟。研究表明：随着集中补给水动力条件的增强,裂隙暂态存储水量呈线性增加趋势,溶质穿透曲线由单峰型向双峰型转变;管道和裂隙中的平均流速呈负相关关系,溶质在管道和裂隙中的滞留时间差决定了穿透曲线的形态;溶质暂态存储引发了穿透曲线的拖尾效应和双峰现象,对岩溶地下水溶质运移过程具有重要的控制作用。     

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

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

多孔介质粒径变异系数对污染物运移的影响

多孔介质污染物运移对于明晰地下水污染很重要,但在多孔介质中粒径变异系数(coefficientofvariation,COV)对内部微观孔隙结构中污染物运移过程影响的研究还存在不足.为此,基于随机算法,提出了一种不同COV且孔隙度一致的多孔介质几何模型构建方法,通过对Navier-Stokes(N-S)方程和对流-扩散方程(advection-diffusion equation,简称ADE)进行耦合求解得到多孔介质地下水流场及污染物浓度场,引入克里斯琴森均匀系数,定量评价流场流速分布的均匀性,基于MIM(mobile and immobile)模型和ADE模型分析耦合求解得到的穿透曲线特征.结果表明：随着粒径COV的增大,流场流速分布的不均性增强,MIM模型中的溶质流动区域占比β、无量纲传质率α^*均增大;MIM模型拟合优度高于ADE模型,且随COV增大,ADE模型的拟合全局误差Ei逐步增大.总体上,粒径COV控制了溶质流动区域和非流动区域的大小及其之间的溶质交换强度,造成了多孔介质内部溶质运移的“非费克”行为,使得ADE模型的误差逐步增大,对于较大COV的多孔介...     

土石混合介质中非反应性阴离子运移试验研究

为探讨化学物质在土石混合介质中的运移过程和机理,采用饱和稳态流下的Cl-1混合置换试验,测定水流和溶质运移过程,分析土石混合介质的溶质穿透曲线特征及碎石组成和含量对运移过程的影响。选用CXTFIT2.1的平衡和物理非平衡对流弥散模型,对参数弥散系数D和滞留因子R进行反求。结果显示:不同土石比的D变异较大:0.258~22.31 cm2/h。R的波动范围为0.6~1.54;碎石含量影响土石介质的溶质运移过程表现为平均孔隙流速、弥散系数、弥散度均与土石比成负指数的幂函数关系。对碎石粒径与溶质运移参数进行相关分析发现,小粒径的碎石含量增加,则孔隙流速和弥散系数有减少的趋势,而大于10 mm的碎石有利于溶质的运移。通过土石介质的非反应性阴离子的混合置换试验研究,可以为非均一介质中化学物质运移提供参考。     

非均匀土壤中溶质运移的两区模型及其解析解

给出了非均匀土壤中考虑水动力弥散尺度效应的一维溶质运移两区模型。然后在初始浓度为零,半无限一维空间内常数通量边界条件下,推导出了可动区和不动区溶质浓度分布的解析表达式。并用Mathematica软件包作了数值计算。计算结果表明,该模型的模拟结果符合溶质运移的一般规律,因而该模型在土壤溶质运移的研究和应用中有一定的价值。     

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

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

死端孔隙对溶质运移影响的实验研究

拖尾分布是溶质运移规律研究的热点和难点。死端孔隙、透镜体等空间结构形成的非均质性以及密度变化产生的指状流会导致溶质运移拖尾分布。论文选择保守性溶质氯化钠作为示踪剂,采用砂柱实验系统地研究了以死端孔隙导致的非均质性以及由不同浓度和温度形成的变密度条件对溶质运移拖尾分布的影响。实验结果表明:死端孔隙的存在会延迟溶质运移穿透曲线出峰时间,降低峰值及出现拖尾分布,随着死端孔隙介质体积或数量的增加这一现象越明显;随着溶液浓度的增加拖尾现象越明显,且温度的升高会加剧这一现象。对比死端孔隙、温度和浓度三者对溶质运移拖尾分布的影响,其中影响最大的是死端孔隙,其次是浓度和温度。     

灰坝工程领域岩体裂隙渗流浸润线及出漏点溶质浓度数值模拟

灰坝工程中灰水渗漏易对地下水产生污染,为预防治理地下水污染,应研究灰场污染物的运移机理,在运移工程中由于岩体裂隙的粗糙性和各向异性,渗流及溶质运移模型实际非常复杂。本文着眼于建立将粗糙裂隙岩体非连续渗流、溶质运移数学模型,应用于陡立裂隙网络岩体渗流浸润线及出漏点溶质浓度数值模拟,为定量分析灰场对地下水的影响提供理论依据。     

Fluid and Solute Transport from a Conduit to the Matrix in a Carbonate Aquifer System

Within carbonate systems, the working hypothesis suggests that when a conduit is flooded fluid and solute migrate from the conduit into the matrix. This flux of fluid and solute into the matrix creates a reservoir that can be slowly released once the flooding recedes. Although hypothesized, these fluxes have never been measured. To quantify the distance that a fluid and solute would move into a matrix, the fluxes of fluid and solute from a conduit into a matrix were simulated for nine different carbonate aquifer systems. Two independent numerical approaches were used to simulate (1) fluid flux into the matrix and (2) solute flux into the matrix during a flooding event. When flooding occurs within the conduit, the volume of water transported into and stored in the matrix with a high porosity and high hydraulic conductivity (Floridan Aquifer) was less than 0.34 m³ along a 1 m length of conduit, resulting in a penetration depth of 7.2×10⁻² m into the matrix. In a low porosity and low hydraulic conductivity matrix (Ozark Plateau), the volume of water transported into and stored in the matrix was less than 6.85×10⁻⁵ m³ along a 1 m length of conduit, resulting in a penetration depth of 2.0×10⁻⁴ m into the matrix. Simulated solute flow shows that less than 0.1% of the solute moves in to the matrix. The two approaches demonstrate that during high flow conditions fluid and solute are forced through the conduits, with very little moving into the carbonate matrix. Once the fluid and solute enter a conduit and are moving, they will remain in the conduit until they are discharged at an outlet. Thus, a carbonate matrix does not become a reservoir for solute and fluid during a high-flow event.     

Effect of source integration on the geochemical fluxes from springs

Geochemical fluxes from watersheds are typically defined using mass-balance methods that essentially lump all weathering processes operative in a watershed into a single flux of solute mass measured in streamflow at the watershed outlet. However, it is important that we understand how weathering processes in different hydrological zones of a watershed (i.e., surface, unsaturated, and saturated zones) contribute to the total geochemical flux from the watershed. This capability will improve understanding of how geochemical fluxes from these different zones may change in response to climate change. Here, the geochemical flux from weathering processes occurring solely in the saturated zone is investigated. This task, however, remains exceedingly difficult due to the sparsity of subsurface sampling points, especially in large, remote, and/or undeveloped watersheds. In such cases, springflow is often assumed to be a proxy for groundwater (defined as water residing in fully saturated geologic formations). However, springflow generation may integrate different sources of water including, but not limited to, groundwater. The authors’ hypothesis is that long-term estimates of geochemical fluxes from groundwater using springflow proxies will be too large due to the integrative nature of springflow generation. Two conceptual models of springflow generation are tested using endmember mixing analyses (EMMA) on observations of spring chemistries and stable isotopic compositions in a large alpine watershed in the San Juan Mountains of southwestern Colorado. In the “total springflow” conceptual model, springflow is assumed to be 100% groundwater. In the “fractional springflow” conceptual model, springflow is assumed to be an integration of different sources of water (e.g., groundwater, unsaturated flow, preferential flow in the soil, etc.) and groundwater is only a fractional component. The results indicate that groundwater contributions in springflow range from 2% to 100% overall and no springs are consistently composed of 100% groundwater; providing support for the fractional springflow conceptual model. Groundwater contributions are not strongly correlated with elevation, spring contributing area, spring discharge, or seasonality. This variability has a profound effect on long-term geochemical fluxes. The geochemical fluxes for total springflow overestimate long-term solute release by 22–48% as compared to fractional springflow. These findings illustrate that springflow generation, like streamflow generation, integrates many different sources of water reflecting solute concentrations obtained along many different geochemical weathering pathways. These data suggest that springs are not always ideal proxies for groundwater. Springs may be integrating very distinct portions of the groundwater flow field and these groundwater contributions may become mixed at the spring emergence with much younger sources of water that have never resided in the groundwater system.     

Geological structure, recharge processes and underground drainage of a glacierised karst aquifer system, Tsanfleuron-Sanetsch, Swiss Alps

The relationships between stratigraphic and tectonic setting, recharge processes and underground drainage of the glacierised karst aquifer system ‘Tsanfleuron-Sanetsch’ in the Swiss Alps have been studied by means of various methods, particularly tracer tests (19 injections). The area belongs to the Helvetic nappes and consists of Jurassic to Palaeogene sedimentary rocks. Strata are folded and form a regional anticlinorium. Cretaceous Urgonian limestone constitutes the main karst aquifer, overlain by a retreating glacier in its upper part. Polished limestone surfaces are exposed between the glacier front and the end moraine of 1855/1860 (Little Ice Age); typical alpine karrenfields can be observed further below. Results show that (1) large parts of the area are drained by the Glarey spring, which is used as a drinking water source, while marginal parts belong to the catchments of other springs; (2) groundwater flow towards the Glarey spring occurs in the main aquifer, parallel to stratification, while flow towards another spring crosses the entire stratigraphic sequence, consisting of about 800 m of marl and limestone, along deep faults that were probably enlarged by mass movements; (3) the variability of glacial meltwater production influences the shape of the tracer breakthrough curves and, consequently, flow and transport in the aquifer.     

Increase of vulnerability of karst aquifers due to leakage from landfills

Karst aquifers are very vulnerable to various pollution sources. Among these, landfills can contribute contaminants over long periods. Groundwater quality monitoring is required to assess the impact of a landfill leachate on underlying aquifer water or spring discharge. Tracer tests are useful for selecting the location and frequency of sampling. Three landfill sites in karst areas were studied with tracer tests. Additional insight was accomplished by the comparison of the tracer transfer and the breakthrough curves obtained from these tests with the results of two other tracer tests carried out in the same period in karst areas without a landfill. To explain the differences observed, the hypothesis of increased permeability below the landfill due to the presence of inorganic acids in leachates was further tested. The results of detailed, long-term monitoring of contaminated drip water in the Postojna Cave were used to verify the hypothesis. The analysis of the simultaneously increased content of calcium, magnesium, and contaminants in the drip water indicates a direct correlation between limestone dissolution and contaminants. Increased dissolution increases secondary porosity and thus permeability of the vadose zone resulting in higher vulnerability of underlying aquifers and springs in the vicinity landfills.     

Application of a solute transport model under variable velocity conditions in a conduit flow aquifer: Olalde karst system, Basque Country, Spain

A model based on numerical solutions, which allows for solving the dispersion equation under variable recharge and velocity conditions, is developed to simulate solute transport in conduit flow aquifers during flow recession periods. As an example, the evolution of a tracer in the little known karst conduit that links the sinking stream of Oma valley to the Olalde spring is investigated in the karstic region of Santa Eufemia-Ereñozar (Basque Country, Spain). The model, with different hypothetical structures, allows for obtaining series of tracer breakthrough curves, which are fitted to experimental data using an optimization algorithm. These results, although they can be used to simulate the tracer evolution between the two points considered, do not allow for determining the internal structure and spatial disposition of contributions in the aquifer.     

长白山火山区温泉温室气体排放通量研究

温泉是深部岩浆活动在地表的直接表现,并且向大气圈排放大量的温室气体.然而,国内尚无火山区温泉排放的温室气体通量研究报道.我国长白山火山区水热活动强烈,主要有湖滨温泉带、聚龙温泉群、锦江温泉以及火山口外围的十八道沟温泉.本文利用数字皂膜流量计测量温泉气体排放通量,并结合前人对长白山火山区温泉气体成分的研究成果,估算了研究区温泉所排放的温室气体通量.结果表明,长白山火山区温泉排放的CO2通量为6.9×104t·a-1,CH4排放通量为428.44t·a-1,与意大利Pantelleria Island火山区温泉排放的温室气体通量规模相当.本文的测试结果表明:数字皂膜流量计在火山区温室气体排放通量估算研究中的应用是可行的.     

Combined use of natural and artificial tracers to determine the hydrogeological functioning of a karst aquifer: the Villanueva del Rosario system (Andalusia,southern Spain)

Analysis of natural responses of karst springs provides information on the behavior of the aquifers they drain. Detailed monitoring and qualitative and quantitative analyses of natural responses, and environmental—total organic carbon (TOC), NO₃ ^?, Cl^? and intrinsic fluorescence—and artificial (fluorescent dye) tracers, in the water drained by Villanueva del Rosario spring (southern Spain), suggest the existence of a conduit flow system with rapid flows and very short transit times of water through the aquifer. This is in agreement with uranine and eosin breakthrough curves and with simple numerical models done using these data. However, due to the low capacity for natural regulation, not all the recharge effects are simultaneously transmitted to the spring water; given a single input, the system modulates and transfers hydrodynamic variations faster than variations of chemical composition and of water temperature. Additionally, time lags between maximum concentrations of natural and artificial tracers show that the global system response (including diffuse infiltration) is faster and more sensitive than that produced from infiltration concentrated at a single point on the surface (sinkholes).     

山西娘子关泉群及其水的来源

娘子关泉域群泉是中国北方最大的岩溶泉,泉域汇水面积达7 436 km2,前人认为:泉域内岩溶水由北、西、南3面向娘子关一带径流汇集,由于娘子关一带下奥陶系燧石团块或条带白云岩相对隔水层隆起,并被桃河侵蚀出露,使岩溶地下水溢出地表成泉群,其主要含水层为中奥陶系含石膏碳酸盐岩。但各泉的水化学、同位素特征有差异,娘子关泉群并不是出自统一源。文章通过水化学、同位素、水文地质剖面等方法研究得出: 娘子关泉域存在两个含水层、三个子系统:中奥陶系灰岩含水层和中上寒武系白云岩含水层;西部奥陶系岩溶水系统、东部奥陶系岩溶水系统和东部中上寒武系岩溶水系统。泉域内城西泉与程家泉出露于中奥陶系下马家沟组泥灰岩之上,含水层为中奥陶系灰岩裂隙、溶隙水,由于区域下马家沟组泥灰岩隆起隔水出露地表成泉,属于东部奥陶系岩溶水系统;坡底泉、五龙泉、河北泉、水帘洞泉、苇泽关泉其补给主要来源于中上寒武系含水岩组,为承压上升泉,属于东部中上寒武系岩溶水系统。