Top Banner
Green’s function retrieval by iterative substitution or inversion (?) of the Marchenko equation Joost van der Neut (Delft University of Technology)
41

Green’s function retrieval by iterative substitution or inversion (

Jan 18, 2018

Download

Documents

Acknowledgements Satyan Singh, Jyoti Behura, Roel Snieder Filippo Broggini, Dirk-Jan van Manen Kees Wapenaar, Evert Slob Jan Thorbecke Ivan Vasconcelos Bowen Guo, Jerry Schuster Andrey Bakulin
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Green’s function retrieval by iterative substitution or inversion (

Green’s function retrieval by iterative substitution or inversion (?) of the Marchenko equation

Joost van der Neut(Delft University of Technology)

Page 2: Green’s function retrieval by iterative substitution or inversion (

Acknowledgements

Kees Wapenaar, Evert SlobJan Thorbecke

Satyan Singh, Jyoti Behura, Roel Snieder

Ivan Vasconcelos

Filippo Broggini, Dirk-Jan van Manen

Bowen Guo, Jerry Schuster

Andrey Bakulin

Page 3: Green’s function retrieval by iterative substitution or inversion (

Introduction

Page 4: Green’s function retrieval by iterative substitution or inversion (

Input for Marchenko Redatuming

Data Multiple-freeWavelet-freeGhost-freeData

SourceFunction

SurfaceReflectivity

Data

1^^^^ ^

(Van Groenestijn & Verschuur, 2009;Lin and Herrmann, 2013)

EPSI output

Page 5: Green’s function retrieval by iterative substitution or inversion (

Input for Marchenko Redatuming

Data Multiple-freeWavelet-freeGhost-freeData

SourceFunction

SurfaceReflectivity

Data

1

2

Focal point

^^^^ ^

(Van Groenestijn & Verschuur, 2009;Lin and Herrmann, 2013)

EPSI output

BackgroundGreen’s function

Page 6: Green’s function retrieval by iterative substitution or inversion (

Marchenkoredatuming

X

gup

gdown

g0

The aim of Marchenko redatuming

Page 7: Green’s function retrieval by iterative substitution or inversion (

gup gdown

= *

Retrieve e.g. by inversion

Xdatum

Redatuming below a complex overburden

Page 8: Green’s function retrieval by iterative substitution or inversion (

Examples

Page 9: Green’s function retrieval by iterative substitution or inversion (

R T0

Conventional image Marchenko image

Example 1

Page 10: Green’s function retrieval by iterative substitution or inversion (

Example 2

Page 11: Green’s function retrieval by iterative substitution or inversion (

Example 2 – Conventional image

Page 12: Green’s function retrieval by iterative substitution or inversion (

Example 2 – Marchenko image

Page 13: Green’s function retrieval by iterative substitution or inversion (

Model

Target area Image

Example 3 - Conventional

Page 14: Green’s function retrieval by iterative substitution or inversion (

Model

Target area Image

Example 3 – Marchenko(with adaptive subtraction)

Page 15: Green’s function retrieval by iterative substitution or inversion (

Example 4 - ConventionalVelocity Density Image

Red=

without multiples

Yellow=

with multiples

Yellow Red

7000m/s

0 4000kg / m3

0

Page 16: Green’s function retrieval by iterative substitution or inversion (

Example 4 - MarchenkoVelocity Density Image

Red=

without multiples

Yellow=

Marchenko result

Yellow Red

7000m/s

0 4000kg / m3

0

Page 17: Green’s function retrieval by iterative substitution or inversion (

Example 4 with erroneous (constant) velocity - ConventionalVelocity Density Image

Red=

without multiples

Yellow=

with multiples

Yellow Red

7000m/s

0 4000kg / m3

0

Page 18: Green’s function retrieval by iterative substitution or inversion (

Example 4 with erroneous (constant) velocity – MarchenkoVelocity Density Image

Red=

without multiples

Yellow=

Marchenko result

Yellow Red

7000m/s

0 4000kg / m3

0

Page 19: Green’s function retrieval by iterative substitution or inversion (

Example 5 - Sigsbee

Behura et al. (2014)

ConventionalMarchenko

Page 20: Green’s function retrieval by iterative substitution or inversion (

Theory

Page 21: Green’s function retrieval by iterative substitution or inversion (

The focusing function (Wapenaar et al., 2014)

Focal point

Focusing functionResponse Earth

Page 22: Green’s function retrieval by iterative substitution or inversion (

Heterogeneous

The focusing function (Wapenaar et al., 2014)

Focal point

Focusing functionResponse Earth

Page 23: Green’s function retrieval by iterative substitution or inversion (

The focusing function (Wapenaar et al., 2014)

Response

Heterogeneous

Focal point

Focusing functionEarth

Page 24: Green’s function retrieval by iterative substitution or inversion (

Data acting on the focusing function

Reflection response

Heterogeneous

Focal point

Focusing function

Heterogeneous

Time-reversedFocusing function& Green’s function

Page 25: Green’s function retrieval by iterative substitution or inversion (

Representation in matrix-vector notation

Time-reversal Focusingfunction

Green’sfunction

Convolutionwith data

Focusingfunction

Page 26: Green’s function retrieval by iterative substitution or inversion (

Representation in matrix-vector notation

Time-reversal Focusingfunction

Green’sfunction

Convolutionwith data

Focusingfunction

Page 27: Green’s function retrieval by iterative substitution or inversion (

Representation in matrix-vector notation

Time-reversal Focusingfunction

Green’sfunction

Convolutionwith data

Focusingfunction

Unknown Unknown

Page 28: Green’s function retrieval by iterative substitution or inversion (

Exploiting causality

Time-reversedFocusing function

Green’s function

Page 29: Green’s function retrieval by iterative substitution or inversion (

Exploiting causalityWindowfunction

10

Time-reversedFocusing function

Green’s function

Page 30: Green’s function retrieval by iterative substitution or inversion (

Exploiting causality

Windowfunction

Focusingfunction

10

Timereversal

Page 31: Green’s function retrieval by iterative substitution or inversion (

Exploiting causality

10

Windowfunction

Green’sfunction

Page 32: Green’s function retrieval by iterative substitution or inversion (

Exploiting causality

10

01

Windowfunction

BackgroundGreen’s function

Green’sfunction

Page 33: Green’s function retrieval by iterative substitution or inversion (

Green’s function representation

Apply

Intial Green’sfunction (model)

Window Convolutionwith data

Focusingfunction

TimeReversal

Page 34: Green’s function retrieval by iterative substitution or inversion (

Green’s function representation

Apply

Intial Green’sfunction (model)

Window Convolutionwith data

Focusingfunction

TimeReversal

Page 35: Green’s function retrieval by iterative substitution or inversion (

Velocity Density

7000m/s

0 4000kg / m3

0

Retrieved Green’s functions - 1D results from spgl1

Iterative solution(50 iterations)

spgl1 inversion(100 iterations)

Page 36: Green’s function retrieval by iterative substitution or inversion (

Applications for Least-Squares Migration?

Page 37: Green’s function retrieval by iterative substitution or inversion (

Least-Squares Migration

perturbation

Data =

Page 38: Green’s function retrieval by iterative substitution or inversion (

Born approximation

perturbation

Data =

Page 39: Green’s function retrieval by iterative substitution or inversion (

Using Marchenko-based Green’s functions?

perturbation

Data =

Marchenko mapping:

Page 40: Green’s function retrieval by iterative substitution or inversion (

Discussion

Page 41: Green’s function retrieval by iterative substitution or inversion (

Revisiting the problem

1. EPSI (inversion):

2. Focusing functionRetrieval (inversion):

3. Green’s functionRetrieval (forward):

4. Least-SquaresMigration(inversion): δδ