Homogénéisation Des Données

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

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

MINISTERE DE LENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE SCIENTIFIQUE

ECOLE NATIONALE SUPERIEURE DE LHYDRAULIQUE

HYDROLOGIE

TRAVAIL DIRIGE N3

THEME :

Homognisation des donnes

Travail tabli par :


SAIDOUN HAKIM

Anne universitaire 2000/2001


Introduction :
En hydrologie les donnes pluviomtriques constituent une base primordiale dans lanalyse et
llaboration des formules qui caractriseront la rgion et sa pluviomtrie, mais faudrait que
ces donnes soient fiables ? Les donnes prisent dans les stations pluviomtriques peuvent
tre fausses par plusieurs types derreurs (report de donnes, lecture, appareillage, utilisation
des eaux de pluie de la station ) do la ncessit de correction de ces derniers.

But : la correction et lextension des donnes pluviomtriques des deux sries appartenant
deux stations voisines. (voir le sujet ci- dessus )

I- Calcul du coefficient de corrlation et la dtection de lerreur


systmatique et la correction des donnes :
1- Calcul du coefficient de corrlation :
Il donne le degr de liaison entre les deux sries, il est donn par la relation suivante :

r = (Xi X) (Yi Y) / (Xi X)2 * (Yi Y)2


les variables sont comme suit :

X=
12 ( Xi ) / K = 412.66 ; 12Y = ( Yi ) / K = 312.5 ; K = 12
Xi Yi Xi - X Yi Y (Xi X) (Yi Y) (Xi X)2 (Yi Y)2

376 332 -36.66 19.5 -714.99 1344.44 380.25

553 446 140.33 133.5 18734.49 19693.44 17822.25

563 502 150.33 189.5 28488.16 22600.10 35910.25

330 300 -82.66 -12.5 1033.33 6833.77 156.25

477 437 64.33 124.5 8009.49 4138.78 15500.25

482 454 69.33 141.5 9810.66 4807.11 20022.25

359 280 -53.66 -32.5 1744.16 2880.10 1056.25

326 179 -86.66 -133.5 11569.99 7511.10 17822.25

382 190 -30.66 -122.5 3756.65 940.44 15006.25

369 230 -43.66 -82.5 3602.49 1906.77 6806.25

454 240 41.33 72.5 -2996.66 1708.44 5256.25

281 160 -131.66 -152.5 20079.15 17336.09 23256.25

Xi Yi (Xi X) (Yi Y) (Xi X)2 (Yi Y)2

4952 3750 103116.94 91700.58 158995


Tableau 1 : les diffrentes sommations pour la dterminationDu coefficient de corrlation
r = (Xi X) (Yi Y) / (Xi X)2 * (Yi Y)2
do
r = 0.85

2- Correction des donnes :


On utilise la mthode des doubles masses des deux sries hydrologiques de stations
pluviomtriques voisines (mme condition climatique).
La station de OUED EL ABTAL a fonctionner dix (10) ans de plus que la station de SIDI
AEK, do la premire cit sera la station de rfrence ou de base quon supposera juste et
lautre sera la station a tudier.
On prend une priode connue qui correspond a la srie de la station de SIDI AEK qui a
fonctionner pendant douze (12) ans cest a dire de 1969/1970 a 1980/1981.
On tabli le graphe des prcipitations cumul de la mthode des doubles masses afin de
dtecter lanomalie qui se traduit par une cassure de la droite. (voir figure 1)

Xi Yi Xi Yi

376 332 376 332

553 446 929 778

563 502 1492 1280

330 300 1822 1580

477 437 2299 2017

482 454 2781 2471

359 280 3140 2751

326 179 3466 2930

382 190 3848 3120

369 230 4217 3350

454 240 4671 3590

281 160 4952 3750

Tableau 2 : cumule des prcipitations

Voir figure 1 : on distingue la cassure de la droite des prcipitations cumules entre les annes
1974/1975 - 1976/1977.

II- correction des prcipitations annuelles :


Suite aux calculs des pentes de droites qui correspondent aux droites avant et aprs la cassure.
La pente m2 correspond la droite avant la cassure et la pente m1 correspond la droite aprs
la cassure, le rapport (m2 / m1 ) sera le rapport de correction des prcipitations de la station
SIDI AEK aprs la cassure (facteur de multiplication).
m1 = 2.25 / 4.0 = 0.563 m2 = 3.8 / 4.3 = 0.884
m2 / m1 = 0.884 / 0.563
m2 / m1 = 1.57

Station OUED EL Station


Anne
ABTAL SIDI AEK
69/70 376 332

70/71 553 446

71/72 563 502

72/73 330 300

73/74 477 437

74/75 482 454

75/76 359 280

76/77 326 281

77/78 382 298

78/79 369 361

79/80 454 377

80/81 281 251

Tableau 3 : prcipitations annuelles corriges.

1- Calcul du coefficient de corrlation aprs correction :

r = (Xi X) (Yi Y) / (Xi X)2 * (Yi Y)2


les variables sont comme suit :

12X= ( Xi ) / K = 412.66 ; 12Y = ( Yi ) / K = 359.92 ; K = 12

Xi Yi Xi - X Yi Y (Xi X) (Yi Y) (Xi X)2 (Yi Y)2

376 332 -36.66 -27.92 1023.55 1344.44 779.53


553 446 140.33 86.08 12079.61 19693.44 7409.77

563 502 150.33 142.08 21358.89 22600.10 20186.73

330 300 -82.66 -59.92 4952.99 6833.77 3590.41

477 437 64.33 77.08 4958.56 4138.78 5941.33

482 454 69.33 94.08 6522.57 4807.11 8851.05

359 280 -53.66 -79.92 4288.51 2880.10 6387.21

326 281 -86.66 -78.92 6839.21 7511.10 6228.37

382 298 -30.66 -61.92 1898.47 940.44 3834.09

369 361 -43.66 1.08 -47.15 1906.77 1.17

454 377 41.33 17.08 705.92 1708.44 291.73

281 251 -131.66 -108.92 14340.41 17336.09 11863.57

Xi Yi (Xi X) (Yi Y) (Xi X)2 (Yi Y)2

4952 4319 91700.58


Tableau 4 : les diffrentes sommations pour la dtermination
Du coefficient de corrlation aprs correction.

r = (Xi X) (Yi Y) / (Xi X)2 * (Yi Y)2


do
r = 0.94
Constatation :
Le coefficient de corrlation avant correction est estim r = 0.85 cest dire que le
coefficient de dtermination R = r2 = 0.72, alors on a 72% des donnes sont justifier, aprs
correction et par la mthode des doubles masses, on atteint un coefficient de corrlation
estim r = 0.94 soit un coefficient de dtermination R = r2 = 0.88, alors on a 88% des
donnes sont justifier, do une amlioration nette des donnes vrifier.
La mthode des doubles masses permet daugmenter la fiabilit des donnes.

III- vrification de lhomognit de la srie corrige :


le test de la mdiane permet de vrifier lhomognit de la srie par les lois suivantes :
Ns > (N+1-U1- /2 N+1 )
Ts < 3.3 [log10(N) + 1]
N = 12
On range la srie corrige par ordre croissant :
251 ; 280 ; 281 ; 298 ; 300 ; 332 ; 361 ; 377 ; 437 ; 446 ; 454 ; 502.
la mdiane Y50% = (332+361)/ 2 = 346.5
pour la variante on a :
Ns = 6
Ts = 6
On calcul les deux quations du test de la mdiane :
(N+1-U1- /2 N+1 ) = 2.96
3.3 [log10(N) + 1] = 6.86
do :
Ns > (N+1-U1- /2 N+1 )
Ts < 3.3 [log10(N) + 1] vrifi
On a alors une srie de donnes homogne.

IV- Extension de la srie :


On tabli le diagramme de dispersion (voir figure 2), on voix que le nuage des points est serr
do lextension des donnes de la station de SIDI AEK est possible.
On utilise la mthode des moindres carrs afin de dterminer la droite de rgression
Y=aX+b

12X= ( Xi ) / K = 412.66 ; 12Y = ( Yi ) / K = 359.92 ; K = 12


Xi Yi X2i Xi * Yi

376 332 141376 124832

553 446 305809 246638

563 502 316969 282626

330 300 108900 99000

477 437 227529 208449

482 454 232324 218828

359 280 128881 100520

326 281 106276 91606

382 298 145924 113836

369 361 131769 133209

454 377 206116 171158

281 251 78961 70531

Xi Yi X2i XiYi

4952 4319 2130834 1861233

la droite de rgression Y = a X + b tel que :


a = [ (XI Yi ) - 12Y XI ] / [ (X2i) 12X XI ] = 0.9
b = 12Y a 12X = -11.47
Do la droite de rgression : Y = (0.9) X (11.47)
on fait lextension des donnes de la station SIDI AEK :
Anne Station OUED EL ABTAL Station SIDI AEK
61/62 266 228
62/63 399 348
63/64 457 400
64/65 212 179
65/66 371 322
66/67 344 298
67/68 472 413
68/69 483 423
69/70 376 332
70/71 553 446
71/72 563 502
72/73 330 300
73/74 477 437
74/75 482 454
75/76 359 280
76/77 326 281
77/78 382 298
78/79 369 361
79/80 454 377
80/81 281 251
81/82 432 377
82/83 446 390
Somme :8834 Somme : 7697
Moyenne :401.55 Moyenne :349.86

V- Moyenne, variance et coefficient de correlation :


1- La moyenne :

Vous aimerez peut-être aussi