Methodology for determining the parameters of drilling mode for directional straight sections of well using screw downhole motors Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

ê Vladimirs. Litvinenko, Mikhail V. Dvoynikov

Methodology for determining the parameters of drilling mode.

UDC 622.24

Methodology for determining the parameters of drilling mode for directional straight sections of well using screw downhole motors

Vladimir S. LITVINENKO, Mikhail V. DVOYNIKOV»

Saint Petersburg Mining University, Saint Petersburg, Russia

Article presents results of study on possibility of increasing the efficiency of drilling directional straight sections of wells using screw downhole motors (SDM) with a combined method of drilling with rotation of drilling string (DS). Goal is to ensure steady-state operation of SDM with simultaneous rotation of DS by reducing the amplitude of oscillations with adjusting the parameters of drilling mode on the basis of mathematical modeling for SDM - DS system.

Results of experimental study on determination of extrema distribution of lateral and axial oscillations of SDM frame depending on geometrical parameters of gerotor mechanism and modes ensuring stable operation are presented.

Approaches to development of a mathematical model and methodology are conceptually outlined that allow determining the range of self-oscillations for SDM - DS system and boundaries of rotational and translational wave perturbations for a heterogeneous rod with an installed SDM at drilling directional straight sections of well. This mathematical model of SDM - DS system's dynamics makes it possible to predict optimal parameters of directional drilling mode that ensure stable operation of borehole assembly.

Key words: well drilling; screw downhole motor; vibration; drilling string

How to cite this article: Litvinenko V.S., Dvoynikov M.V. Methodology for determining the parameters of drilling mode for directional straight sections of well using screw downhole motors. Journal of Mining Institute. 2020. Vol. 241, p. 105-112. DOI: 10.31897/PML2020.1.105

Introduction. During drilling of extended directional and horizontal sections of wells using volumetric principle engines, part of axial load on the bit is not transmitted due to frictional force arising between walls of well and drilling tool [2].

To ensure required load on the bit a combined drilling method is used in production process. Distinctive feature of this method is in the joint operation of drilling string (DS) and screw down-hole motor (SDM) [9]. In the process of their joint work, torsional, lateral and axial oscillations can occur depending on type of SDM, its energy characteristics and DS, which acts as an elastic unbalanced rod [4, 6, 10].

It should be noted that SDM that is located in lower part of DS has its own beating of the frame, nature of occurrence of which is associated with work of its power section, represented by a planetary reductor. Moreover, frequency, amplitude and direction of the frame beats depend on design of gerotor mechanism, hydraulic component of drilling mud flow, as well as load on the bit [3].

To determine parameters of well drilling mode by a combined method, it is necessary to develop a technique that allows providing forecast and control of stable operation of borehole assembly (BHA), based on mathematical modeling of elastic properties of DS stress-strain state and characteristics of SDM [11, 12].

Methodology and results of the research. Stability of SDM operation is characterized by working mode of power section, in which there is no intensive decrease in rotor rotation frequency with increasing torque on motor shaft [1].

It is known that axis of rotor rotates around its own axis, and also makes a transferring movement around axis of stator, directed counterclockwise. Moreover, frequency of transferring (planetary) rotation of rotor's axis relative to stator's axis is higher than rotor rotation frequency around its own axis.

ê Vladimirs. Litvinenko, Mikhail V. Dvoynikov

Methodology for determining the parameters of drilling mode.

Angular rotation velocity of rotor's axis relative to stator's axis, which determines beat frequency of the frame,

ro = -zz ror, (1)

where zz - number of rotor teeth; ror - angular rotation velocity of rotor around its own axis. Motor's frame beats depend on inertial Fin and hydraulic Fh forces acting on rotor,

2

Fin = mzn ro e,

(2)

Fh = Mind/ezi,

(3)

where Mind - indicator moment; e - eccentricity; zi and z2 - number of stator and rotor teeth; m -rotor mass; ro - angular velocity.

During engine start, a skew moment arises, causing instability of rotor rolling along stator teeth and leading to additional beating of SDM frame. Skew moment is

PpDt2

Ms =-

4n

(4)

where D - stator diameter at tooth cavities; Pp - pressure difference; t - rotor pitch.

Experimental study of motor's frame beats is performed at the test bench. Bench is equipped with an automatic control system that provides real-time output of SDM main energy characteristics to panel of a personal computer. To study beats of SDM, oscillation sensors are installed on frame.

Study results for vibration acceleration and amplitude of motor's frame oscillations under different operating conditions are shown in Fig. 1.

Based on experimental study, shaft rotation frequency is determined, which ensures minimal lateral oscillations and optimal axial beats of motor.

Modeling of tool operation is carried out on an advanced mathematical model of E.K.Yunin and V.K.Khegai [8].

At well drilling, it is required to determine the combination of load on the bit along depth P and rotor rotation frequency n0 so that drilling time t of specified interval is minimal under condition of optimal energy costs [7].

DS can be represented as a composite rod, including a section of length Li with outer and inner diameters DLi, dLi, interval of heavy drilling pipes (HDP) with length L2 and outer and inner diameters Dl , dLi and interval, represented by SDM frame and navigation system with length L3 and diameters DL , dL^. Current well depth H = L1 + L2 + L3 in process of drilling a certain interval increases due to deepening of bottomhole. At this, let us assume L2, L3 = const and due to increasing of L1 + AL, H also rises.

1 1 2 2 2 3 3 3 Sensor installation location

■ Longitudinal vibrations of SDM frame

>

400 300 200 100

1 1 1 2 2 2 3 3 3 Sensor installation location

Lateral vibrations of SDM frame

F ig.i Amplitude (a) and vibration acceleration (b) of frame depending on sensor installation location at the motor's frame 1 - upper SDM sub; 2 - middle of working bodies' active part of SDM; 3 - upper spindle sub

b

a

0

1

ê Vladimirs. Litvinenko, Mikhail V. Dvoynikov

Methodology for determining the parameters of drilling mode.

«i(si, t)

9l(Sl, t) U2(S2, t)

92(S2, t) U3(S3, t) 93(S3, t)

N

Si e [0,Li] d1 A-A Di

H

s2 e[0,L2] d2 B-B d2

S3 e[0,L3]

d3 C-C D3

P

Fig.2. Computational scheme for study of rotational and translational oscillations of SDM - DS system operation

Let us consider that sections are made of various materials. Therefore, first, second and third section corresponds to propagation velocity of rotational oscillations, respectively XL, XL , XL ,

propagation velocity of translational oscillations, respectively xL, XL, XL. Computational scheme

for analyzing DS behavior during rotational and translational motion is shown in Fig. 2.

Differential equation of rotational and translational motion of composite heterogeneous rod with boundary and initial conditions [5]:

a 2 9i , , ( a9i 1 ,2 a2 9i

at + ^{lT J = ÀLi as!2"

Si e [0, Li ];

a 2 9 2

at

a 2 9 3

at

+ A,

+ A

a9 2 1 ,2 a 2 9 2

at

Se

L2 2 > *2

aS

[0, L2 ];

ô9 31 ,2 a2 93

at

aS

2 , S3 e 0, L3

■ [0, L3 ].

a2ui ( aui 1 2 a2ui

at + f-i{nr h'-ist-s'ie[°, Li];

a2 u 2 ( du 21 2 a2 u 2

"öT + f X2 lit" l_x L2 "aST

a 2U3 + r ( aU3 1 _ „, 2 a 2U3

at u31 at J_xL as2'S3

e[0, L2 ]; (5)

s, e

[0, L3 ].

Boundary conditions for rotational and translational motion:

si = 0; 9 = n0t, M = GiJi

a9i as,

au

Si = 0; ui = h, N =

as.

Si = Li; S2 = 0; GiJi ^ GJ ^

as,

ôs,,

Si = Li; S2 = 0; EiFI—1 = E2F

aui = E F au2 ---E2F2——

as as0

si = Li; S2 = 0; 9 = 92,

si = Li; S2 = 0; ui = u2,

S

ê Vladimirs. Litvinenko, Mikhail V. Dvoynikov

Methodology for determining the parameters of drilling mode.

S2 = L2; S3 = 0; G2J2 = G3J3 , S2 = L2; S3 = 0; E2F2 = E3F3 ,

ds2 ÖS3 ds2 ÖS3

52 = L3; S3 = 0; 92 = 93, S2 = L3; S3 = 0; U2 = U3,

53 = L3; G3J3 ^ = -mh(p, «h). S3 = L3; E3F3 ^ = p(«h).

dS3 dS3

Initial conditions for rotational motion at t = 0:

9i(si, t = 0) =

_ ^ - / («") 2

2AT

2

S -

/ («0)L1

A2r

+e

Zt, K) L2

- + S

fT2(«0)L3 + MH(P, «H)

V L3

G3 J 3

92(S2, t = 0) = /1(Li) S22 -

2Ar

/T2 («0 )L2

- + s

/T3 («0 )L3 , Mh(P, «H)

v

V L3

+

GJ^,

93(S3, t = 0) = /1(Li) + /2(L2) + ^^ s2 -

2A2

2 J3

/T3 («0 )L3 , Mh(P, «H)

v

V L

- + -

G3 J 3

S,£ [0, Li], S"£ [0, L2], S3 £ [0, L3];

d^, _ d^2 Ö9

dt

= «,

dt

= «C

dt

S2 ;

2

2

S3;

= «0.

Initial conditions for translational motion at t = 0:

u,(s,, t = 0) =

_m_/Ti(«0) 2

"X

2 Ji

S -

/T1(«0) L

xl,

+e

/T2 («0 )L2

X L,

■ + S

^ /T 2 («0) L3 P(«h)A

+

EiFi

/

Uz(S2, t = 0) = /i(Li) + S22 -

2X

/T2(«0) L2

X L,

- + S

/T3(«0) L3 P(«h)

X L,

+

/T2(«0) 2

U,(S3, t = 0) = /i(Li) + /z(Lz) + S3 -

2X L

/T3(«0) L3 P(«h)

X 2,

33

Sie [0, Li], Sze [0, Lz], S3£ [0, L3];

du, _ du2 _ du 3 _

= XL, U1, = XL2 U2 , = X L3 U 3,

dt dt dt

>Si;

S2 ;

S3;

where 91(51, t), ^2(s2, t), 93(s3, t) - rotation angles of current cross-sections of string in corresponding sections; s1, s2, s3 - current position of cross-section; w1(s1, t), u2(s2, t), u3(s3, t) - translational movements of current cross-sections of string in corresponding sections; s1, s2, s3 - current position of cross-section; h - value of translational movement with transfer, lateral to truncation of round;

/J-^11 /J-^21, /x J-^3 J - dissipative members characterizing the resistance to rotation of

ê Vladimirs. Litvinenko, Mikhail V. Dvoynikov

Methodology for determining the parameters of drilling mode.

drilling string in corresponding sections; /.^-^j, /- dissipative members

characterizing the resistance of drilling string translational movement; n0 - rotation velocity for upper end of string; Gi, G2, G3 - shear modules for materials in corresponding sections; E1, E2, E3 -elastic moduli of materials in corresponding sections under tension or compression; Ji, J2, J3 - polar moments of inertia in cross-section of string in corresponding sections; Fi, F2, F3 - cross-sectional area of string in corresponding sections; MH(P, nH) - moment of resistance to rotation of composite rod lower section from the side of rock; P - axial load on end of composite rod lower section;

P(nH) - axial load on end of composite rod lower section; nH =

dt

- rotation frequency for

s, = L,

G J

end of composite rod lower section, 0 = 1 1 - coefficient of moment-force ratio of the first and

G2 J 2

second sections during rotation; s = G2J2 - coefficient of moment-force ratio of the second and

G3 J 3

E F

third sections during rotation; 0 = 1 1 - coefficient of moment-force ratio of the first and second

E2 F2

E F

sections during translational movement; s = 2 2 - coefficient of moment-force ratio of the second

E3 F3

and third sections during translational movement.

This problem is most clearly solved for case in which the values of dissipative members of system are equal to zero. Following formulas are used:

Ar

H =

A

LMn-

Al ^2

A

ln

A

AL, ln AL, AMH + G3 J3nH

Al

^3

ln

Ar

Al ^2

Ar

ln

Ar

Al

^3

ln

A L L3 AM h - G3 J3nH

A L L3 AM h + G3J 3nH

A L L3 AM H - G3 J3nH

A L L3 AM H + G3J 3nH

A L L3 AM H - G3 J3nH

A L L3 AM h + G3 J3nH

+ G2 J 2n0

G2 J 2n0

+ GiJino

H <

Al

ln

Ar

Al ^2

Ar

ln

Ar

A L

—— ln

^3 A L3 AM H " G3J 3nH

A L3 ln A L3 AM H + G3J 3nH

^3 AL3 AMh - G3 J3nH

+ G2J2n0

- G2 J2 n0

- G1 J1 n0

AL3 ln AL3 AMh + G3 J3nH

^3 A L3 AM H - G3J 3nH

+ G2J2

- G2 J2

+ G1 J1 n0

Ar

Al ^2

Ar

ln

Ar

AL3 ln AL3 AMh + G3 J3nH

^3 AL3 AMH - G3 J3nH

AL3 ln AL3 AMh + G3 J3nH

^3 AL3 amh - G3 J3nH

+ G2 J2

- G2 J2

- GJi

(6)

2

2

2

2

ê Vladimirs. Litvinenko, Mikhail V. Dvoynikov

Methodology for determining the parameters of drilling mode.

.3 L3

P = J („02 )1+ke 21L3

x,

.3 L3

1 - ke2X L

P =

P =

.3 L3 1 + ke X L3 P r ch V1+* L2 + I 2Xr 2Xr 2Xr V L1 l2 l3 y + kch .2L, 2Xr 2Xr 2Xr v L1 L 2 L 3 y y

.3 L3 1 - ke L3 n0 .3 L3 1 + ke X L3 P f sh v f sh v L2 + .3L3 ^ 2Xr 2Xr 2Xr v L1 L 2 L3 y L, ^ ^ 2Xr 2Xr 2Xr v L1 L 2 L3 y - ksh + ksh .2L2 .5l5 Y| 2Xr 2Xr 2Xr v L1 L 2 L 3 y y i.L .2L .A Y] 2Xr 2Xr 2Xr L1 L 3 L3 y

.3 L3 1 - ke L3 n0 f ch v L, + ^ 2Xr 2Xr 2Xr L1 L 2 L 3 1 - k2 - kch f.L .2l .3L Y| 2Xr 2Xr 2Xr v L1 L 2 L3 y y

sh

2Xr 2Xr 2À,r

L1 L2 l3 y

+ kch

.1L .2L 2 .3L3

L3 y

(7)

n =

L

2

where Li - DS length; L2 - HDP length; L3 - length of SDM frame and navigation system; k - coefficient of wave reflection, rotational oscillations at the interface of heterogeneous sections of a composite

k, — k- hr G J — hr , GnJn hr GnJn — hr GJ? .... „„

rod, k = , kj = L 1 1-v 2 2, k2 = L3 2 2-; m, ^3 - dissipation coeffi-

k1 + k2 h L2 G1J1 + h L G2 J2 h L3 G2 J2 + h L2 G3 J3

cients in corresponding sections, hl , hl , hl - propagation velocity of rotational oscillations in corresponding sections; PB and PH - axial loads on bottom end of SDM frame corresponding to upper and lower boundaries of self-oscillations; n* - rotor rotation frequency at PB = PH.

Task for case, in which the values of dissipative members of system are equal to zero, and propagation depth of translational oscillations of drilling tool, represented as a composite rod of three heterogeneous sections, is solved by system (6). At the same time G1, G2, G3 are replaced by E1, E2, E3 and J1, J2, J3 by F1, F2, F3, and also hL, hL , hL - propagation velocity of rotational oscillations by propagation velocity of translational oscillations by XL , XL, XL in corresponding sections. Obtained equations determine conditions for occurrence possibility of translational self-oscillations of DS, represented as a composite rod in process of translating to deepen bottomhole of well [13]. Axial loads on lower end of SDM frame, corresponding to upper and lower boundaries of self-oscillations during translational movement of PB and PH, are determined by the equation (7). At the same time propagation velocity of rotational oscillations hL , hL , hL is replaced by

Xl , XL, Xl , elastic modulus G1, G2, G3 and polar moments of inertia in cross-section J1, J2, J3 are

replaced by E1, E2, E3 and F1, F2, F3 respectively.

As a result of calculations based on developed mathematical model (5), range of self-oscillations' onset during rotation and translational movement of SDM - DS system was revealed. Input parameters for calculating rotational and translational movements:

L1 = 1800 m; L2 = 190 m; L3 = 10 m; J1 = 5.84M0-6 m4; J2 = 1.94M0-6 m4; J3 = 4.928-10-6 m4; k = 0.106; G1 = G2 = G3 = 8-1010 Pa; hL, hL, hL = 3200 m/s; «0 = [0; 7] rad/s;

ê Vladimirs. Litvinenko, Mikhail V. Dvoynikov

Methodology for determining the parameters of drilling mode.

^ = 0.1; ^ = 0.2; ^ = 0.3;

Li = 1800 m; L2 = 190 m; L3 = 10 m; F1 = 1.018-10-3 m2; F2 = 1.8M0-3 m2; F3 = 8.042-10-4 m2;

k = 0.106; E1 = E2 = E3 = 2-1010 Pa; xL, h2, = 5320 m/s; n = [0; 7] rad/s;

^ = 0.1; ^ = 0.2; ^3 = 0.3.

Results of mathematical modeling are presented in fig.3. Comparison of obtained study results for SDM frame oscillations in bench conditions with calculated values of boundaries of DS self-oscillations allows determining the range of stable operation for SDM - DS system. Values located under the line indicated by lower boundary of PB self-oscillations mean absence of vibration - uniform translational and rotational movement of tool, between upper PH and lower PB boundaries - a temporary stop (jamming), above the upper PH - braking (no rotation).

Developed methodology of determining required parameters for drilling mode of inclined sections in a well, ensuring stable operation of BHA, is as follows.

SDM is started and pressure drop is determined during its operation in idle mode. Then, required load on the bit (according to work plan and geological and technical schedule) is created and pressure drop is fixed taking into account loading of gerotor mechanism. On the basis of SDM test bench diagram, optimal range of shaft rotation frequency with corresponding pressure drop is graphically determined. At the same time, maximum allowable decrease in rotation frequency of SDM shaft is noted, which corresponds to optimum amplitudes of frame lateral oscillations.

According to developed mathematical model, boundaries of DS self-oscillations' onset are calculated. After constructing the graphical dependencies, required frequency and load on the bit are determined, at which DS is in permissible range of stable operation. Noting modes of DS stable operation, correlation is made with load on the bit, at which SDM will also be in mode of optimal energy characteristics. If rotation frequency of SDM shaft (according to test bench diagram), determined by pressure drop, has decreased by more than 70 %, load on the bit is reduced. Based on graphical dependences (Fig.3) for boundaries range of self-oscillations' onset at given rotation frequencies of DS and load on the bit, rotation frequency of top drive is adjusted to ensure stable operation of the system while maintaining mechanical drilling speed [14, 15].

Conclusion. On the basis of experimental study, amplitudes and frequencies of oscillations of SDM frame were determined along entire length of power section of working bodies and spindle part for various operating modes of motor. It has been established that to reduce amplitude of motor's lateral oscillations and ensure its stable operation, range of shaft rotation frequencies must be maintained within 70 % of SDM rotation frequency in idle mode.

P, kN 430 ■

348 ■

265 -

183 -

101 -

19 0

P, kN 5370 ■

4320 ■

3269 -

2219 -

1168 -

117

1.4 2.8 4.2 5.6 7 n0, rad/s

0

1.4 2.8 4.2 5.6

7 n, rad/s

Fig.3. Boundaries of rotational (a) and translational self-oscillations (b) of SDM - DS system

1 - Ph; 2 - PB

b

a

1

1

ê Vladimirs. Litvinenko, Mikhail V. Dvoynikov

Methodology for determining the parameters of drilling mode.

Mathematical model of SDM - DS system has been developed, which allows predicting the range of DS self-oscillations and boundaries of rotational and translational wave disturbances for case of string modeling as a heterogeneous rod at drilling directionally straight sections of well.

Methodology has been developed for determining required parameters of drilling mode for directional straight sections of well that ensure stable operation of BHA, based on conditions of maintaining stable operation of system, taking into account maximum permissible rotational frequency of SDM and boundaries of DS self-oscillations' onset.

The developed methodology and technical recommendations aimed at ensuring stable operation of SDM with simultaneous rotation of drilling string at drilling directional wells are used in the branch of "LUKOIL-Inzhiniring" LLC - "KogalymNIPIneft".

REFERENCES

1. Dvoynikov M.V. Research on technical and technological parameters of inclined drilling. Zapiski Gornogo instituta. 2017. Vol. 223, p. 86-92. DOI: 10.18454/PML2017.1.86

2. Dvoynikov M.V., Blinov P.A., Morozov V.A. Optimization of drilling parameters for stabilization sections of zenith angle. Izvestiya vysshikh uchebnykh zavedenii. Ne/t' i gaz. 2016. N 6, p. 65-71 (in Russian).

3. Dvoynikov M.V. Designing of well trajectory for efficient drilling by rotary controlled systems. Zapiski Gornogo instituta. 2018. Vol. 231, p. 254-262. DOI: 10.25515/PMI.2018.3.254

4. Dvoynikov M.V., Muraev Yu.D. Technical and technological solutions ensuring stable operation of a screw downhole motor. Zapiski Gornogo instituta. 2016. Vol. 218, p 198-205 (in Russian).

5. Dvoynikov M.V., Blinov P.A., Morozov V.A., Kunshin A.A. Study of dynamics for drilling string and screw downhole motor. Vestnik assotsiatsii burovykh podryadchikov. 2016. N 2, p. 8-12 (in Russian).

6. Morozov V.A., Dvoynikov M.V. Justification for selection of directional drilling mode parameters using screw downhole motors. Stroite'stvo ne/tyanykh i gazovykh skvazhin na sushe i na more. 2019. N 2, p. 15-18. DOI: 10.30713/0130-3872-2019-2-1518 (in Russian).

7. Morozov V.A., Dvoynikov M.V., Blinov P.A. Study of optimal range of stable operation for "drilling bit - screw downhole motor - drilling string" system. Ne/tegazovoe delo. 2018. Vol. 16. N 2, p. 35-43 (in Russian).

8. Yunin E.K., Khegai V.K. Deep Drilling Dynamics. Moscow: OOO «Nedra-Biznestsentr», 2004, p. 286 (in Russian).

9. Liu X.H., Liu Y.H., Feng D. Downhole Propulsion/Steering Mechanism for Wellbore Trajectory Control in Directional Drilling. Applied Mechanics and Materials. 2013. Vol. 318, p. 185-190.

10. Larsen Lena Kyrvestad. Tools and Techniques to Minimize Shock and Vibration to the Bottom Hole Assembly. University of Stavanger. 2014, p. 135. URL: https://brage.bibsys.no/xmlui/handle/11250/220906

11. Leine R.I., van Campen D.H., Keultjes W.J.G. Stick-slip whirl interaction in drillstring dynamics. ASME Journal of Vibration and Acoustics. 2002. Vol. 124(2), p. 209-220.

12. Li Z., Guo B. Analysis of longitudinal vibration of drillstring in air and gas drilling. Rocky Mountain Oil and Gas Technology Symposium. Denver. Colorado. SPE 107697-MS. Society of Petroleum Engineers. 2007. DOI: org/10.2128/107697-MS

13. Samuel R., Robertson J.E. Vibration Analysis and Control with Hole-Enlarging Tools. Annual Technical Conference and Exhibition. Florence. Italy. SPE 134512-MS. Society of Petroleum Engineers. 2010. DOI: org/10.2118/134512-MS

14. Vromen T.G.M. Control of stick-slip vibrations in drilling systems: PhD thesis. Eindhoven University of Technology. Eindhoven, 2015, p. 256.

15. Zhu Xiaohua, Liping Tang, Qiming Yang. A literature review of approaches for stick-slip vibration suppression in oilwell drillstring. Advances in Mechanical Engineering. 2014. N 6. DOI: 10.1155/2014/967952

Authors: Vladimir S. Litvinenko, Doctor of Engineering Sciences, Rector, rectorat@spmi.ru (Saint Petersburg Mining University, Saint Petersburg, Russia), Mikhail V. Dvoynikov, Doctor of Engineering Sciences, Head of Department, dvoynikov@spmi.ru (Saint Petersburg Mining University, Saint Petersburg, Russia). The paper was received on 30 October, 2019. The paper was accepted for publication on 17 December, 2019.