Earthquake location in strongly heterogeneous media

Abstract

Traveltime computation methods for strongly heterogeneous 3-D media developed during recent years are well suited for earthquake location. We present here a new method based on the traveltime algorithm of Podvin-Lecomte, related to the inverse problem formulation of Tarantola & Valette. The Podvin-Lecomte method, based on the Huygens principle, is very robust and allows arbitrary surface topography and station placement even for borehole instruments. First arrival traveltimes are computed for each of the recording stations using a fine 3-D velocity mesh (up to 106 cells on a workstation). The traveltime grid allows the use of the Tarantola & Valette formulation, which enables a full non-linear approach. The solution is given as a 3-D probability density function of hypocentre coordinates, which accounts for the arrival time measurements as well as a priori information for the location, the accuracy of both the arrival time readings and the computation of the theoretical traveltimes. This powerful method called 3dgridloc gives the location of the induced seismicity of the gas field of Lacq (France) using 443 520 cells of a 3-D velocity mesh and the observations from nine recording stations, one of which is located at the bottom of a 3880 m deep borehole. Location of synthetic foci as well as more than 500 actual earthquakes shows the real advantages of this new method versus the classical hypo71. A new insight into the induced seismicity is now possible: induced seismicity may occur as far away as 10 km from the gas reservoir and involve a much greater volume of rock than expected using earlier locations.