A A14 A24 A34 A44 A45 A46
|
|
- Edwin Hunter
- 6 years ago
- Views:
Transcription
1 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 Ltd 5 Claremont Buildings Claremont Bank Shrewsbury SY1 1RJ UK 2 University of oronto 170 College Street oronto Ontario M5S 3E3 Canada *junwei@appliedseismology.com Summary Based on Lax-Friedrichs scheme we developed a fast sweeping method of calculating the direct traveltime field for both quasi-p- and quasi-s-wave in general anisotropic media. he ray vector and the polarization can be obtained from the slowness vector at any given points in the model and the seismic ray path can be traced following the opposite direction of ray vectors. We illustrate the efficiency and accuracy of this algorithm using 2-D and 3-D numerical examples. Introduction A variety of mechanisms can cause seismic anisotropy such as alignments of mineral grains preferential orientation of fractures sequences of fine lithological layers and non-hydrostatic stress. In particular strong layering in shale deposits causes a non-negligible deviation from homogeneous isotropic velocity models when processing seismic data from this type of lithology. Accounting anisotropy in both active and passive seismic data processing has thus become a routine and critical step. he high frequency asymptotics for wave propagation in elastic solids reduces the isotropic seismic wave equation to a type of static Hamilton-Jacobi equation the so-called eikonal equation simply expressed as (1) where the operator takes the gradient of the traveltime field. calculates its norm and c is the propagation speed. In a general anisotropic medium however there is no simple explicit expression for the traveltime field and approximations such as weak anisotropy (homsen 1986) and first-order perturbation (Jech and Pšenčik 1989) must be incorporated to simplify the problem. In addition although it is convex in terms of slowness for qp-wave the anisotropic eikonal equation for one of the qs-waves is non-convex and conventional methods (e.g. sai et al. 2003; Sethian and Vladimirsky 2003; Mensch and Farra 1999) can only be used to solve the traveltime for qp-wave. Kao et al. (2004) presented the first fast numerical method based on discretizing the Hamilton-Jacobi equation directly without assuming convexity of the Hamiltonian and thus can be used to solve both qp- and qs-wave traveltime. In addition unlike seismic ray tracing Lax-Friedrichs scheme based fast sweeping method only needs the eigenvalues (Musgrave 1970) to update the traveltime field iteratively. In the following sections after showing that the analytical expressions of eikonal equations from Christoffel matrix are available only for isotropic and transversely isotropic medium we provide the workflow of solving the general anisotropic eikonal equation using Lax-Friedrichs sweeping scheme. We demonstrate the accuracy and convergence rate of our method on a 2-D isotropic medium and compare the three vectors: ray slowness and polarization in a 3-D homogeneous triclinic medium. Analytical Eikonal Equations he general anisotropic medium is fully characterized by 21 independent elastic parameters. Following the Voigt recipe the four-order density normalized elastic tensor can be expressed by a 6x6 symmetric matrix: A11 A12 A13 A14 A15 A16 A12 A22 A23 A24 A25 A 26. (2) A13 A23 A33 A34 A35 A36 A A14 A24 A34 A44 A45 A46 A15 A25 A35 A45 A55 A 56 A16 A26 A36 A46 A56 A66 he element of the Christoffel matrix Γ 3x3 are given by where p 1 p 2 p 3 are the three component of the slowness vector p i.e. the direction of the phase velocity and normal
2 Computation of Rays raveltimes Slowness and Polarization in General Anisotropic Medium to the wavefront. he matrix has three eigenvalues G m (xp) and eigenvectors g m (xp) with m=1 for qp-wave and m=23 for two qs-waves. If the slowness vector is known at a given point x in the model the ray vector r = [r 1 r 2 r 3 ] generally not in parallel with the slowness vector can be calculated as (Červenỳ 2001). (3) Using the element values of isotropic and transversely isotropic medium A33 A33 2A44 A33 2A A33 2A44 A33 A33 2A (4) A33 2A44 A33 2A44 A A A A A44 A11 A11 2A66 A A11 2A66 A11 A (5) A13 A13 A A A A A66 in Christoffel matrix Γ 3x3 we obtain the following eigenvalues and eigenvectors respectively: P P S S 0 G A p p p g p p p S G A p p p g p 0 p (6a) G A p p p g p p S L N GqP M N gqp p1 p3 p2 p A A L N Gq S1 M N gq S1 p1 p3 p2 p A A G A p A p p g p p qs qs 2 where L ( A33 A44 ) p3 ( A11 A44 ) p1 ( A11 A44 ) p2 M ( A A ) p ( A A ) p ( A A ) p N A A A A p A A A A p A 2A A A p p A 2A A A p (6b) A11 4A13 A33 A44 A11 A33 2A13 A44 p p 3. A11 4A13 A33 A44 A11 A33 2A13 A44 p2 In the isotropic case with letting the square root of the eigenvalue in (6a) equal to one we obtain eq. (1) with c=v p for P-wave and c=v s for S-wave. he ray direction is identical to the slowness direction and the polarization of P-wave is orthogonal to the polarization of S-wave. It is easy to verify that using the analytical expressions in (6b) the angle between the ray vector slowness vector and the polarization for qp-wave is nonzero and form an orthogonal coordinate. We will quantify the angle difference in the numerical example section. Obtaining analytical expressions for other less symmetrical anisotropy (e.g. orthorhombic) becomes significantly challenging. he Lax-Friedrichs sweeping scheme only requires the value of and thus we can numerically solve the general anisotropic eikonal equation (7) with a point source at x 0 and. Lax-Friedrichs Sweeping Algorithm Kao et al. (2004) gave a workflow of solving general static Hamilton-Jacobi equations in 2-D medium. Here we provide a workflow specifically tailored for eq. (7) in 3-D. 1. Initialization. We assign zero value to (x 0 ) at the source grid point and a large positive value to other grid points where the value should be larger than the possible largest arrival time. 2. Alternating Sweeping. At each iteration we perform eight alternating sweeps covering all the inner grid points whereas the grid points at the edge will be calculated in step 3. Each sweep is independent and can be easily parallelized on a shared memory machine. 3. Enforcing computational boundary conditions. After each sweep we enforce computational boundary condition. In our case for a point source inside the model there must be only outflow energy at the computational boundary. 4. Convergence check. Given the convergence criterion ε > 0 we check if is satisfied where n is the iteration number. Otherwise we return to step 2. At step 2 the following formulae are used to update the traveltime at grid point (ijk): i j k i 1 j k i 1 j k i j1 k i j1 k i j k 1 i j k 1 1 Gm x (8a) x y z i 1 j k i 1 j k i j1 k i j1 k i j k 1 i j k 1 x 2x y 2y z 2z n 1 n min (8b) i j k i j k i j k where 1 x y z and σ x σ y σ z are artificial viscosities satisfying x y z. (8c) According to our experience σ x = σ y = σ z = 5.0 is adequate for qp-wave and σ x = σ y = σ z = 3.0 for qs-wave. o enforce the outflow computational boundary condition in 3-D we first update the grid points on boundary planes
3 Computation of Rays raveltimes Slowness and Polarization in General Anisotropic Medium using the inner grid points and then the grid points on the edges using grid points inside the planes. For example assuming N x N y N z are the grid points along X Y and Z axes i=0 1 N x -1 j=0 1 N y -1 and k=0 1 N z -1 on the bottom XY boundary plane (k=0) we impose the following conditions: min max 2 (9a) i j0 i j1 i j2 i j2 i j0 min max 2 0 j0 1 j0 2 j0 2 j0 0 j0 min max 2 N x 1 j0 N x 2 j0 N x 3 j0 N x 2 j0 N x 1 j0 min max 2 i00 i10 i20 i20 i00 i N y 10 i N y 20 i N y 30 i N y 30 i N y 10 min max 2 (9b) (9c) Eq. (9a) ensures the energy at points on XY plane is from inside of the model and eq. (9b) and (9c) ensure the energy at the edge points is from inside of the XY plane. Equations (9) are trivially modified for computational boundary conditions on other planes. he details can be found in Kao et al. (2004). o evaluate its numerical efficiency and stability we apply the algorithm to a large 3-D model with strong triclinic anisotropy. o increase the anisotropy the elastic stiffness tensor is modified from a weak triclinic medium (Mah and Schmitt 2003) (10) A where only the upper part of the symmetric matrix is shown and the plus/minus signs mark the modification of the original weak triclinic elastic tensor. he density is assumed to be 1000 kg/m 3 for simplicity. Numerical Examples We first apply the algorithm to eq. (1) with c=1. A single point source is located at the center (x 0 y 0 ) = ( ) of the 2-D modeling region 1.0 x 1.0. he exact solution can be used to evaluate the convergence order and accuracy of the method. able 1 shows the apparent convergence order and L error as the mesh size approaches zero. able 1: he 2-D isotropic eikonal equation case: the error decreases with mesh size linearly (first-order). he convergence order p is calculated from and. he mesh size in b is selected in all other cases as a trade-off between computation time and accuracy. Mesh size L -error (e) Convergence Order (p) 1/ / / / / / able 2: he numerical efficiency and stability of solving the traveltime field for both qp- and qs-wave in a strong triclinic medium. Wave σ x=σ y=σ z CPU time (sec) Iteration No. qp qs-fast qs-slow Figure 1: he isosurfaces of the traveltime field for both qp- and qs-wave in a 3-D anisotropic medium. Significant elongation of the wavefront and shear wave splitting are consistent with the strong elastic stiffness assigned to the elastic tensor (see eq. 10). Note that the isosurface values for qp- and qs-waves are different for display purpose. able 3: he numerical efficiency and stability of solving the traveltime field for both qp- and qs-wave in a weak triclinic medium. Wave σ x=σ y=σ z CPU time (sec) Iteration No. qp qs-fast qs-slow able 2 presents the CPU time and number of iterations for solving the traveltime on 8 million grid points (200 x 200 x 200) for qp qs-fast and qs-slow waves on a single core. he traveltime will converge after several iterations and the isosurface of the traveltime field are shown in Figure 1. We found that increasing artificial viscosity values would increase the number of iterations and a trade-off should be found between a small viscosity value and the stability of the algorithm. In this study we use transversely isotropic
4 Computation of Rays raveltimes Slowness and Polarization in General Anisotropic Medium reference medium to derive the first-order analytical expression of anisotropic eikonal equation similar to Mensch and Farra (1999) but for both qp- and qs-waves. We can then estimate a feasible artificial viscosity value and choose the smaller one between the user defined (5 for qp-wave and 3 for qs-wave) and the estimated value. Next we solve the traveltime and derive slowness vector ray vector and polarization for a more realistic anisotropic medium i.e. the original weak triclinic medium shown in eq. (10) without modification. he CPU time is much less than the strong triclinic case (see able 3). he angle difference between the slowness ray and polarization is generally nonzero (Figure 2) and better quantified in Figure 3 for qp-wave. In general homogeneous-anisotropic media the angle difference between the slowness and ray and the slowness and polarization can be up to 20 o whereas the angle difference between the ray and the polarization is much smaller i.e. ~5 o. Figure 3: he angle difference as a function of θ and Φ on the isosurface shown in Figure 2b. In this triclinic case the deviation of the ray vector from the polarization is smaller than that of slowness vector. Figure 4: An example showing the irreversibility of the slowness path in a weak triclinic medium. As the source (S1) and the receiver (R1) switches positions the slowness path following the gradient of the traveltime field is irreversible. o trace the seismic ray path the ray vector must be calculated at each step from the slowness vector and as expected the ray path is reversible. Conclusions Figure 2: (a) he isosurfaces of qp qs-fast and qs-slow wavefronts in a weak triclinic medium. (bcd) show the angle difference between the slowness vector (blue arrow) ray vector (black arrow) and the polarization (red arrow) for grid point on an isosurface. In general the angle difference is nonzero. In applications where sources significantly outnumber receivers (e.g. Huang et al. 2013) receivers are usually regarded as virtual sources and the seismic reciprocity is applied to obtain the traveltime and ray path from an arbitrary grid to the receivers as if the source is located at the grid. In anisotropic medium care must be taken when tracing the ray path from the traveltime field since the ray direction must be calculated from the slowness at each step using eq. (3). Figure 4 gives an example of two paths following slowness vector and ray vectors. As the source and receiver position switches ray path is reversible whereas the slowness path is not as it is the phase velocity direction instead of the group velocity (energy) direction. In this paper we present a Lax-Friedrichs scheme-based fast sweeping method to calculate the direct traveltime field for both qp- and qs-waves in general anisotropic media. We give the workflow of updating traveltime and enforcing computational boundary conditions for 3-D cases. he convergence test on a 2-D example verifies its accuracy and the linear convergence rate (first order). he numerical efficiency and stability are demonstrated using a strong triclinic model. We quantify the difference between three vectors: slowness ray and polarization in a weak triclinic model and show that the ray path is reversible with smaller angle deviation from the polarization than the slowness. Acknowledgement We thank the members of staff at Applied Seismology Consultants Ltd for their support in this work. he authors are grateful to Dr. Leon homsen for the communication regarding anisotropy.
5 Computation of Rays raveltimes Slowness and Polarization in General Anisotropic Medium REFERENCE Červenỳ V Seismic ray theory Cambridge University Press 725 pages. Huang J. J. Reyes-Mountes and R.P. Young 2013 Automated microseismic event location using finite difference traveltime calculation and enhanced waveform stacking 75th EAGE Conference & Exhibition incorporating SPE EUROPEC 2013 London UK. Jech J. and Pšenčik I First-order approximation method for anisotropic media Geophys. J. Int Kao Y.C. S. Osher and J. Qian 2004 Lax- Friedrichs sweeping scheme for static Hamilton-Jacobi equations Journal of Computational Physics Mensch. and Farra V Computation of qpwave rays traveltimes and slowness vectors in orthorhombic media. Geophysical Journal International 138: doi: /j x x Mah M. and D.R. Schmitt 2003 Determination of the complete elastic stiffnesses from ultrasonic phase velocity measurements Journal of Geophysical Research Vol 108 No. B doi: /2001jb Musgrave M.J.P Crystal Acoustics: Introduction to the Study of Elastic Waves and Vibrations in Crystals Hen-Day San Francisco 288 pages. Sethian J. and A. Vladimirsky 2003 Ordered Upwind Methods for Static Hamilton--Jacobi Equations: heory and Algorithms SIAM Journal on Numerical Analysis 41: homsen L Weak elastic anisotropy: Geophysics sai Y.-H.R. L.-. Cheng S. Osher H.-K. Zhao 2003 Fast sweeping algorithms for a class of Hamilton Jacobi equations SIAM J. Numer. Anal. 41 (2)
Fermat 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 informationA comparison of exact, approximate, and linearized ray tracing methods in transversely isotropic media
Ray Tracing in TI media A comparison of exact, approximate, and linearized ray tracing methods in transversely isotropic media Patrick F. Daley, Edward S. Krebes, and Laurence R. Lines ABSTRACT The exact
More informationHTI anisotropy in heterogeneous elastic model and homogeneous equivalent model
HTI anisotropy in heterogeneous elastic model and homogeneous equivalent model Sitamai, W. Ajiduah*, Gary, F. Margrave and Pat, F. Daley CREWES/University of Calgary. Summary In this study, we compared
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 informationKey words. Hamilton-Jacobi, Eikonal, viscosity solution, Fast Marching, Fast Sweeping
COMPUTATIONAL STUDY OF FAST METHODS FOR THE EIKONAL EQUATION PIERRE A. GREMAUD AND CHRISTOPHER M. KUSTER Abstract. A computational study of the Fast Marching and the Fast Sweeping methods for the Eikonal
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 informationCOMPUTATIONALLY EFFICIENT RAY TRACING ALGORITHM FOR SIMULATION OF TRANSDUCER FIELDS IN ANISOTROPIC MATERIALS
Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation NDE 2011, December 8-10, 2011 COMPUTATIONALLY EFFICIENT RAY TRACING ALGORITHM FOR SIMULATION OF TRANSDUCER FIELDS IN ANISOTROPIC
More informationAVO for one- and two-fracture set models
GEOPHYSICS, VOL. 70, NO. 2 (MARCH-APRIL 2005); P. C1 C5, 7 FIGS., 3 TABLES. 10.1190/1.1884825 AVO for one- and two-fracture set models He Chen 1,RaymonL.Brown 2, and John P. Castagna 1 ABSTRACT A theoretical
More informationUNIVERSITY OF CALGARY. Traveltime Tomography in Isotropic and Transversely Isotropic Media. Marco Perez
Important Notice This copy may be used only for the purposes of research and private study, and any use of the copy for a purpose other than research or private study may require the authorization of the
More informationEstimating interval shear-wave splitting from multicomponent virtual shear check shots
GEOPHYSICS, VOL. 73, NO. 5 SEPTEMBER-OCTOBER 2008 ; P. A39 A43, 5 FIGS. 10.1190/1.2952356 Estimating interval shear-wave splitting from multicomponent virtual shear check shots Andrey Bakulin 1 and Albena
More informationA note on two-point paraxial travel times
Vlastislav Červený Charles University in Prague, Faculty of Mathematics and Physics, Department of Geophysics, Ke Karlovu 3, 6 Praha, Czech Republic vcerveny@seis.karlov.mff.cuni.cz Received: October,
More informationFoolproof AvO. Abstract
Foolproof AvO Dr. Ron Masters, Geoscience Advisor, Headwave, Inc Copyright 2013, The European Association of Geoscientists and Engineers This paper was prepared for presentation during the 75 th EAGE Conference
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 informationRotation parameters for model building and stable parameter inversion in orthorhombic media Cintia Lapilli* and Paul J. Fowler, WesternGeco.
otation parameters for model building and stable parameter inversion in orthorhombic media Cintia Lapilli* and Paul J Fowler, WesternGeco Summary Symmetry groups commonly used to describe seismic anisotropy
More informationPath Planning and Numerical Methods for Static (Time-Independent / Stationary) Hamilton-Jacobi Equations
Path Planning and Numerical Methods for Static (Time-Independent / Stationary) Hamilton-Jacobi Equations Ian Mitchell Department of Computer Science University of British Columbia Outline Path Planning
More informationAzimuthal Anisotropy Investigations for P and S Waves: A Physical Modelling Experiment
zimuthal nisotropy Investigations for P and S Waves: Physical Modelling Experiment Khaled l Dulaijan, Gary F. Margrave, and Joe Wong CREWES Summary two-layer model was built using vertically laminated
More informationTheoretical Background for OpenLSTO v0.1: Open Source Level Set Topology Optimization. M2DO Lab 1,2. 1 Cardiff University
Theoretical Background for OpenLSTO v0.1: Open Source Level Set Topology Optimization M2DO Lab 1,2 1 Cardiff University 2 University of California, San Diego November 2017 A brief description of theory
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 informationRaypath and traveltime computations Main for 2D transversely isotropic media with Heading dipping symmetry axes Bing Zhou Authors Stewart Greenhalgh
Exploration Geophysics (2006) 37, 150 159 Raypath and traveltime computations Main for 2D transversely isotropic media with Heading dipping symmetry axes Bing Zhou Authors Stewart Greenhalgh Key Words:
More informationNUMERICAL SIMULATION OF IRREGULAR SURFACE ACOUSTIC WAVE EQUATION BASED ON ORTHOGONAL BODY-FITTED GRIDS
- 465 - NUMERICAL SIMULATION OF IRREGULAR SURFACE ACOUSTIC WAVE EQUATION BASED ON ORTHOGONAL BODY-FITTED GRIDS LIU, Z. Q. SUN, J. G. * SUN, H. * LIU, M. C. * GAO, Z. H. College for Geoexploration Science
More informationHigh order fast sweeping methods for static Hamilton-Jacobi equations
High order fast sweeping methods for static Hamilton-Jacobi equations Yong-Tao Zhang, Hong-Kai Zhao and Jianliang Qian 3 ABSTRACT We construct high order fast sweeping numerical methods for computing viscosity
More information3-D traveltime computation using Huygens wavefront tracing
GEOPHYSICS, VOL. 66, NO. 3 MAY-JUNE 2001); P. 883 889, 12 FIGS., 1 TABLE. 3-D traveltime computation using Huygens wavefront tracing Paul Sava and Sergey Fomel ABSTRACT Traveltime computation is widely
More informationMODEL OF ULTRASONIC RAY PATHS IN AN ANISOTROPIC WELD BASED ON RAY METHOD Q. Liu, H. Wirdelius Chalmers Lindholmen University College, Göteborg, Sweden
MODEL OF ULTRASONIC RAY PATHS IN AN ANISOTROPIC WELD BASED ON RAY METHOD Q. Liu, H. Wirdelius Chalmers Lindholmen University College, Göteborg, Sweden Abstract: Experiences from isotropic materials suggest
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 informationAnisotropic 3D Amplitude Variation with Azimuth (AVAZ) Methods to Detect Fracture-Prone Zones in Tight Gas Resource Plays*
Anisotropic 3D Amplitude Variation with Azimuth (AVAZ) Methods to Detect Fracture-Prone Zones in Tight Gas Resource Plays* Bill Goodway 1, John Varsek 1, and Christian Abaco 1 Search and Discovery Article
More informationA Fast Sweeping Method for Static Convex Hamilton Jacobi Equations
Journal of Scientific Computing, Vol. 3, Nos. /2, May 27 ( 27) DOI:.7/s95-6-924-6 A Fast Sweeping Method for Static Convex Hamilton Jacobi Equations Jianliang Qian, ong-tao Zhang, 2 and Hong-Kai Zhao 3
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
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 informationInterpolation within ray tubes - state of the art
Interpolation within ray tubes - state of the art Petr Bulant Department of Geophysics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Praha 2, Czech Republic, http://sw3d.cz/staff/bulant.htm
More informationElastic physical seismic modeling Don C. Lawton,ScottP. Cheadle, EricV. Gallantand MalcolmB. Bertram
Elastic physical seismic modeling Don C. Lawton,ScottP. Cheadle, EricV. Gallantand MalcolmB. Bertram ABSTRACT A physieal seismic modeling system for studying elastic waves in solid models has been established.
More informationJoint Anisotropic Parameter Inversion Based on Nonhyperbolic Moveout and Azimuthal Amplitude Variance
URTe: 2455484 Joint Anisotropic Parameter Inversion Based on Nonhyperbolic Moveout and Azimuthal Amplitude Variance Fangyu Li*, Jie Qi, Bin Lyu, Kurt J. Marfurt, University of Oklahoma. opyright 2 16,
More informationThe Interface Fresnel Zone revisited
The Interface Fresnel Zone revisited Nathalie Favretto-Cristini, Paul Cristini and Eric de Bazelaire ABSTRACT We determine the part of reflectors which actually affects the reflected wavefield, which is
More informationOutline. Level Set Methods. For Inverse Obstacle Problems 4. Introduction. Introduction. Martin Burger
For Inverse Obstacle Problems Martin Burger Outline Introduction Optimal Geometries Inverse Obstacle Problems & Shape Optimization Sensitivity Analysis based on Gradient Flows Numerical Methods University
More informationGEOPHYSICS. Time to depth conversion and seismic velocity estimation using time migration velocities
Time to depth conversion and seismic velocity estimation using time migration velocities Journal: Manuscript ID: Manuscript Type: Date Submitted by the Author: Geophysics draft Letters n/a Complete List
More informationBlood Particle Trajectories in Phase-Contrast-MRI as Minimal Paths Computed with Anisotropic Fast Marching
Blood Particle Trajectories in Phase-Contrast-MRI as Minimal Paths Computed with Anisotropic Fast Marching Michael Schwenke 1, Anja Hennemuth 1, Bernd Fischer 2, Ola Friman 1 1 Fraunhofer MEVIS, Institute
More informationDevelopment of a particle difference scheme for dynamic crack propagation problem
Development of a particle difference scheme for dynamic crack propagation problem *Young-Choel Yoon 1) and Kyeong-Hwan Kim ) 1) Dept. of Civil Engineering, Myongji College, Seoul 10-776, Korea ) Dept.
More informationSeparation of diffracted waves in TI media
CWP-829 Separation of diffracted waves in TI media Yogesh Arora & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines Key words: Diffractions, Kirchhoff, Specularity, Anisotropy ABSTRACT
More informationOrdered upwind methods for static Hamilton Jacobi equations
Ordered upwind methods for static Hamilton Jacobi equations J. A. Sethian* and A. Vladimirsky Department of Mathematics, University of California, Berkeley, CA 94720 Edited by Alexandre J. Chorin, University
More informationThe Level Set Method applied to Structural Topology Optimization
The Level Set Method applied to Structural Topology Optimization Dr Peter Dunning 22-Jan-2013 Structural Optimization Sizing Optimization Shape Optimization Increasing: No. design variables Opportunity
More informationA High-Order Fast-Sweeping Scheme for Calculating First-Arrival Travel Times with an Irregular Surface
Bulletin of the Seismological Society of America, Vol. 103, No. 3, pp. 2070 2082, June 2013, doi: 10.1785/0120120199 A High-Order Fast-Sweeping Scheme for Calculating First-Arrival Travel Times with an
More informationTh SBT3 03 An Automatic Arrival Time Picking Method Based on RANSAC Curve Fitting
3 May June 6 Reed Messe Wien Th SBT3 3 An Automatic Arrival Time Picking Method Based on RANSAC Curve Fitting L. Zhu (Georgia Institute of Technology), E. Liu* (Georgia Institute of Technology) & J.H.
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 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 informationm=[a,b,c,d] T together with the a posteriori covariance
zimuthal VO analysis: stabilizing the model parameters Chris Davison*, ndrew Ratcliffe, Sergio Grion (CGGVeritas), Rodney Johnston, Carlos Duque, Musa Maharramov (BP). solved using linear least squares
More informationP052 3D Modeling of Acoustic Green's Function in Layered Media with Diffracting Edges
P052 3D Modeling of Acoustic Green's Function in Layered Media with Diffracting Edges M. Ayzenberg* (Norwegian University of Science and Technology, A. Aizenberg (Institute of Petroleum Geology and Geophysics,
More informationPrestack Kirchhoff time migration for complex media
Stanford Exploration Project, Report 97, July 8, 998, pages 45 6 Prestack Kirchhoff time migration for complex media Tariq Alkhalifah keywords: time migration, anisotropy ABSTRACT Constructing the seismic
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 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 informationHigher-order schemes for 3-D first-arrival traveltimes and. amplitudes
Higher-order schemes for 3-D first-arrival traveltimes and amplitudes Songting Luo, Jianliang Qian 2 and Hongkai Zhao 3 Department of Mathematics, Michigan State University, East Lansing, MI 48824. Email:
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 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 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 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 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 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 informationThe azimuth-dependent offset-midpoint traveltime pyramid in 3D HTI media
The azimuth-dependent offset-midpoint traeltime pyramid in 3D HTI media Item Type Conference Paper Authors Hao Qi; Stoas Alexey; Alkhalifah Tariq Ali Eprint ersion Publisher's Version/PDF DOI.9/segam3-58.
More informationAnisotropic 3D Amplitude Variation with Azimuth (AVAZ) Methods to Detect Fracture Prone Zones in Tight Gas Resource Plays
Antropic 3D Amplitude Variation with Azimuth (AVAZ) Methods to Detect Fracture Prone Zones in Tight Gas Resource Plays Bill Goodway* EnCana Corp., Calgary, AB, Canada william.goodway@encana.com and John
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 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 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 informationFederal Institute for Materials Research and Testing (BAM), Unter den Eichen 87, Berlin, Germany
Jannis Bulling 1, Jens Prager 1, Fabian Krome 1 1 Federal Institute for Materials Research and Testing (BAM), Unter den Eichen 87, 12205 Berlin, Germany Abstract: This paper addresses the computation of
More informationAdvanced Image Reconstruction Methods for Photoacoustic Tomography
Advanced Image Reconstruction Methods for Photoacoustic Tomography Mark A. Anastasio, Kun Wang, and Robert Schoonover Department of Biomedical Engineering Washington University in St. Louis 1 Outline Photoacoustic/thermoacoustic
More informationModeling lattice structured materials with micropolar elasticity. Accuracy of the micropolar model. Marcus Yoder. CEDAR presentation Spring 2017
Modeling lattice structured materials with micropolar elasticity Accuracy of the micropolar model CEDAR presentation Spring 2017 Advisors: Lonny Thompson and Joshua D. Summers Outline 2/25 1. Background
More informationFast marching methods
1 Fast marching methods Lecture 3 Alexander & Michael Bronstein tosca.cs.technion.ac.il/book Numerical geometry of non-rigid shapes Stanford University, Winter 2009 Metric discretization 2 Approach I:
More informationElastic wave mode separation for TTI media
CWP-628 Elastic wave mode separation for TTI media Jia Yan and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT The separation of wave modes for isotropic elastic wavefields is typically
More informationFracture Quality from Integrating Time-Lapse VSP and Microseismic Data
Fracture Quality from Integrating Time-Lapse VSP and Microseismic Data Mark E. Willis, Daniel R. Burns, Rongrong Lu, M. Nafi Toksöz, Earth Resources Laboratory Dept. of Earth, Atmospheric, and Planetary
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 informationOptimizing Fully Anisotropic Elastic Propagation on 2nd Generation Intel Xeon Phi Processors. 2 Wave propagation in anisotropic elastic media
Optimizing Fully Anisotropic Elastic Propagation on 2nd Generation Intel Xeon Phi Processors 1 Introduction In the last few years, acoustic isotropic wave propagation has been the preferred simulation
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 informationMinimizing Fracture Characterization Uncertainties Using Full Azimuth Imaging in Local Angle Domain
P-237 Minimizing Fracture Characterization Uncertainties Using Full Azimuth Imaging in Local Angle Domain Shiv Pujan Singh*, Duane Dopkin, Paradigm Geophysical Summary Shale plays are naturally heterogeneous
More informationInversion of Fracture Parameters by Using the Artificial Neural Network
IC -AI 99 International Conference 57 Inversion of Fracture Parameters by Using the Artificial Neural Network Hengchang Dai, Xiang-Yang Li, and Colin MacBeth British Geological Survey, Murchison House,
More informationDynamic Programming Algorithms for Planning and Robotics in Continuous Domains and the Hamilton-Jacobi Equation
Dynamic Programming Algorithms for Planning and Robotics in Continuous Domains and the Hamilton-Jacobi Equation Ian Mitchell Department of Computer Science University of British Columbia research supported
More informationA hybrid tomography method for crosswell seismic inversion
hybrid tomography method A hybrid tomography method for crosswell seismic inversion Danping Cao and Wenyuan Liao SUMMARY The waveform tomography, or full waveform inversion (FWI), is in general better
More informationFast Marching Methods (1) are numerical algorithms for
Fast methods for the Eikonal and related Hamilton Jacobi equations on unstructured meshes J. A. Sethian* and A. Vladimirsky Department of Mathematics, University of California, Berkeley, CA 94720 Communicated
More informationSUMMARY INTRODUCTION THEORY. The anisotropic elastic wave equation is
Passive wavefield imaging using the energy norm Daniel Rocha and Paul Sava, Center for Wave Phenomena, Colorado School of Mines, Jeffrey Shragge, Centre for Energy Geoscience, The University of Western
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
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 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 informationComparison of the Parallel Fast Marching Method, the Fast Iterative Method, and the Parallel Semi-Ordered Fast Iterative Method
Procedia Computer Science Volume 0, 06, Pages 7 75 ICCS 06. The International Conference on Computational Science Comparison of the Parallel Fast Marching Method, the Fast Iterative Method, and the Parallel
More informationAn O(N) Level Set Method. Seongjai Kim. Abstract. Level set algorithms have been extensively applied for the problems
An O(N) Level Set Method Seongjai Kim Abstract A propagating interface can develop corners and discontinuities as it advances. Level set algorithms have been extensively applied for the problems in which
More informationSeismic physical modeling II: VVAZ and AVAZ effects observed on reflections from isolated HTI targets
Seismic physical modeling II: VVAZ and AVAZ effects observed on reflections from isolated HTI targets Joe Wong ABSTRACT Seismic physical modeling was used to investigate VVAZ and AVAZ phenomena associated
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 informationProcessing converted-wave data in the tau-p domain: rotation toward the source and moveout correction
τ-p domain converted-wave processing Processing converted-wave data in the tau-p domain: rotation toward the source and moveout correction Raul Cova and Kris Innanen ABSTRACT The asymmetry of the converted-wave
More informationLaurent D. Cohen 2 CEREMADE, Université Paris Dauphine PARIS CEDEX 16 - FRANCE
The shading zone problem in geodesic voting and its solutions for the segmentation of tree structures. Application to the segmentation of Microglia extensions Youssef Rouchdy 1,2 University of Pennsylvania
More informationFourth-order NMO velocity for P-waves in layered orthorhombic media vs. offset-azimuth
Fourt-order NMO velocity for P-waves in layered orrombic media vs. set-azimut Zvi Koren* and Igor Ravve Paradigm Geopysical Summary We derive te fourt-order NMO velocity of compressional waves for a multi-layer
More informationSEG/San Antonio 2007 Annual Meeting
Imaging steep salt flanks by super-wide angle one-way method Xiaofeng Jia* and Ru-Shan Wu, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary The super-wide angle one-way
More informationFull-waveform inversion for reservoir characterization: A synthetic study
CWP-889 Full-waveform inversion for reservoir characterization: A synthetic study Nishant Kamath, Ilya Tsvankin & Ehsan Zabihi Naeini ABSTRACT Most current reservoir-characterization workflows are based
More informationSummary. Introduction
Chris Davison*, Andrew Ratcliffe, Sergio Grion (CGGeritas), Rodney Johnston, Carlos Duque, Jeremy Neep, Musa Maharramov (BP). Summary Azimuthal velocity models for HTI (Horizontal Transverse Isotropy)
More informationTwo Semi-Lagrangian Fast Methods for Hamilton-Jacobi-Bellman Equations
Two Semi-Lagrangian Fast Methods for Hamilton-Jacobi-Bellman Equations Simone Cacace 1, Emiliano Cristiani 2(B), and Maurizio Falcone 1 1 Dipartimento di Matematica, Sapienza Università di Roma, Rome,
More informationdp i = u ray c2 dt with t the travel time.
Three dimensional two point ray tracing using paraxial rays in cartessian coordinates Ronit Strahilevitz*, Paradigm Geophysical; Dan Kosloff, Department of Geophysics, Tel Aviv University, Israel and Paradigm
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. 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 informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationA level set based Eulerian method for paraxial multivalued traveltimes
Journal of Computational Physics 97 (24) 7 736 www.elsevier.com/locate/jcp A level set based Eulerian method for paraxial multivalued traveltimes Jianliang Qian *, Shingyu Leung Department of Mathematics,
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 informationAVO Analysis with Multi-Offset VSP Data
AVO Analysis with Multi-Offset VSP Data Z. Li*, W. Qian, B. Milkereit and E. Adam University of Toronto, Dept. of Physics, Toronto, ON, M5S 2J8 zli@physics.utoronto.ca T. Bohlen Kiel University, Geosciences,
More informationBi-directional seismic vibration control of spatial structures using passive mass damper consisting of compliant mechanism
Bi-directional seismic vibration control of spatial structures using passive mass damper consisting of compliant mechanism Seita TSUDA 1 and Makoto OHSAKI 2 1 Department of Design, Okayama Prefectural
More informationLevel Set Methods and Fast Marching Methods
Level Set Methods and Fast Marching Methods I.Lyulina Scientific Computing Group May, 2002 Overview Existing Techniques for Tracking Interfaces Basic Ideas of Level Set Method and Fast Marching Method
More informationThe Conventional and Non-conventional Seismic Differencing
Summary The Conventional and Non-conventional Seismic Differencing Vanja Vracar* University of Calgary, CREWES Project, Alberta, Canada vmilicev@ucalgary.ca and Robert J. Ferguson University of Calgary,
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 information