CONTEXTE DE L'INVENTION 100011 Une faille peut être considérée comme étant une discontinuité de surface tridimensionnelle complexe finie dans un volume de sol ou de roche. Des fractures, y compris, sans limitation, des diaclases, des filons, des dil(es, des couches de solution sous pression avec des stylolites, etc., peuvent être amenées à se propager de manière intentionnelle, en vue d'augmenter la perméabilité dans des formations telles que les schistes, dans lesquelles l'optimisation du nombre, de l'emplacement, et de la taille des fractures dans la formation augmente le rendement de ressources telles que le gaz de schistes. 10 100021 La contrainte, en mécanique des milieux continus, peut être considérée comme une mesure des forces internes agissant au sein d'un volume. Cette contrainte peut être définie comme une mesure de la force moyenne par unité de surface au niveau d'une surface au sein du volume sur lequel agissent des forces internes. Des forces internes peuvent être produites entre les particules contenues dans le volume en réaction à des forces externes 15 appliquées au volume. 100031 Une compréhension de l'origine et de l'évolution des failles et de l'histoire tectonique des régions faillées peut être établie en reliant des données d'orientation des failles, de direction de glissement, géologiques et géodésiques à l'état de la contrainte dans la croûte terrestre. Dans certains problèmes inverses, les directions des contraintes principales 20 distantes et un rapport de leurs grandeurs sont contraints par une analyse de données de terrain concernant les orientations des failles et les directions de glissement telles qu'elles sont déduites d'artéfacts tels que les striures apparaissant sur des surfaces de failles exposées. Si les failles supportent le frottement de glissement, les résultats de l'inversion de la contrainte pourraient être différents par comparaison à des failles sans frottement. 25 RÉSUMÉ DE L'INVENTION 100041 De manière générale, selon un aspect, l'invention concerne un procédé de prédiction de la contrainte régionale d'un volume souterrain. Le procédé comprend les opérations de : obtenir un modèle du volume souterrain, dans lequel le modèle comprend une matrice de modèle représentant une relation entre un résultat de glissement de faille modélisé généré par le modèle et une condition aux limites appliquée au modèle, dans lequel la condition aux limites comprend un attribut de contrainte régionale et un attribut de frottement de faille ; calculer, au moyen du modèle, le résultat de glissement de faille modélisé sur la base d'une valeur sélectionnée de l'attribut de contrainte et d'une valeur, sélectionnée de l'attribut de frottement ; calculer une fonction de coût représentant une différence entre le résultat de glissement de faille modélisé et une mesure du volume souterrain ; et minimiser la fonction de coût par ajustement itératif d'au moins la valeur sélectionnée de l'attribut de contrainte, dans lequel l'ajustement itératif de la valeur sélectionnée en vue de minimiser la fonction de coût génère une prédiction de la contrainte régionale du volume souterrain. [0005] De manière générale, selon un aspect, l'invention concerne un système de prédiction de l'activité de faille d'un volume souterrain. Le système comprenant : un dispositif de détection configuré pour obtenir une mesure du volume souterrain ; un moteur de modélisation de contrainte et de fracture configuré pour : obtenir un modèle du volume souterrain, dans lequel le modèle comprend une matrice de modèle représentant une relation entre un résultat de glissement de faille modélisé généré par le modèle et une condition aux limites appliquée au modèle, dans lequel la condition aux limites comprend un attribut de contrainte régionale et un attribut de frottement de faille ; calculer, en utilisant le modèle, le résultat de glissement de faille modélisé sur la base d'une valeur sélectionnée de l'attribut de contrainte et d'une valeur sélectionnée de l'attribut de frottement ; calculer une fonction de coût représentant une différence entre le résultat de glissement de faille modélisé et une mesure du volume souterrain ; et minimiser la fonction de coût par ajustement itératif d'au moins la valeur sélectionnée de l'attribut de contrainte et de l'attribut de frottement, dans lequel l'ajustement itératif de la valeur sélectionnée en vue de minimiser la fonction de coût génère une prédiction de la contrainte régionale du volume souterrain ; et un dispositif de commande configuré pour générer, sur la base de la prédiction de la contrainte régionale, un signal de commande d'une opération sur le terrain du volume souterrain. [0006] De manière générale, selon un aspect, l'invention concerne un support lisible par ordinateur non volatil stockant des instructions pour prédire l'activité de faille d'un volume souterrain. Les instructions, lorsqu'elles sont exécutées par un processeur d'ordinateur, comprennent des fonctionnalités pour : obtenir un modèle du volume souterrain, dans lequel le modèle comprend une matrice de modèle qui représente une relation entre un résultat de glissement de faille modélisé généré par le modèle et une condition aux limites appliquée au modèle, dans lequel la condition aux limites comprend un attribut de contrainte régionale et un attribut de frottement de faille ; calculer, en utilisant le modèle, le résultat de glissement de faille modélisé sur la base d'une valeur sélectionnée de l'attribut de contrainte et d'une valeur sélectionnée de l'attribut de frottement ; calculer une fonction de coût représentant une différence entre le résultat de glissement de faille modélisé et une 10 mesure du volume souterrain ; et minimiser la fonction de coût par ajustement itératif d'au moins la valeur sélectionnée de l'attribut de contrainte et de l'attribut de frottement, dans lequel l'ajustement itératif de la valeur sélectionnée pour minimiser la fonction de coût génère une prédiction de la contrainte régionale du volume souterrain. [0007] Le présent résumé est présenté en vue d'introduire un choix de concepts qui sont 15 décrits par ailleurs ci-après dans la description détaillée. Le présent résumé n'est pas destiné à identifier des caractéristiques déterminantes ou essentielles de l'objet revendiqué, ni à être utilisé en tant qu'aide pour délimiter l'étendue de l'objet revendiqué. BRÈVE DESCRIPTION DES DESSINS [0008] Des modes de réalisation de l'inversion de la contrainte régionale en utilisant des 20 failles de frottement sont décrits en référence aux figures suivantes. Des références numériques identiques sont utilisées dans l'ensemble des figures pour désigner des caractéristiques et composants identiques. [0009] La figure 1 est un schéma d'un exemple de système de modélisation de la contrainte et des fractures. 25 [0010] La figure 2 est un schéma de principe d'un exemple de l'environnement informatique destiné à effectuer la modélisation de contrainte et de fracture en utilisant des failles de frottement. [0011] La figure 3 est un schéma de principe d'un exemple de moteur de modélisation de contrainte et de fracture. [0012] La figure 4 est un schéma de principe comparant des techniques de restitution de la paléocontrainte. [0013] La figure 5 est un schéma d'un exemple de procédé appliqué en vue de fracturer et de conjuguer des plans de failles utilisant des ensembles de données sans information de grandeur. [0014] La figure 6 est un organigramme d'un exemple de procédé de modélisation de contrainte et de fracture utilisant des failles de frottement. [0015] La figure 7 et la figure 8 représentent un exemple de modélisation de contrainte, de fracture et d'activité de faille utilisant des failles de frottement. DESCRIPTION DÉTAILLÉE [0016] Dans la description détaillée présentée ci-après de modes de réalisation, de nombreux détails spécifiques sont présentés afin de procurer une compréhension plus approfondie. Cependant, l'homme du métier pourra noter que des modes de réalisation peuvent être mis en oeuvre sans ces détails spécifiques. Dans d'autres cas, des caractéristiques bien connues n'ont pas été décrites en détail afin d'éviter de compliquer inutilement la description. [0017] Des exemples de modes de réalisation fournissent un procédé et un système de modélisation de contrainte et de fracture utilisant des failles de frottement. Un exemple de système simule un modèle de contrainte de champ lointain pour un volume souterrain, en calculant une ou plusieurs contrainte, déformation et/ou déplacement. Dans un ou plusieurs modes de réalisation, la relation reliant les conditions aux limites et la géométrie (par exemple la géométrie des failles) du modèle de contrainte de champ lointain à la distribution des glissements sur les failles est précakulée sous la forme d'une matrice de modèle représentée dans l'équation (1) ci-après. [0018] S = M*be (1) [0019] où S est la distribution des glissements sur les failles, M est la matrice de modèle, * représente la multiplication matricielle/vectorielle, et be est la condition aux limites appliquée au modèle et comprend des paramètres (désignés comme étant l'attribut de contrainte) liés à la contrainte de champ lointain et des paramètres (désignés comme étant l'attribut de frottement) liés au frottement sur les surfaces de faille. Dans un ou plusieurs modes de réalisation, la matrice de modèle comprend une dépendance non linéaire vis-à-vis de l'attribut de frottement. En conséquence, la matrice de modèle précalculée dépend de caractéristiques de frottement des failles dans le volume souterrain. Sur la base de la matrice de modèle précalculée, le système peut générer des résultats relativement rapides 10 compte tenu d'une simulation de modèle de contrainte de champ lointain en ajustant le rapport/l'orientation de champ lointain en vue de restituer des paramètres de paléocontrainte, tels que des paramètres de contrainte, de déformation et de déplacement pour un point quelconque (désigné comme étant des points d'observation, ou des points de données) dans le volume souterrain lorsque l'utilisateur fait varier la valeur de contrainte de 15 champ lointain. Le système peut restituer un ou plusieurs événements tectoniques, ou un tenseur de contrainte représenté par un rapport de grandeurs principales et une orientation associée, en utilisant au moins l'une de la géométrie de faille, du frottement de glissement de faille, des données de puits de forage (y compris des donnés d'orientation de fracture et de plan de faille secondaire), de données de système de positionnement global (GPS), de 20 données de radar à ouverture synthétique interférométrique (InSAR), des horizons plissés et faillés, des inclinomètres, des lignes et des stries de glissement de faille. Le système peut utiliser différents types de données géologiques provenant de l'interprétation sismique, de lectures de puits de forage et d'observations sur le terrain en vue de fournir de nombreux résultats, comme la propagation prédite de fractures sur la base d'un champ de contrainte 25 perturbé. Dans l'ensemble du présent document, le terme de "contrainte de champ lointain" peut également désigner une contrainte régionale ou une contrainte tectonique. Définitions [0020] Dans la description ci-après, certaines variables sont utilisées afin de simplifier la présentation. Le tableau 1 ci-après indique chaque variable pouvant être utilisée et la définition correspondante des variables conformément à l'un ou plusieurs des modes de réalisation. Tableau 1 Variable Définition GR Tenseur de contrainte de champ lointain régionale GH Contrainte principale horizontale maximum (extraite de CYR) OEh Contrainte principale horizontale minimum (extraite de o-R) av Contrainte principale verticale (extraite de aR) R Rapport de contrainte défini par (o-2-a3)/ (ai- a3) et e [0,1] R' Autre rapport de contrainte e [0,3], défini par R' - R pour un régime de faille normale, R'=2 - R pour un régime de faille décrochante, et R'=2 + R pour un régime de faille inverse Matrice de rotation définie par les trois angles d'Euler n Nombre d'éléments triangulaires constituant les surfaces des failles P Nombre de points d'observation ou points de données b Vecteur de glissement sur un élément triangulaire d'une surface de faille 0 Orientation de la contrainte de champ lointain régionale. Dans l'un ou plusieurs des modes de réalisation, l'orientation est définie vers le Nord et en sens horaire ai Valeur de la contrainte principale maximum d'un tenseur de contrainte (72 Valeur de la contrainte principale intermédiaire d'un tenseur de contrainte O-3 Valeur de la contrainte principale minimum d'un tenseur de contrainte R o Matrice de rotation suivant l'axe z d'un angle 0 (voir plus haut) k Paramètre d'échelle défini sous la forme 0-2/ ai Composante (id) d'un tenseur de contrainte ai Paramètres d'optimisation utilisés pour le principe de superposition Ep Tenseur de déformation calculé au point P OEp Tenseur de contrainte calculé au point P 21p Champ de déplacement calculé au point P be Vecteur de glissement calculé sur l'élément triangulaire e A Matrice 3x3 reliant a et OER np Normale à un plan au point P Champ de déplacement calculé au point P, où l'exposant c indique "calculé" Variable Définition Champ de déplacement calculé au point P, où l'exposant m indique "mesuré" /-1 Coefficient de frottement du glissement sur la surface de la faille Exposé général [0021] Le présent document décrit la modélisation de contrainte et de fracture en utilisant des failles de frottement. Connaissant diverses données d'entrée, comme la géométrie des failles avec un coefficient de frottement du glissement, et des ensembles de données sélectionnables ou facultatifs ou des mesures de données, parmi lesquelles un ou plusieurs d'un rejet de faille, de direction de rejet en profondeur ou de stries de glissement, de mesures de contrainte, de données de fractures, d'orientations de plan de faille secondaire, de données GPS, de données InSAR, de données géodésiques provenant d'inclinomètres de surface, de télémétrie laser, etc., l'exemple de système permet de générer ou de reconstituer rapidement de nombreux types de résultats. Les systèmes et procédés décrits ici appliquent de multiples passes de simulation géomécanique à des surfaces de failles présentant une géométrie complexe dans au moins trois dimensions (3D), et les failles sont, par nature, de dimension finie et non infinie ou semi-infinie. Les résultats peuvent comprendre, à titre d'exemple, un ou plusieurs des paramètres de contrainte, de déformation, et/ou de déplacement en réponse à un ou plusieurs d'une requête d'un utilisateur ou d'un paramètre mis à jour, des états de contrainte distants pour de multiples événements tectoniques, la prédiction de la fracturation future voulue, la différenciation de fractures préexistantes par rapport à des fractures induites, etc. Les diverses données d'entrée peuvent être déduites de données de puits de forage, d'interprétation sismique, d'observation sur le terrain, etc. [0022] Les exemples de systèmes et de procédés décrits ci-après peuvent être appliqués à de nombreuses opérations différentes de réservoir et souterraines, y compris, sans limitation, les opérations d'exploration et de production de gaz naturel et autres hydrocarbures, le stockage du gaz naturel, la fracturation hydraulique et la stimulation matricielle en vue d'augmenter la production des réservoirs, la gestion des ressources en eau, y compris le développement et la protection environnementale de nappes phréatiques et autres ressources en eau, le captage et le stockage souterrain du dioxyde de carbone (CO2), etc. [0023] Dans un exemple de réalisation, un système applique une technique des éléments limites 3D en utilisant des failles de frottement pour des milieux isotropes hétérogènes couvrant un espace entier (c'est-à-dire en considérant l'effet de la surface de la terre) ou un demi-espace (ne tenant pas compte de l'effet de la surface de la terre). Sur la base de la matrice de modèle précalculée, l'exemple de système permet d'évaluer une fonction de coût en vue de générer des résultats rapides, comme des paramètres de contrainte, de déformation et de déplacement pour un point quelconque dans un volume souterrain lorsque l'utilisateur fait varier la valeur de contrainte de- champ lointain. Dans une forme de réalisation, le système peut utiliser la géométrie des failles et le coefficient de frottement du glissement de faille et des données de forage, comprenant par exemple l'orientation des fractures, des données de plan de faille secondaire et/ou une mesure de contrainte in situ par fracture hydraulique, pour restituer un ou plusieurs évènements tectoniques ou un tenseur de contrainte représenté par un rapport de grandeur principal et l'orientation associée. Le système peut utiliser de nombreux types différents de données géologiques provenant de l'interprétation sismique, de mesures de forage, et d'observations sur le terrain en vue de fournir divers résultats, tels qu'une propagation prédite des fractures sur la base du champ de contrainte perturbé. [0024] La figure 1 représente un exemple de système de modélisation de contrainte et de fracture 100. L'exemple de système 100 permet de résoudre une variété de problèmes géomécatiques. La géométrie des failles peut être connue (et facultativement, le rejet de faille mesuré et les contraintes d'inégalité imposées, par exemple normale, chevauchante, etc., peuvent être cumules). L'utilisateur peut avoir accès à une ou plusieurs de données provenant de puits de forages (par exemple, orientation de fracture, mesures de contrainte in-situ, plans de faille secondaire), données géodésiques (par exemple, InSAR, GPS, et inclinomètre), et/ou horizons interprétés. Un exemple de moteur de modélisation de contrainte et de fracture 102 et/ou des exemples de procédés correspondants peuvent reconstituer le régime d'état de contrainte et tectonique distant pour un/des événement(s) tectonique(s) pertinent(s), ainsi que la discontinuité de déplacement sur les failles, et estimer le déplacement et les champs de déformation et de contrainte perturbés en tout point au sein du système, par exemple, pour un coefficient de glissement imposé donné sur des failles (c'est-à-dire le coefficient de frottement d'une faille glissante, également désigné coefficient de faille de glissement). [0025] En utilisant la matrice de modèle précalculée, le système 100 ou moteur 102 de modélisation de contrainte et de fracture peut effectuer une simulation en modifiant les conditions aux limites (par exemple la contrainte de champ lointain ou le coefficient de frottement du glissement) sans recalculer le modèle de matrice. Ensuite, comme cela a été introduit plus haut, des applications de l'exemple de système 100 peuvent comprendre un 10 ou plusieurs d'une interpolation de contrainte et d'une modélisation de fracture, d'une restitution d'un ou plusieurs événement(s) tectonique(s), d'un contrôle de qualité concernant des failles interprétées, d'un calcul en temps réel de champs de contrainte et de déplacement perturbés lorsque l'utilisateur effectue un ou plusieurs d'une estimation de paramètres, d'une prédiction de la propagation de fracture, d'une distinction de fractures 15 préexistantes par rapport à des fractures induites, et de nombreuses autres applications. Exemple d'environnement [0026] La figure 2 représente l'exemple de système 100 de la figure 1 dans le contexte d'un environnement informatique, dans lequel une modélisation de contrainte et de fracture en utilisant des failles de frottement peut être effectuée. 20 [0027] Dans l'exemple illustré, le dispositif informatique 200 est relié par voie de communication via des dispositifs de détection (et des dispositifs de commande) à un milieu réel, à titre d'exemple un volume de terre souterrain réel 202, un réservoir 204, un bassin de dépôt, un fond marin, etc., et des puits associés 206 destinés à produire une ressource de pétrole, à la gestion de ressources en eau, ou à des services se rapportant au 25 carbone, etc. [0028] Dans la forme de réalisation représentée, un dispositif informatique 200 met en oeuvre un composant, tel que le moteur de modélisation de contrainte et de fracture 102 et le moteur d'affichage graphique 231. Le moteur de modélisation de contrainte et de fracture 102 et le moteur d'affichage graphique 231 sont illustrés sous la forme de logiciels, mais peuvent être mis en oeuvre sous la forme de matériels ou sous la forme d'une combinaison de matériels et d'instructions logicielles. 100291 Le moteur de modélisation de contrainte et de fracture 102 comprend des fonctionnalités permettant d'effectuer une analyse, en utilisant par exemple une fonction de coût d'une valeur de rapport de contrainte et d'une valeur d'orientation. L'exécution de l'analyse est présentée ci-après. Dans un ou plusieurs modes de réalisation, la valeur du rapport de contrainte peut être définie comme étant (I) dans les équations ci-après. Le moteur de modélisation de contrainte et de fracture 102 peut en outre comprendre des fonctionnalités permettant de calculer, pour un point particulier dans le diagramme du domaine des contraintes, une prédiction de fracture, un champ de contrainte perturbé, et/ou un champ de déplacement. [0030] Lors de son exécution, par exemple sur le processeur 208, le moteur de modélisation de contrainte et de fracture 102 est fonctionnellement connecté à un moteur d'affichage graphique 231. A titre d'exemple, le moteur de modélisation de contrainte et de fracture 102 peut faire partie de la même application logicielle que le moteur d'affichage graphique 231, le moteur d'affichage graphique 231 peut être un module d'extension pour moteur de modélisation de contrainte et de fracture 102, ou un autre procédé peut être utilisé pour connecter le moteur d'affichage graphique 231 au moteur de modélisation de contrainte et de fracture 102. Se référant toujours à la figure 2, le moteur d'affichage graphique 231 comprend en outre des fonctionnalités permettant de recevoir une sélection d'un point, demander une prédiction de fracture, un champ de contrainte perturbé, et/ou un champ de déplacement pour le point particulier au moteur de modélisation de contrainte et de fracture 102 en réponse à la sélection, et présenter la prédiction de fracture, le champ de contrainte perturbé, et/ou le champ de déplacement. Dans un ou plusieurs modes de réalisation, le moteur d'affichage graphique 231 est fonctionnellement connecté à l'unité de commande d'interface utilisateur 230 et à l'afficheur 232. 100311 Le dispositif informatique 200 peut être un ordinateur, un réseau d'ordinateurs, ou un autre dispositif qui comporte un processeur 208, une mémoire 210, une unité de stockage de données 212, et d'autres matériels associés tels qu'une interface réseau 214 et un lecteur de support 216 pour lire et écrire sur un support de stockage amovible 218. Le support de stockage amovible 218 peut par exemple être un disque compact (CD) ; un disque numérique polyvalent / disque vidéo numérique (DVD) ; un disque à mémoire flash, etc. Le support de stockage amovible 218 contient des instructions qui, lorsqu'elles sont exécutées par le dispositif informatique 200, amènent le dispositif informatique 200 à mettre en oeuvre un ou plusieurs exemples -de procédés décrits ici. Par conséquent, le support de stockage amovible 218 peut comprendre des instructions permettant de mettre en oeuvre et d'exécuter l'exemple de moteur de modélisation de contrainte et de fracture 102 et/ou de moteur d'affichage graphique 231. Au moins certaines parties de l'exemple de 10 moteur de modélisation de contrainte et de fracture 102 peuvent être stockées sous la forme d'instructions sur un exemplaire donné du support de stockage amovible 218, du dispositif amovible, ou dans une unité de stockage de données locale 212, en vue d'être chargées dans la mémoire 210 pour leur exécution par le processeur 208. Plus précisément, des instructions logicielles ou un code de programme lisible par ordinateur destiné à mettre en 15 oeuvre des modes de réalisation peuvent être stockés, de manière temporaire ou permanente, en tout ou partie, sur un support lisible par ordinateur non volatil tel qu'un CD, un DVD, un dispositif de stockage local ou distant, une mémoire locale ou distante, une disquette, ou tout autre dispositif de stockage lisible par ordinateur. [0032] Bien que les exemples illustrés de moteur de modélisation de contrainte et de fracture 20 102 et de moteur d'affichage graphique 231 soient représentés sous la forme d'un programme résidant dans la mémoire 210, le moteur de modélisation de contrainte et de fracture 102 et/ou le moteur d'affichage graphique 231 peuvent être mis en oeuvre sous la forme de matériels, tels qu'un circuit intégré spécifique d'applications (ASIC) ou sous la forme d'une combinaison de matériels et de logiciels. 25 [0033] Dans cet exemple de système, le dispositif informatique 200 reçoit des données entrantes 220, comme la géométrie des failles, le frottement de glissement des failles, et de nombreux autres types de données, en provenance de multiples sources, telles que des mesures de puits de forage 222, une observation sur le terrain 224, et une interprétation sismique 226. Le dispositif informatique 200 peut recevoir un ou plusieurs types 30 d'ensembles de données 220 via l'interface réseau 214 qui peut également recevoir des données en provenance d'un réseau (par exemple, l'Intemet 228), telles que des données GPS et des données InSAR. [0034] Le dispositif informatique 200 peut déterminer (ou calculer) et compiler des résultats de modélisation, des résultats de simulateur, et des résultats de commande. L'unité de commande d'affichage 230 peut fournir en sortie des images de modèles géologiques et des images- et des données de simulation à un afficheur 232. Les images peuvent être une simulation 2D ou 3D 234 de résultats de contrainte et de fracture utilisant des failles de frottement. L'exemple de moteur de modélisation de contrainte et de fracture 102 peut également générer une ou plusieurs interfaces utilisateur visuelles (UIs) pour fournir en entrée et/ou afficher des données. [0035] L'exemple de moteur de modélisation de contrainte et de fracture 102 peut également générer, ou finalement produire, des signaux de commande en vue de commander des opérations de terrain associées au volume souterrain. A titre d'exemple les opérations de terrain peuvent être effectuées en utilisant des équipements de forage et d'exploration, des injecteurs et des vannes de commande de puits (ou autres dispositifs de commande lors de la commande réelle du réservoir 204), du réseau de transport et de distribution, des installations de surface, etc. [0036] Par conséquent, un exemple de système 100 peut comprendre un dispositif informatique 200 et une unité d'affichage graphique interactive 232. Le système informatique dans son ensemble peut constituer des simulateurs, des modèles, et l'exemple de moteur de modélisation de contrainte et de fracture 102. Exemple de moteur [0037] La figure 3 représente l'exemple de moteur de modélisation de contrainte et de-fracture 102 de manière plus détaillée que sur la figure 1 et la figure 2. Le mode de réalisation illustré est un exemple de configuration présenté à des fins de description, en vue d'introduire des caractéristiques et des composants d'un moteur mettant en oeuvre l'exemple de modélisation de contrainte et de fracture en utilisant des failles de frottement. Les composants illustrés sont des exemples. Des configurations ou combinaisons de composants différentes de celles qui sont représentées peuvent être utilisées pour mettre en oeuvre les fonctions de modélisation de la contrainte et des fractures, et des composants différents ou supplémentaires peuvent également être utilisés. Comme cela a été introduit plus haut, l'exemple de moteur de modélisation de contrainte et de fracture 102 peut être mis en oeuvre dans des matériels, ou dans des combinaisons de matériels et de logiciels.
Les composants illustrés sont couplés par voie de communication les uns aux autres pour communiquer de la manière souhaitée. Des flèches sont représentées pour suggérer un flux de traitement ou un flux de données, étant donné que les composants peuvent communiquer les uns avec les autres de la manière souhaitée. [0038] L'exemple de moteur de modélisation de contrainte et de fracture 102 illustré sur la figure 3 comprend un tampon destiné aux ensembles de données 302 ou au moins un accès aux ensembles de données 302, un moteur d'initialisation 304, un sélecteur de paramètres de contrainte et de frottement 308, un moteur d'évaluation du coût 310, et un tampon ou une sortie pour les résultats 312. Ces composants sont représentés à des fins de description. D'autres composants, ou un autre agencement des composants, peut conduire à diverses implémentations de l'exemple de moteur de modélisation de contrainte et de fracture 102. Les fonctionnalités de l'exemple de moteur de modélisation de contrainte et de fracture 102 vont être décrites ci-dessous. Fonctionnement de l'exemple de système et de moteur [0039] La figure 4 représente divers procédés de restitution de la paléocontrainte. La figure 4(a) représente une technique classique utilisant des distributions de glissement sur des failles en tant qu'entrée de données, mais-- fournissant des résultats instables. La technique de la figure 4(a) ne permet pas d'augmenter les distributions de glissement avec d'autres formes d'entrées de données, telles que des données GPS, etc. La figure 4(b) représente une méthode de Monte Carlo 2D utilisant ou non des failles de frottement. L'exemple de procédé représenté sur la figure 4(b) peut être mis en oeuvre par l'exemple de - moteur de modélisation de contrainte et de fracture 102. A titre d'exemple, dans une forme de réalisation, le moteur d'initialisation 304 génère la matrice de modèle associée à un volume souterrain 202. Le sélecteur de paramètres de contrainte et de frottement 308 sélectionne de manière itérative des valeurs de paramètres aléatoires (par exemple un rapport de contrainte et un angle d'orientation) pour la contrainte de champ lointain à appliquer au modèle en vue de minimiser un coût au niveau du moteur d'évaluation du coût 310. En outre, dans un ou plusieurs modes de réalisation, une valeur de l'attribut de frottement (par exemple le coefficient de frottement du glissement des failles) est affectée par un utilisateur lors des étapes d'application itérative du modèle. En variante, dans un ou plusieurs modes de réalisation, l'attribut de frottement (par exemple le coefficient de frottement du glissement des failles en association avec le rapport de contrainte et l'orientation) est attribué de manière aléatoire à une valeur par le sélecteur de paramètres de contrainte et de frottement 308 lors des étapes d'application itérative du modèle. Bien que 10 certains procédés d'inversion de la paléocontrainte puissent appliquer un scénario mécanique complet, le moteur de modélisation de contrainte et de fracture 102 améliore certaines techniques en utilisant de multiples types de données et en intégrant un frottement de glissement sur les failles. Les ensembles de données 302 à utiliser sont généralement de deux types : ceux qui fournissent des informations d'orientation (comme les fractures, les 15 plans de faille secondaire avec un angle de frottement interne, et les striures de failles, etc.), et ceux qui fournissent des informations de grandeur (comme le glissement de faille, des données GPS, des données InSAR, etc.). Certains procédés d'inversion de la paléocontrainte sont calculés en utilisant des mesures de glissement sur des plans de faille. [0040] La technique représentée sur la figure 4(a) exécute une opération qui s'effectue en 20 deux stades : (i) résolution du tenseur de contrainte distant initial cyR sur les éléments de faille ne comportant pas de données de déplacement relatif et résolution pour des déplacements relatifs inconnus bj et, (ii) utilisation des déplacements relatifs calculés et connus pour résoudre cyR. On utilise un solveur itératif effectuant un cycle entre les stades (i) et (ii) jusqu'à convergence. 25 [00411 La technique représentée sur la figure 4(b), qui constitue le procédé principal décrit ici, est fondée sur un algorithme de Monte Carlo utilisant des failles de frottement. Dans ce diagramme, n et p sont un nombre d'éléments triangulaires modélisant respectivement les failles et le nombre de points de données. Pour une simulation donnée, une contrainte de champ lointain aléatoire GR est choisie et la discontinuité de déplacement correspondante b sur les failles est calculée. Dans un ou plusieurs modes de réalisation, le coefficient de frottement du glissement est attribué à une valeur constante dans l'ensemble du procédé de minimisation itérative de la fonction de coût. Dans un ou plusieurs modes de réalisation, le coefficient de frottement de glissement est attribué à une valeur aléatoire en association avec la contrainte de champ lointain o-R et les deux sont inversés ensemble via le procédé de minimisation itérative de la fonction de coût. Un post-traitement est ensuite effectué aux points de données et, selon le type de mesures, des fonctions de coût sont calculées en utilisant le champ de déplacement, de déformation, ou de contrainte. Selon ce scénario, pour des centaines de milliers de simulations, le meilleur coût (proche de zéro) est retenu 10 en tant que solution. [0042] Une technique numérique permettant de mettre en oeuvre les exemples de procédés est décrite ci-dessous. Un tenseur distant réduit utilisé pour la simulation et une estimation de la complexité sont également décrits. [0043] Dans une forme de réalisation, une formulation appliquée par le moteur de 15 modélisation de contrainte et de fracture 102 (par exemple le moteur d'initialisation 304) peut être exécutée en utilisant IBEM3D, successeur de POLY3D (POLY3D est décrit par F. Maerten, P. G. Resor, D. D. Pollard, et L. Maerten, Inverting for slip on threedimensional fault surfaces using angular dislocations, Bulletin of the Seismological Society of America, 95:1654-1665, 2005, et par A. L. Thomas, Poly3D: a three 20 dimensional, polygonal element, displacement discontinuity boundary element computer program with applications to fractures, faults, and cavities in the earth's crust, Thèse de Master, Université de Stanford, 1995). IBEM3D est un code à éléments de frontière fondé sur la solution analytique d'une dislocation angulaire dans un espace entier ou un demi-espace élastique homogène ou inhomogène. Un solveur itératif est utilisé pour des raisons 25 de vitesse et à des fins de parallélisation sur des architectures multi-coeurs (Voir par exemple F. Maerten, L. Maerten, et M. Cooke, Solving 3d boundary element problems using constrained iterative approach, Computational Geosciences, 2009). Des contraintes d'inégalité sont ici utilisées pour imposer un frottement de glissement à des surfaces de failles. En d'autre termes, IBEM3D peut être utilisé ainsi pour déterminer (ou calculer) la 30 dépendance de la matrice de modèle M vis-à-vis du coefficient de frottement de glissement. Dans le code sélectionné, les failles sont représentées par des surfaces triangulaires à déplacement et à coefficient de frottement discontinus. Cela offre l'avantage que les surfaces de faille tridimensionnelles approximent plus précisément les surfaces planes incurvées et de lignes d'extrémité incurvées sans introduire de chevauchements ou d' interst ices.
100441 Des conditions aux limites mixtes peuvent être prescrites, et lorsque des conditions aux limites de traction sont spécifiées, le moteur d'initialisation 304 donne une solution pour des composantes de Burgers inconnues. Conformément au présent procédé, des conditions aux limites de traction sont imposées pour les trois axes locaux des éléments triangulaires (axes d'inclinaison, de glissement et normal), étant donné que les failles sont du type à frottement (voir Maerten et al., Solving 3D boundary element problems using constrained iterative approach, Computational Geosciences, 2010). Après que le modèle de contrainte de champ lointain est construit et que la matrice de modèle a été précalculée, il est possible de déterminer (ou calculer) en tout point (au sein de l'espace entier ou du demi- espace, le déplacement), la déformation ou la contrainte en des points d'observation, sous la forme d'un post-traitement. Plus précisément, le champ de contrainte en tout point d'observation est donné par le champ de contrainte perturbé dû aux failles de glissement plus la contribution de la contrainte distante. Par conséquent, l'obtention du champ de contrainte perturbé du fait du glissement sur les failles n'est pas suffisante. En outre, l'estimation du glissement de faille à partir des interprétations sismiques est donnée le long de la direction d'inclinaison. Rien n'est connu le long de la direction de glissement, et un scénario mécanique complet peut être mis en jeu pour reconstituer les composantes inconnues du vecteur de glissement, étant donné qu'il impactera le champ de contrainte perturbé. Le fait de modifier la contrainte de champ lointain (orientation et/ou grandeurs relatives) modifie la distribution du glissement et par conséquent, le champ de contrainte perturbé. Le champ de contrainte perturbé dépend également du coefficient de frottement du glissement qui peut être attribué par l'utilisateur à une valeur fixe tout au long des itérations du procédé de Monte Carlo. Dans un mode de réalisation différent, le coefficient de frottement du glissement, comme l'attribut de contrainte de champ lointain, peut être attribué à une valeur aléatoire dans l'ensemble des itérations de la méthode de Monte Carlo. Il en résulte que la grandeur et l'orientation de la distribution du glissement sont modifiées de manière correspondante. En général, un code tel que IBEM3D est bien adapté au calcul de la totalité des vecteurs de déplacement sur les failles, et peut être optimisé en utilisant une technique à base de matrice H. L'inconnue associée aux fms de la modélisation reste l'estimation de la contrainte de champ lointain devant être imposée en tant que conditions aux limites lors du processus d'inversion. Dans un ou plusieurs modes de réalisation, le rapport de contrainte et l'orientation de la contrainte de champ lointain sont déterminées dans le processus d'inversion sur la base d'un ensemble donné de surfaces de failles en utilisant la méthode de Monte Carlo mentionnée précédemment. Plus précisément, l'ensemble de surfaces de failles est décrit ci-dessous comme étant une géométrie de failles connue alors que la contrainte de champ lointain est modélisée par des formulations mathématiques décrites ci-après. [0045] Dans un exemple de procédé qui peut être mis en oeuvre par le moteur de modélisation de contrainte et de fracture 102, un modèle composé de multiples surfaces de faille est soumis à un coefficient de frottement du glissement et à un tenseur de contrainte de champ lointain constant OER défini dans le système de coordonnées global par l'équation (2) ci-après en supposant une contrainte de champ lointain sub-horizontale (mais la présente méthodologie n'est pas limitée à ce cas). [0046] Dans la formulation présentée ci-après, 0 est défini comme étant une valeur comprise entre 0 et n, plutôt qu'entre -7r/2 et 7r/2. En outre, 0 = 0 correspond au nord et l'angle est défini en sens horaire. RT (2) cYj 13,(p,Ar Etant donné que l'ajout d'une contrainte hydrostatique ne modifie pas a-R, le tenseur de contrainte de champ lointain o-R peut être écrit comme indiqué à l'équation (3), où Ro est la matrice de rotation suivant l'axe vertical (sens horaire), avec 0 e [0,n] : = crR = - av RTe (3) [0047] Lors de l'utilisation du mode de réalisation ci-dessus, la définition d'une contrainte régionale comporte trois inconnues, à savoir (01 - av), (GH - Gy), et O. Lorsqu'on exprime l'équation (2) en utilisant ai, G2 et G3 pour les trois régimes de failles d'Anderson (Anderson, E., The dynamics of faulting. Edinburgh Geol. Soc., 1905, 8(3):387-402), en factorisant avec (G1-453) et en introduisant le rapport de contrainte R = (G2-o3)/(a1-c53) E [0,1], l'équation (4) suivante donne le résultat suivant : - cy3)Ro -1 R1 RT pour un régime de faille normale e [ - a3)R9 [ -R 1- R RT pour un régime de faille décrochante (4) - )R0[ R 1 RT pour un régime de faille chevauchante 10 [0048] En remplaçant R par R' comme indiqué dans l'équation (5) suivante, un paramètre de forme de contrainte unique R' est créé pour les trois régimes de failles combinés : R E [0, 1] pour un régime de faille normale 2- R E [1, 2] pour un régime de faille décrochante (5) 2 + R E [2, 3] pour un régime de faille chevauchante, [0049] Dans un ou plusieurs modes de réalisation, alors que R est désigné comme étant le rapport de contrainte, R' incorpore le rapport de contrainte R et le régime de faille 15 combinés. Lorsqu'on omet le facteur d'échelle (cri - cy3), le tenseur de contrainte régionale dans (23) est défini par deux paramètres, 0 et R' conformément aux limites suivantes : O < Re < 3 O < 0 < 7r (6) [0050] A partir de (0, R'), il est maintenant possible de déterminer les composants de la contrainte de champ lointain : GR = R' 20 (0,R')-40-oo, croi, an) avec : { oG,,, = sin fl - cogi9cos28 au = COe )3 Sin e co. e = sin I? - cos ,8 sin.= e et p étant un angle défini par : ( 13 --= tan-1 pour un régime de faille normale R' y 13 = tan- R-1) pour un régime de faille décrochante p = tan-1 1 pour un régime de faille chevauchante 3- R'j (7) (8) [0051] En outre, l'algorithme permettant de déterminer la contrainte de champ lointain en utilisant une faille de frottement est décrit dans l'ALGORITHME (1) ci-après : ALGORITHME (1): Estimation de la paléocontrainte en utilisant des failles de frottement 10 Entrée : Géométrie des failles + coefficient de frottement du glissement Entrée : Ensemble de données Sortie : o-Ril la paléocontrainte estimée 15 Initialiser: Déterminer (ou calculer) et stocker la matrice de modèle telle quelle ou en utilisant une approximation telle que la matrice H (voir Maerten F. et al., Solving 3D boundary element problems using constrained iterative approach, Computational Geosciences, 2010) Soit c = 1 // coût initial 20 pour i = 1 à np faire Générer aléatoirement 0 e [0,7r] Générer aléatoirement R' e [0,3] Construire la contrainte de champ lointain (Gap, am, ail) en utilisant les équations (7) et (8) Appliquer la contrainte de champ lointain à la matrice de modèle comme indiqué dans l'équation (1) et déterminer (ou calculer) le glissement de faille et la contrainte/la déformation/le déplacement aux points de données // Déterminer (ou calculer) le coût correspondant, en utilisant par exemple la fonction de coût décrite ci-après d = coût(jeu de données) si d < c alors GR= tenseur de contrainte donné par (R',0) c = d fin fin [0052] Bien que la présentation ci-dessus décrive un exemple de formulation, d'autres formulations peuvent être utilisées sans que l'on s'écarte du cadre des revendications. A titre d'exemple, le coefficient de frottement du glissement peut être attribué à une valeur constante (par exemple à 0,4) pour l'Algorithme (1). Dans cet exemple, l'Algorithme (1) est utilisé pour inverser la contrainte de champ lointain GR. Dans un autre exemple, en plus de générer aléatoirement 0 et R', le coefficient de frottement du glissement peut être assigné à une valeur aléatoire pour chaque itération de l'Algorithme (1). Dans cet exemple, l'Algorithme (1) est utilisé pour inverser à la fois la contrainte de champ lointain GR et le coefficient de frottement du glissement. Ensemble de données [0053] La particularité de ce procédé considère que de nombreux types différents d'ensembles de données 302 peuvent être utilisés pour contraindre l'inversion. Certaines données peuvent ne comprendre que des informations d'orientation, tandis qu'un autre ou d'autres type(s) de données peuvent en outre comprendre des informations de déplacement et/ou de grandeur de la contrainte. [0054] Pour des fractures ouvertes (par exemple, des diaclases, des filons, des dikes), l'orientation de la normale au plan de fracture indique la direction de moindre contrainte de compression en 3D (a3). De même, les normales aux couches de solution sous pression et de stylolites indiquent la direction de la contrainte de plus forte compression (ai). L'utilisation de mesures des orientations de fractures, de couches de solution sous pression, et de stylolites en tant qu'ensembles de données 312 permet au moteur de modélisation de contrainte, de fracture et d'activité de faille 102 de reconstituer le régime tectonique ayant généré ces caractéristiques. [0055] La figure 5 représente les plans de fracture et de faille conjuguée. La figure 5(a) indique l'orientation de o-3 par rapport à une fracture ouverte (diaclases, filons, dikes) donnée par sa normale e en 3D. La figure 5(b) est une représentation semblable à la figure 5(a) excepté qu'elle concerne une orientation de cy3par rapport à une diaclase donnée par sa normale projetée e (par exemple, sa trace sur un affleurement). Les figures 5(c) et 5(d) montrent la même représentation que les figures 5(a) et 5(b) excepté pour un stylolite. La figure 5(e) représente l'orientation de cr2 et ai par rapport à des plans de faille conjuguée donnés par l'un de la normalefi,' en 3D et de l'angle de frottement interne 0 (cet angle de frottement n'est pas lié au frottement du glissement sur la faille). Le but à atteindre est de déterminer la meilleure coupe de la contrainte de champ de lointain o-R, et par conséquent les 0 et R', compte tenu de certaines orientations des plans de fracture ouverte pour lesquels les normales coïncident avec les directions de la contrainte de moindre compression o31' en P, ou, de manière équivalente, pour lesquels le plan de la fracture contient la contrainte de plus forte compression (oei), comme illustré sur les figures 5(a) et 5(b). [0056] Des exemples de fonction de coût fondés sur divers types différents d'ensembles de données 302 et géométries des failles illustrée sur la figure 5 sont décrits ci-après. Ut;ilisation de données associées à des diaclases [0057] La fonction de coût pour les diaclases est définie par : (22) où "." est le produit scalaire, 71. est la normale à un plan de fracture, et m est le nombre de points d'observation. [0058] De même, pour des couches de solution sous pression et des stylolites, la fonction de coût est définie comme dans l'équation (22) en utilisant la contrainte de moindre compression o-3 comme dans l'équation (24) : ((ce v (24) Utilisation de données associées à des plans de faille secondaire [0059] L'orientation de plans de faille secondaire qui se développent au voisinage de failles actives plus grandes peut être estimée en utilisant un critère de rupture de Coulomb, défini par l'équation (25): tan(20) = 1/y (25) où 0 est l'angle des plans de rupture par rapport à la contrainte de compression principale maximum ai et y est le coefficient de frottement interne. Deux plans de rupture conjugués se coupent le long de a2 et l'orientation de la faille est influencée par l'orientation des contraintes principales et la valeur du frottement. [0060] La fonction de coût est donc définie par l'équation (26) : (26) où ai est la direction de la contrainte de plus forte compression et o-2 est la direction de la contrainte principale intermédiaire. Le premier terme du membre de droite de l'équation (26) conserve une orthogonalité entre le o-2 calculé et la normale au plan de la faille, tandis que le second terme fait en sorte que l'angle entre le o calculé et le plan de faille est proche de 0 (voir figure 5(E)). Utilisation de données associées à des striures de failles [0061] Dans le cas des striures de failles, la fonction de coût est définie comme indiqué à l'équation (27): !en - (27) où do' and d représentent respectivement le vecteur de glissement normalisé provenant d'une simulation et le vecteur de glissement mesuré. Utilisation de données GPS 10 [0062] Dans le cas d'un ensemble de données GPS, la fonction de coût est définie dans l'équation (29) : (29) où /If est l'élévation mesurée normalisée globalement modifiée au point P par rapport à 15 l'horizon, et up est la modification d'élévation calculée normalisée globalement pour une contrainte de champ lointain donnée. Le premier terme du membre de droite de l'équation (29) représente une minimisation de l'angle entre les deux vecteurs de déplacement, tandis que le second terme représente une minimisation de la différence de la norme.
20 Utilisation de données InSAR [0063] Lorsqu'on utilise un ensemble de données InSAR, il existe deux possibilités. Soit on calcule les vecteurs de déplacement globaux des mesures en utilisant le déplacement u le long de la direction de la ligne de visée du satellite g., auquel cas on utilise l'équation (30) : 25 Ut. = U. (30) et le même processus que celui utilisé pour l'ensemble de données GPS (ci-dessus) est appliqué avec le ep, calculé, soit on calcule les vecteurs de déplacement calculés le long de la ligne de visée du satellite, auquel cas on utilise l'équation (31) : C -T (31) Up = U. S où "." est le produit scalaire. La fonction de coût est par conséquent donnée par l'équation (32) : !'Isar = E( (32) Utilisation de données associées à un horizon aplani [0064] En utilisant le plan moyen d'un horizon sismique donné (horizon aplani), le moteur de modélisation de contrainte et de fracture 102 détermine (ou calcule) tout d'abord la modification de l'élévation pour chaque point constituant l'horizon. On utilise ensuite la fonction de coût GPS pour laquelle la composante uz est obtenue, cela donnant l'équation (33) : (33) [0065] Si un pré-plissement ou un post-plissement de la zone est observé, le plan moyen ne peut plus être utilisé en tant qu'intermédiaire. Par conséquent, une surface de coupe lisse et continue doit être construite qui élimine les déformations génératrices de failles tout en conservant les plissements. Le même processus que pour le plan moyen est ensuite utilisé pour estimer la paléocontrainte. Dans certaines circonstances et avant de définir la surface de coupe continue, un procédé peut filtrer l'horizon d'entrée des bruits à hautes fréquences, tels que les ondulations et les bosses, tout en conservant les déformations naturelles. Utilisation d'informations de rejet en profondeur 100661 Lorsque des données de rejet en profondeur sont utilisées, la fonction de coût est définie comme indiqué dans l'équation (34) : în (34) L où le est la grandeur de rejet en profondeur mesurée pour un élément triangulaire e, et b: est la grandeur de rejet en profondeur calculée.
10 Utilisation de multiples types d'informations disponibles 100671 L'exemple de moteur de modélisation de contrainte et de fracture 102 peut combiner les fonctions de coût précédemment décrites pour mieux contraindre l'inversion de contrainte en utilisant des données disponibles (par exemple toutes les données disponibles, y compris, sans limitation, des données d'orientation de faille et de fracture, 15 des données GPS, des données InSAR, des données d'horizons aplanis, des mesures de rejet en profondeur obtenues par réflexion sismique, des striures de failles, etc.). En outre, les données peuvent être pondérées différemment et un poids peut également être affecté à chaque donnée pour chaque coordonnée. Conclusion et perspectives 20 100681 L'exemple de moteur de modélisation de contrainte et de fracture 102 utilisant des failles de frottement s'applique à des modèles complexes. En outre, la formulation exécutée par l'exemple de moteur de modélisation de contrainte et de fracture 102 permet d'effecteur une inversion de la paléocontrainte en utilisant de multiples types de données telles que orientation de fracture, plans de faille secondaire, GPS, InSAR, rejet de faille, et de stries 25 de glisseeent de faille. Dans une forme de réalisation, en utilisant l'orientation de fracture et/ou les plans de faille secondaire provenant de puits de forages, le moteur de modélisation de contrainte et de fracture 102 restitue un ou plusieurs événements tectoniques, le tenseur de contrainte restitué étant donné par l'orientation et le rapport des grandeurs principales. L'exemple de moteur de modélisation de contrainte et de fracture 102 et les procédés associés peuvent être appliqués dans une large gamme d'applications, comprenant : l'interpolation de la contrainte dans un réservoir faillé de manière complexe, la prédiction de fractures, le contrôle de qualité sur les failles interprétées, le calcul en temps réel de champs de contrainte et de déplacement perturbés tout en effectuant de manière interactive une estimation de paramètres, la prédiction de fracture, le discernement de la fracturation induite de fractures préexistantes, etc.
100691 En variante, une autre application du moteur de modélisation de contrainte et de fracture 102 et des procédés associés consiste en une évaluation du champ de contrainte perturbé (et par conséquent d'unid'événement(s) tectonique(s)) pour la récupération de "gaz de schistes." Comme les schistes présentent une faible perméabilité de la matrice, la production de gaz en quantités commerciales peut faire intervenir des fractures en vue de créer une perméabilité. Cela peut être effectué par fracturation hydraulique afin de créer des fractures artificielles étendues au voisinage de puits de forage, et implique donc une bonne compréhension de la manière dont les fractures se propageront en fonction du champ de contrainte perturbé. Exemples de procédés 100701 Bien que les divers blocs de ces organigrammes soient présentés et décrits de manière séquentielle, l'homme du métier notera que certains ou tous les blocs peuvent être exécutés dans des ordres différents, peuvent être combinés ou omis, et que certains ou tous les blocs peuvent être exécutés en parallèle. En outre, les blocs peuvent être exécutés de manière active ou passive. A titre d'exemple, certains blocs peuvent être exécutés en utilisant une scrutation ou peuvent être pilotés par interruption conformément à un ou plusieurs modes de réalisation. A titre d'exemple, des blocs de détermination peuvent ne pas faire intervenir de processeur pour traiter une instruction sans qu'une interruption ait été reçue pour signifier qu'une condition existe, conformément à un ou plusieurs modes de réalisation. A titre d'autre exemple, des blocs de détermination peuvent être exécutés en effectuant un test, par exemple en vérifiant une valeur de données pour tester si la valeur concorde avec la condition testée conformément à un ou plusieurs modes de réalisation. [0071] La figure 6 représente un exemple de procédé 1300 de modélisation de contrainte et de fracture utilisant des failles de frottement. Dans l'organigramme, les opérations sont résumées sous forme de blocs individuels. L'exemple de procédé 1300 peut être mis en oeuvre par des matériels ou des combinaisons de matériels et de logiciels, par exemple, par l'exemple de moteur de modélisation de contrainte et de fracture illustré sur la figure 3 ci- dessus [0072] Au bloc 1302, la géométrie des failles pour un volume de terre souterrain est reçue. [0073] Au bloc 1304, au moins un ensemble de données associé au volume de terre souterrain est également reçu. Dans un ou plusieurs modes de réalisation, l'ensemble de données comprend une mesure du volume souterrain, par exemple des données d'interprétation sismique, des données de puits de forage et/ou des données d'observation sur le terrain. Plus précisément, la mesure du volume souterrain peut comprendre des données de géométrie de faille, des données d'orientation de fracture, des données d'orientation des stylolites, des données de plan de faille secondaire, des données de rejet de faille, des données de stries de glissement, des données GPS, des données InSAR, des données de télémétrie laser, des données d'inclinomètre, des données de déplacement pour une faille géologique et/des données de grandeur de la contrainte pour la faille géologique. [0074] Au bloc 1306, un modèle du volume souterrain est obtenu. Dans un ou plusieurs modes de réalisation, le modèle comprend une matrice de modèle représentant une relation entre un résultat de glissement de faille modélisé généré par le modèle et une condition aux limites appliquée au modèle. A titre d'exemple, le modèle peut être fondé sur l'équation (1) décrite précédemment, dans laquelle la condition aux limites comprend un attribut de contrainte régionale et un attribut de frottement de faille. [0075] Au bloc 1310, un segment de post-traitement du procédé commence, qui peut déterminer (ou calculer) de nombreux résultats temps réel en utilisant des failles de frottement. [0076] Au bloc 1312, des valeurs de paramètres aléatoires pour l'attribut de contrainte régionale (par exemple, R' et 0 décrits précédemment) sont sélectionnées. Dans un ou plusieurs modes de réalisation, une valeur de l'attribut de frottement de faille (par exemple le coefficient de frottement du glissement) est également sélectionnée. [0077] Au bloc 1314, le résultat de glissement de faille modélisé est calculé au moyen du modèle fondé sur les valeurs sélectionnées aléatoirement de l'attribut de contrainte. A titre d'exemple, les valeurs sélectionnées aléatoirement sont sélectionnées par un processeur d'ordinateur en utilisant la méthode de Monte Carlo. Dans un ou plusieurs modes de réalisation, le résultat de glissement de faille modélisé est calculé sur la base d'une valeur sélectionnée par l'utilisateur de l'attribut de frottement. Dans un ou plusieurs modes de réalisation, le résultat de glissement de faille modélisé est calculé sur la base d'une valeur sélectionnée aléatoirement de l'attribut de frottement, qui est également sélectionnée par le processeur d'ordinateur en utilisant la méthode de Monte Carlo. Dans un ou plusieurs modes de réalisation, le résultat de glissement de faille modélisé comprend un champ de déformation, de contrainte et/ou de déplacements perturbé, ainsi que d'autres attributs dépendant des failles de frottement du volume souterrain. Dans un ou plusieurs modes de réalisation, un coût est calculé pour représenter une différence entre les données modélisées (par exemple le résultat de glissement de faille modélisé) et les données mesurées (par exemple une mesure du volume de terre souterrain). [0078] Au bloc 1316, le coût associé au résultat de glissement de faille modélisé venant d'être calculé est évalué. Si le coût n'est pas satisfaisant, le procédé reboucle ensuite au bloc 1312 pour sélectionner de nouvelles valeurs de paramètres aléatoires. Si le coût est satisfaisant, le procédé se poursuit ensuite au bloc 1318. [0079] Au bloc 1318, l'attribut de contrainte régionale provenant du bloc 1314 est appliqué au volume de terre souterrain en vue de générer des résultats de modélisation, par exemple, par rapport à une requête concernant le volume de terre souterrain ou en réponse à un paramètre mis à jour concernant le volume de terre souterrain. Dans un ou plusieurs modes de réalisation, les résultats de modélisation sont générés en temps réel. En d'autres termes, les résultats de modélisation sont générés en un temps prédéterminé (par exemple une seconde, une minute, une heure, etc.) après réception de la requête ou du paramètre mis à jour concernant le volume de terre souterrain. A titre d'exemple, les résultats de modélisation peuvent comprendre un ou plusieurs d'une inversion de contrainte, d'un champ de contrainte, d'une valeur de contrainte de champ lointain, d'une interpolation de contrainte dans un réservoir faillé complexe, d'un champ de contrainte perturbé, d'un rapport de contrainte et d'une orientation associée, d'un ou plusieurs événements tectoniques, d'une discontinuité de déplacement d'une faille, d'un glissement de faille, d'un déplacement estimé, d'une déformation perturbée, d'une distribution du glissement sur les failles, d'un contrôle de qualité concernant des failles interprétées, d'une prédiction de fracture, d'une prédiction de la propagation de fracture conformément à un champ de contrainte perturbé, d'un calcul en temps réel de champs de contrainte et de déplacement perturbés tout en effectuant une estimation interactive de paramètres, ou un discernement d'une fracture induite d'une fracture préexistante. [0080] Au bloc 1320, une requête ou un paramètre mis à jour concernant le volume de terre souterrain est reçu, qui fait naitre ou initie la génération des résultats de modélisation dans la section de résultats temps réel (1310) du procédé 1300. [0081] La figure 7 et la figure 8 représentent un exemple de modélisation de contrainte, de fracture et d'activité de faille. Plus précisément, la figure 8 représente une vue en perspective du volume de terre souterrain 801 comportant de multiples failles 802 représentées sous la forme de surfaces triangulaires. Aux fins de la modélisation, le volume de terre souterrain 801 est quadrillé en un nombre de points d'observation. Ces points d'observation se situant sur la section transversale 803 sont représentés sous forme de points. Bien que cela ne soit pas explicitement représenté autrement que sur la section transversale 803, les points d'observation peuvent exister dans l'ensemble du volume de terre souterrain 801. Comme décrit ci-dessus, des valeurs de contrainte, de déformation et/ou de déplacement peuvent être calculées pour chacun de ces points d'observation. [0082] En outre, comme illustré sur la figure 8, une vue de dessus A 804 correspond à la section transversale 803 et représente l'orientation calculée de diaclases A 806 qui est représentée sous la forme de lignes du type empreinte digitale dans l'ensemble de la section transversale 803. Comme décrit ci-dessus, pour des fractures ouvertes (par exemple, des diaclases, des filons, des dikes), l'orientation de la normale au plan de la fracture indique la direction de la contrainte de moindre compression (c'est-à-dire a3). Plus précisément, l'orientation calculée de diaclases A 806 est basée sur des valeurs optimisées de 0 et de R' en utilisant le procédé de Monte Carlo, comme décrit en référence à la figure 6 ci-dessus, dans laquelle aucun frottement n'est pris en compte pour les multiples failles 802. En d'autres termes, le coefficient de frottement du glissement y est mis à 0 dans la matrice de modèle. Dans un ou plusieurs modes de réalisation, le coefficient de frottement du glissement y est attribué à 0 par un utilisateur. Dans un ou plusieurs modes de réalisation, le coefficient de frottement du glissement y peut être inclus dans les variables aléatoires de Monte Carlo (en association avec 0 et R') et la valeur 0 attribuée par le sélecteur de paramètres de contrainte et de frottement 308 représenté sur la figure 3 ci- dessus. En outre, la vue de dessus A 804 est superposée avec des segments en gras représentant des orientations observées de diaclases de la diaclase observée A (805-1), de la diaclase observée B (805-2), de la diaclase observée C (805-3), et de la diaclase observée D (805-4). L'ensemble des orientations de ces diaclases observées constitue un exemple des ensembles de données 302 illustrés sur la figure 3 ci-dessus. [0083] Une comparaison de l'orientation calculée de diaclases A 806 et des segments en gras représentant des orientations de diaclases observées montree l'accord entre l'orientation calculée et l'orientation observée pour la diaclase observée A (805-1) et la diaclase observée D (805-4). Cependant, l'orientation calculée et l'orientation observée sont en désaccord grossier pour la diaclase observée B (805-2) et la diaclase observée C (805-3). Conformément à cette comparaison visuelle, le moteur d'évaluation du coût 310 représenté sur la figure 3 ci-dessus détermine que la sélection de la valeur "0" pour le coefficient de frottement du glissement correspondant aux failles 802 ne conduit pas au coût le plus faible, même en choisissant toutes les combinaisons possibles de R' et O. En conséquence, une valeur différente du coefficient de frottement du glissement est sélectionnée, comme illustré sur la figure 8 ci-après. Dans un ou plusieurs modes de réalisation, la valeur différente est choisie manuellement par l'utilisateur. Dans un ou plusieurs modes de réalisation, la valeur différente est choisie par le sélecteur de paramètres de contrainte et de frottement 308 conformément à la méthode de Monte Carlo. La figure 8 montre que la vue de dessus B 808 correspond à la section transversale 803 et représente l'orientation calculée de diaclases B 807 sur la base du coefficient de frottement du glissement qui est mis à la valeur de 0,6. Une comparaison de l'orientation calculée de diaclases B 808 et des segments en gras représentant des orientations de diaclases observées fait apparaître l'accord entre l'orientation calculée et l'orientation observée pour toutes les diaclases, à savoir la diaclase observée A (805-1), la diaclase observée B (805-2), la diaclase observée C (805-3), et la diaclase observée D (805-4). Conformément à cette comparaison visuelle, le moteur d'évaluation du coût 310 détermine que la sélection de la valeur de 0,6 pour le coefficient de frottement du glissement, tout en optimisant les valeurs de 0 et R', conduit au coût le plus faible. En conséquence, les valeurs optimisées de 0 et R' et la sélection particulière de la valeur 0,6 pour le coefficient de frottement du glissement sont utilisées lors de la génération des résultats 312 illustrés sur la figure 3 ci-dessus. Conclusion [0084] Bien que des exemples de systèmes et de procédés aient été décrits selon un langage propre aux caractéristiques structurales et/ou à des actions méthodologiques, il est à noter que l'objet de l'invention défini dans les revendications annexées n'est pas limité aux caractéristiques ou aux actions spécifiques décrites. Les caractéristiques et les actions spécifiques sont en effet décrites en tant qu'exemples de formes de réalisation des systèmes, des procédés et des structures revendiqués. [0085] En outre, bien que la description présentée ci-dessus concerne un nombre limité de modes de réalisation, l'homme du métier tirant bénéfice de la présente description notera qu'il est possible de concevoir d'autres modes de réalisation ne s'écartant pas du cadre des revendications. [0086] Bien que seul un petit nombre d'exemples de réalisation aient été décrits ci-dessus en détail, l'homme du métier notera que de nombreuses variantes sont envisageables dans les exemples de modes de réalisation sans que l'on s'écarte matériellement de la présentation de l'inversion de la contrainte en utilisant des failles de frottement. En conséquence, toutes ces variantes doivent être considérées comme entrant dans le cadre de la présente invention, telle qu'elle est définie dans les revendications annexées. Dans les revendications, des clauses portant sur les moyens et les fonctions doivent être considérées comme couvrant les structures décrites ici lorsqu'elles mettent en oeuvre la fonction citée et non seulement les équivalents structuraux, mais également des structures équivalentes. Par conséquent, bien qu'un clou et une vis puissent ne pas être des équivalents structuraux, en ce sens qu'un clou utilise une surface cylindrique pour fixer les unes aux autres des pièces de bois, tandis qu'une vis utilise une surface hélicoïdale, dans l'environnement de la fixation de pièces en bois, un clou et une vis peuvent être des structures équivalentes.