EP2453106B1 - Method for characterising the network of fractures of a fractured deposit and method for exploiting same - Google Patents
Method for characterising the network of fractures of a fractured deposit and method for exploiting same Download PDFInfo
- Publication number
- EP2453106B1 EP2453106B1 EP11290475.0A EP11290475A EP2453106B1 EP 2453106 B1 EP2453106 B1 EP 2453106B1 EP 11290475 A EP11290475 A EP 11290475A EP 2453106 B1 EP2453106 B1 EP 2453106B1
- Authority
- EP
- European Patent Office
- Prior art keywords
- zone
- fracture
- zones
- fractures
- fin
- 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.)
- Not-in-force
Links
- 238000000034 method Methods 0.000 title claims description 57
- 230000035699 permeability Effects 0.000 claims description 51
- 239000012530 fluid Substances 0.000 claims description 33
- 238000004519 manufacturing process Methods 0.000 claims description 31
- 238000004088 simulation Methods 0.000 claims description 30
- 238000012360 testing method Methods 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 16
- 238000009826 distribution Methods 0.000 claims description 7
- 206010017076 Fracture Diseases 0.000 description 204
- 208000010392 Bone Fractures Diseases 0.000 description 105
- 238000004364 calculation method Methods 0.000 description 19
- 238000013459 approach Methods 0.000 description 7
- 230000008901 benefit Effects 0.000 description 7
- 229930195733 hydrocarbon Natural products 0.000 description 7
- 150000002430 hydrocarbons Chemical class 0.000 description 7
- 238000005259 measurement Methods 0.000 description 7
- 238000002347 injection Methods 0.000 description 6
- 239000007924 injection Substances 0.000 description 6
- 238000012512 characterization method Methods 0.000 description 5
- 230000004907 flux Effects 0.000 description 5
- 238000010276 construction Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000011084 recovery Methods 0.000 description 3
- 230000001052 transient effect Effects 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 241000897276 Termes Species 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- VYMDGNCVAMGZFE-UHFFFAOYSA-N phenylbutazonum Chemical compound O=C1C(CCCC)C(=O)N(C=2C=CC=CC=2)N1C1=CC=CC=C1 VYMDGNCVAMGZFE-UHFFFAOYSA-N 0.000 description 2
- 230000000135 prohibitive effect Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 240000008042 Zea mays Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 229940082150 encore Drugs 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- 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
-
- 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/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/008—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells by injection test; by analysing pressure variations in an injection or production test, e.g. for estimating the skin factor
Definitions
- the present invention relates to the field of the exploitation of underground deposits, such as hydrocarbon deposits comprising a network of fractures.
- the invention relates to a method for characterizing the fracture network and thus to construct a representation of the deposit.
- the invention also relates to a method using this representation to optimize the management of such an operation by means of a prediction of fluid flows likely to occur through this medium, to simulate a production of hydrocarbons according to various production scenarios. .
- the oil industry and more specifically the exploration and exploitation of deposits, especially oil, require the acquisition of the best possible knowledge of the underground geology to efficiently provide an assessment of reserves, a modeling of production , or the management of the operation.
- the determination of the location of a production well or an injection well, the constitution of the drilling mud, the completion characteristics, the choice of a hydrocarbon recovery process (such as the injection of water for example) and the parameters necessary for the implementation of this process (such as the injection pressure, the production rate, etc.) require a good knowledge of the deposit. Knowing the deposit means knowing the petrophysical properties of the subsoil at any point in space.
- Cracked reservoirs are an extreme type of heterogeneous reservoirs with two contrasting media, a matrix medium containing most of the oil in place and having a low permeability, and a cracked medium representing less than 1% of the oil in place. and highly conductive.
- the cracked medium itself can be complex, with different sets of cracks characterized by their respective density, length, orientation, inclination and aperture.
- a "fracture” is a plane discontinuity, very thin in relation to its extension, which represents a plane of rupture of a rock in the deposit.
- the knowledge of the distribution and the behavior of these fractures makes it possible to optimize the localization and the spacing between the wells that one to drill through the oil field.
- the geometry of the fracture network conditions the displacement of the fluids both at the reservoir scale and at the local scale where it determines elementary matrix blocks in which the oil is trapped. Knowing the distribution of fractures, is therefore very useful, also, at a later stage, for the tank engineer who tries to calibrate the models he builds to simulate the deposits in order to reproduce or predict past production curves. or future.
- geoscientists have three-dimensional images of the deposits, making it possible to locate a large number of fractures.
- the tank engineering specialist implements a calculation software, called a “reservoir simulator” (or “flow simulator”), which calculates the flows and the evolution of the pressures within the reservoir represented by the reservoir model. tank.
- a “reservoir simulator” or “flow simulator”
- the results of these calculations allow it to predict and optimize the deposit in terms of flow and / or quantity of recovered hydrocarbons.
- the calculation of the behavior of the reservoir according to a given production scenario constitutes a "reservoir simulation”.
- This simplified representation is called “double-medium approach”, it is proposed by Warren JE et al. in "The Behavior of Naturally Fractured Reservoirs", SPE Journal (September 1963), 245-255 .
- This technique consists in considering the fractured medium as two continua exchanging fluids between them: matrix blocks and fractures. This is called the “double medium” or “double porosity” model.
- the "double medium” modeling of a fractured deposit consists in discretizing this deposit into two sets of meshes (called grids) superimposed, constituting the grid "crack” and the grid “matrix”.
- Each elementary volume of the fractured deposit is thus conceptually represented by two meshes, one "crack” and the other "matrix", coupled together (that is to say, exchanging fluids).
- these two meshes represent the set of matrix blocks delimited by fractures present in this place of the reservoir.
- meshes have lateral dimensions MF (commonly 100 or 200 m) given the size of the fields and limited possibilities of simulation software in terms of capacity and calculation time.
- the fracture reservoir elemental volume contains innumerable fractures forming a complex network delimiting multiple matrix blocks. of variable dimensions and shapes depending on the geological context.
- Each of the constituent real blocks exchanges fluids with the fractures that surround it at a rate (flow) that is specific to it because of the size and shape of this particular block.
- the approach consists, for each elementary volume (mesh) of reservoir, to represent the real fractured medium as a set of matricial blocks all identical, parallelepipedic, delimited by an orthogonal and regular network of fractures oriented along the main directions of flow: for each mesh, the so-called “equivalent” permeabilities of this fracture network are thus determined and defines a "representative" matrix block (of the real (geological) distribution of the blocks), unique and of parallelepipedal shape. It is then possible to formulate and calculate the matrix-crack exchange fluxes for this "representative" block, and to multiply the result by the number of such blocks in the elementary volume (mesh) to obtain the flux at the scale of this mesh.
- Another method, described in the patent application US 2008/0091396 A1 relates to a method of modeling hydraulic fractures, in which zones are built in the region of the well approaches.
- the invention involves a simplification of the fracture networks at the local scale of the well drainage area, in order to be able to simulate the wells tests of fractured reservoirs and thus to calibrate the conductivities of the families of fractures.
- This hydraulic calibration of fractures results in a set of parameters characterizing the fracture network (or fracture model).
- This fracture model is then used to construct a reservoir-scale dual-flow model.
- the statistical parameters can be chosen from the following parameters: fracture density, fracture length, fracture orientation, fracture inclination, fracture opening, and distribution of fractures within the deposit.
- a shape ratio for each zone is determined so as to reproduce an anisotropy of flow around the well, and the areas so as to respect the form report.
- This aspect ratio can be determined by the main values of the permeability tensor.
- a distance between the boundaries between zones is defined, so as to give a weight equal to each zone in terms of the pressure difference observed on each zone in steady-state flow regime.
- This distance can be defined by fixing the lengths of one of the two axes of two successive ellipses to values in geometric progression of constant reason.
- three zones are constructed, a first zone (ZNS) containing the well within which no simplification of the image is made, a second zone (ZP) in contact with the first zone within of which a first simplification of the image is performed, and a third zone (ZL) in contact with the second zone in which a second simplification of the image is performed, the second simplification being more important than the first simplification.
- ZNS first zone
- ZP second zone
- ZL third zone
- modelizations are intended to provide a description of the deposit, characterized by the structure / geometry and the petrophysical properties of the deposits or geological formations that constitute it.
- modelizations are based on a representation of the deposit, in a set of meshes.
- Each of these meshes represents a given volume of the deposit, and constitutes an elementary volume of the deposit.
- the set of meshes is a discrete representation of the deposit, called reservoir model.
- the figure 2 illustrates a two-dimensional view of a tank model. Fractures are represented by lines. Meshes are not represented.
- the specialists in charge of the statistical characterization of the reservoir perform direct and indirect observations (OF) of the reservoir. For this, they have 1) well cores extracted from the reservoir on which a statistical study of the intersected fractures is carried out, 2) characteristic outcrops of the reservoir which have the advantage of providing a large-scale vision of the fracture network 3 ) seismic images allowing them to identify large geological events.
- OF direct and indirect observations
- the objective of the characterization and modeling of the deposit fracture network is to propose a model of fractures validated on the local flows around the wells. This fracture model is then extended to the reservoir scale in order to perform production simulations. To do this, flow properties are associated with each mesh of the reservoir model (MR) (permeability tensor, porosity) of the two media (fracture and matrix).
- MR reservoir model
- PSF statistical parameters
- DNN discrete fracture network
- FIG. 2 illustrates a realization of a network of fractures / faults at the scale of a reservoir.
- Each mesh of the reservoir model thus represents a discrete network of fractures delimiting a set of porous matrix blocks of irregular shapes and sizes delimited by fractures.
- Such an image is represented on the figure 3 .
- This discrete fracture network constitutes a representative image of the real network of fractures delimiting the matrix blocks.
- the next step is to determine the flow properties of the initial fractures (C f , e) and then calibrate these properties by means of well-test simulations on discrete, realistic image-inherited local flow models.
- the actual (geological) fracture network at the reservoir scale Although only covering a limited area (drainage area) around the well, such a well testing simulation model still has a large number of computational nodes if the fracture network is dense. As a result, the size of the systems to be solved and / or the duration of the calculations often remain prohibitive. Hence the need for a procedure to simplify the fracture network.
- the fracture network obtained in the previous step and representative of the actual fractured reservoir can not be used to simulate, i.e. reproduce and / or predict, the local flows around the well.
- the invention is not limited to the definition of three zones. It is also possible to divide the domain into n zones, the simplification of the network being increased from the zone 1 (ZNS) to the zone n (the furthest away from the well). It is thus possible to create a zone ZNS, n1 zones of type ZP, and n2 zones of type ZL.
- fracture network modeling is to simulate well-flow responses (transient or pseudo-permanent flow tests, interference, flowmetry, etc.). This involves simulating, for example, the production of oil via each well drilled through the reservoir.
- one of the two preceding methods may be chosen so as to optimize the precision and speed of the calculations, by applying the method described in the application of EP Patent 2,037,080 , based on the calculation of a connectivity index.
- This technique makes it possible, as a preliminary step, to determine the permeability tensors of a few cells of the reservoir model surrounding the well, and considered representative of the flow of the ZNS.
- the diagonalization of these permeability tensors provides the eigenvectors oriented along the main directions of flow sought.
- the delimitation of the sub-zones is based on the comparison of the equivalent permeabilities calculated on the neighboring "blocks" materializing these arcs.
- the analytical method for calculating equivalent permeabilities is preferably used because it is in this case to determine if a block has the same dynamic behavior as the neighboring block.
- the analytical approach provides erroneous results, the errors are systematic, similar from one block to the other, which allows the comparison of the results between blocks which, let us specify it, does not require a great precision considering the simple objective of definition of zone.
- the analytical approach is quite justified, with the considerable advantage of allowing much faster calculations, guaranteeing a practical feasibility.
- the scaling calculations will allow the fracture network of these sub-zones to be replaced by a simplified network having the same flow properties as the original network. In this case, and contrary to the above, the calculation of the equivalent fracture permeability tensor must be as accurate as possible.
- the fractures of the simplified networks obtained are extended beyond the limits of sub-zones (or “blocks") to ensure sufficient partial “overlap” and thus sufficient horizontal connectivity of the simplified sub-area networks (or “blocks” of ZP) nearby.
- the fractures of the simplified network can thus be extended by a length equal to 60% of the maximum spacing (s 1 large , s 2 large ) fractures of this network.
- the next step is the calibration of the fracture flow properties (conductivity and opening of fractures), locally around the wells. This requires the simulation of well tests. According to the invention, this simulation of well tests is performed on the simplified flow models ( figure 6 ).
- the tank engineer can use the methods described in the following documents, applied to the entire tank this time: FR 2,757,947 ( US 6,023,656 ) FR 2 757 957 ( US 6,064,944 ) and EP 2,037,080 . These methods make it possible to calculate the equivalent fracture permeabilities and the equivalent block dimensions for each of the cells of the reservoir model.
- the reservoir engineer chooses a production process, for example the water injection recovery process, the optimal implementation scenario for the field in question.
- the definition of an optimal water injection scenario will consist, for example, in determining the number and location (position and spacing) of the injection and production wells in order to best take into account the impact of the fractures on the water. progression of fluids within the reservoir.
- the simulator solves all the equations specific to each mesh and each of the two grids of the model (equations involving the matrix-crack exchange formula described above) and thus delivers the solution values of the unknowns S (t) (saturations, pressures, concentrations, temperature, etc.) at this instant t. From this resolution, comes the knowledge of the quantities of oil produced and the state of the deposit (distribution of pressures, saturations, etc %) at the moment considered.
Description
La présente invention concerne le domaine de l'exploitation de gisements souterrains, tels que des gisements d'hydrocarbures comportant un réseau de fractures.The present invention relates to the field of the exploitation of underground deposits, such as hydrocarbon deposits comprising a network of fractures.
En particulier, l'invention concerne une méthode pour caractériser le réseau de fractures et ainsi construire une représentation du gisement. L'invention concerne également une méthode utilisant cette représentation pour optimiser la gestion d'une telle exploitation au moyen d'une prédiction des écoulements de fluides susceptibles de se produire à travers ce milieu, pour simuler une production d'hydrocarbures suivant divers scénarios de production.In particular, the invention relates to a method for characterizing the fracture network and thus to construct a representation of the deposit. The invention also relates to a method using this representation to optimize the management of such an operation by means of a prediction of fluid flows likely to occur through this medium, to simulate a production of hydrocarbons according to various production scenarios. .
L'industrie pétrolière, et plus précisément l'exploration et l'exploitation de gisements, notamment pétroliers, nécessitent d'acquérir une connaissance aussi parfaite que possible de la géologie souterraine pour fournir de façon efficace une évaluation des réserves, une modélisation de la production, ou la gestion de l'exploitation. En effet, la détermination de l'emplacement d'un puits de production ou d'un puits d'injection, la constitution de la boue de forage, les caractéristiques de complétion, le choix d'un procédé de récupération des hydrocarbures (tel que l'injection d'eau par exemple) et des paramètres nécessaires à la mise en oeuvre de ce procédé (tels que la pression d'injection, le débit de production,...) nécessitent de bien connaître le gisement. Connaître le gisement signifie notamment connaître les propriétés pétrophysiques du sous-sol en tout point de l'espace.The oil industry, and more specifically the exploration and exploitation of deposits, especially oil, require the acquisition of the best possible knowledge of the underground geology to efficiently provide an assessment of reserves, a modeling of production , or the management of the operation. In fact, the determination of the location of a production well or an injection well, the constitution of the drilling mud, the completion characteristics, the choice of a hydrocarbon recovery process (such as the injection of water for example) and the parameters necessary for the implementation of this process (such as the injection pressure, the production rate, etc.) require a good knowledge of the deposit. Knowing the deposit means knowing the petrophysical properties of the subsoil at any point in space.
Pour ce faire, depuis longtemps, l'industrie pétrolière allie les mesures sur champ (in situ) aux modélisations expérimentales (réalisées au laboratoire) et/ou numériques (réalisées au moyen de logiciels). Les modélisations des gisements pétroliers constituent donc une étape technique indispensable à toute exploration ou exploitation de gisement. Ces modélisations ont pour but de fournir une description du gisement.To do this, the oil industry has for a long time been combining on-field measurements (in situ) with experimental (laboratory-generated) and / or digital (software-based) models. Modeling oil fields is therefore a technical step essential to any exploration or exploitation of deposits. These models are intended to provide a description of the deposit.
Les réservoirs fissurés constituent un type extrême de réservoirs hétérogènes comportant deux milieux contrastés, un milieu matriciel contenant la plus grande part de l'huile en place et présentant une faible perméabilité, et un milieu fissuré représentant moins de 1 % de l'huile en place et hautement conducteur. Le milieu fissuré lui-même peut être complexe, avec différents ensembles de fissures caractérisés par leur densité, longueur, orientation, inclinaison et ouverture respectives.Cracked reservoirs are an extreme type of heterogeneous reservoirs with two contrasting media, a matrix medium containing most of the oil in place and having a low permeability, and a cracked medium representing less than 1% of the oil in place. and highly conductive. The cracked medium itself can be complex, with different sets of cracks characterized by their respective density, length, orientation, inclination and aperture.
Les spécialistes en charge de l'exploitation de réservoirs fracturés, ont besoin de parfaitement connaître le rôle des fractures. On appelle "fracture", une discontinuité plane, de très faible épaisseur par rapport à son extension, et qui représente un plan de rupture d'une roche du gisement. D'une part, la connaissance de la distribution et du comportement de ces fractures permet d'optimiser la localisation et l'espacement entre les puits que l'on compte forer au travers du gisement pétrolifère. D'autre part, la géométrie du réseau de fractures conditionne le déplacement des fluides tant à l'échelle du réservoir qu'à l'échelle locale où elle détermine des blocs matriciels élémentaires dans lesquels l'huile est piégée. Connaître la distribution des fractures, est donc très utile, aussi, à un stade ultérieur, pour l'ingénieur de réservoir qui cherche à calibrer les modèles qu'il construit pour simuler les gisements afin d'en reproduire ou prédire les courbes de production passées ou futures. A ces fins, les spécialistes de géosciences disposent d'images tridimensionnelles des gisements, permettant de localiser un grand nombre de fractures.The specialists in charge of the operation of fractured reservoirs, need to know perfectly the role of the fractures. A "fracture" is a plane discontinuity, very thin in relation to its extension, which represents a plane of rupture of a rock in the deposit. On the one hand, the knowledge of the distribution and the behavior of these fractures makes it possible to optimize the localization and the spacing between the wells that one to drill through the oil field. On the other hand, the geometry of the fracture network conditions the displacement of the fluids both at the reservoir scale and at the local scale where it determines elementary matrix blocks in which the oil is trapped. Knowing the distribution of fractures, is therefore very useful, also, at a later stage, for the tank engineer who tries to calibrate the models he builds to simulate the deposits in order to reproduce or predict past production curves. or future. For these purposes, geoscientists have three-dimensional images of the deposits, making it possible to locate a large number of fractures.
Ainsi, pour reproduire ou prédire (i.e. "simuler") la production d'hydrocarbures lors de la mise en production d'un gisement suivant un scénario de production donné (caractérisé par la position des puits, la méthode de récupération, ...), le spécialiste en ingénierie de réservoir met en oeuvre un logiciel de calcul, appelé « simulateur de réservoir » (ou « simulateur d'écoulement »), qui calcule les écoulements et l'évolution des pressions au sein du réservoir représenté par le modèle de réservoir. Les résultats de ces calculs lui permettent de prévoir et d'optimiser le gisement en termes de débit et/ou de quantité d'hydrocarbures récupérés. Le calcul du comportement du réservoir suivant un scénario de production donné constitue une « simulation de réservoir ».Thus, to reproduce or predict (ie "simulate") the production of hydrocarbons during the production of a deposit according to a given production scenario (characterized by the position of the wells, the recovery method, ...) , the tank engineering specialist implements a calculation software, called a "reservoir simulator" (or "flow simulator"), which calculates the flows and the evolution of the pressures within the reservoir represented by the reservoir model. tank. The results of these calculations allow it to predict and optimize the deposit in terms of flow and / or quantity of recovered hydrocarbons. The calculation of the behavior of the reservoir according to a given production scenario constitutes a "reservoir simulation".
On connaît une méthode pour optimiser l'exploitation d'un gisement de fluide traversé par un réseau de fractures, dans laquelle on simule des écoulements de fluides dans le gisement au moyen d'une modélisation simplifiée mais réaliste du gisement. Cette représentation simplifiée est appelée "approche double milieu", elle est proposée par
Face à une telle complexité géométrique du milieu réel, la démarche consiste, pour chaque volume élémentaire (maille) de réservoir, à représenter le milieu fracturé réel comme un ensemble de blocs matriciels tous identiques, parallélépipédiques, délimités par un réseau orthogonal et régulier de fractures orientées suivant les directions principales d'écoulement : pour chaque maille, on détermine ainsi les perméabilités dites "équivalentes" de ce réseau de fractures et définit un bloc matriciel, dit « représentatif » (de la distribution réelle (géologique) des blocs), unique et de forme parallélépipédique. Il est alors possible de formuler et calculer les flux d'échange matrice-fissure pour ce bloc "représentatif", et d'en multiplier le résultat par le nombre de tels blocs dans le volume élémentaire (maille) pour obtenir le flux à l'échelle de cette maille.Faced with such a geometrical complexity of the real environment, the approach consists, for each elementary volume (mesh) of reservoir, to represent the real fractured medium as a set of matricial blocks all identical, parallelepipedic, delimited by an orthogonal and regular network of fractures oriented along the main directions of flow: for each mesh, the so-called "equivalent" permeabilities of this fracture network are thus determined and defines a "representative" matrix block (of the real (geological) distribution of the blocks), unique and of parallelepipedal shape. It is then possible to formulate and calculate the matrix-crack exchange fluxes for this "representative" block, and to multiply the result by the number of such blocks in the elementary volume (mesh) to obtain the flux at the scale of this mesh.
Il convient cependant de noter que le calcul des perméabilités équivalentes nécessite de connaître les propriétés d'écoulement (i.e. les conductivités) des fractures discrètes du modèle géologique.It should however be noted that the calculation of the equivalent permeabilities requires knowing the flow properties (ie. The conductivities) discrete fracture of the geological model.
C'est pourquoi, préalablement à la construction de ce modèle de réservoir équivalent (dit "modèle de réservoir double milieu") tel que décrit plus haut, il est d'abord nécessaire de simuler les réponses en écoulement de quelques puits (tests d'écoulement transitoires ou pseudo-permanents, interférences, débitmétrie, etc.) sur des modèles extraits du modèle géologique donnant une représentation discrète (réaliste) des fractures alimentant ces puits. L'ajustement des réponses en pression/débit simulées sur les mesures de champ permet de calibrer les conductivités des familles de fractures. Bien que ne couvrant qu'une aire limitée (aire de drainage) autour du puits, un tel modèle de simulation de tests de puits comporte encore de très nombreux noeuds de calcul si le réseau de fractures est dense. En conséquence, la taille des systèmes à résoudre et/ou la durée des calculs demeurent souvent prohibitives.This is why, prior to the construction of this equivalent reservoir model (called the "double-medium tank model") as described above, it is first necessary to simulate the flow responses of a few wells (test results). transient or pseudo-permanent flow, interferences, flowmetry, etc.) on models extracted from the geological model giving a discrete (realistic) representation of the fractures feeding these wells. Adjusting the simulated pressure / flow responses on the field measurements calibrates the conductivities of the families of fractures. Although only covering a limited area (drainage area) around the well, such a well testing simulation model still has a large number of computational nodes if the fracture network is dense. As a result, the size of the systems to be solved and / or the duration of the calculations often remain prohibitive.
Une autre méthode, décrite dans la demande de brevet
Pour surmonter cette difficulté, l'invention comporte une simplification des réseaux de fractures à l'échelle locale de l'aire de drainage du puits, afin de pouvoir simuler les tests de puits de réservoirs fracturés et calibrer ainsi les conductivités des familles de fractures. Cette calibration hydraulique des fractures aboutit à un jeu de paramètres caractérisant le réseau de fractures (ou modèle de fractures). Ce modèle de fracture est par la suite utilisé pour construire un modèle d'écoulement double milieu à l'échelle du réservoir.To overcome this difficulty, the invention involves a simplification of the fracture networks at the local scale of the well drainage area, in order to be able to simulate the wells tests of fractured reservoirs and thus to calibrate the conductivities of the families of fractures. This hydraulic calibration of fractures results in a set of parameters characterizing the fracture network (or fracture model). This fracture model is then used to construct a reservoir-scale dual-flow model.
Ainsi, l'objet de l'invention concerne une méthode pour optimiser l'exploitation d'un gisement de fluide traversé par un réseau de fractures et par au moins un puits, dans laquelle on construit une représentation dudit gisement de fluide, on discrétise ledit gisement en un ensemble de mailles, et on caractérise les fractures par des paramètres statistiques (PSF) à partir d'observations dudit gisement. La méthode comporte les étapes suivantes :
- a) on déduit desdits paramètres statistiques (PSF) un tenseur de perméabilité équivalente et une ouverture moyenne desdites fractures, à partir desquels on construit une image représentative du réseau de fractures délimitant des blocs poreux et des fractures ;
- b) on déduit dudit tenseur une direction d'écoulement de fluide autour dudit puits ;
- c) on définit autour dudit puits, une première zone de frontière elliptique centrée sur ledit puits et contenant ledit puits, et au moins une seconde zone de frontière elliptique centrée sur ledit puits et formant une couronne elliptique avec la frontière elliptique de ladite première zone, lesdites zones étant orientées selon ladite direction d'écoulement de fluide ;
- d) on simplifie ladite image représentative du réseau de fractures de façon différente dans chacune desdites zones ;
- e) on utilise ladite image simplifiée pour construire ladite représentation du gisement de fluide ;
- f) on utilise ladite représentation du gisement de fluide et un simulateur d'écoulement pour optimiser l'exploitation dudit gisement de fluide.
- a) deducing from said statistical parameters (PSF) an equivalent permeability tensor and an average opening of said fractures, from which a representative image of the fracture network delimiting porous blocks and fractures is constructed;
- b) deducing from said tensor a direction of fluid flow around said well;
- c) defining around said well, a first elliptical boundary zone centered on said well and containing said well, and at least one second elliptical boundary zone centered on said well and forming an elliptical crown with the elliptical boundary of said first zone, said zones being oriented in said direction of fluid flow;
- d) simplifying said representative image of the fracture network differently in each of said zones;
- e) using said simplified image to construct said representation of the fluid reservoir;
- f) using said representation of the fluid reservoir and a flow simulator to optimize the exploitation of said fluid reservoir.
Selon l'invention, les paramètres statistiques (PSF) peuvent être choisis parmi les paramètres suivants : densité de fractures, longueur de fractures, orientation de fractures, inclinaison de fractures, ouverture de fractures, et distribution de fractures au sein du gisement.According to the invention, the statistical parameters (PSF) can be chosen from the following parameters: fracture density, fracture length, fracture orientation, fracture inclination, fracture opening, and distribution of fractures within the deposit.
Selon un mode réalisation, on détermine un rapport de forme pour chaque zone, défini à partir des longueurs des axes de l'ellipse constituant la frontière de la zone, de manière à reproduire une anisotropie d'écoulement autour du puits, et on construit les zones de façon à respecter le rapport de forme. On peut déterminer ce rapport de forme au moyen des valeurs principales du tenseur de perméabilité.According to one embodiment, a shape ratio for each zone, defined from the lengths of the axes of the ellipse constituting the zone boundary, is determined so as to reproduce an anisotropy of flow around the well, and the areas so as to respect the form report. This aspect ratio can be determined by the main values of the permeability tensor.
Selon un mode réalisation, on définit une distance entre les frontières entre zones, de façon à donner un poids égal à chaque zone en termes de différence de pression constatée sur chaque zone en régime d'écoulement permanent. Cette distance peut être définie en fixant les longueurs d'un des deux axes de deux ellipses successives à des valeurs en progression géométrique de raison constante.According to one embodiment, a distance between the boundaries between zones is defined, so as to give a weight equal to each zone in terms of the pressure difference observed on each zone in steady-state flow regime. This distance can be defined by fixing the lengths of one of the two axes of two successive ellipses to values in geometric progression of constant reason.
Selon un autre mode de réalisation, on construit trois zones, une première zone (ZNS) contenant le puits au sein de laquelle aucune simplification de l'image n'est réalisée, une seconde zone (ZP) au contact de la première zone au sein de laquelle une première simplification de l'image est réalisée, et une troisième zone (ZL) au contact de la seconde zone au sein de laquelle une seconde simplification de l'image est réalisée, la seconde simplification étant plus importante que la première simplification.According to another embodiment, three zones are constructed, a first zone (ZNS) containing the well within which no simplification of the image is made, a second zone (ZP) in contact with the first zone within of which a first simplification of the image is performed, and a third zone (ZL) in contact with the second zone in which a second simplification of the image is performed, the second simplification being more important than the first simplification.
De façon avantageuse, les seconde et troisième zones peuvent être découpées en sous-zones en appliquant les étapes suivantes :
- ladite seconde zone est divisée en un nombre de sous-zones égal à un nombre de blocs de mailles présents dans ladite zone, un bloc de mailles désignant un empilement vertical de mailles ;
- ladite troisième zone est divisée en réalisant les étapes suivantes :
- ∘ on découpe tous les degrés la frontière de ladite troisième zone, définissant 360 arcs ;
- ∘ on définit une sous-zone en reliant des points-extrémités de chacun desdits arcs au centre de l'ellipse formant la frontière ;
- ∘ pour chacune desdites sous-zones, on calcule un tenseur de perméabilité de fracture équivalente à partir duquel on détermine une orientation des écoulements dans ladite sous-zone ;
- ∘ on compare les valeurs de perméabilité de fracture équivalente et l'orientation de l'écoulement entre sous-zones voisines ; et
- ∘ on regroupe des sous-zones voisines en une seule sous-zone lorsqu'une différence entre les valeurs de perméabilité est inférieure à un premier seuil et lorsqu'une différence entre les orientations des écoulements est inférieure à un second seuil.
- said second zone is divided into a number of sub-zones equal to a number of mesh blocks present in said zone, a mesh block denoting a vertical stack of meshes;
- said third zone is divided by performing the following steps:
- ∘ the borders of said third zone, defining 360 arcs, are cut off all the degrees;
- A sub-area is defined by connecting end-points of each of said arcs to the center of the ellipse forming the boundary;
- ∘ for each of said sub-zones, an equivalent fracture permeability tensor is calculated from which a direction of the flows in said sub-zone is determined;
- Comparing the values of equivalent fracture permeability and the orientation of the flow between neighboring sub-zones; and
- Grouping neighboring sub-zones into a single sub-zone when a difference between the permeability values is lower than a first threshold and when a difference between the orientations of the flows is lower than a second threshold.
Selon l'invention, on peut simplifier l'image en réalisant les étapes suivantes :
- on construit un réseau de fracture équivalent (RFE) à ladite image, au moyen d'une représentation dite de Warren et Root, dans laquelle le réseau est caractérisé par des espacements de fractures (s1 fin, s2 fin) dans deux directions orthogonales de perméabilité principale, par un paramètre d'ouverture de fractures (efin), par des conductivités de fractures (Cf1 fin et Cf2 fin), et une perméabilité (km fin) d'un milieu matriciel entre fractures ;
- on simplifie ledit réseau de fracture équivalent (RFE) au moyen d'un coefficient d'écartement des fractures (G) du réseau dont la valeur est inférieure à une valeur Gmax-zone définie sur chacune desdites zones afin de garantir une connectivité suffisante entre zones simplifiées et zones non simplifiées.
- constructing an equivalent fracture network (RFE) to said image, by means of a so-called Warren and Root representation, in which the network is characterized by fracture spacings (s 1 end , s 2 end ) in two orthogonal directions main permeability, by a fracture opening parameter (e fin ), by fracture conductivities (C f1 fin and C f2 end ), and a permeability (k m fin ) of a matrix medium between fractures;
- said equivalent fracture network (EFR) is simplified by means of a network gap fracture coefficient (G) whose value is lower than a maximum G -zone value defined on each of said zones in order to guarantee sufficient connectivity between simplified areas and un-simplified areas.
On peut, pour une sous-zone donnée, définir la valeur Gmax-zone de la façon suivante :
- DLM : dimension latérale minimale de la sous-zone donnée ;
- s1 fin, s2 fin : les espacements de fractures dans la représentation dite de Warren et Root.
- DLM: minimum lateral dimension of the given sub-area;
- s 1 end , s 2 fin : fracture spacings in the representation of Warren and Root.
Enfin, l'invention concerne également une méthode pour optimiser la gestion d'un gisement. Elle comporte les étapes suivantes :
- on réitère à l'étape a) en modifiant lesdits paramètres statistiques (PSF) de façon à minimiser une différence entre un résultat de test de puits et un résultat de simulation de test de puits à partir de ladite image simplifiée ;
- on associe à chacune desdites mailles au moins une valeur de perméabilité équivalente et une valeur d'ouverture moyenne desdites fractures, lesdites valeurs étant déterminées à partir desdits paramètres statistiques (PSF) modifiés.
- on simule des écoulements de fluides dans ledit gisement au moyen d'un simulateur d'écoulement et des valeurs de perméabilité équivalente et des valeurs d'ouverture moyenne desdites fractures associées à chacune desdites mailles ;
- on sélectionne un scénario de production permettant d'optimiser la production du gisement à l'aide de ladite simulation des écoulements de fluides; et
- on exploite ledit gisement selon ledit scénario permettant d'optimiser la production du gisement.
- repeating said statistical parameters (PSF) in step a) so as to minimize a difference between a well test result and a well test simulation result from said simplified image;
- each of said meshes is associated with at least one equivalent permeability value and an average opening value of said fractures, said values being determined from said modified statistical parameters (PSF).
- fluid flows are simulated in said reservoir by means of a flow simulator and equivalent permeability values and mean opening values of said fractures associated with each of said meshes;
- a production scenario is selected which makes it possible to optimize the production of the deposit by means of said simulation of the flows of fluids; and
- the said deposit is exploited according to the said scenario making it possible to optimize the production of the deposit.
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.Other characteristics and advantages of the method according to the invention will appear on reading the following description of nonlimiting examples of embodiments, with reference to the appended figures and described below.
-
La
figure 1 illustre les différentes étapes de la méthode selon l'invention.Thefigure 1 illustrates the different steps of the method according to the invention. -
La
figure 2 illustre une réalisation d'un réseau de fractures/failles à l'échelle d'un réservoirThefigure 2 illustrates a realization of a network of fractures / faults at the scale of a reservoir -
La
figure 3 illustre un réseau de fracture discret (DFN) initialThefigure 3 illustrates an initial discrete fracture network (DFN) -
La
figure 4 illustre un réseau de fracture équivalent (RFE), dit de Warren et Root.Thefigure 4 illustrates an equivalent fracture network (RFE), said Warren and Root. -
La
figure 5 illustre la création de zones et sous-zones nécessaires à la simplification du réseau de fracture équivalent (RFE).Thefigure 5 illustrates the creation of zones and sub-areas necessary for the simplification of the equivalent fracture network (RFE). -
La
figure 6 illustre un réseau de fractures équivalent (RFES) simplifié selon l'invention.Thefigure 6 illustrates a simplified fracture network (RFES) simplified according to the invention.
La méthode selon l'invention pour optimiser l'exploitation d'un gisement, utilisant la méthode selon l'invention de caractérisation du réseau de fractures, comporte quatre étapes, comme illustré sur la
- 1- Discrétisation du gisement en un ensemble de mailles (MR)
- 2- Modélisation du réseau de fractures (DFN, RFE, RFES)
- 3- Simulation des écoulements de fluides (SIM) et optimisation des conditions de production du gisement (OPT)
- 4- Exploitation (globale) optimisée du gisement (EXPLO)
- 1- Discretization of the deposit into a set of meshes (MR)
- 2- Modeling the fracture network (DFN, RFE, RFES)
- 3- Simulation of fluid flows (SIM) and optimization of production conditions of the deposit (OPT)
- 4- Optimized Exploitation (Global) of the Deposit (EXPLO)
Depuis longtemps, l'industrie pétrolière allie les mesures sur champ (in situ) aux modélisations expérimentales (réalisées au laboratoire) et/ou numériques (réalisées au moyen de logiciels). Les modélisations des gisements pétroliers, constituent donc une étape technique indispensable à toute exploration ou exploitation de gisement. Ces modélisations ont pour but de fournir une description du gisement, caractérisée par la structure/géométrie et les propriétés pétrophysiques des dépôts ou formations géologiques qui le constituent.For a long time, the oil industry has combined on-field measurements (in situ) with experimental (laboratory-generated) and / or digital (software-based) models. The modeling of the oil fields, therefore constitute a technical step essential to any exploration or exploitation of deposits. These modelizations are intended to provide a description of the deposit, characterized by the structure / geometry and the petrophysical properties of the deposits or geological formations that constitute it.
Ces modélisations se basent sur une représentation du gisement, en un ensemble de mailles. Chacune de ces mailles représente un volume donné du gisement, et constitue un volume élémentaire du gisement. L'ensemble des mailles constitue une représentation discrète du gisement, appelée modèle de réservoir.These modelizations are based on a representation of the deposit, in a set of meshes. Each of these meshes represents a given volume of the deposit, and constitutes an elementary volume of the deposit. The set of meshes is a discrete representation of the deposit, called reservoir model.
Les spécialistes connaissent de nombreux outils logiciels permettant de construire de tel modèle de réservoir, à partir de données (DG) et mesures (MG) relatives au gisement.Specialists are familiar with many software tools for constructing a reservoir model, using data (DG) and measurements (MG) for the deposit.
La
Pour prendre en compte le rôle du réseau de fractures dans la simulation des écoulements au sein du gisement, il est nécessaire d'associer à chacun de ces volumes élémentaires (mailles du modèle de réservoir) une modélisation des fractures.To take into account the role of the fracture network in the simulation of the flows within the deposit, it is necessary to associate with each of these elementary volumes (meshes of the reservoir model) a fracture modeling.
Ainsi, un objet de l'invention concerne une méthode pour construire une représentation d'un gisement de fluide traversé par un réseau de fractures et par au moins un puits. Cette méthode comporte la discrétisation du gisement en un ensemble de mailles (étape 1 précédemment décrite). Puis la méthode comporte les étapes suivantes :
- a. on caractérise les fractures par des paramètres statistiques (PSF) à partir d'observations du gisement,
- b. on déduit de ces paramètres statistiques (PSF) un tenseur de perméabilité équivalente et une ouverture moyenne des fractures, à partir desquels on construit une image représentative du réseau de fractures délimitant des blocs poreux et des fractures ;
- c. on déduit de ce tenseur une direction d'écoulement de fluide autour du puits ;
- d. on définit autour du puits, une première zone de frontière elliptique centrée sur le puits et contenant ce puits, au moins une seconde zone de frontière elliptique centrée sur le puits et de limite intérieure confondue avec la frontière elliptique de la première zone, les zones étant orientées selon la direction d'écoulement de fluide ;
- e. on simplifie l'image représentative du réseau de fractures dans chaque maille appartenant à au moins une zone ;
- f. on réitère à l'étape b) en modifiant les paramètres statistiques (PSF) de façon à minimiser la différence entre le résultat de test de puits et le résultat de simulation de test de puits à partir de l'image simplifiée ; et
- g. on associe à chacune des mailles au moins une valeur de perméabilité équivalente et une valeur d'ouverture moyenne des fractures, ces valeurs étant déterminées à partir des paramètres statistiques (PSF) modifiés.
- at. Fractures are characterized by statistical parameters (PSF) from observations of the deposit,
- b. from these statistical parameters (PSF) is deduced an equivalent permeability tensor and an average fracture opening, from which a representative image of the fracture network delimiting porous blocks and fractures is constructed;
- vs. a tension direction of fluid around the well is deduced from this tensor;
- d. a first elliptical boundary zone centered around the well and containing the well, at least a second elliptical boundary zone centered on the well and an inner boundary coinciding with the elliptical boundary of the first zone, the zones being defined around the well; oriented in the direction of fluid flow;
- e. the representative image of the fracture network is simplified in each mesh belonging to at least one zone;
- f. reiterating in step b) modifying the statistical parameters (PSF) so as to minimize the difference between the well test result and the well test simulation result from the simplified image; and
- boy Wut. each mesh is associated with at least one equivalent permeability value and an average opening value of the fractures, these values being determined from the modified statistical parameters (PSF).
Ces étapes sont détaillées ci-après.These steps are detailed below.
Les spécialistes en charge de la caractérisation statistique du réservoir, réalisent des observations (OF) directes et indirectes du réservoir. Pour cela, ils disposent 1) de carottes de puits extraites du réservoir sur lesquelles une étude statistique des fractures intersectées est effectuée, 2) d'affleurements caractéristiques du réservoir qui présentent l'avantage de fournir une vision à grande échelle du réseau de fractures 3) d'images sismiques leur permettant d'identifier de grands événements géologiques.The specialists in charge of the statistical characterization of the reservoir, perform direct and indirect observations (OF) of the reservoir. For this, they have 1) well cores extracted from the reservoir on which a statistical study of the intersected fractures is carried out, 2) characteristic outcrops of the reservoir which have the advantage of providing a large-scale vision of the fracture network 3 ) seismic images allowing them to identify large geological events.
Ces mesures permettent de caractériser les fissures par des paramètres statistiques (PSF) : leur densité, leur longueur, leur orientation, leur inclinaison et leur ouverture respectives, et bien sur, leur distribution au sein du réservoir.These measurements make it possible to characterize the cracks by statistical parameters (PSF): their density, their length, their orientation, their inclination and their respective opening, and of course, their distribution within the tank.
A l'issue de cette étape de caractérisation des fractures, nous disposons de paramètres statistiques (PSF) décrivant les réseaux de fractures à partir desquels des images réalistes des réseaux réels (géologiques) peuvent être reconstruites (générées) à l'échelle de chacune des mailles (cellules) du modèle de réservoir considéré (domaine de simulation).At the end of this stage of fracture characterization, we have statistical parameters (PSF) describing the fracture networks from which realistic images of the real (geological) networks can be reconstructed (generated) at the scale of each of the meshes (cells) of the reservoir model considered (simulation domain).
L'objectif de la caractérisation et de la modélisation du réseau de fractures de gisement, est de proposer un modèle de fractures validé sur les écoulements locaux autour des puits. Ce modèle de fractures est alors étendu à l'échelle du réservoir afin de réaliser des simulations de production. Pour ce faire, on associe des propriétés d'écoulement à chaque maille du modèle de réservoir (MR) (tenseur de perméabilité, porosité) des deux milieux (fracture et matrice).The objective of the characterization and modeling of the deposit fracture network is to propose a model of fractures validated on the local flows around the wells. This fracture model is then extended to the reservoir scale in order to perform production simulations. To do this, flow properties are associated with each mesh of the reservoir model (MR) (permeability tensor, porosity) of the two media (fracture and matrix).
Ces propriétés peuvent être déterminées, soit directement à partir des paramètres statistiques (PSF) décrivant les réseaux de fractures, soit à partir d'un réseau de fracture discret (DFN) obtenus à partir des paramètres statistiques (PSF).These properties can be determined either directly from the statistical parameters (PSF) describing the fracture networks, or from a discrete fracture network (DFN) obtained from the statistical parameters (PSF).
A partir d'un modèle de réservoir du gisement étudié, on lui associe en chaque maille une représentation détaillée (DFN) de la complexité interne du réseau de fractures aussi fidèle que possible des observations directes et indirectes du réservoir. La
Pour construire un réseau de fracture discret en chaque maille d'un modèle de réservoir, on peut utiliser des logiciels de modélisation, bien connus des spécialistes, tel que le logiciel FRACAFlow® (IFP, France). Ces logiciels utilisent les paramètres statistiques déterminés à l'étape de caractérisation des fractures.To build a discrete fracture network in each cell of a reservoir model, we can use modeling software, known in the art, such as FracaFlow ® software (IFP, France). These programs use the statistical parameters determined in the fracture characterization step.
L'étape suivante consiste à déterminer les propriétés d'écoulement des fractures initiales (Cf, e), puis à calibrer ces propriétés au moyen de simulation de tests de puits sur des modèles d'écoulement locaux discrets, hérités de l'image réalistes du réseau réel (géologique) de fractures à l'échelle du réservoir. Bien que ne couvrant qu'une aire limitée (aire de drainage) autour du puits, un tel modèle de simulation de tests de puits comporte encore de très nombreux noeuds de calcul si le réseau de fractures est dense. En conséquence, la taille des systèmes à résoudre et/ou la durée des calculs demeurent souvent prohibitives. D'où la nécessité de recourir à une procédure de simplification du réseau de fracture.The next step is to determine the flow properties of the initial fractures (C f , e) and then calibrate these properties by means of well-test simulations on discrete, realistic image-inherited local flow models. the actual (geological) fracture network at the reservoir scale. Although only covering a limited area (drainage area) around the well, such a well testing simulation model still has a large number of computational nodes if the fracture network is dense. As a result, the size of the systems to be solved and / or the duration of the calculations often remain prohibitive. Hence the need for a procedure to simplify the fracture network.
En raison de son extrême complexité géométrique, le réseau de fractures obtenu à l'étape précédente et représentatif du réservoir fracturé réel, ne peut être utilisé pour simuler, i.e. reproduire et/ou prédire, les écoulements locaux autour du puits.Because of its extreme geometric complexity, the fracture network obtained in the previous step and representative of the actual fractured reservoir can not be used to simulate, i.e. reproduce and / or predict, the local flows around the well.
Pour contourner cet obstacle, la méthode selon l'invention, utilise un procédé basé sur la division du domaine de simulation (c'est-à-dire le modèle de réservoir) en au moins trois types de zones autour de chaque puits (
- une première zone, dans laquelle aucune simplification du réseau de fracture n'est réalisée. Cette zone contient le puits et son centre. Elle est notée ZNS, pour "Zone Non Simplifiée".
- une seconde zone, au contact de la première zone, dans laquelle une simplification modérée du réseau de fracture est réalisée. Cette zone est notée ZP, pour "Zone à simplifier Proche du puits".
- une troisième zone, au contact de la seconde zone, dans laquelle une simplification importante du réseau de fracture est réalisée. Cette zone est notée ZL, pour "Zone Lointaine du puits".
- a first zone, in which no simplification of the fracture network is performed. This area contains the well and its center. It is noted ZNS, for "Not Simplified Zone".
- a second zone, in contact with the first zone, in which a moderate simplification of the fracture network is performed. This zone is marked ZP, for "Zone to be simplified Near the well".
- a third zone, in contact with the second zone, in which a major simplification of the fracture network is performed. This zone is marked ZL, for "Far Zone of the Well".
L'invention n'est pas limitée à la définition de trois zones. On peut également diviser le domaine en n zones, la simplification du réseau étant croissante de la zone1 (ZNS) à la zone n (la plus éloignée du puits). On peut ainsi créer une zone ZNS, n1 zones de type ZP, et n2 zones de type ZL.The invention is not limited to the definition of three zones. It is also possible to divide the domain into n zones, the simplification of the network being increased from the zone 1 (ZNS) to the zone n (the furthest away from the well). It is thus possible to create a zone ZNS, n1 zones of type ZP, and n2 zones of type ZL.
La modélisation du réseau de fractures a pour but de simuler les réponses en écoulement de puits (tests d'écoulement transitoires ou pseudo-permanents, interférences, débitmétrie, etc.). Il s'agit de simuler par exemple la production d'huile via chaque puits foré à travers le réservoir.The purpose of fracture network modeling is to simulate well-flow responses (transient or pseudo-permanent flow tests, interference, flowmetry, etc.). This involves simulating, for example, the production of oil via each well drilled through the reservoir.
Pour chaque puits, chaque zone est définie selon une limite extérieure formant une ellipse centrée sur le puits. Les trois zones sont donc concentriques et de frontières elliptiques. Les deux zones simplifiées ZP et ZL ont des limites intérieures correspondant à la limite extérieure des zones respectives ZNS et ZP. A l'exception de la zone non simplifiée, les zones sont donc des couronnes elliptiques centrées sur le puits. Pour construire chaque zone, il convient donc de définir :
- l'orientation de l'ellipse, c'est-à-dire la direction du grand axe de l'ellipse (perpendiculaire à la direction du petit axe) ;
- les dimensions de l'ellipse, c'est-à-dire la longueur des axes.
- the orientation of the ellipse, that is to say the direction of the major axis of the ellipse (perpendicular to the direction of the minor axis);
- the dimensions of the ellipse, that is to say the length of the axes.
Leur orientation est déterminée par les directions d'écoulement déduites d'un calcul de perméabilités équivalentes sur la zone ZNS. Ce type de calcul des perméabilités équivalentes est bien connu des spécialistes. On peut par exemple utiliser la méthode numérique de calcul de propriétés équivalentes de milieux fracturés, implémentée dans le logiciel FracaFlow (IFP Energies nouvelles, France) et rappelé ci-après.Their orientation is determined by the flow directions deduced from a calculation of equivalent permeabilities on the ZNS zone. This type of calculation of equivalent permeabilities is well known to those skilled in the art. One can for example use the numerical method for calculating equivalent properties of fractured media, implemented in the FracaFlow software (IFP Energies nouvelles, France) and recalled hereinafter.
Selon cette méthode, un tenseur de perméabilité, représentatif des propriétés d'écoulement du réseau de fractures discrétisé (DFN) peut être obtenu via deux méthodes de mise à l'échelle (méthodes dites d'upscaling).
- La première méthode, analytique, dite "local analytical upscaling", est basée sur une approche analytique décrite dans les documents suivants :
-
Oda M. (1985): Permeability tensor for discontinuous Rock Masses, Geotechnique Vol 35, 483-495 - demande de
brevet EP 2 037 080
-
- The first method, known as "local analytical upscaling", is based on an analytical approach described in the following documents:
-
Oda M. (1985): Permeability tensor for discontinuous Rock Masses, Geotechnical Vol 35, 483-495 - request
EP Patent 2,037,080
-
Elle présente l'avantage d'être très rapide. Son domaine d'application se limite, toutefois, à des réseaux de fractures bien connectés. Dans le cas contraire, des erreurs importantes sur le tenseur de perméabilité peuvent être constatées.
- La seconde méthode, numérique, dite "local numerical upscaling", est décrite dans les documents suivants :
-
Bourbiaux, B., et al., 1998, "A Rapid and Efficient Methodology to Convert Fractured Reservoir Images into a Dual-Porosity Model", Oil & Gas Science and Technology, Vol. 53, No. 6, Nov.-Déc. 1998, 785-799 -
brevet FR 2.757.947 US 6.023.656 FR 2.757.957 US 6.064.944
-
- The second method, digital, called "local numerical upscaling", is described in the following documents:
-
Bourbiaux, B., et al., 1998, "A Rapid and Efficient Methodology to Convert Fractured Reservoir Images into a Dual-Porosity Model," Oil & Gas Science and Technology, Vol. 53, No. 6, Nov.-Dec. 1998, 785-799 -
FR 2,757,947 US 6,023,656 FR 2 757 957 US 6,064,944
-
Elle est basée sur la résolution numérique des équations d'écoulement sur un maillage discret du réseau de fractures pour différentes conditions aux limites du bloc de calcul considéré. Le tenseur de perméabilités équivalentes est obtenu par identification des rapports entre débit et perte de charge aux limites du bloc de calcul. Cette approche, plus coûteuse que la précédente, présente l'avantage de bien caractériser un réseau donné (même peu connecté).It is based on the numerical resolution of the flow equations on a discrete mesh of the fracture network for different boundary conditions of the considered block of computation. The tensor of equivalent permeabilities is obtained by identifying the ratios between flow and pressure drop at the limits of the calculation block. This approach, more expensive than the previous one, has the advantage of well characterizing a given network (even if it is not connected).
Selon un mode de réalisation, on peut choisir l'une ou l'autre des deux méthodes précédentes de manière à optimiser précision et rapidité des calculs, en appliquant la méthode décrite dans la demande de
Cette technique permet, à titre préliminaire, de déterminer les tenseurs de perméabilité de quelques mailles du modèle de réservoir entourant le puits, et considérées représentatives de l'écoulement de la ZNS. La diagonalisation de ces tenseurs de perméabilité fournit les vecteurs propres orientés suivant les directions principales d'écoulement recherchées. On est alors en mesure d'orienter le domaine elliptique qualifié de ZNS et centré sur le puits, suivant le demi-grand axe de l'ellipse de perméabilité déterminé au moyen de ces calculs préliminaires.This technique makes it possible, as a preliminary step, to determine the permeability tensors of a few cells of the reservoir model surrounding the well, and considered representative of the flow of the ZNS. The diagonalization of these permeability tensors provides the eigenvectors oriented along the main directions of flow sought. We are then able to orient the qualified elliptical domain of ZNS and centered on the well, along the semi-major axis of the ellipse of permeability determined by means of these preliminary calculations.
Puis, pour définir les dimensions de l'ellipse, on définit un rapport de forme de cette ellipse ainsi qu'une distance entre ellipses concentriques (distance séparant les limites elliptiques intérieure et extérieure d'une couronne concentrique donnée) :
- rapport de forme des ellipses. En notant Lmax la demi longueur du grand axe de l'ellipse, et Lmin la demi longueur du petit axe de l'ellipse, ce rapport de forme est défini par Lmax /Lmin. II convient de ne pas le fixer de manière arbitraire mais au contraire de le choisir de manière à reproduire l'anisotropie d'écoulement. On peut utiliser le tenseur de perméabilité équivalente calculé pour déterminer l'orientation. Si l'on note Kmin et Kmax les valeurs principales de ce tenseur, alors les ellipses sont orientées suivant les directions principales de perméabilité avec un rapport de forme LmaxlLmin égal à la racine carrée du rapport Kmax/Kmin :
- distance entre ellipses concentriques (frontières entre zones) : un critère de dimensionnement conforme à une précision de modélisation uniformément distribuée sur le domaine de simulation consiste à fixer les demi-longueurs du grand axe Lmax (ou du petit axe Lmin ) de 2 ellipses successives i+1 et i à des valeurs en progression géométrique de raison r (égale à 2 par exemple) constante, soit :
- shape ratio of ellipses. Noting L max the half length of the major axis of the ellipse, and L min the half length of the minor axis of the ellipse, this aspect ratio is defined by L max / L min . It should not be fixed arbitrarily but rather to choose it so as to reproduce the flow anisotropy. The calculated equivalent permeability tensor can be used to determine the orientation. If we denote K min and K max the principal values of this tensor, then the ellipses are oriented according to the principal directions of permeability with a form ratio L max lL min equal to the square root of the ratio K max / K min :
- distance between concentric ellipses (boundaries between zones): a dimensioning criterion conforming to a uniformly distributed modeling precision on the simulation domain consists in setting the half-lengths of the major axis L max (or the minor axis L min ) of 2 ellipses successive i + 1 and i to values in geometric progression of reason r (equal to 2 for example) constant, that is:
Une fois ce dimensionnement global effectué, les procédures de délimitation et de simplification des sous-zones sont mises en oeuvre en s'appuyant également sur les méthodes de calcul des propriétés équivalentes des milieux fracturés.Once this overall dimensioning has been carried out, the delimitation and simplification procedures of the sub-areas are implemented, also relying on the methods for calculating the equivalent properties of the fractured media.
Chacune des zones, sauf la zone ZNS, est ensuite divisée en sous-zones (ssZP, ssZL1, ssZL2), à l'intérieur desquelles une simplification du réseau de fracture est effectuée. Le nombre de sous-zones par zone dépend du type de la zone (ZNS, ZP ou ZL) et de l'hétérogénéité de cette dernière. Ainsi :
- La zone ZNS est une zone à conserver intacte (non simplifiée); aucune sous-zone n'y est créée.
- Les zones à simplifier les plus proches d'un puits (zones de type ZP), nécessitent une attention particulière. Pour modéliser correctement les variations locales des propriétés d'écoulement dans cette zone, les ZP sont divisées en un nombre de sous-zones égal au nombre de blocs de mailles présents dans les zones ZP. Le terme "bloc de mailles" est employé pour désigner un "empilement" de mailles verticales (de type CPG (Corner Point Grid) par exemple) du modèle de réservoir, délimitées par les mêmes poteaux droits sub verticaux. Il ne s'agit donc pas des blocs équivalents.
- Les zones à simplifier les plus éloignées d'un puits (zone de type ZL), sont des couronnes elliptiques concentriques qui couvrent des surfaces de plus en plus importantes au fur et à mesure qu'on s'éloigne du puits. Étant plus éloignées des puits, il est acceptable d'être moins précis dans la détection des hétérogénéités que dans le cas des zones de type ZP. L'hétérogénéité des propriétés d'écoulement d'une zone ZL reste toutefois le facteur principal contrôlant le découpage en sous-zones. Pour échantillonner une zone ZL, un balayage angulaire, centré sur le puits, est effectué sur la limite elliptique extérieure séparant la zone ZL considérée de sa voisine extérieure. Tous les degrés, un bloc de mailles est sélectionné. Ainsi, la zone ZL est caractérisée par au plus 360 blocs. Pour chacun de ces blocs, un calcul du tenseur de perméabilité de fracture équivalente est effectué. Ce tenseur permet de caractériser les propriétés dynamiques du réseau de fractures du bloc étudié. Les valeurs de perméabilité et l'orientation de l'écoulement obtenues sont alors comparées entre blocs voisins. Si les propriétés de deux blocs voisins sont proches, ces deux blocs seront considérés comme appartenant à une même sous-zone. Dans le cas contraire (différence de 20 % sur les valeurs de perméabilité principale ou de 10 degrés sur les directions principales de perméabilité par exemple), les deux blocs sont affectés à des sous-zones différentes ainsi définies. De cette manière, la limite extérieure de la zone ZL considérée est divisée en arcs de forme elliptique (
figure 1 ). Les sous-zones sont alors obtenues en reliant les points-extrémités de chacun de ces arcs au centre (puits) de l'ellipse (lignes en pointillées de lafigure 5 ) : chaque sous-zone est ainsi définie comme l'aire comprise entre le cloisonnement radial intra-zone (traits pointillés) et les limites elliptiques inter-zones (en trait plein).
- The ZNS zone is an area to be kept intact (not simplified); no subzone is created there.
- Areas to be simplified closer to a well (ZP type zones) require special attention. To correctly model the local variations of the flow properties in this area, the PZs are divided into a number of subfields equal to the number of mesh blocks present in the ZP zones. The term "block of meshes" is used to designate a "stack" of vertical meshes (of the CPG type (Corner Point Grid) for example) of the reservoir model, delimited by the same right sub vertical posts. It is not therefore equivalent blocks.
- The zones to be simplified furthest from a well (ZL-type zone) are concentric elliptical crowns that cover increasingly large surfaces as one moves away from the well. Being further away from wells, it is acceptable to be less accurate in detecting heterogeneities than in the case of ZP-type zones. The heterogeneity of the flow properties of a ZL zone, however, remains the main factor controlling the subzone cutting. To sample a zone ZL, an angular sweep, centered on the well, is performed on the outer elliptical boundary separating the zone ZL considered from its outer neighbor. Every degree, a block of meshes is selected. Thus, the ZL area is characterized by at most 360 blocks. For each of these blocks, a calculation of the equivalent fracture permeability tensor is performed. This tensor makes it possible to characterize the dynamic properties of the fracture network of the studied block. The permeability values and the flow orientation obtained are then compared between neighboring blocks. If the properties of two neighboring blocks are close, these two blocks will be considered as belonging to the same sub-zone. In the opposite case (difference of 20% on the principal permeability values or 10 degrees on the principal directions of permeability for example), the two blocks are assigned to different sub-zones thus defined. In this way, the outer limit of the zone ZL considered is divided into elliptic arcs (
figure 1 ). The subfields are then obtained by connecting the endpoints of each of these arcs to the center (well) of the ellipse (dashed lines of thefigure 5 ): each subfield is thus defined as the area between the intra-zone radial partitioning (dashed lines) and the inter-zone elliptical boundaries (solid line).
Comme indiqué précédemment, la délimitation des sous-zones (points-limites des arcs de frontière elliptiques inter-zones) repose sur la comparaison des perméabilités équivalentes calculées sur les "blocs" voisins matérialisant ces arcs. La méthode analytique de calcul des perméabilités équivalentes est utilisée de préférence, car il s'agit, dans ce cas, de déterminer si un bloc a le même comportement dynamique que le bloc voisin. Bien que, pour un réseau faiblement connecté, l'approche analytique fournisse des résultats erronés, les erreurs sont systématiques, similaires d'un bloc à l'autre, ce qui autorise la comparaison des résultats entre blocs qui, précisons-le, n'exige pas une grande précision compte tenu du simple objectif de définition de zone. Ainsi l'approche analytique est tout à fait justifiée, avec l'avantage considérable de permettre des calculs beaucoup plus rapides, garants d'une faisabilité pratique.As indicated above, the delimitation of the sub-zones (boundary points of the inter-zone elliptic boundary arcs) is based on the comparison of the equivalent permeabilities calculated on the neighboring "blocks" materializing these arcs. The analytical method for calculating equivalent permeabilities is preferably used because it is in this case to determine if a block has the same dynamic behavior as the neighboring block. Although, for a weakly connected network, the analytical approach provides erroneous results, the errors are systematic, similar from one block to the other, which allows the comparison of the results between blocks which, let us specify it, does not require a great precision considering the simple objective of definition of zone. Thus the analytical approach is quite justified, with the considerable advantage of allowing much faster calculations, guaranteeing a practical feasibility.
Une fois le découpage en sous zones effectué, les calculs de mise à l'échelle vont permettre de remplacer le réseau de fractures de ces sous-zones par un réseau simplifié possédant les mêmes propriétés d'écoulement que le réseau originel. Dans ce cas, et contrairement à ce qui précède, le calcul du tenseur de perméabilité de fracture équivalente doit être le plus précis possible.Once the sub-zones have been cut, the scaling calculations will allow the fracture network of these sub-zones to be replaced by a simplified network having the same flow properties as the original network. In this case, and contrary to the above, the calculation of the equivalent fracture permeability tensor must be as accurate as possible.
Afin de profiter au mieux des avantages des deux méthodes de mise à l'échelle, le tenseur de perméabilité est déterminé par l'une ou l'autre de ces deux méthodes par exemple, suivant la procédure de sélection décrite dans le document
- Réseau bien connecté : pour ce cas caractérisé par un indice de connectivité voisin ou dépassant 3 (au moins 3 intersections par fracture du réseau en moyenne), la méthode d'upscaling analytique est choisie car sa précision est garantie compte tenu de la bonne connectivité du réseau, avec l'avantage supplémentaire essentiel de la rapidité.
- Réseau peu/mal connecté : dans ce cas, l'indice de connectivité est compris entre 1 et 3 (ce qui correspond à un nombre d'intersections de fractures compris entre une et trois fois le nombre de fractures) et la méthode d'upscaling numérique est utilisée pour calculer de manière fiable le tenseur de perméabilité.
- Lorsque le réseau est très peu voire non connecté (nombre d'intersections voisin ou inférieur au nombre de fractures), les fractures originelles (peu nombreuses) sont conservées, c'est-à-dire que la sous-zone en question n'est pas simplifiée.
- Well connected network: for this case characterized by a connectivity index close to or greater than 3 (at least 3 intersections per network fracture on average), the analytical upscaling method is chosen because its accuracy is guaranteed given the good connectivity of the network. network, with the essential additional benefit of speed.
- Poorly / poorly connected network: in this case, the connectivity index is between 1 and 3 (which corresponds to a number of fracture intersections between one and three times the number of fractures) and the upscaling method digital is used to reliably calculate the permeability tensor.
- When the network is very little or not connected (number of intersections close to or less than the number of fractures), the original fractures (few) are preserved, ie the sub-zone in question is not not simplified.
Une fois ces calculs de perméabilité équivalente effectués pour chaque sous-zone, les autres paramètres d'écoulement équivalents caractérisant les sous-zones simplifiées sont facilement déterminés suivant les méthodes et équivalences suivantes.
- calcul d'un premier réseau équivalent (réseau parallélépipédique dit de Warren et Root) dit "fin", issu directement de la méthode du
brevet FR 2.757.957 US 6.064.944 Figure 4 ) est caractérisé par des espacements des fractures s1 fin, s2 fin dans les 2 directions orthogonales de perméabilité principale 1 et 2 définissant le réseau de fractures.
Un paramètre supplémentaire, l'ouverture des fractures (efin), caractérise les fractures de ce réseau fin. La valeur des ouvertures de fracture est en pratique quasiment toujours négligeable devant l'espacement entre fractures : il est tenu compte de cette hypothèse dans les formules qui suivent, où on suppose également la même valeur d'ouverture pour les 2 familles de fractures. Compte tenu de l'égalité des porosités du réseau DFN initial (φf ) et du réseau équivalent fin, efin se déduit du volume de fractures du réseau DFN initial, Vf init, et du volume total de roche VT comme suit :
Pour mémoire, les perméabilités principales de fracture équivalentes, issues des calculs décrits plus haut, valent k1= keq Max et k2= keq Min suivant les directions principales 1 et 2. On déduit les conductivités des fractures du réseau équivalent fin, Cf1 fin et Cf2 fin, suivant ces deux directions d'écoulement, en écrivant la conservation des flux par unité d'aire de milieu fracturé :
Enfin, le milieu matriciel entre fractures possède une perméabilité km fin. - remplacement de ce réseau fin (
Figure 4 ) par le réseau équivalent dit "grossier" (Figure 6 ), car comportant des fractures plus espacées afin d'accroître le degré de simplification en vue des simulations d'écoulement ultérieures. Les propriétés géométriques et d'écoulement de ce réseau grossier sont les suivantes:- ∘ des espacements de fractures s1 gros et s2 gros tels que :
- G<Gmax-zone où Gmax-zone est tel que
- Max(s1 gros,s2 gros) < (dimension latérale minimale de la sous-zone )/6
- DLM dimension latérale minimale de la sous-zone donnée ;
- s1 fin, s2 fin : les espacements de fractures dans la représentation dite de Warren et Root.
- ∘ des conductivités de fractures Cf1 gros et Cf2 gros, dont les valeurs permettent de conserver les flux par unité d'aire de milieu fracturé, c'est-à-dire également les perméabilités équivalentes, soit : Cf1 gros=s2 gros .k1 et Cf2 gros=s1 gros.k2
ou encore, sachant que Cf1 fin=s2 fin .k1 et Cf2 fin=S1 fin.k2, - ∘ une ouverture de fracture egros permettant de conserver à nouveau la porosité de fracture φf du réseau initial (égale à celle du réseau équivalent grossier) :
- ∘ une perméabilité de matrice km gros conservant la valeur du paramètre d'échange matrice-fissure :
- ∘ des espacements de fractures s1 gros et s2 gros tels que :
soit
- calculation of a first equivalent network (parallelepipedal network called Warren and Root) said "end", directly from the method of
FR 2,757,957 US 6,064,944 Figure 4 ) is characterized by spacings of fractures s 1 end , s 2 end in the 2 orthogonal directions ofmain permeability 1 and 2 defining the network of fractures.
An additional parameter, the opening of fractures (e end ), characterizes the fractures of this fine network. In practice, the value of the fracture openings is practically always negligible compared to the spacing between fractures: this assumption is taken into account in the following formulas, where the same opening value is also assumed for the two families of fractures. Given the equal porosity of the initial DFN network ( φ f ) and the fine equivalent network, e end is deduced from the fracture volume of the initial DFN network, V f init , and the total volume of rock V T as follows:
For memory, the main permeabilities equivalent fracture, resulting from the calculations described above, are worth k 1 = k eq Max and k 2 = k eq Min followingprincipal directions 1 and 2. We deduce the conductivities of the equivalent network end fractures, C f1 end and C f2 end , following these two directions of flow, by writing the conservation of flows per unit area of fractured medium:
Finally, the matrix medium between fractures has a permeability k m fin . - replacing this fine network (
Figure 4 ) by the so-called "coarse" equivalent network (Figure 6 ), because with more spaced fractures to increase the degree of simplification for subsequent flow simulations. The geometric and flow properties of this coarse network are as follows:- ∘ fracture spacings of 1 large and 2 large such as:
- G <G max-zone where G max-zone is such that
- Max (s 1 large , s 2 large ) <(minimal lateral dimension of the sub-area) / 6
- DLM minimum lateral dimension of the given subfield;
- s 1 end , s 2 fin : fracture spacings in the representation of Warren and Root.
- ∘ fracture conductivities C f1 and C large large f2, whose values help maintain flows through the middle of fractured unit area, that is to say also the equivalent permeabilities, ie: C f1 = s big 2 big .k 1 and C f2 big = s 1 big .k 2
or again, knowing that C f1 end = s 2 end .k 1 and C f2 fi n = S 1 end .k 2 , - ∘ a large fracture opening e allowing retaining the fracture porosity φ f of the initial network (equal to that of the coarse equivalent network):
- ∘ a matrix permeability k m large retaining the value of the matrix-crack exchange parameter:
- ∘ fracture spacings of 1 large and 2 large such as:
is
Enfin, une fois, cette opération d'équivalence réseau fin --> réseau grossier effectuée pour chaque sous-zone (ou "bloc" de la ZP), les fractures des réseaux simplifiés obtenus sont prolongées hors des limites de sous-zones (ou "blocs") afin de garantir un "recouvrement" partiel suffisant et donc une connectivité horizontale suffisante des réseaux simplifiés de sous-zones (ou "blocs" de ZP) voisines. A cette fin, et suivant une procédure éprouvée par des tests, les fractures du réseau simplifié peuvent ainsi être prolongées d'une longueur égale à 60% de l'espacement maximum (s1 gros,s2 gros) des fractures de ce réseau.Finally, once this fine-network network equivalence operation performed for each sub-zone (or "block" of the ZP), the fractures of the simplified networks obtained are extended beyond the limits of sub-zones (or "blocks") to ensure sufficient partial "overlap" and thus sufficient horizontal connectivity of the simplified sub-area networks (or "blocks" of ZP) nearby. For this purpose, and according to a procedure tested by tests, the fractures of the simplified network can thus be extended by a length equal to 60% of the maximum spacing (s 1 large , s 2 large ) fractures of this network.
Une telle image est représentée sur la
L'étape suivante est la calibration des propriétés d'écoulement des fractures (conductivité et ouverture des fractures), localement autour des puits. Celle-ci nécessite la simulation de tests de puits. Selon l'invention cette simulation de tests de puits est réalisée sur les modèles d'écoulement simplifiés (
Ce type de calibration est bien connu des spécialistes. On peut par exemple utiliser la méthode décrite dans le
Les résultats de ces simulations permettent de calibrer (estimer) la géométrie et les propriétés d'écoulement des fractures, telles que les conductivités des réseaux de fractures du réservoir étudié et les ouvertures.The results of these simulations make it possible to calibrate (estimate) the geometry and the fracture flow properties, such as the conductivities of the fracture networks of the studied reservoir and the openings.
A ce stade, l'ingénieur de réservoir dispose alors des données requises pour construire le modèle d'écoulement à l'échelle du réservoir. En effet, les simulations de réservoir fracturé adoptent souvent l'approche "double porosité", proposée par exemple par
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 champ considéré. La définition d'un scénario optimal d'injection d'eau consistera, 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 des fractures sur la progression des fluides au sein du réservoir.The reservoir engineer chooses a production process, for example the water injection recovery process, the optimal implementation scenario for the field in question. The definition of an optimal water injection scenario will consist, for example, in determining the number and location (position and spacing) of the injection and production wells in order to best take into account the impact of the fractures on the water. progression of fluids within the reservoir.
En fonction du scénario choisi, de la représentation double milieu du gisement, et de la formule reliant le flux d'échange de masse et/ou d'énergie à la différence de potentiel matrice-fracture, on est alors capable de simuler la production d'hydrocarbures escomptée, au moyen du simulateur d'écoulement (logiciel) dit à double milieu.Depending on the chosen scenario, the double-middle representation of the deposit, and the formula linking the mass and / or energy exchange flux to the matrix-fracture potential difference, one is then able to simulate the production of hydrocarbons, by means of the so-called dual medium flow simulator (software).
A tout instant t de la production simulée, à partir des données d'entrées E(t) (données fixes ou variables au fur et à mesure du temps simulé), et de la formule reliant le flux (f) d'échange à la différence de potentiel (ΔΦ), le simulateur résout l'ensemble des équations propres à chaque maille et chacune des deux grilles du modèle (équations faisant intervenir la formule d'échange matrice-fissure décrite plus haut) et délivre ainsi les valeurs solutions des inconnues S(t) (saturations, pressions, concentrations, température, ...) à cet instant t. De cette résolution, découle la connaissance des quantités d'huile produites et de l'état du gisement (distribution des pressions, saturations, etc...) à l'instant considéré.At any instant t of the simulated production, from the input data E (t) (fixed or variable data as the simulated time), and from the formula connecting the flow (f) of exchange to the potential difference (ΔΦ), the simulator solves all the equations specific to each mesh and each of the two grids of the model (equations involving the matrix-crack exchange formula described above) and thus delivers the solution values of the unknowns S (t) (saturations, pressures, concentrations, temperature, etc.) at this instant t. From this resolution, comes the knowledge of the quantities of oil produced and the state of the deposit (distribution of pressures, saturations, etc ...) at the moment considered.
En sélectionnant 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 permettant d'optimiser la production du gisement fracturé considéré suivant les critères technico-économiques sélectionnés.By selecting various scenarios, characterized for example by various respective locations of the injectors and producers, and by simulating the production of hydrocarbons for each of them according to step 3, one can select the scenario to optimize the production of the deposit. fractured considered according to the selected technical and economic criteria.
Les spécialistes exploitent alors le gisement selon ce scénario permettant d'optimiser la production du gisement.The specialists then exploit the deposit according to this scenario making it possible to optimize the production of the deposit.
Claims (11)
- A method for optimizing the development of a fluid reservoir traversed by a fracture network and by at least one well, wherein a representation of said fluid reservoir is constructed, said reservoir is discretized into a set of cells and the fractures are characterized by statistical parameters (PSF) from observations of said reservoir, characterized in that the method comprises the following stages:a) deducing from said statistical parameters (PSF) an equivalent permeability tensor and an average opening for said fractures, from which an image representative of the fracture network delimiting porous blocks and fractures is constructed,b) deducing from said tensor a direction of fluid flow around said well,c) defining around said well a first elliptical boundary zone centered on said well and containing said well, and at least a second elliptical boundary zone centered on said well and forming an elliptical ring with the elliptical boundary of said first zone, said zones being oriented in said direction of fluid flow,d) simplifying said image representative of the fracture network in a different manner in each one of said zones,e) using said simplified image to construct said representation of the fluid reservoir,f) using said representation of the fluid reservoir and a flow simulator to optimize the development of said fluid reservoir.
- A method as claimed in claim 1, wherein statistical parameters (PSF) are selected from among the following parameters: fracture density, fracture length, fracture orientation, fracture inclination, fracture opening and fracture distribution within the reservoir.
- A method as claimed in any one of the previous claims, wherein an aspect ratio is determined for each zone, defined from the lengths of the axes of the ellipse making up the boundary of said zone, so as to reproduce a flow anisotropy around said well, and the zones are constructed so as to respect said aspect ratio.
- A method as claimed in claim 3, wherein said aspect ratio is determined by means of the principal values of said permeability tensor.
- A method as claimed in any one of the previous claims, wherein a distance is defined between said boundaries between zones, so as to give an equal weight to each zone in terms of pressure difference recorded in each zone under permanent flow regime conditions.
- A method as claimed in claim 5, wherein said distance is defined by setting the lengths of one of the two axes of two successive ellipses at values in geometric progression of constant ratio.
- A method as claimed in any one of the previous claims, wherein three zones are constructed, a first zone (ZNS) containing the well, within which no simplification of said image is performed, a second zone (ZP) in contact with said first zone, within which a first simplification of said image is performed, and a third zone (ZL) in contact with said second zone, within which a second simplification of said image is performed, said second simplification being more significant than said first simplification.
- A method as claimed in claim 7, wherein said second and third zones are divided into sub-zones by applying the following stages:- said second zone is divided into a number of sub-zones equal to a number of blocks of cells present in said zone, a block of cells designating a vertical pile of cells,- said third zone is divided by carrying out the following stages:∘ dividing every degree the boundary of said third zone, defining 360 arcs,∘ defining a sub-zone by connecting end points of each of said arcs to the centre of the ellipse forming the boundary,∘ for each one of said sub-zones, calculating an equivalent fracture permeability tensor from which an orientation of the flows in said sub-zone is determined,∘ comparing the equivalent fracture permeability values and the flow orientation between neighbouring sub-zones, and∘ grouping neighbouring sub-zones together into a single sub-zone when a difference between the permeability values is below a first threshold and when a difference between the flow orientations is below a second threshold.
- A method as claimed in any one of the previous claims, wherein said image is simplified by carrying out the following stages:- constructing a fracture network (RFE) equivalent to said image, by means of a so-called Warren and Root representation, wherein the network is characterized by fracture spacings (s1 fin, s2 fin) in two orthogonal directions of principal permeability, by a fracture opening parameter (efin), by fracture conductivities (Cf1 fin and Cf2 fin) and a permeability (km fin) of a matrix medium between fractures,- simplifying said equivalent fracture network (RFE) by means of a network fracture spacing coefficient (G) whose value is less than a value Gmax-zone defined on each of said zones in order to guarantee sufficient connectivity between simplified zones and non-simplified zones.
- A method as claimed in any one of the previous claims, wherein the following stages are carried out:- repeating stage a) while modifying said statistical parameters (PSF) so as to minimize a difference between a well test result and a well test simulation result from said simplified image,- associating with each one of said cells at least one equivalent permeability value and an average opening value for said fractures, said values being determined from said modified statistical parameters (PSF),- simulating fluid flows in said reservoir by means of a flow simulator, of the equivalent permeability values and the average opening values of said fractures associated with each one of said cells,- selecting a production scenario allowing the reservoir production to be optimized by means of said fluid flow simulation, and- developing said reservoir according to said scenario allowing the reservoir production to be optimized.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1004398A FR2967200A1 (en) | 2010-11-10 | 2010-11-10 | METHOD FOR CHARACTERIZING FRACTURE FRACTURE NETWORK AND METHOD FOR OPERATING IT |
Publications (2)
Publication Number | Publication Date |
---|---|
EP2453106A1 EP2453106A1 (en) | 2012-05-16 |
EP2453106B1 true EP2453106B1 (en) | 2017-01-18 |
Family
ID=44123420
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP11290475.0A Not-in-force EP2453106B1 (en) | 2010-11-10 | 2011-10-14 | Method for characterising the network of fractures of a fractured deposit and method for exploiting same |
Country Status (3)
Country | Link |
---|---|
US (1) | US8983818B2 (en) |
EP (1) | EP2453106B1 (en) |
FR (1) | FR2967200A1 (en) |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2843469A1 (en) * | 2011-07-28 | 2013-01-31 | Schlumberger Canada Limited | System and method for performing wellbore fracture operations |
WO2014148925A1 (en) * | 2013-03-22 | 2014-09-25 | Auckland Uniservices Limted | Method and system for monitoring and/or controlling fracture connectivity |
FR3007165B1 (en) * | 2013-06-13 | 2016-10-28 | Ifp Energies Now | METHOD FOR OPTIMIZING THE OPERATION OF A FLUID DEPOSITION BY TAKING INTO ACCOUNT A TERM OF GEOLOGICAL AND TRANSIENT EXCHANGE BETWEEN MATRIX BLOCKS AND FRACTURES |
US9416630B2 (en) * | 2013-06-18 | 2016-08-16 | Baker Hughes Incorporated | Grid calculation for fluid flow modeling |
US9523275B2 (en) * | 2013-08-26 | 2016-12-20 | Halliburton Energy Services, Inc. | Identifying an axis of a stimulated reservoir volume for a stimulation treatment of a subterranean region |
WO2016080987A1 (en) * | 2014-11-19 | 2016-05-26 | Halliburton Energy Services, Inc. | Junction models for simulating proppant transport in dynamic fracture networks |
CN104462734B (en) * | 2015-01-14 | 2017-05-17 | 中国石油大学(华东) | Combustible ice exploitation experimental model measuring point distribution design method based on numerical simulation |
JP6554828B2 (en) * | 2015-03-06 | 2019-08-07 | 富士通株式会社 | Information processing apparatus, information processing method, and program |
US10415362B1 (en) * | 2015-06-08 | 2019-09-17 | DataInfoCom USA Inc. | Systems and methods for analyzing resource production |
FR3041026B1 (en) | 2015-09-15 | 2017-10-20 | Ifp Energies Now | METHOD FOR CHARACTERIZING THE NETWORK OF FRACTURES OF A FRACTURE SLOT AND METHOD FOR OPERATING IT |
US10337294B2 (en) | 2015-10-01 | 2019-07-02 | Baker Hughes, A Ge Company, Llc | Reservoir permeability upscaling |
FR3045868B1 (en) | 2015-12-17 | 2022-02-11 | Ifp Energies Now | METHOD FOR CHARACTERIZING AND EXPLOITING AN UNDERGROUND FORMATION COMPRISING A NETWORK OF FRACTURES |
US10267132B2 (en) * | 2015-12-21 | 2019-04-23 | Baker Hughes, A Ge Company, Llc | Eliminating discrete fracture network calculations by rigorous mathematics |
GB2562942B (en) * | 2016-02-29 | 2021-10-27 | Landmark Graphics Corp | Hybrid 3D geocellular representation of selected natural fracture network subsets |
US10329904B2 (en) * | 2016-03-18 | 2019-06-25 | Baker Hughes, A Ge Company, Llc | Determining the robustness of discrete fracture network permeability estimates |
US10853533B2 (en) | 2016-05-09 | 2020-12-01 | Schlumberger Technology Corporation | Three-dimensional fracture abundance evaluation of subsurface formation based on geomechanical simulation of mechanical properties thereof |
US10650107B2 (en) | 2016-05-09 | 2020-05-12 | Schlumberger Technology Corporation | Three-dimensional subsurface formation evaluation using projection-based area operations |
WO2017222509A1 (en) * | 2016-06-22 | 2017-12-28 | Schlumberger Technology Corporation | Visualizations of reservoir simulations with fracture networks |
CN106373193B (en) * | 2016-09-22 | 2019-03-01 | 河海大学 | A kind of side slope assembled discrete element analysis generation method based on box cladding process |
CN106567703B (en) * | 2016-10-08 | 2018-10-12 | 中国石油大学(华东) | A kind of cloth hole optimization method of more radial hole auxiliary pressure breaks |
US10866962B2 (en) | 2017-09-28 | 2020-12-15 | DatalnfoCom USA, Inc. | Database management system for merging data into a database |
CN110185426B (en) * | 2018-02-22 | 2022-07-22 | 中国石油化工股份有限公司 | Shale oil reservoir multistage fracture conductivity optimization method |
CN109101773A (en) * | 2018-09-27 | 2018-12-28 | 北京科技大学 | Fine and close grease horizontal well solid seam net cluster network pressure splits optimization method |
CN114973891B (en) * | 2021-02-22 | 2024-04-09 | 中国石油化工股份有限公司 | Three-dimensional visual fracture-crack hole model and manufacturing method thereof |
US11933165B2 (en) | 2021-03-15 | 2024-03-19 | Saudi Arabian Oil Company | Hydraulic fracture conductivity modeling |
CN113236218B (en) * | 2021-05-31 | 2023-02-10 | 中国石油天然气股份有限公司 | Seam spacing determination method and device based on seam control fracturing reserve evaluation |
CN116181324B (en) * | 2023-03-10 | 2024-02-23 | 重庆科技学院 | Method for evaluating equivalent permeability of reservoir after fracturing |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2733073B1 (en) * | 1995-04-12 | 1997-06-06 | Inst Francais Du Petrole | METHOD FOR MODELING A LAMINATED AND FRACTURED GEOLOGICAL ENVIRONMENT |
US5619945A (en) * | 1996-01-16 | 1997-04-15 | The United States Of America As Represented By The Secretary Of The Navy | Drag reduction polymer ejection system for underwater vehicle |
FR2757947B1 (en) | 1996-12-30 | 1999-01-29 | Inst Francais Du Petrole | METHOD FOR DETERMINING THE EQUIVALENT PERMEABILITY OF A FRACTURE NETWORK IN A MULTI-LAYERED UNDERGROUND |
FR2757957B1 (en) | 1996-12-30 | 1999-01-29 | Inst Francais Du Petrole | METHOD FOR SIMPLIFYING THE MODELING OF A POROUS GEOLOGICAL ENVIRONMENT CROSSED BY AN IRREGULAR FRACTURE NETWORK |
FR2787219B1 (en) | 1998-12-11 | 2001-01-12 | Inst Francais Du Petrole | METHOD FOR MODELING FLUID FLOWS IN A CRACKED MULTI-LAYER POROUS MEDIUM AND CORRELATIVE INTERACTIONS IN A PRODUCTION WELL |
FR2809494B1 (en) * | 2000-05-26 | 2002-07-12 | Inst Francais Du Petrole | METHOD FOR MODELING FLOWS IN A FRACTURE MEDIUM CROSSED BY LARGE FRACTURES |
US7925482B2 (en) * | 2006-10-13 | 2011-04-12 | Object Reservoir, Inc. | Method and system for modeling and predicting hydraulic fracture performance in hydrocarbon reservoirs |
FR2918179B1 (en) | 2007-06-29 | 2009-10-09 | Inst Francais Du Petrole | METHOD FOR ESTIMATING THE PERMEABILITY OF A FRACTURE NETWORK FROM A CONNECTIVITY ANALYSIS |
CA2690992C (en) * | 2007-08-24 | 2014-07-29 | Exxonmobil Upstream Research Company | Method for predicting well reliability by computer simulation |
FR2923930B1 (en) | 2007-11-19 | 2009-11-20 | Inst Francais Du Petrole | METHOD FOR MODELING A POROUS GEOLOGICAL ENVIRONMENT CROSSED BY A FRACTURE NETWORK |
FR2925726B1 (en) * | 2007-12-20 | 2010-04-23 | Inst Francais Du Petrole | METHOD FOR OPTIMIZING THE OPERATION OF A FLUID DEPOSITION BY TAKING INTO ACCOUNT A TERM OF GEOLOGICAL AND TRANSIENT EXCHANGE BETWEEN MATRIX BLOCKS AND FRACTURES |
US8498852B2 (en) * | 2009-06-05 | 2013-07-30 | Schlumberger Tehcnology Corporation | Method and apparatus for efficient real-time characterization of hydraulic fractures and fracturing optimization based thereon |
-
2010
- 2010-11-10 FR FR1004398A patent/FR2967200A1/en active Pending
-
2011
- 2011-10-14 EP EP11290475.0A patent/EP2453106B1/en not_active Not-in-force
- 2011-11-03 US US13/288,049 patent/US8983818B2/en active Active
Non-Patent Citations (1)
Title |
---|
None * |
Also Published As
Publication number | Publication date |
---|---|
EP2453106A1 (en) | 2012-05-16 |
US8983818B2 (en) | 2015-03-17 |
FR2967200A1 (en) | 2012-05-11 |
US20120116740A1 (en) | 2012-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2453106B1 (en) | Method for characterising the network of fractures of a fractured deposit and method for exploiting same | |
EP3144468B1 (en) | Method for characterising the network of fractures of a fractured deposit and method for exploiting same | |
EP2581767B1 (en) | Method for constructing a mesh of a fractured reservoir with a limited number of nodes in the matrix environment | |
EP2212722B1 (en) | Method of modelling a porous geological environment through which a network of fractures runs | |
EP2530493B1 (en) | Method for building a fracture network mesh from a Voronoi diagram | |
EP2072752B1 (en) | Method for optimising the exploitation of a fluid reservoir by taking into consideration a geological and transitional exchange term between matrix blocks and fractures | |
EP2037080B1 (en) | Method for determining the permeability of a network of fractures based on a connectivity analysis | |
EP2466338B1 (en) | Method for exploiting an oil deposit based on a flow model defined by setting pressure and saturation maps to scale | |
EP1600897B1 (en) | Method of generating a hybrid 3D grid a heterogenous formation crossed by a geometric discontinuity in order to carry out simulations | |
CA2854085C (en) | Method for optimising the exploitation of a fluid deposit by taking into consideration a geological and transitional exchange term between matrix blocks and fractures | |
EP3199749A1 (en) | Method for exploiting a fluid reservoir having fractures passing through same by means of a flow simulation based on an exchange flow and a corrective factor | |
FR3045868A1 (en) | METHOD FOR CHARACTERIZING AND EXPLOITING UNDERGROUND FORMATION COMPRISING A FRACTURE NETWORK | |
EP2924232A2 (en) | Method for constructing an optimised network for reservoir simulation in a subterranean formation | |
EP2365359B1 (en) | Method of history matching a geological model | |
EP2806101B1 (en) | Method for exploiting a fracture medium from a wedge reservoir model for wells selected using an equivalent transmissivity pattern | |
EP4112875B1 (en) | Method for simulating the flows of a fluid in an underground formation comprising a network of cracks | |
EP2594733A1 (en) | Method for optimising the exploitation of an underground medium by means of a reservoir study comprising optimised scaling | |
EP3252508B1 (en) | Method for exploiting hydrocarbons from an underground formation, by means of optimal scaling |
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): AL 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 RS SE SI SK SM TR |
|
AX | Request for extension of the european patent |
Extension state: BA ME |
|
17P | Request for examination filed |
Effective date: 20121116 |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
INTG | Intention to grant announced |
Effective date: 20160909 |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): AL 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 RS 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: AT Ref legal event code: REF Ref document number: 863036 Country of ref document: AT Kind code of ref document: T Effective date: 20170215 |
|
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: DE Ref legal event code: R096 Ref document number: 602011034441 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: NL Ref legal event code: MP Effective date: 20170118 |
|
REG | Reference to a national code |
Ref country code: LT Ref legal event code: MG4D |
|
REG | Reference to a national code |
Ref country code: NO Ref legal event code: T2 Effective date: 20170118 |
|
REG | Reference to a national code |
Ref country code: AT Ref legal event code: MK05 Ref document number: 863036 Country of ref document: AT Kind code of ref document: T Effective date: 20170118 |
|
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: 20170118 |
|
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: 20170118 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: 20170118 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: 20170419 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: 20170118 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: 20170518 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
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: 20170118 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: 20170118 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: 20170118 Ref country code: RS 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: 20170118 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: 20170418 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: 20170518 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: 20170118 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: 20170118 |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: PLFP Year of fee payment: 7 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R097 Ref document number: 602011034441 Country of ref document: DE |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
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: 20170118 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: 20170118 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: 20170118 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: 20170118 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: 20170118 |
|
RAP2 | Party data changed (patent owner data changed or rights of a patent transferred) |
Owner name: IFP ENERGIES NOUVELLES |
|
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 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
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: 20170118 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: 20170118 |
|
26N | No opposition filed |
Effective date: 20171019 |
|
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: 20170118 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R119 Ref document number: 602011034441 Country of ref document: DE |
|
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: 20170118 |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: PL |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: MM4A |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: DE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20180501 Ref country code: CH Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20171031 Ref country code: LU Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20171014 Ref country code: LI Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20171031 |
|
REG | Reference to a national code |
Ref country code: BE Ref legal event code: MM Effective date: 20171031 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: BE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20171031 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
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: 20170118 |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: PLFP Year of fee payment: 8 |
|
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: 20171014 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: HU Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO Effective date: 20111014 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: CY Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20170118 |
|
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: 20170118 |
|
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: 20170118 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: AL 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: 20170118 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: NO Payment date: 20211020 Year of fee payment: 11 Ref country code: GB Payment date: 20211026 Year of fee payment: 11 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: FR Payment date: 20211027 Year of fee payment: 11 |
|
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: 20221014 |
|
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: 20221031 Ref country code: FR Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20221031 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: GB Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20221014 |