Home

Tarage et calcul des débits des stations hydrométriques non

Contents

1. 4 1441 1409 141 5140 140 JOAN 365 6070142 181 JOAN 366 142 IBIS J0AN 366 TRADUCTION EN DEBITS DES 16 HAUTEURS CONTENUES DANS UNE CARTE TRI DU DEBET MAXIMAL QMAX 2 IF CHANG 109 110 109 109 READ ge 6 ICHA IGRAD IFCICHA 1109 1108 1109 1109 JHIN 1 Ans ICHA IF IGRAD LL10 1111 1110 1110 IF IGRAD IGHA 1112 111201114 1112 IGRAD 60 TO 1113 1111 CALL TRAD CALL TARAG JIN ICHA 1 16 CALL TRAN TO 148 1113 1114 1108 1115 160 155 4155 159 5155 5159 2002 239 240 243 244 249 1190 161 170 162 163 164 165 175 167 166 168 169 180 190 251 252 250 212 271 219 332 CALL TRAD JMIN IGRAD L JHAX ICHA KK KK 1 1111 CALL TRAD CALL JHIN 1 JHAX IGRAD 50 TO 1115 JNIN 1 JHAX IGRAD CALL TRAD KK JHIN IGRAD 1 JHAX 16 CALL TRAD GO TO 148 JMIN 1 JMAX 16 CALLTRAD 50T0 2001 2002 2002 2002 2002 200242002 200212003 CAS TRAITEMENT DU CAS 1 2 1604149 149 9 607011504151 150 152 150 152 150 150 152 150 1523150 215 JMAX 15 5010155 1F11815 153 1544153 JHAX 12 5010155 JMAX 13 5210155 JHAX 14 15070155 16 IFCIUNIV 3 41555155r5155S 0015971 y JMAX 4Q 16 GROUP 1 J CS JO HOLSIZ0E6 Ji GO TO 101 005159 1 1 17 17 24 HGLITI HAUT LA CALL CORDES JO 16 GROUP 1I LA 03610 615
2. 6010101 YRAITEMENT DES CAS 2 8 ICA CAS 1 GOT0 239 240 240 2 41 242 243 244 ICA JHA 15 5078249 JHA 13 5070249 11 5070249 JHA 7 5070249 JHA 8 65070249 J4A 9 IFUIUNIV 3 190 1190 1190 CALL GRAD Ga 161 161 161 162 161 167 1611 ICA DO 170 LA 1 JX IT TT 24 CAS HG LTI HAUTILA CALL CORDES CONTINUE cu 190 DO 175 L 1 JX IF HDD LA 511163 164 163 1 17 25 60 165 IT 11 24 CAS HGIITIZHAUTILA CALL CORDEB CONTINUE 60 10 190 BO 190 LA 1 JX 1F MOO LA 2 1166 168 166 IT T424 CAS 60 169 IT IT 24 ICA HG ITI HAUTILA CALL CORDEB CONTINUE 00 250 J 1 JHA CAS JtICAJ CAS H 0 D 0 HH 0 D0252N 1 CAS K H 1 1 3 251 252 252 1 IF DEB L EQ 20 HMsMM 1 DEB LI 0 D D OEB L IF H LE O AND HH GT 0 QS0 K01 20 Hai 0 010001 035 10 IFIM GT 0 QJO K0 D H CONTINUE 1 5 73271 272 272 CAS 6 ICA CAS 1 50701279 280 281 282 283 1 291 292 291 290 291 250 291 291 290 291 290 291 4015 1F GRGUP 4 58 51 57 IF GROUP 4 58 56 56 IFIGROUP 41584293 293 1815 54 55 54 50701301 302 301 55 301 55 301 301 55 301 55 301 rHOIS IFIGROUP 7 55 51 51 IF GROUP 6 55 311 311 5 53 54 53 6010 321 322 321 320 321 320 321 321 320 321 3204321 815 IF GRQUP 8 54 53 53 1 46900 8 54 51 a GROUP 8 54 5 13 1 332 531 3373314 534331 331 53 331 53 331 8015
3. 100 200 300 200 750 Fig 18 M thode simplifi e de l cart la cote maximale station de Mopti TABLEAU VII STATION DE MOPTI M thode simplifi e de l cart la cote maximale Sens Ann e N H Q Q 29 Q AQ cm m3 s m3 s m s cm cm Wa m s m3 s D 1951 52 3 26 1 52 527 1570 1820 250 137 712 185 160 23 40 1860 A 27 552 232 280 325 45 139 480 160 21 5 330 C 195253 5 8 6 52 126 130 19 171 719 43983 60 11 1 12 D 1953 54 6 19 2 54 310 505 595 90 151 731 421 160 09 5 600 D 1955 56 7 23 3 56 230 290 20 30 94 731 501 160 66 20 340 D 8 195657 9 18 12 56 542 1770 1945 175 90 686 4 160 70 135 2080 10 18 2 57 233 265 325 60 185 453 160 25 10 315 11 8 5 57 100 74 74 0 0 586 160 160 12 86 C 12 6 8 57 473 1600 1425 175 412 2 257 160 37 55 1370 D 1957 58 14 16 1 58 514 1570 1720 150 87 730 216 160 73 125 1845 gt 15 19 2 58 346 530 765 235 307 384 160 147 10 655 16 12 3 58 268 280 430 150 349 462 160 189 80 350 17 46 58 227 360 30 50 162 480 160 02 0 gt 1958 59 18 25 7 58 414 1225 1095 135 123 707 4293 4160 37 40 1055 19 13 8 58 464 1555 1370 185 4135 243 160 25 35 1335 20 6 9 58 580 2400 2250
4. 329 VARIANTE DU 310 POUR LE TRAITEMENT DE BASE DES DONNEES LIMNIGRAPHIQUES INTEGRALES DE STATIONS NON UNIVOQUES CALCUL DU DEBIT PAR LA METHODE OU GRADIENT LIMNIGRAPHIQUE IUNI V 3 INTEGER STATI INTEGER HEU INTEGER 2 JOUR16 MANTU6 9 ICAR 5 0 70001 y HANJOU368 102 70001 IC 13401368 HALT 2017 12 32 03412 31 12 CHOI CL2 DIMENSION MI 16 7 10116 DIMENSION Q 7000 HEURE 7009 e QJ0 366 DEBHOI 12 eFRAG 18 DIHENSION NODEBA 12 COMMON JHINs STATI ISTAT LHAX HP 16 KHAXy 31 25 3 151 146 61 QMAX 11 61 800 s JGRAD r LUNI Ve MZ y 182 HEU 6 HENG 2Y NHB 5 5 5 BECK NODESA 1 30 58 89 119 150 180 211 242 272 303 333 1 10 16 10 15 0 FORMAT 212 14 18 4 FORMAT 212 214 120 64312 E Ze 21119 FORMAT 1X ERREUR DE STATION POUR 16 FORMAT Lie ERREUR POUR LA STATION 13 12 14 15 FORMAT 12 FORHAT 1Xs ERREUR DE SEQUENCE DES CARTES A LA CARTE 416 FORMAT LX ERREUR DANS LE COMPTE DES JOURS POUR 216 1212 214412461312 13 111411 151 FORMAT 212 144212016113 11 14 FORMAT 1 NUMERO DE 57 4 212 18 4 17 4 0 5 1 JOURNALIERS EN H3 S 212 14 212 11 16113 11 13 FORMAT 212 14 13 13413 11 16X 112 2 2 14 212 11 15 13 11 17 FORMAT 1 2
5. 015 1 127 156 127 NODESU O 160 IFCIRISIAN481 80 NODEBUZNODEBA MAIS JOUR 1 6070 169 NODESU NODEBA NOT S 0 int KENNDERU DKCH IFCJUHTYFJ 162416421535 00163 1 1 6 HAUTA II HAUTA T 0 1 NOC NDCAR JFINNC 50001164 726 776 166 1 1 8 1 5 1 IFLICHANGH 1661709 FES REAN S 5 IC PA JMAX JCHA TF NOCAR 5900 400 417 610 CALL TRADA 410 401 170 500 19 180 651 699 702 700 698 190 195 186 TD 471 CALL IF UNOCAR 50091890 819 10 CALL 70 189 CALL pang no 698 J 1s1FIN JF1JUDERI 51 655 651 NEB J DER 3 1000 G IF X 999 51700 702 702 699 2005 HANTIJ3 10008x ICAR J 1 IF FIN 6 190 196 196 JFIN 1F1N 1 00 195 J JFIM 6 JOR JO 18123 1 MANT J 10 1CAR JI 0 HRITE 7 RIKETA KBASy IBYREP JANNe MOT Ss JOURE YY e HELL A ITS MANTES 1 1CAP 3 1 6 JUOEBI NOCART 291 200 1000 1310 1320 1300 1001 1005 1010 1011 1012 1013 1025 1020 1022 1030 1044 1050 1060 1079 1654 1660 1599 1702 1709 1670 1109 3700 3010 3050 3070 3031 3082 3053 3100 MRITE 6 20 MOI So JOURES JI DEBLJI e EPIN 6 2 LHC ART 00200 1 6 Ha 3 200 200 201 JDTKI JOURT I HEURE K 6
6. 272 13 211 2 7 2 23 1 212 2184 12 61312 4 21 211 14 FORMAT IXe ERREUR STATION ONUR T6 FORMAT I1Y ERREUR DUANNEF POUR LA STATION 8 13 12 15 4 4 15 FORMAT T2 FORMAT LX ERREUR DE SEQUENCE DES CARTES LA CARTE NO J6 FORHAT IXY ERREUR DANS LE COMPTE DES JOURS POUR 16 212 214 12 64312 13 11 11 15 212 14 212 16013 11 14 FORMAT 1 NUMERO DE STATION 4 212 14 Xr ANNEE tr 17 H 1 JOURNALIERS FN M3 51 212 14 212 11 16113 113 131 FORMAT 212 14 139 13113 111 16X 11 FORMAT 212 14 1044 FOPMAT 212 214412 61312 14 211 14 FORMAT 212 149212 11 14113 11 17 FORMAT 1 16 14 614 14 13 27 0 7120 17 FDRMAT 1 NUMERO STATION 4X 212 1494Xo ANNEE 01704 Xe DERITS 1 JOURNALIERS EN LITRES SECOHDE 7012 99 650 620 625 630 635 670 5552 100 19 55 7001 51 7002 60 112 8 10113 11 4 18 10 13 11 4 54122 5415 5 5 8 IF JETAJ650 9999 650 1 8 11620 620 630 6 625 BRE FAORMAT 11 7A2 7A2 TRASSIN REPRESENTATIF 23A1 60 TD 670 MRITF 4 635 FTA RASHYDyBRE 1 742 742 BASSIN EXPERIMENTAL 2341 READIS 13 IRVREP NOMBY WRITE 6 5552 IRVREP NOMBV FORMAT1
7. 17 5 18 17 5 0 0 1 510 2 540 2 595 55 2 1 3 2 9 0 8 20 2615 1815 1 980 165 8 3 11 10 7 2 4 45 2 025 120 123 3 2 4 3 2 9 0 5 1 124 1 100 905 195 21 5 28 27 3 5 8 55 850 1 445 1 250 195 15 6 20 19 5 3 9 50 1 200 1735 1 460 275 18 8 11 10 7 8 1 120 1 580 2 020 1 875 145 7 7 6 5 9 1 8 35 1910 2 540 2 625 85 32 A 3 9 7 1 185 2 440 2 645 2 720 75 2 8 3 2 9 0 1 3 2 725 2 400 2415 15 0 6 5 4 9 5 5 135 2 280 2 545 2 640 95 3 6 7 6 8 10 4 275 2 365 725 425 300 70 5 22 21 4 49 1 210 635 465 385 80 20 8 23 22 4 1 6 5 380 1 045 850 195 23 0 17 16 6 6 4 55 905 1 350 1175 175 14 9 13 12 7 2 2 25 1 200 1 665 1 470 195 13 3 11 10 7 2 6 40 1 510 2 005 1 895 110 5 8 8 7 8 2 0 40 1 855 2 205 2 030 175 8 6 2 1 9 6 7 135 2165 1 130 1 460 330 22 6 13 12 7 9 9 145 1315 1 905 1 500 405 27 0 21 20 4 6 6 100 I 600 515 670 155 23 2 32 31 2 8 0 55 725 1 060 870 190 21 8 25 24 4 2 6 25 845 2 120 1 945 175 9 0 4 3 9 5 1 100 2 045 2 260 2150 110 51 9 8 8 3 1 80 2 070 2 290 2 180 110 5 0 6 5 9 0 9 20 2 160 2 330 2 260 70 3 1 3 2 9 0 2 5 2 265 2110 2 220
8. 110 4 9 5 4 9 0 0 2 220 1 980 2 120 140 6 6 6 8 02 5 2 125 1 840 1 930 90 4 7 2 1 9 2 8 5 1925 1710 1865 155 8 3 4 3 9 4 4 10 1855 1 640 1760 120 6 8 7 6 8 0 0 1 760 Ah sur 5 jours Ah sur 6 heures exprim s en centim tres de hauteurs d eau sont Ah les valeurs brutes Ah pendant les intervalles At et non pas les gradients propre ment dits la division tant une op ration suppl mentaire inutile Le tableau montre l efficacit de la m thode en pr sentant les diverses r ductions d carts absolus et relatifs et les pr cisions sur la connaissance du d bit obtenues pour Dir et Mopti Les r ductions d carts sur AQ et sur AQ Q qui sont de l ordre de 50 55 sont galement de 50 pour Boitiek On se souvient que cette r duction pour AQ tait de 75 sur Boitiek en utilisant la m thode it rative La prise en compte de la vitesse U am liore le gain mais on peut consid rer que cette am lioration est com pens e ici par la simplicit plus grande de la m thode graphique du gradient limni graphique est peu pr s certain que pour des cours d eau pente notable et fortes varia tions de pente la vitesse varie galement de mani re non n gligeable et qu elle doive tre prise en compte si la station ne peut tre dot e de 2 chelles 23 24 H N Date cm 1 7 7 65 139 2 29 7 65 275 3 31 7 65 209
9. 16 6 0 24 3 128 21 31 3 60 120 84 96 12 12 5 5 1 9 10 6 10 86 28 22 4 60 92 55 55 0 0 3 1 1 11 1 56 29 11 5 60 74 39 39 0 0 2 0 7 0 7 0 39 30 23 17 60 376 1075 930 145 15 6 44 16 4 0 8 5 925 D 31 9 6 63 140 130 128 2 15 6 3 1 1 16 7 21 149 C 32 24 5 67 52 44 33 11 33 3 1 0 4 33 7 11 44 gt 33 15 10 67 690 3 370 3 285 85 2 6 6 22 0 4 15 3 300 34 17 10 67 693 3 340 3 315 25 0 8 5 1 9 0 9 30 3 285 35 19 10 67 696 3 450 3 350 100 3 0 7 2 6 0 4 15 3 365 36 24 10 67 702 3 465 3 410 55 L 7 2 6 1 0 35 3375 37 29 10 67 708 3 605 3 480 125 3 6 5 1 9 1 7 60 3 540 38 1 11 67 710 3 565 3 500 65 1 9 5 1 9 0 0 3 500 39 4 11 67 713 3 540 3 535 5 0 1 6 22 2 1 75 3 460 40 7 11 67 716 3 585 3570 15 0 4 4 1 5 0 9 30 3 540 41 11 11 67 719 3 615 3 610 5 0 1 4 1 5 14 50 3 560 C 42 13 11 67 720 3 640 3 620 20 0 5 3 1 1 0 6 20 3 600 D 43 16 11 67 720 3 630 3 620 10 0 3 0 0 0 3 10 3 630 44 19 11 67 718 3 610 3 595 15 0 4 3 11 1 5 55 3 650 45 22 11 67 715 3 440 3 560 120 3 4 6 2 2 1 2 45 3515 46 25 11 67 711 3 460 3510 50 1 6 2 2 0 8 30 3 540 47 28 11 67 704 3 390 3435 45 1 3 7 2 6 1 3 45 3 480
10. 40 12 16 14 02 5 3435 gt 48 30 11 67 701 3330 340 70 21 19 13 0 4 15 3385 50 22 10 68 60 2840 2875 35 12 01 13 35 2840 e 3 51 25 10 68 651 2790 2890 100 35 651 0 0 35 100 2790 53 1 11 68 650 2900 2875 25 09 1 01 LO 30 2905 55 8 11 68 641 2760 2780 20 07 10 09 02 5 2785 1968 60 57 14 11 68 629 2610 2670 60 22 22 20 02 5 2665 59 18 11 68 619 2500 2575 75 29 32 29 0 0 2575 61 22 11 68 608 2390 2480 90 36 43 39 03 10 2490 62 25 11 68 602 2310 200 10 45 49 44 01 5 2415 TABLEAU STATION DE DIRE M thode de l cart la cote maximale AQ Ann e N Date Q Oo AQ Go Hx Q 2 m3 s m3 s m3 s Ka cm cm m s m s 1951 52 1 26 8 51 375 1 060 990 70 7 1 601 226 160 830 2 4 4 52 465 1 245 1 430 185 12 9 136 110 1 540 3 11 1 54 601 2 540 2 580 40 1 5 606 5 0 2 580 1953 54 4 23 2 54 542 1815 1 915 100 5 2 64 190 2 105 5 10 6 54 129 120 155 35 22 6 477 5 160 6 9 8 54 353 1 100 890 210 23 6 613 260 5 885 7 24 8 54 419 1 445 1 190 255 21 5 194 15 1175 1954 55 8 5 9 54 456 1735 1 385 350 25 3 157 50 1 435 9 12 10 54 527 2 020 1 790 230 12 9 86 80 1710 10 15 12 54
11. 47 28 11 67 704 3390 3430 40 12 16 1 5 0 3 10 3440 48 30 11 67 701 3330 3400 70 2 1 19 17 04 15 3385 C 50 22 10 68 650 280 285 35 12 1 01 13 40 2835 51 25 10 68 651 2790 2890 100 3 5 651 0 0 35 100 2790 D 53 1 11 68 650 2900 2875 25 0 9 1 0 1 10 30 2905 55 8 11 68 641 2760 2780 20 0 7 10 09 0 2 5 2785 1968 69 57 14 11 68 629 2610 2670 60 22 22 2 0 0 2 5 2665 59 18 11 68 619 2500 2575 75 2 9 32 2 8 0 41 0 2575 61 22 11 68 608 2390 2480 90 3 6 43 37 01 0 2480 62 25 11 68 602 2310 2420 110 4 5 49 42 03 10 2410 s d 42 E 5 4 8 43 5 3 o 3600 38 KE 3000 de hautes eaux 2 3500 36 46 45 7 2000 3300 Hauteur en cm 690 700 710 720 100 Hauteur en cm ee ur en cm 50 100 150 200 250 700 Fig 17 M thode de l cart la cote maximale station de Mopti Jaugeage de crue Point rectifi 3000 rage de eaux Jaugeage de d crue 2 Voir correction Fig 19 2000 400 cm 4 ob 19 4 e 3 u 300 i SS de basses eaux 1000 200 100 Gr Hauteur enem o 150 200 250 Hauteur
12. Fig 27 Sortie de courbe de correction de pente crique Virgile 09 00 00 00 01 01 01 le 312 le seul dont nous ayons eu besoin jusqu ici Ce programme d rive du programme 310 de calcul des d bits limnigraphiques int graux Nous nous contentons donc ici de pr senter l organigramme g n ral de ce programme 312 en insistant seulement sur ses particularit et en renvoyant le lecteur l article descriptif du 310 pour les blocs d instructions communs 3 5 3 1 PR PARATION DES DONN ES D ENTR E Les donn es d entr e sont a Les relev s limnigraphiques int graux relatifs aux deux chelles b Les informations relatives aux talonnages Les relev s limnigraphiques sont introduits sur les cartes 311 ou 511 qu il s agisse de ceux de l chelle principale HAUTA ou de l chelle secondaire HAUTB destin s au calcul de la correction de non univocit On ne proc de l tablissement des relev s HAUTB que si la station est en situation de non univocit ce qui pour certains types de cours d eau ne se produit pas toujours Ces donn es sont pr sent es sur des feuilles de relev s limnigraphiques int graux de contenu identique celui des cartes COH 311 ou 511 afin de faciliter la perforation Pour concentrer au maximum l information sur les cartes et contrairement ce qui avait t dit dans 3 page 103 figure 4 on doit remplir les six
13. 207 809 510 820 830 PERFCRATION DU VECTEUR 0 INTERPOLE DO 600 N 1 22 AC N 00 550 1 17 Ha17 N 1 J 355 1 15 260 561 560 MAX 366 IF M HAX 1562 562 563 0 6010564 X Q M 1 0 TF LX Zelt I I 1 6070509 X 0005 X 21000 HANT J DE ICAR 13 1 MRITEL 5 5 JETA IBAS ISTAT MILL LMANTELO y CARELO 101 173 WRITES 6 1500 FRAG JETAYIBAS ISTAT MILL WRITE 6 1501 0081030 1 31 00801901521 12 VARIMOLSI OJ JO MOIS 90 291802 805 807 IF 420 1 5 5 0 1 803 804 803 6 1511 0 VAR 6010810 WRITE 6r 1512 J0 VAR 5010810 1815 806 803 806 MRITE 6 1513 J U 11 VARCMOIS M01523 12 5070810 IF 1JO0 30 806 806 809 WRITE 1514 JO VAR 1I VARC3 e VAR 5 DE VARAT 9 VAR 8 VARTIOI 121 CONTINUE 0002401 5 1 12 DEBEH HQ S 0 001 FLOAT MAQHO MOIS 10 ICQHO HQ S WRITES 6 1520 0 HOAN 0 008308015 lt 1 12 015 HOAN HOAN 12e 0 001 10 1 WMRITEL6 1550 MDAN QMAX 9995510 520 520 102 IF IUNIV 3 3111 3104 3104 3104 3101 3111 999 1001 1002 1003 1004 1012 1013 Nm 1011 9 1010 20 QHAX 0 READ 5S LI ETAT BASSE STATI ANNEE HOIS GROUP y CAS CHANG IFLETATI103 99 103 0 5 8 STATI y HAAX2 IFISTATI 3103 99 3103 STOP SUBRO
14. quelques grands fleuves bien particuliers la m thode du gradient limnigraphique est g n rale pour tous les autres cours d eau non univoques une chelle mais sa mise en uvre est plus difficile Il faut en effet disposer de la chronique des hauteurs d eau ant rieures la hau teur HAUT LA traduire cette chronique tant dat e par rapport une origine arbitraire des temps laquelle est galement rattach e la date de HAUT LA Cette chronique dat e est appel e le vecteur hauteurs temps HG IT La correction de non univocit requiert le calcul du gradient limnigraphique DH diff rence de cotes entre la hauteur HAUT LA et celle observ e JGRAD aupara vant JGRAD tant le pas de temps du gradient limnigraphique Dans les exemples tudi s cf 4 1 la correction de d bit est proportionnelle au gradient limnigraphique Tout ceci exige un vecteur HG pas de temps constant et suffisamment petit pour que le calcul de DH soit pr cis Le pas de temps du gradient limnigraphique JGRAD conditionne le choix du pas de temps de IT Il cro t avec l importance du cours d eau passant dans les exem ples choisis de 6 heures Boitiek 250 km 5 jours Mopti Dir 2 300 000 km En pratique l exemple de Boitiek est justiciable des relev s limnigraphiques int graux On peut donc consid rer que pour un cours d eau dont les hauteurs n ont pas tre prises plus de 8 fois par jour le pa
15. 150 67 127 1144 47 105 2145 C 22 22 8 59 515 2000 1725 275 15 9 175 160 01 0 175 D 26 14 3 60 142 120 145 25 172 690 548 16 0 12 2 143 1959 60 27 313 60 120 84 101 17 168 570 160 08 1 100 28 22460 92 55 64 9 141 598 160 L9 1 65 29 11 5 60 74 39 47 8 170 616 160 10 Z 0 47 196061 30 23 760 376 1075 90 416 181 691 315 160 21 20 930 D 1962 63 31 9 6 63 140 130 140 10 71 702 52 160 89 2 152 32 24567 52 44 38 158 676 160 02 0 38 gt 33 15 10 67 690 3370 3285 85 26 720 30 27 01 5 3280 34 17 10 67 63 3340 3315 25 08 27 24 16 55 3260 35 19 10 67 696 3450 3350 100 30 24 22 08 25 3375 36 24 10 67 702 3465 3410 5 16 18 16 0 0 3410 37 29 10 67 708 3605 3470 135 39 12 Li 28 95 3565 38 1 11 67 710 3565 3495 70 20 10 09 L 40 3535 39 41167 713 3540 3520 20 06 7 06 0 o 3520 1967 68 40 7 11 67 716 3585 3560 25 07 4 04 403 10 3570 gt 41 11 11 67 719 365 360 5 01 p 4 L Ol 0 0 3610 5 42 13 11 67 70 3640 3635 5 01 0 0 01 5 3640 D 43 16 11 67 720 3630 365 5 01 0 0 01 5 3630 44 19 11 67 718 3610 3595 15 04 2 02 06 20 3615 gt 45 22 11 67 715 3440 3545 105 29 5 05 24 85 3460 46 25 11 67 71 3460 350 40 11 GE 0 8 03 10 3490 47 28 11 67 704 3390 3430
16. 1X NUMERN DE STATION 4X 212 14 1 CALL TAPAG RBEADI 5 19 JUDEBI FNRMATIITLY 151 55 51 6170011 FORMATI T225 DEBITS INSTANTANFS EN 3 5 4 60 HRITELE 7002 FORMATI TITr DERITS INSTANTANES LITRES SECONDE QHAX Q REA 5 9 2NETA JRAS 5 JANN 4 015241 JOUR 1 HEUCTD MINTIS HAUTE STATION 1044 11 1 6 y JUNTTE y ICHANG y NOCAR 40 310 NDCA NIN AR IF NOCA 5000 308 310 310 RBFAD 5 2 NFIA JRAS N0STA JANNH s MOIS 1 I rHINET HAUTAUT 141 1 62 JUNTTE JCHARNG NOCART 125 130 150 127 80 81 169 161 162 163 726 tes 164 165 400 330 1 1 6 HAUTS 1 HAUTI 1 60 ra 320 340 1 1 6 HAUTA 1 HAUTE 1 NOCART NOCAR 111 1000111 DQ 6000 151 6 1 16001 6001 6000 1 1 CONTINUE TFIN TF IBVRFP HNSTA 1126115r112 HRITEL6 SA 9999 FORMATIT22 DERITS MOYENS JNURNALIERS 1 14 7136 3 5 722 10 8175 MOYENS JOURNALIERS er Le 1477734 1 4 525 1 1 13 1 FORMATIA 13 19 3 FORMAT 1X IF NOCART 1 116 120 116 NOCART ICARTF 1 117 118 117 HPITELG 6 COTN113 IFIMIIL JAHN1119 120 119 MRITEL6 4 y JIRAS NOSTA MILL 60710113 MILL JANN IRIS 1J4NN 41 4 JANN IF NDCART 11161 125 161 TF SOUR 1I 1 177 130 127
17. 6 000 6 100 0 100 1 6 0 0 1 6 0 100 6 000 7 800 8 050 0 250 3 1 5 3 7 0 6 0 045 8 095 0 110 0 112 0 002 1 8 2 1 5 0 3 0 0 112 0 0 4 2 1 DOMAINE D APPLICATION La m thode de tarage avec correction par le gradient limnigraphique para t s appliquer tous les types de stations non univoques x r Pour les grands fleuves r gime hydrologique r gulier hydrogramme annuel en forme de cloche d l importance du bassin alimentaire aux faibles pentes du chenal et la pr sence r gularisatrice de grandes plaines d inondations on a envisag Papplication d une m thode simplifi e par rapport celle du gradient En effet pour ces cours d eau la pente de 1 ligne d eau est peu variable en crue et galement en d crue ce qui permet d assimiler le tarage deux courbes pseudo univoques de crue et de d crue r unies par une courbe de raccordement variable chaque ann e selon la cote maximale atteinte Cette relative univocit de crue et de d crue s explique par le fait que le gradient limnigraphique varie peu ou tr s lentement tout au long de la crue comme de la d crue les limnigrammes annuels sont sensiblement homoth tiques cf fig 15 pour le Niger Dir Il semble pour de tels cours d eau que la valeur du gradient limnigraphique d pende de la cote maximale et de l cart cette cote d o l id e de proc der la correction de pente
18. 8 7 WRITEI6 5 IL CIK L K 1 3 L 1e EMAX MRIYE16 111 L HP L 1 10 HP LMAX HP 1 1 0 H 1 HP 1 DEB 1 CA 1 003001 2 1 HITI H I 11 0 1 HT H I IFIHT HP KMAX 170 110 110 L KMAX 6070150 DL Ze IFIFRT HP L 130 130 140 L L 1 GOTC150 CONTINUE a X HT HPIL DER I X C 1sL1 x C 2 L C 3 L IF49EB 11 2S9 309 200 NES11 0 CONTINUE WAJTEL6 8 43210 220 230 MRITEL 40241 TO 210 25 4 9 ng Ts1 10 NO 104J 1y6 1 1 10 1 8 1 11504104 104 CONTINUE 1 71 1000 7 524 1990 50 018 0 801 802 403 8 3 8933 2 6 20018 I 1 27 24 28 3 GE JO 500 RRITELA 200L FLE DERI 1 1 21 2 27 21 TO SEC 11 6 420021 1 4 1 fORVAT 1X F31 2 F26 1 560 9 I HNFBL I 1 527 2 25 01 TINUE 142 18 216141523167 JAL 10 1 AQ IG 160 167 KHP L 0 140 163 KHP 1 x 100 0 1 160 CONTIRUE 2 LIPAY LNAX L 00 146 J LIVAX 12 164 KHP IX 0 HRITE 75 14 2741 RASSI STATE 61 LI Le Le 16 IUNIV 60 19 502 100 LZILEAX y Lalo L MAX 294 oFANtSr L9 O0NP Ls KWAX LMAX 1 DEMAS 20 5 20 QQNIUT
19. 81 IF JUNIV 1UNIVI91 1010 91 91 HRITE 6 1013 23 1013 FORMAT 1X DE METHODE DE TARASE 1010 RETURN END 24 SUBROUTINE TRAD EN INTEGERSTATI COMMON JMIN 3 5 1 161 3 151 0 3 153 1 61 6 QHAX IT HG 800 JGRAD AC UNIV IFIN HZ 2 61 MINT E 2 8 ALPHA SI BETALSI HB 51 2 0013033 MIN JHAX HAUTES DO 115 2 1 831194119 112 110 6010120 112 IF HAU HPIKHAX 115 119 116 115 CONTINUE 116 99 991 118 117 117 117 0284 10 5070130 119 1F HA1 9 99 201 201 119 201 0 8 3 0 65070130 FORMAT II 0 1 12 641 1 1 WRITE 6 31 29 J 1 32 00 27 Ms1 12 JH 54 M 1 1F J 29 17413 15 TFLHOD JAN 41114 17 14 IF M 2 174 253417 2 31214 16 18 GD TO 17 25 17 25 17 25 17 17 25 17 25 171 1 4 25 20 20 LFCKHOMI325 19919 JO LM HI JQ KM M 10 JD GD TO 21 J JQ KJ H J 10 JD JA JK 1 11 IR MINO 2 16 303 DO 24 5 1F 1 2 90 239 22123 IP IP 1 SJACJH T 12 60 TO 24 19 30 10 JD 1 3 1P JAL JM 1 10 10 10 10 0 0 10 CONTINUE GO TO 27 00 26 1 1 5 JA JM 1 11 CONTINUE TF J 31 26 29 30 FFLMOD J 4 51 EQ O HRITE 6 1 HRITE 6 2 J AA JAlI lI sie tie 6 4 34 122 1 1 50 RETURN END SEP
20. CAS ICA SAC eCORQ IT GROUP HOIS 1815 1 JX DOL30J JHINy HAU HAUT LA IFCHAU LT 0 60 G0 TO 129 SUBROUTINE GRAD INTEGER CAS GROUP COMMON JMIN IMAX HAUT 16 HP 161 DEG116 A 3 15 yQMAX KMAXs T COMMONSTATI ISTAT C 3 15 HMAXOp HHAX 1 e HMAXZ 9 KK AA 9 AB e TUNIV COMHONEG 1300 3 CAS ICA JGRAD LT 9 GROUP HOIS IBIS LA 5 JX 60 TO 1 2 3 4 5 5 5 1 3 28 28 8 60 127 9 27 26 27 26 27 27 26 27 26 27 MOIS 00 115 2 1 22110 1108112 JF 1B15 24425 24 IF GROUP 6 25 10 21 L t 1 015 2325 11 25 120 11 IF 181S 23 24 23 IF HAU HP KHAX 115 119 116 IF GROUP 8 24 2441 CONTINUE 99 991 119 117 117 JF GROUP 10 2313921 DEBE 10 13 IF M0IS 2 23 14 23 6070130 14 1 1815021 22 21 L RHAX 5 IF GROUP 15 22915921 H HAU HP LA 15 IF MOLS 2 22 21 22 1 LoLISHYAU 2 1 Tei Set 21 Jx CAS TF IUNIV 3 90 130 130 60 70 29 GO TO 1118 1119 1120 1121 KK 22 JX CAS 2 DHMAX HAU HHAXO 50 TO 29 60 TO 1122 23 5 3 DHHAX 1 60 TO 29 60 TO 1122 24 1 5 4 DHMAX HAU HMAXL TO 29 60 70 1122 25 4 5 5 HKAX2 HAU GO TO 29 CORQ AA ATAN AB DHMAX 26 1 5 6 DEB J DEB J 1 C0RQ 100 60 T0 29 1 0 8 23 17 0 10 8220 27 JXeCAS 7 IF QUAX LToDEB J IQHAX DER 3 TO 29 50 TO 130 28
21. CORQZAC SQRT DH UI 22 CORQ AC DH OEB J DED J1 1 CORQI 1 170 0 IF QMAX LT DEB J1IQMAX DEB J 111 720 100 19 19 IF N 1 100 100 9 112 T1 480 113 1T2 1 11 0 40 J2x113 1T1 IT IT 1 1 0 17 250 11 12 11 Hz Hz 1HZ 24C8MZ 2 HGCITIER 1TN IT 1 00 30 Jl ITN ITI HG J1 9 99 CONTENUE RETURN END SUBROUTINE YHPRI JAN KJLI KH LN INTEGER 2 KJ 12431 1 90125313 KH 121 LH 12 JA160 DIMENSION AA 12 DATA 71 424 434 44 15 46 47 781 49 40 4 FORMAT 1X FORMAT 1 3712 1 5 1 14 776 JAN FEV MAR AVR HAL JUN JUL 1006 18 14 3 11 8 14 4 2 12 5 3 1 4 1 4 21 2 CT NOV DEC 1012 FORMAT 4X 14 2X 12 8E8 5 4 READ S 1001 STATIeLHAXS HP L 1 16 LUNIV IFISTATI ISTAT 1291 1 URITE 6 1003 2 KMAX LMAX 1 IFIKMAX 8 12 12 10 12 00 6 1 1 3 13 READ 5 1012 STATIsIv C J L LEL 14 IFASTATI ISTATID 403 15 3 WRITE 6 1003 16 4 DUGL 1 KMAX 17 6 18 1011 19 10 DO 9 41 3 gt 045 1002 5 1 1 LECSTATI ISTAT 07 8 7 20 HRITE ge 1003 8 00 9 L 1 KMAX 21 9 AllsL C J LU 1011 READ 5 1006 STATI JUNI Vy IGRADr AC INHB LALPUACIA BET A LI HBE I e Tele 2 57 71 15 271 81471 71 WRITEG6 1004 22 1004 FORMAT 1 Xs ERREUR DE STATION POUR LS GRADIENT
22. JOUR JI 0 HEU J 1 HIN J 1 33 10 ICAR JIZO MWRITE 7 8 JETA IBAS y ISTAT HOE Ss COUR 12 3 HEUT IY MINGI p 1 1 edel 9 6 201 200 KRITE 6 20 MOIS 9 JOUR 35 SD MINTSI DEB IFIN KRITEL6 2LINOCART 00200 J 1 6 IFCJOUR J3 1200 200 201 K K 1 JQ K JGURIJ HEURELK GO4HEUV AN Q K DEB J CONTINUE ICARTE NOCART 1000 1319 1320 1320 1001 1005 1010 1011 1012 1013 1025 1020 1022 1030 1049 1050 1060 1070 1650 1660 t 1699 1702 1709 1670 1100 3000 3010 3050 3050 3070 3080 3081 3082 3083 110 INI NODEBU 2 JOR INI 1 INI 1 DO 1300 IO INI K IF JDLU10 JO 10 11 1320 1310 1320 JOR 10 30R IO 1 TO 1300 JORE 1072 084 10 13 1 CONTINUE NRELEV JOR K IF HODEBU 1001 1010 1001 001005 J 1 NODEBU QJ0t3 10 181521011 1012 1011 H90 365 60702013 NJO 366 IFLNRELEV NJDI 1020 10309 1025 WAITEC TIMILL 2000 NREL NRELEV 1 00 1022 J NRE L NJD 920021 10 HeNODEBU 5 0 3 SmS HEURE M LI HEURE CM F4 QUM Q M9 139 20 IFLHEURE M 13 1440 1050 1060 1060 L JORCHI 03061 5 14 0 IFCL NRELEV 1969 1070 1070 0011003 1 90 X QJ0 J 1F XIL650 16609 1660 10E 10 I 0 60701670 Q 1 99951170091702 1702 X 0 1 X 1 1 1 0 01699 00
23. des hauteurs HAUT LA de la carte trait e b L appel CORDEB pour calculer le gradient limnigraphique et le d bit corrig Apr s l ex cution de cette instruction de calcul D il faut d terminer le nombre exact de relev s de hauteurs d eau r els contenus dans la carte lue En effet les 16 colonnes hauteurs de chaque carte sont remplies mais selon l ann e le mois la position dans le mois et le CAS certaines colonnes ne contiennent que des 9999 indiquant des relev s fictifs c est dire compl tant la carte au del du dernier relev r el du mois Or le vecteur hauteurs temps HG doit tre g n r uniquement sur des relev s r els afin que sa base de temps IT reste chronologiquement valable Ce travail est effectu dans le programme principal pour le CAS 1 valeur de JMAX et fait l objet d une sous routine GRAD pour les CAS 2 8 valeur de JX La constitution du vecteur HG hauteur temps intervient ensuite On fait l hypo th se que les observations effectu es dans une journ e sont r guli rement r parties au cours des 24 heures et par cons quent pour chaque observation le pas de temps est pris gal 24 CAS en heures C est dire qu l issue du vecteur cr en d but de programme l aide des relev s ant rieurs on ajoute les nouveaux relev s d une carte traiter un par un en augmentant la base de temps chaque fois de IT CAS Cette hypoth se de r partition journali
24. maximum valeur du compteur lt crue d crue C points valeur IGRAD 0 Si l ann e calendaire commence avant l tiage A la variable prend la valeur 1 si elle commence apr s 2 Le test s effectue sur la valeur du maximum de l ann e ant rieure N 1 appel e HMAXO On place donc les cartes de donn es dans l ordre suivant a Une carte d identification de station b Cinq cartes lues dans la sous routine TARAG et contenant l talonnage univoque pour les quatre premi res les param tres de l quation de correction de 51 52 non univocit et AB pour la cinqui me carte 513 format 2 F 5 2 pour les deux param tres Une carte des hauteurs maximales de trois ann es N 1 N 1 soit HMAX 0 HMAX 1 HMAX 2 suivant indication de la station STATI et de l ann e LAN format 14 12 3 F 42 d Les cartes COH 301 des hauteurs de l ann e N Par convention s 0 est nulle donc inutile l ann e N commence apr s et 2 0 non nulle KK 1 La variable IGRAD est lue dans les m mes conditions que la variable ICHA indicatrice du changement d talonnage Le processus est le m me une perforation d un chiffre quelconque non nul dans la colonne 80 de la carte 301 contenant le jour du passage l tiage ou au maximum annuel entra ne la lecture d une carte ICHA IGRAD en 2 D dans laquelle IGRAD
25. re r guli re n est pas rigoureuse mais 2 l erreur due cette approximation est n gligeable Comme l unit du pas de temps est on d consid rer part les CAS 5 et 7 puisque 24 heures ne sont pas divisibles ni par 5 ni par 7 d o les hypoth ses faites pour ces deux cas 55 56 a Pour le CAS 5 4 pas de temps de 5 heures obtenus en faisant 25 CAS et 1 pas de temps de 4 heures pour clore la journ e 5 Pour le CAS 7 alternativement 1 pas de 3 heures et 1 pas de 4 heures obtenu en faisant 241 On comprend dans ces conditions que ni le calcul de JX dans la sous routine GRAD la constitution du vecteur HG n taient possibles avec des relev s perfor s en CAS 9 puisque ce CAS peut contenir de 9 16 relev s par jour cette impr cision rendait difficile et complexe le traitement de ce CAS Nous ne l avons pas retenu consid rant que des cours d eau justiciables du CAS 9 pouvaient aussi bien tre trait s en relev s limnigraphiques int graux En conclusion on peut montrer que tout le travail s effectue dans une boucle par exemple celle des CAS 2 3 4 et 6 DS 170 LA LIN IT IT 24 CAS HG IT HAUT LA CALL CORDEB 170 CONTINUE La sous routine CORDEB est donc appel e dans la boucle chaque valeur de JX nombre r el de relev s d une carte compris entre 2 et 16 La figure 31 montre l organigramme de CORDEB sous routine qui a pour but essentiel de cal c
26. simplifi e de la AQ S AQ _ cote maximale Q 9ABz gt Q 30 AH Jegen 29 cste 16 0 jusqu AQ cste 22 5 Unit s des symboles AQ Qo 6 en m compt positivement en crue n gativement en z d crue AQ 0 47 20 430 5 E A M thode simplifi e 17 32 7 10 7 M thode normale 11 AHxen cm 600 400 200 200 400 600 800 2 Correction de l talonnage par AHx Jaugeage de crue e Jaugeage de d crue 10 Fig 19 M thode de l cart la cote maximale station de Mopti 33 34 Tarage de hautes eaux 200 1000 Hauteur en cm Hauteur en cm 400 00 100 200 300 400 500 600 Fig 20 M thode de l cart la cote maximale Fig 21 M thode simplifi e de l cart 1 cote station de Dir maximale station de Dir 27 az M thode simplifi e 20 normale Correction de l talonnage Jaugeage de crue O Point rectifi e Jaugeage de d crue 22 M thode de l cart la cote maximale station de Dir 2400 2300 2000 Q enmYs Tarage 1968 196 Courbe Qo H a la cote maximale Courbe Qo H tgradient limnigraphique M thode du gradient limnigraphique M thode de l cart la
27. 15 13 JELMOPCJAN 41114 17 14 14 1 2117 25 17 15 TIFLJ 31114916918 16 60 IO 117 25 17 25 17 25 17 17 25 17 25 17 17 1 4 3 125 20 20 18 JF KM M 25 19 19 19 Anst wini J9 KM M 10 JD Tn 21 20 JD LJ H J J9 xJI J e1n Jn JA JM 1 11 TN MINO 2 s Ce 9 IP IN 1 nn 24 I IM 5 IF11 2 JD 23 22 23 1 2 JA Ju 1 GO Tn 24 23 19 309 10 1 J0 1 3 1P JALJM4J 10 10 110 N IFLJA JM S 0 LIMA 1 10 24 CONTINUE 60 TO 27 25 00 26 14 5 26 28461841111 27 CONTINUE 1 12 31124 29 30 28 1F MNNIIS 6 5 EG 1 4 6 1 29 HRITE 652 9v AALJAT 1 Da Tei e pi sJ W9TTE 6 4HAR JA 1 1 1 601 PFTURN FNA Ne 2 gt 2 67 68 PROGRAMME 328 VARIANTE DU 301 POUR LA TRADUCTION DES HAUTEURS EN DEBITS AUX STATIONS NOM UNIVOQUES 2 SYST MATIQUE DE LA CARTE 301 AVEC LES HEHES CONVENTIONS QUE DANS POH 301 LE CHIFFRE 20 DESIGNE UN DEBIT NL HESURABLE NI ESTIHABLE POUR LA HAUTEUR OBSERVEE CALCUL DES DEBITS LA METHODE DIECART AU MAXIHUH IUHIV 1 OU PAR CU GRADIENT L IMNIGRAPHIQUE CI UNI V 33 LA METHODE DU GRADIENT TRAITE PAS LE CAS 9 DE 9 16 HAUTEURS LUES PAR JOUR INTEGER ETAT BASSI STATI CAS ANNEE GROUP CHANG REAL MOAN COHHONJHINI HAUT 116 L MAX HP 16 0664 161 413 15 OMAX KMAXs I COMMONSTAJ Ly ISTAT ve 157 HHAXD 9 HMAXZ KK AA AB IUNTY COMMONHG 300 e CAS ICAs J
28. 37 11 1 69 509 1 640 1 675 35 2 1 29 6 2 105 1 780 31 TABLEAU STATION DE DIRE M thode simplifi e de l cart la cote maximale Sens Ann e No Date H Q Qo AQ 52 52 cm mis 18 5 m s cm cm m s 13 5 C 1951 52 1 26 8 51 375 1060 990 70 7 1 601 226 225 154 155 835 D 2 4 4 52 465 1245 1430 185 12 9 136 22 5 9 6 135 1565 D 3 11 1 54 601 2540 2580 50 15 606 5 1 5 0 0 2580 1953 54 4 23 2 54 542 1815 1915 100 5 2 64 19 2 14 0 270 2185 5 10 6 54 129 120 155 35 22 6 477 225 01 4 0 155 C 6 9 8 54 353 1 100 890 210 236 613 260 225 11 10 900 7 24 8 54 419 1445 119 255 21 5 194 22 5 10 10 1180 8 5 9 54 456 1735 1385 350 25 3 157 22 5 H 2 8 40 1425 1954 55 9 12 10 54 527 2020 179 230 12 9 86 22 5 96 170 1620 10 15 12 54 603 2540 2595 55 21 0 30 5 1 130 2465 D 11 19 1 55 609 2645 2680 35 1 3 4 12 01 5 2675 C 1955 56 12 23 11 55 588 240 2400 0 0 612 24 72 72 175 2225 gt 13 8 12 55 604 2545 2610 6 2 5 8 24 49 130 2480 1956 57 14 7 8 56 247 725 425 300 70 5 565 318 22 5 48 0 200 625 C 1957 58 15 30 7 57 237 465 380 85 224 61 374 225 01 0 380 16 4 8 58 342 1 045 840 205 244
29. 500 sont millimes elles concernent simplement les formats de lecture des cartes d iden tification et ceux des impressions et sorties sur cartes Les programmes adapt s au traitement des donn es de BRE gardent leur num ro et deviennent 310 version BRE par exemple Au sujet des stations ordinaires et des stations de BRE il est bon de dire un mot des pas de temps de d pouillement des hauteurs observ es a Le pas de temps fixe est applicable tous les cours d eau variation pas trop rapide et pas trop irr guli re du niveau d eau c est le cas de la majorit des cours d eau intertropicaux Le pas de temps variable est indispensable pour les tr s petits cours d eau dont les variations de niveau sont tr s rapides et irr guli res c est dire pour les quelles la prise de 16 relev s intervalle constant par jour est insuffisante pour suivre le limnigramme c est le cas de tous les BRE et de tous les petits cours d eau torrentiels ou de zone aride c Certains cours d eau pour lesquels le pas de temps fixe suffit plus de 80 90 des jours de l ann e sont l objet quelquefois et durant quelques jours par an de crues soudaines et brutales dont la pr hension exige le recours au pas de temps variable C est le cas des cours d eau du bassin m diterran en par exemple On dit que ceux ci n cessitent un pas de temps mixte pour leur d pouillement En r sum on peut dire que les stati
30. C 999995 17141724172 172 1 1 1 VAR 190 1 TQ 170 171 VAOLA VAREA CCONDS TF VAR 185 195 195 195 PANTIKAL VARIAPIONCOO GQ TQ 145 18 165 JCan e L I TE XPAX 9 1 00 182 181 186 199 1 3 190 WRITE 7 15 ETAT BASSI STAYI RETALAK y MANT X UD a PICAR ILE LS EMAX 10 121 nn 191 1 3 RIIET ASJETAT BASSI STATI METAL Ky MARTIAL 11649 KL 1 1 R 191 WRIIFLZ 15 35 5 r RETALK LICAPIK L 9 SE spoe VARSA LCCCCO 66 Sun 312 VFPSINN ROF LAS D UNE STATINN A MENTA VARTARLE DFUX FCHELLES IPALTEMFNT Dr RASE RES DONAFES LIMNIGRAPHIQUES INTEGRALES CAS DU RASSIN REPRESFUTATIF EXPERIMENTAL INTEGER 5 1 TUTEGE RE 901915 MINTE MAHTI y ITAR 6 42004 70001 yBANJ01368 y ICJOL 1348 171 12017 MAT 12 312 10012 411 MANOLI ICROI 121 e JORET 2020 NTMENSTIN NANF 12 017040 HEURE 7200 yn 10 344 12 OTRENSION 5 7 DEMAIN DATA 1 32 8 29 119 150 188 7211 242 272 323 3937 1 10716 19116407 COMMIN IDECHyIPENTE JMIN y JMAXs STAT Le ISTAT LMAX HPL16 KMAX 9 A 39 15 1 813 15 013 15 D13 15 HAUTA 6 6 e DER 6 OMAX LONAX 2 0 14 LRMA Ne 16 HAUT 6 1 2 3 5 6 7 8 9 10 11 12 13 14 18 20 21 Toit
31. D 48 30 11 67 701 3 330 3 400 70 2 1 8 3 0 0 9 30 3 430 C 50 22 10 68 650 2840 2 900 60 2 1 2 0 7 2 8 80 2 820 C 51 25 10 68 651 2790 2 905 115 4 0 2 0 7 135 2770 D 53 1 11 68 650 2 900 2 900 0 0 1 0 4 0 4 10 2 910 55 8 11 68 641 2 760 2 815 55 2 0 7 2 6 0 6 15 2 830 57 14 11 68 629 2 610 2 705 95 3 5 10 3 7 0 2 5 2710 59 18 11 68 619 2 500 2 610 110 4 2 12 4 5 0 3 10 2 620 61 22 11 68 608 2 390 2 515 125 5 0 14 5 2 02 5 2 520 p 62 25 11 68 602 2 310 2 460 150 6 1 13 4 8 1 3 30 2 430 22 Sens N D 1 OX la yn Date 26 8 51 4 4 52 11 1 54 23 2 54 10 6 54 9 8 54 24 8 54 5 9 54 12 10 54 15 12 54 19 1 55 23 1 55 8 12 55 7 8 56 30 7 57 4 8 58 26 8 58 15 9 58 28 10 58 13 11 58 17 3 59 23 9 59 26 3 60 23 8 61 10 11 67 1 12 67 4 12 67 7 13 12 67 20 1 68 30 1 68 18 12 68 2 1 69 11 1 69 H cm 375 465 601 542 129 353 419 456 527 603 609 588 604 247 237 342 405 458 530 548 456 463 303 346 537 563 567 574 TT ant F 560 535 525 F 509 E TABLEAU HI STATION DE DIRE M thode du gradient limnigraphique AQ AQ Q Qo AQ To AH a To AO m s m s m s m3 s m s 1 060 1 015 45 4 4 22 21 4 17 0 170 845 1 245 1 510 265
32. F 4 3 On notera que les champs consacr s IDECH et IPENTE 2 I4 variables du 312 emploi de deux chelles sont conserv s entre le champ IUNIV et celui de AA Dans la m thode de l cart au maximum la place des cartes de relev s ant rieurs il y a une carte et une seule de hauteurs maximales non suivie d une carte blanche et cette carte des hauteurs maximales se retrouve chaque ann e avant la premi re carte des relev s de l ann e Ce programme s applique tous les cas de non univocit et aux cours d eau dont le r gime tr s variable et irr gulier n cessite un d pouillement pas de temps variable La structure de base du programme est identique celle du 310 3 de trai tement syst matique des relev s limnigraphiques int graux d nomination des relev s pas de temps variable des cours d eau tarage univoque La m thode de correction de non univocit est celle du gradient limnigraphique UNIV 3 par cons quent la structure de programmation est comparable celle qui est utilis e ce m me dessein dans le 328 Amalgame pour 80 du 310 et pour le reste d l ments du 328 le programme 329 ne m rite pas de ce fait une description int grale et d taill e La figure 34 repr sente l organigramme simplifi de ce programme dans lequel apparaissent surtout les parties qui le diff rencient du 310 Nous allons d crire ces parties en
33. HG couvrent huit jours pleins 192 heures et sont obligatoirement perfor es soit en Cas 2 donc 2 relev s par jour ce qui tient en une seule carte 301 Cas 4 donc 4 relev s par jour ce qui tient en 2 cartes 301 Cas 8 donc 8 relev s par jour ce qui tient 4 cartes 301 De cette fa on les cartes de relev s ant rieurs sont pleines de 16 hauteurs sans quivoque et le pas de temps est calculable simplement l aide d une seule quation On vite ainsi le recours la sous routine GRAD et les complications des CAS 1 3 5 et 7 En effet il suffit pour introduire ces relev s ant rieurs dans le vecteur HG apr s mise z ro de IT d crire la boucle 65 DS 46 J 1 16 IT IT 24 CAS 46 HG IT HAUT J Un test de carte blanche ach ve l introduction des relev s ant rieurs et envoie l instruction 101 de lecture des hauteurs de l ann e traiter est bon d indiquer que peut ensuite traiter plusieurs ann es successives sans repasser par cette pr paration du vecteur HG lequel reste en fin d une ann e N videmment disponible pour l ann e N 1 S il n y pas compatibilit entre la carte de relev s ant rieurs la premi re carte de relev s de l ann e traiter test sur le num ro de station STATI et l ann e LAN apr s mission de messages d erreurs le programme envoie sur STOP car il est impossible de passer l an
34. IF GROUP 11153 51 51 IF GFOUP 10 53 2336333 333 341 342 51 52 5 54 55 56 57 56 400 401 402 2003 2005 M 411 1 2005 10 11 12 14 15 T3 l 697 660 699 702 700 698 740 741 745 750 160 16 15 22 20 23 24 26 30 31 32 3 4 35 36 1 41815 51 52 51 GOTO 341 34253415 529 3419 529 341 341 3 1 52 381 015 68009 16 52 51 51 IF GROUP 15 52 1 51 KHxX 1 GOTU400 Vue 6070400 3 400 KHX 4 60r0400 KHX 5 5070400 KHX 6 GnTo400 KHX 7 6070400 KHX 6 IF CAS 2 599 59 401 CAS 0 DO 402 KO 1 KHX JD 8 CAS GRUUP 1 KD 04120 015 0 0 6010101 TRA TEMENT O CAS 9 IFI 141 312006 2005 2005 0 JO GROUP Dep MH 0 0n4104 1 16 Metz IF DEB J 411 410 410 JF DEB J E0 20 MM HH 3 DEB J 0 D D DE8 J IF 40 LE202AN0 MH GT 0 0J J0 H01 20 1 0 0 60 020 4530 0153 10 IFLMSGT 03011907 MOT SI D H 6070101 PERFORATION DE LA MATRICE 94 40 4015 RISE EN VECTEUR Q NATRICE INTERPOLATION DES DEBITS MAMQUANTS Ra 0 DOL8HDTS 1 12 6810111 12 11 13 11 13 11 11 13 11 1 3 11 015 HAX 31 507016 1815114 15 1 HAX 28 607016 MAX 29 607016 HAX 30 IMAX 15 IMIN 1 Rief KINzKLN L 00698JQ 144 X 091J0 M01S J J0 15 KIN 1 6 0 100 10 560 0 10 210 10 TF
35. IT HG 800 JGRAD Ate IUNIV IPIN HZ 142 HEULG BINIG 0414 2 NHB ALPHA 5 5 845 HRITELT p 13 IRAS ISTAT JANNO 5 y KINZ ICJOURI 1 1 JUDEBI 50 TO 3230 3250 KD K1 15 KF K2 39 1 4 0 1513235 3231 3235 41 3235 2 3232 WRITEIT 11 JETA 1945 ISTAT JANNO HOIS KINZ HAHJOLH s ICJOLH y HK Ds MI M HEKA 16 JUDEBE GO TO 3230 3231 WRITE 7 11 JETA IBAS ISTAT MOIS 9 KINZ IMANJOSM y IC Met De 1KF JUDEBI 3230 CONTINUE 4000 2 008100 5 1 12 X DEBMOI HOIS 1F1X 4650 4660 4660 4650 IDE 10 1 0 1 60704670 2660 1 0 4669 1F X 9995 4700 4 702 4702 Dw 4702 0 1 5 1 1 1 608 TU 4669 D 4700 X 0005 x 20 IDE 1000 X 7 4670 MAMOI MDIS 10E 4100 ICHOICHOISI 1 1F SUDEB1 4801 4800 801 14 4801 QHAX 1000 0MAX 8 4800 IF QMAX 9995 48104815 4815 13 4815 QMAX 0 1 QHAX GLUES 17 50704800 21 4810 QMAX QMAX 0005 X 1000 QMAX ICHAX IU WRITE T 12 JETA IBAS ISTAT JANNA HAHJ3IIRDIS1 ICMDIIMDISI HOIS 1 1 16 123 MAMA Xe ICHAX JUDEBI 22 WRITE 6 5055 FRAG 6 5552 SETA IBAS ISTAT IFIJUDEB1 4910 4900 4910 15 4900 HAITE 6 50561HILE 19 GOTO 4920 8910 HRITEC6y 5057 HILL 9 4920 CALL YMPRI MILL MAJ ICJ MAHOL 10701 XMAX 0 0001 FLOAT MAMA X 10 FLUAT LCMAX 1 XUAX XHAX 0 01 HRITELG 5100 XHAK 5100 FORMA
36. T O M Cah Hydrol vol V n 3 3 Roche M 1969 a Traitement automatique des relev s limnigraphiques int graux Cah Hydrol vol VI n 2 63 64 4 SIRCOULON J 1970 lt Programme mixte de transformation hauteurs d eau d bits Cah Hydrol vol 1 5 DUBREUIL 1969 lt Premiers mod les de cartes m canographiques du Service Hydrologique de l O R S T O M pour l tude des bassins repr sentatifs ou exp rimentaux impressions du ler semestre 1969 Note technique n 18 O R S T O M Serv Hydrol Paris ANNEXES MISE EN PADABOLE DES TACIS CnUP3ES TARAGE L RVELCPPESDENIVELES HARPALE CORRE CORRECTIONIN NE Station VARIANLE A DEUX FCHELLES ETAJyBASST9STATT DIMENSION Hp 161 9A 61 HJNTI151 D NT11 C 3y15 Ht200 l DEST 200 16 ONNPELE ES 16 D I AT 15 OCMIRT LS y 20 200 01200 EE TIMENSICH 1 MARTI 15 FRAGILA 3 15 12201 FORMAT 12 1668 2 FCRMSTILEFS 3 FORUAT 15F4 2 FARMATIASES 3 FORMAT EN Le 51396 1e 70202207 VALEURS DE 71 FORMATI 19 11 11 16 12 1 16 613 11 FORVATE 1 1X BAREMF HAUTEURS_DERIT51 FORVAT 21 9 HAUTEUR M1 15 NESITIM3 S FORMAT LXrF27e2rF 2803 FOPVAT 1X 1X VALEURS DE 1 26 41445 1 HPEL EN METRES 124X 12 F14 2 11 2
37. X E0 20 1 E 20 0 al 10 1 1 99952700 702 702 X 0 1 x 5010699 X 00054 X 0 1000 IDE DE HANTIJ IDE 1 34 1 127404740 741 LIM 15 GOTO745 LIMEPAX 15 2 7 3 2 4 8 5 15 015 MANT J I ICARIJ J 1 LIMA IFIKIN 21750 760 760 IMAX MAX IHIN 16 6010691 D018JD 1 MAX P H L Q H QJ JO HG1S TNTERPQLATIGN JJ 1 00218 1 JOAN 1 1 22 22 23 1F Q M 20 21 21 HO M Lis 607021 TF Q 4 21 24 24 COP Q M Q MD H HD HDP HD 1 1 HF Q NI QUN L1 COR List CONTINUE CALCUL ET PERFORATION DES DEBITS MOYENS MENSUELS Inte IHAX 0 00379015 1 12 31 32 31 30 31 30 31 31 30 31 30 31 MOIS MAX 30 601035 HAX 3L cornas 1 1815133 34 33 HAX 25 601035 HEH IMAX INAX MAX Se naa6w LMIN IMAX S sS Q M IMIN LFI 5 40 42 Ku 37 43 45 44 1 0 0 0 9995141 62 42 QHO 1 090 I 1 1 GOTO4n GHO 0005 DE 1090 QHO HAQ OLHOIS DE 1 0801 0152 1 1 0 1 8995 44 45 4 5 QMAX 1 QHAX 1 1 1 601043 0005 DE 1C00 QHAX PAQUA DE ICOHA URITE 7 4 JETA 1845 15741 MILL 5 y TCQMOCHOIS I 015 1 1121 CQHA 500 561 560 563 562 564 509 520 510 550 600 801 802 803 804 805 806
38. cote maximale M thode simplifi e i e Jaugeage de d crue Jaugeage de crue 500 520 540 480 Hautepr Tarage 1968 1959 A T C CI LE NIGER A DIR Hauteur en cm 3300 Tarage 1967 1968 2900 Hauteur cm 710 690 700 LE NIGER MOPTI Fig 23 et 24 Trac s sch matiques de tarages de hautes eaux 5 Les programmes 36 de traitement automatique pour stations non univoques 5 1 Le catalogue des programmes On trouve galement sur le tableau d j cit les gains d efficacit des deux modes d emploi de la m thode de l cart au maximum compar s ceux dus la m thode du gradient La m thode normale de l cart la cote maximale para t tre l g rement plus efficace que les deux autres Elle peut donc tre adopt e de pr f rence pour ce type de cours d eau d autant plus que pour le calcul automatique le mode simplifi n est d aucune utilit la fonction Arc tangente existant en biblioth que de programmes Les programmes ont deux fonctions selon qu ils effectuent la mise en quation de la courbe ou des courbes d talonnage ou qu ils traduisent les hauteurs d eau observ es en d bits Pour les stations univoques le catalogue des programmes est assez simple et d j en grande partie bien co
39. d une journ e au minimum il y ait un relev 0 heure et un autre 24 heures afin de faciliter le calcul du d bit moyen journalier Ceci fait que le relev du jour J 24 heures est suivi du relev du jour J 1 0 heure et que ces deux relev s sont iden tiques tant en hauteur qu en date si on convertit celle ci par rapport une origine des temps Afin de ne pas traiter le 2 relev qui est identique au 197 la fois dans le calcul du vecteur HG et dans la recherche du gradient limnigraphique nous avons bas la constitution du vecteur HG sur cette particularit ce qui nous permet de consti tuer la base de temps uniquement l aide de l heure HEU et des minutes MIN qui accompagnent chaque relev HAUT et n gligeant le jour On calcule ITP base de temps journali re provisoire en dixi mes d heures Cette variable varie de 0 240 A chaque passage par z ro on vite et l introduction de la donn e dans HG et la recherche du gradient limnigraphique DH op rations d j faites pour la pr c dente valeur de qui tait 240 et l on augmente d une unit le compteur totalisateur de jours MZ qui d clenche un mouvement r sultant dans le compteur totalisateur de dixi mes d heures IMZ Le bloc d instructions qui correspond ce travail est le suivant emprunt aux relev s ant rieurs mais identique pour tous les EE et alors SH au d but de CORDEB 32 D 37 1 6 ITP 10 H
40. d tail dans la mesure o elles sont elles m mes diff rentes de ce qu elles sont dans le 328 cf 5 4 3 59 60 Le support de l information hauteurs est la carte 311 dite de relev s limnigraphiques int graux La conception particuli re de traitement de la m thode du gradient limnigra phique dans le 329 est comparable celle qui est la base du 328 a Cr ation d un vecteur hauteur temps HG contenant des hauteurs adres s es en temps partir d une origine quelconque b Introduction des relev s limnigraphiques int graux dans ce vecteur c Recherche du gradient limnigraphique et calcul du d bit corrig Les particularit s du 329 sont essentiellement dues au fait que les relev s limnigraphiques int graux sont d j une succession de couples hauteurs temps ce qui modifie les processus de traitement du vecteur HG En dehors des op rations d initialisation du vecteur HG en cours de traitement de relev s toutes les phases b et c pr c dentes sont effectu es dans la sous routine CORDEB l int rieur d une boucle D 100 J 1 IFIN qui traite en une seule fois les IFIN couples hauteurs temps 6 en g n ral moins de 6 parfois en fin de mois d une carte COH 311 La succession des donn es d entr e s apparente celle du 328 Carte d identification de station 5 cartes de tarage lues en TARAG dont la 5 contient le
41. dans le bief entre les deux chelles par exemple Enfin il est important de pr ciser que l tude concerne uniquement les stations non univoques stables Le tarage des stations instables ne pose pas quant son trai tement en calcul automatique de probl mes particuliers puisque l on peut consid rer du entre deux ruptures de stabilit une station instable est alors stable univoque ou non On remarque ce propos que les cours d eau sont instables quand leur lit et leur contr le aval sont sujets des variations brutales et peu pr visibles souvent apr s une forte crue tels que mouvements de fonds par creusement ou exhaussement et que si le ph nom ne causal est de m me nature que l un de ceux qui est responsable de la non univocit de stations stables modification g om trique du contr le il est d intensit bien sup rieure Le probl me hydraulique du tarage d une station hydrom trique deux chelles situ e dans un bief contr le aval variable et dans lequel l coulement n est pas uniforme a t trait depuis longtemps et ne pr sente aucune difficult On en trouve une excellente d monstration d taill e dans l ouvrage de G REMENIERAS 1 2 partie chapitre Nous nous contentons ici d en pr senter un court r sum L hypoth se du r gime uniforme permet l emploi la formule CHEZY pour exprimer la vitesse du courant en fonction des caract ristiques g om trique
42. e c est dire la m thode it rative cf 3 La traduction de l quation 5 sous une forme traitable en programmation conduit Ah 12 Q Q 1 5 8 9 Pour r aliser la correction du d bit l aide de l quation 8 il faut et il suffit D avoir la vitesse moyenne U Ya b De remplacer Ah dans l quation de correction par La vitesse moyenne est li e la hauteur La forme de liaison est parabolique Comme il s agit d une correction de second ordre il n est pas grave de simplifier et par exemple de remplacer les paraboles par une s rie de droites On peut donc op rer comme suit Dans TARAG lire les param tres des diverses droites et les valeurs HB des bornes inf rieures de hauteurs limites d application de chaque droite Cette lecture peut se faire dans le champ libre de la carte 513 Dans CORDEB calculer H et l introduire dans la correction de d bit L emploi d un test ad quat permet de laisser ces instructions dans le 329 et de les utiliser ou non selon que l on recourt ou non l emploi de la vitesse dans l quation de correction du d bit 1 REMENIERAS 1960 L Hydrologie de l Ing nieur Coll Lab Nat Hydraulique Eyrolles diteur Paris 2 RocHE 1968 Traitement automatique des donn es hydrom triques et des donn es pluviom triques au Service Hydrologique de l O R S
43. eau importants pente faible Elles sont d ailleurs toujours cit es en exemple pour caract riser les raccordements des courbes de crue et de d crue Point de rebroussement pour Dir fig 2 Raccordement arrondi pour Mopti fig 3 On a galement fait l application l oued Boitiek mod le du petit cours d eau irr gulier faible pente d j analys par la m thode it rative cf 3 2 Apr s essais les intervalles de temps At retenus ont t 5 jours pour Mopti et Dir 6 heures pour Boitiek L ensemble des calculs d approximation pour le trac des deux courbes consti tuant le tarage est pr sent en d tail dans les tableaux suivants IT pour Mopti pour Dir et IV pour Boitiek Les courbes de tarage et de correction par le gradient font l objet des figures 9 10 pour Mopti 11 12 pour Dir et 13 14 pour Boitiek Dans ces tableaux outre les symboles d j explicit s nous utilisons les abr via tions suivantes jaugeage de crue ou en tale fin de crue D jaugeage de d crue ou en tale d but de d crue N num ro du jaugeage les num ros manquants correspondent des jau geages soit manifestement erron s soit effectu s en un autre site cote l chelle lors du jaugeage d bit mesur AQ Ah a ordonn es de la courbe de r gression ici la droite o k x AQ AQ 09 Qo OI d bit corrig apr s r gres
44. en rempla ant le gradient limnigraphique par l cart la cote maximale Ahx TABLEAU V SYNTHESE DES GAINS COMPARATIFS ENTRE LES METHODES DE TARAGE Station Mopti Dir M thode 1 D I Ecarts sur AQ ZAC 3 605 3 729 3 729 5 003 5 160 5 160 AQ 1 525 1 178 1 379 2156 1 835 2055 Ecarts sur AQ Qo 5 Gi 359 8 396 7 396 7 413 5 462 3 462 3 0 D 2 192 7 116 5 127 2 178 4 151 0 160 1 R duction des carts sur AQ 57 8 68 5 63 2 57 0 64 5 60 1 en 46 5 70 5 67 9 56 9 67 4 65 4 Moyenne des carts AQ m s i 75 78 78 152 156 156 32 25 28 65 52 62 i 7 5 8 3 8 3 12 5 14 0 14 0 TQ f 4 0 2 4 2 7 5 4 4 6 4 9 Ecart pour 50 AQ mis i 60 60 60 145 130 130 des points 17 10 10 40 25 20 AQ i 4 0 4 5 4 5 8 3 11 7 11 7 Qo f 1 3 1 0 1 1 2 8 1 8 1 8 Pourcentage de points 33 31 31 3 12 12 2 62 71 62 37 51 55 M thodes I du gradient limnigraphique 1 initial de l cart la cote maximale f final simplifi e de l cart la cote maximale He 1 AHx He NI sd AHx 4 Hi 3 Exemple LE NIGER A DIRE Aly Ki 2 1854 _ 1955 1967 1868 KE 9 N re minnale 1942 1943 200 100 0 100 200 T eh lours Ten lours Fig 15 Homoth tie des crues des grands fleuves Fig 16 Passage de la m thode du gradient limnigraphique faible pente et plaines d inondation importantes la m t
45. groupes de lt date hautcur gt de chaque ligne ou carte sans aller une nouvelle ligne ou carte chaque jour moins que l on ne soit en d but de mois Les relev s HAUTB quand ils sont n cessaires doivent correspondre exactement la m me date que ceux de HAUTA c est dire la m me heure m me minute du m me jour Les relev s HAUTB sont donc sur des cartes pleines six groupes de lt date hauteur moins que l on ne soit en d but de mois Les cartes relatives HAUTA sont num rot es de 1 partir du ler jour de l ann e pour lequel on poss de des relev s quel que soit ce jour Les cartes relatives HAUTB ont le m me num ro que les cartes HAUTA correspondantes auxquelles elles s appliquent augment de 5000 nombre que n atteint jamais la somme des cartes HAUTA pour une ann e quel que soit le type de cours d eau Les informations relatives aux talonnages sont contenues dans 13 cartes dont les 12 premi r s ont t produites par le programme de tarage 311 Elles con tiennent pour chacune des 3 courbes d talonnage les hauteurs limites 1 carte COH 306 ou 506 et les coefficients des paraboles 3 cartes COH 307 ou 507 La 13 carte est la 513 cf fig 26 sur laquelle sont port es essentiellement les variables IDECH et IPENTE a IDECH diff rence de cote exprim e en millim tres entre les z ros des deux chelles et format 14 IPENTE d
46. identification bassin po 501 Fin de station non ecture carte identification 502 station Lecture talonnage Jud bi Lecture carte 311 Hauteurs dat es 100 Etat univoque Lecture nouvelle carte 311 des hauteurs de l chelle principale 310 INI N DEBU 2__ Traitement des donn es d une ann e P H 310 de fin d ann e non D 6000 1 16 11 Contr les d ann e de Lecture carte Con 311 Hauteurs dat es station de num ro de carte de s quence de carte 310 non Calcul des d bits dates Carte blanche de ja carte lue IE ou gt 0 5000 Call TRAD Mise en vecteur des d bits Etat univoque dates perforation impression 310 des donn es de Call TRAD la carte jue GOTH 110 Fig 28 Organigramme POH 312 version BRE 45 46 l information hauteur de la Ire carte lue est envoy e dans HAUTB 1 Tandis que si le num ro de carte NOCAR est inf rieur 5000 l tat est univoque et l on inscrit l information hauteur de la lre carte dans HAUTA D La traduction en d bits et le contr le de changement d talonnage s effectuent ensuite ensemble Un nouveau test sur le num ro de carte celui de la 17e carte lue dirige selon les cas vers la sous routine TRADA ou la sous routine TR
47. importantes et lorsque le r gime d coulement dans le bief est non permanent parfois les deux causes interf rent contr le variable et non permanence de l coulement 1 1 1 UN CONTROLE AVAL VARIABLE Un contr le aval peut varier soit par modification g om trique soit par modifi cation hydraulique La modification g om trique peut tre naturelle c est le cas du banc rocheux ou de graviers qui s rode lentement c est le cas de la v g tation aquatique du lit des cours d eau faible pente dont le cycle annuel de croissance est plus moins r gl sur le r gime d coulement qu il perturbe de mani re pseudo cyclique chaque ann e mais avec des intensit s variables La modification g om trique peut tre artificielle c est le cas du barrage seuil mobile faisant fonction de contr le La modification hydraulique s observe dans le cas d ennoiement d un seuil de transformation d un coulement torrentiel en coulement fluvial par exhausse ment du niveau cas typique des stations situ es dans le remous d une confluence etc Dans tous ces exemples de variation du contr le le r sultat hydraulique est une modification soit de la vitesse soit de la pente d coulement Pour les cours d eau contr le aval variable le tarage est d autant plus complexe que 1 variabilit du contr le est grande On peut sch matiser a priori un tel tarage en le repr sentant comme une somme juxtapos e de
48. naturelle pour qu en l absence de mouvement du contr le aval il puisse devenir univoque donc admettre une courbe enveloppe Le tarage complet d une station non univoque deux chelles comprend donc au total trois courbes Une courbe enveloppe dans la plupart des cas Une courbe d nivel e normale Une courbe de correction de pente Le tableau suivant explicite les limites d emploi de la courbe enveloppe au lieu et place des deux autres Zone Zone d application Zone sans d inversion du tarage non univoque influence de pente 2 chelles de pente D nivel e D D lt 0 0 lt D lt hm D gt hm Hauteur amont lue hm hm lt hv dZ hv dZ lt hm hm quelconque Hauteur aval lue hy gt dZ hm dZ gt h gt dZ hy lt dZ Courbe d nivel e normale Ccurte Courbe de correction de pente enveloppe 3 Approche th o rique du tarage des stations non univoques dot es d une seule chelle 3 1 Analyse hydraulique Nous avons trait ici le cas g n ral nettement le plus r pandu c est dire celui des cours d eau diminution de pente non univocit partielle de type B mais il va de soi qu un raisonnement analogue et sym trique est applicable aux cours d eau accroissement de pente non univocit de type A par crues h t rog nes par exemple comme le Wabi Shebelli si ses stations taient dot es de deux chelles La courbe envelop
49. ou une surface de tarage Pour certaines stations il faut effectuer des jaugeages pr s du maximum de hauteur chaque ann e pour esp rer avoir un tarage satisfaisant Pour toutes les stations non uni voques la transformation des hauteurs en d bits est laborieuse et impr cise Cet article rassemble le travail de plusieurs chercheurs qui ont eu pour objectifs de trouver des solutions hydrauliques simples pour l utilisation des courbes de tarage non univoques et de permettre un traitement automatique des donn es en harmonie avec celui utilis pour les stations univoques Apr s une pr sentation des divers cas rencontr s de non univocit du tarage la premi re partie d veloppe les m thodes de r solution tudi es pour les stations dot es de deux chelles puis pour celles qui n en ont qu une La seconde partie expose les programmes de traitement automatique appliqu s aux diff rents cas de non univocit tant pour le calcul du tarage que pour celui des d bits Outre les signataires de l article cette tude a b n fici de la contribution de Y BRUNET MORET qui a propos une approche hydraulique th orique de la solution du tarage des stations dot es d une seule chelle Cf chapitre 3 3 1 Les diff rents types de stations non univoques 1 1 Les causes de non univocit Une station hydrom trique est non univoque lorsque le r gime d coulement tant permanent son contr le aval est sujet des variations
50. pas de temps 24 CAS et 24 1CA Appel C RDEB ss debit moyen journalier Fig 30 D tail du 328 calcul du d bit moyen journalier 49 50 Subroutine TRAD D 130 J JMIN JMAX d bit DEB J cart au maximum selon valeur KK CONTINUE Calcul correction de d bit C RQ et d bit corrig Tri du d bit maximal Q MAX Subroutine C RDEB Calcul du moment IG pour calcul du gradient limnigraphique lt 0 gt 0 D 611 115 1 16 Recherche d une hauteur r elle HG au moment IG ou anterieure a IG CONTINUE oui Calcul gradient ant rieur DH1 non DS 6113 HI1 ITG ITI Recherche d une hauteur r elle HG post rieure au moment Calcul de CBRO 6013 Calcul du debit corrig Tri du d bit maximal Calcul du 8012 gradient post rieur DH2 et de DH D 6114 J2 172 171 Changement d origine de temps du vecteur HG IT recul 00 pa CONTINUE 16114 D 6112 J1 IT3 ITI Effacement du vecteur HG entre les pas 201 et 270 C NTINUE IF ITI 270 RETURN 5026 END Fig 31 Sous routines du 328 5 4 2 TRAITEMENT PAR LA M THODE DE L CART AU MAXIMUM UNIV 1 Cette m thode s applique aux cours d eau limnigramme r gulier d aspect sinusoidal La correction non univocit s eff
51. prend comme valeur le num ro de la position dans la carte pr c dente 301 de la veille du changement crue d crue La figure 30 d crit en d tail les circuits d appel aux sous routines TRAD et TARAG dans tous les cas possibles selon que ICHA et IGRAD sont nuls ou non s par ment ou ensemble et selon leur place respective dans ce dernier cas chaque passage pour IGRAD 0 le compteur KK augmente d une unit avant l appel TRAD La sous routine TRAD cf fig 31 fonctionne comme celle du 301 pour le calcul du d bit DEB J relatif une cote HAUT J Ensuite elle calcule l cart au maximum DHMAX selon la valeur de KK puis la correction de non univocit CORQ estim e ATAN AB DHMAX l aide de l quation en arc tangente et des param tres AA et AB lus dans TARAG Le d bit DEB J est alors corrig et le d bit maximal est tri DEB J DEB 1 CORQ 100 Le calcul du d bit moyen journalier qui s effectue ensuite dans le programme principal apr s le 148 est identique celui du 301 quel que soit le CAS 1 9 Les donn es se terminent par Une carte blanche en fin d ann e test sur ETAT et envoi l instruction 10 de d but de traitement annuel des d bits Une carte blanche en fin de station test sur la carte des hauteurs maximales de l ann e suivante Une carte blanche en fin de donn es t
52. tarages pseudo univoques pente constante selon les conditions de contr le le tarage instantan occupe une position quelconque de la surface de tarage et une onde de crue peut la balayer en totalit ou en partie ex fig 1 La non univocit peut se limiter une plage r duite et m me deux courbes pseudo univoques correspondant deux pentes privil gi es des coulements deux tats du contr le aval par banc rocheux par v g tation etc 1 1 2 UN R GIME NON PERMANENT On note galement des modifications de la vitesse et de la pente d coulement dans le cas de non univocit d une station situ e sur un cours d eau r gime non permanent Les conditions varient entre les phases de crue et de d crue et d autant plus que les ph nom nes sont accentu s On rencontre deux groupes diff rents de cours d eau pour lesquels la non univocit est due la non permanence de l coulement Tout d abord les cours d eau tr s faible pente et r gime hydrologique r gulier c est dire pr sentant un hydrogramme annuel compos d une mont e et d une descente pratiquement sans accident C est le cas bien connu des fleuves tropicaux faible pente et plaines d inondations tels que le Niger et le S n gal La r gularit de l hydrogramme annuel fait que l on observe des r gimes d cou lement de crue et de d crue quasi homog nes d une ann e sur l autre c est dire 1 2 Ess
53. 0 C 22 22 8 59 515 2000 1725 275 15 9 175 114 4 5 80 1805 D 26 14 3 60 142 120 145 25 17 5 690 548 16 7 05 0 145 1959 60 27 31 3 60 120 84 101 17 16 8 570 16 8 0 0 101 28 22 4 60 92 55 64 9 14 1 598 16 9 2 8 2 66 29 11 5 60 74 39 47 8 17 0 616 17 0 0 0 47 1960 61 30 23 7 60 376 1 075 910 165 18 1 691 315 14 6 3 5 30 940 D 1962 63 31 9 6 63 140 130 140 10 7 1 702 52 167 9 6 14 154 32 24 5 67 52 44 38 6 158 676 172 14 0 38 33 15 10 67 690 3370 3285 85 2 6 720 30 2 6 0 0 3285 34 17 10 67 693 3340 3315 25 0 8 27 24 16 55 3260 35 19 10 67 696 3450 3350 100 3 0 24 24 0 9 30 3350 gt 36 24 10 67 702 3465 3410 55 1 6 18 1 6 0 0 3410 37 29 10 67 708 3605 3470 135 3 9 12 11 2 8 95 3565 38 1 11 67 710 3565 3495 70 2 0 10 0 9 1 1 40 3535 39 4 11 67 713 3 540 3520 20 0 6 7 07 01 5 3515 1967 68 40 7 11 67 716 3585 3560 25 0 7 4 04 0 3 10 3570 41 11 11 67 719 3615 3610 5 01 1 04 0 0 3610 42 13 11 67 720 3640 365 5 0 1 0 0 01 5 3640 D 43 16 11 67 720 3630 3635 5 01 0 0 0 1 5 3630 44 19 11 67 718 3610 3595 15 04 2 02 0 6 20 3615 45 22 11 67 715 3 440 3545 105 2 9 5 0 5 2 4 85 3460 46 25 11 67 711 3 460 3500 40 11 9 08 03 10 3490
54. 0 HEUt SI4HINT J Q KXI NF81J CONTINUE ICARTESNOCART 110 INI NOOERU 2 JORTINI LI INT 1309 IN INI K 40010 10 1320613100 1320 308 10 08 10 1 60 1300 JOREINI JORLIO 1 1 CONTINUE HRELEV JOR JFINODEBU 1001 10101001 101005 J 1rNODEBU QJ0 31 10 TEUIBIS 1011e 10121011 NJ0 365 0701013 NJ0 386 TF INRELFV NJO 1020 1030 1025 MRITES 6 TIRILL 2000 HRELSRREL EVER DO 1022 J NRFL NJD 93063 10 M N0DEBU S 0 S S HEURE M L HEURELH 0 M 01M 1 172 IFCHEURELM 13 Laag 11050 1060 1060 1 1 JORIM 5 QJOIL S 1 40 IF 1 NRELEV 1049 1070y1070 0811003 1 0 X 0J01 J 1F X 1650 1660 1660 TDE 10 1 0 JF X 997511700 1702 1702 0241 T 1 1 60701699 0005 10 lt 1000 BANJ0 J InE ICJOtJ I JANNOFMILL 1900 JANNA MILL 1000 8 NC NC 1 3050 J 1 16 1 16 3 16 0 1 IC J 1CJOL1 IF 0 22 306043066 2970 WRITELZ 9 lt FJA 1 ICJ Jale 16 e JUDERT 63703110 11 6 252 PO 308N J 1 LIM 1 4 416 3 16 MALJI MANJOL 11 JC J ICJjQ I HRITFITr KRAS IBVREP e 121 SSL SL ai 113081 3082 2081 11616 7011 5 JAVREP MILL G0TA in MPITE ASLA KETA KBA lt IBVREP MILL 841 6 1012 LS TC DT 1 301 004700 01551 12 0
55. 0471042 3120 312043110 3100 3110 3100 31 7 3110 3100 3110 3200 MOIS NJHN 3 GNIN3125 3110 NJMN 30 60703125 2120 513121 7122 3121 3121 NIMO 28 60703125 3122 NJMN 2 3125 Kl K 1 2 K NJHO TEST 0 TOTAL D n03200 J NJMQ HN K J MAJ HOIS J MANJOIMNI IGJ HO S J ICJOUHN YeHANJN MN 10 tICINEMN 3 1F Y 3130 3200 3200 3130 TEST 10 3200 TOTALSTATAL Y 5713210 3220 3220 3210 DERKOI MOIS 104 50 TO 3226 3220 DFAKOL MOTSI TOTAL FI OATINJHOY 3225 DO 3230 XINZ 1 2 113280 3240 3250 3240 0 1 KF KD 14 WRITELT 18 KRAS TBVREP MOTS KJNZ MANJO H 3 1 JUDERT GOTN323nN 3250 KO K1 15 KF K2 KO KE KO IFIKD 15 3235 3231 3235 2235 KAminel 3232 WPITELZ s11 KRAS MOIS KIN y IC s MEKO ATEM e 104 MEKA 16 JUOFBL GATA 3230 3231 WPITEITs1LI KRAS IBVREP e 15 e KINZ y ICJOLM e M KD l JUDERI 3230 CNANTINUE 4000 2 n0410Q 15 1 12 X DEBHNI MOIS 1F1X 4650 4660 4640 4650 INF 10 4660 1 0 4654 1F Y 999514700 4702 4702 4702 Y 0 1ex 60 TQ 2659 4700 x n005 10 1000 467C MAHOI HOIS InF 4100 ICGRDIIMOIS 1 0 136801 4800 4801 4801 QHMAX 1000 QMAX 4890 9995 4810 4815 4815 4815 QMAX 0 1 0OMAX JU 1U 1 4810 QMA
56. 05 IDE 1000 X MANJO 106 DST JAUNO HT LL 1900 JANNA RILL 1000 NC 0 DO 3050 3 1 16 4 I NC 16 J 16 HA J MANJO I IC S 1CJ0 1 IF NC 2233060 3060 3070 7 93 1 5 104 Ai Asil 016 1 e JUDEBI 60703010 LIM NJO 352 DO 3080 3 1 119 I NC 16 J 16 IC J ICJOtI 9 JETA IBAS ISTAT 104 sJ 1 LIM IFCJUDEB1 3081 3082 308 6 7011 IBAS ISTAT MILL 60703083 6 101 I BASS STAT e HILL HATTEN EK KE GO EH K O 004000 HOLS 1 12 0 03100 31203100 3110 3100 3110 3100 3100 3110 3100 3110 31001 1 5 3100 3110 3120 3121 3122 3125 NJMO 31 60703125 NJHO 30 60703125 1F 101513121 3122 3121 NJHO 28 60703125 29 Kl K 1 KZ lt K NJHO TEST 0 TOTAL 0 003200 J 1s M3MO HNak S HAAL RDISe 3 MANJOTHN JC3 RO S J3 ICJO HN 72 YSHANJO HN 10 1CIOCHNI 3 119 LekMAX 331309 3200 3200 120 H HAU HP EL 3130 TEST 10 DEBIJI H ACLoLIFHeA 2 LI IHAL L 3200 TOTAL TOTAL Y 130 CONTINUE IF TEST 3210 3220 3220 RETURN 3210 DEBHOL MOIS 1l0 END GO TO 3225 3220 DESHGI MOLS TOTAL FELOAT NJKGI SUBROUTINE CORDES 3225 DQ 3230 KINZ 1 2 D INTEGER HEU 2 113240 3240 3250 COMMON AR Ne 1 15 1 16 3 153 3 15 3240 KD K1 1162 OEB 61 QHAX
57. 12 14 6 18 5 LL Eelere 6113 340130680333 1 7120 17 FORMAT 1 s NUMERO DE 1 212 14 amp amp 17 068115 1 JOURNALIERS EN LIFRES SECONDE 8 410 13 11 4 3 56 10 1391138 103 READIS 1IJETA IBAS ISTAT FRAG DO 2109 I 1 800 HG 1 9 99 IF JETADLOO 9999 100 STOP CALL TARAG READ 5 19 JUDESI FORMAT 1 Gr 15 IDAS ISTAT FRAG IFCJUDEBT 51 55 51 MRITE 6 70013 4 22 5 INSTANTANES EN H3 S GOTO 60 1 6 1002 7 7370 8215 INSTANTANES LIYRES SECONDE QHAX 00 Mis EMZ 0 De 0 READIES 21 5 NOSTA JANN 5 1 TD MINAT Je HAUT IIIe Liste 6 JUHITEs ICHANGe NOCART L CHETAI30 11030 LECNOSTA SI STATI31 232031 VRITEC69 14 FORHAT LX ERREUR DE STATION POUR LE GRAOIENT TO 9999 00 37 1 1 6 STP LOSHEUCT MINT 13 6 IF 1TP 38435 38 2 1 ts Zant JY ITP IM2 HG IT HAUT I 60 TO 29 READ 5 2IMETA IBAS NOSTA JAHN Ss I I e MINEI Io HAUTE 6 11219 61 9 JUNITE ICHANG y NOGART 111 10004111 DO 8000 I 1 6 00811116001 6001 6000 Ra 60 TO 6002 CONTINUE IFIN 6 IF ISTAT HOSTA 112 115 112 WRIYEt 3 JANN READ 5 2 NETA 113 110 113 FORMAT 1 1844 FOURMAT 1 Xy HUMERO DE STATION 76X 212 14 13 FORMAT T2
58. 2 DEBETS MOYENS JOURNALIERS El 1X 14 T136 M3 5 11 FORMAT T22 DEB TS HOYENS JOURNALIERS 1 19 73 1 178 5 5 4 FORMAT 77757 JANY 713 FEVR 119 HARS 725 AYRI 9 TIL 737 Le JUIN t f ie JUIL rT49 AOUT 155 SEPT Se 761 OCTO T67 e 173 2 FORMAT 1X 13 1X FORMAT Or De Lt PFORMAT 2 Xy SFCHOCART 1 16120 ie 11117 118 117 117 HRITE 6 6YICART 118 119 120 125 6070113 MILL JANN 119 120 119 RRITE 6 HILL 6010113 MILL JANN 18 5 JANN 4 44 JANN LI 161 125 16 300842 1 127 130 127 130 1 09015 12127 150 127 150 NODEBU O GOTO 160 127 1FtIBISIBOP81 80 80 RODESU NODEBA MO1S JDUR 1 8 GOTO 160 HODEBU NODERA HOLS JOUR 13 1 160 K N0DEBU 161 JMIN 1 IF 4UNITEM162 154 162 162 00163 1 1 6 163 164 165 171 HAUT LI HAUT L 0 1 165 170 165 5 5 ICHA JMAX ICHA CALL TRAD CALL TARAG LFCICHA 6 1729 160 180 JHIN ICHA L 170 180 650 655 660 1 699 702 700 670 698 190 195 196 SHAX IFIN CALL TRAD CALL CORDES 00 698 Ae EIN IFISUDES1 50 655 650 8 1000 X DEB J 0 IF X 9995 700 702 702 Ze Ly 1141 GOTO 659 0005 1932100084 ICAR J I IF IFIN 61190 196 196 JFINRIFIN 1 00 195 4 1 1 6
59. 3 5 5 1 20 29 RETURN CONTINUE END RETURN END 6005 6006 6007 6111 6008 6009 6011 6113 6012 6010 6012 6004 6024 6114 6112 6026 SUBROUTINE CORDES INTEGER CAS COHMHONJMIHe JMAX HAUT 16 HP 116 DEB 16 3 15 I COMMONSTATT ISTAT C 3 15 y HHAXO HHAX Lo 2 LUNIV COMHONHG 300 CAS ICA JGRAD AC CORO IT ROUS 4015 15 Ae IX 11 IT IT1 JGRAD 1 1 1 IDE8 DBBILA 100 10 8 6005 6005 6006 0 TO 600 DO 6111 11 1 16 L2 ITG L1 IF1HG L2149 98 6111 6111 6007 21 99 9816008 6111 6111 CONTINUE 60 6009 DH1 HAUT LA HG L2 IF0L1 1 6010 6010 6009 DO 6113 M1 ITG IT1 IFCHG H15 9 98 61136113 601L1 TF HG M1 99 9816012 6113 6113 CONTINUE DH2 HAUTILA HG H1 H2 H1 1TG 1 DH DH1 M2 DH2 L1 1 L1 1 M2 50 TO 6013 DH DHI AC DH OE8 LA 0EB LAI 1 CURQJ IF DEB L2 LT 0 1DEB LA1 0 IF QHAX LT DEB LA1 QMAA DEB LA IF IT1 270 6026 60269 6024 172 171 200 IT 0 DD 6114 J2 IT2 IT1 IT IT 1 R HG J2 HG ITi R CONTINUE IT3 1T 1 00 6112 S1 173 111 HG J1 9 99 RETURN END 2 12 GO 23 21 23 22 23 22 23 23 22 23 22 23 015 nnn 55 7001 51 7002 60 29 30 31 14 32 35 30 37 110 111 6001 6000 6002 112 113 5055 5552 5056 5057 5058 5070 5071 5072 115 116 PROGRAMME
60. 3 24 25 26 97 98 30 34 96 110 120 130 140 150 299 200 220 230 210 194 59 200 2 9 soi 2701 an2 2092 253 2003 FORMAT 212r14r18A4 212 14413411 214458 12 FORMAT 212y 16 13 12 16144 1X 12 FORMATI2I2 4 3 1263116 111 F RVATI 1 12120 1411846 5 X ETAL M I 14 F30Q9AJ112 16F8 31 FORMATILX 1 STATION PENTE VARJARLE PE DEUX ECHFLLES FORSAT LEF4 3 1564 1 FORMAT 1 In 1X BARFMF 5 DENIVELEES_RAPPORTS 5 FORYAT 2X 1Xp VALEUR DE D9NPIL 25X Lt 5X RAPPCRT DC 77 1 24 2 14 3 FORMATI21X RAPPORT D DN 13Xp RAPPORT O QN FORMATA ISX ICOURRF FHVELOPPE FORMAT 15X NENIVELEE NORMALF FORMAT 15X COURBE DE CCRRECTICN OF NENIVELEE READS 12 84551 STATI FRAG 98 40 98 READS 13 ETAT BASSI STATI NETAL y JUNIV IF ETATL3G 97 32 UNIV 1 5 199 93 100 REAC S IPLMAX HPCL PEADISI2 IQATL L L KMAX LMAX L 5 2 HINTIL 1 READ 5 4 QINT L L 1 y KMAX p 95 Lal KMAX NQzQAIL 1 QA L DI FO U 1 HP L DOIN OINT L QAIL DINT HINT L HP L C l1 L DINT D9 DQIN D DI DINT DI DINITI 241 09 1 11 1 01 C 3 L 04 1 HRITEL6 16 RASSI STATI FRAG 18 10 4 999 702 703 HRITE16 24 0 999 WASTE 6 25 AGT 81
61. 4 15 8 65 300 5 18 8 65 153 6 19 8 65 110 7 28 8 65 302 8 29 8 65 255 9 31 8 65 328 10 7 7 66 242 11 9 7 66 110 12 12 7 66 183 13 17 7 66 288 14 19 7 66 123 15 3 9 66 290 16 23 7 67 101 17 4 8 67 273 5 18 25 8 67 287 5 19 8 9 67 106 20 90 4 2 M thode de l cart a la cote maximale TABLEAU IV STATION DE BOITIEK M thode du gradient limnigraphique AQ Q Qo AQ Keng AQ Q m s m s m s Si m s m s 0 200 0 390 0 190 48 7 30 22 0 26 7 0 104 0 286 5 000 6 300 1 300 20 6 3 2 2 22 8 1 140 5 160 0 900 1 085 0 185 17 0 23 16 8 0 2 0 1 085 9 200 9 800 0 600 61 8 59 D 0 020 9 780 0 400 0 515 0 115 22 4 31 22 6 0 2 0 001 0 516 0 140 0 144 0 004 2 8 2 9 0 1 0 0 144 8 300 10 200 1 900 18 6 18 13 2 54 0 450 9 750 3 000 4 000 1 000 25 0 3 2 2 22 8 0 685 3 315 16 000 18 100 2 100 11 6 16 11 7 0 1 0 015 18 115 3 000 2 850 0 150 5 3 8 5 9 0 6 0 020 2 830 0 125 0 144 009 132 4 29 103 0013 0 131 1 250 0 800 0 450 56 3 77 56 3 0 0 0 800 8 200 8 100 01400 1 2 32 234 2922 1820 6 280 0 240 0 252 002 48 13 95 47 0001 0 263 7 500 8 400 0 900 10 7 18 13 2 2 5 0 190 8 590 0 046 0 068 0 022 32 4 2 1 5 30 9 0 014 0 054
62. 5 119 116 115 CONTINUE 116 IF FBAU 99 99 119 11 1 117 117 0 GO 130 119 LEKMAX 120 HeHAU HPIL L etteat Ze LI 11129 1305130 129 OMAX nFB J 130 CONTINUF RETURN END SURPOUTINE TRANB INTEGER STAT I9FCH IPENTF JT Nr y STATT y ISTAT LMAX OL 16 e KMAX A13515 1 843 15 3 15 003715 eHAUTA 6 HAUTB 6 DER 6 OMAX LOMAX KDHAX Ze DONP 16 KBMAX L8MAX HEI 16 s HAUTI 63 330 JHAX NHEHAUTA J 4FLOAT INECHI 1000 IF NH 208 208 219 209 DEI Un 60 rn 330 209 DH FLOAT 1000 DDN NH DN 215 1 2 KDMAX TFLODN PONP L 210 2109212 210 1 1 70 220 212 1 11715 219 219 21 CONTINUE 219 LeKDYAY 220 Qn DDN OnNPfL EE ORTH TER HAUB HAUTA 37 DN 315 1 2 TF HAUR HPB11 310 310 312 3 0 L L 1 60 320 3l2 IFLHAUR HPREKHAY 315 3194316 315 CONTINUE 316 49 49 319 317 317 317 OFA J 10 Tn 330 319 LzKBMAX 320 HB HAUR HPBIL 1 2013301351 0090 4 0229 3309 339 329 330 CONTINUF SURRNUTINE TARAR INTFGER STATI COMHOM JOFGH IPENT E Au Ne JMAX ISTAT e HEI L i e 143 151 l n143 151 CU 15 213 293 HAUTA 6 6 NEB 6 OMAX e LOMAX Ze DD
63. 581 239 225 19 15 855 17 26 8 58 405 1350 1125 225 20 0 176 225 2 5 30 1095 18 15 9 58 458 1665 1390 275 19 8 123 22 5 2 7 40 1350 1958 59 19 28 10 58 530 2005 1815 190 10 5 51 153 48 85 1730 20 13 11 58 548 2205 1975 230 11 7 33 99 18 35 2010 D 21 17 3 59 456 1130 1385 255 18 4 125 22 5 4 1 55 1440 C 1959 60 22 23 9 59 463 195 1420 485 342 573 110 22 5 117 165 1585 D 23 26 3 60 303 515 665 150 22 5 270 22 5 0 0 685 1961 62 24 23 8 61 346 1060 860 200 233 559 213 22 5 08 5 865 C 29 10 11 67 537 2120 1865 255 13 7 583 46 138 01 0 1865 30 1 12 67 563 2260 2130 130 6 1 20 6 0 0 1 0 2130 31 4 12 67 567 2290 2170 120 5 5 16 4 8 0 7 15 2185 1967 68 32 13 12 67 574 2330 2245 85 3 8 9 27 1 1 25 2270 33 20 1 68 571 2110 2210 100 4 5 12 36 09 20 2190 34 30 1 68 560 1980 2100 120 5 7 23 69 1 2 25 2125 35 18 12 68 535 1840 1855 15 08 538 3 0 9 01 0 1855 1968 69 36 2 1 69 525 1710 1780 70 3 9 13 3 9 0 0 1780 37 11 1 69 509 1640 1675 35 24 29 8 7 6 6 110 1785 Le tableau suivant donne les quations des courbes de correction de pente M thode Station de Station de Dir e e de l cart la cote maximale 12 75 Arc tg 0 706 17 8 Arc tg 1 74
64. 603 2 540 2 595 55 2 1 10 135 2 460 il 19 1 55 609 2 645 2 680 35 1 3 4 0 2 680 1955 56 12 23 11 55 588 2 400 2400 0 0 612 24 170 2230 13 8 12 55 604 2 545 2 610 65 2 5 8 130 2 480 1956 57 14 7 8 56 247 125 425 300 70 5 565 318 195 620 1957 58 15 30 7 57 237 465 380 85 22 4 611 374 10 370 16 4 8 58 342 1045 840 205 24 4 581 239 5 845 17 26 8 58 405 1350 1125 225 20 0 176 25 1 100 1958 59 18 15 9 58 458 1 665 1 390 275 19 8 123 5 1385 19 28 10 58 530 2 005 1 815 190 10 5 51 40 1775 20 13 11 58 548 2 205 1 975 230 11 7 33 45 2 020 21 17 3 59 456 1 130 1 385 255 18 4 125 25 1 410 1959 60 22 23 9 59 463 1 905 1 420 485 34 2 573 110 4 8 210 1 630 23 26 3 60 303 515 665 150 22 5 270 1 8 10 675 1961 62 24 23 8 61 346 1 060 860 200 23 3 559 213 0 0 860 29 10 11 67 537 2 120 1 865 255 13 7 583 46 17 30 1 895 30 1 12 67 563 2 260 2 130 130 6 1 20 0 1 0 2130 1967 68 31 4 12 67 567 2290 2170 120 5 5 16 0 6 15 2185 32 13 12 67 574 2330 2245 85 3 8 9 1 0 20 2 265 33 20 1 68 571 2110 2210 100 4 5 12 0 8 20 2190 34 30 1 68 560 1980 2100 120 5 7 23 11 25 2 125 35 18 12 68 535 1 840 1 855 15 0 8 538 3 0 1 0 1 855 1968 69 36 2 1 69 525 1 710 1 780 70 3 9 13 0 1 0 1 780
65. 8 1 chelle Variable POH 302 Gradient 329 Mixte POH 302 limnigraphique 327 329 IUNIV 3 D importance faible Fixe 302 POH 328 moyenne pente 1 chelle Variable POH 302 Gradient POH 329 faible et relative Mixte POH 302 limnigraphique 327 329 vement tr s variable IUNIV 3 A B s e Importants A r gime r gularis par d bor 1 chelle Fixe POH 302 Ecart POH 328 dements Pente faible IUNIV 1 au maximum et relativement tr s variable 5 2 Programme de tarage En avant propos on indique les modifications mineures apport es au P H 302 des stations deux chelles pour le tarage des stations non univoques une chelle 311 5 2 1 MODIFICATIONS DU 302 est cr une variable caract ristique de la non univocit et qui peut prendre les valeurs suivantes 0 pour une station univoque 1 pour une station non univoque trait e par la m thode de P cart au maximum 3 pour une station non univoque trait e par la m thode du gradient limnigra phique 4 5 et 6 pour une station deux chelles cf 5 2 2 Cette variable TUNIV est introduite valeur 1 ou 3 sur la carte HP des hau teurs limites de tron ons de paraboles en colonne 80 Le nouveau format de lecture est 12 16 F 4 2 12 X 12 Cette variable UNIV est alors perfor e en sortie dans la m me colonne 80 de la carte mod le 306 portant les hauteurs limites in
66. ADB La sous routine TRADA identique la sous routine TRAD des programmes 301 ou 310 calcule le d bit correspondant HAUTA D l aide de la courbe enveloppe d talonnage sans influence de pente La sous routine TRADB calcule le d bit l aide des deux courbes de tarage d nivel e normale et de correction de d nivel e On y calcule d abord 1 d nivel e exacte entre les deux hauteurs DH HAUTA 1 HAUTB J FLOAT IDECH 1000 Si la d nivel e est nulle ou n gative le d bit est fait gal z ro Ceci correspond la non prise en compte des inversions d coulement biefs mar e ou autres quand elles sont rares et n gligeables cf 2 3 1 Si la d nivel e DH est positive on calcule nsuite la d nivel e normale DN FLOAT IPENTE 1000 puis le rapport de d nivel es DDN DH DN Les divisions 1000 s expliquent du fait des unit s des calculs m tres et m s et de celles de lecture de IDECH et IPENTE millim tre L introduction de DDN dans la courbe de correction de d nivel e permet le calcul de la correction de d bit QQN On calcule ensuite le d bit QN correspond la hauteur HAUTA J sur la courbe d talonnage d nivel e normale En d finitive le d bit est gal DEB J QQN ON La suite du traitement des d bits est absolument identique celle du 310 Apr s le traitement complet d une ann e de relev s d clench
67. EU 1 MIN 6 IF ITP 38 35 38 35 4 MZ 1 IMZ 240 MZ G T 37 38 IT IMZ HG IT HAU D 37 C NTINUE Un tel bloc d instructions demande une entorse aux r gles de perforation pour le 1er relev du 1 jour des cartes de relev s ant rieurs En effet s il correspondait 0 heure on aurait 0 et le test conduirait IT 240 On vite cela en don nant donc IT la valeur 1 ce qui s obtient ais ment en consid rant que le ler relev a eu lieu 0 heure 6 minutes L alt ration est minime et l incidence n gligeable La sous routine CORDEB est appel e aussit t apr s la sous routine TRAD qui effectue simplement comme dans tous les programmes de traitement des hauteurs le calcul des d bits pour une carte La sous routine CORDEB corrige tous les d bits de la m me carte Apr s introduction de HAUT dans HG selon un bloc d instructions d j d crit ci dessus on y effectue la recherche du gradient limnigraphique DH le calcul de la correction de d bit CORQ la correction du d bit la mise z ro des d bits corrig s n gatifs et le tri du d bit maximal comme dans le 328 La seule diff rence provient des contraintes dues aux r gles de perforation des relev s qui influent sur le changement d origine du vecteur HG et l effacement de portions de vecteurs devenues inutiles Ces deux op rations doivent tre faites lorsque le vecteur HG est rempli exactem
68. Fig 26 Carte COH 513 Bien que cette carte ait t imprim e dans la s rie 500 relative au traitement des BRE elle est utilisable pour n importe quelle station puisque les 8 premi res colonnes concernent galement l identification de la station comme dans toutes les cartes de 1 s rie 300 la seule diff rence des impressions de haut de carte qui ne correspondent plus ETAT BASSIN et STATION Mais le champ tant le m me ces 3 variables peuvent y figurer sans probl me La variable UNIV est perfor e dans a colonne 12 Rep re de Courbe Les autres champs sont utilis s dans les programmes de traduction en d bits et seront explicit s alors cf POH 312 328 et 329 Les sorties sont identiques celles du programme 302 pour les deux pre mi res courbes dont le titre Courbe enveloppe ou D nivel e normale est imprim La sortie relative la troisi me courbe bien que peu diff rente est donn e sur la figure 27 titre d exemple I Comme dans le 302 modifi cf 5 2 1 la variable IUNIV est perfor e dans la colonne 80 des cartes COH 506 en sortie du POH 311 5 3 Programme de calcul Th oriquement une station deux chelles pourrait concerner une rivi re dont des d bits d une station les relev s de hauteurs d eau seraient susceptibles d tre d pouill s selon l un quel deux chelles 312 des trois pas de temps En fait c
69. GRAD AC 080 IT GROUP HOIS IBIS LA Jx DIMENSION MANTIL7 1 177 213661 9 123 y LCOMDC12 DIMENSIONQJO 16 183 0 122 8 122 02131 12 1000 FORMAT 212 14 1844 1001 FORMAT 4X 14 3X 12 16F4 2 1002 4 3 12 868 5713 868 5 1003 FORMAT 1X ERREUR DE STATION 212 14 312 16 4 2 211 FORMAT 212 14 212 11 16413 111 FORMAT 212 14 13 23013 110 212 16 212 170 13111 212 FORMAT 1 2 FORMAT 1X ERREUR OYANHEE 1500 ed STATION NUMERD 2X 212 147 735X DEBIYS HOYENS JOURNALIERS 14 12 1501 FORMAT I X JANVIER 2 FEVRIER HARS 14X AVRIL 55 HAI ee Luut 3X JUILLET 4 4 0 2 5 1 0 2 0 2MBRE 1Xs DECEMBRE 1511 FORMAT 1 12 2 1229 2 1512 1 0 12 2 12 9 2 1513 FORMAT 1X 12 2X F9 2 9X 10F9 21 1514 0 12 2 3 9 2 9 2 9 2 2 9 4 2 1520 FORMAT 5X 12F9 2 1 1550 FORMAT 38X DEBIT HOYEN 1 763 9 291 4 3 547 37 40 LIMAL 763 29 2 1 3 5 Nm my xxx LECTURE DE LA CARTE OYTDENTIFICATION DES HAUTEURS LIMITES INFFRIEURES ET DE CELLE DES COEFFICIENTS D ETALONNAGE 99 READ S LOOD ETA
70. IEK PAR LA METHODE THEORIQUE Ah Date h At Q H 9 06 Y Qo Uo de jaugeage m cm cm h m s cm s m s cm s m3 s cm s 1 1 65 13 1 139 5 0 2 1 53 0 58 0 345 2 64 0 58 0 345 2 64 29 7 65 32 4 275 3 5 0 15 45 0 85 5 9 18 2 0 84 5 95 18 4 31 7 65 22 6 209 5 0 9 3 98 0 58 1 55 6 9 0 58 1 55 6 9 15 8 65 36 4 300 2 9 2 25 25 0 94 9 8 26 9 0 94 9 8 26 9 18 8 65 14 9 153 4 0 4 2 69 0 58 0 69 4 65 0 58 0 69 4 65 19 8 65 9 5 110 0 0 14 1 47 1 47 1 47 28 8 65 36 8 302 A 8 3 22 55 0 87 9 55 25 9 0 857 9 7 26 3 29 8 65 29 3 255 4 3 0 10 25 0 67 4 5 15 3 0 643 4 65 16 0 31 8 65 41 1 328 1 16 0 38 95 0 98 16 3 39 8 0 98 16 3 38 95 7 71 66 27 3 242 2 3 0 11 00 0 865 3 46 12 7 0 854 3 5 12 9 9 7 66 9 5 110 0 0 125 1 31 1 31 1 31 12 7 66 18 9 183 3 1 25 6 62 1 28 0 98 5 17 1 30 0 96 5 10 17 7 66 34 5 288 2 8 2 23 80 1 06 7 75 22 5 1 06 7 75 22 5 19 7 66 11 1 123 0 0 24 2 16 2 16 2 16 3 9 66 34 8 290 3 7 5 21 55 0 90 8 35 24 0 0 89 8 4 24 2 23 7 67 8 4 101 0 0 046 0 55 0 55 0 55 4 8 67 32 2 273 5 0 6 0 18 65 18 65 18 65 25 8 67 34 4 287 5 0 5 7 8 22 70 0 985 7 9 23 0 0 985 7 9 23 0 8 9 67 9 0 106 0 0 11 1 25 1 25 1 25 d bit nul 7 1 90 0 0 0 0 RK 1 4 Ah re 17 approximation 1 1 5 Ah 2 approximation Y 1 Pour une r solution simple on fait appara tre cart du d bit observ au d bit univoque ce qui permet de transf
71. IUNIV 1 ou par le gradient limnigraphique 1UNIV 3 et condition que ses relev s puissent tre seulement au maximum de 16 par jour ou de 8 par jour pour la m thode du gradient limnigraphique qui ne traite pas le CAS 9 c est dire quand il y a de 9 16 relev s par jour Pratiquement le 328 s applique donc la fois aux cours d eau r guliers importants et faible pente grands fleuves type Niger S n gal et aux cours d eau d importance moyenne pente relativement variable La restriction d emploi de la m thode du gradient limnigraphique n est pas grave car cette m thode est surtout celle des petits cours d eau justiciables du pas de temps variable pour le d pouillement des relev s et pour lesquels existe le 329 cf 5 5 Par cons quent les cours d eau rapides dont la non univocit du tarage doit tre corrig e par la m thode du gradient limnigraphique et dont les relev s ne peuvent pas tre pris moins de 8 fois par jour doivent tre d pouill s selon le pas de temps variable relev s limnigraphiques int graux carte 311 et leurs d bits calcul s l aide du 329 328 une structure identique celle du 301 partir du traitement d une ann e enti re apr s lecture de toutes les cartes de hauteurs de celle ci et traductions en d bits Les sorties imprim es et perfor es sont galement identiques Les parties diff rentes
72. NP 161 KRMAXs LOMAX 16 HAUTIS 1001 FORMAY amp X 13 4X 12 16F4 2 1002 4 13 4 12 8 8 5 13 8 8 5 1 1003 FORMAT 1X ERREUR NE STATION 1006 FOnMAT 4X 13 6X 12 16F4 3 1907 8 413 1 13 1 214 5 INALISTATI LMAX HP IL 151 LMAX IFISTATI ISTAT 1 27 WRITE 16 1003 KMAX LMAX 1 00 5 7 1 3 READ 5 1002 5TATIy Ir Cl Sel rl 19 K MAX IFISTATI ISTAT 3 4 3 MA ITE 6 1003 NO 6 L 1 KMAX Siet ist At REAN S 1002 HPB L Lei LBMAXD 45 1 15 7 8 7 6 1003 KRHAX LBHAX 1 DO 12 J 1 3 READ 5 1002 SYATI Is C J L L l KBHAX IFISTAYI ISTAJ 10 11 10 lv WPITE 6 1003 11 12 L 1 KBMAX 12 8 1 1 0 1 READ S 1006 STATI LDHAX pnNPIL1 L lyLDHAX IFISTATI ISTAT 13 l4 LA 13 6 1003 14 KDMAXeLDMAX l PO 18 J 1r3 541002 5 1 161 4241 LE Le 1 lt 15 47316 17 16 16 WRITE 6 1003 17 nn 18 18 DffeL C JL PEAD S 1007 STATI4NETALs LOECH LPENTE RETURN Nm Fur SUBRDUTINF YMPRILIAN K JLI K LM INTEGER 2 12 31 51 12 531 12 14 12 0 6 DIMFHSION 14612 DATA 14 2 43 644 456 6 47 48 49 Or to 1 1 3 1211 541 14 3 FDPMAT 6X 6 JAN FFV MAR AVR Jus JUL aqu SEP DEI NOY 4 4 HOY p121541 1X WRITE 6 3 DQ 29 4 1 32 00 27 M 1 12 JM 5e M 1 1 29317 13
73. Qo Qo Boitiek 0 73 Qo TABLEAU II STATION DE MOPTI M thode du gradient limnigraphique AQ H 27 Sens N Date 5 gt K Qo E Qo cm m3 s m s m s VA m s m s D 3 26 1 52 527 1 570 1 855 285 15 3 30 11 2 4 1 75 1 780 D 4 27 3 52 232 280 295 15 51 20 7 4 2 3 7 288 C 5 8 6 52 126 130 106 24 22 6 19 7 1 15 5 16 122 D 6 19 2 54 310 505 560 55 9 8 24 8 9 0 9 5 555 7 23 3 56 230 290 290 0 0 16 6 0 6 0 17 307 9 18 12 56 542 1770 1980 210 10 6 30 11 2 0 6 10 1990 10 18 2 57 233 265 297 32 10 8 18 6 7 4 1 12 285 il 8 5 57 100 74 67 7 10 4 2 0 7 11 1 7 74 12 6 8 57 473 1 600 1 485 115 7 8 36 13 4 5 6 85 1 400 D 14 16 1 58 514 1 570 1 760 190 10 8 33 12 3 1 5 25 1 785 15 19 2 58 346 530 760 230 30 3 25 9 3 21 0 160 600 D 16 12 3 58 268 280 380 100 26 3 17 6 3 20 0 75 305 17 4 6 58 227 360 284 76 26 8 72 26 8 0 0 284 18 25 7 58 414 1225 1140 85 7 5 15 5 6 1 9 20 1 160 19 13 8 58 464 1555 1435 120 8 3 22 8 2 0 1 0 1 435 20 6 9 58 580 2 400 2275 125 5 5 25 9 3 3 8 85 2190 22 22 8 59 815 2 000 1765 235 13 3 42 15 6 2 3 40 1725 D 26 14 3 60 142 120 131 11 84
74. T 720 D BIT MAXIMAL 5 6 11 2000 0 5 2 JBAS 3 ROSTA JANN HOES JDUR I HEULI 11515 6 JUNITE ICHANG NOCART SF NETA 1111 39 1111 12 1111 MRITE66 15 IBAS ISTAT FRAG TF JUDEB1 8U51 8055 8051 1 8055 MRITEL6 7001 40 60708060 8051 MRITE 6 7002 8060 QMAX 0 30 GO TO 111 100 END SUBROUTINE TARAG INTEGERSTATI COMMON JHIH 1 5 LMAX 1 161 KHAX 3 15 13 15 146 8 6 QMAXs IT 018001 e SGRAD LUNIVS HZ INZ HEW GI HINTS 2 NHBr ALPHAL 5 5 HSI SY 1001 44 14 3 12416 4 2 2 11 1 43 12 8 8 5 13 8 8 5 2 ERREUR DE 5 1 1 3 1003 FORMAT 1X DATA H 0 0H1 0 00 100 J 1 IFIN ITP 10 HEU J RINtJ 6 1F 11P141 39 41 FZ HZ 1 EHZ 240 H2 IT ITP IHZ HG iT HAUT J IT1 IT IYG IT1 JGRAD 1 I0EB DEB J 100 IF 10EB 1 1 2 0 GO TO 15 15 115 1 DO 20 Ll 2 16 12115 11 1F HG1LZ 9 989109 109 3 IF HG1L2 99 9814 10 10 CONTINUE 5 OHISHAUT J HG LZI 14131431445 20 1 176 171 1 61 19 9 98 20 20 6 12 99 98 7 20 20 CONTINUE DH2 HAUT J HGCHL I H2 M1 ITG DH DH1 H24DH2 L1 1 L1 11 M2 60 TO 8 IF NHR113 16 13 DO 17 1 1 NHB IF HAUT J HB 1 17 21 21 CONTINUE U ALPHACTF HAUT BETA I IFIDH LT 0 COR Q2 AC SQRT ABS DH ZU IF DH GE 0
75. T BASSI ISTAT FRAG IFLETATI98 999 98 88 IBAS BASSI DO 2004 I L 300 2004 1 9 9 CALLTARAG LECTURE DE CARTES SPECIALES POUR STATIONS NON UNIVOQUES IF IUNIV 3 111 63 63 63 READ S LIETAT BASSI STATI ANNEE HOES GROUP HAUT 9 CAS CHANG IAN ANNEE IF EJAT 106 100 62 62 IF STATI ISTAT 60 61 60 60 MRITE 67 1004 1004 FORMATI1X ERREUR DE STATION POUR LE GRADIENT GO TO 999 61 IT x 0 65 00 46 J 1 16 TaIT 24 CAS 46 HGUITI2HAUT JI GO TQ 63 111 READ 5y8 STATI LAN HMAXO HMAXL 2 8 FORHAT 4X 014512036402 3103 IFISTATI ISTAT 1104 110111104 1104 WRITE 6 1103 1103 FORMATILX ERREUR DE TARAGE DE STATION TO 105 1101 KK 0 1116 1117 1116 1116 KK KK 1 60 JO 100 1117 KK KK 2 LECTURE 025 16 HAUTEURS CONTENUES DANS UNF CARTE 100 Ven SHAX 0 101 READ S 91 DETAT BASSI STATT ANNER 8015 GROUP HAUT CHANG 112011201106 1120 IF TUNIV 3 1121 68 68 1121 IFULAN ANMEE 1105 1106 1105 1105 HRITE 6 1107 1107 FORRATfLX ERREUR TARAGE 60 TO 105 68 IF TAN ANNEE 110667 67 87 HRITEL6r 1205 1005 FORMAT 1X ERREUR DMANNEE POUR LE GRADIENT 1 GO TO 999 gt 1106 IF ETAT 103 10 103 4 103 Ku MDLS IGRO GROUP K K 1 IFUSTATI ISTATI 1041065104 10 WRITLI67 10033 105 REAO S LIETAT 1 105 102 105 106 1 107 139 107 107 108 139 108 108 WRITE 6 2 6070105 139 HILLSANNEE
76. T de la pente d coulement Si l on appelle U la vitesse moyenne d coulement en r gime non uniforme et pour la hauteur h l chelle et si l on n glige l accroissement de U pour l accroissement Ah de h l accroissement de pente est gal CA U t repr sentant la distance parcourue par l onde de crue pendant le temps t On crit alors I l AI 1 D Al U At L quation d riv e de CHEZY pour le r gime non uniforme et la hauteur h s crit 1 2 0 To U At 13 14 U En posant y er y toujours positif on peut crire l quation pr c dente de la mani re suivante apr s avoir remplac 0 10 47 1 Ah 3 5 3 PRE 6 1 La d riv e de l quation 3 s annule pour 0 578 valeur minimale pour y laquelle correspond un gradient maximal en valeur absolue en d crue de 0 386 lo Uo est int ressant de noter que sur deux exemples analys s la Kayanga en Casamance et Boitiek en Mauritanie les gradients exp ri mentaux ne d passaient pas ce maximum th orique Sachant que y et pour faciliter les calculs l quation 3 peut s crire 1 Ah o At SS 4 En supposant Io ind pendant de h on peut r soudre cette quation en partant d une valeur choisie a priori pour Io et par it ration jusqu ce que la courbe de tarage univoque Qo H pr sente
77. Tarage et calcul des d bits des stations hydrom triques non univoques P DUBREUIL Directeur de recherches Chef du d partement de la Recherche appliqu e au Service hydrologique de 5 J P LAMAGAT et G VUILLAUME Charg s de recherches Hydrologues au Centre O R S T O M de Lom Togo Beaucoup de cours d eau ont un r gime d coulement non uniforme permanent dit graduellement vari si un bief d un cours d eau de ce type a un fond stable et un contr le g n ralement aval invariable toute station hydrom trique implant e dans ce bief est univoque c est dire qu une hauteur d eau donn e correspond toujours le m me d bit f Certains cours d eau ont par contre soit un contr le aval variable soit un r gime d coulement non permanent c est dire rapidement variable M me si un bief d un tel cours d eau a un fond stable la relation hauteur d bit n est univoque en aucun ses points c est dire que pour une hauteur donn e on peut avoir au moins deux valeurs et m me une infinit de valeurs de d bits dans une certaine gamme autour du d bit correspondant au r gime permanent On parle alors de station hydrom trique non univoque Le tarage d une telle station est compliqu il exige d une part un nombre de mesures de d bits bien plus grand que pour une station univoque et d autre part la courbe de tarage univoque c de la place un faisceau de plusieurs courbes
78. UL Io Lolo KMAX ELC L 1 KMAX POQKP QJNP L 1 928P L mp DnNeiL 1i RnqetL SQ9NIK QQNINIIL QQNPIL J JOINT DDNINT L I PORP L CU Le L 00INT DAONP DOQNINEND INT DO DD INTI 2 NOE LLION 2 L CCNPIL WRITELE 1 RASSE STATE FRAG 992 WRITE S 1A WRITE 6 261 301 WRIIF 6 6 WRITE 667 5 5200 CIA leie Ze LE La KMAXD LES 223 CL NP EL Jetel LMAXT TMAR 10 DORP LIAXI DDNP 1 1 0 Dt1 DONPI1 Q 1 00NP 1 MA 700 1 2 1 DUIISN I 11 0 1 NTsD 1 IFLCT NONP KMAX1 7208 710 719 710 LEKYAX 750 720 147 L 2rKMAX IFLOT DONP L 22730 730 740 132 15151 750 74 CONTINUE 750 lt 1 0 1 Q I X C 1 L X C 2 L C 3 L 1F 911 799 700 700 799 I O 700 CONTINUE W9JTF 6 21 MRITEL6 73 nn 850 I 1 14a 1 10 DO 5104 J 1 3 1 7 10 TFIat1 T159509510495104 5104 CONTINUE 5050 001 1 1 912 10 0 5 10 0 60 70 5800 5801 5802 55890 HRITELE 2I0CINL TI OLI TO 859 SAOI WRITE 6 2702 01 001 67 10 850 5802 MRITEL6 2002 011 011 40 CONTINUE S Da 760 1 1 1 LL 1 762 162 7632 763 KANNPILI O 60 Zen 762 41000 169 CONTINUF LIMAX LMAY L 90 764 I LIPAX 16 764 KDNHP I 9 URITELT 14 ETAT BASSI STATI KADNP CL Lie 16 IUNIV 992 03 145 L 19KNAX DS 165 VAR CLKAL I 120 VARIA ABS VAR 17 IF VARIA
79. UTINE TARAG INTEGERSTATI COMKONJKIH JHAX HAUT 161 16 16 3 15 I COMMONSTATI ISTAT 3 15 HMAX Le HMAX2 9 KK 9 AAs TUNIV COMMONHG 300 CAS ICA JGRAD IT GROUP MOIS IBIS LA JX 4 14 3 12 16 4 2 2 11 14 14 3 12 8 8 5 13 8 8 53 FCRHAT 1 ERREUR DE STATION FORMAT 1X ERREUR D STATION POUR COXRECTION TARAGE 4 14 2 11 8 265 2 10 14 31 FORMAT 1X ERREUR DE METHODE DE TARAGE READIS 1001 STATI LMAX HPIL L 1 16 IUNIV IFISTATE ISTATH 1252 BRITE 6r 1003 KMAX LMAX 1 DO6J 1 3 5 1002 5 1 1 2 1 L l KHAX IFCSTATI ISTATI 314 3 MRITEL4 19031 DEGL 1 KMAX ACI LI CII AL READ S 1012 STATE JUNIVe AA 9 AB IF STATI ISTAT 7 591 HRITE 6 1004 TFLIUNIV IUNIVI 9 1019 9 6 1013 RETURN END KXXX XXXX XXXX KXXK XXXXXX XXXX XXXX KXXX XXXXX KXXX 69 70 110 112 115 115 117 115 120 1118 1119 1120 1121 1122 129 130 SUBROUTINE TRAD INTEGERSTATI COMMONJMIN JMAX HAUT 162 LMAX HP Lei DEBI 16 4413415 1 COMHONST ATI gt ISTAT C 13 3 159 HNAXOr HMAX y HMAXZ e KK 9 AA AR y TUNIV COMMDNHG 300
80. X QUAX 0005 HAMAXS 10009 QA x 1CMAX 111 WRITELY 12 KRAS MAMOT MDIS s ICHOL HOTS MOTS Te 112 JUDERI JFINTYP 1 85 686 690 685 HRITE16 625 ETA BASHYD BRE ca TN 695 690 HRTTEL6 635 1FTA BASHYD BRE 695 HAITE 6 5552 KFTA KBAS TDVREP NOMBY 696 IF JUDERI 4910 4900 4910 4900 MRITF 63 5056 MILL GATI 5021 4919 NOITE A SOSTIMILL 5021 CALL YAPRT MILL MAJ ICJ MANDI y ICHDI XMAX 0 0001 FLOATIMAMAX 4RFLOAT ICMAX 1 XMAX XxMAX 0 01 HRITEL6 5100 5109 FORPMAT T20 0F8JT MAXIMAL INSTANTANE 3 F9 1 2000 5 7 JAAS RASTA JANN MAIS JOUR TD HEI 13 BINEI 18 1141 8 JUNJTE ICHANC NOCAR 1111 99 1111 1111 1 1 amp 40 640 652 640 6 525 60 70 661 652 HPITF 6 635 FTARASHYONBRE 661 16 5552 5 8 IF JUNES1 3051 80 5 4951 8055 HPITE 6 7001 60708060 8051 6 7002 8060 GMAXs0 GQ 40 9999 5700 END SUBROUTINE TRADA INTEGER STATI COMMON INFCHISPENTE JM Ny 6 ISTAT HP1 16 AL 3r 1S 13 15 3 15 003 15 HAUTA 6 y SI y QMAXy 2 008 14 KBMAX LBMAY 16 HAUT 61 00 130 J JMIN HAU BANTA 4 115 L 2yKMAX IF IHAU HP L 3110 110 112 110 19 1 TO 120 112 1 11
81. a station qui a servi de mod le cette tude celle de la Crique Virgile en Guyane Il y avait eu pendant moins de 48 heures conjonction d une crue forte sur l Orapu dans lequel se jette la Crique Virgile et dont le remous influence la pente d coulement et d une forte mar e remontant dans le bief de l Orapu 2 3 2 LA CESSATION D EFFET DE LA PENTE La cessation d effet de la pente est possible dans certains cours d eau qui sans une forte modification du contr le aval auraient une pente d coulement assez forte et peu variable pour que le r gime soit alors consid r comme permanent graduelle ment vari et le tarage comme univoque L allure exponentielle croissante de la courbe de correction de pente indique d ailleurs bien que pour de fortes d nivel es Di le rapport Qi Qn tend asymptoti quement vers une limite Crique Virgile nous a encore servi d indicateur de limite d effet de la pente et son tude a permis de faire un choix simple de celle ci On consid re que la pente il QU eri e ib J 1 1 1 Lem Jsugeage D Jaugeage corrig pour n Hi Fig 5 Trac de 1 courbe de tarage d nivel e normale X Ji Jaugesge NIT 101 011 1 1 1 1 1 1 1 Trac de la courbe de correction de pente Fig 7 Courbe enveloppe de tarage sans influence de la pente 12 devient suffisamment lev e pou
82. ai de classification des tarages non univoques pentes et vitesses peu variables pour une cote donn e On peut alors parler de pseudo univocit de crue et de d crue et le tarage comporte une courbe pseudo univoque de crue une courbe pseudo univoque de d crue et une courbe de raccor dement variable elle chaque ann e en fonction de 1 cote maximale Selon que les maximums de hauteur et de pente donc de d bit co ncident ou non le raccordement des courbes de tarage pseudo univoques de crue et de d crue se fait soit en point de rebroussement Dir et Tonka sur le Niger fig 2 soit en courbe douce Mopti sur le Niger Matam Podor Dagana sur le S n gal fig 3 Dans le second groupe de cours d eau la non univocit est due en fait un r gime h t rog ne de l coulement c est dire que les crues d origine diff rente ont des hydrogrammes diff rents et des conditions d coulement pentes et vitesses diff rentes pour une m me cote Si le cours d eau a une faible pente et draine un petit bassin on peut avoir une non univocit g n ralis e toutes les crues tant plus ou moins diff rentes et le sch ma de tarage est celui de la surface de tarage fig 1 mais avec une dispersion peut tre moindre qu avec un contr le variable C est le cas de 1 Oued Boitiek en Mauritanie ou de la Kayanga en Casamance Si le cours d eau est important et de pente non n gligeable on observe un r gime principal d co
83. ands fleuves r guliers auxquels l application de la m thode de l cart au maximum puisse tre conseill e comme donnant de meilleurs r sultats et comme tant plus simple La m thode du gradient limnigraphique offre un inconv nient et un avantage L inconv nient du gradient limnigraphique est tel qu il conduit un hydrogramme en dents de scie autour de l hydrogramme r el inconnu Les figures 23 et 24 illustrant les exemples pr sent s dans le chapitre 4 2 3 sont particuli rement nettes cet gard Cet hydrogramme en dents de scie n offre pas d inconv nient au niveau du calcul du d bit moyen mensuel et m me de celui du volume d une crue les carts de part et d autre de l hydrogramme r el devant se compenser Par contre ces dents de scie sont g nantes pour d un hydrogramme de crue donn e surtout si cette analyse porte sur la forme et les param tres caract ris tiques de l hydrogramme En effet une d formation possible est le d placement du maximum ce qui va perturber le calcul du temps de mont e et des param tres de l hydrogramme standard par exemple Or un lissage manuel peut seulement corriger la forme de mais pas la position du maximum On entend par lis sage le trac de r el vue aux moindres carts au milieu des dents de scie de l hydrogramme brut calcul La solution pour obtenir un lissage parfait et rapide d une
84. ation de l ensemble BRE 501 sert de test pour la fin de traitement des donn es d un bassin c est dire que m me si l on traite successi vement les donn es de deux bassins d un m me ensemble on doit replacer une carte 501 avant les donn es du 26 bassin La carte COH 501 sert aussi contr ler la fin du travail La seconde carte d identification 502 relative au bassin trait suit toujours la carte 501 Elle sert tester que les donn es d talonnage et les donn es de hauteurs d eau appartiennent bien au bassin Par exemple 6002 IBVREP NGSTA 112 115 112 IBVREP est le num ro d identification du bassin lu en 14 sur la carte 502 et NOSTA est le m me num ro lu sur une carte de hauteurs dat es 511 L instruction 112 d livre un message d erreur de station L instruction 100 CALL TARAG suit Dans cette sous routine on lit les param tres des 3 courbes d talonnage et on les inscrit sous diff rents noms de variables LMAX et LL pour la courbe enveloppe b LBMAX et LL pour la courbe d nivel e normale LDMAX et D LL pour la courbe de correction de d nivel e Apr s la lecture de la premi re carte de relev s COH 511 on teste le num ro de la carte S il est sup rieur 5000 est en tat non univoque il faut donc lire une seconde carte 511 pour avoir l information hauteur HAUTA 1 dans cas Lecture carte
85. dient post rieur DH 2 Un poids proportionnel aux distances L 1 1 et M 2 de ces cotes est attribu chaque gradient pour le calcul de DH Et l on a DH DH 1 M 2 DH 2 L 1 D L 1 9 M 2 Comme le pas de temps du gradient JGRAD ne d passera jamais 8 jours pour tout cours d eau on jug inutile de dimensionner le vecteur HG au del de 300 heures afin d conomiser des positions de m moire En fin de CORDEB suivant la valeur de IT 1 on proc de ou non une modifi cation du vecteur HG suppression des valeurs les plus anciennes changement d origine de temps des valeurs conserv es effacement des anciennes valeurs d cal es On a choisi de proc der cette modification quand IT 1 d passe 270 heures et de ne garder que les 200 derni res heures Le changement d origine transporte les 200 derni res valeurs de HG conserv es aux positions de 1 200 il faut donc ensuite effacer les valeurs qui occupaient des positions entre 200 et IT 1 Tout ceci s effectue dans les boucles DG 6114 et DO 6112 Cette limitation du vecteur HG pr sente un inconv nient en cas de lacune d observation durant une p riode sup rieure 300 heures 12 jours le vecteur est vide de cotes r elles Dans cas la derni re valeur de DH 1 ou de DH est conserv e et utilis e jusqu ce qu une nouvelle puisse tre calcul e la premi re observation apr s la phase de lacune on DH 2 0 videmm
86. du PSH 301 sont les instructions de d but et de fin de traitement les blocs relatifs la traduction des d bits et leur correction Mais ces parties de programme nouvelles sont totalement distinctes selon qu il s agit de la m thode de l cart au maximum ou de celle du gradient limnigraphique L organigramme g n ral 328 fig 29 la logique d taill e du calcul du d bit moyen journalier fig 30 et les organigrammes des sous routines TRAD et CORDEB fig 31 permettent de comprendre les articulations des deux BEE au sein du 328 qui les r unit les tests successifs IF 3 venant aiguiller sur l une ou l autre des m thodes aux moments opportuns Il nous a paru plus clair de traiter chaque m thode s par ment plut t que de d crire en une seule fois le programme complet 47 48 Lecture carte identification station non Appel TARAG ecture des hauteurs maximales des ann e N 1 N N 1 de relev s ant rieurs aux donn es ces relev s dans le vecteur hauteurs temps HG IT erecarte 301 Lecture 1 relev s hauteurs non oui Calcul du d bit moyen mensuel annuel Interpolation fupression tous Appel TRAD pour traduction avant changement crue d crue Changement crue d crue KK 1 Lecture carte 301 SE Appel TRAD pour traduct
87. e les l ments de calcul pour les 19 jaugeages disponibles Deux approximations sur la valeur de Io ont suffi pour obtenir une bonne correction de pente comme le montre le trac de la courbe de tarage sur la figure 8 La somme des carts absolus cette courbe a t r duite de 75 ce qui est tr s satisfaisant 4 1 1 MODE D EMPLOI On peut simplifier la m thode it rative d velopp e au chapitre 3 2 La premi re simplification consiste consid rer comme constant le facteur Ab multipliant Ee dans l quation 5 Cette formule montre d ailleurs que la dispersion du d bit par rapport au d bit univoque diminue quand la vitesse cro t ce qui est compr hensible Les courbes traditionnelles de tarage crue d crue des grands fleuves type Niger ou S n gal dans leurs cours d bordements hors du lit majeur montrent un peu la m me chose l cart entre d bit en crue et d bit en d crue pour une m me Ah e hauteur d cro t quand la hauteur done U et cro t pour s annuler quand il a point de rebroussement Si la vitesse U varie peu pour la gamme des hauteurs h le facteur varie encore moins et l on peut admettre qu il est peu pr s constant et gal ak La seconde simplification consiste 2 suivante JE yo Ab L quation 5 avec ces deux simplifications peut s crire Q Q1 k 6 TABLEAU I CALCUL ITERATIF DU TARAGE DE L OUED BOIT
88. ectue partir d une quation Arc tangente de l cart entre la cote consid r e et la cote du maximum annuel de crue Le traitement des donn es s effectue par ann e civile or l ann e hydrologique ne co ncide pas toujours avec celle ci Par cons quent la connaissance de la cote maximale annuelle de l ann e N consid r e ne suffit pas il faut lui adjoindre soit a La cote du maximum de l ann e ant rieure N 1 si l ann e N commence par phase de d crue b La cote du maximum de l ann e suivante N 1 si l ann e N s ach ve par une phase de crue On voit tout de suite la restriction de cette m thode on ne peut parfois traduire en d bits les hauteurs d une ann e que si le maximum de l ann e suivante est connu L cart la cote maximale annuelle intervient en valeur alg brique dans le calcul de la correction de d bit c est dire qu il est positif en d crue et n gatif en crue cas g n ral cf 4 2 2 faut donc le maximum de r f rence pouvant galement varier entre le d but et la fin de l ann e disposer d un indicateur de pas sage par les cote maximales et minimales annuelles c est la variable IGRAD et d un compteur de position crue d crue permettant de choisir la bonne quation de correction c est la variable KK La figure 32 montre le m canisme d emploi de ces deux variables 1 Fig 32 M thode de l cart
89. ee Perforation impressio P H 310 carte releves 311 Lecture oui Fin de station Fig 34 Organigramme du 329 non 5 5 Programime de traitement des relev s limnigraphiques int graux d une station non univoque une chelle POH 329 Apr s la carte d identification et les 5 cartes de TARAG on trouve les cartes de relev s ant rieurs 1 2 ou 4 suivies d une carte blanche puis les cartes de relev s de la premi re ann e traiter toutes cartes sur le m me mod le 301 Une carte blanche s pare ces derni res de celles de la deuxi me ann e s il a lieu et pour laquelle il n y plus videmment de relev s ant rieurs Une deuxi me et une troisi me cartes blanches terminent une station et tout le travail On note que dans les deux m thodes la cinqui me carte de TARAG est une carte mod le COH 513 fig 26 qui contient les param tres des quations de correc tion du d bit Le format de lecture de cette carte est tabli de telle mani re que les champs affect s aux param tres des deux m thodes n interf rent pas espace vide de 10 colonnes AA et AB pour l quation Arc Tg de la m thode de l cart au maximum sont Jus en F 5 2 JGRAD et AC sont lus en I 4 et F 4 3 pour la m thode du gradient Le format complet comprend galement la lecture de STATI et de IUNIV pour les contr les 4X 14 3X Il 8X 2 F 5 2 10X I4
90. elui ci s ex cute sans difficult puisqu il n y a pas d ambiguit dans le choix de la d nivel e normale le syst me it ratif est plus complexe car il n cessite Ah a Une estimation du pas de temps At pour le calcul du gradient b Une estimation de Io pour le d part des calculs 3 D bit en m3 s O Jaugeage effectu Q H e Jaugesge corrige 1 2 Hauteur Fig 8 Tarage univoque par la m thode th orique de l oued Boitiek 15 4 Tarage des stations non univoques une chelle 4 1 M thode du gradient limnigraphique En outre ce syst me n cessite le trac de deux courbes l une du d bit en fonc tion de la hauteur l autre celle de la vitesse U en fonction de la hauteur Bien que la prise en compte de 1 vitesse puisse tre n glig e cf chap 4 la m thode it rative reste assez longue et son passage en calcul automatique peut galement n tre ni simple ni conomique On peut supposer que pour des cours d eau dont le caract re de non univocit n est pas trop accus les it rations seront limit es En tout tat de cause cette m thode est g n rale et applicable On en fait la d monstration sur un exemple pris parmi les petits cours d eau variations relatives importantes de la pente qui est faible en valeur absolue celui de l oued Boitiek du bassin du Ghorfa en Mauritanie Le tableau I rassembl
91. en quation de la courbe de correction de demivel e Impression des resultats Pertoration carte READ S5 1 LMAX HP L 99 Mise en quation des courbes enveloppe et CH 506 d nivel e normale Impression des resultats Perforation carte 506 DS 165 L 1 KMAX Calcul des coefficients des paraboles pour imprimer des cartes 507 en format sp cial 5 chiffres carac t ristiques Fig 25 Logique 311 Carte L d bits pour HP L format 16 E 5 3 ou carte QQNP L rapports de d bits pour DDNP L en format 16 F 4 3 Carte HINT hauteurs aux points interm diaires 15 F 4 2 ou carte DDNINT rapports de d nivel es aux points interm diaires 15 F 4 3 Carte QINT d bits aux points interm diaires 15 E 5 3 ou carte QQNINT rapports de d bits aux points interm diaires 15 F 4 3 Le bloc des 5 derni res cartes c est dire sans la carte d identification est r p t 3 fois autant qu il y a de courbes de tarage et de pr f rence dans l ordre croissant de 4 6 pour IUNIV Ensuite apr s cet ensemble de 1 5 5 5 16 cartes Carte blanche de fin d talonnage Carte blanche de fin de donn es On remarque que par diff rence avec le 302 la carte NETAL 513 pr c de et ne suit pas les 4 cartes relatives aux paraboles de chaque courbe et qu il faut une nouvelle carte 513 avant chaque nou
92. ent Donc DH est pris gal la valeur conserv e d 1 1 de DH 1 corrig e du rapport T 1 Ma qui est alors de l ordre de 0 40 0 55 Du fait de la r duction de DH 1 Perreur sur la correction de d bit ne peut d passer des proportions admissibles cette erreur diminue et s annule quand les observations nouvelles ont dur un temps sup rieur JGRAD puisque alors une vraie valeur de DH 1 est calculable Enfin bien entendu on n op re pas de correction si le d bit est n gatif ou nul et si apr s correction le d bit devient n gatif on lui donne conventionnellement une Valeur nulle La pr sentation des cartes de donn es pour la m thode du gradient limnigra phique est un peu diff rente de celle utilis e pour la m thode de l cart au maximum 57 58 Lecture carte identification station non Lecture etalonnage Mise zero de la base de temps MZ et IT des compteurs de jour de dixiemes d heuresITP IMZ ben Lecture de 2jours de releves D 37 I 16 t anterieurs sur carte 31 dades releves dans le vecteu hauteur_temps HG IT Lecture cartes releves o 311 CONTINUE oui anterjeurs Controle d erreur INJ BU 2 00 Traitement des donn es Traitement de la carte 311 Calcul des j limnigra H 310 Call CARDEB d une ann
93. ent jusqu un temps IT gal un multiple de 240 61 5 6 Remarques sur la m thode 62 du gradient limnigraphique dixi mes d heures et le d calage d origine doit tre galement un multiple de 240 On a adopt l effacement 720 et la conservation d un bloc de 480 dixi mes d heures De la sorte les compteurs MZ et IMZ cumul s en d but d op rations se retrouvent toujours apr s valoir exactement 2 et 480 respectivement Le d clenchement des op rations n a pas lieu au premier passage de par 720 mais au second c est dire non pas lorsque ITP 240 mais lorsqu il gale z ro aussit t apr s Un compteur provisoire N pourvoit cela De la sorte lorsque l on reprend le traitement d une nouvelle hauteur HAUT en d but de CORDEB il s agit toujours du premier relev d un jour suivant le relev de 0 heure et les valeurs de MZ et IMZ r g n r es auparavant ne sont pas modifi es puisque le premier ITP est diff rent z ro ce qui n aurait pas t le cas en effectuant les op rations au premier passage de IT par 720 correspondant ITP 240 Ces deux op rations s effectuent dans les boucles DG 40 J2 IT3 et D 31 ITN II La m thode du gradient limnigraphique appara t comme la m thode de correc tion de non univocit des stations dot es d une seule chelle dont le domaine d appli cation est presque g n ral En effet il n y a exception que pour les gr
94. es dents de scie dans l hydrogramme 27 TABLEAU VI STATION DE MOPTI M thode de l cart la cote maximale A H Q Gs Hx Ca Q Sens Ann e N Date 9 9 cm m s m s m3 s Sa cm cm m3 s m3 s D 1951 52 3 26 1 52 527 1570 1820 250 13 7 712 185 117 20 35 1785 4 27 3 52 232 280 325 45 13 9 480 16 3 2 4 10 335 1952 53 3 8 6 52 126 130 111 19 171 719 593 169 02 Z 0 11 1953 54 6 19 2 54 310 505 595 90 151 731 421 159 08 5 600 D 1955 56 7 23 3 56 230 290 320 30 94 731 501 165 69 20 340 D 8 9 18 12 56 542 1770 1945 175 907 686 144 10 1 1 1 20 1965 1956 57 10 18 2 57 233 265 325 60 18 5 453 161 24 10 315 11 8 5 57 100 74 74 0 0 586 168 168 12 86 12 6 8 57 473 1 600 1 425 175 12 3 257 13 6 1 3 20 1 405 D 14 16 1 58 514 1570 1720 150 87 730 216 126 3 9 65 1785 1957 58 15 19 2 58 346 530 765 235 30 7 384 15 5 15 2 115 655 16 12 3 58 268 280 430 150 34 9 462 162 187 80 350 C 17 4 6 58 227 360 310 50 16 2 480 163 01 0 310 1958 59 18 25 7 58 414 1 225 1 095 135 12 3 707 293 142 19 20 1075 19 13 8 58 464 1555 1370 185 13 5 243 13 3 02 5 1375 20 6 9 58 580 2400 2250 150 6 7 127 9 3 26 60 219
95. est sur la carte d identification de station Le traitement syst matique des hauteurs de grands fleuves comme le Niger et le S n gal n cessite l introduction d une convention suppl mentaire Dans certains cas exemple Dir l importance de la surface mouill e du fleuve est telle ou en tiage la vitesse d coulement est trop faible pour qu un jaugeage soit r alisable Dans ces conditions en dessous d une certaine cote le d bit est r ellement inconnu faut viter que la routine TRAD n extrapole la courbe de tarage vers le bas en dessous de la cote significative limite On introduit un test au d but de TRAD comparant la hauteur la cote inf rieure limite de calcul du d bit Toute hauteur inf rieure entra ne la convention du d bit n gatif gal 20 se diff renciant ainsi du 10 repr sentatif de l absence de relev s Ce d bit 20 est trait ensuite dans le programme comme le d bit 10 impression perforation etc Toutes les instructions correspondantes sont int gr es dans le 328 Vex clusion de l instruction test IF HAU LT 0 60 GOTO 129 qui ne doit tre plac e dans la sous routine TRAD qu en cas de besoin seulement entre D 130 J HAU HAUT et D 115 L 2 KMAX 5 4 3 TRAITEMENT PAR LA M THODE DU GRADIENT LIMNIGRAPHIQUE IUNIV 3 Alors que la m thode de au maximum tait simple mais applicable
96. eux courbes pour effectuer le d une station deux chelles L une dite de tarage d nivel e normale donne Qn en fonction de Dn L autre dite de correction de pente donne Q Q en fonction de D Dn Ces deux courbes remplacent avantageusement le faisceau de courbes sch matis sur la figure 1 dans lequel chaque courbe Q H est trac e pour une certaine d ni vel e Di l inconv nient l usage d un tel faisceau r side dans la n cessit d inter poler souvent pour estimer un d bit entre deux courbes voisines les d nivel es observ es ne correspondant que rarement aux valeurs Di choisies Il n y a pas d obligation quant au choix de l chelle de r f rence amont ou aval On retient cependant g n ralement l chelle amont plus loign e des variations du contr le aval Cette chelle sert de r f rence pour le calcul des jaugeages et des tarages La seconde chelle aval ne sert qu au calcul des d nivel es Di Hmi La proc dure suivre pour tracer les deux courbes est la suivante cf fig 5 et 6 a Choisir la d nivel e normale Dn b Reporter tous les jaugeages Qi Hi Di sur un m me graphique Hi tant la hauteur l chelle de r f rence amont Calculer pour chacun d eux la valeur que l on aurait mesur e si la d ni vel e avait t pour la m me cote Hi soit Qi 22 d Reporter les points calcul s Qn Hi ainsi calcul s
97. f rieures des tron ons de paraboles HP L Ces nouvelles valeurs de IUNIV remplacent celles qui avaient t mentionn es dans la note technique n 18 51 38 La pr sence de IUNIV non nul provoque galement l impression de station non univoque une chelle en sortie juste apr s l identification de la station suivie de l indication de la m thode de tarage utilis e Ainsi modifi le 302 est utilisable pour toutes les stations une chelle qu elles soient ou non univoques 5 2 2 LE PROGRAMME 311 r alise le calcul sur ordinateur des quations des trois courbes de tarage pour une station non univoque dot e de deux chelles dont l talonnage a t d crit dans le chapitre 2 Les trois courbes de ce tarage sont les suivantes a Une courbe enveloppe hors influence de pente pour laquelle on fait IUNIV 4 b Une courbe d nivel e normale ayant IUNIV 5 Une courbe de correction de d nivel e ayant UNIV 6 Les deux premi res courbes sont des courbes de d bits en fonction de hauteurs la troisi me est une courbe de rapports de d bits en fonction de rapports de d nivel es D DN Les trois courbes sont trait es comme des successions de tron ons de paraboles Le traitement des deux premi res courbes pour leur mise en quation est totale ment identique celui du programme 302 Le traitement pour mise en quation de la courbe de c
98. hode de l cart la cote maximale AH 25 26 La figure 16 montre simplement la similitude des m thodes pour des cours d eau limnigrammes homoth tiques et justifie en quelque sorte que l on remplace pour une cote H le jour j le gradient limnigraphique Ahn H RER Hi_n par l cart au maximum Hx Hx 4 2 2 MODE D EMPLOI Le syst me d quations 6 et 7 est remplac ici par le syst me suivant Q Qo 1 f Ahx 8 9 f Ahx tant une fonction de Ah qu il importe de d terminer et dont la forme n est pas connue a priori En fait si l on assimile le gradient un incr ment de la fonction hauteur l cart la cote maximale est assimilable son int grale et l on peut supposer que la liaison 9 sera de la forme AQUA gt 1 Comme l cart Ahx est toujours positif on doit lui affecter un signe positif lorsque l on est en crue et n gatif lors de la d crue Bien entendu l affectation de signe inverse est possible si la station a un type d coulement pente sup rieure en d crue par rapport la crue cas rare de confluence complexe avec crues des affluents d cal es comme on l observe sur le Tinkisso avant qu il se jette dans le Niger la station de Ouaran La fonction simple qui permet le mieux de repr senter les formes de tendance A parabolique sym triques par rapport au z ro des relations Ce f Ah
99. ilit d installer oui non oui non non deux chelles M thode deux Quant A deux k Ecart de tarage chelles graphique chelles graphique graphique au Maximum L tude d crite dans cet article inclut seulement les diff rents types de stations hydrom triques non univoques pr sent s dans ce chapitre Un seul type a t cart il s agit des stations soumises des variations cycliques de la pente c est dire des stations implant es dans des biefs o l influence de la mar e est importante Ce type de station exige deux chelles mais son tarage et son exploitation sont plus compliqu s que ceux des stations des groupes a et b bien qu ils en soient d riv s 2 Tarage des stations dot es de deux chelles sur cours d eau sujets d importantes et brutales variations de pente 2 1 Point de vue th orique En outre des stations du groupe faible pente peuvent tre soumises l in fluence de la mar e Dans tous les cas l inversion de la pente d coulement s observe r guli rement C est un ph nom ne rare hors de l influence de la mar e mais observable affluents de grands fleuves tels que les cours d eau alimentant certains lacs du delta central lacustre du Niger Nous avons pas retenu dans cet article l analyse des tarages de ces cours d eau inversion de pente pour lesquels il faut prendre en compte la variation des volumes d eau accumul s et transitant
100. ion apr s changement crue Calcul du d bit moyen journalier selon valeur de CAS D tail sur fig 30 oui Passage Appel TRAD pour traduction l ann e suivante Appel TRAD pour traduction avant changement de tarage Appel TARAG nouveau tarage Appel TRAD pour traduction apres changement tarage changement tarage et Fig 29 Organigramme g n ral 328 Appel TRAD pour traduction apr s avant changement crue d crue Changement crue d crue KK KK 1 Appel TRAD pour traduction apres changement crue d crue GTA 148 Qu ICHA GTS calcul e CAS Traitement du CAS 1 Traitement CAS 2 8 Calcul JMA selon Calcul JMAX selon valeu valeur CAS MIS IBIS GTS Be 249 AN 5 20 E gt 0 lt 0 Appel GRAD 1190 D 159 1 ss Calcul d bit moyen journalier DS 170 LA 1 JX Calcul vecteur HG IT pas de temps 24 CAS Appel C RDEB 175 LA t JX Calcul vecteur HG pas de temps 25 CAS et 24 CAS Appel D 5159 LA 1 JMA avec pas de temps 24 CHNTINUE 70 C NTINUE Appel C RDEB Calcul d bit moyen journalie CONTINUE 5159 Traitement du CAS 9 Calcul HG
101. ion pour les deux m thodes ensemble sont donn es sur les figures 17 19 pour Mopti et 20 22 pour Dir L ajustement aux moindres carts par approximations successives des trac s des couples de courbes pour les tarages selon l une ou l autre des m thodes pourrait th oriquement donner des trac s l g rement diff rents pour les courbes de tarage proprement dites Qo H En fait les deux m thodes de correction sont tellement proches qu il n en est rien et que les courbes de tarage H sont confondues Par contre on peut voir que la diff rence avec la m thode du gradient limni graphique est notable et que l quation de correction tant bien distincte les courbes de tarage le sont galement Le tarage complet de la station non univoque par l une ou l autre des m thodes gradient ou cart au maximum ne doit cependant pas donner de diff rences syst matiques titre de v rification on a analys les volutions du d bit autour du maximum selon qu il est calcul par l une ou l autre des m thodes sur l exemple des crues de 1967 68 et 1968 69 Mopti et Dir fig 23 24 Compte tenu de la marge d incertitude sur les d bits on peut admettre que tous les r sultats sont concor dants Tout au plus peut on remarquer que la m thode de l cart au maximum respecte plus la r gularit de l volution des d bits que la m thode du gradient qui introduit par exc s de sensibilit quelqu
102. le minimum d carts par rapport aux points mesur s lors des jaugeages On en effet Qo 2 pour chaque jaugeage et par cons quent l application y de P quation 4 aux jaugeages doit permettre par approximations successives d obtenir les valeurs des d bits en r gime uniforme donc univoque Le calcul d un d bit quelconque lors de la traduction des hauteurs se fait ensuite par application de l quation 4 mise sous la forme 1 5 Cette forme de l quation 5 est identique celle de l quation 2 pr sent e au chapitre pr c dent cf 2 1 En effet dans celle ci la d nivel e D est 1 somme de la d nivel e normale D et d un gradient de d nivel e AD Dans ces conditions l quation 2 devient ADN PTa une longueur pr s la d nivel e repr sente la pente puisque D 1 1 ou D I U At pour reprendre la formulation pr c dente On voit donc que Dn est un coefficient pr s assimilable To U At et que les quations 2 et 5 sont de m me forme Ah Le gradient limnigraphique ES est donc parfaitement capable de remplacer un coefficient pr s la d nivel e D dans une quation de correction du d bit 3 2 Avantages et inconv nients La r solution par it rations de l quation 4 ou 5 est donc une op ration iden de la m thode it rative tique celle du tarage par 2 courbes d une station deux chelles Mais alors que c
103. n e suivante qui n aurait pas de vecteur HG ant rieur d fini comme on peut le faire avec le 301 et m me avec la m thode de l cart au maximum pour lesquels le traitement de chaque ann e est ind pendant La traduction des hauteurs en d bits s effectue ensuite dans le programme prin cipal et comme dans POH 301 avec appel TRAD Mais contrairement la m thode de au maximum TRAD n effectue ni la correction du d bit calcul ni le tri du d bit maximal op rations qui sont faites dans la sous routine sp ciale CORDEB Cette sous routine est appel e au moment du calcul du d bit moyen journalier qui est diff rent selon que l on op re en CAS 1 2 8 ou 9 La figure 30 d taille l or ganigramme du programme pour le calcul du d bit moyen journalier selon la valeur de CAS Cette distinction selon CAS est m me plus complexe que dans le 301 ou la m thode de l cart au maximum car il faut s parer les CAS 5 et 7 on a donc les options suivantes CAS 1 CAS 2 8 subdivis e en CAS 2 3 4 6 CAS 5 CAS 7 Le CAS 9 n est pas trait La structure interne de traitement de chaque option est d ailleurs identique aux d tails pr s elle comprend deux phases successives ex cut es dans une m me boucle d instruction D 5155 pour le CAS 1 161 pour les CAS 2 3 4 et 6 162 pour le CAS 5 et 167 pour le CAS 7 a L introduction dans le vecteur
104. nivel e normale en millim tres pour laquelle est calcul e 1 courbe d talonnage dite lt d nivel e normale format 14 43 44 5 3 2 ORGANIGRAMME DU POH 312 VERSION BRE est dessin sur la figure 28 Le listing complet de ce programme est donn en annexe L organigramme est assez sommaire puisque la logique interne pour plus de 90 des instructions est identique celle du 310 Cette logique t d crite en d tail dans 3 pages 106 107 figures 5 1 5 2 Nous allons simplement commenter les particularit s du POH 312 dues la version BRE d une part et la non univocit d autre part Les cartes de donn es sont lues dans l ordre suivant 1 carte d identification de l ensemble de BRE carte COH 501 1 carte d identification de la station dans l ensemble carte 502 13 cartes d talonnage appel sous routine cartes de relev s pour une ann e la carte des relev s HAUTB quand elle existe pr c de la carte des relev s HAUTA c est dire que l on trouve par exemple une s quence de num ros de cartes telle que 5001 1 5002 2 3 etc si seulement les hauteurs HAUT contenues dans les cartes 1 et 2 ont besoin d une correction de non univocit et par celles de la carte 3 etc 1 carte blanche en d ann e 1 carte blanche en fin de station 1 carte blanche en fin de travail La carte d identific
105. nnu 302 de calcul d talonnage b POH 301 de traduction en d bits des hauteurs observ es et d pouill es en pas de temps fixe carte COH 301 c POH 310 de traduction en d bits des hauteurs observ es et d pouill es en pas de temps variable dites relev s limnigraphiques int graux carte 311 Ces programmes ont t d crits dans les Cahiers d Hydrologie D 3 faut y ajouter pour tre complet le programme POH 327 permettant le traite ment mixte d une station normalement d pouill e en pas de temps fixe mais dont quelques crues ont fait l objet d un second d pouillement en pas de temps variable Ce programme t d crit dans le n 1 du volume de 1970 des Cahiers 7 SIR COULON 4 En r alit ce catalogue n est complet que pour les stations de cours d eau exploit zs en r seau ou pour des tudes particuli res et dont la mise en cartes per for es des donn es s effectue sur les mod les de la s rie 300 Il manque ainsi les stations de bassins repr sentatifs et exp rimentaux BRE pour lesquelles si la structure interne des cartes supports de donn es est la m me le codage d identification est diff rent ce qui n cessite des mod les de cartes un peu diff rents ceux de la s rie 500 5 Pratiquement les modifications apporter aux programmes de traitement de la s rie PGH 300 pour qu ils puissent traiter les stations perfor es sur cartes
106. nte qui peuvent tre faibles mais les variations relatives sont toujours tr s importantes On trouve dans ce groupe les cours d eau influence de v g tation les grands fleuves Niger sup rieur et S n gal moyen et inf rieur les petits cours d eau faible pente Mauritanie 1 2 2 CLASSIFICATION D APR S LA M THODE DE TARAGE Cette classification s appuye sur la pr c dente et sur la possibilit d installer ou non deux chelles dans le bief de la station On voit imm diatement que seules les stations pente notable des groupes a et b ci dessus peuvent tre dot es de deux chelles Si cela est possible la m thode de tarage est celle dite deux chelles Si cela n est pas possible la m thode est celle du gradient limnigraphique Les stations du groupe c pr c dant ne peuvent tre dot es de deux chelles Elles peuvent tre tar es par la m thode du gradient limnigraphique Pour celles de ces stations qui ressortent du type grand cours d eau r gulier faible pente une m thode de tarage simplifi e dite de l cart 1 cote maximale peut tre appliqu e Le tableau suivant r capitule ces deux classifications accroissement diminution de la pente de ja pente grande variabilit relative de 1 pente Type de non univocit Type de 1 faible pente cours d eau a pentenotable a pente notable irr gulier r gulier Possib
107. ompte tenu des conditions hydrologiques d installation de deux chelles il est vraisemblable que la rivi re aurait un r gime justiciable du pas de temps variable dans la grande majorit des cas C est pourquoi nous avons seulement crit le programme correspondant ce type de d pouillement 41 60156501 GUYANE ORAPU CRIQUE VIRGILE CRIQUE VI GILE STATION PENTE VARIABLE ETAL NG 1 VALEURS CE COURBE DE CORRECTION DE DENIVELEE 42 C 1 1 CI2 L 013 1 L 1 0 8981483E 02 0 1306111 02 0 0 5 2 0 2243753 01 0 1512500E 01 0 4480000 L 3 0 5733351 00 0 7953336 00 0 6300000 4 0 2666681 00 0 5266671 00 0 7940000 5 0 1599S87E 00 0 40395 lt 7 00 0 9090000 6 0 8000225 01 0 331997GE 00 0 1000000 0 1200426 00 0 2980122 00 0 1077999 8 0 2665913E 01 0 2506642E 00 0 1145000 9 0 4000957 01 0 2260020 00 0 1205999 VALEUR DENP L i L RAPPORT D ON 1 0 0 6 1 000 2 0 090 7 1 250 3 0 250 1 500 4 0 500 9 1 750 5 0 750 10 2 000 BAREME RAPPORTS DENIVELEES_RAPPORTS DEBITS RAPPORT D DN RAPPORT 0 0 0 0 0 0 0 100 0 463 6 200 0 586 0 300 0 669 0 400 0 737 0 500 0 794 0 600 0 844 0 700 0 889 0 800 0 929 0 900 0 966 1 000 1 000 1 100 1 03 1 200 1 06 1 300 1 09 1 400 1 12 1 500 1 14 1 600 1 17 1 700 1 19 1 800 1 22 1 900 1 24 2 000 1 26 2 100 1 28 2 200 1 30 2 300 1 32 2 400 1 34 2 500 1 35 2 600 1 37 2 700 1 38 2 800 1 40 2 900 1 41
108. ons de BRE requi rent toujours le pas de temps variable et les stations ordinaires la plupart du temps le pas de temps fixe mais parfois l un des deux autres variable ou mixte Ainsi des programmes du catalogue pour stations univoques le POH 310 a vraiment besoin d tre modifi pour le traitement des BRE Ces r flexions sur la nature des stations et sur le pas de temps de d pouillement tant faites on peut examiner le catalogue des programmes de traitement automa tique des stations non univoques Pour cela on a repris la classification de ces stations selon la nature de la non univocit et la m thode de tarage applicable qui a t dress e dans le chapitre 1 par 1 2 Le tableau X pr sente un catalogue exhaustif des programmes n cessaires pour le traitement automatique des stations non univoques Les subdivisions introduites tiennent successivement compte du type de non univocit du type de cours d eau de la possibilit d avoir deux chelles du mode de pas de temps de d pouillement et de la m thode de tarage Bien qu il y ait dix subdivisions le nombre des programmes est seulement de six pour le calcul des d bits a Trois pour les stations deux chelles selon les pas de temps 5 Trois pour toutes les autres stations une chelle trait es par la m thode du gradient limnigraphique ou par celle de l cart au maximum et cela quel que soit le type de non univocit Ce nombre est de de
109. ormer l quation 6 en AQ _ Ah GA 0 Cette quation lin aire exprime que l cart relatif du d bit varie proportionnelle ment au gradient limnigraphique Il est alors possible d op rer graphiquement comme il est souhaitable pour effectuer le tarage et ceci en deux temps On trace une courbe peu pr s m diane entre les points de jaugeages Qie mesur s en crue et ceux mesur s d crue Hi courbe de tarage qui est prise comme base de pseudo univocit et donne Qo en fonction de b L cart AQ entre le d bit jaug ou Qia et Qo est rapport et le rap Ah port port dans le graphique de correction de pente en face de 0 On proc de par t tonnements dans le but d obtenir la meilleure r duction des carts sur ce second graphique au choix du pas de temps t Celui ci est d autant plus long que le cours d eau est lent et a un fort d bit 17 18 L ajustement graphique des trac s des deux courbes H et 25 se fait par la m thode des d viations r siduelles en conservant la premi re courbe son allure parabolique et la seconde celle d une droite 4 1 2 APPLICATIONS SUR TROIS EXEMPLES On a appliqu la m thode du gradient limnigraphique d une part aux stations de Mopti et de Dir sur le Niger Ces deux stations ont t choisies parce qu elles illustrent parfaitement le cas particulier des cours d
110. orrection de d nivel e a une logique interne identique celle du programme 302 mais les variables tant exprim es en milli mes qu il s agisse des rapports de d bits variables QQNP et QQNINT ou des rapports de d nivel es variables DDNP et DDNINT les formats d entr e et de sorties sont diff rents ainsi que quelques instructions arith m tiques dans lesquelles entrent en ligne de compte les unit s d expression des variables L articulation des l ments du programme POH 311 s effectue comme le montre le sch ma simplifi ci dessous fig 25 un test sur UNIV envoie vers une logique type 302 identique ou un peu modifi e selon la courbe mettre en quation il y a une partie de programme commune pour l impression des coefficients de paraboles en format sp cial E galement emprunt e au POH 302 Le programme 311 permet le calcul des 3 courbes de tarage d coup es au plus en 15 tron ons chacune pour autant d talonnages ou de stations qu on le d sire Les donn es d entr e doivent tre pr sent es dans l ordre suivant Carte d identification de station 302 Carte sp ciale d talonnage NETAL IUNIV 513 Carte L hauteurs limites 12 16 4 2 ou carte DDNP L rapports limites de d nivel es 12 16 F 4 3 39 40 Lecture carte d identification etcarte COH 513 Contr le de fin de station non Mise
111. par le contr le de fin d ann e 320 NETA 111 1000 111 qui envoie sur l instruction 1000 on lit nouveau une Ire carte 511 Un second test de carte blanche dirige S il est positif sur traitement d une nouvelle station instruction 99 S il est n gatif sur le test IF NGCA 5000 pour reprendre le traitement d une seconde ann e instruction 40 Les sorties du 312 sont identiques celles du 310 qu il s agisse des d bits imprim s ou perfor s dat s journaliers mensuels annuels et maximaux 5 4 Programme de calcul des d bits d une station non univoque une seule chelle relev es partir de hauteurs pas de temps fixe 328 5 4 1 DOMAINE D APPLICATION Ce programme est une variante du 301 programme de traitement syst ma tique des hauteurs relev es pas de temps fixe 1 8 fois par jour CAS 1 8 9 16 fois par jour CAS 9 et contenues dans les cartes mod le 301 Ce programme calcule les d bits correspondant chaque relev puis les d bits moyens journaliers qu il imprime et perfore calcule galement apr s interpolation des d bits manquants les d bits moyens mensuels et annuels qu il imprime et perfore 2 Le programme 328 est tr s complet puisqu il permet de traiter une station non univoque une chelle quelle que soit la m thode de correction de non univocit par l cart au maximum
112. part un r tablissement de la position du maximum d autre part consiste pr lever sur le limnigramme un nombre de points bien sup rieur au nombre n cessaire Cette multiplication des points att nuera les dents de scie et facilitera le lissage Bien entendu cette multiplication des pr l vements est inopportune car consom matrice de temps si l on proc de manuellement au d pouillement des limnigrammes et la perforation des cartes de donn es Il en est tout autrement si l on op re partir d un lecteur de courbes semi automatique ou automatique L introduction en d but de 1971 d un tel lecteur de courbes dans l quipement du Bureau Central Bibliographie en faisant abstraction de At et en faisant k hydrologique de l O R S T O M nous permettra alors de mettre au point une pro c dure de lissage automatique pour le r tablissement de l hydrogramme r el Le grand int r t de la m thode du gradient limnigraphique est dans sa souplesse de correction du d bit En effet celle ci tient en une seule instruction AC DH Mais si l on convenait d une influence non lin aire de gradient sur le d bit toute autre instruction arithm tique pourrait tre substitu e celle ci sans que le reste du programme soit affect la lecture pr s du coefficient AC dans TARAG bien videmment Ainsi pourrait on tout particuli rement penser la m thode du gradient limni graphique non simplifi
113. pe serait inf rieure et il y aurait deux courbes pour les d nivel es sup rieures un certain seuil estimer Pour r soudre le probl me du tarage des stations non univoques munies d une seule chelle il faut trouver un param tre ais ment mesurable d apr s les observations faites cette chelle et qui soit li la pente d coulement Le premier param tre qui vient l esprit est la vitesse de variation de la hauteur Ah d eau lue l chelle que l on peut exprimer par Ah tant la variation de la hau At teur h pendant l intervalle de temps At Cette vitesse est repr sent e par la pente du limnigramme c est pourquoi nous lui avons donn le nom de gradient limnigraphique afin de r server le mot pente la ligne d eau d coulement On peut th oriquement justifier le r le majeur du gradient limnigraphique Cette approche th orique a t trouv e par Y BRUNET MORET qui nous commu niqu e pour l introduire dans cette tude On donne les symboles Q 1 aux valeurs observ es du d bit de la vitesse moyenne et de la pente d coulement et Uo Lo les m mes valeurs en r gime uni forme En supposant le r gime uniforme on peut crire d apr s l quation de CHEZY K I pour une certaine hauteur h l chelle Si le r gime cesse d tre uniforme c est dire si le plan d eau monte de Ah dans l intervalle de temps At on observe une variation A
114. r tre sans effet sur l coulement quand la cote P chelle aval est inf rieure la diff rence d altitude entre les z ros des deux chelles En adoptant les signes suivants cote absolue l chelle amont cote absolue l chelle aval Zm cote absolue du z ro de l chelle amont Zy cote absolue du z ro de l chelle aval hm cote observ e l chelle amont hy cote observ e l chelle aval On peut crire hm Hm Zm hy H Z et la d nivel e D devient si l on pose dZ Zm Ze D hm dZ hw La pente de l coulement repr sent e par D influe l coulement tant que gt dZ et cesse de l influencer quand hv lt dZ On a alors D gt hm Il existe dans ces condi tions ne courbe de tarage vers laquelle tendent asymptotiquement les courbes Qi Hi d nivel es croissantes courbe de tarage correspondant un tat univoque du cours d eau et que nous appelons courbe enveloppe On remarque fig 7 que les courbes Qi Hi Di commencent sur cette courbe enveloppe lorsque la cote d passe la d nivel e Di Cette situation univoque asymptotique peut pratiquement ne pas exister sur un cours d eau contr le aval variable Seuls les jaugeages permettent de le v rifier et le cas ch ant de tracer cette courbe enveloppe En fait on peut penser que le type de cours d eau susceptible de recevoir deux chelles assez de pente
115. s de temps du gradient peut varier de 1 2 journ e 5 jours Un pas de temps horaire para t tout indiqu pour IT Dans cet exemple la hauteur limite inf rieure de calcul du d bit est 0 60 m 53 54 Il faut donc disposer d un vecteur hauteurs temps HG d une longueur d au moins 5 jours et dans lequel il y a une hauteur chaque heure ou un code d ab sence de hauteur Ce code est tant t 99 99 quand il y une absence de relev ainsi cod e sur la carte 301 9 99 quand il n y a pas de hauteur pour une certaine valeur de La programmation pr sente de ce fait un certain nombre de particularit s que nous allons maintenant d crire En d but de programme principal le vecteur HG est rempli dans la totalit de son champ 300 heures de 9 99 valeurs qui ne seront effac es que pour les moments auxquels correspondront des relev s de hauteurs Ensuite il faut initialiser r ellement le vecteur HG avant de traiter la hauteur du Ier janvier de l ann e consid r e On doit donc lire au pr alable au moins 5 jours de relev s ant rieurs Si ceux ci existent pas de probl me 8715 n existent pas cas de la premi re ann e de relev s d une station il faut les estimer de telle sorte que le gradient soit plausible Cette estimation doit au moins porter sur la premi re valeur afin que le gradient soit toujours d fini Par mesure de simplification les donn es ant rieures du vecteur
116. s du bief section pente on applique la formule aux sections des deux chelles et si l on appelle Q l d bit pour une cote Hm l chelle amont et Hy l chelle aval ces deux cotes tant exprim es en valeur absolue c est dire dans le m me syst me de nivelle ment le d bit pour la m me cote amont mais pour une cote H diff rente l chelle aval est gal 0 Les d bits sont entre eux dans le rapport des racines carr es des d nivel es entre les deux chelles D Hm Hy et DI Hm Hie En r alit comme l coulement n est pas uniforme l quation de CHEZY n est qu une approximation et l exposant a du rapport des d nivel es est diff rent de 1 2 1 Dans la pratique on choisit une d nivel e normale Dn aussi proche que possible de ce que l on pense tre la valeur la plus fr quente mode des diverses d nivel es observ es ou observables compte tenu du r gime d coulement et des mouvements du contr le aval Si le cours d eau s coulait en r gime permanent graduellement 9 2 2 Proc dure d tablissement 10 vari avec une pente gale DA entre les deux chelles le tarage serait univoque Appelons le d bit correspondant cette d nivel e normale Un d bit quelconque Q pour une d nivel e quelconque D s exprime alors l aide de l quation 2 e gt 2 Q faut et il suffit de d
117. s param tres et JGRAD du gradient limnigraphique 1 carte JUDEBI d unit de calcul des d bits n cartes de relev s ant rieurs suivies d une carte blanche Les cartes de relev s de l ann e traiter 1 carte blanche en fin d ann e 1 carte blanche en fin de station 1 carte blanche en fin de travail Le pas de temps de calcul du gradient limnigraphique JGRAD doit tre de l ordre de quelques heures pour un petit cours d eau r gime irr gulier On not 6 heures pour l exemple de Boitiek 250 On peut penser que ce pas de temps JGRAD ne descend pas en dessous de 20 30 minutes et qu il ne d passe pas 48 heures On a donc d cid que l unit du pas de temps serait le dixi me d heure la minute paraissant superflue Le vecteur HG re u une dimension maximale de 800 dixi mes d heures 3 jours 1 3 Les relev s ant rieurs couvrent une p riode gale 48 heures 480 dixi mes d heures avant le 1 jour 0 heure du 1 relev traiter Comme dans le 328 le vecteur HG est pr alablement tous calculs rempli de 9 99 cette cote conventionnelle tant galement utilis e en perforation des relev s limnigraphiques int graux pour d signer que le cours d eau est sec d bit nul La r gle de perforation des relev s conditionne la m thode pour leur introduction dans le vecteur HG En effet cette r gle veut que quel que soit le nombre de relev s
118. sion gal Qo AO Den ms Qenm3 s 4 3600 37 IC AC 3000 SE Tarage de hautes eaux 54 3500 0450 Ke 955 a 0457 Yeso 340 2 204 R of H 61 248 2000 33006 2 2000 Gs 9 Hauteur en cm 200 d di 443 O 1 1000 304 A Tarage de basses eaux 95 Hauteur cm 0 100 200 300 400 b M OS 29 oe 32 Jaugeage de crue O Point rectifi Jaugeage de d crue Fig 9 et 10 M thode du gradient limnigraphique station de Mopti 19 20 Hauteur en cm 2000 Hauteur en cm Jaugeage de crue Point rectifi J age de d crue Jaugeage de d 95 AH5 en Correction de l talonnage par A Hn n 5 jours Fig 11 et 12 M thode du gradient limnigraphique station de Dir en m s Tarage de basses eaux 0 5 Hauteur en cm Hauteur en cm o 250 300 Application un petit cours d eau Station de BOITIEK AQ o Qo o Jaugeage de crue Q Point rectifie 50 Jaugeage de d crue AH encm Correction de l talonnage el Fig 13 et 14 M thode du gradient limnigraphique station de Boitiek 21 Le tableau suivant donne les quations des droites de correction par le gradient limnigraphique Dir 29 0 973 Mopti 26 0 372
119. st tr s peu vraisemblable pour le type de cours d eau consid r pente naturelle non n gligeable Une telle inversion serait possible dans une zone de confluence si l un des affluents avait un marnage tr s sup rieur l autre et si la station hydrom trique situ e sur ce dernier tait trop pr s du confluent Une telle implantation de station est tr s mauvaise et mieux vaut l viter L autre cas lui probable d inversion de pente est celui du bief soumis mar e que nous avions cart a priori justement pour cette raison En effet dans l hypoth se d inversion de pente on doit prendre en compte les volumes accumul s dans le bief dans l op ration de correction des d bits mesur s cette prise en compte est assez d licate et compliquerait suffisamment le probl me pour que nous estimions souhaitable de disjoindre le cas des cours d eau inversion de pente bief mar e et certaines confluences de celui trait ici Si l inversion de pente se produit pisodiquement durant une courte p riode et sans forte d nivel e n gative on peut admettre que la proc dure tablie est applicable en supposant Di 0 et par cons quent Qi 0 durant la p riode incrimin e Cela revient n gliger le volume accumul dans le bief et par cons quent l g rement surestimer le d bit mais il s agit d une erreur de second ordre qui reste admissible Nous avons d ailleurs rencontr cette inversion en traitant l
120. sur le m me graphique et tracer la courbe Qn Hi de d nivel e normale Dn Calculer Qi Qn et Di Da pour chaque jaugeage et porter ces points sur un second graphique afin de tracer la courbe de correction de pente Quelques ajustements sur le trac des deux courbes permettent de minimiser les carts L exposant a t estim partir de la valeur th orique 1 2 et tr s grossi rement apr s un ou deux essais sert simplement pr ciser la position du trac de la courbe Hn En r alit comme une courbe de tarage ne peut tre repr sent e par une quation unique l exposant a doit varier avec la hauteur l chelle On rectifie les trac s de basses et hautes eaux de la courbe Qn Hn en s appuyant sur les r sultats bruts de jaugeages 2 3 Limites d emploi La courbe de correction de pente n est pas consid r e comme une quation Q Qn D D mais comme une succession de tron ons de parabole de la m me fa on qu une courbe de tarage Les deux courbes de tarage qui viennent d tre tablies pourraient tre consid r es en th orie comme s appliquant l ensemble des champs de hauteurs d eau aux chelles et de d nivel es En r alit il est bon de prendre en compte deux limites pratiques l une relative l inversion de la pente d coulement l autre la cessation d effet de cette pente 2 3 1 L INVERSION DE LA PENTE L inversion de la pente d coulement e
121. tante de toutes les causes est principalement responsable de la non univocit On verra cf chap 3 que la vitesse peut jouer un certain r le La mesure de la pente doit donc permettre de r soudre le probl me Cette mesure de ia pente de la ligne d eau d coulement est possible dans le bief en amont du contr le aval d une station hydrom trique si trois conditions sont remplies a On trouve une pente suffisante c est dire de plusieurs centim tres au moins 5 10 sont des minimums souhaitables dans un bief de longueur raisonnable b Aucune modification du d bit ne se manifeste dans ce bief affluent c Le bief a une g om trie r guli re Ces trois conditions r unies il faut disposer des moyens mat riels quipement personnel accessibilit implantation permettant la mise en place de deux chelles ou de deux limnigraphes auxquels les observations doivent tre rigoureusement synchronis es Dans l impossibilit de mesurer la pente de la ligne d eau on peut utiliser la station une seule chelle en proc dant la d termination d un param tre li la pente ou suppos varier dans le m me sens qu elle Comme on le verra plusloin chap 4 ce param tre est tant t le gradient limnigraphique tant t l cart la cote maximale 1 2 1 CLASSIFICATION D APR S LA NATURE DE LA NON UNIVOCIT Une classification des types de stations non univoques peut tre maintenant bauch e en essa
122. ulement homog ne tarage univoque sur lequel se surimpose un autre r gime d coulement d des crues diff rentes d importance secondaire qui provoque une relative non univocit Le sch ma th orique de tarage comporte une courbe univoque de base et une surface de tarage d importance r duite se plaquant au dessus ou en dessous de l autre selon les valeurs des pentes surface plus ou moins explor e par chaque crue h t ro g ne fig 4 Les exemples de ce type abondent a Niger Malanville avec son r gime univoque d origine amont hydrogramme r gulier du Niger issu de la Cuvette Lacustre et le placage des crues locales affluents volta ques et dahom ens b Wabi Shebelli d Hamaro Hadad Imi Ethiopie avec son r gime homog ne du haut bassin basaltique recevant des crues violentes d oueds caract re torrentiel dans son cours moyen Si les crues du cours d eau ont seulement deux origines une principale et une secondaire et si chaque type de crue est homog ne la surface de tarage se r duit deux courbes univoques une pour chaque type de crue g n ralement confondues pour les basses eaux Les types de stations non univoques sont nombreux qu ils se diff rencient par leur cause premi re contr le variable coulement non permanent ou par leur cause seconde barrage seuil mobile zone de confluence r gime de crues etc Dans tous les cas la variation de la pente qui est une r sul
123. uler le gradient limnigraphique DH et ensuite la correction d bit le d bit corrig et de faire le tri du d bit maximal L op ration la plus d licate est le calcul du gradient limnigraphique La figure 33 en donne une repr sentation graphique sur un limnigramme fictif dessin HAUT LA I D o ecteuf Cote r elle il Temps IT 1G 1 Fig 33 Recherche du gradient limnigraphique Pour un relev HAUT LA le gradient limnigraphique DH est d nivel e entre cette hauteur et la cote qui a t observ e un certain temps auparavant ce temps JGRAD tant appel le pas de temps du gradient il est lu dans TARAG La recherche de la cote observ e l instant IG IT 1 JGRAD s effectue en explorant le vecteur partir de l instant final atteint par celui ci on fait 1 11 dans la boucle du programme principal Deux cas sont possibles a l instant IG il y a une cote r elle c est dire diff rente de 99 99 et de 9 99 et par cons quent le gradient se calcule instantan ment en faisant DH HAUT LA HG IG b I n y a pas de cote r clle l instant IG et par cons quent il faut chercher en remontant dans le temps puis en le descendant partir de IG la premi re cote r elle La cote r elle ant rieure permet le calcul du gradient ant rieur DH 1 et la cote r elle post rieure celui du gra
124. urs d eau sujets non univocit limit e de type ou ont une pente naturelle ou normale notable permettant souvent naturellement un r gime uni voque repr sent par la courbe pente normale laquelle correspond la valeur Valeurs de la pente Courbe pente gt normale volution possible dune crue gt H gt Fig 1 Surface de tarage d un cours d eau Fig 2 Tarage d un cours d eau r gulier contr le variable r gime non permanent Raccordement point de rebroussement a 7 Surface de tarage des crues h t rog nes Tarage univoque du r gime principal p He Fig 3 Tarage d un cours d eau r gulier r gime Fig 4 Tarage d un cours d eau crues d origines diff rentes non permanent Raccordement arrondi 1 3 Limites de l tude des stations non univoques la plus fr quente des pentes observ es La non univocit y est due une variation absolue importante de la pente qui en provoque une variation relative non n gligeable Non univocit g n ralis e couvrant les plages et C est le cas de tous les cours d eau faible pente n ayant pas la possibilit d un tat univoque mais seule ment d un tat m dian de la pente autour de laquelle se r pandent toutes les valeurs observ es de la pente La non univocit est due des variations absolues de pe
125. ux pour le calcul des talonnages dont l un est le 302 des stations univoques adapt aux stations non univoques une chelle et l autre le POH 311 pour les stations deux chelles Il y a donc un total de sept nouveaux programmes au catalogue Quatre sont pr sent s en d tail dans cet article ce sont les 311 d talonnage 312 328 et 329 de calcul des d bits Les trois autres programmes peuvent assez facilement s crire partir de ceux ci ou de ceux du catalogue des stations univoques on a fait figurer dans la colonne programmes du tableau X et entre parenth ses les programmes composer pour ces critures qui une exception pr s sont relatifs aux pas de temps mixtes Ces trois programmes seront crits ult rieurement quand le besoin s en fera sentir 37 l exception du 312 version BRE tous les programmes tablis et pr sent s ici concernent le traitement des seules stations ordinaires cartes de donn es sur COH s rie 300 i TABLEAU X CATALOGUE DES PROGRAMMES POUR STATIONS NON UNIVOQUES Type de Type de Stations paS Programmes non univocit cours d eau d pouillement de tarage M thode Programmes D importance Fixe POH 311 301 312 quelconque 2 chelles Variable 311 POH 312 Mixte POH 311 327 312 A pente notable TUNIV 4 5 6 A ou B et tr s variable Fixe POH 302 32
126. velle courbe puisqu elle contient et que le type de traitement d pend de la valeur de cette variable La figure 26 montre un fac simil de la carte COH 513 ETAT NTATIF BASSIN COURBE D ETALONNAGE PARAMETRES DE STATION A DOUBLE ECHELLE AA AB GRAD 7777 001010000 000000000000000000000000 0000 000000000000000000000000000000000000 12345 10 11 12 13 14 15 16117 18 19 20 21 22 23 24 29 26 27 28 7830 31 32 33 34 35 36 37 3829 40141 42 43 444546 47 4845 50 51 52 53 54 55 55 57 58 59 60 62 63 64 65 65 67 68 69 70 71 72 73 74 75 76 77 78 18 80 11111111 111111111111111111 11111111111111111111111111111111111111111111111111111111 HYDROGRAPHIQUE DIFFERENCE ENTRE ZEROS PENTE NORMALE ICHA 5 gt ETALONNAGE NO PARTE DE BASSIN BASSIN 8 REPRESE REPERE DE COURBE EXPERIMENTAL 22222222 222222222222 22222222272222222222 22222222222222222222222222222222222222 3333333333333333 3333333333333333335 445444444 4444444 4444444444444444444445 55555555 55555555 55555555555555555555555555 66666566565665665 S 66666886666666 1111111101111111111101111111111111111111 111111111111111111111111111111111111101 E 5 99999999999999999999999999999999999999999999999999999999999999999999999999999 99 234 o 18 301212223 2A atasa 4 4 3 B7 68 s9 70 12 7314256 7018 80 2 3 49505 4 1 IRM 67 046 DRS T Hudrolocia
127. x pour AH positif et n gatif est la fonction arc tangente On proc de comme dans la m thode du gradient limnigraphique en deux temps a Trac de la courbe de tarage au milieu des points de jaugeages de crue et de d crue donnant Qo H b Report des carts relatifs S Gs soit AQ en fonction de et ajustement graphique d une fonction a Arc tg b par la m thode des d viations r siduelles La fonction pr c dente permet une simplification laquelle ne semble pas a priori introduire d erreur grossi re en consid rant que l cart relatif de Qo varie lin airement en fonction de jusqu un maximum au del duquel AQ Qo reste constant Cela revient assimiler le limnigramme une courbe en toile de tente Le choix de cette correction par trois demi droites en lieu et place de la relation en Arc tangente peut tre introduit sous le nom de m thode simplifi e de l cart au maximum AQ 25 On a alors tant l cart maximal 0 4 2 3 APPLICATIONS AU NIGER DIR ET La m thode des d viations r siduelles est appliqu e avec les m mes symboles que pour le gradient limnigraphique cf 4 1 2 Les tableaux r capitulatifs des calculs d ajustement pour Mopti et Dir selon la m thode normale d cart la cote maximale ou selon la m thode simplifi e sont les tableaux VI IX Les courbes de tarage et de correct
128. yant de simplifier et d harmoniser les divers sch mas de tarage correspondant aux diverses causes de non univocit pr sent es au paragraphe 1 1 fig 1 4 On peut partir du raisonnement suivant une station non univoque est une station univoque subissant des contraintes venant d truire en totalit ou en partie son carac t re univoque Par cons quent toute station sans ces contraintes aurait un caract re univoque et une courbe de tarage univoque correspondant la pente d coulement la plus fr quente ou 1 m diane pente que nous qualifions de pente normale Tous les types de stations non univoques appartiennent l un des trois groupes suivants cf fig 1 a Non univocit limit e de type c est dire que les causes provoquent un accroissement de la pente et l apparition d une surface de tarage au dessus de la courbe pseudo univoque pente normale C est le cas des cours d eau modification g om trique naturelle du contr le aval par rosion c est le cas des cours d eau r gime de crues h t rog nes type Wabi Shebelli b Non univocit limit e de type c est dire que les causes provoquent une diminution de la pente et que la surface de tarage est limit e en dessous de la courbe pente normale C est le cas des stations de confluence g n ralement crique Virgile en amont c est le cas de certains r gimes crues h t rog nes Niger moyen Tous les co

Download Pdf Manuals

Related Search

Related Contents

スペース オルガナイザーモニターアーTKLA  Asrock FM2A85X Extreme4-M  spa lotus de luxe manual de uso e instalación  A User-Guide to Archi An Explicit Runge  Stretch DVI    Jabra Motion Office MS  Manuel d`installation Miele  Newstar USB 1.1 data switch, 4-port  symptome - Narr Isoliersysteme  

Copyright © All rights reserved. Failed to retrieve file