Surface wave analysis for P- and S-wave velocity models
|
|
- Damian Newman
- 5 years ago
- Views:
Transcription
1 Distinguished Lectures in Earth Sciences, Napoli, 24 Maggio 2018 Surface wave analysis for P- and S-wave velocity models Laura Valentina Socco, Farbod Khosro Anjom, Cesare Comina, Daniela Teodor POLITECNICO DI TORINO UNIVERSITA DI TORINO
2 What am I going to talk about Surface wave propagation in layered systems and the concept of geometrical dispersion; The inversion of surface wave dispersion curves for 1D VS model estimation: potentials and limitations; Some examples of applications in laterally varying sites; A change in paradigm: from interval velocity to time-average velocity; A new approach to subsurface characterization: from inversion to data transform; Sensitivity to Poisson s ratio and VP estimation.
3 Surface waves: definitions There are many kinds of surface waves (guided waves) that can be exploited in geophysics: Rayleigh waves (1885): land data, ground surface, in line polarization, always generated; Scholte waves (1947): marine data, seabed, in line polarization, always generated; P guided waves: land data, ground surface, in line polarization; Love waves (1911): land/marine data, ground surface/seabed, horizontal cross line polarization; Lamb waves (1904): on solid plates; Stoneley waves (1924): in boreholes, along borehole walls;
4 Surface wave propagation in the homogeneous medium Rayleigh waves travel in a limited layer close to the free surface; The thickness of the layer depends on the wavelength of the propagating harmonics; Hence in a homogeneous medium high frequencies (short l) propagate close to the surface while low frequencies (long l) reach greater depth. Low frequency High frequency
5 Surface wave propagation in the homogeneous medium In a homogeneous halfspace, the different frequency components of the seismic signal produce different wavelengths, which propagates at different depths, at the Rayleigh wave velocity of the medium.
6 Surface wave propagation in the homogeneous medium In a homogeneous halfspace, the different frequency components of the seismic signal produce different wavelengths, which propagates at different depths, at the Rayleigh wave velocity of the medium (about 0.9 VS).
7 Surface wave propagation in the layered medium geometrical dispersion In a layered halfspace, the different frequency components of the seismic signal produce different wavelengths, which propagates at different depths, with different velocities depending from the depth.
8 Surface wave propagation in the layered medium
9 Surface wave propagation in the layered medium Rayleigh waves appear as high energy (amplitude) signals; body waves The effect of geometrical dispersion in time domain is the tipycal cone shape of the wavelength; The velocity of propagation is lower than body waves and varies for the different harmonics. Surface waves
10 Surface wave propagation in the layered medium Surface wave propagation can be modeled by for a stack of linear homogenous elastic layers. Layer parameters are: S-wave velocity, Poisson s ratio (or P-wave velocity), density and thickness. When we solve the forward problem we have to find the zeros of the Rayleigh secular function and they provides couple or f-v (or f-k) values which are possible solutions of the propagation.
11 Surface wave propagation in the layered medium V S =100m/s V P =200m/s V S =300m/s V P =600m/s 300 H = 5m Cut-off f-v frequencies f-k wavenumber [rad/m] frequency [Hz] phase velocity [m/s] Dispersion curve frequency [Hz] Higher modes Fundamental mode
12 Characterization : inversion of geometrical dispersion A brief introduction Vertical particle about motion surface Phase waves velocity and V R their V R use for near Vsurface characterization; S1 V R = l f? V S2 > V S1 Wavelength l V S3 > V S2 Velocity profile Z Z Short wavelength High frequency Long wavelength Low frequency Frequency f Dispersion Curve Experimental INVERSE PROBLEM
13 Seismology Since the early 60 to map the earth crust. Applications Exploration geophysics Only in the last decade for various targets. Geotechnical engineering Since the early 80 to estimate V S profiles. [Yilmaz et al., 2006] [Shapiro and Ritzwoller, 2002] [Brown, Boore and Stokoe, 2002] [Park et al., 2005]
14 The surface wave method ACQUISITION RAW DATA k f DISPERSION PROCESSING V S INVERSION z VELOCITY MODEL
15 The surface wave method ACQUISITION RAW DATA k f DISPERSION PROCESSING V S INVERSION z VELOCITY MODEL
16 The surface wave method time [s] ACQUISITION 2D FFT offset [m] RAW DATA k f DISPERSION PROCESSING 2p f V( f) = k V S INVERSION z VELOCITY MODEL
17 The surface wave method Shear Wave Velocity (m/s) RAW DATA ACQUISITION Depth (m) z k V S f DISPERSION VELOCITY MODEL PROCESSING INVERSION phase velocity, m/s experimental frequency, Hz
18 Summary of SW methods: potentials Surface waves are easy to gather and to recognize; Robust; Available in any seismic gather no matter which is the purpose of acquisition (seismic reflection, P-wave refraction); Allow S-wave velocity to be estimated.
19 Summary of SW methods: limitations Strong non linearity of the invers problem à the result is non unique (no matter if we use global or local search strategies) and the sensitivity to model parameters decreases with depth limiting the investigation depth ; 1D; VP is not estimated.
20 Laterally heterogeneous media The seismogram is windowed in the spatial domain with several asymmetric Gaussian windows whose maxima span the survey line
21 Laterally heterogeneous media The seismogram is windowed in the spatial domain with several asymmetric Gaussian windows whose maxima span the survey line [Bergamo, P., D. Boiero, L. V. Socco, 2012, Retrieving 2D structures from surface wave data by means of a space-varying spatial windowing: Geophysics, 77, 4, EN39-EN51, doi: /GEO ]
22 Laterally heterogeneous media
23 Laterally heterogeneous media
24 Laterally heterogeneous media
25 Spatially Constrained Inversion Full model built up of a number of 1D shear wave velocity models, model parameters are shear wave velocities and depths; Lateral constraints couple the different 1D-models. The constraints consist of the spatially dependent covariance between the model parameters... and can be considered as a priori information on the geological variation in the area; LCI allows for smooth transitions in model parameters along the profile; All data are inverted simultaneously as one system. [AUKEN and CHRISTIANSEN, 2004]; [SOCCO, D. BOIERO, S. FOTI, R. WISÉN, 2009, Laterally constrained inversion of ground roll from seismic reflection records, Geophysics, 74, 6, G35-G45]
26 Field case: shallow bedrock mapping Land-streamer 48 ch 1.25m 55 ch 2.50m Vibroseis Sweep: 7-80 Hz Length: 15 s
27 Field case: shallow bedrock mapping
28 Field case: shallow bedrock mapping
29 Field case: shallow bedrock mapping
30 Field case: shallow bedrock mapping [Boiero, D., L.V. Socco, S. Stocco, R. Wisén, 2013, Bedrock mapping in shallow environments using surface-wave analysis: The Leading Edge, June, ]
31 Summary: SW method in laterally varying media Several strategies can be implemented to mitigate the limitations of surface wave method; 2D sites can be investigated and the results are in agreement with body wave methods; The inversion process remains a critical aspect of the workflow (deterministic methods are non unique, stochastic methods are computationally heavy). still we estimate only VS and assume VP L
32 What am I going to show? We move from interval velocity to time-average velocity We identify a relationship between wavelength and investigation depth and we use it to get VS without inverting; This relationship is sensitive to Poisson s ratio and we can hence use it to estimate directly VP;
33 The concept of time-average velocity Time-average velocity is used in many applications ranging from seismic hazard (VS,30) to static corrections (provides directly the one-way time). V Sz, = å n n å h i h V i Si h 1, V S1 h i, V Si z h n, V Sn
34 The wavelength-depth relationship W/D We search a relationship between the time-average velocity at a given depth z (V S,z ) and the surface wave velocity at a given wavelength.
35 The wavelength-depth relationship W/D We search a relationship between the weighted average velocity at a given depth z (V S,z ) and the surface wave velocity at a given wavelength. Linear regression
36 The meaning of the W/D relationship [Karray, M., and G. Lefebvre, 2008, Significance and Evaluation of Poisson s ratio, Canadian Geotechnical Journal, 45, ] [Pelekis, P.C., and G.A. Athanasopoulos, 2011, An overview of surface wave methods and a reliability study of a simplified inversion technique, Soil Dynamics and Earthquake Engineering,31, ]
37 The W/D sensitivity to Poisson If this is true, also our W/D relationship should be sensitive to Poisson s ratio. Let s see what happens with a constant Poisson with depth (0.33). Depth or Wavelength [m] VS VP SW DC Estimated VSz true VSz Estimated VPz True VPz V [m/s] Wavelength [m] Depth [m]
38 The W/D sensitivity to Poisson The sensitivity to Poisson s has a strong physical meaning and is related to the pattern of vertical displacement with depth. Depth [m] normalized vertical displacement [-] wavelength [m] Poisson depth [m]
39 The W/D sensitivity to Poisson Hence we can use this sensitivity to estimate the Poisson ratio for the different depth and use it to transform the VSz model into a VPz model. Depth or Wavelength [m] VS VP SW DC Estimated VSz true VSz Estimated VPz True VPz V [m/s] Wavelength [m] Wavelength [m] Depth [m]
40 The W/D sensitivity to Poisson Hence we can use this sensitivity to estimate the Poisson ratio for the different depth and use it to transform the VSz model into a VPz model. Depth or Wavelength [m] V [m/s] VS VP SW DC Estimated VSz true VSz Estimated VPz True VPz Wavelength [m] Wavelength [m] Depth [m] [m]
41 I have shown that A relationship exists between VSz at various depth and SW phase velocity at various wavelengths (W/D); The W/D relationship is sensitive to Poisson s ratio; If the Poisson s ratio is constant with depth this sensitivity is enough to estimate the correct Poisson s ratio and transform the VSz into VPz. but what happens if Poisson ratio varies with depth?
42 The apparent Poisson ratio Here we use the same VS model of the previous example but with variable Poisson with depth (0.33 for the first 4 layers and 0.2 for the deeper ones). Depth or Wavelength [m] VS 10 VP true VPz 20 SW DC true VSz V [m/s] ævpz ö ç VSz n z = è ø 2 2 ævpz ö ç - 1 èvsz ø 2 Depth [m] Poisson [-]
43 The apparent Poisson ratio Here we use the same VS model of the previous example but with variable Poisson with depth (0.33 for the first 4 layers and 0.2 for the deeper ones). Depth or Wavelength [m] VS 10 VP true VPz 20 SW DC true VSz V [m/s] Wavelength [m] Depth [m]
44 The apparent Poisson ratio Here we use the same VS model of the previous example but with variable Poisson with depth (0.33 for the first 4 layers and 0.2 for the deeper ones). Depth or Wavelength [m] VS 10 VP true VPz 20 SW DC true VSz V [m/s] Wavelength [m] Wavelength [m] Depth Depth [m] [m]
45 The apparent Poisson ratio Here we use the same VS model of the previous example but with variable Poisson with depth (0.33 for the first 4 layers and 0.2 for the deeper ones). Depth or Wavelength [m] VS 10 VP true VPz 20 SW DC true VSz V [m/s] Wavelength [m] Depth [m]
46 The apparent Poisson ratio Here we use the same VS model of the previous example but with variable Poisson with depth (0.33 for the first 4 layers and 0.2 for the deeper ones). Depth or or Wavelength [m] [m] V [m/s] VS VS 10 VP VP true VPz true VPz data4 SW DC 20 estimated true VPz VSz Estimated VSz true VSz Depth [m] apparent Poisson 10 estimated Poisson Poisson [-]
47 Retrieving the interval velocity If we can estimate the time-average velocity at any depth we can apply a Dix-like equation and retrieve the interval velocity at all depths à we have the V S and the VP profile. V Si = z zi V i - z z - V i-1 i-1 Si, Si, -1
48 I have shown that Using the W/D and the apparent Poisson s ratio we can estimate VSz and VPz profile; By applying a Dix-like relationship we can transform the time-average velocities into interval velocities;. but To estimate W/D we need the VSz model, so we are in a kind of loop and we need to find a strategy to make all this useful in practical applications.
49 Real world applications To make the outlined concepts (W/D and apparent n) of practical use we operates according to the following workflow: 1. Estimate W/D for one reference model along a line 2. Apply it to all the DCs along the line à VSz 3. Retrieve apparent n and transform VSz à VPz 4. Apply dix-like relation à VS and VP
50 Real world example
51 VS and VP benchmark DHT VS VS model from LCI of dispersion curves along the line- DHT VP VP model from travel time tomography [Socco, L.V., D. Boiero, C. Comina, S. Foti, R. Wisén, 2008, Seismic Characterization of an Alpine Site: Near Surface Geophysics, 6, ]
52 The reference model To define the reference model we select one of the broadband DC obtained by merging active and passive data and we invert it with a MC inversion (Poisson free to vary). Wavelength [m] Phase Velocity [m/s] Depth or Wavelength [m] VS MC inversion 10 VSz SW DC V [m/s]
53 The W/D relationship and apparent Poisson s ratio The W/D is computed for the reference model and, using the reference VS model the W/D for different Poisson s ratio values is computed and compared with the experimental one to find the apparent Poisson s ratio Poisson [-] Wavelength [m] Depth [m] Depth [m] 50 60
54 Results: VSz The comparison between VSz from LCI results and the VSz estimated with W/D data transform has difference mostly below 10%. [Socco L.V., C. Comina, F. Khosro Anjom, 2017, Time-average velocity estimation through surface-wave analysis: Part 1 S-wave velocity: GEOPHYSICS, 82, 3, U49 U59]
55 Results: VPz The comparison between VPz from P-wave tomography and VPz estimated with apparent Poisson s ratio applied to the VSz results has difference mostly below 10%. [Socco L.V., C. Comina, 2017, Timeaverage velocity estimation through surface-wave analysis: Part 2 P- wave velocity: GEOPHYSICS, 82, 3, U61 U73]
56 Results: interval VS LCI From VSz The comparison between VS from LCI results and VS estimated by applying Dix-like relation to VSz data has difference mostly below 15% in the upper part of the model, but rech error of 100% in deeper part where data are coarser. Normalized difference
57 Results: interval VP P-wave tomography From VPz Normalized difference The comparison between VP from P-wave tomography and VP estimated by applying Dix-like relation to VPz data has difference mostly below 20% in the upper part of the model, but reach errors of 100% in deeper part where data are coarser.
58 Borehole benchmark About 700 m far from reference model The comparison between VS and VP local model from down-hole test, inversions (LCI and P-wave tomo) and our approximate estimation shows very good agreement for VS, and reasonable agreement for VP with lower difference from DHT at shallow depth and greater difference for deeper layers.
59 I have shown that We can apply the W/D relationship and the apparent Poisson s ratio concepts in a practical workflow and obtain a VS and VP model along a seismic line knowing only one reference VS model.. but The interval velocity has high uncertainties for noisy data. There is a strong assumption: weak lateral variations of VS and negligible lateral variations of Poisson s ratio.
60 Conclusion The SW dispersion curve is extremely sensitive to time-average velocities and a relationship (W/D) exists between V S,z at various depth and SW phase velocity at various wavelengths. This relationship can be used to estimate V S,z and, consequently a smooth approximated V S profiles for a set of similar models. The W/D is sensitive to Poisson s ratio, hence the direct estimation of V P,z and V P becomes possible. The VP and VS models obtained with this workflow have proven to be suitable as initial models for FWI being within cycle skipping at high frequency.
61 Aknowledgements European Union DG12 for SISMOVALP Project - La Salle data. Total E&P for supporting Farbod Khosro PhD grant. Italian Ministry of Research for supporting Daniela Teodor PhD grant. All the PhDs and Post-docs who contributed developing our processing and inversion codes (Claudio Strobbia, Daniele Boiero, Margherita Maraschini, Claudio Piatti, Paolo Bergamo, Flora Garofalo). Thank you for your attention
SUMMARY INTRODUCTION MOTIVATION
Isabella Masoni, Total E&P, R. Brossier, University Grenoble Alpes, J. L. Boelle, Total E&P, J. Virieux, University Grenoble Alpes SUMMARY In this study, an innovative layer stripping approach for FWI
More informationNorthing (km)
Imaging lateral heterogeneity at Coronation Field with surface waves Matthew M. Haney, Boise State University, and Huub Douma, ION Geophysical/GXT Imaging Solutions SUMMARY A longstanding problem in land
More informationGuidelines for the good practice of surface wave analysis: a product of the InterPACIFIC project
Bull Earthquake Eng (2018) 16:2367 2420 https://doi.org/10.1007/s10518-017-0206-7 ORIGINAL RESEARCH PAPER Guidelines for the good practice of surface wave analysis: a product of the InterPACIFIC project
More informationP34 Determination of 1-D Shear-Wave Velocity Profileusing the Refraction Microtremor Method
P34 Determination of 1-D Shear-Wave Velocity Profileusing the Refraction Microtremor Method E. Baniasadi* (University of Tehran), M. A. Riahi (University of Tehran) & S. Chaychizadeh (University of Tehran)
More informationMulticomponent seismic polarization analysis
Saul E. Guevara and Robert R. Stewart ABSTRACT In the 3-C seismic method, the plant orientation and polarity of geophones should be previously known to provide correct amplitude information. In principle
More informationA Dissertation Presented for the Doctor of Philosophy Degree. The University of Memphis
A NEW PROCEDURE FOR ESTIMATION OF SHEAR WAVE VELOCITY PROFILES USING MULTI STATION SPECTRAL ANALYSIS OF SURFACE WAVES, REGRESSION LINE SLOPE, AND GENETIC ALGORITHM METHODS A Dissertation Presented for
More informationGENERAL GUIDELINES FOR APPLICATION OF THE EXTENDED SUBTRACTION METHOD IN SASSI SOIL-STRUCTURE INTERACTION ANALYSIS
Transactions, SMiRT-22 GENERAL GUIDELINES FOR APPLICATION OF THE EXTENDED SUBTRACTION METHOD IN SASSI SOIL-STRUCTURE INTERACTION ANALYSIS C. C. Chin 1, Nan Deng 2, and Farhang Ostadan 3 1 Senior Engineer,
More informationTomostatic Waveform Tomography on Near-surface Refraction Data
Tomostatic Waveform Tomography on Near-surface Refraction Data Jianming Sheng, Alan Leeds, and Konstantin Osypov ChevronTexas WesternGeco February 18, 23 ABSTRACT The velocity variations and static shifts
More informationPolarization Filter by Eigenimages and Adaptive Subtraction to Attenuate Surface-Wave Noise
Polarization Filter by Eigenimages and Adaptive Subtraction to Attenuate Surface-Wave Noise Stephen Chiu* ConocoPhillips, Houston, TX, United States stephen.k.chiu@conocophillips.com and Norman Whitmore
More informationSUMMARY INTRODUCTION GROUP VELOCITY
Surface-wave inversion for near-surface shear-wave velocity estimation at Coronation field Huub Douma (ION Geophysical/GXT Imaging solutions) and Matthew Haney (Boise State University) SUMMARY We study
More informationMultichannel analysis of surface waves
GEOPHYSICS, VOL. 64, NO. 3 (MAY-JUNE 1999); P. 800 808, 7 FIGS. Multichannel analysis of surface waves Choon B. Park, Richard D. Miller, and Jianghai Xia ABSTRACT The frequency-dependent properties of
More informationModule 2 WAVE PROPAGATION (Lectures 7 to 9)
Module 2 WAVE PROPAGATION (Lectures 7 to 9) Lecture 9 Topics 2.4 WAVES IN A LAYERED BODY 2.4.1 One-dimensional case: material boundary in an infinite rod 2.4.2 Three dimensional case: inclined waves 2.5
More informationSouth Africa CO2 Seismic Program
1 South Africa CO2 Seismic Program ANNEXURE B Bob A. Hardage October 2016 There have been great advances in seismic technology in the decades following the acquisition of legacy, limited-quality, 2D seismic
More informationWS15-B02 4D Surface Wave Tomography Using Ambient Seismic Noise
WS1-B02 4D Surface Wave Tomography Using Ambient Seismic Noise F. Duret* (CGG) & E. Forgues (CGG) SUMMARY In 4D land seismic and especially for Permanent Reservoir Monitoring (PRM), changes of the near-surface
More informationSeismic Reflection Method
1 of 25 4/16/2009 11:41 AM Seismic Reflection Method Top: Monument unveiled in 1971 at Belle Isle (Oklahoma City) on 50th anniversary of first seismic reflection survey by J. C. Karcher. Middle: Two early
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,500 108,000 1.7 M Open access books available International authors and editors Downloads Our
More informationMicrotremor Array Measurements and Three-component Microtremor Measurements in San Francisco Bay Area
Microtremor Array Measurements and Three-component Microtremor Measurements in San Francisco Bay Area K. Hayashi & D. Underwood Geometrics, Inc., United States SUMMARY: Microtremor array measurements and
More informationElectromagnetic Induction
Electromagnetic Induction Recap the motivation for using geophysics We have problems to solve Slide 1 Finding resources Hydrocarbons Minerals Ground Water Geothermal Energy SEG Distinguished Lecture slide
More informationApplied Methods MASW Method
Applied Methods MASW Method Schematic illustrating a typical MASW Survey Setup INTRODUCTION: MASW a seismic method for near-surface (< 30 m) Characterization of shear-wave velocity (Vs) (secondary or transversal
More informationThis presentation was prepared as part of Sensor Geophysical Ltd. s 2010 Technology Forum presented at the Telus Convention Center on April 15, 2010.
This presentation was prepared as part of Sensor Geophysical Ltd. s 2010 Technology Forum presented at the Telus Convention Center on April 15, 2010. The information herein remains the property of Mustagh
More informationESTIMATION OF SHEAR WAVE VELOCITY PROFILES USING MICROTREMOR ARRAY EXPLORATIONS IN ISMAILIA CITY, EGYPT
ESTIMATION OF SHEAR WAVE VELOCITY PROFILES USING MICROTREMOR ARRAY EXPLORATIONS IN ISMAILIA CITY, EGYPT Mohamed Maklad MEE16704 Supervisor: Toshiaki Yokoi Takumi Hayashida ABSTRACT The Spatial Autocorrelation
More informationDesign of an Optimal High Pass Filter in Frequency Wave Number (F-K) Space for Suppressing Dispersive Ground Roll Noise from Onshore Seismic Data
Universal Journal of Physics and Application 11(5): 144-149, 2017 DOI: 10.13189/ujpa.2017.110502 http://www.hrpub.org Design of an Optimal High Pass Filter in Frequency Wave Number (F-K) Space for Suppressing
More informationGCM mapping Vildbjerg - HydroGeophysics Group - Aarhus University
GCM mapping Vildbjerg - HydroGeophysics Group - Aarhus University GCM mapping Vildbjerg Report number 06-06-2017, June 2017 Indholdsfortegnelse 1. Project information... 2 2. DUALEM-421s... 3 2.1 Setup
More informationShallow shear wave velocity structure in two sites of Khartoum, Sudan using methods of seismic dispersion and ambient noise.
Norwegian National Seismic Network Technical Report No. 25 Shallow shear wave velocity structure in two sites of Khartoum, Sudan using methods of seismic dispersion and ambient noise. Prepared by Miguel
More informationMULTI-COMPONENT ACTIVE SOURCE RAYLEIGH WAVE ANALYSIS. Gabriel Gribler. A thesis. submitted in partial fulfillment
MULTI-COMPONENT ACTIVE SOURCE RAYLEIGH WAVE ANALYSIS by Gabriel Gribler A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Geophysics Boise State University
More informationSeismic reflection method
Seismic reflection method Seismic reflection method is based on the reflections of seismic waves occurring at the contacts of subsurface structures. We apply some seismic source at different points of
More informationCorresponding Author William Menke,
Waveform Fitting of Cross-Spectra to Determine Phase Velocity Using Aki s Formula William Menke and Ge Jin Lamont-Doherty Earth Observatory of Columbia University Corresponding Author William Menke, MENKE@LDEO.COLUMBIA.EDU,
More informationAmplitude balancing for AVO analysis
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 356 Amplitude balancing for AVO analysis Arnaud Berlioux and David Lumley 1 ABSTRACT Source and receiver amplitude variations can distort
More informationUSE OF RAYLEIGH WAVE METHODS TO DETECT NEAR SURFACE CONCRETE DAMAGE
USE OF RAYLEIGH WAVE METHODS TO DETECT NEAR SURFACE CONCRETE DAMAGE R. Al Wardany 1,2, J. Rhazi 1, G. Ballivy 1, J.L. Gallias 2, K. Saleh 3 1 Research Group on NDT and Instrumentation, Université de Sherbrooke,
More informationTitleApplication of MEMS accelerometer t. AIZAWA, Takao; KIMURA, Toshinori; M Toshifumi; TAKEDA, Tetsuya; ASANO,
TitleApplication of MEMS accelerometer t Author(s) AIZAWA, Takao; KIMURA, Toshinori; M Toshifumi; TAKEDA, Tetsuya; ASANO, Citation International Journal of the JCRM ( Issue Date 2008-12 URL http://hdl.handle.net/2433/85166
More informationAdvanced Ground Investigation Techniques to Help Limit Risk or Examine Failure. Advanced Subsurface Investigations
Advanced Ground Investigation Techniques to Help Limit Risk or Examine Failure Overview Introduction What is geophysics? Why use it? Common Methods Seismic Ground Radar Electrical Case Studies Conclusion
More informationUnderstanding Seismic Amplitudes
Understanding Seismic Amplitudes The changing amplitude values that define the seismic trace are typically explained using the convolutional model. This model states that trace amplitudes have three controlling
More informationPASSIVE ACOUSTIC AND SEISMIC TOMOGRAPHY WITH OCEAN AMBIENT NOISE IN ORION
Proceedings of the International Conference Underwater Acoustic Measurements: Technologies &Results Heraklion, Crete, Greece, 28 th June 1 st July 2005 PASSIVE ACOUSTIC AND SEISMIC TOMOGRAPHY WITH OCEAN
More informationLow wavenumber reflectors
Low wavenumber reflectors Low wavenumber reflectors John C. Bancroft ABSTRACT A numerical modelling environment was created to accurately evaluate reflections from a D interface that has a smooth transition
More information7. Consider the following common offset gather collected with GPR.
Questions: GPR 1. Which of the following statements is incorrect when considering skin depth in GPR a. Skin depth is the distance at which the signal amplitude has decreased by a factor of 1/e b. Skin
More information28th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
SEISMIC SOURCE LOCATIONS AND PARAMETERS FOR SPARSE NETWORKS BY MATCHING OBSERVED SEISMOGRAMS TO SEMI-EMPIRICAL SYNTHETIC SEISMOGRAMS: IMPROVEMENTS TO THE PHASE SPECTRUM PARAMETERIZATION David. Salzberg
More informationAttenuation of high energy marine towed-streamer noise Nick Moldoveanu, WesternGeco
Nick Moldoveanu, WesternGeco Summary Marine seismic data have been traditionally contaminated by bulge waves propagating along the streamers that were generated by tugging and strumming from the vessel,
More informationInfluence of Peak Factors on Random Vibration Theory Based Site Response Analysis
6 th International Conference on Earthquake Geotechnical Engineering 1-4 November 2015 Christchurch, New Zealand Influence of Peak Factors on Random Vibration Theory Based Site Response Analysis X. Wang
More informationInterferometric Approach to Complete Refraction Statics Solution
Interferometric Approach to Complete Refraction Statics Solution Valentina Khatchatrian, WesternGeco, Calgary, Alberta, Canada VKhatchatrian@slb.com and Mike Galbraith, WesternGeco, Calgary, Alberta, Canada
More informationVariable-depth streamer acquisition: broadband data for imaging and inversion
P-246 Variable-depth streamer acquisition: broadband data for imaging and inversion Robert Soubaras, Yves Lafet and Carl Notfors*, CGGVeritas Summary This paper revisits the problem of receiver deghosting,
More informationThe case for longer sweeps in vibrator acquisition Malcolm Lansley, Sercel, John Gibson, Forest Lin, Alexandre Egreteau and Julien Meunier, CGGVeritas
The case for longer sweeps in vibrator acquisition Malcolm Lansley, Sercel, John Gibson, Forest Lin, Alexandre Egreteau and Julien Meunier, CGGVeritas There is growing interest in the oil and gas industry
More informationDirect Imaging of Group Velocity Dispersion Curves in Shallow Water Christopher Liner*, University of Houston; Lee Bell and Richard Verm, Geokinetics
Direct Imaging of Group Velocity Dispersion Curves in Shallow Water Christopher Liner*, University of Houston; Lee Bell and Richard Verm, Geokinetics Summary Geometric dispersion is commonly observed in
More informationTh P6 01 Retrieval of the P- and S-velocity Structure of the Groningen Gas Reservoir Using Noise Interferometry
Th P6 1 Retrieval of the P- and S-velocity Structure of the Groningen Gas Reservoir Using Noise Interferometry W. Zhou* (Utrecht University), H. Paulssen (Utrecht University) Summary The Groningen gas
More informationAnisotropic Frequency-Dependent Spreading of Seismic Waves from VSP Data Analysis
Anisotropic Frequency-Dependent Spreading of Seismic Waves from VSP Data Analysis Amin Baharvand Ahmadi* and Igor Morozov, University of Saskatchewan, Saskatoon, Saskatchewan amin.baharvand@usask.ca Summary
More informationThere is growing interest in the oil and gas industry to
Coordinated by JEFF DEERE JOHN GIBSON, FOREST LIN, ALEXANDRE EGRETEAU, and JULIEN MEUNIER, CGGVeritas MALCOLM LANSLEY, Sercel There is growing interest in the oil and gas industry to improve the quality
More information3/15/2010. Distance Distance along the ground (km) Time, (sec)
GG45 March 16, 21 Introduction to Seismic Exploration and Elementary Digital Analysis Some of the material I will cover today can be found in the book on pages 19-2 and 122-13. 13. However, much of what
More informationDetection of Obscured Targets
Detection of Obscured Targets Waymond R. Scott, Jr. and James Mcclellan School of Electrical and Computer Engineering Georgia Institute of Technology Atlanta, GA 30332-0250 waymond.scott@ece.gatech.edu
More informationDETERMINATION OF 1-D SHEAR WAVE VELOCITIES USING THE REFRACTION MICROTREMOR METHOD. Satish Pullammanappallil*, William Honjas*, and John N.
DETERMINATION OF 1-D SHEAR WAVE VELOCITIES USING THE REFRACTION MICROTREMOR METHOD Satish Pullammanappallil*, William Honjas*, and John N. Louie^ *Optim LLC, UNR-MS174, 1664 N. Virginia St. Reno, NV 89557;satish@optimsoftware.com
More informationABSTRACT INTRODUCTION. different curvatures at different times (see figure 1a and 1b).
APERTURE WIDTH SELECTION CRITERION IN KIRCHHOFF MIGRATION Richa Rastogi, Sudhakar Yerneni and Suhas Phadke Center for Development of Advanced Computing, Pune University Campus, Ganesh Khind, Pune 411007,
More informationSite-specific seismic hazard analysis
Site-specific seismic hazard analysis ABSTRACT : R.K. McGuire 1 and G.R. Toro 2 1 President, Risk Engineering, Inc, Boulder, Colorado, USA 2 Vice-President, Risk Engineering, Inc, Acton, Massachusetts,
More informationCoda Waveform Correlations
Chapter 5 Coda Waveform Correlations 5.1 Cross-Correlation of Seismic Coda 5.1.1 Introduction In the previous section, the generation of the surface wave component of the Green s function by the correlation
More informationERTH3021 Note: Terminology of Seismic Records
ERTH3021 Note: Terminology of Seismic Records This note is intended to assist in understanding of terminology used in practical exercises on 2D and 3D seismic acquisition geometries. A fundamental distinction
More informationSeismic processing workflow for supressing coherent noise while retaining low-frequency signal
Seismic processing for coherent noise suppression Seismic processing workflow for supressing coherent noise while retaining low-frequency signal Patricia E. Gavotti and Don C. Lawton ABSTRACT Two different
More informationEstimation of the Earth s Impulse Response: Deconvolution and Beyond. Gary Pavlis Indiana University Rick Aster New Mexico Tech
Estimation of the Earth s Impulse Response: Deconvolution and Beyond Gary Pavlis Indiana University Rick Aster New Mexico Tech Presentation for Imaging Science Workshop Washington University, November
More informationDesign of Geophysical Surveys in Transportation
Boise State University ScholarWorks CGISS Publications and Presentations Center for Geophysical Investigation of the Shallow Subsurface (CGISS) 1-1-2004 Design of Geophysical Surveys in Transportation
More informationAir-noise reduction on geophone data using microphone records
Air-noise reduction on geophone data using microphone records Air-noise reduction on geophone data using microphone records Robert R. Stewart ABSTRACT This paper proposes using microphone recordings of
More information2D field data applications
Chapter 5 2D field data applications In chapter 4, using synthetic examples, I showed how the regularized joint datadomain and image-domain inversion methods developed in chapter 3 overcome different time-lapse
More informationSeismic interference noise attenuation based on sparse inversion Zhigang Zhang* and Ping Wang (CGG)
Seismic interference noise attenuation based on sparse inversion Zhigang Zhang* and Ping Wang (CGG) Summary In marine seismic acquisition, seismic interference (SI) remains a considerable problem when
More informationDownloaded 09/04/18 to Redistribution subject to SEG license or copyright; see Terms of Use at
Processing of data with continuous source and receiver side wavefields - Real data examples Tilman Klüver* (PGS), Stian Hegna (PGS), and Jostein Lima (PGS) Summary In this paper, we describe the processing
More informationEffect of Frequency and Migration Aperture on Seismic Diffraction Imaging
IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Effect of Frequency and Migration Aperture on Seismic Diffraction Imaging To cite this article: Y. Bashir et al 2016 IOP Conf. Ser.:
More informationEstimation of a time-varying sea-surface profile for receiver-side de-ghosting Rob Telling* and Sergio Grion Shearwater Geoservices, UK
for receiver-side de-ghosting Rob Telling* and Sergio Grion Shearwater Geoservices, UK Summary The presence of a rough sea-surface during acquisition of marine seismic data leads to time- and space-dependent
More informationBorehole vibration response to hydraulic fracture pressure
Borehole vibration response to hydraulic fracture pressure Andy St-Onge* 1a, David W. Eaton 1b, and Adam Pidlisecky 1c 1 Department of Geoscience, University of Calgary, 2500 University Drive NW Calgary,
More information4.4 The transient electromagnetic method (TEM)
4.4 The transient electromagnetic method (TEM) 4.4 The transient electromagnetic method (TEM) 4.4.1 Basic principles and measuring techniques in TEM By the transient electromagnetic method, TEM, the electrical
More informationTh ELI1 08 Efficient Land Seismic Acquisition Sampling Using Rotational Data
Th ELI1 8 Efficient Land Seismic Acquisition Sampling Using Rotational Data P. Edme* (Schlumberger Gould Research), E. Muyzert (Sclumberger Gould Research) & E. Kragh (Schlumberger Gould Research) SUMMARY
More informationMASWaves User manual
MASWaves User manual Version 1 July 11, 2017 Preface/disclaimers... 2 References... 2 Acknowledgements... 2 1. Introduction... 3 2. Quick start guide... 5 3. MASWaves Dispersion... 11 3.1 Read data (MASWaves_read_data)...
More information=, (1) Summary. Theory. Introduction
Noise suppression for detection and location of microseismic events using a matched filter Leo Eisner*, David Abbott, William B. Barker, James Lakings and Michael P. Thornton, Microseismic Inc. Summary
More informationAnalysis of PS-to-PP amplitude ratios for seismic reflector characterisation: method and application
Analysis of PS-to-PP amplitude ratios for seismic reflector characterisation: method and application N. Maercklin, A. Zollo RISSC, Italy Abstract: Elastic parameters derived from seismic reflection data
More informationHIGH FREQUENCY MASW FOR NON-DESTRUCTIVE TESTING OF PAVEMENTS ACCELEROMETER APPROACH. Abstract
HIGH FREQUENCY MASW FOR NON-DESTRUCTIVE TESTING OF PAVEMENTS ACCELEROMETER APPROACH Downloaded 07/03/14 to 129.237.143.20. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/
More informationTHE ELECTROMAGNETIC FIELD THEORY. Dr. A. Bhattacharya
1 THE ELECTROMAGNETIC FIELD THEORY Dr. A. Bhattacharya The Underlying EM Fields The development of radar as an imaging modality has been based on power and power density It is important to understand some
More informationEnhanced subsurface response for marine CSEM surveying Frank A. Maaø* and Anh Kiet Nguyen, EMGS ASA
rank A. Maaø* and Anh Kiet Nguyen, EMGS ASA Summary A new robust method for enhancing marine CSEM subsurface response is presented. The method is demonstrated to enhance resolution and depth penetration
More informationShort Notes Characterization of a Continuous, Very Narrowband Seismic Signal near 2.08 Hz
Bulletin of the Seismological Society of America, 91, 6, pp. 1910 1916, December 2001 Short Notes Characterization of a Continuous, Very Narrowband Seismic Signal near 2.08 Hz by Kelly H. Liu and Stephen
More informationEFFECTS OF RAYLEIGH AND LOVE WAVES ON MICROTREMOR H/V SPECTRA
2232/4/A EFFECTS OF RAYLEIGH AND LOVE WAVES ON MICROTREMOR H/V SPECTRA Hiroshi ARAI 1 and Kohji TOKIMATSU 2 SUMMARY In order to simulate the horizontal-to-vertical (H/V) spectral ratios of microtremors,
More informationAN APPRAISAL OF SURFACE WAVE METHODS FOR SOIL CHARACTERIZATION
AN APPRAISAL OF SURFACE WAVE METHODS FOR SOIL CHARACTERIZATION By KHIEM TAT TRAN A THESIS PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE
More informationSummary. Volumetric Q tomography on offshore Brunei dataset
Success of high-resolution volumetric Q-tomography in the automatic detection of gas anomalies on offshore Brunei data Fatiha Gamar, Diego Carotti *, Patrice Guillaume, Amor Gacha, Laurent Lopes (CGG)
More informationResolution and location uncertainties in surface microseismic monitoring
Resolution and location uncertainties in surface microseismic monitoring Michael Thornton*, MicroSeismic Inc., Houston,Texas mthornton@microseismic.com Summary While related concepts, resolution and uncertainty
More informationSurvey results obtained in a complex geological environment with Midwater Stationary Cable Luc Haumonté*, Kietta; Weizhong Wang, Geotomo
Survey results obtained in a complex geological environment with Midwater Stationary Cable Luc Haumonté*, Kietta; Weizhong Wang, Geotomo Summary A survey with a novel acquisition technique was acquired
More informationLow Frequency Bottom Reflectivity from Reflection
Low Frequency Bottom Reflectivity from Reflection,Alexander Kritski 1 and Chris Jenkins 2 1 School of Geosciences, University of Sydney, NSW, 2 Ocean Sciences Institute, University of Sydney, NSW. Abstract
More informationGeophysical Journal International. Empirically determined finite frequency sensitivity kernels for surface waves
Empirically determined finite frequency sensitivity kernels for surface waves Journal: Manuscript ID: Draft Manuscript Type: Research Paper Date Submitted by the Author: Complete List of Authors: Lin,
More informationA robust x-t domain deghosting method for various source/receiver configurations Yilmaz, O., and Baysal, E., Paradigm Geophysical
A robust x-t domain deghosting method for various source/receiver configurations Yilmaz, O., and Baysal, E., Paradigm Geophysical Summary Here we present a method of robust seismic data deghosting for
More informationResponse spectrum Time history Power Spectral Density, PSD
A description is given of one way to implement an earthquake test where the test severities are specified by time histories. The test is done by using a biaxial computer aided servohydraulic test rig.
More informationGROUND_ROLL ATTENUATION IN THE RADIAL TRACE DOMAIN
GROUND_ROLL ATTENUATION IN THE RADIAL TRACE DOMAIN Bagheri, M. -, Dr.Riahi, M.A. -, Khaxar, Z.O. -, Hosseini, M -, Mohseni.D, R -. Adress: - Institute of Geophysics, University of Tehran Kargar Shomali
More informationUsing long sweep in land vibroseis acquisition
Using long sweep in land vibroseis acquisition Authors: Alexandre Egreteau, John Gibson, Forest Lin and Julien Meunier (CGGVeritas) Main objectives: Promote the use of long sweeps to compensate for the
More informationTu SRS3 07 Ultra-low Frequency Phase Assessment for Broadband Data
Tu SRS3 07 Ultra-low Frequency Phase Assessment for Broadband Data F. Yang* (CGG), R. Sablon (CGG) & R. Soubaras (CGG) SUMMARY Reliable low frequency content and phase alignment are critical for broadband
More informationGCM mapping Gedved - HydroGeophysics Group - Aarhus University
GCM mapping Gedved - HydroGeophysics Group - Aarhus University GCM mapping Gedved Report number 23-06-2017, June 2017 1. INDHOLDSFORTEGNELSE 1. Indholdsfortegnelse... 1 2. Project information... 2 3. DUALEM-421s...
More informationA k-mean characteristic function to improve STA/LTA detection
A k-mean characteristic function to improve STA/LTA detection Jubran Akram*,1, Daniel Peter 1, and David Eaton 2 1 King Abdullah University of Science and Technology (KAUST), Saudi Arabia 2 University
More informationG003 Data Preprocessing and Starting Model Preparation for 3D Inversion of Marine CSEM Surveys
G003 Data Preprocessing and Starting Model Preparation for 3D Inversion of Marine CSEM Surveys J.J. Zach* (EMGS ASA), F. Roth (EMGS ASA) & H. Yuan (EMGS Americas) SUMMARY The marine controlled-source electromagnetic
More informationCDP noise attenuation using local linear models
CDP noise attenuation CDP noise attenuation using local linear models Todor I. Todorov and Gary F. Margrave ABSTRACT Seismic noise attenuation plays an important part in a seismic processing flow. Spatial
More informationGround Penetrating Radar (day 1) EOSC Slide 1
Ground Penetrating Radar (day 1) Slide 1 Introduction to GPR Today s Topics Setup: Motivational Problems Physical Properties - Dielectric Permittivity and Radiowaves - Microwave Example Basic Principles:
More informationMEASUREMENT OF RAYLEIGH WAVE ATTENUATION IN GRANITE USING
MEASUREMENT OF RAYLEIGH WAVE ATTENUATION IN GRANITE USING LASER ULTRASONICS Joseph O. Owino and Laurence J. Jacobs School of Civil and Environmental Engineering Georgia Institute of Technology Atlanta
More informationGEOPIC, Oil & Natural Gas Corporation Ltd, Dehradun ,India b
Estimation of Seismic Q Using a Non-Linear (Gauss-Newton) Regression Parul Pandit * a, Dinesh Kumar b, T. R. Muralimohan a, Kunal Niyogi a,s.k. Das a a GEOPIC, Oil & Natural Gas Corporation Ltd, Dehradun
More informationRadial trace filtering revisited: current practice and enhancements
Radial trace filtering revisited: current practice and enhancements David C. Henley Radial traces revisited ABSTRACT Filtering seismic data in the radial trace (R-T) domain is an effective technique for
More informationTu A D Broadband Towed-Streamer Assessment, West Africa Deep Water Case Study
Tu A15 09 4D Broadband Towed-Streamer Assessment, West Africa Deep Water Case Study D. Lecerf* (PGS), D. Raistrick (PGS), B. Caselitz (PGS), M. Wingham (BP), J. Bradley (BP), B. Moseley (formaly BP) Summary
More informationSite Response from Incident Pnl Waves
Bulletin of the Seismological Society of America, Vol. 94, No. 1, pp. 357 362, February 2004 Site Response from Incident Pnl Waves by Brian Savage and Don V. Helmberger Abstract We developed a new method
More informationEstimating the epicenters of local and regional seismic sources, using the circle and chord method (Tutorial with exercise by hand and movies)
Topic Estimating the epicenters of local and regional seismic sources, using the circle and chord method (Tutorial with exercise by hand and movies) Author Version Peter Bormann (formerly GFZ German Research
More informationSpectral Analysis of Surface Waves for Damage Detection in Layered Concrete Structures
Paper Number 138, Proceedings of ACOUSTICS 2011 Spectral Analysis of Surface Waves for Damage Detection in Layered Concrete Structures ABSTRACT Mitchell Willcocks (1), Martin Veidt (1) and Greg Palmer
More informationMicrotremor measurements interpretations at the BVG station of the Italian accelerometric network
Microtremor measurements interpretations at the BVG station of the Italian accelerometric network Puglia R. 1, Tokeshi K. 2, Picozzi M. 3, d Alema E. 4, Parolai S. 3, Foti S. 2 1. Istituto Nazionale di
More informationOverview ta3520 Introduction to seismics
Overview ta3520 Introduction to seismics Fourier Analysis Basic principles of the Seismic Method Interpretation of Raw Seismic Records Seismic Instrumentation Processing of Seismic Reflection Data Vertical
More information2012 SEG SEG Las Vegas 2012 Annual Meeting Page 1
Full-wavefield, towed-marine seismic acquisition and applications David Halliday, Schlumberger Cambridge Research, Johan O. A. Robertsson, ETH Zürich, Ivan Vasconcelos, Schlumberger Cambridge Research,
More information41ST ANNUAL MARINE / OFFSHORE INDUSTRY OUTLOOK CONFERENCE GEOSCIENCE UPDATE: MR. ALAN M. BREED EVP SALES & MARKETING
41ST ANNUAL MARINE / OFFSHORE INDUSTRY OUTLOOK CONFERENCE GEOSCIENCE UPDATE: MR. ALAN M. BREED EVP SALES & MARKETING AGENDA GEOSCIENCE UPDATE Market Outlook Financial Overview Regulatory Issues Technology
More informationMicrophone Array Measurements for High-speed Train
Microphone Array Measurements for High-speed Train Korea Research Institute of Standards and Science Hyu-Sang Kwon 2016. 05. 31 2 Contents Railway Noise Sound Images Flow Noise Railway Noise Measurement
More information