Deformation Due to a Pair of Infinite and a Pair of infinite-finite planar faults in Fractional Standard Linear Solid

Résumé Technique : Déformation due à une paire de failles planes infinies et une paire de failles infinies-finies dans un solide linéaire standard fractionnaire

Énoncé du problème
Les régions sismiquement actives contiennent fréquemment des failles multiples en interaction où le glissement le long d'une faille influence considérablement l'accumulation et la distribution des contraintes dans les failles voisines. Bien que les études précédentes aient modélisé les interactions entre failles en utilisant des dérivées d'ordre entier au sein de milieux viscoélastiques (par exemple, les modèles de Maxwell ou de solide linéaire standard), ces approches traditionnelles négligent souvent les effets de « mémoire » inhérents aux matériaux géologiques. Cet article traite de la déformation du système lithosphère-asthénosphère lors de la période aseismique en modélisant deux scénarios spécifiques impliquant des failles de décrochement par glissement rampant (creep) en interaction :

1. Cas-I : L'interaction entre deux failles enfouies infinies.
2. Cas-II : L'interaction entre une faille enfouie infinie et une faille de surface finie.

L'étude postule que le milieu terrestre se comporte comme un demi-espace viscoélastique caractérisé par un modèle de solide linéaire standard fractionnaire (FSLS), qui capture mieux le comportement rhéologique dépendant du temps et les effets de mémoire des roches crustales par rapport aux modèles d'ordre entier.

Méthodologie
Les auteurs ont développé des modèles mathématiques théoriques pour dériver des solutions analytiques pour les composantes du déplacement, de la contrainte et de la déformation. La méthodologie comprend les étapes clés suivantes :

• Formulation mathématique : Le milieu est défini par une équation constitutive de solide linéaire standard fractionnaire impliquant l'opérateur de dérivée fractionnaire de Caputo (${}_{0}D^{\alpha}_{t}$) où l'ordre $\alpha \in (0, 1]$. Les relations contrainte-déformation intègrent la rigidité effective ($\mu$) et la viscosité effective ($\eta$).
• Équations de gouvernance : L'analyse suppose une déformation quasi-statique (négligeant les forces inertielles). Les équations de mouvement de la contrainte et les conditions aux limites sont établies pour les deux géométries de failles (infinie et finie).
• Technique de résolution : Les auteurs utilisent la transformation de Laplace pour gérer les équations différentielles fractionnaires dépendant du temps. Les solutions sont dérivées en utilisant :
  • Le Principe de correspondance (Rybicki, 1971).
  • Une technique de fonction de Green modifiée (Maruyama, 1966) pour rendre compte des conditions de dislocation à travers les plans de faille.
  • La fonction de Mittag-Leffler ($E_{\alpha}$), qui apparaît naturellement de la transformation de Laplace inverse des opérateurs fractionnaires.
• Scénarios analysés :
  • Les solutions sont dérivées pour l'état avant le fluage de la faille (accumulation de contrainte tectonique) et après le fluage de la faille (dislocation due au glissement).
  • Des vitesses de fluage ($V_1, V_2$) et des fonctions de fluage dépendantes de la profondeur sont appliquées pour simuler le mouvement des failles $F_1$ et $F_2$.
  • Des calculs numériques sont effectués en utilisant des paramètres dérivés de données géologiques (par exemple, rigidité $\mu = 3,5 \times 10^{10}$ N/m², viscosité $\eta = 5 \times 10^{19}$ Pa·s) pour générer des résultats graphiques pour diverses inclinaisons de faille ($\theta$), ordres fractionnaires ($\alpha$) et vitesses de fluage.

Contributions clés

1. Modélisation fractionnaire de l'interaction des failles : L'article étend les modèles existants d'interaction entre failles en incorporant le calcul fractionnaire, permettant une représentation plus précise du comportement viscoélastique dépendant de la mémoire de la lithosphère.
2. Solutions analytiques pour des géométries complexes : L'étude fournit des expressions analytiques sous forme fermée pour les champs de déplacement, de contrainte et de déformation pour deux configurations distinctes : une paire de failles enfouies infinies et un système mixte de failles de surface infinies et finies.
3. Quantification des effets de mémoire : La recherche démontre explicitement comment l'ordre fractionnaire $\alpha$ influence l'accumulation et la libération de contrainte, liant le paramètre mathématique à la « mémoire » physique du milieu géologique.

Résultats et conclusions numériques
L'analyse numérique, appuyée par des représentations graphiques, produit plusieurs observations spécifiques :

• Magnitudes de déplacement : Le changement de déplacement de surface à proximité des failles est de l'ordre de $10^{-3}$ m/an. Cela concorde bien avec les données géophysiques observées (par exemple, les rapports de l'USGS de 2 à 6 mm/an), alors que les modèles précédents d'ordre entier prédisaient des valeurs de l'ordre de $10^{-5}$, montrant une divergence avec les observations réelles.
• Influence de la géométrie des failles :
  • Dans le Cas-I (deux failles infinies), la magnitude du déplacement est plus élevée près de la seconde faille qu'à proximité de la première.
  • Dans le Cas-II (faille infinie et faille finie), la magnitude du déplacement est plus faible près de la faille finie qu'à proximité de la faille infinie.
  • Les pics de déplacement se situent près de la trace de la faille et diminuent avec la distance, tendant finalement vers zéro.
• Effet de l'ordre fractionnaire ($\alpha$) : L'étude constate que le taux d'accumulation et de libération de contrainte dépend fortement de $\alpha$. Des valeurs plus faibles de $\alpha$ (indiquant des effets de mémoire plus forts) entraînent des magnitudes plus élevées d'accumulation de contrainte de cisaillement et une relaxation plus lente. À mesure que $\alpha$ augmente vers 1 (approchant le comportement entier), le taux de libération de contrainte augmente.
• Force motrice et distribution des contraintes : La force motrice totale à travers les failles montre une libération rapide après le début du fluage. La magnitude de cette libération varie selon l'inclinaison de la faille et la vitesse de fluage. Les cartes de contours révèlent que l'accumulation de contrainte est la plus significative près de la surface libre et des zones de dislocation de la faille, tandis que la réduction de contrainte est proéminente aux profondeurs correspondant aux emplacements des failles.
• Analyse de la déformation : Les variations de la déformation de cisaillement ($e_{12}$) se situent dans la plage de $10^{-4}$, ce qui est cohérent avec les données d'événements réels tels que le séisme de Tohoku-oki de 2011.

Signification et affirmations
Les auteurs affirment que leur modèle fractionnaire offre une représentation plus générale et plus précise de la dynamique d'interaction des failles par rapport aux modèles traditionnels d'ordre entier. Plus précisément :

• L'approche fractionnaire capture avec succès les effets de mémoire et la non-linéarité des milieux viscoélastiques, qui sont souvent négligés dans les modèles standards.
• Les prédictions du modèle concernant les taux de déplacement ($10^{-3}$) montrent une forte concordance avec les données géophysiques observées, suggérant que le calcul fractionnaire est un outil supérieur pour prédire le comportement sismique et la déformation de subsurface.
• Les résultats indiquent que l'interaction entre les failles influence le champ de contrainte de manière plus étendue qu'une faille unique, et que l'ordre fractionnaire $\alpha$ est un paramètre critique déterminant le taux de libération de contrainte.

L'article conclut que ces modèles théoriques améliorent la compréhension de la déformation de subsurface et de son rôle dans l'influence de la dynamique des failles et de la génération de séismes, fournissant un cadre robuste pour l'analyse des périodes aseismiques dans les zones sismiquement actives.