scholarly journals Extended poroelastic impedance

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. N1-N14 ◽  
Author(s):  
Brian H. Russell ◽  
Ken J. Hedlin

Linearized approximations to the P-wave reflectivity as a function of the incidence angle (called amplitude variation with offset) involve the extraction of band-limited reflectivity terms that are a function of changes in the elastic constants of the earth across each lithologic interface. The most common of these extracted reflectivities are the intercept and gradient, usually labeled [Formula: see text] and [Formula: see text], respectively. The extended elastic impedance (EEI) method uses a rotation angle [Formula: see text] to map [Formula: see text] and [Formula: see text] into a new reflectivity corresponding to a particular elastic parameter. The success of EEI depends on finding an optimum value for the angle [Formula: see text]. This value is usually calculated by correlating the EEI result over a range of [Formula: see text] angles with various elastic parameters and then finding the best correlation coefficient. We have developed a new approach for the interpretation of the EEI method, which incorporates the Biot-Gassmann poroelastic theory and attaches a physical meaning to the [Formula: see text] angle. We call this method extended poroelastic impedance (EPI). The main advantage of the EPI method is that the [Formula: see text] angle is now interpreted as a parameter that is dependent on the dry-rock properties of the reservoir, rather than a parameter whose value is estimated empirically. The method is evaluated by numerical and synthetic seismic examples and by application to field data from a gas sand reservoir.

Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. C1-C7 ◽  
Author(s):  
Subhashis Mallick

Amplitude-variation-with-offset (AVO) and elastic-impedance (EI) analysis use an approximate plane P-wave reflection coefficient as a function of angle of incidence. AVO and EI both can be used in a three-term or a two-term formulation. This study uses synthetic data to demonstrate that the P-wave primary reflections at large offsets can be contaminated by reflections from other wave modes that can affect the quality of three-term AVO or EI results. The coupling of P-waves and S-waves in seismic-wave propagation through finely layered media generates the interfering wave modes. A methodology such as prestack-wave-equation modeling can properly account for these coupling effects. Both AVO and EI also assume a convolutional model whose accuracy decreases as incidence angles increase. On the other hand, wave-equation modeling is based on the rigorous solution to the wave equation and is valid for any incidence angle. Because wave interference is minimal at small angles, a two-term AVO/EI analysis that restricts input from small angles is likely to give more reliable parameter estimates than a three-term analysis. A three-term AVO/EI analysis should be used with caution and should be calibrated against well data and other data before being used for quantitative analysis.


Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. N21-N33 ◽  
Author(s):  
Kristopher A. Innanen

Linear amplitude-variation-with-offset (AVO) approximations, which experience a reduction in accuracy as elastic parameter contrasts become large, may be adjusted with second- and higher-order corrections. Corrective terms can be expressed in many ways, but they only serve a meaningful purpose if they provide the same qualitative interpretability as did the linearization. Some aspects of nonlinear AVO can be understood, quantitatively and qualitatively, in terms of coupling — the interdependence of elastic parameter contrasts amongst themselves in their determination of reflection strengths. Coupling, for instance, explains the weak but nonnegligible dependence of the converted wave reflection coefficient on the lower half-space P-wave velocity. This fact can be exposed by expanding the solutions of the Zoeppritz equations in a particular hierarchy of series. Also explainable through this approach is the mathematical importance of what is sometimes referred to as the “Wiggins approximation,” under which [Formula: see text]. This special number is seen to coincide with a full decoupling of density contrasts from [Formula: see text] and [Formula: see text] contrasts at the second order. The decoupling persists across several variations of the nonlinear AVO approximations, including both expressions in terms of the relative changes [Formula: see text], [Formula: see text], and [Formula: see text], and expressions in terms of single-parameter reflectivities.


2018 ◽  
Vol 22 (3) ◽  
pp. 205-213 ◽  
Author(s):  
Saiq Shakeel Abbasi ◽  
Jiangping Liu ◽  
Nayima Hameed ◽  
Muhsan Ehsan

Elastic impedance inversion is the latest development in the field of hydrocarbon exploration and production. The present research focuses on the improvement of the use of elastic impedance inversion, easing exploration of hydrocarbons. The seismic velocities change with variation in geological constraints. Constant K, which is S-wave to P-wave ratio of the nth layer and n+1 layer across the interface, it must be changed accordingly. This research focuses on testing the effects of K as a constant in the elastic impedance equation. As using the same value of K for all types of formations can give rise to severe errors in the interpretation of data. The importance of the value of K for particular Amplitude Variation with Offset AVO type (I-IV) is studied using different Elastic Impedance Equations. The Reflection Coefficient (RC) curves for each AVO class are generated using Zoeppritz approximation and Elastic Impedance equations. The comparison of RC curves shows significant variations at far offsets in each AVO type using the Constant value of K. When K Calculated is used, AVO type I and Type II shows a good match at near, mid and far offsets. Type III does not change due to the changing value of K. Type IV gives good agreement at near and intermediate offsets. This variation in curves, with the change in the value of K, indicates that it is a significant factor of interpretation using elastic impedance. The application of findings on well logs has given a satisfactory confirmation of the present results. This research can be helpful to resolve severe errors in the interpretation due to the constant value of K.


2015 ◽  
Vol 3 (3) ◽  
pp. SZ15-SZ26 ◽  
Author(s):  
Stephanie Tyiasning ◽  
Dennis Cooke

We have developed a tight gas amplitude variation with offset (AVO) case history from the Cooper Basin of Australia that addressed the exploration problem of mapping thin fluvial tight gas sand bodies. In the Cooper Basin, Permian Toolachee and Patchawarra sands are difficult to interpret on seismic data due to strong reflections from adjacent Permian coals. This is not the common AVO problem of distinguishing between coal and gas sand, but a more difficult class-I AVO problem of mapping fluvial sands beneath a sheet coal that varies in thickness. We have reviewed local rock properties and concluded that Poisson’s ratio is probably the most appropriate rock property to solve the above exploration problem. We have compared various seismic attributes made using the extended elastic impedance (EEI) technique and a rotation of near and far partial stacks. In a synthetic modeling study that included random noise and tuning, we compared the noise-discrimination abilities of three competing AVO crossplot techniques and “rotated” the attributes made from them. These three crossplots were as follows: intercept versus gradient (I-G), full-stack versus far-minus-near (Full-FmN), and near-stack versus far-stack (N-F). Previous papers on this subject have found that (I-G) crossplots had a spurious correlation in the presence of noise that did not occur with the (Full-FmN) and (N-F) crossplots. We found that for our class-I AVO case, (1) the advantage of the (Full-FmN) and (N-F) crossplots disappeared in the presence of tuning, (2) if tuning was present, the optimal rotation angle was determined by the “tuning angle,” not by the noise angle or some desired EEI angle, and (3) if the three different crossplots were rotated by their respective “tuning” angles, the results were identical.


Geophysics ◽  
2021 ◽  
pp. 1-145
Author(s):  
Xiaobo Liu ◽  
Jingyi Chen ◽  
Jing Zeng ◽  
Fuping Liu ◽  
Handong Huang ◽  
...  

Amplitude variation with incidence angle (AVA) analysis is an essential tool for discriminating lithology in the hydrocarbon reservoirs. Compared with the traditional AVA inversion using only P-wave information, joint AVA inversion using PP and PS seismic data provides better estimation of rock properties (e.g., density, P- and S-wave velocities). At present, the most used AVA inversions depend on the approximations of Zoeppritz equations (e.g., Shuey and Aki-Richards approximations), which are not suitable for formations with strong contrast interfaces and seismic data with large incidence angles. Based on the previous derivation of accurate Jacobian matrix, we find that the sign of each partial derivative of reflection coefficient with respect to P-, S-wave velocities and density changes across the interface, represents good indicator for the reflection interfaces. Accordingly, we propose an adaptive stratified joint PP and PS AVA inversion using the accurate Jacobian matrix that can automatically obtain the layer information and can be further used as a constraint in the inversion of in-layer rock properties (density, P- and S-wave velocities). Due to the use of the exact Zoeppritz equations and accurate Jacobian matrix, this proposed inversion method is more accurate than traditional AVA inversion methods, has higher computational efficiency and can be applied to seismic wide-angle reflection data or seismic data acquired for formations with strong contrast interfaces. The model study shows that this proposed inversion method works better than the classical Shuey and Aki-Richards approximations at estimating reflection interfaces and in-layer rock properties. It also works well in handling a part of the complex Marmousi 2 model and real seismic data.


2018 ◽  
Vol 16 (1) ◽  
pp. 21
Author(s):  
Handoyo Handoyo ◽  
Fatkhan Fatkhan ◽  
Fourier D. E. Latief ◽  
Harnanti Y. Putri

Modern technique to estimate of the physical properties of rocks can be done by means of digital imagingand numerical simulation, an approach known as digital rock physics (DRP: Digital Rock Physics). Digital rockphysics modeling is useful to understand microstructural parameters of rocks (pores and rock matrks), quite quickly and in detail. In this paper a study was conducted on sandstone reservoir samples in a rock formation. The core of sandstone samples were calculated porosity, permeability, and elasticity parameters in the laboratory. Then performed digital image processing using CT-Scan that utilizes X-ray tomography. The result of digital image is processed and done by calculation of digital simulation to calculate porosity, permeability, and elastic parameter of sandstones. In addition, there are also predictions of p-wave velocity and wave -S using the empirical equations given by Han (1986), Raymer (1990), and Nur (1998). The results of digital simulation (DRP) in this study provide a higher than the calculations in the laboratory. The digital rock physicsmethod (DRP) combined with rock physics modeling can be a practical and rapid method for determining the rock properties of tiny (microscopic) rock fragments


Geophysics ◽  
2010 ◽  
Vol 75 (5) ◽  
pp. 75A3-75A13 ◽  
Author(s):  
Douglas J. Foster ◽  
Robert G. Keys ◽  
F. David Lane

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.


Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. C19-C29 ◽  
Author(s):  
Brian H. Russell ◽  
David Gray ◽  
Daniel P. Hampson

The technique of amplitude variation with offset (AVO) allows geoscientists to extract fluid and lithology information from the analysis of prestack seismic amplitudes. Various AVO parameterizations exist, all of which involve the sum of three weighted elastic-constant terms. In present-day AVO approaches, the weighting terms involve either knowledge of the incidence angle only, or knowledge of both the incidence angle and the in situ VP/VS ratio. We have used the theory of poroelasticity to derive a generalized AVO approximation that provides the estimation of fluid, rigidity, and density parameters. We have combined two previously independent AVO formulations, thus reducing, instead of adding to, the total number of formulations. This new approach requires knowledge of a third parameter to compute the weights: the dry-rock VP/VS ratio. We have derived a new equation and applied it to model and real data sets. The new formulation has allowed us to estimate fluid properties of the reservoir in a more direct manner than previous formulations.


Geophysics ◽  
2000 ◽  
Vol 65 (3) ◽  
pp. 755-765 ◽  
Author(s):  
Xinhua Sun ◽  
Xiaoming Tang ◽  
C. H. (Arthur) Cheng ◽  
L. Neil Frazer

In this paper, a modification of an existing method for estimating relative P-wave attenuation is proposed. By generating synthetic waveforms without attenuation, the variation of geometrical spreading related to changes in formation properties with depth can be accounted for. With the modified method, reliable P- and S-wave attenuation logs can be extracted from monopole array acoustic waveform log data. Synthetic tests show that the P- and S-wave attenuation values estimated from synthetic waveforms agree well with their respective model values. In‐situ P- and S-wave attenuation profiles provide valuable information about reservoir rock properties. Field data processing results show that this method gives robust estimates of intrinsic attenuation. The attenuation profiles calculated independently from each waveform of an eight‐receiver array are consistent with one another. In fast formations where S-wave velocity exceeds the borehole fluid velocity, both P-wave attenuation ([Formula: see text]) and S-wave attenuation ([Formula: see text]) profiles can be obtained. P- and S-wave attenuation profiles and their comparisons are presented for three reservoirs. Their correlations with formation lithology, permeability, and fractures are also presented.


Geophysics ◽  
2008 ◽  
Vol 73 (1) ◽  
pp. E1-E5 ◽  
Author(s):  
Lev Vernik

Seismic reservoir characterization and pore-pressure prediction projects rely heavily on the accuracy and consistency of sonic logs. Sonic data acquisition in wells with large relative dip is known to suffer from anisotropic effects related to microanisotropy of shales and thin-bed laminations of sand, silt, and shale. Nonetheless, if anisotropy parameters can be related to shale content [Formula: see text] in siliciclastic rocks, then I show that it is straightforward to compute the anisotropy correction to both compressional and shear logs using [Formula: see text] and the formation relative dip angle. The resulting rotated P-wave sonic logs can be used to enhance time-depth ties, velocity to effective stress transforms, and low-frequency models necessary for prestack seismic amplitude variation with offset (AVO) inversion.


Sign in / Sign up

Export Citation Format

Share Document