Search NASA⌕ Search

SEARCH · Search NASA

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 289 records · Page 16

Induced seismicity during the 2012 Newberry EGS stimulation: Assessment of two advanced earthquake detection techniques at an EGS site

To improve the microearthquake detection capability during the 2012 Newberry enhanced geothermal system (EGS) demonstration, we applied two advanced earthquake detection algorithms to the continuous seismic data. These techniques have the power to detect microearthquakes within the seismic noise that would otherwise be missed using traditional techniques. We compare the efficacy of a matched field processing algorithm with a multi-channel waveform correlation algorithm. Both methods identified between 60%–80% more events, with each slightly outperforming the other in different seismic noise contexts. Incorporating either template-matching technique into a seismic processing routine could therefore effectively increase the sensitivity of the seismic network without changes to the existing hardware. Additionally, lowering the magnitude of completeness would significantly add to the number of events that could subsequently inform future decision making at the reservoir level.

15 GEOTHERMAL ENERGY↗

Three-dimensional seismic characterization and imaging of the Soda Lake geothermal field

Accurate characterization of subsurface geophysical properties and detection of the fault system are essential for geothermal energy exploration and production. The Soda Lake geothermal field is in western Nevada with a complex fault system. Previous seismic characterization only produced a low-resolution, smooth velocity model along with a simple, conceptual fault model. Using optimized correlation-based full-waveform inversion, wavefield-separation-based reverse-time migration, and automatic fault detection techniques, we present 3D seismic characterization for the Soda Lake geothermal field using 3D surface seismic data acquired with Vibroseis sources. Here, we obtain 3D high-resolution velocity, density, and acoustic impedance models, 3D seismic images with different grid spacings, and a high-resolution fault system. Consistency check between the constructed faults and currently active injection and production geothermal wells verifies that our seismic inversion and imaging results and detected faults are reliable. These results can provide valuable information for optimizing well placement and geothermal energy production at the Soda Lake geothermal field.

15 GEOTHERMAL ENERGY↗

Study of seismic diffractions caused by a fracture zone at In Salah carbon dioxide storage project

The In Salah CO 2 storage project in Algeria has injected over 3 million tonnes of carbon dioxide into a water filled tight sand formation. Interferometric Synthetic Aperture Radar (InSAR) range change data revealed a double lobbed pattern of deformation that has been modeled as the opening of a sub-vertical fracture, or damage, zone. The location and geometry of the linear feature were subsequently verified using a seismic reflection survey. The elastic properties of the fracture zone, including anisotropic Poisson ratio ( ν ), Young's ( E ) and shear ( G ) moduli, were estimated from coupled geomechanical and hydrological modeling of surface deformation and pressure variations in the injection well. The elastic moduli reflect the fracture properties after CO 2 flow through the fracture zone. Thus, the seismic signature of the fracture zone could be used for monitoring the CO 2 plume. Using the estimated fracture model, we built two and three dimensional models consisting of an anisotropic fracture zone embedded within an isotropic background. Finite-difference modeling of seismic shot gathers allows us to estimate the effects of scattering from the fracture zone, potentially further constraining the geomechanical model. From the seismic modeling results, we find diffracted waves, induced by the fracture zone, which behave similar to point source diffractions. This modeling is intended to guide a search for diffraction events in the 3D surface seismic field data. The modeling results indicate that using the moduli estimated from geomechanical modeling, fracture scattered events would be 100 times lower amplitude than the interface reflections, and thus would be hard to detect. While diffracted waves are observed in the field data, which may imply the need for revision of the fracture model, including shape and elastic moduli, we are not able to match the field observation with our modeled events. This article presents a frontier study on the integration of geomechanical and geophysical methods at the In Salah site as a means to test the estimate of the subsurface CO 2 flooded fracture properties.

58 GEOSCIENCES↗

Validating compositional fluid flow simulations using 4D seismic interpretation and vice versa in the SECARB Early Test—A critical review

We report on strengths and weaknesses of 4D seismic interpretation as a technique for monitoring carbon dioxide that was injected as part of a large scale test associated with commercial enhanced oil recovery (EOR) at Cranfield Field, Mississippi, USA. The goals of the monitoring effort are 1) to make measurements to verify that the CO 2 is contained in the reservoir according to operational designs and model predictions, and 2) that if there are deviations, to provide data which can be used to update the earth models and determine if any mitigation is needed. Our work uses a compositional numerical simulation to model CO 2 flow in the reservoir and compares the results with estimates of changes in seismic properties between a pre-injection and a survey after more than 2 million metric tons injected. The complicated physics of the problem in Cranfield field present challenges to seismic interpretations. These findings reflect partial agreement between the results of the numerical simulation and the seismic interpretations. Possible causes of discrepancies among the fluid flow model and multiple interpretations of the same seismic data sets performed with different work flows illuminates the types of uncertainties that should be considered to achieve the goals of monitoring.

4D seismic↗

Elastic-wave sensitivity-guided adaptive seismic survey design for cost-effective monitoring of geological carbon storage

Effective seismic monitoring is essential for verifying CO₂ containment, detecting potential leakage, and optimizing operational decisions in geologic carbon storage. Here, this study presents a time-adaptive, elastic-wave sensitivity-guided framework for designing cost-effective seismic monitoring layouts for tracking CO₂ plume migration. The method is based on elastic-wave sensitivity analysis, which quantifies how variations in subsurface properties impact seismic wavefields. Two complementary design strategies are developed: one based on selecting a fixed number of seismic sources (Method A), and the other based on selecting source–receiver pairs contributing to a fixed fraction of cumulative elastic-wave sensitivity energy (Method B). The optimization workflow to identify source–receiver configurations with the highest detection potential is demonstrated using a hypothetical GCS scenario at the Kimberlina site in California using simulations of elastic-wave sensitivity data at multiple post-injection timesteps. Results show that both strategies adapt to evolving plume geometries and wavefield sensitivities, with Method B offering broader spatial coverage and Method A ensuring simpler deployment. This framework enables site-specific, cost-effective, and risk-informed seismic survey designs, enhancing the ability to monitor CO₂ migration over time in evolving geological environments

58 GEOSCIENCES↗

Diffraction imaging for seal evaluation using ultra high resolution 3D seismic data

Cap rock is any impermeable or low permeability formation that may trap oil, gas or water, preventing it from migrating to the surface. Evaluation of sealing properties of the cap rock is a critical task for seismic exploration. It is important for oil and gas prospecting as well as for locating reservoirs appropriate for carbon dioxide storage. The latter has been proposed as one solution to global climate change caused by heat-trapping anthropogenic gasses in the atmosphere. A fluid escaped from the reservoir fills voids in the overlying strata thereby causing local changes of acoustic properties. We aim to explore such fluid-saturated areas and thus evaluate sealing properties of underlying cap rocks. We assume that local acoustic impedance changes caused by fluid migration generate seismic diffractions. We isolate diffracted waves from P-Cable seismic data acquired over the area of interest and perform diffraction analysis of the near-surface interval (above 200 m). In the diffraction image obtained, we interpret an extensive system of broad sinuous channels, narrow linear faults, and diffuse fracture networks. In addition, we detect high diffractivity anomalies, which extend vertically and laterally. We note no evidence of these anomalies in a conventional seismic image and associated seismic attribute (similarity and fracture density) volumes. We associate this high diffractivity to a variable pore fluid composition and saturation, which we infer to be due to hydrocarbon migration from the underlying strata. Consequently, we conclude poor seal quality for the underlying cap rocks. During data processing, we test diffraction focusing analysis for velocity model building. In conclusion, this option appears quite valuable for P-Cable data processing because of short offset and limitations of the standard velocity analysis.

58 GEOSCIENCES↗

Validation of numerical models for seismic fluid-structure-interaction analysis of nuclear, safety-related equipment

Seismic design, qualification, and risk assessment of nuclear safety-related vessels filled with liquid will have to consider the interaction between the vessel (tank), the contained liquid, and submerged components, if any. Seismic fluid-structure-interaction (FSI) analysis of nuclear vessels will rely on numerical models, which are required to be verified and validated. This study validates previously-verified numerical models using test data generated from earthquake-simulator experiments involving a cylindrical tank. Two solvers in LS-DYNA for FSI simulations are used for the numerical analysis: Arbitrary-Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD). Numerical and test results are compared for responses critical to seismic design of advanced reactor vessels: hydrodynamic pressures on the tank wall, reactions at the support, and wave heights of the contained liquid. Analysis is performed for one-, two-, and three-directional seismic inputs with a range of intensity, and rocking motions. The accuracy of the numerical results and different methods for outputting wave heights are discussed. Recommendations for validation of seismic FSI numerical models of advanced reactors are provided. As a result, the validation exercise presented in this study is broadly applicable to cylindrical tanks, regardless of industry sector.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

The Sensitivity of Joint Inversions of Seismic and Geodynamic Data to Mantle Viscosity

Seismic tomography has revealed the existence of large-scale velocity heterogeneities in the mantle. The interpretation of seismic velocity anomalies in terms of temperature and chemical composition is nonunique. We use geodynamic observations including gravity, plate motions, dynamic topography, and excess ellipticity of the core-mantle boundary combined with seismic observations to investigate the thermo-chemical structure of the mantle through joint inversions. An outstanding issue, however, is the physical connection between mantle density anomalies and the surface geodynamic observations, which requires knowledge of the mantle viscosity structure. Here we perform joint inversions assuming different viscosity profiles and examine the dependence of the results on the viscosity. We first assume that mantle heterogeneity is due to thermal variations, which places a constraint on the relation between seismic velocity and density, and we subsequently relax the constraint to allow for potential nonthermal effects. In all of our joint inversions, a nonthermal origin of density anomalies is required to explain the geodynamic data, though the amount varies with the assumed viscosity structure. A common observation is a high-density chemical signal in the center of the large low-shear-velocity provinces at the base of the mantle resulting in a near neutral or slightly dense overall buoyancy there. Using the derived density models and their corresponding viscosity profiles, we also calculate instantaneous mantle flow fields. The predicted flow fields derived from joint inversions are generally similar but are quite different from flow fields using density models derived from a posteriori scaling of pure seismic tomography models.

58 GEOSCIENCES↗

Seasonal modulation of crustal seismicity in northeastern Japan driven by snow load

Numerous studies have reported that surface hydrological loading can seasonally modulate seismicity rates at crustal depths. For example, substantial winter snow accumulation occurs across the Japanese Islands, and these snowy regions appear to have seasonally modulated the occurrence of previous large inland earthquakes. Therefore, it is important to investigate the impact of seasonal stress changes on crustal seismicity to deepen our understanding of earthquake generation. Here we constrain seasonal changes in the surface load across northeastern Japan using Global Navigation Satellite System surface displacements and evaluate the potential relationship between temporal trends in inland seismicity and estimated seasonal stress changes. The spatial distribution of the seasonal surface load is consistent with snow depth along the Sea of Japan. The inland seismicity beneath northeastern Japan is modestly modulated by the seasonal stress changes that are induced by the annual snow load. However, this seasonal response is weaker than that in other regions. This weak modulation may be due to the small surface-load-induced stress perturbation relative to the long-term-averaged stressing rate and/or the limited presence of crustal fluids to trigger seismicity in Japan.

58 GEOSCIENCES↗

Evaluation of seismic hazard of the Gökova bay in terms of seismotectonics

While discovering the seismicity of our country, knowing the array of earthquake occurrence which reflects the characteristic tectonic features of each region makes vital contributions to the earthquakes that have occurred and to the pursuit of the processes which might occur in the future. When considering the region’s seismic activity, the presence of active faults that create earthquake within the bay is obvious. Many active fault parts in the Gulf of Gökova region continues their seismic activity with the opening effect that is generally prevailing in Western Anatolia. The region has generally been continuing its seismic activity under the control of normal faults. Considering the marine studies that are made and marine continuity of the faults which are on land in addition to the seismological and tectonic studies, the determination of seismic hazard in the Gulf of Gökova and its surroundings is also important in terms of introducing the earthquake scenarios with minimized errors.

58 GEOSCIENCES↗

Seismic Performance of Deeply Embedded Conceptual Advanced Reactors: Three-Dimensional Nonlinear Soil-Structure Interaction Analyses

The U.S. Nuclear Regulatory Commission’s functional containment concept provides advanced nuclear power plant designers with more flexibility in terms of the civil/structural design if the appropriate set of barriers for prevention of radioactive material release exist. Some of the conceptual advanced reactor structures, without the traditional pressure boundaries of large containment structures, are proposed to be deeply embedded or buried into soil. This approach is expected to provide (1) lesser seismic demands on the structures and safety-critical structures, (2) eased regulatory efforts and overall design against other external hazards such as aircraft impact, and (3) overall cost savings. Further, one of the important aspects of assessing the technical and economic viability of deeply embedding advanced reactor buildings is to assess the seismic performance with the understanding of effects with material and geometric nonlinearities. This study investigates the seismic response of deeply embedded or buried advanced reactors by conducting three-dimensional nonlinear soil-structure interaction analyses. Although the results indicate that there is a general trend of decreased seismic response with increased embedment depths, the change in the dynamic environment with different embedment depths and the nonlinear environment under high-intensity seismic inputs may result in increased peak response at increased embedment depths.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗

Connect the Dots: In Situ 4-D Seismic Monitoring of CO 2 Storage With Spatio-Temporal CNNs

4-D seismic imaging has been widely used in CO 2 sequestration projects to monitor the fluid flow in the volumetric subsurface region that is not sampled by wells. Ideally, real-time monitoring and near-future forecasting would provide site operators with great insights to understand the dynamics of the subsurface reservoir and assess any potential risks. However, due to obstacles such as high deployment cost, availability of acquisition equipment, exclusion zones around surface structures, only very sparse seismic imaging data can be obtained during monitoring. That leads to an unavoidable and growing knowledge gap over time. The operator needs to understand the fluid flow throughout the project lifetime and the seismic data are only available at a limited number of times. This is insufficient for understanding reservoir behavior. To overcome those challenges, we have developed spatio-temporal neural-network-based models that can produce high-fidelity interpolated or extrapolated images effectively and efficiently. Specifically, our models are built on an autoencoder, and incorporate the long short-term memory (LSTM) structure with a new loss function regularized by optical flow. We validate the performance of our models using real 4-D post-stack seismic imaging data acquired at the Sleipner CO 2 sequestration field. We employ two different strategies in evaluating our models. Numerically, we compare our models with different baseline approaches using classic pixel-based metrics. We also conduct a blind survey and collect a total of 20 responses from domain experts to evaluate the quality of data generated by our models. Finally, via both numerical and expert evaluation, we conclude that our models can produce high-quality 2-D/3-D seismic imaging data at a reasonable cost, offering the possibility of real-time monitoring or even near-future forecasting of the CO 2 storage reservoir.

4-D seismic imaging↗

Elastic-Waveform Inversion with Compressive Sensing for Sparse Seismic Data

Accurate velocity models of compressional- and shear-waves are essential for geothermal reservoir characterization and microseismic imaging. Elastic-waveform inversion of multi-component seismic data can provide high-resolution inversion results of subsurface geophysical properties. However, the method requires seismic data acquired using dense source and receiver arrays. In practice, seismic sources and/or geophones are often sparsely distributed on the surface and/or in a borehole, such as 3D vertical seismic profiling (VSP) surveys. We develop a novel elastic-waveform inversion method with compressive sensing for inversion of sparse seismic data. We employ an alternating-minimization algorithm to solve the optimization problem of our new waveform inversion method. We validate our new method using synthetic VSP data for a geophysical model built using geologic features found at the Raft River enhanced-geothermal-system (EGS) field. We apply our method to synthetic VSP data with a sparse source array and compare the results with those obtained with a dense source array. Our numerical results demonstrate that the velocity models produced with our new method using a sparse source array are almost as accurate as those obtained using a dense source array.

58 GEOSCIENCES↗

A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage (Task 2 Report Extracting Stress Data from Seismic Data)

This report describes research accomplishments achieved with funding provided through Department of Energy (DOE) Contract DE-FE0031686. The purpose of this DOE funding is to develop methods that can provide key information about stress fields that act on deep rocks without invading the earth to acquire that information. The research objective was to demonstrate methods that extract the azimuth directions of SHmax and SHmin stress in deep rocks from traditional seismic reflection data like the data that are used to explore for deep oil and gas reservoirs. Battelle performed two seismic investigations and achieved estimates of SHmax and SHmin azimuths in both efforts that agreed with local, non-seismic, ground-truth measurements of SHmax orientation. The first procedure utilized vertical seismic profiling (VSP) data; the second effort utilized three-dimensional (3D) seismic data.

58 GEOSCIENCES↗

CO 2 Injection Monitoring with a Scalable, Automated Sparse Seismic Array

Conference presentation at American Institute of Chemical Engineers (AIChE) Annual Meeting, Virtual, November 16–20, 2020. The scalable, automated sparse seismic array (SASSA) method has been ideated as a cost-effective alternative or complementary to conventional 3D seismic methods. Based on a sparse array of receivers deployed strategically on Earth’s surface and one or more fixed seismic sources, the SASSA method can be used to study specific portions of the subsurface to detect the movement of CO 2 in a reservoir with lower environmental footprint and shorter turnaround time to deliver results than conventional 3D seismic methods. Energy & Environmental Research Center (EERC) carried out seismic monitoring of an oilfield undergoing CO 2 enhanced oil recovery in Montana, USA, using the SASSA method.

02 PETROLEUM↗

Numerical Modeling of Anisotropic Elastic-Wave Sensitivity Propagation for Optimal Design of Time-Lapse Seismic Surveys

Reliable subsurface time-lapse seismic monitoring is crucial for many geophysical applications, such as enhanced geothermal system characterization, geologic carbon utilization and storage, and conventional and unconventional oil/gas reservoir characterization, etc. In this work, we develop an elastic-wave sensitivity propagation method for optimal design of cost-effective time-lapse seismic surveys considering the fact that most of subsurface geologic layers and fractured reservoirs are anisotropic instead of isotropic. For anisotropic media, we define monitoring criteria using qP- and qS-wave sensitivity energies after decomposing qP- and qS-wave components from the total elastic-wave sensitivity wavefield using a hybrid time- and frequency-domain approach. Geophones should therefore be placed at locations with significant qP- and qS-wave sensitivity energies for cost-effective time-lapse seismic monitoring in an anisotropic geology setting. Our numerical modeling results for a modified anisotropic Hess model demonstrate that, compared with the isotropic case, subsurface anisotropy changes the spatial distributions of elastic-wave sensitivity energies. Consequently, it is necessary to consider subsurface anisotropies when designing the spatial distribution of geophones for cost-effective time-lapse seismic monitoring. This finding suggests that it is essential to use our new anisotropic elastic-wave sensitivity modeling method for optimal design of time-lapse seismic surveys to reliably monitor the changes in subsurface reservoirs, fracture zones or target monitoring regions.

97 MATHEMATICS AND COMPUTING↗

Seismic Exploration of Small Bodies

As a result of the Phase I study, we have demonstrated that a mission to a smallbody (asteroid, comet) whose objective is to conduct a seismic experiment to understandthe interior structure, can be accomplished with small (<200 kg) spacecraft launched on asmall launch vehicle (Athena IIc).We modeled the seismic response of a small body and calculated that the energynecessary to propagate through the body and be detected by a seismometer. Thesecalculations provide guidance as to the type of energy source that is required. A simpleenergy source similar to a NASA standard initiator (NSI) can be used, although a singleNSI is insufficient. An NSI is an explosive pyrotechnic that is used to sever connectionson spacecraft. Use of such an energy source has illustrated two additional areas of study -anchoring of the source and sensor to the surface and understanding the efficiency ofenergy propagation from the source into the surface.The spacecraft has the ability to carry and deploy a series of source/sensors toconduct the experiment by placing them on the surface. We identified a suite of candidatenear-Earth asteroids as targets and used one 1991VG as the target to calculate the missiontrajectory and ∆V requirements. Spacecraft and launch vehicle performance are launchthat any of the candidates could be reached with appropriate mass and launch margins.Sources and sensors are deployed from arm attached to the spacecraft. Thespacecraft maneuvers next to the target body and presses the sensor against the surfaceand releases it. After emplacing all of the surface packages, the sensors are monitored fora period of time to measure the seismic noise. Finally, the active seismic experiment isconducted.The work conducted during Phase I demonstrates that a small mission can be designed to conduct an active seismic experiment on a small near-Earth body. Whileother targets may require more performance, the basic architecture is viable for anytarget.We have identified a number of specific technical areas that require more detailed study. Those areas are largely focused on the detailed analysis and design of the source mechanism and anchoring it to the surface. Such topics will be part of a Phase IIproposal.

Spacecraft↗