首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
非均质含水层中渗流与溶质运移研究进展   总被引:1,自引:0,他引:1  
Natural aquifer heterogeneity controls the groundwater flow and solute transport, and how to accurately quantify the flow and solute transport in heterogeneous aquifers has received wide attention by many scholars, and has become a hot research topic in earth science. Theoretically, a systematic review is given by the following aspect: flow and solute transport model, moment analysis, multi scale analysis. The resolved and remained issues for scale conversion in hydrogeology research are pointed out. Secondly, recent advances of heterogeneous field test, uncertainty and velocity connectivity are analyzed. Finally, the geophysical inversion of aquifer heterogeneity, stochastic theory and development of stochastic simulation software, scale conversion and uncertainty of velocity connectivity, and the relationship between heterogeneity and hydrogeological condition on the major four aspects of the future research direction is pointed out.  相似文献   

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

3.
多尺度非均质多孔介质中溶质运移的蒙特卡罗模拟   总被引:4,自引:0,他引:4       下载免费PDF全文
探讨了将蒙特卡罗(Monte Carlo)方法应用于多尺度非均质含水层中溶质运移模拟的方法。所研究的含水层由两种具有不同渗透系数统计特征的多孔介质所组成,每一种多孔介质是非均质的,且其渗透系数场符合平稳假设,而整个模拟区的渗透系数是非平稳的。Monte Carlo方法要求参数是平稳的,因此,分别对两种多孔介质产生若干随机渗透系数场后,用两种方法进行组合,并进行溶质运移的模拟计算。通过对计算结果的分析,综合考虑计算精度、计算时间等因素,得出了处理多尺度非均质多孔介质中溶质运移问题的较好方法。  相似文献   

4.
Saturated flow takes place in geological formations of spatially variable permeability which is regarded as a stationary random space function of given statistical moments. The flow is assumed to be uniform in the mean and the Eulerian velocity field has stationary fluctuations. Water carries solutes that react according to the nonlinear equilibrium Freundlich isotherm. Neglecting pore scale dispersion (high Peclet number), we study the behavior of an initially finite pulse injection of constant concentration. Mean flux-averaged concentration is derived in a simple manner by using the previously determined solution of transport in a homogeneous one-dimensional medium and the Lagrangian methodology developed by Cvetkovic and Dagan [5] to model reactive transport in a three-dimensional flow field. The mean breakthrough curves are computed and the combined effect of reactive parameters and heterogeneity upon reduction of the concentration peak is investigated. Furthermore, with the aid of temporal moments, we determine equivalent reaction and macrodispersion coefficients pertinent to a homogeneous medium. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

5.
The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

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.  相似文献   

7.
In porous aquifers, groundwater flow and solute transport strongly depend on the sedimentary facies distribution at fine scale, which determines the heterogeneity of the conductivity field; in particular, connected permeable sediments could form preferential flow paths. Therefore, properly defined statistics, e.g. total and intrinsic facies connectivity, should be correlated with transport features. In order to improve the assessment of the relevance of this relationship, some tests are conducted on two ensembles of equiprobable realizations, obtained with two different geostatistical simulation methods—sequential indicator simulation and multiple point simulation (MPS)—from the same dataset, which refers to an aquifer analogue of sediments deposited in a fluvial point-bar/channel association. The ensembles show different features; simulations with MPS are more structured and characterised by preferential flow paths. This is confirmed by the analysis of transport connectivities and by the interpretation of data from numerical experiments of conservative solute transport with single and dual domain models. The use of two ensembles permits (1) previous results obtained for single realizations to be consolidated on a more firm statistical basis and (2) the application of principal component analysis to assess which quantities are statistically the most relevant for the relationship between connectivity indicators and flow and transport properties.  相似文献   

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

9.
Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport.  相似文献   

10.
Dispersive mass transport processes in naturally heterogeneous geological formations (porous media) are investigated based on a particle approach to mass transport and on its numerical implementation using LPT3D, a Lagrangian Particle Tracking 3D code. We are currently using this approach for studying microscale and macroscale space–time behavior (advection, diffusion, dispersion) of tracer plumes, solutes, or miscible fluids, in 1,2,3-dimensional heterogeneous and anisotropic subsurface formations (aquifers, petroleum reservoirs). Our analyses are based on a general advection-diffusion model and numerical scheme where concentrations and fluxes are discretized in terms of particles. The advection-diffusion theory is presented in a probabilistic framework, and in particular, a numerical analysis is developed for the case of advective transport and rotational flows (numerical stability of the explicit Euler scheme). The remainder of the paper is devoted to the behavior of concentration, mass flux density, and statistical moments of the transported tracer plume in the case of heterogeneous steady flow fields, where macroscale dispersion occurs due to geologic heterogeneity and stratification. We focus on the case of perfectly stratified or multilayered media, obtained by generating many horizontal layers with a purely random transverse distribution of permeability and horizontal velocity. In this case, we calculate explicitly the exact mass concentration field C(x, t), mass flux density field f(x, t), and moments. This includes spatial moments and dispersion variance 2 x (t) on a finite domain L, and temporal moments on a finite time scale T, e.g., the mass variance of arrival times 2 T (x). The moments are related to flux concentrations in a way that takes explicitly into account finite space–time scales of analysis (time-dependent tracer mass; spatially variable flow through mass). The multilayered model problem is then used in numerical experiments for testing different ways of recovering information on tracer plume migration, dispersion, concentration and flux fields. Our analyses rely on a probabilistic interpretation that emerges naturally from the particle approach; it is based on spatial moments (particle positions), temporal moments (mass weighted arrival times), and probability densities (both concentrations and fluxes). Finally, as an alternative to direct estimations of the flux and concentration fields, we formulate and study the Moment Inverse Problem. Solving the MIP yields an indirect method for estimating the space–time distribution of flux concentrations based on observed or estimated moments of the plume. The moments may be estimated from field measurements, or numerically computed by particle tracking as we do here.  相似文献   

11.
肖勋  施文光  王全荣 《地球科学》2020,45(4):1439-1446
径向弥散是指溶质在径向流场下的迁移规律,被广泛用于描述含水层修复领域中污染物的迁移过程.然而,在现有描述径向弥散的模型中,往往忽略了井内混合效应对溶质径向弥散的影响.建立新的注入井附近溶质径向运移动力学模型,同时考虑井内混合效应与弥散度的尺度效应.采用Laplace变换推导该模型的半解析解,利用Stehfest数值逆变换获取溶质在实数空间的解.通过与不考虑混合效应的模型对比研究混合效应对溶质径向弥散的影响,并利用室内渗流槽中的溶质径向弥散实验数据验证模型的合理性与适用性.结果表明:混合效应和尺度效应对注水井附近溶质径向弥散有显著影响.具体地讲,井内的混合效应越显著,在井壁处及含水层中的穿透曲线越低,溶质浓度达到峰值所需时间越长,与不考虑混合效应模型的差异越明显;随尺度效应的增强,溶质提前穿透且扩散范围变大,溶质浓度达到峰值所需时间越长;与前人的模型相比,本研究模型能更好地模拟注水井附近的溶质径向弥散问题.   相似文献   

12.
The combined influence of temporal fluctuations and spatial heterogeneity on non-reactive solute transport mechanisms in porous media can be understood by performing simulations of steady and unsteady flow and transport in heterogeneous media. The study focuses on issues such as the degree of heterogeneity, correlation length, separation of the combined effects of temporal and spatial variations, and ergodicity conditions under unsteady flow conditions. It is shown that the effect of temporal variations on solute transport is masked by the strong effect of spatial heterogeneity. There is no obvious difference in plume shape between steady and unsteady flow conditions; the first and the second spatial moments of the plume of the unsteady-state flow condition fluctuate around the steady-state flow condition with the same period of oscillations as the input signal at small storage coefficient (S????0.001). At a relatively high standard deviation in hydraulic conductivity and a small storage coefficient, the unsteady flow condition sharpens the temporal variations in macrodispersion coefficients. The magnitude of the longitudinal macrodispersion coefficient under unsteady flow condition is almost doubled at the maximum values. However, the transverse macrodispersion coefficient fluctuates around zero. The Kubo number and Peclet number ranges are 1.2?C64 and 10?C250, respectively.  相似文献   

13.
徐绍辉  张佳宝  刘建立 《地球科学》2001,26(5):529-532,550
描述非饱和土壤中溶质运移的对流弥散方程可分成两部分:对流部分用通量校正运移(FCT)算法求解;弥散部分用常规的隐式差分方法求解.FCT算法包括两个阶段,一个是低阶运移阶段,这一阶段的解,可能会引进过量的数值弥散;另一个是高阶通量校正阶段,通过对反扩散通量进行校正(限定),可有效地消除数值弥散和数值振荡.而水体积分数用FUCG方法求得,能保持质量守恒.通过数值例子验证了FCT算法的有效性.  相似文献   

14.
岩溶区地下水数值模拟研究进展   总被引:2,自引:2,他引:0  
岩溶含水介质的不均一性导致岩溶地下水流动、溶质运移和热量迁移的数学模拟研究成为地下水模拟的难点。本文综述了岩溶区地下水流模拟的几种方法,重点阐述了等效多孔介质法、双重连续介质法和三重介质法的定义、发展过程和适用范围,并回顾了这几种方法的研究成果。从等效多孔介质法到三重介质法,模拟精度不断提高,适用范围也逐渐由大区域实际问题向小区域理论研究过渡。介绍了溶质运移模拟和热迁移模拟的研究方法及实例。溶质运移模拟以对流弥散方程为基础,其中尺度效应是溶质运移模拟的重点研究问题;热量迁移模拟应考虑地下热水密度变化对地下热水运动的影响。溶质运移模拟和热量迁移模拟往往是将迁移模型和已经调试成功的地下水流动模型相耦合,从而达到模拟溶质及热量迁移的目的。由于溶质运移和热量迁移的复杂性,现阶段水流模型多数处于等效多孔介质模型阶段。综合理论及实际应用,指出精确刻画裂隙及管道和注重基础数学算法是岩溶水数值模拟进步的关键。   相似文献   

15.
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.  相似文献   

16.
基于典型畦灌施肥试验观测结果及其模拟结果,对比分析利用混合数值解法和Roe有限体积法分别求解一维畦灌施肥地表水流与溶质运移过程控制方程在数值稳定性与收敛性、计算精度与效率上的差异,验证混合数值解法的计算性能与模拟效果.结果表明,混合数值解法比Roe有限体积法表现出更佳的数值稳定性和收敛性,产生的水平衡误差和平均相对误差...  相似文献   

17.
A Eulerian stochastic method is applied to develop a theory of concentration variance for solute transport in a heterogeneous medium. The study focuses on the effects of kinetic sorption and local dispersion on solute dissipation. Spatial distribution of the concentration variance is obtained by scaling the zero local dispersion form of σc 2. The scaling function resulting from the local dispersion and kinetic sorption is derived in a closed integral form. It satisfies the measurement of total concentration variance resulting from the Eulerian mass balance using spatially integrated concentration moments. The spatially integrated moments bypass the need for classical closures applied to joint moments between concentration and velocity fields. The study results indicate that kinetic sorption reduces the total development of concentration variance in comparison with non-reactive solute transport. Kinetic sorption acts as a reduction mechanism, but not as a dissipating mechanism like the local dispersion. Kinetic sorption and local dispersion are not additive processes and their effects on the concentration variance depend on the stage of transport time.  相似文献   

18.
为深入探究含水层中多孔介质的不规则形态和分布对于地下水流场及溶质运移的作用,根据沉积环境中已被发现的倾斜交叉分层地质单元结构,人为构造多孔介质双层交叉鱼骨结构,使其形成宏观各向异性水力传导特征,通过数值模拟研究该结构及其空间位置对于流场和溶质运移的影响。研究发现:宏观各向异性多孔介质引发了螺旋状流动,导致流线的拉伸和折叠,使得溶质羽发生不规则变形,显著提升了非反应性溶质的稀释程度;宏观各向异性结构和溶质羽在三维空间中的垂向相对位置对于溶质运移具有显著影响。  相似文献   

19.
A three-dimensional model for non-reactive solute transport in physically homogeneous subsurface porous media is presented. The model involves solution of the advection-dispersion equation, which additionally considered temporally dependent dispersion. The model also account for a uniform flow field, first-order decay which is inversely proportional to the dispersion coefficient and retardation factor. Porous media with semi-infinite domain is considered. Initially, the space domain is not solute free. Analytical solutions are obtained for uniform and varying pulse-type input source conditions. The governing solute transport equation is solved analytically by employing Laplace transformation technique (LTT). The solutions are illustrated and the behavior of solute transport may be observed for different values of retardation factor, for which simpler models that account for solute adsorption through a retardation factor may yield a misleading assessment of solute transport in ‘‘hydrologically sensitive’’ subsurface environments.  相似文献   

20.
梁越  王俊杰  刘明维 《岩土力学》2015,36(10):3007-3014
垃圾填埋场滤出液、入侵海水、核废物及生产生活废水等污染物随地下水的迁移威胁人类生存。地下水渗流的随机性导致溶质运移问题更加复杂。根据流网特点,利用流线与水头等势线对求解域进行离散,基于质量守恒原理建立流网单元内溶质浓度求解的隐式有限体积差分格式。基于多孔介质孔隙流速分布规律,利用蒙特卡洛法建立流管单元随机流速场进行溶质运移过程的数值模拟,最后根据数值模拟和模型试验的结果对比,验证了数值模拟方法的准确性。基于流网单元的数值模型中沿流线方向的物质交换由对流和扩散共同作用,而流管间的物质交换只有扩散作用,因此,可在不使用弥散系数下进行污染物运移的模拟。引入随机方法确定流管内流速为研究非均匀流场中染污物的优势迁移提供了新的思路。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号