Temporal Integration Of Seismic Traveltime Tomography
|
|
- Nigel Smith
- 6 years ago
- Views:
Transcription
1 Temporal Integration Of Seismic Traveltime Tomography Jonathan B Ajo-Franklin, Earth Resources Laboratory Dept of Earth, Atmospheric, and Planetary Sciences Massachusetts Institute of Technology Cambridge, MA 2142, Jaime Urban, Jerry M Harris, Stanford University June 1, 5 Abstract Time-lapse geophysical measurements and seismic imaging methods in particular are powerful techniques for monitoring changes in reservoir properties Traditional time-lapse processing methods treat each dataset as an independent unit and estimate changes in reservoir state through differencing these separate inversions We present a general least-squares approach to jointly inverting time-varying property models through use of spatio-temporal coupling operators Originally developed within the medical imaging community, this extension of traditional Tikhonov regularization allows us to constrain the way in which models vary in time, thereby reducing artifacts observed in traditional time-lapse imaging formulations The same methodology can also accommodate changes in experiment geometry as a function of time thus allowing inversion of incremental or incomplete surveys In this case, temporal resolution is traded for improved spatial coverage at individual timesteps We use seismic traveltime tomography as a model problem although almost any geophysical inversion task can be posed within this formalism We apply the developed time-lapse inversion algorithm to a synthetic crosswell dataset designed to replicate a CO 2 sequestration monitoring experiment 1 Introduction Time-lapse geophysical measurements provide a powerful toolbox of techniques for monitoring subsurface flow processes including the enhanced recovery of oil and gas (6), CO 2 sequestration (2) (5), and contaminant remediation (4) Seismic methods in particular are increasingly used to generate quantitative maps of variations in fluid saturation or pressure over spatial domains ranging from 1s of meters for environmental problems to 1s of kilometers for field-wide production monitoring Most traditional time-lapse processing techniques treat each acquired dataset as an independent measurement Temporal variations in subsurface properties are typically determined by subtraction of image pairs Naive image subtraction tends to be sensitive to survey-to-survey changes in S/N ratio and variations in acquisition geometry, both of which can generate artifacts in the resulting time-lapse images This conceptual approach also fails to acknowledge a key component of prior knowledge, mainly that the changes in subsurface seismic response in these environments are due to variations in fluid properties or effective stress and not changes in geologic structure Significant artifacts in time-lapse processing results have spurred development of cross-equalization algorithms (8) designed to match the geometry and signal characteristics of repeated surveys Joint inversion of a series of time-lapse surveys seems to be a natural approach to decreasing image artifacts which result from naive differencing Through use of constraints on how the model can change in time, we can suppress non-repeatable noise while providing a more consistent reconstruction of real variations 1
2 in subsurface properties Day-Lewis etal(3) demonstrated a technique for integrating time-lapse imaging based on using smooth basis functions with a time dimension for tomographic inversion They effectively applied their method to dynamically image the flow of a saline tracer through a large aperture fracture using borehole radar We adopt an approach developed within the medical imaging community based on a straight-forward extension of Tikhonov regularization to include a time dimension Brooks etal (1) apply a temporal constraint operator within a least-squares framework for the inversion of electrocardiography data Zhang etal (9) perform a comparison of spatio-temporal inversion approaches including both the Tikhonov approach and Kalman filtering We largely follow the formulation of Brooks etal with the exception of our choice of regularization operator and the techniques used for solving the coupled system 2 Linear Formulation We will initially consider the general linear inverse problem where a linear operator G, maps a model (m) to a dataset (d), G m = d (1) Consider a series of n datasets (d i ) acquired at different times (t i ), possibly with different source/receiver geometries as encapsulated by variations in the kernel (G i ) but equivalent model parameterizations In the naive case, the inverse problem would be solved independently for each dataset by minimizing an objective function of the form Φ(m 1,, m n ) = G i m i d i λ 2 s Dm i 2 2 (2) where D is a weighting operator and λ s is a regularization parameter The first term in equation 2 measures data misfit while the second measures model length as modified by D If D = I than the resulting minimization of equation 2 is simply the damped least-squares solution Neither of the terms in equation 2 couple solutions across multiple time steps since modification to m i does not effect either the data misfit or length of model m j (i j); this approach seems most applicable in cases where models are not correlated in time In the case of our time-lapse seismic problem, models have a strong temporal correlation since the underlying geologic structure is clearly not changing in the relatively short (geologically speaking) time between repeat surveys We can modify equation 2 to include a temporal cross-coupling term which minimizes the time-lapse change in some model attribute in addition to decreasing data misfit and model norm for individual surveys Consider a combined objective function of the form Φ(m 1,, m n ) = G i m i d i (3) λ 2 s n 1 λ 2 t D m i D m i+1 D m i 2 2 t i where t i = t i+1 t i Minimizing equation 3 is equivalent to the least-squares solution of an augmented 2
3 system with operator G c, G c = G 1 G 2 G n λ s D λ s D (4) λ s D λ t D λ t D t 1 t 1 λ t D λ t D t n 1 t n 1 d 1 d 2 m 1 m 2 G c = d n (5) m n The numerator in the third term of equation 3 penalizes the differences between models which are temporal neighbors while the denominator scales this weight by the time elapsed between surveys The second regularization parameter, λ t, controls the strength of the temporal constraint If D = I this component can be thought of as a time damping term For the simple case of two surveys (n = 2) and th order Tikhonov regularization for both the spatial and spatio-temporal terms (D = I), equation 4 can be reduced to G 1 G 2 λ s I λ s I λ t t I λt t I [ m1 Although damping is the simplest type of temporal regularization to add, different types of operators should also be considered Minimizing change in the first or second spatial derivatives of the model are two alternatives We can imagine computing the associated pseudo-inverse inverse of our augmented operator, G g c m 2 ], which allows us to write the model resolution matrix as, R = G g c 3 = d 1 d 2 (6) G c where R can be seen as a filter
4 which shows how the imaging experiment and the choice of G g c modify the true model In this case R gives us insight into not only the spatial aspects of model resolution but also shows us temporal smearing, the process by which information is spread between our multiple experiments This inversion formalism also provides a good approach to incrementally acquired seismic surveys where the survey geometry at any particular time step n is relatively sparse Although independent inversion of a single survey might yield an image with very low spatial resolution, by jointly inverting a series of surveys we can effectively add spatial aperture in exchange for losing temporal resolution 3 Traveltime Tomography Up to this point, our formulation has been relatively general with no assumption regarding the operation which G performs, the model parametrization represented by m, or the type of data stored as d We will now apply our formulation to the concrete example of seismic traveltime tomography with one temporal dimension and two spatial dimensions In this case we choose each m to be a rectalinear mesh of homogeneous slowness cells while d is a vector of picked first-arrival traveltimes and G is the ray-path matrix For the examples presented in this paper, we will use straight rays to preserve the linearity of the coupled inverse problem but extensions to the non-linear case seem quite feasible We use a split laplacian operator (D xx, D zz ) to allow anisotropic regularization with two spatial parameters (λ sx, λ sz ) and two spatio-temporal parameters (λ tx, λ tz ) for the respective terms in equation 4 Regularization parameters are chosen by observation although use of the L-surface technique advocated by Brooks etal (1) would decrease the amount of manual tuning required in the inversion process The resulting coupled systems were solved using the LSQR algorithm (7) 31 A Synthetic CO 2 Monitoring Problem Several field experiments to date have attempted to monitor the subsurface extent of active CO 2 injections including the McElroy (5) and Weyburn (2) EOR projects For our synthetic experiment we have generated four time-lapse images of a CO 2 flood progressing through a permeable layer shown in row (A) of figure 1 Data was synthesized for a crosswell geometry with sources and receivers evenly spaced near the boundary of the model; the resulting arrays allowed tomographic imaging within a x 98 m subdomain of the 115 x 49 m initial model Gaussian noise ( 3%) was added to the traveltime picks for all four synthetic surveys The lower 3 rows of figure 1 show the inversion results for the uncoupled case (B) and the effect of including the spatio-temporal coupling term (C and D) using a 3 x 7 cell mesh In all three cases the same spatial regularization parameters were used [λ sx = 1, λ sz = 8] while the weak and strongly coupled examples used spatio-temporal coupling parameters of [λ tx = 16, λ tz = 12] and [λ tx = 64, λ tz = 48] respectively While all three inversions successfully imaged the target CO 2 flood, visible at approximately 7 m depth, stronger spatio-temporal constraints appear to suppress velocity artifacts Figure 2 shows the result of differencing the tomograms between the four sequential surveys with the same ordering of conditions Similar trends in noise levels are visible in the difference images In the uncoupled case (row B) the artifacts are of the same order of magnitude as the CO 2 flooded region making confident interpretation difficult The use of strong spatio-temporal coupling clearly decreases this interpretation ambiguity at the expense of accurate difference velocity estimates and spatial resolution 4 Conclusion Temporal regularization provides a natural approach for jointly inverting time-lapse datasets A key topic of future research is understanding the trade-offs between spatial and temporal resolution displayed within our synthetic example While we have not discussed the non-linear tomography problem in this paper, the same techniques could be easily integrated within iterative imaging algorithms Recent developments in the medical imaging community suggest several routes to improving these methods including the use 4
5 T1 T2 T3 T4 (A) (B) (C) (D) X (m) X (m) X (m) X (m) Vp (m/s) Figure 1: Tomography results with variable temporal coupling : (A) True model used to generate synthetics, (B) Uncoupled, (C) Weak temporal coupling, (D) Strong temporal coupling (regularization parameters in text) 5
6 T1 d1 T2 d2 T3 d3 T4 (A) (B) (C) (D) X(m) X(m) X(m) Vp (m/s) 5 Figure 2: Difference tomograms with variable temporal coupling corresponding to the tomograms in figure 1 : (A) True model used, (B) Uncoupled tomograms, (C) Weak temporal coupling, (D) Strong temporal coupling 6
7 of L-surface techniques for regularization parameter selection (1) and the incorporation of prior temporal covariance models (9) Our future work will attempt to integrate and extend these techniques to the seismic monitoring problem within the context of full wave-equation tomography References [1] Brooks, D, Ahmad, G, MacLeod, R, and Maratos, G Inverse electrocardiography by simultaneous imposition of multiple constraints IEEE Transactions On Biomedical Engineering 46, 1 (January 1999), 3 18 [2] Davis, T, Terrel, M, R B, Cardona, R, Kendall, R, and Winarsky, R Multicomponent seismic characterization and monitoring of the CO 2 flood at Weyburn Field, Saskatchewan Leading Edge 22 (July 3), [3] Day-Lewis, F, Harris, J, and Gorelick, S Time-lapse inversion of crosswell radar data Geophysics 67, 6 (November-December 2), [4] Hubbard, S, Chen, J, Peterson, J, Majer, E, Williams, K, Swift, D, Mailloux, B, and Rubin, Y Hydrogeological characterization of the South Oyster bacterial transport site using geophysical data Water Resources Research 37, 1 (October 1), [5] Lazaratos, S, and Marion, B Crosswell seismic imaging of reservoir changes caused by CO 2 injection The Leading Edge 16 (September 1997), [6] Lumley, D Time-lapse seismic reservoir monitoring Geophysics 66, 1 (January-February 1), 5 53 [7] Paige, C, and Saunders, M An algorithm for sparse linear equations and sparse least squares ACM Transactions in Mathematical Software 8 (1982), [8] Rickett, J, and Lumley, D Cross-equalization data processing for time-lapse seismic reservoir monitoring: A case study from the Gulf of Mexico Geophysics 66, 4 (1), [9] Zhang, Y, Ghodrati, A, and Brooks, D An analytical comparison of three spatio-temporal regularization methods for dynamic linear inverse problems in a common statistical framework Inverse Problems 21 (5),
Wave-equation MVA applied to 4-D seismic monitoring
Stanford Exploration Project, Report 112, November 11, 2002, pages 15 21 Short Note Wave-equation MVA applied to 4-D seismic monitoring Paul Sava, John Etgen, and Leon Thomsen 1 INTRODUCTION 4-D seismic
More informationP257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies
P257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies. azemi* (University of Alberta) & H.R. Siahkoohi (University of Tehran) SUMMARY Petrophysical reservoir properties,
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 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 informationIssues During the Inversion of Crosshole Radar Data: Can We Have Confidence in the Outcome?
Boise State University ScholarWorks Geosciences Faculty Publications and Presentations Department of Geosciences 12-1-2006 Issues During the Inversion of Crosshole Radar Data: Can We Have Confidence in
More informationA Data estimation Based Approach for Quasi continuous Reservoir Monitoring using Sparse Surface Seismic Data Introduction Figure 1
A Data estimation Based Approach for Quasi continuous Reservoir Monitoring using Sparse Surface Seismic Data Adeyemi Arogunmati* and Jerry M. Harris, Stanford University, California, USA Introduction One
More informationTime-Lapse Seismic Crosswell Monitoring of CO 2 injected in an Onshore Sandstone Aquifer
Time-Lapse Seismic Crosswell Monitoring of CO 2 injected in an Onshore Sandstone Aquifer Jesper Spetzler (1), Zigiu Xue (2), Hideki Saito (3), Dai Nobuoka (3), Hiroyuki Azuma (3), and Osamu Nishizawa (4)
More informationInterval velocity estimation through convex optimization
Stanford Exploration Project, Report 125, January 16, 2007, pages 1?? Interval velocity estimation through convex optimization Ben Witten and Michael Grant ABSTRACT Convex optimization is an optimization
More 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 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 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 informationA MAP Algorithm for AVO Seismic Inversion Based on the Mixed (L 2, non-l 2 ) Norms to Separate Primary and Multiple Signals in Slowness Space
A MAP Algorithm for AVO Seismic Inversion Based on the Mixed (L 2, non-l 2 ) Norms to Separate Primary and Multiple Signals in Slowness Space Harold Ivan Angulo Bustos Rio Grande do Norte State University
More informationChapter 1. Introduction
Chapter 1 Introduction Time-lapse (4D) seismic imaging has become an established technology for monitoring changes in subsurface reservoir properties. In general, time-lapse seismic imaging involves repetition
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 informationTRAVELTIME TOMOGRAPHY (CONT)
30 March 005 MODE UNIQUENESS The forward model TRAVETIME TOMOGRAPHY (CONT di = A ik m k d = Am (1 Data vecto r Sensitivit y matrix Model vector states the linearized relationship between data and model
More informationUltrasonic Multi-Skip Tomography for Pipe Inspection
18 th World Conference on Non destructive Testing, 16-2 April 212, Durban, South Africa Ultrasonic Multi-Skip Tomography for Pipe Inspection Arno VOLKER 1, Rik VOS 1 Alan HUNTER 1 1 TNO, Stieltjesweg 1,
More informationTimelapse ERT inversion approaches and their applications
Timelapse ERT inversion approaches and their applications THOMAS GÜNTHER Leibniz-Institute for Applied Geophysics (LIAG), Hannover, Germany. thomas.guenther@liag-hannover.de Introduction Aim of timelapse
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 informationIterative resolution estimation in Kirchhoff imaging
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 369?? Iterative resolution estimation in Kirchhoff imaging Robert G. Clapp, Sergey Fomel, and Marie Prucha 1 ABSTRACT We apply iterative
More informationFracture Quality from Integrating Time-Lapse VSP and Microseismic Data
Fracture Quality from Integrating Time-Lapse VSP and Microseismic Data Mark E. Willis, Daniel R. Burns, Rongrong Lu, M. Nafi Toksöz, Earth Resources Laboratory Dept. of Earth, Atmospheric, and Planetary
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 informationTotal variation tomographic inversion via the Alternating Direction Method of Multipliers
Total variation tomographic inversion via the Alternating Direction Method of Multipliers Landon Safron and Mauricio D. Sacchi Department of Physics, University of Alberta Summary Geophysical inverse problems
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 informationA Scalable Parallel LSQR Algorithm for Solving Large-Scale Linear System for Seismic Tomography
1 A Scalable Parallel LSQR Algorithm for Solving Large-Scale Linear System for Seismic Tomography He Huang, Liqiang Wang, Po Chen(University of Wyoming) John Dennis (NCAR) 2 LSQR in Seismic Tomography
More informationSUMMARY THEORY INTRODUCTION
Double-Difference Waveform Inversion of 4D Ocean Bottom Cable Data: Application to Valhall, North Sea Di Yang, Michael Fehler and Alison Malcolm, MIT, Faqi Liu and Scott Morton, Hess Corporation SUMMARY
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 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 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 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 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 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 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 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 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 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 informationThe Conventional and Non-conventional Seismic Differencing
Summary The Conventional and Non-conventional Seismic Differencing Vanja Vracar* University of Calgary, CREWES Project, Alberta, Canada vmilicev@ucalgary.ca and Robert J. Ferguson University of Calgary,
More informationSimultaneous joint inversion of refracted and surface waves Simone Re *, Claudio Strobbia, Michele De Stefano and Massimo Virgilio - WesternGeco
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
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 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 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 informationUsing Resolution-Constrained Adaptive Meshes for Traveltime Tomography
Using Resolution-Constrained Adaptive Meshes for Traveltime Tomography Jonathan B. Ajo-Franklin Earth Resources Laboratory, Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute
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 informationP063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA
1 P063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA I. TRINKS 1, S. SINGH 2, C. CHAPMAN 3, P. BARTON 1, M. BOSCH 4, A. CHERRETT 5 Addresses 1 Bullard Laboratories,
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 informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Aliasing Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 16, 2014 Bee Bednar (Panorama Technologies) Seismic
More 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 informationIPTC Title: High Resolution, Simultaneous VSP and Land Seismic Acquisition. Scott Robinson, Qatar Petroleum
Title: High Resolution, Simultaneous VSP and Land Seismic Acquisition Authors Scott Robinson, Qatar Petroleum primary author o P. O. Box 47, Doha, Qatar; (+974) 440-1465; s_robinson@qp.com.qa Salva R Seeni,
More informationRobustness of the scalar elastic imaging condition for converted waves
CWP-830 Robustness of the scalar elastic imaging condition for converted waves Yuting Duan & Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT For elastic reverse-time migration, one
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 informationSimultaneous source separation via robust time variant Radon operators
Simultaneous source separation via robust time variant Radon operators Summary Amr Ibrahim, University of Alberta, Edmonton, Canada amr.ibrahim@ualberta.ca and Mauricio D. Sacchi, University of Alberta,
More informationIntroduction to the traveltime tomography and prestack migration for borehole-data within REFLEXW
Introduction to the traveltime tomography and prestack migration for borehole-data within REFLEXW In the following the use of the traveltime tomography and the prestack-migration for different types of
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 information4D Seismic Inversion on Continuous Land Seismic Reservoir Monitoring of Thermal EOR
4D Seismic Inversion on Continuous Land Seismic Reservoir Monitoring of Thermal EOR Laurene Michou, CGGVeritas, Massy, France, laurene.michou@cggveritas.com Thierry Coleou, CGGVeritas, Massy, France, thierry.coleou@cggveritas.com
More informationSparse wavelet expansions for seismic tomography: Methods and algorithms
Sparse wavelet expansions for seismic tomography: Methods and algorithms Ignace Loris Université Libre de Bruxelles International symposium on geophysical imaging with localized waves 24 28 July 2011 (Joint
More informationSound-speed tomography using first-arrival transmission ultrasound for a ring array
Sound-speed tomography using first-arrival transmission ultrasound for a ring array Youli Quan* a,b and Lianjie Huang b a Department of Geophysics, Stanford University, Stanford, CA 9435-2215 b Mail Stop
More informationMachine-learning Based Automated Fault Detection in Seismic Traces
Machine-learning Based Automated Fault Detection in Seismic Traces Chiyuan Zhang and Charlie Frogner (MIT), Mauricio Araya-Polo and Detlef Hohl (Shell International E & P Inc.) June 9, 24 Introduction
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 informationTravel Time Inversion in Seismic Tomography
Travel Time Inversion in Seismic Tomography S. Padina, D. Churchill, R. P. Bording March 2, 2006 1 Introduction The word tomography finds its origins in the two Greek words tomos, meaning section, and
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 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 informationAdaptive spatial resampling as a Markov chain Monte Carlo method for uncertainty quantification in seismic reservoir characterization
1 Adaptive spatial resampling as a Markov chain Monte Carlo method for uncertainty quantification in seismic reservoir characterization Cheolkyun Jeong, Tapan Mukerji, and Gregoire Mariethoz Department
More informationAmplitude-preserved low velocity noise suppression
Stanford Exploration Project, Report 84, May 9, 2001, pages 197 213 Amplitude-preserved low velocity noise suppression Hector Urdaneta 1 ABSTRACT Low velocity noise can contaminate reflections and distort
More informationLeast-squares migration/inversion of blended data
Least-squares migration/inversion of blended data Yaxun Tang and Biondo Biondi ABSTRACT We present a method based on least-squares migration/inversion to directly image data collected from recently developed
More informationA new take on FWI: Wavefield Reconstruction Inversion
A new take on FWI: Wavefield Reconstruction Inversion T. van Leeuwen 1, F.J. Herrmann and B. Peters 1 Centrum Wiskunde & Informatica, Amsterdam, The Netherlands University of British Columbia, dept. of
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 informationMulticomponent land data pre-processing for FWI: a benchmark dataset
Multicomponent land data pre-processing for FWI: a benchmark dataset Raul Cova, Bernie K. Law and Kris Innanen CRWES/University of Calgary Summary Successful full-waveform inversion (FWI) studies using
More informationIntroduction to Geophysical Modelling and Inversion
Introduction to Geophysical Modelling and Inversion James Reid GEOPHYSICAL INVERSION FOR MINERAL EXPLORERS ASEG-WA, SEPTEMBER 2014 @ 2014 Mira Geoscience Ltd. Forward modelling vs. inversion Forward Modelling:
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 informationPTE 519 Lecture Note Finite Difference Approximation (Model)
PTE 519 Lecture Note 3 3.0 Finite Difference Approximation (Model) In this section of the lecture material, the focus is to define the terminology and to summarize the basic facts. The basic idea of any
More information3D Inversion of Time-Domain Electromagnetic Data for Ground Water Aquifers
3D Inversion of Time-Domain Electromagnetic Data for Ground Water Aquifers Elliot M. Holtham 1, Mike McMillan 1 and Eldad Haber 2 (1) Computational Geosciences Inc. (2) University of British Columbia Summary
More informationTime-Lapse Full Waveform Inversion Methods. Wubshet M. Alemie. A thesis submitted in partial fulfillment of the requirements for the degree of
Time-Lapse Full Waveform Inversion Methods by Wubshet M. Alemie A thesis submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Geophysics Department of Physics
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 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 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 informationWe G Updating the Reservoir Model Using Engineeringconsistent
We G102 09 Updating the Reservoir Model Using Engineeringconsistent 4D Seismic Inversion S. Tian* (Heriot-Watt University), C. MacBeth (Heriot-Watt University) & A. Shams (Heriot-Watt University) SUMMARY
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 informationTarget-oriented wave-equation inversion with regularization in the subsurface-offset domain
Stanford Exploration Project, Report 124, April 4, 2006, pages 1?? Target-oriented wave-equation inversion with regularization in the subsurface-offset domain Alejandro A. Valenciano ABSTRACT A complex
More informationFull waveform inversion of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC)
of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC) Summary (FWI) is a promising technology for next generation model building. However, it faces with many
More informationLecture 2: Seismic Inversion and Imaging
Lecture 2: Seismic Inversion and Imaging Youzuo Lin 1 Joint work with: Lianjie Huang 1 Monica Maceira 2 Ellen M. Syracuse 1 Carene Larmat 1 1: Los Alamos National Laboratory 2: Oak Ridge National Laboratory
More informationTime-lapse inversion of crosswell radar data
GEOPHYSICS, VOL. 67, NO. 6 (NOVEMBER-DECEMBER 2002); P. 1740 1752, 7 FIGS., 2 TABLES. 10.1190/1.1527075 Time-lapse inversion of crosswell radar data Frederick D. Day-Lewis, Jerry M. Harris, and Steven
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 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 informationEfficient Double-Beam Characterization for Fractured Reservoir. Yingcai Zheng, Xinding Fang, Michael C. Fehler and Daniel R. Burns
Efficient Double-Beam Characterization for Fractured Reservoir Yingcai Zheng, Xinding Fang, Michael C. Fehler and Daniel R. Burns Abstract We proposed an efficient target-oriented method to characterize
More informationReverse time migration in midpoint-offset coordinates
Stanford Exploration Project, Report 111, June 9, 00, pages 19 156 Short Note Reverse time migration in midpoint-offset coordinates Biondo Biondi 1 INTRODUCTION Reverse-time migration (Baysal et al., 198)
More 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 informationInterpolation with pseudo-primaries: revisited
Stanford Exploration Project, Report 129, May 6, 2007, pages 85 94 Interpolation with pseudo-primaries: revisited William Curry ABSTRACT Large gaps may exist in marine data at near offsets. I generate
More 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 informationFive Dimensional Interpolation:exploring different Fourier operators
Five Dimensional Interpolation:exploring different Fourier operators Daniel Trad CREWES-University of Calgary Summary Five-Dimensional interpolation has become a very popular method to pre-condition data
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 informationSeisSpace Software. SeisSpace enables the processor to be able focus on the science instead of being a glorified data manager.
SeisSpace Software OVERVIEW Delivering large volumes of data quickly and accurately remains the hallmark of any seismic processing system. With intuitive analysis tools, state-of-the-art geophysical algorithms
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 informationJoint least-squares migration/inversion: theory
Chapter 3 Joint least-squares migration/inversion: theory This chapter presents inversion methods that overcome some short-comings of conventional time-lapse imaging methods. In the previous chapter, I
More informationFull-waveform inversion for reservoir characterization: A synthetic study
CWP-889 Full-waveform inversion for reservoir characterization: A synthetic study Nishant Kamath, Ilya Tsvankin & Ehsan Zabihi Naeini ABSTRACT Most current reservoir-characterization workflows are based
More informationFlattening without picking
Stanford Exploration Project, Report 112, November 11, 2002, pages 141 151 Flattening without picking Jesse Lomask and Jon Claerbout 1 ABSTRACT We introduce an analytical method for integrating dip information
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 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 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 informationWe G High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation
We G103 05 High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation G. Hilburn* (TGS), Y. He (TGS), Z. Yan (TGS) & F. Sherrill (TGS) SUMMARY An approach
More informationA simulated Simultaneous Source Experiment in Shallow waters and the Impact of Randomization Schemes Rolf Baardman, Roald van Borselen, PGS
simulated Simultaneous Source Experiment in Shallow waters and the Impact of Randomization Schemes Rolf aardman, Roald van orselen, PGS Summary In simultaneous source acquisition, seismic data can be recorded
More informationProgress Report on: Interferometric Interpolation of 3D SSP Data
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
More information