共查询到20条相似文献，搜索用时 31 毫秒
1.
A Lagrangian perturbation method is applied to develop a method of moments for reactive solute flux through a three-dimensional, nonstationary flow field. The flow nonstationarity may stem from medium nonstationarity, finite domain boundaries, and/or fluid pumping and injecting. The reactive solute flux is described as a space–time process where time refers to the solute flux breakthrough in a control plane at some distance downstream of the solute source and space refers to the transverse displacement distribution at the control plane. The analytically derived moments equations for solute transport in a nonstationary flow field are too complicated to solve analytically; therefore, a numerical finite difference method is implemented to obtain the solutions. This approach combines the stochastic model with the flexibility of the numerical method to boundary and initial conditions. The approach provides a tool to apply stochastic theory to reactive solute transport in complex subsurface environments. Several case studies have been conducted to investigate the influence of the physical and chemical heterogeneity of a medium on the reactive solute flux prediction in nonstationary flow field. It is found that both physical and chemical heterogeneity significantly affect solute transport behavior in a nonstationary flow field. The developed method is also applied to an environmental project for predicting solute flux in the saturated zone below the Yucca Mountain Project area, demonstrating the applicability of the method in practical environmental projects. 相似文献
2.
Numerical evaluation of multicomponent cation exchange reactive transport in physically and geochemically heterogeneous porous media 总被引：1，自引：0，他引：1
Experimental evidence and stochastic studies strongly show that the transport of reactive solutes in porous media is significantly
influenced by heterogeneities in hydraulic conductivity, porosity, and sorption parameters. In this paper, we present Monte
Carlo numerical simulations of multicomponent reactive transport involving competitive cation exchange reactions in a two-dimensional
vertical physically and geochemically heterogeneous medium. Log hydraulic conductivity, log K, and log cation exchange capacity (log CEC) are assumed to be random Gaussian functions with spherical semivariograms. Random
realizations of log K and log CEC are used as input data for the numerical simulation of multicomponent reactive transport with CORE2D, a general purpose reactive transport code. Longitudinal features of the fronts of reactive and conservative species are
computed from the temporal and spatial moments of depth-averaged concentrations. Monte Carlo simulations show that: (1) the
displacement of reactive fronts increases with increasing variance of log K, while it decreases with the variance of log CEC; (2) second-order spatial moments increase with increasing variances of
log K and log CEC; (3) uncertainties in the mean arrival time are largest (smallest) for negatively (positively) correlated log
K and Log CEC; (4) cations undergoing competitive cation exchange exhibit different apparent velocities and retardation factors
due to both physical and geochemical heterogeneities; and (5) the correlation between log K and log CEC affects significantly apparent cation retardation factors in heterogeneous aquifers. 相似文献
3.
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. 相似文献
4.
This paper considers advective transport of a soluble contaminant through saturated soil with non-linear sorption of the contaminant onto a stationary porous media. The non-linear sorption isotherms considered in the transport analysis are the Langmuir and Freundlich sorption isotherms. A special case of the Freundlich sorption isotherm is the linear sorption isotherm, and it is shown that in this case transport through a homogeneous soil results in the initial concentration profile simply being translated in the direction of the groundwater flow. However, when the sorption isotherm is non-linear the initial concentration profile distorts as it is translated with the groundwater flow, leading to the development of concentration shock fronts and rarefactions. Analytic solutions to the non-linear first-order hyperbolic equations are developed for a number of contaminant transport problems of practical significance. It is shown that in the case of the Langmuir sorption isotherms, shock fronts develop at the leading edge of the concentration profile while for the Freundlich sorption isotherm shock fronts may develop at either the leading or trailing edge of the concentration profile. Copyright © 1999 John Wiley & Sons Ltd. 相似文献
5.
This paper presents an approach conducive to an evaluation of the probability density function (pdf) of spatio-temporal distributions
of concentrations of reactive solutes (and associated reaction rates) evolving in a randomly heterogeneous aquifer. Most existing
approaches to solute transport in heterogeneous media focus on providing expressions for space–time moments of concentrations.
In general, only low order moments (unconditional or conditional mean and covariance) are computed. In some cases, this allows
for obtaining a confidence interval associated with predictions of local concentrations. Common applications, such as risk
assessment and vulnerability practices, require the assessment of extreme (low or high) concentration values. We start from
the well-known approach of deconstructing the reactive transport problem into the analysis of a conservative transport process
followed by speciation to (a) provide a partial differential equation (PDE) for the (conditional) pdf of conservative aqueous
species, and (b) derive expressions for the pdf of reactive species and the associated reaction rate. When transport at the
local scale is described by an Advection Dispersion Equation (ADE), the equation satisfied by the pdf of conservative species
is non-local in space and time. It is similar to an ADE and includes an additional source term. The latter involves the contribution
of dilution effects that counteract dispersive fluxes. In general, the PDE we provide must be solved numerically, in a Monte
Carlo framework. In some cases, an approximation can be obtained through suitable localization of the governing equation.
We illustrate the methodology to depict key features of transport in randomly stratified media in the absence of transverse
dispersion effects. In this case, all the pdfs can be explicitly obtained, and their evolution with space and time is discussed. 相似文献
6.
7.
The Simpevarp area is one of the alternative sites being considered for the deep geological disposal of high level radioactive waste in Sweden. In this paper, a coupled regional groundwater flow and reactive solute transport model of the Simpevarp area is presented that integrates current hydrogeological and hydrochemical data of the area. The model simulates the current hydrochemical pattern of the groundwater system in the area. To that aim, a conceptual hydrochemical model was developed in order to represent the dominant chemical processes. Groundwater flow conditions were reproduced by taking into account fluid-density-dependent groundwater flow and regional hydrogeologic boundary conditions. Reactive solute transport calculations were performed on the basis of the velocity field so obtained. The model was calibrated and sensitivity analyses were carried out in order to investigate the effects of heterogeneities of hydraulic conductivity in the subsurface medium. Results provided by the reactive transport model are in good agreement with much of the measured hydrochemical data. This paper emphasizes the appropriateness of the use of reactive solute transport models when water-rock interaction reactions are involved, and demonstrates what powerful tools they are for the interpretation of hydrogeological and hydrochemical data from site geological repository characterization programs, by providing a qualitative framework for data analysis and testing of conceptual assumptions in a process-oriented approach. 相似文献
8.
随机水文地质学方法,较传统的确定性方法而言,是解决非均质含水层中水流和溶质运移问题的一种更为合理的手段。据以往研究,假设渗透系数场遵循对数正态分布,利用直接傅立叶变换方法来生成渗透系数随机场。应用基于随机理论的蒙特卡罗方法,来研究渗透系数的空间变异性对污染物运移结果的影响。实例研究表明,污染物在含水层中运移过程中污染羽的展布范围(二阶矩)随着渗透系数空间变异方差的增大而扩大,而污染羽在空间上的质心位置(一阶矩)基本不受方差的影响,仅取决于渗透系数随机场的均值大小。另外还分析了污染羽在各点的浓度变化方差和变异系数分别随渗透系数变异方差的变化状况。 相似文献
9.
We have calculated the coefficient of turbulent diffusion in a random flow with time restoration, describing the interstellar medium. Such a flow abruptly loses its memory at random times, forming a Poisson flow of events. The coefficient of turbulent diffusion in the flow is determined by the rms velocity and correlation time, as in mixing-length theory, but the numerical coefficient differs from that predicted by this theory. The closure equation derived by us for the transport of the mean concentration of a passive scalar takes a more complicated form than obtained in standard mean-field theory, but the main properties of the equation retain their validity. The possibility of extending the results of this exactly solved problem to the problem of transport in the turbulent interstellar medium is discussed. 相似文献
10.
以模拟农田尺度降雨入渗—重分布条件下阿特拉津(atrazine)在非饱和土壤中的淋溶动态为数值算例,通过在北京东郊的一块27m×27m的农田采集100个土壤样品,测定其主要理化性质,间接计算获得各土壤样品的van Genuchten型水力学参数、溶质运移的弥散度及atrazine在土壤中吸附的分配系数,进而基于柱模型假设,运用HY DRUS 1D软件得到不同时刻采样区内土壤水负压和atrazine浓度的平均分布。在此基础上,重点探讨了应用幂平均算法计算不同的等价有效参数对农田尺度下土壤水负压和atrazine浓度平均分布的数值模拟精度的影响。就所设计的降雨入渗—重分布条件下atrazine淋溶的算例而言,土壤水力学参数、溶质运移参数和吸附特性参数取几何平均或调和平均作为有效参数值,空间平面上呈正态分布的土壤水负压的初值采用算术平均值作为流场模拟的初始条件,可以获得满意的数值模拟结果。 相似文献
11.
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. 相似文献
12.
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. 相似文献
13.
Due to changes in lithostatic pressure, differential fracturing across bedding planes and irregularities in depositional environments, hydraulic conductivity exhibits heterogeneities and trends at various spatial scales. Using spectral theory, we have examined the effect of trends in hydraulic conductivity on (1) the solution of the mean equation for hydraulic head, (2) the covariance of hydraulic head, (3) the cross-covariances of hydraulic head and log-hydraulic conductivity perturbations and their gradients, and (4) the effective hydraulic conductivity. It is shown that the field of hydraulic head is sensitive to the presence of trends in ways that cannot be predicted by the classical analysis based on stationary hydraulic conductivity fields. The controlling variables for the second moments of hydraulic head are the mean hydraulic gradient, the correlation scale of log-hydraulic conductivity and its variance, and the slope of the trend in log-hydraulic conductivity. The mean hydraulic gradient introduces complications in the analysis since it is, in general, spatially variable. In this respect, our results are approximate, yet indicative of the true role of spatially variable patterns of log-hydraulic conductivity on groundwater flow systems. 相似文献
14.
Integration of Tracer Test Data to Refine Geostatistical Hydraulic Conductivity Fields Using Sequential Self-Calibration Method 总被引：3，自引：0，他引：3
On the basis of local measurements of hydraulic conductivity,geostatistical methods have been found to be useful in heterogeneity characterization of a hydraulic conductivity field on a regional scale. However,the methods are not suited to directly integrate dynamic production data,such as,hydraulic head and solute concentration,into the study of conductivity distribution. These data,which record the flow and transport processes in the medium,are closely related to the spatial distribution of hydraulic conductivity. In this study,a three-dimensional gradient-based inverse method-the sequential self-calibration (SSC) method-is developed to calibrate a hydraulic conductivity field,initially generated by a geostatistical simulation method,conditioned on tracer test results. The SSC method can honor both local hydraulic conductivity measurements and tracer test data. The mismatch between the simulated hydraulic conductivity field and the reference true one,measured by its mean square error (MSE),is reduced through the SSC conditional study. In comparison with the unconditional results,the SSC conditional study creates the mean breakthrough curve much closer to the reference true curve,and significantly reduces the prediction uncertainty of the solute transport in the observed locations. Further,the reduction of uncertainty is spatially dependent,which indicates that good locations,geological structure,and boundary conditions will affect the efficiency of the SSC study results. 相似文献
15.
探讨了将蒙特卡罗(Monte Carlo)方法应用于多尺度非均质含水层中溶质运移模拟的方法。所研究的含水层由两种具有不同渗透系数统计特征的多孔介质所组成,每一种多孔介质是非均质的,且其渗透系数场符合平稳假设,而整个模拟区的渗透系数是非平稳的。Monte Carlo方法要求参数是平稳的,因此,分别对两种多孔介质产生若干随机渗透系数场后,用两种方法进行组合,并进行溶质运移的模拟计算。通过对计算结果的分析,综合考虑计算精度、计算时间等因素,得出了处理多尺度非均质多孔介质中溶质运移问题的较好方法。 相似文献
16.
Conservative tracer experiments can provide information useful for characterizing various subsurface transport properties. This study examines the effectiveness of three different types of transport observations for sensitivity analysis and parameter estimation of a three-dimensional site-specific groundwater flow and transport model: conservative tracer breakthrough curves (BTCs), first temporal moments of BTCs (m 1), and tracer cumulative mass discharge (M d) through control planes combined with hydraulic head observations (h). High-resolution data obtained from a 410-day controlled field experiment at Vandenberg Air Force Base, California (USA), have been used. In this experiment, bromide was injected to create two adjacent plumes monitored at six different transects (perpendicular to groundwater flow) with a total of 162 monitoring wells. A total of 133 different observations of transient hydraulic head, 1,158 of BTC concentration, 23 of first moment, and 36 of mass discharge were used for sensitivity analysis and parameter estimation of nine flow and transport parameters. The importance of each group of transport observations in estimating these parameters was evaluated using sensitivity analysis, and five out of nine parameters were calibrated against these data. Results showed the advantages of using temporal moment of conservative tracer BTCs and mass discharge as observations for inverse modeling. 相似文献
17.
1,2,4-trichlorobenzene (1,2,4-TCB) is one of the most common pollutants in landfill leachate that impacts shallow drinking groundwater quality; accordingly, the transport mechanism and capacity for remediation of this compound are important to the management of local water resource. In this study, a series of indoor experiments and analytical predictions revealed that the transportation of 1,2,4-TCB in groundwater is similar in medium sand and fine sand. Specifically, the peak time increases with distance between monitoring points and the source point, but the peak relative concentrations decrease with distance, indicating that the concentration of 1,2,4-TCB in groundwater is controlled by adsorption and biodegradation. In addition, transportation of water in a fine sand column was lower than that in a medium sand column; therefore, the adsorption and biodegradation of 1,2,4-TCB in groundwater was lower in the medium sand column. These findings demonstrate that 1,2,4-TCB in groundwater could be removed by natural degradation after about 5 years, but that other methods (permeable reactive barriers, air sparging or biosparging) should be applied in the field to shorten the remediation period and enhance water supply safety around landfill. 相似文献
18.
Stochastic analysis of the effect of spatial variability of diffusion parameters on radionuclide transport in a low permeability clay layer 总被引：1，自引：0，他引：1
Most studies that incorporate subsurface heterogeneity in groundwater flow and transport models only analyze and simulate the spatial variability of hydraulic conductivity. Heterogeneity of the other flow and transport parameters are usually neglected. This approach is often justified, but there are, however, cases in which disregarding the heterogeneity of the other flow and transport parameters can be questionable. In low permeability media, for instance, diffusion is often the dominant transport mechanism. It therefore seems logical to incorporate the spatial variability of the diffusion parameters in the transport model. This study therefore analyses and simulates the spatial variability of the effective diffusion coefficient and the diffusion accessible porosity with geostatistical techniques and incorporates their heterogeneity in the transport model of a low permeability formation. The formation studied was Boom clay (Belgium), a candidate host rock for the deep geological disposal of high-level radioactive waste. The calculated output radionuclide fluxes of this model are compared with the fluxes calculated with a homogeneous model and a model with a heterogeneous hydraulic conductivity distribution. This analysis shows that the heterogeneity of the diffusion parameters has a much larger effect on the calculated output radionuclide fluxes than the heterogeneity of hydraulic conductivity in the low permeability medium under study. 相似文献
19.
This paper presents a theoretical approach to analyse coupled, linear thermoporoelastic fields in a saturated porous medium under radial and spherical symmetry. The governing equations account for compressibility and thermal expansion of constituents, heat sink due to thermal dilatation of water and thermal expansion of the medium, and thermodynamically coupled heat–water flow. It has been reported in the literature that thermodynamically coupled heat–water flows known as thermo-osmosis and thermal filtration have the potential to significantly alter the flow fields in clay-rich barriers in the near field of a underground waste containment scheme. This study presents a mathematical model and examines the effects of thermo-osmosis and thermal-filtration on coupled consolidation fields in a porous medium with a cavity. Analytical solutions of the governing equations are presented in the Laplace transform space. A numerical inversion scheme is used to obtain the time-domain solutions for a cylindrical cavity in a homogeneous or a non-homogeneous medium. A closed form time-domain solution is presented for a spherical cavity in a homogeneous medium. Selected numerical solutions for homogeneous and non-homogeneous media show a significant increase in pore pressure and displacements due to the presence of thermodynamically coupled flows and a negligible influence on temperature. © 1998 John Wiley & Sons, Ltd. 相似文献
20.
A. Berez G. Schäfer F. Ayari M. Trabelsi-Ayadi 《International Journal of Environmental Science and Technology》2016,13(7):1625-1640
This work reports new findings on the physisorption of Foron Blue 291 on natural bentonite under both static and dynamic flow conditions of the aqueous solution that may prove important when using the adsorption technique for wastewater treatment. The quantity of dye retained by the solid matrix under both ideal exchange conditions between dyes and adsorbents given in a batch reactor and a pulse injection of dissolved azo dye in a soil column composed of a mixture of mean grain-sized sand and bentonite at different flow rates and at various mass fractions of bentonite was studied. Column experiments involving the use of a non-reactive tracer (Fluorescein) were performed to examine hydrodynamic behaviours of the clay/sand mixture studied. It is shown that advective–dispersive transport across a clay/sand mixture may be characterized by a double porosity medium. The results obtained on the reactive transport of Foron Blue 291 underscore that Foron Blue 291 adsorption depends heavily on the mass fraction of clay in a clay/sand mixture. With a clay mass fraction of 30 %, 77 % of the dissolved Foron Blue 291 mass was irreversibly removed by the adsorbent. At low flow rates, removal ratios obtained from the dynamic reactive system were similar to those obtained through previous static reactor experiments. 相似文献