scholarly journals Brittle mineral prediction based on rock-physics modelling for tight oil reservoir rocks

2021 ◽  
Vol 18 (6) ◽  
pp. 970-983
Author(s):  
Jing Ba ◽  
Peng Hu ◽  
Wenhui Tan ◽  
Tobias M Müller ◽  
Li-Yun Fu

Abstract The reservoir rocks from Chang-7 member of Yanchang Formation of Ordos Basin are characterised with heterogeneous fabric structures at the pore scale, and low porosity/permeability is exhibited at the macro scale. Precise prediction of reservoir brittleness is of great significance to oil production. Ultrasonic experiments are performed on tight sandstones collected from the target formation. A rock-physics model (RPM) is presented based on the Voigt–Reuss–Hill average (VRH), self-consistent approximation (SCA) and differential effective medium (DEM) theory. The brittleness characteristics relying on mineral composition, porosity and microcrack properties are explored by using the RPM. The Young's modulus increases and Poisson ratio decreases with increasing quartz content. Based on experimental, log and seismic data, brittle mineral analysis of rock physical model is performed at multiple scales. The model accuracy is verified by experimental data and well log data. The brittleness distribution is predicted on the basis of log and seismic data, which can be instructive for the reservoir rock fracturing in actual engineering operations.

2021 ◽  
Vol 40 (12) ◽  
pp. 897-904
Author(s):  
Manuel González-Quijano ◽  
Gregor Baechle ◽  
Miguel Yanez ◽  
Freddy Obregon ◽  
Carmen Vito ◽  
...  

The study area is located in middepth to deep waters of the Salina del Istmo Basin where Repsol operates Block 29. The objective of this work is to integrate qualitative and quantitative interpretations of rock and seismic data to predict lithology and fluid of the Early Miocene prospects. The seismic expression of those prospects differs from age-equivalent well-studied analog fields in the U.S. Gulf of Mexico Basin due to the mineralogically complex composition of abundant extrusive volcanic material. Offset well data (i.e., core, logs, and cuttings) were used to discriminate lithology types and to quantify mineralogy. This analysis served as input for developing a new rock-physics framework and performing amplitude variation with offset (AVO) modeling. The results indicate that the combination of intercept and gradient makes it possible to discriminate hydrocarbon-filled (AVO class II and III) versus nonhydrocarbon-filled rocks (AVO class 0 and IV). Different lithologies within hydrocarbon-bearing reservoirs cannot be discriminated as the gradient remains negative for all rock types. However, AVO analysis allows discrimination of three different reservoir rock types in water-bearing cases (AVO class 0, I, and IV). These conclusions were obtained during studies conducted in 2018–2019 and were used in prospect evaluation to select drilling locations leading to two wildcat discoveries, the Polok and Chinwol prospects, drilled in Block 29 in 2020.


Geophysics ◽  
2005 ◽  
Vol 70 (3) ◽  
pp. O1-O11 ◽  
Author(s):  
Alexey Stovas ◽  
Martin Landrø

We investigate how seismic anisotropy influences our ability to distinguish between various production-related effects from time-lapse seismic data. Based on rock physics models and ultrasonic core measurements, we estimate variations in PP and PS reflectivity at the top reservoir interface for fluid saturation and pore pressure changes. The tested scenarios include isotropic shale, weak anisotropic shale, and highly anisotropic shale layers overlaying either an isotropic reservoir sand layer or a weak anisotropic sand layer. We find that, for transverse isotropic media with a vertical symmetry axis (TIV), the effect of weak anisotropy in the cap rock does not lead to significant errors in, for instance, the simultaneous determination of pore-pressure and fluid-saturation changes. On the other hand, changes in seismic anisotropy within the reservoir rock (caused by, for instance, increased fracturing) might be detectable from time-lapse seismic data. A new method using exact expressions for PP and PS reflectivity, including TIV anisotropy, is used to determine pressure and saturation changes over production time. This method is assumed to be more accurate than previous methods.


2005 ◽  
Vol 7 ◽  
pp. 13-16
Author(s):  
Peter Japsen ◽  
Anders Bruun ◽  
Ida L. Fabricius ◽  
Gary Mavko

Seismic data are mainly used to map out structures in the subsurface, but are also increasingly used to detect differences in porosity and in the fluids that occupy the pore space in sedimentary rocks. Hydrocarbons are generally lighter than brine, and the bulk density and sonic velocity (speed of pressure waves or P-wave velocity) of hydrocarbon-bearing sedimentary rocks are therefore reduced compared to non-reservoir rocks. However, sound is transmitted in different wave forms through the rock, and the shear velocity (speed of shear waves or S-wave velocity) is hardly affected by the density of the pore fluid. In order to detect the presence of hydrocarbons from seismic data, it is thus necessary to investigate how porosity and pore fluids affect the acoustic properties of a sedimentary rock. Much previous research has focused on describing such effects in sandstone (see Mavko et al. 1998), and only in recent years have corresponding studies on the rock physics of chalk appeared (e.g. Walls et al. 1998; Røgen 2002; Fabricius 2003; Gommesen 2003; Japsen et al. 2004). In the North Sea, chalk of the Danian Ekofisk Formation and the Maastrichtian Tor Formation are important reservoir rocks. More information could no doubt be extracted from seismic data if the fundamental physical properties of chalk were better understood. The presence of gas in chalk is known to cause a phase reversal in the seismic signal (Megson 1992), but the presence of oil in chalk has only recently been demonstrated to have an effect on surface seismic data (Japsen et al. 2004). The need for a better link between chalk reservoir parameters and geophysical observations has, however, strongly increased since the discovery of the Halfdan field proved major reserves outside four-way dip closures (Jacobsen et al. 1999; Vejbæk & Kristensen 2000).


2015 ◽  
Vol 3 (2) ◽  
pp. SM23-SM35
Author(s):  
Russell W. Carter ◽  
Kyle T. Spikes

Large-scale subsurface injection of [Formula: see text] has the potential to reduce emissions of atmospheric [Formula: see text] and improve oil recovery. Studying the effects of injected [Formula: see text] on the elastic properties of the saturated reservoir rock can help to improve long-term monitoring effectiveness and accuracy at locations undergoing [Formula: see text] injection. We used two vintages of existing 3D surface seismic data and well logs to probabilistically invert for the [Formula: see text] saturation and porosity at the Cranfield reservoir using a double-difference approach. The first step of this work was to calibrate the rock-physics model to the well-log data. Next, the baseline and time-lapse seismic data sets were inverted for acoustic impedance [Formula: see text] using a high-resolution basis pursuit inversion technique. The reservoir porosity was derived statistically from the rock-physics model based on the [Formula: see text] estimates inverted from the baseline survey. The porosity estimates were used in the double-difference routine as the fixed initial model from which [Formula: see text] saturation was then estimated from the time-lapse [Formula: see text] data. Porosity was assumed to remain constant between survey vintages; therefore, the changes between the baseline and time-lapse [Formula: see text] data may be inverted for [Formula: see text] saturation from the injection activities using the calibrated rock-physics model. Comparisons of inverted and measured porosity from well logs indicated quite accurate results. Estimates of [Formula: see text] saturation found less accuracy than the porosity estimates.


Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. C177-C191 ◽  
Author(s):  
Yunyue Li ◽  
Biondo Biondi ◽  
Robert Clapp ◽  
Dave Nichols

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.


2021 ◽  
Author(s):  
S Al Naqbi ◽  
J Ahmed ◽  
J Vargas Rios ◽  
Y Utami ◽  
A Elila ◽  
...  

Abstract The Thamama group of reservoirs consist of porous carbonates laminated with tight carbonates, with pronounced lateral heterogeneities in porosity, permeability, and reservoir thickness. The main objective of our study was mapping variations and reservoir quality prediction away from well control. As the reservoirs were thin and beyond seismic resolution, it was vital that the facies and porosity be mapped in high resolution, with a high predictability, for successful placement of horizontal wells for future development of the field. We established a unified workflow of geostatistical inversion and rock physics to characterize the reservoirs. Geostatistical inversion was run in static models that were converted from depth to time domain. A robust two-way velocity model was built to map the depth grid and its zones on the time seismic data. This ensured correct placement of the predicted high-resolution elastic attributes in the depth static model. Rock physics modeling and Bayesian classification were used to convert the elastic properties into porosity and lithology (static rock-type (SRT)), which were validated in blind wells and used to rank the multiple realizations. In the geostatistical pre-stack inversion, the elastic property prediction was constrained by the seismic data and controlled by variograms, probability distributions and a guide model. The deterministic inversion was used as a guide or prior model and served as a laterally varying mean. Initially, unconstrained inversion was tested by keeping all wells as blind and the predictions were optimized by updating the input parameters. The stochastic inversion results were also frequency filtered in several frequency bands, to understand the impact of seismic data and variograms on the prediction. Finally, 30 wells were used as input, to generate 80 realizations of P-impedance, S-impedance, Vp/Vs, and density. After converting back to depth, 30 additional blind wells were used to validate the predicted porosity, with a high correlation of more than 0.8. The realizations were ranked based on the porosity predictability in blind wells combined with the pore volume histograms. Realizations with high predictability and close to the P10, P50 and P90 cases (of pore volume) were selected for further use. Based on the rock physics analysis, the predicted lithology classes were associated with the geological rock-types (SRT) for incorporation in the static model. The study presents an innovative approach to successfully integrate geostatistical inversion and rock physics with static modeling. This workflow will generate seismically constrained high-resolution reservoir properties for thin reservoirs, such as porosity and lithology, which are seamlessly mapped in the depth domain for optimized development of the field. It will also account for the uncertainties in the reservoir model through the generation of multiple equiprobable realizations or scenarios.


2021 ◽  
pp. 1-59
Author(s):  
Kai Lin ◽  
Xilei He ◽  
Bo Zhang ◽  
Xiaotao Wen ◽  
Zhenhua He ◽  
...  

Most of current 3D reservoir’s porosity estimation methods are based on analyzing the elastic parameters inverted from seismic data. It is well-known that elastic parameters vary with pore structure parameters such as pore aspect ratio, consolidate coefficient, critical porosity, etc. Thus, we may obtain inaccurate 3D porosity estimation if the chosen rock physics model fails properly address the effects of pore structure parameters on the elastic parameters. However, most of current rock physics models only consider one pore structure parameter such as pore aspect ratio or consolidation coefficient. To consider the effect of multiple pore structure parameters on the elastic parameters, we propose a comprehensive pore structure (CPS) parameter set that is generalized from the current popular rock physics models. The new CPS set is based on the first order approximation of current rock physics models that consider the effect of pore aspect ratio on elastic parameters. The new CPS set can accurately simulate the behavior of current rock physics models that consider the effect of pore structure parameters on elastic parameters. To demonstrate the effectiveness of proposed parameters in porosity estimation, we use a theoretical model to demonstrate that the proposed CPS parameter set properly addresses the effect of pore aspect ratio on elastic parameters such as velocity and porosity. Then, we obtain a 3D porosity estimation for a tight sand reservoir by applying it seismic data. We also predict the porosity of the tight sand reservoir by using neural network algorithm and a rock physics model that is commonly used in porosity estimation. The comparison demonstrates that predicted porosity has higher correlation with the porosity logs at the blind well locations.


Author(s):  
A. Ogbamikhumi ◽  
T. Tralagba ◽  
E. E. Osagiede

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.


Sign in / Sign up

Export Citation Format

Share Document