Programming for Geophysics
|
|
- Howard Wilkins
- 5 years ago
- Views:
Transcription
1 Programming for Geophysics Bill Harlan May 21, 2008
2
3 Some ways to program as a geophysicist Avoid it. Borrow code. Find partners. Use Matlab.
4 Some ways to program as a geophysicist Avoid it. Borrow code. Find partners. Use Matlab. Become a guru/sysadmin. Help others publish first.
5 Some ways to program as a geophysicist Avoid it. Borrow code. Find partners. Use Matlab. Become a guru/sysadmin. Help others publish first. Concentrate on numerics. Fortran. One program per dataset.
6 Some ways to program as a geophysicist Avoid it. Borrow code. Find partners. Use Matlab. Become a guru/sysadmin. Help others publish first. Concentrate on numerics. Fortran. One program per dataset. Build a personal library. Generalize your code for reuse.
7 Learn fundamentals deliberately, not as you go Take a course, read books Data structures, algorithms, object-oriented, functional
8 Learn fundamentals deliberately, not as you go Take a course, read books Data structures, algorithms, object-oriented, functional Learn as a branch of math, not engineering. Concentrate on reusable abstractions, not popular toolkits. Master simplicity, not complexity.
9 Learn fundamentals deliberately, not as you go Take a course, read books Data structures, algorithms, object-oriented, functional Learn as a branch of math, not engineering. Concentrate on reusable abstractions, not popular toolkits. Master simplicity, not complexity. Do not get carried away.
10 Learn best software practices Show and share Source control Tests Small changes (refactoring) Appropriate generalization
11 Examples of generalization/abstraction Seismic data objects with flexible dimensions Separate velocity models from ray tracers Different imaging conditions with different extrapolators
12 Typical geophysical inversions Data simulated by series of non-linear operations Inversion is both over- and under-determined No model parameters fit data perfectly Many models fit data equally well Non-linearity is well-behaved
13 Sensitivity of interval velocity to RMS errors
14 Dix inversion Forward equation cannot fit arbitrary data: j Vj rms 1 = j (Vk int)2 k=1 Explicit inverse may not be physical: Vj int = j(v rms j ) 2 (j 1)(V rms j 1 )2 Instead minimize damped least-squares: j (Vrms j ) 2 [ 1 j j ] 1 (Vk int )2 k=1 2 +ǫ k (V int k ) 2
15 Damp interval velocity roughness
16 Defining an inversion Do not define your solution by the way you solve it. Want to improve the solution without redefining the problem. Instead, identify an objective function (or probabilities). E.g., define rays by minimum time.
17 Lomask s flattening, redone Estimate vertical stretch that flattens reflections. Original: Custom regression, phase-unwrapping New version: A few hundred lines of code specific to inversion JTK reused: structure tensors, Gaussian filters, Gauss-Newton
18 Local dips from structure tensors
19 Estimated vertical shifts in color
20 Flattened with vertical shifts
21 The problem, the data, and the solution Flatten seismic structure with vertical shift τ(x,y,t): flat(x,y,t) = structure[x,y,t +τ(x,y,t)]. Data are slopes p x, p y measured from structure tensors. Want x τ(x,y,t) p x(x,y,t) and y τ(x,y,t) p y(x,y,t). min τ(x,y,t) ( x τ(x,y,t) p x(x,y,t) 2 + y τ(x,y,t) p y(x,y,t) 2 +ǫ τ(x,y,t) 2 )dx dy dt
22 Looks like damped least-squares The best model m fits the data d with a function f(d) by minimizing the vector norms d f(m) 2 d + m 2 m or [d f(m)] 1 d [d f(m)]+m C 1 m m. C Optional covariances: C d E(dd ) and C m E(mm ).
23 Gauss-Newton inversion Finds m to minimize [d f(m)] 1 d C [d f(m)] + [m m 0] 1 m [m m 0 ] C Algorithm: 1. Accepts starting model m 0 2. Approximates f(m 0 + m) f(m 0 )+F m 3. Conjugate-gradients minimizes quadratic for m 4. Line search scales perturbation: m 0 +α m 5. Adds perturbation to reference model for new m 0 6. Returns to step 2
24 Required operations Simulate data from model: d = f(m) Perturb data with model perturbation: d = F (m 0 ) m f(m 0 + m) f(m 0 ) Perturb model with transpose: F (m 0 ) d
25 What is that transpose? Use definition: d (F m) (F d) m Discrete: swap summations Continuous: integrate by parts Examples: Smoothing Smoothing Convolution Correlation Derivative Negative derivative Plane-wave modeling Slant stacks Seismic modeling Migration
26 Inversion sees three abstract operations Vector data = forward(vector model) Vector data = linearized(vector model, Vector refmodel) Vector model = transpose(vector data, Vector refmodel)
27 Required operations for both data and model Vector { scale(float scalar) add(vector other) dot(vector other) [required] } multiplyinversecovariance() [optional] applyhardconstraint()
28 Constrained Dix inversion Solve for smooth interval slowness m. Minimize errors in squared stacking slowness d Forward transform: j 1/d j = (1/j) (1/mk 2 ) k=1 Linearized transform: j d j = (2dj 2 /j) (1/mk 3 ) m k k=1 Transpose transform: m k = (1/mk 3 ) (2dj 2 /j) d j j=k
29 Other applications Tomography: reflection, cross-well, diving, amplitude Generalized Radon transforms Surface-consistent deconvolution Normal moveout corrections Automatic moveout picking Coherency, wavelet/phase attributes Tests for simulations
30 Conclusions More time on computer science quickly saves time Look for opportunities to generalize
31
32 Alternative to covariance Insert simplification filter d = m) where f(s S C 1 m Ĩ S If C 1 m m roughens, then m smooths. S Faster than covariance constraint Can change dynamically during optimization
Stanford Exploration Project, Report 120, May 3, 2005, pages
Stanford Exploration Project, Report 120, May 3, 2005, pages 167 179 166 Stanford Exploration Project, Report 120, May 3, 2005, pages 167 179 Non-linear estimation of vertical delays with a quasi-newton
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 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 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 informationSolution steering with space-variant filters
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 225?? Solution steering with space-variant filters Robert G. Clapp, Sergey Fomel, and Jon Claerbout 1 ABSTRACT Most geophysical problem
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 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 informationAccelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method
Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Wenyong Pan 1, Kris Innanen 1 and Wenyuan Liao 2 1. CREWES Project, Department of Geoscience,
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 informationAngle Gathers for Gaussian Beam Depth Migration
Angle Gathers for Gaussian Beam Depth Migration Samuel Gray* Veritas DGC Inc, Calgary, Alberta, Canada Sam Gray@veritasdgc.com Abstract Summary Migrated common-image-gathers (CIG s) are of central importance
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 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 informationTarget-oriented wave-equation inversion
Stanford Exploration Project, Report 120, May 3, 2005, pages 23 40 Target-oriented wave-equation inversion Alejandro A. Valenciano, Biondo Biondi, and Antoine Guitton 1 ABSTRACT A target-oriented strategy
More informationExtending the search space of time-domain adjoint-state FWI with randomized implicit time shifts
Extending the search space of time-domain adjoint-state FWI with randomized implicit time shifts Mathias Louboutin * and Felix J. Herrmann Seismic Laboratory for Imaging and Modeling (SLIM), The University
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 informationWell-log derived step lengths in Full Waveform Inversion Tunde Arenrin Gary Margrave
Well-log derived step lengths in Full Waveform Inversion Tunde Arenrin Gary Margrave Motivation»Being able to make FWI a part of the seismic processing work flow under current computational capability.»incorporating
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 informationRay based tomography using residual Stolt migration
Stanford Exploration Project, Report 11, November 11, 00, pages 1 15 Ray based tomography using residual Stolt migration Robert G. Clapp 1 ABSTRACT In complex areas, residual vertical movement is not an
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 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 informationStanford Exploration Project, Report 103, April 27, 2000, pages
Stanford Exploration Project, Report 103, April 27, 2000, pages 205 324 204 Stanford Exploration Project, Report 103, April 27, 2000, pages 205 324 Data alignment with non-stationary shaping filters James
More informationZ-99 3D Sub-salt Tomography Based on Wave Equation Migration Perturbation Scans
1 Z-99 3D Sub-salt Tomography Based on Wave Equation Migration Perturbation Scans BIN WANG 1, VOLKER DIRKS 1, PATRICE GUILLAUME 2, FRANÇOIS AUDEBERT 1, ANNING HOU 1 AND DURYODHAN EPILI 1 1 CGG Americas;
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 informationSeparation of specular reflection and diffraction images in Kirchhoff depth migration Faruq E Akbar and Jun Ma, SEIMAX Technologies, LP
Separation of specular reflection and diffraction images in Kirchhoff depth migration Faruq E Akbar and Jun Ma, SEIMAX Technologies, LP Summary Seismic diffractions may occur from faults, fractures, rough
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 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 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 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 informationDealing with errors in automatic velocity analysis
Stanford Exploration Project, Report 112, November 11, 2002, pages 37 47 Dealing with errors in automatic velocity analysis Robert G. Clapp 1 ABSTRACT The lack of human interaction in automatic reflection
More informationINTRODUCTION. Model: Deconvolve a 2-D field of random numbers with a simple dip filter, leading to a plane-wave model.
Stanford Exploration Project, Report 105, September 5, 2000, pages 109 123 Short Note Test case for PEF estimation with sparse data II Morgan Brown, Jon Claerbout, and Sergey Fomel 1 INTRODUCTION The two-stage
More informationAdvanced Math Quadratics Review Name: Dec. 2016
Advanced Math Quadratics Review Name: Dec. 2016 Graph the given quadratic by finding the vertex and building a table around it. Identify the axis of symmetry, maximum or minimum value, domain and range
More informationShort Note. DMO velocity analysis with Jacubowicz s dip-decomposition method. David Kessler and Wai-Kin Chan*
GEOPHYSICS, VOL. 58, NO. 10 (OCTOBER 1993); P. 1517-1524,9 FIGS. Short Note DMO velocity analysis with Jacubowicz s dip-decomposition method David Kessler and Wai-Kin Chan* INTRODUCTION Dip-moveout (DMO)
More informationShaping regularization in geophysical-estimation problems
GEOPHYSICS, VOL. 7, NO. MARCH-APRIL 7 ; P. R9 R6, FIGS., TABLE..9/.476 Shaping regularization in geophysical-estimation problems Sergey Fomel ABSTRACT Regularization is a required component of geophysical-estimation
More informationA MATLAB version of Equivalent Offset Migration (EOM)
A MATLAB version of Equivalent Offset Migration (EOM) John C. Bancroft ABSTRACT Equivalent offset migration continues to be a viable method of prestack migration that is used in anisotropic prestack depth
More informationMigration from a non-flat datum via reverse-time extrapolation
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 50 Migration from a non-flat datum via reverse-time extrapolation Gopal Palacharla 1 ABSTRACT Land surveys usually have elevation changes,
More 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 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 informationIt is widely considered that, in regions with significant
Multifocusing-based multiple attenuation Alex Berkovitch and Kostya Deev, Geomage Evgeny Landa, OPERA It is widely considered that, in regions with significant geologic complexity, methods which work directly
More informationConjugate guided gradient (CGG) method for robust inversion and its application to velocity-stack inversion a
Conjugate guided gradient (CGG) method for robust inversion and its application to velocity-stack inversion a a Published in Geophysics, 71, no. 4, R59-R67, (2006) Jun Ji ABSTRACT This paper proposes a
More informationMultichannel deconvolution imaging condition for shot-profile migration
Stanford Exploration Project, Report 113, July 8, 2003, pages 127 139 Multichannel deconvolution imaging condition for shot-profile migration Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT A significant
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 informationPLEASE NOTE THAT THERE ARE QUESTIONS IN SOME OF THESE SECTIONS THAT ARE TO BE TURNED IN AS HOMEWORK.
Seismic Reflection Lab for Geology 594 BEFORE YOU DO ANYTHING INSERT THE DONGLE ON THE LEFT SIDE OF THE KEYBOARD (LIGHT SIDE UP). The little green light will come on showing the dongle is activated. You
More informationPrestack residual migration in the frequency domain
GEOPHYSICS, VOL. 68, NO. (MARCH APRIL 3); P. 634 64, 8 FIGS. 1.119/1.156733 Prestack residual migration in the frequency domain Paul C. Sava ABSTRACT Prestack Stolt residual migration can be applied to
More informationVelocity versus Offset (VVO) Estimation Using Local Event Correlation and Its Application in Seismic Processing and Analysis
J. Eng. Technol. Sci., Vol. 45, No., 013, 15-139 15 Velocity versus Offset (VVO) Estimation Using Local Event Correlation and Its Application in Seismic Processing and Analysis Supriyono, Awali Priyono,
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 informationRay-based tomography with limited picking
Stanford Exploration Project, Report 110, September 18, 2001, pages 1 113 Short Note Ray-based tomography with limited picking Robert G. Clapp 1 INTRODUCTION In ray-based reflection tomography picking
More informationAVO, migration apertures Fresnel zones, stacking, Q, anisotropy, and fast inversions
AVO, migration apertures Fresnel zones, stacking, Q, anisotropy, and fast inversions John C. Bancroft and grad students University of Calgary CREWES 2003 NSERC 1 Grad students Shuang Sun Pavan Elapavuluri
More informationMissing trace interpolation and its enhancement of seismic processes
Missing trace interpolation Missing trace interpolation and its enhancement of seismic processes Wai-kin Chan and Robert R. Stewart ABSTRACT Many multi-channel seismic algorithms assume that the input
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 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 informationCoherent noise attenuation using Inverse Problems and Prediction Error Filters
Stanford Exploration Project, Report 105, September 5, 2000, pages 1 49 Coherent noise attenuation using Inverse Problems and Prediction Error Filters Antoine Guitton 1 ABSTRACT Two iterative methods that
More information3. parallel: (b) and (c); perpendicular (a) and (b), (a) and (c)
SECTION 1.1 1. Plot the points (0, 4), ( 2, 3), (1.5, 1), and ( 3, 0.5) in the Cartesian plane. 2. Simplify the expression 13 7 2. 3. Use the 3 lines whose equations are given. Which are parallel? Which
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 informationEARTH SCIENCES RESEARCH JOURNAL
EARTH SCIENCES RESEARCH JOURNAL Earth Sci. Res. J. Vol. 10, No. 2 (December 2006): 117-129 ATTENUATION OF DIFFRACTED MULTIPLES WITH AN APEX-SHIFTED TANGENT- SQUARED RADON TRANSFORM IN IMAGE SPACE Gabriel
More informationReview of seismic imaging: Prestack
Review of seismic imaging: Prestack John C. Bancroft Published in the CSEG Recorder - November 001 Review of seismic migration: Prestack ABSTRACT The kinematics of prestack data considers an arbitrary
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 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 informationAngle-domain parameters computed via weighted slant-stack
Angle-domain parameters computed via weighted slant-stack Claudio Guerra 1 INTRODUCTION Angle-domain common image gathers (ADCIGs), created from downward-continuation or reverse time migration, can provide
More informationImprovements in time domain FWI and its applications Kwangjin Yoon*, Sang Suh, James Cai and Bin Wang, TGS
Downloaded 0/7/13 to 05.196.179.38. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ Improvements in time domain FWI and its applications Kwangjin Yoon*,
More informationTarget-oriented wavefield tomography: A field data example
Target-oriented wavefield tomography: A field data example Yaxun Tang and Biondo Biondi ABSTRACT We present a strategy for efficient migration velocity analysis in complex geological settings. The proposed
More informationAn Integrated Pre-Stack Depth Migration Workflow Using Model- Based Velocity Estimation And Refinement
An Integrated Pre-Stack Depth Migration Workflow Using Model- Based Velocity Estimation And Refinement Abstract: Andy Furniss * The advent of velocity-sensitive depth imaging techniques such as pre-stack
More informationA Fortran90 implementation of symmetric nonstationary phaseshift extrapolator
A Fortran90 implementation of symmetric nonstationary phaseshift extrapolator Yanpeng Mi and Gary F Margrave ABSTRACT Symmetric nonstationary phase-shift (SNPS) for 2D (SNPS2D) is a 2D prestack shot-gather
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 informationStanford Exploration Project, Report 120, May 3, 2005, pages
Stanford Exploration Project, Report 120, May 3, 2005, pages 137 159 136 Stanford Exploration Project, Report 120, May 3, 2005, pages 137 159 Update on flattening without picking Jesse Lomask, Antoine
More informationLeast squares Kirchhoff depth migration: important details
Least squares Kirchhoff depth migration: important details Daniel Trad CREWES-University of Calgary Summary Least squares migration has been an important research topic in the academia for about two decades,
More informationEquivalence of source-receiver migration and shot-profile migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 109 117 Short Note Equivalence of source-receiver migration and shot-profile migration Biondo Biondi 1 INTRODUCTION At first glance, shot
More informationWe N Converted-phase Seismic Imaging - Amplitudebalancing Source-independent Imaging Conditions
We N106 02 Converted-phase Seismic Imaging - Amplitudebalancing -independent Imaging Conditions A.H. Shabelansky* (Massachusetts Institute of Technology), A.E. Malcolm (Memorial University of Newfoundland)
More 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 informationWhen is anti-aliasing needed in Kirchhoff migration? a
When is anti-aliasing needed in Kirchhoff migration? a a Published in SEP report, 80, 467-476 (1994) Dimitri Bevc and David E. Lumley ABSTRACT We present criteria to determine when numerical integration
More informationMultifocusing. The Breakthrough in Seismic Imaging
Multifocusing The Breakthrough in Seismic Imaging Conventional Processing GeoMage MF Processing Conventional Processing GeoMage MF Processing Time versus Depth Recorded seismic data Wavefield parameter
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 informationCoherent partial stacking by offset continuation of 2-D prestack data
Stanford Exploration Project, Report 82, May 11, 2001, pages 1 124 Coherent partial stacking by offset continuation of 2-D prestack data Nizar Chemingui and Biondo Biondi 1 ABSTRACT Previously, we introduced
More informationA Python solver for out-of-core, fault tolerant inversion
Stanford Exploration Project, Report 117, October 23, 2004, pages 183?? Short Note A Python solver for out-of-core, fault tolerant inversion Robert G. Clapp 1 INTRODUCTION In the last ten years SEP has
More informationStratigraphic coordinates, a coordinate system tailored to seismic interpretation a
Stratigraphic coordinates, a coordinate system tailored to seismic interpretation a a Published in Geophysical Prospecting, v. 63, 1246-1255, (2015) Parvaneh Karimi and Sergey Fomel ABSTRACT In certain
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 informationReconciling processing and inversion: Multiple attenuation prior to wave-equation inversion
Reconciling processing and inversion: Multiple attenuation prior to wave-equation inversion Claudio Guerra and Alejandro Valenciano ABSTRACT Seismic inversion is very sensitive to the presence of noise.
More informationksa 400 Growth Rate Analysis Routines
k-space Associates, Inc., 2182 Bishop Circle East, Dexter, MI 48130 USA ksa 400 Growth Rate Analysis Routines Table of Contents ksa 400 Growth Rate Analysis Routines... 2 1. Introduction... 2 1.1. Scan
More informationVolumetric flattening: an interpretation tool
Volumetric flattening: an interpretation tool Jesse Lomask and Antoine Guitton *Stanford Exploration Project, Mitchell Bldg., Department of Geophysics, Stanford University, Stanford, CA 94305-2215 3DGeo
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 informationKinematics of Reverse Time S-G Migration
Kinematics of Reverse Time S-G Migration William W. Symes TRIP Seminar Rice University September 2003 www.trip.caam.rice.edu 1 Agenda: explore prestack focussing properties of RT S-G migration, proper
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 informationSUMMARY INTRODUCTION THEORY. Objective function of ERWI
Extended Reflection Waveform Inversion via Differential Semblance Optimization Yujin Liu, William W. Symes and Zhenchun Li, China University of Petroleum (Huadong), Rice University SUMMARY Reflection-based
More informationAttenuation of diffracted multiples with an apex-shifted tangent-squared radon transform in image space
Attenuation of diffracted multiples with an apex-shifted tangent-squared radon transform in image space Gabriel Alvarez, Biondo Biondi, and Antoine Guitton 1 ABSTRACT We propose to attenuate diffracted
More informationGenetic Algorithm for Seismic Velocity Picking
Proceedings of International Joint Conference on Neural Networks, Dallas, Texas, USA, August 4-9, 2013 Genetic Algorithm for Seismic Velocity Picking Kou-Yuan Huang, Kai-Ju Chen, and Jia-Rong Yang Abstract
More informationWS 09 Velocity Model Building
Full wave vs Kinematic approaches Goal in mid 80 s: Full Waveform Inversion (FWI) ~ Tomography plus Migration. There seems to exist considerable confusion regarding the sort of information about the subsurface
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 informationSUMMARY METHOD. d(t 2 = τ 2 + x2
Yujin Liu, Zhi Peng, William W. Symes and Wotao Yin, China University of Petroleum (Huadong), Rice University SUMMARY The Radon transform suffers from the typical problems of loss of resolution and aliasing
More informationSummary. Introduction
Yaofeng He and Ru-Shan Wu Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, C 9564 Summary We propose an approach to perform migration and imaging using secondary scattered
More informationUNIVERSITY OF CALGARY. Analysis and application of the Radon transform. Zhihong Cao A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES
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 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 informationLeast squares Kirchhoff depth migration with. preconditioning
Least squares Kirchhoff depth migration with preconditioning Aaron Stanton University of Alberta, Department of Physics, 4-83 CCIS, Edmonton AB T6G E (April, 3) Running head: Least Squares Migration ABSTRACT
More informationMultiple Suppression and Data Quality Enhancement Using Radon Transform: A Case Study
P - 268 Multiple Suppression and Data Quality Enhancement Using Radon Transform: A Case Study Shashank Shekhar Singh, Uma Shankar, Kalachand Sain, NGRI, Hyderabad Summary The ultimate goal of seismic exploration
More informationCLASSIFICATION OF MULTIPLES
Introduction Subsurface images provided by the seismic reflection method are the single most important tool used in oil and gas exploration. Almost exclusively, our conceptual model of the seismic reflection
More informationSeismic-to-well ties by smooth dynamic time warping
Seismic-to-well ties by smooth dynamic time warping Tianci Cui and Gary F. Margrave ABSTRACT Well tying by SDTW Without knowledge of Q or a check-shot/vsp survey, synthetic seismogram has to be manually
More informationMultiple suppression with land data: A case study
Multiple suppression with land data: A case study Gabriel Alvarez 1 ABSTRACT Some of the most important methods for multiple suppression are based on the moveout difference between the hyperbolas corresponding
More informationTarget-oriented wavefield tomography using demigrated Born data
Target-oriented wavefield tomography using demigrated Born data Yaxun Tang and Biondo Biondi ABSTRACT We present a method to reduce the computational cost of image-domain wavefield tomography. Instead
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 informationWave-equation migration velocity analysis by focusing diffractions and reflections
GEOPHYSICS, VOL. 70, NO. 3 (MAY-JUNE 2005); P. U19 U27, 14 FIGS. 10.1190/1.1925749 Wave-equation migration velocity analysis by focusing diffractions and reflections Paul C. Sava 1, Biondo Biondi 2, and
More informationSUMMARY NOISE REDUCTION BY SPATIAL SLOPE ESTIMATION
: steerable and non-linear filter approach Abdulaziz M. Almuhaidib and M. Nafi Toksöz, Earth Resources Laboratory, MIT SUMMARY We present an approach based on local-slope estimation for the separation
More information