FR 2918179 A1 20090102 FR 0704703 A 20070629 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 . 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 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. REVENDICATIONS 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.