Rock-physics machine learning toolkit for joint litho-fluid facies classification and compaction modeling

2021 ◽  
Vol 40 (10) ◽  
pp. 742-750
Author(s):  
Roman Beloborodov ◽  
James Gunning ◽  
Marina Pervukhina ◽  
Kester Waters ◽  
Nick Huntbatch
Keyword(s):  
Machine Learning ◽  
Em Algorithm ◽  
Prior Information ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Borehole Data ◽  
Facies Classification ◽  
Physics Model ◽  
Rock Physics Model ◽  
Fitting Parameters

Correct lithofacies interpretation sourced from wireline log data is an essential source of prior information for joint seismic inversion for facies and impedances, among other applications. However, this information is difficult to interpret or extract manually due to the multivariate and high dimensionality of wireline logs. Facies inference is also challenging for traditional clustering-based approaches because pervasive compaction trends affect a number of petrophysical measurements simultaneously. Another common pitfall in automated clustering approaches is the inability to account for underlying diagenetic processes that correlate with depth. Here, we address these challenges by introducing a rock-physics machine learning toolkit for joint litho-fluid facies classification. The litho-fluid types are inferred from the borehole data within the objective framework of a maximum-likelihood approach for latent facies variables and rock-physics model parameters, explicitly accounting for compaction and depth effects. The inference boils down to an expectation-maximization (EM) algorithm with strong spatial coupling. Each litho-fluid type is associated with an instance of a particular rock-physics model with a unique set of fitting parameters, constrained to a physically reasonable range. These fitting parameters in turn are inferred using bound-constrained optimization as part of the EM algorithm. Outputs produced by the toolkit can be used directly to specify the necessary prior information for seismic inversion, including per-facies rock-physics models and facies proportions. We present an example application of the tool to real borehole data from the North West Shelf of Australia to illustrate the method and discuss its characteristic features in depth.

Rock physics model‐based seismic inversion

2006 ◽  
Author(s):  
Kyle Spikes ◽  
Jack Dvorkin ◽  
Gary Mavko
Keyword(s):  
Rock Physics ◽  
Seismic Inversion ◽  
Model Based ◽  
Physics Model ◽  
Rock Physics Model

Facies — The drivers for modern inversions

2020 ◽  
Vol 39 (2) ◽  
pp. 102-109
Author(s):  
John Pendrel ◽  
Henk Schouten
Keyword(s):  
Bayesian Analysis ◽  
Prior Information ◽  
Reservoir Characterization ◽  
Low Frequency ◽  
Rock Physics ◽  
Well Log ◽  
Data Set ◽  
Prior Probabilities ◽  
Physics Model ◽  
Rock Physics Model

It is common practice to make facies estimations from the outcomes of seismic inversions and their derivatives. Bayesian analysis methods are a popular approach to this. Facies are important indicators of hydrocarbon deposition and geologic processes. They are critical to geoscientists and engineers. The application of Bayes’ rule maps prior probabilities to posterior probabilities when given new evidence from observations. Per-facies elastic probability density functions (ePDFs) are constructed from elastic-log and rock-physics model crossplots, over which inversion results are superimposed. The ePDFs are templates for Bayesian analysis. In the context of reservoir characterization, the new information comes from seismic inversions. The results are volumes of the probabilities of occurrences of each of the facies at all points in 3D space. The concepts of Bayesian inference have been applied to the task of building low-frequency models for seismic inversions without well-log interpolation. Both a constant structurally compliant elastic trend approach and a facies-driven method, where models are constructed from per-facies trends and initial facies estimates, have been tested. The workflows make use of complete 3D prior information and measure and account for biases and uncertainties in the inversions and prior information. Proper accounting for these types of effects ensures that rock-physics models and inversion data prepared for reservoir property analysis are consistent. The effectiveness of these workflows has been demonstrated by using a Gulf of Mexico data set. We have shown how facies estimates can be effectively used to build reasonable low-frequency models for inversion, which obviate the need for well-log interpolation and provide full 3D variability. The results are more accurate probability-based net-pay estimates that correspond better to geology. We evaluate the workflows by using several measures including precision, confidence, and probabilistic net pay.

Two-parameter prestack seismic inversion of porosity and pore structure parameter of fractured carbonate reservoirs: Part 1 — Methods

2018 ◽  
Vol 6 (4) ◽  
pp. SM1-SM8 ◽  
Author(s):  
Tingting Zhang ◽  
Yuefeng Sun
Keyword(s):  
Pore Structure ◽  
Seismic Data ◽  
Rock Physics ◽  
Seismic Inversion ◽  
High Porosity ◽  
Structure Parameter ◽  
Pore Structure Parameter ◽  
Fractured Zones ◽  
Physics Model ◽  
Rock Physics Model

Fractured zones in deeply buried carbonate hills are important because they often have better permeability resulting in prolific production than similar low-porosity rocks. Nevertheless, their detection poses great challenge to conventional seismic inversion methods because they are mostly low in acoustic impedance and bulk modulus, hardly distinguishable from high-porosity zones or mudstones. A proxy parameter of pore structure defined in a rock-physics model, the so-called Sun model, has been used for delineating fractured zones in which the pore structure parameter is relatively high, whereas the porosity is low in general. Simultaneous seismic inversion of the pore structure parameter and porosity proves to be difficult and nontrivial in practice. Although the pore structure parameter is well-defined at locations where density, P-, and S-velocity are known from logs, estimation of P- and S-velocity information, especially density information from prestack seismic data is rather challenging. A three-step iterative inversion method, which uses acoustic, gradient, and elastic impedance from angle-stacked seismic data as input to the rock-physics model for calculating porosity and bulk and shear pore structure parameters simultaneously, is proposed and implemented to solve this problem. The methodology is successfully tested with well logs and seismic data from a deeply buried carbonate hill in the Bohai Bay Basin, China.

A predictive anisotropic rock-physics model for estimating elastic rock properties of unconventional shale reservoirs

2019 ◽  
Vol 38 (5) ◽  
pp. 358-365 ◽  
Author(s):  
Colin M. Sayers ◽  
Sagnik Dasgupta
Keyword(s):  
Accurate Determination ◽  
Rock Physics ◽  
Seismic Inversion ◽  
P Wave ◽  
Rock Properties ◽  
Anisotropic Rock ◽  
Rock Classification ◽  
Shale Reservoirs ◽  
Physics Model ◽  
Rock Physics Model

This paper presents a predictive rock-physics model for unconventional shale reservoirs based on an extended Maxwell scheme. This model accounts for intrinsic anisotropy of rock matrix and heterogeneities and shape-induced anisotropy arising because the dimensions of kerogen inclusions and pores are larger parallel to the bedding plane than perpendicular to this plane. The model relates the results of seismic amplitude variation with offset inversion, such as P- and S-impedance, to the composition of the rock and enables identification of rock classes such as calcareous, argillaceous, siliceous, and mixed shales. This allows the choice of locations with the best potential for economic production of hydrocarbons. While this can be done using well data, prestack inversion of seismic P-wave data allows identification of the best locations before the wells are drilled. The results clearly show the ambiguity in rock classification obtained using poststack inversion of P-wave seismic data and demonstrate the need for prestack seismic inversion. The model provides estimates of formation anisotropy, as required for accurate determination of P- and S-impedance, and shows that anisotropy is a function not only of clay content but also other components of the rock as well as the aspect ratio of kerogen and pores. Estimates of minimum horizontal stress based on the model demonstrate the need to identify rock class and estimate anisotropy to determine the location of any stress barriers that may inhibit hydraulic fracture growth.

scholarly journals Seismic characterization of geothermal sedimentary reservoirs: A field example from the Copenhagen area, Denmark

2020 ◽  
Vol 8 (2) ◽  
pp. T275-T291 ◽  
Author(s):  
Kenneth Bredesen ◽  
Esben Dalgaard ◽  
Anders Mathiesen ◽  
Rasmus Rasmussen ◽  
Niels Balling
Keyword(s):  
Reservoir Characterization ◽  
Energy Resource ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Amplitude Variation With Offset ◽  
Seismic Characterization ◽  
Rock Physics Modeling ◽  
Physics Modeling ◽  
Physics Model ◽  
Rock Physics Model

We have seismically characterized a Triassic-Jurassic deep geothermal sandstone reservoir north of Copenhagen, onshore Denmark. A suite of regional geophysical measurements, including prestack seismic data and well logs, was integrated with geologic information to obtain facies and reservoir property predictions in a Bayesian framework. The applied workflow combined a facies-dependent calibrated rock-physics model with a simultaneous amplitude-variation-with-offset seismic inversion. The results suggest that certain sandstone distributions are potential aquifers within the target interval, which appear reasonable based on the geologic properties. However, prediction accuracy suffers from a restricted data foundation and should, therefore, only be considered as an indicator of potential aquifers. Despite these issues, the results demonstrate new possibilities for future seismic reservoir characterization and rock-physics modeling for exploration purposes, derisking, and the exploitation of geothermal energy as a green and sustainable energy resource.

scholarly journals Monitoring geological storage of CO2 using a new rock physics model

2022 ◽  
Vol 12 (1) ◽  
Author(s):  
Manzar Fawad ◽  
Nazmul Haque Mondol
Keyword(s):  
Waste Product ◽  
Rock Physics ◽  
Time Lapse ◽  
Seismic Inversion ◽  
Storage Site ◽  
Physics Model ◽  
Rock Physics Model ◽  
Pressure Changes ◽  
Prestack Seismic Inversion ◽  
Migration Pathways

AbstractTo mitigate the global warming crisis, one of the effective ways is to capture CO2 at an emitting source and inject it underground in saline aquifers, depleted oil and gas reservoirs, or in coal beds. This process is known as carbon capture and storage (CCS). With CCS, CO2 is considered a waste product that has to be disposed of properly, like sewage and other pollutants. While and after CO2 injection, monitoring of the CO2 storage site is necessary to observe CO2 plume movement and detect potential leakage. For CO2 monitoring, various physical property changes are employed to delineate the plume area and migration pathways with their pros and cons. We introduce a new rock physics model to facilitate the time-lapse estimation of CO2 saturation and possible pressure changes within a CO2 storage reservoir based on physical properties obtained from the prestack seismic inversion. We demonstrate that the CO2 plume delineation, saturation, and pressure changes estimations are possible using a combination of Acoustic Impedance (AI) and P- to S-wave velocity ratio (Vp/Vs) inverted from time-lapse or four-dimensional (4D) seismic. We assumed a scenario over a period of 40 years comprising an initial 25 year injection period. Our results show that monitoring the CO2 plume in terms of extent and saturation can be carried out using our rock physics-derived method. The suggested method, without going into the elastic moduli level, handles the elastic property cubes, which are commonly obtained from the prestack seismic inversion. Pressure changes quantification is also possible within un-cemented sands; however, the stress/cementation coefficient in our proposed model needs further study to relate that with effective stress in various types of sandstones. The three-dimensional (3D) seismic usually covers the area from the reservoir's base to the surface making it possible to detect the CO2 plume's lateral and vertical migration. However, the comparatively low resolution of seismic, the inversion uncertainties, lateral mineral, and shale property variations are some limitations, which warrant consideration. This method can also be applied for the exploration and monitoring of hydrocarbon production.

Sign in / Sign up

Export Citation Format

Share Document