Jean E. Elkhoury

Seismic waves increase permeability

Earthquakes have been observed to affect hydrological systems in a variety of ways—water well levels can change dramatically, streams can become fuller and spring discharges can increase at the time of earthquakes. Distant earthquakes may even increase the permeability in faults. Most of these hydrological observations can be explained by some form of permeability increase. Here we use the response of water well levels to solid Earth tides to measure permeability over a 20-year period. At the time of each of seven earthquakes in Southern California, we observe transient changes of up to 24° in the phase of the water level response to the dilatational volumetric strain of the semidiurnal tidal components of wells at the Piñon Flat Observatory in Southern California. After the earthquakes, the phase gradually returns to the background value at a rate of less than 0.1° per day. We use a model of axisymmetric flow driven by an imposed head oscillation through a single, laterally extensive, confined, homogeneous and isotropic aquifer to relate the phase response to aquifer properties. We interpret the changes in phase response as due to changes in permeability. At the time of the earthquakes, the permeability at the site increases by a factor as high as three. The permeability increase depends roughly linearly on the amplitude of seismic-wave peak ground velocity in the range of 0.21–2.1 cm s-1. Such permeability increases are of interest to hydrologists and oil reservoir engineers as they affect fluid flow and might determine long-term evolution of hydrological and oil-bearing systems. They may also be interesting to seismologists, as the resulting pore pressure changes can affect earthquakes by changing normal stresses on faults.

Changes in permeability caused by transient stresses: Field observations, experiments, and mechanisms

CHANGES IN PERMEABILITY CAUSED BY TRANSIENT STRESSES: FIELD OBSERVATIONS, EXPERIMENTS, AND MECHANISMS Michael Manga, 1 Igor Beresnev, 2 Emily E. Brodsky, 3 Jean E. Elkhoury, 4 Derek Elsworth, 5 S. E. Ingebritsen, 6 David C. Mays, 7 and Chi-Yuen Wang 1 Received 7 November 2011; revised 15 February 2012; accepted 10 March 2012; published 12 May 2012. [ 1 ] Oscillations in stress, such as those created by earth- quakes, can increase permeability and fluid mobility in geo- logic media. In natural systems, strain amplitudes as small as 10 A6 can increase discharge in streams and springs, change the water level in wells, and enhance production from petroleum reservoirs. Enhanced permeability typically recovers to prestimulated values over a period of months to years. Mechanisms that can change permeability at such small stresses include unblocking pores, either by breaking up permeability-limiting colloidal deposits or by mobilizing droplets and bubbles trapped in pores by capillary forces. The recovery time over which permeability returns to the prestimulated value is governed by the time to reblock pores, or for geochemical processes to seal pores. Monitor- ing permeability in geothermal systems where there is abun- dant seismicity, and the response of flow to local and regional earthquakes, would help test some of the proposed mechanisms and identify controls on permeability and its evolution. Citation: Manga, M., I. Beresnev, E. E. Brodsky, J. E. Elkhoury, D. Elsworth, S. E. Ingebritsen, D. C. Mays, and C.-Y. Wang (2012), Changes in permeability caused by transient stresses: Field observations, experiments, and mechanisms, Rev. Geophys., 50, RG2004, doi:10.1029/2011RG000382. INTRODUCTION [ 2 ] The permeability of Earth’s crust is of great interest because it largely governs key geologic processes such as advective transport of heat and solutes and the generation of elevated fluid pressures by processes such as physical com- paction, heating, and mineral dehydration. For an isotropic Department of Earth and Planetary Science, University of California, Berkeley, California, USA. Department of Geological and Atmospheric Sciences, Iowa State University, Ames, Iowa, USA. Department of Earth and Planetary Sciences, University of California, Santa Cruz, California, USA. Department of Civil and Environmental Engineering, University of California, Irvine, California, USA. Department of Energy and Mineral Engineering, Center for Geomechanics, Geofluids, and Geohazards, EMS Energy Institute, Pennsylvania State University, University Park, Pennsylvania, USA. U.S. Geological Survey, Menlo Park, California, USA. Department of Civil Engineering, University of Colorado Denver, Denver, Colorado, USA. Corresponding author: M. Manga, Department of Earth and Planetary Science, University of California, 307 McCone Hall, Berkeley, CA 94720, USA. (manga@seismo.berkeley.edu) material, permeability k is defined by Darcy’s law that relates the fluid discharge per unit area q to the gradient of hydraulic head h, q ¼A kgr rh; m where r is the fluid density, m the fluid viscosity and g is gravity. The permeability of common geologic media varies by approximately 16 orders of magnitude, from values as low as 10 A23 m 2 in intact crystalline rock, intact shales, and fault cores, to values as high as 10 A7 m 2 in well-sorted gravels. Nevertheless, despite being highly heterogeneous, perme- ability can be characterized at the crustal scale in a manner that provides useful insight [e.g., Gleeson et al., 2011]. [ 3 ] The responses of hydrologic systems to deformation provide some insight into controls on permeability, in par- ticular its evolution in time. For example, the water level in wells and discharge in rivers have both been observed to change after earthquakes. Because earthquakes produce stresses that can change hydrogeologic properties of the crust, hydrologic responses to earthquakes are expected, especially in the near field (within a fault length of the Copyright 2012 by the American Geophysical Union. Reviews of Geophysics, 50, RG2004 / 2012 1 of 24 Paper number 2011RG000382 8755-1209/12/2011RG000382 RG2004

Laboratory observations of permeability enhancement by fluid pressure oscillation of in situ fractured rock

Received 5 June 2010; revised 5 November 2010; accepted 22 December 2010; published 24 February 2011. [1] We report on laboratory experiments designed to investigate the influence of pore pressure oscillations on the effective permeability of fractured rock. Berea sandstone samples were fractured in situ under triaxial stresses of tens of megapascals, and deionized water was forced through the incipient fracture under conditions of steady and oscillating pore pressure. We find that short‐term pore pressure oscillations induce long‐term transient increases in effective permeability of the fractured samples. The magnitude of the effective permeability enhancements scales with the amplitude of pore pressure oscillations, and changes persist well after the stress perturbation. The maximum value of effective permeability enhancement is 5 × 10 −16 m 2 with a background permeability of 1 × 10 −15 m 2 ; that is, the maximum enhanced permeability is 1.5 × 10 −15 m 2 . We evaluate poroelastic effects and show that hydraulic storage release does not explain our observations. Effective permeability recovery following dynamic oscillations occurs as the inverse square root of time. The recovery indicates that a reversible mechanism, such as clogging/unclogging of fractures, as opposed to an irreversible one, like microfracturing, is responsible for the transient effective permeability increase. Our work suggests the feasibility of dynamically controlling the effective permeability of fractured systems. The result has consequences for models of earthquake triggering and permeability enhancement in fault zones due to dynamic shaking from near and distant earthquakes.

Fracture Permeability Alteration due to Chemical and Mechanical Processes: A Coupled High-Resolution Model

Reactive fluid-flow experiments in fractures subjected to normal stress suggest the potential for either increased or decreased permeability resulting from fracture-surface dissolution. We present a computational model that couples mechanical deformation and chemical alteration of fractures subjected to constant normal stress and reactive fluid flow. The model explicitly represents micro-scale roughness of the fracture surfaces and calculates elastic deformation of the rough surfaces using a semi-analytical approach that ensures the surfaces remain in static equilibrium. A depth-averaged reactive transport model calculates chemical alteration of the surfaces, which leads to alteration of the contacting fracture surfaces. The mechanical deformation and chemical alteration calculations are explicitly coupled, which is justified by the disparate timescales required for equilibration of mechanical stresses and reactive transport processes. An idealized analytical representation of dissolution from a single contacting asperity shows that under reaction-limited conditions, contacting asperities can dissolve faster than the open regions of the fracture. Computational simulations in fractures with hundreds of contacting asperities show that the transition from transport-limited conditions (low flow rates) to reaction-rate-limited conditions (high flow rates) causes a shift from monotonically increasing permeability to a more complicated process in which permeability initially decreases and then increases as contacting asperities begin to dissolve. These results are qualitatively consistent with a number of experimental observations reported in the literature and suggest the potential importance of the relative magnitude of mass transport and reaction kinetics on the evolution of fracture permeability in fractures subjected to combined normal stress and reactive fluid flow.

Downscaling of fracture energy during brittle creep experiments

We present mode 1 brittle creep fracture experiments along fracture surfaces that contain strength heterogeneities. Our observations provide a link between smooth macroscopic time-dependent failure and intermittent microscopic stress-dependent processes. We find the large-scale response of slow-propagating subcritical cracks to be well described by an Arrhenius law that relates the fracture speed to the energy release rate. At the microscopic scale, high-resolution optical imaging of the transparent material used (PMMA) allows detailed description of the fracture front. This reveals a local competition between subcritical and critical propagation (pseudo stick-slip front advances) independently of loading rates. Moreover, we show that the local geometry of the crack front is self-affine and the local crack front velocity is power law distributed. We estimate the local fracture energy distribution by combining high-resolution measurements of the crack front geometry and an elastic line fracture model. We show that the average local fracture energy is significantly larger than the value derived from a macroscopic energy balance. This suggests that homogenization of the fracture energy is not straightforward and should be taken cautiously. Finally, we discuss the implications of our results in the context of fault mechanics.