GB2178170A - Evaluation of formation invasion and formation permeability - Google Patents

Abstract

Formation invasion and formation permeability are evaluated by conducting a series of formation resistivity loggings in the operation of a measurement while drilling (MWD) logging system. A first resistivity measurement is taken within a short known time of drilling, and that at least one subsequent measurement is taken so as to obtain incremental data of the rate of invasion and hence, the permeability. <IMAGE>

Description

SPECIFICATION Evaluation of formation invasion and formation permeability This invention relates to the field of borehole logging, and in particular to a process for evaluating formation invasion and formation permeability through resistivity logging.

The principal objectives of well logging are to identify and correlate underground formations, to determine the minerology and physical properties of potential reservoir rock and the nature of the fluids they contain. Porosity and fluid saturation are of most interest in the estimation of hydrocarbon reserves. Permeability, which is a measure of the resistance of the formation to fluid flow, including invasion, and which is needed to estimate recoverable reserves, has been the most difficult parameter to determine and quantify by borehole surveys.

One of the primary difficulties in formation analysis from borehole surveys is the need to determine and compensate for the effects of invasion. Invasion takes place in porous permeable zones where the hydrostatic/dymanic pressure of the drilling mud is greater than the formation pore pressure. The invasion of the mud filtrate will cause a radial variation of the formation resistivity. Since one objective of most logging methods is to measure the properties of the undisturbed formation, logging tools are usually designed to be insensitive to invasion or compensated for this fluid replacement by means of departure curves. In the prior art, this is normally done by running wireline resistivity logs which have different response functions, or depths of investigation (i.e., spherically focused (or short-normal) logs, medium induction logs, and deep induction logs).From these measurements, invasion corrections can be derived. (See, for example, Dresser-Atlas, Log Interpretation Charts, (1979)). More recently, a different approach has been used, which is to solve the "forward" problem. In this method, one assumes a lithology, porosity and water saturation (among other parameters) and varies these until a fit is obtained to the suite of wireline logs. (See, for example the following papers at the SPWLA Twenty-fifth Annual Logging Symposium, New Orleans, June 10-13, 1984. M.A. Yuratich and Walter J. Meger, "The Application of Finite Difference Methods to Normal Resistivity Logs," (Paper V). Mark Alberty and Khaled H. Hashmy, "Application of ULTRA to Log Analysis", (Paper Z).Yih-Yih Lin, Stan Gianzero, and Robert Strickland, "Inversion of Induction Logging Data Using the Least Squares Technique" (Paper AA).) However, these wireline methods have serious drawbacks or shortcomings. In either analysis one is looking at "snapshots" of the formation at a particular time (typically many days after drilling), with each measurement having a different response function. These analyses give no information about the time dependence, i.e., the rate, of the invasion process, but are designed only to reconstruct the formation's properties prior to invasion.

These previous methods involve taking a suite of several wireline measurements days or weeks after drilling has occured at a particular depth in a borehole. Because these wireline logs are taken so long after drilling, invasion is often essentially complete and these logs provide no information about the rate of invasion, and, hence, do not provide information about permeability.

Known methods to obtain permeability information include laboratory analyses of core samples and inferences from other suites of wireline logs. These prior methods are cumbersome, time consuming and/or highly unreliable. Indeed, it is often considered to be doing well to get within an order of magnitude of the actual permeability.

Viewed from one broad aspect there is herein disclosed a method of determining a permeability parameter of a borehole, including the steps of: taking a first resistivity measurement of a zone of a formation within a short and known time after drilling of the zone has occurred; taking at least one resistivity measurement of said zone subsequent to said first measurement to obtain incremental data of the rate of invasion of said zone by drilling fluid; and determining a permeability parameter of said zone from said resistivity measurements.

The method disclosed herein, at least in its preferred forms, employs a series of resistivity logs obtained in a measurement while drilling (MWD) system to generate information relating to formation invasion and formation permeability.

The theory of resistivity logging is based on the concept that dry rock is an insulator, and that, therefore, any resistivity in a borehole formation is related to the quantity and type of fluid in the pores of the formation. Thus, for example, since salt water is highly conductive it will cause low resistivity readings; whereas oil, which is a poor conductor, will cause much higher resistivity readings. Since the drilling mud is at a higher pressure than the pressure of the fluids in the formation, it tends to displace the connate fluids in the formation. This distorts or changes the resistivity readings.

The method of a preferred embodiment makes use of the displacemnent of connate fluids by the drilling mud to generate data by resistivity logging which is commensurate with the rate of formation invasion and formation permeability. In this embodiment, first resistivity logs are taken in various zones within minutes after the rock of each zone has been drilled. Subsequent resistivity logs of the zones; i.e., relogs, are taken at intervals when the drill string is being removed from or lowered into the borehole. The data thus generated is then processed to obtain results which are commensurate with formation invasion and - permeability.

Since a measurement-while drilling resistivity logging tool is used, formation parameters may be measured within minutes after the bit has drilled the rock. Under these conditions, the effects of invasion will be minimal for the first log. Furthermore, since it is feasible to relog sections of particular interest several times during later drilling operations before the zone becomes flushed, i.e., before equilibrium is established between the pressure of the mud filtrate and the connate fluid, this relogging technique provides what amounts to a "motion picture" of the invasion process. Since all measurements are made with the same tool, it is not necessary to account for different tool response profiles; any difference observed can be directly attributed to the results of the drilling process, i.e. invasion and/or formation damage.

The logging and relogging aspect of the method of the preferred embodiment generates data which strikingly illustrates the invasion of drilling fluid into the formation as it occurs; and the analytical aspect of the method produces results which provide indications of formation invasion and permeability.

Some embodiments of the present invention will now be described by way of example with reference to the accompanying drawings in which: Figure 1 is a schematic representation of an MWD system, Figures 2 and 3 are charts showing fluid invasion profiles for a first borehole obtained using conventional wireline techniques and a MWD logging and relogging technique, Figure 4 shows a plot of the parameter F(t) (invasion factor) vs. time, Figure 5 is a chart showing fluid invasion profiles for a second borehole obtained using MWD logging and relogging, Figure 6 is a plot of the parameter IR(t)-R(oc)I vs. time, and Figure 7 is a plot of the slope of resistivity curves of the second borehole over the depth zone of interest in Fig. 5.

Referring to Fig. 1, apparatus for MWD resistivity logging is schematically illustrated. The apparatus is a resistivity-gamma-directional (RGD) MWD measurement tool presentiy in commercial use by Teleco Oilfield Services Inc. The apparatus shown in Fig. 1 is a drill coilar having an upper steel section 10 and a lower non-magnetic sub 12, sections 10 and 12 being threaded together, with appropriate electrical interconnections.

As mentioned, the RGD tool of Fig. 1 consists of two sections 10 and 12 and has an overall length of 39 feet. The upper portion 10 of the tool houses Teleco's standard directional sensor 14, turbine/alternator 16, transmitter 18 and microprocessor and electronics 20. The downhole microprocessor 20 processes the logging data and controls the transmission sequence: continuous formation gamma and resistivity data while rotating, directional data while not rotating. The transmitter 18 is a mud pulse transmitter as shown and described in U.S. Patents Nos.

3,982,431; 4,013,945; and 4,021,774.

An auxiliary sensor module is mounted in sub 12. Located in the center of this unit is the gamma ray detector 22 and an electronics package 24. The collar surface has an insulated region 26 with two electrodes 28 and 30 which operates as a 16 inch short normal resistivity tool. Electrical connection of the resistivity electrodes to electronics 24 may be accomplished in any standard manner known in the art. The geometry of the tool places the measure points of the gamma ray and resistivity sensors at the same depth. This facilitates analysis and allows both measurements to be made in a zone of interest at the same time.

The 16 inch short normal MWD resistivity measurement may, in accordance with conventional practices, be used to: o Correlate with offset well logs o Determine lithology o Define thin beds and bed thickness o Identify potential hydrocarbon-bearing intervals o Determine pore pressure trends The method of this preferred embodiment adds another and very important capability to this device. In accordance with this embodiment, a first MWD resistivity log is taken in a borehole in real time and within minutes after a zone is drilled. Thereafter several relogs are taken of one or more zones of interest (which may be determined from the MWD and offset logs). The relogs are taken during subsequent trips of the drill string out of and back down the well or during short trips. The rate of descent of the drill during relogs is controlled to provide the desired bed resolution.

Before further discussion of the analytical method, a brief discussion is in order of actual relogging examples performed in wells to test and prove the method.

Real-time MWD gamma and resistivity logs were recorded on two Gulf Coast wells ("Well A" and "Well B"). Both wells are onshore and drilled delta plain sediments consisting of clastic sequences of sands and shales. Structurally, the wells were located close to down-to-the-basin faults commonly found throughout the Gulf of Mexico region.

On the basis of MWD and offset logs, zones were selected from each well for further study.

During subsequent trips or short trips, the zones of interest were relogged several times. The rate of descent of the drill string during these relogs was held to approximately 60 feet/hour, resulting in a bed resolution of approximately one foot. While this gave a somewhat less detailed log than was obtained in real time, it was adequate for this analysis.

The selected zones of Well A were logged up to 7 times over a period of 12 days; those of Well B up to 6 times in 7 days. The times are summarized in Table 1. For comparison purposes, wireline logs were run approximately two weeks after the initial drilling of the selected zones.

TABLE 1 Summry of MWD Relog Times Well A Zone 4955-5065 feet 6945-7025 feet 7935-8015 feet Times 0* - 0* (hrs). 3 3 37 33 8 62 76 39 110 104 175 136 300 Logs taken while drilling; generally within 30 min.

Well B Zone 6959-7010 feet 7180-7220 feet Times 0* 0* (hrs.) 3 4 28 17 39 73 105 145 136 Logs taken while drilling; generally within 30 min.

An underlying theory of the method is that MWD relogs of a well will provide time dependent data on fluid invasion. The data obtained from Well A proves that theory. In Fig. 2, the data to the right of the depth stations shows resistivity (in ohm-m) for an expanded view of MWD relogs for the depth zone between 6950 and 7120 feet for a series of logs taken over a time period from while drilling to 300 hours after drilling. The data to the left of the depth stations shows spontaneous potential (in millivolts) for a wireline reading of the same zone. The resistivity of the drilling mud (Rm) ws approximately 1 ohm-m. In the MWD resistivity graph, the locale where resistivity readings are essentially constant (i.e., above 6960 feet and below 7000 feet) indicates the presence of shale, while the zone between 6960 and 7000 feet where resistivity varies indicates a sand zone.The wireline readings confirm the presence of the sand and shale zones.

Several features stand out prominently in the MWD resistivity data of Fig. 2. The shale resistivities are essentially constant over time. There is a clear evolution in time in the sand zone resistivity. (Note: the formation and mud resistivities in this example were such that the borehole corrections were negligible.) The fact that the resistivity increases with time reflects the freshness of the mud. The mud filtrate resistivity is higher than that of the connate fluids. The properties of the mud and the differential pressure remained approximately the same during the relogs. The resistivity changes rapidly in the period just following drilling, then more slowly.The similarity of the last passes reflects the fact that once the sensor's zone of investigation is essentially flushed or pressure equilibrium has been reached, there will be little further change in the apparent resistivity.

The MWD data in Fig. 2 confirm the theory that MWD relogs provide incremental data, and hence data about the rate of invasion, before a zone is flushed or pressure equilibrium is reached. This important data has not heretofore been available.

Fig. 3 compares the first (real time, while drilling) and last (300 hours) MWD runs with two wireline resistivity logs, deep induction and spherically focused, run approximately five days after the last MWD relog. This comparison reveals some further interesting and unexpected results.

The real-time MWD resistivity in the sand zone is lower than even the deep induction wireline log. This implies that when using a short normal log in an MWD application, the readings may be less influenced by invasion than those of a deep induction tool, and therefore read close to the true resistivity R,. The first MWD resistivity was measured at about 5 minutes after drilling, during which time the formation was exposed to invasion. Where a longer time has elapsed, a simple "curve extrapolation" of MWD resistivity back to t=O may give a better estimate of R,.

Also, the final MWD pass has a higher peak value than the spherically-focused wireline log. This indicates that the short range characteristics of the MWD normal devices cause it to read closer to the flushed zone resistivity, R,,. Thus, a single sensor, when used in an MWD relog mode can generate the range of apparent resistivity measurements from Rt to R,, as a function of time.

Having established, from the data of Fig. 2, that MWD relogging generates time dependent information about invasion, consideration can be given to the steps for analyzing that information to obtain a reliable permeability parameter.

The apparent resistivity, as seen by the MWD sensor, will change with respect to time if: o the zone being measured is porous and permeable, o the invasion fluid (filtrate) resistivity differs from that of the formation fluid, and o the borehole pressure exceeds the formation pressure.

To simplify the basis for evaluation of formation invasion, it is assumed that the invasion is of a single-phase fluid of uniform viscosity flowing radially outward from the borehole into a porous and permeable formation. Under these conditions, Darcy's relationship for single phase fluid flow states that the flow rate will be proportional to the formation permeability and the pressure differential and inversely proportional to the fluid viscosities.

In looking at the sequence of relogs in Fig. 2, as a function of time, it can be seen that there is a smooth transition from a reading that is close to R, to one that will approximate R,,, if sufficient time has elapsed (i.e., enough time for the invasion process to be completed within the range of the sensor). ("Sufficient" time will be discussed in more detail below). Since the same sensor is used for all the measurements and the borehole conditions have remained about the same, the changes in formation resistivity must result from the movement of mud filtrate from the borehole into the porous and permeable formation. In the absence of formation damage, the principal change will be the replacement of the connate water with mud filtrate.

In attempting to model this process, it is assumed that the formation resistivity will vary between two limits, R(to) (unperturbed) and R(r) (the final measurement), which will depend on the fluid resistivities of the connate fluid (Rw) and the mud filtrate (rum) and the formation porosity (O, It is necessary to account for the differences in these endpoint parameters in order to measure a rate of invasion. This rate will be proportional to the difference between the borehole pressure (Pbh) and the formation pressure (P,), proportional to the permeability (k) and inversely related to the mud filtrate and connate fluid viscosities. Of these, it is assumed that only the pressure difference may undergo significant change with time.

It is further assumed that the initial invasion by mud filtrate does not completely displace the connate fluid in all the pores. Once the invasion front passes a radial zone, the filtrate will then diffuse by tortuous paths into the surrounding pores. Under these conditions, after an initial abrupt change in resistivity, the transition should be approximately exponential. This assumption ignores the fact that a diffusion equation solution in cylindrical coordinates is not a pure exponential, and the fact that the apparent resistivity is a weighted integral over the range of the sensor. It is expected, however, that the exponential behaviour will dominate, and the test data verify this conclusion.

On the basis of the previous assumption, we can write the apparent formation resistivity (after borehole correction) at any time t, R(t), as a function of the initial and final measured resistivities, R(to) and R(x), as follows: R(t)-R()--(R(r0)-R()) exp[--A f(t)] (1) (where A is a constant related to fluid viscosities) According to the assumptions above, f(t), the time-dependent part of the equation, will depend on the permeability, the pressure differential and the viscosities of the fluids (assumed constant). It will be a time integrated function to reflect the fact that the response at any time is a weighted volume integral.An invasion fraction, F(t) can be defined as follows:

where A' is a proportionality constant related inversely to the viscosities of the drilling fluid and the connate fluid, and k is the permeability of the formation. If the borehole pressure, Pbh, is assumed constant, equation (2) can be rewritten as R(t)R(oo) F(t)= -A" k '(PbhPf)(tto) (3) R(t,)-R(oo) Thus, if a plot were made of F(t) vs. time on a semilogarithmic scale, the slope, after an initial surge, would be proportional to the permeability. Furthermore, since the effects of porosity and connate fluid resistivity are normalized out on the left hand side of equation (3), all formations of the same permeability should be on the same line.An alternate method is to plot IR(t)--R(oc)l vs. time. In this case, formations of equal permeability will lie on parallel lines.

In order to plot these curves, it is necessary to have a value of R(oo), that is, it is necessary to allow sufficient time for the invasion process to be complete within the range of the sensor.

"Sufficient" time can be defined as that time after which repeated logs show no further change; the last reading will correspond to R(oo). This will cause a slight problem because a value of zero cannot be plotted on a semilogarithmic scale. By choosing a value for R() slightly beyond the last point, this problem can be avoided.

In cases where the resistivity is changing even at the last relog, and if a reliable estimate of Rw can be obtained from other data, R(cc) can be estimated from the equation R(m)/R(t,)= R,,/R, (4) It can, thus, be seen from equations (1)-(4) and the related discussion, that a method is provided for determining a permeability parameter of the formation under investigation. This permeability parameter is easily obtained by MWD resistivity relogging; and it is accurate, reliable, and easy to derive. It should be noted that the method derives a permeability parameter which is not an absolute value of porosity, but is, rather, a comparative or relative measure of permeability. Nevertheless, the permeability parameter provides valuable information about the formation.However, if the viscosity and pressure of the connate fluid are known (as from wells in the area), and since the drilling fluid viscosity and pressure are known, the absolute value of permeability can be determined from the foregoing equations and analyses.

Further applications of the method of the present embodiment can be seen by referring to Fig.

4, which is an idealized plot of equation (3). (A plot of a change in borehole pressure vs. time is also shown in Fig. 4). In certain cases, there may be a significant change in the borehole pressure, as when mud weight is increased. According to equation (3), this will cause a change in slope of the curve. If P1 and P2 are the (constant) bottomhole pressures during the two time intervals, and m, and m2 are the corresponding slopes of the curves, then: m2/ml=(P2P,)/(PrP,); or (5) P,=(m2P1mlP2)/(m2m1) (6) So, when mud weight is changed, the method can give an estimate of the formation pressure.

The estimate will be most reliable when the slope change is great, i.e., when the original conditions were close to balance, or the weight change is significant when compared to the differential.

Furthermore, once the relog resistivities cease to change, Archie's equation may be used to estimate the effective porosity when Rmf is known or can be estimated.

Referring now to Figs. 5, 6 and 7, another example is presented.

Fig. 5 shows the relogs of one zone from Well B. As in the data from Well A, there are sections of the log which show extensive changes in resistivity (corresponding to a sand zone), and some which show almost none (corresponding to shale zones). In this case, since the sand zone resistivities are still changing even at the last relog, R(x) was estimated from equation (4), using Rum,=1.05 and R,=0.05.

Fig. 6 shows plots of IR(t)--R(x)l against time for several points from Well B including some from a shallower zone. Time was used as the base of this curve because mud weights were held constant. All of the curves show the expected exponential behaviour, but with different slopes. All of the curves also show a much steeper slope for the first few hours, which is consistent with the assumption of the diffusive nature of invasion.

A plot of the slopes of the resistivity curves as a function of depth over the zone of interest is shown in Fig. 7. To eliminate the fast initial invasion, the slope is taken beginning with the second relog. As expexted from looking at the raw data (Fig. 5) there is a wide range in the slopes, with permeable sands separated by impermeable stringers. Note that although the resistivity change at 7207', is greater than that at 7200', the latter zone appears to be more permeable; the difference is apparently a result of slightly higher shale content at 7200'.

In the absence of any quantitative data on the formation pressure, an estimate of the actual formation permeabilities cannot be derived. According to the method of the present embodiment, however, values in Fig. 7 are proportional to the permeabilities.

As can readily be understood by those skilled in the art from the foregoing discussion, the method disclosed herein presents a simple way to achieve a reliable permeability parameter of desired formation zones through the use of an MWD relogging technique and processing of data thus obtained. The method makes it possible; for the first time; to obtain reliable data about the rate of invasion and formation permeability, and thus makes an important contribution to the art of borehole logging.

The importance of the method disclosed above is not in the fact of relogging, per se, in an MWD environment, but in applying the concept of relogging in an MWD environment to a method whereby accurate and reliable permeability information can be generated. It should be noted that the method provides what may be considered to be an index of permeability. This is, the method provides a comparative or relative indication of permeability rather than an absolute permeability value. This relative information is of great value in evaluating borehole formations, especially since the information is obtained while drilling and does not rely on laboratory analyses of cores.

Modifications to the specific embodiments and to any broad aspects thereof referred to or suggested herein may be apparent to those skilled in the art and the disclosure hereof is intended to encompass any such modifications.

Claims (11)

1. A method of determining a permeability parameter of a borehole, including the steps of: taking a first resistivity measurement of a zone of a formation within a short and known time after drilling of the zone has occurred; taking at least one resistivity measurement of said zone subsequent to said first measurement to obtain incremental data of the rate of invasion of said zone by drilling fluid; and determining a permeability parameter of said zone from said resistivity measurements.

2. A method according to claim 1 further including: taking a plurality of resistivity measurements of said zone subsequent to said first measurement, at least some of said subsequent measurements being taken prior to the zone being flushed.

3. A method according to claim 2 wherein: said first resistivity measurement is taken prior to any significant invasion of said zone by said drilling fluid.

4. A method according to claim 2 or 3 wherein: at least one of said resistivity measurements is taken after the resistivity in said zone ceases to change due to effects of invasion of the drilling fluid.

5. A method according to any preceding claim wherein the step of determining a permeability parameter includes determining an invasion factor F(t) where: R(t)-R(x) F(t) = R(to)R(x) where R(t) is formation resistivity at time t, R(to) is the initial or unperturbed formation resistivity, and R(x) is the final or flushed formation resistivity.

6. A method according to claim 5, including the step of plotting F(t) versus time.

7. A method according to claim 6 wherein F(t) is plotted versus time on a semilogarithmic scale to determine the permeability parameter from the slope of the plot.

8. A method according to claim 7 wherein formation pressure is determined from changes in the slope of said plot which are commensurate with changes in the borehole pressure.

9. A method according to claim 5 wherein an absolute value of permeability is determined from said permeability parameters and from the viscosities and pressures of the drilling fluid and the formation fluid.

10. A method according to any preceding claim wherein said first resistivity measurement is taken while drilling the borehole, and said subsequent resistivity measurements are taken at known times after said first measurement.

11. A method of determining a permeability parameter of a borehole substantially as hereinbefore described with reference to the accompanying drawings.