Optimisation sous contraintes d'un modèle hydrologique pour une représentation de la physique des processus

Les modèles hydrologiques sont utilisés sur une base quotidienne par nombre d'experts notamment pour dimensionner des ouvrages, planifier la production d'hydroélectricité ou limiter les dommages liés aux inondations. Ces simulations informatiques sont conçues afin de reproduire le cycle de l'eau par le biais d'équations afin d'arriver à une estimation fiable des débits à l'exutoire d'un bassin versant. Des variables d'optimisation sont présentes au sein de ces équations et la principale difficulté liée à l'utilisation de ces modèles consiste à trouver l'ensemble de valeurs produisant la série de débits la plus près de celle réellement observée. Cette opération se nomme calage ou étalonnage de modèles hydrologiques et offre une paramétrisation réalisée à partir de données historiques. Dans un contexte de réchauffement planétaire, des entreprises comme Hydro-Québec se doivent de produire des projections futures des débits afin de mesurer l'impact des changements climatiques sur leurs activités. Une paramétrisation basée sur le passé soulève un doute quant à son applicabilité en période future, où les événements hydrométéorologiques pourraient être très différents. De plus, la formulation actuelle du problème d'optimisation ne permet pas d'encadrer les variables intermédiaires associées aux processus hydrologiques, comme la neige, l'évapotranspiration ou le débit de base. Ainsi, un ensemble de paramètres est fourni à la simulation et seule la sortie du modèle, donc les débits, est considérée. L'impact est une représentation des processus hydrologiques ne respectant pas les lois de la physique, ce qui renforce l'hypothèse selon laquelle la présente paramétrisation du modèle ne peut, en toute confiance, être transposée en contexte de changements climatiques. L'objectif visé par ce mémoire est donc, par le biais des outils d'optimisation, de présenter une méthode de calage orientée vers la physique des processus et offrant une paramétrisation robuste en contexte de changements climatiques; notre hypothèse de recherche étant qu'un ensemble de paramètres offrant un respect de la physique dans la représentation des processus hydrologiques est préférable a une solution présentant une meilleure valeur de la fonction objectif, mais modélisant les processus en toute liberté. Ainsi, deux stratégies de calage ont été mises de l'avant et sont comparées à la méthode largement utilisée dans la communauté scientifique, notamment chez Hydro-Québec. La première stratégie proposée vient, par le biais d'une contrainte, recentrer le bilan de l'évapotranspiration réelle sur la valeur offerte par l'Atlas du Canada, alors que la seconde, toujours via une contrainte, met à profit l'immense base de données qu'offrent les modèles climatiques en contraignant la forme du cycle de l'évapotranspiration réelle. Cette approche par contraintes a été raffinée par l'utilisation d'un outil d'analyse de sensibilité des contraintes offert en complément de l'implémentation de MADS (NOMAD). Chaque contrainte fut ainsi traitée afin d'être relâchée ou resserrée selon la sensibilité de la fonction objectif et selon la confiance en la qualité des données observées, souvent entachées d'erreurs. Ce sont donc des milliers d'optimisations qui ont été réalisées par l'algorithme MADS afin d'arriver à une comparaison des différentes méthodes sur le plan de la juste représentation des hydrogrammes et des évapogrammes (ou cycle annuel moyen de l'évapotranspiration), ainsi que sur la variabilité de la paramétrisation et la robustesse des solutions. De ces analyses, il est ressorti que le calage parallèle offre la meilleure performance pour plusieurs raisons. Une forte corrélation a été remarquée entre la qualité des solutions produites par l'optimisation et celles produites en simulant les débits d'un ensemble de données n'ayant pas servi à la calibration, dites de validation. Ainsi, il est possible de conclure à une solide robustesse pour une utilisation en contexte de changements climatiques. La meilleure solution obtenue par la stratégie de calage parallèle s'est aussi révélée capable de produire une meilleure représentation des débits en validation que l'ensemble des stratégies, le tout en modélisant le processus d'évapotranspiration de façon à respecter la physique. L'application la plus notoire de la méthode développée dans ce mémoire est l'adoption de la stratégie de calage parallèle pour le projet (CQ)². Ce projet, collaboration d'Hydro-Québec, du consortium Ouranos sur la climatologie régionale et l'adaptation aux changements climatiques, du Centre hydrique du Québec (CEHQ) et de l'entreprise Rio Tinto Alcan, vise à qualifier et à quantifier les impacts des changements climatiques sur le régime hydrologique québécois. Pour y arriver, une paramétrisation unique a été optimisée à l'aide du calage parallèle pour chacun des 324 bassins versants recouvrant l'ensemble du territoire québécois. Au terme de ce projet, ce sont presque 160 000 scénarios hydroclimatiques qui auront été générés, soit 98 simulations climatiques employant 5 différentes méthodes de post-traitements déployées sur les 324 bassins versants. Cette base de données considérable n'a de valeur que si l'étalonnage du modèle est valide et c'est donc une grande marque de confiance des décideurs du projet envers les résultats obtenus dans cette recherche que d'avoir choisi notre méthode pour la calibration du modèle HSAMI. Les contributions de ce mémoire à l'avancement des connaissances en hydrologie sont nombreuses. La plus flagrante est la présentation d'une méthode de calage robuste pouvant être transposée en contexte de changements climatiques et pouvant être étendue à plusieurs autres processus hydrologiques pour une représentation encore plus orientée vers la physique. D'ailleurs, l'intégration de données complémentaires provenant de modèles climatiques dans le processus d'optimisation est proposée pour la première fois à la communauté scientifique. Le regard nouveau porté sur les processus d'optimisation a permis une amélioration non négligeable de la performance de l'algorithme d'optimisation par le biais de l'emploi approprié des options offertes. L'utilisation du module complémentaire d'analyse de sensibilité a aussi ouvert de nouvelles portes en offrant la possibilité de mettre à profit les connaissances des utilisateurs d'un modèle dans la formulation adéquate des contraintes. L'approche de calage par contraintes proposée dans cette recherche s'en est vue bonifiée et de nouvelles avenues de recherche en ont découlé. Finalement, pour le domaine de la recherche opérationnelle, ce projet a offert une belle opportunité d'application des outils de l'optimisation. Non seulement l'algorithme MADS est-il mieux compris par les chercheurs de l'IREQ, mais le module d'analyse de sensibilité s'est fait connaître et continuera d'être employé.