FR2918179A1 - Methode pour estimer la permeabilite d'un reseau de fractures a partir d'une analyse de connectivite - Google Patents
Methode pour estimer la permeabilite d'un reseau de fractures a partir d'une analyse de connectivite Download PDFInfo
- Publication number
- FR2918179A1 FR2918179A1 FR0704703A FR0704703A FR2918179A1 FR 2918179 A1 FR2918179 A1 FR 2918179A1 FR 0704703 A FR0704703 A FR 0704703A FR 0704703 A FR0704703 A FR 0704703A FR 2918179 A1 FR2918179 A1 FR 2918179A1
- Authority
- FR
- France
- Prior art keywords
- permeability
- network
- meshes
- connectivity index
- connectivity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 103
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000004458 analytical method Methods 0.000 title abstract description 10
- 239000012530 fluid Substances 0.000 claims description 10
- 229930195733 hydrocarbon Natural products 0.000 abstract description 16
- 150000002430 hydrocarbons Chemical class 0.000 abstract description 16
- 239000004215 Carbon black (E152) Substances 0.000 abstract description 12
- 238000004519 manufacturing process Methods 0.000 description 24
- 230000006399 behavior Effects 0.000 description 13
- 238000004364 calculation method Methods 0.000 description 11
- 239000006185 dispersion Substances 0.000 description 7
- 238000002347 injection Methods 0.000 description 7
- 239000007924 injection Substances 0.000 description 7
- 238000005259 measurement Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 238000005325 percolation Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 238000012886 linear function Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000006355 external stress Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- -1 injection pressure Chemical class 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH 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
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)
Abstract
- Méthode pour optimiser l'exploitation d'un gisement d'hydrocarbures fracturé, dans laquelle on détermine la perméabilité du réseau en réalisant un compromis fiable, entre méthodes numérique et analytique.- On discrétise le gisement en un ensemble de mailles et on élabore une description géométrique du réseau de fractures dans chacune des mailles. Puis, on en déduit, au sein de chaque maille, un indice de connectivité des fractures. On détermine la perméabilité du réseau de fractures des mailles dont l'indice de connectivité est supérieur à un premier seuil, et l'on affecte une valeur de perméabilité nulle au sein des autres mailles. On peut déterminer d'autres seuils de façon à choisir entre une méthode numérique ou une méthode analytique pour déterminer la perméabilité. On exploite ces perméabilités dans un simulateur d'écoulement de façon à optimiser l'exploitation du gisement.- Application à l'exploitation de gisements pétroliers.
Description
La présente invention concerne le domaine de l'optimisation de
l'exploitation de gisements souterrains, tels que des gisements d'hydrocarbures, notamment lorsque ceux-ci comportent un réseau de fractures. La méthode selon l'invention, convient notamment pour l'étude des propriétés 5 hydrauliques de terrains fracturés, et notamment pour étudier les déplacements d'hydrocarbures dans des gisements souterrains. En particulier, l'invention concerne une méthode destinée à déterminer la perméabilité d'un réseau de fractures, de façon à prédire les écoulements de fluides susceptibles de se produire au travers du gisement. On peut alors simuler une production d'hydrocarbures en 10 fonction de divers scénarios de production. L'industrie pétrolière, et plus précisément l'exploration et l'exploitation de gisements pétroliers, nécessitent d'acquérir une connaissance aussi parfaite que possible de la géologie souterraine pour fournir de façon efficace une évaluation des réserves, une modélisation de la 15 production, ou la gestion de l'exploitation. En effet, la détermination de l'emplacement d'un puits de production ou d'un puits d'injection, la constitution de la boue de forage, les caractéristiques de complétion, les paramètres nécessaires à la récupération optimale des hydrocarbures (tels que la pression d'injection, le débit de production,...) nécessitent de bien connaître le gisement. Connaître le gisement signifie connaître les propriétés pétrophysiques 20 du sous-sol en tout point de l'espace. Pour ce faire, depuis longtemps, l'industrie pétrolière allie les mesures techniques aux modélisations, réalisées en laboratoire et/ou par des logiciels. Les modélisations des gisements pétroliers constituent donc une étape technique indispensable à toute exploration ou exploitation de gisement. Ces modélisations ont pour but de fournir une description du 25 gisement. État de la technique Les ingénieurs en charge de l'exploitation de réservoirs fracturés, ont besoin de parfaitement connaître le rôle des fractures. On appelle 'fracture", une discontinuité plane, de 30 très faible épaisseur par rapport à son extension, et qui représente un plan de rupture d'une roche du gisement.
D'une part, la connaissance de la distribution et du comportement de ces fractures permet d'optimiser la localisation et l'espacement entre les puits que l'on compte forer au travers du gisement pétrolifère. D'autre part, la géométrie du réseau de fractures conditionne le déplacement des fluides tant à l'échelle du réservoir qu'à l'échelle locale où elle détermine des blocs matriciels élémentaires dans lesquels l'huile est piégée. Connaître la distribution des fractures, est donc très utile, aussi, à un stade ultérieur, pour l'ingénieur de réservoir qui cherche à calibrer les modèles qu'il construit pour simuler les gisements afin d'en reproduire ou prédire les courbes de production passées ou futures.
Les ingénieurs en charge de l'exploitation de réservoirs fracturés, ont donc besoin d'estimer la perméabilité à grande échelle (celle du rayon de drainage d'un puits ou de l'espace inter-puits par exemple) des réseaux de fractures, et de prévoir le comportement hydrodynamique (débit, pression,..) de ces réseaux, en réponse à des sollicitations extérieures imposées via des puits.
A ces fins, les spécialistes de géosciences procèdent en premier lieu à une caractérisation du réseau de fractures, sous la forme d'un ensemble de familles de fractures caractérisées par des attributs géométriques. Puis, en vue de la simulation des écoulements au sein du réservoir fracturé, un modèle numérique est le plus souvent utilisé. Ce modèle est appliqué à une représentation discrétisée du gisement, c'est-à-dire que le gisement est découpé en un ensemble de mailles. L' application du modèle numérique, nécessite de connaître les propriétés d'écoulement du réseau de fractures à l'échelle des mailles, habituellement de tailles hectométriques. En particulier, les perméabilités du réseau de fractures doivent être déterminées. Ceci peut être réalisé de manière fiable à partir d'un calcul d'écoulement, effectué sur un 25 modèle géométrique représentatif du réseau de fractures. Une telle méthode est décrite dans le brevet suivant : FR 2.757.947 (US 6.023.656). Toutefois, cette méthode de calcul numérique est coûteuse en temps de calcul, pour des réservoirs complexes et/ou de grande taille. Souvent la discrétisation d'un gisement conduit à la construction d'une grille comportant des millions de mailles. 30 Le spécialiste dispose alors de méthodes alternatives. Il dispose en effet de méthodes analytiques de calcul. On appelle méthode analytique une ou plusieurs équations permettant de déterminer de manière exacte, sans approximation ni recours à des techniques numériques (itératives, etc.) les inconnues d'un problème en fonction des données. Un exemple de méthode analytique est décrit par exemple dans le document suivant : M.Chen, M. Bai, and J-C Roegiers, Permeability Tensors of Anisotropic Fracture Networks, Mathematical Geology, Vol.31, No. 4, 1999 Toutefois, les méthodes analytiques reposent le plus souvent sur des hypothèses simplificatrices du problème physique et ne permettent pas d'obtenir la précision atteinte par les méthodes numériques qui, elles, permettent de prendre en considération la complexité réelle de la physique dans son intégralité. Il est pourtant parfois crucial de préserver une grande précision dans l'estimation des perméabilités des réseaux de fractures, de façon à pouvoir sélectionner les meilleurs scénarios de production, permettant d'optimiser la production d'hydrocarbures. L'objet de l'invention concerne une méthode pour optimiser l'exploitation d'un gisement d'hydrocarbures comportant un réseau de fractures, dans laquelle on détermine la perméabilité du réseau en réalisant un compromis fiable entre méthodes numérique et analytique. La méthode y parvient en réalisant une analyse quantitative des propriétés de connectivité du réseau de fractures, de façon à limiter l'utilisation de méthode numérique.
La méthode selon l'invention L'invention concerne une méthode pour optimiser l'exploitation d'un gisement comportant un réseau de fractures, dans laquelle on discrétise le gisement en un ensemble de mailles. On élabore également une description géométrique du réseau de fractures dans chacune des mailles. La méthode comporte les étapes suivantes : - on détermine, au sein de chaque maille, un indice de connectivité, fonction au moins du nombre d'intersections entre fractures, au moyen de ladite description géométrique ; - on estime la perméabilité du réseau de fractures de mailles dont l'indice de connectivité est supérieur à un seuil ; - on affecte une valeur fixée de perméabilité, au sein des autres mailles dont l'indice de connectivité est inférieur audit seuil, de façon à limiter le nombre d'estimations de perméabilité ; et - on optimise l'exploitation du gisement, en simulant des écoulements de fluides dans ledit gisement, en fonction des perméabilités du réseau de fractures en chaque maille. Pour optimiser davantage l'estimation de la perméabilité du réseau de fractures en chaque maille, on peut sélectionner, pour chaque maille, une méthode d'estimation de la perméabilité du réseau de fractures en fonction de la valeur de l'indice de connectivité.
La sélection de la méthode peut alors être réalisée, en définissant deux seuils de connectivité, correspondant à deux valeurs d'indice de connectivité définissant trois intervalles d'indice de connectivité. On sélectionne alors une méthode différente pour chacun des intervalles, de façon à optimiser l'estimation de la perméabilité en chaque maille. On choisira la méthode la plus simple préservant la précision des résultats.
Dans ce mode de réalisation, les seuils peuvent être définis de façon empirique, ou n réalisant les étapes suivantes : - on dispose d'un ensemble de mailles comportant chacune un réseau de fractures pour lequel on dispose d'une description géométrique ; - on détermine un indice de connectivité pour chacune des mailles ; - on détermine une perméabilité du réseau en chaque maille, à l'aide d'un simulateur d'écoulement ; - on construit une courbe de la perméabilité en fonction de l'indice de connectivité ; - on définit les seuils en fonction de la forme de ladite courbe, de façon à ce que la perméabilité obéisse à la même loi de comportement en fonction de l'indice de connectivité au sein des trois intervalles définis par les seuils. Dans ce mode de réalisation, on peut choisir l'ensemble de mailles pour lequel on dispose d'une description géométrique, en sélectionnant un ensemble de mailles issues de la discrétisation du gisement, dont les indices sont répartis sur l'intervalle des indices de connectivité calculés pour l'ensemble des mailles issues de la discrétisation du gisement.
Les méthodes d'estimation de la perméabilité peuvent être choisies de la façon suivante : - on estime la perméabilité du réseau de fractures au sein des mailles dont l'indice de connectivité est supérieur au second seuil, au moyen d'une formule analytique ; - on estime la perméabilité du réseau au sein de mailles dont l'indice de connectivité est compris entre les deux seuils, au moyen d'un simulateur d'écoulement. Dans ce cas, on peut estimer la perméabilité en fonction de la valeur de l'indice de connectivité. On peut par exemple : - estimer la perméabilité du réseau de fractures au sein des mailles dont l'indice de connectivité est supérieur au second seuil, au moyen d'une formule analytique dans laquelle on considère que la perméabilité du réseau croît linéairement en fonction de l'indice de connectivité ; - et estimer la perméabilité du réseau au sein de mailles dont l'indice de connectivité est compris entre les deux seuils, au moyen d'une méthode dans laquelle on considère que la perméabilité du réseau n'obéit plus à la même relation qu'au-delà du second seuil. 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. Présentation succincte des figures - La figure 1 représente une courbe de perméabilité de réseau K en fonction de l'indice de connectivité Ic, à partir de laquelle on détermine un seuil de percolationIf et un seuil de linéarité le (If 1 et I~ 3). 25 - La figure 2 illustre la discrétisation en deux dimensions d'un gisement en un ensemble de mailles, et indique les mailles pour lesquelles on ne calcule pas la perméabilité (zone 1, en blanc), les mailles pour lesquelles on calcule la perméabilité à l'aide d'un simulateur d'écoulement (zone 2, en gris), et les mailles pour lesquelles on calcule la perméabilité à l'aide d'une formule linéaire (zone 3, en noir). 30 Description détaillée de la méthode La méthode selon l'invention permet d'optimiser l'exploitation d'un gisement d'hydrocarbures, notamment lorsque celui-ci comporte un réseau de fractures. La méthode permet en particulier de minimiser le temps nécessaire à la détermination des perméabilités du réseau de fractures, tout en préservant une bonne précision des résultats. La méthode comporte six étapes : 1- Discrétisation du gisement en un ensemble de mailles 2-Description géométrique du réseau de fractures 3- Analyse de la connectivité du réseau de fractures 4- Détermination de la perméabilité équivalente d'un réseau de fractures 5- Simulation des écoulements de fluides 6- Optimisation des conditions de production du gisement 1-Discrétisation du gisement en un ensemble de mailles Depuis longtemps, l'industrie pétrolière allie les mesures techniques aux modélisations, réalisées en laboratoire et/ou par des logiciels. Les modélisations des gisements pétroliers, constituent donc une étape technique indispensable à toute exploration ou exploitation de gisement. Ces modélisations ont pour but de fournir une description du gisement, via son architecture sédimentaire et/ou ses propriétés pétrophysiques. Ces modélisations se basent sur une représentation du gisement, en un ensemble de mailles. Chacune de ces mailles représentent un volume donné du gisement. L'ensemble des mailles constitue une représentation discrète du gisement. 2- Description géométrique du réseau de fractures naturelles Le spécialiste des géosciences procède à une caractérisation de la géométrie du réseau de fractures naturelles : il élabore une description géométrique du réseau de fractures, dans chacune des mailles, au moyen d'attributs géométriques pertinents. 6 Cette description géométrique nécessite un ensemble de mesures, réalisées sur le terrain par le géologue. Ces mesures permettent de caractériser le réseau de fractures, de façon à aboutir à une description du réseau sous la forme d'un ensemble de N familles de fractures, caractérisées par des attributs géométriques.
Dans tout ce qui suit, pour des raisons de clarté de l'exposé, nous considérons la situation bidimensionnelle d'un réseau de fractures de N familles contenues dans une couche, sans toutefois que cela ne remette en cause la possibilité d'étendre le domaine d'application de l'invention aux situations tridimensionnelles de réservoirs multicouches et/ou de grande épaisseur par rapport à l'extension verticale des fractures.
En deux dimensions les attributs géométriques relatifs à une famille f peuvent être les suivants : - un angle moyen, Of , d'orientation par rapport à une direction de référence et un angle moyen de dispersion, a f , de cette orientation autour de l'angle moyen ef . Ces paramètres d'orientation sont généralement ajustés à une loi statistique (telle que par exemple la loi de Fisher) ; - une longueur moyenne, L f , et une dispersion autour de cette moyenne ; - une densité, d f , définie comme la longueur cumulée de fracture par unité de surface (m/m2)• Cette description géométrique du réseau de fractures, peut également être déterminée de façon probabiliste. On établit alors une description géométrique du réseau de fractures en affectant à chaque famille de fractures f une loi de probabilité pe f des orientations dans le plan des couches par rapport à une direction de référence, ainsi qu'une loi de probabilité des longueurs p, f , et une densité d f . Chaque loi de probabilité pour le paramètre p vérifie la relation : f pn,f dp =1 En fonction des mesures réalisées, ou des lois de probabilité définies, on affecte à chaque maille de la représentation discrète du gisement, une valeur à chacun des attributs géométriques décrivant le réseau de fractures à l'échelle de cette maille. Les documents suivants décrivent un exemple de techniques utilisables pour accomplir 5 cette tâche : FR 2.725.794 (US 5.661.698), FR 2.725.814 (US 5.798.768), FR 2.733.073 (US 5.659.135) A la fin de cette étape, on a construit une représentation du gisement sous forme d'un ensemble de mailles, et l'on a affecté à chacune de ces mailles un ensemble d'attributs géométriques caractérisant le réseau de fractures au sein de chaque maille. Pour une application au cas 3D, les attributs géométriques permettant d'établir une description géométrique du réseau de fractures naturelles, sont les mêmes paramètres que pour le cas 2D ainsi que : - un angle moyen, yrf, d'orientation par rapport à une direction de référence dans le 15 plan vertical et un angle moyen de dispersion, 33f , de cette orientation autour de l'angle moyen yr f . Ces paramètres d'orientation sont généralement ajustés à une loi statistique ; - une hauteur moyenne, Hf , et une dispersion autour de cette moyenne ; - une densité, d f , définie comme la surface cumulée de fracture par unité de volume (m2/m3). 20 On établit une description géométrique du réseau de fractures en affectant à chaque famille de fractures f une loi de probabilité pe f des orientations dans le plan des couches par rapport à une direction de référence, une loi de probabilité des orientations dans le plan vertical p,, f , une loi de probabilité des longueurs p, f ainsi qu'une loi de probabilité des hauteurs pH ,f , et une densité d f . 25 3- Analyse de la connectivité du réseau de fractures Dans l'optique d'optimiser l'exploitation d'un gisement, on prend en compte, non seulement la géométrie du réseau de fractures, mais également le rôle des fractures sur le comportement hydrodynamique du réservoir. Pour déterminer ce rôle, on définit si le réseau 10 de fractures est connecté, de telle sorte qu'il contribue directement aux écoulements et transports à l'échelle du réservoir. La connaissance de ce degré de connectivité est essentielle pour l'ingénieur de réservoir chargé d'estimer/prédire l'exploitation du réservoir. Selon l'invention, avant de calculer la perméabilité du réseau en chaque maille, qui est soit coûteux, lorsque l'on utilise une méthode précise telle qu'une méthode numérique, soit rapide mais imprécis, lorsque l'on utilise une méthode telle qu'une méthode analytique, on évalue la connectivité du réseau de fractures de la maille considérée. En effet, si les fractures d'un réseau, au sein d'une maille, ne sont pas connectées entre elles, la perméabilité du réseau est nulle. A l'inverse, si les fractures d'un réseau, au sein d'une maille, sont toutes connectées, la perméabilité est importante. En effet, un fluide n'a aucune difficulté à traverser la maille dans ce dernier cas. Pour déterminer le degré de connectivité des fractures d'un réseau au sein d'une maille, on calcule, selon l'invention, un indice représentatif du nombre d'intersections entre les fractures du réseau. En effet, plus les fractures d'un réseau comportent d'intersections, plus elles sont connectées. Cet indice est appelé indice de connectivité , et il est noté Ic. L'indice de connectivité le est alors un paramètre fonction du nombre d'intersections entre les fractures du réseau. On le détermine en chaque maille, à partir des informations issues de la description géométrique.
De façon générale, on peut utiliser la formulation suivante pour définir l'indice de connectivité Ic. le = g, (d )g2 (0, a)g3 (L) avec : g, : une fonction linéaire dépendant de la densité de fracturation du réseau ; g2 : une fonction dépendant de la dispersion des orientations (a) dans le plan horizontal g3 : une fonction linéaire dépendant de la longueur moyenne (L) des fractures ; g4 : une fonction dépendant de la dispersion des orientations (yr) dans le plan vertical. On peut évaluer ces fonctions de façon empirique. On peut aussi définir l'indice de connectivité à partir des valeurs moyennes des attributs géométriques définis plus haut, de la façon suivante : Cas de 2 familles (i, j) d'orientations Bi , 9j constantes: d,.d j .L..L.sin(19i -9jl) Icij = d,Lj+djLi Cas d'une famille [dont la dispersion d'orientations est non négligeable et définie par une loi géostatistique Pe f : 211 2H Icf =d f xLf x f f Pe,f(el)Pe,f(02)Isin(9, -02)Id9,d92 e,=0 92>e, Si on considère les lois de distribution statistiques des paramètres géométriques de chacune des familles f, on peut utiliser l'expression suivante de I, pour un cas à N Familles : N 211 211 E die f f po,; (e1)P0,i (e2)Isin(91 û 92 )I d 9, d 92 1+ i=1 e,=0 e2>6, E Edidi Zf 2fpe,i(ei)Pe,j(02)Isin(91 -02)Id9,d92 i=1 j=i+l 9,=0 92=0 4- Détermination de la perméabilité équivalente du réseau de fractures Pour déterminer le comportement hydrodynamique d'un réservoir, il est nécessaire d'évaluer une perméabilité du réseau de fractures à grande échelle. On calcule alors, en 20 chaque maille, une perméabilité, dite équivalente , du réseau de fractures contenu dans cette maille. Deux méthodes existent : l'une numérique, coûteuse en ressources de calcul pour les réservoirs de grandes dimensions (nombreuses mailles), l'autre, analytique, rapide mais15 approximative car fondée sur des hypothèses simplificatrices, relatives à la géométrie du réseau par exemple. Grâce à l'invention, l'ingénieur de réservoir peut optimiser en coût (temps) et qualité (précision) le calcul des perméabilités de fracture.
Le calcul des perméabilités selon l'invention, s'effectue en analysant tout d'abord la valeur de l'indice de connectivité Ic. Sélection de mailles pour lesquelles il est nécessaire de déterminer la perméabilité Le principe est le suivant : - Si l'indice de connectivité le de la maille considérée indique que les fractures sont déconnectées entre elles, on considère la perméabilité du réseau nulle à grande échelle. Hormis le cas de fractures/failles de grande extension, le rôle des fractures sur le comportement hydrodynamique à grande échelle du réservoir est négligeable (perméabilité du réseau nulle). Dans ce cas, il n'est donc pas nécessaire de calculer la perméabilité du réseau.
Ceci peut concerner des millions de mailles dans un modèle de réservoir fracturé, et l'on évite ainsi un très grand nombre de calculs inutiles. -Si l'indice de connectivité le indique que les fractures sont connectées entre elles, on considère que le réseau acquiert une perméabilité à grande échelle. Le rôle hydraulique des fractures risque d'être sensible et ces dernières doivent être intégrées dans l'étude de la dynamique du réservoir. La valeur seuil de l'indice de connectivité, à partir duquel on considère qu'il est nécessaire de calculer la perméabilité, peut être obtenue de façon empirique, ou encore par des simulations. L'homme du métier pourra notamment utiliser un simulateur d'écoulement, logiciel bien connu des spécialistes, pour définir ce seuil. Ce seuil est appelé seuil de percolation. Il est noté IÉ . Ainsi, l'évaluation de la connectivité du réseau de fractures en chaque maille, permet de sélectionner les mailles de la discrétisation du gisement, pour lesquelles il est nécessaire de déterminer la perméabilité de réseau par une méthode de calcul appropriée. Les autres mailles ont une valeur nulle de perméabilité de réseau. Sélection d'une méthode de détermination de la perméabilité du réseau de fractures Selon un mode de réalisation, on peut exploiter davantage l'indice de connectivité ainsi calculé. En effet, en établissant une courbe de perméabilité en fonction de l'indice de connectivité, on peut définir des comportements de la perméabilité, permettant de définir la technique de détermination la mieux adaptée. Selon un mode de réalisation général, on sélectionne la méthode de calcul de la perméabilité, en définissant des seuils de connectivité, correspondant à des valeurs d'indice de connectivité définissant des intervalles d'indice de connectivité. On sélectionne une méthode pour chacun desdits intervalles. Ces seuils peuvent être définis de façon empirique, ou par exemple, en réalisant les étapes suivantes: i- on dispose d'un ensemble de mailles comportant chacune un réseau de fractures pour lequel on dispose d'une description géométrique ; ii- on calcule l'indice de connectivité de chacune de ces mailles ; iii- on détermine une perméabilité du réseau en chaque maille, à l'aide d'un simulateur d'écoulement par exemple ; iv- on construit une courbe de la perméabilité en fonction de l'indice de connectivité ; v- on définit les seuils en fonction de la forme de cette courbe, de façon à ce que la perméabilité varie en fonction de l'indice de connectivité selon un comportement homogène au sein des trois intervalles définis par les seuils.
On entend par comportement homogène, le fait que, sur un intervalle, la courbe de perméabilité obéit à la même loi de comportement en fonction de l'indice de connectivité. La courbe de perméabilité en fonction de l'indice de connectivité peut être alors modélisée par une formule analytique unique (loi linéaire, polynomiale, etc.). En d'autres termes, sur un intervalle, le réseau possède la même loi de comportement à l'écoulement, i.e. la même loi de perméabilité (comportement hydraulique) en fonction de l'indice de connectivité. 5 En pratique, l'ensemble de mailles de l'étape i, peut être défini de la façon suivante : ayant calculé l'indice de connectivité pour toutes les mailles de la discrétisation du gisement, on sélectionne un ensemble de mailles, dont les indices sont répartis sur l'intervalle des indices de connectivité calculés pour l'ensemble des mailles du gisement. Selon un mode de réalisation particulier, on définit deux seuils de connectivité, définissant trois intervalles d'indice de connectivité. La figure 1 illustre une telle approche. Cette figure représente une courbe de perméabilité de réseau, K, en fonction de l'indice de connectivité Ic. On constate qu'il existe 10 un premier seuil. Ce seuil correspond au seuil de percolation If. . Il est défini sur la figure 1 par I~ z- 1. Il existe un second seuil, noté , appelé seuil de linéarité. Il est défini sur la figure 1 par 4 z- 3 . Au-delà de ce seuil, la courbe est une droite. Ces deux seuils définissent trois intervalles sur lesquels la perméabilité varie selon un comportement homogène en fonction de l'indice de connectivité : en dessous de If., la 15 perméabilité est constante (nulle), au dessus de I~ , la perméabilité croît linéairement. Entre les deux seuils, la perméabilité évolue en fonction de l'indice de connectivité selon une relation unique et non linéaire. L'indice de connectivité, calculé pour chaque maille du modèle de champ fracturé, est alors utilisé de la façon suivante : 20 - le Ic K(Ic)=0 - le : sur cet intervalle, la perméabilité du réseau de fractures croît linéairement en fonction de l'indice de connectivité (Ic) ou de la densité de fractures (d). On peut alors utiliser une méthode de calcul, permettant de déterminer la perméabilité à partir d'une formule analytique. On peut par exemple définir la formule suivante : 25 K(Ic)=a.Ic+b Les coefficients a et h peuvent être déterminés par une simple régression linéaire. - If. <_ le : sur cet intervalle de transition, la perméabilité K du réseau évolue suivant une certaine fonction g dépendant de l'indice de connectivité (Ic) ou de la densité de fractures (d) : K(Ic)=g(Ic) La fonction g est une fonction distincte de la fonction linéaire définie sur l'intervalle le >_ le . Elle est fixée pour un type de réseau donné, c'est-à-dire pour les réseaux dont seule la densité varie (à nombre N de familles, orientations et longueurs de fractures de chaque famille fixés). Dans ce cas, une méthode numérique de calcul des perméabilités de fracture permet d'obtenir une valeur précise de perméabilité. Une telle méthode est décrite dans le document suivant : FR 2.757.947 (US 6.023.656). Cependant, une alternative à la méthode numérique peut être adoptée de façon à augmenter la rapidité des calculs de perméabilité. Elle consiste à faire usage d'une approximation, telle qu'une formule analytique donnant l'évolution de la perméabilité en fonction de l'indice de connectivité. En conclusion, l'évaluation de la connectivité du réseau de fractures en chaque maille, permet de sélectionner une méthode de détermination de la perméabilité adaptée au besoin requis au niveau de chaque maille (i.e. une méthode fiable d'une part, rapide et peu coûteuse en temps de calcul d'autre part). Elle permet par la même occasion de définir trois régions du champ (ou ensemble de mailles) possédant chacune un comportement homogène en perméabilité de fracture. 5- Simulation des écoulements de fluides A ce stade, l'ingénieur réservoir dispose d'une représentation discrétisée (ensemble de mailles) du gisement d'hydrocarbures, dont il souhaite extraire les hydrocarbures. Cette représentation est renseignéeen perméabilité de réseau de fractures, c'est-à-dire qu'à chaque maille, est associée une valeur de perméabilité. L'ingénieur réservoir choisit un procédé de production, par exemple le procédé de récupération par injection d'eau, dont il demeure ensuite à préciser le scénario optimal de mise en oeuvre pour le champ considéré. La définition d'un scénario optimal d'injection d'eau consistera par exemple à fixer le nombre et l'implantation (position et espacement) des puits injecteurs et producteurs afin de tenir compte au mieux de l'impact des fractures sur la progression des fluides au sein du réservoir.
En fonction du scénario choisi, et des perméabilités de réseau de fractures, on est alors capable de simuler la production d'hydrocarbures escomptée, au moyen d'un outil bien connu des spécialistes : un simulateur d'écoulement. Un tel logiciel permet de simuler les écoulements de fluides au sein de gisement. 6- Optimisation des conditions de production du gisement En sélectionnant divers scénarios caractérisés par exemple par diverses implantations respectives des puits injecteurs et producteurs, et en simulant la production d'hydrocarbures pour chacun d'eux selon l'étape 5, on peut sélectionner le scénario permettant d'obtenir une 10 production optimale du gisement. On optimise alors l'exploitation du gisement, en mettant en oeuvre, sur le champ, le scénario de production ainsi sélectionné. Exemple d'application. On discrétise un gisement d'hydrocarbures comportant un réseau de fractures. La figure 15 2 illustre le résultat de ce maillage en deux dimensions. On élabore une description géométrique du réseau de fractures dans chacune des mailles, à l'aide d'informations issues de mesures et analyses géologiques. On détermine, au sein de chaque maille, un indice de connectivité, défini par exemple par la formulation suivante (formule fondée sur les attributs moyens de chacune des familles 20 de fractures) : N N d; .d i .I sin(B; ù 0j )I = i =1 i=i+1 le Cet indice définit le nombre d'intersection moyen entre fractures, au sein de chaque maille. On tient alors compte de l'existence des deux seuils de percolation, I" z- 1, et de 25 linéarité, le . 3, pour classer les mailles en fonction de leur indice de connectivité et définir ainsi trois zones (ou régions) du champ caractérisées respectivement par des valeurs de cet5 indice inférieures à 1É, supérieures à et comprises entre I' et l'c.. Cette définition de zones détermine le choix de la méthode de calcul des perméabilités de fracture. La figure 2 illustre, en deux dimensions, les mailles de la représentation du gisement pour lesquelles on ne calcule par la perméabilité (zone 1 où le , en blanc), les mailles pour lesquelles on calcule la perméabilité à l'aide d'un simulateur d'écoulement (zone 2 où le 5 le <_ le , en gris), et les mailles pour lesquelles on calcule la perméabilité à l'aide d'une formule linéaire (zone 2 où le , en gris). Sur la zone 1 on ne calcule pas la perméabilité. On gagne donc un temps de calcul précieux. Sur la zone 3, on effectue un calcul linéaire qui nous donne la même précision qu'une simulation numérique. Sur la zone 2, pour obtenir une précision importante, on utilise un simulateur d'écoulement. On sélectionne ensuite un procédé de production, par exemple l'injection d'eau. Le mode de mise en oeuvre de ce procédé pour le champ considéré demeure toutefois à spécifier, et plus particulièrement encore si ce champ s'avère fracturé. Différents scénarios de mise en oeuvre, différant les uns des autres par la position des puits par exemple, sont alors définis et comparés sur la base de critères quantitatifs de production/récupération des fluides en place. L'évaluation (prévision) de ces critères de production requiert l'usage d'un simulateur de champ à même de reproduire (simuler) chacun des scénarios. Dans le cas de réservoirs fracturés, les perméabilités du réseau de fractures à l'échelle de résolution du simulateur (la maille de réservoir) constituent des informations de base indispensables pour effectuer ces simulations, et déterminantes pour garantir la fiabilité des prévisions de production. Avantages L'invention permet d'estimer la perméabilité à grande échelle (échelle du rayon de 25 drainage d'un puits ou de l'espace inter-puits par exemple) de ces fractures, de façon rapide et précise. Il est alors possible de prévoir le comportement hydrodynamique (débit, pression,..) en réponse à des sollicitations extérieures imposées via des puits lors de la production d'hydrocarbures.
Les ingénieurs en charge de l'exploitation du gisement ont alors un outil leur permettant de rapidement évaluer la performance de différents scénarios de production, et ainsi, de sélectionner celui qui optimise l'exploitation au regard des critères sélectionnés par l'opérateur, comme d'assurer une production d'hydrocarbure optimale.
Ainsi, l'invention trouve une application industrielle dans l'exploitation de gisements souterrains, comportant un réseau de fractures. Il peut s'agir d'un gisement d'hydrocarbures pour lequel on souhaite optimiser la production, ou un gisement de stockage de gaz par exemple, pour lequel on souhaite optimiser l'injection ou les conditions de stockage.
Claims (9)
1) Méthode pour optimiser l'exploitation d'un gisement comportant un réseau de fractures, dans laquelle on discrétise le gisement en un ensemble de mailles et on élabore une description géométrique du réseau de fractures dans chacune des mailles, caractérisée en ce que la méthode comporte les étapes suivantes : - on détermine, au sein de chaque maille, un indice de connectivité, fonction au moins du nombre d'intersections entre fractures, au moyen de ladite description géométrique ; on estime la perméabilité du réseau de fractures de mailles dont l'indice de connectivité est supérieur à un seuil ; - on affecte une valeur fixée de perméabilité, au sein des autres mailles dont l'indice de connectivité est inférieur audit seuil, de façon à limiter le nombre d'estimations de perméabilité ; et on optimise l'exploitation du gisement, en simulant des écoulements de fluides dans ledit gisement, en fonction des perméabilités du réseau de fractures en chaque maille.
2) Méthode selon la revendication 1, dans laquelle on estime la perméabilité du réseau de fractures, en sélectionnant, pour chaque maille, une méthode d'estimation de la perméabilité du réseau de fractures en fonction de la valeur de l'indice de connectivité.
3) Méthode selon la revendication 2, dans laquelle on sélectionne la méthode en définissant deux seuils de connectivité, correspondant à deux valeurs d'indice de connectivité définissant trois intervalles d'indice de connectivité, et l'on sélectionne une méthode différente pour chacun desdits intervalles, de façon à optimiser l'estimation de la perméabilité en chaque maille.
4) Méthode selon la revendication 3, dans laquelle on définit les seuils de façon empirique.
5) Méthode selon la revendication 3, dans laquelle on définit les seuils en réalisant les étapes suivantes : - on dispose d'un ensemble de mailles comportant chacune un réseau de fractures pour lequel on dispose d'une description géométrique ; - on détermine un indice de connectivité pour chacune des mailles ;- on détermine une perméabilité du réseau en chaque maille, à l'aide d'un simulateur d'écoulement ; on construit une courbe de la perméabilité en fonction de l'indice de connectivité ; on définit les seuils en fonction de la forme de ladite courbe, de façon à ce que la perméabilité obéisse à la même loi de comportement en fonction de l'indice de connectivité au sein des trois intervalles définis par les seuils.
6) Méthode selon la revendication 5, dans laquelle l'ensemble de mailles comportant chacune un réseau de fractures pour lequel on dispose d'une description géométrique, est déterminé en sélectionnant un ensemble de mailles issues de la discrétisation du gisement, dont les indices sont répartis sur l'intervalle des indices de connectivité calculés pour l'ensemble des mailles issues de la discrétisation du gisement.
7) Méthode selon l'une des revendications 3 à 6, dans laquelle : - on estime la perméabilité du réseau de fractures au sein des mailles dont l'indice de connectivité est supérieur au second seuil, au moyen d'une formule analytique ; on estime la perméabilité du réseau au sein de mailles dont l'indice de connectivité est compris entre les deux seuils, au moyen d'un simulateur d'écoulement.
8) Méthode selon la revendication 7, dans laquelle on estime la perméabilité en fonction de la valeur de l'indice de connectivité.
9) Méthode selon la revendication 8, dans laquelle : - on estime la perméabilité du réseau de fractures au sein des mailles dont l'indice de connectivité est supérieur au second seuil, au moyen d'une formule analytique dans laquelle on considère que la perméabilité du réseau croît linéairement en fonction de l'indice de connectivité ; - on estime la perméabilité du réseau au sein de mailles dont l'indice de connectivité est compris entre les deux seuils, au moyen d'une méthode dans laquelle on considère que la perméabilité du réseau n'obéit plus à la même relation qu'au-delà du second seuil.
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0704703A FR2918179B1 (fr) | 2007-06-29 | 2007-06-29 | Methode pour estimer la permeabilite d'un reseau de fractures a partir d'une analyse de connectivite |
| EP08290574.6A EP2037080B1 (fr) | 2007-06-29 | 2008-06-17 | Méthode pour estimer la preméabilité d'un réseau de fractures à partir d'une analyse de connectivité |
| US12/146,832 US8078405B2 (en) | 2007-06-29 | 2008-06-26 | Method of estimating the permeability of a fracture network from a connectivity analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0704703A FR2918179B1 (fr) | 2007-06-29 | 2007-06-29 | Methode pour estimer la permeabilite d'un reseau de fractures a partir d'une analyse de connectivite |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2918179A1 true FR2918179A1 (fr) | 2009-01-02 |
| FR2918179B1 FR2918179B1 (fr) | 2009-10-09 |
Family
ID=39148267
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR0704703A Active FR2918179B1 (fr) | 2007-06-29 | 2007-06-29 | Methode pour estimer la permeabilite d'un reseau de fractures a partir d'une analyse de connectivite |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US8078405B2 (fr) |
| EP (1) | EP2037080B1 (fr) |
| FR (1) | FR2918179B1 (fr) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3144468A1 (fr) | 2015-09-15 | 2017-03-22 | IFP Energies nouvelles | Procédé pour caracteriser le réseau de fractures d'un gisement fracture et procédé pour l'exploiter |
| EP3181804A1 (fr) | 2015-12-17 | 2017-06-21 | IFP Énergies nouvelles | Procede pour caracteriser et exploiter une formation souterraine comprenant un reseau de fractures |
| CN115184172A (zh) * | 2022-07-11 | 2022-10-14 | 西安石油大学 | 一种非常规油气储集层岩石力学特征分析装置 |
Families Citing this family (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8437997B2 (en) * | 2008-01-22 | 2013-05-07 | Exxonmobil Upstream Research Company | Dynamic connectivity analysis |
| AU2010271128B2 (en) * | 2009-07-06 | 2015-10-29 | Exxonmobil Upstream Research Company | Method for seismic interpretation using seismic texture attributes |
| WO2011049648A1 (fr) | 2009-10-20 | 2011-04-28 | Exxonmobil Upstream Research Company | Procédé d'évaluation quantitative de la connectivité pour des paires de puits à des fréquences variables |
| CN101942991A (zh) * | 2010-06-30 | 2011-01-12 | 中国石油大学(北京) | 裂缝各向异性油藏注水开发可预测物理模型建立方法 |
| FR2967200A1 (fr) | 2010-11-10 | 2012-05-11 | IFP Energies Nouvelles | Methode pour caracteriser le reseau de fractures d'un gisement fracture et methode pour l'exploiter |
| FR2976099B1 (fr) | 2011-06-01 | 2013-05-17 | IFP Energies Nouvelles | Methode pour construire un maillage d'un reseau de fractures a partir de diagramme de voronoi |
| FR2981475B1 (fr) | 2011-10-12 | 2013-11-01 | IFP Energies Nouvelles | Methode pour construire un maillage d'un reservoir fracture avec un nombre limite de noeuds dans le milieu matrice |
| RU2015123454A (ru) | 2012-12-13 | 2017-01-16 | Лэндмарк Графикс Корпорейшн | Система, способ и компьютерный программный продукт для оценки и классификации геологических тел с использованием эйлеровой характеристики |
| CA2897656C (fr) | 2013-02-08 | 2021-07-27 | Schlumberger Canada Limited | Appareil et methodologie pour mesurer des proprietes de materiau microporeux a des echelles multiples |
| CN103257089B (zh) * | 2013-04-08 | 2015-06-03 | 中国石油天然气股份有限公司 | 压力脉冲测量装置测量基质和裂缝渗透率的方法 |
| CN103323887B (zh) * | 2013-07-09 | 2016-03-02 | 中煤科工集团重庆研究院有限公司 | 一种煤矿采动稳定区煤层气储量的评估方法及系统 |
| CN104634713B (zh) * | 2013-11-13 | 2017-05-17 | 中国石油化工股份有限公司 | 特高含水期油水相渗比值曲线非线性关系表征方法 |
| CN105089650B (zh) * | 2014-04-22 | 2018-03-13 | 中国石油化工股份有限公司 | 评价储盖组合的方法及储层产能预测方法 |
| CN104847339B (zh) * | 2014-12-09 | 2017-09-22 | 中国石油集团川庆钻探工程有限公司 | 利用压化指数评价页岩气储层的方法 |
| US10242136B2 (en) * | 2015-05-20 | 2019-03-26 | Saudi Arabian Oil Company | Parallel solution for fully-coupled fully-implicit wellbore modeling in reservoir simulation |
| US11294095B2 (en) | 2015-08-18 | 2022-04-05 | Schlumberger Technology Corporation | Reservoir simulations with fracture networks |
| EP3338115A4 (fr) * | 2015-08-18 | 2019-04-24 | Services Petroliers Schlumberger | Simulations de réservoir incluant des réseaux de fracture |
| US10302814B2 (en) * | 2015-08-20 | 2019-05-28 | Baker Hughes, A Ge Company, Llc | Mechanisms-based fracture model for geomaterials |
| US10546072B2 (en) | 2016-03-28 | 2020-01-28 | Baker Huges, A Ge Company, Llc | Obtaining micro- and macro-rock properties with a calibrated rock deformation simulation |
| CN106845786B (zh) * | 2016-12-26 | 2021-01-29 | 中国石油天然气股份有限公司 | 注入井和生产井间储层连通性评价方法和装置 |
| EP3568568B1 (fr) * | 2017-01-13 | 2023-10-11 | Baker Hughes Holdings LLC | Réseaux de fractures discrètes à proximité d'un puits de forage |
| CN109374498B (zh) * | 2018-10-29 | 2021-06-25 | 河海大学 | 一种单裂隙岩体渗流应力耦合系统及方法 |
| CN111535849B (zh) * | 2020-05-15 | 2021-06-01 | 中国石油大学(北京) | 一种瓦斯排采时间计算方法、装置及设备 |
| CN112903565B (zh) * | 2021-02-01 | 2022-10-18 | 核工业北京地质研究院 | 考虑岩石裂隙内部几何特征的渗透率测定方法 |
| US11933165B2 (en) | 2021-03-15 | 2024-03-19 | Saudi Arabian Oil Company | Hydraulic fracture conductivity modeling |
| US11525935B1 (en) | 2021-08-31 | 2022-12-13 | Saudi Arabian Oil Company | Determining hydrogen sulfide (H2S) concentration and distribution in carbonate reservoirs using geomechanical properties |
| US11921250B2 (en) | 2022-03-09 | 2024-03-05 | Saudi Arabian Oil Company | Geo-mechanical based determination of sweet spot intervals for hydraulic fracturing stimulation |
| CN118940552B (zh) * | 2024-10-14 | 2024-12-27 | 湖南科技大学 | 路域水体功能连通性阈值确定方法及趋势研判分析系统 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2757947A1 (fr) * | 1996-12-30 | 1998-07-03 | Inst Francais Du Petrole | Methode pour determiner la permeabilite equivalente d'un reseau de fracture dans un milieu souterrain multi-couches |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2725814B1 (fr) | 1994-10-18 | 1997-01-24 | Inst Francais Du Petrole | Methode pour cartographier par interpolation, un reseau de lignes, notamment la configuration de failles geologiques |
| FR2725794B1 (fr) | 1994-10-18 | 1997-01-24 | Inst Francais Du Petrole | Methode pour modeliser la distribution spatiale d'objets geometriques dans un milieu, tels que des failles dans une formation geologique |
| FR2733073B1 (fr) | 1995-04-12 | 1997-06-06 | Inst Francais Du Petrole | Methode pour modeliser un milieu geologique stratifie et fracture |
| US5757663A (en) * | 1995-09-26 | 1998-05-26 | Atlantic Richfield Company | Hydrocarbon reservoir connectivity tool using cells and pay indicators |
| US7584086B2 (en) * | 2003-09-30 | 2009-09-01 | Exxonmobil Upstream Research Company | Characterizing connectivity in reservoir models using paths of least resistance |
| CA2598955C (fr) * | 2004-11-29 | 2014-10-14 | Chevron U.S.A. Inc. | Procede, systeme et dispositif de stockage de programme afin de simuler l'ecoulement de fluides dans un systeme physique au moyen d'une architecture a orientation d'objet extensible a base de composition dynamique |
-
2007
- 2007-06-29 FR FR0704703A patent/FR2918179B1/fr active Active
-
2008
- 2008-06-17 EP EP08290574.6A patent/EP2037080B1/fr active Active
- 2008-06-26 US US12/146,832 patent/US8078405B2/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2757947A1 (fr) * | 1996-12-30 | 1998-07-03 | Inst Francais Du Petrole | Methode pour determiner la permeabilite equivalente d'un reseau de fracture dans un milieu souterrain multi-couches |
Non-Patent Citations (5)
| Title |
|---|
| LAETITIA MACÉ, LAURENT SOUCHE AND JEAN-LAURENT MALLET: "3D fracture modeling integrating geomechanics and geologic data.", AAPG INTERNATIONAL CONFERENCE, 24 October 2004 (2004-10-24) - 27 October 2004 (2004-10-27), Cancun, pages 1 - 6, XP002472551 * |
| LIYONG LI AND SEONG H. LEE: "Efficient field-scale simulation of black oil in a naturally fractured reservoir via discrete fracture networks and homogenized media.", SOCIETY OF PETROLEUM ENGINEERS, vol. SPE, no. 103901, 5 December 2006 (2006-12-05) - 7 December 2006 (2006-12-07), Beijing, pages 1 - 12, XP002472552 * |
| M. MASIHI, P.R. KING AND P.R. NURAFZA: "Fast estimation of performance parameters in fractured reservoirs using percolation theory", SOCIETY OF PETROLEUM ENGINEERS, vol. SPE, no. 94186, 13 June 2005 (2005-06-13) - 16 June 2005 (2005-06-16), Madrid, pages 1 - 12, XP002472553 * |
| R. BASQUET, C.E. COHEN AND B. BOURBIAUX: "Fracture flow property identification: an optimized implementation of discrete fracture network models", SOCIETY OF PETROLEUM ENGINEERS, vol. SPE, no. 93748, 12 March 2005 (2005-03-12) - 15 March 2005 (2005-03-15), Bahrain, pages 1 - 9, XP002472554 * |
| SAIT I. OZKAYA, JOERG MATTNER: "Fracture connectivity from fracture intersections in borehole image logs", COMPUTERS & GEOSCIENCES, no. 29, 2003, pages 143 - 153, XP002472550 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3144468A1 (fr) | 2015-09-15 | 2017-03-22 | IFP Energies nouvelles | Procédé pour caracteriser le réseau de fractures d'un gisement fracture et procédé pour l'exploiter |
| US10288544B2 (en) | 2015-09-15 | 2019-05-14 | IFP Energies Nouvelles | Method for characterizing the fracture network of a fractured reservoir and method for exploiting it |
| EP3181804A1 (fr) | 2015-12-17 | 2017-06-21 | IFP Énergies nouvelles | Procede pour caracteriser et exploiter une formation souterraine comprenant un reseau de fractures |
| CN115184172A (zh) * | 2022-07-11 | 2022-10-14 | 西安石油大学 | 一种非常规油气储集层岩石力学特征分析装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| EP2037080A1 (fr) | 2009-03-18 |
| EP2037080B1 (fr) | 2019-10-23 |
| US20090005996A1 (en) | 2009-01-01 |
| FR2918179B1 (fr) | 2009-10-09 |
| US8078405B2 (en) | 2011-12-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2037080B1 (fr) | Méthode pour estimer la preméabilité d'un réseau de fractures à partir d'une analyse de connectivité | |
| EP3144468B1 (fr) | Procédé pour caracteriser le réseau de fractures d'un gisement fracture et procédé pour l'exploiter | |
| EP2453106B1 (fr) | Methode pour caractériser le réseau de fractures d'un gisement fracture et méthode pour l'exploiter | |
| EP2466338B1 (fr) | Procédé d'exploitation d'un gisement pétrolier à partir d'un modèle d'écoulement calé au moyen d'une mise à l'échelle de cartes de pression et de saturation | |
| EP2212722B1 (fr) | Methode de modelisation d'un milieu géologique poreux traversé par un réseau de fractures | |
| EP2530493B1 (fr) | Méthode pour construire un maillage d'un réseau de fractures à partir de diagramme de voronoi | |
| 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 | |
| FR2823877A1 (fr) | Methode pour contraindre par des donnees dynamiques de production un modele fin representatif de la repartition dans le gisement d'une grandeur physique caracteristique de la structure du sous-sol | |
| FR3045868A1 (fr) | Procede pour caracteriser et exploiter une formation souterraine comprenant un reseau de fractures | |
| FR3043481A1 (fr) | ||
| WO2009034253A1 (fr) | Procede, programme et systeme informatique de mise a l'echelle de donnees de modele de reservoir d'hydrocarbure | |
| EP3719251A1 (fr) | Procede pour exploiter un reservoir petrolier fracture ayant une taille de pores heterogene | |
| FR2992448A1 (fr) | Procede d'exploitation d'un reservoir geologique a partir d'un modele de reservoir cale au moyen d'un parametrage multi-echelles | |
| CA2821099C (fr) | Procede d'exploitation d'un reservoir geologique a partir d'un modele de reservoir cale par le calcul d'une loi analytique de distribution conditionnelle de parametres incertains du modele | |
| CA2854085C (fr) | Methode pour optimiser l'exploitation d'un gisement de fluide par prise en compte d'un terme d'echange geologique et transitoire entre blocs matriciels et fractures | |
| EP2770162B1 (fr) | Procédé d'exploitation d'un réservoir géologique au moyen d'un modèle de réservoir calé et cohérent vis à vis des propriétés d'écoulement | |
| FR2956746A1 (fr) | Methode pour exploiter un milieu souterrain a partir de calage d'historique d'un modele geologique comportant un reseau de failles. | |
| EP2594733A1 (fr) | Procédé pour optimiser l'exploitation d'un milieu souterrain au moyen d'une etude réservoir comportant une mise à l'échelle optimisée | |
| EP2806101B1 (fr) | Procédé d'exploitation d'un milieu fracture à partir d'un modèle de réservoir calé pour des puits choisis au moyen d'un modèle de transmissivité équivalente | |
| EP4112875B1 (fr) | Procédé pour simuler les écoulements d'un fluide dans une formation souterraine comportant un réseau de fractures | |
| FR3037682A1 (fr) | ||
| WO2021063804A1 (fr) | Procede pour determiner une trajectoire d'un puits dans un reservoir petrolier |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| CD | Change of name or company name | ||
| PLFP | Fee payment |
Year of fee payment: 10 |
|
| PLFP | Fee payment |
Year of fee payment: 11 |
|
| PLFP | Fee payment |
Year of fee payment: 12 |
|
| PLFP | Fee payment |
Year of fee payment: 14 |
|
| PLFP | Fee payment |
Year of fee payment: 15 |
|
| PLFP | Fee payment |
Year of fee payment: 16 |
|
| PLFP | Fee payment |
Year of fee payment: 17 |
|
| PLFP | Fee payment |
Year of fee payment: 18 |