Predicting Site Response
Predicting Site Response
• Based on theoretical calculations
– 1-D equivalent linear, non-linear
– 2-D and 3-D non-linear
• Needs geotechnical site properties
Imaging of Near-Surface Seismic
Slowness (Velocity) and Dampingand Damping
Ratios (Q)Ratios (Q)
• Sβ(z)(shear-wave slowness) (=1/velocity)
• Sα(z)(compressional-wave slowness)
• ξβ(z) (shear-wave damping ratio [Qβ])
Image What?Image What?
Why?Why?
• Site amplification
• Site classification for building codes
• Identification of liquefaction and landslide potential
• Correlation of various properties (e.g., geologic units and Vs)
Why Slowness?
• Travel time in layers directly proportional to slowness; travel
time fundamental in site response (e.g., T = 4*s*h = 4*travel
time)
• Can average slowness from several profiles depth-by-depth
• Slowness is the usual regression coefficient in fits of travel
time vs. depth
• Visual comparisons of slowness profiles more
meaningful for site response than velocity
profiles
Why Show Slowness Rather Than Velocity?
Large apparent differences in velocity in deeper layers (usually
higher velocity) become less important in plots of slowness
Focus attention on what contributes most to travel time in the layers
0 2 4 6
0
20
40
60
80
100
Shear-Wave Slowness (sec/km)
Depth(m)
Garner Valley
SASW Testing
Downhole Seismic
File:C:esg2006papergarner_valley_velocity_slowness_4ppt.draw;Date:2006-08-19;Time:09:00:59
0 500 1000 1500
0
20
40
60
80
100
Shear-Wave Velocity (m/sec)
Depth(m)
Imaging Slowness
• Invasive Methods
– Active sources
– Passive sources
• Noninvasive Methods
– Active sources
– Passive sources
Invasive Methods
• Active Sources
– surface source
– downhole source
• Passive sources
– Recordings of earthquake
waves in boreholes---not
covered in this talk
Invasive Method
Surface Source--
Downhole Receiver
(ssdhr)
(receiver can be on SCPT
rod)
One receiver moved up or
down hole
SURFACE SOURCE ---SUBSURFACE RECEIVERS
• downhole profiling
– velocities from surface
– data gaps filled by average velocity
– expensive (requires hole)
– depth range limited (but good to > 250 m)
• seismic cone penetrometer
– advantages of downhole
– inexpensive
– limited range
– not good for cobbly materials, rock
00.2
0 5 10 15 20 25 30 35 40 45 50
TravelTime(sec)
0.4
velTime(sec)
ile:C:coycreekgibbsjimCOYC_f_r_0_100_sideways_4ppt.draw;Date:2006-08-16;Time:17:27:18
Plotting sideways makes it
easier to see slopes changes
by viewing obliquely (an
exploration geophysics
trick)
Create a record
section—opposite
directions of
surface source
(red, blue traces)
Pick arrivals (black)
CCOC
0 50 100 150 200 250
0
0.2
0.4
0.6
TravelTime(sec)
sig = 1
sig = 2
sig = 3
sig = 4
sig = 5
model
CCOC -- 18 layers
-0.002
0
0.002
0.004
Residuial(sec)
ile:C:coycreekgibbsCoys_detail3_tt_resids_4ppt.draw;Date:2006-08-17;Time:08:26:11
Finer layering
in upper 100m
0 1 2 3 4 5 6
0
50
100
150
200
250
Slowness (sec/km)
Depth(m)
CCOC: S-Wave Slowness
Gibbs (Vs(30) = 232 m/s)
More detail (Vs(30) = 235 m/s)
File:C:coycreekgibbsgibbs_detail3_slowness_300m_4ppt.draw;Date:2006-08-19;Time:09:26:14
Two models
from the same
travel time
picks.
0.1 1 10
0.1
0.2
1
2
10
Frequency (Hz)
Amplification
8 layers
18 layers
vertical incidence,
density=2 gm/cc, Q
= 25, and a
halfspace with
V=1200 m/s and
density = 2.4 gm/cc
at 234 m depth.
File:C:coycreekgibbsnrattle_amps_gibbs_few_more_layers.draw;Date:2006-08-17;Time:08:34:45
The increased
resolution
makes little
difference in
site
amplification
SUBSURFACE SOURCE --- SUBSURFACE RECEIVERS
• crosshole
– “point” measurements in depth
– expensive (2 holes)
– velocity not appropriate for site response
• suspension logger
– rapid collection of data (no casing required)
– average velocity over small depth ranges
– can be used in deep holes
– expensive (requires borehole)
– no way of interpolating across data gaps
Cable Head
Head Reducer
Upper Geophone
Lower Geophone
Filter Tube
Source
Source Driver
Weight
Winch
7-Conductor cable
Diskette
with Data
OYO PS-160
Logger/Recorder
Overall Length ~ 25 ft
From Geovision
Downhole source--- P-S suspensionDownhole source--- P-S suspension logging (aka “PS Log”)logging (aka “PS Log”)
Dominant
frequency =
1000 Hz
Example from
Coyote Creek:
note 1) overall
trend; 2)
“scatter”; 3)
results
averaged over
various depth
intervals
reduces
“noise”
0 2 4 6 8
0
50
100
150
200
250
300
slowness (sec/km)
Depth(m)
CCOC: (Steller)
suspension log values
average of slowness over 5 m intervals
average of slowness over 10 m intervals
File:C:coycreekstellersuscoysx_slowness_300m_4ppt.draw;Date:2006-08-30;Time:02:16:10
“Noise”
fluctuations in
both S and P
logs agree with
variations in
lithology! (No
averaging)
Some Strengths of Invasive Methods
• Direct measure of velocity
• Surface source produces a model from the surface,
with depth intervals of poor or missing data replaced
by average layer (good for site amplification
calculations)
• PS suspension logging rapid, can be done soon after
hole drilled, no casing required, not limited in depth
range
Some Weaknesses of Invasive Methods
• Expensive! (If need to drill hole)
• Surface source may have difficulties in deep holes,
requires cased holes, logging must wait
• PS suspension log does not produce model from the
surface (but generally gets to within 1 to 2 m), and
there is no way of interpolating across depth
intervals with missing data.
Noninvasive Methods
• Active Sources
– e.g., SASW and MASW
• Passive sources (usually
microtremors)
– Single station
– Arrays (e.g., fk, SPAC)
• Combined active—passive sources
Overview of SASW and MASW
Method
• Spectral-Analysis-of-Surface-Waves
(SASW—2 receivers); Multichannel
Analysis of Surface Waves (MASW—
multiple receivers)
• Noninvasive and Nondestructive
• Based on Dispersive Characteristics of
Rayleigh Waves in a Layered Medium
SASW Field Procedure
• Transient or
Continuous Sources
(use several per
site)
• Receiver Geometry
Considerations:
– Near Field Effects
– Attenuation
– Expanding
Receiver Spread
– Lateral Variability
(Brown)
SASW & MASW Data Interpretation
80
60
40
20
0
Depth,m
8006004002000
Shear Wave Velocity, VS, m/s
Rinaldi Receiving Station
1
10
100
Wavelength,λRm
6004002000
Surface Wave Velocity, VR, m/s
Experimental Data
Theoretical Dispersion Curve
Rinaldi Receiving Station
(Brown)
Dispersion curve built from a number of subsets (different
source, different receiver spreads)
Some Factors That Influence
Accuracy of SASW & MASW Testing
• Lateral Variability of Subsurface
• Shear-Wave Velocity Gradient and
Contrasts
• Values of Poisson’s Ratio Assumed
in the inversion of the dispersion
curves
• Background Information on Site
Geology Improves the Models
Noninvasive Methods
• Passive sources (usually
microtremors)
– Single station (much work has been
done on this method---e.g., SESAME
project. I only mention it in passing,
using some slides from an ancientancient
paper)
(Boore & Toksöz, 1969)
Ellipticity (H/V) as a function of frequency depends on earth structure
Noninvasive Methods
• Passive sources (usually
microtremors)
– Multiple stations (usually two-
dimensional arrays)
(Hartzell, 2005)
The array of stations at WSP used by Hartzell
(Hartzell, 2005)
Inverting to obtain velocity profile
Noninvasive Methods
• Often active sources are limited in
depth (hard to generate low-
frequency motions)
• Station spacing used in passive
source experiments often too large
for resolution of near-surface
slowness
• Solution: Combined active—passive
sources
(Yoon and Rix, 2005)
An example
from the
CCOC—WSP
experiment
(active: f > 4
Hz; passive:
f<8 Hz)
Comparing Different Imaging Results at the
Same Site
• Direct comparison of slowness profiles
• Site amplification
– From empirical prediction equations
– Theoretical
• Full resonance
• Simplified (Square-root impedance)
Comparison of
slowness profiles:
0 2 4 6
0
20
40
60
80
100
Shear-Wave Slowness (sec/km)
Depth(m)
Garner Valley
PS Log A
PS Log B
SASW Testing
Downhole Seismic
File:C:esg2006papergarner_valley_slowness_4ppt.draw;Date:2006-08-22;Time:15:54:42
Coyote Creek Blind Interpretation Experiment (Asten and
Boore, 2005)
CCOC = Coyote Creek
Outdoor Classroom
The Experiment
• Measurements and interpretations done voluntarily
by many groups
• Interpretations “blind” to other results
• Interpretations sent to D. Boore
• Workshop held in May, 2004 to compare results
• Open-File report published in 2005 (containing a
summary by Asten & Boore and individual reports
from participants)
0 2 4 6 8 10
0
20
40
60
80
100
Slowness (sec/km)
Depth(m)
Shear Wave
Reference model
Reflection (Williams)
SASW (Bay, forward)
SASW (Stokoe, avg lb, ub)
SASW (Kayen, Wave-Eq)
MASW (Stephenson)
WSP: Active Sources
File:C:coycreekpaperwsp_active_s_deep_shallow.draw;Date:2006-08-19;Time:11:45:50
0 2 4 6 8 10
0
10
20
30
40
Slowness (sec/km)
Shear Wave
Reference model
Reflection (Williams)
SASW (Bay, forward)
SASW (Stokoe, avg lb, ub)
SASW (Kayen, Wave-Eq)
MASW (Stephenson)
WSP: Active Sources
Active sources at WSP: note larger near-surface & smaller
deep slownesses than reference for most methods.
0 2 4 6 8 10
0
50
100
150
200
250
300
Slowness (sec/km)
Depth(m)
Shear Wave
Reference model
SPAC (Asten, pkdec2)
SPAC (Hartzell)
H/V (Lang, Oct04)
Remi (Stephenson, mar05)
Remi (Louie)
WSP: Passive Sources
File:C:coycreekpaperwsp_passive_s_deep_shallow.draw;Date:2006-08-19;Time:11:46:18
0 2 4 6 8 10
0
10
20
30
40
Slowness (sec/km)
Shear Wave
Reference model
SPAC (Asten, pkdec2)
SPAC (Hartzell)
H/V (Lang, Oct04)
Remi (Stephenson, mar05)
Remi (Louie)
WSP: Passive Sources
Passive sources at WSP: note larger near-surface & smaller
deep slownesses than reference for most methods. Models
extend to greater depth than do the models from active
sources
0 2 4 6 8 10
0
10
20
30
40
Slowness (sec/km)
Shear Wave
Reference model
MASW+MAM (Hayashi)
MASW+MAM (Rix)
WSP: Active + Passive Sources
File:C:coycreekpaperwsp_both_s_deep_shallow.draw;Date:2006-08-19;Time:11:47:00
0 2 4 6 8 10
0
50
100
150
200
Slowness (sec/km)
Depth(m)
Shear Wave
Reference model
MASW+MAM (Hayashi)
MASW+MAM (Rix)
WSP: Active + Passive Sources
Combined active & passive sources at WSP: note larger
near-surface slownesses than reference
0.01 0.1 1 10
0.8
0.9
1
1.1
Period (s)
Amplification,relativetotheV30fromtheCCOCboreholeaverage
Red: Active Sources; Blue: Passive & Combined Sources
SASW, CCOC (Stokoe, Cl1)
SASW, CCOC (Stokoe, Cl2 avg)
Reflection, WSP (Williams)
SASW, WSP (Kayen)
MASW, WSP (Stephenson)
SASW, WSP (Stokoe, avg)
MASW+MAM, WSP (Hayashi)
MASW+FK, WSP (Rix)
H/V, WSP (Lang, oct04)
SPAC, WSP (Asten, pkdec2)
SPAC, WSP (Hartzell)
ReMi (Stephenson, mar05)
ReMi (Louie)
File:C:coycreekpaperamps_using_v30.draw;Date:2006-08-18;Time:08:40:28
leading to these small differences in empirically-based
amplifications based on V30 (red=active; blue=passive &
combined)
Average slownesses tend to converge near 30 mconverge near 30 m (coincidence?) with
systematic differences shallower and deeper (both types of source give larger
shallow slowness; at 30 m the slowness from active sources is larger than the
reference and on average is smaller than the reference for passive sources.
0 2 4 6 8 10
1
2
10
20
100
200
Slowness (sec/km)
Depth(m)
Active Sources (CCOC & WSP)
reference model
SASW, CCOC (Bay)
SASW, CCOC (Stokoe, CL1)
SASW, CCOC (Stokoe, CL2 avg)
reflection, WSP (Williams)
SASW, WSP (Bay)
SASW, WSP (Kayen)
MASW, WSP (Stephenson)
SASW, WSP (Stokoe, avg)
0 2 4 6 8 10
1
2
10
20
100
200
Slowness (sec/km)
Passive & Combined Sources (WSP)
reference model
SPAC (Asten, pkdec2)
H/V (Lang, oct04)
SPAC (Hartzell)
ReMi (Stephenson, mar05)
ReMi (Louie)
MASW+MAM, WSP (Hayashi)
MASW+FK, WSP (Rix)
File:C:coycreekpaperccoc_wsp_slowness_active_passive.draw;Date:2006-08-23;Time:09:24:19
1 2 10 20
2
3
4
Frequency (Hz)
Amplification(relativeto1500m/s;nodamping)
Active Sources
reference model
SASW, CCOC (Bay)
SASW, CCOC (Stokoe, CL1)
SASW, CCOC (Stokoe, CL2 avg)
Reflection, WSP (Williams)
SASW, WSP (Bay)
SASW, WSP (Kayen)
MASW, WSP (Stephenson)
SASW, WSP (Stokoe, avg)
reference model with damping ( =0.04s)
Kayen, with damping ( =0.04s)
1 2 10 20
2
3
4
Frequency (Hz)
Passive & Combined Sources (WSP)
reference model
SPAC, WSP (Asten, pkdec2)
H/V, WSP (Lang, oct04)
SPAC, WSP (Hartzell)
ReMi, WSP (Stephenson_mar05)
ReMi, WSP (Louie)
MASW+MAM, WSP (Hayashi)
MASW+FK, WSP (Rix)
File:C:coycreekpaperccoc_wsp_amps_active_passive.draw;Date:2006-08-18;Time:08:43:32
But larger differenceslarger differences at higher frequenciesat higher frequencies (up to 40%)
(V30 corresponds to ~ 2 Hz)
Summary (short)
• Many methods available for imaging seismic
slowness
• Noninvasive methods work well, with some
suggestions of systematic departures from borehole
methods
• Several measures of site amplification show little
sensitivity to the differences in models (on the order
of factors of 1.4 or less)
• Site amplifications show trends with V30, but the
remaining scatter in observed ground motions is
large

More Related Content

PPTX
Seismic data processing
PDF
Fundementals of MASW
PDF
Bp sesmic interpretation
PPTX
Seismic Attributes .pptx
PPTX
Seismic refraction method lecture 21
PPTX
Avo ppt (Amplitude Variation with Offset)
PDF
Seismic sources
PPTX
Lecture 23 april29 static correction
Seismic data processing
Fundementals of MASW
Bp sesmic interpretation
Seismic Attributes .pptx
Seismic refraction method lecture 21
Avo ppt (Amplitude Variation with Offset)
Seismic sources
Lecture 23 april29 static correction

What's hot (20)

PPT
ÖNCEL AKADEMİ: INTRODUCTION TO GEOPHYSICS
PPTX
Seismic interpretation - Fluvial Deltaic System
DOCX
Filtering in seismic data processing? How filtering help to suppress noises.
DOCX
Quantitative and Qualitative Seismic Interpretation of Seismic Data
PPT
Seismic Data Processing, Ahmed Osama
PDF
Basics of seismic interpretation
PPTX
Role of Seismic Attributes in Petroleum Exploration_30May22.pptx
PDF
seismic inversion
PPTX
Introduction to Seismic Method
PDF
Reservoir Geophysics
PPTX
Seismic Attributes
PPT
Fundamentals of Seismic Refraction
PPTX
Reservoir Geophysics : Brian Russell Lecture 1
PDF
Insitu stresses and measurement
PDF
Simple seismic processing workflow
PPT
Refraction Seismic
PDF
Principles of seismic data processing m.m.badawy
PPT
Seismic acquisition
PPTX
PROCESSING PHASE-IV.pptx
PDF
Direct hydrocarbon indicators (DHI)
ÖNCEL AKADEMİ: INTRODUCTION TO GEOPHYSICS
Seismic interpretation - Fluvial Deltaic System
Filtering in seismic data processing? How filtering help to suppress noises.
Quantitative and Qualitative Seismic Interpretation of Seismic Data
Seismic Data Processing, Ahmed Osama
Basics of seismic interpretation
Role of Seismic Attributes in Petroleum Exploration_30May22.pptx
seismic inversion
Introduction to Seismic Method
Reservoir Geophysics
Seismic Attributes
Fundamentals of Seismic Refraction
Reservoir Geophysics : Brian Russell Lecture 1
Insitu stresses and measurement
Simple seismic processing workflow
Refraction Seismic
Principles of seismic data processing m.m.badawy
Seismic acquisition
PROCESSING PHASE-IV.pptx
Direct hydrocarbon indicators (DHI)
Ad

Similar to Multichannel analysis of surface waves (20)

PDF
Vs30 measurements for Seismic Site Classification
PDF
Engineering Seismology
PDF
A Comparison Between Shear Wave Velocities
PDF
Crosshole Seismic Reflection: Coal Mine Fields
PPTX
Introduction to Seismic Method
PPT
S24 naldi
PPT
S24 naldi
PPT
S24 naldi
PPTX
Seismic Method Estimate velocity from seismic data.pptx
PPT
Fundamentals of engineering geophysics.ppt
PDF
2D MASW ANALYSIS FOR GEOTECHNICAL ENGINEERING
PDF
Lateral resolution and lithological interpretation of surface wave profi ling
PPT
Seismic Interpretation in geology.pptxxx
PPTX
SEISMIC METHOD
PDF
Seismic Methods
PDF
Engineering geophysical study of unconsolidated top soil using shallow seismi...
PDF
Seismic Traveltime Tomography Of The Crust And Lithosphere N Rawlinson
PDF
Introduction to seismic interpretation
PPTX
Subsurface Exploration Methods in Geotechnical Engineering
PPTX
New introduction to seismic method
Vs30 measurements for Seismic Site Classification
Engineering Seismology
A Comparison Between Shear Wave Velocities
Crosshole Seismic Reflection: Coal Mine Fields
Introduction to Seismic Method
S24 naldi
S24 naldi
S24 naldi
Seismic Method Estimate velocity from seismic data.pptx
Fundamentals of engineering geophysics.ppt
2D MASW ANALYSIS FOR GEOTECHNICAL ENGINEERING
Lateral resolution and lithological interpretation of surface wave profi ling
Seismic Interpretation in geology.pptxxx
SEISMIC METHOD
Seismic Methods
Engineering geophysical study of unconsolidated top soil using shallow seismi...
Seismic Traveltime Tomography Of The Crust And Lithosphere N Rawlinson
Introduction to seismic interpretation
Subsurface Exploration Methods in Geotechnical Engineering
New introduction to seismic method
Ad

Recently uploaded (20)

PPT
Comprehensive Java Training Deck - Advanced topics
DOCX
ENVIRONMENTAL PROTECTION AND MANAGEMENT (18CVL756)
PPTX
CS6006 - CLOUD COMPUTING - Module - 1.pptx
PDF
Beginners-Guide-to-Artificial-Intelligence.pdf
PDF
ECT443_instrumentation_Engg_mod-1.pdf indroduction to instrumentation
PDF
Research on ultrasonic sensor for TTU.pdf
PPTX
INTERNET OF THINGS - EMBEDDED SYSTEMS AND INTERNET OF THINGS
PDF
Mechanics of materials week 2 rajeshwari
PDF
Principles of operation, construction, theory, advantages and disadvantages, ...
PPTX
Wireless sensor networks (WSN) SRM unit 2
PPTX
DATA STRCUTURE LABORATORY -BCSL305(PRG1)
PPT
Programmable Logic Controller PLC and Industrial Automation
PDF
V2500 Owner and Operatore Guide for Airbus
PDF
Engineering Solutions for Ethical Dilemmas in Healthcare (www.kiu.ac.ug)
PDF
IAE-V2500 Engine Airbus Family A319/320
PDF
Software defined netwoks is useful to learn NFV and virtual Lans
PDF
AIGA 012_04 Cleaning of equipment for oxygen service_reformat Jan 12.pdf
PDF
IAE-V2500 Engine for Airbus Family 319/320
PDF
Project_Mgmt_Institute_-Marc Marc Marc .pdf
PDF
Lesson 3 .pdf
Comprehensive Java Training Deck - Advanced topics
ENVIRONMENTAL PROTECTION AND MANAGEMENT (18CVL756)
CS6006 - CLOUD COMPUTING - Module - 1.pptx
Beginners-Guide-to-Artificial-Intelligence.pdf
ECT443_instrumentation_Engg_mod-1.pdf indroduction to instrumentation
Research on ultrasonic sensor for TTU.pdf
INTERNET OF THINGS - EMBEDDED SYSTEMS AND INTERNET OF THINGS
Mechanics of materials week 2 rajeshwari
Principles of operation, construction, theory, advantages and disadvantages, ...
Wireless sensor networks (WSN) SRM unit 2
DATA STRCUTURE LABORATORY -BCSL305(PRG1)
Programmable Logic Controller PLC and Industrial Automation
V2500 Owner and Operatore Guide for Airbus
Engineering Solutions for Ethical Dilemmas in Healthcare (www.kiu.ac.ug)
IAE-V2500 Engine Airbus Family A319/320
Software defined netwoks is useful to learn NFV and virtual Lans
AIGA 012_04 Cleaning of equipment for oxygen service_reformat Jan 12.pdf
IAE-V2500 Engine for Airbus Family 319/320
Project_Mgmt_Institute_-Marc Marc Marc .pdf
Lesson 3 .pdf

Multichannel analysis of surface waves

  • 2. Predicting Site Response • Based on theoretical calculations – 1-D equivalent linear, non-linear – 2-D and 3-D non-linear • Needs geotechnical site properties
  • 3. Imaging of Near-Surface Seismic Slowness (Velocity) and Dampingand Damping Ratios (Q)Ratios (Q)
  • 4. • Sβ(z)(shear-wave slowness) (=1/velocity) • Sα(z)(compressional-wave slowness) • ξβ(z) (shear-wave damping ratio [Qβ]) Image What?Image What? Why?Why? • Site amplification • Site classification for building codes • Identification of liquefaction and landslide potential • Correlation of various properties (e.g., geologic units and Vs)
  • 5. Why Slowness? • Travel time in layers directly proportional to slowness; travel time fundamental in site response (e.g., T = 4*s*h = 4*travel time) • Can average slowness from several profiles depth-by-depth • Slowness is the usual regression coefficient in fits of travel time vs. depth • Visual comparisons of slowness profiles more meaningful for site response than velocity profiles
  • 6. Why Show Slowness Rather Than Velocity? Large apparent differences in velocity in deeper layers (usually higher velocity) become less important in plots of slowness Focus attention on what contributes most to travel time in the layers 0 2 4 6 0 20 40 60 80 100 Shear-Wave Slowness (sec/km) Depth(m) Garner Valley SASW Testing Downhole Seismic File:C:esg2006papergarner_valley_velocity_slowness_4ppt.draw;Date:2006-08-19;Time:09:00:59 0 500 1000 1500 0 20 40 60 80 100 Shear-Wave Velocity (m/sec) Depth(m)
  • 7. Imaging Slowness • Invasive Methods – Active sources – Passive sources • Noninvasive Methods – Active sources – Passive sources
  • 8. Invasive Methods • Active Sources – surface source – downhole source • Passive sources – Recordings of earthquake waves in boreholes---not covered in this talk
  • 9. Invasive Method Surface Source-- Downhole Receiver (ssdhr) (receiver can be on SCPT rod) One receiver moved up or down hole
  • 10. SURFACE SOURCE ---SUBSURFACE RECEIVERS • downhole profiling – velocities from surface – data gaps filled by average velocity – expensive (requires hole) – depth range limited (but good to > 250 m) • seismic cone penetrometer – advantages of downhole – inexpensive – limited range – not good for cobbly materials, rock
  • 11. 00.2 0 5 10 15 20 25 30 35 40 45 50 TravelTime(sec) 0.4 velTime(sec) ile:C:coycreekgibbsjimCOYC_f_r_0_100_sideways_4ppt.draw;Date:2006-08-16;Time:17:27:18 Plotting sideways makes it easier to see slopes changes by viewing obliquely (an exploration geophysics trick) Create a record section—opposite directions of surface source (red, blue traces) Pick arrivals (black) CCOC
  • 12. 0 50 100 150 200 250 0 0.2 0.4 0.6 TravelTime(sec) sig = 1 sig = 2 sig = 3 sig = 4 sig = 5 model CCOC -- 18 layers -0.002 0 0.002 0.004 Residuial(sec) ile:C:coycreekgibbsCoys_detail3_tt_resids_4ppt.draw;Date:2006-08-17;Time:08:26:11 Finer layering in upper 100m
  • 13. 0 1 2 3 4 5 6 0 50 100 150 200 250 Slowness (sec/km) Depth(m) CCOC: S-Wave Slowness Gibbs (Vs(30) = 232 m/s) More detail (Vs(30) = 235 m/s) File:C:coycreekgibbsgibbs_detail3_slowness_300m_4ppt.draw;Date:2006-08-19;Time:09:26:14 Two models from the same travel time picks.
  • 14. 0.1 1 10 0.1 0.2 1 2 10 Frequency (Hz) Amplification 8 layers 18 layers vertical incidence, density=2 gm/cc, Q = 25, and a halfspace with V=1200 m/s and density = 2.4 gm/cc at 234 m depth. File:C:coycreekgibbsnrattle_amps_gibbs_few_more_layers.draw;Date:2006-08-17;Time:08:34:45 The increased resolution makes little difference in site amplification
  • 15. SUBSURFACE SOURCE --- SUBSURFACE RECEIVERS • crosshole – “point” measurements in depth – expensive (2 holes) – velocity not appropriate for site response • suspension logger – rapid collection of data (no casing required) – average velocity over small depth ranges – can be used in deep holes – expensive (requires borehole) – no way of interpolating across data gaps
  • 16. Cable Head Head Reducer Upper Geophone Lower Geophone Filter Tube Source Source Driver Weight Winch 7-Conductor cable Diskette with Data OYO PS-160 Logger/Recorder Overall Length ~ 25 ft From Geovision Downhole source--- P-S suspensionDownhole source--- P-S suspension logging (aka “PS Log”)logging (aka “PS Log”) Dominant frequency = 1000 Hz
  • 17. Example from Coyote Creek: note 1) overall trend; 2) “scatter”; 3) results averaged over various depth intervals reduces “noise” 0 2 4 6 8 0 50 100 150 200 250 300 slowness (sec/km) Depth(m) CCOC: (Steller) suspension log values average of slowness over 5 m intervals average of slowness over 10 m intervals File:C:coycreekstellersuscoysx_slowness_300m_4ppt.draw;Date:2006-08-30;Time:02:16:10
  • 18. “Noise” fluctuations in both S and P logs agree with variations in lithology! (No averaging)
  • 19. Some Strengths of Invasive Methods • Direct measure of velocity • Surface source produces a model from the surface, with depth intervals of poor or missing data replaced by average layer (good for site amplification calculations) • PS suspension logging rapid, can be done soon after hole drilled, no casing required, not limited in depth range
  • 20. Some Weaknesses of Invasive Methods • Expensive! (If need to drill hole) • Surface source may have difficulties in deep holes, requires cased holes, logging must wait • PS suspension log does not produce model from the surface (but generally gets to within 1 to 2 m), and there is no way of interpolating across depth intervals with missing data.
  • 21. Noninvasive Methods • Active Sources – e.g., SASW and MASW • Passive sources (usually microtremors) – Single station – Arrays (e.g., fk, SPAC) • Combined active—passive sources
  • 22. Overview of SASW and MASW Method • Spectral-Analysis-of-Surface-Waves (SASW—2 receivers); Multichannel Analysis of Surface Waves (MASW— multiple receivers) • Noninvasive and Nondestructive • Based on Dispersive Characteristics of Rayleigh Waves in a Layered Medium
  • 23. SASW Field Procedure • Transient or Continuous Sources (use several per site) • Receiver Geometry Considerations: – Near Field Effects – Attenuation – Expanding Receiver Spread – Lateral Variability (Brown)
  • 24. SASW & MASW Data Interpretation 80 60 40 20 0 Depth,m 8006004002000 Shear Wave Velocity, VS, m/s Rinaldi Receiving Station 1 10 100 Wavelength,λRm 6004002000 Surface Wave Velocity, VR, m/s Experimental Data Theoretical Dispersion Curve Rinaldi Receiving Station (Brown) Dispersion curve built from a number of subsets (different source, different receiver spreads)
  • 25. Some Factors That Influence Accuracy of SASW & MASW Testing • Lateral Variability of Subsurface • Shear-Wave Velocity Gradient and Contrasts • Values of Poisson’s Ratio Assumed in the inversion of the dispersion curves • Background Information on Site Geology Improves the Models
  • 26. Noninvasive Methods • Passive sources (usually microtremors) – Single station (much work has been done on this method---e.g., SESAME project. I only mention it in passing, using some slides from an ancientancient paper)
  • 27. (Boore & Toksöz, 1969) Ellipticity (H/V) as a function of frequency depends on earth structure
  • 28. Noninvasive Methods • Passive sources (usually microtremors) – Multiple stations (usually two- dimensional arrays)
  • 29. (Hartzell, 2005) The array of stations at WSP used by Hartzell
  • 30. (Hartzell, 2005) Inverting to obtain velocity profile
  • 31. Noninvasive Methods • Often active sources are limited in depth (hard to generate low- frequency motions) • Station spacing used in passive source experiments often too large for resolution of near-surface slowness • Solution: Combined active—passive sources
  • 32. (Yoon and Rix, 2005) An example from the CCOC—WSP experiment (active: f > 4 Hz; passive: f<8 Hz)
  • 33. Comparing Different Imaging Results at the Same Site • Direct comparison of slowness profiles • Site amplification – From empirical prediction equations – Theoretical • Full resonance • Simplified (Square-root impedance)
  • 34. Comparison of slowness profiles: 0 2 4 6 0 20 40 60 80 100 Shear-Wave Slowness (sec/km) Depth(m) Garner Valley PS Log A PS Log B SASW Testing Downhole Seismic File:C:esg2006papergarner_valley_slowness_4ppt.draw;Date:2006-08-22;Time:15:54:42
  • 35. Coyote Creek Blind Interpretation Experiment (Asten and Boore, 2005) CCOC = Coyote Creek Outdoor Classroom
  • 36. The Experiment • Measurements and interpretations done voluntarily by many groups • Interpretations “blind” to other results • Interpretations sent to D. Boore • Workshop held in May, 2004 to compare results • Open-File report published in 2005 (containing a summary by Asten & Boore and individual reports from participants)
  • 37. 0 2 4 6 8 10 0 20 40 60 80 100 Slowness (sec/km) Depth(m) Shear Wave Reference model Reflection (Williams) SASW (Bay, forward) SASW (Stokoe, avg lb, ub) SASW (Kayen, Wave-Eq) MASW (Stephenson) WSP: Active Sources File:C:coycreekpaperwsp_active_s_deep_shallow.draw;Date:2006-08-19;Time:11:45:50 0 2 4 6 8 10 0 10 20 30 40 Slowness (sec/km) Shear Wave Reference model Reflection (Williams) SASW (Bay, forward) SASW (Stokoe, avg lb, ub) SASW (Kayen, Wave-Eq) MASW (Stephenson) WSP: Active Sources Active sources at WSP: note larger near-surface & smaller deep slownesses than reference for most methods.
  • 38. 0 2 4 6 8 10 0 50 100 150 200 250 300 Slowness (sec/km) Depth(m) Shear Wave Reference model SPAC (Asten, pkdec2) SPAC (Hartzell) H/V (Lang, Oct04) Remi (Stephenson, mar05) Remi (Louie) WSP: Passive Sources File:C:coycreekpaperwsp_passive_s_deep_shallow.draw;Date:2006-08-19;Time:11:46:18 0 2 4 6 8 10 0 10 20 30 40 Slowness (sec/km) Shear Wave Reference model SPAC (Asten, pkdec2) SPAC (Hartzell) H/V (Lang, Oct04) Remi (Stephenson, mar05) Remi (Louie) WSP: Passive Sources Passive sources at WSP: note larger near-surface & smaller deep slownesses than reference for most methods. Models extend to greater depth than do the models from active sources
  • 39. 0 2 4 6 8 10 0 10 20 30 40 Slowness (sec/km) Shear Wave Reference model MASW+MAM (Hayashi) MASW+MAM (Rix) WSP: Active + Passive Sources File:C:coycreekpaperwsp_both_s_deep_shallow.draw;Date:2006-08-19;Time:11:47:00 0 2 4 6 8 10 0 50 100 150 200 Slowness (sec/km) Depth(m) Shear Wave Reference model MASW+MAM (Hayashi) MASW+MAM (Rix) WSP: Active + Passive Sources Combined active & passive sources at WSP: note larger near-surface slownesses than reference
  • 40. 0.01 0.1 1 10 0.8 0.9 1 1.1 Period (s) Amplification,relativetotheV30fromtheCCOCboreholeaverage Red: Active Sources; Blue: Passive & Combined Sources SASW, CCOC (Stokoe, Cl1) SASW, CCOC (Stokoe, Cl2 avg) Reflection, WSP (Williams) SASW, WSP (Kayen) MASW, WSP (Stephenson) SASW, WSP (Stokoe, avg) MASW+MAM, WSP (Hayashi) MASW+FK, WSP (Rix) H/V, WSP (Lang, oct04) SPAC, WSP (Asten, pkdec2) SPAC, WSP (Hartzell) ReMi (Stephenson, mar05) ReMi (Louie) File:C:coycreekpaperamps_using_v30.draw;Date:2006-08-18;Time:08:40:28 leading to these small differences in empirically-based amplifications based on V30 (red=active; blue=passive & combined)
  • 41. Average slownesses tend to converge near 30 mconverge near 30 m (coincidence?) with systematic differences shallower and deeper (both types of source give larger shallow slowness; at 30 m the slowness from active sources is larger than the reference and on average is smaller than the reference for passive sources. 0 2 4 6 8 10 1 2 10 20 100 200 Slowness (sec/km) Depth(m) Active Sources (CCOC & WSP) reference model SASW, CCOC (Bay) SASW, CCOC (Stokoe, CL1) SASW, CCOC (Stokoe, CL2 avg) reflection, WSP (Williams) SASW, WSP (Bay) SASW, WSP (Kayen) MASW, WSP (Stephenson) SASW, WSP (Stokoe, avg) 0 2 4 6 8 10 1 2 10 20 100 200 Slowness (sec/km) Passive & Combined Sources (WSP) reference model SPAC (Asten, pkdec2) H/V (Lang, oct04) SPAC (Hartzell) ReMi (Stephenson, mar05) ReMi (Louie) MASW+MAM, WSP (Hayashi) MASW+FK, WSP (Rix) File:C:coycreekpaperccoc_wsp_slowness_active_passive.draw;Date:2006-08-23;Time:09:24:19
  • 42. 1 2 10 20 2 3 4 Frequency (Hz) Amplification(relativeto1500m/s;nodamping) Active Sources reference model SASW, CCOC (Bay) SASW, CCOC (Stokoe, CL1) SASW, CCOC (Stokoe, CL2 avg) Reflection, WSP (Williams) SASW, WSP (Bay) SASW, WSP (Kayen) MASW, WSP (Stephenson) SASW, WSP (Stokoe, avg) reference model with damping ( =0.04s) Kayen, with damping ( =0.04s) 1 2 10 20 2 3 4 Frequency (Hz) Passive & Combined Sources (WSP) reference model SPAC, WSP (Asten, pkdec2) H/V, WSP (Lang, oct04) SPAC, WSP (Hartzell) ReMi, WSP (Stephenson_mar05) ReMi, WSP (Louie) MASW+MAM, WSP (Hayashi) MASW+FK, WSP (Rix) File:C:coycreekpaperccoc_wsp_amps_active_passive.draw;Date:2006-08-18;Time:08:43:32 But larger differenceslarger differences at higher frequenciesat higher frequencies (up to 40%) (V30 corresponds to ~ 2 Hz)
  • 43. Summary (short) • Many methods available for imaging seismic slowness • Noninvasive methods work well, with some suggestions of systematic departures from borehole methods • Several measures of site amplification show little sensitivity to the differences in models (on the order of factors of 1.4 or less) • Site amplifications show trends with V30, but the remaining scatter in observed ground motions is large

Editor's Notes

  • #23: The Spectral analysis of surface waves method is the successor to the steady state Rayleigh wave method developed in the 1950’s. Much of the development of the modern SASW method was carried out at UT Austin in the early 1980’s. SASW testing is used to obtain a shear wave velocity profile It is non-invasive and non-destructive - testing is performed on the ground surface and strains are in the elastic range Instead of measuring shear wave velocity directly, Rayleigh wave velocities are measured and Vs is inferred.
  • #24: The general testing setup is shown here. A seismic source generates surface waves, which are monitored by two in-line receivers. Both transient and continuous dynamic sources are used to generate surface waves, with the data usually cleaner from continuous swept-sine sources. A vibroseis truck (slide) was used for the long wavelengths and various hand-held hammers (slide) were used for the short wavelengths. Several factors must be considered in receiver geometry. To avoid near field effects associated with Rayleigh waves and body waves, the distance from the source to the receiver, d1, is at least half of the maximum recorded wavelength. Attenuation reduces signal quality if d1 is greater than 4-10 wavelengths, depending on the source. Therefore, an expanding receiver spread is used, with overlap between the wavelengths recorded in each setup. To minimize lateral variability, forward and reverse profiles are taken, usually with a common centerline. The time records from the two geophones are transformed to the frequency domain to generate the dispersion curve. The most important data are the phase of the cross-power spectrum and the coherence. It is important that the frequency domain calculation be done in the field so that the experiment can be modified as needed.
  • #25: From the unwrapped phase of the cross power spectrum, the Rayleigh wave velocity is calculated, given the frequency and interreceiver distance. The dispersion curves from each receiver spacing are combined to generate the composite dispersion curve, which is representative of the site. Several theoretical solutions are used to model the dispersion curve: Fundamental mode Rayleigh waves only, and full stress wave solutions that incorporate higher modes of Rayleigh wave propagation, body wave energy, and receiver location. The full stress wave solution generally gives better results so the results from that model will be shown. The parameters in the layered earth model used to calculate the dispersion curve consist of layer thickness, shear wave velocity, Poisson’s ratio, and mass density. Usually only shear wave velocity and layer thickness are adjusted to match the dispersion curve, since they have the largest influence. The resolvable depth varies from about one half the longest wavelength to one fifth or less, depending on the site.
  • #26: Obviously, the results from SASW testing were closer to those from downhole testing at some sites. There are several possible reasons for this; SASW testing samples a much larger volume of material than downhole testing. The SASW results are averaged across several hundred meters. If the subsurface varies laterally or is non-homogeneous the material sampled in SASW and downhole seismic testing may be different. At sites where the shear wave velocity increases gradually, the SASW data are easiest to interpret and most accurate. Large shear-wave velocity gradients limit the resolvable depth, because the dispersion curve never levels out to a velocity representative of individual layers or a half space. SASW testing does not resolve layer boundaries as well as average properties. Interface resolution also decreases with depth. The value of Poisson’s ratio assumed in the model for Rayleigh wave dispersion has a greater effect than is often thought. A sensitivity study showed that in saturated sediments.