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

TL;DR

This paper introduces a matrix-based seismic imaging method that compensates for wave velocity heterogeneities, enabling high-resolution 3D imaging of the Earth's subsurface even with limited prior velocity information.

Contribution

A novel distortion matrix framework is developed to improve passive seismic imaging resolution without requiring detailed background velocity models.

Findings

Achieved 80 m resolution at 4 km depth in the San Jacinto Fault Zone.

Resolution surpasses the diffraction limit by a factor of eight.

Subsoil heterogeneities enhance imaging by acting as scattering lenses.

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.