Amplitude and kinematic corrections of migrated images for non-unitary imaging operators
|
|
- Ilene Wade
- 6 years ago
- Views:
Transcription
1 Stanford Exploration Project, Report 113, July 8, 2003, pages Amplitude and kinematic corrections of migrated images for non-unitary imaging operators Antoine Guitton 1 ABSTRACT Obtaining true-amplitude migrated images remains a challenging problem. One possible solution to address it is iterative inversion. However, inversion is an expensive process that can be rather difficult and expensive to apply, especially with 3-D data. In this paper, I propose computing an image that is close to the least-squares inverse image by approximating the Hessian, thus avoiding the need for iterative inversion. The Hessian is approximated with non-stationary matching filters. These filters are estimated from two images: one is the migration result (m 1 ) and the second is the migration result of the remodeled data computed from m 1. Tests on the Marmousi dataset show that this filtering approach gives results similar to iterative least-squares inversion at a lower cost. In addition, because no regularization was used in the inversion process, the filtering method produces an image with fewer artifacts. Applying this method in the angle domain yields similar conclusions. INTRODUCTION It is well known that most of our mathematical operators for seismic processing are not unitary. This means that for any linear transformation L, L L = I where ( ) stands for the adjoint and I is the identity matrix. Having non-unitary operators often results from approximations we make when building those operators. For example, we might not take the irregular and finite acquisition geometry of seismic surveys into account. For migration, different approaches exist to correct for these defects of our operators. Bleistein (1987), based on earlier works from Beylkin (1985), derived an inverse operator for Kirchhoff migration. A similar path is followed by Thierry et al. (1999) with the addition of non-linear inversion with approximated Hessian. Alternatively, least-squares migration with regularization has proven efficient with incomplete surface data, e.g., Nemeth et al. (1999) and illumination problems due to complex structures, e.g., Prucha et al. (1999); Kuehl and Sacchi (2001). Hu et al. (2001) introduce a deconvolution operator that approximates the inverse od the Hessian in the least-squares estimate of the migrated image. However, this method assumes a v(z) medium which means that the deconvolution operators are horizontally invariant. Recently, Rickett (2001) proposed estimating weighting functions from reference 1 antoine@sep.stanford.edu 349
2 350 Guitton SEP 113 images to compensate illumination effects for finite-frequency depth migration. This method corrects for amplitude effects only and requires some smoothing that can be rather difficult to estimate. In this paper, I propose a new strategy for approximating the inverse of the Hessian. This approach aims to estimate a bank of non-stationary matching filters (Rickett et al., 2001) between two migrated images that theoretically embed the effects of the Hessian. This approach is implemented after migration and is relatively cheap to apply. It can be applied on images migrated at zero-offset or in the angle domain. I illustrate this method with the Marmousi dataset. I demonstrate that this approach can effectively recover amplitudes similar to the one obtained with least-squares inversion. In addition, these filtered images have less artifacts than the least-squares result without regularization and at a much reduced cost. THEORY In this section, I show how the least-squares estimate of a migrated image can be approximated using non-stationary matching filters. In terms of cost, this approach is comparable to one and a half iterations of a conjugate-gradient method (CG), the first half iteration being the migration. The cost of estimating the non-stationary filters is negligible compared to the total cost of migration. First, given seismic data d and a migration operator L, we seek a model m such that This goal can be rewritten in the following form Lm = d. (1) 0 r d = Lm d (2) and is called the fitting goal. For migration, a model styling goal (regularization) is necessary to compensate for irregular geometry artifacts and uneven illumination (Prucha et al., 1999). I omit this term in my derivations and focus on the data fitting goal only. By estimating m in a least-squares sense, we want to minimize the objective function The least-squares estimate of the model is given by f (m) = r d 2 = Lm d 2. (3) ˆm = (L L) 1 L d (4) where L L is the Hessian of the transformation. My goal in this paper is to approximate the effects of the Hessian L L using non-stationary matching filters. Approximating the Hessian In equation (4), I define L d as the migrated image m 1 after a first migration such that ˆm = (L L) 1 m 1. (5)
3 SEP 113 Approximating the Hessian 351 In equation (5), ˆm and L L are unknown. Since I am looking for an approximate of the Hessian, I need to find two known images that are related by the same expression as in equation (5). This can be easily achieved by remodeling the data from m 1 with L and remigrating them with L as follows: d 1 = Lm 1 (6) m 2 = L d 1 = L Lm 1. (7) Notice a similarity between equations (5) and (7) except that in equation (7), only L L is unknown. Notice that m 2 has a mathematical significance: it is a vector of the Krylov subspace for the model ˆm. Now, I assume that we can write the inverse Hessian as a linear operator B such that ˆm = Bm 1 (8) and m 1 = Bm 2. (9) Equation (9) can be approximated as a fitting goal for a matching filter estimation problem where B is the convolution matrix with a bank of non-stationary filters (Rickett et al., 2001). This choice is rather arbitrary but reflects the general idea that the Hessian is a transform operator between two similar images. My hope is not to perfectly represent the Hessian, but to improve the migrated image at a lower cost than least-squares migration. In addition in equations (8) and (9), the deconvolution process becomes a convolution, which makes it much more stable and easy to apply. Hence, I can rewrite equation (9) such that the matrix B becomes a vector and m 2 becomes a convolution matrix (Robinson and Treitel, 1980): The goal now is to minimize the residual m 1 = M 2 b. (10) 0 r m1 = m 1 M 2 b (11) in a least-squares sense. Because we have many unknown filter coefficients in b, I introduce a regularization term that penalizes differences between filters as follows: 0 r m1 = m 1 M 2 b 0 r b = Rb (12) where R is the Helix derivative (Claerbout, 1998). The objective function for equation (12) becomes f (b) = r m1 2 +ɛ 2 r b 2. (13) where ɛ is a constant. The least-squares inverse is thus given by Once ˆb is estimated, the final image is obtained by computing ˆb = (M 2 M 2 +ɛ 2 R R) 1 M 2 m 1. (14) ˆm = m 1 ˆb (15)
4 352 Guitton SEP 113 where ( ) is the convolution operator. Therefore, I propose computing first a migrated image m 1, then computing a migrated image m 2 (equation (7)), and finally estimating a bank of non-stationary matching filters b, e.g., equation (12). The final improved image is obtained by applying the matching filters to the first image m 1, e.g., equation (8). In the next section, I illustrate this idea with the Marmousi dataset. I show that an image similar to the least-squares migration image can be effectively obtained. MIGRATION RESULTS I illustrate the proposed algorithm with the Marmousi dataset. I use a prestack split-step migration method with one reference velocity. I demonstrate with zero-offset images and angle gathers that the approximation of the Hessian with adaptive filters gives a result comparable to least-squares migration with fewer artifacts. Zero-offset prestack migration results I first show in Figure 1 a comparison between the migration result of the Marmousi dataset (m 1 ) and the remodeled data (m 2 ). We notice that the migration of the remodeled data (Figure 1b) lowers the amplitudes in the upper part of the model. Therefore, we expect the filters to correct for this difference. Figure 2 displays few estimated filters for the Marmousi result. The filters are ten by ten with 40 patches in depth and 80 along the midpoint axis. I shown only a fifth of these filters in both axes. It is interesting to notice that these filters have their highest value at zero-lag, meaning that we have a strong amplitude correction with few kinematic changes. The zero-lag values are also larger at the top of the model, as anticipated. Looking more closely at these filters, we see that the coefficients follow the structure of the Marmousi model (upper right corner). Having estimated the filters b in equation (12), I apply them to m 1 to obtain an improved image. To validate this approach I show in Figure 3a the result of five conjugate gradient (CG) iterations with the Marmousi data. This results show higher amplitudes at the top but with inversion artifacts. This problem should be addressed with a proper regularization scheme (Prucha et al., 1999). In Figure 3b, I show the corrected image with the approximated Hessian B. The amplitude behavior is very similar to Figure 3a, without the inversion artifacts. Additionally, the cost is much lower. I show in Figure 4 the ratio of the envelope of Figure 3b and 1a. This Figure illustrates that the effects of the non-stationary filters, i.e, the Hessian, are stronger on the top of the model. Angle domain results Now, I apply the same method in the angle domain as opposed to the zero-offset image domain. The angle panels are created after migration from the offset panels directly (Stolt and Weglein,
5 SEP 113 Approximating the Hessian 353 Figure 1: (a) Migration result of the Marmousi dataset, i.e., m 1 in equation (9). (b) Migration result of the remodeled data, i.e., m 2 in equation (9). antoine3-m1.m2 [CR,M]
6 354 Guitton SEP 113 Figure 2: Each cell represents a non-stationary filter with its zero-lag coefficient in the middle. A fifth of the filters are actually shown in both directions. Each filter position corresponds roughly to a similar area in the model space (Figure 1a). After close inspection of the filter coefficients, these filters seem to follow the structure of the Marmousi model. They are also stronger at the top of the model, as expected. antoine3-filter-minv0 [CR]
7 SEP 113 Approximating the Hessian 355 Figure 3: (a) Model estimated after five iterations of CG. The model is noisy because no regularization has been applied. (b) Model estimated after applying the adaptive filters to m 1. The amplitude behavior is similar to (a) without the artifacts and with fewer iterations. antoine3-inv.ampcorr [CR,M]
8 356 Guitton SEP 113 Figure 4: Ratio of the envelopes of Figures 1a and 3b. Brighter colors correspond to higher values. The main effect of the filters is clearly visible at the top. antoine3-amplitude0 [CR] 1985; Weglein and Stolt, 1999; Sava and Fomel, 2000; Sava et al., 2001). I show in Figures 5 and 6 the angle gathers at two locations of the model (3050 and 5550 m) for the migration result (m1), the filtering result (ˆb m 1 ) and the inversion result after five CG iterations. Again, the angle gather after filtering has less artifacts than with inversion with a similar amplitude pattern throughout the section. In Figure 6, we notice that the filtering approach improves the continuity of some reflectors (between 500 and 600 m.) and that the amplitude increases for large angles. This behavior is similar to what Sava et al. (2001) observed for wave equation migration with amplitude corrections. I display in Figure 7 the estimated filters for the angle gather in Figure 5. Again, we see that we have higher amplitudes for filters at large aperture angles. In addition, we notice a smoothing effects of the filters along the angle axis. Finally I show in Figure 8 a comparison of stacked images across angle for m 1 and m 2. The bottom reflectors are stronger than in Figure 1 with fewer migration artifacts, as expected from the stack. The differences between the two stacked-data panels in Figure 8 are similar to what we observe in Figure 1, i.e, stronger amplitudes at the top of the image for m 1 and fewer migration artifacts for m 2. In Figure 9, I display a comparison between the inversion and filtering approaches for the stacked data. Again, the amplitude pattern is the same for both images with far less artifacts in the filtering result (Figure 9b).
9 SEP 113 Approximating the Hessian 357 DISCUSSION In this paper, I presented a method that intends to correct migrated images by approximating the Hessian of the imaging operator with non-stationary matching filters. These filters are estimated from two migration results. One migrated image, m 1, corresponds to the first migration result. The second image, m 2, is computed by re-modeling the data from m 1 and then by re-migrating it. It turns out that the relationship between m 1 and m 2 is similar to the relationship that exists between the least-squares inverse ˆm and m 1. In the proposed approach, this relationship is simply captured by matching filters. I demonstrate with the Marmousi dataset that this approach gives a better image than does least-squares without regularization and at a lower cost. In addition this approach can be used on images migrated at zero-offset or in the angle domain. As opposed to Hu et al. (2001), the correction in the image is completely data driven, does not depend on the velocity, and can be applied with any migration operator. It also works in the poststack or prestack domain without any extra effort. Providing the data and the ability to run at least two migrations to estimate m 2, this method would be easy to apply with 3-D migrated images. Compared to Rickett (2001), this proposed approach does not need reference images. In addition, the multi-dimensional filters offer more degrees of freedom for the correction than does a simple zero-lag weight: in that we might also correct for kinematic changes and move energy locally in the image. In the limit case where we choose one filter per point and only one coefficient (zero-lag) per filter, the matching filter approach would theoretically perform better than Rickett s method because the weights would be optimal in a least-squares sense without ad-hoc smoothing. In the future, it would be valuable to go beyond 2-D filters by extending them to 3-D and to test this 3-D approach with more field data. In addition, these filters could be used as preconditioning operators providing faster convergence for iterative inversion. The stability of the inverse filters still needs to be addressed, however. ACKNOWLEDGMENTS I would like to thank Paul Sava for providing the migration codes. REFERENCES Beylkin, G., 1985, Imaging of discontinuities in the inverse scattering problem by inversion of a causal generalized Radon transform: J. Math. Phys., 26, Bleistein, N., 1987, On the imaging of reflectors in the Earth: Geophysics, 52, no. 07, Claerbout, J., 1998, Multidimensional recursive filters via a helix: Geophysics, 63, no. 05,
10 358 Guitton SEP 113 Figure 5: The cmp location is 3050 m. (a) Migration result. (b) Inversion result. (c) Filtering result (ˆb m 1 ). The filtering result is comparable in amplitude to the inversion panel with fewer artifacts. antoine3-comp.ang.100 [CR,M] Hu, J., Schuster, G. T., and Valasek, P., 2001, Poststack migration deconvolution: Geophysics, 66, no. 3, Kuehl, H., and Sacchi, M., 2001, Generalized least-squares DSR migration using a common angle imaging condtion: 71st Ann. Internat. Mtg, Soc. Expl. Geophys., Expanded Abstracts, Nemeth, T., Wu, C., and Schuster, G. T., 1999, Least-squares migration of incomplete reflection data: Geophysics, 64, no. 1, Prucha, M., Biondi, B., and Symes, W., 1999, Angle-domain common image gathers by waveequation migration: 69th Ann. Internat. Mtg, Soc. Expl. Geophys., Expanded Abstracts, Rickett, J., Guitton, A., and Gratwick, D., 2001, Adaptive Multiple Subtraction with Non- Stationary Helical Shaping Filters: 63rd Mtg., Eur. Assn. Geosci. Eng., Expanded Abstracts, Session: P167. Rickett, J., 2001, Model-space vs. data-space normalization for finite-frequency depth migration: 71st Ann. Internat. Mtg, Soc. Expl. Geophys., Expanded Abstracts,
11 SEP 113 Approximating the Hessian 359 Figure 6: The cmp location is 5550 m. (a) Migration result. (b) Inversion result. (c) Filtering result (ˆb m 1 ). antoine3-comp.ang.300 [CR,M] Figure 7: Sample of the estimated matched filters for Figure 5a. Each filter is 10 by 10. The shape of these filters prove that they tend to smooth energy along the main reflectors. antoine3-filter-minve0 [CR]
12 360 Guitton SEP 113 Figure 8: (a) Stack across angle of the migrated image m 1. (b) Stack across angle of the migrated image m 2. antoine3-comp.stack.ang1 [CR,M]
13 SEP 113 Approximating the Hessian 361 Figure 9: (a) Stack across angle of the inverted image. (b) Stack across angle of the filtered image ˆb m 1. antoine3-comp.stack.ang2 [CR,M]
14 362 Guitton SEP 113 Robinson, E. A., and Treitel, S., 1980, Geophysical signal analysis: Prentice-Hall, Inc. Sava, P., and Fomel, S., 2000, Angle-gathers by Fourier Transform: SEP 103, Sava, P., Biondi, B., and Fomel, S., 2001, Amplitude-preserved common image gathers by wave-equation migration: 71st Ann. Internat. Mtg, Soc. Expl. Geophys., Expanded Abstracts, Stolt, R. H., and Weglein, A. B., 1985, Migration and inversion of seismic data: Geophysics, 50, no. 12, Thierry, P., Operto, S., and Lambare, G., 1999, Fast 2-D ray+born migration/inversion in complex media: Geophysics, 64, no. 1, Weglein, A. B., and Stolt, R. H., 1999, Migration-inversion revisited (1999): The Leading Edge, 18, no. 8,
Amplitude 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 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 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 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 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 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 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 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 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 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 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 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 informationRay based tomography using residual Stolt migration
Stanford Exploration Project, Report 11, November 11, 00, pages 1 15 Ray based tomography using residual Stolt migration Robert G. Clapp 1 ABSTRACT In complex areas, residual vertical movement is not an
More informationPrestack residual migration in the frequency domain
GEOPHYSICS, VOL. 68, NO. (MARCH APRIL 3); P. 634 64, 8 FIGS. 1.119/1.156733 Prestack residual migration in the frequency domain Paul C. Sava ABSTRACT Prestack Stolt residual migration can be applied to
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 informationResidual move-out analysis with 3-D angle-domain common-image gathers
Stanford Exploration Project, Report 115, May 22, 2004, pages 191 199 Residual move-out analysis with 3-D angle-domain common-image gathers Thomas Tisserant and Biondo Biondi 1 ABSTRACT We describe a method
More informationStanford Exploration Project, Report 103, April 27, 2000, pages
Stanford Exploration Project, Report 103, April 27, 2000, pages 205 324 204 Stanford Exploration Project, Report 103, April 27, 2000, pages 205 324 Data alignment with non-stationary shaping filters James
More informationCoherent noise attenuation using Inverse Problems and Prediction Error Filters
Stanford Exploration Project, Report 105, September 5, 2000, pages 1 49 Coherent noise attenuation using Inverse Problems and Prediction Error Filters Antoine Guitton 1 ABSTRACT Two iterative methods that
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 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 informationStanford Exploration Project, Report 120, May 3, 2005, pages
Stanford Exploration Project, Report 120, May 3, 2005, pages 167 179 166 Stanford Exploration Project, Report 120, May 3, 2005, pages 167 179 Non-linear estimation of vertical delays with a quasi-newton
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 information(t x) domain, pattern-based multiple separation
Stanford Exploration Project, Report 103, April 27, 2000, pages 63?? (t x) domain, pattern-based multiple separation Robert G. Clapp and Morgan Brown 1 ABSTRACT Pattern-based signal/noise separation is
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 informationInterpolation with pseudo-primaries: revisited
Stanford Exploration Project, Report 129, May 6, 2007, pages 85 94 Interpolation with pseudo-primaries: revisited William Curry ABSTRACT Large gaps may exist in marine data at near offsets. I generate
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 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 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 informationAttenuation of diffracted multiples with an apex-shifted tangent-squared radon transform in image space
Attenuation of diffracted multiples with an apex-shifted tangent-squared radon transform in image space Gabriel Alvarez, Biondo Biondi, and Antoine Guitton 1 ABSTRACT We propose to attenuate diffracted
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 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 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 informationStanford Exploration Project, Report 111, June 9, 2002, pages INTRODUCTION THEORY
Stanford Exploration Project, Report 111, June 9, 2002, pages 231 238 Short Note Speeding up wave equation migration Robert G. Clapp 1 INTRODUCTION Wave equation migration is gaining prominence over Kirchhoff
More informationCoherent partial stacking by offset continuation of 2-D prestack data
Stanford Exploration Project, Report 82, May 11, 2001, pages 1 124 Coherent partial stacking by offset continuation of 2-D prestack data Nizar Chemingui and Biondo Biondi 1 ABSTRACT Previously, we introduced
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 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 informationDealing with errors in automatic velocity analysis
Stanford Exploration Project, Report 112, November 11, 2002, pages 37 47 Dealing with errors in automatic velocity analysis Robert G. Clapp 1 ABSTRACT The lack of human interaction in automatic reflection
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 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 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 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 informationReconciling processing and inversion: Multiple attenuation prior to wave-equation inversion
Reconciling processing and inversion: Multiple attenuation prior to wave-equation inversion Claudio Guerra and Alejandro Valenciano ABSTRACT Seismic inversion is very sensitive to the presence of noise.
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 informationMissing trace interpolation and its enhancement of seismic processes
Missing trace interpolation Missing trace interpolation and its enhancement of seismic processes Wai-kin Chan and Robert R. Stewart ABSTRACT Many multi-channel seismic algorithms assume that the input
More informationSummary. Figure 1: Simplified CRS-based imaging workflow. This paper deals with the boxes highlighted in green.
Smoothing and automated picking of kinematic wavefield attributes Tilman Klüver and Jürgen Mann, Geophysical Institute, University of Karlsruhe, Germany Copyright 2005, SBGf Sociedade Brasiliera de Geofísica
More information3D pyramid interpolation
3D pyramid interpolation Xukai Shen ABSTRACT Seismic data geometries are not always as nice and regular as we want due to various acquisition constraints. In such cases, data interpolation becomes necessary.
More information3-D prestack migration of common-azimuth data
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 126 3-D prestack migration of common-azimuth data Biondo Biondi and Gopal Palacharla 1 ABSTRACT In principle, downward continuation of 3-D
More informationAmplitude-preserved low velocity noise suppression
Stanford Exploration Project, Report 84, May 9, 2001, pages 197 213 Amplitude-preserved low velocity noise suppression Hector Urdaneta 1 ABSTRACT Low velocity noise can contaminate reflections and distort
More informationSeismic data interpolation with the offset continuation equation
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 125?? Seismic data interpolation with the offset continuation equation Sergey Fomel 1 ABSTRACT I propose a finite-difference offset
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 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 informationWave-equation MVA applied to 4-D seismic monitoring
Stanford Exploration Project, Report 112, November 11, 2002, pages 15 21 Short Note Wave-equation MVA applied to 4-D seismic monitoring Paul Sava, John Etgen, and Leon Thomsen 1 INTRODUCTION 4-D seismic
More informationChapter 1. 2-D field tests INTRODUCTION AND SUMMARY
Chapter 1 2-D field tests INTRODUCTION AND SUMMARY The tomography method described in the preceding chapter is suited for a particular class of problem. Generating raypaths and picking reflectors requires
More informationA pattern-based technique for ground-roll and multiple attenuation
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? A pattern-based technique for ground-roll and multiple attenuation Antoine Guitton, Morgan Brown, James Rickett, and Robert Clapp 1 ABSTRACT
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 informationTransformation to dip-dependent Common Image Gathers
Stanford Exploration Project, Report 11, November 11, 00, pages 65 83 Transformation to dip-dependent Common Image Gathers Biondo Biondi and William Symes 1 ABSTRACT We introduce a new transform of offset-domain
More informationInterval velocity estimation through convex optimization
Stanford Exploration Project, Report 125, January 16, 2007, pages 1?? Interval velocity estimation through convex optimization Ben Witten and Michael Grant ABSTRACT Convex optimization is an optimization
More informationLeast-squares joint imaging of primaries and multiples
Stanford Exploration Project, Report 111, June 10, 2002, pages 17 47 Least-squares joint imaging of primaries and multiples Morgan Brown 1 ABSTRACT Multiple reflections provide redundant, and sometimes
More informationStanford Exploration Project, Report 100, April 20, 1999, pages
Stanford Exploration Project, Report 100, April 20, 1999, pages 197 210 196 Stanford Exploration Project, Report 100, April 20, 1999, pages 197 210 Directional smoothing of non-stationary filters Robert
More informationStanford Exploration Project, Report 124, April 4, 2006, pages 49 66
Stanford Exploration Project, Report 124, April 4, 2006, pages 49 66 48 Stanford Exploration Project, Report 124, April 4, 2006, pages 49 66 Mapping of specularly-reflected multiples to image space: An
More informationDSR Migration Velocity Analysis by Differential Semblance Optimization
DSR Migration Velocity Analysis by Differential Semblance Optimization A. Khoury (Total E&P France), W. W. Symes (Rice University), P. Williamson and P. Shen (Total E&P USA Inc.) Society of Exploration
More informationMultidimensional recursive filter preconditioning in geophysical estimation problems
GEOPHYSICS, VOL. 68, NO. 2 (MARCH-APRIL 2003); P. 577 588, 20 FIGS., 1 TABLE. 10.1190/1.1567228 Multidimensional recursive filter preconditioning in geophysical estimation problems Sergey Fomel and Jon
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 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 informationWide-azimuth angle gathers for wave-equation migration
1 Wide-azimuth angle gathers for wave-equation migration 2 3 Paul Sava (Center for Wave Phenomena, Colorado School of Mines) Ioan Vlad (Statoil) 4 5 6 7 (September 14, 2010) GEO-2010-???? Running head:
More informationConjugate guided gradient (CGG) method for robust inversion and its application to velocity-stack inversion a
Conjugate guided gradient (CGG) method for robust inversion and its application to velocity-stack inversion a a Published in Geophysics, 71, no. 4, R59-R67, (2006) Jun Ji ABSTRACT This paper proposes a
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 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 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 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 informationLeast-squares joint imaging of primaries and pegleg multiples: 2-D field data test
Stanford Exploration Project, Report 113, July 8, 2003, pages 17 30 Least-squares joint imaging of primaries and pegleg multiples: 2-D field data test Morgan Brown 1 ABSTRACT In this paper I present an
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 information96 Alkhalifah & Biondi
Stanford Exploration Project, Report 97, July 8, 1998, pages 95 116 The azimuth moveout operator for vertically inhomogeneous media Tariq Alkhalifah and Biondo L. Biondi 1 keywords: AMO, DMO, dip moveout
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 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 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 informationLeast-squares imaging and deconvolution using the hybrid norm conjugate-direction solver
Least-squares imaging and deconvolution using the hybrid norm conjugate-direction solver Yang Zhang and Jon Claerbout ABSTRACT To retrieve a sparse model, we applied the hybrid norm conjugate-direction
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 information=, (1) SEG/New Orleans 2006 Annual Meeting
U. Albertin* +, P. Sava ++, J. Etgen +, and M. Maharramov + + BP EPTG, Houston, Texas, ++ Colorado School of Mines, Goldin, Colorado Summary A methodology for velocity updating using one-way wavefield
More informationAmplitude-preserved multiple suppression
Stanford Exploration Project, Report 82, May 11, 2001, pages 1 45 Amplitude-preserved multiple suppression David Lumley, Dave Nichols, and Thorbjørn Rekdal 1 ABSTRACT Water-bottom and pegleg multiples
More informationMultiple attenuation in the image space
GEOPHYSICS, VOL. 70, NO. 1 (JANUARY-FEBRUARY 2005); P. V10 V20, 15 FIGS. 10.1190/1.1852789 Multiple attenuation in the image space Paul Sava 1 and Antoine Guitton 1 ABSTRACT Multiples can be suppressed
More informationSeismic reflection data interpolation with differential offset and shot continuation a
Seismic reflection data interpolation with differential offset and shot continuation a a Published in Geophysics, 68, 733-744 (2003) Sergey Fomel ABSTRACT I propose a finite-difference offset continuation
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 informationOffset plane waves vs. common-azimuth migration for sub-salt imaging
Stanford Exploration Project, Report, October 5, 999, pages?? Offset plane waves vs. common-azimuth migration for sub-salt imaging Biondo Biondi keywords: depth migration, wave-propagation, wave equation
More informationData interpolation in pyramid domain
Data interpolation in pyramid domain Xukai Shen ABSTRACT Pyramid domain is defined as a frequency-space domain with different spatial grids for different frequencies. Data interpolation in pyramid domain
More informationEstimation of multiple scattering by iterative inversion, Part II: Practical aspects and examples
GEOPHYSICS, VOL. 62, NO. 5 (SEPTEMBER-OCTOBER 1997); P. 1596 1611, 18 FIGS. Estimation of multiple scattering by iterative inversion, Part II: Practical aspects and examples D. J. Verschuur* and A. J.
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 informationData dependent parameterization and covariance calculation for inversion of focusing operators
Stanford Exploration Project, Report 11, September 18, 21, pages 1 91 Data dependent parameterization and covariance calculation for inversion of focusing operators Barbara E. Cox 1 ABSTRACT The Common
More informationData preparation for prestack depth migration
Data preparation for prestack depth migration Hugh D. Geiger Data preparation for prestack migration SUMMARY The goal of prestack depth migration is to reconstruct an image of the subsurface with the highest
More informationAmplitude-preserved wave-equation migration
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? Amplitude-preserved wave-equation migration Paul Sava and Biondo Biondi 1 ABSTRACT We analyze the amplitude variation as a function of
More informationPlane-wave migration in tilted coordinates
Plane-wave migration in tilted coordinates Guojian Shan and Biondo Biondi ABSTRACT Most existing one-way wave-equation migration algorithms have difficulty in imaging steep dips in a medium with strong
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 informationIllumination compensation with Poynting vectors
Illumination compensation using Poynting vectors, with special treatment for multiples Alan Richardson and Alison E. Malcolm, Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute
More informationAutomatic wave equation migration velocity analysis by differential semblance optimization
Automatic wave equation migration velocity analysis by differential semblance optimization Peng Shen, Christiaan Stolk, William W. Symes The Rice Inversion Project Rice University Thanks to Dr. Scott Morton
More informationRay-based tomography with limited picking
Stanford Exploration Project, Report 110, September 18, 2001, pages 1 113 Short Note Ray-based tomography with limited picking Robert G. Clapp 1 INTRODUCTION In ray-based reflection tomography picking
More informationChapter 3. 2-D Field Data Results
Chapter 3 2-D Field Data Results In 1997, WesternGeco distributed a 2-D dataset, acquired in the Mississippi Canyon region of the Gulf of Mexico, for the testing of multiple suppression algorithms. As
More informationSearching the Sea of Galilee The splendors and miseries of iteratively reweighted least squares
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 267 Searching the Sea of Galilee The splendors and miseries of iteratively reweighted least squares Sergey Fomel and Jon F. Claerbout 1 As
More information3-D prestack depth migration of common-azimuth data
Stanford Exploration Project, Report 84, May 9, 200, pages 25 3-D prestack depth migration of common-azimuth data Biondo Biondi and Gopal Palacharla ABSTRACT We extended the common-azimuth 3-D prestack
More informationEstimating a 2D stationary PEF on sparse data
Stanford Exploration Project, Report 117, October 23, 2004, pages 129 140 Estimating a 2D stationary PEF on sparse data Jesse Lomask 1 ABSTRACT A stationary 2D PEF and missing data are simultaneously estimated
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 information