USGS Science⌕ Search

Geology topics

P. A. Hsieh

Publications and source records attributed to P. A. Hsieh.

14 recordsLinked to original sources
USGS2011

Estimating Hydraulic Parameters When Poroelastic Effects Are Significant

For almost 80 years, deformation-induced head changes caused by poroelastic effects have been observed during pumping tests in multilayered aquifer-aquitard systems. As water in the aquifer is released from compressive storage during pumping, the aquifer is deformed both in the horizontal and vertical directions. This deformation in the pumped aquifer causes deformation in the adjacent layers, resulting in changes in pore pressure that may produce drawdown curves that differ significantly from those predicted by traditional groundwater theory. Although these deformation-induced head changes have been analyzed in several studies by poroelasticity theory, there are at present no practical guidelines for the interpretation of pumping test data influenced by these effects. To investigate the impact that poroelastic effects during pumping tests have on the estimation of hydraulic parameters, we generate synthetic data for three different aquifer-aquitard settings using a poroelasticity model, and then analyze the synthetic data using type curves and parameter estimation techniques, both of which are based on traditional groundwater theory and do not account for poroelastic effects. Results show that even when poroelastic effects result in significant deformation-induced head changes, it is possible to obtain reasonable estimates of hydraulic parameters using methods based on traditional groundwater theory, as long as pumping is sufficiently long so that deformation-induced effects have largely dissipated. ?? 2011 The Author(s). Journal compilation ?? 2011 National Ground Water Association.

Ground Water↗
USGS2009

Numerical models of caldera deformation: Effects of multiphase and multicomponent hydrothermal fluid flow

Ground surface displacement (GSD) in large calderas is often interpreted as resulting from magma intrusion at depth. Recent advances in geodetic measurements of GSD, notably interferometric synthetic aperture radar, reveal complex and multifaceted deformation patterns that often require complex source models to explain the observed GSD. Although hydrothermal fluids have been discussed as a possible deformation agent, very few quantitative studies addressing the effects of multiphase flow on crustal mechanics have been attempted. Recent increases in the power and availability of computing resources allow robust quantitative assessment of the complex time-variant thermal interplay between aqueous fluid flow and crustal deformation. We carry out numerical simulations of multiphase (liquid-gas), multicomponent (H 2O-CO2) hydrothermal fluid flow and poroelastic deformation using a range of realistic physical parameters and processes. Hydrothermal fluid injection, circulation, and gas formation can generate complex, temporally and spatially varying patterns of GSD, with deformation rates, magnitudes, and geometries (including subsidence) similar to those observed in several large calderas. The potential for both rapid and gradual deformation resulting from magma-derived fluids suggests that hydrothermal fluid circulation may help explain deformation episodes at calderas that have not culminated in magmatic eruption.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Determining anisotropic transmissivity using a simplified Papadopulos method

The straight-line method presented by Papadopulos requires a minimum of three observation wells for determining the transmissivity tensor of a homogeneous and anisotropic aquifer. A simplification of this method was developed for fractured aquifers where the principal directions of the transmissivity tensor are known prior to implementation, such as when fracture patterns on outcropping portions of the aquifer may be used to infer the principal directions. This new method assumes that observation wells are drilled along the two principal directions from the pumped well, thus reducing the required number of observation wells to two. This method was applied for an aquifer test in the fractured Navajo Sandstone of southwestern Utah and yielded minimum and maximum principal transmissivity values of 70 and 1800 m 2/d, respectively, indicating an anisotropy ratio of ???24 to 1.

Ground Water↗
USGS2005

Assessing an open-well aquifer test in fractured crystalline rock

Use of open wells to conduct aquifer tests in fractured crystalline rock aquifers is potentially problematic, because open wells can hydraulically connect highly permeable fracture zones at different depths within the rock. Because of this effect, it is questionable whether estimates of the hydraulic properties of the rock obtained from analysis of an open-well aquifer test in fractured crystalline rock are representative of the true rock properties. To address this issue, rock properties estimated from such a test are compared with rock properties estimated from an aquifer test at the same location in which packers were installed in wells to restrict vertical fluid movement along the wellbores. Results indicate that when the open-hole test is analyzed with a numerical model that accurately represents the three-dimensional heterogeneity structure in the wellfield, the open-hole test estimates of the rock hydraulic properties are similar to the packer-isolated test estimates. When the open-well test is analyzed using a model having an oversimplified representation of the rock heterogeneity, the estimates of rock transmissivity are similar to the transmissivity estimates from the packer-isolated test, but the storativity estimates are larger than the storativity estimates from the packer-isolated test. These results indicate the importance of accurately representing the rock heterogeneity structure when analyzing an open-hole aquifer test in fractured crystalline rock.Use of open wells to conduct aquifer tests in fractured crystalline rock aquifers is potentially problematic, because open wells can hydraulically connect highly permeable fracture zones at different depths within the rock. Because of this effect, it is questionable whether estimates of the hydraulic properties of the rock obtained from analysis of an open-well aquifer test in fractured crystalline rock are representative of the true rock properties. To address this issue, rock properties estimated from such a test are compared with rock properties estimated from an aquifer test at the same location in which packers were installed in wells to restrict vertical fluid movement along the wellbores. Results indicate that when the open-hole test is analyzed with a numerical model that accurately represents the three-dimensional heterogeneity structure in the wellfield, the open-hole test estimates of the rock hydraulic properties are similar to the packer-isolated test estimates. When the open-well test is analyzed using a model having an oversimplified representation of the rock heterogeneity, the estimates of rock transmissivity are similar to the transmissivity estimates from the packer-isolated test, but the storativity estimates are larger than the storativity estimates from the packer-isolated test. These results indicate the importance of accurately representing the rock heterogeneity structure when analyzing an open-hole aquifer test in fractured crystalline rock.

New Hampshire↗
USGS2005

How good are estimates of transmissivity from slug tests in fractured rock?

Slug tests in fractured rock usually are interpreted with models that assume homogeneous formation properties, even though hydraulic properties of fractures can vary by many orders of magnitude over the length of boreholes. To investigate the impact of heterogeneity on the interpretation of slug tests in fractured rock, slug tests were conducted over large intervals of boreholes in crystalline rock in central New Hampshire, and interpreted using a homogeneous model. The results of the slug tests were then compared with estimates of transmissivity from fluid-injection tests conducted over shorter intervals in the same boreholes. The fluid-injection tests showed transmissivity to vary more than six orders of magnitude over the length of the boreholes; however, the sum of the transmissivities from the fluid-injection tests were within an order of magnitude of the transmissivity estimated from the slug tests. Although the two estimates of transmissivity were within an order of magnitude of each other, the water level responses during the slug tests did not exactly match the responses predicted by the homogeneous model. To investigate the effect of heterogeneity on water level responses during slug tests, a Laplace-transform solution was developed for slug tests conducted in boreholes containing multiple fractures with hydraulic properties that vary over the length of the borehole. A comparison of this solution with the homogeneous model shows no difference between the shape of water level responses in a homogeneous formation and a (layered) heterogeneous formation. Furthermore, the transmissivity estimated using a homogeneous model is within an order of magnitude of the prescribed transmissivity in the heterogeneous model. Thus, differences between responses predicted from a homogeneous model and measured water levels during slug tests can be attributed to phenomena such as nonradial flow in the vicinity of the borehole, and not heterogeneous hydraulic properties over the length of the borehole. The experimental results of this investigation show that even when conditions such as nonradial flow are present in the vicinity of the borehole, interpretations of slug tests using a homogeneous model provided order-of-magnitude estimates of transmissivity in the crystalline rock terrane under consideration.Slug tests in fractured rock are usually interpreted with models that assume homogeneous formation properties, although hydraulic properties of fractures can vary by many order of magnitude over the length of the boreholes. To investigate the impact of heterogeneity on the interpretation of slug tests in fractured rocks, slug tests were conducted over large intervals of boreholes in crystalline rock in central New Hampshire. Tests results were interpreted using a homogeneous model. A Laplace-transform solution was developed for slug tests conducted in boreholes containing multiple fractures with hydraulic properties. Finally, a comparison was made between this solution and the homogeneous model.

Ground Water↗
USGS2005

Seasonal seismicity at western United States volcanic centers

We examine 20-yr data sets of seismic activity from 10 volcanic areas in the western United States for annual periodic signals (seasonality), focusing on large calderas (Long Valley caldera and Yellowstone) and stratovolcanoes (Cascade Range). We apply several statistical methods to test for seasonality in the seismic catalogs. In 4 of the 10 regions, statistically significant seasonal modulation of seismicity (> 90% probability) occurs, such that there is an increase in the monthly seismicity during a given portion of the year. In five regions, seasonal seismicity is significant in the upper 3 km of the crust. Peak seismicity occurs in the summer and autumn in Mt. St. Helens, Hebgen Lake/Madison Valley, Yellowstone Lake, and Mammoth Mountain. In the eastern south moat of Long Valley caldera (LVC) peak seismicity occurs in the winter and spring. We quantify the possible external forcing mechanisms that could modulate seasonal seismicity. Both snow unloading and groundwater recharge can generate large stress changes of > 5 kPa at seismogenic depths and may thus contribute to seasonality. ?? 2005 Elsevier B.V. All rights reserved.

Earth and Planetary Science Letters↗
USGS2002

Crosswell seismic investigation of hydraulically conductive, fracture bedrock near Mirror Lake, New Hampshire

Near Mirror Lake, New Hampshire (USA), hydraulically conductive, fractured bedrock was investigated with the crosswell seismic method to determine whether this method could provide any information about hydraulic conductivity between wells. To this end, crosswell seismic data, acoustic logs from boreholes, image logs from boreholes, and single borehole hydraulic tests were analyzed. The analysis showed that, first, the P-wave velocities from the acoustic logs tended to be higher in schist than they were in granite. (Schist and granite were the dominant rock types). Second, the P-wave velocities from the acoustic logs tended to be low near fractures. Third, the hydraulic conductivity was always low (always less than to 10 −8 m/s) where no fractures intersected the borehole, but the hydraulic conductivity ranged from low to high (from less than to 10 −10 m/s to 10 −4 m/s) where one or more fractures intersected the borehole. Fourth, high hydraulic conductivities were slightly more frequent when the P-wave velocity was low (less than 5200 m/s) than when it was high (greater than or equal to 5200 m/s). The interpretation of this statistical relation was that the fractures tended to increase the hydraulic conductivity and to lower the P-wave velocity. This statistical relation was applied to a velocity tomogram to create a map showing the probability of high hydraulic conductivity; the map was consistent with results from independent hydraulic tests.

New Hampshire↗
USGS2002

Evaluation of longitudinal dispersivity estimates from forced-gradient tracer tests in heterogeneous aquifers

Converging radial-flow and two-well tracer tests are simulated in two-dimensional aquifers to investigate the effects of heterogeneity and forced-gradient test configuration on longitudinal dispersivity (??L) estimates, and to compare ??L estimates from forced-gradient tests with ??L values that characterize solute spreading under natural-gradient flow. Results indicate that in both mildly and highly heterogeneous aquifers, ??L estimates from two-well tests are generally larger than those from radial-flow tests. In mildly heterogeneous aquifers, ??L estimates from two-well tests with relatively large tracer transport distances are similar to ??L values from natural-gradient simulations. In highly heterogeneous aquifers, ??L estimates from two-well tests at all tracer transport distances are typically smaller than ??L values from natural-gradient simulations.

IAHS-AISH Publication↗
USGS2002

Evaluation of longitudinal dispersivity estimates from forced-gradient tracer tests in heterogeneous aquifers

Converging radial-flow and two-well tracer tests are simulated in two-dimensional aquifers to investigate the effects of heterogeneity and forced-gradient test configuration on longitudinal dispersivity (??L) estimates, and to compare ??L estimates from forced-gradient tests with ??L values that characterize solute spreading under natural-gradient flow. Results indicate that in mildly heterogeneous aquifers, for tests with relatively large tracer transport distances, ??L estimates from the two test types are generally similar, and are also similar to ??L values determined from natural-gradient tracer simulations. In highly heterogeneous aquifers, ??L estimates from two-well tests are generally larger than those from radial-flow tests, and the ??L estimates from both test types are typically smaller than the ??L values determined from natural-gradient simulations.

Acta Universitatis Carolinae, Geologica↗
USGS1981

A transient laboratory method for determining the hydraulic properties of 'tight' rocks-I. Theory

Transient pulse testing has been employed increasingly in the laboratory to measure the hydraulic properties of rock samples with low permeability. Several investigators have proposed a mathematical model in terms of an initial-boundary value problem to describe fluid flow in a transient pulse test. However, the solution of this problem has not been available. In analyzing data from the transient pulse test, previous investigators have either employed analytical solutions that are derived with the use of additional, restrictive assumptions, or have resorted to numerical methods. In Part I of this paper, a general, analytical solution for the transient pulse test is presented. This solution is graphically illustrated by plots of dimensionless variables for several cases of interest. The solution is shown to contain, as limiting cases, the more restrictive analytical solutions that the previous investigators have derived. A method of computing both the permeability and specific storage of the test sample from experimental data will be presented in Part II.

International Journal of Rock Mechanics and Mining↗
USGS1981

A transient laboratory method for determining the hydraulic properties of 'tight' rocks-II. Application

In Part I a general analytical solution for the transient pulse test was presented. Part II presents a graphical method for analyzing data from a test to obtain the hydraulic properties of the sample. The general solution depends on both hydraulic conductivity and specific storage and, in theory, analysis of the data can provide values for both of these hydraulic properties. However, in practice, one of two limiting cases may apply in which case it is possible to calculate only hydraulic conductivity or the product of hydraulic conductivity times specific storage. In this paper we examine the conditions when both hydraulic parameters can be calculated. The analyses of data from two tests are presented. In Appendix I the general solution presented in Part I is compared with an earlier analysis, in which compressive storage in the sample is assumed negligible, and the error in calculated hydraulic conductivity due to this simplifying assumption is examined.

International Journal of Rock Mechanics and Mining↗