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

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

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

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

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

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

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

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

含水层非均质性不同刻画方法对地下水流和溶质运移预测的影响

为探讨含水层非均质性不同刻画方法对地下水流和溶质运移预测的影响,基于非均质含水层砂箱实验,分别用传统等效均质模型、克立金插值和水力层析刻画含水层渗透系数场,并探讨了先验信息对水力层析结果的影响.将不同方法估算的渗透系数场用以预测地下水流和溶质运移过程,以此判断不同方法估算结果的优劣,分析含水层非均质性对地下水流和溶质运移的影响.结果表明:与克立金插值法相比,水力层析法可以更好地刻画含水层非均质性,较准确地预测地下水流和溶质运移过程;钻孔岩心渗透系数样本值作为先验信息可以提高水力层析法估算结果的精度;传统等效均质模型无法准确预测地下水流和溶质运移过程.含水层非均质性的增强将导致溶质污染羽分布形态和运移路径的空间变异性增强,并且优势通道直接决定溶质的分布及运移路径.      

穿透曲线反演管壁溶质通量的实验研究

泉口穿透曲线提供了一个便捷有效的方法来研究岩溶含水层中污染物的运移。在旱季,滞留在石灰岩中的污染物通 过岩壁释放入管道,与管道水相混合并运输到下游泉口,因此,泉口穿透曲线能够反演出污染物释放时的管壁溶质通量。 笔者利用一个基于运移稀释模型的解析解,由40条实验室测量的穿透曲线来反演出相应的管壁溶质通量。结果表明,其基 本形状与穿透曲线相似,管壁透水的非均匀性在穿透曲线上明显表现出来,而在管壁溶质通量上表现为强振荡。因此,管 壁溶质通量比泉口穿透曲线对管壁透水的非均匀性有更高的分辨率。该通量是一个连接泉口穿透曲线和污染物在岩石中早 期分布的中间量,本文提供了一种研究管壁污染物释放过程的新方法。     

注入时间和静水压力对孔隙热储层中Cl-运移影响

迄今为止,注入时间和静水压力对溶质在深层承压地热水中的运移规律影响研究少有报道。通过模拟35℃的低温地热环境,开展了注入时间1,2,3,4,5 h以及静水压力0,6,9 MPa条件下Cl^-的运移柱模拟试验。采用CXTFIT 2.1软件进行数值模拟,探讨了孔隙型热储砂土中Cl^-的运移规律和影响因素。结果表明：在模拟的低温孔隙型热储层中,不同注入时间和静水压力下,Cl^-的运移曲线均呈正态对称分布,一维对流弥散(CDE)模型也可较好地表征其穿透曲线,因此溶质扩散过程符合菲克定律。注入时间的不同,会引起Cl^-的穿透曲线、运移参数发生变化,这与不同注入时间条件下溶质注入总量、柱内溶质浓度差以及分子扩散能力不同有关。在不同静水压力条件下,弥散系数从0 MPa的25.22 cm²/h增加到9MPa的36.13 cm²/h,分子扩散系数、机械弥散系数以及弥散度也随之增大,因此溶质的弥散作用随静水压力的增大而增强。研究结果对于丰富地下水的溶质运移理论具有重要意义。     

一种新的土壤冻融特征曲线模型

土壤冻融特征曲线(SFTC)可以描述冻融过程中未冻水含量随负温的变化关系.准确刻画土壤冻融特征曲线对土壤的冻融过程及相关的水热耦合运移研究至关重要.以往研究中土壤冻融特征曲线适用范围较窄,无法满足实际需要.通过类比水分特征曲线,考虑初始含水量和溶质浓度对未冻水含量的影响,提出了一种新的冻融特征曲线模型,经广泛的文献数据...     

土壤水分与溶质运移机制的分形理论研究进展

土壤中水分和溶质运移一直是土壤-水环境系统中的研究热点,也是目前仍未得到很好解决的问题.将分形理论应用于土壤水分溶质运移机制的研究,探讨其领域中的众多复杂问题,是一种全新的思路和方法.在对土壤结构定量化的分形表征进行简要阐述的基础上,重点介绍了分形理论在水分特征曲线和水力传导度等土壤水分运动基本参数的分形模型、土壤水分运移过程模拟、土壤溶质运移的非费克现象、弥散度的尺度效应以及溶质运移机制研究方面所取得的一系列成果,并就分形理论今后在土壤水分、溶质运移研究中的应用作了展望.     

系统维度对变密度溶质运移的影响研究

关于地下水变密度流与溶质运移的研究通常局限于二维系统中,三维系统中的研究非常有限。然而,仍有零星研究表明,三维系统中的对流过程与二维系统不同。文章通过SEAWAT-2000进行数值模拟,系统研究了二维和三维系统中的变密度溶质运移过程,利用舍伍德数、空间矩和稀释指数量化了系统的不稳定性和溶质扩散、稀释程度。结果表明:二维系统产生分散指流,而三维系统中因扩散程度的增强使得分散指流的产生受到抑制,但是,三维系统的不稳定性比二维系统更强,对流入渗更快,与指流的产生与否并无直接关联。另外,三维系统中溶质的稀释程度大于二维系统,但溶质在二维系统中将更快接近稀释的最大值,传统二阶中心矩可能会造成自由对流和不稳定性溶质运移过程中溶质扩散和稀释程度的错误估计。研究结果将有助于正确预测三维自然含水层中的自由对流和溶质运移。     

Modeling preferential water flow and solute transport in unsaturated soil using the active region model

Preferential flow and solute transport are common processes in the unsaturated soil, in which distributions of soil water content and solute concentrations are often characterized as fractal patterns. An active region model (ARM) was recently proposed to describe the preferential flow and transport patterns. In this study, ARM governing equations were derived to model the preferential soil water flow and solute transport processes. To evaluate the ARM equations, dye infiltration experiments were conducted, in which distributions of soil water content and Cl⁻ concentration were measured. Predicted results using the ARM and the mobile–immobile region model (MIM) were compared with the measured distributions of soil water content and Cl⁻ concentration. Although both the ARM and the MIM are two-region models, they are fundamentally different in terms of treatments of the flow region. The models were evaluated based on the modeling efficiency (ME). The MIM provided relatively poor prediction results of the preferential flow and transport with negative ME values or positive ME values less than 0.4. On the contrary, predicted distributions of soil water content and Cl⁻ concentration using the ARM agreed reasonably well with the experimental data, with ME values higher than 0.8. The results indicated that the ARM successfully captured the macroscopic behavior of preferential flow and solute transport in the unsaturated soil.     

Using TDR and modeling tools to investigate effects of interactive factors on preferential flow and transport in field sandy clay soil

Investigation of the extent of preferential flow and transport affected by several factors and their interactions in the vadose zone using the advanced measurement and modeling techniques is crucial for protection of groundwater from agricultural chemicals like fertilizers and pesticides. The objectives of this study were to investigate the interactive effects of factors like soil structure, initial soil water content (SWC), and application rate on preferential flow and transport using the time domain reflectometry (TDR) measurements of SWC and electrical conductivity (EC) in the plots of 12 treatments in a sandy clay field soil, applying the models (HYDRUS-1D and MACRO) to the measured data, performing Tukey test statistical analysis, and relating model parameters to basic soil properties, consequently, flow and transport characteristics. Analysis of response times and the changes in SWC and EC with time during the experiments in the profiles of the treatments confirmed the existence of preferential flow and transport at the site. As long as the other factors or conditions were constant; undisturbed versus disturbed soil, wet versus dry initial SWC, and high versus low application rate caused preferential flow and transport in the soil based on the TDR measurements. Overall HYDRUS-1D had better performance than MACRO in the simulations of the measured data. These results suggest that different scales of these factors in different field soils need to be further studied for better understanding the flow and transport processes in the vadose zone.     

土壤中优势流的几个基本问题研究

优势流是指土壤在整个入流边界上接受补给,但只通过少部分土体的快速运移,优势流是一种普遍存在的现象,而不是一种特例,它受许多因素的控制,如土壤中的大空隙,土壤结构,土壤质地,土壤水分含量,土壤初始水分含量,水和溶质的施加速率及溶质的施加方法等,优势流的产生机理主要有两种,一种是由土壤介质的非均质所驱动的优势流;另一种是湿润锋的不稳定性所驱动的优势流,目前优势流的监测方法主要取土壤原状土,实验室内的土     

多孔介质中非均匀流动特性的染色示踪试验研究

通过土壤染色剂进行的4组试验,对不同介质结构条件下的水流和溶质非均匀运动规律,非均匀流动变异信息分布特征关系以及全局性非均匀流动示踪方法进行了研究。结果表明,即使在相对比较均匀的介质条件下,流动也表现出明显的非均匀特性;对数正态分布能够较好的反映水流的运动分布模式,相比水流运动,溶质的运动和分布规律明显不同,表现出更多的不确定性和变异性。     

Effects of entrapped residual air bubbles on tracer transport in heterogeneous soil: Magnetic resonance imaging study

Magnetic resonance imaging (MRI) was used to study the process of infiltration and solute transport in an undisturbed soil sample of coarse sandy loam. The sample was subjected to the recurrent ponded infiltration (RPI) experiment, which was carried out in order to assess the changes in the entrapped air volume and its impact on steady state flow rates and solute breakthrough. The main stages of the first and second experimental RPI runs were monitored using an MRI sequence that follows both water density and magnetic relaxation. In a steady state stage of each experimental run a nickel nitrate pulse was injected in order to visualize the solute breakthrough. Effluent from the sample was collected for chemical analysis and a breakthrough curve of the nickel was constructed. To obtain information about the soil structure and to reveal potential preferential pathways, the soil sample was scanned using computed tomography. The local nickel ion transport breakthrough was evaluated from MR images in a series of local observation points distributed along the selected preferential pathways.The preferential flow instability phenomenon with the emphasis on air bubble formation was shown by detecting a 60% decrease of the steady state infiltration rate. The detailed analyses of MRI measurements at observation points revealed air bubble formation, producing a flow rate decrease accompanied by redirection of nickel ion transport trajectories. By analyzing M0 maps it was found that the volumetric water content decrease was 2.2%.     

填埋场垃圾中优先流特征的染色示踪试验

生活垃圾具有大孔隙特性,垃圾中水分运动常表现出优先流特征。针对垃圾组成复杂、难以切片的特点,设计了一种能方便地得到垃圾染色剖面的试验方法。使用填埋场不同深度处的垃圾样进行染色示踪试验,对染色后的垃圾柱分别进行水平切片和竖向切片。基于数字图像处理技术,研究了垃圾中优先流特征。垃圾中大孔隙多、基质少,染色模式与土中不同,基质部分的染色比例较高。试验所得的染色面积比和染色深度能反映出垃圾中优先流的明显程度,染色面积比与参与优先流的大孔隙比例有关,染色面积比越大,优先流越明显,而染色深度与基质流比例有关,上部垃圾中基质流比例越低,入渗深度越深。不同入渗强度的试验结果表明,入渗强度越大时垃圾中优先流程度越高;不同初始含水率的试验结果表明,初始含水率高的垃圾试样中优先流更明显;不同埋深垃圾试样的试验结果表明,浅层垃圾中更易产生优先流。竖向切片与水平向切片所得规律完全一致,试验方法在实现垃圾染色示踪和切片研究的同时,未改变垃圾中优先流原有规律,准确揭示了垃圾中优先流特征。     

长江三峡花岗岩林地坡面降雨渗流与水土流失关系研究

在大量的野外观测试验的基础上,根据土壤水分的形态和能态,对长江三峡花岗岩林地坡面降雨渗流机制进行了分析研究,研究发现在试验区当降雨强度大于表层土壤入渗率,即地表产生积水的情况下,长江三峡花岗岩林地坡面优先流在土壤含水量未达到最大毛管持水量时开始产流,早于土壤渗流产流的时间;在降雨强度小于表层土壤入渗率,即地表不产生积水的情况下,土壤渗流产流早于优先流产流;从壤中流产流过程来看,在地表不产生积水的条件下,土壤渗流先于优先流产生;在地表产生积水而土壤水分含量未达到最大毛管持水量的条件下,优先流先于土壤渗流产生。在两种情况下,优先流的峰值都先于土壤渗流。     

Upscaling of Navier–Stokes equations in porous media: Theoretical, numerical and experimental approach

The accurate estimation of hydraulic conductivity is important for many geotechnical engineering applications, as the presence of fluids affects all aspects of soil behaviour, including its strength. Darcy’s law is the key experimental (or phenomenological) equation employed to model ground water flow. Yet, this phenomenological equation can be linked to a more fundamental microscale model of flow through the pore spaces of the porous material. This paper provides an experimental verification of the relationships between Darcy’s law (macroscale) and the Navier–Stokes equations (microscale) for actual complex pore geometries of a granular material. The pore geometries are experimentally obtained through state-of-the-art X-ray computer assisted micro-tomography. From the numerical modelling of the microscale flow based on actual pore geometries, it is possible to quantify and visualize the development of pore-scale fluid preferential flow-paths through the porous material, and to assess the importance of pore connectivity in soil transport properties.     

降水入渗补给过程中优先流的确定

优先流是降水、灌溉水等入渗补给地下水的主要形式之一, 流速快, 流动路径复杂, 难以定量描述.针对优先流难以定量描述的问题, 以郑州地中渗透仪观测资料为基础, 探讨了新乡亚砂土等试筒降水入渗过程及其中的优先流补给量比例.根据土壤的水力性质、气候等资料建立不存在优先流的数值模拟模型来刻画降水入渗补给过程, 通过模拟获得的地下水入渗补给量与实测地下水入渗补给量的历时曲线, 将大于模拟值的实测值视为优先流的量及确定其在总补给量中所占的比例.结果表明, 优先流占总补给量的比例在10%~80%之间; 随着土壤粘性增加, 优先流所占比例呈增加趋势; 随地下水位埋深的增大, 优先流所占比例呈逐渐下降趋势.      

Rapid runoff via shallow throughflow and deeper preferential flow in a boreal catchment underlain by frozen silt (Alaska, USA)

In high-latitude catchments where permafrost is present, runoff dynamics are complicated by seasonal active-layer thaw, which may cause a change in the dominant flowpaths as water increasingly contacts mineral soils of low hydraulic conductivity. A 2-year study, conducted in an upland catchment in Alaska (USA) underlain by frozen, well-sorted eolian silt, examined changes in infiltration and runoff with thaw. It was hypothesized that rapid runoff would be maintained by flow through shallow soils during the early summer and deeper preferential flow later in the summer. Seasonal changes in soil moisture, infiltration, and runoff magnitude, location, and chemistry suggest that transport is rapid, even when soils are thawed to their maximum extent. Between June and September, a shift occurred in the location of runoff, consistent with subsurface preferential flow in steep and wet areas. Uranium isotopes suggest that late summer runoff erodes permafrost, indicating that substantial rapid flow may occur along the frozen boundary. Together, throughflow and deep preferential flow may limit upland boreal catchment water and solute storage, and subsequently biogeochemical cycling on seasonal to annual timescales. Deep preferential flow may be important for stream incision, network drainage development, and the release of ancient carbon to ecosystems.