CN°1

Sismique L3
Cours de sismique destiné au troisième

année licence en géophysique
SAMAI S. [email protected] USTHB 2022-2023

1- Théorie des ondes et du rai sismique
Onde sismique:
L’onde sismique est la propagation dans un milieu continu d’une déformation
brusque produite lors d’un choc ou par une vibration mécanique transitoire. Ce
phénomène est induit par la nature dispersive des milieux élastiques. Les ondes
sismiques peuvent être produites naturellement par un tremblement de Terre
(séisme) ou par ébranlement artificiel à la surface de la Terre
Milieu continue:
C’est une notion fondamentale en géophysique. Il est défini comme une

répartition, homogène ou non, d’un matériau physique, caractérisé entre autres,
par sa densité et ses propriétés élastiques (vitesses sismiques). Le milieu
physique continu peut être considéré à des échelles qui varient de celui de la
planète Terre toute entière jusqu’au modèle géologique réduit (échantillons de
roches).
SAMAI S [email protected] USTHB 2022-2023

1-Théorie des ondes et du rai sismique
Un milieu physique continu est représenté par une répartition dans l’espace
d’une matière renfermée dans une surface fermée (S) et occupant un volume
(V) dans un système lié à un repère.
Schéma d’un milieu continu dans l’espace

Théorie des ondes et du rai sismique
Exemples de milieux continus
- La Terre entière renfermée dans sa

surface extérieure.
- Les constituants de la Terre séparés

par des discontinuités. (Atmosphere,
Lithosphère, Manteau, Noyau)
- Structures géologiques (bassin sédimentaire,

réservoir pétrolier, gisement aurifère, nappes
aquifères.
- Modèle géologique réduit, Echantillon de

roche

1- Élasticité:
Un milieu est dit élastique s’il est capable de retrouver son état initial après
avoir subi une déformation de contrainte auquel il a été soumis. Les lois les
plus simples de l’élasticité ont été établies en considérant un milieu semi-
infini de propriétés élastiques uniformes. Cette propriété est fondamentale
dans le domaine de la sismologie et de la sismique du fait qu’elle contrôle
les mécanismes et l’équation de propagation de l’onde sismique.

# Contrainte
Appelée aussi champ de contrainte. C’est une force par unité de surface. Deux cas
de figures peuvent se présenter, contrainte normale et tangentielle
𝑭𝑁
En SI: 𝝈, 𝝉 = 𝟐
𝑃𝑎 = 10−𝟓 𝐵𝑎𝑟
𝑆𝑚

Déformation
- C’est un changement de forme, appelé aussi champ ou potentiel de
déformation.
- Déplacement relatif des particules au sein d’un milieu. On distingue deux types
de déformations qui correspondent, chacune à un type de contrainte, la
dilatation ou le cisaillement,

Dilatation ou compression:
𝐴′ 𝐵′ − 𝐴𝐵
ε𝑥 =
𝐴𝐵
𝑑𝑥 − 𝑢 + 𝑢 + 𝑑𝑢 − 𝑑𝑥
=
𝑑𝑥
𝑑𝑢
=
𝑑𝑥

Nous reprenons le même raisonnement pour les deux autres dimensions pour
déterminer la dilatation pour le cas d’un volume 3D
Avant déformation
𝑉 = 𝑑𝑥𝑑𝑦𝑑𝑧 (Cube à l’intérieur)

Après déformation
𝑉′ = 𝑑𝑥𝑑𝑦𝑑𝑧(1 + ε𝑥 ) (1 + εy ) (1 + εz )
Nous déterminons la déformation du volume que nous appelons déformation
totale par:
𝑉′ − 𝑉
𝜽=
𝑉
𝑑𝑥𝑑𝑦𝑑𝑧 1 + ε𝑥 )(1+ ε𝑦 1 + ε𝑧 − 𝑑𝑥𝑑𝑦𝑑𝑧
𝜽=
𝑑𝑥𝑑𝑦𝑑𝑧
𝑑𝑥𝑑𝑦𝑑𝑧 1 + ε𝑥 )(1+ ε𝑦 1 + ε𝑧 − 𝑑𝑥𝑑𝑦𝑑𝑧

𝜽=
𝑑𝑥𝑑𝑦𝑑𝑧
𝜽 = ε𝑥 + ε𝑦 + ε𝑧

𝜕𝑢 𝜕𝑣 𝜕𝑤
𝜽= + +
𝜕𝑥 𝜕𝑦 𝜕𝑧
Déformation en cisaillement:
Nous aurons pour ce type de déformation les expressions suivantes:

Il s’agit de déformation infinitésimale, cela veut dire que:
𝛼 𝛼
tg ±𝜑 ~ ±𝜑
2 2
Nous pouvons écrire alors:
𝛼 𝜕𝑢
+𝜑 =
2 𝜕𝑦
𝛼 𝜕𝑣
−𝜑 =
2 𝜕𝑥
La somme de ces deux équations membre à membre donne la
déformation en cisaillement dans le plan xy:
𝜕𝑢 𝜕𝑣
𝛼= +
𝜕𝑦 𝜕𝑥

Le cisaillement est associé à une rotation que nous pourrons tirer

du même système d’équation par une soustraction des ses deux
membres:
1 𝜕𝑢 𝜕𝑣
𝜑= −
2 𝜕𝑦 𝜕𝑥
Exercice:
Retrouvez les expression du cisaillement et de rotation dans les
plans XZ et YZ

Loi de Hook
Une relation linéaire entre la contrainte et la déformation.
E: Module de
d’élasticité
longitudinale appelé
module de Young.
ρ: La densité; λ: Le
premier coefficient de
Lamé; µ: Le deuxième
coefficient de Lamé,
appelé aussi
coefficient de rigidité
ou module d’élasticité
transversale G;

Rapport de déformation longitudinale entre deux sens orthogonaux
Le coefficient de poisson:
𝜺𝒚𝒚
𝝂=−
𝜺𝒙𝒙
𝜺𝒚𝒚 = −𝝂𝜺𝒙𝒙
1
𝜺𝒙𝒙 = 𝝈𝒙𝒙
𝐸
𝝂
𝜺𝒚𝒚 = − 𝝈𝒙𝒙
𝐸
Analyse élargie au 3D:
Dans un système de coordonnées de 3D, XYZ (cartésien) la loi de Hook, pour un

milieu homogène et isotrope (cas le plus simple), est une relation linéaire de
tenseurs. Elle permet de décrire dans l’espace l’évolution de du champ de
déformation sous l’effet du champ de contrainte. Dans l’espace (XYZ), le champ
de contrainte est décomposé en trois composantes normale de dilatation-
compression (𝝈𝒙𝒙 , 𝝈𝒚𝒚 , 𝝈𝒛𝒛 ) à chacune d’elles sont associés deux composantes
de cisaillement
(𝜏𝑥𝑦 , 𝜏𝑥𝑧 ), (𝜏𝑦𝑥 , 𝜏𝑦𝑧 ), (𝜏𝑧𝑥 , 𝜏𝑧𝑦 ).
Le champ de déformation se définie dans ce même espace par des dilatations
compressions dans les directions normales

Potentiel de dilatation 3D:
𝜕𝑢
1 𝜺𝒙𝒙 =
𝜺𝒙𝒙 = 𝝈𝒙𝒙 𝜕𝑥
𝐸 𝜕𝑣
1 𝜺𝒙𝒙 =
𝜺𝒚𝒚 = 𝝈𝒚𝒚 𝜕𝑦
𝐸
1 𝜕𝑤
𝜺𝒛𝒛 = 𝝈𝒛𝒛 𝜺𝒙𝒙 =
𝐸 𝜕𝑧

Potentiel de cisaillement / rotation :
Cisaillement :
𝜕𝑢 𝜕𝑣
𝜺𝒙𝒚 = +
𝜕𝑦 𝜕𝑥
𝜕𝑣 𝜕𝑤
𝜺𝒚𝒛 = +
𝜕𝑧 𝜕𝑦
𝜕𝑤 𝜕𝑢
𝜺𝒛𝒙 = +
𝜕𝑥 𝜕𝑧
𝜕𝑢 𝜕𝑣
α= −
α(𝑥, 𝑦, 𝑧) 𝜕𝑦 𝜕𝑥
𝜕𝑣 𝜕𝑤
Rotation : φ 𝛽(𝑥, 𝑦, 𝑧) , 𝛽= −
𝜕𝑧 𝜕𝑦
γ(𝑥, 𝑦, 𝑧) 𝜕𝑤 𝜕𝑢
γ= −
𝜕𝑥 𝜕𝑧

Tenseur de contrainte tenseur de déformation

Relation entres propriété élastique:
D’après la formule de Poisson, la dilatation dans l’axe XX induit une

compression dans deux autres axes YY et ZZ.
La dilatation selon XX sera alors la somme de la composante

principale due à 𝝈𝒙𝒙 et des deux autres induit par 𝝈𝒚𝒚 et
𝝈𝒛𝒛 :
1 𝝂 𝝂
𝜺𝒙𝒙 = 𝝈𝒙𝒙 − 𝝈𝒚𝒚 − 𝝈𝒛𝒛
𝐸 𝐸 𝐸
Nous reprenons le même raisonnement pour les deux axes
YY et ZZ, nous obtenons le système d’équation:

𝟏 𝝂 𝝂
− −
𝜺𝒙𝒙 𝑬 𝑬 𝑬 𝝈𝒙𝒙
𝝂 𝟏 𝝂
𝜺𝒚𝒚 = − − 𝝈𝒚𝒚 (1)
𝑬 𝑬 𝑬
𝜺𝒛𝒛 𝝂 𝝂 𝟏 𝝈𝒛𝒛
− −
𝑬 𝑬 𝑬
Module de compressibilité K:
Prenons le cas d’égalité du trois composantes en état de
compression:
𝝈𝒙𝒙 = 𝝈𝒚𝒚 = 𝝈𝒛𝒛 = −𝒑

La somme des trois équations du système membre à membre

nous obtenons:
𝐸(ε𝑥x + ε𝑦y + ε𝑧z ) = (1 − 2𝝂)(𝝈𝒙𝒙 + 𝝈𝒚𝒚 + 𝝈𝒛𝒛 )

𝐸𝜽 = (1 − 2𝝂)(−𝟑𝒑)
Le module de compressibilité est lié linéairement à la pression
total p et la dilatation totale par la relation:
−𝐾𝜽 = 𝒑
Il en résulte que:
𝑬
𝑲=
𝟑(𝟏 − 𝟐𝝂)
Relation entre le module de rigidité et le coefficient de Poisson:
𝒅𝟎 = 𝒂 𝟐
Conformément à cette géométrie on peut écrire que :
𝒅𝟐 = 𝒂𝟐 + (𝒂 + 𝒂𝒕𝒈𝝋)𝟐

𝒅𝟐 = 𝒂𝟐 + 𝒂𝟐 (𝟏 + 𝟐𝒕𝒈𝝋 + 𝒕𝒈𝟐 𝝋)
𝟐
𝟏 𝟐
= 𝟐𝒂 𝟏 + 𝒕𝒈𝝋 + 𝒕𝒈 𝝋
𝟐
~𝝋
𝟏
𝒅 = 𝒅𝟎 (𝟏 + 𝝋)𝟐
Que nous pouvons l’écrire :
𝟏
𝒅 ≈ 𝒅𝟎 𝟏+ 𝝋
𝟐

L’allongement relatif selon la diagonal correspondra au

cisaillement :
∆𝒅 𝒅 − 𝒅𝟎
=
𝒅𝟎 𝒅𝟎
∆𝒅 𝒅 − 𝒅𝟎
=
𝒅𝟎 𝒅𝟎
1
𝒅𝟎 + 𝒅𝟎 𝝋 − 𝒅𝟎
= 2
𝒅𝟎
𝟏
= 𝝋
𝟐

Origine du cisaillement ?:
C’est la variation des contraintes dans les direction orthogonales
qui est responsable de la création des forces de cisaillement :
La variation de contrainte ∆p fluctue autour d’une valeur

moyenne p tel que:
𝟏
𝒑 = 𝝈𝒙𝒙 + 𝝈𝒚𝒚
𝟐
La variation ∆p sera alors définie par :
∆𝒑 = 𝒑 − 𝝈𝒙𝒙
= 𝝈𝒚𝒚 − 𝒑
C’est cette même ∆p qui vaut la contrainte de cisaillement
𝝉𝒙𝒚 que nous définissons par:
∆𝒑 = 𝝉𝒙𝒚 = µ𝝋
Cette variation de contrainte ∆𝒑 provoque un déplacement
extension selon l’axe x:
∆𝒙 ∆𝒑
=
𝒙 𝑬

La compression dans l’axe Y due à cette même variation de
contrainte ∆𝒑 contribue elle aussi à une extension selon l’axe X
via la relation de Poisson:
∆𝒙 𝝂∆𝒑
=
𝒙 𝑬
La résultante des deux dernières dilatations relatives est identique
à l’allongement relatif selon la diagonale d
∆𝒅 𝟏 ∆𝒑
= 𝝋= 𝟏+𝝂
𝒅 𝟐 𝑬
Et du fait que: ∆𝒑 = 𝝉𝒙𝒚 = µ𝝋
𝑬 µ Coefficient de rigidité:
On déduit : µ=
𝟐 𝟏+𝝂 µ 1ier coefficient de Lamé
Le deuxième coefficient de Lamé 𝝀:

Prenons la première équation du système (1):
𝑬𝜺𝒙𝒙 = 𝝈𝒙𝒙 − 𝝂𝝈𝒚𝒚 − 𝝂𝝈𝒛𝒛
Nous pouvons l’écrire sous cette forme:
𝑬𝜺𝒙𝒙 = (𝟏 + 𝝂)𝝈𝒙𝒙 − 𝝂𝝈𝒙𝒙 − 𝝂𝝈𝒚𝒚 − 𝝂𝝈𝒛𝒛

= (𝟏 + 𝝂)𝝈𝒙𝒙 − 𝝂(𝝈𝒙𝒙 + 𝝈𝒚𝒚 + 𝝈𝒛𝒛 )
Qui s’écrit aussi:

𝟏+𝝂 𝑬
𝝈𝒙𝒙 + 𝝈𝒚𝒚 + 𝝈𝒛𝒛 = 𝝈𝒙𝒙 − 𝜺𝒙𝒙
𝝂 𝝂

𝐸(ε𝑥x + ε𝑦y + ε𝑧z ) = (1 − 2𝝂)(𝝈𝒙𝒙 + 𝝈𝒚𝒚 + 𝝈𝒛𝒛 )

Qui est équivalente à:
𝐸𝜃
(𝝈𝒙𝒙 + 𝝈𝒚𝒚 + 𝝈𝒛𝒛 ) = 𝐴𝑣𝑒𝑐 𝜃=ε𝑥x + ε𝑦y + ε𝑧z
(1 − 2𝝂)
On déduit alors que:
𝟏+𝝂 𝑬 𝐸𝜃
𝝈𝒙𝒙 − 𝜺𝒙𝒙 =
𝝂 𝝂 (1 − 2𝝂)
𝝂𝑬 𝑬
𝝈𝒙𝒙 = 𝜽+ 𝜺𝒙𝒙
(𝟏 + 𝝂)(𝟏 − 𝟐𝝂) (𝟏 + 𝝂)

C’est cette dernière équation qui nous permet de définir le deuxième
coefficient de Lamé :
𝝂𝑬
𝝀=
(𝟏 + 𝝂)(𝟏 − 𝟐𝝂)
Et la première équation de la loi de Hook pour les contraintes normales
dans un milieu homogène et isotrope s’écrit:
𝝈𝒙𝒙 = 𝝀𝜽 + 𝟐µ𝜺𝒙𝒙
De la même manière, on établi les deux autres équations, puis on
constitue la loi tensorielle de Hook :
𝝀𝜽 + 𝟐µ𝜺𝒊𝒋 𝒊=𝒋
𝝈𝒊𝒋 = ൝
𝟐µ𝜺𝒊𝒋 𝒊≠ 𝒋
𝑂ù: 𝑖, 𝑗 = 𝑥, 𝑦, 𝑧
Propagation d’une déformation:
Pour commencer, intéressons-nous au cas d’une propagation d’une
déformation dans filon de roche représenté en cette figure :
L’état de contrainte dans ce milieu dépendra de deux paramètres: la

distance et le temps. Prenons l’élément de masse dm situé entre x et
dx. Cet élément à instant t et un endroit donné x est soumis au deux
forces F(x) et F(x+dx). Nous appliquons la loi de la dynamique

෍ 𝑭𝒊 = 𝑚𝒂
𝑭 𝒙 + 𝒅𝒙 − 𝑭 𝒙 = 𝒂𝑑𝑚
Fixons maintenant les termes de cette équation:
𝜕𝐹
𝐹 𝑥 + 𝑑𝑥 =F(x)+𝜕𝑥 𝑑𝑥
Nous avons aussi :
𝜕2 𝑢(𝑥,𝑡)
𝐹 = 𝜎𝑆, 𝑑𝑚 = 𝜌𝑑𝑣 = 𝜌𝑆𝑑𝑥, 𝑎=
𝜕𝑡 2
En remplaçant dans l’équation nous obtenons :
𝜕𝐹
F(x)+𝜕𝑥 𝑑𝑥 − 𝐹 𝑥 = 𝜌𝑆𝑑𝑥𝒂
𝜕𝜎 𝜕 2 𝑢(𝑥, 𝑡)
=𝜌 (2)
𝜕𝑥 𝜕𝑡 2
Par l’application de la loi de Hook linéaire:

𝜕𝑢(𝑥,𝑡)
𝜎=𝐸𝜀= 𝐸
𝜕𝑥
En remplaçant dans (2) nous obtenons:
𝜕 2 𝑢(𝑥, 𝑡) 𝜕 2 𝑢(𝑥, 𝑡)
𝐸 2
=𝜌
𝜕𝑥 𝜕𝑡 2
𝜕 2 𝑢(𝑥, 𝑡) 𝜌 𝜕 2 𝑢(𝑥, 𝑡)
2
=
𝜕𝑥 𝐸 𝜕𝑡 2
𝜕 2 𝑢(𝑥, 𝑡) 1 𝜕 2 𝑢(𝑥, 𝑡)
2
= 2 (3)
𝜕𝑥 𝑉 𝜕𝑡 2

𝐸
𝐿𝑒 𝑝𝑎𝑟𝑎𝑚è𝑡𝑟𝑒: 𝑉 = , 𝑒𝑠𝑡 𝑎𝑝𝑝𝑒𝑙é 𝑣𝑖𝑡è𝑠𝑠𝑒 𝑑′ 𝑜𝑛𝑑𝑒 é𝑙𝑎𝑠𝑡𝑖𝑞𝑢𝑒
𝜌
L’équation (3) est une forme d’équation de la physique mathématique

connue sous le nom équation d’onde à une seule dimension:
La même analyse élargie par analogie au cas 3D nous permet

l’établissement des équations d’ondes pour les deux types de
déformations, nous aurons alors:

Onde de compression dilatation:
C’est la cas de l’onde P ayant pour équation
𝜕2𝜽 𝜕2𝜽 𝜕2𝜽 𝜌 𝜕2𝜽

2
+ 2+ 2=
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝝀 + 2µ 𝜕𝑡 2
𝜕2𝜽 𝜕2𝜽 𝜕2𝜽 1 𝜕2𝜽

2
+ 2 + 2 = 2 2 (4)
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝑽𝑷 𝜕𝑡
𝝀+𝟐µ
𝑽𝑷 =
𝝆
𝜽: Est appelé potentiel de déformation totale de compression dilatation

Onde de cisaillement:
C’est la cas de l’onde S ayant pour équation
𝜕 2 𝝋𝒙 𝜕 2 𝝋𝒙 𝜕 2 𝝋𝒙 𝜌 𝜕 2 𝝋
2
+ 2
+ 2
=
𝜕𝑥 𝜕𝑦 𝜕𝑧 µ 𝜕𝑡 2
𝜕 2 𝝋𝒙 𝜕 2 𝝋𝒙 𝜕 2 𝝋𝒙 1 𝜕2𝝋
2
+ 2
+ 2
= 2 2 (5)
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝑽𝑺 𝜕𝑡
µ
𝑽𝑺 =
𝝆
𝝋: Est appelé potentiel de rotation / cisaillement

CN°1

Transféré par

Droits d'auteur :

Formats disponibles

CN°1

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

CN°1

Transféré par

Droits d'auteur :

Formats disponibles

Sismique L3

Cours de sismique destiné au troisième

SAMAI S. [email protected] USTHB 2022-2023

C’est une notion fondamentale en géophysique. Il est défini comme une

SAMAI S [email protected] USTHB 2022-2023

Schéma d’un milieu continu dans l’espace

SAMAI S [email protected] USTHB 2022-2023

- La Terre entière renfermée dans sa

- Les constituants de la Terre séparés

- Structures géologiques (bassin sédimentaire,

- Modèle géologique réduit, Echantillon de

SAMAI S [email protected] USTHB 2022-2023

SAMAI S [email protected] USTHB 2022-2023

SAMAI S [email protected] USTHB 2022-2023

SAMAI S [email protected] USTHB 2022-2023

SAMAI S [email protected] USTHB 2022-2023

𝑉 = 𝑑𝑥𝑑𝑦𝑑𝑧 (Cube à l’intérieur)

SAMAI S [email protected] USTHB 2022-2023

𝑑𝑥𝑑𝑦𝑑𝑧 1 + ε𝑥 )(1+ ε𝑦 1 + ε𝑧 − 𝑑𝑥𝑑𝑦𝑑𝑧

SAMAI S [email protected] USTHB 2022-2023

Nous aurons pour ce type de déformation les expressions suivantes:

SAMAI S [email protected] USTHB 2022-2023

SAMAI S [email protected] USTHB 2022-2023

Le cisaillement est associé à une rotation que nous pourrons tirer

SAMAI S [email protected] USTHB 2022-2023

SAMAI S [email protected] USTHB 2022-2023

Analyse élargie au 3D:

Dans un système de coordonnées de 3D, XYZ (cartésien) la loi de Hook, pour un

SAMAI S [email protected] USTHB 2022-2023

SAMAI S [email protected] USTHB 2022-2023

SAMAI S [email protected] USTHB 2022-2023

Tenseur de contrainte tenseur de déformation

SAMAI S [email protected] USTHB 2022-2023

D’après la formule de Poisson, la dilatation dans l’axe XX induit une

La dilatation selon XX sera alors la somme de la composante

SAMAI S [email protected] USTHB 2022-2023

𝝈𝒙𝒙 = 𝝈𝒚𝒚 = 𝝈𝒛𝒛 = −𝒑

SAMAI S [email protected] USTHB 2022-2023

La somme des trois équations du système membre à membre

𝐸(ε𝑥x + ε𝑦y + ε𝑧z ) = (1 − 2𝝂)(𝝈𝒙𝒙 + 𝝈𝒚𝒚 + 𝝈𝒛𝒛 )

Relation entre le module de rigidité et le coefficient de Poisson:

Conformément à cette géométrie on peut écrire que :

SAMAI S [email protected] USTHB 2022-2023

SAMAI S [email protected] USTHB 2022-2023

L’allongement relatif selon la diagonal correspondra au

SAMAI S [email protected] USTHB 2022-2023

La variation de contrainte ∆p fluctue autour d’une valeur

La variation ∆p sera alors définie par :

SAMAI S [email protected] USTHB 2022-2023

Le deuxième coefficient de Lamé 𝝀:

𝑬𝜺𝒙𝒙 = (𝟏 + 𝝂)𝝈𝒙𝒙 − 𝝂𝝈𝒙𝒙 − 𝝂𝝈𝒚𝒚 − 𝝂𝝈𝒛𝒛

Qui s’écrit aussi:

SAMAI S [email protected] USTHB 2022-2023

𝐸(ε𝑥x + ε𝑦y + ε𝑧z ) = (1 − 2𝝂)(𝝈𝒙𝒙 + 𝝈𝒚𝒚 + 𝝈𝒛𝒛 )

SAMAI S [email protected] USTHB 2022-2023

L’état de contrainte dans ce milieu dépendra de deux paramètres: la

SAMAI S [email protected] USTHB 2022-2023

Par l’application de la loi de Hook linéaire:

SAMAI S [email protected] USTHB 2022-2023

L’équation (3) est une forme d’équation de la physique mathématique