L3Pro - LPG2 - Methodes sismiques

Clara Levy, Anne Sieminski

UJF - OSUG - LGIT

[email protected]
[email protected]

1/67
Cours 1 : Introduction a la sismique refraction

2/67
Methodes de prospection geophysique

3/67
Un analogue des ondes sismiques

4/67
Propagation en milieu homogene

5/67
Propagation en milieu homogene

6/67
Vitesses sismiques
Milieu vP (m/s) vS (m/s) (g/cm3 )
Roche alteree 300-700 100-300 1.7-2.4
Sable sec 400-1200 100-500 1.5-1.7
Sable humide 1500-4000 100-500 1.9-2.1
Argiles 1100-2500 200-800 2.0-2.4
Marnes 2000-3000 750-1500 2.1-2.6
Gres 3000-4500 1200-2800 2.1-2.4
Calcaire 3500-6000 2000-3300 2.4-2.7
Craie 2300-2600 1100-1300 1.8-2.3
Salt 4500-5500 2500-3100 2.1-2.3
Granite 4500-6000 2500-3300 2.5-2.7
Basalte 5000-6000 2800-3400 2.7-3.1
Charbon 2200-2700 1000-1400 1.3-1.8
Eau 1450-1500 - 1
Glace 3400-3800 1700-1900 0.9
Petrole 1200-1250 - 0.6-0.9
7/67
Vitesses sismiques
Milieu vP hors nappe (m/s) vP sous nappe (m/s)
Sable fin 300-700 1450-1700
Graves 500-900 1700-2300
Argile 500-1400 1400-1700
Marnes 1800-2100 2100-3000
Gres 800-3000 2000-4000
Craie fracturee 800-1500 1700-2300
Craie saine 1800-2500 2300-3200
Calcaire fracture 900-2000 1700-3000
Calcaire compact 3000-5000 3500-5000
Granite fracture 1000-2500 2500-4500
Granite sain 4500-5500 4500-5500
Basalte fracture 800-2500 2300-4500
Basalte sain 5500-6000 5500-6000
Schiste altere 500-2000 1700-2500
Schiste sain 2500-4500 3000-4500
8/67
Representation de la propagation

9/67
Les ondes de la sismique refraction

10/67
Preparation a la manip de terrain

Une interface horizontale

Onde directe
x
td (x) =
v0
Droite de pente 1/v0 Hodochrone t(x)

11/67
Preparation a la manip de terrain
r
Une interface horizontale v1 + v0
xd = 2h0
v1 v0

Onde conique
s
x 1 1
tc (x) = + 2h0 2
2
v1 v0 v1

Droite de pente 1/v1 pour

x xc = 2h0 tan ic
Hodochrone t(x)
12/67
Preparation a la manip de terrain

Discussion sur le materiel et ecriture de la check list

13/67
Cours 2 : Retour sur la manip de sismique refraction

14/67
Rappels des episodes precedents

Premier cours
I Principe de la methode : temps darrivee dondes sismiques
artificielles image des proprietes mecaniques du sous-sol
(structure, nature, contraste)
I Ondes sismiques : propagation dune deformation elastique,
ondes P et S (+ ground roll), vitesse selon la nature de la
roche, plusieurs trajets possibles entre source et recepteur
I Sismique refraction (basee sur les ondes coniques, temps de
parcours suivant le rai sismique, hodochrones, quasi-directe)

Terrain (avec Mai-Linh et Clara)

I Acquisition de sismique refraction (manipulation et
deploiement du materiel, experience des difficultes
dacquisition)
15/67
Rappels des episodes precedents

16/67
Programme du jour

Retour sur lexperience de terrain

I Explication et mise au clair du protocole experimental
(notes de terrain)
I Interpretation des donnees avec le logiciel sardine

17/67
Geometrie dacquisition

Profils lineaires de 24 capteurs (milieu 2D sinon ?)

Plusieurs tirs/sources (stack, differents emplacements, source en

onde P et S)

Longueur du profil ?
I Elle doit etre assez grande pour voir londe conique arriver en
premier sur la moitie du dispositif (pour les deux tirs en bout)
I De lordre de 6h0 (a verifier sur le terrain)

Intertrace ? (ni trop petit, ni trop grand echantillonnage spatial)

18/67
Distance du point de brisure

Avec qqles a priori sur le milieu,

on peut estimer xb
r
v1 + v0
xb = 2h0
v1 v0

Exemple :
pour h0 = 10m et v1 = 4500 m/s
v0 [m/s] 400 2000 3000
xb [m]

Conclusion : on choisira une

longueur de profil L = 2xb ' 6h0

19/67
Visualisation

Premier tir en onde P

Stack : justification, origine du bruit, rapport signal sur bruit (SNR)

Ps
SNRN = SNR N avec SNR = 10 log
Pn
Z t2
1
et Ps = |s(t)|2 dt
t2 t1 t1

Pourquoi une plaque sous le marteau ?

Fonction source : forme du signal produit par la source, doit etre

proche dune impulsion pour une interpretation facile.

20/67
Transformee de Fourier

Pour analyser le contenu frequentiel dun signal (i.e. faire ressortir

les differentes frequences que contient le signal)

Frequence (en Hz) = 1/T (periode en s)

TF directe (avec t le temps) :

Z +
TF [s(t)] = S() = s(t)e i2t dt

TF inverse :
Z +
TF 1 [S()] = s(t) = S()e i2t dt

21/67
Decomposition dun signal

22/67
TF dun signal periodique

23/67
Echantillonnage

Pas Te [s] et frequence e = 1/Te [Hz] dechantillonnage

Theoreme de Shannon-Nyquist : e 2max
Filtrage anti-aliasing

24/67
Contenu frequentiel

Le contenu frequentiel dun signal sismique depend de la fonction

source, de la propagation dans le sol (attenuation) et de la
frequence propre du capteur (systeme oscillant amorti)

25/67
Contenu frequentiel
Le contenu frequentiel dun signal sismique depend de la fonction
source, de la propagation dans le sol (attenuation) et de la
frequence propre du capteur (systeme oscillant amorti)

Geophone Mark Products L-28 10 Hz

26/67
Contenu frequentiel

La resolution des images (i.e. precision des images, epaisseur des

couches en sismique refraction) depend du contenu frequentiel

Londe est sensible au milieu sur une zone de taille finie controlee
par la longueur donde [m] avec = vT (londe voit le milieu
au travers de )

La resolution est de lordre de /4

Precision ou profondeur de penetration (les HF sattenuent plus)

27/67
Contenu frequentiel

Filtrage : on enleve du signal certaines frequences (anti-alisaing,

rejection du 50 Hz)

Traitement du signal (TF, echantillonnage, filtrage)

28/67
Tir aller (A)

Hodochrone de londe conique :

x 2h0 v0
tc (x) = + cos ic avec sin ic =
v1 v0 v1
Intercept et delai sismique

Hypotheses ?

Desaccord avec les autres methodes geophysiques ?

29/67
Decouverte de sardine

Visualisation et discussion des traces en temps et en frequence

(bruit, stack)
Pointe des premieres arrivees
Trace des hodochrones, lien avec la theorie vue precedemment
Visualisation dun tir aller uniquement

30/67
Cours 3 : Interpretation des donnees de sismique
refraction

31/67
Episodes precedents

Sismique refraction

I Principe de la methode : source controlee, temps des

premieres arrivees (ondes directes et coniques), hodochrones
I Decouverte du materiel et acquisition des donnees sur le
terrain
I Rappels des conditions experimentales et visualisation des
donnees (sardine)

32/67
Programme

Sismique refraction (fin)

1. Comment choisir le capteur ? (discussion-cours)

2. Interpretation des donnees Modele en vP et vS (TP)
I Pourquoi plusieurs tirs ?
I Discussion du modele final
3. A quoi ca sert ? (cours)

Dernier cours
I Autres methodes sismiques (cours) : sismique reflexion,
tomographie, ondes de surface, rapport H/V

33/67
Sismometrie

Comment choisir le capteur sismique ?

eost.u-strasbg.fr/musee

34/67
Sismometrie
Par sa frequence propre 0
Capteur sismique = systeme oscillant avec 0 k, l

Resonance lorsque la frequence du signal excitant est proche de 0

En sismometrie, on cherche la resonance 0 ' sol
(6= risque sismique !)
35/67
Sismometrie

Le constructeur donne 0 et la fonction instrument F ()

Geophone Mark Products L-28 10 Hz

36/67
Sismometrie
Capteur = filtre F () reponse du capteur en fonction de

F () depend de 0 et de lamortissement h

Pourquoi amortir ? Comment obtenir F () ? (calcul, calibration)

37/67
Sismometrie

F () = |F ()|e i() Le signal de sortie y (t) est decale dans le

temps par rapport au signal dentree x(t) (dephasage)

38/67
Sismometrie

Le bon capteur est celui dont la reponse instrumentale est grande

sur la bande de frequence des mouvements du sol a enregistrer...

Quelle est la bande de frequence des mouvements du sol ?

39/67
Sources actives
Des ondes de frequences differentes selon la source

40/67
Sources actives

Des ondes de frequences differentes selon la source

Comment retrouver le contenu frequentiel de la source ?

41/67
Transformee de Fourier

42/67
Sources actives

Comment choisir la source ?

I Ce qui est possible (cout, technique, autorisation)

I Puissance (profondeur de la cible, longueur du profil)
I Precision de limagerie (HF pour bonne resolution et lecture
des differentes ondes, sismique reflexion, tomographie)

43/67
Propagation

Ne pas oublier les effets de la

propagation : les ondes HF sont
plus attenuees que les ondes BF
Lattenuation depend de la
nature du milieu

44/67
Sources naturelles

Microseismes (bruit oceanique et atmospherique) : 5 s < T < 20 s

Seismes tectoniques (sismologie) :
I a proche distance T ' 1 s
I a grande distance T ' 10 s - 100 s
(pourquoi different ?)
Et pour les deplacements avec T > 1000 s ?

45/67
Capteur electromagnetique
Le signal de sortie dun geophone est la vitesse de deplacement v
(une bobine qui bouge dans un champ magnetique courant
dinduction tel que i v )

46/67
Capteur electromagnetique

du
v= dt Le signal v (t) est plus HF que u(t)

Comment passer de v (t) a u(t) ?

47/67
Capteur large-bande

Bcp delectronique, systeme asservi (feedback), mesure de la force

pour ramener la masse acceleration, plus de 0 large bande
de frequence, mais cher et reglage delicat

48/67
Accelerometre
Capteurs electromagnetiques regles pour supporter de forts
mouvements du sol ( 2g )
Pourquoi ?
Comment ?

0 >> sol + fort amortissement

zsol
(z + 2z + 02 z = zsol )

49/67
En plus du capteur

Sismometre = capteur + horloge + amplificateur-numeriseur +

stockage + alimentation

Bon couplage avec le sol, environnement calme

50/67
Interpretation

Objectifs :
I Modele en vP et vS sous les deux profils
(24 geophones, x=3 m)
I Estimation des incertitudes
I Critique du resultat

Telecharger sardine :
www-lgit.obs.ujf-grenoble.fr/sieminsa/teaching/l3pro
Quels fichiers analyser ?
Interpretation conjointe de plusieurs tirs

51/67
Rappels

Onde directe
x
td (x) =
v0
Onde conique
x 2h0
tc (x) = + cos ic
v1 v0
v0
avec sin ic =
v1

Hyp : une interface

horizontale plane avec v1 > v0

52/67
Plusieurs tirs

Tir aller (A) et tir retour (B)

I Fermeture en temps
I Pendage de linterface

Tir centre (C)

I Confirme A et B + controle de la geometrie

Tirs lointains (O et P)
I Identification de la conique la plus profonde (toit du
substratum ?)

53/67
Fermeture de temps

En premier, verifier que t(xAB ) = t(xBA )

Gros problemes sinon (acquisition, pointe)

54/67
Pendage de linterface

(A) (B)
Si v1 6= v1 Pendage de
linterface
La vitesse mesuree est une
vitesse apparente dependante
du pendage

Ondes coniques montante et descendante :

x x
tup (x) = + Iup tdown (x) = + Idown
vup vdown

vup , vdown v1 , Iup , Idown profondeur sous les sources

55/67
Pendage de linterface

Expressions completes :
x x
tup (x) = + Iup tdown (x) = + Idown
vup vdown
v0 v0
vup = > vdown =
sin(ic ) sin(ic + )
2hup cos ic 2hdown cos ic
Iup = > Idown =
v0 v0
56/67
Pendage de linterface

Comment retrouver v1 et si << 1 ?

1 1 1 1 v1 1 1
' ( + ) ' ( + )
v1 2 vup vdown 2 vup vdown

57/67
Plusieurs interfaces horizontales

Si plusieurs brisures de pentes Plusieurs interfaces avec v

Hodochrone de londe conique sur la neme interface :

n1
(n) x X 2hi (in) (in) vi
tc = + cos ic avec sin ic =
vn vi vn
i

58/67
Couche cachee et non-unicite

Info a priori sur le milieu

(geologie, forage)
Combiner differentes
methodes geophysiques

59/67
Variations laterales de vitesse

Il faut des methodes plus sophistiquees !

I Sismique refraction interpretee avec la methode plus-minus
I Sismique reflexion ?
I Tomographie sismique (3D)

60/67
A quoi ca sert ?

A petite echelle, la sismique refraction sert a caracteriser le

sous-sol avant limplantation douvrages dart
I Geometrie et evolution du socle rocheux (vP )
I Proprietes mecaniques, etat physique des differentes couches
(vP et vS )

Possibilite de detecter et caracteriser des aquiferes (vP )

61/67
Qualite des roches et vP

Le RMQ est dautant plus faible que la roche est fracturee, alteree
et contient peu deau
Relation lineaire entre RMQ et vP

62/67
Rippabilite et vp

Table Caterpillar pour

tracteur D9G (390CV)
avec defonceuse 9B une
dent

Une roche alteree (vP faible) est plus facilement rippable

63/67
Caracterisation des aquiferes

Bon contraste de vitesse sur le toit de laquifere (vP augmente

avec la saturation en eau) Profondeur et geometrie de la nappe
64/67
Qualite des sols et vS

Relation entre compaction (penetrometrie) et vS

La compaction et la fracturation sont capitales pour la stabilite des

sols (glissement de terrain, liquefaction, effets de site)

65/67
Alea sismique

Vitesse de
reference :
P
hi
vS,30 = P i
i i /vi
h

(est-ce bien
pertinent ?)

66/67
Etude de la croute

A grande echelle, etude de la croute terrestre (topographie du

Moho, interfaces intracrustales)

Formation des chanes de montagne, ...

67/67