Education Ph.D. Candidate, Rice University, Houston, TX, USA, 08/ Present. M.S, Shanghai Jiao Tong University, Shanghai, China, 09/ /2009
|
|
- Ernest Benedict Jones
- 5 years ago
- Views:
Transcription
1 Yin Huang Education Ph.D. Candidate, Rice University, Houston, TX, USA, 08/ Present Dissertation Topic: Nonlinear Extended Waveform Inversion M.A. with Master Thesis: Transparency property of one dimensional acoustic wave equations Relevant courses: Numerical Differential Equations; Optimization; High Performance Computing; Geophysical Data Analysis M.S, Shanghai Jiao Tong University, Shanghai, China, 09/ /2009 Dissertation topic: Comparision of many numerical methods for saddle point system arising from the mixed finite element method of elliptic problems with nonsmooth coefficients Research Interests Seismic waveform inversion Forward and inverse problems for non-homogeneous medium High performace computing 1
2 Linearized Extended Waveform Inversion: Reduced Objective Function and Its Gradient Yin Huang Advisor: William Symes The Rice Inversion Project Annual Review Meeting April 18,
3 Introduction Seismic reflection inversion Over-determined: highly redundant in the observed data; (Gauthier et al., 1986; Santosa & Symes, 1989; Virieux & Operto, 2009) Extended model fitting Under-determined: model has more degree of freedom than data. Data misfit + extended modeling + differential semblance smooth objective function of model parameter. (Symes & Kern 1994, Symes, 1999, Stolk & Symes 2003) 3
4 Introduction Seismic reflection inversion Over-determined: highly redundant in the observed data; (Gauthier et al., 1986; Santosa & Symes, 1989; Virieux & Operto, 2009) Extended model fitting Under-determined: model has more degree of freedom than data. Data misfit + extended modeling + differential semblance smooth objective function of model parameter. (Symes & Kern 1994, Symes, 1999, Stolk & Symes 2003) 3
5 Introduction Seismic reflection inversion Over-determined: highly redundant in the observed data; (Gauthier et al., 1986; Santosa & Symes, 1989; Virieux & Operto, 2009) Extended model fitting Under-determined: model has more degree of freedom than data. Data misfit + extended modeling + differential semblance smooth objective function of model parameter. (Symes & Kern 1994, Symes, 1999, Stolk & Symes 2003) 3
6 Waveform Inversion and Extended Modeling Extended Waveform Inversion: Given data d D, find m M so that (Symes, 1986, 1991; Biondi and Almomin, 2012) F [ m] d. M = {m(x)} physical model space: velocity, density, bulk modulus,... M = {m(x, h)} extended model space, M M. D data space. F : M D forward map: acoustic, elastic... F : M D extended forward map. Extended model separation: m m l + δ m. Background model is physical. Reflectivity is extended. 4
7 Waveform Inversion and Extended Modeling Extended Waveform Inversion: Given data d D, find m M so that (Symes, 1986, 1991; Biondi and Almomin, 2012) F [ m] d. M = {m(x)} physical model space: velocity, density, bulk modulus,... M = {m(x, h)} extended model space, M M. D data space. F : M D forward map: acoustic, elastic... F : M D extended forward map. Extended model separation: m m l + δ m. Background model is physical. Reflectivity is extended. 4
8 Linearized Inversion Given data d D, find m l, δ m so that D F [m l ]δ m d F [m l ]. D F is the derivative, or Born approximation. (Symes and Carazzone, 1991; Chauris and Noble, 2001; Mulder and ten Kroode, 2002; Shen and Symes, 2008; Symes 2008.) Objective function: J[m l, δ m] = 1 2 D F [m l ]δ m (d F [m l ]) 2 + α2 Aδ m 2 2 A annihilator, Aδm = 0 for all δm M. 5
9 Linearized Inversion Given data d D, find m l, δ m so that D F [m l ]δ m d F [m l ]. D F is the derivative, or Born approximation. (Symes and Carazzone, 1991; Chauris and Noble, 2001; Mulder and ten Kroode, 2002; Shen and Symes, 2008; Symes 2008.) Objective function: J[m l, δ m] = 1 2 D F [m l ]δ m (d F [m l ]) 2 + α2 Aδ m 2 2 A annihilator, Aδm = 0 for all δm M. 5
10 Reduced Objective Function The value of the reduced objective function is the minimum value of J for fixed m l J[m l ] = min J[m l, δ m]. δ m Let Minimizer of J N[m l ] = D F [m l ] T D F [m l ] + α 2 A T A δ m[m l ] = N[m l ] 1 D F [m l ] T (d F [m l ]). 6
11 Gradient of J[m l ] Directional derivative D J[m l ]dm l = D ml J[m l, δ m]dm l + D δ m J[m l, δ m]d ml δ mdm l. Second term is 0, if δ m is solved exactly. Direct computation gives gradient J[m l ] = D 2 F [m l ] T [δ m, D F [m l ]δ m (d F [m l ])] +DF [m l ] T (D F [m l ]δ m (d F [m l ])). Second term is zero if apply an appropriate cutoff function to the data. D 2 F [ml ] T : M D such that for dml M, q M and φ D D 2 F [m l ][dm l, q], φ = dm l, D 2 F [m l ] T [q, φ]. 7
12 Gradient of J[m l ] Directional derivative D J[m l ]dm l = D ml J[m l, δ m]dm l + D δ m J[m l, δ m]d ml δ mdm l. Second term is 0, if δ m is solved exactly. Direct computation gives gradient J[m l ] = D 2 F [m l ] T [δ m, D F [m l ]δ m (d F [m l ])] +DF [m l ] T (D F [m l ]δ m (d F [m l ])). Second term is zero if apply an appropriate cutoff function to the data. D 2 F [m l ] T : M D such that for dm l M, q M and φ D D 2 F [m l ][dm l, q], φ = dm l, D 2 F [m l ] T [q, φ]. 7
13 Example: Computation of D 2 F [ml ] T For acoustic constant density wave equation: q M, φ D. Let ω solve the zero final value problem ( 1 2 ) t 2 ω = dx r φδ(x x r ). m 2 l D F [m l ] T φ = 2 m l dtω u. Let ω 0 solve the zero final value problem ( 1 2 ) ( ) q t 2 ω 0 = 2 ω. m l m 2 l D 2 F T [m l ][q, φ] = 2 dm l m l dh with ū the Born approximation wave field. dt (( ū + q ) ) u ω + ω 0 u. m l 8
14 Example: Computation of D 2 F [ml ] T For acoustic constant density wave equation: q M, φ D. Let ω solve the zero final value problem ( 1 2 ) t 2 ω = dx r φδ(x x r ). m 2 l D F [m l ] T φ = 2 m l dtω u. Let ω 0 solve the zero final value problem ( 1 2 ) ( ) q t 2 ω 0 = 2 ω. m l m 2 l D 2 F T [m l ][q, φ] = 2 dm l m l dh with ū the Born approximation wave field. dt (( ū + q ) ) u ω + ω 0 u. m l 8
15 Gradient needs a correction term Why need correction term? δ m J[m l, δ m] is not zero numerically. D ml δ m could be large. Neglect of the second term leads to large errors (Symes and Kern, 1994). Iteration method is used to solve for δ m, P k polynomial with degree k. δ m k [m l ] = P k (N[m l ])D F [m l ] T (d F [m l ]). D ml δ m k [m l ]dm l = D(P k (N[m l ]))dm l D F [m l ] T (d F [m l ]) +P k (N[m l ])D 2 F [ml ] T [dm l, d F [m l ]] P k (N[m l ])D F [m l ] T DF [m l ]dm l. 9
16 Gradient needs a correction term Why need correction term? δ m J[m l, δ m] is not zero numerically. D ml δ m could be large. Neglect of the second term leads to large errors (Symes and Kern, 1994). Iteration method is used to solve for δ m, P k polynomial with degree k. δ m k [m l ] = P k (N[m l ])D F [m l ] T (d F [m l ]). D ml δ m k [m l ]dm l = D(P k (N[m l ]))dm l D F [m l ] T (d F [m l ]) +P k (N[m l ])D 2 F [ml ] T [dm l, d F [m l ]] P k (N[m l ])D F [m l ] T DF [m l ]dm l. 9
17 Chebyshev iteration Chebyshev is preferred. Coefficients of P k does not depend on N[m l ], once spectrum bound of N[m l ] is known in advance. Easy to analyze D ml δ m k [m l ]. Coefficients in conjugate gradient depends on N[m l ]. Bound on derivative is not obvious to calculate. For operator A, want P k (A) A 1 I AP k (A) 0 Fix k, find P k with leading coefficient 1, such that P k = arg min p k max 1 λp k(λ) λ with λ [λ min, λmax], i.e. spectrum bound of operator A. NOTE: Chebyshev polynomial is the unique optimal solution. (Varga 1962) 10
18 Chebyshev iteration Chebyshev is preferred. Coefficients of P k does not depend on N[m l ], once spectrum bound of N[m l ] is known in advance. Easy to analyze D ml δ m k [m l ]. Coefficients in conjugate gradient depends on N[m l ]. Bound on derivative is not obvious to calculate. For operator A, want P k (A) A 1 I AP k (A) 0 Fix k, find P k with leading coefficient 1, such that P k = arg min p k max 1 λp k(λ) λ with λ [λ min, λmax], i.e. spectrum bound of operator A. NOTE: Chebyshev polynomial is the unique optimal solution. (Varga 1962) 10
19 Gradient with a correction term P k (N[m l]) is of order k 2. D(P k (N[m l ])) D ml N[m l ]P k (N[m l]) is of order k 2. D(P k (N[m l ]))r k = O(k 2 r k ). Normal residual r k = N[m l ]δ m k D F [m l ] T (d F [m l ]) decrease exponentially in k. D δ m J[m l, δ m]d ml δ mdm l = D ml (P k (N[m l ]))dm l D F [m l ] T (d F [m l ]), r k + P k (N[m l ])D 2 F [ml ] T [dm l, d F [m l ]], r k dm l, D 2 F [ml ] T [P k (N[m l ])r k, d F [m l ]]. Gradient with correction term J[m l ] D 2 F [ml ] T [δ m, D F [m l ]δ m (d F [m l ])] +D 2 F [m l ] T [P k (N[m l ])r k, d F [m l ]]. 11
20 Gradient with a correction term P k (N[m l]) is of order k 2. D(P k (N[m l ])) D ml N[m l ]P k (N[m l]) is of order k 2. D(P k (N[m l ]))r k = O(k 2 r k ). Normal residual r k = N[m l ]δ m k D F [m l ] T (d F [m l ]) decrease exponentially in k. D δ m J[m l, δ m]d ml δ mdm l = D ml (P k (N[m l ]))dm l D F [m l ] T (d F [m l ]), r k + P k (N[m l ])D 2 F [ml ] T [dm l, d F [m l ]], r k dm l, D 2 F [ml ] T [P k (N[m l ])r k, d F [m l ]]. Gradient with correction term J[m l ] D 2 F [ml ] T [δ m, D F [m l ]δ m (d F [m l ])] +D 2 F [m l ] T [P k (N[m l ])r k, d F [m l ]]. 11
21 Restarting Chebyshev iteration Simple scheme: Compute no. of iterations k and Chebyshev coefficients Initialize the iteration, compute Rayleigh Quotient (RQ) to get spectrum bound λ Run Chebyshev iteration until Get bigger λ, restart Chebyshev Or run k steps and return. NOTE: Performance depends on condition number. Good preconditioner is needed. Implementation is available in RVL. Restarting curve for depth-oriented extended waveform inversion (Liu et al. 2013). 12
22 Summary Reduced objective function Smooth in long scale model D 2 F [ml ] T of extended acoustic constant density forward operator Gradient Correction term to the gradient Chebyshev iteration Better than CG for operator with moderate condition number and uniformly distributed spectrum Restarting scheme for operators without prior information of spectrum Available in RVL 13
23 Future work Short term Fill in detailed mathematical analysis Find good preconditioner for shot-coordinate extension Long term Extend to nonlinear waveform inversion (D. Sun s PhD thesis) 14
24 Acknowledgements Great thanks to Current and former TRIP team members. Sponsors of The Rice Inversion Project. Wonderful audience. 15
SUMMARY. (Sun and Symes, 2012; Biondi and Almomin, 2012; Almomin and Biondi, 2012).
Inversion Velocity Analysis via Differential Semblance Optimization in the Depth-oriented Extension Yujin Liu, William W. Symes, Yin Huang,and Zhenchun Li, China University of Petroleum (Huadong), Rice
More informationYin Huang. Education. Research Interests. Rice University, Houston, TX, USA. Shanghai Jiao Tong University, Shanghai, China
Yin Huang Education Rice University, Houston, TX, USA Ph.D. Candidate, Computational and Applied Mathematics 08/2010 - Present Expected Dissertation Topic: Inverse problems for wave propagation for nonsmooth
More informationSource function estimation in extended full waveform inversion
Source function estimation in extended full waveform inversion Lei Fu The Rice Inversion Project (TRIP) May 6, 2014 Lei Fu (TRIP) Source estimation in EFWI May 6, 2014 1 Overview Objective Recover Earth
More informationA nonlinear differential-semblance algorithm for waveform inversion
A nonlinear differential-semblance algorithm for waveform inversion Dong Sun The Rice Inversion Project Annual Review Meeting March 3, 212, 1/15 Agenda 1 Full Waveform Inversion 2 Extension and Differential
More informationA Nonlinear Differential Semblance Strategy for Waveform Inversion: Experiments in Layered Media Dong Sun and William W Symes, Rice University
A Nonlinear Differential Semblance Strategy for Waveform Inversion: Experiments in Layered Media Dong Sun and William W Symes, Rice University SUMMARY This paper proposes an alternative approach to the
More informationOverview Avoiding Cycle Skipping: Model Extension MSWI: Space-time Extension Numerical Examples
Overview Avoiding Cycle Skipping: Model Extension MSWI: Space-time Extension Numerical Examples Guanghui Huang Education University of Chinese Academy of Sciences, Beijing, China Ph.D. in Computational
More informationComparison between least-squares reverse time migration and full-waveform inversion
Comparison between least-squares reverse time migration and full-waveform inversion Lei Yang, Daniel O. Trad and Wenyong Pan Summary The inverse problem in exploration geophysics usually consists of two
More informationObjective functions for full waveform inversion
Objective functions for full waveform inversion William Symes The Rice Inversion Project EAGE 2012 Workshop: From Kinematic to Waveform Inversion A Tribute to Patrick Lailly Agenda Overview Challenges
More informationFull Waveform Inversion via Matched Source Extension
Full Waveform Inversion via Matched Source Extension Guanghui Huang and William W. Symes TRIP, Department of Computational and Applied Mathematics May 1, 215, TRIP Annual Meeting G. Huang and W. W. Symes
More informationTHE RICE INVERSION PROJECT
THE RICE INVERSION PROJECT Mario Bencomo, Lei Fu, Jie Hou, Guanghui Huang, Rami Nammour, and William Symes Annual Report 6 Copyright 5-6 by Rice University i TRIP6 TABLE OF CONTENTS William W. Symes, The
More informationThe Search for a Cycle-Skipping Cure: an Overview
The Search for a Cycle-Skipping Cure: an Overview William Symes Computational and Applied Mathematics & Earth Science, Rice University The Rice Inversion Project IPAM, May 2017 Agenda Fixes within FWI
More informationComments on the Geophysics paper Multiparameter 11norm. waveform fitting: Interpretation of Gulf of Mexico reflection
,..-., Comments on the Geophysics paper Multiparameter 11norm waveform fitting: Interpretation of Gulf of Mexico reflection seismograms by H. Djikp&s6 and A. Tarantola Susan E. Minkoff * t and William
More informationElastic full waveform inversion for near surface imaging in CMP domain Zhiyang Liu*, Jie Zhang, University of Science and Technology of China (USTC)
Elastic full waveform inversion for near surface imaging in CMP domain Zhiyang Liu*, Jie Zhang, University of Science and Technology of China (USTC) Summary We develop an elastic full waveform inversion
More informationAn Approximate Inverse to the Extended Born Modeling Operator Jie Hou, William W. Symes, The Rice Inversion Project, Rice University
Jie Hou, William W. Symes, The Rice Inversion Project, Rice University SUMMARY We modify RTM to create an approximate inverse to the extended Born modeling operator in 2D. The derivation uses asymptotic
More informationAbstract. A Nonlinear Differential Semblance Algorithm for Waveform Inversion. Dong Sun
2 Abstract A Nonlinear Differential Semblance Algorithm for Waveform Inversion by Dong Sun This thesis proposes a nonlinear differential semblance approach to full waveform inversion as an alternative
More informationWenyong Pan and Lianjie Huang. Los Alamos National Laboratory, Geophysics Group, MS D452, Los Alamos, NM 87545, USA
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 019 SGP-TR-14 Adaptive Viscoelastic-Waveform Inversion Using the Local Wavelet
More informationAccelerating. Extended Least Squares Migration. Jie Hou TRIP 2014 Review Meeting May 1st, 2015
Accelerating Extended Least Squares Migration Jie Hou TRIP 2014 Review Meeting May 1st, 2015 Linearized Inverse Problem Fm = d Adjoint Operator m = F T d Approximate Inverse m = F d 2 Linearized Inverse
More informationIntroduction to Seismic Imaging
Introduction to Seismic Imaging Alison Malcolm Department of Earth, Atmospheric and Planetary Sciences MIT August 20, 2010 Outline Introduction Why we image the Earth How data are collected Imaging vs
More informationAttenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation
Attenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation Gaurav Dutta, Kai Lu, Xin Wang and Gerard T. Schuster, King Abdullah University of Science and Technology
More informationReverse Time Shot-Geophone Migration and
Reverse Time Shot-Geophone Migration and velocity analysis, linearizing, inverse problems, extended models, etc. etc. William W. Symes PIMS Geophysical Inversion Workshop Calgary, Alberta July 2003 www.trip.caam.rice.edu
More informationWave-equation tomography for anisotropic parameters
Wave-equation tomography for anisotropic parameters Yunyue (Elita) Li and Biondo Biondi ABSTRACT Anisotropic models are recognized as more realistic representations of the subsurface where complex geological
More informationFull-Waveform Inversion with Gauss- Newton-Krylov Method
Full-Waveform Inversion with Gauss- Newton-Krylov Method Yogi A. Erlangga and Felix J. Herrmann {yerlangga,fherrmann}@eos.ubc.ca Seismic Laboratory for Imaging and Modeling The University of British Columbia
More informationA 27-point scheme for a 3D frequency-domain scalar wave equation based on an average-derivative method
Geophysical Prospecting 04 6 58 77 doi: 0./365-478.090 A 7-point scheme for a 3D frequency-domain scalar wave equation based on an average-derivative method Jing-Bo Chen Key Laboratory of Petroleum Resources
More informationA projected Hessian for full waveform inversion
CWP-679 A projected Hessian for full waveform inversion Yong Ma & Dave Hale Center for Wave Phenomena, Colorado School of Mines, Golden, CO 80401, USA (c) Figure 1. Update directions for one iteration
More informationSUMMARY INTRODUCTION. f ad j (t) = 2 Es,r. The kernel
The failure mode of correlation focusing for model velocity estimation Hyoungsu Baek 1(*), Henri Calandra 2, and Laurent Demanet 1 1 Dept. of Mathematics and Earth Resources Lab, Massachusetts Institute
More informationEstimating a pseudounitary operator for velocity-stack inversion
Stanford Exploration Project, Report 82, May 11, 2001, pages 1 77 Estimating a pseudounitary operator for velocity-stack inversion David E. Lumley 1 ABSTRACT I estimate a pseudounitary operator for enhancing
More informationExtensions and Nonlinear Inverse Scattering
Extensions and Nonlinear Inverse Scattering William W. Symes CAAM Colloquium, September 2004 Data parameters: time t, source location x s, and receiver location x r, (vector) half offset h = x r x s, scalar
More informationModel Extensions and Inverse Scattering: Inversion for Seismic Velocities
Model Extensions and Inverse Scattering: Inversion for Seismic Velocities William W. Symes Rice University October 2007 William W. Symes ( Rice University) Model Extensions and Inverse Scattering: Inversion
More informationFull waveform inversion and the inverse Hessian
Full waveform inversion and the inverse Hessian Gary Margrave, Matt Yedlin and Kris Innanen ABSTRACT FWI and the inverse Hessian Full waveform inversion involves defining an objective function, and then
More informationSource estimation for frequency-domain FWI with robust penalties
Source estimation for frequency-domain FWI with robust penalties Aleksandr Y. Aravkin, Tristan van Leeuwen, Henri Calandra, and Felix J. Herrmann Dept. of Earth and Ocean sciences University of British
More informationSUMMARY REVIEW OF THE FREQUENCY DOMAIN L2 FWI-HESSIAN
Efficient stochastic Hessian estimation for full waveform inversion Lucas A. Willemsen, Alison E. Malcolm and Russell J. Hewett, Massachusetts Institute of Technology SUMMARY In this abstract we present
More informationTwo-Scale Wave Equation Modeling for Seismic Inversion
Two-Scale Wave Equation Modeling for Seismic Inversion Susan E. Minkoff Department of Mathematics and Statistics University of Maryland Baltimore County Baltimore, MD 21250, USA RICAM Workshop 3: Wave
More informationImage-space wave-equation tomography in the generalized source domain
Image-space wave-equation tomography in the generalized source domain Yaxun Tang, Claudio Guerra, and Biondo Biondi ABSTRACT We extend the theory of image-space wave-equation tomography to the generalized
More informationRegistration-guided least-squares waveform inversion
Registration-guided least-squares waveform inversion Hyoungsu Baek 1, Henri Calandra, Laurent Demanet 1 1 MIT Mathematics department, TOTAL S.A. January 15 013 Abstract Full waveform inversion with frequency
More informationModeling, migration, and inversion in the generalized source and receiver domain
Modeling, migration, and inversion in the generalized source and receiver domain Yaxun Tang ABSTRACT I extend the theory of Born modeling, migration and inversion to the generalized source and receiver
More information1.2 Derivation. d p f = d p f(x(p)) = x fd p x (= f x x p ). (1) Second, g x x p + g p = 0. d p f = f x g 1. The expression f x gx
PDE-constrained optimization and the adjoint method Andrew M. Bradley November 16, 21 PDE-constrained optimization and the adjoint method for solving these and related problems appear in a wide range of
More informationMigration with Implicit Solvers for the Time-harmonic Helmholtz
Migration with Implicit Solvers for the Time-harmonic Helmholtz Yogi A. Erlangga, Felix J. Herrmann Seismic Laboratory for Imaging and Modeling, The University of British Columbia {yerlangga,fherrmann}@eos.ubc.ca
More informationApproximate- vs. full-hessian in FWI: 1D analytical and numerical experiments
Approximate- vs. full-hessian in FWI: 1D analytical and numerical experiments Raul Cova and Kris Innanen ABSTRACT Feasibility of using Full Waveform Inversion (FWI) to build velocity models has been increasing
More informationElastic least-squares migration with two-way wave equation forward and adjoint operators
Elastic least-squares migration with two-way wave equation forward and adjoint operators Ke Chen and Mauricio D. Sacchi, Department of Physics, University of Alberta Summary Time domain elastic least-squares
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
More informationFull waveform inversion of shot gathers in terms of poro-elastic parameters
Full waveform inversion of shot gathers in terms of poro-elastic parameters Louis De Barros, M. Dietrich To cite this version: Louis De Barros, M. Dietrich. Full waveform inversion of shot gathers in terms
More informationFWI with Compressive Updates Aleksandr Aravkin, Felix Herrmann, Tristan van Leeuwen, Xiang Li, James Burke
Consortium 2010 FWI with Compressive Updates Aleksandr Aravkin, Felix Herrmann, Tristan van Leeuwen, Xiang Li, James Burke SLIM University of British Columbia Full Waveform Inversion The Full Waveform
More informationElastic least-squares reverse time migration
CWP-865 Elastic least-squares reverse time migration Yuting Duan, Paul Sava, and Antoine Guitton Center for Wave Phenomena, Colorado School of Mines ABSTRACT Least-squares migration (LSM) can produce images
More informationImage amplitudes in reverse time migration/inversion
Image amplitudes in reverse time migration/inversion Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University TRIP Seminar, October 6, 2010
More informationDownloaded 11/08/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Velocity tomography based on turning waves and finite-frequency sensitivity kernels Xiao-Bi Xie*, University of California at Santa Cruz, Jan Pajchel, Statoil, Norway Summary We present a velocity tomography
More informationMathematics of Seismic Imaging Part 1: Modeling and Inverse Problems
Mathematics of Seismic Imaging Part 1: Modeling and Inverse Problems William W. Symes Rice University 0. Introduction How do you turn lots of this... 0 offset (km) -4-3 -2-1 1 2 time (s) 3 4 5 (field seismogram
More informationProgress and Prospects for Wave Equation MVA
Progress and Prospects for Wave Equation MVA William W. Symes TRIP Annual Review, January 007 Agenda MVA, Nonlinear Inversion, and Extended Modeling MVA with two different extensions: the coherent noise
More informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Inverse Scattering Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 30, 2014 Bee Bednar (Panorama Technologies)
More informationWavefield Reconstruction Inversion (WRI) a new take on wave-equation based inversion Felix J. Herrmann
Wavefield Reconstruction Inversion (WRI) a new take on wave-equation based inversion Felix J. Herrmann SLIM University of British Columbia van Leeuwen, T and Herrmann, F J (2013). Mitigating local minima
More informationUncertainty quantification for Wavefield Reconstruction Inversion
Uncertainty quantification for Wavefield Reconstruction Inversion Zhilong Fang *, Chia Ying Lee, Curt Da Silva *, Felix J. Herrmann *, and Rachel Kuske * Seismic Laboratory for Imaging and Modeling (SLIM),
More informationBureau of Economic Geology, The University of Texas at Austin, Austin, TX. Research Associate, Geophysics from 11/2004
Paul C. Sava Bureau of Economic Geology, The University of Texas at Austin, University Station, Box X, Austin, TX 78758, (512) 471-0327 paul.sava@beg.utexas.edu, http://sepwww.stanford.edu/sep/paul Research
More informationFull Waveform Inversion (FWI) with wave-equation migration. Gary Margrave Rob Ferguson Chad Hogan Banff, 3 Dec. 2010
Full Waveform Inversion (FWI) with wave-equation migration (WEM) and well control Gary Margrave Rob Ferguson Chad Hogan Banff, 3 Dec. 2010 Outline The FWI cycle The fundamental theorem of FWI Understanding
More informationAPPLICATION OF OPTIMAL BASIS FUNCTIONS IN FULL WAVEFORM INVERSION
METHODS AND APPLICATIONS OF ANALYSIS. c 2004 International Press Vol. 11, No. 3, pp. 345 352, September 2004 004 APPLICATION OF OPTIMAL BASIS FUNCTIONS IN FULL WAVEFORM INVERSION PING SHENG, GANG SUN,
More informationIntroduction to Seismic Imaging
Introduction to Seismic Imaging Alison Malcolm Department of Earth, Atmospheric and Planetary Sciences MIT August 20, 2010 Outline Introduction Why we image the Earth How data are collected Imaging vs
More informationSeismic inverse scattering by reverse time migration
Seismic inverse scattering by reverse time migration Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University MSRI Inverse Problems Seminar,
More informationVasil Khalidov & Miles Hansard. C.M. Bishop s PRML: Chapter 5; Neural Networks
C.M. Bishop s PRML: Chapter 5; Neural Networks Introduction The aim is, as before, to find useful decompositions of the target variable; t(x) = y(x, w) + ɛ(x) (3.7) t(x n ) and x n are the observations,
More information2C9 Design for seismic and climate changes. Jiří Máca
2C9 Design for seismic and climate changes Jiří Máca List of lectures 1. Elements of seismology and seismicity I 2. Elements of seismology and seismicity II 3. Dynamic analysis of single-degree-of-freedom
More informationNonlinear seismic imaging via reduced order model backprojection
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, Vladimir Druskin 2 and Mikhail Zaslavsky 2 University of Houston, 2 Schlumberger-Doll Research Center Mamonov, Druskin,
More informationMatrix Probing and Simultaneous Sources: A New Approach for Preconditioning the Hessian
Matrix Probing and Simultaneous Sources: A New Approach for Preconditioning the Hessian Curt Da Silva 1 and Felix J. Herrmann 2 1 Dept. of Mathematics 2 Dept. of Earth and Ocean SciencesUniversity of British
More informationExtended isochron rays in prestack depth (map) migration
Extended isochron rays in prestack depth (map) migration A.A. Duchkov and M.V. de Hoop Purdue University, 150 N.University st., West Lafayette, IN, 47907 e-mail: aduchkov@purdue.edu (December 15, 2008)
More informationWe E Multiparameter Full-waveform Inversion for Acoustic VTI Medium with Surface Seismic Data
We E16 4 Multiarameter Full-waveform Inversion for Acoustic VI Medium with Surface Seismic Data X. Cheng* (Schlumberger) K. Jiao (Schlumberger) D. Sun (Schlumberger) & D. Vigh (Schlumberger) SUMMARY In
More informationNon-linear least squares
Non-linear least squares Concept of non-linear least squares We have extensively studied linear least squares or linear regression. We see that there is a unique regression line that can be determined
More informationTieyuan Zhu Postdoctoral Fellow, Jackson School of Geosciences, the University of Texas at Austin Mail address: Telephone: Website:
Postdoctoral Fellow, Jackson School of Geosciences, the University of Texas at Austin Mail address: Telephone: Email: Website: 10100 Burnet Rd. #130, Austin TX 01-650- 308-6506 tzhu@jsg.utexas.edu http://www.jsg.utexas.edu/tyzhu/
More informationCAAM 454/554: Stationary Iterative Methods
CAAM 454/554: Stationary Iterative Methods Yin Zhang (draft) CAAM, Rice University, Houston, TX 77005 2007, Revised 2010 Abstract Stationary iterative methods for solving systems of linear equations are
More informationW011 Full Waveform Inversion for Detailed Velocity Model Building
W011 Full Waveform Inversion for Detailed Velocity Model Building S. Kapoor* (WesternGeco, LLC), D. Vigh (WesternGeco), H. Li (WesternGeco) & D. Derharoutian (WesternGeco) SUMMARY An accurate earth model
More informationGEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS
Methods in Geochemistry and Geophysics, 36 GEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS Michael S. ZHDANOV University of Utah Salt Lake City UTAH, U.S.A. 2OO2 ELSEVIER Amsterdam - Boston - London
More informationP016 Toward Gauss-Newton and Exact Newton Optimization for Full Waveform Inversion
P016 Toward Gauss-Newton and Exact Newton Optiization for Full Wavefor Inversion L. Métivier* ISTerre, R. Brossier ISTerre, J. Virieux ISTerre & S. Operto Géoazur SUMMARY Full Wavefor Inversion FWI applications
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationVollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich
Seminar des FA Ultraschallprüfung Vortrag 1 More info about this article: http://www.ndt.net/?id=20944 Vollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich Thomas
More informationarxiv: v1 [math.na] 1 Apr 2015
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, University of Houston; Vladimir Druskin and Mikhail Zaslavsky, Schlumberger arxiv:1504.00094v1 [math.na] 1 Apr 2015
More informationAn Immersed Boundary Method for Restricted Diffusion with Permeable Interfaces
An Immersed Boundary Method for Restricted Diffusion with Permeable Interfaces Huaxiong Huang Kazuyasu Sugiyama Shu Takagi April 6, 009 Keywords: Restricted diffusion; Permeable interface; Immersed boundary
More informationStanford Exploration Project, Report 115, May 22, 2004, pages
Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 248 Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 First-order lateral interval velocity estimates without picking
More informationSeismic Imaging. William W. Symes. C. I. M. E. Martina Franca September
Seismic Imaging William W. Symes C. I. M. E. Martina Franca September 2002 www.trip.caam.rice.edu 1 0 offset (km) -4-3 -2-1 1 2 time (s) 3 4 5 How do you turn lots of this... (field seismogram from the
More informationAcoustic wavefield imaging using the energy norm
CWP-831 Acoustic wavefield imaging using the energy norm Daniel Rocha 1, Nicolay Tanushev 2 & Paul Sava 1 1 Center for Wave Phenomena, Colorado School of Mines 2 Z-Terra, Inc ABSTRACT Wavefield energy
More informationWave-equation reflection tomography: annihilators and sensitivity kernels
Geophys. J. Int. (2006 167, 1332 1352 doi: 10.1111/j.1365-246X.2006.03132.x Wave-equation reflection tomography: annihilators and sensitivity kernels Maarten V. de Hoop, 1 Robert D. van der Hilst 2 and
More informationEstimation of the Residual Stiffness of Fire-Damaged Concrete Members
Copyright 2011 Tech Science Press CMC, vol.22, no.3, pp.261-273, 2011 Estimation of the Residual Stiffness of Fire-Damaged Concrete Members J.M. Zhu 1, X.C. Wang 1, D. Wei 2, Y.H. Liu 2 and B.Y. Xu 2 Abstract:
More informationA novel full waveform inversion method based on a timeshift nonlinear operator
A novel full waveform inversion method based on a timeshift nonlinear operator Journal: Manuscript ID GJI-S--0.R Manuscript Type: Research Paper Date Submitted by the Author: -Dec-0 Complete List of Authors:
More informationRICE UNIVERSITY. An Alternative Approach to Differential Semblance Velocity Analysis via Normal Moveout Correction. Chao Wang.
RICE UNIVERSITY An Alternative Approach to Differential Semblance Velocity Analysis via Normal Moveout Correction by Chao Wang A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE
More informationWaveform inversion for attenuation estimation in anisotropic media Tong Bai & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines
Waveform inversion for attenuation estimation in anisotropic media Tong Bai & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines SUMMARY Robust estimation of attenuation coefficients remains
More informationScattering-based decomposition of sensitivity kernels for full-waveform inversion
Scattering-based decomposition of sensitivity kernels for full-waveform inversion Daniel Leal Macedo (Universidade Estadual de Campinas) and Ivan asconcelos (Schlumberger Cambridge Research) Copyright
More informationSeismic imaging and multiple removal via model order reduction
Seismic imaging and multiple removal via model order reduction Alexander V. Mamonov 1, Liliana Borcea 2, Vladimir Druskin 3, and Mikhail Zaslavsky 3 1 University of Houston, 2 University of Michigan Ann
More informationON THE CONSTRUCTION OF DUALLY FLAT FINSLER METRICS
Huang, L., Liu, H. and Mo, X. Osaka J. Math. 52 (2015), 377 391 ON THE CONSTRUCTION OF DUALLY FLAT FINSLER METRICS LIBING HUANG, HUAIFU LIU and XIAOHUAN MO (Received April 15, 2013, revised November 14,
More informationApplied Mathematics 205. Unit V: Eigenvalue Problems. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit V: Eigenvalue Problems Lecturer: Dr. David Knezevic Unit V: Eigenvalue Problems Chapter V.4: Krylov Subspace Methods 2 / 51 Krylov Subspace Methods In this chapter we give
More informationIterative linearized migration and inversion
Stanford Exploration Project, Report 123, October 31, 2005, pages 101 128 Iterative linearized migration and inversion Huazhong Wang 1 ABSTRACT The objective of seismic imaging is to obtain an image of
More informationThe Lanczos and conjugate gradient algorithms
The Lanczos and conjugate gradient algorithms Gérard MEURANT October, 2008 1 The Lanczos algorithm 2 The Lanczos algorithm in finite precision 3 The nonsymmetric Lanczos algorithm 4 The Golub Kahan bidiagonalization
More informationTowards full waveform inversion: A torturous path
FWI Towards full waveform inversion: A torturous path J. Helen Isaac and Gary F. Margrave ABSTRACT Full waveform inversion (FWI) can be viewed as an iterative cycle involving forward modelling, pre-stack
More informationEdinburgh Research Explorer
Edinburgh Research Explorer Nonlinear Scattering-based Imaging in Elastic Media Citation for published version: Ravasi, M & Curtis, A 213, 'Nonlinear Scattering-based Imaging in Elastic Media' Paper presented
More informationA timely and necessary antidote to indirect methods and so-called P-wave FWI
A timely and necessary antidote to indirect methods and so-called P-wave FWI ARTHUR B. WEGLEIN, University of Houston Editor s note: The following article brings to light a cautionary concern (and a set
More informationAN EIGENVALUE STUDY ON THE SUFFICIENT DESCENT PROPERTY OF A MODIFIED POLAK-RIBIÈRE-POLYAK CONJUGATE GRADIENT METHOD S.
Bull. Iranian Math. Soc. Vol. 40 (2014), No. 1, pp. 235 242 Online ISSN: 1735-8515 AN EIGENVALUE STUDY ON THE SUFFICIENT DESCENT PROPERTY OF A MODIFIED POLAK-RIBIÈRE-POLYAK CONJUGATE GRADIENT METHOD S.
More informationNonlinear acoustic imaging via reduced order model backprojection
Nonlinear acoustic imaging via reduced order model backprojection Alexander V. Mamonov 1, Vladimir Druskin 2, Andrew Thaler 3 and Mikhail Zaslavsky 2 1 University of Houston, 2 Schlumberger-Doll Research
More informationParallel Numerics, WT 2016/ Iterative Methods for Sparse Linear Systems of Equations. page 1 of 1
Parallel Numerics, WT 2016/2017 5 Iterative Methods for Sparse Linear Systems of Equations page 1 of 1 Contents 1 Introduction 1.1 Computer Science Aspects 1.2 Numerical Problems 1.3 Graphs 1.4 Loop Manipulations
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationThe signature of attenuation and anisotropy on AVO and inversion sensitivities
The signature of attenuation and anisotropy on AVO and inversion sensitivities Shahpoor Moradi and Kristopher Innanen University of Calgary, Department of Geoscience, Calgary, Canada Summary We study the
More informationCompressive sampling meets seismic imaging
Compressive sampling meets seismic imaging Felix J. Herrmann fherrmann@eos.ubc.ca http://slim.eos.ubc.ca joint work with Tim Lin and Yogi Erlangga Seismic Laboratory for Imaging & Modeling Department of
More informationPreface to the Second Edition. Preface to the First Edition
n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More informationSome consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston
ain enu Some consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston Summary When S-wave velocity is absent, approximate Gassmann
More informationYouzuo Lin and Lianjie Huang
PROCEEDINGS, Thirty-Ninth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 24-26, 2014 SGP-TR-202 Building Subsurface Velocity Models with Sharp Interfaces
More information