Spectral-Element and Adjoint Methods in Seismology
|
|
- Bonnie Turner
- 6 years ago
- Views:
Transcription
1 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
2 Introduction to the Spectral-Element Method (SEM)
3 Governing Equations Equation of motion: Boundary condition: Initial conditions: Earthquake source:
4 Weak Form Weak form valid for any test vector Boundary conditions automatically included Source term explicitly integrated Finite-fault (kinematic) rupture:
5 Finite-Elements Mapping from reference cube to hexahedral elements: Volume relationship: Jacobian of the mapping: Jacobian matrix:
6 Lagrange Polynomials and Gauss- Lobatto-Legendre (GLL) Points The 5 degree 4 Lagrange polynomials: Degree 4 GLL points: General definition: GLL points are n+1 roots of: Note that at a GLL point:
7 Interpolation Representation of functions on an element in terms of Lagrange polynomials: Gradient on an element:
8 Integration Integration of functions over an element based upon GLL quadrature: Integrations are pulled back to the reference cube In the SEM one uses: interpolation on GLL points GLL quadrature Degree 4 GLL points:
9 The Diagonal Mass Matrix Representation of the displacement: Degree 4 Lagrange polynomials: Representation of the test vector: Weak form: Diagonal mass matrix: Integrations are pulled back to the reference cube In the SEM one uses: interpolation on GLL points GLL quadrature Degree 4 GLL points:
10 Assembly Need to distinguish: local mesh global mesh Efficient routines available from FEM applications Global equations: Global SEM time-marching is based upon an explicit second-order time scheme
11 Parallel Implementation Global mesh partitioning: Regional mesh partitioning: Cubed Sphere: 6 n2 mesh slices n x m mesh slices
12 Basin Code: SPECFEM3D_BASIN Freely available for non-commercial use from geodynamics.org Manual available at geodynamics.org 3D Attenuation 3D Anisotropy Ocean load Topography & bathymetry Kinematic ruptures Movies Models: Harvard 3D Southern California model SOCAL 1D model Homogeneous half-space model Adjoint capabilities
13 Southern California Simulations A A B B Vp(km/s) n x m mesh slices
14 Periods > 2 s Komatitsch et al. 2004
15 June 12, 2005, M=5.1 Big Bear QuickTime and a YUV420 codec decompressor are needed to see this picture.
16 3D Regional Forward Simulations Periods > 6 s 379 km June 12, 2005, M=5.1 Big Bear 100km Qinya Liu
17 San Andreas Rupture Scenario QuickTime and a YUV420 codec decompressor are needed to see this picture.
18 San Andreas Rupture Scenario: Quantitative Seismic Hazard Assessment Swami Krishnan
19 Near Real-Time Applications Automated near real-time simulations of all M>3.5 events ShakeMovies at Soon: CMT solutions Synthetic seismograms
20 SPECFEM3D_BASIN: Future Plans Switch to a (parallel) CUBIT hexahedral finite-element mesher (Casarotti, Lee) Topography & bathymetry Major geological interfaces Basins Fault surfaces Use ParMETIS or SCOTCH for mesh partitioning & load-balancing Retain the SPECFEM3D_BASIN solver (takes ParMETIS meshes; Komatitsch) Add dynamic rupture capabilities (Ampuero, Lapusta, Kaneko) Taipei Basin 102km N 88km 100km S.-J. Lee
21 Global Code: SPECFEM3D_GLOBE Freely available for non-commercial use from geodynamics.org Manual available at geodynamics.org 3D Attenuation 3D Anisotropy Ocean load Topography & bathymetry Rotation Self-gravitation (Cowling approximation) Kinematic ruptures Movies Models: 1D models: isotropic PREM, transversely isotropic PREM, AK135, IASP91, 1066A S20RTS Crust2.0 Adjoint capabilities
22 Global Simulations PREM benchmarks Cubed sphere mesh Parallel implementation: 2 6 n mesh slices
23 SEM Implementation of Attenuation Anelastic, anisotropic constitutive relationship: Equivalent Standard Linear Solid (SLS) formulation: Attenuation (3 SLSs) Memory variable equation: Unrelaxed modulus: Modulus defect: Physical dispersion (3 SLSs)
24 Effect of Attenuation
25 Attenuation
26 Effect of Anisotropy
27 SEM Implementation of Anisotropy
28 Antipodal Transverse Record
29 Vertical PKP
30 Rotation & Self-Gravitation Wave equation solved by SPECFEM3D_GLOBE in crust, mantle and inner core: Wave equation solved by SPECFEM3D_GLOBE in fluid outer core: SPECFEM3D_GLOBE uses domain decomposition between the fluid outer core and the solid inner core and mantle matching exactly: continuity of traction continuity of the normal component of displacement
31 Full Gravity Versus Cowling
32 SEM Implementation of Gravity
33 Effect of Ocean
34 SEM Implementation of Ocean Modified boundary condition:
35 3D Mantle Models S20RTS (Ritsema et al. 1999)
36 Crustal and Topographic Models Crust 2.0 (Bassin et al. 2000) ETOPO5
37 Great 2004 Sumatra-Andaman Earthquake Finite slip model (Chen et al., 2005) Main shock & aftershocks (Harvard)
38 Sumatra Surface Waves QuickTime and a YUV420 codec decompressor are needed to see this picture.
39 Surface-Wave Fits Vala Hjorleifsdottir
40 SPECFEM3D_GLOBE: Future Plans On-demand TeraGrid applications: Automated, near real-time simulations of all M>6 earthquakes Analysis of past events (more than 20,000 events) Seismology Web Portal (prototype available at this meeting) Petascale simulations: Global simulations at 1-2 Hz New doubling brick (perfect load-balancing)
41 Adjoint Spectral-Element Simulations (ASEM)
42 Adjoint Tomography PDE-constrained waveform tomography: Change in the waveform misfit function:
43 Adjoint Equations Adjoint wavefield: Adjoint equation of motion: Adjoint boundary conditions: Adjoint initial conditions: Adjoint source:
44 Frechet derivative The Frechet derivative may be expressed as: Density and elastic tensor kernels:
45 Isotropic Kernels For isotropic perturbations we have: where and we have defined the strain deviators In terms of wave speeds: where
46 Traveltime Frechet Derivatives Traveltime tomography: Change in the misfit function: Traveltime anomaly in terms of banana-donut kernel (Dahlen): The kernel is a weighted sum of banana-donut kernels :
47 Adjoint Wavefield Kernel in terms of the adjoint wavefield: The adjoint wavefield is generated by the adjoint source Notes: Need simultaneous access to the regular wavefield at time t and the adjoint wavefield at time T t Use of the time-reversed velocity as the source for the adjoint wavefield
48 Numerical Implementation Need simultaneous access to and During calculation of adjoint field, reconstruct by solving the `backward wave equation Need to store from a previous forward simulation: Last snapshot Wavefield absorb on artificial boundaries Challenge: `Undoing attenuation
49 2D Adjoint Tomography
50 Construction of a Banana-Donut Kernel Tape et al. 2006
51 Adjoint Tomography
52 Construction of an Event Kernel Event Kernel: Sum of weighted bananadonut kernels Two simulations per event
53 Construction of the Misfit Kernel Misfit kernel: Sum of all the event kernels Two simulations per event Gradient of misfit function
54 Conjugate Gradient Algorithm
55 Conjugate Gradient Algorithm
56 Joint Structure-Source Inversion
57
58 Southern California Rayleigh Waves 25 events
59 Toward 3D Tomography: SPECFEM3D Adjoint Capabilities
60 3D Sensitivity Kernels S SS PS PS SPECFEM3D_BASIN Liu & Tromp 2006
61 3D Body-Wave Sensitivity Kernels September 3, 2002, M=4.2 Yorba Linda
62 3D Surface-Wave Sensitivity Kernels Rayleigh (HEC) Love (HEC)
63 Our First 3D Event Kernels!
64 Global 3D Body-Wave Sensitivity Kernels Finite-Frequency effects 20 second P wave 9 second P wave SPECFEM3D_GLOBE Qinya Liu
65 Global 3D Body-Wave Sensitivity Kernels Pdiff PKP ScS SKS Qinya Liu
66 PKP Kernel
67 P P Kernel Liu & Tromp 2006
68 Global 3D Transversely Isotropic Kernels α h P α v α h SH α v β h β v β h β v η ρ η ρ Anne Sieminski
69 Transversely Isotropic Parameters: A, C, L, N, F 1 ζ: J, K, M 2 ζ: G, B, H 3 ζ: D 4 ζ: E Anne Sieminski
70 1D Versus 3D Kernels 1D 3D Ying Zhou
71 Time-Reversal Imaging: Glacial Earthquakes Greenland Carene Larmat
72 Measuring all available Phase & Amplitude Anomalies
73
74 Measure all suitable phases RAR F2=0.93 dt=1.00 F2=0.83 dt=1.00 F2=0.83 dt=6.00 F2=0.68 dt=2.00 F2=0.78 dt=-1.00 PP PKiKP ppkikp spkikp SKiKP SKSac S SKKSac psksac ps ssksac ss P PcP pp ScS SPn sp PnS SS SKKSdf SKKSac P P df Time (s) Seismograms ANMO F2=0.91 F2=0.93 F2=0.68 F2=0.92 F2=0.91 F2=0.83 F2=0.84 F2=0.83 F2=0.69 F2=0.85 F2=0.91 F2=0.92 F2=0.51 F2=0.61 F2= Time (s) Alessia Maggi
75 Conclusions Adjoint methods: Choose an observable, e.g., waveforms or cross-correlation traveltimes Choose a measure of misfit, e.g., least-squares Determine the appropriate adjoint source for this observable & measurement Use fully 3D reference models Any arrival suitable for measurement No dependence on the number of stations, components, or measurements 3D sensitivity kernels may be calculated based upon two forward simulations for each earthquake Number of simulations: 3 * (# earthquakes) * (# iterations) Full anisotropy for the same cost Attenuation remains a challenge Regional simulations: One 3 minute forward simulation accurate to 1.5 seconds takes 45 minutes on a 75 node cluster 150 events and 3 iterations would require 1800 simulations, i.e., three weeks of dedicated CPU time on 75 nodes Near real-time simulations Global simulations: One 1 hour forward simulation accurate to 20 seconds takes 4 hours on a 75 node cluster 500 events and 3 iterations would require 6,000 simulations, i.e., 100 days on a 750 node cluster Near real-time simulations On-demand global seismology Petascale application
CO2 sequestration crosswell monitoring based upon spectral-element and adjoint methods
CO2 sequestration crosswell monitoring based upon spectral-element and adjoint methods Christina Morency Department of Geosciences, Princeton University Collaborators: Jeroen Tromp & Yang Luo Computational
More informationSpectral-element Simulations of Elastic Wave Propagation in Exploration and Geotechnical Applications
Spectral-element Simulations of Elastic Wave Propagation in Exploration and Geotechnical Applications Lin Zheng*, Qi Zhao*, Qinya Liu, Bernd Milkereit, Giovanni Grasselli, University of Toronto, ON, Canada
More informationSPECFEM3D_GLOBE. Mostly developed at Caltech (USA) and University of Pau (France) History: v1.0: 1999/2000 ; v3.6: 2005; v4.
SPECFEM3D_GLOBE Min Chen Vala Hjörleifsdóttir Sue Kientz Dimitri Komatitsch Qinya Liu Alessia Maggi David Michéa Brian Savage Bernhard Schuberth Leif Strand Carl Tape Jeroen Tromp The SPECFEM3D source
More informationHigh-Order Finite-Element Earthquake Modeling on very Large Clusters of CPUs or GPUs
High-Order Finite-Element Earthquake Modeling on very Large Clusters of CPUs or GPUs Gordon Erlebacher Department of Scientific Computing Sept. 28, 2012 with Dimitri Komatitsch (Pau,France) David Michea
More informationA Framework for Online Inversion-Based 3D Site Characterization
A Framework for Online Inversion-Based 3D Site Characterization Volkan Akcelik, Jacobo Bielak, Ioannis Epanomeritakis,, Omar Ghattas Carnegie Mellon University George Biros University of Pennsylvania Loukas
More informationA Wavelet-Based Method for Simulation of Seismic Wave Propagation
AGU Fall Meeting, A Wavelet-Based Method for Simulation of Seismic Wave Propagation Tae-Kyung Hong & B.L.N. Kennett Research School of Earth Sciences The Australian National University Abstract Seismic
More informationIntroduction to seismic body waves Tomography
Introduction to seismic body waves Tomography Presented by Han Zhang Prepared by group work of Justin Wilgus and Han Zhang Outlines 1. Velocity derived from wave equation 2. Two commonly used reference
More informationToward Exascale Seismic Imaging & Inversion
Toward Exascale Seismic Imaging & Inversion Jeroen Tromp Department of Geosciences Program in Applied & Computational Mathematics Princeton Institute for Computational Science & Engineering Ebru Bozdağ,
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 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 informationSpectral-element simulations of elastic wave propagation in exploration and geotechnical applications
Earthq Sci (2014) 27(2):179 187 DOI 10.1007/s11589-014-0069-9 RESEARCH PAPER Spectral-element simulations of elastic wave propagation in exploration and geotechnical applications Lin Zheng Qi Zhao Bernd
More information3D Partitioning in SPECFEM3D Internal Mesher
Available online at www.prace-ri.eu Partnership for Advanced Computing in Europe 3D Partitioning in SPECFEM3D Internal Mesher Valentin Pavlov * Rila Solutions EAD, Acad. G. Bonchev str., bl. 27, Sofia
More informationComparison of Numerical Methods for Seismic Wave Propagation and Source Dynamics - the SPICE Code Validation
Third International Symposium on the Effects of Surface Geology on Seismic Motion Grenoble, France, 30 August - 1 September 2006 Paper Number: 62 Comparison of Numerical Methods for Seismic Wave Propagation
More informationComputational Fluid Dynamics - Incompressible Flows
Computational Fluid Dynamics - Incompressible Flows March 25, 2008 Incompressible Flows Basis Functions Discrete Equations CFD - Incompressible Flows CFD is a Huge field Numerical Techniques for solving
More informationOverview of the High-Order ADER-DG Method for Numerical Seismology
CIG/SPICE/IRIS/USAF WORKSHOP JACKSON, NH October 8-11, 2007 Overview of the High-Order ADER-DG Method for Numerical Seismology 1, Michael Dumbser2, Josep de la Puente1, Verena Hermann1, Cristobal Castro1
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 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 information(f) (a) Time (s)
A spectral element method or wave propagation simulation near a uid-solid interace Dimitri Komatitsch (1), Christophe Barnes (2), Jeroen Tromp (1) (1) Harvard University, Mailcode 252-21, Caltech, Pasadena,
More informationComputational Seismology: Simulating Seismic Wavefields for AlpArray
Computational Seismology: Simulating Seismic Wavefields for AlpArray Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich 1 Introduction 2 Goals of lecture Introduction
More information3D PHYSICS-BASED NUMERICAL SIMULATIONS: ADVANTAGES AND LIMITATIONS OF A NEW FRONTIER TO EARTHQUAKE GROUND MOTION PREDICTION.
6th National Conference on Earthquake Engineering & 2nd National Conference on Earthquake Engineering and Seismology June 14-16, 2017 3D PHYSICS-BASED NUMERICAL SIMULATIONS: ADVANTAGES AND LIMITATIONS
More informationComputation of the gravity gradient tensor due to topographic masses using tesseroids
Computation of the gravity gradient tensor due to topographic masses using tesseroids Leonardo Uieda 1 Naomi Ussami 2 Carla F Braitenberg 3 1. Observatorio Nacional, Rio de Janeiro, Brazil 2. Universidade
More informationAn Efficient 3D Elastic Full Waveform Inversion of Time-Lapse seismic data using Grid Injection Method
10 th Biennial International Conference & Exposition P 032 Summary An Efficient 3D Elastic Full Waveform Inversion of Time-Lapse seismic data using Grid Injection Method Dmitry Borisov (IPG Paris) and
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 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 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 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 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 informationImprovements 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 informationDownloaded 10/23/13 to Redistribution subject to SEG license or copyright; see Terms of Use at where P.
eismic wave scattering in D random media: A finite-difference simulation and slowness domain analysis Xiao-Bi Xie* Institute of Geophysics and lanetary hysics, University of California at anta Cruz ummary
More information2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123
2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs
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 informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) doi: 10.1111/j.1365-246X.2012.05676.x Numerical simulation of seismic wave propagation produced by earthquake by using a particle method Junichi
More informationTWO-DIMENSIONAL FEM MODELING OF SEISMIC WAVE PROPAGATION IN COAL SEAM. By Hongliang Wang
EGEE 50 Fall 005 Semester Paper. Instructor: Dr. Derek Elsworth TWO-DIMENSIONAL FEM MODELING OF SEISMIC WAVE PROPAGATION IN COAL SEAM By Hongliang Wang (huw103@psu.edu) Abstract In-seam seismic signals
More informationElastic 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 informationFinite Volume Discretization on Irregular Voronoi Grids
Finite Volume Discretization on Irregular Voronoi Grids C.Huettig 1, W. Moore 1 1 Hampton University / National Institute of Aerospace Folie 1 The earth and its terrestrial neighbors NASA Colin Rose, Dorling
More informationHigh-Frequency Simulations of Global Seismic Wave Propagation Using SPECFEM3D_GLOBE
High-Frequency Simulations of Global Seismic Wave Propagation Using SPECFEM3D_GLOBE Laura Carrington a, Dimitri Komatitsch b,c, Michael Laurenzano a, Mustafa M Tikir a, David Michéa b, Nicolas Le Goff
More information10 years of CyberShake: Where are we now and where are we going
10 years of CyberShake: Where are we now and where are we going with physics based PSHA? Scott Callaghan, Robert W. Graves, Kim B. Olsen, Yifeng Cui, Kevin R. Milner, Christine A. Goulet, Philip J. Maechling,
More informationA Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere
A Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere Ram Nair Computational and Information Systems Laboratory (CISL) National Center for Atmospheric Research
More informationis decomposed into P and S components localized in both space and slowness where P
eismic wave scattering in D random media: A finite-difference simulation and slowness domain analysis Xiao-Bi Xie* Institute of Geophysics and lanetary hysics, University of California at anta Cruz ummary
More informationEfficient data-driven strategy for 3D model-preconditioning FWI
Efficient data-driven strategy for 3D model-preconditioning FWI P.T. Trinh 1,2, R. Brossier 1, L. Métivier 1, L. Tavard 1, J. Virieux 1 and P. Wellington* 1 ¹ ISTerre/LJK/GRICAD Univ. Grenoble Alpes, CNRS
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 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 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 informationNSE 1.7. SEG/Houston 2005 Annual Meeting 1057
Finite-difference Modeling of High-frequency Rayleigh waves Yixian Xu, China University of Geosciences, Wuhan, 430074, China; Jianghai Xia, and Richard D. Miller, Kansas Geological Survey, The University
More informationAVO modeling in TESSERAL package
AVO modeling in TESSERAL package AVO modeling can be performed for complex media: isotropic, anisotropic, absorbing, gradient, thin-layered, for plain and curved reflecting boundaries. Package estimates
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (207) 209, 489 52 Advance Access publication 207 March 0 GJI Seismology doi: 0.093/gji/ggx03 Bessel smoothing filter for spectral-element mesh P.T. Trinh,
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 informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
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 information2012 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
EFFECTS OF 3D SURFACE TOPOGRAPHY ON REGIONAL AND TELESEISMIC SIGNALS FROM UNDERGROUND EXPLOSIONS Thorne Lay 1, Megan Avants 1, Xiao-bi Xie 1, and Arthur J. Rodgers 2 University of California Santa Cruz
More informationTh C 02 Model-Based Surface Wave Analysis and Attenuation
Th C 02 Model-Based Surface Wave Analysis and Attenuation J. Bai* (Paradigm Geophysical), O. Yilmaz (Paradigm Geophysical) Summary Surface waves can significantly degrade overall data quality in seismic
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) 191, 1391 1402 doi: 10.1111/j.1365-246X.2012.05686.x Optimized discrete wavelet transforms in the cubed sphere with the lifting scheme implications
More informationNUMERICAL MODELING OF STRONG GROUND MOTION USING 3D GEO-MODELS
th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 24 Paper No. 92 NUMERICAL MODELING OF STRONG GROUND MOTION USING D GEO-MODELS Peter KLINC, Enrico PRIOLO and Alessandro
More informationDownloaded 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 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 informationMigration from a non-flat datum via reverse-time extrapolation
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 50 Migration from a non-flat datum via reverse-time extrapolation Gopal Palacharla 1 ABSTRACT Land surveys usually have elevation changes,
More informationSMSIM Calibration for the SCEC Validation Project NGA-East Workshop Berkeley 12/11/12
SMSIM Calibration for the SCEC Validation Project NGA-East Workshop Berkeley 12/11/12 SWUS GMC SSHAC Level 3 Carola Di Alessandro GeoPentech with inputs from D.M. Boore OUTLINE Introduce SMSIM (Stochastic
More informationAn explicit feature control approach in structural topology optimization
th World Congress on Structural and Multidisciplinary Optimisation 07 th -2 th, June 205, Sydney Australia An explicit feature control approach in structural topology optimization Weisheng Zhang, Xu Guo
More informationAccelerating 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 informationWaveform tomography in geophysics and helioseismology
26 Waveform tomography in geophysics and helioseismology laura j. cobden, andreas fichtner, and vincent c. h. tong 26.1 What is waveform tomography? Seismic tomography in which we construct images of a
More informationVerification and Validation for Seismic Wave Propagation Problems
Chapter 26 Verification and Validation for Seismic Wave Propagation Problems (1989-2-24-25-28-29-21-211-217-) (In collaboration with Dr. Nima Tafazzoli, Dr. Federico Pisanò, Mr. Kohei Watanabe and Mr.
More informationComparison of 3D and 1D Wave Propagation Effects in the San Francisco Bay Area
Comparison of 3D and 1D Wave Propagation Effects in the San Francisco Bay Area Arben Pitarka SWUS GMC Workshop 2 October, 2013 LLNL-PRES-XXXXXX This work was performed under the auspices of the U.S. Department
More informationA A14 A24 A34 A44 A45 A46
Computation of qp- and qs-wave Rays raveltimes Slowness Vector and Polarization in General Anisotropic Media Junwei Huang 1 * Juan M. Reyes-Montes 1 and R. Paul Young 2 1 Applied Seismology Consultants
More informationThe Finite Element Method
The Finite Element Method A Practical Course G. R. Liu and S. S. Quek Chapter 1: Computational modeling An overview 1 CONTENTS INTRODUCTION PHYSICAL PROBLEMS IN ENGINEERING COMPUTATIONAL MODELLING USING
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 informationHierarchical Hybrid Grids
Hierarchical Hybrid Grids IDK Summer School 2012 Björn Gmeiner, Ulrich Rüde July, 2012 Contents Mantle convection Hierarchical Hybrid Grids Smoothers Geometric approximation Performance modeling 2 Mantle
More informationA COMPARISON OF SCATTERING FROM 2-D AND 3-D ROUGH INTERFACE. Ningya Cheng. Craig A. Schultz. M. N afi Toksoz
A COMPARISON OF SCATTERING FROM 2-D AND 3-D ROUGH INTERFACE Ningya Cheng Earth Resources Laboratory Department of Earth, Atmospheric, and Planetary Sciences Massachusetts Institute of Technology Cambridge,
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationby Chris Scheingraber, Kasra Hosseini, Robert Barsch, and Karin Sigloch
Electronic Supplement to ObsPyLoad - a tool for fully automated retrieval of seismological waveform data doi: 10.1785/0220120103 by Chris Scheingraber, Kasra Hosseini, Robert Barsch, and Karin Sigloch
More informationRegional tomographic inversion of the amplitude and phase of Rayleigh waves with 2-D sensitivity kernels
Geophys. J. Int. (26) 66, 48 6 doi:./j.365-246x.26.2972.x Regional tomographic inversion of the amplitude and phase of Rayleigh waves with 2-D sensitivity kernels Yingjie Yang and Donald W. Forsyth Department
More informationRay and wave methods in complex geologic structures
JOURNAL OF THE BALKAN GEOPHYSICAL SOCIETY, Vol. 5, No 2, May 2002, p. 57-64, figs. 13 Ray and wave methods in complex geologic structures Antonios P. Vafidis 1, Ioannis F. Louis 2 and George A. Karantonis
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 informationA 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 informationSeismic Soil-Structure Interaction Analysis of the Kealakaha Stream Bridge on Parallel Computers
Seismic Soil-Structure Interaction Analysis of the Kealakaha Stream Bridge on Parallel Computers Seung Ha Lee and Si-Hwan Park Department of Civil and Environmental Engineering University of Hawaii at
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 informationDevelopment of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics
Development of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics I. Pantle Fachgebiet Strömungsmaschinen Karlsruher Institut für Technologie KIT Motivation
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 informationSource 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 information2012 Tesseral Technologies
In the past, full-wavefield seismic modeling was not widely utilized in geophysical studies due to the lack of PC computational power. Less computationally intensive and simplified methods of seismic modeling
More informationTu P13 08 A MATLAB Package for Frequency Domain Modeling of Elastic Waves
Tu P13 8 A MATLAB Package for Frequency Domain Modeling of Elastic Waves E. Jamali Hondori* (Kyoto University), H. Mikada (Kyoto University), T.N. Goto (Kyoto University) & J. Takekawa (Kyoto University)
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 informationGeometric 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 informationRapid Processing of Synthetic Seismograms using Windows Azure Cloud
Rapid Processing of Synthetic Seismograms using Windows Azure Cloud Vedaprakash Subramanian, Liqiang Wang Department of Computer Science University of Wyoming En-Jui Lee, Po Chen Department of Geology
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 3, 2017, Lesson 1
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 3, 2017, Lesson 1 1 Politecnico di Milano, February 3, 2017, Lesson 1 2 Outline
More informationRefraction Full-waveform Inversion in a Shallow Water Environment
Refraction Full-waveform Inversion in a Shallow Water Environment Z. Zou* (PGS), J. Ramos-Martínez (PGS), S. Kelly (PGS), G. Ronholt (PGS), L.T. Langlo (PGS), A. Valenciano Mavilio (PGS), N. Chemingui
More informationFermat s principle and ray tracing in anisotropic layered media
Ray tracing in anisotropic layers Fermat s principle and ray tracing in anisotropic layered media Joe Wong ABSTRACT I consider the path followed by a seismic signal travelling through velocity models consisting
More informationMC 1.2. SEG/Houston 2005 Annual Meeting 877
PS-wave moveout inversion for tilted TI media: A physical-modeling study Pawan Dewangan and Ilya Tsvankin, Center for Wave Phenomena, Colorado School of Mines (CSM), Mike Batzle, Center for Rock Abuse,
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More information3D Finite Difference Time-Domain Modeling of Acoustic Wave Propagation based on Domain Decomposition
3D Finite Difference Time-Domain Modeling of Acoustic Wave Propagation based on Domain Decomposition UMR Géosciences Azur CNRS-IRD-UNSA-OCA Villefranche-sur-mer Supervised by: Dr. Stéphane Operto Jade
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 informationSEG/New Orleans 2006 Annual Meeting
Accuracy improvement for super-wide angle one-way waves by wavefront reconstruction Ru-Shan Wu* and Xiaofeng Jia, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary To
More informationLarge-scale workflows for wave-equation based inversion in Julia
Large-scale workflows for wave-equation based inversion in Julia Philipp A. Witte, Mathias Louboutin and Felix J. Herrmann SLIM University of British Columbia Motivation Use Geophysics to understand the
More informationDavid Wagner, Kaan Divringi, Can Ozcan Ozen Engineering
Internal Forces of the Femur: An Automated Procedure for Applying Boundary Conditions Obtained From Inverse Dynamic Analysis to Finite Element Simulations David Wagner, Kaan Divringi, Can Ozcan Ozen Engineering
More informationSimulation of seismic wave propagation in an asteroid based upon an unstructured MPI spectral-element method: blocking and non-blocking communication strategies Roland Martin 1, Dimitri Komatitsch 1,2,
More informationThe goal of full-waveform inversion (WFI) is to derive
SPECIAL Full SECTION: waveform F u l l w ainversion v e f o r m i n v e r s i o n Waveform inversion including well constraints, anisotropy, and attenuation CHAO WANG, DAVID YINGST, JIANYONG BAI, JACQUES
More informationLARGE-SCALE NORTHRIDGE EARTHQUAKE SIMULATION USING OCTREE-BASED MULTIRESOLUTION MESH METHOD
16th ASCE Engineering Mechanics Conference July 16-18, 3, University of Washington, Seattle LARGE-SCALE NORTHRIDGE EARTHQUAKE SIMULATION USING OCTREE-BASED MULTIRESOLUTION MESH METHOD Eui Joong Kim 1 Jacobo
More informationNOISE PROPAGATION FROM VIBRATING STRUCTURES
NOISE PROPAGATION FROM VIBRATING STRUCTURES Abstract R. Helfrich, M. Spriegel (INTES GmbH, Germany) Noise and noise exposure are becoming more important in product development due to environmental legislation.
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 informationSUMMARY NOISE REDUCTION BY SPATIAL SLOPE ESTIMATION
: steerable and non-linear filter approach Abdulaziz M. Almuhaidib and M. Nafi Toksöz, Earth Resources Laboratory, MIT SUMMARY We present an approach based on local-slope estimation for the separation
More informationParallel 3-D simulation of seismic wave propagation in heterogeneous anisotropic media: a grid method approach
Geophys. J. Int. (6) 165, 875 888 doi: 1.1111/j.1365-46X.6.946.x arallel 3-D simulation of seismic wave propagation in heterogeneous anisotropic media: a grid method approach Hongwei Gao and Jianfeng Zhang
More information