automatic selection of reference velocities for recursive depth migration hugh geiger and gary...

68
Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Upload: maude-copeland

Post on 13-Jan-2016

216 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Automatic selection of reference velocities for recursive depth migration

Hugh Geiger and Gary Margrave

CREWES Nov 2004

Page 2: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

POTSI* Sponsors:

*Pseudo-differential Operator Theory in Seismic Imaging

Page 3: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

The problem:• Many recursive wavefield extrapolators require

a limited set of reference velocities for efficient implementation

• How should these reference velocities be chosen?

example velocity profiles

Page 4: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Objectives:• Efficient computation

- a minimum number of reference velocities• Accurate wavefield extrapolation

- reference velocities ‘close’ to model velocities

too few reference velocities?

~3000 ms-1 (~230%)

Page 5: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Some specific requirements• PSPI – lower and upper bounding velocities

(vmin,vmax)- ideally minimize large interpolations (wavefield is a weighted summation)

• Split-step – more accurate focusingwith a slower velocity

too many reference velocities?

~80 ms-1 (~ 5%)~70 ms-1 (~ 4%)~70 ms-1 (~ 4%)

~320 ms-1 (~ 17%)

~40 ms-1 (~ 2%)

Page 6: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Basic approach 1: Linear progression• choose an approximate velocity spacing dV

nV=round((vmax-vmin)/dV)

vstep=(vmax-vmin)/nV• what is a good choice for dV?

- empirical testing required- reasonable for both low and high velocities?

~500 ms-1 (~20%)

~500 ms-1 (~10%)

Page 7: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

(detail forsubsequent figures)

Basic approach 2a: Geometric progression• choose an appropriate percentage step vprcnt

v(i) = (1+vprcnt)*v(i-1)- Kessigner (1992) recommends vprcnt=0.15

- start at vmin for profile?- start at vmin for complete velocity model?

(perhaps if using lookup tables)

~370 ms-1 (~15%)

~750 ms-1 (~15%)

Page 8: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Statistical method of Bagaini et al (1995)• choose a preliminary dv (geometric?)• equally spaced bins over vmin:vmax (e.g. nBtemp = 6)

• bin the velocities to give probability density Pi, Pi = 1

• optimal number of bins by statistical entropy S= PilogPi

nBopt = round(exp(S)+0.5) (e.g. nBopt = 5)

Page 9: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

• calculate cumulative probability distribution Yi, Yi = Pi

• each optimal bin to hold 1/nBopt (e.g. 0.2)- start at vmin

- linearly interpolate from temporary bin boundaries (e.g. at 0.2, 0.4, 0.6, 0.8)

• Is this optimal? - bins not necessarily close to peaks

Page 10: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

New peak search method• cluster velocities

- new cluster where jump exceeds vprcntmax

• Now, within each cluster:

• use Bagaini method for optimal number of bins nBopt

• create a new probability distribution with finer bins

• descending sort of Pi’s, choose all Pi’s where Pi < 0.9

• place vtemp at all Pi’s, include vmin,vmax

• use ‘greedy search’ to combine closely spaced Pi’s

- start search at bin spacing of 1, then 2, etc.- weighted linear average to move vtemp

• stop when at vtemp= nBopt

Page 11: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004
Page 12: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004
Page 13: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi bandlimited reflectivity

Page 14: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

PSPI with velocity clustering algorithm

data: deconpr 50 13 .0002 whiten [4 16 35 60] static –60msshot: ricker fdom 24 ghost array phsrot –68 (to zp) whiten [4 16 35 60]

Page 15: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi bandlimited reflectivity

Page 16: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi shallow reflectivity

Page 17: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Linear

Page 18: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Geometric

Page 19: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Bagaini

Page 20: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Peak Search

Page 21: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Modified Bagaini: clusters

Page 22: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi shallow reflectivity

Page 23: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Linear

Page 24: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Geometric

Page 25: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Bagaini

Page 26: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Peak Search

Page 27: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Modified Bagaini: clusters

Page 28: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Static shifts – affect focusing

Page 29: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

a) source wavefield in (x,z,t) b) reflected wavefield in (x,z,t)

x

z

t

zx

z

t

c) direct + reflected arrival at z=0 d) another perspective of (c)

xz

t

xt

(figures courtesy J. Bancroft)

horizontal reflector (blue)

Page 30: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

z=0

z=2

z=3

z=4

z=1

x

t

reflector

With a static shift of the source and/or receiver wavefield,the extrpolated wavefields will not be time coincidentat the reflector, causing Focusing and positioning errors.

Page 31: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

rec arraysource array

hard water bottom

free surface

Complications for Marmousi imaging:•free-surface and water bottom ghosting and multiples modify wavelet•source and receiver array directivity•two-way wavefield, one-way extrapolators•heterogeneous velocity

Page 32: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

220m

32m28m

0mx=400m

v=1500m/s ρ=1000kg/m3

v=1549m/s ρ=1478kg/m3

v=1598m/s ρ=1955kg/m3

v=1598m/s ρ=4000kg/m3

x=0m

Marmousi source array: 6 airguns at 8m spacing, depth 8m receiver array: 5 hydrophones at 4m spacing, depth 12m

Page 33: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

reflector

downgoing transmitted wave

upgoing reflected wave

receiver array @ 45ºreceiver array @ 0º

Modeled with finite difference code (courtesy Peter Manning) to examine response of isolated reflector at 0º and ~45º degree incidence

Page 34: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

normal incidence reflection

~45 degree incidence reflection

~60ms ~60ms

desired 24 Hzzero-phase Ricker wavelet

Marmousi airgun wavelet

After free-surface ghosting and water-bottom multiples, the Marmousi airgun wavelet propagates as ~24 Hz zero-phase Ricker with 60 ms delay.

Page 35: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

• The deconvolution chosen for the Marmousi data set is a simple spectral whitening followed by a gap deconvolution (40ms gap, 200ms operator)

• this yields a reasonable zero phase wavelet in preparation for depth imaging

Deconvolution

Page 36: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

• the receiver wavefield is then static shifted by -60ms to create an approximate zero phase wavelet

• if the receiver wavefield is extrapolated and imaged without compensating for the 60ms delay, focusing and positioning are compromised, as illustrated using a simple synthetic for a diffractor

diffractor imaging with no delay diffractor imaging with 60ms delay

Page 37: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

reflectivity x: 4000-6000 z: 0-1000

Page 38: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

PSPI whiten [4 16 35 60] cvel .2% clip 6

data: deconpr 50 13 .0002 whiten [4 16 35 60] static 0msshot: ricker fdom 24 ghost array phsrot –68 (to zp) unwhiten

Page 39: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

PSPI whiten [4 16 35 60] cvel .2% clip 6

data: deconpr 50 13 .0002 whiten [4 16 35 60] static –16msshot: ricker fdom 24 ghost array phsrot –68 (to zp) unwhiten

Page 40: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

PSPI whiten [4 16 35 60] cvel .2% clip 6

data: deconpr 50 13 .0002 whiten [4 16 35 60] static –32msshot: ricker fdom 24 ghost array phsrot –68 (to zp) unwhiten

Page 41: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

PSPI whiten [4 16 35 60] cvel .02% clip 6

data: deconpr 50 13 .0002 whiten [4 16 35 60] static –56msshot: ricker fdom 24 ghost array phsrot –68 (to zp) unwhiten

Page 42: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

PSPI whiten [4 16 35 60] cvel .02% clip 6

data: deconpr 50 13 .0002 whiten [4 16 35 60] static –56msshot: ricker fdom 24 ghost array phsrot –45 (to zp) whiten [4 16 35 60]

Page 43: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi bandlimited reflectivity- shifted to match Zhang et al.

(2003)

Page 44: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Zhang et al. (2003) – positioning not accurate

Page 45: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

PSPI reference velocities: peak search- shifted to match Zhang et al.

(2003)

Page 46: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi bandlimited reflectivity (as before)

Page 47: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi shallow reflectivity

Page 48: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi target reservoir

Page 49: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi target reservoir

Page 50: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi target reservoir

Page 51: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi target reservoir

Page 52: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi target reservoir

Page 53: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Marmousi target reservoir

Page 54: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004
Page 55: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004
Page 56: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004
Page 57: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

V=3000m/sV=2000m/s

cos taper 70°-87.5°

PSPI creates discontinuities at boundaries – smoothing may be good!

Page 58: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Conclusions• Preprocessing to zero phase, shot

modeling, and correction of static shifts important for imaging

• Optimal selection of reference velocities desired to maximize accuracy and efficiency of wavefield extrapolation

• Linear or geometric progression does not take into account distribution of velocities

• Bagaini et al. method does not necessarily pick reference velocities close to model velocities

• New peak search method selects reference velocites close to model velocities

Page 59: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Conclusions (cont)• However, Bagaini method performs well

on Marmousi!• Our PSPI implementation provides a good

standard for judging our other algorithms

Page 60: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004
Page 61: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004
Page 62: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

(-60ms static + 16.67m down)

Page 63: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

-56ms static + 8.33m down

Page 64: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

(-56ms static + 8.33m down)

Page 65: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

reflectivity x: 3000-8500 z: 0-3000

Page 66: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004
Page 67: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

Shift down 8.33m

Page 68: Automatic selection of reference velocities for recursive depth migration Hugh Geiger and Gary Margrave CREWES Nov 2004

PSPI whiten [4 16 30 55] cvel .2% clip 4

data: deconpr 50 13 .0002 whiten [4 16 35 60] static –60msshot: ricker fdom 24 ghost array phsrot –68 (to zp) whiten [4 16 35 60]