Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D seismic data”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Enhancing trap and fault seal analyses by integrating observations from HR3D seismic data with well logs and conventional 3D seismic data, Texas inner shelf

An understanding of trap and fault seal quality is critical for assessing hydrocarbon prospectivity. To achieve this, modern analytical techniques leverage well data and conventional industry-standard 3D seismic data to evaluate the trap, and any faults displacing the reservoir and top seal intervals. Above all, geological interpretation provides the framework of trap and fault seal analyses, but can be hindered by the data resolution, quality and acquisition style of the conventional seismic data. Furthermore, limiting the analysis to only the petroleum system at depth may lead to erroneous perceptions because interpreting overburden features, such as shallow faults or gas chimneys, can provide valuable observations with respect to container performance, and can to help validate trap and fault seal predictions. Here, a supplement to conventional 3D data are high-resolution 3D seismic (HR3D) data, which provide detailed images of the overburden geology. This study utilizes an HR3D seismic volume in the San Luis Pass area of the Texas inner shelf, where shallow fault tips and a sizeable gas chimney are interpreted over an unsuccessful hydrocarbon prospect. Static post-drill fault seal and trap analyses suggest that the primary fault displacing the structural closure could have withheld columns of gas c. 100 m high, but disagree with our HR3D seismic interpretations and dry-well analyses. From our results, we hypothesize that tertiary gas migration through fault conduits reduced the hydrocarbon column in the prospective Early Miocene reservoir, and may have resulted from continued movement along the intersecting faults. Overall, this study reinforces the importance of understanding the overburden geology and geohistory of faulted prospects, and demonstrates the utility of pre-drill HR3D acquisition when conducting trap and fault seal analyses.

3D seismic↗
DOE OSTI2018

Utah FORGE: 2D and 3D Seismic Data

This set of data contains raw and Initially-processed 2D and 3D seismic data from the Utah FORGE study area near Roosevelt Hot Springs. Reprocessed versions of these data can be accessed at the linked submission in the Resources section titled "Reprocessed Seismic Reflection Data." The zipped archives numbered from 1-100 to 1001-1122 contain 3D seismic uncorrelated shot gatherers SEG-Y files. The zipped archives numbered from 1-100C to 1001-1122C contain 3D seismic correlated shot gatherers SEG-Y files. Other data have intuitive names.

15 GEOTHERMAL ENERGY↗
DOE OSTI2017

Diffraction imaging for seal evaluation using ultra high resolution 3D seismic data

Cap rock is any impermeable or low permeability formation that may trap oil, gas or water, preventing it from migrating to the surface. Evaluation of sealing properties of the cap rock is a critical task for seismic exploration. It is important for oil and gas prospecting as well as for locating reservoirs appropriate for carbon dioxide storage. The latter has been proposed as one solution to global climate change caused by heat-trapping anthropogenic gasses in the atmosphere. A fluid escaped from the reservoir fills voids in the overlying strata thereby causing local changes of acoustic properties. We aim to explore such fluid-saturated areas and thus evaluate sealing properties of underlying cap rocks. We assume that local acoustic impedance changes caused by fluid migration generate seismic diffractions. We isolate diffracted waves from P-Cable seismic data acquired over the area of interest and perform diffraction analysis of the near-surface interval (above 200 m). In the diffraction image obtained, we interpret an extensive system of broad sinuous channels, narrow linear faults, and diffuse fracture networks. In addition, we detect high diffractivity anomalies, which extend vertically and laterally. We note no evidence of these anomalies in a conventional seismic image and associated seismic attribute (similarity and fracture density) volumes. We associate this high diffractivity to a variable pore fluid composition and saturation, which we infer to be due to hydrocarbon migration from the underlying strata. Consequently, we conclude poor seal quality for the underlying cap rocks. During data processing, we test diffraction focusing analysis for velocity model building. In conclusion, this option appears quite valuable for P-Cable data processing because of short offset and limitations of the standard velocity analysis.

58 GEOSCIENCES↗
DOE OSTI2017

Numerical Simulations of 3D Seismic Data Final Report CRADA No. TC02095.0

This was a collaborative effort between Lawrence Livermore National Security, LLC (formerly The Regents of the University of Califomia)/Lawrence-Livermore National Laboratory (LLNL) and Schlumberger Cambridge Research (SCR), to develop synthetic seismic data sets and supporting codes.

02 PETROLEUM↗
DOE OSTI2017

Gas hydrate characterization from a 3D seismic dataset in the deepwater eastern Gulf of Mexico

Principal component analysis of spectral decomposition results combined with amplitude and frequency seismic attributes derived from 3D seismic data are used for the identification and characterization of gas hydrate deposits in the deepwater eastern Gulf of Mexico. In the central deepwater Gulf of Mexico (GoM), logging while drilling LWD data provided insight to the amplitude response of gas hydrate saturation in sands, which could be used to characterize complex gas hydrate deposits in other sandy deposits. In this study, a large 3D seismic data set from equivalent and distal Plio Pleistocene sandy channel deposits in the deepwater eastern Gulf of Mexico is screened for direct hydrocarbon indicators for gas hydrate saturated sands.

gas hydrate↗
DOE OSTI2017

Gas Hydrate Characterization from a 3D Seismic Dataset in the Eastern Deepwater Gulf of Mexico

The presence of a gas hydrate petroleum system and seismic attributes derived from 3D seismic data are used for the identification and characterization of gas hydrate deposits in the deepwater eastern Gulf of Mexico. In the central deepwater Gulf of Mexico (GoM), logging while drilling (LWD) data provided insight to the amplitude response of gas hydrate saturation in sands, which could be used to characterize complex gas hydrate deposits in other sandy deposits. In this study, a large 3D seismic data set from equivalent and distal Plio-Pleistocene sandy channel deposits in the deepwater eastern Gulf of Mexico is screened for direct hydrocarbon indicators for gas hydrate saturated sands.

gas hydrate↗
DOE OSTI2015

Prestack reverse time migration for 3D marine reflection seismic data

Prestack reverse time migration (RTM) is a method for imaging the subsurface using the inner product of wavefield extrapolation in shot domain and in receiver domain. It is well known that RTM is better for preserving amplitudes and phases than other prestack migrations. Since 3D seismic data is huge data volume and it needs heavy computing works, it requires parallel computing in order to have a meaningful depth image of the 3D subsurface. We implemented a parallelized version of 3D RTM for prestack depth migration. The results of numerical example for 3D SEG/EAGE salt model showed good agreement with the original geological model. We applied RTM to offshore 3D seismic reflection data. The study area is 12 × 25 km with 120 survey lines. Shot and receiver spacing is 25 m and 12.5 m. The line spacing is 100 m. Shot gathers were preprocessed to enhance signal to noise ratio and velocity model was calculated from conventional stack velocity. Both of them were used to obtain 3D image using RTM. The results show reasonable subsurface image.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

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↗
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↗
DOE OSTI2014

Analysis of ancient-river systems by 3D seismic time-slice technique: A case study in northeast Malay Basin, offshore Terengganu, Malaysia

Fluvial sandstones constitute one of the major clastic petroleum reservoir types in many sedimentary basins around the world. This study is based on the analysis of high-resolution, shallow (seabed to 500 m depth) 3D seismic data which generated three-dimensional (3D) time slices that provide exceptional imaging of the geometry, dimension and temporal and spatial distribution of fluvial channels. The study area is in the northeast of Malay Basin about 280 km to the east of Terengganu offshore. The Malay Basin comprises a thick (> 8 km), rift to post-rift Oligo-Miocene to Pliocene basin-fill. The youngest (Miocene to Pliocene), post-rift succession is dominated by a thick (1–5 km), cyclic succession of coastal plain and coastal deposits, which accumulated in a humid-tropical climatic setting. This study focuses on the Pleistocene to Recent (500 m thick) succession, which comprises a range of seismic facies analysis of the two-dimensional (2D) seismic sections, mainly reflecting changes in fluvial channel style and river architecture. The succession has been divided into four seismic units (Unit S1-S4), bounded by basin-wide strata surfaces. Two types of boundaries have been identified: 1) a boundary that is defined by a regionally-extensive erosion surface at the base of a prominent incised valley (S3 and S4); 2) a sequence boundary that is defined by more weakly-incised, straight and low-sinuosity channels which is interpreted as low-stand alluvial bypass channel systems (S1 and S2). Each unit displays a predictable vertical change of the channel pattern and scale, with wide low-sinuosity channels at the base passing gradationally upwards into narrow high-sinuosity channels at the top. The wide variation in channel style and size is interpreted to be controlled mainly by the sea-level fluctuations on the widely flat Sunda land Platform.

58 GEOSCIENCES↗
DOE OSTI2016

The character and amplitude of ‘discontinuous’ bottom-simulating reflections in marine seismic data

Bottom-simulating reflections (BSRs) identified in seismic data are well documented; and are commonly interpreted to indicate the presence of gas hydrates along continental margins, as well as to estimate regional volumes of gas hydrate. A BSR is defined as a reflection that sub-parallels the seafloor but is opposite in polarity and cross-cuts dipping sedimentary strata. BSRs form as a result of a strong negative acoustic impedance contrast. BSRs, however, are a diverse seismic phenomena that manifest in strikingly contrasting ways in different geological settings, and in different seismic data types. We investigate the characteristics of BSRs, using conventional and high resolution, 2D and 3D seismic data sets in three locations: the Terrebonne and Orca Basins in the Gulf of Mexico, and Blake Ridge on the US Atlantic Margin. The acquisition geometry and frequency content of the seismic data significantly impact the resultant character of BSRs, as observed with depth and amplitude maps of the BSRs. Furthermore, our amplitude maps reinforce the concept that the BSR represents a zone, over which the transition from hydrate to free gas occurs, as opposed to the conventional model of the BSR occurring at a single interface. Our results show that a BSR can be mapped in three dimensions but it is not spatially continuous, at least not at the basin scale. Rather, a BSR manifests itself as a discontinuous, or patchy, reflection and only at local scales is it continuous. We suggest the discontinuous nature of BSRs is the result of variable saturation and distribution of free gas and hydrate, acquisition geometry and frequency content of the recorded seismic data. Lastly, the commonly accepted definition of a BSR should be broadened with careful consideration of these factors, to represent the uppermost extent of enhanced amplitude at the shallowest occurrence of free gas trapped by overlying hydrate-bearing sediments.

03 NATURAL GAS↗
DOE OSTI2014

Potosi Reservoir Modeling; History and Recommendations

As a part of a larger project co-funded by the United States Department of Energy (US DOE) to evaluate the potential of formations within the Cambro-Ordovician strata above the Mt. Simon as potential targets for carbon sequestration in the Illinois and Michigan Basins, the Illinois Clean Coal Institute (ICCI) requested Schlumberger to evaluate the potential injectivity and carbon dioxide (CO₂) plume size of the Cambrian Potosi Formation. The evaluation of this formation was accomplished using wireline data, core data, pressure data, and seismic data from two projects: the US DOE-funded Illinois Basin–Decatur Project being conducted by the Midwest Geological Sequestration Consortium in Macon County, Illinois, as well as data from the Illinois – Industrial Carbon Capture and Sequestration (IL-ICCS) project funded through the American Recovery and Reinvestment Act. In 2010, technical performance evaluations on the Cambrian Potosi Formation were performed through reservoir modeling. The data included formation tops from mud logs, well logs from the Verification Well 1 (VW1) and the Injection Well (CCS1), structural and stratigraphic formation from three dimensional (3D) seismic data, and field data from several waste water injection wells for the Potosi Formation. The intention was for two million tonnes per annum (MTPA) of CO₂ to be injected for 20 years into the Potosi Formation. In 2013, updated reservoir models for the Cambrian Potosi Formation were evaluated. The data included formation tops from mud logs, well logs from the CCS1, VW1, and Verification Well 2 (VW2) wells, structural and stratigraphic formation from a larger 3D seismic survey, and field data from several waste water injection wells for Potosi Formation. The objective is to simulate the injection of CO₂ at a rate 3.5 million tons per annum (3.2 million tonnes per annum [MTPA]) for 30 years 106 million tons (96 MT total) into the Potosi Formation. The Potosi geomodeling efforts have evolved from using data from a single well in 2010 to the inclusion of data from three wells in 2013 which largely leverage the porosity and permeability logs plus knowledge of lost circulation zones. The first Potosi model (Potosi Geobody Model 2010) attempted to use the available seismic inversion data to inform the geomodel and predict vugular zones in advance of drilling VW1. Lost circulation zones in VW1 came in as the geologists predicted. The model was not implemented in subsequent simulation work. To date, reservoir models used for flow simulation work have relied predominately on Gaussian distributed properties (porosity and permeability) and have employed a single injection well. Potosi Model 2013b incorporated the new VW2 logs, and exhibited an extra level of sophistication by delineating the vugular intervals. This method added further realism that likely represents the best reservoir approximation to date. Where the 2010 reservoir models were 10 by 10 mi (16 by 16 km) in area, the 2013 models were expanded in size to 30 by 30 mi (48 by 48 km). The latest reservoir simulations show that a minimum of four injectors might be required to meet target injection rates. Still, there is data that requires further scrutiny and modeling methodologies that require testing for the Potosi Formation. This work is currently ongoing, and the next phase of the reservoir modeling intends to implement valuable data like porosity derived from seismic inversion, seismically derived geobodies, or a combination of both to further define vugular zones and the porosity distribution within the Potosi Formation. Understanding the dual porosity, dual permeability character of the Potosi remains the greatest challenge in representing this formation. Further analysis of the FMI* fullbore formation microimager data may aid in assessing this uncertainty. The Potosi Formation is indeed an interesting formation, and recommendations to further characterize it are included in the following list: - Data acquisition to identify the vugs permeability, distribution, and interconnectivity could be considered to perform a more rigorous evaluation of the Potosi Formation injectivity and capacity. This could be achieved by performing an injection test on a vugular interval to determine the vugs permeability, and an interference test between wells to evaluate the local vugs extent and interconnectivity. - A thorough study of the available FMI data may reveal specifics on estimating the vug to matrix ratio. This estimate could be used to further condition the porosity distribution. Porosity logs alone might underestimate the formation’s porosity associated with vugs. Porosity mapping derived from the seismic inversion could also be used in the succeeding task to characterize the lateral porosity distribution within the reservoir. This could involve the geobody methodology previously attempted in 2010. With or without seismic inversion porosity mapping, it is worth exploring whether increased lateral heterogeneity plays a significant role in Potosi injectivity. Investigations on vugular, dolomitic outcrops suggest that there may be significantly greater lateral heterogeneity than what has been modeled here. - The FMI data also reveals the presence of and helps describe open fractures. The presence of fractures will further enhance the formation’s permeability. The task of leveraging this data in the geomodeling effort still remains. Under the best of circumstances, this data describing open fractures may be combined with seismic attributes to delineate fracture corridors. Fracture modeling would certainly add another layer of sophistication to the model. Its contribution and applicability remain to be explored. - Facies modeling within the Potosi has yet to be thoroughly addressed. The carbonates during the time of deposition are believed to be regionally extensive. However, it may be worth delineating the reservoir with other regional wells or modern day analogues to understand the extent of the Potosi. More specifically, the model could incorporate lateral changes or trends if deemed necessary to represent facies transition. - Currently there is no fracture gradient data available for the Potosi in the Decatur project area. The acquisition of the fracture pressure data could be considered to determine an appropriate maximum allowable bottomhole injection pressure. This would allow the evaluation of injectivity and the required number of wells in a more precise manner. - Special core analysis (SCAL) to determine the relative permeability and capillary pressure of the vugs and matrix could be considered to have a better estimation of the reservoir injectivity and plume extent. - Formation water sampling and analysis could be considered for the Potosi to estimate the water salinity and properties. A vertical flow performance evaluation could be considered for the succeeding task to determine the appropriate tubing size, the required injection tubing head pressure, and to investigate whether the corresponding well injection rate falls within the tubing erosional velocity limit. - A simulation using several injectors could also be considered to determine the required number of wells to achieve the injection target while taking into account the pressure interference.

Smith, Valerie↗
DOE OSTI2015

Factors Affecting Seismic Velocity in Alluvium

Yucca Flat at the Nevada National Security Site has been selected as the Source Physics Experiment (SPE) Dry Alluvium Geology Phase II site. The alluvium in this part of Yucca Flat is typical of desert basin fill, with discontinuous beds that are highly variable in clast size and provenance. Detailed understanding of the subsurface geology will be needed for interpretation of the SPE seismic data. A 3D seismic velocity model, created for Yucca Flat using interval seismic velocity data, shows variations in velocity within alluvium near the SPE Phase II site beyond the usual gradual increase of density with depth due to compaction. In this study we examined borehole lithologic logs, geophysical logs, downhole videos, and laboratory analyses of sidewall core samples to understand which characteristics of the alluvium are related to these variations in seismic velocity. Seismic velocity of alluvium is generally related to its density, which can be affected by sediment provenance, clast size, gravel percentage, and matrix properties, in addition to compaction. This study presents a preliminary subdivision of the alluvial strata in the SPE Phase II area into mappable units expected to be significant to seismic modeling. Further refinements of the alluvial units may be possible when seismic data are obtained from SPE Phase II tests.

58 GEOSCIENCES↗
DOE OSTI2019

Utah FORGE: Seismic Reflection Data

This is 2D and 3D seismic reflection data from Utah FORGE reprocessed during Phase 2c. The readme file containing an explanation of the data including data formats, software that can be used, processing, and projection and datum used. The Reprocessing document gives the rationale for reprocessing and shows examples of the improvements that were obtained. For all 3D and 2D data the following datasets were created and output in SEG-Y format: - Unmigrated Time - Prestack Time Migration (PSTM), Unenhanced (UnEnh) and Enhanced (Enh) - Prestack Depth Migration (PSDM), Unenhanced (UnEnh) and Enhanced (Enh) - Velocity Model used for Migration

15 GEOTHERMAL ENERGY↗
DOE OSTI2014

The Potosi Reservoir Model 2013

As a part of a larger project co-funded by the United States Department of Energy (US DOE) to evaluate the potential of formations within the Cambro-Ordovician strata above the Mt. Simon as potential targets for carbon sequestration in the Illinois and Michigan Basins, the Illinois Clean Coal Institute (ICCI) requested Schlumberger to evaluate the potential injectivity and carbon dioxide (CO2) plume size of the Cambrian Potosi Formation. The evaluation of this formation was accomplished using wireline data, core data, pressure data, and seismic data from the US DOE-funded Illinois Basin–Decatur Project (IBDP) being conducted by the Midwest Geological Sequestration Consortium in Macon County, Illinois. In 2010, technical performance evaluations on the Cambrian Potosi Formation were performed through reservoir modeling. The data included formation tops from mud logs, well logs from the VW1 and the CCS1 wells, structural and stratigraphic formation from three dimensional (3D) seismic data, and field data from several waste water injection wells for Potosi Formation. Intention was for two million tons per annum (MTPA) of CO2 to be injected for 20 years. In the preceding, the 2010 Potosi heterogeneous model (referred to as the "Potosi Dynamic Model 2010" in this topical report) was re-run using a new injection scenario; 3.2 MTPA for 30 years. The extent of the Potosi Dynamic Model 2010, however, appeared too small for the new injection target. It was not sufficiently large enough to accommodate the evolution of the plume. The new model, Potosi Dynamic Model 2013a, was built by extending the Potosi Dynamic Model 2010 grid to 30 miles x 30 miles (48.3km x48.3km), while preserving all property modeling workflows and layering. This model was retained as the base case of Potosi Dynamic Model 2013a. The Potosi reservoir model was updated to take into account the new data from the verification well VW2 which was drilled in 2012. The new porosity and permeability modeling was performed to take into account the log data from the new well. Revisions of the 2010 modeling assumptions were also done on relative permeability, capillary pressures, formation water salinity, and the maximum allowable well bottomhole pressure. Dynamic simulations were run using the injection target of 3.2 MTPA for 30 years. This new dynamic model was named Potosi Dynamic Model 2013b. Due to the major uncertainties on the vugs permeability, two models were built; the Pessimistic and Optimistic Cases. The Optimistic Case assumes vugs permeability of 9,000 mD, which is analog to the vugs permeability identified in the pressure fall off test of a waste water injector in the Tuscola site, approx. 40 miles (64.4km) away from the IBDP area. The Pessimistic Case assumes that the vugs permeability is equal to the log data, which does not take into account the permeability from secondary porosity. The probability of such case is deemed low and could be treated as the worst case scenario, since the contribution of secondary porosity to the permeability is neglected and the loss circulation events might correspond to a much higher permeability. It is considered important, however, to identify the range of possible reservoir performance since there are no rigorous data available for the vugs permeability. The Optimistic Case gives an average CO2 injection rate of 0.8 MTPA and cumulative injection of 26 MT in 30 years, which corresponds to 27% of the injection target. The injection rate is approx. 3.2 MTPA in the first year as the well is injecting into the surrounding vugs, and declines rapidly to 0.8 MTPA in year 4 once the surrounding vugs are full and the CO2 start to reach the matrix. This implies that according to this preliminary model, a minimum of four (4) wells could be required to achieve the injection target. This result is lower than the injectivity estimated in the Potosi Dynamic Model 2013a (43 MT in 30 years), since the permeability model applied in the Potosi Dynamic Model 2013b is more conservative. This revision was deemed necessary to treat the uncertainty in a more appropriate manner. As the CO2 follows the paths where vugs interconnection exists, a reasonably large and irregular plume extent was created. For the Optimistic Case, the plume extends 17 miles (27.4km) in E-W and 14 miles (22.5km) in N-S directions after 30 years. After injection is completed, the plume continues to migrate laterally, mainly driven by the remaining pressure gradient. After 100 years post injection, the plume extends 20 miles (32.2km) in E-W and 15.5 miles (24.9km) in N-S directions. Should the targeted cumulative injection of 96 MT be achieved; a much larger plume extent could be expected. For the Optimistic Case, the increase of reservoir pressure at the end of injection is approximately 1200 psia (8,274 kPa) around the injector and gradually decreases away from the well. The reservoir pressure increase is less than 30 psia (206.8 kPa) beyond 14 miles (22.5km) away from injector. Should the targeted cumulative injection of 96 MT be achieved; a much larger areal pressure increase could be expected. The initial reservoir pressure is nearly restored after approximately 100 years post injection. The presence of matrix slows down the pressure dissipations. The Pessimistic Case gives an average CO2 injection rate of 0.2 MTPA and cumulative injection of 7 MT in 30 years, which corresponds to 7% of the injection target. This implies that in the worst case scenario, a minimum of sixteen (16) wells could be required to achieve the injection target. The present evaluation is mainly associated with uncertainty on the vugs permeability, distribution, and interconnectivity. The different results indicated by the Optimistic and Pessimistic Cases signify the importance of vugs permeability characterization. Therefore, injection test and pressure interference test among the wells could be considered to evaluate the local vugs permeability, extent, and interconnectivity. Porosity mapping derived from the seismic inversion could also be used in the succeeding task to characterize the lateral porosity distribution within the reservoir. With or without seismic inversion porosity mapping, it is worth exploring whether increased lateral heterogeneity plays a significant role in Potosi injectivity. Investigations on vugular, dolomitic outcrops suggest that there may be significantly greater lateral heterogeneity than what has been modeled here. Facies modeling within the Potosi has yet to be thoroughly addressed. The carbonates during the time of deposition are believed to be regionally extensive. However, it may be worth delineating the reservoir with other regional wells or modern day analogues to understand the extent of the Potosi. More specifically, the model could incorporate lateral changes or trends if deemed necessary to represent facies transition. Data acquisitions to characterize the fracture pressure gradient, the formation water properties, the relative permeability, and the capillary pressure could also be considered in order to allow a more rigorous evaluation of the Potosi storage performance. A simulation using several injectors could also be considered to determine the required number of wells to achieve the injection target while taking into account the pressure interference.

Adushita, Yasmin↗
DOE OSTI2021

Determining fast-S and slow-S propagation directions with SV-P data produced by buried explosives and recorded with vertical geophones

We have evaluated the concept of practicing S-wave reflection seismology with legacy 3D seismic data generated by a P-wave source and recorded with only vertical geophones. This type of S-wave imaging is based on the principle that seismic P-wave sources not only produce a downgoing illuminating P wavefield, but they also simultaneously produce a downgoing illuminating SV wavefield that, in almost all cases, is suitable for S-wave reflection imaging. The S-mode used in this study is the SV-P, or converted-P, mode. This mode involves a downgoing illuminating SV wavefield and an upgoing reflected P-mode that is recorded by vertical geophones. In flat-layered stratigraphy, the lengths of the SV and P raypaths in SV-P imaging are identical to the lengths of the SV and P raypaths in P-SV imaging with P-sources and 3C geophones. P-SV imaging of deep rocks has been practiced for more than two decades; SV-P imaging is a new concept. SV-P data should provide the same options for investigating deep rocks as do P-SV data. We have determined one of the equivalences between SV-P data extracted from vertical-geophone data and P-SV data extracted from horizontal geophones: that both modes react to azimuth-dependent variations in the S velocity in anisotropic rocks. Azimuthal variations in the SV-P traveltime can be used to define the polarization direction of the fast-S-wave mode, which is also the azimuth of the maximum horizontal stress (SHmax). Our investigation demonstrates a noninvasive method for monitoring changes in the SHmax azimuth across a CO 2 storage reservoir, or any targeted porous rock, as fluids are cycled into, and then out of, that rock’s pore space.

Geochemistry & Geophysics↗
DOE OSTI2014

The Potosi Reservoir Model 2013c, Property Modeling Update

As part of a larger project co-funded by the United States Department of Energy (US DOE) to evaluate the potential of formations within the Cambro-Ordovician strata above the Mt. Simon as potential targets for carbon sequestration in the Illinois and Michigan Basins, the Illinois Clean Coal Institute (ICCI) requested Schlumberger to evaluate the potential injectivity and carbon dioxide (CO2) plume size of the Cambrian Potosi Formation. The evaluation of this formation was accomplished using wireline data, core data, pressure data, and seismic data from this project as well as two other separately funded projects: the US DOE-funded Illinois Basin–Decatur Project (IBDP) being conducted by the Midwest Geological Sequestration Consortium (MGSC) in Macon County, Illinois, and the Illinois Industrial Carbon Capture and Sequestration (ICCS) project funded through the American Recovery and Reinvestment Act. In 2010, technical performance evaluations on the Cambrian Potosi Formation were performed through reservoir modeling. The data included formation tops from mud logs, well logs from the Verification Well #1 (VW1) and the Injection Well (CCS1), structural and stratigraphic formation from three dimensional (3D) seismic data, and field data from several waste water injection wells for Potosi Formation. The intention was for 2.2 million tons per annum (2 million tonnes per annum [MTPA]) of CO2 to be injected for 20 years. In the Task Error! Reference source not found., the 2010 Potosi heterogeneous model (referred to as the "Potosi Dynamic Model 2010") was re-run using a new injection scenario of 3.5 million tons per annum (3.2 MTPA) for 30 years. The extent of the Potosi Dynamic Model 2010, however, appeared too small for the new injection target. The models size was insufficient to accommodate the evolution of the plume. The new model, Potosi Dynamic Model 2013a, was built by extending the Potosi Dynamic Model 2010 grid to 30 by 30 mi (48 by 48 km), while preserving all property modeling workflows and layering. This model was retained as the base case. In the preceding Task [1], the Potosi reservoir model was updated to take into account the new data from the Verification Well #2 (VW2) which was drilled in 2012. The porosity and permeability modeling was revised to take into account the log data from the new well. Revisions of the 2010 modeling assumptions were also done on relative permeability, capillary pressures, formation water salinity, and the maximum allowable well bottomhole pressure. Dynamic simulations were run using the injection target of 3.5 million tons per annum (3.2 MTPA) for 30 years. This dynamic model was named Potosi Dynamic Model 2013b. In this Task, a new property modeling workflow was applied, where seismic inversion data guided the porosity mapping and geobody extraction. The static reservoir model was fully guided by PorosityCube interpretations and derivations coupled with petrophysical logs from three wells. The two main assumptions are: porosity features in the PorosityCube that correlate with lost circulation zones represent vugular zones, and that these vugular zones are laterally continuous. Extrapolation was done carefully to populate the vugular facies and their corresponding properties outside the seismic footprint up to the boundary of the 30 by 30 mi (48 by 48 km) model. Dynamic simulations were also run using the injection target of 3.5 million tons per annum (3.2 MTPA) for 30 years. This new dynamic model was named Potosi Dynamic Model 2013c. Reservoir simulation with the latest model gives a cumulative injection of 43 million tons (39 MT) in 30 years with a single well, which corresponds to 40% of the injection target. The injection rate is approx. 3.2 MTPA in the first six months as the well is injecting into the surrounding vugs, and declines rapidly to 1.8 million tons per annum (1.6 MTPA) in year 3 once the surrounding vugs are full and the CO2 start to reach the matrix. After, the injection rate declines gradually to 1.2 million tons per annum (1.1 MTPA) in year 18 and stays constant. This implies that a minimum of three (3) wells could be required in the Potosi to reach the injection target. The injectivity evaluated in this Task was higher compared to the preceding Task, since the current facies modeling (guided by the porosity map from the seismic inversion) indicated a higher density of vugs within the vugular zones. 5 As the CO2 follows the paths where vugs interconnection exists, a reasonably large and irregular plume extent was created. After 30 years of injection, the plume extends 13.7 mi (22 km) in E-W and 9.7 mi (16 km) in N-S directions. After injection finishes, the plume continues to migrate laterally, mainly driven by the remaining pressure gradient. After 60 years post-injection, the plume extends 14.2 mi (22.8 km) in E-W and 10 mi (16 km) in N-S directions, and remains constant as the remaining pressure gradient has become very low. Should the targeted cumulative injection of 106 million tons (96 MT) be achieved; a much larger plume extent could be expected. The increase of reservoir pressure at the end of injection is approximately 1,200 psia (8,274 kPa) around the injector and gradually decreases away from the well. The reservoir pressure increase is less than 10 psia (69 kPa) beyond 14 mi (23 km) away from injector. Should the targeted cumulative injection of 106 million tons (96 MT) be achieved; a much larger areal pressure increase could be expected. The reservoir pressure declines rapidly during the first 30 years post injection and the initial reservoir pressure is nearly restored after 100 years post-injection. The present evaluation is mainly associated with uncertainty on the vugs permeability and interconnectivity. The use of porosity mapping from seismic inversion might have reduced the uncertainty on the lateral vugs body distributions. However, major uncertainties on the Potosi vugs permeability remains. Therefore, injection test and pressure interference test among the wells could be considered to evaluate the local vugs permeability, extent, and interconnectivity. Facies modeling within the Potosi has yet to be thoroughly addressed. The carbonates during the time of deposition are believed to be regionally extensive. However, it may be worth delineating the reservoir with other regional wells or modern day analogues to understand the extent of the Potosi. More specifically, the model could incorporate lateral changes or trends if deemed necessary to represent facies transition. Data acquisitions to characterize the fracture pressure gradient, the formation water properties, the relative permeability, and the capillary pressure could also be considered in order to allow a more rigorous evaluation of the Potosi storage performance. A simulation using several injectors could also be considered to determine the required number of wells and appropriate spacing to achieve the injection target while taking into account the pressure interference.

Adushita, Yasmin↗