EP2253797A1 - Méthode d'exploitation de milieu poreux au moyen d'une modélisation d'écoulements de fluide - Google Patents
Méthode d'exploitation de milieu poreux au moyen d'une modélisation d'écoulements de fluide Download PDFInfo
- Publication number
- EP2253797A1 EP2253797A1 EP10290174A EP10290174A EP2253797A1 EP 2253797 A1 EP2253797 A1 EP 2253797A1 EP 10290174 A EP10290174 A EP 10290174A EP 10290174 A EP10290174 A EP 10290174A EP 2253797 A1 EP2253797 A1 EP 2253797A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- well
- reservoir
- simulator
- mesh
- flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/12—Methods or apparatus for controlling the flow of the obtained fluid to or in wells
Abstract
Description
- La présente invention concerne le domaine de l'exploitation de milieux souterrains. L'invention permet notamment d'améliorer l'injectivité et la productivité de puits forés à travers un milieu poreux, tel qu'un gisement d'hydrocarbures ou un réservoir de stockage géologique de CO2.
- Les phénomènes locaux qui peuvent se produire autour d'un puits, tels que l'endommagement, ont un impact énorme sur l'injectivité ou la productivité d'un puits. Dans l'industrie pétrolière, il est très important de prédire l'injectivité ou la productivité, surtout lorsqu'il y a des altérations de formation au voisinage de puits, qui changent la capacité d'injection ou de production du puits.
- Depuis des années, des efforts conséquents sont effectués au moyen de techniques expérimentales, en laboratoire, ou de méthodes numériques de modélisations, afin de prendre en compte ces phénomènes locaux autour des puits, et leurs impacts sur l'injectivité ou la productivité.
- Les méthodes numériques qui permettent de modéliser l'écoulement de fluides au sein d'un puits (injectivité et productivité d'un puits), comportent la construction de deux modèles distincts : le modèle de réservoir (« reservoir model ») et le modèle d'abords de puits (« near-wellbore model »).
- Un modèle de réservoir comporte deux éléments :
- un maillage, dit « maillage de réservoir », constitué d'un ensemble de mailles discrétisant spatialement le réservoir.
- un simulateur d'écoulement. On appelle simulateur d'écoulement, un logiciel permettant de modéliser les écoulements de fluides au sein d'un milieu poreux, grâce au maillage de réservoir. Ce logiciel permet ainsi de simuler des données/propriétés dynamiques des fluides (eau, huile, gaz) : pression, flux (quantité de matière traversant une surface), saturation, débits, concentrations. Par exemple, un simulateur permet d'estimer, pour un scénario d'exploitation de puits donné (scénario de production ou scénario d'injection), et pour un intervalle de temps donné : les saturations en eau, huile et gaz, les débits d'huile, de gaz et d'eau, le water-cut (fraction d'eau dans la production liquide), le GOR (rapport de gaz et huile dans la production), les concentrations en polymère absorbé sur la roche du milieu poreux, les débits d'injection de polymère, si une solution de polymère est injectée dans le réservoir par l'intermédiaire de puits d'injection, ...
- Un modèle d'abords de puits comporte deux éléments :
- un maillage, dit « maillage d'abords de puits », constitué d'un ensemble de mailles discrétisant spatialement le puits et ses abords. Ses abords appartiennent donc au milieu poreux dans lequel le puits est foré.
- un simulateur d'écoulement, permettant de simuler, grâce au maillage d'abords de puits, des données/propriétés dynamiques des fluides (eau, huile, gaz).
- Ces deux types de modèles, réservoir et abord de puits, sont en général autonomes et découplés. Les phénomènes locaux étant généralement limités au voisinage immédiat du puits (de quelques centimètres à quelques mètres), des très petites mailles sont nécessaires pour le maillage d'abord de puits, alors que des mailles de tailles plus grandes sont utilisées pour les maillages de réservoir de façon à accélérer les calculs.
- On connaît des techniques permettant d'utiliser un seul et même simulateur d'écoulement de réservoir pour ces deux maillages. On peut par exemple utiliser la technique des maillages dits « hybrides », combinant en un seul maillage, des mailles pour le maillage du réservoir et des mailles pour un maillage localement affiné autour du puits. On associe à ce type de maillage un seul simulateur d'écoulement pour mieux prendre en compte les comportements des écoulements au voisinage du puits dans une simulation du champ.
- Mais les simulations simultanées des écoulements dans le réservoir, qui nécessitent un très grand nombre de mailles, et dans les régions proches du puits avec de petites mailles, qui nécessitent des petits pas de temps pour assurer la stabilité du calcul, posent des problèmes numériques, en particulier, le problème du temps de calcul (temps CPU).
- Ainsi, ont été mises au point des techniques de décomposition de domaine, décrites par exemple dans GAIFFE, S. "Maillages Hybrides et Décomposition de Domaine pour la Modélisation des Réservoirs Pétroliers", Thèse de Doctorat, Université Paris 6, 2000, et des techniques de fenêtrage ("windowing"), décrites par exemple dans le document suivant :MLACNIK, M.J. and HEINEMANN, Z.E. "Using well windows in full field reservoir simulation", paper SPE 66371 presented at the SPE Reservoir Simulation Symposium, Houston, TX, U.S.A., February 2001.
- Mais, certains points délicats comme la convergence, la stabilité ou le temps de calcul posent des problèmes lors d'applications industrielles. De plus, la méthode de décomposition de domaine n'est pas toujours "conservative" (détérioration du bilan de masse dans le modèle en fonction du temps), ce qui ne convient pas à une utilisation pratique de la méthode. Par ailleurs, toutes ces techniques demandent de reformuler les équations mathématiques et les conditions aux limites développées dans les simulateurs d'écoulement, et de nouveaux développements sont nécessaires pour intégrer les solutions proches et éloignées du puits dans un seul modèle, dont la réalisation est un travail long et difficile.
- Un objet de l'invention concerne une méthode mise en oeuvre par ordinateur pour modéliser des écoulements de fluides au sein d'un milieu poreux traversé par au moins un puits. La méthode comporte l'utilisation d'un premier simulateur d'écoulement permettant de simuler l'écoulement des fluides au sein du milieu poreux à partir d'indices de productivité numériques reliant des pressions de fluides à des débits de fluides, et l'utilisation d'un second simulateur d'écoulement pour simuler l'écoulement des fluides aux abords du puits à partir de conditions aux limites. La méthode comporte les étapes suivantes :
- a- on simule les écoulements de fluides au sein du milieu au moyen du premier simulateur sur un intervalle de temps défini entre des temps T0 et T1 , et on en déduit des conditions aux limites mises à jour pour le second simulateur ;
- b- on simule les écoulements de fluides aux abords du puits au moyen du second simulateur sur le même intervalle de temps, en utilisant les conditions aux limites mises à jour, et on en déduit des indices de productivité numériques mis à jour pour le premier simulateur ; et
- c- on modélise les écoulements de fluides au sein du milieu poreux pendant une période de temps entre T0 et Tn , où Tn>T1, en réitérant les étapes a et b, pour des intervalles de temps successifs compris entre T0 et Tn .
- Selon l'invention, chaque intervalle de temps successif peut avoir une longueur qui est fonction d'un pas de temps de calcul du premier simulateur d'écoulement et d'un pas de temps du second simulateur d'écoulement. Par exemple, chaque intervalle de temps successifs peut avoir une longueur égale à un pas de temps du premier simulateur d'écoulement.
- On peut déduire les conditions aux limites par une interpolation linéaire des résultats du premier simulateur entre les temps de début et de fin des intervalles de temps successifs. Quant aux indices de productivité numériques, on peut les déduire en comparant des débits calculés par le premier simulateur et des débits calculés par le second simulateur.
- Selon un mode de réalisation, on simule les écoulements de fluides au sein du milieu au moyen du premier simulateur sur un premier maillage discrétisant le milieu poreux en un ensemble de mailles, et on simule les écoulements de fluides aux abords du puits au moyen du second simulateur sur un second maillage discrétisant le puits et ses abords en un ensemble de mailles. Ce second maillage est généré en contraignant des mailles situées sur le bord du second maillage, de façon à ce que leurs interfaces coïncident avec les interfaces des mailles du premier maillage.
- Dans le cas où l'on modélise des écoulements multiphasiques, on met à jour des multiplicateurs d'indices de productivité numériques, au lieu des indices de productivité numériques eux-mêmes, pour chaque phase, en comparant des débits par phase calculés par le premier simulateur et des débits par phase calculés par le second simulateur.
- L'invention concerne également une méthode pour exploiter un réservoir poreux souterrain au moyen d'au moins un puits le traversant, au moins un fluide circulant entre le réservoir et le puits. Selon cette méthode, on acquiert des données relatives à la géométrie du réservoir poreux, à partir desquelles on construit une discrétisation du réservoir en un ensemble de mailles, appelé maillage de réservoir, et on construit une discrétisation du puits et de ses abords en un ensemble de mailles, appelé maillage d'abord de puits. Cette méthode comporte également les étapes suivantes :
- a- on choisit un scénario d'exploitation du réservoir poreux ;
- b- on associe au maillage de réservoir un premier simulateur d'écoulement permettant de simuler l'écoulement des fluides au sein du réservoir, à partir d'au moins les données suivantes : le scénario de production, des données d'entrées relatives au fluide et au réservoir, des indices de productivité numériques permettant de relier des pressions à des débits, des conditions aux limites ;
- c- on associe au maillage d'abord de puits un second simulateur d'écoulement pour
simuler l'écoulement des fluides aux abords du puits, à partir d'au moins les données suivantes : des données d'entrées relatives au fluide et au réservoir, des conditions aux limites ; - d- on modélise les écoulements de fluides au sein du milieu poreux et aux abords du puits, au moyen de la méthode selon l'une des revendications 1 à 7 ; et
- e- on modifie le scénario d'exploitation et l'on répète l'étape d, jusqu'à obtenir un scénario d'exploitation optimal.
- Selon cette méthode d'exploitation, on peut tenir compte d'un endommagement du puits par un fluide de forage en modélisant une invasion du réservoir poreux par le fluide de forage aux étapes d et e.
- Le scénario d'exploitation peut comporter une injection d'une solution de polymère par le puits, et l'on peut alors modéliser les écoulements pour prévenir une venue d'eau. Le scénario d'exploitation peut également comporter une injection d'une solution acide dans le puits, et l'on peut alors modéliser les écoulements pour évaluer l'impact d'une stimulation acide.
- D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après.
-
- La
figure 1 illustre les principales étapes de la méthode selon l'invention. - La
figure 2 montre le schéma de couplage entre le modèle de réservoir et le modèle abord de puits. L'axe T correspond au temps. - La
figure 3 montre le maillage grossier utilisé pour la simulation du champ dans un modèle de réservoir. - La
figure 4 montre le maillage fin utilisé pour simuler les phénomènes détaillés des écoulements autour de puits dans un modèle abord de puits. - La
figure 5 montre les deux maillages utilisés dans le couplage. La figure de gauche représente le maillage de réservoir pour la simulation de champ, et la figure de droite représente le maillage au voisinage de puits dans le modèle abord de puits. Les mailles aux bords (en gris) dans le modèle abord de puits coïncident avec les mailles de même couleur dans le maillage de réservoir. - La
figure 6 montre le maillage grossier pour la simulation du champ dans le cas de l'endommagement par le fluide de forage. Γ x- et Γ x+ correspondent à deux frontières de ce maillage dans la direction x, et Γ y- et Γ y+ correspondent à deux frontières de ce maillage dans la direction y. - La
figure 7 montre le maillage localement raffiné autour du puits pour simuler la solution de référence dans le cas de l'endommagement par fluide de forage. - Les
figures 8A et 8B montrent les maillages de couplage pour simuler l'endommagement par le fluide de forage. Le maillage de gauche (figure 8A ) correspond au maillage pour la simulation de champ, et le maillage de droite (figure 8B ) correspond au maillage dans le modèle abord de puits. - La
figure 9 montre les perméabilités relatives pendant le forage et la production.
L'axe X est la saturation sans unité. L'axe Y est la perméabilité relative. Il n'y a pas d'unité. La courbe "krw drilling" est la courbe de perméabilité relative de l'eau pendant le forage. La courbe "kro drilling" est la courbe de perméabilité relative de l'huile pendant le forage. Les courbes "krw production" et "kro production" sont des courbes de perméabilité relative de l'eau et de l'huile respectivement, pendant la production. - Les
figures 10A et10B comparent le volume d'invasion du fluide de forage simulé par la méthode de couplage avec celui de la solution de référence. Lafigure 10A montre le débit d'invasion pendant le forage. L'axe X est le temps exprimé en jour. L'axe Y est le débit exprimé en m3/jour. La courbe R est la solution de référence. La courbe CM est la simulation avec la méthode de couplage. Lafigure 10B montre le volume d'invasion en fonction du temps. L'axe X est le temps en jour. L'axe Y est le débit en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation avec la méthode de couplage. - La
figure 11 compare les débits de production d'huile. L'axe X est le temps en jour.
L'axe Y est le débit en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation avec la méthode de couplage. La courbe S est la simulation sans l'endommagement du puits. La courbe CK est la simulation avec l'endommagement uniquement par les boues de forage (« cakes »). - La
figure 12 montre les perméabilités dans la couche 3 sur le maillage grossier du modèle de réservoir dans l'application à la prévention de venue d'eau. Il y a un injecteur et un producteur. - La
figure 13 montre le maillage raffiné autour du producteur pour simuler la solution de référence. - Les
figures 14A et 14B montrent les maillages du couplage. Lafigure 14A est le maillage du modèle de réservoir, et lafigure 14B est le maillage du modèle abord de puits. - La
figure 15 montre le débit d'injection de polymère dans le puits traité. L'axe X est le temps en jour. L'axe Y est le débit en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation du modèle de réservoir avec la méthode de couplage. La courbe S est la simulation directe avec le modèle de réservoir sans couplage. La courbe NW est la simulation du modèle abord de puits avec couplage. - Les
figures 16A à 16E montrent le débit d'injection de polymère dans les couches. Lafigure 16A montre le débit d'injection de polymère dans la couche 1. Lafigure 16B montre le débit d'injection de polymère dans la couche 2. Lafigure 16C montre le débit d'injection de polymère dans la couche 3. Lafigure 16D montre le débit d'injection de polymère dans la couche 4. Lafigure 16E montre le débit d'injection de polymère dans la couche 5. L'axe X est le temps en jour. L'axe Y est le débit en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation du modèle de réservoir avec la méthode de couplage. La courbe S est la simulation directe avec le modèle de réservoir sans le couplage. La courbe NW est la simulation du modèle abord de puits avec (évidemment) le couplage. - La
figure 17 montre le débit d'huile du producteur. L'axe X est le temps en jour. L'axe Y est le débit d'huile en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation du modèle de réservoir avec la méthode de couplage. La courbe S est la simulation directe avec le modèle de réservoir sans le couplage. - La
figure 18 montre le débit d'eau du producteur. L'axe X est le temps en jour. L'axe Y est le débit d'eau en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation du modèle de réservoir avec la méthode de couplage. La courbe S est la simulation directe avec le modèle de réservoir sans le couplage. - La
figure 19 montre la courbe de water-cut du producteur. L'axe X est le temps exprimé en jour. L'axe Y est le water-cut, sans unité. La courbe R est la solution de référence. La courbe CM est la simulation du modèle de réservoir avec la méthode de couplage. La courbe S est la simulation directe avec le modèle de réservoir sans le couplage. - Les
figures 20A à 20D montrent une carte de saturation en eau à 1100 jours. Lafigure 20A correspond à la solution de référence dans le champ. Lafigure 20B correspond à la carte obtenue avec le modèle de réservoir en maillage grossier avec couplage. Lafigure 20C montre la solution de référence au voisinage du puits. Lafigure 20D montre la saturation en eau au voisinage de puits simulée avec le modèle abord de puits. - Les
figures 21A à 21D montrent une carte de pression à 1100 jours. Lafigure 21A montre la solution de référence dans le champ. Lafigure 21B montre la solution obtenue avec le modèle de réservoir en maillage grossier avec couplage. Lafigure 21C montre la solution de référence au voisinage du puits. Lafigure 21D montre la solution avec le modèle abord de puits. - L'invention concerne une méthode pour exploiter un milieu poreux souterrain, en injectant un fluide dans le milieu via au moins un puits, et/ou en produisant un fluide présent dans le milieu au moyen d'au moins un puits également. La méthode comporte une modélisation des écoulements de fluides dans le système constitué par le milieu poreux (réservoir et abords des puits) et le puits. Il s'agit donc, en particulier, de modéliser l'injectivité ou la productivité de puits traversant un milieu poreux.
- La
figure 1 illustre les principales étapes de la méthode : - 1. Choix d'un scénario d'exploitation du milieu poreux, scénario de production et/ou scénario d'injection (SCE) ;
- 2. Choix d'un simulateur d'écoulement (RSIM) compatible avec un maillage de réservoir donné, et choix d'un simulateur d'écoulement (NWSIM) compatible avec un maillage d'abord de puits donné ;
- 3. Au moyen d'un couplage entre les deux simulateurs (EST_CAL, et
figure 2 ), estimation des écoulements de fluides, c'est-à-dire, par exemple, du volume injecté ou du volume produit, sur un intervalle de temps donné, et - 4. Détermination du scénario d'exploitation optimal par modification du scénario d'exploitation et répétition de l'étape 3 (OPT).
- Il peut s'agir d'un scénario de production des hydrocarbures contenus dans le milieu poreux (réservoir), ou d'un scénario d'injection de gaz acide, tel que le CO2, dans un réservoir souterrain en vue du stockage du gaz acide. Un scénario se caractérise par la position des puits, la méthode de récupération ou d'injection, les débits et durée d'injection et/ou de production, les conditions opératoires dans ces puits telles que le débit ou la pression de fond.
- Dans le cadre de la production, l'ingénieur réservoir choisit un procédé de production, par exemple le procédé de récupération par injection d'eau, dont il demeure ensuite à préciser le scénario optimal de mise en oeuvre pour le réservoir considéré. La définition d'un scénario optimal consiste, par exemple, à fixer le nombre et l'implantation (position et espacement) des puits injecteurs et producteurs afin de tenir compte au mieux de l'impact d'hétérogénéités au sein du réservoir, par exemple chenaux de perméabilité, fractures, etc., sur la progression des fluides au sein du réservoir. En fonction du scénario choisi, et de la représentation géométrique du réservoir, on est alors capable de simuler la production d'hydrocarbures escomptée, au moyen de l'outil bien connu des spécialistes : le simulateur d'écoulement.
- Le choix d'un scénario, par la définition de multiples caractéristiques techniques, est une étape bien connue des spécialistes.
- Pour choisir un simulateur d'écoulement, il est nécessaire de connaître le type de maillage sur lequel le simulateur doit fonctionner.
- Le « maillage de réservoir » est constitué d'un ensemble de mailles discrétisant spatialement le réservoir (milieu poreux + puits). Un exemple de maillage de réservoir est illustré sur la
figure 3 , ce maillage est grossier. Certaines mailles correspondent à la partie « milieu poreux », d'autres, correspondent à la partie où le puits est foré. On parle pour ces dernières de mailles de puits du maillage de réservoir. - Le « maillage d'abords de puits » est constitué d'un ensemble de mailles discrétisant spatialement le puits et ses abords. Un exemple de maillage d'abord de puits est illustré sur la
figure 4 , ce maillage est fin pour simuler les phénomènes détaillés autour du puits. Ses abords appartiennent donc au milieu poreux dans lequel le puits est foré. Certaines mailles correspondent à la partie « milieu poreux », d'autres, correspondent à la partie « puits ». On parle pour ces dernières de mailles de puits du maillage d'abords de puits. - La génération de ces maillages, que ce soit le maillage de réservoir ou le maillage d'abord de puits, est une étape bien connue du spécialiste, qui connaît de nombreuses méthodes pour les construire. Par exemple, des techniques de construction d'un maillage d'abord de puits sont décrites dans les documents suivants :
- Boe, O., Flynn, J. et Reiso, E., "On Near Wellbore Modeling and Real Time Reservoir Management", SPE 66369, Houston, Texas, USA, 11-14 Feb. 2001.
- On connaît également des méthodes pour construire des maillages de réservoir, à partir de données relatives à la géométrie du milieu (sismiques, diagraphies...), décrites par exemple dans les documents suivants :
- Flandrin, N., Bennis, C. et Borouchaki, H., "3D Hybrid Mesh Generation for Reservoir Simulation", ECMOR, Canne, France, 30 August - 2 September, 2004.
- Pour définir un modèle de réservoir, il est nécessaire d'associer un simulateur d'écoulement au maillage de réservoir. De même, pour définir un modèle d'abord de puits, il est nécessaire d'associer un simulateur d'écoulement au maillage d'abord de puits.
- Comme il est connu de l'homme de l'art, pour fonctionner, un simulateur d'écoulement a besoin de certaines données, dites données d'entrée :
- caractéristiques géométriques du réservoir, caractéristiques de la roche, caractéristiques des fluides en place et des fluides injectés (masse volumique, viscosité), courbes de perméabilités relatives, courbes de pression capillaire, saturations initiales en fluides,...
- conditions aux limites du domaine simulé et aux puits où sont injectés ou produits des fluides. Les conditions aux limites, sont les valeurs de données dynamiques, telles que la pression, le débit ou le flux, les saturations en fluides, sur les bords du maillage ou dans les mailles qui forment les bords du maillage de réservoir ou d'abord de puits. Un exemple de conditions aux limites peut être : un flux nul sur tous les bords du maillage, ou des saturations et pressions imposées sur les mailles aux bords du maillage.
- éventuellement des Indices de Productivité numériques (IP). La connexion entre la pression dans les mailles traversées par un puits et les pressions dans le puits lui-même est réalisée à l'aide d'un Indice de Productivité numérique (IP). L'IP numérique peut être calculé par une formule analytique dans le code ou donné par l'utilisateur du logiciel (le simulateur). En général, le simulateur calcule un IP numérique par une formule analytique au début de simulation. Mais, si l'utilisateur donne un IP numérique dans le jeu de donnée d'entrée, c'est l'IP numérique de l'utilisateur qui est prise en compte dans la simulation.
- Selon l'invention, il est possible d'utiliser tout type de simulateur d'écoulement, que ce soit pour le modèle de réservoir, ou le modèle d'abord de puits. En effet, un objet de l'invention concerne une méthode de couplage, qui permet de coupler de manière très simple un modèle de réservoir, pour la simulation du réservoir, et un modèle d'abord de puits, qui est un modèle autonome pour simuler les phénomènes détaillés autour du puits.
- Concernant le simulateur du modèle de réservoir, il peut s'agir du logiciel PumaFlow® (IFP, France) par exemple.
- Concernant le simulateur du modèle d'abord de puits, on peut utiliser celui décrit dans le document suivant : DING, Y., RENARD, G. : "Evaluation of Horizontal Well Performance after Drilling Induced Formation Damage", J. of Energy Resources Technology, Vol. 127, Sept., 2005.
- Il s'agit ici de modéliser l'injectivité ou la productivité d'un puits traversant le milieu poreux, et permettant l'exploitation de ce milieu. Cette modélisation se fait sur un intervalle de temps donné, D=[T0 ; Tn ]. Par exemple, on modélise le comportement du système milieu+puits sur 20 ans, en considérant le scénario d'exploitation préalablement choisi.
- La technique utilisée ici, consiste à réaliser un couplage entre les deux simulateurs d'écoulement.
- Un maillage grossier est souvent utilisé pour le modèle de réservoir, et un maillage fin est habituellement nécessaire pour simuler les phénomènes détaillés autour du puits. La
figure 5 montre les deux maillages utilisés dans le couplage. La figure de gauche représente le maillage de réservoir pour la simulation de champ, et la figure de droite représente le maillage au voisinage de puits dans le modèle abord de puits. Les mailles aux bords (en gris) dans le modèle abord de puits coïncident avec les mailles de même couleur dans le maillage de réservoir. La croix indique l'emplacement du puits. - Les pas de temps utilisés dans le modèle abord de puits sont généralement beaucoup plus petits que ceux du modèle de réservoir. Le modèle de réservoir est principalement utilisé pour simuler les écoulements dans le réservoir pris dans son intégralité.
- Le temps T0 est le temps auquel commence le couplage. L'algorithme de couplage dans un cadre général comporte les étapes suivantes, illustrées sur la
figure 2 : - 3a. On initialise les modèles.
- on initialise le modèle de réservoir (RINIT) en affectant aux mailles du maillage de réservoir des valeurs de porosité, perméabilité, pression et saturations en fluide. L'initialisation comporte également la définition de conditions aux limites du modèle de réservoir. Ces conditions peuvent être définies par un flux nul (pas d'échange vers l'extérieur du domaine) ou par un flux ou une pression imposé sur les bords extérieurs des mailles de bord du maillage du modèle de réservoir (échange avec l'extérieur). Les conditions opératoires dans ces puits, telles que le débit ou la pression de fond, sont imposées sous forme d'un historique d'injection pour les injecteurs, et d'un historique de production dans les producteurs.
- on initialise le modèle d'abord de puits (NWINIT) en affectant aux mailles du maillage d'abord de puits des valeurs de porosité, perméabilité, pression et saturations en fluide. Cette affectation est réalisée en utilisant des techniques de mise à l'échelle des résultats du modèle de réservoir. Ces techniques étant connues des spécialistes. L'initialisation comporte également la définition de conditions aux limites du modèle d'abord de puits. Ces conditions peuvent également être définies au moyen des résultats du modèle de réservoir.
- 3b. On définit au moins un pas de temps, noté ΔT, pour échanger des données dynamiques entre le modèle de réservoir et le modèle d'abord de puits, pendant la modélisation sur l'intervalle de temps D.
Ce pas de temps ΔT peut être choisi en fonction du pas de temps ΔTR du simulateur d'écoulement du modèle de réservoir, et le pas de temps ΔTNW du simulateur d'écoulement du modèle d'abord de puits (ΔTR >ΔTNW).
Théoriquement, ΔT doit être le plus petit possible pour assurer la convergence des solutions dans les deux modèles. Mais généralement l'utilisation du pas de temps utilisé pour la simulation du modèle de réservoir est suffisante. Mais, du point de vue pratique, on a parfois besoin de réaliser une simulation d'abord de puits de façon autonome plus longtemps. Ceci se traduit par une réduction de la fréquence de couplage. C'est pourquoi, selon la méthode, le pas de temps ΔT pour les échanges de données entre le modèle de réservoir et le modèle d'abord de puits est un paramètre ajustable.
Selon un mode de réalisation, le pas de temps ΔT peut varier sur l'intervalle de temps D. On peut par exemple utiliser un premier pas de temps entre T0 et Ti, et un second pas de temps entre Ti et Tn. Un exemple d'une telle application est illustré ci-après. Sur lafigure 2 , on note RSIM(T1 ) une simulation du simulateur de réservoir effectuée entre T0 et T1 , et NWSIM(T1 ) une simulation du simulateur d'abord de puits effectuée entre T0 et T1 . - 3c. On réalise une simulation d'écoulements avec le modèle de réservoir entre le temps T0 et le temps T1 = T0 + ΔT.
Les résultats de cette simulation sont :- la pression et les saturations en fluides à la fin du pas de temps dans chaque maille du maillage de réservoir, en particulier dans les mailles qui sont communes avec les mailles du maillage d'abord de puits et qui vont servir de conditions aux limites du modèle d'abord de puits ;
- les débits en fluide (eau, huile, gaz) et les pressions dans les puits d'injection et de production.
- 3d. On met à jour les conditions aux limites (MAJCL) du modèle d'abord de puits en utilisant les résultats de la simulation d'écoulements réalisée avec le modèle de réservoir entre T0 et T1 (étape 3c).
Les conditions aux limites, sont les valeurs de données dynamiques, telles que la pression ou le flux, les saturations, dans les mailles qui forment les limites du maillage de réservoir ou d'abord de puits. Selon un exemple, les conditions aux limites sont définies par un flux nul sur tous les bords du maillage des abords de puits, et une porosité très grande (1000000, par exemple) dans toutes ces mailles.
Ainsi, au cours de cette étape, on utilise les résultats du simulateur d'écoulement du modèle de réservoir, pour déterminer des valeurs que l'on impose comme conditions aux limites pour le simulateur d'écoulement du modèle d'abord de puits au temps T0.
Les conditions aux limites peuvent être calculées à chaque pas de temps du modèle d'abord de puits en interpolant linéairement les résultats de simulation du modèle de réservoir entre T0 et T1. - 3e. On réalise une simulation d'écoulements au voisinage du puits avec le modèle d'abords de puits entre le temps T0 et le temps T1, avec les conditions aux limites mises à jour à l'étape 3d.
Les résultats de cette simulation sont, au moins :- la pression et les saturations en fluides à la fin du pas de temps dans chaque maille du modèle d'abords de puits;
- les débits en fluide (eau, huile, gaz) et les pressions dans le puits d'injection ou de production selon du type de puits modélisé dans le modèle d'abords de puits.
Ces résultats permettent de déterminer un Indice de Productivité numérique (IP).
- 3f. La connexion entre la pression dans les mailles traversées par un puits et les pressions dans le puits lui-même est réalisée à l'aide d'un Indice de Productivité numérique (IP). Les formules de Peaceman sont en général utilisées pour calculer cet indice. On met alors à jours les indices de productivité numériques (MAJIP) du modèle de réservoir, en utilisant les résultats de la simulation d'écoulements réalisée avec le modèle d'abord de puits entre T0 et T1. En fait, si à la fin de la simulation, au temps T1 , les résultats simulés au puits avec le modèle d'abord de puits et avec le modèle de réservoir ne sont pas les mêmes, les indices de productivité numériques dans le modèle de réservoir sont modifiés pour ajuster les résultats de la simulation du modèle de réservoir à ceux du modèle d'abord de puits.
- 3g. On répète les étapes 3c (éventuellement 3b) à 3f, avec un nouvel intervalle de temps (de T1 à T2 , puis de T2 à T3, ..., puis de Tn-1 à Tn )
-
- i :
- numéro de maille de puits dans le maillage (de réservoir ou d'abord de puits)
- p :
- phase du fluide. Les phases p peuvent être l'eau, l'huile ou le gaz
- Q p,i :
- débit de la phase p dans la maille de puits i du maillage (de réservoir ou d'abord de puits)
- λ p,i :
- mobilité de la phase p dans la maille de puits i du maillage (de réservoir ou d'abord de puits); λ p,i dépend essentiellement de la perméabilité relative et de la viscosité de la phase p
- IP¡ :
- indice de productivité numériques dans la maille de puits i du maillage (de réservoir ou d'abord de puits)
- Pp,i :
- pression de la phase p dans la maille de puits i du maillage (de réservoir ou d'abord de puits)
- P wf,i :
- pression dans le puits, au fond, au niveau du réservoir dans la maille de puits i du maillage (de réservoir ou d'abord de puits)
- L'indice de productivité numérique IP prend en compte : l'effet géométrique de la maille de puits i du maillage, la perméabilité du milieu poreux dans la maille du puits et un coefficient de skin. Un coefficient de skin est un coefficient, bien connu de l'homme de l'art, utilisé pour représenter l'endommagement d'un puits dans une maille.
-
- i :
- numéro de maille de puits dans le maillage de réservoir
- j :
- numéro de maille de puits dans le maillage d'abord de puits
- Wi :
- ensemble des mailles de puits du maillage d'abord de puits correspondant à un raffinement de la maille de puits i du maillage de réservoir
- p :
- phase du fluide. Les phases p peuvent être l'eau (w), l'huile (o) ou le gaz (g)
- IP r,i :
- indice de productivité numériques dans la maille de puits i du maillage de réservoir, et utilisé dans le modèle de réservoir
- P nw,p,j :
- pression de la phase p dans la maille de puits j du maillage d'abord de puits, calculée avec le modèle abord de puits
- P r,p,i :
- pression de la phase p dans la maille de puits i du maillage de réservoir, calculée avec le modèle de réservoir
- P wf,j :
- pression dans le puits au niveau du réservoir dans la maille de puits j du maillage d'abord de puits
- IP nw,j :
- indice de productivité numériques dans la maille de puits j du maillage d'abord de puits, et utilisé dans le modèle d'abord de puits
- Les variables IPi , P nw,p,j , P r,p,i et P wf,j sont fonction du temps T.
-
- Q nw,i :
- débit du fluide (unique phase) calculé avec le modèle d'abord de puits dans la section correspondant à la partie du puits dans la maille de puits i du maillage de réservoir
- Q r,i :
- débit du fluide (unique phase) calculé avec le modèle de réservoir dans la même section, correspondant à la partie du puits dans la maille de puits i du maillage de réservoir
- Cette formule montre clairement que la mise à jour de l'indice de productivité numérique correspond à la correction du débit de fluide du modèle de réservoir par rapport au débit de fluide du modèle abord de puits : si les deux modèles donnent le même résultat en terme de débit, alors
- En testant divers scénarios, caractérisés par exemple par diverses implantations respectives des puits injecteurs et producteurs, et en simulant la production d'hydrocarbures pour chacun d'eux selon l'étape 3, on peut sélectionner le scénario optimal. Le scénario optimal est le scénario permettant d'obtenir une production optimale du gisement dans la cadre de la production d'un réservoir, ou le scénario permettant d'obtenir l'injectivité optimale dans le gisement dans le cadre d'injection de fluide dans le réservoir (injection d'eau pour une production améliorée, ou injection de gaz acides).
- Pour tester divers scénarios d'exploitation, on modifie (ΔSCE) le scénario sélectionné à l'étape 1, par exemple en modifiant l'emplacement d'un puits.
- On optimise alors l'exploitation du gisement, en mettant en oeuvre, sur le champ, le scénario de production ainsi sélectionné.
- Selon l'invention, il est tout fait possible de coupler un modèle de réservoir avec plusieurs modèles d'abord de puits.
- Selon un mode de réalisation particulier de l'invention, l'étape 2, au cours de laquelle on construit les maillages, est modifiée.
- La simulation, réalisée au moyen du modèle de réservoir dans l'étape 3c, fournit des propriétés dynamiques des fluides comme la pression ou les saturations dans la période allant de T0 à T1 sur toutes les mailles grossières. Cependant, la détermination des conditions aux limites dans l'étape 3b nécessite l'interpolation de la pression ou du flux sur les bords du modèle abord de puits. Pour réduire les erreurs dans l'interpolation, nous pouvons contraindre, lors de la génération de maillage, les mailles de bord du modèle abord de puits pour qu'elles coïncident avec les interfaces des mailles du modèle de réservoir. De plus, les mailles de bord dans le modèle abord de puits sont aussi contraintes afin qu'elles coïncident avec des mailles du modèle de réservoir (
figure 3 ). De cette façon, le transfert des données dynamiques du modèle de réservoir vers le modèle abord de puits est direct sur ces mailles. Dans le modèle abord de puits lui-même, les conditions aux limites sont de flux nul. Afin de maintenir les propriétés dynamiques aux bords du modèle, des porosités de très grande valeur (1000000, par exemple) sont affectées aux mailles de bord. Ce type de conditions aux limites est consistant avec la plupart des modèles d'écoulement, et son implémentation est facile. - Pour certains problèmes, les changements d'écoulement autour du puits sont liés aux écoulements multiphasiques. Dans ce cas-là, nous pouvons aussi faire la mise à jour des indices de productivité numérique par phase. Pour cela, on reformule la relation de pression/débit, en introduisant un coefficient, appelé multiplicateur de l'indice de productivité :
- M p,i est le multiplicateur de l'indice de productivité pour la phase p dans la maille de puits i.
-
- Q r,p,i (T 1) :
- débit de la phase p calculé par le modèle de réservoir dans la maille de puits i du maillage de réservoir au temps T1
- Q nw,p,i (T 1) :
- débit de la phase p calculé par le modèle d'abord de puits dans le même secteur du puits (voir l'ensemble Wi ) au temps T1
- M p,i (T 0) :
- multiplicateur d'indice de productivité numérique pour la phase p dans le modèle de réservoir aux temps T0 (avant mise à jour du modèle)
- M p,i (T 1) :
- multiplicateur d'indice de productivité numérique pour la phase p dans le modèle de réservoir aux temps T1 (après mise à jour du modèle)
- La méthode de couplage selon l'invention peut être utilisée pour modéliser différents phénomènes détaillés autour du puits, comme par exemple, l'endommagement par fluide de forage ou de complétion, la stimulation acide, l'écoulement non-Darcéen autour du puits, le problème de gaz à condensat, le dépôt d'asphaltène, l'endommagement par l'injection de CO2, la prévention de venue d'eau ou de gaz, la venue de sable, les dépôts minéraux, l'impact des complétions, etc. Ici, nous présentons en particulier un exemple d'application pour l'endommagement de la formation pétrolifère par le fluide de forage lors du forage du puits, et un exemple d'application pour la prévention de venue d'eau lorsqu'un puits en production produit une quantité importante d'eau, et que l'on cherche à réduire cette production d'eau.
- Pour simplifier encore la méthode de couplage, nous faisons la mise à jour des données en utilisant les valeurs au temps Tn , au lieu de l'interpolation linéaire en temps entre Tn et T n+1 , pour la simulation du modèle abord de puits dans la période allant de Tn à T n+1 . Ce choix est intéressant, car il permet de faire des simulations en parallèle sur différentes machines pour le modèle de réservoir et le modèle abord de puits.
- Un modèle de réservoir standard est utilisé pour la simulation du champ. Le modèle abord de puits développé par DING, Y. and RENARD, G. : "Evaluation of Horizontal Well Performance after Drilling Induced Formation Damage" J. of Energy Resources Technology, Vol. 127, Sept., 2005, est utilisé pour simuler l'endommagement de la formation par le forage. Il faut noter que la physique avancée de l'endommagement n'est pas modélisée dans la simulation du champ avec le modèle de réservoir.
- Un réservoir de taille 1000m x 1000m x 10m est considéré. Un maillage cartésien avec 20 mailles dans la direction x, 20 mailles dans la direction y et 1 maille dans la direction z est utilisé pour la simulation du champ (
figure 6 ). Les tailles de mailles sont donc 50m x 50m x 10m. La pression initiale de réservoir est de 200 bars. Un puits producteur doit être foré dans le block (15, 15, 1). Il est représenté par un cercle noir sur lafigure 6 . L'endommagement de ce puits par le fluide de forage est étudié avec la méthode selon l'invention. - Le réservoir est homogène de perméabilité 200 mD et de porosité 0.15. Les conditions aux limites de ce réservoir sont des flux nuls, sauf sur le bord Γ x- (
figure 6 ), où la pression est constante (200 bars). - Pour obtenir la solution de référence, le maillage est raffiné autour du puits (
figure 7 ). Un modèle spécifique, qui prend en compte la physique avancée de l'endommagement, est utilisé sur ce maillage pour simuler la solution de référence. Comme l'endommagement par le fluide de forage est généralement limité à quelques centimètres ou quelques dizaines de centimètres autour de puits, nous avons besoin de très petites mailles dans la zone raffinée (Tableau 1). Le diamètre du puits est de 21.6 cm. Pour que le puits soit inclus dans une maille, la taille de la maille puits est de 22 cm. Les autres mailles autour du puits sont beaucoup plus petites avec une taille de 2 cm. Les maillages utilisés pour le couplage sont illustrés sur lesfigures 8A et 8B . Le maillage du modèle abord de puits (figure 8B ) correspond à la zone raffinée et aux mailles autour dans le maillage de référence. Les mailles aux bords du modèle abord de puits coïncident avec des mailles du modèle de réservoir.Tableau 1-Tailles des mailles autour du puits Taille de maille dans la direction x (m) Taille de maille dans la direction y (m) 50 42.7 30 20 16 8 4 2 1 0.51 0.3 50 42.7 30 20 16 8 4 2 1 0.51 0.3 0.16 0.08 0.04 0.02 0.02 0.02 0.02 0.16 0.08 0.04 0.02 0.02 0.02 0.02 0.02 0.02 0.22 0.22 0.02 0.02 0.02 0.02 0.02 0.04 0.08 0.02 0.02 0.02 0.02 0.02 0.04 0.08 0.16 0.16 0.3 0.51 1 2 4 8 16 20 30 42.7 50 0.3 0.51 1 2 4 8 16 20 30 42.7 50 - Il est supposé que le réservoir est épais, et que ce modèle correspond seulement à la première couche du réservoir. Le temps de contact entre le fluide de forage et le réservoir est de 2 jours. La pression au cours du forage au fond du puits est de 250 bars. La perméabilité et l'épaisseur du « cake » externe formé par la boue de forage sont égales à 0.001 mD et 0.2 cm. L'épaisseur du « cake » interne est de 2 cm avec une perméabilité moyenne réduite à 20 mD durant la période de forage et de 40 mD dans la période de production. La viscosité du fluide de forage est de 30 cPo. L'hystérésis de la perméabilité relative entre les périodes de forage et de mise en production est présenté à la
figure 9 . Une saturation en eau irréductible de 30% liée au filtrat (fluide de forage) qui va envahir la formation pendant la phase de forage va rester bloquée dans le milieu poreux quand le puits va être remis en production. - Les volumes d'invasion de fluide de forage sont comparés à la
figure 10 pour la simulation avec la méthode de couplage et la solution de référence obtenue en utilisant le maillage avec le raffinement local (figure 7 ). Les pas de temps pour la mise à jour des données dans le couplage sont présentés dans le Tableau 2. Lafigure 10 montre que le volume d'invasion de fluide est correctement simulé avec la méthode de couplage. Le petit écart entre la solution de couplage et la solution de référence dans la période entre 0.1 et 0.3 jour peut être amélioré en utilisant des petits pas d'itération en temps pour échanger les données dans le couplage.Tableau 2-Pas de temps pour la mise à jour des données dans le couplage Période (jour) Pas de temps (jour) 0 - 0.01 0.001 0.01 - 0.1 0.01 0.1 - 3 0.1 3 -10 1 10 - 200 10 - Après les 2 jours de forage, le puits est fermé pendant 1 jour pour la mise en place de sa complétion, puis il est mis en production. Le couplage est effectué jusqu'au 10ème jour. Au-delà de 10 jours, l'effet de l'endommagement autour de puits devient stable et les IPs numériques dans le modèle de réservoir ne changent quasiment plus. Nous n'avons plus besoin de couplage pour continuer la simulation du champ avec le modèle de réservoir. La courbe de production d'huile simulée par le modèle de réservoir, qui est couplé avec le modèle abord de puits pendant les 10 premiers jours, est présentée à la
figure 11 . Cette courbe est très proche de la solution de référence. - Si l'endommagement n'est pas pris en compte ou si seulement la présence des « cakes » est considérée dans la simulation, les résultats sont très imprécis avec des erreurs de plus de 20% (
figure 11 ). Prendre en compte des phénomènes autour de puits tel que l'endommagement par le fluide de forage est important pour la gestion de réservoir, et la méthode de couplage proposée convient parfaitement pour simuler ce type de problèmes. - Dans la procédure de prévention de venue d'eau, une solution de polymère est injectée dans un puits producteur pour une courte durée dans le but de réduire la grande quantité d'eau produite en même temps que l'huile. Une partie du polymère est absorbée sur la roche, et une autre partie est dispersée dans l'eau. Le polymère injecté a pour effet de réduire la mobilité de la phase eau par l'augmentation de sa viscosité et par la diminution de la perméabilité relative de cette phase. Donc, dans la méthode de couplage, l'approche la plus adaptée est de mettre à jour le multiplicateur d'IP numérique pour la phase eau.
- Un réservoir de 1000m x 1000m x 25m est considéré comme exemple. Un maillage cartésien avec 20 mailles dans la direction x, 20 mailles dans la direction y et 5 mailles dans la direction z est utilisé pour la simulation du champ. La taille des mailles est de 50m x 50m x 5m. Le réservoir est hétérogène. La perméabilité est présentée à la
figure 12 . Le rapport des perméabilités dans les directions verticale et horizontale est de 0.1. La pression initiale du réservoir est de 200 bars. - Il y a un puits injecteur (INJ) et un puits producteur (PROD) comme montré à la
figure 12 . La pression au puits injecteur est imposée à 300 bars, et la pression au puits producteur est contrainte à 150 bars pendant la production. Après production pendant 1000 jours, le water-cut (débit d'eau rapporté au débit total) du puits producteur atteint 85%. La procédure de prévention de venue d'eau est alors appliquée pour réduire la quantité d'eau produite. Une solution de polymère avec une concentration de 2500 ppm est injectée dans le producteur avec une pression de fond de 300 bars pendant 2 jours. Ensuite, le puits est remis en production. Cette procédure de prévention de venue d'eau est simulée avec la méthode selon l'invention. - Afin d'avoir une solution de référence, un raffinement local autour du puits producteur est utilisé (
figure 13 ). La taille des mailles autour du puits est de 0.617 m dans la direction x. Le maillage pour le couplage est présenté à lafigure 14 . Des mailles aux bords du modèle abord de puits coïncident avec des mailles dans le modèle de réservoir. La physique du polymère peut être considérée dans les deux modèles (modèle abord de puits et modèle de réservoir).Tableau 3 - Pas de temps pendant le couplage Période (jour) Pas de temps (jour) 0 - 950 - 950-970 2 970 - 1000 28 1000 - 1000.1 0.01 1000.1 - 1005 0.1 1005 - 1030 1 1030 - 1100 2 1100-3000 - - Le couplage commence à 950 jours et se termine à 1100 jours, soit une période de 150 jours au total. Les pas de temps pour les échanges des données dans la méthode de couplage sont présentés dans le Tableau 3. Dans les premiers 50 jours (de 950 à 1000 jours) de couplage, il n'y a pas de polymère injecté. Cette période est uniquement utilisée pour assurer une bonne initialisation du modèle abord de puits. Les IPs numériques globaux sont mis à jour au début du couplage (de 950 à 970 jours) pour prendre en compte les effets des maillages entre le modèle de réservoir et le modèle abord de puits. Pendant la période d'injection de polymère (entre 1000 et 1002 jours), les IPs numériques globaux sont encore recalculés pour intégrer l'effet induit par le polymère injecté (on pourrait aussi mettre à jour les multiplicateurs d'IP numérique pour la phase eau). Mais, quand le puits est remis en production (à 1003 jours), les multiplicateurs d'IP numérique pour la phase d'eau sont mis à jour.
- La
figure 15 compare les débits d'injection de polymère dans le puits pour les différentes simulations: la solution de référence, la simulation sur le maillage de réservoir avec couplage, la simulation directe sur le maillage de réservoir sans couplage et la simulation avec le modèle abord de puits (avec couplage). Lesfigures 16A à 16E montrent les mêmes comparaisons couche par couche. Pour la simulation directe avec le maillage de réservoir sans couplage, le volume de polymère injecté est largement surestimé. Quand la simulation du maillage grossier est couplée avec le modèle abord de puits, les résultats sont significativement améliorés. Au début du couplage, le débit d'injection est grand, mais il est rapidement corrigé par la mise à jour de l'IP dû au couplage. Si l'on veut avoir plus de précision sur le débit d'injection de polymère, il suffit de se référer aux résultats de simulation avec le modèle abord de puits. Avec ce modèle, le volume injecté et la distribution du polymère autour du puits sont tous les deux correctement simulés. - Les
figures 17, 18 et19 présentent les courbes de débits d'huile, d'eau et de water-cut pour le modèle de réservoir avec couplage, le modèle de réservoir sans couplage et la solution de référence. Les résultats du modèle de réservoir avec couplage sont globalement satisfaisants. Lafigure 20 présente la carte de saturation d'eau à la fin du couplage (1100 jours), et lafigure 21 montre la carte de pression à 1100 jours. Comparé aux solutions de référence, le couplage donne des résultats globalement satisfaisants.
Claims (11)
- Méthode mise en oeuvre par ordinateur pour modéliser des écoulements de fluides au sein d'un milieu poreux traversé par au moins un puits, dans laquelle on utilise un premier simulateur d'écoulement permettant de simuler l'écoulement des fluides au sein du milieu poreux à partir d'indices de productivité numériques reliant des pressions de fluides à des débits de fluides, et on utilise un second simulateur d'écoulement pour simuler l'écoulement des fluides aux abords du puits à partir de conditions aux limites, caractérisé en ce que:a- on simule les écoulements de fluides au sein du milieu au moyen du premier simulateur sur un intervalle de temps défini entre des temps T0 et T1 , et on en déduit des conditions aux limites mises à jour pour le second simulateur ;b- on simule les écoulements de fluides aux abords du puits au moyen du second simulateur sur le même intervalle de temps, en utilisant les conditions aux limites mises à jour, et on en déduit des indices de productivité numériques mis à jour pour le premier simulateur ; etc- on modélise les écoulements de fluides au sein du milieu poreux pendant une période de temps entre T0 et Tn , où Tn >T1 , en réitérant les étapes a et b, pour des intervalles de temps successifs compris entre T0 et Tn .
- Méthode selon la revendication 1, dans laquelle chaque intervalle de temps successif a une longueur qui est fonction d'un pas de temps de calcul du premier simulateur d'écoulement et d'un pas de temps du second simulateur d'écoulement.
- Méthode selon la revendication 1, dans laquelle chaque intervalle de temps successifs a une longueur égale à un pas de temps du premier simulateur d'écoulement.
- Méthode selon l'une des revendications précédentes, dans laquelle on déduit les conditions aux limites par une interpolation linéaire des résultats du premier simulateur entre les temps de début et de fin desdits intervalles de temps successifs.
- Méthode selon l'une des revendications précédentes, dans laquelle on déduit les indices de productivité numériques en comparant des débits calculés par le premier simulateur et des débits calculés par le second simulateur.
- Méthode selon l'une des revendications précédentes, dans laquelle on simule les écoulements de fluides au sein du milieu au moyen du premier simulateur sur un premier maillage discrétisant le milieu poreux en un ensemble de mailles, et on simule les écoulements de fluides aux abords du puits au moyen du second simulateur sur un second maillage discrétisant le puits et ses abords en un ensemble de mailles, ledit second maillage étant généré en contraignant des mailles situées sur le bord dudit second maillage, de façon à ce que leurs interfaces coïncident avec les interfaces des mailles dudit premier maillage.
- Méthode selon l'une des revendications précédentes, dans laquelle on modélise des écoulements multiphasiques, et on met à jour des multiplicateurs d'indices de productivité numériques, au lieu des indices de productivité numériques eux-mêmes, pour chaque phase, en comparant des débits par phase calculés par le premier simulateur et des débits par phase calculés par le second simulateur.
- Méthode, selon l'une des revendications 1 à 7, pour exploiter un réservoir poreux souterrain au moyen d'au moins un puits traversant la réservoir, au moins un fluide circulant entre le réservoir et le puits, dans laquelle on acquiert des données relatives à la géométrie dudit réservoir poreux, à partir desquelles on construit une discrétisation du réservoir en un ensemble de mailles, appelé maillage de réservoir, et on construit une discrétisation du puits et de ses abords en un ensemble de mailles, appelé maillage d'abord de puits, caractérisé en ce que l'on réalise les étapes suivantes :a- on choisit un scénario d'exploitation du réservoir poreux ;b- on associe au maillage de réservoir un premier simulateur d'écoulement permettant de simuler l'écoulement des fluides au sein du réservoir, à partir d'au moins les données suivantes : le scénario de production, des données d'entrées relatives au fluide et au réservoir, des indices de productivité numériques permettant de relier des pressions à des débits, des conditions aux limites ;c- on associe au maillage d'abord de puits un second simulateur d'écoulement pour simuler l'écoulement des fluides aux abords du puits, à partir d'au moins les données suivantes : des données d'entrées relatives au fluide et au réservoir, des conditions aux limites ;d- on modélise les écoulements de fluides au sein du réservoir et aux abords du puits, au moyen de la méthode selon l'une des revendications 1 à 7 ; ete- on modifie le scénario d'exploitation et l'on répète l'étape d, jusqu'à obtenir un scénario d'exploitation optimal.
- Méthode selon la revendication 8, dans laquelle on tient compte d'un endommagement du puits par un fluide de forage en modélisant une invasion du réservoir poreux par ledit fluide de forage aux étapes d et e.
- Méthode selon la revendication 8, dans laquelle le scénario d'exploitation comporte une injection d'une solution de polymère par le puits, et l'on modélise les écoulements pour prévenir une venue d'eau.
- Méthode selon la revendication 8, dans laquelle le scénario d'exploitation comporte une injection d'une solution acide dans le puits, et l'on modélise les écoulements pour évaluer l'impact d'une stimulation acide.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0902533A FR2945879B1 (fr) | 2009-05-20 | 2009-05-20 | Methode d'exploitation de milieu poreux au moyen d'une modelisation d'ecoulements de fluide |
Publications (2)
Publication Number | Publication Date |
---|---|
EP2253797A1 true EP2253797A1 (fr) | 2010-11-24 |
EP2253797B1 EP2253797B1 (fr) | 2020-02-19 |
Family
ID=41426264
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP10290174.1A Active EP2253797B1 (fr) | 2009-05-20 | 2010-04-01 | Méthode d'exploitation de milieu poreux au moyen d'une modélisation d'écoulements de fluide |
Country Status (4)
Country | Link |
---|---|
US (1) | US8694297B2 (fr) |
EP (1) | EP2253797B1 (fr) |
CA (1) | CA2704060C (fr) |
FR (1) | FR2945879B1 (fr) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2989200A1 (fr) * | 2012-04-10 | 2013-10-11 | IFP Energies Nouvelles | Procede de selection des positions de puits a forer pour l'exploitation d'un gisement petrolier |
FR2997721A1 (fr) * | 2012-11-08 | 2014-05-09 | Storengy | Radonip : nouvelle methodologie de determination des courbes de productivite des puits d'exploitation de stockages et gisements de fluides compressibles |
FR3002270A1 (fr) * | 2013-02-21 | 2014-08-22 | IFP Energies Nouvelles | Procede d'exploitation d'un reservoir geologique au moyen d'un modele de reservoir cale et coherent vis a vis des proprietes d'ecoulement |
WO2021118714A1 (fr) * | 2019-12-11 | 2021-06-17 | Exxonmobil Upstream Research Company | Méthodologie de semi-élimination permettant de simuler des caractéristiques d'écoulement élevées dans un réservoir |
Families Citing this family (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2702965C (fr) | 2007-12-13 | 2014-04-01 | Exxonmobil Upstream Research Company | Partitionnement parallele adaptatif de donnees sur une simulation de reservoir utilisant une grille non structuree |
CN101878434B (zh) | 2008-04-10 | 2014-05-07 | 普拉德研究及开发股份有限公司 | 用于表征被井眼穿过的地质构造的方法 |
US8078328B2 (en) * | 2008-05-03 | 2011-12-13 | Saudi Arabian Oil Company | System, program product, and related methods for performing automated real-time reservoir pressure estimation enabling optimized injection and production strategies |
US9134457B2 (en) * | 2009-04-08 | 2015-09-15 | Schlumberger Technology Corporation | Multiscale digital rock modeling for reservoir simulation |
EP2534605B1 (fr) | 2010-02-12 | 2020-06-17 | Exxonmobil Upstream Research Company | Procédé et système de partitionnement de modèles de simulation parallèle |
US9367564B2 (en) | 2010-03-12 | 2016-06-14 | Exxonmobil Upstream Research Company | Dynamic grouping of domain objects via smart groups |
US8583410B2 (en) * | 2010-05-28 | 2013-11-12 | Ingrain, Inc. | Method for obtaining consistent and integrated physical properties of porous media |
US9540911B2 (en) * | 2010-06-24 | 2017-01-10 | Schlumberger Technology Corporation | Control of multiple tubing string well systems |
WO2012102784A1 (fr) | 2011-01-26 | 2012-08-02 | Exxonmobil Upstream Research Company | Procédé d'analyse des compartiments d'un réservoir en utilisant la structure topologique d'un modèle de terre 3d |
WO2013043531A1 (fr) * | 2011-09-20 | 2013-03-28 | Bp Corporation North America Inc. | Génération automatique d'un affinement local de maillage au niveau de fractures hydrauliques pour simuler des réservoirs compacts de gaz |
GB2512372B (en) * | 2013-03-28 | 2020-07-29 | Total Sa | Method of modelling a subsurface volume |
EP3008281A2 (fr) | 2013-06-10 | 2016-04-20 | Exxonmobil Upstream Research Company | Planification interactive d'un site de puits |
US10480314B2 (en) | 2013-07-26 | 2019-11-19 | Schlumberger Technology Corporation | Well treatment |
US9864098B2 (en) | 2013-09-30 | 2018-01-09 | Exxonmobil Upstream Research Company | Method and system of interactive drill center and well planning evaluation and optimization |
US11414975B2 (en) * | 2014-07-14 | 2022-08-16 | Saudi Arabian Oil Company | Quantifying well productivity and near wellbore flow conditions in gas reservoirs |
US9816366B2 (en) * | 2014-07-14 | 2017-11-14 | Saudi Arabian Oil Company | Methods, systems, and computer medium having computer programs stored thereon to optimize reservoir management decisions |
CN104500040B (zh) * | 2014-10-16 | 2017-06-06 | 西南石油大学 | 水平井酸化过程中井筒多段流体移动界面跟踪方法 |
CN105626007B (zh) * | 2014-11-07 | 2018-06-15 | 中国石油化工股份有限公司 | 基于岩心尺度油藏中不同部位过水倍数计算方法 |
US10280722B2 (en) | 2015-06-02 | 2019-05-07 | Baker Hughes, A Ge Company, Llc | System and method for real-time monitoring and estimation of intelligent well system production performance |
CA2992714A1 (fr) * | 2015-08-21 | 2017-03-02 | Halliburton Energy Services, Inc. | Procede et protocole de modelisation precise de formations de champ proche dans des simulations de puits de forage |
CN111094694B (zh) * | 2017-09-08 | 2022-08-12 | 罗克斯软件解决方案有限公司 | 井裂缝建模 |
CN111712823A (zh) * | 2017-12-14 | 2020-09-25 | 斯伦贝谢技术有限公司 | 用于模拟储层模型的系统和方法 |
CN109882149B (zh) * | 2018-01-29 | 2023-04-07 | 西南石油大学 | 一种模拟缝洞型碳酸盐岩凝析气藏生产动态的实验装置及方法 |
US10845354B2 (en) | 2018-05-21 | 2020-11-24 | Newpark Drilling Fluids Llc | System for simulating in situ downhole drilling conditions and testing of core samples |
CN110778312B (zh) * | 2019-10-09 | 2022-08-30 | 东北石油大学 | 一种模拟气藏边底水侵入的模型以及求取水侵系数的方法 |
CN110924908B (zh) * | 2019-11-08 | 2021-10-08 | 中国石油大学(华东) | 一种水驱油藏注采参数确定方法及计算机可读存储介质 |
US11586790B2 (en) | 2020-05-06 | 2023-02-21 | Saudi Arabian Oil Company | Determining hydrocarbon production sweet spots |
CN112069690B (zh) * | 2020-09-11 | 2024-03-08 | 中海石油(中国)有限公司 | 一种深水断块油藏长水平井多级油嘴测试产能的评价方法 |
CN113065261B (zh) * | 2021-04-25 | 2024-01-02 | 中国长江三峡集团有限公司 | 基于水热耦合模拟的地热资源回收率的评价方法 |
US20230102461A1 (en) * | 2021-09-24 | 2023-03-30 | Saudi Arabian Oil Company | Estimating well downtime factor in field modeling |
CN114001654B (zh) * | 2021-11-01 | 2024-03-26 | 北京卫星制造厂有限公司 | 工件端面位姿评价方法 |
US11613957B1 (en) | 2022-01-28 | 2023-03-28 | Saudi Arabian Oil Company | Method and system for high shut-in pressure wells |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1729153A1 (fr) * | 2005-06-02 | 2006-12-06 | Institut Français du Pétrole | Méthode pour simuler les écoulements de fluides au sein d'un réservoir à l'aide d'une discrétisation de type chimère |
US20080319726A1 (en) * | 2007-06-19 | 2008-12-25 | Schlumberger Technology Corporation | System and method for performing oilfield simulation operations |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2801710B1 (fr) * | 1999-11-29 | 2002-05-03 | Inst Francais Du Petrole | Methode pour generer un maillage hybride permettant de modeliser une formation heterogene traversee par un ou plusieurs puits |
WO2007061618A2 (fr) * | 2005-11-22 | 2007-05-31 | Exxonmobil Upstream Research Company | Systeme et procede de simulation |
EP2038809B1 (fr) * | 2006-06-26 | 2019-07-03 | Exxonmobil Upstream Research Company | Procédé de comparaison et de réallocation de production |
US8775141B2 (en) * | 2007-07-02 | 2014-07-08 | Schlumberger Technology Corporation | System and method for performing oilfield simulation operations |
CA2696638C (fr) * | 2010-03-16 | 2012-08-07 | Exxonmobil Upstream Research Company | Utilisation d'une emulsion dont la phase externe est un solvant pour la recuperation in situ de petrole |
-
2009
- 2009-05-20 FR FR0902533A patent/FR2945879B1/fr not_active Expired - Fee Related
-
2010
- 2010-04-01 EP EP10290174.1A patent/EP2253797B1/fr active Active
- 2010-04-27 US US12/767,866 patent/US8694297B2/en not_active Expired - Fee Related
- 2010-05-17 CA CA2704060A patent/CA2704060C/fr not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1729153A1 (fr) * | 2005-06-02 | 2006-12-06 | Institut Français du Pétrole | Méthode pour simuler les écoulements de fluides au sein d'un réservoir à l'aide d'une discrétisation de type chimère |
US20080319726A1 (en) * | 2007-06-19 | 2008-12-25 | Schlumberger Technology Corporation | System and method for performing oilfield simulation operations |
Non-Patent Citations (11)
Title |
---|
BOE O ET AL: "On Near Wellbore Modeling and Real Time Reservoir Management", SPE RESERVOIR SIMULATION SYMPOSIUM, SPE, US, no. SPE 66369, 11 February 2001 (2001-02-11), pages 1 - 14, XP007911025 * |
BOE, O.; FLYNN, J.; REISO, E.: "On Near Wellbore Modeling and Real Time Reservoir Management", SPE 66369, 11 February 2001 (2001-02-11) |
DATABASE COMPENDEX [online] ENGINEERING INFORMATION, INC., NEW YORK, NY, US; 14 June 2007 (2007-06-14), PHILLIPS P D ET AL: "Identifying reservoir potential in shale-dominated thinly bedded clastic reservoirs with a near-well-bore modeling approach", XP002561635, Database accession no. E20084611698499 * |
DATABASE COMPENDEX [online] ENGINEERING INFORMATION, INC., NEW YORK, NY, US; 28 February 2007 (2007-02-28), KROGSTAD S ET AL: "Multiscale mixed-finite-element modeling of coupled wellbore/near-well flow", XP002561634, Database accession no. E20072210623980 * |
DING, Y.; RENARD, G.: "Evaluation of Horizontal Well Performance after Drilling Induced Formation Damage", J. OF ENERGY RESOURCES TECHNOLOGY, vol. 127, September 2005 (2005-09-01) |
FLANDRIN, N.; BENNIS, C.; BOROUCHAKI, H.: "3D Hybrid Mesh Generation for Reservoir Simulation", ECMOR, 30 August 2004 (2004-08-30) |
GAIFFE, S.: "Maillages Hybrides et Décomposition de Domaine pour la Modélisation des Réservoirs Pétroliers", THÈSE DE DOCTORAT, 2000 |
JOHANSEN ET AL: "Iterative techniques in modeling of multi-phase flow in advanced wells and the near well region", JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING, ELSEVIER, AMSTERDAM, NL, vol. 58, no. 1-2, 24 July 2007 (2007-07-24), pages 49 - 67, XP022166386, ISSN: 0920-4105 * |
MLACNIK M ET AL: "Using Well Windows in Full Field Reservoir Simulation", SPE RESERVOIR SIMULATION SYMPOSIUM, SPE, US, no. SPE 66371, 11 February 2001 (2001-02-11), pages 1 - 8, XP007911024 * |
MLACNIK, M.J.; HEINEMANN, Z.E.: "Using well windows in full field reservoir simulation", SPE RESERVOIR SIMULATION SYMPOSIUM, February 2001 (2001-02-01) |
YUGUANG CHEN ET AL: "Upscaled modeling of well singularity for simulating flow in heterogeneous formations", COMPUTATIONAL GEOSCIENCES, KLUWER ACADEMIC PUBLISHERS, DO, vol. 12, no. 1, 4 January 2008 (2008-01-04), pages 29 - 45, XP019571351, ISSN: 1573-1499 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2989200A1 (fr) * | 2012-04-10 | 2013-10-11 | IFP Energies Nouvelles | Procede de selection des positions de puits a forer pour l'exploitation d'un gisement petrolier |
EP2650471A1 (fr) * | 2012-04-10 | 2013-10-16 | IFP Energies nouvelles | Procédé de sélection des positions de puits a forer pour l'exploitation d'un gisement pétrolier |
US9411915B2 (en) | 2012-04-10 | 2016-08-09 | Ipf Energies Nouvelles | Method of selecting positions of wells to be drilled for petroleum reservoir development |
FR2997721A1 (fr) * | 2012-11-08 | 2014-05-09 | Storengy | Radonip : nouvelle methodologie de determination des courbes de productivite des puits d'exploitation de stockages et gisements de fluides compressibles |
WO2014072627A1 (fr) * | 2012-11-08 | 2014-05-15 | Storengy | Nouvelle methodologie de determination des courbes de productivite des puits d'exploitation de stockages et gisements de fluides compressibles |
CN104981585A (zh) * | 2012-11-08 | 2015-10-14 | 斯多恩吉公司 | 确定开发贮藏和储藏可压缩流体的井生产率曲线的新方法 |
CN104981585B (zh) * | 2012-11-08 | 2018-09-21 | 斯多恩吉公司 | 确定开发贮藏和储藏可压缩流体的井生产率曲线的新方法 |
FR3002270A1 (fr) * | 2013-02-21 | 2014-08-22 | IFP Energies Nouvelles | Procede d'exploitation d'un reservoir geologique au moyen d'un modele de reservoir cale et coherent vis a vis des proprietes d'ecoulement |
EP2770162A1 (fr) * | 2013-02-21 | 2014-08-27 | IFP Energies nouvelles | Procédé d'exploitation d'un réservoir géologique au moyen d'un modèle de réservoir calé et cohérent vis à vis des propriétés d'écoulement |
WO2021118714A1 (fr) * | 2019-12-11 | 2021-06-17 | Exxonmobil Upstream Research Company | Méthodologie de semi-élimination permettant de simuler des caractéristiques d'écoulement élevées dans un réservoir |
Also Published As
Publication number | Publication date |
---|---|
CA2704060C (fr) | 2018-02-27 |
CA2704060A1 (fr) | 2010-11-20 |
US8694297B2 (en) | 2014-04-08 |
FR2945879B1 (fr) | 2011-06-24 |
FR2945879A1 (fr) | 2010-11-26 |
US20100299125A1 (en) | 2010-11-25 |
EP2253797B1 (fr) | 2020-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2253797B1 (fr) | Méthode d'exploitation de milieu poreux au moyen d'une modélisation d'écoulements de fluide | |
EP2072752B1 (fr) | Méthode pour optimiser l'exploitation d'un gisement de fluide par prise en compte d'un terme d'échange géologique et transitoire entre blocs matriciels et fractures | |
EP2212722B1 (fr) | Methode de modelisation d'un milieu géologique poreux traversé par un réseau de fractures | |
EP3144468A1 (fr) | Procédé pour caracteriser le réseau de fractures d'un gisement fracture et procédé pour l'exploiter | |
US10941642B2 (en) | Structure for fluid flowback control decision making and optimization | |
FR2961844A1 (fr) | Ecoulement multiphase dans un puits de forage et fracture hydraulique connectee | |
WO2009129060A1 (fr) | Procédé pour déterminer un ensemble de valeurs nettes actualisées pour influencer le forage d’un puits et augmenter la production | |
WO2009034253A1 (fr) | Procede, programme et systeme informatique de mise a l'echelle de donnees de modele de reservoir d'hydrocarbure | |
EP3358129A1 (fr) | Procede d'exploitation d'un gisement d'hydrocarbures par injection d'un gaz sous forme de mousse | |
EP2963235B1 (fr) | Procede d'exploitation d'un gisement petrolier a partir d'une technique de positionnement des puits a forer | |
WO2009027598A1 (fr) | Procede, programme et systeme informatique de conciliation de donnees de modele de reservoir d'hydrocarbure | |
Bedrikovetsky et al. | Formation-damage evaluation from nonlinear skin growth during coreflooding | |
EP3199749A1 (fr) | Procede pour l'exploitation d'un gisement de fluide traverse par des fractures au moyen d'une simulation d'ecoulement basee sur un flux d'echange et un facteur correctif | |
EP2469306A2 (fr) | Procédé pour localiser des barrières hydrauliques au sein d'une couche géologique de stockage de gaz | |
FR3000579A1 (fr) | Fractures a segments multiples | |
EP2365359B1 (fr) | Méthode de calage d'historique d'un modèle géologique | |
Delshad et al. | A mixed-wet hysteretic relative permeability and capillary pressure model for reservoir simulations | |
Delaplace et al. | Reservoir simulations of a polymer flood pilot in the pelican lake heavy oil field (Canada): A step forward | |
EP3192964B1 (fr) | Procede de production d'hydrocarbures comportant un index de productivite des puits sous effet thermique | |
CN109726450B (zh) | 确定页岩气藏水平井重复压裂时机的方法及设备 | |
EP3763913A1 (fr) | Procede d'exploitation d'un gisement d'hydrocarbures par injection d'un gaz sous forme de mousse | |
FR3040509A1 (fr) | ||
Kazemi Nia Korrani et al. | Upscaling low-salinity benefit from laboratory scale to field scale: An ensemble of models with a relative permeability uncertainty range | |
Gamal et al. | An approach for determination of the economically optimal production controlling parameters from water drive oil reservoirs | |
Peshcherenko et al. | Effective Modeling of Stimulation and Production Decline From Tight Naturally Fractured Carbonate Reservoirs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO SE SI SK SM TR |
|
AX | Request for extension of the european patent |
Extension state: AL BA ME RS |
|
17P | Request for examination filed |
Effective date: 20110524 |
|
RAP1 | Party data changed (applicant data changed or rights of an application transferred) |
Owner name: IFP ENERGIES NOUVELLES |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: EXAMINATION IS IN PROGRESS |
|
17Q | First examination report despatched |
Effective date: 20171204 |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: GRANT OF PATENT IS INTENDED |
|
INTG | Intention to grant announced |
Effective date: 20190919 |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE PATENT HAS BEEN GRANTED |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO SE SI SK SM TR |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D Free format text: NOT ENGLISH |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: EP |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R096 Ref document number: 602010063139 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: AT Ref legal event code: REF Ref document number: 1235173 Country of ref document: AT Kind code of ref document: T Effective date: 20200315 |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: FG4D Free format text: LANGUAGE OF EP DOCUMENT: FRENCH |
|
REG | Reference to a national code |
Ref country code: NO Ref legal event code: T2 Effective date: 20200219 |
|
REG | Reference to a national code |
Ref country code: NL Ref legal event code: MP Effective date: 20200219 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: FI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: FR Payment date: 20200429 Year of fee payment: 11 Ref country code: NO Payment date: 20200422 Year of fee payment: 11 |
|
REG | Reference to a national code |
Ref country code: LT Ref legal event code: MG4D |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IS Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200619 Ref country code: BG Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200519 Ref country code: HR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: SE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: LV Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: GR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200520 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: GB Payment date: 20200429 Year of fee payment: 11 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: NL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: CZ Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: ES Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: LT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: PT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200712 Ref country code: SK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: RO Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: SM Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: EE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: DK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R119 Ref document number: 602010063139 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: AT Ref legal event code: MK05 Ref document number: 1235173 Country of ref document: AT Kind code of ref document: T Effective date: 20200219 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MC Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: PL |
|
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
26N | No opposition filed |
Effective date: 20201120 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: LU Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20200401 Ref country code: LI Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20200430 Ref country code: AT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: CH Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20200430 Ref country code: IT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: DE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20201103 |
|
REG | Reference to a national code |
Ref country code: BE Ref legal event code: MM Effective date: 20200430 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: SI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: BE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20200430 Ref country code: PL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20200401 |
|
REG | Reference to a national code |
Ref country code: NO Ref legal event code: MMEP |
|
GBPC | Gb: european patent ceased through non-payment of renewal fee |
Effective date: 20210401 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: NO Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210430 Ref country code: GB Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210401 Ref country code: FR Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210430 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: TR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: MT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: CY Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |