A compressive sensing framework for seismic source parameter estimation

作者:Rodriguez Ismael Vera*; Sacchi Mauricio; Gu Yu J
来源:Geophysical Journal International, 2012, 191(3): 1226-1236.
DOI:10.1111/j.1365-246X.2012.05659.x

摘要

Simultaneous estimation of origin time, location and moment tensor of seismic events is critical for automatic, continuous, real-time monitoring systems. Recent studies have shown that such systems can be implemented via waveform fitting methods based on pre-computed catalogues of Greens functions. However, limitations exist in the number and length of the recorded traces, and the size of the monitored volume that these methods can handle without compromising real-time response. This study presents numerical tests using a novel waveform fitting method based on compressive sensing, a field of applied mathematics that provides conditions for sampling and recovery of signals that admit a sparse representation under a known base or dictionary. Compressive sensing techniques enable us to determine source parameters in a compressed space, where the dimensions of the variables involved in the inversion are significantly reduced. Results using a hypothetical monitoring network with a dense number of recording stations show that a compressed catalogue of Greens functions with 0.004 per cent of its original size recovers the exact source parameters in more than 50 per cent of the tests. The gains in processing time in this case drop from an estimated 90 days to browse a solution in the uncompressed catalogue to 41.57 s to obtain an estimation using the compressed catalogue. For simultaneous events, the compressive sensing approach does not appear to influence the estimation results beyond the limitations presented by the uncompressed case. The main concern in the use of compressive sensing is detectability issues observed when the amount of compression is beyond a minimum value that is identifiable through numerical experiments. Tests using real data from the 2002 June 18 Caborn Indiana earthquake show that the presence of noise and inaccurate Greens functions require a smaller amount of compression to reproduce the solution obtained with the uncompressed catalogue. In this case, numerical simulation enables the assessment of the amount of compression that provides a reasonable rate of detectability. Overall, the numerical experiments demonstrate the effectiveness of our compressed domain inversion method in the real-time monitoring of seismic sources with dense networks of receivers. As an added benefit of the compression process, the size of the monitored volume can also be increased under specific restrictions while maintaining the real-time response.

  • 出版日期2012-12