WS 09 Velocity Model Building

Size: px
Start display at page:

Download "WS 09 Velocity Model Building"

Transcription

1

2 Full wave vs Kinematic approaches Goal in mid 80 s: Full Waveform Inversion (FWI) ~ Tomography plus Migration. There seems to exist considerable confusion regarding the sort of information about the subsurface which one might expect to extract using the L2 approach to the inverse problem of reflection seismology, and also regarding the quality of the L2 results relative to the output of conventional processing methods. The relation between L2 inversion and migration of both stacked and unstacked data is now well understood, at least in principle [see Lailly (1984) or Tarantola (1984) and Beylkin (1985)]. On the other hand, the possibility of reliable subsurface parameter estimation from realistically band-limited data seems to arise various opinions. One reads several compelling arguments in the literature that extraction of velocity trends (the main out-of-passband components of interest in reflection seismology from band-limited (low-cut) data by L2 inversion is impossible, or so difficult as to be infeasible (Tarantola, 1986). One also encounters convincing arguments in various successful stories of FWI by different groups around the turn of the new century on synthetics and on real data when both reflections and diving waves (refractions) are available. F. Santosa, W.W. Symes, 1989, An analysis of least-squares velocity inversion, 1989, SEG geophyiscal Monograph, number 4.

3 Full wave vs Kinematic approaches However, capturing long-to-intermediate wavelength information from reflection data with a least-square waveform misfit is challenging (may be impossible with a local optimization techniques) Limitations of the least-square waveform inversion formulation have been identified - Cycle skipping issue requiring a decent initial model (kinematically compatible) - Capturing long-wavelength content from reflected waves is ambiguous and difficult - Lack of intermediate wavelength information when only reflected waves are recorded Reformulations of FWI for overcoming these two pitiful features - relaxing cycle skipping - capture reflection travel-time information in a waveform inversion manner (leading to different automatic migration velocity analysis). While promising, these approaches are challenging in practice. We have chosen here to focus on the second aspect, which only start now to appear as an industrial solution.

4 Program Part 1: Overview 9h05-9h40: G. Chavent: Data Space Reflectivity and the Migration based Travel Time approach to FWI 9h40-10h15: B. Symes: Extended waveform inversion 10h15-10h30 : discussion Part 2: Methods/Examples 10h30-10h50 T. Alkhalifah: Optimizing the coefficients of the leading terms of the Born Series: FWI+MVA+more 10h50-11h10 V. Cheverda: FWI for elastic media: macrovelocity reconstruction 11h10-11h30 A. Gomes: Extending the reach of FWI with reflection data: Potential and challenges 11h30-11h50 B. Biondi: Application of tomographic FWI (TFWI) to large-scale field datasets: challenges and insights. 11h50-12h10: discussion 12h10-13h15: Lunch

5 Program Part 3: Applications 13h20-13h40 U. Albertin: Advances in and Challenges for Broadband Multiparameter Full Waveform Inversion 13h40-14h00: M.Warner: Reflection FWI without a good starting model 14h00:14h20: D. Vigh: Updating velocity fields beyond the diving waves 14h20-14h40: P. Nandi: Wave-Equation Migration Velocity Analysis in the Surface Common Offset Domain: Application to Viking Graben Field Data 14h40-15h00: Discussion Part 4: Algorithms/Methods 15h00-15h20: R. Soubaras: Mitigating the gradient artefacts of Migration Velocity Analysis by Gauss-Newton update 15h20-15h40 R.Valensi: Reflection Waveform Inversion method: solutions to the reflectivity-background coupling problem and consequences on the convergence 15h40-16h00 R. Brossier: From RWI to JFWI: including diving waves in reflectionbased velocity model building 16h00:16h20: H. Chauris: Image-domain versus data-domain velocity analysis based on true-amplitude subsurface extended migration 16h20-17h00: discussion

6 Questions we may discuss Misfit function (reformulation): How should we extend the model? In model space or in data space, with or without the Born approximation? Does it really matter? Is the migration step necessary for any reflection based full wave approach? What is feasible in 3D? How to capture both reflection and refraction information? What are the crosstalks between modeled and observed events? How should we compare the data (time pixel, trace, gather, )? What can we borrow from the progresses in conventional FWI?

7 Questions we may discuss What are the main potential benefits? What is the role of phases and amplitudes? Joint or sequential inversion of phases and amplitude? How does it compare with a picking-based approach? What will we need to have solutions as efficient as ray based tomography? Should we combine picking-based approach and waveform approach? Is coherency enhancement and even phase identification crucial (at which stage of the inversion)?

8 Questions we may discuss Which physics should we consider? Is an acoustic approach sufficient? How do we handle AVO effects? What kind of model do we retrieve (background, impedances, ) Should we consider multi-parameter inversion? What data pre-processing? Do we need a data hierarchization (primary only, separation betwen reflection and transmission/refraction, ) What is the role of the low frequencies? Other items?

Full waveform inversion guided migration velocity analysis Thibaut Allemand* and Gilles Lambaré (CGG)

Full waveform inversion guided migration velocity analysis Thibaut Allemand* and Gilles Lambaré (CGG) Thibaut Allemand* and Gilles Lambaré (CGG) Summary While difficulties of full waveform inversion (FWI) using reflected waves only have been well identified, the combination of FWI with migration velocity

More information

Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method

Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Wenyong Pan 1, Kris Innanen 1 and Wenyuan Liao 2 1. CREWES Project, Department of Geoscience,

More information

Improvements in time domain FWI and its applications Kwangjin Yoon*, Sang Suh, James Cai and Bin Wang, TGS

Improvements in time domain FWI and its applications Kwangjin Yoon*, Sang Suh, James Cai and Bin Wang, TGS Downloaded 0/7/13 to 05.196.179.38. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ Improvements in time domain FWI and its applications Kwangjin Yoon*,

More information

We B3 12 Full Waveform Inversion for Reservoir Characterization - A Synthetic Study

We B3 12 Full Waveform Inversion for Reservoir Characterization - A Synthetic Study We B3 12 Full Waveform Inversion for Reservoir Characterization - A Synthetic Study E. Zabihi Naeini* (Ikon Science), N. Kamath (Colorado School of Mines), I. Tsvankin (Colorado School of Mines), T. Alkhalifah

More information

Refraction Full-waveform Inversion in a Shallow Water Environment

Refraction 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 information

Building starting model for full waveform inversion from wide-aperture data by stereotomography

Building starting model for full waveform inversion from wide-aperture data by stereotomography Building starting model for full waveform inversion from wide-aperture data by stereotomography Vincent Prieux 1, G. Lambaré 2, S. Operto 1 and Jean Virieux 3 1 Géosciences Azur - CNRS - UNSA, France;

More information

Seismic Inversion: Progress and Prospects

Seismic 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 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

Model parametrization strategies for Newton-based acoustic full waveform

Model 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 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

SUMMARY INTRODUCTION THEORY. Objective function of ERWI

SUMMARY 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 information

Target-oriented wavefield tomography: A field data example

Target-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 information

Full-waveform inversion for reservoir characterization: A synthetic study

Full-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 information

WS10-A03 Alternative Multicomponent Observables for Robust Full-waveform Inversion - Application to a Near-surface Example

WS10-A03 Alternative Multicomponent Observables for Robust Full-waveform Inversion - Application to a Near-surface Example 1-4 June IFEMA Madrid WS-A3 Alternative Multicomponent Observables for Robust Full-waveform Inversion - Application to a Near-surface Example R. Valensi* (formerly IFSTTAR), R. Brossier (Univ. Grenoble

More information

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

Downloaded 10/26/15 to Redistribution subject to SEG license or copyright; see Terms of Use at Reflection Full-waveform Inversion for Inaccurate Starting Models Abdullah AlTheyab* and G. T. Schuster, King Abdullah University of Science and Technology (KAUST) SUMMARY We propose a seismic-reflection

More information

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

Downloaded 09/16/13 to Redistribution subject to SEG license or copyright; see Terms of Use at Time-domain incomplete Gauss-Newton full-waveform inversion of Gulf of Mexico data Abdullah AlTheyab*, Xin Wang, Gerard T. Schuster, King Abdullah University of Science and Technology Downloaded 9// to

More information

F031 Application of Time Domain and Single Frequency Waveform Inversion to Real Data

F031 Application of Time Domain and Single Frequency Waveform Inversion to Real Data F031 Application of Time Domain and Single Frequency Waveform Inversion to Real Data D Yingst (ION Geophysical), C. Wang* (ION Geophysical), J. Park (ION Geophysical), R. Bloor (ION Geophysical), J. Leveille

More information

Successful application of joint reflection/refraction tomographic velocity inversion in a shallow water marine environment.

Successful application of joint reflection/refraction tomographic velocity inversion in a shallow water marine environment. Successful application of joint reflection/refraction tomographic velocity inversion in a shallow water marine environment. Sergey Birdus 1, Dean Criddle 2, Alexey Artyomov 1, Li Li 1, Qingbing Tang 1,

More information

Target-oriented wavefield tomography using demigrated Born data

Target-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 information

Stochastic conjugate gradient method for least-square seismic inversion problems Wei Huang*, Hua-Wei Zhou, University of Houston

Stochastic conjugate gradient method for least-square seismic inversion problems Wei Huang*, Hua-Wei Zhou, University of Houston Stochastic conjugate gradient method for least-square seismic inversion problems Wei Huang*, Hua-Wei Zhou, University of Houston Summary With the development of computational power, there has been an increased

More information

Source Estimation for Wavefield Reconstruction Inversion

Source Estimation for Wavefield Reconstruction Inversion Source Estimation for Wavefield Reconstruction Inversion Zhilong Fang * and Felix J. Herrmann * * Seismic Laboratory for Imaging and Modeling SLIM), University of British Columbia Abstract Wavefield reconstruction

More information

Seismic Modeling, Migration and Velocity Inversion

Seismic 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 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

A new take on FWI: Wavefield Reconstruction Inversion

A new take on FWI: Wavefield Reconstruction Inversion A new take on FWI: Wavefield Reconstruction Inversion T. van Leeuwen 1, F.J. Herrmann and B. Peters 1 Centrum Wiskunde & Informatica, Amsterdam, The Netherlands University of British Columbia, dept. of

More information

Iterative resolution estimation in Kirchhoff imaging

Iterative 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 information

A projected Hessian matrix for full waveform inversion Yong Ma and Dave Hale, Center for Wave Phenomena, Colorado School of Mines

A projected Hessian matrix for full waveform inversion Yong Ma and Dave Hale, Center for Wave Phenomena, Colorado School of Mines A projected Hessian matrix for full waveform inversion Yong Ma and Dave Hale, Center for Wave Phenomena, Colorado School of Mines SUMMARY A Hessian matrix in full waveform inversion (FWI) is difficult

More information

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

Downloaded 09/03/13 to Redistribution subject to SEG license or copyright; see Terms of Use at Full-waveform inversion in a shallow water environment: A North Sea 3D towed-streamer data example Kathy Zou*, Lars Tore Langlo, Grunde Rønholt, Jaime Ramos-Martinez and Steve Kelly, PGS Summary We apply

More information

SUMMARY. u s(e, x,t)u r(e, x +,t), (1)

SUMMARY. 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 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

SUMMARY. earth is characterized by strong (poro)elasticity.

SUMMARY. 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 information

SUMMARY INTRODUCTION THEORY

SUMMARY 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 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

Th G Surface-offset RTM Gathers - What Happens When the Velocity is Wrong?

Th G Surface-offset RTM Gathers - What Happens When the Velocity is Wrong? Th G103 01 Surface-offset RTM Gathers - What Happens When the Velocity is Wrong? J.P. Montel* (CGG) SUMMARY Surface offset migrated common image gathers (SOCIGs) built by wave equation migration (WEM)

More information

G004 Robust Frequency-domain Multi-parameter Elastic Full Waveform Inversion

G004 Robust Frequency-domain Multi-parameter Elastic Full Waveform Inversion G004 Robust Frequency-domain Multi-parameter Elastic Full Waveform Inversion R. Brossier* (Joseph Fourier University), S. Operto (CNRS/University of Nice-Sophia-Antipolis/Géoazur) & J. Virieux (Joseph

More information

H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography

H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography M. Cavalca* (Schlumberger - Westerngeco) & R.P. Fletcher (Schlumberger - Westerngeco) SUMMARY Estimation of the

More information

Full waveform inversion of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC)

Full waveform inversion of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC) of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC) Summary (FWI) is a promising technology for next generation model building. However, it faces with many

More information

Image-warping waveform tomography

Image-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 information

Building starting models for full waveform inversion from wide-aperture data by stereotomography

Building starting models for full waveform inversion from wide-aperture data by stereotomography Geophysical Prospecting doi: 10.1111/j.1365-2478.2012.01099.x Building starting models for full waveform inversion from wide-aperture data by stereotomography Vincent Prieux 1, Gilles Lambaré 2,Stéphane

More information

Least squares Kirchhoff depth migration with. preconditioning

Least 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 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

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

SUMMARY. amounts to solving the projected differential equation in model space by a marching method.

SUMMARY. 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 information

Z-99 3D Sub-salt Tomography Based on Wave Equation Migration Perturbation Scans

Z-99 3D Sub-salt Tomography Based on Wave Equation Migration Perturbation Scans 1 Z-99 3D Sub-salt Tomography Based on Wave Equation Migration Perturbation Scans BIN WANG 1, VOLKER DIRKS 1, PATRICE GUILLAUME 2, FRANÇOIS AUDEBERT 1, ANNING HOU 1 AND DURYODHAN EPILI 1 1 CGG Americas;

More information

Target-oriented wave-equation inversion

Target-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 information

Wave-equation inversion prestack Hessian

Wave-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 information

Geophysics 1 TOWARDS AN AUTOMATIC FULL WAVE INVERSION: SYNTHETIC STUDY CASES

Geophysics 1 TOWARDS AN AUTOMATIC FULL WAVE INVERSION: SYNTHETIC STUDY CASES Geophysics 1 TOWARDS AN AUTOMATIC FULL WAVE INVERSION: SYNTHETIC STUDY CASES Jean Kormann [1,2]: jean.kormann@unileoben.ac.at, Juan Esteban Rodríguez [1]: juan.rodriguez@bsc.es, Natalia Gutierrez [1]:

More information

GEOPHYS 242: Near Surface Geophysical Imaging. Class 4: First-Arrival Traveltime Tomography Mon, April 11, 2011

GEOPHYS 242: Near Surface Geophysical Imaging. Class 4: First-Arrival Traveltime Tomography Mon, April 11, 2011 GEOPHYS 242: Near Surface Geophysical Imaging Class 4: First-Arrival Traveltime Tomography Mon, April 11, 2011 Wavefront tracing methods - speed versus accuracy Inversion algorithms - any magic approach?

More information

Multicomponent land data pre-processing for FWI: a benchmark dataset

Multicomponent land data pre-processing for FWI: a benchmark dataset Multicomponent land data pre-processing for FWI: a benchmark dataset Raul Cova, Bernie K. Law and Kris Innanen CRWES/University of Calgary Summary Successful full-waveform inversion (FWI) studies using

More information

Building accurate initial models using gain functions for. waveform inversion in the Laplace domain

Building accurate initial models using gain functions for. waveform inversion in the Laplace domain Building accurate initial models using gain functions for waveform inversion in the Laplace domain Wansoo Ha 1 and Changsoo Shin 2 1 Department of Energy Resources Engineering, Pukyong National University

More information

u = v is the Laplacian and represents the sum of the second order derivatives of the wavefield spatially.

u = v is the Laplacian and represents the sum of the second order derivatives of the wavefield spatially. Seismic Elastic Modelling Satinder Chopra, Arcis Corporation, Calgary, Canada Summary Seismic modeling experiments were performed to understand the seismic responses associated with a variety of subsurface

More information

Incorporating FWI velocities in Simulated Annealing based acoustic impedance inversion

Incorporating FWI velocities in Simulated Annealing based acoustic impedance inversion Incorporating FWI velocities in Simulated Annealing based acoustic impedance inversion Date: 27 September 2018 Authors: Nasser Bani Hassan, Sean McQuaid Deterministic vs Stochastic Inversion Deterministic

More information

Post-stack iterative modeling migration and inversion (IMMI)

Post-stack iterative modeling migration and inversion (IMMI) Post-stack iterative modeling migration and inversion (IMMI) Gary F. Margrave SUMMARY Post-stack IMMI The possibility of a post-stack process resembling full-waveform inversion (FWI) is investigated. As

More information

SUMMARY LEAST-SQUARES MIGRATION THEORY

SUMMARY 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 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 High fidelity imaging using reflections, refractions and multiples: North Sea example, Johan Sverdrup Field Grunde Rønholt*, Øystein Korsmo, Barbara Danielsen, Samuel Brown, Sverre Brandsberg-Dahl, Alejandro

More information

Residual move-out analysis with 3-D angle-domain common-image gathers

Residual 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 information

U043 3D Prestack Time Domain Full Waveform Inversion

U043 3D Prestack Time Domain Full Waveform Inversion U043 3D Prestack Time Domain Full Waveform Inversion D.V. Vigh* (WesternGeco), W.E.S. Starr (WesternGeco) & K.D. Kenneth Dingwall (WesternGeco) SUMMARY Despite the relatively high computational demand,

More information

The determination of the correct

The determination of the correct SPECIAL High-performance SECTION: H i gh-performance computing computing MARK NOBLE, Mines ParisTech PHILIPPE THIERRY, Intel CEDRIC TAILLANDIER, CGGVeritas (formerly Mines ParisTech) HENRI CALANDRA, Total

More information

Inversion after depth imaging

Inversion 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 information

Imaging with multiples using LSRTM

Imaging with multiples using LSRTM Chapter 4 Imaging with multiples using LSRTM In this chapter, I present a technique for imaging both primaries and higher-order multiples using joint least-squares reverse-time migration (joint-lsrtm).

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

E044 Ray-based Tomography for Q Estimation and Q Compensation in Complex Media

E044 Ray-based Tomography for Q Estimation and Q Compensation in Complex Media E044 Ray-based Tomography for Q Estimation and Q Compensation in Complex Media M. Cavalca* (WesternGeco), I. Moore (WesternGeco), L. Zhang (WesternGeco), S.L. Ng (WesternGeco), R.P. Fletcher (WesternGeco)

More information

Application of Weighted Early-Arrival Waveform Inversion to Shallow Land Data

Application of Weighted Early-Arrival Waveform Inversion to Shallow Land Data Application of Weighted Early-Arrival Waveform Inversion to Shallow Land Data Han Yu, Dongliang Zhang and Xin Wang King Abdullah University of Science and Technology, Division of Physical Science and Engineering,

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

Geometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT. August Lau and Chuan Yin.

Geometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT. August Lau and Chuan Yin. Geometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT August Lau and Chuan Yin January 6, 2017 Abstract The goal of seismic processing is to convert input data

More information

On the Scattering Effect of Lateral Discontinuities on AVA Migration

On the Scattering Effect of Lateral Discontinuities on AVA Migration On the Scattering Effect of Lateral Discontinuities on AVA Migration Juefu Wang* and Mauricio D. Sacchi Department of Physics, University of Alberta, 4 Avadh Bhatia Physics Laboratory, Edmonton, AB, T6G

More information

Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers

Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers J.W.C. Sherwood* (PGS), K. Sherwood (PGS), H. Tieman (PGS), R. Mager (PGS) & C. Zhou

More information

Application of time-lapse full waveform inversion of vertical seismic profile data for the identification of changes introduced by CO2 sequestration

Application of time-lapse full waveform inversion of vertical seismic profile data for the identification of changes introduced by CO2 sequestration Application of time-lapse full waveform inversion of vertical seismic profile data for the identification of changes introduced by CO2 sequestration Anton Egorov* Andrej Bόna Roman Pevzner Stanislav Glubokovskikh

More information

GEOPHYS 242: Near Surface Geophysical Imaging. Class 5: Refraction Migration Methods Wed, April 13, 2011

GEOPHYS 242: Near Surface Geophysical Imaging. Class 5: Refraction Migration Methods Wed, April 13, 2011 GEOPHYS 242: Near Surface Geophysical Imaging Class 5: Refraction Migration Methods Wed, April 13, 2011 Migration versus tomography Refraction traveltime and wavefield migration The theory of interferometry

More information

Static corrections in challenging cases

Static corrections in challenging cases P-146 Static corrections in challenging cases Seweryn Tlalka, Geofizyka Torun, Poland Summary Seismic data processing in challenging areas such as sand or rock desserts of Asia, North Africa, or Middle

More information

PS wave AVO aspects on processing, inversion, and interpretation

PS wave AVO aspects on processing, inversion, and interpretation PS wave AVO aspects on processing, inversion, and interpretation Yong Xu, Paradigm Geophysical Canada Summary Using PS wave AVO inversion, density contrasts can be obtained from converted wave data. The

More information

Chapter 5. 3D data examples REALISTICALLY COMPLEX SYNTHETIC INVERSION. Modeling generation and survey design

Chapter 5. 3D data examples REALISTICALLY COMPLEX SYNTHETIC INVERSION. Modeling generation and survey design Chapter 5 3D data examples In this chapter I will demonstrate the e ectiveness of the methodologies developed in the previous chapters using 3D data examples. I will first show inversion results of a realistically

More information

Wave-equation MVA applied to 4-D seismic monitoring

Wave-equation MVA applied to 4-D seismic monitoring Stanford Exploration Project, Report 112, November 11, 2002, pages 15 21 Short Note Wave-equation MVA applied to 4-D seismic monitoring Paul Sava, John Etgen, and Leon Thomsen 1 INTRODUCTION 4-D seismic

More information

The impact of the acoustic approximation on time-lapse FWI Bram Willemsen, M.I.T, Jun Cao, ConocoPhillips and Baishali Roy, ConocoPhillips

The impact of the acoustic approximation on time-lapse FWI Bram Willemsen, M.I.T, Jun Cao, ConocoPhillips and Baishali Roy, ConocoPhillips Bram Willemsen, M.I.T, Jun Cao, ConocoPhillips and Baishali Roy, ConocoPhillips SUMMARY Conventional time-lapse analysis methods are adequate for relatively simple geologies, but may degrade in the presence

More information

A MAP Algorithm for AVO Seismic Inversion Based on the Mixed (L 2, non-l 2 ) Norms to Separate Primary and Multiple Signals in Slowness Space

A MAP Algorithm for AVO Seismic Inversion Based on the Mixed (L 2, non-l 2 ) Norms to Separate Primary and Multiple Signals in Slowness Space A MAP Algorithm for AVO Seismic Inversion Based on the Mixed (L 2, non-l 2 ) Norms to Separate Primary and Multiple Signals in Slowness Space Harold Ivan Angulo Bustos Rio Grande do Norte State University

More information

Headwave Stacking in Terms of Partial Derivative Wavefield

Headwave Stacking in Terms of Partial Derivative Wavefield Geosystem Engineering, 7(1), 21-26 (March 2004) Headwave Stacking in Terms of Partial Derivative Wavefield Changsoo Shin School of Civil, Urban and Geosystem Engineering, Seoul National University, San

More information

UNIVERSITY OF CALGARY. Iterative Multiparameter Elastic Waveform Inversion Using Prestack Time Imaging and Kirchhoff approximation.

UNIVERSITY OF CALGARY. Iterative Multiparameter Elastic Waveform Inversion Using Prestack Time Imaging and Kirchhoff approximation. UNIVERSITY OF CALGARY Iterative Multiparameter Elastic Waveform Inversion Using Prestack Time Imaging and Kirchhoff approximation by Hassan Khaniani A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES

More information

Full-waveform inversion using seislet regularization a

Full-waveform inversion using seislet regularization a 1 Full-waveform inversion using seislet regularization a a Published in Geophysics, 82, no. 5, A43-A49, (2017) Zhiguang Xue 1, Hejun Zhu 2 and Sergey Fomel 1 ABSTRACT Because of inaccurate, incomplete

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

Which parameterization for acoustic vertical transverse. isotropic full waveform inversion? - Part 2: synthetic and

Which parameterization for acoustic vertical transverse. isotropic full waveform inversion? - Part 2: synthetic and Which parameterization for acoustic vertical transverse isotropic full waveform inversion? - Part : synthetic and real data case studies from Valhall Yaser Gholami, Romain Brossier, Stéphane Operto, Vincent

More information

An overview of full-waveform inversion in exploration geophysics

An overview of full-waveform inversion in exploration geophysics GEOPHYSICS, VOL. 7, NO. 6 NOVEMBER-DECEMBER 9 ; P. WCC7 WCC5, 5 FIGS., TABLE..9/.67 An overview of full-waveform inversion in exploration geophysics J. Virieux and S. Operto ABSTRACT Full-waveform inversion

More information

3D image-domain wavefield tomography using time-lag extended images

3D 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 information

AVO Analysis with Multi-Offset VSP Data

AVO Analysis with Multi-Offset VSP Data AVO Analysis with Multi-Offset VSP Data Z. Li*, W. Qian, B. Milkereit and E. Adam University of Toronto, Dept. of Physics, Toronto, ON, M5S 2J8 zli@physics.utoronto.ca T. Bohlen Kiel University, Geosciences,

More information

Tu A11 01 Seismic Re-processing and Q-FWI Model Building to Optimise AVO and Resolution in the Shallow Wisting Discovery

Tu A11 01 Seismic Re-processing and Q-FWI Model Building to Optimise AVO and Resolution in the Shallow Wisting Discovery Tu A11 01 Seismic Re-processing and Q-FWI Model Building to Optimise AVO and Resolution in the Shallow Wisting Discovery G. Apeland* (WesternGeco), P. Smith (WesternGeco), O. Lewis (WesternGeco), S. Way

More information

P063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA

P063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA 1 P063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA I. TRINKS 1, S. SINGH 2, C. CHAPMAN 3, P. BARTON 1, M. BOSCH 4, A. CHERRETT 5 Addresses 1 Bullard Laboratories,

More information

Directions in 3-D imaging - Strike, dip, both?

Directions 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 information

Target-oriented wave-equation inversion with regularization in the subsurface-offset domain

Target-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 information

Full-model wavenumber inversion: An emphasis on the appropriate wavenumber continuation

Full-model wavenumber inversion: An emphasis on the appropriate wavenumber continuation GEOPHYSICS, VOL. 81, NO. 3 (MAY-JUNE 2016); P. R89 R98, 12 FIGS. 10.1190/GEO2015-0537.1 Full-model wavenumber inversion: An emphasis on the appropriate wavenumber continuation Tariq Alkhalifah 1 ABSTRACT

More information

Marmousi synthetic dataset

Marmousi synthetic dataset Stanford Exploration Project, Report DATALIB, January 16, 2002, pages 1?? Marmousi synthetic dataset Carmen B. Mora 1 ABSTRACT Marmousi is a 2-D synthetic dataset generated at the Institut Français du

More information

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

Downloaded 11/13/12 to Redistribution subject to SEG license or copyright; see Terms of Use at VTI Waveform Inversion with Practical Strategies: Application to 3D Real Data Chao Wang, David Yingst, Robert Bloor, and Jacques Leveille, ION Geophysical SUMMARY The term of Waveform inversion (WFI) refers

More information

SUMMARY. method to synthetic datasets is discussed in the present paper.

SUMMARY. method to synthetic datasets is discussed in the present paper. Geophysical modeling through simultaneous Joint Inversion of Seismic, Gravity and Magnetotelluric data Michele De Stefano (1), Daniele Colombo (1) WesternGeco EM - Geosystem, via Clericetti 42/A, 20133

More information

=, (1) SEG/New Orleans 2006 Annual Meeting

=, (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 information

Analysis of the traveltime sensitivity kernels for an acoustic transversely isotropic medium with a vertical axis of symmetry

Analysis of the traveltime sensitivity kernels for an acoustic transversely isotropic medium with a vertical axis of symmetry Geophysical Prospecting, 2017, 65, 22 34 doi: 10.1111/1365-2478.12361 Analysis of the traveltime sensitivity kernels for an acoustic transversely isotropic medium with a vertical axis of symmetry Ramzi

More information

Elastic full waveform Inversion for land walkaway VSP data

Elastic full waveform Inversion for land walkaway VSP data Elastic full waveform Inversion for land walkaway VSP data Olga Podgornova, Scott Leaney, arwan Charara, Schlumberger; Eric von Lunen, Nexen Energy ULC Summary Full waveform inversion (FWI) is capable

More information

Well-log derived step lengths in Full Waveform Inversion Tunde Arenrin Gary Margrave

Well-log derived step lengths in Full Waveform Inversion Tunde Arenrin Gary Margrave Well-log derived step lengths in Full Waveform Inversion Tunde Arenrin Gary Margrave Motivation»Being able to make FWI a part of the seismic processing work flow under current computational capability.»incorporating

More information

Wave-equation migration from topography: Imaging Husky

Wave-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

Provides New Tech Solution For Seismic Velocity Modeling

Provides New Tech Solution For Seismic Velocity Modeling JNURY 2009 The etter usiness Publication Serving the Exploration / Drilling / Production Industry Provides New Tech Solution For Seismic Velocity Modeling y Moritz M. Fliedner & Dimitri evc SNT CLR, C.

More information

DSR Migration Velocity Analysis by Differential Semblance Optimization

DSR Migration Velocity Analysis by Differential Semblance Optimization 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

More information

Full waveform inversion by deconvolution gradient method

Full 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 information

Automatic wave equation migration velocity analysis by differential semblance optimization

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 information

Coherent noise attenuation using Inverse Problems and Prediction Error Filters

Coherent noise attenuation using Inverse Problems and Prediction Error Filters Stanford Exploration Project, Report 105, September 5, 2000, pages 1 49 Coherent noise attenuation using Inverse Problems and Prediction Error Filters Antoine Guitton 1 ABSTRACT Two iterative methods that

More information