Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SEISMIC”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Utah FORGE: Composite 3D Seismic Velocity Model

This is a composite 3D seismic velocity that was constructed from compiled information from several local studies regarding seismic velocities and structural information. This seismic velocity model is provided in NonLinLoc format (slow_len), which is readily usable in NonLinLoc software. Other model formats and versions of the model can be produced using the Python script provided with this data set. Details on how the model was created and prior velocity and structural information was used is provided in the accompanying documentation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Utah FORGE: Seismic Velocity Models, February 2021

This dataset contains a map, showing the Utah FORGE seismic stations, and seismic velocity model data. There are 61 1-D velocity models which are in a compressed TAR file. A paper is referenced at the end of this description which discusses the use of these data in 3D modelling. The paper summary follows: We expand the application of spatial autocorrelation (SPAC) from typical 1-D Vs profiles to quasi-3-D imaging via Bayesian Monte Carlo inversion (BMCI) using a dense nodal array (49 nodes) located at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site. Combinations of 4 and 9 geophones in subarrays provide for 36 and 25 1-D Vs profiles, respectively. Profiles with error bars are determined by calculating coherency functions that fit observations in a frequency range of 0.2-5 Hz. Thus, a high-resolution quasi-3-D Vs model from the surface to 2.0 km depth is derived and shows that surface-parallel sedimentary strata deepen to the west, consistent with a 3-D seismic reflection survey. Moreover, the resulting Vs profile is consistent with a Vs profile derived from distributed acoustic sensing (DAS) data located in a borehole at the FORGE site. The quasi-3-D velocity model shows that the base of the basin dips ~22 degrees to the west and topography on the basement interface coincident with the Mag Lee Wash suggests that the bedrock interface is an unconformity. Reference: Zhang, H. and K. L. Pankow (2021). High-resolution Bayesian spatial auto-correlation (SPAC) pseudo-3D Vs model of Utah FORGE site with a dense geophone array, Geophys. Res. Int, https://doi.org/10.1093/gji/ggab049

15 GEOTHERMAL ENERGY↗
DOE OSTI2015

On the use of remote infrasound and seismic stations to constrain the eruptive sequence and intensity for the 2014 Kelud eruption

The February 2014 eruption of Kelud volcano (Indonesia) destroyed most of the instruments near it. We use remote seismic and infrasound sensors to reconstruct the eruptive sequence. The first explosions were relatively weak seismic and infrasound events. A major stratospheric ash injection occurred a few minutes later and produced long-lasting atmospheric and ground-coupled acoustic waves that were detected as far as 11,000 km by infrasound sensors and up to 2300 km away on seismometers. A seismic event followed ~12 minutes later and was recorded 7000 km away by seismometers. We estimate a volcanic intensity around 10.9, placing the 2014 Kelud eruption between the 1980 Mount St. Helens and 1991 Pinatubo eruptions intensities. As a result, we demonstrate how remote infrasound and seismic sensors are critical for the early detection of volcanic explosions, and how they can help to constrain and understand eruptive sequences.

58 GEOSCIENCES↗
DOE OSTI2017

Spatiotemporal changes of seismic attenuation caused by injected CO 2 at the Frio-II pilot site, Dayton, TX, USA

A continuous active source seismic monitoring data set was collected with crosswell geometry during CO 2 injection at the Frio-II brine pilot, near Liberty, TX. Previous studies have shown that spatiotemporal changes in the P wave first arrival time reveal the movement of the injected CO 2 plume in the storage zone. To further constrain the CO 2 saturation, particularly at higher saturation levels, we investigate spatial-temporal changes in the seismic attenuation of the first arrivals. The attenuation changes over the injection period are estimated by the amount of the centroid frequency shift computed by local time-frequency analysis. We observe that (1) at receivers above the injection zone seismic attenuation does not change in a physical trend; (2) at receivers in the injection zone attenuation sharply increases following injection and peaks at specific points varying with distributed receivers, which is consistent with observations from time delays of first arrivals; then, (3) attenuation decreases over the injection time. The attenuation change exhibits a bell-shaped pattern during CO 2 injection. Under Frio-II field reservoir conditions, White's patchy saturation model can quantitatively explain both the P wave velocity and attenuation response observed. We have combined the velocity and attenuation change data in a crossplot format that is useful for model-data comparison and determining patch size. Our analysis suggests that spatial-temporal attenuation change is not only an indicator of the movement and saturation of CO 2 plumes, even at large saturations, but also can quantitatively constrain CO 2 plume saturation when used jointly with seismic velocity.

58 GEOSCIENCES↗
DOE OSTI2021

Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors

Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Developing a Consistent Travel-Time Framework for Comparing Three-Dimensional Velocity Models for Seismic Location Accuracy

Abstract Location algorithms have historically relied on simple, one-dimensional (1D) velocity models for fast seismic event locations. 1D models are generally used as travel-time lookup tables, one for each seismic phase, with travel-times pre-calculated for event distance and depth. These travel-time lookup tables are extremely fast to use and this fast computational speed makes them the preferred type of velocity model for operational needs. Higher-dimensional (i.e., three-dimensional—3D) seismic velocity models are becoming readily available and provide more accurate event locations over 1D models. The computational requirements of these 3D models tend to make their operational use prohibitive. Additionally, comparing location accuracy for 3D seismic velocity models tends to be problematic, as each model is determined using different ray-tracing algorithms. Attempting to use a different algorithm than the one used to develop a model usually results in poor travel-time prediction. We demonstrate and test a framework to create first-P and first-S 3D travel-time correction surfaces using an open-source framework ( PCalc + GeoTess , https://www.sandia.gov/salsa3d/software/geotess ) that easily stores 3D travel-time and uncertainty data. This framework produces fast travel-time and uncertainty predictions and overcomes the ray-tracing algorithm hurdle because the lookup tables can be generated using the exact ray-tracing algorithm that is preferred for a model.

3D velocity models↗
DOE OSTI2018

Water-table and discharge changes associated with the 2016–2017 seismic sequence in central Italy: hydrogeological data and a conceptual model for fractured carbonate aquifers

A seismic sequence in central Italy from August 2016 to January 2017 affected groundwater dynamics in fractured carbonate aquifers. Changes in spring discharge, water-table position, and streamflow were recorded for several months following nine Mw 5.0–6.5 seismic events. Data from 22 measurement sites, located within 100 km of the epicentral zones, were analyzed. The intensity of the induced changes were correlated with seismic magnitude and distance to epicenters. The additional post-seismic discharge from rivers and springs was found to be higher than 9 m{sup 3}/s, totaling more than 0.1 km{sup 3} of groundwater release over 6 months. This huge and unexpected contribution increased streamflow in narrow mountainous valleys to previously unmeasured peak values. Analogously to the L’Aquila 2009 post-earthquake phenomenon, these hydrogeological changes might reflect an increase of bulk hydraulic conductivity at the aquifer scale, which would increase hydraulic heads in the discharge zones and lower them in some recharge areas. The observed changes may also be partly due to other mechanisms, such as shaking and/or squeezing effects related to intense subsidence in the core of the affected area, where effects had maximum extent, or breaching of hydraulic barriers.

58 GEOSCIENCES↗
DOE OSTI2016

TOUGH2-seed: A coupled fluid flow and mechanical-stochastic approach to model injection-induced seismicity

Understanding the injection-induced triggering mechanism is a fundamental step towards controlling the seismicity generated by deep underground exploitation. Here we propose a modeling approach based on coupling the TOUGH2 simulator with a geomechanical-stochastic model. Additionally, the hydro-mechanical-stochastic model provides a good representation of different mechanisms influencing each other during and after the injection phase. Each mechanism affects the induced seismicity in a different way and at different times during the reservoir stimulation, confirming that a complex interaction is in place, and that more sophisticated and physics-based approaches coupled with statistical model are required to explain such a complex interaction. In addition to previous statistical and hybrid models, our approach accounts for a full 3D formulation of both stresses and fluid flow, further including all the TOUGH2 capabilities. Furthermore, it includes interactions between triggered seismic events through calculation of static stress transfer. Here, in this work, we present the main capabilities of TOUGH2-SEED and apply the model to the Basel EGS case, successfully reproducing the injection pressure as well as the evolution of the seismicity.

58 GEOSCIENCES↗
DOE OSTI2016

Effects of Asperity Distribution on Fluid Flow and Induced Seismicity During Deep Geothermal Exploitation

This work investigates the injection-induced seismic response of a heterogeneous fault plane, featuring low-permeability asperities embedded into a high-permeability damage zone. We simulate the pressure evolution with a hydrogeological simulator, accounting for the heterogeneous fault plane. Seismicity occurs then on the asperities, represented as unstable patches reactivating by means of the Mohr-Coulomb criterion. The hydrological and seismic modules are implicitly coupled to account for effects of asperity reactivation on the permeability. Results show that permeability changes may cause at a later time a change in seismicity propagation. Finally, we also investigated such effects by varying the density of asperities.

58 GEOSCIENCES↗
DOE OSTI2021

Review of analytical studies on seismic fluid-structure interaction of base-supported cylindrical tanks

Analytical solutions for seismic fluid-structure interaction (FSI) of tanks are used by analysts and engineers for a number of purposes, including preliminary sizing, risk assessment, and verification of numerical models. Analytical solutions for seismic FSI of tanks were first developed in the 1930s and have been extended and/or modified in the following nine decades to accommodate different seismic inputs and boundary conditions. This report reviews published analytical solutions for rigid and flexible, base-supported, cylindrical tanks subjected to unidirectional horizontal motion of a small amplitude. The studies parsed the fluid-structure (FSI) response into an impulsive and a convective component, and addressed frequencies of lateral motions of the tank (i.e., impulsive frequencies), hydrodynamic pressures, frequencies of waves (i.e., convective frequencies), wave heights, and reactions at the base. The solutions were categorized into exact and approximate solutions based on the methodologies used. Herein, FSI responses are calculated and compared for a range of tank dimensions using the analytical solutions from the different studies. The responses presented here are normalized to be unitless, and can be used for tanks with different dimensions and mechanical properties, and subjected to different input motions. A sample steel water tank is analyzed using the normalized solutions for two earthquake ground motions with different frequency contents. Ignoring tank flexibility may significantly underestimate seismic responses of the tank and its contained fluid. The approximate solutions, which have been widely used for decades, may not be sufficiently accurate in the modern era. Calculation errors in the reviewed studies are identified and corrected, and differences are quantified for the sample tank.

42 ENGINEERING↗
DOE OSTI2021

The Water-Fe-Pressure dependent single-crystal elastic properties of wadsleyite: Implications for the seismic anisotropy in the upper Mantle Transition Zone

Recent seismic studies suggest an anisotropic Mantle Transition Zone (MTZ) in areas adjacent to subducted slabs. Wadsleyite is the main anisotropy contributor in the upper MTZ, therefore the interpretation of these seismic observations requires the knowledge of single-crystal elastic moduli (C ij s) and the deformation-induced lattice preferred orientation (LPO) of wadsleyite. Wadsleyite can host up to 3 wt% water in its crystal structure as point defects in the form of hydroxyl groups, however, the combined effect of water content, Fe content, and pressure on the C ij s of wadsleyite remains unclear. In this study, we measured the high-pressure single-crystal elasticity of a synthetic hydrous Fe-bearing wadsleyite (0.14(5) wt% water, Fe#=9.4, Fe 3+ /ΣFe=0.3) up to 18.2 (2) GPa. In combination with previous experimental data, we separated the effects of pressure, water, and Fe contents on the Cijs and intrinsic elastic anisotropy of wadsleyite. Overall, our results suggest that the intrinsic elastic anisotropy of wadsleyite decreases with pressure, water, and Fe contents. At 15 GPa, increasing the water content by 0.1 wt% or Fe# by 1 decreases the V P and V S anisotropy of wadsleyite by ~1.1-1.3%, and ~0.8-1.3% in average, respectively. Combined with the LPO determined in previous deformation experiments, we modeled the seismic anisotropy in the upper MTZ generated by a sub-vertical mantle flow near cold subducted slabs and a sub-horizontal mantle flow in the ambient mantle. In both scenarios, the LPO of wadsleyite leads to V SV (vertically polarized shear wave velocity) > V SH (horizontally polarized shear wave velocity). Our results suggest that wadsleyite may account for a weak anisotropic MTZ (<1%) on the global scale. Considering the fact that water decreases the elastic anisotropy but promotes LPO of wadsleyite, seismic anisotropy may not be a good water sensor.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A rock physics and seismic reservoir characterization study of the Rock Springs Uplift, a carbon dioxide sequestration site in Southwestern Wyoming

We present a reservoir geophysics study, including rock physics modeling and seismic inversion, of a carbon dioxide sequestration site in Southwestern Wyoming, namely the Rock Springs Uplift, and build a petrophysical model for the potential injection reservoirs for carbon dioxide sequestration. Our objectives include the facies classification and the estimation of the spatial model of porosity and permeability for two sequestration targets of interest, the Madison Limestone and the Weber Sandstone. The available dataset includes a complete set of well logs at the location of the borehole available in the area, a set of 110 core samples, and a seismic survey acquired in the area around the well. The proposed study includes a formation evaluation analysis and facies classification at the well location, the calibration of a rock physics model to link petrophysical properties and elastic attributes using well log data and core samples, the elastic inversion of the pre-stack seismic data, and the estimation of the reservoir model of facies, porosity and permeability conditioned by seismic inverted elastic attributes and well log data. In particular, the rock physics relations are facies-dependent and include granular media equations for clean and shaley sandstone, and inclusion models for the dolomitized limestone. The permeability model has been computed by applying a facies-dependent porosity-permeability relation calibrated using core sample measurements. Finally, the study shows that both formations show good storage capabilities. The Madison Limestone includes a homogeneous layer of high-porosity high-permeability dolomite; the Weber Sandstone is characterized by a lower average porosity but the layer is thicker than the Madison Limestone.

58 GEOSCIENCES↗
DOE OSTI2020

CO 2 induced changes in Mount Simon sandstone: Understanding links to post CO 2 injection monitoring, seismicity, and reservoir integrity

The purpose of this study is to quantify geochemical reactions of CO 2 and brine with subsurface samples taken from the Mt. Simon sandstone and identify any potential alterations of the geomechanical rock properties that could lead to changes observable in seismic monitoring or result in changes of micro seismicity such as those observed at the Illinois Basin-Decatur Project (IBDP) site. Two Mt. Simon sandstone samples from 6919.3 feet (2109.0 m) and 6926.1 feet (2111.1 m) depth were exposed to supercritical CO 2 (scCO 2 ) dissolved in brine at in-situ reservoir conditions for one month. Geochemical, spectral, scanning electron microscopy, and petrophysical methods were used to analyze the samples before and after the one-month exposure. Significant changes were observed. Multiple geomechanical properties were chosen to form the framework against which to interrogate the petrophysical data: Young’s modulus (E), Poisson’s ratio (ν), lambda·rho (λρ), and mu·rho (μρ). In this study we conclude that framework composition, porosity, heterogeneities, effective pressure, and reactive geochemistry are first order controls on trends in the E-μ and λρ-μρ cross plot spaces. Changes in porosity, permeability, dynamic moduli, and brittleness with exposure to these fluids were observed. No change in ultrasonic P-wave attenuation (Q P ) was observed. Geochemical alteration causes a distinct shift in λρ-μρ in both samples as well as changes in E, ν, and P and S wave seismic velocity values. Furthermore, these observations could provide insight into subsurface monitoring using seismic methods including amplitude variation with offset (AVO) classification.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity

Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated with simultaneous modeling of the two failure phenomena because stress and deformation evolve non-linearly in both plasticity and seismicity. However, in reservoirs undergoing long-term production, plastic failure can alter the stress paths of points on a fault such that the onset, location, and magnitude of actual seismic events can no longer be predicted by a poroelastic simulation due to inaccurate stress and deformation history. We present a computational framework for coupled multiphase flow, finite strain poroplastic deformation, and dynamic fault slip and use it to understand the impact of plastic deformation on the onset, location, and magnitude of induced fault slip events. We evaluate the impact of plasticity on reservoir pressure, deformation, induced stress, and fault slip by comparing infinitesimal strain elastic and finite strain poro-elastoplastic models. For real-world applications, we consider different scenarios where the reservoir is either mechanically weaker or stronger than the caprock. We analyze the stress evolution as a function of the change in reservoir pressure to understand the role of contrast in reservoir and caprock elastic moduli on geomechanical stability. The results show that the poroplastic reservoir exhibits larger vertical deformation and delayed slip than the poroelastic reservoir after the same amount of oil production. For the same amount of pressure drop, a reservoir with a smaller modulus than the caprock displays a larger vertical displacement and an earlier onset of both plastic failure and fault slip. For a reservoir with a larger modulus than the caprock, vertical displacement is larger on the reservoir top boundary and smaller on the ground surface, and a higher pressure drop is needed to induce plastic failure and fault slip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Modeling of fault activation and seismicity by injection directly into a fault zone associated with hydraulic fracturing of shale-gas reservoirs

We conducted three-dimensional coupled fluid-flow and geomechanical modeling of fault activation and seismicity associated with hydraulic fracturing stimulation of a shale-gas reservoir. We simulated a case in which a horizontal injection well intersects a steeply dip- ping fault, with hydraulic fracturing channeled within the fault, during a 3-hour hydraulic fracturing stage. Consistent with field observations, the simulation results show that shale-gas hydraulic fracturing along faults does not likely induce seismic events that could be felt on the ground surface, but rather results in numerous small microseismic events, as well as aseismic deformations along with the fracture propagation. The calculated seismic moment magnitudes ranged from about -2.0 to 0.5, except for one case assuming a very brittle fault with low residual shear strength, for which the magnitude was 2.3, an event that would likely go unnoticed or might be barely felt by humans at its epicenter. The calculated moment magnitudes showed a dependency on injection depth and fault dip. We attribute such dependency to variation in shear stress on the fault plane and associated variation in stress drop upon reactivation. Our simulations showed that at the end of the 3-hour injection, the rupture zone associated with tensile and shear failure extended to a maximum radius of about 200 m from the injection well. The results of this modeling study for steeply dipping faults at 1000 to 2500 m depth is in agreement with earlier studies and field observations showing that it is very unlikely that activation of a fault by shale-gas hydraulic fracturing at great depth (thousands of meters) could cause felt seismicity or create a new flow path (through fault rupture) that could reach shallow groundwater resources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Coupled Hydromechanical Modeling of Induced Seismicity From CO 2 Injection in the Illinois Basin

Injection of CO 2 for geologic carbon sequestration into deep sedimentary formations involves fluid pressure increases that engage hydromechanical processes that can cause seismicity by activation of existing faults. In this work, we use a coupled multiphase fluid flow and geomechanical simulator to model spatiotemporal fluid pressure and stress changes in order to study the poroelastic effect of CO 2 injection on faults in crystalline basement rock below the injection zone. The seismicity rate along features interpreted to be basement faults is modeled using Dieterich's rate-and-state earthquake nucleation model. The methodology is applied to microseismicity detected during CO 2 injection into the Mount Simon formation during the Illinois Basin—Decatur Project. The modeling accurately captures an observed reduction in seismicity rate when the injection in the second well was into a slightly shallower zone above the base of the Mount Simon formation. Moreover, the modeling shows that it is important to consider poroelastic stress changes, in addition to fluid pressure changes for accurately modeling of the observed seismicity rate.

58 GEOSCIENCES↗
DOE OSTI2020

Analytical Solutions for Seismic Fluid-Structure Interaction of Head-Supported Cylindrical Tanks

Fluid-structure interaction (FSI) of earthquake shaking is a design consideration in the civil, nuclear, chemical, and mechanical industries. Preliminary sizing of structures and components in these industries will often use analytical solutions because there is insufficient information to warrant the analysis of numerical models. Here, this study extends prior analytical work on seismic FSI of base-supported tanks in the 1980s to accommodate head-supported tanks with geometries similar to those proposed for new-build, liquid (fluid)-filled, advanced reactors. Analytical solutions are derived for flexible, head-supported, cylindrical tanks subjected to small-amplitude, unidirectional, horizontal seismic motion, and address frequencies of lateral motions of the tank, hydrodynamic pressures, frequencies of waves, wave heights, and reactions at the head support. Solutions are presented for a range of tank dimensions and could be used for a) preliminary design of head-supported cylindrical tanks and their supporting structures with idealized seismic input, and b) verification of numerical models to be used for final seismic analysis, design, and probabilistic risk assessment.

42 ENGINEERING↗
DOE OSTI2015

Earthquake location determination using data from DOMERAPI and BMKG seismic networks: A preliminary result of DOMERAPI project

DOMERAPI project has been conducted to comprehensively study the internal structure of Merapi volcano, especially about deep structural features beneath the volcano. DOMERAPI earthquake monitoring network consists of 46 broad-band seismometers installed around the Merapi volcano. Earthquake hypocenter determination is a very important step for further studies, such as hypocenter relocation and seismic tomographic imaging. Ray paths from earthquake events occurring outside the Merapi region can be utilized to delineate the deep magma structure. Earthquakes occurring outside the DOMERAPI seismic network will produce an azimuthal gap greater than 180{sup 0}. Owing to this situation the stations from BMKG seismic network can be used jointly to minimize the azimuthal gap. We identified earthquake events manually and carefully, and then picked arrival times of P and S waves. The data from the DOMERAPI seismic network were combined with the BMKG data catalogue to determine earthquake events outside the Merapi region. For future work, we will also use the BPPTKG (Center for Research and Development of Geological Disaster Technology) data catalogue in order to study shallow structures beneath the Merapi volcano. The application of all data catalogues will provide good information as input for further advanced studies and volcano hazards mitigation.

58 GEOSCIENCES↗