Patent Publication Number: US-7222022-B2

Title: Method of determining properties relating to an underbalanced well

Description:
FIELD OF THE INVENTION 
   The invention relates a method of determining properties relating to an underbalanced well, and in particular deriving properties such as pore pressure, permeability and porosity for fluid-producing formations contributing to fluid output from the well. 
   BACKGROUND TO THE INVENTION 
   Boreholes are sometimes drilled using drilling fluid which has a pressure substantially less than the pressure of the fluid from the formation. This is known as underbalanced drilling. Underbalanced drilling is often used where fluid-bearing formations are known to be delicate and prone to damage, so as to maintain the integrity of the formation. Typically a number of different subterranean structures, with different properties, are drilled through before the actual production formation of interest is reached. The pressure of fluid from the formation will often therefore vary during drilling. It is often important to ensure that the drilling remains underbalanced at all times to minimise formation damage. 
   Underbalanced drilling is also used generally where, for example, faster drill speeds are required or where the life of a drill bit needs to be extended. 
   The formations surrounding the borehole can be characterised by a pore pressure, porosity and permeability. When underbalanced drilling, an estimate of the pore pressure is typically made, and the pressure of the drilling fluid is then chosen in an attempt to ensure that underbalanced drilling is achieved at all times. However, the estimates of pore pressure are generally very inaccurate and as such it is often difficult to perform underbalanced drilling with any degree of reliability or control. 
   The estimate of pore pressure can be used to derive the permeability of the formations, but the estimated pore pressure can be very inaccurate so causing errors in the values of permeability. 
   The present invention aims to provide a method which supplies more information about formations whilst drilling and aims to enable more controlled underbalanced drilling to be achieved. 
   SUMMARY OF THE INVENTION 
   In accordance with the present invention, there is provided a method of calculating properties relating to a subterranean formation. The method comprises drilling a borehole into the subterranean formation; measuring a first pressure in the borehole and a first fluid flow rate when the drilling has progressed to the first location; measuring a second pressure in the borehole and a second fluid flow rate when the drilling has progressed to the second location; and calculating a property of at least a portion of the formation using the first and second pressures and the first and second fluid flow rates. 
   The fluid flow rates can be measured as the fluid exits the borehole at or near the surface, or may be measured in close proximity to the drill bit. The pressure measurements are preferably of annular bottomhole pressures. The step of calculating preferably comprises calculating at least two of the following types of properties: pore pressure, porosity, and permeability. 
   The method can also preferably include a third set of measurements taken when the drilling has progressed to a third location in the formation, and the step of calculating further comprises calculating all three of the following properties: pore pressure, porosity, and permeability. 
   Variations in pressure in the borehole can be induced by various methods, including one or more of the following: altering the flow rate of drilling fluid; placing a tool in the borehole which emits acoustic pulses into fluid within the well; altering the density of drilling fluid used; use of a choke unit. However, variations in pressure can also be caused by unintentional variations in pumping of drilling fluid. 
   The step of calculating preferably comprises using a first relationship between the first flow rate and the first pressure, a second relationship between the second flow rate and the second pressure, and solving the first and second relationships to obtain a value for the property, with the first and second relationships preferably expressing the measured flow rates as a function of drawdown, rate of penetration, and the rate response of a portion of the formation. 
   Advantageously, the step of drilling is preferably not interrupted during the measurement steps. 
   The present invention is also embodied in a system for calculating properties relating to a subterranean formation, comprising: a pressure sensor configured to measure pressures in a borehole in the formation in close proximity to a drill bit used to drill the borehole; a flow sensor configured to measure flow rates of fluid flowing through the borehole; and a processor adapted to calculate a property of at least a portion of the formation using first and second measured pressures and first and second fluid flow rates, wherein the first and second pressures are measured by the pressure sensor when the drilling has progressed to a first and second location respectively, and the first and second flow rates are measured by the flow sensor when the drilling has progressed to a first and second location respectively. 
   As used herein, the term “induce” when referring to pressure variations includes both intentional and unintentional changes in pressure. For example, the induced pressure changes can be caused by uncontrolled variations in the drilling fluid pumping speeds, or other “noise” in the form of unintentional pressure variations induced by the drilling process. 
   The pressure variations cause a change in flow rate from formations along a length of borehole, and as such a change in the production flow rate of the well. The variations in pressure can be used to calculate the pore pressure. 
   The variation in annular bottomhole pressure causes changes in the flow rate from the formations and as such the production flow rate, i.e. the total output of the well, changes. The variations in production flow rate allow analysis of the profile of the formations along the length of the borehole. By monitoring pressure variations over a small distance of drilled borehole over which distance one can assume that properties of the reservoir remain constant, and by correlating the changes in production flow rate with pressure variations, pore pressure for formations over the given length can be determined. The method thus avoids the need to estimate pore pressure, and instead provides a way of calculating a true pore pressure much more accurately. As drilling proceeds, monitoring of the pressure variations continues and with the changes in production flow rate, a profile of pore pressure along the length of the borehole is obtained. Thus, if desired, real time measurements conducted whilst drilling can be used to create a real time profile of the pore pressure and, where desired, also porosity. Permeability may also be derived. 
   By obtaining a profile of pore pressure along the length of the borehole, and not using an estimate or assumption of the pore pressure, the properties of the formation are known with a great deal of resolution. 
   By using a real time profile, a number of advantages are achieved in that the pore pressure of the well is constantly monitored as drilling occurs. Typically formations where underbalanced drilling is required have a pressure of around 10 MPa and thus the induced pressure variations are preferably kept in the range 2 MPa–5 MPa so as to ensure that the drilling is kept underbalanced. Thus it can be guaranteed that drilling is underbalanced at all times. 
   Permeability steering can also be undertaken, and as the properties of the formation are known along its length, the need for testing the well after drilling, and the need to shut down the well, when testing, can be avoided. 
   Productivity steering is also possible where the well is redesigned during drilling based on the measurements obtained, so maximising productivity. The method in accordance with the invention may also be used in a variety of other well operations, including during completion of a well and for targeted, or intelligent, perforating of the well. 
   The invention is also of advantage in that permeable zones and damaged zones are identified with a great degree of accuracy, and as such it is simpler to identify where casings and cement need to be perforated when completing the well. The invention also allows benchmark testing as drilling occurs. 
   In accordance with another aspect of the invention, there is provided apparatus for performing the above described method. 

   
     BRIEF DESCRIPTION OF THE DRAWINGS 
     The invention will now be described by way of example and with reference to the accompanying drawings in which: 
       FIG. 1  shows a schematic view of a section through a borehole, and used for explanation; 
       FIG. 2  shows a graph of pressure variation with time; 
       FIG. 3  shows a series of graphs relating to pressure variation downhole and subsequent calculation of properties of formations using a method in accordance with the present invention; 
       FIG. 4  illustrates a series of graphs showing properties obtained using a prior art method; 
       FIG. 5  shows a system for calculating properties relating to a subterranean formation, according to a preferred embodiment of the invention; and 
       FIG. 6  shows steps involved in calculating properties relating to a subterranean formation, according to a preferred embodiment of the invention. 
   

   DETAILED DESCRIPTION OF THE INVENTION 
     FIG. 1  illustrates a borehole  10  which has been artificially divided along its length into a series of layers z 1 ,z 2 ,z 3  . . . z n . As long as drilling occurs at a suitable rate, one can assume that the change in properties of the formations along the borehole are negligible over small distances. By artificially dividing the borehole along its length into formation layers at successive distances of z 1 ,z 2 ,z 3  . . . z n , one can assume for layer z 1 , drilled at annular bottomhole pressure BHP 1  at a time t 1  and formation z 2 , drilled at annular bottomhole pressure BHP 2  and at time t 2 , that the properties of the two formations z 1 ,z 2  are constant and that any changes in the production flow rate or flux from the well are as a result of the change in pressure. This allows one to solve two simultaneous equations relating to the properties of these artificial layers to deduce the pore pressure of fluid in the reservoir for layers z 1  and Z 2  and also the permeability of the permeable rock forming layers z 1  and Z 2 . 
   By measuring a further pressure change for layer Z 3 , three simultaneous equations are arrived at and these can be solved for the three variables of pore pressure, i.e. pressure of fluid in the reservoir, permeability and porosity. Permeability and porosity are both properties of the permeable rock. By conducting an analysis in this way, and sub-dividing the length of the borehole into a series of impedances, an accurate profile of the true properties of the formations along the length of the borehole is obtained. 
   The method described derives formation pressure for underbalanced drilling of a reservoir where the well is drilled such that the pressure within the wellbore is below the formation pressure of the reservoir. During drilling, a continual influx of formation fluid occurs into the wellbore which results in changes in the production rate of the well as the borehole is drilled. By measuring variations in the annular bottomhole pressure during drilling, the local pore pressure along the length of the borehole can be calculated. According to a preferred embodiment pressure changes can be intentionally induced downhole, such as for a gas reservoir by pulsing the liquid injection rate at surface, and for a liquid reservoir by varying the gas injection rate, and monitoring the changes in pressure, so that the local pore pressure can be derived from the transient response of the well. 
   The variations in pressure can be sinusoidal in nature, as shown in  FIG. 2  which illustrates the pressure downhole P w  as a function of time. Stepped changes in pressure can also be used for analysis by the method, and in general, any type of variation in pressure can be used to practice the present invention so long as the rate of change is sufficient for the resolution required. As shown in  FIG. 2 , if the time Δt between annular bottomhole pressure  1  (BHP 1 ) and annular bottomhole pressure  2  (BHP 2 ) is 1 hour, and the drilling rate is 2 m an hour, then the chosen depth of formation z 1  is 1 m, and formation Z 2  is also 1 m. The more accurately the production flow rate of the well can be measured, the less change in pressure is needed to achieve suitable data for analysis. 
   The resolution of the pore pressure profile obtained will depend on the resolution and accuracy of the measurements of annular bottomhole pressure and production flow rate made whilst drilling, in conjunction with the rate of penetration (ROP) of the drill bit. For a low ROP compared to the sampling rate of the pressure and flow rate data, the spatial resolution of the pore pressure profile will be high. 
   According to a preferred embodiment, pore pressure and permeability profiles are simultaneously derived whilst drilling. Thus the assumption of a fixed pore pressure for the length of the borehole can be dispensed with. Furthermore, direct measurement of local drawdown, i.e. (bottomhole pressure−surface pressure)/(pore pressure−wellbore pressure), enhances the accuracy of the derived permeability values. 
   When analysing the data using artificial layers as shown in  FIG. 1 , the mass flow rate Q from a segment or section of the reservoir at position z 1 , drilled at time t 1 , is described by
 
 Q ( z   1   ,t   1 )=Δ p ( z   1   ,t   1 )ρ U ( z   1   ,t   1 )Δ t   1   ∫Hdt   (1)
 
where Δp(z 1 ,t 1 )is the local drawdown, ρ is the density of the produced fluid, Δt 1  is the duration of drilling the zone, U(z 1 ,t 1 )is the rate of penetration (ROP) of the drill bit (which is a function of time but assumed constant over Δt 1 ), and H is the local rate response of the reservoir to the drawdown, which depends on the local reservoir characteristics.
 
   The mass flow rates dQ from a segment of the reservoir at position z 1 , drilled over time Δt 1 , is described by
 
 dQ ( z   1   ,t   1 )=Δ p ( z   1   t   1 )ρ U ( z   1   t   1 )Δ t   1   ∫Hdt   (2)
 
where Δp(z 1 ,t 1 ) is the local drawdown, equal to [P f (z 1 ,t 1 )−P BHP (z 1 ,t 1 ) 2 ], ρ is the density of the produced fluid, U(z 1 ,t 1 ) is the rate-of-penetration (ROP) which is a function of time but assumed constant over Δt 1 , and H is the local rate response of the reservoir to the drawdown, which depends on the local reservoir characteristics. P f  is the formation pressure and P BHP  is the bottomhole pressure.
 
   For a known reservoir pore pressure and porosity profile, this may be used to derive the formation permeability profile. 
   For a second layer at z 2 , adjacent to the layer at z 1 , but drilled at bottomhole pressure ΔP 2 =P f (Z 2 ,t 2 )−P BHP (Z 2 ,t 2 ), where P BHP (Z 2 ,t 2 )≠P BHP (z 1 ,t 1 ), the production rate will change according to
 
 dQ ( z   2   ,t   2 )=Δ p ( z   2   ,t   2 )ρ U ( z   2   ,t   2 )Δ t   2   ∫Hdt   (3)
 
   For a sufficiently high data acquisition rate compared to the ROP and the heterogeneity of the reservoir itself, we assume that the reservoir characteristics do not vary significantly between z 1  and Z 2 , and that therefore the rate response, H, is constant between the two, i.e. we assume
 
 P   f ( z   1   ,t   1 )= P   f ( z   2   ,t   2 )  (4)
 
and the permeability κ of z 1  and z 2  is such that
 
κ( z   1   ,t   1 )=κ( z   2   ,t   2 )  (5)
 
and the porosity φ is
 
φ( z   1   ,t   1 )=φ( z   2   ,t   2 )  (6)
 
   The volumetric flow rate from any segment k at bottom hole immediately after drilling may be written 
   
     
       
         
           
             
               
                 
                   dq 
                   k 
                 
                 = 
                 
                   Δ 
                   ⁢ 
                   
                       
                   
                   ⁢ 
                   
                     p 
                     k 
                   
                   ⁢ 
                   
                     
                       4 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       π 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         r 
                         w 
                       
                       ⁢ 
                       
                         κ 
                         k 
                       
                     
                     μ 
                   
                   ⁢ 
                   
                     
                       Δ 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         z 
                         k 
                       
                     
                     
                       
                         log 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         
                           ( 
                           
                             Γ 
                             k 
                           
                           ) 
                         
                       
                       - 
                       γ 
                     
                   
                 
               
             
             
               
                 ( 
                 7 
                 ) 
               
             
           
         
       
     
   
   where r w  is the radius of the wellbore, μ is the viscosity of the produced fluid, γ is Euler&#39;s constant which equals 1.78, and 
   
     
       
         
           
             
               
                 
                   Γ 
                   k 
                 
                 = 
                 
                   
                     4 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     
                       κ 
                       k 
                     
                     ⁢ 
                     Δ 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     
                       t 
                       k 
                     
                   
                   
                     ϕ 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     μ 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     
                       c 
                       t 
                     
                     ⁢ 
                     
                       r 
                       w 
                       2 
                     
                   
                 
               
             
             
               
                 ( 
                 8 
                 ) 
               
             
           
         
       
     
   
   for Δt k  the duration of drilling for segment k, where c t  is the compressibility of fluid flowing in the reservoir. 
   Note that the flow rate from any zone at bottomhole might be measured at bottomhole, or estimated from surface. 
   So the volumetric flow rate from the reservoir for segments  1  and  2  immediately on being drilled over equal timescale Δt are 
                     dq   1     ⁡     (     t   1     )       =       (       P   f     -       P   BHP     ⁡     (       t   1     ,     z   1       )         )     ⁢       4   ⁢           ⁢   π   ⁢           ⁢   κ   ⁢           ⁢     r   w       μ     ⁢       Δ   ⁢           ⁢     z   1           log   ⁡     (       4   ⁢           ⁢   κ   ⁢           ⁢   Δ   ⁢           ⁢   t       ϕ   ⁢           ⁢   μ   ⁢           ⁢     c   t     ⁢     r   w   2         )       -   γ                 (   9   )             and                             dq   2     ⁡     (     t   2     )       =       (       P   f     -       P   BHP     ⁡     (       t   2     ,     z   2       )         )     ⁢       4   ⁢           ⁢   π   ⁢           ⁢   κ   ⁢           ⁢     r   w       μ     ⁢       Δ   ⁢           ⁢     z   2           log   ⁡     (       4   ⁢           ⁢   κ   ⁢           ⁢   Δ   ⁢           ⁢   t       ϕ   ⁢           ⁢   μ   ⁢           ⁢     c   t     ⁢     r   w   2         )       -   γ                 (   10   )               
where the permeability κ, pore pressure P f  and porosity φ of the rock are unknown but the same in each equation. Terms involving the permeability and porosity of the rock are eliminated from these two equations to yield the local pore pressure, P f (z), applicable to any two such zones, N-1 and N, drilled at constant time interval (for the expression given) but of varying thickness Δz N-1  and Δz N . This is written
 
   
     
       
         
           
             
               
                 
                   
                     P 
                     f 
                   
                   ⁡ 
                   
                     ( 
                     z 
                     ) 
                   
                 
                 = 
                 
                   
                     ( 
                     
                       
                         
                           
                             dq 
                             
                               N 
                               - 
                               1 
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               t 
                               
                                 N 
                                 - 
                                 1 
                               
                             
                             ) 
                           
                         
                         ⁢ 
                         
                           
                             
                               P 
                               BHP 
                             
                             ⁡ 
                             
                               ( 
                               
                                 t 
                                 N 
                               
                               ) 
                             
                           
                           ⁡ 
                           
                             [ 
                             
                               
                                 Δ 
                                 ⁢ 
                                 
                                     
                                 
                                 ⁢ 
                                 
                                   z 
                                   N 
                                 
                               
                               
                                 Δ 
                                 ⁢ 
                                 
                                     
                                 
                                 ⁢ 
                                 
                                   z 
                                   
                                     N 
                                     - 
                                     1 
                                   
                                 
                               
                             
                             ] 
                           
                         
                       
                       - 
                       
                         
                           
                             dq 
                             N 
                           
                           ⁡ 
                           
                             ( 
                             
                               t 
                               N 
                             
                             ) 
                           
                         
                         ⁢ 
                         
                           
                             P 
                             BHP 
                           
                           ⁡ 
                           
                             ( 
                             
                               t 
                               
                                 N 
                                 - 
                                 1 
                               
                             
                             ) 
                           
                         
                       
                     
                     ) 
                   
                   
                     ( 
                     
                       
                         
                           
                             dq 
                             
                               N 
                               - 
                               1 
                             
                           
                           ⁡ 
                           
                             ( 
                             
                               t 
                               
                                 N 
                                 - 
                                 1 
                               
                             
                             ) 
                           
                         
                         ⁡ 
                         
                           [ 
                           
                             
                               Δ 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               
                                 z 
                                 N 
                               
                             
                             
                               Δ 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               
                                 z 
                                 
                                   N 
                                   - 
                                   1 
                                 
                               
                             
                           
                           ] 
                         
                       
                       - 
                       
                         
                           dq 
                           N 
                         
                         ⁡ 
                         
                           ( 
                           
                             t 
                             N 
                           
                           ) 
                         
                       
                     
                     ) 
                   
                 
               
             
             
               
                 ( 
                 11 
                 ) 
               
             
           
         
       
     
   
   In the case that the porosity of the rock is already known, either of the equations (9) or (10) (defining dq 1 (t 1 ) or dq 2 (t 2 )) can be used to derive the permeability κ. 
   In the event that the porosity is not known then by considering the flow from a third segment, assumed to have similar reservoir characteristics as segments  1  and  2 , drilled at a third bottomhole pressure, three equations with three unknowns are obtained. 
   The solution method for three simultaneous equations can take many forms. For example, if we drill the third segment, segment  3 , over a timescale Δt n , where Δt n ≠Δt, then 
   
     
       
         
           
             
               
                 
                   
                     dq 
                     3 
                   
                   ⁡ 
                   
                     ( 
                     
                       t 
                       3 
                     
                     ) 
                   
                 
                 = 
                 
                   
                     ( 
                     
                       
                         P 
                         f 
                       
                       - 
                       
                         
                           P 
                           BHP 
                         
                         ⁡ 
                         
                           ( 
                           
                             
                               t 
                               3 
                             
                             , 
                             
                               z 
                               3 
                             
                           
                           ) 
                         
                       
                     
                     ) 
                   
                   ⁢ 
                   
                     
                       4 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       π 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       κ 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         r 
                         w 
                       
                     
                     μ 
                   
                   ⁢ 
                   
                     
                       Δ 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         z 
                         3 
                       
                     
                     
                       
                         log 
                         ⁡ 
                         
                           ( 
                           
                             
                               4 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               κ 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               Δ 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               
                                 t 
                                 n 
                               
                             
                             
                               ϕ 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               
                                 c 
                                 t 
                               
                               ⁢ 
                               μ 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               
                                 r 
                                 w 
                                 2 
                               
                             
                           
                           ) 
                         
                       
                       - 
                       γ 
                     
                   
                 
               
             
             
               
                 ( 
                 12 
                 ) 
               
             
           
         
       
     
   
   By using the equation describing the flow from segment  1  to write 
                   log   ⁡     (   ϕ   )       =       log   ⁡     (       4   ⁢           ⁢   κ   ⁢           ⁢   Δ   ⁢           ⁢   t       μ   ⁢           ⁢     c   t     ⁢     r   w   2         )       -   γ   -       (       P   f     -       P   BHP     ⁡     (       t   1     ,     z   1       )         )     ⁢       4   ⁢           ⁢   π   ⁢           ⁢   κ   ⁢           ⁢     r   w             dq   1     ⁡     (     t   1     )       ⁢   μ       ⁢   Δ   ⁢           ⁢     z   1                 (   13   )               
and substituting this into the equation (12) to give
 
   
     
       
         
           
             
               
                 
                   
                     dq 
                     3 
                   
                   ⁡ 
                   
                     ( 
                     
                       t 
                       3 
                     
                     ) 
                   
                 
                 = 
                 
                   
                     ( 
                     
                       
                         P 
                         f 
                       
                       - 
                       
                         
                           P 
                           BHP 
                         
                         ⁡ 
                         
                           ( 
                           
                             
                               t 
                               3 
                             
                             , 
                             
                               z 
                               3 
                             
                           
                           ) 
                         
                       
                     
                     ) 
                   
                   ⁢ 
                   
                     
                       4 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       π 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       κ 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         r 
                         w 
                       
                     
                     μ 
                   
                   ⁢ 
                   
                     
                       Δ 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         z 
                         3 
                       
                     
                     
                       
                         log 
                         ⁡ 
                         
                           ( 
                           
                             
                               Δ 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               
                                 t 
                                 n 
                               
                             
                             
                               Δ 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               t 
                             
                           
                           ) 
                         
                       
                       + 
                       
                         
                           ( 
                           
                             
                               P 
                               f 
                             
                             - 
                             
                               
                                 P 
                                 BHP 
                               
                               ⁡ 
                               
                                 ( 
                                 
                                   
                                     t 
                                     1 
                                   
                                   , 
                                   
                                     z 
                                     1 
                                   
                                 
                                 ) 
                               
                             
                           
                           ) 
                         
                         ⁢ 
                         
                           
                             4 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             π 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             κ 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             
                               r 
                               w 
                             
                           
                           
                             
                               
                                 dq 
                                 1 
                               
                               ⁡ 
                               
                                 ( 
                                 
                                   t 
                                   1 
                                 
                                 ) 
                               
                             
                             ⁢ 
                             μ 
                           
                         
                       
                     
                   
                   ⁢ 
                   Δ 
                   ⁢ 
                   
                       
                   
                   ⁢ 
                   
                     z 
                     1 
                   
                 
               
             
             
               
                 ( 
                 14 
                 ) 
               
             
           
         
       
     
   
   This is re-arranged to give 
                       κ   =       ⁢       μ     4   ⁢           ⁢   π   ⁢           ⁢     r   w         ⁢       log   ⁡     (       Δ   ⁢           ⁢   t       Δ   ⁢           ⁢     t   n         )       [         (       P   f     -       P   BHP     ⁡     (       t   1     ,     z   1       )         )     ⁢       Δ   ⁢           ⁢     t   1           dq   1     ⁡     (     t   1     )           -                           ⁢       (       P   f     -       P   BHP     ⁡     (       t   3     ,     z   3       )         )     ⁢       Δ   ⁢           ⁢     z   3           dq   3     ⁡     (     t   3     )           ]       -   1                   (   15   )               
an expression for the permeability of segments  1 ,  2  and  3 .
 
   Since permeability and pore pressure are now defined, equation (13) gives 
                 ϕ   =     exp   ⁡     [       log   ⁡     (       4   ⁢           ⁢   κ   ⁢           ⁢   Δ   ⁢           ⁢   t       μ   ⁢           ⁢     c   t     ⁢     r   w   2         )       -   γ   -       (       P   f     -       P   BHP     ⁡     (     t   1     )         )     ⁢           ⁢       4   ⁢           ⁢   π   ⁢           ⁢   κ   ⁢           ⁢     r   w             dq   1     ⁡     (     t   1     )       ⁢           ⁢   μ       ⁢   Δ   ⁢           ⁢     z   1         ]               (   16   )               
the porosity appropriate to segments  1 ,  2  and  3 .
 
   This process may be completed for a series of three segments, drilled at three different bottomhole pressures, throughout the entire drilling operation, to yield pore pressure, permeability, and porosity profile of the near wellbore region, with no prior information regarding these characteristics required. 
   Note that where no bottomhole flow rate is possible, the formation pressure becomes 
                   P   f     ≈       (         dq     N   -   1       ⁢         P   BHP     ⁡     (       t   N     ,     z   N       )       ⁡     [     1   +       Δ   ⁢           ⁢     z   N         Δ   ⁢           ⁢     z     N   -   1             ]         -       dq   N     ⁢       P   BHP     ⁡     (       t     N   -   1       ,     z     N   -   1         )           )       (         dq     N   -   1       ⁡     [     1   +       Δ   ⁢           ⁢     z   n         Δ   ⁢           ⁢     z     N   -   1             ]       -     dq   N       )               (   17   )               
where the Q&#39;s are the total volumetric output from the entire reservoir at bottomhole (as measured from surface and suitably corrected for bottomhole conditions).
 
   Again, permeability is calculated easily where the porosity is known from either segment. 
   Alternatively, a similar procedure to that outlined above may be used to determine both permeability and porosity, as well as pore pressure from surface measurements. 
   If the accuracy of the flow meter requires a target change in flow rate from the two individual zones compared to the total volumetric flow rate, then for detectability we have 
                     (       dq   N     -     dq     N   -   1         )         q   s     ⁡     (     t   N     )         =   T           (   18   )               
and writing P BHP (t N )=r N P BHP (t N-1 ), then we find
 
   
     
       
         
           
             
               
                 
                   r 
                   N 
                 
                 = 
                 
                   1 
                   - 
                   
                     
                       
                         
                           
                             q 
                             s 
                           
                           ⁡ 
                           
                             ( 
                             
                               t 
                               N 
                             
                             ) 
                           
                         
                         ⁢ 
                         
                           ( 
                           
                             1 
                             - 
                             T 
                           
                           ) 
                         
                       
                       - 
                       
                         
                           q 
                           s 
                         
                         ⁡ 
                         
                           ( 
                           
                             t 
                             
                               N 
                               - 
                               1 
                             
                           
                           ) 
                         
                       
                       - 
                       
                         
                           4 
                           ⁢ 
                           
                               
                           
                           ⁢ 
                           π 
                           ⁢ 
                           
                               
                           
                           ⁢ 
                           
                             r 
                             w 
                           
                           ⁢ 
                           
                             κ 
                             
                               N 
                               - 
                               1 
                             
                           
                           ⁢ 
                           
                             
                               P 
                               f 
                             
                             ⁡ 
                             
                               ( 
                               
                                 t 
                                 
                                   N 
                                   - 
                                   1 
                                 
                               
                               ) 
                             
                           
                           ⁢ 
                           Δ 
                           ⁢ 
                           
                               
                           
                           ⁢ 
                           
                             z 
                             
                               N 
                               - 
                               1 
                             
                           
                         
                         
                           μ 
                           ⁡ 
                           
                             ( 
                             
                               
                                 log 
                                 ⁡ 
                                 
                                   [ 
                                   
                                     Γ 
                                     
                                       N 
                                       - 
                                       1 
                                     
                                   
                                   ] 
                                 
                               
                               - 
                               γ 
                             
                             ) 
                           
                         
                       
                     
                     
                       
                         
                           P 
                           BHP 
                         
                         ⁡ 
                         
                           ( 
                           
                             t 
                             
                               N 
                               - 
                               1 
                             
                           
                           ) 
                         
                       
                       ⁢ 
                       
                         
                           ∑ 
                           
                             k 
                             = 
                             1 
                           
                           
                             
                               N 
                               ⁡ 
                               
                                 ( 
                                 
                                   t 
                                   N 
                                 
                                 ) 
                               
                             
                             - 
                             1 
                           
                         
                         ⁢ 
                         
                           
                             4 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             π 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             
                               r 
                               w 
                             
                             ⁢ 
                             
                               κ 
                               k 
                             
                             ⁢ 
                             Δ 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             
                               z 
                               k 
                             
                           
                           
                             μ 
                             ⁡ 
                             
                               ( 
                               
                                 
                                   log 
                                   ⁡ 
                                   
                                     [ 
                                     
                                       Γ 
                                       k 
                                     
                                     ] 
                                   
                                 
                                 - 
                                 γ 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                 
               
             
             
               
                 ( 
                 19 
                 ) 
               
             
           
         
       
     
   
   Before using the method in accordance with this invention, the required BHP variation r N , may be needed to obtain the target variation in flow rate, T. This can be obtained by using estimates of the formation pressure and permeability. Thereafter, derived values obtained using real data from the well can be used to update the values of r N  and T. 
   For a spatial resolution R required in the pore pressure profile, a timescale Δt is associated with the annular BHP variations of 
   
     
       
         
           
             
               
                 
                   Δ 
                   ⁢ 
                   
                       
                   
                   ⁢ 
                   
                     t 
                     
                       N 
                       - 
                       1 
                     
                   
                 
                 = 
                 
                   
                     Δ 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     
                       t 
                       N 
                     
                   
                   = 
                   
                     R 
                     
                       
                         U 
                         ⁡ 
                         
                           ( 
                           
                             t 
                             N 
                           
                           ) 
                         
                       
                       + 
                       
                         U 
                         ⁡ 
                         
                           ( 
                           
                             t 
                             
                               N 
                               - 
                               1 
                             
                           
                           ) 
                         
                       
                     
                   
                 
               
             
             
               
                 ( 
                 20 
                 ) 
               
             
           
         
       
     
   
   The method disclosed here is a means of deriving the reservoir pore pressure profile, real time, whilst drilling underbalanced. In this methodology, variations in bottomhole pressure during underbalanced drilling operations, and the subsequent variations in produced flow rates at surface, are interpreted in a manner which allows the local pore pressure to be obtained to high spatial resolution. 
     FIG. 3  shows a series of graphs illustrating simulated data and the pore pressure and permeability which can be derived from such data using an algorithm according to the method of the present invention. The same series of graphs can be achieved for real data, using bottomhole pressure over time, measured depth, surface and standpipe pressures and surface flow meter of gas into and out of the wellbore. 
     FIG. 3(   a ) shows the produced volumetric flux from the reservoir at bottomhole as a function of the distance drilled.  FIG. 3(   c ) shows fluctuations in bottomhole annular pressure as a function of the drilling depth. Using this data, and the expressions derived herein,  FIG. 3(   b ) shows the pore (or formation) pressure in MPa derived using the changes in pressure as a function of distance, with  FIG. 3(   d ) illustrating the permeability profile derived again as a function of distance.  FIG. 3(   e ) shows time as a function of measured depth. 
     FIG. 4  illustrates what is achieved when the same simulated data relating to borehole pressure and production flow is analysed by setting the pore or formation pressure to 10 MPa using a prior art algorithm which derives permeability using an estimated formation pressure.  FIG. 4(   a ) shows produced volumetric flux as a function of distance,  FIG. 4(   b ) shows the fluctuations in bottomhole annular pressure,  FIG. 4(   c ) shows the permeability profile derived using the estimated constant pore pressure of 10 MPa, and  FIG. 4(   d ) shows time as a function of measured depth. The prior art algorithm derives an incorrect permeability profile as shown in  FIG. 4(   c ) even in the case of a fairly homogenous but non-constant formation pressure profile as shown in  FIG. 3(   b ). 
   The measured variations in BHP shown in the example of  FIG. 3  are such that detectable variations in gas flow at surface may be derived over periods of one, to several hours. Accuracy of production rates is facilitated in these cases by adopting a steady injection rate. 
   The present invention can be used by measuring unintentionally caused pressure variations such as from uncontrolled variations in the mud pumping speed or variability of influx from the reservoir. Thus, unintentional variations in pressure can be used, so long as the rate change is sufficient for the resolution required given the particular drilling situation (for example, the rate of penetration, flow measurement accuracy). 
   According to another embodiment of the present invention, the pressure variations can be intentionally induced. According to a preferred embodiment the composition of the drilling fluid can be changed during drilling. This can be accomplished for example by changing the ratio of gas to liquid in the drilling fluid. Pressure variations can also be induced by changing the pumping rates of the drilling fluid, or actuating a moveable constriction in the system either downhole or on the surface. According to another preferred embodiment, the annular pressure of the drilling fluid at the surface can be altered using a choke unit. The variations can also be induced using a specially shaped section of pipe or nozzle that causes a resonance in the fluid pressure. 
     FIG. 5  shows a system for calculating properties relating to a subterranean formation, according to a preferred embodiment of the invention. Although derrick  44  is shown placed on a land surface  42 , the invention is also applicable to offshore and transition zone drilling operations. Borehole  46 , shown in dashed lines, is being formed in the subterranean formation  40  using bit  54  and drill string  58 . The lower portion of drill string  58  comprises a bottom hole assembly (“BHA”)  56 . The BHA  56  in turn, comprises a number of devices, including annular pressure sensor  60 , downhole flowmeter  70  and telemetry subassembly  64 . 
   At the surface  42 , are located the circulating system, not shown, for circulating the drilling fluid (which includes the mud pumps), rotating system, not shown, to rotate the drill string and drill bit, and a hoisting system, not shown, for suspending the drill string with the proper force. 
   According to the invention, data from the pressure sensor  60  and flow meter  70  are transmitted to the telemetry subassembly  64  via a cable, not shown. Telemetry subassembly  64  then converts the data from electrical form to some other form of signals, such as mud pulses. However, Telemetry subassembly  64  could use other types of telemetry such as torsional waves, in drill string  58 , or an electrical connection via a cable. The telemetry signals from subassembly  64  are received by a receiver, not shown, located in surface equipment  66 . The receiver converts the telemetry signals back into electronic form (if necessary) and then transmits the data to a logging unit  68  for recording and further processing. Logging unit  68  comprises a computer/data processor, data storage, display and control logic. 
   Also preferably provided in surface equipment  66  is are surface fluid pressure sensors that (1) measure the pressure of the fluid coming out of the annulus (i.e. the annular region between drillstring  58  and the borehole wall of borehole  46 , and (2) measure the standpipe pressure (i.e. the pressure of the fluid inside the drillstring  58 ). Surface equipment  66  also preferably comprises flow sensors that measure the flow rates of both injection and outflow. According to another embodiment of the invention, a choke unit is provided in surface equipment  66  for altering the pressure of the fluid. 
   According to an alternative embodiment, a coiled tubing drilling arrangement is used instead of derrick  66 , and drillstring  58 . In this case the data from flow meter  70  and pressure sensor  60  is transmitted via a wireline connection to the surface. 
   In operation, the computer located in logging unit  68  is used to calculate the properties such as pore pressure, porosity, and permeability using the data from the various sensors, according to the invention as herein described. 
     FIG. 6  shows steps involved in calculating properties relating to a subterranean formation, according to a preferred embodiment of the invention. Step  100  is the drilling process in which the borehole is formed in the subterranean formation. Although step  100  is shown as the first step in  FIG. 6 , in practice the other steps of the invention (e.g. step  102  to  108  in  FIG. 6 ) are carried out during the drilling step  100 . In step  102  the pressure and fluid flow rates are measured when the drilling has progressed to a certain point, or depth. According to preferred embodiment described above, the pressure in the borehole is measured using an annular pressure sensor located in the bottom hole assembly, and the fluid flow rate is either measured at the surface, or using a downhole flow meter. Additionally, although the drilling process  100  can be stopped during the measurement, according to a preferred embodiment, the measurements are taken as the drilling proceeds. In steps  104  and  106  the same or similar measurements are taken when the drilling has progressed to two other points. Finally, in step  108  the properties of the formation are calculated using the measurements. As has been described above, if only one or two properties are being calculated, then measurement from only two of the three locations are preferably used in the calculation step. 
   The above-described embodiments are illustrative of the invention only and are not intended to limit the scope of the present invention.