Joint simulation of transmissivity and storativity fields conditional to steady-state and transient hydraulic head data

The self-calibrated method has been extended for the generation of equally likely realizations of transmissivity and storativity conditional to transmissivity and storativity data and to steady-state and transient hydraulic head data. Conditioning to transmissivity and storativity data is achieved by means of standard geostatistical co-simulation algorithms, whereas conditioning to hydraulic head data, given its non-linear relation to transmissivity and storativity, is achieved through non-linear optimization, similar to standard inverse algorithms. The algorithm is demonstrated in a synthetic study based on data from the WIPP site in New Mexico. Seven alternative scenarios are investigated, generating 100 realizations for each of them. The differences among the scenarios range from the number of conditioning data, to their spatial configuration, to the pumping strategies at the pumping wells. In all scenarios, the self-calibrated algorithm is able to generate transmissivity–storativity realization couples conditional to all the sample data. For the specific case studied here the results are not surprising. Of the piezometric head data, the steady-state values are the most consequential for transmissivity characterization. Conditioning to transient head data only introduces local adjustments on the transmissivity fields and serves to improve the characterization of the storativity fields.     

Geostatistical inversing for large-contrast transmissivity fields

The estimation of field parameters, such as transmissivity, is an important part of groundwater modeling. This work deals with the quasilinear geostatistical inverse approach to the estimation of the transmissivity fields from hydraulic head measurements. The standard quasilinear approach is an iterative method consisting of successive linearizations. We examine a synthetic case to evaluate the basic methodology and some modifications and extensions. The first objective is to evaluate the performance of the quasilinear approach when applied to strongly heterogeneous (or “high-contrast”) transmissivity fields and, when needed, to propose improvements that allow the solution of such problems. For large-contrast cases, the standard quasilinear method often fails to converge. However, by introducing a derivative-free line search as a polishing step after each Gauss–Newton iteration, we have found that convergence can be practically assured. Another issue is that the quasilinear procedure, which uses linearization about the best estimate to evaluate estimation variances, may lead to inaccurate estimation of the variance of the estimated variable. Our numerical results suggest that this may not be a particularly serious problem, though it is hard to say whether this conclusion will apply to other cases. Nevertheless, since the quasilinear approach is an approximation, we propose a potentially more accurate but computer-intensive Markov Chain Monte Carlo (MCMC) procedure based on conditional realizations generated through the quasilinear approach and accepted or rejected according to the Metropolis–Hastings algorithm. Six transmissivity fields with increasing contrast were generated and one thousand conditional realizations were computed for each studied case. The MCMC procedure proposed in this work gives an overall more accurate picture than the quasilinear approach but at a considerably higher computational cost.     

Transmissivity identification through multi-directional aquifer stimulation

This paper presents a geostatistical approach to multi-directional aquifer stimulation in order to better identify the transmissivity field. Hydraulic head measurements, taken at a few locations but under a number of different steady-state flow conditions, are used to estimate the transmissivity. Well installation is generally the most costly aspect of obtaining hydraulic head measurements. Therefore, it is advantageous to obtain as many informative measurements from each sampling location as possible. This can be achieved by hydraulically stimulating the aquifer through pumping, in order to set-up a variety of flow conditions. We illustrate the method by applying it to a synthetic aquifer. The simulations provide evidence that a few sampling locations may provide enough information to estimate the transmissivity field. Furthermore, the innovation of, or new information provided by, each measurement can be examined by looking at the corresponding spline and sensitivity matrix. Estimates from multi-directional stimulation are found to be clearly superior to estimates using data taken under one flow condition. We describe the geostatistical methodology for using data from multi-directional simulations and address computational issues.     

Transmissivity identification through multi-directional aquifer stimulation

This paper presents a geostatistical approach to multi-directional aquifer stimulation in order to better identify the transmissivity field. Hydraulic head measurements, taken at a few locations but under a number of different steady-state flow conditions, are used to estimate the transmissivity. Well installation is generally the most costly aspect of obtaining hydraulic head measurements. Therefore, it is advantageous to obtain as many informative measurements from each sampling location as possible. This can be achieved by hydraulically stimulating the aquifer through pumping, in order to set-up a variety of flow conditions. We illustrate the method by applying it to a synthetic aquifer. The simulations provide evidence that a few sampling locations may provide enough information to estimate the transmissivity field. Furthermore, the innovation of, or new information provided by, each measurement can be examined by looking at the corresponding spline and sensitivity matrix. Estimates from multi-directional stimulation are found to be clearly superior to estimates using data taken under one flow condition. We describe the geostatistical methodology for using data from multi-directional simulations and address computational issues.     

Stochastically conditioned flow paths and travel time distribution in highly heterogeneous synthesized aquifers / Chemins d’écoulement et distribution des temps de parcours au sein d’aquifères synthétiques hautement hétérogènes sous conditions stochastiques

Abstract

Abstract Characterization of heterogeneity at the field scale generally requires detailed aquifer properties such as transmissivity and hydraulic head. An accurate delineation of these properties is expensive and time consuming, and for many if not most groundwater systems, is not practical. As an alternative approach, stochastic representation of random fields is used and presented in this paper. Specifically, an iterative stochastic conditional simulation approach was applied to a hypothetical and highly heterogeneous pre-designed aquifer system. The approach is similar to the classical co-kriging technique; it uses a linear estimator that depends on the covariance functions of transmissivity (T), and hydraulic head (h), as well as their cross-covariances. A linearized flow equation along with a conditional random field generator constitutes the iterative process of the conditional simulation. One hundred equally likely realizations of transmissivity fields with pre-specified geostatistical parameters were generated, and conditioned to both limited transmissivity and head data. The successful implementation of the approach resulted in conditioned flow paths and travel-time distribution under different degrees of aquifer heterogeneity. This approach worked well for fields exhibiting small variances. However, for random fields exhibiting large variances (greater than 1.0), an iterative procedure was used. The results show that, as the variance of the ln[T] increases, the flow paths tend to diverge, resulting in a wide spectrum of flow conditions, with no direct discernable relationship between the degree of heterogeneity and travel time. The applied approach indicates that high errors may result when estimation of particle travel times in a heterogeneous medium is approximated by an equivalent homogeneous medium.     

Effect of Conditioning Randomly Heterogeneous Transmissivity on Temporal Hydraulic Head Measurements in Transient Two-Dimensional Aquifer Flow

We investigated the effect of conditioning transient, two-dimensional groundwater flow simulations, where the transmissivity was a spatial random field, on time dependent head data. The random fields, representing perturbations in log transmissivity, were generated using a known covariance function and then conditioned to match head data by iteratively cokriging and solving the flow model numerically. A new approximation to the cross-covariance of log transmissivity perturbations with time dependent head data and head data at different times, that greatly increased the computational efficiency, was introduced. The most noticeable effect of head data on the estimation of head and log transmissivity perturbations occurred from conditioning only on spatially distributed head measurements during steady flow. The additional improvement in the estimation of the log transmissivity and head perturbations obtained by conditioning on time dependent head data was fairly small. On the other hand, conditioning on temporal head data had a significant effect on particle tracks and reduced the lateral spreading around the center of the paths.     

Effect of Conditioning Randomly Heterogeneous Transmissivity on Temporal Hydraulic Head Measurements in Transient Two-Dimensional Aquifer Flow

We investigated the effect of conditioning transient, two-dimensional groundwater flow simulations, where the transmissivity was a spatial random field, on time dependent head data. The random fields, representing perturbations in log transmissivity, were generated using a known covariance function and then conditioned to match head data by iteratively cokriging and solving the flow model numerically. A new approximation to the cross-covariance of log transmissivity perturbations with time dependent head data and head data at different times, that greatly increased the computational efficiency, was introduced. The most noticeable effect of head data on the estimation of head and log transmissivity perturbations occurred from conditioning only on spatially distributed head measurements during steady flow. The additional improvement in the estimation of the log transmissivity and head perturbations obtained by conditioning on time dependent head data was fairly small. On the other hand, conditioning on temporal head data had a significant effect on particle tracks and reduced the lateral spreading around the center of the paths.     

Radial basis function neural network for hydrologic inversion: an appraisal with classical and spatio-temporal geostatistical techniques in the context of site characterization

This paper investigates three techniques for spatial mapping and the consequential hydrologic inversion, using hydraulic conductivity (or transmissivity) and hydraulic head as the geophysical parameters of concern. The data for the study were obtained from the Waste Isolation and Pilot Plant (WIPP) site and surrounding area in the remote Chihuahuan Desert of southeastern New Mexico. The central technique was the Radial Basis Function algorithm for an Artificial Neural Network (RBF-ANN). An appraisal of its performance in light of classical and temporal geostatistical techniques is presented. Our classical geostatistical technique of concern was Ordinary Kriging (OK), while the method of Bayesian Maximum Entropy (BME) constituted an advanced, spatio-temporal mapping technique. A fusion technique for soft or inter-dependent data was developed in this study for use with the neural network. It was observed that the RBF-ANN is capable of hydrologic inversion for transmissivity estimation with features remaining essentially similar to that obtained from kriging. The BME technique, on the other hand, was found to reveal an ability to map localized lows and highs that were otherwise not as apparent in OK or RBF-ANN techniques.     

On the geostatistical approach to the inverse problem

The geostatistical approach to the inverse problem is discussed with emphasis on the importance of structural analysis. Although the geostatistical approach is occasionally misconstrued as mere cokriging, in fact it consists of two steps: estimation of statistical parameters (“structural analysis”) followed by estimation of the distributed parameter conditional on the observations (“cokriging” or “weighted least squares”). It is argued that in inverse problems, which are algebraically undetermined, the challenge is not so much to reproduce the data as to select an algorithm with the prospect of giving good estimates where there are no observations. The essence of the geostatistical approach is that instead of adjusting a grid-dependent and potentially large number of block conductivities (or other distributed parameters), a small number of structural parameters are fitted to the data. Once this fitting is accomplished, the estimation of block conductivities ensues in a predetermined fashion without fitting of additional parameters. Also, the methodology is compared with a straightforward maximum a posteriori probability estimation method. It is shown that the fundamental differences between the two approaches are: (a) they use different principles to separate the estimation of covariance parameters from the estimation of the spatial variable; (b) the method for covariance parameter estimation in the geostatistical approach produces statistically unbiased estimates of the parameters that are not strongly dependent on the discretization, while the other method is biased and its bias becomes worse by refining the discretization into zones with different conductivity.     

Combining geostatistics and flow simulators to identify transmissivity

The reconstruction of the transmissivity field from the more numerous experimental hydraulic head data, an inverse problem, remains the focus of continuing stochastic-based research. The difficulty of this problem arises not only from the complexity of the diffusion equation that links the two variables, but also from taking into account the physical aspects of the site under study; e.g. the boundary conditions, the effective recharge, and the geology. In practical applications, the validity of purely analytical techniques proposed to date is limited by certain simplifying assumptions, like the linearization of the flow equation, made in order to obtain a solution. For this reason, a hybrid methodology combining geostatistical techniques with deterministic numerical flow simulators is proposed. This combination allows the numerical calculation of the direct and cross covariances needed to cokrige the transmissivity from both the transmissivity and hydraulic head data. The flexibility of numerical flow simulators takes away the need for the simplifying assumptions of analytical techniques to apply the proposed methodology.     

Stochastic delineation of capture zones: classical versus Bayesian approach

A Bayesian approach to characterize the predictive uncertainty in the delineation of time-related well capture zones in heterogeneous formations is presented and compared with the classical or non-Bayesian approach. The transmissivity field is modelled as a random space function and conditioned on distributed measurements of the transmissivity. In conventional geostatistical methods the mean value of the log transmissivity and the functional form of the covariance and its parameters are estimated from the available measurements, and then entered into the prediction equations as if they are the true values. However, this classical approach accounts only for the uncertainty that stems from the lack of ability to exactly predict the transmissivity at unmeasured locations. In reality, the number of measurements used to infer the statistical properties of the transmissvity field is often limited, which introduces error in the estimation of the structural parameters. The method presented accounts for the uncertainty that originates from the imperfect knowledge of the parameters by treating them as random variables. In particular, we use Bayesian methods of inference so as to make proper allowance for the uncertainty associated with estimating the unknown values of the parameters. The classical and Bayesian approach to stochastic capture zone delineation are detailed and applied to a hypothetical flow field. Two different sampling densities on a regular grid are considered to evaluate the effect of data density in both methods. Results indicate that the predictions of the Bayesian approach are more conservative.     

Image-guided inversion in steady-state hydraulic tomography

In steady-state hydraulic tomography, the head data recorded during a series of pumping or/and injection tests can be inverted to determine the transmissivity distributions of an aquifer. This inverse problem is usually under-determined and ill-posed. We propose to use structural information inferred from a guiding image to constrain the inversion process. The guiding image can be drawn from soft data sets such as seismic and ground penetrating radar sections or from geological cross-sections inferred from the wells and some geological expertise. The structural information is extracted from the guiding image through some digital image analysis techniques. Then, it is introduced into the inversion process of the head data as a weighted four direction smoothing matrix used in the regularizer. Such smoothing matrix allows applying the smoothing along the structural features. This helps preserving eventual drops in the hydraulic properties. In addition, we apply a procedure called image-guided interpolation. This technique starts with the tomogram obtained from the image-guided inversion and focus this tomogram. These new approaches are applied on four synthetic toy problems. The hydraulic distributions estimated from the image-guided inversion are closer to the true transmissivity model and have higher resolution than those computed from a classical Gauss–Newton method with uniform isotropic smoothing.     

Worth of head data in well-capture zone design: deterministic and stochastic analysis

In this study, we examine the effects of conditioning spatially variable transmissivity fields using head and/or transmissivity measurements on well-capture zones. In order to address the challenge posed by conditioning a flow model with spatially varying parameters, an innovative inverse algorithm, the Representers method, is employed. The method explicitly considers this spatial variability.

A number of uniform measurement grids with different densities are used to condition transmissivity fields using the Representers method. Deterministic and stochastic analysis of well-capture zones are then examined. The deterministic study focuses on comparison between reference well-capture zones and their estimated mean conditioned on head data. It shows that model performance due to head conditioning on well-capture zone estimation is related to pumping rate. At moderate pumping rates transmissivity observations are more crucial to identify effects arising from small-scale variations in pore water velocity. However, with more aggressive pumping these effects are reduced, consequently model performance, through incorporating head observations, markedly improves. In the stochastic study, the effect of conditioning using head and/or transmissivity data on well-capture zone uncertainty is examined. The Representers method is coupled with the Monte Carlo method to propagate uncertainty in transmissivity fields to well-capture zones. For the scenario studied, the results showed that a combination of 48 head and transmissivity data could reduce the area of uncertainty (95% confidence interval) in well-capture zone location by over 50%, compared to a 40% reduction using either head or transmissivity data. This performance was comparable to that obtained through calibrating on three and a half times the number of head observations alone.     

利用GCV优化的正则化方法对硬X射线望远镜图象重建研究

图象重建是一个求逆过程具有不适定性。本文介绍的是一种从少数投影数据进行图象重建的线性代数法－Ｐｈｉｌｌｉｐｓ－Ｔｉｋｈｏｎｏｖ正则化方法,同时选择优化工具ＧＣＶ确定最优的正则参数。为验证该方法,我们对日本Ｙｏｈｋｏｈ卫星上塔载的硬Ｘ射线望远镜（ＨＸＴ）拍摄的太阳图象进行重建。     

How to Find Aquifer Statistics Utilizing Pumping Tests? Two Field Studies Using welltestpy

We present a workflow to estimate geostatistical aquifer parameters from pumping test data using the Python package welltestpy . The procedure of pumping test analysis is exemplified for two data sets from the Horkheimer Insel site and from the Lauswiesen site, Germany. The analysis is based on a semi-analytical drawdown solution from the upscaling approach Radial Coarse Graining, which enables to infer log-transmissivity variance and horizontal correlation length, beside mean transmissivity, and storativity, from pumping test data. We estimate these parameters of aquifer heterogeneity from type-curve analysis and determine their sensitivity. This procedure, implemented in welltestpy , is a template for analyzing any pumping test. It goes beyond the possibilities of standard methods, for example, based on Theis' equation, which are limited to mean transmissivity and storativity. A sensitivity study showed the impact of observation well positions on the parameter estimation quality. The insights of this study help to optimize future test setups for geostatistical aquifer analysis and provides guidance for investigating pumping tests with regard to aquifer statistics using the open-source software package welltestpy .     

Incorporating subjective and stochastic uncertainty in an interactive multi-objective groundwater calibration framework

The interactive multi-objective genetic algorithm (IMOGA) combines traditional optimization with an interactive framework that considers the subjective knowledge of hydro-geological experts in addition to quantitative calibration measures such as calibration errors and regularization to solve the groundwater inverse problem. The IMOGA is inherently a deterministic framework and identifies multiple large-scale parameter fields (typically head and transmissivity data are used to identify transmissivity fields). These large-scale parameter fields represent the optimal trade-offs between the different criteria (quantitative and qualitative) used in the IMOGA. This paper further extends the IMOGA to incorporate uncertainty both in the large-scale trends as well as the small-scale variability (which can not be resolved using the field data) in the parameter fields. The different parameter fields identified by the IMOGA represent the uncertainty in large-scale trends, and this uncertainty is modeled using a Bayesian approach where calibration error, regularization, and the expert’s subjective preference are combined to compute a likelihood metric for each parameter field. Small-scale (stochastic) variability is modeled using a geostatistical approach and added onto the large-scale trends identified by the IMOGA. This approach is applied to the Waste Isolation Pilot Plant (WIPP) case-study. Results, with and without expert interaction, are analyzed and the impact that expert judgment has on predictive uncertainty at the WIPP site is discussed. It is shown that for this case, expert interaction leads to more conservative solutions as the expert compensates for some of the lack of data and modeling approximations introduced in the formulation of the problem.     

A comparison of two stochastic inverse methods in a field-scale application

Inverse modeling is a useful tool in ground water flow modeling studies. The most frequent difficulties encountered when using this technique are the lack of conditioning information (e.g., heads and transmissivities), the uncertainty in available data, and the nonuniqueness of the solution. These problems can be addressed and quantified through a stochastic Monte Carlo approach. The aim of this work was to compare the applicability of two stochastic inverse modeling approaches in a field-scale application. The multi-scaling (MS) approach uses a downscaling parameterization procedure that is not based on geostatistics. The pilot point (PP) approach uses geostatistical random fields as initial transmissivity values and an experimental variogram to condition the calibration. The studied area (375 km2) is part of a regional aquifer, northwest of Montreal in the St. Lawrence lowlands (southern Québec). It is located in limestone, dolomite, and sandstone formations, and is mostly a fractured porous medium. The MS approach generated small errors on heads, but the calibrated transmissivity fields did not reproduce the variogram of observed transmissivities. The PP approach generated larger errors on heads but better reproduced the spatial structure of observed transmissivities. The PP approach was also less sensitive to uncertainty in head measurements. If reliable heads are available but no transmissivities are measured, the MS approach provides useful results. If reliable transmissivities with a well inferred spatial structure are available, then the PP approach is a better alternative. This approach however must be used with caution if measured transmissivities are not reliable.