Quantitative prediction of total organic carbon content in shale-gas reservoirs using seismic data: A case study from the Lower Silurian Longmaxi Formation in the Chang Ning gas field of the Sichuan Basin, China

2018 ◽  
Vol 6 (4) ◽  
pp. SN153-SN168 ◽  
Author(s):  
Sheng Chen ◽  
Wenzhi Zhao ◽  
Qingcai Zeng ◽  
Qing Yang ◽  
Pei He ◽  
...  

We present a quantitative prediction of total organic carbon (TOC) content for shale-gas development in the Chang Ning gas field of the Sichuan Basin (China). We have used the rock-physics analysis method to define the geophysical characteristics of the reservoir and the most sensitive elastic parameter to TOC content. We established a quantitative prediction template of the TOC content by rock-physics modeling. Well data and 3D seismic data were combined for prestack simultaneous inversion to obtain the most sensitive elastic parameter data volume. According to the prediction template, we transformed the sensitive elastic parameter data volume to the TOC content volume. The rock-physics analysis indicates that the reservoir with a high TOC content in the Lower Silurian Longmaxi Formation (Fm) of the Chang Ning (CN) gas field is characterized by low density, low P-wave velocity ([Formula: see text]), low S-wave velocity ([Formula: see text]), low Poisson’s ratio (PR), and low ratio of P-wave velocity to S-wave velocity ([Formula: see text]). Density is the most sensitive elastic parameter to TOC content. The rock-physics model suggests that density is negatively correlated with TOC content, and the relationship between them changes under different porosities. The reservoir with high TOC content is mainly distributed at the bottom of the Longmaxi Fm and in the central and east central area of the study field. The quantitative prediction results are in good agreement with the log interpretation and production test. Therefore, it has important implications for the efficient development of the shale-gas reservoir in the basin.

2013 ◽  
Vol 788 ◽  
pp. 701-704
Author(s):  
Xin Gong Tang ◽  
Jian Bin Ma ◽  
Kui Xiang ◽  
Liang Jun Yan ◽  
Wen Bao Hu

Petrophysical study is playing an important role in oil and gas exploration. Shale gas and shale oil is blooming in recent years in many countries. Less rock physics knowledge is known about shale relatively to other rock type such as sandstone and limestone. In this paper, we carried out a rock physical study of shale core sample which is drilled from north China. The plan distribution of permeability, P wave velocity, S wave velocity and complex resistivity were acquired based on AutoScan-IIplatform. The results show that the permeability of the shale sample is basically low with values of 0.1 to several micro Darcy (mD) except some fracture areas in the surface, which has values of about several tens mD. The permeability can basically describe the distribution of the fracture. The complex resistivity has the similar characteristics with permeability, which is also roughly corresponding to the position of the facture. As for the Vp and Vs, although not very good correspondence with the surface, they are still approximately present the high and low velocity feature of the core sample as well. This result is significantly helpful for shale gas exploration and production.


2021 ◽  
Author(s):  
Sheng Chen ◽  
Qingcai Zeng ◽  
Xiujiao Wang ◽  
Qing Yang ◽  
Chunmeng Dai ◽  
...  

Abstract Practices of marine shale gas exploration and development in south China have proved that formation overpressure is the main controlling factor of shale gas enrichment and an indicator of good preservation condition. Accurate prediction of formation pressure before drilling is necessary for drilling safety and important for sweet spots predicting and horizontal wells deploying. However, the existing prediction methods of formation pore pressures all have defects, the prediction accuracy unsatisfactory for shale gas development. By means of rock mechanics analysis and related formulas, we derived a formula for calculating formation pore pressures. Through regional rock physical analysis, we determined and optimized the relevant parameters in the formula, and established a new formation pressure prediction model considering P-wave velocity, S-wave velocity and density. Based on regional exploration wells and 3D seismic data, we carried out pre-stack seismic inversion to obtain high-precision P-wave velocity, S-wave velocity and density data volumes. We utilized the new formation pressure prediction model to predict the pressure and the spatial distribution of overpressure sweet spots. Then, we applied the measured pressure data of three new wells to verify the predicted formation pressure by seismic data. The result shows that the new method has a higher accuracy. This method is qualified for safe drilling and prediction of overpressure sweet spots for shale gas development, so it is worthy of promotion.


2019 ◽  
Vol 38 (10) ◽  
pp. 762-769
Author(s):  
Patrick Connolly

Reflectivities of elastic properties can be expressed as a sum of the reflectivities of P-wave velocity, S-wave velocity, and density, as can the amplitude-variation-with-offset (AVO) parameters, intercept, gradient, and curvature. This common format allows elastic property reflectivities to be expressed as a sum of AVO parameters. Most AVO studies are conducted using a two-term approximation, so it is helpful to reduce the three-term expressions for elastic reflectivities to two by assuming a relationship between P-wave velocity and density. Reduced to two AVO components, elastic property reflectivities can be represented as vectors on intercept-gradient crossplots. Normalizing the lengths of the vectors allows them to serve as basis vectors such that the position of any point in intercept-gradient space can be inferred directly from changes in elastic properties. This provides a direct link between properties commonly used in rock physics and attributes that can be measured from seismic data. The theory is best exploited by constructing new seismic data sets from combinations of intercept and gradient data at various projection angles. Elastic property reflectivity theory can be transferred to the impedance domain to aid in the analysis of well data to help inform the choice of projection angles. Because of the effects of gradient measurement errors, seismic projection angles are unlikely to be the same as theoretical angles or angles derived from well-log analysis, so seismic data will need to be scanned through a range of angles to find the optimum.


Geophysics ◽  
1998 ◽  
Vol 63 (5) ◽  
pp. 1659-1669 ◽  
Author(s):  
Christine Ecker ◽  
Jack Dvorkin ◽  
Amos Nur

We interpret amplitude variation with offset (AVO) data from a bottom simulating reflector (BSR) offshore Florida by using rock‐physics‐based synthetic seismic models. A previously conducted velocity and AVO analysis of the in‐situ seismic data showed that the BSR separates hydrate‐bearing sediments from sediments containing free methane. The amplitude at the BSR are increasingly negative with increasing offset. This behavior was explained by P-wave velocity above the BSR being larger than that below the BSR, and S-wave velocity above the BSR being smaller than that below the BSR. We use these AVO and velocity results to infer the internal structure of the hydrated sediment. To do so, we examine two micromechanical models that correspond to the two extreme cases of hydrate deposition in the pore space: (1) the hydrate cements grain contacts and strongly reinforces the sediment, and (2) the hydrate is located away from grain contacts and does not affect the stiffness of the sediment frame. Only the second model can qualitatively reproduce the observed AVO response. Thus inferred internal structure of the hydrate‐bearing sediment means that (1) the sediment above the BSR is uncemented and, thereby, mechanically weak, and (2) its permeability is very low because the hydrate clogs large pore‐space conduits. The latter explains why free gas is trapped underneath the BSR. The seismic data also indicate the absence of strong reflections at the top of the hydrate layer. This fact suggests that the high concentration of hydrates in the sediment just above the BSR gradually decreases with decreasing depth. This effect is consistent with the fact that the low‐permeability hydrated sediments above the BSR prevent free methane from migrating upwards.


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).


Geophysics ◽  
1994 ◽  
Vol 59 (9) ◽  
pp. 1362-1376 ◽  
Author(s):  
Jan L. Fatti ◽  
George C. Smith ◽  
Peter J. Vail ◽  
Peter J. Strauss ◽  
Philip R. Levitt

The Geostack technique is a method of analyzing seismic amplitude variation with offset (AVO) information. One of the outputs of the analysis is a set of direct hydrocarbon indicator traces called “fluid factor” traces. The fluid factor trace is designed to be low amplitude for all reflectors in a clastic sedimentary sequence except for rocks that lie off the (mudrock line.) The mudrock line is the line on a crossplot of P‐wave velocity against S‐wave velocity on which water‐saturated sandstones, shales, and siltstones lie. Some of the rock types that lie off the mudrock line are gas‐saturated sandstones, carbonates, and igneous rocks. In the absence of carbonates and igneous rocks, high amplitude reflections on fluid factor traces would be expected to represent gas‐saturated sandstones. Of course, this relationship does not apply exactly in nature, and the extent to which the mudrock line model applies varies from area to area. However, it is a useful model in many basins of the world, including the one studied here. Geostack processing has been done on a 3-D seismic data set over the Mossel Bay gas field on the southern continental shelf of South Africa. We found that anomalously high amplitude fluid factor reflections occurred at the top and base of the gas‐reservoir sandstone. Maps were made of the amplitude of these fluid factor reflections, and it was found that the high amplitude values were restricted mainly to the gas field area as determined by drilling. The highest amplitudes were found to be located roughly in the areas of best reservoir quality (i.e., highest porosity) in areas where the reservoir is relatively thick.


2021 ◽  
Vol 11 (4) ◽  
pp. 1809-1822
Author(s):  
Alexander Ogbamikhumi ◽  
Osakpolor Marvellous Omorogieva

AbstractThe application of quantitative interpretation techniques for hydrocarbon prospect evaluation from seismic has become so vital. The effective employment of these techniques is dependent on several factors: the quality of the seismic and well data, sparseness of data, the physics of rock, lithological and structural complexity of the field. This study adopts reflection pattern, amplitude versus offset (AVO), Biot–Gassmann fluid substitution and cross-plot models to understand the physics of the reservoir rocks in the field by examining the sensitivity of the basic rock properties; P-wave velocity, S-wave velocity and density, to variation in lithology and fluid types in the pore spaces of reservoirs. This is to ascertain the applicability of quantitative seismic interpretation techniques to explore hydrocarbon prospect in the studied field. The results of reflection pattern and AVO models revealed that the depth of interest is dominated by Class IV AVO sands with a high negative zero offset reflectivity that reduces with offset. The AVO intercept versus gradient plot indicated that both brine and hydrocarbon bearing sands can be discriminated on seismic. Fluid substitution modelling results revealed that the rock properties will favourably respond to variation in oil saturation, but as little as 5% gas presence will result in huge change in the rock properties, which will remain constant upon further increments of gas saturation, thereby making it difficult to differentiate between economical and sub-economical saturations of gas on seismic data. Rock physics cross-plot models revealed separate cluster points typical of shale presence, brine sands and hydrocarbon bearing sands. Thus, the response of the rock properties to the modelling processes adopted favours the application of quantitative interpretation techniques to evaluate hydrocarbon in the field.


Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. B89-B112 ◽  
Author(s):  
G-Akis Tselentis ◽  
Nikolaos Martakis ◽  
Paraskevas Paraskevopoulos ◽  
Athanasios Lois

We have studied using traveltimes of P- and S-waves and initial seismic-pulse rise-time measurements from natural microearthquakes to derive 3D P-wave velocity VP information (mostly structural) as well as P- and S-wave velocity VP/VS and P-wave quality factor QP information (mostly lithologic) in a known hydrocarbon field in southern Albania. During a 12-month monitoring period, 1860 microearthquakes were located at a 50-station seismic network and were used to obtain the above parameters. The data set included earthquakes with magnitudes ranging from –0.1 to 3.0 R (Richter scale) and focal depths typically occurring between 2 and 10 km. Kohonen neural networks were implemented to facilitate the lithological classification of the passive seismic tomography (PST) results. The obtained results, which agreed with data from nearby wells, helped delineate the structure of the reservoir. Two subregions of the investigated area, one corresponding to an oil field and one to a gas field, were correlated with the PST results. This experiment showed that PST is a powerful new geophysical technique for exploring regions that present seismic penetration problems, difficult topographies, and complicated geologies, such as thrust-belt regions. The method is economical and environmentally friendly, and it can be used to investigate very large regions for the optimal design of planned 2D or 3D conventional geophysical surveys.


2021 ◽  
Author(s):  
Wanbo Xiao ◽  
Siqi Lu ◽  
Yanbin Wang

<p>Despite the popularity of the horizontal to vertical spectral ratio (HVSR) method in site effect studies, the origin of the H/V peaks has been controversial since this method was proposed. Many previous studies mainly focused on the explanation of the first or single peak of the H/V ratio, trying to distinguish between the two hypotheses — the S-wave resonance and ellipticity of Rayleigh wave. However, it is common both in numerical simulations and practical experiments that the H/V ratio exhibits multiple peaks, which is essential to explore the origin of the H/V peaks.</p><p>The cause for the multiple H/V peaks has not been clearly figured out, and once was simply explained as the result of multi subsurface layers. Therefore, we adopted numerical method to simulate the ambient noise in various layered half-space models and calculated the H/V ratio curves for further comparisons. The peak frequencies of the H/V curves accord well with the theoretical frequencies of S-wave resonance in two-layer models, whose frequencies only depend on the S wave velocity and the thickness of the subsurface layer. The same is true for models with varying model parameters. Besides, the theoretical formula of the S-wave resonance in multiple-layer models is proposed and then supported by numerical investigations as in the cases of two-layer models. We also extended the S-wave resonance to P-wave resonance and found that its theoretical frequencies fit well with the V/H peaks, which could be an evidence to support the S-wave resonance theory from a new perspective. By contrast, there are obvious differences between the higher orders of the H/V ratio peaks and the higher orders of Rayleigh wave ellipticity curves both in two-layer and multiple-layer models. The Rayleigh wave ellipticity curves are found to be sensitive to the Poisson’s ratio and the thickness of the subsurface layer, so the variation of the P wave velocity can affect the peak frequencies of the Rayleigh wave ellipticity curves while the H/V peaks show slight change. The Rayleigh wave ellipticity theory is thus proved to be inappropriate for the explanation of the multiple H/V peaks, while the possible effects of the Rayleigh wave on the fundamental H/V peak still cannot be excluded.</p><p>Based on the analyses above, we proposed a new evidence to support the claim that the peak frequencies of the H/V ratio curve, except the fundamental peaks, are caused by S-wave resonance. The relationship between the P-wave resonance and the V/H peaks may also find further application.</p>


Sign in / Sign up

Export Citation Format

Share Document