Geodetic imaging: reservoir monitoring using satellite interferometry

SUMMARY Fluid fluxes within subsurface reservoirs give rise to surface displacements, particularly over periods of a year or more. Observations of such deformation provide a powerful tool for mapping fluid migration within the Earth, providing new insights into reservoir dynamics. In this paper we use Interferometric Synthetic Aperture Radar (InSAR) range changes to infer subsurface fluid volume strain at the Coso geothermal field. Furthermore, we conduct a complete model assessment, using an iterative approach to compute model parameter resolution and covariance matrices. The method is a generalization of a Lanczos-based technique which allows us to include fairly general regularization, such as roughness penalties. We find that we can resolve quite detailed lateral variations in volume strain both within the reservoir depth range (0.4‐2.5 km) and below the geothermal production zone (2.5‐5.0 km). The fractional volume change in all three layers of the model exceeds the estimated model parameter uncertainty by a factor of two or more. In the reservoir depth interval (0.4‐2.5 km), the predominant volume change is associated with northerly and westerly oriented faults and their intersections. However, below the geothermal production zone proper [the depth range 2.5‐5.0 km], there is the suggestion that both north- and northeast-trending faults may act as conduits for fluid flow.

[1]  Huajian Gao,et al.  Dislocations in inhomogeneous media via a moduli perturbation approach: General formulation and two‐dimensional solutions , 1994 .

[2]  Frédérique Rémy,et al.  Ice flow physical processes derived from the ERS-1 high-resolution map of the Antarctica and Greenland ice sheets , 1999 .

[3]  R. W. Ward,et al.  Three‐dimensional Q −1 model of the Coso Hot Springs Known Geothermal Resource Area , 1980 .

[4]  C. Bacon,et al.  Rhyolite thermobarometry and the shallowing of the magma reservoir, Coso volcanic field, California , 2000 .

[5]  W. Ellsworth,et al.  Teleseismic evidence for a low‐velocity body under the Coso Geothermal Area , 1980 .

[6]  Didier Massonnet,et al.  Land subsidence caused by the East Mesa Geothermal Field, California, observed using SAR interferometry , 1997 .

[7]  Takuo Maruyama,et al.  Static elastic dislocation in an infinite and semi-infinite medium , 1964 .

[8]  Didier Massonnet,et al.  Two examples of the use of SAR interferometry on displacement fields of small spatial extent , 1996 .

[9]  R. Kwok,et al.  Greenland Ice Sheet Surface Properties and Ice Dynamics from ERS-1 SAR Imagery , 1993, Science.

[10]  Satoshi Fujiwara,et al.  Coseismic crustal deformations of 1994 Northridge, California, earthquake detected by interferometric JERS 1 synthetic aperture radar , 1996 .

[11]  D. Castillo,et al.  Deep hydraulic fracture imaging: Recent advances in tiltmeter technologies , 1997 .

[12]  K. Feigl,et al.  The displacement field of the Landers earthquake mapped by radar interferometry , 1993, Nature.

[13]  Daniel Dzurisin,et al.  Renewed uplift at the Yellowstone Caldera measured by leveling surveys and satellite radar interferometry , 1999 .

[14]  D. Plouff,et al.  Aeromagnetic and gravity surveys in the Coso Range, California , 1980 .

[15]  K. Feigl,et al.  Discrimination of geophysical phenomena in satellite radar interferograms , 1995 .

[16]  Charles Werner,et al.  On the derivation of coseismic displacement fields using differential radar interferometry: The Landers earthquake , 1994 .

[17]  R. Goldstein,et al.  Topographic mapping from interferometric synthetic aperture radar observations , 1986 .

[18]  A. R. Mitchell,et al.  The Finite Difference Method in Partial Differential Equations , 1980 .

[19]  M. Bruno,et al.  Cost-Effective Monitoring of Injected Steam Migration Using Surface Deformation Analysis , 1994 .

[20]  Jonathan M. Lees,et al.  Three-dimensional P and S wave velocity structures of the Coso Geothermal Area, California, from microseismic travel time data , 1999 .

[21]  R. G. Pratt,et al.  Traveltime tomography in anisotropic media—II. Application , 1992 .

[22]  C. Weaver,et al.  Seismicity of the Coso Range, California , 1980 .

[23]  G. R. Roquemore.  Structure, tectonics, and stress field of the Coso Range, Inyo County, California , 1980 .

[24]  I. Stakgold.  Green's Functions and Boundary Value Problems , 1979 .

[25]  M. A. Chinnery,et al.  Elastic Dislocations in a Layered Half-Space—I. Basic Theory and Numerical Methods , 1974 .

[26]  G. B. Dalrymple,et al.  Late Cenozoic volcanism, geochronology, and structure of the Coso Range, Inyo County, California , 1980 .

[27]  Thatcher,et al.  Migration of fluids beneath yellowstone caldera inferred from satellite radar interferometry , 1998, Science.

[28]  A. Tarantola.  Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[29]  D. Vasco,et al.  Inversion for sources of crustal deformation and gravity change at the Yellowstone caldera , 1990 .

[30]  K. Mogi.  Relations between the Eruptions of Various Volcanoes and the Deformations of the Ground Surfaces around them , 1958 .

[31]  Howard A. Zebker,et al.  Seasonal subsidence and rebound in Las Vegas Valley, Nevada, observed by Synthetic Aperture Radar Interferometry , 2001 .

[32]  M. Walck.  Three‐dimensional Vp /Vs variations for the Coso Region, California , 1988 .

[33]  P. Segall,et al.  Volume strain within The Geysers geothermal field , 1999 .

[34]  Jonathan M. Lees,et al.  P wave anisotropy, stress, and crack distribution at Coso geothermal field, California , 1999 .

[35]  C. Lanczos.  An iteration method for the solution of the eigenvalue problem of linear differential and integral operators , 1950 .

[36]  Don W. Vasco,et al.  Monitoring of Fluid Injection and Soil Consolidation Using Surface Tilt Measurements , 1998 .

[37]  Gene H. Golub,et al.  Matrix computations , 1983 .

[38]  J. P. Henry,et al.  Laboratory investigation of the mechanical behaviour of Tournemire shale , 1997 .

[39]  Mark Simons,et al.  Deformation and seismicity in the Coso geothermal area, Inyo County, California: Observations and modeling using satellite radar interferometry , 2000 .

[40]  D. Wood,et al.  The geometry of a large-scale nitrogen gas hydraulic fracture formed in Devonian shale: an example of fracture mapping with tiltmeters , 1982 .

[41]  R. Goldstein,et al.  Mapping small elevation changes over large areas: Differential radar interferometry , 1989 .

[42]  L. A. Wright.  Late Cenozoic fault patterns and stress fields in the Great Basin and westward displacement of the Sierra Nevada block , 1976 .

[43]  R. Goldstein,et al.  Satellite Radar Interferometry for Monitoring Ice Sheet Motion: Application to an Antarctic Ice Stream , 1993, Science.

[44]  K. Feigl,et al.  Radar interferometry and its application to changes in the Earth's surface , 1998 .

[45]  D. Vasco,et al.  Global Earth structure: inference and assessment , 2002 .

[46]  W. Menke.  Geophysical data analysis : discrete inverse theory , 1984 .

[47]  Jeffrey R. Ridgway,et al.  The development of a deep-towed gravity meter, and its use in marine geophysical surveys of offshore southern California and an airborne laser altimeter survey of Long Valley, California , 1997 .

[48]  J. D. Eshelby.  The determination of the elastic field of an ellipsoidal inclusion, and related problems , 1957, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[49]  R. Clayton,et al.  Applications of attenuation tomography to Imperial Valley and Coso‐Indian Wells Region, southern California , 1988 .

[50]  R. Parker.  Geophysical Inverse Theory , 1994 .

[51]  Freysteinn Sigmundsson,et al.  Crustal Deformation from 1992 to 1995 at the Mid-Atlantic Ridge, Southwest Iceland, Mapped by Satellite Radar Interferometry , 1997, Science.

[52]  Kenneth W. Hudnut,et al.  Detection of aquifer system compaction and land subsidence using interferometric synthetic aperture radar, Antelope Valley, Mojave Desert, California , 1998 .

[53]  H. Kanamori,et al.  Anomalous Shear Wave Attenuation in the Shallow Crust Beneath the Coso Volcanic Region, California , 1988 .

[54]  R. Bamler,et al.  Synthetic aperture radar interferometry , 1998 .

[55]  Kazuaki Nakamura,et al.  Distribution of Quaternary rhyolite domes of the Coso Range, California: implications for extent of the geothermal anomaly. , 1980 .

[56]  D. Massonnet,et al.  Deflation of Mount Etna monitored by spaceborne radar interferometry , 1995, Nature.

[57]  J. Moore,et al.  Integrated mineralogical and fluid inclusion study of the Coso geothermal systems, California , 1996 .

[58]  Leen Weijers,et al.  Surface Tiltmeter Fracture Mapping Reaches New Depths - 10,000 Feet and Beyond? , 1998 .

[59]  James H. Dieterich,et al.  Finite element modeling of surface deformation associated with volcanism , 1975 .

[60]  I. Smith.  Programming the finite element method: With application to geomechanics , 1982 .

[61]  Jonathan M. Lees,et al.  Microseismicity, stress, and fracture in the Coso geothermal field, California , 1998 .

[62]  R. Clayton,et al.  P wave velocity variations in the Coso Region, California, derived from local earthquake travel times , 1987 .

[63]  Don W. Vasco,et al.  Using surface deformation to image reservoir dynamics , 2000 .

[64]  Christophe Delacourt,et al.  Post‐eruptive deformation associated with the 1986–87 and 1989 lava flows of Etna detected by radar interferometry , 1997 .

[65]  N. E. Goldstein,et al.  Using surface displacement and strain observations to determine deformation at depth, with an application to Long Valley Caldera, California , 1988 .

[66]  Hubert Fabriol,et al.  Monitoring and modeling land subsidence at the Cerro Prieto Geothermal Field, Baja California, Mexico, using SAR interferometry , 1999 .

[67]  Kurt L. Feigl,et al.  Estimation of an earthquake focal mechanism from a satellite radar interferogram: Application to the December 4, 1992 Landers aftershock , 1995 .

[68]  Sarva Jit Singh.  Static deformation of a multilayered half‐space by internal sources , 1970 .