Search NASA⌕ Search

SEARCH · Search NASA

Results for “SEISMICITY”

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

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

At least 163 records · Page 9

Seismic Assessment of a Long-Term Lunar Habitat

The establishment of secure earth-independent long-term lunar habitats has been envisioned by numerous government agencies and private companies. Recent advancements in assessing seismic hazards caused by shallow moonquakes have highlighted the importance of incorporating this phenomenon into the design of robust and resilient lunar structures. However, further research is required to explore lunar habitat design that considers seismic loads. This paper proposes assessing the structural response of a lunar habitat made of sulfur concrete covered with a regolith layer. The numerical model of the structure is subjected to gravitational, internal pressure and seismic loads. The seismic analysis of the structure is carried out using spectral and nonlinear time history methods. Conditional mean spectra for shallow moonquakes with return periods of 75, 475, 970 and 2475 years are used in the seismic analysis. The records used for the temporal analyses were ground motions that agree with a preliminary seismic hazard on the Moon. The results of temporal analyses reveal that shallow moonquakes with return periods greater than 475 years can lead to the loss of the global stability of the structure. Consequently, the findings imply that seismic loads have the potential to impose unacceptable demands on lunar structures constructed from in-situ materials like sulfur concrete. Hence, it is imperative to incorporate seismic considerations in the design process for developing resilient and long-term lunar habitats.

Lunar habitat↗

Joint impedance and facies inversion of time-lapse seismic data for improving monitoring of CO 2 incidentally stored from CO 2 EOR

Time-lapse seismic monitoring is an effective and proven technology for mapping the distribution of CO 2 in a subsurface reservoir. When injected CO 2 displaces other reservoir fluids, porous-medium properties are changed and thus the seismic impedance changes, causing time-lapse seismic amplitude differences in the injection zones. The analysis and interpretation of images created from these amplitude differences can provide information about reservoir architecture and the CO 2 migration within the reservoir. Incorporating seismic inversion and rock physics into the interpretation of time-lapse seismic data can considerably improve the modeling and monitoring to detect and assess the location of CO 2 over time. The joint inversion method presented in this paper has an integral representation of the geology in the inversion algorithm using elastic facies, which provides information about the spatial distribution of the geologic heterogeneities controlling the movement of fluids in the reservoir. The method was successfully applied to time-lapse seismic data from a mature oil field undergoing CO 2 enhanced oil recovery. The estimated seismic acoustic impedances and facies reflect the characteristics of individual geologic facies and fluid conditions of the reservoir subject to CO 2 injection. Finally, the probabilities estimated by the joint impedance and facies inversion for the reservoir's litho-fluid facies can be used for forecasting CO 2 saturation and pressure changes within the target reservoir.

4D seismic↗

Aseismic Motions Drive a Sparse Seismicity During Fluid Injections Into a Fractured Zone in a Carbonate Reservoir

An increase in fluid pressure in faults can trigger seismicity and large aseismic motions. Understanding how fluid and faults interact is an essential goal for seismic hazard and reservoir monitoring, but this key relation remains unclear. We developed an in situ experiment of fluid injections at a 10 meter scale. Water was injected at high pressure in different geological structures inside a fault damaged zone, in limestone at 280 m depth in the Low Noise Underground Laboratory (France). Induced seismicity, as well as strains, pressure, and flow rate, was continuously monitored during the injections. Although nonreversible deformations related to fracture reactivations were observed for all injections, only a few tests generated seismicity. Events are characterized by a 0.5-to-4 kHz content and a small magnitude (approximately -3.5). They are located within 1.5 m accuracy between 1 and 12 m from the injections. Comparing strain measurements and seismicity shows that more than 96% of the deformation is aseismic. The seismic moment is also small compared to the one expected from the injected volume. Moreover, a dual seismic behavior is observed as (1) the spatiotemporal distribution of some cluster of events is clearly independent from the fluid diffusion (2) while a diffusion-type pattern can be observed for some others clusters. Furthermore, the seismicity might therefore appear as an indirect effect to the fluid pressure, driven by aseismic motion and related stress perturbation transferred through failure.

58 GEOSCIENCES↗

Seismic Data Gathering and Validation

Three recent earthquakes in the last seven years have exceeded their design basis earthquake values (so it is implied that damage to SSC’s should have occurred). These seismic events were recorded at North Anna (August 2011, detailed information provided in [Virginia Electric and Power Company Memo]), Fukushima Daichii and Daini (March 2011 [TEPCO 1]), and Kaswazaki-Kariwa (2007, [TEPCO 2]). However, seismic walk downs at some of these plants indicate that very little damage occurred to safety class systems and components due to the seismic motion. This report presents seismic data gathered for two of the three events mentioned above and recommends a path for using that data for two purposes. One purpose is to determine what margins exist in current industry standard seismic soil-structure interaction (SSI) tools. The second purpose is the use the data to validated seismic site response tools and SSI tools. The gathered data represents free field soil and in-structure acceleration time histories data. Gathered data also includes elastic and dynamic soil properties and structural drawings. Gathering data and comparing with existing models has potential to identify areas of uncertainty that should be removed from current seismic analysis and SPRA approaches. Removing uncertainty (to the extent possible) from SPRA’s will allow NPP owners to make decisions on where to reduce risk. Once a realistic understanding of seismic response is established for a nuclear power plant (NPP) then decisions on needed protective measures, such as SI, can be made.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗

Geothermal Induced Seismicity National Environmental Policy Act Review

In 2016, the U.S. Bureau of Land Management (BLM) contracted with the National Renewable Energy Laboratory (NREL) to assist the BLM in developing and building upon tools to better understand and evaluate induced seismicity caused by geothermal projects. This review of NEPA documents for four geothermal injection or EGS projects reveals the variety of approaches to analyzing and mitigating induced seismicity. With the exception of the Geysers, where induced seismicity has been observed and monitored for an extended period of time due to large volumes of water being piped in to recharge the hydrothermal reservoir, induced seismicity caused by geothermal projects is a relative new area of study. As this review highlights, determining the level of mitigation required for induced seismic events has varied based on project location, when the review took place, whether the project utilized the International Energy Agency or DOE IS protocols, and the federal agency conducting the review. While the NEPA reviews were relatively consistent for seismic monitoring and historical evaluation of seismic events near the project location, the requirements for public outreach and mitigation for induced seismic events once stimulation has begun varied considerably between the four projects. Not all of the projects were required to notify specific community groups or local government entities before beginning the project, and only one of the reviews specifically stated the project proponent would hold meetings with the public to answer questions or address concerns.

15 GEOTHERMAL ENERGY↗

Time-lapse seismic inversion for CO 2 saturation with SeisCO2Net: An application to Frio-II site

Seismic monitoring of geological CO 2 storage (GCS) involves highly nonlinear seismic inversion and petrophysical inversion, making it challenging to estimate CO 2 volume efficiently and detect possible early CO 2 leakages. Deep learning (DL) using convolutional neural networks (CNNs) has shown promise in solving highly nonlinear seismic inversion problems. However, direct estimation of CO 2 plume extent/saturation from time-lapse seismic gathers using DL is still underexplored, with no reported field applications to date. The investigation of field data is primarily hindered by scarcity of field data for neural network training. Other obstacles include highly nonlinear seismic-petrophysics inverse relationship, and presence of noise in field seismic data. We introduce SeisCO2Net, a deep CNN that predicts CO 2 saturation maps directly from time-lapse full waveform shot gathers. For training, we use site-specific geological information, fluid flow physics, rock physics, and seismic modeling to generate synthetic datasets that closely resemble the CO 2 storage site. Synthetic tests show promising results, inspiring us to apply SeisCO2Net's trained weights on field data collected at Frio-II GCS site by leveraging transfer learning principles. As reference, we compare SeisCO2Net's predicted CO 2 saturation maps with results obtained from physics-based inversion. Our analyses show both methods display similar CO 2 plume shapes, reasonable CO 2 plume characteristics, and comparable saturation values. Our results suggest pre-training CNNs on physics-informed synthetic datasets and then applying the learned weights to field data is a viable approach to estimating field CO 2 saturation. This method effectively addresses the scarcity of field training data, thus encouraging the feasibility of long-term GCS monitoring.

58 GEOSCIENCES↗

Projected Seismic Activity at the Tiger Stripe Fractures on Enceladus, Saturn, From an Analog Study of Tidally Modulated Icequakes Within the Ross Ice Shelf, Antarctica

Abstract Dissipation of tidal energy is expected to generate seismicity on icy‐ocean worlds; however, the distribution and timing of this seismic activity throughout an orbital cycle is not known. We used new observations from an icy‐ocean‐world analog environment on Earth to examine the relationship between tidally driven tensile stress and seismic activity within an ice shell. We investigated a pair of rifts within Antarctica's Ross Ice Shelf which are tidally stressed in a manner analogous to the orbital cycle of tidal stress experienced by Enceladus' Tiger Stripe Fractures. We found that seismic activity at the Antarctic rifts is sensitive to both the amplitude and the rate of tensile stress across the rifts. We combined these findings with calculated stress values along Enceladus' Tiger Stripe Fractures to predict seismic‐activity levels expected along the ice‐shell fractures. We predict a peak in seismicity along the four Tiger Stripe Fractures when Enceladus is 90°–120° past pericenter in its orbit around Saturn, at which point tensile stresses would reach ∼2/3 of their maximum value. We also used the magnitude distribution of icequakes along Antarctic rifts to investigate implications for the likely size of stick‐slip rupture patches along icy faults on Enceladus. Our findings predict that the Tiger Stripe Fractures should produce sustained, low‐magnitude seismic events that involve rupture along discrete portions of each fracture's total length. We predict that seismicity would fall to 50% of peak levels when stresses across the Tiger Stripe Fractures are dominantly compressional.

Olsen, Kira G.↗

Pore Pressure Diffusion and Onset of Induced Seismicity

The Raton Basin of Southern Colorado and Northern New Mexico experienced a significant increase in seismicity over the last several decades, including two M5+ earthquakes. The increase in seismicity started soon after an increase in wastewater disposal associated with coal-bed methane production. Twenty-nine wastewater injection wells have operated in the Basin, with several high-rate injectors located within a kilometer of the most seismically active regions. We developed a pore pressure model to investigate the relationship between pore pressure diffusion and onset of induced seismicity, with an emphasis on the Trinidad fault zones. Our model indicates that sufficient pore pressure increase had accumulated on the northern extent of the Trinidad fault zone to induce seismicity at the time of early-stage seismic activity. The hydrogeologic architecture of the Basin played a key role in controlling pore pressure perturbation. The basal sedimentary interval with moderately low diffusivity restricted vertical pore pressure diffusion, creating a substantial lag time between wells becoming active and the expression of perturbation in the fault zones. The spatiotemporal relationship between early-stage seismicity and wastewater injection in conjunction with modeled pore pressures suggests that the early-stage seismicity in the Raton Basin was likely triggered by pore pressure diffusion.

58 GEOSCIENCES↗

Insights Into Seismicity Associated With Flexibly Operating Enhanced Geothermal System From Real‐Time Distributed Acoustic Sensing

Enhanced Geothermal Systems (EGS) have the capacity to broaden the accessible resource pool for geothermal power generation. Traditionally viewed as a “baseload” resource, their flexible operation might also enable dispatchable load‐following generation and long‐term energy storage, aligning them with the evolving landscape of decarbonized electricity systems. However, increasing permeability and extracting energy during EGS operations can induce microseismic events; for many prior EGS efforts, some associated seismicity has been observed. While energetically beneficial, the flexibility of EGS operations prompts our inquiry into whether new types of operations will yield previously unseen seismicity patterns. We demonstrate the use of distributed acoustic sensing (DAS) with real‐time edge computing to monitor seismicity during a pilot test of a cyclically operated EGS facility at the Blue Mountain geothermal field. Our focus lies in uncovering seismicity insights from the real‐time microseismic catalog, particularly during load‐following dispatchability tests simulating flexible EGS operation. Here, we find that variations in pore pressure consistently correlate with seismicity, and that controlling pressure cycles during flexible operations appears to constrain microseismic activity during subsequent cycles. The spatio‐temporal evolution of microseismic clouds recorded during cyclic injection cycles fits diffusive models over our available observation period. Additionally, seismicity elevation lags behind pore pressure increases, likely due to pressure diffusion to the fracture system boundary. Through real‐time monitoring, we offer novel insights into seismicity associated with flexibly operating EGS. Our findings suggest that leveraging DAS and edge computing can inform EGS operations and help mitigate induced seismicity.

Chamarczuk, Michal [Rice Univ., Houston, TX (Unite↗

Iterative multi-task learning and inference from seismic images

Seismic interpretation aims to extract quantitative and interpretable attributes from a seismic image produced using some migration method to inform characteristics of a subsurface reservoir or target of interest. Current paradigms for computing seismic attributes mostly rely on single-task algorithms. We develop an iterative, multi-task machine learning method to learn and infer multiple attributes from a seismic image. This method is composed of two stages: a multi-task inference stage and a multi-modal, multi-task refinement stage. The basic mechanism of this method is that we train a multi-task inference neural network (NN) to estimate a set of attributes, including a relative geological time (RGT), a denoised higher-resolution (DHR) seismic image, and multiple fault attributes (including probability, dip, and strike), from a low-resolution, noisy seismic image; then we input the inferred attributes to a multi-task refinement NN to enhance the raw inference results iteratively. The two multi-task NNs are trained separately based on synthetic seismic images and associated attributes generated by a geological modeling algorithm. The software we intend to release is a PyTorch implementation of this multi-task learning method for both 2D and 3D cases along with scripts to run the training/validation. The algorithm and software can be a useful tool for automatic seismic interpretation.

Gao, Kai↗

Seismic response to paleo-sand dunes in the Nugget Sandstone Formation, southwestern Wyoming

We have analyzed a 3D seismic survey acquired for a carbon sequestration project on top of the Moxa Arch in southwestern Wyoming. We observed a zone of discontinuous reflectors on vertical slices of seismic amplitude volume, whereas, the northwest–southeast lineations were observed on the time slices. We performed a seismic to well tie that suggested that the lineations occur within the Nugget Sandstone. The Nugget Sandstone is an eolian sandstone deposit of Early Jurassic age, deposited as a subtropical dune field, and equivalent to the Navajo Sandstone of southwestern Utah. Petrophysical analysis indicates that the Nugget Sandstone is dominated by clean sandstone (70%–80%), whereas evaporites, including halite and anhydrite, are present in certain zones. Previous outcrop studies on the Navajo Sandstone indicate the wind direction to be predominantly northeast–southwest. Seismic attributes, including coherence and curvature, displayed on stratal slices within the Nugget Sandstone interval indicate the presence of lineations in the northwest–southeast direction with irregular spacing. These lineations are approximately perpendicular to the inferred dominant wind direction. We computed the dominant wind direction from the average azimuth of the lineations as seen on the curvature attribute in the Nugget Sandstone interval. Geological feature: Eolian sand dunes with interdunal evaporites Seismic appearance: Parallel lineations with irregular spacing on seismic attribute horizon slices Alternative interpretations: Canyons at continental slopes; slope failures Features with a similar appearance: Marine bars; contourites Formation: The Nugget Sandstone — equivalent to the Navajo Sandstone Age: Early Jurassic Location: Moxa Arch, Wyoming Seismic data: Obtained by the University of Wyoming with U.S. DOE funding Contributors: Dhruv Agrawal, Brady Lujan, Sumit Verma, Shuvajit Bhattacharya, and Subhashis Mallick Analysis tools: Coherence and curvature attributes; seismic inversion; petrophysical inversion

Geochemistry & Geophysics↗

Spectral Characteristics of Daily to Seasonal Ground Motion at the Piñon Flats Observatory from Coherence of Seismic Data

We investigate coherences of seismic data recorded during three years (2015–2017) at the Piñon Flats Observatory (PY) array and a collocated 148 m deep borehole station B084, along with oceanic data from a buoy southwest of the PY array. Seismic and barometric recordings at PY stations are analyzed with a multitaper spectral technique. The coherence of signals from seismic sources is >0.6 at 0.05–8 Hz between closely spaced (<65m) surface stations and decreases to ~0.2 in frequency bands in which the wavelengths are smaller than interstation distances. There are several local coherence increases at 1–8 Hz between nearby (<65 m) surface stations, whereas large coherence values between a surface and 148 m deep borehole stations are only present at the secondary microseism (~0.14 Hz). These points to significant modification of seismic recordings in the top crust, and those continual near-surface failures might produce shallow rapidly attenuating signals at surface stations. Incoherent local atmospheric effects induce incoherent seismic signals in low- and high-frequency ranges through different coupling mechanisms. Between 0.003 and 0.05 Hz, atmospheric loadings generate ground tilts that contaminate the two horizontal seismic recordings and decrease their coherence, whereas the vertical component is less affected. At 1–8 Hz, coupling of atmospheric pressure with surface structures transmits incoherent signals into the ground, degrading the seismic coherence in all three components. The two horizontal coherences show seasonal variations with extended coherent frequency bands in winter and spring, likely to be produced by seasonal variations in microseisms and local ground tilts. The coherences also contain high anomalies between 2 and 4 Hz resulting from anthropogenic activities. The results provide useful information on instrument characteristics and variations in the shallow crustal response to earthquakes, seasonal and ambient sources of seismic energy, along with atmospheric pressure–temperature changes and anthropogenic activities.

Geochemistry & Geophysics↗

Modeling of Seismic Waves Through Geologic Metamaterials

This project conducted a modeling study on seismic invisibility cloaks that render geologic targets invisible to seismic waves, using the concept of seismic metamaterials. We present a parametric numerical study on the behaviors of seismic waves through cloaks with different design parameters as well as degrees of geologic heterogeneity. In addition, a seismic cloaking strategy is proposed for a future field-scale experiment at a real-world test bed. This feasibility study will guide future field experiment designs and ultimately allow us to conduct systematic field-scale tests employing Sandia’s existing resources and field expertise. The ultimate goal is to develop methods and design parameters of seismic invisibility cloaks to protect against natural and man-made seismic waves. Seismic cloaking has potential applications in several areas of national security, energy, and natural hazard reduction.

58 GEOSCIENCES↗

​​Hydrogen Transportation and Distributed Energy Systems Seismic Risk Assessment for Cascadia Subduction Zone Airport Facilities​

Portland International Airport in Oregon is exploring operating a fleet of 28 fuel cell electric buses to support airport operations and possibly provide backup power during outages. In this report, we evaluate the risk of hydrogen deployment at an airport considering the potential for seismic activity. We present a seismic risk assessment for pressurized piping for a generic hydrogen refueling station at the airport. Only pressurized piping was considered in the risk analysis because of its vulnerability and availability of fragility information unlike the other components. Seismic capacity of other components should be included as more information becomes available. To characterize the seismic hazard, we considered the initiating event frequency of a Cascadia Subduction Zone earthquake event based on data from the literature. Seismic stressors that translate to the site from an earthquake event are expressed as peak ground acceleration. Site disturbance is characterized as a function of soil conditions, represented by different shear wave velocity values, which affect how seismic waves propagate and impact structures. Pipe fragility curves as a function of seismic capacity were used to correlate ground acceleration and failure probability. These fragility characteristics were combined with site-specific soil conditions to calculate the probability of pipe rupture during a seismic event. It should be noted that multiple simultaneous failures due to an earthquake event as a common cause are not considered. Some epistemic uncertainties such as self-exciting shocks, aftershocks, and the time variant nature of the ground motion are not considered either.

08 HYDROGEN↗

Prairie State Generating Company 2D Seismic Interpretation

The objectives of the Illinois Storage Corridor (ISC) project are to accelerate commercial deployment of carbon capture utilization and storage at two individual sites and receive approvals for Underground Injection Control (UIC) Class VI permits for construction at each site (ISC Project Narrative, 2020). As part of this project, and as part of the subsurface geologic characterization, two-dimensional (2D) seismic data was acquired at both sites. This report summarizes the findings from two phases of 2D acquisition and seismic interpretation at the Prairie State Generating Company site near Marissa, Illinois. The seismic data indicates the presence of three faults that completely transect the storage and confining units. The initial 2021 2D seismic acquisition revealed the presence of a feature two miles east of the Lively Grove #1 characterization well. The feature is a polyphase fault with a component of strike-slip motion, forming a small positive flower structure. The second phase of 2D seismic acquisition in 2022 constrained the maximum extent of the fault, indicating it has a relatively limited length. Other 2022 seismic lines revealed the presence of two other faults that transect the storage and confining units in the southwest and northeast portions of the 2D seismic acquisition area. The 2021 and 2022 2D seismic surveys have identified three specific locations within the project area that may have an elevated risk of out of zone CO 2 migration due to faulting. Further technical work will be needed to quantify this risk and determine how and if this risk will impact the placement of CO 2 injection wells across the project area.

01 COAL, LIGNITE, AND PEAT↗

Mantle Plume Magmatism in Elysium Planitia as Constrained by InSight Seismic Observations

Based on cratering statistics, there has been active volcanism in the last 100 Ma and possibly the last 1 to 10 Ma in both central Elysium Planitia and in the Cerberus Fossae region of south-eastern Elysium. Seismic measurements by the InSight mission show a strong concentration of seismicity in Cerberus Fossae which are interpreted as due to magma transport along dike systems. Together, these observations strongly support the existence of geologically recent and possibly of currently active magmatism in Elysium. InSight seismic observations were also used to construct seismic velocity models for the martian mantle. These models are effectively a regional aver-age of the seismic velocity structure for the Elysium region, where both the InSight lander and most of the measured seismic events are located. These seismic velocity models have been interpreted in terms of the lithospheric thickness and mantle potential temperature. Khan et al. found a lithospheric thickness of 400 to 600 km and a mantle potential temperature of 1325 to 1425 °C. Durán et al. found a lithospheric thickness of ~450 km and a mantle potential tempera-ture of 1375-1475 °C. Drilleau et al. found a thermal lithosphere of 420-660 km and a mantle potential temperature of 1380-1560 °C. In this study, we combine these InSight seismic observations with finite element mantle plume magma production models to place new constraints on the properties of the mantle beneath Elysium Planitia.

Walter S Kiefer↗

Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.

02 PETROLEUM↗