Improving seismic QP estimation using rock-physics constraints

Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. MR187-MR198 ◽  
Author(s):  
Yi Shen ◽  
Jack Dvorkin ◽  
Yunyue Li
Keyword(s):  
Seismic Data ◽  
Closed Form Solution ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Analytical Relation ◽  
Gas Reservoir ◽  
Seismic Section ◽  
Data Set ◽  
Attenuation Model ◽  
Text Model

Our goal is to accurately estimate attenuation from seismic data using model regularization in the seismic inversion workflow. One way to achieve this goal is by finding an analytical relation linking [Formula: see text] to [Formula: see text]. We derive an approximate closed-form solution relating [Formula: see text] to [Formula: see text] using rock-physics modeling. This relation is tested on well data from a clean clastic gas reservoir, of which the [Formula: see text] values are computed from the log data. Next, we create a 2D synthetic gas-reservoir section populated with [Formula: see text] and [Formula: see text] and generate respective synthetic seismograms. Now, the goal is to invert this synthetic seismic section for [Formula: see text]. If we use standard seismic inversion based solely on seismic data, the inverted attenuation model has low resolution and incorrect positioning, and it is distorted. However, adding our relation between velocity and attenuation, we obtain an attenuation model very close to the original section. This method is tested on a 2D field seismic data set from Gulf of Mexico. The resulting [Formula: see text] model matches the geologic shape of an absorption body interpreted from the seismic section. Using this [Formula: see text] model in seismic migration, we make the seismic events below the high-absorption layer clearly visible, with improved frequency content and coherency of the events.

Integrated VTI model building with seismic data, geologic information, and rock-physics modeling — Part 1: Theory and synthetic test

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. C177-C191 ◽  
Author(s):  
Yunyue Li ◽  
Biondo Biondi ◽  
Robert Clapp ◽  
Dave Nichols
Keyword(s):  
Seismic Data ◽  
Seismic Anisotropy ◽  
Model Building ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Additional Information ◽  
Synthetic Test ◽  
Rock Physics Modeling ◽  
Physics Modeling ◽  
Geologic Information

Seismic anisotropy plays an important role in structural imaging and lithologic interpretation. However, anisotropic model building is a challenging underdetermined inverse problem. It is well-understood that single component pressure wave seismic data recorded on the upper surface are insufficient to resolve a unique solution for velocity and anisotropy parameters. To overcome the limitations of seismic data, we have developed an integrated model building scheme based on Bayesian inference to consider seismic data, geologic information, and rock-physics knowledge simultaneously. We have performed the prestack seismic inversion using wave-equation migration velocity analysis (WEMVA) for vertical transverse isotropic (VTI) models. This image-space method enabled automatic geologic interpretation. We have integrated the geologic information as spatial model correlations, applied on each parameter individually. We integrate the rock-physics information as lithologic model correlations, bringing additional information, so that the parameters weakly constrained by seismic are updated as well as the strongly constrained parameters. The constraints provided by the additional information help the inversion converge faster, mitigate the ambiguities among the parameters, and yield VTI models that were consistent with the underlying geologic and lithologic assumptions. We have developed the theoretical framework for the proposed integrated WEMVA for VTI models and determined the added information contained in the regularization terms, especially the rock-physics constraints.

scholarly journals Time-Lapse Seismic Monitoring of Onshore Reservoirs in Niger Delta, Field ‘K’ as a Case Study

2017 ◽  
Vol 1 (1) ◽  
pp. 23-27 ◽  
Author(s):  
A. Ogbamikhumi ◽  
T. Tralagba ◽  
E. E. Osagiede
Keyword(s):  
Seismic Data ◽  
Acoustic Impedance ◽  
Niger Delta ◽  
Rock Physics ◽  
Time Lapse ◽  
Seismic Survey ◽  
Impedance Change ◽  
Data Set ◽  
Reservoir Fluid ◽  
4D Seismic

Field ‘K’ is a mature field in the coastal swamp onshore Niger delta, which has been producing since 1960. As a huge producing field with some potential for further sustainable production, field monitoring is therefore important in the identification of areas of unproduced hydrocarbon. This can be achieved by comparing production data with the corresponding changes in acoustic impedance observed in the maps generated from base survey (initial 3D seismic) and monitor seismic survey (4D seismic) across the field. This will enable the 4D seismic data set to be used for mapping reservoir details such as advancing water front and un-swept zones. The availability of good quality onshore time-lapse seismic data for Field ‘K’ acquired in 1987 and 2002 provided the opportunity to evaluate the effect of changes in reservoir fluid saturations on time-lapse amplitudes. Rock physics modelling and fluid substitution studies on well logs were carried out, and acoustic impedance change in the reservoir was estimated to be in the range of 0.25% to about 8%. Changes in reservoir fluid saturations were confirmed with time-lapse amplitudes within the crest area of the reservoir structure where reservoir porosity is 0.25%. In this paper, we demonstrated the use of repeat Seismic to delineate swept zones and areas hit with water override in a producing onshore reservoir.

Common-focus point-based target-oriented imaging approach for continuous seismic reservoir monitoring

Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. M41-M48 ◽  
Author(s):  
Hongwei Liu ◽  
Mustafa Naser Al-Ali
Keyword(s):  
Seismic Data ◽  
Rock Physics ◽  
Time Lapse ◽  
Velocity Estimation ◽  
Estimation Methods ◽  
Data Sets ◽  
Data Set ◽  
Velocity Models ◽  
Reservoir Monitoring ◽  
Focus Point

The ideal approach for continuous reservoir monitoring allows generation of fast and accurate images to cope with the massive data sets acquired for such a task. Conventionally, rigorous depth-oriented velocity-estimation methods are performed to produce sufficiently accurate velocity models. Unlike the traditional way, the target-oriented imaging technology based on the common-focus point (CFP) theory can be an alternative for continuous reservoir monitoring. The solution is based on a robust data-driven iterative operator updating strategy without deriving a detailed velocity model. The same focusing operator is applied on successive 3D seismic data sets for the first time to generate efficient and accurate 4D target-oriented seismic stacked images from time-lapse field seismic data sets acquired in a [Formula: see text] injection project in Saudi Arabia. Using the focusing operator, target-oriented prestack angle domain common-image gathers (ADCIGs) could be derived to perform amplitude-versus-angle analysis. To preserve the amplitude information in the ADCIGs, an amplitude-balancing factor is applied by embedding a synthetic data set using the real acquisition geometry to remove the geometry imprint artifact. Applying the CFP-based target-oriented imaging to time-lapse data sets revealed changes at the reservoir level in the poststack and prestack time-lapse signals, which is consistent with the [Formula: see text] injection history and rock physics.

Characterization of a carbonate geothermal reservoir using rock-physics-guided deep neural networks

2021 ◽  
Vol 40 (10) ◽  
pp. 751-758
Author(s):  
Fabien Allo ◽  
Jean-Philippe Coulon ◽  
Jean-Luc Formento ◽  
Romain Reboul ◽  
Laure Capar ◽  
...  
Keyword(s):  
Neural Networks ◽  
Seismic Data ◽  
Deep Neural Networks ◽  
Synthetic Data ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Geothermal Reservoir ◽  
Rock Properties ◽  
Hybrid Techniques ◽  
Geothermal Wells

Deep neural networks (DNNs) have the potential to streamline the integration of seismic data for reservoir characterization by providing estimates of rock properties that are directly interpretable by geologists and reservoir engineers instead of elastic attributes like most standard seismic inversion methods. However, they have yet to be applied widely in the energy industry because training DNNs requires a large amount of labeled data that is rarely available. Training set augmentation, routinely used in other scientific fields such as image recognition, can address this issue and open the door to DNNs for geophysical applications. Although this approach has been explored in the past, creating realistic synthetic well and seismic data representative of the variable geology of a reservoir remains challenging. Recently introduced theory-guided techniques can help achieve this goal. A key step in these hybrid techniques is the use of theoretical rock-physics models to derive elastic pseudologs from variations of existing petrophysical logs. Rock-physics theories are already commonly relied on to generalize and extrapolate the relationship between rock and elastic properties. Therefore, they are a useful tool to generate a large catalog of alternative pseudologs representing realistic geologic variations away from the existing well locations. While not directly driven by rock physics, neural networks trained on such synthetic catalogs extract the intrinsic rock-physics relationships and are therefore capable of directly estimating rock properties from seismic amplitudes. Neural networks trained on purely synthetic data are applied to a set of 2D poststack seismic lines to characterize a geothermal reservoir located in the Dogger Formation northeast of Paris, France. The goal of the study is to determine the extent of porous and permeable layers encountered at existing geothermal wells and ultimately guide the location and design of future geothermal wells in the area.

Assessment of gas hydrate using pre-stack seismic inversion in Mahanadi Basin, offshore eastern India

2021 ◽  
pp. 1-42
Author(s):  
Maheswar Ojha ◽  
Ranjana Ghosh
Keyword(s):  
Gas Hydrate ◽  
Rock Physics ◽  
Type Equation ◽  
Seismic Inversion ◽  
Log Analysis ◽  
Well Log ◽  
Seismic Section ◽  
Bottom Simulating Reflector ◽  
Hydrate Saturation ◽  
Well Log Analysis

The Indian National Gas Hydrate Program Expedition-01 in 2006 has discovered gas hydrate in Mahanadi offshore basin along the eastern Indian margin. However, well log analysis, pressure core measurements and Infra-Red (IR) anomalies reveal that gas hydrates are distributed as disseminated within the fine-grained sediment, unlike massive gas hydrate deposits in the Krishna-Godavari basin. 2D multi-channel seismic section, which crosses the Holes NGHP-01-9A and 19B located at about 24 km apart shows a continuous bottom-simulating reflector (BSR) along it. We aim to investigate the prospect of gas hydrate accumulation in this area by integrating well log analysis and seismic methods with rock physics modeling. First, we estimate gas hydrate saturation at these two Holes from the observed impedance using the three-phase Biot-type equation (TPBE). Then we establish a linear relationship between gas hydrate saturation and impedance contrast with respect to the water-saturated sediment. Using this established relation and impedance obtained from pre-stack inversion of seismic data, we produce a 2D gas hydrate-distribution image over the entire seismic section. Gas hydrate saturation estimated from resistivity and sonic data at well locations varies within 0-15%, which agrees well with the available pressure core measurements at Hole 19. However, the 2D map of gas hydrate distribution obtained from our method shows maximum gas hydrate saturation is about 40% just above the BSR between the CDP (common depth point) 1450 and 2850. The presence of gas-charged sediments below the BSR is one of the reasons for the strong BSR observed in the seismic section, which is depicted as low impedance in the inverted impedance section. Closed sedimentary structures above the BSR are probably obstructing the movements of free-gas upslope, for which we do not see the presence of gas hydrate throughout the seismic section above the BSR.

Reservoir Prediction Under Control of Sedimentary Facies

2017 ◽  
Vol 25 (03) ◽  
pp. 1750022
Author(s):  
Xiuwei Yang ◽  
Peimin Zhu
Keyword(s):  
Seismic Data ◽  
Low Frequency ◽  
Prediction Method ◽  
Rock Physics ◽  
Real Data ◽  
Sedimentary Facies ◽  
Seismic Inversion ◽  
Rock Properties ◽  
Reservoir Prediction ◽  
Frequency Model

Acoustic impedance (AI) from seismic inversion can indicate rock properties and can be used, when combined with rock physics, to predict reservoir parameters, such as porosity. Solutions to seismic inversion problem are almost nonunique due to the limited bandwidth of seismic data. Additional constraints from well log data and geology are needed to arrive at a reasonable solution. In this paper, sedimentary facies is used to reduce the uncertainty in inversion and rock physics modeling; the results not only agree with seismic data, but also conform to geology. A reservoir prediction method, which incorporates seismic data, well logs, rock physics and sedimentary facies, is proposed. AI was first derived by constrained sparse spike inversion (CSSI) using a sedimentary facies dependent low-frequency model, and then was transformed to reservoir parameters by sequential simulation, statistical rock physics and [Formula: see text]-model. Two numerical experiments using synthetic model and real data indicated that the sedimentary facies information may help to obtain a more reasonable prediction.

Processing, inversion, and interpretation of a 2D seismic data set from the North Viking Graben, North Sea

Geophysics ◽  
2003 ◽  
Vol 68 (3) ◽  
pp. 837-848 ◽  
Author(s):  
Gislain B. Madiba ◽  
George A. McMechan
Keyword(s):  
Seismic Data ◽  
Dynamic Range ◽  
Velocity Ratio ◽  
Total Porosity ◽  
Shear Velocity ◽  
Seismic Section ◽  
Data Set ◽  
Amplitude Variation With Offset ◽  
Fluid Content ◽  
The North

Simultaneous elastic impedance inversion is performed on the 2D North Viking Graben seismic data set used at the 1994 SEG workshop on amplitude variation with offset and inversion. P‐velocity (Vp), S‐velocity (Vs), density logs, and seismic data are input to the inversion. The inverted P‐impedance and S‐impedance sections are used to generate an approximate compressional‐to‐shear velocity ratio (Vp/Vs) section which, in turn, is used along with water‐filled porosity (Swv) derived from the logs from two wells, to generate fluid estimate sections. This is possible as the reservoir sands have fairly constant total porosity of approximately 28 ± 4%, so the hydrocarbon filled porosity is the total porosity minus the water‐filled porosity. To enhance the separation of lithologies and of fluid content, we map Vp/Vs into Swv using an empirical crossplot‐derived relation. This mapping expands the dynamic range of the low end of the Vp/Vs values. The different lithologies and fluids are generally well separated in the Vp/Vs–Swv domain. Potential hydrocarbon reservoirs (as calibrated by the well data) are identified throughout the seismic section and are consistent with the fluid content estimations obtained from alternative computations. The Vp/Vs–Swv plane still does not produce unique interpretation in many situations. However, the critical distinction, which is between hydrocarbon‐bearing sands and all other geologic/reservoir configurations, is defined. Swv ≤ 0.17 and Vp/Vs ≤ 1.8 are the criteria that delineate potential reservoirs in this area, with decreasing Swv indicating a higher gas/oil ratio, and decreasing Vp/Vs indicating a higher sand/shale ratio. As these criteria are locally calibrated, they appear to be valid locally; they should not be applied to other data sets, which may exhibit significantly different relationships. However, the overall procedure should be generally applicable.

Probabilistic petrophysical-properties estimation integrating statistical rock physics with seismic inversion

Geophysics ◽  
2010 ◽  
Vol 75 (3) ◽  
pp. O21-O37 ◽  
Author(s):  
Dario Grana ◽  
Ernesto Della Rossa
Keyword(s):  
Seismic Data ◽  
Measurement Errors ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Joint Estimation ◽  
Physical Models ◽  
Inversion Method ◽  
Petrophysical Properties ◽  
Reservoir Properties ◽  
Statistical Rock Physics

A joint estimation of petrophysical properties is proposed that combines statistical rock physics and Bayesian seismic inversion. Because elastic attributes are correlated with petrophysical variables (effective porosity, clay content, and water saturation) and this physical link is associated with uncertainties, the petrophysical-properties estimation from seismic data can be seen as a Bayesian inversion problem. The purpose of this work was to develop a strategy for estimating the probability distributions of petrophysical parameters and litho-fluid classes from seismics. Estimation of reservoir properties and the associated uncertainty was performed in three steps: linearized seismic inversion to estimate the probabilities of elastic parameters, probabilistic upscaling to include the scale-changes effect, and petrophysical inversion to estimate the probabilities of petrophysical variables andlitho-fluid classes. Rock-physics equations provide the linkbetween reservoir properties and velocities, and linearized seismic modeling connects velocities and density to seismic amplitude. A full Bayesian approach was adopted to propagate uncertainty from seismics to petrophysics in an integrated framework that takes into account different sources of uncertainty: heterogeneity of the real data, approximation of physical models, measurement errors, and scale changes. The method has been tested, as a feasibility step, on real well data and synthetic seismic data to show reliable propagation of the uncertainty through the three different steps and to compare two statistical approaches: parametric and nonparametric. Application to a real reservoir study (including data from two wells and partially stacked seismic volumes) has provided as a main result the probability densities of petrophysical properties and litho-fluid classes. It demonstrated the applicability of the proposed inversion method.

scholarly journals Seismic Interpretation and Inversion Leading to an Accurate Reservoir Characterization in a Carbonate Environment

2021 ◽  
Vol 2 (12) ◽  
pp. 1229-1230
Author(s):  
Yasir Bashir ◽  
Nordiana Mohd Muztaza ◽  
Nur Azwin Ismail ◽  
Ismail Ahmad Abir ◽  
Andy Anderson Bery ◽  
...  
Keyword(s):  
Seismic Data ◽  
Acoustic Impedance ◽  
Reservoir Characterization ◽  
Rock Physics ◽  
Seismic Inversion ◽  
Seismic Interpretation ◽  
And Inversion

Seismic data acquired in the field show the subsurface reflectors or horizon among the geological strata, while the seismic inversion converts this reflector information into the acoustic impedance section which shows the layer properties based on lithology. The research aims to predict the porosity to identify the reservoir which is in between the tight layer. So, the output of the seismic inversion is much more batter than the seismic as it is closer to reality such as geology. Seismic inversion is frequently used to determine rock physics properties, for example, acoustic impedance and porosity.

Seismic reservoir characterization of the Bone Spring and Wolfcamp Formations in the Delaware Basin: Challenges and uncertainty in characterization using rock physics — A case study: Part 2

2020 ◽  
Vol 8 (4) ◽  
pp. T1057-T1069
Author(s):  
Ritesh Kumar Sharma ◽  
Satinder Chopra ◽  
Larry Lines
Keyword(s):  
Seismic Data ◽  
Reservoir Characterization ◽  
Rock Physics ◽  
Well Log ◽  
Data Set ◽  
Frequency Model ◽  
Log Data ◽  
Delaware Basin ◽  
Seismic Reservoir

The discrimination of fluid content and lithology in a reservoir is important because it has a bearing on reservoir development and its management. Among other things, rock-physics analysis is usually carried out to distinguish between the lithology and fluid components of a reservoir by way of estimating the volume of clay, water saturation, and porosity using seismic data. Although these rock-physics parameters are easy to compute for conventional plays, there are many uncertainties in their estimation for unconventional plays, especially where multiple zones need to be characterized simultaneously. We have evaluated such uncertainties with reference to a data set from the Delaware Basin where the Bone Spring, Wolfcamp, Barnett, and Mississippian Formations are the prospective zones. Attempts at seismic reservoir characterization of these formations have been developed in Part 1 of this paper, where the geologic background of the area of study, the preconditioning of prestack seismic data, well-log correlation, accounting for the temporal and lateral variation in the seismic wavelets, and building of robust low-frequency model for prestack simultaneous impedance inversion were determined. We determine the challenges and the uncertainty in the characterization of the Bone Spring, Wolfcamp, Barnett, and Mississippian sections and explain how we overcame those. In the light of these uncertainties, we decide that any deterministic approach for characterization of the target formations of interest may not be appropriate and we build a case for adopting a robust statistical approach. Making use of neutron porosity and density porosity well-log data in the formations of interest, we determine how the type of shale, volume of shale, effective porosity, and lithoclassification can be carried out. Using the available log data, multimineral analysis was also carried out using a nonlinear optimization approach, which lent support to our facies classification. We then extend this exercise to derived seismic attributes for determination of the lithofacies volumes and their probabilities, together with their correlations with the facies information derived from mud log data.

Sign in / Sign up

Export Citation Format

Share Document