49
III.3. Traitement des données
Nous avons récupéré les données des enregistrements, via la carte flash, sur PC avec le
logiciel ReadCity, et traité les signaux à l’aide du logiciel libre geopsy (figure 23) qui a été
développé en marge du programme européen SESAME (Site Effets Assessment using
Ambiant Excitations), par Wathelet (www.geopsy.org) qui permet de calculer en particulier
les rapports H /V.
Ce logiciel 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être) suivant, par exemple, la technique
de Konno et Ohmachi (1998), avec un paramètre de lissage b égale à 40 ;
5. Calcul de la moyenne quadratique des spectres horizontaux lissés Sh (Ni)
(composantes Nord-Sud 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 ;
III.3. Traitement des données
Nous avons récupéré les données des enregistrements, via la carte flash, sur PC avec le
logiciel ReadCity, et traité les signaux à l’aide du logiciel libre geopsy (figure 23) qui a été
développé en marge du programme européen SESAME (Site Effets Assessment using
Ambiant Excitations), par Wathelet (www.geopsy.org) qui permet de calculer en particulier
les rapports H /V.
Ce logiciel 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être) suivant, par exemple, la technique
de Konno et Ohmachi (1998), avec un paramètre de lissage b égale à 40 ;
5. Calcul de la moyenne quadratique des spectres horizontaux lissés Sh (Ni)
(composantes Nord-Sud 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 ;
