Least squares Kirchhoff depth migration: important details
|
|
- Moris Riley
- 5 years ago
- Views:
Transcription
1 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, but only recently it has attracted interest from the industry. The main reason is that from a practical point of view its ratio of benefit/cost has not been sufficient for its use in seismic exploration. Another problem is that these benefits are mixed with effects from the filtering techniques used to regularize the inversion, which are computationally much cheaper. In this paper, I discuss some challenges with least squares Kirchhoff depth migration. This algorithm, although less precise than the more popular least squares reverse time migration, has the advantage of being fast enough to be applied in a production environment, and flexible enough to be applied without data regularization. This last characteristic makes it a good candidate to understand benefits in terms of footprint acquisition and aliasing. In addition, its limitations in terms of modelling/imaging accuracy make more evident some problems that exist but are often ignored when using reverse time migration with synthetic data. Introduction Despite the existence of more accurate imaging algorithms, the most common migration algorithm for land seismic data is still Kirchhoff migration, mostly because of its flexibility to adapt to irregular sampling and its efficiency. Kirchhoff migration is basically a weighted mapping between data space (the original acquired data), and the image space (migrated section), where samples are weighted and shifted temporally and spatially. These weights come from the mathematics of integration and wave propagation, usually assuming infinite aperture, and taking only partially into account the exact locations of shots and receivers. Hard work over many decades has led to weights that produce a good approximation to true amplitude reflectivity. However, these weights have limitations because they do not contain information about the sampling operator and acquisition illumination. Least squares migration (LSMIG) is a different approach to migration where, rather than mapping and weighting to find the model, a modeling operator is inverted by optimization algorithms (Nemeth et al., 1999; Kuhl and Sacchi, 2003). The inverse of this operator applied to seismic data produces an image that by design can predict the acquired data in a least squares sense. Although the inverse of this modeling operator is closely related to the imaging (mapping) operator, it represents a different approach to the problem and it takes into account the actual geometry (sampling and aperture) of the data. By inverting a modeling operator that contains the acquisition pattern, we can in principle remove acquisition and illumination effects. Probably the most attractive characteristic of least squares migration is that it falls into a wide range of algorithmic techniques used not only in seismic but in physics in general. The mathematics and physics of LSMIG allows one to include concepts like prior information, in the form of mathematical constraints that help to obtain better images. In practice, however, LSMIG is a difficult process because of the high cost of applying the modeling/adjoint operator multiple times, and the challenging goal of creating an image that not only looks reasonable but also can predict the data. Although stacking samples (mapping) is a very robust process, inverting a multidimensional operator is usually an ill-conditioned problem. From the practical point of view the advantages of LSMIG have not been shown in the industry except for specific examples, mostly 2D synthetic examples. Its application to 3D seismic processing is only starting to be attempted as an alternative to traditional methods and it remains as an elusive application of academic research. This is particularly true for land data sets, where noise is complex and correct signal amplitudes are very difficult to predict correctly. However, many possible advantages over conventional GeoConvention
2 migration keep motivating researchers to solve these difficulties. Among these possibilities, we can mention (1) elimination of footprint noise caused by poor sampling and (2) better amplitudes in the image by compensation of illumination problems. In addition, the availability of the LSMIG engine provides a powerful and physical tool for interpolation, and de-noising (Trad, 2015). Theory The LSMIG implementation in this work uses a general philosophy widely used to solve for geophysical problems that originates from the Stanford Exploration Project and is nicely explained in Claerbout (1992). Given a transformation of some model m to some data d through some mathematical operator L, d = Lm. (1) we can approximately solve the inverse problem of recovering the model that produced the data by inverting the operator in a least squares sense. The meaning of the different components of equation (1) depends on the details of the operator and are often associated with well-known transformations. For the case of LSMIG, L is a synthesis or modeling operator, that generates data from reflectivity using Green functions as basis functions. Any transformation could be used to generate the same data by this procedure, independently of whether it is physical or not. What is different and requires some consideration is the meaning we assign to the model. Even for very similar operators, small variations can change what the model represents. For example, by changing the weights in the operator we can switch the model from reflectivity (Kirchhoff) to velocity perturbations (Born). Once a meaningful operator is chosen, we can estimate the model by inverting the operator. To invert operators, we define a cost or objective function, which is a mathematical expression that measures the undesired characteristics of the model. The most common is to find a model that honors the data in a least error sense, measured on some norm, and has a minimum of information not required by the data. This statement of goals is commonly presented as minimize W m m p p subject to W d (d Lm) = φ d (2) where φ d is some estimate of the noise level in the data plus a residual due to the failure of the proposed model to explain the data. p and q indicate that different norms can be applied to measure the norm of vectors. W d could be a matrix or vector of data weights, often a diagonal matrix containing the inverse of the standard deviation of the data, but more generally it could be any kind of filter that leaves out bad data. W m is an operator of model weights that can be customized to enhance our preferences regarding the model. For example, W m might be a gradient (roughening) operator, then minimizing W m m will produce a smooth model. After absorving W d and W m into L by a change of variables (Trad et al., 2003), the solution to this problem results on the usual LS formulation: m=(l H L) -1 L H d (3) were the (L H L) operator, known as the Hessian, is iteratively deconvolved from the migration result (L H d). For LSMIG we usually need to favor, but not enforce, smooth models that lead to continuous structures. A typical choice for Kirchhoff type of migrations, is to use W m to enforce smooth changes of reflectivity with offset or angle. This can be done as in Nemeth et al. (1999) by taking differences between consecutive reflectivities at different offsets, as in Kuhl and Sacchi (2003) by a sparse tau p transform of the reflectivity, or as in Moghaddam and Herrmann (2005) by a sparse curvelet transform of the reflectivity. A generic choice for W m is the inverse of a multidimensional operator that takes a multioffset migrated volume and performs smoothing, for example a triangular filter, a dip filter, a fxy filter or a Radon transform. In Fomel (2007), this formulation is generalized with the name of shaping filter. Weighting It is not trivial to achieve the optimal weighting scheme for LS Kirchhoff migration. A common approach is to use standard migration weights with a migration operator, design its adjoint such that it passes the adjoint test (Claerbout, 1992), and use these two operators into a CG algorithm. Weights in migration algorithms often take the form of deconvolution (DC) weights, which come from the imaging principle GeoConvention
3 (Claerbout, 1971) as a ratio of scattered (u s) over incident (u i) wavefields. This leads to migration weights that are proportional to a ratio of shot and receiver wavefield amplitudes, or, after approximating for constant velocity, to a ratio of traveltime inverses. (Dellinger et al., 1999; Zhang et al., 2000). On the other hand, Nemeth et al. (1999) sets the operator L as a modelling rather than a migration (Casasanta, personal communication). This approach changes the imaging condition from deconvolution (migration) to cross-correlation (modelling). In the cross-correlation imaging condition (XC) the LS weights are proportional to the product of amplitudes, which leads, after using similar approximations, to a product of traveltimes instead of a ratio. Testing with simple flat synthetics supports this interpretation. Figure 1a-b illustrates that deconvolution weights produce the correct amplitudes for migration (a), but as we iterate during the inversion, the amplitudes get distorted (b). Figure 1c-d show that cross-correlation weights produce the wrong image for migration (c), but the amplitudes become correct after a few iterations (d). Another related issue is whether the phase filter operator should be included inside the modelling - migration operator so it would act for every iteration, or should be only pre-applied as a data preconditioner, or not applied at all. The exact formulation for the phase shift operator depends on the dimensionality of the data acquisition, being equal to iω for 3D, i ω for 2.5D or ω for 2D (equations 2 in Bleistein and Gray (2001)). In this paper, I refer to this filter simply as ω filter. Figure 2 compares the amplitude spectra. These tests show that if the ω filter is not included inside the operator, then the spectrum of the image shift towards higher frequencies with iterations. Another factor to consider with data weights is whether to include or not velocity weights. In theory amplitude weights contain a power of velocity (v 2 for 2D, v 3 for 3D). Although in migration this factor can be replaced by a depth gain after migration, in LSMIG it has a different effect because by including the velocities inside the weights the convergence is focused to areas with higher velocity. This is one of many situations where simplifications taken in migration do not carry over to LSMIG. Noise control One of the most important problems in LSMIG is the accumulation of noise with iterations. This noise comes from broken contributions to the model as the algorithm tries to predict small details on the data that are not properly mapped by the migration operator, in other words, when the operator does not approximate the physics correctly. For example, this effect is accentuated with discontinuities or gaps in Kirchhoff depth migration traveltime tables, which make the operator a poor representation of the physics. Conversely, this noise disappears if the data are created by the same numerical approximation used in the migration engine. For this reason, LSMIG noise is more obvious with a Kirchhoff migration than with reverse time migration (RTM). However, RTM may not show this problem if tested with finite difference data, but the noise may appear in real data processing as the numerical approximation deviate from the actual physics that produced the data. As usual in ill-conditioned inversions, this problem requires some form of noise control by means of numerical regularization. Regularization in this work is performed by a multidimensional triangular filter (the W m operator in equations 2). The filter acts on the inline, crossline and offset dimensions. The length of the filter, which is independent across dimensions, permits one to control the noise effectively, but it is limited by the structural complexity of the image. In Figure 3a we see one iteration of LSMIG (equivalent to migration) for the Marmousi data set, Figure 3b we see the same result after 10 iterations when no noise control is applied, and Figure 3c shows the same with a triangle filter along crosslines (cdps) and offsets. The length of the filter was very mild to preserve structure (length along inline was 5, length along offset was 5). Increasing the length would make the result cleaner, but some faults would be smoothed. We need to emphasize that improvements in the image we seek should not be a pure consequence of the regularization filter but of the inversion of the Hessian operator L H L. Convergence A critical aspect for a successful application of LSMIG in an industrial setting is to improve its convergence. One factor that has a large effect on convergence is the dimensionality of the model. When the model is defined in terms of offset gathers instead of their summation, convergence improves GeoConvention
4 dramatically. This is understandable because of less mixing, which makes prediction easier, and therefore convergence faster. Using offset gathers has the effect of subdividing the problem into smaller sub-inversions. Similarly, convergence for a part of the data is much faster than for the whole data. The trade-off, however, is that by using part of the data inversion is only partially accomplished, even if data fitting is better. Conclusions Although straightforward in principle, there are many options related to the operator design whose choice can have significant influence in results and efficiency, and are not obvious from the theoretical point of view. In particular it is important to understand the operator weights in terms of phase and amplitude as well as regularization operators to control the noise during iterations. Acknowledgements My gratitude to Lorenzo Casasanta and Sam Gray for many enlightful discussions about migration during my work at CGG. I also thank CREWES sponsors for contributing to this seismic research. I also gratefully acknowledge support from NSERC (Natural Science and Engineering Research Council of Canada) through the grant CRDPJ a) b) c) d) Figure 1: Top DC weights, a) migration, b) LSMIG, bottom XC weights, c) migration d) LSMIG Figure 2: Effect of phase filter. Top DC weights, a) pre-applied b) inside operator, bottom XC weights, c) pre-applied d) inside operator Figure 2: Noise with iterations. a) 1 iteration b) 9 iterations without noise control c) 9 iterations with filter inside operator. GeoConvention
5 References Bleistein, N., and Gray, S. H., 2001, From the hagedoorn imaging technique to kirchhoff migration and inversion: Geophysical Prospecting, 49, No. 6, Claerbout, J., 1992, Earth sounding analysis, Processing versus inversion: Blackwell Scientific Publications. Claerbout, J. F., 1971, Toward a unified theory of reflector mapping: GEOPHYSICS, 36, No. 3, Dellinger, J., Murphy, G., Etgen, J., Fei, T., and Gray, S., 1999, Efficient 2.5-d trueâa Ramplitude migration:ˇ The Leading Edge, 18, No. 8, Docherty, P., 1991, A brief comparison of some kirchhoff integral formulas for migration and inversion: GEOPHYSICS, 56, No. 8, Fomel, S., 2007, Shaping regularization in geophysical-estimation problems: GEOPHYSICS, 72, No. 2, R29 R36. Kuhl, H., and Sacchi, M. D., 2003, Least squares wave equation migration for avp/ava inversion: GEOPHYSICS, 68, No. 1, Lumley, D. E., Claerbout, J. F., and Bevc, D., 1994, AntiâA Raliased kirchhoff 3⡠A RD migration, ˇ Moghaddam, P. P., and Herrmann, F. J., 2005, Migration preconditioning with Curvelets, Nemeth, T., Wu, C., and Schuster, G. T., 1999, Least squares migration of incomplete reflection data: GEOPHYSICS, 64, No. 1, Trad, D., 2015, Least squares kirchhoff depth migration: implementation, challenges, and opportunities, in 2015 SEG Annual Meeting, Society of Exploration Geophysicists, Trad, D., Ulrych, T., and Sacchi, M., 2003, Latest views of the sparse radon transform: GEOPHYSICS, 68, No. 1, Zhang, Y., Gray, S., and Young, J., 2000, Exact and approximate weights for Kirchhoff migration, GeoConvention
Five 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationAngle-gather time migration a
Angle-gather time migration a a Published in SEP report, 1, 141-15 (1999) Sergey Fomel and Marie Prucha 1 ABSTRACT Angle-gather migration creates seismic images for different reflection angles at the reflector.
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 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 informationLeast. Least squares migration of converted wave seismic data. Introduction FOCUS ARTICLE. Aaron Stanton and Mauricio D. Sacchi
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
More informationTarget-oriented wave-equation inversion
Stanford Exploration Project, Report 120, May 3, 2005, pages 23 40 Target-oriented wave-equation inversion Alejandro A. Valenciano, Biondo Biondi, and Antoine Guitton 1 ABSTRACT A target-oriented strategy
More informationSEG/New Orleans 2006 Annual Meeting
3-D tomographic updating with automatic volume-based picking Dimitri Bevc*, Moritz Fliedner, Joel VanderKwaak, 3DGeo Development Inc. Summary Whether refining seismic images to evaluate opportunities in
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 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 informationProcessing converted-wave data in the tau-p domain: rotation toward the source and moveout correction
τ-p domain converted-wave processing Processing converted-wave data in the tau-p domain: rotation toward the source and moveout correction Raul Cova and Kris Innanen ABSTRACT The asymmetry of the converted-wave
More informationStochastic conjugate gradient method for least-square seismic inversion problems Wei Huang*, Hua-Wei Zhou, University of Houston
Stochastic conjugate gradient method for least-square seismic inversion problems Wei Huang*, Hua-Wei Zhou, University of Houston Summary With the development of computational power, there has been an increased
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 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 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 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 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 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 informationThe power of stacking, Fresnel zones, and prestack migration
The power of stacking, Fresnel zones, and prestack migration John C. Bancroft and Shuang Sun ABSTRACT The stacking of common midpoint (CMP) gathers assumes the presence of specula reflection energy and
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 informationApproximate Constant Density Acoustic Inverse Scattering Using Dip-Dependent Scaling Rami Nammour and William Symes, Rice University
Approximate Constant Density Acoustic Inverse Scattering Using Dip-Dependent Scaling Rami Nammour and William Symes, Rice University SUMMARY This abstract presents a computationally efficient method to
More informationFull waveform inversion by deconvolution gradient method
Full waveform inversion by deconvolution gradient method Fuchun Gao*, Paul Williamson, Henri Houllevigue, Total), 2012 Lei Fu Rice University November 14, 2012 Outline Introduction Method Implementation
More informationTarget-oriented computation of the wave-equation imaging Hessian
Stanford Exploration Project, Report 117, October 23, 2004, pages 63 77 Target-oriented computation of the wave-equation imaging Hessian Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT A target-oriented
More informationSubsalt illumination analysis using RTM 3D dip gathers Zhengxue Li *, Bing Tang and Shuo Ji, CGGVeritas
Zhengxue Li *, Bing Tang and Shuo Ji, CGGVeritas Summary Reverse Time Migration (RTM) is now the preferred option for subsalt imaging in deep water Gulf of Mexico, and its 3D angle gather output plays
More informationMathematics of Multidimensional Seismic Imaging, Migration, and Inversion
N. Bleistein J.K. Cohen J.W. Stockwell, Jr. Mathematics of Multidimensional Seismic Imaging, Migration, and Inversion With 71 Illustrations Springer Contents Preface List of Figures vii xxiii 1 Multidimensional
More informationSource-receiver migration of multiple reflections
Stanford Exploration Project, Report 113, July 8, 2003, pages 75 85 Source-receiver migration of multiple reflections Guojian Shan 1 ABSTRACT Multiple reflections are usually considered to be noise and
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 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 informationRobust Seismic Image Amplitude Recovery Using Curvelets
Robust Seismic Image Amplitude Recovery Using Curvelets Peyman P. Moghaddam Joint work with Felix J. Herrmann and Christiaan C. Stolk UNIVERSITY OF BRITISH COLUMBIA University of Twente References Curvelets
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 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 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 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 informationRemoving blended sources interference using Stolt operator Amr Ibrahim, Department of Physics, University of Alberta,
Removing blended sources interference using Stolt operator Amr Ibrahim, Department of Physics, University of Alberta, amr.ibrahim@ualberta.ca Summary We use Stolt migration/de-migration operator to remove
More informationRobustness of the scalar elastic imaging condition for converted waves
CWP-830 Robustness of the scalar elastic imaging condition for converted waves Yuting Duan & Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT For elastic reverse-time migration, one
More 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 informationGeometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT. August Lau and Chuan Yin.
Geometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT August Lau and Chuan Yin January 6, 2017 Abstract The goal of seismic processing is to convert input data
More informationMulticomponent land data pre-processing for FWI: a benchmark dataset
Multicomponent land data pre-processing for FWI: a benchmark dataset Raul Cova, Bernie K. Law and Kris Innanen CRWES/University of Calgary Summary Successful full-waveform inversion (FWI) studies using
More 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 informationEfficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers
Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers J.W.C. Sherwood* (PGS), K. Sherwood (PGS), H. Tieman (PGS), R. Mager (PGS) & C. Zhou
More 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 informationWhen is anti-aliasing needed in Kirchhoff migration? a
When is anti-aliasing needed in Kirchhoff migration? a a Published in SEP report, 80, 467-476 (1994) Dimitri Bevc and David E. Lumley ABSTRACT We present criteria to determine when numerical integration
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 informationP071 Land Data Regularization and Interpolation Using Azimuth Moveout (AMO)
P071 Land Data Regularization and Interpolation Using Azimuth Moveout (AMO) A.M. Popovici* (3DGeo Development Inc.), S. Crawley (3DGeo), D. Bevc (3DGeo) & D. Negut (Arcis Processing) SUMMARY Azimuth Moveout
More information# of PDE solves flops conventional 2N s N s N hx N hy N hz this paper 2N x N r N s
AVA analysis and geological dip estimation via two-way wave-equation based extended images Rajiv Kumar, Tristan van Leeuwen and Felix J. Herrmann University of British Columbia, Department of Earth, Ocean
More informationModel parametrization strategies for Newton-based acoustic full waveform
Model parametrization strategies for Newton-based acoustic full waveform inversion Amsalu Y. Anagaw, University of Alberta, Edmonton, Canada, aanagaw@ualberta.ca Summary This paper studies the effects
More informationSUMMARY 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 informationSeismic data interpolation beyond aliasing using regularized nonstationary autoregression a
Seismic data interpolation beyond aliasing using regularized nonstationary autoregression a a Published in Geophysics, 76, V69-V77, (2011) Yang Liu, Sergey Fomel ABSTRACT Seismic data are often inadequately
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 informationEquivalent offset migration: the implementation and application update
Equivalent offset migration: the implementation and application update Xinxiang Li, Yong Xu and John C. Bancroft INTRODUCTION We have some improvements about equivalent offset migration (EOM) method. In
More informationSUMMARY INTRODUCTION THEORY
Illumination compensation for subsalt image-domain wavefield tomography Tongning Yang, Center for Wave Phenomena, Colorado School of Mines, Jeffrey Shragge, School of Earth and Environment, University
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 informationG012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry
G012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry D.F. Halliday* (Schlumberger Cambridge Research), P.J. Bilsby (WesternGeco), J. Quigley (WesternGeco) & E. Kragh (Schlumberger
More informationSummary. Introduction. Source separation method
Separability of simultaneous source data based on distribution of firing times Jinkun Cheng, Department of Physics, University of Alberta, jinkun@ualberta.ca Summary Simultaneous source acquisition has
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 informationSeparation of specular reflection and diffraction images in Kirchhoff depth migration Faruq E Akbar and Jun Ma, SEIMAX Technologies, LP
Separation of specular reflection and diffraction images in Kirchhoff depth migration Faruq E Akbar and Jun Ma, SEIMAX Technologies, LP Summary Seismic diffractions may occur from faults, fractures, rough
More 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 informationSeismic data Interpolation in the Continuous Wavenumber Domain, Flexibility and Accuracy
Seismic data Interpolation in the Continuous Wavenumber Domain, Flexibility and Accuracy Ye Zheng Geo-X Exploration Services Inc. Summary Our recently developed algorithm, ASFT (Arbitrarily Sampled Fourier
More informationWe G High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation
We G103 05 High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation G. Hilburn* (TGS), Y. He (TGS), Z. Yan (TGS) & F. Sherrill (TGS) SUMMARY An approach
More informationWave-equation migration from topography: Imaging Husky
Stanford Exploration Project, Report 123, October 31, 2005, pages 49 56 Short Note Wave-equation migration from topography: Imaging Husky Jeff Shragge 1 INTRODUCTION Imaging land seismic data is wrought
More informationCLASSIFICATION OF MULTIPLES
Introduction Subsurface images provided by the seismic reflection method are the single most important tool used in oil and gas exploration. Almost exclusively, our conceptual model of the seismic reflection
More informationA Review of Current Marine Demultiple Techniques with Examples from the East Coast of Canada
A Review of Current Marine De Techniques with Examples from the East Coast of Canada ABSTRACT R.Brooymans*, T.Mojesky and L.Pham CGG Canada Services Ltd., Calgary rbrooymans@ca.cgg.com Introduction Multiple
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 informationImprovements in time domain FWI and its applications Kwangjin Yoon*, Sang Suh, James Cai and Bin Wang, TGS
Downloaded 0/7/13 to 05.196.179.38. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ Improvements in time domain FWI and its applications Kwangjin Yoon*,
More 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 informationAzimuth Moveout (AMO) for data regularization and interpolation. Application to shallow resource plays in Western Canada
Azimuth Moveout (AMO) for data regularization and interpolation. Application to shallow resource plays in Western Canada Dan Negut, Samo Cilensek, Arcis Processing, Alexander M. Popovici, Sean Crawley,
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 informationAngle-domain parameters computed via weighted slant-stack
Angle-domain parameters computed via weighted slant-stack Claudio Guerra 1 INTRODUCTION Angle-domain common image gathers (ADCIGs), created from downward-continuation or reverse time migration, can provide
More informationLessons in Spatial Sampling Or... Does Anybody Know the Shape of a Wavefield?
Lessons in Spatial Sampling Or... Does Anybody Know the Shape of a Wavefield? Norm Cooper - Mustagh Resources Ltd., Calgary, Canada In a homogeneous, isotropic, infinite half space (remember those regimes
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 informationShaping regularization in geophysical-estimation problems
GEOPHYSICS, VOL. 7, NO. MARCH-APRIL 7 ; P. R9 R6, FIGS., TABLE..9/.476 Shaping regularization in geophysical-estimation problems Sergey Fomel ABSTRACT Regularization is a required component of geophysical-estimation
More informationF020 Methods for Computing Angle Gathers Using RTM
F020 Methods for Computing Angle Gathers Using RTM M. Vyas* (WesternGeco), X. Du (WesternGeco), E. Mobley (WesternGeco) & R. Fletcher (WesternGeco) SUMMARY Different techniques can be used to compute angle-domain
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 informationWell-log derived step lengths in Full Waveform Inversion Tunde Arenrin Gary Margrave
Well-log derived step lengths in Full Waveform Inversion Tunde Arenrin Gary Margrave Motivation»Being able to make FWI a part of the seismic processing work flow under current computational capability.»incorporating
More informationIMAGING USING MULTI-ARRIVALS: GAUSSIAN BEAMS OR MULTI-ARRIVAL KIRCHHOFF?
IMAGING USING MULTI-ARRIVALS: GAUSSIAN BEAMS OR MULTI-ARRIVAL KIRCHHOFF? Summary Samuel H. Gray* Veritas DGC Inc., 715 Fifth Ave. SW, Suite 2200, Calgary, AB Sam_gray@veritasdgc.com Carl Notfors Veritas
More informationDeconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC
Deconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC SUMMARY We use the recently introduced multiscale and multidirectional curvelet transform to exploit the
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 informationWe Fast Beam Migration Using Plane Wave Destructor (PWD) Beam Forming SUMMARY
We-02-12 Fast Beam Migration Using Plane Wave Destructor (PWD) Beam Forming A.M. Popovici* (Z-Terra Inc.), N. Tanushev (Z-Terra Inc.), I. Sturzu (Z-Terra Inc.), I. Musat (Z-Terra Inc.), C. Tsingas (Saudi
More informationImplicit and explicit preconditioning for least squares nonstationary phase shift
Implicit and explicit preconditioning for least squares nonstationary phase shift Marcus R. Wilson and Robert J. Ferguson ABSTRACT An implicit preconditioned conjugate gradient scheme is derived to implement
More informationComments on wavefield propagation using Reverse-time and Downward continuation
Comments on wavefield propagation using Reverse-time and Downward continuation John C. Bancroft ABSTRACT Each iteration a of Full-waveform inversion requires the migration of the difference between the
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 informationSeisSpace Software. SeisSpace enables the processor to be able focus on the science instead of being a glorified data manager.
SeisSpace Software OVERVIEW Delivering large volumes of data quickly and accurately remains the hallmark of any seismic processing system. With intuitive analysis tools, state-of-the-art geophysical algorithms
More information3D image-domain wavefield tomography using time-lag extended images
CWP-748 3D image-domain wavefield tomography using time-lag extended images Tongning Yang and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT Image-domain wavefield tomography is
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 informationEARTH SCIENCES RESEARCH JOURNAL
EARTH SCIENCES RESEARCH JOURNAL Earth Sci. Res. J. Vol. 10, No. 2 (December 2006): 117-129 ATTENUATION OF DIFFRACTED MULTIPLES WITH AN APEX-SHIFTED TANGENT- SQUARED RADON TRANSFORM IN IMAGE SPACE Gabriel
More information