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

Title: W15_geowave “3D stochastic full waveform geophysical models”

Abstract

This report describes results for simulated, full waveform geophysical models.

Authors:
 [1];  [1]
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1342871
Report Number(s):
LA-UR-17-20849
DOE Contract Number:
AC52-06NA25396
Resource Type:
Technical Report
Country of Publication:
United States
Language:
English
Subject:
58 GEOSCIENCES; Planetary Sciences

Citation Formats

Larmat, Carene, and Maceira, Monica. W15_geowave “3D stochastic full waveform geophysical models”. United States: N. p., 2017. Web. doi:10.2172/1342871.
Larmat, Carene, & Maceira, Monica. W15_geowave “3D stochastic full waveform geophysical models”. United States. doi:10.2172/1342871.
Larmat, Carene, and Maceira, Monica. Fri . "W15_geowave “3D stochastic full waveform geophysical models”". United States. doi:10.2172/1342871. https://www.osti.gov/servlets/purl/1342871.
@article{osti_1342871,
title = {W15_geowave “3D stochastic full waveform geophysical models”},
author = {Larmat, Carene and Maceira, Monica},
abstractNote = {This report describes results for simulated, full waveform geophysical models.},
doi = {10.2172/1342871},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Fri Feb 03 00:00:00 EST 2017},
month = {Fri Feb 03 00:00:00 EST 2017}
}

Technical Report:

Save / Share:
  • Performance of the MCMC inversion according to the number of cores for the computation. A) 64 cores. B) 480 cores. C) 816 cores. The true model is represented by the black line. Vsv is the wave speed of S waves polarized in the vertical plane, ξ is an anisotropy parameter. The Earth is highly anisotropics; the wavespeed of seismic waves depends on the polarization of the wave. Seismic inversion of the elastic structure is usually limited to isotropic information such as Vsv. Our research looked at the inversion of Earth anisotropy.
  • Adequate path calibrations are crucial for improving the accuracy of seismic event location and origin time, size, and mechanism, as required for CTBT monitoring. There is considerable information on structure in broadband seismograms that is currently not fully utilized. The limitations have been largely theoretical. the development and application to solid earth problems of powerful numerical techniques, such as the Spectral Element Method (SEM), has opened a new era, and theoretically, it should be possible to compute the complete predicted wavefield accurately without any restrictions on the strength or spatial extent of heterogeneity. This approach requires considerable computational power, whichmore » is currently not fully reachable in practice. We propose an approach which relies on a cascade of increasingly accurate theoretical approximations for the computation of the seismic wavefield to develop a model of regional structure for the area of Eurasia located between longitudes of 30 and 150 degrees E, and latitudes of -10 to 60 degrees North. The selected area is particularly suitable for the purpose of this experiment, as it is highly heterogeneous, presenting a challenge for calibration purposes, but it is well surrounded by earthquake sources and, even though they are sparsely distributed, a significant number of high quality broadband digital stations exist, for which data are readily accessible through IRIS (Incorporated Research Institutions for Seismology) and the FDSN (Federation of Digital Seismic Networks). The starting models used will be a combination of a-priori 3D models recently developed for this region, combining various geophysical and seismological data, and a major goal of this study will be to refine these models so as to fit a variety of seismic waveforms and phases.« less
  • Adequate path calibrations are crucial for improving the accuracy of seismic event location and origin time, size, and mechanism, as required for CTBT monitoring. There is considerable information on structure in broadband seismograms that is currently not fully utilized. The limitations have been largely theoretical. the development and application to solid earth problems of powerful numerical techniques, such as the Spectral Element Method (SEM), has opened a new era, and theoretically, it should be possible to compute the complete predicted wavefield accurately without any restrictions on the strength or spatial extent of heterogeneity. This approach requires considerable computational power, whichmore » is currently not fully reachable in practice. We propose an approach which relies on a cascade of increasingly accurate theoretical approximations for the computation of the seismic wavefield to develop a model of regional structure for the area of Eurasia located between longitudes of 30 and 150 degrees E, and latitudes of -10 to 60 degrees North. The selected area is particularly suitable for the purpose of this experiment, as it is highly heterogeneous, presenting a challenge for calibration purposes, but it is well surrounded by earthquake sources and, even though they are sparsely distributed, a significant number of high quality broadband digital stations exist, for which data are readily accessible through IRIS (Incorporated Research Institutions for Seismology) and the FDSN (Federation of Digital Seismic Networks). The starting models used will be a combination of a-priori 3D models recently developed for this region, combining various geophysical and seismological data, and a major goal of this study will be to refine these models so as to fit a variety of seismic waveforms and phases.« less
  • The goal of this project was to develop a method to exploit viscoelastic rock and fluid properties to greatly enhance the sensitivity of surface seismic measurements to the presence of hydrocarbon saturation. To reach the objective, Southwest Research Institute scientists used well log, lithology, production, and 3D seismic data from an oil reservoir located on the Waggoner Ranch in north central Texas. The project was organized in three phases. In the first phase, we applied modeling techniques to investigate seismic- and acoustic-frequency wave attenuation and its effect on observable wave attributes. We also gathered existing data and acquired new datamore » from the Waggoner Ranch field, so that all needed information was in place for the second phase. During the second phase, we developed methods to extract attenuation from borehole acoustic and surface seismic data. These methods were tested on synthetic data constructed from realistic models and real data. In the third and final phase of the project, we applied this technology to a full data set from the Waggoner site. The results presented in this Final Report show that geological conditions at the site did not allow us to obtain interpretable results from the Q processing algorithm for 3D seismic data. However, the Q-log processing algorithm was successfully applied to full waveform sonic data from the Waggoner site. A significant part of this project was technology transfer. We have published several papers and conducted presentations at professional conferences. In particular, we presented the Q-log algorithm and applications at the Society of Exploration Geophysicists (SEG) Development and Production Forum in Austin, Texas, in May 2005. The presentation attracted significant interest from the attendees and, at the request of the SEG delegates, it was placed on the Southwest Research Institute Internet site. The presentation can be obtained from the following link: http://www.swri.org/4org/d15/elecsys/resgeo/ppt/Algorithm.pps In addition, we presented a second application of the Q algorithm at the SEG International Conference in Houston, Texas, in May 2005. The presentation attracted significant interest there as well, and it can be obtained from the following link: http://www.swri.org/4org/d15/elecsys/resgeo/ppt/attenuation.pps.« less
  • We present a detailed study of a novel Cherenkov imaging detector called the Focusing DIRC (FDIRC) with waveform digitizing electronics. In this test study, the FDIRC prototype has been instrumented with seven Hamamatsu H-8500 MaPMTs. Waveforms from ~450 pixels are digitized with waveform sampling electronics based on the BLAB2 ASIC, operating at a sampling speed of ~2.5 GSa/s. The FDIRC prototype was tested in a large cosmic ray telescope (CRT) providing 3D muon tracks with ~1.5 mrad angular resolution and muon energy of E muon greater than 1.6 GeV. In this study we provide a detailed analysis of the tailsmore » in the Cherenkov angle distribution as a function of various variables, compare experimental results with simulation, and identify the major contributions to the tails. We demonstrate that to see the full impact of these tails on the Cherenkov angle resolution, it is crucial to use 3D tracks, and have a full understanding of the role of ambiguities. These issues could not be fully explored in previous FDIRC studies where the beam was perpendicular to the quartz radiator bars. This work is relevant for the final FDIRC prototype of the PID detector at SuperB, which will be tested this year in the CRT setup.« less