skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Do Simple Analytical Models Capture Complex Fractured Bedrock Hydraulics? Oscillatory Flow Tests Suggest Not
Abstract Fractured sedimentary bedrock aquifers represent complex flow systems that may contain fast, fracture‐dominated flow paths and slower, porous media‐dominated flow paths. Thus, characterizing the dynamics of flow and transport through these aquifers remains a fundamental hydrogeologic challenge. Recent studies have demonstrated the utility of a novel hydraulic testing approach, oscillatory flow testing, in field settings to characterize single bedrock fractures embedded in low‐porosity sedimentary bedrock. These studies employed an idealized analytical model assuming Darcian flow through a nondeforming, constant‐aperture, nonleaky fracture for data interpretation, and reported period‐dependent effective fracture flow parameters. Here, we present the application of oscillatory flow testing across a range of frequencies and inter‐well spacings on a fracture embedded in poorly cemented sedimentary bedrock with considerable primary porosity at the Field Site for Research in Fractured Sedimentary Rock. Consistent with previous studies, we show an apparent period‐dependence in returned flow parameters, with hydraulic diffusivity decreasing and storativity increasing with increasing oscillation period, when assuming an idealized fracture conceptual model. We present simple analyses that examine non‐Darcian flow and borehole storage effects as potential test design artifacts and a simple analytical model that examines fluid leakage to the surrounding host rock as a potential hydraulic mechanism that might contribute to the period‐dependent flow parameters. These analyses represent a range of conceptual assumptions about fracture behavior during hydraulic testing, none of which account for the measured responses during oscillatory flow testing, leading us to argue that other hydraulic processes (e.g., aperture heterogeneity and/or fracture hydromechanics) are necessary to accurately represent pressure propagation through fractured sedimentary bedrock.  more » « less
Award ID(s):
1654649
PAR ID:
10397408
Author(s) / Creator(s):
;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Groundwater
Volume:
61
Issue:
6
ISSN:
0017-467X
Format(s):
Medium: X Size: p. 816-833
Size(s):
p. 816-833
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Distributed acoustic sensing (DAS) was originally intended to measure oscillatory strain at frequencies of 1 Hz or more on a fiber optic cable. Recently, measurements at much lower frequencies have opened the possibility of using DAS as a dynamic strain sensor in boreholes. A fiber optic cable mechanically coupled to a geologic formation will strain in response to hydraulic stresses in pores and fractures. A DAS interrogator can measure dynamic strain in the borehole, which can be related to fluid pressure through the mechanical compliance properties of the formation. Because DAS makes distributed measurements, it is capable of both locating hydraulically active features and quantifying the fluid pressure in the formation. We present field experiments in which a fiber optic cable was mechanically coupled to two crystalline rock boreholes. The formation was stressed hydraulically at another well using alternating injection and pumping. The DAS instrument measured oscillating strain at the location of a fracture zone known to be hydraulically active. Rock displacements of less than 1 nm were measured. Laboratory experiments confirm that displacement is measured correctly. These results suggest that fiber optic cable embedded in geologic formations may be used to map hydraulic connections in three‐dimensional fracture networks. A great advantage of this approach is that strain, an indirect measure of hydraulic stress, can be measured without beforehand knowledge of flowing fractures that intersect boreholes. The technology has obvious applications in water resources, geothermal energy, CO2sequestration, and remediation of groundwater in fractured bedrock. 
    more » « less
  2. ABSTRACT:Creation of a fracture network in a hydraulic fracturing process is essential for subsurface energy extraction and CO2 sequestration. It is facilitated by reactivation of pre-existing intersecting weak layers and cemented cracks in the rock. In this study, a poromechanical model is developed for the hydraulic fracturing process in rocks containing such pre-existing weak layers. Based on the mixture theory, the crack band model is used to simulate the growth of a crack system. The governing equations with the parameters for hydromechanical coupling are derived, to describe the evolution of the opening and branching of cracks caused by water injection. Microplane model M7 is adopted to characterize the deformation and fracturing of the solid skeleton of the rock, and the Poiseuille law is used to characterize fluid flow through the hydraulic fractures. Numerical simulations are performed to reproduce and interpret recently published laboratory-scale hydraulic fracturing experiments conducted at Los Alamos National Laboratory (LANL). In these experiments, the rock was represented by confined plaster slabs containing orthogonal intersecting weak layers of higher porosity. Numerical simulations reveal how poromechanical characteristics such as the Biot coefficient and the fluid injection rate lead to various typical fracture modes observed in the experiments. These modes include formation of one dominant planar crack or various orthogonal fracture networks. 
    more » « less
  3. Fractures in Earth's critical zone influence groundwater flow and storage and promote chemical weathering. Fractured materials are difficult to characterize on large spatial scales because they contain fractures that span a range of sizes, have complex spatial distributions, and are often inaccessible. Therefore, geophysical characterizations of the critical zone depend on the scale of measurements and on the response of the medium to impulses at that scale. Using P-wave velocities collected at two scales, we show that seismic velocities in the fractured bedrock layer of the critical zone are scale-dependent. The smaller-scale velocities, derived from sonic logs with a dominant wavelength of ~0.3 m, show substantial vertical and lateral heterogeneity in the fractured rock, with sonic velocities varying by 2,000 m/s over short lateral distances (~20 m), indicating strong spatial variations in fracture density. In contrast, the larger-scale velocities, derived from seismic refraction surveys with a dominant wavelength of ~50 m, are notably slower than the sonic velocities (a difference of ~3,000 m/s) and lack lateral heterogeneity. We show that this discrepancy is a consequence of contrasting measurement scales between the two methods; in other words, the contrast is not an artifact but rather information—the signature of a fractured medium (weathered/fractured bedrock) when probed at vastly different scales. We explore the sample volumes of each measurement and show that surface refraction velocities provide reliable estimates of critical zone thickness but are relatively insensitive to lateral changes in fracture density at scales of a few tens of meters. At depth, converging refraction and sonic velocities likely indicate the top of unweathered bedrock, indicative of material with similar fracture density across scales. 
    more » « less
  4. ABSTRACT Analysis of PRISM and SNOTEL station data paired with USGS streamflow gage data in the western United States shows that, in snow‐dominated mountainous watersheds, streamflow regimes differ between watersheds with karst geology and their non‐karst neighbours. These carbonate aquifers exhibit a spectrum of flow paths encompassing karst conduits, including large fractures or voids that transmit water readily to springs and other surface waters, and matrix flow paths through soils, highly fractured bedrock, or porous media bedrock grains. A well‐connected karst aquifer will discharge a large portion of its accumulated precipitation to surface water via springs and other groundwater flow paths on an annual scale, exhibiting a lagged response to precipitation presenting as a “memory effect” in hydrograph time series. These patterns were observed in the hydrologic records of gaged watersheds with exposed or near‐surface carbonate layers accounting for > 30% of their drainage area. In western snow‐dominated watersheds, where paired streamflow and SNOTEL data are available, analysis of the precipitation and flow time series shows low‐flow volume is strongly related to karst aquifer conditions and winter precipitation when compared to low‐flow volumes present in non‐karst watersheds, which have a complex relationship to multiple driving metrics. Analysis of normalised streamflow and cumulative precipitation in karst watersheds show that low‐flow conditions are highly dependent on the preceding winter precipitation and streamflow in both wet and dry periods. In non‐karst watersheds, increased precipitation primarily impacts high‐flow, spring runoff volumes with no clear relationship to low‐flow periods. When comparing cumulative streamflow and precipitation volumes within each water year and over longer timescales, karst watersheds show the potential filling and draining of large amounts of karst storage, whereas non‐karst watersheds demonstrate a more stable storage regime. Communities in many western US watersheds are dependent on snow‐dominated karst watersheds for their water supply. This analysis, using widely available hydrologic data, can provide insight into the recharge and storage processes within these watersheds, improve our ability to assess current flow regimes, anticipate the impacts of climate change on water availability, and help manage water supplies. 
    more » « less
  5. Abstract Earth's Critical Zone (CZ), the near‐surface layer where rock is weathered and landscapes co‐evolve with life, is profoundly influenced by the type of underlying bedrock. Previous studies employing the CZ framework have focused primarily on landscapes dominated by silicate rocks. However, carbonate rocks crop out on approximately 15% of Earth's ice‐free continental surface and provide important water resources and ecosystem services to ∼1.2 billion people. Unlike silicates, carbonate minerals weather congruently and have high solubilities and rapid dissolution kinetics, enabling the development of large, interconnected pore spaces and preferential flow paths that restructure the CZ. Here we review the state of knowledge of the carbonate CZ, exploring parameters that produce contrasts in the CZ in different carbonate settings and identifying important open questions about carbonate CZ processes. We introduce the concept of a carbonate‐silicate CZ spectrum and examine whether current conceptual models of the CZ, such as the conveyor model, can be applied to carbonate landscapes. We argue that, to advance beyond site‐specific understanding and develop a more general conceptual framework for the role of carbonates in the CZ, we need integrative studies spanning both the carbonate‐silicate spectrum and a range of carbonate settings. 
    more » « less