skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Subsurface fault geometries in Southern California illuminated through Full-3D Seismic Waveform Tomography (F3DT)

ORCiD logo;
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Additional Journal Information:
Journal Volume: 703-704; Journal Issue: C; Related Information: CHORUS Timestamp: 2017-10-08 22:10:24; Journal ID: ISSN 0040-1951
Country of Publication:

Citation Formats

Lee, En-Jui, and Chen, Po. Subsurface fault geometries in Southern California illuminated through Full-3D Seismic Waveform Tomography (F3DT). Netherlands: N. p., 2017. Web. doi:10.1016/j.tecto.2017.03.005.
Lee, En-Jui, & Chen, Po. Subsurface fault geometries in Southern California illuminated through Full-3D Seismic Waveform Tomography (F3DT). Netherlands. doi:10.1016/j.tecto.2017.03.005.
Lee, En-Jui, and Chen, Po. Sat . "Subsurface fault geometries in Southern California illuminated through Full-3D Seismic Waveform Tomography (F3DT)". Netherlands. doi:10.1016/j.tecto.2017.03.005.
title = {Subsurface fault geometries in Southern California illuminated through Full-3D Seismic Waveform Tomography (F3DT)},
author = {Lee, En-Jui and Chen, Po},
abstractNote = {},
doi = {10.1016/j.tecto.2017.03.005},
journal = {Tectonophysics},
number = C,
volume = 703-704,
place = {Netherlands},
year = {Sat Apr 01 00:00:00 EDT 2017},
month = {Sat Apr 01 00:00:00 EDT 2017}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.tecto.2017.03.005

Save / Share:
  • 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 algorithmmore » 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.« less
  • Cited by 3
  • Seismic tomography becomes important tool recently for imaging complex subsurface. It is well known that imaging complex rich fault zone is difficult. In this paper, The application of time domain inverse scattering wave tomography to image the complex fault zone would be shown on this paper, especially an efficient time domain inverse scattering tomography and their run in cluster parallel computer which has been developed. This algorithm is purely based on scattering theory through solving Lippmann Schwienger integral by using Born's approximation. In this paper, it is shown the robustness of this algorithm especially in avoiding the inversion trapped inmore » local minimum to reach global minimum. A large data are solved by windowing and blocking technique of memory as well as computation. Parameter of windowing computation is based on shot gather's aperture. This windowing technique reduces memory as well as computation significantly. This parallel algorithm is done by means cluster system of 120 processors from 20 nodes of AMD Phenom II. Benchmarking of this algorithm is done by means Marmoussi model which can be representative of complex rich fault area. It is shown that the proposed method can image clearly the rich fault and complex zone in Marmoussi model even though the initial model is quite far from the true model. Therefore, this method can be as one of solution to image the very complex mode.« less
  • High-frequency, cross-well seismic data, from the Midale oil field of southeastern Saskatchewan, are analyzed for direct and reflected energy. The goal of the analysis is to produce interpretable sections to assist in enhanced oil recovery activities (CO[sub 2] injection) in this field. Direct arrivals are used for velocity information while reflected arrivals are used for velocity information while reflected arrivals are processed into a reflection image. Raw field data show a complex assortment of wave types that includes direct compressional and shear waves and reflected shear waves. A traveltime inversion technique (layer stripping via ray tracing) is used to obtainmore » P- and S-wave interval velocities from the respective direct arrivals. The velocities from the cross-well inversion and the sonic log are in reasonable agreement. The subsurface coverage of the cross-well geometry is investigated; it covers zones extending from the source well to the receiver well and includes regions above and below the source/receiver depths. Upgoing and downgoing primary reflections are processed, in a manner similar to the vertical seismic profiling/common-depth-point (VSP/CDP) map, to construct the cross-well images. A final section is produced by summing the individual reflection images from each receiver-gather map. This section provides an image with evidence of strata thicknesses down to about 1 m. Synthetic seismograms are used to interpret the final sections. Correlations can be drawn between some of the events on the synthetic seismograms and the cross-well image.« less
  • A rigorous full waveform inversion of seismic data has been a challenging subject partly because of the lack of precise knowledge of the source. Since currently available approaches involve some form of approximations to the source, inversion results are subject to the quality and the choice of the source information used. We propose a new full waveform inversion methodology that does not involve source spectrum information. Thus potential inversion errors due to source estimation can be eliminated. A gather of seismic traces is first Fourier-transformed into the frequency domain and a normalized wavefield is obtained for each trace in themore » frequency domain. Normalization is done with respect to the frequency response of a reference trace selected from the gather, so the complex-valued normalized wavefield is dimensionless. The source spectrum is eliminated during the normalization procedure. With its source spectrum eliminated, the normalized wavefield allows us construction of an inversion algorithm without the source information. The inversion algorithm minimizes misfits between measured normalized wavefield and numerically computed normalized wavefield. The proposed approach has been successfully demonstrated using a simple two-dimensional scalar problem.« less