Variabilité événementielle des paramètres

Cette fonctionnalité permet de faire varier pour chaque événement un ou plusieurs paramètres de production ou de transfert, notamment en fonction des conditions initiales de l' hydrosystème pour chaque épisode : état de saturation des sols, évolution saisonnière de la végétation, travail du sol etc….

Un paramètre de production (ou de transfert) peut ainsi varier pour chaque événement, et cette variabilité est liée à l'utilisation d'une variable externe qui doit être renseignée pour chaque épisode : valeur estimée du paramètre, débit de base, indice d'humidité fourni par Météo-France, mesures d'humidité locales etc… Le paramètre en question peut être directement associé à la variable externe (dans le cas d'une valeur estimée du paramètre), ou être le résultat d'une opération effectuée sur la variable externe (dans le cas d'une relation entre le paramètre et le débit de base par exemple). Les opérations actuellement disponibles sont du type :

Y = aX+b

Y = aX b

Y = aLn(X)+b

où Y désigne la valeur du paramètre pour l'événement i, et X la variable externe associée à l'événement i

En pratique

Faire varier un paramètre d'un événement à l'autre impose de déclarer ce paramètre comme varX (où X peut prendre les valeurs 1 à 6). Cette syntaxe doit être obligatoirement respectée.

Les valeurs associées à varX apparaissent et peuvent être modifiées en activant l'icône  positionnée à la droite de la liste des fichiers de production (ou de transfert).

 

Cette icône ouvre une nouvelle fenêtre, où apparaissent pour chaque événement les valeurs associées aux variables var1 à var6 :

 

Les valeurs var1 à var6 sont initialement fixées à 0. Ces valeurs peuvent être modifiées :

• directement dans la fenêtre, en surlignant la ligne contenant la valeur à modifier. Cette ligne apparaît alors dans le bandeau supérieur, et peut être modifiée. Les modifications sont validées, après avoir activé le bouton "Appliquer"
• à partir d'un fichier externe, contenant les valeurs souhaitées. Ce fichier doit être préparé préalablement, puis appelé directement à l'aide du browser. La structure de ce fichier doit être la suivante :

2 1973 10150200 0 0 0 0 0 0 6 1976 09111500 0 0 0 0 0 0 7 1976 11082000 0 0 0 0 0 0 9 1977 10190500 0 0 0 0 0 0 12 1977 12071200 0 0 0 0 0 0 21 1994 09220100 0 0 0 0 0 0 22 1994 10190600 0 0 0 0 0 0 23 1994 11020900 0 0 0 0 0 0 etc….

où les colonnes 1 à 3 désignent le n°, l'année et la date début( mmjjhhmn) de chaque épisode, et les colonnes 4 à 9 les 6 valeurs de var1 à var6 associées à chaque épisode. Le séparateur des champs est un espace. 

NB : dans un premier temps, ce fichier peut être constitué en déclarant un nom quelconque dans le browser, et en validant l'opération par OK

La figure suivante montre un exemple où les colonnes var1 et var2 ont été complétées avec le débit de base en début de chaque épisode en colonne 1, et la valeur optimale du paramètre 1 en colonne 2.

Les simulations réalisées en déclarant var2 pour le paramètre 1 seront conduites avec des valeurs de ce paramètre différentes pour chaque épisode, et égales à la valeur déclarée dans la colonne var2.

Il est également possible de déterminer la valeur du paramètre, à partir d'opérations transformant les valeurs déclarées de varX, en activant le bouton  situé à la droite de chaque variable. Ce bouton a pour effet de proposer 3 opérations : Y = aX+b, Y = aX b, Y = aLn(X) +b

Dans le cas ci-dessus, l'opération utilisée est Y = aLn(X) +b avec a = -41.8 et b= 328.8, où Y désigne le paramètre qui sera effectivement utilisé dans le modèle, et X la valeur associée à var1. Ce cas correspond à l'utilisation de la régression déterminée entre le modèle S du SCS et le débit de base, dans le cas du Gardon à Anduze. Après avoir activé cette opération, le bouton placé à droite de var1 apparaît en rouge, pour signaler que les valeurs de la colonne var1 ne sont pas substituées directement au paramètre de production, mais transformées avant substitution. Cette option sera désactivée si on revient aux paramètres initiaux , c'est-à-dire Y= aX+b, a=1, b=1.

 

Les paramètres associés à des variables événementielles peuvent faire l'objet de procédures de calibration automatique, par optimisation d'une fonction critère. Les variables événementielles originales ou transformées seront utilisées comme valeur initiale du paramètre pour l'événement considéré.