DETERMINATION OF DRILL BIT RATE OF PENETRATION FROM SURFACE MEASUREMENTS |
|||||||
申请号 | EP93904240.4 | 申请日 | 1993-02-22 | 公开(公告)号 | EP0626032B1 | 公开(公告)日 | 1995-12-13 |
申请人 | SERVICES PETROLIERS SCHLUMBERGER; Anadrill International SA; | 发明人 | BOOER, Anthony K., "Orchard House"; | ||||
摘要 | A method of determining the rate of penetration Δd of a drill bit at the end of a drill string while drilling a well, comprising: (a) measuring the vertical displacement S of the drill string at the surface, (b) determining a state space description of S comprising a state space measurement equation (I), and a state evolution equation (II), wherein S and Δd are as previously defined, Δ is the difference operator for time τ between adjacent samples k and k+1, Μ is the drill string compliance, ς is the noise term associated with the surface displacement measurement, r is the noise term associated with fluctuations in the state, and h is the hookload, and (c) applying a Kalman filter to said equations to obtain an estimate of the state parameters including Δd. | ||||||
权利要求 | |||||||
说明书全文 | The present invention relates to a method of determining the rate of penetration (ROP) of a drill bit from measurements made at the surface while drilling. In the rotary drilling of wells such as hydrocarbon wells, a drill bit is located at the end of a drill string formed from a number of hollow drill pipes attached end to end which is rotated so as to cause the bit to drill into the formation under the applied weight of the drill string. The drill string is suspended from a hook and as the bit penetrates the formation, the hook is lowered so as to allow the drill string to descend further into the well. The ROP has been found to be a useful parameter for measuring the drilling operation and provides information about the formation being drilled and the state of the bit being used. Traditionally, ROP has been measured by monitoring the rate at which the drill string is lowered into the well at the surface. However, as the drill string, which is formed of steel pipes, is relatively long the elasticity or compliance of the string can mean that the actual ROP of the bit is considerably different to the rate at which the string is lowered into the hole. The errors which can be caused by this effect become progressively larger as the well becomes deeper and the string longer, especially if the well is deviated when increased friction between the string and the borehole wall can be encountered. Certain techniques have been proposed to overcome these potential problems. In both US 2,688,871 and US 3,777,560 the drill string is considered as a spring, and, the elasticity of the string being calculated theoretically from the length of the drill string and the Young's modulus of the pipe used to form the string, this information is then used to calculate ROP from the load applied at the hook suspending the drill string and the rate at which the string is lowered into the well. These methods suffer from the problem that no account is taken of the friction encountered by the drill string as a result of contact with the wall of the well. FR 2038700 proposes a method to overcome this problem in which the modulus of elasticity is measured in situ. This is achieved by determining the variations in tension to which the drill string is subjected as the bit goes down the well until it touches the bottom. Since it is difficult to determine exactly when the bit touches the bottom from surface measurements, strain gauges are provided near the bit and a telemetry system is required to relay the information to the surface. This method still does not provide measurements when drilling is taking place and so is inaccurate as well as difficult to implement. By contrast, in FR 2,165,851 (= AU 44,424/72) there is employed a mathematical model describing the drill bit cutting rate - the model necessitates a knowledge of the drill depth, the drill rotational speed, and the weight on bit, and its use involves the application of a Kalman-Bucy filter - to derive an ROP value. This method suffers from the obvious problems of having to know what is really going on at the bit, and the model utilised applies only to roller cone bits. The later GB 2,129,141 A tries to deal with the problem in a related way, applying Kalman filtering to a model that treats the drillstring as an elastic cable. and provides a downhole bit-acceleration measurement device (together with a "motionless tool" sensor necessary for correcting certain errors in depth measurement). Though quite useful, this method, like that of the aforementioned FR 2,165,811. suffers from its requirement for knowledge of downhole conditions. A simpler method is proposed in US 4,843,875 (incorporated herein by reference) in which ROP is measured from surface measurements while drilling is taking place. This method uses the following model: It is an object of the present invention to provide a method of determining ROP from surface measurements which can he used where the approach outlined above is undesirable or inappropriate. In accordance with one aspect of the present invention, there is provided a method of determining the rate of penetration Δd of a drill bit at the end of a drill string while drilling a well, comprising:
The present invention uses the same basic model as that of the aforementioned US 4,843,875 formulated in state space, and uses Kalman filtering as a continuously adaptive technique to solve the state parameters. As is explained in more detail hereinafter, the actual surface vertical displacement s when measured becomes the approximate value S because of various uncertainties and inaccuracies referred to generally as "noise" (the term "ρ"). The present invention will now be described, by way of example, with reference to the accompanying drawings in which:
Referring now to the drawings, the data shown in Figures 1 and 2 are obtained from experimental apparatus designed to provide realistic drilling data in controllable conditions. Such apparatus is described in US 4928521 which is incorporated herein by reference. The two examples from the experimental apparatus demonstrate the difficulties with ROP calculations. Figure 1(a) shows the raw depth measurement from an experiment in a drilling test machine with a PDC bit drilling marble. The derivative of this measurement, calculated by differencing adjacent points, is shown in 1(b). A "noise" level of about ±2 mm is apparent, and totally masks the underlying trend. Smoothing this derivative, as shown in Figure 1(c) (10 second averaging used) yields an indication of the ROP, but the estimate still has a high variance and the averaging has introduced a damped response to sharp changes in weight on bit. Further reduction of the variance by increasing the averaging time will result in a steady state estimate of ROP never being achieved for the finite duration drilling segments. Another example is shown in Figure 2, taken from a test in a different drilling machine. Figure 2(a) shows the depth measurement and 2(b) its derivative. Here again, the derivative calculation is very noisy, but the nature of the noise is different - it is not due to vibrations, but to quantisation (about 0.2 mm steps) in the original depth measurement. Figure 2(c) shows a 2 second average of the depth derivative. The underlying ROP trend is apparent but the variance due to measurement quantisation is still high. Increasing the averaging time would blur the boundaries between the different drilling segments. Both these examples demonstrate the problem with the direct calculation of ROP as a derivative of depth. Vibrations and measurement quantisation noise are also observed in field measurements. An alternative approach to ROP estimation is provided by the present invention by the use of a state-space approach. A state-space model comprises two equations: a measurement equation describing how observable measurements relate to the state vector, and a state evolution equation showing how the components of the state vector evolve in time. The state vector itself is a complete description of the system and contains parameters to be estimated. The state-space model applicable to the ROP problem has a state vector X with components: displacement s, surface ROP Δs. compliance Λ and downhole ROP Δd The observed parameter is displacement s, so the measurement equation (H = measurement matrix) is simply and the state evolves in this manner (Φ = state transition matrix) where τ is the time interval between adjacent sampling instants indicated by subscripts k and k+1. The depth measurement itself will contain noise and the above "model" chosen to represent the system will not be exactly true (e.g. there may be perturbing accelerations). The measurement and state evolution equations can be modified to include additive noise components (ρk, rk) which account for these discrepancies. In a general formulation for state space models, the matrices H and Φ may also be time-varying. Using conventional notation (y = observed output values) we have The measurement noise variance R can be estimated from the depth derivative. In the case of the data shown in Figure 1, the standard deviation of the noise is calculated to be ∼ 1 mm, so R = 1 x 10⁻⁶. For the Figure 2 data, the quantisation step size controls the variance, giving R = 4 x 10⁻⁸. An estimation of Q may be made by considering a perturbing acceleration ak. so the state covariance is where the variance of the acceleration, is chosen on the basis of knowledge of expected ROP variations. The ratio between R and σ Since the Kalman filter is a recursive estimator, initial conditions are required for X̂ and P. In the following examples, the initial conditions The above processing has been applied to the two drilling machine examples previously shown. Figure 3(a) shows the ROP estimate for the Figure 1 data and should be compared with Figure 1(c), the conventional ROP estimate, σ Figure 4 shows the processing applied using the data shown in Figure 2. Here the choice of σ In the following Appendices, Appendix A gives the Kalman filter equations, Appendix B gives a generalised code in Matlab to implement the Kalman filter, and Appesndix C gives an example of use of the code. |