DSR Migration Velocity Analysis by Differential Semblance Optimization
|
|
- Annabelle Gardner
- 5 years ago
- Views:
Transcription
1 DSR Migration Velocity Analysis by Differential Semblance Optimization A. Khoury (Total E&P France), W. W. Symes (Rice University), P. Williamson and P. Shen (Total E&P USA Inc.) Society of Exploration Geophysicists, New Orleans, Oct 6
2 Agenda Automated construction of velocity models based on optimization and prestack imaging can succeed, given appropriate choice of objective function, and kinematically consistent prestack imaging engine Cost: O() migrations. Factors affecting accuracy include () aperture, () coherent noise in image volume originating either in data or in imaging operator.
3 Outline Objective approaches to automated velocity analysis Differential semblance for DSR migration Examples Conclusion
4 Objective approaches to automated velocity analysis 3
5 Objective Velocity Analysis Concept: identify an objective function J[v;d] so that kinematic consistency of data d, velocity v J[v;d] extremal Proposed objectives for velocity estimation from reflection data: Output Least Squares (Tarantola 984) Stack Power / Semblance (Toldi 985) Differential Semblance (Symes 986) All take form J[v; d] = W[v]d ; W[v] is a v-dependent weighting operator. 4
6 Properties of Objectives Size of typical problems must use gradient (Newton-based) optimization. objective must be smooth as function of velocity, data. Theorem: J[v;d] smooth in v,d and quadratic in d differential semblance (Stolk & Symes 3). Numerics: Output Least Squares, Semblance objectives have many spurious local minima, whereas Differential Semblance is (sometimes provably) unimodal. Upshot: Gradient-based optimization for large velocity updates from wideband reflection data Differential Semblance. 5
7 Differential semblance for DSR migration 6
8 Survey Sinking Claerbout: Reflected wavefield extrapolated to depth z = solution p(m, h, z, t) of p z = v (s, z) t + s + v (r, z) t + p r where s = m h is source x, r = m+h is receiver x, and m and h are (subsurface, not survey!) midpoint and half-offset. Initial value at z = is data: p(m, h,, t) = d(m, h,t). Assumes upcoming reflected wavefield, zero phase wavelet,... This paper: D propagation - 3D description similar with common azimuth approximation (Biondi and Palacharla, 996). 7
9 Depth Stepping Operator square roots must be approximated - we used the Generalized Screen Propagator ( GSP ) approach (Wu, 994), following DeHoop and Le Rousseau (). After discretization in depth, extrapolation given by p n+ = H n [v]p n, p = d where H n [v] = GSP approximation to DSR operator, using velocity values in nth depth slab {z n z z n+ }, and p n (m, h, t) = p(m, h,z n,t). 8
10 Imaging and Focusing Claerbout s imaging principle: Image extracted from extrapolated field at zero time, zero offset: I(m, z n ) = p n (m,, ) If velocity model is correct, then energy in zero time field (offset image gather) is focused at zero (subsurface) offset: p n (m, h, ), h [Other possible imaging principles: angle image gathers (Prucha et al 999, Sava and Rickett ), time lag gathers (Sava and Fomel, 5).] 9
11 Differential Semblance for DSR Measure focusing of offset image gathers by () scaling by h, () taking mean square: J DS [v;d] = hp n (m, h, ) m,h,n Presumption: if v kinematically consistent with data p focused at h = J DS [v;d] minimized over v. Theory: in ray-theoretic limit, presumption is correct (De Hoop & Stolk, ; De Hoop et al, 5) - need assume only absence of turning rays, consistent with DSR. [Note: analogous principle fails eg. for Kirchhoff common offset migration if multipathing occurs - imaging artifacts (Stolk & WWS 4).]
12 The Algorithm Because J DS is smooth and (perhaps) unimodal, can use rapidly convergent quasi- Newton algorithm. Limited memory variant of Broyden-Fletcher-Goldfarb-Shanno (Nocedal, 98 - available through Netlib); Velocity parametrization - bicubic splines, sigmoid representation to enforce bounds; Gradient computation - adjoint state method applied to DSR. [See abstract, references for details]
13 Examples
14 Example : Marmousi reflectivity, linear velocity Data (both examples) generated by time domain FD method for Born modeling. Source wavelet = Hz zero phase trapezoidal bandpass filter. Target velocity, used to generate data: linear, =.5 km/s at z =, = 4.5 km/s at z = 3km, represented on bicubic spline grid of 6 nodes in x ( x =.8 km) and 5 nodes in z ( z =.75 km). Initial velocity also linear, = km/s at z =, = 4.5 km/s at z = 3 km. Reflectivity = Marmousi velocity model (Versteeg and Grau, 99) minus m smoothing. Data geometry same as original. 3
15 Example : Marmousi reflectivity, linear velocity Position (km) 4 6 Position (km) 4 6 Left: Migrated image at initial v. Right: after LBFGS v updates. 4
16 Example : Marmousi reflectivity, linear velocity Position (km) 4 6 Position (km) Left: Input reflectivity ( true image). Right: image from DS velocity analysis. Good focusing and geometry in center. Some residual sag from initial velocity error remains on sides. Suggests that larger aperture more accurate v DS. 5
17 Example : Layered reflectivity, smoothed Marmousi velocity Position (km) 5 Position (km) 4 6 Modeling Inputs. Left: Smoothed Marmousi velocity (6 m smoothing width). Right: layered reflectivity. 6
18 Example : Layered reflectivity, smoothed Marmousi velocity Position (km) 5 Position (km) 4 6 Left: Linear velocity, initial guess for optimization. Right: image. 7
19 Example : Layered reflectivity, smoothed Marmousi velocity Position (km) 5 Position (km) 4 6 After LBFGS iterations. Left: velocity. Right: image. 8
20 Example : Layered reflectivity, smoothed Marmousi velocity Position (km) 5 Position (km) 4 6 Modeling Inputs. Left: Smoothed Marmousi velocity (6 m smoothing width). Right: layered reflectivity. 9
21 Example : Layered reflectivity, smoothed Marmousi velocity Offset (km) - Offset (km) - Offset (km) - Plot of image gather scaled by h and squared (sum of all such scaled, squared gathers = J DS ). Left: initial model. Center: LBFGS updates. Right: input ( true ) model. J DS for updated model actually smaller than for input model!
22 Example : Layered reflectivity, smoothed Marmousi velocity How can a wrong model generate better focus than a right model? Cause: coherent noise in image gathers. Possible causes of larger h signal in image gathers for true velocity: coherent noise in data [various remedies, but unlikely to be a factor here]; edge diffractions [remedy: taper on all axes]; mismigration of high angle events - consistent with appearance of gathers. [remedy: better propagator? Two-way RTM?].
23 Conclusions
24 Conclusions Have demonstrated a DSR-based automated VA prototype, extensible to 3D via common azimuth approximation. Cost of velocity analysis a few s of migrations. With present components, algorithm makes large model updates and greatly improves focusing of images. Available aperture affects accuracy (cf. tomography). Coherent noise in image volume, from data or imaging operator, can degrade accuracy. Action items for further research: assess effect of better (more expensive!) extrapolators in reducing operator-induced coherent noise; investigate effect and mitigation / modeling of coherent data noise; quantify aperture influence on velocity resolution. 3
25 Acknowledgements The authors are grateful to Total E&P USA for support of Alexandre Khoury s VIS research at Rice University, and for permission to present these results, to Henri Calandra for his encouragement and advice, and to Maarten de Hoop and Christiaan Stolk for many useful conversations. WWS thanks the sponsors of The Rice Inversion Project for their long-term support of related research on automation of velocity model building. 4
Automatic wave equation migration velocity analysis by differential semblance optimization
Automatic wave equation migration velocity analysis by differential semblance optimization Peng Shen, Christiaan Stolk, William W. Symes The Rice Inversion Project Rice University Thanks to Dr. Scott Morton
More informationSeismic Inversion: Progress and Prospects
Seismic Inversion: Progress and Prospects William W. Symes Rice University SEG 07 William W. Symes ( Rice University) Seismic Inversion: Progress and Prospects 24-09-2007 1 / 18 Introduction Focus: recent
More informationEquivalence of source-receiver migration and shot-profile migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 109 117 Short Note Equivalence of source-receiver migration and shot-profile migration Biondo Biondi 1 INTRODUCTION At first glance, shot
More informationTarget-oriented wavefield tomography: A field data example
Target-oriented wavefield tomography: A field data example Yaxun Tang and Biondo Biondi ABSTRACT We present a strategy for efficient migration velocity analysis in complex geological settings. The proposed
More informationAmplitude and kinematic corrections of migrated images for non-unitary imaging operators
Stanford Exploration Project, Report 113, July 8, 2003, pages 349 363 Amplitude and kinematic corrections of migrated images for non-unitary imaging operators Antoine Guitton 1 ABSTRACT Obtaining true-amplitude
More informationAutomatic Velocity Analysis via Shot Profile. Migration
Automatic Velocity Analysis via Shot Profile Migration Peng Shen 1, William W. Symes 2 1 Shell Internatioal Exploration and Production 2 Department of Computational and Applied Mathematics, Rice University
More informationKinematics of Reverse Time S-G Migration
Kinematics of Reverse Time S-G Migration William W. Symes TRIP Seminar Rice University September 2003 www.trip.caam.rice.edu 1 Agenda: explore prestack focussing properties of RT S-G migration, proper
More informationTarget-oriented wavefield tomography using demigrated Born data
Target-oriented wavefield tomography using demigrated Born data Yaxun Tang and Biondo Biondi ABSTRACT We present a method to reduce the computational cost of image-domain wavefield tomography. Instead
More informationTarget-oriented wave-equation inversion with regularization in the subsurface-offset domain
Stanford Exploration Project, Report 124, April 4, 2006, pages 1?? Target-oriented wave-equation inversion with regularization in the subsurface-offset domain Alejandro A. Valenciano ABSTRACT A complex
More informationSUMMARY. amounts to solving the projected differential equation in model space by a marching method.
Subsurface Domain Image Warping by Horizontal Contraction and its Application to Wave-Equation Migration Velocity Analysis Peng Shen, Shell International E&P, William W. Symes, Rice University SUMMARY
More informationRay based tomography using residual Stolt migration
Stanford Exploration Project, Report 11, November 11, 00, pages 1 15 Ray based tomography using residual Stolt migration Robert G. Clapp 1 ABSTRACT In complex areas, residual vertical movement is not an
More informationResidual move-out analysis with 3-D angle-domain common-image gathers
Stanford Exploration Project, Report 115, May 22, 2004, pages 191 199 Residual move-out analysis with 3-D angle-domain common-image gathers Thomas Tisserant and Biondo Biondi 1 ABSTRACT We describe a method
More informationWave-equation inversion prestack Hessian
Stanford Exploration Project, Report 125, January 16, 2007, pages 201 209 Wave-equation inversion prestack Hessian Alejandro A. Valenciano and Biondo Biondi ABSTRACT The angle-domain Hessian can be computed
More informationStanford Exploration Project, Report 111, June 9, 2002, pages INTRODUCTION THEORY
Stanford Exploration Project, Report 111, June 9, 2002, pages 231 238 Short Note Speeding up wave equation migration Robert G. Clapp 1 INTRODUCTION Wave equation migration is gaining prominence over Kirchhoff
More information3D angle gathers from wave-equation extended images Tongning Yang and Paul Sava, Center for Wave Phenomena, Colorado School of Mines
from wave-equation extended images Tongning Yang and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY We present a method to construct 3D angle gathers from extended images obtained
More informationIMAGING 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 informationLab 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 informationInversion after depth imaging
Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.
More informationAmplitude and kinematic corrections of migrated images for nonunitary imaging operators
GEOPHYSICS, VOL. 69, NO. 4 (JULY-AUGUST 2004); P. 1017 1024, 19 FIGS. 10.1190/1.1778244 Amplitude and kinematic corrections of migrated images for nonunitary imaging operators Antoine Guitton ABSTRACT
More informationSUMMARY INTRODUCTION THEORY. Objective function of ERWI
Extended Reflection Waveform Inversion via Differential Semblance Optimization Yujin Liu, William W. Symes and Zhenchun Li, China University of Petroleum (Huadong), Rice University SUMMARY Reflection-based
More informationSUMMARY. earth is characterized by strong (poro)elasticity.
with illumination compensation Tongning Yang, Center for Wave Phenomena, Colorado School of Mines, Jeffrey Shragge, School of Earth and Environment, University of Western Australia, and Paul Sava, Center
More informationWave-equation angle-domain common-image gathers for converted waves
GEOPHYSICS VOL. 73 NO. 1 JANUARY-FEBRUARY 8; P. S17 S6 17 FIGS. 1.119/1.81193 Wave-equation angle-domain common-image gathers for converted waves Daniel A. Rosales 1 Sergey Fomel Biondo L. Biondi 1 and
More informationAdaptive 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 information3D image-domain wavefield tomography using time-lag extended images
CWP-748 3D image-domain wavefield tomography using time-lag extended images Tongning Yang and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT Image-domain wavefield tomography is
More informationAVO, migration apertures Fresnel zones, stacking, Q, anisotropy, and fast inversions
AVO, migration apertures Fresnel zones, stacking, Q, anisotropy, and fast inversions John C. Bancroft and grad students University of Calgary CREWES 2003 NSERC 1 Grad students Shuang Sun Pavan Elapavuluri
More informationTarget-oriented computation of the wave-equation imaging Hessian
Stanford Exploration Project, Report 117, October 23, 2004, pages 63 77 Target-oriented computation of the wave-equation imaging Hessian Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT A target-oriented
More informationSeismic data interpolation with the offset continuation equation
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 125?? Seismic data interpolation with the offset continuation equation Sergey Fomel 1 ABSTRACT I propose a finite-difference offset
More informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Full Waveform Inversion Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 18, 2014 Bee Bednar (Panorama Technologies)
More informationSUMMARY INTRODUCTION THEORY
Illumination compensation for subsalt image-domain wavefield tomography Tongning Yang, Center for Wave Phenomena, Colorado School of Mines, Jeffrey Shragge, School of Earth and Environment, University
More informationFull waveform inversion by deconvolution gradient method
Full waveform inversion by deconvolution gradient method Fuchun Gao*, Paul Williamson, Henri Houllevigue, Total), 2012 Lei Fu Rice University November 14, 2012 Outline Introduction Method Implementation
More informationWriting 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 informationPrestack residual migration in the frequency domain
GEOPHYSICS, VOL. 68, NO. (MARCH APRIL 3); P. 634 64, 8 FIGS. 1.119/1.156733 Prestack residual migration in the frequency domain Paul C. Sava ABSTRACT Prestack Stolt residual migration can be applied to
More informationAngle-gather time migration a
Angle-gather time migration a a Published in SEP report, 1, 141-15 (1999) Sergey Fomel and Marie Prucha 1 ABSTRACT Angle-gather migration creates seismic images for different reflection angles at the reflector.
More informationLeast squares Kirchhoff depth migration with. preconditioning
Least squares Kirchhoff depth migration with preconditioning Aaron Stanton University of Alberta, Department of Physics, 4-83 CCIS, Edmonton AB T6G E (April, 3) Running head: Least Squares Migration ABSTRACT
More informationWave-equation migration velocity analysis with time-lag imaging
1 Wave-equation migration velocity analysis with time-lag imaging 2 3 Tongning Yang and Paul Sava Center for Wave Phenomena, Colorado School of Mines 4 5 6 (September 30, 2010) Running head: WEMVA with
More informationTarget-oriented wave-equation inversion
Stanford Exploration Project, Report 120, May 3, 2005, pages 23 40 Target-oriented wave-equation inversion Alejandro A. Valenciano, Biondo Biondi, and Antoine Guitton 1 ABSTRACT A target-oriented strategy
More informationRefraction Full-waveform Inversion in a Shallow Water Environment
Refraction Full-waveform Inversion in a Shallow Water Environment Z. Zou* (PGS), J. Ramos-Martínez (PGS), S. Kelly (PGS), G. Ronholt (PGS), L.T. Langlo (PGS), A. Valenciano Mavilio (PGS), N. Chemingui
More informationModel parametrization strategies for Newton-based acoustic full waveform
Model parametrization strategies for Newton-based acoustic full waveform inversion Amsalu Y. Anagaw, University of Alberta, Edmonton, Canada, aanagaw@ualberta.ca Summary This paper studies the effects
More informationOptimal Scaling of Prestack Migration
Optimal Scaling of Prestack Migration William W. Symes and Eric Dussaud Rice University, Total E&P USA SEG 07 William W. Symes and Eric Dussaud ( Rice University, Optimal Scaling Totalof E&P Prestack USA)
More informationReverse time migration in midpoint-offset coordinates
Stanford Exploration Project, Report 111, June 9, 00, pages 19 156 Short Note Reverse time migration in midpoint-offset coordinates Biondo Biondi 1 INTRODUCTION Reverse-time migration (Baysal et al., 198)
More informationFast 3D wave-equation migration-velocity analysis using the prestack exploding-reflector model
Fast 3D wave-equation migration-velocity analysis using the prestack exploding-reflector model Claudio Guerra and Biondo Biondi ABSTRACT In areas of complex geology, velocity-model definition should use
More informationIllumination-based normalization for wave-equation depth migration
Illumination-based normalization for wave-equation depth migration James E. Rickett ChevronTexaco Exploration and Production Technology Company, 6001 Bollinger Canyon Road, San Ramon, CA 94583-2324 formerly
More informationEARTH SCIENCES RESEARCH JOURNAL
EARTH SCIENCES RESEARCH JOURNAL Earth Sci. Res. J. Vol. 10, No. 2 (December 2006): 117-129 ATTENUATION OF DIFFRACTED MULTIPLES WITH AN APEX-SHIFTED TANGENT- SQUARED RADON TRANSFORM IN IMAGE SPACE Gabriel
More informationDirections in 3-D imaging - Strike, dip, both?
Stanford Exploration Project, Report 113, July 8, 2003, pages 363 369 Short Note Directions in 3-D imaging - Strike, dip, both? Marie L. Clapp 1 INTRODUCTION In an ideal world, a 3-D seismic survey would
More informationSUMMARY LEAST-SQUARES MIGRATION THEORY
Making the most out of the least (squares migration) Gaurav Dutta, Yunsong Huang, Wei Dai, Xin Wang, and G.T. Schuster King Abdullah University of Science and Technology SUMMARY Standard migration images
More informationModel-space vs. data-space normalization for recursive depth migration
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? Model-space vs. data-space normalization for recursive depth migration James Rickett 1 ABSTRACT Illumination problems caused by finite-recording
More informationMultichannel 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 informationKinematics of shot-geophone migration
Christiaan C. Stolk, Maarten V. de Hoop, William W. Symes submitted January 22nd, 2009 revised June 18th, 2009 ABSTRACT Recent analysis and synthetic examples have shown that many prestack depth migration
More informationRobustness of the scalar elastic imaging condition for converted waves
CWP-830 Robustness of the scalar elastic imaging condition for converted waves Yuting Duan & Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT For elastic reverse-time migration, one
More information=, (1) SEG/New Orleans 2006 Annual Meeting
U. Albertin* +, P. Sava ++, J. Etgen +, and M. Maharramov + + BP EPTG, Houston, Texas, ++ Colorado School of Mines, Goldin, Colorado Summary A methodology for velocity updating using one-way wavefield
More information96 Alkhalifah & Biondi
Stanford Exploration Project, Report 97, July 8, 1998, pages 95 116 The azimuth moveout operator for vertically inhomogeneous media Tariq Alkhalifah and Biondo L. Biondi 1 keywords: AMO, DMO, dip moveout
More informationTheory of 3-D angle gathers in wave-equation seismic imaging
J Petrol Explor Prod Technol (2011) 1:11 16 DOI 10.1007/s13202-011-000-8 ORIGINAL PAPER - EXPLORATION GEOPHYSICS Theory of 3-D angle gathers in wave-equation seismic imaging Sergey Fomel Received: 26 September
More informationChapter 6. Common Image Gathers
Chapter 6 Common Image Gathers All the migration methods presented in the previous chapters produce a migrated cube function of the three spatial coordinates. To perform velocity analysis and amplitude
More informationIllumination 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 informationJie 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 informationSummary. Introduction
Yaofeng He and Ru-Shan Wu Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, C 9564 Summary We propose an approach to perform migration and imaging using secondary scattered
More informationIterative resolution estimation in Kirchhoff imaging
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 369?? Iterative resolution estimation in Kirchhoff imaging Robert G. Clapp, Sergey Fomel, and Marie Prucha 1 ABSTRACT We apply iterative
More informationAngle 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 informationMitigation 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 informationF020 Methods for Computing Angle Gathers Using RTM
F020 Methods for Computing Angle Gathers Using RTM M. Vyas* (WesternGeco), X. Du (WesternGeco), E. Mobley (WesternGeco) & R. Fletcher (WesternGeco) SUMMARY Different techniques can be used to compute angle-domain
More informationEquivalent offset migration: the implementation and application update
Equivalent offset migration: the implementation and application update Xinxiang Li, Yong Xu and John C. Bancroft INTRODUCTION We have some improvements about equivalent offset migration (EOM) method. In
More informationTransformation to dip-dependent Common Image Gathers
Stanford Exploration Project, Report 11, November 11, 00, pages 65 83 Transformation to dip-dependent Common Image Gathers Biondo Biondi and William Symes 1 ABSTRACT We introduce a new transform of offset-domain
More informationSUMMARY ELASTIC SCALAR IMAGING CONDITION
Robust 3D scalar imaging condition for elastic RTM Yuting Duan, presently at Shell International Exploration and Production Inc., formerly at Center for Wave Phenomena, Colorado School of Mines Paul Sava,
More informationImage-warping waveform tomography
Geophysical Prospecting, 015, 63, 1050 1069 doi: 10.1111/1365-478.15 Image-warping waveform tomography Francesco Perrone, and Paul Sava Center for Wave Phenomena, Colorado School of Mines, Golden, CO 80401,
More informationWave-equation migration velocity analysis II: Subsalt imaging examples. Geophysical Prospecting, accepted for publication
Wave-equation migration velocity analysis II: Subsalt imaging examples Geophysical Prospecting, accepted for publication Paul Sava and Biondo Biondi Stanford Exploration Project, Mitchell Bldg., Department
More informationSUMMARY. «u scat(x, t) = 2δv 2
Reverse time migration-inversion from single-shot data Christiaan C. Stolk, University of Amsterdam, Maarten V. de Hoop, Purdue University and Tim J.P.M. Op t Root, University of Twente SUMMARY Reverse
More informationSUMMARY. u s(e, x,t)u r(e, x +,t), (1)
Image-domain and data-domain waveform tomography: a case study Esteban Díaz 1, Yuting Duan 1, Gerhard Pratt 2, and Paul Sava 1 1 Center for Wave Phenomena, Colorado School of Mines, 2 University of Western
More informationSource-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 informationWave-equation migration from topography: Imaging Husky
Stanford Exploration Project, Report 123, October 31, 2005, pages 49 56 Short Note Wave-equation migration from topography: Imaging Husky Jeff Shragge 1 INTRODUCTION Imaging land seismic data is wrought
More information# of PDE solves flops conventional 2N s N s N hx N hy N hz this paper 2N x N r N s
AVA analysis and geological dip estimation via two-way wave-equation based extended images Rajiv Kumar, Tristan van Leeuwen and Felix J. Herrmann University of British Columbia, Department of Earth, Ocean
More informationAngle-domain common-image gathers for migration velocity analysis by. wavefield-continuation imaging
Angle-domain common-image gathers for migration velocity analysis by wavefield-continuation imaging Biondo Biondi and William Symes 1 Stanford Exploration Project, Mitchell Bldg., Department of Geophysics,
More informationFull-waveform inversion for reservoir characterization: A synthetic study
CWP-889 Full-waveform inversion for reservoir characterization: A synthetic study Nishant Kamath, Ilya Tsvankin & Ehsan Zabihi Naeini ABSTRACT Most current reservoir-characterization workflows are based
More informationWave-equation migration velocity analysis. II. Subsalt imaging examples
Geophysical Prospecting, 2004, 52, 607 623 Wave-equation migration velocity analysis. II. Subsalt imaging examples P. Sava and B. Biondi Department of Geophysics, Stanford University, Mitchell Building,
More informationSEG/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 informationSUMMARY 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 informationAngle-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 informationWide-azimuth angle gathers for wave-equation migration
1 Wide-azimuth angle gathers for wave-equation migration 2 3 Paul Sava (Center for Wave Phenomena, Colorado School of Mines) Ioan Vlad (Statoil) 4 5 6 7 (September 14, 2010) GEO-2010-???? Running head:
More informationMulti-source Least-squares Migration of Gulf of Mexico Data
Multi-source Least-squares Migration of Gulf of Mexico Data Xin Wang. King Abdullah University of Science and Technology, Thuwal 955-69, Kingdom of Saudi Arabia Corresponding author is Xin Wang. E-mail
More informationSEG/San Antonio 2007 Annual Meeting
Imaging steep salt flanks by super-wide angle one-way method Xiaofeng Jia* and Ru-Shan Wu, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary The super-wide angle one-way
More informationP071 Land Data Regularization and Interpolation Using Azimuth Moveout (AMO)
P071 Land Data Regularization and Interpolation Using Azimuth Moveout (AMO) A.M. Popovici* (3DGeo Development Inc.), S. Crawley (3DGeo), D. Bevc (3DGeo) & D. Negut (Arcis Processing) SUMMARY Azimuth Moveout
More informationFast 3D wave-equation migration-velocity analysis using the. prestack exploding-reflector model
Fast 3D wave-equation migration-velocity analysis using the prestack exploding-reflector model Claudio Guerra and Biondo Biondi Formerly Stanford Exploration Project, Geophysics Department, Stanford University,
More informationChapter 1. 2-D field tests INTRODUCTION AND SUMMARY
Chapter 1 2-D field tests INTRODUCTION AND SUMMARY The tomography method described in the preceding chapter is suited for a particular class of problem. Generating raypaths and picking reflectors requires
More informationAttenuation of diffracted multiples with an apex-shifted tangent-squared radon transform in image space
Attenuation of diffracted multiples with an apex-shifted tangent-squared radon transform in image space Gabriel Alvarez, Biondo Biondi, and Antoine Guitton 1 ABSTRACT We propose to attenuate diffracted
More informationSUMMARY METHOD. d(t 2 = τ 2 + x2
Yujin Liu, Zhi Peng, William W. Symes and Wotao Yin, China University of Petroleum (Huadong), Rice University SUMMARY The Radon transform suffers from the typical problems of loss of resolution and aliasing
More informationLeast squares Kirchhoff depth migration: important details
Least squares Kirchhoff depth migration: important details Daniel Trad CREWES-University of Calgary Summary Least squares migration has been an important research topic in the academia for about two decades,
More informationSummary. Figure 1: Simplified CRS-based imaging workflow. This paper deals with the boxes highlighted in green.
Smoothing and automated picking of kinematic wavefield attributes Tilman Klüver and Jürgen Mann, Geophysical Institute, University of Karlsruhe, Germany Copyright 2005, SBGf Sociedade Brasiliera de Geofísica
More information3-D vertical cable processing using EOM
Carlos Rodriguez-Suarez, John C. Bancroft, Yong Xu and Robert R. Stewart ABSTRACT Three-dimensional seismic data using vertical cables was modeled and processed using equivalent offset migration (EOM),
More informationLinearized wave-equation migration velocity analysis by image warping
CWP-708 Linearized wave-equation migration velocity analysis by image warping Francesco Perrone 1, Paul Sava 1, Clara Andreoletti 2, and Nicola Bienati 2 1 Center for Wave Phenomena, Colorado School of
More informationMultiple attenuation in the image space
GEOPHYSICS, VOL. 70, NO. 1 (JANUARY-FEBRUARY 2005); P. V10 V20, 15 FIGS. 10.1190/1.1852789 Multiple attenuation in the image space Paul Sava 1 and Antoine Guitton 1 ABSTRACT Multiples can be suppressed
More informationSUMMARY THEORY INTRODUCTION
Acoustic 3D least-squares reverse time migration using the energy norm Daniel Rocha, Paul Sava & Antoine Guitton Center for Wave Phenomena, Colorado School of Mines SUMMARY We propose a least-squares reverse
More informationMigration 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 informationEnhanced 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 informationWhich Depth Imaging Method Should You Use? A Roadmap In The Maze of 3 D Depth Imaging
Which Depth Imaging Method Should You Use? A Roadmap In The Maze of 3 D Depth Imaging By, Dimitri Bevc, 3DGeo Development Inc., 4633 Old Ironsides Drive, Suite 401, Santa Clara, CA 95054, dimitri@3dgeo.com
More informationLow-rank representation of extended image volumes applications to imaging and velocity continuation
Low-rank representation of extended image volumes applications to imaging and velocity continuation Rajiv Kumar 1, Marie Graff-Kray 2, Tristan van Leeuwen 3, and Felix J. Herrmann 1 1 Georgia Institute
More information3D 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 informationScalar imaging condition for elastic reverse time migration
GEOPHYSICS, VOL. 80, NO. 4 (JULY-AUGUST 2015); P. S127 S136, 17 FIGS. 10.1190/GEO2014-0453.1 Scalar imaging condition for elastic reverse time migration Yuting Duan 1 and Paul Sava 1 ABSTRACT Polarity
More informationReconciling processing and inversion: Multiple attenuation prior to wave-equation inversion
Reconciling processing and inversion: Multiple attenuation prior to wave-equation inversion Claudio Guerra and Alejandro Valenciano ABSTRACT Seismic inversion is very sensitive to the presence of noise.
More informationLeast-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 informationPre-Stack Exploding-Reflector Model
Chapter 2 Pre-Stack Exploding-Reflector Model This chapter introduces the pre-stack exploding-reflector model (PERM). PERM uses exploding reflectors as the initial conditions to synthesize data, in a manner
More informationStacking 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