FR 3048300 A1 20170901 FR 1750680 A 20170127 Représéntation Hybride Géocellulaire En 3D De Sous-Ensembles De Réseau De Fracture Naturelle Choisis DOMAINE DE LA DESCRIPTION La présente divulgation concerne généralement la modélisation des réservoirs d'hydrocarbures et, plus spécifiquement, un système de simulation des réservoirs qui représente un sous-ensemble d'un réseau de fracture naturelle (« RFN ») directement dans la simulation, alors que l'approximation du restant du RFN est réalisée par une formulation multicontinuum. CONTEXTE 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 du réservoir avant la production. Par contre, des modèles de simulation du réservoir peuvent être créés pour simuler le flux de fluide dans le réservoir au cours de sa durée de vie de production. Dans un flux de travail intégré, des modèles statiques, représentés par des grilles structurées ou non-structurées, peuvent également être utilisés dans la simulation des réservoirs, servant souvent de point initial pour la simulation du flux au cours du temps. Avec les modèles de simulation de réservoir, la modélisation des fractures à l'intérieur d'un réservoir peut représenter un défi, 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 décrites comme 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. Au cours de leur histoire géologique, certaines fractures naturelles peuvent également être bouchées par des dépôts à l'intérieur des fissures. La bonne modélisation des fractures est importante étant donné que des propriétés des fractures telles que la distribution spatiale, l'ouverture, la longueur, la hauteur, la conductivité et la 3048300 2 connectivité affectent de façon importante le flux des fluides du réservoir vers le puits de forage. Les fractures naturelles peuvent avoir un impact énorme sur la production d'un réservoir conventionnel ou non-conventionnel. Étant donné que les fractures naturelles peuvent influencer l'efficacité d'une fracturation hydraulique, aussi bien que la productivité du puits et la connectivité entre les puits, la simulation de ces scénarios est très importante pour les équipes d'actifs qui tentent de prendre les bonnes décisions lors de l'exploration et la production des réservoirs de schiste. Des exemples de ces décisions sont l'estimation des réserves, la planification du champ (pour une récupération primaire, secondaire ou tertiaire), l'optimisation de la complétion, l'optimisation de l'espacement entre les puits, l'optimisation des installations en surface, etc. Souvent ces décisions sont complexes, nécessitant à la fois des considérations financières et techniques. Par conséquent, les simulateurs de réservoir ont modélisé des fractures naturelles pendant très longtemps (décennies). Cependant, la pratique la plus répandue actuellement consiste toujours à utiliser une porosité double ou plus généralement une approximation du continuum double pour représenter toutes les fractures naturelles dans le modèle de simulation. Certains produits de simulation utilisent des grilles non-structurées pour le maillage de toutes les fractures naturelles directement dans la grille de simulation. Cependant, ces grilles ont tendance à avoir un nombre de cellules extrêmement grand qui fait que la simulation est trop grande ou trop gourmande en ressources pour la plupart des systèmes informatiques. BRÈVE DESCRIPTION DES DESSINS La figure 1 illustre un organigramme d'un système de simulation hybride 100 selon certains modes de réalisation illustratifs de la présente divulgation ; La figure 2 est un organigramme d'un procédé permettant de générer une grille géocellulaire hybride en 3D, selon certains procédés illustratifs de la présente divulgation ; 3048300 3 La figure 3A est une illustration d'une interface utilisateur graphique d'une trace stéréonet pour la sélection de sous-ensembles de fracture naturelle, selon certains modes de réalisation illustratifs de la présente divulgation ; La figure 3B illustre un schéma de filtrage alternatif pour sélectionner un sous-ensemble de fracture, selon encore un autre mode de réalisation illustratif de la présente divulgation ; La figure 3C est une visualisation illustrative en 3D d'un RFN avec deux sous-ensembles de fracture ; La figure 4 illustre un nombre de fractures raffinées intégrées dans une grille géocellulaire en 3D, selon un exemple de la présente divulgation ; La figure 5 est un organigramme d'un procédé alternatif permettant de générer une grille géocellulaire hybride en 3D, selon certains procédés illustratifs de la présente divulgation ; La figure 6 est un procédé permettant de valider une grille géocellulaire hybride en 3D, selon certains procédés illustratifs de la présente divulgation ; La figure 7 illustre une comparaison des résultats de simulation par rapport à la production gazière cumulative ; et La figure 8 illustre une comparaison des champs de pression dans la matrice à 2000 jours, illustrée pour le réservoir en entier (à gauche) et une vue détaillée de la zone (à droite) mise en évidence par la boîte de Type 2 hybride (la ligne du bas de la première colonne). DESCRIPTION DES MODES DE REALISATION ILLUSTRATIFS Les modes de réalisation illustratifs et les méthodologies apparentées de la présente divulgation sont décrits ci-dessous de la façon dont ils pourraient être utilisés dans un simulateur de réservoir qui génère une représentation géocellulaire hybride en 3D des RFN. Dans un souci de clarté, toutes les caractéristiques d'une implémentation ou d'une 3048300 4 méthodologie ne sont pas décrites dans cette spécification. Il sera, bien sûr, apprécié que dans le développement d'un quelconque mode de réalisation réel, de nombreuses décisions spécifiques à une implémentation doivent être prises afin d'atteindre les objectifs spécifiques des développeurs, tels que la conformité à des contraintes apparentées au système ou aux considérations monétaires, qui varieront d'une concrétisation à une autre. En outre, il sera apprécié qu'un tel effort de développement puisse être complexe et chronophage, mais serait néanmoins une entreprise de routine pour les hommes de métier qui bénéficient de cette divulgation. D'autres aspects et avantages des divers modes de réalisation et des méthodologies apparentées de la divulgation deviendront évidents à la lumière de la description et des figures suivantes. Comme il sera décrit ci-dessous, les procédés et modes de réalisation illustratifs de la présente divulgation modifient une grille géocellulaire en 3D pour représenter un sous-ensemble de RFN directement dans la simulation, alors que l'approximation du restant du RFN est réalisée par une formulation multi-continuum. En règle générale, les modes de réalisation divulgués ici décrivent deux procédés alternatifs par lesquels un utilisateur peut sélectionner, de façon interactive, lequel des sous-ensembles de RFN qui sera directement simulé et dont l'approximation sera réalisée, pour ainsi créer un modèle hybride avec à la fois une modélisation à double continuum et une modélisation directe des RFN basée sur une grille non-structurée. Dans un premier procédé illustratif, une grille géocellulaire, ayant des propriétés à double continuum, est d'abord générée à partir d'un sous-ensemble de fracture choisi d'un RFN. Un autre sous-ensemble de fracture choisi du RFN est ensuite directement maillé, ou intégré, dans la grille. Dans un deuxième procédé illustratif, un sous-ensemble de fracture choisi du RFN est directement intégré dans la grille, par lequel des propriétés à double continuum sont ensuite calculées sur la grille. Ainsi, un utilisateur peut intégrer certaines des fractures naturelles pour une utilisation directe dans la simulation, alors qu'on peut faire l'approximation des autres fractures sous forme de modèles à double porosité ou à double continuum. Par conséquent, 3048300 une approche hybride est rendue possible par laquelle les approximations inhérentes dans le double continuum sont équilibrées contre les avantages de la modélisation directe des fractures naturelles pour donner une simulation équilibrée dans une fourchette de temps raisonnable, même pour des grands RFN (par ex., >100 000 fractures). Par conséquent, le simulateur a une performance améliorée de beaucoup. La production réelle du modèle hybride n'est généralement pas onéreuse en termes de temps écoulé même dans un poste de travail d'ingénierie type (par ex. HP Z800). Cependant, il est possible de se trouver dans des situations dans lesquelles la grille et/ou la taille des RFN est tellement grande que les calculs nécessaires pour créer un modèle hybride deviennent, eux-mêmes, un goulot d'étranglement. Les calculs de mise à l'échelle représentent les parties les plus gourmandes en informatique de ce procédé, et ces calculs peuvent être réalisés en parallèle en distribuant les calculs sur de multiples processus fonctionnant en parallèle sur un ordinateur ou plusieurs ordinateurs, dans un environnement informatique à haute performance (HPC) ou dans un environnement de nuage (avec de multiples ordinateurs virtuels). En outre, les calculs de mise à l'échelle avec une cellule de la grille sont indépendants des calculs avec une autre grille. Ceci rend très pratique la production de modèles hybrides mêmes avec de grands problèmes. Une fois le modèle hybride créé, la simulation elle-même peut être exécutée dans un procédé unique ou de multiples procédés en parallèle, utilisant un logiciel disponible dans le commerce. Dans une description plus généralisée des procédés alternatifs mentionnés ci-dessus, un réseau de fracture ayant une pluralité de cellules, à l'intérieur de celui-ci, est créé en utilisant les données sur les réservoirs. Les RFN sont composés de fractures qui peuvent être très hétérogènes par rapport à leur emplacement spatial, leur forme et leurs dimensions physiques et leurs propriétés. Plusieurs applications logicielles, dont beaucoup sont disponibles dans le commerce (par ex., DecisionSpace® RFN), peuvent être utilisées pour générer des représentations discrètes de ces réseaux. Généralement, dans ces applications, les algorithmes stochastiques sont utilisés en raison de leur 3048300 6 degré élevé d'incertitude associé à la caractérisation des fractures souterraines. Quelquefois, des études de caractérisation de la fracture réalisées sur des excroissances de la surface sont utilisées comme analogues à partir desquels des données sont fournies aux applications logicielles de génération de RFN. Le logiciel de génération de RFN rend possible la génération et la visualisation de multiples réalisations de réseaux de fracture dans le même logiciel utilisé pour visualiser des données géologiques et géophysiques (G&G) souterraines, généralement dans le cadre d'un modèle terrestre partagé. La visualisation elle-même est utile pour améliorer la compréhension géologique dans un réservoir naturellement fractionné, mais une utilisation plus percutante des RFN générés consiste à les utiliser dans des modèles de simulation de réservoir afin de quantifier leur effet sur l'écoulement de fluide et, par conséquent, la productivité du puits dans des réservoirs naturellement fractionnés. Une fois généré, un premier sous-ensemble de fracture est sélectionné parmi le réseau de fractures. Ensuite, un second sous-ensemble de fracture, différent du premier sous-ensemble de fracture, est sélectionné parmi le réseau de fractures. Les cellules de la grille géocellulaire en 3D proche du premier sous-ensemble de fracture sont raffinées, et le second sous-ensemble de fracture est mis à l'échelle. Ensuite, le premier sous-ensemble de fracture et le second sous-ensemble de fracture mis à l'échelle sont intégrés dans la grille géocellulaire en 3D, générant ainsi une grille géocellulaire en 3D hybride. Il est à noter que les grilles géocellulaires en 3D décrites ici peuvent être générées en utilisant des grilles structurées ou non-structurées. Dans un procédé alternatif, le premier et le second sous-ensembles de fracture sont encore une fois sélectionnés comme dans le procédé précédent. Ici, cependant, une grille géocellulaire en 3D à haute résolution est générée à travers un raffinement ou la création directe de la grille à haute résolution. Ensuite, le premier et le second sous-ensembles de fracture sont intégrés dans la grille géocellulaire en 3D. Le premier sous-ensemble de fracture et le second sous-ensemble de fracture sont ensuite mis à l'échelle dans la grille géocellulaire en 3D. 3048300 7 Dans les modes de réalisation illustratifs décrits ici, diverses techniques peuvent être utilisées pour raffiner et intégrer (également appelé « résolution ») les fractures. Le raffinement, tel qu'il est utilisé ici, décrit la division des cellules dans l'espace en 3D, un procédé qui en fin de compte, entraîne la formation d'une grille géocellulaire en 3D à résolution de maillage plus élevée. Le réseau de fracture est ensuite résolu (ou intégré), à l'intérieur de la grille géocellulaire en 3D à haute résolution. Le terme intégration/résolution, tel qu'il est utilisé ici, décrit la création d'un modèle terrestre, utilisant le maillage à haute résolution, dans lequel les zones des fractures et les zones proches des fractures dans le modèle ont une résolution de maillage plus élevée en comparaison à d'autres zones du modèle. Étant donné que le maillage sera être utilisé pour une simulation numérique, en général, les cellules plus petites (en raison du 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 la formule mathématique. La mise à l'échelle décrit le fait de réaliser une approximation de l'effet des fractures naturelles sur l'écoulement à travers une cellule de grille d'une grille d'entrée en calculant les propriétés réelles de la fracture, avec celles qui quantifient l'écoulement de fluide dans la matrice-fracture. Les propriétés calculées, notamment les perméabilités directionnelles de la fracture, la porosité de la fracture et certaines mesures de la fonction de transfert (telles que, par ex., le facteur de forme ou la transmissibilité de l'échange de fluide matrice-fracture), sont appliquées à la cellule de la grille à double continuum qui correspond à la cellule de la grille d'entrée dans la matrice sur laquelle la mise à l'échelle est réalisée. Une diversité de techniques de raffinement peut être utilisée avec les modes de réalisation de la présente divulgation. Une technique de raffinement illustrative est appelée algorithme anisotrope géométriquement adaptatif (« AGAR »). AGAR est un procédé permettant d'intégrer des réseaux de fracture dans des grilles géocellulaires en 3D existantes, qui 3048300 8 existent sous forme d'une grille structurée et des valeurs de propriétés physiques. Le réseau de fracture est sous la forme d'un maillage en surface (avec ou sans collecteur) dans un espace tridimensionnel (« 3D »). La grille structurée est ensuite raffinée de façon anisotrope pour intégrer les fractures et permettre un classement par grade approprié de la cellule au niveau de la région proche de la fracture. En règle générale, la technique AGAR examine des cellules afin de déterminer les cellules qui nécessitent un raffinement, de raffiner ces cellules dans les directions appropriées (par conséquent, l'anisotropie), et ensuite de terminer une fois qu'aucun raffinement de cellules n'est nécessaire. En outre, étant donné que l'AGAR applique une approche anisotrope lors de la manipulation des cellules, le nombre de données et ses exigences informatiques associées sont grandement réduites. Par conséquent, les procédés divulgués sont plus rapides et plus efficaces que les procédés conventionnels qui appliquent une analyse de force brute à toutes les cellules, nécessitant des calculs plus gourmands en ressources et en données. En outre, en raison de la relation de la grille géocellulaire en 3D générée par rapport à la grille géocellulaire en 3D originelle, l'interpolation des caractéristiques est minimisée et la grille géocellulaire en 3D générée respecte les valeurs de la propriété originelle. Ainsi, en sus des autres avantages décrits ici, les 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 aussi les besoins en calculs informatiques. La figure 1 illustre un organigramme d'un système de simulation hybride 100 selon certains modes de réalisation illustratifs de la présente divulgation. Le système de simulation hybride 100 comprend au moins 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 (I/O) 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 3048300 9 instructions logicielles stockées à l'intérieur du moteur hybride 110 conformément aux modes de réalisation illustratifs et les procédés décrits ici, peuvent être stockées dans un stockage 104 ou sur un quelconque autre support de stockage non-transitoire lisible par ordinateur. Même s'il n'est pas explicitement illustré dans la FIG. 1, il sera reconnu que le système de simulation hybride 100 peut être connecté à un ou plusieurs réseaux publics ou privés à travers une ou plusieurs connexions de réseaux appropriées. Il sera également reconnu que les instructions logicielles comprenant le moteur hybride 110 peuvent également être téléchargées dans un stockage 104 à partir d'un support de stockage approprié (par ex., une mémoire portable/un disque dur, un CD-ROM, etc.) par des procédés sans fil ou sur fil. En outre, les hommes du métier apprécieront que cette divulgation puisse être pratiquée avec une diversité de configurations de système informatique, comprenant des dispositifs portatifs, des systèmes multiprocesseurs, et des dispositifs électroniques programmables par l'utilisateur ou à base de microprocesseur, des mini-ordinateurs, 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 description. La description peut être pratiquée dans des environnements informatiques en réseau dans lesquels les tâches sont réalisées par des dispositifs télécommandés qui sont reliés à travers un réseau de télécommunication. Dans un environnement informatique en réseau, les modules de programme peuvent se trouver à la fois sur un support de stockage informatique local et distant, 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. En se rapportant toujours à la Fig. 1, dans certains modes de réalisation illustratifs, le moteur hybride 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 des propriétés pétrophysiques et mécaniques 3048300 pour les fractures, comme il sera compris par les hommes du métier bénéficiant de cette divulgation. Un exemple de module de système de fracture peut être, par ex., les plateformes Fracpro® ou Gohfer®. En se rapportant toujours au mode de réalisation illustratif de la Fig. 1, un moteur hybride 110 comprend également un module de modélisation de la terre 114 pour s'intégrer avec les données de la fracture contenues dans le module du système de fracture 112 afin de procurer une visualisation stratigraphique souterraine des fractures et réservoirs comprenant, par ex., une interprétation géoscience, une modélisation du système pétrolifère, une analyse géochimique, un grillage stratigraphique, le faciès, le volume net de la cellule et la 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é. Des plateformes illustratives de modélisation de la terre comprennent le DecisionSpace®, qui est une application logicielle disponible commercialement et développée et commercialisée par Landmark Graphics Corporation, Houston, Texas. Cependant, les hommes du métier qui bénéficient de cette divulgation réalisent qu'une diversité d'autres plateformes de modélisation de la terre ou d'autres jeux d'applications logicielles peuvent également être utilisés avec la présente divulgation. Le moteur de simulation hybride 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, par ex., des données de diagraphie, des trajectoires de puits, des données sur les propriétés pétrophysiques des roches, des données sur les propriétés mécaniques des roches, des données de surface, des données de faille, des données provenant des puits environnants, des données déduites des géostatistiques, etc. La base de données (non illustrée) qui stocke ces informations peut se trouver 3048300 11 dans le module du système de fracture 112 ou à un emplacement distant. Une plateforme illustrative de base de données est, par ex., le logiciel INSITE®, disponible commercialement chez Halliburton Energy Services Inc., Houston, Texas. Les hommes du métier qui bénéficient de cette divulgation réalisent qu'il existe une diversité de plateformes logicielles et des systèmes associés permettant de récupérer, de stocker et d'intégrer les données apparentées aux puits, tel qu'il est décrit ici. En outre, un moteur hybride 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. En tant que tel, 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 plateformes de flux de travail qui peut être utilisée dans cet objectif. 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/intégrer 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é de sélectionner des sous-ensembles de fracture pour une simulation directe et une mise à l'échelle afin de réduire les exigences informatiques sur le système de traitement. La figure 2 est un organigramme d'un procédé 200 permettant de générer une grille géocellulaire hybride en 3D, selon certains procédés illustratifs de la présente divulgation. Au niveau du bloc 202, le système de simulation hybride 100 est initialisé afin de générer un modèle de réservoir terrestre (c.-à-d., une grille géocellulaire en 3D) ayant un RFN à l'intérieur de celle-ci. Au cours de l'initialisation, le moteur hybride 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., 3048300 12 une grille géocellulaire et des données de fracture. La géométrie pour les fractures saisies peut être fournie de plusieurs façons, et peut comprendre des descriptions de surfaces analytiques ou discrètes. Les propriétés de la fracture nécessaires pour un quelconque procédé illustratif sont définies par le calcul de la propriété à double continuum spécifique utilisé et la simulation de réservoir spécifique qui est exécutée, comme il sera compris par les hommes du métier qui bénéficient de cette divulgation. Au niveau du bloc 204, le système de simulation hybride 100 invite l'utilisateur à sélectionner des sous-ensembles de fracture dans le RFN. Ici, les sous-ensembles comprendront un premier sous-ensemble de fracture et un second sous-ensemble de fracture différent du premier. Ce filtrage défini par l'utilisateur peut être réalisé de différentes façons. Par ex., la sélection peut être réalisée à l'aide d'un stéréonet, d'un attribut de fracture ou d'une quelconque autre visualisation appropriée. La figure 3A est une illustration d'une interface utilisateur graphique d'une tracée stéréonet pour la sélection de sous-ensemble de fracture naturelle, selon certains modes de réalisation illustratif de la présente divulgation. Comme le démontre la figure 3A, un utilisateur peut sélectionner un jeu de fractures de façon interactive. L'utilisateur peut utiliser un dispositif d'entrée graphique, telle qu'une souris ou un stylo, pour sélectionner une zone contiguë du diagramme stéréonet ou pour sélectionner une ou plusieurs barres, l'une à côté de l'autre, dans l'histogramme. L'utilisateur peut ensuite répéter de telles actions avec une action modifiée, de sorte qu'il/qu'elle peut créer une collection de toutes les zones sélectionnées en ajoutant ou en retranchant des zones ou des barres de la collection. Une fois créée, la collection dans son intégralité représente les jeux de fractures sélectionnés par l'utilisateur, c.-à-d., toutes les fractures qui sont représentées dans la zone sélectionnée du stéréonet ou dans les barres sélectionnées de l'histogramme comprennent le jeu de fractures sélectionné par l'utilisateur, parmi toutes les fractures existantes dans les RFN qui sont considérées, pour d'autres actions. La figure 3B illustre un schéma de filtrage alternatif pour sélectionner un sous-ensemble de fracture, selon encore un autre mode de 3048300 13 réalisation illustratif de la présente divulgation. Dans cet exemple, un utilisateur peut sélectionner les sous-ensembles de fracture pour la modélisation directe en utilisant une diversité d'attributs de factures 30 comprenant, par ex., l'azimut, la longueur, la porosité, la perméabilité, la conductivité, l'ouverture ou d'autres attributs définis par l'utilisateur. Des opérations booléennes 32, avec des fourchettes 34, peuvent également être utilisées pour définir les attributs, comme on le souhaite. Dans encore d'autres techniques de sélection, un utilisateur peut filtrer des sous-ensembles de fracture basé sur un facteur de connectivité 36 représentant une distance connectée de chaque fracture à une perforation de puits. Il peut également y avoir une fourchette de facteurs de connectivité acceptables. Par ex., un facteur de connectivité peut être calculé en se basant sur une matrice de proximité qui représenterait un graphique de toutes les fractures. Le facteur de connectivité serait basé sur la « distance » connectée d'une fracture à la perforation du puits. Par ex., un utilisateur pourrait sélectionner comme sous-ensemble toutes les fractures qui sont connectées à une perforation de puits par deux degrés de séparation ou moins. L'utilisateur peut utiliser une valeur de facteur de connectivité comme un seuil contre lequel toutes les fractures doivent être filtrées de sorte que seules les fractures qui ont un facteur de connectivité supérieure au seuil seront directement modélisées ou mises à l'échelle. La figure 3C est une visualisation illustrative en 3D d'un RFN avec deux sous-ensembles de fracture. Une fois le premier et le second sous-ensembles de fracture sélectionnés au niveau du bloc 204, les cellules de la grille proches de/dans le premier sous-ensemble de fracture sont raffinées au niveau du bloc 206. Ici, le moteur hybride 110 analyse toutes les cellules dans/proches du premier sous-ensemble de fracture afin de déterminer les fractures qui nécessitent un raffinement dans une quelconque direction dans l'espace 3D. Comme il a été précédemment mentionné, une diversité de techniques de raffinement peut être utilisée, l'une d'entre elles étant la technique AGAR. En utilisant AGAR, les bords de la cellule sont raffinés jusqu'à « n » fois et aucun bord n'est divisé (et, ainsi, division de la cellule) si la 3048300 14 division produirait un bord plus court que le « targetSize ». « targetSize » représente la taille de maillage désiré, ou la longueur du bord, nécessaire pour 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 grilles géocellulaires en 3D, les procédés doivent déterminer si une cellule doit être raffinée et dans quelle direction ces cellules doivent l'être. Dans certains modes de réalisation illustratifs, il peut y avoir deux types de règles permettant de déterminer si une 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). Dans un exemple généralisé pour résumer le procédé AGAR, 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 « U » peut pointer dans une direction légèrement différente qu'un autre rebord « U ». En considérant la direction U, le moteur hybride 110 analysera toutes les cellules dans le modèle et raffinera les cellules proches des fractures 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 mis à l'échelle croise une fracture ; 2) pour toute paire de rebords opposée dans la direction U, exactement l'un des deux rebords mis à l'échelle 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 mis à l'échelle comportent deux ou plusieurs noeuds suspendus ; ou 5) la cellule non-mis à l'échelle 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 C spécifique d'une implémentation ou fournie par un utilisateur le système de simulation hybride 100 ne divisera aucun rebord (donc, des cellules) qui est plus court que C x targetSize. Dans 3048300 ce procédé illustratif, C est choisi comme V2, mais pourrait ê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. Un « 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 d'un hexaèdre. 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 modes de réalisation illustratifs, l'exigence de la longueur du rebord peut être omise pour les règles 3 et 4 si un respect strict de la directive du noeud suspendu unique est adopté. Une fois que les directions de raffinement ont été déterminées pour une cellule, un moteur hybride 110 réalise le raffinement (c'est à dire, division des bords) pour cette cellule. Les cellules peuvent être raffinées en une diversité de configurations à l'intérieur de l'espace 3D u-v-w. Le moteur hybride 110, réalise de façon itérative les procédés décrits ici en déterminant d'abord si les cellules proches de la fracture doivent être raffinées dans les directions U, V ou W, et ensuite, si tel est le cas, la division des bords dans les directions nécessaires. Les configurations possibles du raffinement des cellules vont d'un raffinement dans aucune direction de raffinement à un raffinement dans toutes les directions de raffinement. É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. La figure 4 illustre un nombre de fractures raffinées 40 intégrées dans une grille géocellulaire en 3D 42, selon un exemple de la présente divulgation. Au niveau du bloc 208, le système de simulation hybride 100 met à l'échelle le second sous-ensemble de fracture. La mise à l'échelle entraîne les propriétés du continuum de fracture, comprenant les propriétés de transfert entre la matrice et les continuums de fracture. Plusieurs techniques existent pour calculer les approximations à double continuum, 3048300 16 aussi bien que des approximations multi-continuums. Par ex., une telle technique d'approximation de continuum est décrite dans Oda, M., 1985, Permeability Tensor for Discontinuous Rock Masses, Geotechnique Vol. 35, pp 483. Par conséquent, au niveau du bloc 208, la grille géocellulaire en 3D contient maintenant les propriétés à double continuum (par ex., une porosité double). Au niveau du bloc 210, le système de simulation hybride 100 intègre maintenant le premier sous-ensemble de fracture et le second sous-ensemble de fracture mis à l'échelle dans la grille géocellulaire en 3D, générant ainsi une grille géocellulaire en 3D hybride. Ici, les propriétés (par ex., la perméabilité) des cellules représentant le premier sous-ensemble de fracture sont ajustées pour tenir compte des propriétés de ces fractures (c.- à-d., directement modélisées), alors que l'approximation des cellules restantes représentant le second sous-ensemble de fracture utilisant la matrice et les propriétés du continuum de fracture, est décrite ci-dessous. Dans certains modes de réalisation illustratifs, les cellules du premier sous-ensemble de fracture comprendront seulement les propriétés de la fracture. Par conséquent, la sortie obtenue est une grille géocellulaire en 3D qui possède une résolution de maillage plus élevée dans les zones entourant le premier sous-ensemble de fracture (c.-à-d., des zones proches de la fracture et les fractures), et une résolution de maillage plus faible dans les zones du second sous-ensemble de fracture. La figure 5 est un organigramme d'un procédé 500 permettant de générer une grille géocellulaire hybride en 3D, selon certains procédés illustratifs de la présente divulgation. Au niveau du bloc 502, 504 et 506, le RFN est généré, le premier et le second sous-ensembles de fracture sont sélectionnés, et les cellules adjacentes au premier sous-ensemble de fracture sont raffinées, comme il est décrit ci-dessus dans le procédé 200. Comme il a été précédemment décrit, le raffinement peut être réalisé en utilisant, par ex., la technique AGAR ou une technique plus simple telle que le Stadia Meshing. Au niveau du bloc 508, le système de simulation hybride 100 intègre le premier sous-ensemble de fracture à l'intérieur de la grille géocellulaire en 3048300 17 3D. Cependant, au niveau du bloc 508, le système n'intègre aucune propriété de continuum de fracture dans le continuum matriciel (comme pour le procédé 200) ; au lieu de cela, le système raffine simplement les cellules selon les besoins. Mais également, la grille géocellulaire pourrait être une grille préexistante qui a été raffinée ou une description de réservoir (par laquelle la grille a été générée sous forme d'une grille à haute résolution). Néanmoins, la sortie obtenue au niveau du bloc 508 est une grille géocellulaire ayant la résolution spatiale souhaitée proche des fractures hydrauliques mais avec des propriétés représentant seulement la matrice. Au niveau du bloc 510, les propriétés de la fracture du premier sous-ensemble de fracture et du second sous-ensemble de fracture sont ensuite mises à l'échelle afin de déterminer les propriétés du continuum de fracture. Toute technique d'approximation souhaitée peut être utilisée ici, comme il a été précédemment décrit. Ici, le système utilise des propriétés du premier sous-ensemble de fracture, du second sous-ensemble de fracture et de la grille raffinée pour faire une approximation des propriétés du continuum de fracture. La sortie obtenue et une grille géocellulaire en 3D qui peut être utilisée pour une diversité de simulations de réservoir. Dans certains procédés alternatifs de la présente divulgation, le raffinement du bloc 506 peut comprendre le remplacement des propriétés matricielles des cellules du premier sous-ensemble de fracture par des propriétés du continuum de fracture des cellules du premier sous-ensemble de fracture. Par conséquent, lorsque la mise à l'échelle du bloc 510 est réalisée, les propriétés du continuum de fracture du premier sous-ensemble de fracture sont déjà intégrées alors que les propriétés du continuum matriciel ont été éliminées des cellules du premier sous-ensemble de fracture. Dans certains autres procédés illustratifs de la présente divulgation, les grilles géocellulaires hybrides en 3D peuvent être validées. Ici, le modèle hybride peut être comparé à un modèle correspondant à toutes les fractures naturelles modélisées directement à des fins d'assurance de la qualité. Ceci pourrait ne pas toujours être possible si les fractures naturelles sont trop nombreuses pour le système informatique dans lequel la 3048300 18 modélisation est réalisée, mais sur des systèmes informatiques suffisamment puissants cette comparaison peut être utilisée pour valider l'approximation réalisée dans le modèle hybride. La validation sera réalisée, par ex., en comparant la production du puits (comprenant, par ex., les débits et les pressions au fond du puits, les pressions à la tête du puits, etc.) aussi bien que la réponse dynamique du réservoir. Cette dernière peut être caractérisée par les propriétés récurrentes telles que, par ex., les pressions de phase, les saturations de phase, etc. Si la validation démontre que le modèle approximatif (hybride) n'est pas suffisamment précis pour l'application courante, l'utilisateur peut utiliser les résultats de la simulation pour manuellement éliminer par filtrage du modèle hybride les fractures naturelles qui n'ont aucun impact ou très peu d'impact. Par ex., un grand nombre de fractures naturelles pourrait n'avoir aucun impact sur les profils de diminution de la pression. Ainsi, leur élimination par filtrage pour le modèle de simulation hybride n'aura aucun impact important sur les résultats de simulation du modèle hybride. À la lumière de ce qui précède, la figure 6 est un procédé 600 permettant de valider une grille géocellulaire hybride en 3D, selon certains procédés illustratifs de la présente divulgation. Après génération des procédés de la grille géocellulaire en 3D 200 ou 500, une seconde grille géocellulaire en 3D comportant un second réseau de fracture dans celui-ci est générée au niveau du bloc 602. Au niveau du bloc 604, les cellules dans la grille proches ou dans le second réseau de fracture sont raffinées. Au niveau du bloc 606, le second réseau de fracture est intégré à l'intérieur de la seconde grille géo cellulaire en 3D. Ensuite, au niveau du bloc 608, la grille géocellulaire hybride en 3D est validée par comparaison avec une seconde grille géocellulaire en 3D. Ensuite, tel qu'il a été présenté ci-dessus, certaines fractures peuvent être éliminées du modèle hybride lorsque ces fractures n'ont aucun impact important sur les résultats de la simulation. Dans encore d'autres procédés illustratifs, les temps d'exécution de la simulation peuvent être davantage écourtés à travers 3048300 19 l'utilisation d'une simulation en phase unique (par ex., gaz seulement), l'utilisation de cellules de grilles plus brutes dans le modèle de simulation, etc. Au cours du test des procédés illustratifs décrits ici, la technique hybride a été comparée aux procédés discrets et doubles utilisant le même cas basé sur un gaz en phase unique. Le procédé discret comprenait une formulation de continuum unique dans laquelle une résolution fine de maillage est fournie à proximité de toutes les fractures. Le procédé double impliquait une formulation dans laquelle l'intégralité de l'ensemble des fractures naturelles a été mise à jour par rapport au continuum de fracture et les fractures hydrauliques ont été modélisées dans le cadre du continuum matriciel. Le cas de comparaison double a été simulé sur la grille obtenue après résolution des fractures hydrauliques contenues dans celle-ci. Le cas discret a été simulé sur une grille à porosité unique qui a résolu les fractures hydrauliques et toutes les fractures naturelles se trouvant dans celles-ci. Deux types de procédés hybrides ont été utilisés : pour le procédé hybride de Type 1, les RFN majeurs (c.-à-d., le premier sous-ensemble de fracture), avec les fractures hydrauliques, ont été d'abord résolus dans la grille du continuum matriciel, et ensuite les RFN mineurs (c.-à-d., le second sous-ensemble de fracture) ont été mis à l'échelle dans la grille de continuum de fracture. Ce dernier a été réalisé à travers le même algorithme de mise à jour que pour les cas doubles. Pour le procédé hybride de Type 2, la mise à l'échelle a été réalisée suivie d'un raffinement : les fractures naturelles mineures ont été mises à l'échelle jusqu'à la grille de continuum de fracture, et ensuite la grille a été raffinée pour intégrer les fractures hydrauliques et les fractures naturelles majeures. Les fractures naturelles majeures sont celles qui sont directement intégrées dans la grille de continuum matriciel, alors que les fractures mineures sont considérées comme une approximation à travers la mise à l'échelle et sont représentées dans la grille de continuum de fracture. Dans ce cas, les fractures naturelles majeures été les cinq grandes fractures qui sont plus longues que les fractures plus courtes qui sont regroupées plus à proximité ensemble. 3048300 Une comparaison des résultats de simulation par rapport à la production gazière cumulative est illustrée dans la 7. En tenant compte du cas discret comme étant celui le plus précis, les procédés double et hybride de Type 2 ont donné de bons résultats, et même le procédé hybride de Type 2 pourrait être approprié pour des études exploratoires. Le cas double était plus rapide à simuler (14 secondes) que le cas discret (67 min en série ou 18 min sur 8 processus). Cependant, les cas hybrides étaient plus proches du cas double en termes de performance (251 sec pour l'hybride de Type 1 et 560 sec pour l'hybride de Type 2). La variation au niveau des résultats entre les deux types hybrides peut être attribuée à sa nature partiellement double, et elle peut donc rencontrer le même problème que la technique double dans laquelle le choix d'un procédé de mise à l'échelle peut engendrer une différence importante dans les résultats, qui est difficile à prédire à priori. La distribution de la pression à partir des résultats de la simulation (Fig. 8) pour les quatre cas ci-dessous, comme on s'y attendait, le procédé discret a la meilleure résolution et le procédé double la pire. Mais les champs de pression pour les deux hybrides sont plus proches en résolution au discret qu'au double. Plus spécifiquement, la figure 8 illustre une comparaison des champs de pression dans la matrice à 2000 jours, illustrée pour le réservoir en entier (à gauche) et une vue détaillée de la zone (à droite) mise en évidence par la boîte de Type 2 hybride (la ligne du bas de la première colonne). Les modes de réalisation et les procédés illustrés décrits ici procurent un certain nombre d'avantages. Les procédés divulgués améliorent l'efficacité d'un simulateur de réservoir lors de la planification du développement des actifs gaziers et pétroliers en améliorant la performance du système tout en générant la grille de simulation, aussi bien que le temps d'exécution des simulations pour ces actifs lorsque des fractures naturelles sont connues pour ou suspectés d'affecter la productivité du puits. La performance du système est améliorée principalement en représentant un sous-ensemble du RFN directement sur la grille, alors que l'approximation du restant des RFN est réalisée par des formulations à double continuum. En 3048300 21 outre, la présente divulgation décrit également un mécanisme de filtrage permettant de sélectionner interactivement lequel des sous-ensembles de RFN doit être maillés directement. Les divers procédés décrits ici peuvent être totalement ou partiellement automatisés. Du point de vue de l'utilisateur, il peut sembler que le procédé soit constitué de moins de trois étapes. Les modèles géocellulaires hybrides en 3D générés peuvent être utilisés dans une diversité d'applications, telles que dans la simulation de réservoir. Une application principale et souvent unique de la simulation de réservoir consiste en l'estimation des réserves récupérables des actifs gaziers et pétroliers. La capacité de créer et de simuler un modèle de simulation hybride rend possible la prédiction plus précise et plus opportune du recouvrement final (EUR) d'un réservoir naturellement fracturé, conventionnel et non-conventionnel, d'une taille typique ou d'une taille très grande. Ensuite, la simulation du réservoir peut être utilisée pour planifier l'exploitation d'un actif gazier et pétrolier naturellement fractionné. Divers scénarios impliquant des procédés multiphasiques complexes dans des formations souterraines, tels que la récupération de pétrole par injection d'eau ou injection de gaz, une fois simulée, fournissent les bases pour l'optimisation des réservoirs naturellement fractionnés. Une autre utilisation des modèles hybrides de simulation de réservoir consiste en la réalisation d'une analyse d'incertitude pour de grands réservoirs naturellement fractionnés. Dans ces flux de travail itératifs, plusieurs variantes (quelquefois des centaines ou des milliers) d'un modèle de réservoir de base doivent être simulées et donc des simulations à exécution rapide sont nécessaires. L'analyse de l'incertitude quantifie la fourchette des résultats (par ex., des facteurs de récupération, la valeur nette présente, la production cumulative à un moment donné dans le futur) qu'on attend de l'actif lorsque les distributions de probabilité fournies des variables dans le modèle (par ex., la distribution des RFN, les perméabilités matricielles, les porosités, les faciès géologiques) qui ne sont pas connues précisément, mais affecteraient de façon importante le résultat. L'optimisation d'un modèle de complétion d'un puits de forage dans le schiste, tel que le Eagle Ford, dans 3048300 22 lequel non seulement la perméabilité réelle du réservoir mais également le schéma de propagation de la fracture hydraulique est significativement influencé par les fractures naturelles, bénéficie de l'utilisation de modèles hybrides du type décrit dans la présente demande. En outre, une étude « rétrospective » peut être réalisée en utilisant un modèle de simulation hybride pour déterminer les causes d'un échec de complétion (y compris la fracturation hydraulique) dans un puits à problèmes et pour améliorer de telles opérations de fond de puits dans des puits futurs. Les modes de réalisation et les procédés décrits ici concernent également l'un ou plusieurs des paragraphes suivants : 1. Un procédé implémenté par ordinateur permettant de générer une grille géocellulaire hybride en 3 dimensions (« 3D »), le procédé comprenant l'utilisation des données provenant d'un réservoir pour générer un réseau de fracture ; la sélection d'un premier sous-ensemble de fracture dans le réseau de fracture ; la sélection d'un second sous-ensemble de fracture, différent du premier sous-ensemble de fracture, dans le réseau de fracture ; le raffinement des cellules d'une grille géocellulaire en 3D à proximité du premier sous-ensemble de fracture ; la mise à l'échelle du second sous-ensemble de fracture et l'intégration du premier sous-ensemble de fracture et du second sous-ensemble de fracture mis à l'échelle à l'intérieur de la grille géocellulaire en 3D, générant ainsi une grille géocellulaire hybride en 3D. 2. Un procédé implémenté par ordinateur tel que défini dans le paragraphe 1, dans lequel le raffinement des cellules comprend la division d'un bord des cellules. 3. Un procédé implémenté par ordinateur tel que défini dans les paragraphes 1 ou 2, dans lequel la mise à l'échelle du second sous-ensemble de fracture comprend la réalisation d'une approximation multicontinuum utilisant le second sous-ensemble de fracture. 4. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 1 à 3, dans lequel le premier sous-ensemble 3048300 23 de fracture est choisi en utilisant un facteur de connectivité représentant une distance connectée de chaque fracture à une perforation de puits. 5. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 1 à 4, dans lequel le premier et le second sous-ensembles de fracture sont choisis par stéréonet, visualisation ou attribut de fracture. 6. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 1 à 5, comprenant également la génération d'une seconde grille géocellulaire en 3D comportant un second réseau de fracture à l'intérieur de celle-ci ; le raffinement des cellules dans la seconde grille géocellulaire en 3D à proximité du second réseau de fracture ; l'intégration du second réseau de fracture à l'intérieur de la seconde grille géocellulaire en 3D ; et la validation de la grille géocellulaire hybride en 3D en comparant la grille géocellulaire hybride en 3D à la seconde grille géocellulaire en 3D. 7. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 1 à 6, comprenant également l'enlèvement des fractures de la grille géocellulaire hybride en 3D basé sur la validation. 8. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 1 à 7, comprenant également la réalisation d'une simulation de réservoir utilisant la grille géocellulaire en 3D. 9. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 1 à 8, dans lequel la simulation du réservoir est utilisée pour planifier, réaliser ou analyser une opération de puits de forage. 10. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 1 à 9, dans lequel la simulation du réservoir est utilisée pour prédire la productivité d'un puits. 11. Un procédé implémenté par ordinateur permettant de générer une grille géocellulaire hybride en trois dimensions (« 3D »), le procédé comprenant l'utilisation des données provenant d'un réservoir pour générer un réseau de fracture comportant une pluralité de cellules dans celui- 3048300 24 ci ; la sélection d'un premier sous-ensemble de fracture dans le réseau de fracture ; la sélection d'un second sous-ensemble de fracture, différent du premier sous-ensemble de fracture, dans le réseau de fracture ; la génération d'une grille géocellulaire en 3D ayant une haute résolution à proximité du premier sous-ensemble de fracture ; l'intégration du premier sous-ensemble de fracture à l'intérieur de la grille géocellulaire en 3D ; et la mise à l'échelle du premier sous-ensemble de fracture et du second sous-ensemble de fracture, générant ainsi une grille géocellulaire hybride en 3D. 12. Un procédé implémenté par ordinateur tel que défini dans le paragraphe 11, dans lequel la génération d'une grille géocellulaire en 3D comprend la division d'un bord des cellules à proximité du premier sous-ensemble de fractures ; ou la génération directe d'une grille géocellulaire en 3D à haute résolution utilisant une description de réservoir. 13. Un procédé implémenté par ordinateur tel que défini dans les paragraphes 11 ou 12, dans lequel la génération de la grille géocellulaire en 3D comprend également le remplacement des propriétés matricielles des cellules par des propriétés de fracture du premier sous-ensemble de fracture. 14. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 11 à 13, dans lequel la mise à l'échelle comprend la réalisation d'une approximation multi-continuum utilisant le premier sous-ensemble de fracture, le second sous-ensemble de fracture et la grille à haute résolution à proximité du premier sous-ensemble de fracture. 15. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 11 à 14, dans lequel le premier sous-ensemble de fracture est choisi en utilisant un facteur de connectivité représentant une distance connectée de chaque fracture à une perforation de puits. 16. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 11 à 15, dans lequel le premier et le second sous-ensembles de fracture sont choisis par stéréonet, visualisation ou attribut de fracture. 3048300 17. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 11 à 16, comprenant également la génération d'une seconde grille géocellulaire en 3D comportant un second réseau de fracture à l'intérieur de celle-ci ; le raffinement des cellules dans le second réseau de fracture ; l'intégration du second réseau de fracture à l'intérieur de la seconde grille géocellulaire en 3D ; et la validation de la grille géocellulaire hybride en 3D en comparant la grille géocellulaire hybride en 3D à la seconde grille géocellulaire en 3D. 18. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 11 à 17, comprenant également l'enlèvement des fractures de la grille géocellulaire hybride en 3D basé sur la validation. 19. Un procédé implémenté par ordinateur tel que défini dans la quelconque des paragraphes 11 à 18, dans lequel la grille géocellulaire en 3D utilisée dans l'intégration du premier sous-ensemble de fracture a été fournie sous forme d'une grille géocellulaire en 3D préexistante ; ou elle a été générée en utilisant une description de réservoir. 20. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 11 à 19, comprenant également la réalisation d'une simulation de réservoir utilisant la grille géocellulaire en 3D. 21. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 11 à 20, dans lequel la simulation du réservoir est utilisée pour planifier, réaliser ou analyser une opération de puits de forage. 22. Un procédé implémenté par ordinateur tel que défini dans l'un quelconque des paragraphes 11 à 22, dans lequel la simulation du réservoir est utilisée pour prédire la productivité d'un puits. En outre, l'un quelconque des procédés illustratifs décrits ici peuvent être implémentés par un système comprenant un circuit de traitement ou un support non-transitoire lisible par ordinateur contenant des instructions qui, lorsqu'elles sont exécutées par au moins un processeur, amène le processeur à réaliser l'un quelconque des procédés décrits ici. 3048300 26 Même si divers modes de réalisation et procédés ont été illustrés et décrits, la divulgation n'est pas limitée à de tels modes de réalisation ou procédés et on comprendra quelle englobe des modifications et des variations qui seront évidentes à un spécialiste du domaine. Par conséquent, il doit être compris que la divulgation n'est pas destinée à être limitée aux formes particulières divulguées. Au lieu de cela, l'invention doit couvrir toutes les modifications, les équivalents et les alternatifs qui sont à l'intérieur de l'esprit et de la portée de la divulgation, telle qu'elle est définie dans les revendications ci-jointes. REVENDICATIONS1. Procédé implémenté par ordinateur pour générer une grille géocellulaire hybride en trois dimensions (« 3D »), le procédé comprenant : l'utilisation des données de réservoir pour générer un réseau de fracture ; la sélection d'un premier sous-ensemble de fracture dans le réseau de fracture ; la sélection d'un deuxième sous-ensemble de fracture, différent du premier sous-ensemble de fracture, dans le réseau de fracture ; le raffinement des cellules d'une grille géocellulaire en 3D à proximité du premier sous-ensemble de fracture ; la mise à l'échelle du second sous-ensemble de fracture ; et l'intégration du premier sous-ensemble de fracture et le deuxième sous-ensemble de fracture mis à l'échelle dans la grille géocellulaire en 3D, générant ainsi une grille géocellulaire hybride en 3D. 2. Procédé implémenté par ordinateur tel que défini dans la revendication 1, dans lequel le raffinement des cellules comprend la division d'un bord des cellules. 3. Procédé implémenté par ordinateur tel que défini dans la revendication 1, dans lequel la mise à l'échelle du second sous-ensemble de fracture comprend la réalisation d'une approximation multi-continuum utilisant le second sous-ensemble de fracture. 4. Procédé implémenté par ordinateur tel que défini dans la revendication 1, dans lequel le premier sous-ensemble de fracture est choisi en utilisant un facteur de connectivité représentant une distance connectée de chaque fracture à une perforation de puits. 5. Procédé implémenté par ordinateur tel que défini dans la revendication 1, dans lequel le premier et le second sous-ensembles 3048300 28 de fracture sont choisis par stéréonet, visualisation ou attribut de fracture. 6. Procédé implémenté par ordinateur tel que défini dans la revendication 1, comprenant également : la génération d'une grille géocellulaire en 3D ayant un second réseau de fracture à l'intérieur de celle-ci ; le raffinement des cellules de la grille géocellulaire en 3D à proximité du second sous-ensemble de fracture ; l'intégration du second réseau de fracture à l'intérieur de la seconde grille géocellulaire en 3D ; et la validation de la grille géocellulaire hybride en 3D en comparant la grille géocellulaire hybride en 3D à la seconde grille géocellulaire en 3D. 7. Procédé implémenté par ordinateur tel que défini dans la revendication 6, comprenant également l'enlèvement des fractures de la grille géocellulaire hybride en 3D basé sur la validation. 8. Procédé implémenté par ordinateur tel que défini dans la revendication 1, comprenant également la réalisation d'une simulation de réservoir utilisant la grille géocellulaire en 3D. 9. Procédé implémenté par ordinateur tel que défini dans la revendication 8, dans lequel la simulation du réservoir est utilisée pour planifier, réaliser ou analyser une opération de puits de forage. 10. Procédé implémenté par ordinateur tel que défini dans la revendication 8, dans lequel la simulation du réservoir est utilisée pour prédire la productivité d'un puits. 11. Support lisible par ordinateur comprenant des instructions qui, lorsqu'elles sont exécutées par au moins un processeur, amènent le processeur à réaliser le procédé tel que décrit dans la revendication 1. 12. Système comprenant un circuit pour réaliser un procédé tel que défini dans la revendication 1. 3048300 29 13. Procédé implémenté par ordinateur pour générer une grille géocellulaire hybride en trois dimensions (« 3D »), le procédé comprenant : l'utilisation des données de réservoir pour générer un réseau de fracture ayant une pluralité de cellules dans celui-ci ; la sélection d'un premier sous-ensemble de fracture dans le réseau de fracture ; la sélection d'un deuxième sous-ensemble de fracture, différent du premier sous-ensemble de fracture, dans le réseau de fracture ; le génération d'une grille géocellulaire en 3D à haute résolution à proximité du premier sous-ensemble de fracture ; l'intégration du premier sous-ensemble de fracture à l'intérieur de la grille géocellulaire en 3D ; et la mise à l'échelle du premier sous ensemble de fracture et le second sous-ensemble de fracture, générant ainsi une grille géocellulaire hybride en 3D. 14. Procédé implémenté par ordinateur tel que défini dans la revendication 13, dans lequel la génération de la grille géocellulaire en 3D comprend : la division d'un rebord de cellules à proximité du premier sous-ensemble de fracture ; ou la génération directe d'une grille géocellulaire en 3D à haute résolution utilisant une description de réservoir. 15. Procédé implémenté par ordinateur tel que défini dans la revendication 14, dans lequel la génération de la grille géocellulaire en 3D comprend également le remplacement des propriétés matricielles des cellules par des propriétés de fracture du premier sous-ensemble de fracture.