Elastic wave propagation in fractured media with spring-type and frictional contact deformation laws

Le sol sous nos pieds est rarement un bloc solide et ininterrompu. Dans de nombreux endroits, des réservoirs de pétrole profonds à la croûte terrestre elle-même, la roche est sillonnée de fractures — de fines fissures planes qui peuvent s'étendre sur des kilomètres. Lorsque l'énergie voyage à travers ce paysage fracturé, qu'il s'agisse d'une onde sismique provenant d'un tremblement de terre ou d'une impulsion sonore utilisée pour inspecter un matériau, les fissures changent l'histoire. Elles dispersent l'énergie, affaiblissent le signal et provoquent parfois même le frottement ou le glissement des parois rocheuses les unes contre les autres. Comprendre exactement comment ces ondes interagissent avec les fissures est vital pour les géologues qui tentent de cartographier ce qui se trouve sous terre et pour les ingénieurs qui assurent la sécurité des structures. Cependant, la physique de cette interaction est notoirement difficile à définir. Une simple fissure n'est pas seulement un vide ; c'est une limite complexe où les parois rocheuses peuvent presser l'une contre l'autre, glisser avec friction, ou s'ouvrir et se refermer brusquement, tout en étant traversées par l'onde en une fraction de seconde.

Pendant des décennies, les scientifiques se sont appuyés sur des modèles simplifiés pour prédire ces comportements, traitant souvent les fissures comme si elles étaient composées de ressorts invisibles qui s'étirent et se compriment. Bien qu'utiles, ces modèles s'effondrent lorsque les ondes sont assez fortes pour forcer les parois rocheuses à se toucher, à glisser ou à se séparer complètement. Ils ne peuvent pas capturer le mouvement de type « stick-slip » (adhérence-glissement) où la friction maintient la roche en place jusqu'à ce que la force devienne trop grande, ni décrire avec précision le durcissement non linéaire qui se produit lorsqu'une fissure est pressée pour se fermer. Une nouvelle étude menée par des chercheurs de l'Université de Bergen, du Centre de recherche norvégien et de l'Université technique de Munich aborde ces lacunes en construisant un cadre de calcul complet qui simule le mouvement des ondes élastiques à travers la roche fracturée avec un réalisme sans précédent. Au lieu de s'appuyer sur des hypothèses simplifiées, l'équipe a développé une méthode capable de gérer quatre niveaux de complexité distincts, allant d'un comportement de base de type ressort à la mécanique de contact complète où la friction et l'ouverture et fermeture physiques des fissures sont calculées en temps réel.

Les chercheurs ont construit un environnement numérique unifié où la roche et les fissures sont traitées comme un système unique et interconnecté. Dans leurs simulations, la roche est représentée comme un volume tridimensionnel, tandis que les fractures sont modélisées comme des surfaces bidimensionnelles intégrées en son sein. Cela permet à l'ordinateur de suivre comment une onde voyageant à travers la roche frappe une fissure, comment la fissure se déforme, et comment cette déformation renvoie des ondulations dans la roche. L'équipe a testé quatre manières spécifiques dont les fissures pourraient se comporter. Les deux premiers modèles utilisaient des lois de type ressort, où la force à travers la fissure est proportionnelle au mouvement de la roche, avec une version permettant une relation linéaire simple et une autre utilisant une règle non linéaire plus complexe qui tient compte du fait que les fissures deviennent plus difficiles à comprimer à mesure qu'elles se ferment. Les deux autres modèles introduisaient la réalité de la friction et du contact. Dans ces scénarios, les faces de la fissure sont autorisées à se toucher et à glisser l'une contre l'autre, régies par des lois de friction similaires à celles qui empêchent les pneus d'une voiture de déraper, tout en imposant une règle stipulant que les parois rocheuses ne peuvent pas se traverser mutuellement.

Pour s'assurer de l'exactitude de leur nouveau cadre, l'équipe a réalisé une série de tests rigoureux. Ils ont simulé des ondes frappant une fissure unique et ont comparé leurs résultats informatiques à des solutions mathématiques connues pour des cas simples. Ils ont constaté que leur méthode pouvait prédire le comportement de l'onde avec une grande précision, correspondant aux attentes théoriques tant pour la vitesse de l'onde que pour la quantité d'énergie qui passait à travers ou rebondissait. Ils ont ensuite poussé le système plus loin, simulant des ondes suffisamment fortes pour provoquer l'ouverture, l'adhérence et le glissement des fissures. Dans ces scénarios plus complexes, où les faces de la fissure frottaient l'une contre l'autre avec un coefficient de friction de 1,0, les simulations ont montré que la méthode restait stable et produisait des résultats cohérents. Les chercheurs ont vérifié que leur approche pouvait gérer les changements soudains de physique qui se produisent lorsqu'une fissure passe de l'adhérence au glissement, une transition qui pose souvent problème aux modèles plus simples.

Lorsque l'équipe a comparé les quatre modèles différents côte à côte, les différences dans les motifs d'ondes résultants étaient frappantes. Les simulations ont révélé que le choix du modèle importe considérablement. Lorsqu'une fissure était autorisée à se fermer et à se durcir selon la loi non linéaire, elle transmettait beaucoup plus d'énergie d'onde qu'une fissure modélisée avec un simple ressort linéaire. C'est parce que le modèle non linéaire tient compte du fait qu'à mesure que la fissure se ferme, les parois rocheuses pressent l'une contre l'autre plus étroitement, rendant la fissure effectivement plus rigide et plus transparente pour l'onde. De plus, l'inclusion de la friction a modifié la forme des fronts d'onde. Dans les modèles frictionnels, l'onde qui émerge de l'autre côté de la fissure était plus large et plus déformée que dans les versions sans friction. Les chercheurs ont observé que la friction faisait que les faces de la fissure restaient collées plus longtemps, retardant la libération de l'énergie et altérant la façon dont l'onde se dispersait.

L'étude ne s'est pas arrêtée aux fissures uniques. Les chercheurs ont appliqué leur cadre à un domaine tridimensionnel contenant six fractures intersectées, créant un petit réseau de fissures aux propriétés différentes. Dans ce cadre complexe, ils ont assigné différentes limites de fermeture maximale à différents ensembles de fractures, ce qui signifie que certaines fissures pouvaient se fermer plus étroitement que d'autres. La simulation a montré que le champ d'ondes devenait hautement hétérogène. Les fractures qui étaient plus rigides et pouvaient se fermer plus étroitement réfléchissaient plus d'énergie et en transmettaient moins, créant des ombres distinctes dans le motif de l'onde. Les chercheurs ont été capables de suivre exactement à quel point chaque fissure s'était fermée et à quel point elle avait glissé, révélant que même de petites différences dans les propriétés physiques des fissures menaient à des réponses de déformation sensiblement différentes. Par exemple, une fissure autorisée à se fermer plus étroitement présentait un comportement d'« adhérence » plus marqué, maintenant ses faces ensemble plus longtemps avant de glisser, ce qui affectait à son tour l'onde qui la traversait.

Les implications de ces découvertes s'étendent au-delà des résultats immédiats. En démontrant qu'un seul cadre de calcul peut gérer tout, du simple comportement de ressort à la mécanique de contact frictionnelle complexe, les chercheurs ont fourni un outil qui peut être utilisé pour interpréter les données sismiques réelles avec une plus grande fidélité. La capacité de simuler la façon dont les ondes se dispersent et s'atténuent en présence de fissures frictionnelles, d'ouverture et de fermeture signifie que les géoscientifiques peuvent mieux distinguer les différents types de formations rocheuses et de fractures remplies de fluides. L'étude confirme que ignorer les aspects non linéaires et frictionnels de la mécanique des fractures peut conduire à des erreurs significatives dans la prédiction de la manière dont l'énergie circule dans le sous-sol. Bien que le travail actuel se concentre sur les interactions purement mécaniques, le cadre est construit sur une base compatible avec les modèles de flux de fluides, suggérant une voie vers des simulations encore plus réalistes où l'eau ou le pétrole circulant à travers les fissures interagit avec l'onde passante. Ce travail démontre qu'en embrassant toute la complexité du comportement des fissures rocheuses, nous pouvons construire une image plus claire du monde caché sous nos pieds.