An improved longitudinal vibration model and dynamic characteristic of sucker rod string

Mingming Xing1 , Shimin Dong2

1, 2School of Mechanical Engineering, Yanshan University, Qinhuangdao 066004, China

2Corresponding author

Journal of Vibroengineering, Vol. 16, Issue 7, 2014, p. 3432-3448.
Received 4 August 2014; received in revised form 15 September 2014; accepted 4 October 2014; published 15 November 2014

Copyright © 2014 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

Considering the influence of the nonlinear characteristics of plunger load and the friction of sucker rod string (SRS) on the SRS’s longitudinal vibration, an improved simulation model of SRS’s longitudinal vibration is derived. In the details, based on the flow characteristic of non-Newtonian power law fluid (NNPLF), a velocity model of NNPLF between pump plunger and pump barrel is established. Then the law of the velocity distribution is solved out with Lagrange multiplier method. Therefore, with the law of the velocity distribution of NNPLF, the computing models of nonlinear friction of pump plunger and clearance leakage between pump plunger and barrel are derived. Taking account of the influence of some parameters on the plunger load, such as plunger friction, hydraulic loss of pump and clearance leakage, an improved simulation model of plunger load is derived. The dynamic response is solved out with fourth order Runge-Kutta method. Comparing experiment results with simulated results, good agreement is found, which shows the simulation model is feasible. The influences of the different parameters on pump pressure and pump plunger load are analyzed, such as stroke number, power law exponent, consistency coefficient and gap between plunger and pump barrel. Simulation result indicates that the opening time of standing valve and traveling valve is affected by the parameters, and the maximum and minimum loads of pump plunger are affected by stroke number. In addition, the influence of SRS absorber on SRS’s longitudinal vibration is analyzed.

Keywords: sucker rod string, longitudinal vibration, non-Newtonian power law fluid, pump plunger load, pump pressure.

1. Introduction

The sucker rod pumping system is one of the mechanical recovery methods in the world’s oil industry and it is by far the predominant artificial lift method used in the world for producing oil well [1-3]. And SRS is a slender rod of several kilometers, which transmits the movement from mechanical rotary motion to plunger linear motion. Meanwhile, the characteristic analysis of SRS’s longitudinal vibration is foundation of system fault diagnosis and oil production forecast. Therefore, the dynamic simulation of SRS draws more and more attention from domestic and foreign scholars. One of the most relevant and successful simulation model describing the sucker rod dynamics is developed by Gibbs [4], then Gibbs [5-8] extends SRS dynamics to computer-based diagnostics of the down-hole well conditions. An improved model of SRPS is built by Dale Russel Doty and Zelimir Schmidt et al. [9], incorporating the dynamics of the liquid columns as well as SRS through a system of partial differential equations. Subsequently, based on the improved model, the dynamic simulation of SRPS is carried out by S. D. L. Lekia and H. A. Tripp et al., and the model is verified by comparing it with the experimental data [10-12]. Then Lukasiewicz, S. A. and J. Xu, et al. [13, 14] apply the dynamic simulation technique of SRPS to a slant-hole or horizontal well, and build a longitudinal vibration model of SRS in a horizontal well. Recently, DaCunha, I. N. Shardakov and Luan Guo-hua [1, 3, 15] analyse boundary conditions effecting rod law of motion and correct the excitation forces of rod longitudinal vibration. In the detail, The SRS’s longitudinal vibration is described with the semi-infinite-spatial-domain wave equation by DaCunha. Shardakov thinks that the multivalued force-velocity relation arises due to the operation of pump valves and Coulomb’s friction, the numerical model of SRS’s longitudinal vibration is built. A new one-dimensional vibration model of SRS for deep well is built by G. H. Luan, considering the structure of the improved side-flow pump.

However, the most of aforementioned approaches have not been studied sufficiently. In the detail, the empirical formula or simplified formula is usually used to build the simulation model of plunger load [16-18]. Obviously, it is hardly possible for empirical formula to meet the precision requirement of plunger load. And the simplified formula is obtained based on the specific oil wells, which does not suit the other wells. For example, there are two methods for computing the friction of pump plunger. The experiment is carried out with water instead of lubricant [19], and the empirical formula is given [17, 19, 20], as follows:

(1)
F p f = 0.94 D p δ - 140 ,

where, Fpf is the liquid friction of pump plunger (N); Dp is the diameter of pump plunger (mm); δ is the gap between pump plunger and pump barrel (mm).

The Eq. (1) was given under some assumptions by Wan B. J. [19]. Wan thought that the produced fluid in oil-water was the unstable emulsion when the oil well had high water content such as more than 50 %. In the oil well, the liquid was easily layered into crude oil and water, the crude oil was dispersed phase and the water is continuous phase. Therefore, Wan thought that the liquid was mainly water between pump plunger and pump barrel. Based on the above analysis, the lubricant of half dry friction of plunger could be replaced with water [19].

As a result, the above formula is obtained under the high water content condition, which only applies to the water-flood or thin oil well. However, when the oil well is a polymer flooding or heavy oil well, the liquid of oil well is Non-Newtonian power law fluid, the Eq. (1) will be inappropriate.

With the second method, the liquid is used as Newtonian power law fluid, and the liquid friction of pump plunger is derived based on gap flow theory of hydromechanics [20], as follows:

(2)
F p f = π L p D p Δ P δ 2 L p + μ v p δ 1 1 - ε 2 ,

where, μ is the instantaneous dynamic viscosity of liquid, kg/(s∙m); vp is the velocity of pump plunger (m/s); Lp is the length of pump plunger (m); ΔP is the pressure difference inside and outside of pump (Pa); ε is the eccentric distance (m).

The Eq. (2) was derived in ideal circumstances [20]. For example, the inner surface of pump barrel as well as the outer surface of pump plunger is defined as an ideal cylindrical surface. However, those surfaces are not absolute straight cylindrical surfaces, and the crevices of the surface are not formed by ideal parallel surfaces. In addition, the Eq. (2) is derived under the Newtonian fluid, which will not apply to the oil well of non-Newtonian power law fluid. According to the above analysis, the computing model of plunger friction needs to be derived for the oil well of non-Newtonian power law fluid.

The current simulation models of plunger load contain the model of hydraulic loss [21, 22] and the model of clearance leakage [23, 24], the formulas are as follows:

(3)
Δ P s = ρ g 2 g γ 2 A p A v 2 v P 2 ,
(4)
q v = π δ 3 D p Δ p i o 12 μ L p 1 + 3 ε 2 2 δ 2 ± π δ D p v p 2 ,

where, ρg is the gravity gradient of liquid (N/m3); Ap is the cross-sectional area of pump plunger (m2); Av is the area of valve’s aperture (m2); γ is the flow coefficient.

From Eq. (3), the hydraulic loss of pump is proportional with the speed squared of pump plunger. However, in order to obtain linear solution, the formula of hydraulic loss is usually linearized at present. Generally, the speed of pump plunger is defined as constant with stroke number and frequency. Therefore, the simulation result will not reflect the influence of the instantaneous speed of pump plunger on hydraulic loss of pump. With Eq. (4), the problem of annular pressure difference-shear leakage direction of pump is discussed, but the influence of the NNPLF on the annular pressure is not considered.

According to the analyses above, the simulation models of plunger nonlinear friction, hydraulic loss and clearance leakage need to be improved with the action of NNPLF. Considering the influence of the time-varying characteristics of some dynamic parameters on the plunger load, such as plunger nonlinear friction, hydraulic loss and clearance leakage, the simulation model of plunger load needs to be corrected. Finally, considering the influence of the nonlinear characteristics of plunger load and SRS friction on the SRS’s longitudinal vibration, the simulation model of SRS’s longitudinal vibration needs to be improved.

In this paper, in section 2, an improved simulation model of mix SRS with absorber is built. The computational model of pump plunger load is corrected. In the details, based on the flow characteristic of NNPLF, a new computational model of NNPLF’s velocity between pump plunger and pump barrel is established. Then the law of the velocity distribution is solved out with Lagrange multiplier method. Therefore, the nonlinear friction between NNPLF and SRS as well as the nonlinear friction between NNPLF and pump plunger is derived. Taking account of the influence of nonlinear characteristics of some parameters on the pump plunger load, such as plunger friction, hydraulic loss and clearance leakage, an improved model of the pump load is derived. In section 3, the dynamic response is solved out with fourth order Runge-Kutta method. In section 4, a wireless indicator is used to verify the feasibility of the proposed dynamics simulation method, and the sensitive parameters of plunger load and dynamic characteristic of mix SRS with absorber are analyzed. In section 5, the signification of dynamic simulation research on SRS and the conclusions are summarized.

2. The dynamic simulation model of SRS’s longitudinal vibration

The number stages of SRS assemblage are supposed to be k, and the polished rod displacement is u*. The relative displacement between polished rod and SRS top is u0. Considering the influence of some factors on the SRS longitudinal vibration, such as SRS absorber between polished rod and SRS top, nonlinear friction between SRS and NNPLF and nonlinear plunger load, a mechanical model of SRS’s longitudinal vibration is established, as shown in Fig. 1.

Fig. 1. The mechanical model of SRS’s longitudinal vibration

 The mechanical model of SRS’s longitudinal vibration

Based on the mechanical model of longitudinal vibration, a mathematical model is built, as follows:

(5)
m 1 u ¨ 1 + k 0 u 1 - u * + k 1 u 1 - u 2 + c 0 u ˙ 1 - u ˙ * + c 1 u ˙ 1 - u ˙ 2 n = 0 , m 2 u ¨ 2 + k 1 u 2 - u 1 + k 2 u 2 - u 3 + c 1 u ˙ 2 - u ˙ 1 n + c 2 u ˙ 2 - u ˙ 3 n = 0 , m i u ¨ i + k i - 1 u i - u i - 1 + k i u i - u i + 1 + c i - 1 u ˙ i - u ˙ i - 1 n + c i u ˙ i - u ˙ i + 1 n = 0 , m n - 1 u ¨ n - 1 + k n - 2 u n - 1 - u n - 2 + k n - 1 u n - 1 - u n + c n - 2 u ˙ n - 1 - u ˙ n - 2 n + c n - 1 u ˙ n - 1 - u ˙ n n = 0 , m n u ¨ n + k n - 1 u n - u n - 1 + c n - 1 u ˙ n - u ˙ n - 1 n - F p t = 0 ,

where:

(6)
F p t = P p t t + F p f .

When {u}=u1, u2, , un-1, un T, the equation is obtained as follows:

(7)
M u ¨ + K u + C u ˙ n = F ,

where, ci is the friction coefficient (Pa.s); Fp is the load of pump plunger (N); ki is the stiffness coefficient of the ith element, m2; Ci is the damping coefficient of the ith element, Pa; ui is the displacement of the ith element, m. Pp(t) is the liquid load of pump plunger (kN); Fpf is the friction between pump plunger and pump barrel (kN).

2.1. The nonlinear friction of pump plunger and clearance leakage

A structure schematic of pump plunger and pump barrel is given in Fig. 2. Lp is the length of pump plunger, Rp is the radius of pump plunger, Rpb is the radius of pump barrel, δ is the gap between pump plunger and pump barrel.

Fig. 2. The structure schematic of pump plunger and pump barrel

 The structure schematic of pump plunger and pump barrel

The liquid of oil well is supposed to be NNPLF with a homogeneous incompressible fluid and the oil well is a vertical well. The pump plunger and pump barrel are defined as rigid body with no elastic deformation and the pump barrel is static.

Based on the constitutive equation and shearing stress equation of NNPLF, the formula is given, as follows:

(8)
τ p l = K p l d v p d r n , d d r r τ p l = r d p d z .

The radial velocity gradient of NNPLF between pump plunger and pump barrel is derived based on the above formula:

(9)
d v p d r = f 2 K p l r + 2 C r 1 n ,

where:

(10)
f = p z ,

where, C is the undetermined parameters; n is the power law exponent; kpl is the consistency coefficient; f is the axial pressure gradient of NNPLF in the annulus.

The velocity of NNPLF in annulus is influenced by pressure difference flow and pulsatile flow, the velocity distribution curve is given in Fig. 3. In Fig. 3, R0 is the distance from a point with zero velocity of NNPLF to the centerline of pump plunger (m); Rm is the distance from a point with maximum velocity of NNPLF to the centerline of pump plunger (m); Rp is the radius of pump plunger (m); Rpb is the radius of pump barrel (m).

Fig. 3. The velocity distribution curve of NNPLF

 The velocity distribution curve of NNPLF

When r=Rm and dvp/dr=0, then:

(11)
d v p d r = f 2 K p l r - R m 2 r 1 n .

Therefore, the velocity distribution law of NNPLF and the boundary conditions are obtained:

(12)
v r , t = r R p b f 2 K p l r - R m 2 r 1 n d r , R p R 0 f 2 K p l r - R m 2 r 1 n d r = d u p l d t , 2 π R p R p b v r d r = π R p l 2 d u p l d t , R 0 R m f 2 K p l r - R m 2 r 1 n d r = R m R p b f 2 K p l r - R m 2 r 1 n d r .

The analytical formulas of the parameters are not derived, such as R0, Rm, f. So the Lagrange multiplier method is used to find the numerical solution. Firstly, the variables are defined, as follow:

(13)
x ^ = x 1       x 2         x 3 = R 0         R m         f .

Combining Eqs. (12) and (13), the constraint and the objective functions are derived:

(14)
m i n F x ^ = R p R p b f 2 K p l r - R m 2 r 1 n d r - v p l , g 1 = R p b - x 2 0 , g 2 = x 2 - x 1 0 , g 3 = x 1 - R p 0 , g 4 = x 2 R p b x 3 2 K p l r - x 2 2 r 1 n d r 0 , h 1 = 2 π R p R p b r r R p b x 3 2 K p l x - x 2 2 x 1 n d x d r - π R p l 2 v p l = 0 , h 2 = R p x 1 x 3 2 K p l r - x 2 2 r 1 n d r - v p l = 0 , h 3 = x 1 x 2 x 3 2 K p l r - x 2 2 r 1 n d r - x 2 R p b x 3 2 K p l r - x 2 2 r 1 n d r = 0 .

The parameters R0, Rm, f are solved out based the Lagrange multiplier method, and the velocity distribution law is expressed fully. Therefore, the friction of pump plunger is given as follows:

(15)
F p f = 2 π R p L p K p l d v p d r r = R p n .

Eq. (15) can be defined as follows:

(16)
F p f = λ p l v p n ,

where:

(17)
λ p l = 2 π K p l L p R p 2 - R m 2 R p R 0 r - R m 2 r 1 n d r n   ,

where, λpl is the resistance coefficient of pump plunger, N·s/mn/m.

The computing method of the nonlinear friction between SRS and NNPLF is the same as for the nonlinear friction between pump plunger and pump barrel, so we will not repeat it here.

Combining Eq. (11), the computing model of clearance leakage between pump plunger and pump plunger is derived, as follows:

(18)
q v m = π D p n 1 + 2 n Δ p i o k p l L p 1 n R m - R p 1 + 2 n n +   R p b - R m 1 + 2 n n - π D p v p R m - R p .

Combining Eqs. (11) and (17), when n= 1, μ=kpl, the equation of velocity of pump plunger is obtained, as follows:

(19)
v p = 1 2 Δ p i o μ L p R m - R p 2 - R p b - R m 2 .

Combining Eqs. (18) and (19), when n= 1, μ=kpl, δ=Rpb-Rp, the computing model of clearance leakage of Newtonian fluid under laminar flow is derived:

(20)
q v = π δ 3 D p Δ p i o 12 μ L p - π δ D p v p 2 .

When ε= 0, the Eq. (4) is the same as above Eq. (20).

2.2. The computational model of liquid load of pump plunger

Currently, the computational model of liquid load of pump plunger is given as follows:

(21)
P p t = A p p d - p - A L r p d ,

where, Ap is the cross area of pump plunger (m2), ALr is the cross area of SRS bottom (m2), pd is the discharge pressure (Pa), p is the instantaneous liquid pressure in pump (Pa).

From Eq. (21), the instantaneous liquid pressure is a major factor which affects liquid load of pump plunger. The current model of instantaneous liquid pressure is built based on the Newtonian fluid. In this paper, a new model of instantaneous liquid pressure is built based on the NNPLF. The liquid pressure law in pump is given in Fig. 4 [25]. In Fig. 4, Vx is the instantaneous volume of movement distance of pump plunger (m3), Vp is the effective volume of movement distance of pump plunger (m3), V0 is the clearance volume of pump (m3).

Fig. 4. The liquid pressure law in pump

 The liquid pressure law in pump

2.2.1. Liquid flowing into pump in the up stroke

When the standing valve is open, the liquid is drawn into the pump. At this point, the instantaneous liquid pressure of pump is influenced by the submergence pressure and the hydraulic loss. The formula of instantaneous liquid pressure is derived:

(22)
P = P s - Δ P s ,

where:

(23)
Δ p i = n v ρ 1 γ 2 d u p d t 2 A p A v 2 3 n + 1 2 n + 1 n + 3 2 5 n + 1 A v A p 2 - A v A p + 3 2 n + 1 2 5 n + 1 ,      
      i = s , d ,

where, ΔPs is the part hydraulic loss of standing valve (Pa); up is the instantaneous displacement of pump plunger (m); ρ is the density of liquid (kg/m3). nv is the number of traveling valve.

2.2.2. Fluid leaving pump in the down stroke

When the traveling valve is open, the liquid is discharged from the pump. At this point, the instantaneous liquid pressure of pump is influenced by the discharge pressure and the hydraulic loss. The formula of instantaneous liquid pressure is derived:

(24)
P = P d + Δ P d ,

where, ΔPd is the part hydraulic loss of traveling valve (Pa), and the formula of hydraulic loss of traveling valve is Eq. (23).

2.2.3. Compressing phase of gas in the down stroke

When the standing valve is close and before the traveling valve is open, the pump plunger moves from top dead to down. Then the gas is compressed in the pump. Based on the polytropic process of gas, the formula of instantaneous liquid pressure is derived:

(25)
p = V g V g - V p + V x - Δ V l d n q P s - Δ P s ,

where:

(26)
V g = r a 1 + r a V 0 + V p , V x = A p u p - A p Δ P L t E t A t ,

where, Vg is the volume of free gas in the pump (m3); Ps is the submergence press of pump (Pa); ΔVl is the volume of clearance leakage between pump plunger and pump barrel (m3); ΔP is the pressure difference between inner and outer pump (Pa); ra is the ratio of gas and liquid (m3/m3); Dpb is the diameter of pump barrel (m); μ is the dynamic viscosity, kg/(s∙m); Lp is the length of pump plunger (m); up is the instantaneous displacement of pump plunger (m); Lt is the length of tube (m); Et is the elasticity modulus of tube (N/m2); At is the cross area of tube (m3).

When the traveling valve is opening, the time is defined as ttv, the volume of movement distance of pump plunger is defined as Vtv, the volume of clearance leakage is defined as ΔVl(ttv). The the equilibrium equation is obtained as follows:

(27)
V g - V p + V t v - Δ V l d t t v = V g P s - Δ P s P d + Δ P d 1 n q .

So, the gas of the clearance volume is derived as follows:

(28)
V o g = V 0 V 0 + V t v V g - V p + V t v - Δ V l d t t v .

2.2.4. Expanding phase of gas in the up stroke

When the traveling valve is close and before the standing valve is open, the pump plunger moves from bottom dead to up. Then the gas is expanded in the pump. With polytropic process of gas, the formula of instantaneous liquid pressure is derived:

(29)
p = V o g V o g + V x - Δ V l s n q P d + Δ P d .

When the standing valve is opening, the time is defined as tsv, the volume of movement distance of pump plunger is defined as Vsv, the volume of clearance leakage is defined as ΔVl(tsv). Then the equilibrium equation is obtained as follows:

(30)
V o g + V s v - Δ V l s t s v = V o g P d + Δ P d P s - Δ P s 1 n q .

3. The numerical model of SRS’s longitudinal vibration

The system equations of nonlinear longitudinal vibration of SRS are solved by fourth-order Runge-Kutta method, and the numerical model of longitudinal vibration of SRS is obtained. In this paper, the SRS is decomposed into five hundred elements, so the element length is four meters. Based on the basic parameters of beam pump system, the dynamic responses of system are obtained. Comparing the plunger displacements of two contiguous periods, if good agreement is found between the two simulation results within the allowed error range, the computing of the corresponding program will be stopped, and the dynamic response will be obtained:

(31)
d y 1 = y 2 , d y 2 = 1 m 1 k 0 t y 1 - u * - k 0 t + k 1 t y 1 + k 1 t y 3 + c 0 y 2 - u ˙ * - c 0 + c 1 y 2 n + c 1 y 4 n , d y 3 = y 4 , d y 4 = 1 m 2 k 1 t y 1 - k 1 t + k 2 t y 3 + k 2 t y 5 + c 1 y 2 n - c 1 + c 2 y 4 n + c 2 y 6 n , d y 2 n - 3 = y n - 2 , d y n - 2 = 1 m n - 1 k n - 2 t y 2 n - 5 - k n - 2 t + k n - 1 t y 2 n - 3 + k n - 1 t y 2 n - 1 + c n - 2 y 2 n - 4 n - c n - 2 + c n - 1 y 2 n - 2 n + c n - 1 y 2 n n , d y 2 n - 1 = y 2 n , d y 2 n = 1 m n k n - 1 t y 2 n - 3 - k n - 1 t y 2 n - 1 + c n - 1 y 2 n - 2 n - c n - 1 y 2 n n + P p t ,

where:

(32)
y 1 = u 1 , y 2 = u ˙ 1 , y 3 = u 2 , y 4 = u ˙ 2 ,                                                       y 2 n - 3 = u n - 1 , y 2 n - 2 = u ˙ n - 1 , y 2 n - 1 = u n , y 2 n = u ˙ n .

In addition, for the first calculation, the initial conditions are the same for all the plots as given below:

(33)
y j = 0 ,         j = 1 ,   2 ,   3 , ,   2 n .

4. Experiment and discussion

The basic parameters are as follows: the well number is #62-39; the outdoor temperature during testing Tho is 28°C; the pumping unit type is YCYJ12-5-89HB; the motor type is YCCH250-8; the rated power of the motor with high rotating variation PH is 30 kW; the length of belt L is 4 m; the linear density of single belt ρ is 0.37 kg/m; the cross-sectional area of belt A is 1.55×10-4 m2; the coefficient of friction of the belt against the pulleys μ is 0.6; the diameter of the small pulley installed on the motor output shaft D1 is 0.3 m; the diameter of the big pulley installed on the gearbox input shaft D2 is 0.9 m; the efficiency of the gearbox system η is 0.95; the transmission ratio of the gearbox system ibox is 35; the Young’s modulus of SRS E1 as well as E2 is 2.1×1011 Pa; the stroke number n is 4 min-1; the stroke length S is 5 m; the density of the rod ρ1 as well as ρ2 is 7850 kg/m3; the outside diameter of tube Dt is 73 mm; the SRS is 22×1200 m+19×800 m; the diameter of the pump Dp is 38 mm; the pump setting depth Lpd is 2000 m; the depth of dynamic liquid level Hdl is 1600 m; the oil pressure Po is 0.3 MPa; the casing pressure Pc is 0.5 MPa; the ratio of gas and oil Rgo is 50 % m3/m3; the density of crude oil; ρo is 870 kg/m3; the middle depth of oil Lz is 2500 m.

4.1. Precision verification

A wireless indicator is installed between polished rod and SRS, in Fig. 5. With the indicator, the surface dynamometer card is obtained in Fig. 6. Based on the simulation model of surface driving units [2] and the improved simulation model in this paper, the dynamic response of beam pump system is obtained with above basic parameters. In the detail, the simulation model of two degree of freedom torsional vibration for surface driving system is built, and the equivalent rotational inertia of two discs are derived with the masses and rotational inertia of every moving part for surface driving system [2]. Combining the simulation models of surface and down-hole system, the numerical model of beam pump system is solved with Runge-Kutta method. And the surface dynamometer card is obtained in Fig. 6.

According to the Fig. 6, the maximum and minimum polished rod loads of the measured curve are 78.14 kN, 38.71 kN, and the maximum and minimum polished rod loads of simulation curve are 76.01 kN, 37.13 kN. The relative errors of measured and simulate curves are 2.73 %, 4.08 %. In the up and down stroke, the maximum and minimum relative errors between the measured curve of the polished rod load and the simulation curve are 4.86 %, 0.71 %. The simulation curves of the dynamometer card based on the new model are found to be in good agreement with the experimental curves. Therefore, the new model in this paper is accurate enough to be used for engineering practice.

Fig. 5. The wireless indicator

 The wireless indicator

Fig. 6. Comparison of the numerical solution with experimental load–displacement curve

 Comparison of the numerical solution with experimental load–displacement curve

4.2. Affecting factors of plunger load

4.2.1. Hydraulic loss

When the clearance leakage between pump plunger and pump barrel Qt= 0, the influence of the hydraulic loss on the plunger load is only considered and the comparison curves of pump plunger load with no hydraulic loss or hydraulic loss are given in Fig. 7(a). In Fig. 7(b), the comparison curves of pump plunger load with current model or new model are given. When other parameters do not change, considering the influence of stroke number and power law exponent on hydraulic loss, the curves of pump plunger load with different stroke number are given in Fig. 7(c), and the curves of pump plunger load with different power law exponents are given in Fig. 7(d).

Fig. 7. The influence of hydraulic loss on plunger load

 The influence of hydraulic loss on plunger load

a) Influence of hydraulic loss

 The influence of hydraulic loss on plunger load

b) Comparison of current and new model

 The influence of hydraulic loss on plunger load

c) Stroke number

 The influence of hydraulic loss on plunger load

d) Power law exponent

As you can see in the above figures, when the hydraulic loss is taken into consideration, the maximum load of pump plunger increases and the minimum load of pump plunger decreases. With increasing stroke number and power law exponent, the maximum load of pump plunger is increasing and the minimum load of pump plunger is decreasing. In addition, the area of pump diagram is decreasing with stroke number and power law exponent increasing.

Fig. 8. The influence of clearance leakage on plunger load

 The influence of clearance leakage on plunger load

a) Influence of clearance leakage

 The influence of clearance leakage on plunger load

b) Comparison of current and new model

 The influence of clearance leakage on plunger load

c) Stroke number

 The influence of clearance leakage on plunger load

d) Power law exponent

 The influence of clearance leakage on plunger load

e) Gap between plunger and pump barrel

 The influence of clearance leakage on plunger load

f) Consistency coefficient

4.2.2. Clearance leakage

When the hydraulic loss of pump valve ΔP= 0, the influence of the clearance loss between pump plunger and pump barrel is only considered and the comparison curves of pump plunger load with no clearance leakage or clearance leakage are given in Fig. 8(a). In Fig. 8(b), the comparison curves of pump plunger load with current model or new model are given. When other parameters do not change, considering the influence of some parameters on clearance leakage, the curves of plunger load with different parameters are given, such as stroke number, power law exponent, consistency coefficient and gap between plunger and pump barrel, in Fig. 8(c), (d), (e), (f).

In Fig. 8, when the clearance leakage between pump plunger and pump barrel is taken into consideration, the falling rate of liquid pressure is decreasing in the up stroke. Therefore, the standing valve opens with short delay. On the other hand, the increasing rate of liquid pressure is increasing in the down stroke, so the traveling valve opens in advance. With stroke number, power law exponent and consistency coefficient decreasing, the lagged time of standing valve opening is increasing and the lead time of traveling valve opening is decreasing. In addition, with the gap between pump plunger and pump barrel increasing, the lagged time of standing valve opening is increasing and the lead time of traveling valve opening is decreasing.

Fig. 9. The influence of hydraulic loss and clearance leakage on plunger load

 The influence of hydraulic loss and clearance leakage on plunger load

a) Influence of hydraulic loss and clearance leakage

 The influence of hydraulic loss and clearance leakage on plunger load

b) Comparison of current and new model

 The influence of hydraulic loss and clearance leakage on plunger load

c) Stroke number

 The influence of hydraulic loss and clearance leakage on plunger load

d) Power law exponent

 The influence of hydraulic loss and clearance leakage on plunger load

e) Gap between plunger and pump barrel

 The influence of hydraulic loss and clearance leakage on plunger load

f) Consistency coefficient

4.2.3. Hydraulic loss and clearance leakage

The influence of hydraulic loss of pump valve on plunger load as well as the influence of clearance leakage between pump plunger and pump barrel is considered, the comparison curves of plunger load with no hydraulic loss or clearance leakage are given in Fig. 9(a). In Fig. 9(b), the comparison curves of pump plunger load with current model or new model are given. When other parameters do not change, with the influence of some parameters on plunger load, the curves of plunger load with different parameters are given, such as stroke number, power law exponent, consistency coefficient and gap between plunger and pump barrel, in Fig. 9(c), (d), (e), (f). As see in Fig. 9, considering the influence of hydraulic loss and clearance leakage on plunger load, the standing valve opens with short delay and the maximum load of pump plunger increases in the up stroke. In the down stroke, the traveling valve opens in advance and the minimum load of pump plunger decreases. With the stroke number decreasing, the lagged time of standing valve opening is increasing and the maximum load is decreasing in the up stroke, the lead time of traveling valve opening is decreasing and the minimum load is increasing in the down stroke. In addition, with the power law exponent and consistency coefficient decreasing, the lagged time of standing valve opening is increasing in the up stroke, the lead time of traveling valve opening is decreasing in the down stroke. When the gap between pump plunger and pump barrel is decreasing, the lagged time of standing valve opening is decreasing in the up stroke, the lead time of traveling valve opening is increasing in the down stroke.

4.3. Impact analysis of SRS absorber

When c0= 6.0×105 N·s/m, k0= 2.0×106 N/m, the displacement of pump plunger and load of polished rod are simulated. In addition, the plunger displacement and polished rod load are also obtained with no SRS absorber. The relative curves of displacement and load are given in Fig. 10. When c0= 0 and other parameters do not change, the excitation frequency and stiffness of SRS absorber are changed simultaneously, the three-dimensional curves of plunger displacement and polished rod load are obtained in Fig. 11. When c0= 6.0×105 N·s/m and other parameters do not change, the excitation frequency and stiffness of SRS absorber are changed simultaneously, the three-dimensional curves of plunger displacement and polished rod load are obtained in Fig. 12.

Fig. 10. The amplitude-frequency curve

 The amplitude-frequency curve

a) Pump plunger displacement

 The amplitude-frequency curve

b) Polished rod load

As you can see in figures, when the SRS absorber is ignored, there are two resonance points in the frequency range of 0 to 20 rad/s. But there are there resonance points in same frequency range with SRS absorber, which is a major characteristic of nonlinear vibration. In addition, the polished rod load decreases in non-resonant region when the SRS absorber is considered. The influences of excitation frequency and absorber stiffness on the plunger displacement and polished rod load are intuitively expressed in three-dimensional curves. The above conclusions are very important to the design of the SRS absorber in practical applications.

Fig. 11. The amplitude-frequency with no damped absorber

 The amplitude-frequency with no damped absorber

a) Pump plunger displacement

 The amplitude-frequency with no damped absorber

b) Polished rod load

Fig. 12. The amplitude-frequency with damped absorber

 The amplitude-frequency with damped absorber

a) Pump plunger displacement

 The amplitude-frequency with damped absorber

b) Polished rod load

5. Conclusions

The following conclusions can be summarized from the theoretical formulation and the numerical studies in this paper:

1) Considering the influence of nonlinear characteristic of SRS friction and plunger load on the SRS’s longitudinal vibration, an improved simulation model of longitudinal vibration for mix absorber with SRS is built. In the detail, selecting NNPLF as the research object, some new models are established, such as the nonlinear friction between pump plunger and NNPLF, the hydraulic loss of pump and the clearance leakage between pump plunger and pump barrel. The dynamic responses are obtained with zero initial conditions based on fourth order Ronge-Kutta method.

2) In order to verify the feasibility of the proposed dynamics simulation model, the wireless indicator is used. It shows that the simulated and experimental results of polished rod load have a good match.

3) Considering the influence of hydraulic loss and clearance leakage on the instantaneous liquid pressure of pump, the influences of some parameters on the pump plunger load are analyzed, such as stroke number, power law exponent, consistency coefficient and the gap between pump plunger and pump barrel. The results show that when the stroke number is increasing, the maximum load of pump plunger is increasing and minimum load of pump plunger is decreasing. With the stroke number, the power law exponent and consistency coefficient decreasing, the lagged time of standing valve opening is increasing in the up stroke, the lead time of traveling valve opening is decreasing in the down stroke. In addition, when the influence of SRS absorber on SRS’s longitudinal vibration is taken into consideration, there is an obvious characteristic of nonlinear longitudinal vibration. The above conclusions are very important to system fault diagnosis and the design of the SRS absorber in practical applications.

Acknowledgements

This project is supported by National Natural Science Foundation of China (No. 50974108) and National Natural Science Foundation of China (No. 51174175).

References

  1. Luan G. H., He Sh. L., Yang Zh. A prediction model for a new deep-rod pumping system. Journal of Petroleum Science and Engineering, Vol. 80, Issue 1, 2011, p. 75-80. [CrossRef]
  2. Xu G. M., Zhou F., Wu H. X. Dynamic coupled model for oil-pumping system and electrical machine matching energy-save simulation. Journal of Harbin Institute of Technology, Vol. 43, Issue 5, 2011, p. 126-130, (in Chinese). [CrossRef]
  3. Shardakov I. N., Wasserman I. N. Numerical modelling of longitudinal vibrations of a sucker rod string. Journal of Sound and Vibration, Vol. 329. Issue 3, 2010, p. 317-327. [CrossRef]
  4. Gibbs S. G. Predicting the behavior of sucker rod pumping systems. Journal of Petroleum Technology, Vol. 15, Issue 7, 1963, p. 769-778. [CrossRef]
  5. Gibbs S. G., Neely A. B. Computer diagnosis of down-hole conditions in sucker rod pumping wells. Journal of Petroleum Technology, Vol. 18, Issue 1, 1966, p. 91-98. [CrossRef]
  6. Gibbs S. G. A Review of methods for design and analysis of rod pumping installations. Journal of Petroleum Technology, Vol. 34, Issue 12, 1982, p. 2931-2940. [CrossRef]
  7. Gibbs S. G. Design and diagnosis of deviated tod-pumped wells. Journal of Petroleum Technology, 1992, p. 774-781. [CrossRef]
  8. Everitt T. A., Jennings J. W., Texas A. An improved finite-difference calculation of downhole dynamometer cards for sucker-rod pumps. SPE Production Engineering, 1992. [CrossRef]
  9. Doty D. R., Schmidt Z. An improved model for sucker rod pumping. Society of Petroleum Engineers Journal, Vol. 23, Issue 1, 1983, p. 33-41. [CrossRef]
  10. Lekia S. D. A coupled rod and fluid dynamic model for predicting the behavior of sucker-rod pumping systems part 1: model theory and solution methodology. SPE Production and Facilities, Vol. 10, Issue 1, 1995, p. 26-33. [CrossRef]
  11. Lekia S. D. A coupled rod and fluid dynamic model for predicting the behavior of sucker-rod pumping systems – part 2: parametric study and demonstration of model capabilities. SPE Production and Facilities, Vol. 10, Issue 1, 1995, p. 34-40. [CrossRef]
  12. Tripp H. A., Kilgore J. J. A Comparison between predicted and measured walking beam pump parameters. SPE Technical Conference and Exhibition held in New Orleans, 1990. [CrossRef]
  13. Lukasiewicz S. A. Dynamic Behavior of the sucker rod string in the inclined well. SPE Production Operations Symposium held in Oklahoma City, 1991. [CrossRef]
  14. Xu J. A method for diagnosing the performance of sucker rod string in straight inclined wells. American/Caribbean Pelroleum Engineering Conference held In Buenos AIres, 1994. [CrossRef]
  15. DaCunha J. J., Gibbs. Modeling a finite-length sucker rod using the semi-infinite wave equation and a proof to Gibbs’ conjecture. Society of Petroleum Engineers Journal, Vol. 14, Issue 1, 2009, p. 112-119. [CrossRef]
  16. Podio A. L., Gomez J., Mansure A. J. Laboratory-instruments sucker-rod pump. SPE Production and Facilities, Vol. 18, Issue 2, 2003, p. 104-113. [CrossRef]
  17. Liu X. F., Qi Y. G. A modern approach to the selection of sucker rod pumping systems in CBM wells. Journal of Petroleum Science and Engineering, Vol. 76, Issue 3-4, 2011, p. 100-108. [CrossRef]
  18. Dotson B., Nunez-Paclibon E. Gas well liquid loading from the power perspective. SPE Production and Operation, 2011. [CrossRef]
  19. Wan B. J. The Design and Calculation of Mechanical Oil Production. Petroleum Industry Press, 1988, (in Chinese). [CrossRef]
  20. Zhao H. J. A discussion of the friction between pump plunger and barrel. China Petroleum Machinery, Vol. 21, Issue 2, 1993, p. 34-37, (in Chinese). [CrossRef]
  21. Wang Z. B., Yingchuan L., et al. A simple numerical model for the prediction of multiphase mass flow rate through chokes. Petroleum Science and Technology, Vol. 29, 2011, p. 2545-2553. [CrossRef]
  22. Enfis M. S., Ahmed R. M. The hydraulic effect of tool-joint on annular pressure loss. SPE Production Operations Symposium held in Oklahoma City, 2011. [CrossRef]
  23. Wu Y. Q., Wu X. D. Research on the problem of annular pressure difference-shear leakage direction of rod pump. China Petroleum Machinery, Vol. 40, Issue 8, 2012, p. 107-110, (in Chinese). [CrossRef]
  24. Jeong Y. T., Shah S. N. Analysis of tool joint effects for accurate friction pressure loss calculations. IADC/SPE Conference held in Dallas, 2004. [CrossRef]
  25. Gilbert W. E. An oil-well pump dynagraph. American Petroleum Institute (API), 1936. [CrossRef]