SUPRPUMP: An Interactive Program for Well Test Analysis and Design

Abstract. We have developed a program which aids in the design and analysis of pumping tests and slug tests. In design mode, the program emphasizes calculation and plotting of the sensitivities of drawdown (or head) to well function parameters. In analysis mode, the program can analyze a given set of experimental data. For pumping tests, the program allows multiple observation wells and multiple variable-rate pumping wells. The program is written in a modular fashion, allowing easy addition of well functions to the currently existing library. An example based on a hypothetical pumping test illustrates the utility of sensitivity analysis for well test design.     

ANALYSIS OF AQUIFER TESTS CONDUCTED IN FRACTURED ROCK: A REVIEW OF THE PHYSICAL BACKGROUND AND THE DESIGN OF A COMPUTER PROGRAM FOR GENERATING TYPE CURVES

Abstract. Assessing the hydrogeological character of fractured bedrock between two or more wells is usually accomplished through an analysis of the results of traditional pumping tests. Because these tests are subject to field boundary conditions that may be different from those conducted in porous media, alternative testing techniques or analytical methods are required. To provide an adaptable interpretive package for analyzing interwell hydraulic tests in fractured rock, a FORTRAN program was developed to generate type curves for pumping tests and pulse interference tests. A review of the physical conditions which most influence the results of these tests is conducted to provide the background for the development and use of the program. The analytical solutions included in the program consider inner boundary conditions such as wellbore storage, infinitesimally thin skin and finite-thickness skin. In addition, a solution is included to analyze results influenced by wellbore storage and infinitesimally thin skin at both the pumping and observation wells. The solutions are given in the Laplace domain and are numerically inverted to generate data for the type curves. The curves can be plotted using any standard plotting program, or they can be plotted manually. The program is menu driven, easily adaptable to include additional solutions, and can be executed on a personal computer.     

A COMPUTERIZED TECHNIQUE FOR ESTIMATING THE HYDRAULIC CONDUCTIVITY OF AQUIFERS FROM SPECIFIC CAPACITY DATA

Abstract. Specific capacity data obtained from well construction reports can provide useful estimates of hydraulic conductivity (K). A simple computer program has been developed which can correct specific capacity data for partial penetration and well loss and, using an iterative technique, provide rapid estimates of K at hundreds of data points. The program allows easy data handling and is easily linked with existing statistical programs or contour mapping routines. The method was tested at two field sites in Wisconsin, one underlain by a sandy outwash aquifer, the other by fractured dolomite. In both areas, estimates of K from corrected specific capacity data agree reasonably well with data from pumping tests.     

A CONSTANT DISCHARGE DEVICE FOR AQUIFER TESTS ON LARGE-DIAMETER WELLS

Abstract. The discharge in pumping tests conducted on large-diameter wells declines with the fall of water level. This violates the basic requirement of pumping tests that the discharge should remain constant. A simple device to keep the discharge constant has therefore been developed. This device is also useful in keeping the discharge much lower than the rated capacity of the farmer's pump, and thus in increasing the total test period.
The discharge from the farmer's pump is released in a barrel having a V-notch at the top and an arrangement for returning a part of the discharge to the well. The return flow is gradually reduced through a control valve, and thus the net outflow is kept constant.
This device has been used for conducting pumping tests on several dug wells.     

A new method to characterize hydraulic short-circuits in defective borehole seals

A new approach has been developed to detect, characterize, and quantify hydraulic short-circuits in boreholes with faulty seals. The methodology, applicable to an aquifer-aquitard-aquifer system, involves a series of successive, constant-rate pumping tests in the lower aquifer while determining the leakage rate with a simultaneous nonreactive tracer test. During each pumping step, the tracer is injected under constant concentration and constant hydraulic head from a piezometer in the upper aquifer. If a seal defect exists, the tracer will follow the leak and will be recovered from the pumped water. The theoretical equations relate the leakage rate, the pumping rate, the concentration of the injected tracer, and the recovered concentration. Leakage rates can be determined for any pumping rate. The theory is tested using numerical analysis and a full-scale field test.     

SIMPLE ANALYTICAL SOLUTIONS FOR THE HP41CV PROGRAMMABLE CALCULATOR

Abstract. Two useful programs have been developed for the Hewlett Packard HP41CV programmable calculator. The THEIS program is designed to simulate a well pumping from a confined or unconfined aquifer. Drawdown, residual drawdown, t/t₁ and t/r² are calculated. The BOUN program is designed to solve for drawdown in a well pumping from an aquifer bounded by two parallel impermeable barriers. The programs can be used in aquifer pumping test design, pumping test analysis, and aquifer response predictions.     

A GENERAL PURPOSE MICROCOMPUTER AQUIFER TEST EVALUATION TECHNIQUE

Abstract. Although determination of aquifer characteristics from pumping test data is generally carried out using type curves or other graphical techniques, a number of computer methods have been developed recently for this purpose. Based on the principle of least squares, these methods of nonlinear regression analysis can be applied to any flow system for which analytical expressions of the drawdown distribution are known. In view of the growing general interest in the application of microcomputers in ground-water hydrology, a BASIC routine has been developed for estimating any number of aquifer parameters. The least squares solution is calculated by Marquardt's algorithm, using the singular-value decomposition of the Jacobian matrix. The robust computing method obtained can be applied to all kinds of pumping tests. Aquifer characteristics as well as their standard deviations are computed with optimal speed and accuracy. The technique is demonstrated by a simple application to steady flow in a leaky aquifer and an example is provided. Other applications are easily implemented and programs for unsteady-state aquifer tests, recovery tests and multiple aquifer tests are available.     

天津王4井抽水观测水汞影响因素的探讨

王4井因受井区地热开发的影响,于2008年12月断流停测。2010年对原观测井口装置进行改造,开始进行抽水观测水汞试验研究。在新的观测方式下,为确保观测数据稳定、可靠,本文从抽水观测的观测条件方面,对不同取样位置和抽水时间对水汞测值的影响进行了试验。试验结果表明,不同抽水时间和不同取样位置抽水观测对水汞测值影响不大,这一结论与水汞本身的化学特性相一致。认为在没有自流井的地区,只要选择灵敏的观测点使用抽水方式采样,固定观测条件,依然可以进行水汞监测。     

Flow to a well in a five-layer system with application to the Oxnard Basin

Nearly 40 years ago, Neuman (1968) developed an analytical solution for drawdown in a system of three aquifers separated by two aquitards when one of the aquifers is pumped at a constant rate. Whereas the simpler case of two aquifers separated by one aquitard has been presented by Neuman and Witherspoon (1969a), the full five-layer solution has not been previously evaluated. We do so here using numerical inversion of its Laplace-transformed version and present selected results graphically in dimensionless form. The solution demonstrates that the effect of pumping propagates across all five layers, adding emphasis to a question previously raised by Neuman and Witherspoon about the validity of leaky aquifer theories that disregard drawdowns in unpumped aquifers. A large-scale, long-term pumping test spanning three aquifers separated by two aquitards near Oxnard, California, has been conducted and analyzed by Neuman and Witherspoon (1972). They evaluated the vertical hydraulic diffusivities of the aquitards using the Neuman-Witherspoon ratio method and their specific storage values on the basis of laboratory consolidation tests. We reinterpret the Oxnard pumping test by coupling the five-layer analytical solution of Neuman (1968) with the parameter estimation code PEST (Doherty 2002) and validate our results against drawdowns from a subsequent pumping test at the site. Our parameter estimates compare favorably with those of Neuman and Witherspoon (1972).     

Practical issues in imaging hydraulic conductivity through hydraulic tomography

Hydraulic tomography has been developed as an alternative to traditional geostatistical methods to delineate heterogeneity patterns in parameters such as hydraulic conductivity (K) and specific storage (S(s)). During hydraulic tomography surveys, a large number of hydraulic head data are collected from a series of cross-hole tests in the subsurface. These head data are then used to interpret the spatial distribution of K and S(s) using inverse modeling. Here, we use the Sequential Successive Linear Estimator (SSLE) of Yeh and Liu (2000) to interpret synthetic pumping test data created through numerical simulations and real data generated in a laboratory sandbox aquifer to obtain the K tomograms. Here, we define "K tomogram" as an image of K distribution of the subsurface (or the inverse results) obtained via hydraulic tomography. We examine the influence of signal-to-noise ratio and biases on results using inverse modeling of synthetic and real cross-hole pumping test data. To accomplish this, we first show that the pumping rate, which affects the signal-to-noise ratio, and the order of data included into the SSLE algorithm both have large impacts on the quality of the K tomograms. We then examine the role of conditioning on the K tomogram and find that conditioning can improve the quality of the K tomogram, but can also impair it, if the data are of poor quality and conditioning data have a larger support volume than the numerical grid used to conduct the inversion. Overall, these results show that the quality of the K tomogram depends on the design of pumping tests, their conduct, the order in which they are included in the inverse code, and the quality as well as the support volume of additional data that are used in its computation.     

Redundant and Nonredundant Information for Model Calibration or Hydraulic Tomography

Drawdown data from independent pumping tests have widely been used to validate the estimated hydraulic parameters from inverse modeling or hydraulic tomography (HT). Yet, the independent pumping test has not been clearly defined. Therefore, the goal of this paper is to define this independent pumping test concept, based on the redundant or nonredundant information about aquifer heterogeneity embedded in the observed heads during cross-hole pumping tests. The definition of complete, moderate redundancy and high nonredundancy of information are stipulated using cross-correlation analysis of the relationship between the head and heterogeneity. Afterward, data from numerical experiments and field sequential pumping test campaigns reinforce the concept and the definition.     

Characterization of High-K Pathways by Borehole Flowmeter and Tracer Tests

This paper reviews different borehole flowmeter analysis methods and evaluates their applicability to a test site composed of fluvial deposits. Results from tracer and aquifer tests indicate that the aquifer is highly heterogeneous and that low-K skin effects exist at the wells. Borehole flowmeter tests were performed at 37 wells. An appropriate method for calculating borehole flowmeter K values was developed based on results from multiwell pumping tests, single-well pumping tests, and slug tests. The flowmeter data produced 881 K values. The trends and the magnitude of the K values are consistent with results from geologic investigations, recirculating tracer tests, and large-scale multiwell pumping tests. The field tests illustrate that high-K deposits can significantly affect ground-water flows in some heterogeneous fluvial aquifers.     

地层测试评价仪(FET)及其在中国海域的应用

地层测试评价仪(Formation Evaluation Tool)是中海油田服务股份有限公司引进澳大利亚 CROCKER公司技术而自行开发、制造的新一代地层测试评价系统。该仪器已经在中国近海渤海、南海等海域取得了成功应用,其精确地层压力测试能帮助判断多套油水系统,尤其在一些薄层和边水油藏也取得了可靠的测压资料,很好的解决了油水、气水界面确定问题;此外在目的层取样时所具有的污染流体抽排和地层流体特性实时检测功能保证获取不受泥浆污染的真实地层流体样品,对确定流体性质和分析流体成分起到了重要作用。本文主要介绍地层测试评价仪(FET)的原理和用途,以及在中国海域的成功应用实例。     

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.     

Single Well Tracer Tests in Aquifer Characterization

When the purpose of aquifer testing is to yield data for modeling aqueous mass transport, pumping tests and gradient measurement can only partially satisfy characterization requirements. Effective porosity, ground water flow velocity, and the vertical distribution of hydraulic conductivity within the aquifer are left as unknowns. Single well tracer methods, when added to the testing program, can be used to estimate these parameters. A drift, and pumpback test yields porosity and velocity, and point-dilution testing yields depth-discrete hydraulic information, A single emplacement of tracer into a test well is sufficient to conduct both tests. The tracer tests are facilitated by a simple method for injecting and evenly distributing the tracer solution into a wellbore, and by new ion-selective electrode instrumentation, specifically designed for submersible service, for monitoring the concentration of tracers such as bromide.     

A Method to Determine the Formation Constants of Leaky Aquifers, and Its Application to Pumping Test Data

Abstract. A method to calculate aquifer transmissivity, storage coefficient, and the leakage coefficient from pumping test data for a leaky aquifer is presented. The method is carried out by a computer program and is based on a minimization of the sum of squares of differences between drawdown in the observation well and the theoretical values from the Hantush and Jacob formula. No user defined starting points are necessary. Random error estimates for the parameters are given. Applications of the method are illustrated using data from pumping tests performed in leaky aquifers at the Cauca River Valley, Colombia.     

Temporal sampling and role of flux measurements for subsurface heterogeneous characterization in groundwater basins using hydraulic tomography

Accurate characterization of heterogeneity in groundwater basins is crucial to the sustainable management of groundwater resources. This study explores the temporal sampling issues and the role of flux measurements in the characterization of heterogeneity in groundwater basins using numerical experiments. The experiments involve a digital basin imitating the groundwater basin of the North China Plain (NCP), where the groundwater exploitation reduction program is ongoing. Using the experiments, we champion that the reduction program could collect groundwater level information induced by operational variations of existing pumping wells at different locations in the basin. Such a dataset could serve as a basin-scale hydraulic tomography (HT) to characterize the basin-scale heterogeneity cost-effectively. Both steady-state and transient-state inversion experiments demonstrate the advantage of HT surveys in characterizing basin-scale heterogeneity over conventional pumping tests at fixed well locations. Additionally, head data at the early, intermediate, and late time from well hydrographs should be selected for the HT analysis to maximize HT's power and save computational costs. When accurate geological zones are incorporated in prior information, flux measurements significantly improve parameter estimates based on conventional pumping tests. However, their effects are less noticeable for long-term HT surveys in such basin-scale aquifers without fissures or fractures. This basin-scale tomographic survey example serves a guide for field data collection and optimization of the analysis of future basin-scale HT.     

Optimization of groundwater abstraction from a coastal aquifer

ABSTRACT

An analytical mathematical model, based on Jacobian elliptic functions, has been used to identify feasible wellfield locations and pumping rates for large-scale abstraction from an unconfined coastal aquifer. The choice of optimum feasible wellfield strategy has been made using a simple economic model which calculates the cost of the pipelines required to transport the abstracted groundwater to a large coastal city which forms the demand centre. Results indicate that the cheapest wellfield design would be a single large wellfield. However, because of the need to maintain at least a minimum supply in the city until a new surface water source is developed, a better solution may well be to develop two smaller wellfields pumping a greater total abstraction.     

An Application of Hydraulic Tomography to a Large‐Scale Fractured Granite Site,Mizunami, Japan

While hydraulic tomography (HT) is a mature aquifer characterization technology, its applications to characterize hydrogeology of kilometer‐scale fault and fracture zones are rare. This paper sequentially analyzes datasets from two new pumping tests as well as those from two previous pumping tests analyzed by Illman et al. (2009) at a fractured granite site in Mizunami, Japan. Results of this analysis show that datasets from two previous pumping tests at one side of a fault zone as used in the previous study led to inaccurate mapping of fracture and fault zones. Inclusion of the datasets from the two new pumping tests (one of which was conducted on the other side of the fault) yields locations of the fault zone consistent with those based on geological mapping. The new datasets also produce a detailed image of the irregular fault zone, which is not available from geological investigation alone and the previous study. As a result, we conclude that if prior knowledge about geological structures at a field site is considered during the design of HT surveys, valuable non‐redundant datasets about the fracture and fault zones can be collected. Only with these non‐redundant data sets, can HT then be a viable and robust tool for delineating fracture and fault distributions over kilometer scales, even when only a limited number of boreholes are available. In essence, this paper proves that HT is a new tool for geologists, geophysicists, and engineers for mapping large‐scale fracture and fault zone distributions.     

A review of thermal response test analysis using pumping test concepts

The design of ground-coupled heat pump systems requires knowledge of the thermal properties of the subsurface and boreholes. These properties can be measured with in situ thermal response tests (TRT), where a heat transfer fluid flowing in a ground heat exchanger is heated with an electric element and the resulting temperature perturbation is monitored. These tests are analogous to standard pumping tests conducted in hydrogeology, because a system that is initially assumed at equilibrium is perturbed and the response is monitored in time, to assess the system's properties with inverse modeling. Although pumping test analysis is a mature topic in hydrogeology, the current analysis of temperature measurements in the context of TRTs is comparatively a new topic and it could benefit from the application of concepts related to pumping tests. The purpose of this work is to review the methodology of TRTs and improve their analysis using pumping test concepts, such as the well function, the superposition principle, and the radius of influence. The improvements are demonstrated with three TRTs. The first test was conducted in unsaturated waste rock at an active mine and the other two tests aimed at evaluating the performance of thermally enhanced pipe installed in a fully saturated sedimentary rock formation. The concepts borrowed from pumping tests allowed the planning of the duration of the TRTs and the analysis of variable heat injection rate tests accounting for external heat transfer and temperature recovery, which reduces the uncertainty in the estimation of thermal properties.