首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 687 毫秒
1.
A methodology for identifying the geometry of different materials in highly heterogeneous porous media in discrete inverse problems (DIP) is described. It applies a geostatistical approach within the differential system method (DSM). DSM calculates conductivity values along an integration path beginning at a point with known conductivity. In aquifers with zero source terms, DSM completely describes the conductivity field through a spatially distributed parameter depending on hydraulic head gradients and integration path. A factor analysis of the structural components of this parameter (i.e. coregionalisation analysis) was carried out to identify the geometry of different materials, corresponding to distinct statistically homogeneous areas. The equivalent conductivity values for homogeneous areas were estimated.This approach was applied for a synthetic aquifer. The identification of geometry was accurate and the estimates of equivalent parameters were good, compared with reference values. The accuracy of the results depended on errors in hydraulic gradients, compared with conductivity gradients.  相似文献   

2.
This work presents a stochastic diagrammatic theory for the calculation of the effective hydraulic conductivity of heterogeneous media. The theory is based on the mean-flux series expansion of a log-normal hydraulic conductivity medium in terms of diagrammatic representations and leads to certain general results for the effective hydraulic conductivity of three-dimensional media. A selective summation technique is used to improve low-order perturbation analysis by evaluating an infinite set of diagrammatic terms with a specific topological structure that dominates the perturbation series. For stochastically isotropic media the selective summation yeilds the anticipated exponential expression for the effective hydraulic conductivity. This expression is extended to stochastically anisotropic media. It is also shown that in the case of non homogeneous media the uniform effective hydraulic conductivity is replaced by a non-local tensor kernel, for which general diagrammatic expressions are obtained. The non-local kernel leads to the standard exponential behavior for the effective hydraulic conductivity at the homogeneous limit.  相似文献   

3.
This work presents a stochastic diagrammatic theory for the calculation of the effective hydraulic conductivity of heterogeneous media. The theory is based on the mean-flux series expansion of a log-normal hydraulic conductivity medium in terms of diagrammatic representations and leads to certain general results for the effective hydraulic conductivity of three-dimensional media. A selective summation technique is used to improve low-order perturbation analysis by evaluating an infinite set of diagrammatic terms with a specific topological structure that dominates the perturbation series. For stochastically isotropic media the selective summation yeilds the anticipated exponential expression for the effective hydraulic conductivity. This expression is extended to stochastically anisotropic media. It is also shown that in the case of non homogeneous media the uniform effective hydraulic conductivity is replaced by a non-local tensor kernel, for which general diagrammatic expressions are obtained. The non-local kernel leads to the standard exponential behavior for the effective hydraulic conductivity at the homogeneous limit.  相似文献   

4.
This study investigates and quantifies the influence of physical heterogeneity in granular porous media, represented by materials with different hydraulic conductivity, on the migration of nitrate, used as an amendment to enhance bioremediation, under an electric field. Laboratory experiments were conducted in a bench‐scale test cell under a low applied direct current using glass bead and clay mixes and synthetic groundwater to represent ideal conditions. The experiments included bromide tracer tests in homogeneous settings to deduce controls on electrokinetic transport of inorganic solutes in the different materials, and comparison of nitrate migration under homogeneous and heterogeneous scenarios. The results indicate that physical heterogeneity of subsurface materials, represented by a contrast between a higher‐hydraulic conductivity and lower‐hydraulic conductivity material normal to the direction of the applied electric field exerts the following controls on nitrate migration: (1) a spatial change in nitrate migration rate due to changes in effective ionic mobility and subsequent accumulation of nitrate at the interface between these materials; and (2) a spatial change in the voltage gradient distribution across the hydraulic conductivity contrast, due to the inverse relationship with effective ionic mobility. These factors will contribute to higher mass transport of nitrate through low hydraulic conductivity zones in heterogeneous porous media, relative to homogeneous host materials. Overall electrokinetic migration of amendments such as nitrate can be increased in heterogeneous granular porous media to enhance the in situ bioremediation of organic contaminants present in low hydraulic conductivity zones.  相似文献   

5.
A data assimilation method is developed to calibrate a heterogeneous hydraulic conductivity field conditioning on transient pumping test data. The ensemble Kalman filter (EnKF) approach is used to update model parameters such as hydraulic conductivity and model variables such as hydraulic head using available data. A synthetical two-dimensional flow case is used to assess the capability of the EnKF method to calibrate a heterogeneous conductivity field by assimilating transient flow data from observation wells under different hydraulic boundary conditions. The study results indicate that the EnKF method will significantly improve the estimation of the hydraulic conductivity field by assimilating continuous hydraulic head measurements and the hydraulic boundary condition will significantly affect the simulation results. For our cases, after a few data assimilation steps, the assimilated conductivity field with four Neumann boundaries matches the real field well while the assimilated conductivity field with mixed Dirichlet and Neumann boundaries does not. We found in our cases that the ensemble size should be 300 or larger for the numerical simulation. The number and the locations of the observation wells will significantly affect the hydraulic conductivity field calibration.  相似文献   

6.
According to common understanding, the advective velocity of a conservative solute equals the average linear pore-water velocity. Yet direct monitoring indicates that the two velocities may be different in heterogeneous media. For example, at the Camp Dodge, Iowa, site the advective velocity of discrete Cl- plumes was less than one tenth of the average pore-water velocity calculated from Darcy's law using the measured hydraulic gradient, effective porosity, and hydraulic conductivity (K) from large-scale three-dimensional (3D) techniques, e.g., pumping tests. Possibly, this difference reflects the influence of different pore systems, if the K relevant to transient solute flux is influenced more by lower-K heterogeneity than a steady or quasi-steady water flux. To test this idea, tracer tests were conducted under controlled laboratory conditions. Under one-dimensional flow conditions, the advective velocity of discrete conservative solutes equaled the average pore-water velocity determined from volumetric flow rates and Darcy's law. In a larger 3D flow system, however, the same solutes migrated at approximately 65% of the average pore-water velocity. These results, coupled with direct observation of dye tracers and their velocities as they migrated through both homogeneous and heterogeneous sections of the same model, demonstrate that heterogeneity can slow the advective velocity of discrete solute plumes relative to the average pore-water velocity within heterogeneous 3D flow sytems.  相似文献   

7.
An ensemble Kalman filter (EnKF) is developed to identify a hydraulic conductivity distribution in a heterogeneous medium by assimilating solute concentration measurements of solute transport in the field with a steady‐state flow. A synthetic case with the mixed Neumann/Dirichlet boundary conditions is designed to investigate the capacity of the data assimilation methods to identify a conductivity distribution. The developed method is demonstrated in 2‐D transient solute transport with two different initial instant solute injection areas. The influences of the observation error and model error on the updated results are considered in this study. The study results indicate that the EnKF method will significantly improve the estimation of the hydraulic conductivity field by assimilating solute concentration measurements. The larger area of the initial distribution and the more observed data obtained, the better the calculation results. When the standard deviation of the observation error varies from 1% to 30% of the solute concentration measurements, the simulated results by the data assimilation method do not change much, which indicates that assimilation results are not very sensitive to the standard deviation of the observation error in this study. When the inflation factor is more than 1.0 to enlarge the model error by increasing the forecast error covariance matrix, the updated results of the hydraulic conductivity by the data assimilation method are not good at all. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

8.
Hydraulic conductivity distribution and plume initial source condition are two important factors affecting solute transport in heterogeneous media. Since hydraulic conductivity can only be measured at limited locations in a field, its spatial distribution in a complex heterogeneous medium is generally uncertain. In many groundwater contamination sites, transport initial conditions are generally unknown, as plume distributions are available only after the contaminations occurred. In this study, a data assimilation method is developed for calibrating a hydraulic conductivity field and improving solute transport prediction with unknown initial solute source condition. Ensemble Kalman filter (EnKF) is used to update the model parameter (i.e., hydraulic conductivity) and state variables (hydraulic head and solute concentration), when data are available. Two-dimensional numerical experiments are designed to assess the performance of the EnKF method on data assimilation for solute transport prediction. The study results indicate that the EnKF method can significantly improve the estimation of the hydraulic conductivity distribution and solute transport prediction by assimilating hydraulic head measurements with a known solute initial condition. When solute source is unknown, solute prediction by assimilating continuous measurements of solute concentration at a few points in the plume well captures the plume evolution downstream of the measurement points.  相似文献   

9.
The process of attempting to model ground-water systems requires a good understanding of the spatial variation of aquifer hydraulic properties. The capabilities of the more recent innovative flowmeters such as the electromagnetic and heat pulse flowmeters provide the sensitivity to measure ambient flows and pump-induced flows. These flowmeters provide the measurements of pump-induced vertical flows which are analyzed to obtain vertical variations in horizontal hydraulic conductivity, K(z). With discrete areal K-values, K(x, y), and vertical profiles of K, provided by multiwell testing, the essential elements are present to produce a three-dimensional hydraulic conductivity field. The advent of these new flow measuring devices has contributed much to the motivation behind this paper. This paper presents the results of applying deterministic and stochastic methodology to the three-dimensional interpolation of hydraulic properties, specifically, hydraulic conductivity, K. Three of the approaches applied in this paper are deterministic in nature, inverse-distance weighting, inverse-distance-squared weighting, and ordinary kriging, while the fourth is a stochastic approach based on self-affine fractals. All of the methods are applied to measured data collected from 14 wells at a site in the United States near Mobile, Alabama. The three-dimensional K-distributions generated by each of the methods are used as inputs to an advective based transport model with the resulting model output compared to a two-well tracer study run previously at the same site.  相似文献   

10.
Wang F  Bright J 《Ground water》2004,42(5):760-766
The influence on solute transport of the small-scale spatial variation of aquifer hydraulic conductivity (K) was analyzed by comparing results from fine-grid (2 m by 2 m) simulations of a synthetic heterogeneous aquifer to those from coarse-grid (8 m by 4 m) simulations of an equivalent homogeneous aquifer. Realizations of the K field of the heterogeneous aquifer were generated, using the Monte Carlo approach, from a lognormal distribution with mean log K of 2 (K in m/d) and three levels of log K variance of 0.1, 0.5, and 1.0. Numerical simulation results show that the average standard deviation of point concentrations increased from 1.21 to 5.78 when the value of log K variance was increased from 0.1 to 1.0. The average discrepancy between modeled concentrations (obtained from a coarse-grid deterministic numerical simulation) and the actual mean point concentrations (obtained from fine-grid Monte Carlo numerical simulations) increased from 0.91 to 4.23 with the increase in log K variance. The results from this study illustrate the uncertainty in predictions from contaminant transport models due to their inability to simulate the effects of heterogeneities at scales smaller than the model grid.  相似文献   

11.
A comprehensive analysis of steady flow patterns in saturated and unsaturated, possibly heterogeneous, isotropic soils is presented. It is shown that, at any point, the divergence of the unit tangent vector field to the streamlines is equal to the directional derivative along the streamlines of the orthogonal cross-sectional area of an infinitesimal stream tube divided by that area and also equal to the mean curvature of the surface of constant total head. Expressions are derived for the distribution of the flux, the water content, the velocity, the hydraulic conductivity, the total head, and the pressure head along a stream line or an infinitesimal, stream tube. Among the results is a simpler derivation, further interpretation, and extension of earlier work on calculating the hydraulic conductivity distribution from detailed measurements of the total head distribution in combination with measurements of the hydraulic conductivity at a few locations. In the last section, the jumps of various quantities at an interface are discussed.  相似文献   

12.
Analytical solutions for the water table and lateral discharge in a heterogeneous unconfined aquifer with time-dependent source and fluctuating river stage were derived and compared with those in an equivalent homogeneous aquifer. The heterogeneous aquifer considered consists of a number of sections of different hydraulic conductivity values. The source term and river stage were assumed to be time-dependent but spatially uniform. The solutions derived is useful in studying various groundwater flow problems in a horizontally heterogeneous aquifer since the spatially piecewise-constant hydraulic conductivity and temporally piecewise-constant recharge and lateral discharge can be used to quantify variations in these processes commonly observed in reality. Applying the solutions derived to an aquifer of three sections of different hydraulic conductivity values shown that (1) the aquifer heterogeneity significantly increases the spatial variation of the water table and thus its gradient but it has little effect on lateral discharge in the case of temporally and spatially uniform recharge, (2) the time-dependent but spatially uniform recharge increases the temporal variation of groundwater table over the entire aquifer but its effect on lateral discharge is limited in the zone near the river, and (3) the effect of river stage fluctuation on the water table and lateral discharge is limited in the zone near the river and the effect of the heterogeneity is to increase lateral discharge to or recharge from the river.  相似文献   

13.
A localized ensemble Kalman filter (EnKF) method is developed to assimilate transient flow data to calibrate a heterogeneous conductivity field. To update conductivity value at a point in a study domain, instead of assimilating all the measurements in the study domain, only limited measurement data in an area around the point are used for the conductivity updating in the localized EnKF method. The localized EnKF is proposed to solve the problems of the filter divergence usually existing in a data assimilation method without localization. The developed method is applied, in a synthetical two dimensional case, to calibrate a heterogeneous conductivity field by assimilating transient hydraulic head data. The simulations by the data assimilation with and without localized EnKF are compared. The study results indicate that the hydraulic conductivity field can be updated efficiently by the localized EnKF, while it cannot be by the EnKF. The covariance inflation and localization are found to solve the problem of the filter divergence efficiently. In comparison with the EnKF method without localization, the localized EnKF method needs smaller ensemble size to achieve stabilized results. The simulation results by the localized EnKF method are much more sensitive to conductivity correlation length than to the localization radius. The developed localized EnKF method provides an approach to improve EnKF method in conductivity calibration.  相似文献   

14.
Radar determination of the spatial structure of hydraulic conductivity   总被引：2，自引：0，他引：2  
Spatial variability of hydraulic conductivity exerts a predominant control on the flow of fluid through porous media. Heterogeneities influence advective pathways, hydrodynamic dispersion, and density-dependent dispersion; they are, therefore, a key concern for studies of ground water resource development, contaminant transport, and reservoir engineering. Ground-penetrating radar contributes to the remote, geophysical characterization of the macroscale variability of natural porous media. On a controlled excavation of a glacial-fluvial sand and gravel deposit in the Fanshawe Delta area (Ontario, Canada), the hydraulic conductivity field of a 45 x 3 m vertical exposure was characterized using constant-head permeameter measurements performed on undisturbed horizontal sediment cores. Ground-penetrating radar data were collected along the excavation face in the form of both reflection and common midpoint surveys. Comparison of geostatistical analyses of the permeameter measurements and the radar data suggests thatthe horizontal correlation structure of radar stack velocity can be used to directly infer the horizontal correlation structure of hydraulic conductivity. The averaging nature of the common midpoint survey is manifest in the vertical correlation structure of stack velocity, making it less useful. Radar reflection data do not exhibit a spatial structure similar to that of hydraulic conductivity possibly because reflections are a result of material property contrasts rather than the material properties themselves.  相似文献   

15.
The analytic element method is well suited for the Gardner hydraulic conductivity function, but is limited in describing real soils. Therefore, parameter equivalence between the van Genuchten and Gardner hydraulic conductivity functions is explored for the case of steady vertical flow through a homogeneous medium with a single inclusion, i.e., a binary soil. The inclusion has different hydraulic parameters than the background medium. Equivalence is established using three methods: (1) effective capillary drive; (2) capillary length; (3) and a least-squares optimization method that aims to fit a Gardner function to a corresponding van Genuchten function by minimizing the difference in log conductivity over a specified pressure range. Comparisons between hydraulic models are made based on scatterplots of pressure head and the vertical Darcian flux obtained using a finite-element numerical solution with both constitutive relations. For applicability of an equivalent Gardner function over a broad range of pressure heads, the crossover pressure must be maintained between the two parametric functions. The crossover pressure is defined as the pressure in which the hydraulic conductivity of the inclusion is equal to the background. It can be shown that a hybrid methodology of preserving the crossover pressure exactly and using the effective capillary drive will result in hydraulic parameters that are easily obtained and provide good agreement between the conductivity functions of the GR model to the VG model.  相似文献   

16.
We consider heterogeneous media whose properties vary in space and particularly aquifers whose hydraulic conductivity K may change by orders of magnitude in the same formation. Upscaling of conductivity in models of aquifer flow is needed in order to reduce the numerical burden, especially when modeling flow in heterogeneous aquifers of 3D random structure. Also, in many applications the interest is in average values of the dependent variables over scales larger or comparable to the conductivity length scales. Assigning values of the conductivity Kb to averaging domains, or computational blocks, is the topic of a large body of literature, the problem being of wide interest in various branches of physics and engineering. It is clear that upscaling causes loss of information and at best it can render a good approximation of the fine scale solution after averaging it over the blocks.The present article focuses on upscaling approaches dealing with random media. It is not meant to be a review paper, its main scope being to elucidate a few issues of principle and to briefly discuss open questions. We show that upscaling can be usually achieved only approximately, and the result may depend on the particular upscaling scheme adopted. The typically scarce information on the statistical structure of the fine-scale conductivity imposes a strong limitation to the upscaling problem. Also, local upscaling is not possible in nonuniform mean flows, for which the upscaled conductivity tensor is generally nonlocal and it depends on the domain geometry and the boundary conditions. These and other limitations are discussed, as well as other open topics deserving further investigation.  相似文献   

17.
Probabilistic analysis by Monte Carlo Simulation method (MCSM) is a computationally prohibitive task for a reactive solute transport involving coupled PDEs with nonlinear source/sink terms in 3-D heterogeneous porous media. The perturbation based stochastic finite element method (SFEM) is an attractive alternative method to MCSM as it is computationally efficient and accurate. In the present study SFEM is developed for solving nonlinear reactive solute transport problem in a 3-D heterogeneous medium. Here the solution of the biodegradation problem involving a single solute by a single class of microorganisms coupled with dynamic microbial growth is attempted using this method. The SFEM here produces a second-order accurate solution for the mean and a first-order accurate solution for the standard deviation of concentrations. In this study both the physical parameters (hydraulic conductivity, porosity, dispersivity and diffusion coefficient) and the biological parameters (maximum substrate utilization rate and the coefficient of cell decay) are considered as spatially varying random fields. A comparison between the MCSM and SFEM for the mean and standard deviation of concentration is made for 1-D and 3-D problem. The effects of heterogeneity on the degradation of substrate and growth of biomass concentrations for a range of variances of input parameters are discussed for both 1-D and 3-D problems.  相似文献   

18.
Three-dimensional numerical simulations using a detailed synthetic hydraulic conductivity field developed from geological considerations provide insight into the scaling of subsurface flow and transport processes. Flow and advective transport in the highly resolved heterogeneous field were modeled using massively parallel computers, providing a realistic baseline for evaluation of the impacts of parameter scaling. Upscaling of hydraulic conductivity was performed at a variety of scales using a flexible power law averaging technique. A series of tests were performed to determine the effects of varying the scaling exponent on a number of metrics of flow and transport behavior. Flow and transport simulation on high-performance computers and three-dimensional scientific visualization combine to form a powerful tool for gaining insight into the behavior of complex heterogeneous systems.Many quantitative groundwater models utilize upscaled hydraulic conductivity parameters, either implicitly or explicitly. These parameters are designed to reproduce the bulk flow characteristics at the grid or field scale while not requiring detailed quantification of local-scale conductivity variations. An example from applied groundwater modeling is the common practice of calibrating grid-scale model hydraulic conductivity or transmissivity parameters so as to approximate observed hydraulic head and boundary flux values. Such parameterizations, perhaps with a bulk dispersivity imposed, are then sometimes used to predict transport of reactive or non-reactive solutes. However, this work demonstrates that those parameters that lead to the best upscaling for hydraulic conductivity and head do not necessarily correspond to the best upscaling for prediction of a variety of transport behaviors. This result reflects the fact that transport is strongly impacted by the existence and connectedness of extreme-valued hydraulic conductivities, in contrast to bulk flow which depends more strongly on mean values. It provides motivation for continued research into upscaling methods for transport that directly address advection in heterogeneous porous media.An electronic version of this article is available online at the journal's homepage at http://www.elsevier.nl/locate/advwatres or http://www.elsevier.com/locate/advwatres (see “Special section on vizualization”. The online version contains additional supporting information, graphics, and a 3D animation of simulated particle movement.©1998 Elsevier Science Limited. All rights reserved  相似文献   

19.
The hydraulic conductivity of heterogeneous porous media depends on the distribution function and the geometry of local conductivities at the smaller scale. There are various approaches to estimate the effective conductivity Keff at the larger scale based on information about the small scale heterogeneity. A critical geometric property in this ‘upscaling’ procedure is the spatial connectivity of the small-scale conductivities. We present an approach based on the Euler-number to quantify the topological properties of heterogeneous conductivity fields, and we derive two key parameters which are used to estimate Keff. The required coefficients for the upscaling formula are obtained by regression based on numerical simulations of various heterogeneous fields. They are found to be generally valid for various different isotropic structures. The effective unsaturated conductivity function Keff (ψm) could be predicted satisfactorily. We compare our approach with an alternative based on percolation theory and critical path analysis which yield the same type of topological parameters. An advantage of using the Euler-number in comparison to percolation theory is the fact that it can be obtained from local measurements without the need to analyze the entire structure. We found that for the heterogeneous field used in this study both methods are equivalent.  相似文献   

20.
Non-unique solutions of inverse problems arise from a lack of information that satisfies necessary conditions for the problem to be well defined. This paper investigates these conditions for inverse modeling of water flow through multi-dimensional variably saturated porous media. It shows that in order to obtain a unique estimate of hydraulic parameters, along each streamline of the flow field (1) spatial and temporal head observations must be given; (2) the number of spatial and temporal head observations required should be greater or equal to the number of unknown parameters; (3) the flux boundary condition or the pumping rate of a well must be specified for the homogeneous case and both boundary flux and pumping rate are a must for the heterogeneous case; (4) head observations must encompass both saturated and unsaturated conditions, and the functional relationships for unsaturated hydraulic conductivity/pressure head and for the moisture retention should be given, and (5) the residual water content value also need to be specified a priori or water content measurements are needed for the estimation of the saturated water content.For field problems, these necessary conditions can be collected or estimated but likely involve uncertainty. While the problems become well defined and have unique solutions, the solutions likely will be uncertain. Because of this uncertainty, stochastic approaches are deemed to be appropriate for inverse problems as they are for forward problems to address uncertainty. Nevertheless, knowledge of these necessary conditions is critical to reduce uncertainty in both characterization of the vadose zone and the aquifer, and prediction of water flow and solute migration in the subsurface.  相似文献   

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

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