Chapter I, Introduction To The Theory of Poro-Elasticity and Two Applications
Chapter I, Introduction To The Theory of Poro-Elasticity and Two Applications
poro-elasticity
and two applications.
Lecture 1
24 March 2009
Introduction
Poro-elasticity theory
References
Part I : Introduction
Objectives
A porous medium.
Pression solution : a deformation mechanism (Houseknecht and
Hathon, 1987).
Part II : Poro-elasticity theory
Introduction
• Hypothesis : isotropy.
Part II : Poro-elasticity theory
Recall Hooke’s law
1
∆σ = (3K − 2G) δ tr() + 2G , (1)
3
K : compressibility and G : shear modulus
• The French notation !
a is a vector of coordinate ai ,
A is a second-order tensor of coordinate Aij
δ second-order identity tensor δij (Kronecker delta 1 if i = j, zero
otherwise)
tr() the trace operator ii
Repeated indices imply summation : ii = 11 + 22 + 33
Note : δii = δ11 + δ22 + δ33 = 3
Part II : Poro-elasticity theory
Recall Hooke’s law
1
∆σ = (3K − 2G) δ tr() + 2G ,
3
• Interpreting the compressibility modulus K
3
2
1
1 1 +ε
33
1 + ε 11
1
1 1 + ε 22
Part II : Poro-elasticity theory
Recall Hooke’s law
∆σ = C
I: (3)
... there exists a fourth-order tensor : C
I = 3KJI + 2GIK
• Expressed in terms of two fundamental tensors
1
JI = δ⊗δ , IK = 1IIIS − JI (4)
3
JI : JI = JI , IK : IK = IK et IK : JI = O
I , (5)
E E
K = G= . (6)
3(1 − 2ν) 2(1 + ν)
∆σ = −αδ∆pf , (9)
α : Biot’s constant
The sign convention.
∆ pf
Part II : Poro-elasticity theory
Extension of Hooke’s law to saturated media
1 1 1 1 α
=( − ) δtr(∆σ) + ∆σ + δ∆pf . (11)
3K 2G 3 2G 3K
Part II : Poro-elasticity theory
Extension of Hooke’s law to saturated media
a) b)
−∆ pf n −∆ pf n
∆p
∆p f
f
=
Ks
• Consider a porous solid which skeleton is composed of an elastic,
compressible material Ks .
Pore network is connected and fluid saturated.
• Apply fluid pressure change ∆pf and a density of external forces
T d = −∆pf n
• Two analyses of the same test :
(1) at the level of the micro-structure
(2) from a macroscopic level (ignore pore level description)
Part II : Poro-elasticity theory
Interpretation of Biot’s coefficient α
∆σ = −∆pf δ , (13)
• The solid phase local deformation is purely spherical
1
=− ∆pf δ , (14)
3Ks
• The fluid phase sustains the same deformation, thus homogeneous
deformation for the whole.
Part II : Poro-elasticity theory
Interpretation of Biot’s coefficient
1 α
= δtr(∆σ) + ∆pf δ . (15)
9K 3K
• Approaches (1), eq. (14) and (2), eq. (15) are identical :
K
α=1− . (16)
Ks
a) b)
−P n −P n
∆p
f
• First test is drained, force density on boundary T d = −Pn and
change in fluid pressure ∆pf .
1 α
=− Pδ + ∆pf δ . (17)
3K 3K
• Second test :
1
=− Pδ . (18)
3Ku
• Adjust the pressure change of first test so that the deformations
(17) and (18) are identical :
Ku − K
∆pf = P . (19)
αKu
Part II : Poro-elasticity theory
Change in fluid pressure for undrained tests
• Need to generalize
B Ku − K
∆pf = −tr(∆σ) with B = , (20)
3 αKu
Skempton’s coefficient B.
Between 0 et 1.
Tends to 0 for very compressible fluids, tends to 1 if solid and fluid
have the same compressibility.
Part II : Poro-elasticity theory
Reciprocity or Maxwell-Betti theorem extended to poro-elasticity
T 1d T d2
b b2
∆ pf 1
1
∆ pf
2
• Theorem
T d1 · u 2 dS +
R R
∂ΩT Ω
b1 · u 2 + α∆p1 tr(2 )dV =
T d2
R R
∂ΩT
· u 1 dS + Ω
b2 · u 1 + α∆p2 tr(1 )dV . (21)
I : 1 − α´p1 tr(2 )
Material response : ∆σ 1 : 2 = 2 : C
The position of the shear zone in California with respect to the San
Andreas fault (Yeats et al., 1997).
Part III : Application to post-seismic rebound
The earthquake of Landers, 1992
Interferogram with two SAR images before and after the Landers
earthquake , June 1992 (Yeats, Sieh and Allen 1997, after Massonet
et al. 1993)
Part III : Application to post-seismic rebound
Following Peltzer et al. (1996)
F IG .:
Three sets of SAR images are used corresponding to (1) 7/08/92 and
24/09/95 (red on Figure 29), (2) 27/09/92 and 23/01/96 (green) and
(3) 10/01/93 and 23/05/95 (blue).
Part III : Application to post-seismic rebound
Following Peltzer et al., (1996)
w = (νu − ν)1 h .
l =5 km, h = 4 km, δl = 3 m, the two Poisson’s ratios differ by 0.03, so
delayed displacement by as much as 7 cm !
Part V : References