首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
We present a critical analysis of experimental findings on vegetation–flow–sediment interactions obtained through both laboratory and field experiments on tidal and coastal environments. It is well established that aquatic vegetation provides a wide range of ecosystem services (e.g. protecting coastal communities from extreme events, reducing riverbank and coastal erosion, housing diverse ecosystems), and the effort to better understand such services has led to multiple approaches to reproduce the relevant physical processes through detailed laboratory experiments. State-of-the-art measurement techniques allow researchers to measure velocity fields and sediment transport with high spatial and temporal resolution under well-controlled flow conditions, yielding predictions for hydrodynamic and sediment transport scenarios that depend on simplified or bulk vegetation parameters. However, recent field studies have shown that some simplifications on the experimental setup (e.g. the use of rigid elements, a single diameter, a single element height, regular or staggered layout) can bias the outcome of the study, by either hiding or amplifying some of the relevant physical processes found in natural conditions. We discuss some observed cases of bias, including general practices that can lead to compromises associated with simplified assumptions. The analysis presented will identify potential pathways to move forward with laboratory and field measurements, which could better inform predictors to produce more robust, universal and accurate predictions on flow–vegetation–sediment interactions. © 2020 John Wiley & Sons, Ltd.  相似文献   

2.
Bacteriophages have been used in soil column studies for the last several decades as surrogates to study the fate and transport behavior of enteric viruses in groundwater. However, recent studies have shown that the transport behavior of bacteriophages and enteric viruses in porous media can be very different. The next generation of virus transport science must therefore provide more data on mobility of enteric viruses and the relationship between transport behaviors of enteric viruses and bacteriophages. To achieve this new paradigm, labor intensity devoted to enteric virus quantification method must be reduced. Recent studies applied quantitative polymerase chain reaction (qPCR) to column filtration experiments to study the transport behavior of human adenovirus (HAdV) in porous media under a variety of conditions. A similar approach can be used to study the transport of other enteric viruses such as norovirus. Analyzing the column samples with both qPCR and culture assays and applying multiplex qPCR to study cotransport behavior of more than one virus will provide information to under‐explored areas in virus transport science. Both nucleic acid extraction kits and one‐step lysis protocols have been used in these column studies to extract viral nucleic acid for qPCR quantification. The pros and cons of both methods are compared herein and solutions for overcoming problems are suggested. As better understanding of the transport behavior of enteric viruses is clearly needed, we strongly advocate for application of rapid molecular tools in future studies as well as optimization of protocols to overcome their current limitations.  相似文献   

3.
4.
Abstract

Many advances in simulating single and two-phase fluid flow and heat transport in porous media have recently been made in conjunction with geothermal energy research. These numerical models reproduce system thermal and pressure behaviour and can be used for heat-transport problems other than those associated with geothermal energy development, such as high-level radioactive waste disposal and heat-storage projects. Although these models are general, additional research is necessary before they can be applied to certain site-specific problems that are concerned with additional processes, such as mass transport and flow in fractured media.  相似文献   

5.
Qinhong Hu  Jean E. Moran 《水文研究》2005,19(14):2671-2687
An analytical method that employs ion chromatography has been developed to exploit the use of fluorobenzoic acids (FBAs) and halides more fully as hydrologic tracers. In a single run, this reliable, sensitive, and robust method can simultaneously separate and quantify halides (fluoride, chloride, bromide, and iodide) and up to seven FBAs from other common groundwater constituents (e.g. nitrate and sulphate). The usefulness of this analytical method is demonstrated in both field and laboratory tracer experiments. The field study examines the hydrologic response of fractures and the matrix to different flow rates and the contribution of matrix diffusion in chemical transport. Laboratory tracer experiments with eight geologic media from across the USA—mostly from Department of Energy facilities where groundwater contamination is prevalent and where subsurface characterization employing tracers has been ongoing or is in need—reveal several insights about tracer transport behaviour: (1) bromide and FBAs are not always transported conservatively; (2) the delayed transport of these anionic tracers is likely related to geologic media characteristics, such as organic matter, pH, iron oxide content, and clay mineralogy; (3) use of iodine as a hydrologic tracer should take into account the different sorption behaviours of iodide and iodate and the possible conversion of iodine's initial chemical form; (4) the transport behaviour of potential FBA and halide tracers under relevant geochemical conditions should be evaluated before beginning ambitious, large‐scale field tracer experiments. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

6.
In a flume experiment with steady flow conditions, H. A. Einstein recognised the transport of bedload particles as consisting of steps of rolling, sliding, or saltation with intermittent rest periods, and introduced the concept of an average, ‘virtual’ transport velocity. This virtual velocity then has also been derived from tracer studies in the field by dividing the travelled distance of a tracer by the duration of competent flow. As a consequence, the virtual velocity in the field is represented by one single value only, despite the unsteady flow variables. Tracer measurements in a river have not been yet used to express transport velocity as a direct function of these actual variables, and insights from tracer measurements into the processes of sediment transport remain limited. In particular, the unsteady conditions for bedload in the field have impeded the derivation of sediment transport characteristics as determined from laboratory experiments, as well as the transfer of laboratory insights to a field setting. We introduce a method of data regression for the derivation of an ‘unsteady’ virtual velocity from repeated surveys of tracer positions. The regression program called graVel (provided as supplementary material) relates the integral of an excess flow variable term to measured travel distances, yielding the most probable threshold value for entrainment and the coefficient of linear and non‐linear formulas. An extended regression allows additional fitting of the exponent in non‐linear formulas. Application to published tracer data from the Mameyes River, Puerto Rico, shows that the unsteady virtual velocity is more likely governed by non‐linear relations to excess Shields stress, similar to bedload transport, than by relations linking the particle velocity linearly to excess shear velocity. Partial agreements with non‐dimensional results derived from the larger, non‐wadeable Mur River encourage the establishment of a generalised formula for the unsteady virtual velocity. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

7.
A new lumped-parameter approach to simulating unsaturated flow processes in dual-porosity media such as fractured rocks or aggregated soils is presented. Fluid flow between the fracture network and the matrix blocks is described by a non-linear equation that relates the imbibition rate of the local difference in liquid-phase pressure between the fractures and the matrix blocks. Unlike a Warren-Root-type equation, this equation is accurate in both the early and late time regimes. The fracture/matrix interflow equation has been incorporated into an existing unsaturated flow simulator, to serve as a source/sink term for fracture gridblocks. Flow processes are then simulated using only fracture gridblocks in the computational grid. This new lumped-parameter approach has been tested on two problems involving transient flow in fractured/porous media, and compared with simulations performed using explicit discretisation of the matrix blocks. The new procedure seems to accurately simulate flow processes in unsaturated fractured rocks, and typically requires an order of magnitude less computational time than do simulations using fully-discretised matrix blocks.  相似文献   

8.
Infiltration into frozen soil plays an important role in soil freeze–thaw and snowmelt-driven hydrological processes. To better understand the complex thermal energy and water transport mechanisms involved, the influence of antecedent moisture content and macroporosity on infiltration into frozen soil was investigated. Ponded infiltration experiments on frozen macroporous and non-macroporous soil columns revealed that dry macroporous soil produced infiltration rates reaching 103 to 104 mm day−1, two to three orders of magnitude larger than dry non-macroporous soil. Results suggest that rapid infiltration and drainage were a result of preferential flow through initially air-filled macropores. Using recorded flow rates and measured macropore characteristics, calculations indicated that a combination of both saturated flow and unsaturated film flow likely occurred within macropores. Under wet conditions, regardless of the presence of macropores, infiltration was restricted by the slow thawing rate of pore ice, producing infiltration rates of 2.8 to 5.0 mm day−1. Reduced preferential flow under wet conditions was attributed to a combination of soil swelling, due to smectite-rich clay (that reduced macropore volume), and pore ice blockage within macropores. In comparison, dry soil column experiments demonstrated that macropores provided conduits for water and thermal energy to bypass the frozen matrix during infiltration, reducing thaw rates compared with non-macroporous soils. Overall, results showed the dominant control of antecedent moisture content on the initiation, timing, and magnitude of infiltration and flow in frozen macroporous soils, as well as the important role of macropore connectivity. The study provides an important data set that can aid the development of hydrological models that consider the interacting effects of soil freeze–thaw and preferential flow on snowmelt partitioning in cold regions.  相似文献   

9.
Tracer experiments conducted in the laboratory on undisturbed core samples (<7.3-cm-diameter) have been a standard method for estimating hydraulic and transport properties of fractured till since the 1980s. This study assesses the relationship between visible fractures on the top and bottom of core samples and the resulting hydraulic and mass transport properties of the core. We hypothesized that more visible fractures would indicate the presence of a well-connected fracture network, leading to greater hydraulic conductivity (K) values and earlier chemical breakthrough times. To test this hypothesis, water flow and bromide (Br-) tracer experiments were performed on 10, 16-cm diameter, 16-cm-tall samples of fractured Dows Formation till from central Iowa. Visually identifiable fractures were present on the top and bottom of every sample. Results indicate that the visual identification of fractures does not predict a connected fracture network, as some samples produced breakthrough curves showing rapid first arrival times and shapes characteristic of solute transport in a fractured medium, while others appeared similar to an unfractured medium. No correlation was found between the number of visible fractures and K (Pearson's r = 0.25), or Br- first arrival time (r = −0.33), but a strong negative correlation between K and first arrival time (r = −0.92). Results indicate that the sample volume was not large enough to reliably contain a connected fracture network. Thus, testing large volumes of till at the field scale coupled with fracture-flow modeling likely represents the best approach for estimating hydraulic and mass transport properties for fractured till.  相似文献   

10.
Numerical simulation of flows in shallow reservoirs has to be checked for its consistency in predicting real flow conditions and sedimentation patterns. Typical flow patterns may exhibit flow separation at the inlet, accompanied by several recirculation and stagnation areas all over the reservoir surface. The aim of the present research project is to study the influence of the geometry of a reservoir on sediment transport and deposition numerically and experimentally, focusing on a prototype reservoir depth between 5 and 15 m as well as suspended sediment transport.
A series of numerical simulations is presented and compared with scaled laboratory experiments, with the objective of testing the sensitivity to different flow and sediment parameters and different turbulence closure schemes. Different scenarios are analyzed and a detailed comparison of preliminary laboratory tests and some selected simulations are presented.
The laboratory experiments show that suspended sediment transport and deposition are determined by the initial flow pattern and by the upstream and downstream boundary conditions. In the experiments, deposition in the rectangular basin systematically developed along the left bank, although inflow and outflow were positioned symmetrically along the centre of the basin. Three major horizontal eddies developed influencing the sediment deposition pattern. Although asymmetric flow patterns are privileged, a symmetric pattern can appear from time to time. This particular behaviour could also be reproduced by a two-dimensional depth-averaged flow and sediment transport model (CCHE2D). The paper presents numerical simulations using different turbulence closure schemes (k-ε and eddy viscosity models). In spite of the symmetric setup, these generally produced an asymmetric flow pattern that can easily switch sides depending on the assumptions made for the initial and boundary conditions. When using the laboratory experiment as a reference, the most reliable numerical results have been obtai  相似文献   

11.
An Eulerian perturbation approach was applied to develop a method of moment for solute transport in a nonstationary, fractured medium. The conceptualized fractured medium is described through a dual-porosity model. Stochastic governing equations for mean concentration and concentration covariance were analytically derived to the first-order accuracy of log-conductivity variance and solved with a numerical method––a finite difference method. The developed method is called a numerical Eulerian method of moment (NEMM). This method was compared with the stationary transport theory [Water Resour. Res. 36(7) (2000) 1665] for predicting mean concentration and its spatial moments. The comparison indicated that the two methods matched very well in predicting first and second spatial moments. NEMM solutions were also compared with Monte Carlo simulations for solute transport in stationary fractured media. The results of the two methods were consistent for calculating small log conductivity variance. The theory was then used to study effects of various parameters and nonstationarity of the medium on flow and transport processes. Results indicated that medium nonstationarity would significantly influence the solute transport process. The nonstationary transport theory relaxes many assumptions adopted in stationary theories and paves the way for applying the NEMM to many environmental projects, especially in analyzing uncertainty of solute transport.  相似文献   

12.
A three-dimensional, reactive numerical flow model is developed that couples chemical reactions with density-dependent mass transport and fluid flow. The model includes equilibrium reactions for the aqueous species, kinetic reactions between the solid and aqueous phases, and full coupling of porosity and permeability changes that result from precipitation and dissolution reactions in porous media. A one-step, global implicit approach is used to solve the coupled flow, transport and reaction equations with a fully implicit upstream-weighted control volume discretization. The Newton–Raphson method is applied to the discretized non-linear equations and a block ILU-preconditioned CGSTAB method is used to solve the resulting Jacobian matrix equations. This approach permits the solution of the complete set of governing equations for both concentration and pressure simultaneously affected by chemical and physical processes. A series of chemical transport simulations are conducted to investigate coupled processes of reactive chemical transport and density-dependent flow and their subsequent impact on the development of preferential flow paths in porous media. The coupled effects of the processes driving flow and the chemical reactions occurring during solute transport is studied using a carbonate system in fully saturated porous media. Results demonstrate that instability development is sensitive to the initial perturbation caused by density differences between the solute plume and the ambient groundwater. If the initial perturbation is large, then it acts as a “trigger” in the flow system that causes instabilities to develop in a planar reaction front. When permeability changes occur due to dissolution reactions occurring in the porous media, a reactive feedback loop is created by calcite dissolution and the mixed convective transport of the system. Although the feedback loop does not have a significant impact on plume shape, complex concentration distributions develop as a result of the instabilities generated in the flow system.  相似文献   

13.
Phast4Windows is a Windows® program for developing and running groundwater‐flow and reactive‐transport models with the PHAST simulator. This graphical user interface allows definition of grid‐independent spatial distributions of model properties—the porous media properties, the initial head and chemistry conditions, boundary conditions, and locations of wells, rivers, drains, and accounting zones—and other parameters necessary for a simulation. Spatial data can be defined without reference to a grid by drawing, by point‐by‐point definitions, or by importing files, including ArcInfo® shape and raster files. All definitions can be inspected, edited, deleted, moved, copied, and switched from hidden to visible through the data tree of the interface. Model features are visualized in the main panel of the interface, so that it is possible to zoom, pan, and rotate features in three dimensions (3D). PHAST simulates single phase, constant density, saturated groundwater flow under confined or unconfined conditions. Reactions among multiple solutes include mineral equilibria, cation exchange, surface complexation, solid solutions, and general kinetic reactions. The interface can be used to develop and run simple or complex models, and is ideal for use in the classroom, for analysis of laboratory column experiments, and for development of field‐scale simulations of geochemical processes and contaminant transport.  相似文献   

14.
Fractures in porous media have been documented extensively. However, they are often omitted from groundwater flow and mass transport models due to a lack of data on fracture hydraulic properties and the computational burden of simulating fractures explicitly in large model domains. We present a MATLAB toolbox, FracKfinder, that automates HydroGeoSphere (HGS), a variably saturated, control volume finite-element model, to simulate an ensemble of discrete fracture network (DFN) flow experiments on a single cubic model mesh containing a stochastically generated fracture network. Because DFN simulations in HGS can simulate flow in both a porous media and a fracture domain, this toolbox computes tensors for both the matrix and fractures of a porous medium. Each model in the ensemble represents a different orientation of the hydraulic gradient, thus minimizing the likelihood that a single hydraulic gradient orientation will dominate the tensor computation. Linear regression on matrices containing the computed three-dimensional hydraulic conductivity (K) values from each rotation of the hydraulic gradient is used to compute the K tensors. This approach shows that the hydraulic behavior of fracture networks can be simulated where fracture hydraulic data are limited. Simulation of a bromide tracer experiment using K tensors computed with FracKfinder in HGS demonstrates good agreement with a previous large-column, laboratory study. The toolbox provides a potential pathway to upscale groundwater flow and mass transport processes in fractured media to larger scales.  相似文献   

15.
The delineation of wellhead protection areas (WHPAs) under uncertainty is still a challenge for heterogeneous porous media. For granular media, one option is to combine particle tracking (PT) with the Monte Carlo approach (PT‐MC) to account for geologic uncertainties. Fractured porous media, however, require certain restrictive assumptions under this approach. An alternative for all types of media is the capture probability (CP) approach, which is based on the solution of the standard advection‐dispersion equation in a backward mode, making use of the analogy between forward and backward transport processes. Within this context, we review the current controversy about the correct form of the conceptual model for transport, finding that the advection‐diffusion model, which represents the diffusive interchange between streamtubes with differing velocities, is more physically realistic than the conventional advection‐dispersion model. For mildly to moderately heterogeneous materials, stochastic theories and simulation experiments show that this process converges at the field scale to an effective advection‐dispersion process that can be simulated with conventional transport models using appropriate macrodispersivity values. For highly heterogeneous materials, stochastic theories do not yet exist but there is no reason why the process should not converge naturally as well. Macrodispersivities appear to be formation‐specific. The advection‐dispersion model can be used for capture zone delineation in heterogeneous granular media. For fractured porous systems, hybrid equivalent porous medium and discrete fracture network or CP‐based approaches may have potential. In general, capture zones delineated by PT without MC will always be too small and should not be used as a basis for land‐use decisions.  相似文献   

16.
Quantifying preferential flow in soils: A review of different techniques   总被引:19,自引:0,他引:19  
Preferential flow (PF) in soil has both environmental and human health implications since it favours contaminant transport to groundwater without interaction with the chemically and biologically reactive upper layer of soil. PF is, however, difficult to measure and quantify. This paper reviews laboratory and field techniques, such as breakthrough curves, dye tracing, and scanning techniques, for evaluating PF in soil at different scales. Advanced technologies, such as scanning techniques, have increased our capability to quantify transport processes within the soil with minimal soil disturbance. Important issues with respect to quantifying PF concern large-scale studies, frozen soil conditions, tracing techniques for particles and gases, a lack of simple mathematical tools for interpreting field data, and the lack of a systematic approach for comparing PF data resulting from different measurement techniques. Also, more research is required to quantify the relative importance of the various PF processes that occur in soil rather than the integrated result of all PF processes in soils.  相似文献   

17.
Discrete-fracture and dual-porosity models are infrequently used to simulate solute transport through fractured unconsolidated deposits, despite their more common application in fractured rock where distinct flow regimes are hypothesized. In this study, we apply four fracture transport models--the mobile-immobile model (MIM), parallel-plate discrete-fracture model (PDFM), and stochastic and deterministic discrete-fracture models (DFMs)--to demonstrate their utility for simulating solute transport through fractured till. Model results were compared to breakthrough curves (BTCs) for the conservative tracers potassium bromide (KBr), pentafluorobenzoic acid (PFBA), and 1,4-piperazinediethanesulfonic acid (PIPES) in a large-diameter column of fractured till. Input parameters were determined from independent field and laboratory methods. Predictions of Br BTCs were not significantly different among models; however, the stochastic and deterministic DFMs were more accurate than the MIM or PDFM when predicting PFBA and PIPES BTCs. DFMs may be more applicable than the MIM for tracers with small effective diffusion coefficients (De) or for short timescales due to differences in how these models simulate diffusion or incorporate heterogeneities by their fracture networks. At large scales of investigation, the more computationally efficient MIM and PDFM may be more practical to implement than the three-dimensional DFMs, or a combination of model approaches could be employed. Regardless of the modeling approach used, fractures should be incorporated routinely into solute transport models in glaciated terrain.  相似文献   

18.
Diffusive mass exchange into immobile water regions within heterogeneous porous aquifers influences the fate of solutes. The percentage of immobile water is often unidentified in natural aquifers though. Hence, the mathematical prediction of solute transport in such heterogeneous aquifers remains challenging. The objective of this study was to find a simple analytical model approach that allows quantifying properties of mobile and immobile water regions and the portion of immobile water in a porous system. Therefore, the Single Fissure Dispersion Model (SFDM), which takes into account diffusive mass exchange between mobile and immobile water zones, was applied to model transport in well‐defined saturated dual‐porosity column experiments. Direct and indirect model validation was performed by running experiments at different flow velocities and using conservative tracer with different molecular diffusion coefficients. In another column setup, immobile water regions were randomly distributed to test the model applicability and to determine the portion of immobile water. In all setups, the tracer concentration curves showed differences in normalized maximum peak concentration, tailing and mass recovery according to their diffusion coefficients. These findings were more pronounced at lower flow rates (larger flow times) indicating the dependency of diffusive mass exchange into immobile water regions on tracers' molecular diffusion coefficients. The SFDM simulated all data with high model efficiency. Successful model validation supported the physical meaning of fitted model parameters. This study showed that the SFDM, developed for fissured aquifers, is applicable in porous media and can be used to determine porosity and volume of regions with immobile water. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

19.
One of the more advanced approaches for simulating groundwater flow in fractured porous media is the discrete-fracture approach. This approach is limited by the large computational overheads associated with traditional modeling methods. In this work, we apply the Lanczos reduction method to the modeling of groundwater flow in fractured porous media using the discrete-fracture approach. The Lanczos reduction method reduces a finite element equation system to a much smaller tridiagonal system of first-order differential equations. The reduced system can be solved by a standard tridiagonal algorithm with little computational effort. Because solving the reduced system is more efficient compared to solving the original system, the simulation of groundwater flow in discretely fractured media using the reduction method is very efficient. The proposed method is especially suitable for the problem of large-scale and long-term simulation. In this paper, we develop an iterative version of Lanczos algorithm, in which the preconditioned conjugate gradient solver based on ORTHOMIN acceleration is employed within the Lanczos reduction process. Additional efficiency for the Lanczos method is achieved by applying an eigenvalue shift technique. The “shift” method can improve the Lanczos system convergence, by requiring fewer modes to achieve the same level of accuracy over the unshifted case. The developed model is verified by comparison with dual-porosity approach. The efficiency and accuracy of the method are demonstrated on a field-scale problem and compared to the performance of classic time marching method using an iterative solver on the original system. In spite of the advances, more theoretical work needs to be carried out to determine the optimal value of the shift before computations are actually carried out.  相似文献   

20.
This review focuses on investigations of groundwater flow and solute transport in karst aquifers through laboratory scale models (LSMs). In particular, LSMs have been used to generate new data under different hydraulic and contaminant transport conditions, testing of new approaches for site characterization, and providing new insights into flow and transport processes through complex karst aquifers. Due to the increasing need for LSMs to investigate a wide range of issues, associated with flow and solute migration karst aquifers this review attempts to classify, and introduce a framework for constructing a karst aquifer physical model that is more representative of field conditions. The LSMs are categorized into four groups: sand box, rock block, pipe/fracture network, and pipe-matrix coupling. These groups are compared and their advantages and disadvantages highlighted. The capabilities of such models have been extensively improved by new developments in experimental methods and measurement devices. Newer technologies such as 3D printing, computed tomography scanning, X-rays, nuclear magnetic resonance, novel geophysical techniques, and use of nanomaterials allow for greater flexibilities in conducting experiments. In order for LSMs to be representative of karst aquifers, a few requirements are introduced: (1) the ability to simulate heterogeneous distributions of karst hydraulic parameters, (2) establish Darcian and non-Darcian flow regimes and exchange between the matrix and conduits, (3) placement of adequate sampling points and intervals, and (4) achieving some degree of geometric, kinematic, and dynamic similitude to represent field conditions.  相似文献   

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

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