Research Article Mathematical Modeling and

Hindawi Publishing Corporation Advances in Mechanical Engineering Volume 2013, Article ID 125936, 7 pages http://dx.doi.org/10.1155/2013/125936

Research Article Mathematical Modeling and Validation on a New Valve Core of the Throttle Valve in MPD Y. K. Fu,1 Q. Y. Liu,2 G. R. Wang,1 and S. Y. Tao1 1 2

School of Mechatronic Engineering, Southwest Petroleum University, Chengdu 610500, China State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu 610500, China

Correspondence should be addressed to Q. Y. Liu; [email protected] Received 29 May 2013; Revised 19 August 2013; Accepted 19 August 2013 Academic Editor: M. Affan Badar Copyright © 2013 Y. K. Fu et al. 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. To achieve different demands for pressure control in managed pressure drilling (MPD), the relationship between the pressure difference and the valve opening was analyzed in the working process of MPD throttle. The throttle flow area equation and the flow area isosurface equation at different openings were derived. Then, the mathematical model of the throttle valve surface shape was established by solving the envelope curve of isosurface cluster. Furthermore, the shape curve equation of the valve core surface was derived and modified by applying this mathematical model. Finally, to verify the performance of this investigation, the comparison of the results between simulation and design values was performed. The results showed that the correlation between the throttle pressure difference and the valve opening is linearly related, which proves the correctness of the mathematical model for the throttle valve surface.

1. Introduction With the growing exhaustion of the oil and gas resources, the exploration of the formation is becoming deeper and more complex. The complex drilling problems, such as the kick, well leakage, and the leakage of harmful gas, are increasing rapidly. As a new technology of oil drilling, the managed pressure drilling (MPD) can realize precise control in the whole drilling, pipe connection, and tripping processes by controlling the surface casing pressure, and it has become an effective way to alleviate these complicated situations previously mentioned [1, 2]. As one of the key equipment to achieve precise control of the wellhead back pressure in MPD, the throttle valve can be used to control the bottom hole pressure by adjusting the opening size of throttle valve to form different throttle pressures [3, 4]. With the development of MPD, the requirements of the wellhead back pressure control have gradually been transformed form high-pressure difference to high-precision control, and the study of the variation between the throttle pressure and the valve opening size has become the core

problem [5]. As the main component of adjusting the throttle pressure in MPD, the shape curve of the valve core surface is directly related to the precision of the wellhead back pressure control [6]. At present, the valve core that is used in the conventional drilling consists mainly of needle valve core, cylinder valve core, and wedge valve core, for which numerous studies have been carried out and has already formed a perfect theory system. On account of the throttle valve that is used in MPD, some manufacturers suggest the use of throttle valve in conventional drilling and save the cost to a certain extent, but the throttle pressure with different openings present a nonlinear relationship that is prone to failure and hence cannot meet the requirements of MPD [7, 8]. Therefore, a special throttle valve used in MPD has been developed by T3 Energy Services and Expro International Group Ltd. [9]. The valve products of these two companies have been widely used by Schlumberger, Weatherford, Optimal, and other oil companies. In these throttle valves, the flow coefficient (Cv) changes with a linear relationship in the range of 30%–70% valve openings, but the throttle pressure changes nonlinearly which cannot meet all

Downloaded from ade.sagepub.com by guest on August 26, 2015

2

Advances in Mechanical Engineering

the conditions in MPD [10, 11]. Therefore, it is necessary to develop a new throttle valve, in which the throttle pressure changes linearly with the valve opening. This paper studies the relationship between the throttle pressure difference and the valve opening in the working process of throttle valve in MPD and derives the throttle flow area equation and the flow area isosurface equation. Then, the mathematical model of the throttle valve surface shape is established by solving the envelope curve of isosurface cluster. Furthermore, the shape curve equation of the valve core surface is derived and modified by applying the mathematical model. Finally, to verify the performance of this investigation, a comparison of simulation and design values is performed.

6 5 7 4 8 9 3 10 2 11 1

2. Working Principle of the Throttle Valve As a key equipment of controlling the wellhead back pressure accurately in MPD, the regulating performance of throttle valve plays an important role in the control precision of the wellhead back pressure. Figure 1 shows the sectional view of the throttle valve. The throttle valve is mainly composed of a valve body, valve cover, screw drive assembly, valve seat, and a valve core assembly. The screw drive assembly is the central drive of the throttle valve to perform an action, and the valve core assembly is the key component to regulate the wellhead back pressure precisely. In the working process, the hydraulic motor or a hand wheel rotates the screw drive assembly which gives the spool assembly a vertical reciprocating motion, thus changing the valve opening. This provides adjustment in the throttle pressure and ultimately controls the wellhead back pressure.

Figure 1: The sectional view of the throttle valve. (1) valve body, (2) valve bonnet, (3) union nut, (4) connection flange, (5) screw drive assembly, (6) handwheel, (7) Coupling, (8) packing gland, (9) combined seal, (10) valve core assembly, and (11) valve seat.

From the boundary conditions 𝑙 = 𝑙max 󳨀→ Δ𝑝 = Δ𝑝min , 𝑙 = 𝑙min 󳨀→ Δ𝑝 = Δ𝑝max .

Substituting the boundary conditions into (1), the values of the 𝐾 and 𝐶 can be figured out as 𝐾 = 1,

3. Development of the Mathematical Model

Substituting (3) into (1), we can rewrite (1) as follows: Δ𝑝 𝑙 =1− . Δ𝑝max 𝑙max

𝑄 = 𝐾V √ 𝐾V = 5.09

(1)

where Δ𝑝 is the chocking pressure at one valve opening, Δ𝑝max is the throttle pressure at the minimum opening, Δ𝑝/Δ𝑝max is the relative pressure, 𝑙 is the displacement of valve core at one opening, 𝑙max is the displacement of valve core at maximum opening, and 𝑙/𝑙max is the relative displacement of the valve core.

(4)

3.2. The Flow Area Equation. According to the flow equation of throttle valve

3.1. Pressure Characteristic Equation. According to the throttle pressure characteristics of throttle valve in MPD, the throttle pressure changes linearly with valve opening and is given by the following equation: Δ𝑝 𝑙 =𝐾 + 𝐶, Δ𝑝max 𝑙max

(3)

𝐶 = −1.

Being the main component for adjusting throttle pressure in MPD, the valve core design is directly related to precisely control the wellhead back pressure. Therefore, the study and development of the mathematical model for the shape curve of the throttle valve core surface becomes the key point [12– 16].

(2)

Δ𝑝 , 𝜌

𝜋𝐷2 4√𝐴2 /𝐴2𝑟 − 1

(5) ,

(6)

where 𝑄 is the volume flow rate, 𝐾V is the coefficient of flow, 𝐴 is the flow area at the maximum opening, and 𝜌 is the destiny of fluid. Solving (5)-(6) simultaneously, the flow area equation of throttle valve can be written as 𝐴𝑟 = √

𝜌𝜋2 𝑄2 𝐷4 . 5.092 𝜋2 𝐷4 × (1 − 𝑙/𝑙max ) Δ𝑝max + 16𝑄2 𝜌


(7)

3

Ar

1

y


N1

N a

M1

M

d

X

D

Figure 2: Schematic diagram of the flow area.

3.3. Mathematical Model of the Valve Core Surface. As shown in Figure 2, 𝐿 𝑀𝑁 is the minimum distance between the valve core and the valve seat, and the flow area 𝐴 𝑟 is the lateral area of pyramidal face 𝑀𝑁𝑁1 𝑀1 at one opening. From the lateral area equation of pyramidal face, the flow area 𝐴 𝑟 can be calculated as 𝜋 𝐴 𝑟 = (𝐷 + 𝑑) 𝐿 𝑀𝑁, (8) 2 where 𝑑 = 𝐷 − 2𝑥, 𝐿 𝑀𝑁

where 𝐷 is the maximum diameter, 𝛼 is the included angle, and 𝑥 and 𝑦 are, respectively, the transverse and longitudinal coordinates of 𝑁 at the valve core curve. Using 𝑀 as the reference point at one opening, the level set curve is formed by the circular truncated cone at different 𝛼 values, as shown in Figure 2. With the increase of the opening, the equation for the level set curve cluster can be written as 𝐹 = 𝑥2 − 𝐷𝑥

(9)

𝑥 = . cos 𝛼

(10)

Solving (9)-(10) simultaneously, the flow area equation becomes 𝐴 (11) 𝑥2 − 𝐷𝑥 + 𝑟 cos 𝛼 = 0. 𝜋 Substituting (7) into (11), the flow area equation takes the form

+

𝜌𝜋2 𝑄2 𝐷4 cos 𝛼 √ . 2 2 4 𝜋 5.09 𝜋 𝐷 × (1 − 𝑙/𝑙max ) Δ𝑝max + 16𝑄2 𝜌 (13)

As an envelope curve cluster of 𝐹, the curve of the valve core must meet the following conditions:

2

𝑥 − 𝐷𝑥 +

2

2

4

𝜌𝜋 𝑄 𝐷 cos 𝛼 √ = 0, 2 2 4 𝜋 5.09 𝜋 𝐷 × (1 − 𝑙/𝑙max ) Δ𝑝max + 16𝑄2 𝜌 (12)

𝐹 = 0,

(14a)

𝜕𝐹 = 0. 𝜕𝛼

(14b)

According to (14a), the level set curve cluster at different openings can be written as


4


𝑥2 − 𝐷𝑥 +

Table 1: The data of calculation example.

𝜌𝜋2 𝑄2 𝐷4 cos 𝛼 √ = 0. 𝜋 5.092 𝜋2 𝐷4 × (1 − 𝑙/𝑙max ) Δ𝑝max + 16𝑄2 𝜌 (15)

According to (14b), 𝜕𝐹 cos 𝛼 = 𝜕 (𝑥2 − 𝐷𝑥 + 𝜕𝛼 𝜋 ×√ × (𝜕𝛼)−1 . (16) According to Figure 2, the displacement of valve core at one opening 𝑙 can be written as (17)

Equations (9), (15), (16), and (17) represent the mathematical model of the valve core which can be solved simultaneously 𝜕 (𝑥2 − 𝐷𝑥 +

×√

cos 𝛼 𝜋

𝜌𝜋2 𝑄2 𝐷4 ) 5.092 𝜋2 𝐷4 × (1 − 𝑙/𝑙max ) Δ𝑝max + 16𝑄2 𝜌

× (𝜕𝛼)−1 = 0, 𝑥2 − 𝐷𝑥 + ×√

Values 1700 kg/m3 20 L/s 38.1 mm 28 mm 10 MPa 1 mm

Table 2: The coordinates of the valve core at different openings.

𝜌𝜋2 𝑄2 𝐷4 ) 5.092 𝜋2 𝐷4 × (1 − 𝑙/𝑙max ) Δ𝑝max + 16𝑄2 𝜌

𝑙 = 𝑦 + 𝑥 tan 𝛼.

Parameters Density 𝜌 Flow rate 𝑄 Diameter 𝐷 The maximum opening 𝑙max The maximum throttle pressure Δ𝑝max Opening interval Δ𝑙

cos 𝛼 𝜋

𝜌𝜋2 𝑄2 𝐷4 = 0, 5.092 𝜋2 𝐷4 × (1 − 𝑙/𝑙max ) Δ𝑝max + 16𝑄2 𝜌 𝑙 = 𝑦 + 𝑥 tan 𝛼,

Openings (mm) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27

Coordinate 𝑥 1.6 1.6 1.61 1.62 1.66 1.69 1.72 1.78 1.84 1.9 1.98 2.06 2.15 2.26 2.38 2.53 2.7 2.92 3.2 3.52 3.74 4.33 6 6.12 6.3 6.54 6.71

𝑦 −0.5044 −1.206 −1.963 −2.72 −3.489 −4.189 −4.933 −5.706 −7.199 −8.657 −10.18 −11.66 −13.12 −14.62 −16.08 −17.58 −19.02 −20.44 −21.95 −23.15 −23.75 −24.57 −26.97 −27.16 −27.48 −28 −28.5

𝑑 = 𝐷 − 2𝑥. (18)

4. Solution and Curve Fitting

The curve fit equation of the valve core surface is 𝑦 = −0.03736𝑥5 + 0.92553𝑥4 − 9.23562𝑥3

4.1. The Valve Core Curve Equation. To verify the validity of the mathematical model of throttle valve core, a project from Sinopec of China is used as a calculation example; the parameters of this example are shown in Table 1. A computer program was written for the mathematical model of throttle valve core and its coordinates were calculated at different openings. The results are shown in Table 2. Fitting the coordinate values in Table 2, the shape curve of the valve core surface is shown in Figure 3.

+ 46.52065𝑥2 − 119.84944𝑥 + 103.3261.

(19)

4.2. Modification of the Fitting Curve Equation. By calculating the distance between the valve seat surface and valve core curve, the minimum flow area was calculated using the equation 𝐴 𝑟 = 𝜋𝐿 𝑀𝑁(𝐷−𝑥) and the correctness of the fitting curve equation was verified in comparison with the design values.



5

Radial coordinate (mm) 0

0

1

2

3

4

5

6

7

0

1

2

Radial coordinate (mm) 3 4 5

6

7

−5 Axial coordinate (mm)

−5 Axial coordinate (mm)

0

−10 −15 −20

−10 −15 −20

−25

−25

−30

−30

Equation

y = intercept + B1∗x1 + B2∗x2 + B3∗x3 + B4∗x4 + B5∗x5

Adj. R2

0.99892

B

Intercept B1 B2 B3 B4 B5

Equation

Value Standard error 103.3261 10.08794 −119.84944 15.77034 9.23296 46.52065 −9.23562 2.53198 0.92553 0.32642 0.01595 −0.03736

y = intercept + B1∗x1 + B2∗x2 + B3∗x3

+ B4∗x4 + B5∗x5 + B6∗x6

Adj. R2

1

B

Intercept B1 B2 B3 B4 B5 B6

Figure 3: The curve fit for the valve core using values from Table 2.

Value Standard error 139.71239 0.59961 −184.9369 1.21714 92.70609 0.98986 −25.89161 0.41264 4.1516 0.09298 −0.3564 0.01074 4.96246E−4 0.01264

Figure 4: The modification of the valve core fitting curve.

Fitting the coordinate values in Table 3, the modified shape curve of the valve core surface is shown in Figure 4. From Figure 4, the modified curve fit equation of the valve core can be written as (20). Using this equation, a new throttle valve core can be designed and manufactured: 5

4

𝑦 = 0.01264𝑥 − 0.3564𝑥 + 4.1516𝑥 − 25.89161𝑥 + 92.70609𝑥2 − 184.9369𝑥 + 139.71239.

3

(20)

4.3. Relationship of Throttle Pressure with Different Openings. One has 𝑄 = 5.09 ×

𝐴 √𝐴2 /𝐴2𝑟 − 1

×√

Δ𝑃 . 𝜌

6

4

2

(21) 0

From the flow equation of throttle valve, the throttle pressure is derived using (21) resulting in (22) 𝑄2 𝜌 𝐴2 − 1) . Δ𝑃 = ( 5.092 𝐴2 𝐴2𝑟

8 Throttle pressure (MPa)

6

10

0

3

6

9

12 15 18 Opening (mm)

21

24

27

30

Figure 5: The throttle pressure as a function of valve opening.

(22)

The throttle pressure can be calculated by (22) after the flow area at different openings is calculated by (7), and this relationship of throttle pressure with valve opening is shown in Figure 5. As shown in Figure 5, with the increase of the valve opening, the throttle pressure decreases linearly, thus making it easier to control the wellhead back pressure.

5. Comparison of Simulation Results with Design Values The flow model at different openings is analyzed using computer simulation. The throttle pressure from the simulation results is compared with the design values, and the feasibility of the valve core is verified [17–20].


6

Advances in Mechanical Engineering Table 3: The modified coordinates.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27

Coordinate 𝑥 1.591 1.625 1.659 1.696 1.734 1.775 1.817 1.862 1.910 1.960 2.014 2.072 2.134 2.202 2.275 2.356 2.446 2.545 2.657 2.784 2.929 3.096 3.287 3.504 3.741 3.983 4.217

𝑦 −0.949 −1.961 −2.930 −3.937 −4.924 −5.937 −6.921 −7.918 −8.921 −9.901 −10.892 −11.883 −12.864 −13.854 −14.827 −15.806 −16.783 −17.739 −18.690 −19.624 −20.534 −21.411 −22.235 −22.986 −23.631 −24.141 −24.519

Valve core

Inflow

Outflow

Figure 6: The three-dimensional flow model. 10

8 Throttle pressure (MPa)

Openings (mm)

6

4

2

0 0

5.2. Comparison and Validation. To verify the feasibility of this design, the simulation results of throttle pressure are compared with the design values. Figure 8 shows the variation of the designed throttle pressure with simulation results at different openings. The simulation results of throttle pressure are a little smaller than the design values when the opening is between 0 and 18 mm, they are almost equal when the opening is between 21 and 24 mm, and they are little smaller at an opening of 27 mm. The numerical simulation accuracy is affected by number of cells in the computational grid quality and other factors, and the small deviation of throttle pressure between the simulation

6

9


21

24

27

30

Figure 7: Simulation results of throttle pressure at different openings. 12 10 Throttle pressure (MPa)

5.1. Simulation Modeling. Due to the fact that throttle pressure changes at different valve openings, the numerical simulation of three-dimensional flow model of Figure 6 is studied at ten different valve openings [21, 22], and the results are shown in Figure 7. The maximum throttle pressure is found to be 8.825 MPa at an opening of 0 mm, and the minimum throttle pressure is 0.722 MPa at an opening of 27 mm. The throttle pressure decreases as the valve is opened gradually and this decrease is approximately linear.

3

8 6 4 2 0 0

3

6

9


21

24

27

30

Simulation results Design values

Figure 8: The comparison of the design values with simulation results at different openings.



7

and design values can hardly be avoided. Therefore, the simulation results of throttle pressure with the design values are nearly the same, which verifies the correctness of the design of the valve core.

6. Conclusions

[11]

[12]

(1) The pressure characteristics of throttle valve in MPD are fully considered, and the mathematical model of the valve core surface is derived, solved, and modified. (2) The throttle pressure of the throttle valve at different openings is solved with the modified equations of valve core surface, and the variation of throttle pressure with valve opening is approximately linear. (3) Based on the simulation analysis of flow model at different openings, the variation of simulated throttle pressure with valve opening is approximately linear. (4) The simulation results of throttle pressure with the design values are nearly the same, which verifies the correctness of the design of the valve core surface.

[13]

[14]

[15]

[16]

Acknowledgment [17]

The work was supported by the Natural Science Foundation of China (Grant no. 51134004). [18]

References [1] Y. C. Zhou, M. Cui, and Y. J. Zha, “Discussion and prospect of managed pressure drilling technology,” Petroleum Drilling Techniques, vol. 36, no. 4, pp. 1–4, 2008. [2] B. Huang, X. B. Shi, and Z. L. Li, “Research and application of MPD technique,” Drilling & Production Technology, vol. 33, no. 5, pp. 1–4, 2010. [3] J. Kozicz, “Managed-pressure drilling—recent experience, potential efficiency gains, and future opportunities,” in Proceedings of the IADC/SPE Asia Pacific Drilling Technology Conference and Exhibition, pp. 391–395, November 2006. [4] A. K. Das, J. R. Smith, and P. J. Frink, “Simulations comparing different initial responses to kicks taken during managed pressure drilling,” in Proceedings of the SPE/IADC Drilling Conference and Exhibition, pp. 1231–1243, March 2008. [5] T. Zhang, G. H. Liu, and J. Li, “Research on multi-level parallel choke manifold system,” Petroleum Drilling Techniques, vol. 40, no. 2, pp. 98–103, 2012. [6] Y. Q. Jin, G. Liu, and Z. Q. Sun, “Research on the relation between throttle valve opening and throttle pressure in pressure control drilling,” China Petroleum Machinery, vol. 40, no. 10, pp. 11–14, 2012. [7] Q. Zhao, K. Wang, and Q. F. Guo, “A new choke manifold for the throttle control pressure drilling,” China, 200920246846.8[P].2010-08-11. [8] S. L. Fang, X. W. Yang, and K. Wang, “A controlled pressure drilling throttling device with two-stage throttle function,” China, 201020270446.3[P]. 2011-03-30. [9] M. E. Parsons, “Erosion resistant wear sleeve for high pressure valve: USA,” US6446664B1 [P]. 2002-09-10. [10] A. Grace and P. Frawley, “Experimental parametric equation for the prediction of valve coefficient (Cv) for choke valve trims,”

[19]

[20]

[21]

[22]

International Journal of Pressure Vessels and Piping, vol. 88, no. 2-3, pp. 109–118, 2011. H. Chattopadhyay, A. Kundu, B. K. Saha, and T. Gangopadhyay, “Analysis of flow structure inside a spool type pressure regulating valve,” Energy Conversion and Management, vol. 53, no. 1, pp. 196–204, 2012. B. Kan, M. W. Jin, and J. N. Ding, “Mathematical modeling of the interference of seal pair in triple-offset butterfly valve,” Proceedings of the Institution of Mechanical Engineers C, vol. 226, no. 12, pp. 3026–3031, 2012. S. H. Somashekhar, M. Singaperumal, and R. Krishna Kumar, “Mathematical modelling and simulation of a jet pipe electrohydraulic flow control servo valve,” Proceedings of the Institution of Mechanical Engineers I, vol. 221, no. 3, pp. 365–382, 2007. J. Benajes, E. Reyes, and J. M. Lujân, “Modelling study of the scavenging process in a turbocharged diesel engine with modified valve operation,” Proceedings of the Institution of Mechanical Engineers C, vol. 210, no. 4, pp. 383–393, 1996. M. Y. Salloom and Z. Samad, “Design and modeling magnetorheological directional control valve,” Journal of Intelligent Material Systems and Structures, vol. 23, no. 2, pp. 155–167, 2012. C. A. Kitio Kwuimy and C. Nataraj, “Modeling and dynamic analysis of a magnetically actuated butterfly valve,” Nonlinear Dynamics, vol. 70, no. 1, pp. 435–451, 2012. N. D. Vaughan, D. N. Johnston, and K. A. Edge, “Numerical simulation of fluid flow in poppet valves,” Proceedings of the Institution of Mechanical Engineers C, vol. 206, no. 2, pp. 119– 127, 1992. S.-P. Rwei, H.-I. Pai, and I.-C. Wang, “Fluid simulation of the airflow in interlacing nozzles,” Textile Research Journal, vol. 71, no. 7, pp. 630–634, 2001. H.-X. Wan, J. Fang, and H. Huang, “Numerical simulation on a throttle governing system with hydraulic butterfly valves in a marine environment,” Journal of Marine Science and Application, vol. 9, no. 4, pp. 403–409, 2010. J. Venkatesan, G. Nagarajan, R. V. Seeniraj, and R. Murugan, “Experimental validation of a mathematical model of a reedvalve reciprocating air compressor from an automotive-braking system,” International Journal of Automotive Technology, vol. 11, no. 3, pp. 317–322, 2010. Y. K. Fu, J. Liu, and J. Wang, “3D flow simulation and erosion mechanism analysis of the high-pressure throttle valves,” Oil Field Equipment, vol. 39, no. 7, pp. 5–8, 2010. Q. Y. Liu, D. W. Shan, and L. X. Chen, “Study on oscillation of valve stem induced by fluid flow in popper valve,” Oil Field Equipment, vol. 36, no. 5, pp. 24–27, 2007.
