Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SEISMICITY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Fault Detection on Seismic Structural Images Using a Nested Residual U-Net

Automatic identification of faults on seismic structural images is a challenging yet crucial task in quantitative seismic interpretation. Human picking or attribute-based fault detection methods may misidentify faults on noisy, complex seismic images. In this work, we develop a new automatic fault detection method using a nested residual U-shaped convolutional neural network. Each of the encoders and decoders in this neural network is a residual U-Net, leading to a nested architecture. The final fault map results from the fusion of three fault maps with low, medium, and high fault resolutions. We demonstrate the excellent fault-detection capability of our nested neural network using a series of synthetic and field seismic images. We find that our approach produces clearer and more interpretable fault maps than the current state-of-the-art U-Net fault detection method, particularly on noisy seismic images. Our new automatic fault detection method can facilitate reliable quantitative seismic interpretation on field seismic images.

58 GEOSCIENCES↗
DOE OSTI2023

Coordinate-Based Seismic Interpolation in Irregular Land Survey: A Deep Internal Learning Approach

Physical and budget constraints often result in irregular sampling, which complicates accurate subsurface imaging. Preprocessing approaches, such as missing trace or shot interpolation, are typically employed to enhance seismic data in such cases. Recently, deep learning has been used to address the trace interpolation problem at the expense of large amounts of training data to adequately represent typical seismic events. Nonetheless, most research in this area has focused on trace reconstruction, with little attention having been devoted to shot interpolation. Furthermore, existing methods assume regularly spaced receivers/sources failing in approximating seismic data from real (irregular) surveys. This work presents a novel shot gather interpolation approach which uses a continuous coordinate-based representation of the acquired seismic wavefield parameterized by a neural network. The proposed unsupervised approach, which we call coordinate-based seismic interpolation (CoBSI), enables the prediction of specific seismic characteristics in irregular land surveys without using external data during neural network training. Importantly, experimental results on real and synthetic 3-D data validate the ability of the proposed method to estimate continuous smooth seismic events in the time-space and frequency-wavenumber domains, improving sparsity or low-rank-based interpolation methods.

58 GEOSCIENCES↗
DOE OSTI2022

High-Precision Characterization of Seismicity from the 2022 Hunga Tonga-Hunga Ha'apai Volcanic Eruption

The earthquake swarm accompanying the January 2022 Hunga Tonga-Hunga Ha'apai (HTHH) volcanic eruption includes a large number of posteruptive moderate-magnitude seismic events and presents a unique opportunity to use remote monitoring methods to characterize and compare seismic activity with other historical caldera-forming eruptions. We compute improved epicentroid locations, magnitudes, and regional moment tensors of seismic events from this earthquake swarm using regional to teleseismic surface-wave cross correlation and waveform modeling. Precise relative locations of 91 seismic events derived from 59,047 intermediate-period Rayleigh- and Love-wave cross-correlation measurements collapse into a small area surrounding the volcano and exhibit a southeastern time-dependent migration. Regional moment tensors and observed waveforms indicate that these events have a similar mechanism and exhibit a strong positive compensated linear vector dipole component. Precise relative magnitudes agree with regional moment tensor moment magnitude ($M_w$) estimates while also showing that event sizes and frequency increase during the days after the eruption followed by a period of several weeks of less frequent seismicity of a similar size. The combined information from visual observation and early geologic models indicate that the observed seismicity may be the result of a complex series of events that occurred after the explosive eruption on 15 January, possibly involving rapid resupply of the magma chamber shortly after the eruption and additional faulting and instability in the following weeks. In addition, we identify and characterize an $M_w$ 4.5 event five days before the paroxysmal explosion on 15 January, indicating that additional seismic events preceding the main eruption could have been identified with improved local monitoring. As a result, our analysis of the HTHH eruption sequence demonstrates the value of potentially utilizing teleseismic surface-wave cross correlation and waveform modeling methods to assist in the detailed analysis of remote volcanic eruption sequences.

58 GEOSCIENCES↗
DOE OSTI2012

Seismic Imaging and Inversion: Application of Linear Theory (2012), Cambridge University Press, co-authored with Bob Stolt

Extracting information from seismic data requires knowledge of seismic wave propagation and reflection. The commonly used method involves solving linearly for a reflectivity at every point within the Earth, but this book follows an alternative approach which invokes inverse scattering theory. By developing the theory of seismic imaging from basic principles, the authors relate the different models of seismic propagation, reflection and imaging - thus providing links to reflectivity-based imaging on the one hand and to nonlinear seismic inversion on the other. The comprehensive and physically complete linear imaging foundation developed presents new results at the leading edge of seismic processing for target location and identification. This book serves as a fundamental guide to seismic imaging principles and algorithms and their foundation in inverse scattering theory and is a valuable resource for working geoscientists, scientific programmers and theoretical physicists.

58 GEOSCIENCES↗
DOE OSTI2019

Work Plan for the Site-wide SSHAC Level 3 Probabilistic Seismic Hazard Analysis at INL

This work plan provides the overall framework, key activities, key tasks, and schedule for the execution of a probabilistic seismic hazard analysis at Idaho National Laboratory (INL). The probabilistic seismic hazard analysis will follow processes and procedures appropriate for the Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 study provided in U. S. Nuclear Regulatory Commission guidance, NUREG-2213. The SSHAC Level 3 probabilistic seismic hazard analysis (PSHA) is being conducted as an INL site-wide effort. The work plan was developed and revised by participants who will perform the SSHAC Level 3 PSHA. The work plan also includes the quality assurance requirements commensurate with American Society of Mechanical Engineers Nuclear Quality Assurance NQA 1 2008/2009 and 2011 Addendum for design nuclear facilities under the U.S. Department of Energy and U. S. Nuclear Regulatory Commission. The scope of the SSHAC Level 3 PSHA is designed to produce the mean-centered definition of the seismic hazard with the appropriate treatment of uncertainties supported by the associated technical justifications. The SSHAC Level 3 study will properly implement and completely document the activities of evaluation, integration, and documentation. The evaluation process will consider data, models, and methods proposed by the larger technical community that are relevant to the PSHA at INL. The integration process will build seismic-source and ground-motion characterization models that represent the center, body, and range of technically defensible interpretations, taking into consideration the outcomes of the evaluation process. Key tasks and activities—such as a project kickoff meeting, three workshops, four formal working meetings, and a briefing meeting with the Participatory Peer Review Panel—will occur for the seven-member technical integration teams to construct seismic-source and ground-motion-characterization models. Both preliminary and final hazard-input documents will be developed to perform sensitivity analyses and the final PSHA, respectively. The entire study will be thoroughly documented in draft and final reports. The Participatory Peer Review Panel will review the SSHAC Level 3 study from start to finish, including the draft and final INL Site-wide PSHA reports, so they can provide a closure letter documenting the study’s adherence to SSHAC processes and technical adequacy of the final products.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2011

Seismic Fragility Analysis of a Degraded Condensate Storage Tank

The Korea Atomic Energy Research Institute (KAERI) and Brookhaven National Laboratory are conducting a collaborative research project to develop seismic capability evaluation technology for degraded structures and components in nuclear power plants (NPPs). One of the goals of this collaboration endeavor is to develop seismic fragility analysis methods that consider the potential effects of age-related degradation of structures, systems, and components (SSCs). The essential part of this collaboration is aimed at achieving a better understanding of the effects of aging on the performance of SSCs and ultimately on the safety of NPPs. A recent search of the degradation occurrences of structures and passive components (SPCs) showed that the rate of aging related degradation in NPPs was not significantly large but increasing, as the plants get older. The slow but increasing rate of degradation of SPCs can potentially affect the safety of the older plants and become an important factor in decision making in the current trend of extending the operating license period of the plants (e.g., in the U.S. from 40 years to 60 years, and even potentially to 80 years). The condition and performance of major aged NPP structures such as the containment contributes to the life span of a plant. A frequent misconception of such low degradation rate of SPCs is that such degradation may not pose significant risk to plant safety. However, under low probability high consequence initiating events, such as large earthquakes, SPCs that have slowly degraded over many years could potentially affect plant safety and these effects need to be better understood. As part of the KAERI-BNL collaboration, a condensate storage tank (CST) was analyzed to estimate its seismic fragility capacities under various postulated degradation scenarios. CSTs were shown to have a significant impact on the seismic core damage frequency of a nuclear power plant. The seismic fragility capacity of the CST was developed for five cases: (1) a baseline analysis where the design condition (undegraded) is assumed, (2) a scenario with degraded stainless steel tank shell, (3) a scenario with degraded anchor bolts, (4) a scenario with anchorage concrete cracking, and (5) a perfect correlation of the above three degradation scenarios. This paper will present the methodology for the time-dependent fragility calculation and discuss the insights drawn from this study. To achieve a better understanding of the effects of aging on the performance of structures and passive components (SPCs) in nuclear power plants (NPPs), the Korea Atomic Energy Research Institute (KAERI) and Brookhaven National Laboratory (BNL) are collaborating to develop seismic fragility analysis methods that consider age-related degradation of SPCs. The rate of age-related degradation of SPCs was not found to be significantly large, but increasing as the plants get older. The slow but increasing rate of degradation of SPCs can potentially affect the safety of the older plants and become an important factor in decision making in the current trend of extending the operating license period of the plants (e.g., in the U.S. from 40 years to 60 years, and even potentially to 80 years). In this paper, a condensate storage tank (CST) was analyzed to estimate its seismic fragility capacities under various postulated degradation scenarios. This paper will present the methodology for the time-dependent fragility calculation and discuss the insights drawn from this study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Factors Affecting Seismic Velocity in Alluvium

Yucca Flat at the Nevada National Security Site has been selected as the Source Physics Experiment (SPE) Dry Alluvium Geology Phase II site. The alluvium in this part of Yucca Flat is typical of desert basin fill, with discontinuous beds that are highly variable in clast size and provenance. Detailed understanding of the subsurface geology will be needed for interpretation of the SPE seismic data. A 3D seismic velocity model, created for Yucca Flat using interval seismic velocity data, shows variations in velocity within alluvium near the SPE Phase II site beyond the usual gradual increase of density with depth due to compaction. In this study we examined borehole lithologic logs, geophysical logs, downhole videos, and laboratory analyses of sidewall core samples to understand which characteristics of the alluvium are related to these variations in seismic velocity. Seismic velocity of alluvium is generally related to its density, which can be affected by sediment provenance, clast size, gravel percentage, and matrix properties, in addition to compaction. This study presents a preliminary subdivision of the alluvial strata in the SPE Phase II area into mappable units expected to be significant to seismic modeling. Further refinements of the alluvial units may be possible when seismic data are obtained from SPE Phase II tests.

58 GEOSCIENCES↗
DOE OSTI2016

Fault structure, stress, or pressure control of the seismicity in shale? Insights from a controlled experiment of fluid-induced fault reactivation

It is clear that clay formations are present in reservoirs and earthquake faults, but questions remain on their mechanical behavior, as they can vary from ductile (aseismic) to brittle (seismic). An experiment, at a scale of 10 m, aims to reactivate a natural fault by fluid pressure in shale materials. The injection area was surrounded by a dense monitoring network comprising pressure, deformation, and seismicity sensors, in a well-characterized geological setting. Thirty-two microseismic events were recorded during several injection phases in five different locations within the fault zone. Their computed magnitude ranged between -4.3 and -3.7. Their spatiotemporal distribution, compared with the measured displacement at the injection points, shows that most of the deformation induced by the injection is aseismic. Whether the seismicity is controlled by the fault architecture, mineralogy of fracture filling, fluid, and/or stress state is then discussed. The fault damage zone architecture and mineralogy are of crucial importance, as seismic slip mainly localizes on the sealed-with-calcite fractures which predominate in the fault damage zone. As no seismicity is observed in the close vicinity of the injection areas, the presence of fluid seems to prevent seismic slips. The fault core acts as an impermeable hydraulic barrier that favors fluid confinement and pressurization. Therefore, the seismic behavior seems to be strongly sensitive to the structural heterogeneity (including permeability) of the fault zone, which leads to a heterogeneous stress response to the pressurized volume.

58 GEOSCIENCES↗
DOE OSTI2025

Toward more-robust, AI-enabled subsurface seismic imaging for geotechnical applications

Non-invasive seismic imaging has the potential to cost-effectively evaluate large volumes of subsurface material to inform geotechnical site investigation. However, seismic imaging using full waveform inversion (FWI) requires significant computational time and is dependent on an initial starting model. As a result, FWI has not yet been widely adopted into geotechnical practice. Previous efforts, on relatively simple two-layered models, indicate that data-driven artificial intelligence (AI) models may be as effective as FWI at predicting 2D images of shear wave velocity (V s ). Furthermore, the AI model predictions can be made almost instantaneously after data acquisition and do not require an initial starting model. We examine the generality of these findings by developing a new AI model for subsurface seismic imaging, whereby we make several notable contributions. First, we architect a multimodal AI model that combines time- and frequency-domain representations of the seismic wavefield to predict a 50 m by 20 m subsurface image of V s . Second, we developed a new diverse dataset of 100,000 images with their corresponding seismic wavefields to train the AI model. Third, we propose four physics-informed data augmentations for data-driven seismic imaging. Fourth, we develop two prediction consistency tests to evaluate the model’s performance when the true subsurface is unknown. Our final model, which has been made publicly available, is capable of predicting a subsurface V s image from a single seismic wavefield with an average, mean absolute percent error (MAPE) of 24 %. The predictive model is applied to a field dataset and shown to be consistent with local geology and shear-wave refraction measurements from the same location.

Artificial intelligence↗
DOE OSTI2015

Induced seismicity in geothermal reservoirs: A review of forecasting approaches

In order to reach Europe's 2020 and 2050 targets in terms of greenhouse gas emissions, geothermal resources will have to contribute substantially to meeting carbon-free energy needs. However, public opinion may prevent future large-scale application of deep geothermal power plants, because induced seismicity is often perceived as an unsolicited and uncontrollable side effect of geothermal development. In the last decade, significant advances were made in the development of models to forecast induced seismicity, which are either based on catalogues of induced seismicity, on the underlying physical processes, or on a hybrid philosophy. In this paper, we provide a comprehensive overview of the existing approaches applied to geothermal contexts. This overview will outline the advantages and drawbacks of the different approaches, identify the gaps in our understanding, and describe the needs for geothermal observations. Most of the forecasting approaches focus on the stimulation phase of enhanced geothermal systems which are most prone to generate seismic events. Besides the statistical models suited for real-time applications during reservoir stimulation, the physics-based models have the advantage of considering sub-surface characteristics and estimating the impact of fluid circulation on the reservoir. Thus, to mitigate induced seismicity during major hydraulic stimulations, application of hybrid methods in a decision support system seems the best available solution. So far, however, little attention has been paid to geochemical effects on the failure process and to production periods. Quantitative modelling of induced seismicity still is a challenging and complex matter. Appropriate resources remain to be invested for the scientific community to continue its research and development efforts to successfully forecast induced seismicity in geothermal fields. This is a prerequisite for making this renewable energy resource sustainable and accessible worldwide.

58 GEOSCIENCES↗
DOE OSTI2018

Seismicity induced by simultaneous abrupt changes of injection rate and well pressure in Hutubi Gas Field

Hutubi gas field, the largest gas storage field in China, has been operated on annual injection/extraction cycles since 9 June 2013. Here, we study the seismicity near the gas field from 9 June 2013 to 22 October 2015, a time span that the gas field has experienced three injection periods and two extraction periods, and explore its physical mechanism based on the relationship between seismicity and field operation. We identify 273 events (ML >1) in the region within 10 km of the gas field, with 97% of those occurring in the first two injection periods, 0.4% in the third injection period and 1% in the two extraction periods. Seismicity in the first two injection periods occurs mostly as shallow clusters (focal depth < 2 km) at two locations: with one along the fault that marks the southern boundary of the gas field and the other about 2 km south to the southeastern tip of the gas field with the seismicity distributed along north-south direction. The seismicity does not correlate with total gas injection volume, injection rate or well pressure. It instead occurs 11–17 hours after simultaneous abrupt increases/decreases of gas injection rate and well pressure in the field operation in the first two injection periods when some accumulative injection has reached. Such relationship is consistent with a physical mechanism that the seismicity near Hutubi gas field is induced on pore-pressured faults with a rate- and state-dependent friction law through an abrupt change of stress in elastic and undrained poroelastic responses to simultaneous abrupt changes of injection rate and well pressure. Our study also points to the possibility that induced seismicity may be controllable in some practical field operations.

58 GEOSCIENCES↗
DOE OSTI2020

A New Crater Near InSight: Implications for Seismic Impact Detectability on Mars

A new 1.5 m diameter impact crater was discovered on Mars only ~40 km from the InSight lander. Context camera images constrained its formation between 21 February and 6 April 2019; follow-up High Resolution Imaging Science Experiment images resolved the crater. During this time period, three seismic events were identified in InSight data. In this paper, we derive expected seismic signal characteristics and use them to evaluate each of the seismic events. However, none of them can definitively be associated with this source. Atmospheric perturbations are generally expected to be generated during impacts; however, in this case, no signal could be identified as related to the known impact. Using scaling relationships based on the terrestrial and lunar analogs and numerical modeling, we predict the amplitude, peak frequency, and duration of the seismic signal that would have emanated from this impact. The predicted amplitude falls near the lowest levels of the measured seismometer noise for the predicted frequency. Hence, it is not surprising this impact event was not positively identified in the seismic data. Finding this crater was a lucky event as its formation this close to InSight has a probability of only ~0.2, and the odds of capturing it in before and after images are extremely low. We revisit impact-seismic discriminators in light of real experience with a seismometer on the Martian surface. Using measured noise of the instrument, we revise our previous prediction of seismic impact detections downward, from ~a few to tens, to just ~2 per Earth year, still with an order of magnitude uncertainty.

58 GEOSCIENCES↗
DOE OSTI2020

Clustering earthquake signals and background noises in continuous seismic data with unsupervised deep learning

The continuously growing amount of seismic data collected worldwide is outpacing our abilities for analysis, since to date, such datasets have been analyzed in a human-expert intensive, supervised fashion. Moreover, analyses that are conducted can be strongly biased by the standard models employed by seismologists. In response to both of these challenges, we develop a new unsupervised machine learning framework for detecting and clustering seismic signals in continuous seismic records. Our approach combines a deep scattering network and a Gaussian mixture model to cluster seismic signal segments and detect novel structures. To illustrate the power of the framework, we analyze seismic data acquired during the June 2017 Nuugaatsiaq, Greenland landslide. We demonstrate the blind detection and recovery of the repeating precursory seismicity that was recorded before the main landslide rupture, which suggests that our approach could lead to more informative forecasting of the seismic activity in seismogenic areas.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Seismic data are rich in information about subsurface formations and fluids

Seismic attributes are defined as any measured or computed information derived from seismic data. Throughout the last decades extensive work has been done in developing variety of mathematical approaches to extract maximum information from seismic data. Nevertheless, geoscientists found that seismic is still mature and rich in information. In this paper a new seismic attribute is introduced. Instantaneous energy seismic attribute is an amplitude based attribute that has the potential to emphasize anomalous amplitude associated with hydrocarbons. Promising results have been obtained from applying the attribute on seismic section traversing hydrocarbon filled sand from Alberta, Canada.

58 GEOSCIENCES↗
DOE OSTI2021

Profiling the Quito basin (Ecuador) using seismic ambient noise

SUMMARY Quito, the capital of Ecuador, with more than 2.5 M inhabitants, is exposed to a high seismic hazard due to its proximity to the Pacific subduction zone and active crustal faults, both capable of generating significant earthquakes. Furthermore, the city is located in an intermontane piggy-back basin prone to seismic wave amplification. To understand the basin’s seismic response and characterize its geological structure, 20 broad and medium frequency band seismic stations were deployed in Quito’s urban area between May 2016 and July 2018 that continuously recorded ambient seismic noise. We first compute horizontal-to-vertical spectral ratios to determine the resonant frequency distribution in the entire basin. Secondly, we cross-correlate seismic stations operating simultaneously to retrieve interstations surface-wave Green’s functions in the frequency range of 0.1–2 Hz. We find that Love waves travelling in the basin’s longitudinal direction (NNE–SSW) show much clearer correlograms than those from Rayleigh waves. We then compute Love wave phase-velocity dispersion curves and invert them in conjunction with the HVSR curves to obtain shear-wave velocity profiles throughout the city. The inversions highlight a clear difference in the basin’s structure between its northern and southern parts. In the centre and northern areas, the estimated basin depth and mean shear-wave velocity are about 200 m and 1800 ms−1, respectively, showing resonance frequency values between 0.6 and 0.7 Hz. On the contrary, the basement’s depth and shear-wave velocity in the southern part are about 900 m and 2500 ms−1, having a low resonance frequency value of around 0.3 Hz. This difference in structure between the centre-north and the south of the basin explains the spatial distribution of low-frequency seismic amplifications observed during the Mw 7.8 Pedernales earthquake in April 2016 in Quito.

Pacheco, D. (ORCID:0000000282355096)↗
DOE OSTI2015

Discovering new events beyond the catalogue—application of empirical matched field processing to Salton Sea geothermal field seismicity

Using empirical matched field processing (MFP), we compare 4 yr of continuous seismic data to a set of 195 master templates from within an active geothermal field and identify over 140 per cent more events than were identified using traditional detection and location techniques alone. In managed underground reservoirs, a substantial fraction of seismic events can be excluded from the official catalogue due to an inability to clearly identify seismic-phase onsets. Empirical MFP can improve the effectiveness of current seismic detection and location methodologies by using conventionally located events with higher signal-to-noise ratios as master events to define wavefield templates that could then be used to map normally discarded indistinct seismicity. Since MFP does not require picking, it can be carried out automatically and rapidly once suitable templates are defined. In this application, we extend MFP by constructing local-distance empirical master templates using Southern California Earthquake Data Center archived waveform data of events originating within the Salton Sea Geothermal Field. We compare the empirical templates to continuous seismic data collected between 1 January 2008 and 31 December 2011. The empirical MFP method successfully identifies 6249 additional events, while the original catalogue reported 4352 events. The majority of these new events are lower-magnitude events with magnitudes between M0.2–M0.8. Here, the increased spatial-temporal resolution of the microseismicity map within the geothermal field illustrates how empirical MFP, when combined with conventional methods, can significantly improve seismic network detection capabilities, which can aid in long-term sustainability and monitoring of managed underground reservoirs.

15 GEOTHERMAL ENERGY↗
DOE OSTI2016

Case History : Using time-lapse vertical seismic profiling data to constrain velocity-saturation relations: the Frio brine pilot CO 2 injection: Time-lapse VSP to constrain VSRs

CO 2 sequestration projects benefit from quantitative assessment of saturation distribution and plume extent for field development and leakage prevention. In this work, we carry out quantitative analysis of time-lapse seismic by using rock physics and seismic modelling tools. We investigate the suitability of Gassmann's equation for a CO 2 sequestration project with 1600 tons of CO 2 injected into high-porosity, brine-saturated sandstone. We analyze the observed time delays and amplitude changes in a time-lapse vertical seismic profile dataset. Both reflected and transmitted waves are analyzed qualitatively and quantitatively. To interpret the changes obtained from the vertical seismic profile, we perform a 2.5D elastic, finite-difference modelling study. The results show a P-wave velocity reduction of 750 m/s in the proximity of the injection well evident by the first arrivals (travel-time delays and amplitude change) and reflected wave amplitude changes. These results do not match with our rock physics model using Gassmann's equation predictions even when taking uncertainty in CO 2 saturation and grain properties into account. We find that time-lapse vertical seismic profile data integrated with other information (e.g., core and well log) can be used to constrain the velocity–saturation relation and verify the applicability of theoretical models such as Gassmann's equation with considerable certainty. The study shows that possible nonelastic factors are in play after CO 2 injection (e.g., CO 2 –brine–rock interaction and pressure effect) as Gassmann's equation underestimated the velocity reduction in comparison with field data for all three sets of time-lapse vertical seismic profile attributes. Here, our work shows the importance of data integration to validate the applicability of theoretical models such as Gassmann's equation for quantitative analysis of time-lapse seismic data.

58 GEOSCIENCES↗
DOE OSTI2020

Code for Calculating Regional Seismic Travel Time

The RSTT software computes predictions of the travel time of seismic energy traveling from a source to a receiver through 2.5D models of the seismic velocity distribution within the Earth. The two primary applications for the RSTT library are tomographic inversion studies and seismic event location calculations. In tomographic inversions studies, a seismologist begins with number of source-receiver travel time observations and an initial starting model of the velocity distribution within the Earth. A forward travel time calculator, such as the RSTT library, is used to compute predictions of each observed travel time and all of the residuals (observed minus predicted travel time) are calculated. The Earth model is then modified in some systematic way with the goal of minimizing the residuals. The Earth model obtained in this way is assumed to be a better model than the starting model if it has lower residuals. The other major application for the RSTT library is seismic event location. Given an Earth model, an initial estimate of the location of a seismic event, and some number of observations of seismic travel time thought to have originated from that event, location codes systematically modify the estimate of the location of the event with the goal of minimizing the difference between the observed and predicted travel times. The second application, seismic event location, is routinely implemented by the military as part of its effort to monitor the Earth for nuclear tests conducted by foreign countries.

BALLARD, SANFORD↗