Université des Sciences et de la Technologie Houari Boumediene

F S T G A T- Faculté de Géologie. L2-A (2023/2024)

Prospection Sismique

I) INTRODUCTION:
La prospection sismique est une méthode géophysique qui permet de
visualiser les structures géologiques en profondeur grâce à l'analyse des échos
d'ondes sismiques. La sismique a été mise en place grâce à la sismologie, qui
est l'étude des ondes sismiques générées lors de la survenance des séismes.

Les ondes sismiques étudiées peuvent avoir des causes naturelles

(séismes) ou artificielles (Camion vibreur, explosif, canon à air, etc.). Dans tous
les cas, elles suivent les mêmes lois de propagation que les ondes lumineuses
(lois de Snell et Descartes). Le comportement des sous-sols prospectés
repondent aussi aux criteres de la mécanique des milieux élastiques
(propagation d'ondes dans les milieux élastiques).

C’est la prospection géophysique la plus importante vu le nombre de

géophysiciens et géologues qui lui font appel. De grandes entreprises
mondiales (CGG, BGP, Western-Geco, Schlumberger, Total, Sonatrach, etc...)
octroient des budgets important pour améliorer cette prospection qui est
essentiellement utilisée dans l’exploration pétrolière, génie civil…etc.

II) Etapes de la prospection Sismique:

II-1) Acquisition des données sismiques:

 Sismique à terre: on utilise des camions vibreurs ou de la dynamite (Fig.

01) pour créer un front d'ode qui va se propager en profondeur. Le signal
émis va être réfléchie/réfracté et capté par un ensemble de géophones
placés en surface selon un profil rectiligne à des distances égales (Fig.
01). Tout l'ensemble est raccordé à une station de calcul et de traitement
préliminaire sur le terrain.

__________________________________________________________________________________
Prospection sismique 1
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

Géophone

Ondes réfléchies

Fig. 01: Compagne de sismique en terre

 Sismique en mer: on procède de la même manière avec un bateau

possédant des canons à air comprimé, créant ainsi des ondes qui vont se
propager à travers la couche d'eau et les fonds marins. Ces ondes vont
êtres enregistrées par hydrophones (Fig. 02). Les données de la sismique
en mer vont subir le même traitement que ceux à terre.

Fig. 02: Compagne de sismique en mer

__________________________________________________________________________________
Prospection sismique 2
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

 Sismique de puits (VSP): La nécessité de mieux contraindre les réservoirs

pour les exploiter de façon très optimale a conduit les géophysiciens à
développer la technique de sismique de puits à plus haute résolution. On
intervient dans un puits avec des géophones placés le long de la paroi du
puits et des sources en surface (profil sismique vertical: Vertical Seismic
Profil), ou vis versa (Fig. 03). Le VSP vient pour réconforter les données de
sismique en surface.

Fig. 03: Compagne de sismique de puits

II-2) Traitement des données: Apres la récolte des données vient l'étape
de traitement où les données sismiques subissent des traitements
informatiques très complexes qui permettront d'obtenir une image très nette
du sous-sol prospecté. Les ondes émises se propagent suivant les lois de
réflexion et de réfraction de Snell (appelées aussi lois de Descartes) et sont en
partie réfléchies à chaque changement de facies géologique. Ce changement de
facies est détecté par un changement de vitesse de propagation.

II-3) Interprétation des données sismique : La vitesse serait le paramètre

physique calculé en fonction de la distance entre la source et le géophone
(mesure de distance et mesure de temps). Ces vitesses de propagation des
ondes sismiques obtenues seront l'outil de l'interprétation. En effet, les vitesses
sont en fonction de l'état et de la densité de la matière. Certains types d'ondes
se propagent autant dans les liquides, les solides et les gaz, alors que d'autres
types ne se propagent que dans les solides.

__________________________________________________________________________________
Prospection sismique 3
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

On admet la présence en profondeur d'une discontinuité importante

séparant deux milieux dont les propriétés sont très différentes (pas la même
densité). Au niveau de cette discontinuité, les ondes sismiques doivent subir
une réflexion/réfraction importante. Cette déviation des ondes sismiques à la
discontinuité implique une variation de vitesse de propagation des ondes (Fig.
04).

O: Point d'ébranlement
AA', BB', CC': Fronts d'ondes
R: Rayon sismique
Demi-hyperbole en rouge: 1ere discontinuité
Demi-hyperbole en bleu: 2eme discontinuité

Fig. 04: Propagation de fronts d'ondes

L'interprétation des données sismique est régie par les formules de

propagation des ondes élastiques (mécanique des milieux continus). Suite à un
ébranlement provoqué ou naturel, des ondes élastiques prennent naissance à
partir de la source et se propagent dans le sol avec des vitesses sismiques
propre à elles. Les ondes sismiques sont très rapides, avec des vitesses de
l'ordre de quelques kilomètres par seconde. Ces vitesses dépendent du
matériau traversé, de la pression et de la température (Tableau 1).

III) Propagation d'une onde élastique : Dans un milieu homogène, isotrope

deux sortes d'ondes peuvent principalement se propager:

III-1) Ondes de volume : Les deux premières ondes qui arriveront aux
géophones sont les ondes de volume. Ce sont des ondes qui vont déformer le
volume excité.
__________________________________________________________________________________
Prospection sismique 4
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

La 1ère onde qui va arrivée est appelée onde P, onde longitudinale (Fig.
05). C'est une onde qui va déformer le volume sans le cisailler car sa direction
de propagation est parallèle à la direction de la trajectoire de l’onde. Elle
provoque un changement de volume (en le dilatant et en le comprimant).

Fig. 05: Onde première P (Onde longitudinale). La direction d'oscillation est //

à la direction de propagation

La 2éme onde qui arrivera aux géophones sera l'onde S, ondes transversale
(Fig. 06). C'est une onde qui va déformer le volume qu'elle traverse en le
cisaillant car sa direction de propagation est perpendiculaire à la direction de la
trajectoire de l’onde. Elle provoque un changement de forme en le cisaillant.

Fig. 06: Onde secondaire S (Onde transversale). Sa direction d'oscillation est 

à la direction de propagation.

Remarques: Seules ces deux ondes feront l'objet de l'interprétation des

données sismiques.

__________________________________________________________________________________
Prospection sismique 5
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

Nous donnons les formules avec lesquelles on calcule les deux vitesses P
et S :

√
Vp= 2 µ+ ¿ ¿---------------------------(01)
❑

Vs=
√ µ --------------------------------(02)
❑

Avec  et  (les coefficients de Lamé) sont des paramètres caractérisant

les propriétés élastiques du milieu traversé par ces ondes. Ce milieu peut être
solide ou liquide.

 défini le changement de forme du volume (dilatation ou compression).

Il est défini par la formule suivante:

.E
¿
¿¿
------------(03)

 défini le changement de forme par le cisaillement. Il est donné par la

formule suivante:

E
¿
2¿ ¿
------------------(04)

Avec E est le module d'élasticité (module de Young) défini par le rapport

entre la contrainte et la déformation.

 est le coefficient de Poisson, il défini la déformation transversale et

longitudinale.

Dans les solides, μ et λ sont très voisins mais dans les fluides, il n’y a pas
de cisaillement (μ=0).

Les ondes P se propagent dans les solides, les liquides et les gaz, alors
que les ondes S ne se propagent que dans les solides. On sait aussi que la
vitesse de propagation des ondes sismiques est proportionnelle à la densité du
matériel dans lequel elles se propagent. L'augmentation progressive de la
vitesse des ondes P et S indique une augmentation de densité du matériel à

__________________________________________________________________________________
Prospection sismique 6
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)
mesure qu'on s'enfonce en profondeur (Tab. 1). La chute subite de la vitesse
des ondes P est reliée au changement d'état de la matière (de solide à liquide).

La brusque interruption de propagation des ondes S indique la présence

d'un matériau dans un état liquide. Les ondes S, ondes transversales, ne se
propagent pas dans les milieux liquides car le coefficient de cisaillement  est
nul dans les liquides (Tab. 1). Leurs disparitions, indiquent qu'au niveau de
cette discontinuité, on passe d'un milieu solide à un milieu liquide. Ainsi, grâce
à la sismologie, on révèle l'existence de deux enveloppes terrestres très
différentes (Tab. 2).

Tableau 1: Ordre de grandeur des vitesses P et S pour certaines roches

__________________________________________________________________________________
Prospection sismique 7
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

Tableau 2: Ordre de grandeur des vitesses P et S à l’intérieur de la terre

III-2) Ondes de surface : Ces ondes sont appelés ondes de surface car
elles se propagent le long de la surface de la terre sans atteindre les grandes
profondeurs. Elles sont générées par l’arrivée des ondes de volumes à la
surface du globe suite à une très énergétique source comme dans le cas d'un
séisme ou bien une explosion nucléaire. Dans le cas du séisme, plus celui-ci
est fort plus les ondes de surface sont importantes, et elles peuvent faire
plusieurs fois le tour de la terre. On distingue deux types d'ondes:

Onde de Love (L): Les ondes de Love provoquent un ébranlement

horizontal comme les ondes S, mais elles ne se propagent que dans les solides
hétérogènes. (Fig. 07).

Fig. 07 : Propagation des ondes de Love (L)

Onde de Rayleigh (R) : Se déplaçant dans le plan vertical avec un

mouvement elliptique. Elle possède donc à la fois une composante
longitudinale et une composante transversale, assimilable à la houle de mer
(Fig. 08).

__________________________________________________________________________________
Prospection sismique 8
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

Fig. 08: Propagation des ondes de Rayleigh (R)

IV) Principe de propagation des ondes:

Les méthodes sismiques utilisent un champ élastique. A la suite d'une

excitation provoquée, des ondes élastiques prennent naissance et se propagent
dans le sol avec des vitesses sismiques propre à elle. Le sous-sol est composé
de couches géologiques différentes Chaque changement de couche est défini
par une interface. Le front d'onde sera ou réfléchit ou transmis/réfracté par les
différentes interfaces rencontrées (Fig. 09).

Fig. 09: Schéma représentatif des rais sismique réfléchis et transmis

IV-1) Onde directe: Ce sont des ondes longitudinales qui vont se propager
dans le premier milieu sans aucune réflexion ni réfraction. Elles vont être
enregistrées en premier par les géophones proches du point de tir. Leur
hodochrone est de la forme f(x) = ax (une droite qui passera par l'origine). La
pente de cette droite donnera l'inverse de la vitesse de la première couche (1/
V 1).

__________________________________________________________________________________
Prospection sismique 9
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)
IV-2) Onde réfléchie: Lorsque une onde sismique atteint, sous un angle
d’incidence (i 1, i 2, i3, i '1, i '2 , i '3) la surface de séparation entre la formation de
vitesse V1 et épaisseur h1, elle sera réfléchie par le réflecteur et renvoyée selon
une direction symétrique par rapport à la normale du point d'incidence (avec le
même angle d'incidence (i 1 ¿ . Mais une bonne partie de son énergie sera
transmise au deuxième milieu de vitesse V2 et d'épaisseur h2 avec un angle i 2
supérieur à i1 On parle alors d'onde transmise (Fig. 09). Cette dernière se
propage dans le second milieu et donne à son tour une onde réfléchie et une
onde réfractée lorsqu’elle rencontre la surface de séparation entre les milieux.
Ce processus est répété jusqu’à dissipation complète de l’énergie sismique.

Pour une formation à deux couches, la loi de Snell (ou de Descartes)

s'écrit par rapport au premier réflecteur comme suit:

sin i1 sin i2
V1
=
V2
----------------(05),

avec la condition que V 1< V 2

La loi s'écrit par rapport au 3eme réflecteur comme suit:

sin i1 sin i2 sin i 3

V1
=
V2
=
V3
------------------(06),

avec la condition que V 1<V 2<V 3

L’onde réfléchie de la 1ere interface est enregistrée par un récepteur

(géophone G2 sur Fig. 09) situé à une distance x de la source. Le temps de
parcours aller-retour entre la source et le récepteur est donné par la formule:

1
t=
V1
√ x 2 + 4 h2----------------------- (07)

__________________________________________________________________________________
Prospection sismique 10
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

Cette courbe qui décrit la relation entre le temps de parcours de l'onde

et la distance source-récepteur est connue sous le nom d'HODOCHRONE. La
caractéristique de l'hodochrone d'une onde réfléchie est une hyperbole de la
forme f(x) = ax² +b (Fig. 10).

L'analyse de l'hodochrone nous fournis des informations sur la vitesse de

la couche et son épaisseur. Seul le premier terme de l'équation dépend de la
distance x. Le second terme ne dépend que des paramètres de la couche. Ce
2eme terme n'est que le temps aller-retour de la source à la base de la couche.
Ce terme va nous donner l'épaisseur de la 1ere couche:

2h
Pour x=0 dans ( 01 ) on aura t 0=
V1

Finalement, tous les géophones qui sont proche du point de tir vont
capter en premier les ondes directes, elles vont être suivies par les ondes
réfléchies (Fig. 10).

Fig. 10:Hodochrones des ondes réfléchies (hyperbole) et les ondes directes

(droite)

On peut généraliser la loi de Snell et Descartes pour un terrain composé

de n couches par la formule:

sin i1 sin i2 sin i 3 sin i n

= = =… −−−−(08)
V1 V2 V3 Vn

Avec V 1 <V 2 <V 3 <…< V n

__________________________________________________________________________________
Prospection sismique 11
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

Cependant, plus on s'éloignera de la source plus l'angle d'incidence

augmentera, ce qui va donner naissance aux ondes réfractées.

IV-3) Onde réfractée: Pour enregistrer les ondes réfractées, Il existe une
valeur critique de l'angle d'incidence " i 1" pour le quel "i 2"= 90°. Cette de l'angle
i 2 va permettre au rayon sismique de se propager dans le 2 ème milieu le long du
réflecteur de vitesse V2 (Fig. 11).

Fig. 11: Rai réfracté par la discontinuité à partir d'une valeur limite de l'angle
d'incidence

sin i1 sin i2
Pour V1 < V2, la loi de Snell sera : V1
=
V2

sin i1 sin 90 °
i 1 va atteindre une valeur critique i c i 2=90 ° = V
V1 2

V1
sin i c = −−−−−−( 8 )
V2

V1
i c ¿ arcsin −−−−−−(9)
V2

A) Cas d'une formation à 02 couches: Le trajet de l'onde réfractée (conique)

pour aller de la source jusqu'au géophone B, pour une formation à 2 couches
(Fig. 12), est donné par la formule suivante:

x 2 h. cos i c
t conique= + --------------- 10
V2 V1

__________________________________________________________________________________
Prospection sismique 12
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

V2

Fig. 12: Trajet de l'onde conique [ACDB] pour une formation à deux couches

L'équation (10) est du 1er degré de la forme f ( x )=a x+ b. La pente de cette

1
droite est V (Fig. 13).
2

L'intercepte t 0 (pour x = 0) de la droite (formule 10) permet de

calculer l’épaisseur h de la 1ere couche :

2 h. cos i c t 0 .V 1
t 0=
V1
h=
2cos ic
---------(11)

V1
avec i c =arcsin ( V )
2

Il existe une distance critique x c entre la source 'S' et les géophones à

partir de laquelle on commence à enregistrer la réfractée (Fig. 13). Cette
distance peut être déterminée graphiquement (projection sur l'axe des
distance) ou bien par la formue: x c =2 h . tan i c

__________________________________________________________________________________
Prospection sismique 13
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

À la distance x b, l'onde directe et l'onde réfractée arrivent au même

moment au même géophone : t directe =tr é fract é e . Cette distance peut être
déterminée graphiquement (Fig. 13), ou bien par la formule suivante:

2 h . √ V 2 ²−V 1 ²
x b= -------------(11)
V 2−V 1

Cependant, tous les géophones situés à gauche de ce point x b, le

temps du trajet de l'onde directe est inférieur au temps du trajet de la réfracté
et inversement pour les géophones situés à droite du point x b (Fig. 13).

Fig. 13:La droite de pente 1/V1 (Ondes directes tracées en rouge sur le
graphe). La droite de pente 1/V2 (ondes réfractées tracées en bleu sur le
graphe)

B) Cas d'une formation à 03 couches: De la même manière on peut calculer

le trajet des différents rais réfractés qui parviennent des différents réflecteurs,
en exprimant de proche en proche le temps de trajets réfractés totalement à la
surface des 2, 3, 4, ……et n couches. On obtient les équations de droites qui
correspondent à une formation à 03 couches tabulaires (Fig. 14) :

__________________________________________________________________________________
Prospection sismique 14
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)
Hodochrone de la réfractée de 2ème couche:
x 2 h2 cos ic 2 2 h1 cos i c1
t 3= + +
V3 V2 V1

Hodochrone de la réfractée de 3ème couche:

x 2 h3 cos i c3 2 h2 cos i c 2 2h 1 cos i c 1
t 4= + + +
V4 V3 V2 V1

La formule générale de l'hodochrone de la niemme-1 couche serait:

n −1
x 2 h cos i ci
t n= + ∑ i
V n n =1 Vi

avec la loi généralisée de Snell-Descartes :

sin i c1 sini c 2 Sini c3 1

V1
=
V2
= V3
=…
Vn

Connaissant les vitessesV 1, V 2, V 3, ….et V n et les angles critiques à chaque

interface (par leurs Sinus) on peut alors calculer les épaisseurs h1, h2, h3 , ………, h n
(Fig. 14).

V4
Fig. 14: Trajet des ondes coniques pour une formation à 3 couches

__________________________________________________________________________________
Prospection sismique 15
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)

Fig. 15: Hodochrone d'une formation à 03 couches

V) Conclusion:

La prospection sismique est essentiellement utilisée dans la recherche

pétrolière, plus particulièrement la sismique réflexion car elle permet
d'identifier les structures géologiques à grandes profondeurs. En Algérie, des
champs d'hydrocarbures ont été découverts essentiellement à l'aide de la
prospection sismique, tel que le champ pétrolier de Hassi-Messaoud, qui
contenaient des réserves estimées à plusieurs Milliards de Tonnes d'huiles, et
le champ gazier de Hassi R'mel avec des réserves approximatives de 2 Trillons
de mètres cubes de Gaz.

En géologie structurale, la sismique réfraction est très demandée

dans la recherche minière, en Génie Civil (étude pour les fondations des
bâtisses, optimisation d'implantation de structures).

__________________________________________________________________________________
Prospection sismique 16
 Université des Sciences et de la Technologie Houari Boumediene
F S T G A T- Faculté de Géologie. L2-A (2023/2024)
En hydrogéologie la sismique réfraction est largement utilisée pour
identifier les positions et les niveaux des nappes, les zones fracturées/altérés
qui peuvent servir de réservoirs.

Dans le milieu aquatique la sismique est utilisée pour le dragage,

ensouillage de câble, pipeline.

VI) References:

Benaissa Z. Cours de sismique, FSTGAT, Dpt de Géophysique

Dubois J. & Diament M. - Géophysique, Cours et exercices corrigés - Editions

Dunod

Lavergne M. - Méthodes Sismiques - Ed. Technip

Yilmaz, O. (2000). Seismic data Analysis. Society of Exploration Geophysicists

__________________________________________________________________________________
Prospection sismique 17