Fractional-Order Viscoelastic Modeling of Seismic P-, S-, and Rayleigh-Wave Attenuation: Marmara Seismotectonic Context, Fault Geometry, and a Reproducible Strong-Motion Framework

Bien en dessous de la surface de la Terre, le sol ne se comporte pas comme un simple bloc de roche solide. Lorsqu'un séisme survient, l'énergie se propage vers l'extérieur sous forme d'ondes, mais à mesure que ces ondes traversent la croûte, elles perdent de l'énergie et changent de forme. Ce processus, connu sous le nom d'atténuation, est la raison pour laquelle un sol qui tremble est ressenti différemment à une certaine distance de la source qu'immédiatement à côté de la rupture. Pendant des décennies, les scientifiques ont tenté de décrire cette perte d'énergie en utilisant les règles standards de la physique, en traitant la croûte terrestre soit comme un ressort parfaitement élastique, soit comme un fluide visqueux simple. Cependant, le monde réel est souvent plus complexe que ces deux extrêmes le permettent. Le sol est un mélange désordonné de roche, de fluide et de fissures qui se souvient de ses mouvements passés d'une manière que la physique standard peine à capturer. C'est ici que une branche des mathématiques appelée calcul fractionnaire devient utile. Au lieu d'utiliser des nombres entiers pour décrire la réaction des matériaux, le calcul fractionnaire permet d'utiliser des valeurs intermédiaires, offrant un moyen de modéliser la « mémoire » de la roche avec beaucoup moins d'hypothèses.

Dans une nouvelle étude axée sur la mer de Marmara, une région d'activité géologique intense située juste au sud d'Istanbul, les chercheurs Taylan Demir et Atakan Koçyiğit ont construit un cadre mathématique plus précis pour décrire la façon dont les ondes sismiques s'estompent. La mer de Marmara repose sur la faille nord-anatolienne, une immense fracture dans la croûte terrestre où les plaques tectoniques glissent l'une contre l'autre. Comme cette faille passe directement sous une métropole majeure, comprendre exactement comment les ondes sismiques voyagent et s'affaiblissent dans cette zone spécifique est crucial pour concevoir des bâtiments plus sûrs et prédire l'intensité des secousses futures. Les chercheurs ont cherché à tester si un modèle fractionnaire, qui rend compte du comportement complexe et dépendant de l'historique de la roche, pouvait expliquer l'atténuation des ondes sismiques mieux que les modèles traditionnels. Ils n'ont pas tenté de prédire quand le prochain grand séisme se produirait ; ils se sont plutôt concentrés sur la manière de reconstruire avec précision les secousses qui surviennent après qu'un événement a été enregistré.

L'équipe a commencé par cartographier la scène géologique. Ils ont examiné les segments spécifiques de la faille qui passent sous la mer, notamment les branches de Tekirdağ, du bassin central et de Kumburgaz, et ont étudié l'historique des séismes dans la région, tels que le séisme majeur de l'İzmit en 1999 et l'événement plus petit de Silivri en 2019. Ce contexte géologique était essentiel car les mathématiques qu'ils développaient devaient être ancrées dans la forme et la structure réelles de la zone de faille. Ils ont ensuite introduit une nouvelle façon de décrire le comportement de la roche. Imaginez la croûte terrestre non pas comme un simple ressort qui reprend sa forme instantanément, mais comme un matériau qui se relâche lentement au fil du temps, conservant la mémoire de la façon dont il a été étiré ou comprimé. Les chercheurs ont utilisé un outil mathématique appelé dérivée de Caputo pour décrire cet effet de mémoire. Cet outil leur a permis de créer un modèle qui relie les propriétés physiques de la roche à la façon dont les ondes perdent de l'énergie lors de leur déplacement.

En utilisant ce nouveau modèle, les chercheurs ont dérivé le comportement de trois types différents d'ondes sismiques : les ondes P, qui sont les plus rapides et arrivent en premier ; les ondes S, qui font vibrer le sol latéralement ; et les ondes de Rayleigh, qui roulent à la surface. Ils ont prouvé mathématiquement que leur modèle entraîne toujours une perte d'énergie, ce qui est une condition nécessaire pour toute description réaliste d'un séisme. Si un modèle suggérait que les ondes pouvaient gagner de l'énergie en voyageant, cela serait physiquement impossible. Leurs équations ont montré que l'approche fractionnaire produit une courbe de perte d'énergie lisse et continue qui correspond à la façon dont les roches réelles se comportent sur une large gamme de fréquences. Il s'agit d'une amélioration significative par rapport aux anciens modèles, qui peinent souvent à s'ajuster aux données à travers différentes fréquences sans devenir excessivement compliqués.

Pour tester si cette théorie fonctionne réellement, l'équipe a mené une expérience contrôlée utilisant des données générées par ordinateur. Ils ont créé un scénario de séisme synthétique où ils connaissaient exactement les règles de l'atténuation des ondes. Ils ont ensuite demandé à leur modèle fractionnaire de découvrir ces règles en se basant uniquement sur les motifs des ondes. Les résultats ont été frappants. Le modèle fractionnaire a réussi à récupérer les paramètres corrects avec une très grande précision, réduisant l'erreur de ses prédictions de près de 87 % par rapport à un modèle traditionnel utilisant des règles à nombres entiers. Dans le modèle traditionnel, les chercheurs avaient dû forcer les chiffres pour qu'ils correspondent, ce qui entraînait une mauvaise adéquation avec le comportement réel des ondes. Le modèle fractionnaire, en permettant cet effet de mémoire « intermédiaire », a capturé la physique de la perte d'énergie de manière beaucoup plus naturelle.

L'étude a également examiné comment ce cadre pourrait être appliqué aux données réelles de la région de Marmara. Les chercheurs ont identifié le séisme de Silivri de 2019 comme un candidat idéal pour un premier test en conditions réelles. Ils ont exposé un plan clair sur la manière dont les ingénieurs et les sismologues pourraient télécharger les enregistrements réels de cet événement et utiliser leurs nouvelles équations pour reconstruire l'intensité des secousses. Cela permettrait de calculer des mesures précises du mouvement du sol, telles que l'accélération ou la vitesse maximale, qui sont les chiffres utilisés par les architectes pour concevoir des bâtiments résistants aux séismes. Les auteurs ont pris soin de préciser que cet outil sert à analyser ce qui s'est déjà produit, et non à prédire l'avenir. Il ne nous dit pas quand la prochaine grande rupture se produira, mais il offre un prisme bien plus net pour comprendre la violence des secousses une fois qu'elles surviennent.

En combinant une preuve mathématique rigoureuse avec une carte géologique claire et un plan d'analyse de données réelles, cet article fait passer la conversation de la théorie abstraite à l'application pratique. Il démontre que l'approche fractionnaire n'est pas seulement une curiosité mathématique, mais une méthode viable et plus précise pour décrire la réalité complexe des ondes sismiques en mer de Marmara. Ce travail comble le fossé entre la géologie désordonnée et stratifiée de la zone de faille et les équations propres et prévisibles nécessaires à la sécurité de l'ingénierie. Bien que le calibrage complet par rapport à des milliers d'enregistrements sismiques réels soit une tâche future, cette étude a posé les bases nécessaires, prouvant que la méthode fonctionne en théorie et en simulation, et qu'elle est prête à être testée contre le sol réel qui se trouve sous Istanbul.