A simple approximation to the P‐wave reflection coefficient and its implication in the inversion of amplitude variation with offset data

Geophysics ◽  
1993 ◽  
Vol 58 (4) ◽  
pp. 544-552 ◽  
Author(s):  
Subhashis Mallick
Keyword(s):  
Reflection Coefficient ◽  
Wave Reflection ◽  
P Wave ◽  
S Wave ◽  
Linear Inversion ◽  
Amplitude Variation With Offset ◽  
Shear Moduli ◽  
Amplitude Versus Offset ◽  
Class 1 ◽  
Two Parameters

I derive an approximate formula for the plane P‐wave reflection coefficient as a function of ray‐parameter. The approximation shows that the behavior of the P‐wave reflection coefficient at nonnormal angles of incidence is mainly controlled by two parameters: (1) [Formula: see text], the fluid‐fluid reflection coefficient (i.e., the reflection coefficient when the S‐wave velocities in both media are set to zero) and (2) Δμ/ρ, the ratio of the contrast in shear moduli to the average bulk density. I also show that the other formulas for the P‐wave reflection coefficient given by R. Bortfeld and R. Shuey can be approximately derived from my formula. I give numerical examples to demonstrate the accuracy of the formula. Using a least‐squares inversion of theoretical values for the reflection coefficients, I demonstrate that, in a linear inversion of amplitude versus offset data, Δμ/ρ is better estimated than the contrast in Poisson’s ratio, Δσ. Finally, comparing the exact reflection coefficient with Shuey’s approximation for typical shale and Class 1 gas‐sand models, I show that Shuey’s approximation under‐estimates the value of |Δσ| for such reflections.

Two-step joint PP- and PS-wave three-term amplitude-variation with offset inversion

2016 ◽  
Vol 4 (4) ◽  
pp. T613-T625 ◽  
Author(s):  
Qizhen Du ◽  
Bo Zhang ◽  
Xianjun Meng ◽  
Chengfeng Guo ◽  
Gang Chen ◽  
...  
Keyword(s):  
Reflection Coefficient ◽  
Wave Reflection ◽  
Coefficient Matrix ◽  
P Wave ◽  
Inversion Method ◽  
S Wave ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
Avo Inversion ◽  
Pp Wave

Three-term amplitude-variation with offset (AVO) inversion generally suffers from instability when there is limited prior geologic or petrophysical constraints. Two-term AVO inversion shows higher instability compared with three-term AVO inversion. However, density, which is important in the fluid-type estimation, cannot be recovered from two-term AVO inversion. To reliably predict the P- and S-waves and density, we have developed a robust two-step joint PP- and PS-wave three-term AVO-inversion method. Our inversion workflow consists of two steps. The first step is to estimate the P- and S-wave reflectivities using Stewart’s joint two-term PP- and PS-AVO inversion. The second step is to treat the P-wave reflectivity obtained from the first step as the prior constraint to remove the P-wave velocity related-term from the three-term Aki-Richards PP-wave approximated reflection coefficient equation, and then the reduced PP-wave reflection coefficient equation is combined with the PS-wave reflection coefficient equation to estimate the S-wave and density reflectivities. We determined the effectiveness of our method by first applying it to synthetic models and then to field data. We also analyzed the condition number of the coefficient matrix to illustrate the stability of the proposed method. The estimated results using proposed method are superior to those obtained from three-term AVO inversion.

Application of instantaneous rotations to S‐wave vertical seismic profiling

Geophysics ◽  
1997 ◽  
Vol 62 (5) ◽  
pp. 1365-1368
Author(s):  
M. Boulfoul ◽  
Doyle R. Watts
Keyword(s):  
High Amplitude ◽  
Seismic Profile ◽  
P Wave ◽  
Petroleum Exploration ◽  
S Wave ◽  
Amplitude Variation With Offset ◽  
Seismic Profiling ◽  
Vertical Seismic Profiling ◽  
Seismic Models ◽  
Low Amplitude

The petroleum exploration industry uses S‐wave vertical seismic profiling (VSP) to determine S‐wave velocities from downgoing direct arrivals, and S‐wave reflectivities from upgoing waves. Seismic models for quantitative calibration of amplitude variation with offset (AVO) data require S‐wave velocity profiles (Castagna et al., 1993). Vertical summations (Hardage, 1983) of the upgoing waves produce S‐wave composite traces and enable interpretation of S‐wave seismic profile sections. In the simplest application of amplitude anomalies, the coincidence of high amplitude P‐wave reflectivity and low amplitude S‐wave reflectivity is potentially a direct indicator of the presence of natural gas.

Low‐cost geophysical investigations of a paleochannel aquifer in the Eastern Goldfields, Western Australia

Geophysics ◽  
2002 ◽  
Vol 67 (3) ◽  
pp. 690-700 ◽  
Author(s):  
Josef Holzschuh
Keyword(s):  
Water Table ◽  
Western Australia ◽  
Seismic Reflection ◽  
Wave Reflection ◽  
Gravity Data ◽  
P Wave ◽  
Gravity Modeling ◽  
S Wave ◽  
Sand And Gravel ◽  
Gravel Aquifer

Compressional (P) wave and shear (S) wave seismic reflection techniques were used to delineate the sand and gravel aquifer within a highly saline clay‐filled paleochannel in the Eastern Goldfields of Western Australia. The seismic refraction and gravity methods were also used to investigate the paleochannel. The unsaturated loose fine‐grained sand up to 10 m in depth at the surface is a major factor in degrading subsurface imaging. The seismic processing needed to be precise, with accurate static corrections and normal moveout corrections. Deconvolution enhanced the aquifer and other paleochannel reflectors. P‐wave reflection and refraction layer depths had good correlation and showed a total of six boundaries: (1) water table, (2) change in velocity (compaction) in the paleochannel sediments, (3) sand and gravel aquifer, (4) red‐brown saprolite and green saprolite boundary, (5) weathered bedrock, and (6) unweathered bedrock. P‐wave explosive and hammer sources were found to have similar signal characteristics, and the aquifer and bedrock were both imaged using the hammer source. The deep shots below the water table have the most broadband frequency response for reflections, but stacking clear reflections was difficult. The S‐wave reflection results showed high lateral and vertical resolution of the basal saprolite clay, the sand and gravel aquifer, and very shallow clays above the aquifer. The S‐wave reflection stacking velocities were 10–20% of the P‐waves, increasing the resolution of the S‐wave section. The gravity data were modelled to fit the known drilling and P‐wave seismic reflection depths. The refraction results did not identify the top of bedrock, so refraction depths were not used for the gravity modeling in this highly weathered environment. The final gravity model mapped the bedrock topography beyond the lateral extent of the seismic and drilling data.

Three-parameter prestack seismic inversion based on L1-2 minimization

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. R753-R766 ◽  
Author(s):  
Lingqian Wang ◽  
Hui Zhou ◽  
Yufeng Wang ◽  
Bo Yu ◽  
Yuanpeng Zhang ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Seismic Inversion ◽  
Thin Layers ◽  
P Wave ◽  
S Wave ◽  
Reservoir Prediction ◽  
Linear Inversion ◽  
Prestack Inversion ◽  
The Difference ◽  
Predictive Filtering

Prestack inversion has become a common approach in reservoir prediction. At present, the critical issue in the application of seismic inversion is the estimation of elastic parameters in the thin layers and weak reflectors. To improve the resolution and the accuracy of the inversion results, we introduced the difference of [Formula: see text] and [Formula: see text] norms as a nearly unbiased approximation of the sparsity of a vector, denoted as the [Formula: see text] norm, to the prestack inversion. The nonconvex penalty function of the [Formula: see text] norm can be decomposed into two convex subproblems via the difference of convex algorithm, and each subproblem can be solved efficiently by the alternating direction method of multipliers. Compared with the [Formula: see text] norm regularization, the [Formula: see text] minimization can reconstruct reflectivities more accurately. In addition, the [Formula: see text]-[Formula: see text] predictive filtering was introduced to guarantee the lateral continuity of the location and the amplitude of the reflectivity series. The generalized linear inversion and [Formula: see text]-[Formula: see text] predictive filtering are combined for stable elastic impedance inversion results, and three parameters of P-wave velocity, S-wave velocity, and density can be inverted with the Bayesian linearized amplitude variation with offset inversion. The inversion results of synthetic and real seismic data demonstrate that the proposed method can effectively improve the resolution and accuracy of the inversion results.

Frequency-dependent spherical-wave reflection coefficient inversion in acoustic media: Theory to practice

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. R425-R435
Author(s):  
Binpeng Yan ◽  
Shangxu Wang ◽  
Yongzhen Ji ◽  
Xingguo Huang ◽  
Nuno V. da Silva
Keyword(s):  
Reflection Coefficient ◽  
Wave Reflection ◽  
Incident Angle ◽  
Spherical Wave ◽  
Frequency Dependent ◽  
Amplitude Variation With Offset ◽  
Reflection Data ◽  
Alternative Approach ◽  
Acoustic Media ◽  
The Time Domain

As an approximation of the spherical-wave reflection coefficient (SRC), the plane-wave reflection coefficient does not fully describe the reflection phenomenon of a seismic wave generated by a point source. The applications of SRC to improve analyses of seismic data have also been studied. However, most of the studies focus on the time-domain SRC and its benefit to using the long-offset information instead of the dependency of SRC on frequency. Consequently, we have investigated and accounted for the frequency-dependent spherical-wave reflection coefficient (FSRC) and analyzed the feasibility of this type of inversion. Our inversion strategy requires a single incident angle using reflection data for inverting the density and velocity ratios, which is distinctly different from conventional inversion methods using amplitude variation with offset. Hence, this investigation provides an alternative approach for estimating media properties in some contexts, especially when the range of aperture of the reflection angles is limited. We apply the FSRC theory to the inversion of noisy synthetic and field data using a heuristic algorithm. The multirealization results of the inversion strategy are consistent with the feasibility analysis and demonstrate the potential of the outlined method for practical application.

Bayesian closed-skew Gaussian inversion of seismic AVO data for elastic material properties

Geophysics ◽  
2010 ◽  
Vol 75 (1) ◽  
pp. R1-R11 ◽  
Author(s):  
Omid Karimi ◽  
Henning Omre ◽  
Mohsen Mohammadzadeh
Keyword(s):  
Elastic Material ◽  
Material Properties ◽  
P Wave ◽  
S Wave ◽  
High Dimensions ◽  
Seismic Amplitude ◽  
The North ◽  
Avo Inversion ◽  
Amplitude Versus Offset ◽  
The North Sea

Bayesian closed-skew Gaussian inversion is defined as a generalization of traditional Bayesian Gaussian inversion, which is used frequently in seismic amplitude-versus-offset (AVO) inversion. The new model captures skewness in the variables of interest; hence, the posterior model for log-transformed elastic material properties given seismic AVO data might be a skew probability density function. The model is analytically tractable, and this makes it applicable in high-dimensional 3D inversion problems. Assessment of the posterior models in high dimensions requires numerical approximations, however. The Bayesian closed-skew Gaussian inversion approach has been applied on real elastic material properties from a well in the Sleipner field in the North Sea. A comparison with results from traditional Bayesian Gaussian inversion shows that the mean square error of predictions of P-wave and S-wave velocities are reduced by a factor of two, although somewhat less for density predictions.

Seismic wave attenuation and modulus dispersion in sandstones

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. D211-D231 ◽  
Author(s):  
James W. Spencer ◽  
Jacob Shine
Keyword(s):  
Wave Attenuation ◽  
High Viscosity ◽  
P Wave ◽  
S Wave ◽  
Local Flow ◽  
Shear Moduli ◽  
Low Viscosity ◽  
Confining Pressures ◽  
Wave Induced ◽  
Dispersion And Attenuation

We have conducted laboratory experiments over the 1–200 Hz band to examine the effects of viscosity and permeability on modulus dispersion and attenuation in sandstones and also to examine the effects of partial gas or oil saturation on velocities and attenuations. Our results have indicated that bulk modulus values with low-viscosity fluids are close to the values predicted using Gassmann’s first equation, but, with increasing frequency and viscosity, the bulk and shear moduli progressively deviate from the values predicted by Gassmann’s equations. The shear moduli increase up to 1 GPa (or approximately 10%) with high-viscosity fluids. The P- and S-wave attenuations ([Formula: see text] and [Formula: see text]) and modulus dispersion with different fluids are indicative of stress relaxations that to the first order are scaling with frequency times viscosity. By fitting Cole-Cole distributions to the scaled modulus and attenuation data, we have found that there are similar P-wave, shear and bulk relaxations, and attenuation peaks in each of the five sandstones studied. The modulus defects range from 11% to 15% in Berea sandstone to 16% to 26% in the other sandstones, but these would be reduced at higher confining pressures. The relaxations shift to lower frequencies as the viscosity increased, but they do not show the dependence on permeability predicted by mesoscopic wave-induced fluid flow (WIFF) theories. Results from other experiments having patchy saturation with liquid [Formula: see text] and high-modulus fluids are consistent with mesoscopic WIFF theories. We have concluded that the modulus dispersion and attenuations ([Formula: see text] and [Formula: see text]) in saturated sandstones are caused by a pore-scale, local-flow mechanism operating near grain contacts.

Frequency-dependent spherical-wave nonlinear AVO inversion in elastic media

2020 ◽  
Vol 223 (2) ◽  
pp. 765-776
Author(s):  
Guangsen Cheng ◽  
Xingyao Yin ◽  
Zhaoyun Zong
Keyword(s):  
Reflection Coefficient ◽  
Wave Reflection ◽  
Fresnel Zone ◽  
Spherical Wave ◽  
Inversion Method ◽  
Elastic Media ◽  
Amplitude Variation With Offset ◽  
Avo Inversion ◽  
Well Location ◽  
The Stability

SUMMARY The plane-wave reflection coefficient (PRC) plays a remarkable role in conventional amplitude variation with offset (AVO) analysis and inversion. Compared with the widely exploited PRC that breaks down at the near- and supercritical incidence angles, the spherical-wave reflection coefficient (SRC) can overcome the influence of wide-angle reflection and give an accurate description of the actual seismic wave reflection phenomenon based on spherical-wave fronts. However, SRC is not widely used in AVO inversion due to its nonlinearity and computational complexity. In our study, the characteristics of frequency–depth-dependent monochromatic SRC are discussed and a novel three-parameter SRC is derived. Compared with the conventional six-parameter SRC, the novel three-parameter SRC improves the stability of spherical-wave AVO inversion. In addition, the concept of SRC within the Fresnel zone is proposed, and the accuracy of SRC within the Fresnel zone in the deep subsurface is tested. Finally, a nonlinear spherical-wave AVO inversion method for elastic media is proposed, which can make full use of all frequency components of wavelet. The robustness of the proposed method is verified by the application on synthetic seismogram with white Gaussian noise. The feasibility and practicability of this method are verified by comparing the spherical-wave AVO inversion results with the filtered well logs at the known well location.

The effects of pore‐fluid salinity on ultrasonic wave propagation in sandstones

Geophysics ◽  
1998 ◽  
Vol 63 (3) ◽  
pp. 928-934 ◽  
Author(s):  
Simon M. Jones ◽  
Clive McCann ◽  
Timothy R. Astin ◽  
Jeremy Sothcott
Keyword(s):  
Wave Propagation ◽  
Ultrasonic Wave ◽  
Wave Attenuation ◽  
Clay Content ◽  
Seismic Wave Propagation ◽  
Ultrasonic Pulse ◽  
Pore Fluid ◽  
P Wave ◽  
S Wave ◽  
Shear Moduli

Petrophysical interpretation of increasingly refined seismic data from subsurface formations requires a more fundamental understanding of seismic wave propagation in sedimentary rocks. We consider the variation of ultrasonic wave velocity and attenuation in sandstones with pore‐fluid salinity and show that wave propagation is modified in proportion to the clay content of the rock and the salinity of the pore fluid. Using an ultrasonic pulse reflection technique (590–890 kHz), we have measured the P-wave and S-wave velocities and attenuations of 15 saturated sandstones with variable effective pressure (5–60 MPa) and pore‐fluid salinity (0.0–3.4 M). In clean sandstones, there was close agreement between experimental and Biot model values of [Formula: see text], but they diverged progressively in rocks containing more than 5% clay. However, this effect is small: [Formula: see text] changed by only 0.6% per molar change in salinity for a rock with a clay content of 29%. The variation of [Formula: see text] with brine molarity exhibited Biot behavior in some samples but not in others; there was no obvious relationship with clay content. P-wave attenuation was independent of pore‐fluid salinity, while S-wave attenuation was weakly dependent. The velocity data suggest the frame bulk and shear moduli of sandstones are altered by changes in the pore‐fluid salinity. One possible mechanism is the formation damage caused by clay swelling and migration of fines in low‐molarity electrolytes. The absence of variation between the attenuation in water‐saturated and brine‐saturated samples indicates the attenuation mechanism is relatively unaffected by changes in the frame moduli.

Sign in / Sign up

Export Citation Format

Share Document