Effect of fault zone architecture on injection-induced seismicity.
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
Seismic velocities and elastic moduli of rocks are known to vary significantly with applied stress, which indicates that these materials exhibit nonlinear elasticity. Monochromatic waves in nonlinear elastic media are known to generate higher harmonics and combinational frequencies. Such effects have the potential to be used for broadening the frequency band of seismic sources, characterization of the subsurface, and safety monitoring of civil engineering infrastructure. However, knowledge on nonlinear seismic effects is still scarce, which impedes the development of their practical applications. To explore the potential of nonlinear seismology, we performed three experiments: two in the field and one in the laboratory. The first field experiment used two vibroseis sources generating signals with two different monochromatic frequencies. The second field experiment used a surface orbital vibrator with two eccentric motors working at different frequencies. In both experiments, the generated wavefield was recorded in a borehole using a fiber-optic distributed acoustic sensing cable. Both experiments showed combinational frequencies, harmonics, and other intermodulation products of the fundamental frequencies both on the surface and at depth. Laboratory experiments replicated the setup of the field test with vibroseis sources and showed similar nonlinear combinations of fundamental frequencies. Amplitudes of the nonlinear signals observed in the laboratory showed variation with the saturating fluid. These results confirm that nonlinear components of the wavefield propagate as body waves, are likely to generate within rock formations, and can be potentially used for reservoir fluid characterization.
Abstract not provided.
The seismic design of an advanced nuclear reactor must consider the interaction of vessel internal components with the surrounding coolant: fluid–structure interaction (FSI). Available analytical solutions for FSI of submerged components do not accommodate multiple-component, intense seismic inputs and complex reactor and internal geometries. Physical testing of reactor vessels and internals to inform seismic design is impractical and cost-prohibitive, leaving the use of verified and validated, robust numerical models as the only plausible option for analysis and design. Physical data that could be used for validating such numerical models for multi-component shaking are not available. This article describes an experimental program performed on a 6-degree-of-freedom earthquake simulator to generate data that could support validation of seismic FSI numerical models for submerged components in commercial finite element codes. A scaled model of a base-supported reactor vessel with simplified representations of submerged internals was tested to generate submerged-component response histories for a range of seismic inputs. The generated data were used to validate numerical models in the finite element code LS-DYNA. Numerical models were validated for calculation of hydrodynamic pressure responses on internals, in-water frequencies of internals, and acceleration responses of internals. The generated data and the analysis recommendations could aid engineering analysts designing submerged components and systems for seismic effects.
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
The spectacular Lumpur Sidoarjo (Lusi) eruption started in northeast Java on the 29 of May 2006 following a M6.3 earthquake striking the island [1,2]. Initially, several gas and mud eruption sites appeared along the reactivated strike-slip Watukosek fault system [3] and within weeks several villages were submerged by boiling mud. The most prominent eruption site was named Lusi. The Lusi seismic experiment is a project aims to begin a detailed study of seismicity around the Lusi area. In this initial phase we deploy 30 seismometers strategically distributed in the area around Lusi and along the Watukosek fault zone that stretches between Lusi and the Arjuno Welirang (AW) complex. The purpose of the initial monitoring is to conduct a preliminary seismic campaign aiming to identify the occurrence and the location of local seismic events in east Java particularly beneath Lusi.This network will locate small event that may not be captured by the existing BMKG network. It will be crucial to design the second phase of the seismic experiment that will consist of a local earthquake tomography of the Lusi-AW region and spatial and temporal variations of vp/vs ratios. The goal of this study is to understand how the seismicity occurring along the Sunda subduction zone affects to the behavior of the Lusi eruption. Our study will also provide a large dataset for a qualitative analysis of earthquake triggering studies, earthquake-volcano and earthquake-earthquake interactions. In this study, we will extract Green’s functions from ambient seismic noise data in order to image the shallow subsurface structure beneath LUSI area. The waveform cross-correlation technique will be apply to all of recordings of ambient seismic noise at 30 seismographic stations around the LUSI area. We use the dispersive behaviour of the retrieved Rayleigh waves to infer velocity structures in the shallow subsurface.
This report summarizes the results of experimental and numerical studies conducted to evaluate the effect of accident thermal loads on the seismic response of steel-plate composite (SC) walls and reinforced concrete (RC) walls and beams. The research was sponsored by U.S. Department of Energy through the Nuclear Energy University Program (NEUP), Award Number DE-NE0008256. Summary and conclusions of the research conducted are presented in this report. Based on the observations, analysis and design recommendations for SC and RC structures are also summarized
This paper discusses a new approach for investigating the seismic wave velocity of rock, specifically carbonates, as affected by their pore structures. While the conventional routine of seismic velocity measurement highly depends on the extensive laboratory experiment, the proposed approach utilizes the digital rock physics view which lies on the numerical experiment. Thus, instead of using core sample, we use the thin section image of carbonate rock to measure the effective seismic wave velocity when travelling on it. In the numerical experiment, thin section images act as the medium on which wave propagation will be simulated. For the modeling, an advanced technique based on artificial neural network was employed for building the velocity and density profile, replacing image's RGB pixel value with the seismic velocity and density of each rock constituent. Then, ultrasonic wave was simulated to propagate in the thin section image by using finite difference time domain method, based on assumption of an acoustic-isotropic medium. Effective velocities were drawn from the recorded signal and being compared to the velocity modeling from Wyllie time average model and Kuster-Toksoz rock physics model. To perform the modeling, image analysis routines were undertaken for quantifying the pore aspect ratio that is assumed to represent the rocks pore structure. In addition, porosity and mineral fraction required for velocity modeling were also quantified by using integrated neural network and image analysis technique. It was found that the Kuster-Toksoz gives the closer prediction to the measured velocity as compared to the Wyllie time average model. We also conclude that Wyllie time average that does not incorporate the pore structure parameter deviates significantly for samples having more than 40% porosity. Utilizing this approach we found a good agreement between numerical experiment and theoretically derived rock physics model for estimating the effective seismic wave velocity of rock.
Abstract The primary phase of the Earth’s lower mantle, (Al, Fe)‐bearing bridgmanite, transitions to the post‐perovskite (PPv) phase at Earth’s deep mantle conditions. Despite extensive experimental and ab initio investigations, there are still important aspects of this transformation that need clarification. Here, we address this transition in (Al 3+ , Fe 3+ )‐, (Al 3+ )‐, (Fe 2+ )‐, and (Fe 3+ )‐bearing bridgmanite using ab initio calculations and validate our results against experiments on similar compositions. Consistent with experiments, our results show that the onset transition pressure and the width of the two‐phase region depend distinctly on the chemical composition: (a) Fe 3+ ‐, Al 3+ ‐, or (Al 3+ , Fe 3+ )‐alloying increases the transition pressure, while Fe 2+ ‐alloying has the opposite effect; (b) in the absence of coexisting phases, the pressure‐depth range of the Pv‐PPv transition is likely too broad to cause a sharp D” discontinuity (<30 km); (c) the average Clapeyron slope of the two‐phase regions are consistent with previous measurements, calculations in MgSiO 3 , and inferences from seismic data. In addition, (d) we observe a softening of the bulk modulus in the two‐phase region. The consistency between our results and experiments gives us the confidence to proceed and examine this transition in aggregates with different compositions computationally, which will be fundamental for resolving the most likely chemical composition of the D" region by analyses of tomographic images.
Explore the source record for details and available documents.
Not Available
Not Available
Abstract not provided.
Most earth materials are anisotropic with regard to seismic wave-speeds, especially materials such as shales, or where oriented fractures are present. However, the base assumption for many numerical simulations is to treat earth materials as isotropic media. This is done for simplicity, the apparent weakness of anisotropy in the far field, and the lack of well-characterized anisotropic material properties for input into numerical simulations. One approach for addressing the higher complexity of actual geologic regions is to model the material as an orthorhombic medium. We have developed an explicit time-domain, finite-difference (FD) algorithm for simulating three-dimensional (3D) elastic wave propagation in a heterogeneous orthorhombic medium. The objective of this research is to investigate the errors and biases that result from modeling a non-isotropic medium as an isotropic medium. This is done by computing “observed data” by using synthetic, anisotropic simulations with the assumption of an orthorhombic, anisotropic earth model. Green’s functions for an assumed isotropic earth model are computed and then used an inversion designed to estimate moment tensors with the “observed” data. One specific area of interest is how shear waves, which are introduced in an anisotropic model even for an isotropic explosion, affect the characterization of seismic sources when isotropic earth assumptions are made. This work is done in support of the modeling component of the Source Physics Experiment (SPE), a series of underground chemical explosions at the Nevada National Security Site (NNSS).
Multiple U.S. national laboratories, universities and industrial collaborators are conducting collaborative research under the EGS Collab project supported by the U.S. Department of Energy, to understand the fracture creation and imaging during fracturing in enhanced geothermal systems. Microearthquake hypocenter locations and focal mechanisms are used to monitor hydraulic fracturing growth and characterization at the EGS Collab experimental site at the Sanford Underground Research Facility using seismic receivers in multiple monitoring wells. In this study, we develop a methodology for optimal design a 3D borehole seismic array for cost-effective seismic monitoring in anisotropic media using not only the relationship between receiver distributions and standard deviation errors of microearthquake hypocenter locations, but also that between receiver distributions and focal-mechanism inversion errors. Our results indicate that microearthquake hypocenter locations and focal mechanisms can be reasonably well reconstructed for the EGS Collab Experiment I using six monitoring wells, including four fracture-parallel monitoring wells and two orthogonal wells. Lastly, eight seismic receivers evenly distributed in four parallel monitoring wells or twelve receivers in all six monitoring wells are required for hypocenter location, and twelve receivers evenly distributed in six wells or sixteen receivers in four wells are needed for focal-mechanism inversion.
A classification of the equipment important for the safety of nuclear power plants is proposed in terms of its dynamic behavior under seismic loading. An extended bank of data from dynamic tests over the entire range of thermal and mechanical equipment in generating units with VVER-1000 and RBMK-1000 reactors is analyzed. Results are presented from a study of the statistical behavior of the distribution of vibrational frequencies and damping decrements with the “small perturbation” factor that affects the measured damping decrements taken into account. A need to adjust the regulatory specifications for choosing the values of the damping decrements with specified inertial loads on equipment owing to seismic effects during design calculations is identified. Minimum values of the decrements are determined and proposed for all types of equipment as functions of the directions and natural vibration frequencies of the dynamic interactions to be adopted as conservative standard values in the absence of actual experimental data in the course of design studies of seismic resistance.
Currently the Department of Energy (DOE) and the nuclear industry perform seismic soil-structure interaction (SSI) analysis using equivalent linear numerical analysis tools. For lower levels of ground motion, these tools should produce reasonable in-structure response values for evaluation of existing and new facilities. For larger levels of ground motion these tools likely overestimate the in-structure response (and therefore structural demand) since they do not consider geometric nonlinearities (such as gaping and sliding between the soil and structure) and are limited in the ability to model nonlinear soil behavior. The current equivalent linear SSI (SASSI) analysis approach either joins the soil and structure together in both tension and compression or releases the soil from the structure for both tension and compression. It also makes linear approximations for material nonlinearities and generalizes energy absorption with viscous damping. This produces the potential for inaccurately establishing where the structural concerns exist and/or inaccurately establishing the amplitude of the in-structure responses. Seismic hazard curves at nuclear facilities have continued to increase over the years as more information has been developed on seismic sources (i.e. faults), additional information gathered on seismic events, and additional research performed to determine local site effects. Seismic hazard curves are used to develop design basis earthquakes (DBE) that are used to evaluate nuclear facility response. As the seismic hazard curves increase, the input ground motions (DBE’s) used to numerically evaluation nuclear facility response increase causing larger in-structure response. As ground motions increase so does the importance of including nonlinear effects in numerical SSI models. To include material nonlinearity in the soil and geometric nonlinearity using contact (gaping and sliding) it is necessary to develop a nonlinear time domain methodology. This methodology will be known as, NonLinear Soil-Structure Interaction (NLSSI). In general NLSSI analysis should provide a more accurate representation of the seismic demands on nuclear facilities their systems and components. INL, in collaboration with a Nuclear Power Plant Vender (NPP-V), will develop a generic Nuclear Power Plant (NPP) structural design to be used in development of the methodology and for comparison with SASSI. This generic NPP design has been evaluated for the INL soil site because of the ease of access and quality of the site specific data. It is now being evaluated for a second site at Vogtle which is located approximately 15 miles East-Northeast of Waynesboro, Georgia and adjacent to Savanna River. The Vogtle site consists of many soil layers spanning down to a depth of 1058 feet. The reason that two soil sites are chosen is to demonstrate the methodology across multiple soil sites. The project will drive the models (soil and structure) using successively increasing acceleration time histories with amplitudes. The models will be run in time domain codes such as ABAQUS, LS-DYNA, and/or ESSI and compared with the same models run in SASSI. The project is focused on developing and documenting a method for performing time domain, non-linear seismic soil structure interaction (SSI) analysis. Development of this method will provide the Department of Energy (DOE) and industry with another tool to perform seismic SSI analysis.
Executive Summary The overall objective of this work was the development of surface and borehole seismic methodologies using both compressional and shear waves for characterizing faults and fractures in Enhanced Geothermal Systems. We used both surface seismic and vertical seismic profile (VSP) methods. We adapted these methods to the unique conditions encountered in Enhanced Geothermal Systems (EGS) creation. These conditions include geological environments with volcanic cover, highly altered rocks, severe structure, extreme near surface velocity contrasts and lack of distinct velocity contrasts at depth. One of the objectives was the development of methods for identifying more appropriate seismic acquisition parameters for overcoming problems associated with these geological factors. Because temperatures up to 300º C are often encountered in these systems, another objective was the testing of VSP borehole tools capable of operating at depths in excess of 1,000 m and at temperatures in excess of 200º C. A final objective was the development of new processing and interpretation techniques based on scattering and time-frequency analysis, as well as the application of modern seismic migration imaging algorithms to seismic data acquired over geothermal areas. The use of surface seismic reflection data at Brady's Hot Springs was found useful in building a geological model, but only when combined with other extensive geological and geophysical data. The use of fine source and geophone spacing was critical in producing useful images. The surface seismic reflection data gave no information about the internal structure (extent, thickness and filling) of faults and fractures, and modeling suggests that they are unlikely to do so. Time-frequency analysis was applied to these data, but was not found to be significantly useful in their interpretation. Modeling does indicate that VSP and other seismic methods with sensors located at depth in wells will be the most effective seismic tools for getting information on the internal structure of faults and fractures in support of fluid flow pathway management and EGS treatment. Scattered events similar to those expected from faults and fractures are seen in the VSP reported here. Unfortunately, the source offset and well depth coverage do not allow for detailed analysis of these events. This limited coverage also precluded the use of advanced migration and imaging algorithms. More extensive acquisition is needed to support fault and fracture characterization in the geothermal reservoir at Brady's Hot Springs. The VSP was effective in generating interval velocity estimates over the depths covered by the array. Upgoing reflection events are also visible in the VSP results at locations corresponding to reflection events in the surface seismic. Overall, the high temperature rated fiber optic sensors used in the VSP produced useful results. Modeling has been found useful in the interpretation of both surface reflection seismic and VSP data. It has helped identify possible near surface scattering in the surface seismic data. It has highlighted potential scattering events from deeper faults in the VSP data. Inclusion of more detailed fault and fracture specific stiffness parameters are needed to fully interpret fault and fracture scattered events for flow properties (Pyrak-Nolte and Morris, 2000, Zhu and Snieder, 2002). Shear wave methods were applied in both the surface seismic reflection and VSP work. They were not found to be effective in the Brady's Hot Springs area. This was due to the extreme attenuation of shear waves in the near surface at Brady's. This does not imply that they will be ineffective in general. In geothermal areas where good shear waves can be recorded, modeling suggests they should be very useful for characterizing faults and fractures.