共查询到20条相似文献，搜索用时 500 毫秒
1.
This paper presents a new numerical tool to model the two-dimensional contaminant transport through saturated porous media
using a meshfree method, called radial point interpolation method (RPIM) with polynomial reproduction. In RPIM, an approximate
solution is constructed entirely in terms of a set of nodes and no characterisation of the interrelationship of the nodes
is needed. The advection–dispersion equation with sorption is considered to illustrate the applicability of the RPIM. The
Galerkin weak form of the governing equation is formulated using 2D meshfree shape functions constructed using thin plate
spline radial basis functions. MATLAB code is developed to obtain the numerical solution. Three numerical examples are presented
and the results are compared with those obtained from the finite element method and analytical solutions. In order to test
the practical applicability and performance of the RPIM, two case studies of contaminant transport through landfill liners
are presented. A good agreement is obtained between the results of the RPIM and the field investigation data. 相似文献
2.
Determination of groundwater solute transport parameters in finite element modelling using tracer injection and withdrawal testing data 总被引：1，自引：0，他引：1 
下载免费PDF全文
下载免费PDF全文 Nguyen Van Hoang 《地下水科学与工程》2021,9(4):292-303
The groundwater tracer injection and withdrawal tests are often carried out for the determination of aquifer solute transport parameters. However, the parameter analyses encounter a great difficulty due to the radial flow nature and the variability of the temporal boundary conditions. An adaptive methodology for the determination of groundwater solute transport parameters using tracer injection and withdrawal test data had been developed and illustrated through an actual case. The methodology includes the treatment of the tracer boundary condition at the tracer injection well, the normalization of tracer concentration, the groundwater solute transport finite element modelling and the method of least squares to optimize the parameters. An application of this methodology was carried out in a field test in the South of Hanoi city. The tested aquifer is Pleistocene aquifer, which is a main aquifer and has been providing domestic water supply to the city since the French time. Effective porosity of 0.31, longitudinal dispersivity of 2.2 m, and hydrodynamic dispersion coefficients from D = 220 m2/d right outside the pumping well screen to D =15.8 m2/d right outside the tracer injection well screen have been obtained for the aquifer at the test site. The minimal sum of squares of the differences between the observed and model normalized tracer concentration is 0.00119, which is corresponding to the average absolute difference between observed and model normalized concentrations of 0.035 5 (while 1 is the worst and 0 is the best fit). 相似文献
3.
One-dimensional unsteady solute transport along unsteady flow through inhomogeneous medium 总被引：1，自引：0，他引：1
SANJAY K YADAV ATUL KUMAR DILIP K JAISWAL NAVEEN KUMAR 《Journal of Earth System Science》2011,120(2):205-213
The one-dimensional linear advection–diffusion equation is solved analytically by using the Laplace integral transform. The
solute transport as well as the flow field is considered to be unsteady, both of independent patterns. The solute dispersion
occurs through an inhomogeneous semi-infinite medium. Hence, velocity is considered to be an increasing function of the space
variable, linearly interpolated in a finite domain in which solute dispersion behaviour is studied. Dispersion is considered
to be proportional to the square of the spatial linear function. Thus, the coefficients of the advection–diffusion equation
are functions of both the independent variables, but the expression for each coefficient is considered in degenerate form.
These coefficients are reduced into constant coefficients with the help of a new space variable, introduced in our earlier
works, and new time variables. The source of the solute is considered to be a stationary uniform point source of pulse type. 相似文献
4.
To better understand the movement and transport of water and pollution through the coarse gravel unsaturated zone, the presented
research was conducted to estimate water flow and transport processes with a tracing experiment in a lysimeter in the Selniška
Dobrava. A combined tracing experiment was performed with deuterated water and the fluorescent dye—uranine. The fastest and
dominant flow velocities were calculated based on injection time, the first tracer appearance time and the time of highest
concentration. Mean flow velocity and vertical dispersion were estimated by an analytical best-fit method using one-dimensional
convection–dispersion model. Deuterium was confirmed as an ideal conservative tracer and a more suitable tracer than dye (uranine)
for the study of water flow in the unsaturated zone of a coarse gravel aquifer. The retardation factor of the dye as compared
with deuterium was 1.13–1.75, which is in agreement with previously published results. Artificial tracers, especially deuterated
water, were also identified as a very useful tool to assess other properties and differences in water flow in the unsaturated
zone of a coarse gravel aquifer such as velocity and dispersion. 相似文献
5.
选择北京平原区水文地质、环境地质等方面都比较典型的沙子营垃圾堆放场,建立了其水文地质模型。在充分收集资料、分析得出部分计算所需参数后,采用现场弥散实验、勘查取样测试等方法,求得了该含水层的弥散系数等参数;用二维非稳定流溶质运移方程对污染物在此含水层中的迁移扩散规律、速度和污染范围等进行了模拟计算;采用现场钻探、取样测试分析等方法,评价了该含水层的实际污染状况。实验模拟计算和现场调查结果表明:污染物在潜水含水层中的运移规律遵循二维非稳定流场中的溶质运移方程,污染物在潜水含水层中的运移速度约为86.25m/a,迁移扩散主要发生在地下水流向上,侧向扩散宽度极小,是地下水流向上的1/17。 相似文献
6.
为确定研究区含水层的弥散系数,进一步开展咸水体运移的数值模拟,根据研究区的地质背景和水文地质条件,针对含水层中咸水体分布格局设计并开展了示踪试验。利用试验获得的数据初步分析了咸水体垂向和横向运移特征,提出了依赖于单尺度示踪试验结果推求多尺度下弥散系数的内推-外延法,据此计算了试验尺度上的弥散系数。示踪试验结果直观地反映了咸水体对抽水井中氯离子质量浓度的影响,确定了试验尺度上的横向弥散度为7.42 m,垂向弥散度为41.4 m。基于弥散方程的解析解,采用半程拟合的方法进行尺度效应的分析,确立了横向及垂向弥散度与试验尺度关系分别为y=0.007 5·e0.129 3x和y=0.102 9e0.049 2x。根据上述定量关系可进行同类地区相似条件下的弥散系数计算和尺度效应分析,从而为进一步开展咸淡水交互地带水盐运移的模拟和预测提供参数基础。 相似文献
7.
Tomás Morales-Juberías Martín Olazar Jose María Arandes Pedro Zafra Iñaki Antigüedad Félix Basauri 《Environmental Geology》1997,30(3-4):143-151
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. 相似文献
8.
弥散度是刻画孔隙介质中溶质运移和扩散的重要参数,对于污染物的预测和修复至关重要,但野外示踪试验往往会选择忽略真实存在的井内混合效应。通过室内砂槽实验方法,模拟具有水平分层结构的含水层,该含水层主要由3种介质充填而成。采用埋藏传感器和井中布设传感器2种监测方式,对比在有/无混合效应情况下,穿透曲线的形态差异,进而探究井内混合效应对弥散尺度依赖性的影响情况。实验结果表明,井内混合效应会使穿透曲线呈现阶梯式增长,并伴有显著的拖尾现象;当使用对流弥散方程进行计算时,混合效应会导致弥散度被高估;观测到的弥散度与真实弥散度的差异会随着注入井和观测井间距离的增加而增大;此外,2种观测方式（埋藏/井内）均能发现弥散的尺度依赖性,且井内混合效应显著增强了弥散尺度效应,该实验结果可为污染物运移的评价和预测提供参考。 相似文献
9.
This paper describes a natural-gradient field tracer test to characterise solute-transport properties in a sand and gravel
aquifer in the Hebei Province, northern China. Some laboratory-scale column tests on aquifer material and a local-scale field
borehole-dilution test have been conducted previously, but the field test reported herein represents the only large-scale
tracer test in the aquifer, which is the sole water supply to the city of Shi Jiazhuang and which is threatened by urban pollution.
The aim of the study was to quantify the transport behaviour of nonreactive pollutants in this aquifer. Little quantitative
data are available concerning its solute-transport properties; thus, the results of the tracer test are significant and critical
for understanding pollutant transport and fate. The in-situ tracer test was carried out in the aquifer using a slug injection
of the geochemically conservative, radioactive iodine tracer 131I. The longitudinal (α
L
) and transverse (α
T
) hydrodynamic dispersivities for solute transport in the field are 1.72 and 0.0013 m, respectively. The ratio of longitudinal
dispersivity α
L
and the flow length at the field scale is 1:10. The ratio between α
L
and α
T
from the in-situ test (~1,300:1) demonstrates a dominant longitudinal dispersion in this fluvial sand and gravel aquifer.
The tracer test further indicates a relatively short transit time for the aquifer (linear velocities ~13 m/d) under natural-gradient
conditions.
Electronic Publication 相似文献
10.
The transport of chemically reactive solutes (e.g. surfactants, CO2 or dissolved minerals) is of fundamental importance to a wide range of applications in oil and gas reservoirs such as enhanced
oil recovery and mineral scale formation. In this work, we investigate exponential time integrators, in conjunction with an
upwind weighted finite volume discretisation in space, for the efficient and accurate simulation of advection–dispersion processes
including non-linear chemical reactions in highly heterogeneous 3D oil reservoirs. We model sub-grid fluctuations in transport
velocities and uncertainty in the reaction term by writing the advection–dispersion–reaction equation as a stochastic partial
differential equation with multiplicative noise. The exponential integrators are based on the variation of constants solution
and solve the linear system exactly. While this is at the expense of computing the exponential of the stiff matrix representing
the finite volume discretisation, the use of real Léja point or the Krylov subspace technique to approximate the exponential
makes these methods competitive compared to standard finite difference-based time integrators. For the deterministic system,
we investigate two exponential time integrators, the second-order accurate exponential Euler midpoint (EEM) scheme and exponential
time differencing of order one (ETD1). All our numerical examples demonstrate that our methods can compete in terms of efficiency
and accuracy compared with standard first-order semi-implicit time integrators when solving (stochastic) partial differential
equations that model mixing and chemical reactions in 3D heterogeneous porous media. Our results suggest that exponential
time integrators such as the ETD1 and EEM schemes could be applied to typical 3D reservoir models comprising tens to hundreds
of thousands unknowns. 相似文献
11.
山东淄博市大武水源地裂隙岩溶水中污染物运移的数值研究 总被引：11，自引：0，他引：11
在分析研究淄博市大武水源地裂隙岩溶含水层的水力性质和污染物运移特征的基础上 ,对裂隙岩溶水的水头和污染物运移进行数值研究。目前国内外对裂隙岩溶水进行数值计算时 ,通常用等价多孔介质模型 ,但裂隙岩溶介质和多孔介质有很大不同。裂隙岩溶介质的储水和导水空间为裂隙网络 ,导水系数大 ,地下水的实际平均流速比孔隙水大得多 ,但给水度和贮水系数小。当用等价多孔介质模型进行模拟时应考虑这些特点。对于污染物运移的模拟 ,要同时求解水头方程和对流弥散方程 ,可采用MODFLOW和MT3D软件进行模拟。研究区裂隙岩溶水水头的数值计算表明 ,等效多孔介质模型水头的拟合误差能满足国标GB/T144 97- 93的要求。各时段地下水水量均衡计算的精度也满足要求。对流弥散方程的数值计算 ,由于Peclet数高达 95 .6 7,对流占绝对优势 ,可能存在数值弥散和数值振荡 ,因而采用多种方法进行了比较。对于同一问题 ,同时采用上游有限差分法 (UFDM) ,混合的欧拉拉格朗日方法 (特征线法MOC、改进特征线法MMOC和混合特征线法HMOC) ,总变异消减法(TVD)进行计算 ,并比较其结果。结果表明 ,混合特征线法 (HMOC)和总变异消减法 (TVD)比较适合于对流占优势的运移问题计算。由于渗透系数K和有效孔隙度θ对溶质运移结果的影响很大 ,? 相似文献
12.
Papers presented at a two-day jointly sponsored IAHS/AGU symposium on groundwater contamination are briefly summarized. This international symposium was held 11–12 May, 1989, in Baltimore, Maryland. Presentations encompassed recent research developments in three general areas: abiotic and biotic processes governing contaminant transport; aquifer rehabilitation; and the influence of agricultural practices and nonpoint sources on aquifer quality. Contributions offered an interesting mixture of theoretical, mathematical, laboratory, and field studies. In the first session, transport processes explored ranged from dispersion and fingering to nonequilibrium sorption, metals complexation, and bacteria migration. The use of optimization modeling in the design of remediation strategies was the focus of another session. Here theoretical studies were presented alongside case histories of aquifer rehabilitation. In a final session, a number of models for agricultural management were described. These presentations were complemented by case studies of actual aquifer degradation resulting from land-use and management practices. 相似文献
13.
定量研究污染物在地下水中的运移过程通常采用数值模拟方法。MT3D是一套基于有限差分方法的污染物运移模拟软件，近年来在国外水文地质和水环境模拟等领域的研究中已经得到较为广泛的认可。MT3D比较全面地考虑了污染物在地下水中的对流、弥散和化学反应等过程，可以灵活处理各种复杂的源汇项和边界条件，能够准确模拟承压、无压和越流含水层中的污染物运移过程。MT3D具有模块化的程序结构、灵活的求解方法以及全面的模拟功能，非常适合实际问题的研究，值得在国内推广使用。 相似文献
14.
Aidan Foley Gerd Cachandt Janine Franklin Fergus Willmore Tim Atkinson 《Hydrogeology Journal》2012,20(3):483-498
The Corallian limestone of northern England (UK) is widely exploited for water supplies and exhibits the karstic phenomena of sinking rivers, conduit development and groundwater velocities of several kilometres per day. To test a number of model-derived source protection zones and elucidate contaminant transport mechanisms in the aquifer, three tracer tests were conducted from a set of swallow-holes draining the River Derwent toward public water supply wells in the eastern part of the aquifer. Tracers used included: Enterobacter cloacae (bacteriophage), Photine C (optical brightener), sodium fluorescein (fluorescent dye) and sulphur hexafluoride (dissolved gas), the varying properties of which make them suitable analogues for different types of potential contaminant. Observed tracer transport times and arrival patterns indicate that tracer transport occurs through karstic channels embedded in a network of primary fissures which exert control over tracer concentrations once initial tracer plumes have passed. A dipole flow system is observed between the swallow-holes and the closest abstraction well, whilst previously modelled source protection zones do not accurately reflect either groundwater velocity or those areas of the aquifer supplying the wells. These findings imply that managing such aquifers for potential contamination should rely upon empirical tracer evidence for source-protection zone modelling. 相似文献
15.
Experimental determination of three-dimensional dispersivities in homogeneous porous medium 总被引：1，自引：1，他引：0
Laboratory tracer experiments were conducted to investigate longitudinal dispersivity (α
x
) as well as the transversal (α
y
) and vertical (α
z
) dispersivities in homogeneous 3–5 mm sandy aquifer. The experiments were carried out in a channel 12-m long, 1.35-m wide
and 0.60-m high which was built in the Hydraulics Laboratory of Civil Engineering Department in Dokuz Eylul University. NaCl
was used as a tracer and conductivity values were measured at 220 measurement points. Mass Transport 3 Dimensional (Zheng
and Wang in SERDP-99-1, US Army Engineer Research and Development Center, Vicksburg, MS, 1999; MT3DMS code) which is a three-dimensional solute transport simulation model incorporating finite differences solution option
was used to solve the three-dimensional advective–dispersive transport equation. The estimated dispersivity values were modified
until an acceptable compatibility between the observed and calculated concentrations at measurement points was reached. The
best match was obtained for α
x
= 12 cm, α
y
/α
x
= 0.2 and α
z
/α
x
= 0.05. These values are compatible with those encountered in the literature. 相似文献
16.
Effective porosity value was analyzed from the tritium concentration of sampled groundwater using a three-dimensional groundwater-flow and advection-dispersion code based on the finite element method. The effective porosity value was about 10%. Porosity values measured from core samples were 7–15%. The groundwater flow velocity estimated from the tritium concentrations was about 1 × 10–5 cm s–1. Therefore, during the low groundwater flow velocity condition, effective porosity and porosity values were the same. At the same test site, a 0.48% effective porosity value, determined by another tracer test injecting Br– solution into the aquifer during groundwater level change, was smaller than the porosity value when the flow velocity was 1.8×10–2 cm s–1. Thus the effective porosity value is concluded to be due to groundwater flow velocity. The specific yield value was calculated to be 0.6% by the total volume of tunnel seepage water and the total volume of the rock unsaturated during tunnel construction. However, as pore water continued to be drained after the groundwater level change was completed, the specific yield value became larger than 0.6%. Thus specific yield value is concluded to be due to drainage time. 相似文献
17.
径向弥散是指溶质在径向流场下的迁移规律,被广泛用于描述含水层修复领域中污染物的迁移过程.然而,在现有描述径向弥散的模型中,往往忽略了井内混合效应对溶质径向弥散的影响.建立新的注入井附近溶质径向运移动力学模型,同时考虑井内混合效应与弥散度的尺度效应.采用Laplace变换推导该模型的半解析解,利用Stehfest数值逆变换获取溶质在实数空间的解.通过与不考虑混合效应的模型对比研究混合效应对溶质径向弥散的影响,并利用室内渗流槽中的溶质径向弥散实验数据验证模型的合理性与适用性.结果表明:混合效应和尺度效应对注水井附近溶质径向弥散有显著影响.具体地讲,井内的混合效应越显著,在井壁处及含水层中的穿透曲线越低,溶质浓度达到峰值所需时间越长,与不考虑混合效应模型的差异越明显;随尺度效应的增强,溶质提前穿透且扩散范围变大,溶质浓度达到峰值所需时间越长;与前人的模型相比,本研究模型能更好地模拟注水井附近的溶质径向弥散问题. 相似文献
18.
Sergio E. Serrano 《Mathematical Geology》1996,28(6):701-721
A large body of existing theories of flow and contaminant transport in aquifers ignore the presence of recharge, eliminate the boundary conditions, neglect transient conditions in groundwater flow, conceive hydraulic gradients as linear, and require parameter variability to be stationary and Gaussian. The most outstanding and difficult to justify assumption is the subjective small size of the stochastic terms (i.e., small perturbation methods), which usually is forced by considering the logarithm of the hydraulic conductivity. Several problems in flow and contaminant subsurface hydrology, such as the enhanced dispersion parameters with plume size or time after injection, remain to be observed in the light of a stochastic theory that allows a more realistic consideration of physical and hydrologic properties. In this article, an attempt is made to reformulate a contaminant transport equation (the variable dispersion equation, VDE) with transport parameters in terms of regional hydrologic and aquifer hydraulic properties, such as recharge rate, spatially random transmissivity, hydraulic gradient, aquifer thickness, and soil porosity. Subsequently, a general analytic procedure, the method of decomposition, is used to derive a solution to the VDE. This procedure does not require small perturbation, logarithmic transformations, or specific probability law assumptions. Comparison tests with existing theoretical and field results are given. The tests illustrate the enhanced dispersion and shifting concentration effects produced by the variable dispersion equation. Finally a generalization of the method to nonstationary dispersion in three-dimensional domains is proposed. 相似文献
19.
山东烟台夹河中、下游地区海水入侵三维水质数值模拟研究 总被引：17，自引：0，他引：17
建立了三维变密度对流弥散水质数学模型来研究山东省烟台夹河中、下游地区咸淡水界面的运移规律。以四面体为基本离散单元 ,推导出三维海水入侵变密度水质模型求解的数值方法 ,其中水流方程求解时运用了迦辽金有限单元法。溶质运移方程求解时运用了欧拉拉格朗日混合方法 ,将对流项与弥散项分离 ,用传统迦辽金有限元方法求解弥散项 ;采用自适应MOC MMOC法求解对流项 ,以消除人工过量和数值弥散。根据地下水的潮汐效应观测信息 ,确定了含水系统的海底延伸边界 ;利用该地区地下水水头及水质长观资料识别了模型的水文地质参数 ,探讨了夹河地区海水入侵的原因 :认为夹河下游地区滨海地带地下水过量开采是造成烟台地区海水入侵的主要原因。此外 ,海水随潮定期地倒灌进入夹河 ,通过局部岩性天窗侵入淡水含水层加剧了沿夹河河床两侧地下水的咸化。同时还预测了几种情况下地下水的水质演化趋势 ,为防止和减轻夹河地区海水入侵提供合理、科学的依据。 相似文献
20.
E. César S. Wildemeersch Ph. Orban S. Carrière S. Brouyère A. Dassargues 《Hydrogeology Journal》2014,22(5):1087-1100
Models are the only tools capable of predicting the evolution of groundwater systems at a regional scale, by taking into account a large amount of information. This study presents the association of a water balance model (WetSpass) with a groundwater flow and solute transport model (SUFT3D, saturated and unsaturated flow and transport in 3D) in order to simulate the present and future groundwater quality in terms of nitrate in the Upper Dyle basin (439 km2) Belgium. The HFEMC (hybrid finite element mixing cell) method implemented in the SUFT3D code is used to model groundwater flow and nitrate transport. Spatially distributed recharge, modelled with WetSpass, is considered for prescribing the recharge to the groundwater flow model. The feasibility of linking the WetSpass model with the finite-elements SUFT3D code is demonstrated. Time evolution and distribution of nitrate concentration are then simulated using the calibrated model. Nitrate inputs are spatially distributed according to land use. The spatial simulations and temporal trends are compared with previously published data on this aquifer and show good results. 相似文献