共查询到20条相似文献，搜索用时 476 毫秒
1.
2.
3.
4.
Modelling adsorptive solute transport in soils needs a number of parameters to describe its reaction kinetics and the values of these parameters are usually determined from batch and displacement experiments. Some experimental results reveal that when describing the adsorption as first-order kinetics, its associated reaction rates are not constants but vary with pore water velocity. Explanation of this varies but an independent verification of each explanation is difficult because simultaneously measuring the spatiotemporal distributions of dissolved and adsorbed solutes in soils is formidable. Pore-scale modelling could play an important role to address this gap and has received increased attention over the past few years. This paper investigated the transport of adsorptive solute in a simple porous medium using pore-scale modelling. Fluid flow through the void space of the medium was assumed to be laminar and in saturated condition, and solute transport consisted of advection and molecular diffusion; the sorption and desorption occurring at the fluid–solid interface were modelled as linear first-order kinetics. Based on the simulated spatiotemporal distribution of dissolved and adsorbed solutes at pore scale, volumetric-average reaction kinetics at macroscopic scale and its associated reactive parameters were measured. Both homogeneous adsorption where the reaction rates at microscopic scale are constant, and heterogeneous adsorption where the reaction rates vary from site to site, were investigated. The results indicate that, in contrast to previously thought, the macroscopic reaction rates directly measured from the pore-scale simulations do not change with pore velocity under both homogeneous and heterogeneous adsorptions. In particular, we found that for the homogeneous adsorption, the macroscopic adsorption remains first-order kinetic and can be described by constant reaction rates, regardless of flow rate; whilst for the heterogeneous adsorption, the macroscopic adsorption kinetics continues not to be affected by flow rate but is no longer first-order kinetics that can be described by constant reaction rates. We discuss how these findings could help explain some contrary literature reports over the dependence of reaction rates on pore water velocity. 相似文献
5.
Flow and transport simulation in karst aquifers remains a significant challenge for the ground water modeling community. Darcy's law–based models cannot simulate the inertial flows characteristic of many karst aquifers. Eddies in these flows can strongly affect solute transport. The simple two-region conduit/matrix paradigm is inadequate for many purposes because it considers only a capacitance rather than a physical domain. Relatively new lattice Boltzmann methods (LBMs) are capable of solving inertial flows and associated solute transport in geometrically complex domains involving karst conduits and heterogeneous matrix rock. LBMs for flow and transport in heterogeneous porous media, which are needed to make the models applicable to large-scale problems, are still under development. Here we explore aspects of these future LBMs, present simple examples illustrating some of the processes that can be simulated, and compare the results with available analytical solutions. Simulations are contrived to mimic simple capacitance-based two-region models involving conduit (mobile) and matrix (immobile) regions and are compared against the analytical solution. There is a high correlation between LBM simulations and the analytical solution for two different mobile region fractions. In more realistic conduit/matrix simulation, the breakthrough curve showed classic features and the two-region model fit slightly better than the advection-dispersion equation (ADE). An LBM-based anisotropic dispersion solver is applied to simulate breakthrough curves from a heterogeneous porous medium, which fit the ADE solution. Finally, breakthrough from a karst-like system consisting of a conduit with inertial regime flow in a heterogeneous aquifer is compared with the advection-dispersion and two-region analytical solutions. 相似文献
6.
Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based on the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This study provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods. 相似文献
7.
Charles J. Paradis Emma R. Dixon Lauren M. Lui Adam P. Arkin Jack C. Parker Jonathan D. Istok Edmund Perfect Larry D. McKay Terry C. Hazen 《Ground water》2019,57(2):292-302
The breakthrough curve obtained from a single-well push-pull test can be adjusted to account for dilution of the injection fluid in the aquifer fluid. The dilution-adjusted breakthrough curve can be analyzed to estimate the reaction rate of a solute. The conventional dilution-adjusted method assumes that the ratios of the concentrations of the nonreactive and reactive solutes in the injection fluid vs. the aquifer fluid are equal. If this assumption is invalid, the conventional method will generate inaccurate breakthrough curves and may lead to erroneous conclusions regarding the reactivity of a solute. In this study, a new method that generates a dilution-adjusted breakthrough curve was theoretically developed to account for any possible combination of nonreactive and reactive solute concentrations in the injection and aquifer fluids. The newly developed method was applied to a field-based data set and was shown to generate more accurate dilution-adjusted breakthrough curves. The improved dilution-adjusted method presented here is simple, makes no assumptions regarding the concentrations of the nonreactive and reactive solutes in the injection and aquifer fluids, and easily allows for estimating reaction rates during push-pull tests. 相似文献
8.
《Advances in water resources》2005,28(5):479-492
The transient storage model (TSM) has been widely used in studies of stream solute transport and fate, with an increasing emphasis on reactive solute transport. In this study we perform sensitivity analyses of a conservative TSM and two different reactive solute transport models (RSTM), one that includes first-order decay in the stream and the storage zone, and a second that considers sorption of a reactive solute on streambed sediments. Two previously analyzed data sets are examined with a focus on the reliability of these RSTMs in characterizing stream and storage zone solute reactions. Sensitivities of simulations to parameters within and among reaches, parameter coefficients of variation, and correlation coefficients are computed and analyzed. Our results indicate that (1) simulated values have the greatest sensitivity to parameters within the same reach, (2) simulated values are also sensitive to parameters in reaches immediately upstream and downstream (inter-reach sensitivity), (3) simulated values have decreasing sensitivity to parameters in reaches farther downstream, and (4) in-stream reactive solute data provide adequate data to resolve effective storage zone reaction parameters, given the model formulations. Simulations of reactive solutes are shown to be equally sensitive to transport parameters and effective reaction parameters of the model, evidence of the control of physical transport on reactive solute dynamics. Similar to conservative transport analysis, reactive solute simulations appear to be most sensitive to data collected during the rising and falling limb of the concentration breakthrough curve. 相似文献
9.
Pressure pulsing technology is an innovative method that has been developed with the aim of overcoming preferred flow paths associated with remediation techniques that rely on the injection of reagents. Numerical and field experiments were conducted to assess how pressure pulsing affects groundwater flow and solute transport during reagent injection. A series of field experiments were performed at two field sites where a monitoring network designed to capture the breakthrough of solutes delivered from an injection well was installed. Pressure pulsing and conventional injection methods were used at each site. One site was comprised of fine sand with low heterogeneity, while the other was moderately heterogeneous with discrete layers varying from fine sand to silt. The data suggest that breakthrough was more uniform for the pressure pulsing injections; however, this difference was minor and complicated by sorption of some of the tracers employed. The groundwater flow and solute transport modeling exercise simulated the rapid boundary pressure modulation that occurs in association with pressure pulsing. Two‐dimensional (2D) simulations revealed that repeated sudden onset of injection cessation produces brief periods of gradient reversal and the development of a mixing zone near the injection well. The spatial extents of this mixing zone were found to be highly dependent upon the hydraulic diffusivity of the medium, with medium heterogeneity and pulsing frequency playing secondary roles. Three‐dimensional (3D) numerical simulations were used to benchmark the observations from one of the field sites. The results from the modeling effort showed that solute breakthrough from a pressure pulsing injection is more dispersed relative to a conventional injection as a result of the mixing zone phenomenon; however, we were unable to directly observe this mixing zone using the instrumentation deployed at the two field sites. 相似文献
10.
11.
《Advances in water resources》2007,30(6-7):1593-1607
An experimental and numerical investigation was conducted to study the colonization dynamics of a bioluminescent bacterium, Pseudomonas fluorescens HK44, during growth in a porous medium under steady, variably saturated flow conditions. Experiments were conducted in a thin-slab light transmission chamber filled with uniform, translucent quartz sand. Steady, variably saturated flow conditions were established using drip emitters mounted on the top of the chamber, with glucose applied through a central dripper located directly above an inoculated region of the chamber. Periodic pulses of salicylate and a dye tracer were applied to induce bioluminescence of the bacterium to monitor colony expansion and to track changes in the hydraulic and transport properties of the sand. Changes in the apparent water saturation of the sand were quantified by monitoring light transmission through the chamber with a CCD camera. The colonized region expanded laterally by about 15 cm, and upward against the flow by 7–8 cm during the 6-day experiment while apparent saturations in the colonized region decreased by 7–9% and the capillary fringe dropped by ∼5 cm. The observed data were reproduced approximately using a numerical model that accounted for the processes of water flow, solute and bacterial transport, cell growth and accumulation, glucose and oxygen consumption, and gas diffusion and exchange. The results of this study illustrate some of the complexities associated with coupled flow, reactive transport, and biological processes in variably saturated porous media, such as localized desaturation, capillary fringe lowering effects, and upstream movement of bacterial colonization, that may not readily observable using other experimental techniques. 相似文献
12.
We used the 3D continuum-scale reactive transport models to simulate eight core flood experiments for two different carbonate rocks. In these experiments the core samples were reacted with brines equilibrated with pCO2 = 3, 2, 1, 0.5 MPa (Smith et al., 2013 [27]). The carbonate rocks were from specific Marly dolostone and Vuggy limestone flow units at the IEAGHG Weyburn-Midale CO2 Monitoring and Storage Project in south-eastern Saskatchewan, Canada. Initial model porosity, permeability, mineral, and surface area distributions were constructed from micro tomography and microscopy characterization data. We constrained model reaction kinetics and porosity–permeability equations with the experimental data. The experimental data included time-dependent solution chemistry and differential pressure measured across the core, and the initial and final pore space and mineral distribution. Calibration of the model with the experimental data allowed investigation of effects of carbonate reactivity, flow velocity, effective permeability, and time on the development and consequences of stable and unstable dissolution fronts.The continuum scale model captured the evolution of distinct dissolution fronts that developed as a consequence of carbonate mineral dissolution and pore scale transport properties. The results show that initial heterogeneity and porosity contrast control the development of the dissolution fronts in these highly reactive systems. This finding is consistent with linear stability analysis and the known positive feedback between mineral dissolution and fluid flow in carbonate formations. Differences in the carbonate kinetic drivers resulting from the range of pCO2 used in the experiments and the different proportions of more reactive calcite and less reactive dolomite contributed to the development of new pore space, but not to the type of dissolution fronts observed for the two different rock types. The development of the dissolution front was much more dependent on the physical heterogeneity of the carbonate rock. The observed stable dissolution fronts with small but visible dissolution fingers were a consequence of the clustering of a small percentage of larger pores in an otherwise homogeneous Marly dolostone. The observed wormholes in the heterogeneous Vuggy limestone initiated and developed in areas of greater porosity and permeability contrast, following pre-existing preferential flow paths.Model calibration of core flood experiments is one way to specifically constrain parameter input used for specific sites for larger scale simulations. Calibration of the governing rate equations and constants for Vuggy limestones showed that dissolution rate constants reasonably agree with published values. However the calcite dissolution rate constants fitted to the Marly dolostone experiments are much lower than those suggested by literature. The differences in fitted calcite rate constants between the two rock types reflect uncertainty associated with measured reactive surface area and appropriately scaling heterogeneous distribution of less abundant reactive minerals. Calibration of the power-law based porosity–permeability equations was sensitive to the overall heterogeneity of the cores. Stable dissolution fronts of the more homogeneous Marly dolostone could be fit with the exponent n = 3 consistent with the traditional Kozeny–Carman equation developed for porous sandstones. More impermeable and heterogeneous cores required larger n values (n = 6–8). 相似文献
13.
Chaotic advection at the pore scale: Mechanisms,upscaling and implications for macroscopic transport
The macroscopic spreading and mixing of solute plumes in saturated porous media is ultimately controlled by processes operating at the pore scale. Whilst the conventional picture of pore-scale mechanical dispersion and molecular diffusion leading to persistent hydrodynamic dispersion is well accepted, this paradigm is inherently two-dimensional (2D) in nature and neglects important three-dimensional (3D) phenomena. We discuss how the kinematics of steady 3D flow at the pore scale generate chaotic advection—involving exponential stretching and folding of fluid elements—the mechanisms by which it arises and implications of microscopic chaos for macroscopic dispersion and mixing. Prohibited in steady 2D flow due to topological constraints, these phenomena are ubiquitous due to the topological complexity inherent to all 3D porous media. Consequently 3D porous media flows generate profoundly different fluid deformation and mixing processes to those of 2D flow. The interplay of chaotic advection and broad transit time distributions can be incorporated into a continuous-time random walk (CTRW) framework to predict macroscopic solute mixing and spreading. We show how these results may be generalised to real porous architectures via a CTRW model of fluid deformation, leading to stochastic models of macroscopic dispersion and mixing which both honour the pore-scale kinematics and are directly conditioned on the pore-scale architecture. 相似文献
14.
Thomas Willingham Changyong Zhang Charles J. Werth Albert J. Valocchi Mart Oostrom Thomas W. Wietsma 《Advances in water resources》2010
A key challenge for predictive modeling of transverse mixing and reaction of solutes in groundwater is to determine values of transverse dispersivity (αT) in heterogeneous flow fields that accurately describe mixing and reaction at the pore scale. We evaluated the effects of flow focusing in high permeability zones on mixing enhancement using experimental micromodel flow cells and pore-scale lattice-Boltzmann-finite-volume model (LB-FVM) simulations. Micromodel results were directly compared to LB-FVM simulations using two different pore structures, and excellent agreement was obtained. Six different flow focusing pore structures were then systematically tested using LB-FVM, and both analytical solutions and a two-dimensional (2D) continuum-scale model were used to fit αT values to pore-scale results. Pore-scale results indicate that the overall rate of mixing-limited reaction increased by up to 40% when flow focusing occurred, and it was greater in pore structures with longer flow focusing regions and greater porosity contrast. For each pore structure, αT values from analytical solutions of transverse concentration profiles or total product at a given longitudinal location showed good agreement for nonreactive and reactive solutes, and values determined in flow focusing zones were always smaller than those downgradient after the flow focusing zone. Transverse dispersivity values from the 2D continuum model were between values within and downgradient from the flow focusing zone determined from analytical solutions. Also, total product and transverse concentration profiles along the entire pore structure from the 2D continuum model matched pore scale results. These results indicate that accurate quantification of pore-scale flow focusing with transverse dispersion coefficients is possible only when the entire flow and concentration fields are considered. 相似文献
15.
We review the analysis of the dynamics of reactive transport in disordered media, emphasizing the nature of the chemical reactions and the role of small-scale fluctuations induced by the structure of the porous medium. We are motivated by results and interpretations of laboratory-scale experiments, for which detailed characterization of the system is possible. Modeling approaches based on continuum and particle tracking (PT) schemes are examined critically, highlighting how fluctuations are incorporated. The continuum approach spans a large literature. Traditional formats of reactive transport equations, such as the advection–dispersion–reaction equation (ADRE), are based on a series of assumptions related mainly to scale separation and relative magnitude of time scales involved in the reactive transport setting. These assumptions as well as further developments are assessed in depth. PT methods offer an alternative means of accounting for pore-scale dynamics, wherein space–time transitions are drawn from appropriate probability distributions that have been tested to account for anomalous transport. While PT methods have been employed for many years to describe conservative transport, their application to laboratory-scale reactive transport problems in the context of both Fickian and non-Fickian regimes is relatively recent. We concentrate on experimental observations of different types of reactions in disordered media: (1) the dynamics of a bimolecular reactive transport (A + B → C) in passive (non-reactive) media, and (2) a multi-step chemical reaction, as exemplified in the process of dedolomitization involving both dissolution and precipitation. The fluctuations in a number of the key variables controlling the processes prove to have a dominant role; elucidation of this role forms the basis of the present study and the comparison of methods. 相似文献
16.
The finite-element method based on a Galerkin technique was used to formulate the problem of simulating the two-dimensional (cross-sectional) transient movement of water and solute in saturated or partially saturated nonuniform porous media. The numerical model utilizes linear triangular elements. Nonreactive, as well as reactive solutes whose behaviour can be described by a distribution coefficient or first-order reaction term were considered. The flow portion of the model was tested by comparison of the model results with experimental and finite-difference results for transient flow in an unsaturated sand column and the solute transport portion of the model was tested by comparison with analytical solution results. The model was applied to a hypothetical case involving movement of water and solutes in tile-drained soils. The simulation results showed the development of distinct solute leaching patterns in the soil as drainage proceeded. Although applied to a tile drainage problem in this study, the model should be equally useful in the study of a wide range of two-dimensional water and solute migration problems. 相似文献
17.
《Advances in water resources》2007,30(6-7):1583-1592
A major problem preventing widespread implementation of microbial injection strategies for bioremediation and/or microbially enhanced oil recovery is the tendency of bacteria to strongly adhere to surfaces in the immediate vicinity of the injection point. Long term (weeks to months) nutrient starvation of bacteria prior to injection can decrease attachment and enhance transport through porous media. This paper summarizes results of starvation-enhanced transport experiments in sand columns of 30 cm, 3 m, and 16 m in length. The 16 m column experiments compared transport, breakthrough and distribution of adhered cells for starved and vegetative cultures of Klebsiella oxytoca, a copious biofilm producer. Results from these experiments were subsequently used to design and construct a field-scale biofilm barrier using starved Pseudomonas fluorescens. The 30 cm and 3 m sand columns experiments investigated starvation-enhanced transport of Shewanella algae BrY, a dissimilatory metal-reducing bacterium. In both cases the vegetative cells adsorbed onto the sand in higher numbers than the starved cells, especially near the entrance of the column. These results, taken together with studies cited in the literature, indicate that starved cells penetrate farther (i.e. higher breakthrough concentration) and adsorb more uniformly along the flow path than vegetative cells. 相似文献
18.
This study formulates and analyzes continuous time random walk (CTRW) models in radial flow geometries for the quantification of non-local solute transport induced by heterogeneous flow distributions and by mobile–immobile mass transfer processes. To this end we derive a general CTRW framework in radial coordinates starting from the random walk equations for radial particle positions and times. The particle density, or solute concentration is governed by a non-local radial advection–dispersion equation (ADE). Unlike in CTRWs for uniform flow scenarios, particle transition times here depend on the radial particle position, which renders the CTRW non-stationary. As a consequence, the memory kernel characterizing the non-local ADE, is radially dependent. Based on this general formulation, we derive radial CTRW implementations that (i) emulate non-local radial transport due to heterogeneous advection, (ii) model multirate mass transfer (MRMT) between mobile and immobile continua, and (iii) quantify both heterogeneous advection in a mobile region and mass transfer between mobile and immobile regions. The expected solute breakthrough behavior is studied using numerical random walk particle tracking simulations. This behavior is analyzed by explicit analytical expressions for the asymptotic solute breakthrough curves. We observe clear power-law tails of the solute breakthrough for broad (power-law) distributions of particle transit times (heterogeneous advection) and particle trapping times (MRMT model). The combined model displays two distinct time regimes. An intermediate regime, in which the solute breakthrough is dominated by the particle transit times in the mobile zones, and a late time regime that is governed by the distribution of particle trapping times in immobile zones. These radial CTRW formulations allow for the identification of heterogeneous advection and mobile-immobile processes as drivers of anomalous transport, under conditions relevant for field tracer tests. 相似文献
19.
《Advances in water resources》2007,30(4):838-850
Treatment wetlands play an important role in reducing nutrient content and heavy metals in wastewater and run-off water. The treatment efficiency strongly depends on flow pattern and residence times of the water. Here, we study the impact of different factors on water flow patterns based on a tracer experiment with tritiated water in a 2.6 ha constructed wetland pond. A 2D flow and inert transport model was used to evaluate the relative importance of bottom topography, vegetation distribution, water exchange with stagnant zones and dispersion.Results from computer simulations and independent measurements of friction losses as well as wetland geometry showed that variations in bottom topography, formed by several deep zones, decreased the variance in water residence times to a minor extent. Heterogeneity in vegetation, on the other hand, significantly contributed to the spread in water residence times and explained the multiple peaks observed in the breakthrough curves. Analyses showed that in the Ekeby treatment wetland, basin shape explained about 10% of the variance in the observed residence times, whereas vegetation explained about 60–80%. To explain all variance secondary factors were needed, such as dispersion and water exchange with stagnant zones. These were shown to contribute to the spread of residence times and primarily to the long tail of the observed breakthrough curves. 相似文献
20.
Kaili Wang Guanhua Huang 《Advances in water resources》2011,34(6):671-683
The effect of aquifer heterogeneity on flow and solute transport in two-dimensional isotropic porous media was analyzed using the Monte Carlo method. The two-dimensional logarithmic permeability (ln K) was assumed to be a non-stationary random field with its increments being a truncated fractional Lévy motion (fLm). The permeability fields were generated using the modified successive random additions (SRA) algorithm code SRA3DC [http://www.iamg.org/CGEditor/index.htm]. The velocity and concentration fields were computed respectively for two-dimensional flow and transport with a pulse input using the finite difference codes of MODFLOW 2000 and MT3DMS. Two fLm control parameters, namely the width parameter (C) and the Lévy index (α), were varied systematically to examine their effect on the resulting permeability, flow velocity and concentration fields. We also computed the first- and second-spatial moments, the dilution index, as well as the breakthrough curves at different control planes with the corresponding concentration fields. In addition, the derived breakthrough curves were fitted using the continuous time random walk (CTRW) and the traditional advection-dispersion equation (ADE). Results indicated that larger C and smaller α both led to more heterogeneous permeability and velocity fields. The Lévy-stable distribution of increments in ln K resulted in a Lévy-stable distribution of increments in logarithm of the velocity (ln v). Both larger C and smaller α created sharper leading edges and wider tailing edges of solute plumes. Furthermore, a relatively larger amount of solute still remained in the domain after a relatively longer time transport for smaller α values. The dilution indices were smaller than unity and increased as C increased and α decreased. The solute plume and its second-spatial moments increased as C increased and α decreased, while the first-spatial moments of the solute plume were independent of C and α values. The longitudinal macrodispersivity was scale-dependent and increased as a power law function of time. Increasing C and decreasing α both resulted in an increase in longitudinal macrodispersivity. The transport in such highly heterogeneous media was slightly non-Gaussian with its derived breakthrough curves being slightly better fitted by the CTRW than the ADE, especially in the early arrivals and late-time tails. 相似文献