TRAVELTIME TOMOGRAPHY (CONT)
|
|
- Roxanne Lloyd
- 5 years ago
- Views:
Transcription
1 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 via a mapping or sensitivity matrix, where i = 1,,, N (data and k = 1,,, M (model parameters. In the case of GRT traveltime tomography, A can be thought of as a sparse sampling matrix whose nonzero elements are the lengths of segments of the raypath within each discretized voxel of the Earth. l 11 l 1M ray 1 A M O M M l N 1 l NM ray N Normally, N >> M the number of data (typically ~ is far greater than the number of model parameters (typically ~ making (1 nominally overdetermined. Consequently, A is normally not square. We can invert for m by first left-multiplying (1 by A T T T A Am = A d and then, since A T A is square, we complete the inversion to yield the classic least squares solution, as follows: T 1 T m = ( A A A d. ( The penalty function that is satisfied by ( is ε m ( T 1 T A A d A (3 where m (plugged in for m in (3 minimizes ε in the least squares sense. Aside: penalty metrics other than the squared norm in (3 can be used; the absolute norm, for example. As stated above, N >> M Æ the system is nominally overdetermined. But A is often rank-deficient mathematically, rank(a << M and therefore the system behaves as though it is underdetermined. In other words, many of the N observations are linearly dependent and so do not provide complementary information on the subsurface. To understand this, consider the worst-case scenario where we have 10 6 observations that travel along the exact same raypath like the figure shown at right. Only cells traversed by the ray are actually sampled, meaning that unsampled cells are completely unconstrained. From an inverse point of view, the unconstrained parameter cells make the problem underdetermined there are more parameters than S R 1
2 30 March 005 there are complementary pieces of information to constrain the system. When an inverse problem is underdetermined, all the following statements are true and/or synonymous: {1} the problem is ill-posed; {} A T A is singular/non-invertible/rank-deficient; {3} one or more rows or columns of A and A T A are linearly dependent; {4} one or more parameters are unconstrained; {5} A and A T A both have nontrivial null spaces. All parts of a discretized model that are unconstrained or, equivalently, that are not mapped by A from model space into data space reside in the null space of A. This means that ε is independent of unconstrained parameters; an unconstrained parameter can take on any value without affecting ε. So the least squares solution is not unique if one or more parameters are unconstrained! In fact, there is an infinite set of solutions that equally minimize ε. Namely, any prospective solution that contains the true solution (which we don t know linearly combined with any of the vectors in null (A will equally minimize ε. This is termed nonuniqueness. REGUARIZATION AND DAMPING Model regularization essentially adds additional constraints to the inverse problem to address the solution nonuniqueness that arises as a consequence of ill-posedness. For example, the expression a + b = is ill-posed if we are asked to find values for both a and b. There are an infinite number of combinations of a and b that equally satisfy this expression, so the inverse problem is nonunique. But if we add an extra condition or constraint that specifies λ 1 ( a b = 0, where λ 1 is a user-specified weight (a.k.a. agrange multiplier, we now have a well-posed system of two equations in two unknowns (provided λ 1 0. The unique solution to this system is a = b = 1. This second expression is a form of model regularization. It assumes that a functional relationship exists between the model parameters (in this case, it basically assumes that a = b that can be used to add information and thereby constrain the original problem. This idea can be applied to (3 by augmenting it with an additional regularization condition: ( T ε= m A A 1 A T d + λ 1 m (4 where, for example, can defined as the gradient or first-difference operator M O O O O that basically reduces the difference between adjacent model parameters in m. This type of constraint requires that neighboring model parameters be of similar magnitude. For this reason, the second term on the RHS of (4 is often referred to as a smoothness constraint. The strength
3 30 March 005 of the weight λ 1 specifies how important parameter smoothness should be in relation to the main objective (the first term on the RHS, which is to minimize the residual in the least squares sense. Regularization can be implemented by augmenting the system of equations. In block matrix form, the regularized forward problem is written A d m =, λ 1 0 (5 and m can be found in the usual way by (. Another constraint customarily employed in inversion is norm damping. For unconstrained parameters, it is often desirable to at least constrain them to be physically reasonable. In fact, if we don t have any data that constrains particular parameters, it s just as well to make them zero, or at least something small. Going back to the simple two-parameter problem at the beginning of this section, we might additionally constrain the problem by a second constraint: λ a =0 λ b =0 again, where λ is a weight (agrange multiplier specifying the importance of this objective. Applying this idea to (4, we create a new penalty function, ( T ε= m A A 1 A T d + m +λ m. (6 Squared error λ 1 Smoothness Norm constraint Constraint Damping applies a norm constraint to the problem. That is to say, it constrains the parameters to be small; large values in m would make its norm large and would therefore make ε big. In the limit as λ Æ, all the parameters would be zero. The augmented forward problem is expressed as A d λ1 m = 0, λ I 0 and, as with (5, m can be found by employing (. (7 RESOUTION Equation (1 establishes a forward mapping from model to data space (cartooned by the top arrow in the figure at right and Equation ( establishes the inverse mapping from data to model space Data Space A T T ( A A 1 A Model Space 3
4 30 March 005 (cartooned by the bottom arrow. We would especially like to know how well we resolve model parameters via the inverse mapping. Combining (1 and (, we get m ( A A ( A A m, (8 and we define the resolution matrix as so that R ( T 1 T T 1 T A A ( A m Rm. A (9 Clearly, if R = I then all the parameters are equally and perfectly resolved. In practice, R will often be diagonally dominant, but the main diagonals will not all be one and the off-diagonals will often be different from zero. By analyzing the diagonals of R, we gather information on how well resolved the parameters are (resolution increases as the diagonal values approach unity. The off-diagonal terms tell us about smearing between parameters. However, in seismic tomography, A frequently contains billions or even trillions of elements, so formally computing R is impractical (impossible?. There are a number of ways to circumnavigate this difficulty. The two methods described here are the checkerboard resolution test and the hypothesis test. The checkerboard resolution test basically utilizes a synthetic Earth model wherein seismic velocities alternate between slow and fast, in a checkerboard pattern: m syn is mapped to synthetic data by m syn checkerboard. d syn = Am syn, and then d syn is inverted to create the inversion model We can then consider the resolution metric T 1 T syn m = ( A A A d. m m syn syn or we can simply visually inspect m syn to see how well it tallies with m syn. The following observations can be made about checkerboard resolution tests: The greatest model resolution, not surprisingly, occurs over continental landmasses Many regions over open oceans apparently cannot be resolved at all. Deep mantle model parameters are generally better resolved than ones nearer the surface because there is better data coverage provided by deep diving waves. Model regularizations used in the inversions tend to provide poorer amplitude resolution than spatial resolution. These test types assume A is known. But A is an approximation of physical reality. syn 4
5 30 March 005 The hypothesis testing method is less arbitrary than the checkerboard method. Rather than test the resolution of an unrealistic Earth model, which is only marginally informative at best, this method tests a null hypothesis. The null hypothesis can be any Earth model that is similar to the current inversion model, but that differs in some way that would be important from an interpretive viewpoint. et the current and null hypothesis models be denoted, respectively, m 0 and m null. The null hypothesis Earth model is mapped to data space via A and then inverted back -1 T -1 A so we re considering (A T A (A T A m null. If m null is similar to to model space via (A T A -1 (A T A (A T A m null at those places where m null was purposely made to differ from m 0, then there -1 is sufficient resolution in these regions and the null hypothesis is rejected. If (A T A (A T A m null looks like m 0 then the null hypothesis is accepted the data do not constrain the parameters enough to distinguish between m 0 and m null. So, when the null hypothesis is rejected, we have resolution in regions of interest; but if it is accepted, then both the null hypothesis and the original inversion model can equally explain the data, which presents a problem as far as interpretation is concerned. Hypothesis testing is easily visualized using the block diagram shown in Figure 1. Tomographic model with anomaly ( m 0 Null hypothesis model (m null with dashed outline of m 0 for comparison Hypothesis test Figure 1 Block diagram cartooning hypothesis testing for resolution. The top left panel depicts a tomographic image with one anomaly. We would T 1 T ( like to know, for example, whether the A A ( A A m null data provides enough resolution to unequivocally interpret this anomaly as a contiguous body or perhaps as two separate bodies (as seen in the top right nu ll hypothesis panel. To do this, apply (A T A -1 A T A to m null and check to see whether the two-anomaly or hypothesis is preserved. If it is, the hypothesis is rejected; there is enough resolution to distinguish between the two models. If not, the hypothesis is Reject Resolution is sufficient Accept Resolution is accepted; we do not have enough to distinguish between models. insufficient to distinguish resolution to reasonably distinguish between models. between the null hypothesis and the original image. 5
6 30 March 005 FRESNE ZONES AND MUTIPATH/SCATTERING We have d = Gδµ dv Volume but we only need to worry about those portions of Gδµ that are nonzero. Geometric Ray Theory employs a high-frequency approximation that allows us to assume that Gδµ is only nonzero along the ray path. So the preceding integral can be reduced to δ t = s dl. ray This simplification breaks down when we consider finite frequencies. For finite frequencies, we employ the analytical tools of Fresnel zones and multipathing/scattering. These concepts provide more realistic 3D kernels applicable to Earth models. Fresnel Zones S R Rays traveling from S to R will follow not only the GRT raypath, but also any number of excursive paths. If the lengths of these excursive paths,, do not differ from (the GRT raypath length by more λ /4, they will constructively interfere. In other words: < λ 4 constructive interference T T < T 4 = s + w + ( s + w w S λ 1/ w S - S 8 R w denotes the radius of a Fresnel Zone. It specifies a cross-sectional area through which constructively interfering rays travel from S to R. So the waveform observed at R is no longer dependent only on the GRT raypath, but the volume described by the Fresnel Zone. The arrival at R has now sampled a volume, as opposed to a line. ~ 1000km Æ w ~ 30km ~10,000km Æ w ~ 100km The excursive raypaths arrive later than the GRT raypath, creating constructive interference and the arrival coda. 6
7 30 March 005 The Fresnel Zone is shaped like a banana and is frequently referred to as the banana donut: In cross-section, the banana donut is represented by a series of concentric contours that correspond to the spatial sensitivity of the banana donut at each point. Interestingly, for finite frequencies, the sensitivity is identically zero at the center of the cross-section, which provides the counterintuitive insight that the GRT raypath is a path of zero sensitivity! Banana Donut Sensitivity Kernel Zero sensitivity The banana donut kernel (BDK replaces the elements in A, as described earlier in these notes. Because the BDK is a volume and not a raypath, A is no longer nearly as sparse as it is in the case of GRT. Moreover, because sensitivity is now distributed in a volume, the BDK A is less ill-posed than the GRT A. Hence, banana donut kernels will tend to increase resolution, as compared to GRT kernels. Multipathing The preceding ideas also naturally give rise to multipathing, which acknowledges that more than one raypath can satisfy Fermat s Principle. For example, a scattered ray might travel through a faster region and therefore arrive at the same time as the GRT ray. Consequently, the integral fast d = Volume Gδµ dv can be simplified to the sum of two integrals and this is the central idea of multipathing. δt = s dl ray1 ray + s dl, 7
(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 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 informationModel parametrization strategies for Newton-based acoustic full waveform
Model parametrization strategies for Newton-based acoustic full waveform inversion Amsalu Y. Anagaw, University of Alberta, Edmonton, Canada, aanagaw@ualberta.ca Summary This paper studies the effects
More informationIntroduction to seismic body waves Tomography
Introduction to seismic body waves Tomography Presented by Han Zhang Prepared by group work of Justin Wilgus and Han Zhang Outlines 1. Velocity derived from wave equation 2. Two commonly used reference
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 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 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 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 informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More informationTravel Time Tomography using Neural Networks
Travel Time Tomography using Neural Networks Yoshiya Oda Tokyo Metropolitan University, Japan Tomohisa Ishiyama Tokyo Metropolitan University, Japan Shinya Yokono Tokyo Metropolitan University, Japan SUMMARY:
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 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 informationA simple method for determining the spatial resolution of a
A simple method for determining the spatial resolution of a general inverse problem Meijian An Institute of Geomechanics, CAGS, Beijing 100081, China Email: meijianan@yahoo.com.cn ABSTRACT The resolution
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 informationModelling, migration, and inversion using Linear Algebra
Modelling, mid., and inv. using Linear Algebra Modelling, migration, and inversion using Linear Algebra John C. Bancroft and Rod Blais ABSRAC Modelling migration and inversion can all be accomplished using
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More informationApplications of the tau-p Inversion Method. Shoba Maraj
Applications of the tau-p Inversion Method Shoba Maraj Final Report Theoretical Seismology Geop 523 Due 6th May 2011 Introduction Generally, velocity increases with depth in the Earth. They are related
More informationFermat s principle and ray tracing in anisotropic layered media
Ray tracing in anisotropic layers Fermat s principle and ray tracing in anisotropic layered media Joe Wong ABSTRACT I consider the path followed by a seismic signal travelling through velocity models consisting
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 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 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 informationSource Estimation for Wavefield Reconstruction Inversion
Source Estimation for Wavefield Reconstruction Inversion Zhilong Fang * and Felix J. Herrmann * * Seismic Laboratory for Imaging and Modeling SLIM), University of British Columbia Abstract Wavefield reconstruction
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 informationAmplitude within the Fresnel zone for the zero-offset case
Amplitude within the Fresnel zone for the zero-offset case Shuang Sun and John C. Bancroft ABSTRACT Reflection energy from a linear reflector comes from an integrant over an aperture often described by
More informationAdvanced Operations Research Techniques IE316. Quiz 1 Review. Dr. Ted Ralphs
Advanced Operations Research Techniques IE316 Quiz 1 Review Dr. Ted Ralphs IE316 Quiz 1 Review 1 Reading for The Quiz Material covered in detail in lecture. 1.1, 1.4, 2.1-2.6, 3.1-3.3, 3.5 Background material
More informationShortest-path calculation of first arrival traveltimes by expanding wavefronts
Stanford Exploration Project, Report 82, May 11, 2001, pages 1 144 Shortest-path calculation of first arrival traveltimes by expanding wavefronts Hector Urdaneta and Biondo Biondi 1 ABSTRACT A new approach
More informationFull waveform inversion by deconvolution gradient method
Full waveform inversion by deconvolution gradient method Fuchun Gao*, Paul Williamson, Henri Houllevigue, Total), 2012 Lei Fu Rice University November 14, 2012 Outline Introduction Method Implementation
More informationIterative Algorithms I: Elementary Iterative Methods and the Conjugate Gradient Algorithms
Iterative Algorithms I: Elementary Iterative Methods and the Conjugate Gradient Algorithms By:- Nitin Kamra Indian Institute of Technology, Delhi Advisor:- Prof. Ulrich Reude 1. Introduction to Linear
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 informationA projected Hessian matrix for full waveform inversion Yong Ma and Dave Hale, Center for Wave Phenomena, Colorado School of Mines
A projected Hessian matrix for full waveform inversion Yong Ma and Dave Hale, Center for Wave Phenomena, Colorado School of Mines SUMMARY A Hessian matrix in full waveform inversion (FWI) is difficult
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 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 information1. Techniques for ill-posed least squares problems. We write throughout this lecture = 2. Consider the following least squares problems:
ILL-POSED LEAST SQUARES PROBLEMS. Techniques for ill-posed least squares problems. We write throughout this lecture = 2. Consider the following least squares problems: where A = min b Ax, min b A ε x,
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 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 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 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 informationTh E An Interferometric Interpretation of Marchenko Redatuming
Th E0 An Interferometric Interpretation of Marchenko Redatuming J. van der Neut (Delft University of Technology), I. Vasconcelos (chlumberger Gould Research) & K. Wapenaar (Delft University of Technology)
More informationIntroduction to Geophysical Inversion
Introduction to Geophysical Inversion Goals Understand the non-uniqueness in geophysical interpretations Understand the concepts of inversion. Basic workflow for solving inversion problems. Some important
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 informationContents. I Basics 1. Copyright by SIAM. Unauthorized reproduction of this article is prohibited.
page v Preface xiii I Basics 1 1 Optimization Models 3 1.1 Introduction... 3 1.2 Optimization: An Informal Introduction... 4 1.3 Linear Equations... 7 1.4 Linear Optimization... 10 Exercises... 12 1.5
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 informationLECTURE 13: SOLUTION METHODS FOR CONSTRAINED OPTIMIZATION. 1. Primal approach 2. Penalty and barrier methods 3. Dual approach 4. Primal-dual approach
LECTURE 13: SOLUTION METHODS FOR CONSTRAINED OPTIMIZATION 1. Primal approach 2. Penalty and barrier methods 3. Dual approach 4. Primal-dual approach Basic approaches I. Primal Approach - Feasible Direction
More informationThe interpolation of a 3-D data set by a pair of 2-D filters
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 275 The interpolation of a 3-D data set by a pair of 2-D filters Matthias Schwab and Jon F. Claerbout 1 ABSTRACT Seismic 3-D field data is
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 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 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 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 informationWave-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 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 informationIMAGING USING MULTI-ARRIVALS: GAUSSIAN BEAMS OR MULTI-ARRIVAL KIRCHHOFF?
IMAGING USING MULTI-ARRIVALS: GAUSSIAN BEAMS OR MULTI-ARRIVAL KIRCHHOFF? Summary Samuel H. Gray* Veritas DGC Inc., 715 Fifth Ave. SW, Suite 2200, Calgary, AB Sam_gray@veritasdgc.com Carl Notfors Veritas
More 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 informationConvolution Product. Change of wave shape as a result of passing through a linear filter
Convolution Product Change of wave shape as a result of passing through a linear filter e(t): entry signal (source signal) r(t): impulse response (reflectivity of medium) (a) The spikes are sufficiently
More informationWave-equation migration velocity analysis II: Subsalt imaging examples. Geophysical Prospecting, accepted for publication
Wave-equation migration velocity analysis II: Subsalt imaging examples Geophysical Prospecting, accepted for publication Paul Sava and Biondo Biondi Stanford Exploration Project, Mitchell Bldg., Department
More informationThree critical concerns for marine seismics with portable systems. Source strength and tuning Streamer length 2D vs. 3D
Three critical concerns for marine seismics with portable systems Source strength and tuning Streamer length 2D vs. 3D Airgun Source Single airgun Multiple airguns signal strength volume 1/3 1 x (200 in
More 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 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 informationCopyright 2007 Pearson Addison-Wesley. All rights reserved. A. Levitin Introduction to the Design & Analysis of Algorithms, 2 nd ed., Ch.
Iterative Improvement Algorithm design technique for solving optimization problems Start with a feasible solution Repeat the following step until no improvement can be found: change the current feasible
More informationTemporal Integration Of Seismic Traveltime Tomography
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,
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 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 informationMulti-source Least-squares Migration of Gulf of Mexico Data
Multi-source Least-squares Migration of Gulf of Mexico Data Xin Wang. King Abdullah University of Science and Technology, Thuwal 955-69, Kingdom of Saudi Arabia Corresponding author is Xin Wang. E-mail
More informationShort Note. Non-stationary PEFs and large gaps. William Curry 1 INTRODUCTION
Stanford Exploration Project, Report 120, May 3, 2005, pages 247 257 Short Note Non-stationary PEFs and large gaps William Curry 1 INTRODUCTION Prediction-error filters (PEFs) may be used to interpolate
More informationData Acquisition. Chapter 2
Data Acquisition Chapter 2 1 st step: get data Data Acquisition Usually data gathered by some geophysical device Most surveys are comprised of linear traverses or transects Typically constant data spacing
More informationAmplitude and kinematic corrections of migrated images for non-unitary imaging operators
Stanford Exploration Project, Report 113, July 8, 2003, pages 349 363 Amplitude and kinematic corrections of migrated images for non-unitary imaging operators Antoine Guitton 1 ABSTRACT Obtaining true-amplitude
More 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 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 informationA Probabilistic Approach to the Hough Transform
A Probabilistic Approach to the Hough Transform R S Stephens Computing Devices Eastbourne Ltd, Kings Drive, Eastbourne, East Sussex, BN21 2UE It is shown that there is a strong relationship between the
More informationFourier Transforms. Laboratory & Computational Physics 2
Fourier Transforms Laboratory & Computational Physics 2 Last compiled August 8, 2017 1 Contents 1 Introduction 3 1.1 Prelab questions................................ 3 2 Background theory 5 2.1 Physical
More informationFourier Transforms. Laboratory & Comp. Physics 2
Fourier Transforms Laboratory & Comp. Physics 2 Last compiled August 8, 2017 Contents 1 Introduction 3 1.1 Prelab questions............ 3 2 Background theory 6 2.1 Physical discussion of Fourier transforms.................
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 informationIrradiance Gradients. Media & Occlusions
Irradiance Gradients in the Presence of Media & Occlusions Wojciech Jarosz in collaboration with Matthias Zwicker and Henrik Wann Jensen University of California, San Diego June 23, 2008 Wojciech Jarosz
More informationContrast Optimization A new way to optimize performance Kenneth Moore, Technical Fellow
Contrast Optimization A new way to optimize performance Kenneth Moore, Technical Fellow What is Contrast Optimization? Contrast Optimization (CO) is a new technique for improving performance of imaging
More informationLessons in Spatial Sampling Or... Does Anybody Know the Shape of a Wavefield?
Lessons in Spatial Sampling Or... Does Anybody Know the Shape of a Wavefield? Norm Cooper - Mustagh Resources Ltd., Calgary, Canada In a homogeneous, isotropic, infinite half space (remember those regimes
More informationSummary. Figure 1: Simplified CRS-based imaging workflow. This paper deals with the boxes highlighted in green.
Smoothing and automated picking of kinematic wavefield attributes Tilman Klüver and Jürgen Mann, Geophysical Institute, University of Karlsruhe, Germany Copyright 2005, SBGf Sociedade Brasiliera de Geofísica
More informationIntroduction to Geometric Algebra
Introduction to Geometric Algebra Lecture 6 Intersection and Union of Subspaces Professor Leandro Augusto Frata Fernandes laffernandes@ic.uff.br Lecture notes available in http://www.ic.uff.br/~laffernandes/teaching/2013.1/topicos_ag
More informationVidale, J., 1988, Finite-difference calculation of travel times, BSSA, Vo. 78, No. 6,
Class 4: First-Arrival Traveltime Tomography Mon, Sept 21, 2009 Wavefront tracing methods Inversion algorithms Model regularization Review of tomography case histories, presentation by Dr. Xianhuai Zhu
More informationHeadwave Stacking in Terms of Partial Derivative Wavefield
Geosystem Engineering, 7(1), 21-26 (March 2004) Headwave Stacking in Terms of Partial Derivative Wavefield Changsoo Shin School of Civil, Urban and Geosystem Engineering, Seoul National University, San
More informationREGULAR GRAPHS OF GIVEN GIRTH. Contents
REGULAR GRAPHS OF GIVEN GIRTH BROOKE ULLERY Contents 1. Introduction This paper gives an introduction to the area of graph theory dealing with properties of regular graphs of given girth. A large portion
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 informationLecture 9 - Matrix Multiplication Equivalences and Spectral Graph Theory 1
CME 305: Discrete Mathematics and Algorithms Instructor: Professor Aaron Sidford (sidford@stanfordedu) February 6, 2018 Lecture 9 - Matrix Multiplication Equivalences and Spectral Graph Theory 1 In the
More informationMath 308 Autumn 2016 MIDTERM /18/2016
Name: Math 38 Autumn 26 MIDTERM - 2 /8/26 Instructions: The exam is 9 pages long, including this title page. The number of points each problem is worth is listed after the problem number. The exam totals
More informationWave-equation migration velocity analysis. II. Subsalt imaging examples
Geophysical Prospecting, 2004, 52, 607 623 Wave-equation migration velocity analysis. II. Subsalt imaging examples P. Sava and B. Biondi Department of Geophysics, Stanford University, Mitchell Building,
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 informationSimulation of Multipoint Ultrasonic Flowmeter
Simulation of Multipoint Ultrasonic Flowmeter Jakub Filipský 1,*, Jiří Nožička 2 1,2 CTU in Prague, Faculty of Mechanical Engineering, Department of Fluid Mechanics and Thermodynamics, Technická 4, 166
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 informationUNIVERSITY OF CALGARY. Traveltime Tomography in Isotropic and Transversely Isotropic Media. Marco Perez
Important Notice This copy may be used only for the purposes of research and private study, and any use of the copy for a purpose other than research or private study may require the authorization of the
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 informationImplementing non-stationary filtering in time and in Fourier domain
Stanford Exploration Project, Report 111, June 9, 2002, pages 347 361 Implementing non-stationary filtering in time and in Fourier domain Gabriel Alvarez 1 ABSTRACT Non-stationary filtering of seismic
More informationThe power of stacking, Fresnel zones, and prestack migration
The power of stacking, Fresnel zones, and prestack migration John C. Bancroft and Shuang Sun ABSTRACT The stacking of common midpoint (CMP) gathers assumes the presence of specula reflection energy and
More informationIn other words, we want to find the domain points that yield the maximum or minimum values (extrema) of the function.
1 The Lagrange multipliers is a mathematical method for performing constrained optimization of differentiable functions. Recall unconstrained optimization of differentiable functions, in which we want
More informationThe Interface Fresnel Zone revisited
The Interface Fresnel Zone revisited Nathalie Favretto-Cristini, Paul Cristini and Eric de Bazelaire ABSTRACT We determine the part of reflectors which actually affects the reflected wavefield, which is
More informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Full Waveform Inversion Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 18, 2014 Bee Bednar (Panorama Technologies)
More informationScanning Real World Objects without Worries 3D Reconstruction
Scanning Real World Objects without Worries 3D Reconstruction 1. Overview Feng Li 308262 Kuan Tian 308263 This document is written for the 3D reconstruction part in the course Scanning real world objects
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 information25 The vibration spiral
25 The vibration spiral Contents 25.1 The vibration spiral 25.1.1 Zone Plates............................... 10 25.1.2 Circular obstacle and Poisson spot.................. 13 Keywords: Fresnel Diffraction,
More informationDOWNLOAD PDF BIG IDEAS MATH VERTICAL SHRINK OF A PARABOLA
Chapter 1 : BioMath: Transformation of Graphs Use the results in part (a) to identify the vertex of the parabola. c. Find a vertical line on your graph paper so that when you fold the paper, the left portion
More informationChapter 15 Introduction to Linear Programming
Chapter 15 Introduction to Linear Programming An Introduction to Optimization Spring, 2015 Wei-Ta Chu 1 Brief History of Linear Programming The goal of linear programming is to determine the values of
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 informationProcessing converted-wave data in the tau-p domain: rotation toward the source and moveout correction
τ-p domain converted-wave processing Processing converted-wave data in the tau-p domain: rotation toward the source and moveout correction Raul Cova and Kris Innanen ABSTRACT The asymmetry of the converted-wave
More information