Next Article in Journal
Prediction of Phase Equilibrium Conditions and Thermodynamic Stability of CO2-CH4 Gas Hydrate
Next Article in Special Issue
Research on the Outburst–Rockburst Coupling Disaster Law Based on True Triaxial Unloading Tests
Previous Article in Journal
Research and Development of Inland Green and Smart Ship Technologies in China
Previous Article in Special Issue
Experimental Study on Dynamic and Static Combined Dynamics of Temperature–Water-Coupled Sandstone and Energy Consumption Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on the Vertical Stability of Drilling Wellbore under Optimized Constraints

1
School of Civil Engineering and Architecture, Anhui University of Science and Technology, Huainan 232001, China
2
Engineering Research Center of Underground Mine Engineering, Ministry of Education, Anhui University of Science and Technology, Huainan 232001, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(6), 2317; https://doi.org/10.3390/app14062317
Submission received: 8 February 2024 / Revised: 5 March 2024 / Accepted: 6 March 2024 / Published: 9 March 2024
(This article belongs to the Special Issue Mechanics, Damage Properties and Impacts of Coal Mining)

Abstract

:
With the development of coal resource extraction and wellbore construction proceeding towards deeper depths, the stability of drilling wellbore structures has become increasingly severe, even posing a barrier to the use of drilling method technology in deep wellbore construction. To address this issue, this study raised an optimized constraints method involving pre-throwing cement slurry to the bottom before wellbore decent, altering bottom constraints. Firstly, the critical depth and instability criterion of this optimized method was derived by catastrophe theory. Subsequently, the role of single-factor and multi-factor sensitivity analyses on critical depth was discussed. The engineering effects of optimized constraint methods were contrasted and examined in several drilling projects. Finally, the characteristic values of real engineering were computed using numerical techniques and ABAQUS2020 software, and the efficacy of optimization approaches was examined and validated. The results revealed that the critical depth increased by 41.39 ± 5%. The influence factors described in order of the degree were the counterweight water height, the elastic modulus, the thickness of the wellbore, and the self-weight of the wellbore, sequentially. The conclusion on structural stability between the numerical calculation solution and theoretical calculation solution was completely the same. The optimized constraints method can effectively improve the stability of the wellbore structure.

1. Introduction

As one of the major fossil energy sources, coal has brought a huge supply to power, steel, chemical, and other industries with the growing global demand for energy [1,2]. Coal mine wellbore construction remains an important project at present. However, during the process of drilling method construction, after the drilling rig has drilled the well hole, the builders will measure the deviation of the well hole and then place the prefabricated bottom of the well wall and each section of the well wall into the well hole, then connect them with flanges in turn. The resulting well wall structure is a round tubular structure, which is referred to as the wellbore structure by many scholars [3,4,5]. The wellbore structure is prone to longitudinal instability during the construction process in drilling methods. Once the wellbore structure loses its stability, it can lead to collapse accidents, causing significant personnel and economic losses [6,7,8,9,10]. Moreover, the depth of coal resource extraction in western China is gradually increasing, and the construction depth of new shafts is about to reach a historic height [11,12]. Therefore, the stability of drilling wellbore structures is an urgent engineering problem to be solved.
The drilling method is an advanced directional drilling technology that is widely used in energy extraction, mineral exploration, and rescue [13,14,15,16]. In general, the drilling method is divided into three stages: floating, floated to the bottom but not filled and fixed, and filled and fixed. The second stage is the weakest stage, and it is easy for the wellbore to experience instability once there is any minor interference. Many scholars have conducted numerous studies on the stability characteristics of this stage. Hong [17,18] first proposed the problem of structural instability in the wellbore and analyzed structural stability treating it as a slender rod with hinged ends. Niu et al. [19,20] established the critical depth of the equal cross-section wellbore structure on full water and non-full water by the principle of minimum potential energy and energy method. Cheng et al. [21] further established the critical depth of the variable cross-section drilling wellbore structure. Liu et al. [22,23] established the critical depth of the equal cross-section drilling wellbore structure using catastrophe theory considering the sudden and irreversible characteristics of the wellbore instability process. Rong et al. [24] used numerical calculation methods to analyze the vertical stability of the entire drilling process. Cheng et al. [25] monitored the lateral displacement during construction on site. In summary, a relatively large amount of research on drilling wellbore structural stability has been accumulated. However, these studies only focus on the structural stability of traditional construction methods without proposing effective optimization measures to improve structural instability and conducting targeted research on optimization measures.
In this study, an optimization measure will be raised, and a targeted theoretical and numerical study on optimization measures will be carried out. Considering that wellbore structure is most prone to instability during the moment of suspended sinking and touching the bottom, the contact conditions at the bottom of the wellbore will directly affect structural stability controlling. So, we raise the idea of pre-throwing cement slurry to the bottom before wellbore falls to the bottom. In this optimized method, cement mortar is injected into the bottom by grouting pipe before the wellbore sinks to the bottom. The injection amount and the speed of counterweight water control the speed of the wellbore sinking. When the wellbore sinks to the bottom under the action of counterweight water, the cement has already been pre-thrown. At this time, the bottom of wellbore structure is placed on the pre-casted cement and firmly restricts the displacement in all directions, forming a fixed contact at the bottom. This optimized construction measure improves the structural stability of wellbore structure at the moment when the suspension sinks to the bottom by changing the constraint state from hinged to fixed.
In order to further study and verify the rationality and effectiveness of the optimized constraints method, a series of theoretical and numerical analysis studies were conducted. In terms of theoretical research, the mechanical characteristics of the wellbore structure under optimization measures were analyzed and obtained the total potential energy function and catastrophic instability model. Based on catastrophe theory, the sufficient and necessary conditions for wellbore instability were analyzed and obtained the critical depth. Based on current drilling engineering, the changes in critical depth before and after optimization measures were compared. The influencing factors and degree of critical depth after using optimization measures were also analyzed. In terms of numerical research, a numerical calculation method for vertical stability characteristic values was established, and the finite element model was established. Optimization measures were calculated and implemented according to actual engineering. In order to make it easier for readers to quickly understand the technical terms and abbreviations used in the text and to refer to them easily, a list of symbols has been added here to explain what the symbols represent. For more details, please refer to Table 1.

2. Theoretical Methodology

2.1. Force Analysis of Wellbore under Optimized Construction Method

Before the wellbore is suspended and sinks to the bottom, the grouting pipe is used to fill the cement slurry on the bedrock surface so that the bottom of wellbore can fall on the bedrock surface. The wellbore’s external side is wrapped by cement slurry, and the bottom of the wellbore is in contact with the bedrock surface. It can be assumed that the wellbore is regarded as a slender pressure rod with an upper hinge support and a lower end fixed support. For the convenience of this discussion, the wellbore is regarded as the whole of the continuity of a homogeneous nature, and we made the following basic assumptions:
(1)
The material of the wellbore is completely elastic, obeying the law of Huke;
(2)
The displacement and deformation of the wellbore during flexion is small;
(3)
The effects of friction and bonding between the wall and mud are not taken into account.
The drilling wellbore is regarded as a slender rod hinged at the upper end and solidly supported at the lower end, so the deflection curve function of the drilling wellbore y can be expressed as [26]
y = n = 1 C n cos 2 n + 1 2 H π x cos 2 n 1 2 H π x
where n is the number of the half wave number of the flexural curve function; C n is a coefficient to be determined; x is the height at any point on the wellbore; and H is the length of the wellbore. Based on the assumption (2), taking n equal to 1, the deflection curve deformation function of the drilling wellbore y can be obtained as
y = δ cos 3 π x 2 H cos π x 2 H
where δ is the maximum deformation value of the deflection curve.
When the wellbore sinks to the bottom, the forces on the wellbore structure include the counterforce of wellhead support R B , the vertically downward self-weight of the wellbore q c , the lateral pressure of mud on the external surface of the wellbore P m , the lateral pressure of counterweight water on the internal surface of the wellbore P w , and the counterforce on the bottom of the wellbore R A . These forces are shown in Figure 1.
The simplified force analysis of the wellbore is shown in Figure 2.
For the wellbore, its curvature equation is
1 k = y 1 + y 2 3 2
where y is the first derivative of y and y is the second derivative of y . y 2 is much less than 1, so it can be negligible, and the curvature equation is
1 k y
Therefore, the in-plane bending moment M is
M = 1 k = y 1 + y 2 3 2 y
From the assumption (2), it can be seen that the positive angle between the outer normal of any point on the shaft axis and the X-axis α is approximately equal to zero, and the bending deformation is small, so in the simplified calculation, it can be considered that
sin α α d y d x , cos α 1 , α 2 0
The wellbore weight per unit length q c is
q c = π 4 D 2 d 2 γ c 1 ρ + γ s ρ
where D is the outer diameter of the wellbore; d is the inner diameter of the wellbore; ρ is the ratio of reinforcement; and γ c and γ s are the weight of the concrete and rebar, respectively. The direction of q c is vertically downward.
The lateral pressure of mud per unit length on the outside of the wellbore P m is
P m = π 4 D 2 γ m sin α = q m sin α
where γ m is the weight of mud and q m is the gravity of mud per unit length on the outside of the wellbore. The direction of P m is inward to the wellbore along the normal direction of a point.
The lateral pressure of counterweight water per unit length on the inside of the wellbore P w is
P w = π 4 d 2 γ w sin α = q w sin α
where γ w is the weight of counterweight water and q w is the gravity of counterweight water per unit length on the inside of the wellbore. The direction of P w is outward to the wellbore along the normal direction of a point.

2.2. Building the Total Potential Energy Function

From the force analysis of the wellbore in Section 2.1, the total potential energy equation of the drilled wellbore structure is obtained [27]
        Π = U + V
where   U is the elastic strain energy released by the vertical flexural deformation of the structure and V is the total external potential energy, which consists of the following components: W c , the change in the external potential energy of the wellbore structure caused by the self-weight;   W m , the change in the external potential energy due to mud on the outside of the wellbore; and   W w , the change in the external potential energy due to counterweight water on the inside of the wellbore. Since no displacement occurs under the function of forces R A and R B , the potential energy generated by the function of forces R A and R B is zero. Therefore, the total potential energy equation for the wellbore structure   Π can be written as
    Π = U + W c + W m + W w
The elastic strain energy of wellbore cylinder is
U = 1 2 0 H E I y 2 d x = 41 E I π 4 32 H 3 δ 2
where E represents the elastic modulus of the material used for the wellbore structure. Since the wellbore structure can be either fully reinforced concrete or steel plate concrete, the elastic modulus of the wellbore structure needs to be determined based on the structural form of the wellbore. For commonly encountered reinforced concrete wellbore structures, the elastic modulus of the material used for the wellbore structure E is the composite elastic modulus of the steel bar and concrete converted according to the reinforcement ratio; I is the sectional moment of inertia of the wellbore.
The external potential energy change due to self-weight is
W c = 1 2 0 H q c H x y 2 d x = 1 5 16 π 2 q c δ 2
The change in external potential energy due to the lateral pressure of mud is
W m = W m x + W m y
W m x = 0 H P m x cos α d x = 0 H q m sin α cos α d x = 0
( W m ) y = 0 H P m sin α x 1 2 y 2 d x = 0 H q m ( sin α ) 2 x 1 2 y 2 d x = q m δ 4 π 2 183 π 2 256 H 2 + 308 75
The change in external potential energy under the lateral pressure of counterweight water is
W w = W w x + W w y
W w x = 0 H w P w cos α y d x = 0 H w q w sin α cos α y d x = 2 q w sin 2 π H w 2 H sin 2 π H w H δ 2
( W w ) y = 0 H w P w sin α x 1 2 y 2 d x = 0 H w q w ( sin α ) 2 x 1 2 y 2 d x = π 4 q w δ 4 32 H 2 1 π 2 [ 17,657 1200 cos π H w H 9913 600 cos 2 π H w H + 9 ( 1267 120 + 549 1200 cos 3 π H w H + 1 100 cos 4 π H w H + 1 8 cos 5 π H w H ) ] sin 2 π H w 2 H 15 π ( 59 15 sin π H w H + 83 48 sin 2 π H w H 37 30 sin 3 π H w H + 27 80 sin 4 π H w H 9 50 sin 5 π H w H + 9 80 sin 6 π H w H ) H w H + 183 8 ( H w H ) 2
set
A = 1 π 2 [ 17,657 1200 cos π H w H 9913 600 cos 2 π H w H + 9 ( 1267 120 + 549 1200 cos 3 π H w H + 1 100 cos 4 π H w H + 1 8 cos 5 π H w H ) ] sin 2 π H w 2 H 15 π ( 59 15 sin π H w H + 83 48 sin 2 π H w H 37 30 sin 3 π H w H + 27 80 sin 4 π H w H 9 50 sin 5 π H w H + 9 80 sin 6 π H w H ) H w H + 183 8 ( H w H ) 2
Then, Equation (19) can then be simplified as
( W w ) y = A δ 4 π 4 32 H 2 q w
According to Equations (12)–(21), the total potential energy function of the system can be obtained
Π = A π 4 q w 32 H 2 + q m π 2 183 π 2 256 H 2 + 308 75 δ 4 + 41 E I π 4 32 H 3 + 1 5 16 π 2 q c 2 q w sin 2 π H w 2 H sin 2 π H w H δ 2 + 0 δ
An equivalent transformation of the Equation (22) yields
A 1 = 0
A 2 = 41 E I π 4 32 H 3 + 1 5 16 π 2 q c 2 q w sin 2 π H w 2 H sin 2 π H w H
A 4 = A π 4 q w 32 H 2 + q m π 2 183 π 2 256 H 2 + 308 75
Then, set
x = 4 4 A 4 δ
m = A 2 A 4
n = A 1 4 A 4 4
Then, the total potential energy function of the system can be written as
Π = 1 4 x 4 + 1 2 m x 2 + n x
which shows that the form of the potential function satisfies the standard form of the cusp catastrophe model [28,29]. In Equation (29), x , m , and n are parameters in the standard form of the total potential energy function of the cusp catastrophe model, x is the state variable of the system, and m and n are two control variables of the system. Set the first-order derivative of the total potential energy function equal to zero, then
Π x = x 3 + m x + n = 0
Those points whose first-order derivative of the total potential energy function is zero are singularities, also known as critical points. The surface formed by the set of critical points is called the equilibrium surface, as shown in Figure 3.
However, critical points do not necessarily guarantee that the system is stable. Only the point at which the total energy function takes a unique extreme value is stable, so the set of non-isolated singularities also requires that the second-order derivative of Equation (29) is also zero, so
Π x = 3 x 2 + m = 0
From Equations (30) and (31), we can obtain the equation of bifurcation points set B
Δ = 4 m 3 + 27 n 2 = 0

2.3. Derivation of Equations for Critical Depth

2.3.1. Necessary Condition

It has been seen from Figure 3 that the structure is in an unstable state when m 0 , which means that the points appear inside the bifurcation points set B. For Equation (22), when Δ = 0 , it has two identical roots out of three roots, which is the critical state for the system to mutate. When Δ > 0 , it has only one real root and two complex roots, at which time we consider the system to be stable. When Δ < 0 , it has three real roots, at which time the system will exceed the critical state and then mutate. Using the discriminant Δ , we can also determine whether the system is stable or not. From Equations (27), (28) and (32), it can be determined that the destabilization condition of the wellbore structure is
Δ = 4 m 3 + 27 n 2 = 4 A 2 A 4 3 0
from which it can be determined that
A 2 0
By substituting Equation (24) into Equation (34), we can obtain
H 41 E I π 4 ( 10 π 2 32 ) q c + 64 q w sin 2 π H w 2 H sin 2 π H w H 3
Equation (35) is a necessary condition for wellbore structure instability. The results show that vertical wellbore structure instability will occur as long as the depth of the wellbore is within its range. In addition, when Equation (35) obtains an equal sign, it is the critical condition for shaft wall structure instability, so the critical depth of shaft wall instability is
H c r = 41 E I π 4 ( 10 π 2 32 ) q c + 64 q w sin 2 π H w 2 H sin 2 π H w H 3

2.3.2. Sufficient Condition

If a wellbore is destabilized, its total potential energy function equation satisfies the standard form of the cusp catastrophe model. To meet the condition
Δ = 4 m 3 + 27 n 2 0
m must be negative. Combined with Equation (27), it can be seen that
A 2 0
by substituting Equation (24) into Equation (38), we determine that the calculation result is the same as Equation (35), so Equation (35) is a sufficient condition for the occurrence of the vertical instability of the wellbore. The critical depth of wellbore destabilization is also the same as Equation (36).

2.3.3. Equation for Critical Depth of Wellbore Destabilization

According to Section 2.3.1 and Section 2.3.2, the sufficient and necessary conditions for vertical instability to occur are that the depth of the wellbore is within the range shown in Equation (35), and Equation (36) is the critical depth for the occurrence of vertical instability of the wellbore.

3. Analysis and Discussion

3.1. Comparison Analyses

In this part, several representative drilling wellbore projects were selected, and the critical depth values before and after using optimization measures were calculated. The critical depth before adopting optimization measures was calculated based on reference [23]. The critical depth after adopting optimization measures was calculated by the Equation (36) in this article. The selected representative projects included the Panji west wind wellbore, the Keke gai coal mine wind wellbore, the Banji coal mine auxiliary wellbore, and the Banji coal mine main wellbore. The comparison calculation basis and process are shown in Table 2 and Table 3. The calculation results are shown in Figure 4.
By comparing the calculation results of the critical depth of instability before and after the optimized construction process in Table 3, it can be seen that the wellbore stability changed significantly due to the different constraints under different construction processes. The critical depth of the wellbore instability calculated under the optimized construction technology was significantly higher than that under traditional construction technology.
In the theoretical analysis, the critical depths of instability of the Panji west wind wellbore, Keke gai coal mine wind wellbore, the Banji coal mine auxiliary wellbore, and the Banji coal mine main wellbore were previously increased by 45.29%, 46.99%, 42.08%, and 45.06%, respectively, after adopting the optimized measures. The differences are 3.9%, 5.6%, 0.69%, and 3.67%, respectively, compared with the results of the literature [33], which shows that the optimized construction measures are reasonable and reliable to improve the stability of wellbores by around 41.39 ± 5%.

3.2. Discussion on Critical Depth Influence Factors

From Equation (29), it can be seen that the factors affecting the critical depth of wellbores are elastic modulus E , the self-weight of the wellbore q c , the thickness of the wellbore , and the height of counterweight water H w . In this section, a two-by-two combination of single-factor analysis of the effects of the elastic modulus, the self-weight of the wellbore, the thickness of wellbore, and the height of counterweight water on the critical depth was performed by using the control variable method.
Taking the Banji coal mine auxiliary wellbore as an example, the depth of the wellbore is 638.084 m, the minimum height of the counterweight water is 333.3 m, the external diameter is 9.3 m, the internal diameter is 7.6 m, the elastic modulus is 38.15 Gpa, the moment of inertia of the cross-section is 202.65 m4, the self-weight per unit length is 6.6 × 105 N/m, and the weight of counterweight water per unit length is 4.44 × 105 N/m. We use the above parameters as a baseline to distribute the change in the elastic modulus of the wellbore E from 33.65 Gpa to 43.15 Gpa at intervals of 0.5 Gpa, the self-weight of the wellbore q c from 570,000 N/m to 760,000 N/m at intervals of 10,000 Gpa, the thickness of the wellbore from 760 mm to 950 mm at intervals of 10 mm, and the height of counterweight water H w from 333.3 m to 618.3 m at intervals of 15 m to calculate and analyze the change of the critical instability depth of the deep wellbore structure after adopting the optimized construction measure. To briefly analyze the trend of any two influencing factors on the critical depth of wellbore instability H c r , the analytical diagram of influencing factors is arranged and drawn according to Table 4.
As can be seen in Table 4, the combinations  E - E , q c - q c , H w - H w , and  -  in the table are considered duplicates because they are the same influencing factors. In the blank section above the duplicates, the combinations of any two influencing factors are duplicated with all those below the duplicates in this section. It is sufficient to consider analyzing only the following six factor combinations:  E - q c , E - H w , E - H w , q c - H w , q c - , and H w - . The calculation results are shown in Figure 5.
From Figure 5a–c, it can be seen that the values of the critical depths of the wellbore increase with the increase in elastic modulus regardless of any self-weight of the wellbore, the height of counterweight water, and the thickness of the wellbore. Under a different elastic modulus, the critical depths of the wellbore instability decrease with the self-weight of wellbore and increase with the increase in the thickness of the wellbore, and the degree of increase tends to be slow. When the injection heights of counterweight water are 483.3 m, 453.3 m, 468.3 m, 483.3 m, and 498.3 m, the critical depths of the wellbore instability increase with the increase of counterweight water height, and the degree of increase tends to become urgent.
From Figure 5d,e, it can be visualized that under different heights of counterweight water, the effects of the self-weight of the wellbore on the critical depths are all the same, and the critical depths of the wellbore decrease with the increase in the self-weight of the wellbore. Under different self-weights of the wall, the effects of the thickness of the wellbore on the critical depths are also the same, and the critical depths of the wellbore increase with the growth in the thickness of the wellbore.
As can be seen from Figure 5f, the critical depth of wellbore instability first decreases and then increases with the increase in the height of counterweight water, starting from the minimum counterweight water of 333.3 m. When the height of counterweight water is between 375 m and 400 m, the critical depth reaches the lowest point. When the height of counterweight water is injected above the lowest point, the critical depth increases sharply at first and then slowly with the increase in the height of counterweight water.

3.3. Multi-Factor Sensitivity Discussion

The sensitivity coefficient represents the relative measure of the degree of influence of uncertainty factors on evaluation indicators [34,35]. The formula for calculating the sensitivity coefficient is s A F = Δ A A Δ F F , where Δ A A is the corresponding rate of change of the evaluation indicator A , and Δ F F is the rate of change of the uncertainty factor F . The larger the absolute value of the sensitivity coefficient, the more sensitive the evaluation indicator A is to the uncertainty factor F . It can be seen from Figure 5 that the critical depth of the wellbore structural instability after the adoption of optimized construction measures is affected by several factors to varying degrees. The critical depth of instability of the wellbore structure after the optimization measures was affected by some factors to different degrees. The sensitivity coefficients of critical depth to each factor are shown in Figure 6.
Figure 6 shows that in descending order, the effect degree of each factor on the critical depth of the wellbore instability is the elastic modulus (0.333) > the height of counterweight water (0.278) > the self-weight of the wellbore (0.248) > the thickness of the wellbore (0.118). Additionally, the sensitivity coefficients of critical depth of the wellbore instability to the elastic modulus, height of counterweight water, and thickness of the wellbore are all positive, while the sensitivity coefficient of the self-weight of the wellbore is negative, indicating that the critical depth of the wellbore instability changes in the same direction as the elastic modulus, height of counterweight water, and thickness of the wellbore and changes in the opposite direction of the self-weight of the wellbore.

4. Numerical Calculation Verification

4.1. Principles of Stability Numerical Calculation and Numerical Models

Euler stability was determined to be the cause of the wellbore’s stability issue. The generalized eigenvalue solution of Euler stability in numerical analysis and the incremental solution of extreme stability were used to determine this stability. Literature [36] provides a detailed explanation of the numerical calculation process for the stability of wellbore structures. Based on numerical calculation principles, the finite element software ABAQUS can be used to calculate numerical simulation eigenvalue and analyze the structural stability. When the eigenvalue λ satisfies λ < 0 or λ > 1, the wellbore is vertically stable, and when 0 < λ < 1, the wellbore will be unstable.
The numerical simulation section of this paper is based on the Banji coal mine auxiliary wellbore as the engineering basis and establishes a physical model of the second stage of the wellbore construction process, which is the stage where the wellbore floated to the bottom but is not filled and fixed. Since a wellbore is cylindrical in shape, with constant inner and outer diameters, it can be regarded as an axisymmetric model. The establishment of the entire wellbore structure model is done through the input file. According to the characteristics of the wellbore structure, the entire wellbore is divided into two layers: the lower layer is the bottom of wellbore, and the upper layer is the joint of wellbore. The establishment process of the wellbore is shown below.
(1)
Definition of nodes and node sets. The node types are defined by *Node and Nset. A total of 15,433 nodes are defined in the entire wellbore model, and the height of wellbore model is 638.084 m. The nodes at the bottom of the wellbore are named “BOTTOM”, and the nodes at the top of the wellbore are named “TOP”.
(2)
Definition of elements and elements sets. In this structure, the elements in the bottommost of the wellbore are defined using three-node triangular linear film strain linear shell elements, referred to as S3R. All other elements are specified to use four-node quadrilateral linear reduced integration shell elements, referred to as S4R.
(3)
Define surfaces. Surfaces are defined by the *Surface command, which describes the faces composed of elements. For all elements, their inner surfaces are named “In”, and their outer surfaces are named Out.
(4)
Define materials. Materials are defined by the Materials command. The wellbore structure is made of reinforced concrete. The concrete is graded C70. The elastic modulus of concrete is 38.15 GPa, and its Poisson’s ratio is 0.285. The reinforcing steel has an elastic modulus of 210 GPa, and its Poisson’s ratio is also 0.285.
(5)
Define shell parameters. Defined by the Shell section command, the thickness of the shell is 0.85 m.
(6)
Definition of loads. The height of counterweight water is 333.3 m, and the height of mud is 638.084 m.
(7)
Define boundaries. The *Boundary command is used to define boundary conditions. We define the corresponding constraints on BOTTOM and TOP.
(8)
Define buckling. Buckling is defined by the Buckle Step command. We set the number of eigenvalues requested to five.
The finite element calculation model of the equal cross-sectional wellbore, established according to above definitions, is shown in Figure 7a. The types of elements used in the finite element model are depicted in Figure 7b,c.

4.2. Comparative Analysis of Numerical Models

The stabilities of the wellbore structure before and after the optimization of construction measures are compared by numerical calculation techniques in ABAQUS 2020. The two finite element models are consistent in the definition of material parameters, section characteristics, analysis steps and load structures, but the constraint conditions at the bottom of the wellbore have changed. The constraint comparison of the wellbore models before and after optimization is shown in Table 5.
The boundary conditions of the finite element model of the drilling wellbore established under the two constraints in Table 5 are shown in Figure 8.
The cloud diagram obtained after the model was submitted for analysis is shown in Figure 9.
As can be seen from Figure 9, the stress conditions of wellbore structures before and after optimization are similar. The maximum displacement values for both cases are 1 m. However, due to the different locations of maximum deformation, the eigenvalue of the wellbore structure under the traditional method is 0.9978, which is less than 1, indicating that the wellbore structure is unstable. In contrast, the eigenvalue of the optimized wellbore structure is 2.7471, which is more than 1, indicating that the structure is stable. This is in agreement with the theoretical calculations of the wellbore structure’s stability. Therefore, the optimization measure can significantly enhance the stability of wellbore structures.

5. Conclusions

(1)
In this paper, it is proposed that cement mortar can fill the the bottom of the well in advance to improve the stability of the wellbore before the wall is suspended and sunk to the bottom. According to the contact characteristics of the bottom of the wall after optimization, the wall cylinder is regarded as an equal-section compression rod hinged at the upper end and solidly supported at the lower end, and based on the mechanism of cusp catastrophe characteristics, the critical depth of instability of an equal-section drilling wall in a non-full-water state is deduced from catastrophe theory.
(2)
Compared with the traditional method, the critical depth of the wellbore structure after the optimization method can be increased by 41.39 ± 5%. Through numerical simulation using a sub-well as an engineering example, the eigenvalue of the traditional method is 0.9978, which is smaller than 1, indicating that the wall is unstable, while the eigenvalue of the optimization method is 2.7471, which is larger than 1, indicating that the wall is stable, which is in line with the theoretical calculation. Both the theoretical and numerical analysis reveal that the optimization method is reliable and applicable.
(3)
The critical depth of the wellbore structure by using optimization measures is influenced by multiple factors. It increases with the increase in the elastic modulus and the thickness of the wellbore, decreases with the increase in the self-weight of the wellbore, and shows a trend of first decreasing and then increasing with the increase in the height of counterweight water.
(4)
The critical depth of the wellbore structure by using optimization measures is influenced by multiple factors to varying degrees. The degree to which they are affected in descending order is elastic modulus, counterweight water height, self-weight of the wellbore, and thickness of the wellbore.
(5)
In practical engineering, it is of great practical significance to use high-strength steel with a higher modulus of elasticity or increase the height of appropriate counterweight water and reduce the thickness of internal and external steel plates to improve the stability of the wellbore.

Author Contributions

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

Funding

This research was funded by Anhui Province Key Laboratory of Building Structure and Under-ground Engineering (Hefei), Anhui Jianzhu University (KLBSUE-2022-01).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author (privacy).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhao, S.; Chen, X.; Zhao, X. Analysis of current status and development direction of coal mining technology in China. In Proceedings of the 4th International Conference on Energy Science and Applied Technology (ESAT 2018), Chongqing, China, 29–30 December 2018. [Google Scholar]
  2. Gustavo, A.; Peter, G.; Jürgen, K.; Jaime, R.O. Sustainability of coal mining. is Germany a post-mining model for Colombia? Resour. Policy 2023, 81, 103358. [Google Scholar]
  3. Fares, A.; Kirill, P.; Vitaly, Z.; Sergey, Z. Development of a Software Tool for Visualizing a Mine (Wellbore) in the Industrial Drilling of Oil Wells. Processes 2023, 11, 624. [Google Scholar]
  4. Zou, Y.; Zhu, X.; Wu, X. Experimental Study on Enhancing Wellbore Stability of Coal Measures Formation with Surfactant Drilling Fluid. Chem. Technol. Fuels Oil 2021, 57, 179–187. [Google Scholar] [CrossRef]
  5. Zhang, G.; He, S.; Tang, M.; Kong, L. The mechanism and countermeasures of inclined well wellbore instability in Dibei deep coal seam. J. Pet. Explor. Prod. Technol. 2022, 12, 2833–2848. [Google Scholar] [CrossRef]
  6. Liu, M. Technical study on horizontal coalbed methane wells in high dip strata in the Urumqi Mining Area. China Coalbed Methane 2020, 17, 24–28. [Google Scholar]
  7. Wu, B.; Wang, J.; Zhong, M.; Xu, C.; Qu, B. Multidimensional Analysis of Coal Mine Safety Accidents in China—70 Years Review. Mining Metall. Explor. 2023, 40, 253–262. [Google Scholar] [CrossRef]
  8. Zhang, X.; Zhang, S. Research on the Path of Constructive News of Disaster Report in the Era of Convergence-Media—Taking the Coal Mine Collapse Accident in Alxa Left Banner, Internal Mongolia as an Example. Acad. J. Humanit. Soc. Sci. 2023, 6, 116–122. [Google Scholar]
  9. Ashena, R.; Elmgerbi, A.; Rasouli, V.; Ghalambor, A.; Rabiei, M.; Bahrami, A. Severe wellbore instability in a complex lithology formation necessitating casing while drilling and continuous circulation system. J. Pet. Explor. Prod. Technol. 2020, 10, 1511–1532. [Google Scholar] [CrossRef]
  10. Wang, Y.; Song, L.; He, Y.; Wang, G.; Wang, Q.; He, Y.; Gao, J. Similar simulation test study on the influence of mining on shaft deflection of deep soil based on hydraulic permeability model. Energy Explor. Exploit. 2023, 41, 19–43. [Google Scholar] [CrossRef]
  11. Liu, Z.; Song, C.; Cheng, S.; Hong, W.; Jing, G.; Li, X.; Wang, L.; Zhao, J. Development history and status quo of raise boring technologies and equipment in China. IJCST 2021, 49, 32–65. [Google Scholar]
  12. Jia, D.; Lan, S.; Wang, X. Evolution Characteristics of Concrete Support Strength of the Inclined Shaft Wall in Xiaobaodang Coal Mine. Energy Technol. Manag. 2024, 49, 141–145. [Google Scholar]
  13. Łukasz, B. The analysis of shaft sinking progress as a function of technical and organizational parameters. Multidiscip. Asp. Prod. Eng. 2019, 2, 203–212. [Google Scholar]
  14. Eldert, J.V.; Schunnesson, H.; Saiang, D.; Funehag, J. Improved filtering and normalizing of Measurement-While-Drilling (MWD) data in tunnel excavation. Tunn. Undergr. Space Technol. 2020, 103, 103467. [Google Scholar] [CrossRef]
  15. Zheng, Y.; Zhao, Z.; Zhao, T.; Ma, C.; Zhang, K.; Qi, Y. A new method for monitoring coal stress while drilling process: Theoretical and experimental study. Energy Source Part A 2023, 45, 3980–3993. [Google Scholar] [CrossRef]
  16. Zhang, W.; Li, C.; Jin, J.; Qu, X.; Fan, S.; Xin, C. A new monitoring-while-drilling method of large diameter drilling in underground coal mine and their application. MEAS 2020, 173, 108840. [Google Scholar] [CrossRef]
  17. Hong, B. Revisiting “Axial Stabilization of Drilling Walls in Mud”. IJCST 2008, 33, 121–125. [Google Scholar]
  18. Hong, B. Axial stabilization of the wall of a deep well drilled by the drilling method in mud. IJCST 1980, 9, 22–25+21+68. [Google Scholar]
  19. Niu, X.; Hong, B.; Yang, R. Theoretical study of the vertical structural stability of a borehole wall barrel filled with counterweighted water drilling in mud. IJCST 2005, 30, 463–466. [Google Scholar]
  20. Niu, X.; Hong, B.; Yang, R. Theoretical study of the vertical structural stability of wall barrels in non-full counterweight water drilled wells. JRMGE 2006, 27, 1897–1901. [Google Scholar]
  21. Cheng, H.; Liu, J.; Rong, C.; Yao, Z. Variable cross section shaft drilling lining’s vertical stability in thick alluvium. J. Coal Sci. 2008, 33, 1351–1357. [Google Scholar]
  22. Liu, J.; Liu, Y.; Fu, X. Mechanism analysis of catastrophic instability of drilling shaft lining based on Python language. Coal Sci. Technol. 2022, 50, 75–81. [Google Scholar]
  23. Liu, J.; Cheng, H.; Rong, C.; Wang, C. Analysis of Cusp Catastrophic Model for Vertical Stability of Drilling Shaft Lining. Adv. Civ. Eng. 2020, 2020, 8891751. [Google Scholar] [CrossRef]
  24. Rong, C.; Wang, Y.; Liu, J.; Cheng, H. Numerical simulation on vertical stability critical height of variable cross section drilling shaft in deep alluvium. J. Coal Sci. 2007, 32, 1277–1281. [Google Scholar]
  25. Cheng, H.; Wang, M.; Nie, Q.; Liang, J. Study on mechanical behavior of shaft drilling liner during floating and dropping. J. Beijing Jiaotong Univ. 2012, 36, 48–51. [Google Scholar]
  26. Yang, S. Influence Analysis of Buckling Stability of Super-Long Piles under Different Working Conditions Based on Energy Method. Master’s Thesis, China University of Geosciences, Wuhan, China, 2019. [Google Scholar]
  27. Luo, L.; Zhang, Y. A new method for establishing the total potential energy equations of steel members based on the principle of virtual work. Structures 2023, 52, 904–920. [Google Scholar] [CrossRef]
  28. Yuan, Y.; Li, J. A Cusp Catastrophe Theory Model for Evaluation of Rock Slope Stability. Geol. Explor. 2021, 57, 183–189. [Google Scholar]
  29. Xiao, Y.; Yang, S.; Li, M. A Cusp Catastrophe Model for Alluvial Channel Pattern and Stability. Water 2020, 12, 780. [Google Scholar] [CrossRef]
  30. Bai, W. Experience of drilling mud wall protection in PanSanxifeng wellbore. JCM 1986, 1, 18–20. [Google Scholar]
  31. Fan, J.; Feng, H.; Song, C.; Ren, H.; Ma, Y.; Wang, Q.; Tan, J.; Liu, Q.; Li, C. Key technology and equipment of intelligent mine construction of whole mine mechanical rock breaking in Keke gai Coal Mine. J. Coal Sci. 2022, 47, 499–514. [Google Scholar]
  32. Sun, J. An analysis of vertical shaft forces in drilling well wall cylinders. JCM 2008, 29, 27–30. [Google Scholar]
  33. Zhang, Q.; Liu, Y.; Xu, J.; Gao, L.; An, C. Influence of axially homogenized loads and constraints on the destabilizing length of oil drilling columns. MIPT 2016, 46, 135–142. [Google Scholar]
  34. Gulhan, H.; Dizaji, R.F.; Hamidi, M.N.; Abdelrahman, A.M.; Basa, S.; Cingoz, S.; Koyuncu, I.; Guven, H.; Ozgun, H.; Ersahin, M.E. Modelling of high-rate activated sludge process: Assessment of model parameters by sensitivity and uncertainty analyses. Sci. Total Environ. 2024, 915, 170102. [Google Scholar] [CrossRef] [PubMed]
  35. Takeda, K.; Shimizu, K.; Sato, M.; Katayama, M.; Nakayama, S.M.M.; Tanaka, K.; Ikenaka, Y.; Hashimoto, T.; Minato, R.; Oyamada, Y. Sensitivity assessment of diphacinone by pharmacokinetic analysis in invasive black rats in the Bonin (Ogasawara) Archipelago, Japan. Pestic. Biochem. Phys. 2024, 199, 105767. [Google Scholar] [CrossRef]
  36. Li, Z. Theoretical Analysis of Vertical Stability of Drilling Shaft Lining Based on Catastrophe Theory. Master’s Thesis, Anhui University of Science and Technology, Huainan, China, 2019. [Google Scholar]
Figure 1. Force analysis model of drilling wellbore under non-full water conditions. (a) Longitudinal section force diagram. (b) Top view of force diagram.
Figure 1. Force analysis model of drilling wellbore under non-full water conditions. (a) Longitudinal section force diagram. (b) Top view of force diagram.
Applsci 14 02317 g001
Figure 2. Simplified force analysis diagram of wellbore structure.
Figure 2. Simplified force analysis diagram of wellbore structure.
Applsci 14 02317 g002
Figure 3. Equilibrium planes and control variable planes of the cusp catastrophe model.
Figure 3. Equilibrium planes and control variable planes of the cusp catastrophe model.
Applsci 14 02317 g003
Figure 4. Analysis diagram of critical depth growth after optimization of each wellbore.
Figure 4. Analysis diagram of critical depth growth after optimization of each wellbore.
Applsci 14 02317 g004
Figure 5. Effect of each factor on the critical depth of wellbore. (a) Effect of E and q c on critical depths. (b) Effect of E and H w on critical depths. (c) Effect of E and on critical depths. (d) Effect of q c and H w on critical depths. (e) Effect of and q c on critical depths. (f) Effect of H w and on critical depths.
Figure 5. Effect of each factor on the critical depth of wellbore. (a) Effect of E and q c on critical depths. (b) Effect of E and H w on critical depths. (c) Effect of E and on critical depths. (d) Effect of q c and H w on critical depths. (e) Effect of and q c on critical depths. (f) Effect of H w and on critical depths.
Applsci 14 02317 g005aApplsci 14 02317 g005b
Figure 6. The sensitivity coefficient of the critical depth of wellbore to each influencing factor.
Figure 6. The sensitivity coefficient of the critical depth of wellbore to each influencing factor.
Applsci 14 02317 g006
Figure 7. Finite element model of the entire wellbore and the types of elements used. (a) Numerical model of the complete wellbore. (b) Element type of S3R. (c) Element type of S4R.
Figure 7. Finite element model of the entire wellbore and the types of elements used. (a) Numerical model of the complete wellbore. (b) Element type of S3R. (c) Element type of S4R.
Applsci 14 02317 g007
Figure 8. Comparative boundary condition analysis of wellbore finite element models before and after optimization. (a) Boundary condition of wellbore before optimization. (b) Boundary condition of wellbore after optimization.
Figure 8. Comparative boundary condition analysis of wellbore finite element models before and after optimization. (a) Boundary condition of wellbore before optimization. (b) Boundary condition of wellbore after optimization.
Applsci 14 02317 g008
Figure 9. The displacement cloud diagram of the finite element model for wellbore structure. (a) Under conventional construction method. (b) Under optimized construction method.
Figure 9. The displacement cloud diagram of the finite element model for wellbore structure. (a) Under conventional construction method. (b) Under optimized construction method.
Applsci 14 02317 g009
Table 1. Glossary of symbols and their meanings.
Table 1. Glossary of symbols and their meanings.
SymbolsMeaningsSymbolsMeanings
y The deflection curve function of wellbore   Π The total potential energy of wellbore
y The first derivative of y   U Elastic strain energy
y The second derivative of y   V Total external potential energy
δ The maximum deformation value of the deflection curve   E Elastic modulus of the material used for wellbore
D The outer diameter of wellbore I Sectional moment of inertia of wellbore
d The inner diameter of wellbore A 1 The coefficient of δ
k Radius of curvature A 2 The coefficient of δ 2
M The in-plane bending moment A 4 The coefficient of δ 4
ρ The ratio of reinforcement x The state variable of the system
q c Self-weight of wellbore m A control variable of the system
P m Lateral pressure of mud on the external surface of wellbore n A control variable of the system
P w Lateral pressure of counterweight water on the internal surface of wellbore Π x The first derivative of   Π
R A Counterforce on the bottom of wellbore Π x The second derivative of   Π
R B The counterforce of wellhead
support
  W c The external potential energy change due to self-weight
H Depth of wellbore   W m The change in external potential energy due to the lateral pressure of mud
γ c The weight of concrete   W w The change in external potential energy under the lateral pressure of counterweight water
γ s The weight of steel W m x The horizontal component of   W m
γ m The weight of mud W m y The vertical component of   W m
q m Gravity of mud per unit length on the outside of wellbore. W w x The horizontal component of   W w
γ w The weight of counterweight water W w y The vertical component of   W w
q w Gravity of counterweight water per unit length on the inside of wellbore. Δ Stability criterion
H w Height of counterweight water Thickness of wellbore
Table 2. Critical depth equations for different construction methods.
Table 2. Critical depth equations for different construction methods.
Traditional MethodsOptimized Methods
Restrictive conditionsWith both ends hingedUpper end hinged
lower end fixed support
Sources of informationLiterature [23]Equation (36)
Critical depth formulas H c r 1 = 2 E I π 4 q c π 2 + 4 q w sin 2 π H w H 3 H c r 2 = 41 E I π 4 ( 10 π 2 32 ) q c + 64 q w sin 2 π H w 2 H sin 2 π H w H 3
Table 3. Comparison calculation process of critical depth before and after optimization measures for representative engineering use.
Table 3. Comparison calculation process of critical depth before and after optimization measures for representative engineering use.
Shaft WallsPanji West Wind WellboreKekegai Coal Mine Wind WellboreBanji Coal Mine Auxiliary Wellbore Banji Coal Mine Main Wellbore
SourcesLiterature [30]Literature [31]Literature [32]Literature [23]
E (GPa)35.538.8938.1538.85
I (m4)126.568.26203.33139.66
D (m)97.29.38.3
d (m)767.66.5
q c (N/m)457,008328,159660,000469,000
q w (N/m)376,957276,948444,348351,000
H w (m)222.2215.2333.3291.703
Depths of wellbore H (m)508542.5638.084659.675
Critical depths under the traditional method H c r 1 (m)527.28495.98566.95563.58
Stability under the
traditional method
StableUnstableUnstableUnstable
Critical depths under
the optimized method H c r 2 (m)
766.11729.06805.52817.55
Stability under the
optimized method
StableStableStableStable
Table 4. Influence factor combination analysis matrix.
Table 4. Influence factor combination analysis matrix.
E q c H w
E E - E
q c E - q c q c - q c
H w E - H w q c - H w H w - H w
E - q c - H w - -
Table 5. Comparison of constraints of wellbore models before and after optimization.
Table 5. Comparison of constraints of wellbore models before and after optimization.
Construction TechniquesConstraints on BOTTOMConstraints on TOP
Traditional construction technologyU1, U2, U3U1, U3, UR2
Optimized construction technologyU1, U2, U3, UR1, UR2, UR3U1, U3, UR2
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

Pan, R.; Liu, J.; Cheng, H.; Fan, H. Study on the Vertical Stability of Drilling Wellbore under Optimized Constraints. Appl. Sci. 2024, 14, 2317. https://doi.org/10.3390/app14062317

AMA Style

Pan R, Liu J, Cheng H, Fan H. Study on the Vertical Stability of Drilling Wellbore under Optimized Constraints. Applied Sciences. 2024; 14(6):2317. https://doi.org/10.3390/app14062317

Chicago/Turabian Style

Pan, Ruixue, Jimin Liu, Hua Cheng, and Haixu Fan. 2024. "Study on the Vertical Stability of Drilling Wellbore under Optimized Constraints" Applied Sciences 14, no. 6: 2317. https://doi.org/10.3390/app14062317

APA Style

Pan, R., Liu, J., Cheng, H., & Fan, H. (2024). Study on the Vertical Stability of Drilling Wellbore under Optimized Constraints. Applied Sciences, 14(6), 2317. https://doi.org/10.3390/app14062317

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