Introduction to seismic body waves Tomography
|
|
- Milo Morgan
- 5 years ago
- Views:
Transcription
1 Introduction to seismic body waves Tomography Presented by Han Zhang Prepared by group work of Justin Wilgus and Han Zhang
2 Outlines 1. Velocity derived from wave equation 2. Two commonly used reference 1- D models 3. Naming seismic body waves at global scale 4. Basic ideas of body wave tomography 5. Possible improvement on inversion
3 Seismic wave equation we got in last class, ρu = λ + 2μ u μ u where dot above vector u denotes partial differential on time (. ), and two dots./ means do it twice (.0 ). And =.,.,. as we discussed. λ, μ and ρ are two./ lame parameters and density of media respectively. By dot timing operator from left side, we can get (using v = 0 ), ρ u = λ + 2μ 7 u which is compression wave equation since u denotes volume change of media. V 9 = λ + 2μ ρ Similarly, by taking the curl of two side, we can get, ρ u = μ 7 u for S wave since u is shear component of displacement. V ; = μ ρ
4 Seismic wave equation we got in last class, Despite the fact that density of materials inside Earth is ρu = λ + 2μ u μ u increasing with depth, lame parameters λ, μ are also increasing with depth and their enlargement ), factor and two is dots greater than ). that And of = density,.,.,. which as we discussed. finally leads λ, μ and to ρan are two./ increase function of wave velocity with depth inside mantle. where dot above vector u denotes partial differential on time (../ means do it twice (.0 lame parameters and density of media respectively. By dot timing operator from left side, we can get (using v = 0 ), ρ u = λ + 2μ 7 u which is compression wave equation since u denotes volume change of media. V 9 = λ + 2μ ρ Similarly, by taking the curl of two side, we can get, ρ u = μ 7 u for S wave since u is shear component of displacement. V ; = μ ρ
5 Reference Models PREM reference model km 400 km km Vs ρ Vp Depth (km) 2400 CMB ICB Velocity (km/s) or Density (g/cm ) 14
6 PREM reference model AK135 reference model km 400 km km km ρ Vs Vp km ρ Vs Vp Depth (km) CMB Depth (km) CMB ICB 5200 ICB Velocity (km/s) or Density (g/cm 3 ) Velocity (km/s) or Density (g/cm 3 )
7 PREM reference model AK135 PREM reference model km 400 km km 400 km km ρ Vs Vp km ρ Vs Vp Depth (km) CMB Depth (km) CMB ICB 5200 ICB Velocity (km/s) or Density (g/cm 3 ) Velocity (km/s) or Density (g/cm 3 )
8 Datasets of PREM and AK135 PREM 1. Geodetic data 2. Free oscillation and long- period surface wave 3. Body wave (primarily, direct P and S) Starting model of density derived from state equation. AK Body wave traveltime data (absolute traveltime of 18 seismic body waves and traveltime differences between them) Starting model can be regarded as IASP91
9 Phases Naming (global) source Generally, upgoing waves start with lower case letter while downgoing phases using upper case. Same rule for reflection/conversion points. Considering difference of sp and SP to see how to name different phases. (Shearer, Intro. to Seismology)
10 Phases Naming (exercise) Blue for P wave Red for S wave 60 Ray 1 Ray 2 ppcp 410 s 60 ScS 660 ScS CMB ICB 0 Animation available here
11 Body Wave Tomography Common assumptions of body wave traveltime tomography, 1. One good reference model available to make sure velocity perturbation of model region is small enough for linear inversion; 2. All blocks in our model region are well sampled by observations to reduce uncertainty of inversion result. Idea of body wave traveltime tomography based on ray theory, 1. Traveltime misfit between observed data and prediction of one certain ray is contributed by the blocks lie on this ray path only; 2. Total traveltime perturbation along the ray path can be summed from the product of traveltime in each block (calculated using reference model) with the fractional velocity perturbation within the block.
12 Body Wave Tomography r H = r 7 = r P = r Q = 0.01 sec; r S = r T = r U = r V = r W = 0 sec; a b c d e f 4 5 For each ray path, we can get one equation looks like this, r = = BCDEFG t =@ where r = is traveltime misfit between observation and prediction of ray path i, t =@ is absolute traveltime in block j of ray path i, and is velocity perturbation of block j. 8 g h i 6 7 For ray path 1 in the setting showed left, the equation should be, r H = P I + P J + P K by assuming t HI = t HJ = t HK = 1 sec. Rays are numbered aside their arrows and blocks are marked by low- case characters inside them. Simple traveltime (1 sec or 2 sec for one block) is assumed for practice.
13 Body Wave Tomography a b c d e f 4 5 r H = r 7 = r P = r Q = 0.01 sec; r S = r T = r U = r V = r W = 0 sec; 8 g h i 6 7 P I = P Y = 1%; P B = P E = P J = P [ = P K = P \ = P = = 0. White blocks have zero perturbation relative to reference model while black ones have perturbation of - 1%.
14 Body Wave Tomography Rewrite residual equation using matrix to get general form, d = Gm Residuals = Rays Model where d (1 x i) represents residuals of all the ray path, G (j x i) is related to ray path configuration and sensitive kernels adopted and m (1 x j) denotes perturbation model. In the case of i < j, which means the number of independent equations is less than parameters need to be recovered, we can get general solution only (infinitely many solutions). In the case of i = j, which means we have same amount of independent equations and model parameters, we will get exactly one solution. In the case of i > j, which means we have more independent equations than model parameters, theoretically we cannot get solution from above form. However, we are able to get least- square solution (L2 norm) of this kind of problem by solving G T d = G T Gm which theoretically have exactly one solution.
15 GAP_P4 Model (Fukao and Obayashi, 2013) Slab atop at 660 km discontinuity
16 GAP_P4 Model (Fukao and Obayashi, 2013) Slab cross 660 km discontinuity
17 Possible Improvements Body wave tomography works pretty well in helping us understand deep Earth structures, however, there are some defects of it, which may be overcame by joint inversion with other dataset. Some are listed here for checking, 1. L2 norm tend to results in a smooth model, which will average small- scale heterogeneity to its surroundings. Combining with converted wave modelling (such as receiver function). 2. Limited by the fact that most seismometers are deployed on continents and islands, which leads to poor ray path cover at crust and uppermost mantle beneath ocean. Combing with data from OBSs (ocean bottom seismometers) and/or electrical resistivity tomography (which is also expensive in ocean but is able to get very good shallow structures).
18
19 Finite Frequency Kernels Banana- doughnut kernels showing the sensitivity of P- wave travel times at 60 epicentral distance to velocity perturbations in the mantle. (Hung et al, 2000)
A NUMERICAL EXPERIMENT ON NONLINEAR IMAGE RECONSTRUCTION FROM FIRST-ARRIVAL TIMES FOR TWO-DIMENSIONAL ISLAND ARC STRUCTURE
Letter J. Phys. Earth, 34, 195-201, 1986 A NUMERICAL EXPERIMENT ON NONLINEAR IMAGE RECONSTRUCTION FROM FIRST-ARRIVAL TIMES FOR TWO-DIMENSIONAL ISLAND ARC STRUCTURE Ichiro NAKANISHI and Kaoru YAMAGUCHI
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 informationTRAVELTIME TOMOGRAPHY (CONT)
30 March 005 MODE UNIQUENESS The forward model TRAVETIME TOMOGRAPHY (CONT di = A ik m k d = Am (1 Data vecto r Sensitivit y matrix Model vector states the linearized relationship between data and model
More informationGEOPHYS 242: Near Surface Geophysical Imaging. Class 8: Joint Geophysical Inversions Wed, April 20, 2011
GEOPHYS 4: Near Surface Geophysical Imaging Class 8: Joint Geophysical Inversions Wed, April, 11 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical,
More informationAn imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults
Earth Planets Space, 52, 3 11, 2000 An imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults Takumi Murakoshi 1, Hiroshi Takenaka
More informationSeismic modelling with the reflectivity method
Yongwang Ma, Luiz Loures, and Gary F. Margrave ABSTRACT Seismic modelling with the reflectivity method Numerical seismic modelling is a powerful tool in seismic imaging, interpretation and inversion. Wave
More informationSuccessful 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 informationA Scalable Parallel LSQR Algorithm for Solving Large-Scale Linear System for Seismic Tomography
1 A Scalable Parallel LSQR Algorithm for Solving Large-Scale Linear System for Seismic Tomography He Huang, Liqiang Wang, Po Chen(University of Wyoming) John Dennis (NCAR) 2 LSQR in Seismic Tomography
More information3D S wave statics from direct VSP arrivals Michael J. O'Brien, Allied Geophysics, and Paritosh Singh, Colorado School of Mines
3D S wave statics from direct VSP arrivals Michael J. O'Brien, Allied Geophysics, and Paritosh Singh, Colorado School of Mines Summary We observe that direct S wave arrivals in 3D VSP data contain the
More informationWe N Converted-phase Seismic Imaging - Amplitudebalancing Source-independent Imaging Conditions
We N106 02 Converted-phase Seismic Imaging - Amplitudebalancing -independent Imaging Conditions A.H. Shabelansky* (Massachusetts Institute of Technology), A.E. Malcolm (Memorial University of Newfoundland)
More information(x, y, z) m 2. (x, y, z) ...] T. m 2. m = [m 1. m 3. Φ = r T V 1 r + λ 1. m T Wm. m T L T Lm + λ 2. m T Hm + λ 3. t(x, y, z) = m 1
Class 1: Joint Geophysical Inversions Wed, December 1, 29 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical, statistical Deal with data in the same
More informationSparse wavelet expansions for seismic tomography: Methods and algorithms
Sparse wavelet expansions for seismic tomography: Methods and algorithms Ignace Loris Université Libre de Bruxelles International symposium on geophysical imaging with localized waves 24 28 July 2011 (Joint
More informationGEOPHYS 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 informationTravel Time Tomography using Neural Networks
Travel Time Tomography using Neural Networks Yoshiya Oda Tokyo Metropolitan University, Japan Tomohisa Ishiyama Tokyo Metropolitan University, Japan Shinya Yokono Tokyo Metropolitan University, Japan SUMMARY:
More informationVelocity (Vs) Display (1-D Profile)
PS User Guide Series 2015 Velocity (Vs) Display (1-D Profile) Prepared By Choon B. Park, Ph.D. January 2015 Table of Contents Page 1. Overview 2 2. Files 6 3. View 7 4. Edit 9 5. Chart 11 6. Layer 12 7.
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 informationP063 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 informationChapter 2 Methodology of Seismic Tomography
Chapter Methodology of Seismic Tomography Abstract In this chapter we introduce the tomographic methods which are widely used to study three-dimensional (3-D) seismic velocity attenuation and anisotropy
More informationSUMMARY. 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 informationSeismic velocity tomography with co-located soft data
Seismic velocity tomography with co-located soft data Bob Clapp Stanford University SCRF Affiliate Meeting May 2009 Motivation soft data constrains flow simulation model 2 Outline Introduction Seismic
More informationVidale, J., 1988, Finite-difference calculation of travel times, BSSA, Vo. 78, No. 6,
Class 4: First-Arrival Traveltime Tomography Mon, Sept 21, 2009 Wavefront tracing methods Inversion algorithms Model regularization Review of tomography case histories, presentation by Dr. Xianhuai Zhu
More informationUltrasonic Multi-Skip Tomography for Pipe Inspection
18 th World Conference on Non destructive Testing, 16-2 April 212, Durban, South Africa Ultrasonic Multi-Skip Tomography for Pipe Inspection Arno VOLKER 1, Rik VOS 1 Alan HUNTER 1 1 TNO, Stieltjesweg 1,
More informationReflectivity modeling for stratified anelastic media
Reflectivity modeling for stratified anelastic media Peng Cheng and Gary F. Margrave Reflectivity modeling ABSTRACT The reflectivity method is widely used for the computation of synthetic seismograms for
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 informationTotal variation tomographic inversion via the Alternating Direction Method of Multipliers
Total variation tomographic inversion via the Alternating Direction Method of Multipliers Landon Safron and Mauricio D. Sacchi Department of Physics, University of Alberta Summary Geophysical inverse problems
More informationEfficient 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 informationThe geometry of reflection and refraction Wave conversion and reflection coefficient
7.2.2 Reflection and Refraction The geometry of reflection and refraction Wave conversion and reflection coefficient The geometry of reflection and refraction A wave incident on a boundary separating two
More informationComparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media
Raytracing in VTI media Comparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media Jianli Yang and Don C. Lawton ABSTRACT Determination of the conversion
More informationComparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media
Comparison of two P-S conversion-point mapping approaches for Vertical Transversely Isotropic (VTI) media Jianli Yang* and Don C. Lawton CREWES, University of Calgary, 5 UniversityDR.NW, Calgary, AB, TN
More informationInterval velocity estimation through convex optimization
Stanford Exploration Project, Report 125, January 16, 2007, pages 1?? Interval velocity estimation through convex optimization Ben Witten and Michael Grant ABSTRACT Convex optimization is an optimization
More informationApplications of the tau-p Inversion Method. Shoba Maraj
Applications of the tau-p Inversion Method Shoba Maraj Final Report Theoretical Seismology Geop 523 Due 6th May 2011 Introduction Generally, velocity increases with depth in the Earth. They are related
More informationSUMMARY THEORY INTRODUCTION
Double-Difference Waveform Inversion of 4D Ocean Bottom Cable Data: Application to Valhall, North Sea Di Yang, Michael Fehler and Alison Malcolm, MIT, Faqi Liu and Scott Morton, Hess Corporation SUMMARY
More informationLecture 2: Seismic Inversion and Imaging
Lecture 2: Seismic Inversion and Imaging Youzuo Lin 1 Joint work with: Lianjie Huang 1 Monica Maceira 2 Ellen M. Syracuse 1 Carene Larmat 1 1: Los Alamos National Laboratory 2: Oak Ridge National Laboratory
More informationFinite-difference elastic modelling below a structured free surface
FD modelling below a structured free surface Finite-difference elastic modelling below a structured free surface Peter M. Manning ABSTRACT This paper shows experiments using a unique method of implementing
More informationSeparation of specular reflection and diffraction images in Kirchhoff depth migration Faruq E Akbar and Jun Ma, SEIMAX Technologies, LP
Separation of specular reflection and diffraction images in Kirchhoff depth migration Faruq E Akbar and Jun Ma, SEIMAX Technologies, LP Summary Seismic diffractions may occur from faults, fractures, rough
More informationPEAT SEISMOLOGY Lecture 16: Seismic Tomography I
PEAT8002 - SEISMOLOGY Lecture 16: Seismic Tomography I Nick Rawlinson Research School of Earth Sciences Australian National University Introduction Seismic data represent one of the most valuable resources
More informationSpectral-Element and Adjoint Methods in Seismology
Spectral-Element and Adjoint Methods in Seismology Jeroen Tromp Min Chen, Vala Hjorleifsdottir, Dimitri Komatitsch, Swami Krishnan, Qinya Liu, Alessia Maggi, Anne Sieminski, Carl Tape & Ying Zhou Introduction
More informationSummary. Introduction
Dmitry Alexandrov, Saint Petersburg State University; Andrey Bakulin, EXPEC Advanced Research Center, Saudi Aramco; Pierre Leger, Saudi Aramco; Boris Kashtan, Saint Petersburg State University Summary
More informationTh LHR5 08 Multi-modal Surface Wave Inversion and Application to North Sea OBN Data
Th LHR5 08 Multi-modal Surface Wave Inversion and pplication to North Sea ON Data S. Hou (CGG), D. Zheng (CGG), X.G. Miao* (CGG) & R.R. Haacke (CGG) SUMMRY Surface-wave inversion (SWI) for S-wave velocity
More informationChapter 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 informationSEG/San Antonio 2007 Annual Meeting
Yaofeng He* and Ru-Shan Wu Modeling and Imaging Laboratory, Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, CA, 95064, USA Summary A one-way and one-return boundary
More informationSeismic refraction surveys
Seismic refraction surveys Seismic refraction surveys generate seismic waves that are refracted back to Earth s surface from velocity and density discontinuities at depth. Uses : Small-scale: geotechnical
More information2010 SEG SEG Denver 2010 Annual Meeting
Localized anisotropic tomography with checkshot : Gulf of Mexico case study Andrey Bakulin*, Yangjun (Kevin) Liu, Olga Zdraveva, WesternGeco/Schlumberger Summary Borehole information must be used to build
More informationTravel-Time Sensitivity Kernels in Long-Range Propagation
Travel-Time Sensitivity Kernels in Long-Range Propagation Emmanuel Skarsoulis Foundation for Research and Technology Hellas Institute of Applied and Computational Mathematics P.O. Box 1385, GR-71110 Heraklion,
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 informationPS 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 informationElectromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah
Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Summary In this paper we present a new approach to the interpretation
More informationRecovering 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 informationGG450 4/5/2010. Today s material comes from p and in the text book. Please read and understand all of this material!
GG450 April 6, 2010 Seismic Reflection I Today s material comes from p. 32-33 and 81-116 in the text book. Please read and understand all of this material! Back to seismic waves Last week we talked about
More informationMultimode surface wave tomography for the Australian region using a three-stage approach incorporating finite frequency effects
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2002jb002254, 2004 Multimode surface wave tomography for the Australian region using a three-stage approach incorporating finite frequency effects
More informationDownloaded 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 informationE044 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 informationAnisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS
Anisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS Summary Anisotropic depth model building using surface seismic data alone is non-unique and
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 informationTh E An Interferometric Interpretation of Marchenko Redatuming
Th E0 An Interferometric Interpretation of Marchenko Redatuming J. van der Neut (Delft University of Technology), I. Vasconcelos (chlumberger Gould Research) & K. Wapenaar (Delft University of Technology)
More informationG021 Subsalt Velocity Analysis Using One-Way Wave Equation Based Poststack Modeling
G021 Subsalt Velocity Analysis Using One-Way Wave Equation Based Poststack Modeling B. Wang* (CGG Americas Inc.), F. Qin (CGG Americas Inc.), F. Audebert (CGG Americas Inc.) & V. Dirks (CGG Americas Inc.)
More informationIntroduction to ANSYS DesignXplorer
Lecture 4 14. 5 Release Introduction to ANSYS DesignXplorer 1 2013 ANSYS, Inc. September 27, 2013 s are functions of different nature where the output parameters are described in terms of the input parameters
More informationPressure and shear wave separation of ocean bottom seismic data
Pressure and shear wave separation of ocean bottom seismic data Ohad Barak ABSTRACT PZ summation is a common method for separating the upgoing wavefield from the downgoing wavefield in data acquired by
More informationFull 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 informationWS10-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 informationA MODELING STUDY OF LOW-FREQUENCY CSEM IN SHALLOW WATER
A MODELING STUDY OF LOW-FREQUENCY CSEM IN SHALLOW WATER Michael A. Frenkel and Sofia Davydycheva EMGS ASA ABSTRACT The applicability and resolution power of low-frequency Controlled-Source Electromagnetic
More informationImage-guided 3D interpolation of borehole data Dave Hale, Center for Wave Phenomena, Colorado School of Mines
Image-guided 3D interpolation of borehole data Dave Hale, Center for Wave Phenomena, Colorado School of Mines SUMMARY A blended neighbor method for image-guided interpolation enables resampling of borehole
More informationMachine-learning Based Automated Fault Detection in Seismic Traces
Machine-learning Based Automated Fault Detection in Seismic Traces Chiyuan Zhang and Charlie Frogner (MIT), Mauricio Araya-Polo and Detlef Hohl (Shell International E & P Inc.) June 9, 24 Introduction
More informationH003 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 informationThe Development of a Volume-of-Fluid Interface Tracking Method for Modeling Problems in Mantle Convection
The Development of a Volume-of-Fluid Interface Tracking Method for Modeling Problems in Mantle Convection Jonathan Robey 1 September 12, 2016 1 Funded in part by a GAANN Fellowship 1 Motivation Mantle
More informationInteractive Receiver Function Forward Modeller-2 (IRFFM2 v1.2) General Information & User Manual
Interactive Receiver Function Forward Modeller-2 (IRFFM2 v1.2) Hrvoje Tkalčić 1, Debdeep Banerjee, Jason Li and Chris Tarlowski 2 Research School of Earth Sciences The Australian National University Canberra,
More informationNoise-thresholding sparse-spike inversion with global convergence: calibration and applications
Noise-thresholding sparse-spike inversion with global convergence: calibration and applications Douglas S. Sassen and Michael Glinsky, Ion Geophysical 2013 SEG Douglas Sassen 9/25/13 1 Illustrative Abstract!
More informationHigh definition tomography brings velocities to light Summary Introduction Figure 1:
Saverio Sioni, Patrice Guillaume*, Gilles Lambaré, Anthony Prescott, Xiaoming Zhang, Gregory Culianez, and Jean- Philippe Montel (CGGVeritas) Summary Velocity model building remains a crucial step in seismic
More informationSUMMARY. computationally more expensive but suitable for currently available parallel computers.
algorithms on parallel computers Suhas Phadke*, Dheeraj Bhardwaj and Sudhakar Yerneni Centre for Development of Advanced Computing, Pune University Campus, GaneshKhind, Pune 411007, India SUMMARY This
More informationResearch Collection. Localisation of Acoustic Emission in Reinforced Concrete using Heterogeneous Velocity Models. Conference Paper.
Research Collection Conference Paper Localisation of Acoustic Emission in Reinforced Concrete using Heterogeneous Velocity Models Author(s): Gollob, Stephan; Vogel, Thomas Publication Date: 2014 Permanent
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 informationWe 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 informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Aliasing Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 16, 2014 Bee Bednar (Panorama Technologies) Seismic
More informationCAPTURING UNCERTAINTIES IN SOURCE-SPECIFIC STATION CORRECTIONS DERIVED FROM A 3-D MODEL. Sponsored by Defense Threat Reduction Agency
CAPTURING UNCERTAINTIES IN SOURCE-SPECIFIC STATION CORRECTIONS DERIVED FROM A 3-D MODEL N.M. Shapiro, 1 M.H. Ritzwoller, 1 M.P. Barmin, 1 P.B. Weichman, 2 L.F. Tenorio, 3 W.C. Navidi 3 University of Colorado
More informationTu Crossline Reconstruction Using Aliased 3D Deghosted Up- and Downgoing Wavefields
Tu-04-10 Crossline Reconstruction Using Aliased 3D Deghosted Up- and Downgoing Wavefields D.-J. van Manen* (Schlumberger), M. Vassallo (WesternGeco), A.K. Özdemir (WesternGeco), A. Özbek (Schlumberger)
More informationShortest-path calculation of first arrival traveltimes by expanding wavefronts
Stanford Exploration Project, Report 82, May 11, 2001, pages 1 144 Shortest-path calculation of first arrival traveltimes by expanding wavefronts Hector Urdaneta and Biondo Biondi 1 ABSTRACT A new approach
More informationTying well information and seismic data
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 319 Tying well information and seismic data Arnaud Berlioux 1 ABSTRACT Well log information and seismic data for a given horizon may not tie
More informationThe. Reverse-time depth migration in elastic media. Introduction FOCUS ARTICLE. Zaiming Jiang. Coordinated by Rob Holt
FOCUS ARTICLE Coordinated by Rob Holt Reverse-time depth migration in elastic media Zaiming Jiang KEY SEISMIC SOLUTIONS, CALGARY, ALBERTA, CANADA The main objective of this paper is to sketch a reverse-time
More informationM. 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 informationRay tracing in geophysics
Ray tracing in geophysics Seminar Author: Klemen Kunstelj Adviser: dr. Andrej Gosar and dr. Jure Bajc 18.12.2002 Abstract Alp 2002 is an international seismological project. The main issue is to construct
More informationAn illustration of adaptive Marchenko imaging
An illustration of adaptive Marchenko imaging Joost van der Neut 1, Kees Wapenaar 1, Jan Thorbecke 1, Evert Slob 1, and Ivan Vasconcelos 2 Abstract In Marchenko imaging, wavefields are retrieved at specified
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 informationu = 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 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 informationTesting Joint Inversion Code with Geologically Realistic Models
Testing Joint Inversion Code with Geologically Realistic Models Angela Carter-McAuslan, Peter Lelièvre, Gary Blades, Colin Farquharson GAC-MAC Joint Annual Meeting May 27 th, 2012 Background Joint Inversion
More informationRealtime Water Simulation on GPU. Nuttapong Chentanez NVIDIA Research
1 Realtime Water Simulation on GPU Nuttapong Chentanez NVIDIA Research 2 3 Overview Approaches to realtime water simulation Hybrid shallow water solver + particles Hybrid 3D tall cell water solver + particles
More informationJoint seismic traveltime and TEM inversion for near surface imaging Jide Nosakare Ogunbo*, Jie Zhang, GeoTomo LLC
Jide Nosaare Ogunbo*, Jie Zhang, GeoTomo LLC Summary For a reliable interpretation of the subsurface structure, the joint geophysical inversion approach is becoming a viable tool. Seismic and EM methods
More informationLOTOS Code for Local Earthquake Tomographic Inversion: Benchmarks for Testing Tomographic Algorithms
Bulletin of the Seismological Society of America, Vol. 99, No. 1, pp. 194 214, February 2009, doi: 10.1785/0120080013 LOTOS Code for Local Earthquake Tomographic Inversion: Benchmarks for Testing Tomographic
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 informationGEOPHYS 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 informationInterpolation by Spline Functions
Interpolation by Spline Functions Com S 477/577 Sep 0 007 High-degree polynomials tend to have large oscillations which are not the characteristics of the original data. To yield smooth interpolating curves
More informationHeadwave 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 informationVelocity Model Sampling and Interpolation: The Gradient Problem
January 11, 2013 1 2 3 4 5 Kirchhoff migration uses raytracing equations for prestack depth migration. d x dt = v( x)2 p d p dt = v( x) v( x) These equations estimate the paths of seismic energy through
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 informationDeconvolution in the radial trace domain
R-T domain deconvolution Deconvolution in the radial trace domain David C. Henley ABSTRACT The radial trace (R-T) domain has been shown to be useful for coherent noise attenuation and other seismic wavefield
More informationQ-compensation in complex media ray-based and wavefield extrapolation approaches Maud Cavalca, Robin Fletcher and Marko Riedel, WesternGeco.
ray-based and wavefield extrapolation approaches Maud Cavalca, Robin Fletcher and Marko Riedel, WesternGeco. Summary We apply and compare three model-based Q-compensation approaches. The first two approaches
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 informationREDUCED ORDER MODELING IN MULTISPECTRAL PHOTOACOUSTIC TOMOGRAPHY
REDUCED ORDER MODELING IN MULTISPECTRAL PHOTOACOUSTIC TOMOGRAPHY Arvind Saibaba Sarah Vallélian Statistical and Applied Mathematical Sciences Institute & North Carolina State University May 26, 2016 OUTLINE
More informationMulticomponent 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 informationChapter N/A: The Earthquake Cycle and Earthquake Size
This sheet provides a general outline of the topics covered on your upcoming exam. This is not an exhaustive list, but is simply a general list of key points. In other words, most exam questions will be
More information