首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
A calibration method to solve the groundwater inverse problem under steady- and transient-state conditions is presented. The method compares kriged and numerical head field gradients to modify hydraulic conductivity without the use of non-linear optimization techniques. The process is repeated iteratively until a close match with piezometric data is reached. The approach includes a damping factor to avoid divergence and oscillation of the solution in areas of low hydraulic gradient and a weighting factor to account for temporal head variation in transient simulations. The efficiency of the method in terms of computing time and calibration results is demonstrated with a synthetic field. It is shown that the proposed method provides parameter fields that reproduce both hydraulic conductivity and piezometric data in few forward model solutions. Stochastic numerical experiments are conducted to evaluate the sensitivity of the method to the damping function and to the head field estimation errors.  相似文献   

2.
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.  相似文献   

3.
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.  相似文献   

4.
A physically based inverse method is developed using hybrid formulation and coordinate transform to simultaneously estimate hydraulic conductivity tensors, steady‐state flow field, and boundary conditions for a confined aquifer under ambient flow or pumping condition. Unlike existing indirect inversion techniques, the physically based method does not require forward simulations to assess model‐data misfits. It imposes continuity of hydraulic head and Darcy fluxes in the model domain while incorporating observations (hydraulic heads, Darcy fluxes, or well rates) at measurement locations. Given sufficient measurements, it yields a well‐posed inverse system of equations that can be solved efficiently with coarse grids and nonlinear optimization. When pumping and injection are active, well rates are used as measurements and flux sampling is not needed. The method is successfully tested on synthetic aquifer problems with regular and irregular geometries, different hydrofacies and flow patterns, and increasing conductivity anisotropy ratios. All problems yield stable inverse solutions under increasing head measurement errors. For a given set of observations, inversion accuracy is strongly affected by the conductivity anisotropy ratio. Conductivity estimation is also affected by flow pattern: within a hydrofacies, when Darcy flux component is very small, the corresponding directional conductivity perpendicular to streamlines becomes less identifiable. Finally, inversion is successful even if the location of aquifer boundaries is unknown. In this case, the inversion domain is defined by the location of the measurements.  相似文献   

5.
Including geophysical data in ground water model inverse calibration   总被引：1，自引：0，他引：1  
Dam D  Christensen S 《Ground water》2003,41(2):178-189
A nonlinear regression method is developed that can be used to estimate parameters of a ground waterflow model from a combination of observations of hydrological variables and observations of geophysical properties that are functionally related with the hydraulic conductivity. The procedure estimates: parameters characterizing the hydraulic conductivity field (e.g., zonal or pilot point values); geophysical properties that have been observed and that are functionally related with the hydraulic conductivity parameters; and a few parameters of the function that relates the hydraulic conductivity parameters with the geophysical properties (the type of function is assumed known). A fidelity factor, sigma(r)2, of a term of the minimized objective function reflects the faith one has in the validity of this functional relationship. The estimation methodology has been tested by means of synthetic models. The experimental results demonstrate that the number of estimated hydraulic conductivity parameters can be increased by adding geophysical observations to the set of hydrological observations that are traditionally used for model calibration. The improvement of the estimated hydraulic conductivity field and the simulated hydraulic head field can be significant but is dependent on the number, the locations, and the uncertainty of geophysical observations. The sensitivity of the estimation results to the value of sigma(r) is small for the studied problems except when the uncertainty of geophysical observations is high. In the latter case, a large sigma(r) value was found to be optimal to avoid that hydraulic conductivity estimates are closely tied to corresponding but highly uncertain geophysical observations.  相似文献   

6.
Ye Zhang 《Ground water》2014,52(3):343-351
Modeling and calibration of natural aquifers with multiple scales of heterogeneity is a challenging task due to limited subsurface access. While computer modeling plays an essential role in aquifer studies, large uncertainty exists in developing a conceptual model of an aquifer and in calibrating the model for decision making. Due to uncertainties such as a lack of understanding of subsurface processes and a lack of techniques to parameterize the subsurface environment (including hydraulic conductivity, source/sink rate, and aquifer boundary conditions), existing aquifer models often suffer nonuniqueness in calibration, leading to poor predictive capability. A robust calibration methodology is needed that can address the simultaneous estimations of aquifer parameters, source/sink, and boundary conditions. In this paper, we propose a multistage and multiscale approach that addresses subsurface heterogeneity at multiple scales, while reducing uncertainty in estimating the model parameters and model boundary conditions. The key to this approach lies in the appropriate development, verification, and synthesis of existing and new techniques of static and dynamic data integration. In particular, based on a given set of observation data, new inversion techniques can be first used to estimate aquifer large‐scale effective parameters and smoothed boundary conditions, based on which parameter and boundary condition estimation can be refined at increasing detail using standard or highly parameterized estimation techniques.  相似文献   

7.
The groundwater inverse problem of estimating heterogeneous groundwater model parameters (hydraulic conductivity in this case) given measurements of aquifer response (such as hydraulic heads) is known to be an ill-posed problem, with multiple parameter values giving similar fits to the aquifer response measurements. This problem is further exacerbated due to the lack of extensive data, typical of most real-world problems. In such cases, it is desirable to incorporate expert knowledge in the estimation process to generate more reasonable estimates. This work presents a novel interactive framework, called the ‘Interactive Multi-Objective Genetic Algorithm’ (IMOGA), to solve the groundwater inverse problem considering different sources of quantitative data as well as qualitative expert knowledge about the site. The IMOGA is unique in that it looks at groundwater model calibration as a multi-objective problem consisting of quantitative objectives – calibration error and regularization – and a ‘qualitative’ objective based on the preference of the geological expert for different spatial characteristics of the conductivity field. All these objectives are then included within a multi-objective genetic algorithm to find multiple solutions that represent the best combination of all quantitative and qualitative objectives. A hypothetical aquifer case-study (based on the test case presented by Freyberg [Freyberg DL. An exercise in ground-water model calibration and prediction. Ground Water 1988;26(3)], for which the ‘true’ parameter values are known, is used as a test case to demonstrate the applicability of this method. It is shown that using automated calibration techniques without using expert interaction leads to parameter values that are not consistent with site-knowledge. Adding expert interaction is shown to not only improve the plausibility of the estimated conductivity fields but also the predictive accuracy of the calibrated model.  相似文献   

8.
We present a geostatistically based inverse model for characterizing heterogeneity in parameters of unsaturated hydraulic conductivity for three-dimensional flow. Pressure and moisture content are related to perturbations in hydraulic parameters through cross-covariances, which are calculated to first-order. Sensitivities needed for covariance calculations are derived using the adjoint state sensitivity method. Approximations of the conditional mean parameter fields are then obtained from the cokriging estimator. Correlation between parameters and pressure – moisture content perturbations is seen to be strongly dependent on mean pressure or moisture content. High correlation between parameters and pressure data was obtained under saturated or near saturated flow conditions, providing accurate estimation of saturated hydraulic conductivity, while moisture content measurements provided accurate estimation of the pore size distribution parameter under unsaturated flow conditions.  相似文献   

9.
We present a geostatistically based inverse model for characterizing heterogeneity in parameters of unsaturated hydraulic conductivity for three-dimensional flow. Pressure and moisture content are related to perturbations in hydraulic parameters through cross-covariances, which are calculated to first-order. Sensitivities needed for covariance calculations are derived using the adjoint state sensitivity method. Approximations of the conditional mean parameter fields are then obtained from the cokriging estimator. Correlation between parameters and pressure – moisture content perturbations is seen to be strongly dependent on mean pressure or moisture content. High correlation between parameters and pressure data was obtained under saturated or near saturated flow conditions, providing accurate estimation of saturated hydraulic conductivity, while moisture content measurements provided accurate estimation of the pore size distribution parameter under unsaturated flow conditions.  相似文献   

10.
The hydraulic gradient between aquifers and rivers is one of the most variable properties in a river/aquifer system. Detailed process understanding of bank storage under hydraulic gradients is obtained from a two‐dimensional numerical model of a variably saturated aquifer slice perpendicular to a river. Exchange between the river and the aquifer occurs first at the interface with the unsaturated zone. The proportion of total water exchanged through the river bank compared to the river bed is a function of aquifer hydraulic conductivity, partial penetration, and hydraulic gradient. Total exchange may be estimated to within 50% using existing analytical solutions provided that unsaturated zone processes do not strongly influence exchange. Model‐calculated bank storage is at a maximum when no hydraulic gradient is present and increases as the hydraulic conductivity increases. However, in the presence of a hydraulic gradient, the largest exchange flux or distance of penetration does not necessarily correspond to the highest hydraulic conductivity, as high hydraulic conductivity increases the components of exchange both into and out of an aquifer. Flood wave characteristics do not influence ambient groundwater discharge, and so in large floods, hydraulic gradients must be high to reduce the volume of bank storage. Practical measurement of bank storage metrics is problematic due to the limitations of available measurement technologies and the nested processes of exchange that occur at the river‐aquifer interface. Proxies, such as time series concentration data in rivers and groundwater, require further development to be representative and quantitative.  相似文献   

11.
This paper presents an in situ falling-head method for measuring hydraulic conductivity of beach sediments in tidal environment. A polyvinyl chloride (PVC) standpipe was vertically pushed into the submerged beach sediments so that its lower part was filled by a sediment column. During the experiment, the sediments were submerged by sea water and the standpipe top was higher than the sea level. The pipe was fully filled with sea water at the beginning of the experiment. Then the water level time series inside and outside the standpipe were recorded. Analytical solutions were derived to describe the relation among the sediment's hydraulic conductivity and the water levels inside and outside the standpipe and used to analyze the experiment data obtained from the intertidal zone of Puqian Bay, Haikou, Hainan Province, China. The water levels predicted by the analytical solution agreed very well with all the experiment data. Experiments for horizontal hydraulic conductivity estimation were also conducted using L-shaped standpipe which bends from vertical to horizontal in the beach sediments. The averaged hydraulic conductivity anisotropy ratio at the study area is about 2.9. After each in situ experiment, the sediments in the standpipe were stored in a plastic box and transported to university laboratory to measure the hydraulic conductivity using falling-head method. It is found that the in situ hydraulic conductivity averages one order of magnitude greater than the laboratory one, indicating that the original beach surface sediments were loose due to tidal and wave actions and that the samples were significantly compacted during the transportation to laboratory.  相似文献   

12.
AN EXERCISE IN GROUND-WATER MODEL CALIBRATION AND PREDICTION   总被引：1，自引：0，他引：1  
Abstract. For a classroom exercise, nine groups of graduate students calibrated a numerical ground-water flow model to a set of perfectly observed hydraulic head data for a hypothetical phreatic aquifer. All groups used exactly the same numerical model and identical sets of observed data. After calibration, the students predicted the hydraulic head distribution in the aquifer resulting from a modification in one boundary condition. A quantitative analysis of the results of this calibration-prediction exercise vividly demonstrates some of the difficulties in parameter identification for ground-water flow models. Group predictions differed significantly. Successful prediction was strongly correlated with successful estimation of conductivity values, and was essentially unrelated to successful estimation of aquifer bottom elevations or with the number of trial-and-error simulations required for calibration. Most importantly, success in prediction was unrelated to success in matching observed heads under premodification conditions. In this sense, good calibration did not lead to good prediction.  相似文献   

13.
This paper presents the Kalman Filtered Double Constraint Method (DCM‐KF) as a technique to estimate the hydraulic conductivities in the grid blocks of a groundwater flow model. The DCM is based on two forward runs with the same initial grid block conductivities, but with alternating flux‐head conditions specified on parts of the boundary and the wells. These two runs are defined as: (1) the flux run, with specified fluxes (recharge and well abstractions), and (2) the head run, with specified heads (measured in piezometers). Conductivities are then estimated as the initial conductivities multiplied by the fluxes obtained from the flux run and divided by the fluxes obtained from the head run. The DCM is easy to implement in combination with existing models (e.g., MODFLOW). Sufficiently accurate conductivities are obtained after a few iterations. Because of errors in the specified head‐flux couples, repeated estimation under varying hydrological conditions results in different conductivities. A time‐independent estimate of the conductivities and their inaccuracy can be obtained by a simple linear KF with modest computational requirements. For the Kleine Nete catchment, Belgium, the DCM‐KF yields sufficiently accurate calibrated conductivities. The method also results in distinguishing regions where the head‐flux observations influence the calibration from areas where it is not able to influence the hydraulic conductivity.  相似文献   

14.
A groundwater model characterized by a lack of field data about hydraulic model parameters and boundary conditions combined with many observation data sets for calibration purpose was investigated concerning model uncertainty. Seven different conceptual models with a stepwise increase from 0 to 30 adjustable parameters were calibrated using PEST. Residuals, sensitivities, the Akaike information criterion (AIC and AICc), Bayesian information criterion (BIC), and Kashyap's information criterion (KIC) were calculated for a set of seven inverse calibrated models with increasing complexity. Finally, the likelihood of each model was computed. Comparing only residuals of the different conceptual models leads to an overparameterization and certainty loss in the conceptual model approach. The model employing only uncalibrated hydraulic parameters, estimated from sedimentological information, obtained the worst AIC, BIC, and KIC values. Using only sedimentological data to derive hydraulic parameters introduces a systematic error into the simulation results and cannot be recommended for generating a valuable model. For numerical investigations with high numbers of calibration data the BIC and KIC select as optimal a simpler model than the AIC. The model with 15 adjusted parameters was evaluated by AIC as the best option and obtained a likelihood of 98%. The AIC disregards the potential model structure error and the selection of the KIC is, therefore, more appropriate. Sensitivities to piezometric heads were highest for the model with only five adjustable parameters and sensitivity coefficients were directly influenced by the changes in extracted groundwater volumes.  相似文献   

15.
This paper presents the analytical properties of the sensitivity of the two-dimensional, steady-state groundwater flow equation to the flow parameters and to the boundary conditions, based on the perturbation approach. These analytical properties are used to provide guidelines for model design, model calibration and monitoring network design. The sensitivity patterns are shown to depend on the nature of both the perturbed parameter and the variable investigated. Indeed, the sensitivity of the hydraulic head to the hydraulic conductivity extends mainly in the flow direction, while the sensitivity to the recharge spreads radially. Besides, the sensitivity of the flow longitudinal velocity to the hydraulic conductivity propagates in both the longitudinal and transverse directions, whereas the sensitivity of the flow transverse velocity propagates in the diagonal directions to the flow. The analytical results are confirmed by application examples on idealized and real-world simulations. These analytical findings allow some general rules to be established for model design, model calibration and monitoring network design. In particular, the optimal location of measurement points depends on the nature of the variable of interest. Measurement network design thus proves to be problem-dependent. Moreover, adequate monitoring well network design may allow to discriminate between the possible sources of error.  相似文献   

16.
Sedimentary units generally present anisotropy in their hydraulic properties, with higher hydraulic conductivity along bedding planes, rather than perpendicular to them. This common property leads to a modeling challenge if the sedimentary structure is folded. In this paper, we show that the gradient of the geological potential used by implicit geological modeling techniques can be used to compute full hydraulic conductivity tensors varying in space according to the geological orientation. For that purpose, the gradient of the potential, a vector normal to the bedding, is used to construct a rotation matrix that allows the estimation of the 3D hydraulic conductivity tensor in a single matrix operation. A synthetic 2D cross section example is used to illustrate the method and show that flow simulations performed in such a folded environment are highly influenced by this rotating anisotropy. When using the proposed method, the streamlines follow very closely the folded formation. This is not the case with an isotropic model.  相似文献   

17.
Hydraulic tomography (HT) is a method for resolving the spatial distribution of hydraulic parameters to some extent, but many details important for solute transport usually remain unresolved. We present a methodology to improve solute transport predictions by combining data from HT with the breakthrough curve (BTC) of a single forced‐gradient tracer test. We estimated the three dimensional (3D) hydraulic‐conductivity field in an alluvial aquifer by inverting tomographic pumping tests performed at the Hydrogeological Research Site Lauswiesen close to Tübingen, Germany, using a regularized pilot‐point method. We compared the estimated parameter field to available profiles of hydraulic‐conductivity variations from direct‐push injection logging (DPIL), and validated the hydraulic‐conductivity field with hydraulic‐head measurements of tests not used in the inversion. After validation, spatially uniform parameters for dual‐domain transport were estimated by fitting tracer data collected during a forced‐gradient tracer test. The dual‐domain assumption was used to parameterize effects of the unresolved heterogeneity of the aquifer and deemed necessary to fit the shape of the BTC using reasonable parameter values. The estimated hydraulic‐conductivity field and transport parameters were subsequently used to successfully predict a second independent tracer test. Our work provides an efficient and practical approach to predict solute transport in heterogeneous aquifers without performing elaborate field tracer tests with a tomographic layout.  相似文献   

18.
SCE-UA方法在新安江模型参数优化中的应用   总被引：9，自引：0，他引：9  
以前在使用新安江模型时人们遇到的最大困难可归因于缺乏有效的参数全局优化的数学方法,事实上对于一个缺乏经验的人来说,模型参数的人工试错计算的过程是一个相当不容易的过程,并且耗时颇多,为此,近些年来研究者们正在探索把概念性水文模型中的专家经验与自动优化计算相结合的方法或者数学优化中的全局优化方法,如SEC-UA方法,本文首先简述新安江模型,而后采用3个大小和气候条件各不相同的流域对SCE-UA算法就在新安江模型计算的参数优化进行了研究,研究结果表明,SCE-UA算法用来进行新安江模型的参数优化所取得的效果是好的,从率定和检验的结果来看,SCE-UA算法可以使得率定的新安江模型的参数达到全局最优并且从概念上也合理。  相似文献   

19.
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.  相似文献   

20.
Over the past several decades, different groundwater modeling approaches of various complexities and data use have been developed. A recently developed approach for mapping hydraulic conductivity (K) and specific storage (Ss) heterogeneity is hydraulic tomography, the performance of which has not been compared to other more “traditional” methods that have been utilized over the past several decades. In this study, we compare seven methods of modeling heterogeneity which are (1) kriging, (2) effective parameter models, (3) transition probability/Markov Chain geostatistics models, (4) geological models, (5) stochastic inverse models conditioned to local K data, (6) hydraulic tomography, and (7) hydraulic tomography conditioned to local K data using data collected in five boreholes at a field site on the University of Waterloo (UW) campus, in Waterloo, Ontario, Canada. The performance of each heterogeneity model is first assessed during model calibration. In particular, the correspondence between simulated and observed drawdowns is assessed using the mean absolute error norm, (L1), mean square error norm (L2), and correlation coefficient (R) as well as through scatterplots. We also assess the various models on their ability to predict drawdown data not used in the calibration effort from nine pumping tests. Results reveal that hydraulic tomography is best able to reproduce these tests in terms of the smallest discrepancy and highest correlation between simulated and observed drawdowns. However, conditioning of hydraulic tomography results with permeameter K data caused a slight deterioration in accuracy of drawdown predictions which suggests that data integration may need to be conducted carefully.  相似文献   

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

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