Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression
Abstract
Full-3D seismic waveform tomography (F3DT) is the latest seismic tomography technique that can assimilate broadband, multi-component seismic waveform observations into high-resolution 3D subsurface seismic structure models. The main drawback in the current F3DT implementation, in particular the scattering-integral implementation (F3DT-SI), is the high disk storage cost and the associated I/O overhead of archiving the 4D space-time wavefields of the receiver- or source-side strain tensors. The strain tensor fields are needed for computing the data sensitivity kernels, which are used for constructing the Jacobian matrix in the Gauss-Newton optimization algorithm. In this study, we have successfully integrated a lossy compression algorithm into our F3DT SI workflow to significantly reduce the disk space for storing the strain tensor fields. The compressor supports a user-specified tolerance for bounding the error, and can be integrated into our finite-difference wave-propagation simulation code used for computing the strain fields. The decompressor can be integrated into the kernel calculation code that reads the strain fields from the disk and compute the data sensitivity kernels. During the wave-propagation simulations, we compress the strain fields before writing them to the disk. To compute the data sensitivity kernels, we read the compressed strain fields from the disk and decompress themmore »
- Authors:
-
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Univ. of Wyoming, Laramie, WY (United States)
- National Cheng Kung Univ., Tainan City (Taiwan)
- Publication Date:
- Research Org.:
- Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR)
- OSTI Identifier:
- 1258513
- Alternate Identifier(s):
- OSTI ID: 1358859
- Report Number(s):
- LLNL-JRNL-678678
Journal ID: ISSN 0098-3004
- Grant/Contract Number:
- AC52-07NA27344
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Computers and Geosciences
- Additional Journal Information:
- Journal Volume: 93; Journal ID: ISSN 0098-3004
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING
Citation Formats
Lindstrom, Peter, Chen, Po, and Lee, En-Jui. Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression. United States: N. p., 2016.
Web. doi:10.1016/j.cageo.2016.04.009.
Lindstrom, Peter, Chen, Po, & Lee, En-Jui. Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression. United States. https://doi.org/10.1016/j.cageo.2016.04.009
Lindstrom, Peter, Chen, Po, and Lee, En-Jui. Thu .
"Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression". United States. https://doi.org/10.1016/j.cageo.2016.04.009. https://www.osti.gov/servlets/purl/1258513.
@article{osti_1258513,
title = {Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression},
author = {Lindstrom, Peter and Chen, Po and Lee, En-Jui},
abstractNote = {Full-3D seismic waveform tomography (F3DT) is the latest seismic tomography technique that can assimilate broadband, multi-component seismic waveform observations into high-resolution 3D subsurface seismic structure models. The main drawback in the current F3DT implementation, in particular the scattering-integral implementation (F3DT-SI), is the high disk storage cost and the associated I/O overhead of archiving the 4D space-time wavefields of the receiver- or source-side strain tensors. The strain tensor fields are needed for computing the data sensitivity kernels, which are used for constructing the Jacobian matrix in the Gauss-Newton optimization algorithm. In this study, we have successfully integrated a lossy compression algorithm into our F3DT SI workflow to significantly reduce the disk space for storing the strain tensor fields. The compressor supports a user-specified tolerance for bounding the error, and can be integrated into our finite-difference wave-propagation simulation code used for computing the strain fields. The decompressor can be integrated into the kernel calculation code that reads the strain fields from the disk and compute the data sensitivity kernels. During the wave-propagation simulations, we compress the strain fields before writing them to the disk. To compute the data sensitivity kernels, we read the compressed strain fields from the disk and decompress them before using them in kernel calculations. Experiments using a realistic dataset in our California statewide F3DT project have shown that we can reduce the strain-field disk storage by at least an order of magnitude with acceptable loss, and also improve the overall I/O performance of the entire F3DT-SI workflow significantly. The integration of the lossy online compressor may potentially open up the possibilities of the wide adoption of F3DT-SI in routine seismic tomography practices in the near future.},
doi = {10.1016/j.cageo.2016.04.009},
journal = {Computers and Geosciences},
number = ,
volume = 93,
place = {United States},
year = {Thu May 05 00:00:00 EDT 2016},
month = {Thu May 05 00:00:00 EDT 2016}
}
Web of Science
Works referenced in this record:
Seismic Tomography
journal, October 1984
- Anderson, Don L.; Dziewonski, Adam M.
- Scientific American, Vol. 251, Issue 4
Adaptive Local Cosine transform for Seismic Image Compression
conference, December 2006
- Aparna, P.; David, Sumam
- 2006 International Conference on Advanced Computing and Communications
Bandwidth Reduction through Multithreaded Compression of Seismic Images
conference, May 2011
- Aqrawi, Ahmed A.; Elster, Anne C.
- Distributed Processing, Workshops and Phd Forum (IPDPSW), 2011 IEEE International Symposium on Parallel and Distributed Processing Workshops and Phd Forum
Low bit-rate efficient compression for seismic data
journal, December 2001
- Averbuch, A. Z.; Meyer, F.; Stromberg, J. -O.
- IEEE Transactions on Image Processing, Vol. 10, Issue 12
Finite-Moment Tensor of the 3 September 2002 Yorba Linda Earthquake
journal, June 2005
- Chen, P.
- Bulletin of the Seismological Society of America, Vol. 95, Issue 3
Full 3D Tomography for the Crustal Structure of the Los Angeles Region
journal, August 2007
- Chen, P.; Zhao, L.; Jordan, T. H.
- Bulletin of the Seismological Society of America, Vol. 97, Issue 4
Full three-dimensional tomography: a comparison between the scattering-integral and adjoint-wavefield methods
journal, July 2007
- Chen, Po; Jordan, Thomas H.; Zhao, Li
- Geophysical Journal International, Vol. 170, Issue 1
The TeraShake Computational Platform for Large-Scale Earthquake Simulations
book, January 2009
- Cui, Yifeng; Olsen, Kim; Chourasia, Amit
- Advances in Geocomputing
Seismic Data Compression Using 2D Lifting-Wavelet Algorithms
journal, January 2015
- Fajardo, Carlos; Reyes, Oscar Mauricio; Ramirez, Ana
- Ingeniería y Ciencia, Vol. 11, Issue 21
Full seismic waveform tomography for upper-mantle structure in the Australasian region using adjoint methods
journal, December 2009
- Fichtner, Andreas; Kennett, Brian L. N.; Igel, Heiner
- Geophysical Journal International, Vol. 179, Issue 3, p. 1703-1725
Group testing for image compression
journal, August 2002
- Hong, E. S.; Ladner, R. E.
- IEEE Transactions on Image Processing, Vol. 11, Issue 8
Seismic data compression
conference, January 1990
- Jonsson, S. B.; Spanias, A. S.
- Ninth Annual International Phoenix Conference on Computers and Communications. 1990 Conference Proceedings
Full-3-D tomography for crustal structure in Southern California based on the scattering-integral and the adjoint-wavefield methods: F3DT FOR SOUTHERN CALIFORNIA CRUST
journal, August 2014
- Lee, En-Jui; Chen, Po; Jordan, Thomas H.
- Journal of Geophysical Research: Solid Earth, Vol. 119, Issue 8
Testing Waveform Predictions of 3D Velocity Models against Two Recent Los Angeles Earthquakes
journal, October 2014
- Lee, E. -J.; Chen, P.; Jordan, T. H.
- Seismological Research Letters, Vol. 85, Issue 6
An optimized parallel LSQR algorithm for seismic tomography
journal, December 2013
- Lee, En-Jui; Huang, He; Dennis, John M.
- Computers & Geosciences, Vol. 61
Fast and Efficient Compression of Floating-Point Data
journal, September 2006
- Lindstrom, Peter; Isenburg, Martin
- IEEE Transactions on Visualization and Computer Graphics, Vol. 12, Issue 5
Fixed-Rate Compressed Floating-Point Arrays
journal, December 2014
- Lindstrom, Peter
- IEEE Transactions on Visualization and Computer Graphics, Vol. 20, Issue 12
i PACKMAN: high-quality, low-complexity texture compression for mobile phones
conference, January 2005
- Ström, Jacob; Akenine-Möller, Tomas
- Proceedings of the ACM SIGGRAPH/EUROGRAPHICS conference on Graphics hardware - HWWS '05
Seismic tomography of the southern California crust based on spectral-element and adjoint methods
journal, January 2010
- Tape, Carl; Liu, Qinya; Maggi, Alessia
- Geophysical Journal International, Vol. 180, Issue 1
Comparison of wavelet image coding schemes for seismic data compression
conference, October 1997
- Vassiliou, Anthony A.; Wickerhouser, Mladen V.
- Optical Science, Engineering and Instrumentation '97, SPIE Proceedings
Seismic data compression based on integer wavelet transform
journal, November 2004
- Wang, Xi-zhen; Teng, Yun-tian; Gao, Meng-tan
- Acta Seismologica Sinica, Vol. 17, Issue S1
Seismic data compression by an adaptive local cosine/sine transform and its effects on migration: Adaptive local cosine/sine transform
journal, November 2000
- Wang, Yongzhong; Wu, Ru-Shan
- Geophysical Prospecting, Vol. 48, Issue 6
Seismic data Compression Methods
journal, August 1974
- Wood, Lawrence C.
- GEOPHYSICS, Vol. 39, Issue 4
Adaptive Seismic Data Compression Using Wavelet Packets
conference, July 2006
- Wu, W.; Yang, Z.; Qin, Q.
- 2006 IEEE International Symposium on Geoscience and Remote Sensing
Fast Lossless Compression of Seismic Floating-Point Data
conference, May 2009
- Xie, Xing; Qin, Qianqing
- 2009 International Forum on Information Technology and Applications (IFITA)
Frechet Kernels for Imaging Regional Earth Structure Based on Three-Dimensional Reference Models
journal, December 2005
- Zhao, L.
- Bulletin of the Seismological Society of America, Vol. 95, Issue 6
Strain Green's Tensors, Reciprocity, and Their Applications to Seismic Source and Structure Studies
journal, October 2006
- Zhao, L.; Chen, P.; Jordan, T. H.
- Bulletin of the Seismological Society of America, Vol. 96, Issue 5
A fast compression algorithm for seismic data from non-cable seismographs
conference, October 2012
- Zheng, Fan; Liu, Shufen
- 2012 World Congress on Information and Communication Technologies (WICT)
Works referencing / citing this record:
Data Reduction Techniques for Simulation, Visualization and Data Analysis: Survey on Scientific Data Reduction Techniques
journal, March 2018
- Li, S.; Marsaglia, N.; Garth, C.
- Computer Graphics Forum, Vol. 37, Issue 6
Compression Challenges in Large Scale Partial Differential Equation Solvers
journal, September 2019
- Götschel, Sebastian; Weiser, Martin
- Algorithms, Vol. 12, Issue 9
An Efficient Encoding Algorithm Using Local Path on Huffman Encoding Algorithm for Compression
journal, February 2019
- Erdal, Erdal; Ergüzen, Atilla
- Applied Sciences, Vol. 9, Issue 4