Lab 5: 2D FFT (Fortran version)
|
|
- Basil Harrington
- 6 years ago
- Views:
Transcription
1 Due Date: 17:00, Tuesday, November 1, 2011 TA: Xukai Shen Lab 5: 2D FFT (Fortran version) Conrad Schlumberger 1 ABSTRACT In this lab you will modify programs to filter data based on the dip in the frequency domain, and use it to process a VSP and remove multiples from a cmp gather. This lab is, once again, available online. Download the tar file at labs:lab5.tar.gz Then untar it with tar -xvzf Lab5.tar.gz VSP wavefield separation Figure 1 shows a common shot gather from a VSP (vertical seismic profile) shot by Schlumberger somewhere in Italy. A VSP is a seismic survey where the shots are near the surface and the geophones are down a well. Processing VSP s is done differently than conventional surface seismic. This shot gather was recorded by three component receivers, hence the three panels in the Figure. Each represents a different direction of Earth motion. A simple way to process VSP s is to separate the wavefield into upgoing and downgoing waves, apply a linear moveout correction, and stack to give an estimate of the reflectivity profile. For this part of the exercise, you just need to separate the upgoing and downgoing wavefields. You need to do two things, first you need to open fourier.f90 and implement the function fth, input and output definitions are given, along with some intermediate variables, dont change the input or the output, add variables if necessary, make sure you implement the correct forward and adjoint, then you need to modify Dip filter.f90 to apply the 2D FFT correctly and to apply an appropriate mute, 1 conrad@schlum.com
2 BEI - Lab 5 2 2D FFT Figure 1: VSP shot gather preferably with a taper.. You will need to put the parameters you use to describe dip in the parameter files up.p and down.p. Then compile the code by typing gmake Dip filter.x. You can then run gmake upgoing.view and gmake downgoing.view to create and view your results. Once you like your results, find the plot commands in paper.tex and un-comment them to include your results in the paper. You can view the results in the F k domain with the command gmake up fk domain.view and gmake down fk domain.view. Include these in the paper by un-commenting the appropriate plot commands. F k multiple suppression You have just written a program to separate upgoing and downgoing wavefields. You can now use this program to suppress multiples in the CMP gather displayed in Figure 2. Create the file intervel.h containing an RMS velocity field somewhere between the primary and multiple moveout velocities by typing gmake intervel.h. This velocity function will be used to apply normal moveout. The command gmake nomultiples.h will moveout the data cmp.hh with the SEPlib program NMO so that the multiples will now be dipping downwards and the primaries upwards, and will then apply your wave-field separation (dip-filter) program to remove the multiple energy. Put your dip limitations in the parameter file multiple.p. Be conservative - it is better to keep a bit of multiple energy rather than throw out any signal. Finally, we apply inverse NMO (NMO inv=1) to transform the
3 BEI - Lab 5 3 2D FFT data back to a normal cmp gather. You can create and view your results with gmake nomultiples.view. To see the results of your filtering in the F k domain, you can type gmake cmp fk domain.view. Include this result by un-commenting the plot commands in paper.tex. Next, run gmake subsampled.view to subsample the data along the x-axis, and try to remove the multiples using the parameter file sub.p. Having fewer spatial offsets leads to data-aliasing and means the primaries and multiples are poorly separated in the Fourier domain. You can view the Fourier domain result with gmake sub fk domain.view. Include these results in paper.tex also. Figure 2: Multiple contaminated cmp. VSP QUESTIONS 1. Which VSP panel do you think is measuring the vertical component of Earth motion? Why? 2. VSP s are usually higher frequency than surface seismic. Why do you think this is the case?
4 BEI - Lab 5 4 2D FFT 3. When you look at the VSP result in the frequency domain, you see that there is a blank strip in the middle. What caused this? Multiple suppression 1. Seismic amplitudes can often be inverted to estimate rock properties (impedance, etc). Multiples can corrupt the amplitudes though, making this kind of analysis impossible. Do you think F k multiple suppression makes the amplitudes more or less reliable? 2. How did you choose the mute range for the multiples? How did it change for the subsampled case? Extra Credit: Who was Conrad Schlumberger, what was he/she famous for? Who was his/her brother? Extra Credit: With a green pen mark all spelling and gramar error s in this Lab. Also comment on confusing statements in the Lab or other broken and outdated features/questions.! Fourier transform and apply dip f i l t e r program D i p f i l t e r use f o u r i e r i m p l i c i t none complex, a l l o c a t a b l e, dimension ( :, : ) : : data i n t e g e r n1, n2 i n t e g e r hetch, putch c a l l i n i t p a r ( ) c a l l doc ( /homes/ sep / c l a u d i o /sepwww/lab5/lab5/ D i p f i l t e r. r s ) i f (0>= hetch ( n1, d, n1 ) ) then c a l l e r e x i t ( Trouble reading n1 from h i s t o r y ) i f (0>= hetch ( n2, d, n2 ) ) then c a l l e r e x i t ( Trouble reading n2 from h i s t o r y ) a l l o c a t e ( data ( n1, n2 ) )
5 BEI - Lab 5 5 2D FFT c a l l d i p f i l t ( n1, n2, data ) end program D i p f i l t e r s ubroutine d i p f i l t ( n1, n2, data ) i n t e g e r nw, nk, n1, n2, iw, ik, i3, n3 r e a l dw, dk, d1, d2, w, k, w0, k0, vel, wt complex data ( n1, n2 )! Un comment and put your dip parameters here! from param : SOME THAT YOU DECIDE ON i n t e g e r hetch, putch, auxputch i f (0>= hetch ( n1, d,nw ) ) then c a l l e r e x i t ( Trouble reading n1 from h i s t o r y ) i f (0>= hetch ( n2, d, nk ) ) then c a l l e r e x i t ( Trouble reading n2 from h i s t o r y ) i f (0>= hetch ( n3, d, n3 ) ) then c a l l e r e x i t ( Trouble reading n3 from h i s t o r y ) i f (0>= hetch ( d1, f, d1 ) ) then c a l l e r e x i t ( Trouble reading d1 from h i s t o r y ) i f (0>= hetch ( d2, f, d2 ) ) then c a l l e r e x i t ( Trouble reading d2 from h i s t o r y ) i f ( 0. ne. auxputch ( n1, d, n1, fk domain ) ) then c a l l e r e x i t ( Trouble w r i t i n g n1 to aux fk domain ) i f ( 0. ne. auxputch ( n2, d, n2, fk domain ) ) then c a l l e r e x i t ( Trouble w r i t i n g n2 to aux fk domain ) i f ( 0. ne. auxputch ( e s i z e, d,8, fk domain ) ) then c a l l e r e x i t ( Trouble w r i t i n g e s i z e to aux fk domain ) dw = 1/((nw 1) d1 ) dk = 1/(( nk 1) d2 ) w0 = 1/(2 d1 ) k0 = 1/(2 d2 ) i f ( 0. ne. auxputch ( d1, f,dw, fk domain ) ) then c a l l e r e x i t ( Trouble w r i t i n g d1 to aux fk domain ) i f ( 0. ne. auxputch ( d2, f, dk, fk domain ) ) then c a l l e r e x i t ( Trouble w r i t i n g d2 to aux fk domain ) i f ( 0. ne. auxputch ( o1, f,w0, fk domain ) ) then c a l l e r e x i t ( Trouble w r i t i n g o1 to aux fk domain ) i f ( 0. ne. auxputch ( o2, f, k0, fk domain ) ) then c a l l e r e x i t ( Trouble w r i t i n g o2 to aux fk domain ) do i 3 =1,n3 c a l l s r e e d ( in, data, n1 n2 8)
6 BEI - Lab 5 6 2D FFT! Un comment and f i x t h e s e 2 c a l l s to FFT! c a l l f t 1 a x i s (?,?,?,?,?)! c a l l f t 2 a x i s (?,?,?,?,?) do iw = 1, n1 w = ( iw 1) dw + w0 do i k = 1, n2 k = ( i k 1) dk + k0!#####! Fix t h i s so t h e r e i s a dip f i l t e r, p r e f e r a b l y with a taper!##### wt = 1. 0 data ( iw, i k ) = data ( iw, i k ) wt! Un comment and f i x t h e s e 2 c a l l s to i n v e r s e FFT c a l l s r i t e ( fk domain, data, n1 n2 8)! c a l l f t 1 a x i s (?,?,?,?,?)! c a l l f t 2 a x i s (?,?,?,?,?) c a l l s r i t e ( out, data, n1 n2 8) return module f o u r i e r i m p l i c i t c o n t a i n s none s ubroutine f t 1 a x i s ( adj, sign1, n1, n2, cx ) i n t e g e r i2, adj, n1, n2 complex cx ( n1, n2 ) r e a l s i g n 1 do i 2= 1, n2 c a l l f t h ( adj, sign1, 1, n1, cx ( 1, i 2 ) ) s ubroutine f t 2 a x i s ( adj, sign2, n1, n2, cx ) i n t e g e r i1, adj, n1, n2 complex cx ( n1, n2 ) r e a l s i g n 2 do i 1= 1, n1 c a l l f t h ( adj, sign2, n1, n2, cx ( i1, 1 ) ) s ubroutine f t h ( adj, sign, m1, n12, cx ) i n t e g e r i, adj, m1, n12 r e a l s i g n
7 BEI - Lab 5 7 2D FFT complex cx (m1, n12 ) complex temp ( n12 )! implement f t h using f t u s ubroutine f t u ( s i g n i, nx, cx )! complex f o u r i e r transform with unitary s c a l i n g!! 1 nx s i g n i 2 pi i ( j 1) (k 1)/nx! cx ( k ) = sum cx ( j ) e! s q r t ( nx ) j=1 for k = 1, 2,..., nx=2 i n t e g e r! i n t e g e r nx, i, j, k, m, i s t e p r e a l s i g n i, s c a l e, arg complex cx ( nx ), cmplx, cw, cdel, ct c a l l pad2 ( nx ) s c a l e = 1. / s q r t ( 1. nx ) do i= 1, nx cx ( i ) = cx ( i ) s c a l e j = 1 k = 1 do i= 1, nx i f ( i<=j ) then ct = cx ( j ) cx ( j ) = cx ( i ) cx ( i ) = ct m = nx/2 do while ( j>m. and. m>1) j = j m m = m/2! && means.and. j = j+m do i s t e p = 2 k cw = 1. arg = s i g n i / k c d e l = cmplx ( cos ( arg ), s i n ( arg ) ) do m= 1, k do i= m, nx, i s t e p ct=cw cx ( i+k ) cx ( i+k)=cx ( i ) ct cx ( i )=cx ( i )+ ct cw = cw c d e l k = i s t e p i f ( k>=nx ) then
8 BEI - Lab 5 8 2D FFT e x i t = s ubroutine pad2 ( nx ) i n t e g e r : : nx, ni ni=1 do while ( nx > ni ) ni=ni 2 nx=ni end module
Lab 4: Kirchhoff migration (Matlab version)
Due Date: Oktober 29th, 2012 TA: Mandy Wong (mandyman@sep.stanford.edu) Lab 4: Kirchhoff migration (Matlab version) Robert U. Terwilliger 1 ABSTRACT In this computer exercise you will modify the Kirchhoff
More informationOil and Gas Research Institute Seismic Analysis Center Faults Detection Using High-Resolution Seismic Reflection Techniques
Oil and Gas Research Institute Seismic Analysis Center Faults Detection Using High-Resolution Seismic Reflection Techniques Ghunaim T. Al-Anezi (KACST) March 2013 1 Objectives The objective of the survey
More informationStanford Exploration Project, Report SERGEY, November 9, 2000, pages 683??
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 683?? 682 Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 683?? Velocity continuation by spectral methods Sergey
More informationLab 1-1-D Gradient Operator
Due Date: 17:00, Friday, October 7, 2011 TA: Yunyue Elita Li (yunyue.li@sep.stanford.edu) Lab 1-1-D Gradient Operator William of Orange 1 ABSTRACT The gradient of a continuous function of more than one
More informationReflection Seismic Method
Reflection Seismic Method Data and Image sort orders; Seismic Impedance; -D field acquisition geometries; CMP binning and fold; Resolution, Stacking charts; Normal Moveout and correction for it; Stacking;
More informationAnalysis of multicomponent walkaway vertical seismic profile data
Analysis of multicomponent walkaway vertical seismic profile data Bona Wu, Don C. Lawton, and Kevin W. Hall ABSTRACT A multicomponent walkaway VSP data processed for PP and PS imaging as well to study
More informationP Wave Reflection and Refraction and SH Wave Refraction Data Processing in the Mooring, TN Area
P Wave Reflection and Refraction and SH Wave Refraction Data Processing in the Mooring, TN Area Abstract: Author: Duayne Rieger Home Institution: Slippery Rock University of Pennsylvania REU Institution:
More informationThe 2D/3D i-cube Workflow for Subsurface Imaging of Complex Structures
The 2D/3D i-cube Workflow for Subsurface Imaging of Complex Structures Öz Yilmaz CTO, GeoTomo LLC With the i-cube workflow, you circumvent the velocity uncertainty in imaging complex structures, you preserve
More informationSubsalt imaging by common-azimuth migration
Stanford Exploration Project, Report 100, April 20, 1999, pages 113 125 Subsalt imaging by common-azimuth migration Biondo Biondi 1 keywords: migration, common-azimuth, wave-equation ABSTRACT The comparison
More informationStanford Exploration Project, Report 97, July 8, 1998, pages
Stanford Exploration Project, Report 97, July 8, 1998, pages 337 342 336 Stanford Exploration Project, Report 97, July 8, 1998, pages 337 342 Short Note Factorization of cross spectra Jon Claerbout 1 keywords:
More informationWave-equation tomography for anisotropic parameters
Wave-equation tomography for anisotropic parameters Yunyue (Elita) Li and Biondo Biondi ABSTRACT Anisotropic models are recognized as more realistic representations of the subsurface where complex geological
More informationPART A: Short-answer questions (50%; each worth 2%)
PART A: Short-answer questions (50%; each worth 2%) Your answers should be brief (just a few words) and may be written on these pages if you wish. Remember to hand these pages in with your other exam pages!
More informationMultiple realizations using standard inversion techniques a
Multiple realizations using standard inversion techniques a a Published in SEP report, 105, 67-78, (2000) Robert G Clapp 1 INTRODUCTION When solving a missing data problem, geophysicists and geostatisticians
More informationEstimating a pseudounitary operator for velocity-stack inversion
Stanford Exploration Project, Report 82, May 11, 2001, pages 1 77 Estimating a pseudounitary operator for velocity-stack inversion David E. Lumley 1 ABSTRACT I estimate a pseudounitary operator for enhancing
More informationStanford Exploration Project, Report 105, September 5, 2000, pages 41 53
Stanford Exploration Project, Report 105, September 5, 2000, pages 41 53 40 Stanford Exploration Project, Report 105, September 5, 2000, pages 41 53 Short Note Multiple realizations using standard inversion
More informationMoveout, velocity, and stacking
Moveout, velocity, and stacking Jon Claerbout In this chapter we handle data as though the earth had no dipping reflectors. The earth model is one of stratified layers with velocity a (generally increasing)
More informationP S-wave polarity reversal in angle domain common-image gathers
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? P S-wave polarity reversal in angle domain common-image gathers Daniel Rosales and James Rickett 1 ABSTRACT The change in the reflection
More informationERTH3021: Exploration and Mining Geophysics
Course Profile ERTH3021: Exploration and Mining Geophysics Semester 2, 2014 Course Outline: This course builds on the general introduction provided by ERTH2020, and examines common applied-geophysical
More informationVector, Matrix, and Tensor Derivatives
Vector, Matrix, and Tensor Derivatives Erik Learned-Miller The purpose of this document is to help you learn to take derivatives of vectors, matrices, and higher order tensors (arrays with three dimensions
More informationResidual Statics using CSP gathers
Residual Statics using CSP gathers Xinxiang Li and John C. Bancroft ABSTRACT All the conventional methods for residual statics analysis require normal moveout (NMO) correction applied on the seismic data.
More informationA Petroleum Geologist's Guide to Seismic Reflection
A Petroleum Geologist's Guide to Seismic Reflection William Ashcroft WILEY-BLACKWELL A John Wiley & Sons, Ltd., Publication Contents Preface Acknowledgements xi xiii Part I Basic topics and 2D interpretation
More informationStolt residual migration for converted waves
Stanford Exploration Project, Report 11, September 18, 1, pages 1 6 Stolt residual migration for converted waves Daniel Rosales, Paul Sava, and Biondo Biondi 1 ABSTRACT P S velocity analysis is a new frontier
More informationStanford Exploration Project, Report 115, May 22, 2004, pages
Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 248 Stanford Exploration Project, Report 115, May 22, 2004, pages 249 264 First-order lateral interval velocity estimates without picking
More informationEAS 535 Laboratory Exercise Weather Station Setup and Verification
EAS 535 Laboratory Exercise Weather Station Setup and Verification Lab Objectives: In this lab exercise, you are going to examine and describe the error characteristics of several instruments, all purportedly
More informationBaseline VSP processing for the Violet Grove CO 2 Injection Site
Baseline VSP processing for Violet Grove Baseline VSP processing for the Violet Grove CO 2 Injection Site Marcia L. Couëslan, Don C. Lawton, and Michael Jones * ABSTRACT Injection of CO 2 for enhanced
More informationReversible Stolt migration
Reversible Stolt migration William A. Burnett Robert J. Ferguson ABSTRACT The nonstationary equivalent of the Fourier shift results in a general one-dimensional integral transform that applies to many
More informationNear-surface velocity characterization via RVSP and multicomponent seismic refraction experiments
Near-surface multicomponent RVSP and refraction survey Near-surface velocity characterization via RVSP and multicomponent seismic refraction experiments Darryl G. Party and Don C. Lawton Abstract A shallow
More informationMandatory Assignment 2013 INF-GEO4310
Mandatory Assignment 2013 INF-GEO4310 Deadline for submission: 12-Nov-2013 e-mail the answers in one pdf file to vikashp@ifi.uio.no Part I: Multiple choice questions Multiple choice geometrical optics
More informationInvitation to Futterman inversion
Invitation to Futterman inversion Jon Claerbout ABSTRACT A constant Q earth model attenuates amplitude inversely with the number of wavelengths propagated, so the attenuation factor is e ω (z/v)/q. We
More informationarxiv: v1 [physics.geo-ph] 23 Dec 2017
Statics Preserving Sparse Radon Transform Nasser Kazemi, Department of Physics, University of Alberta, kazemino@ualberta.ca Summary arxiv:7.087v [physics.geo-ph] 3 Dec 07 This paper develops a Statics
More informationChapter 7. Seismic imaging. 7.1 Assumptions and vocabulary
Chapter 7 Seismic imaging Much of the imaging procedure was already described in the previous chapters. An image, or a gradient update, is formed from the imaging condition by means of the incident and
More informationSub Basalt Imaging Using Low Frequency Processing and Angle stack In Saurashtra Region, Western India
5th Conference & Exposition on Petroleum Geophysics, Hyderabad-2004, India PP 90-94 Sub Basalt Imaging Using Low Frequency Processing and Angle stack In Saurashtra Region, Western India A.K.Sarkar & R.Mohan
More informationApplications of finite-difference modelling to coalscale seismic exploration
Applications of finite-difference modelling to coalscale seismic exploration Shaun Strong 1,2, Troy Peters 1 1. Velseis Pty Ltd 2 University of Queensland Introduction Geological environments with significant
More informationMadagascar tutorial: Field data processing
Madagascar tutorial: Field data processing Maurice the Aye-Aye ABSTRACT In this tutorial, you will learn about multiple attenuation using parabolic Radon transform (Hampson, 1986). You will first go through
More informationH Seismic imaging in and around salt bodies: problems and pitfalls. Ian Jones 1, Ian Davison 2. Summary
H1-2-3 Seismic imaging in and around salt bodies: problems and pitfalls Ian Jones 1, Ian Davison 2 1 ION, Chertsey, UK, 2 EarthMoves, Egham, UK Summary Salt movement often results in steeply-dipping complex
More informationThe effect of anticlines on seismic fracture characterization and inversion based on a 3D numerical study
The effect of anticlines on seismic fracture characterization and inversion based on a 3D numerical study Yungui Xu 1,2, Gabril Chao 3 Xiang-Yang Li 24 1 Geoscience School, University of Edinburgh, UK
More informationElastic versus acoustic inversion for marine surveys
Elastic versus acoustic inversion for marine surveys Item Type Article Authors Mora, Peter; Wu, Zedong Citation Mora P, Wu Z (2018) Elastic versus acoustic inversion for marine surveys. Geophysical Journal
More informationAVO inversion in V (x, z) media
Stanford Exploration Project, Report 97, July 8, 998, pages 75 94 AVO inversion in V (x, z) media Yalei Sun and Wenjie Dong keywords:.5-d Kirchhoff integral, AVO inversion, fluid-line section ABSTRACT
More informationThe course web site s notes entitled Switches, Gates and Circuits which can be found at will be useful to you throughout the lab.
Circuits Lab Names: Objectives Learn how circuits work, different properties of circuits, and how to derive truth tables, circuits from functions and simply circuits using Boolean circuit equivalence..
More informationTime-migration velocity analysis by velocity continuation
GEOPHYSICS, VOL. 68, NO. 5 (SEPTEMBER-OCTOBER 2003); P. 1662 1672, 15 FIGS. 10.1190/1.1620640 Time-migration velocity analysis by velocity continuation Sergey Fomel ABSTRACT Time-migration velocity analysis
More informationChapter 1. Introduction EARTH MODEL BUILDING
Chapter 1 Introduction Seismic anisotropy in complex earth subsurface has become increasingly important in seismic imaging due to the increasing offset and azimuth in modern seismic data. To account for
More informationINTRODUCTION. The depth-variable-velocity prestack migration in offset-midpoint coordinates (Yilmaz, 1979) is formulated as
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 58 Short Note Prestack migration by split-step DSR Alexander M. Popovici 1 INTRODUCTION The DSR prestack migration equation, though defined
More informationProcessing and interpretation of 2d land seismic reflection survey over bishop wood Yorkshire, England
Vol. 7(8), pp. 172-184, July, 2013 DOI: 10.5897/JGMR2013.0186 ISSN 2006 9766 2013 Academic Journals http://www.academicjournals.org/jgmr Journal of Geology and Mining Research Review Processing and interpretation
More informationGEO 365N/384S Seismic Data Processing Computational Assignment 5
GEO 365N/384S Seismic Data Processing Computational Assignment 5 Team: Longhorns ABSTRACT In this assignment, you will experiment with noise removal using different kinds of the Radon transform: 1. Linear
More informationAmplitude preserving AMO from true amplitude DMO and inverse DMO
Stanford Exploration Project, Report 84, May 9, 001, pages 1 167 Amplitude preserving AMO from true amplitude DMO and inverse DMO Nizar Chemingui and Biondo Biondi 1 ABSTRACT Starting from the definition
More informationDeconvolution imaging condition for reverse-time migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 83 96 Deconvolution imaging condition for reverse-time migration Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT The reverse-time
More informationStructural Geology, GEOL 330 Spring 2012
Developing the Magic Eye for folds This lab exercise is designed to get you thinking about the chronology of structural processes and and the resultant map patterns in areas with flat topography. You may
More informationInterval anisotropic parameters estimation in a least squares sense Case histories from West Africa
P-263 Summary Interval anisotropic parameters estimation in a least squares sense Patrizia Cibin*, Maurizio Ferla Eni E&P Division (Milano, Italy), Emmanuel Spadavecchia - Politecnico di Milano (Milano,
More informationTh P7 02 A Method to Suppress Salt-related Converted Wave Using 3D Acoustic Modelling
Th P7 0 A Method to Suppress Salt-related Converted Wave Using 3D Acoustic Modelling J. Kumar* (Petroleum Geo-Services), M. Salem (ENI E&P), D.E. Cegani (ENI E&P) Summary Converted waves can be recorded
More informationSeismic reflection notes
1 of 13 20/02/2006 3:29 PM Seismic reflection notes Introduction This figure was modified from one on the Lithoprobe slide set. LITHOPROBE is "probing" the "litho" sphere of our continent, the ground we
More informationTrace balancing with PEF plane annihilators
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 333 Short Note Trace balancing with PEF plane annihilators Sean Crawley 1 INTRODUCTION Land seismic data often shows large variations in energy
More informationStepwise Multiple Elimination using Linear Transforms: An alternative approach to SRME for stacking multiple-free near offsets in the Barents Sea
: An alternative approach to SRME for stacking multiple-free near offsets in the Barents Sea Anthony Hardwick*, Carl Lang and Bent Kjølhamar, TGS. Summary A novel, de-aliased method called Stepwise Multiple
More informationSummary. Introduction
: Red Lodge, Montana Jingqiu Huang*, University of Houston, and Joe Wong, University of Calgary. Summary We report on a suite of geophysical surveys conducted on glacial sediments near Red Lodge, Montana.
More informationTh Guided Waves - Inversion and Attenuation
Th-01-08 Guided Waves - Inversion and Attenuation D. Boiero* (WesternGeco), C. Strobbia (WesternGeco), L. Velasco (WesternGeco) & P. Vermeer (WesternGeco) SUMMARY Guided waves contain significant information
More informationSolving with Absolute Value
Solving with Absolute Value Who knew two little lines could cause so much trouble? Ask someone to solve the equation 3x 2 = 7 and they ll say No problem! Add just two little lines, and ask them to solve
More informationANGLE-DEPENDENT TOMOSTATICS. Abstract
ANGLE-DEPENDENT TOMOSTATICS Lindsay M. Mayer, Kansas Geological Survey, University of Kansas, Lawrence, KS Richard D. Miller, Kansas Geological Survey, University of Kansas, Lawrence, KS Julian Ivanov,
More informationBorehole Seismic Monitoring of Injected CO 2 at the Frio Site
Borehole Seismic Monitoring of Injected CO 2 at the Frio Site * Daley, T M (tmdaley@lbl.gov), Lawrence Berkeley National Lab., 1 Cyclotron Rd, Berkeley, CA 94720 Myer, L (lrmyer@lbl.gov), Lawrence Berkeley
More informationDipping Layer Refraction Problem Moveout and Coincident Source-Receiver Format
Environmental and Exploration Geophysics II Dipping Layer Refraction Problem Moveout and Coincident Source-Receiver Format tom.h.wilson tom.wilson@mail.wvu.edu Department of Geology and Geography West
More informationLab 3 Acceleration. What You Need To Know: Physics 211 Lab
b Lab 3 Acceleration Physics 211 Lab What You Need To Know: The Physics In the previous lab you learned that the velocity of an object can be determined by finding the slope of the object s position vs.
More informationThe seismic response to strong vertical velocity change
10 th Biennial International Conference & Exposition P-073 The seismic response to strong vertical velocity change Ian F. Jones, ION GX Technology Summary Conventional seismic data processing, whether
More informationMatrix formulation of adjoint Kirchhoff datuming
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 377 Matrix formulation of adjoint Kirchhoff datuming Dimitri Bevc 1 ABSTRACT I present the matrix formulation of Kirchhoff wave equation datuming
More informationTime domain passive seismic processing at Valhall
Stanford Exploration Project, Report 125, January 16, 2007, pages 1?? Time domain passive seismic processing at Valhall Brad Artman ABSTRACT Passive recordings from an array of 2500 four component receiver
More informationAn Alternative Approach to Process the Wide-Angle Reflection Data By Pre-stack Compositing of Gathers for Sub-basalt Imaging
5th Conference & Exposition on Petroleum Geophysics, Hyderabad-2004, India PP 288-293 An Alternative Approach to Process the Wide-Angle Reflection Data By Pre-stack Compositing of Gathers for Sub-basalt
More informationExamples of prestack depth migration in TI media
Eamples of prestack depth migration in TI media Robert J. Ferguson and Gary F. Margrave ABSTRACT Wave field etrapolation by nonstationary phase shift can be formulated to allow velocity variation with
More informationGeorge Mason University ECE 201: Introduction to Signal Analysis Spring 2017
Assigned: March 20, 2017 Due Date: Week of April 03, 2017 George Mason University ECE 201: Introduction to Signal Analysis Spring 2017 Laboratory Project #6 Due Date Your lab report must be submitted on
More informationProbabilistic Graphical Models Homework 2: Due February 24, 2014 at 4 pm
Probabilistic Graphical Models 10-708 Homework 2: Due February 24, 2014 at 4 pm Directions. This homework assignment covers the material presented in Lectures 4-8. You must complete all four problems to
More informationLecture Notes (Math 90): Week IX (Tuesday)
Lecture Notes (Math 90): Week IX (Tuesday) Alicia Harper November 13, 2017 Existence theorems and monotonic functions 1 The Story 1.1 Two more existence theorems! Recall the intermediate value theorem:
More informationSatish Singh* (IPG Paris, France, Tim Sears (British Gas, UK), Mark Roberts (IPG Paris, Summary. Introduction P - 92
P - 92 Fine-Scale P- and S-Wave Velocities From Elastic Full Waveform Inversion of Multi-Component and Time-Lapse Data: Future of Quantitative Seismic Imaging Satish Singh* (IPG Paris, France, singh@ipgp.jussieu.fr),
More informationDownloaded 11/08/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Velocity tomography based on turning waves and finite-frequency sensitivity kernels Xiao-Bi Xie*, University of California at Santa Cruz, Jan Pajchel, Statoil, Norway Summary We present a velocity tomography
More informationVollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich
Seminar des FA Ultraschallprüfung Vortrag 1 More info about this article: http://www.ndt.net/?id=20944 Vollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich Thomas
More informationDownloaded 05/01/17 to Redistribution subject to SEG license or copyright; see Terms of Use at
Mapping Imbricate Structures in the Thrust Belt of Southeast Turkey by Large-Offset Seismic Survey Oz Yilmaz*, Anatolian Geophysical, Istanbul, Turkey; and Serdar Uygun, Ali Ölmez, and Emel Çalı, Turkish
More informationData examples of logarithm Fourier domain bidirectional deconvolution
Data examples of logarithm Fourier domain bidirectional deconvolution Qiang Fu and Yi Shen and Jon Claerbout ABSTRACT Time domain bidirectional deconvolution methods show a very promising perspective on
More information22.4 Choosing a Method for Solving Quadratic Equations
Name Class Date.4 Choosing a Method for Solving Quadratic Equations Essential Question: How can you choose a method for solving a given quadratic equation? Resource Locker Explore 7 x - 3x - 5 = 0 Comparing
More informationDrift time estimation by dynamic time warping
Drift time estimation by dynamic time warping Tianci Cui and Gary F. Margrave CREWES, University of Calgary cuit@ucalgary.ca Summary The drift time is the difference in traveltime at the seismic frequency
More informationElements of 3D Seismology Second Edition
Elements of 3D Seismology Second Edition Copyright c 1993-2003 All rights reserved Christopher L. Liner Department of Geosciences University of Tulsa August 14, 2003 For David and Samantha And to the memory
More informationThere is no pure P- or S-wave land seismic source André J.-M. Pugin*, Geological Survey of Canada, and Oz Yilmaz, Anatolian Geophysical
There is no pure P- or S-wave land seismic source André J.-M. Pugin*, Geological Survey of Canada, and Oz Yilmaz, Anatolian Geophysical Summary: We conducted a field experiment at a soil site near Ottawa
More informationHigh Resolution Seismic for Minerals
High Resolution Seismic for Minerals Greg Turner HiSeis Pty Ltd MEGWA April 2013 WELCOME HiSeis Pty Ltd, Suite 4 Enterprise Unit 3, 9 DeLaeter Way, Bentley, WA 6102 Ph: +61 8 9470 9866 www.hiseis.com Used
More informationLab 4: Introduction to Signal Processing: Fourier Transform
Lab 4: Introduction to Signal Processing: Fourier Transform This laboratory requires the following equipment: Matlab The laboratory duration is approximately 3 hours. Although this laboratory is not graded,
More informationLAB 5: THE EARTHQUAKE CYCLE
NAME: LAB TIME: LAB 5: THE EARTHQUAKE CYCLE This lab will introduce you to the basic quantitative concepts of the earthquake cycle for a vertical strike-slip fault. Most of your time will be spent calculating
More informationA modeling and migration study of fault shadows: A synthetic example
Fault shadow imaging A modeling and migration study of fault shadows: A synthetic example Sitamai W. Ajiduah and Gary F. Margrave This report presents a study that uses 2D finite difference modeling and
More informationComplex-beam Migration and Land Depth Tianfei Zhu CGGVeritas, Calgary, Alberta, Canada
Page 1 of 10 Home Articles Interviews Print Editions Complex-beam Migration and Land Depth Tianfei Zhu CGGVeritas, Calgary, Alberta, Canada DECEMBER 2012 FOCUS ARTICLE Summary Gaussian-beam depth migration
More informationEE3210 Lab 3: Periodic Signal Representation by Fourier Series
City University of Hong Kong Department of Electronic Engineering EE321 Lab 3: Periodic Signal Representation by Fourier Series Prelab: Read the Background section. Complete Section 2.2(b), which asks
More informationFull Waveform Inversion (FWI) with wave-equation migration. Gary Margrave Rob Ferguson Chad Hogan Banff, 3 Dec. 2010
Full Waveform Inversion (FWI) with wave-equation migration (WEM) and well control Gary Margrave Rob Ferguson Chad Hogan Banff, 3 Dec. 2010 Outline The FWI cycle The fundamental theorem of FWI Understanding
More informationInvestigating the Barnett Shale
Investigating the Barnett Shale An Integrated Workflow from Petrophysics to Visualisation and Seismic Decomposition Vision for Energy The Barnett Shale 1981: Barnett Shale discovery 8,000+ wells to date
More informationChapter 7: Reflection Seismology Homework Solutions (Jan. 2010)
Chapter 7: eflection Seismology Homework Solutions (Jan. 200). Why do marine seismic reflection surveys not record (a) S waves? (b) refracted rays? 2 μ a) For ideal fluid, μ=0, thus, v s = = 0 ρ b) eflection
More informationReciprocity of tensor wavefields
Chapter 4 Reciprocity of tensor wavefields This chapter consists of two parts: one shows the validity of the reciprocity principle for anisotropic heterogeneous media using a properly designed numerical
More informationVelocity structure of the methane hydrate formation at the Blake Outer Ridge
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 123 Velocity structure of the methane hydrate formation at the Blake Outer Ridge Christine Ecker 1 ABSTRACT Seismic analysis of data from the
More informationEdinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains
Frequency-dependent AVO attribute: theory and example Xiaoyang Wu, 1* Mark Chapman 1,2 and Xiang-Yang Li 1 1 Edinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains Road, Edinburgh
More informationProject: Equilibrium Simulation
UNIT 02: BCLN CHEMISTRY 12 - Rev. July, 2015 Project: Simulation Potential Credits: /20 Name: Goal: In this project you will explore the concept of reversible reactions and calculating K eq using some
More informationSynthesized shot gathers from Valhall
Chapter 3 Synthesized shot gathers from Valhall SUMMARY BP made available assive recordings from an array of 2300 four component receiver stations at the Valhall reservoir in the Norwegian North Sea from
More informationDesigning Information Devices and Systems I Spring 2018 Lecture Notes Note Introduction to Linear Algebra the EECS Way
EECS 16A Designing Information Devices and Systems I Spring 018 Lecture Notes Note 1 1.1 Introduction to Linear Algebra the EECS Way In this note, we will teach the basics of linear algebra and relate
More informationOn non-stationary convolution and inverse convolution
Stanford Exploration Project, Report 102, October 25, 1999, pages 1 137 On non-stationary convolution and inverse convolution James Rickett 1 keywords: helix, linear filtering, non-stationary deconvolution
More informationStatics preserving projection filtering Yann Traonmilin*and Necati Gulunay, CGGVeritas
Yann Traonmilin*and Necati Gulunay, CGGVeritas Summary Projection filtering has been used for many years in seismic processing as a tool to extract a meaningful signal out of noisy data. We show that its
More informationLecture 17: Floyd-Hoare Logic for Partial Correctness
Lecture 17: Floyd-Hoare Logic for Partial Correctness Aims: To look at the following inference rules Page 1 of 9 sequence; assignment and consequence. 17.1. The Deduction System for Partial Correctness
More informationImprovement of stacking image by anisotropic velocity analysis using P-wave seismic data
P-46 Improvement of stacking image by anisotropic velocity analysis using P-wave seismic data Laxmidhar Behera*, and Prakash Khare, National Geophysical Research Institute (NGRI), Hyderabad Summary Anisotropy
More informationPartial Fraction Decomposition Honors Precalculus Mr. Velazquez Rm. 254
Partial Fraction Decomposition Honors Precalculus Mr. Velazquez Rm. 254 Adding and Subtracting Rational Expressions Recall that we can use multiplication and common denominators to write a sum or difference
More informationEOS 350 MIDTERM OCT 4, 2013 STUDENT NAME: TEAM #:
EOS 350 MIDTERM OCT 4, 2013 STUDENT NAME: TEAM #: Some equations which may, or may not, be useful: Distance from sensor to a dipole z ~ x ½, Distance to line of dipoles z ~ 0.75x ½ B = μh, M = κh Seismic
More informationJumping on a scale. Data acquisition (TI 83/TI84)
Jumping on a scale Data acquisition (TI 83/TI84) Objective: In this experiment our objective is to study the forces acting on a scale when a person jumps on it. The scale is a force probe connected to
More informationLecture 20: Discrete Fourier Transform and FFT
EE518 Digital Signal Processing University of Washington Autumn 2001 Dept of Electrical Engineering Lecture 20: Discrete Fourier Transform and FFT Dec 10, 2001 Prof: J Bilmes TA:
More informationDownloaded 09/17/13 to Redistribution subject to SEG license or copyright; see Terms of Use at Log data.
Downloaded 9/17/13 to 99.186.17.3. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ Extracting polar antropy parameters from mic data and well logs Rongrong
More information