FR3058447A1

FR3058447A1 - Inclusion de diffusion selective pour une simulation de reservoir pour la recuperation des hydrocarbures

Info

Publication number: FR3058447A1
Application number: FR1759381A
Authority: FR
Inventors: Saeedeh Mohebbinia; Terry Wayne Wong
Original assignee: Landmark Graphics Corp
Current assignee: Landmark Graphics Corp
Priority date: 2016-11-08
Filing date: 2017-10-06
Publication date: 2018-05-11

Abstract

Un procédé comprend la création d'un modèle de diffusion pour une simulation de 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 diffusion. La création du modèle de diffusion pour une maille parmi plusieurs mailles comprend la détermination d'un rapport de flux d'un flux de convection à un flux de diffusion estimé pour la maille, la détermination du fait que le rapport de flux est inférieur ou non à un seuil, et en réponse à la 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. Le procédé comprend également la réalisation de la simulation de la récupération d'hydrocarbures à partir du réservoir sur la base du modèle de diffusion.

Description

DESCRIPTION [0009J La description qui suit comprend des exemples de systèmes, de procédés, de techniques et de flux de programme qui concrétisent des aspects de la divulgation. Cependant, il est entendu que cette divulgation peut être 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.

[00101 Divers modes de réalisation se rapportent à des simulations de récupération de 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, La diffusion moléculaire peut jouer un rôle important dans la récupération des hydrocarbures à partir des réservoirs fractionnés. La diffusion moléculaire peut comprendre tout mouvement thermique de particules (y compris le liquide et le gaz) à 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 être engendré en raison des gradients 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 peut entraîner une sous-estimation de la récupération d'hydrocarbures. Sans la diffusion, le gaz injecté s'écoule principalement à travers les 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

2016-IPM-100753-U1-FR 3 peuvent interagir entre elles au travers de leurs interfaces de maille au cours de la réalisation des calculs de diffusion.

[00111 Les modèles existants pour la diffusion sont soit des modèles à coefficient constant simples qui ne peuvent appréhender correctement les éléments physiques, soit 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.

[001.2J Les calculs de diffusion peuvent impliquer de nombreuses opérations matricielles, comprenant les calculs des coefficients de diffusion, des flux et des dérivés 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 celtes 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 hétérogènes avec différentes propriétés de matrice et de fracture, il peut exister des zones avec des régimes d’écoulement différents allant des régimes dominés par la diffusion (principalement proches des fractures) jusqu'aux régimes dominés par la convection (éloignés des fractures). La réalisation des calculs de diffusion pour les mailles à 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 au cours de la simulation.

[0013J 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 diffusion négligeable en comparaison à la convection. La zone la plus active au niveau de laquelle le transfert massique par diffusion maximal se produit peut se situer autour du front de fluide, où les forces d'entraînement sont à 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

2016-IPM-100753-U1-FR 4 transition, dans laquelle la composition passe de la composition du fluide injecté à celle du fluide en place. L’étendue de la zone 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 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 diffusion devient insignifiant. Par conséquent, l'extension de la zone de transition peut être corrélée au rapport du flux de convection sur le flux de diffusion, ce qui représente la définition d'un nombre de Péclet.

[0014] 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éterminés 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 de récupération d’hydrocarbures peuvent être basés sur le résultat de la simulation de réservoir.

Exemple de système [0015] La 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 1 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 forage. En outre, bien qu’il soit représenté avec un seul puits 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.

[0016] 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 production 103 ou à d'autres emplacements. Le 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

2016-IPM-100753-U1-FR 5 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 differentes, 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, [0017J La zone souterraine 104 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 granité, 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. La zone souterraine 104 peut comprendre des formations étanches au gaz qui comprennent des roches à faible perméabilité (par ex., le schiste, le charbon, ou autres).

[0018J Le système de puits 100 illustré dans la figure 1 comprend un système d'injection 108. Le système d'injection 108 peut être utilisé pour réaliser un traitement d'injection par 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 être appliqué au niveau d'un emplacement d'injection unique ou à de multiples emplacements d'injection dans une 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 de multiples emplacements différents d'injection de gaz dans un seul puits de forage, de multiples emplacements d'injection de gaz dans de multiples puits de forage différents, ou une quelconque combinaison 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 forage inclinés, des puits de forage horizontaux, des puits de forage courbés, ou des combinaisons de ceux-ci ou d'autres.

[0019} En réponse à l'injection de gaz de traitement dans le puits d'injection 102 et en 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

2016-IPM-100753-U1-FR 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 provenant 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 stockage, etc., pour un nouveau traitement, stockage et transport de fluides de réservoir récupérés du puits de production 103. En 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.

10020) Le gaz peut être fourni à partir d'un, camion muni d’un compresseur, ou à partir 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 d’injection 112 installé dans le puits d'injection 102. Le conduit de production 162 et le conduit d’injection 112 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 du puits d'injection. 102, sans cuvelage. Le 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.

(0021j Le système d'injection 108 peut également comprendre des capteurs en surface et en. fond de puits pour mesurer la pression, le débit, la température ou d'autres paramètres de traitement ou de production. Par exemple, le système d'injection 108 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 108 peut également comprendre des commandes de pompe et/ou d'autres types de commandes pour démarrer, arrêter, augmenter, diminuer ou commander 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

2016-IPM-100753-U1-FR 7 communiquer avec l'équipement afin, de surveiller et commander le traitement par injection.

[0022j La zone souterraine 104 peut inclure les fractures naturelles 151 à 156. En variante ou en plus, un traitement d’injection de fluide peut également créer des fractures dans la zone souterraine 104, ou stimuler davantage les fractures naturelles 151 à 156. 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 être 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 nonfraeturée, ou les deux. Les fractures peuvent également être reliées au système de production et peuvent constituer le conduit principal pour la production à partir du réservoir vers les puits de forage de production.

[0023J Dans certaines réalisations, le sous-système informatique 110 peut exécuter des 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 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 soussystème informatique 110. Comme autre exemple, le sous-système de commande de traitement d'injection peut modifier, actualiser ou générer un plan d'injection de gaz basé sur des simulations réalisées par le sous-système informatique 110 en temps réel au cours 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 nouveaux puits.

[0024j 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 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 compositionnelle en se basant sur les données

2016-IPM-100753-U1-FR 8 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 de gaz, des densités de gaz, des compositions de gaz, et d'autres paramètres de commande en se basant sur les données fournies par les simulations. Dans certains cas, les données fournies par les simulations peuvent être affichées en temps réel au cours du traitement d'injection, par ex., à un ingénieur ou à un autre exploitant du système de puits 100.

10025j 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 de plus grande diffusion 241. à 249.

10026J Chacune des rangées du système de mailles 200 comprend sept mailles. De haut en bas dans l'ordre croissant, les rangées du système de mailles 200 peuvent être 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 200 comprend les mailles de moindre diffusion 208 à 209, les mailles de diffusion moyenne 228 à 230 et les mailles de moindre diffusion 210 à 211. La troisième rangée du système de mailles 200 comprend les mailles de moindre diffusion 212 à 213, la maille de diffusion 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 2.34 à 235, les mailles de plus grande diffusion 242 à 244 et la maille de moindre diffusion 216. La cinquième rangée du système de mailles 200 comprend les mailles de moindre diffusion 217 à. 219, la maille de diffusion moyenne 236, les mailles de plus grande diffusion 245à 246 et la maille de moindre diffusion 220. La sixième rangée du système de mailles 200 comprend les mailles de moindre diffusion 221 à 224, la maille de diffusion moyenne 237 et les mailles de plus grande diffusion 247 à 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.

2016-IPM-100753-U1-FR 9 [0027J 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 à une valeur fixe du flux. Le flux de diffusion des mailles de moindre diffusion 201 à 227 peut être inférieur au flux de diffusion de l’une quelconque des mailles de diffusion moyenne 228 à 240 ou des mailles de plus grande diffusion 241 à 249. 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.

[0028j 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 (kg/m²-s). Le flux de diffusion de chacune des mailles de diffusion moyenne 228 à 240 peut être dans une plage de 10 à moins de 100 kg/m²-s. Le flux de diffusion de chacune des mailles de plus grande diffusion 241 à 249 peut être un flux égal ou supérieur à 100kg/m²-s.

(0029J Dans certains modes de réalisation, une maille peut être représentée en deux dimensions en estimant la troisième dimension (par exemple, la hauteur) pour être égale pour chaque maille. Dans de tels modes de réalisation bidimensionnels, chacune des mailles peut être rectangulaire et de dimensions égales. 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 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 mailles dodécaédriques, etc.

2016-IPM-100753-U1-FR 10

Exemple d'opérations (0030J Des opérations sont maintenant décrites pour réduire le temps de calcul de la diffusion molé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.

(0031J La figure 3 illustre un schéma fonctionnel pour contourner les calculs de diffusion dans une simulation de diffusion moléculaire pour les mailles d’un réservoir qui sont dominées par le flux de convection, selon certains modes de 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 peut ê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 300 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 dispositif informatique représenté sur la figure 4 (décrit plus en détail ci-dessous), un processeur peut exécuter des instructions pour effectuer les opérations du diagramme de 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.

[0032j 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. Une variable de suivi peut être utilisée pour déterminer quelle maille est une prochaine maille. Par exemple, la simulation de réservoir 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 d'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ées à la maille correspondante à cette valeur d'indice sui vante. 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 être é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 référence à

2016-IPM-100753-U1-FR 11 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 maille devient la maille décrite ci-dessous pour les opérations aux blocs 302 à 320.

10033j Au bloc 302, on. détermine si la maille est entourée par des mailles dominées par 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 maille 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 à tout flux de convection à travers la surface de la maille dominée par la 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, en 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 inchangés, dans lequel un coefficient de 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.

[0034j Par exemple, en référence à la figure 2, la maille de moindre diffusion 222 comporte quatre mailles voisines. .Les quatre mailles voisines sont les mailles de moindre diffusion 218,223, 226 et 221. 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 (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.

2016-IPM-100753-U1-FR [0035} Au bloc 303, 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 être déterminés par référence croisée entre des composés et des mélanges connus dans le système de réservoir avec une table de données. En variante, ou de plus, les 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.

[0036} Le flux de diffusion peut être déterminé à l’aide de l'équation 1 basée sur les coefficients de diffusion SM auxiliaires, B^, où φ est la porosité, l'indice a désigne une phase, S_a est la saturation de phase, c est la densité molaire, Xj est la composition du composant i, est le potentiel chimique du composant i, R est la constante de gaz. universelle, et T est la température absolue de la maille :

Cl)

[0037} 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 et/ou 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, etc.

[0038} 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 paramètres 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 d'injection peut être déterminée en se basant sur une condition limite de flux de gaz é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 peut ê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

2016-IPM-100753-U1-FR 13 température et la saturation de phase peuvent changer en raison d’un flux de gaz d'injection au niveau de la maille pour simuler une opération d’injection. Les paramètres opérationnels peuvent être basés sur les réglages des pompes ou d'autres équipements dans un réservoir. Alternativement, ou de plus, les paramètres opérationnels peuvent être déterminés en fonction des résultats de la corrélation des paramètres avec une EOS appropriée, des données de capteur, des paramètres basés sur les résultats d'une simulation, etc.

(00391 Dans certains modes de réalisation, les éléments de la matrice complète des coefficients de diffusion SM auxiliaires peuvent être déterminés à l'aide des équations 2 et

3, où B™ est un coefficient de diffusion SM auxiliaire hors diagonale, B™ est un coefficient de diffusion SM auxiliaire diagonal, n_c est le nombre de composants, EU- est un coefficient de diffusion SM du composant i dilué dans le composant k, xt est la composition du composant k, et D_inc est le coefficient de diffusion du composant i dilué dans le gaz injecté :

(2)	nSU _ ^Xi o_ti — — D.
(3)	_Dsu _ ,, z 1

Σ* i=i _ni*k U.

îk f0040j Dans certains modes de réalisation, les coefficients de diffusion SM peuvent être déterminés en utilisant l'équation de Vignes généralisée pour les mélanges multicomposants 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ù D_Sj est le coefficient de diffusion SM de la paire de composants i-j, est le coefficient de diffusion infinie du composant i dilué dans la composante j, Dp est le coefficient de dilution infinie du composant] diluée dans le composant i, xi est la concentration du composant i et xj est la concentration du composant j :

(4) n:., «çcsr’ ^V k*l,j

10041j 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

2016-IPM-100753-U1-FR. 14 l’équation Darken ou par l'équation. Caldwell et Babb. Dans certains modes de réalisation, les coefficients de diffusion SM peuvent être remplacés par des coefficients de diffusion déterminés par des résultats de simulation, tels que des résultats de simulation, basés sur des simulations de dynamique moléculaire, [00421 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 S.M 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 due au gradient de concentration de ce composant seul, de sorte que seuls les éléments diagonaux de la matrice des coefficients de diffusion peuvent être utilisés à la place de la matrice entière. La seconde prémisse dépend de deux hypothèses, la première étant que le 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 diffusion se produit lors d'une phase gazeuse avec une composition qui est approximativement celle de la composition de gaz injecté.

[0043J 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 Bff est le coefficient de diffusion SM auxiliaire pour toute valeur de i et j, est le coefficient de diffusion infinie du composant n dilué dans le gaz injecté, les valeurs de D“_r sont les coefficients de diffusion infinie du composant i dilué dans le gaz injecté pour toute valeur de i, et x, est la concentration du composant /' :

2016-IPM-100753-U1-FR 15 (5) nSM °ii

ΙΐϊΤΙ

Xi-41.0 n$M ^ϋ1,η-1

-,-1 — dSM dSM r<i - «i-<d rjCO ^2,1 n°° ^υη-2,1 n^œ ^un-l,l [0044J 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 (EOS) découlant de la concordance des expériences des fluides. En supposant que les coefficients de diffusion de dilution infinie entre les gaz injectés et les autres hydrocarbures dans le pétrole sont dans te même ordre de grandeur, (D“i — D_t“ ) 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 négligeable, le calcul du flux à l'aide des éléments non diagonaux peut être ignoré. 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 :

D™ (6) lim

SM

-i ri,i

CO n°° ^υη-2,ΐ n® [0045J L'équation 6 peut être utilisée pour fournir une estimation de t'inverse d'un coefficient de diffusion SM auxiliaire. L'équation 1 peut être simplifiée en reconnaissant que lorsqu'on utilise la simplification des coefficients de diffusion SM, tous les éléments de la matrice des coefficients de diffusion SM auxiliaires peuvent être considérés comme étant nuis sauf dans le cas où ï est égal àj.

[0046> 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 Pé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

2016-IPM-100753-U1-FR 16 diffusion estimé peut être déterminé pour chaque maille en. utilisant l'équation 7 cidessous :

_Kî __ Convective flux Difflisive flux [0047J Dans certains modes de réalisation, 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 maille 249, une condition limite pourrait définir le flux de convection du côté droit à 30 kg/m²-s. Dans certains modes de réalisation, les conditions limites peuvent être définies au bloc 308. Par ailleurs, la condition limite peut être définie au bloc 306 en fonction d'un paramètre opérationnel. Par ailleurs ou en outre, un flux de convection pourrait être déterminé lors de la détermination de flux dans le bloc 312 ou le bloc 314, décrit plus en détail ci-dessous.

[00481 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 considéré comme é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.

[0049J .Au bloc 311,1a 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 2, 3 et 4.

[0050J Au bloc 312, un. flux de diffusion complet est déterminé en fonction de la matrice complète de coefficients de diffusion SM, Le flux de diffusion complet peut être déterminé en fonction d'une matrice complète de coefficients de diffusion SM 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

2016-IPM-100753-U1-FR Π déterminé en utilisant l'équation 1 en fonction des éléments de la matrice complète de coefficients de diffusion SM auxiliaires.

[OOSIJ Au bloc 314, le calcul de flux de diffusion complet est contourné et la maille est marquée comme étant dominée par la convection. Dans certains modes de réalisation, les calculs de flux 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 SM, 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 à 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 marquage de « CONVECTION ».

[0052j Au 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 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 349, 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 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 figure 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 diagramme de flux 300 se poursuivent au bloc 301. Sinon, des opérations du diagramme de flux 300 se poursuivent au bloc 320.

[0053j Au bloc 320, toutes les mailles sont mises à jour 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 du bilan des matériaux (de masse) et les équations de pression peuvent être résolues pour fournir des flux de convection de toutes les mailles et/ou d'autres paramètres de réservoir. Par exemple, les équations du bilan des matériaux (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

2016-IPM-100753-U1-FR 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 1, le changement des conditions d’une maille pour refléter une pluralité de réglages de paramètres opérationnels tels que le débit d’injection ou la composition de gaz 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.

(0054] 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 diffusion de dilution infinie et/ou les coefficients de diffusion 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.

|0055] Dans certains modes de réalisation, les opérations et les résultats du diagramme de flux 300 peuvent ê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. gaz injecté, peut modifier les conditions limites de flux de convection de mailles. La modification des conditions limites de flux de convection de mailles peut entraîner une récupération de pétrole estimée differente. 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 paramè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. 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 peut être utilisée pour optimiser les emplacements au niveau desquels il. faut forer au moins l’un d’un puits de production et d’un puits d'injection de gaz. Par exemple, la sélection de différentes mailles comme sites de puits de production ou sites de puits d'injection peut fournir

2016-IPM-100753-U1-FR 19 differentes estimations de la récupération de pétrole en fonction des paramètres géologiques au niveau ou près des mailles. Un site de puits de production optimal ou un site de puits d'injection optimal peut être déterminé en fonction d'une récupération de pétrole estimée maximale.

[0056J 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 être effectuées en parallèle ; et les opérations peuvent être effectuées dans un ordre différent.

Il sera compris qu'au moins certains des blocs des illustrations du diagramme de flux et/ou des diagrammes de bloc, et des combinaisons de blocs dans les illustrations du diagramme de flux et/ou des diagrammes de blocs peuvent être réalisées par le code de programme. Le code de programme peut être fourni à un processeur ou à un ordinateur polyvalent, à un ordinateur spécialisé ou à une autre machine ou un autre appareil programmable.

ΙΑΓΜΙΊ.Ι. DF DISPOSITIF INFORMATIQUE (00571 La figure 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 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 cœurs, de multiples nœuds 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, 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 N'RAM, 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.

[0058J Le dispositif informatique 400 comprend également un stockage permanent de données 409. Le stockage permanent de données 409 peut être un lecteur de disque dur, tel qu’un dispositif de stockage magnétique. Le dispositif informatique 400 comprend également un bus 403 (par ex., PCI, ISA, PCI-Express, HyperTransport® bus,

2016-IPM-100753-U1-FR. 20

InfiniBandK 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 SONET, une interface sans fil, etc.),

I0Ü59J Le dispositif informatique 400 comprend également un simulateur 411. Le 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 fonctionnalité peut être réalisée avec un circuit intégré spécifique à une application, dans la logique mise en œuvre 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 ex,, des cartes vidéo, des cartes audio, d'autres interfaces réseau, des dispositifs périphériques, etc.). Le processeur 401, l’interface de réseau 405 et le stockage permanent de données 409 sont couplés au bus 403. 'Même si elle est illustrée comme étant couplée au bus 403, la mémoire 407 peut être couplée au processeur 401.

[0O6OJ 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, microcode, etc,), ou une combinaison d'aspects logiciels et matériels qui peuvent généralement être désignés ici comme « circuit », « module » ou « système ». La fonctionnalité présentée sous forme de modules/d’unités individuels dans les exemples d'illustrations peut être organisée différemment selon une quelconque plateforme (système d'exploitation et/ou matériel), écosystème d'application, interfaces, préférences de programmeur, langue de programmation, préférences d'administrateur, etc.

[0061 j Une quelconque combinaison d'un ou de plusieurs supports lisibles par ordinateur peut être utilisée ici. Le support lisible par ordinateur peut être un support de signal lisible par ordinateur ou un support de stockage lisible par ordinateur. Un support de stockage lisible par ordinateur peut être, par ex., mais sans limitation, un système, un appareil ou un dispositif) qui utilise une quelconque technologie ou combinaison de technologie électronique, magnétique, optique, électromagnétique, infrarouge ou semi-conductrice pour stocker un code de programme. D'autres exemples plus spécifiques (une liste non

2016-IPM-100753-U1-FR. 21 exhaustive) du support de stockage lisible par ordinateur comprendrait les éléments suivants : une disquette d’ordinateur portable, un disque dur, une mémoire vive (RAM), une mémoire morte (ROM), une mémoire morte programmable effaçable (EPROM ou une mémoire Flash), une mémoire morte sur disque compact à lecture seule (CD-ROM), un dispositif de stockage optique, un dispositif de stockage magnétique, ou une quelconque combinaison appropriée des éléments précédents. Dans le contexte de ce document, un support de stockage lisible par ordinateur peut être un quelconque support tangible qui peut contenir ou stocker un programme pour une utilisation par ou en relation avec un système, un appareil ou un dispositif d'exécution d'instructions. Un support de stockage lisible par ordinateur n'est pas un support de signal lisible par ordinateur.

[0062J Un support de signal lisible par ordinateur peut comprendre un signal de données propagé avec un code de programme lisible par ordinateur inscrit dans celui-ci, par ex., en bande de base ou en tant qu’onde porteuse. Un tel. signal propagé peut prendre une quelconque diversité de formes, comprenant, mais sans limitation, une forme électromagnétique, optique, ou une quelconque combinaison appropriée de celles-ci. Un. support de signal lisible par ordinateur peut être un quelconque support lisible par ordinateur qui n'est pas un support de stockage lisible par ordinateur et qui peut communiquer, propager ou transporter un. programme pour être utilisé par ou en relation avec un. système d’exécution d’instructions, un appareil, ou un. dispositif.

10063) Le code de programme inscrit sur un support lisible par ordinateur peut être transmis en utilisant un quelconque support approprié, comprenant, mais sans limitation, le Wi-Fi, une ligne câblée, un câble de fibre optique, la RF, etc., ou une quelconque combinaison appropriée des éléments précédents. Le code de programme informatique permettant de réaliser les opérations pour les aspects de cette divulgation peut être écrit dans une quelconque combinaison d'un ou de plusieurs langages de programmation, comprenant un langage de programmation orienté objet tel que le langage de programmation Java®, C++, ou semblable ; un langage de programmation, dynamique tel que Python ; un. langage de script tel que le langage de programmation Perl ou le langage de script PowerShell ; et des langages de programmation, procéduraux classiques, tels que le langage « C » ou des langages de programmation semblables. Le code de programme peut s'exécuter entièrement sur un ordinateur autonome, peut s'exécuter d'une façon distribuée à travers de multiples ordinateurs et peut s'exécuter sur un ordinateur tout en

2016-IPM-100753-U1-FR 22 fournissant des résultats ou en. acceptant des données d'entrées provenant d'une autre machine.

[0064} Le codes de programme/instructions peuvent également être stockés dans un support lisible par ordinateur qui peut ordonner à une machine de fonctionner d’une façon particulière, de sorte que les instructions stockées dans le support lisible par ordinateur produisent un article de fabrication comprenant des instructions qui réalisent la fonction/l'action mentionnée dans le bloc ou les blocs du diagramme de flux et/ou du diagramme en blocs.

[0065} Sachant que les aspects de cette divulgation sont décrits en référence à diverses réalisations et exploitations, il sera compris que ces aspects sont de nature illustrative et que le cadre des revendications n'est pas limitée par ceux-ci. De nombreuses variations, modifications, additions et améliorations sont possibles.

10066} Plusieurs cas peuvent être fournis pour les composants, opérations ou structures décrits ici sous forme de cas unique. Finalement, les limites entre les divers composants, opérations et stockages des données sont quelque peu arbitraires, et des opérations particulières sont illustrées dans le contexte de configurations illustratives spécifiques. D’autres attributions de fonctionnalités sont envisagées et peuvent entrer dans le champ d’application de la divulgation. En. règle générale, les structures et les fonctionnalités présentées sous forme de composants distincts dans les exemples de configurations peuvent être réalisées sous forme d'une structure ou d'un composant combiné. De la même façon, les structures et les fonctionnalités présentées sous forme d'un composant unique peuvent être réalisées sous forme de composants distincts. Celles-ci et d’autres variations, modifications, additions et améliorations peuvent entrer dans le champ d’application de la divulgation.

Exemple de modes de réalisation [0067} Les exemples de réalisation comprennent les éléments suivants :

[0068} Mode de réalisation 1 : Un 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 diffusion.

2016-IPM-100753-U1-FR. 23 dans lequel la création du modèle de diffusion pour une maille parmi une pluralité de mailles comprend, la détermination d'un rapport de flux de convection sur un flux de diffusion estimé pour la maille ; la détermination, du fait que le rapport de flux est inférieur à un seuil ; et en réponse à la 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 la réalisation de la simulation de la récupération d'hydrocarbures à partir du réservoir sur la base du modèle de diffusion.

10069} Mode de réalisation 2 : Le procédé selon le mode de réalisation 1, dans lequel la création du modèle de diffusion pour la maille comprend : déterminer si la maille est entourée 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 par d'autres mailles ayant le rapport de flux supérieur au seuil, contourner la détermination du flux de diffusion complet pour la maille pour le modèle de diffusion.

[0070} Mode de réalisation 3 : Le procédé selon les modes de réalisation 1 ou 2, dans lequel, en réponse à la détermination que la maille est entourée par d'autres mailles ayant le rapport de flux supérieur au seuil, l’application de valeurs de coefficient de diffusion préalablement déterminées pour la maille.

[0071J Mode de réalisation 4 : Le procédé selon les modes de réalisation 1 à 3, dans lequel la création du modèle de diffusion pour les mailles 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 complet pour la maille pour le modèle de diffusion.

[0072j Mode de réalisation 5 : Le procédé selon l’un quelconque des modes de réalisation 1 à 4, dans lequel, en. réponse à la détermination que le rapport de flux est supérieur ou égal au seuil, l'application de valeurs de coefficient de diffusion préalablement déterminées pour la maille.

[0073} Mode de réalisation 6 : Le procédé selon l'un quelconque des modes de réalisation 1 à 5, dans lequel la création du modèle de diffusion comprend la création d'un modèle de diffusion moléculaire.

[0074} Mode de réalisation 7 : Le procédé selon l'un quelconque des modes de réalisation 1 à 6, dans lequel la création du modèle de diffusion pour la maille comprend

2016-IPM-100753-U1-FR 24 la détermination du flux de diffusion, estimé, dans lequel la détermination du flux de diffusion estimé comprend : la détermination d'un coefficient de diffusion de dilution infinie pour chaque paire de composants de gaz non injecté et de gaz injecté dans la maille du volume simulé du réservoir.

[0075J Mode de réalisation. 8 : Le procédé selon l'un quelconque des modes de réalisation 1 à 7, dans lequel la détermination du flux de diffusion estimé comprend : l'établissement de la différence deux par deux entre les coefficients de diffusion de dilution infinie déterminés pour chaque paire de gaz non injecté et de gaz injecté à zéro.

(00761 Mode de réalisation 9 : Le procédé selon l’un quelconque des modes de réalisation 1 à 8, comprenant en outre : 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 simulation.

(0077J Mode de réalisation 10 : Le procédé selon l'un quelconque des modes de réalisation 1 à 9, comprenant en outre : la récupération d’hydrocarbures à partir du réservoir via un puits de production en fonction de l'injection du gaz injecté dans un puits d'injection à un. débit qui est basé, au moins en partie, sur un résultat de la simulation.

[00781 Mode de réalisation 11 : Le procédé selon l’un, quelconque des modes de réalisation 1 à 10, comprenant en outre : la récupération d'hydrocarbures à partir du réservoir via un puits de production, en. fonction de l’injection du gaz injecté dans un. puits d'injection, dans lequel une composition du gaz injecté est basée, au moins en partie, sur un résultat de la simulation.

[00791 Mode de réalisation. 12 : Un ou plusieurs supports de stockage non. transitoires lisibles par ordinateur comprenant un code the programme pour une simulation, le code de programme permettant : de 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 dans la pluralité de fractures, dans lequel Îe réservoir est subdivisé en une pluralité de mailles dans le modèle de diffusion, dans lequel le code de programme pour créer le modèle de diffusion pour une maille de la pluralité de mailles comprend le code de programme pour, déterminer si la maille est entourée par d'autres mailles ayant un rapport de flux de convection sur un flux de diffusion estimé supérieur à un. seuil ; en réponse à une détermination selon laquelle la maille est entourée par d'autres mailles

2016-IPM-100753-U1-FR 25 ayant le rapport de flux supérieur au seuil, contourner la détermination d'un flux de diffusion complet de la maille pour le modèle de diffusion ; déterminer le rapport de flux pour la maille ; déterminer si le rapport de flux pour la maille est inférieur au seuil ; et en. réponse à une détermination selon laquelle le rapport de flux pour la maille est inférieur au seuil, déterminer le flux de diffusion complet pour la maille aux fins d'inclusion dans le modèle de diffusion ; et effectuer la simulation de la récupération d’hydrocarbures à partir du réservoir en fonction du modèle de diffusion.

[0080J Mode de réalisation 13 : Un ou plusieurs supports de stockage non transitoires lisibles par ordinateur selon le mode de réalisation 12, dans lequel le code de programme comprend un code de programme pour : en réponse à la détermination que la maille est entourée par d'autres mailles ayant un rapport de flux supérieur au seuil, appliquer des valeurs de coefficients de diffusion préalablement déterminées pour la maille.

[008lj Mode de réalisation 14 : Un ou plusieurs supports de stockage non transitoires lisibles par ordinateur selon les modes de réalisation 12 ou 13, dans lesquels le code de programme pour créer le modèle de diffusion pour la maille comprend un code de programme pour : en réponse à la détermination que le rapport de flux est supérieur ou égal au seuil, contourner la détermination du flux de diffusion complet pour la maille pour le modèle de diffusion.

[0082j Mode de réalisation 15 : Un ou plusieurs supports de stockage non transitoires lisibles par ordinateur selon les modes de réalisation 12 à 14, dans lesquels le code de programme comprend un code de programme pour : en réponse à la détermination que le rapport de flux est supérieur ou égal au seuil, appliquer des valeurs de coefficients de diffusion préalablement déterminées pour la maille.

[0083J Mode de réalisation 16 : Un système comprenant : un processeur ; et un support lisible par machine ayant un code de programme exécutable par le processeur pour amener le processeur à créer un modèle de diffusion pour une simulation de 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 diffusion, dans lequel le code de programme permettant au processeur de créer le modèle de diffusion pour une maille de la pluralité de mailles comprend un code de programme exécutable par le processeur pour amener le

2016-IPM-100753-U1-FR. 26 processeur à déterminer un rapport de flux d'un flux de convection sur un flux de diffusion estimé pour la maille ; déterminer si le rapport de flux est inférieur à un seuil ; et en. réponse à la détermination que le rapport de flux est inférieur au seuil, déterminer un flux de diffusion complet pour la maille aux fins d'inclusion dans le modèle de diffusion ;

et effectuer la simulation de la récupération d’hydrocarbures à partir du réservoir en fonction du modèle de diffusion.

[00841 Mode de réalisation 17 : Le système selon le mode de réalisation 16, dans lequel le code de programme 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 amener le processeur : à déterminer si la maille est entourée 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 par d'autres mailles ayant le rapport de flux supérieur au seuil, contourner la détermination du flux de diffusion complet pour la maille pour le modèle de diffusion.

[0085j Mode de réalisation 18 : Le système selon les modes de réalisation 16 ou 17, dans lequel le code de programme 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 amener le processeur à : en réponse à la détermination que le rapport de flux est supérieur ou égal au seuil, contourner la détermination du flux de diffusion complet pour la maille pour le modèle de diffusion. ; et appliquer des valeurs de coefficients de diffusion préalablement déterminées pour la maille.

[0086J Mode de réalisation 19 : Le système selon l’un quelconque des modes de réalisation 16 à 18 comprenant également : une pompe d'injection pour pomper le gaz injecté à un débit dans un. puits d'injection pour produire des hydrocarbures à partir du réservoir, et dans lequel le débit est basé, au moins en partie, sur un résultat de la simulation.

[0087J Mode de réalisation 20 : Le système selon l'un quelconque des modes de réalisation 16 à 19, dans lequel une composition, du gaz injecté à pomper par la pompe d'injection est basée, au moins en partie, sur un résultat de la simulation.

2016-IPM-100753-U1-FR 27

Claims

REVENDICATIONS

1. 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 diffusion, dans lequel la création du modèle de diffusion pour une maille parmi plusieurs mailles comprend, déterminer un rapport de flux d’un flux de convection à un flux de diffusion estimé pour la maille ;

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 de diffusion complet pour la maille aux fins 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, le contournement de la détermination du 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

2016-IPM-100753-U1-FR. 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.
4. Procédé selon la revendication 1, dans lequel la création du modèle de diffusion pour la maille comprend la détermination du flux de diffusion 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 injecté dans la maille du volume simulé du réservoir.
5. Procédé selon la revendication. 1, comprenant également :

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 simulation.
6. Procédé selon la revendication 1, comprenant également :

récupérer les hydrocarbures à partir du réservoir par le biais d’un puits de production en fonction, de l'injection du gaz injecté dans un puits 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 :

un processeur ; et un support lisible par ordinateur ayant un code de programme exécutable par le processeur pour amener le processeur à, 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 diffusion, dans lequel le code de programme qui amène le processeur à créer le modèle de diffusion pour une maille

2016-IPM-100753-U1-FR 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 à 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 récupération d’hydrocarbures à partir du réservoir sur la base du modèle de diffusion.
8. 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 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 ayant le rapport de flux 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.
9. Système scion la revendication 7, dans lequel le code de programme 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 :

en réponse à une détermination que le rapport de flux est supérieur ou égal au seuil, le contournement de la détermination du flux de diffusion complet pour la maille pour le modèle de diffusion, et

2016-IPM-100753-U1-FR. 30 l’application, de valeurs de coefficient de diffusion préalablement déterminées pour la maille.
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 à 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 gaz injecté à pomper par la pompe d'injection est basée, au moins en partie, sur un. résultat de la simulation.

2016-IPM-100753-U1-FR

1/4

2016-IPM-100753-U1-FR

2/4

201 2·2 203 — γ\ρ L· 206 71 21 19 S 21 21 222 23 Îsix! X/SZA 7/77 J* 232 6 '777 y 223 Λ 21 215 . 7777/ 21 217 η 21 21 221 Μ 222 £231 Z οο* §TMa KjçsMS 225 2â 22 7//// <zz/i 7//4. ssçsæ (TC rfp M