FR 3041026 A1 20170317 FR 1558653 A 20150915 1 La présente invention concerne le domaine de l'industrie pétrolière, et plus particulièrement l'exploration et l'exploitation de gisements d'hydrocarbures ou de sites de stockage géologique de gaz. En particulier, l'invention concerne un procédé pour construire une représentation d'une formation souterraine traversée par un réseau de fractures quelconque, et l'utilisation de cette représentation pour la simulation des écoulements de fluide dans la formation considérée. L'exploration et l'exploitation de gisements, notamment pétroliers, nécessitent d'acquérir une connaissance aussi parfaite que possible de la géologie souterraine, et ce, afin de fournir de façon efficace une évaluation des réserves, une modélisation de la production, ou de la gestion de l'exploitation. En effet, la détermination de l'emplacement d'un puits de production ou d'un puits d'injection au sein d'un gisement d'hydrocarbures, la constitution de la boue de forage, les caractéristiques de complétion, le choix d'un procédé de récupération des hydrocarbures (tel que l'injection d'eau par exemple) et des paramètres nécessaires à la mise en oeuvre de ce procédé (tels que la pression d'injection, le débit de production,...) nécessitent de bien connaltre le gisement. La connaissance d'un gisement signifie disposer d'une description aussi précise que possible de la structure, des propriétés pétrophysiques, des propriétés des fluides, etc, du gisement étudié. Pour acquérir cette connaissance, l'industrie pétrolière allie les mesures sur champ (réalisées in situ, lors de campagnes sismiques, de mesures dans des puits, de carottages etc) aux modélisations expérimentales (réalisées au laboratoire) ainsi qu'aux simulations numériques (réalisées au moyen de logiciels). La formalisation de cette connaissance passe ensuite par l'établissement d'une maquette du sous-sol, connue sous le terme de « modèle géologique », qui permet de rendre compte de ces aspects de façon approchée. Généralement, ce type de maquette est représenté sur un ordinateur, et l'on parle alors de modèle numérique. Un modèle de réservoir comporte une représentation maillée (sous la forme d'une grille régulière, ou sous la forme d'un maillage de façon plus générale), généralement tridimensionnelle, au moins une propriété pétrophysique (porosité, perméabilité, saturation...) étant attribuée à chacune des mailles de la représentation maillée. 3041026 2 Afin de reproduire ou prédire (i.e. "simuler") la production d'hydrocarbures réelle, le spécialiste en ingénierie de réservoir met en oeuvre un logiciel de calcul, appelé « simulateur de réservoir ». Le simulateur de réservoir est un simulateur d'écoulement, qui calcule les écoulements et l'évolution des pressions au sein du réservoir représenté 5 par un « modèle de réservoir ». Si la puissance informatique disponible pour réaliser les simulations d'écoulement le permet, le modèle de réservoir peut se confondre avec le modèle géologique. Dans le cas contraire, le modèle de réservoir peut être obtenu à l'issue d'une technique d'« upscaling » (mise à l'échelle), qui permet de passer du modèle géologique (modèle aux mailles plus fines) au modèle réservoir (modèle aux mailles plus 10 grossières). Cette étape d'upscaling est bien connue du spécialiste en ingénierie de réservoir et peut être réalisée par exemple à l'aide du logiciel CobraFlow (IFP Energies nouvelles, France). Les résultats de ces calculs permettent alors de prévoir et d'optimiser dès schémas d'exploitation (définition du nombre de puits à implanter, de leur position, du mode de 15 récupération assistée etc) du gisement étudié afin d'améliorer les débits et/ou les quantités d'hydrocarbures récupérés. Le calcul du comportement du réservoir suivant un scénario de production donné constitue une simulation de réservoir. Les documents suivants seront cités au cours de la description : 20 Bourbiaux, B., Cacas, M. C., Sarda, S. and Sabathier J. C., 1998, "A Rapid and Efficient Methodology to Convert Fractured Reservoir Images into a Dual-Porosity Model", Oil & Gas Science and Technology, Vol. 53, No. 6, Nov.-Déc. 1998, 785-799. Delorme, M., Atfeh, B., Allken, V. and Bourbiaux, B., 2008, Upscaling Improvement 25 for Heterogeneous Fractured Reservoir Using a Geostatistical Connectivity Index, edited in Geostatistics 2008, VIII International Geostatistics Congress, Santiago, Chile. Oda, M., 1985, Permeability tensor for discontinuous Rock Masses, Geotechnique Vol 35, Issue 4, 483-495. Warren, J.E. et Root, P.J., "The Behavior of Naturally Fractured Reservoirs", SPE 30 Journal (septembre 1963), 245-255. Fourno, A., Grenier, C., Benabderrahmane, A., Delay, F., 2013, A continuum voxel approach to model flow in 3D fault networks: A new way to obtain up-scaled hydraulic conductivity tensors of grid cells, Journal of Hydrology, 493, 68-80. 3041026 3 Les réservoirs fracturés constituent un type extrême de réservoirs géologiques hétérogènes, comportant deux milieux contrastés: un milieu matriciel contenant la plus grande part de l'huile en place et présentant une faible perméabilité, et un milieu fracturé représentant moins de 1 % de l'huile en place et hautement conducteur. Le milieu fracturé 5 lui-même peut être complexe, composé de différentes familles de fractures caractérisées par leur densité, et des distributions portant sur la longueur des fractures de chaque famille, leur orientation dans l'espace et leur ouverture. On appelle « fracture », une discontinuité plane, de très faible épaisseur par rapport à son extension, et qui représente un plan de rupture d'une roche du gisement. 10 Les spécialistes en charge de l'exploitation de réservoirs fracturés, ont besoin de parfaitement connaIre le rôle des fractures. D'une part, la connaissance de la distribution et du comportement de ces fractures permet d'optimiser les schémas de production d'un gisement d'hydrocarbures, c'est-à-dire le nombre de puits à forer, leur localisation, leur géométrie, le type de fluide de récupération injecté etc. D'autre part, la géométrie du 15 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Êre la distribution des fractures est également très utile pour l'ingénieur de réservoir qui cherche à calibrer les modèles qu'il construit, de sorte à ce que les modèles construits reproduisent les courbes de production passées, et ce, afin de 20 prédire les productions futures de façon fiable. A ces fins, les spécialistes de géosciences disposent d'images tridimensionnelles des gisements, permettant de localiser un grand nombre de fractures. Confronté à la complexité des milieux fracturés, les spécialistes utilisent bien 25 souvent une approche de type « double milieu » pour représenter ce type de milieu. Cette approche, décrite par exemple dans (Warren J.E. et al., 1963), consiste à considérer le milieu fracturé comme deux continuums échangeant des fluides entre eux : des blocs matriciels et des fractures. On parle alors de modèle « double milieu » ou « double porosité ». Ainsi, la modélisation en « double milieu » d'un gisement fracturé consiste à 30 discrétiser ce gisement en deux ensembles de mailles (appelés grilles) superposés, constituant la grille "fissure" et la grille "matrice". Chaque volume élémentaire du gisement fracturé est ainsi conceptuellernent représenté par deux mailles, l'une "fissure" et l'autre "matrice", couplées entre elles (c'est-à-dire échangeant des fluides). Dans la réalité du champ fracturé, ces deux mailles représentent l'ensemble des blocs matriciels délimités par des fractures présentes en cet endroit du réservoir. En effet, le plus souvent, les mailles ont des dimensions latérales hectométriques (couramment 100 ou 200 m) compte 3041026 4 tenu de la taille des champs et des possibilités limitées des logiciels de simulation en termes de capacité et temps de calcul. Il en résulte que, pour la plupart des champs fracturés, le volume élémentaire (maille) de réservoir fracturé renferme d'innombrables fractures formant un réseau complexe délimitant de multiples blocs matriciels de 5 dimensions et formes variables suivant le contexte géologique. Chacun des blocs réels constitutifs échange des fluides avec les fractures qui l'entourent à un rythme (débit) qui lui est propre car dépendant des dimensions et de la forme de ce bloc particulier. Face à une telle complexité géométrique du milieu réel, la démarche du spécialiste en ingénierie de réservoir consiste, pour chaque volume élémentaire (maille) de réservoir, 10 à représenter le milieu fracturé réel comme un ensemble de blocs matriciels tous identiques, parallélépipédiques, délimités par un réseau orthogonal et régulier de fractures orientées suivant les directions principales d'écoulement : pour chaque maille, on détermine ainsi les perméabilités dites "équivalentes" de ce réseau de fractures et on définit un bloc matriciel, dit « représentatif » (de la distribution réelle (géologique) des 15 blocs), unique et de forme parallélépipédique. Il est alors possible de formuler et calculer les flux d'échange matrice-fissure pour ce bloc "représentatif', et d'en multiplier le résultat par le nombre de tels blocs dans le volume élémentaire (maille) pour obtenir le flux à l'échelle de cette maille. Il convient cependant de noter que le calcul des perméabilités équivalentes 20 nécessite de connaître les propriétés d'écoulement (Le. les conductivités) des fractures discrètes du modèle de réservoir. C'est pourquoi, préalablement à la construction de ce modèle de réservoir équivalent (dit "modèle de réservoir double milieu") tel que décrit plus haut, il est d'abord nécessaire de simuler les réponses en écoulement de quelques puits (tests d'écoulement transitoires ou pseudo-permanents dits aussi tests de puits, tests 25 d'interférences, tests de débit, etc.) sur des modèles extraits du modèle de réservoir donnant une représentation discrète (réaliste) des fractures alimentant ces puits. L'ajustement des réponses en pression/débit simulées sur les mesures de champ permet de calibrer les conductivités des familles de fractures. Bien que ne couvrant qu'une aire limitée (aire de drainage) autour du puits, un tel modèle de simulation de tests de puits 30 comporte encore de très nombreux noeuds de calcul si le réseau de fractures est dense. En conséquence, la taille des systèmes à résoudre et/ou la durée des calculs demeurent souvent prohibitives. 35 3041026 5 Etat de la technique On connait la demande de brevet FR 2967200 (US 8983818) qui décrit un procédé permettant une simplification des réseaux de fractures à l'échelle locale de l'aire de 5 drainage du puits, afin de pouvoir simuler, en un temps de calcul raisonnable, les tests de puits de réservoirs fracturés et calibrer ainsi les conductivités des familles de fractures. Cette calibration hydraulique des fractures aboutit à un jeu de paramètres caractérisant le réseau de fractures (ou modèle de fractures). Ce modèle de fractures est par la suite utilisé pour construire un modèle d'écoulement double milieu à l'échelle du réservoir. 10 Toutefois, cette méthode est uniquement adaptée aux réseaux de fractures sensiblement orthogonales aux couches et continus d'une couche à une autre, c'est-à-dire à des réseaux de fractures dont le comportement de l'écoulement est assimilable à un comportement 2D (cf notamment les demandes de brevet EP 2530493 (US 9103194) et EP 2581767 (US 13/644,479), ainsi que le document (Fourno et al., 2013)). Un exemple 15 d'un tel réseau de fractures est donné en figure 2a. Dans cet exemple, la formation comporte deux couches, le réseau de fractures est continu d'une couche à une autre (assurant ainsi la connectivité d'une couche à une autre) et les fractures sont sensiblement orthogonales aux limites des couches. 20 On conne également les demandes de brevet EP 2530493 (US 9103194) et EP 2581767 (US 2013/0096889) qui permettent de prendre en compte des réseaux de fractures qui ne sont pas nécessairement sub-orthogonales et/ou continues d'une couche à une autre, tout en limitant le nombre de noeuds de calcul. Pour y parvenir, la demande EP 2530493 (US 9103194) est basée sur l'utilisation d'un diagramme de Voronoï sur 25 chaque plan de fracture, afin de construire un maillage du milieu fracturé qui donne une bonne estimation des volumes en présence à chaque noeud et des lignes de courant représentant le déplacement des fluides de noeuds à noeuds. Ainsi, les maillages ne comportent que 3 ou 4 noeuds aux intersections, permettant de traiter des millions de fractures en un temps de calcul raisonnable. La demande EP 2581767 (US 30 2013/0096889) permet de s'affranchir de l'hypothèse restrictive d'avoir autant de noeuds fracture que de noeuds matrice tout en préservant les volumes en place et la physique d'écoulement grâce à un schéma décentré. L'utilisation des structures type arbre (et en particulier de l'octree) permet d'accélérer la construction du maillage du milieu fracture (limitation du nombre de calcul d'intersections), de simplifier les calculs de termes 35 d'échange matrice/fracture ainsi que de bien estimer les volumes matrice associées aux 3041026 6 noeuds fractures. Cette méthode, appliquée ici au maillage d'un milieu dual en 3D, s'applique à tout problème dual impliquant de fortes hétérogénéités. Ainsi, les demandes EP 2530493 (US 9103194) et EP 2581767 (US 2013/0096889) permettant de représenter les réseaux de fractures non nécessairement sub-orthogonales et/ou continues d'une 5 couche à une autre, des simulations d'écoulement 3D sur ces réseaux 3D de fractures, réalisées à l'aide d'un simulateur 3D d'écoulement, ont permis d'accéder à une meilleure connaissance des écoulements 3D dans les gisements étudiés. Toutefois, ces méthodes restent aujourd'hui prohibitives en termes de temps de calcul, notamment lorsqu'un nombre important de simulations d'écoulement (à l'échelle du réservoir ou à l'échelle 10 d'une zone d'influence d'un ou plusieurs puits) est nécessaire (cas de l'optimisation de schéma de production, du calage d'historique de production, ou encore de la calibration des paramètres statistiques des familles de fractures par des procédures de sensibilité aux paramètres). Un des objets de la présente invention concerne un procédé pour représenter de 15 façon appropriée une formation souterraine comportant un réseau de fractures dont au moins une famille peut être caractérisée par des pendages quelconques et/ou par une continuité variable des fractures d'une couche à une autre. En particulier, la présente invention établit un paramètre permettant de qualifier le type de fractures et adapte, en fonction de ce paramètre, la simplification du milieu fracturé. De plus, cette simplification 20 n'est pas localisée à l'échelle locale de l'aire de drainage du puits, mais peut directement être réalisée à l'échelle du réservoir. Ce modèle de fractures simplifié, adapté à la nouvelle génération de simulateur de réservoir 3D, est ensuite utilisé pour optimiser les schémas de production du fluide contenu dans la formation étudiée, et ce, en des temps de calculs abordables. 25 Le procédé selon l'invention Ainsi, la présente invention concerne un procédé pour l'exploitation d'un fluide au 30 sein d'une formation souterraine comprenant au moins une couche traversée par un réseau de fractures et par au moins un puits, dans lequel on construit une représentation maillée de ladite au moins une _couche à partir de mesures de propriétés relatives à ladite au moins une couche, et on définit, au moins pour ladite couche, une zone d'intérêt regroupant un ensemble de mailles de ladite représentation maillée relative à ladite couche et comprenant 3041026 7 au moins ledit puits. Le procédé comporte au moins les étapes suivantes pour chacune desdites couches : A. à partir desdites mesures, on caractérise ledit réseau de fractures par des paramètres statistiques, et on construit, à partir desdits paramètres, un modèle 5 de réseau de fractures discret; B. pour chacun desdits puits contenus dans ladite zone d'intérêt, on définit au moins trois zones de simplification (ZS1, ZS2, ZS3), une première zone de simplification (ZS1) contenant ledit puits, une seconde zone de simplification (ZS2) ayant comme frontière interne la frontière externe de ladite première 10 zone (ZS1), et une troisième zone de simplification (ZS3) ayant comme frontière interne la frontière externe de ladite seconde zone (ZS2), l'ensemble desdites zones de simplification de l'ensemble desdits puits couvrant l'ensemble des mailles de ladite zone d'intérêt ; C. à partir d'au moins lesdits paramètres statistiques, on détermine, pour au 15 moins chacune desdites zones de simplification (ZS1, ZS2, ZS3) de ladite zone d'intérêt, un tenseur de perméabilité équivalente et une valeur d'un paramètre P caractérisant l'orientation et la continuité vers les autres couches des fractures dudit réseau ; D. on détermine une simplification dudit modèle en chacune desdites zones de 20 simplification, en fonction de ladite zone, dudit tenseur de perméabilité équivalente et de ladite valeur dudit paramètre P ; Puis, à partir de ladite représentation maillée, dudit modèle de réseau de fractures simplifié, et d'un simulateur d'écoulement, on définit un schéma d'exploitation optimal dudit fluide de ladite formation et on exploite ledit fluide de ladite formation en fonction dudit schéma 25 d'exploitation optimal. Avantageusement, lesdits paramètres statistiques peuvent être choisis parmi les paramètres suivants : densité de fractures, longueur de fractures, orientation des fractures dans l'espace, ouverture de fractures, et distribution de fractures au sein de ladite couche. 30 Préférentiellement, lesdites zones de simplification peuvent être des cylindres d'axe vertical et de section elliptique, centrées autour dudit puits. Selon un mode de réalisation de la présente invention, le grand axe d'une desdites zones de simplification peut être orienté selon une direction préférentielle d'écoulement, 3041026 8 ladite direction préférentielle étant déterminée à partir d'une diagonalisation dudit tenseur de perméabilité équivalente. Avantageusement, la distance séparant deux desdites zones peut être déterminée en fonction d'une progression géométrique de raison constante. 5 Selon un mode de réalisation de la présente invention, on peut ne réaliser aucune simplification dudit modèle de réseau de fractures dans ladite première zone de simplification (ZS1), réaliser une première simplification dans ladite seconde zone (ZS2), et une seconde simplification dans ladite troisième zone (ZS3), ladite seconde simplification étant plus importante que ladite première simplification. 10 Selon un mode de mise en oeuvre de l'invention, lesdites seconde et/ou troisième zones de simplification (ZS2, ZS3) peuvent être découpées en sous-zones de simplification annulaires. Avantageusement, lesdites secondes zones de simplification (ZS2) peuvent être découpées en sous-zones de simplification en colonne. 15 Préférentiellement, lesdites troisièmes zones de simplification (ZS3) peuvent être découpées en sous-zones de simplification localement homogènes. Selon un mode de mise en oeuvre de l'invention, on peut déterminer ledit paramètre P pour au moins chacune desdites zones de simplification (ZS1, ZS2, ZS3), au moyen d' : - une valeur d'un paramètre PO mesurant l'orientation dans l'espace desdites 20 fractures de ladite couche selon une formule du type : dev2 dey 180 PO = - 27-0 + 1 avec dey = arccos(dtr3. * (Pli, et/ou 904 - une valeur d'un paramètre PC mesurant la continuité desdites fractures de ladite couche vers les autres couches selon une formule du type : k3 25 où kl, k2, et k3 sont les valeurs propres dudit tenseur de perméabilité après diagonalisation. Selon un mode de réalisation de la présente invention, ladite simplification peut consister en le calcul d'au moins une valeur d'espacement entre lesdites fractures, une ouverture desdites fractures, et une valeur de perméabilité. 30 Selon un mode de mise en oeuvre de l'invention, ladite simplification dudit modèle peut être réalisée selon trois directions de l'espace si le paramètre P caractérise des PC = (kl+ke 3041026 9 fractures non sensiblement orthogonales aux couches et/ou continues vers les autres couches. Selon un autre mode de mise en oeuvre de l'invention, ladite simplification dudit modèle peut être réalisée selon deux directions de l'espace si le paramètre P caractérise 5 des fractures sensiblement orthogonales aux couches et continues vers les autres couches. Préférentiellement, ledit schéma d'exploitation optimal peut être défini en déterminant un procédé de récupération dudit fluide, ainsi qu'un nombre, une implantation et une géométrie de puits injecteurs et producteurs permettant de satisfaire des critères 10 technico-économiques prédéfinis. Avantageusement, ladite exploitation dudit fluide en fonction dudit schéma d'exploitation optimal peut consister au moins à forer lesdits puits injecteurs et producteurs et à produire ledit fluide selon ledit procédé de récupération. En outre, l'invention concerne un produit programme d'ordinateur téléchargeable 15 depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, comprenant des instructions de code de programme pour la mise en oeuvre du procédé selon la description ci-dessus, lorsque ledit programme est exécuté sur un ordinateur. 20 D'autres caractéristiques et avantages du procédé selon l'invention, apparaltront à 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. 25 Présentation succincte des figures - La figure 1 présente un exemple de réseau de fractures selon une section horizontale réalisée dans une couche d'une formation souterraine. - La figure 2a illustre un exemple de réseau de fractures 2.5D dans le cas d'une 30 formation souterraine composée de deux couches. - La figure 2b présente un exemple de réseau de fractures 3D dans le cas d'une formation souterraine composée de quatre couches. 3041026 10 La figure 2c illustre un exemple de réseau de fractures 3D dans le cas d'une formation souterraine composée d'une seule couche. La figure 3 illustre la création de zones et sous-zones de simplification autour d'un puits traversant une couche d'une formation souterraine. 5 - La figure 4 illustre un réseau de fractures simplifié DFNs sur une grille fine. - La figure 5 illustre un modèle de réseau de fractures simplifié DFNs selon l'invention. La figure 6 présente le recouvrement entre modèle de réseaux de fractures simplifiés DFNs de sous-zones de simplification voisines. - La figure 7 illustre un modèle de réseau de fractures simplifié DFNs résultant de la 10 mise en oeuvre de la présente invention dans le cas d'une couche d'une formation souterraine traversée par cinq puits. Description détaillée du procédé 15 Les définitions suivantes sont utilisées au cours de la description de l'invention : - Tests de puits : il s'agit d'une injection ou d'une production d'un fluide dans le sous-sol au niveau d'un puits donné. Le test de puits est caractérisé par une courbe de production (quantité de fluide injecté ou produit) et d'une courbe de pression (courbe qui caractérise l'évolution de la pression dans la 20 structure en fonction de la courbe de production). Les tests de puits sont utiles pour calibrer les modèles d'écoulement/de réservoir utilisés pour les simulations d'écoulement. Tests d'interférence : il s'agit d'une association de plusieurs tests de puits. Chaque puits est caractérisé par sa courbe de production et sa courbe de 25 pression. La différence avec un test de puits unique par puits réside dans le fait que chaque puits peut être perturbé par le fonctionnement d'un autre. ... Le procédé selon l'invention concerne l'exploitation d'un fluide au sein d'une formation souterraine comprenant au moins une couche traversée par un réseau de 30 fractures et par au moins un puits. De façon non limitative, la formation souterraine est constituée d'au moins une couche de roche réservoir (c'est-à-dire d'une roche au moins poreuse et préférentiellement perméable). L'exploitation du fluide est alors appelée 3041026 11 exploitation des ressources, et vise à extraire le fluide de cette roche réservoir. Non limitativement, le fluide en question est de type hydrocarbures (huile, gaz, bitume etc). Une formation souterraine peut être constituée de plusieurs couches géologiques, chaque couche géologique pouvant comprendre un réseau de fractures propre, et chaque 5 réseau pouvant être constitué de plusieurs familles de fractures. Par la suite et à des fins de simplification, on appellera « réseau de fractures 3D » un réseau de fractures constitué d'au moins une famille de fractures de pendages (orientations) quelconques (figure 2c) et/ou caractérisée par une continuité variable des fractures d'une couche à une autre (figure 2b) ; et on appellera « réseaux de fractures 2.5D » un réseau de fractures 10 constitué de familles de fractures sensiblement perpendiculaires aux couches et dont les fractures sont continues d'une couche à une autre (figure 2a). La présente invention concerne la prise en compte des réseaux de fractures quelconques dans une simulation d'écoulement, et ce en des temps de calcul raisonnables, même dans le cas d'un réseau de fractures 3D. Pour ce faire, la présente 15 invention consiste à modéliser le réseau de fractures de la couche considérée (voir l'étape 1 décrite ci-dessous), puis à le simplifier (voir l'étape 2 décrite ci-dessous), afin de permettre des simulations numériques d'écoulement (voir l'étape 3 décrite ci-dessous) en des temps de calcul raisonnables. 20 La présente invention requiert de disposer : - de mesures de propriétés relatives au moins à la couche étudiée : il peut s'agir de mesures de propriétés pétrophysiques réalisées in situ, en différents points de la formation étudiée, telles que la porosité, la perméabilité, et la lithologie (c'est-à-dire le 25 type de roche), la perméabilité relative ou encore la pression capillaire. Ces mesures peuvent avoir été obtenues par exemple par carottage, via des diagraphies réalisées dans des puits, par des campagnes d'acquisition sismique. Mais il s'agit aussi de mesures (débits d'huile, débits d'eau, variations de pression par exemple) liées aux écoulements dans la couche étudiée, obtenues par exemple par des mises en production du fluide en 30 certains puits traversant la formation, au cours de tests de puits ou encore de tests d'interférence ; - d'une représentation maillée représentative d'au moins la couche étudiée : appelée aussi « modèle de réservoir », il s'agit d'une sorte de maquette du sous-sol construite dans le but de décrire aussi précisément que possible la structure, les 35 propriétés pétrophysiques, les propriétés des fluides, de la formation étudiée, et d'au 3041026 12 moins la couche étudiée. Cette maquette est généralement représentée sur un ordinateur, et consiste en un maillage ou grille, chacune des mailles de cette grille comportant une ou plusieurs valeurs de propriétés relatives à la formation étudiée (porosité, perméabilité, saturation, etc). Un modèle de réservoir se doit de vérifier autant que possible les 5 propriétés relatives à la couche étudiée collectées sur le terrain : les données de diagraphie mesurées le long des puits, les mesures réalisées sur des échantillons de roche prélevés dans les puits, les données déduites de campagnes d'acquisition sismique, les données de production comme les débits d'huile, d'eau, les variations de pression etc. Le spécialiste en simulation de réservoir a pleine connaissance de 10 méthodes pour construire une telle représentation maillée d'une formation. A noter que le modèle de réservoir peut se confondre avec le modèle géologique lorsque la puissance informatique est suffisante pour permettre des calculs numériques de simulation d'écoulement sur une grille à mailles fines. Dans les autres cas, le spécialiste pourra avoir recours à une technique d' « upscaling » (mise à l'échelle) afin de passer d'un modèle aux 15 mailles fines (le modèle géologique) à un modèle aux mailles plus grossières (le modèle de réservoir). Cette étape d'upscaling peut être réalisée par exemple à l'aide du logiciel CobraFlow (IFP Energies nouvelles, France). - une zone d'intérêt (Z01), regroupant un ensemble de mailles de ladite représentation maillée relatives à la couche étudiée et comprenant au moins le puits 20 traversant la couche concernée : les dimensions et la géométrie de cette zone peuvent être variables selon le mode de mise en oeuvre de la présente invention considéré. Dans le cas le plus général, c'est-à-dire lors de simulations de production du fluide de la formation étudiée mises en oeuvre pour tester différents schémas d'exploitation, la ZOI peut correspondre à l'ensemble du réservoir, et peut comprendre au moins deux puits, et 25 préférentiellement l'ensemble des puits. Dans le cas où préalablement à cette simulation, on procédera à une calibration du modèle utilisé pour la simulation, via des tests de puits, la ZOI considérée correspondra avantageusement à une partie du réservoir, d'une dimension kilométrique, incluant un seul puits. Dans le cas de tests d'interférence, la ZOI peut alors correspondre à une portion du réservoir, et comprend plusieurs puits. A noter 30 que la ZOI peut être discontinue et être par exemple formée d'une somme de sous-zones d'intérêt définies localement autour d'un nombre de donné de puits. - d'un simulateur d'écoulement : un simulateur d'écoulement est un programme numérique, exécuté sur un ordinateur, qui est utilisé pour prédire l'écoulement de fluides au sein des différentes couches d'une formation. La simulation d'écoulement, dite aussi 35 simulation de réservoir, consiste à prédire numériquement la production d'un fluide piégé dans une ou plusieurs couches d'une formation souterraine, la production nécessitant 3041026 13 l'existence d'au moins un puits injecteur (dans lequel on injecte un autre fluide, qui va chasser le fluide piégé) et d'un puits producteur (vers lequel le fluide piégé, chassé va se déplacer et duquel il va pouvoir être extrait). Avantageusement, on aura recours à un simulateur d'écoulement dit « double milieu », qui permet de modéliser les échanges 5 entre les blocs de roches non fracturés et le réseau de fracture, sans -nécessiter la discrétisation de ces blocs. Un exemple d'un tel simulateur est le logiciel PUMAFLOW (I FP Energies nouvelles, France). 10 Le procédé selon l'invention comporte au moins les étapes suivantes, les étapes 1 et 2 étant répétées pour chacune des couches de la formation étudiée : 1. Modélisation du réseau de fractures 1.1. Détermination de paramètres statistiques 15 1.2. Construction d'un modèle de réseau de fractures discret 2. Simplification du modèle de réseau de fractures 2.1. Définition de zones de simplification 2.2. Détermination d'un tenseur de perméabilité équivalente 2.3. Détermination d'un paramètre P caractérisant la géométrie des fractures 20 2.4. Simplification du modèle de réseau de fractures discret 3- Exploitation du fluide de la formation Les principales étapes de la présente invention sont détaillées ci-après. Les étapes 25 1 et 2 seront déclinées pour une seule couche de la formation étudiée et sont à répéter pour l'ensemble des couches de la formation. La mise en oeuvre de l'invention nécessite l'existence d'au moins un puits traversant chacune des couches considérées, mais la mise en oeuvre de l'invention ne nécessite pas de prendre en compte tous les puits traversant une couche donnée. 30 3041026 14 1. Modélisation du réseau de fractures Au cours de cette étape, il s'agit de déterminer un modèle du réseau de fractures de la couche étudiée, respectant les mesures de propriétés relatives à la couche réalisées in situ. 5 La réalisation de cette étape peut comporter au moins les sous-étapes suivantes : 1.1. Détermination de paramètres statistiques Au cours de cette sous-étape, il s'agit de déterminer des paramètres statistiques permettant de caractériser le réseau de fractures de la couche étudiée. Cette sous-étape 10 requiert des mesures, permettant une caractérisation directe ou indirecte de la fracturation. Des informations sur la fracturation d'une couche géologique peuvent être obtenues via des : carottes, extraites de la formation étudiée, et à partir desquelles une étude 15 statistique des fractures intersectées peut être effectuée ; affleurements, qui présentent l'avantage de fournir une vision à grande échelle du réseau de fractures ; images sismiques, permettant d'identifier de grands événements géologiques. 20 Selon l'invention, on utilise des paramètres statistiques (notés PSF par la suite) afin de caractériser les fractures observées lors des mesures. Selon un mode de mise en oeuvre de l'invention, on caractérise la densité des fractures observées, leur longueur, leur orientation dans l'espace, leur ouverture, et leur distribution au sein de la couche étudiée. 25 On dispose ainsi de paramètres statistiques (PSF) décrivant le réseau de fractures de la couche étudiée, à partir desquels un modèle de ce réseau peut être construit, à l'échelle de chacune des mailles de la représentation maillée. 1.2. Construction d'un modèle de réseau de fractures discret 30 Selon l'invention, on construit ensuite une image réaliste du réseau de fractures caractérisé par les paramètres statistiques (PSF) déterminés à la sous-étape précédente 3041026 15 réelles, en passant par un modèle de réseau de fractures discret (connu sous l'acronyme DFN, pour « Discrete Fracture Network » en anglais). Pour ce faire, à partir d'une représentation maillée de la formation étudiée (et au moins de la couche considérée), on associe en chaque maille, une représentation 5 détaillée (DFN) de la complexité interne du réseau de fractures, aussi fidèle que possible des observations directes et indirectes du réservoir. La figure 1 présente une section horizontale prise dans une couche d'une formation, et le réseau de fractures, observable à l'échelle de la formation, pour cette section. La figure 2 présente le réseau discret de fractures visibles à l'échelle de la maille tridimensionnelle dont les limites 2D sont 10 présentées par le carré en figure 1. On peut observer en figure 2 un ensemble de blocs matriciels poreux, de formes et tailles irrégulières, délimités par des fractures. Ce réseau de fractures discret constitue une image représentative du réseau réel de fractures délimitant les blocs matriciels. Pour construire un modèle de réseau de fractures discret (DFN) en chaque maille 15 de la représentation maillée, on peut utiliser des logiciels de modélisation, bien connus des spécialistes, tel que le logiciel FRACAFlowe (IFP Energies nouvelles, France). Ces logiciels utilisent les paramètres statistiques (PSF) déterminés précédemment. 20 2. Simplification du modèle de réseau de fractures En raison de son extrême complexité géométrique, le modèle de réseau de fractures (DFN) obtenu à l'étape précédente et représentatif du réservoir fracturé réel, ne peut être utilisé pour simuler, i.e. reproduire et/ou prédire, les écoulements locaux autours des puits. En effet, le DFN obtenu à l'étape précédente présente bien souvent un grand 25 nombre de fractures et de nombreux points d'intersection entre fractures. Le maillage sera donc trop complexe ou présentera trop de noeuds de calcul, rendant les simulations inabordables en temps de calcul. L'objectif de cette étape est de déterminer un modèle de réseau de fractures simplifié, respectant néanmoins les mesures d'écoulement réalisées in situ. 30 Ainsi, au cours de cette étape, il s'agit de simplifier le modèle de réseau de fractures discret déterminé à l'étape précéderge pour la couche considérée. Selon l'invention, cette simplification est différente en fonction des zones de simplification décrites autour de chacun des puits considérés pour la couche en question (voir l'étape 2.1). Chaque zone de simplification est caractérisée par 2 paramètres qui lui sont propres : un tenseur de 35 perméabilité équivalente (voir l'étape 2.2) et un paramètre (P) caractérisant à la fois 3041026 16 l'orientation et la continuité (vers les autres couches de la formation étudiée) du réseau de fractures de cette couche (voir l'étape 2.3). 2.1. Définition de zones de simplification 5 Au cours de cette sous-étape, il s'agit de définir des zones de simplification pour chacun des puits considérés. Selon l'invention, on définit autour de chaque puits, au moins trois zones de simplification (ZS1, ZS2, ZS3), la première desdites zones de simplification (ZS1) contenant ledit puits, la seconde desdites zones de simplification (ZS2) ayant comme frontière interne la frontière externe de ladite première zone (ZS1), et 10 la troisième desdites zones de simplification (ZS3) ayant comme frontière interne la frontière externe de ladite seconde zone (ZS2). De plus, selon l'invention, l'ensemble des zones de simplification pour l'ensemble des puits considérés doit couvrir l'ensemble des mailles de la ZOI, et ce quel que soit le nombre de puits. L'invention n'est pas limitée à la définition de trois zones de 15 simplification pour chacun des puits considérés, et une mise en oeuvre de l'invention peut être réalisée pour NZ zones de simplification, avec NZk3. Selon un mode préférentiel de mise en oeuvre de l'invention, les zones de simplification correspondent à des cylindres de section elliptique, l'axe de révolution des cylindres pouvant être orienté selon un axe vertical (voir par exemple les sections 20 horizontales des zones ZS1, ZS2 et ZS3 créées autour d'un puits W en figure 3). Les sections dans un plan horizontal des zones de simplification sont alors des ellipses (ZS1) ou des couronnes elliptiques (ZS2 et ZS3) concentriques, centrées autour d'un puits. Selon un autre mode de réalisation de l'invention, les zones de simplification correspondent à des cylindres de section elliptique, l'axe de révolution des cylindres 25 pouvant être orienté selon un axe défini par la perpendiculaire au pendage (ligne de plus grande pente) de la couche considérée, ou encore un axe défini par l'axe du puits considéré. Selon un mode de mise en oeuvre de l'invention, une zone de simplification ZS2 et/ou une zone ZS3 est découpée en N sous-zones de simplification en rajoutant, entre 30 les frontières internes et externes de la zone considérée, N-1 frontières intermédiaires, la géométrie de la frontière intermédiaire n (avec n compris entre 1 et N-1) pouvant résulter d'une combinaison linéaire (par exemple fonction de n et N) des courbes décrivant les frontières internes et externes de la zone de simplification considérée. Par exemple, dans le cas d'une zone de simplification cylindrique de section elliptique, le découpage de cette 35 zone de simplification en deux sous-zones de simplification revient à rajouter une frontière 3041026 17 cylindrique intermédiaire (voir par exemple la frontière ann_ZS3 en figure 3), de section elliptique. On parlera par la suite de « sous-zone annulaire » pour ce type de sous-zones de simplification. Ce mode de mise en oeuvre sera préféré dans le cas de propriétés d'écoulement relativement homogènes. 5 Selon un mode de réalisation dans lequel les zones de simplification ou sous-zones de simplification sont des cylindres de section horizontale à géométrie elliptique, on peut très préférentiellement orienter les ellipses (direction du grand axe de chaque ellipse) en fonction des directions d'écoulement préférentielles déduites d'un calcul de perméabilités équivalentes. Un tel calcul est précisément décrit à la sous-étape suivante. De même, le 10 rapport de forme de l'ellipse (en notant Lmax la demi-longueur du grand axe de l'ellipse, et Lm, la demi-longueur du petit axe de l'ellipse, ce rapport de forme est défini par Lmax/Lmin) et la distance entre ellipses concentriques (distance séparant les limites elliptiques intérieure et extérieure d'une couronne concentrique donnée) peuvent être déterminées en fonction de l'anisotropie d'écoulement dans la couche étudiée. 15 Selon un mode particulier de mise en oeuvre de l'invention, on peut définir le rapport de forme selon une formule du type : Lmax _ iimax(ki,k2) Lm~nmin(ki,k2) où k1 et k2 sont respectivement les valeurs propres, obtenues après diagonalisation du tenseur de perméabilité équivalente, dans les directions principales les plus proches d'un 20 plan horizontal. Selon un autre mode particulier de mise en oeuvre de l'invention, afin de définir la distance entre ellipses concentriques (frontières entre zones), on opte pour un critère de dimensionnement conforme à une précision de modélisation uniformément distribuée sur le domaine de simulation, qui consiste à fixer les demi-longueurs du grand axe Lmax (ou du 25 petit axe Lin) de 2 ellipses successives i et i+1 (avec i variant de 1 à I, et I au moins égal à 3) à des valeurs en progression géométrique de raison r constante (égale à 2 par exemple), selon une formule du type : Lmax(i+1)= Lmin(i+lj r Lmax(i) Ltnin(i) Cette règle permet de donner un poids égal à chaque zone i (i=1 à I, en termes de 30 différence de pression constatée sur chaque couronne en régime d'écoulement permanent. 3041026 18 Selon un mode préféré de mise en oeuvre de l'invention, une zone de simplification ZS2 est découpée en un nombre de sous-zones égal au nombre de colonnes de mailles présents dans la zone ZS2. Le terme "colonne de mailles" est employé pour désigner un "empilement" de mailles verticales de la représentation maillée Ces sous-zones seront 5 alors des quadrilatères dont les longueurs des cotés correspondront à la discrétisation suivant les axes x et y de la représentation maillée (sous-zones ss_ZS2 en figure 3). On parlera par la suite de « sous-zone de simplification en colonne » pour ce type de sous-zones de simplification. Selon un mode préféré de mise en oeuvre de l'invention, une zone de simplification 10 ZS3 est découpée en sous-zones de simplification, en tenant compte des hétérogénéités des propriétés d'écoulement. On peut par exemple, de façon non limitative, délimiter des sous-zones de simplification d'une zone ZS3 en regroupant des sous-ensembles de mailles ayant des propriétés d'écoulement voisines. Par exemple dans le cas d'une zone ZS3 correspondant à un cylindre d'axe vertical et de section elliptique, un balayage 15 angulaire, centré sur le puits, est effectué. Les mailles contenues par degré échantillonné sont regroupées en une colonne de mailles. Ainsi, la zone ZS3 est caractérisée par au plus 360 colonnes de mailles. Pour chacune de ces colonnes de mailles, un calcul du tenseur de perméabilité de fracture équivalente est effectué (cf. étape suivante). Ce tenseur permet de caractériser les propriétés dynamiques du réseau de fractures de la 20 colonne de mailles étudiée. Les valeurs de perméabilité et l'orientation de l'écoulement obtenues sont alors comparées entre colonnes de mailles voisines. Si les propriétés équivalentes de deux colonnes de mailles voisines sont proches, ces deux colonnes de mailles seront considérées comme appartenant à une même sous-zone. Dans le cas contraire (différence de 20 % sur les valeurs de perméabilité principale ou de 10 degrés 25 sur les directions principales de perméabilité par exemple), les deux colonnes de mailles sont affectées à des sous-zones différentes. De cette manière, la limite extérieure, dans un plan horizontal, de la zone ZS3 considérée est divisée en arcs de forme elliptique (cf figure 3). Les limites dans un plan horizontal des sous-zones sont alors obtenues en reliant les points-extrémités de chacun de ces arcs au centre (puits W) de l'ellipse (lignes 30 en pointillés ss_ann_ZS3 et ss_ZS3 en figure 3) : chaque sous-zone à simplifier est ainsi définie comme la réunion des sous-ensembles de mailles constituant l'aire comprise entre le cloisonnement radial infra-zone (lignes en traits pointillés ss_ann_ZS3 et ss_ZS3) et les limites elliptiques interzones (lignes en trait plein ann_ZS3 et ZS3). On obtient de cette manière des sous-zones de simplification de la zone ZS3 sensiblement homogènes quant 35 aux propriétés d'écoulement, chacune étant caractérisée par une perméabilité équivalente 3041026 19 propre. On parlera par la suite de « sous-zone de simplification localement homogène » pour ce type de sous-zones de simplification. Selon un mode de mise en oeuvre de l'invention, un découpage en sous-zones de simplification localement homogènes peut être appliqué, en sus, à des sous-zones de 5 type annulaires pour une zone de simplification ZS3 (sous-zones ss_ann_ZS3 et ss_ZS3 en figure 3). Selon un mode de mise en oeuvre de l'invention, un découpage en sous-zones de simplification en colonne peut être appliqué, en sus, à des sous-zones de type annulaires pour une zone de simplification ZS2. 10 2.2. Détermination d'un tenseur de perméabilité équivalente Au cours de cette sous-étape, il s'agit de calculer un tenseur de perméabilité équivalente pour au moins chacune des zones de simplification déterminées à la sous-étape précédente, l'objectif étant, dans la sous-étape suivante, de remplacer le réseau de 15 fractures de chaque zone de simplification par un réseau simplifié possédant les mêmes propriétés d'écoulement que le réseau originel. Très préférentiellement, dans le cas où des sous-zones de simplification ont été définies pour une zone de simplification donnée à la sous-étape précédente, on calcule un tenseur de perméabilité équivalente en chacune de ces sous-zones de simplification. 20 Le calcul de perméabilités équivalentes est effectué à partir du modèle de réseau de fractures DFN, ou directement à partir des PSF. Ce type de calcul des perméabilités équivalentes est bien connu des spécialistes. On peut par exemple utiliser la méthode numérique de calcul de propriétés équivalentes de milieux fracturés, implémentée dans le 25 logiciel FracaFlow (IFP Energies nouvelles, France) et rappelé ci-après. Selon un mode de mise en oeuvre de l'invention, un tenseur de perméabilité, représentatif des propriétés d'écoulement du réseau de fractures discrétisé (DFN) peut être obtenu par exemple par l'une des deux méthodes de mise à l'échelle (méthodes dites d'upscaling) suivantes : 30 - une première méthode, analytique, dite "local analytical upscaling", est basée sur une approche analytique décrite par exemple dans le document (Oda, 1985) ou encore dans la demande de brevet FR 2918179 (US 8078405). Cette méthode présente l'avantage d'être très rapide. Son domaine d'application se limite, toutefois, à des réseaux 3041026 20 de fractures bien connectés. Dans le cas contraire, des erreurs importantes sur le tenseur de perméabilité peuvent être constatées. - une seconde méthode, numérique, dite "local numerical upscaling", est décrite dans le document (Bourbiaux, et al., 1998), ou dans la demande de brevet FR 2757947 5 (US 6023656) pour le calcul des perméabilités équivalentes. Cette méthode est basée sur la résolution numérique des équations d'écoulement sur un maillage discret du réseau de fractures pour différentes conditions aux limites appliquées à la colonne de mailles de calcul considéré. Le tenseur de perméabilités équivalentes est obtenu par identification des rapports entre débit et perte de charge aux limites du domaine de calcul. Cette 10 approche, plus coûteuse que la précédente, présente l'avantage de bien caractériser un réseau donné (même peu connecté). Selon un mode de réalisation de la présente invention, afin de profiter au mieux des avantages des deux méthodes de mise à l'échelle, le tenseur de perméabilité est déterminé par l'une ou l'autre de ces deux méthodes, par exemple, suivant la procédure 15 de sélection décrite dans le document FR 2918179 (US 8078405), fondée sur la valeur de l'indice de connectivité du réseau de fractures. Cet indice, représentatif du rapport entre le nombre d'intersections entre fractures et le nombre de fractures, est calculé pour chaque unité de la colonne de mailles considérée. Sa valeur permet de considérer le réseau comme très bien connecté, peu/mal connecté ou non connecté. Le choix de la méthode 20 de mise à l'échelle est alors effectué, comme décrit dans (Delorme, et al., 2008), et résumé ci-après : - Réseau bien connecté : pour ce cas caractérisé par un indice de connectivité voisin de 3 ou dépassant 3 (au moins 3 intersections par fracture du réseau en moyenne), la méthode d'upscaling analytique est choisie car sa précision est garantie compte tenu 25 de la bonne connectivité du réseau, avec l'avantage supplémentaire essentiel de la rapidité. - Réseau peu/mal connecté : dans ce cas, l'indice de connectivité est compris entre 1 et 3 (ce qui correspond à un nombre d'intersections de fractures compris entre une et trois fois le nombre de fractures) et la méthode d'upscaling numérique est utilisée pour 30 calculer de manière fiable le tenseur de perméabilité. - Lorsque le réseau est très peu voire non connecté (nombre d'intersections voisin ou inférieur au nombre de fractures), les fractures originelles (peu nombreuses) sont conservées, et il n'y a pas lieu de simplifier le réseau de fractures. 3041026 21 Selon un mode de mise en oeuvre de l'invention, on diagonalise le tenseur de perméabilité équivalente ainsi calculé afin de déterminer les trois directions principales de l'écoulement, que l'on notera par la suite clin, dir2 et dir3, la direction diri (respectivement dir2, et dir3) étant choisie de sorte à ce que son angle avec l'axe des abscisses 5 (respectivement l'axe des ordonnées, l'axe vertical) soit le plus faible. On notera alors kt k2, k3 les valeurs propres de ce tenseur, associées respectivement aux trois directions principales diri, dir2 et dir3. Selon un mode de mise en oeuvre de l'invention, on calcule en outre une porosité de fractures représentative du volume de fractures de la zone ou sous-zone. Cette porosité Of 10 est obtenue via une formule du type : of =Volume de fractures Volume total où le volume de fractures présentes dans le volume considéré peut être obtenu à partir du PSF ou du DFN déterminés à l'étape 1, et Volume total est le volume total de la zone ou sous-zone de simpification. 15 Selon un mode de mise en oeuvre de l'invention, on calcule en outre les dimensions caractéristiques d'un bloc matriciel représentatif (dit équivalent) pour chaque zone de simplification, et préférentiellement pour chaque sous-zone de simplification. On notera a, b et c les dimensions caractéristiques de ce bloc matriciel respectivement dans les directions clin, dir2 et dir3. Un tel calcul peut être réalisé selon la technique décrite dans 20 les demandes de brevet FR 2757957 (US 6064944) ou FR2923930 (US 8688424). 2.3. Détermination d'un paramètre P caractérisant la géométrie des fractures A partir d'au moins les paramètres statistiques (PSF) déterminés à l'étape 25 précédente, on détermine une valeur d'un paramètre, noté P; caractérisant à la fois l'orientation dans l'espace et la continuité vers les autres couches des fractures du réseau de fractures contenus au moins dans chacune des zones de simplification déterminées à la sous-étape 2.1. Ainsi, la valeur de ce paramètre P va permettre de déterminer si le réseau de fractures d'intérêt est de type 2.5D ou de type 3D dans chacune des zones de 30 simplification. Selon un mode préféré de mise en oeuvre de l'invention, un paramètre P est déterminé, le cas échéant, pour chaque sous-zone de simplification définie à l'étape 2.1, 3041026 22 en particulier dans le cas d'hétérogénéités des propriétés d'écoulement, et très préférentiellement dans le cas où des sous-zones de simplification dites localement homogènes ont été prédéfinies. 5 Selon un mode de mise en oeuvre de la présente invention, on définit un paramètre binaire P valant 0 ou 1 pour chaque (sous-)zone de simplification. La valeur 1 du paramètre P indique que les familles de fractures du réseau de fractures contenu dans une sous-zone de simplification donnée sont à la fois sensiblement orthogonales aux limites de la couche considérée et sensiblement continues vers les autres couches de la 10 formation (figure 2a). La valeur 0 indique que le réseau de fractures de la (sous-)zone de simplification considérée comprend au moins une famille de fractures de pendages quelconques (figure 2c) et/ou caractérisée par une continuité variable des fractures d'une couche à une autre (figure 2b). Autrement dit, la valeur 1 du paramètre P indique que le réseau de fractures de la (sous-)zone de simplification considérée est 2.5D, alors que la 15 valeur 0 du paramètre P indique que le réseau de fractures de la (sous-)zone de simplification considérée est 3D Selon un mode de mise en oeuvre de l'invention, on détermine un paramètre PC qualifiant la continuité de la fracturation d'une couche donnée vers les autres couches de la formation. Le calcul de ce paramètre peut être effectué par zone de simplification ou 20 par sous-zone de simplification. Pour ce faire, on considère un tenseur de perméabilité équivalente calculée sur la hauteur de la couche considérée pour la (sous-)zone de simplification considérée (cf sous-étape 2.2), et on diagonalise ce tenseur. Si k3 (perméabilité sub-verticale) est très différente de ki-Fk2 (somme des perméabilités subhorizontales), alors la fracturation n'est pas continue sur la hauteur de la couche 25 considérée pour la (sous-)zone considérée. Plus précisément, on peut définir un paramètre PC selon une formule du type : k3 PC = (kl+k2) Lorsque PC REVENDICATIONS1. Procédé pour l'exploitation d'un fluide au sein d'une formation souterraine comprenant au moins une couche traversée par un réseau de fractures et par au moins un puits, dans lequel, on construit une représentation maillée de ladite au moins une couche à partir de mesures de propriétés relatives à ladite au moins une couche, et on définit, au moins pour ladite couche, une zone d'intérêt regroupant un ensemble de mailles de ladite représentation maillée relative à ladite couche et comprenant au moins ledit puits, caractérisé en ce qu'on réalise au moins les étapes suivantes pour chacune desdites couches : A. à partir desdites mesures, on caractérise ledit réseau de fractures par des paramètres statistiques, et on construit, à partir desdits paramètres, un modèle de réseau de fractures discret; B. pour chacun desdits puits contenus dans ladite zone d'intérêt, on définit au moins trois zones de simplification (ZS1, ZS2, ZS3), une première zone de simplification (ZS1) contenant ledit puits, une seconde zone de simplification (ZS2) ayant comme frontière interne la frontière externe de ladite première zone (ZS1), et une troisième zone de simplification (ZS3) ayant comme frontière interne la frontière externe de ladite seconde zone (ZS2), l'ensemble desdites zones de simplification de l'ensemble desdits puits couvrant l'ensemble des mailles de ladite zone d'intérêt ; C. &partir d'au moins lesdits paramètres statistiques, on détermine, pour au moins chacune desdites zones de simplification (ZS1, ZS2, ZS3) de ladite zone d'intérêt, un tenseur de perméabilité équivalente et une valeur d'un paramètre P caractérisant l'orientation et la continuité vers les autres couches des fractures dudit réseau ; D. on détermine une simplification dudit modèle en chacune desdites zones de simplification, en fonction de ladite zone, dudit tenseur de perméabilité équivalente et de ladite valeur dudit paramètre P ; et en ce que, à partir de ladite représentation maillée, dudit modèle de réseau de fractures simplifié, et d'un simulateur d'écoulement, on définit un schéma d'exploitation optimal dudit fluide de ladite formation et on exploite ledit fluide de ladite formation en fonction dudit schéma d'exploitation optimal. 3041026 37 2- Procédé selon la revendication 1, dans lequel lesdits paramètres statistiques sont choisis parmi les paramètres suivants : densité de fractures, longueur de fractures, orientation des fractures dans l'espace, ouverture de fractures, et distribution de fractures au sein de ladite couche. 5 3- Procédé selon l'une des revendications précédentes, dans lequel lesdites zones de simplification sont des cylindres d'axe vertical et de section elliptique, centrées autour dudit puits. 4- Procédé selon la revendication 3, dans lequel le grand axe d'une desdites zones de simplification est orienté selon une direction préférentielle d'écoulement, ladite 10 direction préférentielle étant déterminée à partir d'une diagonalisation dudit tenseur de perméabilité équivalente. 5- Procédé selon l'une des revendications précédentes 3 ou 4, dans lequel la distance séparant deux desdites zones est déterminée en fonction d'une progression géométrique de raison constante. 15 6- Procédé selon l'une des revendications précédentes, dans lequel aucune simplification dudit modèle de réseau de fractures n'est réalisée dans ladite première zone de simplification (ZS1), une première simplification est réalisée dans ladite seconde zone (ZS2), et une seconde simplification est réalisée dans ladite troisième zone (ZS3), ladite seconde simplification étant plus importante que 20 ladite première simplification. 7- Procédé selon l'une des revendications précédentes, dans lequel lesdites seconde et/ou troisième zones de simplification (ZS2, ZS3) sont découpées en sous-zones de simplification annulaires. 8- Procédé selon l'une des revendications précédentes, dans lequel lesdites 25 secondes zones de simplification (ZS2) sont découpées en sous-zones de simplification en colonne. 9- Procédé selon l'une des revendications précédentes, dans lequel lesdites troisièmes zones de simplification (ZS3) sont découpées en sous-zones de simplification localement homogènes. 30 10- Procédé selon l'une des revendications précédentes, dans lequel on détermine ledit paramètre P pour au moins chacune desdites zones de simplification (ZS1, ZS2, ZS3), au moyen d' : 3041026 38 une valeur d'un paramètre PO mesurant l'orientation dans l'espace desdites fractures de ladite couche selon une formule du type : dev2 180 PO = - 2 -dey + 1 avec dey = arccos(chr3. z) * - et/ou 9cfr 0 (PI)' - une valeur d'un paramètre PC mesurant la continuité desdites fractures de 5 ladite couche vers les autres couches selon une formule du type : k3 PC = (k1+k2)' où kl, k2, et k3 sont les valeurs propres dudit tenseur de perméabilité après diagonalisation. 11- Procédé selon l'une des revendications précédentes, dans lequel ladite 10 simplification consiste en le calcul d'au moins une valeur d'espacement entre lesdites fractures, une ouverture desdites fractures, et une valeur de perméabilité. 12- Procédé selon l'une des revendications précédentes, dans lequel ladite simplification dudit modèle est réalisée selon trois directions de l'espace si le paramètre P caractérise des fractures non sensiblement orthogonales aux 15 couches et/ou continues vers les autres couches. 13- Procédé selon l'une des revendications 1 à 11, dans lequel ladite simplification dudit modèle est réalisée selon deux directions de l'espace si le paramètre P caractérise des fractures sensiblement orthogonales aux couches et continues vers les autres couches. 20 14- Procédé selon l'une des revendications précédentes, dans lequel ledit schéma d'exploitation optimal est défini en déterminant un procédé de récupération dudit fluide, ainsi qu'un nombre, une implantation et une géométrie de puits injecteurs et producteurs permettant de satisfaire des critères technico-économiques prédéfinis. 15- Procédé selon la revendication 14, dans lequel ladite exploitation dudit fluide en 25 fonction dudit schéma d'exploitation optimal consiste au moins à forer lesdits puits injecteurs et producteurs et à produire ledit fluide selon ledit procédé de récupération. 16- Produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou 30 exécutable par un processeur, comprenant des instructions de code de programme pour la mise en oeuvre du procédé selon l'une des revendications précédentes, lorsque ledit programme est exécuté sur un ordinateur.