SEG/New Orleans 2006 Annual Meeting. Non-orthogonal Riemannian wavefield extrapolation Jeff Shragge, Stanford University
|
|
- Tracey Wilkins
- 5 years ago
- Views:
Transcription
1 Non-orthogonal Riemannian wavefield extrapolation Jeff Shragge, Stanford University SUMMARY Wavefield extrapolation is implemented in non-orthogonal Riemannian spaces. The key component is the development of a dispersion relationship appropriate for propagating wavefields on generalized non-orthogonal meshes. This wavenumer contains a numer of mixed-domain fields in addition to velocity that represent coordinate system geometry. An extended split-step Fourier approximation of the extrapolation wavenumer is developed that provides accurate results when multiple reference parameters sets are used. Three examples are presented that demonstrate the validity of the theory. An important consequence is that greater emphasis can e placed on generating smoother computational meshes, rather than satisfying restrictive semi-orthogonal criteria, leading to more accurate and efficient Riemannian wavefield extrapolation. INTRODUCTION Riemannian wavefield extrapolation (RWE) (Sava and Fomel, 005) generalizes wavefield extrapolation to non-cartesian coordinate systems. The original implimentation assumed that coordinate systems are at least semi-orthogonal and characterized y an extrapolation direction orthogonal to the other two axes. This resulted in a wave-equation dispersion relationship for the extrapolation wavenumer containing mixed-domain fields additional to velocity that encode coordinate system geometry. However, semi-orthogonal geometry can e an overly restrictive assertion ecause many computational meshes have greatly varying mixed-domain coefficients that cause numerical instaility during wavefield extrapolation. Initially, RWE was designed for dynamic applications where wavefields are extrapolated on ray-ased coordinate systems oriented in the wave propagation direction. This approach generally generates high-quality Green s functions; however, numerical instaility (i.e. zero-division) occurs wherever the ray coordinate system triplicates. Sava and Fomel (005) address this issue y iteratively smoothing the velocity model until these triplications vanish. This solution, though, somewhat counters the original purpose of RWE: coordinate systems conformal to the propagation direction. A more geometric RWE application is performing wavefield extrapolation to and from surfaces of irregular geometry. Shragge and Sava (005) formulate a wave-equation migration from topography strategy that poses wavefield extrapolation directly in locally orthogonal meshes conformal to the acquisition surface. Although successful in areas with longer wavelength and lower amplitude topography, imaging results degrade in situations involving more rugged acquisition topography. However, a more general oservation is the genetic link etween degraded image quality and the grid compression/extension demanded y semi-orthogonal geometry. A solution to these prolems is eliminating the semi-orthogonal constraint y allowing wavefield propagation in non-orthogonal Riemannian coordinate spaces. Non-orthogonality introduces two additional terms in the RWE dispersion relationship and makes existing coefficient terms more complex. Non-stationary coefficients in the resulting extrapolation wavenumer can e handled with an extended split-step Fourier approach (Stoffa et al., 990). Importantly, this solution affords greater flexiility in coordinate system design while facilitating more rapid mesh generation. Furthermore, greater emphasis can e placed on optimizing grid quality y controlling grid clustering and generating smoother coefficient fields. This paper develops a D wave-equation dispersion relationship for performing non-orthogonal RWE. I first discuss generalized Riemannian geometry and show how the acoustic wave-equation can e formulated in a non-orthogonal Riemannian space. Susequently, I develop an expression for a one-way wavefield extrapolation wavenumer and present the corresponding split-step Fourier approximation. I then present two analytic -D non-orthogonal coordinate systems to help validate the developed extrapolation wavenumer expressions. The paper concludes with the example of RWE-generated Green s function estimates through a SEG- EAGE model section. GENERALIZED RIEMANNIAN GEOMETRY Geometry in a generalized D Riemannian space is descried y a symmetric metric tensor, g ij g ji, that relates the geometry in a non-orthogonal coordinate system, {x, x, x }, to an underlying Cartesian mesh, {ξ,ξ,ξ }. The matrix-form metric tensor is, gij g g g g g g g g g g g g, () g g g g where g, g, g, g, g and are functions linking the two coordinate systems through, ( e.g. () g ij x k x k g x ) k x k ξ ξ (Summation notation - e.g. g ii g + g + - is used throughout this paper.) The associated (or inverse) metric tensor, g ij, is defined y g ij g g ij, where g is metric tensor matrix determinant. The associated metric tensor is given y, g ij g g g g g g g g g g g g g g g g g g g, g g g g g g g g g g g () and has the following metric determinant, g (g g g ) g g + g g g g g (g g g ). (4) Weighted metric tensor, m ij g g ij, is a useful definition for the following development. GENERALIZED -D ACOUSTIC WAVE-EQUATION The acoustic wave-equation for wavefield, U, in a generalized Riemannian space is, U ω s (x)u, (5) where the ω is frequency, s is the propagation slowness, and is the Laplacian operator, U ( ) ij U m. (6) g Sustituting equation 6 into 5 generates a Helmholtz equation appropriate for propagating waves through a D space, ( ) ij U m ω s U. (7) g Expanding derivative terms and multiplying y g yields, m ij ξ i U ξ j + m ij U g ω s U. (8) 6
2 Defining n j as, n j mij ξ i m j ξ + m j ξ leads to a more compact notation of equation 8, + m j ξ, (9) U n j + m ij U g ω s U. (0) ξ j The -D wave-equation dispersion relation is developed y replacing the partial differential operators acting on wavefield U with their Fourier-domain duals, (m ij k ξi in j ) k ξj g ω s, () where k ξi is the dual of differential operator ξ. Equation represents the dispersion relationship for wavefield propagation on a i generalized -D Riemannian space. One-way wavefield extrapolation Developing an expression for the extrapolation wavenumer requires isolating one of the wavenumers in equation (herein assumed to e coordinate ξ ). Introducing indicies i, j,,, into equation and rearranging terms yields, m kξ + (m k ξ + m ) k ξ i n k ξ ) g ω s + i (n k ξ + n k ξ m kξ m k ξ m k ξ k ξ. () An expression for wavenumer k ξ can e otained y completing the square and isolating the k ξ contriutions, k ξ a k ξ a k ξ + ia () ± a4 ω a5 k ξ a 6 k ξ a 7 k ξ k ξ + a 8 i k ξ + a 9 i k ξ a0, where a i are non-stationary coefficients given y, a g g n m s g ( g ) g ( g g g n ( ) m m n n ( m ) m m n ( m ) g )... T n m. (4) where i i (ξ ) are reference values of a i a i (ξ,ξ,ξ ). The split-step approximation is developed y performing a Taylor expansion aout each coefficient a i and evaluating the results at stationary reference values i. The stationary values of k ξ and k ξ are assumed to e zero. This leads to a split-step correction of, kξ SSF k ξ ( ) ( ) k a a + ξ ( ) 0 a a a a0 0, (6) 0 0 where 0 denotes "with respect to a reference medium". The partial differential expressions in equation 6 are, a, 0 a 4 ω, ω 0 a 0, ω 0 (7) resulting in the following split-step Fourier correction, kξ SSF i (a ) + 4 ω (a 4 4 ) 4 ω 0 0 (a 0 0 ). (8) 4 ω 0 Note that we could use many reference media followed y interpolation similar to the phase-shift plus interpolation (PSPI) technique of Gazdag and Sguazzero (984). Importantly, even though there are additional a i coefficients in the dispersion relationship, these can e made smooth through mesh regularization such that fewer sets of reference parameters are needed to accurately represent wavenumer k ξ. In addition, situations exist where some coefficients are zero or negligile. For example, four coefficients adequately descrie a weakly non-orthogonal coordinate system (i.e. max {g, g, g } << max{g, g, }) within a kinematic (i.e. no imaginary terms) approximation. EXAMPLE - -D SHEARED CARTESIAN COORDINATES An instructive example is a coordinate system is a -D sheared Cartesian mesh defined y, x cosθ x 0 sinθ ξ ξ, (9) where θ is the shear angle of the coordinate system (θ 90 is Cartesian). This system has a metric tensor g ij given y, xk x k gij ξ ξ x k ξ x k ξ x k x k ξ ξ x k x k ξ ξ cosθ cosθ, (0) Note that the coefficients contain a mixture of m ij and g ij terms, and that positive definite terms, a 4,a 5,a 6 and a 0 in equation, are squared such that the familiar Cartesian split-step Fourier correction is recovered elow. Split-Step Fourier Approximation The extrapolation wavenumer defined in equations and 4 cannot e implemented purely in the Fourier domain due to the presence of mixed-domain fields (i.e. a function of oth ξ and k ξ simultaneously). This can e addressed using an extended version of the split-step Fourier approximation (Stoffa et al., 990) that uses Taylor expansions to separate k ξ into two parts: k ξ kξ P S +k SSF ξ. Wavenumers k P S ξ and k SSF ξ represent a pure Fourier domain phase-shift and a mixed ω x domain split-step correction, respectively. The phase-shift term is given y, k P S ξ k ξ k ξ + i (5) ± 4 ω 5 k ξ 6 k ξ 7 k ξ k ξ + 8 i k ξ + 9 i k ξ 0, with determinant g sin θ and associated metric tensor g ij, g ij sin θ cosθ cosθ. () Note that ecause the tensor in equation is spatially invariant, equation 0 simplifies to, U g ij and generates the following dispersion relation, U ω s U, () g ij k ξi k ξj ω s. () Expanding terms leads to an expression for wavenumer k ξ, ( k ξ g k ξ ± s ω g ( ) g ) kξ. (4) 7
3 Sustituting the values of the associated metric tensor in equation into equation 4 yields, k ξ cosθ k ξ ± sinθ s ω k ξ. (5) A numerical test using a Cartesian coordinate system sheared at 5 from vertical is shown in figure. The ackground velocity model is 500 ms and the zero-offset data consist of 4 flat plane-waves t 0., 0.4, 0.6 and 0.8 s. As expected, the zero-offset migration results image reflectors at depths z00, 600, 900, and 00 m. Note that the propagation has created oundary artifacts: those on the left are reflections due to a truncated coordinate system while those on the right are hyperolic diffractions caused y truncated plane-waves. Depth (m) Distance (m) Figure : Polar ellipsoidal coordinate system example. Figure : Sheared Cartesian coordinate system test. Coordinate system shear angle and velocity are θ 5 and 500 ms, respectively. Zero-offset data consist of 4 flat plane-wave impulses at t 0.,0.4,0.6 and 0.8 s that are correctly imaged at depths z 00, 600, 900, and 00 m. EXAMPLE - POLAR ELLIPSOIDAL A second instructive example is a stretched polar coordinate system (see figure ). A polar ellipsoidal coordinate system is specified y, x a(ξ )ξ cosξ. (6) x a(ξ )ξ sinξ Parameter a a(ξ ) is a smooth function controlling coordinate system ellipticity and has curvature parameters ξ a and c a ξ. The metric tensor g ij is, gij a ξ a ξ a ξ ( + a ), (7) with determinant g a 4 ξ. The associated metric and weighted associated metric tensors are given y, g ij +a a 4 a ξ a ξ a ξ and m ij ξ ( +a ) a a a ξ (8) Tensors g ij and m ij specify a wavenumer appropriate for extrapolating wavefields on a -D non-orthogonal mesh. However, ecause the coordinate system is spatially variant, we must also. Figure : Ellipsoidal polar coordinate system test example. Upper left: v(z) z velocity function overlain y a polar ellipsoidal coordinate system defined y parameter a + 0.ξ 0.05ξ. Upper right: velocity model in the RWE domain. Bottom right: Imaged reflectors in RWE domain. Bottom left: the RWE domain image mapped to a Cartesian mesh. compute the n i fields: n a + ac a and n 0. Inserting these values yields the following extrapolation wavenumer k ξ, k ξ ξ ( a a k ξ ± a ξ s ω ξ k ξ ik ξ ξ + ) ac a. (9) Note that the kinematic version of equation 9 is, k ξ ξ a k ξ ± a s ω k ξ, (0) Figure shows an ellipsoidal polar coordinate system defined y a + 0.ξ 0.05ξ. The upper left panel shows a v(z) z velocity function overlain y the coordinate system mesh. The upper right panel presents velocity model as mapped into the RWE domain. The data used in this test consisted of 4 flat plane-waves. Given this experimental setup, propagated flat plane-waves should not end in the Cartesian domain ecause of the v v(z) velocity model, even though there is velocity variation across each extrapolation step in the RWE domain. Hence, the impulses have curvature in the RWE domain (lower right panel). The lower left panel shows the RWE domain imaging results mapped ack to a Cartesian domain. Consistent with theory, the flat planewaves are imaged as flat reflectors. Note that the edge effects are again caused y coordinate system and plane-wave truncation. 8
4 EXAMPLE - RWE GREEN S FUNCTION GENERATION The final test uses RWE to generate Green s function estimates. The test velocity model is a slice of the SEG-EAGE salt model (see figure 4). Velocity contrasts etween the salt ody and sediment cause complex wavefield propagation including significant wavefield triplication and multipathing. The upper left panel shows the velocity model with an overlain coordinate mesh generated y differential meshing (Liseikin, 004). The mesh is a traveltime-ased coordinate system where the first and last extrapolation steps are formed y the 0.04 s and.5 s Eikonal equation isochrons. The velocity model in the RWE domain is illustrated in the upper right. The lower right panel shows the impulse response test in the RWE domain. The 7 impulsive waves conform fairly well to the traveltime steps, except where they enter the salt ody to the lower left of the image. The migration results mapped ack to Cartesian space are shown in the lower left panel. The wavefield to the left of the shot point is fairly complicated and the energy in the salt ody and the corresponding upward refracted (and perhaps reflected?) wavefields are strongly present. Figure 4: Example of RWE generated Green s Functions on structured non-orthogonal mesh for a slice through the SEG-EAGE Salt data set velocity model. Top left: Salt model in physical space with an overlain ray-ased mesh. Top right: Velocity model in the RWE domain. Bottom right: Wavefield propagated in ray-coordinates through velocity model shown in the upper right. Bottom left: Wavefield in ottom right interpolated ack to Cartesian space. Figure 5 presents a comparison test etween RWE and Cartesian extrapolation. Beneath and right of the shot point the wavefields are fairly similar except for a phase-change. However, they are significantly different to the left ecause Cartesian-ased extrapolation neither propagates energy laterally with the same accuracy nor upward at all. Hence, this energy is asent from the wavefield in the lower panel. CONCLUDING REMARKS This paper discusses Riemannian wavefield extrapolation in D non-orthogonal coordinate systems. The extrapolation wavenumer is decoupled from the other wavenumers allowing for an extended split-step Fourier approximation. Examples indicate that wavefields can e extrapolated on non-orthogonal coordinate meshes. Hence, users can concentrate more on mesh design to optimizing mixed-domain field smoothness leading to more accurate RWE. ACKNOWLEDGEMENTS Thanks to Paul Sava and Sergey Fomel for helpful discussions. REFERENCES Gazdag, J. and P. Sguazzero, 984, Migration of seismic data y phase-shift plus interpolation: Geophysics, 49, 4. Liseikin, V., 004, A Computational Differential Geometry Approach to Grid Generation: Springer-Verlag, Berlin. Sava, P. and S. Fomel, 005, Riemannian wavefield extrapolation: Geophysics, 70, T45 T56. Shragge, J. and P. Sava, 005, Wave-equation migration from topography: 75th Ann. Internat. Mtg., Soc. of Expl. Geophys., Expanded Astracts, Stoffa, P. L., J. T. Fokkema, R. M. de Luna Freire, and W. P. Kessinger, 990, Split-step Fourier migration: Geophysics, 55, no. 04, Figure 5: Results of a comparison test etween generalized Riemannian wavefield extrapolation (top panel) and Cartesian-ased extrapolation (ottom panel). 9
5 EDITED REFERENCES Note: This reference list is a copy-edited version of the reference list sumitted y the author. Reference lists for the 006 SEG Technical Program Expanded Astracts have een copy edited so that references provided with the online metadata for each paper will achieve a high degree of linking to cited sources that appear on the We. REFERENCES Gazdag, J., and P. Sguazzero, 984, Migration of seismic data y phase-shift plus interpolation: Geophysics, 49, 4. Liseikin, V., 004, A Computational Differential Geometry Approach to Grid Generation: Springer-Verlag. Sava, P., and S. Fomel, 005, Riemannian wavefield extrapolation: Geophysics, 70, T45 T56. Shragge, J., and P. Sava, 005, Wave-equation migration from topography: 75th Annual International Meeting, SEG, Expanded Astracts, Stoffa, P. L., J. T. Fokkema, R. M. de Luna Freire, and W. P. Kessinger, 990, Split-step Fourier migration: Geophysics, 55,
Mesh generation using differential methods
Stanford Exploration Project, Report 24, April 4, 2006, pages 45?? Mesh generation using differential methods Jeff Shragge ABSTRACT This paper examines a differential gridding method for generating computational
More informationImplicit 3-D depth migration by wavefield extrapolation with helical boundary conditions
Stanford Exploration Project, Report 97, July 8, 1998, pages 1 13 Implicit 3-D depth migration by wavefield extrapolation with helical boundary conditions James Rickett, Jon Claerbout, and Sergey Fomel
More informationSeismic imaging using Riemannian wavefield extrapolation
Stanford Exploration Project, Report 114, October 14, 2003, pages 1 31 Seismic imaging using Riemannian wavefield extrapolation Paul Sava and Sergey Fomel 1 ABSTRACT Riemannian spaces are described by
More informationINTRODUCTION. The depth-variable-velocity prestack migration in offset-midpoint coordinates (Yilmaz, 1979) is formulated as
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 58 Short Note Prestack migration by split-step DSR Alexander M. Popovici 1 INTRODUCTION The DSR prestack migration equation, though defined
More informationMorse, P. and H. Feshbach, 1953, Methods of Theoretical Physics: Cambridge University
Bibliography Albertin, U., D. Yingst, and H. Jaramillo, 2001, Comparing common-offset Maslov, Gaussian beam, and coherent state migrations: 71st Annual International Meeting, SEG, Expanded Abstracts, 913
More informationFast wavefield extrapolation by phase-shift in the nonuniform Gabor domain
Fast wavefield extrapolation by phase-shift in the nonuniform Gabor domain Jeff P. Grossman* and Gary F. Margrave Geology & Geophysics, University of Calgary 2431 22A Street NW, Calgary, AB, T2M 3X8 grossman@geo.ucalgary.ca
More informationRiemannian wavefield extrapolation
GEOPHYSICS, VOL. 70, NO. 3 MAY-JUNE 2005); P. T45 T56, 14 FIGS. 10.1190/1.1925748 Riemannian wavefield extrapolation Paul Sava 1 and Sergey Fomel 2 ABSTRACT Riemannian spaces are described by nonorthogonal
More informationP S-wave polarity reversal in angle domain common-image gathers
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? P S-wave polarity reversal in angle domain common-image gathers Daniel Rosales and James Rickett 1 ABSTRACT The change in the reflection
More informationStable wide-angle Fourier finite-difference downward extrapolation of 3-D. wavefields
Stable wide-angle Fourier finite-difference downward extrapolation of 3-D wavefields Biondo Biondi Stanford Exploration Project, Mitchell Bldg., Department of Geophysics, Stanford University, Stanford,
More informationPlane-wave migration in tilted coordinates
Stanford Exploration Project, Report 124, April 4, 2006, pages 1 16 Plane-wave migration in tilted coordinates Guojian Shan and Biondo Biondi ABSTRACT Plane-wave migration in tilted coordinates is powerful
More informationSUMMARY ANGLE DECOMPOSITION INTRODUCTION. A conventional cross-correlation imaging condition for wave-equation migration is (Claerbout, 1985)
Comparison of angle decomposition methods for wave-equation migration Natalya Patrikeeva and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Angle domain common image gathers offer
More informationFast wavefield extrapolation by phase-shift in the nonuniform Gabor domain
Gabor wavefield extrapolation Fast wavefield extrapolation by phase-shift in the nonuniform Gabor domain Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureux ABSTRACT Wavefield extrapolation for
More informationMatrix formulation of adjoint Kirchhoff datuming
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 377 Matrix formulation of adjoint Kirchhoff datuming Dimitri Bevc 1 ABSTRACT I present the matrix formulation of Kirchhoff wave equation datuming
More informationA Padé approximation to the scalar wavefield extrapolator for inhomogeneous media
A Padé approimation A Padé approimation to the scalar wavefield etrapolator for inhomogeneous media Yanpeng Mi, Zhengsheng Yao, and Gary F. Margrave ABSTRACT A seismic wavefield at depth z can be obtained
More informationMain Menu SUMMARY INTRODUCTION
- a new method of solution Lasse Amundsen, Børge Arntsen, Arne Reitan, Eirik Ø Dischler, and Bjørn Ursin StatoilHydro Reseach Center, and NTNU SUMMARY We present a wave propagation method rigorous in one-way
More informationTOM 2.6. SEG/Houston 2005 Annual Meeting 2581
Oz Yilmaz* and Jie Zhang, GeoTomo LLC, Houston, Texas; and Yan Shixin, PetroChina, Beijing, China Summary PetroChina conducted a multichannel large-offset 2-D seismic survey in the Yumen Oil Field, Northwest
More informationOne-step extrapolation method for reverse time migration
GEOPHYSICS VOL. 74 NO. 4 JULY-AUGUST 2009; P. A29 A33 5 FIGS. 10.1190/1.3123476 One-step extrapolation method for reverse time migration Yu Zhang 1 and Guanquan Zhang 2 ABSTRACT We have proposed a new
More informationFourier finite-difference wave propagation a
Fourier finite-difference wave propagation a a Published in Geophysics, 76, T123-T129, (2011) Xiaolei Song and Sergey Fomel ABSTRACT We introduce a novel technique for seismic wave extrapolation in time.
More informationAn acoustic wave equation for orthorhombic anisotropy
Stanford Exploration Project, Report 98, August 1, 1998, pages 6?? An acoustic wave equation for orthorhombic anisotropy Tariq Alkhalifah 1 keywords: Anisotropy, finite difference, modeling ABSTRACT Using
More informationMAS. Summary. (Popovici. 1996), and semi-recursive Kirchoff (Bevc, 1997).
Prestack depth migration for complex 2D structure using phase-screen propagators 9 Peter Roberts *, Lian-Jie Huang, Los AIamos National Laboratory, Charles Burch. Conoco Inc.. Michael ehlt ri i d Ste1,en
More informationSUMMARY INTRODUCTION THEORY
Least-squares reverse-time migration using one-step two-way wave extrapolation by non-stationary phase shift Junzhe Sun, Sergey Fomel, and Jingwei Hu, The University of Texas at Austin Downloaded 11/06/14
More informationComplex-beam Migration and Land Depth Tianfei Zhu CGGVeritas, Calgary, Alberta, Canada
Page 1 of 10 Home Articles Interviews Print Editions Complex-beam Migration and Land Depth Tianfei Zhu CGGVeritas, Calgary, Alberta, Canada DECEMBER 2012 FOCUS ARTICLE Summary Gaussian-beam depth migration
More informationCompensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS
Compensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS SUMMARY Anelastic properties of the earth cause frequency dependent energy
More informationTraveltime sensitivity kernels: Banana-doughnuts or just plain bananas? a
Traveltime sensitivity kernels: Banana-doughnuts or just plain bananas? a a Published in SEP report, 103, 61-68 (2000) James Rickett 1 INTRODUCTION Estimating an accurate velocity function is one of the
More informationStolt residual migration for converted waves
Stanford Exploration Project, Report 11, September 18, 1, pages 1 6 Stolt residual migration for converted waves Daniel Rosales, Paul Sava, and Biondo Biondi 1 ABSTRACT P S velocity analysis is a new frontier
More informationAcoustic anisotropic wavefields through perturbation theory
GEOPHYSICS, VOL. 78, NO. 5 (SEPTEMBER-OCTOBER 2013); P. WC41 WC50, 22 FIGS. 10.1190/GEO2012-0391.1 Acoustic anisotropic wavefields through perturbation theory Tariq Alkhalifah 1 ABSTRACT Solving the anisotropic
More informationAmplitude preserving AMO from true amplitude DMO and inverse DMO
Stanford Exploration Project, Report 84, May 9, 001, pages 1 167 Amplitude preserving AMO from true amplitude DMO and inverse DMO Nizar Chemingui and Biondo Biondi 1 ABSTRACT Starting from the definition
More informationStanford Exploration Project, Report 105, September 5, 2000, pages
Stanford Exploration Project, Report 5, September 5, 2, pages 57 66 56 Stanford Exploration Project, Report 5, September 5, 2, pages 57 66 Short Note Simple factorization of implicit finite-difference
More informationG022 Multi-azimuth Seismic Data Imaging in the Presence of Orthorhombic Anisotropy
G0 Multi-azimuth Seismic Data Imaging in the Presence of Orthorhomic Anisotropy Y. Xie* (CGGVeritas), S. Birdus (CGGVeritas), J. Sun (CGGVeritas) & C. Notfors (CGGVeritas) SUMMARY The presence of orthorhomic
More informationP137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients
P137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients M. Ayzenberg (StatoilHydro), A. Aizenberg (Institute of Petroleum Geology and Geophysics),
More informationJ.A. Haugen* (StatoilHydro ASA), J. Mispel (StatoilHydro ASA) & B. Arntsen (NTNU)
U008 Seismic Imaging Below "Dirty" Salt J.A. Haugen* (StatoilHydro ASA), J. Mispel (StatoilHydro ASA) & B. Arntsen (NTNU) SUMMARY Base and sub salt seismic imaging is still an unresolved issue. To solve
More informationLowrank finite-differences and lowrank Fourier finite-differences for seismic wave extrapolation in the acoustic approximation a
Lowrank finite-differences and lowrank Fourier finite-differences for seismic wave extrapolation in the acoustic approximation a a Published in Geophysical Journal International, 193, 960-969, (2013) Xiaolei
More informationRobert W. Vestrum and R. James Brown
Group versus phase velocity in laboratory measurements on orthorhombic samples Robert W. Vestrum and R. James Brown ABSTRACT In laboratory measurements of traveltimes across anisotropic materials, there
More informationStanford Exploration Project, Report 115, May 22, 2004, pages
Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 248 Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 First-order lateral interval velocity estimates without picking
More informationSubsalt imaging by common-azimuth migration
Stanford Exploration Project, Report 100, April 20, 1999, pages 113 125 Subsalt imaging by common-azimuth migration Biondo Biondi 1 keywords: migration, common-azimuth, wave-equation ABSTRACT The comparison
More informationDirect nonlinear traveltime inversion in layered VTI media Paul J. Fowler*, Alexander Jackson, Joseph Gaffney, and David Boreham, WesternGeco
Paul J. Fowler*, Alexander Jackson, Joseph Gaffney, and David Boreham, WesternGeco Summary We present a scheme for direct nonlinear inversion of picked moveout traveltimes in block-layered isotropic or
More informationWave-equation tomography for anisotropic parameters
Wave-equation tomography for anisotropic parameters Yunyue (Elita) Li and Biondo Biondi ABSTRACT Anisotropic models are recognized as more realistic representations of the subsurface where complex geological
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (04) 98 653 66 GJI Seismology doi: 0.093/gji/ggu9 Effective orthorhombic anisotropic models for wavefield extrapolation Wilson Ibanez-Jacome Tariq Alkhalifah
More informationOn non-stationary convolution and inverse convolution
Stanford Exploration Project, Report 102, October 25, 1999, pages 1 137 On non-stationary convolution and inverse convolution James Rickett 1 keywords: helix, linear filtering, non-stationary deconvolution
More informationWave field extrapolation using a new multiplication/convolution method. work in progress
Wave field extrapolation using a new multiplication/convolution method work in progress Christiaan C. Stolk ABSTRACT In this paper we consider numerical methods for the acoustic equation, based on evolution
More informationThe Application of Depth Migration for Processing GPR Data
The Application of Depth Migration for Processing GPR Data Dang Hoai Trung 1, *, Nguyen Van Giang, and Nguyen Thanh Van 1 1 Faculty of Physics and Engineering Physics, VNUHCM University of Science Vietnam
More informationSUMMARY INTRODUCTION THEORY
Seismic modeling using the frozen Gaussian approximation Xu Yang, University of California Santa Barbara, Jianfeng Lu, Duke University, and Sergey Fomel, University of Texas at Austin SUMMARY We adopt
More informationPrestack exploding reflector modelling and migration for anisotropic media
Geophysical Prospecting, 015, 63, 10 doi: 10.1111/1365-478.1148 Prestack exploding reflector modelling and migration for anisotropic media T. Alkhalifah Physical Sciences and Engineering, King Abdullah
More informationExamples of prestack depth migration in TI media
Eamples of prestack depth migration in TI media Robert J. Ferguson and Gary F. Margrave ABSTRACT Wave field etrapolation by nonstationary phase shift can be formulated to allow velocity variation with
More informationStanford Exploration Project, Report 123, October 31, 2005, pages 25 49
Stanford Exploration Project, Report 123, October 31, 2005, pages 25 49 24 Stanford Exploration Project, Report 123, October 31, 2005, pages 25 49 Residual moveout in anisotropic angle-domain common image
More informationExact wavefield extrapolation in 2D for v(x)
Exact extrapolation Exact wavefield extrapolation in D for v(x) Gary F. Margrave, Michael P. Lamoureux, Peter Gibson, Richard A. Bale, and Jeff Grossman ABSTRACT We pose the wavefield extrapolation problem
More informationSolving physics pde s using Gabor multipliers
Solving physics pde s using Michael P. Lamoureux, Gary F. Margrave, Safa Ismail ABSTRACT We develop the mathematical properties of, which are a nonstationary version of Fourier multipliers. Some difficulties
More informationElastic wavefield separation for VTI media
CWP-598 Elastic wavefield separation for VTI media Jia Yan and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT The separation of wave modes from isotropic elastic wavefields is typically
More informationStanford Exploration Project, Report SERGEY, November 9, 2000, pages 683??
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 683?? 682 Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 683?? Velocity continuation by spectral methods Sergey
More informationLowrank RTM for converted wave imaging
Lowrank RM for converted wave imaging Lorenzo Casasanta * (CGG), Zhiguang Xue (U Austin), Sam Gray (CGG) Summary his paper is an attempt to fill the technology gap existing between pure P- and PS-wave
More informationExtended isochron rays in prestack depth (map) migration
Extended isochron rays in prestack depth (map) migration A.A. Duchkov and M.V. de Hoop Purdue University, 150 N.University st., West Lafayette, IN, 47907 e-mail: aduchkov@purdue.edu (December 15, 2008)
More informationMulticomponent imaging with distributed acoustic sensing Ivan Lim Chen Ning & Paul Sava, Center for Wave Phenomena, Colorado School of Mines
Multicomponent imaging with distributed acoustic sensing Ivan Lim Chen Ning & Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY The usage of Distributed Acoustic Sensing (DAS in geophysics
More informationAzimuth moveout for 3-D prestack imaging
Azimuth moveout for 3-D prestack imaging Biondo Biondi, Sergey Fomel, and Nizar Chemingui Stanford Exploration Project, 360 Mitchell Bldg., Stanford, CA 94305-15 (July 30, 1997) ABSTRACT We introduce a
More informationNonstationary filters, pseudodifferential operators, and their inverses
Nonstationary filters, pseudodifferential operators, and their inverses Gary F. Margrave and Robert J. Ferguson ABSTRACT An inversion scheme for nonstationary filters is presented and explored. Nonstationary
More informationSTANFORD EXPLORATION PROJECT
STANFORD EXPLORATION PROJECT Gabriel Alvarez, Brad Artman, James Berryman, Biondo Biondi, Robert Clapp, Antoine Guitton, Jeff Shragge, Alejandro Valenciano, and Huazhong Wang Report Number 123, October
More informationDownloaded 08/20/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Viscoacoustic modeling and imaging using low-rank approximation Junzhe Sun 1, Tieyuan Zhu 2, and Sergey Fomel 1, 1 The University of Texas at Austin; 2 Stanford University SUMMARY A constant-q wave equation
More informationAn efficient wave extrapolation method for tilted orthorhombic media using effective ellipsoidal models
An efficient wave extrapolation method for tilted orthorhombic media using effective ellipsoidal models Item Type Conference Paper Authors Waheed Umair bin; Alkhalifah Tariq Ali Eprint version Pre-print
More informationEstimation of fracture parameters from reflection seismic data Part II: Fractured models with orthorhombic symmetry
GEOPHYSICS, VOL. 65, NO. 6 (NOVEMBER-DECEMBER 2000); P. 1803 1817, 5 FIGS., 1 TABLE. Estimation of fracture parameters from reflection seismic data Part II: Fractured models with orthorhomic symmetry Andrey
More informationElastic wave-equation migration for laterally varying isotropic and HTI media. Richard A. Bale and Gary F. Margrave
Elastic wave-equation migration for laterally varying isotropic and HTI media Richard A. Bale and Gary F. Margrave a Outline Introduction Theory Elastic wavefield extrapolation Extension to laterally heterogeneous
More informationThe i-stats: An Image-Based Effective-Medium Modeling of Near-Surface Anomalies Oz Yilmaz*, GeoTomo LLC, Houston, TX
The i-stats: An Image-Based Effective-Medium Modeling of Near-Surface Anomalies Oz Yilmaz*, GeoTomo LLC, Houston, TX Summary Near-surface modeling for statics corrections is an integral part of a land
More informationSection 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)).
Difference Equations to Differential Equations Section 8.5 Applications: Pendulums Mass-Spring Systems In this section we will investigate two applications of our work in Section 8.4. First, we will consider
More informationFull waveform inversion in the Laplace and Laplace-Fourier domains
Full waveform inversion in the Laplace and Laplace-Fourier domains Changsoo Shin, Wansoo Ha, Wookeen Chung, and Ho Seuk Bae Summary We present a review of Laplace and Laplace-Fourier domain waveform inversion.
More information3D VTI traveltime tomography for near-surface imaging Lina Zhang*, Jie Zhang, Wei Zhang, University of Science and Technology of China (USTC)
Downloaded 01/03/14 to 16.01.198.34. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ 3D VTI traveltime tomography for near-surface imaging Lina Zhang*, Jie
More informationEdinburgh Research Explorer
Edinburgh Research Explorer Nonlinear Scattering-based Imaging in Elastic Media Citation for published version: Ravasi, M & Curtis, A 213, 'Nonlinear Scattering-based Imaging in Elastic Media' Paper presented
More informationImage-space wave-equation tomography in the generalized source domain
Image-space wave-equation tomography in the generalized source domain Yaxun Tang, Claudio Guerra, and Biondo Biondi ABSTRACT We extend the theory of image-space wave-equation tomography to the generalized
More informationFinite difference modelling of the full acoustic wave equation in Matlab
Acoustic finite difference modelling Finite difference modelling of the full acoustic wave equation in Matlab Hugh D Geiger and Pat F Daley ABSTRACT Two subroutines have been added to the Matlab AFD (acoustic
More informationTOM 1.7. Sparse Norm Reflection Tomography for Handling Velocity Ambiguities
SEG/Houston 2005 Annual Meeting 2554 Yonadav Sudman, Paradigm and Dan Kosloff, Tel-Aviv University and Paradigm Summary Reflection seismology with the normal range of offsets encountered in seismic surveys
More informationElastic wave-mode separation for VTI media
GEOPHYSICS, VOL. 74, NO. 5 SEPTEMBER-OCTOBER 2009 ; P. WB19 WB32, 19 FIGS. 10.1190/1.3184014 Elastic wave-mode separation for VTI media Jia Yan 1 and Paul Sava 1 ABSTRACT Elastic wave propagation in anisotropic
More informationMultidimensional partitions of unity and Gaussian terrains
and Gaussian terrains Richard A. Bale, Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureu ABSTRACT Partitions of unity play an important rôle as amplitude-preserving windows for nonstationary
More informationTime-migration velocity analysis by velocity continuation
GEOPHYSICS, VOL. 68, NO. 5 (SEPTEMBER-OCTOBER 2003); P. 1662 1672, 15 FIGS. 10.1190/1.1620640 Time-migration velocity analysis by velocity continuation Sergey Fomel ABSTRACT Time-migration velocity analysis
More informationv v Downloaded 01/11/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
The pseudo-analytical method: application of pseudo-laplacians to acoustic and acoustic anisotropic wae propagation John T. Etgen* and Serre Brandsberg-Dahl Summary We generalize the pseudo-spectral method
More informationAttenuation compensation in viscoacoustic reserve-time migration Jianyong Bai*, Guoquan Chen, David Yingst, and Jacques Leveille, ION Geophysical
Attenuation compensation in viscoacoustic reserve-time migration Jianyong Bai*, Guoquan Chen, David Yingst, and Jacques Leveille, ION Geophysical Summary Seismic waves are attenuated during propagation.
More informationP125 Method for Calibrating Seismic Imaging Velocities
P125 Method for Calibrating Seismic Imaging Velocities P.C. Docherty* (Fairfield Industries Inc.) SUMMARY Anisotropy can cause large depth errors in images obtained with isotropic velocity analysis and
More informationFast Kirchhoff migration in the wavelet domain
Exploration Geophysics (2002) 33, 23-27 Fast Kirchhoff migration in the wavelet domain Valery A. Zheludev 1 Eugene Ragoza 2 Dan D. Kosloff 3 Key Words: Kirchhoff migration, wavelet, anti-aliasing ABSTRACT
More informationStanford Exploration Project, Report 97, July 8, 1998, pages
Stanford Exploration Project, Report 97, July 8, 998, pages 7 6 Stanford Exploration Project, Report 97, July 8, 998, pages 7 The offset-midpoint traveltime pyramid in transversely isotropic media Tariq
More informationSUMMARY INTRODUCTION LOCALIZED PHASE ESTIMATION
Local similarity with the envelope as a seismic phase detector Sergey Fomel, The University of Texas at Austin, and Mirko van der Baan, University of Alberta SUMMARY We propose a new seismic attribute,
More informationApproximate Methods for Time-Reversal Processing of Large Seismic Reflection Data Sets. University of California, LLNL
147th Annual Meeting: Acoustical Society of America New York City, May 25, 2004 Approximate Methods for Time-Reversal Processing of Large Seismic Reflection Data Sets Speaker: James G. Berryman University
More informationA modeling and migration study of fault shadows: A synthetic example
Fault shadow imaging A modeling and migration study of fault shadows: A synthetic example Sitamai W. Ajiduah and Gary F. Margrave This report presents a study that uses 2D finite difference modeling and
More informationStress induced seismic velocity anisotropy study Bin Qiu, Mita Sengupta, Jorg Herwanger, WesternGeco/Schlumberger
anisotropy study Bin Qiu, Mita Sengupta, Jorg Herwanger, WesternGeco/Schlumberger Summary A method of calculating stress induced seismic velocity anisotropy is presented and studied. In a first step an
More informationShot-profile migration of multiple reflections
Stanford Exploration Project, Report 111, June 10, 2002, pages?? 33 Shot-profile migration of multiple reflections Antoine Guitton 1 ABSTRACT A shot-profile migration algorithm is modified to image multiple
More informationDeconvolution imaging condition for reverse-time migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 83 96 Deconvolution imaging condition for reverse-time migration Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT The reverse-time
More informationSummary. 3D potential field migration. Various approximations of the Newton method result in the imaging condition for the vector of densities, :
24 A case study from Broken Hill, Australia, with comparison to 3D regularized inversion Michael S. Zhdanov, University of Utah and TechnoImaging, Xiaojun Liu, University of Utah, Le Wan*, Martin Čuma,
More informationDownloaded 05/01/17 to Redistribution subject to SEG license or copyright; see Terms of Use at
Mapping Imbricate Structures in the Thrust Belt of Southeast Turkey by Large-Offset Seismic Survey Oz Yilmaz*, Anatolian Geophysical, Istanbul, Turkey; and Serdar Uygun, Ali Ölmez, and Emel Çalı, Turkish
More informationVelocity analysis using AB semblance a
Velocity analysis using AB semblance a a Published in Geophysical Prospecting, v. 57, 311-321 (2009) Sergey Fomel ABSTRACT I derive and analyze an explicit formula for a generalized semblance attribute,
More informationExploring three-dimensional implicit wavefield extrapolation with the helix transform a
Exploring three-dimensional implicit wavefield extrapolation with the helix transform a a Published in SEP report, 95, 43-60 (1997) Sergey Fomel and Jon F. Claerbout 1 ABSTRACT Implicit extrapolation is
More informationAVAZ inversion for fracture orientation and intensity: a physical modeling study
AVAZ inversion for fracture orientation and intensity: a physical modeling study Faranak Mahmoudian*, Gary F. Margrave, and Joe Wong, University of Calgary. CREWES fmahmoud@ucalgary.ca Summary We present
More informationAn Improved Dual Sensor Summation Method with Application to Four-Component (4-C) Seafloor Seismic Data from the Niger Delta
Earth Science Research; Vol. 4, No. 2; 2015 ISSN 1927-0542 E-ISSN 1927-0550 Published by Canadian Center of Science and Education An Improved Dual Sensor Summation Method with Application to Four-Component
More informationBureau of Economic Geology, The University of Texas at Austin, Austin, TX. Research Associate, Geophysics from 11/2004
Paul C. Sava Bureau of Economic Geology, The University of Texas at Austin, University Station, Box X, Austin, TX 78758, (512) 471-0327 paul.sava@beg.utexas.edu, http://sepwww.stanford.edu/sep/paul Research
More informationTh Guided Waves - Inversion and Attenuation
Th-01-08 Guided Waves - Inversion and Attenuation D. Boiero* (WesternGeco), C. Strobbia (WesternGeco), L. Velasco (WesternGeco) & P. Vermeer (WesternGeco) SUMMARY Guided waves contain significant information
More informationA phase-space formulation for elastic-wave traveltime tomography
A phase-space formulation for elastic-wave traveltime tomography Eric Chung 1, Jianliang Qian 2, Gunther Uhlmann 3, Hong-Kai Zhao 4 1 Applied and Computational Mathematics, California Institute of Technology,
More informationPseudo-acoustic wavefield propagation for anisotropic media
Anisotropic wavefield propagation Pseudo-acoustic wavefield propagation for anisotropic media Ben D. Wards, Gary F. Margrave, and Michael P. Lamoureux ABSTRACT Reverse-time migration (RTM) is a powerful
More informationNonlinear seismic imaging via reduced order model backprojection
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, Vladimir Druskin 2 and Mikhail Zaslavsky 2 University of Houston, 2 Schlumberger-Doll Research Center Mamonov, Druskin,
More informationVelocity Update Using High Resolution Tomography in Santos Basin, Brazil Lingli Hu and Jianhang Zhou, CGGVeritas
Lingli Hu and Jianhang Zhou, CGGVeritas Summary The exploration interest in the Santos Basin offshore Brazil has increased with the large deep water pre-salt discoveries, such as Tupi and Jupiter. As the
More informationTraveltime approximations for inhomogeneous transversely isotropic media with a horizontal symmetry axis
Geophysical Prospecting, 2013, 61, 495 503 doi: 10.1111/j.1365-2478.2012.01067.x Traveltime approximations for inhomogeneous transversely isotropic media with a horizontal symmetry axis Tariq Alkhalifah
More informationWave-equation prestack depth migration for sub-basalt P and converted wave imaging
Stanford Exploration Project, Report 110, September 18, 2001, pages 1 31 Wave-equation prestack depth migration for sub-basalt P and converted wave imaging Morgan Brown, Biondo Biondi, and Clement Kostov
More informationConstrained inversion of P-S seismic data
PS Inversion Constrained inversion of P-S seismic data Robert J. Ferguson, and Robert R. Stewart ABSTRACT A method to estimate S-wave interval velocity, using P-S seismic data is presented. The method
More informationStatic Corrections for Seismic Reflection Surveys
Static Corrections for Seismic Reflection Surveys MIKE COX Volume Editors: Series Editor: Eugene F. Scherrer Roland Chen Eugene F. Scherrer Society of Exploration Geophysicists Tulsa, Oklahoma Contents
More informationInvestigating fault shadows in a normally faulted geology
Investigating fault shadows in a normally faulted geology Sitamai Ajiduah* and Gary Margrave CREWES, University of Calgary, sajiduah@ucalgary.ca Summary Fault shadow poses a potential development risk
More informationRobot Position from Wheel Odometry
Root Position from Wheel Odometry Christopher Marshall 26 Fe 2008 Astract This document develops equations of motion for root position as a function of the distance traveled y each wheel as a function
More informationAnisotropic tomography for TTI and VTI media Yang He* and Jun Cai, TGS
media Yang He* and Jun Cai, TGS Summary A simultaneous anisotropic tomographic inversion algorithm is developed. Check shot constraints and appropriate algorithm preconditioning play an important role
More information