共查询到20条相似文献，搜索用时 31 毫秒
1.
《Advances in water resources》1998,21(6):487-498
Transport of a sorbing solute in a two-dimensional steady and uniform flow field is modeled using a particle tracking random walk method. The solute is initially introduced from an instantaneous point source. Cases of linear and nonlinear sorption isotherms are considered. Local pore velocity and mechanical dispersion are used to describe the solute transport mechanisms at the local scale. The numerical simulation of solute particle transport yields the large scale behavior of the solute plume. Behavior of the plume is quantified in terms of the center-of-mass displacement distance, relative velocity of the center-of-mass, mass breakthrough curves, spread variance, and longitudinal skewness. The nonlinear sorption isotherm affects the plume behavior in the following way relative to the linear isotherm: (1) the plume velocity decreases exponentially with time; (2) the longitudinal variance increases nonlinearly with time; (3) the solute front is steepened and tailing is enhanced 相似文献
2.
A comprehensive numerical study was undertaken to investigate transport of a variable-density, conservative solute plume in an unconfined coastal aquifer subject to high and low frequency oceanic forcing. The model combined variable-density saturated flow for groundwater and solute transport, and wave hydrodynamics from a 2D Navier–Stokes solver. A sinusoidal tidal signal was specified by implementing time-varying heads at the seaward boundary. The solute plume behavior was investigated under different oceanic forcing conditions: no forcing, waves, tide, and combined waves and tide. For each forcing condition, four different injected solute densities (freshwater, brackish water, seawater, brine) were used to investigate the effects of density on the transport of the injected plume beneath and across the beach face. The plume’s low-order spatial moments were computed, viz., mass, centroid, variance and aspect ratio. The results confirmed that both tide- and wave-forcing produce an upper saline plume beneath the beach face in addition to the classical saltwater wedge. For the no-forcing and tide-only cases (during rising tides), an additional small circulation cell below the beach face was observed. Oceanic forcing affects strongly the solute plume’s flow path, residence time and discharge rate across the beach face, as well as its spreading. For the same oceanic forcing, solute plumes with different densities follow different trajectories from the source to the discharge location (beach face). The residence time and plume spreading increased with plume density. It was concluded that simulations that neglect the effect of waves or tides cannot reproduce accurately solute plume dispersion and also, in the case of coasts with small waves or tides, the solute residence time in the aquifer. 相似文献
3.
Pore-scale dispersion (PSD), aquifer heterogeneity, sampling volume, and source size influence solute concentrations of conservative tracers transported in heterogeneous porous formations. In this work, we developed a new set of analytical solutions for the concentration ensemble mean, variance, and coefficient of variation (CV), which consider the effects of all these factors. We developed these models as generalizations of the first-order solutions in the log-conductivity variance of point concentration proposed by [Fiori A, Dagan G. Concentration fluctuations in aquifer transport: a rigorous first-order solution and applications. J Contam Hydrol 2000;45(1–2):139–163]. Our first-order solutions compare well with numerical simulations for small and moderate formation heterogeneity and from small to large sampling and source volumes. However, their performance deteriorates for highly heterogeneous formations. Successively, we used our models to study the interplay among sampler size, source volume, and PSD. Our analysis shows a complex and important interaction among these factors. Additionally, we show that the relative importance of these factors is also a function of plume age, of aquifer heterogeneity, and of the measurement location with respect to the mean plume center of gravity. We found that the concentration moments are chiefly controlled by the sampling volume with pore-scale dispersion playing a minor role at short times and for small source volumes. However, the effect of the source volume cannot be neglected when it is larger than the sampling volume. A different behavior occurs for long periods, which may be relevant for old contaminations, or for small injection volumes. In these cases, PSD causes a significant dilution, which is reflected in the concentration statistics. Additionally, at the center of the mean plume, where high concentrations are most likely to occur, we found that sampling volume and PSD are attenuating mechanisms for both concentration ensemble mean and coefficient of variation, except at very large source and sampler sizes, where the coefficient of variation increases with sampler size and PSD. Formation heterogeneity causes a faster reduction of the ensemble mean concentrations and a larger uncertainty at the center of the mean plume. Therefore, our results highlight the importance of considering the combined effect of formation heterogeneity, exposure volume, PSD, source size, and measurement location in performing risk assessment. 相似文献
4.
We analyze the impact of conditioning to measurements of hydraulic transmissivity on the transport of a conservative solute. The effects of conditioning on solute transport are widely discussed in the literature, but most of the published works focuses on the reduction of the uncertainty in the prediction of the plume dispersion. In this study both ensemble and effective plume moments are considered for an instantaneous release of a solute through a linear source normal to the mean flow direction, by taking into account different sizes of the source. The analysis, involving a steady and spatially inhomogeneous velocity field, is developed by using the stochastic finite element method. Results show that conditioning reduces the ensemble moment in comparison with the unconditioned case, whereas the effective dispersion may increase because of the contribution of the spatial moments related to the lack of stationarity in the flow field. As the number of conditioning points increases, this effect increases and it is significant in both the longitudinal and transverse directions. Furthermore, we conclude that the moment derived from data collected in the field can be assessed by the conditioned second-order spatial moment only with a dense grid of measured data, and it is manifest for larger initial lengths of the plume. Nevertheless, it seems very likely that the actual dispersion of the plume may be underestimated in practical applications. 相似文献
5.
Transport of sorbing solutes in 2D steady and heterogeneous flow fields is modeled using a particle tracking random walk technique. The solute is injected as an instantaneous pulse over a finite area. Cases of linear and Freundlich sorption isotherms are considered. Local pore velocity and mechanical dispersion are used to describe the solute transport mechanisms at the local scale. This paper addresses the impact of the degree of heterogeneity and correlation lengths of the log-hydraulic conductivity field as well as negative correlation between the log-hydraulic conductivity field and the log-sorption affinity field on the behavior of the plume of a sorbing chemical. Behavior of the plume is quantified in terms of longitudinal spatial moments: center-of-mass displacement, variance, 95% range, and skewness. The range appears to be a better measure of the spread in the plumes with Freundlich sorption because of plume asymmetry. It has been found that the range varied linearly with the travelled distance, regardless of the sorption isotherm. This linear relationship is important for extrapolation of results to predict behavior beyond simulated times and distances. It was observed that the flow domain heterogeneity slightly enhanced the spreading of nonlinearly sorbing solutes in comparison to that which occurred for the homogeneous flow domain, whereas the spreading enhancement in the case of linear sorption was much more pronounced. In the case of Freundlich sorption, this enhancement led to further deceleration of the solute plume movement as a result of increased retardation coefficients produced by smaller concentrations. It was also observed that, except for plumes with linear sorption, correlation between the hydraulic conductivity and the sorption affinity fields had minimal effect on the spatial moments of solute plumes with nonlinear sorption. 相似文献
6.
In karst aquifers with significant matrix permeability, water and solutes are exchanged between the conduits and carbonate matrix. Transport through the matrix increases the spread of solutes and increases travel times. This study numerically evaluates advective solute transport in synthetic karst systems that contain 3D branching conduit networks. Particle tracking is performed to analyze the spatial and temporal transport history of solute that arrives at the conduit outlet. Three measures of transport connectivity are used to quantify the solute migration behavior: the skewness of the particle arrival time distribution, the normalized fifth percentile of arrival times, and the fraction of the total travel time that occurs within conduits. All three of these metrics capture the influence of conduit network geometry on solute transport. A more tortuous network leads to enhanced conduit-matrix mixing, which reduces the transport connectivity and yields a broader distribution of solute arrival times. These results demonstrate that the conduit network geometry is an important control on solute transport in karst systems with a permeable matrix. 相似文献
7.
8.
V. D. Cvetkovic G. Dagan A. M. Shapiro 《Stochastic Environmental Research and Risk Assessment (SERRA)》1991,5(1):45-54
The problem of one-dimensional transport of passive solute by a random steady velocity field is investigated. This problem is representative of solute movement in porous media, for example, in vertical flow through a horizontally stratified formation of variable porosity with a constant flux at the soil surface. Relating moments of particle travel time and displacement, exact expressions for the advection and dispersion coefficients in the Focker-Planck equation are compared with the perturbation results for large distances. The first- and second-order approximations for the dispersion coefficient are robust for a lognormal velocity field. The mean Lagrangian velocity is the harmonic mean of the Eulerian velocity for large distances. This is an artifact of one-dimensional flow where the continuity equation provides for a divergence free fluid flux, rather than a divergence free fluid velocity. 相似文献
9.
The problem of one-dimensional transport of passive solute by a random steady velocity field is investigated. This problem is representative of solute movement in porous media, for example, in vertical flow through a horizontally stratified formation of variable porosity with a constant flux at the soil surface. Relating moments of particle travel time and displacement, exact expressions for the advection and dispersion coefficients in the Focker-Planck equation are compared with the perturbation results for large distances. The first- and second-order approximations for the dispersion coefficient are robust for a lognormal velocity field. The mean Lagrangian velocity is the harmonic mean of the Eulerian velocity for large distances. This is an artifact of one-dimensional flow where the continuity equation provides for a divergence free fluid flux, rather than a divergence free fluid velocity. 相似文献
10.
Previous studies of hyporheic zone focused largely on the net mass transfer of solutes between stream and streambed. Solute transport within the bed has attracted less attention. In this study, we combined flume experiments and numerical simulations to examine solute transport processes in a streambed with periodic bedforms. Solute originating from the stream was subjected to advective transport driven by pore water circulation due to current–bedform interactions as well as hydrodynamic dispersion in the porous bed. The experimental and numerical results showed that advection played a dominant role at the early stage of solute transport, which took place in the hyporheic zone. Downward solute transfer to the deep ambient flow zone was controlled by transverse dispersion at the later stage when the elapsed time exceeded the advective transport characteristic time tc (= L/uc with L being the bedform length and uc the characteristic pore water velocity). The advection-based pumping exchange model was found to predict reasonably well solute transfer between the overlying water and streambed at the early stage but its performance deteriorated at the later stage. With dispersion neglected, the pumping exchange model underestimated the long-term rate and total mass of solute transfer from the overlying water to the bed. Therefore both advective and dispersive transport components are essential for quantification of hyporheic exchange processes. 相似文献
11.
Z. J. Kabala 《Stochastic Environmental Research and Risk Assessment (SERRA)》1997,11(4):331-348
Under the assumption that local solute dispersion is negligible, a new general formula (in the form of a convolution integral)
is found for the arbitrary k-point ensemble moment of the local concentration of a solute convected in arbitrary m spatial
dimensions with general sure initial conditions. From this general formula new closed-form solutions in m=2 spatial dimensions
are derived for 2-point ensemble moments of the local solute concentration for the impulse (Dirac delta) and Gaussian initial
conditions. When integrated over an averaging window, these solutions lead to new closed-form expressions for the first two
ensemble moments of thevolume-averaged solute concentration and to the corresponding concentration coefficients of variation (CV). Also, for the impulse (Dirac delta) solute concentration
initial condition, the second ensemble moment of thesolute point concentration in two spatial dimensions and the corresponding CV are demonstrated to be unbound. For impulse initial conditions the CVs
for volume-averaged concentrations axe compared with each other for a tracer from the Borden aquifer experiment. The point-concentration
CV is unacceptably large in the whole domain, implying that the ensemble mean concentration is inappropriate for predicting
the actual concentration values. The volume-averaged concentration CV decreases significantly with an increasing averaging
volume. Since local dispersion is neglected, the new solutions should be interpreted as upper limits for the yet to be derived
solutions that account for local dispersion; and so should the presented CVs for Borden tracers. The new analytical solutions
may be used to test the accuracy of Monte Carlo simulations or other numerical algorithms that deal with the stochastic solute
transport. They may also be used to determine the size of the averaging volume needed to make a quasi-sure statement about
the solute mass contained in it. 相似文献
12.
The concentration fluctuations resulting from hazardous releases in the subsurface are modeled through the concentration moments. The local solute exposure concentration, resulting from the heterogeneous velocity field and pore scale dispersion in the subsurface, is a random function characterized by its statistical moments. The approximate solution to the exact equation that describes the evolution of concentration standard moments in the aquifer transport is proposed in a recursive form. The expressions for concentration second, third and fourth central moments are derived and evaluated for various flow and transport conditions. The solutions are sought by starting from the exact upper bound solution with the zero pore scale dispersion and introducing the physically based approximation that allows the inclusion of the pore scale dispersion resulting in simple closed-form expressions for the concentration statistical moments. The concentration moments are also analyzed in the relative and absolute frame of reference indicating their combined importance in the practical cases of the subsurface contaminant plume migration. The influence of pore scale dispersion with different source sizes and orientations are analyzed and discussed with respect to common cases in the environmental risk assessment problems. The results are also compared with the concentration measurements of the conservative tracer collected in the field experiments at Cape Cod and Borden Site. 相似文献
13.
A three-dimensional finite-volume ELLAM method has been developed, tested, and successfully implemented as part of the U.S. Geological Survey (USGS) MODFLOW-2000 ground water modeling package. It is included as a solver option for the Ground Water Transport process. The FVELLAM uses space-time finite volumes oriented along the streamlines of the flow field to solve an integral form of the solute-transport equation, thus combining local and global mass conservation with the advantages of Eulerian-Lagrangian characteristic methods. The USGS FVELLAM code simulates solute transport in flowing ground water for a single dissolved solute constituent and represents the processes of advective transport, hydrodynamic dispersion, mixing from fluid sources, retardation, and decay. Implicit time discretization of the dispersive and source/sink terms is combined with a Lagrangian treatment of advection, in which forward tracking moves mass to the new time level, distributing mass among destination cells using approximate indicator functions. This allows the use of large transport time increments (large Courant numbers) with accurate results, even for advection-dominated systems (large Peclet numbers). Four test cases, including comparisons with analytical solutions and benchmarking against other numerical codes, are presented that indicate that the FVELLAM can usually yield excellent results, even if relatively few transport time steps are used, although the quality of the results is problem-dependent. 相似文献
14.
Numerical experiments are conducted to examine the effects of gravity on monodisperse and polydisperse colloid transport in water-saturated fractures with uniform aperture. Dense colloids travel in water-saturated fractures by advection and diffusion while subject to the influence of gravity. Colloids are assumed to neither attach onto the fracture walls nor penetrate the rock matrix based on the assumptions that they are inert and their size is larger than the pore size of the surrounding solid matrix. Both the size distribution of a colloid plume and colloid density are shown to be significant factors impacting their transport when gravitational forces are important. A constant-spatial-step particle-tracking code simulates colloid plumes with increasing densities transporting in water-saturated fractures while accounting for three forces acting on each particle: a deterministic advective force due to the Poiseuille flow field within the fracture, a random force caused by Brownian diffusion, and the gravitational force. Integer angles of fracture orientation with respect to the horizontal ranging from ±90° are considered: three lognormally distributed colloid plumes with mean particle size of 1 μm (averaged on a volumetric basis) and standard deviation of 0.6, 1.2 and 1.8 μm are examined. Colloid plumes are assigned densities of 1.25, 1.5, 1.75 and 2.0 g/cm3. The first four spatial moments and the first two temporal moments are estimated as functions of fracture orientation angle and colloid density. Several snapshots of colloid plumes in fractures of different orientations are presented. In all cases, larger particles tend to spread over wider sections of the fracture in the flow direction, but smaller particles can travel faster or slower than larger particles depending on fracture orientation angle. 相似文献
15.
An empirical hyperbolic scale-dependent dispersion model, which predicts a linear growth of dispersivity close to the origin and the attainment of an asymptotic dispersivity at large distances, is presented for deterministic modelling of field-scale solute transport and the analysis of solute transport experiments. A simple relationship is derived between local dispersivity, which is used in numerical simulations of solute transport, and effective dispersivity, which is estimated from the analysis of tracer breakthrough curves. The scale-dependent dispersion model is used to interpret a field tracer experiment by nonlinear least-squares inversion of a numerical solution for unsaturated transport. Simultaneous inversion of concentration-time data from several sampling locations indicates a linear growth of the dispersion process over the scale of the experiment. These findings are consistent with the results of an earlier analysis based on the use of a constant dispersion coefficient model at each of the sampling depths. 相似文献
16.
A Eulerian analytical method is developed for nonreactive solute transport in heterogeneous, dual-permeability media where the hydraulic conductivities in fracture and matrix domains are both assumed to be stochastic processes. The analytical solution for the mean concentration is given explicitly in Fourier and Laplace transforms. Instead of using the fast fourier transform method to numerically invert the solution to real space (Hu et al., 2002), we apply the general relationship between spatial moments and concentration (Naff, 1990; Hu et al., 1997) to obtain the analytical solutions for the spatial moments up to the second for a pulse input of the solute. Owing to its accuracy and efficiency, the analytical method can be used to check the semi-analytical and Monte Carlo numerical methods before they are applied to more complicated studies. The analytical method can be also used during screening studies to identify the most significant transport parameters for further analysis. In this study, the analytical results have been compared with those obtained from the semi-analytical method (Hu et al., 2002) and the comparison shows that the semi-analytical method is robust. It is clearly shown from the analytical solution that the three factors, local dispersion, conductivity variation in each domain and velocity convection flow difference in the two domains, play different roles on the solute plume spreading in longitudinal and transverse directions. The calculation results also indicate that when the log-conductivity variance in matrix is 10 times less than its counterpart in fractures, it will hardly influence the solute transport, whether the conductivity field is matrix is treated as a homogeneous or random field. 相似文献
17.
Solute transport in heterogeneous media: The impact of anisotropy and non-ergodicity in risk assessment 总被引：1，自引：1，他引：0
X. Sánchez-Vila J. Solís-Delfín 《Stochastic Environmental Research and Risk Assessment (SERRA)》1999,13(5):365-379
Conceptual model selection is a key issue in risk assessment studies. We analyze the effect of a number of conceptual aspects
related to solute transport in two-dimensional heterogeneous media. The main issues addressed are non-ergodicity, anisotropy
in the correlation structure of the transmissivity field, and dispersion at the local scale. In particular, we study the development
of a solute plume when mean flow is oriented at an angle with respect to the principal directions of anisotropy. The study
is carried out in a Lagrangian framework using Monte Carlo analysis.
Of special interest is the evolution of individual plumes. A number of aspects are analyzed, namely the location of the center
of mass for each plume and the different ways to compute the angles that the main axes of the plume develop with respect to
the direction of the mean flow. Stochastic theories based upon ergodicity conclude that the plume gets oriented in the mean
flow direction. In our non-ergodic simulations, the mean of the offset angles, for each individual plume in each particular
realization, is offset from the mean flow direction towards the direction of maximum anisotropy. If, instead, the analysis
is performed on the ensemble plume (superposition of all different simulations), it is then found oriented closer to the direction
of the mean flow than the average offset angle for the different plumes considered separately. This last result adds an extra
word of caution to the use of ensemble averaged values in solute transport studies.
Serious implications for risk assessment follow from the conceptual model adopted. First, in any single realization there
will a large uncertainty in locating the plume at any given time; second, real dilution would be less than what would be expected
if the macrodispersion values obtained for ergodic conditions were applied; third, the volume that is affected by a non-zero
concentration is smaller than that predicted from macrodispersion concepts; fourth, the orientation of the plume does not
correspond to that of the mean flow; and fifth, accounting for local dispersion helps reducing uncertainty. 相似文献
18.
Matheron and de Marsily [Matheron M, de Marsily G. Is the transport in porous media always diffusive? A counter-example. Water Resour Res 1980;16:901–17] studied transport in a perfectly stratified infinite medium as an idealized aquifer model. They observed superdiffusive solute spreading quantified by anomalous increase of the apparent longitudinal dispersion coefficient with the square root of time. Here, we investigate solute transport in a vertically bounded stratified random medium. Unlike for the infinite medium at asymptotically long times, disorder-induced mixing and spreading is uniquely quantified by a constant Taylor dispersion coefficient. Using a stochastic modeling approach we study the effective mixing and spreading dynamics at pre-asymptotic times in terms of effective average transport coefficients. The latter are defined on the basis of local moments, i.e., moments of the transport Green function. We investigate the impact of the position of the initial plume and the initial plume size on the (highly anomalous) pre-asymptotic effective spreading and mixing dynamics for single realizations and in average. Effectively, the system “remembers” its initial state, the effective transport coefficients show so-called memory effects, which disappear after the solute has sampled the full vertical extent of the medium. We study the impact of the intrinsic non-ergodicity of the confined medium on the validity of the stochastic modeling approach and study in this context the transition from the finite to the infinite medium. 相似文献
19.
We analyze the impact of a linear trend in the mean log-conductivity on the transport of a conservative tracer in a bounded domain. The effects of such a linear trend on solute transport were analyzed in depth for unbounded domains (Rubin and Seong, Water Resour Res 30(11):2901–2911, 1994; Indelman and Rubin, Water Resour Res 31(5):1257–1265, 1995; Water Resour Res 32(5):1257–1265, 1996), whereas studies concerning this special case of medium nonstationarity in finite domains usually focus on head or flow statistics (Guadagnini et al., Stoch Environ Res Risk Assess, 17:394–407, 2003). In this study both ensemble and effective plume moments are provided for an instantaneous release of a solute through a linear source normal to the mean flow direction, by taking into account different sizes of the source. The analysis involving a steady velocity field spatially nonstationary is developed by using the stochastic finite element method. Results show that ensemble moments are affected by increasing trends both parallel and normal to the mean flow direction, but the impact on effective plume moments is very different. A parallel trend does not seem to influence the effective second moments; while a normal trend, although modifies the transverse effective moment only weakly, strongly increases the longitudinal one, especially for large initial sizes of the source. Furthermore, the increase of the particle displacement variance produced by a parallel trend in the finite domain disagrees with the results obtained in an unbounded domain, due to the boundary conditions here considered making both head and velocity moments nonstationary and nonsymmetric. 相似文献
20.
Natural Gradient Tracer Test to Evaluate Natural Attenuation of MTBE Under Anaerobic Conditions 总被引：1，自引：0，他引：1
A natural gradient tracer test using perdeuterated MTBE was conducted in an anaerobic aquifer to determine the relative importance of dispersion and degradation in reducing MTBE concentrations in ground water. Preliminary ground water chemistry and hydraulic conductivity data were used to place the tracer within an existing dissolved MTBE plume at Port Hueneme, California. Following one year of transport, the tracer plume was characterized in detail.
Longitudinal dispersion was identified as the dominant mechanism for lowering the perdeuterated MTBE concentrations. The method of moments was used to determine the longitudinal and lateral dispersion coefficients (0.85 m2 /day and 0.08 m2 /day, respectively). A mass-balance analysis, carried out after one year of transport, accounted for 110% of the injected mass and indicated that no significant mass loss occurred. The plume structure created by zones of higher and lower hydraulic conductivity at the site was complex, consisting of several localized areas of high tracer concentration in a lower concentration plume. This is important because the aquifer has generally been characterized as exhibiting fairly minor heterogeneity. In addition, the tracer plume followed a curved flowpath that deviated from the more macroscopic direction of ground water flow inferred from local ground water elevation measurements and the behavior of the existing plume. Understanding the mass balance, plume structure, curvature of the tracer plume, and consequently natural attenuation behavior required the detailed sampling approach employed in this study. These data imply that a detailed understanding of site hydrogeology and an extensive sampling network may be critical for the correct interpretation of monitored natural attenuation of MTBE. 相似文献
Longitudinal dispersion was identified as the dominant mechanism for lowering the perdeuterated MTBE concentrations. The method of moments was used to determine the longitudinal and lateral dispersion coefficients (0.85 m