A method to allow temporal variation of velocity in travel-time tomography using microearthquakes induced during hydraulic fracturing

Abstract Hydraulic injections produce fluid-filled fractures that reduce the seismic velocity of the rock compared to intact rock. The travel times of microearthquakes induced by the injections may be used to discern changes in the rock velocities, as well as locating the microearthquakes. Determining the volumes of rock where the velocities have changed provides indirect evidence for the location of the injected fluid, and the character of the changes produced in the fractured rock. Available data are generally insufficient to resolve both the spatial and temporal changes within the rock. To extract information about temporal changes, and to obtain an improved image of the velocity structure, we chose a parameterization scheme in which the velocities of each block are allowed to change from the background velocity only after a threshold number of microearthquakes have occurred in the block. Regularizing by constraining the velocity of all the altered blocks to be similar helps stabilize the inversion. The regularization can be relaxed somewhat to allow the velocity of an altered block to be different from other altered blocks if the travel-time data are compelling. The parameterization scheme is justified since observations show that the volume of the seismically stimulated rock increases linearly with the volume of the injected fluid. We applied the method to data collected in a region of Precambrian crystalline rock that was injected with 21,600 m 3 of water. We use travel times from a total of 3886 microearthquakes that were induced by the injection. The mean RMS travel-time residual decreases about 7%. The velocity structure contains a low-velocity zone located near the injection region. Other distinct low-velocity zones are identified. The pattern of microearthquake locations found using our method appears to contain more structure than the pattern found in locations determined using a homogeneous velocity structure. Two clear low-velocity regions found near the point where water was injected into the rock are separated by a region whose velocity did not change. The region of unaltered velocity had a large number of microearthquakes.

[1]  James N. Albright,et al.  Interpretation of seismic data from hydraulic fracturing experiments at the Fenton Hill, New Mexico, hot dry rock geothermal site , 1982 .

[2]  Clifford H. Thurber,et al.  A fast algorithm for two-point seismic ray tracing , 1987 .

[3]  Leigh House,et al.  Detailed joint structure in a geothermal reservoir from studies of induced microearthquake clusters , 1997 .

[4]  Chuen Hon Arthur Cheng,et al.  Seismic imaging using microearthquakes induced by hydraulic fracturing , 1994 .

[5]  Gene Simmons,et al.  The effect of saturation on velocity in low porosity rocks , 1969 .

[6]  Michael Fehler,et al.  Stress control of seismicity patterns observed during hydraulic fracturing experiments at the Fenton Hill hot dry rock geothermal energy site, New Mexico , 1987 .

[7]  G. Pavlis,et al.  The mixed discrete‐continuous inverse problem: Application to the simultaneous determination of earthquake hypocenters and velocity structure , 1980 .

[8]  F. Harlow,et al.  A theoretical study of geothermal energy extraction , 1972 .

[9]  M. Fehler,et al.  Changes in compressional and shear wave velocities and dynamic moduli during operation of a hot dry rock geothermal system , 1983 .

[10]  Alberto Michelini,et al.  Seismological studies at Parkfield. I. Simultaneous inversion for velocity structure and hypocenters using cubic B-splines parameterization , 1991, Bulletin of the Seismological Society of America.

[11]  Leigh House,et al.  Locating microearthquakes induced by hydraulic fracturing in crystalline rock , 1987 .

[12]  M. Fehler Changes in P wave velocity during operation of a hot dry rock geothermal system , 1981 .

[13]  P. Dawson,et al.  Nonlinear teleseismic tomography at Long Valley Caldera, using three‐dimensional minimum travel time ray tracing , 1995 .