660 likes | 875 Views
Making CMP’s. From chapter 16 “Elements of 3D Seismology” by Chris Liner. Outline. Convolution and Deconvolution Normal Moveout Dip Moveout Stacking. Outline. Convolution and Deconvolution Normal Moveout Dip Moveout Stacking. Convolution means several things:.
E N D
Making CMP’s From chapter 16 “Elements of 3D Seismology” by Chris Liner
Outline • Convolution and Deconvolution • Normal Moveout • Dip Moveout • Stacking
Outline • Convolution and Deconvolution • Normal Moveout • Dip Moveout • Stacking
Convolution means several things: • IS multiplication of a polynomial series • IS a mathematical process • IS filtering
Convolution means several things: • IS multiplication of a polynomial series A * B = C E.g., A= 0.25 + 0.5 -0.25 0.75]; B = [1 2 -0.5]; C = [0.2500 1.0000 0.6250 0 1.6250 -0.3750]
Convolutional Model for the Earth output input Reflections in the earth are viewed as equivalent to a convolution process between the earth and the input seismic wavelet.
Convolutional Model for the Earth output input SOURCE * Reflection Coefficient = DATA (input) (earth) (output) where * stands for convolution
Convolutional Model for the Earth SOURCE * Reflection Coefficient = DATA (input) (earth) (output) where * stands for convolution (MORE REALISTIC) SOURCE * Reflection Coefficient + noise = DATA (input) (earth) (output) s(t) * e(t) + n(t) = d(t)
Convolution in theTIME domainis equivalent toMULTIPLICATION in theFREQUENCYdomain s(t) * e(t) + n(t) = d(t) FFT FFT FFT s(f,phase) x e(f,phase) + n(f,phase) = d(f,phase) Inverse FFT d(t)
CONVOLUTION as a mathematical operator signal has 3 terms (j=3) -1 2 -1/2 earth Reflection Coefficient has 4 terms (k=4) 1/4 1/4 time 1/2 z 1/2 -1/4 3/4 -1/4 3/4 Reflection Coefficients with depth (m)
0 0 0 -1/2 2 1 0 0 0 0 0 0 0 0 0 0 x x x x x = = = = = 0 0 0 1/4 1/2 -1/4 3/4 0 0 0 +
0 0 0 -1/2 2 -1 0 0 0 0 0 0 0 0 0 0 0 x x x x x = = = = = = 0 0 0 1/4 1/2 -1/4 3/4 0 0 0 +
0 0 0 -1/2 2 1 0 0 0 0 0 0 0 0 0 0 0 0 x x x x x x x = = = = = = = 0 0 0 1/4 1/2 -1/4 3/4 0 0 0 +
0 0 0 -1/2 2 1 0 0 0 0 0 0 0 1/4 0 0 0 0 1/4 x x x x x x x x = = = = = = = = 0 0 0 1/4 1/2 -1/4 3/4 0 0 0 +
0 0 0 1/2 1/2 0 0 0 0 1 0 0 0 -1/2 2 1 0 0 0 0 x x x x x x x x x = = = = = = = = = 0 0 0 1/4 1/2 -1/4 3/4 0 0 0 +
0 0 0 -1/8 1 -1/4 0 0 0 0 5/8 x x x x x x x x x x = = = = = = = = = = 0 0 0 -1/2 2 1 0 0 0 0 0 0 0 1/4 1/2 -1/4 3/4 0 0 0 +
0 0 0 0 -1/4 -1/2 3/4 0 0 0 0 x x x x x x x x x x = = = = = = = = = = 0 0 0 1/4 1/2 -1/4 3/4 0 0 0 0 0 0 -1/2 2 1 0 0 0 0 +
0 0 0 1/8 1 1/2 0 0 0 1 5/8 x x x x x x x x = = = = = = = = 0 0 0 1/4 1/2 -1/4 3/4 0 0 0 + 0 0 0 -1/2 2 1 0 0 0 0
0 0 0 -3/8 0 0 0 -3/8 x x x x x x x = = = = = = = 0 0 0 1/4 1/2 -1/4 3/4 0 0 0 + 0 0 0 -1/2 2 1 0 0 0 0
0 0 0 0 0 0 0 x x x x x x = = = = = = 0 0 0 1/4 1/2 -1/4 3/4 0 0 0 + 0 0 0 -1 2 -1/2 0 0 0 0
MATLAB %convolution a = [0.25 0.5 -0.25 0.75]; b = [1 2 -0.5]; c = conv(a,b) d = deconv(c,a) c = 0.2500 1.0000 0.6250 0 1.6250 -0.3750 matlab
Outline • Convolution and Deconvolution • Normal Moveout • Dip Moveout • Stacking
Normal Moveout Hyperbola: x T
Normal Moveout x T “Overcorrected” Normal Moveout is too large Chosenvelocity for NMO is too (a) large (b) small
Normal Moveout x T “Overcorrected” Normal Moveout is too large Chosenvelocity for NMO is too (a)large (b) small
Normal Moveout x T “Under corrected” Normal Moveout is too small Chosenvelocity for NMO is (a) too large (b) too small
Normal Moveout x T “Under corrected” Normal Moveout is too small Chosenvelocity for NMO is (a) too large (b) too small
Vinterval from Vrms Dix, 1955
Vrms V1 V2 Vrms < Vinterval V3
Multiples and Primaries x M1 T M2
Conventional NMO before stacking x M1 NMO correction V=V(depth) e.g., V=mz + B T M2 “Properly corrected” Normal Moveout is just right Chosenvelocity for NMO is correct
Over-correction (e.g. 80% Vnmo) x x M1 M1 NMO correction V=V(depth) e.g., V=0.8(mz + B) T T M2 M2
f-k filtering before stacking (Ryu) x x M1 NMO correction V=V(depth) e.g., V=0.8(mz + B) T T M2 M2
Correct back to 100% NMO x x M1 M1 NMO correction V=V(depth) e.g., V=(mz + B) T T M2 M2
Outline • Convolution and Deconvolution • Normal Moveout • Dip Moveout • Stacking
Outline • Convolution and Deconvolution • Normal Moveout • Dip Moveout • Stacking
Dip Moveout (DMO) (Ch. 19; p.365-375) How do we move out a dipping reflector in our data set? m Offset (m) TWTT (s) z
Dip Moveout • A dipping reflector: • appears to be faster • its apex may not be centered Offset (m) For a dipping reflector: Vapparent = V/cos dip TWTT (s) e.g., V=2600 m/s Dip=45 degrees, Vapparent = 3675m/s
CONFLICTING DIPS Different dips CAN NOT be NMO’d correctly at the same time Offset (m) 3675 m/s TWTT (s) 2600 m/s Vrms for dipping reflector too low & overcorrects Vrms for dipping reflector is correct but undercorrects horizontal reflector
DMO Theoretical Background (Yilmaz, p.335) (Levin,1971) is layer dip “NMO”
DMO Theoretical Background (Yilmaz, p.335) (Levin,1971) “DMO”
Three properties of DMO “DMO” “NMO” (1) DMO effect at 0 offset = ? (2) As the dip increases DMO (a) increases (B) decreases (3) As velocity increases DMO (a) increases (B) decreases
Three properties of DMO “DMO” “NMO” (1) DMO effect at 0 offset = 0 (2) As the dip increases DMO (a) increases(B) decreases (3) As velocity increases DMO (a) increases(B) decreases
Application of DMO aka “Pre-stack partical migration” • (1) DMOafter NMO (applied to CDP/CMP data) • but before stacking • DMO is applied to Common-Offset Data • Is equivalent to migration of stacked data • Works best if velocity is constant
DMO Implementation before stack -I Offset (m) (1) NMO using background Vrms TWTT (s)