Split-step complex Padé-Fourier depth migration
|
|
- Joy Howard
- 6 years ago
- Views:
Transcription
1 Geophys. J. Int. (007) 7, doi: 0./j X x Split-step complex Padé-Fourier depth migration Linbin Zhang,,, James W. Rector III, 3 G. Michael Hoversten and Sergey Fomel 4 Department of Material Science and Engineering, University of California, Berkeley, CA 9470, USA. linbin.zhang@total.com One Cyclotron Road, Lawrence Berkeley National Laboratory, MS 90-6, Berkeley, CA 9470, USA 3 Department of Civil and Environmental Engineering, University of California, Berkeley, CA 9470, USA 4 Bureau of Economic Geology, The University of Texas at Austin, University Station, Box X, Austin, TX , USA Accepted 007 September 7. Received 007 August 7; in original form 006 November 8 SUMMARY We present a split-step complex Padé-Fourier migration method based on the one-way wave equation. The downward-continuation operator is split into two downward-continuation operators: one operator is a phase-shift operator and the other operator is a finite-difference operator. A complex treatment of the propagation operator is applied to mitigate inaccuracies and instabilities due to evanescent waves. It produces high-quality images of complex structures with fewer numerical artefacts than those obtained using a real approximation of a square-root operator in the one-way wave equation. Tests on zero-offset data from the SEG/EAGE salt data show that the method improves the image quality at the cost of an additional 0 per cent computational time compared to the conventional Fourier finite-difference method. Key words: depth migration, Fourier finite-difference method, Padé approximation. GJI Seismology INTRODUCTION Wave equation migration is a rapidly growing tool for complex structure imaging (Claerbout 985; Stoffa et al. 990; Lee et al. 99; Ristow & Ruhl 994; Wu & Jin 997; Bonomi & Cazzola 999). Wave equation migration methods can produce accurate images provided that the background velocity field is well defined. Most widely used wave equation methods are based on solving the one-way acoustic wave equation. There are many different methods to numerically solve the one-way equation. They can be classified into three main approaches: Fourier methods (Gazdag & Sguazzero 984; Stoffa et al. 990) solved in the wavenumber domain and space domain; finite-difference methods (FD) (Claerbout 985; Hale 99) and mixed methods that are a linear combination of spectral and FD methods (Ristow & Rulh 994; Biondi 00). The Fourier methods, such as Split-Step Fourier method (SSF) (Stoffa et al. 990) and Phase Shift Plus Interpolation (PSPI) (Gazdag & Sguazzero 984; Bagani et al. 995) can handle both wide-angle propagation and lateral velocity variations, but the imaging becomes less accurate when the lateral velocity variations are strong. The FD methods can easily handle strong lateral velocity variations. When implicit extrapolators are used in FD methods, the one-way wave equation is solved in the frequency-space domain, and the square root of the wave operator is approximated using continued fractions expansions (Claerbout 985). The essential shortcoming Now at: Total E&P USA Inc., 800 Gessner, Suite 700, Houston, TX 7704, USA. of this method is that it is not efficient for wide-angle propagation. High-order approximations to the square root operators are required for imaging of steep dips. When explicit extrapolators are used (Holberg 988; Hale 99), the exponential operator is approximated by finite length spatial filters. The disadvantage of explicit methods is that their stability is not guaranteed if the filter is not designed carefully. In 994, Ristow and Ruhl proposed a Fourier finite-difference method (FFD) which combined a phase-shift method in the frequency-wavenumber domain and an FD method in the frequencyspace domain. The FFD method combines the advantages of the FD methods and the phase shift methods and is more accurate than SSF, PSPI and FD methods. Huang & Fehler (00) proposed a SSF Padé migration method that allows for the use of a larger grid spacing than conventional FD migration methods. Huang and Fehler s method is equivalent to applying the Crank Nicholson scheme to the FD part in the FFD method of Ristow and Ruhl. In both FD and FFD migration methods, the real Padé approximation or a Taylor expansion is usually used to approximate the square-root operator. For example, the real [n/n] Padé approximation to the square-root operator + X is given by: + X + n j= a j X + b j X, where X = V (x,z), and a x j and b j are the real-valued coefficients. If X < (evanescent waves), the right-hand side of the above equation is real and the left-hand side is complex. The real approximation of the square-root operators is inconsistent and will cause accuracy and stability problems. It cannot handle evanescent ω 308 C 007 The Authors
2 Split-step complex Padé-Fourier depth migration 309 waves properly (Milinazzo et al. 997; Yevick & Thomson 000). In FFD method, usually the minimum velocity is used as the reference velocity. In this case, most evanescent wave energy still exists in the wavefield. If the real complex Padé approximation is used, we will have inaccuracy and instability problem. To deal with this problem, a complex Padé approximant technique (Milinazzo et al. 997) can be employed. Before the Padé approximations, one can rewrite + X as: + X = e iθ/ + [(X + )e iθ ] = e iθ/ + X, where X = (X + )e iθ, and then apply the Padé approximation to the operator + X. This procedure maps the poles of the Padé approximant along the real axis for X < into the positive imaginary half-space. This method improved the stability significantly (Milinazzo et al. 997; Lingevitch & Collins 998; Yevick & Thomson 000). Zhang et al. (003) have applied this technique to the FD migration method. However, the FD method is not efficient for wide-angle propagation because using a [/] Padé approximant restricts the maximum migration angle to about 45. In this paper, we apply complex Padé approximant technique to FFD method to mitigate the inaccuracies and instabilities caused by evanescent waves. In the first part of the paper, we derive a split-step complex Padé-Fourier solution for the one-way wave equation using a Padé approximant. Then we present zero-offset migration results for the SEG/EAGE salt model. A numerical comparison of FFD, split-step complex Padé migration and our method highlights the advantages of the proposed method. METHOD The -D, one-way acoustic wave equation can be expressed in frequency domain as: P ω z = i V + x P = i ω + XP, () V where P = P(z, x, ω) is the wavefield, i =, and ω is the circular frequency. X = V ω x, () and V = V (z, x) is the velocity. The solution for eq. () is extrapolating in depth: P(z + z, x,ω) = e iδ +X P(z, x,ω), (3) where δ = ω z and z is the spacing in the z-direction. V Let δ 0 = ω V ref z and X 0 = V ref, where V ω x ref denotes the minimum velocity at the layer [z, z + z]. Multiplying and dividing by the approximate exponential operator e iδ 0 +X0, the exact exponential operator e iδ +X can be written as (Ristow & Ruhl 994; Huang & Fehler 00): e iδ +X = e iδ 0 +X0 +iδ 0 ( +X m + X m ), (4) where m = V V ref. Finally, we express the wave field depth extrapolation equation as (see Appendix A for the derivation) : P(z + z, x,ω) e iδ 0 +X0 e iδ 0( m ) ( ) a + b X P(z, x,ω). a + b X (5) The first two exponential terms are implemented using SSF algorithm (Stoffa et al. 990) and the third term is implemented by FD method. When θ = 0, eq. (5) is identical to FFD method with the Crank Nicholson scheme applied to solve the third term of eq. (5) (Huang & Felher 00). Therefore, the maximum migration angle of eq. (5) is the same as that of Huang and Felher s and FFD method. However, eq. (5) can properly handle the propagation of evanescent waves. To obtain more accurate results with FD method, we employ the Douglas operator to approximate the second derivative (Tsuchiya et al. 00): x x δ x x, (6) + /δx where δ x P(x) = P(x + x) P(x) + P(x x), x is the spacing in x-direction. One can also apply [/] Padé approximation to the exponential function E = e iδ 0 A to obtain: E + iδ 0 A + (iδ 0 A) iδ 0 A +. (7) (iδ 0 A) Using eq. (7), eq. (5) becomes: P(z + z, x,ω) e iδ 0 +X0 e iδ 0( m ) + iδ 0 A + (iδ 0 A) iδ 0 A + P(z, x,ω), (8) (iδ 0 A) Using the multistep method (Hadley 99) to solve the equation: P(z + z, x,ω) = + iδ 0 A + (iδ 0 A) iδ 0 A + P(z, x,ω), (9) (iδ 0 A) we get the following two-step algorithm: P(z + k z/, x,ω) = + a k(iδ 0 A) P(z + (k ) z/, x,ω), for k =, (0) a k (iδ 0 A) where a k = (3 + ( ) k 3i)/. The truncation error of [/] Padé approximant of eq. (5) is O( z ) while the overall truncation error of eq. (0) is O( z 4 ). So eq. (0) allows us to use a larger grid spacing z than that of eq. (5). NUMERICAL RESULTS We test our method with post-stack data for the SEG/EAGE salt model. The velocity model is shown in Fig.. The data set consists of 90 traces, each trace contains 66 time samples and a sampling interval of 8 ms (Fig. ). The velocity model grid is with x = z =.9 m. No additional filters were applied to the migration results. Figs 3(a) and (b) show the images using eq. (5) with θ = 0 (FFD) and θ = 5, respectively. Fig. 3(b) gives an image with fewer evanescent-energy artefacts than those produced by the FFD migration (Fig. 3a) while the accuracy of both migrations is similar to each other. The computational time of the Split-step complex Padé-Fourier migration increases 0 per cent compared to the FFD method and the memory usage is the same for both methods. Fig 3(c) shows the image using complex Padé migration method for the exploding reflector data set. Comparing Figs 3(b) and (c) we see that our method gives better images of the top and base interfaces of the salt body (Location A and B in those figures). Moreover, the image in Fig. 3(c) has fewer numerical artefacts than Fig. 3(a). Figs 4(a) (c) show the images with θ = 0,0 and 30, respectively. These images show that the artefacts decrease as θ increases, C 007 The Authors, GJI, 7,
3 30 L. Zhang et al. Figure. SEG/EAGE -D velocity model. Figure. Exploding reflector modelled data (zero offset equivalent). but the events become weaker, for example at the location pointed to by the arrow. Thus it appears that there is a trade-off in image quality and θ. Experiments show that θ should be less than 5. Figs 5 and 6 show the images using eq. (0) with θ = 5 and vertical grid spacing.9 and m, respectively. Comparing Fig. 5 with Figs 3(a) and (b), the [/] Padé approximant yields improved image quality. The reflectors under the salt body are much clearer than those in Fig. 3. However it took double CUP times of that using eq. (8). With a coarse spacing, we can also obtain a good image (Fig. 6) but it is not as good as Fig. 5. This is because the true velocity model cannot be represented well by the coarse spacing of m. Therefore, we lacked necessary resolution in migration. CONCLUSIONS We have developed a split-step complex Padé-Fourier technique for seismic migration. Our method makes use of the advantages of both FD methods and spectral methods. Because the complex treatment of the propagation operator is used, it treats the evanescent waves with more accuracy, which improves the stability of the Figure 3. Comparison of migration images of an exploding-reflector data set for the -D model shown in Fig.. (a) Split-step complex Padé-Fourier migration with θ = 0 (FFD). (b) Split-step complex Padé-Fourier migration with θ = 5. (c) Split-step complex Padé migration with θ = 5. The images shown in (b) and (c) have less artefacts than that in (a). Both images in (a) and (b) give better images of the top and base interfaces of the salt body than that in (c) (Locations A and B in those figures). C 007 The Authors, GJI, 7,
4 Split-step complex Padé-Fourier depth migration 3 Figure 5. Split-step complex Padé-Fourier migration (eq. 8) with θ = 5 and vertical grid spacing.9 m. Figure 6. Split-step complex Padé-Fourier migration (eq. 8) with θ = 5 and vertical grid spacing m. numerical approximation. Synthetic data migration results show that the proposed method can produce higher quality images. ACKNOWLEDGMENT This work was supported by DOE. We are grateful to the editors and two reviewers for their useful comments and suggestions. REFERENCES Figure 4. Images with different θ: (a) θ = 0 ; (b) θ = 0 and (c) θ = 30. As θ increases the event that the arrow pointed to is weaker but the image contains less noise. Bagani, C., Bonomi, E. & Pieroni, E., 995. Data parallel implementation of 3D PSPI, in 65th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, pp Biondi, B., 00. Stable wide-angle Fourier finite-difference downward extrapolation of 3-D wavefields, Geophysics, 67, Bonomi, E. & Cazzola, L., 999. Prestack imaging of compressed seismic data: a Monte Carlo approach, in 69th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, pp Claerbout, J. F., 985. Imaging the Earth s interior, Blackwell Scientific Publications, Inc. C 007 The Authors, GJI, 7,
5 3 L. Zhang et al. Gazdag, J. & Sguazzero, P., 984. Migration of seismic data by phase shift plus interpolation, Geophysics, 49, 4 3. Hale, D., 99. Stable explicit depth extrapolation of seismic wavefields, Geophysics, 56, Hadley, G.R., 99. Multistep method for wide angle beam prpagation, Opt. Lett., 7, Holberg, O., 988. Towards optimum one-way wave propagation, Geophys. Prospect., 36, Huang, L. & Fehler C.F., 00. Split-step Fourier Padé migration, in 7th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, pp Lee, M., Mason, L.M. & Jackson, G.M., 99. Split-step Fourier shotrecord migration with deconvolution imaging, Geophysics, 56, Lingevith, J.F. & Collins, M.D., 998. Wave propagation in range-dependent poro-acoustic waveguides, J. Acoust. Soc. Am., 04, Milinazzo, F.A., Zala, C.A. & Brooke, G.H., 997. Rational square-root approximations for parabolic equation algorithms, J. Acoust. Soc. Am., 0, Ristow, D. & Ruhl, T., 994. Fourier finite-difference migration, Geophysics, 59, Stoffa, P.L., Fokkema, J.T, Freir, R.M. & Kessinger, W.P., 990. Split-step Fourier migration, Geophysics, 55, Tsuchiya, T., Anada, T., Endoh, N., Nakamura, T., Aoki, T. & Kaihou, I., 00. An efficient method combined the Douglas operator scheme tosplitstep Pade approximation of higher-order parabolic equation, in Ultrasonics Symposium, 00 IEEE, pp Wu, R. & Jin, S., 997. Windowed GSP(Generalized Screen Propagators) migration applied to SEG-EAGE salt model data, in 67th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, pp Yevick, D. & Thomson, 000. Complex Padé approximants for wide-angle acoustic propagators, J. Acoust. Soc. Am., 08, Zhang, L., Rector, J.W. & Hoversten, G.M, 003. Split-step complex Padé migration, J. Seismic Explor.,, APPENDIX A To handle evanescent waves, we first rotate the square-root operators in eq. (4) by θ (Milinazzo et al. 997), + X = e iθ/ [( + X)e iθ ] + == e iθ/ + X, (A) [( ) ] + Xm = eiθ/ + Xm e iθ + = e iθ/ + X, (A) where X = ( + X)e iθ and X = ( + X m )e iθ. Eq. (4) becomes: e iδ +X = e iδ 0 +X0 +iδ 0( m )+iδ 0[ m ( +X ) ( +X )], (A3) ω V ref ze iθ/. Let s set E = e iδ 0 [ m ( +X ) ( +X )] = e iδ 0 A, where A = m ( + X ) ( + X ). Then we apply where δ 0 = Taylor expansion to A, A ( ) X m X 8 ( ) X m X + = ( ) 8 (e iθ ) m (5 e iθ ) + ( ) e iθ (3 e iθ )X 4m m e iθ 8m ( m 3 ) X +, using the approximation: ax ax( bx +...). + bx We obtain: A 8 (e iθ ) ( ) m (5 e iθ ) + ( 4m Applying Padé approximation to e iδ 0 A,wehave: E + iδ 0 A iδ 0 A = a + b X a + b X, m ) e iθ (3 e iθ )X + (m +m+) m e iθ X. (A4) (3 e iθ ) C 007 The Authors, GJI, 7,
6 Split-step complex Padé-Fourier depth migration 33 where ( ) a = + iδ 0 8 (e iθ ) m (5 e iθ ), b = (m + m + ) m { e iθ ( m 3 (3 e iθ ) + ) (e iθ )(5 e iθ )e iθ iδ 0 6m 3 (3 e iθ ) ( ) a = iδ 0 8 (e iθ ) m (5 e iθ ), b = (m + m + ) m { e iθ ( m 3 (3 e iθ ) ) (e iθ )(5 e iθ )e iθ iδ 0 6m 3 (3 e iθ ) + m } 4m e iθ (3 e iθ ), + m } 4m e iθ (3 e iθ ). (A5) C 007 The Authors, GJI, 7,
= 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 information3-D parallel shot-gather prestack depth migration
3-D parallel shot-gather prestack depth migration Wensheng Zhang 1, Guanquan Zhang 1 and Xingfu Cui 2 1 LSSC, Institute of Computational Mathematics and Scientific/Engineering Computing, Academy of Mathematics
More informationModeling 32D scalar waves using the Fourier finite2difference method
50 6 007 11 CHINESE JOURNAL OF GEOPHYSICS Vol. 50, No. 6 Nov., 007,,..,007,50 (6) :1854 186 Zhang J H,Wang W M, Zhao L F,et al. Modeling 3D scalar waves using the Fourier finitedifference method. Chinese
More informationSEG/New Orleans 2006 Annual Meeting
Accuracy improvement for super-wide angle one-way waves by wavefront reconstruction Ru-Shan Wu* and Xiaofeng Jia, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary To
More informationSEG/San Antonio 2007 Annual Meeting
Imaging steep salt flanks by super-wide angle one-way method Xiaofeng Jia* and Ru-Shan Wu, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary The super-wide angle one-way
More 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 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 informationChapter 1 Downward-continuation methods
Chapter 1 Downward-continuation methods The numerical solution of the one-way wave equation is the cornerstone of all downward-continuation migration methods. Over the years a wide variety of solutions
More information3D PRESTACK DEPTH MIGRATION WITH FACTORIZATION FOUR-WAY SPLITTING SCHEME
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Supp, Pages 183 196 c 5 Institute for Scientific Computing and Information 3D PRESTACK DEPTH MIGRATION WITH FACTORIZATION FOUR-WAY SPLITTING
More 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 informationWAVE EQUATION DEPTH MIGRATION USING COMPLEX PADÉ APPROXIMATION
Revista Brasileira de Geofísica (2009) 27(1): 95-102 2009 Sociedade Brasileira de Geofísica ISSN 0102-261X www.scielo.br/rbg WAVE EQUATION DEPTH MIGRATION USING COMPLEX PADÉ APPROXIMATION Reynam C. Pestana
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 informationMigration from a non-flat datum via reverse-time extrapolation
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 50 Migration from a non-flat datum via reverse-time extrapolation Gopal Palacharla 1 ABSTRACT Land surveys usually have elevation changes,
More informationThe Conventional and Non-conventional Seismic Differencing
Summary The Conventional and Non-conventional Seismic Differencing Vanja Vracar* University of Calgary, CREWES Project, Alberta, Canada vmilicev@ucalgary.ca and Robert J. Ferguson University of Calgary,
More informationDownward continuation with explicit extrapolators and common-azimuth migration
Downward continuation with explicit extrapolators and common-azimuth migration Eric Dussaud ABSTRACT We present a review of explicit extrapolation schemes in the context of 3-D poststack migration, and
More informationOffset-domain pseudoscreen prestack depth migration
GEOPHYSICS, VOL. 67, NO. 6 (NOVEMBER-DECEMBER 2002); P. 895 902, 7 FIGS. 0.90/.527089 Offset-domain pseudoscreen prestack depth migration Shengwen Jin, Charles C. Mosher, and Ru-Shan Wu ABSTRACT The double
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 information3D plane-wave migration in tilted coordinates
Stanford Exploration Project, Report 129, May 6, 2007, pages 1 18 3D plane-wave migration in tilted coordinates Guojian Shan, Robert Clapp, and Biondo Biondi ABSTRACT We develop 3D plane-wave migration
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 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 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 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 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 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 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 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 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 informationGEOPHYS 242: Near Surface Geophysical Imaging. Class 5: Refraction Migration Methods Wed, April 13, 2011
GEOPHYS 242: Near Surface Geophysical Imaging Class 5: Refraction Migration Methods Wed, April 13, 2011 Migration versus tomography Refraction traveltime and wavefield migration The theory of interferometry
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 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 informationRiemannian Wavefield Extrapolation: Non-orthogonal coordinate systems
Riemannian Wavefield Extrapolation: Non-orthogonal coordinate systems Jeff Shragge Stanford Exploration Project, Mitchell Bldg., Department of Geophysics, Stanford University, Stanford, CA 94305-2215 (June
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 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 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 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 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 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 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 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 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 informationRiemannian Wavefield Extrapolation: Non-orthogonal. coordinate systems
Riemannian Wavefield Extrapolation: Non-orthogonal coordinate systems Jeff Shragge (October 1, 27) Running head: ABSTRACT Riemanninan wavefield extrapolation (RWE) is used to model one-way wave propagation
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 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 informationMIGRATION BY EXTRAPOLATION OF TIME- DEPENDENT BOUNDARY VALUES*
Geophysical Prospecting 31,413-420, 1983. MIGRATION BY EXTRAPOLATION OF TIME- DEPENDENT BOUNDARY VALUES* G.A. McMECHAN** ABSTRACT MCMECHAN, G.A. 1983, Migration by Extrapolation of Time-Dependent Boundary
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 informationA full-wave equation based seismic illumination analysis method
A full-wave equation based seismic illumination analysis method Xiao-Bi Xie and Hui Yang University of California, Santa Cruz, CA 95064, USA Summary We propose a full-wave equation based method for seismic
More informationCrosswell Imaging by 2-D Prestack Wavepath Migration
Crosswell Imaging by 2-D Prestack Wavepath Migration Hongchuan Sun ABSTRACT Prestack wavepath migration (WM) is applied to 2-D synthetic crosswell data, and the migrated images are compared to those from
More informationWhich Depth Imaging Method Should You Use? A Roadmap In The Maze of 3 D Depth Imaging
Which Depth Imaging Method Should You Use? A Roadmap In The Maze of 3 D Depth Imaging By, Dimitri Bevc, 3DGeo Development Inc., 4633 Old Ironsides Drive, Suite 401, Santa Clara, CA 95054, dimitri@3dgeo.com
More informationAzimuth Moveout Transformation some promising applications from western Canada
Azimuth Moveout Transformation some promising applications from western Canada Satinder Chopra and Dan Negut Arcis Corporation, Calgary, Canada Summary Azimuth moveout (AMO) is a partial migration operator
More 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 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 informationFast Gabor Imaging with Spatial Resampling
Fast Gabor Imaging with Spatial Resampling Yongwang Ma* University of Calgary, Calgary, AB yongma@ucalgary.ca and Gary Margrave and Chad Hogan University of Calgary, Calgary, AB, Canada Summary The Gabor
More informationMany-core and PSPI: Mixing fine-grain and coarse-grain parallelism
Many-core and PSPI: Mixing fine-grain and coarse-grain parallelism Ching-Bih Liaw and Robert G. Clapp ABSTRACT Many of today s computer architectures are supported for fine-grain, rather than coarse-grain,
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 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 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 informationSeismic reflection data interpolation with differential offset and shot continuation
GEOPHYSICS, VOL. 68, NO. 2 (MARCH-APRIL 2003; P. 733 744, 12 FIGS. 10.1190/1.1567243 Seismic reflection data interpolation with differential offset and shot continuation Sergey Fomel ABSTRACT I propose
More informationShot-profile migration of GPR data
Shot-profile migration of GPR data Jeff Shragge, James Irving, and Brad Artman 1 ABSTRACT Multi-offset ground-penetrating radar data possess a number of important advantages over constant-offset data,
More informationWEI: Wave-Equation Imaging Library
Stanford Exploration Project, Report 111, June 9, 2002, pages 381 392 WEI: Wave-Equation Imaging Library Paul Sava and Robert G. Clapp 1 ABSTRACT This paper introduces WEI, a new library in the SEP library
More informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Aliasing Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 16, 2014 Bee Bednar (Panorama Technologies) Seismic
More informationGEOPHYSICS. Time to depth conversion and seismic velocity estimation using time migration velocities
Time to depth conversion and seismic velocity estimation using time migration velocities Journal: Manuscript ID: Manuscript Type: Date Submitted by the Author: Geophysics draft Letters n/a Complete List
More 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 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 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 informationReverse time migration with random boundaries
Reverse time migration with random boundaries Robert G. Clapp ABSTRACT Reading wavefield checkpoints from disk is quickly becoming the bottleneck in Reverse Time Migration. We eliminate the need to write
More information13 Recursive prestack depth migration using CFP-gathers 1
13 Recursive prestack depth migration using CFP-gathers 1 Jan Thorbecke 2 A.J. Berkhout 3 The Common Focus-Point Technology (CFP) describes prestack migration by two focusing steps: focusing in emission
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 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 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 informationDownloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at
Recovering the Reflectivity Matrix and Angle-dependent Plane-wave Reflection Coefficients from Imaging of Multiples Alba Ordoñez PGS/UiO*, Walter Söllner PGS, Tilman Klüver PGS and Leiv J. Gelius UiO Summary
More 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 informationShort Note. DMO velocity analysis with Jacubowicz s dip-decomposition method. David Kessler and Wai-Kin Chan*
GEOPHYSICS, VOL. 58, NO. 10 (OCTOBER 1993); P. 1517-1524,9 FIGS. Short Note DMO velocity analysis with Jacubowicz s dip-decomposition method David Kessler and Wai-Kin Chan* INTRODUCTION Dip-moveout (DMO)
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 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 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 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 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 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 information3-D traveltime computation using Huygens wavefront tracing
GEOPHYSICS, VOL. 66, NO. 3 MAY-JUNE 2001); P. 883 889, 12 FIGS., 1 TABLE. 3-D traveltime computation using Huygens wavefront tracing Paul Sava and Sergey Fomel ABSTRACT Traveltime computation is widely
More 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 informationConverted wave dip moveout
Stanford Exploration Project, Report 111, June 9, 2002, pages 47 59 Converted wave dip moveout Daniel Rosales 1 ABSTRACT Dip moveout (DMO) introduces a dip-dependent correction for a more appropiate transformation
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 informationAFRL-HE-BR-SR A STUDY OF INVERSE METHODS FOR PROCESSING OF RADAR DATA
AFRL-HE-BR-SR-2006-0010 A STUDY OF INVERSE METHODS FOR PROCESSING OF RADAR DATA Mohamed F. Chouikha, Ph.D. Howard University Department of Electrical and Computer Engineering 2400 Sixth Street NW Washington,
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 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 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 informationEfficient traveltime compression for 3D prestack Kirchhoff migration
Geophysical Prospecting, 2011, 59, 1 9 doi: 10.1111/j.1365-2478.2010.0886.x Efficient traveltime compression for 3D prestack Kirchhoff migration Tariq Alkhalifah The Institute for Geophysical Research,
More informationSEG Houston 2009 International Exposition and Annual Meeting
Yueming Ye *1, Ru-Shan Wu and Zhenchun Li 2 Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz, CA 95064 Summary Migration with data acquired on surface with irregular topography
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 informationStratigraphic coordinates, a coordinate system tailored to seismic interpretation a
Stratigraphic coordinates, a coordinate system tailored to seismic interpretation a a Published in Geophysical Prospecting, v. 63, 1246-1255, (2015) Parvaneh Karimi and Sergey Fomel ABSTRACT In certain
More informationSummary. Introduction
Yaofeng He and Ru-Shan Wu Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, C 9564 Summary We propose an approach to perform migration and imaging using secondary scattered
More informationPyramid-shaped grid for elastic wave propagation Feng Chen * and Sheng Xu, CGGVeritas
Feng Chen * and Sheng Xu, CGGVeritas Summary Elastic wave propagation is elemental to wave-equationbased migration and modeling. Conventional simulation of wave propagation is done on a grid of regular
More information3-D depth migration via McClellan transformations
GEOPHYSICS, VOL. 56. NO. I I (NOVEMBER 1991): P. 177X-17X5. I? FIGS. 3-D depth migration via McClellan transformations Dave Hale* ABSTRACT Three-dimensional seismic wavefields may be extrapolated in depth,
More informationDownloaded 05/09/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Elastic converted-wave path migration for subsalt imaging Ru-Shan Wu*, Rui Yan, Xiao-Bi Xie, Modeling and Imaging Laboratory, Earth and Planetary Sciences/IGPP, University of California, Santa Cruz, David
More 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 informationRecovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples
Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples A. Ordoñez* (PGS), W.F. Sollner (PGS), T. Klüver (PGS) & L.G. Gelius (UiO) SUMMARY A joint
More informationTheory of 3-D angle gathers in wave-equation seismic imaging
J Petrol Explor Prod Technol (2011) 1:11 16 DOI 10.1007/s13202-011-000-8 ORIGINAL PAPER - EXPLORATION GEOPHYSICS Theory of 3-D angle gathers in wave-equation seismic imaging Sergey Fomel Received: 26 September
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 informationk are within the range 0 k < k < π dx,
Acoustic Wave Equation by Implicit Spatial Derivative Operators Dan Kosloff*, Department of Geophysics, Tel-Aviv University and Paradigm, Reynam Pestana, Department of Geophysics, Federal University of
More information