Chapitre 3 - Collecte Et Analyse Statistique Des Données

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 48

CHAP III : Collecte et analyse statistique des données

1. Bathymétrie et morphologie

2. Niveau d’eau projet

3. Calcul des extrêmes

4. Analyse fréquentielle de la houle


Collecte et analyse statistique des données

2
Collecte et analyse statistique des données

BATHYMÉTRIE ET MORPHOLOGIE

La topographie est la description des caractéristiques géométriques de la surface du sol au-dessus de l’eau. Trois
méthodes sont habituellement utilisées

Nivellement classique

Photogrammétrie est une technique d’analyse de


photographies aériennes qui permet de quantifier la
géométrie des objets dans les clichés par rapport à un
système de coordonnées fixes

Scannage laser en 3D : le balayage au laser


permet d’obtenir rapidement des données en 3D
(10 000 points par seconde) avec une densité
relativement élevée (espacement entre les points
d'environ 0.1 m voire mieux).
3
Collecte et analyse statistique des données

BATHYMÉTRIE ET MORPHOLOGIE

La bathymétrie est la description de la surface du sol sous l’eau. Il s’agit d’une condition aux limites particulièrement
importante pour le dimensionnement géométrique et structurel des ouvrages. Elle a une influence significative sur le volume
de dragage, le volume d’enrochement nécessaire ainsi que sur les actions hydrauliques. À la côte, par exemple, la hauteur
d’eau peut limiter la hauteur de la houle. Sur les cartes marines, le niveau du fond de la mer est défini par rapport à un niveau
de référence (p. ex. cote marine), généralement égal au niveau de plus basse mer astronomique du site.

La bathymétrie d’une zone est normalement déterminée à l’aide d’un bateau équipé d’un sondeur acoustique mono faisceau
(travaillant point par point) ou multifaisceaux (travaillant profil par profil).

4
Collecte et analyse statistique des données

BATHYMÉTRIE ET MORPHOLOGIE

Dans la pratique, ces levés sont effectués à partir d’un bateau, à l’exception des zones intertidales et en eau peu profonde
pour lesquelles il est possible d’employer des méthodes topographiques.

• la position horizontale (X, Y) du bateau, qui est généralement obtenue à l’aide d’instruments topographiques
classiques tels que le GPS différentiel (DGPS), le laser ou des dispositifs optiques;
• la hauteur d’eau (d), donnée par le sondeur ;
• la position verticale (Z) du niveau de l’eau qui est obtenue soit à l’aide d’équipements topographiques classiques (tels
que le DGPS, le laser, les dispositifs optiques), soit par la mesure de la marée ou à l’aide d’échelles référencées.

NB: Il convient de prêter une attention particulière à la


densité et à l’espacement des lignes de levé. Elles
doivent être suffisamment proches pour donner une juste
représentation des caractéristiques et du niveau du fond
sans que le coût soit excessif.

5
Collecte et analyse statistique des données

BATHYMÉTRIE ET MORPHOLOGIE

La morphodynamique est la description des processus de sédimentation et d’érosion qui se produisent sur le fond. À cause
de l’érosion et de la sédimentation, il est possible que la bathymétrie évolue avec le temps. Ces variations sont
problématiques lorsque l’on effectue des mesures pour déterminer la bathymétrie à des fins de dimensionnement ou de
construction.

La vitesse d’évolution morphologique du fond dépend de l'intensité du transport sédimentaire. Des taux de transport
élevés entraînent généralement une évolution relativement rapide du fond. Ces variations morphologiques sont toutefois
très lentes comparées aux variations des conditions aux limites hydrauliques ; en règle générale, elles s'étalent sur plusieurs
années voire des dizaines d'années, à l'exception de l'affouillement local près des ouvrages.

6
Collecte et analyse statistique des données

BATHYMÉTRIE ET MORPHOLOGIE LIÉES AUX OUVRAGES MARITIMES

Il est fondamental de connaître la bathymétrie et la morphologie du fond de la mer pour dimensionner des ouvrages
côtiers, par exemple dans le cas de hauteurs de houle limitées par la profondeur d'eau. Dans bien des cas, la connaissance
des variations morphologiques du fond au cours du temps est tout aussi essentielle au dimensionnement que la
détermination du niveau moyen du fond, dans la mesure où il faut déterminer le niveau du fond le plus bas devant
l'ouvrage, pour son dimensionnement.

De fait, que le fond de la mer soit en équilibre dynamique ou instable sur le plan morphologique, l’amplitude des variations
du niveau du fond à prendre en compte dans le dimensionnement doit être identifiée

7
Collecte et analyse statistique des données

Collecte et analyse des données (Houle et


niveaux d’eau)

Données calculées Données Mesurées

Modèle transformation Houlographes (Houle)


Marégraphe (Marée)

Wavewatch III
WAM : Puertos del esatados

8
Collecte et analyse statistique des données

Niveau d’eau projet

Le niveau d'eau de projet est la somme du niveau de la marée la plus haute, des décotes et des surcotes dues à des
conditions météorologiques particulières et de la valeur prévisible pour une période de retour de 100 ans de
remontée lente du niveau de la mer due aux altérations climatiques :

a) Niveau d’eau projet haut

b) Niveau d’eau projet bas

9
Collecte et analyse statistique des données

Niveau d’eau projet

a) Niveau d’eau projet haut

1. Moyenne des pleines mers des vives eaux (m/Zh)

2. Dans des conditions atmosphériques exceptionnelles, le vent peut provoquer des surélévations de niveau de
l’ordre de 0,20 à 0,30 m.

3. Dans la majorité des cas un vent soufflant à 30 m/s en provenance de la mer est associé à une dépression
barométrique de 990 mbar (donnant naissance une surcote de +0,25 m)

4. A partir de l'analyse systématique de la longue série d'observations marégraphiques disponibles pour certains
sites du Globe, il a été conclu que le niveau de la mer est dans une phase d'ascension, avec une tendance
d'environ 0.5 à 0.6 mm / an.

10
Collecte et analyse statistique des données

Niveau d’eau projet

a) Niveau d’eau projet bas

1. Moyenne des basses mers des vives eaux (m/Zh)

2. Dans des conditions atmosphériques exceptionnelles, le vent peuvent provoquer des abaissements de niveau de
l’ordre de 0,20 à 0,30 m.

11
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Les ouvrages en mer, construites généralement pour protéger les ports de la houle, sont soumises à son action. La
connaissance des valeurs extrêmes est donc un élément clé pour leur conception. La houle projet appelée aussi la
houle de dimensionnement est la hauteur de houle significative caractérisée par une période de retour X (dont X est la
durée de vie de projet).

La période de retour, ou temps de retour, caractérise le temps statistique entre deux occurrences d'un événement
naturel d'une intensité donnée. Ce terme est très utilisé pour caractériser les risques naturels.

Le concept de période de retour quand on traite de probabilités combinées est moins simple que lorsqu’on a affaire
à une seule variable. Un extrême de probabilité combinée peut être défini comme la probabilité qu’une hauteur de
la houle spécifique soit dépassée simultanément au dépassement d’un niveau d’eau donné.

12
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet


Cordonnées la fréquence Direction Temps
cartésiennes relative des vagues

Densité d’action des vagues

l'énergie donnée aux vagues par le vent

Fréquence relative

Opérateur différentiel

Groupe de la phase des houles

Vitesses de propagation Transfert d’énergie Transfert d’énergie


dans l’espace (x,y) en fréquence Angulaire
13
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Capteur de pression

Datawell

14
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Les houles, un phénomène aléatoire.

Analyse vague par Analyse spectrale


vague.
15
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

16
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Etape 2 : Analyse statistiques

Période de retour

Calcul de Hs Projet

Détermination de la période de pic et


de la Direction associée à HS

17
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

L’estimation statistique des événement tempétueux peut se faire de deux


manières

la méthode des la méthode des valeurs


maxima annuels supérieures à un seuil

la première, dénommée block maxima method (BM) et la seconde approche est désignée sous le terme Peaks-Over-
Threshold method (POT). En pratique, on se ramène dans les deux cas à un cadre de statistique paramétrique et
d’ajustement de paramètres sur les données. Selon le contexte, l’une ou l’autre des approches peut se révéler mieux
adaptée, mais il est le plus souvent utile de les mettre toutes deux en œuvre et d’en comparer les résultats

18
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

la méthode des maxima annuels

19
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

la méthode des maxima annuels

Cette méthode d’analyse que les valeurs maximales sur un intervalle de temps est donné souvent pour un an (maxima
annuels)

20
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Le modèle de dépassement de seuil (POT)

21
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Le modèle de dépassement de seuil (POT)

22
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Le modèle de dépassement de seuil (POT)

La théorie des valeurs extrêmes par méthode POT, propose différentes méthodes pour estimer un seuil au-delà
duquel une observation sera considérée comme valeur extrême. La sélection du seuil pour l'application de la
méthode Peaks Over Threshold (POT) à un seul échantillon est un sujet de longue date qui n'est toujours pas résolu
[Davison et Huser, 2015; Cavanaugh et al., 2015]. Bien qu'il y ait eu diverses propositions pour automatiser la
sélection des seuils [par exemple, Dupuis, 1998; Thompson et al., 2009; Solari et Losada, 2012, entre autres], la
pratique la plus courante consiste à utiliser des méthodes graphiques (voir Coles [2001].

La fonction moyenne des excès (FME), qui permet de décrire la prédiction du dépassement du seuil u lorsqu’un
excès se produit, est définie par :

23
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Le modèle de dépassement de seuil (POT)

Cette fonction moyenne des excès est estimée par la somme des excès dépassant un certain seuil élevé u, divisé par le
nombre d’observations qui dépassent ce seuil.

Dans le travail de (DAVISON & SMITH, 1990), l’étude de sensibilité au changement du seuil de sélection des évènements a
montré qu’à partir du seuil optimal, l’espérance résiduelle des excès devient approximativement linéaire. Cette approche
pratique visant à choisir le seuil u consiste à tracer l’estimateur empirique de l’espérance résiduelle et à choisir u de manière
à ce que ên(u) soit approximativement linéaire.
24
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Le modèle de dépassement de seuil (POT)

• la détermination d’un seuil à partir duquel un événement est considéré comme


atypique est basée sur la moyenne des excès qui représente la somme des excès
au dessus du seuil µ devisé par le nombre N des données

25
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

• L’un des inconvénients de la première procédure réside dans le choix d’un


événement non signifiant lorsqu’aucune tempête importante n’est enregistrée
durant une année et l’autre concerne la non prise en compte des événements
importantes enregistrées dans la même année.

• Dans le cas de séries courtes, d’une longueur comprise entre 10 et 20 ans, la


méthode des valeurs supérieures à un seuil est préconisée (Javelle, 2001).

26
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

• L’étape suivante consiste à constituer une série comportant suffisamment de


données pour réaliser l’ajustement d’une loi statistique (Dautrebande et al.,
2006).

• La prédétermination des événements tempétueux selon une loi de probabilité


est choisie parmi une diversité des lois (Meylan et al, 2008).

• Le choix de la loi et de la méthode d’estimation les mieux adaptées pour


l’ajustement de la distribution des évènement se fait graphiquement puis il est
confirmé numériquement par des critères statistiques.

27
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

28
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet Lois AIC BIC

Log.normal
La méthode graphique est cependant insuffisante pour choisir les modèles
411.4272 414.3586
adéquats à l’estimation. Il est nécessaire donc de faire intervenir les deux critères
statistiques AIC (Akaike Information Criterion) (Akaike, 1973) et BIC (Bayesian

Gamma Normal
Information Criterion) (Schwarz, 1978). 438.3114 441.2429

416.8131 419.7446
Ces deux critères permettent de classifier les lois selon une performance d’estimation.
La loi qui obtient les plus faibles valeurs des deux critères est la plus performante

Exp
431.2351 432.7009
pour la prédétermination des événements extrêmes.

Weibull
422.3308 425.2623
Ces critères se définissent de la façon suivante :
AIC = -2*logLik+2*K
BIC = -2*logLik+K*log(N)

Weibull
400.01 398.85
Où LogLik représente la log-vraisemblance des paramètres associés aux données,
K désigne le nombre de paramètres indépendants du modèle.
N est le nombre d’individus composant l’échantillon. 29
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Période de retour

Par définition, la période de retour T de la valeur x est l’espérance mathématique de la variable aléatoire définie
comme la durée séparant deux dépassements successifs (dépassements vers les valeurs supérieures) de cette
valeur (associée à une probabilité P. Cette définition rigoureuse permet pour un processus donné d’attribuer à
toute valeur x sa période de retour T(x) et, inversement, d’attribuer une valeur x(T) à toute période de retour T.

Cette notion est primordiale dans la détermination des houles extrêmes : si X désigne une hauteur significative
de houle, x(T = 100) désignera la valeur de la houle centennale, à condition que la période de retour soit
exprimée en années.

Cette remarque nous amène d’ailleurs à introduire un paramètre d’importance dans l’analyse des échantillons de
valeurs extrêmes par méthode de renouvellement (POT). Le taux moyen des évènements de tempêtes, ou mean
rate. Soit NT le nombre de tels évènements sur une période de K années. Alors, le taux moyen λT est défini
comme suit :

30
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Période de retour

La période de retour correspondant à fonction de répartition de


la loi F(x) dans le cas de la méthode des maximas annuels (Coles
2001) :

1
𝑇(𝑥) =
(1 − 𝐹 𝑥 )
La période de retour correspondant à fonction de répartition de
la loi F(x) dans le cas de la méthode de renouvellement (ou Sur une durée de vie planifiée de N années pour
méthode POT : Peaks Over Threshold-Pics au-dessus du seuil) est une structure (pas nécessairement la même que la
donnée par l’équation (Coles 2001) : période de retour de projet), la probabilité de
rencontrer la condition d'une houle caractérisée
1
𝑇(𝑥) = par une période de retour TR, au moins une fois,
𝜆𝑡(1 − 𝐹(𝑥)) est donnée par la distribution de Poisson:
Où 𝝀𝒕 est le nombre moyen d’évènements supérieurs au seuil
par an.

31
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Le modèle de dépassement de seuil (POT)

Il existe quelques logiciels pour pouvoir travailler sur les valeurs extrêmes sans tout fois programmer soi-même.
Exemple (HYFRAN).

Période de retour Proba au non dépassement Hs (m) Intervalle de confiance

100 0,99 9,5 8,27 - 10,7

50 0,98 8,79 7,73 - 9,84

20 0,95 7,85 7,01 - 8,70

10 0,9 7,13 6,44 - 7,82

5 0,8 6,38 5,84 - 6,91

3 0,6667 5,78 5,35 - 6,20

2 0,5 5,24 4,89 - 5,59

1 0,0001 2,64 2,12 - 3,16

32
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Après détermination des périodes de retour des Hs (Hauteur significative de la houle), il est nécessaire de tracer la
correspondance entre Hs/Tp (Période de pic) et Hs/Direction
20.00

18.00

16.00

14.00

12.00
Tp (s)

10.00

8.00

6.00

4.00

2.00

0.00
0 1 2 3 4 5 6 7 8 9
Hs (m)

33
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

Après détermination des périodes de retour des Hs (Hauteur significative de la houle), il est nécessaire de tracer la
correspondance entre Hs/Tp (Période de pic) et Hs/Direction

350

300

250
Direction (°N)

200

150

100

50

0
0 1 2 3 4 5 6 7 8 9
Hs(m)

34
Collecte et analyse statistique des données

Calcul des extrêmes : Houle projet

TP : ANALYSE DES EXTREMS : A partir des données ci-dessous, représentant les valeurs de Hs maximale annuelle
enregistrées sur la période 2005-2018, calculer les Hs associées aux période de retour :

• Tr = 10 ans
• Tr=50 ans
• Tr=100 ans

Nb: La loi adoptée pour le calcul des Hs Projets doit être justifiée

35
Collecte et analyse statistique des données

Analyse fréquentielle des données de la houle

Analyse statistiques

Analyse fréquentielle

Caractérisation des états de mer dans un milieu donnée

36
Collecte et analyse statistique des données

Analyse fréquentielle des données de la houle

La statistique constitue un outil indispensable pour effectuer un calcul de fréquence afin de tirer profit au mieux de
toutes les informations concernant les fréquences d’occurrences (Hs/Direction, Hs/Tp, Tp/Direction et
Direction/Hs/Tp.

En statistiques, le concept de fréquence absolue est le nombre d'apparitions (occurrences) d'une valeur dans une série
statistique. Il suffit de prendre la fréquence absolue et de la rapporter au nombre total d'éléments de la série. En
langage statistique, la fréquence relative s'obtient en divisant la taille de la classe par la taille de la population.

37
Collecte et analyse statistique des données

Analyse fréquentielle des données de la houle

L’analyse de la fréquence d’occurrence passe par 4 étapes.

1. Préparation d’un ensemble de données (encore appelé échantillon ou population). Dans notre
cas, ça concerne les paramètres Hs, Tp et direction.
2. Choix d’une valeur particulière ou plutôt d’un intervalle de valeur (ex : 1>Hs>1.25)
3. Calcul de la fréquence absolue ; à savoir le nombre d'apparitions (appelées aussi
occurrences) de la valeur en question ou plutôt des valeurs à l’intérieur d’un intervalle.
4. Calcul de la fréquence relative dans la série. Il suffit de diviser la fréquence absolue par la taille de
la série qui corresponde au nombre total d'éléments de la série.

NB : En additionnant les fréquences relatives de tous les éléments d'une série, on doit obtenir une
valeur de 100%.

38
Collecte et analyse statistique des données

Analyse fréquentielle des données de la houle

La largeur Δx (ΔHs pour la hauteur significative de la houle, ΔTp pour la période de pic et Δ𝜃° pour la direction de la
houle) de l'intervalle utilisé pour calculer les fréquences est choisi de manière à ce qu'un nombre suffisant d'observations se
retrouve à l'intérieur de chacun des intervalles de telle manière que les variations de l'histogramme soient suffisamment
douces.

Si le nombre d'observations "ni" est compris dans l'intervalle "i" des données [xi − Δx , xi ] , et si le nombre total
d'observations est "n", alors la Fonction de Fréquence Relative s'écrit : ni/N

39
Collecte et analyse statistique des données

Analyse fréquentielle des données de la houle

"ni" : le nombre d'observations


compris dans l'intervalle "i" des
données [xi − Δx , xi ]

" N " : le nombre total des


observations

Le pourcentage d’occurrence = ni/N

L’objectif est de déterminer les


condition les plus fréquentes
ni= 3 et N =14

40
Collecte et analyse statistique des données

Analyse fréquentielle des données de la houle

1. Code sur Excel :


1. Fréquence d’occurrence Hs/Direction :

=(NB.SI.ENS(B:B;">20";B:B;"<=40";C:C;"<="&Lx;C:C;">"&Ky)/$M$z)*100
Où ;

o B:B : Colonne de la direction de la houle.


o C:C : Colonne de la hauteur significative de la houle
o &Lx &Ky : Largeur de l’intervalle de la variable Hs
o $M$z: Nombre total d'éléments de la série
42
Collecte et analyse statistique des données

Analyse fréquentielle des données de la houle

1. Code sur Excel :


1. Fréquence d’occurrence Tp/Direction :

=(NB.SI.ENS(B:B;">20";B:B;"<=40";C:C;"<="&Lx;C:C;">"&Ky)/$M$z)*100
Où ;

o B:B : Colonne de la direction de la houle.


o C:C : Colonne de la période de pic.
o &Lx &Ky : Largeur de l’intervalle de la variable Tp
o $M$z: Nombre total d'éléments de la série
43
Collecte et analyse statistique des données

Analyse fréquentielle des données de la houle

1. Code sur Excel :


1. Fréquence d’occurrence Direction/Hs/Tp :
o B:B : Colonne de la direction de la houle.
Nombre total
o C:C : Colonne de la hauteur significative . d’observation

o D:D :Colonne de la période de pic.

Intervalle direction Intervalle Hauteur significative Intervalle période de pic

44
Collecte et analyse statistique des données

Analyse fréquentielle des données de la houle

45
Collecte et analyse statistique des données

Analyse fréquentielle des données de la houle

Le principal paramètre de l’état de mer est la hauteur


significative qui nous donne des renseignements sur
l’énergie des vagues et sa distribution au large, outre la Hs
la distribution de la période de pic nous renseigne sur le
système de vague qui prédomine. La direction est aussi un
paramètre essentiel qui nous donne les renseignements sur
la provenance des vagues.

46
Collecte et analyse statistique des données

Analyse fréquentielle des données de la houle

47
Collecte et analyse statistique des données

TP: Analyse fréquentiels (TP1 analyse fréquentielle de la houle)

 Calculer le pourcentage d’occurrence des houles provenant de la direction 40-50 degrés par rapport au
Nord avec (40 et 50 incluses).

 Calculer le pourcentage d’occurrence des houles provenant de la direction 100-180 degrés ayant des Hs
entre 0 et 5m (intervalles incluses).

 Calculer le pourcentage d’occurrence des houles provenant de la direction 0-80 degrés ayant des Tp entre
6 et 12s (intervalles incluses).

 Calculer le pourcentage d’occurrence des houles provenant de la direction 5-40 degrés ayant des Hs entre
0 et 5m et des période entre 6s et 16s (intervalles incluses).

48

Vous aimerez peut-être aussi