FR 2538849 A1 19840706 FR 8222100 A 19821230 PROCEDE ET DISPOSITIF POUR DETERMINER LES CARACTERISTIQUES D'ECOULEMENT D'UN FLUIDE DANS UN PUITS A PARTIR DE MESURES DE TEMPERATURE L'invention concerne l'étude des caractéristiques d'écoulement d'un fluide dans un puits, plus particulièrement dans les zones du puits o une circulation de fluide se produit entre le puits et les formations géologiques traversées par le puits, telles que les zones productrices dans les puits produc- teurs d'hydrocarbures. Le débit du fluide s'écoulant dans un puits peut être mesuré au moyen d'une sonde de débit suspendue à l'extrémité d'un câble En déplaçant la sonde le long du puits, on obtient une diagraphie de débit indiquant le débit en fonction de la profondeur Dans le cas d'un puits producteur, l'examen de cette diagraphie permet de déterminer le débit des fluides, le débit des fluides entrant dans le puits dans les différentes zones productrices, et donc d'évaluer les rendements respectifs de ces différentes zones. Les sondes de débit utilisées habituellement sont des dispositifs à turbine qui ne fournissent des mesures précises que si l'écoulement est stable Or l'entrée dans le puits des fluides provenant de la formation perturbe considérablement l'écoulement, les perforations pratiquées dans le tubage du puits pour permettre l'entrée des fluides étant dirigées radia- lement par rapport à l'axe du puits On constate de façon générale que les sondes de débit ne réagissent qu'avec un certain retard à l'entrée de fluides dans le puits. Lorsque la séparation entre zones productrices consé- cutives est suffisante pour que l'écoulement puisse se stabili- ser, la sonde de débit fournit des mesures précises, et le débit des différentes zones productrices peut être évalué correctement Mais si la séparation entre zones productrices est insuffisante pour permettre à l'écoulement de se stabili- ser, les mesures fournies par la sonde de débit sont imprécises et ne permettent pas en tout cas d'évaluer individuellement le rendement des différentes zones productrices. Une autre approche repose sur la relation qui existe entre le débit et la température du fluide dans le puits. Considérons le cas d'un puits producteur de pétrole La tempéra- ture des formations, et donc du liquide de formation, décroît avec la profondeur selon le gradient géothermique Si l'on suppose que l'équilibre thermique est réalisé entre le liquide dans le puits et les formations environnantes, la température dans le puits suit la même variation géothermique Mais du fait de l'inertie thermique du fluide entrant dans le puits et s'écoulant vers la surface, la température dans le puits au-dessus du point d'entrée est portée à une température supé- rieure à la température correspondant à la variation géothermi- que la décroissance de la température est donc plus lente que le gradient géothermique et une certaine distance est nécessai- re pour qu'elle revienne à la valeur géothermique - le brevet US 3 913 398 (Curtis) utilise un modèle selon lequel la courbe de température au-dessus d'une zone de production se rapproche de façon exponentielle d'une asymptote parallèle au gradient géothermique la distance de relaxation intervenant dans l'exponentielle est un paramètre qui dépend de la densité, de la chaleur spécifique et du débit massique du fluide ainsi que du temps pendant lequel le puits a produit. l'examen du profil de température permet d'identifier les intervalles présentant l'allure d'une exponentielle On ajuste par itération la distance de relaxation jusqu'à ce que le profil de température déterminé sur la base du modèle corresponde au profil mesuré La valeur obtenue pour la distance de relaxation est alors utilisée pour calculer le débit massique et, connaissant la densité du fluide, le débit volumétrique. 3 2538849 Cette technique est intéressante en combinaison avec une mesure de débit en particulier pour déceler -les écoulements de fluide à l'extérieur du tubage du puits La sonde de débit n'est en effet sensible qu'à l'écoulement dans le puits, alors que la mesure de température est influencée par tous les types d'écoulement, qu'ils aient lieu à l'intérieur ou à l'extérieur du puits. Toutefois, cette technique, repose sur des hypothèses qui s'écartent souvent de façon notable des conditions réelles. Elle suppose en particulier que la séparation entre zones productrices soit assez grande, au minimum d'une trentaine de mètres, et que sur l'intervalle qui part d'un point d'entrée de fluide, le gradient géothermique et les propriétés de transfert thermique soient constants En fait, le gradient géothermique ne reste constant que sur de faibles intervalles correspondant à une formation uniforme, et les caractéristiques de transfert thermique peuvent également varier avec l'état du puits (diamètre du tubage, qualité de -la cimentation entre le tubage et la formation. Il est fréquent en pratique que plusieurs zones productives d'épaisseur limitée se succèdent avec des interval- les très courts, de l'ordre de quelques mètres Dans ce cas, la technique du brevet cité fournira le débit total de l'ensemble des zones productives, mais ne permettra pas d'évaluer les contributions respectives des différentes zones. Selon l'invention, on mesure avec une sonde de diagra- phie la température le long d'un intervalle de profondeur comprenant des zones de circulation de fluide entre le puits et les formations, on mesure le débit aux extrémités de cet inter- valle, et on détermine pour chaque profondeur, à partir des températures mesurées et desdits débits, une valeur de débit, l'ensemble des valeurs ainsi obtenues étant enregistré sous la forme d'un profil de débit synthétique. L'invention-repose sur une relation entre la tempéra- ture T à une profondeur z et le débit à la 'même profondeur. Cette relation est déduite de l'analyse des transferts thermiques sur un petit intervalle dz, et consiste en une équation différentielle du premier ordre reliant le débit et sa dérivée à des grandeurs pouvant être déduites du profil de 4 25388494 température, à savoir la température et sa dérivée, et à la température Te -du fluide dans la formation, cette température Te variant comme le gradient géothermique. De façon appropriée, on résout par itération l'équation différentielle, la grandeur itérante étant la température Te, en imposant des solutions au sommet et au bas de l'intervalle de profondeur considéré égales aux débits mesurés à ces profondeurs On obtient ainsi une valeur de débit à chaque profondeur dudit intervalle La valeur de convergence 1 o de la température Te permet en outre de déduire le gradient géothermique, à titre d'information complémentaire. le profil de débit synthétique ainsi obtenu permet d'analyser avec précision les contributions des différentes zones productrices dans le cas d'un puits producteur Il remédie par conséquent aux insuffisances de la sonde de débit dans les zones perturbées. Une manière simple de mesurer le débit est de combi- ner une sonde de débit à la sonde de température. L'invention sera bien comprise, à la lecture de la description ci-après, faite en référence aux dessins annexés. Dans les dessins: la figure 1 représente schématiquement un appareil -de diagraphie dans un puits producteur de pétrole, la figure 2 illustre sous forme de diagramme le traitement permettant d'obtenir un profil de débit synthétique à partir des mesures fournies par l'appareil de la figure 1, les figures 3 a, 3 b donnent des exemples de tels profils de débit synthétiques comparés aux profils de débit mesurés. La figure 1 représente un forage 10 traversant des formations 11 et dans lequel on a placé un tubage 12 L'espace entre le tubage et les parois du forage a été rempli avec du ciment 13 Une série de perforations telles que 15 ont été pratiquées dans une zone productrice 16 pour permettre au fluide présent dans cette zone de s'écouler dans le puits De même, on a représenté une seconde zone productrice 17 et la série de perforations 18 correspondante. Au-dessus des zones productrices est disposée à l'intérieur du tubage une colonne de production 20, avec un 2538849, bouchon (packer) 21 fermant l'espace 22 entre la colonne de production et le tubage pour contraindre le fluide produit à s'écouler à l'intérieur de la colonne de production. La figure 1 montre schématiquement un appareil de diagraphie comprenant un dispositif de fond 30 suspendu à un câble électrique 31 et reliée par ce câble à un équipement de surface 32 L'équipement de surface comprend un treuil, non représenté, pour l'enroulement du câble, et fournit la puissan- ce électrique et les signaux de commande nécessaires au lo fonctionnement du dispositif de fond Les signaux de mesure fournis par le dispositif de fond sont enregistrés au moyen d'un enregistreur à bande magnétique 33 et d'un enregistreur graphique 34 associés à l'équipement de surface De façon classique, des signaux indicatifs de la profondeur du dispositif de fond sont fournis par un dispositif 29 associé au treuil, qui mesure le déplacement du câble Ces signaux permettent d'associer une profondeur à chaque valeur de mesure produite par le dispositif de fond. Le dispositif de fond comprend un ensemble de sondes placées bout à bout et sensibles à différentes caractéristi- ques, et une cartouche de télémétrie 35 assurant l'interface entre les sondes et le câble 31. Un ensemble de sondes typique comprend une sonde de débit 36 sensible à la vitesse d'écoulement du fluide, par exemple un débitmètre à turbine tel que décrit dans le brevet US 3 954 006, une sonde 37 sensible à la densité du fluide telle que le gradiomanomètre décrit dans le brevet US 3 455 157, une sonde de température 38 et un dispositif 39 détectant les joints du tubage. Il est connu que les sondes de débit ne fournissent des mesures précises que si l'écoulement est relativement stable En face des perforations, les mesures de débit sont faussées par les perturbations dues à l'entrée radiale des fluides dans le puits Des exemples en sont fournis par les diagrammes de débit des figures 3 a, 3 b Sur toutes ces figures, la profondeur est en ordonnée et le débit est en abscisse, le profil de débit tel qu'il est fourni par la sonde de débit étant représenté en trait plein. 6 2538849, Sur tous ces diagrammes, on a indiqué les profondeurs en pieds, et chaque graduation représente deux pieds, soit environ 63 cm. Sur la figure 3 a, on note deux zones productives A 1 et A 2 La zone A 1 s'étend de la cote 9670 à la cote 9600, c'est-à-dire sur environ 22 m La zone A 2 s'étend de la cote 9550 à la cote 9590, donc sur 6,5 m On note que la sonde de débit ne réagit qu'avec un retard très net à l'entrée de fluide par les perforations Dans le cas de la zone A 2, on n'observe un accroissement significatif du débit que près du sommet de cette zone,-environ 5 m au-dessus du bas de la zone. La figure 3 b montre un autre exemple, avec trois zones productives B 1, B 2, B 3 séparées par des intervalles très courts La production de la zone inférieure Bl n'est prise en compte qu'une fois dépassé le sommet de cette zone Ce retard peut être imputable en partie à un calage en profondeur défec- tueux du profil de débit Mais quoi qu'il en soit, le profil de débit-de la figure 36 ne permet pas de déterminer les contribu- tions respectives des zones B 1, B 2 et B 3. La relation, indiquée plus haut, entre le débit et la température du fluide qui s'écoule dans le puits est mise à profit -de la manière expliquée ci-après pour produire un profil de débit synthétique reflétant de façon plus précise et sans retard les échanges de fluide entre le puits et les formations. Considérons un volume dans le puits oh un échange de fluide se produit entre le puits et la formation. Appelons z la profondeur, T la température du fluide qui s'écoule, Te la température du fluide dans la formation, Tf la température du fluide entrant dans le puits, Th la tempéra- ' turenà l'interface entre la formation et le ciment entourant le tubage du puits, et u le débit. Le principe de conservation de l'énergie conduit à la relation: d T 4 U T -Th d T 1 du =Dh+ (T Tf)- dz Dpop u dzf U () Dans cette relation U signifie le coefficient de transfert thermique entre le puits et la formation 7 2538849 D est le diamètre intérieur du tubage, pest la densité du fluide Cp est sa capacité calorifique Le terme 4 U T Th Dpc u correspond à l'échange thermique entre le fluide qui s'écoule et la formation, dû à la différence entre la température T du fluide et la température Th- d T h Le terme dépend de la phase du fluide. dzf Son expression complète est d T l 1 d P d F Vdjz/ = TL gp dz +Udû ( 2) Le premier terme est imputable aux effets de dilation du fluide Il est négligeable pour un liquide, et peut être déduit de l'équation d'état pour un gaz. Le second terme d F/dz représente la chaleur résul- tant du frottement du fluide sur les parois du tubage Ce terme vaut d P 1 fu 2 dz 2 À D ( 3) oh f est le coefficient de frottement et D le diamètre inté- rieur du tubage. Le troisième terme (T Tf)1 du u dz exprime la variation de température due à l'échange de fluide entre la formation et le puits Ce terme est nul dans les zones stationnaires oh aucun échange de fluide ne se produit, car alors du/dz est nul. Dans le cas o le fluide produit est un liquide, la température d'entrée Tf peut être remplacée par la température Te dans la formation - Dans le cas o il-s'agit d'un gaz, les températures Te et Tf sont habituellement différentes, car le gaz se dilate à l'entrée dans le puits du fait de la différence AP entre la pression de la formation et la pression du puits, et cette dila- tation s'accompagne d'une variation de température par effet de JouleThomson On a alors Tf T = e P ( 4) le coefficient de Joule-Thompson Y A est habituelle- ment positif, auquel cas le gaz se refroidit en entrant dans le puits. la résolution de l'équation différentielle ( 1) en vue de déterminer, pour chaque profondeur une valeur de débit u à partir de la température mesurée T, exige la connaissance d'un grand nombre de paramètres dont beaucoup sont en fait difficiles à évaluer avec une précision suffisante. Ainsi, le coefficient de transfert thermique U, exprime la résistance au transfert thermique offerte par le fluide dans le puits, le tubage et la couche de -ciment Il est donc lui-même fonction de nombreux paramètres (conductivité thermique et épaisseur des différentes couches). La température Th à l'interface ciment-formation diffère de la température Te du fluide dans la formation entre le fluide et la formation L'équilibre thermique est très long à s'établir après la mise en production du puits Par consé- quent, la différence Th Te dépend du temps pendant lequel le puits a produit, et également de la diffusivité et de la conduc- tivité thermiques de la formation, quantités mal connues. Une autre difficulté réside dans la détermination de l'équation d'état, lorsque le fluide est un gaz ou un mélange gaz-liquide Les équations réellement représentatives pour les hydrocarbures sont des équations empiriques qui exigent la connaissance de la composition du mélange. Une approche qui permet d'éliminer tous ces para- mètres difficiles à estimer consiste à considérer globalement 2538849, les deux premiers termes de l'équation ( 1) comme la partie (d T/dz), du gradient de température total non imputable aux échanges de fluide Dans une zone stationnaire, ainsi qu'on l'a vu, le troisième terme de l'équation ( 1) est nul (du/dz = o); la somme des deux premiers termes est donc égale au gradient de température total dans une zone stationnaire Cette quantité peut être déduite facilement du profil de température mesuré, une fois identifié une zone o le gradient de température est constant sur un intervalle de profondeur suffisant. L'équation ( 1)devient alors, dans le cas d'une production de liquide (Tf = Te) (d PT (d T du -; -u T ( 5) dz T e Avec cette approche, la seule grandeur qui ne se déduise pas directement du profil de température est la tempé- rature Te du fluide dans la formation. e La température T e suit le gradient géothermique mais même si on peut supposer que celui-ci est constant sur l'inter- valle de profondeur considéré, sa valeur est mal connue La technique utilisée pour faire face à cette difficulté consiste 2 o à intégrer l'équation ( 5) par itération, la grandeur itérante étant Tes en imposant des solutions aux extrémités de l'inter- valle égales aux valeurs que fournit la sonde de débit pour les mêmes profondeurs. En résumé le traitement conduisant à un profil de débit synthétique nécessite en plus du profil de température, la connaissance des valeurs de débit aux extrémités de l'inter- valle de profondeur que l'on considère. Ce traitement, exécuté sur un calculateur numérique programmé, comprend les phases ci-après, schématisées sur la 3 e figure 2. Le profil de température mesuré et enregistré se compose d'une série de valeurs discrètes, avec un point de mesure tous les 15 cm. On procède tout d'abord à un filtrage (bloc 50) du profil de température pour éliminer des petites variations dues au bruit l Ie filtrage peut consister de façon appropriée à convoluer le profil de température chaque fois sur une fenêtre donnée, par exemple représentant 20 points de mesure, avec une fonction triangle. L'étape suivante est la dérivation (bloc 51) du profil de température filtré, qui fournit en chaque point la grandeur d T/dz (gradient de température total). lo L'analyse-du gradient permet de délimiter les zones stationnaires, o ce gradient est constant au bruit près et a une valeur nettement plus faible que devant les entrées de fluide Les zones d'entrée de fluide sont définies comme les zones non stationnaires. Pour cette analyse, on considère tout d'abord la zone s'étendant d'une profondeur ZA située à une certaine distance (par exemple 1,5 m) audessus de la plus haute des perforations jusqu'à l'extrémité supérieure, de profondeur zh' de l'intervalle dans lequel la température a été mesurée Les profondeurs des différentes perforations étant connues, la cote ZA est facile à déterminer. On calcule la moyenne M de la dérivée d T/dz sur cette zone (z A, Zh) (bloc 52). On détermine ensuite (bloc 53) l'écart-type S de la dérivée d T/dz sur tout l'intervalle de mesure, soit de la cote zb à la cote Zh, pour avoir une estimation du bruit sur l'ensem- ble de l'intervalle, et on calcule la somme D = M + S (bloc 54). On calcule la moyenne m de la dérivée d T/dz (bloc ) sur N points de mesure consécutifs (avec par exemple N = 4), et on compare cette moyenne m à la valeur D (bloc 56) Si m est inférieur à D, la zone correspondante est considérée comme appartenant à une zone stationnaire On procède à cette compa- raison pour tous les groupes de N points. La comparaison est faite en suivant l'ordre des pro- fondeurs, du haut vers le bas Pour chaque groupe d'ordre i, après la comparaison de sa moyenne mi ave D, on compare le signe de la différence mi D avec le signe de la différence mi i D relative au groupe précédent d'ordre i-i (bloc 57). 2538849, Si les signes sont les mêmes, les deux groupes appartiennent à la même zone (stationnaire ou non stationnai- re) S'ils sont différents, les deux groupes appartiennent à des zones différentes, et le groupe de signe positif est identi- fié comme une limite d'une zone non stationnaire Si le groupe i-1 donne un signe négatif et le groupe i un signe positif, la limite définira le front de montée d'un pic -sur la courbe f(z) = d T/dz, un tel pic étant imputable A une entrée de fluide. Dans le cas inverse, il s'agira du front descendant d'un tel pic Les limites ainsi déterminées, complétées d'un code "front de montée" ou "front de descente", sont enregistrées dans le bloc 58. Une zone d'entrée de fluide est identifiée lorsqu'un front de montée, puis un front de descente ont été trouvés (bloc 59). -Lorsqu'on analyse de cette façon la zone à l'extrémité inférieure de l'intervalle, on peut trouver un front de montée indiquant la fin d'une zone stationnaire, mais aucun front de descente Cela signifie que la zone en question suit le gradient géothermique, qui est plus important que le gradient dans les zones stationnaires situées entre des zones d'entrée de fluide ou au-dessus de ces zones Dans ce cas, on calcule (bloc 60) la moyenne de la dérivée d T/dz dans cette zone Cette moyenne peut être utilisée dans les phases ultérieu- res du traitement en tant que valeur du gradient géothermique. Lorsque ces opérations ont été effectuées sur l'ensemble de l'intervalle, qui se trouve ainsi fractionné en zones d'entrée de fluide et zones stationnaires (avec éventuel- lement une zone géothermique au bas de l'intervalle), on élimi- ne (bloc 61) les zones stationnaires trop fines, par exemple qui représentent moins de deux points de mesure ( 30 cm). L'identification des zones d'entrée de fluide et des zones stationnaires permet ensuite de calculer (bloc 70) pour chaque profondeur l'expression (de (=d Tz)d(T rdze = rdzl ja rdl 2538849, figurant au numérateur dans l'équation ( 5), expression qui correspond au gradient de température imputable aux seules entrées de fluide. On examine tout d'abord si la profondeur considérée se trouve dans une zone stationnaire ou dans une zone d'entrée de fluide, par comparaison avec les limites des zones telles que déterminées ci-dessus Si l'on se trouve dans une zone stationnaire, (d T/dz)e est en principe voisin de zéro On calcule sa valeur en prenant pour valeur de (d T/dz) la moyenne, sur la zone stationnaire'coneidér 4 e, du gradient (d T/dz) résultant de la mesure. Dans les zones d'entrée de fluide, on préfère prendre en-considération les moyennes de (d T/dz) dans les deux zones stationnaires entourant la zone d'entrée de fluide considérée, et prendre pour valeur de (d T/dz)a l'interpolation linéaire entre ces deux moyennes. Une fois calculé le terme (d T/dz)e sur l'ensemble de l'intervalle, il reste à intégrer l'équation différentielle ( 5). Comme expliqué plus haut, la façon dont la tempéra- ture Te entrant dans l'équation ( 5) varie avec la profondeur est mal connue, et on intègre l'équation différentielle ( 5) de façon itérative, avec Te comme grandeur itérante- L'intégration (bloc 72) fournit un ensemble de valeurs de débit u(z) pour chaque profondeur z correspondant à, une mesure de température. L'ensemble de valeurs retenu est celui pour lequel la valeur de débit uh au sommet de l'intervalle de profondeur considéré -est égale à la valeur U O mesurée par la sonde de débit au même niveau (ou du moins ne diffère de 1 Uo que d'une quantité inférieure à une valeur déterminée). La température Te peut être exprimée par Te = Tb az Tb étant la température du fluide dans la formation au bas de l'intervalle considéré et a le gradient géothermique dans cet intervalle. L'itération peut porter soit sur la température Tb, soit sur le gradient géothermique a. 2538849, On prend comme estimation Tb de Tb la température au sommet de la zone géothermique, identifiée par l'analyse expo- sée ci-dessus Si 1 tanalyse n'a pas permis de déceler une zone géothermique, on prend, à défaut, la température au bas de l'intervalle. Pour le gradient géothermique, on prend comme estima- tion la valeur déduite de l'analyse des zones ci-dessus, ou à défaut, une valeur donnée au départ qui repose sur les informa- tions géologiques disponibles. Lorsque l'itération porte sur la température Tb, on * fixe le gradient à la valeur a et Te s'écrit. * * Te = Tb D Tb a z l'itérant étant D Tb. Lorsque l'itération porte sur le gradient, Te s'écrit Te = Tb (a* + Da)z l'itérant étant Da. Le choix du mode d'itération est à fixer au départ. Si l'on appelle IT l'itérant, qu'il s'agisse de l'itération sur la température ou sur le gradient, le processus itératif consiste à arriver à la valeur I Tf telle que P(I Tf) = UO, o P est le résultat de l'intégration de l'équation ( 5) au sommet de l'intervalle de profondeur considéré. Il faut noter tout d'abord que la valeur ub au bas de l'intervalle ne doit pas être nulle, car le modèle décrit plus haut n'autorise pas une valeur de débit nulle Si, comme il est habituel, le bas de l'intervalle est plus profond que la zone de production la plus profonde, ou correspond au tout début de cette zone, la valeur mesurée par la sonde de débit sera nulle. Dans ce cas, on choisira une valeur ub très faible, égale par exemple à 1 % de la valeur U mesurée au sommet. Une étape préalable du processus itératif (bloc 71) consiste à déterminer les intervalles de valeurs autorisées pour IT Comme Te est présent au dénominateur dans l'équation ( 5), il faut interdire toutes les valeurs de IT qui annulent le 2538849, dénominateur Ces valeurs, pour des valeurs données des paramè- tres Tb et a, dépendent des valeurs prises par la température mesurée T Ces valeurs constituent un ensemble discret de valeurs qu'on peut assimiler à un intervalle continu On en déduit l'intervalle interdit pour IT (I Ta, I Tb) et ses complé- ments, à savoir les deux intervalles autorisés I I Ta) et I 2 (I Tb, +c Q). Ensuite, le processus itératif proprement dit peut se dérouler On sait que l'équation P(IT) U O = O a une solution unique I Tf, cette solution pouvant se trouver d'un côté ou de l'autre de l'intervalle interdit. Le processus itératif se déroule en fait en deux étapes, la première permettant d'identifier l'intervalle o se situe la solution et la seconde permettant d'arriver à la convergence. La première étape prend fin lorsqu'on a détecté le premier changement de signe de la fonction (IT) UO par rapport à un signe de référence Le signe de référence correspond à une valeur de l'itérant ITO située dans l'un des intervalles autorisés et très éloignée de la limite Par exemple on prend IT, = I Ta + K l I Ta I Tbl, K étant très grand. Les valeurs successives de l'itérant sont alors définies comme suit (bloc 73) On effectue une itération avec une valeur IT 1 située dans l'autre intervalle 12, également éloignée de la limite I Tb, et on poursuit en changeant chaque fois d'intervalle avec des valeurs IT 2, etc plus proches des limites I Ta et I Tb, par exemple IT 2 = I Ta + K/n l I Ta I Tbl, avec N 10 A chaque itération, le signe de la valeur F(I Ti) Uo est comparé ausigne de référence (bloc 74) La première valeur d'itérant donnant un changement de signe appartient à l'intervalle oh se situe la solution. 2538849. Dans la seconde étape, l'intégration se fait de la même façon (bloc 75) mais l'itérant est défini (bloc 76) selon l'un quelconque des procédés classiques applicables à une fonction convexe sur un intervalle l'itération est poursuivie jusqu'à la convergence, c'est-à-dire jusqu'à ce que la différence F(IT) U O devienne inférieure à une quantité déterminée E(bloc de comparaison 77). L'ensemble des valeurs de débit u(z) correspondant à la convergence sont enregistrées (bloc 78) en vue de l'obtention d'un profil de débit synthétique (bloc 79). * Sur la figure 3 a, le profil de température est représenté en trait pointillé, et le profil de débit synthéti- que résultant du traitement exposé ci-dessus est en trait mixte On note que le profil synthétique donne une image bien plus précise dos entrées de fluide que le profil mesuré En particulier, les entrées dans le zone supérieure A 2 sont prises en compte sans retard. L'exemple de la figure 3 b est encore plus net le profil synthétique permet d'évaluer les contributions indivi- -20 duelles des trois zones B 1, B 2, B 3 comme étant respectivement de 70 %, 15 % et 15 % Une telle évaluation n'était pas possible à partir du profil mesuré. Il est à noter que le procédé décrit ci-dessus est capable de fournir une valeur du gradient géothermique, si l'itération porte sur ce gradient Pour cela, on modifiera l'estimation a avec la -valeur A a correspondant à la convergence. Une variante au procédé décrit plus haut consiste à diviser l'intervalle total de mesure en plusieurs zones pour l'exécution du processus itératif On prend alors comme valeurs de débit imposées les valeurs mesurées par la sonde de débit au bas et au sommet de chaque zone. Le traitement commence par la zone la plus profonde et continue dans le sens des profondeurs décroissantes Pour la zone Z 1 la plus profonde, le traitement est identique à ce qui a été décrit, à ceci près que la valeur imposée U est la valeur de débit mesurée au sommet de cette zone. 2538849, Pour la zone suivante Z 2, on doit tenir compte du fait que la température mesurée au bas de la zone n'est plus du tout une estimation appropriée de la température T du fluide dans la formation On prend alors, comme estimation Tb pour la zone Z 2, la valeur de Te au sommet de la zone Z 1, cette valeur étant calculée avec la valeur finale I Tf de l'itérant. le même traitement s'applique successivement pour les zones suivantes. Dans le cas d'un puits producteur de gaz, il faut appliquer la correction liée à l'effet Joule-Thompson Le coefficient de Joule-Thompson Y et la différence de pression AP sont des paramètres que l'on peut considérer comme constants sur toute la hauteur de l'intervalle de mesure, et qu'il n'est pas difficile de déterminer. Cette correction est également à appliquer lorsqu'une entrée de gaz se produit dans un puits producteur de pétrole. les zones d'entrée de gaz peuvent être identifiées facilement par le profil de densité fourni par le gradiomanomètre En effet, l'entrée de gaz dans le puits se traduit par une forte chute de la densité. 2538849, REVENDICATIONS 1 Procédé pour déterminer le débit de fluide dans un puits, le long d'un intervalle comprenant des zones d'échange de fluide entre le puits et les formations, caractérisé par le fait que l'on mesure au moyen d'une sonde de diagraphie la température du fluide le long de cet intervalle, on mesure le débit au moins aux extrémités de cet intervalle, et on détermi- ne pour chaque profondeur, à partir des températures mesurées et des débits mesurés, une valeur de débit, l'ensemble des valeurs obtenues étant enregistré sous la -forme d'un profil de débit synthétique. 2 Procédé selon la revendication 1, caractérisé par le fait que le profil de débit synthétique est obtenu par inté- gration sur le dit intervalle d'une équation différentielle reliant, pour chaque profondeur, le débit à la température du fluide et à la température Te du fluide dans la formation, la résolution étant faite par itération, avec la température Te en tant que grandeur itérante, le processus d'itération reposant sur la comparaison de la valeur résultant de l'intégration sur l'ensemble de l'intervalle avec la valeur de débit mesurée au sommet de l'intervalle. 3 Procédé selon la revendication 2, caractérisé par le fait que, si le débit mesuré est nul au bas de l'intervalle, la valeur de départ pour l'intégration est choisie comme une petite fraction de la valeur de débit mesurée au sommet de l'intervalle. 4 Procédé selon l'une des revendications 2 et 3, caractérisé par le fait que, pour un liquide, l'équation diffé- rentielle s'écrit sous la forme (d T ld T \-ai/ tdz/e du__ _ _ _ _ dz T Te dz u étant le débit, T la température du fluide mesurée, Te la température du fluide dans la formation et (d T/dz)O le gradient de température dû au transfert thermique entre le puits et la formation, ce gradient étant déterminé à partir des mesures de température dans les zones du dit intervalle o aucun échange notable de fluide ne se produit.