EP0125164B1 - Méthode de détermination des caractéristiques d'une formation souterraine produisant un fluide - Google Patents

Méthode de détermination des caractéristiques d'une formation souterraine produisant un fluide Download PDF

Info

Publication number
EP0125164B1
EP0125164B1 EP84400781A EP84400781A EP0125164B1 EP 0125164 B1 EP0125164 B1 EP 0125164B1 EP 84400781 A EP84400781 A EP 84400781A EP 84400781 A EP84400781 A EP 84400781A EP 0125164 B1 EP0125164 B1 EP 0125164B1
Authority
EP
European Patent Office
Prior art keywords
theoretical
experimental
graph
function
evolution
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired
Application number
EP84400781A
Other languages
German (de)
English (en)
Other versions
EP0125164A1 (fr
Inventor
Dominique Bourdet
Timothy Whittle
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Flopetrol Services Inc
Original Assignee
Flopetrol Services Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Flopetrol Services Inc filed Critical Flopetrol Services Inc
Publication of EP0125164A1 publication Critical patent/EP0125164A1/fr
Application granted granted Critical
Publication of EP0125164B1 publication Critical patent/EP0125164B1/fr
Expired legal-status Critical Current

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • E21B49/008—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells by injection test; by analysing pressure variations in an injection or production test, e.g. for estimating the skin factor

Definitions

  • the present invention relates to the testing of oil wells, making it possible to determine the physical characteristics of the system formed by a well and an underground formation (also called “reservoir”) producing hydrocarbons through the well. More specifically, the invention relates to a method according to which the flow of fluid produced by the well is modified by closing or opening a valve which is located on the surface or in the well. The resulting pressure variations are measured and recorded at the bottom of the well as a function of the time elapsed since the start of the tests, that is to say since the modification of the flow rate. The characteristics of the underground well-formation system can be deduced from these experimental data.
  • the experimental data from the well tests are analyzed by comparing the response of the underground formation to a change in the flow rate of the fluid produced, with the behavior of theoretical models having well-defined characteristics and subjected to the same change in flow rate as the formation studied.
  • pressure variations as a function of time characterize the behavior of the well-formation system and the removal at constant flow rate of fluids, by opening a valve in the initially closed well. is the test condition that is applied to the training and the theoretical model.
  • their behaviors are identical, it is assumed that the system studied and the theoretical model are identical both quantitatively and qualitatively. In other words, these tanks are assumed to have the same physical characteristics.
  • the characteristics obtained from this comparison depend on the theoretical model: the more complicated the model, the more numerous the characteristics that can be determined.
  • the basic model is represented by a homogeneous formation with impermeable upper and lower limits and with an infinite radial extension. The flow in the formation is then radial, directed towards the well.
  • the wall effect is defined by a coefficient S which characterizes the damage or stimulation of the part of the formation adjacent to the well.
  • the compression or decompression effect of the fluid in the well is characterized by a coefficient C which results from the difference in the flow of fluid produced by the well, between the underground formation and the wellhead, when a valve located at the head well is either closed or open.
  • the coefficient C is usually expressed in barrels per psi, one barrel being equal to 0.16 m 3 and one psi equal to 0.069 bar.
  • each curve is characterized by one or more dimensionless numbers which each represent a characteristic (or a combination of characteristics) of the theoretical system formed by a well and a reservoir.
  • a dimensionless parameter is defined by the actual parameter (pressure for example) multiplied by an expression that includes certain characteristics of the well-reservoir system so as to make the dimensionless parameter independent of these characteristics.
  • the coefficient S characterizes only the parietal effect but is independent of the other characteristics of the reservoir and of the experimental conditions such as the flow rate, the viscosity of the fluid, the permeability of the formation. etc.
  • the experimental curve and one of the standard curves represented with the same coordinate scales have an identical shape but are offset with respect to each other.
  • the offsets along the two axes, on the ordinate for pressure and on the abscissa for time, are proportional to the characteristic values of the well-reservoir system which can thus be determined.
  • Qualitative information on the underground formation is obtained by the identification of the different regimes on the graph in logarithmic scale representing the experimental data. Knowing that a particular characteristic of the well-reservoir system, such as for example a vertical fracture, is characterized by a particular regime, all the different regimes appearing on the graph of the experimental data are identified in order to select the appropriate model of well-reservoir system. . Specialized graphs taking into account only part of the experimental data allow the more precise determination of the characteristics of the system. We then return to the graph on a logarithmic scale taking into account all the data to confirm the choice of the system and the quantitative determination of the characteristics of the training. The latter are obtained by selecting a standard curve having the same shape as the experimental curve and by determining the offset of the coordinate axes of the experimental curve with respect to the theoretical curve.
  • Patent 4,328,705 also describes a method in which the standard curves are represented using dimensionless pressure P o for the ordinate axis and the ratio td / Cd for the abscissa axis, t D being the dimensionless time and C d the dimensionless coefficient characterizing the compression or decompression effect of the fluid in the well.
  • P o for the ordinate axis
  • t D being the dimensionless time
  • C d the dimensionless coefficient characterizing the compression or decompression effect of the fluid in the well.
  • the present invention relates to a method for determining the characteristics of a well-reservoir system allowing better identification between the experimental behavior of the analyzed system formed by the well and the underground formation and the behavior of a theoretical model.
  • This model is general, namely that the formation can be homogeneous or heterogeneous and that it takes into account the parietal effect and the effect of compression or decompression of the fluid and possibly the double porosity of the reservoir and fractures of the well.
  • the method according to the present invention allows a global and unique analysis of the behavior of the well-reservoir system, without resorting to specialized analyzes.
  • the invention also allows the analysis of experimental data when the condition imposed on the system is the closing of the well, thanks to a judicious choice of parameters.
  • the method according to the present invention can also be advantageously combined with a method of the prior art.
  • Said theoretical evolution can advantageously be that of the logarithm of the product P ' D ⁇ t D / C D as a function of the logarithm of t D / C D and said experimental evolution is that of the logarithm of the product ⁇ P' ⁇ ⁇ t as a function of the logarithm of ⁇ t .
  • Said theoretical evolution can also be a function of an index representing a quantity characteristic of the product C D e 2s .
  • the change in fluid flow corresponds to the closing of the well, we can advantageously compare said theoretical evolution with the experimental evolution of the logarithm of the expression: as a function of the logarithm of the time intervals ⁇ t, tp being the time during which the well was put into production.
  • Certain steps of the present invention in particular the identification of experimental data with the behavior of a theoretical model having very precise characteristics, can be implemented using a computer.
  • these steps are advantageously implemented by drawing a theoretical graph in Cartesian coordinates and on a logarithmic scale, said graph representing the theoretical evolution of the derivative P ' o as a function of t D / C D or even the theoretical evolution of the product P ' D - t D / C D as a function of t D / C D.
  • the experimental curve correspond with one of the standard curves of the theoretical graph and we can determine certain physical characteristics of the underground well-formation system.
  • the subject of the invention is also the theoretical graphs obtained as indicated above.
  • measurements are generally made in order to determine the physical characteristics of the underground formation producing these hydrocarbons. This preliminary stage before production is very important because it makes it possible to define the most appropriate conditions to produce these hydrocarbons and possibly to improve production.
  • One of these measures consists in varying the flow rate of the fluid produced, by opening or closing a valve placed at the well head or in the well itself. and recording the resulting variations in pressure as a function of the time elapsed from the change in flow rate of the fluid produced. We can for example completely close the well and record the resulting pressure increase (we will then obtain an experimental curve called the “build-up” curve). We can also have a well produced whose production had been stopped and record the corresponding pressure reduction (the experimental curve obtained is then called "drawdown curve").
  • Pressure variations as a function of time can be monitored using a probe lowered into the well at the end of a cable.
  • the cable can be electric and in this case transmit the pressure information directly to a recorder located on the surface.
  • the pressure variations are recorded in memories located in the probe. These memories are then read on the surface.
  • One can also install a pressure gauge in a side pocket of the production column of the well, near the producing formation.
  • a conductive cable located in the annular space between the production column and the casing connects the pressure gauge to a recorder located on the surface.
  • the values measured by the pressure probes generally do not correspond to the pressure itself, but to a quantity characteristic of the pressure such as, for example, a difference of two frequencies. Thereafter, we will use the expression "value of pressure for convenience and clarity, while keeping in mind that the experimental data can correspond to a quantity pressure characteristic.
  • the graph in Figure 1 characterizes the behavior of a homogeneous reservoir model and a well exhibiting the wall effect and the effect of compression or decompression of the fluid in the well.
  • the curves of FIG. 1 are characterized by three distinct parts: the left part of the graph corresponding to short times and being characteristic of the effect of decompression of the fluid from the well (this effect is most important at the opening of the valve) ; the right part of the graph corresponding to a pure radial flow from the reservoir and an intermediate part between the left and right parts corresponding to a transient flow regime between the two preceding limit flows.
  • This intermediate flow is a function of the decompression effect of the fluid and of the wall effect.
  • the method of determining physical characteristics using the graph in Figure 1 has been improved by following the evolution, no longer of the mathematical derivative of the dimensionless pressure, but by following the evolution, as a function of t D / C D , of the product of the derivative P ' o of the dimensionless pressure (derivative with respect to t D / C D ) by the ratio t D / C D.
  • This new method is illustrated in FIG. 2 by a graph representing the behavior of a homogeneous formation presenting the parietal effect and the effect of decompression of the fluid from the well.
  • the data axis corresponds to P ' D.
  • t D / C D and the abscissa axis corresponds to t D / C D
  • P ' D being the derivative of P D with respect to t D / C D.
  • Equation (7) can be written:
  • FIG. 3 illustrates the use of the graph of the standard curves of FIG. 2. This graph has been reproduced in FIG. 3 with. on the ordinate, P ' D.
  • the value of Y makes it possible to determine the product kh.
  • the value of the fluid flow q is generally known by measurements having been carried out beforehand with a flow meter or a separator, and the values of the coefficient B of expansion of the fluid and of its viscosity ⁇ are determined by the analysis of fluid samples (usually called “PVT”). Consequently, the value of the permeability-thickness product kh can be determined by knowing the value Y measured.
  • the value of the coefficient S of the wall effect is determined by the correspondence of the experimental curve with one of the standard curves, the correspondence of the two curves leading to the value of C D e 2s .
  • the value of C o is determined by the value of C by the following equation: in which ⁇ bc, h represents the porosity-compressibility-thickness product. known from geological studies (such as, for example, analysis of samples or electrical logs) and r is the radius of the well.
  • the value of the coefficient S can therefore be calculated from the value of C D e 2s ,
  • the standard curves represented in FIGS. 1 and 2 correspond to the behavior of a theoretical model of homogeneous formation when suddenly an increase in the flow rate of the fluid is imposed. produced by the formation, and more especially when a valve is opened on the surface of the well to make it produce at constant flow when it was previously closed (“drawdown • curve).
  • the experimental curve is plotted on a logarithmic scale by plotting the time intervals ⁇ t on the x axis and on the y axis: tp representing the duration during which the training was put into production.
  • tp representing the duration during which the training was put into production.
  • the analysis of well tests can then be done by comparing this experimental curve with the standard curves of the graph in Figure 2.
  • FIG. 4 shows an example of application to a formation having a double porosity.
  • the fluid produced by the formation is contained in the matrix, that is to say in the rock composing the formation. and in the interstices or cracks contained in the matrix. We therefore have a system in which the fluid contained in the matrix first flows through the cracks before passing into the well.
  • the coefficient ⁇ characterizes the ratio of the volume of fluid produced by the cracks to the volume of fluid produced by the total system (matrix + crack).
  • the coefficient ⁇ characterizes the delay of the matrix to produce the fluid in the cracks compared to the production of the cracks themselves.
  • the graph in Figure 4 corresponds to a theoretical training model with double porosity. On this graph, we have represented, in solid lines, the standard curves corresponding to the homogeneous model, identical to those of FIG. 2, in dotted lines of the standard curves by choosing as index and semi-dotted standard curves by choosing as index
  • a typical experimental curve characterizing a double porosity formation has also been represented by points.
  • the use of the graph in FIG. 4 makes it possible to determine the values of the coefficients w and ⁇ , in addition to the values of kh, C and S. It is noted that the curves characterizing the behavior of a heterogeneous model have a very marked shape in applying the method according to the present invention.
  • the present invention also makes it possible to plot on the same theoretical graph the standard curves of FIG. 2.
  • P ' D. t D / C D as a function of t D / C D but also the standard curves P o as a function of t D / C D described in US Pat. No. 4,328,705.
  • the juxtaposition of these two series of standard curves on the same graph is shown in Figure 5.
  • the method for determining the characteristics of an underground formation which has just been described has many advantages.
  • the analysis of well tests can be carried out using a single graph, whereas the methods of the prior art call for a general graph on a logarithmic scale using all the experimental data and a specialized graph in semi-logarithmic scale and taking into account only part of the experimental data.
  • the calibration of the experimental curve with the standard curves can be done without ambiguity.
  • the combination of the standard curves of the prior art with the standard curves of the present invention on the same graph has a certain advantage.
  • defining a new time. given by equation (15) makes it possible to analyze the well tests carried out by closing the well.

Landscapes

  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Fluid Mechanics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Sampling And Sample Adjustment (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Description

  • La présente invention concerne les essais de puits d'hydrocarbures, permettant de déterminer les caractéristiques physiques du système formé d'un puits et d'une formation souterraine (appelée aussi « réservoir ») produisant des hydrocarbures à travers le puits. De façon plus précise, l'invention se rapporte à une méthode selon laquelle le débit de fluide produit par le puits est modifié en fermant ou en ouvrant une vanne qui se trouve à la surface ou dans le puits. Les variations de pression résultantes sont mesurées et enregistrées au fond du puits en fonction du temps écoulé depuis le début des essais, c'est-à-dire depuis la modification du débit. Les caractéristiques du système puits-formation souterraine peuvent être déduites de ces données expérimentales. Les données expérimentales des essais de puits sont analysées en comparant la réponse de la formation souterraine à un changement de débit du fluide produit, avec le comportement de modèles théoriques ayant des caractéristiques bien définies et soumis au même changement de débit que la formation étudiée. Habituellement, les variations de pression en fonction du temps caractérisent le comportement du système puits-formation et l'enlèvement à débit constant de fluides, par l'ouverture d'une vanne dans le puits initialement fermé. est la condition d'essai qui est appliquée à la formation et au modèle théorique. Lorsque leurs comportements sont identiques, on suppose que le système étudié et le modèle théorique sont identiques aussi bien au point de vue quantitatif que qualitatif. En d'autres termes, ces réservoirs sont supposés avoir les mêmes caractéristiques physiques.
  • Les caractéristiques obtenues de cette comparaison dépendent du modèle théorique : plus le modèle est compliqué, plus les caractéristiques qui peuvent être déterminées sont nombreuses. Le modèle de base est représenté par une formation homogène avec des limites supérieure et inférieure imperméables et avec une extension radiale infinie. Le débit dans la formation est alors radial, dirigé vers le puits.
  • Cependant, le modèle théorique le plus couramment utilisé est plus compliqué. Il comporte les caractéristiques du modèle de base auquel on ajoute des conditions internes telles que l'effet pariétal («skin effect en anglais) et l'effet de compression ou de décompression du fluide dans le puits (« wellbore storage » en anglais). L'effet pariétal est défini par un coefficient S qui caractérise l'endommagement ou la stimulation de la partie de la formation adjacente au puits. L'effet de compression ou de décompression du fluide dans le puits est caractérisé par un coefficient C qui résulte de la différence de débit de fluide produit par le puits, entre la formation souterraine et la tête de puits, lorsqu'une vanne située en tête de puits est soit fermée soit ouverte. le coefficient C est exprimé habituellement en baril par psi, un baril étant égal à 0,16 m3 et un psi égal à 0,069 bar.
  • Le comportement d'un modèle théorique est représenté de façon commode par un graphe de courbes-type qui représentent les variations de pression, en fonction du temps, du fluide au fond du puits. Ces courbes sont habituellement tracées en coordonnées cartésiennes et en échelle logarithmique, la pression sans dimension étant portée en ordonnée et le temps sans dimension étant porté en abscisse. De plus, chaque courbe est caractérisée par un ou plusieurs nombres sans dimension qui représentent chacun une caractéristique (ou une combinaison de caractéristiques) du système théorique formé par un puits et un réservoir. Un paramètre sans dimension est défini par le paramètre réel (la pression par exemple) multiplié par une expression qui inclut certaines caractéristiques du système puits-réservoir de façon à rendre le paramètre sans dimension indépendant de ces caractéristiques.
  • C'est ainsi que le coefficient S caractérise uniquement l'effet pariétal mais est indépendant des autres caractéristiques du réservoir et des conditions expérimentales telles que le débit, la viscosité du fluide, la perméabilité de la formation. etc. Lorsque le modèle théorique et le système étudié puits-formation correspondent, la courbe expérimentale et l'une des courbe-types représentées avec les mêmes échelles de coordonnées ont une forme identique mais sont décalées l'une par rapport à l'autre. Les décalages suivant les deux axes, en ordonnée pour la pression et en abscisse pour le temps, sont proportionnels à des valeurs de caractéristiques du système puits-réservoir qui peuvent ainsi être déterminées.
  • Des informations qualitatives sur la formation souterraine, telle que par exemple la présence d'une fracture, sont obtenues par l'identification des différents régimes sur le graphe en échelle logarithmique représentant les données expérimentales. Sachant qu'une caractéristique particulière du système puits-réservoir, telle que par exemple une fracture verticale, se caractérise par un régime particulier, tous les différents régimes apparaissant sur le graphe des données expérimentales sont identifiées pour sélectionner le modèle de système puits-réservoir approprié. Des graphes spécialisés ne tenant compte que d'une partie des données expérimentales permettent la détermination plus précise des caractéristiques du système. On revient ensuite au graphe en échelle logarithmique prenant en compte toutes les données pour confirmer le choix du système et la détermination quantitative des caractéristiques de la formation. Ces dernières sont obtenues en sélectionnant une courbe-type ayant la même forme que la courbe expérimentale et en déterminant le décalage des axes de coordonnées de la courbe expérimentale par rapport à la courbe théorique.
  • A un même modèle théorique correspondent plusieurs graphes de courbes-type. Ceci dépend des paramètres sans dimension choisis pour la représentation des axes de coordonnées du graphe, ainsi que d'un ou plusieurs index», Un « index n'est autre qu'un paramètre supplémentaire (ou une combinaison de paramètres) choisi pour la représentation des courbes, en complément des paramètres sans dimension des axes de coordonnées. La comparaison des différentes méthodes utilisées est donnée dans l'article intitulé « A Comparison Between Different Skin and Wellbore Storage Type Curves for Early-Time Transient Analysis » (Comparaison entre les différentes courbes-type avec effet pariétal et effet de compression ou de décompression pour l'analyse des transitoires précoces) par A.C. Gringarten & al., publié par « Society of Petroleum Engineers of AIME », (n° SPE 8205). Le brevet des Etats-Unis d'Amérique 4,328,705 décrit également une méthode selon laquelle les courbes-type sont représentées en utilisant la pression sans dimension Po pour l'axe des ordonnées et le rapport td/Cd pour l'axe des abscisses, tD étant le temps sans dimension et Cd le coefficient sans dimension caractérisant l'effet de compression ou de décompression du fluide dans le puits. L'inconvénient de la méthode décrite dans ce brevet est que les courbes-type ont des formes variant relativement lentement l'une par rapport à l'autre. Il en résulte une certaine incertitude dans le choix de la courbe-type correspondant à la courbe expérimentale. On remarque également que pour effectuer une analyse complète, on est obligé de faire appel non seulement à un graphe en échelle logarithmique représentant l'ensemble des données expérimentales, mais d'utiliser également des graphes spécialisés, en échelle semi-logarithmique par exemple, pour n'analyser qu'une partie des données mais d'une façon plus précise.
  • On a déjà songé à utiliser la dérivée mathématique de la pression sans dimension, P'D, au lieu de la pression sans dimension Pd. C'est ainsi que dans l'article intitulé « Application of the P'p Function to Interference Analysis» (Application de la fonction P'D à l'analyse des tests d'interférence) publié dans « Journal of Petroleum Technology " d'août 1980, page 1465, l'évolution de la dérivée P'D (dérivée par rapport à tD) en fonction de to est utilisée pour analyser des essais d'interférence entre un puits de production et un puits d'observation. On enregistre les variations de pression dans le puits d'observation lorsqu'on modifie le débit de fluide produit par le puits de production. Dans ce cas. l'effet pariétal et l'effet de compression ou de décompression du fluide n'interviennent pas. C'est donc un cas très simple dans lequel on analyse la réponse de la formation souterraine dans un puits éloigné du puits producteur. Il en résulte qu'il n'y a pas une famille de courbes-type mais une seule courbe.
  • La dérivée de la' pression P'D (dérivée par rapport à tD) a été également mise à profit pour caractériser les réservoirs contenant deux failles étanches parallèles autour du réservoir, dans l'article intitulé « Detection and Location of Two Parallel Sealing Faults Around a Well » (Détection et localisation de deux failles étanches parallèles autour d'un puits) paru dans « Journal of Petroleum Technology d'octobre 1980, page 1701. Cet article ne traite que d'un problème particulier.
  • Le comportement en pression d'un puits produisant un fluide légèrement compressible à travers un plan unique d'une fracture verticale dans un réservoir infini a été analysé à l'aide de la dérivée mathématique de la pression sans dimension P'o, (dérivée par rapport à un temps sans dimension tDf) dans l'article intitulé « Application of P'D Function to Vertically Fractured Wells " (Application de la fonction P'o au puits fracturé verticalement) paru dans « Society of Petroleum Engineers of AIME, SPE 11028, du 26-29 septembre 1982.
  • Cet article ne concerne qu'un cas particulier dans lequel la courbe-type est unique et pour lequel les avantages de l'utilisation de la dérivée de la pression ne sont pas évidents par rapport aux méthodes classiques. De plus, l'effet pariétal et l'effet de compression ou de décompression du fluide n'interviennent pas.
  • La présente invention a pour objet une méthode de détermination des caractéristiques d'un système puits-réservoir permettant une meilleure identification entre le comportement expérimental du système analysé formé par le puits et la formation souterraine et le comportement d'un modèle théorique. Ce modèle est général, à savoir que la formation peut être homogène ou hétérogène et qu'il tient compte de l'effet pariétal et de l'effet de compression ou de décompression du fluide et éventuellement de la double porosité du réservoir et des fractures du puits. La méthode selon la présente invention permet une analyse globale et unique du comportement du système puits-réservoir, sans recourir à des analyses spécialisées. L'invention permet également l'analyse des données expérimentales lorsque la condition imposée au système est la fermeture du puits, grâce à un choix judicieux des paramètres. La méthode selon la présente invention peut également être avantageusement combinée à une méthode de l'art antérieur.
  • De façon plus précise, la présente invention concerne une méthode de détermination des caractéristiques physiques d'un système formé d'un puits et d'une formation souterraine contenant un fluide et communiquant avec ledit puits, ladite formation présentant de l'effet pariétal et/ou le fluide se comprimant ou se décomprimant dans le puits et ladite formation étant homogène ou hétérogène. Selon la méthode, on provoque un changement de débit dudit fluide, on mesure une grandeur caractéristique de la pression P du fluide à des intervalles de temps successifs Δt et on compare.
    • - d'une part, à partir d'un modèle théorique de système de puits-réservoir, l'évolution théorique du logarithme de la dérivée P'D de la pression sans dimension en fonction du logarithme de tD/CD, ladite dérivée P'D étant par rapport à tD/CD, to représentant le temps sans dimension et CD le coefficient sans dimension de l'effet de compression ou de décompression du fluide dans le puits, avec
    • - d'autre part, l'évolution expérimentale du logarithme de la dérivée ΔP' de la pression en fonction du logarithme des intervalles de temps correspondants Δt, ladite dérivée ΔP' étant par rapport au temps t. et on détermine de la comparaison desdites évolutions théorique et expérimentale au moins une caractéristique du système puits-formation, choisie parmi le produit kh de la perméabilité k par l'épaisseur de ladite formation h, le coefficient CD et le coefficient S de l'effet pariétal.
  • Ladite évolution théorique peut avantageusement être celle du logarithme du produit P'D ◆ tD/CD en fonction du logarithme de tD/CD et ladite évolution expérimentale est celle du logarithme du produit ΔP' · Δt en fonction du logarithme de Δt.
  • Ladite évolution théorique peut également être fonction d'un index représentant une grandeur caractéristique du produit CDe2s. Lorsque le changement de débit du fluide correspond à la fermeture du puits, on peut avantageusement comparer ladite évolution théorique avec l'évolution expérimentale du logarithme de l'expression :
    Figure imgb0001
    en fonction du logarithme des intervalles de temps △t, tp étant le temps pendant lequel le puits a été mis en production.
  • Certaines étapes de la présente invention, notamment l'identification des données expérimentales avec le comportement d'un modèle théorique ayant des caractéristiques bien précises, peuvent être mises en oeuvre à l'aide d'un ordinateur. Cependant, ces étapes sont avantageusement mises en oeuvre en traçant un graphe théorique en coordonnées cartésiennes et en échelle logarithmique, ledit graphe représentant l'évolution théorique de la dérivée P'o en fonction de tD/CD ou encore l'évolution théorique du produit P'D - tD/CD en fonction de tD/CD.
  • On peut également tracer une courbe expérimentale à l'aide des données expérimentales avec la même échelle logarithmique que ledit graphe théorique, la courbe expérimentale représentant soit l'évolution expérimentale de ΔP' en fonction de △t, soit l'évolution expérimentale du produit ΔP' · Δt en fonction de Δt. On peut alors faire correspondre la courbe expérimentale avec l'une des courbes-type du graphe théorique et on peut en déterminer certaines caractéristiques physiques du système puits-formation souterraine.
  • L'invention a également pour objet les graphes théoriques obtenus tel qu'indiqué précédemment.
  • L'invention sera mieux comprise à l'aide de la description qui va suivre de modes de réalisation de l'invention donnés à titre d'exemples explicatifs mais nullement limitatifs. La description se rapporte aux dessins qui l'accompagnent dans lesquels :
    • la figure 1 représente en échelle logarithmique un graphe de courbe-type représentant P'o en fonction de tD/CD, l'index représentant les valeurs de CDe2S ;
    • la figure 2 montre un graphe de courbes-type en échelle logarithmique représentant P'D . tD/CD en fonction de tD/CD, l'index étant CDe2S ;
    • la figure 3 illustre la méthode selon la présente invention pour la détermination des caractéristiques physiques d'une formation souterraine produisant un fluide ;
    • la figure 4 représente en échelle logarithmique un graphe de courbes-type représentant P' D . tD/CD en fonction de tD/CD pour une formation souterraine à double porosité ; et
    • la figure 5 représente deux séries de courbes-type en échelle logarithmique, l'une montrant des courbes-type connues et l'autre montrant des courbes-type selon la présente invention.
  • Avant de mettre un puits d'hydrocarbures en production, on effectue généralement des mesures de façon à déterminer les caractéristiques physiques de la formation souterraine produisant ces hydrocarbures. Cette étape préliminaire avant la production est très importante car elle permet de définir les conditions les plus appropriées pour produire ces hydrocarbures et éventuellement pour améliorer la production. L'une de ces mesures consiste à faire varier le débit du fluide produit, en ouvrant ou fermant une vanne placée en tête de puits ou dans le puits lui-même. et à enregistrer les variations résultantes de pression en fonction du temps écoulé à partir de la modification de débit du fluide produit. On peut par exemple fermer complètement le puits et enregistrer l'augmentation de pression résultante (on obtiendra alors une courbe expérimentale dite courbe de « build-up •). On peut aussi faire produire un puits dont la production avait été arrêtée et enregistrer la diminution de pression correspondante (la courbe expérimentale obtenue s'appelle alors « courbe de drawdown •).
  • Les variations de pression en fonction du temps peuvent être suivies à l'aide d'une sonde descendue dans le puits à l'extrémité d'un câble. Le câble peut être électrique et dans ce cas transmettre les informations de pression directement à un enregistreur situé à la surface. Lorsque le câble n'est pas conducteur, les variations de pression sont enregistrées dans des mémoires disposées dans la sonde. Ces mémoires sont ensuite lues en surface. On peut également installer une jauge de pression dans une poche latérale de la colonne de production du puits, à proximité de la formation productrice. Un câble conducteur situé dans l'espace annulaire compris entre la colonne de production et le tubage relie la jauge de pression jusqu'à un enregistreur situé à la surface. Un tel dispositif est décrit par exemple dans les brevets des Etats-Unis n° 3,939.705 et 4.105.279.
  • Les valeurs mesurées par les sondes de pression ne correspondent généralement pas à la pression elle-même, mais à une grandeur caractéristique de la pression telle que, par exemple, une différence de deux fréquences. Par la suite, on utilisera l'expression « valeur de la pression par commodité et clarté, tout en gardant à l'esprit que les données expérimentales peuvent correspondre à une grandeur caractéristique de la pression.
  • La figure 1 représente un graphe de nouvelles courbes-type, en échelle logarithmique, représentant la dérivée mathématique P'D de la pression sans dimension PD en fonction du rapport tD/CD, tD représente le temps sans dimension et CD représentant la valeur du coefficient sans dimension de l'effet de compression ou de décompression du fluide dans le puits. La dérivée mathématique P'D est faite par rapport à tD/CD. De plus, les variations de la dérivée de la pression P'D sont représentées par rapport à un index CDe2s, qui n'est autre qu'une combinaison de deux caractéristiques physiques CD et S du système puits-réservoir analysé. Il est à remarquer que l'index CDe2s peut prendre une valeur quelconque, qui n'est pas nécessairement entière, La valeur de la pression sans dimension PD est donnée par l'équation suivante, en employant le système d'unités utilisé couramment dans l'industrie pétrolière et appelé « Oil field units page 185 du livre intitulé « Advances in Well Test Analysis publié par « Society of Petroleum Engineers of AIME» ― 1977 :
    Figure imgb0002
    dans laquelle :
    • k représente la perméabilité de la formation souterraine,
    • h est l'épaisseur de la formation,
    • ΔP est la variation de pression.
    • 5q est le débit du fluide en surface,
    • B est le coefficient de dilatation du fluide entre le réservoir et la surface (appelé en anglais « formation volume factor ») et
    • µ est la viscosité du fluide.
  • La dérivée mathématique P'D de la pression sans dimension PD, par rapport à tD/CD, est donnée par l'équation suivante :
    Figure imgb0003
    dans laquelle ΔP' est la dérivée (par rapport au temps t) de la variation de pression ΔP en fonction de l'intervalle de temps Δt qui représente l'intervalle de temps s'étant écoulé depuis le début de l'essai de la formation, c'est-à-dire l'intervalle de temps entre l'instant de mesure et l'instant de modification de débit du fluide.
  • La valeur du rapport tD/CD, dans le même système d'unités que pour les équations précédentes, est donnée par:
    Figure imgb0004
    dans laquelle C est l'effet de compression ou de décompression du fluide dans le puits.
  • Le graphe de la figure 1 caractérise le comportement d'un modèle de réservoir homogène et d'un puits présentant l'effet pariétal et l'effet de compression ou de décompression du fluide dans le puits.
  • Ce graphe est obtenu à partir de l'équation (A.2) de l'article intitulé « Determination of fissure volume and block size in fractured reservoirs by type-curve analysis » (Détermination du volume des fissures et de la taille des blocs dans les réservoirs fissurés par analyse par courbes-type) publié par « Society of Petroleum Engineers » en septembre 1980. n° SPE 9293. Cette équation est donnée dans le domaine de Laplace. L'inversion dans le domaine des temps réels est obtenue à l'aide d'un algorithme d'inversion, tel que celui décrit par exemple par H. Stehfest dans « Communications of the ACM, D-5 du 13.01.70. n° 1 page 47.
  • Les courbes de la figure 1 se caractérisent par trois parties distinctes : la partie gauche du graphe correspondant aux temps courts et étant caractéristique de l'effet de décompression du fluide du puits (cet effet est le plus important à l'ouverture de la vanne) ; la partie droite du graphe correspondant à un écoulement pur radial du réservoir et une partie intermédiaire entre les parties gauche et droite correspondant à un régime d'écoulement transitoire entre les deux écoulements-limite précédents. Cet écoulement intermédiaire est fonction de l'effet de décompression du fluide et de l'effet pariétal.
  • Sur la partie gauche du graphe, les courbes tendent vers une asymptote correspondant à une dérivée égale à 1. En effet, au tout début des essais, le phénomène prédominant est l'effet de décompression du puits lequel est caractérisé par l'équation :
    Figure imgb0005
    • La dérivée de la pression sans dimension par rapport à tdCD peut s'écrire :
      Figure imgb0006
  • On voit que la dérivée P'D pour ce type d'écoulement est égal à 1 et que les courbes-type se réduisent à une droite de pente nulle. La partie droite du graphe de la figure 1. qui correspond à un écoulement radial infini dans une formation homogène, est caractérisée par l'équation :
    Figure imgb0007
    1 n représentant le logarithme néperien.
  • En dérivant PD par rapport à tD/CD, on obtient :
    Figure imgb0008
    et en passant en échelle logarithmique :
    Figure imgb0009
  • On remarque que la courbe représentée par l'équation (8) est une droite de pente égale à - 1. Pour les temps courts et les temps longs, les courbes sont rectilignes et indépendantes de CDe2s, ce qui est un avantage considérable par rapport aux méthodes de l'art antérieur, Entre les deux asymptotes, pour les temps intermédiaires, chaque courbe d'index CDe2s a une forme différente bien contrastée.
  • Si dP représente la différence de deux mesures successives de la pression du fluide dans le puits et si dt représente l'intervalle de temps (court) séparant ces deux mesures successives, on calcule pour tous les couples de mesures successives les valeurs ΔP' = dP/dt. Ce calcul permet de déterminer de façon pratique les valeurs successives de la dérivée mathématique ΔP' qui par définition est égale au rapport dP/dt lorsque dt tend vers zéro. En traçant la courbe ΔP' en fonction de Δt (Δt étant l'intervalle de temps entre l'instant de la mesure considérée et l'instant de modification de débit du fluide) de façon à former un graphe expérimental, et en prenant les mêmes échelles logarithmiques que celles adoptées pour tracer les courbes-type de la figure 1, on peut déterminer les caractéristiques physiques du système puits-formation souterraine. En effet, le décalage des ordonnées de la courbe expérimentale et des courbes-type permet de déterminer la valeur de C (ce qui est évident au vu de l'équation (2), en effectuant log P'D - log ΔP' et en connaissant les valeurs de q et B). Le décalage des abscisses de la courbe expérimentale par rapport à la courbe-type choisie permet de déterminer la valeur kh (connaissant C et µ, ce qui est évident au vu de l'équation (3) en effectuant log tD/CD - log Δt). Enfin le choix de la courbe-type correspondant à la courbe expérimentale permet de déterminer le coefficient S (par le calcul préalable de Co à partir de l'équation (14) comme il sera montré par la suite). Le graphe théorique de la figure 1 étant utilisé de la même façon que celui de la figure 2, par comparaison avec la courbe expérimentale, seule l'utilisation du graphe de la figure 2 est illustrée (figure 3).
  • . La méthode de détermination des caractéristiques physiques par utilisation du graphe de la figure 1 a été améliorée en suivant l'évolution, non plus de la dérivée mathématique de la pression sans dimension, mais en suivant l'évolution, en fonction de tD/CD, du produit de la dérivée P'o de la pression sans dimension (dérivée par rapport à tD/CD) par le rapport tD/CD. Cette nouvelle méthode est illustrée à la figure 2 par un graphe représentant le comportement d'une formation homogène présentant l'effet pariétal et l'effet de décompression du fluide du puits. L'axe des données correspond à P'D . tD/CD et l'axe des abscisses correspond à tD/CD, P'D étant la dérivée de PD par rapport à tD/CD.
  • De plus, l'index CDe2s a été choisi pour représenter les courbes-type. Comme dans le cas de la figure 1, l'effet prédominant au début de l'essai de puits est l'effet de-décompression du fluide. Cet effet correspond aux équations (4) et (5). A partir de l'équation (5), on peut écrire :
    Figure imgb0010
  • On remarque dans cette dernière équation, que pour les temps courts, les courbes-type tendent vers une asymptote de pente égale à 1.
  • Pour les temps longs, correspondant à la partie droite du graphe de la figure 2, les équations (6) et (7) demeurent valables puisqu'on se trouve en fin de test en écoulement radial infini pour une formation homogène. L'équation (7) peut s'écrire :
    Figure imgb0011
  • Il en résulte que pour les temps longs, la valeur du produit P'D . tD/CD est égale à 0.5 et que les courbes-type tendent vers une asymptote de pente nulle.
  • On remarque que pour le régime d'écoulement intermédiaire situé au centre du graphe de la figure 2, les courbes-type sont de forme très contrastée, ce qui permet une identification beaucoup plus précise de la courbe expérimentale avec l'une des courbes-type que par les méthodes de l'art antérieur. Par rapport au graphe de la figure 1, on peut dire que le graphe de la figure 2 correspond, en première approximation. à une rotation de 45° du graphe de la figure 1. Cependant, les courbes-type ont un relief plus accentué et la présentation du graphe de la figure 2 est plus pratique. Les valeurs de l'index CDe2s sont indiquées sur les courbes-type. La figure 3 illustre l'utilisation du graphe des courbes-type de la figure 2. Ce graphe a été reproduit sur la figure 3 avec. en ordonnée, P'D . tD/CD et en abscisse tD/CD. Les différences de pression dP mesurées dans le puits, pour des différences de temps successifs dt, sont utilisées pour calculer les valeurs ΔP' = dP/dt comme indiqué précédemment. On multiplie les valeurs successives de ΔP' par les intervalles de temps correspondants Δt et on trace ensuite un graphe expérimental représentant le produit ΔP' . At en ordonnée en fonction de Δt en abscisse. Les valeurs de ΔP sont en psi (1 psi = 0,068 bar) et les valeurs de At sont en heure. Les deux graphes théorique et expérimental ont la même échelle logarithmique. On commence par superposer la partie droite, qui est rectiligne. de la courbe expérimentale tracée sur la figure 3 à l'aide de points, à la partie rectiligne des courbes-type à droite du graphe. Ceci est facile à réaliser puisque cette partie des courbes est une droite de pente nulle. On décale ensuite le graphe expérimental le long de l'axe des temps de façon à faire correspondre sa partie gauche avec la partie gauche des courbes-type. Ceci est également facile puisque cette partie des courbes-type est une droite de pente égale à 1. Si la formation souterraine étudiée a un comportement homogène, la courbe expérimentale doit alors se superposer parfaitement, aux erreurs de précision des mesures près, à une courbe-type. Dans l'exemple montré sur la figure 3, cette courbe-type correspond à CDe2S = 1010. Le décalage des axes de coordonnées de la courbe expérimentale avec les axes des courbes-type permet de déterminer les valeurs du produit kh et la valeur de l'effet de décompression du fluide dans le puits. En effet. en combinant les équations (2) et (3). on obtient :
    Figure imgb0012
    qui s'écrit :
    Figure imgb0013
  • Le membre de gauche de cette dernière équation correspond au décalage des ordonnées représenté par Y sur la figure 3.
  • La valeur de Y permet de déterminer le produit kh. En effet, la valeur du débit de fluide q est généralement connue par des mesures ayant été préalablement effectuées avec un débitmètre ou un séparateur, et les valeurs du coefficient B de dilatation du fluide et de sa viscosité µ sont déterminées par l'analyse d'échantillons du fluide (analyse appelée habituellement « PVT •). En conséquence, la valeur du produit perméabilité-épaisseur kh peut être déterminée en connaissant la valeur Y mesurée.
  • De la même façon, l'équation (3) peut s'écrire
  • Figure imgb0014
  • Le membre de gauche de cette équation correspond au décalage X des abscisses de la courbe-type choisi et de la courbe expérimentale. Connaissant la valeur de ce décalage X ainsi que les valeurs de la viscosité µ et du produit kh, on déduit de l'équation (13) la valeur du coefficient C de l'effet de décompression du fluide du puits.
  • La valeur du coefficient S de l'effet pariétal est déterminé par la correspondance de la courbe expérimentale avec l'une des courbes-type, la correspondance des deux courbes conduisant à la valeur de CDe2s. La valeur de Co est déterminée par la valeur de C par l'équation suivante :
    Figure imgb0015
    dans laquelle <b c,h représente le produit porosité-compressibilité-épaisseur. connu à partir d'études géologiques (telles que par exemple, l'analyse d'échantillons ou les diagraphies électriques) et r est le rayon du puits. La valeur du coefficient S peut donc être calculée à partir de la valeur de CDe2s,
  • Les courbes-type représentées sur les figures 1 et 2 correspondent au comportement d'un modèle théorique de formation homogène lorsqu'on impose soudainement une augmentation du débit du fluide produit par la formation, et plus spécialement lorsqu'on ouvre une vanne en surface du puits pour le faire produire à débit constant alors qu'il était fermé auparavant (courbe de « drawdown •).
  • Selon l'une des caractéristiques de la présente invention, pour les analyses d'essais de puits correspondant à la fermeture du puits, la courbe expérimentale est tracée en échelle logarithmique en portant en abscisse les intervalles de temps Δt et en ordonnée :
    Figure imgb0016
    tp représentant la durée pendant laquelle la formation a été mise en production. L'analyse des essais de puits peut alors être faite par la comparaison de cette courbe expérimentale avec les courbes-type du graphe de la figure 2.
  • La représentation des courbes-type, avec en ordonnée P'D . tD/CD et en abscisse tD/CD, est utilisable non seulement pour les formations souterraines homogènes mais également pour les formations non homogènes présentant par exemple une double porosité. La figure 4 montre un exemple d'application à une formation ayant une double porosité. Dans ce cas, le fluide produit par la formation est contenu dans la matrice, c'est-à-dire dans la roche composant la formation. et dans les interstices ou fissures contenues dans la matrice. On a donc un système dans lequel le fluide contenu dans la matrice s'écoule d'abord dans les fissures avant de passer dans le puits. Le coefficient ω caractérise le rapport du volume de fluide produit par les fissures au volume de fluide produit par le système total (matrice + fissure). Le coefficient λ caractérise le retard de la matrice à produire le fluide dans les fissures par rapport à la production des fissures elles-mêmes. Le graphe de la figure 4 correspond à un modèle théorique de formation ayant une double porosité. Sur ce graphe, on a représenté, en traits pleins, les courbes-type correspondant au modèle homogène, identiques à celles de la figure 2, en pointillés des courbes-type en choisissant comme index
    Figure imgb0017
    et en semi-pointillé des courbes-type en choisissant comme index
    Figure imgb0018
  • Les courbes en pointillés représentent l'équation :
    Figure imgb0019
  • Les courbes en semi-pointillés représentent l'équation :
    Figure imgb0020
  • On a également représenté par des points une courbe expérimentale typique caractérisant une formation à double porosité. L'utilisation du graphe de la figure 4 permet de déterminer les valeurs des coefficients w et λ, en plus des valeurs de kh, C et S. On remarque que les courbes caractérisant le comportement d'un modèle hétérogène ont une forme très marquée en appliquant la méthode selon la présente invention.
  • La présente invention permet également de tracer sur un même graphe théorique les courbes-type de la figure 2. P'D . tD/CD en fonction de tD/CD mais également les courbes-type Po en fonction de tD/CD décrites dans le brevet des Etats-Unis d'Amérique n° 4.328,705. La juxtaposition de ces deux séries de courbes-type sur un même graphe est montrée sur la figure 5. On peut en effet réaliser cette superposition sur un même graphe car pour passer de P'D . tD/CD aux données expérimentales qui sont ΔP' · Δt, il faut multiplier ces dernières par un coefficient qui est donné par l'équation (11 ). Pour passer de PD aux données expérimentales ΔP, dans le cas des courbes-type du brevet cité, il faut multiplier ces dernières par le même coefficient que précédemment. On peut donc superposer les deux séries de courbes-type et porter en ordonnée, avec une même échelle, PD et P'D . tD/CD. Pour utiliser le graphe théorique de la figure 5. on utilise alors un même graphe expérimental comportant deux courbes représentant en ordonnée les variations de pression ΔP pour l'une et ΔP'D . tD/CD pour l'autre, Δt étant porté en abscisse pour les deux courbes. Le graphe combiné de la figure 5 permet une comparaison plus précise des deux courbes expérimentales avec les courbes-type.
  • La méthode de détermination des caractéristiques d'une formation souterraine qui vient d'être décrite présente de nombreux avantages. Ainsi, l'analyse des essais de puits peut-être effectuée à l'aide d'un seul graphe, alors que les méthodes de l'art antérieur font appel à un graphe général en échelle logarithmique en utilisant toutes les données expérimentales et à un graphe spécialisé en échelle semi-logarithmique et ne prenant en compte qu'une partie des données expérimentales. Du fait du comportement des modèles de systèmes formation-puits au début et en fin d'essai de puits (temps courts et temps longs sur les graphes) qui se traduit par des lignes droites de pentes bien définies pour les deux extrémités des courbes-type, le calage de la courbe expérimentale avec les courbes-type peut se faire sans ambiguïté. La combinaison des courbes-type de l'art antérieur avec les courbes-type de la présente invention sur un même graphe présente un avantage certain. De plus. la définition d'un nouveau temps. donné par l'équation (15), permet d'analyser les essais de puits effectués par la fermeture du puits.

Claims (16)

1. Méthode de détermination des caractéristiques physiques d'un système formé d'un puits et d'une formation souterraine contenant un fluide et communiquant avec ledit puits, ladite formation présentant de l'effet pariétal et/ou le fluide se comprimant ou se décomprimant dans le puits et ladite formation étant homogène ou hétérogène, méthode selon laquelle on provoque un changement de débit dudit fluide et on mesure une grandeur caractéristique de la pression P du fluide à des intervalles de temps successifs Δt. ladite méthode étant caractérisée en ce que l'on compare,
- d'une part, à partir d'un modèle théorique de système de puits-réservoir, l'évolution théorique du logarithme de la dérivée P'D de la pression sans dimension en fonction du logarithme de tD/CD, ladite dérivée P'D étant par rapport à tD/CD, to représentant le temps sans dimension et Co le coefficient sans dimension de l'effet de compression ou de décompression du fluide dans le puits, avec
- d'autre part, l'évolution expérimentale du logarithme de la dérivée ΔP' de la pression en fonction du logarithme des intervalles de temps correspondants Δt, ladite dérivée ΔP' étant par rapport au temps t, et en ce que l'on détermine de la comparaison desdites évolutions théorique et expérimentale au moins une caractéristique du système puits-formation, choisie parmi le produit kh de la perméabilité k par l'épaisseur de ladite formation h, le coefficient C.
2. Méthode selon la revendication 1, caractérisée en ce que ladite évolution théorique est celle du logarithme du produit P'D . tD/CD en fonction du logarithme de tD/CD et ladite évolution expérimentale est celle du logarithme du produit ΔP'. Δt en fonction du logarithme de Δt.
3. Méthode selon la revendication 2, caractérisée en ce qu'en plus, on compare l'évolution théorique du logarithme de PD en fonction du logarithme de tD/CD avec l'évolution expérimentale du logarithme de ΔP en fonction du logarithme de Δt.
4. Méthode selon l'une des revendications 1, 2 ou 3, caractérisée en ce que ladite évolution calculée est également fonction d'une grandeur caractéristique du produit CDe2S et en ce qu'on détermine la valeur du coefficient S de l'effet pariétal.
5. Méthode selon l'une des revendications précédentes caractérisée en ce que, lorsque la formation étudiée présente une double porosité, ladite évolution théorique est en plus fonction des index
Figure imgb0021
et
Figure imgb0022
dans lesquelles λ caractérise le retard de production de fluide par la roche de la formation souterraine comparé à la production de fluide par les fissures de la formation souterraine et ω représente le rapport du volume de fluide produit par lesdites fissures par le volume de fluide produit par le système total et en ce que l'on détermine de la comparaison des évolutions théorique et expérimentale les valeurs de λ et w.
6. Méthode selon l'une des revendications 2 à 5 caractérisée en ce que lorsque ledit changement de débit du fluide correspond à la fermeture du puits, on compare ladite évolution calculée avec l'évolution expérimentale du logarithme de l'expression :
Figure imgb0023
en fonction du logarithme des intervalles de temps Δt, tp étant le temps pendant lesquels le puits a été mis en production.
7. Méthode selon les revendications 1 et 4 caractérisée en ce que l'on trace un graphe théorique de courbes-type en coordonnées cartésiennes et en échelle logarithmique, ledit graphe représentant ladite évolution théorique de la dérivée P'D, en fonction de tD/CD,
8. Méthode selon les revendications 5 et 7 caractérisée en ce qu'on trace en plus deux familles de courbes-type correspondant aux index •
Figure imgb0024
et
Figure imgb0025
9. Méthode selon l'une des revendications 7 et 8 caractérisée en ce qu'on trace une courbe expérimentale en coordonnées cartésiennes et avec la même échelle logarithmique que ledit graphe théorique, ladite courbe expérimentale représentant ladite évolution expérimentale de ΔP' en fonction de Δt, en ce qu'on fait correspondre ladite courbe expérimentale avec l'une des courbes-type dudit graphe théorique et en ce qu'on détermine au moins l'une des caractéristiques kh, C. S, λ et ω par le décalage des axes de coordonnées du graphe théorique et de la courbe expérimentale et par le choix de l'une des courbes-type.
10. Méthode selon la revendication 9, caractérisée en ce qu'on détermine le coefficient C par le décalage des axes d'ordonnées de la courbe expérimentale et du graphe théorique, kh par le décalage des axes d'abscisses de la courbe expérimentale et du graphe théorique et S, w et λ par le choix de la courbe-type du graphe théorique correspondant à la courbe expérimentale.
11. Méthode selon les revendications 2 et 4 caractérisée en ce que l'on trace un graphe théorique de courbes-type en coordonnées cartésiennes et en échelle logarithmique, ledit graphe représentant ladite évolution théorique du produit P' . tD/CD, en fonction de tD/CD.
12. Méthode selon les revendications 3 et 11, caractérisée en ce que l'on superpose audit graphe théorique, un deuxième graphe théorique représentant l'évolution théorique de PD en fonction de tD/CD.
13. Méthode selon les revendications 5 et 11 caractérisée en ce qu'on trace en plus deux familles de courbes-type correspondant aux index
Figure imgb0026
et
Figure imgb0027
14. Méthode selon la revendication 11 ou 13 caractérisée en ce qu'on trace une courbe expérimentale en coordonnées cartésiennes et avec la même échelle logarithmique que ledit graphe théorique, ladite courbe expérimentale représentant ladite évolution expérimentale du produit ΔP' · Δt, en fonction de Δt, en ce qu'on fait correspondre ladite courbe expérimentale avec l'une des courbes-type dudit graphe théorique et en ce qu'on détermine au moins l'une des caractéristiques kh, C. S, \ et ω par le décalage des axes de coordonnées du graphe théorique et de la courbe expérimentale et par le choix de l'une des courbes-type.
15. Méthode selon les revendications 12 et 14, caractérisée en ce qu'en plus on trace une deuxième courbe expérimentale représentant l'évolution expérimentale de AP en fonction de Δt, en ce qu'on fait correspondre ladite deuxième courbe expérimentale avec l'une des courbes-type dudit deuxième graphe théorique et en ce qu'on détermine au moins l'une des caractéristiques kh, C et S, par le décalage des axes de coordonnées du deuxième graphe théorique et de la deuxième courbe expérimentale et par le choix de l'une des courbes-type.
16. Méthode selon les revendications 14 ou 15, caractérisée en ce qu'on détermine le coefficient kh par le décalage des axes d'ordonnées de la courbe expérimentale et du graphe théorique, C par le décalage des axes d'abscisses de la courbe expérimentale et du graphe théorique et S par le choix de la courbe-type du graphe théorique correspondant à la courbe expérimentale.
EP84400781A 1983-04-22 1984-04-19 Méthode de détermination des caractéristiques d'une formation souterraine produisant un fluide Expired EP0125164B1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR8307075 1983-04-22
FR8307075A FR2544790B1 (fr) 1983-04-22 1983-04-22 Methode de determination des caracteristiques d'une formation souterraine produisant un fluide

Publications (2)

Publication Number Publication Date
EP0125164A1 EP0125164A1 (fr) 1984-11-14
EP0125164B1 true EP0125164B1 (fr) 1986-12-30

Family

ID=9288354

Family Applications (1)

Application Number Title Priority Date Filing Date
EP84400781A Expired EP0125164B1 (fr) 1983-04-22 1984-04-19 Méthode de détermination des caractéristiques d'une formation souterraine produisant un fluide

Country Status (6)

Country Link
US (1) US4597290A (fr)
EP (1) EP0125164B1 (fr)
CA (1) CA1209699A (fr)
DE (1) DE3461844D1 (fr)
FR (1) FR2544790B1 (fr)
NO (1) NO841473L (fr)

Families Citing this family (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8418429D0 (en) * 1984-07-19 1984-08-22 Prad Res & Dev Nv Estimating porosity
FR2569762B1 (fr) * 1984-08-29 1986-09-19 Flopetrol Sa Etu Fabrications Procede d'essai de puits d'hydrocarbures
FI75631C (fi) * 1984-11-20 1988-07-11 Reijonen Veli Oy Foerfarande foer dimensionering av grundvattensbrunn.
NO170037C (no) * 1985-07-23 1992-09-02 Flopetrol Services Inc Fremgangsmaate for maaling av stroemningshastigheter i en borebroenn.
FR2585404B1 (fr) * 1985-07-23 1988-03-18 Flopetrol Procede de determination des parametres de formations a plusieurs couches productrices d'hydrocarbures
EP0429078A1 (fr) * 1986-05-15 1991-05-29 Soletanche Procédé pour la mesure de la perméabilité d'un terrain
FR2613418B1 (fr) * 1987-04-02 1995-05-19 Schlumberger Cie Dowell Procede de traitement matriciel dans le domaine petrolier
US4797821A (en) * 1987-04-02 1989-01-10 Halliburton Company Method of analyzing naturally fractured reservoirs
US4843878A (en) * 1988-09-22 1989-07-04 Halliburton Logging Services, Inc. Method and apparatus for instantaneously indicating permeability and horner plot slope relating to formation testing
US5050674A (en) * 1990-05-07 1991-09-24 Halliburton Company Method for determining fracture closure pressure and fracture volume of a subsurface formation
US5005643A (en) * 1990-05-11 1991-04-09 Halliburton Company Method of determining fracture parameters for heterogenous formations
US5517593A (en) * 1990-10-01 1996-05-14 John Nenniger Control system for well stimulation apparatus with response time temperature rise used in determining heater control temperature setpoint
US5247829A (en) * 1990-10-19 1993-09-28 Schlumberger Technology Corporation Method for individually characterizing the layers of a hydrocarbon subsurface reservoir
FR2710687B1 (fr) * 1993-09-30 1995-11-10 Elf Aquitaine Procédé d'évaluation de l'endommagement de la structure d'une roche entourant un puits.
US5501273A (en) * 1994-10-04 1996-03-26 Amoco Corporation Method for determining the reservoir properties of a solid carbonaceous subterranean formation
FR2747729B1 (fr) * 1996-04-23 1998-07-03 Elf Aquitaine Methode d'identification automatique de la nature d'un puits de production d'hydrocarbures
US6832515B2 (en) 2002-09-09 2004-12-21 Schlumberger Technology Corporation Method for measuring formation properties with a time-limited formation test
US7178392B2 (en) 2003-08-20 2007-02-20 Schlumberger Technology Corporation Determining the pressure of formation fluid in earth formations surrounding a borehole
US7031841B2 (en) * 2004-01-30 2006-04-18 Schlumberger Technology Corporation Method for determining pressure of earth formations
DE602005013483D1 (de) * 2005-02-28 2009-05-07 Schlumberger Technology Bv Verfahren zur Messung von Formationseigenschaften mit einem Formationstester
US7793718B2 (en) * 2006-03-30 2010-09-14 Schlumberger Technology Corporation Communicating electrical energy with an electrical device in a well
US7712524B2 (en) * 2006-03-30 2010-05-11 Schlumberger Technology Corporation Measuring a characteristic of a well proximate a region to be gravel packed
US8056619B2 (en) 2006-03-30 2011-11-15 Schlumberger Technology Corporation Aligning inductive couplers in a well
US7735555B2 (en) * 2006-03-30 2010-06-15 Schlumberger Technology Corporation Completion system having a sand control assembly, an inductive coupler, and a sensor proximate to the sand control assembly
US20080230221A1 (en) * 2007-03-21 2008-09-25 Schlumberger Technology Corporation Methods and systems for monitoring near-wellbore and far-field reservoir properties using formation-embedded pressure sensors
US7580797B2 (en) * 2007-07-31 2009-08-25 Schlumberger Technology Corporation Subsurface layer and reservoir parameter measurements
US8121790B2 (en) * 2007-11-27 2012-02-21 Schlumberger Technology Corporation Combining reservoir modeling with downhole sensors and inductive coupling
US8136395B2 (en) 2007-12-31 2012-03-20 Schlumberger Technology Corporation Systems and methods for well data analysis
US8078402B2 (en) * 2008-07-16 2011-12-13 Schlumberger Technology Corporation Method of ranking geomarkers and compositional allocation of wellbore effluents
US20100147066A1 (en) * 2008-12-16 2010-06-17 Schlumberger Technology Coporation Method of determining end member concentrations
US8839850B2 (en) * 2009-10-07 2014-09-23 Schlumberger Technology Corporation Active integrated completion installation system and method
US20110192596A1 (en) * 2010-02-07 2011-08-11 Schlumberger Technology Corporation Through tubing intelligent completion system and method with connection
US9249559B2 (en) 2011-10-04 2016-02-02 Schlumberger Technology Corporation Providing equipment in lateral branches of a well
US9644476B2 (en) 2012-01-23 2017-05-09 Schlumberger Technology Corporation Structures having cavities containing coupler portions
US9175560B2 (en) 2012-01-26 2015-11-03 Schlumberger Technology Corporation Providing coupler portions along a structure
US9938823B2 (en) 2012-02-15 2018-04-10 Schlumberger Technology Corporation Communicating power and data to a component in a well
US20130282286A1 (en) * 2012-04-20 2013-10-24 Chevron U.S.A. Inc. System and method for calibrating permeability for use in reservoir modeling
US10036234B2 (en) 2012-06-08 2018-07-31 Schlumberger Technology Corporation Lateral wellbore completion apparatus and method
GB2499523B (en) * 2013-02-15 2014-04-09 Petroleum Experts Ltd Modelling of transient hydrocarbon reservoirs
CN114060018B (zh) * 2020-08-04 2024-05-28 中国石油天然气股份有限公司 一种储层动态储量确定方法、系统、设备及可读存储介质

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3550445A (en) * 1968-01-19 1970-12-29 Exxon Production Research Co Method for testing wells for the existence of permeability damage
US3604256A (en) * 1969-01-31 1971-09-14 Shell Oil Co Method for measuring the average vertical permeability of a subterranean earth formation
US3636762A (en) * 1970-05-21 1972-01-25 Shell Oil Co Reservoir test
US4328705A (en) * 1980-08-11 1982-05-11 Schlumberger Technology Corporation Method of determining characteristics of a fluid producing underground formation

Also Published As

Publication number Publication date
US4597290A (en) 1986-07-01
FR2544790B1 (fr) 1985-08-23
CA1209699A (fr) 1986-08-12
NO841473L (no) 1984-10-23
FR2544790A1 (fr) 1984-10-26
DE3461844D1 (en) 1987-02-05
EP0125164A1 (fr) 1984-11-14

Similar Documents

Publication Publication Date Title
EP0125164B1 (fr) Méthode de détermination des caractéristiques d&#39;une formation souterraine produisant un fluide
CN107701180B (zh) 一种基于密闭取心的原始油藏含水饱和度计算方法
EP0113285B1 (fr) Procédé et dispositif pour déterminer les caractéristiques d&#39;écoulement d&#39;un fluide dans un puits à partir de mesures de températures
EP0698722B1 (fr) Méthode pour tester de formations à faible perméabilité
US5247830A (en) Method for determining hydraulic properties of formations surrounding a borehole
EP1167948B1 (fr) Méthode pour évaluer des paramètres physiques d&#39;un gisement souterrain à partir de débris de roche qui y sont prélevés
EP1963886B1 (fr) Methode d&#39;evaluation quantitative des pressions de fluides et de detection des surpressions d&#39;un milieu souterrain
FR2853071A1 (fr) Methode et dispositif pour evaluer des parametres physiques d&#39;un gisement souterrain a partir de debris de roche qui y sont preleves
BRPI0821665B1 (pt) método para determinar confiança em medidas tomadas por uma ferramenta de teste
CA2624305A1 (fr) Procedes et systemes permettant de determiner des proprietes de reservoir propres aux formations souterraines
CA2517244A1 (fr) Methode de determination de parametres d&#39;ecoulement multiphasique d&#39;un milieu poreux avec prise en compte de l&#39;heterogeneite locale
EP1506388A1 (fr) Methode d evaluation de la courbe de pression capillaire des roches d un gisement souterrain a partir de mesures sur des debris de roche
Novakowski Analysis of pulse interference tests
EP0174890B1 (fr) Procédé d&#39;essai de puits d&#39;hydrocarbures
EP1522850A2 (fr) Méthode et dispositif pour mesurer l&#39;anisotropie de résistivité d&#39;échantillons de roche présentant des litages
EP3327435A1 (fr) Procede d&#39;exploitation et/ou de surveillance d&#39;un aquifere comportant au moins un gaz dissous
EP3626929A1 (fr) Procede d&#39;exploitation d&#39;un gisement d&#39;hydrocarbures par injection d&#39;un polymere
EP2128653A2 (fr) Procédé pour localiser l&#39;origine spatiale d&#39;un évènement sismique se produisant au sein d&#39;une formation souterraine
EP2315055B1 (fr) Méthode pour interpréter des enregistrements sismiques répétitives
US20210148227A1 (en) Relating To Injection Wells
WO2009153443A1 (fr) Methode pour evaluer des pressions de fluides et detecter des surpressions dans un milieu souterrain
FR2831917A1 (fr) Procede de determination de la variation de la permeabilite relative a au moins un fluide d&#39;un reservoir contenant des fluides en fonction de la saturation en l&#39;un d&#39;entre eux
EP3502751A1 (fr) Procédé de modelisation d&#39;un bassin sédimentaire
FR2585404A1 (fr) Procede de determination des parametres de formations a plusieurs couches productrices d&#39;hydrocarbures
FR2836227A1 (fr) Methode pour evaluer des parametres physiques d&#39;un gisement souterrain a partir de debris de roche qui y sont preleves

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

AK Designated contracting states

Designated state(s): DE GB IT NL

17P Request for examination filed

Effective date: 19850214

17Q First examination report despatched

Effective date: 19860220

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): DE GB IT NL

REF Corresponds to:

Ref document number: 3461844

Country of ref document: DE

Date of ref document: 19870205

ITF It: translation for a ep patent filed

Owner name: DOTT. FRANCO CICOGNA

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed
ITTA It: last paid annual fee
PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: NL

Payment date: 19910430

Year of fee payment: 8

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DE

Payment date: 19910529

Year of fee payment: 8

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NL

Effective date: 19921101

NLV4 Nl: lapsed or anulled due to non-payment of the annual fee
PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DE

Effective date: 19930101

REG Reference to a national code

Ref country code: GB

Ref legal event code: IF02

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20020417

Year of fee payment: 19

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GB

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20030419

GBPC Gb: european patent ceased through non-payment of renewal fee

Effective date: 20030419