Quantitative characterization of hydrocarbon reservoir using integrated seismic rock physics analysis: an integrated approach using seismic data, seismic rock physics of well-log and core

2009 ◽  
Vol 2009 (1) ◽  
pp. 1
Author(s):  
B. Nurhandoko ◽  
M. Choliq ◽  
K. Triyoso ◽  
I. Soemantri ◽  
S. Praptono ◽  
...  
2020 ◽  
Vol 8 (4) ◽  
pp. T1057-T1069
Author(s):  
Ritesh Kumar Sharma ◽  
Satinder Chopra ◽  
Larry Lines

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.


2018 ◽  
Vol 6 (2) ◽  
pp. T325-T336 ◽  
Author(s):  
Ritesh Kumar Sharma ◽  
Satinder Chopra ◽  
James Keay ◽  
Hossein Nemati ◽  
Larry Lines

The Utica Formation in eastern Ohio possesses all the prerequisites for being a successful unconventional play. Attempts at seismic reservoir characterization of the Utica Formation have been discussed in part 1, in which, after providing the geologic background of the area of study, the preconditioning of prestack seismic data, well-log correlation, and building of robust low-frequency models for prestack simultaneous impedance inversion were explained. All these efforts were aimed at identification of sweet spots in the Utica Formation in terms of organic richness as well as brittleness. We elaborate on some aspects of that exercise, such as the challenges we faced in the determination of the total organic carbon (TOC) volume and computation of brittleness indices based on mineralogical and geomechanical considerations. The prediction of TOC in the Utica play using a methodology, in which limited seismic as well as well-log data are available, is demonstrated first. Thereafter, knowing the nonexistence of the universally accepted indicator of brittleness, mechanical along with mineralogical attempts to extract the brittleness information for the Utica play are discussed. Although an attempt is made to determine brittleness from mechanical rock-physics parameters (Young’s modulus and Poisson’s ratio) derived from seismic data, the available X-ray diffraction data and regional petrophysical modeling make it possible to determine the brittleness index based on mineralogical data and thereafter be derived from seismic data.


Geophysics ◽  
2001 ◽  
Vol 66 (4) ◽  
pp. 988-1001 ◽  
Author(s):  
T. Mukerji ◽  
A. Jørstad ◽  
P. Avseth ◽  
G. Mavko ◽  
J. R. Granli

Reliably predicting lithologic and saturation heterogeneities is one of the key problems in reservoir characterization. In this study, we show how statistical rock physics techniques combined with seismic information can be used to classify reservoir lithologies and pore fluids. One of the innovations was to use a seismic impedance attribute (related to the [Formula: see text] ratio) that incorporates far‐offset data, but at the same time can be practically obtained using normal incidence inversion algorithms. The methods were applied to a North Sea turbidite system. We incorporated well log measurements with calibration from core data to estimate the near‐offset and far‐offset reflectivity and impedance attributes. Multivariate probability distributions were estimated from the data to identify the attribute clusters and their separability for different facies and fluid saturations. A training data was set up using Monte Carlo simulations based on the well log—derived probability distributions. Fluid substitution by Gassmann’s equation was used to extend the training data, thus accounting for pore fluid conditions not encountered in the well. Seismic inversion of near‐offset and far‐offset stacks gave us two 3‐D cubes of impedance attributes in the interwell region. The near‐offset stack approximates a zero‐offset section, giving an estimate of the normal incidence acoustic impedance. The far offset stack gives an estimate of a [Formula: see text]‐related elastic impedance attribute that is equivalent to the acoustic impedance for non‐normal incidence. These impedance attributes obtained from seismic inversion were then used with the training probability distribution functions to predict the probability of occurrence of the different lithofacies in the interwell region. Statistical classification techniques, as well as geostatistical indicator simulations were applied on the 3‐D seismic data cube. A Markov‐Bayes technique was used to update the probabilities obtained from the seismic data by taking into account the spatial correlation as estimated from the facies indicator variograms. The final results are spatial 3‐D maps of not only the most likely facies and pore fluids, but also their occurrence probabilities. A key ingredient in this study was the exploitation of physically based seismic‐to‐reservoir property transforms optimally combined with statistical techniques.


Solid Earth ◽  
2016 ◽  
Vol 7 (3) ◽  
pp. 943-958 ◽  
Author(s):  
Xènia Ogaya ◽  
Juan Alcalde ◽  
Ignacio Marzán ◽  
Juanjo Ledo ◽  
Pilar Queralt ◽  
...  

Abstract. Hontomín (N of Spain) hosts the first Spanish CO2 storage pilot plant. The subsurface characterization of the site included the acquisition of a 3-D seismic reflection and a circumscribed 3-D magnetotelluric (MT) survey. This paper addresses the combination of the seismic and MT results, together with the available well-log data, in order to achieve a better characterization of the Hontomín subsurface. We compare the structural model obtained from the interpretation of the seismic data with the geoelectrical model resulting from the MT data. The models correlate well in the surroundings of the CO2 injection area with the major structural differences observed related to the presence of faults. The combination of the two methods allowed a more detailed characterization of the faults, defining their geometry, and fluid flow characteristics, which are key for the risk assessment of the storage site. Moreover, we use the well-log data of the existing wells to derive resistivity–velocity relationships for the subsurface and compute a 3-D velocity model of the site using the 3-D resistivity model as a reference. The derived velocity model is compared to both the predicted and logged velocity in the injection and monitoring wells, for an overall assessment of the computed resistivity–velocity relationships. The major differences observed are explained by the different resolution of the compared geophysical methods. Finally, the derived velocity model for the near surface is compared with the velocity model used for the static corrections in the seismic data. The results allowed extracting information about the characteristics of the shallow unconsolidated sediments, suggesting possible clay and water content variations. The good correlation of the velocity models derived from the resistivity–velocity relationships and the well-log data demonstrate the potential of the combination of the two methods for characterizing the subsurface, in terms of its physical properties (velocity, resistivity) and structural/reservoir characteristics. This work explores the compatibility of the seismic and magnetotelluric methods across scales highlighting the importance of joint interpretation in near surface and reservoir characterization.


Sign in / Sign up

Export Citation Format

Share Document