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.     

Conventional and combined pump-and-treat systems under nonuniform background flow

We investigate the performance of vertical hydraulic barriers in combination with extraction wells for the partial hydraulic isolation of contaminated aquifer areas. The potential advantage of such combinations compared to a conventional pump-and-treat system has already been demonstrated in a previous study. Here we extend the scope of the performance analysis to the impact of uncertainty in the regional flow direction as well as to highly heterogeneous aquifer transmissivity distributions. In addition, two new well-barrier scenarios are proposed and analyzed. The hydraulic efficiency of the scenarios is rated based on the expected (mean) reduction of the pumping rate that is required to achieve downgradient contaminant capture. The uncertain spatial distribution of aquifer transmissivity is considered by means of unconditioned Monte Carlo simulations. The significance of uncertain background flow conditions is incorporated by computing minimized pumping rates for deviations of the regional flow direction up to 30 degrees from a normative base case. The results give an answer on how pumping rates have to be changed for each barrier-well combination in order to achieve robust systems. It is exposed that in comparison to installing exclusively wells, the barrier-supported approach generally yields savings in the (average) pumping rate. The particular efficiency is shown to be highly dependent on the interaction of variance and integral scale of transmissivity distribution, well and barrier position, as well as direction of background flow.     

Interpretation of Pumping Tests in Heterogeneous Aquifers with Constant Head Boundary

This paper investigates the impact of heterogeneity of the transmissivity field on the interpretation of steady-state pumping test data from aquifer systems delimited by constant head boundaries such as aquifers adjacent to lakes or rivers. Spatially variable transmissivity fields are randomly generated and used to simulate the drawdown due to a pumping well located at different distances from a constant head boundary. The steady-state drawdown simulated at different observation wells are then interpreted using the Hantush method (Hantush 1959). The numerical simulations show that, in contrast to the case of infinite aquifer domains, the interpreted transmissivity varies depending on well locations and the separation distance between pumping well and boundary relative to the correlation length. The ensemble-averaged estimated transmissivity varies between the geometric mean and the arithmetic mean, and can even exceed the arithmetic mean in a narrow domain adjacent to the boundary. It approaches the geometric mean of the underlying transmissivity field only if the distance between the pumping well is more than 20 times the characteristic length of the transmissivity field.     

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.     

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.     

Optimal implicit staggered‐grid finite‐difference schemes based on the sampling approximation method for seismic modelling

We propose new implicit staggered‐grid finite‐difference schemes with optimal coefficients based on the sampling approximation method to improve the numerical solution accuracy for seismic modelling. We first derive the optimized implicit staggered‐grid finite‐difference coefficients of arbitrary even‐order accuracy for the first‐order spatial derivatives using the plane‐wave theory and the direct sampling approximation method. Then, the implicit staggered‐grid finite‐difference coefficients based on sampling approximation, which can widen the range of wavenumber with great accuracy, are used to solve the first‐order spatial derivatives. By comparing the numerical dispersion of the implicit staggered‐grid finite‐difference schemes based on sampling approximation, Taylor series expansion, and least squares, we find that the optimal implicit staggered‐grid finite‐difference scheme based on sampling approximation achieves greater precision than that based on Taylor series expansion over a wider range of wavenumbers, although it has similar accuracy to that based on least squares. Finally, we apply the implicit staggered‐grid finite difference based on sampling approximation to numerical modelling. The modelling results demonstrate that the new optimal method can efficiently suppress numerical dispersion and lead to greater accuracy compared with the implicit staggered‐grid finite difference based on Taylor series expansion. In addition, the results also indicate the computational cost of the implicit staggered‐grid finite difference based on sampling approximation is almost the same as the implicit staggered‐grid finite difference based on Taylor series expansion.     

An improved straight-line fitting method for analyzing pumping test recovery data

Theis (1935) derived an exact solution for the residual drawdown in a well after the cessation of a pumping test by summing two drawdowns: one (s1), caused by imaginary continuation of the original pumping and the other (s2), due to an imaginary injection at the same constant rate. We approximated the Theis solution to obtain a simple linear relation for determining the transmissivity and storage coefficient from recovery data. Unlike other existing straight-line fitting methods, in our method, we applied different approximations to the well functions in the solutions of s1 and s2. We used the well-known Cooper-Jacob approximation for s1, truncating the expansion of the well function in s2 to its first three terms. For the same level of truncation errors, while the Cooper-Jacob approximation requires the argument u1相似文献   

A method for evaluating horizontal well pumping tests

Predicting the future performance of horizontal wells under varying pumping conditions requires estimates of basic aquifer parameters, notably transmissivity and storativity. For vertical wells, there are well-established methods for estimating these parameters, typically based on either the recovery from induced head changes in a well or from the head response in observation wells to pumping in a test well. Comparable aquifer parameter estimation methods for horizontal wells have not been presented in the ground water literature. Formation parameter estimation methods based on measurements of pressure in horizontal wells have been presented in the petroleum industry literature, but these methods have limited applicability for ground water evaluation and are based on pressure measurements in only the horizontal well borehole, rather than in observation wells. This paper presents a simple and versatile method by which pumping test procedures developed for vertical wells can be applied to horizontal well pumping tests. The method presented here uses the principle of superposition to represent the horizontal well as a series of partially penetrating vertical wells. This concept is used to estimate a distance from an observation well at which a vertical well that has the same total pumping rate as the horizontal well will produce the same drawdown as the horizontal well. This equivalent distance may then be associated with an observation well for use in pumping test algorithms and type curves developed for vertical wells. The method is shown to produce good results for confined aquifers and unconfined aquifers in the absence of delayed yield response. For unconfined aquifers, the presence of delayed yield response increases the method error.     

Issues in characterizing heterogeneity and connectivity in non-multiGaussian media

The performances of kriging, stochastic simulations and sequential self-calibration inversion are assessed when characterizing a non-multiGaussian synthetic 2D braided channel aquifer. The comparison is based on a series of criteria such as the reproduction of the original reference transmissivity or head fields, but also in terms of accuracy of flow and transport (capture zone) forecasts when the flow conditions are modified. We observe that the errors remain large even for a dense data network. In addition some unexpected behaviours are observed when large transmissivity datasets are used. In particular, we observe an increase of the bias with the number of transmissivity data and an increasing uncertainty with the number of head data. This is interpreted as a consequence of the use of an inadequate multiGaussian stochastic model that is not able to reproduce the connectivity of the original field.     

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.     

Random errors in the hydrogeological parameters of a two-layer stratum evaluated by multi-well pumping tests

Random errors in estimates of flow and storage parameters of a two-layer stratum, derived from multi-well pumping test data are examined. It was found that the relative random error with a confidence probability of 0.954 is appreciably less than 100% only for the estimates of aquifer transmissivity, while the estimates of the hydraulic conductivity and the specific yield of the overlying deposits show relative random errors with the same confidence probability in excess of 100%. Therefore, the random errors of parameters are to be given as multiples of the established values.     

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.     

Fractional conservation of mass

The traditional conservation of mass equation is derived using a first-order Taylor series to represent flux change in a control volume, which is valid strictly for cases of linear changes in flux through the control volume. We show that using higher-order Taylor series approximations for the mass flux results in mass conservation equations that are intractable. We then show that a fractional Taylor series has the advantage of being able to exactly represent non-linear flux in a control volume with only two terms, analogous to using a first-order traditional Taylor series. We replace the integer-order Taylor series approximation for flux with the fractional-order Taylor series approximation, and remove the restriction that the flux has to be linear, or piece-wise linear, and remove the restriction that the control volume must be infinitesimal. As long as the flux can be approximated by a power-law function, the fractional-order conservation of mass equation will be exact when the fractional order of differentiation matches the flux power-law. There are two important distinctions between the traditional mass conservation, and its fractional equivalent. The first is that the divergence term in the fractional mass conservation equation is the fractional divergence, and the second is the appearance of a scaling term in the fractional conservation of mass equation that may eliminate scale effects in parameters (e.g., hydraulic conductivity) that should be scale-invariant.     

Static and dynamic conceptual model of a complexly fractured crystalline rock aquifer

A conceptual model of anisotropic and dynamic permeability is developed from hydrogeologic and hydromechanical characterization of a foliated, complexly fractured, crystalline rock aquifer at Gates Pond, Berlin, Massachusetts. Methods of investigation include aquifer‐pumping tests, long‐term hydrologic monitoring, fracture characterization, downhole heat‐pulse flow meter measurements, in situ extensometer testing, and earth tide analysis. A static conceptual model is developed from observations of depth‐dependent and anisotropic permeability that effectively compartmentalizes the aquifer as a function of foliation intensity. Superimposed on the static model is dynamic permeability as a function of hydraulic head in which transient bulk aquifer transmissivity is proportional to changes in hydraulic head due to hydromechanical coupling. The dynamic permeability concept is built on observations that fracture aperture changes as a function of hydraulic head, as measured during in situ extensometer testing of individual fractures, and observed changes in bulk aquifer transmissivity as determined from earth tides during seasonal changes in hydraulic head, with higher transmissivity during periods of high hydraulic head, and lower transmissivity during periods of relatively lower hydraulic head. A final conceptual model is presented that captures both the static and dynamic properties of the aquifer. The workflow presented here demonstrates development of a conceptual framework for building numerical models of complexly fractured, foliated, crystalline rock aquifers that includes both a static model to describe the spatial distribution of permeability as a function of fracture type and foliation intensity and a dynamic model that describes how hydromechanical coupling impacts permeability magnitude as a function of hydraulic head fluctuation. This model captures important geologic controls on permeability magnitude, anisotropy, and transience and therefor offers potentially more reliable history matching and forecasts of different water management strategies, such as resource evaluation, well placement, permeability prediction, and evaluating remediation strategies.     

Dynamic impedances of a weakly anelastic medium

Approximate formulas are derived to evaluate in the frequency domain the dynamic impedances of a weakly anelastic medium based on its pure elastic behaviour. The correspondence principle is applied to the elastic solution of a boundary-value problem followed by the expansion of the anelastic solution in a Taylor series about the elastic state. Taking the magnitude of material damping into account (small damping) only the first two terms of the Taylor series have been used. The derivatives of first order in the expansion can be determined by the central difference approximation; this requires only the evaluation of differences between neighbouring elastic solutions.     

Cautions and Suggestions for Geochemical Sampling in Fractured Rock

Collecting water samples for geochemical analyses in open bedrock boreholes or in discrete intervals of boreholes intersected by multiple fractures is likely to yield ambiguous results for ground water chemistry because of the variability in the transmissivity, storativity, and hydraulic head of fractures intersecting the borehole. Interpreting chemical analyses of water samples collected in bedrock boreholes requires an understanding of the hydraulic conditions in the borehole under the ambient flow regime in the aquifer as well as during sampling. Pumping in open boreholes, regardless of the pumping rate and the location of the pump intake, first draws water from the borehole and then from fractures intersecting the borehole. The time at which the volumetric rate of water entering the borehole from fractures is approximately equal to the pumping rate can be identified by monitoring the logarithm of drawdown in the borehole as a function of the logarithm of time. Mixing of water entering the borehole from fractures with water in the borehole must be considered in estimating the time at which the pump discharge is representative of aquifer water. In boreholes intersected by multiple fractures, after the contribution from the borehole volume has diminished, the contribution of fractures to the pump discharge will be weighted according to their transmissivity, regardless of the location of the pump intake. This results in a flux-averaged concentration in the pump discharge that is biased by the chemical signature of those fractures with the highest transmissivity. Under conditions where the hydraulic head of fractures varies over the length of the borehole, open boreholes will be subject to ambient flow in the water column in the borehole. In some instances, the magnitude of the ambient flow may be similar to the designated pumping rate for collecting water samples for geochemical analyses. Under such conditions, the contributions to the pump discharge from individual fractures will be a function not only of the transmissivity of the fractures, but also of the distribution of hydraulic head in fractures intersecting the borehole. To reduce or eliminate the deleterious effects of conducting geochemical sampling in open boreholes, a straddle-packer apparatus that isolates a single fracture or a series of closely spaced fractures is recommended. It is also recommended that open boreholes be permanently outfitted with borehole packers or borehole liners in instances where maintaining the hydraulic and chemical stratification in the aquifer is of importance. In a field example, a comparison of results from sampling in an open borehole and in discrete intervals of the same borehole showed dramatic differences in the concentrations of chemical constituents in the water samples, even though chemical field parameters stabilized prior to both open borehole and discrete interval sampling.     

基于改进模态参数灵敏度法的结构损伤识别研究

推导了模态参数对于损伤构件的一阶和二阶灵敏度矩阵,并对在推导一阶和二阶振型灵敏度的过程中产生的模态截尾误差进行了改进。根据泰勒级数展开的原理分别建立了一阶和二阶的灵敏度方程。考虑到一阶灵敏度方程求解速度快和二阶灵敏度方程求解精度高的特点,本文提出了一种用于结构损伤识别的混合迭代算法,该算法用二阶非线性的解析解作为算法的第一次迭代值,用一阶灵敏度方程的求解值对该算法的第一次迭代值进行关于泰勒级数截尾误差的修正。研究表明,本文提出的混合迭代算法由于采用了精确度较高的二阶非线性解析解作为迭代修正的初值,因此,迭代修正精度更高,收敛性更好。