Brevet FR3058447A1 — publié le 20180511

2016-WM- 00753-1.1 - INCLUSION DE DIFFUSION SÉLECTI E POUR UNE SIMULATION D RÉSERVOIR POUR LA RÉCUPÉRA TION DES HYDROCARBURES CON EX FIL. 100011 La divulgation porte en général sur le traitement des données pour la récupération des hydrocarbures, et plus particulièrement la simulation des réservoirs pour la récupération des hydrocarbures. 100021 Les simulations tic réservoirs peuvent prévoir l'écoulement et le comportement de phase des fluides dans un réservoir pour prévoir les quantités de production et d'injection qui incluent la récupération progressive d'hydrocarbures, les besoins en solvants miscibles, l'efficacité de l'utilisation des solvants. etc. 1.es simulations peuvent aider à évaluer les procédés de récupération assistée de pétrole t LOR). y compris l'injection de gaz et l'inondation chimique. L'injection de gd/ en général peut etre une bonne option pour l'EOR, en particulier pour les réservoirs fracturés, où il existe une esande zone de 15 contact entre le gaz injecté et le fluide en place. 100031 Au cours des procédés de récupération par injection de gaz, le gaz ou le fluide pressurisé peut être transmis à partir d'un puits de forage dans le réservoir à pression élevée, et le gaz pressurisé deplace le pétrole à l'intérieur de la roche de réservoir. Le pétrole déplacé et le gaz injecté peuvent ensaite être produits dans des puits de forage de 20 production. La création d'un modèle compositionnel fiable pour .repre,enter ces procédés nécessite des caractéristiques spéciales dans fa simulation, telles que des modèles pour la diffusion moléculaire entre le gaz injecté à l'intérieur des fractures et les fluides stockés dans la roche de réservoir. Cependant. le calcul précis de la diffusion moléculaire est généralement un procédé coûteux et compliqué au cours de la simulation. 25 MT VE DESCRIPTION DES FIGURES 100041 Les modes de réalisation rie la divulgation peuvent être mieux compris en faisant référence aux illustrations ci-jointes. 100051 La figure 1 est un diagramme schématique d'un exemple de système de puits et d'un sous-système informatique, selon certains modes de réalisation. 3058447 2016-1PM-100753-U l-hR 2 100061 Li figure L. décrit une représentation visuelle des mailles d'un réservoir dans une simulation de diffusion souterraine, selon certains modes de réalisation. 100071 La ligure 3 décrit un diagramme de flux pour contourner les calculs de fliffUSion dans une simulation de diffusion moketilaire pour les mailles d'un réservoir qui sont 5 dominées par le flux de convection, selon certains modes de réalisation. 100081 La figure 4illustre un exemple d'llfl dispositif informatique, selon certains modes de réalisation. DESC.i'RIPTION 1(101)91 La description qui suit comprend des exemples de systèmes. de procédés, de I0 techniques et de flux de programme qui concrétisent des aspects de la divulgation. Cependant, il est entendu que cette divulgation peut etre pratiquée en absence de ces détails spécifiques. Dans d'autres cas, des structures et des techniques bien connues n'ont pas été illustrées en détail afin de ne pas obscurcir la description. 100101 Divers modes de réalisation se rapportent à des simulations de récupération de 15 pétrole à partir de réservoirs de pétrole. Certains modes de réalisation comprennent des opérations pour réduire le temps de calcul de la diffusion moléculaire lors d'une simulation d'injection de gaz dans un réservoir en milieu fracturé, afin de réduire le temps d'exécution totale de la simulation. I.a diffusion moléculaire peut jouer un rôle important dans la récupération des hydrocarbures à partir dés réservoirs fractionnés. La diffusion 20 moléculaire peut comprendre tout mouvement thermique de particules ( y compris le liquide et le gai), à des températures au-dessus du zéro absolu. La diffusion moléculaire peut être un facteur important dans la récupération d'hydrocarbures lors des injections de gaz lorsque le drainage gravitationnel est inefficace. Un écoulement efficace entre des fractures et des blocs matriciels peut etre engendré en raison des gradients 25 compositionnels. Mais également, la diffusion de composants légers dans la phase de pétrole peut entraîner des changements favorables dans les propriétés des fluides, tels que la réduction de la viscosité et le gonflement du pétrole dans la matrice. Le fait d'ignorer la diffusion dans la simulation petit entraîner une sous-estimation de hi récupération d'hydrocarbures. Sans la diffusion, le gaz injecté s'écoule principalement à travers les 30 fractures, ce qui entraîne une percée précoce. Dans certains modes de réalisation, un volume simulé du réservoir peut être subdivisé en un certain nombre de mailles qui 3058447 2016-1PM-100753-U - peuvent interagir entre elles au travers de leurs interfaces de maille au cours de la réalisation des calculs de diffusion. 100111 Tes modèles existants pour la diffusion sont soit des modèles t coefficient constant simples qui ne peuvent appréhender correctement les éléments physiques, soit 5 des modèles sophistiqués coûteux qui peuvent avoir des coûts de calcul élevés pour des problèmes réels. Certains modes de réalisation fournissent une procédure permettant de réduire le temps pour réaliser des calculs de diffusion sans perte de précision. Des opérations pour accélérer les calculs de diffusion peuvent être particulièrement utiles pour des cas de sites utilisant un grand nombre de composants dans les simulations. Le contournement ou l'omission des calculs de diffusion inutiles est particulièrement importante parce qu'il peut y avoir des zones dans lesquelles il n'y a pas de grand contraste de composition entre des mailles voisines. 100121 Les calculs de diflUsion peuvent impliquer de nombreuses opérations matricielles, comprenant les calculs des coefficients de diffusion, des flux et des dérivés 15 des termes de flux. La réalisation de ces opérations pour chaque maille d'un réservoir à simuler peut être très chronophage. La diffusion peut être négligeable pour certaines mailles qui ne sont pas en contact avec le fluide injecté ou pour celles qui sont complètement balayées par le solvant et pour lesquelles il ne reste rien pour l'échange de masse. En outre, pour les réservoirs lietérogenes avec (titi-Ci-entes propriétés de matrice et 20 de fracture. il peut exister des zones avec des 1-écimes d'écoulement différents allant des régimes dominés par la diffusion (principalement proches des fractures) jusqu'aux régimes dominés par ia convection (éloignés des fractures). La réalisation (les calculs de diffusion pour les mailles a diffusion négligeable ou pour les mailles dominées seulement par la convection entraîne une plus grande complexité et un temps d'exécution plus long 25 au cours de la simulation. 100131 Dans certains modes de réalisation, une simulation de réservoir suit dynamiquement les zones au niveau desquelles la diffusion est un mécanisme de transfert massique efficace et évite les calculs de diffusion inutiles dans les zones de di frasion négligeable en comparaison à la convection. 1...a zone la plus active au niveau de laquelle 30 le transfert massique par diffusion maximal se produit peut se situer autour du front de fluide. où les forces d'entraînement sont a la puissance maximale en raison de la grande différence dans les propriétés du fluide. Le front de fluide comporte une zone de 3058447 2016-IPM-100753-U - 4 transition, dans laquelle la composition passe de la composition du fluide injecté à celle du fluide en place. L'étendue de la 7011C de transition dépend du fluide et des propriétés de milieu poreux. Les grandes diffusivités peuvent entraîner une propagation plus éloignée de molécules et créer ainsi une zone de transition plus large. Inversement, si la vitesse 5 globale est très élevée, le temps disponible pour la propagation de molécules peut être faible et, par conséquent_ l'effet de cii fusion devient insignifiant. Par conséquent, l'extension de la zone de transition peut être corrélée au rapport du flux cle convection sur le flux de diffusion, ce qui représente la définition d'un nombre de Péclet. 100141 Dans certains modes de réalisation, un résultat de la simulation de réservoir peut être utilisé pendant une opération réelle de récupération d'hydrocarbures. Par exemple, les emplacements de l'un ou des deux d'un puits d'injection de gaz ou d'un puits de production peuvent être détermines en fonction d'un résultat de la simulation de réservoir. De plus, l'un ou l'autre ou les deux d'un débit et d'une composition d'un gaz d'injection à injecter dans le puits d'injection de gaz pendant l'opération (le récupération. 15 d'hydrocarbures peuvent être basés sur le résultat de la simulation de réservoir. Exemple de système 100151 figure 1 est un schéma d'un exemple de système de puits et d'un sous-système informatique, selon certains modes de réalisation. Un système de puits 100 comprend un puits d'injection 102 et un puits de production 103 dans une zone souterraine 104 en dessous de la surface de sol 106. Le puits d'injection 102 et le puits de production 103 illustrés dans la figure I sont représentés comme des puits verticaux. Cependant, certains modes de réalisation peuvent être incorporés dans des systèmes de puits qui comprennent toute combinaison d'orientation horizontale. verticale, inclinée, courbée ou d'autres orientations de puits de foraue. outre, bien qu'il soit représente avec un seul puits 25 d'injection et un puits de production, le système de puits 100 peut comprendre un ou plusieurs puits de traitement supplémentaires, des puits d'observation, des puits de production, etc. 100161 Le sous-système informatique 110 peut comprendre un ou plusieurs dispositifs ou systèmes informatiques situés au niveau du puits d'injection 102 et du puits de 30 production 103 ou à d'autres emplacements. I_e sous-système informatique 110, ou un quelconque de ses composants, peut être situé à l'écart des autres composants illustrés dans la figure 1. Par exemple, le sous-système informatique 110 peut être situé dans un 3058447 2016-IPM-100753-U - centre de traitement de données, dans une installation informatique, ou dans un autre emplacement approprié. Le système de puits 100 peut comprendre des caractéristiques additionnelles ou différentes, et les caractéristiques du système de puits 100 peuvent être disposées telles comme le montre la figure 1 ou dans une autre configuration. 5 10017] La zone souterraine 10-1- peut comprendre un réservoir qui contient des ressources d'hydrocarbures telles que du pétrole, du gaz naturel, ou autres. Par exemple, la zone souterraine 104 peut comprendre toute ou une partie d'une formation rocheuse (par ex., le schiste, le charbon, le grès, le granite, ou autres) qui contient du gaz naturel. La zone souterraine 104 peut comprendre des formations de roches naturellement fracturées ou de roches naturelles qui ne sont pas fracturées à un quelconque degré significatif. a zone souterraine 104 peut comprendre des formations étanches au gaz qui comprennent des roches à faible perméabilité (par ex., é schiste. le charbon, ou autres). 100181 Le système de puits 100 illustré dans la ligure 1 comprend un système d'injection 108.11...e système d'injection 108 peut être utilisé pour réaliser un traitement d'injection par 15 lequel un gaz est injecté dans la zone souterraine 104 dans le puits d'injection 102. Par exemple, le système d'injection 108 peut comprendre une pompe d'injection pour injecter un gaz de traitement dans la zone souterraine 104 dans le puits d'injection 102. Par exemple, un procédé de déplacement de gaz d'injection peut etre appliqué na niveau d'un emplacement d'injection unique ou à de multiples emplacements d'injection dans une 20 zone souterraine. et le gaz peut être injecté au cours d'une période de temps unique ou au cours de multiples périodes de temps différentes. Dans certains cas. un schéma de récupération d'injection de gaz peut utiliser (le multiples emplacements différents d'injection de gaz dans un seul puits de tbrage, de multiples emplacements d'injection de gaz dans de multiples puits de forage différents, ou une quelconque combinaison 25 appropriée. En outre, le schéma de récupération d'injection de gaz peut injecter du gaz à travers un quelconque type de puits de forage approprié tel que, par ex., des puits de forage verticaux. des puits de forae inclinés, des puits de fora - horizontaux, des puits de forage courbés, ou des combinaisons de ceux-ci ou d'autres. 10019] En réponse à l'injection de gaz de traitement dans le puits d'injection 102 et en 30 raison d'un gradient potentiel, des fluides de réservoir peuvent s'écouler dans le puits de production 103 au travers d'un conduit de production 162. Les fluides de réservoir peuvent être récupérés du puits de production 103. Même s'il n'est pas illustré, le système 3058447 2016-1PM-100753-U - 6 de puits 100 peut également comprendre des systèmes de commande de production et des installations en surface permettant de récupérer et de traiter les fluides de réservoir pro\.cnant du puits de production 103. Le système de puits 100 peut également comprendre des installations de séparation en surface, des pipelines, des unités de 5 stockage, etc., pour un nouveau traitement, stockage et transport de fluides de réservoir récupérés du puits de production 103. là) outre, le système de puits 100 peut produire des fluides de réservoir et injecter du gaz à partir de multiples emplacements dans la zone souterraine. En outre, la production peut être réalisée à un quelconque point avant, pendant et après l'injection de gaz de traitement. La production peut également être réalisée à partir de multiples zones à l'intérieur du même puits de forage. De plus, bien que le système de puits 100 illustre un puits unique de production, la production peut également être réalisée à partir d'une quelconque combinaison de puits verticaux, déviés et horizontaux. 100201 Le gaz peut 1etre fourni à partir d'un camion muni d'un compresseur, ou à partir 15 d'un pipeline de gaz et d'installations de compresseur en surface. Le gaz de traitement peut être transmis à travers le puits d'injection 102 à partir de la surface du sol 106 par un conduit crin_jcction 112 installé dans le puits d'injection 102. Le conduit de production 162 et le conduit d'injection 1 12 peuvent inclure un cuvelage cimenté à la paroi du puits d'injection 102. Dans certaines réalisations, on peut laisser Ouverte la totalité ou une partie 20 du puits d'injection 102, sans cuvelage. 1.e conduit de production 162 et le conduit d'injection 112 peuvent inclure une colonne de travail, un tubage enroulé, un tuyau articulé, ou d'autres types de conduit. 100211 Le système d'injection 108 peut également comprendre des capteurs en surface et en tond de puits pour mesurer la pression, le débit, la température ou d'autres paramètres 25 de traitement ou de production. Par exemple, le système d'injection 105 peut comprendre des manomètres ou un autre équipement qui mesure la pression dans le puits d'injection. 102 au niveau de ou proche de la surface du sol 106 ou à d'autres emplacements. Le système d'injection 105 peut également comprendre des commandes de pompe et ou d'autres types de commandes pour démarrer, arrêter, augmenter, diminuer ou commander. 30 autrement le pompage, ainsi que des commandes pour sélectionner ou commander autrement des gaz pompés au cours du traitement d'injection. Le système d'injection 108 peut comprendre un sous-système de commander de traitement par injection pour 3058447 2016-1PM-100753-U 1-hR 7 communiquer avec l'équipement afin de surveiller et commander le traitement par injection. 100221 La zone souterraine 104 peut inclure les fractures naturelles 151 à 156. Ln variante ou en plus, un traitement d'injection de fluide peut également créer des fractures 5 dans la zone souterraine 104, ou stimuler davantage les fractures naturelles 1: I ti 15(i. Généralement, les fractures peuvent comprendre des fractures de tout type quelconque, nombre, longueur, forme, géométrie ou ouverture. Les fractures peuvent se prolonger dans une quelconque direction ou orientation, et elles peuvent etre formées à de multiples étapes ou intervalles, à des moments différents ou simultanément. Les fractures peuvent se prolonger à travers la roche naturellement fracturée, ou dans des zones de roche non-fracturée, ou les deux. 1.es fractures peuvent également etre reliées au système de production et peuvent constituer le conduit principal pour ta production à partir du réservoir vers les puits de forage de production. 100231 Dans certaines réalisations, le sous-système informatique 110 peut exécuter des 15 instructions pour simuler le réservoir de pétrole dans le système de puits 100 au cours des opérations d'injection de gaz. Le sous-système informatique 110 peut réaliser des simulations avant, pendant ou après le traitement d'injection. Dans certaines réalisations, le sous-système de commande de traitement d'injection commande le traitement d'injection en se basant sur les simulations réalisées par le sous-système informatique 20 110. Par exemple, un schéma de pompage ou d'autres aspects d'un plan d'injection de gaz peuvent être générés à l'avance en se basant sur des simulations réalisées par le sous-système informatique 110. Coinme autre exemple, le sous-système de commande dc. traitement d'injection peut modifier, actualiser ou générer un plan d'injection de gaz basi, sur des simulations réalisées par le sous-système informatique 110 en temps réel au cours 25 du traitement d'injection. Dans certaines réalisations, le sous-système de commande de production peut commander la production de puits existants, le traitement de reconditionnement de puits existants et le forage de nouveaut puits. 100241 Dans certains cas, les simulations sont basées sur les données provenant du système de puits 100. Par exemple, des manomètres, des débitmètres. un équipement 30 microsismique, des inclinomètres ou d'autres équipements peuvent réaliser des mesures avant, pendant ou après un traitement d'injection et le sous-système informatique 110 peut réaliser la simulation de réservoir eompositionnelle en se basant sur les données 3058447 2016-1PM-100753-U - mesurées. Dans certains cas, le sous-système de commande de traitement d'injection peut sélectionner ou modifier (par ex., augmenter ou diminuer) des pressions Lie gaz, des densités de gaz. dos compositions de gaz, et d'autres paramètres de commande co se basant sur les données fournies par les simulations. Dans certains cas, les données 5 fournies par les simulations peuvent être affichées en temps réel au cours du traitement d'injection, par ex., à uan. ingénieur ou à un autre exploitant du système de puits 100. 100251 La figure 2 présente une représentation visuelle de mailles d'un réservoir dans une simulation de la diffusion souterraine, selon certains modes de réalisation. Un système de mailles 200 est rectangulaire et comprend 49 mailles rectangulaires avec sept rangées et sept colonnes. Les mailles comprennent des mailles de moindre diffusion 201227, des mailles de diffusion moyenne 228 à 240 et des mailles (Ie plus grande diffusion 241 à 249. 100261 Chacune des rangées du système de mailles 200 comprend sept mailles. De haut en bas dans Tordre croissant, les rangées du système de mailles 200 peuvent être 15 ordonnées en première rangée, deuxième rangée, troisième rangée, quatrième rangée, cinquième rangée, sixième rangée et septième rangée. La première rangée du système de mailles 200 comprend les mailles de moindre diffusion 201 à 207. La deuxième rangée du système de mailles :on comprend les mailles de moindre diffusion 208 à 209, les mailles de diffusion moyenne 228 à 230 ei les mailles de moindre diffusion 210 à 21 ! . 20 troisième rangée du système de mailles 200 comprend les mailles de moindre diffusion 2 à 213, la maille de di itusion moyenne 231, la maille de plus grande diffusion 241, les mailles de diffusion moyenne 232 à 233, et la maille de moindre diffusion 214. La quatrième rangée du système de mailles 200 comprend la maille de moindre diffusion 215, les mailles de diffusion moyenne 234 à 235, les mailles de plus grande diffusion 242 25 à 244 et ia maille de moindre diffusion 2 i6. La cinquième rangée du système de mailles 200 comprend les mailles de moindre diffusion 217 a 219, la maille de diffusion moyenne 236, les mailles de plus grande diffusion 245à 246 et la maille de moindre diffusion 220. I,a sixième rangée du système de mailles 200 comprend les mailles de moindre diffusion 221 à 224, la maille de diffusion moyenne 2217 et les mailles de plus grande diffUsion 247 30 à 248. La septième rangée du système de mailles 200 comprend les mailles de moindre diffusion 225 à 227, les mailles de diffusion moyenne 238 à 240 et la maille de plus grande diffusion 249. 3058447 ")016-1PM-100753-U 1-hR 9 100271 Le flux de diffusion de chacune des mailles peut être décrit comme étant moindre. moyen ou plus grand l'un par rapport à l'autre et sans limiter aucun flux de diffusion inférieur à, approximativement égal, ou supérieur a une valeur fixe du flux. Le flux de diffusion des mailles de moindre diffitsion 201 à 227 peut être inférieur au flux de 5 diffusion de l'une quelconque des mailles de diffusion moyenne 228 à 240 ou des mailles de plus grande diffusion 241 à 149. Le flux de diffusion de chacune des mailles de diffusion moyenne 228 à 240 est supérieur au flux de diffusion de l'une quelconque des mailles de moindre diffusion 201 à 227 et inférieur au flux de diffusion de l'une quelconque des mailles de plus grande diffusion 241 à 249. Le flux de diffusion de chacune des mailles de plus grande diffusion 241 à 249 est supérieur au flux de diffusion de l'une quelconque des mailles de moindre diffusion 201 à 227 ou des mailles de diffusion moyenne 228 à 240. 100281 Par exemple, le flux de diffusion de chacune des mailles de moindre diffusion 201 à 227 peut être dans une plage de moins de 10 kilos par mètre carré par seconde 15 (kg 'm2-s). Le flux de diffusion de chacune des mailles de diffusion moyenne 228 à 240 peut être dans une plage de IO a moins de 100 kg;m2-s. Le flux de diffusion de chacune des mailles de plus grande diffusion 241 à 249 peut être un flux égal ou supérieur à 100 100291 Dans certains modes de réalisation, une maille peut être représentée en deux 20 dimensions en estimant la troisième dimension (par exemple, lu hauteur) pour être égale pour chaque maille. Dans de tels modes de réalisation bidimensionnels, chacune c es mailles peut être rectangulaire et de dimensions égaies. Alternativement, les mailles peuvent être des quadrilatères non équilatéraux de différentes tailles. Les mailles d'autres simulations de réservoir peuvent être représentées par d'autres polygones 25 bidimensionnels, des ellipsoïdes bidimensionnels, des mailles Voronoï, etc. Dans certains modes de réalisation, une maille peut être représentée en trois dimensions par des formes polyédriques. des ellipsoïdes tridimensionnels, des formes combinant des surfaces courbes et planes, etc. Par exemple, un système de simulation de réservoir tridimensionnel peut inclure des mailles en cube, des mailles en prisme hexagonal, des 30 .mailles dodécaédriques, etc. 3058447 2016-1 PM-100753-LJ l - Exemple d'opérations 100301 Des opérations sont maintenant décrites pour réduire le temps de calcul de la diffusionmoléculaire dans la simulation compositionnelle du réservoir lors des injections de gaz dans le milieu fracturé, afin de réduire le temps d'exécution total de la simulation. 5 100311 La f !n.ire 3 illustre un schéma fonctionnel pour contourner les calculs de tirusion dans une SHImlilation de diffusion moléculaire pour les mailles d-un réservoir qui sont dominées par le flux de convection, selon certains modes Oc réalisation. Les opérations d'un diagramme de flux 300 peuvent dynamiquement suivre les zones de transition, dans lesquelles la diffusion est plus efficace ou dominante, et éviter des calculs dans des zones dominées par la convection. Diminuer le temps total de simulation en. ignorant certaines opérations petit être particulièrement important pour les cas de sites comprenant des zones où la diffusion est dominée par la convection. Les opérations du diagramme de flux 301) peuvent être effectuées par un logiciel, un matériel, un microprogramme ou une combinaison de ceux-ci. Par exemple, en référence 'à un exemple de dispositi r i1tormatiqt.te représente sur la figure 4 (deerit plus en détail ci-dessous), un processeur peut exécuter (les instructions pour effectuer les opérations du diagramme rte flux 300. Les opérations du diagramme de flux 300 peuvent être effectuées pour n'importe quel nombre d'intervalles de temps dans un système de simulation de réservoir. Les opérations du diagramme de flux 300 commencent à 301. 100321 Au bloc 301, la prochaine maille est déterminée. Dans certains modes de réalisation. une simulation de réservoir peut comprendre une pluralité de mailles. variable de suivi peut être utilisée pour déterminer quelle maille est une prochaine maille. Par exemple, la simulation de résenoir peut inclure un indice de maille. dans lequel chaque valeur d'indice de l'indice de maille correspond à une maille parmi une pluralité de mailles. La détermination d'une prochaine maille peut inclure la modification (l'une variable de suivi d'une valeur d'indice antérieure à une valeur d'indice suivante. Les opérations aux blocs 302 à 320 peuvent ensuite être appliquée; à la maille correspondante à cette valeur d'indice suivante. Par exemple, en référence à la figure 2, si une variable de suivi était initialement égale à 201. la détermination d'une prochaine maille peut augmenter la variable de suivi pour etre égale à 202. Dans certains modes de réalisation, si la variable de suivi était à l'origine vide ou inexistante, la variable de suivi peut passer à une valeur initiale dans la plage d'indice de la maille. Par exemple, en reference à 3058447 2016-1PM- 00753-1.1 - l'exemple de la figure 2, pour une première itération. la valeur initiale de la variable de suivi peut être définie sur 201. Une fois les opérations au bloc 301 terminées, la prochaine devient la maille décrite ci-dessous pour les opérations aux blocs 302 a 320. 100331 Au bloc 302, on détermine si la maille est entourée par des mailles dominées par 5 la convection. Dans certains modes de réalisation. une maille peut. être entourée par des mailles voisines. Les mailles voisines peuvent être identifiées comme étant dans un régime d'écoulement dominé par la convection (c'est-à-dire identifiée comme des mailles dominées par la convection). Par exemple, une maille voisine peut inclure une variable de marquage qui marque la !n'aille voisine en tant que maille dominée par la convection. Alternativement, ou de plus, une maille voisine peut être définie comme ayant un nombre de Péclet élevé (c'est-à-dire supérieur au seuil de transition) et être définie comme une maille dominée par la convection. Dans certains modes de réalisation, chaque flux de diffusion à travers une surface d'une maille dominée par la convection est négligeable par rapport i tout flux de convection à travers la surface de la maille dominée par la 15 convection. Dans certains modes de réalisation, si une maille est entourée par des mailles dominées par la convection, on suppose que toutes les faces de cette maille ont un flux de diffusion nul et peuvent rester inchangées pour l'opération ce cours. Si toutes les faces de cette maille ont un flux de diffusion nul, il n'est pas nécessaire de calculer le flux de diffusion et les coefficients initiaux de diffusion de dilution infinie peuvent rester 20 inchangés, dans lequel un coefficient (le diffusion de dilution infinie d'un composant i est défini comme un coefficient de diffusion du composant i dilué à l'infini dans un deuxième gaz de composition. 100341 Par exemple, en référence à la figure 2, la maille de moindre diffusio comporte quatre mailles voisines. Les quatre mailles voisines sont les mailles (le moindre 25 diffusion 21 S. 223,126 et . Si chacune de ces mailles de moindre diffusion est marquée comme étant une maille dominée par la convection, alors l'opération peut déterminer que la maille de moindre diffusion 222 est déterminée comme étant entourée par des mailles dominées par la convection. Si la maille est entourée de mailles dominées par la convection, les opérations du diagramme de flux 300 se poursuivent au bloc 314 30 (qui est décrit plus en détail ci-dessous). Si la maille n'est pas entourée de mailles dominées par la convection, les opérations du diagramme de flux 300 se poursuivent au bloc 303. 3058447 201 6-1PM- I 00753-U - 100351 Au bloc 3Ct 3, on détermine les coefficients de diffusion de dilution infinie de la maille. Dans certains modes de réalisation, les coefficients de diffusion de dilution infinie peuvent etre déterminés par référence croisée entre cfes composés et des mélanoes connus dans le système de réservoir avec une table de données. Un variante, ou de plus. les 5 coefficients de diffusion de dilution infinie peuvent être déterminés en combinant des équations d'état avec des simulations de fluide expérimental. Dans certains modes de réalisation, les coefficients de diffusion de dilution infinie peuvent être déterminés en les rendant égaux à une valeur initiale, telle qu'un ou zéro. 100361 Le flux de diffusion peut être déterminé à l'aide de l'équation 1 basée sur les coefficients de diffusion SM auxiliaires, , où est la porosité. l'indice a désigne une phase, S, est la saturation de phase, c est la densité molaire, xi est la composition du composant i, tif est le potentiel chimique du composant 1, R est la constante de gaz universelle, et Test la température absolue de la maille : Pa,/ ,n, - 100371 Dans certains modes de réalisation, des paramètres tels que la porosité ou la 15 température peuvent être des paramètres géologiques. Un paramètre géologique peut être lié à une propriété de la formation des substrats géologiques, des roches crou des hydrocarbures. Un paramètre géologique peut être déterminé en fonction des résultats d'un tableau de données, en application de relations empiriques entre des quantités mesurables et des propriétés géologiques, des données de capteurs, des simulations. 20 100381 Dans certains modes de réalisation, des paramètres tels que la densité molaire ou la saturation en phase au niveau d'une maille d'injection ou d'une face de maille peuvent être des paramètres opérationnels. Les plinttuën-es opérationnels peuvent être déterminés en fonction d'une donnée d'entrée qui correspond à une opération physique dans un puits. Par exemple. la densité molaire au niveau d'une maille qui correspond au site d'un puits 25 d'injection peut être déterminée en se basant sur une condition limite de flux de uaz établie pour un débit. d'injection de gaz au niveau du puits d'injection. Dans certains modes de réalisation, un paramètre peut être à la fois un paramètre opérationnel et un paramètre géologique. Par exemple, une maille petit être définie par sa température initiale et sa saturation de phase sur la base des données de paramètres géologiques. La 3058447 2016-1PM-100753-U - température et la saturation de phase peuvent changer en raison d'un flux de gal d'injection au niveau de la maille pour simuler une opération d'injection. Les paramètres opérationnels peuvent être basés SUI: les réglages des pompes ou d'autres équipements dans un réservoir. Alternativement, ou de plus, les paramètres opérationnels peuvent être 5 déterminés en fonction des résultats de la corrélation des paramètres avec une LOS appropriée, des données de capteur, des paramètres bases sur les résultats d'une simulation, etc. 100391 Dans certains modes de réalisation, les éléments dc la matrice complète des coefficients de diffusion SM auxiliaires peuvent etre déterminés à l'aide des équations 2 et 3, où Bisti est un coefficient de diffusion SM auxiliaire hors diaizonale, Br est un coefficient de diffusion SM auxiliaire diagonal, n, est le nombre de composants, Do est art coefficient de diffusion SM du composant i diltta dans le composant est la composition du composant f. et Di', est le coefficient de diffusion du composant i dilué dans le gaz injecté : (2) B''' 15 100401 Dans certains modes de réalisation, les coefficients dc diffusion 5M peuvent être déterminés en utilisant l'équation de Vignes généralisée pour les mélanges multi- composants en fonction des coefficients de dilution infinie. Une forme de l'équation Vignes généralisée est présentée ci-dessous dans l'équation 4 où est le coefficient de diffusion SM de la paire Lie composants i ) est le coefficient de diffusion infinie du 20 composant i dilué dans la composante j, est k coefficient de dilution infinie du composant j diluée dans le composant i, xi est la concentration du composant i et xj est a concentration du composant j : (4) L) )' 1_.)j1 I 100411 Dans certains modes de réalisation, les coefficients de diffusion SM peuvent être remplacés par d'autres coefficients provenant d'autres équations empiriques, telles que 3058447 2016-WM-100753-LJ l-hR 14 l'équation Darkcn ou par l'équation Caldwell et Babb. Dans certains modes de réalisation, les coefficients de diffusion SM peuvent etre remplacés par des coefficients de diffusion determinés par des resultals de simulation, tels que des résultats de simulation. basés sur des simulations de dynamique moléculaire. 5 100421 Au bloc 306, un flux de diffusion estimé est déterminé en utilisant une simplification des coefficients de diffusion. des paramètres géologiques et des paramètres opérationnels. Dans certains modes de réalisation, le flux de diffusion peut être déterminé en se basant en partie sur une matrice de coefficients de diffusion de Stefan-Maxwell (SM) auxiliaires. Deux prémisses peuvent être mises à profit pour simplifier les calculs des coefficients de diffusion SM auxiliaires afin de fournir un flux de diffusion estimé. Dans la première prémisse, on suppose que la diffusion dominante d'un composant est chie au gradient de concentration de ce composant seul, de sorte que seuls les éléments diagonaux de !a matrice des coefficients de diffusion peuvent être utilisés à la place de la matrice entière. seconde prémisse dépend de deux hypothèses. la première étant que le 15 coefficient de diffusion de la phase gazeuse est plus élevé que le coefficient de diffusion. de la phase liquide de sorte que l'écoulement par diffusion maximal se produit dans la phase gazeuse et la seconde étant que la composition du gaz se rapproche rapidement de celle du gaz injecté. Ces deux hypothèses peuvent conduire à supposer que la majorité de la diffitsion se produit lors d'une phase gazeuse avec une composition qui est 20 approximativement celle de la composition de gaz injecté. 100431 L'hypothèse ci-dessus permet de déterminer les coefficients de diffusion SM auxiliaires en utilisant l'approximation représentée dans l'équation 5 ci-dessous, où chaque valeur de Bir est le coefficient de diffusion SM auxiliaire pour toute valeur de i et j. D'if 1 est le coefficient de diffusion infinie du composant n dilué dans le gaz injecté, 25 les valeurs de Du sont les coefficients de diffusion infinie du composant i dilué dans le gaz injecté pour toute valeur de est la concentration du composant 3058447 100441 Les éléments de la matrice des coefficients de diffusion de dilution infinie peuvent être déterminés à partir des propriétés de l'équation d'état (FOS) découlant de la concordance des expériences des fluides. supposant que les coefficients de diffusion lilution infinie entre les gai injectés et les autres livdrocarbitres dans le pétrole sont 5 dans le même ordre de grandeur. D;,-)_ Do) dans la première rangée de l'équation 1. peut recevoir un flux de diffusion négligeable. Si les éléments non diagonaux donnent un flux de diffusion rlégligeable, le calcul du flux a l'aide des éléments non diagonaux peut être ignore. Cela peut aboutir à une approximation de l'équation 5 sous la forme de l'équation 6_ dans laquelle à la fois la matrice comprenant les coefficients de diffusion SM auxiliaires et la matrice comprenant des coefficients de diffusion de dilution infinie sont des matrices diagonales 2016-1PM-100753-U - 15 (5) D271 (Dri°71 - 2,1 D .;r-- Dr,1 0 2,l 0 (6) 11m 100451 L'équation 6 peut être utilisée pour fournir estimation de l'inverse d'un coefficient de diffusion SM au.\ iliaire. L'équation 1 peut être simplifiée en reconnaissant que lorsqu'on utilise la simplification des coefficients de diffusion SM. tous les éléments 15 de la matrice des coefficients de diffusion SM auxil lires peuvent être considérés comme étant nuls sauf dans le cas où i est égal al 100461 Au bloc 308, le nombre de Péclet est déterminé en partie en se basant sur le flux de diffusion estimé. Le nombre de l'éclet peut être défini comme un rapport du flux de convection sur le flux de diffusion estimé. Le rapport du flux de convection sur le flux de 3058447 2016-1PM-10075340 - 16 diffusion estimé peut être déterminé pour chaque maille en utilisant l'équation 7 ci-dessous Convective flux (7) Diffusivc flux 10047] Dans certains modes de réalisanon, le flux de convection peut être déterminé à partir d'une ou de plusieurs conditions limites de la maille. Par exemple, au niveau de la 5 maille 249, une condition limite pourrait définir le flux de convection du côté droit à 30 kg.. in'-s. Dans certains modes de réalisation_ les conditions limites peuvent être définies au bloc 308. Par ailleurs, la condition limite peut etre définie au bloc 306 en t'onction d'un paramètre opérationnel. Par ailleurs ou en outre, un flux de convection pourrait être déterminé lors de !a détermination de flux dans le bloc 312 ou le bloc 314. décrit plus en détail ci-dessous. 100481 Au bloc 310, on détermine si le nombre de Péclet est inférieur à un seuil de transition. Le seuil de transition peut être d'une quelconque valeur supérieure à 1. Dans certains modes de réalisation le seuil de transition peut être de cinq. Si le nombre Péclet de la maille est inférieur au seuil de transition, l'écoulement au niveau de la maille est 15 considéré connue étant dominé par la diffusion. Si le nombre de Péclet de la maille est supérieur ou égal au seuil de transition, l'écoulement au niveau de la maille est considéré comme étant dominé par la convection. Si le nombre de Péclet de la maille est inférieur au seuil de transition, des opérations du diagramme de flux 300 se poursuivent au bloc 311. Sinon, des opérations du diagramme de flux 300 se poursuivent au bloc 314. 20 10049] Au bloc 31 1, la matrice complète des coefficients de diffusion SM est déterminée en fonction des coefficients de diffusion de dilution infinie. Dans certains modes de réalisation, la matrice complète des coefficients de diffusion SM peut être déterminée en utilisant une combinaison des équations et 4. 100501 Au bloc 312, un flux de diffusion conlplet est déterminé en .fonction de la matrice 25 complète de coefficients de diffusion 5M. Le flux de diffusion complet peut être déterminé en fonction d'une matrice complète de coefficients de diffusion 5M auxiliaires. Les éléments de la matrice complète de coefficients de diffusion SM auxiliaires peuvent être déterminés à l'aide des équations 2 et 3. Le flux de diffusion complet peut ensuite être 3058447 2016-1PM-100753-U l-hR 17 déterminé en utilisant l'équation 1 en fonction des éléments de la matrice complète de coefficients de diffusion SNI auxiliaires. 100511 Au bloc 3 I-1, le calcul de flux de diffusion coin') et est contourné et la maille est marquée comme étant dominée par la convection. Dans certains modes de réalisation, les 5 calculs de Ilux de diffusion complets peuvent être contournés en ne déterminant pas le flux de diffusion basé sur la matrice complète de coefficients de diffusion SNI. Dans certains modes de réalisation, le marquage de la maille comme étant dominée par la convection peut inclure la définition d'une variable de marquage associée à la maille qui marquera la maille en tant que maille dominée par la convection. Par exemple.en référence ài la figure 2. le marquage de la maille de moindre diffusion 208 comme étant dominée par la convection peut inclure la définition d'une variable de marqune de « CONVECTION ». 100521 An bloc 316, on détermine si la maille est la dernière maille pour les calculs de flux. Dans certains modes de réalisation, la détermination qu'une maille est la dernière 15 maille pour des calculs de flux peut inclure la détermination qu'une limite de la matrice d'indexation des mailles a été atteinte. Par exemple, en référence à la figure 2, une limite de la matrice d'indexation des mailles peut être 3-P). Si la variable de suivi est 349, l'opération peut déterminer qu'une limite de la matrice d'indexation des mailles a été atteinte. Dans certains modes de réalisation, la détermination d'une maille comme étant la 20 dernière maille pour des calculs de flux peut inclure la détermination qu'une limite de maille prédéterminée a été atteinte. Par exemple, en référence à la ligure 2. la maille 220 peut être une limite de maille prédéterminée. Si la variable de suivi est 220, alors la limite de maille prédéterminée a été atteinte. Si la maille n'est pas la dernière maille pour des calculs de flux, des opérations du dia2ramme de flux 300 se poursuivent nu bloc 301. 25 Sinon, des opérations du diagramme de flux 300 se poursuivent au bloc 320. 100531 Au bloc 320, toutes les mailles sont mises iijour en résolvant les équations du bilan des matériaux (de masse) et les équations de pression. Dans certains modes de réalisation, les équations dit bilan des matériaux (de masse) et les équations de pression peuvent être résolues pour fournir des Ilux de con \'Ceii0Ildc toutes les mailles eLott 30 d'autres paramètres de réservoir. Par exemple, les équations du bilan des matériau (de masse) et les équations de pression peuvent être résolues pour fournir des paramètres tels que la température, la pression, les compositions de phase, leurs modifications, etc. Dans 3058447 2016-1PM-100753-U - 18 certains modes de réalisation_ les résultats des équations du bilan des matériaux (de masse) et des équations de pression peuvent être utilisés pour déterminer une quantité d'hydrocarbures récupérés par un ou plusieurs puits de production. Par exemple. en référence à la figure I le changement des conditions d'une maille pour refléter une 5 pluralité de réglages de paramètres opérationnels tels que le débit d'injection ou la composition de ai commandée par le système d'injection 108 peut entraîner une pluralité d'estimations des hydrocarbures récupérés à partir du puits de production 103. 100541 Dans certains modes de réalisation, les coefficients de diffusion de dilution infinie et ou les coefficients de diffusion SM peuvent être définis pour une itération future en fonction des solutions aux équations du bilan des matériaux (de masse) et aux équations de pression lors d'une procédure de mise à jour. En outre, les coefficients de diffirsion de dilution infinie et/ou les coefficients de diffitsion SM des mailles qui sont entourées par des mailles dominées par la convection peuvent être ignorés par la procédure de mise à jour. Les opérations du diagramme de flux 300 sont complètes. 15 100551 Dans certains modes de réalisation, les opérations et les résultats du diagramme de flux 300 peinent être utilisés dans des opérations d'injection de gaz réelles ou dans une planification d'emplacement de puits. Par exemple, le changement d'un paramètre opérationnel, tel qu'un débit d'injection de gaz du ga7 injecté, peut modifier les conditions limites de flux de convection de mailles. La modification des conditions limites de fluN, de 20 convection de mailles peut entraîner une récupération de pétrole estimée différente. Dans certains modes de réalisation. la corrélation entre les paramètres opérationnels et la récupération de pétrole estimée peut être utilisée pour déterminer des paratriétres opérationnels optimaux. Par exemple, des paramètres opérationnels tels que le schéma d'injection de gaz. le débit d'injection de gaz, la composition d'un gaz injecté, etc. 25 peuvent être optimisés. Le paramètre ou l'ensemble des paramètres qui entraînent une récupération de pétrole estimée maximale ou le débit de récupération de pétrole estimée le plus élevé peuvent être sélectionnés comme paramètres opérationnels optimaux pour la récupération de pétrole dans le puits. Dans certains modes de réalisation, la corrélation entre les paramètres géologiques et la récupération de pétrole estimée petit être utilisée 30 pour optimiser les emplacements au niveau (lesquels il faut forer au moins l'un d'un puits de production et d'un puits d'injection de 2az. Par exemple, la sélection de différentes mailles comme sites de puits de production ou sites de puits d'injection peut fournir 3058447 2016-1PM-100753-U - 19 différentes estimations de la récupération de pétrole en fonction des paramètres géologiques au niveau ou près des mailles. Lin site de puits de production optimal ou un site de puits d'injection optimal peut etre &tenni. le en fonction d'une récupération de pétrole estimée maximale. 5 100561 Le diagramme de flux est fourni pour faciliter la compréhension des illustrations et ne doit pas être utilisé pour limiter le champ d'application des revendications. Le diagramme de flux illustre des exemples d'opérations qui peuvent varier dans le champ d'application des revendications. Des opérations supplémentaires peuvent être effectuées ; moins d'opérations peuvent être effectuées ; les opérations peuvent etre effectuées en parallèle ; et les opérations peuvent être effectuées dans un ordre différent. lisera compris qu'au moins certains des blocs des illustrations du diagramme de flux ct/ou des diagrammes de bloc, et des combinaisons de blocs dans les illustrations du diagramme de flux clou des diagrammes de blocs peuvent Cire réalisées par le code de programme. 1.e code de programme peut être fourni à un processeur ou à un ordinateur 15 polyvalent, à un ordinateur spécialisé ou à une autre machine ou un autre appareil. prograMMable. EXEMPLE DE DISPOSITIF INFORMXFIQUE 100571 La ligure 4 illustre un exemple d'un dispositif infOrmatique selon certains modes de réalisation. Le dispositif informatique 400 représenté dans la figure 4 peut être un 20 exemple d'au moins une partie du sous-système informatique 110 illustré dans la figure 1. Le dispositif informatique 400 comprend un processeur 401 (comprenant éventuellement de multiples processeurs, de multiples coeurs, de multiples noeuds et. ou réalisant un traitement multiprocessus, etc.). Le dispositif informatique 400 une mémoire 407. La mémoire 407 peut être une mémoire de système (par ex., un ou plusieurs parmi un cache, 25 une SRAM, une DRAM. une RAM à condensateur de calage à zéro. une RAM à double transistor, une eDRAM, une EDO RAM. une DDR RAM. une EEPROM, une NRAM. une RRAM, une SONOS, une PRAM. etc.) et l'une quelconque ou plusieurs des réalisations possibles précédemment décrites du support lisible par ordinateur. 100581 Le dispositif informatique 400 comprend également un stockage permanent de 30 données 409. Le stockage permanent de données 409 peut Cire un lecteur de disque dur, tel qu'un dispositir de stockage magnétique. I.e dispositif in ii-irmatique 400 comprend également un bits 403 (par ex., PCI, ISA. PCI- uess, Flyper-franspolt4:' bus, 3058447 2(116-11)M- 100753-01-1; 20 IntiniBand".K.. bus. NuBus, etc.) et une interface de réseau 405 par ex., une interface de canal de fibres, une interface Ethernet, une interface de système d'Internet pour petit ordinateur, une interface SON T, une interface sans fil, etc.). 100591 Le dispositif informatique 400 comprend également un simulateur 411. Le 5 simulateur 411 peut réaliser des opérations de simulation, comme il est décrit ci-dessus. L'une quelconque des fonctionnalités précédemment décrites peut être partiellement (ou intégralement) réalisée dans le matériel et ou sur le processeur 401. Par exemple, la tbnetionnalite peut être réalisée avec un circuit intégré spécifique à une application, dans la logique mise en oeuvre dans le processeur 401, dans un coprocesseur sur un dispositif ou une carte périphérique, etc. En outre, des réalisations peuvent également comprendre plus ou moins de composants supplémentaires qui ne sont pas illustrés dans la figure 4 (par .'N' des cartes vidéo. des cartes audio, d'antres interfaces réseau, des dispositifs périphériques, etc.). Ise processeur 401, l'interfuee de réseau 405 el le stockage permanent c données 409 sont couplés au bus 403. Même si elle est illustrée comme étant couplée 15 au bus 403, la mémoire 407 peut être couplée au processeur 401. 100601 Comme il sera compris, des aspects de cette divulgation peuvent être concrétisés sous forme d'un système, d'un procédé ou d'un code de programme d'instructions stockés sur un ou plusieurs supports lisibles par ordinateur. Par conséquent, les aspects peuvent prendre la forme de matériel. de logiciels (y compris micrologiciel, logiciel résident, 20 microcode, etc.), ou une combinaison d'aspects logiciels et matériels qui peuvent généralement être désignés ici comme « circuit », REVENDICATIONS1. Procédé comprenant créer un modèle de diffusion pour une simulation de la récupération d'hydrocarbures à partir d'un réservoir ayant une pluralité de fractures lors de l'injection d'un gaz injecté dans la pluralité de fractures, dans lequel le réservoir est subdivisé en une pluralité de mailles dans le modèle de di rrusion. dans lequel la création du modèle de diffusion pouf une maille parmi plusieurs mailles comprend, déterminer un rapport de flux d'un flux de convection a un flux de diffusion estimé pour la ; déterminer si le rapport de flux est inférieur au seuil ou non ; et en réponse à la détermination que le rapport de flux est inférieur au seuil, la détermination d'un flux clé diffusion complet pour la maille aux lins d'inclusion dans le modèle de diffusion ; et réaliser la simulation de la récupération d'hydrocarbures à partir du réservoir sur. la base du modèle de diffusion. 2. Procédé selon la revendication 1, dans lequel la création du modèle de diffusion pour la maille comprend : déterminer le fait que la maille est entourée ou non par d'autres mailles ayant le rapport de flux supérieur au seuil ; et en réponse à la détermination que la maille est entourée d'autres mailles ayant le rapport de flux supérieur au seuil. k contournement de la détermination dit flux de diffusion complet pour la maille dans le modèle de diffusion, et l'application de valeurs de coefficient de diffusion préalablement déterminées pour la maille. 3. Procédé selon la revendication 1, dans lequel la création du modèle de diffusion pour la maille comprend : en réponse la détermination que le rapport de flux est supérieur ou égal au seuil. le contournement de la détermination du flux de diffusion 3058447 2016-1PM-100753-U - 28 complet pour la maille pour le modèle de diffusion et l'application de valeurs de coefficient de diffusion préalablement déterminées pour la maille. 5 4. Procédé selon la revendication 1, dans lequel la création du modèle de diffusion pour maille comprend la détermination du fluN (1C di Ifusion estimé, dans lequel. la détermination du flux de diffusion estimé comprend : déterminer un coefficient de diffusion de dilution infinie pour chaque paire de composants de gaz non injecté et de gaz, injecte dans la maille du volume simulé du réservoir. 5. Procédé selon la revendication 1. comprenant éttalement le forage d'au moins l'un d'un puits d'injection et d'un puits de production dans le réservoir à un endroit qui repose, au moins en partie, sur un résultat de la 15 simulation. 6. Procédé selon la revendication 1, comprenant également : récupérer les hydrocarbures d partir du réservoir par le biais d'un puits de production en fonction de l'injection du gaz injecté dans un puits 20 d'injection avec un débit basé au moins en partie sur un résultat de la simulation, et dans lequel une composition du gaz injecté est basée, au moins en partie, sur un résultat de la simulation. 7. Système comprenant : 25 un processeur ; et un support lisible par ordinateur avant un code de programme exécutable par le processeur pour amener le processeur a. créer un modèle de diffusion pour une simulation de la récupération d'hydrocarbures à partir d'un réservoir ayant 30 une pluralité de fractures lors de l'injection d'un ga7 injecté dans ln pluralité de fractures, dans lequel le réservoir est subdivisé en une pluralité de mailles dans le modèle de diffusion, dans lequel le code de programme qui amène le processeur à créer le modèle de diffusion pour une maille 3058447 2016-IPM-100753-U l-hR 29 parmi plusieurs mailles comprend le code de programme exécutable par le processeur pour faire en sorte que le processeur, détermine un rapport de flux d'un flux de convection à 5 un flux de diffusion estimé pour la maille ; détermine si le rapport de flux est inférieur à un seuil ou non ; et en réponse à une détermination que le rapport de flux est inférieur au seuil, la détermination d'un flux de diffusion complet pour la maille aux fins d'inclusion dans le modèle de diffusion ; et réalise la simulation de la recupération d'hydrocarbures à. partir du réservoir sur la base du modèle de di lfusion. 15 8. Système selon la revendication 7, dans lequel le code de prouramme exécutable par le processeur pour amener le processeur à créer le modèle de diffusion pour la maille comprend un code de programme exécutable par le processeur pour faire en sorte que le processeur détermine si la maille est entourée par d'autres mailles a. ait le rapport de flux 20 supérieur au seuil ; et en réponse à une détermination que la maille est entourée d'autres mailles ayant le rapport de flux supérieur au seuil. le contournement de la détermination du flux de diffusion complet pour la maille pour le modèle de diffusion. 25 9. Système selon la revendication 7, dans lequel le code de programme exécutable par le processeur pour amener le processeur à créer le modèle de difibsion pour la maille comprend un code de programme exécutable par le processeur pour faire en sorte que le processeur : 30 en réponse à une déten iination que le rapport de flux est supérieur ou égal u se[ e contournement de la détermination du flux de diffUsion complet pour la maille pour le modèle de diffusion, et 3058447 2016-1PM-100753-U - 30 l'application de valeurs de coefficient de diffusion préalablement déterminées pour la maillc. 10. Système de la revendication 7, comprenant également 5 une pompe d'injection pour pomper le gaz injecté à un débit d'un puits d'injection pour produire des hydrocarbures a partir du réservoir, dans lequel le débit est basé, au moins en partie, sur un résultat de la simulation, et dans lequel une composition du Luiz injecté à pomper par la pompe d'injection est basée, au moins en partie, sur un, résultat de la simulation.