FR3028333A1 - Creation de maillage de reservoir utilisant un raffinement anisotrope prolonge, adapte a la geometrie, d'un polyedre - Google Patents
Creation de maillage de reservoir utilisant un raffinement anisotrope prolonge, adapte a la geometrie, d'un polyedre Download PDFInfo
- Publication number
- FR3028333A1 FR3028333A1 FR1559673A FR1559673A FR3028333A1 FR 3028333 A1 FR3028333 A1 FR 3028333A1 FR 1559673 A FR1559673 A FR 1559673A FR 1559673 A FR1559673 A FR 1559673A FR 3028333 A1 FR3028333 A1 FR 3028333A1
- Authority
- FR
- France
- Prior art keywords
- fracture
- cells
- fractures
- mesh
- refined
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 230000002035 prolonged effect Effects 0.000 title 1
- 238000000034 method Methods 0.000 claims abstract description 141
- 238000001125 extrusion Methods 0.000 claims abstract description 14
- 238000003860 storage Methods 0.000 claims description 35
- 238000004088 simulation Methods 0.000 claims description 19
- 230000015654 memory Effects 0.000 claims description 17
- 238000007670 refining Methods 0.000 claims description 7
- 230000004044 response Effects 0.000 abstract description 4
- 238000004422 calculation algorithm Methods 0.000 description 38
- 238000012360 testing method Methods 0.000 description 24
- 230000008569 process Effects 0.000 description 21
- 238000004891 communication Methods 0.000 description 20
- 238000010586 diagram Methods 0.000 description 14
- 230000006870 function Effects 0.000 description 12
- 238000004590 computer program Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 8
- 230000015572 biosynthetic process Effects 0.000 description 6
- 239000012530 fluid Substances 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 208000035126 Facies Diseases 0.000 description 2
- 241000080590 Niso Species 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000013100 final test Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G01V20/00—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
Abstract
Des systèmes et des procédés sont proposés pour la création des maillages de réservoir utilisant un raffinement anisotrope étendu, adapté à la géométrie, d'un polyèdre. Dans un exemple, un procédé comprend l'identification, basée sur la spécification du réservoir (204), d'un jeu de fractures (206, 208) comprenant des fractures (206) permettant la dimension 2,5 (2,5D), et d'autres fractures (208). Le procédé comprend également la génération d'un modèle de réservoir intermédiaire (212) comprenant un maillage d'extrusion qui modélise les fractures (206) permettant le 2,5D dans un espace en 3D. En réponse à la détermination que les cellules dans le maillage doivent être raffinées dans une direction à l'intérieur de l'espace en 3D, le procédé raffine les cellules de façon anisotrope dans le maillage correspondant à d'autres fractures (208). Le procédé comprend également la résolution d'un réseau de fracture à l'intérieur du modèle de réservoir intermédiaire (212) en utilisant les cellules raffinées, et ensuite en générant un modèle de réservoir de Terre à l'aide du réseau de fracture.
Description
Domaine technique [0001] La présente divulgation concerne généralement la modélisation des réservoirs d'hydrocarbures et, plus particulièrement, des techniques permettant de créer des maillages de réservoir utilisant un raffinement anisotrope prolongé, adapté à la géométrie, d'un polyèdre. Arrière-plan [0002] Dans l'industrie du gaz et du pétrole, la modélisation des réservoirs implique la construction d'un modèle informatique d'un réservoir de pétrole afin d'améliorer l'estimation des réserves et la prise de décision concernant le développement du champ. Par ex., des modèles géologiques peuvent être créés pour obtenir une description statique d'un réservoir avant la production. Des modèles de simulation du réservoir peuvent être créés pour simuler et prédire le flux de fluide dans un réservoir au cours de sa durée de vie de production. [0003] Un problème rencontré avec les modèles de simulation de réservoir est le défi de la modélisation des fractures à l'intérieur d'un réservoir, qui demande une compréhension approfondie des caractéristiques du flux matriciel, de la connectivité du réseau de la fracture et de l'interaction matrice-fracture. Les fractures peuvent être des fissures ouvertes ou des vides à l'intérieur de la formation, et elles peuvent être naturelles ou générées artificiellement à partir d'un puits de forage. La modélisation précise des fractures est importante étant donné que des propriétés de la fracture telles que la distribution spatiale, l'ouverture, la longueur, la hauteur, la conductivité et la connectivité affectent de façon importante le flux des fluides du réservoir vers le puits de forage. [0004] Des techniques de génération de maillage sont utilisées dans la modélisation des réservoirs. Deux techniques traditionnelles de génération de maillage pour la simulation des réservoirs en 3D sont le maillage basé sur la structure et le maillage basé sur l'extrusion. Dans les techniques structurées, des hexaèdres sont connectés dans un espace logique 3D i-j-k dans lequel chaque noeud de maillage intérieur est adjacent à 8 hexaèdres. Les prolongements des techniques structurées comprennent un raffinement local du maillage dans lequel des régions locales d'une grille originale sont remplacées par des grilles plus fines. Ceci peut être chronophage, coûteux en ressources informatiques et un fardeau prohibitif lorsqu'on traite des géométries générales de réservoir, telles que des surfaces de fracture en 3D arbitraires. En raison de la nature 2,5 dimensionnelle (2,5 D) inhérente des techniques d'extrusion existantes, des limitations semblables s'appliquent à ces techniques. Par ailleurs, il existe des techniques de maillage complètement non-structurées, comme par ex., les schémas de maillage par tétraédralisation et polyédriques. La complexité accrue de ces techniques entraîne souvent une diminution de la robustesse en comparaison à des techniques structurées, particulièrement, en présence d'une entrée de géométrie imparfaite (c.-à-d., une « géométrie sale »). Présentation Dans un ou plusieurs mode(s) de réalisation du présent exposé, un procédé implémenté par ordinateur comprend : la réception, au niveau d'un dispositif informatique, d'une spécification de réservoir ; l'identification, basée sur la spécification du réservoir, d'un jeu de fractures comprenant : des fractures permettant la 2,5 dimension (2,5D) ; et d'autres fractures ; la génération d'un modèle de réservoir intermédiaire comprenant un maillage d'extrusion qui modélise les fractures permettant la 2,5D dans un espace tridimensionnel (3D) répondant à la détermination que les cellules dans le maillage doivent être raffinées dans une direction à l'intérieur de l'espace 3D, le raffinement des cellules de façon anisotrope dans le maillage correspondant aux autres fractures ; la résolution, par le dispositif informatique, d'un réseau de fracture à l'intérieur du modèle de réservoir intermédiaire en utilisant les cellules raffinées ; et la génération d'un modèle de réservoir de Terre à l'aide du réseau de fracture. Dans un ou plusieurs mode(s) de réalisation du présent exposé, les fractures permettant la 2,5D possèdent une géométrie qui a été discrétisée dans un plan en 2D par une collection de segments de ligne, et la génération du modèle de réservoir intermédiaire comprend, pour chaque segment de ligne associé à chaque fracture dans les fractures permettant la 2,5D : la génération d'un jeu de stades à un rayon spécifié provenant du segment de ligne ; la génération des boucles fermées autour des segments de ligne de la fracture associés à une fracture permettant la 2,5D ; et la génération des éléments de forme à l'intérieur des boucles fermées du segment de la ligne, le procédé comprenant également : la génération du maillage sous forme de maillage contraint autour des boucles fermées des fractures permettant la 2,5D pour remplir un espace restant dans le plan en 2D. Dans un ou plusieurs mode(s) de réalisation du présent exposé, le procédé comprend également : l'utilisation d'un modèle de réservoir de Terre dans une simulation de réservoir ; et l'affichage de la simulation de réservoir sur un écran du dispositif informatique. Dans un ou plusieurs mode(s) de réalisation du présent exposé : les fractures permettant le 2,5D peuvent comprendre des fractures représentées par un ou plusieurs plans sensiblement verticaux ; et les autres fractures comprennent des fractures représentées par un ou plusieurs plans sensiblement horizontaux.
Dans un ou plusieurs mode(s) de réalisation du présent exposé, le raffinement anisotrope des cellules dans le maillage comprend : l'identification d'une direction à l'intérieur de l'espace 3D dans laquelle les cellules doivent être raffinées ; et la division d'un rebord des cellules, le rebord étant dans la direction à l'intérieur de l'espace 3D. Dans un ou plusieurs mode(s) de réalisation du présent exposé, le rebord est divisé le long d'un plan qui est perpendiculaire à un axe d'une fracture à l'intérieur d'un réseau de fracture.
Dans un ou plusieurs mode(s) de réalisation du présent exposé, le raffinement anisotrope des cellules comprend : la détermination d'un emplacement d'une fracture dans un modèle de réservoir intermédiaire ; l'identification des cellules dans la région proche de la fracture ; la détermination d'une taille de cible des rebords des cellules proches de la fracture ; et la détermination si les rebords des cellules proches de la fracture dans la direction à l'intérieur de l'espace 3D sont plus longs que la taille de la cible multipliée par un variable scalaire, et les cellules proches de la fracture ayant des rebords dans la direction à l'intérieur de l'espace 3D qui sont plus longs que la taille de la cible multipliée par le variable scalaire sont parmi les cellules qui sont raffinées. Dans un ou plusieurs mode(s) de réalisation du présent exposé, la détermination de l'emplacement de la fracture comprend : l'échelonnement des rebords des cellules proches de la fracture qui sont plus longs que la taille de la cible multipliée par la variable scalaire ; et la détermination si les rebords échelonnés croisent la fracture, dans lequel les cellules proches de la fracture qui possèdent des rebords échelonnés qui croisent la fracture sont des cellules qui sont raffinées.
Dans un ou plusieurs mode(s) de réalisation du présent exposé, la détermination de l'emplacement de la fracture comprend également : l'échelonnement des rebords des cellules proches de la fracture qui sont plus longs que la taille de la cible multipliée par la variable scalaire ; et pour une quelconque paire opposée des rebords échelonnés, la détermination si exactement un des deux rebords échelonnés croisent une fracture, dans lequel les cellules proches de la fracture qui ont exactement un rebord échelonné croisant la fracture sont les cellules qui sont raffinées.
Dans un ou plusieurs mode(s) de réalisation du présent exposé, la détermination de l'emplacement de la fracture comprend également : la détermination si les rebords de cellules proches de la fracture qui sont plus longs que la taille de la cible multipliée par la variable scalaire comportent deux ou plusieurs noeuds suspendus, dans lequel les cellules proches de la fracture ayant deux ou plusieurs noeuds suspendus sont raffinées.
Dans un ou plusieurs mode(s) de réalisation du présent exposé, le raffinement anisotrope des cellules dans le maillage comprend également, pour une quelconque paire opposée des rebords des cellules proches de la fracture qui sont plus longs que la taille de la cible multipliée par la variable scalaire, la détermination si exactement un des deux rebords comportent deux ou plusieurs noeuds suspendus, et dans lequel les cellules proches de la fracture qui comportent deux ou plusieurs noeuds suspendus sont raffinées. Dans un ou plusieurs mode(s) de réalisation du présent exposé, le raffinement anisotrope des cellules dans le maillage comprend également : l'échelonnement des rebords des cellules proches de la fracture qui sont plus longs que la taille de la cible multipliée par la variable scalaire ; et la détermination si les cellules proches de la fracture ayant des rebords échelonnés croisent une fracture alors qu'aucune cellule proche de la fracture ne croise une fracture, dans lequel les cellules proches de la fracture qui possèdent des rebords échelonnés qui croisent la fracture sont raffinées. Dans un ou plusieurs mode(s) de réalisation du présent exposé, le raffinement anisotrope des cellules dans le maillage comprend également : la détermination d'un emplacement d'une fracture dans un modèle de réservoir intermédiaire ; l'identification des cellules dans la région proche de la fracture ; l'identification d'une direction à l'intérieur de l'espace 3D dans laquelle les cellules proches de la fracture doivent être raffinées ; et la détermination si les cellules proches de la fracture ayant un rebord dans la direction identifiée ont également pas plus d'un noeud suspendu, dans lequel les cellules proches de la fractures qui n'ont pas plus d'un noeud suspendu sont raffinées. Dans un ou plusieurs mode(s) de réalisation du présent exposé, le raffinement anisotrope des 30 cellules dans le maillage comprend également : la détermination d'un emplacement d'une fracture dans un modèle de réservoir intermédiaire ; l'identification des cellules dans la région proche de la fracture ; et pour une quelconque paires opposée des rebords de cellule proche de la fracture, la détermination si exactement l'un des deux rebords opposés possède deux ou plusieurs noeuds suspendus, dans lequel les cellules proches de la fracture qui possèdent les deux ou plusieurs noeuds suspendus sont raffinées. Dans un ou plusieurs mode(s) de réalisation du présent exposé, un support de stockage lisible par ordinateur a des instructions exécutables stockées sur ledit support de stockage, qui, lorsqu'elles sont exécutées par un dispositif informatique, permet au dispositif informatique de réaliser des opérations, les instructions comprenant : des instructions permettant de recevoir une spécification de réservoir ; des instructions pour l'identification, basée sur la spécification du réservoir, d'un jeu de fractures comprenant : des fractures permettant la 2,5 dimension (2,5D) ; et d'autres fractures ; des instructions pour la génération d'un modèle de réservoir intermédiaire comprenant un maillage d'extrusion qui modélise les fractures permettant la 2,5D dans un espace tridimensionnel (3D) ; répondant à la détermination que les cellules dans le maillage doivent être raffinées dans une direction à l'intérieur de l'espace 3D, instructions pour le raffinement anisotrope des cellules dans le maillage correspondant aux autres fractures ; instructions pour la résolution d'un réseau de fracture à l'intérieur du modèle de réservoir intermédiaire en utilisant les cellules raffinées ; et des instructions pour la génération d'un modèle de réservoir de Terre à l'aide du réseau de fracture. Dans un ou plusieurs mode(s) de réalisation du présent exposé, les instructions pour le raffinement anisotrope des cellules dans le maillage comprennent : des instructions pour l'identification d'une direction à l'intérieur de l'espace 3D dans laquelle les cellules doivent être raffinées ; et instructions pour la division d'un rebord des cellules, le rebord étant dans la direction à l'intérieur de l'espace 3D.
Dans un ou plusieurs mode(s) de réalisation du présent exposé : les fractures permettant le 2,5D peuvent comprendre des fractures représentées par un ou plusieurs plans sensiblement verticaux ; et les autres fractures comprennent des fractures représentées par un ou plusieurs plans sensiblement horizontaux. Dans un ou plusieurs mode(s) de réalisation du présent exposé, un système comprend: un processeur ; un dispositif d'affichage ; et une mémoire ayant des instructions stockées sur celle-ci, qui, lorsqu'elles sont exécutées par un processeur, permet au processeur d'exécuter des opérations comprenant : la réception d'une spécification de réservoir ; l'identification, basée sur la spécification du réservoir, d'un jeu de fractures comprenant : des fractures permettant la 2,5 dimension (2,5D) ; et d'autres fractures ; la génération d'un modèle de réservoir intermédiaire comprenant un maillage d'extrusion qui modélise les fractures permettant la 2,5D dans un espace tridimensionnel (3D) répondant à la détermination que les cellules dans le maillage doivent être raffinées dans une direction à l'intérieur de l'espace 3D, le raffinement des cellules de façon anisotrope dans le maillage correspondant aux autres fractures ; la résolution d'un réseau de fracture à l'intérieur du modèle de réservoir intermédiaire en utilisant les cellules raffinées ; et la génération d'un modèle de réservoir de Terre à l'aide du réseau de fracture. Dans un ou plusieurs mode(s) de réalisation du présent exposé: les fractures permettant le 2,5D comprennent des fractures représentées par un ou plusieurs plans verticaux ; et les autres fractures comprenant des fractures représentées par un ou plusieurs plans sensiblement horizontaux.
Dans un ou plusieurs mode(s) de réalisation du présent exposé, le système comprend également : l'affichage, sur le dispositif d'affichage, du modèle de réservoir de Terre.
Brève description des dessins [0005] La FIG. 1 est un organigramme d'un système de raffinement anisotrope prolongé, adapté à la géométrie, permettant la création de maillages de réservoir selon certains exemples illustratifs de la présente divulgation. [0006] La FIG. 2 est un schéma de travail illustrant un procédé permettant de créer des maillages de réservoir utilisant un raffinement anisotrope prolongé, adapté à la géométrie, d'un polyèdre, selon certains exemples illustratifs de la présente divulgation. [0007] La FIG. 3 illustre des exemples de configuration de raffinement pour une cellule d'hexaèdre unique à l'intérieur d'un espace en 3D U-V-W, selon certains exemples illustratifs de la présente divulgation. [0008] Les Figures 4 et 5 illustrent des cellules de prisme triangulaire représentatives montrant un système coordonné logique en 2D U-W et un schéma de numéro pour les rebords des prismes, selon certains exemples illustratifs de la présente divulgation. [0009] La FIG. 6 illustre des configurations de raffinement pour une cellule de prisme unique utilisant un paradigme à direction U unique à l'intérieur du système coordonné logique en 2D U-W, selon certains exemples illustratifs de la présente divulgation. [0010] La FIG. 7 illustre des configurations de raffinement dans 4 directions utilisant un paradigme à multiples directions U à l'intérieur du système coordonné logique en 2D U-W, selon certains exemples de la présente divulgation. [0011] La FIG. 8 illustre une interface utilisateur interactive permettant d'afficher un modèle contenant des prismes hexaèdres et triangulaires, selon certains aspects de la présente divulgation. [0012] La FIG. 9 illustre une interface utilisateur interactive permettant d'afficher des fractures raffinées dans un modèle utilisant un paradigme à direction U unique, selon certains exemples illustratifs de la présente divulgation. [0013] La FIG. 10 illustre un prisme non-triangulaire représentatif (par ex., un élément prismatique général) montrant le système coordonné logique en 2D, selon certaines caractéristiques de la présente divulgation. [0014] La FIG. 11 illustre un raffinement possible d'un élément prismatique général en subdivisant l'élément en hexaèdres, selon certains exemples illustratifs de la présente divulgation. [0015] La FIG. 12 illustre une cellule représentative montrant le système coordonné local en 3D et un schéma de numéro, selon certains aspects illustratifs de la présente divulgation. [0016] Les Figures 13 et 14 illustrent l'application des fonctions de mise à l'échelle à une cellule et à ses rebords, respectivement, selon certaines caractéristiques illustratives de la présente divulgation. [0017] Les Figures 15 et 16 illustrent deux raffinements différents générés en utilisant les procédés illustratifs de la présente divulgation. [0018] Les Figures 17 et 18 sont des exemples d'images de maillages générées utilisant les aspects illustratifs de la présente divulgation. [0019] Les Figures 19 et 20 sont des modèles de Terre complets générés utilisant certains des exemples de procédés de la présente divulgation. [0020] Les Figures 21 et 22 sont des organigrammes illustrant un procédé de raffinement de cellule anisotrope prolongé, adapté à la géométrie, selon certains exemples illustratifs de la présente divulgation. [0021] La FIG. 23 illustre une image de fracture en 3D modélisée à l'aide de certains exemples illustratifs de la présente divulgation. [0022] La FIG. 24 illustre un exemple d'un jeu de fractures discrétisées en 3D croisant un plan en 2D, selon certains aspects illustratifs de la présente divulgation. [0023] La FIG. 25 est un organigramme d'un procédé permettant de modéliser des fractures en 3D, selon certains exemples illustratifs de la présente divulgation. [0024] La FIG. 26 est une illustration de la génération d'un maillage informatisé autour d'un segment d'une ligne de fracture, selon certaines caractéristiques illustratives de la présente divulgation. [0025] La FIG. 27 est une illustration de la génération d'un maillage informatisé autour des segments d'une ligne de fracture croisée, selon certains exemples illustratifs de la présente divulgation. [0026] La FIG. 28 illustre un exemple de maillages informatisés autour d'un réseau complexe de segments d'une ligne de fracture, selon certains exemples illustratifs de la présente divulgation. [0027] La FIG. 29 illustre une grille non-structurée générée autour de géométries complexe qui comprennent des segments croisés de ligne de fracture, selon certains exemples illustratifs de la présente divulgation. [0028] La FIG. 30 est un diagramme d'un exemple d'un système informatique dans lequel certains modes de réalisation de la présente divulgation peuvent être implémentés. Description détaillée [0029] Certains aspects et caractéristiques concernent la génération de maillages informatisés hybrides autour de fractures complexes et discrètes dans le but de simuler des réservoirs. Des exemples illustratifs et des méthodologies apparentées de la présente divulgation sont décrits ci-dessus tels qu'ils peuvent être utilisés dans un système qui applique un raffinement anisotrope prolongé, adapté à la géométrie, d'un polyèdre (AGAR), pour créer un maillage de réservoir. Dans un mode de réalisation, un algorithme AGAR prolongé est utilisé en association avec des mailleurs en 2,5 D, tels qu'un mailleur stadimétrique. Des exemples de mailleurs stadimétrique sont décrits par rapport aux Figs. 24-28 Dans un mode de réalisation, un algorithme AGAR utilise un modèle de la Terre composé entièrement d'hexaèdres comme structure de base dans laquelle il résout les réseaux de fracture qui sont modélisés. L'extension de l'algorithme AGAR permet aux systèmes et procédés décrits ici d'utiliser des topologies plus générales dans un modèle basé sur l'entrée. L'extension AGAR permet d'avoir un flux de travail qui résout les réseaux de fracture en 2 phases : une phase de maillage en 2,5D ; et une phase de raffinement en 3D. Cette combinaison peut maintenir les avantages de la vitesse de calcul et de l'orientation de commande d'un mailleur à 2,5D lorsqu'on le souhaite, tout en invoquant également un algorithme AGAR pour résoudre les réseaux de fracture généraux en 3D. [0030] Un algorithme AGAR prolongé est utilisé de sorte que les modes de réalisation ne soient pas limités à l'utilisation d'un maillage structuré. En particulier, certains modes de réalisation utilisent initialement un mailleur stadimétrique en 2,5D qui procure un mailleur non-structuré, et ensuite l'utilisation d'une version prolongée, améliorée de l'algorithme AGAR pour raffiner les fractures en 3D. Certains exemples utilisent les techniques AGAR pour mettre à l'échelle les rebords et les cellules. Ces techniques peuvent être appliquées directement à des types d'éléments additionnels, tels que, par ex., des prismes triangulaires et d'autres éléments prismatiques généraux. Selon certains modes de réalisation, un jeu de règles est défini pour raffiner les éléments plus généraux et pour s'assurer que le raffinement de la cellule permette une topologie de maillage conforme. [0031] Dans un mode de réalisation, une extension AGAR permet toujours un raffinement anisotrope pour la majorité des éléments dans la plupart des modèles tout en éliminant la restriction qui fait que les éléments puissent être des hexaèdres. Une caractéristique d'un exemple d'un flux de travail combiné (voir, par ex., le flux de travail illustré dans la figure 2), est que l'utilisateur peut obtenir le maillage en 2,5D relativement rapidement et efficacement pour le comptage des cellules à chaque fois que les fractures le permettent sans besoin spécifique de savoir lesquels des fractures permettent ou pas une telle solution. Souvent, les fractures les plus conductrices sont représentées en utilisant une géométrie permettant le 2,5 D, et dans ces cas, les exemples de techniques présentées ici peuvent produire un maillage hautement raffiné dans ces régions, tout en ressoudant toujours toutes les autres fractures, bien qu'avec une résolution plus faible. [0032] Certains modes de réalisation peuvent maintenir les avantages du maillage à 2,5D, lorsque c'est possible, tout en permettant également l'utilisation des fractures complètement générales comme entrée. Ces exemples peuvent améliorer les algorithmes AGAR non- prolongés en éliminant le besoin que les éléments soient des hexaèdres. Certains exemples peuvent également améliorer les techniques qui dépendent seulement d'un maillage stadimétrique à 2,5D en éliminant le besoin que les fractures utilisées permettent un maillage à 2,5D. Tout en libérant les utilisateurs de ces restrictions, la prolongation AGAR divulguée et son association avec un maillage stadimétrique à 2,5D libèrent les utilisateurs du besoin de comprendre, et d'être impliqués de façon importante dans, le procédé de maillage. [0033] Des procédés illustratifs décrits ici peuvent modifier une grille structurée d'un modèle de la Terre pour ensuite créer un maillage quantitatif qui 1) peut être utilisé dans une simulation de réservoir, 2) respecte les propriétés et la géométrie du modèle de la Terre originel et 3) résout les systèmes de fracturation hydraulique naturelle additionnels. Les divers exemples décrits ici peuvent utiliser un algorithme AGAR prolongé pour insérer des réseaux de fracture dans un modèle de Terre existant, qui existe sous forme d'une grille structurée ayant des valeurs de propriétés physiques associées à celui-ci. Le réseau de fracture peut être sous la forme d'un maillage de surface (manifold ou non) dans un espace en 3D. La grille structurée du modèle de la Terre est ensuite raffiné de façon anisotrope pour résoudre les fractures et permettre un classement par grade approprié de la cellule au niveau de la région proche de la fracture. [0034] En règle générale, les cellules peuvent être examinées afin d'identifier celles qui ont besoin d'un raffinement. Ces cellules peuvent être raffinées dans les directions appropriées (ainsi, anisotrope). Lorsqu'aucune cellule n'a besoin de raffinement, on peut mettre fin au procédé. Le raffinement, tel qu'il est utilisé ici, décrit la division des cellules dans l'espace en 3D qui, en fin de compte, entraîne la formation d'un modèle de Terre à résolution de maillage plus élevé. Un réseau de fracture peut ensuite être résolu à l'intérieur du modèle de la Terre en utilisant les cellules raffinées. La résolution d'un réseau de fracture peut inclure la génération d'un modèle de la Terre et l'utilisation d'un maillage à haute résolution, dans lequel les fractures et les zones proches des fractures dans le modèle possèdent une résolution de maillage plus élevé que les autres zones du modèle. Étant donné que le maillage peut être utilisé pour une simulation numérique, en général, les cellules plus petites (après raffinement) peuvent donner des solutions plus précises en comparaison à des cellules plus grandes. Étant donné que les simulateurs de réservoir génèrent des approximations des formules mathématiques, plus les cellules sont petites, plus l'approximation mime de façon rapprochée les formules mathématiques. [0035] Étant donné que les procédés appliquent une approche anisotrope prolongée lors de la manipulation des cellules, le comptage des données et ses besoins informatiques associés peuvent être grandement réduits. Par conséquent, le procédé peut être plus rapide et plus efficace que les procédés conventionnels qui appliquent une analyse de force brute, nécessitant des calculs plus gourmands en ressources et en données. En outre, en raison de la relation du modèle de la Terre généré par rapport au modèle de la Terre originel, l'interpolation des caractéristiques est minimisée et le modèle de la Terre générée respecte également les valeurs des caractéristiques originelles. Ainsi, certains modes de réalisation de la présente divulgation permettent une modélisation avec une résolution de maillage plus élevée tout en réduisant les besoins en calculs informatiques. [0036] La FIG. 1 est un organigramme d'un système de simulation et de modélisation de réservoir 100. Le système 100 comprend un processeur 102, un stockage non-transitoire lisible par ordinateur 104, un module de communication émetteur/récepteur/réseau 105, des dispositifs entrée/sortie (1/0) optionnels 106 et un dispositif d'affichage optionnel 108 (par ex., un affichage pour une interface utilisateur), tous interconnectés à travers un système de bus 109. Les instructions logicielles qui peuvent être exécutées par un processeur 102 pour implémenter les instructions logicielles stockées à l'intérieur du moteur AGAR 110 peuvent être stockées dans un stockage 104 ou sur un quelconque support de stockage non-transitoire lisible par ordinateur. De la même façon, les instructions exécutables par un processeur 102 pour implémenter le module de maillage stadimétrique à 2,5D 116 peuvent être stockées dans un stockage 104. Comme illustrée dans la Fig. 1, le module de maillage stadimétrique à 2,5D 116 peut implémenter un algorithme de maillage en 2,5D ou être incarné sous forme d'un maillage stadimétrique en 2,5D. Même s'il n'est pas explicitement illustré dans la FIG. 1, le système 100 peut être connecté à un ou plusieurs réseaux publics ou privés à travers une ou plusieurs connexions de réseaux appropriées. Comme il est présenté en référence au système informatique 2900 illustré dans la FIG. 29, les instructions logicielles qui comprennent le moteur AGAR prolongé 110 et le module de maillage stadimétrique à 2,5D 116 peuvent également être téléchargées dans le stockage 104 à partir d'un CD-ROM ou d'un autre support de stockage approprié à travers des connexions câblées ou non. [0037] En outre, les techniques divulguées ici peuvent être réalisées avec une diversité de configuration des systèmes informatiques, comprenant des appareils portatifs, des systèmes multiprocesseurs, des composants électroniques grand public basés sur des multiprocesseurs ou programmables, des mini ordinateurs, des serveurs, des ordinateurs centraux, etc. Un nombre quelconque de systèmes informatiques et de réseaux informatiques sont acceptables pour une utilisation avec la présente divulgation. Les techniques divulguées ici peuvent être réalisées dans des environnements informatiques distribués dans lesquels les tâches sont réalisées par des dispositifs de traitement à distance qui sont reliés à travers des réseaux de communication. Par ex., les techniques divulguées ici peuvent être réalisées par un regroupement de dispositifs informatiques fonctionnant en groupe sur une ferme de serveurs.
Dans un environnement informatique distribué, les modules de programme peuvent être localisés à la fois sur un support de stockage informatique local ou à distance, y compris des dispositifs de stockage à mémoire. La présente divulgation peut donc être implémentée en relation avec divers matériels, logiciels ou une combinaison de ceux-ci dans un système informatique ou un autre système de traitement. [0038] Dans certains modes de réalisation, le moteur AGAR prolongé 110 comprend un module de système de fracture 112 et un module de modélisation de la terre 114. Le module du système de la fracture 112 procure la définition géométrique du réseau de fracture, avec les distributions de propriété pétrophysique et mécanique pour les fractures. Des exemples de modules de systèmes de fracture peuvent être implémentés en utilisant, par ex., les plates- formes Fracpro® ou Gohfer®. D'autres plates-formes et systèmes logiciels peuvent, cependant, être utilisés pour implémenter le module de systèmes de fracture 112. [0039] Le moteur AGAR prolongé 110 comprend également un modèle de modélisation de la terre 114. Le modèle de modélisation de la terre 114 s'intègre aux données de la fracture contenues dans le module du système de fracture 112 pour permettre une visualisation stratigraphique sous la surface des fractures et du réservoir. La visualisation peut comprendre, par ex., une interprétation géo-scientifique, une modélisation d'un système pétrolifère, une analyse géochimique, un maillage stratigraphique, un faciès, un volume net de la cellule et une modélisation de la propriété pétrophysique. En outre, chaque module de modélisation de la terre 114 modélise des trajets de puits, en sus d'une coupe à travers le faciès et les données sur la porosité. Les plates-formes de modélisation de la terre comprennent des plates-formes telles que, par ex., DecisionSpace®, qui est commercialement disponible chez Halliburton Energy Services Inc., Houston, Texas. Une diversité d'autres plates-formes de modélisation de la terre, peut, cependant, également être utilisée, avec ou à la place, avec la présente divulgation. [0040] Le moteur AGAR prolongé 110 peut également comprendre des circuits de traitement permettant une récupération robuste des données et l'intégration des données historiques ou en temps réel, apparentées au réservoir, qui couvrent tous les aspects de la planification des puits, les procédés de construction et d'achèvement tel que, par ex., le forage, le bétonnage, la diagraphie par câble, le test et la simulation de puits. En outre, de telles données peuvent comprendre, des données sur le forage, les trajectoires des puits, des données sur la propriété pétrophysique des roches, les données sur la propriété mécanique des roches, les données en surface, les données sur les failles, les données provenant des puits environnants, les données déduits des géostatistiques, etc. La base de données (non illustrée) qui stocke ces informations peut se trouver dans le module du système de fracture 112 ou à un emplacement distant. Une plate-forme de base de données est, par ex., le logiciel INSITE®, disponible commercialement chez Halliburton Energy Services Inc., Housto, Texas. Une diversité d'autres plates-formes de base de données, de plates-formes logicielles, et des systèmes associés peut être utilisée pour récupérer, stocker et intégrer les données apparentées aux puits, tel qu'il est décrit ici. [0041] En outre, un moteur AGAR prolongé 110 peut également comprendre des capacités d'automatisation du flux de travail multi-domaines, qui peut connecter une diversité d'applications techniques souhaitées. Ainsi, la sortie provenant d'une application, ou d'un module, peut devenir l'entrée pour une autre, procurant ainsi la capacité à analyser la façon dont divers changements affectent le placement du puits et/ou la conception de la fracture. Les gens de métier qui prennent connaissance de cette divulgation réalisent qu'il existe une diversité de plates-formes de flux de travail qui peut être utilisée dans cet objectif. [0042] Comme il est généralement décrit ci-dessus, les procédés et les modes de réalisation de la présente divulgation décrivent des algorithmes pour le raffinement automatique du maillage afin de résoudre des caractéristiques dans un réservoir pour des applications qui comprennent, sans limitation, la modélisation et la simulation de réservoir.
Un aspect important de la présente divulgation est la capacité d'appliquer sélectivement un raffinement anisotrope afin de réduire les degrés de liberté pour des analyses ultérieures utilisant le modèle de Terre généré. Exemple de flux de travail [0043] La FIG. 2 illustre un exemple d'un flux de travail 200 pour les créations des maillages de réservoir utilisant un raffinement anisotrope prolongé, adapté à la géométrie, d'un polyèdre. L'entrée 202 au flux de travail 200 comprend un modèle de réservoir (par ex., des données techniques du réservoir 204), et deux jeux de fractures 206 et 208. Un exemple de modèle d'entrée pour l'entrée 202 est donné dans la FIG. 8. Dans un mode de réalisation, le flux de travail 200 commence lors de la réception des données techniques d'un réservoir fournies par un utilisateur 204. Au niveau du bloc 214, un algorithme AGAR prolongé est utilisé en association avec un mailleur stadimétrique en 2,5D 210. Les données techniques du réservoir 204 est l'entrée nécessaire pour le mailleur stadimétrique en 2,5D 210. Comme illustré, le mailleur stadimétrique en 2,5D 210 peut être un mailleur par extrusion. Comme un composant additionnel de l'entrée 202 dans le flux de travail 200, un réseau de fracture fourni par l'utilisateur peut être reçu. Ce réseau de fracture peut être représenté sous forme de propriétés de géométrie et de fracture. Dans l'exemple du mode de réalisation de la FIG. 2, le jeu de fracture est divisé en 2 groupes : ceux qui doivent être résolus dans le maillage en 2,5D (par ex., les fractures permettant le 2,5D 206) ; et ceux qui ne doivent être résolus que dans le maillage en 3D final (par ex., d'autres fractures 208). Les caractéristiques des fractures permettant le 2,5D 206 sont ainsi appelées parce qu'elles permettent un maillage en 2,5D. Tel qu'illustrées, les caractéristiques des fractures permettant le 2,5D 206 peuvent comprendre des fractures représentées par des plans verticaux. Les caractéristiques des fractures permettant le 2,5D 206 peuvent également comprendre des fractures représentées par des plans sensiblement verticaux par rapport à l'horizon. Les autres fractures 208 peuvent comprendre des fractures représentées par des plans horizontaux. Les autres fractures 208 peuvent également comprendre des fractures représentées par des plans sensiblement horizontaux (par ex., des plans presque horizontaux). Les jeux de fractures 206 et 208 ne sont pas nécessairement exclusifs. C.-à-d., certains sous-ensembles de fractures peuvent se trouver dans les deux jeux de fractures 206 et 208. Dans le flux de travail 200, cependant, toutes les fractures d'entrée sont catégorisées comme étant dans au moins l'un des deux jeux de fractures 206, 208. [0044] Dans le flux de travail 200, le maillage en 2,5D est d'abord créé en utilisant un mailleur stadimétrique en 2,5D 210 afin de résoudre les fractures permettant le 2,5D 206. Le flux de travail 200 peut utiliser le mailleur stadimétrique en 2,5D 210 pour produire un maillage extrudé en 2,5 D, qui est illustré comme le modèle de réservoir intermédiaire 212 dans la FIG. 2. Le mailleur stadimétrique en 2,5D 210 peut générer le modèle de réservoir intermédiaire 212 de sorte qu'il comprenne un maillage d'extrusion qui modélise les fractures permettant le 2,5D 206 dans un espace en 3D. Le maillage d'extrusion modélise un certain volume du réservoir (par ex., le réservoir correspondant aux données techniques du réservoir 204), y compris les fractures permettant le 2,5D 206. Des exemples de mailleurs stadimétriques en 2,5D sont décrits en détail en référence aux FIG. 24-28 suivants. [0045] Ensuite, le flux de travail 200 applique l'algorithme AGAR prolongé au niveau du bloc 214 pour résoudre les autres fractures 208. Comme illustré, le modèle de réservoir intermédiaire 212 et le résultat de la résolution des autres fractures 208 sont associés pour donner un modèle de réservoir rempli 216. Le modèle de réservoir rempli 216 peut être compris dans la sortie 218 du flux de travail 200. Dans un exemple, la sortie 218 peut être présentée sous forme de graphique dans une interface utilisateur, tel que l'exemple d'interface donné dans la FIG. 9. Dans un exemple alternatif, si seulement les fractures permettant le 2,5D 206 sont comprises dans une entrée 202, le bloc 214 peut être sauté, et le modèle de réservoir intermédiaire 212 est utilisé comme le modèle de réservoir plein 216. Exemples de techniques de raffinement [0046] En général, les algorithmes de raffinement illustratifs décrits ici raffinent les rebords des cellules jusqu'à «n » fois et ne divise aucun rebord (et, ainsi, divise la cellule) qui produirait un rebord plus court que le `targetSize.' Tel que défini ici, le largetSize' représente la taille de maille souhaitée, ou la longueur de rebord, permettant de résoudre la largeur de la fracture. Dans certains modes de réalisation, la valeur du targetSize est fournie au système sur la base d'un niveau souhaité de précision et le délai souhaité jusqu'à la solution. Afin de générer les modèles de Terre, les procédés peuvent déterminés si une cellule doit être raffinée et dans quelle direction ces cellules doivent être raffinées. Dans certains modes de réalisation, il peut y avoir deux types de règles permettant de déterminer si un cellule doit être raffinée : 1) les règles de gradation (les règles qui assurent des transitions lentes du côté élément) et 2) les règles de croisement (les règles qui assurent que les fractures sont adéquatement représentées). [0047] Dans un exemple pour résumer le procédé, considérons les directions U (une direction dans un espace 3D). Il est à noter qu'U est une direction dans un sens topologique et non dans un sens Euclidien. La "direction U" pour chaque cellule est indépendante de la direction U pour une cellule avoisinante. En outre, dans une cellule, un rebord 'Il' peut pointer dans une direction légèrement différente qu'un autre rebord 'Il'. Comme il sera décrit plus en détail ci-dessous, en considérant la direction U, le moteur AGAR prolongé 110 analysera toutes les cellules dans un modèle et raffinera des cellules presque fracturées dans la direction U si tous les rebords U sont plus longs que C x targetSize et qu'au moins une des règles suivantes est satisfaite : 1) un rebord U échelonné croise une fracture ; 2) pour toute paire de rebords opposée dans la direction U, exactement l'un des deux rebords échelonnés croise une fracture ; 3) un rebord U comporte deux ou plusieurs noeuds 'suspendus' ; 4) pour toute paire de rebords opposée dans la direction U, exactement l'un des deux rebords échelonnés comportent deux ou plusieurs noeuds suspendus ; ou 5) la cellule non-échelonnée croise une fracture mais aucun rebord échelonné de la cellule ne croise une fracture. En se référant au C x targetSize, pour un variable scalaire spécifique à une implémentation ou fourni par un utilisateur, C, le système 100 ne divisera aucun rebord (donc, des cellules) qui est plus court que C x targetSize. Dans ce procédé illustratif, C est choisi comme -\/2, mais peut être d'autres variables tels que, par exemple, 4/3 ou toute valeur finie. Ce procédé est répété dans chaque direction dans l'espace 3D. [0048] Tel qu'il est utilisé ici, le terme "noeud suspendu" est généralement utilisé pour décrire un noeud créé au cours du raffinement d'une cellule avoisinante, qui n'est pas nécessaire pour maintenir la géométrie sous-jacente de la cellule d'intérêt. Par exemple, la géométrie sous-jacente peut être un hexaèdre (tel que dans les exemples des FIG. 3, 8 et 12- 13), un prisme (tel que le démontre les exemples des FIG. 4-8 et 10-11), ou d'autres types d'éléments extrudés. Considérons une cellule H, lorsqu'une cellule avoisinante est raffinée, un noeud supplémentaire est ajouté à la cellule H en forme d'hexaèdre. Ce noeud est "suspendu" par rapport à la cellule H. Dans encore d'autres exemples, l'exigence de la longueur du rebord peut être omise pour les règles 3 et 4 si une adhérence stricte à la directive du noeud suspendu est adoptée. [0049] La FIG. 3 illustre un jeu 300 de configurations de raffinement possibles 306, 312 et 316 pour une cellule hexaèdre unique 304 à l'intérieur d'un espace en 3D U-V-W 302. En particulier, la FIG. 3 illustre comment, pour un hexaèdre, certains modes de réalisation peuvent utiliser un espace local, logique U-V-W 302, et chaque cellule est analysée et raffinée dans un certain sous-ensemble de ces trois dimensions, U, V et W comme le démontre les étapes 308, 310 et 314 de la FIG. 3. [0050] Une fois que les directions de raffinement ont été déterminées pour une cellule, un moteur AGAR prolongé 110 réalise le raffinement (c'est à dire, diviser les bords) pour cette cellule. La FIG. 3 illustre huit configurations de raffinements possibles pour une cellule unique à l'intérieur d'un espace en 3D U-V-W 302. Le moteur AGAR prolongé 110, via l'algorithme, peut réaliser un raffinement de façon itérative en déterminant d'abord si les cellules proches de la fracture doivent être raffinées dans les dimensions/directions U, V ou W (étapes, 308, 310, et 314 illustrées dans la FIG. 3), et ensuite, si tel est le cas, la division des rebords dans les directions illustrées pour créer des raffinements de cellules 306, 312 et/ou 316. Les huit configurations possibles de raffinement de cellules 304, 306, 312 et 314 vont du raffinement sans direction de raffinement (hexaèdre 304) au raffinement dans toutes les 3 directions de raffinement (le raffinement le plus à droite 316). Étant donné que ces cellules seront celles qui seront positionnées dans la région proche de la fracture, la résolution du maillage des images ainsi obtenues sera grandement améliorée. [0051] La FIG. 4 illustre une cellule de prisme triangulaire représentative 404 et un système de coordonnée logique à 2D U-W 402 et un schéma de numéro pour les rebords du prisme 404. Certains modes de réalisation comprennent deux alternatives pour étendre le concept de raffinement susmentionné aux prismes. L'une consiste en l'option « direction U unique » illustrée dans la FIG. 4. Comme illustrée dans la Fig. 4, la direction U unique englobe tous les rebords de face triangulaire du prisme 404. Selon cette définition des directions, le raffinement, lorsqu'il se produit sur les faces triangulaires du prisme 404, est toujours isotrope par rapport à la face. Dans l'exemple de la FIG. 4, en utilisant l'option direction U unique 400, chacun des rebords 1 à 6 est considéré comme étant dans la direction U. Le raffinement dans U raffinera tous les rebords et le raffinement peut se faire dans l'un quelconque des sous-ensembles de {U, W} . Une autre option est appelée ici l'option « direction U multiples ». Celle-ci est illustrée dans la FIG. 5. La FIG. 5 illustre une cellule de prisme triangulaire 504 et un système de coordonnée logique à 2D U-W 502 et un schéma de numéro pour les rebords du prisme 504. En utilisant l'option « direction U multiples » 500 de la FIG. 5, le raffinement serait effectué dans les paires de direction U. Ainsi, le raffinement pourrait être réalisé dans un sous ensemble de directions 1U1U2, U1U3, U2U3, W} 502 illustré dans la FIG. 5. En règle générale, dans les règles appliquées par certains modes de réalisation, on n'autorise pas toutes les 16 permutations dans les 4 directions (par ex., W, Ul, U2, et U3) au cours de l'étape de raffinement. Dans l'exemple de la FIG. 5, les multiples directions U sont considérées en paires. Par ex., le raffinement dans U1U2 raffinera les rebords 1, 2, 4, 5 du prisme 504. [0052] La FIG. 6 illustre un exemple du paradigme à direction unique U 600. Suite à la notation générale et à l'algorithme prolongé AGAR présenté ci-dessus, et l'utilisation de la numérotation canonique du prisme triangulaire dans la FIG. 6, étiqueter les rebords de la cellule, Ci, comme teid) pour J = 1...9. Considérons maintenant une fonction Booléenne abstraite sur les rebords, Test(( ). Par ex., la cellule peut être raffinée dans : -[6.,,1 , ei,2 ,ei.3, et. e r , U si un quelconque rebord provenant de l'ensemble satisfait Test(teid}) (voir l'étape 606 de la FIG. 6), et -. W si un quelconque rebord provenant de l'ensemblete e 1,7' satisfait Test(teJ) (voir étape 610 de la FIG. 6). [0053] La condition susmentionnée est appelée la « condition principale » pour la gradation. Nous allons maintenant définir les « conditions opposantes ». Pour les quadrilatères dans le prisme 604, la règle est la même pour les hexaèdres (par ex., comme il a été présenté ci-dessus en référence à la FIG. 3). Si, disons deux rebords opposés sur la face, exactement un de ces rebords satisfait Test(tei/}), alors la direction qui divise ces 2 rebords seront également raffinés. Par exemple, les rebords tei,4 ,11 sont des rebords opposés sur une face de quadrilatère. Si le rebord tei,41 satisfait Test(tei,i}) mais pas le rebord teial, alors cette condition sera satisfaite, et la direction W serait ajoutée à l'ensemble de raffinement. De la même façon, si le rebord satisfait satisfait Test(te1,/}) mais pas le rebord [e.) la la direction W serait ajoutée au jeu de raffinement. En ce qui concerne la face triangulaire, si exactement un satisfait Test(tei,/}), alors U est ajouté au jeu de raffinement. Par exemple, si le rebord tei,4) t te e satisfait Test( id ) mais pas les rebords i,5' i,6), la direction W serait ajoutée au jeu de raffinement. [0054] En ce qui concerne les règles de croisement de l'algorithme AGAR prolongé, nous définissons Test(teid}) comme vrai si et seulement si te } a été raffiné au moins 2 fois. En ce qui concerne les règles de gradation de l'algorithme AGAR, nous définissons Test(teid)) comme vrai si le rebord échelonné croise une fracture d'intérêt. La FIG. 6 illustre les combinaisons possibles de raffinements 608 et 612 utilisant le paradigme à direction U unique 600. [0055] La FIG. 7 illustre la façon dont on peut, dans une alternative de la direction U unique présentée ci-dessus en référence à la FIG. 6, sélectionner un paradigme à multiples directions U 700. En particulier, la FIG. 7 illustre les raffinements 708, 712, 716, 720 et 722 dans les 4 dimensions/directions logiques 702 pour un paradigme à multiples directions U700 d'un prisme 704. En utilisant le paradigme à multiples directions U 700, la cellule sera raffinée dans - U1U2 et U1U3 si un quelconque rebord provenant de l'ensemble tei,1 satisfait Test(teEd}) (vers les étapes 706 et 710 de la FIG. 7), - te ,r satisfait U1U2 et U2U3 si un quelconque rebord de l'ensemble E,2 Test(t i,./}) (voir les étapes 706 et 714 de la FIG. 7), - U1U3 et U2U3 si un quelconque rebord de l'ensemble te e i,3 satisfait 15te 1 ei, satisfait Test(tei,/}) voir Test( td- ) (voir les étapes 710 et 714 de la FIG. 7), et -te W si un quelconque rebord de l'ensemble l'étape 718 de la FIG. 7). [0056] En sus de ces conditions principales, nous avons des conditions opposées. Pour les quadrilatères, les deux paradigmes U sont très semblables. Si, prenant deux rebords 20 opposés sur la face, exactement l'un de ces rebords satisfait Test(teid)), alors les directions qui divise les deux rebords seront également raffinées. Si le rebord te[,7) satisfait Test(te 1) mais pas le rebord tet,81, la direction W serait ajoutée au jeu de raffinement. Dans un exemple, pour la règle de la face triangulaire, si exactement un rebord satisfait Test(t JI), alors U1U2 et U1U3 peuvent être ajoutés au jeu de raffinement. Dans ce paradigme, nous n'ajoutons aucune 25 condition opposée pour les faces triangulaires. Cependant, des règles alternatives peuvent être déterminées et testées dans le cadre des paradigmes au-delà des paradigmes 600 et 700 illustrés dans les FIG. 6 et 7, respectivement. [0057] Les Figures 8 et 9 sont des exemples d'interfaces utilisateur. En particulier, la FIG. 8 illustre une interface utilisateur 800 affichant un modèle d'entrée contenant à la fois 30 des prismes hexaèdres 802 et triangulaires 804. La FIG. 9 illustre une interface utilisateur affichant un jeu de fracture 902 raffiné dans le modèle de la FIG. 8 utilisant le paradigme à direction U unique. Comme illustrée dans la Fig. 9, les avantages anisotropes sont grandement réduits pour les prismes, utilisant le paradigme à direction U unique, mais seulement dans les directions horizontales. La majeure partie du comportement anisotrope est maintenue dans la direction verticale. La raison en est que la direction verticale n'est pas nécessairement raffinée. [0058] La FIG. 10 illustre un prisme non-triangulaire représentatif 1004 (par ex., un élément prismatique général) montrant le système coordonné logique en 2D 1002. La FIG. 11 illustre un raffinement possible d'un élément prismatique général 1104 en subdivisant l'élément en hexaèdres. Comme il est décrit ci-dessus, un algorithme AGAR prolongé peut être utilisé pour implémenter des procédés qui prolongent le schéma AGAR basé sur un hexaèdre pour également inclure des prismes triangulaires, tels que ceux illustrés dans les FIG. 4-7. Pour les maillages quadrilatères et triangulaires prolongés, cette combinaison est suffisante. En assumant une topologie d'extrusion cohérente, les paradigmes présentés ci-dessus en référence aux FIG. 4-7 entraîneront la formation d'un maillage valable même si certains des éléments sont dégénérés. Les Figures 10 et 11 illustrent la façon dont certains modes de réalisation peuvent directement prolonger le paradigme à direction U unique vers des prismes non triangulaires, utilisant un algorithme de subdivision mi-point. Pour les prismes non-triangulaires, tels que le prisme non-triangulaire 1004 de la FIG. 10, on peut maintenir l'anisotropie dans l'espace U-W 1002. Comme on peut le voir dans le raffinement 1104 de la FIG. 11, en utilisant cet algorithme de raffinement, les éléments prismatiques non- hexaèdres 1004 sont subdivisés en hexaèdres. Étant donné que les hexaèdres ont 3 dimensions logiques, comparé aux deux illustrées dans les FIG. 10 et 11 pour les éléments prismatiques non-hexaèdres, l'opportunité pour un raffinement anisotrope est augmentée après le raffinement de l'élément originel. [0059] Pour généraliser davantage, nous pouvons appliquer un algorithme de subdivision mi-point pour toutes les cellules polyèdres convexes, même les cellules non-extrudées. Il est à noter que ceci modifiera l'algorithme de raffinement pour les prismes triangulaires qui sont présentés ci-dessus en référence aux FIG. 4-7. Cependant, pour les cellules non-extrudées, nous utilisons un espace U à dimension unique, écartant les raffinements anisotropes pour ces cellules générales. Dans les cas où la majorité des cellules sont extrudées (cellules hexaèdres ou d'autres cellules prismatiques), l'anisotropie pour ces cellules procurera toujours un avantage important par rapport au raffinement totalement isotrope. En outre, le raffinement anisotrope sur une cellule régulière avoisinant une cellule irrégulière, ne viole pas l'exigence de conformité sur le maillage. La face polygonale de la cellule avoisinante irrégulière est tout simplement divisée. [0060] Maintenant qu'une description sommaire d'un exemple de procédé de raffinement a été donnée, une discussion plus détaillée d'exemples de techniques pour la création de maillage de réservoir sera maintenant présentée. Afin de générer les modèles de réservoir de la Terre, un exemple de procédé détermine d'abord quelles cellules doivent être raffinées et dans quelle direction ces cellules doivent l'être. Pour réaliser ceci dans cet exemple, un moteur AGAR prolongé 110 commence d'abord avec une grille hexaèdre structurée, uniforme et axialement alignée. Cependant, dans d'autres exemples, il n'est pas nécessaire de commencer avec une grille uniforme, axialement alignée. Néanmoins, lorsque des données sont entrées dans l'algorithme de raffinement, un moteur AGAR prolongé 110 nécessite : 1) M = Cd, le jeu initial des cellules de maillage devant être raffiné, 2) F = {FI}, le jeu de facettes triangulaires représentant les fractures, 3) co, la « largeur » souhaitée des cellules pour résoudre adéquatement les fractures, et 4) a, la taille du maillage souhaité « le long des » fractures. [0061] Parce que la grille initiale, M, est uniforme, il existe une taille de grille caractéristique, TT. Ainsi, le moteur AGAR prolongé 110 calcul deux nombres : niso = logi Éq.(1), et T 2 Éq.(2), où ntot = g1- -T 2 rani= ntot - niso. qui définira le nombre de niveaux de raffinement réalisé par le moteur AGAR prolongé 110. Le moteur AGAR prolongé 110 effectuera nn veaux Iso n.sode raffinement isotrope pour obtenir ÎX , proche des fractures. Le moteur AGAR prolongé 110 un maillage de taille, effectuera nn veaux an n an/de raffinement isotrope pour obtenir un maillage de « largeur », ntot 2 , proche des fractures. Il est à noter que nous avons a=-'--1 a et [0062] Par conséquent, dans un procédé illustratif, l'algorithme général appliqué à un moteur AGAR prolongé 110 peut être exprimé sous la forme de : raffiner(n, doIsotropic, targetSize) co =25 pour itération = 1... ,-n-min(n,itération) sizeLimit = targetSize * z pour C in M //boucle au-dessus de chaque cellule dans le maillage pour cette itération directions = determineDirections(C,F, iteration,sizeLimit,n,doIsotropic) refineInGivenDirections(directions,C) si aucun changement n'a été effectué durant cette itération routine de sortie [0063] Dans cette routine 1, l'algorithme « Raffiner » est appelé avec un nombre entier pour indiquer le nombre de niveaux de raffinement à exécuter et un Booléen pour indiquer si c'est un raffinement isotrope ou non. Il est à noter que, dans cet exemple, la routine est appelée deux fois. Le moteur AGAR prolongé 110 appelle raffiner(niso,true,â) suivi de nam, f aise, Essentiellement, pour exécuter n niveaux de raffinement, le moteur AGAR prolongé 110 fait une boucle au-dessus des cellules jusqu'à ce qu'aucun autre changement n'est nécessaire. Étant donné que le moteur AGAR prolongé 110, via l'algorithme, interdit la division de tout rebord qui est plus court que C x targetSize au cours d'une itération, en fin de compte, l'algorithme se terminera. L'une des parties les plus complexes de l'algorithme est encapsulée dans le procédé « determineDirections », donné dans la routine 2 ci-dessous qui retourne un sous-ensemble de direction locale en 3D {U, V, W} dépendamment des directions qui doivent être raffinées. Il est à noter que « sizeLimit », utilisé pour définir le rebord le plus petit lors d'une itération, est grand pour les itérations précoces. Par le l'itération, et pour toutes les itérations subséquentes, sizeLimit est égal à targetSize. La Routine 2 peut être exprimée sous forme de : determineDirections(C,F,iteration,sizeLimit,n,doIsotropic) directions = 0 // Ajouter une quelconque direction de raffinement nécessaire à cause du croisement de la fracture/cellule //Ceci n'est permis que pendant les n premières itérations si l'itération <= n // dans le stade de croisement directions = getDirectionsFromIntersections(C,F) //ajouter une quelconque direction de raffinement nécessaire en raison de la gradation //ceci ne sera nécessaire que pour l'itération > 1 Si l'itération > 1 directions = directions U getDirectionsFromGradation(C)C) // si nous voulons isotrope et nous trouvons une quelconque direction, effectuer toutes les directions si doisotrope & (directions 0) directions = ALL // supprimer toutes les directions interdites, indépendamment de ce qu'on a trouvé directions = directions - getDisallowedDirections(C,F,sizeLimit) directions de retour [0064] Dans cette Routine 2 illustrative, `determineDirections' utilise principalement des sous-routines getDirectionsFromGradation et getDirectionsFRomIntersections. Comme le démontre le `determineDirections' à partir de la Routine 2, c'est seulement au cours des premières n itérations, que le moteur AGAR prolongé 110 raffine les cellules proches de la fracture en fonction de l'interaction avec les fractures, qui sera décrit plus en détail ci-dessous. Au cours d'une première itération après que le moteur AGAR prolongé 110 raffine les cellules en se basant sur la 'règle d'un noeud suspendu'. C'est à dire, si un rebord de la cellule a été divisé plus d'une fois (en raison d'un voisin qui est raffiné), la cellule peut être raffinée. Si le moteur AGAR prolongé 110 réalise un raffinement "isotrope", alors, à chaque fois qu'il a besoin de raffiner une cellule dans une quelconque direction, il raffine la cellule dans toutes les directions. [0065] Toutes les directions identifiées pour le raffinement doivent passer un test final avant d'être acceptées. En faisant le raffinement dans une direction donnée, 4 rebords peuvent être divisés. Si l'un quelconque de ces 4 rebords est plus court que C x targetSize, alors la direction est rejetée par le moteur AGAR prolongé 110. [0066] Étant donné que le maillage initial peut comprendre un hexaèdre, chaque cellule peut être représentée par un système de numérotation canonique décrit ci-dessous en référence à la FIG. 12. La FIG. 12 illustre une cellule représentative 1204 montrant le système de coordonnée en 3D local 1202 et le schéma de numérotation 1200. Dans la FIG. 12, les directions U-V-W sont définies par le mappage des noeuds sur le système de numérotation (à noter qu'il est dessiné comme un système gaucher). Dans le schéma de numérotation 1200, les numéros représentant les rebords sont encerclés, alors que les numéros représentant les noeuds ne le sont pas. Le mappage du système U-V-W 1202 sur le système x-y-z peut différer des cellules avoisinantes. Dans cet exemple, prenons les rebords {0,2,4,6} comme les rebords U, les rebords {1,3,5,7} comme les rebords V, et les rebords {8,9,10,11} comme les rebords W. [0067] En utilisant le schéma de numérisation canonique dans la FIG. 12, les rebords de la cellule, Ci, sont étiquetés comme { ev} pour j = 0 ... 11. Le jeu de retour provenant de getDirectionsFromGradation(C1) comprend : 1) U si un quelconque rebord de l'ensemble {e,0, e1,2, e1,4, e1,6} a été raffiné au moins 2 fois, 2) V si un quelconque rebord de l'ensemble { e'i, Co, e,5, a été raffiné au moins 2 fois, et 3) W si un quelconque rebord de l'ensemble { e,8, e1,9, e1,10, e1,11} a été raffiné au moins 2 fois. Cette condition est appelée la « condition principale » pour la gradation. En outre, si nous prenons deux rebords opposés sur une face, exactement un de ces rebords a été raffiné 2 fois, alors la direction qui divise ces 2 rebords seront également raffinés par le moteur AGAR prolongé 110. Par exemple, les rebords { e,8, Co} sont des rebords opposés sur une face. Si le rebord {e,8} a été raffiné 2 fois mais le rebord { e,9} a été raffiné une fois, alors cette condition serait satisfaite, et la direction U serait ajouté à l'ensemble de raffinement. Ce test est appelé la « condition opposée » pour la gradation. À noter que pour tous les tests dans getDirectionsFromGradation, le moteur AGAR prolongé 110 considère seulement le nombre de fois que le rebord a été divisé (c.-à-d., raffiné) avant l'itération actuelle. De cette façon, le moteur AGAR prolongé 110 diminue sa dépendance sur l'ordre dans lequel les cellules peuvent être visualisées dans l'algorithme.
Test du croisement cellule-facette [0068] Un aspect de l'algorithme de raffinement utilisé dans certains modes de réalisation est la définition d'un test pour déterminer si une cellule donnée, Ci, croise une facette de fracture donnée, Fk. Voir les FIG. 4, 5 et 12 par ex., les schémas de numérotation utilisée par le moteur AGAR prolongé 110 pour les noeuds et les rebords de la cellule. Afin de continuer de décrire ce procédé illustratif, le moteur AGAR prolongé 110 peut définir la transformée d'échelonnement, Sfl(), qui échelonne la distance entre le centre d'un objet (le noeud) et chaque noeud par la valeur scalaire, fi. Dans cet algorithme illustratif, fi est défini comme 2, mais la détermination de la valeur la plus appropriée de fi est un choix de conception et peut être modifié, comme le comprendra un homme de métier qui a connaissance de cette divulgation. À noter que le moteur AGAR prolongé 110 peut appliquer cette transformée à la cellule, produisant une cellule échelonnée, ou aux rebords, donnant un jeu de rebords échelonnés. Ces deux applications donnent des résultats différents, comme le démontre les FIG. 13 et 14, qui illustre la fonction échelonnée, Sfl(), appliquée à une cellule, et la fonction d'échelonnage, Sfl(), appliquée aux rebords de la cellule, respectivement. [0069] Toujours en référence aux FIG. 13 et 14, la fonction Boolénne Tfl(C'F) est vraie si et seulement s'il existe une fracture Fk E F de sorte que le croisement de la cellule échelonnée (Se(C3 n Fk 0). et de la fracture n'est pas vide Pour les rebords, la fonction est légèrement altérée de sorte qu'elle sera fausse à chaque fois que le rebord se trouve dans le plan de la facette. Ainsi, pour les rebords, Tfl(e,J,F) si et seulement s'il existe une fracture Fk E F de sorte qu'à la fois le croisement du rebord échelonné et de la fracture n'est pas vide la normale de Fk n'est pas perpendiculaire (ou transversal dans d'autres exemples) au rebord (R.S'e(e,) n Fk # 0)&(normal(F k) e1> E) Essentiellement, la fonction est vraie lorsque l'entité échelonnée croise les facettes, excepté dans le cas où le rebord se trouve dans le plan de la facette. Ainsi, ceci entraîne la division du rebord le long d'un plan qui est transversal (par ex. perpendiculaire) à un axe de la fracture à l'intérieur du réseau de fracture. [0070] Maintenant, un algorithme illustratif utilisé dans la Routine 2 qui définit getDirectionsFromIntersections(C'F) est spécifié. À noter que, comme dans le cas des tests de gradation ci-dessus, les tests de croisements comprennent des « conditions principales » et des « conditions opposées » analogues. Cette routine 3 peut être exprimée comme, par ex. : Set<Directions> getDirectionsFromIntersections(C'F) directionSet = // ne rien faire si aucune facette ne croise la cellule échelonnée Si Tfl(C'F) == faux retourner directionSet // la condition « principale » //Si un rebord échelonné croise le jeu de fracture, nous effectuons un raffinement perpendiculairement à celui-ci T j3(e F)IT F)IT 13(e 0, F)IT directionSet = directionSet U U (e 0, F)1T (e , F)IT i3(e i,s, F)IT i,7 directionSet = directionSet U V pour un certain c petit et spécifié). si (T(e48, F) I Tp,(eo, F) I Tp,(e00, F) I T /3(e directionSet = directionSet U W //la condition « opposée » //Si les rebords « opposés » ont ainsi différents, effectuer un raffinement de la direction de division 'T (e8, F) # T f3( F) T e(e 00, F) # T f3(e F) si T 13(e F) T 15,(e F) I T 13(e i,5, F) # T 13(e, 7, F) directionSet = directionSet U F) # T e(e F) I T /3(e F) # T e(e 00, F) si T 13(e F) # T #(e 42, F) I T i3(e F) # T 13(e 46, F) directionSet = directionSet U V ,a(e 43, F) # T e(e 47, F) I T (e F) # T )5,(e F) si T (e 40, F) # T /3(e n I T /3(e i,2, F) # T (e 46, n directionSet = directionSet U W //Si le jeu est vide mais qu'une facette croise la cellule non-échelonnée, raffiner tout Si(directionSet , 0 & T0 (C, F)) directionSet =UUVUW [0071] À noter que l'algorithme « Raffiner » peut être invoqué ou appelé par un moteur AGAR prolongé 110 avec un nombre entier pour indiquer le nombre de niveaux de raffinement à exécuter et un Bouléen pour indiquer si c'est un raffinement isotrope ou non. [0072] Les Figures 15 et 16 illustrent les images 1500 des maillages AGAR prolongés générés avec des caractéristiques illustratives de la présente divulgation. Dans certains exemples de procédés, les « conditions opposées » permettant de déterminer les directions du raffinement, peuvent être utilisées pour restreindre le raffinement anisotrope dans une région localisée (par ex., proche de la fracture). Par ex., prenons la FIG. 15, qui démontre un raffinement AGAR proche de l'extrémité d'une fracture plane 1500, et la FIG. 16, qui montre le même modèle de la FIG. 15 utilisant un prolongement de l'AGAR, les « conditions opposées » étant désactivées. Dans les exemples des FIG. 15 et 16, la fracture 1500 est illustrée comme étant entourée de diverses cellules 1502 dans la région proche de la fracture. En absence des « conditions opposées » pour la sélection de la direction du raffinement tel que le démontre la FIG. 16, l'extrémité de la fracture 1500 n'aurait pas été résolue et le raffinement anisotrope aurait continué à travers le haut de l'image 1500. [0073] Les Figures 17 et 18 sont des images additionnelles des maillages AGAR générées utilisant des aspects illustratifs de la présente divulgation. La FIG. 17 illustre une extrémité de la fracture 1500 avec fi = 2, donnant ainsi un maillage contenant 1500 cellules. La FIG. 18 illustre une extrémité de la fracture 1500 avec fi = 4, donnant ainsi un maillage avec 2260 cellules. [0074] Les Figures 19 et 20 montrent un modèle de réservoir rempli avec à la fois des maillages de fractures hydrauliques et naturelles utilisant les procédés illustratifs décrits ici. Dans la FIG. 19, une tranche à travers un exemple de maillage généré avec l'algorithme de raffinement est illustrée, comprenant la fracture hydraulique 1900 et les fractures naturelles 1902. Dans l'exemple de la FIG. 19, la taille de maillage caractéristique du maillage originel était de 250 pieds et les tailles cibles étaient w = 3,125 pieds, a = 250 pieds. C.-à-d., pour ce maillage, aucun des niveaux n'utilise le raffinement isotrope, qui donne flani = 6. Comme illustrée dans la Fig. 19, la résolution du maillage permise par les cellules proches de la fracture est plus élevée dans les zones entourant les fractures 1900 et 1902 en comparaison aux zones qui n'entourent pas les fractures 1900 ou 1902. Dans la FIG. 20, le modèle de la FIG. 19 est illustré avec certaines fractures 2000, 2002 superposées. Exemples de procédés [0075] À la lumière de la description détaillée précédente, un procédé illustratif sera maintenant décrit en référence aux organigrammes des FIG. 21 et 22. Les procédés 2100 et 2200 illustrés dans les FIG. 21 et 22, respectivement, sont décrits en référence aux exemples de caractéristiques de la FIG. 1. Cependant, les procédés 2100 et 2200 ne sont pas limités à ces exemples de caractéristiques. Comme illustré dans la Fig. 21, le procédé 2100 commence au niveau du bloc 2102. Dans un exemple, le bloc 2102 comprend un système d'initialisation 100 afin de générer un modèle de réservoir de Terre. Au cours de l'initialisation, le moteur AGAR prolongé 110 télécharge des données à partir du module du système de fracture 112 et du module de modélisation de la Terre 114. De telles données peuvent comprendre, par ex., une grille géocellulaire et des données de fracture. Au niveau du bloc 2104, le moteur AGAR prolongé 110 analyse toutes les cellules dans un modèle de Terre en 3D afin de déterminer si les cellules dans un maillage doivent être raffinées de façon anisotrope (par ex., raffinées dans une quelconque direction dans l'espace 3D). Si un moteur AGAR prolongé 110 détermine que la cellule qui est analysée ne nécessite pas un raffinement, la commande est transférée au bloc 2106. [0076] Au niveau du 2106, le moteur AGAR prolongé 110 détermine ensuite s'il existe d'autres cellules à analyser. S'il est déterminé qu'il n'existe pas d'autres cellules à analyser, la commande est transférée au bloc 2116 et le procédé 2100 se termine. Sinon, s'il est déterminé qu'il existe une autre cellule qui reste à analyser, la commande est retransmise au bloc 2104 où cette cellule sera maintenant analysée pour déterminer si le raffinement pourrait être nécessaire. [0077] Si, au niveau du bloc 2104, le moteur AGAR prolongé 110 détermine que cette cellule nécessite un raffinement, l'algorithme repasse au niveau du bloc 2108 dans lequel cette cellule est raffinée. Comme il a été précédemment décrit, les règles de gradation et de croisement peuvent être utilisées pour déterminer si une cellule doit être raffinée. Une fois qu'un moteur AGAR prolongé 110 raffine la cellule dans les directions déterminées comme étant nécessaires, l'algorithme passe au niveau du bloc 2109 où une détermination est réalisée afin de voir s'il y a d'autres cellules à analyser. S'il est déterminé qu'il existe d'autres cellules à analyser, la commande est retransmise vers le bloc 2106. Sinon, s'il est déterminé qu'il n'existe pas d'autres cellules à analyser, la commande est retransmise vers le bloc 2110. Au niveau du bloc 2110, le moteur AGAR prolongé 110 est utilisé pour résoudre la fracture à l'intérieur du modèle Terre en utilisant les cellules raffinées. Le modèle ainsi obtenu possédera un niveau plus élevé de résolution de maillage dans ces zones entourant les fractures (par ex., zones proches de la fracture et les fractures) et une résolution de maillage moins élevée dans ces zones proches de la fracture. [0078] La FIG. 22 est un organigramme d'un procédé de raffinement illustratif 2200 utilisé dans le bloc 2104 présenté ci-dessus en référence à la FIG. 21. Pour commencer cette analyse, il faut noter que l'emplacement des fractures à l'intérieur du système de fracture est connu à partir des données provenant du module du système de fracture 112. En utilisant ces données, aussi bien qu'avec les règles de gradation et le croisement précédemment décrites, le moteur AGAR prolongé 110 peut localiser ces cellules dans les régions proches de la fracture du modèle qui nécessite un raffinement. Pour les cellules qui nécessitent un raffinement, tous les rebords (et, ainsi, les cellules) seront divisés en deux dans chaque direction de raffinage. Les cellules nouvellement créées seront ensuite de nouveau analysées dans la prochaine itération, augmentant, en fin de compte, la résolution du maillage. Le maillage à résolution élevée ainsi obtenue procure des résultats numériques plus précis dans les régions proches de la fracture et les régions de la fracture. La C x targetSize (où C = -\/2, par ex.) est déterminée comme précédemment décrit ci-dessus dans la Section A. [0079] Néanmoins, au niveau du bloc 2104(i), le moteur AGAR prolongé 110 détermine si tous les rebords dans les directions U-V-W pour la cellule (chaque direction est analysée séparément) sont plus longs que C x targetSize. Si ces rebords sont plus longs que C x targetSize, aucun raffinement n'est requis, et la commande est transmise au bloc 2106. Si, cependant, tous les rebords sont plus longs que C x targetSize, les 5 règles peuvent ensuite être appliquées. À ce point, au niveau des blocs 2104(ii)-(vi), le moteur AGAR prolongé 110 détermine si un ou plusieurs des éléments suivants sont satisfaits : ii) un rebord directionnel échelonné croise une fracture ; iii) pour une quelconque paire de rebords opposés dans la direction 3D, exactement un des deux rebords échelonnés croise une fracture ; iv) un rebord directionnel possède deux ou plusieurs noeuds « suspendus » ; v) pour une quelconque paire de rebords opposés dans la direction 3D, exactement un des deux rebords possède deux ou plusieurs noeuds suspendus ; et vi) la cellule non-échelonnée croise une fracture mais aucun rebord échelonné de la cellule ne croise une fracture. [0080] Au niveau du bloc 2104(h), en référence à la FIG. 13, le moteur AGAR prolongé 110 détermine si la version échelonnée de la cellule 1300 (c.-à-d., 1300S) possèdent un rebord dans la direction analysée qui croise une fracture. Au niveau du bloc 2104(iii), le moteur AGAR prolongé 110 (par exemple) peut analyser une paire de rebords dans une direction donnée (par ex., les rebords 5 et 1 dans la FIG. 12), afin de déterminer si seulement un des rebords échelonnés croise la fracture. Au niveau du bloc 2104(iv), le moteur AGAR déterminera si un rebord dans la direction analysée possède 2 ou plusieurs noeuds suspendus. Dans la FIG. 22, l'élément 2209 est donné comme une illustration de ce principe. L'élément 2209 est décrit ci-dessous. [0081] Comme illustré dans la Fig. 22, l'élément 2209 comprend deux faces de cellules adjacentes 2210 et 2220 (illustrées en 2D pour des raisons de simplicité) ayant deux noeuds 2212 et 2214 positionnés entre eux. Dans l'exemple de l'élément 2209, le premier raffinement aurait divisé la cellule 2210 en deux en utilisant le noeud 2212. Le deuxième raffinement aurait divisé l'une des nouvelles cellules en 2 en utilisant le noeud 2214. Par conséquent, les noeuds 2212 et 2214 sont « suspendus » par rapport à la cellule 2220 parce que, même si la cellule 2210 a été raffinée (dans une itération précédente), la cellule 2220 n'a pas été raffinée. Étant donné que la cellule 2220 possède au moins deux noeuds suspendus, le moteur AGAR prolongé 110 le raffinera également dans des itérations postérieures. [0082] Au niveau du bloc 2104(v) du procédé 2200, le moteur AGAR prolongé 110 analyse des paires de rebords opposés dans la direction analysée afin de déterminer si seulement l'un des rebords possède 2 ou plusieurs noeuds suspendus. Au niveau du bloc 2104(vi), le moteur AGAR prolongé 110 détermine si une cellule non-échelonnée croise une fracture mais qu'aucun rebord échelonné de la cellule ne croise une fracture. [0083] Si la détermination est « oui » à l'une quelconque des 5 règles décrites ci-dessus en référence au bloc 2104(i-v), le procédé passe la commande au bloc 2108 où cette cellule sera raffinée. Sinon, si la détermination est « non » à l'une quelconque des 5 règles, la commande est retransmise au bloc 2106. Ce procédé continuera de façon itérative jusqu'à ce que chaque cellule dans le modèle ait été analysée. Une fois l'analyse complétée, le moteur AGAR prolongé 110 va générer un modèle de la Terre en utilisant les cellules raffinées. Le modèle de la Terre généré peut être utilisé dans une diversité d'applications, telle que dans la simulation de réservoir. Ensuite, la simulation de réservoir peut être utilisée pour concevoir et réaliser un achèvement de puits de forage et pour implémenter, de façon conséquente, les opérations de fond de puits. Exemples de techniques de maillage stadimétrique [0084] La FIG. 23 illustre une image 2300 des fractures en 3D qui peuvent être modélisées avec les caractéristiques divulguées. Comme on peut le voir dans l'image 2300, les couches de la formation de la Terre comprennent des fractures à l'intérieur de la formation. Ces fractures peuvent être décrites comme des fissures ouvertes ou des vides à l'intérieur de la formation, et elles peuvent être soit naturelles soit générées artificiellement à partir d'un puits de forage. La compréhension et la modélisation des caractéristiques propres à ces fractures sont importantes étant donné que les fractures permettent et influencent le flux des fluides du réservoir vers le puits de forage. Des images telles que l'image 2300 peuvent être obtenues ou générées en utilisant des catalogues d'image. Les catalogues d'image peuvent utiliser un transducteur rotatif pour mesurer l'impédance acoustique à travers toute la paroi du trou de forage afin d'identifier la présence et les directions des fractures rocheuses, aussi bien que de comprendre la direction de l'inclinaison de la stratigraphie. [0085] La FIG. 24 illustre une perspective large d'un exemple d'un jeu de fractures en 3D 2410 qui croise, au niveau des segments de ligne 2420, un plan en 2D 2400. Comme on peut le voir dans la FIG. 24, un jeu de stades en boucle fermé peut être généré autour de chacun des segments de ligne 2420 et un maillage contraint remplit l'espace restant du plan en 2D 2400. Le jeu de stades en boucle fermé peut être généré en utilisant le procédé 2500 décrit ci-dessous en référence à la FIG. 25. [0086] La FIG. 25 est un organigramme illustrant un procédé 2500 permettant la modélisation des fractures en 3D en utilisant un mailleur stadimétrique en 2,5D. Le procédé 2500 commence au niveau du bloc 2502 en recevant un jeu de fractures avec une géométrie qui a été discrétisée dans un plan en 2D par une collection de segments de ligne. Par ailleurs, le procédé 2500 peut commencer par la réalisation de la discrétisation d'un jeu de fractures dans un plan en 2D par une collection de segments de ligne. La collection de segments de ligne peut représenter le croisement entre le plan en 2D et les fractures géologiques en 3D. Comme le démontre le bloc 2503, chaque fracture peut être représentée par une collection de segments en ligne droite pour faire l'approximation d'une courbure de la fracture. [0087] Conformément à un exemple, au niveau du bloc 2504, pour chaque segment de la ligne de fracture dans le plan en 2D, le procédé 2500 réalise le bloc 2506 pour générer un jeu de stades avec un rayon spécifié provenant du segment de la ligne de fracture respectif. Ensuite, au niveau du bloc 2508, le procédé 2500 génère ensuite des boucles fermées autour de tous les segments de ligne d'une fracture. Dans certains modes de réalisation, le bloc 2508 comprend un procédé permettant de générer les boucles fermées autour des segments de ligne de la fracture. Ceci peut comprendre le calcul d'un croisement de tous les côtés des stades pour chaque rayon spécifié pour chaque segment de ligne de la fracture, comme le démontre le bloc 2509 dans la FIG. 25, et le rejet des segments contenus pour chaque segment de ligne dans chaque segment de ligne de fracture qui est totalement contenu par le stade d'autres segments de ligne dans le segment de la ligne de fracture, comme le démontre le bloc 2511. [0088] Après l'achèvement du bloc 2508 (et les blocs optionnels 2509 et 2511 dans certains modes de réalisation), le procédé 2500 passe au bloc 2510. Au niveau du bloc 2510, des éléments de forme peuvent être générés à l'intérieur des boucles fermées du segment de la ligne droite. Par ex., dans un mode de réalisation, le procédé 2500 génère des segments paramétriques le long d'une longueur et d'un rayon de chaque segment de ligne droite, tel que le montre le bloc 2513. [0089] Le procédé 2500 forme ensuite des éléments quadrilatères là où c'est possible à l'intérieur de la région structurée au niveau du bloc 2515, et forme des polygones à l'intérieur des régions restantes des boucles fermées au niveau du bloc 2517. [0090] Après avoir généré les éléments de forme, le procédé 2500 peut générer un maillage contraint autour des boucles fermées du jeu de segments de ligne de fracture, remplissant le restant du plan en 2D au niveau du bloc 2512). Dans un mode de réalisation, un algorithme de triangulation de Delaunay peut être utilisé pour générer le maillage contraint autour des boucles fermées du jeu de segments de la ligne de fracture. Ainsi, le plan en 2D consiste maintenant entièrement d'éléments de cellule des segments de la ligne de fracture et du maillage contraint. [0091] À ce point, le procédé 2500 peut extruder chacune des cellules dans le plan en 2D en une 3e dimension permettant de créer une ou plusieurs couches de cellules en 3D. Les cellules à l'intérieur d'une boucle fermée d'un segment de la ligne de fracture peuvent représenter une fracture en 3D, alors que les cellules à l'intérieur du maillage peuvent représenter des couches rocheuses englobant ces fractures. Ainsi, le procédé 2500 peut attribuer des propriétés de réservoir telles que, sans limitation, la porosité et la perméabilité, à chacune des cellules en 3D afin de modéliser le flux de fluide du réservoir, tel que le démontre le bloc 2516. Finalement, le procédé 2500 peut saisir le modèle cellulaire en 3D dans un programme de simulation, tel que, sans limitation le logiciel de simulation de réservoir Nexus®, pour la réalisation d'une simulation numérique et pour évaluer le flux de fluide, tel que le démontre le bloc 2518). [0092] La FIG. 26 donne une vue plus rapprochée du procédé 2500 en illustrant un exemple pour la génération d'un maillage informatique autour d'un segment de ligne de fracture unique en conformité avec les modes de réalisation divulguée. En commençant avec le diagramme 2602, un jeu de stades est généré autour d'un segment de ligne 2600. Comme le démontre le diagramme 2602, chaque stade dans le jeu de stades est composé de deux côtés linéaires connectés par deux arcs afin d'entourer complètement le segment de ligne droite. La distance à partir de chaque côté vers le segment à ligne droite est un rayon constant. Dans certains modes de réalisation, la distance du rayon peut être une valeur de variable modifiable par l'utilisateur. [0093] Dans le diagramme 2604, les segments paramétriques le long d'une longueur et du rayon de chaque segment de ligne droite sont générés en conformité avec le bloc 2513 du procédé 2500. Des éléments quadrilatères peuvent ensuite être formés là où c'est possible à l'intérieur de la région structurée comme il est référencé par le bloc 2513 du procédé 2500. Le diagramme 2608 illustre le maillage contraint généré autour des boucles fermées de la ligne de segment 2600. [0094] La FIG. 27 procure une autre vue rapprochée du procédé 2500 en illustrant un exemple de la génération de maillages informatiques autour des segments de ligne de fracture croisée. Par ex., le diagramme 2702 illustre un jeu de stade généré autour de 3 segments de ligne de fracture. Le résultat du diagramme 2702 nécessite que le procédé calcule un croisement de tous les côtés des stades pour chaque rayon spécifié pour chacun du segment de la ligne de fracture croisée telle que référencé dans le bloc 2509 et de rejeter des segments contenus pour chaque segment de ligne de fracture qui sont complètement contenues par des stades d'autres 27 segments de ligne de fracture tel que référencé dans le bloc 2511. [0095] Le diagramme 2704 illustre les résultats de la génération des éléments de forme à l'intérieur des boucles fermées des segments de la ligne de fracture tel que référencé dans le bloc 2510 de la FIG. 25. Comme on peut le voir, les segments paramétriques le long d'une longueur et d'un rayon de chaque segment de ligne de fracture sont générés en conformité avec le bloc 2513 de la FIG. 25. Des éléments quadrilatères peuvent ensuite être formés là où c'est possible à l'intérieur de la région structurée comme il est référencé par le bloc 2515 illustré dans le diagramme 2706. En outre, des polygones peuvent être formés à l'intérieur des régions restantes des boucles fermées des segments de la ligne de fracture tel qu'il est indiqué au niveau du bloc 2517 de la FIG. 25. Le diagramme 2708 illustre un maillage contraint généré autour des boucles fermées des segments de la ligne de fracture tel qu'il est référencé au niveau du bloc 2512 du procédé 2500. [0096] La FIG. 28 illustre un exemple de maillages informatiques générés autour d'un réseau complexe de segments de ligne de fracture. En particulier, le diagramme 2802 indique un jeu de fractures avec une géométrie qui a été discrétisée dans un plan en 2D par une collection de segments de ligne. Le diagramme 2804 illustre les résultats d'un jeu de stades qui est généré autour de chacun des segments de la ligne de fracture. Le diagramme 2806 illustre une vue explosée des segments de la ligne de fracture découlant de l'exécution du processus restant 2500 décrite ci-dessus en référence à la FIG. 25. [0097] Comme illustrée dans la Fig. 28, l'algorithme divulgué peut rapidement générer des grilles non-structurées en utilisant des éléments structurés autour de géométries complexes. Les cellules dans le plan en 2D illustrées dans la FIG. 28 peuvent ensuite être extrudées dans une 3e dimension pour former des couches de cellules en 3D. On peut attribuer aux cellules en 3D des propriétés de réservoir pour permettre une simulation numérique. [0098] Comme autre exemple, la FIG. 29 présente une illustration 2900 de géométrie complexe impliquant une pluralité de segments de ligne de fracture croisée dans laquelle certains modes de réalisation peuvent rapidement générer une cellule de grille en 2D qui peut être extrudée dans des éléments en 3D pour réaliser des simulations numériques. Exemples d'implémentation de système informatique [0099] Même si des exemples de modes de réalisation ont été décrits en termes d'appareils, de systèmes, de services et de procédés, il est envisagé que certaines fonctionnalités décrites ici puissent être implémentées dans des logiciels sur des microprocesseurs, tels que la puce de microprocesseur intégrée dans des dispositifs informatiques, tels que le système informatique 3000 illustré dans la FIG. 30. Dans divers modes de réalisation, une ou plusieurs des fonctions des divers composants peuvent être implémentées dans un logiciel qui commande un dispositif informatique, tel que le système informatique 3000, qui est décrit ci-dessous en référence à la FIG. 30. [0100] Afin d'implémenter les diverses caractéristiques et fonctions décrites ci-dessus, certains ou tous les éléments des systèmes (par ex., le système 100 de la FIG. 1), les flux de travail et les procédés (par ex., les procédés 200, 2100 et 2500 de la FIG. 2, 21, et 25) peuvent être implémentés en utilisant des éléments du système informatique de la FIG. 30. Plus particulièrement, la FIG. 30 illustre un exemple d'un système informatique 3000 permettant d'implémenter les techniques en conformité avec la présente divulgation. [0101] Les aspects de la présente divulgation illustrés dans les FIG. 1-28, ou une quelconque ou plusieurs parties et une ou plusieurs fonctions de ceux-ci, peuvent être implémentés en utilisant un matériel, des modules logiciels, un micrologiciel, un support tangible lisible par ordinateur ayant une logique ou des instructions stockées sur celui-ci, ou une combinaison de ceux-ci et peuvent être implémentés dans un ou plusieurs systèmes informatiques ou d'autres systèmes de traitement. [0102] La FIG. 30 illustre un exemple d'un système informatique 3000 dans lequel les modes de réalisation de la présente divulgation, ou des parties de ceux-ci, peuvent être implémentés sous forme d'instructions ou de code lisible par ordinateur. Par ex., certaines fonctionnalités réalisées par le système 100 illustré dans la FIG. 1, peuvent être implémentées dans le système informatique 3000 en utilisant un matériel, un logiciel, un micrologiciel, un support non-transitoire lisible par ordinateur ayant des instructions stockées sur celui-ci, ou une combinaison de ceux-ci et peut être implémenté dans un ou plusieurs systèmes informatiques ou d'autres systèmes de traitement. Un matériel, un logiciel, ou toute combinaison de ceux-ci peut incarner certains modules et composants utilisés pour implémenter le système 100 et le flux de travail 200 illustrés dans les FIG. 1 et 2 présentés ci- dessus. De la même façon, un matériel, un logiciel, ou toute combinaison de ceux-ci peut incarner certains modules et composants utilisés pour implémenter les étapes dans les organigrammes illustrés dans les FIG. 21, 22 et 25 présentés ci-dessus. [0103] Si une logique programmable est utilisée, une telle logique peut s'exécuter sur une plate-forme de traitement disponible dans le commerce ou sur un dispositif spécialisé. Un homme de métier peut apprécié que les modes de réalisation de l'objet de la divulgation peut être pratiqué avec diverses configurations de systèmes informatiques, comprenant des systèmes multiprocesseurs multi-coeurs, des minis ordinateurs, des ordinateurs centraux, des ordinateurs liés ou regroupés avec des fonctions distribuées, aussi bien que des ordinateurs omniprésents ou miniatures qui peuvent être intégrées dans virtuellement tous les dispositifs. [0104] Par ex., au moins un dispositif de processeur et de mémoire peut être utilisé pour implémenter les modes de réalisation décrits ci-dessus. Un processeur peut être un processeur unique, une pluralité de processeurs, ou des combinaisons de ceux-ci. Les processeurs peuvent comporter un ou plusieurs « coeurs » de processeur. [0105] Certains modes de réalisation sont décrits en termes de cet exemple de système informatique 3000. Après la lecture de cette description, il sera évident à un spécialiste du domaine pertinent comment implémenter les exemples divulgués utilisant d'autres systèmes informatiques et/ou d'autres architectures informatiques. Même si les opérations peuvent être décrites sous forme d'un procédé séquentiel, certaines des opérations peuvent, en fait, être réalisées en parallèle, simultanément et/ou dans un environnement distribué, et avec le code de programme stocké localement ou à distance pour un accès par les appareils à processeur unique ou à multiples processeurs. En outre, dans certains modes de réalisation l'ordre des opérations peut être réarrangé sans s'écarter de l'esprit de l'objet divulgué. [0106] Le processeur 3004 peut être un processeur spécialisé ou un processeur polyvalent. Comme il sera compris par les spécialistes du domaine pertinent, le processeur 3004 peut également être un système à processeur unique dans un système multiprocesseur multicceurs, un tel système fonctionnant seul ou dans un regroupement de dispositifs informatiques ou dans un regroupement ou une ferme de serveurs. Le processeur 3004 est relié à une infrastructure de communication 3006, par ex., un bus, une file de message, un réseau ou un schéma de transfert de messages multicceurs. Dans un mode de réalisation, le système de bus 109 du système 100 décrit ci-dessus en référence à la FIG. 1 peut être incarné sous forme d'une infrastructure de communication, 3006 illustrée dans la FIG. 30. Selon un mode de réalisation, le processeur 102 du système 100 décrit ci-dessus en référence à la FIG. 1 peut être incarné sous forme du processeur 3004 illustré dans la FIG. 30. [0107] Le système informatique 3000 comprend également une mémoire principale 3008, par exemple, une mémoire RAM et peut également comprendre une mémoire secondaire 3010. La mémoire secondaire 3010 peut comprendre, par ex., un disque dur 3012, et un disque de stockage amovible 3014. Le disque de stockage amovible 3014 peut comprendre une disquette, un disque à bande magnétique, un disque optique, une mémoire flash, etc. Dans un mode de réalisation non limitant, le stockage 104 du système 100 de la FIG. 1 peut être incarné sous forme de la mémoire principale 3008 illustrée dans la FIG. 30. [0108] Le disque de stockage amovible 3014 lit et/ou écrit sur une unité de stockage amovible 3018 d'une façon bien connue. L'unité de stockage amovible 3018 peut comprendre une disquette, un disque à bande magnétique, un disque optique, etc., qui est lu par le disque de stockage amovible 3014 qui peut également enregistrer sur celui-ci. Comme il sera compris par des spécialistes du domaine pertinent, l'unité de stockage amovible 3018 comprend un support de stockage non-transitoire lisible par ordinateur ayant stocké, sur celui-ci, un logiciel informatique et/ou des données. [0109] Dans des implémentations alternatives, une mémoire secondaire, 3010 peut comprendre d'autres moyens semblables pour permettre au programme informatique ou à d'autres instructions d'être téléchargées dans le système informatique 3000. De tels moyens peuvent comprendre, par ex., une unité de stockage amovible 3022 et une interface 3020. Des exemples de tels moyens peuvent comprendre un module logiciel et un module d'interface (tel que celui trouvé dans les appareils de jeux vidéo), une puce de mémoire amovible (telle qu'une EPROM ou EEPROM) et la prise associée, et d'autres unités de stockage amovibles 3022 et des interfaces 3020 qui permettent le transfert de logiciels et de données de l'unité de stockage amovible 3022 vers le système informatique 3000. [0110] Le système informatique 3000 peut comprendre une interface de communication 3024. L'interface de communication 3024 permet le transfert de logiciel et de données entre un système informatique 3000 et des dispositifs externes. L'interface de communication 3024 peut comprendre un modem, une interface de réseau (telle qu'une carte Éthernet), un port de communication, un emplacement et une carte PCMCIA, ou des éléments semblables. Le logiciel et les données transférés à travers l'interface de communication 3024 peuvent être sous la forme de signaux, qui peuvent être électroniques, électromagnétiques, optiques, ou d'autres signaux capables d'être reçus par l'interface de communication 3024. Ces signaux peuvent être fournis à l'interface de communication 3024 à travers une voie de communication 3026. La voie de communication 3026 transporte des signaux et peut être implémenté avec un fil ou un câble, la fibre optique, une ligne de téléphone, une liaison de téléphone cellulaire ou une liaison RF ou d'autres canaux de communication. [0111] Tels qu'ils sont utilisés ici, les termes « support lisible par ordinateur » et « support non-transitoire lisible par ordinateur » sont généralement utilisés pour décrire des supports tels que des mémoires, telles que la mémoire principale 3008 et une mémoire secondaire 3010, qui peuvent être des mémoires à semi-conducteurs (par ex., DRAM, etc.). Un support lisible par ordinateur et un support non-transitoire lisible par ordinateur peut également décrire une unité de stockage amovible 3018, une unité de stockage amovible 3022, et un disque dur installé dans le disque dur 3012. Les signaux transportés à travers la voie de communication 3026 peuvent également incarner la logique décrite ici. Ces produits de programme informatique peuvent également constituer un moyen de fournir un logiciel au système informatique 3000. [0112] Les programmes informatiques (également appelés logique de commande d'ordinateur) peuvent être stockés dans la mémoire principale 3008 et/ou dans la mémoire secondaire 3010. Les programmes informatiques peuvent également être reçus à travers l'interface de communication 3024. De tels programmes informatiques, lorsqu'ils sont exécutés, permettent au système informatique 3000 d'implémenter les présents exemples tels que décrits ici. En particulier, les programmes informatiques, lorsqu'ils sont exécutés, permettent au processeur 3004 d'implémenter les procédés de la présente divulgation, tels que les étapes dans le flux de travail illustré dans la FIG. 2 et les procédés illustrés par les organigrammes des FIG. 21, 22 et 25, présentés ci-dessus. Par conséquent, de tels programmes informatiques représentent les commandes du système informatique 3000. Lorsque des caractéristiques sont implémentées en utilisant un logiciel, celui-ci peut être stocké dans un programme informatique et téléchargé dans le système informatique 3000 en utilisant un disque de stockage amovible 3014, une interface 3020 et un disque dur 3012, ou une interface de communication 3024. [0113] Dans certains modes de réalisation, un ou plusieurs des dispositifs d'affichage 108 d'un système 100 et des écrans utilisés pour afficher les interfaces utilisateur illustrées dans les FIG. 8 et 9 pourraient être moniteurs d'ordinateur 3030 illustré dans la FIG. 30. Le moniteur d'ordinateur 3030 du système informatique 3000 peut être implémenté sous la forme d'un dispositif tactile (c.-à-d., un écran tactile). Le moniteur d'ordinateur 3030 peut se connecter à l'infrastructure de communication à travers une interface d'affichage 3002 afin d'afficher le contenu électronique reçu. Par ex., le moniteur d'ordinateur 3030 peut être utilisé pour afficher des modèles d'entrée et des fractures raffinées. Mais également, par ex., le moniteur d'ordinateur 3030 peut être utilisé pour afficher les interfaces utilisateur 800 et 900 illustrées dans les FIG. 8 et 9. [0114] Certains modes de réalisation peuvent être orientés également vers un programme informatique comprenant un logiciel stocké sur un quelconque support utilisable par un ordinateur. Un tel logiciel, lorsqu'il est exécuté sur un ou plusieurs dispositifs de traitement de données, entraîne le fonctionnement d'un ou des dispositifs de traitement de données tel qu'il est décrit ici. Certains modes de réalisation utilisent un support lisible ou utilisable par un ordinateur. Des exemples de supports utilisables par un ordinateur comprennent, sans limitation, des dispositifs de stockages principaux (par ex., tout type de mémoire à accès aléatoire), des dispositifs de stockage secondaires (par ex., des disques durs, des disquettes, des CD ROM, des disques ZIP, des bandes, des dispositifs de stockage magnétique et des dispositifs de stockage optique, des MEMS, des dispositifs de stockage nanotechnologique, etc.), et des supports de communication (par ex., des réseaux de communication câblés ou non, des réseaux locaux, des réseaux étendus, et l'intranet, etc.). [0115] Bien que l'objet de la présente divulgation ait été décrit en détail par rapport aux modes de réalisation spécifiques de celui-ci, il sera apprécié que les spécialistes du domaine, lorsqu'ils auraient acquis une compréhension de la description précédente puissent facilement apporter des altérations à, et des variations de, et des équivalents de tels modes de réalisation. Par conséquent, il doit être compris que la présente divulgation a été présentée dans un objectif d'exemple plutôt que de limitation, et n'écarte pas l'inclusion de telles modifications, variations et/ou ajouts à l'objet de la présente divulgation comme il serait facilement évident à un homme de métier. En effet, les procédés et les systèmes décrits ici peuvent être incarnés dans une variété d'autres formes ; en outre, diverses omissions, substitutions et changements dans la forme des procédés et de systèmes décrits ici peuvent être apportés sans s'écarter de l'esprit de la présente divulgation. Les revendications ci-jointes et leurs équivalents sont destinés à couvrir de telles formes ou modifications qui font partie de la portée et de l'esprit de la présente divulgation. [0116] Sauf en cas de mention contraire spécifique, il est apprécié qu'à travers toute cette spécification les discussions utilisant les termes tels que « en traitement », « calculer », « calcul », « détermination » et « identification » et des termes semblables décrivent des actions ou des procédés d'un dispositif informatique, tel qu'un ou plusieurs ordinateurs ou un ou des autres dispositifs de calcul électronique semblables, qui manipulent ou transforment des données représentées sous la forme de quantités électroniques ou magnétiques à l'intérieur des mémoires, des registres ou d'autres dispositifs de stockage d'informations, des dispositifs de transformation ou des dispositifs d'affichage de la plate-forme informatique. [0117] Le ou les systèmes présentés ici ne sont pas limités à une quelconque architecture ou configuration matérielle. Un dispositif informatique peut comprendre un quelconque agencement de composants qui procure un résultat dépendant d'une ou de plusieurs entrées. Des dispositifs informatiques appropriés comprennent des systèmes informatiques basés sur des microprocesseurs polyvalents qui accèdent à un logiciel stocké qui programme ou configure le système informatique d'un appareil de calcul polyvalent vers un appareil de calcul spécialisé implémentant l'un ou plusieurs des modes de réalisation de l'objet de la présente divulgation. Tout type de programmation, de script ou tout autre type de langage ou de combinaison de langage peut être utilisée pour implémenter les enseignements contenus ici dans un logiciel qui doit être utilisé dans la programmation ou la configuration d'un dispositif informatique. [0118] Des exemples des procédés décrits ici peuvent être réalisés dans le fonctionnement de tels dispositifs informatiques L'ordre des blocs présentés dans les exemples susmentionnés peut être varié, par ex., les blocs peuvent être réaménagés, combinés et/ou divisés en sous-blocs. Certains blocs ou procédés peuvent être exécutés en parallèle. [0119] Le langage conditionnel utilisé ici, tel que, entre autres « peut », « pourrait », « pouvait », « peut », « par ex. », et des mots semblables, sauf en cas de mention contraire spécifique, ou autrement compris dans le contexte utilisé, est généralement destiné à transmettre que certains exemples comprennent, alors que d'autres exemples ne comprennent pas, certaines caractéristiques, éléments et/ou étapes. Ainsi, un tel langage conditionnel, n'est pas généralement destiné à impliquer que des caractéristiques, des éléments et/ou des étapes sont d'une façon quelconque nécessaires pour un ou plusieurs exemples ou qu'un ou plusieurs des exemples comprennent nécessairement une logique pour la prise de décision, avec ou sans une entrée de l'auteur ou une invite, que ses caractéristiques, éléments et/ou étapes soient incluses ou non ou qu'elles doivent être exécutées dans un exemple donné. [0120] Les termes « comprenant », « incluant », « ayant » et des termes semblables sont synonymes et sont utilisés inclusivement, d'une façon ouverte et n'excluent pas des éléments, caractéristiques, actions, opérations additionnelles, etc. Mais également, le terme « ou » est utilisé dans son sens inclusif (et non pas dans son sens exclusif) de sorte que lorsqu'il est utilisé, par ex., pour relier une liste d'éléments, le terme « ou » veut dire un, quelques-uns ou tous les éléments dans la liste. En outre, l'utilisation de « basé sur » est destiné à être ouvert et inclusif, en ce qu'un procédé, une étape, un calcul ou tout autre action « basée sur » une ou plusieurs conditions ou valeurs décrites peuvent, en pratique, être basé sur des conditions ou des valeurs additionnelles à par celle décrites. [0121] Les diverses caractéristiques et procédés décrits ci-dessus peuvent être utilisés indépendamment l'un de l'autre, ou peuvent être associé de diverses façons. Toutes les combinaisons possibles et les sous-combinaisons sont destinées à faire partie de la portée de la présente divulgation. En outre, certains procédés ou blocs de procédé peuvent être omis dans certaines implémentations. Les procédés et les processus décrits ici ne sont également pas limités à une quelconque séquence donnée, et les blocs ou les états apparentés à ceux-ci peuvent être réalisés dans d'autres séquences qui sont appropriées. Par ex., les blocs et les états décrits peuvent être réalisés dans un ordre autre que celui qui est spécifiquement divulgué, ou de multiples blocs ou états peuvent être combinés en un bloc ou état unique. Les exemples de blocs ou d'états peuvent être réalisés en série, en parallèle, ou d'une autre façon. Des blocs ou des états peuvent être ajoutés ou enlevés des exemples divulgués. De la même façon, les exemples de systèmes et de composants décrits ici peuvent être configurés différemment de la façon décrite. Par ex., des éléments peuvent être ajoutés à, ou retranchés de, ou réarrangés en comparaison aux exemples divulgués. [0122] La description précédente des modes de réalisation, y compris les modes de réalisation illustrés, a été présentée seulement dans un but illustratif et il n'est pas envisagé que la description soit exhaustive ou qu'elle se limite aux formes précises divulguées. De nombreuses modifications, adaptations et utilisations de celle-ci seront apparentes à un spécialiste du domaine. [0123] Les modes de réalisation et les procédés décrits ici concernent également l'un ou plusieurs des paragraphes suivants. Tel qu'il est utilisé ci-dessus, toute référence à une série 25 d'exemples doit être comprise comme une référence à chacun de ces exemples de façon disjonctive (par ex., « Exemples 1-4 » doit être compris comme « Exemples 1, 2, 3 ou 4 »). [0124] L'Exemple 1 est un procédé implémenté par ordinateur qui comprend la réception, au niveau d'un dispositif informatique, d'une spécification de réservoir. Le procédé identifie, en fonction de la spécification du réservoir, un jeu de fractures comprenant : Des fractures 30 (2,5D) permettant la 2,5D ; et d'autres fractures. Le procédé comprend également la génération d'un modèle de réservoir intermédiaire comprenant un maillage d'extrusion qui modélise les fractures permettant le 2,5D dans un espace en 3D. En réponse à la détermination que les cellules dans le maillage doivent être raffinées dans une direction à l'intérieur de l'espace en 3D, le procédé raffine les cellules de façon anisotrope dans le maillage correspondant à d'autres fractures. Le procédé comprend également la résolution, par le dispositif informatique, d'un réseau de fracture à l'intérieur du modèle de réservoir intermédiaire en utilisant les cellules raffinées ; et en générant un modèle de réservoir de Terre à l'aide du réseau de fracture. [0125] L'Exemple 2 est le procédé de l'Exemple 1, dans lequel les fractures permettant la 2,5D possèdent une géométrie qui a été discrétisée dans un plan en 2D par une collection de segments de ligne, et dans lequel la génération du modèle de réservoir intermédiaire comprend, pour chaque segment de ligne associé avec chaque fracture dans les fractures permettant la 2,5D : la génération d'un jeu de stades avec un rayon spécifié à partir du segment de ligne ; la génération de boucles fermées autour de tous les segments de ligne associés à une fracture permettant la 2,5D ; et la génération des éléments de forme à l'intérieur des boucles fermées de la ligne de segment. Dans l'Exemple 2, le procédé comprend également la génération de maillage sous forme de maillage contraint autour des boucles fermées des fractures permettant la 2,5D pour remplir un espace restant dans le plan en 2D. [0126] L'Exemple 3 est le procédé des Exemples 1 ou 2 dans lequel le procédé comprend également l'utilisation du modèle de réservoir de terre dans une simulation de réservoir, et ensuite l'affichage de la simulation de réservoir sur un écran du dispositif informatique. [0127] L'Exemple 4 est le procédé des Exemples 1-3, dans lequel : les fractures permettant la 2,5D comprennent des fractures représentées par un ou plusieurs plans sensiblement verticaux ; et les autres fractures comprennent des fractures représentées par un ou plusieurs plans sensiblement horizontaux. [0128] L'Exemple 5 est le procédé des Examples 1-4, dans lequel le raffinement anisotrope des cellules dans le maillage comprend l'identification d'une direction à l'intérieur de l'espace 3D dans laquelle les cellules doivent être raffinées, et la division d'un rebord des cellules, le rebord étant dans la direction à l'intérieur de l'espace 3D. [0129] L'Exemple 6 est le procédé de l'Exemple 5, dans lequel le rebord est divisé le long d'un plan qui est perpendiculaire à un axe d'une fracture à l'intérieur d'un réseau de fracture. [0130] L'Exemple 7 est le procédé des Exemples 1-6, dans lequel le raffinement anisotrope des cellules comprend : la détermination d'un emplacement d'une fracture à l'intérieur d'un modèle de réservoir intermédiaire ; l'identification de cellules dans la région proche de la fracture ; la détermination d'une taille de cible des rebords des cellules proches de la fracture ; et la détermination si les rebords des cellules proches de la fracture dans la direction à l'intérieur de l'espace 3D sont plus longs que la taille de la cible multipliée par une variable scalaire. Les cellules proches de la fracture possèdent des rebords dans la direction à l'intérieur de l'espace 3D qui sont plus longs que la taille de la cible multipliée par la variable scalaire sont parmi les cellules qui sont raffinées. [0131] L'Exemple 8 est le procédé de l'Exemple 7, dans lequel la détermination d'un emplacement de la fracture comprend l'échelonnement des rebords de cellules proches de la fracture qui sont plus longs que la taille de la cible multipliée par la variable scalaire, et la détermination si les rebords échelonnés croisent la fracture. Les cellules proches de la fracture qui possèdent des rebords échelonnés qui croisent la fracture sont des cellules qui sont raffinées. [0132] L'Exemple 9 est le procédé de l'Exemple 7 ou 8, dans lequel la détermination d'un emplacement de la fracture comprend également l'échelonnement des rebords des cellules proches de la fracture qui sont plus longs que la taille de la cible multipliée par la variable scalaire, et pour une quelconque paire opposée des bords échelonnés, la détermination si les rebords échelonnés croisent la fracture. Les cellules proches de la fracture qui ont exactement un rebord 10 qui croise la fracture sont les cellules qui sont raffinées. [0133] L'Exemple 10 est le procédé de l'Exemple 7-9, dans lequel la détermination de l'emplacement de la fracture comprend également la détermination si les rebords de cellules proches de la fracture qui sont plus longs que la taille de la cible multipliée par la variable scalaire, comportent deux ou plusieurs noeuds suspendus, et dans lequel les cellules proches de la fracture qui comportent deux ou plusieurs noeuds suspendus sont raffinées. [0134] L'Exemple 11 est le procédé de l'Exemple 7-10, dans lequel le raffinement anisotrope des cellules dans le maillage comprend également, pour une quelconque paire opposée des rebords des cellules proches de la fracture qui sont plus longs que la taille de la cible multipliée par la variable scalaire, la détermination si exactement un des deux rebords comportent deux ou plusieurs noeuds suspendus, et dans lequel les cellules proches de la fracture qui comportent deux ou plusieurs noeuds suspendus sont raffinées. [0135] L'Exemple 12 est le procédé de l'Exemple 11, dans lequel dans lequel le raffinement anisotrope des cellules dans le maillage comprend l'échelonnement des rebords de cellules proches de la fracture qui sont plus longs qu'une taille de la cible multipliée par la variable scalaire, et la détermination si les cellules proches de la fracture ayant des rebords échelonnés croisent une fracture alors qu'aucune cellules proche de la fracture ne croise une fracture. Les cellules proches de la fracture qui possèdent des rebords échelonnés qui croisent la fracture sont raffinées. [0136] L'Exemple 13 est le procédé des Exemples 7-11, dans lequel le raffinement anisotrope des cellules dans le maillage comprend également : la détermination d'un emplacement d'une fracture à l'intérieur d'un modèle de réservoir intermédiaire ; l'identification des cellules dans une région proche de la fracture ; l'identification d'une direction à l'intérieur de l'espace 3D dans laquelle les cellules proches de la fracture doivent être raffinées ; et la détermination si les cellules proches de la fracture ayant un rebord dans la direction identifiée possèdent également pas plus d'un noeud suspendu. Les cellules proches de la fractures qui n'ont pas plus d'un noeud suspendu sont raffinées. [0137] L'Exemple 14 est le procédé des Exemples 7-11, dans lequel le raffinement anisotrope des cellules comprend également : la détermination d'un emplacement d'une fracture à l'intérieur d'un modèle de réservoir intermédiaire ; l'identification de cellules dans la région proche de la fracture ; et pour une quelconque paire opposée des rebords de cellule proche de la fracture, la détermination si exactement l'un des deux rebords opposés possède deux ou plusieurs noeuds suspendus. Les cellules proches de la fracture qui possèdent deux ou plusieurs noeuds suspendus sont raffinées. [0138] L'Exemple 15 est un support de stockage lisible par ordinateur qui a des instructions exécutables stockées sur celui-ci. Si les instructions sont exécutées par un dispositif informatique, les instructions permettent aux dispositifs informatiques de réaliser les opérations. Les instructions comprennent : des instructions permettant de recevoir une spécification de réservoir. Les instructions comprennent également des instructions permettant d'identifier, en se basant sur la spécification du réservoir, un jeu de fractures comprenant des fractures (2,5D) permettant la 2,5D, et d'autres fractures. Les instructions comprennent également des instructions pour générer un modèle de réservoir intermédiaire comprenant un maillage d'extrusion qui modélise les fractures permettant le 2,5D dans un espace en 3D. En réponse à la détermination que les cellules dans le maillage doivent être raffinées dans une direction à l'intérieur de l'espace 3D, les instructions comprennent des instructions pour raffiner les cellules de façon anisotrope dans le maillage correspondant aux autres fractures. Les instructions comprennent également les instructions pour la résolution d'un réseau de fracture à l'intérieur du modèle de réservoir intermédiaire en utilisant les cellules raffinées, et des instructions pour générer un modèle de réservoir de Terre à l'aide du réseau de fracture. [0139] L'Exemple 16 est le support de stockage de l'Exemple 15, dans lequel les instructions pour le raffinement anisotrope des cellules dans le maillage comprennent des instructions permettant d'identifier une direction à l'intérieur de l'espace 3D dans lequel les cellules doivent être raffinées, et des instructions pour la division d'un rebord d'une cellule, le rebord étant dans la direction à l'intérieur de l'espace 3D. [0140] L'Exemple 17 est le support de stockage des Exemples 15ou 16, où les fractures permettant la 2,5D comprennent des fractures représentées par un ou plusieurs plans sensiblement verticaux, et où les autres fractures comprennent des fractures représentées par un ou plusieurs plans sensiblement horizontaux. [0141] L'Exemple 18 est un système qui comprend un processeur, un dispositif d'affichage, et une mémoire ayant des instructions stockées sur celle-ci, qui, lorsqu'elles sont exécutées par un processeur, entraîne le processeur à exécuter des opérations. Les opérations comprennent la réception d'une spécification de réservoir et l'identification, basée sur la spécification du réservoir, d'un jeu de fractures comprenant des fractures (2,5D) permettant la 2,5D, et d'autres fractures. Les opérations comprennent également la génération d'un modèle de réservoir intermédiaire comprenant un maillage d'extrusion qui modélise les fractures permettant le 2,5D dans un espace en 3D. En réponse à la détermination que les cellules dans le maillage doivent être raffinées dans une direction à l'intérieur de l'espace 3D, les opérations comprennent le raffinement des cellules de façon anisotrope dans le maillage correspondant aux autres fractures. Les opérations comprennent également les instructions pour la résolution d'un réseau de fracture à l'intérieur du modèle de réservoir intermédiaire en utilisant les cellules raffinées, et ensuite de générer un modèle de réservoir de Terre à l'aide du réseau de fracture. [0142] L'Exemple 19 est le système de l'Exemple 18, où les fractures permettant la 2,5D comprennent des fractures représentées par un ou plusieurs plans verticaux, et les autres fractures comprennent des fractures représentées par un ou plusieurs plans horizontaux. [0143] L'Exemple 20 est le système des Exemples 18 ou 19, dans lequel les opérations comprennent également l'affichage du modèle de réservoir de la Terre sur le dispositif d'affichage.30
Claims (17)
- REVENDICATIONS1. Procédé implémenté par ordinateur, caractérisé en ce qu'il comprend : la réception, au niveau d'un dispositif informatique (3000), d'une spécification de réservoir (204) ; l'identification, basée sur la spécification du réservoir (204), d'un jeu de fractures (206, 208) comprenant : des fractures (206) permettant la 2,5D; et d'autres fractures (208) ; la génération d'un modèle de réservoir intermédiaire (212) comprenant un maillage d'extrusion qui modélise les fractures (206) permettant la 2,5D dans un espace 3D; répondant à la détermination que les cellules dans le maillage doivent être raffinées dans une direction à l'intérieur de l'espace 3D, le raffinement des cellules de façon anisotrope dans le maillage correspondant aux autres fractures (208) ; la résolution, par le dispositif informatique, d'un réseau de fracture à l'intérieur du modèle de réservoir intermédiaire (212) en utilisant les cellules raffinées ; et la génération d'un modèle de réservoir de Terre à l'aide du réseau de fracture.
- 2. Procédé implémenté par ordinateur selon la revendication 1, dans lequel les fractures (206) permettant la 2,5D possèdent une géométrie qui a été discrétisée dans un plan en 2D par une collection de segments de ligne, et dans lequel la génération du modèle de réservoir intermédiaire (212) comprend, pour chaque segment de ligne associé à chaque fracture dans les fractures (206) permettant la 2,5D : la génération d'un jeu de stades à un rayon spécifié provenant du segment de ligne ; la génération des boucles fermées autour des segments de ligne de la fracture associés à une fracture (206) permettant la 2,5D ; et la génération des éléments de forme à l'intérieur des boucles fermées du segment de la ligne, le procédé comprenant également : la génération du maillage sous forme de maillage contraint autour des boucles fermées 30 des fractures permettant la 2,5D pour remplir un espace restant dans le plan en 2D.
- 3. Procédé implémenté par ordinateur selon la revendication 1 ou 2, comprenant également :l'utilisation d'un modèle de réservoir de Terre dans une simulation de réservoir ; et l'affichage de la simulation de réservoir sur un écran (108, 3030) du dispositif informatique (3000).
- 4. Procédé implémenté par ordinateur selon l'une quelconque des revendications 1 à 3, dans lequel : les fractures (206) permettant le 2,5D peuvent comprendre des fractures représentées par un ou plusieurs plans sensiblement verticaux ; et les autres fractures comprennent des fractures représentées par un ou plusieurs plans sensiblement horizontaux.
- 5. Procédé implémenté par ordinateur selon l'une quelconque des revendications 1 à 4, dans lequel le raffinement anisotrope des cellules dans le maillage comprend : l'identification d'une direction à l'intérieur de l'espace 3D dans laquelle les cellules doivent être raffinées ; et la division d'un rebord des cellules, le rebord étant dans la direction à l'intérieur de l'espace 3D.
- 6. Procédé implémenté par ordinateur selon la revendication 5, dans lequel le rebord est divisé le long d'un plan qui est perpendiculaire à un axe d'une fracture à l'intérieur d'un réseau de fracture.
- 7. Procédé implémenté par ordinateur selon l'une quelconque des revendications 1 à 6, dans lequel le raffinement anisotrope des cellules comprend : la détermination d'un emplacement d'une fracture (206) dans un modèle de réservoir intermédiaire (212) ; l'identification des cellules dans la région proche de la fracture ; la détermination d'une taille de cible des rebords des cellules proches de la fracture ; et la détermination si les rebords des cellules proches de la fracture (206) dans la direction à l'intérieur de l'espace 3D sont plus longs que la taille de la cible multipliée par un variable scalaire,dans lequel les cellules proches de la fracture ayant des rebords dans la direction à l'intérieur de l'espace 3D qui sont plus longs que la taille de la cible multipliée par le variable scalaire sont parmi les cellules qui sont raffinées.
- 8. Procédé implémenté par ordinateur selon la revendication 7, dans lequel la détermination de l'emplacement de la fracture (206) comprend : l'échelonnement des rebords des cellules proches de la fracture qui sont plus longs que la taille de la cible multipliée par la variable scalaire ; et la détermination si les rebords échelonnés croisent la fracture (206), dans lequel les cellules proches de la fracture qui possèdent des rebords échelonnés qui croisent la fracture (206) sont des cellules qui sont raffinées.
- 9. Procédé implémenté par ordinateur selon la revendication 7, dans lequel la détermination de l'emplacement de la fracture (206) comprend également : l'échelonnement des rebords des cellules proches de la fracture qui sont plus longs que la taille de la cible multipliée par la variable scalaire ; et pour une quelconque paire opposée des rebords échelonnés, la détermination si exactement un des deux rebords échelonnés croisent une fracture (206), dans lequel les cellules proches de la fracture (206) qui ont exactement un rebord échelonné croisant la fracture sont les cellules qui sont raffinées.
- 10. Procédé implémenté par ordinateur selon la revendication 7, dans lequel la détermination de l'emplacement de la fracture (206) comprend également : la détermination si les rebords de cellules proches de la fracture (206) qui sont plus longs que la taille de la cible multipliée par la variable scalaire comportent deux ou plusieurs noeuds suspendus, dans lequel les cellules proches de la fracture (206) ayant deux ou plusieurs noeuds suspendus sont raffinées.
- 11. Procédé implémenté par ordinateur selon la revendication 7, dans lequel le raffinement anisotrope des cellules dans le maillage comprend également, pour une quelconque paire opposée des rebords des cellules proches de la fracture (206) qui sont plus longs que la taille de la cible multipliée par la variable scalaire, la détermination si exactementun des deux rebords comportent deux ou plusieurs noeuds suspendus, et dans lequel les cellules proches de la fracture (206) qui comportent deux ou plusieurs noeuds suspendus sont raffinées.
- 12. Procédé implémenté par ordinateur selon la revendication 11, dans lequel le raffinement anisotrope des cellules dans le maillage comprend également : l'échelonnement des rebords des cellules proches de la fracture (206) qui sont plus longs que la taille de la cible multipliée par la variable scalaire ; et la détermination si les cellules proches de la fracture ayant des rebords échelonnés croisent une fracture (206) alors qu'aucune cellule proche de la fracture ne croise une fracture, dans lequel les cellules proches de la fracture (206) qui possèdent des rebords échelonnés qui croisent la fracture sont raffinées.
- 13. Procédé implémenté par ordinateur selon la revendication 7, dans lequel le raffinement anisotrope des cellules dans le maillage comprend également : la détermination d'un emplacement d'une fracture (206) dans un modèle de réservoir intermédiaire (212) ; l'identification des cellules dans la région proche de la fracture (206) ; l'identification d'une direction à l'intérieur de l'espace 3D dans laquelle les cellules proches de la fracture (206) doivent être raffinées ; et la détermination si les cellules proches de la fracture (206) ayant un rebord dans la direction identifiée ont également pas plus d'un noeud suspendu, dans lequel les cellules proches de la fractures (206) qui n'ont pas plus d'un noeud suspendu sont raffinées.
- 14. Procédé implémenté par ordinateur selon la revendication 7, dans lequel le raffinement anisotrope des cellules dans le maillage comprend également : la détermination d'un emplacement d'une fracture (206) dans un modèle de réservoir intermédiaire (212) ; l'identification des cellules dans la région proche de la fracture (206) ; et pour une quelconque paires opposée des rebords de cellule proche de la fracture (206), la détermination si exactement l'un des deux rebords opposés possède deux ou plusieurs noeuds suspendus,dans lequel les cellules proches de la fracture (206) qui possèdent les deux ou plusieurs noeuds suspendus sont raffinées.
- 15. Support de stockage lisible par ordinateur ayant des instructions exécutables stockées sur ledit support de stockage, qui, lorsqu'elles sont exécutées par un dispositif informatique (3000), permet au dispositif informatique (3000) de réaliser des opérations, caractérisé en ce que les instructions comprennent des instructions pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications 1 à 14.
- 16. Système (100) caractérisé en ce qu'il comprend: un processeur (102) ; un dispositif d'affichage (108); et une mémoire ayant des instructions stockées sur celle-ci, qui, lorsqu'elles sont exécutées par le processeur (102), permet au processeur d'exécuter des opérations comprenant des opérations pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications 1 à 14.
- 17. Système (100) selon la revendication 16, comprenant également : l'affichage, sur le dispositif d'affichage (108), du modèle de réservoir de Terre.20
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/US2014/065177 WO2016076847A1 (fr) | 2014-11-12 | 2014-11-12 | Création de maillages de réservoirs par raffinement étendu anisotrope géométriquement adaptatif de polyèdres |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| FR3028333A1 true FR3028333A1 (fr) | 2016-05-13 |
Family
ID=55806035
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR1559673A Withdrawn FR3028333A1 (fr) | 2014-11-12 | 2015-10-12 | Creation de maillage de reservoir utilisant un raffinement anisotrope prolonge, adapte a la geometrie, d'un polyedre |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US10564317B2 (fr) |
| AR (1) | AR102183A1 (fr) |
| CA (1) | CA2963928C (fr) |
| FR (1) | FR3028333A1 (fr) |
| GB (1) | GB2545608B (fr) |
| WO (1) | WO2016076847A1 (fr) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| MX2016000151A (es) * | 2013-07-02 | 2016-03-01 | Landmark Graphics Corp | Mallado de estadios de 2 , 5d. |
| CA2996269C (fr) | 2015-09-24 | 2021-01-12 | Landmark Graphics Corporation | Simulation de reservoirs fractures en utilisant de multiples mailles |
| US20170184760A1 (en) * | 2015-12-23 | 2017-06-29 | Paradigm Sciences Ltd. | Device, system and method for a structure and stratigraphy preserving transformation of a geological model |
| US10255720B1 (en) * | 2016-10-13 | 2019-04-09 | Bentley Systems, Incorporated | Hybrid mesh from 2.5D and 3D point data |
| CN106991244B (zh) * | 2017-04-13 | 2019-06-28 | 河海大学 | 一种基于图论的裂隙网络连通性及渗流计算的方法 |
| CN108287371A (zh) * | 2018-01-31 | 2018-07-17 | 中南大学 | 直流电阻率无单元法中的背景网格自适应剖分方法 |
| US11156744B2 (en) | 2019-01-10 | 2021-10-26 | Emerson Paradigm Holding Llc | Imaging a subsurface geological model at a past intermediate restoration time |
| US10520644B1 (en) | 2019-01-10 | 2019-12-31 | Emerson Paradigm Holding Llc | Imaging a subsurface geological model at a past intermediate restoration time |
| US11353622B2 (en) * | 2020-01-06 | 2022-06-07 | Saudi Arabian Oil Company | Systems and methods for hydrocarbon reservoir three dimensional unstructured grid generation and development |
| CN111401451B (zh) * | 2020-03-17 | 2023-06-30 | 中国石油化工股份有限公司 | 油藏地质模型中储层构型界面自动识别方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5740342A (en) * | 1995-04-05 | 1998-04-14 | Western Atlas International, Inc. | Method for generating a three-dimensional, locally-unstructured hybrid grid for sloping faults |
| EP1107187A1 (fr) * | 1999-12-10 | 2001-06-13 | Institut Francais Du Petrole | Méthode pour générer un maillage sur une formation hétérogène traversée par une ou plusieurs discontinuités géométriques dans le but de réaliser des simulations |
Family Cites Families (35)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2801710B1 (fr) * | 1999-11-29 | 2002-05-03 | Inst Francais Du Petrole | Methode pour generer un maillage hybride permettant de modeliser une formation heterogene traversee par un ou plusieurs puits |
| AU2001271633A1 (en) * | 2000-06-29 | 2002-01-14 | Object Reservoir, Inc. | System and method for defining and displaying a reservoir model |
| US7369973B2 (en) * | 2000-06-29 | 2008-05-06 | Object Reservoir, Inc. | Method and system for representing reservoir systems |
| US7248259B2 (en) * | 2001-12-12 | 2007-07-24 | Technoguide As | Three dimensional geological model construction |
| FR2870621B1 (fr) * | 2004-05-21 | 2006-10-27 | Inst Francais Du Petrole | Methode pour generer un maillage hybride conforme en trois dimensions d'une formation heterogene traversee par une ou plusieurs discontinuites geometriques dans le but de realiser des simulations |
| US7809537B2 (en) * | 2004-10-15 | 2010-10-05 | Saudi Arabian Oil Company | Generalized well management in parallel reservoir simulation |
| FR2891383B1 (fr) * | 2005-09-26 | 2008-07-11 | Inst Francais Du Petrole | Methode pour simuler des ecoulements de fluides au sein d'un milieu discretise par un maillage hybride |
| US7925482B2 (en) * | 2006-10-13 | 2011-04-12 | Object Reservoir, Inc. | Method and system for modeling and predicting hydraulic fracture performance in hydrocarbon reservoirs |
| US8190414B2 (en) * | 2008-03-26 | 2012-05-29 | Exxonmobil Upstream Research Company | Modeling of hydrocarbon reservoirs containing subsurface features |
| BRPI0919207B1 (pt) * | 2008-09-19 | 2020-08-18 | Chevron U.S.A. I.N.C | método e sistema implementado por computador para uso na modelagem de um sistema de reservatório geomecânico, e, método implementado por computador para uso na modelagem de uma fratura |
| US9026416B2 (en) * | 2008-12-03 | 2015-05-05 | Chevron U.S.A. Inc. | System and method of grid generation for discrete fracture modeling |
| GB2477901A (en) * | 2008-12-03 | 2011-08-17 | Chevron Usa Inc | System and method for predicting fluid flow characteristics within fractured subsurface reservoirs |
| FR2948215B1 (fr) * | 2009-07-16 | 2011-06-24 | Inst Francais Du Petrole | Methode pour generer un maillage hexa-dominant d'un milieu souterrain faille |
| US8494827B2 (en) * | 2009-09-25 | 2013-07-23 | Exxonmobil Upstream Research Company | Method of predicting natural fractures and damage in a subsurface region |
| US9085957B2 (en) * | 2009-10-07 | 2015-07-21 | Exxonmobil Upstream Research Company | Discretized physics-based models and simulations of subterranean regions, and methods for creating and using the same |
| CN102612682B (zh) * | 2009-11-12 | 2016-04-27 | 埃克森美孚上游研究公司 | 用于储层建模和模拟的方法和设备 |
| US8727017B2 (en) * | 2010-04-22 | 2014-05-20 | Exxonmobil Upstream Research Company | System and method for obtaining data on an unstructured grid |
| US8731873B2 (en) * | 2010-04-26 | 2014-05-20 | Exxonmobil Upstream Research Company | System and method for providing data corresponding to physical objects |
| US8463586B2 (en) * | 2010-06-22 | 2013-06-11 | Saudi Arabian Oil Company | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5D unstructured grids |
| WO2012015521A1 (fr) * | 2010-07-29 | 2012-02-02 | Exxonmobil Upstream Research Company | Procédé et système de modélisation d'un réservoir |
| CA2820942C (fr) * | 2010-12-10 | 2016-12-20 | Baker Hughes Incorporated | Procede d'amelioration de simulation de reservoir et d'extraction a partir de reservoirs fractures |
| CN102110183B (zh) * | 2010-12-30 | 2011-09-28 | 中国石油化工股份有限公司胜利油田分公司地质科学研究院 | 一种反映流体沿储层大裂缝窜流的数值模拟方法 |
| US9051825B2 (en) * | 2011-01-26 | 2015-06-09 | Schlumberger Technology Corporation | Visualizing fluid flow in subsurface reservoirs |
| US9194222B2 (en) * | 2011-04-19 | 2015-11-24 | Halliburton Energy Services, Inc. | System and method for improved propped fracture geometry for high permeability reservoirs |
| EP2756382A4 (fr) * | 2011-09-15 | 2015-07-29 | Exxonmobil Upstream Res Co | Opérations matricielles et vectorielles optimisées dans des algorithmes à instructions limitées qui effectuent des calculs eos |
| FR2987138A1 (fr) * | 2012-02-22 | 2013-08-23 | Schlumberger Technology Corp | Titre non renseigne. |
| AU2012376793B2 (en) * | 2012-04-10 | 2015-02-12 | Halliburton Energy Services, Inc. | Method and apparatus for generating seismic pulses to map subterranean fractures |
| NO339744B1 (no) * | 2012-08-10 | 2017-01-30 | Schlumberger Technology Bv | Hybrid, lokal ikke-avstemmende metode for flerfasestrømningssimuleringer i heterogent sprukne medier |
| US20140136171A1 (en) * | 2012-11-13 | 2014-05-15 | Chevron U.S.A. Inc. | Unstructured Grids For Modeling Reservoirs |
| US20140236559A1 (en) * | 2013-02-18 | 2014-08-21 | Saudi Arabian Oil Company | Systems, methods, and computer-readable media for modeling complex wellbores in field-scale reservoir simulation |
| US20160170085A1 (en) * | 2013-07-02 | 2016-06-16 | Landmark Graphics Corporation | 2.75d meshing algorithm |
| MX2016000151A (es) * | 2013-07-02 | 2016-03-01 | Landmark Graphics Corp | Mallado de estadios de 2 , 5d. |
| CA2951235C (fr) | 2014-07-11 | 2021-04-13 | Landmark Graphics Corporation | Affinement adaptatif en geometrie anisotrope pour la creation de maillage de reservoirs |
| US10571601B2 (en) * | 2014-10-09 | 2020-02-25 | Schlumberger Technology Corporation | Geologic fault representation for flow simulation |
| US10961834B2 (en) * | 2015-11-09 | 2021-03-30 | Halliburton Energy Services, Inc. | Fracture network fluid flow simulation with junction area modeling |
-
2014
- 2014-11-12 US US15/517,006 patent/US10564317B2/en active Active
- 2014-11-12 GB GB1705594.8A patent/GB2545608B/en active Active
- 2014-11-12 WO PCT/US2014/065177 patent/WO2016076847A1/fr active Application Filing
- 2014-11-12 CA CA2963928A patent/CA2963928C/fr active Active
-
2015
- 2015-10-06 AR ARP150103218A patent/AR102183A1/es unknown
- 2015-10-12 FR FR1559673A patent/FR3028333A1/fr not_active Withdrawn
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5740342A (en) * | 1995-04-05 | 1998-04-14 | Western Atlas International, Inc. | Method for generating a three-dimensional, locally-unstructured hybrid grid for sloping faults |
| EP1107187A1 (fr) * | 1999-12-10 | 2001-06-13 | Institut Francais Du Petrole | Méthode pour générer un maillage sur une formation hétérogène traversée par une ou plusieurs discontinuités géométriques dans le but de réaliser des simulations |
Also Published As
| Publication number | Publication date |
|---|---|
| US20170299770A1 (en) | 2017-10-19 |
| US10564317B2 (en) | 2020-02-18 |
| CA2963928A1 (fr) | 2016-05-19 |
| GB201705594D0 (en) | 2017-05-24 |
| CA2963928C (fr) | 2019-06-25 |
| WO2016076847A1 (fr) | 2016-05-19 |
| AR102183A1 (es) | 2017-02-08 |
| GB2545608B (en) | 2020-06-17 |
| GB2545608A (en) | 2017-06-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR3028333A1 (fr) | Creation de maillage de reservoir utilisant un raffinement anisotrope prolonge, adapte a la geometrie, d'un polyedre | |
| FR3041803A1 (fr) | ||
| Pellerin et al. | Automatic surface remeshing of 3D structural models at specified resolution: A method based on Voronoi diagrams | |
| FR2987903A1 (fr) | Structures de failles geologiques contenant des non-conformites. | |
| EP2924232A2 (fr) | Procédé de construction d'un maillage optimise pour la simulation de reservoir dans une formation souterraine | |
| FR3055723A1 (fr) | Modelisation basee sur un point-vecteur des proprietes de reservoir de petrole pour un modele de simulation de reservoir sans grille | |
| FR3048300A1 (fr) | Representation hybride geocellulaire en 3d de sous-ensembles de reseau de fracture naturelle choisis | |
| US10324228B2 (en) | Gridless simulation of a fluvio-deltaic environment | |
| Alfaleh et al. | Topological data analysis to solve big data problem in reservoir engineering: Application to inverted 4D seismic data | |
| FR3027134A1 (fr) | Utilisation d'un volume elementaire representatif pour determiner un volume de sous-ensemble dans un modele terrestre d'une zone d'interet | |
| AU2013393306B2 (en) | 2.5D stadia meshing | |
| EP2873989A1 (fr) | Procédé de construction d'une grille représentative de la distribution d'une propriété physique d'une formation souterraine par simulation statistique multipoint conditionnelle | |
| CA2965871C (fr) | Definition de facies petrographique non lineaire pour un modele de simulation de reservoir | |
| GB2529957B (en) | 2.75D meshing algorithm | |
| EP3017430B1 (fr) | Algorithme de génération de volume pour maillage de réseau discret | |
| US10310116B2 (en) | Method for determining a stacked model describing architectural elements | |
| Katzmayr et al. | An iterative algorithm for generating constrained Voronoi grids | |
| US10605055B2 (en) | Integrated hydrocarbon fluid distribution modeling | |
| Aydin et al. | Exploring structural uncertainty using a flow proxy in the depositional domain | |
| EP3058398B1 (fr) | Procédé et dispositif pour la détermination contrôlée de canaux | |
| WO2015055905A1 (fr) | Procédé et dispositif de détermination de groupes de faciès géologiques | |
| Tarrahi et al. | Markov-Chain Monte Carlo with Locally Varying Mean Kriging for Improved Reservoir Model Calibration and Uncertainty Assessment | |
| Bahar et al. | Integrated lithofacies and petrophysical properties simulation | |
| FR2983615A1 (fr) | Rotations a partir de directions par gradient |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PLFP | Fee payment |
Year of fee payment: 2 |
|
| PLFP | Fee payment |
Year of fee payment: 3 |
|
| PLSC | Publication of the preliminary search report |
Effective date: 20180302 |
|
| PLFP | Fee payment |
Year of fee payment: 4 |
|
| PLFP | Fee payment |
Year of fee payment: 5 |
|
| ST | Notification of lapse |
Effective date: 20210605 |