Search NASA⌕ Search

SEARCH · Search NASA

Results for “acoustic log”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Construction of the damped Ly⁢𝛼 absorber catalog for DESI DR2 Ly⁢𝛼 BAO

We present the Damped Ly⁢𝛼 Toolkit for automated detection and characterization of damped Ly⁢𝛼 absorbers (DLAs) in quasar spectra. Our method uses quasar spectral templates with and without absorption from intervening DLAs to reconstruct observed quasar forest regions. The best-fitting model determines whether a DLA is present while estimating the redshift and HI column density. With an optimized quality cut on detection significance (Δ⁢𝜒$^{2}_{𝑟}$ >0.03), the technique achieves an estimated 80% purity and 79% completeness when evaluated on simulated spectra with S/N>2 that are free of broad absorption lines (BALs). We provide a catalog containing candidate DLAs from the DLA Toolkit detected in DESI DR1 quasar spectra, of which 21 719 were found in S/N>2 spectra with predicted log 10 ⁡(𝑁 𝙷𝙸 )>20.3 and detection significance Δ⁢𝜒$^{2}_{𝑟}$ >0.03. We compare the Damped Ly⁢𝛼 Toolkit to two alternative DLA finders based on a convolutional neural network and Gaussian process models. We present a strategy for combining these three techniques to produce a high-fidelity DLA catalog from DESI DR2 for the Ly⁢𝛼 forest baryon acoustic oscillation measurement. The combined catalog contains 41 152 candidate DLAs with log 10 ⁡(𝑁 𝙷𝙸 )>20.3 from quasar spectra with S/N>2. We estimate this sample to be approximately 85% pure and 79% complete when BAL quasars are excluded.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Detecting fractures and monitoring hydraulic fracturing processes at the first EGS Collab testbed using borehole DAS ambient noise

Enhanced geothermal systems (EGS) require cost-effective monitoring of fracture networks. We validate the capability of using borehole distributed acoustic sensing (DAS) ambient noise for fracture monitoring using core photos and core logs. The EGS Collab project has conducted 10 m scale field experiments of hydraulic fracture stimulation using 50–60 m deep experimental wells at the Sanford Underground Research Facility (SURF) in Lead, South Dakota. The first EGS Collab testbed is located at 1616.67 m (4850 ft) depth at SURF and consists of one injection well, one production well, and six monitoring wells. All wells are drilled subhorizontally from an access tunnel called a drift. The project uses a single continuous fiber-optic cable installed sequentially in the six monitoring wells to record DAS data for monitoring hydraulic fracturing during stimulation. We analyze 60 s time records of the borehole DAS ambient noise data and compute the noise root-mean-square (rms) amplitude on each channel (points along the fiber cable) to obtain DAS ambient noise rms amplitude depth profiles along the monitoring wellbore. Our noise rms amplitude profiles indicate amplitude peaks at distinct depths. We compare the DAS noise rms amplitude profiles with borehole core photos and core logs and find that the DAS noise rms amplitude peaks correspond to the locations of fractures or lithologic changes indicated in the core photos or core logs. We then compute the hourly DAS noise rms amplitude profiles in two monitoring wells during three stimulation cycles in 72 h and find that the DAS noise rms amplitude profiles vary with time, indicating the fracture opening/growth or closing during the hydraulic stimulation. Our results demonstrate that borehole DAS passive ambient noise can be used to detect fractures and monitor fracturing processes in EGS reservoirs.

58 GEOSCIENCES↗
DOE OSTI2022

HFTS-2 (Final Report)

This report delves into the research program, identified tasks, various studies, observations, results, and conclusions under the Hydraulic Fracturing Test Site (HFTS-2) project. This test site is located in Loving County close to the Texas, New Mexico border. The program targeted the Wolfcamp formation in Delaware Basin. The aim of the project was to improve our understanding of hydraulic fracturing processes to not only optimize unconventional resource development in US Basins but to minimize environmental footprint as well. Various new diagnostic technologies were demonstrated at this site, including distributed acoustic microseismic, high resolution strain using Rayleigh frequency band (DSS-RFS), and in-fill proppant logging tool. Significant conclusions from this study include: 1) lateral and vertical fracture growth mapping as well as proppant transport behavior, 2) Use of DSS-RFS as a diagnostic tool to understand cluster level behavior, 3) Validation of proppant log, both at core scale as well as in-fill well implementation, and 4) New modeling tools to map far-field strain response using distributed strain sensing data. Specifically, far-field fracture distribution is influenced by linear fracture corridors defined by the stress-state as well as clustering behavior influenced by spatio-temporal effects. Depletion impacts on new stimulation was demonstrated and technically quantified. Significant vertical fracture height growth was observed. However, our understanding is limited by available gauge data availability. Numerical strain modeling techniques have been developed and demonstrated to help understand observed strain response behavior. Finally, proppant log has been further validated as a useful tool to map propped zones within the SRV. It has also shown very strong correlation with observed drawdown behavior at this test site. This report will detail many of these studies as well as technical results and conclusions from said studies. These include general subsurface characterization work, completion designs as implemented, studies looking at fracture geometry, high resolution microseismic study including source mechanisms, advanced DSS-RFS studies including drainage characterization, core characterization results, proppant analysis (core and in-fill child well), geochemistry, etc.

02 PETROLEUM↗
DOE OSTI2021

Persistent, “Mysterious” Seismoacoustic Signals Reported in Oklahoma State during 2019

Here, we report on the source of seismoacoustic pulses that were observed across the state of Oklahoma (OK) during summer of 2019, and the subject of national media coverage and speculation. Seismic network data collected across four U.S. states and interviews with witnesses to the pulse’s effect on residential structures demonstrate that they were triggered by routine ammunition disposal operations conducted by McAlester Army Ammunition Plant (McAAP). During these operations, conventional explosives destroy obsolete munitions stored in pits through a controlled sequence of electronically timed shots that occur over tens of minutes. Despite noise-abatement efforts that reduce coupling of acoustic energy with air, some lower frequency, subaudible (infrasonic) sound radiates from these shots as discrete pulses. We use nine months of blast log documents, seismic network records, analyst picks, and physical modeling to demonstrate that seismic stations as far as 640 km from McAAP sample these pulses, which record seasonal patterns in stratospheric and tropospheric winds, as well as the dynamic formation of waveguides and shadow zones. Digital short-term average to long-term average detectors that we augment with dynamic thresholds and time-binning operations identify these pulses with a fair probability, when compared with visual observations. Our analyses thereby provide estimates of observation rates for both partial and full sequences of these pulses, as well as single shots. We suggest that disposal operations can exploit existing, composite seismic networks to predict where residents are likely to witness blasting. Crucially, our data also show that dense seismic networks can record multiscale atmospheric processes in the absence of infrasound arrays.

58 GEOSCIENCES↗
DOE OSTI2020

Visualization of hydraulic fracture using physics-informed clustering to process ultrasonic shear waves

Ultrasonic transmission is sensitive to the spatial variation in mechanical properties of materials due to the presence of cracks/fractures. Wave propagation through fractured media introduces changes in the frequency content, travel time and transmission coefficient of the wave. A workflow based on physics-informed unsupervised learning is developed to process the transmitted ultrasonic-shear waveforms to non-invasively visualize the geomechanical alterations due to hydraulic fracturing. Novelty of the work involves the assignment of both statistically consistent and physically consistent clusters to the measurements of shear waveforms acquired across the one axial and two frontal planes. Physically consistent/relevant information is incorporated by considering the travel time of the peak of spectral energy and transmission coefficient of the transmitted waveform. The proposed workflow generates maps of geomechanical alterations across the frontal and axial planes of the sample. The outputs of the workflow are in good agreement with independent techniques viz. acoustic emission and X-ray computed tomography. Finally, the proposed workflow can be adapted for improved fracture characterization in the subsurface when processing sonic-logging, cross-wellbore seismic or surface seismic waveform data.

42 ENGINEERING↗
DOE OSTI2020

Oscillating scalar fields and the Hubble tension: A resolution with novel signatures

Here, we present a detailed investigation of a subdominant oscillating scalar field [“early dark energy” (EDE)] in the context of resolving the Hubble tension. Consistent with earlier work, but without relying on fluid approximations, we find that a scalar field frozen due to Hubble friction until log 10 ($z_c$) ~ 3.5, reaching $ρ_{\text{EDE}}(z_c)/ρ_{\text{tot}}$ ~ 10 % and diluting faster than matter afterwards, can bring cosmic microwave background (CMB), baryonic acoustic oscillations, supernovae luminosity distances, and the late-time estimate of the Hubble constant from the SH0ES Collaboration into agreement. A scalar field potential that scales as $V (\phi) \propto \phi ^{2n}$ with 2 ≲ $\textit{n}$ ≲ 3.4 around the minimum is preferred at the 68% confidence level, and the Planck polarization places additional constraints on the dynamics of perturbations in the scalar field. In particular, the data prefer a potential that flattens at large field displacements. A Markov-chain Monte Carlo analysis of mock data shows that the next-generation CMB observations (i.e., CMB-S4) can unambiguously detect the presence of the EDE at a very high significance. This projected sensitivity to the EDE dynamics is mainly driven by improved measurements of the E-mode polarization. We also explore new observational signatures of EDE scalar field dynamics: (i) We find that depending on the strength of the tensor-to-scalar ratio, the presence of the EDE might imply the existence of isocurvature perturbations in the CMB. (ii) We show that a strikingly rapid, scale-dependent growth of EDE field perturbations can result from parametric resonance driven by the anharmonic oscillating field for $\textit{n}$ ≈ 2. This instability and ensuing potentially nonlinear, spatially inhomogeneous, dynamics may provide unique signatures of this scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗
DOE OSTI2025

Validation of the DESI 2024 Lyα forest BAO analysis using synthetic datasets

The first year of data from the Dark Energy Spectroscopic Instrument (DESI) contains the largest set of Lyman-α (Lyα) forest spectra ever observed. This data, collected in the DESI Data Release 1 (DR1) sample, has been used to measure the Baryon Acoustic Oscillation (BAO) feature at redshift z = 2.33. In this work, we use a set of 150 synthetic realizations of DESI DR1 to validate the DESI 2024 Lyα forest BAO measurement presented in [1]. The synthetic data sets are based on Gaussian random fields using the log-normal approximation. We produce realistic synthetic DESI spectra that include all major contaminants affecting the Lyα forest. The synthetic data sets span a redshift range 1.8 < z < 3.8, and are analyzed using the same framework and pipeline used for the DESI 2024 Lyα forest BAO measurement. To measure BAO, we use both the Lyα auto-correlation and its cross-correlation with quasar positions. We use the mean of correlation functions from the set of DESI DR1 realizations to show that our model is able to recover unbiased measurements of the BAO position. We also fit each mock individually and study the population of BAO fits in order to validate BAO uncertainties and test our method for estimating the covariance matrix of the Lyα forest correlation functions. Finally, we discuss the implications of our results and identify the needs for the next generation of Lyα forest synthetic data sets, with the top priority being to simulate the effect of BAO broadening due to non-linear evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

CoLoRe-2LPT: Lyman-$α$ mock catalogues for the validation of DESI cosmological analyses

The Lyman-$α$ (Ly$α$) forest has become a crucial probe for studying the large-scale structure of the universe at high redshift ($z > 2$), providing powerful constraints on Baryon Acoustic Oscillations (BAO) and the full-shape (FS) clustering of matter. As a key ingredient for upcoming BAO and FS analyses, we present a new generation of fast cosmological Ly$α$ mocks based on second-order Lagrangian perturbation theory (2LPT). These new mocks significantly improve upon previous log-normal approaches, both at accurately capturing small scale clustering and at recovering the non-linear broadening of the BAO peak. They are able to reproduce Ly$α$ statistics within $10\%$ of the latest DESI measurement; including the Ly$α$ bias and the redshift-space distortion $β$ parameter, mean transmitted flux, and 1D power spectrum. The corresponding quasar (QSO) clustering is also improved with respect to previous approaches, calibrated against high-resolution Abacus simulations, recovering the observational QSO linear bias to less than $5\%$ and improving redshift-space distortions via 2LPT velocities and the addition of Fingers-of-God effects. Furthermore, these mocks incorporate high column density systems and metal lines, allowing us to explore the effects and systematics induced by these astrophysical contaminants. This new set of mocks has been key for enhancing the modeling and validation of the DESI DR2 Ly$α$ full shape cosmological analysis. This work provides a physically motivated and computationally efficient tool for simulating current and next-generation Ly$α$ surveys and validating FS and BAO analysis.

Bernal, M.F. Ruiz-Herrera [Madrid, CIEMAT] (ORCID:↗
DOE OSTI2020

Completion design improvement using a deep convolutional network

Maximizing stimulated natural and hydraulic fracture network is one of the primary hydraulic fracturing concerns for economic production from a horizontal shale gas well. Geomechanical facies and preexisting fractures in each stage are identified based on similarities in formation characteristics to optimize the locations of perforation clusters. This often requires analyzing large volumes of drilling, Logging While Drilling (LWD) and Measurement While Drilling (MWD) data. In this paper, we develop a methodology that calculates the mechanical specific energy (MSE) using real-time drill string acceleration signals directly from its definition. High resolution vibration signals have been collected using a tri-axial accerlometer, which was an auxiliary tool included in acoustic borehole imager. This technique provides a cost-efficient solution for engineered completion design. Furthermore, we adopt deep Convolutional Neural Network (CNN) with signal processing to build a data pipeline that effectively extracts patterns from dynamic acceleration signals for rock lateral MSE classification. First, we apply discrete wavelet transform and Short-Time Fourier Transform (STFT) for signal denoising and pattern recognition. Then we construct an image dataset using multi-scale image fusion at pixel level from 3 sensor channels, including axial, lateral acceleration spectrograms and zero-padded revolutions per minute (RPM). The resulted RGB image dataset includes 4,000 images of 5 MSE ranges with various rock strength conditions. Our results demonstrate that the proposed deep learning model can achieve more than 90% classification accuracy. The deep learning results, as a reference source, were applied in selected Marcellus Shale Energy and Environmental Lab (MSEEL) wells engineered completion located in the Marcellus shale gas site.

03 NATURAL GAS↗
DOE OSTI2025

Modeling of the high column density systems in the Lyman-Alpha forest

The Lyman-α forests observed in the spectra of high-redshift quasars can be used as a tracer of the cosmological matter density to study baryon acoustic oscillations (BAO) and the Alcock-Paczynski effect. Extraction of cosmological information from these studies requires modeling of the forest correlations. While the models depend most importantly on the bias parameters of the intergalactic medium (IGM), they also depend on the numbers and characteristics of high-column-density systems (HCDs) ranging from Lyman-limit systems with column densities log N HI /1cm -2 > 17 to damped Lyman-α systems (DLAs) with log N HI /1cm -2 > 20.2. These HCDs introduce broad damped absorption characteristic of a Voigt profile. Consequently they imprint a component on the power spectrum whose modes in the radial direction are suppressed, leading to a scale-dependent bias. Using mock data sets of known HCD content, we test a model that describes this effect in terms of the distribution of column densities of HCDs, the Fourier transforms of their Voigt profiles and the bias of the halos containing the HCDs. Our results show that this physically well-motivated model describes the effects of HCDs with an accuracy comparable to that of the ad-hoc models used in published forest analyses. We also discuss the problems of applying the model to real data, where the HCD content and their bias is uncertain.

Lyman alpha forest↗
DOE OSTI2023

Well proppants

Deployable smart acoustic resonance particles have dense cores, compliant matrixes surrounding the cores and stiff outer shells surrounding the matrixes. The particles have mechanical stress sensitivities that provide unique band gap shifts when compressed. Groups of similar particles with similar stress sensitivities and similar band gap shifts are added at different times to hydraulic fluids, as circulated through wells with the fluid and pushed into fractures. A plural, sonic monopole well logging tool is lowered into the well to determine locations and depth of fractures and local pressures by distinct resonance of individual groups.

Pollock, Jacob↗
DOE OSTI2020

Vitality models found useful in modeling tag-failure times in acoustic-tag survival studies

Abstract Acoustic telemetry studies often rely on the assumption that premature tag failure does not affect the validity of inferences. However, in some cases this assumption is possibly or likely invalid and it is necessary to apply a correction to estimation procedures. The question of which approaches and specific models are best suited to modeling acoustic tag failures has received little research attention. In this short communication, we present a meta-analysis of 42 acoustic tag-life studies, originally used to correct survival studies involving outmigrating juvenile salmonids in the Columbia/Snake river basin. We compare the performance of nine alternative parametric models including common failure–time/survival models and the vitality models of Li and Anderson Theor Popul Biol 76:118–131, (2009) and Demogr Res 28:341–372, (2013). The tag-life studies used acoustic tags from three different tag manufacturers, had expected lifetimes between 12 and 61 days, and had dry weights ranging from 0.22 to 1.65 g. In 57% of the cases, the vitality models of Li and Anderson Theor Popul Biol 76:118–131, (2009) and Demogr Res 28:341–372, (2013) fit the tag-failure times best. The vitality models were also the second-best choices in 17% of the cases. Together, the vitality models, log-logistic, (19%), and gamma models (14%) accounted for 90% of the models selected. Unlike more traditional failure–time models (e.g., Weibull, Gompertz, gamma, and log-logistic), the vitality models are capable of characterizing both the early onset of tag failure due to manufacturing errors and the anticipated battery life. We provide further guidance on appropriate sample sizes (50–100 tags) and procedures to be considered when applying precise tag-life corrections in release–recapture survival studies.

Skalski, John R. (ORCID:0000000270702505)↗
DOE OSTI2023

Slow dynamics in a single bead with mechanical conditioning and transient heating

The contact stiffness of an aluminum bead confined between two slabs diminishes upon mechanical conditioning, and then recovers like log(t) after the conditioning ceases. Here that structure is evaluated for its response to transient heating and cooling, with and without accompanying conditioning vibrations. It is found that, under heating or cooling alone, stiffness changes are mostly consistent with temperature dependent material moduli; there is little or no slow dynamics. Hybrid tests in which vibration conditioning is followed by heating or cooling lead to recoveries that begin like log(t) and then become more complex. On subtracting the known response to heating or cooling alone we discern the influence of higher or lower temperatures on slow dynamic recovery from vibrations. It is found that heating accelerates the initial log(t) recovery, but by an amount more than predicted by an Arrhenius model of thermally activated barrier penetrations. Furthermore, transient cooling has no discernable effect, in contrast to the Arrhenius prediction that it slows recovery.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Gas hydrate quantification in Walker Ridge block 313, Gulf of Mexico, from full-waveform inversion of ocean-bottom seismic data

The Gulf of Mexico (GOM) Joint Industry Project Leg 2 logging-while-drilling data in Walker Ridge lease block 313 (WR313) in the GOM detected gas hydrate in coarse- and fine-grained sediments at sites WR313-G and WR313-H. The coarse-grained units are thin ([Formula: see text]) and highly saturated, whereas the fine-grained unit is thick (approximately 200 m) with low saturation and fracture-filling gas hydrate. Unlike its coarse-grained counterpart, the seismic character of the fine-grained unit does not clearly indicate the presence of gas hydrate, which would likely have remained undiscovered in the absence of drilling. In this paper, through frequency-domain acoustic full-waveform inversion (FWI) of ocean-bottom seismometer data along a 2D multichannel seismic transect near sites WR313-G and WR313-H, we detect and quantify gas hydrate in the fine-grained unit. Key results are as follows: First, the base of the gas hydrate stability zone, which is not obvious in the reflection profile, can be discerned in the FWI results. Second, the gas hydrate in the fine-grained unit is mainly confined to the area between two sets of opposite-dipping normal faults implying that the fault architecture may be partially responsible for this gas hydrate accumulation and distribution.

Geochemistry & Geophysics↗
DOE OSTI2023

Time-lapse VSP integration and calibration of subsurface stress field utilizing machine learning approaches: A case study of the morrow B formation, FWU

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 a region around the injector well located within Farnsworth field unit (FWU), where there is an ongoing CO 2 -enhanced oil recovery (EOR) operation. In our study, a site-specific rock physics model was created from extensive geological, geophysical, and geomechanical characterization through 3D seismic data, well logs, and core assessed as part of the 1D MEM conducted on the characterization well within the study area. Here, the Biot-Gassmann workflow was utilized to combine the rock physics and reservoir simulation outputs to determine the seismic velocity change due to fluid substitution. Modeled seismic velocities attributed to mean effective stress were determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from ultrasonic seismic velocity measurements. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer (PSO) was utilized to minimize a penalty function created between the modeled seismic velocities and the observed time-lapse VSP dataset. 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. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.

58 GEOSCIENCES↗
DOE OSTI2020

Evaluation of a pair-based, joint-likelihood association approach for regional infrasound event identification

SUMMARY A Bayesian framework for the association of infrasonic detections is presented and evaluated for analysis at regional propagation scales. A pair-based, joint-likelihood association approach is developed that identifies events by computing the probability that individual detection pairs are attributable to a hypothetical common source and applying hierarchical clustering to identify events from the pair-based analysis. The framework is based on a Bayesian formulation introduced for infrasonic source localization and utilizes the propagation models developed for that application with modifications to improve the numerical efficiency of the analysis. Clustering analysis is completed using hierarchical analysis via weighted linkage for a non-Euclidean distance matrix defined by the negative log-joint-likelihood values. The method is evaluated using regional synthetic data with propagation distances of hundreds of kilometres in order to study the sensitivity of the method to uncertainties and errors in backazimuth and time of arrival. The method is found to be robust and stable for typical uncertainties, able to effectively distinguish noise detections within the data set from those in events, and can be made numerically efficient due to its ease of parallelization.

58 GEOSCIENCES↗
DOE OSTI2021

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↗