Constrained inversion of P-S seismic data
|
|
- Adam Lamb
- 5 years ago
- Views:
Transcription
1 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 is composed of three steps. First, the P-S data are converted to a relative change in S-velocity section. The section is then linearized. Constrained inversion is used to convert this last section to S-interval velocity. To illustrate the method, the S velocity of a field data set was derived, and a Vp/Vs section was computed. The constrained inversion was found to provide reasonable velocities and Vp/Vs. INTRODUCTION Accurate measurements of S-wave velocity (β), P-wave velocity (α) and Vp/Vs (α/β) can be used to describe the lithologic properties of rock. Such relationships can be complex, but when used in conjunction with well-log or other independent measurements, they can provide a coherent interpretation of reservoir properties (Tatham and McCormick, 1991). Until recently, seismic exploration has consisted mainly of the interpretation of P-P data, and P velocity inversion has long been part of the interpretation process (Oldenburg et al., 1983; Smith and Gidlow, 1987; Russell and Hampson, 1991). However, since the early 1980 s, P-S seismic surveying has developed rapidly, and P-S data is increasingly available (Stewart, 1994). Because more P-S data is being acquired, and because S interval velocity is an important indicator of rock lithology, it is desirable to have practical S inversion methods for P-S data. In this paper, one method will be described which uses constrained linearinversion to give an S interval velocity estimate. Description begins with an overview of the required input, followed by a chart of the basic inversion flow. The governing equations are described, and an inversion example is presented. INVERSION OVERVIEW The S velocity inversion presented here requires four sources of data: a P interval velocity in P-S time, an S interval velocity model in P-S time, an average α/β ratio, and seismic traces gathered into common-conversion points. Output is an S velocity section, where the amplitude of each sample is the S interval velocity at that point. The P-S reflectivity data should be scaled to represent true reflectivity, and be as broad-band, and noise-free as possible. NMO correction should be applied. The P velocity should be obtained from P velocity inversion, and should be stretched to tie the P-S data. Thus it is converted from a log in P-P time, to a log in P-S time. The average α/β is used during the inversion process to model angles at conversion points. If the input CCP gathers are asymptotically approximated, the Vp/Vs used in the approximation should be used in the inversion. The inversion flow is summarized in Figure 1. CREWES Research Report Volume 7 (1995) 18-1
2 Ferguson and Stewart P-SV Inversion Flow Chart P-S seismic section P-wave Velocity Correlate S-wave velocity model S-wave velocity model in P-S time P-wave velocity model in P-S time NMO corrected CCP gathers Inversion S-WAVE VELOCITY TRACES FIG. 1. Flow-chart for P-S inversion. RELATIVE CHANGE IN S-WAVE VELOCITY Estimation of S velocity from P-S seismic data requires a relationship between β and P-S reflectivity. A simple relationship, is the ratio of velocity change across an interface to the average velocity across the interface ( β/β), or normalized change in velocity. This value can be estimated using the following equations. At each conversion point: where R ps = 4c ρ ρ + d β β, (1) c = - αtan φ 8β 1-2β2 α 2 sin2 θ + 2β α cos θ cos φ, d = αtan φ 2β 4β 2 α 2 sin2 θ - 4β α cos θ cos φ, where Rps is the P-S reflection coefficient. The angles θ and φ are calculated as averages across the interface. The values α/α, β/β, and ρ/ρ are the relative changes in P-wave velocity, S-wave velocity, and density (Aki and Richards, 1980). Using an assumed relationship between ρ and α such as ρ = κ α 0.25 (Gardner et al, 1974), the P- S coefficient equation is: 18-2 CREWES Research Report Volume 7 (1995)
3 PS Inversion R ps = c α α + d β β. (2) Equation (2) can be cast as a least squares problem and solved for β/β (Stewart, 1990). The sum of the squares of the error at a single interface is: ε = R - R ps 2, (3) where R and R ps are the recorded and modeled P-S reflectivity respectively. Summation is over trace offset. Expansion of (3) gives the result: ε = R 2-2 α α Rc - 2 β β Rd + α2 α 2 c 2 + β2 β 2 d α β αβ cd (4) Differentiation of (4) with respect to β/β will give a value of β/β that minimizes the error function: ε β = - 2 Rd + 2 β β β d α α cd = 0. (5) Rearranging (5) to solve for β/β gives: β β = n offsets i = 1 R i ps d i - α α n offsets i = 1 n offsets i = 1 d i 2 c i d i. (6) CREWES Research Report Volume 7 (1995) 18-3
4 Ferguson and Stewart The angles θ and φ (required by coefficients c and d ) are approximated using RMS velocities and Snell s law (Dix, 1955; Zheng and Stewart, 1991; Slawinski and Slawinski, 1994). The RMS velocities and α/α can both be derived from the P-wave interval velocity model. Having established a relationship between P-S reflection data and the relative change in S velocity (equation 6), constrained linear-inversion can be used to convert β/β to a straight-forward velocity. CONSTRAINED LINEAR INVERSION Velocity estimates from P-P seismic inversion have often been compromised by a lack of low-frequency information in typical seismic data (Hendrick and Hearn, 1993). That is, velocity estimates exhibit low-frequency instability. In the P-P case, linear inversion can be stabilized by constraining the velocity solution to be close to an initial guess (Russell and Hampson, 1991). The result is that velocity information not resolved by the data is provided by the initial velocity estimate. Similarly, a constrained algorithm can be used to estimate actual velocity from the normalized velocity changes. The first step is to linearize β/β from equation (6). This can be done by scaling β/β with a macro-velocity model (7b). The change in S velocity ( β) across an interface can be represented by: β i = β i+1 - β i (7A) or, using β/β from equation (6), β i = β β -1 i α i α β mod (7B) Here, β/β has been modified to approximate β using the model P velocity and the model α/β. Equation 7a is the model response and the modified form of β/β (equation 7b) is the observed data. The model response (7a) is a linear function of β and as such, can be represented by the following first-order Taylor expansion (underlined symbols indicate vectors): β = β o + Z β - β o (8) 18-4 CREWES Research Report Volume 7 (1995)
5 PS Inversion where βο and βο are initial guesses of the change in S velocity across an interface and the S interval velocity respectively, β is the model interval velocity and Z is an n x n+1 matrix of the partial derivatives of 7a: Z = (9) where a single element of Z is: Z ij = β i β j Setting δ = β βο gives the change in S velocity vector: β = β o + Z δ (10) Let error vector e represent the difference between (10) and the observed data from equation (7b): e = β obs - β o - Zδ (11) The sum of the squares of (11) is: e T e = β obs - β o - Zδ T β obs - β o - Zδ (12) At this point a constraint λ is introduced to the squared error function in the form of a Lagrange multiplier problem (Lines and Treitel, 1984). CREWES Research Report Volume 7 (1995) 18-5
6 Ferguson and Stewart S = e T e + λ δ T δ - δ o 2 (13) where, S (eq. 13) is minimized subject to the constraint that the absolute value of δ is constant (Lines and Treitel, 1984). The use of λ has a twofold effect. Singularities in matrix inversion (equation 16) are avoided, and the resulting S velocity estimates are constrained to be close to an initial guess (Lines and Treitel, 1984). Differentiation of S (eq. 13) with respect to δ, and setting the result equal to zero, gives a modified form of the least-squared error: S δ = δ δt Z T T Zδ - β obs - β o Zδ - δ δt Z T β obs - β o + δ β T obs - β o βobs - β o + δ λ δt 2 δ - δ o (14) Solving for δ : δ = Z T Z + λi -1 Z T β obs - β o (15) Replacing δ with β βo: β = Z T Z + λi -1 Z T β obs - β o + β o (16) 18-6 CREWES Research Report Volume 7 (1995)
7 PS Inversion Equation (16) is the constrained linear-inversion equation. It ensures that the input model provides as much low-frequency information to the velocity estimates as the user thinks is required. P-S INVERSION APPLICATIONS The P-S inversion method developed above was used to provide an estimate of the S velocity and Vp/Vs for line EKW-002 of the Lousana 3-C data set. This line was acquired as part of a two line, 3-C survey by Unocal (Miller et al, 1994). Available to this study were the vertical and radial component seismic data, and two sonic logs (P sonics only). Also available were horizon interpretations and macro Vp/Vs ratios from the work of Miller et. al The seismic data sets were processed to preserve true amplitude, and a broad frequency-band. The first step in the inversion process was to derive an estimate of the P velocity for the entire line. Wells and (Figure 2) were used as initial P velocity models, and P-P inversion was accomplished using the Promax processing package, using a frequency-domain post-stack inversion (Figure 3). Then, the horizon picks, shown in Figure 4, were used to correlate the P velocity to P-S time (Figure 5); this step employed a sinc function interpolator to resample the stretched P velocity. Correlation of P-wave Sonics to P-P Stack VIK MISS WAB IRET FIG. 2. P sonics correlated to P-P seismic data, well 2-20 is on the left and well is on the right. The correlated logs were used as initial velocities for P-P inversion. Macro Vp/Vs ratios (See Miller et al, 1994, pg. 7-21) were used to scale the P-sonic logs to S-sonic logs; these were used as initial S-velocity input for the P-S inversion (Figure 6). The inversion input used for line 002 is summarized: 1) P-velocity from P-P inversion CREWES Research Report Volume 7 (1995) 18-7
8 Ferguson and Stewart 2) Correlation of P-wave velocity using horizon picks and sinc function resampling 3) Initial guess S-wave velocity from P-sonic logs scaled with macro Vp/Vs Figure 7 is the P-S inversion of the data. At this stage, issues regarding the proper scaling of the seismic data, and how tightly the S-velocity solution should have been constrained to the initial velocity, were not uniquely known. The P velocity from P-P inversion, and the S velocity from P-S inversion were then used to estimate a continuos Vp/Vs for the line (Figure 8). P-wave Velocity From P-P Inversion Viking Miss Wabmun Nisku Cooking L. FIG. 3. P interval velocity from post-stack inversion. Initial velocities from P sonics are overlain at and The time axis is ms, the cdp range is FIG. 4. Horizon picks. Picks were used to correlate the P-velocity section of Figure 3 to P-S time CREWES Research Report Volume 7 (1995)
9 PS Inversion P-wave Velocity Correlated to P-SV Time FIG. 5. P velocity from Figure 3 correlated to P-S time using the horizon picks. The time axis is ms, the cdp range is Correlation of Shear Sonics to PSV Stack VIK MISS WAB IRET Fig.6. Correlation of S-sonic logs to P-S data. The correlated logs were used as initial S velocities for P-S inversion. CREWES Research Report Volume 7 (1995) 18-9
10 Ferguson and Stewart S-wave Velocity Frome P-SV Inversion Viking Miss Wabmun Nisku Cooking L. FIG. 7. S velocity from P-S inversion. The S sonics, which were used as initial S- velocity estimates, are indicated at and VP/VS Ratio Nisku Cooking Lake VP/VS Ratio FIG. 8. Vp/Vs derived by dividing the P velocity (Figure 5) by the S (Figure 7). The time range is 1850 to 2100 ms to show greater detail. The Vp/Vs near 1 at the Nisku level is probably due to multiple interference originating in the Manville coals (Miller et al., 1994). An anomalous Vp/Vs increase to the right of can be seen. The cdp range is CREWES Research Report Volume 7 (1995)
11 PS Inversion DISCUSSION Both of the P and S velocity estimates (Figure 3 and 7 respectively) show a fairly consistent lateral velocity, with no sudden velocity changes. The S velocity estimate may be somewhat over-constrained to the initial guess velocity. This means that the ratio of reflection data to background velocity is to highly biased to the back ground velocity. The effect can be seen in the smooth nature of the S velocity estimate (Figure 7), compared with the P estimate (Figure 3). The Vp/Vs section of Figure 8 is included as an example of the continuous Vp/Vs display, that is possible with both P and P-S inversion estimates. Examination of Figure 8 shows the effects of multiple contamination at the Nisku level (Miller et al., 1994), the Vp/Vs is near 1, and a Vp/Vs decrease corresponding with the zone just above the Cooking Lake level. The multiple interference at the Nisku level probably originates in the Manville coals (Miller et al., 1994), and the Vp/Vs increase to the right of is as yet unexplained. CONCLUSION A method is presented which can be used to estimate S interval velocity from P- S seismic data. The method is composed of three steps: P-S data is converted to a relative change in velocity section, this section is scaled by a macro-velocity model to give a simple-change in velocity section, this last section is converted to S interval velocity using constrained linear-inversion. A field data set was inverted to derive S interval velocity and a Vp/Vs section. Contamination by multiples resulted in a Vp/Vs estimate of 1 at the Nisku level, and an anomalous increasing Vp/Vs was found above the Cooking lake level. These results are tempered by the non-unique nature, with regards to scaling and initial velocity constraint, of the inversion solution. More work must be done to determine if the results are valid, and how to set parameter boundaries accordingly. The authors wish to thank Henry C. Bland for his contributions to this work. We gratefully appreciate the continued support of the CREWES Project Sponsors. REFERENCES Aki, K. and Richards, P. G., 1980, Quantitative seismology: Theory and methods, v.1: W. H. Freeman and Co. Dix, H. C., 1955, Seismic velocities from surface measurements: Geophysics, 20, Gardener, G.H.F., Gardener, L. W., and Gregory, A. R., 1974, Formation velocity and density-the diagnostic basics of stratigraphic traps: Geophysics, 39, Hendrick, N., and Hearn, S., 1993, Evaluation of seismic trace inversion techniques: Exploration Geophysics, 24, Lines, L. R., and Treitel, S., 1984, A review of least-squares inversion and its application to geophysical problems: Geoph. Prosp., 32, Miller, S. L. M., Harrison, M. P., Lawton, D. C., Stewart, R. R., and Szata, K. J., Analysis of P-P and P-SV seismic data from Lousana Alberta: CREWES Project Research Report, v.7, Oldenburg, D. W., Scheuer, T., and Levy, S., 1983, Recovery of the acoustic impedance from reflection seismograms: Geophysics, 48, Russell, B., and Hampson, D., 1991, Comparison of post-stack seismic inversion methods: 61st Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, Slawinski, M. A, and Slawinski, R. A., 1994, Energy partition at the boundary between anisotropic media; part one: generalized Snell's law: CREWES Project Research Report, v.6, CREWES Research Report Volume 7 (1995) 18-11
12 Ferguson and Stewart Smith, G. C., and Gidlow, P. M., 1987, Weighted stacking for rock property estimation and detection of gas: Geophysical Prospecting, 35, Stewart, R. R., 1990, Joint P and P-SV inversion: CREWES Project Research Report, v.2, Stewart, R. R., 1991, Rapid map and inversion of P-SV waves: Geophysics, 56, Stewart, R. R., 1994, The present and promise of P-S seismic exploration: CREWES Project Research Report, v.6, Tatham, R. H., and McCormak, M. D, 1991, Multicomponent seismology in petroleum exploration: Society of Exploration Geophysicists. Zheng, C., and Stewart, R. R., 1991, Using P-SV waves to improve conventional AVO estimates: A synthetic study: CREWES Project Research Report, v.3, CREWES Research Report Volume 7 (1995)
A 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 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 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 informationModelling 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 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 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 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 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 informationFigure 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.
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
More informationRobert R. Stewart, Henry C. Bland, Qi Zhang*, and Felix Guthoff**
Average versus interval Vp/Vs Average versus interval Vp/Vs Robert R. Stewart, Henry C. Bland, Qi Zhang*, and Felix Guthoff** *Matrix Geoservices Ltd., Calgary; **Institute für Geophysik, Münster, Germany
More informationPost-stack inversion of the Hussar low frequency seismic data
Inversion of the Hussar low frequency seismic data Post-stack inversion of the Hussar low frequency seismic data Patricia E. Gavotti, Don C. Lawton, Gary F. Margrave and J. Helen Isaac ABSTRACT The Hussar
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 informationDelineating a sandstone reservoir at Pikes Peak, Saskatchewan using 3C seismic data and well logs
Delineating a sandston reservoir at Pikes Peak Delineating a sandstone reservoir at Pikes Peak, Saskatchewan using 3C seismic data and well logs Natalia L. Soubotcheva and Robert R. Stewart ABSTRACT To
More informationBandlimited impedance inversion: using well logs to fill low frequency information in a non-homogenous model
Bandlimited impedance inversion: using well logs to fill low frequency information in a non-homogenous model Heather J.E. Lloyd and Gary F. Margrave ABSTRACT An acoustic bandlimited impedance inversion
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 (AVO) and post-stack inversion of the Hussar low frequency seismic data
Pre-stack (AVO) and post-stack inversion of the Hussar low frequency seismic data A.Nassir Saeed, Gary F. Margrave and Laurence R. Lines ABSTRACT Post-stack and pre-stack (AVO) inversion were performed
More informationInterpretation of baseline surface seismic data at the Violet Grove CO 2 injection site, Alberta
Violet Grove seismic interpretation Interpretation of baseline surface seismic data at the Violet Grove CO 2 injection site, Alberta Fuju Chen and Don Lawton ABSTRACT Time-lapse seismic technology has
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 informationAnalysis of multicomponent walkaway vertical seismic profile data
Analysis of multicomponent walkaway vertical seismic profile data Bona Wu, Don C. Lawton, and Kevin W. Hall ABSTRACT A multicomponent walkaway VSP data processed for PP and PS imaging as well to study
More informationSensitivity of interval Vp/Vs analysis of seismic data
Sensitivity of interval Vp/Vs analysis of seismic data Rafael J. Asuaje* and Don C. Lawton CREWES, University of Calgary Summary The objective of this paper is to set guidelines for an interpreter to accurately
More informationMulticomponent seismic surveys at Sibbald Flats, Alberta
Multicomponent seismic surveys at Sibbald Flats, Alberta Robert R. Stewart and Don C. Lawton ABSTRACT We conducted a 3C seismic survey, using vibrators, over a 5 km line in the Sibbald Flats area some
More informationAn investigation of the free surface effect
An investigation of the free surface effect Nasser S. Hamarbitan and Gary F. Margrave, An investigation of the free surface effect ABSTRACT When P and S seismic waves are incident on a solid-air interface
More informationFeasibility and design study of a multicomponent seismic survey: Upper Assam Basin
P-276 Summary Feasibility and design study of a multicomponent seismic survey: Upper Assam Basin K.L.Mandal*, R.K.Srivastava, S.Saha, Oil India Limited M.K.Sukla, Indian Institute of Technology, Kharagpur
More informationMapping the conversion point in vertical transversely isotropic (VTI) media
Mapping the conversion point in vertical transversely isotropic (VTI) media Jianli Yang and Don. C. Lawton Conversion-point mapping ABSTRACT The important aspect of converted-wave (P-S) seismology is that
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 informationDetecting fractures using time-lapse 3C-3D seismic data
data Zimin Zhang, Don C. Lawton and Robert R. Stewart ABSTRACT This report presents the interpretation of time-lapse 3C-3D seismic data for fracture detection in a Saskatchewan potash mine. Seismic interpretation
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 informationInversion and interpretation of multicomponent seismic data: Willesden Green, Alberta
Willesden Green Inversion and interpretation of multicomponent seismic data: Willesden Green, Alberta Roxana M. Varga, Robert R. Stewart and John C. Bancroft ABSTRACT In this paper, 2D multicomponent seismic
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 informationSUMMARY INTRODUCTION LOCALIZED PHASE ESTIMATION
Local similarity with the envelope as a seismic phase detector Sergey Fomel, The University of Texas at Austin, and Mirko van der Baan, University of Alberta SUMMARY We propose a new seismic attribute,
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 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 informationStructure-constrained relative acoustic impedance using stratigraphic coordinates a
Structure-constrained relative acoustic impedance using stratigraphic coordinates a a Published in Geophysics, 80, no. 3, A63-A67 (2015) Parvaneh Karimi ABSTRACT Acoustic impedance inversion involves conversion
More informationProcessing and interpretation of 3C-3D seismic data: Clastic and carbonate numerical models
3C-3D Interpretation: Numerical Model Processing and interpretation of 3C-3D seismic data: Clastic and carbonate numerical models Glenn A. Larson 1 and Robert R. Stewart ABSTRACT Three-dimensional (3-D)
More informationPenn West Pembina Cardium CO 2 EOR seismic monitoring program
Penn West Pembina Cardium CO 2 EOR seismic monitoring program Don Lawton Marcia Coueslan, Fuju Chen Henry Bland, Abdullah Alshuhail University of Calgary Calgary, Alberta, Canada Penn West Petroleum CO
More informationRock Physics and Quantitative Wavelet Estimation. for Seismic Interpretation: Tertiary North Sea. R.W.Simm 1, S.Xu 2 and R.E.
Rock Physics and Quantitative Wavelet Estimation for Seismic Interpretation: Tertiary North Sea R.W.Simm 1, S.Xu 2 and R.E.White 2 1. Enterprise Oil plc, Grand Buildings, Trafalgar Square, London WC2N
More informationA Petroleum Geologist's Guide to Seismic Reflection
A Petroleum Geologist's Guide to Seismic Reflection William Ashcroft WILEY-BLACKWELL A John Wiley & Sons, Ltd., Publication Contents Preface Acknowledgements xi xiii Part I Basic topics and 2D interpretation
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 informationVelocity and VTI anisotropy scanning in multicomponent seismic data
PS multi-parameter scan Velocity and VTI anisotropy scanning in multicomponent seismic data Christopher O. Ogiesoba *, James E. Gaiser **, and Robert R. Stewart * ABSTRACT We present a prestack method
More informationThe Estimation of PS Wave Statics By Using SCDCC + Stack Power Method
The Estimation of PS Wave Statics By Using SCDCC + Stack Power Method ABSTRACT Summary Aimin Xue * Petrol Sound Geoservices. 8, MacEwan Ridge Close, N.W., Calgary, Alberta T3K3A7, Canada. aiminx@telus.net
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 informationRobust one-step (deconvolution + integration) seismic inversion in the frequency domain Ivan Priezzhev* and Aaron Scollard, Schlumberger
Robust one-step (deconvolution + integration) seismic inversion in the frequency domain Ivan Priezzhev and Aaron Scollard, Schlumberger Summary Seismic inversion requires two main operations relative to
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 informationP S-wave polarity reversal in angle domain common-image gathers
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? P S-wave polarity reversal in angle domain common-image gathers Daniel Rosales and James Rickett 1 ABSTRACT The change in the reflection
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 informationAutomatic time picking and velocity determination on full waveform sonic well logs
Automatic time picking and velocity determination on full waveform sonic well logs Lejia Han, Joe Wong, John C. Bancroft, and Robert R. Stewart ABSTRACT Full waveform sonic logging is used to determine
More informationCase History. 3-D AVO analysis and modeling applied to fracture detection in coalbed methane reservoirs. Antonio C. B. Ramos and Thomas L.
GEOPHYSICS, VOL. 62, NO. 6 (NOVEMBER-DECEMBER 1997); P. 1683 1695, 23 FIGS., 1 TABLE. Case History 3-D AVO analysis and modeling applied to fracture detection in coalbed methane reservoirs Antonio C. B.
More informationResidual Statics using CSP gathers
Residual Statics using CSP gathers Xinxiang Li and John C. Bancroft ABSTRACT All the conventional methods for residual statics analysis require normal moveout (NMO) correction applied on the seismic data.
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 informationABSTRACT INTRODUCTION
Migration velocity analysis for P-S data Miqrationvelocity analysis for P-S data Wai-Kin Chan and Robert R. Stewart ABSTRACT Migration velocity analysis (MVA) in depth (Reshef, 1992) has been successfully
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 informationAVAZ inversion for fracture orientation and intensity: a physical modeling study
AVAZ inversion for fracture orientation and intensity: a physical modeling study Faranak Mahmoudian*, Gary F. Margrave, and Joe Wong, University of Calgary. CREWES fmahmoud@ucalgary.ca Summary We present
More informationReservoir Characterization of Plover Lake Heavy-Oil Field
Reservoir Characterization of Plover Lake Heavy-Oil Field Larry Lines* Chorus, University of Calgary, Calgary, Alberta, Canada lrlines@ucalgary.ca Joan Embleton, Mathew Fay, Steve Larter, Tony Settari
More informationA031 Porosity and Shale Volume Estimation for the Ardmore Field Using Extended Elastic Impedance
A31 Porosity and Shale Volume Estimation for the Ardmore Field Using Extended Elastic Impedance A.M. Francis* (Earthworks Environment & Resources Ltd) & G.J. Hicks (Earthworks Environment & Resources Ltd)
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 informationTime lapse view of the Blackfoot AVO anomaly
Time lapse view of the Blackfoot AVO anomaly Han-xing Lu, Gary F. Margrave and Colin C. Potter Time lapse view of the Blackfoot AVO SUMMARY In the Blackfoot field, southeast of Calgary there is an incised
More informationShallow P and S velocity structure, Red Deer, Alberta
Shallow P and S velocity structure, Red Deer, Alberta P & S velocity structure Don C. Lawton, Meredith A. McArthur, Rachel T. Newrick and Sarah E. Trend ABSTRACT A multioffset vertical seismic profile
More informationSEG/New Orleans 2006 Annual Meeting
Carmen C. Dumitrescu, Sensor Geophysical Ltd., and Fred Mayer*, Devon Canada Corporation Summary This paper provides a case study of a 3D seismic survey in the Leland area of the Deep Basin of Alberta,
More informationINVERSION OF SEISMIC DATA FOR ASSESSING FLUID REPLACEMENT IN THE NISKU FORMATION PATRICIA GAVOTTI DON C. LAWTON
INVERSION OF SEISMIC DATA FOR ASSESSING FLUID REPLACEMENT IN THE NISKU FORMATION PATRICIA GAVOTTI DON C. LAWTON TABLE OF CONTENTS Objective Introduction Project Pioneer data Seismic processing overview
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 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 informationMapping the P-S Conversion Point in VTI Media. * Jianli Yang Don C. Lawton
Mapping the P-S Conversion Point in VTI Media * Jianli Yang Don C. Lawton Outline! Introduction! Theory! Numerical modeling methodology and results! NORSARD anisotropy ray mapping! Discussion and conclusions!
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 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 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 informationBaseline VSP processing for the Violet Grove CO 2 Injection Site
Baseline VSP processing for Violet Grove Baseline VSP processing for the Violet Grove CO 2 Injection Site Marcia L. Couëslan, Don C. Lawton, and Michael Jones * ABSTRACT Injection of CO 2 for enhanced
More informationTime-lapse seismic modelling for Pikes Peak field
Time-lapse seismic modelling for Pikes Peak field Ying Zou*, Laurence R. Bentley and Laurence R. Lines University of Calgary, 2500 University Dr, NW, Calgary, AB, T2N 1N4 zou@geo.ucalgary.ca ABSTRACT Predicting
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 informationPorosity prediction using cokriging with multiple secondary datasets
Cokriging with Multiple Attributes Porosity prediction using cokriging with multiple secondary datasets Hong Xu, Jian Sun, Brian Russell, Kris Innanen ABSTRACT The prediction of porosity is essential for
More informationThe relationship between elastic-wave velocities and density in sedimentary rocks: A proposal
260 The relationship between elastic-wave velocities and density in sedimentary rocks: A proposal Susan L.M. Miller and Robert R. Stewart ABSTRACT When information about density (0) is unavailable, it
More informationPluto 1.5 2D ELASTIC MODEL FOR WAVEFIELD INVESTIGATIONS OF SUBSALT OBJECTIVES, DEEP WATER GULF OF MEXICO*
Pluto 1.5 2D ELASTIC MODEL FOR WAVEFIELD INVESTIGATIONS OF SUBSALT OBJECTIVES, DEEP WATER GULF OF MEXICO* *This paper has been submitted to the EAGE for presentation at the June 2001 EAGE meeting. SUMMARY
More informationImpact of Phase Variations on Quantitative AVO Analysis
Impact of Phase Variations on Quantitative AVO Analysis Summary A.K. Srivastava, V. Singh*, D.N. Tiwary and V. Rangachari GEOPIC, ONGC, Dehradun Email: ak_sri3@rediffmail.com Effectiveness of AVO techniques
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 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 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 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 informationSeismic Inversion on 3D Data of Bassein Field, India
5th Conference & Exposition on Petroleum Geophysics, Hyderabad-2004, India PP 526-532 Seismic Inversion on 3D Data of Bassein Field, India K.Sridhar, A.A.K.Sundaram, V.B.G.Tilak & Shyam Mohan Institute
More informationInterpretation of PP and PS seismic data from the Mackenzie Delta, N.W.T.
Interpretation of PP and PS seismic data from the Mackenzie Delta, N.W.T. Carlos E. Nieto* CREWES. The University of Calgary nieto@geo.ucalgary.ca and Robert R. Stewart CREWES. The University of Calgary
More informationSimultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait
Simultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait Osman Khaled, Yousef Al-Zuabi, Hameed Shereef Summary The zone under study is Zubair formation of Cretaceous
More informationThe reason why acoustic and shear impedances inverted
SPECIAL The Rocky SECTION: Mountain The Rocky region Mountain region Comparison of shear impedances inverted from stacked PS and SS data: Example from Rulison Field, Colorado ELDAR GULIYEV, Occidental
More informationRadiation pattern in homogeneous and transversely isotropic attenuating media
Radiation pattern in homogeneous and transversely isotropic attenuating media Satish Sinha*, Sergey Abaseyev** and Evgeni Chesnokov** *Rajiv Gandhi Institute of Petroleum Technology, Rae Bareli, UP 229010
More informationSeismic methods in heavy-oil reservoir monitoring
Seismic methods in heavy-oil reservoir monitoring Duojun A. Zhang and Laurence R. Lines ABSTRACT Laboratory tests show that a significant decrease in acoustic velocity occurs as the result of heating rock
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 informationPorosity prediction using attributes from 3C 3D seismic data
Porosity prediction Porosity prediction using attributes from 3C 3D seismic data Todor I. Todorov, Robert R. Stewart, and Daniel P. Hampson 1 ABSTRACT The integration of 3C-3D seismic data with petrophysical
More informationMinimum entropy deconvolution with frequency-domain constraints
GEOPHYSICS, VOL. 59, NO. 6 (JUNE 1994); P. 938-945, 9 FIGS., 1 TABLE. Minimum entropy deconvolution with frequency-domain constraints Mauricio D. Sacchi*, Danilo R. Velis*, and Alberto H. Cominguez ABSTRACT
More informationAn Improved Dual Sensor Summation Method with Application to Four-Component (4-C) Seafloor Seismic Data from the Niger Delta
Earth Science Research; Vol. 4, No. 2; 2015 ISSN 1927-0542 E-ISSN 1927-0550 Published by Canadian Center of Science and Education An Improved Dual Sensor Summation Method with Application to Four-Component
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 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 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 informationAcoustic impedance inversion and CO 2 flood detection at the Alder Flats ECBM project
ECBM CO 2 flood detection Acoustic impedance inversion and CO 2 flood detection at the Alder Flats ECBM project Jason M. McCrank and Don C. Lawton ABSTRACT The 3D post-stack vertical component data from
More informationAn Improvement of Velocity Variation with Offset (VVO) Method in Estimating Anisotropic Parameters
DOI 10.7603/s40875-014-0003-0 GSTF Journal of Geological Sciences (JGS) Vol. No.1, July 015 An Improvement of Velocity Variation with Offset (VVO) Method in Estimating Anisotropic Parameters Ida Herawati,
More informationHampsonRussell. A comprehensive suite of reservoir characterization tools. cgg.com/geosoftware
HampsonRussell A comprehensive suite of reservoir characterization tools cgg.com/geosoftware HampsonRussell Software World-class geophysical interpretation HampsonRussell Software is a comprehensive suite
More informationRefraction analysis of the Blackfoot 2D-3C data
Refraction analysis of the Blackfoot 2D-3C data Refraction analysis of the Blackfoot 2D-3C data Jocelyn Dufour, and Don C. Lawton ABSTRACT A three-component 2D seismic line was acquired in the Blackfoot
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 informationSynthetic Seismogram A Tool to Calibrate PP & PS Seismic Data
P-475 Summary Synthetic Seismogram A Tool to Calibrate PP & PS Seismic P. Sugadha*, M. K. Jain, M. Singh, ONGC Conventional P wave technology will not meet all the requirements of the industry. Hence,
More informationOTC Paper Number 15248
OTC Paper Number 15248 The Vp/Vs Inversion Procedure: A Methodology for Shallow Water Flow (SWF) Prediction from Seismic Analysis of Multicomponent Data. Carlos Moreno* / The University of Oklahoma, John
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 informationDeterministic and stochastic inversion techniques used to predict porosity: A case study from F3-Block
Michigan Technological University Digital Commons @ Michigan Tech Dissertations, Master's Theses and Master's Reports 2015 Deterministic and stochastic inversion techniques used to predict porosity: A
More informationFull Waveform Modelling for Converted Waves Seismic Reflections in Mountainous and Marine Environment
Advances in Petroleum Exploration and Development Vol. 7, No. 2, 2014, pp. 21-29 DOI:10.3968/5089 ISSN 1925-542X [Print] ISSN 1925-5438 [Online] www.cscanada.net www.cscanada.org Full Waveform Modelling
More informationIdentification of Velocity Variations in a Seismic Cube Using Neural Networks
Identification of Velocity Variations in a Seismic Cube Using Neural Networks Dario Sergio Cersósimo, Claudia Ravazoli, Ramón García-Martínez PhD Program & Geophisic Group of the Astronomic & Geophisic
More information