Next Article in Journal
Assessment of Hydrogen and Volatile Fatty Acid Production from Fruit and Vegetable Waste: A Case Study of Mediterranean Markets
Previous Article in Journal
Comparison of Different Topologies of Thermal Management Subsystems in Multi-Stack Fuel Cell Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanism of Casing String Curvature Due to Displacement of Surface Strata

1
Faculty of Civil and Environmental Engineering and Architecture, Bydgoszcz University of Science and Technology, Kaliskiego 7, 85-796 Bydgoszcz, Poland
2
Laboratory of Modeling of Damping Systems, Pidstryhach-Institute for Applied Problems in Mechanics and Mathematics of the National Academy of Sciences of Ukraine, Mykytynetska Str. 3, 76002 Ivano-Frankivsk, Ukraine
3
Department of Oil and Gas Well Drilling, Ivano-Frankivsk National Technical University of Oil and Gas, Karpatska Str. 15, 76019 Ivano-Frankivsk, Ukraine
4
Department of Geoinformation Systems and Technologies, Ivano-Frankivsk National Technical University of Oil and Gas, Karpatska Str. 15, 76019 Ivano-Frankivsk, Ukraine
*
Author to whom correspondence should be addressed.
Energies 2022, 15(14), 5031; https://doi.org/10.3390/en15145031
Submission received: 4 June 2022 / Revised: 2 July 2022 / Accepted: 6 July 2022 / Published: 10 July 2022

Abstract

:
One of the main problems of well operation is the risk of uncontrolled leakage of hydrocarbons into the environment. This problem is especially relevant for the long-term operation of wells. The idea for this study was inspired by a real industrial problem that the authors of the article were involved in solving. At several operating gas wells, an abnormal slope of the production tree occurred, which raised the question of the safety of their further operation. An analysis of known studies and current regulatory documents did not allow us to assess the safety of using a gas well based on the measured kinematic parameters of production tree deviations. A mathematical model for the deformation of a package of casing strings when the surface layer of the rock is displaced is developed in the article. A boundary value problem is formulated for differential equations of pipe bending on an elastic foundation. Based on the results of solving this problem, an unambiguous relationship was established between the maximum bending stress in the surface pipe and the angle of inclination of the production tree. The quantitative characteristics of the connection depend on the geometric and mechanical properties of the pipes and on the thickness and mechanical parameters of the formations. It was established that the existing inclination of the production tree can be achieved due to the beginning of the plastic bending of the surface pipe under the slickenside, which does not exclude the exhaustion of the safety margin of the surface pipe and indicates the possible operation of the casing string in a pre-emergency state. In general, the obtained results develop analytical approaches to assessing the behavior of underground structures of a production well in unstable soil bodies.

1. Introduction

In the mining industry, significant attention is paid to safety and production efficiency. In particular, this is reflected in the annual growth of efforts and costs aimed at ensuring the integrity of oil and gas wells [1,2,3]. The life cycle of a well includes several stages, such as design, construction, operation and abandonment. One of the main problems at the stage of well operation is the risk of uncontrolled leakage of hydrocarbons into the environment. This problem is especially relevant for the long-term operation of wells. The human factor has been identified as the most common cause of catastrophic accidents in the oil and gas industry [4]. However, there are threats that are relatively independent of human impact, for example, the destruction of casing strings due to landslides [5,6] or the collapse of deep wells built in unstable formations such as mudstones, bischofites, shales and layers of plastic mineral salt [7,8,9,10,11].
The idea of the presented study was the result of the emergence of a real industrial problem, in the solution of which the authors of the article were directly involved. Several active gas wells experienced an abnormal inclination of the production tree while the wells continued to produce. Because of this, the issue of the safety of further operation of the mentioned wells became extremely acute. An analysis of known studies and current regulatory documents did not allow us to assess the safety of using a gas well according to the parameters of the production tree deviations measured on the daylight surface. It turned out that there are no methods for such an assessment. Our study aims to fill the identified theoretical gap.
The problem of casing string failure due to collapse or shear exists in many oil and gas fields in a number of countries and can lead to emergencies, premature shut-in of wells and significant economic damage [12]. Ensuring the reliability and durability of casing strings is closely related to their strength calculations. Across the range of issues for assessing the strength of casing strings, the most difficult are the issues of calculating collapse resistance to external pressure. Research on this problem is carried out simultaneously in two main areas: theoretical and experimental studies to determine the critical pressure of casing collapse as well as determination of the design external pressure for casing strings [13,14]. It should be recognized that the problem of determining the design external pressure for casing strings has not found a final solution to date. Difficulties in solving this problem lie in the fact that at present there are almost no means for measuring the actual pressures in the casing string at different stages of well construction and operation. Therefore, in the methods of calculating casing strings, when determining the design external pressure, there are always a number of assumptions, and design schemes do not always correspond to the real object of study [15].
In [16], based on experimental and numerical results, expressions were obtained for calculating the shear strength of the casing pipe at various irregularity coefficients. The authors of [17] obtained regularities of the effect of axial length on the casing string failure. Probabilistic estimates of the strength of casing strings are studied in the work [18]. Finite element modeling of various scenarios of static and dynamic loads of casing strings was carried out in [19,20]. A comparative analysis of widely used methods for predicting the destruction of casing pipes under combined loads is presented in [21].
One of the key factors positively influencing the life of a well is the centering of casing strings, which ensures the quality of annulus cementing [22,23,24,25]. The greatest success is achieved through the use of centralizers with non-linear rigidity characteristics, as described in [26,27,28]. Other problems of centering of long objects are discussed in [29,30].
To assess the stress–strain state and limit equilibrium in zones of displacement of the surface strata and in areas of unstable or damaged foundation, it is appropriate to use the methods developed for analyzing the integrity of underground pipelines as extended rod-shell systems [31,32,33,34,35]. In particular, for pipelines built in mountainous areas, such calculations of strength parameters can be carried out according to the kinematic characteristics of the movements of the damaged foundation [36,37]. On the basis of simple schemes for the interaction of rods with the surrounding rock, many other studies involving innovative methods of well construction can be carried out [38,39,40,41,42]. Models of rod [43,44,45] and plate [46,47,48] structures resting on elastic or inelastic Winkler foundations are effective for such engineering simulations. The results of limit state calculations based on 1D models can be refined, taking into account the stress concentration in a pipe thread [49,50,51,52,53] and through the application of fracture mechanics methods to closable crack-like defects in bent and compressed cylindrical shells [54,55,56,57,58,59]. In this case, the model of contact of crack edges along a line has proven itself well [60,61,62,63].
Therefore, the initiation of this study is the inclination of the production tree at the production well Bytkiv-40, located in the Carpathians on the territory of the Bytkiv gas condensate field (Figure 1). The directly measured deviation of the upper part of the surface string from the vertical reaches 6° and is associated with shear processes. Within the area where the well is located, a planned-altitude geodetic survey, ground-based geophysical surveys using electrometric methods [64], as well as a detailed complex interpretation of well logging data for the upper part of the section [65] were performed. The geophysical prerequisites for landslide hazard include anomalies of geophysical parameters, which, by their physical content, are unconditionally associated with the activation of landslides, in particular, zones of low electrical resistance, lithologically consistent with clays, which are potential slickensides with a depth of up to 7 m.
The interpretation of the described problem does not fit into the typical algorithms that are offered to the manufacturers by the current directive documents. Bibliosemantic analysis showed that today there is no method for analytical evaluation of the safe use of underground structures of a production well by the parameters of production tree inclines measured on the daylight surface.
The aim of this work is to explain the possible mechanism of casing string curvature in the surface displacement zone and to develop a model for estimating the maximum stresses from the kinematic parameters of the production tree inclinations.

2. Materials and Methods

2.1. Output Data and Problem Statement

We study the formulated problem of estimating the stress state of casing pipes according to ground measurements of the incline of the production tree on the basis of a well-tested model of the interaction of rods with the base through the Winkler interlayer [66,67].
Direct measuring (Figure 2) set the angle of inclination of the surface casing to the vertical:
θ 0 = arcsin 113 1180 0.096 0.1 rad 5.5
According to geophysical studies, the probable displacement of the upper rock layer occurred above the slickenside at depth b = 7   m .
Let us consider the behavior of a casing string in surface rock layers. We direct the axis O z with the beginning on the daylight surface along the axis of the column deep into the Earth (Figure 3).
A cemented surface casing with a conductor casing is modeled by a piecewise homogeneous rod with tubular cross sections, which in the shear section ( 0 z b ) bends by the forces of plastic resistance of the rock, and below the slickenside ( z > b ) interacts with bedrock according to Winkler’s law with an elastic bed coefficient k . The reinforcing effect of the cement stone on the mechanical behavior of the casing in the cemented areas is not taken into account. At small angles of production tree inclination, the interaction of the conductor with the internal uncemented technical, production and flow strings from above through the wellhead equipment flanges and through contact at depth is neglected. Such assumptions significantly simplify the solution of the problem, transferring it to the category of linear ones.
Let us investigate the distribution of forces, moments, angle of rotation and transverse movement of the column and estimate the magnitude of the maximum stresses.

2.2. Boundary Value Problem of String Bend

The boundary value problem for the system of differential equations of string bending has the following form:
(a) equation in the region:
E J y ′′′′ = m q , z ( 0 , a ) ;
E J 2 y ′′′′ = q , z ( a , b ) ;
E J 2 y ′′′′ + k y = 0 , z ( b , ) ;
(b) boundary and matching conditions:
M ( 0 ) = 0 , Q ( 0 ) = 0 ;
y ( a 0 ) = y ( a + 0 ) , y ( a 0 ) = y ( a + 0 ) ,
M ( a 0 ) = M ( a + 0 ) , Q ( a 0 ) = Q ( a + 0 ) ;
y ( b 0 ) = y ( b + 0 ) , y ( b 0 ) = y ( b + 0 ) ,
M ( b 0 ) = M ( b + 0 ) , Q ( b 0 ) = Q ( b + 0 ) ;
M ( ) = 0 , Q ( ) = 0 .
where
y is a horizontal displacement of the string axis;
q = σ Y D 2 , σ Y is a plastic resistance of the rock;
m = D 1 / D 2 , D 1 , D 2 are outer diameters of the conductor and the surface casing;
a is a conductor length;
b is a depth of the slickenside;
E J 1 π 16 D 1 3 t 1 , E J 2 π 16 D 2 3 t 2 are the bending rigidity of the pipes of the conductor and surface casing, respectively;
t 1 and t 2 are their wall thickness;
E is Young’s modulus of pipe material;
E J = E J 1 + E J 2 = E J 2 ( 1 + n ) is a total rigidity of the conductor with a coaxial surface pipe ( n = E J 1 E J 2 ( D 1 D 2 ) 3 t 1 t 2 , if the pipes are connected perfectly, n = 0 , if the pipes slide mutually, and only the windage of the conductor is taken into account);
M is a bending moment;
Q is a transverse force;
k is a bedrock coefficient.
An analytical study of the boundary value problem (Equations (1)–(7)) was performed using the method of initial parameters [68], which turned out to be effective in simulating casing rod centralizers [69].

3. Results and Analysis

Using the technique of initial parameters [60], solutions of Equations (1)–(3) satisfying conditions (4), (5), and (7) were constructed in the following form:
when z [ 0 , b ]
Q ( z ) = q ( m z ( m 1 ) ( z a ) H ( z a ) ) ,
M ( z ) = q ( m z 2 2 ! ( m 1 ) ( z a ) 2 2 ! H ( z a ) ) ,
y ( z ) = θ 0 q E J 2 ( m 1 + n z 3 3 ! + m n 1 + n z 3 a 3 3 ! H ( z a ) ( m 1 ) ( z a ) 3 3 ! H ( z a ) ) ,
y ( z ) = Δ + θ 0 z q E J 2 ( m 1 + n z 4 4 ! + m n 1 + n ( z 4 a 4 4 ! a 3 ( z a ) ) H ( z a ) ( m 1 ) ( z a ) 3 3 ! H ( z a ) ) ;
when z [ b , )
Q ( z ) = 2 γ q b e z b γ ( ( A + B ) cos z b γ + ( B A ) sin z b γ ) ,
M ( z ) = 2 q b 2 e z b γ ( B cos z b γ + A sin z b γ ) ,
y ( z ) = q b 2 γ E J 2 e z b γ ( ( B A ) cos z b γ ( A + B ) sin z b γ ) ,
y ( z ) = q b 2 γ 2 E J 2 e z b γ ( A cos z b γ + B sin z b γ ) .
where
γ = ( 4 E J 2 k D 2 ) 1 / 4 is a characteristic length parameter that decreases with an increase in the degree of sticking of the string in the bedrock;
k is a coefficient of subgrade reaction (Pa/m);
θ 0 , Δ , A , B are arbitrary constants;
H ( z ) is the Heaviside function.
Satisfying the matching conditions of (6), we obtain
A = 1 4 ( m ( 1 + 2 γ b ) ( m 1 ) ( 1 α ) ( 1 α + 2 γ b ) ) ,
B = 1 4 ( m ( m 1 ) ( 1 α ) 2 ) ,
θ 0 = q b 3 3 ! E J 2 ( 3 ! γ b ( B A ) + m 1 + n + m n 1 + n ( 1 α 3 ) ( m 1 ) ( 1 α ) 3 ) ,
Δ b = θ 0 q b 3 4 ! E J 2 ( 4 ! γ 2 b 2 A + m 1 + n + m n 1 + n ( 1 α 4 4 ! α 3 ( 1 α ) ) ( m 1 ) ( 1 α ) 4 ) ,
and α = a / b .
Thus, Equations (8)–(10) represents solution of the problem i at a given load q .
Let us investigate the function of the bending moment to the extremum. From the condition Q = 0 we find a stationary point z = b + γ a r c t g A + B A B .
Therefore,
| M | max = | M | ( z ) = q b 2 2 ( A 2 + B 2 ) e a r c t g A + B A B .
This moment corresponds to the maximum stress in the surface pipe σ max = | M | max D 2 2 J 2 or
σ max = q b 2 D 2 2 J 2 2 ( A 2 + B 2 ) e a r c t g A + B A B .
Excluding from relations of Equations (10) and (11) the value q b 2 , we obtain the relationship between the highest stress in the surface string and the angle of inclination of the production tree θ 0 , which is the only parameter available for direct measurement:
σ max = 3 2 E θ 0 D 2 b F ,
F = 2 2 ( A 2 + B 2 ) e a r c t g A + B A B 3 ! γ b ( B A ) + m 1 + n + m n 1 + n ( 1 α 3 ) ( m 1 ) ( 1 α ) 3 .
In a similar way, a relationship was established between the horizontal displacement of the top of the surface pipe and its angle of rotation:
Δ = b θ 0 G ,
G = 1 1 4 4 ! γ 2 b 2 A + m 1 + n + m n 1 + n ( 1 α 4 4 ! α 3 ( 1 α ) ) ( m 1 ) ( 1 α ) 4 3 ! γ b ( B A ) + m 1 + n + m n 1 + n ( 1 α 3 ) ( m 1 ) ( 1 α ) 3 .
Numerical studies of the solution were performed for specific parameters of the casing pipes of the Bytkiv-40 well: E 2 10 11 Pa , D 1 = 529   mm , D 2 = 325   mm , t 1 t 2 10   mm , a = 5 m , b = 7   m . Then n 1.628 , m 4.312 , α = 5 / 7 .
According to the results of Equations (8)–(10), the distribution of the transverse force, bending moment, angle of rotation and lateral displacement were studied (Figure 4). The graphical results of the force and deformation analysis of the interaction of underground structures of a borehole with a soil body are presented in a dimensionless form, since in this case one solved dimensionless option corresponds to a whole group of possible dimensional problems. Therefore, in Figure 4, going from left to right, we have: the results of the force analysis—the distribution of the shear force (Figure 4a) and the distribution of the bending moment (Figure 4b) along the length of the underground well structure; the results of the deformation analysis—the distribution of the angles of rotation of the cross sections of the system (Figure 4c) and the distribution of transverse displacements (Figure 4d). The transverse force is equal to zero on the daylight surface, as a whole it changes nonmonotonically, changes symbol once and has a kink in the area of the slickenside. The bending moment behaves nonmonotonically and has an extremum outside the slickenside. The rotation angle graph has linear and non-linear sections, which are associated with the design feature of the underground part of the well. The angles of rotation over the entire interval of the study retain a constant symbol, which determines the interval of monotonous change in the graph of transverse displacements of the underground structure of the well. The maximum values of the angle of rotation and transverse displacements are reached on the daylight surface.
At different values of the coefficient of subgrade reaction of the layered bedrock, corresponding to the values γ / b = 2 ; 1 ; 0.5 , we found the value of the maximum stress in the surface pipe and the lateral displacement of the surface pipe at the wellhead, achieved at a given angle of inclination: θ 0 0.1 rad 5.5 (Table 1).
It is difficult to establish the value of the coefficient of subgrade reaction of the layered bedrock with a cement layer of long-term operation. Therefore, in the calculations, the value k varied over a wide range. A change of value k by a factor of 16 corresponds to a change of value γ by a factor of two.
As can be seen from Table 1, increasing the degree of sticking of casing pipes in the bedrock at a given θ 0 leads to a significant increase in maximum stresses in the wall of the surface pipe. For real θ 0 0.1 rad , Δ 1   m and k 10 5 10 6 Pa / m stresses within the string may reach up to σ max 350 500 MPa . At the same time, according to [18], the yield strength for casing pipes of strength group D is σ Y = 373 552 MPa .
Thus, production tree inclination θ 0 0.1 rad can be achieved due to the beginning of the plastic bending of the surface pipe under the slickenside, which does not exclude the exhaustion of the margin of safety of the surface pipe and indicates the possible operation of the casing in the limit state, and this is, in turn, a sign of a pre-emergency situation.

4. Conclusions

In order to explain the mechanism of casing string curvature as a result of the displacement of the surface strata and to assess the safety of pipes in unstable ground, a calculation scheme for determining the stress in the casing pipe based on the observed deflection of the ground part of the well structure is proposed for the first time.
A mathematical model was developed for the deformation of a package of casing strings when the upper rock layer is displaced. Based on the results of solving the boundary value problem for the differential equation of bending, a relationship was established between the maximum bending stress in the surface pipe and the slope angle of the production tree. The quantitative characteristics of this connection depend on the geometric and mechanical properties of the pipes and the thickness and mechanical parameters of the formations. It was established that the existing inclination of the production tree can be achieved due to the beginning of the plastic bending of the surface pipe under the slickenside, determined at a depth of 7 m, which does not exclude the exhaustion of the margin of safety of the surface pipe and indicates the possible operation of the casing string in the limit state, and this is, in turn, a sign of pre-emergency situations.
Passive waiting can lead to further plastic bending of the surface casing and the destruction of the integrity of the strings. Therefore, it is recommended to continue instrumental observations of the inclination of the string, perform another survey near the well in order to determine the halo of gas contamination, and carry out engineering work to reduce the slope load on the well. In the future, when designing a well, it will be necessary to provide for the study of the upper part of the section in order to prevent the exploitation of the territory in complicated geological, hydrogeological and geomorphological conditions.
Future studies of the described problem should be developed, taking into account the inelastic properties of the unstable foundation.

Author Contributions

Conceptualization, M.D. and I.S.; methodology, I.S. and M.D.; software, O.M.; validation, M.D., O.M., E.K.; formal analysis, M.D., O.M.; investigation,. I.S., E.K., M.D.; resources, E.K., M.D.; writing—original draft preparation, O.M., I.S.; writing—review and editing, M.D., I.S.; visualization, O.M.; supervision, M.D.; project administration, M.D., E.K.; funding acquisition, M.D., I.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schultz, R.A.; Mutlu, U.; Bere, A. Critical Issues in Subsurface Integrity. In Proceedings of the 50th U.S. Rock Mechanics/Geomechanics Symposium, Houston, TX, USA, 26–29 June 2016. paper number ARMA-2016-037. [Google Scholar]
  2. Khalifeh, M.; Saasen, A. General Principles of Well Barriers. In Introduction to Permanent Plug and Abandonment of Wells; Ocean Engineering & Oceanography; Springer: Cham, Switzerland, 2020; Volume 12. [Google Scholar] [CrossRef] [Green Version]
  3. Yakoot, M.S.; Elgibaly, A.A.; Ragab, A.M.S. Well integrity management in mature fields: A state-of-the-art review on the system structure and maturity. J. Pet. Explor. Prod. Technol. 2021, 11, 1833–1853. [Google Scholar] [CrossRef]
  4. Nwankwo, C.D.; Arewa, A.O.; Theophilus, S.C.; Esenowo, V.N. Analysis of accidents caused by human factors in the oil and gas industry using the HFACS-OGI framework. Int. J. Occup. Saf. Ergon. 2022, 28, 1642–1654. [Google Scholar] [CrossRef]
  5. Abdideh, M.; Hedayati Khah, S. Analytical and Numerical Study of Casing Collapse in Iranian Oil Field. Geotech. Geol. Eng. 2018, 36, 1723–1734. [Google Scholar] [CrossRef]
  6. Ghodusi, F.; Jalalifar, H.; Jafari, S. Analysis of the Casing Collapse in Terms of Geomechanical Parameters and Solid Mechanics. J. Chem. Pet. Eng. 2019, 53, 211–225. [Google Scholar] [CrossRef]
  7. Zhang, B.; Guan, Z.; Zhang, Q.; Han, D. Prediction of sustained annular pressure and the pressure control measures for high pressure gas wells. Pet. Explor. Dev. 2015, 42, 567–572. [Google Scholar] [CrossRef]
  8. Lin, Y.; Deng, K.; Sun, Y.; Zeng, D.; Xia, T. Through-wall yield collapse pressure of casing based on unified strength theory. Pet. Explor. Dev. 2016, 43, 506–513. [Google Scholar] [CrossRef]
  9. Zachary, C.; Bello, O.; Teodoriu, C. Calculation and prediction of casing collapse strength based on a new yield strength acquisition method. J. Nat. Gas Sci. Eng. 2021, 95, 104149. [Google Scholar] [CrossRef]
  10. Kovbasiuk, I.M.; Martsynkiv, O.B.; Femiak, Y.M.; Vytvytskyi, I.I.; Zhdanov, Y.D. Research of the stressed state of saline rocks of Prykarpattia deposits under the influence of thermobaric conditions. In Proceedings of the XIV International Scientific Conference “Monitoring of Geological Processes and Ecological Condition of the Environment”, Kyiv, Ukraine, 10–13 November 2020; European Association of Geoscientists & Engineers: Bunnik, The Netherlands, 2020; Volume 2020, pp. 1–5. [Google Scholar] [CrossRef]
  11. Kovbasiuk, I.M.; Martsynkiv, O.B.; Vytvytskyi, I.I.; Pavlyshyn, L.V.; Zhdanov, Y.D.; Paliichuk, I.I. Problems of casing wells in potassium-magnesium salts of the Dnieper-Donets Rift and ways of their solution. In Proceedings of the 15th International Conference Monitoring of Geological Processes and Ecological Condition of the Environment, Kyiv, Ukraine, 17–19 November 2021; European Association of Geoscientists & Engineers: Bunnik, The Netherlands; Volume 2021, pp. 1–5. [Google Scholar] [CrossRef]
  12. Zhao, W.; Ge, J.; Pathegama Gamage, R.; Li, Y.; Song, Z.; Wang, T. Study of Deformation Law of Casing Local Lateral Collapse Based on the Principle of Virtual Work. Energies 2019, 12, 3717. [Google Scholar] [CrossRef] [Green Version]
  13. Lin, B.; Coe, D.; Harris, R.; Timothy, T. On Some Factors Affecting Casing Collapse Resistance under External Pressure. In Proceedings of the SPE Annual Technical Conference and Exhibition, Denver, CO, USA, 5–7 October 2020. SPE-201352-MS. [Google Scholar] [CrossRef]
  14. Zhang, H.; Sun, B.; Yan, G.; Wang, Z.; Huang, M. Distribution laws and effects analysis of casing external pressure taking elastic parameters matching into account. Petroleum 2016, 2, 108–115. [Google Scholar] [CrossRef] [Green Version]
  15. Deng, K.; Lin, Y.; Liu, W.; Zeng, D.; Sun, Y.; Li, K. Equations to Calculate Casing Collapse Strength Under Nonuniform Load Based on New ISO Model. J. Press. Vessel Technol. 2016, 138, 054501. [Google Scholar] [CrossRef]
  16. Li, D.; Yu, F.; Fan, H.; Wang, R.; Yang, S.; Yan, X. Casing Collapse Strength Analysis under Nonuniform Loading Using Experimental and Numerical Approach. Math. Probl. Eng. 2021, 2021, 5577250. [Google Scholar] [CrossRef]
  17. Kuanhai, D.; Bing, L.; Yuanhua, L.; Yongxing, S.; Jiping, Y.; Zuwen, T. Experimental study on the influence of axial length on collapse properties of N80 casings under non-uniform load. Thin-Walled Struct. 2019, 138, 137–142. [Google Scholar] [CrossRef]
  18. Gouveia, L.P.; Lima Junior, E.T.; Santos, J.P.L.; Lira, W.W.M.; Anjos, J.L.R.; Oliveira, F.L. Probabilistic assessment of API casing strength in serviceability limit state. J. Pet. Explor. Prod. Technol. 2020, 10, 2089–2104. [Google Scholar] [CrossRef] [Green Version]
  19. Cao, Y.; Dou, Y.; Li, M.; Suo, H. Theoretical and Simulation Study on the Collapse Strength of a Curved Casing in Horizontal Wells. Procedia Struct. Integr. 2019, 22, 84–92. [Google Scholar] [CrossRef]
  20. Dashtaki, B.B.; Lashkaripour, G.R.; Ghafoori, M.; Moghaddas, N.H. Numerical modeling of casing collapse in Gachsaran formation in Sirri-E oilfield in Persian Gulf. J. Pet. Sci. Eng. 2021, 196, 108009. [Google Scholar] [CrossRef]
  21. Brechan, B.; Teigland, A.; Dale, S.; Sangesland, S.; Kornberg, E. Collapse prediction of pipe subjected to combined loads. J. Pet. Sci. Eng. 2020, 191, 107158. [Google Scholar] [CrossRef]
  22. ISO 10427-2:2004; Petroleum and Natural Gas Industries—Equipment for Well Cementing—Part 2: Centralizer Placement and Stop-Collar Testing. American National Standards Institute: Washington, DC, USA, 2007.
  23. Casing Accessories: Cement Integrity Starts Here; Weatherford: Sandnes, Norway, 2016; p. 56.
  24. Liu, G.; Weber, L.D. Centralizer Selection and Placement Optimization. In Proceedings of the SPE Deepwater Drilling and Completions Conference, Galveston, TX, USA, 20 June 2012. SPE-150345-MS. [Google Scholar] [CrossRef]
  25. Kotskulych, Y.S.; Seniushkovych, M.V.; Martsynkiv, O.B.; Vytvytskyi, I.I. Casing string centeringin directional wells. Nauk. Visnyk Natsionalnoho Hirnychoho Univ. 2016, 2, 23–30. [Google Scholar]
  26. Vytvytskyi, I.I.; Seniushkovych, M.V.; Shatskyi, I.P. Calculation of distance between elastic-rigid centralizers of casing. Nauk. Visnyk Natsionalnoho Hirnychoho Univ. 2017, 5, 28–35. [Google Scholar]
  27. Shatskyi, I.; Vytvytskyi, I.; Senyushkovych, M.; Velychkovych, A. Modelling and improvement of the design of hinged centralizer for casing. IOP Conf. Ser. Mater. Sci. Eng. 2019, 564, 12073. [Google Scholar] [CrossRef]
  28. Shatskyi, I.; Velychkovych, A.; Vytvytskyi, I.; Senyushkovych, M. Modeling of nonlinear properties of casing centralizers equipped with axial thrust. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1018, 012003. [Google Scholar] [CrossRef]
  29. Vlasiy, O.; Mazurenko, V.; Ropyak, L.; Rogal, O. Improving the aluminum drill pipes stability byoptimizing the shape of protector thickening. East.-Eur. J. Enterp. Technol. 2017, 1, 25–31. [Google Scholar] [CrossRef]
  30. Balalayeva, E.; Kukhar, V.; Artiukh, V.; Filatov, V.V.; Simonova, O. Computer Aided Simulation of Behavior of Extrusion Press Ram Jointly with Rotary Elastic Compensator. Adv. Intell. Syst. Comput. 2019, 983, 475–488. [Google Scholar] [CrossRef]
  31. Velychkovych, A.S.; Andrusyak, A.V.; Pryhorovska, T.O.; Ropyak, L.Y. Analytical model of oil pipeline overground transitions, laid in mountain areas. Oil Gas Sci. Technol. 2019, 74, 65. [Google Scholar] [CrossRef] [Green Version]
  32. Kryzhanivs’kyi, E.I.; Rudko, V.P.; Shats’kyi, I.P. Estimation of admissible loads upon a pipeline in the zone of sliding ground. Mater. Sci. 2004, 40, 547–551. [Google Scholar] [CrossRef]
  33. Orynyak, I.V.; Bogdan, A.V. Problem of large displacements of buried pipelines. Part 1. Working out a numerical procedure. Strength Mater. 2007, 39, 257–274. [Google Scholar] [CrossRef]
  34. Zhang, J.; Liang, Z.; Han, C.J. Finite element analysis of wrinkling of buried pressure pipeline under strike-slip fault. Mechanika 2015, 21, 31–36. [Google Scholar] [CrossRef] [Green Version]
  35. Zapuklyak, V.; Melnichenko, Y.; Poberezhnyi, L.; Kyzymyshyn, Y.; Grytsuliak, H.; Komada, P.; Amirgaliyev, Y.; Kozbakova, A. Development of main gas pipeline deepening method for prevention of external effects. In Mechatronic Systems 1: Applications in Transport, Logistics, Diagnostics, and Control; Routledge: London, UK, 2021; pp. 75–87. [Google Scholar] [CrossRef]
  36. Shats’kyi, I.P.; Struk, A.B. Stressed state of pipeline in zones of soil local fracture. Strength Mater. 2009, 41, 548–553. [Google Scholar] [CrossRef]
  37. Shatskyi, I.; Perepichka, V.; Vaskovskyi, M. Longitudinal waves in an elastic rod caused by sudden damage to the foundation. Theor. Appl. Mech. 2021, 48, 29–37. [Google Scholar] [CrossRef]
  38. Velichkovich, A.S. Design features of shell springs for drilling dampers. Chem. Petrol. Eng. 2007, 43, 458–461. [Google Scholar] [CrossRef]
  39. Moisyshyn, V.; Levchuk, K. The impact of vibration mechanism’installation place on the process of retrieving stuck drill pipe. Min. Miner. Depos. 2016, 10, 65–76. [Google Scholar] [CrossRef] [Green Version]
  40. Ropyak, L.Y.; Pryhorovska, T.O.; Levchuk, K.H. Analysis of materials and modern technologies for PDC drill bit manufacturing. Prog. Phys. Met. 2020, 21, 274–301. [Google Scholar] [CrossRef]
  41. Bazaluk, O.; Slabyi, O.; Vekeryk, V.; Velychkovych, A.; Ropyak, L.; Lozynskyi, V. A technology of hydrocarbon fluid production intensification by productive stratum drainage zone reaming. Energies 2021, 14, 3514. [Google Scholar] [CrossRef]
  42. Bazaluk, O.; Velychkovych, A.; Ropyak, L.; Pashechko, M.; Pryhorovska, T.; Lozynskyi, V. Influence of Heavy Weight Drill Pipe Material and Drill Bit Manufacturing Errors on Stress State of Steel Blades. Energies 2021, 14, 4198. [Google Scholar] [CrossRef]
  43. Shatskii, I.P.; Perepichka, V.V. Shock-wave propagation in an elastic rod with a viscoplastic external resistance. J. Appl. Mech. Tech. Phys. 2013, 54, 1016–1020. [Google Scholar] [CrossRef]
  44. Shatskyi, I.; Perepichka, V. Problem of dynamics of an elastic rod with decreasing function of elastic-plastic external resistance. In Dynamical Systems in Applications, Proceedings of the DSTA 2017, Lodz, Poland, 11–14 December 2017; Awrejcewicz, J., Ed.; Springer: Cham, Switzerland, 2018; Volume 249, pp. 335–342. [Google Scholar] [CrossRef]
  45. Pelekhan, B.; Dutkiewicz, M.; Shatskyi, I.; Velychkovych, A.; Rozhko, M.; Pelekhan, L. Analytical Modeling of the Interaction of a Four Implant-Supported Overdenture with Bone Tissue. Materials 2022, 15, 2398. [Google Scholar] [CrossRef]
  46. Shats’kyi, I.P.; Makoviichuk, M.V. Contact interaction of the crack edges in the case of bending of a plate with elastic support. Mater. Sci. 2003, 39, 371–376. [Google Scholar] [CrossRef]
  47. Nobili, A.; Radi, E.; Lanzoni, L. A cracked infinite Kirchhoff plate supported by a two-parameter elastic foundation. J. Eur. Ceram. Soc. 2014, 34, 2737–2744. [Google Scholar] [CrossRef]
  48. Aydin, E.; Dutkiewicz, M.; Ozturk, B.; Sonmez, M. Optimization of elastic spring supports for cantilever beams. Struct. Multidiscip. Optim. 2020, 62, 55–81. [Google Scholar] [CrossRef]
  49. Shats’kyi, I.; Lyskanych, O.; Kornuta, V. Combined deformation conditions for fatigue damage indicator and well–drilling tool joint. Strength Mater. 2016, 48, 469–472. [Google Scholar] [CrossRef]
  50. Pryhorovska, T.; Ropyak, L. Machining Error Influnce on Stress State of Conical Thread Joint Details. In Proceedings of the IEEE 8th International Conference on Advanced Optoelectronics and Lasers (CAOL), Sozopol, Bulgaria, 6–8 September 2019; Volume 9019544, pp. 493–497. [Google Scholar] [CrossRef]
  51. Tutko, T.; Dubei, O.; Ropyak, L.; Vytvytskyi, V. Determination of Radial Displacement Coefficient for Designing of Thread Joint of Thin-Walled Shells. In Advances in Design, Simulation and Manufacturing IV; Lecture Notes in Mechanical Engineering Book Series; Springer: Cham, Switzerland, 2021; pp. 153–162. [Google Scholar] [CrossRef]
  52. Kopei, V.; Onysko, O.; Panchuk, V.; Pituley, L.; Schuliar, I. Influence of Working Height of a Thread Profile on Quality Indicators of the Drill-String Tool-Joint. In Proceedings of the Grabchenko’s International Conference on Advanced Manufacturing Processes, Odessa, Ukraine, 7–10 September 2021; Springer: Cham, Switzerland, 2021; pp. 395–404. [Google Scholar] [CrossRef]
  53. Ropyak, O.Y.; Vytvytskyi, V.S.; Velychkovych, A.S.; Pryhorovska, T.O.; Shovkoplias, M.V. Study on grinding mode effect on external conical thread quality. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1018, 012014. [Google Scholar] [CrossRef]
  54. Rong, L.; Tie, Z.; Wu, X.J.; Wang, C.H. Crack closure effect on stress intensity factors of an axially and a circumferentially cracked cylindrical shell. Int. J. Fract. 2004, 125, 227–248. [Google Scholar] [CrossRef]
  55. Shats’kyi, I.P. Closure of a longitudinal crack in a shallow cylindrical shell in bending. Mater. Sci. 2005, 41, 186–191. [Google Scholar] [CrossRef]
  56. Shats’kyi, I.P.; Makoviichuk, M.V. Contact interaction of crack lips in shallow shells in bending with tension. Mater. Sci. 2005, 41, 486–494. [Google Scholar] [CrossRef]
  57. Shats’kyi, I.P.; Makoviichuk, M.V. Analysis of the limiting state of cylindrical shells with cracks with regard for the contact of crack lips. Strength Mater. 2009, 41, 560–565. [Google Scholar] [CrossRef]
  58. Shatskii, I.P.; Makoviichuk, N.V. Effect of closure of collinear cracks on the stress-strain state and the limiting equilibrium of bent shallow shells. J. Appl. Mech. Tech. Phys. 2011, 52, 464–470. [Google Scholar] [CrossRef]
  59. Dovbnya, K.; Shevtsova, N. Studies on the stress state of an orthotropic shell of arbitrary curvature with the through crack under bending loading. Strength Mater. 2014, 46, 345–349. [Google Scholar] [CrossRef]
  60. Shatsky, I.P. Bending of the plate weakened by the crack with contacting edges. Dopovidi Akademii Nauk Ukrainskoi Rsr Seriya A-Fiziko-Matematichni Ta Technichni Nauki 1988, 7, 49–51. [Google Scholar] [CrossRef]
  61. Shatskii, I.P. Contact of the edges of the slit in the plate in combined tension and bending. Mater. Sci. 1989, 25, 160–165. [Google Scholar] [CrossRef]
  62. Young, M.J.; Sun, C.T. Influence of crack closure on the stress intensity factor in bending plates—A classical plate solution. Int. J. Fract. 1992, 55, 81–93. [Google Scholar] [CrossRef]
  63. Shatskii, I.P. Model for contact of crack boundaries in a bending plate. J. Math. Sci. 2001, 103, 357–362. [Google Scholar] [CrossRef]
  64. Kuzmenko, E.D.; Maksumchuk, V.Y.; Bagriy, S.M.; Sapuzhak, O.Y.; Chepurnyi, I.V.; Dzoba, U.O. On the relevance of using a complex combination of NIEMFE and EM methods in forecasting rock deformation. In Proceedings of the 17th International Conference on Geoinformatics—Theoretical and Applied Aspects, Kyiv, Ukraine, 14–17 May 2018; Volume 2018, pp. 1–6. [Google Scholar] [CrossRef]
  65. Kuzmenko, E.D.; Bezsmertyi, A.F.; Vdovyna, O.P. Study of Landslide Processes by Geophysical Methods; IFNTUOG: Ivano-Frankivsk, Ukraine, 2009. (In Ukrainian) [Google Scholar]
  66. Barber, J.R.; Ciavarella, M. Contact mechanics. Int. J. Solids Struct. 2000, 37, 29–43. [Google Scholar] [CrossRef]
  67. Popov, V.L. Contact Mechanics and Friction: Physical Principles and Applications, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar] [CrossRef]
  68. Beer, F.; Johnston, E.; DeWolf, J.; Mazurek, D. Mechanics of Materials; McGraw-Hill Education: New York, NY, USA, 2015. [Google Scholar]
  69. Shatskyi, I.; Velychkovych, A.; Vytvytskyi, I.; Seniushkovych, M. Analytical models of contact interaction of casing centralizers with well wall. Eng. Solid Mech. 2019, 7, 355–366. [Google Scholar] [CrossRef]
Figure 1. Production tree incline.
Figure 1. Production tree incline.
Energies 15 05031 g001
Figure 2. Angle θ 0 measurement scheme.
Figure 2. Angle θ 0 measurement scheme.
Energies 15 05031 g002
Figure 3. Calculation scheme of the problem: 1—surface casing, 2—conductor casing.
Figure 3. Calculation scheme of the problem: 1—surface casing, 2—conductor casing.
Energies 15 05031 g003
Figure 4. Diagrams of transverse force, bending moment, angle of rotation and lateral displacement at γ / b = 1 ; z = b = 7 m —slickenside. (a) the distribution of the shear force; (b) the distribution of the bending moment; (c) the distribution of the angles of rotation of the cross sections of the system; (d) the distribution of transverse displacements.
Figure 4. Diagrams of transverse force, bending moment, angle of rotation and lateral displacement at γ / b = 1 ; z = b = 7 m —slickenside. (a) the distribution of the shear force; (b) the distribution of the bending moment; (c) the distribution of the angles of rotation of the cross sections of the system; (d) the distribution of transverse displacements.
Energies 15 05031 g004
Table 1. The results of the calculation of mechanical parameters.
Table 1. The results of the calculation of mechanical parameters.
Coefficient of Subgrade Reaction
k , Pa/m
γ b Maximum Stresses
σ max , MPa
Angle of Inclination
θ 0 , rad
Lateral Displacement on the Surface
Δ , m
4.3 10 3 21500.11.82
6.9 10 4 12930.11.20
1.1 10 6 0.55400.10.95
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dutkiewicz, M.; Shatskyi, I.; Martsynkiv, O.; Kuzmenko, E. Mechanism of Casing String Curvature Due to Displacement of Surface Strata. Energies 2022, 15, 5031. https://doi.org/10.3390/en15145031

AMA Style

Dutkiewicz M, Shatskyi I, Martsynkiv O, Kuzmenko E. Mechanism of Casing String Curvature Due to Displacement of Surface Strata. Energies. 2022; 15(14):5031. https://doi.org/10.3390/en15145031

Chicago/Turabian Style

Dutkiewicz, Maciej, Ivan Shatskyi, Oleh Martsynkiv, and Eduard Kuzmenko. 2022. "Mechanism of Casing String Curvature Due to Displacement of Surface Strata" Energies 15, no. 14: 5031. https://doi.org/10.3390/en15145031

APA Style

Dutkiewicz, M., Shatskyi, I., Martsynkiv, O., & Kuzmenko, E. (2022). Mechanism of Casing String Curvature Due to Displacement of Surface Strata. Energies, 15(14), 5031. https://doi.org/10.3390/en15145031

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