Corps de l’article

1. Introduction

La température de l’eau intervient de façon indirecte mais déterminante dans la pollution des eaux de surface. Par exemple, l’eau possède une capacité d’auto-épuration, qui est contrôlée essentiellement par sa température. L’auto-épuration de l’eau est limitée par la quantité d’oxygène dissoute. L’eau étant très polaire, elle ne peut dissoudre que des composés chimiques ayant une polarité différente de zéro, en établissant des liaisons de forte énergie de type hydrogène. Ainsi, l’oxygène est peu soluble dans l’eau, à cause de sa faible polarité par rapport à celle de l’eau. En eau douce, la solubilité de l’oxygène dans l’eau est de 14,4 mg/L à 0 °C. Elle diminue rapidement avec une élévation de la température. À 30 °C, elle n’est plus que de 7,53 mg/L.

Lorsque la température augmente, la quantité d’oxygène diminue, ce qui provoque le développement d’autres micro-organismes capables de vivre en anaérobiose (absence d’oxygène) avec pour conséquence :

  1. la formation de gaz putrides malodorants;

  2. la disparition de certains organismes modifiant la diversité des espèces; et

  3. la réduction de la capacité de charge et d’auto-épuration de l’eau.

En somme, la température de l’eau est un paramètre très important pour toute étude qualitative d’un cours d’eau. Ainsi la connaissance, la prédiction et la prévision de la température de l’eau de rivière représentent un atout capital.

Les applications des réseaux de neurones dans le domaine hydrométéorologique, à des fins de prévision, concernent généralement la prévision des débits des rivières, comme l’estimation des précipitations, la prévision des apports naturels aux réservoirs d’irrigation. Dans le domaine de la prévision de la qualité de l’eau, les études réalisées portaient sur la prévision de la salinité (Maier et Dandy, 1996), l’acidité et la conductivité (Bastaracheet al., 1997) et sur le pouvoir de dissolution de l’eau des produits chimiques à base de carbone et d’azote (Clair et Etherman, 1998).

Dans ce travail, nous proposons un modèle de prédiction des variations de la température d’un cours d’eau en fonction de variables climatiques telles que la température de l’air ambiante, le débit d’eau et la quantité des précipitations se déversant dans le cours d’eau. Des réseaux de neurones ont été utilisés pour approximer la relation entre ces différentes variables. Nous avons considéré trois paramètres, disponibles dans nos données historiques relatives à la température de l’eau relevée sur le site étudié. Ces paramètres sont jugés pertinents : la température de l’air (Cluis, 1972; Mohseni et Stefan, 1999; Piperet al., 1982; Stefan et Preud’homme, 1993), le débit d’eau et la quantité de précipitations enregistrés au cours de la journée (Caissieet al., 2001).

Dans un premier temps, nous proposons un modèle purement statique, considérant que la température de l’eau ne dépend que des variables mesurées au cours de la journée. Ce modèle est simple car il est basé sur l’estimation de la relation entre les différentes variables, qui sont mesurées à chaque instant. D’autre part, les processus de variation de la température de l’eau et de l’air sont des systèmes évolutifs dans le temps. Il est alors judicieux d’en tenir compte en introduisant d’autres variables liées à la dynamique de ces régimes. Ainsi, nous proposons un autre type de modèle pour prédire la température de l’eau à différents horizons dans le futur jusqu’au jour (j + 7).

2. Description du site

Les données disponibles proviennent du site du ruisseau du Panola Mountain. La région est située à 25 km au sud-ouest d’Atlanta (E.-U.). Le climat de la région est considéré comme subtropical tempéré avec une température moyenne de 16 °C et une moyenne de précipitation de l’ordre de 1,24 m (Hopper, 2001). Les données historiques relevées sur le site (la température et le débit d’eau de la rivière, la température de l’air et la quantité de précipitations) couvrent la période de 1998 à 2003. Notons la présence de données manquantes sur l’ensemble de cette période. Les valeurs horaires des trois premières variables relevées pendant la journée sont converties en moyenne journalière et la quantité de précipitations est transformée en une variable cumulative de la journée. La journée commence à 1 heure du matin, et se termine à minuit. La figure 1 représente l’évolution de la température durant la période considérée. Quelques paramètres statistiques de ces données sont indiqués dans le tableau 1.

Figure 1

Variation de la température de l’eau durant la période de l’étude.

Variation of the water temperature during the period under study.

Variation de la température de l’eau durant la période de l’étude.

-> Voir la liste des figures

Tableau 1

Quelques paramètres statistiques descriptifs relatifs à la température de l’eau en fonction des années de l’étude.

Statistical parameters of the data set corresponding to the water temperature during the period studied.

Année

Variables

1998

1999

2000

2001

2002

Nb de jours

365

360

301

238

236

Tmoyenne

15,3

15

15,8

16,2

16,2

Tmin-max

(8,7‑20,7)

(5,5‑21,6)

(7,7‑21,1)

(7,0‑19,9)

(7‑20)

-> Voir la liste des tableaux

3. Méthodologie des réseaux de neurones

3.1 Réseaux de neurones : introduction

Le problème que se posent les scientifiques en recherchant des modèles déterministes de type « cause à effet » est le suivant : à partir d’un ensemble de mesures de variables d’un processus de nature quelconque (physique, chimique ou biologique) et des résultats de ce processus, comment trouver une relation déterministe entre ces variables et ce résultat, représentée par une équation mathématique, valable dans le domaine où les mesures ont été effectuées. Dans notre étude, nous utilisons une technique qui a fait ses preuves en calcul statistique d’approximation fonctionnelle : les réseaux de neurones. Elle permet de faire des prévisions, d’élaborer des modèles, de reconnaître des formes ou des signaux, etc. (Dreyfuset al., 2002).

Dans notre cas, les données disponibles forment une série chronologique dont on peut effectuer des prévisions des événements futurs, mais dont les facteurs influençant les événements passés doivent rester les mêmes. À notre connaissance, aucun travail n’a été fait dans ce sens. Par contre, la problématique de la prévision en hydrologie a été largement étudiée, surtout le cas de la relation débit-précipitation (Coulibalyet al., 2000; Luket al., 2000; Tothet al., 2000; Zealandet al., 1999).

On distingue deux grands types d’architecture des réseaux de neurones pour résoudre ce type de problématique : les réseaux de neurones non bouclés et les réseaux de neurones bouclés. Dans notre cas, nous utiliserons une variante de la deuxième catégorie en introduisant comme variables d’entrée le résultat du processus de l’événement passé dans une nouvelle fonction. Une seule architecture de type réseaux à couches (MLP : un réseau de neurones de type perceptron multicouche) sera utilisée (Figure 2) : une couche de neurones est connectée à une couche de neurones « cachés », qui est reliée à une couche de neurones de sortie. L’activité des neurones d’entrée code l’information qui est présentée au réseau. Ce sont les variables explicatives externes citées précédemment. L’activité de chaque neurone caché est déterminée par l’activité des neurones d’entrée et les valeurs des coefficients de connexion ou poids. L’activité des neurones de sortie dépend de l’activité des neurones cachés et des coefficients de connexion reliant ces neurones à ceux de sortie. La couche de sortie contient un seul neurone, représentant la température de l’eau.

Figure 2

Représentation d’une simple architecture à couches de type MLP des réseaux de neurones.

Representation of a simple MLP architecture of the neural network.

Représentation d’une simple architecture à couches de type MLP des réseaux de neurones.

-> Voir la liste des figures

Un neurone est avant tout un opérateur mathématique, dont la valeur numérique se calcule par quelques lignes d’un programme. Un neurone réalise une somme pondérée suivie d’une fonction non linéaire f. Cette fonction f doit être bornée, continue et dérivable. Elle peut aussi avoir la forme d’une fonction seuil si le résultat recherché est de type booléen (soit 0 ou 1). Les fonctions les plus fréquemment utilisées sont les fonctions sigmoïdes.

3.2 Principe d’un modèle de réseaux de neurones

Le principe est le même que pour n’importe quel processus statistique caractérisé par une relation déterministe entre des causes et des effets. Des exemples d’apprentissage, c’est-à‑dire des jeux de valeurs de neurones d’entrée et de valeurs des neurones de sortie correspondants, sont fournis au réseau de neurones. L’apprentissage consiste tout simplement à calculer les coefficients des connexions (poids) entre les différentes couches de telle manière que les sorties du réseau de neurones soient, pour les exemples utilisés, aussi proches que possible des sorties désirées. Pour ce faire, on utilise des algorithmes d’optimisation qui minimisent une fonction de coût constituée d’une mesure de l’écart entre la réponse réelle du réseau et la réponse désirée. Cette optimisation se fait de manière itérative, en modifiant les poids en fonction du gradient de la fonction du coût : le gradient est estimé par une méthode spécifique au réseau de neurones, dite méthode de rétropropagation de l’erreur. Les poids sont initialisés aléatoirement avant l’apprentissage, puis modifiés itérativement jusqu’à obtention d’un compromis satisfaisant entre la précision de l’approximation sur un ensemble de validation disjoint au précédent et celle obtenue sur l’ensemble des données d’apprentissage.

3.3 Mise en oeuvre des modèles de réseaux de neurones

Pour réaliser un modèle de réseaux de neurones, trois étapes successives sont nécessaires :

  • Tout d’abord, on choisit l’architecture du réseau, c’est‑à‑dire le nombre de neurones de la couche cachée. Le réseau ne doit être ni trop souple, ni trop rigide. Il existe maintenant plusieurs méthodes adaptées à ces considérations, telles que la régularisation bayésienne. Nous disposons de plus d’un millier de données pour un nombre de neurones d’entrée ne dépassant pas six variables. La fonction dans le modèle neuronal pourrait être définie sans surajustement, car le nombre de données est assez élevé et la température de l’eau ne varie pas de manière significative d’un jour à l’autre.

  • Les poids sont calculés en minimisant l’erreur d’approximation sur les données de l’ensemble d’apprentissage, de telle manière que le réseau réalise la tâche désirée.

  • Enfin, la qualité du réseau obtenu est estimée à partir d’exemples ne faisant pas partie de l’ensemble d’apprentissage.

Dans notre étude, les deux fonctions d’activation utilisées sont les suivantes : une fonction de type sigmoïde dans la couche intermédiaire et une fonction linéaire dans la couche de sortie. Des travaux précédents ont montré que ce couple de fonctions permet d’approximer quasiment tous les types de relations non linéaires (Dreyfuset al., 2002).

Pour enseigner une tâche à un réseau, il faut ajuster les coefficients de chaque neurone et minimiser la différence entre la sortie désirée et la sortie effective. Ce procédé impose de calculer la dérivée d’une quantité nommée coût J par rapport aux coefficients de connexion, c’est‑à‑dire estimer la variation de l’erreur en fonction des variations de chaque coefficient de connexion. La méthode la plus utilisée pour déterminer ces dérivées est la méthode de rétropropagation de l’erreur.

Une fois que les activités de tous les neurones de sortie ont été déterminées pour l’ensemble des données d’apprentissage, le calcul du coût J est défini par l’expression suivante :

où y est la température observée, g(x,w) est la température calculée, x et w représentent respectivement l’ensemble des variables explicatives et les coefficients de connexion ou poids, N est le nombre d’exemples et k est l’indice des données de l’ensemble d’apprentissage.

Évaluer le gradient consiste à évaluer le gradient du coût partiel Jk(w) relatif à l’observation k, et de faire ensuite la somme sur tous les exemples. L’algorithme de rétropropagation consiste essentiellement en l’application répétée de la règle des dérivées composées.

δkj désigne la valeur du gradient du coût partiel par rapport au potentiel du neurone i, et ce pour l’exemple k. Il peut être exprimé de la manière suivante pour les neurones de la couche cachée :

Le même type de calcul est appliqué aux neurones de la couche cachée. Finalement, les paramètres du réseau sont modifiés par la formule suivante à l’itération i de l’apprentissage.

La quantité µi représente le pas du gradient ou pas d’apprentissage.

Cette méthode représente des inconvénients quand le gradient de la fonction de coût devient trop petit. Pour éviter ce problème, des techniques basées sur la méthode de Newton sont utilisées en considérant que le pas d’apprentissage est égal à l’inverse de la dérivée de second ordre de la fonction de coût. Parmi ces techniques, citons l’algorithme de Levenberg-Marquardt.

Avant tout apprentissage, toutes les variables sont normalisées et centrées pour éviter la surestimation de l’effet des variables ayant des grandeurs importantes. Seules les variables représentant la température de l’air, le débit et la quantité des précipitations sont retenues dans une première phase. La corrélation entre les valeurs calculées (prédites) et observées (expérimentales) est examinée au moyen d’une relation linéaire :

Tobs = a*Tcalc + b. T désignant la température de l’eau, a et b représentent la pente de la droite et la valeur de température à l’origine.

Les paramètres statistiques sont calculés à partir des relations ci-dessous :

r2 et s représentent respectivement le coefficient de corrélation et l’écart-type de l’erreur et n est le nombre d’exemples traités.

4. Réalisation du modèle neuronal pour la prédiction et la prévision de la températion de l’eau

Le modèle consiste en une série de trois couches de neurones. La première couche représente les variables retenues précédemment : la température de l’air, le débit du cours d’eau et la quantité des précipitations pendant la journée. La dernière couche représente la température de l’eau. Les paramètres du réseau (les valeurs des poids) sont initialisés à l’aide d’une fonction qui fournit des valeurs aléatoires. L’apprentissage des données, c’est-à-dire la minimisation de la fonction de coût, est réalisé à l’aide de l’algorithme de Lenvenberg-Marquardt qui donne de meilleurs résultats que les algorithmes usuels à pas constant, cités précédemment. La couche intermédiaire est composée de sept neurones. L’architecture retenue est donc de type 3‑7‑1. Pour un nombre de neurones dans la couche cachée supérieur à 7, le gain n’est pas significatif.

4.1 Généralisation par validation croisée

Notre objectif est de choisir un modèle capable d’apprendre les données qu’on lui propose, mais aussi qui offre une bonne généralisation à des fins de prédiction, en d’autres termes qui évite le surajustement lors d’un test. Pour ce faire, une des techniques utilisées consiste à décomposer l’ensemble de données en trois sous-ensembles. Le premier sert à l’apprentissage du modèle et le deuxième à ajuster le modèle et éviter le surapprentissage en effectuant un arrêt prématuré quand son erreur augmente lors de l’apprentissage. Le troisième sous-ensemble sert à évaluer les performances du modèle. Cette technique est assez lourde à réaliser et pose la problématique du choix du deuxième sous-ensemble pour le calibrage du modèle. Dans notre cas, nous avons utilisé une méthode basée sur la technique de régularisation bayésienne par modération des poids. Elle consiste à pénaliser les valeurs élevées des poids en modifiant la fonction de coût (Dreyfuset al., 2002). Cette technique force les paramètres (les poids) à ne pas prendre des valeurs trop élevées, et, par conséquent, à éviter le surajustement.

Nous avons testé le pouvoir prédictif du modèle neuronal par la technique de la validation croisée (leave-20%-out cross-validation method). Dans cette procédure, nous retirons 20 % de l’échantillon de base qui servira de test et le reste à l’apprentissage. Les données sont classées par ordre chronologique. Pour les cinq décompositions testées, le choix des valeurs initiales de poids n’a presque pas d’effet et la convergence de la fonction de coût est observée après 300 itérations en moyenne.

Les résultats présentés au tableau 2 sont plus que satisfaisants. Seuls 14 cas présentent une erreur supérieure à 3 °C et un seul cas avec une erreur supérieure à 4 °C. Il est important de signaler que ces paramètres statistiques proviennent d’une validation croisée dont le test porte sur l’ensemble des données. La figure 3 présente l’évolution de la température prédite en fonction de la température observée pour les cinq ensembles de test.

Tableau 2

Résultats statistiques de la validation croisée pour une architecture 3-7-1.

Statistical results of the cross validation for the architecture 3-7-1.

Coefficient de détermination R2

0,92

Déviation standard (RMSE, s)

1 °C

Histogramme de l’erreur (res = obs-calc)

‑1<res<1

‑2<res<2

3<res

 

70 %

93 %

0,9 %

-> Voir la liste des tableaux

Figure 3

Températures observées et calculées pour l’ensemble des données en validation croisée.

Observed and calculated temperatures for the whole set by means of the cross validation technique.

Températures observées et calculées pour l’ensemble des données en validation croisée.

-> Voir la liste des figures

La nature de la décomposition de l’ensemble des données permet d’évaluer la performance du modèle. Lorsque les ensembles du test sont composés à partir de données prélevées une par une à intervalles égaux, les paramètres statistiques sont nettement améliorés (R2 = 0,95, s = 0,8 °C), ceci est dû au fait que la température ne change pas de manière significative d’un jour à l’autre. La figure 4 présente l’erreur résiduelle pour l’ensemble des données en validation croisée.

Figure 4

Répartition des résidus pour l’ensemble des données en validation croisée.

Distribution of the residuals for the whole set by using cross validation.

Répartition des résidus pour l’ensemble des données en validation croisée.

-> Voir la liste des figures

L’analyse des données historiques de l’évolution de la température de l’eau en rivière montre que celle-ci ne varie pas beaucoup à court terme, c’est-à-dire sur un ou deux jours. La température de l’air ambiante est la variable explicative, dont l’effet est le plus important sur ce processus. Une simple analyse statistique entre ces deux variables montre que la température de l’eau dépend non seulement de la température de l’air ambiant du même jour, mais aussi des jours précédents. Dans la perspective d’augmenter la capacité d’adaptabilité du modèle neuronal, nous avons considéré les températures de l’air ambiant des jours précédents comme des variables explicatives supplémentaires. Pour une meilleure flexibilité du réseau, seuls les deux jours précédents sont pris en considération. Le nombre de neurones dans la couche intermédiaire est maintenu égal à 7. Les paramètres statistiques de la validation croisée sont nettement meilleurs. Le coefficient de détermination est passé à 0,96, l’écart-type de l’erreur est de l’ordre de 0,7 °C. Comparativement à la figure 4, la figure 5 montre une nette amélioration de la répartition des résidus, seul un point présente une erreur entre -1 et -1,5.

Figure 5

Résidus pour l’ensemble des données en validation croisée en considérant les valeurs des températures de l’air aux jours précédents.

Cross validation residuals for the whole set by considering the air temperature for the previous days.

Résidus pour l’ensemble des données en validation croisée en considérant les valeurs des températures de l’air aux jours précédents.

-> Voir la liste des figures

D’une manière générale, les résultats sont très satisfaisants et justifient le recours à l’approche par les réseaux de neurones dans la prédiction de la température de l’eau en rivière.

4.3 Modélisation dynamique par les réseaux de neurones

Les données historiques consistent en une série chronologique de l’évolution de la température de l’eau en fonction du temps. Les modélisations s’inscrivant dans la stratégie de Box et Jenkins (Box et Jenkins, 1976) AR, MA, ARMA, ARIMA et SARIMA se prêtent bien à ce type de série, puisqu’elles considèrent que les données peuvent être fortement corrélées, comme le sont habituellement les séries successives d’une manière générale. L’observation est liée d’une manière ou d’une autre aux observations précédentes. Dans notre cas, un modèle de type ARMA se réduit tout simplement à une relation linéaire. La valeur de la température de l’eau ŷ(j) est reliée aux valeurs précédentes ŷ(j‑k), les variables explicatives précédemment retenues (température de l’air, précipitation et débit) et les erreurs de chaque estimation e (t‑k) suivant la relation ci-dessous :

ak, bl et cm représentent les nombres de jours précédents pris en considération et na, nb et nc représentent le nombre des jours précédents pour chaque variable pris en compte dans le modèle. La température de l’air, les précipitations et le débit du cours d’eau sont représentés par l’élément x(t–l), les valeurs estimées de la température de l’eau sont représentées par ŷ(t‑k). En général, la régression linéaire ne convient pas puisque les données sont supposées dépendantes, leur contribution étant alors arbitraire et ne reflétant pas la réalité de la relation (Zhang, 2001).

Dans la première partie, nous nous sommes intéressés à la prédiction de la température de l’eau en rivière sans faire référence à un temps spécifique. Nous nous proposons dans cette partie de tenir compte de la chronologie des données. Il s’agit d’estimer la température de l’eau des jours à venir (j + i, i = 1,2,..). Notons que l’architecture du modèle neuronal se complique notoirement à chaque étape, puisqu’avec notre technique, le nombre d’entrées augmente. Afin de préserver la parcimonie du modèle neuronal, nous avons jugé utile de réduire le nombre de paramètres d’entrée, en éliminant les deux variables représentant le cumul des précipitations et le débit ainsi que les variables représentant les températures de l’air ambiant des jours précédents. Les paramètres statistiques d’évaluation des modèles restent les mêmes.

Pour réaliser une multitude d’estimations par les réseaux de neurones, nous utilisons deux méthodes. La première procède par itération, elle utilise la valeur estimée du jour j pour estimer la valeur de la température de l’eau du jour j + 1. Nous réalisons ainsi autant de relations que de valeurs à estimer, en utilisant les relations suivantes :

x(j) représente la température de l’air du jour j, les fonctions f0, f1, f2…..f7 sont générées par les modèles neuronaux de chaque étape. La figure 6a illustre cet ensemble de fonctions.

Figure 6

(a) Architecture neuronale pour la prévision par itération, (b) Architecture neuronale pour la prévision par bloc.

Architecture for the forecast with iteration (a) and with blocks (b).

(a)

(b)

-> Voir la liste des figures

La deuxième méthode est beaucoup plus simple, elle consiste à réaliser une estimation de tous les jours en une seule fois. Une seule fonction est utilisée (Figure 6b) :

Le nombre de neurones dans la couche intermédiaire a été fixé à sept, la technique de minimisation de la fonction de coût étant celle de Levenberg-Marquardth. Pour tester le modèle, une simple décomposition de la base de données initiale a été effectuée. En effet, l’ensemble des trois premières années a servi à l’apprentissage du modèle neuronal, les deux dernières années servant au test. Les paramètres statistiques sont résumés au tableau 3.

Tableau 3

Coefficient de détermination et erreur standard. (a) Prévision en procédant par itération, (b) prévision par bloc.

Coefficient of determination and standard error: (a) forecast by iteration, (b) forecast by block.

Jour j

0

1

2

3

4

5

6

7

(a) R2 (s)

0,99(0,39)

0,95(0,83)

0,89(1,20)

0,86(1,36)

0,84(1,43)

0,83(1,47)

0,83(1,50)

0,82(1,53)

erreur_max

1,38

3,37

3,99

4,90

5,41

6,14

5,53

5,50

nb (res>2 °C)

0 cas

12

47

64

72

79

88

86

(b) R2 (s)

0,99(0,40)

0,95(0,82)

0.89(1,16)

0,86(1,34)

0,84(1,43)

0,83(1,47)

0,82(1,51)

0,84(1,53)

erreur_max

1,38

3,38

4,10

5,08

5,00

5,68

5,09

5,03

nb (res>2 °C)

15 cas

15

48

71

76

75

82

87

res =Température calculée - Température observée

nb (res >2°C) : nombre d’observations avec une erreur supérieure à 2

R2 et s sont respectivement le coefficient de détermination et l’écart-type

-> Voir la liste des tableaux

En général, le coefficient de détermination diminue subitement dès que l’on prédit au-delà du premier jour. Néanmoins, les valeurs restent acceptables. Les deux méthodes conduisent à des résultats semblables. L’erreur standard augmente aussi. La figure 7 présente la dispersion de l’erreur pour les quatre jours avec la méthode (b). Nous remarquons qu’elle est presque la même d’un jour à l’autre. Cette dispersion ne reflète pas la suite de la chronologie des événements. Les données correspondant aux intervalles, allant du mois de mai jusqu’au mois de décembre 2001 et du mois de janvier jusqu’au mois d’août 2002, sont manquantes. Les valeurs des résiduels sont faibles au cours de la période chaude. En effet, la température de l’eau en rivière varie faiblement en période chaude.

Figure 7

Résidus de la température de l’eau pour les quatre jours (méthode b).

Residuals of the temperature for the four days (method b).

Résidus de la température de l’eau pour les quatre jours (méthode b).

-> Voir la liste des figures

Dans la bibliographie, aucune étude n’a été réalisée pour expliquer la différence entre les deux méthodes. Les avis sont contradictoires quant aux performances d’une méthode par rapport à l’autre. Dans notre cas, la méthode (a) ne présente aucun avantage, malgré le fait qu’on procède par itération. La brusque diminution des paramètres statistiques est due principalement à la propagation de l’erreur, puisque nous nous servons de la valeur calculée de la température de l’eau du jour j pour réaliser celle du jour j + 1. De surcroît, la nature du modèle connexionniste amplifie cette erreur, en raison de la liaison qui existe entre les paramètres.

5. Conclusion

Dans ce travail, nous avons montré que le modèle connexionniste basé sur les réseaux de neurones est approprié à la modélisation de la température de l’eau d’un ruisseau. La température de l’air ambiante et les températures de l’eau des jours précédents suffisent pour construire ce modèle.

D’autre part, nous avons utilisé deux techniques pour réaliser des prévisions dans le futur. Elles ont conduit à des résultats identiques jugés très satisfaisants. La technique qui procède par itération est prometteuse si l’amplification de la propagation de l’erreur est contrôlée.

Dans la pratique, le modèle serait nettement meilleur si l’on disposait d’un modèle de prédiction de la température de l’air des jours j + 1, j + 2…..j + 6. Les fonctions qui permettent de prévoir la température de l’eau auraient la forme suivante :

x̂(j) représente le modèle de prédiction de la température de l’air.