350 likes | 766 Views
Prestack Migration. Intuitive. Least Squares. Green’s Theorem. Migration. 3D Prestack Diffraction Stack Migration. Motivation: ZO only good if no lateral vel change. s. g. x. . =. d(x’ , + ). x g. s x. s, g. 3D Prestack Diffraction Stack Migration.
E N D
Prestack Migration Intuitive Least Squares Green’s Theorem Migration
3D Prestack Diffraction Stack Migration Motivation: ZO only good if no lateral vel change s g x
= d(x’, + ) xg sx s,g 3D Prestack Diffraction Stack Migration m(x) = Trial image pt x s g x
Prestack DS Migration Theory Outline RTM RTM vsPoststackvsPrestack MATLAB Code
Question: Why Prestack when poststack migration seems good enough? Prestack Migration Answer: Stacking to get stacked section assumes layered medium assumption. Solution: Migrate shot gathers so no layer assumption needed. This is prestack migration.
+ sx xg + T(s,g) = s,g sx xg g s sx xg d(s,g, ) x Diffraction Stack Migration: Prestack Up time Down time Where is scatterer? Narrow band case: direct wave correlated with data
m = L d T x ~ ~ W( ) m(x) 2 w 115. Diffraction Stack Modeling: Prestack d= L m i i ~ sx e xg e d(s,g) = A(s,x) A(g,x)
m = L d T d d s,g s,g ~ ~ W( ) W( )=1 ò ò Broadband case d(s,g, + ) = sx xg A(x,g) A(s,x) 2 w 115. Diffraction Stack Migration: Prestack - - i i * sx ~ e xg e ~ m(x) = d(s,g) A(s,x) A(g,x) .. m(x) Narrow band case: direct wave correlated with data
Prestack DS Migration Theory Outline RTM RTM vsPoststackvsPrestack MATLAB Code
m(x) = d (g, ) xg g traces ZO Diffraction Stack Migration Trial image pt x g x
m(x) = d (g, ) xg g traces ZO Diffraction Stack Migration Trial image pt x 2D dot product of migration Operator and d(g,t) g Migration Image x
m(x) = d (g, ) xg g traces ZO Reverse Time Migration Trial image pt x ~Scattered RTM g Super-resolution x x
m(x) = d (g, ) xg g traces ZO Reverse Time Migration Trial image pt x ~Scattered RTM g Super-resolution x x
Prestack DS Migration Theory Outline RTM RTM vsPoststackvsPrestack MATLAB Code
Shortest path ray Maximum energy ray Shortest traveltime ray Shortest Traveltime or Shortest Raypath Maximum Energy Traveltimes Types of Traveltimes
Prestack DS Migration Theory Outline RTM RTM vsPoststackvsPrestack DS MATLAB Code
MATLAB Inefficient Prestack Migration for isx=1:nx % Loop over shot for igx=1:nx % Loop over receivers for ix=1:nx % Loop over model x for iz=1:nx % Loop over model z t=timer(ix,iz,isx)+timer(ix,iz,igx) sample=gather(isx,igx,t) % Shot gather has 2 time derivatives mig(ix,iz)=mig(ix,iz)+sample end end end end
1. No assumption about velocity model Prestack Migration 2. More sensitive to velocity model errors compared to poststack migration 2 6 3. More than 10 – 10 times slower than poststack migration 4. More sensitive to velocity model than time migration