Next Article in Journal
Essential Norm of t-Generalized Composition Operators from F(p, q, s) to Iterated Weighted-Type Banach Space
Previous Article in Journal
The Blow-Up of the Local Energy Solution to the Wave Equation with a Nontrivial Boundary Condition
Previous Article in Special Issue
Investigating Behavior of Slider–Crank Mechanisms with Bearing Failures Using Vibration Analysis Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Method for Predicting the Behavior of a Sucker Rod Pumping Unit Based on the Polished Rod Velocity

College of Science, China University of Petroleum (East China), Qingdao 266580, China
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(9), 1318; https://doi.org/10.3390/math12091318
Submission received: 20 March 2024 / Revised: 16 April 2024 / Accepted: 17 April 2024 / Published: 25 April 2024

Abstract

:
Fault dynamometer cards are the basis of the diagnosis technique for sucker rod pumping systems. Predicting fault cards with a pumping condition model is an economical and effective method. The usual model is described by a mixed function of the pump displacement and pump load, and it is difficult to use in the prediction method based on the analytical solution of the sucker rod string wave equation. In this paper, a normal pumping condition model described by a function of polished rod velocity is proposed. For the analytical solution of the sucker rod wave equation, an iterative prediction algorithm with pumping condition models is proposed, its convergence is analyzed, and then it is validated by classical finite difference method simulated cards and measured surface dynamometer cards. The results show that the proposed algorithm is accurate. The algorithm has a maximum relative error of 0.10% for the classical method simulated card area and 1.45% for the measured card area. The research of this paper provides an effective scheme for the design, prediction, and fault diagnosis of a sucker rod pumping system with an analytical solution.

1. Introduction

A sucker rod pumping system is an artificial lift instrument that is commonly installed worldwide [1]. It comprises a surface unit, a sucker rod string, and a subsurface pump [2,3]. The subsurface pump consists of a standing valve at the bottom of the well and a traveling valve attached to a rod [4]. A sucker rod pumping system is usually set up in an open-air environment and requires a long operation time, so monitoring its working conditions is very important for oil production [5]. However, its working condition is difficult to monitor directly because it operates in a small-diameter oil tube thousands of meters underground. In production practice, the pumping condition is usually identified by analyzing the surface dynamometer card [6], which is called the diagnosis technology. Many advanced analytical methods have been applied in diagnosis technology based on fault dynamometer cards [7,8,9,10,11]. It is impossible to have all kinds of fault dynamometer cards in a real oil well. One of the most economical and effective methods for obtaining the card is simulation based on the wave equation of the sucker rod string with the establishment of a pumping condition model, which is called prediction technology [12,13].
As early as 1963, a one-dimensional wave equation of a sucker rod string under normal and gas interference pumping conditions was built by Gibbs [14]. Then, for different pumping conditions, a two-dimensional wave equation [15,16], a three-dimensional wave equation [17,18] and other wave equations [19,20,21,22,23] were established. In 2020, Xiaoxiao et al. [24] simulated the working process of a sucker rod pumping system under fault conditions based on the three-dimensional wave equation of a sucker rod string by adopting operating characteristic models of the pump valves.
Usually, pumping condition models are described by the mix function of the pump displacement (or velocity) and pump load. Due to unknown pump positions at the beginning of pumping, the finite difference method is employed to solve the wave equation [14,15,16,17,18,19,20,24]. It should first be determined where the pump is in its stroke, and then the pump load should be determined within one time step [14,25]. However, the finite difference method should satisfy the stability conditions [26]. The analytical solution method based on the Fourier series for the one-dimensional wave equation of a sucker rod string can overcome this shortcoming, but it requires the pump load–time function within one pumping cycle as the boundary condition [26,27]. It is difficult to use traditional pumping condition models.
To overcome this shortcoming, a normal pumping condition model with the polished rod velocity as a function is established. An iterative prediction algorithm is proposed. The novelty of this work is that the algorithm can use the analytical solution of the wave equation to predict the behavior of pumping units and is based only on the polished rod velocity.
The remainder of the paper is organized as follows: Section 2 describes the mathematical problem of this work. Section 3 describes the approach to the problem solution. Section 4 presents the results and discussion. Section 5 presents the conclusions of the paper.

2. Mathematical Problem

2.1. One-Dimensional Sucker Rod String Wave Equation

In the actual production of oilfields, multi-tapered rod string is the majority used. The wave equation of the i-th stage rod string, along with its boundary and continuity conditions, is provided below [26]:
2 u i t 2 c i 2 2 u i x 2 + v i u i t = 0       l L i < x < l u 0 , t = u a ( t ) E r K A r K u x , t x x = L = P p ( t ) u i l , t = u i + 1 l , t E r i A r i u i x , t x x = l = E r i + 1 A r i + 1 u i + 1 x , t x x = l
where c i = E r i g / ρ i , l = k = 1 i L k , Er represents Young’s modulus of the rod, g represents the gravitational constant, ρ represents the density of the rod material, and L represents the length of the rod string.
When the weight of the sucker rod in the fluid is considered separately, the pump load is described as follows [21,23]:
P p t = A p p d p t A r e p d
where Ap is the plunger area, pd is the pump discharge pressure, p is the pump pressure, and Are is the cross area of the last section of the rod string.

2.2. Traditional Normal Pumping Conditions

Under normal pumping conditions with the tubing string anchored, only the sucker rod string undergoes elastic changes during the plunger movement caused by polished rod movement. Based on the assumption of the normal pumping condition and its movement rule [28], the relation between the pump displacement and load can be expressed as follows:
u p t = 0 , P p t < W 0 , 0 t t m P p t = W 0 , P p t W 0 , 0 t t m u p t = u m , P p t > 0 , t m < t T P p t = 0 , P p t 0 , t m < t T
where W0 = Ap(pdps), tm is the downstroke start time, um is the plunger position at the downstroke start time, and T is the pumping period. Equation (3) can be described as the famous Robin boundary condition [14,29]. Obviously, the time-varying mixed boundary of displacement and load is a complex, which switches automatically according to the pump load.

2.3. Analytical Solution of the One-Dimensional Sucker Rod String Wave Equation

After, the polished rod displacement and the pump load are approximated by the truncated Fourier series as follows [26]:
u a ( t ) = ν 0 2 + n = 1 N ( ν n cos n ω t + δ n sin n ω t ) P p ( t ) = σ 0 2 + n = 1 N ( σ n cos n ω t + τ n sin n ω t )
where N is the number of Fourier series; ν0, νn, δn, σ0, σn, and τn are the Fourier coefficients; and ω = 2πn/T, n = 1, 2, 3, …
Equation (1) is solved by a matrix expression of analytical solutions in predictive analysis by our team. Further details are given in Reference [26]. The main results are as follows:
u i x , t = ξ i + η i x + n = 1 N Q i n x cos n ω t + P i n x sin n ω t D i x , t = E r i A r i η i + E r i A r i n = 1 N Q i n x cos n ω t + P i n x sin n ω t
where
Q i n ( x ) P i n ( x ) Q i n ( x ) P i n ( x ) T = M F i n x 4 × 4 C i n 4 × 1 M F i n x 4 × 4 = M 1 i n x M 2 i n x F 1 i n x F 2 i n x C i n 4 × 1 = ε i n π i n κ i n μ i n T
M 1 i n x = cosh β i n x cos α i n x sinh β i n x sin α i n x sinh β i n x sin α i n x cosh β i n x cos α i n x i n x M 2 i n x = cosh β i n x sin α i n x sinh β i n x cos α i n x sinh β i n x cos α i n x cosh β i n x sin α i n x i n x F 1 i n x = M 1 i n x x , F 2 i n x = M 2 i n x x β i n = n ω c i 2 1 + 1 + v i n ω 2 , α i n = n ω c i 2 1 + 1 + v i n ω 2
Considering the boundary conditions, the coefficients of C n 4 K × 1 can be obtained as follows:
C n 4 K × 1 = T M n 4 K × 4 K 1 U D n 4 K × 1 C n 4 K × 1 = C 1 n 4 × 1 C i n 4 × 1 C K n 4 × 1 T U D n 4 K × 1 = U n 2 × 1 0 2 × 1 0 4 × 1 0 2 × 1 D n 2 × 1 T U n 2 × 1 = v n δ n T ,   D n 2 × 1 = σ n τ n T / E r K A r K
Thus, the pump displacement and the polished rod load can be given as follows:
u p t = ν 0 2 σ 0 2 k = 1 K L r k E r k A r k n = 1 N Q K n L cos n ω t + P K n L sin n ω t P R L t = σ 0 2 + E r 1 A r 1 n = 1 N Q 1 n 0 cos n ω t + P 1 n 0 sin n ω t
where
Q K n ( L ) P K n ( L ) T = M 1 K n L M 2 K n L C K n 4 × 1 Q 1 n ( 0 ) P 1 n ( 0 ) T = F 11 n 0 F 21 n 0 C 1 n 4 × 1
The analytical solution is mainly based on the complex theory and separation of variables to solve the wave equation proposed by Gibbs in 1967 [30]. Obviously, this analytical solution is a frequency domain method that requires a periodic pump load-time function and cannot be used to solve the definite solution problem of Equations (1) and (3).

3. Approach to the Problem Solution

To solve the definite solution problem using the analytic method, an iterative algorithm, which ranges from static to dynamic, can be employed. In the algorithm, the sucker rod string is equivalent to a spring, and a periodic pump load-time function of polished rod velocity is obtained. Taking the polished rod displacement and pump load as the initial values, the pump displacement can be calculated according to the analytical solution. According to the calculated pump displacement, the new polished rod velocity can be obtained according to the modified method. Thus, the new periodic pump load-time function of polished rod velocity is obtained. After several iterations, the pump load will become a dynamic load to meet the requirements.

3.1. Model of Normal Pumping Conditions

There are four stages in a pumping cycle: the first stage is the loading portion of the upstroke, the second is the fully loaded portion of the upstroke, the third is the unloading portion of the downstroke, and the last is the unloaded portion of the downstroke [25].
The normal pumping condition model should consider the anchored state of the tubing string. With the tubing string anchored, only the sucker rod string undergoes elastic changes during the plunger movement; meanwhile, in the anchored state, the sucker rod string and the tubing string also exhibit elastic changes caused by polished rod movement [28]. Considering the elastic movement, the sucker rod string and the tubing string can be equivalent to springs. A schematic diagram during the loading portion of the upstroke is shown in Figure 1.
As shown in Figure 1a, when the pumping speed is low enough, it has a pump load (i.e., pump plunger load) as follows:
P p t = u a t / k r
where ua(t) is the polished rod displacement, kr is the flexibility of the equivalent spring of the sucker rod string (i.e., the derivative of stiffness), and kr = Lr/(ErAr). The pump load in the unloading portion of the downstroke differs only by a constant from Equation (11).
According to Figure 1b, when the pump displacement (i.e., the position of the pump plunger) derived from the elastic movement of the tubing string is up(t),
P p t = u a t / k e
u a t = u p t k e / k t
where ke = δtkt + kr, kt is the flexibility of the tubing string, kt = Lt/(EtAt), and Lt, Et, and At are the length, Young’s modulus, and cross area of the tubing string, respectively. The δt = 0 represents the state of tubing anchored while δt = 0 represents the state of tubing unanchored. Obviously, the equivalent springs of the sucker rod string and the tubing string are in series. It is a static model of the sucker rod string system. When considering the dynamic characteristics, the spring-mass-damper system [22] can be considered as an alternative, which is under study.
Considering Equation (12) and differentiating Equation (2) yields
v a t = A p k e d p t d t
where va(t) is the velocity of the polished rod. Equation (14) is suitable for both the loading portion and the unloading portion.
When the pump pressure is less than the intake pressure, which is ps, the standing valve will open. When the pump pressure is greater than the discharge pressure, the traveling valve will open. Considering these valve opening conditions and Equation (14), the pump pressure can be given by the recurrent equation as follows:
p f d t i = p t i 1 Δ t A p k e v a t i 1 p t i = p f d t i , p f d t i > p s , 0 t i < t m p s , p f d t i p s , 0 t i < t m p f d t i , p f d t i < p d , t m t i T p d , p f d t i p d , t m t i T
where △t is the time increment.
u p t i = u a t i k r A p p d p t i
The initial value is
p t 0 = p d v p t 0 = 0
The advantage of the recurrent Equation (15) is that it does not consider the constant um in the downstroke of Equation (3), and it is easy to include other fault conditions. For example, considering the leakage state of the pump valve, Equation (15) can be modified as follows:
p f d t i = p t i 1 Δ t A p k e v a t i 1 v p l t i 1
where vpl is the plunger velocity due to leakage of the pump valve. For the standing valve and traveling valve, the plunger velocity is given as follows [24]:
v p l s = ζ s ζ p A p p t i 1 p s e s v p l t = ζ t ζ p A p p d p t i 1 e t ζ p = π D p δ 3 12 μ l p
where Dp represents the diameter of the plunger, lp represents the length of the plunger, δ represents the clearance between the plunger and pump barrel, ζs and ζt are the leakage coefficients of the standing valve and traveling valve, respectively, and es and et are the leakage exponents of the standing valve and traveling valve, respectively.

3.2. Iterative Algorithm

The pumping condition model proposed in the paper is a static model. If it is used in the pump displacement and the polished rod load resolution of Equation (9), an error will be introduced. The higher the pumping speed is, the greater the error.
To solve the definite solution problem by the analytic method, the iterative algorithm should consist of the following steps:
Step 1: Calculate the pump pressure and pump load according to the pumping model based on the polished rod velocity.
p 0 t = f v p t P p 0 t = A p p d p 0 t A r e p d
Step 2: Approximate the polished rod displacement and the pump load by the truncated Fourier series according to Equation (4).
Step 3: Calculate the pump displacement up0(t) and the polished rod load according to Equation (9).
Step 4: Modify the polished rod displacement as follows:
u a 0 t = u p 0 t + k r A p p d p 0 t
Step 5: Calculate the tolerance error as follows:
Δ u a = max u a ( t ) u a 0 ( t )
Step 6: Evaluate the tolerance error.
If △ua > ε, update the polished rod displacement as follows, and return to step 1.
u a t = u a t + u a 0 t 2
Otherwise, if △uaε, stop the program, and obtain the results of the surface dynamometer card, i.e., [ua0(t), PRL(t)], and the downhole card, i.e., [up0(t), Pp0(t)].
The flowchart of the algorithm is shown in Figure 2.
The computer program used in this paper is MATLAB R2021a, and the pseudocode of the algorithm is shown in Algorithm 1.
Algorithm 1: Prediction algorithm with iteration
Input: np, K, and a series system parameter
Calculate: ci, pd, ps, W0, ks, kt, T (i.e., 60/np), ω, and vi
Set: J = 400, t = linspace(0, T, J + 1), N = J/2
Calculate: ua(t) according to kinematic equation of pumping unit’s movement
Approximated by Fourier series with trpaz function: ua(t)→ν0, νn, and δn refer to Equation (4)
Mathematics 12 01318 i001
Output: ua0(t), PRL(t), up0(t), Pp0(t)
Note: The equation for vi is given in Reference [26].

3.3. Theoretical Analysis of the Prediction Algorithm

To prove the algorithm, the normal pumping condition with the tubing anchored and a single rod is considered as an example. During the loading portion, the pump is stationary [14], that is, up(t) = 0. Next, it needs be proven that the pump displacement approaches 0 during the loading portion after iteration.

3.3.1. Theoretical Basis

There is another analytical solution of the single rod wave equation, which is based on separating variables without using the complex method and has a clearer physical significance [27]. It is used to illustrate the mechanism of iteration. The pump displacement based on the single rod wave equation can be given as follows:
u p r t = u a t + Δ 1 u a t k r P p t k r Δ 2 P p t Δ 1 u a t = m = 0 u a m t sin ω m l Δ 2 P p t = k r m = 0 P p m t sin ω m l
where
u a m t = n = 1 P s m n cos η n t + Q s m n sin η n t P p m t = n = 1 P p m n cos η n t + Q p m n sin η n t
A s m n B s m n = 2 l ω m η n 2 ν η n ν η n η n 2 ν n δ n P s m n Q s m n = ω m c 2 η n 2 v η n v η n ω m c 2 η n 2 1 A s m n B s m n
A p m n B p m n = 2 ( 1 ) m l ω m 2 η n 2 ν η n ν η n η n 2 σ n τ n ω m c 2 η n 2 v η n v η n ω m c 2 η n 2 ] P p m n Q p m n = A p m n B p m n
where l is the length of the rod string; ωm is the eigenfrequency; and ωm = (2m + 1)π/(2l), m = 0, 1, 2, …
Obviously, in the static state, △1ua(t) = △2Pp(t) = 0, and Equation (20) is reduced to Equation (16).

3.3.2. Results of the Iterative Process

In step 1 of the iterative process,
P p 0 t = u a t k r
Thus, the pump displacement can be calculated according to step 3, i.e., Equation (20), as follows:
u p 0 t = Δ 1 u a t Δ 2 u a t
The pump displacement is not zero when it is not in the static state. The aim of the iteration is to ensure that the pump displacement is 0 during the loading portion, i.e., up(t) = 0.
According to step 4, the modified polished rod displacement is as follows:
u a 0 t = u a t + Δ 1 u a t Δ 2 u a t
Then, the new pump load according to step 1 is as follows:
P p 1 t = u a 0 t k r
The new pump displacement is calculated according to step 3, i.e., Equation (20), as follows:
u p 1 t = Δ 2 u p 0 t
According to step 4, the modified polished rod displacement is as follows:
u a 1 t = u a t + Δ 1 u a t Δ 2 u p 0 t Δ 2 u a t
Then, the new pump load is obtained according to step 1.
P p 2 t = u a 1 t k r
It has the new calculated pump displacement according to step 3, i.e., Equation (20), as follows:
u p 2 t = Δ 2 2 u p 0 t
where △2 is an operator. The △22 = △22, i.e., △22up0(t) = △2[△2up0(t)].
Therefore, by generalizing this derivation, the following expression can be obtained after the i-th iteration.
u p i t = 1 i Δ 2 i u p 0 t
If the operator △2 can reduce the value of up0(t), the value of upi(t) will approach zero after the i-th iteration, i.e., the iteration is convergence. Next, is the discussion of the convergence of iterations.

3.3.3. Convergence of Iterations

Considering the tubing anchored state, i.e., Equations (12) and (16), and according to Equations (21)–(23), the iterative matrix of the operator △2 can be written as the following expression:
M n = k r k e σ 0 2 E r K A r K m = 0 2 ( 1 ) m l ω m 2 sin ω m l A m n B m n B m n A m n
where
A m n = ω m c 2 η n 2 η n 2 v η n 2 ω m c 2 η n 2 2 + v η n 2 B m n = ω m c 2 η n 2 v η n + v η n η n 2 ω m c 2 η n 2 2 + v η n 2
The role of σ0/(2ErKArK) is to adjust the order of magnitude.
Obviously, when the maximum value of Mn is less than 1, the iteration converges. When only m = 0, the value of ωm−2 is the largest. Equations (25) and (26) show the convergence mechanism, which is that the maximum value of Mn occurs when m = 0 and η n = ω 0 c 2 v 2 2 , i.e., harmonic resonance. The disadvantage is that Equations (25) and (26) are only applicable to the case of a single rod. Further study of the case of multi-tapered rods is ongoing.
To visualize the relation between the iterative matrix elements Mn and the Fourier series number n, especially for multi-tapered rods, the iterative matrix in the matrix expression of the analytical solution [26] is given as follows:
M n = k r k e σ 0 2 E r K A r K M n L 2 × 4 M n 0 2 × 4 F n L 2 × 4 1   for   single   rod ;
M n = k r k e σ 0 2 E r K A r K M n L 2 × 4 T M i n v 1 4 × 4   for   multi-tapered   rod
where T M i n v 1 4 × 4 = T M i n v 4 K 3 : 4 K , 4 K 3 : 4 K , T M i n v = T M n 4 K × 4 K 1 .
Thus, the convergence curve between the iterative matrix elements Mn and Fourier series number n can be drawn when solving the wave equation with Equations (4)–(10).

4. Results and Discussion

4.1. Validation Study

Two methods are used to validate the prediction algorithm. One is to validate the prediction algorithm by comparison with the simulated results of the classical finite differential solution [14,25,26], and the other is based on the measured surface dynamometer card. The reason for the comparison with the finite difference method is that it is a well-known and proven solution for the predictive analysis of sucker rod pumping systems. The further details of the finite differential solution are given in reference [26] and are not repeated here.
The measured data are from two different kinds of wells belonging to an oilfield of the Sinopec Oilfield Company of China. The basic parameters of the wells are listed in Table 1.
In the analytical solution, the time increment is T/400, the Fourier series number is 200, and the tolerance error is 0.1%. In the finite differential solution, the time increment is 0.95Lmin/c, which can satisfy the stability conditions, and the tolerance error is 0.1%.

4.1.1. Comparison with the Finite Differential Solution

Well 1 is used here as an example for comparative study. After two iterations, the finite difference solution satisfies the accuracy requirement, and the simulation results are shown in Figure 3. After the first iteration, the results simulated by the analytical solution are also shown in Figure 3. The first iteration of the analytical solution means that the pumping condition model is a static model. Figure 3 shows that the surface dynamometer cards and downhole cards simulated by both solutions are consistent. The relative area error of the surface card and downhole card is 0.02%. This demonstrates the feasibility of the proposed prediction algorithm using only the static pumping condition model at a speed of 1.4 min−1. In this paper, the relative area error is defined as (AfAa)/Af × 100%, where Af is the card area of the finite difference solution, and Aa is that of the analytical solution.
To further demonstrate the applicability of the algorithm proposed in this paper, the cards at a pumping speed of 5.0 min−1 are simulated by both solutions. The simulation results of the analytical solution after the first iteration are shown in Figure 4. As illustrated in Figure 4, for the downhole cards, the loading and unloading portions of the analytical solution are significantly smaller than those of the finite difference solution. In the surface dynamometer card, the amplitude of the fluctuation in the upper stroke and downstroke of the analytical solution is slightly smaller than that of the finite difference solution. Hence, the use of the static pump simulation model at this pumping speed will result in a larger error.
After eight iterations, the accuracy requirement of the analytical solution is met, and the results are shown in Figure 5. As shown in Figure 5, the surface dynamometer cards and downhole cards simulated by both solutions are consistent, and the area relative errors of the surface card and downhole card are −0.10% and 0.01%, respectively, which demonstrates that the iterative algorithm proposed in this paper can eliminate the error of the static model and achieve the same accuracy as the finite difference solution.

4.1.2. Comparison with the Measured Card

The surfaces and downhole cards simulated by the analytical solution and the measured surface dynamometer cards of well 1 and well 2 are shown in Figure 6 and Figure 7, respectively. Table 2 compares the data taken from the simulated surface dynamometer card and the measured data, which are recorded by the dynamometer.
As illustrated in Figure 6, the simulated surface dynamometer card has good consistency with the measured card. As shown in Table 2, the loads and area predicted by the analytical solution closely matched the actual loads, especially for the card area, and the relative area error to the measured card area is 1.30%. The difference may come from errors due to dynamometer resolution.
According to the simulated surface dynamometer card by the analytical solution of well 2, the buoyant rod weight was adjusted by −0.91 kN. This error may be due to calculation method of the rod buoyancy weight different from the actual. As illustrated in Figure 7, there is little difference between the two surface dynamometer cards. Table 2 shows that the loads and area predicted by the analytical solution also closely match the actual loads, and the relative area error to the measured card area is 1.45%. These results indicate that the prediction algorithm is feasible. These differences may be due to dynamometer resolution or the neglect of hydrodynamic effects. Further research is ongoing.

4.2. Convergence Study

The convergence characteristics of the algorithm are studied for wells 1 and 2 at different pumping speeds. For well 1 at pumping speeds of 1.4 min−1 and 5.0 min−1, iterative matrix elements Mn with the increasing Fourier series are calculated according to Equation (28) and are shown in Figure 8 and Figure 9. The relative characteristic data are given in Table 3. The resonant frequency is calculated according to 2nπ/T.
Figure 8 and Figure 9 show that the values of the iterative matrix elements are less than one; so, the iterations converge. A comparison of Figure 8 and Figure 9 also shows that the number of harmonics increases, and the resonance position decreases, as the pumping speed increases. Table 3 shows that despite the different pumping speeds, the frequency of the resonance and the maximum value of the iterative matrix element are almost constant. The mechanism is shown in Equations (25) and (26), i.e., the convergence of the algorithm depends on the material of the rod, its length, the fluid viscosity, and the tube anchoring state. Table 3 also shows that the number of iterations at a pumping speed of 5.0 min−1 is greater than that at a pumping speed of 1.4 min−1. The number of iterations of the algorithm depends on the initial value and the structure of the algorithm. A comparison of Figure 4 and Figure 5 shows that the initial value at a pumping speed of 5.0 min−1 is far from the true solution. This may be the reason for the greater number of iterations at a speed of 5.0 min−1.
Figure 10 and Figure 11 show the convergence curves of well 2 at pumping speeds of 4.1 min−1 and 6.0 min−1, which are calculated according to Equation (27). The relative characteristic data are also shown in Table 3. The results show that there is the same convergence law between well 1 and well 2. The difference is that the algorithm of well 2 is more likely to converge because the maximum value of the iterative matrix elements is smaller. The resonance frequency calculated from ηn is 8.1267, which is very close to the resonance frequency in Table 3. The error is due to too few Fourier series.

5. Conclusions

This paper establishes a normal pump condition model based on the polished rod velocity. The iterative prediction algorithm proposed in this paper predicts the behavior of a sucker rod pumping unit based on an analytical solution. The algorithm is validated by classical finite difference method simulated cards and measured surface dynamometer cards. The convergence of the algorithm is analyzed theoretically and numerically. The following conclusions can be drawn:
(1)
In the normal pump condition model, the recursive equation for pump pressure is based on the polished rod velocity, which can easily provide the pump load–time function within one pumping cycle, naturally consider the anchoring state of the tubing, and include other fault conditions.
(2)
The algorithm can use the analytical solution of the wave equation to predict the behavior of the pumping unit only based on the polished rod velocity. Comparison with the simulated cards of the classical finite difference method shows that the maximum area relative error is 0.10%, and the proposed algorithm can achieve the same level of accuracy as the classical finite difference method. When compared with the measured surface cards, the area relative error is 1.45%, indicating that the algorithm is feasible.
(3)
The convergence of the algorithm is analyzed theoretically. An expression for the iteration matrix is given, which can be applied to both single rod and multi-tapered rods. The expression shows that the convergence of the algorithm depends on the material of the rod, its length, the fluid viscosity, and the tube anchoring state. Numerical results show that the algorithm converges in the two wells given in this paper. The smaller the maximum value of the iterative matrix elements, the easier it is for the algorithm to converge. The convergence analysis provides assurance of the accuracy and reliability of the algorithm.
Based on the limitations of this study, the authors would like to propose the following directions for future research: (a) a more widely applicable iterative algorithm; (b) a theoretical study of the convergence of multi-tapered sucker rod pumping systems; and (c) an iterative algorithm study of multi-fault pumping conditions.

Author Contributions

Conceptualization, J.Y.; methodology, J.Y.; software, J.Y.; validation, J.Y.; formal analysis, J.Y.; investigation, J.Y.; resources, J.Y.; data curation, J.Y.; writing—original draft preparation, J.Y.; writing—review and editing, J.Y.; visualization, J.Y.; supervision, J.Y.; project administration, H.M.; funding acquisition, H.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Shandong province Natural Science Foundation of China (Grant No. ZR2021MD067), the Fundamental Research Funds for the Central Universities (Grant No. 22CX03011A).

Data Availability Statement

The data presented in this study are available on request.

Acknowledgments

Comments and suggestions from the editor and reviewers are very much appreciated.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Takacs, G. A critical analysis of power conditions in sucker-rod pumping systems. J. Petrol. Sci. Eng. 2022, 210, 110061. [Google Scholar] [CrossRef]
  2. Ramos, A.J.A.; Almeida Júnior, D.S.; Santos, A.R.; Araújo, E.A.; Aum, P.T. Energy dissipation analysis of oil wells sucker rod string. J. Braz. Soc. Mech. Sci. Eng. 2020, 42, 108. [Google Scholar] [CrossRef]
  3. Langbauer, C.; Langbauer, T.; Fruhwirth, R.; Mastobaev, B. Sucker rod pump frequency-elastic drive mode development–from the numerical model to the field test. Liq. Gas. Energy Resour. 2021, 1, 64–85. [Google Scholar] [CrossRef]
  4. Bangert, P. Diagnosing and Predicting Problems with Rod Pumps using Machine Learning. In Proceedings of the SPE Middle East Oil and Gas Show and Conference, Manama, Bahrain, 18–21 March 2019. [Google Scholar] [CrossRef]
  5. Wang, X.; He, Y.; Li, F.; Dou, X.; Wang, Z.; Xu, H.; Fu, L. A Working Condition Diagnosis Model of Sucker Rod Pumping Wells Based on Big Data Deep Learning. In Proceedings of the International Petroleum Technology Conference, Beijing, China, 26–28 March 2019. [Google Scholar] [CrossRef]
  6. Xu, J. A Method for Diagnosing the Performance of Sucker Rod String in Straight Inclined Wells. In Proceedings of the SPE Latin America/Caribbean Petroleum Engineering Conference, Buenos Aires, Argentina, 27–29 April 1994. [Google Scholar] [CrossRef]
  7. Zheng, B.; Gao, X. Sucker rod pumping diagnosis using valve working displacement and parameter optimal continuous hidden Markov model. J. Process Cont. 2017, 59, 1–12. [Google Scholar] [CrossRef]
  8. Zheng, B.; Guo, X.; Ki, X. Diagnosis of Sucker Rod Pump based on generating dynamometer cards. J. Process Cont. 2019, 77, 76–88. [Google Scholar] [CrossRef]
  9. Zhang, A.; Gao, X. Supervised dictionary-based transfer subspace learning and applications for fault diagnosis of sucker rod pumping systems. Neurocomputing 2019, 338, 293–306. [Google Scholar] [CrossRef]
  10. He, Y.; Zang, C.; Zeng, P.; Wang, M.; Wan, G.; Dong, Q. Automatic Recognition of Sucker-Rod Pumping System Working Conditions Using Few-Shot Indicator Diagram Based on Meta-learning. In Proceedings of the International Conference on Intelligent Automation and Soft Computing, Chicago, IL, USA, 26–28 May 2021. [Google Scholar] [CrossRef]
  11. Han, Y.; Li, K.; Ge, F.; Xu, W. Online fault diagnosis for sucker rod pumping well by optimized density peak clustering. ISA Trans. 2022, 120, 222–234. [Google Scholar] [CrossRef] [PubMed]
  12. Takacs, G. Calculation of Operational Parameters. In Sucker-Rod Pumping Handbook: Production Engineering Fundamentals and Long-Stroke Rod Pumping; Gulf Professional Publishing: Houston, TX, USA, 2015; Chapter 4; pp. 247–376. [Google Scholar] [CrossRef]
  13. Moreno, G.A.; Garriz, A.E. Sucker rod string dynamics in deviated wells. J. Petrol. Sci. Eng. 2019, 184, 106534. [Google Scholar] [CrossRef]
  14. Gibbs, S.G. Predicting the behavior of sucker-rod pumping systems. SPE J. 1963, 15, 769–778. [Google Scholar] [CrossRef]
  15. Doty, D.R.; Schmidt, Z. An Improved Model for Sucker Rod Pumping. SPE J. 1983, 23, 1. [Google Scholar] [CrossRef]
  16. Lekia, S.D.L.; Evans, R.D. A Coupled Rod and Fluid Dynamic Model for Predicting the Behavior of Sucker-Rod Pumping Systems. SPE Prod. Fac. 1995, 10, 26–33. [Google Scholar] [CrossRef]
  17. Yu, G.A.; Wu, Y.J.; Wang, G.Y. Three dimensional vibration in a sucker rod beam pumping system. Acta Petrolei Sin. 1989, 10, 76–83. (In Chinese) [Google Scholar] [CrossRef]
  18. Lollback, P.A.; Wang, G.Y.; Rahman, S.S. An alternative approach to the analysis of sucker-rod dynamics in vertical and deviated wells. J. Pet. Sci. Eng. 1997, 17, 313–320. [Google Scholar] [CrossRef]
  19. Khodabandeh, A.; Miska, S. A Simple Method for Predicting the Performance of a Sucker-Rod Pumping System. In Proceedings of the SPE Eastern Regional Meeting, Lexington, KY, USA, 22–25 October 1991. [Google Scholar] [CrossRef]
  20. Khodabandeh, A.; Miska, S. A New Approach for Modeling Fluid Inertia Effects on Sucker-Rod Pump Performance and Design. In Proceedings of the SPE Rocky Mountain Regional Meeting, Casper, WY, USA, 18–21 May 1992. [Google Scholar] [CrossRef]
  21. Xing, M.; Dong, S. An improved longitudinal vibration model and dynamic characteristic of sucker rod string. J. Vibroeng. 2014, 16, 3432–3448. [Google Scholar]
  22. Xing, M. Response analysis of longitudinal vibration of sucker rod string considering rod buckling. Adv. Eng. Softw. 2016, 99, 49–58. [Google Scholar] [CrossRef]
  23. Li, W.; Dong, S.; Sun, X. An Improved Sucker Rod Pumping System Model and Swabbing Parameters Optimized Design. Math. Probl. Eng. 2018, 2018, 4746210. [Google Scholar] [CrossRef]
  24. Lv, X.; Wang, H.; Liu, Y.; Chen, S.; Lan, W.; Sun, B. A novel method of output metering with dynamometer card for SRPS under fault condition. J. Pet. Sci. Eng. 2020, 192, 107098. [Google Scholar] [CrossRef]
  25. Schafer, D.J.; Jennings, J.W. An investigation of analytical and numerical sucker rod pumping mathematical models. In Proceedings of the SPE Annual Technical Conference and Exhibition, Dallas, TX, USA, 24 September 1987. [Google Scholar] [CrossRef]
  26. Yin, J.; Sun, D.; Yang, Y. Predicting multi-tapered sucker-rod pumping systems with the analytical solution. J. Pet. Sci. Eng. 2021, 197, 108115. [Google Scholar] [CrossRef]
  27. Yin, J.; Sun, D.; Yang, Y. A Novel Method for Diagnosis of Sucker-Rod Pumping Systems Based on the Polished-Rod Load Vibration in Vertical Wells. SPE J. 2020, 25, 2470–2481. [Google Scholar] [CrossRef]
  28. Takacs, G. The Analysis of Sucker-Rod Pumping Installations. In Sucker-Rod Pumping Handbook: Production Engineering Fundamentals and Long-Stroke Rod Pumping; Gulf Professional Publishing: Houston, TX, USA, 2015; Chapter 6; pp. 423–503. [Google Scholar] [CrossRef]
  29. DaCunha, J.J.; Gibbs, S.G. Modeling a Finite-Length Sucker Rod Using the Semi-Infinite-Wave Equation and a Proof to Gibbs’ Conjecture. Soc. Petrol. Eng. 2009, 14, 112–119. [Google Scholar] [CrossRef]
  30. Gibb, S.G. Method of Determining Sucker Rod Pump Performance. U.S. Patent 3343409, 26 September 1967. [Google Scholar]
Figure 1. Schematic diagram of the normal pumping condition model.
Figure 1. Schematic diagram of the normal pumping condition model.
Mathematics 12 01318 g001
Figure 2. Flowchart of the prediction algorithm.
Figure 2. Flowchart of the prediction algorithm.
Mathematics 12 01318 g002
Figure 3. Simulated cards with a pumping speed of 1.4 min−1.
Figure 3. Simulated cards with a pumping speed of 1.4 min−1.
Mathematics 12 01318 g003
Figure 4. Simulated cards with a pumping speed of 5.0 min−1 and the first iteration of the analytic solution.
Figure 4. Simulated cards with a pumping speed of 5.0 min−1 and the first iteration of the analytic solution.
Mathematics 12 01318 g004
Figure 5. Simulated cards with the pumping speed of 5.0 min−1.
Figure 5. Simulated cards with the pumping speed of 5.0 min−1.
Mathematics 12 01318 g005
Figure 6. Measured surface dynamometer card and simulated cards of well 1.
Figure 6. Measured surface dynamometer card and simulated cards of well 1.
Mathematics 12 01318 g006
Figure 7. Measured surface dynamometer card and simulated cards of well 2.
Figure 7. Measured surface dynamometer card and simulated cards of well 2.
Mathematics 12 01318 g007
Figure 8. Convergence curve of well 1 at a pumping speed of 1.4 min−1.
Figure 8. Convergence curve of well 1 at a pumping speed of 1.4 min−1.
Mathematics 12 01318 g008
Figure 9. Convergence curve of well 1 at a pumping speed of 5.0 min−1.
Figure 9. Convergence curve of well 1 at a pumping speed of 5.0 min−1.
Mathematics 12 01318 g009
Figure 10. Convergence curve of well 2 at pumping speed of 4.1 min−1.
Figure 10. Convergence curve of well 2 at pumping speed of 4.1 min−1.
Mathematics 12 01318 g010
Figure 11. Convergence curve of well 2 at pumping speed of 6.0 min−1.
Figure 11. Convergence curve of well 2 at pumping speed of 6.0 min−1.
Mathematics 12 01318 g011
Table 1. Basic parameters of oil well.
Table 1. Basic parameters of oil well.
ItemsValuesValues
Pumping unitLong-stroke pumping unit
(Rotaflex)
Beam pumping unit
(CYJ14-4.8-73HB)
Pump stroke, m6.04.2
Pumping speed, min−11.44.1
Sucker rod string, mm×m25 × 372.5 + 22 × 518.222 × 986.2
Tubing string, mm×m76 × 856.7, unanchored62 × 980.6, unanchored
Pump diameter, mm6357
Pump depth, m900.6997.3
Fluid density, kg/m3998.46990.30
Dynamic liquid level, m629579
Oil pressure, MPa0.90.1
Casing pressure, MPa00
Fluid viscosity, mPa.s800747.5
Gas/oil ratio, m3/m300
Rod and tube’s density, kg/m37850 (Steel)7850 (Steel)
Rod and tube’s Young’s modulus, GPa210 (Steel)210 (Steel)
Table 2. Comparison of the data predicted from the algorithm with the dynamometer recorded data.
Table 2. Comparison of the data predicted from the algorithm with the dynamometer recorded data.
Well NumberMin Load (kN)Max Load (kN)Area (kNm)
Well 1simulated20.9051.90140.35
measured23.7750.81142.19
Well 2simulated17.4648.2892.11
measured18.3847.2590.80
Table 3. Comparison of the data convergence results.
Table 3. Comparison of the data convergence results.
Well
Number
Pumping
Speed (min−1)
Iteration
Number
Resonance
Position
Resonance
Frequency (Hz)
Max |MnA| (m)Max |MnB| (m)
Well 11.43679.82270.30640.5386
5.08199.94840.29330.5159
Well 24.13198.15770.12280.2009
6.03138.16810.11740.2010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yin, J.; Ma, H. A Novel Method for Predicting the Behavior of a Sucker Rod Pumping Unit Based on the Polished Rod Velocity. Mathematics 2024, 12, 1318. https://doi.org/10.3390/math12091318

AMA Style

Yin J, Ma H. A Novel Method for Predicting the Behavior of a Sucker Rod Pumping Unit Based on the Polished Rod Velocity. Mathematics. 2024; 12(9):1318. https://doi.org/10.3390/math12091318

Chicago/Turabian Style

Yin, Jiaojian, and Hongzhang Ma. 2024. "A Novel Method for Predicting the Behavior of a Sucker Rod Pumping Unit Based on the Polished Rod Velocity" Mathematics 12, no. 9: 1318. https://doi.org/10.3390/math12091318

APA Style

Yin, J., & Ma, H. (2024). A Novel Method for Predicting the Behavior of a Sucker Rod Pumping Unit Based on the Polished Rod Velocity. Mathematics, 12(9), 1318. https://doi.org/10.3390/math12091318

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop