Seismic reservoir characterization of Utica-Point Pleasant Shale with efforts at quantitative interpretation — A case study: Part 1

2018 ◽  
Vol 6 (2) ◽  
pp. T313-T324 ◽  
Author(s):  
Satinder Chopra ◽  
Ritesh Kumar Sharma ◽  
Hossein Nemati ◽  
James Keay

The Utica Shale is one of the major source rocks in Ohio, and it extends across much of the eastern United States. Its organic richness, high content of calcite, and development of extensive organic porosity make it a perfect unconventional play, and it has gained the attention of the oil and gas industry. The primary target zone in the Utica Play includes the Utica Formation, Point Pleasant Formation, and Trenton Formation intervals. We attempt to identify the sweet spots within the Point Pleasant interval using 3D seismic data, available well data, and other relevant data. This has been done by way of organic richness and brittleness estimation in the rock intervals. The organic richness is determined by weight % of total organic carbon content, which is derived by transforming the inverted density volume. Core-log petrophysical modeling provides the necessary relationship for doing so. The brittleness is derived using rock-physics parameters such as the Young’s modulus and Poisson’s ratio. Deterministic simultaneous inversion along with a neural network approach are followed to compute the rock-physics parameters and density using seismic data. The correlation of sweet spots identified based on the seismic data with the available production data emphasizes the significance of integration of seismic data with all other relevant data.

Geophysics ◽  
2021 ◽  
Vol 86 (6) ◽  
pp. M197-M209
Author(s):  
Kun Luo ◽  
Zhaoyun Zong ◽  
Xingyao Yin ◽  
Hong Cao ◽  
Minghui Lu

A Gaussian mixture Hamiltonian Monte Carlo (HMC) Bayesian method has been developed for the inversion of petrophysical parameters such as pyrolysis parameter S1, which is driven by a statistical shale rock-physics model. Pyrolysis parameter S1 can be used to indicate the content of free or adsorbed hydrocarbons in source rock, and it is an important indicator to evaluate the production of shale oil reservoirs. However, most studies on pyrolysis parameters are based on pyrolysis experiments and there is no relevant study to inverse pyrolysis parameter S1 from seismic data. In addition, compared to the total organic carbon content, pyrolysis S1 is more accurate for evaluating gas and oil in shale. In particular, high values of pyrolysis S1 can directly indicate the content of shale oil. We have developed a strategy for assessing shale oil sweet spots through estimating pyrolysis S1 and other petrophysical parameters. Based on the Gaussian mixture assumptions for the prior distribution of the model, we build a joint distribution to link the pyrolysis parameter S1 with elastic attributes, and then we derive a formulation to inverse S1 with the Bayesian model. Due to the components of the Gaussian mixture, the HMC method has been used to sample the posterior distribution. Our study finds that the HMC method for sampling can improve the efficiency and allow a more robust quantification of the uncertainty; also, application to real seismic data sets indicates that the delineation of sweet spots is more accurate combined with pyrolysis S1.


2016 ◽  
Vol 35 (2) ◽  
pp. 147-171 ◽  
Author(s):  
Sheng Chen ◽  
Wenzhi Zhao ◽  
Yonglin Ouyang ◽  
Qingcai Zeng ◽  
Qing Yang ◽  
...  

W4 block of Sichuan Basin is a pioneer in shale gas exploration and development in China. But geophysical prospecting is just at its beginning and thus has not provided enough information about how sweet spots distribute for the deployment of horizontal well. This paper predicted sweet spots based on logging and 3D seismic data. Well logging interpretation method was used to get the key evaluation parameters of shale reservoir and determine the distribution of sweet spots in vertical direction. Rock physics analysis technology was used to define the elastic parameters that were sensitive to the key evaluation parameters, such as TOC and gas content of shale gas reservoir. At the same time the quantitative relationships between them were established. Based on the result of seismic rock physics analysis, prestack inversion was carried out to predict the transverse plane distribution of the key evaluation parameters of shale reservoir. These research results are integrated to determine the distribution of sweet spots. The results show that sweet spots in this area were characterized by high TOC content, high gas content, high GR, high Young’s modulus, low Poisson’s ratio, low density, and low P-wave velocity. Density was the most sensitive elastic parameters to TOC of the reservoir. The optimal combination for predicting the gas content is composed of six parameters include density, Poisson’s ratio, and so on. Sweet spots in this block vertically concentrate within 30 m above the bottom of Longmaxi Formation. Two classes of sweet spots have been predicted in this area, class I sweet spots are recommended to be prioritized for development. This study effectively predicted the spatial distribution of sweet spots, which provide important guidance for the development of the area.


2017 ◽  
Vol 5 (4) ◽  
pp. T641-T652 ◽  
Author(s):  
Mark Sams ◽  
Paul Begg ◽  
Timur Manapov

The information within seismic data is band limited and angle limited. Together with the particular physics and geology of carbonate rocks, this imposes limitations on how accurately we can predict the presence of hydrocarbons in carbonates, map the top carbonate, and characterize the porosity distribution through seismic amplitude analysis. Using data for a carbonate reef from the Nam Con Son Basin, Vietnam, the expectations based on rock-physics analysis are that the presence of gas can be predicted only when the porosity at the top of the carbonate is extremely high ([Formula: see text]), but that a fluid contact is unlikely to be observed in the background of significant porosity variations. Mapping the top of the carbonate (except when the top carbonate porosities are low) or a fluid contact requires accurate estimates of changes in [Formula: see text]. The seismic data do not independently support such an accurate estimation of sharp changes in [Formula: see text]. The standard approach of introducing low-frequency models and applying rock-physics constraints during a simultaneous inversion does not resolve the problems: The results are heavily biased by the well control and the initial interpretation of the top carbonate and fluid contact. A facies-based inversion in which the elastic properties are restricted to values consistent with the facies predicted to be present removes the well bias, but it does not completely obviate the need for a reasonably accurate initial interpretation in terms of prior facies probability distributions. Prestack inversion improves the quality of the facies predictions compared with a poststack inversion.


Geophysics ◽  
2008 ◽  
Vol 73 (3) ◽  
pp. C13-C21 ◽  
Author(s):  
Arild Buland ◽  
Odd Kolbjørnsen ◽  
Ragnar Hauge ◽  
Øyvind Skjæveland ◽  
Kenneth Duffaut

A fast Bayesian inversion method for 3D lithology and fluid prediction from prestack seismic data, and a corresponding feasibility analysis were developed and tested on a real data set. The objective of the inversion is to find the probabilities for different lithology-fluid classes from seismic data and geologic knowledge. The method combines stochastic rock physics relations between the elastic parameters and the different lithology-fluid classes with the results from a fast Bayesian seismic simultaneous inversion from seismic data to elastic parameters. A method for feasibility analysis predicts the expected modification of the prior probabilities to posterior probabilities for the different lithology-fluid classes. The feasibility analysis can be carried out before the seismic data are analyzed. Both the feasibility method and the seismic lithology-fluid probability inversion were applied to a prospect offshore Norway. The analysis improves the probability for gas sand from 0.1 to about 0.2–0.4 with seismic data.


Author(s):  
Nina Skaarup ◽  
James A. Chalmers

NOTE: This article was published in a former series of GEUS Bulletin. Please use the original series name when citing this article, for example: Skaarup, N., & Chalmers, J. A. (1998). A possible new hydrocarbon play, offshore central West Greenland. Geology of Greenland Survey Bulletin, 180, 28-30. https://doi.org/10.34194/ggub.v180.5082 _______________ The discovery of extensive seeps of crude oil onshore central West Greenland (Christiansen et al. 1992, 1994, 1995, 1996, 1997, 1998, this volume; Christiansen 1993) means that the central West Greenland area is now prospective for hydrocarbons in its own right. Analysis of the oils (Bojesen-Koefoed et al. in press) shows that their source rocks are probably nearby and, because the oils are found within the Lower Tertiary basalts, the source rocks must be below the basalts. It is therefore possible that in the offshore area oil could have migrated through the basalts and be trapped in overlying sediments. In the offshore area to the west of Disko and Nuussuaq (Fig. 1), Whittaker (1995, 1996) interpreted a few multichannel seismic lines acquired in 1990, together with some seismic data acquired by industry in the 1970s. He described a number of large rotated fault-blocks containing structural closures at top basalt level that could indicate leads capable of trapping hydrocarbons. In order to investigate Whittaker’s (1995, 1996) interpretation, in 1995 the Geological Survey of Greenland acquired 1960 km new multichannel seismic data (Fig. 1) using funds provided by the Government of Greenland, Minerals Office (now Bureau of Minerals and Petroleum) and the Danish State through the Mineral Resources Administration for Greenland. The data were acquired using the Danish Naval vessel Thetis which had been adapted to accommodate seismic equipment. The data acquired in 1995 have been integrated with the older data and an interpretation has been carried out of the structure of the top basalt reflection. This work shows a fault pattern in general agreement with that of Whittaker (1995, 1996), although there are differences in detail. In particular the largest structural closure reported by Whittaker (1995) has not been confirmed. Furthermore, one of Whittaker’s (1995) smaller leads seems to be larger than he had interpreted and may be associated with a DHI (direct hydrocarbon indicator) in the form of a ‘bright spot’.


Geology ◽  
2011 ◽  
Vol 39 (12) ◽  
pp. 1167-1170 ◽  
Author(s):  
Helge Løseth ◽  
Lars Wensaas ◽  
Marita Gading ◽  
Kenneth Duffaut ◽  
Michael Springer

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):  
Rick Schrynemeeckers

Abstract Current offshore hydrocarbon detection methods employ vessels to collect cores along transects over structures defined by seismic imaging which are then analyzed by standard geochemical methods. Due to the cost of core collection, the sample density over these structures is often insufficient to map hydrocarbon accumulation boundaries. Traditional offshore geochemical methods cannot define reservoir sweet spots (i.e. areas of enhanced porosity, pressure, or net pay thickness) or measure light oil or gas condensate in the C7 – C15 carbon range. Thus, conventional geochemical methods are limited in their ability to help optimize offshore field development production. The capability to attach ultrasensitive geochemical modules to Ocean Bottom Seismic (OBS) nodes provides a new capability to the industry which allows these modules to be deployed in very dense grid patterns that provide extensive coverage both on structure and off structure. Thus, both high resolution seismic data and high-resolution hydrocarbon data can be captured simultaneously. Field trials were performed in offshore Ghana. The trial was not intended to duplicate normal field operations, but rather provide a pilot study to assess the viability of passive hydrocarbon modules to function properly in real world conditions in deep waters at elevated pressures. Water depth for the pilot survey ranged from 1500 – 1700 meters. Positive thermogenic signatures were detected in the Gabon samples. A baseline (i.e. non-thermogenic) signature was also detected. The results indicated the positive signatures were thermogenic and could easily be differentiated from baseline or non-thermogenic signatures. The ability to deploy geochemical modules with OBS nodes for reoccurring surveys in repetitive locations provides the ability to map the movement of hydrocarbons over time as well as discern depletion affects (i.e. time lapse geochemistry). The combined technologies will also be able to: Identify compartmentalization, maximize production and profitability by mapping reservoir sweet spots (i.e. areas of higher porosity, pressure, & hydrocarbon richness), rank prospects, reduce risk by identifying poor prospectivity areas, accurately map hydrocarbon charge in pre-salt sequences, augment seismic data in highly thrusted and faulted areas.


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.


Sign in / Sign up

Export Citation Format

Share Document