Geophysical Methods for Geotechnical Site Characterization

Short Course – 23rd of February 2014
Introduction to Geophysical Tests
VS (m/s)
Sebastiano Foti
(ITALY)
Email: [email protected]
www.soilmech.polito.it/people/foti_sebastiano
Geophysical methods for geotechnical site characterization
Outline
• Geophysical methods
– Scope and potential for geotechnical and
geoenvironmental characterization
– Use of seismic velocities
– Significance of other geophysical parameters
– In-hole vs surface methods
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Geophysical Methods
Geophysical methods are indirect surveying techniques based on
measurements carried out on the ground surface or in holes. They
allow the distribution of physical properties of the subsurface to be
estimated and correlated with engineering information.
They are based on the excitation of an object with an energy field (artificial
or natural) and on the measurement of the object response.
The interpretation of the object response allows the object to be
characterised.
\
?
Geocongress 2014 23rd February 2014
?
?
?
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Geophysical parameters
•
•
•
•
•
•
Density
Electrical Conductivity (or Resistivity)
Electrical Permittivity
Magnetic Suscettibility
Chargeability
Seismic velocities
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Geophysical parameters
•
•
•
•
•
•
Density
Electrical Conductivity (or Resistivity)
Electrical Permittivity
Magnetic Suscettibility
Chargeability
Seismic velocities
Direct relationship to mechanical
parameters of the medium
(Elastic Moduli)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Applicability of in situ tests
(Mayne et al, 2002)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Geotechnical and geoenvironmental site characterization
In the context of site characterization for engineering purposes,
the role of geophysical methods is twofold:
• evaluation of geometrical boundaries to model subsoil
conditions (e.g. stratigraphy but also physical inclusions or
hydrogeological features);
• evaluation of physical/mechanical parameters of direct use
for geotechnical modeling.
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Identification of stratigraphic sequence
Pugin et al., 2009
Seismic methods:
e.g. seismic reflection to identify an acquifer
In combination with conventional investigation:
e.g. boreholes logs allow calibration / identification of litography
geophysical surveys allow for 2D/3D extension
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Lateral variations (shallow faults)
e.g. seismic methods: surface wave tests
Geological model
(expected)
2D VS model from
surface wave
analysis
Updated geological
model
[Ivanov et al., 2006]
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Cavity detection
Example 1: void detection in a minerary area in
canada with pseudo-2D VS sections from surface
wave analysis
Example 2: (ERT) Electrical Resistivity
Tomography and (GPR) Ground Penetrating Radar
surveys reveal a sinkhole beneath a house
ERT
VS
Xu et al., 2008
GPR
Dobecki and Upchurch, 2006
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Identification of stratigraphic sequence / local litography
e.g. electrical methods to identify clays below sands
Turesson and Lind, 2005
Non-seismic methods:
Powerful tools to investigate lateral variations at the site
(e.g. for assessing the potential for differential settlements)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Hydro - geophysics
Sea coast-line
S1
Depth [ m ]
0
100
Land
Distance [ m ]
200
300
400
500
600
700
800
900
1000
1100
1200
-20
-40
-60
-80
-100
Salt water
intrusion
Macchiareddu - Cagliari (Italy)
50 150 250 350 450 550 650 750 850
[Ohm m]
Apparent resistivity pseudosection
Profile: n.3 Cardiga
2D rendering of time domain EM vertical 1D profiles for salt water
intrusion in coastal aquifer.
Courtesy of Alberto Godio
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Pollutants and waste detection
Buried waste disposal
Courtesy of Valentina Socco
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Hydrogeological / environmental applications
Electrical Resistivity Tomography (ERT)
(Martìnez-Pagàn et al., 2009)
Resistivity is sensitive to:
• pore fluid content
• pore fluid conductivity
Geocongress 2014 23rd February 2014
Saturated vs unsaturated(for coarse materials)
Identification and monitoring of plumes
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Monitoring in environmental applications
Example:
3D resistivity tomography on lab
soil samples for diffusion of
conductive plume monitoring.
(Comina et al., 2011).
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
site characterization
for engineering
purposes
Geotechnical
and geoenvironmental
site characterization
In the context of site characterization for engineering purposes,
the role of geophysical methods is twofold:
• evaluation of geometrical boundaries to model subsoil
conditions (e.g. stratigraphy but also physical inclusions or
hydrogeological features);
• evaluation of physical/mechanical parameters of direct use
for geotechnical modeling.
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Indagini geofisiche per la caratterizzazione geotecnica
Bulk waves
Longitudinal wave
(Primary wave – P)
VP =
λ + 2µ
=
ρ
M
ρ
Shear wave
(Secondary wave – S)
VS =
µ
=
ρ
G
ρ
(Animation courtesy of prof.Braile)
CISM Udine 7 Febbraio 2014
Sebastiano Foti
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Seismic methods
In a linear elastic medium
Shear wave propagation
G = ρ V S2
In soils
G 0 = ρ V S2
Animation courtesy of
Dr. L. Braile, Purdue University
Gsec
G0
τ
G0
1.0
Gsec
G0
Gsec
γc
γ
γc
γ
Strain range of
geophysical test
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Role of G0 in geotechnical engineering
• Evaluation of seismic site response
• Foundation vibrations
• Dynamic soil structure interaction
• Vibrations (e.g. railroads, industrial activities, …)
• Liquefaction suscettivity assessment
• Monitoring of ground improvement projects
• Correlation to operative values of G at medium strains
• Numerical simulations with advanced constitutive laws
• Evaluation of disturbance of soil samples
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Typical strain ranges for geotechnical problems
(Atkinson, 2000)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Site vs Lab (Pisa)
G0 lab (MPa)
120
90
60
30
0
0
30
60
90
120
G0 sito (MPa)
G 0 = ρ V S2
(Cross-Hole Test)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
G0 site vs lab
G
0
150
(m/s)
300
600
0.25
/G
0, sito
0.50
0.80 1.00
1.50
2.00
valori progetto
ROSRINE
campioni
ricostituiti
da sabbie
cementate
XX
X
X
X
X
campioni indisturbati
di sabbie cementate
V
S, sito
450
0.10
0,lab
750
900
Pisa
1050
0.0
valori usuali per
campioni di roccia
0.5
1.0
V
S,lab
G 0 = ρ V S2
Geocongress 2014 23rd February 2014
/V
1.5
S,sito
(Stokoe e Santamarina, 2000)
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Assessment of sample quality
The ratio
VS(lab) / VS(field)
Gives an indication of
sample quality
it can be used also for
coarse grained soils
DeGroot et al (2011).
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Site characterization for seismic projects
120
30
100
G (MPa)
80
25
G ( γ ) lab
20
60
15
40
10
20
5
0
0
0.0001
0.001
0.01
0.1
D (%) ∆ u/σ '0 (%)
G 0(sito ) = ρ ⋅ VS 2
1
Shear
strain γ (%)
deformazione
tangenziale,
• G0 from in situ testing (geophysics)
• G/G0(γ) e D(γ) from laboratory tests
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Site characterization for seismic projects
120
30
 G(γ) 
G ( γ ) = (G 0 ) sito ⋅

 G 0  lab
25
80
20
60
15
40
10
20
5
0
0
0.0001
0.001
0.01
0.1
D (%) ∆ u/σ '0 (%)
G (MPa)
100
1
deformazione
tangenziale,
Shear
strain γ (%)
In situ tests investigate a large volume of soil
whereas laboratory testing concerns small samples
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Body Waves
Compressional wave
VP =
M
ρ
(after Bolt, 1976)
Shear wave
VS =
G
ρ
Vertically polarized SV
or
Horizontally polarized SH
Geocongress 2014 23rd February 2014
(after Bolt, 1976)
Direction of Propagation
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
In a linear elastic isotropic homogeneous medium
VP =
M
ρ
VS: shear wave velocity
VP: dilational wave velocity
ρ: density
G: shear modulus
VS =
G
ρ
M: laterally constrained modulus
(oedometric conditions)
Note: In saturated soils VP is strongly influenced by the
compressibility of the pore fluid (water)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Biot Theory
Macroscopic approach: the medium is modeled as a binary continuum
arising from the superposition of a fluid and a solid phase occupying
simultaneously the same regions of space. The porosity is the link between
the two.
Hypothesis:
- isotropic, linear elastic soil skeleton
- a non-dissipative compressible fluid saturates all voids
- no relative motion between the solid and the fluid phases
(valid for low frequency range)
Writing the equations of motion for the porous media and applying the
Helmholtz decomposition, it is possible to show the existence of two
different compressional waves and of a unique shear waves.
The fastest compressional wave is called of the first kind or P-wave, the
slowest is called of the second kind or Biot wave.
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Biot solution
Under the hypothesis of grain incompressibility, the velocity of propagation of body
waves in porous media can be written as:
VP =
VS =
(K
SK
F
K
+ 43 ⋅ G ) +
n
ρ
G
ρ
where
ρ = (1 − n) ⋅ ρ S + n ⋅ ρ F
Geocongress 2014 23rd February 2014
ρS grain density
ρF water density
KF
KSK
G
n
water bulk modulus
soil skeleton bulk modulus
shear modulus
porosity
νSK Poisson ratio of the (evacuated) soil skeleton
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Biot solution
Under the hypothesis of grain incompressibility, the velocity of propagation of body
waves in porous media can be written as:
VP =
VS =
(K
SK
F
K
+ 43 ⋅ G ) +
n
ρ
G
ρ
ρ = (1 − n) ⋅ ρ S + n ⋅ ρ F
where
ρS grain density
ρF water density
KF
KSK
G
n
water bulk modulus
soil skeleton bulk modulus
shear modulus
porosity
νSK Poisson ratio of the (evacuated) soil skeleton
4 ⋅ (ρ S − ρ F ) ⋅ K F
ρ − (ρ ) −
 1 −ν SK  2
2
 ⋅VS
VP − 2 ⋅ 
SK 
 1 − 2ν 
n=
2 ⋅ (ρ S − ρ F )
S
S 2
Geocongress 2014 23rd February 2014
ρS,ρF,KF: standard values
VP & VS: measured
νSK : range 0.1÷0.4
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Soil porosity from seismic velocities
Leaning Tower of Pisa site
Velocity of Propagation [m/s]
0
500
0
1000
1500
2000
2500
Porosity
0
0.2
0.4
0.6
0.8
1
10
Depth [m]
20
30
40
50
60
70
80
LAB (Laval)
LAB (Osterberg)
Vp
Vs
Cross-Hole
test
(Foti et al., 2002)
Geocongress 2014 23rd February 2014
4 ⋅ (ρ S − ρ F ) ⋅ K F
ρ − (ρ ) −
 1 −ν SK  2
2
 ⋅VS
VP − 2 ⋅ 
SK 
ν
1
−
2


n=
2 ⋅ (ρ S − ρ F )
S
S 2
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Degree of saturation
Also very limited desaturation has a strong effect on the VP
Valle-Molina (2006)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Influence of degree of salutarion on liquefaction resistance
saturation degree strongly affect liquefaction resistance
VP can be used to monitor saturation and esclude liquefaction
Tsukamoto et al (2006)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Example: Zelasny Most tailing dam
West dam
Geocongress 2014 23rd February 2014
Jamiolkowski, 2012
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Non seismic methods
Quantitative use of geophysical parameters other than
seismic velocities is less straightforward and typically require
the use of empirical correlations with geotechnical parameters
Example: electrical conductivity of soils
Trasport parameter related to:
-
fluid properties (solubility of ionic species, concentration);
-
mineralogy and specific surface of the solid grains;
-
porosity and fabric
σw : pore fluid conductivity
Archie
σt = σw nm Srp
n: porosity
Bruggeman
σt = σw n3/2
m = 3/2 : theoretical
Waxman & Smits
σt = X (σw + σs)
σs : clay surface conductivity
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
S: saturation
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Example at Lab scale
Polito – 2D ERT (Borsic et al., 2005)
Identification of zones with different compaction
levels in sand
Coarse Matrix
n ≈ 0.48
Dense Inclusion
n ≈ 0.43
1
mS/cm
2
Estimated values with
Bruggeman equation
Matrix n ≈ 0.46
3
Inclusion n ≈ 0.42
3. Tomographic reconstruction
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
In-hole vs surface methods
(Invasive vs Non-invasive methods)
2
1
D
Cross-Hole Test (CHT)
Down-Hole Test (DHT)
Seismic Cone (SCPT)
Seismic Dilatometer (SDMT)
P-S Suspension Logging
Vertical Seismic Profiling (VSP)
Geocongress 2014 23rd February 2014
X
3
n
X
Surface Waves Methods SWM
(SASW, MASW, microtremors)
Seismic Refraction
(P-waves or SH-waves)
Seismic Reflection
(P-waves or SH-waves)
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
In-hole vs surface methods
Invasive Tests
Direct measurements: simple
and accurate interpretation
Costs and flexibility (in time and
space)
Good resolution also at great
depth
Non-intrusive (e.g. important for
waste landfills)
Easier standardization
Average properties (dynamic
behaviour of the whole soil
deposit)
Additional information from
borehole logging or the
penetration of the cone
Disadvantages
Non-Invasive Tests
Costs and necessity of planning
well in advance
Local measurement
Geocongress 2014 23rd February 2014
Large volumes are investigated
Complex interpretation (indirect
measurements based on inversion
procedures or heavy processing)
Accuracy and resolution at depth
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
In-hole vs surface methods
\
VS1
VS2
VS3
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Combined use of geophysical methods
Synergies between different techniques can be
exploited at different level of integration:
• Level 1: comparison for validation / calibration
• Level 2: data integration and data fusion (combining
different information on the same medium)
• Level 3: a priori info (one method help the other)
• Level 4: joint inversion (simultaneous interpretation of
different dataset)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Level 1: Comparison In-Hole methods vs SASW
Vs (m/s)
0
400
800
0
0
400
5
Saluggia
5
Depth (m)
Depth (m)
Cross Hole
25
Cross Hole
25
SASW-fk
SASW-fk
30
30
Vs (m/s)
400
Vs (m/s)
Vs (m/s)
800
0
400
800
0
1200
400
800
1200
1600
0
0
0
Pontremoli
Pontremoli
Pontremoli
site 1
site 2
site 3
Depth (m)
5
10
5
Depth (m)
0
Pisa
5
20
Down Hole
15
300
15
20
SASW-fk
200
10
10
10
100
0
15
Depth (m)
0
800
0
Castelnuovo
Depth (m)
Vs (m/s)
Vs (m/s)
10
5
10
Down Hole
Down Hole
Down Hole
SASW-fk
SASW-fk
SASW-fk
15
Geocongress 2014 23rd February 2014
15
15
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Level 2: Data integration and data fusion
Electrical resistivity tomography
Pugin et al., 2009
SH-wave seismic reflection
resistivity
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Combined use
•
•
•
•
Level 1: comparison for validation
Level 2: data fusion
Example: synergies of
Level 3: a priori info
seismic refraction and
Level 4: joint inversions
surface wave analysis
(SWM)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Example of synergy: SW + VP refraction
Same testing setup and
equipment
P-waves
VP1
Receivers (geophones)
VP2
VP3 ≈ VP2
Rayleigh waves
Experimental data contain both
surface waves and
direct/refracted P waves
VS1
VS2
VS3
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
P-WAVE REFRACTION
70
travel time [ms]
60
50
VP1
40
VP2
30
VP3 ≈ VP2
20
10
0
10
20
30
40
50
spacing [m]
60
KF
( K + ⋅ G) +
n
VP =
S
(1 − n) ⋅ ρ + n ⋅ ρ F
SK
4
3
Shallow water table masks
variation of the mechanical
properties of the solid skeleton
(influence of the pore fluid)
Geocongress 2014 23rd February 2014
70 0
0
80
200
400
600
Vp [m/s]
800
1000
1200
1400
1600
Water table
1
Depth [m]
0
2
3
4
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Experimental Data
Shear Wave Velocity (m/s)
Depth (m)
5
0
200
400
600
800 1000
starting profile
inversion #1
inversion #2
inversion #3
cross-hole test
10
15
20
25
Hp#1 Water table from P-wave refraction
Hp#2 No water table
Hp#3 Water table deeper than Hp #1
700
phase velocity, m/s
0
experimental
inversion #1
inversion #2
inversion #3
600
500
400
300
200
100
0
10
20
30
40
50
60
70
frequency, Hz
30
Geocongress 2014 23rd February 2014
(Foti and Strobbia, 2002)
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Level 4: joint inversion
(Piatti et al., 2012b)
A single inversion problems is solved considering all the available
experimental information: the best fit parameters for both VP and VS
models are obtained
A single misfit parameter include misfit on Rayleigh wave dispersion
curve and P-wave travel times
T -1
1


+
L=
d
-g
m
C
d
-g
m
(
)
(
)
(
)
(
)
obs
obs
obs

 N + M + A 
dobs = ( log (VR1 ) , log (VR2 ) , ...., log (VRN ' ) ) ( log ( t1 ) , log ( t2 ) , ...., log ( tN '' ) ) 
gSW ( m) 
g ( m) = 

gPR ( m) 
m = [( log ( h1 ) , log ( h2 ) , ...., log ( hn ) ) ( log (VS1 ) , log (VS 2 ) , ...., log (VS n+1 ) )
Geocongress 2014 23rd February 2014
( log (V ) ,
P1
SEBASTIANO FOTI
log (VP2 ) , ...., log (VPn+1 ) )]
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Example on synthetic data
(Piatti et al., 2012b)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Example on synthetic data
(Piatti et al., 2012b)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Experimental data
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Case History #1
Combination of seismic and electrical methods for
the assessment of site conditions for seepage
analysis along an embankment
• Combination of several methods for reliable evaluation of cover
thickness
• Joint inversion to improve accuracy
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
The PO river
LENGTH: 650 km
DISCHARGE
ave.= 1450 m3/s
max.= (nov 2000): 13000m3/s
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Seepage potential
Floods very often start with localized seepage that can
degenerate causing inundations
10 extreme events each 100 years
Levees for a total length over 2400 km
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Seepage potential
Geology: alluvial deposits: recent sands, gravel, clay
TARGET: clayey layer: continuity, thickness
Water level can reach 10 m
above the ground surface
Anthropic soil
?
Thickness of low permeability layer?
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Geophysical investigation
large extension of the areas
Interest in fast geophysical tests from the surface
At a test site several
methods have been
tested and compared
VES ERT
HEP
SWM
Prefr SHrefr
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Combinations MASW + VES
ρ
app
Processing
VR
Apparent
resistivity
AB
dispersion
curve
ω
Inversion
VS
ρ 0
G0
VS profile
Z
resistivity
profile
Geocongress 2014 23rd February 2014
Z
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Joint inversion VES + MASW
Physical parameters: shear velocity and resistivity
Assumed parameter distribution: stack of homogeneous isotropic layers
VS, ρ
MODEL PARAMETERS:
n
ρ
n
VS
n-1 H
VS, ρ
LINK BETWEEN THE TWO MODELS:
geometry, thickness of the layers
(same position of interfaces: independent
variations of the two parameters, a variation
of resistivity does not imply a variation of
seismic shear velocity )
VS, ρ
VS, ρ
From 4n-2 to 3n-1 unkowns
with the same experimental information
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Joint inversion VES + MASW
ρ
Processing
app
AB
VR
Joint Inversion
ρ
0
VS
dispersion
curve
ω
G0
VS profile
Z
Geocongress 2014 23rd February 2014
Z
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Field test results
(Comina et al., 2004)
Resistivity curves
3
10
Electric stratigraphy
Seismic stratigraphy
0
0
joint
joint
2
10
1
10 0
10
10 1
AB/2 [m]
10
2
-2
-2
-4
-4
-6
-6
-8
-8
Best estimate
single
of the -10
clay layer
thickness
-12
z [m]
joint
-10
Dispersion curves
350
Vr [m/s]
single
-12
300
single
250
joint
200
z [m]
ro [Ohm.m]
single
-14
-14
-16
-16
150
-18
1000
5
10
15
freq [Hertz]
20
Geocongress 2014 23rd February 2014
25
30
50 100 150 200 250 300
ro [Ohm.m]
SEBASTIANO FOTI
-18
100 200 300 400 500
Vs [m/s]
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
The joint inversion algorithm
(Garofalo, 2014)
Surface-wave
propagation
Body-wave
propagation
Electrical
current flow
Dispersion
curves
P-wave
travel-time
Apparent
resistivity
Joint
Inversion
Physical Laws
Physical
inversion
Thickness h
S-wave velocity VS
m
P-wave velocity Vp
Resistivity Rho
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
dobs
reference model
Geometric
regularization
Structural
inversion
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
The model and the physical links
Unsaturated – skeleton
h1
VS,1
VP,1
Rho1
Saturated – mixture
of fluid and grains
h2
VS,2
VP,2
Rho2
VS,3
VP,3
Rho3
Half Space – clay
Model Parameters m:
Geocongress 2014 23rd February 2014
thickness (h)
density (ρ)
S-wave velocity (VS)
P-wave velocity (VP)
Resistivity (Rho)
SEBASTIANO FOTI
Poisson’s ratio
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
The model and the physical links
Unsaturated – skeleton
h1
VS,1
VP,1
Rho1
Saturated – mixture
of fluid and grains
h2
VS,2
VP,2
Rho2
VS,3
VP,3
Rho3
Half Space – clay
Model Parameters m:
Seismic wave Velocities
Porosity
Resistivity
Geocongress 2014 23rd February 2014
thickness (h)
density (ρ)
S-wave velocity (VS)
P-wave velocity (VP)
Resistivity (Rho)
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Field example: Results
(Garofalo, 2014)
Saturated layer
4 ⋅ ( ρ S − ρ F )⋅ K F
ρ − (ρ ) −
 1− ν SK  2
2
⋅VS
VP − 2 ⋅ 
SK 
 1− 2ν 
φ=
2 ⋅ (ρ S − ρ F )
S
S 2
φ = 0.42 – 0.52
Bellotti and Selleri (1969)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Field example: Results
(Garofalo, 2014)
ν = 0.13
Bergamo and Socco
(2013)
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Case history #2
Investigation of volcanoclastic slopes
• Combination of several in situ geophysical tests to increase the
reliability of the results
• Combination of laboratory and in situ testing for the assessment of
saturation conditions
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Flowslides of 1998 in Campania
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Sarno
Air-fall pyroclastic deposits
(Cascini et al., 2008)
flowslides occurred in May 1998
(Cascini et al., 2008)
Cover soils formed by volcanic ashes from the Vesuvio
(few meters thick) over a carbonatic bedrock
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Site characterization
Objectives
• Quantification of potential volume of the flow (for the
design of mitigation infrastructures): thickness of the
soil cover
• Prevision of onset of the flowslide: assessment and
monitoring of saturation condition of the soil cover
Critical issues
• Very difficult site logistics with steep and vegetated
slopes poses strong limitations in the use of conventional
site tests (boreholes and penetration testing)
• Necessity of investigating large areas
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Combination of different geophysical approaches
Surface wave method (MASW)
(Cosentini et al., 2012)
Seismic tomography (VP)
Electrical resistivity tomography
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Comments
• Electical and seismic (VP) tomography show that the
assumption of a layered medium in MASW is reasonable
• Inversion of MASW shows the relevance of higher
modes at this site: surface wave analysis is not a simple
and straightforward task
• The estimated thickness of the cover material is
comparable with different methods
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
(Cosentini et al., 2012)
Laboratory calibration of Archie’s law for unsat materials
σt = σw nm Srp
n: porosity
S: saturation
σw : pore fluid conductivity
The two exponet m and p are found by fitting laboratory data
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
(Cosentini et al., 2012)
Mapping resistivity into degree of saturation
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Closing Remarks
• Geophysical test provide useful tools for
geotechnical site characterization
– evaluation of geometrical boundaries to model subsoil
conditions (e.g. stratigraphy but also physical
inclusions or hydrogeological features);
– evaluation of physical/mechanical parameters of
direct use for geotechnical modeling.
• VS G0; sample quality
• VP saturation; porosity (+M0 ν for dry soils)
• Surface wave methods are cost and time
effective but their interpretation is not simple
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Closing remarks
• Importance of choosing the right technique
for the specific application
• Integration of different techniques reduces
uncertainties
• Laboratory experimental can provide a
framework and calibration for quantitative
interpretation of field tests
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO
Geophysical methods for geotechnical site characterization
Thank you for your attention
Acknowledgments
Prof. Laura Valentina Socco (DIATI - Politecnico di Torino)
Dr Cesare Comina (University of Torino)
Prof. Guido Musso (DISEG – Politecnico di Torino)
Dr Renato Cosentini (Politecnico di Torino)
Ms Flora Garofalo (PhD student at Politecnico di Torino)
Dr Paolo Bergamo (now at Queen’s University, Belfast – UK)
Dr Margherita Maraschini (now at Fugro - UK)
Dr Daniele Boiero (now at Western-Gico - UK)
Dr Claudio Piatti (now at D’Apollonia - Italy)
Dr Claudio Strobbia (now at Western-Gico - UK)
Additional material available at
www.soilmech.polito.it/download
Geocongress 2014 23rd February 2014
SEBASTIANO FOTI
POLITECNICO DI TORINO