Reconstruction de l'onde non-linéaire à une bouée à partir des données de rivage et applications au problème inverse du tsunami pour une bathymétrie à pentes par morceaux

Informations fondamentales

• ID de l'article: 2510.14177
• Titre: Reconstruction of the non-linear wave at a buoy from shoreline data and applications to the tsunami inverse problem for piece-wise sloping bathymetry
• Auteurs: Oleksandr Bobrovnikov¹, Madison Jones², Shriya Prasanna³, Josiah Smith², Alexei Rybkin¹, Efim Pelinovsky⁴'⁵
• Affiliations: ¹University of Alaska Fairbanks, ²University of Colorado Boulder, ³University of Washington, ⁴HSE University, ⁵Institute of Applied Physics
• Classification: math.AP (Analyse mathématique appliquée)
• Date de publication: 17 octobre 2025
• Lien de l'article: https://arxiv.org/abs/2510.14177

Résumé

Cet article traite d'un problème inverse important : étant donné les données de remontée d'une onde de tsunami, peut-on récupérer sa forme initiale ? La recherche est menée dans le cadre des équations non-linéaires d'eau peu profonde, un modèle largement utilisé pour l'étude de la propagation et de l'inondation des tsunamis. Les auteurs considèrent une topographie de fond marin à pente finie et démontrent qu'il est possible de récupérer les conditions aux limites (déplacement du niveau d'eau et vitesse) à une bouée virtuelle à partir des données de rivage. L'article discute également de l'assemblage des équations d'eau peu profonde avec les équations de Boussinesq dans une topographie de fond marin à pentes par morceaux plus complexe, afin de récupérer les conditions initiales et d'introduire les effets dispersifs dans le modèle.

Contexte et motivation de la recherche

Importance du problème

1. Menaces de catastrophes: Le séisme de 1964 en Alaska, le tsunami de l'océan Indien de 2004, le grand tremblement de terre du Japon oriental de 2011 et d'autres catastrophes majeures de tsunami ont causé d'énormes pertes en vies humaines et dégâts économiques, soulignant l'importance de l'alerte précoce aux tsunamis et de la réduction des risques.
2. Besoins des systèmes d'alerte: Avec le développement du système DART (Deep-ocean Assessment and Reporting of Tsunamis), un cadre théorique est nécessaire pour guider le déploiement des bouées afin d'améliorer la précision et la rapidité de l'alerte aux tsunamis.
3. Défis de la modélisation mathématique: La propagation des tsunamis implique des effets non-linéaires complexes, particulièrement dans les zones côtières, nécessitant des modèles mathématiques précis pour décrire le comportement des ondes.

Limitations des méthodes existantes

Les méthodes de réversion de tsunami existantes comprennent principalement quatre approches:

1. Méthodes combinant données sismiques et de forme d'onde
2. Méthodes de reconstruction utilisant uniquement les données de forme d'onde
3. Méthodes de réversion utilisant la distribution des sédiments
4. Méthodes de réversion basées sur les données d'inondation

Cependant, ces méthodes présentent des limitations dans le traitement d'une topographie de fond marin à pente finie, particulièrement en considérant les effets combinés des non-linéarités et des effets dispersifs.

Motivation de la recherche

La motivation de cet article est d'étendre les travaux antérieurs des auteurs dans le cas d'une topographie de fond marin à pente infinie, en développant une théorie de réversion applicable à une topographie de fond marin à pente finie plus réaliste, et en fournissant un cadre théorique pour guider le déploiement des bouées.

Contributions principales

1. Extension théorique: Extension de la théorie de réversion des tsunamis du fond marin à pente infinie à une topographie de fond marin à pente finie, fournissant une solution analytique complète.
2. Développement d'algorithmes: Proposition d'un processus algorithmique complet pour récupérer les conditions aux limites à une bouée virtuelle à partir des données de remontée du rivage.
3. Assemblage de modèles: Développement d'une méthode pour assembler les équations non-linéaires d'eau peu profonde (NSWE) avec les équations de Boussinesq dans une topographie de fond marin à pentes par morceaux, considérant simultanément les effets non-linéaires et dispersifs.
4. Vérification numérique: Validation de l'efficacité de la méthode par plusieurs scénarios physiquement réalistes, incluant différents types de formes d'onde initiales.
5. Efficacité computationnelle: Fourniture d'une implémentation algorithmique efficace avec une complexité O(N³), supportant le calcul rapide des caractéristiques des tsunamis.

Détails de la méthode

Définition de la tâche

Entrée: Données de remontée du rivage R(t) Sortie: Déplacement du niveau d'eau η(L,t) et vitesse u(L,t) à la bouée virtuelle Contraintes: Condition d'onde non-déferlante, validité de la transformation de Carrier-Greenspan

Cadre mathématique fondamental

1. Système d'équations non-linéaires d'eau peu profonde

En unités sans dimension, le système (1) est:

∂tη + u∂x(x + η) + (x + η)∂xu = 0  (conservation de la masse)
∂tu + u∂xu + ∂xη = 0              (conservation de la quantité de mouvement)

2. Transformation de Carrier-Greenspan

La transformation de linéarisation clé (5):

φ(σ,τ) = u(x,t)
σ = x + η(x,t)
ψ(σ,τ) = η(x,t) + u²(x,t)/2
τ = t - u(x,t)

Cette transformation convertit le système non-linéaire (1) en système linéaire (6):

∂τψ + σ∂σφ + φ = 0
∂τφ + ∂σψ = 0

3. Solution du problème direct

En introduisant les nouvelles variables ρ² = σ/σL et la fonction θ(ρ,τ) = ψ(σ,τ) - ψb(τ), en utilisant la décomposition de Fourier:

θ(τ,ρ) = Σ cn(τ)J₀(jnρ)

où J₀ est la fonction de Bessel de première espèce et jn est le n-ième zéro de J₀.

4. Formules centrales de l'algorithme de réversion

Équation du rivage (20):

ψsh(τ) = Σ (bn/√an) ∫₀τ sin(√an(τ-ξ))ψ''b(ξ)dξ + ψb(τ)

Formule de réversion (24):

ψb(τ) = L⁻¹[Lψsh/(s²Σ bn/(an+s²) + 1)⁻¹]

Processus complet de l'algorithme de réversion

1. Conversion des données: Utilisation de la CGT pour convertir R(t) en ψsh(τ)
2. Récupération des limites: Calcul de ψb(τ) à partir de ψsh(τ) via la formule (24)
3. Reconstruction du champ: Utilisation des formules (15,16) pour calculer ψ(σ,τ) et φ(σ,τ) dans tout le domaine
4. Transformation inverse: Exécution de la CGT inverse pour récupérer les conditions aux limites physiques u(L,t) et η(L,t)

Points d'innovation technique

1. Méthode de solution analytique: Contrairement aux méthodes numériques, fournit une solution de réversion complètement analytique, avec une efficacité computationnelle élevée et une précision contrôlable.
2. Technique de transformation de Laplace: Application astucieuse de la transformation de Laplace et de la transformée de Fourier rapide inverse (iFFT) pour améliorer l'efficacité computationnelle.
3. Stratégie d'assemblage de modèles: Assemblage innovant des équations NSWE et de Boussinesq, traitant simultanément les effets non-linéaires et dispersifs.

Configuration expérimentale

Scénarios de vérification numérique

Les auteurs ont conçu quatre scénarios de déplacement initial différents pour la vérification:

1. Forme d'onde bimodale: η(x,0) = 0.005 sech(x-16) - 0.003e^(-(x-13)²)
2. Forme d'onde soliton: η(x,0) = 0.005 cosh⁻²(2x-6)
3. Forme d'onde gaussienne: η(x,0) = 0.005e^(-(x-7)²)
4. Forme d'onde multi-pics: η(x,0) = 0.005e^(-2(x-6)²) + 0.003e^(-(x-10)²)

Indicateurs d'évaluation

• Comparaison de précision: Comparaison du ψb(τ) récupéré avec la valeur exacte ψᵉb(τ)
• Récupération des conditions aux limites: Vérification de la précision de récupération de u(L,t) et η(L,t)
• Efficacité computationnelle: Mesure du temps d'exécution et de la complexité de l'algorithme

Détails d'implémentation

• Configuration des limites: L = 1 (unités sans dimension)
• Nombre de points de données: N = 1500 points de données
• Paramètre de troncature: Zéros de la fonction de Bessel tronqués à n = 500
• Plateforme de calcul: Environnement Google Colab

Résultats expérimentaux

Résultats principaux

1. Vérification de la précision de réversion

Les figures 4 et 5 montrent que pour tous les quatre scénarios de test, la cohérence entre les conditions aux limites récupérées ψb(τ) et φb(τ) et la solution exacte est excellente, validant l'efficacité de l'algorithme.

2. Performance computationnelle

• Temps d'exécution: Le calcul de 1500 points de données prend 76 secondes sur Google Colab
• Vérification de complexité: La complexité temporelle est conforme à l'analyse théorique O(N³)
• Résultat d'ajustement: La pente de l'ajustement des moindres carrés entre le temps et le nombre de points de données est 2,82, proche de la valeur théorique 3

3. Récupération des conditions aux limites physiques

La figure 7 montre les conditions aux limites physiques récupérées u(L,t) et η(L,t) pour les quatre scénarios, indiquant que l'algorithme peut reconstruire avec précision les informations du champ d'onde à la bouée virtuelle.

Expériences d'assemblage des équations de Boussinesq

1. Vérification de la solution de double soliton de Hirota

Utilisant la solution de double soliton de Hirota comme condition aux limites, validation réussie de la méthode d'assemblage des équations NSWE et de Boussinesq. Paramètres expérimentaux:

• q₁ = 0.1, q₂ = √0.1
• Différentes combinaisons de paramètres temporels t₁, t₂

2. Reconstruction des conditions initiales

La figure 13 montre les résultats de reconstruction des conditions initiales au moment de l'événement de tsunami (tevent = -100s):

• Position exacte des pics: x₁ = 119.55, x₂ = 123.63
• Position des pics récupérés: x₁ = 119.50, x₂ = 123.64
• Erreur minimale, validant la précision de la méthode

Comparaison avec les équations linéaires d'eau peu profonde

La figure 16 montre la comparaison entre le modèle de Boussinesq et le modèle linéaire d'eau peu profonde sous les mêmes conditions aux limites, mettant en évidence l'impact des effets dispersifs.

Travaux connexes

Directions de recherche en réversion de tsunami

1. Réversion conjointe sismique-forme d'onde: Yokota et al. (2011), Yue et al. (2014)
2. Réversion de données de forme d'onde pure: Ho et al. (2019), Yamanaka and Tanioka (2024)
3. Réversion de sédiments: Tang and Weiss (2015), Ioki and Tanioka (2016)
4. Réversion de données d'inondation: Lee et al. (2021), Piatanesi et al. (1996)

Développement des méthodes numériques

• Différences finies/Éléments finis: Kounadis and Dougalis (2020), Xing (2017)
• Réseaux de neurones informés par la physique: Brecht et al. (2025)
• Modélisation de topographie par morceaux: Synolakis (1987), Varsoliwala and Singh (2021)

Avantages relatifs de cet article

Comparé aux travaux existants, cet article fournit une méthode de réversion complètement analytique, évitant les erreurs cumulatives des méthodes numériques, et traite systématiquement pour la première fois le problème de réversion pour une topographie à pente finie.

Conclusions et discussion

Conclusions principales

1. Complétude théorique: Extension réussie de la théorie de réversion des tsunamis à une topographie de fond marin à pente finie, fournissant une solution analytique complète.
2. Valeur pratique: L'algorithme peut être utilisé pour guider le déploiement optimal des bouées de systèmes d'alerte aux tsunamis tels que DART, améliorant l'efficacité de l'alerte.
3. Intégration de modèles: Réalisation réussie de l'assemblage des équations NSWE et de Boussinesq, fournissant une nouvelle voie pour la modélisation des tsunamis dans des topographies plus complexes.

Limitations

1. Hypothèse d'onde non-déferlante: L'algorithme n'est applicable que dans le cas d'ondes non-déferlantes; la méthode échoue lorsque le déterminant jacobien est nul.
2. Limitation du modèle unidimensionnel: Actuellement, seul le cas unidimensionnel est considéré, tandis que la propagation réelle des tsunamis est un phénomène tridimensionnel.
3. Simplification de la topographie: Bien que les topographies par morceaux soient considérées, elles restent très simplifiées par rapport à la topographie réelle du fond marin.
4. Portée des effets dispersifs: L'applicabilité des équations de Boussinesq reste controversée dans certains cas.

Directions futures

1. Extension multidimensionnelle: Développement de théories de réversion bidimensionnelles et tridimensionnelles
2. Topographie complexe: Traitement du problème de réversion pour une topographie de fond marin à pentes multiples
3. Critères d'onde déferlante: Établissement de critères clairs pour les données de rivage acceptables
4. Applications pratiques: Combinaison avec les données d'observation réelles de tsunamis pour la vérification

Évaluation approfondie

Points forts

1. Rigueur mathématique: Basée sur une théorie mathématique stricte, fournissant une solution analytique complète, évitant l'accumulation d'erreurs des méthodes numériques.
2. Efficacité computationnelle: La complexité O(N³) est acceptable compte tenu des exigences de temps réel de l'alerte aux tsunamis, avec des avantages de vitesse de calcul particulièrement évidents pour les modèles unidimensionnels.
3. Innovation forte: Première résolution systématique du problème de réversion des tsunamis pour une topographie à pente finie, avec assemblage innovant de différentes équations.
4. Vérification suffisante: Vérification numérique par plusieurs scénarios de forme d'onde, démontrant la stabilité et la précision de la méthode.

Insuffisances

1. Limitations de la portée d'application: L'hypothèse d'onde non-déferlante limite l'application de la méthode dans les cas fortement non-linéaires.
2. Simplification du modèle de topographie: La topographie réelle du fond marin est bien plus complexe que les pentes par morceaux, l'efficacité pratique de la méthode reste à vérifier.
3. Manque de vérification avec données réelles: Seules les données synthétiques sont utilisées pour la vérification, sans comparaison avec les observations réelles de tsunamis.
4. Compatibilité des conditions aux limites: Bien que les conditions de compatibilité soient mentionnées, leur impact dans les applications pratiques n'est pas discuté en détail.

Impact

1. Contribution théorique: Fournit une base mathématique importante pour la théorie de réversion des tsunamis, particulièrement dans les méthodes analytiques.
2. Valeur pratique: Fournit des conseils théoriques pour le déploiement de bouées dans les systèmes d'alerte aux tsunamis, avec une importance significative pour la réduction des risques de catastrophes.
3. Signification méthodologique: L'application innovante de la transformation de Carrier-Greenspan fournit des références pour d'autres problèmes d'ondes non-linéaires.

Scénarios d'application

1. Conception de systèmes d'alerte aux tsunamis: Guidage du déploiement optimal des équipements de surveillance tels que les bouées DART
2. Reconstruction de tsunamis historiques: Reconstruction d'événements historiques basée sur les données de sédiments de paléo-tsunamis
3. Recherche théorique: Fourniture d'outils mathématiques pour l'étude des mécanismes de propagation des tsunamis
4. Applications d'ingénierie: Références pour la conception de travaux de protection côtière

Références

L'article cite une littérature riche et pertinente, comprenant principalement:

• Carrier and Greenspan (1958): Travail original de la transformation CGT
• Rybkin et al. (2021, 2023, 2024, 2025): Série de travaux de l'équipe d'auteurs
• Synolakis (1987): Travail classique de modélisation de topographie par morceaux
• Hirota (1973): Contribution importante aux solutions soliton des équations de Boussinesq

Évaluation globale: Cet article est un travail de haute qualité en physique mathématique, apportant des contributions importantes à la théorie de réversion des tsunamis. La méthode possède une base mathématique stricte et une bonne efficacité computationnelle, avec une valeur théorique et pratique importante pour l'alerte aux tsunamis et la réduction des risques de catastrophes. Bien qu'il existe certaines limitations, il jette une base solide pour le développement ultérieur du domaine.