A Distortion Matrix Framework for High-Resolution Passive Seismic 3-D Imaging: Application to the San Jacinto Fault Zone, California

Abstract

Reflection seismic imaging usually suffers from a loss of resolution and contrast because of the fluctuations of the wave velocities in the Earth's crust. In the literature, phase distortion issues are generally circumvented by means of a background wave velocity model. However, it requires a prior tomography of the wave velocity distribution in the medium, which is often not possible, especially in depth. In this paper, a matrix approach of seismic imaging is developed to retrieve a three-dimensional image of the subsoil, despite a rough knowledge of the background wave velocity. To do so, passive noise cross-correlations between geophones of a seismic array are investigated under a matrix formalism. More precisely, the detrimental effect of wave velocity fluctuations on imaging is overcome by introducing a novel mathematical object: The distortion matrix. This operator essentially connects any virtual source inside the medium with the distortion that a wavefront, emitted from that point, experiences due to heterogeneities. A time reversal analysis of the distortion matrix enables the estimation of the transmission matrix that links each real geophone at the surface and each virtual geophone in depth. Phase distortions can then be compensated for any point of the underground. Applied to seismic data recorded along the Clark branch of the San Jacinto fault zone, the present method is shown to provide an image of the fault until a depth of 4 km with a transverse resolution of 80 m. Strikingly, this resolution is almost one eighth below the diffraction limit imposed by the geophone array aperture. The heterogeneities of the subsoil play the role of a scattering lens and of a transverse wave guide which increase drastically the array aperture.

24 pages, 12 figures

References in corpus (6)

• Non-invasive real-time imaging through scattering layers and around corners via speckle correlations
• Detection and imaging in a random medium: A matrix method to overcome multiple scattering and aberration
• A random matrix approach to detect defects in a strongly scattering polycrystal: how the memory effect can help overcome multiple scattering
• Matrix Approach of Seismic Imaging: Application to the Erebus Volcano, Antarctica
• Ultrasound Matrix Imaging-Part II: The distortion matrix for aberration correction over multiple isoplanatic patches
• Ultrasound Matrix Imaging-Part I: The focused reflection matrix, the F-factor and the role of multiple scattering

Cited by in corpus (4)

• Three-Dimensional Ultrasound Matrix Imaging
• Distribution of seismic scatterers in the San Jacinto Fault Zone, southeast of Anza, California, based on passive matrix imaging
• Matrix imaging as a tool for high-resolution monitoring of deep volcanic plumbing systems with seismic noise
• Source and defect localization in thin elastic plates of arbitrary geometry using eigenmodes