Calculation of a Synthetic Gather Using the Aki-Richards Approximation to the Zoeppritz Equations

Total Page:16

File Type:pdf, Size:1020Kb

Calculation of a Synthetic Gather Using the Aki-Richards Approximation to the Zoeppritz Equations University of New Orleans ScholarWorks@UNO University of New Orleans Theses and Dissertations Dissertations and Theses Fall 12-15-2012 Calculation of a Synthetic Gather using the Aki-Richards Approximation to the Zoeppritz Equations Graham Ganssle [email protected] Follow this and additional works at: https://scholarworks.uno.edu/td Part of the Geophysics and Seismology Commons Recommended Citation Ganssle, Graham, "Calculation of a Synthetic Gather using the Aki-Richards Approximation to the Zoeppritz Equations" (2012). University of New Orleans Theses and Dissertations. 1566. https://scholarworks.uno.edu/td/1566 This Thesis is protected by copyright and/or related rights. It has been brought to you by ScholarWorks@UNO with permission from the rights-holder(s). You are free to use this Thesis in any way that is permitted by the copyright and related rights legislation that applies to your use. For other uses you need to obtain permission from the rights- holder(s) directly, unless additional rights are indicated by a Creative Commons license in the record and/or on the work itself. This Thesis has been accepted for inclusion in University of New Orleans Theses and Dissertations by an authorized administrator of ScholarWorks@UNO. For more information, please contact [email protected]. Calculation of a Synthetic Gather using the Aki-Richards Approximation to the Zoeppritz Equations A Thesis Submitted to the Graduate Faculty of the University of New Orleans in partial fulfillment of the requirements for the degree Master of Science In Applied Physics By Graham R. Ganssle B.S. University of New Orleans, 2009 December, 2012 Table of Contents List of Figures ............................................................................................................................................. iii Abstract ........................................................................................................................................................ iv Chapter 1: Introduction ................................................................................................................................ 1 Chapter 2: Background ................................................................................................................................. 3 Seismic Gathers ........................................................................................................................................ 3 Zoeppritz Equations and Linear Approximations .................................................................................... 5 Wireline Logging ....................................................................................................................................... 7 Chapter 3: Techniques and Methods ............................................................................................................ 9 Data ......................................................................................................................................................... 10 Wavelet Choice ....................................................................................................................................... 10 Normal Incidence Synthetic Trace ......................................................................................................... 11 Offset Incorporation and AVO ............................................................................................................... 13 Chapter 4: Analysis ...................................................................................................................................... 17 Algorithm Differences ............................................................................................................................ 18 1) Vertical Scale ............................................................................................................................... 18 2) Wavelets ..................................................................................................................................... 21 3) Sample Rates ............................................................................................................................... 22 4) Approximations ........................................................................................................................... 22 Comparison of the two synthetic gathers ............................................................................................. 23 Chapter 5: Conclusions and Future Work ................................................................................................... 27 Conclusions ............................................................................................................................................. 27 Future Work ............................................................................................................................................ 27 References .................................................................................................................................................. 28 Appendix ..................................................................................................................................................... 29 Vita .............................................................................................................................................................. 57 ii List of Figures Figure 1: Selected ray paths from source to receiver ................................................................................... 1 Figure 2: (left) wireline log data and (right) seismic data ............................................................................. 3 Figure 3: Seismic gather with offset angle in degrees on the horizontal axis and time in seconds on the vertical axis ................................................................................................................................................... 4 Figure 4: Graphical representation of an acoustic wave striking a boundary between two media with different acoustic impedances ...................................................................................................................... 5 Figure 5: Sonic logging tools. (left) short spaced, (right) long spaced .......................................................... 8 Figure 6: Map of southeast Louisiana with the well located at the tip of the red arrow ............................. 9 Figure 7: (left) Ricker wavelet r(t) with f = 1, (right) its Fourier transform R(ω) ........................................ 11 Figure 8: The normal trace plotted by using the ratio of reflected amplitude to incident amplitude for the reflection coefficient series of eq. 7 convolved with the Ricker wavelet of eq. 5. ..................................... 12 Figure 9: A zero offset trace plotted using the Aki-Richards equation ....................................................... 14 Figure 10: A synthetic gather for SEBB #39ST1 plotted using the authors algorithm. Horizontal scale is offset angle, with one trace for every five degrees. Vertical scale is data index number increasing downward ................................................................................................................................................... 16 Figure 11: A synthetic gather of SEBB #39ST1 plotted using IHS Kingdom ................................................ 17 Figure 12: A data index/depth map from the author's algorithm .............................................................. 19 Figure 13: The discrete derivative of the data in figure 12 ......................................................................... 19 Figure 14: A time/depth map from IHS Kingdom ....................................................................................... 20 Figure 15: The discrete derivative of the data in figure 14 ......................................................................... 21 Figure 16: (left) an Ormsby wavelet and (right) its Fourier transform ....................................................... 22 Figure 17: Synthetic gathers, (left) created using IHS Kingdom, (right) created using the author's algorithm ..................................................................................................................................................... 23 Figure 18: A cartoon of a class 1 AVO response (Hampson-Russell, 2005) ................................................ 25 Figure 19: A cartoon of a class 4 AVO response (Hampson-Russell, 2005) ................................................ 25 iii Abstract A synthetic seismic gather showing amplitude versus offset can be analyzed by the interpretive geophysicist to predict rock properties useful in oil exploration. For the research reported in this thesis, reflection coefficients derived from measured well log data are convolved with a Ricker wavelet to create a synthetic seismic trace. The Zoeppritz equations describe the propagation of an acoustic wave across an interface between two viscous media of different acoustic impedances with respect to increasing offset angle. The Aki-Richards linear approximation to the Zoeppritz equations is applied here to create a synthetic seismic gather with offset angles up to fifty degrees. The Aki-Richards approximation has not been used in this fashion prior to this research. The resulting gather is compared to a corresponding synthetic gather created using commercially available approximations and
Recommended publications
  • AVO Principles, Processing and Inversion
    AVO principles, processing and inversion AVO principles, processing and inversion Hong Feng and John C. Bancroft ABSTRACT In this paper, we will review the basic principles of AVO and introduce some amplitude preserving algorithms. Seismic data processing sequences for AVO analysis and some algorithms for AVO inversion will be discussed. Topics related to AVO, such as rock physics, AVO modeling, AVO interpretation, and AVO pitfalls are beyond the scope of this report. INTRODUCTION The AVO (amplitude variation with offset) technique assesses variations in seismic reflection amplitude with changes in distance between shot points and receivers. AVO analysis allows geophysicists to better assess reservoir rock properties, including porosity, density, lithology and fluid content. Around 1900, Knott and Zoeppritz developed the theoretical work necessary for AVO theory (Knott, 1899; Zoeppritz, 1919), given the P-wave and S-wave velocities along with the densities of the two bounding media. They developed equations for plane-wave reflection amplitudes as a function of incident angle. Bortfeld (1961) simplified Zoeppritz’s equation, making it easier to understand how reflection amplitudes depend on incident angle and physical parameters. Koefoed (1955) described the relationship of AVO to change in Poisson’s ratio across a boundary. Koefoed’s results were based on the exact Zoeppritz equation. The conclusions drawn by Koefoed are the basis of today’s AVO interpretation. Rosa (1976) and Shuey (1985) did breakthrough work related to predicting lithology from AVO, inspired by Koefoed’s article. Ostrander (1982) illustrated the interpretation benefits of AVO with field data, accelerating the second era of amplitude interpretation. Allen and Peddy (1993) collected seismic data, well information, and interpretation from the Gulf Coast of Texas and published a practical look at AVO that documented not only successes but also the interpretational pitfalls learned from dry holes.
    [Show full text]
  • Feasibility of Multi-Component Seismic for Mineral Exploration
    Western Australian School of Mines: Minerals, Energy and Chemical Engineering Feasibility of Multi-component Seismic for Mineral Exploration Pouya Ahmadi This thesis is presented for the Degree of Doctor of Philosophy of Curtin University December 2018 Declaration To the best of my knowledge and belief, this thesis contains no material previously published by any other person except where due acknowledgement has been made. This thesis contains no material, which has been accepted for the award of any other degree or diploma in any university. Signature: …………………………………………. Date: …………………………07/12/2018 To my wife, parents and sister. iii Abstract The mineral industry traditionally uses potential field and electrical geophysical methods to explore for subsurface mineral deposits. These methods are very effective when exploring for shallow targets. As shallow deposits are now largely depleted, alternative technologies need to be introduced to discover deeper targets or to extend known resources deeper. Only the seismic reflection method can reach the desired depths as well as provide the high-definition images of the subsurface that are required. Conventional seismic reflection methods make use of only one component (vertical) of the wave field (P-wave) which is very effective in resolving complex structural shapes but less so for rock characterisation. The proposition of this work is that the utilisation of the full vector field (3 component seismic data (3C)) could provide additional information that will help in rock identification, classification and detection of subtle structures, alteration and lithological boundaries that are of crucial importance for mineral exploration. Shear wave splitting, polarisation changes and velocity differences between polarised modes have at present undefined potential in resolving lithological changes, internal composition of shear zones, subtle faults and fractured zones.
    [Show full text]
  • Avo Analysis
    CHAPTER ONE 1.1 INTRODUCTION TO AVO ANALYSIS 1.1.1BACKGROUND Seismic amplitude-versus offset (AVO) analysis has been a powerful geophysical method in aiding the direct detection of gas from seismic records. This method of seismic reflection data is widely used to infer the presence of hydrocarbons. The conventional analysis is based on such a strategy: first, one estimates the effective elastic parameters of a hydrocarbon reservoir using elastic reflection coefficient formulation; which models the reservoir as porous media and infers the porous parameters from those effective elastic parameters. Strictly speaking, this is an approximate approach to the reflection coefficient in porous media. The exact method should be based on the rigorous solutions in porous media. The main thrust of AVO analysis is to obtain subsurface rock properties using conventional surface seismic data. The Pooisson‟s ratio change across an interface has been of particular interest. These rock properties can then assist in determining lithology, fluid saturants, and porosity. It has been shown through solution of the Knott energy equations (or Zoeppritz equations) that the energy reflected from an elastic boundary varies with the angle of incidence of the incident wave (Muskat and Meres, 1940). This behavior was studied further by Koefoed (1955, 1962): He established in 1955 that, the change in reflection coefficient with the incident angle is dependent on the Poisson‟s ratio difference across an elastic boundary. Poisson‟s ratio is defined as the ratio of transverse strain to longitudinal strain (Sheriff, 1973), and is related to the P-wave and S-wave velocities of an elastic medium by equation 1 below.
    [Show full text]
  • Msc Thesis Justin Lo 1047357.Pdf (12.23Mb)
    A Modern Reinvestigation of the 1987 Moutere Depression Seismic Survey Justin Chen Lo A thesis submitted for the degree of Master of Science At the University of Otago, Dunedin, New Zealand. Table of Contents Abstract .................................................................................................................................................. iv Acknowledgements ................................................................................................................................ vi 1.0 Introduction ...................................................................................................................................... 1 1.1 Project Aim .................................................................................................................................... 1 1.2 Geological Background ................................................................................................................. 2 1.2.1 Moutere Depression .............................................................................................................. 2 1.2.2 Moutere Gravels & Glen Hope Formation ............................................................................. 2 1.2.3 Western Province ....................................................................................................................... 3 1.2.4 Eastern Province .................................................................................................................... 3 1.2.5 South of Nelson .....................................................................................................................
    [Show full text]