FR 2824651 A1 20021115 FR 0106216 A 20010509 en mode d'adressage actif. DESIGNATION DU DOMAINE TECHNIQUE La présente invention concerne une méthode pour déterminer par simulation numérique les conditions optimales à imposer dans un puits horizontal (ou complexe) foré s au travers d'un gisement souterrain, pour éliminer progressivement (restaurer) par balayage naturel au moyen des fluides de production issus du gisement, des dépôts ou cakes qui se sont formés dans au moins une zone à la périphérie du puits, suite aux opérations de forage et de complétion. Les gens de l'art savent que l'on distingue les cakes dits internes formés par 0 envahissement des pores de la roche par la boue et les cakes dits externes constitués d'un revêtement de boue sur la paroi externe du puits. ETAT DE LA TECHNIQUE Les dommages causés aux formations environnant les puits horizontaux (ou complexes), souvent des trous ouverts équipés pour la production, constitue un point critique pour les champs pétroliers en offshore profond o, pour obtenir des coûts de développement acceptables, on ne met en production qu'un nombre limité de puits très productifs. Les essais que l' on peut faire pour caractériser l'endommagement des formations au voisinage d'un puits sont d'une importance primordiale. Ils permettent de choisir le fluide de forage le mieux approprié pour minimiser ou réduire la détérioration de la perméabilité au voisinage des puits et aussi d'optimiser les techniques de nettoyage des puits. Au cours des cinq dernières années, le demandeur a mis au point un équipement spécifique d'essai en laboratoire et des procédures destinées à caractériser l' endommagement des formations du au forage au cours d' opérations en surpression et à 2s quantifier les performances de différentes techniques de nettoyage mises en _uvre dans l'industrie comme en témoignent les publications suivantes: - Alfenore, J. et al, " What Really Matters in our Quest of Minimizing Formation Damage in Open Hole Horizontal Wells ", 1999, SPE 54731; - Longeron, D. et al, " Experimental Approach to Characterize Drilling Mud Invasion, Formation Damage and Cleanup Efficiency in Horizontal Wells with Openhole Completions " 2000, SPE 58737; ou - Longeron, D. et al, " an Integrated Experimental Approach for Evaluating s Formation Damages due to Drilling and Completion Fluids ", l995, SPE 30089. Cependant les études menées en laboratoire s'avèrent souvent insuffisantes à elles seules, pour modéliser de façon réaliste, les conditions de production à appliquer dans les puits pour restaurer au mieux la perméabilité des formations environnantes, sans pour autant provoquer des venues de sable. La modélisation des processus de restauration des 0 formations environnant un puits présentent un grand intérêt économique pour la mise en exploitation des champs. LA METHODE SELON L'INVENTION La méthode selon l' invention permet de simuler au mieux les conditions optimales à imposer dans un puits foré au travers d'un gisement souterrain à une trajectoire lS quelconque, pour l'élimination progressive par des fluides issus du gisement, des dépôts ou cakes qui se sont formés dans au moins une zone à la périphérie du puits, suite aux opérations de forage. Elle comporte l' acquisition de données initiales obtenues par des mesures en laboratoire des valeurs de perméabilité initiale (ki) des formations entourant le puits, de l'épaisseur des cakes ainsi que des valeurs de perméabilité endommagée (kd) et de perméabilité restaurée (kf) de la zone, en fonction de la distance (r) à la paroi du puits, la discrétisation de la zone endommagée par un maillage cylindrique en 3D formant des blocs d'épaisseur radiale faible relativement au diamètre du puits et la résolution dans ce maillage, d'équations de diffusivité modélisant les écoulements des fluides au travers des 2s cakes en tenant compte des données initiales mesurées et en modélisant l'évolution de la perméabilité en fonction des débits (Q) de fluides s'écoulant au travers des cakes, et on en déduit les conditions optimales à appliquer pour la mise en production du puits. On modélise la restauration de la perméabilité en tout point à distance (r) de la paroi en considérant par exemple que la perméabilité varie proportionnellement à l'écart entre la perméabilité endommagée (kd) et la perméabilité restaurée (kf), le coefficient de proportionnalité dépendant d'une loi empirique de variation de la perméabilité en fonction de la quantité de fluides au travers des cakes. La simulation opérée suivant la méthode permet aux ingénieurs de réservoir de mieux prévoir le meilleur schéma d' exploitation du gisement en évitant les inconvénients tels que les venues de sable. Elle permet aussi aux foreurs, compte tenu de données de perméabilité connues ou estimées, de choisir des fluides convenant plus particulièrement pour le forage des puits et la mise en place des équipements. PRESENTATION SOMMAIRE DES FIGURES 0 D'autres caractéristiques et avantages de la méthode et du dispositif selon l'invention, appara^tront à la lecture de la description ci-après d'un exemple non limitatif de réalisation, en se référant aux dessins annexés o: - la Fig.1 montre les courbes de variation en fonction de la distance r à la paroi du puits endommagé, d'un premier coefficient multiplicateur c(r) de la perméabilité d' endommagement et d' un deuxième coefficient multiplicateur c2(r) de la perméabilité restaurée; - la Fig.2 montre une loi empirique de variation d'un coefficient de variation de la perméabilité à la distance r de la paroi du puits endommagé, en fonction du débit de fluides Qs au travers des cakes; - la Fig.3 montre un exemple de maillage radial pour la résolution des équations de diffusivité; - la Fig.4 illustre le calcul du flux F avec un maillage radial; - les Fig.Sa et Sb illustrent le calcul de l'IP numérique respectivement sans cake externe et avec cake externe Cext, autravers d'une maille Wcell; - la Fig.6 montre schématiquement une portion de puits de longueur L et de rayon r\v comportant 4 zones de profondeur r centrées autour du puits, avec des perméabilités k différentes l00mD ou l000mD et un cake interne d'épaissseur rin; - les Fig.7, 8 montrent les variations en fonction de la distance d au puits, des coefficients multiplicateurs respectivement de perméabilité endommagée C1(r) et de perméabilité restaurée ou de retour c2(r), qui ont été mesurées au laboratoire dans différentes zones et utilisées dans les exemples; - la Fig.9 montre la courbe de variation de la perméabilité c(r) dans le cake interne en fonction du volume cumulé q de fluide par unité de surface offerte à l'écoulement, mesuré au laboratoire et utilisée dans les exemples; - les Fig. lOa à lOc montrent respectivement les variations en fonction du temps t(d) exprimé en jours, des déLits d'huile FR (en m3/j) dans différentes zones perforces le long o du puits, correspondant à 3 simulations différentes SM1 à SM3, dans l'exemple 1 (cas a); - les Fig. lla et llb montrent les variations en fonction du temps t(d) exprimé en jours, du coefficient de perméabilité c(r) du cake interne dans deux zones différentes le long du puits (exemple 1); - la Fig. 12 montre la variation en fonction du temps du débit total FR(m3/d) dans le cas c de l'exemple 1, pour trois simulations différentes SM1 à SM3; - la Fig. 13 montre la distribution du cake externe le long de la portion de puits, dans l'exemple 2; - les Fig. 14a à 14c montrent respectivement, dans l'exemple 2, la distribution sur la longueur L(m) du puits, du cake externe (Fig.14a) et du débit FR le long du puits au temps t=0.5 j (Fig.14b) et au temps t=5j (Fig.14c); - les Fig. 15a à 15f montrent respectivement, dans l'exemple 2, la distribution sur la longueur L(m) du puits, du cake externe (Fig.lSa) et du débit FR le long du puits, respectivement au temps t = 0.1j (Fig. 15b), t = 0.3j (Fig.15c), t = 0.5j (Fig.15d), t = lj (Fig.15e), et t = 5j (Fig.15f); - la Fig. 16 montre le débit total du puits FR en fonction du temps exprimé en jours, dans l'exemple 2, pour les deux cas cl et c2;; la Fig.17 est un tableau montrant un exemple de maillage avec NX mailles réparties le long du puits, progressivement plus épaisses en s'éloignant radialement de la paroi du puits (direction r(m)); et s - la Fig.18 est un tableau montrant la durée d'application t(d) exprimée en jours, d'une pression de fond de puits P(bar) imposée. DESCRIPTION DETAILLEE DE LA METHODE 1) Acquisition de données de laboratoire s Les essais d'endommagement des formations sont d'une importance primordiale pour minimiser ou pour réduire la détérioration de la perméabilité au voisinage des puits en choisissant le fluide de forage le mieux approprié et en optimisant les techniques de nettoyage des puits. Au cours des cinq dernières années, le demandeur a mis au point un équipement spécifique d'essai en laboratoire et des procédures destinées à caractériser l'endommagement des formations dû au forage au cours d'opérations en surpression et à quantifier les performances de différentes techniques de nettoyage mises en _uvre dans l'industrie comme le témoignent les publications suivantes: Alfenore, J. et al, "What Really Matters in our Quest of Minimizing Formation Damage in Open Hole Horizontal Wells ", 1999, SPE 54731; Longeron, D. et al, " Experimental Approach to Characterize Drilling Mud Invasion, Formation Damage and Cleanup Efficiency in Horizontal Wells with Openhole Completions " 2000, SPE 58737; ou Longeron, D. et al, " an Integrated Experimental Approach for Evaluating Formation Damages due to Drilling and Completion Fluids ", 1995, SPE 30089. Les essais de pression de fuite sont effectués avec une cellule de filtration dynamique qui peut recevoir des carottes d'un diamètre de 5 cm et d'une longueur pouvant atteindre cm. La cellule est équipée par exemple de cinq prises de pression situces à 5, 10, 15, 20 et 25 cm de la face d'entrée de la carotte. Les prises de pression permettent de surveiller les chutes de pression à travers six sections de la carotte pendant que l'on fait circuler la boue 2s et circuler en retour l'huile afin de simuler la production du puits. En vue de reproduire le processus dynamique d' invasion de la boue et du filtrat de la boue, les essais de laboratoire sont réalisés dans des conditions de puits représentatives (température, surpression et taux de cisaillement appliqué à la boue, carottes saturées en huile et en eau connce, etc.). De l'huile est ensuite injectée en sens inverse (courant de retour) à débit constant afin de simuler la production du puits. L'évolution des perméabilités de retour est calculée, pour chaque section, en fonction du volume cumulé d'huile injectée. La valeur finale stabilisée de la perméabilité de retour est ensuite comparée à la perméabilité initiale non détériorée afin d'évaluer la détérioration résiduelle en fonction de la distance par rapport à la face s d'entrée de la carotte. D'une manière générale, on a observé qu'une quantité totale de l0 à PV (maximum une centaine de PV) d'huile injectée était suffisante pour obtenir une valeur stabilisée de la perméabilité de retour après endommagement avec une boue à base d'huile. II - Modèle numérique simplifié pour supprimer les dommages au voisinage du i0 puits Considérant un puits foré dans la zone pétrolière avec une boue à base d'huile, on admet que les propriétés de l'huile dans le réservoir sont identiques à celles observées dans le filtrat. L'équation de l'écoulement au voisinage du puits est ainsi régie par une équation monophasique qui s'écrit comme suit: | - div(-Vp) = c) 3t (1) o p représente la pression, k, la perméabilité absolue, Il. la viscosité, c, la compressibilité et l, la porosité. La viscosité 11 et la compressibilité c dans le filtrat sont considérées comme identiques à celles observoes dans l'huile qui sature le réservoir. On considère que la pression initiale dans le réservoir est hydrostatique au début de la production. II 1 Modélisation du cake de filtration interne Le cake de filtration interne réduit la perméabilité du réservoir au voisinage du puits. Comme nous l'avons indiqué plus haut, les réductions de la perméabilité à l'issue de la période de forage et à la fin d'un nettoyage complet peuvent étre obtenues à partir de mesures effectuées en laboratoire. Pour la modélisation, nous utilisons le facteur de réduction de la perméabilité sous forme adimensionnelle pour représenter la variation de la perméabilité. L'utilisation de la forme adimensionnelle présente l'avantage de permettre de regrouper les données par zones géologiques. Soient kj la perméabilité initiale, kd la perméabilité après endommagement et kf la perméabilité finale de retour; la perméabilité après endommagement et la perméabilité s finale de retour sont généralement fonction de r la distance au puits. c(r) = 4() et c2(r) =-étant les courbes du facteur de réduction de la perméabilité kj ki en fonction de r avant le nettoyage et après le retour de fluide respectivement (Figure 1), la variation de la perméabilité au voisinage du puits est généralement limitée par ces deux courbes au cours de la période de retour de fluide. c,(r) correspond à la courbe de 0 perméabilité après endommagement et c2(r) correspond à la courbe de perméabilité de retour stabilisée. Comme nous l'avons indiqué plus haut, la variation de la perméabilité dans la zone occupée par le cake de filtration interne au cours de la période de retour de fluide dépend de la quantité d'huile produite qui s'écoule vers le puits. Nous utilisons la forme S adimensionnelle suivante pour décrire cette variation (Figure 2): C (Q)= k Q) k (2) o Q est le débit d'écoulement total à travers le milieu poreux dans le sens de l'écoulement divisé par la surface poreuse (surface des pores offerte à l'écoulement). Cette courbe représente la variation de la perméabilité par rapport à l'écoulement à travers une unité de surface poreuse. Elle correspond généralement à un sens d'écoulement donné. En pratique, le sens de l'écoulement est la direction radiale vers le puits. Lorsque Q = 0, il n'y a aucun écoulement permettant de nettoyer le cake de filtration, la perméabilité correspond à la perméabilité après endommagement avec k(O)=kd. Nous avons donc c0(0) = 0. Lorsque Q est très grand, le cake de filtration est entièrement nettoyé, la perméabilité correspond à 2s la perméabilité de retour finale avec k(+oo)=kf. Dans ce cas, nous avons Co(+oo)=1. La courbe de variation de la perméabilité peut étre mesurce à partir de données de laboratoire et elle peut être considérce comme indépendante de l'emplacement dans une carotte. Ainsi, une courbe est utilisce pour chaque zone géologique. Cette courbe est monotone. Son maximum est généralement atteint pour plusieurs m3 (ou plusieurs dizaines de m3) de fluide traversé par unité de surface poreuse. La perméabilité k à la distance r du puits au cours de la période de retour de fluide peut s'écrire sous la forme triviale suivante: k(r,Q)=(kf(r)-k,(r)) rk () k () +k,l(r) (3) En utilisant les courbes adimensionnelles définies ci-dessus et en prenant en compte lO l'équation (2), le facteur de réduction de la perméabilité c(r,Q) est exprimé par: c(r, Q) = (C2 (r) - c (r))CO (Q) + c (r) (4) Initialement, Q = 0, la réduction de la perméabilité correspond à celle obtenue après invasion du filtrat (perméabilité après endommagement) : c(r,O) = c,(r) (5) 1S A la fin du retour de fluide, lorsque la quantité de fluide en écoulement Q est très importante avec co(Q) 1, la réduction de la perméabilité correspond à l'état restauré avec la perméabilité de retour stabilisée: c(r,Q) = c2(r) (6) La variation de la perméabilité dans la zone occupée par le cake de filtration interne est modélisce avec l'équation (3). Contrairement au cake de filtration interne, l'impact du cake de filtration externe décrit ci-après est modélisé sous la forme d'un coefficient de paroi dans le modèle numérique discrétisé. II.2 Maillace et schémas numériques Un maillage cylindrique rHx est utilisé pour la modélisation de l'écoulement du fluide au voisinage d'un puits horizontal (Figure 3): r est la direction radiale, perpendiculaire à l'axe du puits, H est la direction angulaire et x est la direction le long du s puits. Avec ce maillage, les limites du puits sont discrétisces et des mailles de très petite taille peuvent être utilisées pour discrétiser la zone occupée par le cake de filtration interne. En général, le rayon du puits est de l'ordre de quelques centimètres, et l'épaisseur du cake de filtration interne varie entre quelques centimètres et quelques décimètres. En vue d' obtenir une bonne description du phénomène d' élimination du cake de filtration, les 0 mailles utilisées au voisinage du puits varient entre quelques millimètres et quelques centimètres. Pour des mailles cylindriques, un schéma standard numérique pour l' approximation du flux entre deux points peut être utilisé pour modéliser l'écoulement. Par exemple, le flux entre deux mailles voisines i et i+1 dans la direction radiale est calculé par (Figure 4): Fi+l,2 = T,+l,2(pi+l - Pi) () avec: Ti+l,2= 1 l l jk (8) ln ri+l'2 + ln ri+1 krti ri krti+l ri+l/2 o j et k sont les indices des mailles considérées dans les directions O et x, r représente la distance de la maille i au puits, rj+"2 est la distance de l'interface des mailles considérées au puits, krj est la perméabilité de la maille i dans la direction radiale, AU et x sont les longueurs des mailles dans les directions (3 et x et Ti la transmissivité entre mailles. Le terme " mailles de puits " désignant les mailles qui discrétisent les limites du puits, les conditions aux limites du puits sont traitées dans les mailles de puits. La pression interne du puits pw et le débit du puits q; sur une maille i considérée peuvent être reliés par la formule de discrétisation suivante (Figure Sa): q; = IPj(pi - pw) () avec: Plj= r jxk (10) ln.v krl r o rw est le rayon du puits. Cette discrétisation aux limites du puits est similaire à l' approximation de l'écoulement du fluide entre deux mailles. Cependant, pour la discrétisation des limites du puits, le coefficient de discrétisation est désigné par l' index de productivité numérique IP et non par la transmissivité T. et le flux F est remplacé par le déhit du puits qi. Cette notation est cohérente par rapport au modèle numérique de puits couramment utilisé, et le coefficient de paroi peut être intégré dans le terme de l'index de productivité numérique IP. La perméabilité kr,; varie au cours du retour de fluide dans la zone occupée par le 0 cake de filtration interne selon la formule présentée dans la section précédente. Ainsi, la transmissivité et l' index de productivité numérique IP varient également dans la simulation au cours de la période de retour de fluide. II.3 Modélisation du cake de filtration externe La présence du cake de filtration externe peut étre prise en compte dans la formule de discrétisation via l' indice numérique IP. Dans le cas de la présence d' un cake de filtration externe d'épaisseur de et de perméabilité ke' la pression du puits Pw correspond à la pression sur le rayon rw - de et non pas sur le rayon rw. La chute de pression est élevée à travers le cake de filtration externe qui se trouve dans la zone situce entre rw - de et rw. En utilisant une nouvelle fois l'équation (9) pour relier la pression du puits pw, la pression des mailles de puits pi et le débit du puits qi, le coefficient de discrétisation IP devrait intégrer l' impact du cake de filtration externe comme suit (Figure Sb): -ln-±ln w (11) k, i ri ke rwde On admet que le cake de filtration externe est élirniné si la différence de pression au travers de son épaisseur est supérieure à un seuil donné. Ainsi, au début du retour de fluide, le coefficient numérique IP est calculé en utilisant l'équation (11) qui intogre la présence du cake externe si ce dernier existe. Une fois que la différence de pression au travers du cake de filtration est supérieure au seuil donné, l'index de productivité numérique IP est calculé avec l'équation (10). La perméabilité ke du cake de filtration externe pourrait généralement être très s inférieure à la perméabilité au sein du réservoir ou dans la zone occupée par le cake de filtration interne. Ainsi, en présence du cake de filtration externe, le coefficient numérique IP est très petit. Pour réaliser les simulations on peut utiliser un outil de modélisation d'écoulements tel que le modèle ATHOS par exemple (ATHOS est un modèle de modélisation numérique 0 développé par 1'IFP). Le schéma de discrétisation utilisé est un schéma à 5-points classique pour modéliser l'équation de diffusivité en maillage cylindrique. Dans les mailles au voisinage immédiat du puits, un IP numérique est utilisé pour relier la pression dans ces mailles, la pression au fond du puits et le débit d'écoulement vers le puits. Comme la perméabilité au voisinage du puits change pendant la période de dégorgement, les transmissivités autour du puits et 1'IP changent aussi en fonction de la variation des perméabilités. Les courbes qui définissent les coefficients multiplicateurs de perméabilités en fonction de la distance au puits, cj(r) et c2(r), sont entrées dans le simulateur sous forme de tableaux de valeurs. Les valeurs correspondantes dans chaque maille sont calculées à partir de ces courbes en utilisant une interpolation linéaire comme expliqué précédemment. Le volume poreux cumulé de fluide traversant une interface entre deux mailles dans la direction radiale r est utilisé pour calculer le coefficient multiplicateur de transmissivité entre ces deux mailles à chaque instant considéré. III Résultats numériques Nous présentons deux exemples pour illustrer les capacités de la méthode qui a été développée: le premier concerne le dégorgement d'un cake interne sans cake externe et le deuxième le dégorgement en présence d'un cake interne et également d'un cake externe. Exemple 1: Dégorgement en présence du cake interne seul On considère une partie d'un puits horizontal d'une longueur de 20 m, qui traverse 4 zones représentatives alternativement de deux types d'hétérogénéités différentes (Figure 6). Les perméabilités k des milieux correspondants, initialement sans aucun endommagement, sont de 1000 et 100 mD. La longueur de chaque milieu traversé est de 5 s m. Les valeurs de la perméabilité dans les mailles o s'est formé le cake interne du à l'endommagement sont entrces manuellement dans le jeu de donnces. Les courbes, par zones, du coefficient multiplicateur de la perméabilité endommagée en fonction de la distance à la paroi du puits c(r) sont données dans la Figure 7. Les courbes de la perméabilité retour c2(r) sont présentées dans la Figure 8. Ces courbes sont discontinues 0 car les données fournies par laboratoire sur seulement quelques points. Plus le nombre de points est important, et mieux la courbe de laboratoire est représentée. La variation de la perméabilité pendant le nettoyage en fonction de la quantité de fluide écoulé à l'unité de surface poreuse, c0(V), est présentée à la Figure 9. En pratique, le plateau maximum peut étre atteint avec quelques mètres cubes de fluide passé par unité surface. Comme on l'a vu, un maillage cylindrique est utilisé pour les simulations. Le réservoir est très grand dans la direction radiale avec un rayon extérieur de 1750 m o la condition limite est de flux nul. Sur les frontières aux deux extrémités du puits, la condition est aussi de flux nul. Les nombres et les tailles de mailles dans les directions r et x sont donnés dans la Fig.17 (0=360 ). Le puits est discrétisé en 80 mailles suivant sa longueur. Chaque zone de perméabilité constante est ainsi discrétisée en 20 mailles de 0.25 m. La pression initiale dans le réservoir à la cote du puits est sensiblement de 320 bar. Deux simulations ont été faites avec différentes conditions imposées au puits: a). Un débit de 20m3/j est imposé au puits pendant 1.5 jours. L'écoulement au voisinage du puits simulé avec la méthode présentée cidessus tenant compte de la 2s variation de la perméabilité en cours du temps est notée SM1. Cette simulation est comparée à deux autres simulations en utilisant le modèle d'écoulement classique avec des perméabilités inchangées, égales d'une part aux perméabilités après endommagement c:(r), et d'autre part aux perméabilités retour c2(r). Ces deux simulations sont notées SM2 et SM3. Les résultats de simulations sont présentés pour les mailles des puits 31 et 40 situces au milieu et à la limite d'une des zones peu perméables, et pour les mailles 41 et 50 situées en limite et au milieu de la zone plus perméable suivante. La Figure 10 montre les déhits d'huile au niveau de ces mailles pour les trois scénarios simulés: SM1, SM2 et s SM3. Les simulations avec des perméabilités fixées, SM2 et SM3, donnent des débits constants par maille, ce qui est normal puisque la limite dans la direction r n'est pas atteinte pour le temps court simulé (1.5 jours). Par contre, les débits varient quand on modélise les variations de perméabilités dans le cake interne pendant la remise en production. Au temps 0, ces déhits sont identiques à ceux des simulations avec les 0 perméabilités consécutives à l'endommagement. Ensuite, ils se différencient en raison des perméabilités qui augmentent dans le cake interne suite au nettoyage par l'huile de la formation. Très rapidement, au bout d'une journée, ces débits rejoignent ceux simulés avec les perméabilités retour. Les variations de perméabilité dans les mailles 31 et 50 sont présentées aux Fig.lla, llb respectivement. Ces variations correspondent à celles dans les deux zones. Les perméabilités dans les états d'endommagement et de retour sont aussi représentées. La variation de perméabilité au cours du nettoyage est bornée par ces valeurs. Au bout d'un jour, la perméabilité dans la zone la plus perméable (maille 50) rejoint presque la valeur de perméabilité retour, et celle dans la zone moins perméable (maille 31) ne change pas beaucoup. Mais, comme la variation entre la perméabilité d'endommagement et celle de retour est très faible dans la zone peu perméable, les résultats de simulation dépendent principalement de la variation de perméabilité dans la zone la plus perméable. Sur les résultats présentés à la Fig.10, les détits augmentent dans les zones plus perméables et ils rejoignent très rapidement ceux de la simulation SM3. Les débits dans les zones peu 2s perméables diminuent car les simulations sont faites avec un débit total imposé au puits. Avec cette modélisation, nous pouvons également obtenir la variation de vitesse locale due au dégorgement du cake. b). Une différence de pression de 1 bar est imposée pendant 1.5 jour. La Fig.12 montre la variation en fonction du temps t exprimé en jours, des débits simulés FR (exprimés en m3/j) correspondants au puits. Dans le cas de perméabilité inchangée (SM2 et SM3), les débits diminuent avec le temps. Par contre, la modélisation d'un dogorgement progressif donne un débit croissant jusqu'à un jour environ, avant de diminuer. L' augmentation de déhit pendant la période initiale est due à l'augmentation de perméabilité dans le cake interne au cours de la remise en production. Les résultats dans les mailles de puits 31, 40, 41 et 50 sont trèssimilaire aux ceux du cas a. Les débits avec la modélisation du nettoyage de cake au temps t=0 sont égaux à ceux simulés avec les perméabilités d'endommagement, et ensuite ils varient et rejoignent les débits simulés avec les perméabilité de retour. Dans cet exemple, nous observons que le nettoyage du puits est assez rapide quel que 0 soit le scénario modélisé. Dans tous les cas, les résultats de simulation du dégorgement progressif SM1 sont très proches au bout d'un jour de ceux obtenus avec les perméabilités retour SM3. I1 est possible de fournir les détails des résultats aux temps courts, comme par exemple les déhits le long du puits, les pressions et les vitesses au voisinage du puits, pour mieux connâître ce qui se passe au cours du dégorgement. Toutefois, les performances du puits aux temps longs, après quelques jours, sont quasi identiques, quelles que soient les configurations étudiées, sachant que les aspects géomécaniques ne sont pas pris en compte. Dans cette hypothèse, il semble donc que les impacts du cake interne sur la performance du puits soient très limités en temps et qu' il suffise en général d' étudier cette performance en considérant la perméabilité restaurée, c'est-à-dire celle de la configuration notée SM3. Exemple 2. Présence d'un cake externe non uniforme le long du drain On considère la même géométrie de puits que dans l'exemple précédent. Dans cet exemple, le réservoir est homogène avec une perméabilité de 1000mD dans le milieu poreux. Le cake externe n'a pas une présence homogène le long du puits. Dans certains endroits, il n'y pas de cake externe, et dans les endroits o le cake externe est présent, 2s celui-ci a une perméabilité kext de 1 mD et une épaisseur reX de 4 mm comme dans l'exemple précédent. La distribution de la présence du cake externe est donnée dans la Fig.13. La différence de pression nécessaire pour l'arrachement du cake externe est toujours fixée à 0.5 bar. Deux types de conditions aux limites sont utilisés dans les simulations. Pour le premier cas, une pression de 318.2 bar est imposce au fond du puits, soit une différence de pression entre le réservoir et le puits de 1.8 bar. Pour le deuxième cas, nous imposons plusieurs paliers consécutifs de pression pour arriver à une chute totale de pression de 1.8 s bar (Tableau 2). Les Figures 14 et 15 présentent la distribution du cake externe et la répartition du débit le long du puits pour ces deux cas à différents temps de production. Dans le premier cas, les débits sont uniformes le long du puits, car les cakes externes sont entièrement arrachés dès le début. Dans le deuxième cas, la répartition de déhit varie en fonction du 0 temps, car les cakes externes sont arrachés non uniformément à des temps différents. De plus, il y a touj ours des cakes externes qui ne peuvent pas être arrachés au bout de 5 jours. La Fig.16 montre la production du puits pour ces deux cas. Dans le premier cas, la production du puits est plus forte, car tous les cakes externes sont arrachés dès le début. Mais, le maximum du débit local le long du puits est toujours inférieur à 3m3/m jour. Dans le deuxième cas, le débit du puits est plus faible, mais le détit local peut être très fort avec un maximum de 4.5m3/mjour. Les cakes ne peuvent toujours pas étre arrachés à certains endroits. La performance du puits est largement réduite dans ce cas. Cet exemple montre que les procébures de dégorgement peuvent influencer la performance du puits même dans un réservoir homogène. Bien qu'il paraisse tentant de vouloir imposer une forte différence de pression entre le puits et la formation, puisque c'est la procédure qui permet d'arracher le plus vite et le plus uniformément le cake externe qui limite le débit du puits, il peut étre dangereux pour l'intégrité du puits d'agir ainsi si la formation est non consolidée et que des venues de sable risquent de se produire avec comme conséquence le bouchage du puits. C' est l'un 2s des intérêts de la présente invention que de pouvoir définir la meilleure procédure de dégorgement du puits sans provoquer le risque évoqué dès lors que l'on connâît la vitesse du fluide à partir de laquelle le sable perd sa cohésion. REVENDICATIONS l) Méthode pour simuler les conditions optimales à imposer dans un puits foré au travers d'un gisement souterrain à une trajectoire quelconque, pour l'élimination progressive par des fluides issus du gisement, des dépôts ou cakes qui se sont formés dans au moins une zone à la périphérie du puits, suite aux opérations de forage et de complétion, caractérisée en ce qu'elle comporte: - l' acquisition de données initiales obtenues par des mesures de laboratoire de l'épaisseur et des cakes ainsi que des valeurs de perméabilité endommagée (kd) et de perméabilité restaurée (kf) de la zone entourant le puits, en fonction de la distance (r) à la 0 paroi du puits, suivant la valeur de la perméabilité initiale (ki) de la formation autour du puits; - la discrétisation de la zone endommagée par un maillage cylindrique en 3D formant des blocs d'épaisseurs radiales petites relativement au diamètre du puits; et - la résolution dans ce maillage de l'équation de diffusivité modélisant les lS écoulements des fluides au travers des cakes en tenant compte des données initiales mesurées et en modélisant l'évolution de la perméabilité en fonction des débits (Q) de fluides s'écoulant au travers des cakes, pour en déduire les conditions optimales à appliquer pour la mise en production du puits. 2) Méthode selon la revendication 1, caractérisée en ce que l'on modélise la restauration de la perméabilité en tout point à distance (r) de la paroi en considérant que la perméabilité varie proportionnellement à l'écart entre la perméabilité endommagée (kd) et la perméabilité restaurce (kf), le coefficient de proportionnalité dépendant d'une loi empirique de variation de la perméabilité en fonction de la quantité (Q) de fluides au