Traitement et analyse de séries chronologiques continues de ...

Partie 6 – Chapitre 17 : Test des modèles Accumulation-Erosion-Transfert
17.16). D’autres essais (non présentés ici) ont également mis en évidence la sensibilité du
calage à la valeur du paramètre Ms ini . Pour éviter ce problème, nous avons fixé dans la suite des
tests la valeur de Ms ini égale à la valeur de Ms max . Nous faisons donc l’hypothèse qu’en début de
simulation l’intégralité du stock de sédiments est disponible sur la surface du bassin.
12
ks
Ms ini
Statistique de Gelman et Rubin
10
8
6
4
Ms max
ARRA
a
Kq1
Kq2
no q
std 0
2
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Nombre d'évaluation du DREAM
x 10 4
Figure 17.16. Diagnostic de convergence des paramètres (critère de Gelman et Rubin), pour les paramètres de
calage (modèle qualité et modèle d’erreur) ; résultats pour le modèle qualité (niveau de complexité NC1), site de
Chassieu, 1 ère itération ; les valeurs du critère pour les paramètres K Q1 et K Q2 ne parviennent pas en dessous de la
valeur seuil
Pour la deuxième itération, nous avons adopté des bornes a priori plus appropriées pour les
paramètres K Q1 et K Q2 (cf. Tableau 17.8). Les distributions de ces derniers sont cette fois mieux
identifiées (Figure 17.17) mais montrent néanmoins une corrélation linéaire significative. Pour
le paramètre k s , la distribution a posteriori est nettement tronquée à la valeur 2, suggérant que la
distribution a priori pour ce paramètre n’est pas assez large. Ce paramètre est relatif à
l’accumulation des sédiments de temps sec. Cela suggère donc que lors des épisodes de temps
sec le modèle sous-estime la vitesse de reconstitution du stock de sédiments. Une corrélation
significative est également observée entre les paramètres ARRA et Ms max , ce qui n’est pas
vraiment étonnant étant donné que par construction ils peuvent se compenser. Enfin, si le
paramètre to Q du réservoir linéaire double est bien identifié, sa valeur la plus probable est faible
(proche de 1 minute).
301