Sismique réfraction
Nous pouvons retrouver les variations de
vitesse et les épaisseurs des couches d’un
milieu tabulaire grâce aux tempsdes premières
arrivées sismiques.
1
Rai sismique
Sismique réfraction
On réduit l’information d'une trace à seulement un temps d’arrivée, en pointant la position de la première
ondelette visible.
3
Rai sismique
Identification des premières arrivées
La précision des pointés esttrès importante, car elle influence directement la précision des
résultats.
Des données trop bruitéesrendent difficile, voire
impossible, un pointage précis!
4
Rai sismique
Sismique réfraction
Bruits ou arrivées ?
Pour y parvenir, il faut des données de qualitécomportant des
premières arrivées franches.
5
Rai sismique
Front d’onde
L’identification des premièresarrivées
consiste à tracer dans l’espace le
déplacement du front d’onde parvenant
aux récepteurs. Il faut donc identifier la
même phase sur toutes les traces!
Front d’onde: Positions spatiales d’une
onde dont laphase est équivalente (par
exemple le minimum ou le maximum
d’une onde).
6
Rai sismique
Front d’onde
L’identification des premièresarrivées
consiste à tracer dans l’espace le
déplacement du front d’onde parvenant
aux récepteurs. Il faut donc identifier la
même phase sur toutes les traces!
Front d’onde: Positions spatiales d’une
onde dont laphase est équivalente (par
exemple le minimum ou le maximum
d’une onde).
7
Rai sismique
Principe de Huygens
Pour estimer les vitesses sismiquesdu sol, il faut
estimer quel a été le parcours du front d’onde qui
est arrivé à la surface. Comment ?
Principe de Huygens: Chaque point du front
d’onde agit comme une source d’une onde
sphérique, de même amplitude et de même phase.
L’interface constructive et destructive des sources
secondaires produit un front d’onde continu et quise
déplace dans le temps.
8
Rai sismique
Principe de Huygens
Pour estimer les vitesses sismiquesdu sol, il faut estimer
quel a été le parcours du front d’onde qui est arrivé à la
Front de l’onde Front de l’onde
surface. Comment ? réfléchie incidente
Principe de Huygens: Chaque point du front
d’onde agit comme une source d’une onde
sphérique, de même amplitude et de même phase.
L’interface constructive et destructive des sources
secondaires produit un front d’onde continu et qui
se déplace dans le temps.
nt de l’onde
réfractée
9
Rai sismique
Rai sismique
Il est plus commode d’utiliser le rai sismique pour
retracer le parcours etla direction du front d’onde.
Front de l’onde Front de l’onde
Le rai est perpendiculaire au front d’onde.
réfléchie incidente
Rai
Rai sismique: Vecteur perpendiculaire au front incident
d’onde etpointant dans sa direction de propagation.
Rai
réfléchi
Rai
réfracté
Front de l’onde
réfractée
10
Rai sismique
Lois de Snell-Descartes
Il est possible de prédire le parcoursdes rais sismiques grâce à la loi de Snell-
Descartes:
1. Loi de la réflexion: Rai
incident
Rai
réfléchi
2. Loi de la réfraction:
Rai
réfracté
11
Rai sismique
Réfractions en sismique
Attention! En sismique, nous avons plusieurs types d’ondes (P et S). Uneinterface produira
une conversion demode, et chaque type d’onde répond à l’équation de réfraction
Onde SV
généralisée: Onde P réfléchie
incidente
Onde P
réfléchie
Onde SV
réfractée
Onde P
réfractée
où p est le paramètre de rai.
12
Rai sismique
Réfractions SH
Les ondes SH (polarisée horizontalement, parallèle
au plan deréfraction) ne produisent pas de
conversion d’onde!
Onde SH
incidente
La loi de la réfraction restecependant la même!
Onde SH
réfléchie
Onde SH
réfractée
13
Rai sismique
Gradient de vitesse
Qu’arrive-t-il si la vitesse augmentegraduellement ?
Aucune réflexion n’est produite, maisle paramètre de Rai
rai reste toujours constant! incident
Rai
Les rais deviennent alors courbes etpeuvent remonter
réfléchi
à la surface!
Rai
réfracté
14
Rai sismique
Réfraction critique
Comment varie l’angle en fonctiondu contraste de vitesse entre les
milieux 1 et 2 ?
Prenons V2/V1 = 5, alors:
Réfraction
• critique
•
15
Rai sismique
Réfraction critique
L’angle critique est donné par:
• La réfraction critique ne survientque si
V2>V1
• L’angle se propage à la vitesse V2 lelong de Réfraction
l’interface critique
• Selon le principe de Huygens, uneonde est
réémise vers le haut, à unangle égal à l’angle
critique.
16
La réfraction
• Modèle à deuxcouches
• Un gradient de vitesseest présent
dans chaque couche, ce qui
courbe les rais
• Propagation purementacoustique,
pas de conversion de mode, ni
d’ondes de surface
17
Jean Virieux after Anne Obermann
La réfraction
18
La réfraction
• L’angle critique estatteint
• La réfraction est maintenant
parallèle àl'interface
Angle critique
• Une réflexion totale interne est
maintenantatteinte, c’est-à-dire
que le coefficient de réflexion est
de 1.
19
La réfraction
20
La réfraction
• La réflexion atteint maintenant Onde conique
la surface
• On discerne maintenant une
onde « conique » causéepar la
réfraction critique montante
21
La réfraction Onde
Distance
• Les ondes réfractéespar le
critique
gradient de vitesse (les rais
courbes jaunes) arrivent
maintenant avant l’onde
directe
• À ce moment, la distance
critique dugradient est
atteinte
22
La réfraction
23
La réfraction
24
La réfraction
• La réfraction critique rejoint
maintenant l’onde réfractée par
legradient de vitesse
• Un brusque changement de
penteest observé, car l’onde
conique se propage à la vitesse
du milieu 2
• La distance de croisement de
la réfraction au milieu 2est
alors atteinte.
25
La réfraction
26
Jean Virieux after Anne Obermann
La réfraction
27
La réfraction
• À long déport, les ondes
réfléchies et lesondes courbées
par legradient de vitesse ont le Onde
guidée
même temps de propagation
• Une seule arrivée est alors
observée sur lesismogramme
28
La réfraction
29
Milieu tabulaire à 2 couches horizontales
30
2 couches horizontales
Parcours des rais sismiques
L’énergie peut voyager de la S R
source au receveur selon plusieurs parcours:
• Onde directe
• Onde réfléchie
• Onde réfractée
31
2 couches horizontales
Réciprocité
Réciprocité en sismique: Les S R
temps d’arrivée entre deux points
fixes sont équivalents, peu importe siles positions
du tir et du receveur sont interverties.
• Ceci permet de vérifier la précision du pointage
des arrivées (tSR = tRS)
• Permet de vérifier des erreurs sur le
positionnement du temps 0 lors del’acquisition
• Les méthodes d’interprétation sebasent sur
le principe de réciprocité
32
2 couches horizontales
Dromochronique
Temps
La dromochronique est le graphique des temps d’arrivées
en fonction de ladistance. Elle permet d’interpréter les
levés de sismique réfraction.
Distance
S R R R
33
2 couches horizontales
Onde directe
Temps
L’onde directe voyage sur un parcours direct entre la source et
les receveurs. Son temps d’arrivée est:
Distance
S R R R
34
2 couches horizontales
Onde réfléchie
Temps
L’onde réfléchie subit une réflexion aupoint milieu de S et
R. Son temps d’arrivée est:
Distance
S R R
M
35
2 couches horizontales
Onde réfractée
Temps
L’onde réfractée suit une droite dont lapente est donnée par
l’inverse de la vitesse du milieu 2.
Distance
S R R
36
2 couches horizontales
Temps d’intercepte
Temps
L’ordonnée à l’origine de la droite dela réfraction
s’appelle le temps d’intercepte:
Distance
S R R
37
2 couches horizontales
Distance critique
Temps
Attention! Aucune réfraction ne peutarriver avant la
distance critique!
L’ordonnée à l’origine ne représentepas une arrivée
physique!
Distance
S R R R
38
2 couches horizontales
Distance de croisement
Temps
Une autre distance importante est la distance de
croisement, le point où l’onde réfractée rejoint l’onde directe:
Distance
S R
39
2 couches horizontales
Interprétation deux couches horizontales
1. Calculer les vitesses des milieux 1et 2 grâce à
l’inverse des pentesdes ondes directes et
Temps
réfractées
2. Estimer l’épaisseur soit:
• par les temps d’intercepte:
Distance
• par les distances de croisement
40
2 couches horizontales
Exemple: Interprétation à deux couches
Exemple: Interprétez les temps d’arrivées suivants selon un modèle à deux couches horizontales.
Distance (m) Temps (ms)
0 0.0
5 16.7
10 33.3
15 49.2
20 52.5
25 55.9
30 59.2
35 62.5
41
2 couches horizontales
Exemple: Interprétation à deux couches
Exemple: Interprétez les temps d’arrivées suivants selon un modèle à deuxcouches horizontales.
1. Tracez la dromochronique
42
2 couches horizontales
Exemple: Interprétation à deux couches
Exemple: Interprétez les temps d’arrivées suivants selon un modèle à deuxcouches horizontales.
1. Tracez la dromochronique
43
2 couches horizontales
Exemple: Interprétation à deux couches
Exemple: Interprétez les temps d’arrivées suivants selon un modèle à deuxcouches horizontales.
1. Tracez la dromochronique
2. Identifiez les brisures de pente et assignez les
arrivées à un réfracteur
44
2 couches horizontales
Exemple: Interprétation à deux couches
Exemple: Interprétez les temps d’arrivées suivants selon un modèle à deuxcouches horizontales.
1. Tracez la dromochronique
2. Identifiez les brisures de pente et assigner les
arrivées à un réfracteur
3. Obtenir les vitesses par l’inverse despentes V2 = 1500 m/s
V1 = 300 m/s
45
2 couches horizontales
Exemple: Interprétation à deux couches
Exemple: Interprétez les temps d’arrivées suivants selon un modèle à deuxcouches horizontales.
1. Tracez la dromochronique
2. Identifiez les brisures de pente et assignez les
arrivées à un réfracteur
t0 = 39 ms
3. Obtenir les vitesses par l’inverse despentes
4. Déterminez le temps d’intercepte et ladistance de V2 = 1500 m/s
croisement
V1 = 300 m/s
xi = 14.7 m
46
2 couches horizontales
Exemple: Interprétation à deux couches
Exemple: Interprétez les temps d’arrivées suivants selon un modèle à deuxcouches horizontales.
1. Tracez la dromochronique
2. Identifiez les brisures de pente et assignez les
arrivées à un réfracteur
t0 = 39 ms
3. Obtenir les vitesses par l’inverse despentes
4. Déterminez le temps d’intercepte et ladistance de V2 = 1500 m/s
croisement
5. Déterminez la profondeur sous le tir:
V1 = 300 m/s
xi = 14.7 m
47
2 couches horizontales
Interprétation de n couches horizontales
Temps
Dans le cas de n couches horizontales, le temps d’arrivée est
donné par:
On voit que la pente dépend de la vitesse de la couche N. Distance
S R
48
2 couches horizontales
Interprétation de n couches horizontales
1. Identifier chaque bris de pente sur les dromochroniques, et en tirer les vitessesdes n couches ainsi que leur
temps et leur distance d’intercepte.
2. Trouver itérativement les épaisseurs de chaque couche:
Distances de croisement Temps d’intercepte
49
2 couches horizontales
Interprétation de n couches horizontales
3. Trouver le déplacement sous le tir pour lequel la profondeur est valide
S R
50
Milieu tabulaire avec une couche inclinée
51
Couche inclinée
Couches inclinées
Temps
Les vitesses apparentes du réfracteur dépendent
maintenant de la direction du tir!
Il faut toujours un tir direct et un tir inverse en
—
sismique réfraction ! t0
xi— Distance
A
B
52
Couche inclinée
Couches inclinées
Le temps de propagation entre A et B:
A B
hA
hB
53
Couche inclinée
Couches inclinées
Le temps de propagation entre A et B:
54
Couche inclinée
Couches inclinées
Le temps de propagation entre A et B:
55
Couche inclinée
Couches inclinées
Le temps de A à B est bien sur égal au temps de B à A.
Cependant, la vitesse apparente sera différente!
A x B1 B2 B3 B4
hA
hB
Le parcours réfracteur
56
surface se rallonge en x
Couche inclinée
Couches inclinées
Le temps de A à B est bien sur égal au temps de B à A.
Cependant, la vitesse apparente sera différente!
A1 A1 A2 A3 x B
hA
hB
Le parcours réfracteur
57 surface se raccourcit en x
Couche inclinée
Couches inclinées
Pour trouver les vitesses, il faut exprimer les hauteurs en A et B en fonction de x:
A1 x B
hA
hB
58
Couche inclinée
Couches inclinées
La vitesse est la dérivée de la distance en fonction du temps. Sachant:
La vitesse pour un tir en A:
59
Couche inclinée
Couches inclinées
Temps
Les vitesses apparentes sontdonnées par:
t—
0
Distance
A B
60
Couche inclinée
Couches inclinées
Comment déterminer le pendage et la vitesse vraie ?
61
Couche inclinée
Temps d’intercepte
Temps
Quels sont les temps d’intercepte ?
Le temps de parcours était:
t—
0
—
Pour les tirs en A et en B: xi Distance
62
Couche inclinée
Distances de croisement
Temps
Quelles sont les distances de croisement?
Le temps de parcours était:
t—
0
Pour les tirs en A et en B:
xi— Distance
63
Couche inclinée
Couches inclinées
Attention! Les hauteurs hA et hB sont les longueurs perpendiculaires à l’interfacesous les tirs en A et B! Pour retrouver les
profondeurs au droit des tirs, il faut:
hA
hB
64
Couche inclinée
Interprétation d’un réfracteur incliné
1. Déterminer les vitesses V1, V2+ et V2-
directement sur les dromochroniques
2. Déterminer l’angle critique et lependage à l’aide
de:
3. Déterminer la vitesse vraie du
deuxième milieu:
4. Déterminer les profondeurs grâce
aux distances de croisement ou aux
temps d’intercepte:
65
Corrections et limitations
66
Corrections et limitations
Correction topographique
La topographie change le parcours des rais. Pour la
source:
Après un peu d’algèbre:
e: Élévation
h: Profondeur de la sourced: Élévation de
référence
67
Corrections et limitations
Correction topographique
La correction topographique totale est la
somme des corrections de la source et du
géophone:
Connaissant:
es: Élévation de la source er: Élévation du
géophone
h: Profondeur de la sourced: Élévation
de référence
68
Corrections et limitations
Couche cachée
Que se passe-t-il lorsqu’une couche intermédiaire s’amincit ?
Temps
Distance
69
Corrections et limitations
Couche cachée
Que se passe-t-il lorsqu’une couche intermédiaire s’amincit ?
Temps
Distance
70
Corrections et limitations
Couche cachée
Que se passe-t-il lorsqu’une couche intermédiaire s’amincit ?
Temps
Distance
71
Corrections et limitations
Couche cachée
Temps
Une couche cachée est une couche dont laréfraction ne devient
jamais une première arrivée!
Les causes:
• Contraste de vitesse trop faible
• Épaisseur trop faible
• Espacement des géophones trop
grands
Distance
Dans ce cas, la profondeur du réfracteurvisible sera erronée!
72
Corrections et limitations
Couche cachée: correction
On peut corriger pour une couche cachéedans le cas simple
d’un milieu à 3 couches horizontales.
Facteur de correction
Pour cela, il faut connaître la vitesse V2 etune des épaisseurs.
V2 et hr = h1 + h2 V2 et h1 V2 et h2
73
Corrections et limitations
Inversion de vitesse
Temps
En présence d’une inversion de vitesse(V2 < V1):
• Aucune réfraction critique n’estpossible
• La dromochronique ne peut donc montrer un bris de
pente pour cettecouche
• L’interprétation sera un modèle à 2
couches
Distance
• La profondeur du réfracteur seraexagérée
La sismique réfraction n’est pasapplicable dans
ce cas!
74
Corrections et limitations
Faille
Temps
Une faille ou une dépression brusque duroc causera des
discontinuités de la dromochronique.
Lorsque la hauteur de la discontinuité est faible par rapport à la
profondeur du réfracteur, la hauteur peut être approximée par: Distance
75
observées sur les dromochroniques:
A. Variations topographiques rapides ou
changement localde vitesse de la
couche 1
B. Lentille de vitesse anormale
C. Topographie localisée duréfracteur
D. Changement de vitesselocale du
réfracteur
76
Corrections et limitations
Corrections et limitations
Bref, les dromochroniques sont rarement des droites parfaites en pratique, cequi rend l’interprétation non triviale!
Effet combiné de la
topographie de surface
et du roc!
77
