PROCESS FOR MEASURING FLOW AND DETERMINING THE PARAMETERS OF MULTILAYER HYDROCARBON-PRODUCING FORMATIONS

Abstract: The invention concerns a process for measuring the flow as a function of time of the several layers of a subterranean multilayer hydrocarbon-producing formation through which a well is drilled, and for determining the parameters of formation by comparison of trends in well behavior shown by actual measurements with trends in behavior established theoretically. The flow rate contributions of individual layers (1 to 5) are determined from cumulative measurements made above successive layers by a flowmeter (13) moved vertically within a wellbore. Experimental curves are derived that represent the relative variations in time of the flowrates of layers ( q j ) or groups of layers in the formation using flowrate measurement points obtained when a total flowrate variation ( Q) has been imposed on the well at a given time between two given constant values. These experimental curves then compared to theoretical model curves (G, H) established for various values of the characteristic parameters of a subterranean formation, to determine values of the parameters of the formation.

Description: This invention involves a process for measuring the flow as a function of time of the several layers of a subterranean multilayer hydrocarbon-producing formation through which a well is drilled, and of the formation. Measurements of pressure in oil wells as a function of time in order to determine the characteristics of the productive subterranean formations through which the wells are drilled, have long been known. Although such measurements make it possible to determine a considerable number of parameters characterizing subterranean formations in general, they are insufficient in the case of complex reservoirs such as multilayer formations. A single pressure curve cannot in effect supply the data necessary for the determining the characteristics specific to the various layers, such as their permeability and skin coefficient. Patent FR-A-2,544,790 describes a method for evaluating a single layer reservoir. Downhole pressure measurements are performed after a step change of the surface flowrate. The logarithm of the pressure derivative is determined and compared to type curves derived from a single layer model to determine characteristic parameters such as the skin and permeability of the total reservoir. A process for testing multilayer systems was proposed by Gao ("The Crossflow Behavior and the Determination of Reservoir Parameters by Drawdown Tests In Multilayer Reservoirs", SPE paper No. 12580, submitted for publication September 29, 1983). Using the semipermeable wall model published by Deans and Gao in SPE paper No. 11966 presented at the 58th Annual Conference and Exposition at San Francisco, October 5-8, 1983, this process consists of testing each layer individually and recording a series of pressure curves. Such a process involves at least three inconveniences. First, it takes a long time. Second, the interpretation of the curves is tricky if there is any transfer flow between formation layers. Finally, during testing, the well is never in an activity mode similar to a real production situation. Another method of investigating multilayer systems which is described in French patent FR-A-2,434,923 is to use variations of flow and pressure as a function of depth in a stabilized well, i.e., a well In which production is at a constant surface pressure and flowrate. This type of measurement leads to a "snapshot" of the flow and pressure at each layer for a given surface flowrate and pressure. The data obtained can be presented for various successive surface flowrates in the form of a series of pressure/flow curves for each layer. Here, there are two inconveniences. First, not all wells reach a stabilized flow situation. In addition, it was shown (Lefkovits, H.C., Hzebroek, P., Allen, E.E. and Matthews, C.S.: "A Study of the Behavior of Bounded Reservoirs Composed of Stratified Layers", J. Pet. Tech., March 1961) that the respective flowrates of the layers vary with time. Thus, this process is applicable only to wells which actually reach a steady state. Based on the state of the art thus recalled, the purpose of the invention is an original process for determining the characteristic parameters of a subterranean hydrocarbon-producing formation through which a well is drilled, as set forth in claim 1. Preferred embodiments of the process are set forth in the dependent claims. Such a process makes it possible to determine the characteristic parameters of a subterranean multilayer formation, based mainly on the known fact that variations of the respective flowrates of the layers of the formation are, during an initial period immediately following the well flowrate change, sensitive to wall skin and layer permeability effects and, in a later period, to interlayer fluid transfer effects. Other characteristics and advantages of the invention will become more clearly apparent from the following description and attached drawings of a non-limitative example. Figure 1 is a schematic vertical cross-section of an oil well drilled in a multilayer formation into which a flowmeter has been lowered. Figure 2 represents a curve obtained by moving the flowmeter in the well. Figures 3 and 4 show two sets of flow curves prepared based on curves such as those in Figure 2. Figure 5 represents an experimental well pressure/time curve and the derivative curve. Figures 6 through 8 represent flowrate relative variation curves for layers relative to total well flow, distinct groups of layers or zones relative to total well flow and layers relative to the total flowrate of the zone to which the layers belong, respectively. Figure 1 shows an oil well 10 drilled in a formation containing several oil-bearing layers, i.e., five layers, 1, 2, 3, 4 and 5. The intermediate layers 12, 34, 45 separating layers 1 and 2, 3 and 4, and 4 and 5 respectively, have a certain vertical permeability such that there may be oil flow through these intermediate layers. On the other hand, layer 23 between layers 2 and 3 is impermeable and there is not oil flow between these two layers. Each group of layers between which vertical oil flow can occur and which is isolated by impermeable layers is called a "zone". In this example, the formation includes two zones Z1 and Z2, Z1 being composed of layers 1 and 2 and Z2 of layers 3, 4 and 5. By definition, there cannot be any vertical transfer of oil between two zones. This division of the subterranean formation into layers and zones has proven quite advantageous for the interpretation of the results obtained and is one of the characteristics of this invention. The various layers and zones can be identified by making a recording of the flowrate as a function of depth through the formation. Previously made logs of the well can also be used. When the well is placed in production, it delivers a total oil flow to the surface through its production string 11. The annular space between the casing 10 and the production string 11 is sealed off by the packer 9. The partial flowrates q1 to q5 of layers 1 through 5 make up the total flowrate. The total flowrate Q measured above the formation is the sum of flowrates q1 to q5. It should be noted that the flowrate measured on the surface may differ due to the wellbore storage effect. These flowrates are identical if the effect is zero. The process according to the invention uses the variations in time of the pressure p in the well and of the partial flowrates q1 to q5 of the various layers resulting from a modification made to the total well flowrate Q. The flow measurements are made using a flowmeter 13 (for example, as described in French patent No. 74 / 22 391) lowered into the well at the end of a cable 14, then moved vertically several times in a sweeping movement throughout the total depth of the formation. When it is immediately above a layer, the flowmeter measures the cumulative flowrate of that layer and those below it. With each passage, a curve such as the one shown in Figure 2 is recorded, indicating the flowrate measured as a function of depth, from which the cumulative flowrates of the various layers 1 through 5, i.e., q5, q5 + q4, q5 + q4 + q3, etc., recorded in front of intermediate layers 45, 34, 23, etc., can be deduced. The times at which these flow measurements are made are also recorded. This makes it possible to trace the curves representing the variations in cumulative flow as a function of time. Figure 3 is a semi-logarithmic representation of such a set of curves, with the flowrates expressed in barrels per day (1 barrel = 158.98 litres). Using these curves and simple subtraction, it is possible to trace the set of curves in Figure 4, which represent the variations of the flowrate specific to each layer 1 through 5 as a function of time. In Figures 3 and 4, the measurement times seem to be the same for all the layers. In fact, because of the sweeping movement of the flowmeter 13, these points are offset from one curve to the next. This obviously has no effect on the curve tracing operations. The cable used may or may not be electrical. If it is, the data from the flowmeter are transmitted through the cable to the surface to be recorded and processed. When the cable used is a simple "piano string", the data are recorded by a downhole recorder with memories. Such a recorder is, as an example, described in British patent application No. 82 31560. It may be helpful to use several flowmeters connected end-to-end so as to record the respective flowrates of several layers at once or that of a single layer at very short intervals of time. Pressure measurements are made using a pressure gauge (for example, as described in the French patent published under the No. 2 496 884) which can be installed stationary either at the well-head or (as represented at 16) at the top of the formation or be connected to the flowmeter 13 (at 16'). In the latter case, it must be kept in mind that the pressure gauge is subjected to the pressure of an oil column of variable height. Measurements of the pressure p in the well as a function of time are obtained in this way. Like the flowrate measurements, the pressure measurements are transmitted by electrical cable to the surface to be recorded or are recorded in the well using a recorder. In its measurement phase, the process according to the invention consists essentially of varying the well flowrate Q by a quantity ΔQ at a time t1, and measuring the pressure and the flowrates of the respective layers of the formation just prior to time t1, then for a certain period thereafter. The measurements taken after time t1 make it possible to prepare the sets of flow curves in Figures 3 and 4 and the pressure curve in Figure 5 as a function of the time Δt lapsed after time t1. More precisely, Figure 5 represents the variations of the quantity Δp x (Q/ΔQ) with Δp designating the pressure differential measured between times t1 and t1 + Δt. The pressure scale is expressed in psi (1 psi = 6.9 kPa approx.). The variations of the derivative of the aforesaid quantity Δp are also represented. The abscissa and ordinate scales are logarithmic. The table at the end of the description gives an example of simulated values of variations of pressure Δp (in psi) as a function of time Δt (Δt being expressed in days and counted from time t1). The values of the variations of the derivative (Δp)' are also indicated, calculated as explained below, as are the flowrate measurements q1 to q5 (expressed in barrels/day and represented on Figure 4) of the five layers of the formation. The derivative (Δp)' is calculated as a function of the logarithm of Δt, i.e.: The method of calculation and the interpretation of the derivative data are described in the published French patent application No. 83 07075 dated 22 April 1983. The negative values obtained for Δp and (Δp)' can be explained by the overlay principle which is well known to specialists. Briefly, in order for the measurements to be usable and meaningful, the well flow time prior to the flowrate change must be very long compared to the period of time during which measurements are made after the flowrate change. In order to plot the curves in Figure 5, the valves of Δp and (Δp)' from the table are multiplied by the ratio Q/ΔQ of the flowrate Q before the change at time t1 to the variation in the flowrate ΔQ before and after time t1. In the example in Figure 5, Q = 500 and ΔQ = 500 - 200 = 300, as the flowrate after time t1 was reduced from 500 to 200 barrels/day. This is equivalent to normalizing the pressure values after time t1 with the values that would have been obtained before time t1. The normalization of the curves is important as regards the pressure measurements as well as the flowrate measurements because it makes possible the use of the values measured just before time t1 as asympototic values for very long time periods Δt. This characteristic of the invention, which is important in practice, will be explained in connection with Figure 7 (points P1 and P2). For the interpretation of the experimental pressure data, a classical analysis, well known to specialists, is made, consisting of plotting various logarithmic and semilogarithmic graphs to diagnose the wellbore storage effect; the oil flow regime in the reservoir, which can be radial and considered to be infinite at the scale of the well; the presence of several productive layers and the presence of possible reservoir limits. Thus, in logarithmic scales, the wellbore storage effect is shown by a slope equal to 1 for the pressure and pressure derivative curves for short time periods (beginnings of curves) and the presence of a limit or boundary to the reservoir is shown by an increased in the pressure and pressure derivative values for long periods (ends of curves). These diagnostic methods are commonly used in the petroleum industry and are described, for example, in U.S. patent 4 328 705 and published French patent application No. 83 07075. The convolution of the total flow measured at the bottom of the well with the pressure can also be used to eliminate the wellbore storage effect from the pressure measurements (in this case, the pressure is called rate convolved pressure). This technique is published in "Interpretation of Pressure Build-up Test using In-situ Measurement of Afterflow", Journal of Petroleum Technology, January 1985. In its measurement interpretation phase, the process according to the invention includes the following parameter determination operations: A) kh (average product of formation permeability k x thickness h for the overall formation); B) k j and s j (horizontal permeability and skin coefficient of layer j, with j varying from 1 to 5 in this example); C) type and position of the external limit or boundary of the formation (which determines the extent and type of the formation); D) vertical permeability between layers. Based on pressure measurements, using the following formula: in which: ΔQ is the change made in the well flowrate (expressed in barrels/day) at the time t₁; B is the relative volume factor of the oil in the formation and at the surface (equal to the ratio between the volumes of oil in the formation and at the surface); µ is the viscosity of the oil expressed in centipoises; (Δp)' M is the value of the derivative of the pressure p as a function of the logarithm of time in the flat part of the derivative curve (Figure 5). This flat portion reflects an infinite action radial flow. In this example, Figure 5 shows that: (Δp)' M . (Q/ΔQ) = 5 Thus, (Δp)' M = 3. In addition, other measurements revealed that: B = 1.2 and µ = 1. Thus, kh = 8.472 md.ft = 25.42 x 10² µm².m. B. Determination of k j and s j For each layer j, the curve representing as a function of time the fraction of the variation of total flow attributable to the layer involved, i.e., the quantity Δq j /ΔQ, based on the values from the table and the curves in Figure 4. This results in five series of points in semilogarithmic representation (Figure 6) for the five layers 1 through 5 of the formation, respectively, as a function of the time Δt lapsed after the time t1. Each series of points is then compared with a theoretical model to determine which of the curves suitably fits the series of points involved, at least during the initial period following the time t1 of the flowrate change. It has been acknowledged that during this period the model used may correspond to the absence of flow between layers in the formation, with an infinite external boundary. After the initial period, deviations may appear between the measurement points and the theoretical curve due to interlayer flow and external boundary effects, or overlay effects. The theoretical model was established based on the following formula: in which: q jD is the Laplace transform of the dimensionless flowrate of layer j; s j is the skin coefficient of layer j; K₀ and K₁ are the modified Bessel functions of the first and second types; P wD is the Laplace transform of the dimensionless pressure in the well; while σ j is given by: σ = (ω j z/k j ) 1/2 in which: with: in which (h) j designates the product porosity x height of layer j, n the number of layers in the formation and z the Laplace space variable. Equation (1) does not give flowrate as a function of time. To obtain it, the inverse Laplace transform given by the Stehfest algorithm is applied (see "Numerical inversion of Laplace transforms", D-5, Communications of the ACM, January 1970, No. 1, pages 47 to 49). When fitting is achieved with a given curve of the theoretical model, the skin coefficient s j of the layer j involved, appearing in formula (1), can be deduced from it, as can its permeability, which also appears in formula (1) as the product (kh) j of the permeability and height of said layer, the latter parameter being known by previously made log measurements, while the product kh is determined using the pressure measurements explained above. To illustrate cases which could be encountered in practice, two theoretical curves G and H (dotted lines) are shown in Figure 6 which do not correctly fit the series of measurement points for formation layers 1 and 2 on the left side of the figure, while there is a good fit on the right side (significant deviations do, however, occur at the extreme right due to boundary and interlayer flow effects). The examination of the position of curve G shows, for example, that the skin coefficient selected for it is too low and should be increased. For curve H, the opposite is true: the skin effect must be reduced, even though a modification of the value of the skin effect of curve G has an influence on the other curves. These operations make it possible to determine the horizontal permeability k j and the skin coefficient s j of each layer of the formation. Using a different operational mode of this invention, these parameters can be determined as explained below: It is known to specialists that the result of the mathematical operation of convolution of the derivative of the variations of flow q with dimensionless pressure P D (the pression that would be obtained if no other parameters intervened in the formation and the well to influence the pressure value and if the flow rate were constant) represents the variations of the pressure P sf effectively measured in the well in front of the formation. This is expressed by the following equation: P sf (T) being the value of the pressure variation measured in the well at time T. To obtain P D , which is the pressure value being sought, requires the mathematical deconvolution between the effectively measured pressure P sf and the flow. However, the results obtained by deconvolution can be sprinkled with significant errors if the experimental data include some noise. Convolution is thus the preferred operation. This is why, within the framework of this invention, the flow variations for each layer and the pressure variations in the well were measured. It was then shown that the convolution of the flow variations for each layer with the pressure variations in the well provides the pressure response of the layer as if it were the only one producing a fluid, provided, however, that there is no interlayer flow. Thus, armed with the pressure response of each layer, it is possible to use the classical methods of interpretation for each, specifically the pressure/time curves plotted on semilogarithmic scales which make it possible to determine the permeability and skin effect. C. Determination of the external boundary of each zone The point is to determine the boundary type of each zone: seemingly infinite boundary; no-flow boundary, behaving like an impermeable seal, with all the liquid flowing into the well coming from the formation zone located inside this boundary; Constant pressure boundary. In the first instance, it is as though there were no boundary. In the other two cases, the radius of the boundary must also be specified. For this determination, a graph (Figure 7) similar to Figure 6 is prepared, in which each series of points corresponds to a zone i of the formation based on the above definition, and no longer to a given layer (of course, a zone may contain only one layer). In this example, there are two zones Z1 and Z2 (Figure 1) and the two series of points represent, respectively, the following quantities: as a function of the time Δt. The values from the table are used to plot the experimental curves in Figure 7. In addition, theoretical models corresponding to the above formula (1), as well as the following formulas, were used: In these formulas, formula (2) being already known, the quantities are defined as follows: in which I₀, I₁, K₀ and K₁ are modified Bessel functions of the first and second type and r eD is the dimensionless external radius of the formation. Formula (1) refers to the case of a boundary behaving as if it were infinite, formula (2) a no-flow boundary and formula (3) a constant pressure boundary. Figure 7 shows the fitting achieved (in the initial period) between the series of points corresponding to zones Z1 and Z2 and the curves defined based on formula (2). It can be concluded that the boundary of the zones under study is of the "no-flow" type. If the fitting is achieved with a curve ending with a horizontal, such as curve K corresponding to formula (3), sketched at the top right of Figure 7, the boundary is of the "constant pressure" type. If the curve ends with a slight downward bend, like curve L resulting from formula (1), the boundary is of the apparently infinite type. Indeed, in these operations relative to the various zones of the formation, there is no need to envision a vertical transfer flow situation, since such transfers are by definition inexistent between zones. In the case represented in Figure 7, showing a no-flow boundary, i.e., a situation in which the production volume of the formation is limited, the straight portion of the curve for each zone tends towards a value equal to the product h for the zone, i.e., in this example, 0.4 for zone Z1 and 0.6 for zone Z2. These values are in addition known due to previous logging operations. Figure 7 also shows that the straight portion of the curves deviates from the experimental points. This is because of an overlay effect due to the fact that the well in question was placed in production for 200 hours (8.33 days) and then its production flow rate was reduced (from 500 to 200 barrels/day) at time t1, for another 200 hours. However, if a measurement is made at the end of the first 200-hours period, just prior to time t1, the results obtained can be recorded on the figure (points P1 and P2) and considered as measurement points obtained after the initial period, after time t1, without overlay effect. As can be seen, the theoretical curves are very close to these points. It results from the foregoing remark that, in practice, there is no need to make measurements at times distant from time t1, since measurements made just prior to that time can advantageously replace them. Thus, the measurements at the points located approximately between 10⁰ and 10¹ days after t1 need not be made, and this considerably shortens the total time required for well measurements. As a result, it is possible to determine using only the measurements taken just before t1 whether the boundary is of the no-flow type (if the measurement points correspond to the known values of h) or not (if the measurement points do not correspond), since this effect depends solely on boundary conditions. If the boundary is recognized to be non-infinite, its radius is determined by finding the radius value which leads to the best fit between the model curves and the measurement curves in the period following the initial period. The measurement points recorded just prior to the change in the well's flow rate are also very useful in this phase of determination of formation parameters. As in the case of the determination of permeability and skin effect, there is a second possible method using convolution operations. It has, in fact, been shown that the convolution of the flow variations of each zone with the pressure variations in the well provides the pressure response of the zone involved. As in the case of the individual layers, this harks back to a classical well test interpretation, particularly for the determination of the boundary of each zone. D. Determination of interlayer permeability Having determined the horizontal permeability and the skin coefficient of each layer of the formation, the type and location of the zone boundaries, the vertical permeability between layers remains to be determined. This is done by means of an analysis in each zone of the formation of the flow rates of the layers of the zone as compared to the zone's total flow rate. As is shown in Figure 8, we show as a function of Δt, still using a semilogarithmic representation and based on the table data, the quantities Δq i / j/ΔQ i in which ΔQ i designates the flow rate variation of zone i and Δqij the flow rate variation of layer j belonging to zone i. Figure 8, as an example, is limited to the measurements for zone Z1, composed of layers 1 and 2, the values for which are indicated in the table (page 19). The theoretical model used, established for the case in which there are transfer flows between layers, is derived from the following formula: in which: q jD is the dimensionless flow rate of layer j of zone i, which contains m i layers; C D is the dimensionless wellbore storage constant; K j = ratio of the product permeability x height for layer j to the average product kh permability x height for zone i; σ k1 represents the quantity σ k for zone i; σ k are the roots of the equation γ n = 0 in which γ j is a polynomial defined by recurrence by: γ i = γ j-1 a jj - γ j-2 a j,j-1 a j-1,j for j = 2, ..., n with γ₀ = 1 and γ₁ = a₁₁, a jk designating the elements of the matrix [a jk ]: α ki / j is a coefficient relative to layer j, root k, for zone i, defined by the formula: b ki is an external boundary condition coefficient defined by the formulas: b ki = 0 for an apparently infinite boundary for a no-flow boundary for a constant pressure boundary, while B ki / 1 is related to A ki / 1 by the equation: A ki / 1 being determined based on well conditions. The skin effect coefficient values obtained in interpretation phase B) are used, keeping in mind the type and location of the formation's external boundary as determined in phase C). Finally, a set of values is sought for the parameters λ j of interlayer permeability between layers j and j + 1 of each zone i such as to achieve good fit of the curves for all the Δq j / ΔQ i ratios considered. More precisely, it can be noted that the appearance of the curves in the left-hand portion of the figure depens on permeability and skin effect, while the right-hand side of Figure 8 depens also on the type of boundary and transfer flows. Since the permeability, skin effect and boundary type are known, the only remaining parameter is transfer flow, for which different values are tried until a good curve fit is achieved. These operations are repeated for each of the zones of the formation, so as to determine the transfer flow parameters for all layers. The set of calculations and curve fitting operations just described as part of the process according to the invention can be done my hand or, preferably, by a digital calculator. In the first instance, sets of typical curves are traced using the equations given above. These sets of curves are a graphic representation of the behavior of the theoretical models. A digital calculator can also be used to select the values of the parameters being sought which correspond to a perfect fit between the theoretical and experimental variations of the various functions of the pressure and flow rates (variation of pressure, of the derivative of pressure, of the fraction of the variation of the total flow rate for a given layer and for a given zone and of the fraction of the variation of flow rate of a layer as compared to the flow rate of the zone to which it belongs, all as a function of time).

Claims: A process for determining characteristic parameters such as permeability and skin coefficient of a subterranean hydrocarbon-producing formation through which a well is drilled, comprising the steps of: determining relative variations of pressure vs. time, in response to changes imposed on the surface total flowrate of the well, from pressure data measured in the well before a change in flowrate and at a plurality of times starting immediately after such change, and computing a composite value (kh) of said parameters, for the formation analyzed as a single layer, from a comparison of the experimental pressure variations thus determined with theoretical variations derived from a model of the system made up of a single producing formation and a well traversing the formation, by deducing said composite value from the parameter value of the model which best fits the experimental variations, characterized in that, in order to determine the parameters of individual layers or groups of layers constituting the formation, it further comprises the step of: determining the relative variations vs. time, in response to said changes of the surface total flowrate, of the flowrate contributed by each layer or group of layers from measurements of the flowrate above each layer or group of layers taken before a change and at a plurality of times starting immediately after said change; and computing individual values (k j ,s j ) of said parameters, for each layer or group of layers, from a comparison between the experimental flowrate variations thus determined and theoretical flowrate variations corresponding to a model made up of several layers or several groups of layers and of a well traversing said layers, by deducing said individual values from the parameter values for which a best fit is obtained with the experimental flowrate variations. A process according to claim 1, characterized in that: the flowrate determining step comprises, for each layer j of the formation, determining as a function of time the variations of the fraction Δq j /ΔQ representing the ratio of the flowrate variation of layer j to the variation of the total downhole well flowrate; and the step of computing said individual values comprises a first comparison between the variations of the fraction Δq j /ΔQ and theoretical flowrate variations corresponding to the model to deduce the horizontal permeability and skin coefficient parameters k j and s j of each layer j. A process according to claim 2, characterized in that said first comparison is made between the variations of the fraction Δq j /ΔQ during an early period following a surface flowrate change and theoretical layer flowrate variations corresponding to a model associated with various values of permeability and skin coefficient. A process according to anyone of claims 2 to 3, characterized by further comprising the steps of: for each zone i of the formation representative of a group of productive layers confined between two impermeable intermediate layers, determining as a function of time the variations of the fraction Δq i /ΔQ representing the ratio of the composite flowrate variations of said zone i to the variation of the total well flowrate; and making a second comparison between the variations of the fraction Δq i /ΔQ and theoretical flowrate variations of the model to deduce the type and position of the external boundary of said zone i from said second comparison. A process according to claim 4, characterized in that said second comparison is made between the time variations of the fraction Δq i /ΔQ during a late period following said surface flowrate change and theoretical flowrate variations of zones having various types of external boundary. A process according to claim 4 or 5, characterized by further comprising the steps of: for each layer j of a zone i of the formation, determining as a function of time the variations of the fraction Δq j /ΔQ i representing the ratio of the flowrate variations of said layer j to the composite flowrate variations of the zone i; and making a third comparison between the variations of the fraction Δq j /ΔQ i with corresponding theoretical flowrate variations of the model to deduce the interlayer permeability parameters of the zone i from said third comparison. A process according to claim 6, characterized in that said third comparison is made between the variations of the fraction Δq j /ΔQ i during a late period following said surface flowrate change and the theoretical flowrate variations of layers of the zone i which have various interlayer permeability parameters.

Inventors: EHLIG-ECONOMIDES, CHRISTINE
Applicants: FLOPETROL SERVICES, INC.
Application year: 1986
Main classification: E21B47
IPC codes: E21B49, E21B47