LA GEOPHYSIQUE POUR LES

GEOLOGUES
Tome 4
LES METHODES SISMIQUES

INTRODUCTION :
Cet ouvrage traitant la géophysique pour les géologues, est le premier du
genre, publie par l'O.P.U. Nous sommes tres heureux de le présenter
aujourd'hui ; longtemps souhaite et attendu de la part de notre communaute
universitaire: (Enseignants ,Etudiants et Chercheurs).Nous esperons a travers
ces ouvrages:Enrichir d'une part le rayon des Sciences de la Terre de notre
bibliothèque Nationale et d'autre part percer les Mysteres de l'Univers par la
connaissance de la Geophysique et de ses differentes methodes (deja publies
par l'O.P.U :La Magnétotellurique, Bases physiques des méthodes sismiques et De la Pesanteur à
la Gravimetrie).Les géologues trouveront dans cet ouvrage les bases theoriques
et pratiques de la géophysique. Il interessera également le grand public,
curieux de s'avoir ce qui se passe ou ce que cache le sous-sol; dans le
domaine des hydrocarbures, dans la recherche minière, en hydrogéologie,
dans les travaux publiques (Batiments, Ponts et Chausses) ; les Chemins de
fer,.etc...
Dans ce livre, nous presentrons les bases theoriques et pratiques de la
géophysique:
1- La géophysique de surface: Méthodes électriques et électromagnétiques ,
les méthodes gravimétriques et magnétiques, les méthodes radiométriques,
les méthodes sismiques et seismiques (seismologie).
2 - La géophysique des puits ou diagraphies.
 2

Tout d'abord quelques notions sur la théorie du potentiel (distribution des

champs: champ électrique, magnétique, électromagnétique, gravimétrique, et
radiométrique ). Nous présenterons ensuite les caracteristiques "géo-
physiques ou pétro-physiques" des roches: L'élasticite; les propriétés
électriques et électromagnétiques, les propriétés magnétiques et
gravimétriques ainsi que les propriétés radiométriques des roches.
Nous donnerons un aperçu sur les différentes techniques d'énrégistrement de
traitement et d'interprétation de ces paramètres.

TABLE DES MATIERES

CHAPITRE I INTRODUCTION A LA SISMIQUE

I-1 Elements de base de la sismique

I-1-1 La sismique reflexion
I-1-2 La sismique refraction
I-1-3 La sismique codee
I-1-3-1 Etude de quelques procedes en sismique codee
I-1-3-1-1 Procede vibroseis
I-1-3-1-2 Procede sosie
I-1-3-1-3 Procede mini-sosie

CHAPITRE II BASES PHYSIQUES DE LA PROSPECTION

SISMIQUE

II-1- Principe de la theorie de l'elasticite

II-1-1 Le module d'allongement (Module de Young E )
II-1-2- Le coefficient de contraction (Poisson )
II-1-3 Le module de regidite (Coulomb)
II-1-4 Le module d'incompressibilite (K)

2
 3

II-2 Ondes sismiques

II-2-1 Ondes longitudinales
II-2-2 Ondes transversales ou de cisaillement
II-2-3 Ondes superficielles (Ondes de Rayleigh
et Ondes de Love)

CHAPITRE III PARAMETRES CARACTERISANT UNE

OSCILLATION

III-1 Introduction
III-2 Oscillations harmoniques et periodiques
III-3 Superposition des oscillations
III-4 Decomposition des oscillations
III-5 Decomposition spectrale des oscillations
III-6 Representation des spectres
III-6-1 Proprietes de conformites reciproques
III-6-2 Addition des spectres
III-6-3 Similitude des oscillations et des spectres
III-6-4 Spectre d'une oscillation decalee en temps par
rapport a une autre oscillation
III-6-5 Spectre d'une oscillation serree ou comprimee
sur l'axe des temps
III-6-6 Spectres de certaines impulsions
III-6-6-1 Cas d'une oscillation sinusoidale
III-6-6-2 Cas d'une oscillation sinusoidale amortie
III-6-6-3 Cas d'une fonction rectangulaire
III-6-6-4 Cas d'une impulsion rectangulaire de duree t
III-6-6-5 Cas d'une Impulsion de Dirac
III-6-6-6 Cas d'une Impulsion de Berlape
III-6-6-7 Impulsion en Cloche
III-6-7 Decomposition par des Fonctions Rectangulaires
III-6-8 Convolution
III-6-8-1 Pratique de la Convolution

CHAPITRE IV PROPAGATION DES ONDES SISMIQUES

3
 4

IV -1 Introduction
IV-2 Principe de la Sismique Geometrique
IV-2-1 Principe de Huyguens
IV-2-2 Principe de Fermat
IV-2-3 Principe de Reciprocite
IV-2-4 Principe de Superposition
IV-3 Hodochrone Vitesse Apparente
IV-4 L'Equation de l'Onde et les Solutions Particulieres
IV-5 Reflexion et Refraction des Ondes Sismiques
IV-5-1 Generalites sur les Hodochrones,sur les Champs de
Temps des Ondes Reflechies et Refractees et Ondes
Coniques
IV-6 Impedance Accoustiques,Coefficient de Reflexion
et de Refraction
IV-6-1 Analyse de quelques cas Particuliers
IV-6-1-1 Incidence Normale
IV-6-2 Variation du Coefficient de Reflexion avec l'Angle
d'Incidence
IV-6-3 Relation entre le Coefficient de Reflexion
et de Transmission
IV-6-4 Coefficient de Reflexion et Proprietes Physiques des
Sediments
IV-6-4-1 Coefficient de Reflexion et Diametre des Grains
IV-6-4-2 Coefficient de Reflexion et Porosite
IV-6-4-3 Coefficient de Reflexion et Frequence
IV-6-4-4 Mesure du Coefficient de Reflexion
IV-6-4-4-1 Incidence Sur la Surface Libre (surface de la terre)
IV-6-4-4-2 Incidence sur la Surface du Milieu de Regidite
Absolue
IV-6-4-4-3 L'Angle d'Incidence est Petit
IV-6-5 Reflexion et Refraction des Ondes dans
les Milieux a Plusieurs Couches
IV-7 Les Amplitudes en Sismique Reflexion
IV-7-1 Etudes des Amplitudes
IV-7-2 Les Facteurs Dependant de la Subsurface
IV-7-3 Les Facteurs Lies a la Reflectivite des Miroirs

4
 5

Sismiques
IV-7-4 Les Differents Facteurs qui Affectent les
Amplitudes de Reflexion
IV-8 Les Effets de la Nappe Sismographiques
IV-9 Les Effets de la Transmission
IV-9-1 Distorsion de l'Ondelette
IV-10 Principe de la Deconvolution Multicanal
IV-10-1 Coefficient de Reflexion
(Calcules aux divers Interfaces)
IV-10-2 Formule de Wyllie
IV-11 Influence de la Zone Alteree
IV-11-1 Correction d'Amplitude (Surface Consistante)
IV-11-2 Amortissement des Ondes Sismiques
IV-12 Diffraction des Ondes
IV-13 Ondes Utiles et Parasites

CHAPITRE V ACQUISITION DES DONNEES SISMIQUES

V-1 Introduction
V-2 Sismographe ou Geophone
V-2-1 Introduction
V-2-2 Etude d'un Sismographe(Sismographe Electromagnetique
a Flux Coupe)
V-2-3 Description d'un Sismographe
V-2-4 Multiplication des Aismographes ou de Geophones
V-3 Aspect Aleatoire des Sismographes
V-4 Classification des Ondes Parasites
V-4-1 Les Parasites Horizontaux
V-4-1-1 Les Ondes de Rayleigh
V-4-1-2 Les Ondes de Love
V-4-1-3 Les Ondes Fantomes
V-4-1-4 Les Ondes Guidees
V-4-1-5 Les Bruits Naturels
V-4-1-6 Les Bruits Provoques par le Tir
V-4-1-6-1 L'Onde Directe
V-4-1-6-2 L'Onde Refractee

5
 6

V-4-1-6-3 L'Onde Aerienn

V-4-1-6-4 Les Reflexions Multiples
V-4-1-6-5 Les Ondes Diffractees
V-4-1-6-6 Les Oscillations Parasites

CHAPITRE VI PARAMETRES GEO - SISMIQUES

VI - 1 Introduction
VI-1-1 Notions sur les Vitesses Aismiques
VI-1-2 Expression Mathematique des Vitesses des Ondes
Elastiques Dans les Milieux Solides non Poreux
VI-1-3 Les Vitesses Rencontrees en Sismique
VI-1-3-1 Vitesse Instantanee
VI-1-3-2 Vitesse Moyenne
VI-1-3-3 Vitesse d'Intervalle
VI-1-3-4 Vitesse quadratique
VI-1-3-5 Vitesse par Tranche
VI-1-4 Mesure des Differents Types de Vitesses
VI-1-4-1 La Vitesse Instantanee
VI-1-4-2 La Vitesse Moyenne
VI-1-4-3 La Vitesse d'Intervalles
VI-1-4-3-1 Determination Indirecte
VI-1-4-4 Autres Moyens
VI-1-4-4-1 Interets des Vitesses d'Intervalles
VI-1-4-5 La Vitesse Quadratique Moyenne
VI-1-4-6 La Vitesse de Sommation
VI-2 rincipaux Facteurs Modifiant la Vitesse dans les
differentes Roches
VI-2-1 Introduction
VI-2-2 La Lithologie
VI-2-3 La Nature de la Roche
VI-2-4 Age de la Roche
VI-2-5 Pression - Profondeur
VI-2-6 Anisotropie
VI-2-7 La Porosite
VI-2-7-1 La Porosite Primaire

6
 7

VI 2-7-2 La Porosite Secondaire

VI-2-7-3 Determination de la Porosite a Partir des Diagraphies
VI-2-7-3-1 Diagraphie Sonique
VI-2-7-3-2 Diagraphie Neutronique(Gamma-Gamma
pour un Sable a Gaz)
VI-2-7-3-3 Diagraphie de Densite (Gamma-Gamma
pour un Sable a Gaz)
VI-2-7-3-4 Diagraphie Electrique
VI-2-7-4 Porosite dans les Formations (Consolidees et non
Consolidees)
VI-2-7-5 Role de la Densite
VI-3 Zone Alteree (WZ)
CHAPITRE VII REPRESENTATION GRAPHIQUE
DES DIFFERENTES ONDES

VII-1 Hodochrones des Differentes Ondes

VII-2 Champ des Temps et Indicatrices des Ondes Directes
VII-2-1 La Source se Trouve a la Surface de la Terre
VII-2-2 L'Indicatrice est Lineaire
VII-2-3 Indicatrices des Ondes Reflechies
VII-2-3-1 Indicatrice Pour un Reflecteur Plan
VII-2-3-2 Indicatrice Longitudinale
VII-2-3-3 Indicatrice non Longitudinale
VII-2-3-4 Indicatrice des Ondes Reflechies
VII-2-3-5 Indicatrice des Ondes Diffractees
VII-2-3-6 Indicatrice des Ondes Multiples
VII-2-3-6-1 Les Multiples de Surface
VII-2-3-6-2 Les Multiples Internes
VII-2-3-7 Dromochronique des Ondes Refractees
VII-2-3-7-1
VII-3 Hodochrones des Ondes Sismiques
(Dans le Milieu a Plusieurs Couches)
VII-3-1 Indicatrices des Ondes Reflechies
VII-3-2 Rapport entre les Hodochrones des Ondes Differentes
VII-3-2-1 Rapport entre les Indicatrices des Ondes Directes et
reflechies

7
 8

VII-3-2-2 Rapport entre les Hodochrones des Ondes

Reflechies et Refractees

CHAPITRE VIII FILTRAGE SPATIAL

VIII-1 Introduction
VIII-2 Qu'appelle-t-on Vitesse Apparente
VIII-2-1 Rapport entre la Vitesse Apparente et la Vitesse Reelle
de Propagation en Fonction de l'Angle d'imergence
VIII-3 Definition des Bruits
VIII-3-1 Bruits Aleatoirs
VIII-3-2 Bruits Semi-Aleatoirs
VIII-3-3 Bruits de Diaphonie
VIII-4 Filtrage Spatial ou ( en Nombre d'Onde , K)
VIII-4-1 Calcul de la Fonction de Transfert
VIII-4-2 Etude de la Fonction de Transfert
VIII-4-3 Calcul et Representation Graphique des Courbes
d'Attenuation.

8
 9

CHAPITRE I INTRODUCTION A LA SISMIQUE

I - 1 - Elements de Base de la Sismique:

La sismique ou seismique s'emploie pour designer d'une maniere generale les

vibrations et les secousses et d'une maniere particuliere aux tremblements de terre
ou seismes. C'est une methode geophysique ,sa theorie est construite a partir des
connaissances de l'elasticite ,de l'optique geometrique,de la physique vibratoire ,
du traitement du signal , et enfin de la geologie .
Pour connaître la constitution géologique d'une region,les géophysiciens
procèdent a diverses mesures de surface. Pour préciser en particulier les courbes
de niveau des couches geologiques,ils emploient les methodes de la prospection
sismique:Une explosion est provoquee au voisinage des couches superficielles du
terrain etudie.L'onde de choc,provoque artificiellement se propage dans le sous-
sol.
Les couches geologiques etant de densite et de nature differentes,lorsque le front
d'onde franchit la frontiere separant deux couches,une partie de l'energie
transportee est :reflechie et refractee vers la surface du fait de la discontinuite des
constantes elastiques des couches.
En surface des sismographes enregistrent les arrivees de l'energie liberee par
l'explosion durant les 4 ou 5 secondes qui la suivent.
En pointant ces arrivees sur les sismogrammes ( section temps) et en recoupant
leurs indications avec celles procurees par d'autres mesures,le sismicien precise
sa connaissance du sous-sol:position des discontinuites,vitesse de propagation ,
parfois coefficient d'absorption

9
 10

On constate depuis quelques annees que les methodes sismiques sont de plus en
plus utilisees a des fins de reconnaissance.Il ya a cette evolution des raisons
techniques lies aux conditions geologiques et aussi des raisons economiques.
En Algerie,comme on le sait,la province triasique est l'une des regions les plus
riches en hydrocarbures du Sahara Algerien.
Longtemps consideree comme presentant peu de perspectives,ce n'est qu'en 1952
qu'on a entrepris des travaux geologiques et geophysiques intensifs ayant pour but
la recherche d'accumulations commerciales d'hydrocarbures.
En 1956,des champs geants y ont ete decouverts essentiellement a l'aide de la
jprospection sismique;par exemple :Hassi-Messaoud,avec des reserves en place
estimees a 5 Milliards de Tonnes d'Huiles et Hassi-R'mel avec des reserves
approximatives de 2 Trillons de metres cubes de Gaz et 500 Millions de
Condensat.
Donc,on peut constater que le role de la prospection sismique dans la recherche
petroliere est tres grand.
On distingue deux methodes de prospection sismique:la sismique reflexion et la
sismique refraction.Elles ont toutes deux pour but de determiner la profondeur et la
forme des discontinuites geologiques,qui constituent le sous-sol.

I - 1 - 1 La Sismique Reflexion:

Largement employee pour la prospection industrielle,est responsable de la

decouverte de tres nombreux gisements de petrole.
Le principe de cette methode est simple:des ondes sismiques
(elastiques)engendrees par l'explosion de charges de quelques Kg d'explosion se
propagent a travers les terrains sedimentaires ou metamorphiques,se reflechissent
sur les surfaces de separation des formations de nature differentes.
Les ondes reflechies sont enregistrees a la surface de la terre par une serie de
sismographes etales sur une surface appelee Nappe de geophones ou trace (le
plus souvent:24,48,96,192 Traces).
Si l'on connait la vitesse de propagation des ondes sismiques dans les differents
niveaux,la mesure des temps de propagation permet de calculer la profondeur de
ces niveaux.

I - 1 - 2 La Sismique Refraction:

10
 11

La sismique refraction n'est pas en mesure de detailler comme la methode par

reflexion,la morphologie de toutes les formations profondes.Elle permet,en
revanche,de determiner les vitesses de propagation des ondes sismiques,ainsi que
la profondeur des differentes interfaces.L'onde sismique se refracte sur l'interface
et revient aux sismographes qui sont situes a la surface du sol.Le temps de
parcours de l'onde refractee permet de determiner la profondeur des
interfaces.Notons que,loin de ce concurrencer ces deux methodes se completent et
il est souhaitable de faire suivre une etude par refraction,que l'on doit considerer
comme une etude de reconnaissance,d'une etude par reflexion. Au cours de son
application la prospection sismique a connu des mutations,par exemple
l'enregistrement photographique a ete remplacee par l'enregistrement analogique
sur les bandes magnetiques et a l'heure actuelle on l'a abondonnee totalement au
profit de l'enregistrement digital. L'Acquisition et le traitement des donnees
sismiques sont realises par ordinateur en temps reel ou legerement differe.

I-1-3 La sismique codee

La sismique codee selon (T.Bernhardt et J.H Peacock) est redevable des

techniques de codage developpees dans la theorie de la communication.En
effet,dans cette theorie,la transmission des informations se fait selon un code etabli
auparavant;code qui permet a un fil unique de transmettre pendant un temps
donne,plusieurs informations a differents recepteurs. Pour en venir a la sismique
codee,on sait que l'energie necessaire pour obtenir des informations sur le sous-
sol peut etre injectee en une emission faible de longue duree (7 secondes).
Contrairement a la sismique conventionnelle qui est caracterisee par l'injection
d'une seule emission puissante (exemple: l'explosion de dynamite). Le signal emis en
sismique "codee" est fonction du procede employe: Il est bref et discontinu dans le
cas du "sosie" et du "mini-sosie";tandis qu'il est continue et variable en frequence
dans le temps,dans le cas du "vibroseis".
L'emission represente un ensemble de sequences codees determinees par un
calculateur electronique.La sequence est une succession d'instant d'emission ou
une succession d'equidistqnce entre sources d'emission dans le cas des procedes
utilisant plusieurs sources.
Si nous prenons le cas de la sequence qui est constituee d'une succession
d'instant d'emission;on peut remarquer que les intervalles de temps entre
emissions successives sont au plus egaux au plus long temps aller et retour de

11
 12

l'onde.L'enregistrement des Informations se fait dans le domaine temporel ou dans

le domaine frequentiel.En sismique conventionnelle,au cours de l'enregistrement
des informations,il est interdit de declencher une autre impulsion.
Les enregistrements obtenus en sismique codee sont illisibles;en effet nous
enregistrons un melange d'informations du aux differents miroirs rencontres par
l'onde au cours de sa propagation;qu'ils soient obtenus en vibrosismique
(vibroseis,seiscode,mini-sosie) ou en sismique codee impulsionnelle (sosie).

I-1-3-1 But de la sismique codee

Les principaux buts de la sismique codee sont:L'obtension d'une quantite

importante d'informations,l'amelioration du rapport su bruit,la reduction du cout des
operations tant a terre qu'en mer.

I-1-3-2 Principe de base de la sismique codee

Un emission sismique donne naissance a un signal S(t) qui se propage dans le sol
ou dans l'eau (sismique marine),se reflechi ou se refracte;on enregistre S(t)
convolue avec les differents coefficients que l'onde emise rencontre au cours de sa
propagation;ce principe qui vient d'etre enonce caracterise bien la sismique
conventionnelle;par contre en sismique codee on determine auparavant ou au
cours de l'emission une sequence codee y(t).
S(t) etant le signal caracterisant la source utilisee;ce signal se propage et arrive
aux hydrophones ou geophones,on enregistre alors un signal f(t) qui est la
convolution de S(t) par y(t),c'est-a-dire:

f(t) = S(t) * y(t) * RC log (t)

ou:
* est le signe representant la convolution
y(t) la sequence codee
S(t) etant le signal caracterisant la source utilisee
f(t) est la convolution de S(t) par y(t)
RClog(t) parametres caracterisant la reflexion ou la refraction.

I-1-3-3 Enregistrement:

12
 13

L'equipement utilise pour les enregistrements terrains differe peu selon les
procedes utilises.
On a tout d'abord des recepteurs (hydrophones ou geophones),des amplificateurs
qui amplifient le signal recu par les recepteurs,d'un calculateur qui doit realiser
l'operation de crosscorrelation entre la sequence d'emission et de la trace recue.Le
calculateur numerise toute les traces qui y sont introduite (code binaire),indiquant
ainsi la presence ou l'absence d'une emission acoustique au temps
correspondant;la visualisation est imediate.

I-1-4 Etude de quelques procedes en sismique codee

I-1-4-1 Le procede Vibroseis

Le vibroseis est un des procedes les plus anciens de la prospection sismique.Ce

procede utilise comme source d'emission un signal de reference genere par des
vibrateurs (La source d'emission est constitue d'une masse de quelques
tonnes,actionnee par un vibrateur hydraulique generant ainsi une onde sinusoidale
de frequence variant de facon monotone dans le temps).Le signal sinusoidal est
"choisi" en fonction de la bande de ferquence sismique utilisee.

"Le choix du signal de reference dans le procede Vibroseis a fait l'objet de plusieurs
recherches;ainsi plusieurs auteurs ont contribuer a leur determination;parmi eux on peut
citer:Klander et All (1960);Cook et Paolillo (1964);Nigel Anstey (1964);Robert Geyer (1972).
La gamme de frequence et la duree du signal doivent etre choisies de telle maniere que la fonction
d'autocorrelation presente un pic central bien developpe et plus etroit par rapport aux pics
secondaires qui l'entourent.Il faut aussi noter qu'une attenuation de l'amplitude des pics secondaires
peut etre obtenue a la suite d'une addition successive de signaux appartenant a la meme bande de
frequence,emis en phase."

I-1-4-2 Le procede Sosie

La technique sosie (technique de sismique codee impulsionnelle),mise au point par la

SNPA;est une technique basee sur l'utilisation des sources non explosives
hautement repetitives:telles que:Etincelleur;canon a gaz,a air ou a vapeur,utilisee
dans la prospection maritime.

13
 14

L'utilisation du procede dans le cas des exploitations a faibles profondeur d'eau,les

petites sources sont plus pratiques a cause d'une emission riche en basses
frequences,une maniabilite facile et un couplage satisfaisant.
Le procede sosie est caracterise par l'energie declenchee dans l'eau,sous la forme
d'une impulsion de duree tres reduite,de meme gain et de meme amplitude;la
longueur de l'enregistrement traite;l'ordre de couverture;la distance entre les points
communs miroir ( CDP:Commun Depth Points).

I-1-4-3 Le procede Mini-sosie

Le mini-sosie est un procede de sismique codee,caracterise par l'utilisation d'une

dameuse comme source sismique et surtout par le changement du taux d'emission
par variation de vitesse de "frappe" de la source.
Le principe du procede est simple,la source "dameuse" ou "pilon"comme celui
utilise pour les travaux en genie-civil;est constitue a sa base par une plaque qui
frappe le sol,sur ce dernier est place un capteur qui transmet le temps exact de
l'emission de l'impulsion S(t),caracterisant la source.
Pour traiter les donnees enregistrees,on applique la crosscorrelation de la
sequence codee et du signal recu f(t).

F(t)=f(t)*y(t)=AFC y(t)*S(t)*RC log(t)

ou :
* est le signe representant la convolution
y(t) la sequence codee
S(t) etant le signal caracterisant la source utilisee
f(t) est la convolution de S(t) par y(t)
AFC est la fonction d'autocorrelation
RClog(t) parametres caracterisant la reflexion ou la refraction.
F(t) est l'information utile

CHAPITRE II - BASES PHYSIQUES DE LA PROSPECTION

SISMIQUE:

II - 1- Principe de la Theorie de l'Elasticite (Deformation)

14
 15

Il n'existe pas de corps indeformable.Un solide,du fait qu'il est soumis a des
tensions,se deforme.Il joue le role de ressort,il est elastique,au sens vulgaire du
terme.L'elasticite est une propriete des corps a reprendre leur forme initiale ou leur
volume apres les avoir perdus a la suite d'une compression ou d'une tension.
En theorie d'elasticite on admet aue forces et deformations sont infinement
petites.Cela revient a eriger le principe que les relations entre forces et
deformations sont reversibles c'est-a-dire qu'en supprimant les tensions on
supprime ipso-facto des deformations.
Cela a pour consequence que,si deux systemes quelconques de forces F,F1
engendrent respectivement deux systemes D,D1 de deformations,l'action
simultanee F+F1 des deux systemes de forces engendrent le systeme D+D1 de
deformations. Afin de mieux comprendre les methodes sismiques,il est necessaire
de rappeler ou de preciser quelques points essentiels sur la propagation des
ondes elastiques dans les solides homogenes et isotropes. Tout d'abord notons
que l'on considere approximativement le milieu geologique comme un corps
elastique.Considerons une force F appliquee sur une petite surface d'un materiel
homogene et isotrope suspendu:Trois sortes de deformations elastiques peuvent y
resulter:des deformqtions de compression,des deformation de tensions (on les
appelles encore des deformations de volume) et des deformations de cisaillement
(ou des deformations de forme).En effet les deformations de compression et de
tension sont de meme nature,seul le sens de la force appliquee
perpendiculairement a la surface ayant varie.Il s'agit,dans les deux cas,de
deformations longitudinales(de volume).Dans le cas d'une deformation de
cisaillement,la force est parallele a la surface.Il en resulte que nous n'avons plus
qu'une deformation angulaire (de forme).On peut definir les deformations d'un
solide isotrope a l'aide de deux parametres elastiques,choisis parmi les quatres
suivants:
II -1-1 Le Module d'Allongement (Module de Young E )

Dans le cas d'une tension ou d'une compression,l'experience montre que,si la

tension ou la compression est faible les deformations d'un corps sont directement
proportionnelles a la force appliquee F.Cette loi s'appelle la loi de HOOKE.
Si le corps a un diametre ''d'' et en tenant compte que l'allongement relatif est
proportionnel a la force appliquee,on peut ecrire:

15
 16

Dl / l = 1 .F / E.p.r2 = T / E = 1.4.F / E.p.d2

ou bien:

E = T / (Dl / l )

ou Dl / l : est l'allongement relatif du corps,'' l '' represente la longueur initiale du

materiel ( Dl ) l'allongement consecutif a l'application de la force ''F'' ou bien de la
tension '' T ''.

La loi de Hooke: En presence de faibles contraintes,il existe une proportionnalite

entre contraintes et deformation resultantes.Ainsi sous l'action d'une faible
contrainte un corps elastique change de dimensions et de forme.

II-1-2 Le Coefficient de Contraction ( Poisson )

Le schema ( ) montre que lorsqu'une tension est appliquee,on a une allongement

du corps,;ais en meme temps on a un amincissement dans la direction
perpendiculaire a celle de la tension.Dans le cas d'une compression ( fig: ),on
aura,au contraire une diminution de la longueur de l'element considere et un
epaississement suivant une direction perpendiculaire a celle de la compression.Le
coefficientit de contraction ou coefficient de Poisson,designe par le terme s est
egal a:

s=(Dd/d) / (Dl / l)

Ce coefficient caracterise le rapport entre la compression relave et l'allongement

relatif.Il ne depasse jamais 0.5 et il est de l'ordre de 0.25 pour la plupart des
roches.s est une grandeur sans unite.
Les coefficients de Poisson s et le module de Young E sont des caracteristiques
principales de l'elasticite des materiels isotropes.Ils ne dependent ni de la forme,ni
de la dimension des corps.
Ces coefficients sont mesures pour de nombreux materiaux,dont les valeurs sont
presentees ci-dessous:

16
 17

Materiel E.10 s Materiel E.10 s

Fer 19.6 0.27 Cal caire 2.5 - 6.0 0.22 - 0.35

Verre 6.0 0.27 Marnes 1.5 - 4.5 0.3 - 0.4
Granite 3.5 0.3-0.25 Gres 1.5 - 3.0 0.2 - 0.35

II-1-3 Le Module de Regidite ( de Coulomb)

Si une force F agit parallelement a la surface S de l'element solide considere ci-

dessus (fig ) la pression P exercee sera egale a:

P=F/S

Si l'amplitude de la deformation est faible on poura ecrire la relation suivante:

P=m.f

avec: m=E / (2.(1+s))

ou m : constante de proportionalite,s'appelle ;odule de rigidite.La valeur de m est

egale generalement a peu-pres a la moiyie de celle de E.

II-1-4 Le Module d'Incompressibilite K.

Si un corps est soumis a l'effet d'une pression P uniforme dans toutes les
directions,son volume diminuera de DV et on aura:

K=P / (DV/V)

La valeur K ainsi definie est appelee module d'incompressibilite,son inverse est le

module de compressibilite.Il existe encore les coefficients m et l appeles les
constantes de LAME.
Le coefficient l est exprime par :

17
 18

l= E / (1+s)(1-2s)

Dans le cas d'un corps homogene,isotrope et parfaitement elastique,les relations

suivantes relient les parametres precedement mentionnes:

m=1/2.E/(1+s)=3.K(1-2s)/2(1+s)

K=1/3.E/(1-2s)=2/3.(m(1+s)/(1-2s))

E=3K(1-2s)=2m(1+s)

s=1 / 2-E / 6K=(1-(2m / 3K)) / 2+(2m/3K)

Notons que,dans le cas d'un liquide parfait m et E = 0 et s=1/2

Le tableau 1 resume les relations entre les differentes constantes.

II-2 Ondes Sismiques:

Si un point quelconque d'un milieu isotrope est soumis a l'action d'une force F
variant dans le temps fig 3 celle-ci provoque des deformations qui varient suivant la
meme loi que la force appliquee.Les deformations se transmettent a travers le
corps dans toutes les directions avec une vitesse qui est fonction des proprietes
elastiques du milieu.Les particules du milieu se trouvant a certaine distance de ce
point se mettent a osciller successivement en transmettant les oscillations
(deformation) d'une particule a l'autre.Donc,on aura ce que l'on nomme ondes
elastiques ou sismique.Si l'on considere un milieu solide et homogene,plusieurs
types d'ondes elastiques peuvent s'y propager a s'avoir :des ondes
longitudinales,des ondes transversales et ondes dites superficielles.

II-2-1 Les Ondes Logitudinales:

C'est S.D. Poisson qui,en 1830,a montre qu'un corps homogene et isotrope de
dimensions illimitees peut transmettre des ondes de deux vitesses differentes et
que a une grande distance de la source de perturbation,le mouvement transmis par

18
 19

l'onde la plus rapide est longitudinalem c'est-a-dire que le deplacement des

particules est parallele a la direction de propagation de l'onde,le mouvement
transmis par l'onde la plus lente est transversale,c'est-a-dire que le deplacement
est a angle droit suivant la direction de propagation;en d'autre terme,pendant la
propagation d'une onde longitudinale les particules du milieu sont soumises a des
deformations de volume sans rotation coincidant avec le sens avec le sens de
propagation de l'onde fig .Donc la deformation de volume donne naissance aux
ondes longitudinales.L'onde longitudinale P s'accompagne d'un changement de
volume;elle represente une succession rapide des espaces de compression et de
dilatation a travers le milieu.
La vitesse de propagation de l'onde longitudinale est liee a la densite du milieu et a
des constantes elastiques par la relation suivante:

Vp = ÖE(1-s)/r(1+s)(1-2s)

Vp = Ö(K+4/3m)r = Ö3k(1-s)/r(1+s)

S'il s'agit d'un liquide parfait m=0 :

Vp=ÖK/r

II-2-2 Ondes Transversales ou de Cisaillement:

Lorsqu'une onde transversale se propage dans un milieu,le mouvement des

particules qui subissent des deformations sans changement de volume est
perpendiculaire a la direction de propagation de l'onde,donc tangent a l'onde Fig
La celerite est egalement liee a la densite de la roche et a pour expression

Vs= ÖE / r.2(1+s) = Öm/r

Si on compare les expressions (1,5) a (1,7) et (1,6) a (1,7) relatives aux vitesses
des ondes longitudinales et celles des ondes transversales on obtient:

Vp / Vs = Ö2(1-s) / (1-2s) = ÖK/m + 4/3

19
 20

Etant donne que le coefficient de Poisson est generalement voisin de 0.25 il en

resulte que:

Vp/Vs = Ö3 = 1.732

Donc,les ondes longitudinales sont plus rapides que les ondes transversales.Les
ondes transversales ne peuvent exister que dans les solides.

II-2-3 Ondes Superficielles: Ondes de Rayleigh et Ondes de Love:

En 1885 Rayleigh a montre que le long d'une surface libre limitant un demi-espace
homogene,isotrope et elastique,des ondes peuvent etre transmises.Les particules
du milieu subissent un mouvement elliptique et retrograde par rapport a la direction
de propagation,le grand axe de l'ellipse etant habituellement vertical.
Donc le passage de l'onde de Rayleigh provoque le deplacement a la fois
horizontal et vertical dans un plan vertical,parallele a la direction de propagation de
l'onde.La vitesse de ces ondes aux quelles on a donne le nom d'onde de Rayleigh
est faible,comparee a celle des ondes longitudinales et transversales.Pour un
milieu donne,la vitesse de l'onde de Rayleigh est a peu pres de 0.9 fois celle de
l'onde transversale:

Vr = 0.9.Vs

On a egalement montre que de tels types d,ondes peuvent prendre naissance

meme si la limite n'est pas libre ou si le milieu n'est pas homogene.
Des ondes de ce genre se propagent,par exemple,le long de l'interface separant
deux milieux epais.On a donne a ces ondes le nom d'onde de Stoneley (Howell
1969)
Il existe encore un autre type d'onde superficielle,a s'avoir les ondes de Love ou
les ondes superficielles transversales.Il s'agit d'onde de cisaillement,les
mouvements des particules etant paralleles a la surface et a l'angle droit avec la
direction de propagation.
Les ondes de Love:

20
 21

Elles prennent naissance dans le cas ou Vs1 < Vs2,le passage de l'onde de Love
provoque un deplacement le meme que celui des ondes S sans mouvement
vertical
Les ondes de surface constituent des parasites en prospection sismique,mais
utilisees en seismologie.Leur generation est attribuee a l'heterogeneite des
terrains.Leur passage limite au voisinage de la surface,les particules sont
soumises a une combinaison des mouvements longitudinaux et transversaux dont
l'amplitude diminue en fonction de la profondeur.Ces ondes constituent une
source de bruits indesirables dans l'enregistrement sismique et qu'on tente
generalement d'attenuer pour une meilleur exploitation des donnees sismiques.On
pense que les ondes de Rayleigh constituent le principal composant du ground-
roll*,parasite souvent tres genant qui fort heureusement n'existe pas en sismique
marine.Le ground-roll,se propage a vitesse relativement faible (150 a 1500 m/s)
avec un caractere basse frequence,dont on ne peut total ement se debarasser.
L'enregistrement des vibrations provenant d'un tremblement de terre eloigne
permet de mettre en evidence des ondes superficielles.Les mouvements
ondulatoires reveles par les sismogrammes montrent que leur periode est de
l'ordre de 10 s,d'ou le nom general d'ondes longues qu'on leur a donne.

II -3 Facteurs affectant l'energie des ondes acoustiques

La redistribution de l'energie transmise par une source sismique dans les milieux
favorables a la propagation des ondes s'effectue de la maniere suivante:
1- Soit par reflexion et transmission au niveau des interfaces
2- Soit par expansion geometrique du front d'onde,ce phenomene etant appele
"divergence spherique".

- La dispersion est due aux irregularites du milieu traverse

- L'absorption est une conversion de l'energie acoustique en chaleur.

II-3-1 La divergence spherique

Au cours de l'expansion du front d'onde incident (divergence spherique) l'energie emise

est repartie sur une surface croissante.Par consequent,la densite d'energie par
unite de surface decroit proportionnellement au rayon du front d'onde.
L'attenuation provoquee resulte de phenomenes purement geometriques.

21
 22

appelee " la divergence spherique".

Pour corriger cette decroissance d'energie,plusieurs methodes sont proposees :
Gutenberg (1963);O'brien et Lucas (1971) ;Newman (1973).
Elles tiennent compte de la stratification des terrains,de la structure des reflecteurs
et de la distance separant le recepteur de la source (offset).

II-3-2 Attenuation par reflexion et transmission

L'energie incidente est egalement redistribuee a chaque interface* sous forme

d'energie reflechie ou refractee.
* Interface: C'est la limite separant deux milieux de proprietes physiques differentes.En
sismique,le facteur influent connu est l'impedance acoustique,produit de la vitesse par la
densite;une variation de l'un de ces deux parametres suffit alors a creer une interface et a
engendrer cette repartition d'energie.

Dans le cas d'une onde plane a incidence normale et d'un modele bicouche le
rapport de l'energie reflechie (Er) sur l'energie incidente (Ei) est donnee par la
relation suivante:

| Er/Ei|i=0 = [( Z2 - Z1 )/(Z2+Z1)]2 = [(V2.r2 - V1.r1)/(V2.r2 + V1.r1)]2

ou : Z1 est l'impedance acoustique du premier milieu

V1 sa vitesse et r1 sa densite ; Z1= V1.r1

Z2 est l'impedance acoustique du dernier milieu

V2 sa vitesse et r2 sa densite ; Z2= V2.r2

En terme d'qmplitudes cette relation devient:

| Ar/Ai|i=0 = ( Z2 - Z1 )/(Z2+Z1) = (V2.r2 - V1.r1)/(V2.r2 + V1.r1)

CHAPITRE III - PARAMETRES CARACTERISANT

UNE OSCILLATION:

III - 1 - Introduction :

22
 23

La source sismique engendre une impulsion tres breve de l'ordre de la

milliseconde,a partir de laquelle on pourrait esperer enregistrer une trace
"impulsionnelle",ou les arrivees des ondes seraient detectees avec une grande
precision.En effet ,l'onde sismique subit des filtrages,tant par terrain que par
l'appareillage dont l'effet est de modifier l'impulsion initiale par un train d'ondes
plus ou moins long.L'onde sismiaue se propage dans un espace a trois
dimensions.Sa representation graphique tridimensionnelle s'avere difficile,c'est
pourquoi,on se borne a la representation de l'amplitude de vibration "u" d'un point
du milieu traverse par l'onde sismique en fonction du temps.
La courbe presentee sur la fig c'est-a-dire :

u=f(t)

Pour "x"=Cste est appelee "trace sismique ou sismogramme".La trace sismique a

l'aspect d'une impulsion amortie avec le temps;le temps "t1" du commencement de
la vibration est appele "impetus".Aux temps t2,t4,t6 etc..;on observe des bosses(les
"maximums "et les "minimums").
Le temps:

T=t6 - t2

Est appele "periode".Elle indique en unite de temps le nombre de fois ou le corps

oscillant passe par le meme etat.
La valeur inverse de la periode est la frequence lineaire:

f=1/T ; w=2pf

w:etant la frequence circulaire,l'unite de la frequence lineaire est le Hertz celle de

la frequence circulaire est rad/seconde.Il y a une autre representation des
caracteres des oscillations des differentes particules du milieu en temps donne
t=Const. c'est-a-dire on etudie la fonction u=f(x),t=const,ou u etant les valeurs des
deplacements des differents points donnes.Cette fonction est nommee "profil
d'onde".Les differents points du milieu sont deplaces de leur position stable aux
differentes valeurs au moment donne,comme il est presente sur la fig .La distance

23
 24

l entre les bosses de meme signe est appelee "la longueur d'onde".Elle est liee a
la periode par la relation suivante:

l = V / f = V.T [m]

Pour les ondes longitudinales,l'espace correspondant a deux mouvements de

dilatation ou a deux mouvements de compression caracterise la longueur d'onde.

K = 1 / l = f / V = 1 / V.T [m-1]

En prospection sismique on etudie encore la valeur inverse a la longueur d'onde

appelee "nombre d'onde".

III - 2 - Oscillations Harmoniques et Periodiques

Il est evident que les ondes elastiques d'une impulsion de courte duree ,sont
classees parmi les oscillations dites "non-periodiques",c'est-a-dire parmi celles
dont la forme ne se repete jamais.L'etude de ce type d'sentees par la relation
suivante:

F(t)=F(t+nT); (n=0,1,2,..)

Les formes d'impulsion se repetent periodiquement fig...Pour les oscillations

harmoniques les variations d'un parametre physique s'effectuent suivant les lois de
"sinus ou cosinus".Si l'on designe par u(t) la deviation d'une particule quelconque a
un instant "t" on a :

u(t)=a0sin(wt+j)=a0sin(2pf.t+j)

ou : a0 etant l'amplitude initiale

j phase initiale
w frequence circulaire.
lorsque j=p/2 on a une oscillation harmonique cosinusoidale.Les oscillations
harmoniques se caracterisent par les parametres suivants:Amplitude;Frequence (f)
ou Periode (T);Phase (j).

24
 25

III-3 Superposition des Oscillations

Le resultat de superposition (addition) des oscillations harmoniques engendrees

simultanement presente un grand interet en sismique.Soient deux oscillations
harmoniques u1(t) et u2(t) de meme freqoscillations est assez difficile;c'est
pourquoi on tache de ramener l'etude des oscillations nom-periodiques a l'etude
des oscillations periodiques et le plus souvent aux oscillations harmoniques.Ces
dernieres ne representent qu'un type,le plus simple parmi les oscillations
periodiques.Les oscillations periodiques peuvent etre preuence,mais distinguees
par leurs amplitude a1 et a2 et leur phases initiales j1 et j2:

u1=a1sin(2pf.t +j1)

u2=a2.sin(2pf.t+j2)

L'oscillation resultante sera la superposition des deux oscillations considerees:

U(t) = u1(t)+u2(t)
= a1sin (2pf.t+j1) + a2.sin(2pf.t+j2)
= a1sin 2pf.t.cosj1+ a1.sinj1.cos2pf.t+a2sin 2pf.t.cosj2+a2.sinj2.cos2pf.t
=sin 2pf.t(a1.cosj1+a2.cosj2) + cos2pf.t(a1.sinj1 + a2.sinj2)

En designant:

a.cosj=a1.cosj1+a2.cosj2 et a.sinj=a1.sinj1 + a2.sinj2

Þ U(t) = a.sin(2pf.t +j)

a et j peuvent etre determines en resolvant les equations precedentes

mentionnees ci-dessus:

a=Öa12+a22+2a1a2.cos(j1-j2)

j=arctg (a1.sinj1 + a2.sinj2) / (a1.cosj1+a2.cosj2)

25
 26

Donc,on peut constater,que le resultat de la superposition de deux oscillations

harmoniques differentes par les amplitudes et les phases initiales est une
oscillation harmonique de meme frequence,mais avec l'amplitude et la phase
initiale exprimees par les relations precedentes.Il est a noter que,le resultat sera
plus complique si on considere la superposition des oscillations de frequences
differentes.Dans ce cas le resultat somme ne sera pas harmonique.

III - 4 - Decomposition des Oscillations

Theoriquement,un sismogramme a un contenu certain:l'ensemble des reflexions

consecutives a l'explosion,provenant de couches geologiques d'emplacement
determine.Mais si l'on effectue deux tirs an meme endroit,les deux sismogrammes
different notablement,du fait de nombreux facteurs aleatoires.Les impulsions S(t)
engendrees par les explosions ne sont pas identiques,elles subissent une
deformation et un allongement avec la profondeur et l'absorption.Enfin,plusieurs
types de bruits s'ajoutent aux signaux sismiques.C'est pourqoi il est difficile de
comparer les oscillations engendrees par differentes sources.Pour eliminer les
ondes parasites,les sismiciens procedent a des analyses frequentielles.
On utilise le plus souvent pour la decomposition des oscillations:la Transformee de
Fourier et l'integral de Duhamel.Afin de mieu comprendre les principes generaux
de decomposition des oscillations un rappel mathematique serait
utile.Notemment,on se contentera de mentionner les formules les plus courantes et
principes generaux utilises en sismique.

III -5- Decomposition Spectrale des Oscillations:

Si "t" est une variable reelle,f(t) une fonction aux valeurs relles ou complexes
definie en tous les points et indefinement derivable sur la droite relle,et de plus a
decroissance rapide a l'infini ( conditions de Dyrechlet) alors on peut lui faire
correspondre une fonction que l'on nomme " Transformee de Fourier ".

F(t)=1/2pòs(w) ejwt dw

ou bien

26
 27

S(w)=1/2pòF(t)e-jwt dt

La fonction S(w) est une fonction complexe de variable w (frequence circulaire ou

pulsation).Une fonction f(t) peut etre representee comme la somme de la fonction
paire F1(t) et impaire F2(t).En appliquant la formule d'Euler on peut obtenir la
relation suivante:

S(w)=ò[F1(t)+F2(t)]e-jwt dt = ò[F1(t)coswt-jF(t)sinwt+F2(t)coswt-jF1(t)sinwt]

En tenant compte que la fonction "coswt" est paire

et la fonction "sinwt" est impaire on peut ecrire:

S(w) = 2òF1(t)coswtdt -2jòF2(t)sinwtdt = A(w) + Jb(w)

F(t) n'existe que dans le domaine t³0

ou :

Re{S(w)}=A(w)=2òF1(t)coswtdt

Im{S(w)}=B(w)= -2òF2(t)sinwtdt

ou : Re{S(w)} = A , Partie reelle

Im{S(w)} = B , Partie imaginaire de la fonction S(w)
S(w) peut etre represente sous la forme suivante:

S(w) = S(w).ejj(w)

S(w): etant le module

j(w):variable de la fonction complexe S(w).
Ils sont liee aux fonctions A(w) et B(w) par les expressions suivantes:

S(w) = ÖA2(w) + B2(w)

j(w) = Arctg B(w) / A(w)

27
 28

La fonction S(w) est nommee le spectre d'amplitude,j(w) etant le spectre de

phase.La partie relle A(w) et le module S(w) du spectre complexe S(w) sont les
fonctions paires de la frequence w,par contre la partie imaginaire B(w) et la variable
j(w) sont les fonctions impaires:

A(-w) = A( w) ;B(-w)=-B(w) ;S(-w) = S(w) ;j(-w) = -j(w).

En effet la determination du spectre d'une oscillation consiste en la determination

des deux fonctions relles S(w) et j(w) .Enfin,on nomme "Transformee en cosinus "
la fonction :

C(w) = òF(t)coswt dt = 2òF(t)coswt dt

En vertu de cette expression toute fonction (ou oscillation) f(t) peut etre consideree
comme la somme d'un certain nombre de cosinusoides de freauences
differentes,chaque cosinusoide etant caracterisee par sa propre amplitude et sa
phase.Autrement dit, chaque oscillation complexe (nom-periodique) peut etre
decomposee en oscillations harmoniques simples. fig

III - 6- Representation des Spectres:

La representation des oscillations periodiques ou non-periodiques par la somme

des oscillations harmoniques est appelee decomposition spectrale.Cette
decomposition est unique,c'est-a-dire,toute oscillation ne peut etre presentee qu'a
l'aide d'une seule combinaison des oscillations harmoniques.Le resultat de la
decomposition spectrale d'une oscillation periodique se presente un spectre de
raies fig ,ou sur l'axe d'abscisse on porte la frequence lineaire "f" et sur l'axe des
ordonnees les amplitudes S(w).Le spectre des oscillations periodiques et
harmoniques est compose de lignes equidistantes Le spectre d'amplitude des
oscillations non periodiques ou complexes est une courbe continue,dans le
systeme des coordonnees "frequence-amplitude" fig ,puisque les intervalles entre
des lignes spectrales diminuent infinement,comme l'oscillation decomposee etant
la fonction non-periodique.Il comprend un grand nombre d'oscillations harmoniques
ayant des amplitudes et des frequences differentes.Enfin pour construire les
spectres des phases,on porte sur l'axe des ordonnees les phases initiales et sur
l'axe des abscisses les frequences fig .Maintenant,considerons les proprietes

28
 29

principales des Transformees de Fourier qui sont tres importantes en "pratique"

dans la decomposition spectrale:

III-6-1- Propriete de Conformite Reciproque et Synonyme:

A toute impulsion F(t) correspond un spectre S(w) et inversement a tout spectre

S(w) correspond une fonction F(t) c'est-a-dire la decomposition spectrale de toute
oscillation est unique.Si les spectres complexes de deux oscillations coincident,on
peut dire qu'elles sont identiques.

III-6-2- Addition des Spectres:

On sait que les Transformees de Fourier sont lineaires,c'est pourquoi on peut dire:
si S1(w) et S2(w) sont respectivement les spectres complexes des fonctions F1(t)
et F2(t) la somme sera:

S(w) = S1(w) + S2(w)

ou : S(w) est le spectre complexe de la fonction sommaire:

F(t) = F1(t)+F2(t)

III-6-3- Similitude des oscillations et des spectres:

Soit deux fonctions semblables:

F1(t) et F2(t) = kF1(t)

ou k etant une constante arbitraire ,qui ne depend pas de "t" en sachant que les
deux Transformees de Fourier sont lineaires,les spectres complexes S1(w) et S2(w)
de ces fonctions,sont generalement semblables,autrement dit:

S2(w) = kS1(w).

III - 6 - 4 - Spectre d'une Oscillation Decalee en Temps

29
 30

par rapport a une autre Oscillation

Soit deux oscillations :

F1(t) et F2(t)=F(t+t)

ou: t est une constante,dans ce cas:

S2(w) =S1(w) ejwt

Les spectres d'amplitude de ces oscillations sont identiques

S1(w) =S2(w) ; mais le spectre de phase de l'oscillation decalee sera:

j2(w)=j1(w)+wt

III -6 - 5 - Spectre d'une Oscillation Serree

ou Comprimee sur l'Axe des Temps

Soit deux oscillations :

F1(t) et F2(t) =F1(pt) ou: "p" = const.

On distingue trois cas:

1- p=1 le signal est constant

1- p>1 le signal se comprime
2- p<1 le signal se dilate

S2(w)=ò[F(pt)e-jwt dt = 1/pò[F1(t1)e-j(w/p)t1 dt1 = 1/p(w/p)

ou: t1=pt ; t=t1/p ; dt=d(t1/p) = (1/p)dt1.

Par consequent,lorsque la duree d'une oscillation augmente de "p" fois sur l'axe
des temps,le spectre d'amplitude se serrera de la meme proportion sur l'axe des
frequences et inversement.

30
 31

III -6 - 6 Spectre de Certaines Impulsions.

Examinons quelques impulsions theoriques et leurs spectres d'amplitude

III-6-6-1- Cas d'une Oscillation Sinusoidale de la Forme F(t)=asinw0t

Le spectre d'amplitude de cette oscillation de duree infinie fig

III-6-6-2- Cas d'un Segment de la Sinusoide ayant l'expression suivante:

0 t£0

{
F(t)= sin w0t 0 < t < 2pn/w0
0 t³ 2pn/w0

ou "n" nombre de periodes;le spectre complexe de l'oscillation est:

S(w)=òsinw0te-jwt dt = [2w0sinnp(w/w0)] /(w20 - w2) ej(npw/w0 - p/2

Le spectre d'amplitude sera:

S(w)=2w0sinnp(w/w0)] /(w20 - w2)

Considerons les spectres d'amplitude pour n=1,2,3

III-6-6-3- Une Oscillation Sinusoidale Amortie

Elle est definie par l'equation suivante:

F(t) = {0 quand t£0 et e-atsinw0t quand t>0 ; son spectre complexe est:
S(w)=òe-(a+jw)t sinw0tdt = w0/(a+jw)2w20.

Le spectre d'amplitude s'exprime par :

S(w ) = w0/Ö(a2+w20 - w2)2 + 4a2w2

III-6-6-4- Une fonction Rectangulaire Unique F(t)

31
 32

Est Donnee par la Formule Suivante:

1(t) = {0 quand t<0 et 1 quand t ³0 ; son spectre complexe est:

S(w)=òe-jwt dt=1/jw

le spectre d'amplitude est egal a:

S(w ) = 1/w

III-6-6-5- Impulsion Rectangulaire de Duree t:

Qui s'exprime par l'equation suivante: F(t)={0 quand t£0 et A quand 0<t<t et 0
t³t}On peut considerer cette fonction comme la difference de deux fonctions
rectangulaires uniques:

F(t)=A[1((t) - 1(t-t)].

Le spectre conplexe S(w ) de l'impulsion rectangulaire sera egal a:

S(w ) = [2Atsin(wt/2)/wt]ejwt

Le spectre d'amplitude de cette impulsion peut etre calcume d'apres l'expression

suivante:

S(w ) =2Csin(wt/2)/wt

ou C=At etant l'aire de l'impulsion consideree.

III-6-6-6- Impulsion de Dirac (unitaire) F(t)=d(t).

L'impulsion dite de Dirac d(t), joue le role de l'impulsion unite,dans les problemes
de physique impulsionnelle.La forme de cette impulsion est indifferente a condition

32
 33

que la base soit infinement petite,elle peut prendre 2 formes classiques:1- la forme
d'un rectangle de largeur zet de hauteur 1/z dans les problemes ou l'on doit
echantillonner des courbes. 2- La forme plus "physique"d'une cloche tres etroite et
d'une hauteur assez grande dans les autres problemes.Impulsion de Dirac est
exprimee de la maniere suivante:
ì¥ t<0
F(t)=íòF(t)dt = {1 si t1<0<t2 ou bien 0 sit1t2>0
î0 t¹0

Le spectre complexe est:

S(w)=1=S(w)

et coincide avec le spectre d'amplitude:S(w)=1 C'est-a-dire que le spectre

comprend toutes les frequences d'egale mesure.
Nous avons deja note que des impulsions sismiques ne durent que pendant un
temps tres bref et elles ont la forme d'une impulsion croissante-amortie.C'est
pourqoi il est important de considerer les spectres de quelques impulsions qui
rssemblent aux impulsions sismiques.L'impulsion de Berlague et l'impulsion de
cloche sont de cette serie d'impulsions.

III-6-6-7- Impulsion de Berlague:

Elle s'exprime par la relation suivante:

F(t)=ít.e-2pat.sinf0t t>0 et 0 pout t£0

ou t.e-2pat est l'amplitude d'impulsion.son spectre complexe peut etre trouve a

l'aide de l'expressio suivante:

S(f)=2f0/(2p)2{a[(a2+f 2)2 - f 4+4a2f 2]+jf[(f02-f 2)2+a4+4a2f 2]/[(a2+ f 2 -f 2 )2 4a2 f 2]2

III-6-6-8- Impulsion en Cloche:

Nous avons examine les spectres de plusieurs impulsions et on peut constater que
les formes des spectres obtenus ne coincident pas avec celles des

33
 34

impulsions.C'est bien le cas de dire que les formes des spectres et des impulsions
different toujours.Mais il y a une exception de cette regle generale et a titre
d'exemple on peut citer l'impulsion ayant la forme de cloche et qui est nommee
"impulsion de cloche": a.e-(2pat)2 Cette impulsion est symetrique par rapport a
l'origine de systeme de coordonees et son spectre complexe ne contient que la
partie reelle:

S(f) = a.e-(f/2a)2/2aÖp

Il est commode d'avoir un multiplicateur egal a une et c'est pourquoi nous

remplacons " a" par: 2aÖp :

a = 2aÖp

on obtient alors:

F(t) =2aÖpe-(2pat)2

Þ S(f) = e-(bf)2

ou b=1/2a. Nous avons examine certaines impulsions theoriques et leurs

spectres.Maintenant notons quelques proprietes importantes de la decomposition
spectrale des fonctions:
- Le maximum de spectre d'amplitude coincide avec la frequence apparente
(visible) de l'oscillation.
- La forme du spectre d'impulsion depend de sa duree c'est-a-dire il y a une liaison
inversement proportionnelle entre elles.a)- Lorsque la duree de l'impulsion
croit,son spectre d'amplitude devient plus etroit; b)- Lorsque la duree de l'impulsion
diminue,son spectre s'elargit et la frequence maximum fmax ne coincide plus avec
la frequence apparente (visible).

III -6 -7 - Decomposition par des Fonctions Rectangulaires

- Integrale de Duhamel - Convolution

34
 35

On peut presenter une fonction quelconque comme la somme des impulsions

rectangulaires (notons que dans le cas des transformees de Fourier une fonction
peut etre presentee par la somme des oscillations harmoniques).Prenons une
impulsion rectangulaire determinee par la relation suivante:f(t-t)= 0 si t<t et 1 si
t³tConsiderons,a present une fonction quelconque F(t) presentee sur la fig et
representons cette fonction comme la somme des fonctions rectangulaires.La ligne
brisee OABCDE etc.. designons la par le terme F1,on peut ecrire alors:

F1(t)=OA.1(t)+BC.1(t-dt)+DE.1(t-2dt)+...;

ou : OA=F(0) , BC=dF(dt) , DE=dF(2dt)

dF(tk) etant l'accroissement de la fonction dans le segment tk et tk+dt;on peut

reecrire l'equation precedente comme suit:

F1(t)=F(0).1(t)+ådF(kdt)1(t-kdt) ; t=(n-1)dt

La fonction F1(t) coincide avec la fonction initiale F(t) aux points A,C... En prenant
dt infinement petit on peut avoir la coincidence satsfaisante des fonctions F1(t) et
F(t)

F(t)=F(0)1(t)+ò1(t-t)dF(t)

si F(0)=0 :

F(t)=ò1(t-t)dF(t) ; ou t = kdt

Cette equation est nommee Integrale de Duhamel ou de Convolution.

III -6 - 8 - Convolution:

Soient deux fonctions S(t) et G(t),dont les Transformees de Fourier sont

respectivement S(w) et G(w);soit :

R(w)=G(w).S(w)

35
 36

le produit de ces deux transformees,lui meme transformee de Fourier de R(t);on

demontre que:

R(t)=òS(t)G(t-t).dt=òS(t-t)G(t).dt

que l'on nomme produit de convolution,et que l'on ecrit:

R(t)=S(t)*G(t).

Dans le plan reel,on nomme fonction d'auto-correlation (en abrege FAC) le produit
de convolution:

òG(t).G(t-x)dt = G(x)*G(-x)

considerons un exemple:L'impulsion de Dirac d(t) a pour transformee de Fourier la

fonction D(w)=1,et inversement on peut demontrer que:

S(t)*d(t) = S(t)

d(t) est donc l'unite de convolution,qui constitue un element neutre.Si l'impulsion se

produit a l'instant t elle a pour expression:

d(t-t)

et òd(t-t)dt =1

d(t-t) est la transformee de Fourier e-j2pft et S(t)*d(t-t)=S(t-t)

La convolution est equivaut a un changement d'origine ou translation.

III-6-8-1 - Pratique de la Convolution

Correlations:Une correlation entre deux traces consiste a rechercher dans quelle

mesure elles se ressemblent.En un temps t donne le produit de leurs amplitudes

36
 37

G1(t) .G2(t) fournit par sa valeur et son signe ,une mesure de cette
ressemblance.Dans l'intervalle de t1 t2,la ressemblance des deux traces est
exprimee par l'integrale:

òG1(t).G2(t)dt

Une tres forte valeur positive indique que les traces se ressemblent beaucoup;une
valeur nulle signifie que la correlation n'a decelee aucune ressemblance entre les
traces;on obtient une forte valeur negative pour des traces qui se ressemblent,mais
qui sont en opposition de phase.Si l'on veut etudier les variations de cette
ressemblance entre les traces lorsaue celles-ci se deplacent l'une par rapport a
l'autre,on doit considerer,le produit de convolution des deux traces:

òG1(t).G2(t-t)dt = G1(t)*G2(-t)

que l'on nomme fonction de cross-correlation.L'autocorrelation est la comparaison

de la trace a elle meme:

òG(t).G(t-t)dt

c'est une courbe symetrique par rapport a une valeur maximum correspondant a la
superposition de la trace a elle meme.Considerons un exemple:supposons au'il
faut trouver la convolution de deux series X et K representees par des successions
suivantes:

X Þ x0 x1 x2
K Þ k0 k1 k2

On trouve la convolution de la maniere suivante:

1ere Etape:

x0 x1 x2
k2 k1 k0

y0 = x0.k0

37
 38

2eme etape:
x0 x1 x2
k2 k1 k0

y1 = k1x0 + k0x1

3eme etape:
x0 x1 x2
k2 k1 k0

y2 = k2x0 + k1x1 + k0x2

4eme etape:
x0 x1 x2
k2 k1 k0
y3 = k2x1 + k1x2

5eme etape:
x0 x1 x2
k2 k1 k0

y4 = k2x2

Le resultat de la convolution nous donne la succession yt:

yt=k0x0;(k1x0+k0x1);(k2x0+k1x1+k0x2);(k2x1+k1x2);k2x2

CHAPITRE IV - PROPAGATION DES ONDES SISMIQUES

(Front,Rayon,Isochrone)

IV-1 Introduction:

Considerons un milieu homogene,isotrope et une source ponctuelle d'energie S

engendrant une onde elastique de duree tres breve dt (telle l'explosion d'une

38
 39

petite charge de dynamite) fig ;La force engendrant des ondes elastiques au
point S n'agit que pendant un temps bref.Par consequent,les oscillations ne
peuvent subsister qu'aux points de coordonnees suivantes :

Vp.(t1-dt) £r1£Vp.t1

ou Vp:vitesse de propagation des ondes longitudinales;t1:temps considere;r1:est

la largeur du domaine d'existance des oscillations.Si l'on considere la position des
oscillations en un temps t2 superieur a t1 (t2>t1),alors la condition suivante doit
etre satisfaite:

Vp(t2-dt) £r2£Vpt2.

En analysant ces deux inegalites on peut conclure que les oscillations se

propagent en s'eloignant de la source S.Au moment t1 les oscillations n'existent
qu'entre deux spheres de rayon:

r11=Vp(t1-dt) et r22=Vpt1

Mais pour le moment t2,le rayon de ces spheres sont:

r12=Vp(t2-dt) ; r22=Vpt2.

On peut conclure,que si dans un milieu une impulsion de duree dt se propage a

tous les temps t>dt,on peut distinguer deux domaines:

1 - Le domaine (sphere) interne de rayon Rin:

Rin=Vp.(t-dt)

dans lequel les oscillations engendrees par la sources se sont deja amorties.

2 - La couche spherique de rayon interne (Rin) et d'epaisseur "r":

39
 40

r=Vp.dt

dans laquelle les oscillations subsistent.

3-Le domaine exterieur enveloppe par la sphere de rayon Rex:

Rex = Vp.t

Les oscillations engendrees par la source S n'ont pas atteint ce domaine.La

surface limitant a l'exterieur de couche dans laquelle les oscillations existent et
appelee "front primaire" de l'onde (front avant).La surface limitant de l'interieur
cette couche est appelee (front arriere).En un mot on appelle Front d'Onde la limite
entre le domaine ou les particules du milieu ont deja mises en mouvement et celui
ou elles sont encore au repos.Le front d'onde est mobile et a chaque instant du
temps :t0,t1,t2,t3,... tm;il occupe une nouvelle position.Pour determiner la position
d'un front a un temps "t" on a introduit la notion d'isochrone.On appelle
"isochrone"une surface dont tous les points ont le meme temps t=const.Les
isochrones sont des surfaces qui coincident avec le front a des momemts
differents.Chaque isochrone est caracterisee par son propre temps.fig .On
distingue quatre types d'ondes suivant la forme du front d'onde a s'avoir:
1-onde spherique
2-onde plane
3-onde cylindrique
4-onde conique.
Sur les fronts de ces ondes les deviations des particules sont pareilles et en
phases.Dans le milieu homogene l'onde issue d'une source ponctuelle est
spherique.Si une onde spherique a un rayon important il est possible d'assimiler
cette onde a une onde plane. On peut alors plus simplement considerer les rais
perpendiculaires au front d'onde,rais a qui l'on donne le nom de rayon
sismique.Les rayons des ondes elastiques dans le milieu homogene et isotrope
sont les lignes droites ;dans le milieu heterogene et anistrope les rayons sont
curvilignes.Mais il est a noter que la conception de propagation par rayon n'a pas
ici un caractere physique,la seule realite physique etant le front d'onde.L'ensemble
d'isochrones est appele "champ de temps" si Dt est l'accroissement du temps
t=t(x,y,z) le long du rayon sismique on peut ecrire:

40
 41

Dn=V(x,y,z)Dt

ou Dn est la distance entre deux isochrones;v(x,y,z) etant la vitesse.En passant a

la limite Dn ® 0 on peut determiner le gradient de temps :

t = [ grad t ] = dt/dn=1/v(x,y,z)

et en remplacant "n" par ses composants (x,y,z) on peut ecrire:

(¶t/¶x)2(¶t/¶y)2(¶t/¶z)2 = 1/V2(x,y,z)

Cette equation est nommee Equation de Champ de Temps ou bien l'equation

d'Hamilton.La solution est une fonction t=t(x,y,z),qu'on appelle fonction de
temps.Pour l'ensemble des isochrones on a:

t(xi,yi,zi) = ti

ou ti:etant le temps correspondant a une isochrone fixe choisie parmi un ensemble

d'isochrones.

IV - 2 Principes de la Sismique Geometrique:

La sismique geometrique se conforme aux meme lois que celles de l'optique

geometrique,on examine les indices geometriques caracterisant la propagation des
ondes.
telle que:temps de parcours ,vitesse,trajectoires,vitesse apparentes,front d'onde,
hodochrone.
Ce sont les caracteristiques dites cinematiques.La sismique geometrique se fonde
sur quatre principes:

IV 2-1- Principe de Huyguens :

41
 42

En vertu de ce principe,chaque point du front d'onde "f1"peut etre considere

comme la source d'une nouvelle onde spherique fig .Le nouveau front d'onde "f2"
apres un certain laps de temps "t"est alors l'enveloppe mathematique de toutes les
nouvelles petites ondes elementaires appelees ondelettes.Les vibrations qu'elles
engendrent se detruisent partout mutuellement par interference destructive sauf
sur une surface qui est l'enveloppe mathematique des ondelettes et qui constitue
elle meme un front d'onde.Dans un milieu homogene et isotrope (vi=const)les
fronts elementaires sont spheriques et de rayon :

r = vt

Si le milieu est heterogene (ou la vitesse est variable (v¹cste les rayons des
ondelettes seront differents pour chaque point du front F1

rk = vk.t

Dans ce cas l'enveloppe des ondelettes constitue egalement un front initial F1 fig
Donc,suivant le principe de Huyguens le front d'onde est la surface que l'onde
sismique a atteint a un instant donne.

IV-2-2- Principe de Fermat:

Tout d'abord rappelons que les trajectoires orthogonales aux fronts d'ondes sont
les rayons ou rais sismiques.Le temps de parcours de l'onde le long du rayon
sismique est inferieur au temps de parcours par tout autre chemin.Autrement
dit,l'onde s'efforce de parcourir d'un point a l'autre pendant un temps le plus court
fig C'est justement le principe de Fermat.Il est evident que,dans un milieu
homogene les rayons sismiques sont les lignes droites,mais dans des milieux
heterogenes ce sont les lignes brisees ou curvilignes.

IV-2-3- Principe de reciprocite:

Le temps de parcours,le trajet et le caractere d'ebranlement ne subissent pas de

changement si l'on permute la source d'onde (le point de tir) et le sismographe fig

42
 43

IV-2-4- Principe de superposition:

Si dans un milieu deux ondes se propagent simultanement,chacune de ces ondes

se propage comme si l'autre n'existe pas (elles se propagent
independement).Autrement dit ,on peut construire les isochrones de chaque onde
sans tenir compte de la propagation de l'autre.

IV-3- Hodochrone . Vitesse Apparente:

En Sismique les hodochrones jouent un grand role dans le traitement et

l'interpretation On appelle hodochronne une courbe "temps-distance"c'est-a-dire
une courbe qui est construite d'apres le temps d'arrivee des ondes sismiques aux
differents points de reception.(aux sismographes situes sur la surface de la
terre).Supposons que l'on observe les temps d'arrivee d'une onde sismique a partir
de deux sismographes S1 et S2 situes sur la surface de la terre fig L'onde arrive
au point S1 au temps t1,elle atteint le point S2 au moment:

t2=t1+Dt.

La vitesse vraie de propagation de l'onde sismique est:

V=DS/Dt.

Mais a la surface de la terre l'onde atteint le sismographe S2 en parcourant la

distance Dx=S1®S2 a la vitesse Va qui est nommee vitesse apparente:

Va=Dx/Dt.

La vitesse apparente est une vitesse de propagation du front d'onde le long de la

ligne d'observation (ou bien du profil).Il n'est pas tres difficile de remarquer que
dans le triangle "QS1S2" on a :

DS=Dx.sina=Dx.cose

ou a:est l'angle entre le raie sismique et la normale qui s'appelle "angle

d'incidence"de l'onde (e:angle d'emergence).On peut ecrire que:

43
 44

V=Dx.sina/Dt.=Dx.cose/Dt=Va.sina=V.cose

d'ou:

Va=V/sina=V/cose.

En s'achant que sina<1 et cose<1,on peut constater que la vitesse apparente est
en general superieure a la vitesse vraie:

Va>V pour a<p/2 ; si a=p/2 on a: Va=V

C'est-a-dire que le front d'onde est orthogonal(le raie sismique est parallele a la
ligne d'observation fig )

si a=0,on a : Va®¥

Autrement dit,le front de l'onde est parallele a la ligne d'observation(le raie

sismique est perpendiculaire),la vitesse apparente a une valeur infinement grande
fig .La vitesse apparente peut etre negative ou positive en dependance de la
direction d'emergence de l'onde par rapport au point de tir fig . On peut calculer
la vitesse apparente a partir de l'hodochrone en effet :

Va=Dx/Dt [m/sec].

En sismique,on utilise encore la valeur inverse a la vitesse apparente que l'on

appelle "gradient de Temps":

t=[grad t]=[1/Va];[sec/m].

IV-4 L'Equation de l'Onde et les Solutions Particulieres

Pour engendrer les ondes sismiques il faut appliquer une force;lorsqu'un milieu
elastique initialement au repos est soumis a des forces de volume,l'equilibre initiale
est rompu et un deplacement se produit de facon a retablir un nouvel equilibre.En

44
 45

sismique on etudie generalement,des actions concentrees et les ondes sismiques

se propagent dans le milieu a partir d'une source que l'on peut souvent considerer
comme ponctuelle.Le mouvement des particules du milieu homogene et elastique
pendant la propagation de l'onde sismique peut etre exprime par l'equation
differentielle d'ordre deux avec les derivees partielles:

¶2U/¶t2=V2(¶2U/¶x2+¶2U/¶y2+¶2U/¶z2)=V2 DU

ou U=U(x,y,z,t) est la fonction d'onde qui determine la loi de deviation des

particules dans tous les points et a tout moment "t";V etant la vitesse de
l'onde.Considerons quelques types de solutions de l'equation de l'onde:a) Tout
d'abord etudions la solution de l'equation pour une onde plane:La solution partielle
est:

U=F1(t-X/V)+F2(t+X/V).

Cette solution verifie,donc l'equation ci-dessus;en effet la derivee partielle de U par

rapport a t s'ecrit:

¶U/¶t=F'1(t-x/V)+F'2(t+x/V)

La derivee seconde de U est:

¶2U/¶t2=F''1(t-x/V)

La derivee partielle de U par rapport a "x" s'ecrit:

¶U/¶x=-1/VF'1(t-x/V) +1/VF'2(t+x/V).

La derivee seconde de U est:

¶2U/¶x2=1/V2F''1(t-x/V)+1/V2F''2(t+x/V);

En remplacant ¶2U/¶x2 et ¶2U/¶t2 par ses valeurs trouvees on peut demontrer que
:

45
 46

¶2U/¶t2+V2¶2U/¶t2

Exemple:Considerons une onde plane sinusoidale:

U=sinw(t-x/V)

¶U/¶t = wcosw(t-x/V) ;
¶U/¶x = - w/Vcosw(t-x/V) ;
¶2U/¶t2 = - w2sinw(t-x/V) ;
¶2U/¶x2 = - w2/V2sinw(t-x/V)

d'ou on peut ecrire:

¶2U/¶t2=V2¶2U/¶t2

A present examinons le cas general:La propagation d'une onde dans l'espace.La

solution de l'equation de l'onde est de la forme suivante:

U=F(mx+ny+kz+pt)

ou m,n,k,p etant les constantes qui determinent la position du front plan dans
l'espace.
Elles sont liees a la vitesse V et aux cosinus directeurs.
Exprimons V par les coefficients :m,n,k,p.

¶U/¶t = pF'(mx + ny + kz + pt)

¶U/¶x = mF'(mx + ny + kz + pt)
¶U/¶y = nF'(mx+ny+kz+pt)
¶U/¶z = kF'(mx+ny+kz+pt)

Les derivees secondes de U sont :

¶2U/¶t2 = p2F''(mx+ny+kz+pt)

46
 47

¶2U/¶x2 = m2 F''(mx+ny+kz+pt)
¶2U/¶y2 = n2F''(mx+ny+kz+pt)
¶2U/¶z2 = k2F''(mx+ny+kz+pt)

En injectant ces valeurs dans l'equation d'onde,on peut trouver :

p2F''(mx+ny+kz+pt) = V2(m2+n2+k2)F''(mx+ny+kz+pt)

V=p/Ö(m2+n2+k2)

La solution elementaire de l'equation d'onde s'ecrit egalement sous la forme

suivante :
U=F(xcosNx+ycosny+zcosnz+Vt)
ou :
cosNx=m/Ö(m2+n2+k2)
cosNy=n/Ö(m2+n2+k2);
cosNz=k/Ö(m2+n2+k2).

Sont les cosinus directeurs du raie sismique qui est perpendiculaire au front plan.
On peut transformer la formule en la representant de la maniere suivante:

U=F[x/VcosNx+y/VcosNy+z/VcosNz+t]

si l'onde sismique a la meme direction que celle de l'axe des abscisses,alors;

y=0;z=0;cosNx=1

il en resulte que:

U=F(t-x/V)

il est evident que pour une onde plane se propageant dans la direction de "-x" on
a:U=F(t-x/V).

47
 48

Considerons la solutionn de l'equation d'onde pour une onde spherique:elle

s'exprime comme suit:

U=1/r{F1(t-x/V)+F2(t+x/V)}

ou "r" est le rayon du front spherique.Cette solution de l'equation d'onde comprend

deux parties:-cas ou l'onde se propage suivant la direction "+x" et egalement por
l'onde se propageant dans le sens inverse "-x".L'energie de l'onde spherique
diminue avec la distance,parce que son energie est donc,repartie sur une surface
spherique,qui croit comme le carre du rayon de la sphere. Il en resulte que
l'energie par unite de surface est inversement proportionnelle au carre de la
distance de la source.L'elargissement du front spherique s'appelle "dilatation".La
perte de l'energie de l'onde spherique a cause de dilatation est importante au
voisinage de la source.Aux points plus eloignes de la source le front se redresse et
devient presque plan.De cette dispersion d'energie avec la distance ne sera pas la
seule cause d'attenuation de l'amplitude de l'onde,cette amplitude etant
proportionnelle a la racine carre de l'energie:

A=ÖE

Il faut encore tenir compte des pertes par absorption notamment de la dissipation
de l'energie elastique(frottement des particules du milieu). L'amplitude de l'onde
sismique peut etre calculee a l'aide de la formule suivante:

A=A0.e-ar/rn

ou "a" est le coefficient d'absorption;"n" etant l'exposant de dilatation Lorsque

a=0,l'absorption n'existe pas;n=1 pour l'onde spherique ; n=2 pour l'onde
cylindrique et n=0 pour l'onde plane.On voit donc,que le poids des charges
d'explosif deviendrea rapidement important principalement en sismique
refraction,ou la distance entre le point d'excitation et le point d'enregistrement croit
frequement de 1 a 100 km

IV - 5 - Reflexion et Refraction des Ondes Sismiques:

48
 49

Les ondes sismiques provoquees par une explosion ou par un ebranlement

artificiel se propageent dans tous les sens du point d'excitation;si le milieu est
homogene et isotrope les ondes sismiques suivent leurs parcours en profondeur
sans aucune perturbation jusaqu'a l'amortissement complet.Mais pour prevoir la
propagation des ondes dans un milieu heterogene et stratifie il faut rechercher les
phenomenes produits par la presence de frontieres limitant ce milieu.Tout d'abord
considerons le passage d'une onde sismique a travers la surface de separation de
deux ;ilieux de proprietes physiques differentes.Sur cette surface on a des
phenomenes de reflexion et de refraction fig .Les phenomenes de reflexion et de
refraction des ondes sismiques sont plus complexes qu'en optique,du fait qu'une
onde incidente longitudinale ou transversale donne naissance a deux ondes
reflechies et deux ondes refractees.L'une des deux ondes reflechie sera
longitudinale,tandis que l'autre sera transversale.Considerons deux milieux:fig
;"il" et "it" sont les angles de reflexion des ondes L et T.Les ondes sismiques se
reflechissent et se refractent suivant la loi de Snell-Descartes,qui s'exprime comme
suit:- Tous les rayons:incidents,reflechis et refractes sont situes dans le meme plan
normal a la surface de separation des milieux,qui s'appelle "plan des rayons"; -
L'angle d'incidence,les angles de reflexion et de refraction sont lies aux vitesses
par la relation suivante :

sini(l) / V1 = sini(i) / V1
= sini(r) / V1
= sina(l) / V2
= sina(r) / W2
= sini(r) / W1
= 1 / Va fig

Les valeurs inverses aux vitesses apparentes de toutes les ondes au point "A" sont
egales.Cela veut dire,que sur l'interface le front d'onde d'incidence ne se decroche
pas,seul la direction du front etant changee.On peut exprimer la loi de reflexion
dans le cas d'incidence d'une onde longitudinale Considerons tout d'abord l'onde
reflechie -longitudinale,qui se propage dans le milieu I Nous voyons que:

sini(l)/V1=sini'(l)/V1 ;i=i'.

Donc:l'angle d'incidence est egal a l'angle de reflexion.

49
 50

Maintenant,considerons l'onde reflechie-transversale.Nous voyons que:

sini(l)/V1=sini'(r)/W1 , V1¹W1 Þ i(l)¹i'(r)

Dans le cas d'incidence d'une onde transversale sur l'interface,pour l'onde

reflechie-transversale on a:

sini(r)/W1=sini'r/ W1 Þi(t)=i'(t).

En meme temps nous pouvons demontrer que pour l'onde reflechie longitudinale:

sini(t) / W1 =sini'(l)/V1 ; W1¹V1;i(t)¹i'(l).

La loi de refraction pour les ondes incidentes et refractees de la meme nature (soit
longitudinale soit transversale) peut etre exprimee de la maniere suivante:

sini/sina=V1/V2

Lors de la refraction des ondes sismiques,quatre cas peuvent avoir lieu:

1- Lorsque la vitesse V1 est superieur a la vitesse V2 (V1>V2),l'angle de

refraction devient inferieur a l'angle d'incidence i;

2-Lorsque V1<V2,l'angle a est toujours superieur a l'angle i;

3- Il y a un cas particulier lorsque V2>>V1 l'angle d'incidence i=il devient tel que
l'angle de refraction a=90°,on a:

sinil=V1/V2

Cet angle d'incidence il est nomme angle limite.Dans ce cas le rayon refracte
glisse le long d'interface;

50
 51

4- Si l'angle d'incidece est superieur a l'angle limite "i1" on a le phenomene de

reflexion total fig .En sismique,le cas exceptionnel ou le rayon refracte glisse le
long d'interface,represente un interet important.

IV 5-1 - Generalites sur les Hodochrones,sur les Champs de Temps

des Ondes Reflechies et Refractees;Ondes Coniques

Considerons deux formations homogenes I et II paralleles,dont les vitesses

sismiques sont V1 et V2 telles que V2>V1 fig ;au point O trois sortes d'ondes
pouront etre detectees:* Une onde qui atteint les points de reception S1,S2,...en
ligne droite ou onde directe;* Une onde reflechie par le toit de la formation II c'est
l'onde designee Rx. *Une onde refractee qui a glissee le long de l'interface I-II c'est
l'onde designee par Rf. Considerons,a present,les courbes "temps-distance".Sur
un graphe portons en abscisse les distances "x" separant le point d'enregistrement
S de la source d'ebranlement O et en ordonnee,les temps de propagation des
ondes sismiques.Nous obtiendrons trois courbes:

1-L'indicatrice (hodochrone) de temps de l'onde directe;

2-L'indicatrice de la forme d'une hyperbole-temps de l'onde reflechie Rx

3- La dromochronique de l'onde refractee Rf.

Etant donne que les ondes refractees ne peuvent etre presentees que si la
distance est superieur a la distance limite "x1".La dromochronique ne recoupe
jamais l'axe des temps.Remarquons que,pour cette distance limite "x1",les
hodochrones des ondes reflechies et refractees sont toujours tangentes.Il est a
noter qu'a partir d'une certaine distance,dite distance critique "xc",l'onde refractee
arrive aux sismographes avant les ondes directes.Nous avons note,que lorsqu'un
front d'onde rencontre une surface ou les proprietes physiques changent
brusquement,une partie de l'energie revient a la surface de la terre dans le premier
milieu,(reflexion),une autre partie est transmise dans le deuxieme milieu
(refraction).Le front d'onde,issue d'une source O ponctuelle,est une
sphere,jusqu'au moment ou l'onde incidente atteint l'interface.A partir de ce
moment un front spherique reflechi ayant pour centre l'image O' de O apparait fig
;et en meme temps du point A une onde refractee se propage dans le deuxieme

51
 52

milieu donr le front est une surface de revolution(pqrqllele a un ellipsoide) limitee

au meme cercle.Du fait que le temps s'ecoule et le point de reflexion s'eloigne du
point d'excitation,"i" augmente jusqu'a l'angle limite "i1" :

sinil=V1/V2 (a=p/2).

A partir de ce moment le front refracte etant devenu normal a l'interface,lui

demeure toujours normal.L'onde refractee glisse le long de l'interface avec la
vitesse V2;(V2>>V1).Des que l'onde atteint le point L,le front d'onde refracte
devient perpendiculaire a l'interface et se decroche,donc,de celui des ondes
incidentes (directes) Le front refracte introduit sur le trajet L un etat vibratoire(sur
l'interface),qui a son tour induit un etat vibratoire dans le milieu I.L'application du
principe de Huyguens permet de trouver le nouveau front conique d'onde:C'est
l'enveloppe de toutes les ondelettes crees dans le milieu superieur pqr le passage
du front refracte.Le front d'onde se propage dans ce milieu avec une vitesse qui ne
peut qu'etre egal a V1.D'autre part,ce front coupe la surface de la terre
(interface)sous l'angle limite "i1",c'est-a-dire la vitesse apparente de cette onde
sera egal a V2,puisque:

Varf=V1/sinil=V2.

Donc,nous voyons qu'il est possible,grace a l'onde conique de detecter,a partir de

la surface du milieu I,un ebranlement refractee par le milieu II.Pour cela,il suffit de
placer les geophones au dela d'une certaine distance de la source O;cette distance
est egale a x1: x1 s'appelle distance limite (ou encore la zone morte).Elle est
fonction de la valeur de l'angle limite "il" ( contrastes des vitesses ) ainsi que de
l'epaisseur "h" du milieu (superpose) superieur.

IV-6- Impedance Accoustique:Coefficient de Reflexion et de

Refraction

Lors de l'incidence d'une onde sismique sur une interface plane les ondes
secondaires (longitudinales et transversales) se forment.Les ondes secondaires de
meme type que l'onde incidente s'appellent "ondes monotypes" et celles d'autres
types "ondes d'echanges";l'intensite relative des ondes reflechies et refractees
depend de la repartition de l'energie de l'explosion,des conditions de reflexion et

52
 53

de la vitesse etc..Il a ete demontre qu'il y a les rapports suivants entre les
amplitudes des ondes secondaires et celle des ondes primaires (cas de l'incidence
d'une onde longitudinale sur la surface) :

- Al.sin(il) + At.cos(it) - Bl.sin(al)-Bt.cos(at) = - sin(i) ;

-Al.cos(il) - At.sin(lt) + Bl.cos(al) - Bt.sin(at) = cosi ;

-Al.sin(2al) +At.(V1/W1).cos(2it)+Bl.(r2/r1)(W2/W1)(V1/V2).sin(2al)+Bt.(r2/r1).(V1.W2/W12).cos(2at)= sin(2i) ;

- Al.cos(2il) - At.(W1/V1).sin(2it) - Bl.(r2/r1).(V2/V1).cos(2at) + Bt.(r2/r1).(W2/W1).sin(2at) = - cos(2i) .

ou :Al=al/a;At=at/a;Bl=bl/a;Bt=bt/a.

Les coefficients Al et At determinent les relations entre des amplitudes des ondes
reflechies (Longitidinale et transversale) et de l'onde incidente et s'appelle
"coefficient de reflexion".Physiquement,le coefficient de reflexion represente le
rapport de l'amplitude de l'onde reflechie a celle de l'onde incidente.Le coefficient
de reflexion est generalement de l'ordre du centieme au contact de diverses
formations geologiques,mais peut atteindre des valeurs beaucoup plus grandes en
presence de certains materiaux.A titre d'exemple les couches de charbon
reflechissent 30% de l'energie du signal incident. Les coefficients :Bl et Bt sont les
coefficients de refraction et montrent les relations entre des amplitudes des ondes
refractees (longitudinales et transversales) et de l'onde incidente.

IV-6-1 Analyse de Quelques Cas Particuliers:

IV-6-1-1 Incidence Normale:

Considerons tout d'abord le cas exceptionnel de l'incidence normale d'une onde

longitudinale sur la surface de separation de deux milieux.D'apres l'equation de
Snell nous avons:il=i'l=i't=al=at=0 . Dans ce cas les 4 equations precitees peuvent
etre ecrites de la maniere suivante:

Qt-Bt=0 ;
At+Bt.r2.W2/r1.W1=0 ;
-Al+Bl=1 ;

53
 54

-Al+Bl.r2.V2/r1.V1=-1 ;

d'ou,on a :

At=Bt=0 ;

Al = - (r2.V2/r1.V1).Bl+1
= - (1+Al).(r2.V2/r1.V1)+1
= - (r2.V2/r1.V1)-(r2.V2/r1.V1)Al+1

Al.(1+r2.V2/r1.V1) = (r1.V1 - r2.V2)r1.V1

Þ Al= (r2.V2 - r1.V1)/(r1.V1+ r2.V2)

Bl = 2r1.V1/(r1.V1+ r2.V2)

Donc,lors de l'incidence normale seules les ondes monotypes (logitudinales) se

forment,les ondes d'echange (transversales) etant inexistantes.On se souviendrea
que (la resistivite acoustiques) l'impedance acoutique d'un milieu donne est egal
au produit de la vitesse dans ce milieu par sa densite.Pour les ondes
longitudinales l'impedance acoustique est rV;pour les ondes transversales est
egale a rW.Designons par le terme "m" le rapport r2.V2/r1.V1=m.Alors pour les
coefficients de reflexion et de refraction nous avons les rapports suivants:

Al = (1-m)/(1+m)

Bl = 2/(1+m).

IV-6-2 Variation du Coefficient de Reflexion avec l'Angle d'Incidence:

Pour un angle d'incidence a' quelconque le rapport de la pression de l'onde

reflechie a celle de l'onde incidente est determine a l'aide du coefficient de
Rayleigh (Officier 1958)

R=[r2/r1 - (v12/V2-sin2i)1/2/(1-sin2a')1/2 ]/[r2/r1 + (v12/V2-sin2i)1/2/(1-sin2a')1/2 ]

54
 55

Lorsque le contraste d'impedance acoustique: r2.V2/r1.V1 est faible,seule une

petite partie de l'energie est reflechie par la surface de separation des deux
milieux.
Le coefficient de refraction aller-retour ( en incidence normale ) est :

B = 4.(r1.V1.r2.V2) / (r1V1+r2V2)2 = 1 - A21-2

L'amplitude A du signal reflechi par le second miroir est:

A = A2-3(1-A21-2).

Et pour le troisieme miroir on a:

A = A3-4(1-A21-2).(1-A22-3)

La perte par transmission n'est pas affectee par le signe du coefficient de

reflexion.En negligeant les pertes par divergence spherique et par absorption,il
ressort que l'amplitude d'une reflexion sur un miroir determine est le produit de son
propre coefficient par les produits des coefficients de transmission (aller-retour)
correspondant aux interfaces situees au dessus de ce miroir.Pour un coefficient
A1-2 assez significatif 0.4 par exemple,on peut obtenir une belle reflexion et donc
la perte par transmission associee a un seul reflecteur est
insignifiant.Toutefois,lorsque le nombre d'interfaces devient grand,la perte par
transmission devient plus importante.
1er cas: V2>V1:Lorsque la vitesse V2 est superieure a V1,le coefficient de
reflexion variant avec l'angle d'incidence croit a partir de l'incidence normale
jusqu'a une certaine valeur egale a l'unite a un angle d'incidence " i' " ; tel que : sin
( i ) = V1 / V2 Au dela de cet angle appele critique,le radical dans le nu;erateur du
second terme de l'expression :

T = 2.(r1.V1.r2.V2) / (r1V1+r2V2)

devient imaginaire.Ce radicale peut etre exprime par :

55
 56

± iÖsin2i - V12/V22.

La quantite : - iÖsin2i - V12/V22

est la representation correcte de ce probleme physique donnant a l'onde transmise
une amplitude qui diminue d'une maniere exponentielle dans le milieu sou-
jacent.Pour un angle d'incidence superieure a l'angle critique l'expression devient:

A=[r2/r1 + i(sin2i -V12/V22)1/2/(1-sin2i)1/2] / [r2/r1-i(sin2i-V12/V22)1/2/(1- sin2i)1/2

expression qui peut etre ecrite sous la forme:

A=eiE (Officier 1958)

ou:
E=2tg-1.r1(sin2i - V21/V22) / (r2(1 - sin2i)1/2.

2eme cas:V1>V2:

Dans ce cas le coefficient de reflexion decroit lorsque l'angle d'incidence croit a

partir de i=0 puisque la quantite:

(sin2i-V12/V22)1/2

decroit moins vite que :

( 1 - sin2i )1/2

Le coefficient de reflexion devient nul a une incidence telle que :

r2/r1 = (V12/V22 - sin2i)1/2/(1- sin2i)1/2

En resume;lorsque V2>V1,le coefficient de reflexion augmente d'une maniere

exponentielle a partir de i=0 jusqu'a une valeur egale a l'unite;lorsque l'angle
d'incidence est egal a l'angle critique.Dans le cas ou V1>V2,le coefficient de

56
 57

reflexion diminue de i=0 et devient nul a un angle d'incidence tel qu'il est defini par
l'equation ci-dessus.

IV-6-4 Coefficient de Reflexion et Proprietes Physiques des

Sediments:

IV-6-4-1- Coefficient de Reflexion et Diametre des Grains:

Les formules precedentes montrent que le coefficient de reflexion depend de la

densite et de la vitesse des formations geologiques.Dans le domaine de la
sismique marine,de nombreux travaux portant sur l'etude de la reflexion des ondes
sismiques sur les fonds marins,ont fait ressortir que le coefficient de reflexion
depend de la granilometrie des sediments qui constituent le fond marin.Ces memes
travaux montrent que le coefficient de reflexion croit avec le diametre des grains.La
synthese des travaux de divers auteurs(Mc.Kenzi - 1960 ; Jones et al -1964 ;
Zhitkovski et Lysanov - 1967 ; Song - 1977) fait ressortir que les coefficients de
reflexion proposes pour les fonds marins varient de 0.2 a 0.4 pour les vases
sableuses et les argiles;de 0.4 a 0.6 pour les sables vaseux et de 0.4 a 0.85 pour
les sables.Mis a part les conditions d'experimentation,la dispersion des valeurs du
coefficient de reflexion pour un meme type de sediment serait due a l'aspect
rugueux du fond de la mer:Car l'interface eau de mer-sediment se comporte en
realite non pas comme une discontinuite nette,mais comme une zone de transition
graduelle de l'impedance acoustique.

IV-6-4-2- Coefficient de Reflexion et Porosite:

Des travaux montrent que le coefficient de reflexion est lie lineairement a la

porosite des sediments marins.Engeneral,il diminue lorsque la porosite
augmente.(Faas - 1969) a defini une equation mettant en relation le coefficient de
reflexion a la porosite (f);elle s'exprime par:

f = 0.9795 - 1.4649.A.

IV-6-4-3- Coefficient de Reflexion et Frequence:

57
 58

La dependance entre le coefficient de reflexion et la frequence peut etre attribuee

a l'effet de stratification du fond marin,aux irregularites dans la topographie du fond
marin et a la dependance frequence (perte de frequence sur le fond
marin).Neanmoins cette dependance entre le coefficient de reflexion et la
frequence est differente suivant la bande de frequence utilisee.Zhitkovski et
Lysanov 1967;analysant les travaux realises par plusieurs auteurs,sont parvenus
aux conclusions suivantes:

1 - Ades basses frequences,les proprietes reflectrices du fond marin sont

determinees en premier lieu par les parametres physiques des couches
constituants le fond marin;l'epaisseur,la vitesse,l'absorption,la densite,la porosite.

2 - A des frequences superieures a quelques Khz,les proprietes reflectrices du

fond dependent dans la majorite des cas de la morphologie du fond,les parametres
acoustiques du fond marin n'ayant,dans ce cas aucun effet significatif.

IV-6-4-4- Mesure du Coefficient de Reflexion:

De nombreux auteurs ont tente de determiner les proprietes des milieux stratifies a
partir de l'information vehiculee par l'onde reflechie ou refractee.On utilise la
formule de base du coefficient de reflexion en faisant intervenir soit les amplitudes
soit les vitesses et les densites soit encore en utilisant les techniques de traitement
du signal.

IV-6-4-4-1 Incidence sur la Surface Libre (Surface de la terre):

Il est evident que la densite de l'air est egale a zero,si on la compare avec celle
des roches.Donc,on peut considerer la surface de la Terre comme une surface
libre sur laquelle les deplacements des particules ne sont limites par rien.Dans ce
cas les ondes refractees n'existent plus.Alors les coefficients de reflexion sont :

Al = (M - 1) / (M + 1)

At = 2r.ctg( 2it ) / (1 + M)

58
 59

ou:
M=tgi(l).ctg22it / tgi(t) ;
et
r = V1 / V2

Si l'onde longitudinale atteint perpendiculairement la surface libre,on donc: i (l) = 0

et A l = - 1 ; A t = 0 . Cela veut dire que,dans ce cas nous avons une onde
reflechie de la meme amplitude que l'onde incidente et l'onde transversale n'existe
plus.

IV-6-4-4-2 Incidence sur la Surface du Milieu de Regidite Absolue:

Le coefficient de reflexion peut etre positif et negatif selon le signe de la difference:

r1V1 - r2V2

Cela veut dire que la phase de l'onde reflechie peut etre la meme que celle de
l'onde incidente ou elles se different de 180°.Dans le cas considere les coefficients
de reflexion sont:

Al = cos(i+it)/cos(i-it)

At = - sin2i/cos(i-it)

Le signe (-) montre que dans ce cas il y a un changement de phase de l'onde.

IV-6-4-4-3 L'Angle d'Incidence est Petit:

Dans ce cas on peut ecrire:sini=i et cosi=i ; Les coefficients de reflexion et de

refraction pour les ondes longitudinales sont:

Al = (1 - M ) / ( M + 1 )

Bl = 2/(1+M)

59
 60

Ces formules coincident avec celles de l'incidence normale.Les coefficients de

reflexion et de refraction pour les ondes transversales sont :

At = Al.2i

Bt=Bl i(V2/V1 -1) = Bl i(r-1)

ou r = V2/V1

Autrement dit,lorsque l'angle d'incidence est petit les ondes transversales ont de
faibles intensites,mais augmentent proportionnellement avec l'accroissement de
l'angle d'incidence.Enfin,notons que l'intensite des ondes reflechies et refractees
depend de la difference(contraste)des impedances acoustiques.Si cette differences
est grande,le coefficient de reflexion (l'intensite des ondes reflechies) est aussi
grande et dans ce cas il s'agit de "la surface reflechissante forte".Ils peuvent se
presenter 3 cas:

1 er cas: r1V1 = r2V2 ; Al=0 ,c'est-a-dire l'onde reflechie n'existe pas.

2 eme cas: V1 ³ V2 ; l'onde refractee conique ne se forme pas.

3eme cas: r1V1 ¹ r2V2 ; l'intensite d'onde reflechie est autant plus grande que la
difference entre r1V1 et r2V2 est plus forte.

IV-6-5- Reflexion et Refraction des Ondes dans les Milieux a

Plusieurs Couches:

Dans un milieu stratifie sur chaque surface de separation des couches,les ondes
reflechies et refractees monotypes et d'echange prennent naissance.Nous avons
deja mentionne que si une onde incidente atteint un interface plan,elle donne
naissance a 4 ondes secondaires (2 ondes reflechies et 2 ondes refractees).Il est
evident que chacune de ces 4 ondes donnera naissance a quatre autres ondes
sur l'interface suivant et ainsi de suite.Les milieux a plusieurs couches grace aux
phenomenes de reflexions et de refractions sont les lieux de la dispersion
d'energie qui n'est d'ailleurs pas la seule cause d'attenuation de l'amplitude de

60
 61

l'onde incidente.Afin de mieux comprendre la repartition des energies des ondes

reflechies et refractees en fonction des differences d'impedances acoustiques des
couches et de leur nombre;considerons un milieu stratifie a trois couches fig
Examinons l'intensite de l'onde reflechie sur le troisieme interface lors de
l'incidence normale d'une onde longitudinale.Soit:

a¯ = Amplitude de l'onde incidente

b2¯ = Amplitude de l'onde passante dans la II couche
b3¯ = Amplitude dans de l'onde passanate dans la III couche
a3- = Amplitude de l'onde reflechie sur la III interface
b2- = Amplitude de l'onde passante d'en bas a II couche
b1- =Amplitude de l'onde passante d'en bas a I couche.

Il est evident que,l'amplitude de l'onde refractee traversant dans la IIeme couche

peut etre exprimee par la relation suivante:

b2¯ =a¯.B1¯

ou : B1 egale au coefficient de refraction sur l'interface I. On peut le trouver par

l'expression:

B1¯ = b2¯ /a¯ = 2.r1V1/(r1V1+ r2V2)

= 2 / ( 1 + r2V2 / r1V1)
= 2 / ( 1+a )

a = r2V2/r1V1.

De la meme maniere;le coefficient de refraction sur l'interface II peut etre determine

par la relation suivante:

B2¯=b3¯ /b2¯ = 2.r2V2/(r2V2+ r3V3)

= 2 / ( 1 + r3V3 / r2V2)
= 2 / ( 1+b )

b = r3V3/r2V2.

61
 62

Le coefficient de reflexion sur l'interface III est:

A3=a3-/b3¯ = (r3V3-r4V4)/(r3V3+r4V4)
= [1-r4V4/r3V3]/[1+r4V4/r3V3]
= (1-c)/(1+c)

c = r4V4/r3V3

Les coefficients de refraction sur l'interface II:Lorsque l'onde traverse de bas vers
le haut on:

B'2-= b2-/a3- = 2.r3V3/(r2V2+ r3V3)

= 2 / ( 1 + r2V2 / r3V3)
= 2 / ( 1+b' )

b'=r2V2 / r3V3 .

Le coefficient de refraction sur l'interface I :

B'1-=b1-/b2- = 2.r2V2/(r2V2+ r1V1)

= 2 / ( 1 + r1V1 / r2V2)
= 2 / ( 1+a' )

a' = 1/a=r1V1 / r2V2.

Le rapport entre des amplitudes de l'onde reflechie et de l'onde incidente s'exprime

comme suit:

b1-/a¯=B1¯B2¯A3 B'2-B'1-=A

donc le rapport A=b1-/a¯ est egal a:

A=2.2.2.2(1-c)/[(1+a)(1+b)(1+c)(1+b')(1+a')]

62
 63

remplacons (1+b')(1+a') = (1+1/b) (1+1/b) = [(1+a)(1+b)]/a.b

alors : A = 24ab(1-c)/[(1+a)2(1+b)2(1+c)2]

ou bien encore :

A = b1-/a¯ = 24(r3V3/r1V1)(1-(r4V4/r3V3))/[(1+(r2V2/r1V1))2(1+(r3V3/r2V2))2](1+(r4V4/r3V3))

On peut passer du cas de trois couches au cas general c'est-a-dire a plusieurs

couches.L'amplitude de l'onde reflechie par la derniere interface dans le milieu a
plusieurs couches peut etre determine par l'expression suivante:
b1-=A.a¯=a¯24(rkVk / r1V1)[(1 - (rk-1Vk-1 / rkVk))2] / Pnk =1[(1 + (rk+1Vk+1 / rkVk))2 ;

ou:Pnk=1 etant produit des (1 - (rk-1Vk-1/rkVk))2

pour k,variant entre 1 et n

IV-7 - Les Amplitudes en Sismique Reflexion:

Depuis un certain temps,les geophysiciens s'interessent a la detection directe des

hydrocarbures par l'analyse des amplitudes des ondes reflechies,les spectres des
ondes sismiques reflechies ou encore par la variation des vitesses sismiques.

IV-7-1 Etude des Amplitudes:

On sait depuis longtemps que les amplitudes des ondes sismiques varient ou plus
exactement,elles sont etroitement liees a la geologie du sous-sol.Une ondelette
emise dans le sous-sol est soumise,au cours de sa propagation a une serie de
facteurs qui l'affecte et que l'on peut resumer comme suit: 1- Facteurs independant
de la subsurface:-La geometrie de la nappe sismographique;-La variation de
l'energie emise par la source sismique;- L e couplage de la source et des
recepteurs sismographiques avec le sol;- La directivite des dispositifs d'emission et
de reception .

IV-7-2- Les Facteurs Dependants de la Subsurface :

(independant de la reflectivite des miroirs sismiques )

63
 64

- La divergence spherique ; - L'absorption ; - La focalisation et la defacolisation des

miroirs sismiques; - La difraction; - L a variation de la vitesse avec la profondeur; -
La distorsion de l'ondelette sismique.

IV-7-3- Les Facteurs Lies a la Reflectivite des Miroirs Sismiques:

- Le coefficient de reflexion;- Les pertes par transmission;- Les interferences; - La

conversion d'ondes; - La refraction.

IV-7-4- Les Differents Facteurs qui Affectent les Amplitudes de Reflexion: (Sheriff-
1973).

On peut egalement citer les facteurs tels que les bruits,l'attenuation due aux
appareils d'enregistrement. La trace sismique,la distance inter-geophones et la
longueur d'onde du signal.Une trace sismique comportant N geophones disposes
suivant une nappe lineaire,fournit une reponse ayant pour expression:

|R|=sin NpE/N.sin pE

ou E:represente une constante sans dimension il peut prendre des valeurs

comprises entre 0 et 1 : (0 < E < 1) ;

E = K.Dx = Dx / l = Dx/V.T = f.Dx /V = f.Dt

T:La periode ;Dx:Distance inter-geophone;l:Longueur d'onde V:Vitesse de

propagation.

IV-9 - Les Effets de la Transmission:

La transmissivite joue le role d'attenuateur;les geophones mesurent la vitesse du

deplassement des particules du sol.Lorsque l'impedance acoustique(produit de la
vitesse par le densite) augmente,la vitesse des particules diminue et l'amplitude du
signal devient:A(1-R) ou A etant l'amplitude du signal incident;R etant le coefficient
de reflexion.(1-R) represente la transmissivite pour la vitesse des particules;elle

64
 65

sera:(1-R)(1+R) sur le chemin du retour.D'ou a l'arrivee aux geophones

l'attenuation due a la trnsmissivite sera:

(1-R2).

Considerons le schema suivant:On a:

A2=(1-R1)A0; A3=R2(1-R1)A0 ; A4=R2(1-R2)(1+R1)A0 = R2(1-R22)A0.

Pour un milieu donne,la vitesse des particules est liee a l'energie emise par la
source sismique,mais elle depend surtout de l'impedance acoustique des
formations geologiques.L'effet attenuateur de la transmissivite peut etre
important,notamment sous une multi-couches renfermant du gaz,dans une serie
composee par une alternance de sel anhydrite ou encore lorsqu'il y a existance
d'une serie de bancs minces a tres forte impedance acoustique (cas des bancs de
dolomie).Il ressort que l'effet de la transmissivite a un effet attenuateur car,au fur et
a mesure que le front d'onde sismique atteint les differents reflecteurs,l'energie
qu'il vehicule decroit en fonction des coefficients de reflexion.

IV-9-1 - Distorsion de l'Ondelette:

L'ondelette emise,subit au cours de son trajet,des modifications.Elle subit une

distorsion d'amplitude et une distorsion de phase en meme temps.La meilleur
technique pour y pallir a ce phenomene est d'appliquer une deconvolution
notamment la deconvolution multicanal (Deconvolution surface consistante).Celle-
ci aura donc pour objectif d'effectuer la correction due a la distorsion de
l'ondelette.Cette technique consiste a calculer un filtre pour chaque collection du
point de tir (PT) ,Offset ,position sismo,et point miroir (CDP).Ainsi,les traces
constituant une meme collection,ont le meme filtre.Le filtre final a appliquer a
chaque trace est la multiplication des filtres relatifs au meme PT,meme off-
set,meme recepteur et meme CDP.

IV-10 - Principe de la Deconvolution Multicanal:

S1(t) . F1(t) . T1(t) .

S1(t) . F1(t) . T1(t) .

65
 66

S1(t) . F1(t) . T1(t) . ÅÞ T(t) ß

S1(t) . F1(t) . T1(t) . ß
. Ý
Sn(t) . Fn(t) . Tn(t) . ÆÜ
. D(t) .Ý

Sn(t) : represente la trace sismique correspondant au nieme canal.

Fn(t) : represente la reponse impulsionnelle du neme canal
Tn(t): represente la trace de sortie du neme canal avec: T(t)=S(t)*F(t).
La trace du monocanal de sortie est definie par: T(t) = åTn(t)
D(t): represente la trace de sortie desiree du canal
e(t): represente l'erreur au sens des moindres carres.
Dans le domaine discret,le filtre desire au sens des moindres carres est calcules a
partir de l'expression suivante:

jss.F = jsd

jss Autocorrelation multicanal des traces d'entreejsd ;cross-correlation entre les

traces d'entree et la trace sismique desiree ; F: le filtre desiree est calculee comme
suit:

F=jsd[ jss ]-1

L'application de la convolution a surface constante en traitement en amplitudes

reelles,permet une amelioration notammemt du point de vue
energie,frequence,continuite,coherence.Lorsque une onde sismique atteint une
surface de discontinuite(interface)separant deux couches geologiques
d'impedances acoustiques differentes,l'energie de l'onde incidente se divise en
deux parties:

- Une energie reflechie

- Une energie refractee (traverse l'interface)

Si Ar :l'impedance de l'onde reflechie

At : L'amplitude de l'onde refractee

66
 67

Ar et At sont liees a l'amplitude de l'onde incidente A0 par :

Ar = R.A0 ; At = (1 - R).A0 ;

ou :
R=(r2.V2 - r1.V1)/(r1.V1+ r2.V2) ; (cas d'une incidense normale) ;
r1.V1 et r2.V2 sont respectivement les densites et les vitesses de la premiere et la
seconde couche geologique
R:represente le coefficient de reflexion en incidence normale:

R = (V2.cosa1 - V1.cosa2 ) / (V2.cosa1 + V1.cosa2)

T = 2V2.cosa1 /(V2.cosa1 + V1.cosa2)

ou: R est le coefficient de reflexion a incidence inclinee

T est le coefficient de transmission en incidence inclinee

Plus le contraste des vitesses V1 et V2 est important plus la reflexion est

importante et le coefficient de reflexion est grand.Les variations lentes du
coefficient de reflexion impliqueront une variation lente des amplitudes. -

IV-10-1 - Coefficients de Reflexion - Calcules aux Divers Interfaces

Cas d'un Anticlinal: 1 - Couverture Argileuse ; 2 - Reservoir Sableux a Gaz ; 3 - Reservoir Sableux a Huile ; 4 -
Reservoir Sableux a Eau ; 5 - Couche Sableuse. Pour determiner les coefficients de reflexion sur les differents
interfaces,il est necessaire de connaitre les densites du sable envahi par le Gaz,l'huile et l'eau en utilisant la
formule de Wyllie. fig

IV-10-2 - Formule de Wyllie:

La formule de Wyllie permet d'aboutir au calcul de la vitesse apparente et la

densite apparente d'une couche poreuse :

dapp = df.f + dm..( 1 - f ) Va = Vf.Vm/ [Vf.+ f(Vm. - Vf )] .

67
 68

La relation entre la porosite et la vitesse est exprimee par la formule de Wyllie :

1/Vmo = ( 1 - f ) .1/Vm + Sw. f/Vw + (1- Sw ) f/Vhc .

Vm : Vitesse de la matrice ; Vf. : Vitesse du fluide ; f : Porosite ; df. : Densite du

fluides ; dm.: Densite de la matrice ; Vmo : Vitesse moyenne de la roche ; Vhc :
Vitesse des hydrocarbures ; Vw : Vitesse de l'eau ; Sw : Saturation.

A des coefficients de reflexion tres differentes correspondraient des amplitudes

tres differentes.

IV-11 - Influence de la Zone Alterree:

IV-11-1 - Correction d'Amplitude ( Surface Consistante)

Pour aboutir a la realisation des corrections d'amplitude "Surface Consistante" le

Sismicien doit faire certaines hypotheses notamment:- a) Reciprocite pour la
propagation de la source au miroir plan et de celui-ci vers le recepteur (on admet
que l'emetteur et le recepteur se situent a la meme distance du miroir plan). b) On
suppose une propagation directe de la source au capteur a travers un milieu que
l'on suppose homogene. c) On admet aussi que la zone alteree "WZ" a le meme
effet tant pour l'onde descendante que pour l'onde montante; d) On admet
egalement que la reponse de la "WZ", le couplage de la source et des recepteurs
(y compris leur sensibilite) ainsi que l'attenuation provoquee par la "WZ",sont
constants durant le temps d'enregistrement ; e) L'energie emise par la source
affecte l'ensemble des traces enregistrees d'une maniere analogue; f) Le couplage
recepteur est analogue sur toutes les traces sismiques pour differentes positions
de la source d'emission; g) Les traces sismiques correspondant a un meme
"CDP",renferment la meme information du sous sol.
On admettant ces hypotheses,il est possible de separer les facteurs de subsurface
consistante pour les facteurs les plus importantes soient :

a) - EPT(f) : Reponse frequentielle de la source;

b) - GS(f) : Reponse frequentielle du recepteur;
c) - CSS (f) : La reponse frequentielle du sous-sol,c'est le facteur porteur
de l'information utile (geologique utile) ;

68
 69

d) - O(f): Reponse frequentielle a la position Offset;

e) - CD(f): La reponse frequentielle du canal d'enregistrement.

En prenant en compte les hypotheses simplistes ci-dessus,on peut definir la trace

sismique dans le domaine frequentielle par l'expression suivante:

TS(f) = EPT(f).GS(f).CSS (f).O(f).CD(f).

Le but de la prospection sismique est de determiner "CSS (f)";Les influences de

"EPT(f);GS(f);O(f) , sur le signal utile peuvent etre suprimees. En general,on obtient
un systeme d'equation dont le nombre d'inconnus est inferieur au nombre
d'equations.La resolution peut etre effectuee par differentes methodes (Methode
de Gauss-Seidell).En prenant le Logarithme neperien,on obtient:

Ln TS(f) =Ln EPT(f) + Ln GS(f) + Ln CSS (f) + Ln O(f)

Avec: TS(f) =| TS(f) | ej(f) ; Ln T S(f) =| TS(f) | + ij(f) ;

ou :| TS(f) | : Represente le spectre d'amplitude de la trace sismique;

j(f) : L e spectre de Phase de la trace sismique.

Pendant la sequence de traitement,les deux systemes sont determines

separemment:Le premier relatif aux spectre d'amplitudemse rapporte au traitement
en amplitude surface consistante tandis que le deuxieme systeme d'equations
relatif a la phase a pour but de realiser le calcul des corrections statiques
residuelles.

IV-11-2 - Amortissement des Ondes Sismiques:

Les ondes reflechies et refractees,arrivant a la surface,ont une tres faible

energie.Les deplacements des particules du milieu ne depassent pas quelques
microns.L'intensite des ondes enregistrees a une grande importance pour les
methodes sismiques.On sait que dans un milieu homogene et isotrope,l'intensite
des ondes spheriques diminue inversement avec le trajet parcouru a cause de la
dilatation.Cette loi est vraie pour les ondes reflechies,mais l'intensite des ondes
refractees grace a une autre forme du front (conique,et a une grande distance"x"

69
 70

entre le point de tir et les sismos) est inversement proportionnelle a x2.Il en resulte
que l'intensite des ondes refractees diminue plus vite que celle des ondes
reflechies.Et ce,en fonction du trajet parcouru.Theoriquement l'amplitude d'une
onde en tout point du milieu peut etre calculee d'apres l'expression suivante:

AX = A0 e-ax

ou: A0 : est l'amplitude initiale ; a: etant le coefficient d'absorption; x: distance.

Sachant que : -ax = ln e-ax

On trouve: a= ln (AX / A0)/x; [m-1]

Le coefficient d'absorption depend de la frequence (f) des oscillations:

Pour les roches solides : aS = b.f

pour les roches meubles :am = n.f 2

ou : b et n sont les coefficients de proportionnalite.L'absorption des roches

meubles est plus importantes.Les coefficients :am :depend de la lithologie et de
l'humidite des roches;am ,est d'autant plus petit que l'humidite est grande.En tenant
compte des pertes par absorption,on peut conclure,que les ondes utiles qui
arrivent aux sismographes en provenance des grandes profondeurs possedent une
tres faible intensite par comparaison avec celles venant de faible profondeur.

IV-12 - Diffraction des Ondes:

Si une Onde Sismique rencontre un obstacle dont les dimensions sont de l'ordre
de la longueur d'onde de propagation,cette onde sera perturbee.Le phenomene est
appele la diffraction.Lorsqu'une onde sismique atteint une proeminence dans le
milieu geologique,par exemple le compartiment superieur d'une faible couche,la
diffraction pourrait perturber gravement les enregistrements.Cette proeminence
joue le role d'une nouvelle source,donnant naissance aux ondes qui se propagent
dans toutes les directions du point de diffraction.On peut construire des isochrones
de l'onde diffractee suivant le principe de Huyguens,en considerant le point de

70
 71

diffraction comme une nouvelle source elementaire.En sismique,la notion de

diffraction est associee a celle de Faille,flancs de pente raide,de
monoclinaux,synclinaux et anticlinaux.Sur une section sismique,une diffraction se
traduit par des evenements pentes.

IV-13 - Ondes Utiles et Parasites:

Pour bien comprendre la repartition (division) des ondes: en ondes utiles et ondes
parasites,nous donnons ici une breve discription de la sismique sans entrer dans le
detail:On observe les ondes elastiques qui reviennent a la surface du sol apres
s'etre propagees dans le sous-sol a la suite d'un ebranlement.Dans le sous-sol,les
discontinuites des constantes elastiques ou des densites entrainent des reflexions
,refractions ou diffractions de l'onde provenant de l'ebranlement.Leur etude doit
nous reveler la position des discontinuites,les vitesses de propagation et les
coefficients d'absorption.

CHAPITRE V ACQUISITION DES DONNEES SISMIQUES

V-1 Introduction :

Generalement,on procede a une explosion de dynamite dans les couches

superficielles du terrain.Des sismographes disposes en surface sont relies par des
cables electriques a des modules d'amplification et d'enregistrement a plusieurs
canaux;Dans les 4 ou 5 secondes qui suivent un tir,on enregistre sur bandes
magnetiques ou disque les retours d'energie de l'explosion sous forme de
signaux.Les groupes de sismographes disposes en lignes par rapport au ligne de
l'explosion donnent naissance a des traces (de 24;48 ;96 sismos).On applique
ensuite a ces traces des corrections de temps correspondant aux differences des

71
 72

temps de propagations du point de tir aux sismographes,ainsi au'aux differences

eventuelles d'altitude des sismographes.La detection des reflexions est ensuite
faite par des sismiciens qui identifient a vue les reflexions sur l'enregistrement
photographique (sismogrammes) par des analogies d'amplitude,de forme et de
mise en phase constatees sur toutes les traces et permettant meme de suivre le
profil d'une couche sur les sismogrammes des tirs voisins.

V-2- Sismographe ou Geophone:

V-2-1 Introduction:

En Prospection sismique,le sismographe terrestre ou geophone est l'appareil

charge de traduire en tension electrique le mouvement du sol,le plus souvent sa
composante verticale.(le signal electrique varie dans les limites suivantes : 1/10
mV a 1 Volt).
Il est suppose traduire fidelement le mouvement du sol,c'est-a-dire s'assurer que le
mouvement que prend le boitier du geophone en contact avec le sol est le meme
que celui qu'aurait le sol en son absence,cela se traduit egalement par la fidelite
du couplage sol-geophone dans la bande des signaux sismiques.
Le plus souvent la trace sismique n'est pas la reponse d'un sel geophone,mais
celle d'un groupement de geophone permettant des l'enregistrement une
amelioration du rapport signal sur bruit,par filtrage en nombre d'onde.
Il existe plusieurs modeles de sismographes bases sur des principes differents et
sont choisis a des buts divers.On distingue les capteurs suivants:
a ) - Capteurs a ceramique,reposant sur l'effet piezo-electrique.
b ) - Sismographes a pression (utilises pour les carottages grandes sondes.)
c ) - Sismographes electromagnetique a flux coupe.

V-2-2 Etude d'un Sismographe:(Sismographe electromagnetique a flux coupe)

Le Sismographe electromagnetique a flux coupe est base sur les lois de l'induction
magnetique et de la force electromotrice d'induction (F.E.M.).En effet ,si on
approche un aimant permanent "a" d'une bobine "b",reliee a un millivoltmetre a "0"
central. A chaque mouvement de "a" on note le passage d'un courant "i" dans le
sens indique par le millivoltmetre. [Ce passage de courant existe seulement
lorsque le mouvement existe] Si l'aimant permanent reste fixe mais que la bobine

72
 73

se deplace par rapport a celui-ci il existe egalement un courant "i"; ce courant "i"
est appele courant induit.Dans tous les cas il apparait une F.E.M. induite
puisque,un courant a circule dans cette bobine.C'est la loi de Faraday qu'on peut
rappeler ici:Toute variation de flux a travers un circuit electrique donne naissance a
un F.E.M induite.De nombreux models de sismographes sont une application
directe de cette loi.

V-2-3 Description d'un Sismographe:

(Sismographe electromagnetique a flux coupe)

Soit un barreau aimante solidaire du corps d'un sismographe,et par consequent du

sol sur lequel celui-ci est pose.Autour de cet aimant on a place une bobine
suspendue par 2 "spieders" (inferieur et superieur).Au moment de l'arrivee d'une
onde sismique ou d'une vibration,le sismographe(solidaire du sol) se deplace par
rapport a la bobine qui par inertie a tendance a rester immobile.Il se forme donc un
courant induit dans cette bobine.C'est ce courant induit qui est enregistre sur des
supports (bandes magnetiques ou disque).L'ensemble aimant-bobine-spieders est
enferme dans un boitier etanche.Sur ce boitier on a place 2 bornes a souder
isolees de la masse par des perles de verre.Les 2 bornes etant polarisees
permettent le branchement exterieur du cable sur le sismographe.Ce boitier se
trouve place a son tour dans une enceinte en matiere plastique compose de 2
parties;une partie superieure formant le couvercle qui est fixee au corps inferieur
par deux baques emmanchees;de force ou visees.Une partie inferieure qui recoit le
sismographe;enfin une pointe amovible vissee sur le corps du sismographe facilite
l'implantation de celui-ci dans des terrains de nature differente.

V-2-4 Multiplication des Sismographes ou Geophones:

Quand on observe un enregistrement sismique,on s'apercoit que les signaux utiles

sont plus ou moins noyes dans des bruits parasites.On constate egalement que
ces signaux sont dephases d'un sismographe a l'autre.La premiere idee consiste
alors a augmenter le nombre de geophones,afin d'ameliorer statistiquement la
qualite du signal;le rendement etait tres faible;mais ce qui compte le plus c'est
l'idee,car elle a ouvert la voie
plus tard a une meilleure etude du signal sur bruit.Donc pour ameliorer le rapport
signal/bruit on effectue une multiplication de geophones.

73
 74

Cette multiplication peut s'effectuer soit a l'emission,soit a la reception ou aux deux

a la fois.A l'emission c'est une multiplication de charges et a la reception c'est une
multiplication de geophones.
Pour ameliorer la mise en phase des signaux de deux traces consecutives on
effectue ce que l'on appelle un Mixage. fig

V - 3 - Aspect Aléatoire des Sismogrammes:

En sismique,les ondes sont nommees "utiles" si on les utilise pour l'etude de la

coupe geologique.Par exemple:En sismique - reflexion;les ondes reflechies sont
des ondes utiles,toutes les autres ondes sont considerees comme parasites
(bruits).En sismique refraction ce sont les ondes refractees qui sont utiles;alors que
les ondes reflechies deviennent des parasites,qu'il faut les
éliminer.Theoriquement,un sismogramme a un contenu certain,representant
l'ensemble des reflexions consecutives a l'explosion provenant de couches
geologiques d'emplacement determine.Mais a priori ces emplacements sont
inconnus.De plus,si l'on effectue deux tirs au meme endroit,les deux
sismogrammes different notablement,du fait de nombreux facteurs aleatoires.Les
impulsions S(t) engendrees par les explosions ne sont pas identiques.Elles
subissent une deformation et un allongement avec la profondeur et l'absorption.

V-4 - Classification des Ondes Parasites:

En prospection sismique reflexion,on designe par "bruit",tout ce qui

perturbe,oblitere ou cache l'information sismique;c'est-a-dire,tout ce qui ne se
rapporte pas a un signal reflechi (dans le cas de la sismique reflexion bien
entendu).Donc seules les ondes reflechies sont utiles;les autres doivent etre
eliminees dans la mesure ou elles se distinguent par leurs frequences,leur mode
de propagation etc...La frequence des ondes reflechies est en general,plus grande
que celle des ondes parasites,mais cette regle n'est pas absolue.Le bruit constitue
l'ensemble des oscillations qui sur le "film sismique" ou (monitor) tendent a
masquer les ondes utiles qui nous interessent.
On classe les ondes Parasites en deux Groupes :

1 - Les ondes parasites Horizontaux

74
 75

2 - Les ondes Parasites Verticaux.

Les ondes Parasites Horizontaux sont:

- Les Ondes de Rayleigh

- Les Ondes de Love
- Les Ondes Directes
- Les Ondes Aeriennes
- L'onde Aerienne Couplee a l'Onde de Rayleigh
- Le Bruit Industriel , Naturel et Ambiant
- Le Bruit Aleatoire.

Les ondes Parasites Verticaux sont:

- Les Ondes Multiples

- Les Ondes Fontomes
- Les Reverberations
- Les Dffractions.

V-4-1 - Les Parasites Horizontaux:

V-4-1-1- Les Ondes de Rayleigh:

L'onde de Rayleigh fait partie des ondes parasites appelees "Ground-roll".Ce sont
des ondes superficielles qui prennent naissance a proximite de l'air libre.L'onde de
Rayleigh represente un mellange de l'onde longitudinale et de celle transversale
dont le mouvement a une composante verticale et une composante horizontale
parallele a sa propagation.Son mouvement est elliptique et retrograde fig .Elle
constitue le principal bruit organise que l'on enregistre en prospection sismique et
que les sismiciens appellent Ground-roll.Les geophones sont generalement
implantes verticalement ce qui ne permet pas d'eviter l'enregistrement de sa
composante verticale.Les ondes de Rayleigh ont ete etudiees depuis 1885,et
depuis cette date des progres ont ete realises dans leur etude.Les ondes de
Rayleigh se propagent dans la partie superficielle de la croute terrestre (quelques
kilometres),que l'on assimile a une couche de faible Vitesse.L'evolution des
vitesses des ondes de Rayleigh en fonction de la longueur d'onde,de la

75
 76

frequence;le nombre d'onde est etroitement liee aux vitesses des couches
superficielles,leurs epaisseurs,des parametres elastiques,du rapport des vitesses
des ondes transversales,ainsi que le rapport des densites etc...
Cas d'une seule couche:Si on suppose que la couche a une grande epaisseur
comparee a la longueur d'onde,Rayleigh a montre que cette onde est non
dispersive et se propage dans la couche superficielle avec une vitesse constante
equivalente a : 0.9194 de l'onde n des transversale de la couche superficielle.

VR = 0.9194 VS (Pour un coefficient de Poisson = 0.25)

Si on designe par UX,la composante horizontale du mouvement des ondes de

Rayleigh,celle-ci aura pour expression:

UX=D.[e-0.8475.K.H - 0.5773.e (-0.3933.K.H) ].sinK(Vt-x).

La composante verticale s'exprime par :

UR=D.[-0.8475e-0.8475.K.H + 1.4679.e(-0.3933.K.H) ].cosK(Vt-x).

V:vitesse;t:temps;X:la direction horizontale de la source;H:la profondeur;K:le

nombre d'onde;D:une fonction de K .Les ondes de Rayleigh se caracterisent par un
mouvement elliptique et retrograde ,on determine que sa composante verticale est
environ 1.5 fois plus grande que sa composante horizontale.Cette derniere
composante disparait a une profondeur de 0.192l;ou l est sa longueur
d'onde.Dans le cas de deux couches superficielles,des etudes font ressortir que les
ondes de Rayleigh sont dispersives.

V-4-1-2 - L'Onde de Love:

L'onde de Love a un mouvement horizontal par rapport a la direction de

propagation Elle n'est pas generalement enregistree par les geophones poses
verticalement sur la surface du sol.C'est donc une onde purement transversale.Elle
est liee a la frequence de la zone alteree dans la mesure ou cette derniere est
limitee par la discontinuite Sol/Air d'une part et une couche a caracteristiques
elastiques differentes sur laquelle elle repose L'equation de l'onde de "LOVE" est
donnee par la relation:

76
 77

tg(KH.ÖV2Phase/V2Cisail-(1)-1)=m2/m1Ö(1-V2Love/V2Cisail-(2))/Ö(V2Love/ V2Cisail-(1)-1)

ou :
m1 et m2 : represente respectivement les coefficients de LAME ou modules de
cisaillement des milieux (1) et (2)
H:epaisseur de la couche superficielle
VCisaillement(1) et VCisaillement(2) : represente la vitesse de propagation de l'onde de
cisaillement dans les milieux (1) et (2)
VPhase :Vitesse de phase de l'onde de LOVE en fonction de la frequence

V-4-1-3 - Les Ondes Fantomes:

Ils prennent naissance dans deux cas:1er cas:Lorsque la source sismique se

trouve dans la zone alteree,caracterisee par une grande epaisseur; 2eme
cas:Lorsque la source sismique est disposee au dessous de la base d'une zone
alteree composee de deux couches dont la premiere est caracterisee par une
faible Vitesse.

V-4-1-4 - Les Ondes Guidees:

Les Ondes guidees prennent naissance,la ou les ondes longitudinales sont

coincees dans les premieres couches superficielles a cause du grand contraste
entre le sol et l'air et les couches alterees et les couches situees au dessous . Les
ondes guidees se caracterisent par des vitesses plus grandes que les vitesses des
ondes de Rayleigh et les Ondes directes.

V-4-1-5 - Les Bruits Naturels:

Ces bruits sont caracterises par la frequence un peu superieur a 50 Hz,due a

l'agitation existante a la surface du sol (microseismes).
Les microseismes sont les ebranlements desordonnes du sol,formes par de
differentes causes exterieures tels que:
vents,pluies,circulation,transports,et bruits industriels etc...
Ils etouffent surtout les signaux faibles de la fin d'un enregistrement.

77
 78

V-4-1-6 - Les Bruits Provoques par le Tir:

Ils contiennent des reflexions multiples,des diffractions,des ondes de surface,onde

directe,ondes refractee par la base de la WZ,ondes aeriennes,etc...

V-4-1-6-1 - L'onde Directe:

Ce sont des ondes longitudinales qui se propagent dans la premiere couche de la

coupe geologique sans aucune reflexion et refraction;elle est toujours enregistree
en premier au voisinage du point de tir.

V-4-1-6-2 - L'onde Refractee:

C'est l'onde refractee par la base du terrain altere ou de la WZ . Au voisinage du

point de Tir on l'enregistre apres l'onde directe,mais elle arrive aux sismographes
les plus eloignes avant l'onde directe.

V-4-1-6-3 - L'onde Aerienne:

En cas de debourrage du trou ou de Tir en Surface,l'explosion provoque un souffle

qui apparqitra sur les enregistrements avec une forte amplitude et une haute
frequence (25 - 35 Hz) .Cette onde se propage dans l'air avec une vitesse de 340
m/s environ.L'onde aerienne apparait sur les films de production ou sur l'ecran par
un cone ,appele cone de l'onde aerienne.Celle-ci presente un inconvenient assez
significatif sur les traces sismiques proches du point de Tir.

V-4-1-6-4 - Les Reflexions Multiples :

Tres frequement,on observe sur les enregistrements des réflexions semblant

prévenir d'un reflecteur situe a une certaine profondeur sous un premier
miroir,mais qui ne sont,en effet,que des phenomenes de repetitions a qui l'on
donne le nom de multiples.On peut affirmer que les reflexions multiples sont
presques toujours presentes et qu'elles constituent un probleme insidieux et parfois
difficile a resoudre.Sur les sections sismiques,les multiples genantes,se presentent

78
 79

generalement avec une courbure plus forte que les reflexions reelles.Elles restent
donc,pentees sur la section corrigee.

V-4-1-6-5 - Ondes diffractees:

Ils prennent naissance sur des failles,des corbures ou rupture des couches,dome
de sel etc...,en un mot,par des proemineuses du milieu,elles peuvent perturber les
enregistrement sismiques.

V-4-1-6-6 - Des Oscillations Parasites:

Dues a un mauvais couplage du sismgraphe au sol,ou des reglages mal faits,ainsi

que le bruit de l'appareillage brouillent enfin le signal.Donc,les ondes utiles et
parasites arrivent a la surface de la terre presque simultanement interferent et
defiguent gravement les enregistrements. Pour eviter l'influence des ondes
parasites sur les signaux utiles et ameliorer le rapport signal sur bruit;on applique
tres frequement les differentes techniques:Appareillage et systeme d'observation
(multiplication des geophones et des charges , couverture multiple , les differents
filtres,etc...Enfin,preciser que la notion de "signal" et "bruit" sont tres utilises en
sismique:Une arrivee d'energie,avec son caractere constitue ce que l'on nomme un
signal.L'ensemble des phenomenes parasites est qualifie de bruit,on distingue les
bruits organises et non organises.

CHAPITRE VI - PARAMETRES " GEO-SISMIQUES "

VI-1 Introduction:

79
 80

Les roches composants l'ecorce terrestre sont tres differentes et possedent

defferentes compositions petrographiques,voila pourquoi elles ont differentes
proprietes elastiques. Certaines roches,par exemple,les roches eruptives et
sedimentaires fortement consolidees se rapprochent des corps parfaitement
elastiques,mais des formations geologiques des zones de fissures,ainsi que les
sediments meubles different notablement de milieux parfaitement elastiques;elles
ont des coefficients d'absorption tres grandes.

VI-1-1 Notions sur les Vitesses Sismiques:

Depuis plus de 50 ans l'exploration sismique dans le domaine de la recherche des

hydrocarbures a ete la methode la plus utilisee en raison des immenses
possibimites qu'elle offre dans l'interpretation des donnees.De plus,des
ameliorations constantes sont apportees a cette methode qui est de nos
jours,parvenue a un tres haut niveau de perfectionnement. Toutefois,jusqu'a ces
dernieres annees,l'objectif essentiel des recherches consistait a determiner
l'existence des pieges structuraux tels que anticlinaux,domesmfailles,etc...
Les methodes Sismiques reflexion sont basees sur la mesure des temps d'arrivee
des ondes elastiques engendrees par une source d'ebranlement adequate et
reflechie sur les differents niveaux geologiques.Ces temps d'arrivee en provenance
des differents horizons reflecteurs dependent de la vitesse moyenne a laquelle se
deplace l'onde sismique a travers les differentes formations geologiques.Si l'on
desire supprimer l'effet d'obliquite " NMO" resultant des variations de la distance
separant le point de Tir du point de detection et determiner la profondeur des
principaux horizons (conversion - temps-profondeur).Il faut connaitre la loi de
vitesse qui regit la propagation des ondes sismiques dans le sous-sol.
Les recents progres realises dans l'amelioration des enregistrements et dans le
traitement des donnees sismiques ont conduit a une interpretation specifique liee a
la notion de niveaux d'energie.Ainsi,les roches sedimentaires ont pu etre classees
en depots de forte et de faible energie,ce qui permet d'associer les facies
sismiques aux facies sedimentaires.L'etude detaillee des vitesses de propagation
des ondes sismiques, combinees a d'autres parametres physiques permet de
recueillir des informations de nature lithologique,de rechercher des proprietes
physiques des roches (densites,porosite contenus en fluides etc..),de determiner
quantitativement le pourcentage sable-argile,de detecter les zones sous-
composees etc...De telles etudes detaillees permettent,dans une certaine

80
 81

mesure,de voir les variations rapides de facies et ainsi faciliter la decouverte des
pieges lithologiques et stratigraphiques comme les recifs,les biseaux,les lentilles
sableuses,les chenaux enterres, etc...

VI-1-2 Expressions Mathematiques des Vitesses des Ondes

Elastiques dans les Milieux Solides non Poreux

Le principe de l'etude de la propagation des ondes elastiques dans un milieu

homogene et isotrope derive de la loi de Newton:

F = m.g

ou : F : Forces exterieures qui engendrent des tensions

g : l'acceleration,elle represente la derivee seconde des deformations.
Les deformations permettent d'etablir les equations fondamentales de propagation
des ondes dans un milieu elastique.Pour une onde longitudinale,l'equation de
propagation s'ecrit:

(l+2m).¶2u / ¶x2 - z/g .¶2u / ¶t2 = 0

ou l et m: les coefficients de Lame ;z: poids specifique apparent du milieu ; g:

acceleration de la pesanteur.
La theorie d'elasticite (voir chapitre I),nous enseigne que pour un milieu homogene
et isotrope,la vitesse des ondes longitudinales Vp; et transversales Vs ;sont liees
aux constantes elastiques et a la densite par les relations suivantes:

Vp = | E(1-d) / d . [(1+d)(1-2d)|1/2 = |(k+4/3.m)|1/2

Vs = | E/d.(1+d) / 2|1/2 =(m/d)1/2

VI-1-3 - Vitesses Rencontrees en Sismique:

On denombre 5 types de Vitesses en prospection Sismique:

1- La Vitesse instantanee Vi
2 - La Vitesse Moyenne Vm

81
 82

3 -La vitesse par Tranche Vtr

4 - La Vitesse Quadratique Vq.m
5 - La Vitesse de Sommation VS.

VI-1-3-1- Vitesse Instantanee:Vs

C'est la vitesse du front d'onde qui se deplace dans la direction du rayon

sismique.Elle est la vitesse reelle de propagation des ondes sismiques
Pratiquement,il est impossible de la determiner,car on ne connait jamais
exactement la direction du rayon sismique (a cause des heterogenietes du milieu
geologique).

VI-1-3-2- Vitesse Moyenne:Vm

La vitesse moyenne peut etre determinee a partir d'un carottage,en se basant sur
les donnees d'une analyse.Pour un milieu horizontalement stratifie avec des
couches d'epaisseurs hi et des vitesses vi ;on peut ecrire l'equation suivante pour
la vitesse moyenne Vm :

Vm = å hi / å ti =å hi / å hi /vi = 2å Pi/ T0.

ou Pi : profondeur d'intervalle ; T0 : temps de parcours double. Pour un milieu

continue,il faut remplacer le signe somme "å "dans la formule par le signe integrale
"ò" on obtient la formule suivante donnant la vitesse moyenne en fonction de la
profondeur "z" :

Vm(Z) = Z / tv = Z / ò dZ / Vs(Z) .

En sismique on construit le plus souvent la courbe:

Vm = f(tv)

ou: tv : etant le temps vertical. Enfin,on peut calculer la vitesse instantanee Vs(Z) si
l'on connait la vitesse moyenne Vm(Z) :

82
 83

Vs(Z) =V2m(Z) / (Vm(Z) -Z.d Vm(Z)/dZ ).

VI-1-3-3 Vitesse d'Intervalle :Vi(P)

C'est une vitesse qui correspond a un intervalle geologique determine.Connaissant

les temps " tn e t tn-1 " correspondant aux profondeurs donnees " pn e t pn-1 " :
On calcule la vitesse d'intervalle (on l'appelle encore vitesse par tranche) :

Vi(P) = pn - pn-1 / tn - tn-1 ;Vi(P) = Pi / ti .

VI-1-3-4 Vitesse Quadratique :VRMS

Cette vitesse est calculee d'apres les "Indicatrices des ondes".

VRMS= ÖåV2i .t2i / å t2i .

VI-1-3-5- Vitesse par Tranche :Vp.tr .

Si la vitesse moyenne est a peu pres constante dans une couche qulconque,alors
on la considere comme etant tres proche de la vitesse instantanee et on presente
cette couche comme une couche homogene et,enfin,on l'appelle "tranche
sismique".Cette vitesse moyenne est nommee "vitesse par tranche" Vp.tr Elle
caracterise les proprietes d'une couche sismique. fig .

VI -1-4 La Mesure des differents types de Vitesses.

Il y a plusieurs methodes pour mesurer les vitesses sismiques ,on les divise en
deux grands ensembles.

-1 Les methodes Directes: (mesures dans les puits)

83
 84

Elles se mesurent a l'aide des diagraphies acoustiques (carottage continu de

vitesse) ou a l'aide du carottage sismique (sismosondage) et des pseudologs
d'impedance acoustiques.

-2 Les Methodes Indirectes:

Elles se realisent a partir des methodes statistiques ou a partir des tirs speciaux
comme par exemple :

la methode du : "X2 - T2"

La methode des: "T - Dt "

La methode de la meilleure coincidence etc... et enfin a partir des films par

traitement numerique des donnees sismiques du type "Analyse de Vitesse" sur
ordinateur.

Les Vitesses mesures par carottage sismique et carottage continu de vitesse

necessitent un forage;on mesure les Vitesse suivantes:

VI - 1 - 4-1 La Vitesse Instantanee:

C'est la vitesse d'une couche geologique situee a une profondeur "z" lorsque son
epaisseur tend ( Þ 0 ) vers zero ;

Vi = dz/dt

ou dt: etant le temps simple de propagation du trajet "dz".La vitesse instantanee se

mesure a l'aide du carottage continu de vitesse (CCV) au moyen d'un outil
denomme "sonde sonique" dont le principe est schematise sur la figure .Cette
outil consiste a mesurer d'une maniere continue la valeur de la vitesse instantanee
(Dt = 1/DVi ) le long du trou de forage. Les plus recents outils soniques comportent
deux emetteurs "E1 et E2 " d'impulsions sonores de hautes frequences (20 a 40
KHz) et 4 recepteurs " R1,R2,R3,R4" .Les ondes sonores transmises par le recepteur
E1 sont enregistrees par R1,R3 et celles emise par E2 sont enregistrees par R2,R4 .

84
 85

Cette technique permet de reduire les erreurs par le calcul de la moyenne des
temps entre les recepteurs R1 et R3 et les recepteurs R2 et R4 .Les temps des
trajets mesures sont exprimes en "m.sec / pied".Les resultats de mesure de "t" sont
representes sur un film (log sonique).Les temps de parcours des ondes sonores
dans les formations geologiques sont lies aux vitesses instantanees par la relation:

V(ft / sec) = 106 / Dt ("m.sec / pied").

Les methodes acoustiques permettent de mesurer no seulement le parametre

instantane,mais aussi de determiner d'autres parametres ou de mettre en evidence
certaines caracteres des milieux tels que:
- Attenuation des ondes sonores dans les formations geologiques;
- Obtention des informations sur la porosite;
- Identification de la lithologie et de la zone a pression anormale;
- Mesures des proprietes elastiques des roches;
- Detection des fractures et des cavites des roches ;
- Controler le casing;
- Estimation de la consolidation des roches De plus , outre l'interet du "CCV" pour
le calage des donnees sismiques,il est encore utilise pour realiser les films
synthetiques et effectuer des correlations entre les differents puits.

VI -1- 4-2 La Vitesse Moyenne :

Le rapport "Z / T", ou Z: represente la profondeur et T : le temps simple de

propagation correspondant.Ce rapport s'appelle "Vitesse Moyenne".

Vm = Z / T = [ z1+z2+z3+z4+......Zn] / [T1 +T2+T3+T4 ......Tn]

La vitesse moyenne peut etre egalement definie comme la moyenne temporelle

des vitesses instantanees Vi par l'expression:

Vm = 1/T0 òVi(t)dt = 1/TSVk.tk

ou T0 =Stk = 2SHk.Vk
Hk : epaisseur de la couche "k"

85
 86

T0 : temps double a la verticale

La vitesse moyenne verticale se mesure uniquement dans les forages par la

technique du sismosondage.Elle est utilisee,en prospection sismique,dans les
differents processus de l'interpretation geologique et non pour realiser les
corrections "dynamiques",car elle donnerait des valeurs de Dt trop grandes.Le
principe de mesure de la vitesse moyenne a l'aide de la technique du carottage
sismique consiste fig a mesurer des vitesses moyennes a des niveaux fixes a
l'avance par le sismicien.Pour chaque profondeur,on descend un geophone dans
le puits a la cote desire et on effectue un tir a l'explosif situe a une distance ne
depassant pas 0.5 kilometres de la sonde.Les trajets mesures sont corriges de
l'effet d'obliquite,puis normalise par rapport a un plan de reference.Les vitesses
moyennes depuis le plan de reference jusqu'a la profondeur du sismographe sont
calculees selon la relation:

Vm = HS / T

HS : profondeur plan de reference -sismographe ; T : temps vertical de propagation

de l'onde suivant le trajet HS :

VI-1-4-3 Vitesse d'Intervalle:

VI-1-4-3-1 Determination Indirecte:

Parmi les 5 types de Vitesses,la vitesse par tranche est la plus importante du point
de vue de l'interpretation sismostratigraphique C'est pratiquement la seule vitesse
qui est consideree comme parametre "diagnostic" de la variation de la lithologie,de
la porosite,de la presence de fracturation et de fluide dans les formations
geologiques.Elle est un parametre cle dans l'interpretation stratigraphique des
donnees sismiquesmen temps que parametre interne dans les analyses des facies
sismiques.La vitesse d'intervalle est determinee a partir des analyses de vitesses
au cours du processus de traitement a l'aide de la formule de "Dix".

V1-2 = [( V2RMS2T2 - V2RMS1T1 ) / (T2 -T1 )]1/2

86
 87

VRMS2 ; VRMS1 : Les vitesses quadratiques moyennes liees aux temps T2 et T1 .Cette
formule n'est reellement applicable que pour les couches paralleles et horizontales
et a condition d'assimiler les valeurs des vitesses de sommation deduites des
analyses de vitesses a des vitesses quadratiques moyennes.De grandes
precautions sont a prendre en presence d'un pendage des horizons,car la vitesse
VRMS deduite est generalement plus forte que les valeurs reelles.C'est pourquoi,on
multiplie habituellement les resultats par le cosinus de l'angle de pendage. Lorsque
les reflecteurs ne sont pas paralleles,les erreurs peuvent devenir plus importantes
a cause de differentes perturbations telles que les interferences,les diverses sortes
de bruits,les distorsions provoquees par des anomalies de vitesses de la zone
alteree.Pour y remedier,seule la modelisation peut apporter une precision
acceptable sur les vitesses par tranche.Dans tous les cas de faibles erreurs
provoquees pqr la mesure des courbes normales peuvent produire des erreurs
importantes dans les vitesses d'addition,provoquant ainsi un calcul errone des
vitesses d'intervalles.

VI -1- 4- 4 Autres moyens:

La vitesse par tranche peut etre calculee egalement a partir de la vitesse moyenne
deduite d'un carottage sismique selon la formule:

V1-2 (TK ,TN) =(VMnTn -VMkTk) / (Tn-Tk )

ou : VMn et VMk sont les vitesses moyennes verticales au toit et au mur de la

couche et Tn et Tk les temps correspondants.Il esiste d'autres procedes plus
rapides qui permettent de detecter les vitesses par tranches dont l'un d'entre eux
(procede GSI) consiste a representer les points maxima des analyses de vitesses,
dans un systeme d'axes du type:

V2.T = f(T).

VI - 1- 4-4-1 Interets des Vitesses d'Intervalle:

L'analyse et l'etude des Vitesses par tranche permettent dans des conditions
favorables de nous renseigner sur les variations lithologiques dues a la
porosite.Elles peuvent egalement nous donner de precieux renseignements sur

87
 88

l'effet de la fracturation et de la fissuration des roches,sur l'influence de la

cimentation et sa variation,ou de la mineralisation.Le champ d'application de la
vitesse par tranche,ne cesse de s'accroitre. Au stade actuel,elle est utilisee pour la
delimitation de la couche reservoir renfermant les hydrocarbures.Dans le plan
vertical,la vitesse par tranche peut egalement nous renseigner sur les variations
lithologiques a condition de proceder a la correction due au phenomene de
compaction. En Interpretation des donnees sismiques,la reconnaissance d'une
zone epaisse et sourde ou les analyses de vitesse s'averent efficaces, pour mettre
en evidence les zones sous-compactees,est un facteur favorable a la presence des
hydrocarbures. La vitesse par tranche permet,en outre,de determiner le type de
lithologie,en particulier dans l'identification des sables satures en Gaz (ceux-ci se
caracterisent par une faible valeur des vitesses d'intervalle par rapport aux sables
satures en huile) et enfin de calculer le rapport sable / argile ou calcaire / gres
etc..., Notons que la combinaison des vitesses par tranche avec les autres
parametres,comme par exemple,les amplitudes et les analyses spectrales,permet
de tracer des cartes de variation laterale de facies,d'identifier et de delimiter les
recifs etc...

VI -1-4-5 - La Vitesse Quadratique Moyenne

Si la vitesse d'intervalle et la vitesse moyenne dependent de la profondeur et du

temps de propagation,la vitesse quadratique moyenne depend en plus de la
distance "X" entre l'emetteur et le recepteur et affecte donc la corbure des
indicatrices . L'expression mathematique des trajets theoriques selon la loi de
l'optique est la suivante:

TX,n=T0,n +X2/V2n+C4X4 +C6X6 +......C2mX2m

T0,n :temps double vertical surface-miroir ; X:distance emetteur-recepteur.

Habituellement la serie est tronquee au terme en X2 ce qui fournit une precision
suffisante.La vitesse dans l'expression ci-dessus est appelee vitesse moyenne
quadratique.Elle est definie par:

V2RMS = SV2K.tK / StK

VRMS =(1/T0SV2K.tX )1/2.

88
 89

C'est cette vitesse qui rentre dans l'equation des indicatrices.

TX=[T0 +X2/V2RMS] 1/2

Elle est mesuree a partir des indicatrices des enregistrements sismiques,elle sert
essentiellement a ramener les trajets reflechis a des trajets verticaux.La
vitesseVRMS peut etre definie a partir de la vitesse instantanee comme la racine
carree de la moyenne temporelle du carre de la vitesse instantanee.

V2RMS = 1/T0 òV2i(t).dt

La determination des vitesses VRMS a partir des donnees du carottage continu de

vitesse permettent em appliquant cette formule,de deduire des valeurs de Dt avec
une bonne precision.

VI - 1-4-6 Vitesse de Sommation:

C'est une vitesse qui ne peut etre mise en equation en raison de son lien avec les
trajets sismiques.Elle se determine uniquement a partir des analyses de
vitesses,ce qui,en quelques sorte,optimise le mieux possible l'amplitude de la trace
sommee lors du processus de traitement des donnees sismiques acquises en
couverture multiple.
VI-2 - Principaux Facteurs Modifiant La Vitesse dans les Differentes Roches:

VI -2-1 Introduction:

La methode permettant de remonter d'une reflexion pointee sur un film,au miroir qui
lui a donne naissance,comporte toujours une part d'approximation et
d'hypothese,notamment sur la repartition des vitesses dans le sous sol.La
conversion des temps en distance,s'opere a l'aide de la loi de vitesse;cependant
on n'est jamais a l'abri d'une variation de vitesse dans le sens horizontal et
vertical,car on ne connait pas la vitesse en tout point.Le rapport entre la vitesse de
propagation des ondes sismiques et les coefficients (modules)d'elasticite du milieu
;les modules d'elasticite et la densite des roches varient en fonction de leurs

89
 90

compositions lithologique,de pressions et d'autres facteurs.La densite des roches

varie relativement dans de petites limites: de 1.5 - 1.7 a 3.0 - 3.1 gr/cm3 .D'ou on
peut conclure,que la variation de densite n'influe pas beaucoup sur la vitesse de
propagation des ondes elastiques.Mais le module de Young (E) varie dans de
larges limites;pour differentes roches,il differe de 100 fois.La vitesse de
propagation des ondes sismiques depend en premier lieu de la lithologie des
roches.Elle depend aussi de plusieurs facteurs;parmi ceux-ci on peut citer:la nature
de la roche,son age,la pression,la profondeur,la porosite, etc...

VI-2-2 - La lithologie:

La composition mineralogique des roches varie largement pour chaque type;mais

on peut indiquer les limites possibles des vitesses : Air (310 a 360 m/s) ; Sables
(100 a 1800 m/s) ; Argiles (1200 a 2700 m/s) ; Gres (1700 a 4000 m/s) ;
Calcaires,dolomite (2400 a 6000 m/s) ; Granite (4000 a 5700 m/s) ; Basalte (5000
a 6500 m/s).Il est a noter que pour les roches argileuses,les sables et pour
certaines roches sedimentaires il y a une relation lineaire entre la vitesse et la
densite:

V = 6.r-11

ou V est exprimee em Km/sec.

VI-2-3 La Nature de la Roche:

La porosite ainsi que la teneur en silice diminue la vitesse.La teneur en eau dans
les roches meubles (sables) l'augmente alors qu'elle diminue dans les roches
consolidees.

VI-2-4- Age de la Roche:

L'Age geologique de la roche augmente la probabilite d'action des pressions

(profondeur,forces orogeniques).La vitesse depend donc du passe de la roche.
Dans les memes conditions,les roches anciennes ont des chances d'etre plus
rapides,on note la relation suivante entre l'age absolue de la roche et la vitesse de
propagation des ondes sismiques:

90
 91

V = K(h.T)1/6

ou K: est le coefficient,qui depend de la lithologie de la roche ; h: est la profondeur

; T: etant l'age absolue.

VI-2-5 - Pression-Profondeur:

Ces deux parametres sont etroitement lies,puisque la pression augmente avec la

profondeur,mais l'effet de la pression est le plus important.La vitesse augmente
avec l'accroissement de la pression.

VI-2-6 - Anisotropie:

Les vitesses sont plus grandes suivant la direction des strates que dans la
direction perpendiculaire.Ce fait est important,car on voit que les vitesses
horizontales deduites des dromochroniques sont superieures aux vitesses
verticales deduites par exemple d'un carottage.

VI-2-7 La Porosité:

La Vitesse sismique est essentiellement soumise a l'influence du facteur

porosite.Toute augmentation de ce facteur amenuise la vitesse.Pour les formations
clastiques,la porosite depend de la pression differentielle d'une maniere
irreversible avec la profondeur d'enfouissement.Elle devient generalement plus
faible a mesure que diminue la separation des elements de la roche et
qu'augmente la cimentation.A des valeurs elevees de porosite sont generalement
associees des faibles valeurs de vitesses, et les faibles valeurs de porosite a des
hautes valeurs de vitesse.

VI-2-7-1 La Porosite Primaire

Les pores primaires se forment en meme temps que les roches.Ils changent leurs
tailles,leurs formes dans le temps grace aux differents processus tels que la
cimentation,la consolidation,le metamorphisme etc...,Ils trouvent leur origine dans:-
les pores structuraux entre les grains des formations

91
 92

clastiques(sables,gres,calcaires detritiques, etc...);-les espaces entre les surfaces

d'accumulation des depots des roches sedimentaires;-les pores de formations
enrichies par les squelettes organiques;-les pores intergranulaires et
intercristallins(pores entre les cristaux et les oolithes dans les calcaires et dolomies
etc...);- les pores formes par les bulles de gaz dans les cristaux des roches
magmatiques.
VI-2-7-2 La Porosite Secondaire

Les pores secondaires correspondent aux vides dont l'origine est liee a des
processus de diagenese.Ils peuvent prendre naissance a la suite des phenomenes
geologiques de lixiviation,de deformation des roches sous l'effet de la pression de
cristallisation primaire et secondaire,de deshydratation,de dolomitisation,d'erosion
etc...;La porosite est etroitement liee a la taille des pores qui peuvent etre classes
d'apres leurs dimensions en 4 groupes: - Les "megapores": Sont des pores de
diametre superieures a 10 mm.Ils sont typiques pour les calcaires et les gypses. -
Les pores "supercapillaires" (0.1 - 10 mm) sont rencontres dans les graviers,les
sables a gros grains et a grains moyens,les calcaires a dolomies caverneux etc...-
Les pores "capilaires" (2.10-4 a 0.1 mm) sont caracteristiques des sables et des
gres moyens et fins,des calcaires peu cimentes etc... - Les pores "subcapillaires"
(2.10-4 <) sont typiques pour les roches bien dispersees telles que les argiles ,les
roches argileuses,les calcaires aux cristaux fins a porosite primaire etc...
D'un point de vue petrolier,on considere comme constituant un tres bon
reservoir,les roches dont la porosite atteint 20% a 25% et plus,bons reservoirs
celle ou le pourcentage des vides est compris entre 15 % e t 20 % ,comme
reservoir moyens les roches ou la porosite varie entre 10 % et 15% et negligeable
pour les porosites inferieures a 5%.La plus forte porosite correspond aux roches
clastiques peu cimentees et la plus faible s'observe dans les roches detritiques
compactes,magmatiques et metamorphiques.La classification ci-dessus des
reservoirs ne constitue toutefois que des indications d'ordre general,car de
nombreuses roches peuvent presenter des pourcentages de vides importants sans
etre pour cela de bonne roches reservoirs ou inversement.(En Algerie,la porosite
des deux principaux gisements de petrole et de gaz est de : 15% a 20% dans les
gres triassiques du gisement de gaz de Hassi-R'mel; et 5% a 15% dans les gres
quartzitiques du cambrien du gisement de Hassi-Messaoud.) La porosite des
roches detritiques est etroitement liee a :

92
 93

- L'influence de la forme et de la dimension des grains

- L'influence du coefficient du triage et du degres de cimentation.
- L'influence de la consolidation des grains.
- L'influence de l'arrangement des grains, etc...

La porosite des roches corbonatees presente un interet tout particulier du point de

vue accumulation des hydrocarbures etant donne qu'elles occupent une place de
tout premier plan avec 50% des reserves mondiales.La porosite primaire des
roches cabonatees resulte principalement de sediments chimiques et
biochimiques.Elle peut avoir au debut une valeur considerable et depend du
degres de consolidation.Elle decroit avec le temps sous l'effet de la pression.Aussi
estime-t-on que les roches carbonatees ont souvent une porosite secondaire qui
prend naissance a la suite des processus posterieurs de diagenese tels que la
fissuration,la dissolution et la dolomitisation accompagnee d'une contraction du
volume des pores.On concoit que la porosite des roches carbonatees soit
influencee par de nombreux facteurs contemporains ou posterieurs a leur periode
de depot et que les reservoirs carbonates presentent par rapport aux reservoirs
detritiques terrigenes des caracteres originaux et specifiques.

VI -2-7-3 Determination de la Porosite a Partir des Diagraphies

De ce qui vient d'etre dit plus haut,il est evident que le facteur porosite joue le role
essentiel dans les variations des vitesses sismiques.(Wyllie et All 1951)
mentionnaient que l'influence de la pression sur la vitesse diminue avec la
profondeur d'enfouissement. Autrement dit,au-dela d'une certaine profondeur,les
couches geologiques se comportent comme des roches bien consolidees.La
pression n'a plus alors un grand role a jouer et la vitesse depend en grande partie
de la porosite.Les principales techniques de determination de la porosite restent du
domaine des proprietes geophysiques du puits. Ce sont les methodes
soniques,neutroniques,diagraphie gamma-gamma (densite) et les diagraphies
electriques qui sont largement employees dans la determination de la porosite des
formations geologiques.

VI -2-7-3-1 Diagraphie Sonique

93
 94

Parmi les methodes enumerees ci-dessus,la methode sonique semble etre la plus
interessante du point de vue de la determination de la porosite a partir des
donnees de vitesses dans les roches reservoirs a porosite intergranulaire et
intercristalline.C'est a Wyllie que revient le merite de formuler l'etude de la vitesse
en fonction de la porosite. Sa formule est bien justifiee dans les formations
consolidees et dures telles que les gres et les calcaires.Elle est:

f = (DTtotal - DTmatrice ) / ( DTfluide - DTmatrice ) .

Elle exprime la porosite a partir du temps de parcours de l'onde dans la roche

reservoir. Celle-ci est generalement deduite du log sonique.
En supposant que pour les milieux poreux le temps de parcours total "DTtotal" soit
la somme du temps de parcours dans la matrice et du temps de parcours dans le
fluide qui remplit les pores,on a :

DTtotal = f. DTfluide + DTmatrice .(1-f).

La formule precedente dite "relation temps-moyen" peut etre generalisee pour une
roche constituee par plusieurs composants de differentes porosites :f1; f2; f3 ;........
fn.
On aura:

DTtotal = f1. DT1 + f2DT2 +...........+ fn.DTn +

avec : f1+ f2+...........+ fn.=1

Pour le calcul des porosites dans les roches reservoirs granulaires a l'aide des
donnees de diagraphies sonique,on utilise le plus frequemment la relation
empirique appelee relation temps-moyen.:

Vp = Vf .Vma / (f.(Vma - Vf ) + Vf ) .

C'est une relation qui s'applique essentiellement a la determination des porosites a

partir des mesures de vitesses dans les roches reservoirs de type granulaires
(calcaires gres etc..).Elle suppose que lq duree du trajet par unite de longueur

94
 95

dans une roche poreuse contenant un fluide est la moyenne des durees des trajets
par unite de longueur dans la matrice et dans le fluide 1/Vma et 1/Vf . La figure
illustre la relation vitesse-porosite experimentee sur de roches reservoirs saturees
en eau en appliquant la formule ci-dessus.

VI -2-7-3-2 Diagraphies Neutroniques

(Gamma-Gamma ,pour un sable a gaz)

Selon la societe Schlumberger,il existe une relation entre la porosite lue sur les
diagraphies neutroniques dans un sable et la porosite reelle de ce sable.

freelle =[1/ fN+a.SGN +(1-SGN )]

ou :fN porosite lue sur un abaque etabli pour une saturation a 100% d'un liquide
protonique.freelle :porosite reelle ; SGN :saturation en gaz dans la zone exploree par
les neutrons ; a : rapport des concentrations en hydrogene de l'eau et du gaz dans
les conditions de reservoir.

VI -2-7-3-3 Diagraphies de Densite (Gamma-Gamma ,pour un sable a gaz)

Pour un sable a gaz,la porosite reelle est exprimee par la formule suivante:

fR = (de - dma ) / [(Sgd.(dg-1)+(1- dma )]

dg:densite du gaz ; Sgd :saturation en gaz dans la zone exploree par le gamma-
gamma.En outre,la diagraphie gamma,gamma permet d'evaluer la porosite totale
des formations,au'elle soit intercristalline , intergranulaire ou vacuolaire.

VI -2 -7-3-4 Diagraphie Electrique

Les methodes electriques (resistivite) permettent egalement d'evaluer la porosite si

l'on connait la resistivite Rw de l'eau contenue dans la roche et R0 .D'apres la
formule de D'Archie la porosite est liee au facteur de formation par:

F= R0 / Rw = a/fm

95
 96

ou :R0 :resistivite de la roche saturee par la saumure ;fm :coefficient de

cimentation. Si l'on connait Rx0 ,on a :

F=Rx0 /Rmf = a/fm .

Le coefficient "a" varie d'une couche a l'autre.Pour une roche compacte,elle est
egale a:

F = 0.81/f2

Pour le sable :

F = 1/f2.

En outre,la connaissance de la resistivite des roches clastiques la ou il n'existe pas

de donnees sur les vitesses,permet de deduire celles-ci a l'aide de la formule de
Faust:

Vms= 1948.(Z.R)1/6

ou : R resistivite ; Z:profondeur

VI -2-7-4 Porosite dans les Formations Consolidees

et non Consolidees

Comme nous le savons,le petrole est toujours lie aux roches sedimentaires.Les
sediments meubles se transforment progressivement en roches coherentes sous
l'effet de la diagenese.Ce phenomene reunit (d'apres GUILLOT) le tassement,le
cimentation,et la recristallisation.Des travaux sur les mesures de densite et de
porosite sur des sediments detritiques terrigenes essentiellement argilo-sableux
montrent que ceux-ci subissent au cours de leur enfouissement une compaction
qui se traduit par une augmentation reguliere de la densite et une diminution de la
porosite.Les travqux experimentaux de Y.LACAZE (Elf-Acquitaine) ont conduit
l'auteur a proposer une relation mathematique entre la porosite et la profondeur:

96
 97

f = a.ec.z

ou :z La profondeur. L'experimentation de cette formule sur de nombreux puits a

conduit a formule suivante:

f = 0.4477.e(-6.745.(10-5).zf)

Par analogie,le tassement des sediments sous l'effet de l'enfouissement se traduit

par une augmentation reguliere de la vitesse des ondes sismiques;on obtient:

Pour le sable: V = 1966.e(5.54(10-5).zf)

Pour les argiles V = 1431.e(6.627.(10-5).zf)

Contrairement aux formations non consolidees,les caracteristiques physiques des

formations consolidees sont moins sensibles a l'effet de l'enfouissement,et la
nature des fluides contenus dans les roches affecte beaucoup moins la vitesse
sismique.
VI - 2-7-5 Role de la Densite

La densite des roches est un parametre aussi important que la vitesse.Elle est une
fonction de la densite des divers mineraux qui composent la roche.Elle est mise en
evidence essentiellement par les methodes gravimetriques,par les diagraphies
gamma-gamma et dans une certaine mesure,par la sismique reflexion sensible a
l'impedance acoustique des roches (produit de la densite par la vitesse).La densite
des roches croit generalement avec la profondeur d'enfouissement en raison de
l'expulsion de l'eau des pores sous l'effet de l'augmentation de la compaction.La
relation qui lie la vitesse a la densite,la plus couramment utilisee est la formule de
Gardner et All - 1974:

d = 0.31 V 0.25 (m/s)

d = 0.23 V 0.25 (ft/s)

97
 98

Cette relation est representee sur la fig le sel et l'anhydrite et le charbon

n'obeissent pas a cette loi.Par analogie avec les relations precedentes liant la
porosite,la vitesse et la profondeur;on a pour des sables:

d = 2.65 - 7.14.10 3 V -1.21 (V en cm/s)

d = 2.65 - 0.7387.e(-67045.(10-5).zf)

Pour les roches renfermant des hydrocarbures,la densite en fonction de la porosite

est :

d=dma(1-f).dhf+Sw(dw-dh)

la compressibilite est:

K=K1.f+Kma.(1-f)

ou :dh densite des hydrocarbures

Sw saturation de l'eau ;
dw densite de l'eau ;
K1 et Kma. la compressibilite du liquide et de la matrice.

VI-3- Zone Alteree - WZ :

La couche superficielle de la coupe geologique est presentee par des roches

meubles pour la plupart des regionsmou on effectue des travaux sismiques.La
couche superficielle meuble peu consideree est caracterisee par des petites
valeurs des vitesses des ondes sismiques est nommee "zone alteree" ou "zone de
la WZ en Anglais".Dans cette zone la vitesse Vp.tr .peut etre petite,meme moindre
que la vitesse du son dans l'air.Les parametres de la zone alteree ( son epaisseur
et la vitesse de propagation des ondes sismiques dans cette zone) sont
etroitement lies a l'humidite des roches et a la position du niveau des eaux
souterraines (niveau hydrostatique),car l'accroissement de l'humidite des roches
mene a l'augmentation de la vitesse et cell-ci se rapproche de la valeur
normale,d'habitude,plus que 1500 a 1700 m/sec. C'est pourquoi pour la plupart des

98
 99

zones d'etude,la base de la zone alteree coincide avec le niveau hydrostatique.Par

exemple,pour une region de jijel on a presente la zone alteree fig ,dont la base
est en parfaite coincidence avec le niveau hydrostatique.Mais il y a des regions,ou
on ne peut etablir aucune liaison entre le niveau hydrostatique et la base de la
zone alteree.Comme par exemple,nous donnons la zone alteree d'une region de
Hassi-Messaoud. fig . Pour les terrains,ou les roches de meme composition
affleurent a la surface de la terre,on peut dire que la base de la zone alteree suit la
forme du relief.L'epaisseur de la zone WZ varie dans de larges limites: de 1 - 2m a
80 - 100m et plus;le plus souvent elle est de 8-15m. Dans nos regions
principalement dans le sahara Algerien,cette zone est surtout liee aux dunes de
sable sec,tres defavorable et son epaisseur est souvent superieur a 100m. La zone
alteree n'existe pas dans les regions de la cogelation eternelle.Les vitesses
peuvent varier dans les sens horizontal et vertical.Le caractere de variation de la
vitesse dans le sens vertical depend tout d'abord de l'epaisseur de la zone
alteree,c'est pourquoi ayant en vue l'accroissement du gradient de vitesse en
profondeur,on distingue la zone alteree d'une couche,de deux couches et parfois
de trois couches.Notons que la presence de la zone alteree mene a la refraction
forte des ondes sismiques revenant des grandes profondeurs et les rayons des
ondes sismiques se rapprochent du sens vertical.De plus la zone alteree se
caracterise par un grand coefficient d'amortissement.En passant par cette zone les
ondes sismiques changent sensiblement leurs compositions spectrales;leurs
composantes hautes frequences etant absorbees par la zone alteree. La zone
alteree constitue un obstacle,car elle est tres perturbantes en prospection sismique
reflexion haute resolution.Car cette zone decale dans le temps les trajets
sismiques les unes par rapport aux autres,en raison de sa faible vitesse et joue le
role d'un filtre coupe haut des frequences sismiques.D'ou il est indispensable,si
l'on veut pratiquer une sismique haute resolution,de disposer les points de Tir et
ces geophones sous la zone alteree.Chose pratiquement irrealisable;On effectue
alors des Corrections,dites "Corrections Statiques".

99
 100

CHAPITRE VII REPRESENTATION GRAPHIQUE DES

DIFFERENTES ONDES

VII-1 - Hodochrones de Differentes Ondes :

Les hodochrones sont les donnees de depart essentielles lors du traitement et de

l'interpretation des resultats de la prospection sismique.On distingue les
indicatrices lineaires et superficielles (tridimentionnelles).Dans la pratique de la
sismique les indicatrices superficielles sont presentees par des cartes en
isochrones.il est a noter que dans la theorie des indicatrices on suppose que les
observations sont toujours realisees sur le plan horizontal.En sismique on presente
les indicatrices lineaires sous la forme des graphes des temps d'arrivees des
ondes en fonction des coordonnees des points d'observations.On distingue encore
les indicatrices obtenues le long des profils longitudinaux et non-longitudinaux.Lors
des observations le long des profils longitudinaux les points d'excitation et les
sismographes sont disposes sur la meme ligne et comme origine du systeme de
coordonnees,on prend le point d'excitation ou son epicentre fig .Pour le profil non-
longitudinal comme origine du systeme de coordonnees,on prend la projection du
point d'excitation sur le profil fig .Lors des observations le long des profils
circulaires on peut prendre comme origine un point quelconque fixe fig .A
present,considerons les indicatrices des ondes sismiques pour un mlieu a deux
couches.Les vitesses de propagation des ondes sismiques dans ces milieux
differents.Le point d'excitation est au point O.

VII-2 Champ de Temps et Indicatrices des Ondes Directes :

100
 101

Examinons une onde directe issue du point de Tir O et se propageant dans le

premier milieu avec la vitesse v.Placons l'origine du systeme de coordonnees
(x,y,z) au point de tir O. fig .L'axe Z est dirige vers le bas (verticalement).Prenons
egalement un parametre "r"(rayon),qui determine la distance entre un point
quelconque du milieu et la source O.Dans ce cas nous avons:

¶t/¶x=dt/dr.(x/r)

[cos(x,r)=x/r];¶t/¶y=dt/dr.(y/r)
[cos(y,r)=y/r];¶t/¶z=dt/dr.(z/r)
[cos(z,r)=z/r]

En les portant dans l'equation de champ de temps nous trouvons que:

( dt/dr)2.(x/r)2 + (dt/dr)2.(y/r)2 + ( dt/dr )2.(z/r)2 = 1/V2

( dt/dr)2[(X2 + Y2 + Z2 ) / r2] = 1/V2

ou : r2 = X2 + Y2 + Z2 ; dt/dr = 1/V .

En integrant l'equation,nous trouvons l'expression suivante pour la fonction de

champ de temps:

t=r/V+C.

La constante "C" doit etre determinee d'apres les conditions initiales.En supposant
que l'onde sismique est engendree par la source "O" a l'instant t=0;et en
considerant le champ de temps pour le front d'onde on obtient C=0.Alors la
fonction de temps est egale a:

t = r /V = 1/V = Ö(X2 + Y2 + Z2 )

Il en resulte que les isochrones de l'onde directe sont l'ensemble des spheres
concentriques avec les rayons :

101
 102

ri = vi.ti

ayant le centre O.Il est tout a fait evident que les rayons sismiques coincident avec
ceux des spheres concentriques.Considerons les cas particuliers suivants:

VII-2-1 - Lorsque la Source se trouve a la Surface de la Terre:

Lorsque la source se trouve a la surface de la terre nous avons pour z=0 :

t=1/V(ÖX2 + Y2)

ou bien

X2/V2 + Y2 /V2 - t2/ 1 = 0

L'indicatrice superficielle de l'onde directe est un cone dont le sommet se trouve a

l'origine du systeme de coordonnees fig .

VII-2-2 - Lorsque l'Indicatrice est Lineaire:

L'indicatrice lineaire peut etre trouvee en supposant que : (y = 0 , z = 0 ):

t=±X/V

avec : t = ± K.X ou K = 1/V.

Elle consiste en deux segments des droites issues de l'origine du systeme de

coordonnees fig .La pente de l'indicatrice depend de la vitesse de l'onde
directe.Lorsque la source O se trouve au voisinage de la surface terrestre,les
ondes superficielles et aeriennes se propagent dans tous les sens du point de
tir.Leurs Indicatrices coincident avec celles de l'onde directe.Mais pour les
indicatrices de ces ondes,il est necessaire de prendre leur vitesse correspondante.

VII-2-3 - Indicatrice des Ondes Reflechies.

VII-2-3-1 Indicatrice Superficielles pour un Reflecteur Plan :

102
 103

Ce cas est d'une grande importance dans la pratique de la sismique car : * - Les
discontinuites geologiques reelles sont planes dans les limites des profils : * - La
solution mathematique du probleme est tres simple.Admettons que dans un
milieu,nous avons un reflecteur de pendage j.Ce reflecteur divise le milieu
considere en deux couches de vitesse V et V1,c'est-a-dire la vitesse de
propagation des ondes sismiques entre le reflecteur et la surface de la terre est
egale a V.Dirigeons l'axe "x" du systeme de coordonnees dans le sens de l'amont
du reflecteur et l'axe "y" perpendiculairement a l'axe "x". fig .L'axe "z" est dirige
verticalement (vers le bas).En construisant le point de Tir imaginaire
O'(x0,y0,z0),on peut trouver que: t(OMS) = t(O'S); et les points O',M et S se
trouvent sur la meme droite.Les coordonnees de la source Imaginaire (fictive) sont:

x0 = 2h.sin jy0=0

z0 = 2h.cosj.

ou :h est la profondeur du reflecteur suivant la normale passant par la source;j:

l'angle du pendage reel. Donc,d'apres la figure ,on peut constater que tous les
rayons des ondes reflechies sont issus de la source imaginaire et coincident avec
les rayons de la sphere ayant pour centre le point O'.Cela nous permet d'ecrire
l'equation de champ de temps de l'onde reflechie sous la forme suivante:

( X - X0 )2 + ( Y - Y0 )2 + ( Z - Z0 )2 = V2. t2

Pour trouver l'equation de l'indicatrice superficielle de l'onde reflechie il est

necessaire de resoudre cette equation avec celle du plan d'observation.Pour ce
cas :z=0 ;L'equation devient.
En portant les expressions:

x0=2h.sin jy0=0 ; z0= 2h.cosj

dans l'equation:

( X - X0 )2 + ( Y - Y0 )2 + ( Z - Z0 )2 = V2. t2

103
 104

on obtient :

(X-2h.sinj)2 + Y2 =V2. t2 - 4h2.cos2 j

Dans le systeme de coordonnees (x,y) et avec le parametre "t",c'est une equation

de l'ensemble des cercles concentriques ayant le centre deplace vers l'amont de
valeur :

x0=2h.sinj

suivant l'axe "X".Au point (x0,0) le temps sera minimal :

tmin = 2h.cosj/V

ri = ÖV2. t2i - 4h2.cos2j.

Dans le systeme de coordonnees (x,y,t) l'equation de l'indicatrice a la forme

suivante (z=0) :

(X - X0)2 + Y2 - V2. t2 = -Z02

( X - X0 )2/ Z02+Y2/ Z02- t2/ (Z0/V)2 = -1

(X-2hsinj)2/4h2.cos2j.+ Y2/4h2.cos2j - t2 / (2h.cosj/V)2 = - 1

C'est une equation d'une hyperboloide dont la partie negative n'a pas le caractere
physique (temps "t" etant toujours positif). fig .

VII-2-3-2 - Indicatrice Longitudinale (le long du profil longitudinal):

A present essayons de trouver l'equation de l'indicatrice longitudinale pour un profil

arbitrairement oriente par rapport a la direction de l'amont-aval du reflecteur
.L'equation de la ligne d'observation (du profil) peut etre ecrite par l'expression
suivante :

104
 105

y = X.tgY

ou Y: est l'angle entre l'axe "x" (qui passe par la direction de l'amant-aval) et le
profil sismique.En portant l'expression: y = X.tgY ; dans l'equation :

( X - 2h.sinj ) 2 + Y2 = V2. t2 - 4h2.cos2 j

On trouve:

X2.sec2Y - 4hx.sinj+4h2 =V2. t2.

Si l'on prend comme variable la distance "l" entre la source et les sismographes
(comme c'est l'usage en sismique);on :

l = X/cosY = X.secY

V2. t2 = l2 - 4h.l.sinj.cosY + 4h2; V2. t2 = l2 - 4h.l.sinja+ 4h2

ou :

sinja= sinj.cosY

L'angle ja: est nomme "angle apparent" du reflecteur dans le plan normal passant
par le profil.L'equation:

V2. t2 = l2 - 4h.l.sinja+ 4h2

est une equation de l'hyperbole

( l - l0 )2/Z20 - t2/ t2min = -1

ou :

Z0 = 2h.cosja ; l0 =2h.sinja ; tmin = 2h.cosja / V .

105
 106

Si le profil sismique a la meme direction que celle de l'amont-aval du reflecteur ( Y

= 0 ; ou 180° ) ; on a : l = X ; ja = j ;

V2. t2 = X2±4h.X.sinj+ 4h2

t = (ÖX2±4h.X.sinj+ 4h2) / V

ou : les signes "±" correspondent respectivement a l'amont (-) et a l'aval (+) du

reflecteur.Donc,l'indicatrice longitudinale d'une onde reflechie a la forme d'une
hyperbole,dont les parametres sont en fonction de l'angle apparent ja , de la
profondeur "h" et de la vitesse "V".Le minimum de l'hyperbole est toujours deplace
dans la direction de l'amont a la valeur 2h.sinja . Il convient de dire que le
deplacement le plus grand sera observe sur le profil ayant la direction de l'amont -
aval du reflecteur et le plus petit (egal a 0) sera observe sur le profil,qui passe
perpendiculairement a cette direction.La plus petite valeur du temps minimal (tmin)
de l'hyperbole sera observee dans la direction en amont du reflecteur : (Y=0) et
inversement,la plus grande valeur de ce temps sera observee dans la direction
perpendiculaire (Y =90° ).On peut analyser la forme de l'indicatrice longitudinale
en fonction de la variation de la vitesse,de la profondeur et de l'angle du
pendage.Si la vitesse de propagation des ondes sismiques augmente,la courbure
de l'indicatrice diminue et le minimum se deplace vers le bas fig .Plus l'angle de
pendage est grand plus le minimum se deplace en bas et vers l'amont du miroir fig
.La courbure de l'indicatrice diminue et le minimum se deplace en haut avec
l'augmentation de la profondeur;l'equation des asymptotes de l'hyperbole peut etre
trouvee d'apres la formule indiquee auparavant et qui s'ecrit:

( l - l0 )2/Z20 - t2/ t2min = -1

En egalisant a zero le terme de droite,on trouve l'equation des asymptotes de

l'hyperbole longitudinale:

V2. t2 = (l - 2h.sinja )2 ; t = 1/V(l - 2h.sinja )

Les asymptotes de l'hyperbole ont comme origine la projection du point minimum

de l'indicatrice sur l'axe "x ou l"; et leurs pentes dependent de la vitesse vraie du

106
 107

premier milieu.Par consequent elles sont paralleles a l'indicatrice des ondes

directes.

VII-2-3-3 - Indicatrice Non-Longitudinale (le long du profil non-longitudinale):

Soit un point d'ebranlement "O" a la distance "R" du profil sismique.L'angle entre

ce profil et l'amont-aval du reflecteur est egal a Y fig et fig .En systeme de
coordonnee xoy l'equation du profil s'ecrit de la forme:

y=x.tgY + R.secY

En portant cette expression dans l'expression prcedemment ecrite et qui est:

(X-2h.sinj)2 + Y2 =V2. t2 - 4h2.cos2 j

On obtient :

(X-2h.sinj)2 +(R.secY +XtgY)2 =V2. t2 - 4h2.cos2 j

Mais en realite,on sait,que comme origine du systeme de coordonnees,on prend la

projection du point de tir sur le profil "E".Pour trouver les relations entre les
nouvelles coordonnees et les anciennes,projections le segment ES=1 sur
l'axe;alors:

l.cosY = X+R.sinY Þ X = l.cosY- R.sinY

En portant cette expression dans :

(X-2h.sinj)2 +(R.secY +XtgY)2 =V2. t2 - 4h2.cos2 j

on obtient :

(l.cosY- R.sinY -2h.sinj)2 +[(R.secY +(2.cosY-R.sinY)tgY]2 =V2. t2 - 4h2.cos2 j .

Apres transformation on arrive a :

107
 108

( 1 - 2h.sinj.cosY )2 + q2 =V2. t2

ou : q = Ö(R+2h.sinj.sinY)2 + 4h2.cos2 j
l0 = 2h.sinj.cosY ; (l - l0)2 - V2. t2 =q2
(l - l0)2 /q2 - t2/(q/V)2= -1 .

Comme on voit,c'est l'equation de l'hyperbole.Dans le cas particulier si R=0,nous

obtenons l'equation de l'indicatrice pour le profil longitudinal.Donc,les indicatrices
lineaires des ondes reflechies seront toujours des hyperboles.La courbure de ces
indicatrices diminue avec l'accroissement de "R",(c'est-a-dire,de la distance entre
le profil et la source).Si le profil non-longitudinal passe dans la direction de l'amont-
aval du reflecteur,on a:

l0 = 2h.sinj; q = ÖR2 + 4h2.cos2 j .

Cela veut dire que "q" est l'hypothenus du triangle rectangle avec les cotes :
Z0 = 2h.cosj (la profondeur du point de tir imaginaire) et R.

VII-2-3-4 - Indicatrices des Ondes Reflechies: ( cas du miroir curviligne )

Chaque fois qu'on etudie la forme de l'indicatrice des ondes reflechies,on

supposait que le miroir considere etait plan.Il convient de dire,que quand la surface
reflechissante (miroir) est curviligne,l'indicatrice de l'onde reflechie a une forme
complexe,c'est-a-dire la forme de l'indicatrice peut differer notablement de la forme
hyperbolique.Examinons quelques exemples theoriques qui le confirment.Par
exemple,si la source se trouve dans la focale (distance) du miroir
parabolique,l'indicatrice a la forme d'une droite. fig Dans le cas ou le point
d'excitation est au centre d'un cercle,l'indicatrice se transforme en un point ;fig .
Un noeud peut se former lorsque la courbure du miroir varie dans les limites du
profil ,fig .Enfin la figure nous montre la dependance de la forme de l'indicatrice
en fonction de la courbure du miroir.

VII-2-3-5 - Indicatrice de Ondes Diffractees:

108
 109

On sait,que l'onde diffractee se propage dans tous les sens du point de

diffraction.Les fronts d'une onde diffractee peuvent etre construits d'apres le
principe de Huyguens.L'onde diffractee a ete consideree comme une onde
parasite.Mais ces dernieres annees on a commence a l'utiliser pour etudier des
elements structuraux,tels que les failles,ruptures de couches,dome de sel etc...
Soit une onde directe (incidente) qui se propage dans le milieu homogene.Lorsque
cette onde atteint un obstacle (une faille),au point "D" elle diffracte fig .Le temps
d'arrivee de l'onde diffractee au point d'observation:S(x,0) est:

t = OD/V + DS/V

ou bien

t=tdirect+ tdiffractee

ou : tdirect: est le temps mis par l'onde pour atteindre le point de diffraction "D" et il
est egal a :

tdirect= (ÖX20 + Z20) /V

OD = (ÖOD21 + DD21)

ou OD1= X0 Þ abscisse du point de diffraction; DD1 = Z0 Þprofondeur du point

de diffraction.De la meme maniere on peut trouver que:

tdiffractee=DS/V;DS=Ö(DD1)2+(D1S)2 DD1 =Z0

D1S =X - X0 ;DS=ÖZ02+(X-X0)2

tdiffractee= Ö(X-X0)2+Z02 /V

donc,l'equation de l'indicatrice d'une onde diffractee sera :

t = 1/V [ Ö X02+Z02 + Ö ( X - X0 )2 + Z02] .

109
 110

C'est une hyperbole symetrique par rapport a l'axe vertical passant par le point de
diffraction;c'est-a-dire;le minimum de l'indicatrice se trouve toujours au-dessus du
point de diffraction et le temps d'arrivee au point considere ( D1 ) est :

tD1= tdirecte +Z0/V = 1/V( Ö (X02+Z02 ) +Z0).

La courbure de cette indicatrice depend aussi de la vitesse et de la profondeur du

point de diffraction.

VII-2-3-6 - Indicatrice des Ondes Multiples:

Les ondes multiples sont des ondes qui se reflechissent ou se refractent plusieurs
fois sur le meme interface fig .Lorsqu'une onde ne se reflechit qu'une seule fois
sur un miroir,on dit qu'il s'agit d'une reflexion primaire ou simple.En revanche,si
cette onde primaire,apres s'etre reflechie sous la surface ou l'on a pratique
l'ebranlement,se reflechit de nouveau une ou plusieurs fois sur ce meme miroir,on
parle alors de reflexions multiples fig .Les trajets des reflexions multiples peuvent
etre tres complexes.C'est ainsi qu'en plus des multiples a tajets complets fig , on
rencontre des reflexions dont les multiples ne se sont produits qu'entre deux
miroirs situes sous la surface fig .On distingue encore deux types de reflexions
multiples:

VII-2-3-6-1 - Multiples de surface:

Il s'agit alors des reflexions entre un nouveau reflecteur profond et la surface du sol
fig

VII-2-3-6-2 -Multiples internes:

Dans ce cas les deux niveaux reflechissants sont profonds fig .Considerons la
deduction de l'indicatrice d'une onde reflechie-multiple.Soit "R" un miroir pente
(incline) sous l'angle j;qui donne naissance aux reflexions multiples et soit "h" la
profondeur du miroir suivant la normale fig .Tracons le trajet d'une onde multiple
"OABEC".Si le miroir fictif R' de profondeur h1 avec un point fictif de reflexion B'. Il
est evident,que AB=AB';BE=B'E et l'angle de pendage du miroir R' est egal a 2j. Il

110
 111

en resulte,que nous pouvons trouver l'equation de l'indicatrice des reflexions

multiples en utilisant l'equation des reflexions simples :

t=1/Vmul. Ö4h12 ± 4h12.X.sin2j+ X2.

ou Vmul : est la vitesse de propagation des ondes sismiques dans le milieu se

trouvant entre la surface terrestre et le miroir "R".Le rapport des profondeurs:

h / h1.= sinj/sin2j

ou h / sinj = h1 / sin2j
h1 =h.sin2j/ sinj

En portant cette derniere expression dans :

t=1/Vmul Ö4h12 ± 4h12.X.sin2j+ X2

On obtient l'equation de l'indicatrice d'une onde reflechie multiple 2 fois:

t=1/Vmul Ö4h2.sin2.2j/ sin2j + 4hX.sin2.2j /sinj+ X2

t=Öt02+4h.X/V2mulsin22j/sinj+ X2 / V2mul
ou :
tt0= 2h/Vmul.(sin2j/sinj).

Dans le cas general,l'equation de l'indicatrice d'une onde reflechie multiple (m fois)

s'ecrit de la maniere suivante:

t=1/Vmul Ö4h2m.±4hmsin mj + X2

ou :
hm = h.sin mj /sinj
m: est le nombre des multiples.Pour un miroir horizontal ( j = 0 )

On obtient :

111
 112

t=1/Vmul Ö4h2.m2.+X2 = Öt20 + X2 / V2mul .

Si l'on compare l'indicatrice d'une onde reflechie primaire et celle d'une onde
multiple pour le meme temps t0 on remarque,que la courbure de l'indicatrice d'une
onde multiple est plus grande que celle de l'indicatrice de l'onde reflechie.La
difference entre les deux indicatrices augmente avec l'accroissement du rapport
des vitesses de ces ondes.Dans les milieux geologiques reels la difference entre
les vitesses de ces ondes n'est pas importante et c'est pourquoi il est pratiquement
impossible de separer les reflexions multiples et simples d'apres les donnees de la
sismique par la couverture simple.

VII-2-3-7 - Dromochronique des Ondes Refractees:

En sismique refraction,on etudie les ondes refractees qui sont formees sur la
surface de separation de deux milieu dite "Marqueur" dont les vitesses sont en
relation de la maniere suivante:V1>>V;ou V1 est la vitesse du 2eme milieu. Soit OX la
ligne d'observation fig ;O le point de tir;j l'angle de pendage du marqueur;h la
profondeur du marqueur suivant la normale;V la vitesse du premier milieu.On sait
que les ondes refractees ne se forment que dans le cas suivant

sin il = V/V1

il : etant l'angle limite. Donc,le premier rayon d'une onde refractee cree par la
premiere discontinuite n'arrive qu'a la ligne d'observation a partir (apres) des points
M2 et M1. La distance M1M2 s'appelle "zone morte",ou bien la zone d'absence des
ondes refractees.Pour trouver l'equation de dromochronique,designons le point de
Tir imaginaire par O' et tracons les lignes O'B1 et O'B2. Il est evident
que:OA=O'A;OB2=O'B2. Il n'est pas difficile de trouver le temps d'arrivee de l'onde
refractee au point M2 :

tM2 = ( OB2 + B2M2 ) / V = ( O'B2 + B2M2 ) / V = O'M2 / V

112
 113

de D(triangle) O'O2M2:O'M2 = O'O2/cos (i+j)

de D(triangle) O'O2O: O'O2»2h cosj
d'ou:

O'M2 = 2h cosj/cos (i+j)

et:
tM2 =[2h cosj/cos (i+j) ] /V

XM2 = O2B2 - OO2

O2M2 = 2h cosj/tg (i+j)

OO2 =2h sinj

XM2 = 2h cosj.tg (i+j) - 2h sinj

apres transformation on obtient:

XM2 = 2h sin i / cos (i+j) .

L'equation d'une droite passant par le point F2(tM2 , XM2) ,peut etre ecrite de la
maniere suivante:

(t - tM2 ) / (X - XM2) = Dt / Dx = 1 / Va =sin (i+j) / V ;

En prenant les expressions tM2 , XM2 de :

tM2 =[2h cosj/cos (i+j) ] /V et XM2 = 2h cosj.tg (i+j) - 2h sinj

on trouve:

[t - 2h cosj./Vcos (i+j)] / [X - 2h sin i / cos (i+j)] = sin (i+j) / V

t = X. sin (i+j) / V + 2h / Vcos (i+j).[cosj - sin i .sin (i+j)

= X.sin (i+j) / V + 2h.cosi / V .

113
 114

Apres avoir fait les memes operations nous pouvons obtenir l'equation de la
deuxieme branche de la dromochronique :

t = X.sin (i - j) / V + 2h.cosi / V

Donc , l'equation de la dromochronique de l'onde refractee est :

t = X.sin (i ± j) / V + 2h.cosi / V

La dromochronique contient deux lignes droites appelees "branches" si l'on

prolonge les branches de la dromochroniques,elles se croisent ou se rapprochent
presque sur l'axe "t" en un meme point (I):

I = 2h.cosi / V.

Le temps "I" est appele "INTERCEPT". On peut dire c'est une solution particuliere de
l'equation de dromochronique pour X=0.Dans le cas d'un marqueur horizontal
:j=0,on a :

t=Xsini/V+2h.cosi /V

ou bien :

t = X / V1 + I ; (V1 = V / sin i ).

En analysant les equations de t= f(x) on peut conclure que: - avec l'augmentatiom

de la profondeur du marqueur et avec la diminution du contraste de vitesses V et
V1 (de deux milieux) la zone "morte" devient plus grande. - La difference entre
OM1 et OM2 croit avec l'augmentation de l'angle de pendage du marqueur.

VII-2-3-7-1 -Dromochronique:(cas des marqueurs curvilignes):

Si le marqueur est curviligne,on obtient le changement de la forme des

dromochroniques des ondes refractees et celles ne seront plus lineaires.Lorsque

114
 115

les vitesses de propagation des ondes dans deux milieux (V et V1) sont
constantes,l'angle limite est toujours constant,car:

sin il = V/V1

Va =V/sin a = V / sin (il +j) ; a = (il ±j)

si j=0 ; a = il ; a = const Þ Va =const

si j¹0 alors : a = (il ±j) , a ¹ const Þ Va ¹ const .

Il en resulte,que si le pendage du marqueur varie le long du profil,nous avons le

changement de l'angle d'incidence de l'onde refractee,puisque
a = (il ±j)

la vitesse apparente sera variable:

Va = V / sin (il ±j).

Les differents exemples des dromochroniques sont representes sur les figures

VII-3 Hodochrones des Ondes Sismiques dans le Milieu a plusieurs Couches:

Les milieux geologiques reels contiennent toujours de nombreuses strates ayant

des proprietes physiques differentes.En sismique on les considere comme des
milieux a deux ou plusieurs couches.Mais dans la plupart des cas on considere les
milieux stratifies comme des milieux gradients,c'est-a-dire,dont les proprietes
physiques varient continuellement.Lors de l'analyse des champs de temps et des
hodochrones des differentes ondes sismiques,on distingue les milieux
horizontalement-stratifies,si les pendages des miroirs ne depassent pas 10° a 15°,
et les milieux verticalement-stratifies si le pendage depasse 10° a 15°. Ici nous
nous bornons a analyser des cas simples permettant d'etablir les parametres
(proprietes) cinematiques des ondes sismiques le plus utilisees en pratique.

VII-3-1 - Indicatrice des Ondes Reflechies:

115
 116

Soit une coupe geologique a plusieurs couches.Chaque couche se caracterise par

sa propre epaisseur h1,h2,.....hn et par sa vitesse V1,V2,....Vn fig En prenant en
consideration le phenomene de refraction sur les interfaces intermediaires,tracons
un rayon sismique de l'onde incidente jusqu'a la n eme couche et egalement le
rayone de l'onde reflechie.Les angles de refractions sur les interfaces peuvent etre
trouves d'apres la loi de Snellius :

sin a / V1 = sin a2 / V2 = sin an-1 /Vn-1

Il n'est pas difficile de trouver le trajet de l'onde dans chaque couche et en sachant
la vitesse,on peut determiner le temps de propagation de l'onde reflechie;qui est
egal a:

t = 2.Shl / Vl.cosa2

X = 2..Shl tg al

ou : X est l'abscisse du point de reception.Si l'on decompose les expressions ci

dessus en series de Taylor,on obtient:

t=2.Shl /Vl.+X2/4(Shl/Vl.)-1-X4/64(Shl .V3l.).(Shl.Vl.)-4+..

Si l'on analyse un milieu a trois couches (de deux interfaces) avec le rapport des
vitesses V3>V2>>V1 fig .On aura un point d'intersection de ces deux
indicatrices.La distance entre la source et ce point depend de l'epaisseur h1 et h2 et
les vitesses dans chaque couche V1;V2;V3 respectivement,c'est-a-dire ,dans la
coupe il existe deux marqueurs : R1;R2 ; et V3>V2 V2>V1 ; fig .L'equation de
dromochronique pour un milieu homogene avec un marqueur est la suivante:

t=X/V2+2hl.cos i1,2 / V1

ou l'angle :

i1,2=arcsin V2/V1

116
 117

Pour trouver la dromochronique pour le deuxieme marqueur,on trace le rayon

sismique de facon que l'angle d'incidence soit egal a l'angle limite i2,3.Dans ce cas
l'angle sur la surface intermediaire peut etre trouve d'apres la loi de Snellius-
Descartes :

sin i1,3 / V1 = sin a2,3 /V2 = 1 /V3

ou : l'angle a2,3 est egal a l'angle limite i2,3 sur le deuxieme interface.La vitesse
apparente demeure toujours constante.car l'angle i1,3 ne varie pas (i1,3 etant l'angle
d'incidence sur le premier interface) et elle est egale a :

Va= V1 / sin i1,3 = V3

On peut aussi trouver Va , autrement:

Va =V3=DX/Dt

Il en resulte que la dromochronique est une droite dont le coefficient angulaire est
determine par V3. D'apres la fig ,on peut trouver les coordonnees de point initial
de cette droite:

Xm=2hl.tg i1,3 + 2h2.tg i2,3

tm= 2hl / V1.cos i1,3 + 2h2/V2.cos i2,3

Connaissant le coefficient angulaire,nous pouvons ecrire l'equation de la droite

passant par le point:Xm; tm:

DX/Dt = (X - Vm) / (t - tm).

En integrant les expressions:

Va=V3=DX/Dt

et: X m=2hl.tg i1,3 + 2h2.tg i2,3

117
 118

tm=(2hl /V1)cos i1,3 + (2h2/V2)cos i2,3

Dans la derniere expression qui est:

DX/Dt = (X - Vm) / (t - tm)

On obtient l'equation de la dromochronique des ondes refractees pour le marqueur

R2:

t=X/V3+ (2hl / V1).cos i 1,3 + (2h2/V2)cos i 2,3

De la meme maniere on peut formuler l'equation de la dromochronique pour le

marqueur Rm , quand la coupe se compose de "m" couches avec h1,h2,h3,..hm et
V1,V2,..Vm ;

t=X/Vm+ 2S(hk / Vk).cos i k,m ;m=1,2..n.

On peut avoir une dromochronique pour chaque interface,si la condition suivante

est satisfaite:Vk+1 > Vk ; Dans ce cas l'intercept sera egal a:

I=2S(hk / Vk).cos i k,m

VII-3-2 - Rapport entre les Hodochrones des Ondes Differentes:

VII-3-2-1 - Rapport entre les Indicatrices des Ondes Directes

et Reflechies:

Comme il a ete dit plus haut,pour un milieu homogene et isotrope les indicatrices
des ondes directes sont paralleles aux asymptotes des indicatrices des ondes
reflechies.Mais dans les milieux reels la vitesse de l'onde directe (V0) differe de la
vitesse moyenne (Vm) des ondes sismiques jusqu'au miroir et c'est pourquoi il est
possible que les indicatrices de ces deux types des ondes se croisent;dans le cas
d'un miroir parallele a la surface de la terre leur intersection aura lieu a la distance
"X" , satisfaisant la condition suivante :

118
 119

X21 /V20 = (X21 / V2m)+4h2 /V2m

d'ou:

X1 = 2h / Ö(X2m /V20 ) - 1.

Par exemple,si h=500 m;Vm = 2500 m/s;V0 = 1500 m/s ,les indicatrices se croisent
a la distance » 750 m du point d'excitation.Tres frequement,a ces
distances,l'intensite des ondes directes diminue notablement.

VII-3-2-2 - Rapport entre les Hodochrones des Ondes Reflechies

et Refractees:

En cas d'un interface plan (h,j) les hodochrones sont determinees par les
expressions ci-dessous: pour une onde refractee:

V1.ta =2h.cos i1,2 +X.sin (i - j)

pour une onde reflechie :

V21t2r = (X2 - 4h.sinj +4h2

ou : sin i = V1 /V2 ; ( V1 ;V2 : les vitesses de propagation des ondes sismiques

dans deux milieux).Comme d'habitude dans la direction de l'amont on prend le
signe de l'angle j positif. si x=0 :

ta0 =(2h/V1)cos i = I

tr0 = 2h / V1 .

En d'autres termes l'intercept de la dromochronique est toujours inferieur au temps

de l'indicatrice de l'onde reflechie au-dessus de point de tir.La difference entre ces
deux temps augmente avec la diminution de la difference entre V1 et V2 . Les deux
indicatrices ont un point commun coincident avec le debut de l'onde conique sur le

119
 120

profil.L'abscisse de ce point peut etre determine d'apres l'expression suivante fig

.

Xra = 2h. sin i / cos(i - j) .

On peut observer l'interference des ondes reflechies et refractees,si le marqueur et

le miroir ne coincident pas.

CHAPITRE VIII FILTRAGE SPATIAL

VIII - 1 Introduction:

Le filtrage spatial ou filtrage en nombre d'onde est une operation qui consiste a
disposer sur le terrain des geophones ou de sismographes d'une certaine facon

120
 121

afin de pouvoir selectionner les signaux utils,et par consequent eviter

l'enregistrement des autres signaux supposes indesirable,cette operation se
realise en choisissant correctement les caracteristiques reliants la distances entre
geophones,le nombre de geophones et la vitesse apparentes des ondes de
propagation (signaux utils et bruits).

VIII -2 Qu'appelle-t-on Vitesse Apparente ( Va ).

C'est la vitesse avec laquelle l'onde (signal ou bruit organise),parait se propager le

long du profil.On assimile,dans la mesure du possiblemcette onde a des arches de
sinusoides se propageant a vitesse apparente constante.Frequence (F) et Nombre
d'Ondes (K) sont definis a partir d'une onde plane sinusoidale,du temps (t) pour
une distance (x) donnee et sinusoidale de (x) pour un (t) donne.La figure ,montre
ce qui se passe pour un ebranlement,emis en "A",au temps t=0,et recueillie par les
traces "S et B",au bout des temps " tS et tB " respectivement.OA=T : est la periode
en fonction du temps;f=1/T frequence.OB=la est la periode en fonction de la
distance ou longueur d'onde apparente ;

ka=1/la

Nombre d'onde apparent.L'onde sinusoidale parait se propager le long de Ox a la

vitesse Va = distance parcourue / temps de parcours = OB / OA = la / T Nous ne
pouvons pas ecrire Va directement en fonction de la pente de la droite OM, dans le
diagramme ( Ot , Ox ),a cause de la difference des parametres utilises sur chaque
axe.Cependant dans un systeme ou les vecteurs unitaires "Ot et Ox"restent
identiques a eux meme,on remarque que,plus est forte l'inclinaison de OM sur
Ox,plus la vitesse apparente Va de l'onde correspondante est grande.

VIII - 2 - 1 Rapport entre la Vitesse Apparente Va et la Vitesse

Reelle de Propagation V0 en Fonction de l'Angle
d'Emergence "i" du Rayon Sismique

Les trajectoires arrivant en S1 et S2 sont considerees comme paralleles a cause

de Dx petit: on a par consequent :

sin i =V0 .Dt/Dx = IS2 / S1S2

121
 122

avec:

Va = Dx/Dt

d'ou la formule : Va =V0 / sin i Þ ce qui entraine Va ³ V0 .

Le signal correspond generalement aux reflexions.Pour un emetteur et un

recepteur donne,le schema ci-contre montre que l'angle d'emergence "i" du signal
sera d'autant plus faible que le miroir,sur lequel ce signal se sera reflechi,sera
plus profond.Donc Va tend vers l'infini pour des reflexions sur des miroirs tres
profonds.Au contraire,la plupart des bruits organises se propagent au voisinage de
la surface du sol,en ondes dispersives:Ondes de Love et surtout Ondes de
Rayleigh.Leurs angles d'emergence seront par consequent tres grands,leurs
vitesses apparentes faibles,presque egales a leur vitesses reelles
respectives.Nous vous presentons dans cette etude quelques types de bruits
connus dans le domaine petrolier.

VIII - 3 - Definition des Bruits

Un Bruit est un signal indesirable,son identification est une des operations les plus
delicates,car il se manifeste en premier lieu aleatoirement,et non organise,quand il
est previsible et coherent son etude devient facile pour le geophysicien.

VIII-3-1- Bruits Aleatoires:

Un bruit est dit aleatoire lorsqu'il n'est pas predictible,dans l'espace ou dans le
temps,en un mot,il ne presente aucune coherence avec les autres signaux.
Ce sont les microseismes dus aux manifestations: naturels (vent,pluie
,mouvements de la vegetation ou du terrain,ruissellement);activites humaines (train
, camion , pompe , influence des lignes electriques promeneurs ...); ou bien bruits
d'appareillage ou d'instrumentation ;sur lesquels on ne peut avoir d'action directe
sauf arreter toute activite a proximite du dispositif d'enregistrement,et aussi le
fonctionnement de toutes appareils.

122
 123

VIII-3-2 - Bruits Semi - Aleatoires

Ce sont des bruits consecutifs crees consecutivement a l'explosion , mais non

systematiques en particulier les debourages de trous,retombee de tube,fuite
d'animaux .

VIII -3-3 - Bruits de Diaphonie

L'influence par induction des fils vehiculant les signaux sismique (ou signaux
electriques) lorsque ces derniers sont soumis a des oscillations dans le champ
electromagnetique naturel , ou bien lorsqu'ils sont mutuellement proches les uns
des autres; est appelee "diaphonie".
Il est possible de dresser le tableau suivant,qui donne une premiere approximation
;distingue les signaux des bruits.

f Va k

BRUITS inferieur a presque egal a superieur a 10 c/Km

50 c/s la V reelle de 25 c/Km
50 a 500 m/s

SIGNAUX inferieur a plusieurs inferieur a 10 c/km

150 c/s milliers de m/s

Mesurant les caracteristiques des bruits et des signaux sur le terrain,le tir de bruit
permet d'ameliorer le rapport signal/bruit par un dispositif adequat
d'enregistrement.

VIII - 4 - Filtrage Spatial ou ( Filtrage en Nombre d'Onde ).

La notion du filtrage spatial est liee a la multiplication de geophone ,en effet la

repartition des geophones sur une surface donnee permet d'augmenter la
selectivite des signaux ayants meme directivite et meme phase et ellimine ainsi les
signaux non coherents et non organises.Le filtrage spatial ou en nombre d'onde
consiste a determiner les caracteristiques d'une nappe de geophones disposee sur
la surface du sol.Ces caracteristiques sont :la distance entre geophone,le nombre
de geophone,enfin la geometrie de la nappe.

123
 124

Considerons un cas simple ou l'on a deux recepteurs R1 et R2 espaces de "d";soit

un bruit sinusoidale de pulsation w ou de frequence "f" et d'amplitude "A" se
propageant horizontalement dans le milieu superficiel,suppose homogene et
isotrope,de Vitesse "V" La pulsation w et la frequence f sont lies par la relation :
w=2pf ; Le signal resultant enregistre s'ecrit :

y = A.cos(wt-j)+A.cos(w(t-Dt)-j)

avec :Dt=d/V

(Dt represente le retard en temps entre les recepteurs R1 et R2 ;

y = A.cos(wt-j)+A.cos(wt-wd/V)-j)

Þ y = 2A.cos(wd/2V).cos(wt-wd/2V-j)

etant donne que:

cosa+cosb=2cos(a+b)/2.cos(a-b)/2

L'amplitude du signal resultant "(2A.cos(wd/2V)) " est fonction de la distance "d"

entre les deux recepteurs R1 et R2. Nous cherchons par cette formule les lieux ou
le signal enregistre par les deux recepteurs est fortement attenue,donc l'amplitude
du signal resultant de l'expression precedente doit etre egale a zero:

wd/2V = p/2 + mp Þ d =pV/w + (2V/w).mp

etant donne que:

w=2pf et mÎN ;d = pV/2pf + (2V/2pf).mp

Þ d = V/2f + V.m / f = V/f(m+1/2)

Etant donne que :

124
 125

l = V/f Þ d = l(m+1/2) .

Pour un bruit de longueur d'onde l,il sera completement attenue si la distance entre les deux
recepteurs est de:

l/2;3l/2;5l2;7l/2;9l/2;11l/2..........(2m+1)/2l

ou: m=0,1,2,3,...
Dans le cas ou on cherche une attenuation optimale d'un bruit ;le terme ou l'argument wd/2V =
m'p ; d'ou d = m'l ; c'est-a-dire la distance entre les 2 recepteurs doit etre egale a
:0;l;2l;3l;5l;...etc...;.Cette exemple ellustre bien le principe de filtrage spatial ou en nombre
d'onde.

VIII - 4 - 1 - Calcul de la Fonction de Transfert

Si on admet que le signal S(t),resultant de la somme des differents signaux recus au niveau
des recepteurs repartis sur la surface du sol.L'expression mathematique de ce signal s'ecrit de
la facon suivante:

S(t)=s1(t) + s2(t) +s3(t) +.... +si(t) +..... +SN(t) .

Pour 2 recepteur elementaires voisins distants de "d":

Sj+1(t) =Sj(t -Dt) = S(t - jDt) ;jÎ[1,N-1]

ou :

Dt = d.sin e / V = d / Va

Dt represente le retard entre deux recepteurs voisins: La relation :

Sj+1(t) =Sj(t -Dt) Þ S(t) = s(t)+s(t-Dt)+s(t-2Dt)+...... + s(t-(n-1)Dt)

Þ S(t) = s(t)*[d(t)+d(t-Dt)+d(t-2Dt)+..............+d(t-(n-1)Dt)]

125
 126

Þ S(t) = s(t)*Sd(t-iDt)

De la " meme maniere" on represente le signal emis :

Þ R(t) = r(t)*Sd(t-jDt')

*:represente l'operation de deconvolution ; d(t):represente la fonction de Dirac .

Ce qui signifie que dans les deux cas ,les signaux S(t) et R(t) ont ete filtres respectivement par
un systeme dont la reponse impulsionnelle est:

h(t) = Sd(t-iDt) et g(t) = Sd(t-jDt')

Le filtrage spatial est donc un filtrage lineaire qui se reduit a une convolution dans l'espace. Soit
h(t) la reponse impulsionnelle de la trace;par transformee de Fourier on obtient:

HN( f ) = 1 + e(- 2pjfDt)+ e(- 2pjf2Dt)+....+e(- 2pjf(n-1)Dt)

On remarque que : HN( f ) est une progression geometrique de raison "q" tel que :

q=e(- 2pjfDt)

donc :

HN( f ) = (1 - qN ) / (1 - q ) = 1 - e( - 2pjfNDt) / 1 - e( - 2pjfDt)

ÞH (f)=[e( -pjfNDt)/e( - pjfDt) ][(e(pjfNDt)-e( -pjfNDt))/(e(pjfDt) -e( -pjfDt) )]
N

Þ HN( f ) = e - pjfNDt / e- pjfDt . Sin NpfDt /Sin pfDt

La reponse frequentielle Pour une nappe contenant N geophones,s' ecrit:

HN( f ) = e - pjf(N-1)Dt . Sin NpfDt /Sin pfDt

126
 127

La fonction de transfert du dispositif s'ecrit:

| HN( f ) | = | Sin NpfDt /Sin pfDt |

ou :

| H(ka) | = | Sin Npdka /Sin pdka |

Le spectre de phase est egal a:

j = - p(N-1)f.Dt

et s'exprime en fonction du nombre d'onde ka

par:

j(ka) = - p(N-1)d.ka

le dephasage j correspond au retard temporel Dt tel que :

j = wDt = - p(N-1)d.ka = - p(N-1)d.f/Va

d'ou :

Dt = -(N-1)d/2Va

la distance :

-(N-1)d/2 : correspond au centre de gravite de la ligne. Dt : represente alors le temps mis pour
parcourir la moitie de la longueur de la ligne a la vitesse Va. Cela veut dire que le signal
resultant enregistre etait le meme que s'il avait ete enregistre au barycentre de la ligne.

VIII - 4 -2 Etude de la Fonction de Transfert:

" | H(ka) | = 1/N | Sin Npdka /Sin pdka | "

127
 128

Soit N le nombre de Geophones identiques,places a intervalle "d" constante dont la reponse en

Amplitude est affectee au point milieu ou barycentre de la trace sismique. Dans ce cas on a
l'expression :

| H(ka) | = | Sin Npdka /Sin pdka |

de periode 1/d ,maximum pour k=n/e ; et nul pour k=n/Nd ou : n entier ; la phase,de periode
1/d pour N impaire ou 2/d pour N pair,est constamment egale a 0 ou p dans chaque intervalle
n/Nd et n+1/Nd. Le rapport presente une indetermination qu'on peut lever par les fonctions
equivalentes en utilisant la relation " sin x » x " quand x ® 0 ; c'est-a-dire en utilisant : l'
approximativement par continuite de " sin x / x " en zero. " Sin Npdka / Sin pdka " est
equivalent a : "Npdka / pdka = N ";On obtient : l'expression du module maximal de la courbe de
reponse en amplitude,pour :ka= 0 ,ou ka =1/d .
Le module de la fonction de transfert norme est :

| H(ka) | = 1/N | Sin Npdka/Sin pdka |.

128