Seismic wavefield reconstruction using a pre-conditioned wavelet–curvelet compressive sensing approach
- Creators
- Muir, Jack B.
- Zhan, Zhongwen
Abstract
The proliferation of large seismic arrays have opened many new avenues of geophysical research; however, most techniques still fundamentally treat regional and global scale seismic networks as a collection of individual time-series rather than as a single unified data product. Wavefield reconstruction allows us to turn a collection of individual records into a single structured form that treats the seismic wavefield as a coherent 3-D or 4-D entity. We propose a split processing scheme based on a wavelet transform in time and pre-conditioned curvelet-based compressive sensing in space to create a sparse representation of the continuous seismic wavefield with smooth second-order derivatives. Using this representation, we illustrate several applications, including surface wave gradiometry, Helmholtz–Hodge decomposition of the wavefield into irrotational and solenoidal components, and compression and denoising of seismic records.
Additional Information
© The Author(s) 2021. Published by Oxford University Press on behalf of The Royal Astronomical Society. This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model). Accepted 2021 June 3. Received 2021 May 23; in original form 2021 February 21. Published: 07 June 2021. The authors would like to thank the editor, Andrew Valentine, and an anonymous reviewer, for their insightful reviews that clarified several points. They would also like to thank Charles Langston for a particularly enthusiastic review. JBM acknowledges the financial support of the Origin Energy Foundation and the General Sir John Monash Foundation during his PhD studies. ZZ thanks the support from NSF CAREER award EAR 1848166 and the NSF/IUCRC GMG Center. Seismic data were processed using Obspy (Beyreuther et al. 2010), and figures were created using Matplotlib (Hunter 2007) using Cartopy for mapping (Met Office 2010). Data Availability: Data for the Mw 7.0 Loyalty Islands Earthquake were obtained through the obspy FDSN service using the Southern California Earthquake Data Center as a provider (SCEDC 2013). CSN data for the Mw 7.1 Ridgecrest Earthquake are available through the CSN website at http://csn.caltech.edu/data/. Code and data to recreate the Helmholtz–Hodge decomposition are available at https://github.com/jbmuir/HelmholtzHodgeCSN.Attached Files
Published - ggab222.pdf
Supplemental Material - ggab222_supplemental_files.zip
Files
Name | Size | Download all |
---|---|---|
md5:97d83f93833b3c7b9f1aee80a1bd45e3
|
4.1 MB | Preview Download |
md5:118c586b67317d07a2d44d4a0a27384c
|
20.5 MB | Preview Download |
Additional details
- Eprint ID
- 111716
- Resolver ID
- CaltechAUTHORS:20211102-203603665
- Origin Energy Foundation
- General Sir John Monash Foundation
- EAR-1848166
- NSF
- Created
-
2021-11-02Created from EPrint's datestamp field
- Updated
-
2021-11-02Created from EPrint's last_modified field
- Caltech groups
- Seismological Laboratory, Division of Geological and Planetary Sciences