Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir properties”

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

Deep Direct-Use Feasibility Study Estimation of Reservoir Properties for the Tuscarora Sandstone

This dataset contains black-and-white photographic images of individual core segments from well Preston 119 and photomicrographs taken of petrographic thin sections from well Clay 513. Minipermeameter measurement results of 2,260 samples for fracture and matrix permeability on core segments for well Preston 119. These results are summarized in figures, spreadsheets, and PowerPoints. Microsoft Excel will be necessary to examine the permeability data and software capable of reading a .PDF file will be needed to view the file of annotated core CT scans. Microsoft Powerpoint will be necessary to view the photomicrograph presentation.

15 GEOTHERMAL ENERGY↗

Reservoir Properties of Barremian – Aptian Urgonian Limestones, Se France, Part 1: Influence of Structural History on Porosity-Permeability Variations

Upper Barremian – Lower Aptian inner platform “Urgonian” limestones in the Mont de Vaucluse region, SE France, consist of alternating metre-scale microporous and tight intervals. This paper focuses on the influence of structural deformation on the reservoir properties of the Urgonian limestone succession in a study area near the town of Rustrel. Petrographic, petrophysical and structural data were recovered from five fully-cored boreholes, from the walls of a 100 m long underground tunnel, and from a 50 m long transect at a nearby outcrop. The data allowed reservoir property variations in the Urgonian limestones to be studied from core to reservoir scale. Eleven Reservoir Rock Types (RRTs) were identified based on petrographic features (texture, grain size), reservoir properties (porosity, permeability), and the frequency of structural discontinuities such as fractures, faults and stylolites. Tight and microporous reservoir rock types were distinguished. Tight reservoir rock types were characterised by early cementation of intergranular pore spaces and by the presence of frequent structural discontinuities. By contrast microporous reservoir rock types contained preserved intragranular microporosity and matrix permeability, but had very few structural discontinuities. Observed vertical alternations of microporous and tight rock types are interpreted to have been controlled by the early diagenesis of the Urgonian carbonates. Deformation associated with regional-scale tectonic phases, including Albian – Cenomanian “Durancian” uplift (~105 to 96 Ma) and Pyrenean compression (~55 to 25 Ma), resulted in the modification of the initial petrophysical properties of the Urgonian limestones. An early diagenetic imprint conditioned both the intensity of structural deformations and the associated circulations of diagenetic and meteoric fluids. Evolution of the Reservoir Rock Types is therefore linked both to the depositional conditions and to subsequent phases of structural deformation.

Cochard, J.↗

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

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

15 GEOTHERMAL ENERGY↗

Oil Reservoir Properties Estimation Using Neutal Networks

This paper investigates the applicability as well as the accuracy of artificial neural networks for estimating specific parameters that describe reservoir properties based on seismic data. Our approach relies on JPL's adjoint operators general purpose neural network code to determine the best suited architecture. We believe that results presented in this work demonstrate that artificial neural networks produce surprisingly accurate estimates of the reservoir parameters.

oil reservoir neural networks geoscience data oil ↗

Numerical simulation of gas production from natural gas hydrate deposits with multi-branch wells: Influence of reservoir properties

Vast amounts of natural gas hydrate are buried in subseafloor sediments without impermeable boundaries, which is recognized as an essential energy source for the future. Previously, the multi-branch well was proposed to enhance the recovery efficiency of natural gas hydrate, and the gas production rate has been dramatically improved comparing with the vertical well. However, the multi-branch well shows a terrible performance in gas production duration. Here, as a continuation of the previous study, numerical simulations were conducted to investigate the influence of hydrate reservoir properties on the gas production potential. Results indicate that it is hard to extract hydrate commercially for hydrate accumulations without impermeable boundaries. A high initial hydrate saturation leads to a long gas production duration but a low gas production rate. An increase in the intrinsic permeability of isotropic reservoirs would shorten the gas production duration and result in a low gas recovery ratio. Permeability anisotropy shows a noticeable effect on enhancing the gas recovery ratio and the gas production duration due to the improved pressure propagation pattern. Therefore, in the upcoming field tests, reservoir reconstructions that enhance permeability anisotropy are strongly suggested to obtain better outcomes.

03 NATURAL GAS↗

Inversion of Time-Lapse Seismic Reservoir Monitoring Data Using CycleGAN: A Deep Learning-Based Approach for Estimating Dynamic Reservoir Property Changes

Carbon capture and storage is being pursued globally as a geoengineering measure for reducing the emission of anthropogenic CO 2 the atmosphere. Comprehensive monitoring, verification, and accounting programs must be established for demonstrating the safe storage of injected CO 2 . One of the most commonly deployed monitoring techniques is time-lapse seismic reservoir monitoring (also known as 4-D seismic), which involves comparing 3-D seismic survey data taken at the same study site but over different times. Analyses of 4-D seismic data volumes can help improve the quality of storage reservoir characterization, track the movement of injected CO 2 plume, and identify potential CO 2 spillover/leakage from the storage reservoirblue. However, the derivation of high-resolution CO 2 saturation maps from 4-D seismic data is a highly nonlinear and ill-posed inverse problem, often requiring significant computational effort. In this research, we apply a physics-based deep learning method to facilitate the solution of both the forward and inverse problems in seismic inversion while honoring physical constraints. A cycle generative adversarial neural network (CycleGAN) model is trained to learn the bidirectional functional mappings between the reservoir dynamic property changes and seismic attribute changes, such that both forward and inverse solutions can be obtained efficiently from the trained model. We show that our CycleGAN-based approach not only improves the reliability of 4-D seismic inversion but also expedites the quantitative interpretation. Our deep learning-based workflow is generic and can be readily used for reservoir characterization and reservoir model updates involving the use of 4-D seismic data.

58 GEOSCIENCES↗

Controls on reservoir properties in organic-matter-rich shales: Insights from MICP analysis

Mercury injection capillary pressure (MICP) analysis was used to study organic-matter-rich Devonian shales varying in maturity from marginally mature to postmature to discuss controls on drainage and imbibition in shale reservoirs. Six samples come from the New Albany Shale; their total organic carbon (TOC) ranged from 3.3–8.4 % and vitrinite reflectance (R o ) from 0.54–1.42 %. One Marcellus Shale sample with R o 2.50 % was used to extend the maturity range to the advanced dry gas window. Three coal samples of varying rank were also included for comparison. Porosity of the shales ranges from 1.5–5.0 %, and MICP-derived permeability ranges from 0.17–9.9 mD. Injection capillary pressure curves have a range of entry pressures from 0.055 to0.421 MPa, and no relationship of entry pressure with maturity was revealed. Maturity, however, has an influence on pore size distribution, with the higher-maturity samples having significantly lower volumes of pores with throats less than 10 nm in diameter. Withdrawal capillary pressure curves document a wide range of withdrawal efficiencies from 55 % to almost 100 %, and pores having throats between 3–50 nm in diameter play a major role in the ability of mercury to be released from the pore system. The obtained data provide valuable information about shale properties and their behavior not only in the air-mercury system but also under reservoir conditions.

02 PETROLEUM↗

Simple Statistical Models for Predicting Overpressure Due to CO2 and Low-Salinity Waste-Fluid Injection into Deep Saline Formations

Deep saline aquifers have been used for waste-fluid disposal for decades and are the proposed targets for large-scale CO2 storage to mitigate CO2 concentration in the atmosphere. Due to relatively limited experience with CO2 injection in deep saline formations and given that the injection targets for CO2 sometimes are the same as waste-fluid disposal formations, it could be beneficial to model and compare both practices and learn from the waste-fluid disposal industry. In this paper, we model CO2 injection in the Patterson Field, which has been proposed as a site for storage of 50 Mt of industrial CO2 over 25 years. We propose general models that quickly screen the reservoir properties and calculate pressure changes near and far from the injection wellbore, accounting for variable reservoir properties. The reservoir properties we investigated were rock compressibility, injection rate, vertical-to-horizontal permeability ratio, average reservoir permeability and porosity, reservoir temperature and pressure, and the injectant total dissolved solids (TDS) in cases of waste-fluid injection. We used experimental design to select and perform simulation runs, performed a sensitivity analysis to identify the important variables on pressure build-up, and then fit a regression model to the simulation runs to obtain simple proxy models for changes in average reservoir pressure and bottomhole pressure. The CO2 injection created more pressure compared to saline waste-fluids, when similar mass was injected. However, we found a more significant pressure buildup at the caprock-reservoir interface and lower pressure buildup at the bottom of the reservoir when injecting CO2 compared with waste-fluid injection.

Ansari, Esmail↗

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗

Illitization in the Mt. Simon Sandstone, Illinois Basin, USA: Implications for carbon dioxide storage

Clay mineral cementation is one of the most important controls on sandstone reservoir properties. Here, in this paper, the diagenetic history of the Mt. Simon reservoir complex is studied to reveal the origin, timing, and controls of illitization. Samples of sandstone and shale from Mt. Simon reservoir complex were acquired from the Illinois Basin–Decatur Project (IBDP), a CO 2 storage demonstration project in the central Illinois Basin. Petrographic, SEM, TEM, XRD analyses and K/Ar age dating were completed to identify the major detrital and diagenetic components of the samples and reveal illite to be the major clay component in all samples. Illitic clay coatings in the lower Mt. Simon reservoir are identified as a major control on reservoir properties by inhibiting major precipitation of authigenic quartz during illitization, resulting in highly permeable sandstone, essential for CO 2 storage. The coating box-work morphology and mineralogy are indicative of a detrital smectite origin with subsequent illite growth associated with feldspar dissolution and kaolinite alteration. The mineralogy of bulk material and separate grain size fractions (2–0.6 μm; 0.6–0.2 μm; < 0.2 μm) with illite polytypes 2M 1 , 1M, and 1M d were quantified and age dated via 40 K- 40 Ar methods. The shales or reservoir seals contain the highest proportions of detrital illite with illite in the lower Mt. Simon Sandstone reservoir identified as solely diagenetic. Two major events of illitization are identified throughout the Mt. Simon with more porous reservoir rock exhibiting the older event from approximately 360 to 315 Ma and tighter sandstone units with reservoir properties too low to be considered reservoir exhibiting illite dates from 250 to 220 Ma. This partitioning of illitization is attributed to varied mineralogy controlled by depositional changes and evolution of the greater basin.

58 GEOSCIENCES↗

Machine learning for geophysical characterization of brittleness: Tuscaloosa Marine Shale case study

Brittleness is one of the most important reservoir properties for unconventional reservoir exploration and production. Better knowledge about the brittleness distribution can help to optimize the hydraulic fracturing operation and lower costs. However, there are very few reliable and effective physical models to predict the spatial distribution of brittleness. We have developed a machine learning-based method to predict subsurface brittleness by using multidiscipline data sets, such as seismic attributes, rock physics, and petrophysics information, which allows us to implement the prediction without using a physical model. The method is applied on a data set from Tuscaloosa Marine Shale, and the predicted rock physics template is close to the calculated value from conventional inverted elastic parameters. Therefore, the proposed method helps determine areas of the reservoir that have optimal geomechanical properties for successful hydraulic fracturing.

Geochemistry & Geophysics↗

Data Arrays for Microearthquake (MEQ) Monitoring using Deep Learning for the Newberry EGS Sites

The 'Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties' project looks to apply machine learning (ML) methods to Microearthquake (MEQ) data for imaging geothermal reservoir properties and forecasting seismic events, in order to advance geothermal exploration and safe geothermal energy production. As part of the project, this submission provides data arrays for 149 microearthquakes between the year 2012 and 2013 at the Newberry EGS Site for use with the Deep Learning Algorithm that has been developed. The data provided includes raw waveform data, location data, normalized waveform data, and processed waveform data. Penn State Geothermal Team has shared the following files from the project: - 149 microearthquakes (MEQs) between 2012 and 2013 at Newberry EGS sites, 'Normalized Waveform Inputs.npz' are normalized waveforms. - labels of 149 MEQs: Processed Waveform Inputs.npz - location labels of 149 MEQs: Location Data.npz Note: .npz is the python file format by NumPy that provides storage of array data.

15 GEOTHERMAL ENERGY↗

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↗

Evaluating Probability of Containment Effectiveness at a GCS Sites using integrated assessment modeling approach with Bayesian decision Networks

Improved scientific and engineering understanding of the behavior of geologic CO2 storage together with established regulatory framework and incentive structures raise the prospects for accelerated, large-scale deployment of this greenhouse gas emissions reduction approach. Incentive structures call for the establishment of appropriate verification and accounting approaches to support claims of the integrity of a geologic storage complex and to justify taking credit for long-term storage. In this study, we present a framework for assessing the probability of containment effectiveness over the lifetime of a geologic carbon storage site (e.g., after 70 years of injection and post-injection site performance) using forward stochastic model realizations based on site characterization data and using a monitoring-informed Bayesian network based on hypothetical detectability from surface seismic surveys over the site injection and post-injection phases. The National Risk Assessment Partnership’s open-source Integrated Assessment Model (NRAP-Open-IAM) was utilized to develop an ensemble of 10,000 a priori stochastic forecasts of CO2 containment. Those simulations were used to train the Bayesian network model to estimate the prior probabilities of the CO2 leakage mass into overlying, monitorable aquifers considering the uncertainties in the reservoir properties, permeability of potentially leaky wells and the overlying aquifers. The conditional probabilities in the Bayesian network were either learned from the NRAP-Open-IAM simulations or derived from the predefined detection thresholds for the monitoring method. Observations obtained from monitoring, over time during the site operation phases were then used to generate updated posterior probabilities of containment (and any loss from containment) in the Bayesian network by propagating the prior probabilities through the conditional probabilities. We demonstrate how to construct and use the Bayesian network for verifying the long-term storage complex effectiveness informed by monitoring based on the NRAP-Open-IAM simulations previously developed for the FutureGen 2.0 site. This approach may have relevance for stake holders to demonstrate secure geologic storage, provide a defensible, probabilistic approach to claim credit for geologic storage, and to estimate the likelihood that any fraction of the claimed credit may need to be refunded to the creditor based on available monitoring information.

Bayesian network, Risk assessment, Monitoring, car↗

Assessment of Enhanced Oil Recovery Using Carbon Dioxide in Michigan Basin Silurian Pinnacle Reefs

The Michigan team at Western Michigan University has provided technical analysis and support of the enhanced oil recovery (EOR) using carbon dioxide (CO2) in the Northern Michigan Silurian Pinnacle reefs. This report discussed the relevant data about reservoir properties, including lithologic and depositional characteristics of the reservoir formations and porosity and fluid flow properties through different compartments in the reservoir. The team has also developed models that represent the geological and physical characteristics of the reservoir and seal system.

EOR↗