GB2274667A

GB2274667A - Methods for analysis of drillstring vibration using torsionally induced frequency modulation

Info

Publication number: GB2274667A
Application number: GB9401264A
Authority: GB
Inventors: John Duncan Macpherson
Original assignee: Baker Hughes Inc
Current assignee: Baker Hughes Holdings LLC
Priority date: 1993-02-01
Filing date: 1994-01-24
Publication date: 1994-08-03
Anticipated expiration: 2014-01-24
Also published as: NO940208L; NO307845B1; GB9401264D0; NO940208D0; US5321981A; GB2274667B

Description

2274667 METHODS FOR ANALYSIS OF DRILLSTRING VIBRATION USING TORSIONALLY

INDUCED FREQUENCY MODULATION This invention relates generally to methods and techniques for the analysis of vibration data in oil and gas well drilling. More particularly, this invention relates to the use of frequency modulation sidebands observed from detected frequency domain vibration data for (1) discriminating between surface and downhole vibratory sources and (2) the determination of the rotary speed of bottom-hole-assembly (BHA) components.

During drilling, various sources excite the drillstring. The amplitude of the resultant drillstring vibrations will depend on the level (severity) of the excitation, the system damping and the proximity of the excitation frequency to a natural frequency of the drillstring. When the frequency of any of the excitation sources is a natural frequency of the drillstring (axial, torsional or lateral) then the string resonates. Vibration levels are generally highest at resonance, but high level vibrations may exist in the drillstring, independent of drillstring resonance, whenever a high level of excitation is present.

Torsional vibration often manifests itself as a stick/slip action of the bottomhole assembly (BHA). Indicators of these downhole vibrations can be monitored at or near the surface. This stick/slip phenomenon is described in a paper entitled "A Study of Slip-Stick Motion of the BiC by A. Kyllingstad and G.W. Halsey, Society of Petroleum Engineers (SPE) Paper 16659, Sept. 1987. As discussed in that paper, torsional oscillations are caused by alternating slipping and sticking of the bottom-hole-assembly WHA) as it rotates in the borehole. This phenomenon is associated with a large amplitude, sinusoidal and often saw-tooth like variation in the applied torque. The term slip-stick motion refers to the belief that the amplitude of the torsional oscillations becomes so large that the 2 drillcollar section periodically comes to a complete stop and does not come free until enough torque is built up in the drillstring to overcome the static friction.

Drilling with large amplitude vibrations will result in accelerated drillstring fatigue. A recent study of drillstring failures (see Hill, T. H.; Seshadri, P.V.; Durham, K.S., "A Unified Approach to Drillstring Failure Prevention, " SPElIA0C 22002, March 1991, Amsterdam). Most current efforts aimed at understanding and controlling drillstring vibrations focus on the failure of drillstring components. However, during drilling, the drillstring transfers power from the surface to the bit and high amplitude drillstring vibrations may represent a loss, or waste, of drilling energy. Therefore, high levels of vibration not only result in drillstring component failures but can also result in suboptimum drill rates.

The avoidance of high vibration levels can be attempted in two ways: (a) the BHA can he modelled and a harmonic analysis performed to predict the operating conditions, weight-on-bit (WOB) and rotary speed (RPM), which avoid resonant conditions or (b) the vibrations can be directly monitored while drilling to determine the optimum operating conditions (WOB, RPM and pump rate).

The use of modelling analysis is limited by the number of unknowns that occur in real-time drilling operation. In addition, most models focus on drillstring resonance frequencies, and estimate bit rotary speeds which will avoid exciting these frequencies ("critical spee& analysis).

There is a continuing need for improved methods of using drillstring vibration data to enhance or optimize drilling performance andlor drillstring performance. For example, there is a need for a method to locate the source (e.g. surface or downhole source) of drillstring vibration so that corrective action may be taken 3 to prevent drillstring fatigue and less than optimum drilling rates. In addition, there is a continuing need for methods of using drillstring vibration data to determine the rotary speed of BHA components to enhance or optimize the life of the BHA components and the drilling performance.

The above-discussed and other problems and deficiencies of the prior art are overcome or alleviated by the methods for using vibration data for the optimization or enhancement of drilling andlor drillstring performance of the present invention as set out in any of the independent claims. Torsional oscillations of the drillstring will lead to frequency modulation WM) of the signal from a vibratory source (e.g., the bit). This results, in the frequency domain, in sidebands being present around a detected excitation frequency. These sidebands may be used in advantageous methods for optimizing drillstring and drilling performance.

In a first embodiment of this invention, these sidebands are used to discriminate between downhole and surface vibrational sources caused by torsionally induced frequency modulation. The number of sidebands depends on the ratio of the max angular frequency due to the oscillatory motion and the frequency of these oscillations (termed the modulation index). Since the oscillation frequency is a constant for a given drillstring length, drillstring configuration and wellbore, the number of sidebands is directly proportional to the maximum angular frequency. there is a maximum of these sidebands at the bit (the 'working end') of the pendulum, and a minimum (or zero) of these sidebands at the surface. Therefore, surface vibratory sources will be distinguishable by the absence of sidebands (zero modulation index). Downhole vibratory sources will have sidebands, the number of sidebands (and the modulation index) increasing from the surface to the distal end of the pendulum.

4 Once the location of the drillstring vibration is determined, appropriate action such as changing drilling parameters, for example weight-on-bit (WOB) or mud properties may be taken to optimize drilling and drillstring performance.

In a second embodiment of this invention, the sidebands are used to determine the rotary speed of BHA components. Using the method of this second embodiment, minimum and maximum rotary speeds of a given BHA component is determined as a function of the excitation of the frequency, the frequency of torsional oscillation and the modulation index. Once the minimum and maximum rotary speeds of the BHA components are determined, adjustments can be made to alter the rotary speed of the drillstring and thereby enhance or optimize drilling and drillstring performance. This method is particularly well suited for use in those applications where torsional oscillations are not recognizable in the time domain, but are better recognized in the frequency domain.

The above-described and other features and advantages of the present invention will be appreciated and understood by those skilled in the art from the following detailed description and drawings, in which:

Fig. 1 is a graphical depiction for a drillstring, showing frequency modulation of the signal from a vibration source; Fig. 2 is a graphical depiction for a drillstring showing, in the frequency domain, sidebands being present around the detected excitation frequency; Fig. 3 is a plot of amplitude vs frequency for a drillstring for discriminating between downhole and surface vibrations; Fig. 4 is a cross-sectional schematic side elevation view of a drillstring vibration measurement sub:

Fig. 5 is a graph depicting frequency vs amplitude for data set at 25 Hz; and Fig. 6 is a graph depicting amplitude vs frequency for the time domain data set of Fig. 5.

While drilling wells it is possible for the drilistring to act as a torsional pendulum (also referred to as 'stick-slip' motion). This will result in frequency modulation WM) of any excitation source (e.g., the bit), the degree of modulation being dependent upon the location of the source in the string. Fig. 1 is a graphical depiction showing frequency modulation of the torsional signal from a vibratory source in a drillstring. Sources at the bit will experience a high degree of modulation; sources at the surface should not experience modulation.

Frequency modulation is apparent in the time domain as a series of periodic 'beats'. However, beats can also be generated by closely spaced excitation sources or by amplitude modulation. Thus, under certain circumstances, such periodic 'beats' in the time domain caused by torsional oscillations (e.g., stickslip phenomenon) do not provide a recognizable signature in the time domain. In contrast, in the frequency domain, FM is readily distinguishable since it will generate several sidebands around the modulated frequency. Fig. 2 is a graphical depiction showing, in the frequency domain, sidebands being present around the detected excitation frequency. The spacing between the sidebands is governed by the modulating frequency (the frequency of oscillation of the torsional pendulum - the periodicity of the 'stick-slip' motion). The number of sidebands depends on the ratio of the max angular frequency due to the oscillatory motion and the frequency of these oscillations (termed the modulation index). Since the oscillation frequency is a constant for a given drillstring length, drillstring 6 configuration and wellbore, the number of sidebands is directly proportional to the maximum angular frequency. This is a maximum at the bit (the 'working end') of the pendulum, and a minimum (or zero) at the surface.

Therefore, in accordance with a first embodiment of this invention, surface vibratory sources are distinguishable by the absence of sidebands (zero modulation index). In contrast, downhole vibratory sources will have sidebands, the number of sidebands (and the modulation index) increasing from the surface to the distal end of the pendulum. Thus, referring to Fig. 3, a plot of amplitude vs frequency is depicted for a drillstring. The data is from a vertical onshore well at a depth of approximately 10,900 ft (3370 M), drilling with a rollercone bit and mud motor surface rotary speed of 68, WOB of 26,000 Ibs (116 kN), pump rate of 105 with triplex pumps. As indicated on Fig. 3, downhole vibration sources show sidebands downhole with such sidebands decreasing along the drillstring towards the surface where there are no apparent sidebands.

Once the source of the vibration is determined (e.g., discriminated as originating from downhole or at the surface) using the method of the first embodiment of this invention, corrective action may be taken to optimize drilling performance and preclude drillstring fatigue. Such corrective action may include changing drilling parameters, for example, weight-on-bit (WOB) or mud properties.

Measurement of drillstring vibrations in order to obtain the data for the graphs of Fig. 3 may be obtained from any number of known vibration measurement systems (located downhole or at the surface). Preferably, the drillstring vibration data is collected from a surface measurement system of the type described in Besaisow, A.A., Jan, Y.M., Schuh, F.J., "Detection of Various Drilling Phenomena Utilizing High Frequency Surface Measurements% 1985, SPE 14327 Las Vegas and U.S. Patent No. 4,715,451, all of the contents of which are 7 incorporated herein by reference. A commercial version of this vibration measurement system is known as the ADAMS Mark 3 sub which is used in vibration monitoring services offered by EXLOG, INC. under the servicemark DynaByte. Referring to Fig. 4, this sub is shown at 10 and contains a suite of sensors, consisting of strain gauges 12 and accelerometers 14, mounted on a 4145H modified steel sub 16. The sub 16 fits into the drillstring below the kelly swivel.

Full bridge semiconductor strain gauges 12 are used to measure dynamic axial force and dynamic torsional moment; foil type gauges are used to measure string weight and torque. Paired accelerometers 14 are used to measure axial and torsional acceleration. Several ancillary sensors are also used; these include a magnetometer for surface rotary speed, a pressure transducer for pump pressure and operational sensors such as battery voltage and temperature.

A specialized data acquisition system 18 contained within the sub 16 housing digitizes and encodes data from the sensors. Sample rates are programmable and typically 2083 samples per second per channel is used. Anti-aliasing filters are used prior to sampling to ensure a non-aliased measurement with a 500 Hz bandwidth. The measured data is transferred to an on-site unit, at a rate of 250 kbits/sec, utilizing a microwave transmitter and omnidirectional antenna mounted within the sub housing. In addition to efficiently transferring data from a rotating medium, microwave telemetry also alleviates the need to run data and power lines to the rig floor. Power to the system is provided by a removable and rechargeable battery pack.

In accordance with a second embodiment of the present invention, the sidebands of Fig. 2 are used to determine the rotary speeds of a B11A component. The spacing between sidebands depends upon the periodicity of the torsional 8 oscillations. In accordance with this second embodiment, it has been determined that the minimum and maximum rotary speeds of the BHA component can be found from the following:

RPM min 60 We - a Fm] 1 (a) ROM max 60 [Fe + a Fm] 1 (b) where, Fe = the excitation frequency Fm = the frequency of torsional oscillation a = the modulation index The excitation frequency Fe is the center peak in the 'bundle'.

The frequency of torsional oscillations Fm is given by the frequency spacing of the sidebands.

The modulation index a can be found as follows:

The amplitude of each peak in the 'bundle' is given by the following:

Main Peak (excitation): AO x JO (a) 1st sideband: AO x J1 (a) 2nd sideband: AOX J2 (a) Nth sideband: AO X JN (a) where, AO is the original (unmodulated) amplitude of the excitation; a is the modulation index; and 9 JO J,, are Bessel functions of the first kind evaluated at a.

It will be appreciated from the above that the amplitude ratios of adjacent peaks are given by the ratios of the Bessel functions evaluated at the modulation index. There is, therefore, a direct correspondence between the ratios of the Bessel functions and the peak amplitude ratios. For example, Table 1 shows the amplitude ratio of the first sideband to the center peak versus the modulation index (the table is only one-to-one for values of the modulation index less than approximately 2.4). If the amplitude ratio of the peaks is calculated, then the modulation index can be easily read from the table. Of course, it will be appreciated that the table method may be replaced by computerized software using a more enhanced scheme taking into account ratios of other sidebands.

The method of the second embodiment of this invention may be summarized as follows:

(1) Divide the amplitude of the first sideband by the amplitude of the center peak.

(2) Use the provided table to find the modulation index (or a software programme).

(3) Use Equation la and 1b to find the minimum rotary speed of the drillstring component.

TABLE 1

AMP MOD AMP MOD AMP MOD PATIO INDEX PATIO INDEX RATIO INDEX 0.005 0.010 0.269 0.520 0.598 1.030 0.010 0.020 0.275 0.530 0.606 1.040 0.015 0.030 0.280 0.540 0.614 1.050 0.020 0.040 0.286 0.550 0.622 1.060 0.025 0.050 0.292 0.560 0.630 1.070 0.030 0.060 0.297 0.570 0.638 1.080 0.035 0.070 0.303 0.580 0.646 1.090 0.040 0.080 0.309 0.590 0.654 1.100 0.045 0.090 0.314 0.600 0.633 1.110 0.050 0.100 0.320 0.610 0.671 1.120 0.055 0.110 0.326 0.620 0.680 1.130 0.060 0.120 0.332 0.630 0.688 1.140 0.065 0.130 0.338 0.640 0.697 1.150 0.070 0.140 0.343 0.650 0.706 1.160 0.075 0.150 0.349 0.660 0.715 1.170 0.080 0.160 0.355 0.670 0.724 1.180 0.085 0.170 0.361 0.680 0.733 1.190 0.090 0.180 0.367 0.690 0.742 1.200 0.095 0.190 0.373 0.700 0.752 1.210 0.101 0.200 0.379 0.710 0.761 1.220 0.106 0.210 0.386 0.720 0.771 1.230 0.111 0.220 0.392 0.730 0.781 1.240 0.116 0.230 0.398 0.740 0.791 1.250 0.121 0.240 0.404 0.750 0.801 1.260 0.126 0.250 0.410 0.760 0.811 1.270 0.131 0.260 0.417 0.770 0.821 1.280 0.136 0.270 0.423 0.780 0.831 1.290 0.141 0.280 0.429 0.790 0.842 1.300 0.147 0.290 0.436 0.800 0.853 1.310 0.152 0.300 0.442 0.810 0.863 1.320 0.157 0.310 0.449 0.820 0.874 1.330 0.162 0.320 0.455 0.830 0.886 1.340 0.167 0.330 0.462 0.840 0.897 1.350 0.173 0.340 0.469 0.850' 0.908 1.360 0.178 0.350 0.475 0.860 0.920 1.370 0.183 0.360 0.482 0.870 0.932 1.380 0.188 0.370 0.489 0.880 0.944 1.390 0.194 0.380 0.496 0.890 0.956 1.400 0.199 0.390 0.503 0.900 0.968 1.410 0.204 0.400 0.510 0.910 0.981 1.420 0.209 0.410 0.517 0.920 0.994 1.430 0.215 0.420 0.524 0.930 1.007 1.440 0.220 0.430 0.531 0.940 1.020 1.450 0.226 0.440 0.538 0.950 1.034 1.460 0.231 0.450 0.545 0.960 1.047 1.470 0.236 0.460 0.553 0.970 1.061 1.480 0.242 0.470 0.560 0.980 1.076 1.490 0.247 0.480 0.568 0.990 1.090 1.500 0.253 0.490 0.575 1.000 1.105 1.510 0.258 0.500 0.583 1.010 1.120 1.520 0.264 0.510 0.590 1.020 1.135 1.530 Table 1 (Continued) AMP MOD AMP MOD RATIO INDEX RATIO INDEX 1.151 1.540 2.936 2.050 1.167 1.550 3.021 2.060 1.183 1.560 3.110 2.070 1.200 1.570 3.204 2.080 1.216 1.580 3.304 2.090 1.234 1.590 3.411 2.100 1.251 1.600 3.525 2.110 1.269 1.610 3.646 2.120 1.288 1.620 3.777 2.130 1.307 1.630 3.916 2.140 1.326 1.640 4.067 2.150 1.346 1.650 4.230 2.160 1.366 1.660 4.406 2.170 1.387 1.670 4.598 2.180 1.408 1.680 4.808 2.190 1.429 1.690 5.038 2.200 1.452 1.700 5.291 2.210 1.475 1.710 5.571 2.220 1.498 1.720 5.883 2.230 1.522 1.730 6.233 2.240 1.547 1.740 6.627 2.250 1.572 1.750 7.075 2.260 1.598 1.760 7.590 2.270 1.625 1.770 8.187 2.280 1.653 1.780 8.887 2.290 1.681 1.790 9.720 2.300 1.710 1.800 10.729 2.310 1.741 1.810 11.975 2.320 1.772 1.820 13.553 2.330 1.804 1.830 15.617 2.340 1.837 1.840 18.433 2.350 1.872 1.850 22.505 2.360 1.907 1.860 28.913 2.370 1.944 1.870 40.483 2.380 1.982 1.880 67.655 2.390 2.021 1.890 207.437 2.400 2.062 1.900 2.105 1.910 2.149 1.920 2.195 1.930 2.242 1.940 2.292 1.950 2.344 1.960 2.398 1.970 2.455 1.980 2.514 1.990 2.576 2.000 2.641 2.010 2.709 2.020 2.781 2.030 2.857 2.040 12 As discussed, torsional oscillation produces a frequency modulated signal; this modulation ranges from zero (if the bit comes to a complete stop) the maximum bit speed as the bit "unwinds". Axial vibration levels will vary with the bit speed. Therefore, when the bit is rotating at lower speeds the axial vibration level is lower, and conversely when the bit rotates at higher speeds the vibration levels are higher. This produces visible amplitude modulation (or 'beats') of the axial signals detected at the surface; it may also produce frequency modulation of the monitored axial signals. In the example of Fig. 5, however, the characteristic 'beats' are not apparent. The data is from a vertical onshore well at a depth of approximately 10,900 ft (3370 m), drilling with a rollercone bit and mud motor surface rotary speed of 68, WOB of 26,000 Ibs (116 kN), pump rate of 105 with triplex pumps. Significantly, in the amplitude spectra of this data set, shown in Fig. 6, the frequency modulation of the bit RPM signal is easily identified by the presence of multiple sidebands. These sidebands have a frequency spacing equal to the modulating frequency of 0.18 Hz (period of 5.56 seconds). The surface rotary signal is also modulated; however, the degree of modulation is less than that imposed on the bit RPM signal.

Once frequency modulation is established, the minimum and maximum rotary speeds may be estimated using Equations 1 (a) and 1 (b) as described above. For the example shown, the RPM range experienced by the bit is 261 - 279 RPM. It will be appreciated such a small oscillation is difficult to detect on static torsional channels.

Bit Fe 4.50 Hz Fm 0.18 Hz Amplitude Ratio - 4.3 Modulation Index = 0.85 13 RPM min = 60 [4.5 - 0.8 (0.18fl = 261 RPM RPM max = 60 [4.5 + 0.8 (0.18fl = 279 RPM RPM ave = 270 RPM.

Once the rotary speed of the B11A components is determined using the method of the second embodiment of this invention, the drilling operation can be optimized or improved by, for example, altering the rotary speed of the drillstring.

Claims

14 CLAIMS

1. A method for analyzing drillstring vibration and optimizing drillstring or drilling performance comprising the steps of:

(1) detecting vibratory signals in a drillstring; (2) converting said detected vibratory signals to the frequency domain wherein sidebands are exhibited around at least one or more detected excitation frequencies due to frequency modulation; (3) using said sidebands to analyze drillstring vibration; and (4) optimizing drillstring or drilling performance based on the analysis of step (3).

2. The method of Claim 1 wherein step (3) further comprises the step of:

discriminating between downhole drillstring vibrational sources and surface drillstring vibrational sources using torsionally induced frequency modulation wherein the number of sidebands decrease from a downhole location towards the surface.

3. The method of Claim 2 wherein:

surface drillstring vibrations exhibit an absence of sidebands; and downhole drillstring vibrations exhibit sidebands which increase from the surface toward a distal end of the drillstring.

4. The method of Claim 1 wherein step (3) further comprises the step of:

determining the rotary speed of at least one bottom hole assembly (BHA) component as a function of the excitation frequency, the frequency of torsional oscillation and the modulation index.

5. The method of Claim 4 wherein step (3) further comprises determining is the minimum rotary speed of BHA component (RPM min) and the maximum rotary speed of a BHA component (RPM max) as follows:

RPM min 60 We - a Fm] RPM max 60 [Fe + a Fm] where, Fe the excitation frequency Fm the frequency of torsional oscillation a = the modulation index.

6. A method of discriminating between downhole drillstring vibrational sources and surface drillstring vibrational sources and optimizing drillstring or drilling performance comprising the steps of:

(1) detecting vibratory signals in a drillstring; (2) converting said detected vibratory signals to the frequency domain wherein sidebands are exhibited around at least one or more detected excitation frequencies due to frequency modulation; (3) using said sidebands to discriminate between downhole drillstring vibrational sources and surface drillstring vibrational sources using torsionally induced frequency modulation wherein the number of sidebands decrease from a downhole location towards the surface; and (4) optimizing drillstring or drilling performance based on the discrimination of step (3).

7. The method of claim 2 wherein:

16

8. A method for determining the rotary speed of at least one bottom hole assembly (BHA) component and optimizing drillstring or drilling performance comprising the steps of:

(1) detecting vibratory signals in a drillstring; (2) converting said vibratory signals to the frequency domain wherein sidebands are exhibited around at least one or more detected excitation frequencies due to frequency modulation; (3) using said sidebands to determine the rotary speed of at lest one bottom hole assembly (BHA) component as a function of the excitation frequency, the frequency of torsional oscillation and the modulation index; and (4) optimizing drillstring or drilling performance based on the determination of rotary speed of step (3).

9. The method of Claim 8 wherein step (3) further comprises determining the minimum rotary speed of BHA component (RPM min) and the maximum rotary speed of a BHA component (RPM max) as follows:- RPM min 60 [Fe - a Fml RPM max 60 [Fe + a Fml where, Fe = the excitation frequency Fm = the frequency of torsional oscillation a = the modulation index.

10. A method for analyzing drillstring vibration substantially as hereinbefore described, with reference to the drawings.

11. A method discriminating between downhole drillstring vibration sources and surface drillstring vibration sources substantially as hereinbefore described, 4 17 with reference to the drawings.

12. A method for determining the rotary speed of a component substantially as hereinbefore described, with reference to the drawings.