Brevet FR3041803A1 — publié le 20170331

SIMULATION DE RÉSERVOIRS FRACTURÉS À L'AIDE DE MAILLAGES MULTIPLES DOMAINE TECHNIQUE La présente invention se rapporte généralement à un système et à un procédé de génération d'une grille pouvant être utilisée pour construire un modèle de simulation d'un réservoir de subsurface, et plus particulièrement, à un système et à un procédé configurés pour modéliser des fractures géologiques. ARRIÈRE-PLAN 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 fournir une description statique d'un réservoir avant la production. Des modèles de simulation du réservoir peuvent être créés pour simuler le flux de fluide à l'intérieur du réservoir au cours de sa durée de vie de production. 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. 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 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. 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 1 3041803 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 en dimension 2,5 (2,5D) inhérente des techniques d'extrusion existantes, des limitations semblables s'appliquent à ces techniques. Par ailleurs, il existe des techniques 5 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. En conséquence, une simulation des réservoirs avec de grands systèmes de fracture, 10 présentant une géométrie et une orientation arbitraires, est difficile, et il faut faire des compromis entre une précision suffisante et un temps informatique raisonnable. À titre d'exemple, une simulation d'un réservoir de shale peut générer un réseau de fractures naturelles constitué de dizaines de milliers de fractures définies géométriquement par des centaines de milliers de triangles. Des systèmes complexes de ce genre ne sont pas facilement 15 résolus à l'aide d'un système de maillage d'un réservoir en 2,5D et lorsqu'ils sont résolus en trois dimensions, produisent souvent des modèles de simulation de grandeur prohibitive qui prennent un temps important à simuler. BRÈVE DESCRIPTION DES DESSINS Des modes de réalisation sont décrits en détail ci-après en faisant référence aux figures 20 annexées, dans lesquelles : Les figures 1 A à 1C sont des diagrammes de flux de haut niveau qui décrivent un procédé à mailles multiples pour modéliser un réservoir d'hydrocarbures fracturé à l'aide d'un système à mailles multiples selon un ou plusieurs modes de réalisation ; La figure 2 est une perspective de la surface supérieure d'une définition de maillage structuré 25 ou extrudé d'un réservoir à titre d'exemple fourni comme un modèle de la terre structuré, selon un ou plusieurs modes de réalisation ; La figure 3 est une vue en perspective d'un modèle de fracture hydraulique à titre d'exemple défini à l'aide de surfaces, selon un mode de réalisation ; 2 3041803 La figure 4 est une vue en perspective d'un modèle de fracture naturelle à titre d'exemple défini à l'aide de surfaces, selon un mode de réalisation ; La figure 5A est une vue en plan d'un modèle de réservoir structuré initial à titre d'exemple référencé dans le diagramme de flux de la figure 1 A ; 5 La figure 5B est une vue de dessus en coupe transversale du modèle de réservoir de la figure 5A prise selon les lignes 5B-5B de la figure 5A ; La figure 6A est une vue en plan du modèle de réservoir de la figure 5A, représenté avec une partie du maillage retiré pour y insérer ensuite un maillage de fracture non structuré en 2,5D ; La figure 6B est une vue de dessus en coupe transversale du modèle de réservoir de la figure 10 6A prise selon les lignes 6B-6B de la figure 6A ; La figure 7A est une vue en plan du modèle de réservoir de la figure 6A, représenté avec un maillage inséré de fracture non structuré en 2,5D dans la couche supérieure ; La figure 7B est une vue de dessus en coupe transversale du modèle de réservoir de la figure 7A prise selon les lignes 7B-7B de la figure 7A ; 15 La figure 7C est une vue de dessus en coupe transversale du modèle de réservoir de la figure 7B, représenté avec un maillage de fracture non structuré extrudé vers le bas ; La figure 8 est un diagramme de flux d'un processus pour la modélisation en 2,5D de certaines fractures selon la figure 1 A et certains exemples à titre d'illustration de la présente invention ; 20 La figure 9 décrit des exemples de configurations de raffinement pour une cellule hexaédrique unique à l'intérieur d'un espace tridimensionnel U-V-W, selon certains aspects à titre d'illustration de la présente invention, qui peuvent être utilisés dans le diagramme de flux de la figure lA ; Les figures 10 et 11 décrivent des cellules de prisme triangulaires représentatives représentant 25 un système de coordonnées U-W bidimensionnel logique et un schéma numérique pour les rebords des prismes, selon certains aspects à titre d'illustration de la présente invention, qui peuvent être utilisés dans le diagramme de flux de la figure lA ; 3 3041803 La figure 12 décrit des configurations de raffinement pour une seule cellule de prisme à l'aide d'un seul paradigme de direction U à l'intérieur du système de coordonnées U-W bidimensionnel logique, selon certains exemples à titre d'illustration de la présente invention, qui peuvent être utilisés dans le diagramme de flux de la figure 1 A ; 5 La figure 13 décrit des configurations d'affinement dans quatre directions à l'aide d'un paradigme à directions U multiples à l'intérieur du système de coordonnées U- W bidimensionnel logique, selon certains aspects à titre d'illustration de la présente invention, qui peuvent être utilisés dans le diagramme de flux de la figure 1 A ; La figure 14 décrit un prisme non triangulaire représentatif (par exemple, un élément 10 prismatique général) représentant le système de coordonnées bidimensionnel logique, selon certaines caractéristiques de la présente invention, qui peuvent être utilisée dans le diagramme de flux de la figure lA ; La figure 15 illustre un raffinement possible d'un élément prismatique général en sous-divisant l'élément en hexaèdres, selon certains exemples donnés à titre d'illustration de la 15 présente invention, qui peuvent être utilisés dans le diagramme de flux de la figure 1 A ; Les figures 16 et 17 sont des diagrammes de flux décrivant un processus pour le raffinement de la cellule à géométrie adaptative, anisotrope et prolongé, selon la figure 1 A et certains exemples à titre d'illustration de la présente invention ; La figure 18 est une représentation schématique en 2D d'un schéma permettant de générer des 20 connexions de maillage réservoir - fracture selon la figure 1C et certains exemples à titre d'illustration de la présente invention ; Les figures 19A à 19C illustrent un maillage de fracture naturelle à titre d'exemple qui a été traitée conformément au procédé des figures 1 A à 1C selon un mode de réalisation ; et La figure 20 est un schéma synoptique fonctionnel d'un système informatique à titre 25 d'exemple selon un mode de réalisation dans lequel le procédé des figures 1 A à 1C peut être mis en oeuvre. 4 3041803 DESCRIPTION DÉTAILLÉE La présente divulgation peut comporter une répétition des références numériques et/ou des lettres dans la présente invention. Cette répétition est faite pour des raisons de simplicité et de clarté et ne suppose pas en elle-même un rapport entre les divers modes de réalisation et/ou 5 les configurations abordés dans la discussion. Par ailleurs, les termes relatifs à l'espace, tels que « en dessous », « dessous », « inférieur », « au-dessus », « supérieur », « tête de puits », « fonds de puits », « amont », « aval » et des expressions analogues, peuvent être utilisés ici afin de faciliter la description pour décrire le rapport d'un élément ou d'une caractéristique avec un autre ou d'autres éléments ou avec une autre ou d'autres caractéristiques comme 10 illustré dans les figures. Les termes relatifs à l'espace sont destinés à englober différentes orientations de l'appareil pendant l'utilisation ou le fonctionnement qui s'ajoutent à l'orientation décrite dans les figures. La présente invention présente de manière détaillée un procédé et un système pour simuler des réservoirs avec de grands systèmes de fracture, d'une géométrie et d'une orientation 15 arbitraires, avec une précision suffisante et dans un temps raisonnable à l'aide de systèmes de simulateurs que l'on trouve dans le commerce en résolvant de manière sélective des données dans des modèles de simulation de taille notablement réduite. Plusieurs modèles de fractures et de réservoirs sont combinés pour former un modèle de réservoir unique qui peut être aisément simulé avec une demande de calcul plus faible et qui a toutefois une résolution 20 suffisante pour capturer des modèles d'écoulement à proximité des fractures. Par conséquent, les modes de réalisation divulgués fournissent un système, un procédé et un produit de programme informatique permettant de générer des maillages de calcul hybrides autour de fractures complexes et discrètes pour les besoins d'une simulation de réservoir. Les figures 1 A à 1C forment ensemble un diagramme de flux de haut niveau en trois parties 25 définissant un procédé à maillage multiple 100 pour modéliser un réservoir d'hydrocarbures fracturé selon un ou plusieurs modes de réalisation. La figure 1 A définit une première partie d'un procédé 100 pour créer un maillage volumétrique du réservoir fracturé qui résout de manière appropriée des fractures souhaitées. Plus particulièrement, la première partie du procédé 100 décrite dans la figure 1A révèle un processus initial en deux phases selon un 30 mode de réalisation dans lequel un algorithme de mailleur (2,5D) en 2,5 dimension, tel qu'un mailleur de stade, est utilisé conjointement avec un algorithme prolongé de raffinement à géométrie adaptative, anisotrope (AGAR) pour produire, à partir d'un modèle de la terre en 5 3041803 grille structurée, un modèle de réservoir raffiné. Cette combinaison peut maintenir les avantages d'une vitesse de calcul et d'une maîtrise de l'orientation d'un mailleur en 2,5D, lorsqu'on le souhaite, tout en faisant appel également à un algorithme AGAR prolongé pour résoudre des réseaux généraux de fractures en 3D. Un réseau de fractures peut être représenté 5 comme des propriétés géométriques et de fracture. Le jeu de fractures peut être divisé en deux groupes : Ceux devant être résolus dans le maillage en dimension 2,5 et ceux devant être résolus uniquement dans le maillage tridimensionnel final. Les fractures permettant la dimension 2,5 comprennent des fractures représentées par des plans verticaux ou des plans sensiblement verticaux par rapport à l'horizon. Les autres fractures peuvent comprendre des 10 fractures représentées par des plans horizontaux ou des plans sensiblement horizontaux. Un maillage en dimension 2,5 est d'abord créé à l'aide du mailleur de stade 2,5 pour résoudre les fractures permettant une dimension en 2,5 dans un espace tridimensionnel (3D). Un algorithme AGAR prolongé résout ensuite anisotropiquement les autres fractures pour produire un modèle de réservoir raffiné. Le raffinement, tel qu'utilisé dans la présente 15 invention, désigne la division des cellules dans l'espace tridimensionnel qui aboutit à un modèle supérieur de résolution de maillage. L'algorithme AGAR prolongé peut être directement appliqué à des topologies comprenant des prismes hexaèdres, triangulaires et d'autres éléments prismatiques généraux pour fournir une gradation de cellule appropriée dans la région proche de la fracture. La figure 1B définit une seconde partie d'un procédé 100 20 pour modéliser les fractures comme un maillage volumétrique qui est indépendant du maillage du réservoir. Enfin, comme illustré à la figure 1C, les deux maillages chevauchants indépendants sont combinés pour former un modèle unique pour permettre une simulation rapide d'un réseau de fractures extrêmement complexe avec des résultats suffisamment précis. Toutefois, en fonction du type de données saisies et de la simulation souhaitée, il est possible 25 d'omettre certaines étapes illustrées du procédé 100. De plus, dans d'autres mises en oeuvre, les étapes définies aux figures 1 A à 1C peuvent se produire en dehors de l'ordre décrit dans les figures. Par exemple, deux blocs présentés dans un ordre successif peuvent, en fait, être exécutés sensiblement concomitamment, ou les blocs peuvent parfois être exécutés dans l'ordre inverse, selon la fonctionnalité concernée. 30 Le diagramme de flux des figures 1 A à 1C illustre une architecture, une fonctionnalité et une utilisation de mises en oeuvre possibles de systèmes, de procédés et de produits de programme informatique selon divers modes de réalisation de la présente invention. Dans un 6 3041803 ou plusieurs modes de réalisation, chaque bloc et/ou combinaison de blocs du diagramme de flux illustré peuvent être mis(e) en oeuvre par des systèmes matériels spéciaux qui réalisent les fonctions ou les actions indiquées, ou des combinaisons de matériels informatiques spéciaux avec des instructions informatiques. 5 En référence à la figure 1A, un modèle de réservoir 200 fournit initialement un détail stratigraphique de subsurface d'un réservoir, comprenant, par exemple, une interprétation géoscientifique, une modélisation de système pétrolier, une analyse géochimique, un maillage stratigraphique, le faciès, un volume net de cellule et une modélisation des propriétés pétrophysiques. Les plates-formes de modélisation de la terre comprennent des plateformes 10 telles que, par exemple, DecisionSpace®, qui est disponible commercialement chez Halliburton Energy Services Inc. de Houston au Texas. Une diversité d'autres plates-formes de modélisation de la terre, peut, cependant, également être utilisée avec la présente description ou à la place de cette dernière. De plus, le modèle de réservoir 200 peut permettre une récupération robuste des données et l'intégration des données historiques ou en temps 15 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. De telles données peuvent comprendre, des données sur la diagraphie, 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 20 surface, les données sur les failles, les données provenant des puits environnants, les données déduites des géostatistiques, etc. Une base de données (non illustrée) peut stocker ces informations. Une plate-forme de base de données est, par ex., le logiciel INSITEO, disponible commercialement chez Halliburton Energy Services Inc., Houston, Texas. Une diversité d'autres plates-formes de base de données, de plates-formes logicielles, et des 25 systèmes associés peut être utilisée pour récupérer, stocker et intégrer les données apparentées aux puits. Les utilisateurs de la simulation de réservoir peuvent entrer la description du réservoir sous plusieurs formes. Par exemple, les utilisateurs peuvent commencer par un modèle structuré de la terre ou par une définition qui peut être convertie facilement en un modèle structuré de la 30 terre. Les modèles initiaux de réservoir contiennent uniquement des colonnes à maillage vertical. Cependant, le maillage structuré de réservoir 200 peut être obtenu par tout procédé adapté et fourni sous une forme adaptée. La figure 2 illustre sous forme graphique la 7 3041803 définition d'un maillage structuré ou extrudé d'un réservoir fourni comme un modèle structuré de la terre (EM). L'ombrage des mailles peut indiquer un contenu en hydrocarbures, une densité, une porosité, une combinaison de celles-ci, ou une caractéristique analogue. Le modèle initial de fractures 220 fournit la définition géométrique du réseau de fractures, 5 avec les distributions de propriétés pétrophysique et mécanique pour les fractures, ce qui est également décrit à la figure 1A. Des modules de système de fractures à titre d'exemple peuvent être mis en oeuvre, à l'aide, par exemple, des plates-formes Fracpro® ou Gohfer®. D'autres plates-formes et systèmes logiciels peuvent, cependant, être utilisés pour mettre en oeuvre le module de systèmes de fracture 220. 10 Le modèle de fracture 220 peut être séparé à l'étape 105 en un premier jeu de fractures 223, généralement des fractures hydrauliques, qui sont bien adaptées pour une modélisation en dimension 2,5 (2.5D) et un second jeu de fractures 227, généralement des fractures naturelles, qui sont moins adaptées pour une modélisation en 2,5D. Les deux jeux de fractures ne s'excluent pas nécessairement l'un l'autre. Autrement dit, certains sous-ensembles de fractures 15 peuvent se trouver dans les deux jeux de fractures. Des illustrations graphiques d'un modèle de fracture hydraulique et un modèle de fracture naturelle courant à titre d'exemple sont présentés respectivement dans les figures 3 et 4, qui peuvent être définies comme des surfaces. L'ombrage des surfaces peut représenter la dimension ou la résistance des fluides des fractures, par exemple. 20 Les modèles de réservoir et de fractures des figures 2 à 4 comprennent respectivement les données initiales pour former le modèle final de réseaux de fractures du procédé à mailles multiples 100. En référence à la figure 1A, le modèle de réservoir structuré 200 peut être modifié pour remplacer une partie du maillage de réservoir structuré existant par un ou plusieurs modèles de fracture, tel qu'indiqué à l'étape 110. Se rappeler qu'à l'étape 105, on 25 détermine un sous-ensemble 223 des fractures 220 (particulièrement des fractures hydrauliques) à résoudre dans le modèle de réservoir structuré 200 par un maillage en 2.5D. Généralement, un sous-ensemble 223 est une simple sélection de fractures qui sont connectées au puits de forage par des chemins de haute conductivité. Le sous-ensemble 223 peut être intégré dans un modèle de réservoir structuré 200 à l'étape 110 à l'aide d'un 30 algorithme de maillage en 2,5D, tel qu'un mailleur de stade, pour résoudre les fractures planes, verticales - une géométrie souvent utilisée pour représenter des fractures hydrauliques - pour créer un maillage de réservoir extrudé non structuré 230. 8 3041803 Par exemple, la figure 5A est une vue en plan d'un modèle de réservoir structuré initial à titre d'exemple 200. La figure 5B est une vue de dessus en coupe transversale du modèle initial de réservoir structuré 200 de la figure 5A prise selon les lignes 5B-5B. On identifie une région 202 associée à une fracture hydraulique unique. Les figures 6A et 6B sont respectivement des 5 vues de dessus en plan et en coupe transversale, du modèle de réservoir structuré 200 représenté à une phase intermédiaire de l'étape 110 avec la région de la fracture 202 retirée. Bien qu'une région rectangulaire 202 ait été identifiée et retirée pour y insérer les fractures 220, toute forme de polygone appropriée peut être utilisée en fonction de ce que les conditions imposent. 10 Les figures 7A et 7B sont respectivement des vues de dessus en plan et en coupe transversale, d'un modèle de réservoir non structuré 230 avec un maillage de fracture bidimensionnelle non structurée 204 projeté initialement sur la surface supérieure d'un maillage autrement structuré. Un schéma de maillage bidimensionnel (2D) peut être utilisé pour ajouter une résolution à proximité du sous-ensemble de la fracture 223 pour créer un maillage de fracture 204. 15 Comme représenté dans la vue de dessus en coupe transversale de la figure 7C, le maillage de fracture 204 est alors projeté sur la surface du maillage inférieur suivant et on répète le processus jusqu'au remplissage de la région 202 par de nouveaux éléments. Vous noterez que dans la vue en coupe de la figure 7C, le maillage de fracture en 2,5 D 204 apparaît structuré. Le maillage du modèle initial de réservoir structuré 200 peut être modifié pour retirer la 20 section 202 à remplacer par le maillage de fracture non structuré 204, et les couches du maillage du modèle initial de réservoir structuré 200 peuvent être utilisées pour guider l'extrusion du maillage de fracture 204. Les propriétés du modèle initial de réservoir structuré 200 sont mises en correspondance avec la nouvelle topologie du maillage du modèle de réservoir non structuré 230. 25 Bien qu'il soit possible d'utiliser n'importe quel algorithme de maillage en 2,5 D approprié pour l'étape 110, dans un ou plusieurs modes de réalisation, un mailleur de stades, qui peut être utilisé pour générer rapidement des grilles non structurées à l'aide d'éléments structurés autour de géométries complexes, peut être utilisé comme suit : Un sous-ensemble de fractures 223 peut être discrétisé dans un plan bidimensionnel par une collection de segments de ligne. 30 La collection des segments de ligne représente le croisement entre le plan bidimensionnel et les fractures géologiques tridimensionnelles. Chaque fracture du sous-ensemble 223 peut être 9 3041803 représentée par une collection de segments linéaires pour faire l'approximation d'une courbure de la fracture. Pour chaque segment de la ligne de fracture dans le plan en 2D, un jeu de stades est généré avec un rayon spécifié provenant du segment de ligne de fracture respectif. Ensuite, des 5 boucles fermées sont générées autour de tous les segments de ligne d'une fracture. Le processus de génération de la boucle fermée autour des segments de ligne de la fracture peut comporter 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 et la mise à l'écart des segments contenus pour chaque segment linéaire dans chaque segment de la ligne de fracture qui sont entièrement 10 contenus par des stades d'autres segments de ligne dans le segment de la ligne de fracture. À la suite de la création de boucles fermées, des éléments de forme peuvent être générés à l'intérieur des boucles fermées du segment linéaire. Par ex., des segments paramétriques peuvent être générés sur une longueur et un rayon de chaque segment linéaire. Des éléments quadrilatères sont formés là où c'est possible à l'intérieur de la région structurée, des 15 polygones étant formés à l'intérieur des régions restantes des boucles fermées. Après avoir généré les éléments de forme, un maillage contraint est généré autour des boucles fermées du jeu de segments de la ligne de fracture, remplissant le reste du plan en 2D. Dans un ou plusieurs modes de réalisation, un algorithme de triangulation de Delaunay est utilisé pour générer le maillage contraint autour des boucles fermées du jeu de segments de la ligne 20 de fracture. Ainsi, le plan en 2D est maintenant entièrement constitué d'éléments de cellule des segments de la ligne de fracture et du maillage contraint. À ce point, le procédé peut extruder chacune des cellules dans le plan en 2D en une troisième 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 représentent 25 une fracture en 3D, alors que les cellules à l'intérieur du maillage représentent des couches rocheuses englobant les fractures. Ainsi, le procédé 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. Un algorithme de maillage de stades en 2,5D fera l'objet d'une discussion plus détaillée ci-après. 30 La figure 8 est un diagramme de flux illustrant un procédé 2500 permettant la modélisation des fractures en 3D en utilisant un mailleur stadimétrique en 2,5D, selon un ou plusieurs 10 3041803 modes de réalisation. Le procédé 2500 commence au niveau du bloc 2502 par recevoir un jeu de fractures 223 (figure 1A) 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 223 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 linéaires pour faire l'approximation d'une courbure de la fracture. Au niveau du bloc 2504, pour chaque segment de la ligne de fracture dans le plan en 2D, le 10 procédé 2500 réalise le bloc 2506 pour générer un jeu de stades avec un rayon spécifié à partir 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. Ce processus peut 15 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, tel que représenté au niveau du bloc 2509 de la figure 8, et le mise à l'écart des segments contenus pour chaque segment linéaire dans chaque segment de la ligne de fracture qui sont entièrement contenus par des stades d'autres segments de ligne dans le segment de la ligne de fracture, tel que représenté au niveau du bloc 2511. 20 Après l'achèvement du bloc 2508 (et des 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 linéaire. Par ex., dans un mode de réalisation, le procédé 2500 génère des segments paramétriques sur une longueur et un rayon de chaque segment linéaire, tel que le montre le bloc 2513. 25 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. 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 la ligne de fracture, remplissant le restant du plan en 2D au niveau du bloc 30 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 11 3041803 ligne de fracture. Ainsi, le plan en 2D est constitué maintenant entièrement d'éléments de cellule des segments de la ligne de fracture et du maillage contraint. À ce stade, le procédé 2500 peut extruder chacune des cellules dans le plan en 2D en une troisième dimension permettant de créer une ou plusieurs couches de cellules en 3D. Les 5 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 les 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 10 démontre le bloc 2516. En revenant à la figure 1A, après la création du maillage du modèle non structuré 230 pour capturer le champ de pression à proximité du premier sous-ensemble de fractures 223, à l'étape 115, un algorithme prolongé de mailleur de raffinement à géométrie adaptative anisotrope (AGAR) peut être utilisé pour capturer le flux à proximité du second sous- 15 ensemble de fractures restantes 227. Autrement dit, un algorithme non structuré AGAR peut être utilisé pour raffiner le maillage de modèle non structuré 230 pour résoudre les fractures restantes 227, créant ainsi un maillage de modèle de réservoir raffiné 235. Dans un ou plusieurs modes de réalisation, le raffinement peut être grossier pour toutes les fractures restantes 227. Souvent, une partie importante des fractures du sous-ensemble 227 20 contribuera de manière moins significative au flux du système de réservoir et produira des gradients de pression moins sévères. Lorsque l'on s'attend à ce qu'il en soit ainsi, ces zones ne nécessitent pas une résolution de maillage aussi importante que les fractures 223 gérées principalement avec le schéma de maillage en 2,5 D. En conséquence, il est possible d'obtenir une réduction dans le décompte d'éléments par rapport à une utilisation dans l'ensemble du 25 réservoir de l'algorithme AGAR non structuré. Cependant, l'algorithme AGAR non structuré peut être affiné pour fournir une résolution plus élevée à proximité de fractures particulières dont on s'attend à ce qu'elles contribuent de manière plus significative au flux du réservoir. Un algorithme prolongé de mailleur AGAR adapté pour être utilisé à l'étape 115 selon un ou plusieurs modes de réalisation va maintenant être décrit. En général, les algorithmes de 30 mailleur AGAR raffinent les rebords des cellules jusqu'à « n » fois et ne divisent aucun rebord (c'est-à-dire en divisant la cellule) qui produirait un rebord plus court qu'un 12 3041803 TargetSize. Tel que défini ici, le TargetSize représente la taille de maillage 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 du délai souhaité jusqu'à la solution. Il peut être déterminé si une cellule donnée 5 doit ou ne doit pas être raffinée et dans quelle direction les raffinements doivent se produire. 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). 10 Pour résumer la mise en oeuvre AGAR de l'étape 115, il faut prendre en compte la direction U (une direction parmi trois directions cardinales dans un espace en 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 U peut pointer dans une direction légèrement différente qu'un autre 15 rebord U. Comme il sera décrit plus en détail ci-dessous, en considérant la direction U, l'algorithme AGAR prolongé analysera toutes les cellules dans un modèle et raffinera les cellules proches de la fracture dans la direction U si tous les rebords U sont plus longs que C x TargetSize et qu'au moins une des cinq 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, 20 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 comporte 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 une variable scalaire spécifique 25 à une mise en oeuvre ou fournie par un utilisateur, C, l'étape 115 (figure 1A) ne divisera aucun rebord (donc, des cellules) qui est plus court que C x TargetSize. Dans ce procédé illustratif, C est choisi comme .Ni2, mais peut être d'autres variables tels que, par exemple, 4/3 ou toute valeur finie. Ce procédé est répété pour chaque direction dans l'espace 3D. Tel qu'il est utilisé ici, le terme « noeud suspendu » est généralement utilisé pour désigner un 30 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 ne peut être un hexaèdre (tel que représenté dans les exemples de la figure 9), un 13 3041803 prisme (tel que représenté dans les exemples des figures 10 à 13), 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 5 peut être omise pour les règles 3 et 4 si une adhérence stricte à la directive du noeud suspendu est adoptée. La figure 9 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 figure 9 illustre comment, pour un hexaèdre, certains modes de réalisation 10 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émontrent les sous-étapes 308, 310 et 314 de la figure 9. Une fois que les directions de raffinement souhaitées ont été déterminées pour une cellule, l'étape 115 de l'algorithme AGAR prolongé (figure 1A) réalise le raffinement (c'est-à-dire, en divisant les rebords) pour cette cellule. La 15 figure 9 illustre huit configurations de raffinements possibles pour une cellule unique à l'intérieur d'un espace en 3D U- V- W 302. L'étape de l'algorithme AGAR prolongé 115 (figure 1A) 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 (sous-étapes, 308, 310, et 314 illustrées dans la figure 9), et ensuite, si c'est le cas, en divisant les 20 rebords dans les directions tel que représenté pour créer des raffinements de cellule 306, 312 et/ou 316. Les huit configurations possibles de raffinement de cellules 304, 306, 312 et 314 vont du raffinement dans des directions de non-raffinement (hexaèdre 304) au raffinement dans l'ensemble des 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 25 fracture, la résolution du maillage des images ainsi obtenues sera grandement améliorée. La figure 10 illustre une cellule de prisme triangulaire représentative 404 et un système de coordonnées logique à 2D U-W 402 et un schéma numérique pour les rebords du prisme 404. Certains modes de réalisation comprennent deux alternatives pour étendre le concept de raffinement susmentionné aux prismes. L'un consiste en l'option « direction U unique » 30 illustrée dans la figure 10. 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 14 3041803 de la figure 10, 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 « directions U multiples », illustrée dans la figure 11. 5 La figure 11 illustre une cellule de prisme triangulaire 504 et un système de coordonnées logique à 2D U-W 502 et un schéma numérique pour les rebords du prisme 504. En utilisant l'option des directions U multiples 500 de la figure 11, le raffinement se ferait dans des paires de direction U. Ainsi, le raffinement pourrait se faire dans un sous-ensemble des directions { U1 U2, U1U3, U2 U3, W} 502 représentées dans la figure 11. De manière générale, dans les 10 règles appliquées par certains modes de réalisation, on n'autorise pas l'ensemble des 16 permutations dans les quatre directions (par exemple, W, U1, U2, et U3) au cours de l'étape de raffinement AGAR 115 (figure 1A). Dans l'exemple de la figure 11, les directions U multiples sont considérées en paires. Par ex., le raffinement dans U1U2 raffinera les rebords 1, 2, 4, 5 du prisme 504. 15 La figure 12 illustre un exemple de raffinement AGAR 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 figure 11, les rebords de la cellule, Ci, sont étiquetés comme {e,j} pour j = 1 ... 9. Considérons maintenant une fonction Booléenne abstraite sur les rebords, Test ({e,j}). Par ex., la cellule peut être raffinée 20 dans : U si n'importe quel rebord à partir du jeu {e,,i, e,,2, ei,6} satisfait Test({e,j}) (voire sous-étape 606 de la figure 12), et W si n'importe quel rebord à partir du jeu {ei,7, e,,9} satisfait Test({e,,}) (voir la sous-étape 610 de la figure 12). La condition susmentionnée est appelée la condition principale pour la gradation. Nous allons maintenant définir des conditions opposées. Pour les quadrilatères dans le prisme 604, la 25 règle est la même que pour les hexaèdres (par exemple, comme il a été présenté ci-dessus en référence à la figure 9). Si, disons deux rebords opposés sur la face, exactement un de ces rebords satisfait Test({e, j}), alors la direction qui divise ces deux rebords sera également raffinée. Par exemple, les rebords {e,,4, e,,1} sont des rebords opposés sur une face de quadrilatère. Si le rebord {e,,4} satisfait Test({e,,j}) mais pas le rebord {e,,1}, alors cette 30 condition serait satisfaite, et la direction W serait ajoutée à l'ensemble de raffinement. De même, si le rebord {e,,7} satisfait Test({e,,j}) mais pas le rebord {e,,8}, la direction W serait ajoutée au jeu de raffinement. En ce qui concerne la face triangulaire, si exactement un 15 3041803 satisfait Test({e, j}), alors U est ajouté au jeu de raffinement. Par exemple, si le rebord {ei.4} satisfait Test({ e,,}) mais pas les rebords {ei,5, e,,6}, la direction W serait ajoutée au jeu de raffinement. Pour les règles de gradation de l'algorithme prolongé AGAR, Test({ei j}) est défini comme 5 vrai si et seulement si {e11} a été défini au moins deux fois. Pour les règles de croisement de l'algorithme prolongé AGAR, Test({eid}) est défini comme vrai si le rebord échelonné {e,1} croise une fracture d'intérêt. La figure 12 illustre les combinaisons possibles de raffinements 608 et 612 utilisant le paradigme à direction U unique 600. La figure 13 illustre la façon dont on peut, dans une alternative de la direction U unique 10 présentée ci-dessus en référence à la figure 12, obtenir un paradigme de multiples directions U 700. En particulier, la figure 13 illustre des raffinements 708, 712, 716, 720 et 722 dans les quatre dimensions/directions logiques 702 pour un paradigme de multiples directions U 700 d'un prisme 704. En utilisant le paradigme des multiples directions U 700, la cellule sera raffinée dans : U/U2 et U1U3 si n'importe quel rebord à partir du jeu {e,,i, e,,4} satisfait 15 Test({ev }) (voir sous-étapes 706 et 710 de la figure 13) ; U1U2 et U2U3 si n'importe quel rebord à partir du jeu {e,,2, eo} satisfait Test({e,,,}) (voir sous-étapes 706 et 714 de la figure 13) ; U1U3 et U2U3 si n'importe quel rebord à partir du jeu{e,,3, e,,6} satisfait Test({e, j}) (voir sous-étapes 710 et 714 de la figure 13) ;et W si n'importe quel rebord à partir du jeu eo} satisfait Test({e,'}) (voir la sous-étape 718 de la figure 13). 20 En plus de ces conditions primaires, il existe des conditions opposées. Pour des quadrilatères, les deux paradigmes U sont très similaires. Si, disons deux rebords opposés sur une face, exactement un de ces rebords satisfait Test({e,j}), alors la direction qui divise ces deux rebords sera également raffinée. Si le rebord {e,,7} satisfait Test({e,,,,}) mais pas le rebord {e0}, la direction W serait ajoutée au jeu de raffinement. Dans un exemple, en ce qui 25 concerne la règle de la face triangulaire, si exactement un rebord satisfait Test({e,j}), alors U1 U2 et U/ U3 peuvent être ajoutés au jeu de raffinement. Dans ce paradigme, aucune condition opposée n'est ajoutée pour les faces triangulaires. Cependant, des règles alternatives peuvent être déterminées et testées dans le cadre de paradigmes au-delà des paradigmes 600 et 700 illustrés respectivement dans les figures 13 et 13. 30 La figure 14 décrit un prisme non triangulaire représentatif 1004 (par exemple, un élément prismatique général) représentant le système de coordonnées bidimensionnel logique 1002. 16 3041803 La figure 15 illustre un raffinement possible d'un élément prismatique général 1104 en sous-divisant l'élément en hexaèdres. Comme décrit ci-dessus, l'étape de l'algorithme prolongé AGAR 115 (figure 1A) peut être utilisée pour mettre en oeuvre des procédés qui prolongent le schéma AGAR basé sur un hexaèdre pour inclure également des prismes triangulaires, tels 5 que ceux illustrés dans les figures 10 à 13. Pour des maillages triangulaires et quadrilatéraux extrudés, cette combinaison est suffisante. En supposant une topologie extrudée cohérente, les paradigmes présentés ci-dessus en référence aux figures 10 à 13 permettront d'obtenir un maillage valide si certains des éléments sont dégénérés. Les figures 14 et 15 illustrent la façon dont certains modes de réalisation peuvent prolonger directement le paradigme de 10 direction unique U en des prismes non triangulaires, en utilisant un algorithme de sous- division du point médian. Pour des prismes non triangulaires, tels que le prisme non triangulaire 1004 de la figure 14, l'anisotropie peut être maintenue dans l'espace U-W 1002. Comme on peut le voir dans le raffinement 1104 de la figure 14, en utilisant cet algorithme de raffinement, les éléments prismatiques non hexaédriques 1004 sont subdivisés en hexaèdres. 15 Puisque les hexaèdres ont trois dimensions logiques, par rapport aux deux illustrées dans les figures 14 et 15 pour les éléments prismatiques non hexaédriques, la possibilité d'un raffinement anisotrope augmente après le raffinement de l'élément d'origine. Pour généraliser davantage, un algorithme de sous-division du point médian peut être appliqué pour toutes les cellules polyédriques convexes, y compris les cellules non extrudées. 20 Il est à noter que cela changerait l'algorithme de raffinement pour des prismes triangulaires de celui présenté ci-dessus en référence aux figures 10 à 13. Cependant, pour les cellules non extrudées, un espace U dimensionnel unique peut être utilisé en é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 25 procurera encore un avantage important sur le raffinement entièrement isotrope. En outre, un raffinement anisotrope sur une cellule régulière avoisinant une cellule irrégulière ne viole pas l'exigence de façonnage sur le maillage. La face polygonale de la cellule voisine irrégulière est simplement divisée. Compte tenu de la description détaillée précédente, un procédé illustratif pour l'étape 115 de 30 la figure 1 A va maintenant être décrit en référence aux diagrammes de flux des figures 16 et 17. En faisant référence à la figure 16, le procédé 2100 commence au niveau du bloc 2102. Dans un exemple, le bloc 2102 comprend des données provenant du modèle de réservoir 17 3041803 extrudé non structuré 230. Au niveau du bloc 2104, une cellule dans le modèle de réservoir extrudé non structuré 230 peut être analysée pour déterminer si la cellule dans un maillage doit être raffinée de manière anisotrope (c'est-à-dire raffinée dans n'importe quelle direction particulière dans un espace en 3D). S'il est déterminé que la cellule en cours d'analyse ne 5 nécessite pas de raffinement, la commande passe au bloc 2109. Au niveau du bloc 2109, il est déterminé s'il y a d'autres cellules à analyser. S'il n'y a pas d'autres cellules à analyser, la commande passe au bloc 2110. Autrement, s'il est déterminé qu'il y a une autre cellule restante à analyser, la commande retourne au bloc 2104 où cette cellule sera ensuite analysée pour déterminer si un raffinement peut être nécessaire. 10 Si, au niveau du bloc 2104, il est déterminé que la cellule nécessite un raffinement, l'algorithme se déplace au bloc 2108 où cette cellule est raffinée. Tel que décrit précédemment, les règles de gradation et de croisement peuvent être employées pour déterminer si une cellule doit être raffinée. Une fois que la cellule est raffinée dans toutes les directions déterminées comme étant nécessaires, l'algorithme passe au bloc 2109 où une 15 détermination a lieu pour savoir s'il y a d'autres cellules à analyser. Le processus se répète jusqu'à ce que toutes les cellules aient été analysées. Lorsque l'analyse de toutes les cellules est achevée, au niveau du bloc 2110, le réseau de fractures à l'intérieur du modèle de réservoir est résolu en utilisant les cellules raffinées. Le modèle de réservoir raffiné ainsi obtenu possèdera un niveau de résolution de maillage plus élevé dans les zones encerclant les 20 fractures (par exemple, des zones proches de la fracture et les fractures), et une résolution de maillage plus faible dans les zones non proches de la fracture. La figure 17 est un diagramme de flux d'un procédé de raffinement illustratif 2200 utilisé dans le bloc 2104 de la figure 16. L'emplacement des fractures à l'intérieur du système de fractures est connu sur la base des données de fracture 220 (figure 1A). En utilisant ces 25 données, ainsi que les règles de gradation et de croisement décrites précédemment, ces cellules dans les régions proches de la fracture du modèle qui nécessitent un raffinement peuvent être localisées. Pour les cellules nécessitant un raffinement, tous les rebords (et, ainsi, les cellules) peuvent être coupés en deux dans chaque direction de raffinement. Les cellules nouvellement créées seront ensuite analysées à nouveau dans l'itération suivante, 30 augmentant de ce fait finalement la résolution de maillage. Le maillage à haute résolution ainsi obtenu procure des résultats numériques plus précis dans les régions proches de la 18 3041803 fracture/de la fracture. Le C x TargetSize (par exemple, où C = i2) peut être déterminé comme décrit précédemment. Au niveau du bloc 2104(i), il peut être déterminé 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 5 TargetSize. Si aucun des rebords n'est plus long que C x TargetSize, aucun raffinement n'est requis et la commande passe au bloc 2109. Cependant, si un ou plusieurs rebords sont plus longs que C x TargetSize, les cinq règles illustrées ci-dessus peuvent alors s'appliquer. À ce stade, au niveau des blocs 2104 (ii) à (vi), il peut être déterminé si une ou plusieurs des conditions suivantes sont satisfaites : ii) un rebord directionnel échelonné croise une fracture ; 10 iii) pour toute paire de rebords opposée dans la direction 3D, exactement l'un des deux rebords échelonnés croise une fracture ; iv) un rebord directionnel comporte deux ou plusieurs noeuds « suspendus » ; v) pour toute paire de rebords opposée dans la direction 3D, exactement l'un des deux rebords comporte deux ou plusieurs noeuds suspendus ; ou vi) la cellule non échelonnée croise une fracture mais aucun rebord non échelonné de la cellule ne 15 croise une fracture. Au niveau du bloc 2104(ii), il peut être déterminé si la version échelonnée de la cellule comporte un rebord dans la direction analysée qui croise une fracture. Au niveau du bloc 2104(iii), une paire de rebords opposée dans une direction donnée peut être analysée pour déterminer si seulement un des rebords échelonnés croise une fracture. Au niveau du bloc 20 2104(ii), il peut être déterminé si un rebord dans la direction analysée comporte deux noeuds suspendus ou davantage. Dans la figure 17, l'élément 2209 est fourni comme illustration de ce principe. L'élément 2209 est décrit ci-dessous. Comme représenté à la figure 17, l'élément 2209 comprend deux faces de cellule adjacentes 2210 et 2220 (représentées en 2D pour des raisons de simplicité) ayant deux noeuds 2212 et 25 2214 positionnés entre elles. Dans l'exemple de l'élément 2209, le premier raffinement aurait divisé la cellule 2210 en deux en utilisant le noeud 2212. Le second raffinement aurait divisé l'une des nouvelles cellules en deux en utilisant le noeud 2214. Par conséquent, les noeuds 2212 et 2214 sont suspendus par rapport à la cellule 2220 parce que, bien que la cellule 2210 ait été raffinée (dans une itération précédente), la cellule 2220 n'a pas été raffinée. Puisque la 30 cellule 2220 comporte au moins deux noeuds suspendus, elle sera également raffinée dans des itérations ultérieures. 19 3041803 Au niveau du bloc 2104(v) du procédé 2200, des paires de rebords opposées dans la direction analysée peuvent être examinées pour déterminer si seulement un des rebords comporte deux noeuds suspendus ou davantage. Au niveau du bloc 2104(iv), il peut être déterminé si une cellule non échelonnée croise une fracture mais aucun rebord non échelonné de la cellule ne 5 croise une fracture. Si la décision est positive pour n'importe laquelle des cinq règles illustrées ci-dessus en référence aux blocs 2104(i-v), le procédé 2200 passe la commande au bloc 2108 où cette cellule est raffinée. Autrement, si la décision est négative en ce qui concerne l'ensemble des cinq règles, la commande passe au bloc 2109. Ce procédé continuera de manière itérative 10 jusqu'à ce que chaque cellule dans le modèle ait été analysée, tel que présenté ci-dessus en rapport avec la figure 16. En faisant référence maintenant à la figure 1B, outre la création d'un maillage de réservoir raffiné 235, le procédé 100 comprend la création d'un maillage pour le système de fractures. Le maillage de la fracture peut être représenté comme un maillage volumétrique ou un 15 maillage de surface avec une propriété d'épaisseur (ouverture). Dans certains cas, les données de fracture peuvent être fournies initialement comme une surface discrétisée. Si cette surface est discrétisée avec une résolution suffisante pour une simulation, elle peut être utilisée directement, conjointement à une valeur d'ouverture fournie. Dans d'autres cas, les surfaces devront être maillées ou raffinées pour atteindre la résolution souhaitée. 20 En conséquence, un maillage de surface de fracture approprié 240 peut être construit à partir des données de fracture à l'étape 120. Pour des fractures 220 déjà définies comme un maillage de surface, ce maillage peut être utilisé tel quel, ou des algorithmes de raffinement ou de grossissement peuvent être utilisés à l'étape 120 pour créer la résolution souhaitée. À l'étape 125, chaque élément du maillage de fracture de la surface de la fracture 240 peut alors être 25 épaissi par une valeur d'ouverture de la fracture pour former un maillage volumétrique de la fracture 245. Bien que généralement le volume du réseau de fractures dans le maillage de fracture 245 ne soit pas une fraction suffisante de n'importe quelle cellule de réservoir dans le maillage de réservoir raffiné 235, à ce stade, on peut déterminer quel volume de fracture est situé à 30 l'intérieur de n'importe quelle cellule de réservoir. Si on le souhaite, les propriétés de réservoir du maillage de réservoir raffiné 235 peuvent être modifiées à l'étape 130 pour tenir 20 3041803 compte du volume qui est présent dans le maillage volumétrique de la fracture 245. Les maillages de réservoir et de fracture raffinés 235, 245 se chevaucheront, comme présenté ci-après, le volume de pores dans un maillage de réservoir raffiné 235 peut être réduit pour tenir compte du volume maintenant représenté par le maillage chevauchant du volume de fracture 5 245. En d'autres termes, les volumes des blocs matriciels encerclant les fractures dans le maillage de réservoir raffiné 235 peuvent être modifiés à l'étape 130 pour maintenir le volume de pores exact en retirant le volume de pores des volumes de contrôle de la matrice qui connectent des fractures, obtenant ainsi un maillage de réservoir intermédiaire 250. Le volume de pores retiré dépend du nombre et de la taille des fractures auxquelles le volume de 10 contrôle de la matrice est connecté. Le maillage de la fracture 245 et le maillage du réservoir intermédiaire 250 définissent deux maillages volumétriques discrets qui se chevauchent dans l'espace. En référence maintenant à la figure 1C, selon un ou plusieurs modes de réalisation, le procédé 100 comprend la combinaison du maillage du volume de la fracture 245 avec le maillage du réservoir 15 intermédiaire 250 pour former un modèle de réservoir final 260. À l'étape 135, trois jeux distincts de connexions entre le maillage du volume de la fracture 245 et le maillage du réservoir intermédiaire 250 sont formés : Des connexions matrice-matrice, des connexions fracture-fracture et des connexions matrice-fracture. Des connexions d'un premier jeu, des connexions matrice-matrice, peuvent être formées à 20 l'étape 135a entre des paires de cellules adjacentes dans le maillage de réservoir intermédiaire 250. Ces connexions peuvent être calculées d'une manière classique en utilisant des approximations de flux à deux points (TPFA), des approximations de flux multipoints (MPFA) ou d'autres approximations. Le schéma TPFA utilise le potentiel de pression aux centres de deux cellules adjacentes pour permettre de calculer le flux entre lesdites cellules. 25 Le schéma MPFA utilise de manière complémentaire le potentiel de pression provenant de toutes les cellules avoisinantes de deux cellules adjacentes pour calculer le flux entre des paires de cellules. Avec des grilles orthogonales (la grille peut être orthogonale par rapport au champ de perméabilité - souvent désigné sous le terme « k-orthogonal »), le schéma TPFA peut s'avérer suffisamment précis. Cependant pour des grilles non orthogonales, comme c'est 30 souvent le cas avec des grilles non structurées, le schéma TPFA perd en précision et le schéma MPFA peut être préféré pour maintenir la précision. 21 3041803 Un second jeu de connexions, connexions fracture-fracture, peut être formé à l'étape 135b entre des paires de cellules adjacentes dans le maillage volumique de la fracture 245. Dans un ou plusieurs modes de réalisation, une technique de discrétisation à différence finie de volume de contrôle utilisant une approximation de flux à deux points peut être utilisée. À la 5 fois pour le flux multiphase et unique, la compensation de matériau pour chaque volume de contrôle requiert de connaître les volumes de contrôle avoisinants-une liste de connectivité- et le débit associée à chaque connexion, qui peut être déterminée en utilisant la transmissibilité pour un support poreux fracturé, tel que décrit ci-dessous. Pour n'importe quelle forme du volume de contrôle et dimension du problème, le débit peut 10 être indiqué sous la forme : Q12 = T122(P2 -P1) , (Éq. 1) où p est la pression, Q12 est le débit de la cellule 1 à la cellule 2, T12 est la partie géométrique de la transmissibilité, et X, représente la mobilité du fluide en utilisant les informations en amont. Dans le cas d'un flux multiphase, différents débits, pressions et mobilités s'appliquent 15 pour chaque phase. L'élément de mobilité de la transmissibilité, qui est différent pour chaque phase, peut être calculé de la façon habituelle. Le composant géométrique de la transmissibilité, T12, qui est le même pour chaque phase, peut être exprimé sous la forme : I 2 avec D, T a1a2 a 12 a ± a (Éq. 2) 20 où, A, est la zone de l'interface entre deux volumes de contrôle (en utilisant les informations de CV' où CV, désigne le volume de contrôle ith), k, étant la perméabilité de CV' D, est la distance entre le centroïde de l'interface et le centroïde de CV' n, est le vecteur d'unité normal à l'interface à l'intérieur de CV' et f, est le vecteur d'unité le long de la direction de la ligne joignant le centroïde du volume de contrôle avec le centroïde de l'interface. L'ensemble des 25 informations géométriques nécessaires pour calculer T12 est défini dans le domaine de la grille. Le calcul de la transmissibilité ci-dessus peut être utilisé directement soit pour des problèmes en 2D soit pour des problèmes en 3D. Dans une configuration en 2D, l'interface est un segment, alors que dans une configuration en 3D, l'interface est un polygone. 22 3041803 La connexion numérique entre deux fractures peut être accomplie par l'utilisation d'un volume de contrôle intermédiaire (CV0). L'objet de ce volume de contrôle intermédiaire est de permettre une redirection du flux et une variation de l'épaisseur entre deux fractures. Cependant, dans un ou plusieurs modes de réalisation, il peut être avantageux de ne pas 5 introduire CV0, ni ses inconnues associées, dans le modèle numérique, car un volume de contrôle intermédiaire peut introduire des problèmes numériques en raison de la petite taille de CV0 par rapport aux volumes de contrôle des deux cellules, CV/ et CV2. La transmissibilité entre CV/ et CV2 peut être exprimée, tout en tenant compte implicitement du volume de contrôle intermédiaire CV0, comme suit : 10 T = 10T 02 (Éq. 3) T '12 T10 ± TO2 où T10 est la transmissibilité entre CV/ et CV0 et T02 est la transmissibilité entre CV0 et CV2. En conséquence, T12 est simplement une moyenne harmonique de T10 et T02, qui est appropriée pour des cellules en série. La définition de T12 impose de connaître T10 et T02, ce qui rend nécessaire une définition 15 géométrique du volume de contrôle intermédiaire. Afin d'éviter l'introduction de cette définition géométrique, qui peut être complexe dans certaines configurations, il est possible de faire la simplification suivante. La taille du volume de contrôle intermédiaire étant généralement faible en comparaison des volumes de contrôles adjacents, et le volume de contrôle intermédiaire étant généralement caractérisé par une perméabilité semblable à celle 20 des volumes de contrôle de fracture environnants, on peut supposer que : DI » Do a1 REVENDICATIONS1. Procédé pour modéliser des fractures géologiques tridimensionnelles à l'intérieur d'un réservoir, comprenant : la génération d'un maillage volumétrique de matrice dudit réservoir ; la génération d'un maillage volumétrique de fracture desdites fractures géologiques ; l'établissement de connexions matrice-matrice à l'intérieur dudit maillage de matrice ; l'établissement de connexions fracture-fracture à l'intérieur dudit maillage de fracture ; et l'établissement de connexions matrice-fracture entre lesdits maillages de matrice et de 10 volume de fracture. 2. Procédé selon la revendication 1, dans lequel : ledit maillage de matrice est orthogonal ; et ledit établissement de connexions matrice-matrice comprend en outre la discrétisation 15 à l'aide d'approximations de flux à deux points du potentiel de pression entre des cellules adjacentes dans ledit maillage de matrice. 3. Procédé selon la revendication 1, dans lequel : ledit maillage de matrice comprend des cellules non orthogonales ; et 20 ledit établissement de connexions matrice-matrice comprend en outre la discrétisation à l'aide d'approximations de flux multipoints du potentiel de pression de toutes les cellules avoisinantes des cellules adjacentes dans ledit maillage de matrice. 4. Procédé selon la revendication 1, dans lequel ledit établissement de connexions 25 fracture-fracture comprend en outre : l'approximation d'une transmissibilité entre des cellules adjacentes dans ledit maillage de fracture ; la détermination d'un débit entre lesdites cellules adjacentes à l'aide de ladite transmissibilité ; et 30 la discrétisation à l'aide d'approximations de flux à deux points ou multipoints dudit débit. 34 3041803 5. Procédé selon la revendication 1, dans lequel ledit établissement de connexions matrice-fracture comprend en outre : la détermination de faces se croisant géographiquement entre des cellules dudit maillage de matrice et des cellules dudit maillage de fracture ; et 5 pour chaque face sécante, le calcul d'une valeur de transmissibilité de matrice au niveau de la face dudit maillage de matrice, le calcul d'une valeur de transmissibilité de fracture au niveau de la face dudit maillage de fracture, et le calcul d'une transmissibilité à travers la face entre ledit maillage de matrice et ledit maillage de fracture comme une moyenne harmonique desdites valeurs de transmissibilité de matrice et de fracture. 6. Procédé selon la revendication 5, dans lequel ledit établissement de connexions matrice-fracture comprend en outre : le calcul de ladite valeur de transmissibilité de matrice à l'aide d'une approximation de flux à deux points ; et le calcul de ladite valeur de transmissibilité de fracture à l'aide d'une approximation de flux à deux points. 7. Procédé selon la revendication 5, dans lequel : ladite génération d'un maillage volumétrique de fracture comprend l'extrusion d'une face d'origine de chaque cellule à l'intérieur d'un maillage de surface de fracture par une largeur d'ouverture correspondante ; et ledit établissement de connexions matrice-fracture comprend en outre, pour chaque face, le calcul desdites valeurs de transmissibilité de fracture et de matrice comme fonctions du croisement de la face d'origine correspondante dudit maillage de surface de fracture avec la face dudit maillage de matrice. 8. Procédé selon la revendication 7, dans lequel ledit établissement de connexions matrice-fracture comprend en outre : pour chaque face, le calcul desdites valeurs de transmissibilité de fracture et de 30 matrice comme fonctions du croisement de la face d'origine correspondante dudit maillage de surface de fracture avec la face dudit maillage de matrice. 9. Procédé selon la revendication 8, dans lequel ledit établissement de connexions matrice-fracture comprend en outre : 35 3041803 la détermination d'une distance représentative entre le centroïde d'une cellule dans ledit maillage de matrice et la face sécante de ladite cellule ; et le calcul de ladite valeur de transmissibilité de matrice pour ladite face de ladite cellule en fonction de ladite distance représentative. 5 10. Procédé selon la revendication 1, dans lequel ladite génération dudit maillage volumétrique de matrice comprend en outre : la réception d'une spécification de réservoir ; l'identification, sur la base de la spécification de réservoir, d'un jeu de fractures 10 comprenant des fractures autorisant une dimension 2,5 et d'autres fractures ; la génération d'un modèle de réservoir non structuré comprenant un maillage d'extrusion qui modélise les fractures autorisant une dimension 2,5 en trois dimensions ; le raffinement anisotrope d'une ou plusieurs cellules dans le modèle de réservoir non structuré correspondant aux autres fractures ; 15 la résolution d'un réseau de fractures dans le modèle de réservoir non structuré à l'aide de cellules raffinées ; et la génération d'un modèle de réservoir raffiné à l'aide du réseau de fractures. 11. Procédé selon la revendication 1, dans lequel : 20 les fractures autorisant une dimension 2,5 ont une géométrie qui a été discrétisée dans un plan bidimensionnel par une collection de segments de ligne ; et ladite génération dudit maillage volumétrique de matrice comprend en outre, pour chaque segment de ligne associé à chaque fracture dans les fractures autorisant une 2,5 D, la génération d'un ensemble de stades à un rayon spécifié à partir 25 du segment de ligne, la génération de boucles fermées autour du segment de ligne, et la génération d'éléments de forme à l'intérieur des boucles fermées du segment de ligne, et la génération du maillage en tant que maillage contraint autour de boucles fermées des fractures autorisant une 2,5 D pour remplir un espace restant du plan bidimensionnel. 12. Procédé selon la revendication 10, dans lequel ledit raffinement anisotrope de cellules comprend en outre : 36 3041803 l'identification d'une direction dans les trois dimensions dans lesquelles les cellules doivent être raffinées ; et la division d'un rebord des cellules, le rebord étant dans la direction dans les trois dimensions. 5 13. Procédé selon la revendication 1, dans lequel : le procédé pour modéliser des fractures géologiques tridimensionnelles à l'intérieur du réservoir est exécuté par un système informatique. 10 14. Support de stockage lisible par un ordinateur sur lequel sont stockées des instructions exécutables qui, si elles sont exécutées par un dispositif informatique, entraîne la réalisation, par le dispositif informatique, d'opérations, les instructions comprenant : la génération d'un maillage volumétrique de matrice dudit réservoir ; la génération d'un maillage volumétrique de fracture desdites fractures géologiques ; 15 l'établissement de connexions matrice-matrice à l'intérieur dudit maillage de matrice ; l'établissement de connexions fracture-fracture à l'intérieur dudit maillage de fracture ; et l'établissement de connexions matrice-fracture entre lesdits maillages de matrice et de volume de fracture. 20 15. Support de stockage lisible par un ordinateur selon la revendication 14, dans lequel : ledit maillage de matrice est orthogonal ; et lesdites instructions comprennent en outre la discrétisation à l'aide d'approximations de flux à deux points du potentiel de pression entre des cellules adjacentes dans ledit maillage 25 de matrice pour établir lesdites connexions matrice-matrice. 16. Support de stockage lisible par un ordinateur selon la revendication 14, dans lequel : ledit maillage de matrice comprend des cellules non orthogonales ; et lesdites instructions comprennent en outre la discrétisation à l'aide d'approximations 30 de flux multipoints du potentiel de pression de toutes les cellules avoisinantes des cellules adjacentes dans ledit maillage de matrice pour établir lesdites connexions matrice-matrice. 17. Support de stockage lisible par un ordinateur selon la revendication 14, dans lequel, pour établir lesdites connexions fracture-fracture, lesdites instructions comprennent en outre : 37 3041803 l'approximation d'une transmissibilité entre des cellules adjacentes dans ledit maillage de fracture ; la détermination d'un débit entre lesdites cellules adjacentes à l'aide de ladite transmissibilité ; et 5 la discrétisation à l'aide d'approximations de flux à deux points ou multipoints dudit débit. 18. Support de stockage lisible par un ordinateur selon la revendication 14, dans lequel, pour établir lesdites connexions matrice-fracture, lesdites instructions comprennent en outre : 10 la détermination de faces se croisant géographiquement entre des cellules dudit maillage de matrice et des cellules dudit maillage de fracture ; et pour chaque face sécante, le calcul d'une valeur de transmissibilité de matrice au niveau de la face dudit maillage de matrice, le calcul d'une valeur de transmissibilité de fracture au niveau de la face dudit maillage de fracture, et le calcul d'une transmissibilité à 15 travers la face entre ledit maillage de matrice et ledit maillage de fracture comme moyenne harmonique desdites valeurs de transmissibilité de matrice et de fracture. 19. Support de stockage lisible par un ordinateur selon la revendication 18, dans lequel, pour établir lesdites connexions matrice-fracture, lesdites instructions comprennent en outre : 20 le calcul de ladite valeur de transmissibilité de matrice à l'aide d'une approximation de flux à deux points ; et le calcul de ladite valeur de transmissibilité de fracture à l'aide d'une approximation de flux à deux points. 25 20. Support de stockage lisible par un ordinateur selon la revendication 18, dans lequel : pour générer un maillage volumétrique de fracture, lesdites instructions comprennent en outre l'extrusion d'une face d'origine de chaque cellule à l'intérieur d'un maillage de surface de fracture par une largeur d'ouverture correspondante ; et pour établir lesdites connexions matrice-fracture, lesdites instructions comprennent en 30 outre, pour chaque face, le calcul desdites valeurs de transmissibilité de fracture et de matrice comme fonctions du croisement de la face d'origine correspondante dudit maillage de surface de fracture avec la face dudit maillage de matrice. 38 3041803 21. Support de stockage lisible par un ordinateur selon la revendication 20, dans lequel, pour établir lesdites connexions matrice-fracture, lesdites instructions comprennent en outre : pour chaque face, le calcul desdites valeurs de transmissibilité de fracture et de matrice comme fonctions du croisement de la face d'origine correspondante dudit maillage de 5 surface de fracture avec la face dudit maillage de matrice. 22. Support de stockage lisible par un ordinateur selon la revendication 21, dans lequel, pour établir lesdites connexions matrice-fracture, lesdites instructions comprennent en outre : la détermination d'une distance représentative entre le centroïde d'une cellule dans 10 ledit maillage de matrice et la face sécante de ladite cellule ; et le calcul de ladite valeur de transmissibilité de matrice pour ladite face de ladite cellule en fonction de ladite distance représentative. 23. Support de stockage lisible par un ordinateur selon la revendication 14, dans lequel, 15 pour générer ledit maillage volumétrique de matrice, lesdites instructions comprennent en outre : la réception d'une spécification de réservoir ; l'identification, sur la base de la spécification de réservoir, d'un jeu de fractures comprenant des fractures autorisant une dimension 2,5 et d'autres fractures ; 20 la génération d'un modèle de réservoir non structuré comprenant un maillage d'extrusion qui modélise les fractures autorisant une dimension 2,5 en trois dimensions ; le raffinement anisotrope d'une ou plusieurs cellules dans le modèle de réservoir non structuré correspondant aux autres fractures ; la résolution d'un réseau de fractures dans le modèle de réservoir non structuré à l'aide 25 des cellules raffinées ; et la génération d'un modèle de réservoir raffiné à l'aide du réseau de fractures. 24. Support de stockage lisible par un ordinateur selon la revendication 14, dans lequel : les fractures autorisant une dimension 2,5 ont une géométrie qui a été discrétisée dans 30 un plan bidimensionnel par une collection de segments de ligne ; et pour générer ledit maillage volumétrique de matrice, lesdites instructions comprennent en outre, pour chaque segment de ligne associé à chaque fracture dans les fractures autorisant une 2,5 D, la génération d'un ensemble de stades à un rayon spécifié à partir 39 3041803 du segment de ligne, la génération de boucles fermées autour du segment de ligne, et la génération d'éléments de forme à l'intérieur des boucles fermées du segment de ligne, et la génération du maillage en tant que maillage contraint autour de boucles 5 fermées des fractures autorisant une 2,5 D pour remplir un espace restant du plan bidimensionnel. 25. Support de stockage lisible par un ordinateur selon la revendication 23, dans lequel, pour raffiner de manière anisotrope des cellules, lesdites instructions comprennent en outre : 10 l'identification d'une direction dans les trois dimensions dans lesquelles les cellules doivent être raffinées ; et la division d'un rebord des cellules, le rebord étant dans la direction dans les trois dimensions. 15 26. Système pour modéliser des fractures géologiques tridimensionnelles à l'intérieur d'un réservoir, comprenant : un processeur ; et et une mémoire sur laquelle sont stockées des instructions, qui, si elles sont exécutées par le processeur, entraîne la réalisation par le processeur d'opérations comprenant 20 la génération d'un maillage volumétrique de matrice dudit réservoir, la génération d'un maillage volumétrique de fracture desdites fractures géologiques, l'établissement de connexions matrice-matrice à l'intérieur dudit maillage de matrice, 25 l'établissement de connexions fracture-fracture à l'intérieur dudit maillage de fracture, et l'établissement de connexions matrice-fracture entre lesdits maillages de matrice et de volume de fracture. 30 27. Système selon la revendication 26, dans lequel : ledit maillage de matrice est orthogonal ; et lesdites instructions comprennent en outre la discrétisation à l'aide d'approximations de flux à deux points du potentiel de pression entre des cellules adjacentes dans ledit maillage de matrice pour établir lesdites connexions matrice-matrice. 40 3041803 28. Système selon la revendication 26, dans lequel : ledit maillage de matrice comprend des cellules non orthogonales ; et lesdites instructions comprennent en outre la discrétisation à l'aide d'approximations 5 de flux multipoints du potentiel de pression de toutes les cellules avoisinantes de cellules adjacentes dans ledit maillage de matrice pour établir lesdites connexions matrice-matrice. 29. Système selon la revendication 26, dans lequel, pour établir lesdites connexions fracture-fracture, lesdites instructions comprennent en outre : 10 l'approximation d'une transmissibilité entre des cellules adjacentes dans ledit maillage de fracture ; la détermination d'un débit entre lesdites cellules adjacentes à l'aide de ladite transmissibilité ; et la discrétisation à l'aide d'approximations de flux à deux points ou multipoints dudit 15 débit. 30. Système selon la revendication 26, dans lequel, pour établir lesdites connexions matrice-fracture, lesdites instructions comprennent en outre : la détermination de faces se croisant géographiquement entre des cellules dudit 20 maillage de matrice et des cellules dudit maillage de fracture ; et pour chaque face sécante, le calcul d'une valeur de transmissibilité de matrice au niveau de la face dudit maillage de matrice, le calcul d'une valeur de transmissibilité de fracture au niveau de la face dudit maillage de fracture, et le calcul d'une transmissibilité dans la face entre ledit maillage de matrice et ledit maillage de fracture comme moyenne 25 harmonique desdites valeurs de transmissibilité de matrice et de fracture. 31. Système selon la revendication 30, dans lequel, pour établir lesdites connexions matrice-fracture, lesdites instructions comprennent en outre : le calcul de ladite valeur de transmissibilité de matrice à l'aide d'une approximation de 30 flux à deux points ; et le calcul de ladite valeur de transmissibilité de fracture à l'aide d'une approximation de flux à deux points. 32. Système selon la revendication 30, dans lequel : 41 3041803 pour générer un maillage volumétrique de fracture, lesdites instructions comprennent en outre l'extrusion d'une face d'origine de chaque cellule à l'intérieur d'un maillage de surface de fracture par une largeur d'ouverture correspondante ; et pour établir lesdites connexions matrice-fracture, lesdites instructions comprennent en 5 outre, pour chaque face, le calcul desdites valeurs de transmissibilité de fracture et de matrice comme fonctions du croisement de la face d'origine correspondante dudit maillage de surface de fracture avec la face dudit maillage de matrice. 33. Système selon la revendication 32, dans lequel, pour établir lesdites connexions 10 matrice-fracture, lesdites instructions comprennent en outre : pour chaque face, le calcul desdites valeurs de transmissibilité de fracture et de matrice comme fonctions du croisement de la face d'origine correspondante dudit maillage de surface de fracture avec la face dudit maillage de matrice. 15 34. Système selon la revendication 33, dans lequel, pour établir lesdites connexions matrice-fracture, lesdites instructions comprennent en outre : la détermination d'une distance représentative entre le centroïde d'une cellule dans ledit maillage de matrice et la face sécante de ladite cellule ; et le calcul de ladite valeur de transmissibilité de matrice pour ladite face de ladite 20 cellule en fonction de ladite distance représentative. 35. Système selon la revendication 26, dans lequel, pour générer ledit maillage volumétrique de matrice, lesdites instructions comprennent en outre : la réception d'une spécification de réservoir ; 25 l'identification, sur la base de la spécification de réservoir, d'un jeu de fractures comprenant des fractures autorisant une dimension 2,5 et d'autres fractures ; la génération d'un modèle de réservoir non structuré comprenant un maillage d'extrusion qui modélise les fractures autorisant une dimension 2,5 en trois dimensions ; le raffinement anisotrope d'une ou plusieurs cellules dans le modèle de réservoir non 30 structuré correspondant aux autres fractures ; la résolution d'un réseau de fractures dans le modèle de réservoir non structuré à l'aide des cellules raffinées ; et la génération d'un modèle de réservoir raffiné à l'aide du réseau de fractures. 42 3041803 36. Système selon la revendication 26, dans lequel : les fractures autorisant une dimension 2,5 ont une géométrie qui a été discrétisée dans un plan bidimensionnel par une collection de segments de ligne ; et pour générer ledit maillage volumétrique de matrice, lesdites instructions 5 comprennent en outre, pour chaque segment de ligne associé à chaque fracture dans les fractures autorisant une 2,5 D, la génération d'un ensemble de stades à un rayon spécifié à partir du segment de ligne, la génération de boucles fermées autour du segment de ligne, et la génération d'éléments de forme à l'intérieur des boucles fermées du segment de 10 ligne, et la génération du maillage en tant que maillage contraint autour de boucles fermées des fractures autorisant une 2,5 D pour remplir un espace restant du plan bidimensionnel. 15 37. Système selon la revendication 35, dans lequel, pour raffiner de manière anisotrope des cellules, lesdites instructions comprennent en outre : l'identification d'une direction dans les trois dimensions dans lesquelles les cellules doivent être raffinées ; et la division d'un rebord des cellules, le rebord étant dans la direction dans les trois 20 dimensions. 43