Least. Least squares migration of converted wave seismic data. Introduction FOCUS ARTICLE. Aaron Stanton and Mauricio D. Sacchi
|
|
- Jesse Wilkins
- 5 years ago
- Views:
Transcription
1 Coordinated by Rob Holt Least squares migration of converted wave seismic data Aaron Stanton and Mauricio D. Sacchi DEPARTMENT OF PHYSICS, UNIVERSITY OF ALBERTA, EDMONTON, ALBERTA, CANADA Least squares migration permits the estimation of subsurface models that honor the recorded primary seismic wavefield. In addition, it permits us to include regularization constraints that reduce sampling and illumination artifacts. Least squares migration utilizing one-way wave equation migration has been mainly restricted to imaging via scalar wave propagation operators. We present an extension of wave equation least squares migration for vector wave fields. Synthetic results demonstrate that converted wave least squares migration can improve the resolution of PP and PS images, and can also be used for interpolation and wavefield separation. Introduction Converted waves are a useful tool for characterizing the elastic properties of a reservoir for a number of reasons. The comparison of PP and PS reflectivities can help us identify fluid contacts, and S-wave ray paths provide a tool to illuminate targets hidden beneath gas filled zones. Perhaps most importantly, P and S waves provide complementary information that is helpful when inverting for elastic parameters (Stewart, 1990). Prestack migration is a critical step in the processing of converted waves. For P-wave seismology it is common to approximate wave propagation using the scalar wave equation. This approximation is often sufficient to produce reasonable P-wave structural images. Extending this approximation to converted wave imaging requires some troubling assumptions about amplitude and polarization. Migration algorithms based on the elastic wave equation more accurately model wave propagation and allow multiple wavefields to be imaged simultaneously. There are a variety of approaches to elastic migration. Kuo and Dai (1984) developed an elastic implementation of Kirchhoff migration. Chang and McMechan (1987) applied Reverse Time Migration (RTM) to multicomponent data using an elastic finite difference algorithm, resulting in horizontal and vertical component images. Dellinger and Etgen (1990) proposed the application of Helmholtz decomposition to elastic data using a Fourier domain operator. Later implementations of elastic RTM use this idea to provide distinct PP and PS images (see for example Yan and Sava, 2008). Extensive work has also been done on wave equation based methods to migrate elastic data. Xie and Wu (2005) used an elastic version of split step migration to downward continue elastic data, while Bale (2006) modified the phase shift migration to handle anisotropic elastic wavefields. Wave equation based methods are an attractive option because they perform more accurately than ray-based methods, and are much less computationally expensive than methods that employ finite differences. For example, in time domain RTM many wavefield snapshots for all time samples must be generated and saved prior to imaging, whilst in wave equation migration individual frequency slices are treated independently and the wavefield is recursively updated with depth. The migration of converted waves comes with some additional complications when compared with P-wave migration. Typically PS data have a lower signal to noise ratio than PP data, and a slower S-leg gives converted waves an asymmetric ray path. While this implies that converted wave data have more restrictive aliasing criteria, it also presents an opportunity: given the same acquisition geometry converted waves can image the subsurface with a wider source side aperture than P-wave data. Furthermore, ignoring the effects of attenuation, a given frequency of converted wave data images the earth with a higher resolution than P-wave data. Realistically, attenuation is always a factor (especially in the near surface), which challenges the resolution of both PP and PS depth images (Bale and Stewart, 2002). Some of these challenges can be dealt with before migration. For example, regularization can be applied to reduce the effects of acquisition footprint on the migrated image (Cary, 2011). In this article we aim to address some of the challenges facing converted wave migration via least squares inversion. Least squares migration, in its various forms, has been an active field of research for many years. Lambaré et al. (1992) used an iterative data fitting approach to solve for reflectivity. Nemeth et al. (1999) applied least squares migration with a Kirchhoff operator to image in the presence of poor spatial sampling, while Kuehl and Sacchi (2003) used regularized least squares wave equation migration to generate angle gathers that are unaffected by acquisition footprint. Wang (2005) and Wang et al. (2005) explore the application of 3D wave equation least squares migration with constraints to a real data set of the Western Canadian Sedimentary Basin. In a subsequent contribution, Wang and Sacchi (2007) also explore the incorporation of a sparsity constraint to least square migration to the ubiquitous problem of vertical resolution enhancement. Least squares migration of structurally complex data has been also being investigated via the utilization of preconditioning 48
2 operators that were synthetized with prediction error filters (Wang and Sacchi, 2009). A common theme of all least squares migration algorithms is data fitting which depends greatly on the ability of the migration operator to accurately propagate energy from source to receiver. Extending least squares migration from the acoustic to the elastic case is a natural progression in data fitting. Recently, Stanton (2014) proposed using Least squares prestack time migration to image and regularize converted waves. In this article, we tackle the more complex problem of utilizing Least squares migration with wave equation prestack depth imaging operators. We implement least squares migration using an elastic shot profile wave equation migration operator with an angle gather imaging condition. We use a synthetic data example to demonstrate that the algorithm can improve the resolution of PP and PS images. The algorithm can also be used for interpolation and wavefield separation. Methodology Least squares migration aims to find a migrated image that minimizes the difference between observed and predicted data. The cost function for least squares migration of acoustic data with quadratic regularization is J = TTTTTT dd + µμ mm (1) where T is a sampling operator, L is a forward (demigration) operator, m is the migrated image we wish to solve for, and d are the observed data. We minimize equation 1 using Conjugate Gradients. Each iteration of CG involves applying the forward operator to predict data from the migrated image, computing a data residual, and migrating the residual using the adjoint operator to compute an update to the image. Before adapting equation 1 to the elastic case we must first define our elastic migration operator. We form an elastic migration operator by combining Helmholtz decomposition with a wave equation migration operator that extrapolates P and S scalar potentials independently. After wavefield separation we extrapolate P and S wave scalar potentials using a shot profile wave equation migration operator with an angle gather imaging condition. To generate images m pp and m ps simultaneously, we can write LL = LL 0 0 LL " where L T pp and L T ps implement split step wave equation migration (Stoffa et al., 1990) of PP and PS images respectively. Split step migration consists of downward continuing the wavefield by depth increments using a phase shift determined using a background velocity, and performing a correction for the perturbation of the background velocity at each lateral position. The source is extrapolated using the P wave velocity. On the receiver side, operator L T pp extrapolates using the P wave velocity, while operator L T ps extrapolates using the S-wave velocity. Typically in wave equation migration the source and receiver are moved to a coincident position in the subsurface (migration to zero offset). To do this we extrapolate source and receiver wavefields and cross correlate them once for each lateral position we wish to image. After migration, all shots are summed to form the image. Forming angle gathers requires a small adjustment to this procedure. Instead of simply cross correlating the extrapolated source and receiver wavefields once at each lateral position, the cross correlation is repeated several times for different spatial lags. These lags can be interpreted as the subsurface offset between the two wavefields. Individual shots are stacked after migration and the offset axis is converted to angle using the relation tan γγ = (5) where γ is the half opening angle, and and are are the wavenumbers for the offset and the depth axes respectively (Rickett and Sava, (4) Isotropic elastic data can be decomposed into P and S-wave potentials by taking the divergence and curl of the wavefield components respectively: P = +, S =. (2) This is commonly referred to as Helmholtz decomposition. Etgen (1988) evaluated these derivatives by Fourier transform. Writing this as a linear operator: ik ik (v ) HH = FF FF ik (v ) ik (3) where F x and F x 1 are the forward and inverse Fourier transform on the spatial axis respectively. In this equation the lateral wavenumber is given by k x, and the vertical wavenumbers k z (v p ) and k z (v s ) depend on the P and S wave velocities of the near surface. Technically this upper layer should have a constant velocity, but practically we find that the average velocity is generally a good approximation. Figure 1. P and S wave velocities for a normal fault. The S wave velocities were computed as vv 3. Continued on Page 50 DECEMBER 2014 CSEG RECORDER 49
3 Continued from Page ). For a more detailed derivation of the forward and adjoint operators in least squares shot profile migration, we refer the interested reader to Kaplan et al. (2010). Using the wavefield separation and migration operators defined above, we can combine them to form the forward and adjoint operators needed for inversion. To forward model two component data we use dd = HH LL dd 0 mm 0 LL " mm (6) " 3 shows the result of migration for a constant opening angle of -42. Because of the missing data the wavefield separation did not perform well and some cross-talk persists in both images. Figure 4 shows the result of 5 iterations of least squares migration. Clearly the wavefield separation has been improved as the cross-talk artifacts have been reduced. Another improvement comes in the form of resolution. Figure 5 shows the amplitude spectrum for the migrated (blue) and least squares migrated data (red). For both m pp and m ps the amplitude have Conversely, to migrate two component data we take the adjoint of the forward operator to get mm = LL 0 mm " 0 LL " HH dd dd. (7) Finally, for elastic least squares migration with quadratic regularization the cost function is written J = TTHH LL 0 0 LL " mm mm " dd dd + µμ mm mm " (8) Again, this cost function can be minimized using CG. To improve convergence we can add preconditioning to suppress crosstalk artifacts in the angle domain using an FK fan filter. To demonstrate our algorithm we created a model consisting of bedding transected by a normal fault. The P and S wave velocities for this model are shown in Figure 1. We generated 115 two component shot gathers using elastic finite difference modeling using a pressure source with a 40Hz peak frequency Ricker wavelet. A representative shot gather is shown in Figure 2. Notice the mixture of PP and PS wave events between components. We decimated 50% of the receivers randomly to simulate a more realistic acquisition scenario. Under perfect sampling our Helmholtz decomposition operator, H, performs well in decomposing the vector wavefield to P and S potentials. Figure Figure 3. Migrated PP and PS images for a constant opening angle of -42. Notice some cross-talk artifacts on both sections as a result of mperfect wavefield separation. Figure 2. X and Z components for a common shot gather generated using elastic finite difference modeling with the model shown in Figure 1. 50% of the data have been decimated to simulate irregular receiver acquisition geometry. Figure 4. Least squares migrated PP and PS images for a constant opening angle of -42. Notice that many of the cross-talk artifacts have been reduced and the resolution has been improved. 50
4 Figure 7. Predicted fully sampled P waves on the X and Z components of one shot gather. Figure 5. Amplitude spectra for the PP (top) and PS (bottom) depth sections. The migrated amplitude spectra are plotted in blue, while the least squares migrated amplitude spectra are plotted in red. Least squares migration improves the bandwidth of the migrated image. been broadened for low and high frequencies. Figure 6 shows the misfit for 5 iterations of CG. The algorithm required just a few iterations to converge to a satisfactory result. Figure 7 demonstrates two interesting applications of elastic least squares migration: interpolation and wavefield separation. After finding an acceptable image that best fits the data, the forward operator can be applied to predict data: "#$ dd "#$ dd = LL HH 0 mm 0 LL " mm (9) " This prediction can be used for interpolation and denoising of the data. For wavefield separation one of the images m pp or m ps can be set to zero allowing for the prediction of P waves on both components (or vice versa). In Figure 7 we see the result of predicting P waves at all trace locations for one shot. This approach can be used to predict entirely new source locations, or as shown in this case, to infill missing receivers for a previously existing source. Figure 6. Data misfit for 5 iterations of the inversion. The misfit is computed as the sum of the squared differences between observed and predicted data for both data components only at observed trace locations. Conclusions Least squares migration algorithms attempt to fit recorded data with predictions generated from a migrated image. By improving the accuracy of the migration operator to include elastic wave propagation we expect to improve the ability of least squares migration to fit seismic amplitudes. Our example demonstrates that least squares migration can improve the resolution of PP and PS images. The algorithm can also be used for data reconstruction and wavefield separation. Continued on Page 52 DECEMBER 2014 CSEG RECORDER 51
5 Continued from Page 51 Acknowledgements We thank the sponsors of the Signal Analysis and Imaging Group for their continued support, and Westgrid for computational resources. The figures in this article were generated using the open source software Madagascar. References Bale, R., 2006, Elastic wave-equation depth migration of seismic data for isotropic and azimuthally anisotropic media: PhD thesis, University of Calgary. Bale, R., and Stewart, R., 2002, The impact of attenuation on the resolution of multicomponent seismic data: CREWES research report. Cary, P., 2011, True-Amplitude PS Prestack Time Migration via 5D Interpolation: Presented at the CSEG conference and exhibition. Chang, W.F., and McMechan, G.A., 1987, Elastic reverse-time migration: Geophysics, 52, Dellinger, J., and Etgen, J., 1990, Wavefield separation in two-dimensional anisotropic media: Geophysics, 55, Etgen, J.T., 1988, Elastic prestack migration of two component data: Stanford Exploration Project, 57, Kaplan, S., Routh, P., and Sacchi, M.D., 2010, Derivation of forward and adjoint operators for least-squares shot-profile split-step migration: Geophysics, 75(6), S225-S235. Kuehl, H., and Sacchi, M.D., 2003, Least-squares wave-equation migration for AVP/ AVA inversion: Geophysics, 68, Kuo, J. T., and Dai, T., 1984, Kirchhoff elastic wave migration for the case of noncoincident source and receiver: Geophysics, 49, Lambaré, G., Virieux, J., Madariaga, R., and Jin, S., 1992, Iterative asymptotic inversion in the acoustic approximation: Geophysics, 57, no. 9, , Nemeth, T., Wu, C., and Schuster, G.T., 1999, Least-squares migration of incomplete data: Geophysics, 64, Rickett, J. E., and Sava, P.C., 2002, Offset and angle-domain common image-point gathers for shot-profile migration: Geophysics, 67, Stanton, A., 2014, Least-squares migration of converted waves, GeoConvention, expanded abstracts, 6 pages. Stewart, R., 1990 Joint P and P-SV inversion: CREWES research report. Stoffa, P.L., Fokkema, J.T., de Luna Freire, R.M., and Kessinger, W.P., 1990: Splitstep Fourier migration, Geophysics 55, Wang, J., 2005, 3-D Least-squares Wave-equation AVP/AVA Migration of common Azimuth data: PhD Thesis, Department of Physics, University of Alberta. Wang, J., Kuehl, H., and Sacchi, M.D., 2005, High-resolution wave-equation AVA imaging: Algorithm and tests with a data set from the Western Canadian Sedimentary Basin: Geophysics, 70 (5), Wang, J., and Sacchi, M.D., 2007, High-resolution wave equation AVP imaging with sparseness constraints: Geophysics, 72 (1), S11-S18. Wang, J., and Sacchi, M.D., 2009, Structure constrained least-squares migration: SEG Expanded Abstracts, 28, Xie, X.B., and Wu, R.S., 2005, Multicomponent prestack depth migration using elastic screen method: Geophysics, 70, S30-S37. Yan, J., and Sava, P., 2008, Isotropic angle-domain elastic reverse-time migration: Geophysics, 73, S229 S239. Aaron Stanton completed a B.Sc. in Geophysics from the University of Alberta in From 2007 to 2011 he worked as a processing project leader for CGG. He is now pursuing a Ph.D. in Geophysics at the University of Alberta where his Ph.D. work is focused on the regularization and imaging of multicomponent seismic data. He can be reached at aaron.stanton@ualberta.ca. Mauricio Sacchi received a diploma in geophysics from The National University of La Plata, Argentina, in 1988 and a Ph.D. in geophysics from UBC, Canada, in He joined the Department of Physics at the University of Alberta (Edmonton, Canada) in His research interests are in the area of signal analysis and imaging methods. He has developed and taught short courses for the industry and for SEG, CSEG, and EAGE in seismic signal theory, transform methods for signal enhancement, seismic inversion, and multidimensional data reconstruction. With Tad Ulrych, he wrote the book Information-based processing and inversion with applications (Elsevier). He is the recipient of the 2012 CSEG Medal, and is the 2014 SEG Honorary Speaker for Latin America and the Caribbean. 52 DECEMBER 2014 CSEG RECORDER
Robustness of the scalar elastic imaging condition for converted waves
CWP-830 Robustness of the scalar elastic imaging condition for converted waves Yuting Duan & Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT For elastic reverse-time migration, one
More informationSUMMARY ELASTIC SCALAR IMAGING CONDITION
Robust 3D scalar imaging condition for elastic RTM Yuting Duan, presently at Shell International Exploration and Production Inc., formerly at Center for Wave Phenomena, Colorado School of Mines Paul Sava,
More informationLeast squares Kirchhoff depth migration: important details
Least squares Kirchhoff depth migration: important details Daniel Trad CREWES-University of Calgary Summary Least squares migration has been an important research topic in the academia for about two decades,
More informationLeast squares Kirchhoff depth migration with. preconditioning
Least squares Kirchhoff depth migration with preconditioning Aaron Stanton University of Alberta, Department of Physics, 4-83 CCIS, Edmonton AB T6G E (April, 3) Running head: Least Squares Migration ABSTRACT
More information3D angle decomposition for elastic reverse time migration Yuting Duan & Paul Sava, Center for Wave Phenomena, Colorado School of Mines
3D angle decomposition for elastic reverse time migration Yuting Duan & Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY We propose 3D angle decomposition methods from elastic reverse
More informationSUMMARY. denoise the original data at each iteration. This can be
A comparison of D reconstruction methods Aaron Stanton*, Nadia Kreimer, David Bonar, Mostafa Naghizadeh, and Mauricio Sacchi, Department of Physics, University of Alberta SUMMARY A comparison is made between
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 informationApex Shifted Radon Transform
Apex Shifted Radon Transform Daniel Trad* Veritas DGC Inc., Calgary, Alberta, Canada dtrad@veritasdgc.com ABSTRACT The Apex Shifted Radon transform (ASRT) is an extension of the standard hyperbolic RT,
More informationLeast Squares Kirchhoff Depth Migration: potentials, challenges and its relation to interpolation
Least Squares Kirchhoff Depth Migration: potentials, challenges and its relation to interpolation Daniel Trad CGG Summary Least Squares Migration (LSM), like interpolation, has the potential to address
More informationTarget-oriented wavefield tomography: A field data example
Target-oriented wavefield tomography: A field data example Yaxun Tang and Biondo Biondi ABSTRACT We present a strategy for efficient migration velocity analysis in complex geological settings. The proposed
More informationDownloaded 10/23/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
ACQUISITION APERTURE CORRECTION IN ANGLE-DOMAIN TOWARDS THE TRUE- REFLECTION RTM Rui Yan 1*, Huimin Guan 2, Xiao-Bi Xie 1, Ru-Shan Wu 1, 1 IGPP, Earth and Planetary Sciences Department, University of California,
More informationA least-squares shot-profile application of time-lapse inverse scattering theory
A least-squares shot-profile application of time-lapse inverse scattering theory Mostafa Naghizadeh and Kris Innanen ABSTRACT The time-lapse imaging problem is addressed using least-squares shot-profile
More informationScalar imaging condition for elastic reverse time migration
GEOPHYSICS, VOL. 80, NO. 4 (JULY-AUGUST 2015); P. S127 S136, 17 FIGS. 10.1190/GEO2014-0453.1 Scalar imaging condition for elastic reverse time migration Yuting Duan 1 and Paul Sava 1 ABSTRACT Polarity
More informationWriting Kirchhoff migration/modelling in a matrix form
Writing Kirchhoff migration/modelling in a matrix form Abdolnaser Yousefzadeh and John C. Bancroft Kirchhoff migration matrix ABSTRACT Kirchhoff prestack migration and modelling are linear operators. Therefore,
More informationInterpolation using asymptote and apex shifted hyperbolic Radon transform
Interpolation using asymptote and apex shifted hyperbolic Radon transform Amr Ibrahim, Paolo Terenghi and Mauricio D. Sacchi Department of Physics, University of Alberta PGS Abstract The asymptote and
More informationMulti-source Least-squares Migration of Gulf of Mexico Data
Multi-source Least-squares Migration of Gulf of Mexico Data Xin Wang. King Abdullah University of Science and Technology, Thuwal 955-69, Kingdom of Saudi Arabia Corresponding author is Xin Wang. E-mail
More informationRecovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples
Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples A. Ordoñez* (PGS), W.F. Sollner (PGS), T. Klüver (PGS) & L.G. Gelius (UiO) SUMMARY A joint
More informationThe. Reverse-time depth migration in elastic media. Introduction FOCUS ARTICLE. Zaiming Jiang. Coordinated by Rob Holt
FOCUS ARTICLE Coordinated by Rob Holt Reverse-time depth migration in elastic media Zaiming Jiang KEY SEISMIC SOLUTIONS, CALGARY, ALBERTA, CANADA The main objective of this paper is to sketch a reverse-time
More informationAmplitude and kinematic corrections of migrated images for nonunitary imaging operators
GEOPHYSICS, VOL. 69, NO. 4 (JULY-AUGUST 2004); P. 1017 1024, 19 FIGS. 10.1190/1.1778244 Amplitude and kinematic corrections of migrated images for nonunitary imaging operators Antoine Guitton ABSTRACT
More informationMulticomponent f-x seismic random noise attenuation via vector autoregressive operators
Multicomponent f-x seismic random noise attenuation via vector autoregressive operators Mostafa Naghizadeh and Mauricio Sacchi ABSTRACT We propose an extension of the traditional frequency-space (f-x)
More information3D angle gathers from wave-equation extended images Tongning Yang and Paul Sava, Center for Wave Phenomena, Colorado School of Mines
from wave-equation extended images Tongning Yang and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY We present a method to construct 3D angle gathers from extended images obtained
More informationLeast-squares Wave-Equation Migration for Broadband Imaging
Least-squares Wave-Equation Migration for Broadband Imaging S. Lu (Petroleum Geo-Services), X. Li (Petroleum Geo-Services), A. Valenciano (Petroleum Geo-Services), N. Chemingui* (Petroleum Geo-Services),
More informationInversion after depth imaging
Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.
More informationDownloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at
Recovering the Reflectivity Matrix and Angle-dependent Plane-wave Reflection Coefficients from Imaging of Multiples Alba Ordoñez PGS/UiO*, Walter Söllner PGS, Tilman Klüver PGS and Leiv J. Gelius UiO Summary
More informationProgress Report on: Interferometric Interpolation of 3D SSP Data
Progress Report on: Interferometric Interpolation of 3D SSP Data Sherif M. Hanafy ABSTRACT We present the theory and numerical results for interferometrically interpolating and extrapolating 3D marine
More informationAmplitude and kinematic corrections of migrated images for non-unitary imaging operators
Stanford Exploration Project, Report 113, July 8, 2003, pages 349 363 Amplitude and kinematic corrections of migrated images for non-unitary imaging operators Antoine Guitton 1 ABSTRACT Obtaining true-amplitude
More informationTarget-oriented wavefield tomography using demigrated Born data
Target-oriented wavefield tomography using demigrated Born data Yaxun Tang and Biondo Biondi ABSTRACT We present a method to reduce the computational cost of image-domain wavefield tomography. Instead
More informationAnisotropy-preserving 5D interpolation by hybrid Fourier transform
Anisotropy-preserving 5D interpolation by hybrid Fourier transform Juefu Wang and Shaowu Wang, CGG Summary We present an anisotropy-preserving interpolation method based on a hybrid 5D Fourier transform,
More informationPressure and shear wave separation of ocean bottom seismic data
Pressure and shear wave separation of ocean bottom seismic data Ohad Barak ABSTRACT PZ summation is a common method for separating the upgoing wavefield from the downgoing wavefield in data acquired by
More informationLab 3: Depth imaging using Reverse Time Migration
Due Wednesday, May 1, 2013 TA: Yunyue (Elita) Li Lab 3: Depth imaging using Reverse Time Migration Your Name: Anne of Cleves ABSTRACT In this exercise you will familiarize yourself with full wave-equation
More informationU043 3D Prestack Time Domain Full Waveform Inversion
U043 3D Prestack Time Domain Full Waveform Inversion D.V. Vigh* (WesternGeco), W.E.S. Starr (WesternGeco) & K.D. Kenneth Dingwall (WesternGeco) SUMMARY Despite the relatively high computational demand,
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 informationRefraction Full-waveform Inversion in a Shallow Water Environment
Refraction Full-waveform Inversion in a Shallow Water Environment Z. Zou* (PGS), J. Ramos-Martínez (PGS), S. Kelly (PGS), G. Ronholt (PGS), L.T. Langlo (PGS), A. Valenciano Mavilio (PGS), N. Chemingui
More informationEffects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines
Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic
More informationSUMMARY LEAST-SQUARES MIGRATION THEORY
Making the most out of the least (squares migration) Gaurav Dutta, Yunsong Huang, Wei Dai, Xin Wang, and G.T. Schuster King Abdullah University of Science and Technology SUMMARY Standard migration images
More informationModel-space vs. data-space normalization for recursive depth migration
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? Model-space vs. data-space normalization for recursive depth migration James Rickett 1 ABSTRACT Illumination problems caused by finite-recording
More informationAdaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited
Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears
More informationAngle-domain elastic reverse-time migration
1 Angle-domain elastic reverse-time migration 2 3 4 Jia Yan and Paul Sava Center for Wave Phenomena Colorado School of Mines 5 6 7 8 (January 4, 2008) GEO-2008-???? Running head: Elastic reverse-time migration
More informationReverse time migration in midpoint-offset coordinates
Stanford Exploration Project, Report 111, June 9, 00, pages 19 156 Short Note Reverse time migration in midpoint-offset coordinates Biondo Biondi 1 INTRODUCTION Reverse-time migration (Baysal et al., 198)
More informationAngle-dependent reflectivity by wavefront synthesis imaging
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 477 Angle-dependent reflectivity by wavefront synthesis imaging Jun Ji 1 ABSTRACT Elsewhere in this report, Ji and Palacharla (1994) show
More informationWave-equation inversion prestack Hessian
Stanford Exploration Project, Report 125, January 16, 2007, pages 201 209 Wave-equation inversion prestack Hessian Alejandro A. Valenciano and Biondo Biondi ABSTRACT The angle-domain Hessian can be computed
More informationPractical implementation of SRME for land multiple attenuation
Practical implementation of SRME for land multiple attenuation Juefu Wang* and Shaowu Wang, CGGVeritas, Calgary, Canada juefu.wang@cggveritas.com Summary We present a practical implementation of Surface
More informationOn the Scattering Effect of Lateral Discontinuities on AVA Migration
On the Scattering Effect of Lateral Discontinuities on AVA Migration Juefu Wang* and Mauricio D. Sacchi Department of Physics, University of Alberta, 4 Avadh Bhatia Physics Laboratory, Edmonton, AB, T6G
More informationSummary. Introduction
Yaofeng He and Ru-Shan Wu Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, C 9564 Summary We propose an approach to perform migration and imaging using secondary scattered
More informationHigh Resolution Imaging by Wave Equation Reflectivity Inversion
High Resolution Imaging by Wave Equation Reflectivity Inversion A. Valenciano* (Petroleum Geo-Services), S. Lu (Petroleum Geo-Services), N. Chemingui (Petroleum Geo-Services) & J. Yang (University of Houston)
More informationSUMMARY THEORY INTRODUCTION
Acoustic 3D least-squares reverse time migration using the energy norm Daniel Rocha, Paul Sava & Antoine Guitton Center for Wave Phenomena, Colorado School of Mines SUMMARY We propose a least-squares reverse
More informationTarget-oriented wave-equation inversion with regularization in the subsurface-offset domain
Stanford Exploration Project, Report 124, April 4, 2006, pages 1?? Target-oriented wave-equation inversion with regularization in the subsurface-offset domain Alejandro A. Valenciano ABSTRACT A complex
More informationTh G Surface-offset RTM Gathers - What Happens When the Velocity is Wrong?
Th G103 01 Surface-offset RTM Gathers - What Happens When the Velocity is Wrong? J.P. Montel* (CGG) SUMMARY Surface offset migrated common image gathers (SOCIGs) built by wave equation migration (WEM)
More informationFive Dimensional Interpolation:exploring different Fourier operators
Five Dimensional Interpolation:exploring different Fourier operators Daniel Trad CREWES-University of Calgary Summary Five-Dimensional interpolation has become a very popular method to pre-condition data
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 informationA Nuclear Norm Minimization Algorithm with Application to Five Dimensional (5D) Seismic Data Recovery
A Nuclear Norm Minimization Algorithm with Application to Five Dimensional (5D) Seismic Data Recovery Summary N. Kreimer, A. Stanton and M. D. Sacchi, University of Alberta, Edmonton, Canada kreimer@ualberta.ca
More informationMultichannel deconvolution imaging condition for shot-profile migration
Stanford Exploration Project, Report 113, July 8, 2003, pages 127 139 Multichannel deconvolution imaging condition for shot-profile migration Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT A significant
More informationAngle-domain elastic reverse-time migration
CWP-597 Angle-domain elastic reverse-time migration Jia Yan and Paul Sava Center for Wave Phenomena Colorado School of Mines ABSTRACT Multi-component data are not usually processed with specifically designed
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 informationLeast-squares RTM: Reality and possibilities for subsalt imaging Ping Wang*, Adriano Gomes, Zhigang Zhang, and Ming Wang, CGG
Least-squares RTM: Reality and possibilities for subsalt imaging Ping Wang*, Adriano Gomes, Zhigang Zhang, and Ming Wang, CGG Summary We investigated how current least-squares reverse time migration (LSRTM)
More informationPredicting rugged water-bottom multiples through wavefield extrapolation with rejection and injection
Predicting rugged water-bottom multiples through wavefield extrapolation with rejection and injection Benxi Ke ABSTRACT Although convolution-based and WEM-modeling-based methods for predicting surface-related
More informationReverse-time migration imaging with/without multiples
Reverse-time migration imaging with/without multiples Zaiming Jiang, John C. Bancroft, and Laurence R. Lines Imaging with/without multiples ABSTRACT One of the challenges with reverse-time migration based
More informationSeismic Attributes on Frequency-enhanced Seismic Data
Seismic Attributes on Frequency-enhanced Seismic Data Satinder Chopra* Arcis Corporation, Calgary, Canada schopra@arcis.com Kurt J. Marfurt The University of Oklahoma, Norman, US and Somanath Misra Arcis
More information= 0) is the 2-D Fourier transform of the field (1),, z = 0). The phase ) is defined in the dispersion relation as where
GEOPHYSICS, VOL. 61, NO. 5 (SEPTEMBER-OCTOBER 1996); P. 1412 1416, 4 FIGS. Short Note Prestack migration by split-step DSR. Alexander Mihai Popovici* INTRODUCTION The double-square-root (DSR) prestack
More information1D internal multiple prediction in a multidimensional world: errors and recommendations
1D internal multiple prediction 1D internal multiple prediction in a multidimensional world: errors and recommendations Pan Pan and Kris Innanen ABSTRACT Internal multiples are more difficult to estimate
More informationCommon-angle processing using reflection angle computed by kinematic pre-stack time demigration
Common-angle processing using reflection angle computed by kinematic pre-stack time demigration Didier Lecerf*, Philippe Herrmann, Gilles Lambaré, Jean-Paul Tourré and Sylvian Legleut, CGGVeritas Summary
More informationAVA analysis based on RTM angle-domain common image gather Rui Yan* and Xiao-Bi Xie, IGPP, University of California, Santa Cruz
AVA analysis based on RTM angle-domain common image gather Rui Yan* and Xiao-Bi Xie, IGPP, University of California, Santa Cruz Summary We propose an alternative approach of AVA analysis based on RTM angle-domain
More informationApplication of 3D source deghosting and designature to deep-water ocean bottom node data Xu Li*, Jing Yang, Hui Chen, Melanie Vu, and Ping Wang (CGG)
Xu Li*, Jing Yang, Hui Chen, Melanie Vu, and Ping Wang (CGG) Summary Compared to towed-streamer data, deep-water ocean bottom node (OBN) data by nature have a broader bandwidth; however, the presence of
More informationEquivalence of source-receiver migration and shot-profile migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 109 117 Short Note Equivalence of source-receiver migration and shot-profile migration Biondo Biondi 1 INTRODUCTION At first glance, shot
More informationAngle Gathers for Gaussian Beam Depth Migration
Angle Gathers for Gaussian Beam Depth Migration Samuel Gray* Veritas DGC Inc, Calgary, Alberta, Canada Sam Gray@veritasdgc.com Abstract Summary Migrated common-image-gathers (CIG s) are of central importance
More informationMigration velocity analysis of incomplete/irregular data using least squares Kirchhoff migration
Migration velocity analysis of incomplete/irregular data using least squares Kirchhoff migration Abdolnaser Yousefzadeh, John C. Bancroft, and Gary F. Margrave ABSTRACT In order to attenuate the migration
More informationSimultaneous source separation via robust time variant Radon operators
Simultaneous source separation via robust time variant Radon operators Summary Amr Ibrahim, University of Alberta, Edmonton, Canada amr.ibrahim@ualberta.ca and Mauricio D. Sacchi, University of Alberta,
More informationF031 Application of Time Domain and Single Frequency Waveform Inversion to Real Data
F031 Application of Time Domain and Single Frequency Waveform Inversion to Real Data D Yingst (ION Geophysical), C. Wang* (ION Geophysical), J. Park (ION Geophysical), R. Bloor (ION Geophysical), J. Leveille
More informationMigration from a non-flat datum via reverse-time extrapolation
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 50 Migration from a non-flat datum via reverse-time extrapolation Gopal Palacharla 1 ABSTRACT Land surveys usually have elevation changes,
More informationIllumination-based normalization for wave-equation depth migration
Illumination-based normalization for wave-equation depth migration James E. Rickett ChevronTexaco Exploration and Production Technology Company, 6001 Bollinger Canyon Road, San Ramon, CA 94583-2324 formerly
More informationCFP migration; practical aspects
13 CFP migration; practical aspects Jan Thorbecke 1 13.1 Introduction In the past year the Common Focus Point (CFP) technology has become an important paradigm in the DELPHI research program. The interpretation
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 informationWe N Converted-phase Seismic Imaging - Amplitudebalancing Source-independent Imaging Conditions
We N106 02 Converted-phase Seismic Imaging - Amplitudebalancing -independent Imaging Conditions A.H. Shabelansky* (Massachusetts Institute of Technology), A.E. Malcolm (Memorial University of Newfoundland)
More informationA low rank based seismic data interpolation via frequencypatches transform and low rank space projection
A low rank based seismic data interpolation via frequencypatches transform and low rank space projection Zhengsheng Yao, Mike Galbraith and Randy Kolesar Schlumberger Summary We propose a new algorithm
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 informationMinimum weighted norm wavefield reconstruction for AVA imaging
Geophysical Prospecting, 25, 53, 787 81 Minimum weighted norm wavefield reconstruction for AVA imaging Mauricio D. Sacchi and Bin Liu Department of Physics, University of Alberta, Abadh Bhatia Physics
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 informationSuccesses and challenges in 3D interpolation and deghosting of single-component marinestreamer
Successes and challenges in 3D interpolation and deghosting of single-component marinestreamer data James Rickett*, Schlumberger Gould Research Summary Combining deghosting with crossline interpolation
More informationDownloaded 05/09/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Elastic converted-wave path migration for subsalt imaging Ru-Shan Wu*, Rui Yan, Xiao-Bi Xie, Modeling and Imaging Laboratory, Earth and Planetary Sciences/IGPP, University of California, Santa Cruz, David
More information3D PRESTACK DEPTH MIGRATION WITH FACTORIZATION FOUR-WAY SPLITTING SCHEME
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Supp, Pages 183 196 c 5 Institute for Scientific Computing and Information 3D PRESTACK DEPTH MIGRATION WITH FACTORIZATION FOUR-WAY SPLITTING
More informationBuilding starting model for full waveform inversion from wide-aperture data by stereotomography
Building starting model for full waveform inversion from wide-aperture data by stereotomography Vincent Prieux 1, G. Lambaré 2, S. Operto 1 and Jean Virieux 3 1 Géosciences Azur - CNRS - UNSA, France;
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 informationy Input k y2 k y1 Introduction
A regularization workflow for the processing of cross-spread COV data Gordon Poole 1, Philippe Herrmann 1, Erika Angerer 2, and Stephane Perrier 1 1 : CGG-Veritas, 2 : OMV Data regularization is critical
More informationMitigation of the 3D Cross-line Acquisition Footprint Using Separated Wavefield Imaging of Dual-sensor Streamer Seismic
Mitigation of the 3D Cross-line Acquisition Footprint Using Separated Wavefield Imaging of Dual-sensor Streamer Seismic A.S. Long* (PGS), S. Lu (PGS), D. Whitmore (PGS), H. LeGleut (PGS), R. Jones (Lundin
More informationStacking angle-domain common-image gathers for normalization of illumination a
Stacking angle-domain common-image gathers for normalization of illumination a a Published in Geophysical Prospecting, 59, 244-255 (2011) Guochang Liu, China University of Petroleum-Beijing and The University
More informationDownloaded 01/20/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
D Multi-source Least-squares Reverse Time Migration Wei Dai, C. Boonyasiriwat and Gerard T. Schuster, King Abdullah University of Science and Technology Downloaded 1//1 to 19.171.17.1. Redistribution subject
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 informationSeismic Data Preconditioning for Improved Reservoir Characterization (Inversion and Fracture Analysis)
Seismic Data Preconditioning for Improved Reservoir Characterization (Inversion and Fracture Analysis) Darren Schmidt *, Hampson-Russell, Calgary, AB, Canada Darren.Schmidt@cggveritas.com Alicia Veronesi,
More informationSUMMARY. These two projections are illustrated in the equation
Mitigating artifacts in Projection Onto Convex Sets interpolation Aaron Stanton*, University of Alberta, Mauricio D. Sacchi, University of Alberta, Ray Abma, BP, and Jaime A. Stein, Geotrace Technologies
More informationAzimuth Moveout Transformation some promising applications from western Canada
Azimuth Moveout Transformation some promising applications from western Canada Satinder Chopra and Dan Negut Arcis Corporation, Calgary, Canada Summary Azimuth moveout (AMO) is a partial migration operator
More informationA full-wave equation based seismic illumination analysis method
A full-wave equation based seismic illumination analysis method Xiao-Bi Xie and Hui Yang University of California, Santa Cruz, CA 95064, USA Summary We propose a full-wave equation based method for seismic
More informationMachine-learning Based Automated Fault Detection in Seismic Traces
Machine-learning Based Automated Fault Detection in Seismic Traces Chiyuan Zhang and Charlie Frogner (MIT), Mauricio Araya-Polo and Detlef Hohl (Shell International E & P Inc.) June 9, 24 Introduction
More informationIterative resolution estimation in Kirchhoff imaging
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 369?? Iterative resolution estimation in Kirchhoff imaging Robert G. Clapp, Sergey Fomel, and Marie Prucha 1 ABSTRACT We apply iterative
More informationLeast-squares migration/inversion of blended data
Least-squares migration/inversion of blended data Yaxun Tang and Biondo Biondi ABSTRACT We present a method based on least-squares migration/inversion to directly image data collected from recently developed
More informationG021 Subsalt Velocity Analysis Using One-Way Wave Equation Based Poststack Modeling
G021 Subsalt Velocity Analysis Using One-Way Wave Equation Based Poststack Modeling B. Wang* (CGG Americas Inc.), F. Qin (CGG Americas Inc.), F. Audebert (CGG Americas Inc.) & V. Dirks (CGG Americas Inc.)
More informationDownloaded 09/20/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
Joint SRME and model-based water-layer demultiple for ocean bottom node Hui Huang*, Ping Wang, Jing Yang, Hui Chen (CGG); Pierre-Olivier Ariston, Imtiaz Ahmed, and Nick Bassett (BP) Summary Ocean bottom
More informationAttenuation of water-layer-related multiples Clement Kostov*, Richard Bisley, Ian Moore, Gary Wool, Mohamed Hegazy, Glenn Miers, Schlumberger
Clement Kostov*, Richard Bisley, Ian Moore, Gary Wool, Mohamed Hegazy, Glenn Miers, Schlumberger Summary We present a method for modeling and separation of waterlayer-related multiples in towed streamer
More informationDirections in 3-D imaging - Strike, dip, both?
Stanford Exploration Project, Report 113, July 8, 2003, pages 363 369 Short Note Directions in 3-D imaging - Strike, dip, both? Marie L. Clapp 1 INTRODUCTION In an ideal world, a 3-D seismic survey would
More informationG009 Scale and Direction-guided Interpolation of Aliased Seismic Data in the Curvelet Domain
G009 Scale and Direction-guided Interpolation of Aliased Seismic Data in the Curvelet Domain M. Naghizadeh* (University of Alberta) & M. Sacchi (University of Alberta) SUMMARY We propose a robust interpolation
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 information