Chao Ma!! San Antonio, TX! 2013 M-OSRP Annual Meeting! May 2, 2013!

Size: px
Start display at page:

Download "Chao Ma!! San Antonio, TX! 2013 M-OSRP Annual Meeting! May 2, 2013!"

Transcription

1 One-dimensional analytic analysis of the effects of! treating s in the data as subevents! in the leading-order ISS multiple-removal algorithms:! comparison between free-surface and cases Chao Ma!! San Antonio, TX! 2013 M-OSRP Annual Meeting! May 2, 2013! 1

2 Key points! ISS leading-order attenuation:!! Strengths (data driven):! Predicts all first-order s at all depths at once;! The prediction has correct time and an approximate, wellunderstood amplitude of first-order s;! A limitation & Solution:! A new higher-order attenuation algorithm is developed, retaining the strengths of the current algorithm and addressing one of its limitations.!! 2

3 Key points & Outline! ISS Leading-order attenuation algorithm! A limitation and solution (higherorder modification)! 3

4 Key points & Outline! ISS Leading-order attenuation algorithm! A limitation and solution (higherorder modification)! Two examples to help understanding! 4

5 Key points & Outline! ISS Leading-order attenuation algorithm! A limitation and solution (higherorder modification)! Simplest circumstance (a three reflector example) the limitation can occur! Two examples to help understanding! 5

6 Key points & Outline! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and solution (higherorder modification)! Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification and When this limitation can be serious and needs higher-order modification! 6

7 Key points & Outline! Introduction! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and solution (higherorder modification)! Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification and When this limitation can be serious and needs higher-order modification! Summary! 7

8 Introduction! M-OSRP is concerned with the seismic data processing.! Green s Theorem! 1. Wavefield separation! 2. Wavelet estimation! 3. Deghosting! Inverse Scattering Series! (ISS)! 4. Free-surface multiple removal! 5. Internal-multiple removal! 6. Imaging! 7. Inversion! 8

9 Introduction! Green s Theorem! Data without reference and ghosts! Inverse scattering series(iss)! ISS free surface multiple removal ISS internal multiple removal ISS imaging ISS inversion! attenuation ISS leading-order attenuation subseries for firstorder! 9

10 Introduction! Green s Theorem! Data without reference and ghosts! Inverse scattering series(iss)! ISS free surface multiple removal ISS internal multiple removal ISS imaging ISS inversion! attenuation ISS leading-order attenuation subseries for firstorder! 10

11 Introduction! Green s Theorem! Data without reference and ghosts! Inverse scattering series(iss)! ISS free surface multiple removal ISS internal multiple removal ISS imaging ISS inversion! attenuation b IM = b 1 + b 3 + b 5 +. ISS leading-order attenuation subseries for firstorder! 11

12 Introduction! Green s Theorem! Data without reference and ghosts! Inverse scattering series(iss)! ISS free surface multiple removal ISS internal multiple removal ISS imaging ISS inversion! attenuation Stolt migration! of input data! (P+IM)! b IM = b 1 + b 3 + b 5 +. Leading-order prediction of secondorder s! ISS leading-order attenuation subseries for firstorder! D(t) D(ω ) D( 2ω ) c 0 b 1 (k) b 1 (z) Leading-order prediction of first-order s! k = 2ω ; z c 0t c

13 Introduction! Green s Theorem! Data without reference and ghosts! Inverse scattering series(iss)! ISS free surface multiple removal ISS internal multiple removal ISS imaging ISS inversion attenuates internal multiples of all orders! ISS! leading-order! attenuation Stolt migration! of input data! (P+IM)! b IM = b 1 + b 3 + b 5 +. Leading-order prediction of secondorder s! ISS leading-order attenuation subseries for firstorder! Leading-order prediction of first-order s! 13

14 Introduction! Green s Theorem! Data without reference and ghosts! Inverse scattering series(iss)! ISS free surface multiple removal ISS internal multiple removal ISS imaging ISS inversion attenuates internal multiples of all orders! ISS! leading-order! attenuation Stolt migration! of input data! (P+IM)! b IM = b 1 + b 3 + b 5 +. Leading-order prediction of secondorder s! attenuates internal multiples of first-order! ISS leading-order attenuation subseries for firstorder! Leading-order prediction of first-order s! b 1 + b 3 14

15 Introduction! Green s Theorem! Data without reference and ghosts! Inverse scattering series(iss)! ISS free surface multiple removal ISS internal multiple removal ISS imaging ISS inversion attenuates internal multiples of all orders! ISS! leading-order! attenuation Stolt migration! of input data! (P+IM)! b IM = b 1 + b 3 + b 5 +. Leading-order prediction of secondorder s! attenuates internal multiples of first-order! ISS leading-order attenuation subseries for firstorder! Leading-order prediction of first-order s! b 1 + b 3 15

16 Key points & Outline! Introduction! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and its higher-order modification! Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification and When this limitation can be serious and needs higher-order modification! Summary! 16

17 The ISS leading-order attenuation! The leading-order prediction of the first-order s, in 1D, is! (Araújo et al, 1994; Weglein et al., 1997): b 3 (k) = dz 1 e ikz 1 b (z ) dz e ikz 2 b (z ) dz e ikz 3 b (z ) 1 3 z 1 ε z 2 +ε 17

18 The ISS leading-order attenuation! The leading-order prediction of the first-order s, in 1D, is! (Araújo et al, 1994; Weglein et al., 1997): b 3 (k) = dz 1 e ikz 1 b (z ) dz e ikz 2 b (z ) dz e ikz 3 b (z ) 1 3 z 1 ε z 2 +ε Data after the removal of reference wave, wavelet, ghosts, free-surface multiples.! D(t) D(ω ) D( 2ω c 0 ) b 1 (k) b 1 (z) k = 2ω ; z c t 0 c

19 The ISS leading-order attenuation! The leading-order prediction of the first-order s, in 1D, is! (Araújo et al, 1994; Weglein et al., 1997): Primaries + Internal multiples! Primaries + Internal multiples! b 3 (k) = dz 1 e ikz 1 b (z ) dz e ikz 2 b (z ) dz e ikz 3 b (z ) 1 3 z 1 ε z 2 +ε Primaries + Internal multiples! Data after the removal of reference wave, wavelet, ghosts, free-surface multiples.! D(t) D(ω ) D( 2ω c 0 ) b 1 (k) b 1 (z) k = 2ω ; z c t 0 c

20 The ISS leading-order attenuation! The leading-order prediction of the first-order s, in 1D, is! (Araújo et al, 1994; Weglein et al., 1997): Primaries + Internal multiples! Primaries + Internal multiples! b 3 (k) = dz 1 e ikz 1 b (z ) dz e ikz 2 b (z ) dz e ikz 3 b (z ) 1 3 Air-water/! Primaries + Internal multiples! Air-land! A! D! C! G! z 1 ε z 2 +ε E! F! ABC (+)! DEC (-)! DFG (+)! B! First-order prediction through three primary subevents 20

21 The ISS leading-order attenuation! The leading-order prediction of the first-order s, in 1D, is! (Araújo et al, 1994; Weglein et al., 1997): Primaries + Internal multiples! Primaries + Internal multiples! b 3 (k) = dz 1 e ikz 1 b (z ) dz e ikz 2 b (z ) dz e ikz 3 b (z ) 1 3 Air-water/! Primaries + Internal multiples! Air-land! A! D! C! G! z 1 ε z 2 +ε E! F! ABC (+)! DEC (-)! DFG (+)! B! First-order prediction through three primary subevents 21

22 What leading-order term can achieve (among others)! ISS leading-order attenuation subseries for firstorder.! b 1 (k) + b 3 (k) 1. Predicts all first-order s at all depths at once;! 2. The correct time and an approximate, well-understood amplitude of first-order s;! Data driven; avoids assumptions that many other methods make, e.g., subsurface information.! 22

23 The ISS leading-order attenuation! The leading-order prediction of the first-order s, in 1D, is! (Araújo et al, 1994; Weglein et al., 1997): Primaries + Internal multiples! Primaries + Internal multiples! b 3 (k) = dz 1 e ikz 1 b (z ) dz e ikz 2 b (z ) dz e ikz 3 b (z ) 1 3 Air-water/! Primaries + Internal multiples! Air-land! A! D! C! G! z 1 ε z 2 +ε E! F! ABC (+)! DEC (-)! DFG (+)! B! First-order prediction through three primary subevents 23

24 Possible subevents combinations! b1 = P + I b 3 = (P + I)(P + I)(P + I) = PPP + PPI + IPP + IPI + III + PIP + PII + IIP Where P stands for primaries,! I stands for s. First-order s! 24

25 Possible subevents combinations! b1 = P + I (Two reflector example)! Higher-order s! e.g., second-order, third-order *! b 3 = (P + I)(P + I)(P + I) = PPP + PPI + IPP + IPI + III + PIP + PII + IIP Where P stands for primaries,! I stands for s. First-order s! * Zhang and Shaw (2010)! 25

26 Prediction of a second-order! (an first-order acts as one subevent)! 1! 2! Pimary_2! Pimary_1! A first-order! A second-order internal multiple! 26

27 ISS attenuation! attenuates internal multiples of all orders! attenuates internal multiples of first-order! ISS internal multiple removal ISS leading-order attenuation ISS leading-order attenuation subseries for firstorder! Stolt migration! of input data! (P+IM)! b IM = b 1 + b 3 + b 5 +. Leading-order prediction of first-order internal multiples! 27

28 ISS attenuation! attenuates internal multiples of all orders! attenuates internal multiples of first-order! ISS internal multiple removal ISS leading-order attenuation ISS leading-order attenuation subseries for firstorder! Stolt migration! of input data! (P+IM)! Leading-order prediction of second-order internal multiples! b IM = b 1 + b 3 + b 5 +. Leading-order prediction of first-order internal multiples! + Higher-order internal multiples! 28

29 Key points & Outline! Introduction! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and solution (higherorder modification)! 1. Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification When this limitation can be serious and needs higher-order modification! Free-surface multiple removal! Summary! 29

30 1-D free-surface multiple elimination! Green s Theorem! Data without reference and ghosts! ISS free surface multiple removal Data without reference and ghosts! Prediction of second-order free-surface multiple! R(ω ) = R FS (ω ) + R 2 FS(ω ) + R 3 FS(ω ) + R 4 FS(ω ) + Prediction of first-order free-surface multiple! (-1)! (-1)! (-1)! Air-water! 1! R 1! R 1! R 1! R 1! R 1! R 1! Water bottom! Pimary_1! A first-order free-surface multiple! A second-order free-surface multiple! 30

31 1-D free-surface multiple elimination! Green s Theorem!! term! Data without reference and ghosts! data! ISS free surface multiple removal Primary! Data without reference and ghosts! R(ω ) = R FS (ω ) + R 2 FS(ω ) + R 3 FS(ω ) + R 4 FS(ω ) + First-order! Free-surface multiples! Prediction of first-order free-surface multiple! Prediction of second-order free-surface multiple! Second-order! Free-surface multiples! R FS! 1! 1! 1! R 2 FS! 0! -1! -2! R 3 FS! 0! 0! 1! 31

32 1-D free-surface multiple elimination! 1! R 1! Air-water! Water bottom! Pimary_1! Primary_1! Actual amplitude: (R 1 )(-1)(R 1 )! Predicted amplitude: [R 1 ][R 1 ]! A first-order free-surface multiple! 1! R 1! Actual amplitude: (R 1 )(-1)(R 1 )(-1)(R 1 )! Predicted amplitude: [R 1 (-1)(R 1 )][R 1 ]! Air-water! Water bottom! 32

33 1-D free-surface multiple elimination! Green s Theorem!! term! Data without reference and ghosts! data! ISS free surface multiple removal Primary! Data without reference and ghosts! R(ω ) = R FS (ω ) + R 2 FS(ω ) + R 3 FS(ω ) + R 4 FS(ω ) + First-order! Free-surface multiples! Prediction of first-order free-surface multiple! Prediction of second-order free-surface multiple! Second-order! Free-surface multiples! R FS! 1! 1! 1! R 2 FS! 0! -1! -2! R 3 FS! 0! 0! 1! 33

34 Key points & Outline! Introduction! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and solution (higherorder modification)! Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification When this limitation can be serious and needs higher-order modification! Free-surface multiple removal! To remove second-order free-surface multiples, more than one terms are needed (collective works are needed).! Summary! 34

35 Key points & Outline! Introduction! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and solution (higherorder modification)! Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification When this limitation can be serious and needs higher-order modification! Free-surface multiple removal! Internal multiple attenuation (two reflector example)! Summary! 35

36 ISS attenuation! attenuates internal multiples of all orders! attenuates internal multiples of first-order! ISS internal multiple removal ISS leading-order attenuation ISS leading-order attenuation subseries for firstorder! Stolt migration! Leading-order prediction of second-order internal multiples! b IM = b 1 + b 3 + b 5 +. Leading-order prediction of first-order internal multiples! + Higher-order internal multiples! 36

37 1-D internal-multiple attenuation (two-reflector case)! 1! 2! R 1! T 01! R 2! T 10! Pimary_1! Pimary_2! A first-order! A second-order! 37

38 1-D internal-multiple attenuation (two-reflector case)! Data without reference, ghosts and free-surface multiples! ISS free surface multiple removal ISS leading-order attenuation Stolt migration! Leading-order prediction of secondorder s! b IM (k) = b 1 (k) + b 3 (k) + b 5 (k) +. Leading-order prediction of first-order s!! term! data! Primary! First-order! internalmultiples! Second-order! Internal multiples! b 1! 1! 1! 1! b 3! 0! -1(T 01 T 10 )! -2(T 01 T 10 )+ (T 01 T 10 R 1 ) 2! b 5! 0! 0! 1(T 01 T 10 ) 2! 38

39 1-D internal-multiple attenuation (two-reflector case)! T 01! 1! R 2! 2! T 10! R 1! Pimary_2! Pimary_1! Pimary_2! A first-order! Actual amplitude: T 01 R 2 (-R 1 )R 2 T 12! Predicted amplitude: [T 01 R 2 T 10 ][R 1 ][T 01 R 2 T 10 ]=(-1T 01 T 10 )[T 01 R 2 (-R 1 )R 2 T 12 ]! 39

40 1-D internal-multiple attenuation (two-reflector case)! 1! 2! A second-order! Pimary_1! Pimary_2! A second-order! Actual amplitude: T 01 R 2 (-R 1 )R 2 (-R 1 )R 2 T 10! Predicted amplitude: [T 01 R 2 (-R 1 )R 2 T 10 ][R 1 ][T 01 R 2 T 10 ]! =(-1T 01 T 10 )[T 01 R 2 (-R 1 )R 2 (-R 1 )R 2 T 10 ]! 40

41 1-D internal-multiple attenuation (two-reflector case)! 1! 2! Pimary_2! Pimary_1! Pimary_2! Pimary_1! Pimary_2! A second-order! Actual amplitude: T 01 R 2 (-R 1 )R 2 (-R 1 )R 2 T 10! Predicted amplitude: [T 01 R 2 T 10 ][R 1 ][T 01 R 2 T 10 ] [R 1 ][T 01 R 2 T 10 ]! =(T 01 T 10 ) 2 [T 01 R 2 (-R 1 )R 2 (-R 1 )R 2 T 10 ]! 41

42 Key points & Outline! Introduction! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and solution (higherorder modification)! Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification When this limitation can be serious and needs higher-order modification! To attenuate second-order s, more than one terms are needed (collective works are needed).! Internal multiple attenuation (two reflector example)! Summary! 42

43 Key points & Outline! Introduction! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and solution (higherorder modification)! Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification When this limitation can be serious and needs higher-order modification! Terms within each subseries work collectively to achieve the corresponding task! (e.g., free-surface multiple removal and attenuation)! Summary! 43

44 Key points & Outline! Introduction! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and solution (higherorder modification)! Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification When this limitation can be serious and needs higher-order modification! Summary! 44

45 A limitation of the leading-order algorithm! b1 = P + I (Two reflector example)! Higher-order s! e.g., second-order, third-order *! b 3 = (P + I)(P + I)(P + I) = PPP + PPI + IPP + IPI + III + PIP + PII + IIP Where P stands for primaries,! I stands for s. First-order s! * Zhang and Shaw (2010)! 45

46 A limitation of the leading-order algorithm! b1 = P + I (Two reflector example)! Higher-order s! e.g., second-order, third-order! b 3 = (P + I)(P + I)(P + I) = PPP + PPI + IPP + IPI + III + PIP + PII + IIP Where P stands for primaries,! I stands for s. First-order s! (three reflector example)! Higher-order s! + Spurious events! 46

47 A three reflector model! 1! 2! R 1! T 01! R 2! T 10! T 12! T 21! 3! Pimary_1! Pimary_2! A first-order! R 3! Pimary_3! 47

48 A limitation of the leading-order algorithm! P + I! P + I! P + I! z ikz ε ikz ikz 3 = 1 1 z + ε 1 b ( k ) dze b ( z ) dze b ( z ) dz e b ( z ) 48

49 A limitation of the leading-order algorithm! P I! P z ikz ε ikz ikz 3 = 1 1 z + ε 1 b ( k ) dze b ( z ) dze b ( z ) dz e b ( z ) P I P Primary-Internal multiple-primary! 49

50 A limitation of the leading-order algorithm! 1! 2! 3! Pirmary_3! An internal mutiple_212! Pirmary_3! A spurious event! Assumption: t 3 >(2t 2 -t 1 )! 50

51 Key points & Outline! Introduction! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and solution (higherorder modification)! Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification When this limitation can be serious and needs higher-order modification! Terms within each subseries work collectively to achieve the corresponding task! (e.g., free-surface multiple removal and attenuation)! Summary! 51

52 Key points & Outline! Introduction! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and solution (higherorder modification)! Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification and when this limitation can be serious and needs higher-order modification! Summary! 52

53 Higher-order terms introduced to deal with spurious prediction! P I! P z ikz ε ikz ikz 3 = 1 1 z + ε 1 b ( k ) dze b ( z ) dze b ( z ) dz e b ( z ) P I P 53

54 Higher-order terms introduced to deal with spurious prediction! P I! P z ikz ε ikz ikz 3 = 1 1 z + ε 1 b ( k ) dze b ( z ) dze b ( z ) dz e b ( z ) b 5 PIP (k) = dze ikz b 1 (z) d z e ik z b ( z ) ik d z e z b 3 1 ( z ) z ε z +ε P I pred P 54

55 Higher-order terms introduced to deal with spurious prediction! P I! P z ikz ε ikz ikz 3 = 1 1 z + ε 1 b ( k ) dze b ( z ) dze b ( z ) dz e b ( z ) b 5 PIP (k) = dze ikz b 1 (z) d z e ik z b ( z ) ik d z e z b 3 1 ( z ) z ε z +ε PIP b 5 is derived from!(g d 0V 1 G d 0V 3 G d 0V 1 G d 0) m.s. P I pred P 55

56 When the limitation is significant! When is this spurious prediction important?!! If N=2,!! No spurious prediction;! If N>=3 & Amplitude of primary Amplitude of s,! e.g., Off-shore Brazil, Middle East, Western Canada,!! Higher-order modification should be included.! 56

57 Higher-order terms introduced to deal with spurious prediction! Model parameters! V 1 = 1500m / s; ρ 1 = 1.0g / cm 3 V 2 = 1700m / s; ρ 2 = 1.8g / cm 3 Z 1 = 500m Z 2 = 900m V 3 = 1750m / s; ρ 3 = 1.0g / cm 3 Z 3 = 1530m V 4 = 5000m / s; ρ 4 = 4.0g / cm 3 57

58 Higher-order terms introduced to deal with spurious prediction! Traces P 1! P 1! P 2! P 3! 1 1! Z 1 = 500m Time (s) 2 P 2! I 212! 2! Z 2 = 900m x P 3! 3! Z 3 = 1530m 4 Input data! 58

59 Higher-order terms introduced to deal with spurious prediction! Traces P 1! I 212! 1 1! Z 1 = 500m Time (s) 2 P 2! I 212! 2! Z 2 = 900m x P 3! 3! Z 3 = 1530m 4 Input data! 59

60 Higher-order terms introduced to deal with spurious prediction! Traces Trace P 1! P 2! Time (s) 2 I 212! Time(s) 2 I 2-1-2! x P 3! I ! 4 4 Spurious event! (P 3 -I 212 -P 3 )! Input data! Leading-order prediction b 3! 60

61 Higher-order terms introduced to deal with spurious prediction! Trace Trace Time(s) 2 3 I 212! I 212! Time(s) Spurious event! (P 3 -I 212 -P 3 )! I 21212! I 312! 4 Higher-order modification! Leading-order prediction b 3! Higher-order modification b 5 PIP! 61

62 Key points & Outline! Introduction! ISS Leading-order attenuation algorithm! Two examples to help understanding! A limitation and solution (higherorder modification)! Simplest circumstance (a three reflector example) the limitation can occur! Higher-order modification and when this limitation can be serious and needs higher-order modification! Free-surface multiple removal! Internal multiple attenuation (two reflector example)! Summary! 62

63 Summary! 1. We analyze the effects of treating s as subevents in the ISS leading-order attenuation algorithm;! 2. When the current algorithm can produce spurious prediction;! 3. Higher-order modifications are identified within the ISS framework to accommodate the limitation of the current leading-order algorithm.! Number of reflectors! (N)! ISS attenuation algorithm! N=2! b 1 +b 3! N=3! b 1 +b 3 +b 5 PIP! N>=4 (Hong Liang s talk)! 63

64 Key points! ISS leading-order attenuation:!! Strengths (data driven):! Predicts all first-order s at all depths at once;! The prediction has correct time and an approximate, wellunderstood amplitude of first-order s;! A limitation & Solution:! A new higher-order attenuation algorithm is developed, retaining the strengths of the current algorithm and addressing one of its limitations.!! 64

65 65

66 Reference!! Araujo, F. V., A. B. Weglein, P. M. Carvalho, and R. H. Stolt, 1994, Inverse scattering series for multiple attenuation: An example with surface and s: SEG Technical Program Expanded Abstracts, !! Weglein, A. B., F. A. Gasparotto, P. M. Carvalho, and R. H. Stolt, 1997, An inversescattering series method for attenuating multiples in seismic reflection data:! Geophysics, 62, !! Zhang, H., and S. Shaw, 2010, 1-d analytical analysis of higher order internal multiples predicted via the inverse scattering series based algorithm: SEG Expanded Abstracts, 29, !!! 66

SUMMARY THE ISS INTERNAL-MULTIPLE-ATTENUATION AL- GORITHM

SUMMARY THE ISS INTERNAL-MULTIPLE-ATTENUATION AL- GORITHM Comparing the new Inverse Scattering Series (ISS) internal-multiple-elimination algorithm and the industry-standard ISS internal-multiple-attenuation algorithm plus adaptive subtraction when primaries

More information

1D internal multiple prediction in a multidimensional world: errors and recommendations

1D internal multiple prediction in a multidimensional world: errors and recommendations 1D internal multiple prediction 1D internal multiple prediction in a multidimensional world: errors and recommendations Pan Pan and Kris Innanen ABSTRACT Internal multiples are more difficult to estimate

More information

1.5D internal multiple prediction on physical modeling data

1.5D internal multiple prediction on physical modeling data 1.5D internal multiple prediction on physical modeling data Pan Pan*, Kris Innanen and Joe Wong CREWES, University of Calgary Summary Multiple attenuation is a key aspect of seismic data processing, with

More information

Arthur B. Weglein. versus use multiples. The premise behind that versus phrasing speaks to a competing and adversarial relationship.

Arthur B. Weglein. versus use multiples. The premise behind that versus phrasing speaks to a competing and adversarial relationship. Direct and indirect inversion and a new and comprehensive perspective on the role of primaries and multiples in seismic data processing for structure determination and amplitude analysis Arthur B. Weglein

More information

1.5D internal multiple prediction on physical modeling data

1.5D internal multiple prediction on physical modeling data Internal multiple prediction on physical modeling data 1.5D internal multiple prediction on physical modeling data Pan Pan, Kris Innanen and Joe Wong ABSTRACT Multiple attenuation is a key aspect of seismic

More information

Internal Multiple Attenuation on Radial Gathers With Inverse- Scattering Series Prediction

Internal Multiple Attenuation on Radial Gathers With Inverse- Scattering Series Prediction Internal Multiple Attenuation on Radial Gathers With Inverse- Scattering Series Prediction Frederico Xavier de Melo, Clément Kostov, James Wu Schlumberger Summary We present a novel workflow for pre-stack

More information

Despite the explosion of ne w, innovative technologies

Despite the explosion of ne w, innovative technologies SPECIAL Multiple SECTION: attenuation M u l tiple attenuation Elimination of land internal multiples based on the inverse scattering series YI LUO, PANOS G. KELAMIS, QIANG FU, SHOUDONG HUO, and GHADA SINDI,

More information

An illustration of adaptive Marchenko imaging

An illustration of adaptive Marchenko imaging An illustration of adaptive Marchenko imaging Joost van der Neut 1, Kees Wapenaar 1, Jan Thorbecke 1, Evert Slob 1, and Ivan Vasconcelos 2 Abstract In Marchenko imaging, wavefields are retrieved at specified

More information

Tu N Internal Multiple Attenuation on Radial Gathers With Inverse-scattering Series Prediction

Tu N Internal Multiple Attenuation on Radial Gathers With Inverse-scattering Series Prediction Tu N114 05 Internal Multiple Attenuation on Radial Gathers With Inverse-scattering Series Prediction F. Xavier de Melo* (Schlumberger), C. Kostov (Schlumberger) & J. Wu (Schlumberger) SUMMARY We present

More information

Surface-related multiple elimination with multiple prediction on GPU

Surface-related multiple elimination with multiple prediction on GPU P-68 Surface-related multiple elimination with Summary Ying Shi*, Hong Liu and Bo Li Institute of Geology and Geophysics, Chinese Academy of Sciences Surface-related multiple prediction by SRME call for

More information

Adaptive de-ghosting by kurtosis maximisation. Sergio Grion, Rob Telling, Janet Barnes, Dolphin Geophysical. Summary

Adaptive de-ghosting by kurtosis maximisation. Sergio Grion, Rob Telling, Janet Barnes, Dolphin Geophysical. Summary Adaptive de-ghosting by kurtosis maximisation Sergio Grion, Rob Telling, Janet Barnes, Dolphin Geophysical Summary We discuss the reasons for adaptive de-ghosting and its advantages. The success of broadband

More information

Interpolation using asymptote and apex shifted hyperbolic Radon transform

Interpolation using asymptote and apex shifted hyperbolic Radon transform Interpolation using asymptote and apex shifted hyperbolic Radon transform Amr Ibrahim, Paolo Terenghi and Mauricio D. Sacchi Department of Physics, University of Alberta PGS Abstract The asymptote and

More information

Practical implementation of SRME for land multiple attenuation

Practical implementation of SRME for land multiple attenuation Practical implementation of SRME for land multiple attenuation Juefu Wang* and Shaowu Wang, CGGVeritas, Calgary, Canada juefu.wang@cggveritas.com Summary We present a practical implementation of Surface

More information

We N Converted-phase Seismic Imaging - Amplitudebalancing Source-independent Imaging Conditions

We N Converted-phase Seismic Imaging - Amplitudebalancing Source-independent Imaging Conditions We N106 02 Converted-phase Seismic Imaging - Amplitudebalancing -independent Imaging Conditions A.H. Shabelansky* (Massachusetts Institute of Technology), A.E. Malcolm (Memorial University of Newfoundland)

More information

Jie Hou. Education. Research Interest True Amplitude Seismic Imaging and Inversion Acceleration of Least Squares Migration Inversion Velocity Analysis

Jie Hou. Education. Research Interest True Amplitude Seismic Imaging and Inversion Acceleration of Least Squares Migration Inversion Velocity Analysis Jie Hou Education Rice University Ph.D. Candidate in Geophysics, Earth Science 09/2012 Present China University of Petroleum(East China) 09/2008 06/2012 B.S. in Exploration Geophysics Thesis: High Order

More information

Integrated receiver deghosting and closed-loop surface-multiple elimination

Integrated receiver deghosting and closed-loop surface-multiple elimination Delft University of Technology Integrated receiver deghosting and closed-loop surface-multiple elimination Vrolijk, Jan-Willem; Verschuur, D.J.; Lopez Angarita, Gabriel DOI.9/GEO6-44. Publication date

More information

SEG Houston 2009 International Exposition and Annual Meeting

SEG Houston 2009 International Exposition and Annual Meeting Wave-equation based residual multiple prediction and elimination in migration depth domain as an aid to seismic interpretation in Wang, Manhong Guo, huck Mason, Jun ai, Sampath Gajawada, and Duryodhan

More information

Modelling of Time-varying Rough Sea Surface Ghosts and Source Deghosting by Integral Inversion

Modelling of Time-varying Rough Sea Surface Ghosts and Source Deghosting by Integral Inversion Modelling of Time-varying Rough Sea Surface Ghosts and Source Deghosting by Integral Inversion E. Cecconello* (Petroleum Geo-Services / UiO), E.G. Asgedom (Petroleum Geo-Services), O.C. Orji (Petroleum

More information

UNIVERSITY OF CALGARY. Internal multiple prediction: an application on synthetic data, physical modeling data and field. data

UNIVERSITY OF CALGARY. Internal multiple prediction: an application on synthetic data, physical modeling data and field. data UNIVERSITY OF CALGARY Internal multiple prediction: an application on synthetic data, physical modeling data and field data by Melissa Judith Hernandez Quijada A THESIS SUBMITTED TO THE FACULTY OF GRADUATE

More information

An adaptive method for source-receiver Marchenko redatuming

An adaptive method for source-receiver Marchenko redatuming An adaptive method for source-receiver Marchenko redatuming Myrna Staring (Delft University of Technology), Roberto Pereira (CGG), Huub Douma (CGG), Joost van der Neut (Delft University of Technology)

More information

Th N Deghosting by Echo-deblending SUMMARY. A.J. Berkhout* (Delft University of Technology) & G. Blacquiere (Delft University of Technology)

Th N Deghosting by Echo-deblending SUMMARY. A.J. Berkhout* (Delft University of Technology) & G. Blacquiere (Delft University of Technology) Th N103 11 Deghosting by Echo-deblending A.J. Berkhout* (Delft University of Technology) & G. Blacquiere (Delft University of Technology) SUMMARY Because of the strong sea surface reflectivity, a marine

More information

Introduction. Surface and Interbed Multtple Elimination

Introduction. Surface and Interbed Multtple Elimination Pre-stack Land Surface and Interbed Demultiple Methodology An Example from the Arabian Peninsula Roald van Borselen, Grog Fookes, Michel Schonewille, Constantine Tsingas, Michael West PGS Geophysical;

More information

Source-receiver migration of multiple reflections

Source-receiver migration of multiple reflections Stanford Exploration Project, Report 113, July 8, 2003, pages 75 85 Source-receiver migration of multiple reflections Guojian Shan 1 ABSTRACT Multiple reflections are usually considered to be noise and

More information

Stanford Exploration Project, Report 124, April 4, 2006, pages 49 66

Stanford Exploration Project, Report 124, April 4, 2006, pages 49 66 Stanford Exploration Project, Report 124, April 4, 2006, pages 49 66 48 Stanford Exploration Project, Report 124, April 4, 2006, pages 49 66 Mapping of specularly-reflected multiples to image space: An

More information

Chapter N/A: The Earthquake Cycle and Earthquake Size

Chapter N/A: The Earthquake Cycle and Earthquake Size This sheet provides a general outline of the topics covered on your upcoming exam. This is not an exhaustive list, but is simply a general list of key points. In other words, most exam questions will be

More information

Illumination compensation with Poynting vectors

Illumination compensation with Poynting vectors Illumination compensation using Poynting vectors, with special treatment for multiples Alan Richardson and Alison E. Malcolm, Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute

More information

Reverse time migration of multiples: Applications and challenges

Reverse time migration of multiples: Applications and challenges Reverse time migration of multiples: Applications and challenges Zhiping Yang 1, Jeshurun Hembd 1, Hui Chen 1, and Jing Yang 1 Abstract Marine seismic acquisitions record both primary and multiple wavefields.

More information

Modelling, migration, and inversion using Linear Algebra

Modelling, migration, and inversion using Linear Algebra Modelling, mid., and inv. using Linear Algebra Modelling, migration, and inversion using Linear Algebra John C. Bancroft and Rod Blais ABSRAC Modelling migration and inversion can all be accomplished using

More information

SUMMARY INTRODUCTION NEW METHOD

SUMMARY INTRODUCTION NEW METHOD Reverse Time Migration in the presence of known sharp interfaces Alan Richardson and Alison E. Malcolm, Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology SUMMARY

More information

Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines

Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic

More information

Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited

Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears

More information

A least-squares shot-profile application of time-lapse inverse scattering theory

A least-squares shot-profile application of time-lapse inverse scattering theory A least-squares shot-profile application of time-lapse inverse scattering theory Mostafa Naghizadeh and Kris Innanen ABSTRACT The time-lapse imaging problem is addressed using least-squares shot-profile

More information

Estimation of multiple scattering by iterative inversion, Part II: Practical aspects and examples

Estimation of multiple scattering by iterative inversion, Part II: Practical aspects and examples GEOPHYSICS, VOL. 62, NO. 5 (SEPTEMBER-OCTOBER 1997); P. 1596 1611, 18 FIGS. Estimation of multiple scattering by iterative inversion, Part II: Practical aspects and examples D. J. Verschuur* and A. J.

More information

Challenges and strategies of interbed multiple attenuation in the Asia-Pacific region

Challenges and strategies of interbed multiple attenuation in the Asia-Pacific region Challenges and strategies of interbed multiple attenuation in the Asia-Pacific region Min Wang 1*, Barry Hung 1, Xiang Li 1 and Stephane Fintz 1 demonstrate how to address the issues of interbed multiples

More information

PLEASE NOTE THAT THERE ARE QUESTIONS IN SOME OF THESE SECTIONS THAT ARE TO BE TURNED IN AS HOMEWORK.

PLEASE NOTE THAT THERE ARE QUESTIONS IN SOME OF THESE SECTIONS THAT ARE TO BE TURNED IN AS HOMEWORK. Seismic Reflection Lab for Geology 594 BEFORE YOU DO ANYTHING INSERT THE DONGLE ON THE LEFT SIDE OF THE KEYBOARD (LIGHT SIDE UP). The little green light will come on showing the dongle is activated. You

More information

Autofocusing imaging: Imaging with primaries, internal multiples and free-surface multiples

Autofocusing imaging: Imaging with primaries, internal multiples and free-surface multiples CWP-797 Autofocusing imaging: Imaging with primaries, internal multiples and free-surface multiples Satyan Singh 1, Roel Snieder 1, Jyoti Behura 1,3, Joost van Der Neut 2, Kees Wapenaar 2, and Evert Slob

More information

Attenuation of water-layer-related multiples Clement Kostov*, Richard Bisley, Ian Moore, Gary Wool, Mohamed Hegazy, Glenn Miers, Schlumberger

Attenuation of water-layer-related multiples Clement Kostov*, Richard Bisley, Ian Moore, Gary Wool, Mohamed Hegazy, Glenn Miers, Schlumberger Clement Kostov*, Richard Bisley, Ian Moore, Gary Wool, Mohamed Hegazy, Glenn Miers, Schlumberger Summary We present a method for modeling and separation of waterlayer-related multiples in towed streamer

More information

It is widely considered that, in regions with significant

It is widely considered that, in regions with significant Multifocusing-based multiple attenuation Alex Berkovitch and Kostya Deev, Geomage Evgeny Landa, OPERA It is widely considered that, in regions with significant geologic complexity, methods which work directly

More information

Towards velocity independent collapsing of diffractions

Towards velocity independent collapsing of diffractions Towards velocity independent collapsing of diffractions Early stage concepts and approximations Kris Innanen M-OSRP Annual Meeting, June 5-7, 2007 Annual Report pgs. 279-303 Key Points Liu et al., (2005-2007)

More information

Machine-learning Based Automated Fault Detection in Seismic Traces

Machine-learning Based Automated Fault Detection in Seismic Traces Machine-learning Based Automated Fault Detection in Seismic Traces Chiyuan Zhang and Charlie Frogner (MIT), Mauricio Araya-Polo and Detlef Hohl (Shell International E & P Inc.) June 9, 24 Introduction

More information

Summary. Introduction

Summary. Introduction Dmitry Alexandrov, Saint Petersburg State University; Andrey Bakulin, EXPEC Advanced Research Center, Saudi Aramco; Pierre Leger, Saudi Aramco; Boris Kashtan, Saint Petersburg State University Summary

More information

Angle Gathers for Gaussian Beam Depth Migration

Angle Gathers for Gaussian Beam Depth Migration Angle Gathers for Gaussian Beam Depth Migration Samuel Gray* Veritas DGC Inc, Calgary, Alberta, Canada Sam Gray@veritasdgc.com Abstract Summary Migrated common-image-gathers (CIG s) are of central importance

More information

Downloaded 09/20/16 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 09/20/16 to Redistribution subject to SEG license or copyright; see Terms of Use at Joint SRME and model-based water-layer demultiple for ocean bottom node Hui Huang*, Ping Wang, Jing Yang, Hui Chen (CGG); Pierre-Olivier Ariston, Imtiaz Ahmed, and Nick Bassett (BP) Summary Ocean bottom

More information

Main Menu. providing relatively fast and extremely high-quality and high-resolution performance.

Main Menu. providing relatively fast and extremely high-quality and high-resolution performance. Full-Aimuth Angle Domain Imaging Zvi Koren, Igor Ravve, Evgeny Ragoa, Allon Bartana, Paradigm Geophysical, Dan Kosloff, Tel Aviv University and Paradigm Geophysical Summary This work presents a new seismic

More information

Comments on wavefield propagation using Reverse-time and Downward continuation

Comments on wavefield propagation using Reverse-time and Downward continuation Comments on wavefield propagation using Reverse-time and Downward continuation John C. Bancroft ABSTRACT Each iteration a of Full-waveform inversion requires the migration of the difference between the

More information

Predicting rugged water-bottom multiples through wavefield extrapolation with rejection and injection

Predicting rugged water-bottom multiples through wavefield extrapolation with rejection and injection Predicting rugged water-bottom multiples through wavefield extrapolation with rejection and injection Benxi Ke ABSTRACT Although convolution-based and WEM-modeling-based methods for predicting surface-related

More information

AVA analysis based on RTM angle-domain common image gather Rui Yan* and Xiao-Bi Xie, IGPP, University of California, Santa Cruz

AVA analysis based on RTM angle-domain common image gather Rui Yan* and Xiao-Bi Xie, IGPP, University of California, Santa Cruz AVA analysis based on RTM angle-domain common image gather Rui Yan* and Xiao-Bi Xie, IGPP, University of California, Santa Cruz Summary We propose an alternative approach of AVA analysis based on RTM angle-domain

More information

Downloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at Recovering the Reflectivity Matrix and Angle-dependent Plane-wave Reflection Coefficients from Imaging of Multiples Alba Ordoñez PGS/UiO*, Walter Söllner PGS, Tilman Klüver PGS and Leiv J. Gelius UiO Summary

More information

Mitigation of the 3D Cross-line Acquisition Footprint Using Separated Wavefield Imaging of Dual-sensor Streamer Seismic

Mitigation of the 3D Cross-line Acquisition Footprint Using Separated Wavefield Imaging of Dual-sensor Streamer Seismic Mitigation of the 3D Cross-line Acquisition Footprint Using Separated Wavefield Imaging of Dual-sensor Streamer Seismic A.S. Long* (PGS), S. Lu (PGS), D. Whitmore (PGS), H. LeGleut (PGS), R. Jones (Lundin

More information

The art of well tying with new MATLAB tools

The art of well tying with new MATLAB tools Art of well tying The art of well tying with new MATLAB tools Heather J.E. Lloyd and Gary F. Margrave ABSTRACT To tie seismic data to well logs, four steps are required. The first step involves editing

More information

Introduction modeling, migration, imaging, and inversion

Introduction modeling, migration, imaging, and inversion 1 Introduction modeling, migration, imaging, and inversion 1.1 Seismic data contains information The fundamental assumption of this book is that seismic data contain decipherable information about Earth

More information

Writing Kirchhoff migration/modelling in a matrix form

Writing Kirchhoff migration/modelling in a matrix form Writing Kirchhoff migration/modelling in a matrix form Abdolnaser Yousefzadeh and John C. Bancroft Kirchhoff migration matrix ABSTRACT Kirchhoff prestack migration and modelling are linear operators. Therefore,

More information

We N Depth Domain Inversion Case Study in Complex Subsalt Area

We N Depth Domain Inversion Case Study in Complex Subsalt Area We N104 12 Depth Domain Inversion Case Study in Complex Subsalt Area L.P. Letki* (Schlumberger), J. Tang (Schlumberger) & X. Du (Schlumberger) SUMMARY Geophysical reservoir characterisation in a complex

More information

Lab 3: Depth imaging using Reverse Time Migration

Lab 3: Depth imaging using Reverse Time Migration Due Wednesday, May 1, 2013 TA: Yunyue (Elita) Li Lab 3: Depth imaging using Reverse Time Migration Your Name: Anne of Cleves ABSTRACT In this exercise you will familiarize yourself with full wave-equation

More information

M. Warner* (S-Cube), T. Nangoo (S-Cube), A. Umpleby (S-Cube), N. Shah (S-Cube), G. Yao (S-Cube)

M. Warner* (S-Cube), T. Nangoo (S-Cube), A. Umpleby (S-Cube), N. Shah (S-Cube), G. Yao (S-Cube) Tu A12 15 High-Resolution Reflection FWI M. Warner* (S-Cube), T. Nangoo (S-Cube), A. Umpleby (S-Cube), N. Shah (S-Cube), G. Yao (S-Cube) Summary We demonstrate reflection FWI on a less-than-ideal 3D narrow-azimuth

More information

Angle-domain parameters computed via weighted slant-stack

Angle-domain parameters computed via weighted slant-stack Angle-domain parameters computed via weighted slant-stack Claudio Guerra 1 INTRODUCTION Angle-domain common image gathers (ADCIGs), created from downward-continuation or reverse time migration, can provide

More information

Edinburgh Research Explorer

Edinburgh Research Explorer Edinburgh Research Explorer Multi-dimensional Free-surface Multiple Elimination and Source Deblending of Volve OBC Data Citation for published version: Ravasi, M, Vasconcelos, I, Curtis, A & Kritski, A

More information

Pressure and shear wave separation of ocean bottom seismic data

Pressure and shear wave separation of ocean bottom seismic data Pressure and shear wave separation of ocean bottom seismic data Ohad Barak ABSTRACT PZ summation is a common method for separating the upgoing wavefield from the downgoing wavefield in data acquired by

More information

Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples

Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples A. Ordoñez* (PGS), W.F. Sollner (PGS), T. Klüver (PGS) & L.G. Gelius (UiO) SUMMARY A joint

More information

Coherent partial stacking by offset continuation of 2-D prestack data

Coherent partial stacking by offset continuation of 2-D prestack data Stanford Exploration Project, Report 82, May 11, 2001, pages 1 124 Coherent partial stacking by offset continuation of 2-D prestack data Nizar Chemingui and Biondo Biondi 1 ABSTRACT Previously, we introduced

More information

Migration from a non-flat datum via reverse-time extrapolation

Migration from a non-flat datum via reverse-time extrapolation Stanford Exploration Project, Report 84, May 9, 2001, pages 1 50 Migration from a non-flat datum via reverse-time extrapolation Gopal Palacharla 1 ABSTRACT Land surveys usually have elevation changes,

More information

3-D prestack migration of common-azimuth data

3-D prestack migration of common-azimuth data Stanford Exploration Project, Report 80, May 15, 2001, pages 1 126 3-D prestack migration of common-azimuth data Biondo Biondi and Gopal Palacharla 1 ABSTRACT In principle, downward continuation of 3-D

More information

3D angle decomposition for elastic reverse time migration Yuting Duan & Paul Sava, Center for Wave Phenomena, Colorado School of Mines

3D angle decomposition for elastic reverse time migration Yuting Duan & Paul Sava, Center for Wave Phenomena, Colorado School of Mines 3D angle decomposition for elastic reverse time migration Yuting Duan & Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY We propose 3D angle decomposition methods from elastic reverse

More information

Imaging and Full Waveform Inversion of seismic data from the CO 2 gas cloud at Sleipner

Imaging and Full Waveform Inversion of seismic data from the CO 2 gas cloud at Sleipner Imaging and Full Waveform Inversion of seismic data from the CO 2 gas cloud at Sleipner Espen Birger Raknes, Børge Arntsen, and Wiktor Weibull Norwegian University of Science and Technology (NTNU) Department

More information

= R G SUMMARY. The Santos basin in Brazil suffers from strong internal multiples that overlap with primaries from the pre-salt reservoirs.

= R G SUMMARY. The Santos basin in Brazil suffers from strong internal multiples that overlap with primaries from the pre-salt reservoirs. Marchenko redatuming by adaptive double-focusing on 2D and 3D field data of the Santos basin Myrna Staring, Joost van der Neut and Kees Wapenaar (Delft University of Technology) SUMMARY The Santos basin

More information

Data reconstruction with shot-profile least-squares migration

Data reconstruction with shot-profile least-squares migration Data reconstruction with shot-profile least-squares migration Sam T. Kaplan, Mostafa Naghizadeh and Mauricio D. Sacchi Department of Physics, Room 238 CEB, 11322-89 Avenue, University of Alberta, Edmonton,

More information

P262 Limited-aperture Acquisition Compensation for Shot Profile Imaging

P262 Limited-aperture Acquisition Compensation for Shot Profile Imaging P262 Limited-aperture Acquisition Compensation for Shot Profile Imaging C.E. Yarman* (WesternGeco-Schlumberger), R. Fletcher (WesternGeco- Schlumberger) & X. Xiao (WesternGeco-Schlumberger) SUMMARY We

More information

G009 Multi-dimensional Coherency Driven Denoising of Irregular Data

G009 Multi-dimensional Coherency Driven Denoising of Irregular Data G009 Multi-dimensional Coherency Driven Denoising of Irregular Data G. Poole* (CGGVeritas Services (UK) Ltd) SUMMARY Many land and ocean bottom datasets suffer from high levels of noise which make the

More information

Time in ms. Chan. Summary

Time in ms. Chan. Summary Marine source signature estimation with dual near-field hydrophones Rob Telling*, Sergio Grion Stuart Denny & R. Gareth Williams, Shearwater GeoServices Summary We derive marine seismic signatures using

More information

Reconstruction of irregularly sampled and aliased data with linear prediction filters

Reconstruction of irregularly sampled and aliased data with linear prediction filters Reconstruction o irregularly sampled and aliased data with linear prediction ilters Mostaa Naghizadeh and Mauricio Sacchi Signal Analysis and Imaging Group (SAIG) Department o Physics University o Alberta

More information

A full-wave equation based seismic illumination analysis method

A full-wave equation based seismic illumination analysis method A full-wave equation based seismic illumination analysis method Xiao-Bi Xie and Hui Yang University of California, Santa Cruz, CA 95064, USA Summary We propose a full-wave equation based method for seismic

More information

Edinburgh Research Explorer

Edinburgh Research Explorer Edinburgh Research Explorer Marchenko Imaging of Volve Field, North Sea Citation for published version: Ravasi, M, Vasconcelos, I, Kritski, A, Curtis, A, Da Costa Filho, CA & Meles, G 5, 'Marchenko Imaging

More information

APPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA.

APPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA. APPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA. BY: ISAAC KUMA YEBOAH. Department of Engineering, Regent University

More information

SEG/New Orleans 2006 Annual Meeting

SEG/New Orleans 2006 Annual Meeting 3-D tomographic updating with automatic volume-based picking Dimitri Bevc*, Moritz Fliedner, Joel VanderKwaak, 3DGeo Development Inc. Summary Whether refining seismic images to evaluate opportunities in

More information

Multichannel deconvolution imaging condition for shot-profile migration

Multichannel deconvolution imaging condition for shot-profile migration Stanford Exploration Project, Report 113, July 8, 2003, pages 127 139 Multichannel deconvolution imaging condition for shot-profile migration Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT A significant

More information

Tu STZ1 06 Looking Beyond Surface Multiple Predictions - A Demultiple Workflow for the Culzean High Density OBC Survey

Tu STZ1 06 Looking Beyond Surface Multiple Predictions - A Demultiple Workflow for the Culzean High Density OBC Survey Tu STZ1 06 Looking Beyond Surface Multiple Predictions - A Demultiple Workflow for the Culzean High Density OBC Survey S. Gupta (Schlumberger), A. Merry (Maersk Oil), L.P. Jensen (Maersk Oil), A. Clarke

More information

Progress Report on: Interferometric Interpolation of 3D SSP Data

Progress Report on: Interferometric Interpolation of 3D SSP Data Progress Report on: Interferometric Interpolation of 3D SSP Data Sherif M. Hanafy ABSTRACT We present the theory and numerical results for interferometrically interpolating and extrapolating 3D marine

More information

Postmigration multiple prediction and removal in the depth domain

Postmigration multiple prediction and removal in the depth domain GEOPHYSICS. VOL. 76, NO. 5 (SEPTEMBER-OCTOBER 2011); P. WB217 WB223, 10 FIGS. 10.1190/GEO2011-0010.1 Postmigration multiple prediction and removal in the depth domain Bin Wang 1, Jun Cai 1, Manhong Guo

More information

G012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry

G012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry G012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry D.F. Halliday* (Schlumberger Cambridge Research), P.J. Bilsby (WesternGeco), J. Quigley (WesternGeco) & E. Kragh (Schlumberger

More information

Shallow Reverberation Prediction Methodology with SRME

Shallow Reverberation Prediction Methodology with SRME Shallow Reverberation Prediction Methodology with SRME S.R. Barnes* (Petroleum Geo-Services), R.F. Hegge (Petroleum Geo- Services), R. van Borselen (Petroleum Geo-Services) & J. Owen (Petroleum Geo-Services)

More information

Th ELI1 12 Joint Crossline Reconstruction and 3D Deghosting of Shallow Seismic Events from Multimeasurement Streamer Data

Th ELI1 12 Joint Crossline Reconstruction and 3D Deghosting of Shallow Seismic Events from Multimeasurement Streamer Data Th ELI1 12 Joint Crossline Reconstruction and 3D Deghosting of Shallow Seismic Events from Multimeasurement Streamer Data Y.I. Kamil* (Schlumberger), M. Vassallo (Schlumberger), W. Brouwer (Schlumberger),

More information

Seismic Reflection Method

Seismic Reflection Method Seismic Reflection Method 1/GPH221L9 I. Introduction and General considerations Seismic reflection is the most widely used geophysical technique. It can be used to derive important details about the geometry

More information

Stacking angle-domain common-image gathers for normalization of illumination a

Stacking angle-domain common-image gathers for normalization of illumination a Stacking angle-domain common-image gathers for normalization of illumination a a Published in Geophysical Prospecting, 59, 244-255 (2011) Guochang Liu, China University of Petroleum-Beijing and The University

More information

A Review of Current Marine Demultiple Techniques with Examples from the East Coast of Canada

A Review of Current Marine Demultiple Techniques with Examples from the East Coast of Canada A Review of Current Marine De Techniques with Examples from the East Coast of Canada ABSTRACT R.Brooymans*, T.Mojesky and L.Pham CGG Canada Services Ltd., Calgary rbrooymans@ca.cgg.com Introduction Multiple

More information

Azimuthal binning for improved fracture delineation Gabriel Perez*, Kurt J. Marfurt and Susan Nissen

Azimuthal binning for improved fracture delineation Gabriel Perez*, Kurt J. Marfurt and Susan Nissen Azimuthal binning for improved fracture delineation Gabriel Perez*, Kurt J. Marfurt and Susan issen Abstract We propose an alternate way to define azimuth binning in Kirchhoff prestack migration. This

More information

Downward propagation: using true near-surface velocity model, one-way wave- equation extrapolator:

Downward propagation: using true near-surface velocity model, one-way wave- equation extrapolator: Class 9: Beyond Seismic Statics Corrections Wed, Oct 7, 2009 Dynamic corrections: wave-equation and wavefield datuming methods Skipping the near-surface by interferometry (no model needed) Common-Focus

More information

Q-compensation in complex media ray-based and wavefield extrapolation approaches Maud Cavalca, Robin Fletcher and Marko Riedel, WesternGeco.

Q-compensation in complex media ray-based and wavefield extrapolation approaches Maud Cavalca, Robin Fletcher and Marko Riedel, WesternGeco. ray-based and wavefield extrapolation approaches Maud Cavalca, Robin Fletcher and Marko Riedel, WesternGeco. Summary We apply and compare three model-based Q-compensation approaches. The first two approaches

More information

Source deghosting for synchronized multi-level source streamer data Zhan Fu*, Nan Du, Hao Shen, Ping Wang, and Nicolas Chazalnoel (CGG)

Source deghosting for synchronized multi-level source streamer data Zhan Fu*, Nan Du, Hao Shen, Ping Wang, and Nicolas Chazalnoel (CGG) Source deghosting for synchronized multi-level source streamer data Zhan Fu*, Nan Du, Hao Shen, Ping Wang, and Nicolas Chazalnoel (CGG) Summary Ghost wavefield elimination is pivotal for improving the

More information

Downloaded 10/23/13 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 10/23/13 to Redistribution subject to SEG license or copyright; see Terms of Use at ACQUISITION APERTURE CORRECTION IN ANGLE-DOMAIN TOWARDS THE TRUE- REFLECTION RTM Rui Yan 1*, Huimin Guan 2, Xiao-Bi Xie 1, Ru-Shan Wu 1, 1 IGPP, Earth and Planetary Sciences Department, University of California,

More information

Challenges of pre-salt imaging in Brazil s Santos Basin: A case study on a variable-depth streamer data set Summary

Challenges of pre-salt imaging in Brazil s Santos Basin: A case study on a variable-depth streamer data set Summary Challenges of pre-salt imaging in Brazil s Santos Basin: A case study on a variable-depth streamer data set Jeremy Langlois, Bing Bai, and Yan Huang (CGGVeritas) Summary Recent offshore discoveries in

More information

Y015 Complementary Data-driven Methods for Interbed Demultiple of Land Data

Y015 Complementary Data-driven Methods for Interbed Demultiple of Land Data Y015 Complementary Data-driven Methods for Interbed Demultiple of Land Data S. Sonika* (WesternGeco), A. Zarkhidze (WesternGeco), J. Heim (WesternGeco) & B. Dragoset (WesternGeco) SUMMARY Interbed multiples

More information

Data preparation for prestack depth migration

Data preparation for prestack depth migration Data preparation for prestack depth migration Hugh D. Geiger Data preparation for prestack migration SUMMARY The goal of prestack depth migration is to reconstruct an image of the subsurface with the highest

More information

Seismic Time Processing. The Basis for Modern Seismic Exploration

Seismic Time Processing. The Basis for Modern Seismic Exploration The Future of E&P Seismic Time Processing The Basis for Modern Seismic Exploration Fusion is a leading provider of Seismic Processing for the oil and gas industry from field tapes through final migration.

More information

We LHR5 06 Multi-dimensional Seismic Data Decomposition for Improved Diffraction Imaging and High Resolution Interpretation

We LHR5 06 Multi-dimensional Seismic Data Decomposition for Improved Diffraction Imaging and High Resolution Interpretation We LHR5 06 Multi-dimensional Seismic Data Decomposition for Improved Diffraction Imaging and High Resolution Interpretation G. Yelin (Paradigm), B. de Ribet* (Paradigm), Y. Serfaty (Paradigm) & D. Chase

More information

Summary. Introduction

Summary. Introduction Multiple Prediction by Wavefield Extrapolation in Common-P Domain. Wang, Y. Kim, H. Guan, S. Sen, M. Guo, and K. Yoon TGS, 2500 CityWest oulevard, Suite 2000, Houston, TX 77042, US Summary Wavefield extrapolation

More information

IMAGING USING MULTI-ARRIVALS: GAUSSIAN BEAMS OR MULTI-ARRIVAL KIRCHHOFF?

IMAGING USING MULTI-ARRIVALS: GAUSSIAN BEAMS OR MULTI-ARRIVAL KIRCHHOFF? IMAGING USING MULTI-ARRIVALS: GAUSSIAN BEAMS OR MULTI-ARRIVAL KIRCHHOFF? Summary Samuel H. Gray* Veritas DGC Inc., 715 Fifth Ave. SW, Suite 2200, Calgary, AB Sam_gray@veritasdgc.com Carl Notfors Veritas

More information

Least-squares Wave-Equation Migration for Broadband Imaging

Least-squares Wave-Equation Migration for Broadband Imaging Least-squares Wave-Equation Migration for Broadband Imaging S. Lu (Petroleum Geo-Services), X. Li (Petroleum Geo-Services), A. Valenciano (Petroleum Geo-Services), N. Chemingui* (Petroleum Geo-Services),

More information

Application of 3D source deghosting and designature to deep-water ocean bottom node data Xu Li*, Jing Yang, Hui Chen, Melanie Vu, and Ping Wang (CGG)

Application of 3D source deghosting and designature to deep-water ocean bottom node data Xu Li*, Jing Yang, Hui Chen, Melanie Vu, and Ping Wang (CGG) Xu Li*, Jing Yang, Hui Chen, Melanie Vu, and Ping Wang (CGG) Summary Compared to towed-streamer data, deep-water ocean bottom node (OBN) data by nature have a broader bandwidth; however, the presence of

More information

Enhanced Angular Illumination from Separated Wavefield Imaging (SWIM)

Enhanced Angular Illumination from Separated Wavefield Imaging (SWIM) Enhanced Angular Illumination from Separated Wavefield Imaging (SWIM) S. Lu* (Petroleum Geo-Services), N.D. Whitmore (Petroleum Geo- Services), A.A. Valenciano (Petroleum Geo-Services) & N. Chemingui (Petroleum

More information

A Novel 3-D De-multiple Workflow for Shallow Water Environments - a Case Study from the Brage field, North Sea

A Novel 3-D De-multiple Workflow for Shallow Water Environments - a Case Study from the Brage field, North Sea A Novel 3-D De-multiple Workflow for Shallow Water Environments - a Case Study from the Brage field, North Sea J. Oukili* (PGS), T. Jokisch (PGS), A. Pankov (PGS), B. Farmani (PGS), G. Ronhølt (PGS), Ø.

More information