Complex-beam Migration and Land Depth Tianfei Zhu CGGVeritas, Calgary, Alberta, Canada
|
|
- Theodore Fox
- 5 years ago
- Views:
Transcription
1 Page 1 of 10 Home Articles Interviews Print Editions Complex-beam Migration and Land Depth Tianfei Zhu CGGVeritas, Calgary, Alberta, Canada DECEMBER 2012 FOCUS ARTICLE Summary Gaussian-beam depth migration is a useful alternative to Kirchhoff and wave-equation migrations in imaging geologically complex land areas. It overcomes the limitations of Kirchhoff migration in imaging multipathing arrivals and suffers no dip limits from oneway wave equation approximation. The accuracy of the Gaussian-beam method in solving the wave equation is, however, parameter dependent, limiting its applications in areas with strong near-surface structural variations. We describe in this study a complex-beam method which represents an improvement over the Gaussian-beam migration: It retains the advantages of the Gaussian-beam method over Kirchhoff and wave-equation migration, and is more robust in imaging complex structures. Applications of this method to field and synthetic data show that the new method produces depth images superior to those from Kirchhoff migration. It compares well to the reverse-time migration in imaging steep to overturned structures. Introduction Depth imaging of geologically complex land areas is challenging because of their steeply folded and faulted structures and because of their large lateral variations in topography and near-surface velocity. These difficulties are often further compounded by irregular data-acquisition geometries found in the seismic surveys from these areas. Prestack Gaussianbeam migration provides a useful alternative to Kirchhoff and wave-equation migrations for imaging these areas: It overcomes the limitations of Kirchhoff migration in imaging multipathing arrivals while retaining its flexibility with input geometry and topographic variations (Hill, 2001; Nowack et al., 2003; Gray, 2005; Zhu et al., 2007; Ting and Wang, 2008). Moreover, unlike wave-equation migration, which is based on a one-way wave-equation approximation, Gaussian beams are two-way solutions for the outgoing waves and hence have no problems in imaging steep to overturned structures. Gaussian-beam migration also differs from Kirchhoff migration in that it takes into account the local slope of recorded seismic events and is a direction-oriented migration. While Kirchhoff migration treats seismic records as a collection of discrete samples and back spreads each sample onto its entire diffraction traveltime surface, Gaussian-beam migration models data as a superposition of local plane waves and back propagates each plane wave along a narrow beam in the initial direction specified by its local slope. As a result, Gaussian-beam migration spreads a recorded local plane wave only along a small beam front near where it is reflected, reducing migration swinging artifacts due to the imperfect cancellation of the spread energy by destructive interference. The accuracy of the Gaussian-beam method in solving the wave equation is, however, controlled by a beam parameter, or, equivalently, by an initial beam width at the surface (e.g., Červený et al., 1982; Weber, 1988; Hill, 2001). This parameter introduces the needed smoothing in a geometric-ray solution so that the resulting beam solution will be regular everywhere. But in the Gaussian-beam formulation this parameter also determines how the width of a beam will evolve along its central ray. In the cases where the initial width is small, the width of a beam will increase rapidly along its central ray and becomes excessively large at deeper parts of the beam (Červený et al., 1982), which will in turn lead to inaccuracy in traveltime extrapolation along the beam fronts. This can be a problem for land data where near-surface velocity structures have strong lateral variations. The initial beam width in such cases must be sufficiently small in order to accurately image the shallow structures. We have developed a complex-beam method for shot-domain prestack depth migration which overcomes the difficulties encountered by the Gaussian-beam migration. In this article, we describe this method and show its applications on both synthetic and field data.
2 Page 2 of 10 Complex-beam migration Similar to the Gaussian-beam migration method (e.g., Hill, 2001; Gray, 2005; Zhu et al., 2007), the complex-beam migration algorithm consists of two main components: local planewave decomposition of shot records and beam propagation and imaging of the decomposed plane waves. It differs from Gaussian-beam migration mainly in that complex beams, instead of Gaussian beams, are used for wavefield back propagation and imaging. Local plane-wave decomposition Figure 1 is a schematic 2D representation of local plane-wave decomposition of a shot record where S is shot location and R represents receiver locations. The decomposition begins by partitioning the record into overlapping Gaussian-windowed traces, each being centered on a regularly spaced grid point. These grid points are often referred as beam centres. Slant stack is then used to decompose each window of the tapered traces into local plane waves with a range of slopes, corresponding different initial propagation directions. The number of slopes is determined by a formula derived from wavefield sampling theory. Figure 1 shows one of these plane waves propagating from a Gaussian window centered on L. The width of the Gaussian windows is specified by an initial beam width at the surface, and the spacing between the beam centres is determined by the partition of unity of the overlapping windows (Hill, 2001). Our algorithm for the local slant stack is flexible and capable of dealing with both regular and irregular trace spacings. Thus the algorithm provides an additional benefit in performing a uniform spatial gridding or binning of irregular data, reducing the effects of acquisition footprint. Moreover, as antialiasing filtering is carried out as part of decomposition on the slant stacked data based on the slope of each local plane wave, operator anti-aliasing during the migration is no longer needed which is required by Kirchhoff migration. Standard formulation for migration operator anti-aliasing assumes horizontal reflectors and has been found to be detrimental to amplitude preserving of steep reflections. Beam propagation and imaging After data decomposition, each local plane wave from a beam centre is back propagated in our algorithm by a complex beam constructed around a central ray traced from the beam centre with a given initial direction (the blue solid line on right side of Figure 1). Similar to Gaussian beams, complex beams are local solutions to the wave equation and provide complex valued traveltimes and amplitudes for imaging. The two types of beam are, however, constructed in different manners and have different beam properties. A Gaussian beam is obtained by solving a parabolic wave equation along a given central ray (Červený et al., 1982); its beam width expands hyperbolically at a rate inverse proportional to the initial beam width. This property can, as pointed out earlier, lead to excessive expansion for a beam with small initial beam width and result in inaccuracy in imaging the areas with strong near-surface variations. A complex beam, on the other hand, is constructed by solving the wave equation with a complex source location in a ray-centered coordinate system (Zhu and Chun, 1994; Zhu, 2009). The small imaginary component of this complex location ensures that the beam solution is regular everywhere, and the real component is used in a perturbation scheme for calculating the wavefield in the neighbourhood of the central ray. The beam width is determined by the first Fresnel zone and increases along the central ray at a rate proportional to the square root of the initial beam width. Complex beams thus have no difficulties in handling rapidly varying near-surface structures since a smaller initial beam width will result in a narrower beam along the entire central ray and more accurate traveltime extrapolation along the beam fronts. Figure 1. Schematic sketch of the main components of complex-beam migration for a shot record. S indicates the source location and R the receiver locations. L represents a beam center.
3 Page 3 of 10 Similar to the wavefield at the receivers, the point source wavefield at S is also decomposed into local plane waves and propagated by complex beams with different initial propagation directions (Figure 1). A subsurface image is generated from the overlapping area between a pair of source and receiver beams (shaded area in Figure 1). Accumulating the contributions from all shot-receiver beam pairs for a given beam centre produces a beam-center image, and summing over all beam centers from a shot record gives a common-shot image. The final subsurface image is formed by stacking together all individual common-shot images.
4 Page 4 of 10 Figure 2. Depth migrated images from a 2D Canadian Foothills survey: (a) Gaussian-beam migration. (b) Complex beam migration. Seismic waves are often multipathed in a complex velocity structure and the number of the arrivals in general varies from location to location. Retaining such multiple arrivals is difficult in a typical Kirchhoff implementation which computes traveltimes from each shot or receiver to subsurface points and stores these traveltimes in a table for later image accumulation. To accommodate the multiple arrivals, the traveltime tables from all shots and receivers must contain not only multiple entries of varying size at each subsurface point, but also information to ensure that traveltime interpolations are carried out between points on the same
5 Page 5 of 10 traveltime branch. As a result, Kirchhoff migration is typically single valued, selecting at most one traveltime that links shot or receiver to a subsurface point. By replacing the single Kirchhoff traveltime table from each shot or receiver location with multiple beam solutions, the full beam accumulation process described above provides a natural way for retaining all multipathing arrivals in a subsurface image. This, together with the robust behaviour of complex beams, makes our migration algorithm an accurate and flexible tool for depth imaging in geologically complex areas.
6 Page 6 of 10 Figure 3. Depth migrated images from the Foothills model dataset. (a) Velocity model used for generating the synthetic data. (b) Reverse-time migration. (c) Kirchhoff migration. (d) Wave-equation migration. (e) Complex-beam migration. Examples The complex-beam migration has been tested with both field and synthetic data. The first example shows the results from a 2D Canadian Foothills survey which covers a thrust-belt area. The velocity model from this survey is anisotropic (TTI) and has strong lateral variation near the surface. The data was migrated using both Gaussian-beam and complex-beam methods and the results are displayed in Figures 2a and 2b, respectively. Both images are superior to that produced by Kirchhoff migration, as the Kirchhoff method was unable to image the multipathing arrivals generated by the complex velocity model (Zhu et al., 2009). A comparison of the two beam images, however, shows that the shallow part of the survey is better imaged by the complex-beam migration. The dipping events highlighted by the ellipse in Figure 2b, for example, appear more continuous and better defined than those in Figure 2a, showing that complex-beam is more robust for imaging structures with strong near-surface lateral variations. The second example compares complex-beam migration with Kirchhoff, wave-equation, and reverse-time methods using a synthetic dataset from the 2D Canadian Foothills model shown in Figure 3a. The model is about 25km long and 10km deep and is characterized by large lateral variations in both velocity and topography. The velocity ranges from 3600 to 6000 m/s and the maximum elevation difference reaches 1340m. It also contains steep and overturned folds as highlighted by the ellipses in the figure. Based on finite-difference solution to a two-way wave equation, reverse time is the most accurate migration method and its image in Figure 3b is used here as a benchmark for the comparison. For most parts of the model, all four methods give similar results. They all image, for example, the dipping basal reflector at the bottom of the model and the faulted fold about the reflector accurately. The Kirchhoff result (Figure 3c), however, shows more migration swing artifacts because of its wide-spread migration operators and because of its non-recursive implementation. It also failed to image the steep and overturned structures. The waveequation image (Figure 3d) appears to be cleaner than other three, but it was also unable to image the steep structures due to its one-way approximation. Complexbeam migration (Figure 3e), on the other hand, compares well to the reverse time method and imaged both the steep and overturned structures correctly. The final example is from a 3D Alberta Foothills survey. The survey covers a geologically complex area with a large topographic variation. The maximum elevation difference reaches 1240m. The data was migrated using both Kirchhoff and complex-beam algorithms with a TTI velocity model and the migrated sections on one cross-line are shown in Figures 4a and 4b, respectively. Compared to the Kirchhoff section, the beam image is in general cleaner and has less migration swinging artifacts. The solid yellow ellipse in Figure 4a, for example, highlights an area where migration swinging artifacts (indicated by black arrows) are relatively strong, and because of the interference of these artifacts, the events are poorly defined. The corresponding area on Figure 4b, on the other hand, has no visible swinging artifacts and the events are better defined and more interpretable. The dashed ellipses show another shallow area where events are more continuous and better resolved on the beam section. It is interesting to note that the Kirchhoff section below the Nordegg horizon appears to be noisier compared to the shallow part of the section, making the deep reflectors even less traceable than those on the beam section. One possible explanation is that multipathing arrivals are more prevalent below the complex overburden; not properly imaged by the single-valued Kirchhoff migration, some of these arrivals appear as noise in the section.
7 Page 7 of 10 Figure 4. Depth migrated images from a 3D Alberta Foothills survey : (a) Kirchhoff migration. (b) Complex-beam migration. The ellipses highlight areas where beam is superior to Kirchhoff. The black arrows indicate the migration swinging artifacts which are clearly identifiable on Figure 4a but appear to be absent on Figure 4b. The yellow arrows indicate the CMP locations where the commonimages gathers in Figure 5 were generated.
8 Page 8 of 10 Figure 5. Common-offset common-image gathers from three CMP locations indicated in Figure 4. (a) Kirchhoff migration. (b) Complex-beam migration. The ellipses highlight areas where beam is superior to Kirchhoff.
9 Page 9 of 10 One of important applications of the Kirchhoff and beam methods is to generate various surface and subsurface common image gathers (CIGs) for migration velocity analysis. These gathers can be generated much more efficiently by these two methods than by wave-equation and reverse-time as ray and beam solutions contain all the local information needed for computing these gathers (Gray, 2007). Figure 5 shows one type of these gathers common-offset CIGs from three different CMPs near the location indicated by the yellow arrows on the tops of Figures 4a and 4b. Events are picked from these gathers and used as input to a tomographic inversion scheme for velocity determination (e.g., Zhou et al. 2003). A comparison between the complex-beam (Figure 5b) and Kirchhoff (Figure 5a) results shows that reflection events are much better defined on both shallow and deep parts of the complex-beam gathers. The solid ellipses, for example, highlight some shallow events which are well imaged by the beam method but are poorly defined on the Kirchhoff gathers. Similar to the migrated sections in Figure 4, the improvement of the complexbeam over the Kirchhoff method in signal to noise ratio is more noticeable on the deep part of gathers. The dashed ellipses, for example, indicate the deep events which can be clearly identified in Figure 5b, but are buried in noise and barely visible in Figure 5a. Thus more events can be picked up from the complex-beam gathers and picked up more robustly, which will in turn lead to more accurate velocity determination. Conclusions We have developed a complex-beam method for shot-domain prestack depth migration. Based on a complex-beam solution to the wave equation, this method represents an improvement over the Gaussian-beam method. It retains the advantages of Gaussianbeam over Kirchhoff and wave-equation migration, and is more robust in imaging areas with complex structures. Test results with both synthetic and field data show that the complex-beam images are in general superior to those produced by the Kirchhoff method. They are cleaner and have less migration swinging artifacts, and both shallow and deeper structures are better imaged on the beam sections than those on the Kirchhoff sections. The CIGs generated by complex-beam migration have much higher signal to noise ratio than those from Kirchhoff migration, providing more accurate information for velocity determination. The test results also show that the complex-beam method does not suffer the dip limitations of wave-equation migration and compares well to the reverse time in imaging steep to overturned structures. Thus the complexbeam migration provides an accurate and flexible tool for seismic depth imaging in geologically complex land areas. Acknowledgements We thank Sam Gray for discussions about Gaussian-beam and complex-beam migration and for the Foothills synthetic dataset. We thank Talisman Energy Inc. and Calgary Data library of CGGVeritas for the field data examples, Chuck Ursenbach for the synthetic data examples, and Sam Gray and Pam Nagarajappa for reading the paper. References Gray, S. H., 2005, Gaussian beam migration of common-shot records: Geophysics, 70, S71-S77. Gray, S. H., 2007, Angle gathers for Gaussian beam depth migration: 69th Annual International Conference and Exhibition, EAGE, Expanded Abstract, C018. Hill, N. R., 2001, Prestack Gaussian-beam depth migration: Geophysics, 66, Červený, V., M. M. Popov, and I. Pšen ík, 1982, Computation of wave fields in inhomogeneous media Gaussian beam approach, Geophys. J. R. astr. Soc., 70, Nowack, R. L., M. K. Sen, and P. L. Stoffa, 2003, Gaussian beam migration for sparse common-shot and common-receiver data: 73rd Annual International Meeting, SEG, Expanded Abstracts, Ting, C., and D. Wang, 2008, Controlled beam migration applications to Gulf of Mexico: 78th Annual International Meeting, SEG, Expanded Abstracts, Weber, M., 1988, Computation of body-wave seismograms in absorbing 2-D media using the Gaussian beam method: comparison with exact methods: Geophysical Journal International, 92, Zhou, H., Gray, S., Yong, J., and Zhang, Y., 2003, Tomographic residual curvature analysis: The process and its components, 73th Annual international Mtg., Soc. Expl. Geophys., Expanded Abstracts, Zhu, T., 2009, A complex-ray Maslov formulation of beam migration: 79th Annual International Meeting, SEG, Expanded Abstracts, Zhu, T., and K. Chun, 1994, Complex rays in elastic and anelastic media: Geophysical Journal International, 119, Zhu, T., S. H. Gray, and D. Wang, 2007, Prestack Gaussian-beam depth migration in anisotropic media: Geophysics, 72, S133-S138. Zhu, T., T. Krishnasamy, and D. Trad, 2009, Beam migration of Canadian Foothills Datasets, CSEG Convention, Expanded Abstracts, Article Tags: beam, imaging, migration, wave equation Share this Article: Other Reading That May Interest You: Introduction to the Processing Migration / Imaging Focus Articles - 3 tag match(es) (75%) Overview on the Fundamentals of Imaging - 2 tag match(es) (50%)
10 Page 10 of 10 Shear-wave sourced 3-D VSP Depth Imaging of Tight Gas Sandstones in Rulison Field, Colorado - 2 tag match(es) (50%) Non-hyperbolic MultiFocusing Improves Depth Imaging - 1 tag match(es) (25%) Deconvolution After Migration - 1 tag match(es) (25%)
Daniele Colombo* Geosystem-WesternGeco, Calgary, AB M.Virgilio Geosystem-WesternGeco, Milan, Italy.
Seismic Imaging Strategies for Thrust-Belt Exploration: Extended Offsets, Seismic/Gravity/EM simultaneous Joint-Inversion and Anisotropic Gaussian Beam Pre-Stack Depth Migration Daniele Colombo* Geosystem-WesternGeco,
More informationSUMMARY ANGLE DECOMPOSITION INTRODUCTION. A conventional cross-correlation imaging condition for wave-equation migration is (Claerbout, 1985)
Comparison of angle decomposition methods for wave-equation migration Natalya Patrikeeva and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Angle domain common image gathers offer
More informationTu N Fault Shadow Removal over Timor Trough Using Broadband Seismic, FWI and Fault Constrained Tomography
Tu N118 05 Fault Shadow Removal over Timor Trough Using Broadband Seismic, FWI and Fault Constrained Tomography Y. Guo* (CGG), M. Fujimoto (INPEX), S. Wu (CGG) & Y. Sasaki (INPEX) SUMMARY Thrust-complex
More informationDepth Imaging through Surface Carbonates: A 2D example from the Canadian Rocky Mountains
Depth Imaging through Surface Carbonates: A 2D example from the Canadian Rocky Mountains Yong Hua* University of Calgary, Calgary, AB yhua@ucalgary.ca and Don Lawton University of Calgary, Calgary, AB,
More informationNoise suppression and multiple attenuation using full-azimuth angle domain imaging: case studies
first break volume 33, June 2015 special topic Noise suppression and multiple attenuation using full-azimuth angle domain imaging: case studies Aleksander Inozemtsev 1*, Zvi Koren 1 and Alexander Galkin
More informationP S-wave polarity reversal in angle domain common-image gathers
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? P S-wave polarity reversal in angle domain common-image gathers Daniel Rosales and James Rickett 1 ABSTRACT The change in the reflection
More informationH005 Pre-salt Depth Imaging of the Deepwater Santos Basin, Brazil
H005 Pre-salt Depth Imaging of the Deepwater Santos Basin, Brazil Y. Huang* (CGGVeritas), D. Lin (CGGVeritas), B. Bai (CGGVeritas), S. Roby (CGGVeritas) & C. Ricardez (CGGVeritas) SUMMARY Several discoveries,
More informationVelocity Update Using High Resolution Tomography in Santos Basin, Brazil Lingli Hu and Jianhang Zhou, CGGVeritas
Lingli Hu and Jianhang Zhou, CGGVeritas Summary The exploration interest in the Santos Basin offshore Brazil has increased with the large deep water pre-salt discoveries, such as Tupi and Jupiter. As the
More information3D beam prestack depth migration with examples from around the world
A Publication of Petroleum Geo-Services Vol. 8 No. 8 August 2008 3D beam prestack depth migration with examples from around the world Introduction In 1999 AGS specialized in 2D seismic depth processing.
More informationSummary. Introduction
Detailed velocity model building in a carbonate karst zone and improving sub-karst images in the Gulf of Mexico Jun Cai*, Hao Xun, Li Li, Yang He, Zhiming Li, Shuqian Dong, Manhong Guo and Bin Wang, TGS
More informationAnisotropic tomography for TTI and VTI media Yang He* and Jun Cai, TGS
media Yang He* and Jun Cai, TGS Summary A simultaneous anisotropic tomographic inversion algorithm is developed. Check shot constraints and appropriate algorithm preconditioning play an important role
More informationTOM 2.6. SEG/Houston 2005 Annual Meeting 2581
Oz Yilmaz* and Jie Zhang, GeoTomo LLC, Houston, Texas; and Yan Shixin, PetroChina, Beijing, China Summary PetroChina conducted a multichannel large-offset 2-D seismic survey in the Yumen Oil Field, Northwest
More informationPlane-wave migration in tilted coordinates
Stanford Exploration Project, Report 124, April 4, 2006, pages 1 16 Plane-wave migration in tilted coordinates Guojian Shan and Biondo Biondi ABSTRACT Plane-wave migration in tilted coordinates is powerful
More informationSubsalt imaging by common-azimuth migration
Stanford Exploration Project, Report 100, April 20, 1999, pages 113 125 Subsalt imaging by common-azimuth migration Biondo Biondi 1 keywords: migration, common-azimuth, wave-equation ABSTRACT The comparison
More informationImproving subsurface imaging in geological complex area: Structure PReserving INTerpolation in 6D (SPRINT6D)
Improving subsurface imaging in geological complex area: Structure PReserving INTerpolation in 6D (SPRINT6D) Dan Negut* and Mark Ng, Divestco Inc. Copyright 2008, ACGGP This paper was selected for presentation
More informationP137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients
P137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients M. Ayzenberg (StatoilHydro), A. Aizenberg (Institute of Petroleum Geology and Geophysics),
More informationObservations of Azimuthal Anisotropy in Prestack Seismic Data
Observations of Azimuthal Anisotropy in Prestack Seismic Data David Gray* CGGVeritas, Calgary, AB, Canada dave.gray@cggveritas.com Summary A method for displaying prestack seismic data that highlights
More informationCompensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS
Compensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS SUMMARY Anelastic properties of the earth cause frequency dependent energy
More informationWe A10 12 Common Reflection Angle Migration Revealing the Complex Deformation Structure beneath Forearc Basin in the Nankai Trough
We A10 12 Common Reflection Angle Migration Revealing the Complex Deformation Structure beneath Forearc Basin in the Nankai Trough K. Shiraishi* (JAMSTEC), M. Robb (Emerson Paradigm), K. Hosgood (Emerson
More informationAnisotropic Depth Migration and High-Resolution Tomography in Gulf of Mexico: A Case History
Anisotropic Depth Migration and High-Resolution Tomography in Gulf of Mexico: A Case History Gary Rodriguez, Sherry Yang, Diane Yang, Quincy Zhang, Steve Hightower, TGS Summary We present a case study
More informationPART A: Short-answer questions (50%; each worth 2%)
PART A: Short-answer questions (50%; each worth 2%) Your answers should be brief (just a few words) and may be written on these pages if you wish. Remember to hand these pages in with your other exam pages!
More informationJ.A. Haugen* (StatoilHydro ASA), J. Mispel (StatoilHydro ASA) & B. Arntsen (NTNU)
U008 Seismic Imaging Below "Dirty" Salt J.A. Haugen* (StatoilHydro ASA), J. Mispel (StatoilHydro ASA) & B. Arntsen (NTNU) SUMMARY Base and sub salt seismic imaging is still an unresolved issue. To solve
More informationA Petroleum Geologist's Guide to Seismic Reflection
A Petroleum Geologist's Guide to Seismic Reflection William Ashcroft WILEY-BLACKWELL A John Wiley & Sons, Ltd., Publication Contents Preface Acknowledgements xi xiii Part I Basic topics and 2D interpretation
More informationDip-constrained tomography with weighting flow for paleo-canyons: a case study in Para- Maranhao Basin, Brazil Guang Chen and Lingli Hu, CGG
Dip-constrained tomography with weighting flow for paleo-canyons: a case study in Para- Maranhao Basin, Brazil Guang Chen and Lingli Hu, CGG Summary Para-Maranhao Basin offshore Brazil is well-known for
More informationDeconvolution imaging condition for reverse-time migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 83 96 Deconvolution imaging condition for reverse-time migration Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT The reverse-time
More informationStanford Exploration Project, Report 115, May 22, 2004, pages
Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 248 Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 First-order lateral interval velocity estimates without picking
More informationFinite difference elastic modeling of the topography and the weathering layer
Finite difference elastic modeling of the topography and the weathering layer Saul E. Guevara and Gary F. Margrave ABSTRACT Finite difference 2D elastic modeling is used to study characteristics of the
More information3D VTI traveltime tomography for near-surface imaging Lina Zhang*, Jie Zhang, Wei Zhang, University of Science and Technology of China (USTC)
Downloaded 01/03/14 to 16.01.198.34. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ 3D VTI traveltime tomography for near-surface imaging Lina Zhang*, Jie
More informationChapter 7. Seismic imaging. 7.1 Assumptions and vocabulary
Chapter 7 Seismic imaging Much of the imaging procedure was already described in the previous chapters. An image, or a gradient update, is formed from the imaging condition by means of the incident and
More informationDepth Imaging for Unconventional Reservoir Characterization: Canadian Plains Case Study
Depth Imaging for Unconventional Reservoir Characterization: Canadian Plains Case Study Bill Goodway 1, Greg Purdue 1, Shiang Yong Looi 2, Lijuan (Kathy) Du 2, Mark Rowland 2 1 Apache Canada, 2 Schlumberger
More informationSUMMARY INTRODUCTION THEORY
Seismic modeling using the frozen Gaussian approximation Xu Yang, University of California Santa Barbara, Jianfeng Lu, Duke University, and Sergey Fomel, University of Texas at Austin SUMMARY We adopt
More informationVelocity model building in complex media by multi-layer non-linear slope tomography. Summary Introduction Non-linear slope tomography
Velocity model building in complex media by multi-layer non-linear slope tomography. Patrice Guillaume, Jean-Philippe Montel*, Xiaoming Zhang, Gilles Lambaré, Anthony Prescott, Matthieu Reinier, Richard
More informationFrom PZ summation to wavefield separation, mirror imaging and up-down deconvolution: the evolution of ocean-bottom seismic data processing
From PZ summation to wavefield separation, mirror imaging and up-down deconvolution: the evolution of ocean-bottom seismic data processing Sergio Grion, CGGVeritas Summary This paper discusses present
More informationThe Deconvolution of Multicomponent Trace Vectors
The Deconvolution of Multicomponent Trace Vectors Xinxiang Li, Peter Cary and Rodney Couzens Sensor Geophysical Ltd., Calgary, Canada xinxiang_li@sensorgeo.com Summary Deconvolution of the horizontal components
More informationMatrix formulation of adjoint Kirchhoff datuming
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 377 Matrix formulation of adjoint Kirchhoff datuming Dimitri Bevc 1 ABSTRACT I present the matrix formulation of Kirchhoff wave equation datuming
More informationDownloaded 11/08/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Velocity tomography based on turning waves and finite-frequency sensitivity kernels Xiao-Bi Xie*, University of California at Santa Cruz, Jan Pajchel, Statoil, Norway Summary We present a velocity tomography
More informationThe i-stats: An Image-Based Effective-Medium Modeling of Near-Surface Anomalies Oz Yilmaz*, GeoTomo LLC, Houston, TX
The i-stats: An Image-Based Effective-Medium Modeling of Near-Surface Anomalies Oz Yilmaz*, GeoTomo LLC, Houston, TX Summary Near-surface modeling for statics corrections is an integral part of a land
More informationConverted-wave beam migration with sparse sources or receivers
Geophysical Prospecting, 2015 doi: 10.1111/1365-2478.12226 Converted-wave beam migration with sparse sources or receivers Lorenzo Casasanta and Samuel H. Gray CGG, 715 5th Avenue SW, Calgary,, AB T2P 5A2
More informationElements of 3D Seismology Second Edition
Elements of 3D Seismology Second Edition Copyright c 1993-2003 All rights reserved Christopher L. Liner Department of Geosciences University of Tulsa August 14, 2003 For David and Samantha And to the memory
More informationAcoustic anisotropic wavefields through perturbation theory
GEOPHYSICS, VOL. 78, NO. 5 (SEPTEMBER-OCTOBER 2013); P. WC41 WC50, 22 FIGS. 10.1190/GEO2012-0391.1 Acoustic anisotropic wavefields through perturbation theory Tariq Alkhalifah 1 ABSTRACT Solving the anisotropic
More informationImplicit 3-D depth migration by wavefield extrapolation with helical boundary conditions
Stanford Exploration Project, Report 97, July 8, 1998, pages 1 13 Implicit 3-D depth migration by wavefield extrapolation with helical boundary conditions James Rickett, Jon Claerbout, and Sergey Fomel
More informationTomography for Static Corrections and Prestack Depth Imaging
Tomography for Static Corrections and Prestack Depth Imaging Xianhuai Zhu, Ph.D. Fusion Petroleum Technologies Inc. 25231 Grogan's Mill Road, Suite 175 The Woodlands, TX 77380, USA Summary Turning-ray
More informationTOM 1.7. Sparse Norm Reflection Tomography for Handling Velocity Ambiguities
SEG/Houston 2005 Annual Meeting 2554 Yonadav Sudman, Paradigm and Dan Kosloff, Tel-Aviv University and Paradigm Summary Reflection seismology with the normal range of offsets encountered in seismic surveys
More informationDownloaded 05/01/17 to Redistribution subject to SEG license or copyright; see Terms of Use at
Mapping Imbricate Structures in the Thrust Belt of Southeast Turkey by Large-Offset Seismic Survey Oz Yilmaz*, Anatolian Geophysical, Istanbul, Turkey; and Serdar Uygun, Ali Ölmez, and Emel Çalı, Turkish
More informationAn empirical method for estimation of anisotropic parameters in clastic rocks
An empirical method for estimation of anisotropic parameters in clastic rocks YONGYI LI, Paradigm Geophysical, Calgary, Alberta, Canada Clastic sediments, particularly shale, exhibit transverse isotropic
More informationWave-equation tomography for anisotropic parameters
Wave-equation tomography for anisotropic parameters Yunyue (Elita) Li and Biondo Biondi ABSTRACT Anisotropic models are recognized as more realistic representations of the subsurface where complex geological
More informationFast Kirchhoff migration in the wavelet domain
Exploration Geophysics (2002) 33, 23-27 Fast Kirchhoff migration in the wavelet domain Valery A. Zheludev 1 Eugene Ragoza 2 Dan D. Kosloff 3 Key Words: Kirchhoff migration, wavelet, anti-aliasing ABSTRACT
More informationComparison of two physical modeling studies of 3D P-wave fracture detection
Comparison of two physical modeling studies of 3D P-wave fracture detection ZHONGPING QIAN 1,2, XIANG-YANG LI 1 AND SHANGXU WANG 3 1 British Geological Survey, Murchison House, West Mains Road, Edinburgh
More informationImprovement of stacking image by anisotropic velocity analysis using P-wave seismic data
P-46 Improvement of stacking image by anisotropic velocity analysis using P-wave seismic data Laxmidhar Behera*, and Prakash Khare, National Geophysical Research Institute (NGRI), Hyderabad Summary Anisotropy
More informationSEG Houston 2009 International Exposition and Annual Meeting
TTI/VTI anisotropy parameters estimation by focusing analysis, Part I: theory Jun ai*, Yang He, Zhiming Li, in Wang, Manhong Guo TGS-Nopec Geophysical ompany, itywest lvd. Suite, Houston, TX 7742, US Summary
More informationAngle-domain common-image gathers from anisotropic Gaussian beam migration and its application to anisotropy-induced imaging errors analysis
J. Earth Syst. Sci. (217) 126: 1 c Indian Academy of Sciences DOI 1.17/s124-16-783-x Angle-domain common-image gathers from anisotropic Gaussian beam migration and its application to anisotropy-induced
More informationComplex structural imaging of transition zones in Bohai Bay, China, by OBC technology
Complex structural imaging of transition zones in Bohai Bay, China, by OBC technology Xiaogui Miao 1, Wei Yan 1, Yongxia Liu 1, Joe Zhou 2, and Kunlun Yang 2 Abstract A series of ocean-bottom cable (OBC)
More informationLowrank RTM for converted wave imaging
Lowrank RM for converted wave imaging Lorenzo Casasanta * (CGG), Zhiguang Xue (U Austin), Sam Gray (CGG) Summary his paper is an attempt to fill the technology gap existing between pure P- and PS-wave
More informationStatics preserving projection filtering Yann Traonmilin*and Necati Gulunay, CGGVeritas
Yann Traonmilin*and Necati Gulunay, CGGVeritas Summary Projection filtering has been used for many years in seismic processing as a tool to extract a meaningful signal out of noisy data. We show that its
More informationStatic Corrections for Seismic Reflection Surveys
Static Corrections for Seismic Reflection Surveys MIKE COX Volume Editors: Series Editor: Eugene F. Scherrer Roland Chen Eugene F. Scherrer Society of Exploration Geophysicists Tulsa, Oklahoma Contents
More informationF-K Characteristics of the Seismic Response to a Set of Parallel Discrete Fractures
F-K Characteristics of the Seismic Response to a Set of Parallel Discrete Fractures Yang Zhang 1, Xander Campman 1, Samantha Grandi 1, Shihong Chi 1, M. Nafi Toksöz 1, Mark E. Willis 1, Daniel R. Burns
More informationResidual Statics using CSP gathers
Residual Statics using CSP gathers Xinxiang Li and John C. Bancroft ABSTRACT All the conventional methods for residual statics analysis require normal moveout (NMO) correction applied on the seismic data.
More informationUNIVERSITY OF CALGARY. Prestack depth migration methods for isotropic and polar anisotropic media. Xiang Du A THESIS
UNIVERSITY OF CALGARY Prestack depth migration methods for isotropic and polar anisotropic media by Xiang Du A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILMENT OF THE REQUIREMENTS
More informationTh Guided Waves - Inversion and Attenuation
Th-01-08 Guided Waves - Inversion and Attenuation D. Boiero* (WesternGeco), C. Strobbia (WesternGeco), L. Velasco (WesternGeco) & P. Vermeer (WesternGeco) SUMMARY Guided waves contain significant information
More informationX040 Buried Sources and Receivers in a Karsted Desert Environment
X040 Buried Sources and Receivers in a Karsted Desert Environment C. Berron* (CGGVeritas Services S.A.), E. Forgues (CGGVeritas Services S. A.), M. Jervis (Saudi Aramco), A. Bakulin (Saudi Aramco) & R.
More informationSEG/New Orleans 2006 Annual Meeting. Non-orthogonal Riemannian wavefield extrapolation Jeff Shragge, Stanford University
Non-orthogonal Riemannian wavefield extrapolation Jeff Shragge, Stanford University SUMMARY Wavefield extrapolation is implemented in non-orthogonal Riemannian spaces. The key component is the development
More informationReflection Seismic Method
Reflection Seismic Method Data and Image sort orders; Seismic Impedance; -D field acquisition geometries; CMP binning and fold; Resolution, Stacking charts; Normal Moveout and correction for it; Stacking;
More informationApplication of Seismic Reflection Surveys to Detect Massive Sulphide Deposits in Sediments-Hosted Environment
IOSR Journal of Applied Geology and Geophysics (IOSR-JAGG) e-issn: 2321 0990, p-issn: 2321 0982.Volume 3, Issue 4 Ver. I (Jul - Aug. 2015), PP 46-51 www.iosrjournals.org Application of Seismic Reflection
More information( ) ( ), 1, 1. Downloaded 10/14/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Chao Zhang*, Jie Zhang, University of Science and Technology of China (USTC); Zhikun Sun, GeoTomo LLC Summary Refraction traveltimes have long been applied for deriving long-wavelength statics solutions.
More informationSynthetic seismic modelling and imaging of an impact structure
Modelling and imaging of an impact structure Synthetic seismic modelling and imaging of an impact structure Matteo Niccoli ABSTRACT A geologic depth and velocity model of an impact crater was created and
More informationVollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich
Seminar des FA Ultraschallprüfung Vortrag 1 More info about this article: http://www.ndt.net/?id=20944 Vollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich Thomas
More informationInvestigating fault shadows in a normally faulted geology
Investigating fault shadows in a normally faulted geology Sitamai Ajiduah* and Gary Margrave CREWES, University of Calgary, sajiduah@ucalgary.ca Summary Fault shadow poses a potential development risk
More informationSummary. Introduction
Multi-layer tomography and its application for improved depth imaging Patrice Guillaume, Steve Hollingworth*, Xiaoming Zhang, Anthony Prescott, Richard Jupp, Gilles Lambaré, Owen Pape, CGGVeritas Summary
More informationNear Surface Anomalies and Solutions
Near Surface Anomalies and Solutions Todd Mojesky* CGG Canada Services, Ltd., Calgary, Alberta, Canada tmojesky@cgg.com and David LeMeur and Ali Karagul, CGG Canada Services, Ltd., Calgary, Alberta, Canada
More informationFUNDAMENTALS OF SEISMIC EXPLORATION FOR HYDROCARBON
FUNDAMENTALS OF SEISMIC EXPLORATION FOR HYDROCARBON Instructor : Kumar Ramachandran 10 14 July 2017 Jakarta The course is aimed at teaching the physical concepts involved in the application of seismic
More informationOBS wavefield separation and its applications
P-088 OBS wavefield separation and its applications Sergio Grion*, CGGVeritas Summary This paper discusses present trends in ocean-bottom seismic (OBS) data processing. Current industrial practices for
More informationTh P7 02 A Method to Suppress Salt-related Converted Wave Using 3D Acoustic Modelling
Th P7 0 A Method to Suppress Salt-related Converted Wave Using 3D Acoustic Modelling J. Kumar* (Petroleum Geo-Services), M. Salem (ENI E&P), D.E. Cegani (ENI E&P) Summary Converted waves can be recorded
More informationH Seismic imaging in and around salt bodies: problems and pitfalls. Ian Jones 1, Ian Davison 2. Summary
H1-2-3 Seismic imaging in and around salt bodies: problems and pitfalls Ian Jones 1, Ian Davison 2 1 ION, Chertsey, UK, 2 EarthMoves, Egham, UK Summary Salt movement often results in steeply-dipping complex
More informationPluto 1.5 2D ELASTIC MODEL FOR WAVEFIELD INVESTIGATIONS OF SUBSALT OBJECTIVES, DEEP WATER GULF OF MEXICO*
Pluto 1.5 2D ELASTIC MODEL FOR WAVEFIELD INVESTIGATIONS OF SUBSALT OBJECTIVES, DEEP WATER GULF OF MEXICO* *This paper has been submitted to the EAGE for presentation at the June 2001 EAGE meeting. SUMMARY
More informationA synthetic example of anisotropic P -wave. processing for a model from the Gulf of Mexico
A synthetic example of anisotropic P -wave processing for a model from the Gulf of Mexico Baoniu Han, Tagir Galikeev, Vladimir Grechka, Jérôme Le Rousseau and Ilya Tsvankin Center for Wave Phenomena, Department
More informationIMAGING WITH REVERSE VERTICAL SEISMIC PROFILES USING A DOWNHOLE, HYDRAULIC, AXIAL VIBRATOR
IMAGING WITH REVERSE VERTICAL SEISMIC PROFILES USING A DOWNHOLE, HYDRAULIC, AXIAL VIBRATOR Mary Krasovec and Roger Turpening Earth Resources Laboratory Department of Earth, Atmospheric, and Planetary Sciences
More informationSeismic tests at Southern Ute Nation coal fire site
Seismic tests at Southern Ute Nation coal fire site Sjoerd de Ridder and Seth S. Haines ABSTRACT We conducted a near surface seismic test at the Southern Ute Nation coal fire site near Durango, CO. The
More informationDirty salt velocity inversion: The road to a clearer subsalt image
GEOPHYSICS. VOL. 76, NO. 5 (SEPTEMBER-OCTOBER 2011); P. WB169 WB174, 8 FIGS. 10.1190/GEO2010-0392.1 Dirty salt velocity inversion: The road to a clearer subsalt image Shuo Ji 1, Tony Huang 1, Kang Fu 2,
More informationImaging complex structure with crosswell seismic in Jianghan oil field
INTERPRETER S CORNER Coordinated by Rebecca B. Latimer Imaging complex structure with crosswell seismic in Jianghan oil field QICHENG DONG and BRUCE MARION, Z-Seis, Houston, Texas, U.S. JEFF MEYER, Fusion
More informationP Wave Reflection and Refraction and SH Wave Refraction Data Processing in the Mooring, TN Area
P Wave Reflection and Refraction and SH Wave Refraction Data Processing in the Mooring, TN Area Abstract: Author: Duayne Rieger Home Institution: Slippery Rock University of Pennsylvania REU Institution:
More informationPost-stack inversion of the Hussar low frequency seismic data
Inversion of the Hussar low frequency seismic data Post-stack inversion of the Hussar low frequency seismic data Patricia E. Gavotti, Don C. Lawton, Gary F. Margrave and J. Helen Isaac ABSTRACT The Hussar
More informationChapter 1. Introduction EARTH MODEL BUILDING
Chapter 1 Introduction Seismic anisotropy in complex earth subsurface has become increasingly important in seismic imaging due to the increasing offset and azimuth in modern seismic data. To account for
More informationANGLE-DEPENDENT TOMOSTATICS. Abstract
ANGLE-DEPENDENT TOMOSTATICS Lindsay M. Mayer, Kansas Geological Survey, University of Kansas, Lawrence, KS Richard D. Miller, Kansas Geological Survey, University of Kansas, Lawrence, KS Julian Ivanov,
More informationRough sea estimation for phase-shift de-ghosting Sergio Grion*, Rob Telling and Seb Holland, Dolphin Geophysical
Rough sea estimation for phase-shift de-ghosting Sergio Grion*, Rob Telling and Seb Holland, Dolphin Geophysical Summary This paper discusses rough-sea de-ghosting for variabledepth streamer data. The
More informationINTEGRATED GEOPHYSICAL INTERPRETATION METHODS FOR HYDROCARBON EXPLORATION
INTEGRATED GEOPHYSICAL INTERPRETATION METHODS FOR HYDROCARBON EXPLORATION Instructor : Kumar Ramachandran 31 July 4 August 2017 Jakarta COURSE OUTLINE The course is aimed at imparting working knowledge
More informationEdinburgh Research Explorer
Edinburgh Research Explorer Nonlinear Scattering-based Imaging in Elastic Media Citation for published version: Ravasi, M & Curtis, A 213, 'Nonlinear Scattering-based Imaging in Elastic Media' Paper presented
More informationThe Estimation of PS Wave Statics By Using SCDCC + Stack Power Method
The Estimation of PS Wave Statics By Using SCDCC + Stack Power Method ABSTRACT Summary Aimin Xue * Petrol Sound Geoservices. 8, MacEwan Ridge Close, N.W., Calgary, Alberta T3K3A7, Canada. aiminx@telus.net
More informationSeismic Reflection Views of the 1994 Northridge Earthquake Hypocentral Region Using Aftershock Data and Imaging Techniques
Seismic Reflection Views of the 1994 Northridge Earthquake Hypocentral Region Using Aftershock Data and Imaging Techniques SERGIO CHÁVEZ-PÉREZ and JOHN N. LOUIE SEISMOLOGICAL LABORATORY/174, MACKAY SCHOOL
More informationAnalysis of multicomponent walkaway vertical seismic profile data
Analysis of multicomponent walkaway vertical seismic profile data Bona Wu, Don C. Lawton, and Kevin W. Hall ABSTRACT A multicomponent walkaway VSP data processed for PP and PS imaging as well to study
More informationImproving Resolution with Spectral Balancing- A Case study
P-299 Improving Resolution with Spectral Balancing- A Case study M Fatima, Lavendra Kumar, RK Bhattacharjee, PH Rao, DP Sinha Western Offshore Basin, ONGC, Panvel, Mumbai Summary: The resolution limit
More informationSeismic Waves in Complex 3 D Structures, 26 (2016), (ISSN , online at
Kirchhoff prestack depth migration in simple orthorhombic and triclinic models with differently rotated elasticity tensor: comparison with zero-offset travel-time perturbations Václav Bucha Department
More informationVáclav Bucha. Department of Geophysics Faculty of Mathematics and Physics Charles University in Prague. SW3D meeting June 6-7, 2016 C OM S TR 3 D
Kirchhoff prestack depth migration in simple orthorhombic and triclinic models with differently rotated elasticity tensor: comparison with zero-offset travel-time perturbations Václav Bucha Department
More informationTowards full waveform inversion: A torturous path
FWI Towards full waveform inversion: A torturous path J. Helen Isaac and Gary F. Margrave ABSTRACT Full waveform inversion (FWI) can be viewed as an iterative cycle involving forward modelling, pre-stack
More informationFloatSeis Technologies for Ultra-Deep Imaging Seismic Surveys
FloatSeis Technologies for Ultra-Deep Imaging Seismic Surveys 25 th January, 2018 Aleksandr Nikitin a.nikitin@gwl-geo.com Geology Without Limits Overview 2011-2016 GWL Acquired over 43000 km 2D seismic
More informationComputation of refraction static corrections using first-break traveltime differences
GEOPHYSICS, VOL. 54, NO. IO (OCTOBER 1989); P. 1289-1296, 7 FIGS., 2 TABLES. Computation of refraction static corrections using first-break traveltime differences Don C. Lawton* ABSTRACT Differences in
More informationTitle: Exact wavefield extrapolation for elastic reverse-time migration
Title: Exact wavefield extrapolation for elastic reverse-time migration ummary: A fundamental step of any wave equation migration algorithm is represented by the numerical projection of the recorded data
More informationAn Improved Dual Sensor Summation Method with Application to Four-Component (4-C) Seafloor Seismic Data from the Niger Delta
Earth Science Research; Vol. 4, No. 2; 2015 ISSN 1927-0542 E-ISSN 1927-0550 Published by Canadian Center of Science and Education An Improved Dual Sensor Summation Method with Application to Four-Component
More informationElastic full waveform inversion for near surface imaging in CMP domain Zhiyang Liu*, Jie Zhang, University of Science and Technology of China (USTC)
Elastic full waveform inversion for near surface imaging in CMP domain Zhiyang Liu*, Jie Zhang, University of Science and Technology of China (USTC) Summary We develop an elastic full waveform inversion
More informationTime vs depth in a North Sea case study of the ugly truth Abstract
Time vs depth in a North Sea case study of the ugly truth Thomas Massip 1, Lauren Braidwood 1, Juergen Fruehn 1, Owen Isaac 1, Jonathan Denly 1, Robert Richardson 2, Phil Mollicone 3 1 ION Geophysical;
More informationAVO inversion in V (x, z) media
Stanford Exploration Project, Report 97, July 8, 998, pages 75 94 AVO inversion in V (x, z) media Yalei Sun and Wenjie Dong keywords:.5-d Kirchhoff integral, AVO inversion, fluid-line section ABSTRACT
More information