FR 3046810 A1 20170721 FR 1650309 A 20160115 1 La présente invention concerne le domaine de l'exploitation de milieux souterrains traversés par des puits forés, et en particulier l'exploitation de gisements souterrains contenant 5 des hydrocarbures. Pour exploiter au mieux les milieux souterrains, les industries pétrolières définissent des modèles numériques des zones réservoirs afin de piloter le processus d'exploitation des champs pétroliers. Le procédé de production selon l'invention met en oeuvre un modèle 10 optimisé du point de vue de la détermination des indices, ou index, de productivité. La technique qui consiste à prendre en compte les effets thermiques au niveau de la perméabilité prise en compte dans les modèles de réservoir, en mettant à jour ladite perméabilité en fonction de la pression, température, variation de la porosité et de contraintes 15 volumiques est connue. D'ailleurs, une telle technique a été proposée par Ojagbohunmi et al. (2012) où la perméabilité est fonction de la porosité, qui est elle-même fonction de la pression et de la température (Tran et al., 2004). D'autre part, il existe des méthodes de calcul des indices de productivité des puits (IP). En particulier, ces indices de productivité peuvent être estimés en fonction de la perméabilité des 20 mailles sous-jacentes par les méthodes de Peaceman (1983) ou Mochizuki (1995) pour ne citer que les plus connues. S. Ojagbohunmi, R. Chalaturnyk, J. Leung, Coupling of stress dependent relative permeability and reservoir simulation. SPE 154083, 2012. D. Tran, A. Settari, L. Nghiem, New iterative coupling between a reservoir simulator 25 and a geomechanics module. SPE 2004. D.W. Peaceman, Interpretation of well-block pressures in numerical reservoir simulation with nonsquare grid blocks and anisotropic permeability. SPE, 1983. S. Mochizuki, Well productivity for arbitrarily inclined well. SPE 029133, 1995. 3046810 2 Dans la littérature, on trouve des méthodes qui combinent les deux techniques et prennent en compte les effets de la pression sur les indices de productivité, en les faisant dépendre de la perméabilité, dépendante de la porosité, elle-même fonction de la pression. Toutefois, on ne trouve pas dans la littérature de méthode où les index de productivité des puits sont fonction de la perméabilité qui serait elle-même fonction de la pression, et aussi plus particulièrement de la température. En effet, la porosité n'est généralement pas dépendante de la température. Cela sous-entend que dans tous les cas, si les effets thermiques peuvent être pris en compte dans la perméabilité au niveau du champ, ils ne sont pas pris en compte spécifiquement au niveau des puits forés, ce qui induit une importante incertitude sur la production d'hydrocarbures des puits par l'intermédiaire du calcul de l'index de production IP. La présente invention propose de remédier à ce problème. Cette invention prend mieux en compte les effets thermiques, et plus particulièrement les effets thermiques dans l'environnement des puits. En intégrant ces effets thermiques pour la détermination d'indices de productivité des puits, on peut obtenir un modèle de réservoir plus précis et aboutir à des résultats, et des estimations de production plus proches de la réalité. Ainsi la présente invention concerne un procédé de production d'hydrocarbures contenus dans une formation souterraine traversée par au moins un puits. Pour ce procédé, on réalise les étapes suivantes : - on construit à l'aide d'un ordinateur un modèle numérique de ladite formation souterraine comportant un maillage, et prenant en compte au moins des données pétrophysiques appliquées aux mailles, et les caractéristiques d'au moins un puits traversant ladite formation souterraine, - on détermine une fonction de dépendances en pression et température des perméabilités prises en compte dans ledit modèle, - on détermine un modèle d'indice de productivité fonction de la perméabilité desdites mailles, - on détermine pour au moins une maille dudit modèle un indice de productivité en dépendance de la pression et de la température, en prenant en compte lesdites 3046810 3 fonctions de dépendances en pression et température des perméabilités, ainsi que ledit modèle d'indice de productivité, et - on exploite ladite formation souterraine en fonction dudit indice de productivité. 5 Selon un mode de réalisation de l'invention, lesdites fonctions de dépendances en pression et température des perméabilités prises en compte dans ledit modèle comportent des multiplicateurs de perméabilité Mult,(T , P). Conformément à un mode de réalisation de l'invention, lesdits multiplicateurs sont 10 déterminés en utilisant les modifications de la porosité en fonction de la température et la pression. Selon une mise en oeuvre de l'invention, le calcul de l'indice de productivité est effectué selon : \ 21TV Kx KY h IP(1(eff eff eff ln(;)+ s 15 avec l'épaisseur perpendiculaire à l'écoulement h, rayon de puits rw, rayon de drainage ro et skin factor s. Selon un mode de réalisation de l'invention, le procédé est appliqué à un procédé de drainage thermique. Avantageusement, on peut déterminer l'indice de productivité pour les mailles dudit 20 modèle numérique impactés par ledit puits. En outre, l'invention concerne un produit programme d'ordinateur téléchargeable 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 25 en oeuvre du procédé selon l'une des revendications précédentes, lorsque ledit programme est exécuté sur un ordinateur. L'invention sera mieux comprise à la lecture qui suit de la description et d'un exemple, 30 nullement limitatif, et illustré par les figures ci-après annexées, notamment : 3046810 4 - la figure 1 montre les courbes de production cumulée en fonction du temps dans les trois cas suivant : - sans les effets géomécaniques, - avec effet de la pression seul sur les perméabilités et sur l'IP, - avec les effets thermiques et de pression sur les perméabilités et sur 5 l'IP ; - les figures 2a et 2b donnent les débits de production pour les mêmes trois cas et leurs cumuls de production. La présente invention concerne le domaine de l'exploitation de milieux souterrains. 10 Cette invention permet de mieux prendre en compte les effets thermiques et plus particlulièrement les effets thermiques proches puits en intégrant ces effets aux index de productivité IP des puits. Cela permet donc d'obtenir un modèle de réservoir plus précis. En particulier, cette invention permet d'obtenir des résultats de production plus précis. On rappelle qu'un indice (ou index) de productivité IP définit la capacité de production 15 d'un puits. Cet indice de productivité peut permettre de prévoir l'équipement d'un puits pour un débit déterminé ou la perte de charge qu'il faut imposer à un puits. Afin d'exploiter au mieux les milieux souterrains, les industries pétrolières définissent des modèles numériques permettant de mieux piloter et exploiter les champs pétroliers. 20 La première étape de ces études d'exploitation consiste à créer un modèle numérique qui représente au mieux les données géologiques et pétrophysiques du milieu poreux, par exemple les facies, les perméabilités, ou les porosités de la roche réservoir. On définit aussi dans ce modèle, la position des puits forés ainsi que leur mode de fonctionnement : débit, pression, fluides injectés, nature des produits injectés,... 25 La deuxième étape consiste à définir les effets thermiques. Classiquement, on ajoute aux perméabilités des dépendances en température et pression : K (T , P). Afin de ne pas devoir changer les perméabilités définies à l'étape première, on prend en compte les effets thermiques en définissant un multiplicateur de perméabilité Mult(T, P) afin d'avoir une perméabilité effective Keff (T, P) = Mult(T, P) * K(T, P). De cette manière les effets thermiques sont 30 pris en compte sur toutes les mailles du modèle, mais ne sont pas pris en compte sur les index 3046810 5 de productivité des puits. Les facteurs Mult(T, P) peuvent être obtenus par des lois physiques ou empiriques, par exemple comme selon Touhidi-B aghini (1998). -Touhidi-Baghini, A.: "Absolute Permeability of McMurray Formation Od Sands at 5 Low Confining Stress, Departement of Civil and Environmental Engineering, University of Alberta, Ph.D. Dissertation, 1998 La présente invention permet de prendre en compte ces effets sur les indices de productivité (IP) en choisissant tout d'abord une méthode de calcul d'IP fonction de la 10 perméabilité, par exemple, la méthode connue de « Peaceman », ou la méthode de projection (basée sur la méthode de « Peaceman »), ou autre. En effet, physiquement, la perméabilité de la maille perforée par un puits a obligatoirement un impact sur l'indice de productivité du puits. On calcule alors les indices de productivité des puits IP(T, P) pour au moins une 15 maille, et de préférence pour chaque maille perforée par les puits en utilisant la perméabilité mise à jour dans la deuxième étape : Keff (T, P), c'est-à-dire en tenant compte du multiplicateur de perméabilité, ou de la dépendance en température et pression de la perméabilité. La méthode de « Peaceman » citée plus haut peut notamment être retrouvée décrite 20 dans les documents suivants : - Peaceman, D.W.: "Interpretation of Well-Block Pressures in Numerical Reservoir Simulation", paper SPE 6893, presented at the SPE-AIME 52nd Annual Fall Technical Conference and Exhibition, Denver, Oct 9-12, 1977 - Peaceman, D.W.: "Interpretation of Well-Block Pressures in Numerical Reservoir 25 Simulation With Nonsquare Grid Blocks and Anisotropic Permeability", paper SPE 10528, presented at the 1982 SPE Symposium on Reservoir Simulation, New Orleans, Jan 31-Feb 3 - Peaceman, D.W.: "Interpretation of Wellblock Pressures in Numerical Reservoir Simulation: Part3 - Off-Center and Multiple Wells Within a Wellblock", paper SPE 16976, presented at the 1987 SPE Annual Technical Conference and Exhibition, Dallas, Sept 27-30 30 3046810 6 Selon la présente invention, le procédé de production d'hydrocarbures contenus dans une formation souterraine, par exemple une roche réservoir, traversée par au moins un puits foré consiste en une méthode comportant les cinq étapes suivantes : 5 1 Génération d'un modèle numérique du réservoir ; 2 Définition des dépendances en température de la perméabilité ; 3 Choix d'une modèle des indices de productivité des puits avec dépendance à la perméabilité de la maille perforée ; 4 Calcul des indices de productivité des puits en utilisant les perméabilités 10 dépendantes de la température obtenues à l'étape 2 ; 5 Exploitation du réservoir 1) Génération d'un modèle géologique fin. Cette première étape consiste à construire à l'aide d'un ordinateur et d'un logiciel de 15 modélisation, un modèle de réservoir comprenant un maillage, des propriétés pétrophysiques appliquées aux mailles, ainsi que des positions et trajectoires des puits. Ainsi, le modèle de réservoir représente la formation souterraine. Les propriétés physiques peuvent comprendre notamment la perméabilité et porosité. 20 2) Définition des dépendances en pression et température de la perméabilité. La deuxième étape consiste à définir les effets thermiques. Classiquement, on ajoute aux perméabilités des dépendances en température et pression : K,(T , P) ou i représente la direction (x, y ou z). La perméabilité K,(T , P) est obtenue dans le modèle de réservoir construit à l'étape précédente. 25 Selon un mode de réalisation de l'invention, afin de modifier les perméabilités initiales définies à l'étape 1 dans le modèle géologique, on peut prendre en compte les effets thermiques en introduisant un multiplicateur de perméabilité Mult,(T , P) afin d'avoir une perméabilité effective IÇ 11(T , P) = Mult,(T , P) * Conformément à une variante de réalisation de l'invention, les facteurs Mult(T, P) peuvent 30 être obtenus par des lois physiques ou empiriques (par exemple Touhidi-Baghini (1998)). 3046810 7 Selon une variante de réalisation de l'invention, les facteurs Mult(T, P) peuvent être obtenus au moyen d'un simulateur de réservoir, par exemple le logiciel PumaFlow (IFP Energies nouvelles, France). 5 3) Choix d'un modèle des indices de productivité des puits avec dépendance à la perméabilité de la maille perforée. Lors de cette étape on détermine une relation (appelé également modèle ou méthode de calcul) entre les indices de productivité et la perméabilité. On fait le choix d'un modèle faisant 10 intervenir la perméabilité dite effective. Selon un mode de réalisation de l'invention, l'indice de productivité pour un puits vertical peut être estimé par la formule de Peaceman - 2' c.\ IK fK eYff h InM+ s Avec : l'épaisseur perpendiculaire à l'écoulement h, rayon de puits rw, rayon de drainage ro 15 et « skin factor » (qui est un facteur sans dimension qui traduit les pertes de charge engendrées par une restriction de débit à proximité des puits) s. On remarque que ro est une fonction de la perméabilité : IP(Kff = (1) ro = 0.28 - (r,, /4 + ( )1/4 K:ff (K:fr )1/ 2A,x 2 (K4f )1/ 2 Ay 2 + içff 1/ 2 (2) 20 4) Calcul des indices de productivité des puits en utilisant les perméabilités dépendantes de T et P obtenues à l'étape 2. On calcule alors au moins un indice de productivité des puits IP(T, P) pour au moins une maille du modèle de réservoir, et de préférence pour chaque maille du modèle de réservoir perforée par le(s) puits en utilisant les perméabilités effectives mises à jour dans l'étape 2 : 25 Kjff (T, P) et elfe, P), c'est-à-dire en tenant compte du multiplicateur de perméabilité, ou plus généralement de la fonction de dépendance en pression et température de la perméabilité. Si f(x,....) représente la formule utilisée pour calculer l'IP, alors : 3046810 8 IP (T , P) = f (IQ f(T , PI K: fi(T , P) ....) En utilisant les formules de « Peaceman » (1, 2) et l'expression de Keff (T, P), on obtient l'indice de productivité en fonction de T et P IP(T' P) = 27 1-h i \ V Mult x(T , P)* Kexff * Mult y(T , P)* KeYff ln Iir°I+ s 5 avec ro*: . I Mult Y (T , P)* KY Ax2 + Mult (T P)* Kx Ay2 ro* =0.28 eff VMUity (T, P)* KeYff ± VMUitx (T, P)* Kexff La présente invention est basée sur la prise en compte des effets thermiques sur les index de productivité des puits en tenant compte des effets thermiques sur la perméabilité par 10 l'exhibition du multiplicateur de perméabilité dans la formule d'IP. 5) Exploitation du réservoir Lors de cette étape, on exploite la formation souterraine, c'est-à-dire on produit des hydrocarbures, en fonction des indices de productivité de puits déterminés à l'étape 15 précédente. Selon un mode de réalisation de l'invention, à partir des indices de productivité des puits, les spécialistes peuvent déterminer plusieurs schémas d'exploitation correspondant à différentes configurations possibles d'exploitation du réservoir souterrain : emplacement des puits producteurs et/ou injecteurs, valeurs cibles pour les débits par puits et/ou pour le 20 réservoir, le type d'outils utilisés, les fluides utilisés, injectés et/ou récupérés... Ces schémas de production peuvent être simulés au moyen d'un simulateur d'écoulement et du modèle de réservoir. A partir des indices de productivité des puits, et des éventuelles simulations d'écoulement, les spécialistes peuvent par comparaison choisir le schéma d'exploitation qui 25 leur semble le plus pertinent. Par exemple : 3046810 9 - en comparant le maximum du volume d'huile récupéré, on peut déterminer le schéma de production susceptible de fournir le maximum de récupération ou d'être le plus rentable. 5 - en comparant l'écart type du volume d'huile récupéré, on peut déterminer le schéma de production le moins risqué. On exploite alors le réservoir selon le schéma d'exploitation défini par exemple en forant de nouveaux puits (producteur ou injecteur), en modifiant les outils utilisés, en modifiant les débits et/ou la nature de fluides injectés... 10 L'invention concerne par ailleurs un produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur. Ce programme comprend des instructions de code de programme pour la mise en oeuvre du procédé tel que décrit ci-dessus, lorsque le programme 15 est exécuté sur un ordinateur. Exemples Afin de montrer l'intérêt de cette méthode pour la simulation de production, on 20 l'applique sur un cas simplifié de production de réservoir. On utilise un modèle de simulation de SAGD (Steam-assisted gravity drainage) avec un puits producteur et un puits injecteur. D'abord, on préchauffe les puits fermés, ensuite on ouvre les puits pour équilibrer la pression. On injecte la vapeur dans le puits injecteur et on observe la production d'huile dans le puits producteur. Les variations de la température et de la pression 25 sont très importantes dans ce processus ce qui permettra de mettre en évidence leurs rôles dans les effets géomécaniques. On va démontrer l'influence de la prise en compte de la variation de perméabilité avec la température et la pression sur la production. Pour cela on compare des résultats des simulations : 30 - sans les effets géomécaniques, selon l'art antérieur (noté AA), 3046810 10 - seulement avec les effets de la pression sur les perméabilités et sur l'IP, selon l'art antérieur (noté Keff (P) et IP (P)) ; et - avec les effets thermiques et de pression sur les perméabilités et sur l'IP selon 5 l'invention (noté Keff(T,P) et IP (T, P)). Les perméabilités initiales sont Ky, = Ky = 2000 mD et Kz=1000 mD. Les multiplicateurs sont calculés en utilisant les modifications de la porosité en fonction de la température et la pression 0 = 00exPtcp(P Po)+ -To)} 10 et en utilisant une approximation de la loi empirique de Touhidi-Baghini (1998) où 00, 7'0 et Po sont respectivement : la porosité, la température et la pression initiales, cp compressibilité de roche, cg, expansion thermique K eff MUlti (P, = = exp Ki 00 15 avec i désignant l'un des axes x, y et z, et c, = cy = 2,c, = 5 , et avec c.' , cy et cz paramètres du matériau et ev une contrainte volumique calculée en considérant les variations de la porosité ainsi que de la température, ev = 00c r(P - po+ 3ad (T -To) avec ad coefficient d'expansion thermique drainée. 20 Leurs valeurs sont données dans les tableaux 1 et 2, les multiplicateurs sont dépendants de la direction considérée horizontale ou verticale. Pour cet exemple illustratif et non limitatif, les facteurs Mult(T, P) sont obtenus par des lois physiques ou empiriques, telles que décrites dans le document Touhidi-B aghini (1998). - Touhidi-Baghini, A.: "Absolute Permeability of McMurray Formation Oil Sands at 25 Low Confining Stress, Departement of Civil and Environmental Engineering, University of Alberta, Ph.D. Dissertation, 1998 La figure 1 montre les courbes de production cumulée d'huile (COS) en fonction du temps (T) pour tous ces cas. 3046810 11 En considérant une perméabilité et des IP variables en fonction de la température et de la pression, Keff (T, P) et IP(T, P), par la méthode selon l'invention définie ci-dessus (courbe (-^-)), on obtient l'estimation de la production significativement différente par rapport au cas 5 sans les effets thermiques Keff (P) et IP(P) (courbe (--^--)) selon l'art antérieur. Dans un premier temps, l'effet de la hausse de la température augmente les perméabilités et les IP. Ce qui amène à produire plus rapidement que dans la cas sans effet de la température. L'huile étant produite plus rapidement, on assiste à une baisse de production après 700 jours. La figure 2a montre le débit d'huile en surface en fonction du temps. On observe que la 10 prise en compte des effets géomécaniques pour les indices de productivité influence le débit d'huile dans le puits producteur. En conséquence, la production d'huile cumulée (COS) à la fin de période est modifié (figure 2b). En effet, la température a tendance à augmenter la perméabilité du milieu, par exemple, par l'activation de micro-fractures et par voie de conséquence, elle augmente aussi la valeur de l'index de productivité. La prise en comptes des 15 IP thermique (courbe Keff (T, P) et IP(T, P)) selon l'invention accélère l'arrivée d'eau au producteur, et donc une chute de la productivité d'huile (Figure 2a) aux alentours de 750 jours. Sur la figure 2b, on remarque que la chute de productivité d'huile arrive plus tôt si l'on tient compte des effets thermiques sur les IPs. Ainsi, les IP (T , P) obtenus selon cette présente invention permettent de mieux estimer la production et d'optimiser un procédé de production. 20 En conséquence, il est important de tenir compte des variations d' IP provoquées par des variations de température pour bien estimer la production, notamment dans les cas où la température change significativement. 25 Tableau 1 : Multiplicateurs Mu/tx,y(T, P) définies pour la pression P (col) et la température T (lignes) P (bar) \ 0 5 10 15 20 25 30 35 40 45 50 55 60 T (°C) 0 0.9418 0.9512 0.9608 0.9704 0.9802 0.9900 1.0000 1.0101 1.0202 1.0305 1.0408 1.0513 1.0618 25 0.9656 0.9753 0.9851 0.9950 1.0050 1.0151 1.0253 1.0356 1.0460 1.0565 1.0672 1.0779 1.0887 50 0.9900 1.0000 1.0101 1.0202 1.0305 1.0408 1.0513 1.0618 1.0725 1.0833 1.0942 1.1052 1.1163 3046810 12 75 1.0151 1.0253 1.0356 1.0460 1.0565 1.0672 1.0779 1.0887 1.0997 1.1107 1.1219 1.1331 1.1445 100 1.0408 1.0513 1.0618 1.0725 1.0833 1.0942 1.1052 1.1163 1.1275 1.1388 1.1503 1.1618 1.1735 125 1.0672 1.0779 1.0887 1.0997 1.1107 1.1219 1.1331 1.1445 1.1560 1.1677 1.1794 1.1912 1.2032 150 1.0942 1.1052 1.1163 1.1275 1.1388 1.1503 1.1618 1.1735 1.1853 1.1972 1.2092 1.2214 1.2337 175 1.1219 1.1331 1.1445 1.1560 1.1677 1.1794 1.1912 1.2032 1.2153 1.2275 1.2399 1.2523 1.2649 200 1.1503 1.1618 1.1735 1.1853 1.1972 1.2092 1.2214 1.2337 1.2461 1.2586 1.2712 1.2840 1.2969 225 1.1794 1.1912 1.2032 1.2153 1.2275 1.2399 1.2523 1.2649 1.2776 1.2905 1.3034 1.3165 1.3298 250 1.2092 1.2214 1.2337 1.2461 1.2586 1.2712 1.2840 1.2969 1.3100 1.3231 1.3364 1.3499 1.3634 275 1.2399 1.2523 1.2649 1.2776 1.2905 1.3034 1.3165 1.3298 1.3431 1.3566 1.3703 1.3840 1.3979 300 1.2712 1.2840 1.2969 1.3100 1.3231 1.3364 1.3499 1.3634 1.3771 1.3910 1.4049 1.4191 1.4333 Tableau 2 : Multiplicateurs Mu/tz(T, P) définies pour la pression P (col) et la température T P (bar) \ T (°C) 0 5 10 15 20 25 30 35 40 45 50 55 60 0 0.8607 0.8825 0.9048 0.9277 0.9512 0.9753 1.0000 1.0253 1.0513 1.0779 1.1052 1.1331 1.1618 25 0.9162 0.9394 0.9632 0.9876 1.0126 1.0382 1.0645 1.0914 1.1191 1.1474 1.1764 1.2062 1.2368 50 0.9753 1.0000 1.0253 1.0513 1.0779 1.1052 1.1331 1.1618 1.1912 1.2214 1.2523 1.2840 1.3165 75 1.0382 1.0645 1.0914 1.1191 1.1474 1.1764 1.2062 1.2368 1.2681 1.3002 1.3331 1.3668 1.4014 100 1.1052 1.1331 1.1618 1.1912 1.2214 1.2523 1.2840 1.3165 1.3499 1.3840 1.4191 1.4550 1.4918 125 1.1764 1.2062 1.2368 1.2681 1.3001 1.3331 1.3668 1.4014 1.4369 1.4733 1.5106 1.5488 1.5880 150 1.2523 1.2840 1.3165 1.3499 1.3840 1.4191 1.4550 1.4918 1.5296 1.5683 1.6080 1.6487 1.6905 175 1.3331 1.3668 1.4014 1.4369 1.4733 1.5106 1.5488 1.5880 1.6282 1.6695 1.7117 1.7551 1.7995 200 1.4191 1.4550 1.4918 1.5296 1.5683 1.6080 1.6487 1.6905 1.7333 1.7771 1.8221 1.8682 1.9155 225 1.5106 1.5488 1.5880 1.6286 1.6695 1.7117 1.7551 1.7995 1.8450 1.8917 1.9396 1.9887 2.0391 250 1.6080 1.6487 1.6905 1.7333 1.7771 1.8221 1.8682 1.9155 1.9640 2.0138 2.0647 2.11700 2.1706 275 1.7117 1.7551 1.7995 1.8450 1.8917 1.9396 1.9887 2.0391 2.0907 2.1436 2.1979 2.2535 2.3106 300 1.8221 1.8682 1.9155 1.9640 2.0138 2.0647 2.1170 2.1706 2.2255 2.2819 2.3396 2.3989 2.4596 REVENDICATIONS1) Procédé de production d'hydrocarbures contenus dans une formation souterraine traversée par au moins un puits, caractérisé en ce qu'on réalise les étapes suivantes : - on construit à l'aide d'un ordinateur un modèle numérique de ladite formation souterraine comportant un maillage, et prenant en compte au moins des données pétrophysiques appliquées aux mailles, et les caractéristiques d'au moins un puits traversant ladite formation souterraine, - on détermine une fonction de dépendances en pression et température des perméabilités prises en compte dans ledit modèle, - on détermine un modèle d'indice de productivité fonction de la perméabilité desdites mailles, - on détermine pour au moins une maille dudit modèle un indice de productivité en dépendance de la pression et de la température, en prenant en compte lesdites fonctions de dépendances en pression et température des perméabilités, ainsi que ledit modèle d'indice de productivité, et - on exploite ladite formation souterraine en fonction dudit indice de productivité. 2) Procédé de production d'hydrocarbures selon la revendication 1, dans lequel lesdites fonctions de dépendances en pression et température des perméabilités prises en compte dans ledit modèle comportent des multiplicateurs de perméabilité Mult,(T , P). 3) Procédé de production d'hydrocarbures selon la revendication 2, dans lequel lesdits multiplicateurs sont déterminés en utilisant les modifications de la porosité en fonction de la température et la pression. 4) Procédé de production d'hydrocarbures selon l'une des revendications 1 à 3, dans lequel le calcul de l'indice de productivité est effectué selon : 3046810 14 \ 21-cifrx KY h 113(Keff) = eff eff ln(1r)+s avec l'épaisseur perpendiculaire à l'écoulement h, rayon de puits rw, rayon de drainage ro et skin factor s. 5 5) Procédé de production d'hydrocarbures selon l'une des revendications 1 à 4, appliqué à un procédé de drainage thermique. 6) Procédé de production d'hydrocarbures selon l'une des revendications précédentes, 10 dans lequel on détermine l'indice de productivité pour les mailles dudit modèle numérique impactés par ledit puits. 7) Produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un 15 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.