Anisotropic seismic-waveform inversion: Application to a seismic velocity model from Eleven-Mile Canyon in Nevada
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
A system and method of characterizing properties of a medium from a non-linear interaction are include generating, by first and second acoustic sources disposed on a surface of the medium on a first line, first and second acoustic waves. The first and second acoustic sources are controllable such that trajectories of the first and second acoustic waves intersect in a mixing zone within the medium. The method further includes receiving, by a receiver positioned in a plane containing the first and second acoustic sources, a third acoustic wave generated by a non-linear mixing process from the first and second acoustic waves in the mixing zone; and creating a first two-dimensional image of non-linear properties or a first ratio of compressional velocity and shear velocity, or both, of the medium in a first plane generally perpendicular to the surface and containing the first line, based on the received third acoustic wave.
Abstract not provided.
Presentation to be presented during the Site-Directed Research and Development (SDRD) program FY 2021 review meeting (Webex), September 22–23, 2021.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Seismic isolation is a mature technology that has been implemented in many civil structures, including buildings, bridges, liquid natural gas tanks, and off shore oil platforms, both in the United States and other countries, to mitigate the damaging effects of earthquakes. Seismic isolation has also been used in nuclear structures in France and South Africa, but not yet in the United States: neither in the Department of Energy facilities nor in commercial nuclear power plants (NPPs). This is primarily due to a lack of guidelines, codes and standards for the analysis, design and construction specific to seismically isolated nuclear structures. However, seismic isolation of nuclear structures has seen increased research interest in the recent years and the recently published national consensus standard, America Society of Civil Engineers (ASCE) Standard 4-16 'Seismic analysis of safety related nuclear structures', now incorporates language and commentary (Chapter 12) for seismically isolating surface or near-surface-mounted nuclear facilities, including NPPs. Seismic isolation substantially reduces horizontal seismic loads (demands) on structures, systems, and components. Reduction in demand results in four potential benefits: (1) economic: reduction in capital cost, (2) increased safety: reduction in the mean annual frequency of unacceptable performance, (3) insurance: protection against increases in the known seismic hazard after construction by minimizing the effort to re-qualify and re-certify structures, systems and components, and (4) recertification: the opportunity to certify an existing NPP design for a region of higher seismic hazard. Item (2) above, an increase in safety (reduction of seismic risk), was explored in Huang et al., wherein it was demonstrated that the implementation of seismic isolation reduced seismic risk in nuclear power plant structures. Studies that assess item (1) have not been performed although the use of isolation may lead to large reductions in the capital cost of safety-related nuclear facilities. Funding provided by Nuclear Safety Research and Development in the Department of Energy (DOE) allowed the authors to develop a framework for assessing the economic benefits and reductions in seismic risk afforded by the use of seismic (base) isolation. The framework includes probabilistic risk assessment and estimation of overnight capital costs (OCC's) of a sample generic DOE nuclear facility (GDNF). The project is documented in detail in Bolisetti et al. and Yu et al., and key results are presented here. The GDNF considered in this study is based on a facility constructed at the Idaho National Laboratory. The GDNF is a two-story reinforced concrete building with an embedded basemat and is founded on soft rock at the boundary of Site Class B and Site Class C (B/C boundary) as characterized by ASCE/SEI Standard 7. The GDNF is assumed to handle radiological materials, and its structure, systems and components (SSC's) are designed to effectively confine these materials in the event of an internal failure. The GDNF is populated with components that are generic to safety-related nuclear structures. A single system, which is intended to confine a loss of materials-at-risk (MAR) to the interior of the structure, is chosen for this study. Accordingly, a hypothetical event tree and fault tree are assumed for risk assessment. The isolation system for the GDNF is designed per Chapter 12 of ASCE/SEI 4-16 but is not optimized. The isolation system consists of 38 lead rubber bearings installed in one horizontal plane below the basemat. Details of the GDNF, including the SSC's, event tree and fault tree, and isolation system, are presented in the project report and Yu et al.. This study includes seismic risk assessment and cost estimation for four cases: (1) the GDNF is located at a site of low to moderate seismic hazard and is constructed on a conventional (non-isolated) foundation; (2) the GDNF is located at the same site as case 1, but is seismically isolated; (3) the GDNF is located at a site of high seismic hazard and is constructed on a conventional foundation; and (4) the GDNF is located at the same site of high seismic hazard as case 3, but is seismically isolated. In this study, the DOE site at Idaho National Laboratory (INL) is chosen as the low to moderate seismic hazard site for cases 1 and 2, and the Los Alamos National Laboratory (LANL) is chosen as the high seismic hazard site for cases 3 and 4.
Time-lapse surface seismic surveys have been widely used at carbon sequestration sites for site characterization, monitoring subsurface CO 2 plume migration, and detecting potential CO 2 leakage from a storage reservoir. Monitoring in the first permeable unit directly above the primary seal is important for early detection of CO 2 leakage. Forward modeling of time-lapse seismic data can be used to assess the utility of the seismic method for CO 2 leakage detection. We develop a workflow for forward modeling of time-lapse seismic data, including constructing seismic velocity models using flow simulation outputs, modeling of pre-stack and post-stack synthetic seismic data following seismic data processing sequence and analysis of processed synthetic time-lapse seismic data. We apply the forward modeling and analysis workflow to assessing the detectability and the earliest detection time of seismic monitoring using the hypothetical CO 2 leakage scenarios for a model geologic storage site with realistic geology. We then derive the detection thresholds using the simulated normalized root-mean-square (NRMS) differences for the no-leakage case at a range of signal-to-noise ratios, representing the background noise levels in seismic data. We then compare NRMS differences triggered by the CO 2 leakage to the detection thresholds at each time step to quantify the detectability and the earliest detection time of seismic monitoring. We analyze the effects of the acquisition parameters and elastic parameters on the produced synthetic seismic data and earliest detection time. Our modeling results indicate that high signal-to-noise ratio is needed to detect the CO 2 leakage at the model site. Minimizing the background noise in seismic data is crucial for improving the detectability of the seismic method. Increasing the shot density or increasing the dominant frequency of the source wavelet is likely to increase the possibility of the leakage detection and reduce the time needed for the detection. The elastic parameters used in the rock physics modeling have significant effects on the resultant seismic velocity models and synthetic seismic data, highlighting their critical roles in understanding and interpreting time-lapse seismic reflection data associated with CO 2 monitoring, verification and accounting activities.
Seismic techniques are the dominant geophysical techniques for the characterization of subsurface structures and stratigraphy. The seismic techniques also dominate the monitoring and mapping of reservoir injection and production processes. Borehole seismology, of all the seismic techniques, despite its current shortcomings, has been shown to provide the highest resolution characterization and most precise monitoring results because it generates higher signal to noise ratio and higher frequency data than surface seismic techniques. The operational environments for borehole seismic instruments are however much more demanding than for surface seismic instruments making both the instruments and the installation much more expensive. The current state-of-the-art borehole seismic instruments have not been robust enough for long term monitoring compounding the problems with expensive instruments and installations. Furthermore, they have also not been able to record the large bandwidth data available in boreholes or having the sensitivity allowing them to record small high frequency micro seismic events with high vector fidelity. To reliably achieve high resolution characterization and long term monitoring of Enhanced Geothermal Systems (EGS) sites a new generation of borehole seismic instruments must therefore be developed and deployed. To address the critical site characterization and monitoring needs for EGS programs, US Department of Energy (DOE) funded Paulsson, Inc. in 2010 to develop a fiber optic based ultra-large bandwidth clamped borehole seismic vector array capable of deploying up to one thousand 3C sensor pods suitable for deployment into ultra-high temperature and high pressure boreholes. Tests of the fiber optic seismic vector sensors developed on the DOE funding have shown that the new borehole seismic sensor technology is capable of generating outstanding high vector fidelity data with extremely large bandwidth: 0.01 – 6,000 Hz. Field tests have shown that the system can record events at magnitudes much smaller than M-2.6 at frequencies up to 2,000 Hz. The sensors have also proved to be about 100 times more sensitive than the regular coil geophones that are used in borehole seismic systems today. The fiber optic seismic sensors have furthermore been qualified to operate at temperatures over 300°C (572°F). Simultaneously with the fiber optic based seismic 3C vector sensors we are using the lead-in fiber to acquire Distributed Acoustic Sensor (DAS) data from the surface to the bottom of the vector array. While the DAS data is of much lower quality than the vector sensor data it provides a 1 m spatial sampling of the downgoing wavefield which will be used to build the high resolution velocity model which is an essential component in high resolution imaging and monitoring.
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.
We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.
The passive seismic network deployed on the Moon during the Apollo missions operated for eight years and allowed the observation of seismic activity. More than 12500 seismic events were registered, where 28 were classified as shallow moonquakes with moment magnitudes up to 4.1. Seismic events of this nature pose a significant risk to future long-term lunar habitats; thus, these events must be carefully studied and considered in the seismic design of these structures. This paper proposes a preliminary seismic hazard assessment imposed by shallow moonquakes. The hazard assessment is performed using the Probabilistic Seismic Hazard Analysis (PSHA) methodology, considering previous studies and theories regarding the seismic environment of the Moon. The study zone covers ~860 km2 of the Taurus-Littrow Valley, containing the Apollo 17 landing site and the Lee-Lincoln lobate scarp as the considered seismic source. The seismic hazard is quantified in terms of peak ground acceleration (PGA) and spectral acceleration (5% damped pseudo-acceleration, PSA). Seismic hazard deaggregation scenarios, Uniform Hazard Spectra (UHS) for different hazard levels, and a Conditional Mean Spectrum (CMS) for a target period of 0.2 s are obtained to quantify the seismic hazard on a specific site on the Moon. The developed seismic hazard assessment provides a preliminary approach for realistic scenarios to conduct structural designs that ensure the seismic performance of fully operational long-term lunar structures.
The present ]day internal structure of the Moon provides insight not only into its own formation and evolution, but also that of all rocky planetary bodies. The most direct way to probe a planet fs interior structure is through seismology. As part of the Apollo lunar missions, four seismometers were deployed on the nearside surface of the Moon between the years 1969 and 1972. These instruments operated continuously until 1977, forming the only substantial extraterrestrial seismic data set in existence. These data have been used to constrain various aspects of the seismic velocity and density structure of the Moon. Typical 1-D models recognize a 30-60 km thick crust overlying a nearly constant ]velocity mantle, and extend to a depth of approximately 1000 km, below which the lack of penetrating moonquake ray ]paths precludes the seismic determination of deeper structure. Previously, the lack of observed moonquakes from the far side of the Moon has been used to infer the presence of a highly attenuating (possibly molten) core. Indirect geophysical measurements such as moment of inertia, magnetic induction, lunar laser ranging, and elemental abundances of mare basalts also place varying constraints on core size and state. In combination with seismic studies, these indirect measurements have been used to arrive at a commonly accepted model of the Moon's deepest interior that includes a solid inner and fluid outer core, overlain by a partial melt boundary layer. We recently applied modern array seismology techniques to the Apollo data and revealed detailed core structure, including the first direct confirmation of the presence of a solid inner core. Our study focused on the identification of core ]reflected phases in deep moonquake seismograms. The resulting model of the Moon fs innermost structure was found to be consistent with the commonly accepted model. However, the modeled layer radii may vary by tens of kilometers, as is expected when accounting for uncertainties such as moonquake location, timing errors, and potential seismic heterogeneities. In addition, the modeled velocities may vary with a 1-to-1 trade ]off with the modeled reflector depth. The GRAIL (Gravity Recovery and Interior Laboratory) mission, launched in Sept. 2011, placed two nearly identical spacecraft in lunar orbit. The two satellites make extremely high-resolution measurements of the lunar gravity field, which can be used to constrain the interior structure of the Moon using a "crust to core" approach. GRAIL fs constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten) will complement seismic investigations in several ways. Here we present a progress report on our efforts to advance our knowledge of the Moon fs internal structure using joint gravity and seismic analyses. We will focus on methodology, including 1) refinements to the seismic core constraint accomplished through array processing of Apollo seismic data, made by applying a set of travel time corrections based on GRAIL structure estimates local to each Apollo seismic station; 2) modeling deep lunar structure through synthetic seismograms, to test whether the seismic core model can reproduce the core reflections observed in the Apollo seismograms; and 3) a joint seismic and gravity inversion in which we attempt to fit a family of seismic structure models with the gravity constraints from GRAIL, resulting in maps of seismic velocities and densities that vary from a nominal model both laterally and with depth.
Geological fault detection and characterization are crucial for understanding subsurface dynamics across scales. While methods for fault delineation based on either seismicity location analysis or seismic image reflector discontinuity are well-established, a systematic approach that integrates both data types remains absent. We develop a novel machine learning model that unifies seismic reflector images and seismicity location information to automatically identify geological faults and characterize their geometrical properties. The model encodes a seismic image and a seismicity location image separately, and fuses the encoded features with a spatial-channel attention fusion module to improve the learning of important features in both inputs. We design an automated strategy to generate high-quality synthetic training data and labels. To improve the realism of the seismicity location image, we include random seismicity noise and missing seismicity location associated with some of the faults. We validate the model’s efficacy and accuracy using synthetic data examples and two field data examples. Moreover, we show that fine-tuning the trained model with a small, domain-specific dataset enhances its fidelity for field data applications. The results demonstrate that integrating seismicity location and seismic images into a unified framework allows the end-to-end neural network to achieve higher fidelity and accuracy in delineating subsurface faults and their geometrical properties compared with image-only fault detection methods. Our approach offers an adaptive data-driven tool for geological fault characterization and seismic hazard mitigation, bridging the gap between seismicity location and image-based fault detection methods.
Risk calculations should focus on providing best estimate results, and associated insights, for evaluation and decision-making. Specifically, seismic probabilistic risk assessments (SPRAs) are intended to provide best estimates of the various combinations of structural and equipment failures that can lead to a seismic induced core damage event. However, in some instances the current SPRA approach has large uncertainties, and potentially masks other important events (for instance, it was not the seismic motions that caused the Fukushima core melt events, but the tsunami ingress into the facility). SPRA’s are performed by convolving the seismic hazard (this is the estimate of all likely damaging earthquakes at the site of interest) with the seismic fragility (the conditional probability of failure of a structure, system, or component given the occurrence of earthquake ground motion). In this calculation, there are three main pieces to seismic risk quantification, 1) seismic hazard and nuclear power plants (NPPs) response to the hazard, 2) fragility or capacity of structures, systems and components (SSC), and 3) systems analysis. Two areas where NLSSI effects may be important in SPRA calculations are, 1) when calculating in-structure response at the area of interest, and 2) calculation of seismic fragilities (current fragility calculations assume a lognormal distribution for probability of failure of components). Some important effects when using NLSSI in the SPRA calculation process include, 1) gapping and sliding, 2) inclined seismic waves coupled with gapping and sliding of foundations atop soil, 3) inclined seismic waves coupled with gapping and sliding of deeply embedded structures, 4) soil dilatancy, 5) soil liquefaction, 6) surface waves, 7) buoyancy, 8) concrete cracking and 9) seismic isolation The focus of the research task presented here-in is on implementation of NLSSI into the SPRA calculation process when calculating in-structure response at the area of interest. The specific nonlinear soil behavior included in the NLSSI calculation presented in this report is gapping and sliding. Other NLSSI effects are not included in the calculation. The results presented in this report document initial model runs in the linear and nonlinear analysis process. Final comparisons between traditional and advanced SPRA will be presented in the September 30th deliverable.
The primary goal of the project was met, as improved resolution (from kilometers to tens/hundreds of meters and to a depth of at least 3 km) was added to the seismic interferometry-estimated parameters relative to the work reported in AltaRock Energy Inc. (2014) for the Dixie Valley Geothermal Wellfield. A procedure for initial fault identification, using seismic interferometry, was developed. The project team further improved, and applied a new, cost-effective, and non-invasive exploration method, developed and validated since 2010, during two previous experiments, at Dixie Valley and Soda Lake, Nevada. The method, initially designed for hydrothermal exploration, was adapted to EGS exploration conditions, which require deeper crustal information. Regions of elevated energy for both the 2-7 Hz and 5-15 Hz frequency bands are located above intersections of multiple faults. Thus, an unprecedented, innovative seismic field deployment and ambient seismic noise and signal processing were successfully used to estimate crustal structure and to identify fault locations at depths up to 3 km, at a cost per kilometer of survey that in some cases may be as much as ten times less than the cost of a controlled-source survey. Based on valid and useful seismic data collected and analyzed during Phase I, in Phase II the team has further quantified statistically the advantage of using seismic analysis, in combination with other geological and geophysical techniques, for EGS target identification. A power spectral density analysis of the passive-seismic empirical Green’s functions reveals significantly higher power recorded over the exploration area at the power plant, and in the drilling area, near the BB station 10201, especially on the horizontal components. In this study, the most important relationship observed during the prior report, namely the one between velocity and temperature was explored. The areas that show anomalously high predicted temperature in the geostatistics occur where multiple faults intersect. Not surprisingly, such areas span between the Dixie Valley range front fault and the basin-ward Piedmont faults, in the vicinity of numerous wells. These locations also fall in the area of high EGS favorability from the prior 2014 study. The anomalies also coincide with the high power spectral density values. The result of altering all favorability criteria for the new data reached a maximum of 1/8th of a point and were not large enough to change any favorability value reported in the 2014 study. The new data reinforce the results obtained from the previous 2014 study. The current high-resolution ambient seismic noise project provided excellent results in defining faults in the wellfield. As such the identification of known faults by the project demonstrated the utility of this method for identifying faults at an exploration stage in geothermal development. Such fault identification, while extremely useful in exploration investigation, has not added sufficient new data for the generation of significantly different DVGW favorability maps in the current project. We believe the cost savings of our passive-seismic technique for EGS, over common industry active-source survey techniques, may have a minimum of 30%, and a maximum of 90%. There are many significant non-cost benefits of passive seismic surveys. Our passive-source method provides an environmentally-friendly seismic survey with minimal to negligible impacts. The economic benefit is greater EGS project feasibility.
We report nuclear energy has a key role to play in global decarbonization. Impediments to the widespread deployment of reactors are their projected high capital cost and levelized cost of energy, and time required to analyze, design, license, construct, and commission them. The earthquake load case is a key cost driver for a new build nuclear plant, because near-surface soils and seismic hazard are different at each site, requiring site-specific analysis, design, engineering, qualification, licensing, and regulatory review, essentially making every design First-of-a-Kind (FoaK). To enable deployment at the scale needed for deep decarbonization, the cost and time impact of the seismic load case must be significantly mitigated, and plants must be standardized. Seismic base isolation has been proven to considerably reduce the earthquake response of structures and equipment but has yet to be applied to a nuclear power plant in the United States, in part because the financial impacts, positive or negative, are not known. Because there are no recent non-proprietary data to characterize the influence of the seismic load case on capital cost, it is difficult to confidently quantify the financial benefits of seismic isolation. Scheme-level designs of two fundamentally different advanced reactor buildings were developed to assemble cost data on the influence of the seismic load case. Both buildings were equipped with three bespoke pieces of safety-related equipment and analyzed for incremented levels of earthquake shaking to quantify the seismic penalty on equipment, in terms of vessel weights and horizontal accelerations. Using analysis results, a questionnaire was developed and transmitted to nuclear utilities, reactor developers, engineers, and equipment suppliers to collect cost data on engineering and fabrication costs for these unique pieces of safety-class equipment. Synthesis of the cost data showed that the seismic load case significantly affects the capital cost (sum of engineering and fabrication cost) of safety-class equipment, with engineering costs being comparable to fabrication costs. Standardization of safety-class equipment is made possible by seismic isolation, that is, equipment designed for minimal seismic robustness can resist earthquake shaking at a site of much higher seismic hazard. The average reduction in the capital cost of the safety-related equipment, enabled by seismic isolation, is a factor of two for FoaK equipment and a factor of five for standardized equipment.
Abstract The carbon capture, utilization, and storage (CCUS) framework is an essential component in reducing greenhouse gas emissions, with its success hinging on the comprehensive knowledge of subsurface geology and geomechanics. Passive seismic event relocation and fault detection offer vital insights into subsurface structures and the ability to monitor fluid migration pathways. Accurate identification and localization of seismic events, however, face significant challenges, including the necessity for high‐quality seismic data and advanced computational methods. To address these challenges, we introduce a novel deep learning method, , specifically designed for passive seismic source relocation and fault delineating for passive seismic monitoring projects. By leveraging data domain‐adaptation, allows us to train a neural network with labeled synthetic data and apply it directly to field data. Using , the passive seismic sources are automatically clustered based on their recording time and spatial locations, and subsequently, faults and fractures are delineated accordingly. We demonstrate the efficacy of on a field case study involving injection related microseismic data from Decatur, Illinois area. Our approach accurately and efficiently relocated passive seismic events, identified faults and could aid in potential damage induced by seismicity. Our results highlight the potential of as a valuable tool for passive seismic monitoring, emphasizing its role in ensuring CCUS project safety. This research bolsters the understanding of subsurface characterization in CCUS, illustrating machine learning’s capacity to refine these methods. Ultimately, our work has significant implications for CCUS technology deployment, an essential strategy in combating climate change. Plain Language Summary In our quest to tackle climate change, we use a strategy known as carbon capture, utilization, and storage (CCUS) to keep greenhouse gases out of the atmosphere. This strategy relies heavily on our ability to understand what's happening deep under the earth's surface. To make sure we store super critical safely, we need to accurately map out the geological structure, especially faults, but this is tough without high‐quality data and complex computer programs. We've developed a new tool called “DeFault,” which uses advanced machine learning to improve how we find and map these underground features. “DeFault” is smart enough to learn from numerically simulated data and then apply what it’s learned to real‐world situations. It groups together seismic activity—tiny tremors and shifts in the earth—based on when and where they happen, which helps us spot where there might be cracks or faults. We tested “DeFault” in Illinois, where CO 2 is injected underground, and it successfully pinpointed where these tremors occurred and mapped out the faults, helping to prevent accidents accurately in the future. Our study shows that “DeFault” will be a powerful ally in making CCUS safer and more effective, especially for the Illinois Basin Decatur Project. Key Points Faults and fractures introduced by carbon storage can be monitored by passive seismicity DeFault algorithm enables an automatic process for accurate and efficient passive seismic event locating and clustering