Attenuation of Multiples in Marine Seismic Data. Geology /23/06
|
|
- Priscilla Wilson
- 5 years ago
- Views:
Transcription
1 Attenuation of Multiples in Marine Seismic Data Geology /23/06
2 Marine Seismic Surveying Marine seismic surveying is performed to support marine environmental and civil engineering projects in the coastal zone. Also performed to evaluate and/or explore for new oil and gas reserves in the continental shelf.
3 Types of Surveys Single Channel Profiling: The source and receiver is contained in the same unit. Typically utilized in the upper 10 2 meters for high resolution engineering and environmental surveys. Unit can be deployed quickly and retrofitted to most any vessel. Measures the reflectivity of a wave with an angle of incidence of 0 degrees. Normal to a horizontal layer.
4 Types of Surveys Towed source and hydrophone arrays Utilized for exploration of oil and gas. Strong sources and hundreds of hydrophones peer deeper into the continental shelf.
5 Types of Surveys Multibeam acoustic seafloor mapping High frequency source sweeps the seafloor in a fan shape below the vessel. Not used for seafloor penetration. Used for seafloor contours.
6 Types of Surveys Side Scan Sonar peers laterally to determine relief of a feature, or the seafloor. Towfish flown in the water column. Try to keep the fish at a constant height above the seafloor.
7 Problem Multiples occur in marine seismic data due to the reflectivity of the water layer
8 Review of Reflection Phenomena Reflection Coefficient (R): Ratio of reflected and incident wave amplitudes R = Ar/Ai Ar = Amplitude of the reflected wave Ai = Amplitude of the incident wave Relates to the magnitude of reflection from the interface between two media with different physical properties
9 Acoustic Impedance (Z) Acoustic Impedance is a ratio of acoustic pressure to flow Z = ρv ρ=density of the material (kg/m 3 ) v=speed of the acoustic wave (m/s) Units of Rayles (kg/m 2 /s)
10 Reflection Coefficient (R) The full expression for the reflection coefficient : R=(Z 2 /Z 1 )- (1-(n-1)tan 2 α i ) (Z 2 /Z 1 )+ (1-(n-1)tan 2 α i ) Reflection occurs to some extent at every interface. n = the ratio of (v 2 /v 1 )^2, α is the angle of incidence of the wave ray
11 Reflection Coefficients in the Water Layer Let s take a look at the simple case of the reflection coefficient at the air/water interface and water/rock interface for a wave of normal incidence. air: ρ=1.3 kg/m 3 v=350 m/s seawater: ρ= 1027 kg/m 3 v=1500 m/s rock: ρ= 2650 kg/m 3 v=3000 m/s
12 Quick Computation of the Reflection Coefficients Air/Water Interface: 2 1 Consider the cartoon: The source is located in the water layer, so the wave is propagating toward the air/water interface upward R = (Z2-Z1)/(Z2+Z1) = (1.3*350)-(1027*1500) (1.3*350)+(1027*1500) =
13 Quick Computation of the Reflection Coefficients Water/Rock Interface: 1 2 Consider the cartoon: The source is located in the water layer, so the wave is propagating toward the rock/water interface downward R = (Z2-Z1)/(Z2+Z1) = (2650*3000)-(1027*1500) (2650*3000)+(1027*1500) = 0.67
14 Characteristics of the Reflection Coefficients The air/water boundary can be considered a perfect reflector. The -1 coefficient indicates that the reflection reverses polarity at the interface. The water/rock boundary has a coefficient of 0.67 indicating that the boundary allows transmission of energy, but can also be considered a relatively efficient reflector. The energy does not reverse polarity at this interface. So for any given shot, it is clear that there is energy in the water layer that represents primary reflections and multiple reflections.
15 Types of Multiples Multiples- Events that have undergone more than one reflection Short Path Multiples- Arrives so soon after the primary reflection that it adds a tail to the primary reflection. Long Path Multiples- Reflection between surface and water bottom over a great distance. Looks like a separate event/layer on the seismogram.
16 Short Path Multiples Peg-Leg Multiples- Short path multiples that have been reflected from the top and base of thin reflectors. Delays part of the energy from the primary reflection which effectively lengthens the wavelet. Have the same polarity as the primary. Path: Intuitively you can see that the path is slightly longer, so there is a delay in the arrival.
17 Short Path Multiples Ghosts- Occur in marine surveys when a submerged source is used. Energy travels up from source to the air/water interface, reflects, and continues down as a secondary source. As previously mentioned, the reflection is reversed in polarity from the original source.
18 Short Path Multiples Water Reverberation- Due to large reflection coefficients at the air/water and water/rock interface. This type of multiple also occurs as a long path multiple over long distances.
19 Filtering Basics Transforms- Operations that transform a function of certain variables to a related function of different variables (i.e. from time domain to the frequency domain by the Fourier transform). Transforms a vector to a different vector space for easy computation.
20 Filtering Basics Convolution- A mathematical operator which takes as inputs two functions f and g and produces a third function that in a sense represents the amount of overlap between f and a reversed and translated version of g in the time domain. Convolution is the multiplication of two functions in the Z domain or frequency domain.
21 Filtering Basics Correlation- Indicates the strength and direction of a linear relationship between the inputs of two functions, producing a third output function.
22 Transforms Z Transform- Converts a discretely sampled time function to a polynomial in Z the Z domain. Consider the continuous time function:
23 Z Transform
24 Z Transform The continuous function has been discretely sampled in time steps of 1 unit. The Z transform takes the amplitude of the value of the function as the coefficient of the term of a term in the polynomial. The time unit becomes the exponent of the polynomial term. x t
25 Z Transform The discretely sampled function can be viewed as: A vector: b t = (0.5000, , , , , , , , , , ) A polynomial in the Z domain: B(Z)= Z Z Z Z Z Z Z Z Z Z 10 General form of the Z transform for a discretely sampled function b(t): B(Z)= t b t Z t
26 Features of the Z Transform The Z transform is also known as the unit delay operator. If B(Z)=1+2 Z +0Z 2 -Z 3 -Z 4 A time delay of n time units can be applied to the function by multiplying the polynomial by Z n.
27 Features of the Z Transform B(Z)=1+2Z+0Z 2 -Z 3 -Z 4 Time delayed function Delay by 1 time unit (Z 1 ) ZB(Z)=Z+2Z 2 +0Z 3 -Z 4 -Z 5
28 Features of the Z Transform The Z transform may also be used to create more complicated time functions from simple ones. Say B(Z) represents a source that you can see on a seismogram. B(Z) is said to be the impulse response of the source. We can add another source 5 time units later by simply creating a new function: C(Z)=B(Z)+Z 5 B(Z)
29 Features of the Z Transform If C(Z)=B(Z)+Z 5 B(Z) with B(Z)=0+2Z-1Z 2-1Z 3 +0Z 4 +Z 5 +Z 6 C(Z)= 0+2Z-1Z 2-1Z 3 +0Z 4 +Z 5 +3Z 6-1Z 7-1Z 8 +0Z 9 +Z 10 +Z 11 Notice that the signal has essentially repeated itself. However, look at the 3Z 6 term. Two terms were combined here, Z 5 (2Z)+Z 6. If pulses overlap with each other, the waveforms add together. This is called the linearity assumption.
30 Features of the Z Transform The waves do not interfere with each other despite traveling through heterogeneous materials. Nonlinearity arises from very large amplitude disturbances- not in the case of engineering and exploration seismology.
31 Features of the Z Transform The Taylor Series expansion of f(x) is given by the inverse of the Z transform: The Taylor series:
32 Linear System General form of a linear system: Y(Z) = X(Z)B(Z) where Y(Z) is the output, X(Z) is the input, and B(Z) is the impulse response (if the input is the unit impulse). Input Impulse Response Output
33 Convolution A convolution is an integral (or sum, if the function is discretely sampled) that expresses the amount of overlap of one function A as it is reversed and shifted over another function B. Numerically, the calculation for convolution of two functions can be performed by time reversing one function and shifting it relative to the other, summing the incident values. Take a look at an example:
34 Convolution Convolve A with B, or C = A*B A = (a 0, a 1, a 2, a 3 ) B = (b 0, b 1, b 2 ) a 0 a 1 a 2 a 3 b 2 b 1 b 0 c 0 = a 0 b 0
35 Convolution a 0 a 1 a 2 a 3 b 2 b 1 b 0 c 1 =a 0 b 1 +a 1 b 0 a 0 a 1 a 2 a 3 b 2 b 1 b 0 c 3 =a 1 b 2 +a 2 b 1 +a 3 b 0 a 0 a 1 a 2 a 3 b 2 b 1 b 0 c 2 =a 0 b 2 +a 1 b 1 +a 2 b 0 a 0 a 1 a 2 a 3 b 2 b 1 b 0 c 4 =a 2 b 2 +a 3 b 1
36 Convolution a 0 a 1 a 2 a 3 b 2 b 1 b 0 c 5 =a 3 b 2 The function C is: C = (c 0 c 1 c 2 c 3 c 4 c 5 ) c 0 = a 0 b 0 c 1 = a 0 b 1 +a 1 b 0 c 2 = a 0 b 2 +a 1 b 1 +a 2 b 0 c 3 = a 1 b 2 +a 2 b 1 +a 3 b 0 c 4 = a 2 b 2 +a 3 b 1 c 5 = a 3 b 2
37 Convolution Return to the Z transform briefly: A = a 0 + a 1 Z+ a 2 Z 2 + a 3 Z 3 B = b 0 +b 1 Z+ b 2 Z 2 Multiply the two polynomials: C(Z) = A(Z)B(Z) C=a 0 b 0 +a 0 b 1 Z+a 0 b 2 Z 2 +a 1 b 0 Z+a 1 b 1 Z 2 +a 1 b 2 Z 3 +a 2 b 0 Z 2 +a 2 b 1 Z 3 +a 2 b 2 Z 4 +a 3 b 0 Z 3 +a 3 b 1 Z 4 +a 3 b 2 Z 5
38 Convolution Collection of like terms: c 0 = (a 0 b 0 ) c 1 (Z) = (a 0 b 1 +a 1 b 0 )Z c 2 (Z 2 )= (a 0 b 2 +a 1 b 1 +a 2 b 0 )Z 2 c 3 (Z 3 )= (a 1 b 2 +a 2 b 1 +a 3 b 0 )Z 3 c 4 (Z 4) = (a 2 b 2 +a 3 b 1 )Z 4 c 5 (Z 5 )= (a 3 b 2 )Z 5
39 Convolution From the product of the multiplication of the two functions in the Z domain, we can see that convolution in the time domain is equivalent to multiplication in the Z domain.
40 Convolution Properties: Commutativity f*g = g*f Associativity f*(g*h) = (f*g)*h Associativity with scalar multiplication a(f*g) = (af)*g = f*(ag) Distributivity f*(g+h) = f*g + f*h Differentiation (f*g) = f *g = f*g
41 Correlation The correlation, or cross-correlation, is a measure of similarity of two signals as a function of the relative time shift between them, commonly used to find features in an unknown signal by comparing it to a known one. Auto-correlation is the correlation of a data set with itself. Numerically, the correlation procedure is similar to convolution, except the signals are not time reversed with respect to each other. One signal is shifted with respect to the second signal. The amplitudes of the two signals are scaled with respect to each other and summed. Take a look at an example:
42 Correlation Correlate A with B, or Φ = A B A = (a 0, a 1, a 2, a 3 ) B = (b 0, b 1, b 2 ) Continuing to shift the function by the specified unit of time yields: a 0 a 1 a 2 a 3 b 0 b 1 b 2 Φ ab (0) = a 0 b 0 +a 1 b 1 +a 2 b 2 Φ ab (2) = a 2 b 0 +a 3 b 1 Φ ab (3) = a 3 b 0 a 0 a 1 a 2 a 3 b 0 b 1 b 2 Φ ab (1) = a 1 b 0 +a 2 b 1 +a 3 b 2
43 Correlation The general form of the cross correlation of two data sets, x t and y t : Φ xy (τ) = k x k y k+τ Φ is a function of the relative time between two signals Τ (tau) is the relative shift of the signal y K is the element of each signal that we are looking at
44 Correlation Intuitively, we can see that two similar data sets will have positive products with large sums. Conversely, two dissimilar data sets will have both positive and negative products with small sums. If two data sets are correlated by having a large negative sum, it means the data sets would be similar if one was inverted (the two data sets are similar, but out of phase with respect to each other).
45 Multiple Removal Deterministic Predictive Filtering
46 Deterministic Predictive Filtering Deghosting with predictive filters Recall the geometry of the ghost:
47 Deterministic Predictive Filtering The ghost creates the equivalent of a second apparent shot. The second source follows the initial shot by τ milliseconds, which is a function of the shot depth. Since the ghost is reflected off the free surface (air/water interface) it exhibits 180 degree shift in polarity.
48 Deterministic Predictive Filtering The simplified impulse response of the ghost is shown to the right. The first element is the unit impulse. The second element is the ghost. τ is the time lag of the ghost k is the amplitude of the reflected energy
49 Deterministic Predictive Filtering The objective is to define a filter such that the input is reduced to the desired output, which is the unit spike of energy with no multiples. Filter 1 0
50 Deterministic Predictive Filtering x is the input y is the desired output f is the filter that gets from x to y In the time domain we convolve the input with the filter to get the desired output x*f = y
51 Deterministic Predictive Filtering Recall that convolution in the time domain is multiplication in the Z domain X(z)F(z) = Y(z) Filter design: F(z) = Y(z)/X(z)
52 Deterministic Predictive Filtering F(Z) = Y(Z)/ X(Z) X(Z) = 1 kz τ Y(Z) = 1 F(Z) = 1/(1 kz τ ) Note that the Taylor series is the inverse of the Z transform, so Filter 0
53 Deterministic Predictive Filtering F(Z) = 1/(1 kz τ ) F(Z) = 1 + kz τ k 2 z 2τ +k 3 z 3τ This filter has infinite length, so let s truncate it after the first 2 points and convolve it with X(Z)
54 Deterministic Predictive Filtering Y(Z) = F(Z)*X(Z) f o x o = 1 f t x o + f o x t = k + -k = 0 The multiple has been attenuated, but there is still one more term left to deal with
55 Deterministic Predictive Filtering The ghost has been eliminated but a new false event has been created at 2τ. The false event has amplitude of k 2. By increasing the number of terms in the filter, the amplitude of the false event gets smaller and smaller. Last term of the convolution
56 Deterministic Predictive Filtering Predictive Filtering to Attenuate Reverberation Effects After the source stops putting energy into the water layer, the reflections continue and decrease in amplitude. Due to the large acoustic impedances at the interfaces (i.e. air/water and water/rock).
57 Deterministic Predictive Filtering What is actually happening: The energy from the source propagates toward the air/water interface, reflects with a 180 degree phase shift, propagates toward the seafloor, transmits some energy and reflects again. The energy in the water layer now has an amplitude of k, which is a function of the reflection coefficient and impulse amplitude. The wave now reflects off the sea surface boundary and is phase shifted again. This continues for several cycles of the water layer.
58 Deterministic Predictive Filtering Reverberation has the form: 1,-2k,3k 2 at the receiver All of the terms with the exception of the first one represent multiples. Design a filter similar to that used in the deghosting technique.
59 Deterministic Predictive Filtering Up until now we have only looked at multiples from the source point of view. Multiples occur at both the source and receiver locations. First order reverberation (two types).
60 Deterministic Predictive Filtering Second Order Reverberation (three types). The goal of the filter is to take out the effect of the water layer and maintain the information of the trace below the water/rock interface. As can be seen in the accompanying cartoon, the energy travels through the water column two round trips. 2 nd Order Reverb
61 Deterministic Predictive Filtering Besides the two round trips through the water column, the energy also goes into the earth. The information from the sub-bottom areas is preserved. 2nd Order Reverb
62 Deterministic Predictive Filtering The term for the second order reverberation becomes 3k 2 z 2τ. The coefficient 3 because there are three options for the wave energy to follow. 2 nd Order Reverb
63 Deterministic Predictive Filtering The effect of the filter term 3k 2 z 2τ is: Which removes the information from the red boxes only and preserves information about the sub-bottom.
64 Deterministic Predictive Filtering Also third, fourth, etc. Need to account for these in filter design because they rescale (make the coefficients of the k term larger) the amplitudes of the subsequent reverberations
65 Deterministic Predictive Filtering The signal that we will have to filter is actually: X(Z) = 1-2kz τ +3k 2 z 2τ
66 Deterministic Predictive Filtering The first unit spike still represents the primary energy The filter for dereverberation effects: X(Z) F(Z) = Y(Z) X(Z) = 1-2kz τ +3k 2 z 2τ F(Z) = 1/ (1-2kz τ +3k 2 z 2τ ) F(Z) = 1+2kz τ +k 2 z 2τ
67 Deterministic Predictive Filtering F(Z) = 1+2kz τ +k 2 z 2τ f is: This is called the three point Backus filter
68 Deterministic Predictive Filtering The complete water reverberation function X(Z) that we just defined can be viewed from a different standpoint: Since the energy goes up and down through the water layer: Convolve the unit spike with the simple water layer model that we defined earlier for the downgoing wave. Convolve again for the upgoing wave energy to the surface.
69 Deterministic Predictive Filtering Y*X = (1)*(1,-k,k 2,-k 3,k 4,-k 5,k 6 ) Convolve the result with itself, since it is going back through the water column. (1,-k,k 2,-k 3 )* (1,-k,k 2,-k 3 ) and transform to the Z domain We get 1-2kz τ +3k 2 z 2τ -4k 3 z 3τ which is the complete effect of reverberation in the water layer.
70 Deterministic Predictive Filtering Similarly, if we take the filter in the Z domain for the simple water layer and square it, we get the filter for the complete water layer. Remember that multiplication in the Z domain is equivalent to convolution in the t domain.
71 Performance of the Predictive Filter To check the performance of the filter, convolve the reverberation water layer model with the filter in steps of τ, the time reverberation time lag. X*F=Y Y= (1,-2k,3k 2,-4k 3,5k 4,-6k 5 )*(1,2k,k 2 )
72 Performance of the Predictive Filter y0 = 1(1) = 1 y1 = 1(-2k)+2k(1) = 0 y2 = 1(3k 2 )+2k(-2k)+k 2 (1) = 0 y3 = 1(-4k 3 )+2k(3k 2 )+k 2 (-2k) = 0.. So we see that the filter removes the effect of the reverberation related multiples.
73 Determining Filter Parameters To determine the time lag, τ, between the primary and the ghost, correlate a trace with respect to itself. This form of correlation is called auto correlation. The auto correlation of a trace will determine the time lag, τ, of the reverberation. After performing the operation, the correlation will exhibit a large trough. Remember that correlation is the multiplication and summation of the amplitudes of two signals, while one function is being time shifted past the other. The time lag between the primary and the reverberation coincides with the time shift that occurred to get the large trough.
74 Determining Filter Parameters To determine k, the amplitude of the reflection, you must determine the reflection coefficient of the multiple energy at the seafloor. The inverted amplitude of the first reflection can be divided by the amplitude of the impulse spike to get the reflection coefficient. k = (amplitude of the impulse)(r)
75 Determining Filter Parameters For simplicity, consider a system with a time step of 1 unit. Consider the trace to the right representing a very simple signal.
76 Determining Filter Parameters Find the time lag of the reverberation by autocorrelation of the simple signal seen in the previous slide. ( ) ( ) It appears that the time lag is 7 time units between the source and the reverberation. ( ) ( )
77 Determining Filter Parameters Now determine k Amplitude of the impulse is 3, the inverted value of the first multiple is 2 R = 2/3 k = 3R
78 Determining Filter Parameters Since this is deterministic filter prediction we can find the parameters by the methods above, but this gets to be time consuming to evaluate from trace to trace. We can get best fit values of the two parameters by implementing a least squares inversion and applying them to the stacked data.
79 Assignment Define a filter and implement on the following stacked profile:
The Analysis of Data Sequences in the Time and Frequency Domains 1
The Analysis of Data Sequences in the Time and Frequency Domains 1 D. E. Smylie 1 c D. E. Smylie, 2004. Contents 1 Time Domain Data Sequences 1 1.1 Discrete, Equispaced Data Sequences....... 1 1.2 Power
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 informationEstimating received sound levels at the seafloor beneath seismic survey sources
Proceedings of ACOUSTICS 016 9-11 November 016, Brisbane, Australia Estimating received sound levels at the seafloor beneath seismic survey sources Alec J Duncan 1 1 Centre for Marine Science and Technology,
More informationChapter 7: Reflection Seismology Homework Solutions (Jan. 2010)
Chapter 7: eflection Seismology Homework Solutions (Jan. 200). Why do marine seismic reflection surveys not record (a) S waves? (b) refracted rays? 2 μ a) For ideal fluid, μ=0, thus, v s = = 0 ρ b) eflection
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 informationSurface Waves and Free Oscillations. Surface Waves and Free Oscillations
Surface waves in in an an elastic half spaces: Rayleigh waves -Potentials - Free surface boundary conditions - Solutions propagating along the surface, decaying with depth - Lamb s problem Surface waves
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 informationLast Time. GY 305: Geophysics. Seismology (Marine Surveys) Seismology. Seismology. Other Seismic Techniques UNIVERSITY OF SOUTH ALABAMA
UNIVERSITY OF SOUTH ALABAMA Last Time GY 305: Geophysics Lecture 12: Introduction to (resolution versus penetration) Techniques (marine versus terrestrial) (Marine Surveys) http://www.glossary.oilfield.slb.com/displayimage.cfm?id=236
More informationDeconvolution imaging condition for reverse-time migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 83 96 Deconvolution imaging condition for reverse-time migration Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT The reverse-time
More informationResearch Project Report
Research Project Report Title: Prediction of pre-critical seismograms from post-critical traces Principal Investigator: Co-principal Investigators: Mrinal Sen Arthur Weglein and Paul Stoffa Final report
More informationMagnetotelluric (MT) Method
Magnetotelluric (MT) Method Dr. Hendra Grandis Graduate Program in Applied Geophysics Faculty of Mining and Petroleum Engineering ITB Geophysical Methods Techniques applying physical laws (or theory) to
More informationMandatory Assignment 2013 INF-GEO4310
Mandatory Assignment 2013 INF-GEO4310 Deadline for submission: 12-Nov-2013 e-mail the answers in one pdf file to vikashp@ifi.uio.no Part I: Multiple choice questions Multiple choice geometrical optics
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu.50 Introduction to Seismology Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. .50 Introduction to Seismology
More informationCompensating for attenuation by inverse Q filtering. Carlos A. Montaña Dr. Gary F. Margrave
Compensating for attenuation by inverse Q filtering Carlos A. Montaña Dr. Gary F. Margrave Motivation Assess and compare the different methods of applying inverse Q filter Use Q filter as a reference to
More informationDownloaded 09/17/13 to Redistribution subject to SEG license or copyright; see Terms of Use at Log data.
Downloaded 9/17/13 to 99.186.17.3. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ Extracting polar antropy parameters from mic data and well logs Rongrong
More informationSeismic applications in coalbed methane exploration and development
Seismic applications in coalbed methane exploration and development Sarah E. Richardson*, Dr. Don C. Lawton and Dr. Gary F. Margrave Department of Geology and Geophysics and CREWES, University of Calgary
More information7.2.1 Seismic waves. Waves in a mass- spring system
7..1 Seismic waves Waves in a mass- spring system Acoustic waves in a liquid or gas Seismic waves in a solid Surface waves Wavefronts, rays and geometrical attenuation Amplitude and energy Waves in a mass-
More informationLecture 04: Discrete Frequency Domain Analysis (z-transform)
Lecture 04: Discrete Frequency Domain Analysis (z-transform) John Chiverton School of Information Technology Mae Fah Luang University 1st Semester 2009/ 2552 Outline Overview Lecture Contents Introduction
More informationSeismic Imaging. William W. Symes. C. I. M. E. Martina Franca September
Seismic Imaging William W. Symes C. I. M. E. Martina Franca September 2002 www.trip.caam.rice.edu 1 0 offset (km) -4-3 -2-1 1 2 time (s) 3 4 5 How do you turn lots of this... (field seismogram from the
More informationPART A: Short-answer questions (50%; each worth 2%)
PART A: Short-answer questions (50%; each worth 2%) Your answers should be brief (just a few words) and may be written on these pages if you wish. Remember to hand these pages in with your other exam pages!
More informationPEAT SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II
PEAT8002 - SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II Nick Rawlinson Research School of Earth Sciences Australian National University Waveform modelling P-wave first-motions
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 informationBasic principles of the seismic method
Chapter 2 Basic principles of the seismic method In this chapter we introduce the basic notion of seismic waves. In the earth, seismic waves can propagate as longitudinal (P) or as shear (S) waves. For
More informationAmplitude, Frequency and Bandwidth and their relationship to Seismic Resolution
Environmental and Exploration Geophysics II Amplitude, Frequency and Bandwidth and their relationship to Seismic Resolution tom.h.wilson wilson@geo.wvu.edu Department of Geology and Geography West Virginia
More informationSeismic reflection notes
1 of 13 20/02/2006 3:29 PM Seismic reflection notes Introduction This figure was modified from one on the Lithoprobe slide set. LITHOPROBE is "probing" the "litho" sphere of our continent, the ground we
More informationEOS 350 MIDTERM OCT 4, 2013 STUDENT NAME: TEAM #:
EOS 350 MIDTERM OCT 4, 2013 STUDENT NAME: TEAM #: Some equations which may, or may not, be useful: Distance from sensor to a dipole z ~ x ½, Distance to line of dipoles z ~ 0.75x ½ B = μh, M = κh Seismic
More informationConvolution and Linear Systems
CS 450: Introduction to Digital Signal and Image Processing Bryan Morse BYU Computer Science Introduction Analyzing Systems Goal: analyze a device that turns one signal into another. Notation: f (t) g(t)
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 informationUpscaling Aspects of Spatial Scaling
Aspects of Spatial Scaling 288 We need to relate measurements at different scales Lab Logs Crosswell VSP Surface Seismic How does laboratory rock physics apply to the field? frequency differences sample
More informationAcoustic Velocity, Impedance, Reflection, Transmission, Attenuation, and Acoustic Etalons
Acoustic Velocity, Impedance, Reflection, Transmission, Attenuation, and Acoustic Etalons Acoustic Velocity The equation of motion in a solid is (1) T = ρ 2 u t 2 (1) where T is the stress tensor, ρ is
More informationOBS wavefield separation and its applications
P-088 OBS wavefield separation and its applications Sergio Grion*, CGGVeritas Summary This paper discusses present trends in ocean-bottom seismic (OBS) data processing. Current industrial practices for
More informationElements of 3D Seismology Second Edition
Elements of 3D Seismology Second Edition Copyright c 1993-2003 All rights reserved Christopher L. Liner Department of Geosciences University of Tulsa August 14, 2003 For David and Samantha And to the memory
More informationPhysics General Physics. Lecture 25 Waves. Fall 2016 Semester Prof. Matthew Jones
Physics 22000 General Physics Lecture 25 Waves Fall 2016 Semester Prof. Matthew Jones 1 Final Exam 2 3 Mechanical Waves Waves and wave fronts: 4 Wave Motion 5 Two Kinds of Waves 6 Reflection of Waves When
More informationLand seismic sources
Seismic Sources HOW TO GENERATE SEISMIC WAVES? Exploration seismology mostly artificial sources à active technique Natural sources can also be used (e.g. earthquakes) usually for tectonic studies (passive
More informationTan K. Wang National Taiwan Ocean University, Keelung, Taiwan, R.O.C.
SEISMIC IMAGING IN THE OCEANS Tan K. Wang National Taiwan Ocean University, Keelung, Taiwan, R.O.C. Keywords: Converted wave, multi-channel seismic, ocean-bottom seismometer, travel-time inversion, virtual
More informationSEISMIC WAVE PROPAGATION IN FRACTURED CARBONATE ROCK
Proceedings of the Project Review, Geo-Mathematical Imaging Group (Purdue University, West Lafayette IN), Vol. 1 (2010) pp. 211-220. SEISMIC WAVE PROPAGATION IN FRACTURED CARBONATE ROCK WEIWEI LI AND LAURA
More informationAn Introduction to Geophysical Exploration
An Introduction to Geophysical Exploration Philip Kearey Department of Earth Sciences University of Bristol Michael Brooks Ty Newydd, City Near Cowbridge Vale of Glamorgan Ian Hill Department of Geology
More informationChapter 5: Application of Filters to Potential Field Gradient Tensor Data
Chapter 5: Filters 98 Chapter 5: Application of Filters to Potential Field Gradient Tensor Data 5.1 Introduction One of the objectives of this research is to investigate the use of spatial filters on potential
More informationFrom PZ summation to wavefield separation, mirror imaging and up-down deconvolution: the evolution of ocean-bottom seismic data processing
From PZ summation to wavefield separation, mirror imaging and up-down deconvolution: the evolution of ocean-bottom seismic data processing Sergio Grion, CGGVeritas Summary This paper discusses present
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 informationFloatSeis Technologies for Ultra-Deep Imaging Seismic Surveys
FloatSeis Technologies for Ultra-Deep Imaging Seismic Surveys 25 th January, 2018 Aleksandr Nikitin a.nikitin@gwl-geo.com Geology Without Limits Overview 2011-2016 GWL Acquired over 43000 km 2D seismic
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 informationP137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients
P137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients M. Ayzenberg (StatoilHydro), A. Aizenberg (Institute of Petroleum Geology and Geophysics),
More informationAPPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY
APPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY Timur TEZEL Supervisor: Takuo SHIBUTANI MEE07169 ABSTRACT In this study I tried to determine the shear wave velocity structure in the crust
More informationInvitation to Futterman inversion
Invitation to Futterman inversion Jon Claerbout ABSTRACT A constant Q earth model attenuates amplitude inversely with the number of wavelengths propagated, so the attenuation factor is e ω (z/v)/q. We
More informationFUNDAMENTALS OF SEISMIC EXPLORATION FOR HYDROCARBON
FUNDAMENTALS OF SEISMIC EXPLORATION FOR HYDROCARBON Instructor : Kumar Ramachandran 10 14 July 2017 Jakarta The course is aimed at teaching the physical concepts involved in the application of seismic
More informationApproximate- vs. full-hessian in FWI: 1D analytical and numerical experiments
Approximate- vs. full-hessian in FWI: 1D analytical and numerical experiments Raul Cova and Kris Innanen ABSTRACT Feasibility of using Full Waveform Inversion (FWI) to build velocity models has been increasing
More informationReview paper: Virtual sources and their responses, Part II: data-driven single-sided focusing
Geophysical Prospecting, 017, 65, 1430 1451 doi: 10.1111/1365-478.1495 Review paper: Virtual sources and their responses, Part II: data-driven single-sided focusing Kees Wapenaar 1, Jan Thorbecke 1, Joost
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 informationAdding Value with Broadband Seismic and Inversion in the Central North Sea Seagull Area
H2-2-1 Adding Value with Broadband Seismic and Inversion in the Central North Sea Seagull Area Marnix Vermaas, Andy Lind Apache North Sea, Aberdeen, UK Introduction The merits of modern broadband seismic
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 informationSeismic processing of numerical EM data John W. Neese* and Leon Thomsen, University of Houston
Seismic processing of numerical EM data John W. Neese* and Leon Thomsen, University of Houston Summary The traditional methods for acquiring and processing CSEM data are very different from those for seismic
More informationTu SRS3 02 Data Reconstruction and Denoising of Different Wavefield Components Using Green s Theorem
Tu SRS3 02 Data Reconstruction and Denoising of Different Wavefield Components Using Green s Theorem N. Kazemi* (University of Alberta) & A.C. Ramirez (Statoil) SUMMARY Multicomponent technology is likely
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 informationCompensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS
Compensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS SUMMARY Anelastic properties of the earth cause frequency dependent energy
More informationNumerical Study: Time-Reversed Reciprocal Method and Damage Detection Method for Weld Fracture
Chapter 4 Numerical Study: Time-Reversed Reciprocal Method and Damage Detection Method for Weld Fracture A numerical study is performed to gain insight into applying the proposed method of detecting high-frequency
More informationZ - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.
Z - Transform The z-transform is a very important tool in describing and analyzing digital systems. It offers the techniques for digital filter design and frequency analysis of digital signals. Definition
More informationTowards Modelling Elastic and Viscoelastic Seismic Wave Propagation in Boreholes
Towards Modelling Elastic and Viscoelastic Seismic Wave Propagation in Boreholes NA WANG, DONG SHI, BERND MILKEREIT Department of Physics, University of Toronto, Toronto, Canada M5S 1A7 Summary We are
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 informationDesigning Information Devices and Systems I Fall 2018 Lecture Notes Note Positioning Sytems: Trilateration and Correlation
EECS 6A Designing Information Devices and Systems I Fall 08 Lecture Notes Note. Positioning Sytems: Trilateration and Correlation In this note, we ll introduce two concepts that are critical in our positioning
More informationReflection Seismic Method
Reflection Seismic Method Data and Image sort orders; Seismic Impedance; -D field acquisition geometries; CMP binning and fold; Resolution, Stacking charts; Normal Moveout and correction for it; Stacking;
More informationNonlinear seismic imaging via reduced order model backprojection
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, Vladimir Druskin 2 and Mikhail Zaslavsky 2 University of Houston, 2 Schlumberger-Doll Research Center Mamonov, Druskin,
More informationThe z-transform Part 2
http://faculty.kfupm.edu.sa/ee/muqaibel/ The z-transform Part 2 Dr. Ali Hussein Muqaibel The material to be covered in this lecture is as follows: Properties of the z-transform Linearity Initial and final
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 informationReceiver. Johana Brokešová Charles University in Prague
Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate
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 information3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE
3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE 3.0 INTRODUCTION The purpose of this chapter is to introduce estimators shortly. More elaborated courses on System Identification, which are given
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 informationA particle/collision model of seismic data
Kris Innanen ABSTRACT A particle/collision model of seismic data It can occasionally be a valuable exercise to take a familiar phenomenon and observe it from a different point of view. It turns out, for
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 informationPART 1. Review of DSP. f (t)e iωt dt. F(ω) = f (t) = 1 2π. F(ω)e iωt dω. f (t) F (ω) The Fourier Transform. Fourier Transform.
PART 1 Review of DSP Mauricio Sacchi University of Alberta, Edmonton, AB, Canada The Fourier Transform F() = f (t) = 1 2π f (t)e it dt F()e it d Fourier Transform Inverse Transform f (t) F () Part 1 Review
More informationTim Carr - West Virginia University
Tim Carr - West Virginia University Understanding Seismic Data Resolution (Vertical and Horizontal) Common Depth Points (CDPs) Two way time (TWT) Time versus depth Interpretation of Reflectors 2 Able to
More informationA comparison of the imaging conditions and principles in depth migration algorithms
International Journal of Tomography & Statistics (IJTS), Fall 2006, Vol. 4, No. F06; Int. J. Tomogr. Stat.; 5-16 ISSN 0972-9976; Copyright 2006 by IJTS, ISDER A comparison of the imaging conditions and
More informationSTRATEGIES FOR DETECTING POOR COUPLING IN AN OBS EXPERIMENT
i STRATEGIES FOR DETECTING POOR COUPLING IN AN OBS EXPERIMENT A Thesis by FITRIX PRIMANTORO PUTRO Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements
More informationFinite difference elastic modeling of the topography and the weathering layer
Finite difference elastic modeling of the topography and the weathering layer Saul E. Guevara and Gary F. Margrave ABSTRACT Finite difference 2D elastic modeling is used to study characteristics of the
More informationAn equivalent fluid representation of a layered elastic seafloor for acoustic propagation modelling Matthew W Koessler (1)
An equivalent fluid representation of a layered elastic seafloor for acoustic propagation modelling Matthew W Koessler (1) ABSTRACT (1) Marshall Day Acoustics,6/448 Roberts Rd Subiaco, Australia Modelling
More informationEarth in 2-D, 3-D & 4-D
Earth in 2-D, 3-D & 4-D We will consider the scientific tools and techniques used to map surface features, reconstruct the layered structure of Earth, and interpret Earth history, including the origin
More information56 CHAPTER 3. POLYNOMIAL FUNCTIONS
56 CHAPTER 3. POLYNOMIAL FUNCTIONS Chapter 4 Rational functions and inequalities 4.1 Rational functions Textbook section 4.7 4.1.1 Basic rational functions and asymptotes As a first step towards understanding
More informationDesigning Information Devices and Systems I Fall 2018 Lecture Notes Note Introduction to Linear Algebra the EECS Way
EECS 16A Designing Information Devices and Systems I Fall 018 Lecture Notes Note 1 1.1 Introduction to Linear Algebra the EECS Way In this note, we will teach the basics of linear algebra and relate it
More informationReceiver Function (RF) Estimation Using Short Period Seismological Data
Receiver Function (RF) Estimation Using Short Period Seismological Data Shantanu Pandey Department of Earth Sciences, Kurukshetra University, Kurukshetra-136119 Abstract The advent of the receiver function
More informationCommunication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi
Communication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture - 3 Brief Review of Signals and Systems My subject for today s discussion
More informationRough sea estimation for phase-shift de-ghosting Sergio Grion*, Rob Telling and Seb Holland, Dolphin Geophysical
Rough sea estimation for phase-shift de-ghosting Sergio Grion*, Rob Telling and Seb Holland, Dolphin Geophysical Summary This paper discusses rough-sea de-ghosting for variabledepth streamer data. The
More informationTh Guided Waves - Inversion and Attenuation
Th-01-08 Guided Waves - Inversion and Attenuation D. Boiero* (WesternGeco), C. Strobbia (WesternGeco), L. Velasco (WesternGeco) & P. Vermeer (WesternGeco) SUMMARY Guided waves contain significant information
More information11th Biennial International Conference & Exposition. Keywords Sub-basalt imaging, Staggered grid; Elastic finite-difference, Full-waveform modeling.
Sub-basalt imaging using full-wave elastic finite-difference modeling: A synthetic study in the Deccan basalt covered region of India. Karabi Talukdar* and Laxmidhar Behera, CSIR-National Geophysical Research
More informationInversion for Geoacoustic Model Parameters in Range-Dependent Shallow Water Environments
DISTRIBUTION STATEMENT A: Approved for public release; distribution is unlimited. Inversion for Geoacoustic Model Parameters in Range-Dependent Shallow Water Environments N. Ross Chapman School of Earth
More informationarxiv: v1 [math.na] 1 Apr 2015
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, University of Houston; Vladimir Druskin and Mikhail Zaslavsky, Schlumberger arxiv:1504.00094v1 [math.na] 1 Apr 2015
More informationDesigning Information Devices and Systems I Fall 2018 Lecture Notes Note Positioning Sytems: Trilateration and Correlation
EECS 6A Designing Information Devices and Systems I Fall 08 Lecture Notes Note. Positioning Sytems: Trilateration and Correlation In this note, we ll introduce two concepts that are critical in our positioning
More informationTu 23 A12 Multi-frequency Seafloor Characterization Using Seismic Sources of Opportunity
Tu 23 A12 Multi-frequency Seafloor Characterization Using Seismic Sources of Opportunity M.N. Banda* (University of Bath/Seiche Ltd), Ph. Blondel (University of Bath), M. Burnett (Seiche Ltd), R. Wyatt
More informationStochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore
Stochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore Lecture No. # 33 Probabilistic methods in earthquake engineering-2 So, we have
More informationLet us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a).
7.1. Low-Coherence Interferometry (LCI) Let us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a). The light is split by the beam splitter (BS) and
More information2. Intersection Multiplicities
2. Intersection Multiplicities 11 2. Intersection Multiplicities Let us start our study of curves by introducing the concept of intersection multiplicity, which will be central throughout these notes.
More informationGG710 Remote Sensing in Submarine Environments Sidescan Sonar
GG710 Remote Sensing in Submarine Environments Sidescan Sonar Harold Edgerton, a professor of electrical engineering at the Massachusetts Institute of Technology, developed sidescan sonar technology for
More informationBorehole Seismic Monitoring of Injected CO 2 at the Frio Site
Borehole Seismic Monitoring of Injected CO 2 at the Frio Site * Daley, T M (tmdaley@lbl.gov), Lawrence Berkeley National Lab., 1 Cyclotron Rd, Berkeley, CA 94720 Myer, L (lrmyer@lbl.gov), Lawrence Berkeley
More informationCHALLENGE! (0) = 5. Construct a polynomial with the following behavior at x = 0:
TAYLOR SERIES Construct a polynomial with the following behavior at x = 0: CHALLENGE! P( x) = a + ax+ ax + ax + ax 2 3 4 0 1 2 3 4 P(0) = 1 P (0) = 2 P (0) = 3 P (0) = 4 P (4) (0) = 5 Sounds hard right?
More information3. Magnetic Methods / 62
Contents Preface to the Second Edition / xv Excerpts from Preface to the FirstEdition / xvii Mathematical Conventions / xix 1. Introduction / 1 Reference / 5 2. Gravity Methods / 6 2. I. Introduction /
More informationDE128. Seismic Inversion Techniques. H.H. Sheik Sultan Tower (0) Floor Corniche Street Abu Dhabi U.A.E
DE128 Seismic Inversion Techniques H.H. Sheik Sultan Tower (0) Floor Corniche Street Abu Dhabi U.A.E www.ictd.ae ictd@ictd.ae Course Introduction: Seismic inversion is now commonly used on post-stack and
More informationEffects of Surface Geology on Seismic Motion
4 th IASPEI / IAEE International Symposium: Effects of Surface Geology on Seismic Motion August 23 26, 2011 University of California Santa Barbara VELOCITY STRUCTURE INVERSIONS FROM HORIZONTAL TO VERTICAL
More informationGeophysical Applications Seismic Reflection Processing
Seismic reflection data are routinely acquired for multiple purposes such as exploration, mining, or engineering problems. The seismic data are generally acquired in shot-gathers, i.e. the data is sorted
More informationRAYFRACT IN MARINE SURVEYS. The data provided by the contractor needed total re-picking and reinterpretation.
RAYFRACT IN MARINE SURVEYS As part of a geotechnical assessment and feasibility planning of channel improvement a shallow marine seismic refraction survey was undertaken. The data was initially processed
More informationAdvanced Optical Communications Prof. R. K. Shevgaonkar Department of Electrical Engineering Indian Institute of Technology, Bombay
Advanced Optical Communications Prof. R. K. Shevgaonkar Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture No. # 15 Laser - I In the last lecture, we discussed various
More information