EP2253797B1 - Méthode d'exploitation de milieu poreux au moyen d'une modélisation d'écoulements de fluide - Google Patents
Méthode d'exploitation de milieu poreux au moyen d'une modélisation d'écoulements de fluide Download PDFInfo
- Publication number
- EP2253797B1 EP2253797B1 EP10290174.1A EP10290174A EP2253797B1 EP 2253797 B1 EP2253797 B1 EP 2253797B1 EP 10290174 A EP10290174 A EP 10290174A EP 2253797 B1 EP2253797 B1 EP 2253797B1
- Authority
- EP
- European Patent Office
- Prior art keywords
- well
- mesh
- reservoir
- simulator
- fluids
- 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.)
- Active
Links
- 239000012530 fluid Substances 0.000 title claims description 85
- 238000000034 method Methods 0.000 title claims description 50
- 238000004519 manufacturing process Methods 0.000 claims description 44
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 40
- 238000002347 injection Methods 0.000 claims description 38
- 239000007924 injection Substances 0.000 claims description 38
- 238000005553 drilling Methods 0.000 claims description 35
- 229920000642 polymer Polymers 0.000 claims description 24
- 239000000243 solution Substances 0.000 claims description 14
- 239000002253 acid Substances 0.000 claims description 8
- 230000009545 invasion Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000000638 stimulation Effects 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims 1
- 238000004088 simulation Methods 0.000 description 75
- 238000010168 coupling process Methods 0.000 description 57
- 230000008878 coupling Effects 0.000 description 40
- 238000005859 coupling reaction Methods 0.000 description 40
- 238000013459 approach Methods 0.000 description 30
- 239000012088 reference solution Substances 0.000 description 24
- 230000035699 permeability Effects 0.000 description 19
- 239000007789 gas Substances 0.000 description 14
- 230000004907 flux Effects 0.000 description 10
- 230000015572 biosynthetic process Effects 0.000 description 7
- 239000003921 oil Substances 0.000 description 7
- 230000002265 prevention Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 241001080024 Telles Species 0.000 description 4
- 229930195733 hydrocarbon Natural products 0.000 description 4
- 150000002430 hydrocarbons Chemical class 0.000 description 4
- 240000008042 Zea mays Species 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 235000021183 entrée Nutrition 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000706 filtrate Substances 0.000 description 1
- 238000001595 flow curve Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/12—Methods or apparatus for controlling the flow of the obtained fluid to or in wells
Definitions
- the present invention relates to the field of the exploitation of underground environments.
- the invention makes it possible in particular to improve the injectivity and the productivity of wells drilled through a porous medium, such as a hydrocarbon deposit or a geological CO 2 storage tank.
- the numerical methods which make it possible to model the flow of fluids within a well involve the construction of two distinct models: the model of reservoir (“ reservoir model ”) and the model d ' near-wellbore model'.
- domain decomposition techniques have been developed, described for example in GAIFFE, S. "Hybrid Meshes and Domain Decomposition for the Modeling of Petroleum Reservoirs", Doctoral Thesis, University of Paris 6, 2000 And windowing techniques ( “windowing") as described for example in the following document: MLACNIK, MJ and HEINEMANN, ZE "Using well windows in full field reservoir simulation ", paper SPE 66371 presented at the SPE Reservoir Simulation Symposium, Houston, TX, USA, February 2001 .
- each successive time interval can have a length which is a function of a time step of calculation of the first flow simulator and of a time step of the second flow simulator.
- each successive time interval can have a length equal to a time step of the first flow simulator.
- the fluid flows within the medium are simulated by means of the first simulator on a first mesh discretizing the porous medium into a set of meshes
- the fluid flows around the well are simulated by means of the second simulator on a second mesh discretizing the well and its surroundings into a set of meshes.
- This second mesh is generated by constraining meshes situated on the edge of the second mesh, so that their interfaces coincide with the interfaces of the meshes of the first mesh.
- digital productivity index multipliers are updated, instead of the digital productivity indices themselves, for each phase, by comparing flow rates per phase calculated by the first simulator and flow rates per phase calculated by the second simulator.
- damage to the well can be taken into account by a drilling fluid by modeling an invasion of the porous reservoir by the drilling fluid in steps d and e.
- the operating scenario can include an injection of a polymer solution through the well, and we can then model the flows to prevent water coming in.
- the operating scenario can also include an injection of an acid solution into the well, and we can then model the flows to assess the impact of acid stimulation.
- the invention relates to a method for exploiting an underground porous medium, by injecting a fluid into the medium via at least one well, and / or by producing a fluid present in the medium by means of at least one well.
- the method includes a modeling of the flows of fluids in the system constituted by the porous medium (reservoir and surroundings of the wells) and the well. It is therefore, in particular, to model the injectivity or the productivity of wells passing through a porous medium.
- a scenario can be a scenario of production of hydrocarbons contained in the porous medium (reservoir), or a scenario of injection of acid gas, such as CO2, in an underground reservoir for the storage of acid gas .
- a scenario is characterized by the position of the wells, the recovery or injection method, the injection and / or production flow rates and duration, the operating conditions in these wells such as the flow rate or the bottom pressure.
- the reservoir engineer chooses a production process, for example the recovery process by water injection, of which it then remains to specify the optimal implementation scenario for the reservoir considered.
- the definition of an optimal scenario consists, for example, in fixing the number and the location (position and spacing) of the injector and producer wells in order to best take into account the impact of heterogeneities within the reservoir, for example permeability channels, fractures, etc., on the progression of fluids within the reservoir.
- the flow simulator we are then able to simulate the expected production of hydrocarbons, using the tool well known to specialists: the flow simulator.
- the “reservoir mesh” consists of a set of meshes spatially discretizing the reservoir (porous medium + well).
- An example of a tank mesh is illustrated on the figure 3 , this mesh is coarse. Some meshes correspond to the “porous medium” part, others correspond to the part where the well is drilled. One speaks for these last of meshs of well of the mesh of reservoir.
- the “well access grid” consists of a set of meshes spatially discretizing the well and its surroundings.
- An example of first well mesh is illustrated on the figure 4 , this mesh is fine to simulate the detailed phenomena around the well. Its surroundings therefore belong to the porous medium in which the well is drilled. Some meshes correspond to the “porous medium” part, others correspond to the “well” part.
- an object of the invention relates to a coupling method, which makes it possible to couple in a very simple way a reservoir model, for the simulation of the reservoir, and a well model first, which is an autonomous model for simulating the phenomena detailed around the well.
- tank model simulator it can be the Puma Flow ® software (IFP, France) for example.
- the technique used here consists in coupling between the two flow simulators.
- a coarse mesh is often used for the reservoir model, and a fine mesh is usually necessary to simulate the detailed phenomena around the well.
- the figure 5 shows the two meshes used in the coupling.
- the figure on the left represents the reservoir mesh for the field simulation, and the figure on the right represents the mesh in the vicinity of wells in the approach well model.
- the meshes at the edges (in gray) in the well approach model coincide with the meshes of the same color in the reservoir mesh.
- the cross indicates the location of the well.
- the time steps used in the first well model are generally much smaller than those of the reservoir model.
- the reservoir model is mainly used to simulate flows in the entire reservoir.
- the digital productivity index IP takes into account: the geometric effect of the well mesh i of the mesh, the permeability of the porous medium in the well mesh and a skin coefficient.
- a skin coefficient is a coefficient, well known to those skilled in the art, used to represent the damage of a well in a mesh.
- the variables IP i , P nw, p, j , P r, p, i and P wf, j are a function of time T.
- the optimal scenario can be selected.
- the optimal scenario is the scenario allowing to obtain an optimal production of the deposit within the framework of the production of a reservoir, or the scenario allowing to obtain the optimal injectivity in the deposit within the framework of injection of fluid in the tank (injection of water for improved production, or injection of acid gases).
- the scenario selected in step 1 is modified ( ⁇ SCE ), for example by modifying the location of a well.
- step 2 during which the meshes are constructed, is modified.
- the simulation carried out using the reservoir model in step 3c, provides dynamic properties of fluids such as pressure or saturation in the period from T 0 to T 1 on all coarse meshes.
- the determination of the boundary conditions in step 3b requires the interpolation of the pressure or the flow on the edges of the well approach model.
- the edge meshes in the approach well model are also constrained so that they coincide with meshes of the reservoir model ( figure 3 ). In this way, the transfer of dynamic data from the reservoir model to the approach well model is direct on these meshes.
- the boundary conditions are of zero flux. In order to maintain the dynamic properties at the edges of the model, very high porosities
- M p, i is the multiplier of the productivity index for phase p in the mesh of wells i .
- the coupling method according to the invention can be used to model various detailed phenomena around the well, such as, damage by drilling or completion fluid, acid stimulation, non-Darcean flow around the well, the problem of condensate gas, asphaltene deposition, damage by injection of CO 2 , prevention of water or gas, sand, mineral deposits, the impact of completions, etc.
- damage by drilling or completion fluid acid stimulation, non-Darcean flow around the well
- the problem of condensate gas asphaltene deposition
- damage by injection of CO 2 prevention of water or gas, sand, mineral deposits
- the impact of completions etc.
- we present in particular an example of application for the damage of the oil formation by the drilling fluid during the drilling of the well and an example of application for the prevention of water coming when a well in production produces a significant amount of water, and we are trying to reduce this production of water.
- a standard reservoir model is used for the field simulation.
- a 1000m x 1000m x 10m tank is considered.
- a Cartesian mesh with 20 meshes in the x direction, 20 meshes in the y direction and 1 mesh in the z direction is used for the simulation of the field ( figure 6 ).
- the mesh sizes are therefore 50m x 50m x 10m.
- the initial tank pressure is 200 bar.
- a producing well must be drilled in the block (15, 15, 1). It is represented by a black circle on the figure 6 . The damage to this well by the drilling fluid is studied with the method according to the invention.
- the tank is homogeneous with permeability 200 mD and porosity 0.15.
- the boundary conditions of this reservoir are zero flows, except on the edge ⁇ x- ( figure 6 ), where the pressure is constant (200 bars).
- the mesh is refined around the well ( figure 7 ).
- a specific model which takes into account the advanced physics of the damage, is used on this mesh to simulate the reference solution. Since damage by drilling fluid is generally limited to a few centimeters or a few tens of centimeters around wells, we need very small meshes in the refined zone (Table 1).
- the diameter of the well is 21.6 cm.
- the size of the well mesh is 22 cm.
- the other meshes around the well are much smaller with a size of 2 cm.
- the meshes used for the coupling are illustrated on the Figures 8A and 8B .
- the mesh of the well approach model ( figure 8B ) corresponds to the refined zone and to the meshes around in the reference mesh.
- the meshes at the edges of the approaching well model coincide with meshes of the reservoir model.
- the contact time between the drilling fluid and the reservoir is 2 days.
- the pressure during drilling at the bottom of the well is 250 bars.
- the permeability and the thickness of the external “cake” formed by the drilling mud are equal to 0.001 mD and 0.2 cm.
- the thickness of the internal cake is 2 cm with an average permeability reduced to 20 mD during the drilling period and 40 mD in the production period.
- the viscosity of the drilling fluid is 30 cPo.
- the hysteresis of the relative permeability between the drilling and production periods is presented in the figure 9 . An irreducible water saturation of 30% linked to the filtrate (drilling fluid) which will invade the formation during the drilling phase will remain blocked in the porous medium when the well is returned to production.
- the drilling fluid invasion volumes are compared to the figure 10 for the simulation with the coupling method and the reference solution obtained using the mesh with local refinement ( figure 7 ).
- the time steps for updating the data in the coupling are presented in Table 2.
- the figure 10 shows that the volume of fluid invasion is correctly simulated with the coupling method.
- the small difference between the coupling solution and the reference solution in the period between 0.1 and 0.3 days can be improved by using small iteration steps in time to exchange the data in the coupling.
- the well After the 2 days of drilling, the well is closed for 1 day for the establishment of its completion, then it is put into production. The coupling is effected until 10 th day. Beyond 10 days, the effect of damage around wells becomes stable and the digital IPs in the reservoir model hardly change any more. We no longer need coupling to continue simulating the field with the reservoir model.
- the oil production curve simulated by the reservoir model, which is coupled with the first well model for the first 10 days, is presented in the figure 11 . This curve is very close to the reference solution.
- a polymer solution is injected into a producing well for a short time in order to reduce the large amount of water produced together with the oil. Part of the polymer is absorbed on the rock, and another part is dispersed in water.
- the injected polymer has the effect of reducing the mobility of the water phase by increasing its viscosity and by reducing the relative permeability of this phase. Therefore, in the coupling method, the most suitable approach is to update the digital IP multiplier for the water phase.
- a 1000m x 1000m x 25m tank is considered as an example.
- a Cartesian mesh with 20 meshes in the x direction, 20 meshes in the y direction and 5 meshes in the z direction is used for the simulation of the field.
- the mesh size is 50m x 50m x 5m.
- the reservoir is heterogeneous. Permeability is presented to the figure 12 .
- the ratio of permeabilities in the vertical and horizontal directions is 0.1.
- the initial tank pressure is 200 bar.
- the pressure at the injector well is imposed at 300 bars, and the pressure at the producing well is constrained to 150 bars during production.
- the water-cut (water flow compared to the total flow) of the producing well reaches 85%.
- the water ingress prevention procedure is then applied to reduce the amount of water produced.
- a polymer solution with a concentration of 2500 ppm is injected into the producer with a bottom pressure of 300 bars for 2 days. Then, the well is put back into production. This water prevention procedure is simulated with the method according to the invention.
- a local refinement around the producing well is used ( figure 13 ).
- the mesh size around the well is 0.617 m in the x direction.
- the mesh for coupling is presented in the figure 14 .
- Meshes at the edges of the first well model coincide with meshes in the reservoir model.
- the physics of the polymer can be considered in both models (first well model and reservoir model).
- Table 3 No time during coupling Period (day) No time (day) 0 - 950 - 950-970 2 970 - 1000 28 1000 - 1000.1 0.01 1000.1 - 1005 0.1 1005 - 1030 1 1030 - 1100 2 1100 - 3000 -
- the coupling begins at 950 days and ends at 1100 days, a period of 150 days in total.
- the time steps for exchanging data in the coupling method are presented in Table 3.
- the global digital IPs are updated at the start of the coupling (from 950 to 970 days) to take into account the effects of the meshes between the reservoir model and the approaching well model.
- the global digital IPs are further recalculated to integrate the effect induced by the injected polymer (we could also update the digital IP multipliers for the water phase ).
- the digital IP multipliers for the water phase are updated.
- the figure 15 compare the injection rates of polymer in the well for the different simulations: the reference solution, the simulation on the reservoir mesh with coupling, the direct simulation on the reservoir mesh without coupling and the simulation with the well approach model (with coupling).
- the Figures 16A to 16E show the same layer-by-layer comparisons.
- the volume of polymer injected is greatly overestimated.
- the results are significantly improved.
- the injection rate is high, but it is quickly corrected by the update of the IP due to the coupling. If we want to have more precision on the polymer injection rate, it is enough to refer to the simulation results with the first model of well. With this model, the volume injected and the distribution of the polymer around the well are both correctly simulated.
- the figures 17, 18 and 19 show the oil, water and water-cut flow curves for the tank model with coupling, the tank model without coupling and the reference solution.
- the results of the coupled reservoir model are generally satisfactory.
- the figure 20 presents the water saturation card at the end of the coupling (1100 days), and the figure 21 shows the pressure card at 1100 days. Compared to the reference solutions, the coupling gives globally satisfactory results.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
- La présente invention concerne le domaine de l'exploitation de milieux souterrains. L'invention permet notamment d'améliorer l'injectivité et la productivité de puits forés à travers un milieu poreux, tel qu'un gisement d'hydrocarbures ou un réservoir de stockage géologique de CO2.
- Les phénomènes locaux qui peuvent se produire autour d'un puits, tels que l'endommagement, ont un impact énorme sur l'injectivité ou la productivité d'un puits. Dans l'industrie pétrolière, il est très important de prédire l'injectivité ou la productivité, surtout lorsqu'il y a des altérations de formation au voisinage de puits, qui changent la capacité d'injection ou de production du puits.
- Depuis des années, des efforts conséquents sont effectués au moyen de techniques expérimentales, en laboratoire, ou de méthodes numériques de modélisations, afin de prendre en compte ces phénomènes locaux autour des puits, et leurs impacts sur l'injectivité ou la productivité.
- Les méthodes numériques qui permettent de modéliser l'écoulement de fluides au sein d'un puits (injectivité et productivité d'un puits), comportent la construction de deux modèles distincts : le modèle de réservoir (« reservoir model ») et le modèle d'abords de puits (« near-wellbore model »).
- Un modèle de réservoir comporte deux éléments :
- un maillage, dit « maillage de réservoir », constitué d'un ensemble de mailles discrétisant spatialement le réservoir.
- un simulateur d'écoulement. On appelle simulateur d'écoulement, un logiciel permettant de modéliser les écoulements de fluides au sein d'un milieu poreux, grâce au maillage de réservoir. Ce logiciel permet ainsi de simuler des données/propriétés dynamiques des fluides (eau, huile, gaz) : pression, flux (quantité de matière traversant une surface), saturation, débits, concentrations. Par exemple, un simulateur permet d'estimer, pour un scénario d'exploitation de puits donné (scénario de production ou scénario d'injection), et pour un intervalle de temps donné : les saturations en eau, huile et gaz, les débits d'huile, de gaz et d'eau, le water-cut (fraction d'eau dans la production liquide), le GOR (rapport de gaz et huile dans la production), les concentrations en polymère absorbé sur la roche du milieu poreux, les débits d'injection de polymère, si une solution de polymère est injectée dans le réservoir par l'intermédiaire de puits d'injection, ...
- Un modèle d'abords de puits comporte deux éléments :
- un maillage, dit « maillage d'abords de puits », constitué d'un ensemble de mailles discrétisant spatialement le puits et ses abords. Ses abords appartiennent donc au milieu poreux dans lequel le puits est foré.
- un simulateur d'écoulement, permettant de simuler, grâce au maillage d'abords de puits, des données/propriétés dynamiques des fluides (eau, huile, gaz).
- Ces deux types de modèles, réservoir et abord de puits, sont en général autonomes et découplés. Les phénomènes locaux étant généralement limités au voisinage immédiat du puits (de quelques centimètres à quelques mètres), des très petites mailles sont nécessaires pour le maillage d'abord de puits, alors que des mailles de tailles plus grandes sont utilisées pour les maillages de réservoir de façon à accélérer les calculs.
- On connaît des techniques permettant d'utiliser un seul et même simulateur d'écoulement de réservoir pour ces deux maillages. On peut par exemple utiliser la technique des maillages dits « hybrides », combinant en un seul maillage, des mailles pour le maillage du réservoir et des mailles pour un maillage localement affiné autour du puits. On associe à ce type de maillage un seul simulateur d'écoulement pour mieux prendre en compte les comportements des écoulements au voisinage du puits dans une simulation du champ.
- Mais les simulations simultanées des écoulements dans le réservoir, qui nécessitent un très grand nombre de mailles, et dans les régions proches du puits avec de petites mailles, qui nécessitent des petits pas de temps pour assurer la stabilité du calcul, posent des problèmes numériques, en particulier, le problème du temps de calcul (temps CPU).
- Ainsi, ont été mises au point des techniques de décomposition de domaine, décrites par exemple dans GAIFFE, S. "Maillages Hybrides et Décomposition de Domaine pour la Modélisation des Réservoirs Pétroliers", Thèse de Doctorat, Université Paris 6, 2000, et des techniques de fenêtrage ("windowing"), décrites par exemple dans le document suivant :MLACNIK, M.J. and HEINEMANN, Z.E. "Using well windows in full field reservoir simulation", paper SPE 66371 presented at the SPE Reservoir Simulation Symposium, Houston, TX, U.S.A., February 2001.
- Mais, certains points délicats comme la convergence, la stabilité ou le temps de calcul posent des problèmes lors d'applications industrielles. De plus, la méthode de décomposition de domaine n'est pas toujours "conservative" (détérioration du bilan de masse dans le modèle en fonction du temps), ce qui ne convient pas à une utilisation pratique de la méthode. Par ailleurs, toutes ces techniques demandent de reformuler les équations mathématiques et les conditions aux limites développées dans les simulateurs d'écoulement, et de nouveaux développements sont nécessaires pour intégrer les solutions proches et éloignées du puits dans un seul modèle, dont la réalisation est un travail long et difficile.
- Un objet de l'invention concerne une méthode mise en oeuvre par ordinateur pour modéliser des écoulements de fluides au sein d'un milieu poreux souterrain traversé par au moins un puits, en vue de l'exploitation dudit milieu réalisée en injectant un desdits fluides dans ledit milieu via au moins un desdits puits, et/ou en produisant un desdits fluides présent dans ledit milieu au moyen d'au moins un desdits puits. La méthode comporte l'utilisation d'un premier simulateur d'écoulement permettant de simuler l'écoulement des fluides au sein du milieu poreux à partir d'indices de productivité numériques reliant des pressions de fluides à des débits de fluides, et l'utilisation d'un second simulateur d'écoulement pour simuler l'écoulement des fluides aux abords du puits à partir de conditions aux limites. La méthode comporte les étapes suivantes :
- a- on simule les écoulements de fluides au sein du milieu au moyen du premier simulateur sur un intervalle de temps défini entre des temps T0 et T1, et on déduit des conditions aux limites mises à jour pour le second simulateur par une interpolation linéaire des résultats du premier simulateur entre les temps T0 et T1 ;
- b- on simule les écoulements de fluides aux abords du puits au moyen du second simulateur sur le même intervalle de temps, en utilisant les conditions aux limites mises à jour, et on en déduit des indices de productivité numériques mis à jour pour le premier simulateur en comparant des débits calculés par le premier simulateur et des débits calculés par le second simulateur ; et
- c- on modélise les écoulements de fluides au sein du milieu poreux pendant une période de temps entre T0 et Tn, où Tn >T1, en réitérant les étapes a et b, pour des intervalles de temps successifs compris entre T0 et Tn.
- Selon l'invention, chaque intervalle de temps successif peut avoir une longueur qui est fonction d'un pas de temps de calcul du premier simulateur d'écoulement et d'un pas de temps du second simulateur d'écoulement. Par exemple, chaque intervalle de temps successifs peut avoir une longueur égale à un pas de temps du premier simulateur d'écoulement.
- Selon un mode de réalisation, on simule les écoulements de fluides au sein du milieu au moyen du premier simulateur sur un premier maillage discrétisant le milieu poreux en un ensemble de mailles, et on simule les écoulements de fluides aux abords du puits au moyen du second simulateur sur un second maillage discrétisant le puits et ses abords en un ensemble de mailles. Ce second maillage est généré en contraignant des mailles situées sur le bord du second maillage, de façon à ce que leurs interfaces coïncident avec les interfaces des mailles du premier maillage.
- Dans le cas où l'on modélise des écoulements multiphasiques, on met à jour des multiplicateurs d'indices de productivité numériques, au lieu des indices de productivité numériques eux-mêmes, pour chaque phase, en comparant des débits par phase calculés par le premier simulateur et des débits par phase calculés par le second simulateur.
- L'invention concerne également une méthode pour exploiter un réservoir poreux souterrain au moyen d'au moins un puits le traversant, au moins un fluide circulant entre le réservoir et le puits. Selon cette méthode, on acquiert des données relatives à la géométrie du réservoir poreux, à partir desquelles on construit une discrétisation du réservoir en un ensemble de mailles, appelé maillage de réservoir, et on construit une discrétisation du puits et de ses abords en un ensemble de mailles, appelé maillage d'abord de puits. Cette méthode comporte également les étapes suivantes :
- a- on choisit un scénario d'exploitation du réservoir poreux ;
- b- on associe au maillage de réservoir un premier simulateur d'écoulement permettant de simuler l'écoulement des fluides au sein du réservoir, à partir d'au moins les données suivantes: le scénario de production, des données d'entrées relatives au fluide et au réservoir, des indices de productivité numériques permettant de relier des pressions à des débits, des conditions aux limites ;
- c- on associe au maillage d'abord de puits un second simulateur d'écoulement pour simuler l'écoulement des fluides aux abords du puits, à partir d'au moins les données suivantes : des données d'entrées relatives au fluide et au réservoir, des conditions aux limites ;
- d- on modélise les écoulements de fluides au sein du milieu poreux et aux abords du puits, au moyen de la méthode selon l'une des revendications 1 à 7 ; et
- e- on modifie le scénario d'exploitation et l'on répète l'étape d, jusqu'à obtenir un scénario d'exploitation optimal.
- Selon cette méthode d'exploitation, on peut tenir compte d'un endommagement du puits par un fluide de forage en modélisant une invasion du réservoir poreux par le fluide de forage aux étapes d et e.
- Le scénario d'exploitation peut comporter une injection d'une solution de polymère par le puits, et l'on peut alors modéliser les écoulements pour prévenir une venue d'eau. Le scénario d'exploitation peut également comporter une injection d'une solution acide dans le puits, et l'on peut alors modéliser les écoulements pour évaluer l'impact d'une stimulation acide.
- D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après.
-
- La
figure 1 illustre les principales étapes de la méthode selon l'invention. - La
figure 2 montre le schéma de couplage entre le modèle de réservoir et le modèle abord de puits. L'axe T correspond au temps. - La
figure 3 montre le maillage grossier utilisé pour la simulation du champ dans un modèle de réservoir. - La
figure 4 montre le maillage fin utilisé pour simuler les phénomènes détaillés des écoulements autour de puits dans un modèle abord de puits. - La
figure 5 montre les deux maillages utilisés dans le couplage. La figure de gauche représente le maillage de réservoir pour la simulation de champ, et la figure de droite représente le maillage au voisinage de puits dans le modèle abord de puits. Les mailles aux bords (en gris) dans le modèle abord de puits coïncident avec les mailles de même couleur dans le maillage de réservoir. - La
figure 6 montre le maillage grossier pour la simulation du champ dans le cas de l'endommagement par le fluide de forage. Γx- et Γ x+ correspondent à deux frontières de ce maillage dans la direction x, et Γy- et Γ y+ correspondent à deux frontières de ce maillage dans la direction y. - La
figure 7 montre le maillage localement raffiné autour du puits pour simuler la solution de référence dans le cas de l'endommagement par fluide de forage. - Les
figures 8A et 8B montrent les maillages de couplage pour simuler l'endommagement par le fluide de forage. Le maillage de gauche (figure 8A ) correspond au maillage pour la simulation de champ, et le maillage de droite (figure 8B ) correspond au maillage dans le modèle abord de puits. - La
figure 9 montre les perméabilités relatives pendant le forage et la production. L'axe X est la saturation sans unité. L'axe Y est la perméabilité relative. Il n'y a pas d'unité. La courbe "krw drilling" est la courbe de perméabilité relative de l'eau pendant le forage. La courbe "kro drilling" est la courbe de perméabilité relative de l'huile pendant le forage. Les courbes "krw production" et "kro production" sont des courbes de perméabilité relative de l'eau et de l'huile respectivement, pendant la production. - Les
figures 10A et10B comparent le volume d'invasion du fluide de forage simulé par la méthode de couplage avec celui de la solution de référence. Lafigure 10A montre le débit d'invasion pendant le forage. L'axe X est le temps exprimé en jour. L'axe Y est le débit exprimé en m3/jour. La courbe R est la solution de référence. La courbe CM est la simulation avec la méthode de couplage. Lafigure 10B montre le volume d'invasion en fonction du temps. L'axe X est le temps en jour. L'axe Y est le débit en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation avec la méthode de couplage. - La
figure 11 compare les débits de production d'huile. L'axe X est le temps en jour. L'axe Y est le débit en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation avec la méthode de couplage. La courbe S est la simulation sans l'endommagement du puits. La courbe CK est la simulation avec l'endommagement uniquement par les boues de forage (« cakes »). - La
figure 12 montre les perméabilités dans la couche 3 sur le maillage grossier du modèle de réservoir dans l'application à la prévention de venue d'eau. Il y a un injecteur et un producteur. - La
figure 13 montre le maillage raffiné autour du producteur pour simuler la solution de référence. - Les
figures 14A et 14B montrent les maillages du couplage. Lafigure 14A est le maillage du modèle de réservoir, et lafigure 14B est le maillage du modèle abord de puits. - La
figure 15 montre le débit d'injection de polymère dans le puits traité. L'axe X est le temps en jour. L'axe Y est le débit en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation du modèle de réservoir avec la méthode de couplage. La courbe S est la simulation directe avec le modèle de réservoir sans couplage. La courbe NW est la simulation du modèle abord de puits avec couplage. - Les
figures 16A à 16E montrent le débit d'injection de polymère dans les couches. Lafigure 16A montre le débit d'injection de polymère dans la couche 1. Lafigure 16B montre le débit d'injection de polymère dans la couche 2. Lafigure 16C montre le débit d'injection de polymère dans la couche 3. Lafigure 16D montre le débit d'injection de polymère dans la couche 4. Lafigure 16E montre le débit d'injection de polymère dans la couche 5. L'axe X est le temps en jour. L'axe Y est le débit en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation du modèle de réservoir avec la méthode de couplage. La courbe S est la simulation directe avec le modèle de réservoir sans le couplage. La courbe NW est la simulation du modèle abord de puits avec (évidemment) le couplage. - La
figure 17 montre le débit d'huile du producteur. L'axe X est le temps en jour. L'axe Y est le débit d'huile en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation du modèle de réservoir avec la méthode de couplage. La courbe S est la simulation directe avec le modèle de réservoir sans le couplage. - La
figure 18 montre le débit d'eau du producteur. L'axe X est le temps en jour. L'axe Y est le débit d'eau en m3/j. La courbe R est la solution de référence. La courbe CM est la simulation du modèle de réservoir avec la méthode de couplage. La courbe S est la simulation directe avec le modèle de réservoir sans le couplage. - La
figure 19 montre la courbe de water-cut du producteur. L'axe X est le temps exprimé en jour. L'axe Y est le water-cut, sans unité. La courbe R est la solution de référence. La courbe CM est la simulation du modèle de réservoir avec la méthode de couplage. La courbe S est la simulation directe avec le modèle de réservoir sans le couplage. - Les
figures 20A à 20D montrent une carte de saturation en eau à 1100 jours. Lafigure 20A correspond à la solution de référence dans le champ. Lafigure 20B correspond à la carte obtenue avec le modèle de réservoir en maillage grossier avec couplage. Lafigure 20C montre la solution de référence au voisinage du puits. Lafigure 20D montre la saturation en eau au voisinage de puits simulée avec le modèle abord de puits. - Les
figures 21A à 21D montrent une carte de pression à 1100 jours. Lafigure 21A montre la solution de référence dans le champ. Lafigure 21B montre la solution obtenue avec le modèle de réservoir en maillage grossier avec couplage. Lafigure 21C montre la solution de référence au voisinage du puits. Lafigure 21D montre la solution avec le modèle abord de puits. - L'invention concerne une méthode pour exploiter un milieu poreux souterrain, en injectant un fluide dans le milieu via au moins un puits, et/ou en produisant un fluide présent dans le milieu au moyen d'au moins un puits également. La méthode comporte une modélisation des écoulements de fluides dans le système constitué par le milieu poreux (réservoir et abords des puits) et le puits. Il s'agit donc, en particulier, de modéliser l'injectivité ou la productivité de puits traversant un milieu poreux.
- La
figure 1 illustre les principales étapes de la méthode : - 1. Choix d'un scénario d'exploitation du milieu poreux, scénario de production et/ou scénario d'injection (SCE) ;
- 2. Choix d'un simulateur d'écoulement (RSIM) compatible avec un maillage de réservoir donné, et choix d'un simulateur d'écoulement (NWSIM) compatible avec un maillage d'abord de puits donné ;
- 3. Au moyen d'un couplage entre les deux simulateurs (EST_CAL, et
figure 2 ), estimation des écoulements de fluides, c'est-à-dire, par exemple, du volume injecté ou du volume produit, sur un intervalle de temps donné, et - 4. Détermination du scénario d'exploitation optimal par modification du scénario d'exploitation et répétition de l'étape 3 (OPT).
- Il peut s'agir d'un scénario de production des hydrocarbures contenus dans le milieu poreux (réservoir), ou d'un scénario d'injection de gaz acide, tel que le CO2, dans un réservoir souterrain en vue du stockage du gaz acide. Un scénario se caractérise par la position des puits, la méthode de récupération ou d'injection, les débits et durée d'injection et/ou de production, les conditions opératoires dans ces puits telles que le débit ou la pression de fond.
- Dans le cadre de la production, l'ingénieur réservoir choisit un procédé de production, par exemple le procédé de récupération par injection d'eau, dont il demeure ensuite à préciser le scénario optimal de mise en oeuvre pour le réservoir considéré. La définition d'un scénario optimal consiste, par exemple, à fixer le nombre et l'implantation (position et espacement) des puits injecteurs et producteurs afin de tenir compte au mieux de l'impact d'hétérogénéités au sein du réservoir, par exemple chenaux de perméabilité, fractures, etc., sur la progression des fluides au sein du réservoir. En fonction du scénario choisi, et de la représentation géométrique du réservoir, on est alors capable de simuler la production d'hydrocarbures escomptée, au moyen de l'outil bien connu des spécialistes : le simulateur d'écoulement.
- Le choix d'un scénario, par la définition de multiples caractéristiques techniques, est une étape bien connue des spécialistes.
- Pour choisir un simulateur d'écoulement, il est nécessaire de connaître le type de maillage sur lequel le simulateur doit fonctionner.
- Le « maillage de réservoir » est constitué d'un ensemble de mailles discrétisant spatialement le réservoir (milieu poreux + puits). Un exemple de maillage de réservoir est illustré sur la
figure 3 , ce maillage est grossier. Certaines mailles correspondent à la partie « milieu poreux », d'autres, correspondent à la partie où le puits est foré. On parle pour ces dernières de mailles de puits du maillage de réservoir. - Le « maillage d'abords de puits » est constitué d'un ensemble de mailles discrétisant spatialement le puits et ses abords. Un exemple de maillage d'abord de puits est illustré sur la
figure 4 , ce maillage est fin pour simuler les phénomènes détaillés autour du puits. Ses abords appartiennent donc au milieu poreux dans lequel le puits est foré. Certaines mailles correspondent à la partie « milieu poreux », d'autres, correspondent à la partie « puits ». On parle pour ces dernières de mailles de puits du maillage d'abords de puits. - La génération de ces maillages, que ce soit le maillage de réservoir ou le maillage d'abord de puits, est une étape bien connue du spécialiste, qui connaît de nombreuses méthodes pour les construire. Par exemple, des techniques de construction d'un maillage d'abord de puits sont décrites dans les documents suivants :
Boe, O., Flynn, J. et Reiso, E., "On Near Wellbore Modeling and Real Time Reservoir Management", SPE 66369, Houston, Texas, USA, 11-14 Feb. 2001. - On connaît également des méthodes pour construire des maillages de réservoir, à partir de données relatives à la géométrie du milieu (sismiques, diagraphies...), décrites par exemple dans les documents suivants :
Flandrin, N., Bennis, C. et Borouchaki, H., "3D Hybrid Mesh Génération for Reservoir Simulation", ECMOR, Canne, France, 30 August - 2 September, 2004. - Pour définir un modèle de réservoir, il est nécessaire d'associer un simulateur d'écoulement au maillage de réservoir. De même, pour définir un modèle d'abord de puits, il est nécessaire d'associer un simulateur d'écoulement au maillage d'abord de puits.
- Comme il est connu de l'homme de l'art, pour fonctionner, un simulateur d'écoulement a besoin de certaines données, dites données d'entrée :
- caractéristiques géométriques du réservoir, caractéristiques de la roche, caractéristiques des fluides en place et des fluides injectés (masse volumique, viscosité), courbes de perméabilités relatives, courbes de pression capillaire, saturations initiales en fluides,...
- conditions aux limites du domaine simulé et aux puits où sont injectés ou produits des fluides. Les conditions aux limites, sont les valeurs de données dynamiques, telles que la pression, le débit ou le flux, les saturations en fluides, sur les bords du maillage ou dans les mailles qui forment les bords du maillage de réservoir ou d'abord de puits. Un exemple de conditions aux limites peut être : un flux nul sur tous les bords du maillage, ou des saturations et pressions imposées sur les mailles aux bords du maillage.
- éventuellement des Indices de Productivité numériques (IP). La connexion entre la pression dans les mailles traversées par un puits et les pressions dans le puits lui-même est réalisée à l'aide d'un Indice de Productivité numérique (IP). L'IP numérique peut être calculé par une formule analytique dans le code ou donné par l'utilisateur du logiciel (le simulateur). En général, le simulateur calcule un IP numérique par une formule analytique au début de simulation. Mais, si l'utilisateur donne un IP numérique dans le jeu de donnée d'entrée, c'est l'IP numérique de l'utilisateur qui est prise en compte dans la simulation.
- Selon l'invention, il est possible d'utiliser tout type de simulateur d'écoulement, que ce soit pour le modèle de réservoir, ou le modèle d'abord de puits. En effet, un objet de l'invention concerne une méthode de couplage, qui permet de coupler de manière très simple un modèle de réservoir, pour la simulation du réservoir, et un modèle d'abord de puits, qui est un modèle autonome pour simuler les phénomènes détaillés autour du puits.
- Concernant le simulateur du modèle de réservoir, il peut s'agir du logiciel PumaFlow® (IFP, France) par exemple.
- Concernant le simulateur du modèle d'abord de puits, on peut utiliser celui décrit dans le document suivant : DING, Y., RENARD, G. : "Evaluation of Horizontal Well Performance after Drilling Induced Formation Damage", J. of Energy Resources Technology, Vol. 127, Sept., 2005.
- Il s'agit ici de modéliser l'injectivité ou la productivité d'un puits traversant le milieu poreux, et permettant l'exploitation de ce milieu. Cette modélisation se fait sur un intervalle de temps donné, D=[T0 ; Tn ]. Par exemple, on modélise le comportement du système milieu+puits sur 20 ans, en considérant le scénario d'exploitation préalablement choisi.
- La technique utilisée ici, consiste à réaliser un couplage entre les deux simulateurs d'écoulement.
- Un maillage grossier est souvent utilisé pour le modèle de réservoir, et un maillage fin est habituellement nécessaire pour simuler les phénomènes détaillés autour du puits. La
figure 5 montre les deux maillages utilisés dans le couplage. La figure de gauche représente le maillage de réservoir pour la simulation de champ, et la figure de droite représente le maillage au voisinage de puits dans le modèle abord de puits. Les mailles aux bords (en gris) dans le modèle abord de puits coïncident avec les mailles de même couleur dans le maillage de réservoir. La croix indique l'emplacement du puits. - Les pas de temps utilisés dans le modèle abord de puits sont généralement beaucoup plus petits que ceux du modèle de réservoir. Le modèle de réservoir est principalement utilisé pour simuler les écoulements dans le réservoir pris dans son intégralité.
- Le temps T0 est le temps auquel commence le couplage. L'algorithme de couplage dans un cadre général comporte les étapes suivantes, illustrées sur la
figure 2 : - 3a. On initialise les modèles.
- on initialise le modèle de réservoir (RINIT) en affectant aux mailles du maillage de réservoir des valeurs de porosité, perméabilité, pression et saturations en fluide. L'initialisation comporte également la définition de conditions aux limites du modèle de réservoir. Ces conditions peuvent être définies par un flux nul (pas d'échange vers l'extérieur du domaine) ou par un flux ou une pression imposé sur les bords extérieurs des mailles de bord du maillage du modèle de réservoir (échange avec l'extérieur). Les conditions opératoires dans ces puits, telles que le débit ou la pression de fond, sont imposées sous forme d'un historique d'injection pour les injecteurs, et d'un historique de production dans les producteurs.
- on initialise le modèle d'abord de puits (NWINIT) en affectant aux mailles du maillage d'abord de puits des valeurs de porosité, perméabilité, pression et saturations en fluide. Cette affectation est réalisée en utilisant des techniques de mise à l'échelle des résultats du modèle de réservoir. Ces techniques étant connues des spécialistes. L'initialisation comporte également la définition de conditions aux limites du modèle d'abord de puits. Ces conditions peuvent également être définies au moyen des résultats du modèle de réservoir.
- 3b. On définit au moins un pas de temps, noté ΔT, pour échanger des données dynamiques entre le modèle de réservoir et le modèle d'abord de puits, pendant la modélisation sur l'intervalle de temps D.
Ce pas de temps ΔT peut être choisi en fonction du pas de temps ΔTR du simulateur d'écoulement du modèle de réservoir, et le pas de temps ΔTNW du simulateur d'écoulement du modèle d'abord de puits (ΔTR>ΔTNW).
Théoriquement, ΔT doit être le plus petit possible pour assurer la convergence des solutions dans les deux modèles. Mais généralement l'utilisation du pas de temps utilisé pour la simulation du modèle de réservoir est suffisante. Mais, du point de vue pratique, on a parfois besoin de réaliser une simulation d'abord de puits de façon autonome plus longtemps. Ceci se traduit par une réduction de la fréquence de couplage. C'est pourquoi, selon la méthode, le pas de temps ΔT pour les échanges de données entre le modèle de réservoir et le modèle d'abord de puits est un paramètre ajustable.
Selon un mode de réalisation, le pas de temps ΔT peut varier sur l'intervalle de temps D. On peut par exemple utiliser un premier pas de temps entre T0 et Ti , et un second pas de temps entre Ti et Tn. Un exemple d'une telle application est illustré ci-après. Sur lafigure 2 , on note RSIM(T1 ) une simulation du simulateur de réservoir effectuée entre T0 et T1, et NWSIM(T1 ) une simulation du simulateur d'abord de puits effectuée entre T0 et T1 . - 3c. On réalise une simulation d'écoulements avec le modèle de réservoir entre le temps T0 et le temps T1 = T0 + ΔT.
Les résultats de cette simulation sont :- la pression et les saturations en fluides à la fin du pas de temps dans chaque maille du maillage de réservoir, en particulier dans les mailles qui sont communes avec les mailles du maillage d'abord de puits et qui vont servir de conditions aux limites du modèle d'abord de puits ;
- les débits en fluide (eau, huile, gaz) et les pressions dans les puits d'injection et de production.
- 3d. On met à jour les conditions aux limites (MAJCL) du modèle d'abord de puits en utilisant les résultats de la simulation d'écoulements réalisée avec le modèle de réservoir entre T0 et T1 (étape 3c).
Les conditions aux limites, sont les valeurs de données dynamiques, telles que la pression ou le flux, les saturations, dans les mailles qui forment les limites du maillage de réservoir ou d'abord de puits. Selon un exemple, les conditions aux limites sont définies par un flux nul sur tous les bords du maillage des abords de puits, et une porosité très grande (1000000, par exemple) dans toutes ces mailles.
Ainsi, au cours de cette étape, on utilise les résultats du simulateur d'écoulement du modèle de réservoir, pour déterminer des valeurs que l'on impose comme conditions aux limites pour le simulateur d'écoulement du modèle d'abord de puits au temps T0.
Les conditions aux limites peuvent être calculées à chaque pas de temps du modèle d'abord de puits en interpolant linéairement les résultats de simulation du modèle de réservoir entre T0 et T1. - 3e. On réalise une simulation d'écoulements au voisinage du puits avec le modèle d'abords de puits entre le temps T0 et le temps T1, avec les conditions aux limites mises à jour à l'étape 3d.
Les résultats de cette simulation sont, au moins :- la pression et les saturations en fluides à la fin du pas de temps dans chaque maille du modèle d'abords de puits;
- les débits en fluide (eau, huile, gaz) et les pressions dans le puits d'injection ou de production selon du type de puits modélisé dans le modèle d'abords de puits.
Ces résultats permettent de déterminer un Indice de Productivité numérique (IP). - 3f. La connexion entre la pression dans les mailles traversées par un puits et les pressions dans le puits lui-même est réalisée à l'aide d'un Indice de Productivité numérique (IP). Les formules de Peaceman sont en général utilisées pour calculer cet indice. On met alors à jours les indices de productivité numériques (MAJIP) du modèle de réservoir, en utilisant les résultats de la simulation d'écoulements réalisée avec le modèle d'abord de puits entre T0 et T1 . En fait, si à la fin de la simulation, au temps T1, les résultats simulés au puits avec le modèle d'abord de puits et avec le modèle de réservoir ne sont pas les mêmes, les indices de productivité numériques dans le modèle de réservoir sont modifiés pour ajuster les résultats de la simulation du modèle de réservoir à ceux du modèle d'abord de puits.
- 3g. On répète les étapes 3c (éventuellement 3b) à 3f, avec un nouvel intervalle de temps (de T1 à T2 , puis de T2 à T3, ..., puis de Tn-1 à Tn )
-
- i
- : numéro de maille de puits dans le maillage (de réservoir ou d'abord de puits)
- p
- : phase du fluide. Les phases p peuvent être l'eau, l'huile ou le gaz
- Qp,i
- : débit de la phase p dans la maille de puits i du maillage (de réservoir ou d'abord de puits)
- λp,i
- : mobilité de la phase p dans la maille de puits i du maillage (de réservoir ou d'abord de puits); λp,i dépend essentiellement de la perméabilité relative et de la viscosité de la phase p
- IPi
- : indice de productivité numériques dans la maille de puits i du maillage (de réservoir ou d'abord de puits)
- Pp,i
- : pression de la phase p dans la maille de puits i du maillage (de réservoir ou d'abord de puits)
- Pwf,i
- : pression dans le puits, au fond, au niveau du réservoir dans la maille de puits i du maillage (de réservoir ou d'abord de puits)
- L'indice de productivité numérique IP prend en compte : l'effet géométrique de la maille de puits i du maillage, la perméabilité du milieu poreux dans la maille du puits et un coefficient de skin. Un coefficient de skin est un coefficient, bien connu de l'homme de l'art, utilisé pour représenter l'endommagement d'un puits dans une maille.
-
- i
- : numéro de maille de puits dans le maillage de réservoir
- j
- : numéro de maille de puits dans le maillage d'abord de puits
- Wi
- : ensemble des mailles de puits du maillage d'abord de puits correspondant à un raffinement de la maille de puits i du maillage de réservoir
- p
- : phase du fluide. Les phases p peuvent être l'eau (w), l'huile (o) ou le gaz (g)
- IPr,i
- : indice de productivité numériques dans la maille de puits i du maillage de réservoir, et utilisé dans le modèle de réservoir
- Pnw,p,j
- : pression de la phase p dans la maille de puits j du maillage d'abord de puits, calculée avec le modèle abord de puits
- Pr,p,i
- : pression de la phase p dans la maille de puits i du maillage de réservoir, calculée avec le modèle de réservoir
- Pwf,j
- : pression dans le puits au niveau du réservoir dans la maille de puits j du maillage d'abord de puits
- IPnw,j
- : indice de productivité numériques dans la maille de puits j du maillage d'abord de puits, et utilisé dans le modèle d'abord de puits
- Les variables IPi, Pnw,p,j, Pr,p,i et Pwf,j sont fonction du temps T.
-
- Qnw,i : débit du fluide (unique phase) calculé avec le modèle d'abord de puits dans la section correspondant à la partie du puits dans la maille de puits i du maillage de réservoir
- Qr,i : débit du fluide (unique phase) calculé avec le modèle de réservoir dans la même section, correspondant à la partie du puits dans la maille de puits i du maillage de réservoir
- IPr,i (T 1) et IPr,i (T 0) sont les indices de productivité numérique aux temps T1 et T0 respectivement, c'est-à-dire, avant et après la mise à jour.
- Cette formule montre clairement que la mise à jour de l'indice de productivité numérique correspond à la correction du débit de fluide du modèle de réservoir par rapport au débit de fluide du modèle abord de puits : si les deux modèles donnent le même résultat en terme de débit, alors
- En testant divers scénarios, caractérisés par exemple par diverses implantations respectives des puits injecteurs et producteurs, et en simulant la production d'hydrocarbures pour chacun d'eux selon l'étape 3, on peut sélectionner le scénario optimal. Le scénario optimal est le scénario permettant d'obtenir une production optimale du gisement dans la cadre de la production d'un réservoir, ou le scénario permettant d'obtenir l'injectivité optimale dans le gisement dans le cadre d'injection de fluide dans le réservoir (injection d'eau pour une production améliorée, ou injection de gaz acides).
- Pour tester divers scénarios d'exploitation, on modifie (ΔSCE) le scénario sélectionné à l'étape 1, par exemple en modifiant l'emplacement d'un puits.
- On optimise alors l'exploitation du gisement, en mettant en oeuvre, sur le champ, le scénario de production ainsi sélectionné.
- Selon l'invention, il est tout fait possible de coupler un modèle de réservoir avec plusieurs modèles d'abord de puits.
- Selon un mode de réalisation particulier de l'invention, l'étape 2, au cours de laquelle on construit les maillages, est modifiée.
- La simulation, réalisée au moyen du modèle de réservoir dans l'étape 3c, fournit des propriétés dynamiques des fluides comme la pression ou les saturations dans la période allant de T0 à T1 sur toutes les mailles grossières. Cependant, la détermination des conditions aux limites dans l'étape 3b nécessite l'interpolation de la pression ou du flux sur les bords du modèle abord de puits. Pour réduire les erreurs dans l'interpolation, nous pouvons contraindre, lors de la génération de maillage, les mailles de bord du modèle abord de puits pour qu'elles coïncident avec les interfaces des mailles du modèle de réservoir. De plus, les mailles de bord dans le modèle abord de puits sont aussi contraintes afin qu'elles coïncident avec des mailles du modèle de réservoir (
figure 3 ). De cette façon, le transfert des données dynamiques du modèle de réservoir vers le modèle abord de puits est direct sur ces mailles. Dans le modèle abord de puits lui-même, les conditions aux limites sont de flux nul. Afin de maintenir les propriétés dynamiques aux bords du modèle, des porosités de très grande valeur - (1000000, par exemple) sont affectées aux mailles de bord. Ce type de conditions aux limites est consistant avec la plupart des modèles d'écoulement, et son implémentation est facile.
- Pour certains problèmes, les changements d'écoulement autour du puits sont liés aux écoulements multiphasiques. Dans ce cas-là, nous pouvons aussi faire la mise à jour des indices de productivité numérique par phase. Pour cela, on reformule la relation de pression/débit, en introduisant un coefficient, appelé multiplicateur de l'indice de productivité :
- Mp,i est le multiplicateur de l'indice de productivité pour la phase p dans la maille de puits i.
-
- Qr,p,i (T 1)
- : débit de la phase p calculé par le modèle de réservoir dans la maille de puits i du maillage de réservoir au temps T1
- Qnw,p,i (T 1)
- : débit de la phase p calculé par le modèle d'abord de puits dans le même secteur du puits (voir l'ensemble Wi ) au temps T1
- Mp,i (T 0)
- : multiplicateur d'indice de productivité numérique pour la phase p dans le modèle de réservoir aux temps T0 (avant mise à jour du modèle)
- Mp,i (T 1)
- : multiplicateur d'indice de productivité numérique pour la phase p dans le modèle de réservoir aux temps T1 (après mise à jour du modèle)
- La méthode de couplage selon l'invention peut être utilisée pour modéliser différents phénomènes détaillés autour du puits, comme par exemple, l'endommagement par fluide de forage ou de complétion, la stimulation acide, l'écoulement non-Darcéen autour du puits, le problème de gaz à condensat, le dépôt d'asphaltène, l'endommagement par l'injection de CO2, la prévention de venue d'eau ou de gaz, la venue de sable, les dépôts minéraux, l'impact des complétions, etc. Ici, nous présentons en particulier un exemple d'application pour l'endommagement de la formation pétrolifère par le fluide de forage lors du forage du puits, et un exemple d'application pour la prévention de venue d'eau lorsqu'un puits en production produit une quantité importante d'eau, et que l'on cherche à réduire cette production d'eau.
- Pour simplifier encore la méthode de couplage, nous faisons la mise à jour des données en utilisant les valeurs au temps Tn , au lieu de l'interpolation linéaire en temps entre Tn et T n+ 1, pour la simulation du modèle abord de puits dans la période allant de Tn à T n+ 1. Ce choix est intéressant, car il permet de faire des simulations en parallèle sur différentes machines pour le modèle de réservoir et le modèle abord de puits.
- Un modèle de réservoir standard est utilisé pour la simulation du champ. Le modèle abord de puits développé par DING, Y. and RENARD, G. : "Evaluation of Horizontal Well Performance after Drilling Induced Formation Damage" J. of Energy Resources Technology, Vol. 127, Sept., 2005, est utilisé pour simuler l'endommagement de la formation par le forage. Il faut noter que la physique avancée de l'endommagement n'est pas modélisée dans la simulation du champ avec le modèle de réservoir.
- Un réservoir de taille 1000m x 1000m x 10m est considéré. Un maillage cartésien avec 20 mailles dans la direction x, 20 mailles dans la direction y et 1 maille dans la direction z est utilisé pour la simulation du champ (
figure 6 ). Les tailles de mailles sont donc 50m x 50m x 10m. La pression initiale de réservoir est de 200 bars. Un puits producteur doit être foré dans le block (15, 15, 1). Il est représenté par un cercle noir sur lafigure 6 . L'endommagement de ce puits par le fluide de forage est étudié avec la méthode selon l'invention. - Le réservoir est homogène de perméabilité 200 mD et de porosité 0.15. Les conditions aux limites de ce réservoir sont des flux nuls, sauf sur le bord Γx- (
figure 6 ), où la pression est constante (200 bars). - Pour obtenir la solution de référence, le maillage est raffiné autour du puits (
figure 7 ). Un modèle spécifique, qui prend en compte la physique avancée de l'endommagement, est utilisé sur ce maillage pour simuler la solution de référence. Comme l'endommagement par le fluide de forage est généralement limité à quelques centimètres ou quelques dizaines de centimètres autour de puits, nous avons besoin de très petites mailles dans la zone raffinée (Tableau 1). Le diamètre du puits est de 21.6 cm. Pour que le puits soit inclus dans une maille, la taille de la maille puits est de 22 cm. Les autres mailles autour du puits sont beaucoup plus petites avec une taille de 2 cm. Les maillages utilisés pour le couplage sont illustrés sur lesfigures 8A et 8B . Le maillage du modèle abord de puits (figure 8B ) correspond à la zone raffinée et aux mailles autour dans le maillage de référence. Les mailles aux bords du modèle abord de puits coïncident avec des mailles du modèle de réservoir.Tableau 1-Tailles des mailles autour du puits Taille de maille dans la direction x (m) Taille de maille dans la direction y (m) 50 42.7 30 20 16 8 4 2 1 0.51 0.3 50 42.7 30 20 16 8 4 2 1 0.51 0.3 0.16 0.08 0.04 0.02 0.02 0.02 0.02
0.020.16 0.08 0.04 0.02 0.02 0.02 0.02
0.020.22 0.22 0.02 0.02 0.02 0.02 0.02 0.04 0.08
0.160.02 0.02 0.02 0.02 0.02 0.04 0.08
0.160.3 0.51 1 2 4 8 16 20 30 42.7 50 0.3 0.51 1 2 4 8 16 20 30 42.7 50 - Il est supposé que le réservoir est épais, et que ce modèle correspond seulement à la première couche du réservoir. Le temps de contact entre le fluide de forage et le réservoir est de 2 jours. La pression au cours du forage au fond du puits est de 250 bars. La perméabilité et l'épaisseur du « cake » externe formé par la boue de forage sont égales à 0.001 mD et 0.2 cm. L'épaisseur du « cake » interne est de 2 cm avec une perméabilité moyenne réduite à 20 mD durant la période de forage et de 40 mD dans la période de production. La viscosité du fluide de forage est de 30 cPo. L'hystérésis de la perméabilité relative entre les périodes de forage et de mise en production est présenté à la
figure 9 . Une saturation en eau irréductible de 30% liée au filtrat (fluide de forage) qui va envahir la formation pendant la phase de forage va rester bloquée dans le milieu poreux quand le puits va être remis en production. - Les volumes d'invasion de fluide de forage sont comparés à la
figure 10 pour la simulation avec la méthode de couplage et la solution de référence obtenue en utilisant le maillage avec le raffinement local (figure 7 ). Les pas de temps pour la mise à jour des données dans le couplage sont présentés dans le Tableau 2. Lafigure 10 montre que le volume d'invasion de fluide est correctement simulé avec la méthode de couplage. Le petit écart entre la solution de couplage et la solution de référence dans la période entre 0.1 et 0.3 jour peut être amélioré en utilisant des petits pas d'itération en temps pour échanger les données dans le couplage.Tableau 2-Pas de temps pour la mise à jour des données dans le couplage Période (jour) Pas de temps (jour) 0 - 0.01 0.001 0.01 - 0.1 0.01 0.1 - 3 0.1 3 - 10 1 10 - 200 10 - Après les 2 jours de forage, le puits est fermé pendant 1 jour pour la mise en place de sa complétion, puis il est mis en production. Le couplage est effectué jusqu'au 10ème jour. Au-delà de 10 jours, l'effet de l'endommagement autour de puits devient stable et les IPs numériques dans le modèle de réservoir ne changent quasiment plus. Nous n'avons plus besoin de couplage pour continuer la simulation du champ avec le modèle de réservoir. La courbe de production d'huile simulée par le modèle de réservoir, qui est couplé avec le modèle abord de puits pendant les 10 premiers jours, est présentée à la
figure 11 . Cette courbe est très proche de la solution de référence. - Si l'endommagement n'est pas pris en compte ou si seulement la présence des « cakes » est considérée dans la simulation, les résultats sont très imprécis avec des erreurs de plus de 20% (
figure 11 ). Prendre en compte des phénomènes autour de puits tel que l'endommagement par le fluide de forage est important pour la gestion de réservoir, et la méthode de couplage proposée convient parfaitement pour simuler ce type de problèmes. - Dans la procédure de prévention de venue d'eau, une solution de polymère est injectée dans un puits producteur pour une courte durée dans le but de réduire la grande quantité d'eau produite en même temps que l'huile. Une partie du polymère est absorbée sur la roche, et une autre partie est dispersée dans l'eau. Le polymère injecté a pour effet de réduire la mobilité de la phase eau par l'augmentation de sa viscosité et par la diminution de la perméabilité relative de cette phase. Donc, dans la méthode de couplage, l'approche la plus adaptée est de mettre à jour le multiplicateur d'IP numérique pour la phase eau.
- Un réservoir de 1000m x 1000m x 25m est considéré comme exemple. Un maillage cartésien avec 20 mailles dans la direction x, 20 mailles dans la direction y et 5 mailles dans la direction z est utilisé pour la simulation du champ. La taille des mailles est de 50m x 50m x 5m. Le réservoir est hétérogène. La perméabilité est présentée à la
figure 12 . Le rapport des perméabilités dans les directions verticale et horizontale est de 0.1. La pression initiale du réservoir est de 200 bars. - Il y a un puits injecteur (INJ) et un puits producteur (PROD) comme montré à la
figure 12 . La pression au puits injecteur est imposée à 300 bars, et la pression au puits producteur est contrainte à 150 bars pendant la production. Après production pendant 1000 jours, le water-cut (débit d'eau rapporté au débit total) du puits producteur atteint 85%. La procédure de prévention de venue d'eau est alors appliquée pour réduire la quantité d'eau produite. Une solution de polymère avec une concentration de 2500 ppm est injectée dans le producteur avec une pression de fond de 300 bars pendant 2 jours. Ensuite, le puits est remis en production. Cette procédure de prévention de venue d'eau est simulée avec la méthode selon l'invention. - Afin d'avoir une solution de référence, un raffinement local autour du puits producteur est utilisé (
figure 13 ). La taille des mailles autour du puits est de 0.617 m dans la direction x. Le maillage pour le couplage est présenté à lafigure 14 . Des mailles aux bords du modèle abord de puits coïncident avec des mailles dans le modèle de réservoir. La physique du polymère peut être considérée dans les deux modèles (modèle abord de puits et modèle de réservoir).Tableau 3 - Pas de temps pendant le couplage Période (jour) Pas de temps (jour) 0 - 950 - 950-970 2 970 - 1000 28 1000 - 1000.1 0.01 1000.1 - 1005 0.1 1005 - 1030 1 1030 - 1100 2 1100 - 3000 - - Le couplage commence à 950 jours et se termine à 1100 jours, soit une période de 150 jours au total. Les pas de temps pour les échanges des données dans la méthode de couplage sont présentés dans le Tableau 3. Dans les premiers 50 jours (de 950 à 1000 jours) de couplage, il n'y a pas de polymère injecté. Cette période est uniquement utilisée pour assurer une bonne initialisation du modèle abord de puits. Les IPs numériques globaux sont mis à jour au début du couplage (de 950 à 970 jours) pour prendre en compte les effets des maillages entre le modèle de réservoir et le modèle abord de puits. Pendant la période d'injection de polymère (entre 1000 et 1002 jours), les IPs numériques globaux sont encore recalculés pour intégrer l'effet induit par le polymère injecté (on pourrait aussi mettre à jour les multiplicateurs d'IP numérique pour la phase eau). Mais, quand le puits est remis en production (à 1003 jours), les multiplicateurs d'IP numérique pour la phase d'eau sont mis à jour.
- La
figure 15 compare les débits d'injection de polymère dans le puits pour les différentes simulations: la solution de référence, la simulation sur le maillage de réservoir avec couplage, la simulation directe sur le maillage de réservoir sans couplage et la simulation avec le modèle abord de puits (avec couplage). Lesfigures 16A à 16E montrent les mêmes comparaisons couche par couche. Pour la simulation directe avec le maillage de réservoir sans couplage, le volume de polymère injecté est largement surestimé. Quand la simulation du maillage grossier est couplée avec le modèle abord de puits, les résultats sont significativement améliorés. Au début du couplage, le débit d'injection est grand, mais il est rapidement corrigé par la mise à jour de l'IP dû au couplage. Si l'on veut avoir plus de précision sur le débit d'injection de polymère, il suffit de se référer aux résultats de simulation avec le modèle abord de puits. Avec ce modèle, le volume injecté et la distribution du polymère autour du puits sont tous les deux correctement simulés. - Les
figures 17, 18 et19 présentent les courbes de débits d'huile, d'eau et de water-cut pour le modèle de réservoir avec couplage, le modèle de réservoir sans couplage et la solution de référence. Les résultats du modèle de réservoir avec couplage sont globalement satisfaisants. Lafigure 20 présente la carte de saturation d'eau à la fin du couplage (1100 jours), et lafigure 21 montre la carte de pression à 1100 jours. Comparé aux solutions de référence, le couplage donne des résultats globalement satisfaisants.
Claims (9)
- Méthode mise en œuvre par ordinateur pour modéliser des écoulements de fluides au sein d'un milieu poreux souterrain traversé par au moins un puits en vue de l'exploitation dudit milieu réalisée en injectant un desdits fluides dans ledit milieu via au moins un desdits puits, et/ou en produisant un desdits fluides présent dans ledit milieu au moyen d'au moins un desdits puits, dans laquelle on utilise un premier simulateur d'écoulement permettant de simuler l'écoulement des fluides au sein du milieu poreux à partir d'indices de productivité numériques reliant des pressions de fluides à des débits de fluides, et on utilise un second simulateur d'écoulement pour simuler l'écoulement des fluides aux abords du puits à partir de conditions aux limites, caractérisé en ce que :a- on simule les écoulements de fluides au sein du milieu au moyen du premier simulateur sur un intervalle de temps défini entre des temps T0 et T1, et on déduit des conditions aux limites mises à jour pour le second simulateur par une interpolation linéaire des résultats du premier simulateur entre les temps T0 et T1 ;b- on simule les écoulements de fluides aux abords du puits au moyen du second simulateur sur le même intervalle de temps, en utilisant les conditions aux limites mises à jour, et on en déduit des indices de productivité numériques mis à jour pour le premier simulateur en comparant des débits calculés par le premier simulateur et des débits calculés par le second simulateur; etc- on modélise les écoulements de fluides au sein du milieu poreux pendant une période de temps entre T0 et Tn, où Tn >T1 , en réitérant les étapes a et b, pour des intervalles de temps successifs compris entre T0 et Tn.
- Méthode selon la revendication 1, dans laquelle chaque intervalle de temps successif a une longueur qui est fonction d'un pas de temps de calcul du premier simulateur d'écoulement et d'un pas de temps du second simulateur d'écoulement.
- Méthode selon la revendication 1, dans laquelle chaque intervalle de temps successifs a une longueur égale à un pas de temps du premier simulateur d'écoulement.
- Méthode selon l'une des revendications précédentes, dans laquelle on simule les écoulements de fluides au sein du milieu au moyen du premier simulateur sur un premier maillage discrétisant le milieu poreux en un ensemble de mailles, et on simule les écoulements de fluides aux abords du puits au moyen du second simulateur sur un second maillage discrétisant le puits et ses abords en un ensemble de mailles, ledit second maillage étant généré en contraignant des mailles situées sur le bord dudit second maillage, de façon à ce que leurs interfaces coïncident avec les interfaces des mailles dudit premier maillage.
- Méthode selon l'une des revendications précédentes, dans laquelle on modélise des écoulements multiphasiques, et on met à jour des multiplicateurs d'indices de productivité numériques, au lieu des indices de productivité numériques eux-mêmes, pour chaque phase, en comparant des débits par phase calculés par le premier simulateur et des débits par phase calculés par le second simulateur.
- Méthode, selon l'une des revendications 1 à 5, pour exploiter un réservoir poreux souterrain au moyen d'au moins un puits traversant la réservoir, au moins un fluide circulant entre le réservoir et le puits, dans laquelle on acquiert des données relatives à la géométrie dudit réservoir poreux, à partir desquelles on construit une discrétisation du réservoir en un ensemble de mailles, appelé maillage de réservoir, et on construit une discrétisation du puits et de ses abords en un ensemble de mailles, appelé maillage d'abord de puits, caractérisé en ce que l'on réalise les étapes suivantes :a- on choisit un scénario d'exploitation du réservoir poreux ;b- on associe au maillage de réservoir un premier simulateur d'écoulement permettant de simuler l'écoulement des fluides au sein du réservoir, à partir d'au moins les données suivantes : le scénario de production, des données d'entrées relatives au fluide et au réservoir, des indices de productivité numériques permettant de relier des pressions à des débits, des conditions aux limites ;c- on associe au maillage d'abord de puits un second simulateur d'écoulement pour simuler l'écoulement des fluides aux abords du puits, à partir d'au moins les données suivantes : des données d'entrées relatives au fluide et au réservoir, des conditions aux limites ;d- on modélise les écoulements de fluides au sein du réservoir et aux abords du puits, au moyen de la méthode selon l'une des revendications 1 à 5 ; ete- on modifie le scénario d'exploitation et l'on répète l'étape d, jusqu'à obtenir un scénario d'exploitation optimal.
- Méthode selon la revendication 6, dans laquelle on tient compte d'un endommagement du puits par un fluide de forage en modélisant une invasion du réservoir poreux par ledit fluide de forage aux étapes d et e.
- Méthode selon la revendication 6, dans laquelle le scénario d'exploitation comporte une injection d'une solution de polymère par le puits, et l'on modélise les écoulements pour prévenir une venue d'eau.
- Méthode selon la revendication 6, dans laquelle le scénario d'exploitation comporte une injection d'une solution acide dans le puits, et l'on modélise les écoulements pour évaluer l'impact d'une stimulation acide.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0902533A FR2945879B1 (fr) | 2009-05-20 | 2009-05-20 | Methode d'exploitation de milieu poreux au moyen d'une modelisation d'ecoulements de fluide |
Publications (2)
Publication Number | Publication Date |
---|---|
EP2253797A1 EP2253797A1 (fr) | 2010-11-24 |
EP2253797B1 true EP2253797B1 (fr) | 2020-02-19 |
Family
ID=41426264
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP10290174.1A Active EP2253797B1 (fr) | 2009-05-20 | 2010-04-01 | Méthode d'exploitation de milieu poreux au moyen d'une modélisation d'écoulements de fluide |
Country Status (4)
Country | Link |
---|---|
US (1) | US8694297B2 (fr) |
EP (1) | EP2253797B1 (fr) |
CA (1) | CA2704060C (fr) |
FR (1) | FR2945879B1 (fr) |
Families Citing this family (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2702965C (fr) | 2007-12-13 | 2014-04-01 | Exxonmobil Upstream Research Company | Partitionnement parallele adaptatif de donnees sur une simulation de reservoir utilisant une grille non structuree |
US9581723B2 (en) | 2008-04-10 | 2017-02-28 | Schlumberger Technology Corporation | Method for characterizing a geological formation traversed by a borehole |
EP2277065B1 (fr) * | 2008-05-03 | 2013-09-18 | Saudi Arabian Oil Company | Système, produit de programme et procédés associés pour effectuer une estimation de pression de gisement en temps réel et automatisée permettant des stratégies d'injection et de production optimisées |
US9134457B2 (en) * | 2009-04-08 | 2015-09-15 | Schlumberger Technology Corporation | Multiscale digital rock modeling for reservoir simulation |
IN2012DN05167A (fr) | 2010-02-12 | 2015-10-23 | Exxonmobil Upstream Res Co | |
WO2011112221A1 (fr) | 2010-03-12 | 2011-09-15 | Exxonmobil Upstream Research Company | Groupage dynamique d'objets de domaine au moyen de groupes intelligents |
US8583410B2 (en) * | 2010-05-28 | 2013-11-12 | Ingrain, Inc. | Method for obtaining consistent and integrated physical properties of porous media |
US9540911B2 (en) * | 2010-06-24 | 2017-01-10 | Schlumberger Technology Corporation | Control of multiple tubing string well systems |
US10318663B2 (en) | 2011-01-26 | 2019-06-11 | Exxonmobil Upstream Research Company | Method of reservoir compartment analysis using topological structure in 3D earth model |
WO2013043531A1 (fr) * | 2011-09-20 | 2013-03-28 | Bp Corporation North America Inc. | Génération automatique d'un affinement local de maillage au niveau de fractures hydrauliques pour simuler des réservoirs compacts de gaz |
FR2989200B1 (fr) | 2012-04-10 | 2020-07-17 | IFP Energies Nouvelles | Procede de selection des positions de puits a forer pour l'exploitation d'un gisement petrolier |
FR2997721B1 (fr) * | 2012-11-08 | 2015-05-15 | Storengy | Radonip : nouvelle methodologie de determination des courbes de productivite des puits d'exploitation de stockages et gisements de fluides compressibles |
FR3002270B1 (fr) * | 2013-02-21 | 2015-02-27 | Ifp Energies Now | Procede d'exploitation d'un reservoir geologique au moyen d'un modele de reservoir cale et coherent vis a vis des proprietes d'ecoulement |
GB2512372B (en) * | 2013-03-28 | 2020-07-29 | Total Sa | Method of modelling a subsurface volume |
AU2014278645B2 (en) | 2013-06-10 | 2016-07-28 | Exxonmobil Upstream Research Company | Interactively planning a well site |
US10480314B2 (en) | 2013-07-26 | 2019-11-19 | Schlumberger Technology Corporation | Well treatment |
US9864098B2 (en) | 2013-09-30 | 2018-01-09 | Exxonmobil Upstream Research Company | Method and system of interactive drill center and well planning evaluation and optimization |
US11414975B2 (en) * | 2014-07-14 | 2022-08-16 | Saudi Arabian Oil Company | Quantifying well productivity and near wellbore flow conditions in gas reservoirs |
US9816366B2 (en) * | 2014-07-14 | 2017-11-14 | Saudi Arabian Oil Company | Methods, systems, and computer medium having computer programs stored thereon to optimize reservoir management decisions |
CN104500040B (zh) * | 2014-10-16 | 2017-06-06 | 西南石油大学 | 水平井酸化过程中井筒多段流体移动界面跟踪方法 |
CN105626007B (zh) * | 2014-11-07 | 2018-06-15 | 中国石油化工股份有限公司 | 基于岩心尺度油藏中不同部位过水倍数计算方法 |
US10280722B2 (en) | 2015-06-02 | 2019-05-07 | Baker Hughes, A Ge Company, Llc | System and method for real-time monitoring and estimation of intelligent well system production performance |
CA2992714A1 (fr) * | 2015-08-21 | 2017-03-02 | Halliburton Energy Services, Inc. | Procede et protocole de modelisation precise de formations de champ proche dans des simulations de puits de forage |
CN111094694B (zh) * | 2017-09-08 | 2022-08-12 | 罗克斯软件解决方案有限公司 | 井裂缝建模 |
EP3724793A4 (fr) * | 2017-12-14 | 2021-10-20 | Services Pétroliers Schlumberger | Système et procédé de simulation de modèles de réservoir |
CN109882149B (zh) * | 2018-01-29 | 2023-04-07 | 西南石油大学 | 一种模拟缝洞型碳酸盐岩凝析气藏生产动态的实验装置及方法 |
WO2019226149A1 (fr) | 2018-05-21 | 2019-11-28 | Newpark Drilling Fluids Llc | Système de simulation de conditions de forage en fond de trou in situ et de test d'échantillons carottés |
CN110778312B (zh) * | 2019-10-09 | 2022-08-30 | 东北石油大学 | 一种模拟气藏边底水侵入的模型以及求取水侵系数的方法 |
CN110924908B (zh) * | 2019-11-08 | 2021-10-08 | 中国石油大学(华东) | 一种水驱油藏注采参数确定方法及计算机可读存储介质 |
WO2021118714A1 (fr) * | 2019-12-11 | 2021-06-17 | Exxonmobil Upstream Research Company | Méthodologie de semi-élimination permettant de simuler des caractéristiques d'écoulement élevées dans un réservoir |
US11586790B2 (en) | 2020-05-06 | 2023-02-21 | Saudi Arabian Oil Company | Determining hydrocarbon production sweet spots |
US20220075803A1 (en) * | 2020-09-04 | 2022-03-10 | Saudi Arabian Oil Company | Graph framework (database methods) to analyze trillion cell reservoir and basin simulation results |
CN112069690B (zh) * | 2020-09-11 | 2024-03-08 | 中海石油(中国)有限公司 | 一种深水断块油藏长水平井多级油嘴测试产能的评价方法 |
CN113065261B (zh) * | 2021-04-25 | 2024-01-02 | 中国长江三峡集团有限公司 | 基于水热耦合模拟的地热资源回收率的评价方法 |
US12020192B2 (en) * | 2021-09-24 | 2024-06-25 | Saudi Arabian Oil Company | Estimating well downtime factor in field modeling |
CN114001654B (zh) * | 2021-11-01 | 2024-03-26 | 北京卫星制造厂有限公司 | 工件端面位姿评价方法 |
US11613957B1 (en) | 2022-01-28 | 2023-03-28 | Saudi Arabian Oil Company | Method and system for high shut-in pressure wells |
US12024985B2 (en) | 2022-03-24 | 2024-07-02 | Saudi Arabian Oil Company | Selective inflow control device, system, and method |
CN114839130B (zh) * | 2022-05-12 | 2024-07-19 | 西南石油大学 | 一种高温高压大尺度剖面模型束缚水实验条件的建立方法 |
CN117828732B (zh) * | 2024-01-02 | 2024-05-31 | 中国恩菲工程技术有限公司 | 基于数字孪生的边坡稳定性确定方法及系统、介质、终端 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2801710B1 (fr) * | 1999-11-29 | 2002-05-03 | Inst Francais Du Petrole | Methode pour generer un maillage hybride permettant de modeliser une formation heterogene traversee par un ou plusieurs puits |
FR2886743B1 (fr) * | 2005-06-02 | 2007-07-27 | Inst Francais Du Petrole | Methode pour simuler les ecoulements de fluides au sein d'un reservoir a l'aide d'une discretisation de type chimere |
CA2628721A1 (fr) * | 2005-11-22 | 2007-05-31 | Exxonmobil Upstream Research Company | Systeme et procede de simulation |
EP2038809B1 (fr) * | 2006-06-26 | 2019-07-03 | Exxonmobil Upstream Research Company | Procédé de comparaison et de réallocation de production |
US20080319726A1 (en) * | 2007-06-19 | 2008-12-25 | Schlumberger Technology Corporation | System and method for performing oilfield simulation operations |
US8775141B2 (en) * | 2007-07-02 | 2014-07-08 | Schlumberger Technology Corporation | System and method for performing oilfield simulation operations |
CA2696638C (fr) * | 2010-03-16 | 2012-08-07 | Exxonmobil Upstream Research Company | Utilisation d'une emulsion dont la phase externe est un solvant pour la recuperation in situ de petrole |
-
2009
- 2009-05-20 FR FR0902533A patent/FR2945879B1/fr not_active Expired - Fee Related
-
2010
- 2010-04-01 EP EP10290174.1A patent/EP2253797B1/fr active Active
- 2010-04-27 US US12/767,866 patent/US8694297B2/en not_active Expired - Fee Related
- 2010-05-17 CA CA2704060A patent/CA2704060C/fr not_active Expired - Fee Related
Non-Patent Citations (1)
Title |
---|
ED NENNIE ET AL: "SPE 110 316 - An Investigation into the need of a dynamic coupled well-reservoir simulator", COPYRIGHT, 1 November 2007 (2007-11-01), pages 11 - 14, XP055423977, Retrieved from the Internet <URL:https://www.onepetro.org/download/conference-paper/SPE-110316-MS?id=conference-paper/SPE-110316-MS> [retrieved on 20171110] * |
Also Published As
Publication number | Publication date |
---|---|
EP2253797A1 (fr) | 2010-11-24 |
FR2945879A1 (fr) | 2010-11-26 |
FR2945879B1 (fr) | 2011-06-24 |
US20100299125A1 (en) | 2010-11-25 |
CA2704060C (fr) | 2018-02-27 |
CA2704060A1 (fr) | 2010-11-20 |
US8694297B2 (en) | 2014-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2253797B1 (fr) | Méthode d'exploitation de milieu poreux au moyen d'une modélisation d'écoulements de fluide | |
EP3144468B1 (fr) | Procédé pour caracteriser le réseau de fractures d'un gisement fracture et procédé pour l'exploiter | |
EP2072752B1 (fr) | Méthode pour optimiser l'exploitation d'un gisement de fluide par prise en compte d'un terme d'échange géologique et transitoire entre blocs matriciels et fractures | |
EP2212722B1 (fr) | Methode de modelisation d'un milieu géologique poreux traversé par un réseau de fractures | |
FR2961844A1 (fr) | Ecoulement multiphase dans un puits de forage et fracture hydraulique connectee | |
EP2813668B1 (fr) | Méthode pour optimiser l'exploitation d'un gisement de fluide par prise en compte d'un terme d'échange géologique et transitoire entre blocs matriciels et fractures | |
WO2002006794A1 (fr) | Methode pour modeliser des deplacements de fluides dans un milieu poreux tenant compte d'effets d'hysteresis | |
EP2375352A1 (fr) | Méthode pour exploiter un réservoir pétrolier à partir d'un calage d'historique optimisé | |
EP2963235B1 (fr) | Procede d'exploitation d'un gisement petrolier a partir d'une technique de positionnement des puits a forer | |
Al-Mudhafar | From coreflooding and scaled physical model experiments to field-scale enhanced oil recovery evaluations: Comprehensive review of the gas-assisted gravity drainage process | |
EP2924232A2 (fr) | Procédé de construction d'un maillage optimise pour la simulation de reservoir dans une formation souterraine | |
FR3005766A1 (fr) | Fractures a segments multiples | |
CA2956285A1 (fr) | Procede pour l'exploitation d'un gisement de fluide traverse par des fractures au moyen d'une simulation d'ecoulement basee sur un flux d'echange et un facteur correctif | |
EP2365359B1 (fr) | Méthode de calage d'historique d'un modèle géologique | |
EP1889999B1 (fr) | Méthode pour optimiser la récupération assistée d'un fluide en place dans un milieu poreux par suivi de front | |
Olalotiti-Lawal et al. | Post-combustion CO2 WAG pilot in a mature field: model calibration and optimization | |
Delshad et al. | A mixed-wet hysteretic relative permeability and capillary pressure model for reservoir simulations | |
Zhao et al. | Morphology and growth of fractures in unconventional reservoirs | |
FR3040509A1 (fr) | ||
EP3192964B1 (fr) | Procede de production d'hydrocarbures comportant un index de productivite des puits sous effet thermique | |
EP4112875B1 (fr) | Procédé pour simuler les écoulements d'un fluide dans une formation souterraine comportant un réseau de fractures | |
FR2950993A1 (fr) | Methode pour modeliser des ecoulements de fluides au sein d'un milieu poreux a partir d'indices de productivite optimises | |
EP3252508B1 (fr) | Procédé d'exploitation des hydrocarbures d'une formation souterraine, au moyen d'une mise à l'échelle optimisée | |
Gamal et al. | An approach for determination of the economically optimal production controlling parameters from water drive oil reservoirs | |
Ghadami et al. | Robust chemical EOR modelling from coreflood to full field scale in a Brown field, offshore |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO SE SI SK SM TR |
|
AX | Request for extension of the european patent |
Extension state: AL BA ME RS |
|
17P | Request for examination filed |
Effective date: 20110524 |
|
RAP1 | Party data changed (applicant data changed or rights of an application transferred) |
Owner name: IFP ENERGIES NOUVELLES |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: EXAMINATION IS IN PROGRESS |
|
17Q | First examination report despatched |
Effective date: 20171204 |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: GRANT OF PATENT IS INTENDED |
|
INTG | Intention to grant announced |
Effective date: 20190919 |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE PATENT HAS BEEN GRANTED |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO SE SI SK SM TR |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D Free format text: NOT ENGLISH |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: EP |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R096 Ref document number: 602010063139 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: AT Ref legal event code: REF Ref document number: 1235173 Country of ref document: AT Kind code of ref document: T Effective date: 20200315 |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: FG4D Free format text: LANGUAGE OF EP DOCUMENT: FRENCH |
|
REG | Reference to a national code |
Ref country code: NO Ref legal event code: T2 Effective date: 20200219 |
|
REG | Reference to a national code |
Ref country code: NL Ref legal event code: MP Effective date: 20200219 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: FI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: FR Payment date: 20200429 Year of fee payment: 11 Ref country code: NO Payment date: 20200422 Year of fee payment: 11 |
|
REG | Reference to a national code |
Ref country code: LT Ref legal event code: MG4D |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IS Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200619 Ref country code: BG Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200519 Ref country code: HR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: SE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: LV Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: GR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200520 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: GB Payment date: 20200429 Year of fee payment: 11 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: NL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: CZ Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: ES Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: LT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: PT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200712 Ref country code: SK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: RO Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: SM Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: EE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: DK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R119 Ref document number: 602010063139 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: AT Ref legal event code: MK05 Ref document number: 1235173 Country of ref document: AT Kind code of ref document: T Effective date: 20200219 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MC Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: PL |
|
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
26N | No opposition filed |
Effective date: 20201120 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: LU Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20200401 Ref country code: LI Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20200430 Ref country code: AT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: CH Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20200430 Ref country code: IT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: DE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20201103 |
|
REG | Reference to a national code |
Ref country code: BE Ref legal event code: MM Effective date: 20200430 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: SI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: BE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20200430 Ref country code: PL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20200401 |
|
REG | Reference to a national code |
Ref country code: NO Ref legal event code: MMEP |
|
GBPC | Gb: european patent ceased through non-payment of renewal fee |
Effective date: 20210401 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: NO Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210430 Ref country code: GB Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210401 Ref country code: FR Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210430 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: TR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: MT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 Ref country code: CY Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20200219 |