Page | 34
IV.7. Traitement des données
Nous avons récupéré les données des enregistrements, via la carte flash, sur PC afin de traiter les
signaux à l'aide du logiciel geopsy (www.geopsy.org) développé par Marc Wathelet, en marge
du programme européen SESAME (Site Effects Assessment using Ambiant Excitations),ce
programme permet, en particulier,de calculer les rapports H/V (Figure 18). L’analyse des signaux
a été faite sur des fenêtres stationnaires de 32 secondes.
Ce logiciel (Figure 3) traite les données de la façon suivante :
1. sélection des trois composantes (verticale v, nord-sud n, et est-ouest e) d’un signal de bruit de
fond sismique ;
2. recherche des Ni fenêtres stationnaires simultanément sur les trois composantes, à l'aide d'un
anti-trigger LTA/STA (LTA : Long Term Average ; STA : Short Term Average).
3. calcul individuel des spectres de chaque composante Sv (Ni.), Sn (Ni) et Se (Ni) sur chacune
des Ni fenêtres par FFT (Fast Fourier Transform);
4. lissage des spectres (pour chacune des Ni fenêtres) suivant, par exemple, la technique de
Konno et Ohmachi (1998), avec un paramètre de lissage b égal à 40;
5. calcul de la moyenne quadratique des spectres horizontaux lissés Sh (Ni) (composantes nordsud et est-ouest) pour chaque fenêtre :
6. calcul des rapports spectraux H/V (Ni) pour chaque fenêtre:
7. calcul de la moyenne géométrique du rapport spectral H/V sur toutes les fenêtres:
Sh(Ni)=
=
Sn(Ni)
2 + Se(Ni)
2
2
H / V (Ni) =
Sv(Ni)
Sh(Ni)
H/V= 10
∑ Iog 10 (H / V (Ni))
Ni