Réseau de Koopman gaussien conditionnel : prévision et assimilation de données dans un seul modèle ✛ ✎

Original: "CGKN: A Deep Learning Framework for Modeling Complex Dynamical Systems and Efficient Data Assimilation"
arXiv:2410.20072 · 2024-10-26 · arXiv non-exclusive · analyse indépendante · 4 min
Un réseau de neurones qui prédit un système chaotique tout en assimilant instantanément les observations.

Contexte

Les systèmes dynamiques complexes, du climat à la neurobiologie, sont souvent décrits par des équations trop coûteuses à résoudre, voire inconnues. Les réseaux de neurones profonds ont appris à prédire ces systèmes à partir de données, mais leur boîte noire entrave l'assimilation de données, une technique clé qui intègre les observations dans le modèle et estime les variables cachées. Les méthodes d'ensemble d'assimilation nécessitent des milliers d'exécutions du modèle et souffrent d'erreurs d'échantillonnage. Ce travail se situe à l'intersection de l'apprentissage automatique scientifique et de la théorie de Koopman : les approches voisines soit linéarisent complètement la dynamique, perdant la non-linéarité, soit restent des boîtes noires avec une assimilation inefficace.

Méthodes

CGKN construit une équation différentielle stochastique neuronale en trois étapes. D'abord, un encodeur projette les variables non observées dans un espace latent où la dynamique devient conditionnellement linéaire : étant donné les variables observées, les variables latentes suivent une équation linéaire avec un bruit gaussien. Ensuite, des réseaux de neurones approchent les quatre fonctions de cette équation gaussienne conditionnelle. Enfin, un décodeur ramène les variables latentes dans l'espace d'origine. L'idée clé est une généralisation de la théorie de Koopman : au lieu d'une linéarisation complète, on conserve la non-linéarité par rapport aux variables observées, mais on introduit une linéarité conditionnelle. Cela donne des formules analytiques fermées pour la moyenne et la covariance de l'état latent, qui sont intégrées par un solveur d'équations différentielles ordinaires classique. L'assimilation de données devient analytique et est incluse dans la fonction de perte lors de l'entraînement, éliminant ainsi le besoin de réglage empirique des ensembles.

Résultats

Sur les trois systèmes de test, CGKN a montré une précision de prévision comparable à celle du modèle vrai et des meilleures boîtes noires, tout en effectuant l'assimilation de données de manière analytique. Pour l'équation de Burgers–Sivashinsky, l'erreur de prévision était de 0.284 contre 0.246 pour le modèle vrai, et l'erreur d'assimilation de 0.728 contre 0.695. Dans le système de Lorenz-96 à 40 variables, CGKN a atteint une erreur de prévision de 0.064 et d'assimilation de 0.095, tandis qu'un réseau de Koopman standard donnait respectivement 0.306 et 0.273. Pour El Niño, le modèle a prédit la température de surface de l'océan avec une corrélation supérieure à 0.5 sur un horizon de 12 mois, surpassant la prévision persistante. Dans tous les cas, CGKN a surpassé le réseau de Koopman standard et la régression avec bibliothèque gaussienne conditionnelle, et dans la tâche d'assimilation de données, il a également surpassé la boîte noire DNN, pour laquelle l'assimilation analytique est impossible.

Implications

Ce travail se distingue des approches voisines en ce qu'il ne sacrifie pas la non-linéarité au profit d'une assimilation analytique. Le réseau de Koopman standard linéarise complètement la dynamique latente et perd la contribution des variables observées, tandis que la régression CG-Reg est limitée à une bibliothèque de fonctions prédéfinie. CGKN apprend des transformations non linéaires par des réseaux de neurones et conserve la linéarité conditionnelle, offrant à la fois expressivité et efficacité. C'est la première démonstration qu'une structure gaussienne conditionnelle peut être intégrée dans l'apprentissage profond et utilisée pour l'assimilation analytique de données dans des systèmes turbulents fortement non linéaires.

Développements futurs

Ce thème pourrait évoluer vers une évaluation quantitative rigoureuse de l'incertitude pour les variables non observées d'origine. Actuellement, l'incertitude est estimée par post-traitement via l'analyse des résidus, mais le décodeur non linéaire rend difficile la projection précise de la covariance depuis l'espace latent. Des méthodes d'optimisation multi-objectifs pourraient équilibrer la qualité de la transformation, de la prévision et de l'assimilation de données, ainsi que l'extension à des problèmes de contrôle de systèmes partiellement observés et à des problèmes inverses.

Impact

La méthode est applicable en climatologie, géophysique, science des matériaux et neurobiologie, où des prévisions rapides et l'estimation d'états cachés à partir de données incomplètes sont nécessaires.

Prochaines étapes

La prochaine étape consistera à développer des méthodes efficaces d'estimation de l'incertitude pour les variables d'origine, par exemple via le jacobien du décodeur ou la transformée unscented pour les cas faiblement non linéaires.

Problèmes ouverts clés

Ce travail est directement lié au problème de longue date de l'assimilation de données dans les systèmes non linéaires : comment concilier l'expressivité des modèles profonds avec l'efficacité computationnelle. Il propose une voie à travers des contraintes structurelles, évitant la malédiction de la dimensionnalité caractéristique des méthodes d'ensemble.

Le regard de la machine du savoir

Cette section s'adresse à l'auteur du travail et à ses collègues. Les travaux voisins ont été trouvés par une recherche sémantique dans notre archive, et non par correspondance de mots. Nous n'évaluons pas : nous montrons le champ autour.

Ce que fait ce travail

Les auteurs proposent le modèle CGKN, qui transforme un système dynamique non linéaire en une forme conditionnellement gaussienne à l'aide de réseaux de neurones, permettant ainsi d'obtenir des formules analytiques pour l'assimilation de données. Le modèle est testé sur trois systèmes turbulents : l'équation stochastique de Burgers–Sivashinsky, le système de Lorenz-96 et El Niño. Les résultats montrent une précision de prévision au niveau des meilleures boîtes noires, avec une accélération radicale de l'assimilation de données.

Où est sa force

La force clé de ce travail réside dans la préservation de la non-linéarité pour les variables observées tout en gardant une dynamique cachée conditionnellement linéaire, ce qui offre à la fois l'expressivité des réseaux de neurones et l'efficacité analytique de l'assimilation. C'est la première démonstration de l'intégration d'une structure conditionnellement gaussienne dans l'apprentissage profond pour des systèmes turbulents fortement non linéaires.

Vers où avancer

  • Il pourrait être intéressant d'examiner comment CGKN se comporte dans des problèmes où les données observées sont fortement bruitées ou incomplètes, par exemple dans les mesures astrophysiques, où les réseaux de neurones sont déjà utilisés pour reconstruire des signaux à partir d'observations limitées. À proximité se trouve l'expérience du réseau de neurones sculpteur, qui reconstruit une lentille gravitationnelle à partir de quatre points : là aussi, il faut composer avec une information incomplète et des paramètres cachés.appui : Sculpteur neuronal : comment interpréter une lentille gravitationnelle à partir de quatre points arXiv:2607.22812v1
  • Une suite possible est l'adaptation de CGKN à des systèmes présentant des fluctuations chaotiques rapides, où l'assimilation est importante mais aussi la prévision d'événements rares. Dans les archives, il existe un travail sur la domestication du chaos, où des algorithmes prévoient l'avenir dans des séries chaotiques ; comparer les approches pourrait suggérer des améliorations pour la prévision à long terme.appui : Dompter le chaos : quand les algorithmes prédisent l'avenir arXiv:2607.11344
  • Une autre direction est l'application de CGKN à des problèmes avec des données distribuées spatialement, par exemple en géophysique ou en science des matériaux, où les réseaux de neurones ont déjà appris à extraire des propriétés physiques à partir de signaux complexes. L'expérience avec le réseau de neurones qui « entend » l'élasticité des roches pourrait être utile lors de l'adaptation du modèle à des champs multidimensionnels.appui : La mélodie de la pierre : quand un réseau de neurones « entend » l'élasticité de la roche arXiv:2606.25275

À essayer

  • Tester CGKN sur des données avec des valeurs manquantes : entraîner le modèle en masquant une partie des variables observées et comparer la qualité de l'assimilation avec celle obtenue sur des données complètes.
  • Utiliser l'augmentation de données en ajoutant du bruit aux variables observées pendant l'entraînement afin d'améliorer la robustesse de l'assimilation analytique.
  • Essayer d'appliquer CGKN à la prévision à court terme de signaux chaotiques, où la vitesse de calcul est cruciale, et comparer avec les méthodes d'ensemble.

Pourquoi cela peut compter

Si CGKN peut être étendu à une classe plus large de problèmes avec des données incomplètes et bruitées, cela pourrait fournir un outil universel pour une assimilation rapide de données dans les sciences de la Terre, l'astrophysique et l'ingénierie, où les méthodes d'ensemble coûteuses dominent actuellement.

L'équation de Burgers–Sivashinsky, utilisée dans les tests, décrit non seulement la turbulence, mais aussi l'instabilité du front de flamme, et sa version stochastique génère des événements extrêmes semblables à des rafales soudaines dans un écoulement turbulent.

Chiffres clés

  • Erreur de prévision pour l'équation de Burgers–Sivashinsky: 0.284 (CGKN) contre 0.246 (modèle vrai)
  • Erreur d'assimilation de données pour Lorenz-96: 0.095 (CGKN) contre 0.273 (réseau de Koopman standard)
  • Corrélation de prévision El Niño à 12 mois: supérieure à 0.5
  • Dimension de l'espace latent pour Lorenz-96: 120
  • Longueur de la série d'entraînement pour El Niño: 2000 ans de temps modèle
Notions
réseau neuronal de graphe généralisation processus gaussien coupled general circulation model modélisation physique par réseaux de neurones xgboost koopman operator réseaux antagonistes génératifs filtre de kalman nonlinear laplacian spectral analysis systèmes complexes apprentissage contrastif regroupement régression par processus gaussien grands modèles de langage dynamique des systèmes complexes systèmes dynamiques modèle génératif apprentissage flots normalisants