Figure 1. P wave speed, Vp, S wave speed, Vs, and density, ρ, for the different layers of Ostrander s gas sand model shown in SI units.
|
|
- Camron Jordan
- 5 years ago
- Views:
Transcription
1 Ambiguity in Resolving the Elastic Parameters of Gas Sands from Wide-Angle AVO Andrew J. Calvert - Simon Fraser University and David F. Aldridge - Sandia National Laboratories Summary We investigate the range of elastic parameters that will reproduce the amplitude of P waves reflected at angles of incidence up to 5 o from the top-sand and base-sand interfaces of Ostrander s gas sand model. In the case of the top-sand reflector, which lies at the top of a P wave low velocity layer, it is not possible to resolve either density or S wave speed information given P-to-P amplitude data alone. The density contrast is, however, well resolved for the base-sand reflector, but the S wave information remains ambiguous. Including converted P-to-S reflections in the analysis significantly reduces the range of elastic parameters that can reproduce the reflection amplitudes predicted by Ostrander s model. We suggest that differences in the ability to resolve density are related to the existence and location of the P wave critical point, which permits the determination of the P wave speed contrast across an interface. 1. Introduction Following the early work of Ostrander (1984), who demonstrated that P wave reflections from gas sands could exhibit an anomalous amplitude variation with increasing source-receiver offset (AVO), i.e. angle of incidence, considerable effort has been devoted to the extraction of subsurface elastic parameters from the amplitudes of reflected seismic waves. The majority of this work has focussed on the use of P-to-P (P-P) reflections because this wave type is most commonly utilised in conventional exploration surveys. For two homogeneous isotropic elastic media in welded contact, the amplitude of the P-P reflection can be determined at any angle of incidence from the Zoeppritz equations if the six parameters that completely characterise the elastic properties of the two media, for example the densities and wave speeds, are known. A number of linearised approximations to the expression for P wave reflection amplitude, which offer some insight into the factors that control its variation with changing angle of incidence, have been developed (e.g. Aki and Richards, 1980; Shuey, 1985; Smith and Gidlow, 1987; Fatti et al., 1994). These approximations typically break down at angles of incidence greater than 0 degrees, and are only valid for small changes in the elastic parameters across the interface. Nevertheless, the relatively simple form of these expression has facilitated the estimation of parameters, such as P and S reflectivity, from seismic field data in a computationally efficient fashion using least squares methods. The influence of the density contrast across an interface on the P wave reflection amplitude only becomes marked at large angles of incidence, making it difficult to extract density information with approximations that are restricted to small angles. The extraction of density information and the desire to constrain better the S wave speeds has motivated recent interest in the AVO inversion of wide-angle data (Roberts, 000). In this paper, we investigate the range of elastic parameters that will reproduce the amplitude of P waves reflected at wide-angle from the topsand and base-sand interfaces of Ostrander s Class III gas sand model (Ostrander, 1984). We disregard multiples and converted waves that would, in reality, be generated by an embedded layer. For this study, we consider the two interfaces separately and for each interface calculate the amplitude of the P wave reflection over a wide range of angles of incidence using the exact Zoeppritz equations. We then calculate P wave reflection amplitudes for a broad range of material parameters, and determine a root mean square (RMS) misfit function that we use to characterise those parameters that will reproduce the variation of reflection amplitude with angle of incidence for both interfaces of the gas sand model. In an analysis based on a linearised solution to the Zoeppritz equations, Lines (1999) showed that density is poorly resolved by P wave reflection amplitudes recorded over a limited range of angles of incidence from the top of a low impedance contrast gas sand. Our approach allows us to determine the range of elastic parameters that can be resolved at wide angle when the exact Zoeppritz equations are employed. This range of material parameters can be substantially reduced by incorporating additional information from P-to-S (P-S) converted reflections, but we infer that a full 4-parameter AVO inversion may only be possible when the P wave critical angle occurs close to the range of angles over which amplitudes are available. Figure 1. P wave speed, Vp, S wave speed, Vs, and density, ρ, for the different layers of Ostrander s gas sand model shown in SI units. 1
2 . Methodology For both interfaces of the Ostrander model, synthetic P-to-P reflection amplitude data were calculated using the exact Zoeppritz equations as a function of angle of incidence every o from 0 o to 5 o, which is just less than the P wave critical angle of the lower interface. No noise was added to these amplitude data. The values of these synthetic amplitude data thus indicate the reflection coefficient at a given angle of incidence. Using the parameterisation of the Zoeppritz equations discussed below, amplitude values were calculated for a wide range of trial parameter values, and a least-squares-based objective, i.e. misfit, function was calculated. The objective function is the logarithm to the base 10 of the root mean square average of the differences between the synthetic and trial amplitude data; thus, for example, a RMS misfit of 0.01 corresponds to a value of the objective function of. This formulation of the objective function was chosen primarily to facilitate display of the function s values. Low amplitude values of the objective function indicate the ranges of material parameters in model space that are able to reproduce well the original synthetic amplitude data. We have arbitrarily chosen a value of the objective function of, i.e. a RMS misfit of 0.01, as the cutoff point that indicates the amplitude data have been well fit. A misfit as low as this would only be possible in high quality field data but is not completely unrealistic..1 Parameterisation of the Zoeppritz Equations P wave reflection amplitudes can be computed from the P and S wave speeds plus the densities of the two media in contact, i.e. six parameters. However, the material properties can be specified in the Zoeppritz equations in terms of only four independent variables, which describe the ratios of various elastic properties (Koefoed, 196), for example V P/V P1, V S1/V P1, V S/V P1, and ρ /ρ 1, where V P, V S and ρ refer to the P wave speed, S wave speed, and mass density of the two media, and 1 and refer to the upper and lower layers respectively. Consequently doubling the densities and wave speeds of two layers in contact does not change the amplitude of the reflected, or indeed transmitted, waves. As a result, any approach to invert the reflection amplitudes can only recover these four ratios, or some function thereof, rather than individual parameters such as density or S wave speed, and any background wave speed cannot be resolved from amplitude data alone. The above parameter ratios can be rearranged to provide several other 4-element parameter sets. For our study we have chosen to calculate the objective function as a function of (i) the normal incidence P wave reflection coefficient, Rp; (ii) the Vs/Vp ratio in the upper layer; (iii) the Vs/Vp ratio in the lower layer; (iv) a density ratio, ρ /ρ 1, which is the ratio of the density of the lower layer to the density of the upper layer. As will be seen later this parameter choice has the merit in some situations of restricting ambiguity in material properties to the hyper-plane in model space defined by the two Vs/Vp ratios. Furthermore, what is strictly speaking a four-parameter problem can be reduced to three parameters by assuming that the value of Rp is known. This is a reasonable assumption since, the zero-offset reflection coefficient is the parameter that is best constrained by linearised AVO inversion methods. We are, of course, assuming that the seismic data have been well calibrated to indicate reflection coefficient. Data calibration is not an issue that we address here, but we note that this can, for example, be accomplished in deep-water marine data using the seafloor multiple sequence.. Restrictions on the Values of Elastic Parameters In order to sample the model space as completely as possible, it is necessary to determine the range of parameter values that are consistent with the theory of linear elasticity. In the absence of a vacuum, mass density is always positive. The strain energy of a linear elastic solid is also positive. For an isotropic solid, positive strain energy implies that the shear modulus, µ, and the bulk modulus, κ, which is also known as the incompressibility, are positive (e.g. Debski and Tarantola, 1995). Since the incompressibility is a function of the first Lamé parameter, λ: κ = λ + µ, we must have λ µ >. The Lamé parameter, λ can thus be negative, and these negative values correspond to materials with negative Poisson s ratios such as α- cristobalite (Yeganeh-Haeri et al., 199). Poisson s ratio must lie between 1 and These restrictions in turn limit the ratio of the S wave speed to the P wave speed: VS 0 < V P These inequalities define the regions of a hyper-plane in model space, in which real-world materials can lie. In determining values of the objective function, the Vs/Vp ratios were varied between the above limits, and the density ratio was allowed to vary between 0.5 and Geometry of the Objective Function A measure of the ambiguity of the solution to the nonlinear inversion of these P wave reflection amplitudes for the material parameters can be obtained by calculating hyper-planes passing through the point in model space that represents the parameters used to calculate the original synthetic data. We show two key hyper-planes for the top-sand and base-sand reflectors in Figure. The density and P wave speed ratios appear to be relatively well resolved against one another (Figure a, c); the parameters used to calculate the synthetic amplitudes are Rp = , ρ /ρ 1 = 0.8 for the top-sand reflector and Rp = , ρ /ρ 1 = 1.5 for the base-sand reflector. In contrast, the Vs/Vp ratios are poorly resolved against one another. The Vs/Vp ratios of the shale and gas sand layers are and respectively. In the case of the base-sand reflector, a single, but elongated, minimum of the objective function is apparent. For the top sand reflector, multiple minima exist in the objective function along a trajectory that extends across much of the range of physically realisable values; non-physical parameter combinations, i.e. negative S wave speed in the lower layer satisfy the amplitude data, but are excluded from Figure b by the bounds imposed on the parameters used to calculate the objective function, as discussed previously..
3 4. Elastic Parameters Derived from P-P Amplitudes In order to characterise the range of elastic parameters that fit the synthetic amplitude data, we have fixed the value of Rp to be the actual normal-incidence reflection coefficient, and then calculated values of the objective function. As described previously, we have arbitrarily chosen a RMS misfit of less than 0.01, i.e. values of the objective function less than, to be the set of parameter values that reproduce well the synthetic amplitude data. We have then plotted this volume in the three-dimensional space defined by the density ratio and the two Vs/Vp ratios. The set of points in model space that are able to reproduce the P wave AVO of the base-sand interface lies parallel to the plane defined by the two Vs/Vp ratios, and forms a tube restricted to a very narrow range of density values. Thus the density contrast across the base-sand reflector can be well resolved by nonlinear inversion using the exact Zoeppritz equation for the P wave reflection amplitude. The Vs/Vp ratios are, however, poorly resolved, and a broad range of values is possible. In the case of the top-sand reflector, the volume of points that reproduce the synthetic amplitudes, which approximates a sheet-like structure, extends across almost the entire model space used to calculate values of the objective function. Part of the volume of admissible points occurs at values where the Vs/Vp ratio in the lower, i.e. gas sand, layer are greater than 0.75, indicating a material with negative Poisson s ratio. Clearly for the top-sand interface, the density and the two Vs/Vp ratios are all poorly resolved. We ascribe these significant differences between the two interfaces in the regions of model space that satisfy the synthetic amplitude data to the existence of a P wave critical point, which allows the density contrast to be determined, in the case of the base-sand reflector. At a decrease in P wave speed, there is no critical point and the density cannot be resolved. Figure. Values of the objective function, i.e. log 10 of the RMS amplitude error, along two hyper-planes through the point in model space that represents the parameters used to calculate the synthetic amplitude data. a) Rp-ρ /ρ 1 for top sand, b) Vs/Vp upper layer Vs/Vp lower layer for top sand, c) Rp-ρ /ρ 1 for base-sand, d) Vs/Vp upper layer Vs/Vp lower layer for base-sand.
4 Figure. Representation of the regions of model space that can reproduce well the P-P amplitude data, for a correctly specified Rp. a) Top-sand reflector, b) Base-sand reflector. 5. Elastic Parameters Derived from P-P and P-S Amplitudes To reduce the size of the set of parameter values that will reproduce the synthetic data, and within which a nonlinear inversion scheme must identify a solution, we have modified the objective function to include both P-P and converted P-S reflections. Repeating the preceding analysis, produces the displays shown in Figure 4. For the base-sand interface, Figure 4b indicates that the tradeoff between the two Vs/Vp ratios is eliminated, and that a global minimum can be isolated by an AVO inversion. The ambiguity in solutions for the top-sand interface is greatly reduced by inclusion of P-S reflections, but a broad range of material properties is still possible, and other a priori constraints are required to identify a meaningful solution. Figure 4. Representation of the regions of model space that can reproduce well both the P-P and P-S amplitude data, for a correctly specified Rp. a) Top-sand reflector, b) Base-sand reflector. 4
5 Acknowledgements Jon Downton provided helpful insights insight into linearised inversion methods. This project was funded through the Natural Sciences and Engineering Council of Canada with support from Paradigm Geophysical and Landmark Graphics. References Aki, K., and Richards, P.G., 1980, Quantitative seismology: W.H. Freeman and Co. Debski, W., and Tarantola, A., 1995, Information on elastic parameters obtained from the amplitudes of reflected waves: Geophysics, 60, Fatti, J.L., Smith, G.C., Vail, P.J., Strauss, P.J., and Levitt, P.R., 1994, Detection of gas in sandstone reservoirs using AVO analysis: A -D seismic case history using the Geostack technique: Geophysics, 59, Koefoed, O., 196, Reflection and transmission coefficients for plane longitudinal incident waves: Geophys. Prosp., 10, Lines, L., 1999, Density and AVO, Can. J. Expl. Geophys., 5, -5. Ostrander, W.J., 1984, Plane-wave reflection coefficients for gas sands at nonnormal angles of incidence: Geophysics, 49, Roberts, G., 000, Wide-angle AVO, 70th Ann. Internat. Mtg. Soc. Expl. Geophys., Expanded Abstracts, abstract AVO.. Shuey, R.T., 1985, A simplification of the Zoeppritz equations: Geophysics, 51, Smith, G.C., and Gidlow, P.M., 1987, Weighted stacking for rock property estimation and detection of gas: Geophys. Prosp., 5, Yeganeh-Haeri, A., Weidner, D.J., and Parise, J.B., 199, Elasticity of α-cristobalite: A silicon dioxide with negative Poisson s ratio: Science, 57,
Modelling of linearized Zoeppritz approximations
Modelling of linearized Zoeppritz approximations Arnim B. Haase Zoeppritz approximations ABSTRACT The Aki and Richards approximations to Zoeppritz s equations as well as approximations by Stewart, Smith
More informationJoint PP and PS AVO inversion based on Zoeppritz equations
Earthq Sci (2011)24: 329 334 329 doi:10.1007/s11589-011-0795-1 Joint PP and PS AVO inversion based on Zoeppritz equations Xiucheng Wei 1,2, and Tiansheng Chen 1,2 1 SINOPEC Key Laboratory of Seismic Multi-Component
More informationAn overview of AVO and inversion
P-486 An overview of AVO and inversion Brian Russell, Hampson-Russell, CGGVeritas Company Summary The Amplitude Variations with Offset (AVO) technique has grown to include a multitude of sub-techniques,
More informationUseful approximations for converted-wave AVO
GEOPHYSICS VOL. 66 NO. 6 NOVEMBER-DECEMBER 001); P. 171 1734 14 FIGS. 3 TABLES. Useful approximations for converted-wave AVO Antonio C. B. Ramos and John P. Castagna ABSTRACT Converted-wave amplitude versus
More information= (G T G) 1 G T d. m L2
The importance of the Vp/Vs ratio in determining the error propagation and the resolution in linear AVA inversion M. Aleardi, A. Mazzotti Earth Sciences Department, University of Pisa, Italy Introduction.
More informationLinearized AVO and Poroelasticity for HRS9. Brian Russell, Dan Hampson and David Gray 2011
Linearized AO and oroelasticity for HR9 Brian Russell, Dan Hampson and David Gray 0 Introduction In this talk, we combine the linearized Amplitude ariations with Offset (AO) technique with the Biot-Gassmann
More informationA simple way to improve AVO approximations
A simple way to improve AVO approximations Charles P. Ursenbach A simple way to improve AVO approximations ABSTRACT Some twenty years ago it was suggested that the average angle, = ( + )/, in the Aki-Richards
More informationConstrained inversion of P-S seismic data
PS Inversion Constrained inversion of P-S seismic data Robert J. Ferguson, and Robert R. Stewart ABSTRACT A method to estimate S-wave interval velocity, using P-S seismic data is presented. The method
More informationSimultaneous Inversion of Pre-Stack Seismic Data
6 th International Conference & Exposition on Petroleum Geophysics Kolkata 006 Summary Simultaneous Inversion of Pre-Stack Seismic Data Brian H. Russell, Daniel P. Hampson, Brad Bankhead Hampson-Russell
More informationAzimuthal AVO and Curvature. Jesse Kolb* David Cho Kris Innanen
Azimuthal AVO and Curvature Jesse Kolb* David Cho Kris Innanen Azimuthal AVO What is azimuthal AVO?: Analysis of incidence angle and azimuthal amplitude variations of reflection coefficients; Measures
More informationA Case Study on Simulation of Seismic Reflections for 4C Ocean Bottom Seismometer Data in Anisotropic Media Using Gas Hydrate Model
A Case Study on Simulation of Seismic Reflections for 4C Ocean Bottom Seismometer Data in Anisotropic Media Using Gas Hydrate Model Summary P. Prasada Rao*, N. K. Thakur 1, Sanjeev Rajput 2 National Geophysical
More informationReservoir Characterization using AVO and Seismic Inversion Techniques
P-205 Reservoir Characterization using AVO and Summary *Abhinav Kumar Dubey, IIT Kharagpur Reservoir characterization is one of the most important components of seismic data interpretation. Conventional
More informationIntegrating reservoir flow simulation with time-lapse seismic inversion in a heavy oil case study
Integrating reservoir flow simulation with time-lapse seismic inversion in a heavy oil case study Naimeh Riazi*, Larry Lines*, and Brian Russell** Department of Geoscience, University of Calgary **Hampson-Russell
More informationP125 AVO for Pre-Resonant and Resonant Frequency Ranges of a Periodical Thin-Layered Stack
P125 AVO for Pre-Resonant and Resonant Frequency Ranges of a Periodical Thin-Layered Stack N. Marmalyevskyy* (Ukrainian State Geological Prospecting Institute), Y. Roganov (Ukrainian State Geological Prospecting
More informationSpherical wave AVO-modelling in elastic isotropic media
Spherical wave AVO-modelling in elastic isotropic media Arnim B. Haase and Charles P. Ursenbach ABSTRACT Spherical wave AVO The AVO-response of two-layer isotropic models for AVO-Classes 1, 2, 3 and 4
More informationP-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry
GEOPHYSICS, VOL. 62, NO. 3 (MAY-JUNE 1997); P. 713 722, 7 FIGS., 1 TABLE. P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry Andreas Rüger ABSTRACT
More informationCOMPARISON OF OPTICAL AND ELASTIC BREWSTER S ANGLES TO PROVIDE INVUITIVE INSIGHT INTO PROPAGATION OF P- AND S-WAVES. Robert H.
COMPARISON OF OPTICAL AND ELASTIC BREWSTER S ANGLES TO PROVIDE INVUITIVE INSIGHT INTO PROPAGATION OF P- AND S-WAVES Robert H. Tatham Department of Geological Sciences The University of Texas at Austin
More informationFramework for AVO gradient and intercept interpretation
GEOPHYSICS, VOL. 63, NO. 3 (MAY-JUNE 1998); P. 948 956, 10 FIGS., 2 TABLES. Framework for AVO gradient and intercept interpretation John P. Castagna, Herbert W. Swan, and Douglas J. Foster ABSTRACT Amplitude
More informationGEOPHYSICAL PROSPECTING: DYNAMIC RESERVOIR CHARACTERIZATION AND TIME-LAPSE MULTICOMPONENT SEISMOLOGY FOR RESERVOIR MONITORING UNESCO EOLSS
GEOPHYSICAL PROSPECTING: DYNAMIC RESERVOIR CHARACTERIZATION AND TIME-LAPSE MULTICOMPONENT SEISMOLOGY FOR RESERVOIR MONITORING Steven L. Roche CGGVeritas, Multicomponent Processing & Technology Group Thomas
More informationReservoir connectivity uncertainty from stochastic seismic inversion Rémi Moyen* and Philippe M. Doyen (CGGVeritas)
Rémi Moyen* and Philippe M. Doyen (CGGVeritas) Summary Static reservoir connectivity analysis is sometimes based on 3D facies or geobody models defined by combining well data and inverted seismic impedances.
More informationA New AVO Attribute for Hydrocarbon Prediction and Application to the Marmousi II Dataset*
A New AVO Attribute for Hydrocarbon Prediction and Application to the Marmousi II Dataset* Changcheng Liu 1 and Prasad Ghosh 2 Search and Discovery Article #41764 (2016) Posted January 25, 2016 *Adapted
More informationSEG/San Antonio 2007 Annual Meeting. Summary
A comparison of porosity estimates obtained using post-, partial-, and prestack seismic inversion methods: Marco Polo Field, Gulf of Mexico. G. Russell Young* and Mrinal K. Sen, The Institute for Geophysics
More information23855 Rock Physics Constraints on Seismic Inversion
23855 Rock Physics Constraints on Seismic Inversion M. Sams* (Ikon Science Ltd) & D. Saussus (Ikon Science) SUMMARY Seismic data are bandlimited, offset limited and noisy. Consequently interpretation of
More informationAn empirical method for estimation of anisotropic parameters in clastic rocks
An empirical method for estimation of anisotropic parameters in clastic rocks YONGYI LI, Paradigm Geophysical, Calgary, Alberta, Canada Clastic sediments, particularly shale, exhibit transverse isotropic
More informationP235 Modelling Anisotropy for Improved Velocities, Synthetics and Well Ties
P235 Modelling Anisotropy for Improved Velocities, Synthetics and Well Ties P.W. Wild* (Ikon Science Ltd), M. Kemper (Ikon Science Ltd), L. Lu (Ikon Science Ltd) & C.D. MacBeth (Heriot Watt University)
More informationAn empirical study of hydrocarbon indicators
An empirical study of hydrocarbon indicators Brian Russell 1, Hong Feng, and John Bancroft An empirical study of hydrocarbon indicators 1 Hampson-Russell, A CGGVeritas Company, Calgary, Alberta, brian.russell@cggveritas.com
More informationEdinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains
Frequency-dependent AVO attribute: theory and example Xiaoyang Wu, 1* Mark Chapman 1,2 and Xiang-Yang Li 1 1 Edinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains Road, Edinburgh
More informationElastic impedance inversion from robust regression method
Elastic impedance inversion from robust regression method Charles Prisca Samba 1, Liu Jiangping 1 1 Institute of Geophysics and Geomatics,China University of Geosciences, Wuhan, 430074 Hubei, PR China
More informationApproximations to the Zoeppritz equations and their use in AVO analysis
GEOPHYSICS, VOL. 64, NO. 6 (NOVEMBER-DECEMBER 1999; P. 190 197, 4 FIGS., 1 TABLE. Approximations to the Zoeppritz equations their use in AVO analysis Yanghua Wang ABSTRACT To efficiently invert seismic
More informationRC 1.3. SEG/Houston 2005 Annual Meeting 1307
from seismic AVO Xin-Gong Li,University of Houston and IntSeis Inc, De-Hua Han, and Jiajin Liu, University of Houston Donn McGuire, Anadarko Petroleum Corp Summary A new inversion method is tested to directly
More informationEffects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses
PROC. ITB Eng. Science Vol. 38 B, No. 2, 2006, 159-170 159 Effects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses Fatkhan Program Studi Teknik Geofisika FIKTM-ITB
More informationComparative Study of AVO attributes for Reservoir Facies Discrimination and Porosity Prediction
5th Conference & Exposition on Petroleum Geophysics, Hyderabad-004, India PP 498-50 Comparative Study of AVO attributes for Reservoir Facies Discrimination and Porosity Prediction Y. Hanumantha Rao & A.K.
More informationIntroduction. Theory. GEOHORIZONS December 2007/22. Summary
Seismic amplitude versus angle modeling of a bottom-simulating reflector Maheswar Ojha and Kalachand Sain National Geophysical Research Institute, Uppal Road, Hyderabad - 500 007, India * Corresponding
More informationComments on the Geophysics paper Multiparameter 11norm. waveform fitting: Interpretation of Gulf of Mexico reflection
,..-., Comments on the Geophysics paper Multiparameter 11norm waveform fitting: Interpretation of Gulf of Mexico reflection seismograms by H. Djikp&s6 and A. Tarantola Susan E. Minkoff * t and William
More informationP-wave impedance, S-wave impedance and density from linear AVO inversion: Application to VSP data from Alberta
Linear AVO inversion P-wave impedance, S-wave impedance and density from linear AVO inversion: Application to VSP data from Alberta Faranak Mahmoudian and Gary F. Margrave ABSTRACT In AVO (amplitude variation
More informationSeismic reservoir characterization in offshore Nile Delta.
Seismic reservoir characterization in offshore Nile Delta. Part II: Probabilistic petrophysical-seismic inversion M. Aleardi 1, F. Ciabarri 2, B. Garcea 2, A. Mazzotti 1 1 Earth Sciences Department, University
More informationUnphysical negative values of the anelastic SH plane wave energybased transmission coefficient
Shahin Moradi and Edward S. Krebes Anelastic energy-based transmission coefficient Unphysical negative values of the anelastic SH plane wave energybased transmission coefficient ABSTRACT Computing reflection
More informationAVO modeling and the locally converted shear wave
GEOPHYSICS, VOL. 59, NO. 9 (AUGUST 1994); P. 1237-1248, 11 FIGS., 4 TABLES. AVO modeling and the locally converted shear wave James L. Simmons, Jr.* and Milo M. Backus* ABSTRACT The locally converted shear
More informationAVO and AVA inversion challenges: a conceptual overview
AVO and AVA inversion challenges: a conceptual overview Jeff P. Grossman ABSTRACT Inversion of seismic data for earth parameters involves two main steps: (1) estimate the reflectivity as a function of
More informationInversion of seismic AVA data for porosity and saturation
Inversion of seismic AVA data for porosity and saturation Brikt Apeland Thesis for the degree Master of Science Department of Earth Science University of Bergen 27th June 213 2 Abstract Rock physics serves
More informationSensitivity Analysis of Pre stack Seismic Inversion on Facies Classification using Statistical Rock Physics
Sensitivity Analysis of Pre stack Seismic Inversion on Facies Classification using Statistical Rock Physics Peipei Li 1 and Tapan Mukerji 1,2 1 Department of Energy Resources Engineering 2 Department of
More informationRelationships among elastic-wave values (R pp, R ps, R ss, Vp, Vs, ρ, σ, κ)
Relationships among elastic-wave values (R pp, R ps, R ss, Vp, Vs, ρ, σ, κ) Robert R. Stewart, Qi Zhang, Felix Guthoff * ABSTRACT The average Vp/Vs value of a set of layers is a weighted sum of the interval
More informationERTH2020 Introduction to Geophysics The Seismic Method. 1. Basic Concepts in Seismology. 1.1 Seismic Wave Types
ERTH2020 Introduction to Geophysics The Seismic Method 1. Basic Concepts in Seismology 1.1 Seismic Wave Types Existence of different wave types The existence of different seismic wave types can be understood
More informationThe signature of attenuation and anisotropy on AVO and inversion sensitivities
The signature of attenuation and anisotropy on AVO and inversion sensitivities Shahpoor Moradi and Kristopher Innanen University of Calgary, Department of Geoscience, Calgary, Canada Summary We study the
More informationFluid-property discrimination with AVO: A Biot-Gassmann perspective
Fluid-property discrimination with AVO: A Biot-Gassmann perspective Brian H. Russell, Ken Hedlin 1, Fred J. Hilterman, and Laurence R. Lines ABSTRACT This paper draws together basic rock physics, AVO,
More informationOptimal Zoeppritz approximations
harles. Ursenbach ABSTRAT New approximations are developed for R and R S that are optimal in the sense that they preserve as much accuracy as possible with as much simplicity as possible. In particular,
More informationRock physics integration of CSEM and seismic data: a case study based on the Luva gas field.
Rock physics integration of CSEM and seismic data: a case study based on the Luva gas field. Peter Harris*, Zhijun Du, Harald H. Soleng, Lucy M. MacGregor, Wiebke Olsen, OHM-Rock Solid Images Summary It
More informationStochastic vs Deterministic Pre-stack Inversion Methods. Brian Russell
Stochastic vs Deterministic Pre-stack Inversion Methods Brian Russell Introduction Seismic reservoir analysis techniques utilize the fact that seismic amplitudes contain information about the geological
More informationAN AVO METHOD TOWARD DIRECT DETECTION OF LITHOLOGIES COMBINING P-P AND P-S REFLECTION DATA. A Thesis JUAN RAMON DE JESUS CARCUZ JEREZ
AN AVO METHOD TOWARD DIRECT DETECTION OF LITHOLOGIES COMBINING P-P AND P-S REFLECTION DATA A Thesis by JUAN RAMON DE JESUS CARCUZ JEREZ Submitted to the Office of Graduate Studies of Texas A&M University
More informationSENSITIVITY ANALYSIS OF AMPLITUDE VARIATION WITH OFFSET (AVO) IN FRACTURED MEDIA
SENSITIVITY ANALYSIS OF AMPLITUDE VARIATION WITH OFFSET AVO) IN FRACTURED MEDIA Mary L. Krasovec, William L. Rodi, and M. Nafi Toksoz Earth Resources Laboratory Department of Earth, Atmospheric, and Planetary
More informationMultiple Scenario Inversion of Reflection Seismic Prestack Data
Downloaded from orbit.dtu.dk on: Nov 28, 2018 Multiple Scenario Inversion of Reflection Seismic Prestack Data Hansen, Thomas Mejer; Cordua, Knud Skou; Mosegaard, Klaus Publication date: 2013 Document Version
More informationAVAZ and VVAZ practical analysis to estimate anisotropic properties
AVAZ and VVAZ practical analysis to estimate anisotropic properties Yexin Liu*, SoftMirrors Ltd., Calgary, Alberta, Canada yexinliu@softmirrors.com Summary Seismic anisotropic properties, such as orientation
More informationPre-Stack Seismic Inversion and Amplitude Versus Angle Modeling Reduces the Risk in Hydrocarbon Prospect Evaluation
Advances in Petroleum Exploration and Development Vol. 7, No. 2, 2014, pp. 30-39 DOI:10.3968/5170 ISSN 1925-542X [Print] ISSN 1925-5438 [Online] www.cscanada.net www.cscanada.org Pre-Stack Seismic Inversion
More informationAmplitude variation with offset AVO. and. Direct Hydrocarbon Indicators DHI. Reflection at vertical incidence. Reflection at oblique incidence
Amplitude variation with offset AVO and Direct Hydrocarbon Indicators DHI Reflection at vertical incidence Reflection coefficient R(p) c α 1 S wavespeed β 1 density ρ 1 α 2 S wavespeed β 2 density ρ 2
More informationAmplitude Variation with Offset: Case History in Water Table and Lithological Reflection
IOSR Journal of Applied Geology and Geophysics (IOSR-JAGG) e-issn: 31 0990, p-issn: 31 098.Volume, Issue 5 (Jul-Aug. 014), PP 51-58 Amplitude Variation with Offset: Case History in Water Table and Lithological
More informationA simplification of the Zoeppritz equations. R. T. Shuey* (a) When the underlying medium has the greater longitudinal
GEOPHYSICS. VOL. 50, NO, 4 (APRIL 1985); P. 609-614, 8 FIGS., 7, TABLES. A simplification of the Zoeppritz equations R. T. Shuey* ABSTRACT The compressional wave reflection coefficient R(B) given by the
More informationPitfall in AVO Anisotropic Modeling: Plane Wave vs Spherical Wave. Qing Li May, 2003
Pitfall in AVO Anisotropic Modeling: Plane Wave vs Spherical Wave Qing Li May, 2003 A Client s Model A client uses the following model to model anisotropic AVO effects and raised the question of discrepancy
More informationPractical aspects of AVO modeling
Practical aspects of AVO modeling YONGYI LI, Paradigm Geophysical, Calgary, Canada JONATHAN DOWNTON, Veritas, Calgary, Canada, YONG XU, Arcis Corporation, Calgary, Canada AVO (amplitude variation with
More informationLinearized AVO in viscoelastic media Shahpoor Moradi,Kristopher A. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada
Shahpoor Moradi,Kristopher. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada SUMMRY Study of linearized reflectivity is very important for amplitude versus offset VO) analysis.
More informationEstimation of density from seismic data without long offsets a novel approach.
Estimation of density from seismic data without long offsets a novel approach. Ritesh Kumar Sharma* and Satinder Chopra Arcis seismic solutions, TGS, Calgary Summary Estimation of density plays an important
More informationA simple way to improve PP and PS AVO approximations. Chuck Ursenbach CREWES Sponsors Meeting Thursday, December 1, 2005
A simple way to improve PP and PS AVO approximations Chuck Ursenbach CREWES Sponsors Meeting Thursday, December 1, 2005 Overview Notes on spherical-wave modeling Reflectivity Explorer observations Theoretical
More informationDownloaded 11/20/12 to Redistribution subject to SEG license or copyright; see Terms of Use at
AVO crossplot analysis in unconsolidated sediments containing gas hydrate and free gas: Green Canyon 955, Gulf of Mexico Zijian Zhang* 1, Daniel R. McConnell 1, De-hua Han 2 1 Fugro GeoConsulting, Inc.,
More informationRecent advances in application of AVO to carbonate reservoirs: case histories
Recent advances in application of AVO to reservoirs: case histories Yongyi Li, Bill Goodway*, and Jonathan Downton Core Lab Reservoir Technologies Division *EnCana Corporation Summary The application of
More informationAVO responses for varying Gas saturation sands A Possible Pathway in Reducing Exploration Risk
P 420 AVO responses for varying Gas saturation sands A Possible Pathway in Reducing Exploration Risk Naina Gupta*, Debajyoti Guha, Ashok Yadav, Santan Kumar, Ashutosh Garg *Reliance Industries Limited,
More informationEstimation of Elastic Parameters Using Two-Term Fatti Elastic Impedance Inversion
Journal of Earth Science, Vol. 6, No. 4, p. 556 566, August 15 ISSN 1674-487X Printed in China DOI:.7/s158-15-564-5 Estimation of Elastic Parameters Using Two-Term Fatti Elastic Impedance Inversion Jin
More informationWe apply a rock physics analysis to well log data from the North-East Gulf of Mexico
Rock Physics for Fluid and Porosity Mapping in NE GoM JACK DVORKIN, Stanford University and Rock Solid Images TIM FASNACHT, Anadarko Petroleum Corporation RICHARD UDEN, MAGGIE SMITH, NAUM DERZHI, AND JOEL
More informationSummary. Seismic Field Example
Water-saturation estimation from seismic and rock-property trends Zhengyun Zhou*, Fred J. Hilterman, Haitao Ren, Center for Applied Geosciences and Energy, University of Houston, Mritunjay Kumar, Dept.
More informationReservoir properties inversion from AVO attributes
Reservoir properties inversion from AVO attributes Xin-gang Chi* and De-hua Han, University of Houston Summary A new rock physics model based inversion method is put forward where the shaly-sand mixture
More informationNonlinear Bayesian joint inversion of seismic reflection
Geophys. J. Int. (2) 142, Nonlinear Bayesian joint inversion of seismic reflection coefficients Tor Erik Rabben 1,, Håkon Tjelmeland 2, and Bjørn Ursin 1 1 Department of Petroleum Engineering and Applied
More informationWe Challenges in shale-reservoir characterization by means of AVA and AVAZ
We-06-15 Challenges in shale-reservoir characterization by means of AVA and AVAZ N.C. Banik* (WesternGeco), M. Egan (WesternGeco), A. Koesoemadinata (WesternGeco) & A. Padhi (WesternGeco) SUMMARY In most
More informationEstimating shear-wave velocity logs using P-S seismic inversion. Robert J. Ferguson, Robert R. Stewart,and Henry C. Bland
P-S SeismicInversion Estimating shear-wave velocity logs using P-S seismic inversion. Robert J. Ferguson, Robert R. Stewart,and Henry C. Bland ABSTRACT A method is developed to estimate S-wave velocity
More informationQUANTITATIVE INTERPRETATION
QUANTITATIVE INTERPRETATION THE AIM OF QUANTITATIVE INTERPRETATION (QI) IS, THROUGH THE USE OF AMPLITUDE ANALYSIS, TO PREDICT LITHOLOGY AND FLUID CONTENT AWAY FROM THE WELL BORE This process should make
More informationThe Implementation of AVO Formulations to the Water Saturated Sandstone and Lithological Reflection in Shallow Gound Water
J. Basic. Appl. Sci. Res., 6(5)6-33, 016 016, TextRoad Publication ISSN 090-4304 Journal of Basic and Applied Scientific Research www.textroad.com The Implementation of AVO Formulations to the Water Saturated
More informationA simple algorithm for band-limited impedance inversion
Impedance inversion algorithm A simple algorithm for band-limited impedance inversion Robert J. Ferguson and Gary F. Margrave ABSTRACT This paper describes a seismic inversion method which has been cast
More informationAVO Crossplotting II: Examining Vp/Vs Behavior
AVO Crossplotting II: Examining Vp/Vs Behavior Heath Pelletier* Talisman Energy, Calgary, AB hpelletier@talisman-energy.com Introduction The development of AVO crossplot analysis has been the subject of
More informationDerived Rock Attributes Analysis for Enhanced Reservoir Fluid and Lithology Discrimination
IOSR Journal of Applied Geology and Geophysics (IOSR-JAGG) e-issn: 2321 0990, p-issn: 2321 0982.Volume 5, Issue 2 Ver. I (Mar. - Apr. 2017), PP 95-105 www.iosrjournals.org Derived Rock Attributes Analysis
More informationEFFECTS OF FRACTURED RESERVOIR ELASTIC PROPERTIES ON AZIMUTHAL AVO. Feng Shen, Xiang Zhu, and M. Nafi Toksoz
EFFECTS OF FRACTURED RESERVOIR ELASTIC PROPERTIES ON AZIMUTHAL AVO Feng Shen, Xiang Zhu, and M. Nafi Toksoz Earth Resources Laboratory Department of Earth, Atmospheric, and Planetary Sciences Massachusetts
More informationJoint inversion of AVA data for elastic parameters by bootstrapping
ARTICLE IN PRESS Computers & Geosciences 33 (2007) 367 382 www.elsevier.com/locate/cageo Joint inversion of AVA data for elastic parameters by bootstrapping Hu lya Kurt Istanbul Technical University, Department
More informationMatrix forms for the Knott-Zoeppritz equations
Kris Innanen ABSTRACT In this note we derive convenient matrix forms for the Knott-Zoeppritz equations. The attempt is to recapture results used (quoted not derived) by Levin and Keys and to extend these
More informationQuantitative Interpretation
Quantitative Interpretation The aim of quantitative interpretation (QI) is, through the use of amplitude analysis, to predict lithology and fluid content away from the well bore. This process should make
More informationElastic wave-equation migration for laterally varying isotropic and HTI media. Richard A. Bale and Gary F. Margrave
Elastic wave-equation migration for laterally varying isotropic and HTI media Richard A. Bale and Gary F. Margrave a Outline Introduction Theory Elastic wavefield extrapolation Extension to laterally heterogeneous
More informationSEISMIC INVERSION OVERVIEW
DHI CONSORTIUM SEISMIC INVERSION OVERVIEW Rocky Roden September 2011 NOTE: Terminology for inversion varies, depending on the different contractors and service providers, emphasis on certain approaches,
More informationDetermination of reservoir properties from the integration of CSEM and seismic data
Determination of reservoir properties from the integration of CSEM and seismic data Peter Harris, 1 Rock Solid Images, and Lucy MacGregor, 2 Offshore Hydrocarbons Mapping, discuss the advantages in reservoir
More informationMethane hydrate rock physics models for the Blake Outer Ridge
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 307 Methane hydrate rock physics models for the Blake Outer Ridge Christine Ecker 1 ABSTRACT Seismic analyses of methane hydrate data from
More informationPrediction of Shale Plugs between Wells in Heavy Oil Sands using Seismic Attributes
GRAY, F. DAVID *, PAUL F. ANDERSON, and JAY A. GUNDERSON *, * Veritas DGC, Calgary, AB, Canada, Apache Canada Ltd., Calgary, AB, Canada Prediction of Shale Plugs between Wells in Heavy Oil Sands using
More informationAnisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements
Anisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements K. Key (Scripps Institution of Oceanography), Z. Du* (PGS), J. Mattsson (PGS), A.
More informationIntroducción a la Geofísica
Introducción a la Geofísica 2010-01 TAREA 7 1) FoG. A plane seismic wave, travelling vertically downwards in a rock of density 2200 kg m -3 with seismic velocity 2,000 m s -1, is incident on the horizontal
More informationHydrogeophysics - Seismics
Hydrogeophysics - Seismics Matthias Zillmer EOST-ULP p. 1 Table of contents SH polarized shear waves: Seismic source Case study: porosity of an aquifer Seismic velocities for porous media: The Frenkel-Biot-Gassmann
More informationANISOTROPIC PRESTACK DEPTH MIGRATION: AN OFFSHORE AFRICA CASE STUDY
Copyright 000 by the Society of Exploration Geophysicists ANISOTROPIC PRESTACK DEPTH MIGRATION: AN OFFSHORE AFRICA CASE STUDY Philippe Berthet *, Paul Williamson *, Paul Sexton, Joachim Mispel * * Elf
More information6298 Stress induced azimuthally anisotropic reservoir - AVO modeling
6298 Stress induced azimuthally anisotropic reservoir - AVO modeling M. Brajanovski* (Curtin University of Technology), B. Gurevich (Curtin University of Technology), D. Nadri (CSIRO) & M. Urosevic (Curtin
More informationSatish Singh* (IPG Paris, France, Tim Sears (British Gas, UK), Mark Roberts (IPG Paris, Summary. Introduction P - 92
P - 92 Fine-Scale P- and S-Wave Velocities From Elastic Full Waveform Inversion of Multi-Component and Time-Lapse Data: Future of Quantitative Seismic Imaging Satish Singh* (IPG Paris, France, singh@ipgp.jussieu.fr),
More informationIntroduction: Simultaneous AVO Inversion:
Implementation of AVO, AVOAz Inversion and Ant Tracking Techniques in Wembley Valhalla Integrated Merge 3D Seismic Survey, Alberta Homayoun Gerami, Patty Evans WesternGeco Introduction: The Wembley Valhalla
More informationModelling Class 1 AVO responses of a three layer system
Modelling Class 1 AVO responses of a three layer system Arnim B. Haase ABSTRACT Class 1 three-layer model AVO-responses are computed by a method developed from the Ewing-algorithm for point sources in
More informationP191 Bayesian Linearized AVAZ Inversion in HTI Fractured Media
P9 Bayesian Linearized AAZ Inversion in HI Fractured Media L. Zhao* (University of Houston), J. Geng (ongji University), D. Han (University of Houston) & M. Nasser (Maersk Oil Houston Inc.) SUMMARY A new
More informationAVAZ inversion for fracture orientation and intensity: a physical modeling study
AVAZ inversion for fracture orientation and intensity: a physical modeling study Faranak Mahmoudian and Gary F Margrave ABSTRACT AVAZ inversion We present a pre-stack amplitude inversion of P-wave data
More informationIntegration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties
Integration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties Amaro C. 1 Abstract: The main goal of reservoir modeling and characterization is the inference of
More informationFluid Substitution and AVO Modeling for Hydrocarbon and Lithology Discrimination in X-Field, Offshore Niger Delta
Standard Scientific Research and Essays Vol2 (8): 333-344, August 2014 (ISSN: 2310-7502) http://www.standresjournals.org/journals/ssre Research Article Fluid Substitution and AVO Modeling for Hydrocarbon
More informationApplication of AVO Attributes and Fluid Replacement in Determining AVO Anomaly; Case Study of Aboozar Oil Field, Iran
Application of AVO Attributes and Fluid Replacement in Determining AVO Anomaly; Case Study of Aboozar Oil Field, Iran Sadeq Rahim Bakhtiari 1, Ehsan Farahbakhsh 2, Rasoul Hamidzadeh Moghadam 1 1- Department
More informationExact elastic impedance in orthorhombic media
Exact elastic impedance in orthorhombic media F. Zhang (hina University of Petroleum), X.Y. Li (hina University of Petroleum, British Geological Survey) SUMMARY onventional elastic/ray impedance approximations
More informationSUPPLEMENTARY INFORMATION
Inability of additional parameters to resolve the Rayleigh-Love discrepancy Radial anisotropy is introduced to resolve the Rayleigh-Love misfit discrepancy that exists across large regions of the western
More information