Converted-wave model building and imaging based on common-focus-point methodology

Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. U139-U149
Author(s):  
Hongwei Liu ◽  
Mustafa Naser Al-Ali ◽  
Yi Luo
Keyword(s):  
Elastic Wave ◽  
Wave Velocity ◽  
Seismic Data ◽  
Model Building ◽  
P Wave ◽  
Rock Properties ◽  
S Wave ◽  
Converted Wave ◽  
Seismic Images ◽  
Focus Point

Seismic images can be viewed as photographs for underground rocks. These images can be generated from different reflections of elastic waves with different rock properties. Although the dominant seismic data processing is still based on the acoustic wave assumption, elastic wave processing and imaging have become increasingly popular in recent years. A major challenge in elastic wave processing is shear-wave (S-wave) velocity model building. For this reason, we have developed a sequence of procedures for estimating seismic S-wave velocities and the subsequent generation of seismic images using converted waves. We have two main essential new supporting techniques. The first technique is the decoupling of the S-wave information by generating common-focus-point gathers via application of the compressional-wave (P-wave) velocity on the converted seismic data. The second technique is to assume one common VP/ VS ratio to approximate two types of ratios, namely, the ratio of the average earth layer velocity and the ratio of the stacking velocity. The benefit is that we reduce two unknown ratios into one, so it can be easily scanned and picked in practice. The PS-wave images produced by this technology could be aligned with the PP-wave images such that both can be produced in the same coordinate system. The registration between the PP and PS images provides cross-validation of the migrated structures and a better estimation of underground rock and fluid properties. The S-wave velocity, computed from the picked optimal ratio, can be used not only for generating the PS-wave images, but also to ensure well registration between the converted-wave and P-wave images.

Useful approximations for converted‐wave AVO

Geophysics ◽  
2001 ◽  
Vol 66 (6) ◽  
pp. 1721-1734 ◽  
Author(s):  
Antonio C. B. Ramos ◽  
John P. Castagna
Keyword(s):  
Wave Velocity ◽  
Random Noise ◽  
P Wave ◽  
Rock Properties ◽  
Negative Slope ◽  
High Contrast ◽  
Fractional Change ◽  
S Wave ◽  
Converted Wave ◽  
Low Contrast

Converted‐wave amplitude versus offset (AVO) behavior may be fit with a cubic relationship between reflection coefficient and ray parameter. Attributes extracted using this form can be directly related to elastic parameters with low‐contrast or high‐contrast approximations to the Zoeppritz equations. The high‐contrast approximation has the advantage of greater accuracy; the low‐contrast approximation is analytically simpler. The two coefficients of the low‐contrast approximation are a function of the average ratio of compressional‐to‐shear‐wave velocity (α/β) and the fractional changes in S‐wave velocity and density (Δβ/β and Δρ/ρ). Because of its simplicity, the low‐contrast approximation is subject to errors, particularly for large positive contrasts in P‐wave velocity associated with negative contrasts in S‐wave velocity. However, for incidence angles up to 40° and models confined to |Δβ/β| < 0.25, the errors in both coefficients are relatively small. Converted‐wave AVO crossplotting of the coefficients of the low‐contrast approximation is a useful interpretation technique. The background trend in this case has a negative slope and an intercept proportional to the α/β ratio and the fractional change in S‐wave velocity. For constant α/β ratio, an attribute trace formed by the weighted sum of the coefficients of the low‐contrast approximation provides useful estimates of the fractional change in S‐wave velocity and density. Using synthetic examples, we investigate the sensitivity of these parameters to random noise. Integrated P‐wave and converted‐wave analysis may improve estimation of rock properties by combining extracted attributes to yield fractional contrasts in P‐wave and S‐wave velocities and density. Together, these parameters may provide improved direct hydrocarbon indication and can potentially be used to identify anomalies caused by low gas saturations.

Estimating shear‐wave velocities from P-wave and converted‐wave data

Geophysics ◽  
1999 ◽  
Vol 64 (2) ◽  
pp. 504-507 ◽  
Author(s):  
Franklyn K. Levin
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
P Wave ◽  
S Wave ◽  
Converted Wave ◽  
Wave Data ◽  
Wave Velocities ◽  
Shear Wave Velocities ◽  
P Wave Velocity ◽  
Growing Body

Tessmer and Behle (1988) show that S-wave velocity can be estimated from surface seismic data if both normal P-wave data and converted‐wave data (P-SV) are available. The relation of Tessmer and Behle is [Formula: see text] (1) where [Formula: see text] is the S-wave velocity, [Formula: see text] is the P-wave velocity, and [Formula: see text] is the converted‐wave velocity. The growing body of converted‐wave data suggest a brief examination of the validity of equation (1) for velocities that vary with depth.

Interpretation of AVO anomalies

Geophysics ◽  
2010 ◽  
Vol 75 (5) ◽  
pp. 75A3-75A13 ◽  
Author(s):  
Douglas J. Foster ◽  
Robert G. Keys ◽  
F. David Lane
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
P Wave ◽  
Rock Properties ◽  
Light Hydrocarbons ◽  
S Wave ◽  
Amplitude Variation With Offset ◽  
P Wave Velocity ◽  
Fluid Properties ◽  
Porosity Changes

We investigate the effects of changes in rock and fluid properties on amplitude-variation-with-offset (AVO) responses. In the slope-intercept domain, reflections from wet sands and shales fall on or near a trend that we call the fluid line. Reflections from the top of sands containing gas or light hydrocarbons fall on a trend approximately parallel to the fluid line; reflections from the base of gas sands fall on a parallel trend on the opposing side of the fluid line. The polarity standard of the seismic data dictates whether these reflections from the top of hydrocarbon-bearing sands are below or above the fluid line. Typically, rock properties of sands and shales differ, and therefore reflections from sand/shale interfaces are also displaced from the fluid line. The distance of these trends from the fluid line depends upon the contrast of the ratio of P-wave velocity [Formula: see text] and S-wave velocity [Formula: see text]. This ratio is a function of pore-fluid compressibility and implies that distance from the fluid line increases with increasing compressibility. Reflections from wet sands are closer to the fluid line than hydrocarbon-related reflections. Porosity changes affect acoustic impedance but do not significantly impact the [Formula: see text] contrast. As a result, porosity changes move the AVO response along trends approximately parallel to the fluid line. These observations are useful for interpreting AVO anomalies in terms of fluids, lithology, and porosity.

Using Seismic Data to Predict Shale Pore Pressure and Overpressure Sweet Spots: A Case Study from the Lower Silurian Longmaxi Formation in Sichuan Basin, China

2021 ◽  
Author(s):  
Sheng Chen ◽  
Qingcai Zeng ◽  
Xiujiao Wang ◽  
Qing Yang ◽  
Chunmeng Dai ◽  
...  
Keyword(s):  
Prediction Model ◽  
Wave Velocity ◽  
Shale Gas ◽  
Seismic Data ◽  
P Wave ◽  
Formation Pressure ◽  
S Wave ◽  
P Wave Velocity ◽  
New Formation ◽  
Sweet Spots

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.

Elastic reflectivity vectors and the geometry of intercept-gradient crossplots

2019 ◽  
Vol 38 (10) ◽  
pp. 762-769
Author(s):  
Patrick Connolly
Keyword(s):  
Elastic Property ◽  
Elastic Properties ◽  
Wave Velocity ◽  
Seismic Data ◽  
Measurement Errors ◽  
Rock Physics ◽  
P Wave ◽  
S Wave ◽  
P Wave Velocity ◽  
Well Log Analysis

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.

Simultaneous inversion of PP and PS seismic data

Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. R1-R10 ◽  
Author(s):  
Helene Hafslund Veire ◽  
Martin Landrø
Keyword(s):  
Seismic Data ◽  
Reservoir Characterization ◽  
Model Building ◽  
Joint Inversion ◽  
P Wave ◽  
Inversion Method ◽  
System Of Equations ◽  
S Wave ◽  
Data Set ◽  
Simultaneous Inversion

Elastic parameters derived from seismic data are valuable input for reservoir characterization because they can be related to lithology and fluid content of the reservoir through empirical relationships. The relationship between physical properties of rocks and fluids and P-wave seismic data is nonunique. This leads to large uncertainties in reservoir models derived from P-wave seismic data. Because S- waves do not propagate through fluids, the combined use of P-and S-wave seismic data might increase our ability to derive fluid and lithology effects from seismic data, reducing the uncertainty in reservoir characterization and thereby improving 3D reservoir model-building. We present a joint inversion method for PP and PS seismic data by solving approximated linear expressions of PP and PS reflection coefficients simultaneously using a least-squares estimation algorithm. The resulting system of equations is solved by singular-value decomposition (SVD). By combining the two independent measurements (PP and PS seismic data), we stabilize the system of equations for PP and PS seismic data separately, leading to more robust parameter estimation. The method does not require any knowledge of PP and PS wavelets. We tested the stability of this joint inversion method on a 1D synthetic data set. We also applied the methodology to North Sea multicomponent field data to identify sand layers in a shallow formation. The identified sand layers from our inverted sections are consistent with observations from nearby well logs.

Wavelet distortion correction due to domain conversion

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. V77-V87 ◽  
Author(s):  
Rishi Bansal ◽  
Mike Matheney
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
Seismic Velocity ◽  
Distortion Correction ◽  
P Wave ◽  
Converted Wave ◽  
Shaping Filter ◽  
Seismic Wavelet ◽  
Attribute Analysis ◽  
P Wave Velocity

Converted-wave (PS) data, when converted to PP time, develop time- and location-varying compression of the seismic wavelet due to a variable subsurface [Formula: see text] [Formula: see text]. The time-dependent compression distorts the wavelet in a seismic trace. The lack of a consistent seismic wavelet in a domain-converted PS volume can eventually lead to an erroneous joint PP/PS inversion result. Depth-converted seismic data also have wavelet distortion due to velocity-dependent wavelet stretch. A high value of seismic velocity produces more stretch in a seismic wavelet than a low value. Variable wavelet stretch renders the depth data unsuitable for attribute analysis. A filtering scheme is proposed that corrects for distortion in seismic wavelets due to domain conversions (PS to PP time and time-to-depth) of seismic data in an amplitude-preserving manner. The method uses a Fourier scaling theorem to predict the seismic wavelet in the converted domain and calculates a shaping filter for each time/depth sample that corrects for the distortion in the wavelet. The filter is applied to the domain-converted data using the method of nonstationary filtering. We provide analytical expressions for the squeeze factor [Formula: see text] that is used to predict the wavelet in the converted domain. The squeeze factor [Formula: see text] for PS to PP time conversion is a function of the subsurface [Formula: see text] whereas for PP time-to-depth conversion [Formula: see text] is dependent on subsurface P-wave velocity. After filtering, the squeezed wavelets in domain-converted PS data appear to have resulted from a constant subsurface [Formula: see text], which we denote as [Formula: see text]. Similarly, the filtered depth-converted data appear to have resulted from a constant subsurface P-wave velocity [Formula: see text].

Sediments with gas hydrates: Internal structure from seismic AVO

Geophysics ◽  
1998 ◽  
Vol 63 (5) ◽  
pp. 1659-1669 ◽  
Author(s):  
Christine Ecker ◽  
Jack Dvorkin ◽  
Amos Nur
Keyword(s):  
Internal Structure ◽  
Wave Velocity ◽  
Seismic Data ◽  
Pore Space ◽  
Rock Physics ◽  
P Wave ◽  
S Wave ◽  
High Concentration ◽  
Amplitude Variation With Offset ◽  
Bottom Simulating Reflector

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.

Predicting reservoir quality in the Bakken Formation, North Dakota, using petrophysics and 3C seismic data

2020 ◽  
Vol 8 (4) ◽  
pp. T851-T868
Author(s):  
Andrea G. Paris ◽  
Robert R. Stewart
Keyword(s):  
North Dakota ◽  
Seismic Data ◽  
Rock Physics ◽  
P Wave ◽  
Well Logs ◽  
Reservoir Quality ◽  
Bakken Formation ◽  
S Wave ◽  
Converted Wave ◽  
Wave Velocities

Combining rock-property analysis with multicomponent seismic imaging can be an effective approach for reservoir quality prediction in the Bakken Formation, North Dakota. The hydrocarbon potential of shale is indicated on well logs by low density, high gamma-ray response, low compressional-wave (P-wave) and shear-wave (S-wave) velocities, and high neutron porosity. We have recognized the shale intervals by cross plotting sonic velocities versus density. Intervals with total organic carbon (TOC) content higher than 10 wt% deviate from lower TOC regions in the density domain and exhibit slightly lower velocities and densities (<2.30 g/cm3). We consider TOC to be the principal factor affecting changes in the density and P- and S-wave velocities in the Bakken shales, where VP/ VS ranges between 1.65 and 1.75. We generate the synthetic seismic data using an anisotropic version of the Zoeppritz equations, including estimated Thomsen’s parameters. For the tops of the Upper and Lower Bakken, the amplitude shows a negative intercept and a positive gradient, which corresponds to an amplitude variation with offset of class IV. The P-impedance error decreases by 14% when incorporating the converted-wave information in the inversion process. A statistical approach using multiattribute analysis and neural networks delimits the zones of interest in terms of P-impedance, density, TOC content, and brittleness. The inverted and predicted results show reasonable correlations with the original well logs. The integration of well log analysis, rock physics, seismic modeling, constrained inversions, and statistical predictions contributes to identifying the areas of highest reservoir quality within the Bakken Formation.

scholarly journals Identification of hydrocarbons in chalk reservoirs from surface seismic data: South Arne field, North Sea

2005 ◽  
Vol 7 ◽  
pp. 13-16
Author(s):  
Peter Japsen ◽  
Anders Bruun ◽  
Ida L. Fabricius ◽  
Gary Mavko
Keyword(s):  
Wave Velocity ◽  
North Sea ◽  
Seismic Data ◽  
Pore Space ◽  
Sedimentary Rocks ◽  
Rock Physics ◽  
Acoustic Properties ◽  
P Wave ◽  
S Wave ◽  
Reservoir Rocks

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

Sign in / Sign up

Export Citation Format

Share Document