Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface stress”

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 37 records · Page 2

A Non-Invasive Approach for Elucidating the Spatial Distribution of In-Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage (Task 3 Report Laboratory Characterization of Stress Dependent Wavespeed)

This report presents a summary of Task 3 (Laboratory Experimental Characterization of Stress Dependent Wavespeed in Rocks and Application for In-Situ Stress Estimation) of the project titled “A Non-Invasive Approach for Elucidating the Spatial Distribution of In-Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage”. The main objective of this component of the project was to employ non-destructive techniques to determine how compressional (P)-wave and shear (S)-wave propagation velocities are related to triaxial-stress conditions in laboratory test samples of targeted reservoir rocks. The overarching principle is that such laboratory measurements can provide a basis whereby field measurements of wave-propagation velocities in rocks can be used to ascertain the stress state that would be expected to generate the observed wave-propagation behavior.

58 GEOSCIENCES↗

Determining Stress in Metallic Conducting Layers of Microelectronics Devices Using High Resolution Electron Backscatter Diffraction and Finite Element Analysis

Delayed failure due to stress voiding is a concern with some aging microelectronics, as these voids can grow large enough to cause an open circuit. Local measurements of stress in the metallic layers are crucial to understanding and predicting this failure, but such measurements are complicated by the fact that exposing the aluminum conducting lines will relieve most of their stress. In this study, we instead mechanically thin the device substrate and measure distortions on the thinned surface using high resolution electron backscatter diffraction (HREBSD). These measurements are then related to the stresses in the metallic layers through elastic simulations. This study found that in legacy components that had no obvious voids, the stresses were comparable to the theoretical stresses at the time of manufacture (≈300 MPa). Distortion fields in the substrate were also determined around known voids, which may be directly compared to stress voiding models. In conclusion, the technique presented here for stress determination, HREBSD coupled with finite element analysis to infer subsurface stresses, is a valuable tool for assessing failure in layered microelectronics devices.

HREBSD↗

Ultrasonic measurement of stress in 2219-T87 aluminum plate

The basic relationship of ultrasonic signal velocity to directional subsurface stress is reviewed. Inappropriateness of dependency on a single correlative value of constant for a three dimensional stress field in metallic materials is discussed. Implementation of conventional ultrasonic nondestructive testing capabilities integrated to provide a composite technique for the measurement of orthogonal stress components is described, and the procedures for performing the preparatory calibration and subsequent stress field measurements are presented. In conclusion, the prime effect of stress on ultrasonic signal velocity occurs only in the direction of material excitation or particle motion.

Clotfelter, W. N.↗

Stresses and deformations in elliptical contacts

Topics presented deal with defining conformal and nonconformal surfaces, curvature sum and difference, and surface and subsurface stresses in elliptical contacts. Load-deflection relationships for nonconformal contacts are developed. The deformation within the contact is, among other things, a function of the ellipticity parameter and elliptic integrals of the first and second kinds. Simplified expressions that allow quick calculations of the deformation to be made simply from a knowledge of the applied load, the material properties, and the geometry of the contacting elements are presented.

Hamrock, B. J.↗

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

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

58 GEOSCIENCES↗

Bearing tester data compilation, analysis, and reporting and bearing math modeling

The Shaberth bearing analysis computer program was developed for the analysis of jet engine shaft/bearing systems operating above room temperature with normal hydrocarbon lubricants. It is also possible to use this tool to evaluate the shaft bearing systems operating in cryogenics. Effects such as fluid drag, radial temperature gradients, outer race misalignments and clearance changes were simulated and evaluated. In addition, the speed and preload effects on bearing radial stiffness was evaluated. The Shaberth program was also used to provide contact stresses from which contact geometry was calculated to support other analyses such as the determination of cryogenic fluid film thickness in the contacts and evaluation of surface and subsurface stresses necessary for bearing failure evaluation. This program was a vital tool for the thermal analysis of the bearing in that it provides the heat generation rates at the rolling element/race contacts for input into a thermal model of the bearing/shaft assembly.

Source record↗

Simplified solution for stresses and deformation

Conventional contact deformation analysis for ball bearings, gears, and cams involves tedious iterative procedures or the use of design charts. A simplified approach that makes the elastic deformation at the center of contact easy to calculate was previously reported. The range of validity in which these equations can be used is extended. A simplified approach to the calculation of the location and magnitude of subsurface stresses developed in machine element applications is included.

Hamrock, B. J.↗

Enhanced Microseismicity During Production Pumping Cessation at the San Emidio Geothermal Field (Nevada, USA) in December 2016

Abstract Tectonic activity, geothermal fluids, and microseismic events (MSEs) tend to occur in similar locations as a result of spatiotemporal changes in the subsurface stress state. To quantify this association, we analyze data from a dense seismic array deployed at the San Emidio geothermal field, Nevada for 1 week in December 2016 to coincide with a 19.45‐hr shutdown of all injection and production pumping operations. 123 MSEs were detected, of which 101 occurred during the shutdown. The spatial association of the MSEs with the production wells suggests a causal relationship between the production cessation and the MSEs. Here we performed a detailed analysis to investigate reservoir material properties, distribution of seismically activated faults, and local stress state. We determined the hypocenters, magnitudes, and focal mechanisms for the MSEs, P‐wave tomographic velocity model, and local stress tensor. The results show that most MSEs occurred near the production wells. Magnitudes fall between −2.2 and 0.0 with larger events located closer to the production wells. Most MSEs occurred within a westward‐dipping normal fault zone in the reservoir associated with anomalously low P‐wave velocity values. The focal mechanism and stress inversion results show predominantly normal faulting with the maximum horizontal stress oriented north‐south. We suggest that the MSEs during shutdown were triggered on pre‐existing, small‐scale, critically stressed fault patches in the reservoir as the pore pressure increased around the production wells when the production pumping ceased. We interpret the larger MSE magnitudes closer to the production wells as a result of higher pore pressure increase.

15 GEOTHERMAL ENERGY↗

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗

Unique Aeolian Transport Mechanisms on Mars: Respective Roles of Percussive and Repercussive Grain Populations in the Sediment Load

Experiments show that when sand-size grains impact a sediment surface with energy levels commensurate for Mars, small craters are formed by the ejection of several hundred grains from the bed. The experiments were conducted with a modified crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism. Individual grains of sand could be fired at loose sand targets to observe ballistic effects unhindered by aerodynamic mobilization of the bed. Impact trajectories simulated the saltation process on dune surfaces. Impact craters were not elongated despite glancing (15 deg.) bed impact; the craters were very close to being circular. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 um-diameter grains into similar material. This behavior is explained by deposition of elastic energy in the bed by the "percussive" grain. Impact creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains to open packing and they consequently become forcefully ejected from the site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. A stress model based on repercussive bed dilatancy and interparticle adhesive forces (for smaller grains) predicts, to first order, the observed crater volumes for various impact conditions. On earth, only a few grains are mobilized by a percussive saltating grain; some grains are "knudged" along the ground, and some are partly expelled on short trajectories. These motions constitute reptation transport. On Mars, saltation and reptation become indistinct: secondary or "repercussive" trajectories have sufficient vertical impulse to create a dense saltation population of many tens or hundreds of grains for each single high-speed saltation percussion of the bed. Impact cascading will lead to near-surface distortion of the boundary layer, and choked flow formed by a dense "slurry" of sand, with the majority of grains mobilized by repercussive forces rather than by aerodynamic lift. This proceeds until a fully-matured transport layer imposes self-limitations as grain-population density constrains the free-path motion of individial grains.

Marshall, John R.↗

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗

Unique Aeolian Transport Mechanisms on Mars: Respective Roles of Percussive and Repercussive Grain Populations in the Sediment Load

Experiments show that when sand-size grains impact a sediment surface with energy levels commensurate for Mars, small craters are formed by the ejection of several hundred grains from the bed. The experiments were conducted with a modified crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism. Individual grains of sand could be fired at loose sand targets to observe ballistic effects unhindered by aerodynamic mobilization of the bed. Impact trajectories simulated the saltation process on dune surfaces. Impact craters were not elongated despite glancing (15 deg.) bed impact; the craters were very close to being circular. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. This behavior is explained by deposition of elastic energy in the bed by the "percussive" grain. Impact creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains to open packing and they consequently become forcefully ejected from the site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. A stress model based on repercussive bed dilatancy and interparticle adhesive forces (for smaller grains) predicts, to first order, the observed crater volumes for various impact conditions. On earth, only a few grains are mobilized by a percussive saltating grain; some grains are "knudged" along the ground, and some are partly expelled on short trajectories. These motions constitute reptation transport. On Mars, saltation and reptation become indistinct: secondary or "repercussive" trajectories have sufficient vertical impulse to create a dense saltation population of many tens or hundreds of grains for each single high-speed saltation percussion of the bed. Impact cascading will lead to near-surface distortion of the boundary layer, and choked flow formed by a dense "slurry" of sand, with the majority of grains mobilized by repercussive forces rather than by aerodynamic lift. This proceeds until a fully-matured transport layer imposes self- limitations as grain-population density constrains the free-path motion of individual grains.

Marshall, John R.↗

Nondestructive Superresolution Imaging of Defects and Nonuniformities in Metals, Semiconductors, Dielectrics, Composites, and Plants Using Evanescent Microwaves

We have imaged and mapped material nonuniformities and defects using microwaves generated at the end of a microstripline resonator with 0.4 micrometer lateral spatial resolution at 1 GHz. Here we experimentally examine the effect of microstripline substrate permittivity, the feedline-to-resonator coupling strength, and probe tip geometry on the spatial resolution of the probe. Carbon composites, dielectrics, semiconductors, metals, and botanical samples were scanned for defects, residual stresses, subsurface features, areas of different film thickness, and moisture content. The resulting evanescent microwave probe (EMP) images are discussed. The main objective of this work is to demonstrate the overall capabilities of the EMP imaging technique as well as to discuss various probe parameters that can be used to design EMPs for different applications.

Tabib-Azar, M.↗

The 2002 NASA Faculty Fellowship Program Research Reports

Contents include the following: System Identification of X-33. Neural Network Advanced Ceramic Technology for Space Applications at NASA MSFC. Developing a MATLAB-Based Tool for Visualization and Transformation. Subsurface Stress Fields in Single Crystal (Anisotropic). Contacts Our Space Future: A Challenge to the Conceptual Artist Concept Art for Presentation and Education. Identification and Characterization of Extremophile Microorganisms. Significant to Astrobiology. Mathematical Investigation of Gamma Ray and Neutron. Absorption Grid Patterns for Homeland Defense-Related Fourier Imaging Systems. The Potential of Microwave Radiation for Processing Martian Soil. Fuzzy Logic Trajectory Design and Guidance for Terminal Area.

Bland, J.↗

Software Developed for Analyzing High- Speed Rolling-Element Bearings

COBRA-AHS (Computer Optimized Ball & Roller Bearing Analysis--Advanced High Speed, J.V. Poplawski & Associates, Bethlehem, PA) is used for the design and analysis of rolling element bearings operating at high speeds under complex mechanical and thermal loading. The code estimates bearing fatigue life by calculating three-dimensional subsurface stress fields developed within the bearing raceways. It provides a state-of-the-art interactive design environment for bearing engineers within a single easy-to-use design-analysis package. The code analyzes flexible or rigid shaft systems containing up to five bearings acted upon by radial, thrust, and moment loads in 5 degrees of freedom. Bearing types include high-speed ball, cylindrical roller, and tapered roller bearings. COBRA-AHS is the first major upgrade in 30 years of such commercially available bearing software. The upgrade was developed under a Small Business Innovation Research contract from the NASA Glenn Research Center, and incorporates the results of 30 years of NASA and industry bearing research and technology.

Fleming, David P.↗

Small Volcanic Vents of the Tharsis Volcanic Province, Mars

Distributed-style volcanism is an end member of terrestrial volcanism that produces clusters of small volcanoes when isolated magma bodies ascend from broad magma source regions. Volcano clusters can develop over millions of years, one volcano at a time, and can be used to infer unobserved geologic phenomena, including subsurface stresses and cracks during eruption periods. The Tharsis Volcanic Province covers approximately one quarter of the Martian surface and hosts a large concentration of small volcanoes that formed from distributed volcanism. We present a catalog of 1106 small volcanic vents identified within Tharsis Volcanic Province. This catalog includes morphologic measurements for each cataloged vent. Vent lengths range from 71 m to 51 km, widths range from 40 m to 3.1 km, and 90% of vents have lengths at least 1.5 times their widths. Additionally, 90% of edifices associated with vents have topographic prominences <100 m. Vents are found throughout Tharsis, though they generally form clusters near large volcanoes or among large graben sets. Older regions with volcanic eruption ages of >1 Ga are found at the Tharsis periphery in the Tempe-Mareotis region and Syria Planum. Vents in the Tharsis interior have reported ages <500 Ma. Regional trends in vent orientation and intervent alignment are dependent on nearby central volcanoes and fossae. We use these findings to hypothesize that within the most recent 500 Ma, magma was present under and to the east of the Tharsis Montes and that some of this magma erupted and built hundreds of small volcanoes in this region.

Distributed volcanism↗

Geothermal Hydraulic Stimulation: Overview of Methods and Best Practices

This paper outlines important information related to the hydraulic stimulation of geothermal wells and considerations for regulators responsible for evaluating applications to conduct these operations. Hydraulic stimulation applied to geothermal wells is an evolving technology that is different from the hydraulic fracturing process widely used to complete unconventional oil and gas wells. Important differences relate to operating pressures, variations in rock and fracture processes, and characteristics and disposal of hydraulic fracture fluids. Hydraulic stimulation is not widely used in the geothermal industry; however, future development of enhanced geothermal systems (EGS) will require hydraulic stimulation in order to enhance and create fracture permeability to allow wells to deliver sufficient heat and fluid to power plants. Gigawatts of EGS potential have been identified (U.S. Department of Energy 2019), so it is expected that the use of hydraulic stimulation of geothermal wells will be more common in the future. At present, there are no formal regulations guiding drilling programs or sundry notices that propose hydraulic stimulation of geothermal wells. On federal lands, important constraints and oversight are embedded in the process for obtaining Geothermal Drilling Permits outlined in the Code of Federal Regulations (43 CFR Section 3262.11); in Geothermal Resource Order 2 (GRO 2) guidance for assessment and mitigation of impacts of geothermal operations as well as construction and testing of geothermal wells; and in a Bureau of Land Management (BLM) Induced Seismicity Instruction Memorandum (BLM 2018). Additionally, local knowledge from BLM field offices and expertise of BLM engineers, or that of state regulators for operations on private and state lands, will be applied to fully assess operator applications to conduct geothermal hydraulic stimulation activities. Significant research effort has been directed toward EGS, so geothermal hydraulic stimulation techniques, impacts, and results are evolving and becoming better understood, including adaptation of decades of oil and gas industry experience with hydraulic fracturing of unconventional reservoirs. Recent EGS activities in the United States provide details about hydraulic stimulation with respect to oversight, stimulation design, execution, and results. Based on experience to date, the important issues to address when hydraulically stimulating a geothermal well include the following: 1. Wellbore construction and integrity must be appropriate to protect groundwater and manage stimulation pressures. 2. Understanding of lithology, faults, fractures, and subsurface stress state is necessary to design stimulation plans and predict results. 3. Seismic monitoring allows for observation and mitigation of induced seismicity.

geothermal↗

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↗