Progress Report on: Interferometric Interpolation of 3D SSP Data
|
|
- Sabina Mosley
- 6 years ago
- Views:
Transcription
1 Progress Report on: Interferometric Interpolation of 3D SSP Data Sherif M. Hanafy ABSTRACT We present the theory and numerical results for interferometrically interpolating and extrapolating 3D marine SSP data. For the interpolation of SSP data we use the combination of a natural Green s function (SSP shot gathers) and a model-based Green s function for the water-layer model. Synthetic results show that the aliased SSP data with large receiver intervals can be accurately interpolated to smaller intervals. The resulting virtual shot gather contains some artifacts, a 2D multi-channel local image matching filter is used after interferometric interpolation to remove these artifacts. SSP data extrapolation also employs a natural Green s function and a model-based Green s function. This suggests that a sparse marine SSP survey can yield much more information about the reflectors if data are interpolated and extrapolated by interferometry. This assumes that the sources are located both far outside and within the recording aperture. Limitations of this theory are: interpolated traces are limited to 4 virtual traces between each two real ones and with extrapolation, near offsets virtual traces show some artifacts. INTRODUCTION Typical marine surface seismic profiles (SSP) are ideally designed for a regular recording grid, but in reality it suffers from irregularities in the recording geometry. Coarse receiver spacing and narrow recording apertures are common problem in marine acquisition, especially in the crossline direction. The result can be inadequate subsurface illumination and distortions in the migration image. To alleviate this problem, various algorithms were suggested to fill in the missing traces in marine data (Abma and Kabir, 2006; van Dedem and Verschuur, 2005; Muijs et al., 2007). Most of them require certain assumptions 27 such as, the data have sparse representation in a certain domain, or rely on a priori information, such as knowledge of seismic velocities (Baumstein et al., 2005). In all cases, there is no use of the redundant information available in the free-surface multiples. In this paper we show how to interpolate marine SSP data by interferometry. This method transforms surface and seabed related multiples into primaries recorded at virtual receiver both inside and outside the receiver array. In this proposed technique, no velocity model for the sediments is needed. The interpolated data can be used for migration, velocity analysis, and tomography. The multiples also can be used to illuminate much wider areas of the subsurface. This work is the 3D extension of the 2D interferometric method of Dong and Hanafy (2008) This paper is divided into three parts. First, the theory of interferometric interpolation and extrapolation is presented. Then, tests are carried out on 3D synthetic data generated from a multi-layer model. Finally, a brief conclusion and discussion is given. THEORY Wang et al. (2009) applied the interferometric interpolation method to surface seismic profile (SSP) data. In that method, the surface related multiples are used to predict missing traces in gaps. Sparse SSP data can also be interpolated where an SSP SSP correlation transform is used. This transform is similar in spirit to the VSP VSP correlation transform (Schuster, 2009), except now we use a two-layer model-based Green s function with a band limited point source, as shown in Figure 1. The diagrams show how SSP traces (with both sources and receivers just below the free surface) are correlated with sparsely distributed SSP traces to generate a dense distribution of SSP traces. This correlation operation is required by the acoustic reciprocity equation of correlation type for a two-state system, where one state is the acoustic field associated with the multi-layered model shown in
2 28 Hanafy Figure 1 and the other state is associated with a sea-floor model. The reciprocity equation of correlation type can be described as (Appendix A): G(B A) G 0 (A B) = [G 0 (x B) G(x A) (1) S s n x G(x A) G 0(x B) n x ]d 2 x, where, G 0 (x B) is the Green s function for a water layer model and G(x A) is the Green s function for the actual earth model (Dong and Hanafy, 2008). Here, S s is the boundary just below the sea surface and the integration along the free surface vanishes because both Green s functions are zero there. The contributions from the vertical boundaries at infinity to the left and to the right of the boat will be ignored (Wapenaar and Fokkema, 2006). Here, A and B are just below the free surface and below the horizontal source line S o The above equation is a reciprocity equation of correlation type for two different states, which can be used for interpolation or extrapolation of traces. The far-field approximation to equation 1 yields the SSP SSP far-field transform SSP {}}{ G(B A) = 2ik S s SSP SSP SSP {}}{{}}{{}}{ G(x A) G 0 (x B) dx 2 + G 0 (A B), (2) To implement this equation, the SSP data are used to estimate G(x A) and a finite-difference solution to the wave equation is used to estimate G 0 (x B) for the ocean-layer model that only contains the free surface and ocean bottom interfaces. This FD calculation is possible because the sea floor topography is well known beneath any exploration survey. The key idea for interpolation is that the free surface acts as a perfectly reflecting mirror so that 2 nd and 3 rd views, i.e., free-surface related multiples, of the subsurface can be used to fill in the trace gaps, as indicated by Figure 1. NUMERICAL TEST Multi-layer Velocity Model To test the proposed approach, a 3D acoustic multi-layer velocity model (Figure 2) is generated and finite-difference (FD) solutions of the acoustic wave equation are used to simulate SSP traces in that model. The velocity model is 3000 x 3000 x 1400 m 3 with a 300 by 300 receiver points at receiver interval of 10 m in both X and Y directions. The water layer has a thickness of 750 m and the water bottom is a flat layer. At this stage, one common shot gather is generated, where the source is located at location (10,10,30) m for x, y, and z directions. In the following section we show the 2D test results and then the 3D results. Details for implementing the 2D interpolation and extrapolation method are given in Dong and Hanafy (2008) 2D Interpolation Results A 2D line is extracted from the CSG s to test the proposed approach. The extracted 2D line is located at X = 10 m and has 300 traces in the Y direction with a trace interval of 10 m. Figure 3a shows the original CSG and Figure 3b shows the sparse version of the extracted CSG, where only 75 traces are kept with a trace interval of 40 m. Equation 2 is used to interpolate the missing traces and Figure 3c shows the virtual CSG after one iteration before application of the matching filter. To eliminate the artifacts generated during interpolation process we use the following 2-steps; Use a 2D multi-channel local image matching filter. This matching filter can be repeated several times to minimize the artifacts. Repeat the interpolation - matching filter process several times, where the input CSG to each new iteration is the virtual CSG from the previous iteration. In this work we used 8 interpolation iterations, each interpolation iteration is followed by 10 iterations of matching filter. The final virtual CSG is shown in Figure 3d. Comparing Figures 3a and 3d shows a very good correlation between the virtual CSG and the original CSG. Artifacts at the far offsets are not completely eliminated due to the paucity of true traces at the far offsets that could be used for interpolation. 3D Interpolation Results The proposed approach is tested on a 3D CSG, where Figure 4a shows a part of the original CSG. The sparse CSG contains 60 parallel lines with 100 receivers per line (Figure 4b). The inline offset is 30 m and the crossline offset is 50 m. The sparse version of the CSG contains 6000 traces, and our target is to interpolate missing traces and get a virtual CSG with 10 m offsets in both inline and crossline directions. Equation 2 is used to interpolate the missing traces, where we used 3 interpolation iterations, each interpolation iteration is followed by 8 iterations of matching filter. Figure 4c shows the final virtual CSG with matching filter. Comparing the original (Figure 4a) and virtual (Figure 4c) CSGs shows a very good correlation. SEG/EAGE Velocity Model Another velocity model, the SEG/EAGE velocity model, is used to test the proposed approach. In this test we simulated a near azimuth marine survey with 8 streamers each one has 170 receiver. The crrossline interval is 30 m and the inline interval is 12 m. Near offset is at 0 m, while far offset is at 2028 m. Our target is to interpolate virtual traces in both the inline and crossline directions,
3 3D SSP Interpolation 29 so that the final results has 22 streamers and 805 receiver per streamer. The cross line interval is 10 m and the inline interval is 4 m. Figure 5 shows the SEG/EAGE velocity model we used in this test, Figure 6a shows a sample of the input CSG (1 st and 2 nd streamers), Figure 6b shows the virtual CSG (1 st to 4 th streamers), and Figure 6c shows the real CSG for comparison. Here, the artifacts are highly eliminated and the virtual CSG (Figure 6b) shows a very good correlation with the real CSG (Figure 6c). CONCLUSIONS An interferometric method for interpolating and extrapolating marine SSP data is tested on synthetic multi-layer and SEG/EAGE velocity models. The results show that this method can kinematically interpolate the sparse SSP data to a dense receiver distribution data and can widen the receiver aperture significantly. The least squares image matching filter is shown to suppress artifacts and can partly correct for amplitudes. No velocity model is required to apply this approach, however, the thickness of the water layer and the velocity of seismic wave in the water is required. The interpolated data can be used for migration, velocity analysis, and tomography. The multiples also can be used to illuminate much wider areas of the subsurface. The virtual CSGs are very good correlated with the real CSGs. With interpolation, the number of receivers can be increased 4 times. We do not recommend increase the number of receivers more than 4 times, since we got more artifacts with increasing number of virtual traces. The extrapolation is in principle similar to the interpolation, our extrapolated traces (not included in this work) show higher artifacts at the near offset zone. A future challenge is improve the extrapolated traces and eliminate the artifacts. ACKNOWLEDGMENTS We would like to thank the members of the 2008 University of Utah Tomography and Modeling/Migration (UTAM) Consortium ( APPENDIX A ACOUSTIC RECIPROCITY EQUATION OF CORRELATION TYPE WITH TWO STATES Consider two states, one is the acoustic field associated with the multi-layered model shown in Figure 1 where G(x A) is interpreted as the acoustic wave field excited by an interior harmonic point source at A and recorded at x. And the other state is the acoustic field in the seafloor model, which only consists of a water layer, a free surface, and a sea floor below which is a homogeneous medium with the same velocity and density of layer 1. The Green s function associated with this new state is defined as G 0 (x B), and does not contain reflections from any interface below the sea floor. The Helmholtz equations satisfied by these two Green s functions are ( 2 + k 2 )G(x A) = δ(x A), (A-1) ( 2 + k 2 0)G 0 (x B) = δ(x B), (A-2) where k = ω/v(x) for the multi-layered model and k 0 = ω/v 0 (x) for the sea-floor model. The Laplacians can be weighted by Green s functions to give the following identities: G(x A) 2 G 0 (x B) = [G(x A) G 0 (x B)] G(x A) G 0 (x B), G 0 (x B) 2 G(x A) = [G 0 (x B) G(x A)] G 0 (x B) G(x A). (A-3) Instead of defining the integration volume over the entire multi-layered model the volume is restricted to the ocean layer where both the sea-floor model and multi-layer model agree. Subtracting and integrating the above equations over the ocean volume yields the reciprocity equation of correlation type for two different states: G(B A) G 0 (A B) = [G 0 (x B) G(x A) S s n x G(x A) G 0(x B) ]d 2 x, (A-4) n x where S s is the boundary along the sea floor and the integration along the free vanishes because both Green s functions are zero there. The contributions from the infinite vertical boundaries to the left and right of the boat will be ignored. Equation A-4 is a reciprocity equation of correlation type for two different states, which can be used for the interpolation or extrapolation of traces. REFERENCES Abma, R. and N. Kabir, 2006, 3d interpolation of irregular data with a pocs algorithm: Geophysics, 71, E91 E97. Baumstein, A., M. T. Hadidi, L. Hinkley, and W. S. Ross, 2005, A practical procedure for application of 3d srme to conventional marine data: TLE, Dong, S. and S. M. Hanafy, 2008, Interferometric interpolation and extrapolation of sparse obs data: UTAM Annual Meeting, Muijs, R., J. O. A. Robertsson, and K. Holliger, 2007, Prestack depth migration of primary and surfacerelated multiple reflections: Part 1 - imaging: Geophysics, 72, S59 S69. Schuster, G. T., 2009, Seismic interferometry: Cambridge Press, under press. van Dedem and D. Verschuur, 2005, 3d surface related multiple prediction: A sparse inversion approach: Geophysics, 70, V31 V43. Wang, Y., Y. Luo, and G. Schuster, 2009, Interferometric interpolation of missing seismic data: Geophysics, in press.
4 30 Hanafy Figure 1: Ray diagrams for transforming SSP data to SSP data. Here, the open geophones indicate the locations of virtual geophones where traces are created from the original SSP data recorded at the filled geophone positions. Both G(x A) and G o (x B) can be data based Green s functions, but in this case G o (x B) is model-based and so is computed for the two-layer sea-floor model. Figure 2: The 3D multi-layer velocity model used to test interpolation technique
5 3D SSP Interpolation 31 Figure 3: The 2D test results. a) The original common shot gather with a receiver interval of 10 m. b) The sparse SSP CSG used for the interpolation test, where the receiver interval is 40 m. c) The interpolation results before using the matching filter, where the interpolated data have a receiver interval of 10 m. d) The interpolation results after using the matching filter. Most artifacts are removed and the wavelet is consistent with the true data.
6 32 Hanafy Figure 4: The 3D test results. a) One of the original common shot gathers, here the receiver interval is 10 m. b) The sparse SSP CSG used for the interpolation test, where 60 parallel lines with 100 receivers per line are used. The inline receiver interval is 30 m, the crossline spacing is 50 m, and total number of receivers in the input data are c) The interpolation results after using the matching filter. Most artifacts are removed and the wavelet is consistent with the true data. The virtual CSGs contains 300 lines with 300 receiver per line. The inline receiver interval is 10 m, the crossline spacing is 10 m, and total number of receivers in the output data are 90,000.
7 3D SSP Interpolation 33 Figure 5: The SEG/EAGE velocity model used to test interpolation technique Wapenaar, K. and J. Fokkema, 2006, Green s function representations for seismic interferometry: Geophysics, 71, S133 S146.
8 34 Hanafy Figure 6: The 1 st and 2 nd streamers of the original common shot gather, here the receiver interval is 12 m, streamers offset is 30 m, and each streamer has 170 receivers.
9 3D SSP Interpolation 35 Figure 7: The interpolation results of the SEG/EAGE velocity model. a) The interpolation results after using the matching filter. Most of the artifacts are removed and the wavelet is consistent with the true data. The virtual CSGs contains 22 crosslines with 508 receiver per crossline. The inline receiver interval is 4 m, the crossline spacing is 10 m, and total number of receivers in the output data are 11,176. b) Synthetically generated CSG with the same geometry as in (a).
Interferometric interpolation of sparse marine data
Geophysical Prospecting, 2014, 62, 1 16 doi: 10.1111/1365-2478.12066 Interferometric interpolation of sparse marine data Sherif M. Hanafy 1,2 and Gerard T. Schuster 1 1 King Abdullah University of Science
More informationAPPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA.
APPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA. BY: ISAAC KUMA YEBOAH. Department of Engineering, Regent University
More informationDownward propagation: using true near-surface velocity model, one-way wave- equation extrapolator:
Class 9: Beyond Seismic Statics Corrections Wed, Oct 7, 2009 Dynamic corrections: wave-equation and wavefield datuming methods Skipping the near-surface by interferometry (no model needed) Common-Focus
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 informationPractical implementation of SRME for land multiple attenuation
Practical implementation of SRME for land multiple attenuation Juefu Wang* and Shaowu Wang, CGGVeritas, Calgary, Canada juefu.wang@cggveritas.com Summary We present a practical implementation of Surface
More informationPredicting rugged water-bottom multiples through wavefield extrapolation with rejection and injection
Predicting rugged water-bottom multiples through wavefield extrapolation with rejection and injection Benxi Ke ABSTRACT Although convolution-based and WEM-modeling-based methods for predicting surface-related
More informationCLASSIFICATION OF MULTIPLES
Introduction Subsurface images provided by the seismic reflection method are the single most important tool used in oil and gas exploration. Almost exclusively, our conceptual model of the seismic reflection
More informationEnhanced Angular Illumination from Separated Wavefield Imaging (SWIM)
Enhanced Angular Illumination from Separated Wavefield Imaging (SWIM) S. Lu* (Petroleum Geo-Services), N.D. Whitmore (Petroleum Geo- Services), A.A. Valenciano (Petroleum Geo-Services) & N. Chemingui (Petroleum
More informationG012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry
G012 Scattered Ground-roll Attenuation for 2D Land Data Using Seismic Interferometry D.F. Halliday* (Schlumberger Cambridge Research), P.J. Bilsby (WesternGeco), J. Quigley (WesternGeco) & E. Kragh (Schlumberger
More informationSummary. Introduction
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 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 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 informationMulti-source Least-squares Migration of Gulf of Mexico Data
Multi-source Least-squares Migration of Gulf of Mexico Data Xin Wang. King Abdullah University of Science and Technology, Thuwal 955-69, Kingdom of Saudi Arabia Corresponding author is Xin Wang. E-mail
More informationDownloaded 09/20/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
Joint SRME and model-based water-layer demultiple for ocean bottom node Hui Huang*, Ping Wang, Jing Yang, Hui Chen (CGG); Pierre-Olivier Ariston, Imtiaz Ahmed, and Nick Bassett (BP) Summary Ocean bottom
More information3D SRME: An Acquisition Independent Approach
Summary 3D SRME: An Acquisition Independent Approach A.Pica, P. Lanfranchi*, Ph. Herrmann CGG, Paris, France Conventional data-driven SRME techniques do not require any a priori knowledge of the subsurface
More informationAutofocusing imaging: Imaging with primaries, internal multiples and free-surface multiples
CWP-797 Autofocusing imaging: Imaging with primaries, internal multiples and free-surface multiples Satyan Singh 1, Roel Snieder 1, Jyoti Behura 1,3, Joost van Der Neut 2, Kees Wapenaar 2, and Evert Slob
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 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 informationY015 Complementary Data-driven Methods for Interbed Demultiple of Land Data
Y015 Complementary Data-driven Methods for Interbed Demultiple of Land Data S. Sonika* (WesternGeco), A. Zarkhidze (WesternGeco), J. Heim (WesternGeco) & B. Dragoset (WesternGeco) SUMMARY Interbed multiples
More informationAnisotropy-preserving 5D interpolation by hybrid Fourier transform
Anisotropy-preserving 5D interpolation by hybrid Fourier transform Juefu Wang and Shaowu Wang, CGG Summary We present an anisotropy-preserving interpolation method based on a hybrid 5D Fourier transform,
More informationDownloaded 01/20/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
D Multi-source Least-squares Reverse Time Migration Wei Dai, C. Boonyasiriwat and Gerard T. Schuster, King Abdullah University of Science and Technology Downloaded 1//1 to 19.171.17.1. Redistribution subject
More 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 informationP. Bilsby (WesternGeco), D.F. Halliday* (Schlumberger Cambridge Research) & L.R. West (WesternGeco)
I040 Case Study - Residual Scattered Noise Attenuation for 3D Land Seismic Data P. Bilsby (WesternGeco), D.F. Halliday* (Schlumberger Cambridge Research) & L.R. West (WesternGeco) SUMMARY We show that
More informationReverse time migration of multiples: Applications and challenges
Reverse time migration of multiples: Applications and challenges Zhiping Yang 1, Jeshurun Hembd 1, Hui Chen 1, and Jing Yang 1 Abstract Marine seismic acquisitions record both primary and multiple wavefields.
More informationAttenuation of water-layer-related multiples Clement Kostov*, Richard Bisley, Ian Moore, Gary Wool, Mohamed Hegazy, Glenn Miers, Schlumberger
Clement Kostov*, Richard Bisley, Ian Moore, Gary Wool, Mohamed Hegazy, Glenn Miers, Schlumberger Summary We present a method for modeling and separation of waterlayer-related multiples in towed streamer
More informationH003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography
H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography M. Cavalca* (Schlumberger - Westerngeco) & R.P. Fletcher (Schlumberger - Westerngeco) SUMMARY Estimation of the
More informationA least-squares shot-profile application of time-lapse inverse scattering theory
A least-squares shot-profile application of time-lapse inverse scattering theory Mostafa Naghizadeh and Kris Innanen ABSTRACT The time-lapse imaging problem is addressed using least-squares shot-profile
More informationAttenuation of diffracted multiples with an apex-shifted tangent-squared radon transform in image space
Attenuation of diffracted multiples with an apex-shifted tangent-squared radon transform in image space Gabriel Alvarez, Biondo Biondi, and Antoine Guitton 1 ABSTRACT We propose to attenuate diffracted
More informationCFP migration; practical aspects
13 CFP migration; practical aspects Jan Thorbecke 1 13.1 Introduction In the past year the Common Focus Point (CFP) technology has become an important paradigm in the DELPHI research program. The interpretation
More 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 informationWe ELI1 02 Evaluating Ocean-bottom Seismic Acquisition in the North Sea - A Phased Survey Design Case Study
We ELI1 02 Evaluating Ocean-bottom Seismic Acquisition in the North Sea - A Phased Survey Design Case Study M. Branston* (Schlumberger Geosolutions), R. Campbell (Schlumberger Geosolutions), M. Rowlands
More informationThe ups and downs of ocean-bottom seismic processing: Applications of wavefield separation and up-down deconvolution
SPECIAL M u l t SECTION: i c o m p o Mn ue ln t t i cs oe mi s p m o i nc e n t s e i s m i c The ups and downs of ocean-bottom seismic processing: Applications of wavefield separation and up-down deconvolution
More informationImaging with multiples using LSRTM
Chapter 4 Imaging with multiples using LSRTM In this chapter, I present a technique for imaging both primaries and higher-order multiples using joint least-squares reverse-time migration (joint-lsrtm).
More informationSuccesses and challenges in 3D interpolation and deghosting of single-component marinestreamer
Successes and challenges in 3D interpolation and deghosting of single-component marinestreamer data James Rickett*, Schlumberger Gould Research Summary Combining deghosting with crossline interpolation
More informationU043 3D Prestack Time Domain Full Waveform Inversion
U043 3D Prestack Time Domain Full Waveform Inversion D.V. Vigh* (WesternGeco), W.E.S. Starr (WesternGeco) & K.D. Kenneth Dingwall (WesternGeco) SUMMARY Despite the relatively high computational demand,
More informationShallow Reverberation Prediction Methodology with SRME
Shallow Reverberation Prediction Methodology with SRME S.R. Barnes* (Petroleum Geo-Services), R.F. Hegge (Petroleum Geo- Services), R. van Borselen (Petroleum Geo-Services) & J. Owen (Petroleum Geo-Services)
More informationMultiple attenuation for shallow-water surveys: Notes on old challenges and new opportunities
Multiple attenuation for shallow-water surveys: Notes on old challenges and new opportunities Clement Kostov 1, Frederico Xavier de Melo 1, Abhishek Raj 1, Alexander Zarkhidze 1, Alex Cooke 1, Glenn Miers
More informationSubsalt Marchenko imaging: A Gulf of Mexico example
CWP-888 Subsalt Marchenko imaging: A Gulf of Mexico example Xueyi Jia, Antoine Guitton, Satyan Singh, and Roel Snieder ABSTRACT Marchenko redatuming allows one to use surface seismic reflection data to
More informationInversion after depth imaging
Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.
More 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 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 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 informationLeast Squares Kirchhoff Depth Migration: potentials, challenges and its relation to interpolation
Least Squares Kirchhoff Depth Migration: potentials, challenges and its relation to interpolation Daniel Trad CGG Summary Least Squares Migration (LSM), like interpolation, has the potential to address
More informationStanford Exploration Project, Report 124, April 4, 2006, pages 49 66
Stanford Exploration Project, Report 124, April 4, 2006, pages 49 66 48 Stanford Exploration Project, Report 124, April 4, 2006, pages 49 66 Mapping of specularly-reflected multiples to image space: An
More informationE044 Ray-based Tomography for Q Estimation and Q Compensation in Complex Media
E044 Ray-based Tomography for Q Estimation and Q Compensation in Complex Media M. Cavalca* (WesternGeco), I. Moore (WesternGeco), L. Zhang (WesternGeco), S.L. Ng (WesternGeco), R.P. Fletcher (WesternGeco)
More informationPressure and shear wave separation of ocean bottom seismic data
Pressure and shear wave separation of ocean bottom seismic data Ohad Barak ABSTRACT PZ summation is a common method for separating the upgoing wavefield from the downgoing wavefield in data acquired by
More informationImaging Salt Dome Flank and Dipping Sediments using Time Reversed Acoustics
Imaging Salt Dome Flank and Dipping Sediments using Time Reversed Acoustics Rongrong Lu, Mark Willis, Xander Campman, Jonathan Ajo-Franklin, Nafi Toksöz Earth Resources Laboratory Dept. of Earth, Atmospheric
More informationShot-based pre-processing solutions for wide azimuth towed streamer datasets
first break volume 25, March 2007 focus on azimuth Shot-based pre-processing solutions for wide azimuth towed streamer datasets Philippe Herrmann, 1* Gordon Poole, 2 Antonio Pica, 1 Sylvain Le Roy, 1 and
More informationSubsalt steep dip imaging study with 3D acoustic modeling Lei Zhuo* and Chu-Ong Ting, CGGVeritas
Lei Zhuo* and Chu-Ong Ting, CGGVeritas Summary We present a 3D acoustic wave equation modeling study with the objective of understanding imaging challenges for steep dips (faults and three-way closure)
More informationImproved Imaging through Pre-stack Trace Interpolation for missing offsets of OBC data A case study from North Tapti area of West Coast, India
P-256 Improved Imaging through Pre-stack Trace Interpolation for missing offsets of OBC data A case study from North Tapti area of West Coast, India M Lal, CPS Rana, Ramji Pathak, BN Bhatta, DP Sinha,
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 informationHigh Resolution Imaging by Wave Equation Reflectivity Inversion
High Resolution Imaging by Wave Equation Reflectivity Inversion A. Valenciano* (Petroleum Geo-Services), S. Lu (Petroleum Geo-Services), N. Chemingui (Petroleum Geo-Services) & J. Yang (University of Houston)
More 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 informationSUMMARY LEAST-SQUARES MIGRATION THEORY
Making the most out of the least (squares migration) Gaurav Dutta, Yunsong Huang, Wei Dai, Xin Wang, and G.T. Schuster King Abdullah University of Science and Technology SUMMARY Standard migration images
More informationSummary. Introduction
Multivessel coil shooting acquisition with simultaneous sources Nick Moldoveanu¹, Ying Ji², Craig Beasley¹ ¹WesternGeco, ²Schlumberger Cambridge Research Summary Multivessel coil shooting is a towed-streamer
More informationRefraction Full-waveform Inversion in a Shallow Water Environment
Refraction Full-waveform Inversion in a Shallow Water Environment Z. Zou* (PGS), J. Ramos-Martínez (PGS), S. Kelly (PGS), G. Ronholt (PGS), L.T. Langlo (PGS), A. Valenciano Mavilio (PGS), N. Chemingui
More 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 informationDirections in 3-D imaging - Strike, dip, both?
Stanford Exploration Project, Report 113, July 8, 2003, pages 363 369 Short Note Directions in 3-D imaging - Strike, dip, both? Marie L. Clapp 1 INTRODUCTION In an ideal world, a 3-D seismic survey would
More informationThree critical concerns for marine seismics with portable systems. Source strength and tuning Streamer length 2D vs. 3D
Three critical concerns for marine seismics with portable systems Source strength and tuning Streamer length 2D vs. 3D Airgun Source Single airgun Multiple airguns signal strength volume 1/3 1 x (200 in
More informationWe Survey Design and Modeling Framework for Towed Multimeasurement Seismic Streamer Data
We-04-12 Survey Design and Modeling Framework for Towed Multimeasurement Seismic Streamer Data K. Eggenberger* (Schlumberger), P. Christie (Schlumberger), M. Vassallo (Schlumberger) & D.J. van Manen (Schlumberger)
More informationDownloaded 09/03/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Full-waveform inversion in a shallow water environment: A North Sea 3D towed-streamer data example Kathy Zou*, Lars Tore Langlo, Grunde Rønholt, Jaime Ramos-Martinez and Steve Kelly, PGS Summary We apply
More informationAn illustration of adaptive Marchenko imaging
An illustration of adaptive Marchenko imaging Joost van der Neut 1, Kees Wapenaar 1, Jan Thorbecke 1, Evert Slob 1, and Ivan Vasconcelos 2 Abstract In Marchenko imaging, wavefields are retrieved at specified
More informationChapter 1. 2-D field tests INTRODUCTION AND SUMMARY
Chapter 1 2-D field tests INTRODUCTION AND SUMMARY The tomography method described in the preceding chapter is suited for a particular class of problem. Generating raypaths and picking reflectors requires
More informationQuest for subsalt steep dips Chu-Ong Ting*, Lei Zhuo, and Zhuoquan Yuan, CGGVeritas
Chu-Ong Ting*, Lei Zhuo, and Zhuoquan Yuan, CGGVeritas Summary Advances in wide azimuth (WAZ) acquisition, TTI model building, and Reverse Time Migration (RTM) have brought significant improvements in
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 informationA Novel 3-D De-multiple Workflow for Shallow Water Environments - a Case Study from the Brage field, North Sea
A Novel 3-D De-multiple Workflow for Shallow Water Environments - a Case Study from the Brage field, North Sea J. Oukili* (PGS), T. Jokisch (PGS), A. Pankov (PGS), B. Farmani (PGS), G. Ronhølt (PGS), Ø.
More information1D internal multiple prediction in a multidimensional world: errors and recommendations
1D internal multiple prediction 1D internal multiple prediction in a multidimensional world: errors and recommendations Pan Pan and Kris Innanen ABSTRACT Internal multiples are more difficult to estimate
More 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 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 informationFlexi Binning : A Solution for Irregular Acquisition Geometry
5th Conference & Exposition on Petroleum Geophysics, Hyderabad-2004, India PP 199-203 Flexi Binning : A Solution for Irregular Acquisition Geometry M. S. Rawat 1, M. Goswamy 1 & M. Das 2 1 GEOPIC, ONGC,
More information3D predictive deconvolution for wide-azimuth gathers
3D predictive deconvolution for wide-azimuth gathers P.Hugonnet(1), J.L.Boelle(2), P.Herrmann(1), F.Prat(1), S.Navion(1) (1) CGGVeritas, (2) Total E&P Summary Usual pre-stack predictive deconvolution solutions
More informationTh E An Interferometric Interpretation of Marchenko Redatuming
Th E0 An Interferometric Interpretation of Marchenko Redatuming J. van der Neut (Delft University of Technology), I. Vasconcelos (chlumberger Gould Research) & K. Wapenaar (Delft University of Technology)
More informationAngle-dependent reflectivity by wavefront synthesis imaging
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 477 Angle-dependent reflectivity by wavefront synthesis imaging Jun Ji 1 ABSTRACT Elsewhere in this report, Ji and Palacharla (1994) show
More informationWave-equation inversion prestack Hessian
Stanford Exploration Project, Report 125, January 16, 2007, pages 201 209 Wave-equation inversion prestack Hessian Alejandro A. Valenciano and Biondo Biondi ABSTRACT The angle-domain Hessian can be computed
More informationChallenges and Opportunities in 3D Imaging of Sea Surface Related Multiples Shaoping Lu*, N.D. Whitmore and A.A. Valenciano, PGS
Challenges and Opportunities in 3D Imaging of Sea Surface Related Multiples Shaoping Lu*, N.D. Whitmore and A.A. Valenciano, PGS Summary Conventional shot domain migration constructs a subsurface image
More informationEffect of Structure on Wide Azimuth Acquisition and Processing
P - 71 Effect of Structure on Wide Azimuth Acquisition and Processing Bruce VerWest and Dechun Lin, CGGVeritas, Houston, TX Summary This model study shows that narrow azimuth acquisition may be adequate
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 informationChallenges of pre-salt imaging in Brazil s Santos Basin: A case study on a variable-depth streamer data set Summary
Challenges of pre-salt imaging in Brazil s Santos Basin: A case study on a variable-depth streamer data set Jeremy Langlois, Bing Bai, and Yan Huang (CGGVeritas) Summary Recent offshore discoveries in
More informationSeismic data Interpolation in the Continuous Wavenumber Domain, Flexibility and Accuracy
Seismic data Interpolation in the Continuous Wavenumber Domain, Flexibility and Accuracy Ye Zheng Geo-X Exploration Services Inc. Summary Our recently developed algorithm, ASFT (Arbitrarily Sampled Fourier
More informationTh ELI1 12 Joint Crossline Reconstruction and 3D Deghosting of Shallow Seismic Events from Multimeasurement Streamer Data
Th ELI1 12 Joint Crossline Reconstruction and 3D Deghosting of Shallow Seismic Events from Multimeasurement Streamer Data Y.I. Kamil* (Schlumberger), M. Vassallo (Schlumberger), W. Brouwer (Schlumberger),
More informationMitigation of the 3D Cross-line Acquisition Footprint Using Separated Wavefield Imaging of Dual-sensor Streamer Seismic
Mitigation of the 3D Cross-line Acquisition Footprint Using Separated Wavefield Imaging of Dual-sensor Streamer Seismic A.S. Long* (PGS), S. Lu (PGS), D. Whitmore (PGS), H. LeGleut (PGS), R. Jones (Lundin
More informationEstimation of multiple scattering by iterative inversion, Part II: Practical aspects and examples
GEOPHYSICS, VOL. 62, NO. 5 (SEPTEMBER-OCTOBER 1997); P. 1596 1611, 18 FIGS. Estimation of multiple scattering by iterative inversion, Part II: Practical aspects and examples D. J. Verschuur* and A. J.
More informationSeismic Attributes on Frequency-enhanced Seismic Data
Seismic Attributes on Frequency-enhanced Seismic Data Satinder Chopra* Arcis Corporation, Calgary, Canada schopra@arcis.com Kurt J. Marfurt The University of Oklahoma, Norman, US and Somanath Misra Arcis
More 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 informationConsiderations in 3D depth-specific P-S survey design
Considerations in 3D depth-specific P-S survey design Don C. Lawton and Peter W. Cary 3D P-S survey design ABSTRACT A new sparse-shot design for 3D P-S surveys is introduced. In the sparse shot design
More informationPerforming depth imaging is an essential part of deepwater
SPECIAL Seismic SECTION: modeling S e i s m i c m o d e l i n g The Tempest Project Addressing challenges in deepwater Gulf of Mexico depth imaging through geologic models and numerical simulation ADAM
More informationAVO Analysis with Multi-Offset VSP Data
AVO Analysis with Multi-Offset VSP Data Z. Li*, W. Qian, B. Milkereit and E. Adam University of Toronto, Dept. of Physics, Toronto, ON, M5S 2J8 zli@physics.utoronto.ca T. Bohlen Kiel University, Geosciences,
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 informationu = v is the Laplacian and represents the sum of the second order derivatives of the wavefield spatially.
Seismic Elastic Modelling Satinder Chopra, Arcis Corporation, Calgary, Canada Summary Seismic modeling experiments were performed to understand the seismic responses associated with a variety of subsurface
More informationWe C 07 Least-squares Inversion for Imaging the Full Seismic Wavefield
We C 07 Least-squares Inversion for Imaging the Full Seismic Wavefield S. Lu* (Petroleum Geo-Services), F. Liu (Petroleum Geo-Services), N. Chemingui (Petroleum Geo-Services), M. Orlovich (Petroleum Geo-Services)
More informationSUMMARY THEORY INTRODUCTION
Acoustic 3D least-squares reverse time migration using the energy norm Daniel Rocha, Paul Sava & Antoine Guitton Center for Wave Phenomena, Colorado School of Mines SUMMARY We propose a least-squares reverse
More 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 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 informationMarchenko imaging: Imaging with primaries, internal multiples, and free-surface multiples
GEOPHYSICS, VOL. 8, NO. 5 (SEPTEMBER-OCTOBER 215); P. S165 S174, 16 FIGS., 3 TABLES. 1.119/GEO214-494.1 Downloaded 9/21/15 to 138.67.128.57. Redistribution subject to SEG license or copyright; see Terms
More informationG017 Beyond WAZ - A Modeling-based Evaluation of Extensions to Current Wide Azimuth Streamer Acquisition Geometries
G017 Beyond WAZ - A Modeling-based Evaluation of Extensions to Current Wide Azimuth Streamer Acquisition Geometries M. Cvetkovic* (ION Geophysical), Z. Zhou (ION Geophysical / GXT Imaging Solutions) &
More informationA Review of Current Marine Demultiple Techniques with Examples from the East Coast of Canada
A Review of Current Marine De Techniques with Examples from the East Coast of Canada ABSTRACT R.Brooymans*, T.Mojesky and L.Pham CGG Canada Services Ltd., Calgary rbrooymans@ca.cgg.com Introduction Multiple
More informationTu N Internal Multiple Attenuation on Radial Gathers With Inverse-scattering Series Prediction
Tu N114 05 Internal Multiple Attenuation on Radial Gathers With Inverse-scattering Series Prediction F. Xavier de Melo* (Schlumberger), C. Kostov (Schlumberger) & J. Wu (Schlumberger) SUMMARY We present
More informationInternal Multiple Attenuation on Radial Gathers With Inverse- Scattering Series Prediction
Internal Multiple Attenuation on Radial Gathers With Inverse- Scattering Series Prediction Frederico Xavier de Melo, Clément Kostov, James Wu Schlumberger Summary We present a novel workflow for pre-stack
More informationExamples of GLOBE Claritas Processing
V6.0 Examples of GLOBE Claritas Processing Refraction Statics Removal of Noise (Land, 3D) Removal of Swell Noise Interpolation : shots/receivers Interpolation : 5D (STITCH) Demultiple : High Resolution
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 information