Menu 4 : optimisation des paramètres

  • Généralités
  • Cliquez sur les différents champs de l'image pour avoir plus d'information

menu4

 

Le menu 4 propose un ensemble de méthodes numériques basées sur la minimisation d'un critère d'erreur, pour déterminer les paramètres "optimaux" du modèle. Les spécifications à fournir concernent :

 

Sélection des paramètres à optimiser: le nombre de paramètres à optimiser n'est pas limité. La sélection (ou désélection) d'un paramètre à optimiser est réalisée par double-clic gauche sur le champ correspondant à la valeur du paramètre. Dans la fenêtre ci-dessus, 3 paramètres ont été sélectionnés : param1 production de la classe granites, param3 production de la classe schistes, param1 transfert de la classe rivière.

 

Définition des bornes inférieures et supérieures: la recherche des paramètres optimaux n'est autorisée que dans l'intervalle des bornes inférieure et supérieure déclarées pour chaque paramètre. Ces bornes sont définies par paramètre, et sont donc identiques pour l'ensemble des différentes classes associées à ce paramètre.

 

Station hydrométrique de référence: la station hydrométrique de référence est utilisée pour le calage automatique des paramètres de production et de transfert. Les critères de calage sont calculés par rapport à cette référence. Choisir la station de référence à l'aide de l'icône : image003
On ne peut choisir qu'une seule station hydrométrique de référence. Les procédures de calage automatique des paramètres ne peuvent donc etre appliquées que pour un exutoire, correspondant à cette station de référence.

 

Type d'optimisation: les paramètres de production et de transfert peuvent être calés successivement (d'abord la production, ensuite le transfert) ou simultanément. Dans le premier cas, les paramètres de production sont déterminés par minimisation d'un critère d'erreur sur les volumes écoulés, puis les paramètres de transfert sont ensuite déterminés par minimisation d'un critère d'erreur sur les débits. Dans le second, les paramètres de production et de transfert sont calés ensemble, par minimisation d'un critère d'erreur sur les débits.

 

Evénements groupés ou non: les paramètres de production et de transfert peuvent être calés pour chacun des épisodes averse-crue sélectionnés, et dans ce cas, peuvent différer d'un événement à l'autre (non groupés). Ils peuvent également être déterminés pour l'ensemble des épisodes averse-crue sélectionnés, et sont dans ce cas identiques pour chacun des épisodes (groupés).

 

Domaine de calage: le calcul de la fonction critère peut être effectué sur différentes parties des hydrogrammes de crue, hydrogramme complet, débit maximum de crue, débits compris entre deux valeurs Q1 et Q2 et/ou entre deux instants t1 et t2.

 

Choix du critère: plusieurs critères d'erreur sont disponibles: EQM, EAM, Nash, CREC. En cas d'optimisation séparée des paramètres de production et de transfert, deux critères différents peuvent être utilisés. Le nombre d'itérations maximal doit être spécifié, et sera fixé en fonction du nombre de paramètres à optimiser, des premiers résultats obtenus, etc..

 

NB La procédure d'optimisation est désactivée dans les cas suivants :
- aucun paramètre n'est sélectionné pour l'optimisation,
- le nombre d'itérations est égal à 1.

Exemple: lorsque la procédure d'optimisation est activée, les calculs suivants apparaissent dans le fichier listing :

Optimisation simultanee des parametres de production et de transfert en cours....

Evenement no 1 du 15- 9-1994
Evenement no 2 du 18-10-1994
Evenement no 3 du 01-11-1994

No --Critere------- Seuil----parametre 1--parametre 2--parametre no 3

1 -0.298230E+00 0.131848E-01 0.100010E+03 0.253928E+00 0.100982E+01
2 -0.270198E+00 0.100364E-01 0.100014E+03 0.240179E+00 0.100786E+01
3 -0.270198E+00 0.101377E-01 0.100009E+03 0.237560E+00 0.101113E+01
4 -0.270198E+00 0.167181E-01 0.100004E+03 0.227957E+00 0.101135E+01
5 -0.235492E+00 0.257048E-01 0.100002E+03 0.207224E+00 0.101222E+01
6 -0.189040E+00 0.236732E-01 0.100005E+03 0.174706E+00 0.101899E+01
7 -0.189040E+00 0.408304E-01 0.100005E+03 0.146553E+00 0.102030E+01
8 -0.876965E-01 0.343303E-01 0.999998E+02 0.874091E-01 0.102794E+01
9 -0.876965E-01 0.282667E-01 0.999968E+02 0.396866E-01 0.103848E+01
etc.........

La variable seuil désigne la distance euclidienne entre deux jeux de paramètres testés successivement. La procédure d'optimisation s'interrompt quand ce seuil descend en dessous de 10-6.