共查询到20条相似文献，搜索用时 453 毫秒
1.
《地学前缘(英文版)》2016,7(6)
A new method, the characteristic finite element method(CFEM), was developed to simulate solute transport in a cross-fracture. The solution of this mathematical model for solute transport considered that the contribution of convection and dispersion terms was deduced using the single-step, trace-back method and routine finite element method(FEM). Also, experimental models were designed to verify the reliability and validity of the CFEM. Results showed that experimental data from a single fracture model agreed with numerical simulations obtained from the use of the CFEM. However, routine FEM caused numerical oscillation and dispersion during the calculation of solute concentration. Furthermore, in this cross-fracture model, CFEM simulation results predicted that the arrival time of concentration peak values decreased with increasing flux. Also, the second concentration peak value was obvious with the decrease of flux, which may have resulted from the convergence of solute concentrations from main, and branch, fractures. 相似文献
2.
The hydrodynamic dispersion is an important factor influencing the reactive solute transport in the porous media, and many previous studies assumed that it linearly varied with the average velocity of the groundwater flow. Actually, such linear relationship has been challenged by more and more experimental observations, even in homogeneous media. In this study, we aim to investigate the relationship between hydrodynamics dispersion and the flow velocity in different types of porous media through a laboratory-controlled experiment. The results indicate that (1) the dispersion coefficient should not be a linear function of the flow velocity when the relationship between the flow velocity and the hydraulic gradient can be described by Darcy’s law satisfactorily; (2) Power function works well in describing the dispersion coefficient changing with the flow velocity for different types of porous media, and the power value is between 1.0-2.0 for different particle sizes. 相似文献
3.
溶质暂态存储是岩溶地下水溶质运移过程中的普遍现象。为揭示岩溶管道与裂隙介质间溶质暂态存储机制,本文构建室内管道-裂隙物理模型,开展集中补给条件下的定量示踪试验,运用双区对流弥散模型实现溶质运移过程模拟。研究表明:随着集中补给水动力条件的增强,裂隙暂态存储水量呈线性增加趋势,溶质穿透曲线由单峰型向双峰型转变;管道和裂隙中的平均流速呈负相关关系,溶质在管道和裂隙中的滞留时间差决定了穿透曲线的形态;溶质暂态存储引发了穿透曲线的拖尾效应和双峰现象,对岩溶地下水溶质运移过程具有重要的控制作用。 相似文献
4.
对交叉裂隙渗流传质特性的定量描述是研究整个裂隙网络渗透传质特性的基础。为真实模拟水流及溶质在三维交叉裂隙中的运移过程,首先通过三维轮廓仪获取天然岩石裂隙表面的形貌数据,再应用三维重构技术生成相应的三维交叉裂隙模型,随后求解Navier-Stokes方程,假定溶质运移满足Fick定律,模拟水流和溶质在三维交叉裂隙中的运移过程。通过对比粗糙裂隙模型与平行平板模型的模拟结果发现：粗糙度对流体的分布及流动状态存在显著的影响;不同进、出口工况下的流体流动及溶质运移状态亦表明：裂隙交叉的几何形貌会显著地影响溶质混合行为。这些结果表明,目前被广泛采用的平行平板模型在评估岩体内特别是交叉口的物质运移特性时将导致较大的偏差,在将来的研究中有必要针对裂隙交叉口的几何特征建立修正的模型以提高评估的准确性。 相似文献
5.
径向弥散是指溶质在径向流场下的迁移规律,被广泛用于描述含水层修复领域中污染物的迁移过程.然而,在现有描述径向弥散的模型中,往往忽略了井内混合效应对溶质径向弥散的影响.建立新的注入井附近溶质径向运移动力学模型,同时考虑井内混合效应与弥散度的尺度效应.采用Laplace变换推导该模型的半解析解,利用Stehfest数值逆变换获取溶质在实数空间的解.通过与不考虑混合效应的模型对比研究混合效应对溶质径向弥散的影响,并利用室内渗流槽中的溶质径向弥散实验数据验证模型的合理性与适用性.结果表明:混合效应和尺度效应对注水井附近溶质径向弥散有显著影响.具体地讲,井内的混合效应越显著,在井壁处及含水层中的穿透曲线越低,溶质浓度达到峰值所需时间越长,与不考虑混合效应模型的差异越明显;随尺度效应的增强,溶质提前穿透且扩散范围变大,溶质浓度达到峰值所需时间越长;与前人的模型相比,本研究模型能更好地模拟注水井附近的溶质径向弥散问题. 相似文献
6.
垃圾填埋场滤出液、入侵海水、核废物及生产生活废水等污染物随地下水的迁移威胁人类生存。地下水渗流的随机性导致溶质运移问题更加复杂。根据流网特点,利用流线与水头等势线对求解域进行离散,基于质量守恒原理建立流网单元内溶质浓度求解的隐式有限体积差分格式。基于多孔介质孔隙流速分布规律,利用蒙特卡洛法建立流管单元随机流速场进行溶质运移过程的数值模拟,最后根据数值模拟和模型试验的结果对比,验证了数值模拟方法的准确性。基于流网单元的数值模型中沿流线方向的物质交换由对流和扩散共同作用,而流管间的物质交换只有扩散作用,因此,可在不使用弥散系数下进行污染物运移的模拟。引入随机方法确定流管内流速为研究非均匀流场中染污物的优势迁移提供了新的思路。 相似文献
7.
以基岩渗流方向与裂隙轴向呈45度角为例,探讨了当裂隙轴向与基岩水流斜交时,裂隙流的非达西程度对流场及溶质运移的影响.使用Comsol Multiphysics多物理场仿真软件构建了一个在中部包含水平单裂隙的正方形多孔介质模型,裂隙中的非达西流用Izbash方程去刻画,裂隙的上游基岩中存在持续的溶质源.随着裂隙水流的非达西程度逐渐增强,流场及污染物分布表现出如下特征:（1）裂隙中水流流速逐渐增大;流线在裂隙与基岩界面处的折射逐渐偏离折射定律;（2）裂隙水流的流向逐渐偏向基岩水流的渗流方向;（3）溶质羽宽度变宽但对称性逐渐降低;（4）溶质在水平方向上的浓度峰值逐渐降低,右侧浓度逐渐升高;（5）裂隙产生的回弥散对溶质运移作用逐渐增强,使裂隙中更多的溶质运移到了上层基岩中.总体而言,裂隙流的非达西程度对流场及污染物分布有着显著的影响. 相似文献
8.
迄今为止,注入时间和静水压力对溶质在深层承压地热水中的运移规律影响研究少有报道。通过模拟35℃的低温地热环境,开展了注入时间1,2,3,4,5 h以及静水压力0,6,9 MPa条件下Cl-的运移柱模拟试验。采用CXTFIT 2.1软件进行数值模拟,探讨了孔隙型热储砂土中Cl-的运移规律和影响因素。结果表明：在模拟的低温孔隙型热储层中,不同注入时间和静水压力下,Cl-的运移曲线均呈正态对称分布,一维对流弥散(CDE)模型也可较好地表征其穿透曲线,因此溶质扩散过程符合菲克定律。注入时间的不同,会引起Cl-的穿透曲线、运移参数发生变化,这与不同注入时间条件下溶质注入总量、柱内溶质浓度差以及分子扩散能力不同有关。在不同静水压力条件下,弥散系数从0 MPa的25.22 cm2/h增加到9MPa的36.13 cm2/h,分子扩散系数、机械弥散系数以及弥散度也随之增大,因此溶质的弥散作用随静水压力的增大而增强。研究结果对于丰富地下水的溶质运移理论具有重要意义。 相似文献
9.
裂隙岩体中非饱和渗流与运移的概念模型及数值模拟 总被引：12，自引：2，他引：12
探讨了裂隙岩体中非饱和地下水渗流与溶质运移的几种概念模型的构造及数值模拟问题 ,如裂隙网络模型、连续体模型、等效连续体模型、双孔隙度 (单渗透率 )模型、双渗透率模型、多组份连续体模型等。在裂隙岩体中 ,非饱和地下水的渗流可能只局限于岩体中的岩石组份、或裂隙网络 ,也可能在裂隙和岩石中同时发生 ;对前一种情形只需考虑单一连续体中的流动 ,而后一种情况则需要包括地下水在岩石和裂隙之间的交换。岩体中的裂隙网络往往是溶质运移的主要通道 ;但当溶质在裂隙与岩石之间的渗透和扩散是重要的运移机制时 ,就需要考虑岩石与裂隙界面处的溶质交换。为了模拟岩石与裂隙之间地下水和溶质的交换 ,就需要了解岩石与裂隙之间相互作用的模式和范围 ,使得这类问题的概念模型较单一连续体模型多了一层不确定性、其数值模拟也变得更为困难。因为在实际问题中不易、甚至根本不能判别非饱和渗流的实际形态 ,具体采用哪种模型主要取决于分析的目的和对现场数据的掌握程度。不论哪种模型都会受到模型及参数不确定性的影响 ,因此必须考虑与其他辅助模型的比较. 相似文献
10.
岩溶区地下水数值模拟研究进展 总被引：2，自引：2，他引：0
岩溶含水介质的不均一性导致岩溶地下水流动、溶质运移和热量迁移的数学模拟研究成为地下水模拟的难点。本文综述了岩溶区地下水流模拟的几种方法,重点阐述了等效多孔介质法、双重连续介质法和三重介质法的定义、发展过程和适用范围,并回顾了这几种方法的研究成果。从等效多孔介质法到三重介质法,模拟精度不断提高,适用范围也逐渐由大区域实际问题向小区域理论研究过渡。介绍了溶质运移模拟和热迁移模拟的研究方法及实例。溶质运移模拟以对流弥散方程为基础,其中尺度效应是溶质运移模拟的重点研究问题;热量迁移模拟应考虑地下热水密度变化对地下热水运动的影响。溶质运移模拟和热量迁移模拟往往是将迁移模型和已经调试成功的地下水流动模型相耦合,从而达到模拟溶质及热量迁移的目的。由于溶质运移和热量迁移的复杂性,现阶段水流模型多数处于等效多孔介质模型阶段。综合理论及实际应用,指出精确刻画裂隙及管道和注重基础数学算法是岩溶水数值模拟进步的关键。 相似文献
11.
E. Vetrimurugan M. Senthilkumar L. Elango 《International Journal of Environmental Science and Technology》2017,14(5):1055-1070
Groundwater modelling is an important management tool to study the behaviour of aquifer system under various hydrological stresses. Present study was carried out in deltaic regions of the Cauvery river, with an objective of estimating the minimum river flow required to improve the groundwater quality by numerical modelling. Cauvery river delta is the most productive agricultural plains of south India, but the agricultural activities during the last few decades have decreased due to limited flow in the river and increasing concentration of solutes in groundwater in the eastern parts. In order to understand the causes for increasing concentration, a three-layered finite-difference flow model was formulated to simulate the groundwater head and solute transport. The model was used to simulate the groundwater flow and solute transport for 5 years from July 2007 to June 2012. There was a fairly good agreement between the computed and observed groundwater heads. The chloride and nitrate ions were considered for solute transport modelling. Observed and simulated temporal variation in chloride and nitrate concentrations were comparable. The simulated solute concentrations from July 2007 to June 2012 showed an accumulation of solutes in groundwater of coastal part of the study area. The model was used to find the flow to be maintained in the river and rainfall recharge required to flush the ions into the sea. This can be achieved by maintaining minimal flow in the river and through regulation of fertilizer use as well as by creating awareness of sustainable use of groundwater in this area. 相似文献
12.
地下水污染问题日益严重,研究溶质运移的弥散理论开始应用于实际问题。建立地下水溶质运移模型,对地下水中污染物的运移及发展趋势进行准确预测,是对地下水进行保护、对地下水污染进行控制的基础。而弥散参数的确定则是地下水溶质运移模型建立的关键环节之一,直接影响着模型预测结果的精度和准确性。 对西宁市贵德县地下水污染的水质运移规律进行分析,在贵德县河滨公园林场采用径向收敛流水动力弥散理论方法进行了第四系含水层现场弥散试验,计算了试验场地潜水含水层的弥散度,获得纵向弥散度(aL)为0.843~0.998 cm,横向弥散度(aT)经验推断值为0.17~0.20 cm,为进一步建立该地区的地下水溶质运移模型、预测地下水污染的发展趋势和评价该地区地下水环境质量提供了数据参考。 相似文献
13.
通过对焦作地区浅层地下水中铬(Ⅵ)污染物分布特征进行调查,分析了研究区浅层地下水中铬(Ⅵ)的污染机理,并运用Visual MODFLOW建立地下水流模型及溶质运移模型,模拟预测了浅层地下水中铬(Ⅵ)的迁移规律。结果表明:研究区浅层地下水铬(Ⅵ)污染严重,污染源是位于老君庙西南方向的焦作某电厂堆灰场,主要原因是露天堆放的粉煤灰中的铬(Ⅵ)污染物在长期淋滤作用下下渗污染含水层。气候条件、包气带岩性、地下水化学环境以及人为因素等也间接使浅层地下水铬(Ⅵ)浓度升高;模拟结果显示在未来的五年时间内,受地形和地下水动力场的影响,浅层地下水中铬(Ⅵ)的迁移方向与地下径流方向一致,沿大沙河水流方向上扩散速度更快,污染区域面积逐渐增大。 相似文献
14.
Sensitivity Analysis of a Combined Groundwater Flow and Solute Transport Model Using Local-Grid Refinement: A Case Study 总被引：1，自引：0，他引：1
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. 相似文献
15.
Ali Salim Joodi Stanislas Sizaret Stéphane Binet Ary Bruand Patrick Alberic Michel Lepiller 《Hydrogeology Journal》2010,18(2):295-309
Darcy’s law is the equation of reference widely used to model aquifer flows. However, its use to model karstic aquifers functioning with large pores is problematic. The physics occurring within the karstic conduits requires the use of a more representative macroscopic equation. A hydrodynamic model is presented which is adapted to the karstic aquifer of the Val d’Orléans (France) using two flow equations: (1) Darcy’s law, used to describe water flow within the massive limestone, and (2) the Brinkman equation, used to model water flow within the conduits. The flow equations coupled with the transport equation allow the prediction of the karst transfer properties. The model was tested by using six dye tracer tests and compared to a model that uses Darcy’s law to describe the flow in karstic conduits. The simulations show that the conduit permeability ranges from 5?×?10?6 to 5.5?×?10?5?m2 and the limestone permeability ranges from 8?×?10?11 to 6?×?10?10?m2. The dispersivity coefficient ranges from 23 to 53 m in the conduits and from 1 to 5 m in the limestone. The results of the simulations carried out using Darcy’s law in the conduits show that the dispersion towards the fractures is underestimated. 相似文献
16.
17.
Groundwater flow is an important control on subsurface evaporite (salt) dissolution. Salt dissolution can drive faulting and associated subsidence on the land surface and increase salinity in groundwater. This study aims to understand the groundwater flow system of Gypsum Canyon watershed in the Paradox Basin, Utah, USA, and whether or not groundwater-driven dissolution affects surface deformation. The work characterizes the groundwater flow and solute transport systems of the watershed using a three-dimensional (3D) finite element flow and transport model, SUTRA. Spring samples were analyzed for stable isotopes of water and total dissolved solids. Spring water and hydraulic conductivity data provide constraints for model parameters. Model results indicate that regional groundwater flow is to the northwest towards the Colorado River, and shallow flow systems are influenced by topography. The low permeability obtained from laboratory tests is inconsistent with field observed discharges, supporting the notion that fracture permeability plays a significant role in controlling groundwater flow. Model output implies that groundwater-driven dissolution is small on average, and cannot account for volume changes in the evaporite deposits that could cause surface deformation, but it is speculated that dissolution may be highly localized and/or weaken evaporite deposits, and could lead to surface deformation over time. 相似文献
18.
基于CFP的岩溶管道流溶质运移数值模拟研究 总被引：2，自引：0，他引：2
多重岩溶含水介质的复杂性导致岩溶地下水流动及溶质运移的数学模拟成为地下水研究难点之一。为了探讨岩溶多重含水介质中地下水流溶质运移特征,文章构建了管道流CFP水流模型和MT3DMS溶质运移三维耦合数值模型。在阐述管道流CFP和MT3DMS基本原理的基础上,通过建立水文地质概念模型算例（1个落水洞、4个直管道）,探讨岩溶管道水流及溶质运移规律,分析讨论不同水文地质参数对浓度穿透曲线的影响。研究结果表明:管道流CFP模型能够刻画岩溶管道与基岩裂隙水流交换特征,MT3DMS模型能够模拟穿透曲线的拖尾现象,符合实际岩溶区特征。随着水力梯度、管道直径及管道渗透系数增大,孔隙度减小,浓度曲线峰值越大,峰值到达时间越快,浓度穿透曲线越对称。得出结论:耦合CFP水流模型和MT3DMS溶质运移模型能够刻画岩溶管道流溶质运移规律,为研究岩溶复杂介质污染物运移特征提供一种思路和途径。 相似文献
19.
为探讨含水层非均质性不同刻画方法对地下水流和溶质运移预测的影响,基于非均质含水层砂箱实验,分别用传统等效均质模型、克立金插值和水力层析刻画含水层渗透系数场,并探讨了先验信息对水力层析结果的影响.将不同方法估算的渗透系数场用以预测地下水流和溶质运移过程,以此判断不同方法估算结果的优劣,分析含水层非均质性对地下水流和溶质运移的影响.结果表明:与克立金插值法相比,水力层析法可以更好地刻画含水层非均质性,较准确地预测地下水流和溶质运移过程;钻孔岩心渗透系数样本值作为先验信息可以提高水力层析法估算结果的精度;传统等效均质模型无法准确预测地下水流和溶质运移过程.含水层非均质性的增强将导致溶质污染羽分布形态和运移路径的空间变异性增强,并且优势通道直接决定溶质的分布及运移路径. 相似文献
20.
Groundwater aquifer vulnerability has been assessed by incorporating the major geological and hydrogeological factors that affect and control the groundwater contamination using GIS-based DRASTIC model along with solute transport modeling. This work demonstrates the potential of GIS to derive a vulnerability map by overlying various spatially referenced digital data layers (i.e., depth to water, net recharge, aquifer media, soil media, topography, the impact of vadose zone and hydraulic conductivity) that portrays cumulative aquifer sensitivity ratings in Kishangarh, Rajasthan. It provides a relative indication of groundwater aquifer vulnerability to contamination. The soil moisture flow and solute transport regimes of the vadose zone associated with specific hydrogeological conditions play a crucial role in pollution risk assessment of the underlying groundwater resources. An effort has been made to map the vulnerability of shallow groundwater to surface pollutants of thestudy area, using soil moisture flow and contaminant transport modeling. The classical advection-dispersion equation coupled with Richard’s equation is numerically simulated at different point locations for assessing the intrinsic vulnerability of the valley. The role of soil type, slope, and the land-use cover is considered for estimating the transient flux at the top boundary from daily precipitation and evapotranspiration data of the study area. The time required by the solute peak to travel from the surface to the groundwater table at the bottom of the soil profile is considered as an indicator of avulnerability index. Results show a high vulnerability in the southern region, whereas low vulnerability is observed in the northeast and northern parts. The results have recognized four aquifer vulnerability zones based on DRASTIC vulnerability index (DVI), which ranged from 45 to 178. It has been deduced that approximately 18, 25, 34, and 23% of the area lies in negligible, low, medium and high vulnerability zones, respectively. The study may assist in decision making related to theplanning of industrial locations and the sustainable water resources development of the selected semi-arid area. 相似文献