Simultaneous joint inversion of refracted and surface waves Simone Re *, Claudio Strobbia, Michele De Stefano and Massimo Virgilio - WesternGeco
|
|
- Anabel Simon
- 5 years ago
- Views:
Transcription
1 Simultaneous joint inversion of refracted and surface waves Simone Re *, Claudio Strobbia, Michele De Stefano and Massimo Virgilio - WesternGeco Summary In this paper, we review the near-surface challenges and present a new technique to address issues of shallow velocity model building. The method is based on simultaneous joint inversion of surface waves and firstbreak picks. We show its application to a 2D line acquired in Egypt. Introduction The shallow subsurface often exhibits large and rapid vertical and horizontal variations due to structural and stratigraphic reasons. The low compaction and cementation of shallow formations can generate extreme velocity changes, especially in arid areas without a shallow water table. From a geological and geophysical point of view, the near surface can be extremely complex, with multiple velocity inversion and sharp lateral velocity variations. It is well known that the impact of near-surface perturbations on seismic reflection data must be removed to obtain the correct geometrical image of deeper horizons and representative reservoir attributes. Additionally, the near surface is an important, but challenging, portion of velocity model building in depth imaging. The near-surface characterization is, therefore, an important part of seismic data processing, in particular for land data. Conventional approaches involve the use of refracted waves, or diving waves, for P-wave velocity model estimation. The use of Rayleigh waves can be an alternative when data are acquired properly. Refraction tomography (RT) and surface-wave inversion (SWI) are based on different physical principles, make use of different components of the wavefield, and have different limitations and strengths. Refraction tomography often provides deeper models and estimates directly the compressional wave velocity. But, with land data, refraction techniques can be challenging in areas with a complex near surface. The data quality can be critical; for instance, picking the near offset can be difficult. Often, only a limited offset range can be picked reliably. Moreover, velocity inversions and hidden layers can produce ambiguities, and they might not be resolved uniquely. Surface waves have a very high resolution in the shallow near surface and are very robust versus model complexity and data quality. However, even if the low-frequency penetration reaches hundreds of meters, the vertical resolution decreases with depth and the resolution at the final investigation target may be insufficient in some areas. Finally, conversion from shear to compressional velocities requires calibration with P-wave events. When RT and SWI are used, integration and reconciliation of the two models can be done in different ways. A robust framework for this task is simultaneous joint inversion (SJI). The two measurements are input into a single inverse problem, where a unique model with multiple properties is estimated, minimizing the data prediction error and a link between the two domains. The synergies between the two techniques can be exploited. Different portions of the model are resolved by the interacting contributions of the measurements. Surface waves and refractions in land data Land seismic acquisition usually involves deployment of sources and receivers at the surface. Such a configuration implies the generation and propagation of different wave types, including several events propagating in the near surface, primarily refracted waves and surface waves. Rayleigh waves, the so-called ground roll, constitute the largest component of the source radiation. Refracted P- waves are far less energetic, but often can be easily identified in gathers as first arrivals. Historically, surface waves were usually regarded only as coherent noise, and in conventional land acquisition, receiver arrays attempted to attenuate Rayleigh waves in the field. On the other hand, point-receiver acquisition with low-frequency sources and receivers provides a correct sampling of the Rayleigh wave modes and enables implementation of methods based on the analysis and inversion of surface waves for near-surface characterization. Raw point-receiver data are processed to extract the propagation properties, usually the dispersion curve, which is finally inverted to get a velocity profile. A comprehensive workflow for surface-wave analysis based on continuous and adaptive surface-wave processing, has been developed for land 3D data (Strobbia et al., 2009). This surface-wave processing workflow extracts the local Rayleigh wave properties, in particular the modal dispersion curves. For each given location within a survey, the Rayleigh wave properties are estimated considering a frequency-dependent receiver aperture to honor the frequency-dependent near-field effects, intrinsic attenuation, and variable lateral resolution. The result for an individual receiver line can be plotted as 2D dispersion pseudosection (Figure 1). Each column of the pseudosection represents the local dispersion curve for the considered mode. Multiple receiver and source lines are processed to get a dispersion volume with 3D geometries. The surface-wave propagation properties are then inverted to a shear-wave velocity (V s ) model, matching modeled properties to the measured ones. The vertical and lateral 1914
2 resolution is high and the method can solve velocity inversions and complex near-surface structures, but the investigation depth is generally limited to 100 to 200 m. Figure 1: Modal dispersion pseudosection for a receiver line. Refraction-based techniques, on the other hand, make use of the part of the body wave energy that is (continuously) refracted in the near surface and is observed in records as first arrivals. Refraction tomography tries to determine the subsurface velocity by means of fastest raypaths associated with first-break arrivals. First arrival traveltimes are picked and then inverted to a compressional wave velocity (V P ) model. Given a source-receiver pair and an interval velocity model, it is possible to compute first-arrival time by integrating the slowness along the trajectory of the refracted (i.e., fastest) raypath: t s( x, y, z) dr (1), = r where t is the total traveltime, r is the ray trajectory, and s is the slowness. V 1 V 2 V 3 V 4 V 5 V 6 V 7 V 8 V 9 V 10 V 11 V 12 Figure 2: Example of a refracted ray for a given source-receiver pair and the formulation of RT problem for a simple 2D case. The subsurface is divided into cells: V P is constant inside the cell and the size of the cells defines the spatial resolution of tomography. Making the subsurface discrete leads to the discretization of equation 1 that can be expressed as: N t i = s j j= 1 l (2), where t i is the traveltime for the i-th source-receiver pair, s j is the slowness of the j-th cell of the model, and l is the length of the i-th raypath inside the j-th cell. Equation 2 is written for the 2D case shown in Figure 2, but the extension to the 3D case is straightforward. Collecting equations like equation 2 for each sourcereceiver pair allows building a linear system as: ij ij t = Ls (3), with t being the vector containing all traveltimes, L the matrix of all raypaths in each model cell, and s the slowness vector. Therefore, given a guess velocity (slowness) model, linear system 3 provides traveltimes for all source-receivers pairs (forward modeling). The difference between traveltimes computed using expression 3 and traveltimes coming from first-break picking can be used to invert the linear system and then to update the initial guess velocity model (inversion process). The tomography must then iterate to try to converge on the best estimate of the true model by minimizing the differences between the observed data (refracted pick times) and those computed by ray tracing for the current guess of the model. In areas characterized by shallow dipping reflectors and/or high lateral velocity contrasts, the process of Dix inversion, being based on a flat-layered earth model, cannot give an accurate description of the subsurface. In the same context, refraction tomography is a powerful tool for improving both time and depth processing. On the time-processing side, the interval velocity coming from tomography can be used for accurate static correction computation, while on the depth-processing side, the same velocity model is usually the first step in the model-building workflow by means of the depth migration process and CMP residuals analysis. However, due to issues related to inversion stability, the resolution of a velocity model achieved by refraction tomography is lower than that needed for depth imaging techniques. In addition, some geological configurations, such as cases with velocity inversions, cannot be resolved without independent information. These limitations can be overcome by solving refraction and surface-wave problems simultaneously. RT and SWI present several synergies. They have different depths of investigation, different intrinsic limitations, and different distribution of sensitivity and resolution with depth, but the two measured physical properties (V P and V S ) have a strong correlation. The joint use of the two techniques has, therefore, several advantages. The high lateral resolution and robustness of surface waves allow resolving shallow complex structures such as velocity inversions. Surface-wave data can fill the information gap due to missing offset ranges in the picked traveltimes, for instance, the noisy near offset. On the other hand, the joint inversion of the two measurements allows robust V S to V P conversion, and refraction data provide deeper velocity models. The final model can be used to compute near-surface corrections (for instance, statics) for depth imaging and other data processing applications. Additionally, it provides useful information for near-surface geological modeling in which the V P /V S ratio can play an important role. 1915
3 The simultaneous joint inversion method The term, joint inversion, is commonly used in the oil and gas industry to indicate a wide range of technologies and workflows that aim to integrate different measurements for geophysical exploration. Dell Aversana (2001) integrated seismic and electromagnetic data for structural imaging, Li and Oldenburg (1996b) used borehole and surface magnetic data to invert for susceptibility, and De Stefano and Colombo (2007) inverted linked data within a single cost function. This approach is called simultaneous joint inversion, given that the workflow integrates the measurements in the inversion phase and it is not simply an alternating sequence of single measurement inversions. SJI uses the same regularization and preconditioning that the single domains use for standard inversions. This is one of the key benefits provided by SJI; each domain is regularized (and/or preconditioned) as in single-domain inversions. Most importantly, SJI solves simultaneously for different scales, so that we can handle separately the denser V S unknowns to define precisely the very shallow portions of the models and the coarser V P unknowns to invert down to the main refractors of the model (Figure 4). Within this approach, SJI maintains the single-domain operators and workflows, reaching the benefits of a common inversion as described in the following section. Figure 4: Left: Grid used for refraction tomography. Right: Grid used for surface-wave inversion: short wavelengths (blue) have less penetration than long wavelengths (green). Different measurements invert data on different grids. Example Figure 3: Left: Single-domain approach to invert for V S and V P. Right: The SJI workflow to invert simultaneously for the two different measurements. For the first time, we propose SJI as a technology to invert simultaneously any combination of refracted, reflection, and surface-wave seismic data. In particular, in this paper we focus on real data applications in West Gharib (Egypt) using refracted and surface waves only, and we use a single objective function that is defined and minimized (De Stefano and Colombo, 2007), in contrast to an approach in which multiple objective functions are inverted in separate domains. As shown in Figure 3, SJI acts at the inversion phase where a kernel of the objective function is built from three different elements: 1. Residuals from different domains (misfits in the dispersion curves for V S and delays for first arrivals for V P ). 2. Single-domain constraints (second-order regularizations for lateral and vertical consistency). 3. Inter-domain constraints (empirical law linking V S and V P and some geometrical constraint forcing the models to have the same anomalies). To show the benefits of the presented method, we discuss the results of a 2D line, about 10-km long, extracted from a point-receiver 3D survey. The site is a gravel plain between the Gulf of Suez and the Red Sea Mountains in Egypt. Figure 5: Shear-wave velocity section from surface-wave inversion. Point-receiver data were acquired with 7.5-m receiver spacing and a vibroseis source. For the whole line, data show the presence of strong Rayleigh waves dominated by a fundamental mode, with large lateral velocity variations. Both the surface-wave and refracted-wave properties are estimated on point-receiver data. The extracted dispersion pseudosection and the first breaks are first inverted with a single-domain approach, independently. The result of surface-wave single-domain inversion is plotted in Figure 5. The velocity model does not indicate 1916
4 the presence of major velocity inversion in the near surface. There are, however, sharp lateral velocity variations. It is important to notice that the low-velocity zones identified by the V S inversion are extremely relevant from a geological point of view. They correspond to fault zones, they are observed throughout the whole survey area, and they are consistent with the structural setting of the region (Laake et al., 2010). They are spatially consistent and evident in the raw dispersion volume before inversion. In Figure 6, a wavelength slice is presented, with the location of the 2D line. The final results show good agreement with the shallow stack section, plotted in Figure 8. Figure 8: Shallow stack section of the considered 2D line. Conclusions Figure 6: Principle of surface-wave characterization and dispersion section. The simultaneous joint inversion of the two measurements is run using a link defined on the basis of the two singledomain inversions. The result of the simultaneous joint inversion is shown in Figure 7. Figure 7: Final V P section from SJI of refracted first arrivals and Rayleigh wave dispersion. This final V P section is obtained using the information from both measurements. The deeper part is controlled only by the refracted wave data, while in the shallow portion, the SW data provides high resolution. The final traveltime misfit of the SJI inversion is not larger than the one of the single domain refraction tomography. It is important to remark that, in this case, surface-wave data have confirmed the absence of critical features, such as velocity inversions. The presented method allows construction of a near-surface model, merging in a rational manner, the contribution of refracted waves and surface waves, usually treated as noise and simply removed from seismic records. Surface and refracted waves have different resolution, penetration, and therefore, depth of investigation; they have also different intrinsic limitations. Merging refracted and surface waves properly means taking benefits from both of them, and thus, deriving a better definition of the shallow subsurface parameters. Surface-wave methods are robust and offer high vertical and lateral resolution in the shallow near surface without suffering from the intrinsic limitations of refraction methods (hidden layer, velocity inversions, and others). At the same time, the investigation depth may not be sufficient for some applications, and the refraction techniques provide deeper models, giving directly the P-wave information that can be used in the data processing workflow for statics and velocity modeling, for example, without the need of conversion from V S to V P. The combination of the two techniques has, therefore, several synergies that can be exploited using simultaneous joint inversion. The joint inversion will also contribute to improve the shear-wave model reliability with Poisson s ratio anomalies. In these cases, the a priori assumption of the Poisson ratio may lead to model errors (Foti and Strobbia, 2002), and the use of the refraction data avoids this pitfall. Acknowledgments The authors acknowledge TransGlobe Energy and Dara Petroleum Company for the authorization to show the data. The authors also acknowledge WesternGeco for the authorization to publish this work. 1917
5 EDITED REFERENCES Note: This reference list is a copy-edited version of the reference list submitted by the author. Reference lists for the 2010 SEG Technical Program Expanded Abstracts have been copy edited so that references provided with the online metadata for each paper will achieve a high degree of linking to cited sources that appear on the Web. REFERENCES De Stefano, M., and D. Colombo, 2007, Pre-Stack Depth Imaging via Simultaneous Joint Inversion of Seismic, Gravity, and Magnetotelluric Data: Presented at the 69th EAGE Conference and Exhibition, EAGE. Dell Aversana, P., 2001, Integration of Seismic, MT, and Gravity Data in a Thrust Belt Interpretation: First Break, 19, no. 6, , doi: /j x. Gallardo, L. A., and M. A. Meju, 2004, Joint two-dimensional DC resistivity and Seismic travel time Inversion with Cross-Gradients constraints: Journal of Geophysical Research, 109, doi: /2003jb Hu, W., A. Abubakar, and T. M. Habashy, 2009, Joint Electromagnetic and Seismic Inversion Using Structural Constraints: Geophysics, 74, no. 6, R99 R109, doi: / Li, Y., and D. W. Oldenburg, 1996 Joint Inversion of Surface and Three Component Borehole Magnetic Data: 66th Annual International Meeting, SEG, Strobbia, C., A. Laake, P. L. Vermeer, and A. Glushchenko, 2009, Surface waves: use them, then lose them: EAGE. Foti S., C. Strobbia, 2002, Some Notes on Model Parameters for Surface Wave Data Inversion, Proceeds of SAGEEP Laake A., M. Sheneshen, C. Strobbia, L. Velasco, A. Cutts, 2010, Surface-subsurface Integration Reveals Faults in Gulf of Suez Oilfields: EAGE Virgilio, M., D. Colombo, and A. Dyke, 2008, Seismic Imaging Strategies for Thrust-Belt Exploration: Extended Offsets, Seismic/Gravity/EM Simultaneous Joint-Inversion and Anisotropic GBM Pre- Stack Depth Migration: CSPG CSEG CWLS Joint Convention. 1918
Th LHR5 08 Multi-modal Surface Wave Inversion and Application to North Sea OBN Data
Th LHR5 08 Multi-modal Surface Wave Inversion and pplication to North Sea ON Data S. Hou (CGG), D. Zheng (CGG), X.G. Miao* (CGG) & R.R. Haacke (CGG) SUMMRY Surface-wave inversion (SWI) for S-wave velocity
More information2010 SEG SEG Denver 2010 Annual Meeting
Localized anisotropic tomography with checkshot : Gulf of Mexico case study Andrey Bakulin*, Yangjun (Kevin) Liu, Olga Zdraveva, WesternGeco/Schlumberger Summary Borehole information must be used to build
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 informationRayleigh wave inversion for the near-surface characterization of shallow targets in a heavy oil field in Kuwait
Rayleigh wave inversion for the near-surface characterization of shallow targets in a heavy oil field in Kuwait Near surface complexity can be particularly challenging in the case of shallow targets, such
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 informationWe Solutions for Scattered Surface-wave Attenuation in the Western Desert of Egypt
We-14-11 Solutions for Scattered Surface-wave Attenuation in the Western Desert of Egypt D. Yanchak* (Apache Corporation), D. Monk (Apache Corporation), A.V. Zarkhidze (WesternGeco), P. Blair (WesternGeco),
More informationMAPPING POISSON S RATIO OF UNCONSOLIDATED MATERIALS FROM A JOINT ANALYSIS OF SURFACE-WAVE AND REFRACTION EVENTS INTRODUCTION
MAPPING POISSON S RATIO OF UNCONSOLIDATED MATERIALS FROM A JOINT ANALYSIS OF SURFACE-WAVE AND REFRACTION EVENTS Julian Ivanov, Choon B. Park, Richard D. Miller, and Jianghai Xia Kansas Geological Survey
More informationSuccessful application of joint reflection/refraction tomographic velocity inversion in a shallow water marine environment.
Successful application of joint reflection/refraction tomographic velocity inversion in a shallow water marine environment. Sergey Birdus 1, Dean Criddle 2, Alexey Artyomov 1, Li Li 1, Qingbing Tang 1,
More informationChapter 5. 3D data examples REALISTICALLY COMPLEX SYNTHETIC INVERSION. Modeling generation and survey design
Chapter 5 3D data examples In this chapter I will demonstrate the e ectiveness of the methodologies developed in the previous chapters using 3D data examples. I will first show inversion results of a realistically
More information(x, y, z) m 2. (x, y, z) ...] T. m 2. m = [m 1. m 3. Φ = r T V 1 r + λ 1. m T Wm. m T L T Lm + λ 2. m T Hm + λ 3. t(x, y, z) = m 1
Class 1: Joint Geophysical Inversions Wed, December 1, 29 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical, statistical Deal with data in the same
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 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 informationSummary. Introduction
Application of Surface-wave modeling and inversion in Cordova Embayment of northeastern British Columbia Antoun Salama*, Schlumberger-WesternGeco, Houston, Texas, USA ASalama@slb.com and Niranjan Banik,
More informationSEG/New Orleans 2006 Annual Meeting
3-D tomographic updating with automatic volume-based picking Dimitri Bevc*, Moritz Fliedner, Joel VanderKwaak, 3DGeo Development Inc. Summary Whether refining seismic images to evaluate opportunities in
More informationWave-equation migration from topography: Imaging Husky
Stanford Exploration Project, Report 123, October 31, 2005, pages 49 56 Short Note Wave-equation migration from topography: Imaging Husky Jeff Shragge 1 INTRODUCTION Imaging land seismic data is wrought
More informationSUMMARY. method to synthetic datasets is discussed in the present paper.
Geophysical modeling through simultaneous Joint Inversion of Seismic, Gravity and Magnetotelluric data Michele De Stefano (1), Daniele Colombo (1) WesternGeco EM - Geosystem, via Clericetti 42/A, 20133
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 informationAnisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS
Anisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS Summary Anisotropic depth model building using surface seismic data alone is non-unique and
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 informationAdaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited
Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears
More informationElastic full waveform Inversion for land walkaway VSP data
Elastic full waveform Inversion for land walkaway VSP data Olga Podgornova, Scott Leaney, arwan Charara, Schlumberger; Eric von Lunen, Nexen Energy ULC Summary Full waveform inversion (FWI) is capable
More informationGEOPHYS 242: Near Surface Geophysical Imaging. Class 8: Joint Geophysical Inversions Wed, April 20, 2011
GEOPHYS 4: Near Surface Geophysical Imaging Class 8: Joint Geophysical Inversions Wed, April, 11 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical,
More informationM. Warner* (S-Cube), T. Nangoo (S-Cube), A. Umpleby (S-Cube), N. Shah (S-Cube), G. Yao (S-Cube)
Tu A12 15 High-Resolution Reflection FWI M. Warner* (S-Cube), T. Nangoo (S-Cube), A. Umpleby (S-Cube), N. Shah (S-Cube), G. Yao (S-Cube) Summary We demonstrate reflection FWI on a less-than-ideal 3D narrow-azimuth
More informationEfficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers
Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers J.W.C. Sherwood* (PGS), K. Sherwood (PGS), H. Tieman (PGS), R. Mager (PGS) & C. Zhou
More 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 informationMain Menu. Summary. Survey Design
3D VSP acquisition and 3C processing on a deep subsalt prospect in the Gulf of Mexico John Graves, Steve Checkles, Jacques Leveille, Hess Corporation, Houston; Allan Campbell*, Scott Leaney, C. Peter Deri,
More informationGG450 4/5/2010. Today s material comes from p and in the text book. Please read and understand all of this material!
GG450 April 6, 2010 Seismic Reflection I Today s material comes from p. 32-33 and 81-116 in the text book. Please read and understand all of this material! Back to seismic waves Last week we talked about
More informationcv R z design. In this paper, we discuss three of these new methods developed in the last five years.
Nick Moldoveanu, Robin Fletcher, Anthony Lichnewsky, Darrell Coles, WesternGeco Hugues Djikpesse, Schlumberger Doll Research Summary In recent years new methods and tools were developed in seismic survey
More informationGeogiga Seismic Pro 8.0 Release Notes
Geogiga Seismic Pro 8.0 Release Notes Copyright 2015, All rights reserved. Table of Contents Introduction...1 Part 1 General Enhancements...3 Trace Display...4 Color Section Display...6 Curve Plotting...8
More informationSeismic Reflection Method
Seismic Reflection Method 1/GPH221L9 I. Introduction and General considerations Seismic reflection is the most widely used geophysical technique. It can be used to derive important details about the geometry
More informationSummary. Introduction
Chris Davison*, Andrew Ratcliffe, Sergio Grion (CGGeritas), Rodney Johnston, Carlos Duque, Jeremy Neep, Musa Maharramov (BP). Summary Azimuthal velocity models for HTI (Horizontal Transverse Isotropy)
More informationHigh definition tomography brings velocities to light Summary Introduction Figure 1:
Saverio Sioni, Patrice Guillaume*, Gilles Lambaré, Anthony Prescott, Xiaoming Zhang, Gregory Culianez, and Jean- Philippe Montel (CGGVeritas) Summary Velocity model building remains a crucial step in seismic
More informationBuilding starting model for full waveform inversion from wide-aperture data by stereotomography
Building starting model for full waveform inversion from wide-aperture data by stereotomography Vincent Prieux 1, G. Lambaré 2, S. Operto 1 and Jean Virieux 3 1 Géosciences Azur - CNRS - UNSA, France;
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 informationSeismic Time Processing. The Basis for Modern Seismic Exploration
The Future of E&P Seismic Time Processing The Basis for Modern Seismic Exploration Fusion is a leading provider of Seismic Processing for the oil and gas industry from field tapes through final migration.
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 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 informationMaximizing the value of the existing seismic data in Awali field Bahrain, by utilizing advanced 3D processing technology.
Maximizing the value of the existing seismic data in Awali field Bahrain, by utilizing advanced 3D processing technology. Eduard Maili* (OXY - Tatweer), Scott Burns (OXY), Neil Jones (Consultant, OXY)
More informationA case study for salt model building using CFP full azimuth data W. Gao*, Z. Guo, M. Guo, Q. Zhang, S. Hightower, G. Cloudy Jr. and Z.
case study for salt model building using CFP full azimuth data W. Gao*, Z. Guo, M. Guo, Q. Zhang, S. Hightower, G. Cloudy Jr. and Z. Li, TGS Summary We present a case study of the salt model building for
More informationMain Menu. Well. is the data misfit vector, corresponding to residual moveout, well misties, etc. The L matrix operator contains the d / α
Application of steering filters to localized anisotropic tomography with well data Andrey Bakulin, Marta Woodward*, Yangjun (Kevin) Liu, Olga Zdraveva, Dave Nichols, Konstantin Osypov WesternGeco Summary
More informationTh SRS3 07 A Global-scale AVO-based Pre-stack QC Workflow - An Ultra-dense Dataset in Tunisia
Th SRS3 07 A Global-scale AVO-based Pre-stack QC Workflow - An Ultra-dense Dataset in Tunisia A. Rivet* (CGG), V. Souvannavong (CGG), C. Lacombe (CGG), T. Coleou (CGG) & D. Marin (CGG) SUMMARY Throughout
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 informationMultichannel analysis of surface waves (MASW) active and passive methods
Multichannel analysis of surface waves (MASW) active and passive methods CHOON B. PARK, RICHARD D. MILLER, JIANGHAI XIA, AND JULIAN IVANOV, Kansas Geological Survey, Lawrence, USA Downloaded 01/20/13 to
More informationTu-P05-05 Multi-azimuth Anisotropic Tomography and PreSDM of a North Sea Streamer Survey
Tu-P05-05 Multi-azimuth Anisotropic Tomography and PreSDM of a North Sea Streamer Survey D. Sekulic* (ION Geophysical), O. Matveenko (Total E&P Norge), J.K. Fruehn (ION GXT) & G. Mikkelsen (Total E&P Norge)
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 informationMulti-azimuth velocity estimation
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 87 Multi-azimuth velocity estimation Robert G. Clapp and Biondo Biondi 1 ABSTRACT It is well known that the inverse problem of estimating interval
More informationAzimuthal binning for improved fracture delineation Gabriel Perez*, Kurt J. Marfurt and Susan Nissen
Azimuthal binning for improved fracture delineation Gabriel Perez*, Kurt J. Marfurt and Susan issen Abstract We propose an alternate way to define azimuth binning in Kirchhoff prestack migration. This
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 informationGeogiga Seismic Pro 9.0 Release Notes
Geogiga Seismic Pro 9.0 Release Notes Copyright 2018, All rights reserved. Table of Contents Introduction...1 Part 1 New Modules...3 Modeling2D...3 Surface RT...4 Part 2 General Enhancements...5 Part 3
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 informationMitigating Uncertainties in Towed Streamer Acquisition and Imaging by Survey Planning
Mitigating Uncertainties in Towed Streamer Acquisition and Imaging by Survey Planning M.T. Widmaier* (Petroleum Geo-Services) SUMMARY Uncertainties in seismic images or reservoir characterisation can very
More informationObstacles in the analysis of azimuth information from prestack seismic data Anat Canning* and Alex Malkin, Paradigm.
Obstacles in the analysis of azimuth information from prestack seismic data Anat Canning* and Alex Malkin, Paradigm. Summary The azimuth information derived from prestack seismic data at target layers
More informationThe determination of the correct
SPECIAL High-performance SECTION: H i gh-performance computing computing MARK NOBLE, Mines ParisTech PHILIPPE THIERRY, Intel CEDRIC TAILLANDIER, CGGVeritas (formerly Mines ParisTech) HENRI CALANDRA, Total
More informationElectromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah
Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Summary In this paper we present a new approach to the interpretation
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 informationJoint seismic traveltime and TEM inversion for near surface imaging Jide Nosakare Ogunbo*, Jie Zhang, GeoTomo LLC
Jide Nosaare Ogunbo*, Jie Zhang, GeoTomo LLC Summary For a reliable interpretation of the subsurface structure, the joint geophysical inversion approach is becoming a viable tool. Seismic and EM methods
More informationImage-guided 3D interpolation of borehole data Dave Hale, Center for Wave Phenomena, Colorado School of Mines
Image-guided 3D interpolation of borehole data Dave Hale, Center for Wave Phenomena, Colorado School of Mines SUMMARY A blended neighbor method for image-guided interpolation enables resampling of borehole
More information2012 SEG SEG Las Vegas 2012 Annual Meeting Page 1
Survey design for high-density, high-productivity, broadband vibroseis point-source and receiver acquisition a case study Peter van Baaren, Roger May, Alexander Z,arkhidze, David Morrison, and John Quigley,
More informationDeconvolution in the radial trace domain
R-T domain deconvolution Deconvolution in the radial trace domain David C. Henley ABSTRACT The radial trace (R-T) domain has been shown to be useful for coherent noise attenuation and other seismic wavefield
More informationTu A11 01 Seismic Re-processing and Q-FWI Model Building to Optimise AVO and Resolution in the Shallow Wisting Discovery
Tu A11 01 Seismic Re-processing and Q-FWI Model Building to Optimise AVO and Resolution in the Shallow Wisting Discovery G. Apeland* (WesternGeco), P. Smith (WesternGeco), O. Lewis (WesternGeco), S. Way
More informationGeometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT. August Lau and Chuan Yin.
Geometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT August Lau and Chuan Yin January 6, 2017 Abstract The goal of seismic processing is to convert input data
More informationResidual Moveout Analysis in a 3D dataset from Amplitude Variations with Offfset
P-270 Residual Moveout Analysis in a 3D dataset from Amplitude Variations with Offfset Summary Subhendu Dutta, ONGC The amplitude variations with offset (AVO) is a very useful tool which compliments the
More informationSummary. Introduction
Dmitry Alexandrov, Saint Petersburg State University; Andrey Bakulin, EXPEC Advanced Research Center, Saudi Aramco; Pierre Leger, Saudi Aramco; Boris Kashtan, Saint Petersburg State University Summary
More informationChapter N/A: The Earthquake Cycle and Earthquake Size
This sheet provides a general outline of the topics covered on your upcoming exam. This is not an exhaustive list, but is simply a general list of key points. In other words, most exam questions will be
More information3D inversion of marine CSEM data: A feasibility study from the Shtokman gas field in the Barents Sea
3D inversion of marine CSEM data: A feasibility study from the Shtokman gas field in the Barents Sea M. S. Zhdanov 1,2, M. Čuma 1,2, A. Gribenko 1,2, G. Wilson 2 and N. Black 2 1 The University of Utah,
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 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 informationMinimizing Fracture Characterization Uncertainties Using Full Azimuth Imaging in Local Angle Domain
P-237 Minimizing Fracture Characterization Uncertainties Using Full Azimuth Imaging in Local Angle Domain Shiv Pujan Singh*, Duane Dopkin, Paradigm Geophysical Summary Shale plays are naturally heterogeneous
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 informationIntegral equation method for anisotropic inversion of towed streamer EM data: theory and application for the TWOP survey
Integral equation method for anisotropic inversion of towed streamer EM data: theory and application for the TWOP survey Michael S. Zhdanov 1,2, Masashi Endo 1, Daeung Yoon 1,2, Johan Mattsson 3, and Jonathan
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 informationGravity Gradients in PreStack Depth Migration
Index Table of contents Gravity Gradients in PreStack Depth Migration Ed. K. Biegert Shell International Exploration and Production Summary Subsurface de-risking is presently almost exclusively done by
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 informationG042 Subsalt Imaging Challenges - A Deepwater Imaging Analysis
G042 Subsalt Imaging Challenges - A Deepwater Imaging Analysis M. Cogan* (WesternGeco), J. Gardner (WesternGeco) & N. Moldoveanu (WesternGeco) SUMMARY Upon completion of the final reverse-time migration
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 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 informationMASW Horizontal Resolution in 2D Shear-Velocity (Vs) Mapping
MASW Horizontal Resolution in 2D Shear-Velocity (Vs) Mapping by Choon B. Park Kansas Geological Survey University of Kansas 1930 Constant Avenue, Campus West Lawrence, Kansas 66047-3726 Tel: 785-864-2162
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 informationZero offset VSP processing for coal reflections
Zero offset VSP processing for coal reflections Salman K. Bubshait and Don C. Lawton ABSTRACT Shlumberger has acquired five VSP surveys on behalf of EnCana in Alberta. The primary goal of the VSP surveys
More informationGeogiga Seismic Pro 8.3 Release Notes
Geogiga Seismic Pro 8.3 Release Notes Copyright 2017, All rights reserved. Table of Contents Introduction...1 Part 1 Utility Modules...2 Part 2 Reflection Modules...4 Updates in SF Imager...5 Updates in
More information3-D vertical cable processing using EOM
Carlos Rodriguez-Suarez, John C. Bancroft, Yong Xu and Robert R. Stewart ABSTRACT Three-dimensional seismic data using vertical cables was modeled and processed using equivalent offset migration (EOM),
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 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 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 information=, (1) SEG/New Orleans 2006 Annual Meeting
U. Albertin* +, P. Sava ++, J. Etgen +, and M. Maharramov + + BP EPTG, Houston, Texas, ++ Colorado School of Mines, Goldin, Colorado Summary A methodology for velocity updating using one-way wavefield
More 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 informationCommon-angle processing using reflection angle computed by kinematic pre-stack time demigration
Common-angle processing using reflection angle computed by kinematic pre-stack time demigration Didier Lecerf*, Philippe Herrmann, Gilles Lambaré, Jean-Paul Tourré and Sylvian Legleut, CGGVeritas Summary
More informationHIGH RESOLUTION STACKING OF SEISMIC DATA. Marcos Ricardo Covre, Tiago Barros and Renato da Rocha Lopes
HIGH RESOLUTION STACKING OF SEISMIC DATA Marcos Ricardo Covre, Tiago Barros and Renato da Rocha Lopes School of Electrical and Computer Engineering, University of Campinas DSPCom Laboratory, Department
More informationWe LHR5 06 Multi-dimensional Seismic Data Decomposition for Improved Diffraction Imaging and High Resolution Interpretation
We LHR5 06 Multi-dimensional Seismic Data Decomposition for Improved Diffraction Imaging and High Resolution Interpretation G. Yelin (Paradigm), B. de Ribet* (Paradigm), Y. Serfaty (Paradigm) & D. Chase
More informationAdvances in radial trace domain coherent noise attenuation
Advances in radial trace domain coherent noise attenuation ABSTRACT David C. Henley* CREWES, Department of Geology and Geophysics University of Calgary, Calgary, AB henley@crewes.org The radial trace transform,
More information2D Inversions of 3D Marine CSEM Data Hung-Wen Tseng*, Lucy MacGregor, and Rolf V. Ackermann, Rock Solid Images, Inc.
2D Inversions of 3D Marine CSEM Data Hung-Wen Tseng*, Lucy MacGregor, and Rolf V. Ackermann, Rock Solid Images, Inc. Summary A combination of 3D forward simulations and 2D and 3D inversions have been used
More informationSelection of an optimised multiple attenuation scheme for a west coast of India data set
P-391 Selection of an optimised multiple attenuation scheme for a west coast of India data set Summary R Pathak*, PC Kalita, CPS Rana, Dr. S. Viswanathan, ONGC In recent years a number of new algorithms
More informationClosing the Loop via Scenario Modeling in a Time-Lapse Study of an EOR Target in Oman
Closing the Loop via Scenario Modeling in a Time-Lapse Study of an EOR Target in Oman Tania Mukherjee *(University of Houston), Kurang Mehta, Jorge Lopez (Shell International Exploration and Production
More informationPS wave AVO aspects on processing, inversion, and interpretation
PS wave AVO aspects on processing, inversion, and interpretation Yong Xu, Paradigm Geophysical Canada Summary Using PS wave AVO inversion, density contrasts can be obtained from converted wave data. The
More informationGEOPHYS 242: Near Surface Geophysical Imaging. Class 4: First-Arrival Traveltime Tomography Mon, April 11, 2011
GEOPHYS 242: Near Surface Geophysical Imaging Class 4: First-Arrival Traveltime Tomography Mon, April 11, 2011 Wavefront tracing methods - speed versus accuracy Inversion algorithms - any magic approach?
More information2011 SEG SEG San Antonio 2011 Annual Meeting 3938
Depth imaging Coil data: Multi azimuthal tomography earth model building and depth imaging the full azimuth Tulip coil project Michele Buia 1, Peter Brown 2, Bakhrudin Mansyur 2, Michelle Tham 3, Suyang
More informationAnatomy of common scatterpoint (CSP) gathers formed during equivalent offset prestack migration (EOM)
Anatomy of CSP gathers Anatomy of common scatterpoint (CSP) gathers formed during equivalent offset prestack migration (EOM) John C. Bancroft and Hugh D. Geiger SUMMARY The equivalent offset method of
More informationFoolproof AvO. Abstract
Foolproof AvO Dr. Ron Masters, Geoscience Advisor, Headwave, Inc Copyright 2013, The European Association of Geoscientists and Engineers This paper was prepared for presentation during the 75 th EAGE Conference
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 informationDaniel R.H. O Connell, Ph.D. Geophysicist. Fugro Consultants, Inc. 25 March Diablo Canyon SSHAC SSC Working Meeting March 2014
2011-2012 AB1632 Onshore Seismic Data Daniel R.H. O Connell, Ph.D. Geophysicist Fugro Consultants, Inc. 25 March 2014 Diablo Canyon SSHAC SSC Working Meeting March 2014 2011 and 2012 CCCSIP Onshore Seismic
More information