Geology topics

J.J. Butler

Publications and source records attributed to J.J. Butler.

23 records · Page 2Linked to original sources

USGS2001

1r2dinv: A finite-difference model for inverse analysis of two dimensional linear or radial groundwater flow

We have developed a program for inverse analysis of two-dimensional linear or radial groundwater flow problems. The program, 1r2dinv, uses standard finite difference techniques to solve the groundwater flow equation for a horizontal or vertical plane with heterogeneous properties. In radial mode, the program simulates flow to a well in a vertical plane, transforming the radial flow equation into an equivalent problem in Cartesian coordinates. The physical parameters in the model are horizontal or x-direction hydraulic conductivity, anisotropy ratio (vertical to horizontal conductivity in a vertical model, y-direction to x-direction in a horizontal model), and specific storage. The program allows the user to specify arbitrary and independent zonations of these three parameters and also to specify which zonal parameter values are known and which are unknown. The Levenberg-Marquardt algorithm is used to estimate parameters from observed head values. Particularly powerful features of the program are the ability to perform simultaneous analysis of heads from different tests and the inclusion of the wellbore in the radial mode. These capabilities allow the program to be used for analysis of suites of well tests, such as multilevel slug tests or pumping tests in a tomographic format. The combination of information from tests stressing different vertical levels in an aquifer provides the means for accurately estimating vertical variations in conductivity, a factor profoundly influencing contaminant transport in the subsurface. ?? 2001 Elsevier Science Ltd. All rights reserved.

Computers & Geosciences↗

USGS1999

Pumping tests in networks of multilevel sampling wells: Motivation and methodology

The identification of spatial variations in hydraulic conductivity ( K ) on a scale of relevance for transport investigations has proven to be a considerable challenge. Recently, a new field method for the estimation of interwell variations in K has been proposed. This method, hydraulic tomography, essentially consists of a series of short‐term pumping tests performed in a tomographic‐like arrangement. In order to fully realize the potential of this approach, information about lateral and vertical variations in pumping‐induced head changes (drawdown) is required with detail that has previously been unobtainable in the field. Pumping tests performed in networks of multilevel sampling (MLS) wells can provide data of the needed density if drawdown can accurately and rapidly be measured in the small‐diameter tubing used in such wells. Field and laboratory experiments show that accurate transient drawdown data can be obtained in the small‐diameter MLS tubing either directly with miniature fiber‐optic pressure sensors or indirectly using air‐pressure transducers. As with data from many types of hydraulic tests, the quality of drawdown measurements from MLS tubing is quite dependent on the effectiveness of well development activities. Since MLS ports of the standard design are prone to clogging and are difficult to develop, alternate designs are necessary to ensure accurate drawdown measurements. Initial field experiments indicate that drawdown measurements obtained from pumping tests performed in MLS networks have considerable potential for providing valuable information about spatial variations in hydraulic conductivity.

Water Resources Research↗

USGS1994

The use of slug tests to describe vertical variations in hydraulic conductivity

Multilevel slug tests provide one means of obtaining estimates of hydraulic conductivity on a scale of relevance for contaminant transport investigations. A numerical model is employed here to assess the potential of multilevel slug tests to provide information about vertical variations in hydraulic conductivity under conditions commonly faced in field settings. The results of the numerical simulations raise several important issues concerning the effectiveness of this technique. If the length of the test interval is of the order of the average layer thickness, considerable error may be introduced into the conductivity estimates owing to the effects of adjoining layers. The influence of adjoining layers is dependent on the aspect ratio (length of test interval/well radius) of the tesy interval and the flow properties of the individual layers. If a low-permeability skin is present at the well, the measured vertical variations will be much less than the actual variations, owing to the influence of the skin conductivity on the parameter estimates. A high-permeability skin can also produce apparent vertical variations that are much less than the actual, owing to water flowing vertically along the conductive skin. In cases where the test interval spans a number of layers, a slug test will yield an approximate thickness-weighted average of the hydraulic conductivities of the intersected layers. In most cases, packer circumvention should not be a major concern when packers of 0.75 m or longer are employed. Results of this study are substantiated by recently reported field tests that demonstrate the importance of well emplacement and development activities for obtaining meaningful estimates from a program of multilevel slug tests.

Journal of Hydrology↗

USGS1991

Pumping tests in non-uniform aquifers - The linear strip case

Many pumping tests are performed in geologic settings that can be conceptualized as a linear infinite strip of one material embedded in a matrix of differing flow properties. A semi-analytical solution is presented to aid the analysis of drawdown data obtained from pumping tests performed in settings that can be represented by such a conceptual model. Integral transform techniques are employed to obtain a solution in transform space that can be numerically inverted to real space. Examination of the numerically transformed solution reveals several interesting features of flow in this configuration. If the transmissivity of the strip is much higher than that of the matrix, linear and bilinear flow are the primary flow regimes during a pumping test. If the contrast between matrix and strip properties is not as extreme, then radial flow should be the primary flow mechanism. Sensitivity analysis is employed to develop insight into the controls on drawdown in this conceptual model and to demonstrate the importance of temporal and spatial placement of observations. Changes in drawdown are sensitive to the transmissivity of the strip for a limited time duration. After that time, only the total drawdown remains a function of strip transmissivity. In the case of storativity, both the total drawdown and changes in drawdown are sensitive to the storativity of the strip for a time of quite limited duration. After that time, essentially no information can be gained about the storage properties of the strip from drawdown data. An example analysis is performed using data previously presented in the literature to demonstrate the viability of the semi-analytical solution and to illustrate a general procedure for analysis of drawdown data in complex geologic settings. This example reinforces the importance of observation well placement and the time of data collection in constraining parameter correlation, a major source of the uncertainty that arises in the parameter estimation procedure.

Journal of Hydrology↗

USGS1988

Pumping tests in nonuniform aquifers - The radially symmetric case

Traditionally, pumping-test-analysis methodology has been limited to applications involving aquifers whose properties are assumed uniform in space. This work attempts to assess the applicability of analytical methodology to a broader class of units with spatially varying properties. An examination of flow behavior in a simple configuration consisting of pumping from the center of a circular disk embedded in a matrix of differing properties is the basis for this investigation. A solution describing flow in this configuration is obtained through Laplace-transform techniques using analytical and numerical inversion schemes. Approaches for the calculation of flow properties in conditions that can be roughly represented by this simple configuration are proposed. Possible applications include a wide variety of geologic structures, as well as the case of a well skin resulting from drilling or development. Of more importance than the specifics of these techniques for analysis of water-level responses is the insight into flow behavior during a pumping test that is provided by the large-time form of the derived solution. The solution reveals that drawdown during a pumping test can be considered to consist of two components that are dependent and independent of near-well properties, respectively. Such an interpretation of pumping-test drawdown allows some general conclusions to be drawn concerning the relationship between parameters calculated using analytical approaches based on curve-matching and those calculated using approaches based on the slope of a semilog straight line plot. The infinite-series truncation that underlies the semilog analytical approaches is shown to remove further contributions of near-well material to total drawdown. In addition, the semilog distance-drawdown approach is shown to yield an expression that is equivalent to the Thiem equation. These results allow some general recommendations to be made concerning observation-well placement for pumping tests in nonuniform aquifers. The relative diffusivity of material on either side of a discontinuity is shown to be the major factor in controlling flow behavior during the period in which the front of the cone of depression is moving across the discontinuity. Though resulting from an analysis of flow in an idealized configuration, the insights of this work into flow behavior during a pumping test are applicable to a wide class of nonuniform units.

Journal of Hydrology↗