Analysis of the Influence of Downhole Drill String Vibration on Wellbore Stability

: Most studies related to aspects of wellbore stability, such as wellbore breakage, block dropping, and wellbore expansion, revolve around the physicochemical interaction between drilling ﬂuid and surrounding rock, but relevant studies show that drill string vibration during drilling also has a crucial and even decisive inﬂuence on wellbore stability. In order to thoroughly explore the inﬂuence mechanism of drill string vibration on wellbore stability, our research group established a ﬁnite element ﬂexible simulation model of drill string dynamics and used a storage downhole vibration measurement device to collect downhole real drilling vibration data to verify the correctness of the simulation model. Then, based on the critical conditions of wellbore breakage, a wellbore stability evaluation method was established, and the wellbore stability under different drilling parameters and drilling tool combination conditions was evaluated and analyzed. The research results play an important role in revealing the inﬂuence mechanism of drill string vibration on wellbore stability and can provide theoretical guidance for engineering problems such as wellbore instability risk assessment.


Introduction
In the petroleum industry, with the development of drilling technology, the proportion of deep wells and ultra-deep wells in oil drilling has gradually increased, and the phenomenon of drilling collapse and wellbore falling blocks leading to well leakage or jamming occurs from time to time and will lead to serious safety problems.According to a survey, the global oil industry spends more than USD 6 billion annually on wellbore stabilization, and the time lost due to wellbore instability accounts for 40% of the non-production time of drilling [1,2].It can be seen that the stability of the wellbore during the drilling process is not only related to the drilling efficiency but also has an important impact on the economy and safety of drilling [3,4].Therefore, it is of great engineering significance and economic value to study the mechanism of wellbore instability and the influencing factors and to propose measures to ensure the stability of the wellbore.
During the drilling process, the causes of the wellbore instability are usually divided into two categories [5]: one is the physicochemical interaction between the rock and the drilling fluid, which is particularly prominent in shale; the second is mechanical reasons, including the collision of drill strings on the wellbore, the difference in rock mechanics, and the distribution change in formation pressure.In 2003, Moheb A. Fam et al. [6] studied the coupling of physical and chemical properties leading to changes in the mechanical properties of mudstone, analyzed the process of mudstone diffusion under coupling, and speculated to explain the complex influencing factors of mud shale borehole stability under actual working conditions.In the same year, Oort et al. [7] conducted a study of shale wellbores, revealing complex links between transport processes (e.g., hydraulic flow, that the initial configuration of BHA has a great influence on the vortex trajectory of the drill string.In 2009, Karkoub et al. [25] combined the proportional integral derivative and the lead-lag controllers with a genetic algorithm to design and optimize the controller, and the results showed that the designed controller could adjust the speed of the drilling column system to the safe range in time and shorten the drilling column vibration settling time, thus achieving the effect of reducing the vibration of the drilling column.In 2011, Liao et al. [26] explored the influence of speed and wellbore friction coefficient on the whirl of BHA and BHA vortex, established a BHA concentrated parameter model, and gave the optimal friction coefficient of the drill string during steady-state motion at the contact point.In 2012, Liao et al. [27] investigated the kinematic parameters of the drill column through a research method of numerical simulation and experimental mutual verification and found that the constructed reduced-order model could identify the characteristics of the collision between the drill column and the wellbore, and a small change in the rotational speed would have a large impact on the kinematic state of the drill column in the state of rotational speed and mass imbalance.From 2013 to 2014, Zhu et al. [28,29] established the collision model between the drill string and the wellbore and found that the impact force of the drill string had a greater influence on the stability of the wellbore.In 2017, Vijayan et al. [30] connected two disks with concentrated BHA mass through massless springs to establish a discrete mass-spring system to analyze BHA reverse whirl.The results show that the reverse whirl of BHA is greatly affected by the speed and axial force, and adjusting the ground speed and large hook load can weaken the reverse whirl of BHA, thereby improving the stability of the wellbore.In the same year, Khaled [31,32] and others developed a new model for predicting wellbore stability, and they found that whirl is directly related to wellbore instability, and collision with wellbore during drill string vibration is an important factor affecting wellbore stability.In 2018, Kapitaniak [33] used a new experimental rig to study the forward and reverse whirl of BHA and, for the first time, experimentally demonstrated the coexistence of the forward and reverse whirl of BHA.In addition, they established a double-degree-of-freedom mathematical model to describe the motion of the bottom plane, calibrated it through experimental results, and further studied it by numerical analysis.In 2020, Zheng et al. [34] proposed an observer-based control scheme for the continuous pole configuration of time-lag systems, which was verified to be effective in suppressing the viscous-slip vibration of the drilling column system.In 2023, Li et al. [35] proposed a set of partial differential equations considering the dynamic effects of continuous tubing drilling based on the beam unit bending theory for longitudinal loads to investigate the contact force between continuous tubing and wellbore in horizontal wells, and the results showed that the contact force between deformed continuous tubing and wellbore was related to the deformation compression rate.
In this paper, the influence of drill string vibration on wellbore stability is mainly studied, and the simulation mathematical model and the theoretical basis of wellbore stability are mainly described in Section 2. In Section 3, the modeling process and boundary conditions of drilling tool combinations are explained.In Section 4, according to the evaluation method of wellbore stability, the influence of different drilling tool combinations on wellbore stability is discussed.The main conclusions are presented in Section 5.

Force Analysis Model of Wellbore under Drill String Vibration
Due to the strong vibration of the downhole drill string, the movement state of the drill string in the downhole is more complicated, so the influence of drill string vibration (mainly lateral vibration) should be considered in the wellbore force model.Through analysis, it is found that the wellbore has constraints on the drill string (see Figure 1), and the current common constraint assumptions are as follows: (1) the wellbore is regarded as a rigid body and does not deform; (2) the wellbore is regarded as an elastomer; (3) the wellbore is regarded as a viscoelastic, and factors such as intrusion depth and damping are taken into account.Due to the variety of motion states of the drill string and the frequent nonlinear contact collision with the wellbore, the first two hypotheses are not suitable, so the third hypothesis method is adopted.rigid body and does not deform; (2) the wellbore is regarded as an elastomer; (3) the wellbore is regarded as a viscoelastic, and factors such as intrusion depth and damping are taken into account.Due to the variety of motion states of the drill string and the frequent nonlinear contact collision with the wellbore, the first two hypotheses are not suitable, so the third hypothesis method is adopted.The expression of the normal force is as follows: where the recovery factor expression is as follows: =− In the formula, -Pre-crash velocity, m/s; -Post-collision velocity, m/s.
The damping factor expression is as follows: The tangential friction expression is as follows: The expression of the normal force F n is as follows: In the formula, where the recovery factor c r expression is as follows: In the formula, The damping factor c f expression is as follows: The tangential friction expression is as follows: In the formula, The expressions of frictional moment M o and frictional resistance F z in the Z direction are as follows:

Solving Method for Wellbore Instability
Professor Mohamed Shafik [32] of Cairo University summarized the previous research by comparing and analyzing the collision force between drilling tools and the wellbore and the strength of the wellbore rock mass, considering the influence of continuous vibration of the drill string on the rock mass strength, which will lead to the failure of the wellbore rock, divided into the following three situations: (1) The drilling tool collides with the wellbore, and the collision stress is large and exceeds the bearing capacity of the wellbore rock mass, resulting in the failure of the wellbore rock mass (indoor cyclic load test (marble)).( 2) The drilling tool collides with the wellbore, and the collision stress is small and does not exceed the bearing capacity of the rock mass of the wellbore, but the multiple collisions between the drill string and the wellbore (the number of collisions exceeds 100 times) lead to a decrease in the strength of the rock mass of the wellbore and failure (thick-walled hollow cylinder cyclic loading (fatigue) test).(3) The collision stress between the drilling tool and the wellbore does not exceed the bearing capacity of the rock mass of the wellbore, but the multiple collisions between the drill string and the wellbore (the number of collisions exceeded 10,000 times) lead to a decrease in the strength of the rock mass of the wellbore and failure (cyclic load experiment (Buria sandstone)).
The criteria for determining whether a rock mass has failed are as follows: F N 0.85 × 0.862S (7) In the formula, [σ]-The rock strength of the wellbore, Pa; S-The contact area when the drill string collides with the wellbore, m 2 ; F N -The minimum force required to destroy a rock, N.
The size of the contact area S can be expressed by the following equation [36]: In the formula, π = 3.14; α-Surrounding angle of drill string; L-Drill string length, m.
In order to take advantage of Equation ( 8), it is necessary to know the size of α and the law of its changes.It has been determined that the analytic formula for the surrounding angle α is affected by certain factors, and its values are calculated as follows: In the formula, l-The chord length of the arc corresponding to the central angle α, m.
Machines 2023, 11, 762 6 of 24 The starting and ending points of the chord are determined by the arc length α in contact with the wellbore relative to the drill string, which can be expressed by Equation (10): In the formula, δ-Clay layer thickness, m; ∆-The width of the ring gap, m.
The instability evaluation coefficient of the wellbore can be obtained using the following equation, which calculates the minimum force required to break the rock due to the collision of the drill string, drill bit, and wellbore, as calculated by Equation (7): In the formula, n-The instability evaluation coefficient of the wellbore; F-The maximum collision force per meter, N; F N -The minimum force required to break the rock, N.
Due to the limitation of simulation time, there is no significant change in collision frequency in a short period.Therefore, it is necessary to make a reasonable prediction of the collision frequency.
In the formula, N-Predict the number of collisions between drill string and wellbore per unit length; N 1 -The number of collisions per unit time per unit contact area; T-The pure drilling time for drilling through the length of the BHA, s. where In the formula, N 2 -The number of collisions between the drill string and the wellbore per unit length in the simulation time; T 1 -Simulation time, s; A-The number of units per contact area; V h -Mechanical drilling speed, m/s; l 1 -The length of the BHA, m. where In the formula, S 1 -Surface area of the inner wall of the wellbore per unit length, m 2 .
According to the contact area calculation formula, the contact area between the drill pipe and the drill collar and the wellbore is calculated, and the minimum force F N required to destroy the rock is calculated according to Equation (7).By combining the derivation of the first section, it can be seen that by obtaining the radial displacement and velocity of the center of mass at any moment on the cross-section of the drilling column, the collision force between the drilling column and the wellbore at that moment can be solved F, which can be further combined with the collision frequency to derive the coefficient of the wellbore destabilization, so the establishment of a nonlinear dynamic model of the drilling column and its solution is especially critical.
In this paper, the nonlinear dynamic equation of the drill string is established based on the theory of finite element beam elements, and the Newmark method is used to solve the equation because the equation and solution method are more classic, and the derivation process will not be repeated.

Basic Assumptions
During the drilling process under the combined action of drilling pressure and torque, the drill pipe is deformed under external loads and randomly collides with the wellbore, so the model is assumed as follows: (1) The deformation of the drilling string is limited, and when the combination of drilling tools is deformed, it is guaranteed to have random nonlinear contact collision with the wellbore; (2) The geometric dimensions and material properties of drill pipes and collars remain constant; (3) The wellbore is a viscoelastic body, the axis of the wellbore and the axis of the drill pipe should coincide in the initial state, there is a ring gap between the wellbore and the drill string, and the wellbore cross-section is circular; (4) Treat the surrounding rock within 1 m of the drill string as a whole; (5) Ignore the influence of drilling fluid on the drill string.

Establishing a Simulation Model
To avoid too long of a simulation time, we made a reasonable simplification of the drilling tool combination (see Figure 2), and then established the drill pipe and drill collar model.Firstly, the soild185 unit was used to establish the drill pipe and drill collar (see Table 1).According to the actual motion of the drill string on site, the degree of freedom of each member of the model was analyzed by finite element analysis software, and constraints were added at the appropriate position of the model (see Table 2).Finally, the other high-order modes were closed, and the first eight-order modes that play a major role in the deformation of the flexible body were taken into effect, and the establishment of the simulation model was completed (see Figure 3).

Added parts
Wellbore and ground Drill pipe and wellbore Drill bits and drill pipes Drill pipe and drill collar

Establishing a Simulation Model
To avoid too long of a simulation time, we made a reasonable simplification of the drilling tool combination (see Figure 2), and then established the drill pipe and drill collar model.Firstly, the soild185 unit was used to establish the drill pipe and drill collar (see Table 1).According to the actual motion of the drill string on site, the degree of freedom of each member of the model was analyzed by finite element analysis software, and constraints were added at the appropriate position of the model (see Table 2).Finally, the other high-order modes were closed, and the first eight-order modes that play a major role in the deformation of the flexible body were taken into effect, and the establishment of the simulation model was completed (see Figure 3).In the actual drilling conditions on site, the drilling operation is mainly controlled by adjusting the drilling pressure and speed.Therefore, in this model, the drilling pressure is mainly applied through the weight of the drilling tool, and a one-way force is added to the top to control the drilling pressure.The top drive is equivalent to a rotary drive with single degrees of freedom to provide the rotational speed.The specific boundary conditions are given in Tables 3 and 4.  In the actual drilling conditions on site, the drilling operation is mainly controlled by adjusting the drilling pressure and speed.Therefore, in this model, the drilling pressure is mainly applied through the weight of the drilling tool, and a one-way force is added to the top to control the drilling pressure.The top drive is equivalent to a rotary drive with single degrees of freedom to provide the rotational speed.The specific boundary conditions are given in Tables 3 and 4. According to the needs of simulation analysis, we created a single-centralizer drilling tool combination and a double-centralizer drilling tool combination concerning the above method and changed the position and size of the centralizer for simulation analysis, which will be explained in Section 4.

Model Reliability Validation
In a high-temperature and high-pressure downhole environment, based on fully investigating the research work of predecessors, our research group designed the highfrequency vibration measurement short section of this experiment (Figure 4) according to the construction situation on-site and used a high-precision triaxial accelerometer and gyroscope to measure and record the vibration of drilling tools.To minimize the influence of the shocker and torsion impactor on their measurement results, the measuring device was chosen to be positioned closer to the drill bit in this example, as shown in Figure 5.To minimize the influence of the shocker and torsion impactor on their measurement results, the measuring device was chosen to be positioned closer to the drill bit in this example, as shown in Figure 5.To minimize the influence of the shocker and torsion impactor on their measurement results, the measuring device was chosen to be positioned closer to the drill bit in this example, as shown in Figure 5.When the vibration measurement device measures the vibration of the drilling tool, the two methods of high-frequency and low-frequency measurement are used at the same time, and the filter setting is carried out considering the influence of clutter generated by the downhole vibration of the drilling tool.The field data in this paper were collected by the No. 3 sensor high-pass filter (see Figure 6).The specific installation parameters are shown in Tables 5 and 6.
When the vibration measurement device measures the vibration of the drilling tool, the two methods of high-frequency and low-frequency measurement are used at the same time, and the filter setting is carried out considering the influence of clutter generated by the downhole vibration of the drilling tool.The field data in this paper were collected by the No. 3 sensor high-pass filter (see Figure 6).The specific installation parameters are shown in Table 5 and 6.

Number of Experiments Sensor Number Sampling Frequency Filter Frequency Error
First trial To verify the accuracy and effectiveness of the simulation model, we used the real  To verify the accuracy and effectiveness of the simulation model, we used the real drilling data to confirm the accuracy of the simulation model.First, we obtained the triaxial acceleration of the real drilling data, and then used the simulation model to simulate the vibration of the drilling tool according to the corresponding mark point on the model according to the position of the sensor installation (see Figure 7).Finally, we compared the simulation results and the real drilling data to analyze the similarity of the two.To verify the accuracy and effectiveness of the simulation model, we used the real drilling data to confirm the accuracy of the simulation model.First, we obtained the triaxial acceleration of the real drilling data, and then used the simulation model to simulate the vibration of the drilling tool according to the corresponding mark point on the model according to the position of the sensor installation (see Figure 7).Finally, we compared the simulation results and the real drilling data to analyze the similarity of the two.The arrangement of the sensor in the measuring device is as follows: Y-axis measurement of longitudinal acceleration, X-axis measurement of tangential acceleration, and Z-axis measurement of radial acceleration.The gyroscope direction of the sensor coincides with the positive direction of X, with a drilling pressure of 16 t and speed of 50 rpm.The intercepted radial and tangential acceleration and simulation data were compared and analyzed, as shown in Figures 8 and 9.
From the measured and simulated radial and tangential acceleration curves, the overall trend is similar, the measured radial acceleration fluctuation range is −6~6 g, and the tangential acceleration fluctuation range is −8~10 g.The simulation radial fluctuation interval is −0.5~0.5 g, and the tangential acceleration fluctuation interval is −0.5~0.5 g.The Pearson correlation calculation is used to show that the correlation coefficient between the measured and simulated radial acceleration is 0.6569 and the tangential acceleration correlation coefficient is 0.6982, which can be seen in the Pearson correlation coefficient table; the two are strongly correlated.
The arrangement of the sensor in the measuring device is as follows: Y-axis measurement of longitudinal acceleration, X-axis measurement of tangential acceleration, and Zaxis measurement of radial acceleration.The gyroscope direction of the sensor coincides with the positive direction of X, with a drilling pressure of 16 t and speed of 50 rpm.The intercepted radial and tangential acceleration and simulation data were compared and analyzed, as shown in Figures 8 and 9.  From the measured and simulated radial and tangential acceleration curves, the overall trend is similar, the measured radial acceleration fluctuation range is −6~6 g, and the tangential acceleration fluctuation range is −8~10 g.The simulation radial fluctuation interval is −0.5~0.5 g, and the tangential acceleration fluctuation interval is −0.5~0.5 g.The Pearson correlation calculation is used to show that the correlation coefficient between the measured and simulated radial acceleration is 0.6569 and the tangential acceleration correlation coefficient is 0.6982, which can be seen in the Pearson correlation coefficient table; the two are strongly correlated.
In the acceleration solution, the actual downhole rotational angular velocity of the drill string can be solved by the tangential acceleration measured by the three-axis accelerometer, and the movement trajectory of the drill string during drilling at the measurement point position of the drill string can be calculated according to the three-axis acceleration.Through the comparison of the trajectory diagram (Figure 10), the analysis results of the measured data and the simulation show that whirl occurs at this point, and the whirl trajectories of the two are similar.The arrangement of the sensor in the measuring device is as follows: Y-axis measurement of longitudinal acceleration, X-axis measurement of tangential acceleration, and Zaxis measurement of radial acceleration.The gyroscope direction of the sensor coincides with the positive direction of X, with a drilling pressure of 16 t and speed of 50 rpm.The intercepted radial and tangential acceleration and simulation data were compared and analyzed, as shown in Figures 8 and 9.  From the measured and simulated radial and tangential acceleration curves, the overall trend is similar, the measured radial acceleration fluctuation range is −6~6 g, and the tangential acceleration fluctuation range is −8~10 g.The simulation radial fluctuation interval is −0.5~0.5 g, and the tangential acceleration fluctuation interval is −0.5~0.5 g.The Pearson correlation calculation is used to show that the correlation coefficient between the measured and simulated radial acceleration is 0.6569 and the tangential acceleration correlation coefficient is 0.6982, which can be seen in the Pearson correlation coefficient table; the two are strongly correlated.
In the acceleration solution, the actual downhole rotational angular velocity of the drill string can be solved by the tangential acceleration measured by the three-axis accelerometer, and the movement trajectory of the drill string during drilling at the measurement point position of the drill string can be calculated according to the three-axis acceleration.Through the comparison of the trajectory diagram (Figure 10), the analysis results of the measured data and the simulation show that whirl occurs at this point, and the whirl trajectories of the two are similar.In the acceleration solution, the actual downhole rotational angular velocity of the drill string can be solved by the tangential acceleration measured by the three-axis accelerometer, and the movement trajectory of the drill string during drilling at the measurement point position of the drill string can be calculated according to the three-axis acceleration.Through the comparison of the trajectory diagram (Figure 10), the analysis results of the measured data and the simulation show that whirl occurs at this point, and the whirl trajectories of the two are similar.Combined with the comparison of tangential and radial acceleration, the accuracy of the nonlinear dynamic simulation model of the drill string and the reliability of the simulation results are verified.Combined with the comparison of tangential and radial acceleration, the accuracy of the nonlinear dynamic simulation model of the drill string and the reliability of the simulation results are verified.

Analysis of the Impact of Borehole Risk Assessment Instability under Drill String Vibration
Using the calculation method introduced in Section 2.2, we calculated the contact area between the drill pipe, the drill collar, and the wellbore, and further calculated the damage force of the wellbore at different frequencies, as shown in Table 7.The dynamic model of the drill string in Figures 1-3 is simulated by a drilling pressure of 16 t and rotational speeds of 20, 50, and 80 rpm, respectively, and the collision force is shown in Figure 11.As shown in Figure 11, the collision between the drilling tool and the wellbore has the situation of one-time failure of rock mass and multiple collisions of rock mass failure.When the speed is 80 rpm, multiple positions can destroy the rock mass at one time, and the collision force of the drilling position (0-157.6m from the drill bit) does not change much overall because the diameter of the drill collar is large, the stiffness is strong, and the collision force changes little in different positions.At the interface between BHA and the drill pipe, the collision force is as high as 113,853 N, far more than the force required for a failure, and the drill pipe part (157.6-657.6 m from the drill bit) generally changes greatly in the collision force and collision frequency, because the diameter of the drill pipe is smaller than the borehole, and the drill pipe accounts for most of the drilling tool combination during the drilling process, and the length is longer.When the speed is 50 rpm, the collision force is reduced compared to 80 rpm, and two positions can cause damage to the rock mass.When the speed is 20 rpm, only local collisions occur, and the collision force is small.
By analyzing the collision between the drilling tool and the wellbore, it is found that As shown in Figure 11, the collision between the drilling tool and the wellbore has the situation of one-time failure of rock mass and multiple collisions of rock mass failure.When the speed is 80 rpm, multiple positions can destroy the rock mass at one time, and the collision force of the drilling position (0-157.6m from the drill bit) does not change much overall because the diameter of the drill collar is large, the stiffness is strong, and the collision force changes little in different positions.At the interface between BHA and the drill pipe, the collision force is as high as 113,853 N, far more than the force required for a failure, and the drill pipe part (157.6-657.6 m from the drill bit) generally changes greatly in the collision force and collision frequency, because the diameter of the drill pipe is smaller than the borehole, and the drill pipe accounts for most of the drilling tool combination during the drilling process, and the length is longer.When the speed is 50 rpm, the collision force is reduced compared to 80 rpm, and two positions can cause damage to the rock mass.When the speed is 20 rpm, only local collisions occur, and the collision force is small.
By analyzing the collision between the drilling tool and the wellbore, it is found that the collision force between the drilling tool and the wellbore will increase with the increase in speed.However, we can observe that the collision force is much more than the force required for a failure, but when the collision force is lower than the force required for a failure, the influence of collision frequency should be specifically considered; thus, it is extremely necessary to explore the collision between the drilling tool and the wellbore.
The dynamic model of the above drill string is further simulated and analyzed with a drilling pressure of 16 t and speed of 50 rpm, and the collision force and collision frequency are shown in Figures 12 and 13.
much overall because the diameter of the drill collar is large, the stiffness is strong, and the collision force changes little in different positions.At the interface between BHA and the drill pipe, the collision force is as high as 113,853 N, far more than the force required for a failure, and the drill pipe part (157.6-657.6 m from the drill bit) generally changes greatly in the collision force and collision frequency, because the diameter of the drill pipe is smaller than the borehole, and the drill pipe accounts for most of the drilling tool combination during the drilling process, and the length is longer.When the speed is 50 rpm, the collision force is reduced compared to 80 rpm, and two positions can cause damage to the rock mass.When the speed is 20 rpm, only local collisions occur, and the collision force is small.
By analyzing the collision between the drilling tool and the wellbore, it is found that the collision force between the drilling tool and the wellbore will increase with the increase in speed.However, we can observe that the collision force is much more than the force required for a failure, but when the collision force is lower than the force required for a failure, the influence of collision frequency should be specifically considered; thus, it is extremely necessary to explore the collision between the drilling tool and the wellbore.
The dynamic model of the above drill string is further simulated and analyzed with a drilling pressure of 16 t and speed of 50 rpm, and the collision force and collision frequency are shown in Figures 12 and 13.According to the schematic diagram of the collision force, it can be seen that the collision force at the drill pipe part is significantly greater than the collision force at the BHA, the collision force at the drill pipe part fluctuates greatly, and the maximum collision force of 57,126 N occurs at a distance of 450 m from the drill bit (on the drill pipe), reaching the According to the schematic diagram of the collision force, it can be seen that the collision force at the drill pipe part is significantly greater than the collision force at the BHA, the collision force at the drill pipe part fluctuates greatly, and the maximum collision force of 57,126 N occurs at a distance of 450 m from the drill bit (on the drill pipe), reaching the force required for one failure.The collision force at BHA fluctuated smoothly between 0 and 34,816.5 N, and the maximum collision force at BHA was 34,816.5 N at 155 m from the drill bit.
Using Equation ( 13), the number of collisions at the same point is reasonably predicted (see Figure 14), and the collision frequency in each contact area is obtained.According to the number of collisions (the number of collisions between the drill collar and the wellbore is generally higher than the number of collisions between the drill pipe and the wellbore), the drilling collar part is more likely to collide with the wellbore in the simulation time than the drill pipe part, which is consistent with the research results that the drilling collar part has a larger diameter and more violent whirl than the upper drill pipe.Among them, the frequency of collision between drilling and the wellbore at 22.6 m from the bottom of the well reached 925 times, indicating that the collision between the drilling collar and the wellbore at the bottom of the well was fierce.According to the schematic diagram of the collision force, it can be seen that the collision force at the drill pipe part is significantly greater than the collision force at the BHA, the collision force at the drill pipe part fluctuates greatly, and the maximum collision force of 57,126 N occurs at a distance of 450 m from the drill bit (on the drill pipe), reaching the force required for one failure.The collision force at BHA fluctuated smoothly between 0 and 34,816.5 N, and the maximum collision force at BHA was 34,816.5 N at 155 m from the drill bit.
Using Equation ( 13), the number of collisions at the same point is reasonably predicted (see Figure 14), and the collision frequency in each contact area is obtained.According to the number of collisions (the number of collisions between the drill collar and the wellbore is generally higher than the number of collisions between the drill pipe and the wellbore), the drilling collar part is more likely to collide with the wellbore in the simulation time than the drill pipe part, which is consistent with the research results that the drilling collar part has a larger diameter and more violent whirl than the upper drill pipe.Among them, the frequency of collision between drilling and the wellbore at 22.6 m from the bottom of the well reached 925 times, indicating that the collision between the drilling collar and the wellbore at the bottom of the well was fierce.
Combined with the collision frequency and collision force, the wellbore instability evaluation table is established by Formula ( 12), as shown in Table 8.Combined with the collision frequency and collision force, the wellbore instability evaluation table is established by Formula ( 12), as shown in Table 8.According to the wellbore instability evaluation method, the safety evaluation of the wellbore is carried out (see Figure 15).It can be seen that the instability coefficient of 1 is exceeded at distances of 240, 265, 380, 375, 450, and 455 m from the drill bit, indicating that the wellbore has been damaged and dropped at these six points.

≥
According to the wellbore instability evaluation method, the safety evaluation of the wellbore is carried out (see Figure 15).It can be seen that the instability coefficient of 1 is exceeded at distances of 240, 265, 380, 375, 450, and 455 m from the drill bit, indicating that the wellbore has been damaged and dropped at these six points.

Influence of Single Centralizer Placement on Wellbore Stability
In the actual drilling conditions, the appropriate drilling tool combination will be selected for different formations, which shows that the influence of drilling tool combination on the stability of the wellbore is more prominent, so we changed the drilling tool combination model and simulation analysis.The schematic diagram of drilling tool combinations is shown in Figure 16.In the actual drilling conditions, the appropriate drilling tool combination will be selected for different formations, which shows that the influence of drilling tool combination on the stability of the wellbore is more prominent, so we changed the drilling tool combination model and simulation analysis.The schematic diagram of drilling tool combinations is shown in Figure 16.Numerical simulation was performed on the placement positions of three types of packed hole centralizers, shown in Figure 16, using a drilling pressure of 16 tons and a rotation speed of 50 rpm, and the force and frequency curves are shown in the following figures.
According to the schematic diagram of the collision frequency curves in Figure 17, on the whole, when the packed hole centralizer is installed 28.55 m from the drill bit, the collision frequency of the whole well is generally lower than that of the other two installation positions.At a distance of 105 m from the drill bit, the maximum collision frequencies of the three centralizers are 775, 538, and 452 times, showing that the collision frequency can be effectively reduced when the centralizer is installed 28.55 m from the drill bit, which is 42% and 16% lower than those installed 10.55 m and 19.55 m from the drill bit, respectively.Numerical simulation was performed on the placement positions of three types of packed hole centralizers, shown in Figure 16, using a drilling pressure of 16 tons and a rotation speed of 50 rpm, and the collision force and frequency curves are shown in the following figures.
According to the schematic diagram of the collision frequency curves in Figure 17, on the whole, when the packed hole centralizer is installed 28.55 m from the drill bit, the collision frequency of the whole well is generally lower than that of the other two installation positions.At a distance of 105 m from the drill bit, the maximum collision frequencies of the three centralizers are 775, 538, and 452 times, showing that the collision frequency can be effectively reduced when the centralizer is installed 28.55 m from the drill bit, which is 42% and 16% lower than those installed 10.55 m and 19.55 m from the drill bit, respectively.
According to the collision force curves of different installation positions of the full eye centralizer (see Figure 18), when the packed hole centralizer is installed 10.55 m from the drill bit, the collision force of the whole well does not exceed the force required for a failure, and the maximum collision force is 62,139 N at 60 m from the drill bit.When the packed hole centralizer is installed 19.55 m from the drill bit, five collision forces exceed the required force for one failure, and the maximum collision force of 84,158 N occurs 90 m from the drill bit.When the packed hole centralizer is installed 28.55 m from the drill bit, four collision forces exceed the required force for one failure, and the maximum collision force of 94,161 N occurs 230 m from the drill bit, which is far more than the required force for one failure.By comparing the three centralizer positions, it is found that the collision force of the packed hole centralizer installed 10.55 m from the drill bit is the smallest, and compared with the other two positions, its maximum collision force is reduced by 26% and 34%.
on the whole, when the packed hole centralizer is installed 28.55 m from the drill bit, the collision frequency of the whole well is generally lower than that of the other two installation positions.At a distance of 105 m from the drill bit, the maximum collision frequencies of the three centralizers are 775, 538, and 452 times, showing that the collision frequency can be effectively reduced when the centralizer is installed 28.55 m from the drill bit, which is 42% and 16% lower than those installed 10.55 m and 19.55 m from the drill bit, respectively.According to the collision force curves of different installation positions of the full eye centralizer (see Figure 18), when the packed hole centralizer is installed 10.55 m from the drill bit, the collision force of the whole well does not exceed the force required for a failure, and the maximum collision force is 62,139 N at 60 m from the drill bit.When the packed hole centralizer is installed 19.55 m from the drill bit, five collision forces exceed the required force for one failure, and the maximum collision force of 84,158 N occurs 90 m from the drill bit.When the packed hole centralizer is installed 28.55 m from the drill bit, four collision forces exceed the required force for one failure, and the maximum collision force of 94,161 N occurs 230 m from the drill bit, which is far more than the required force for one failure.By comparing the three centralizer positions, it is found that the collision force of the packed hole centralizer installed 10.55 m from the drill bit is the smallest, and compared with the other two positions, its maximum collision force is reduced by 26% and 34%.By observing the schematic diagram of the wellbore instability coefficient (see Figure 19), it is found that when the packed hole centralizer is installed 19.55 m from the drill bit, a total of six places exceed the wellbore instability coefficient of 1.When the packed hole centralizer is installed 28.55 m from the drill bit, a total of eight places exceed the wellbore instability coefficient of 1, and when it is installed 10.55 m from the drill bit, there are no points exceeding the instability coefficient of 1 in the whole well.Compared with the placement position of the three types of centralizers, the centralizer is placed 10.55 m away from the drill bit, which can effectively ensure the stability and safety of the wellbore.By observing the schematic diagram of the wellbore instability coefficient (see Figure 19), it is found that when the packed hole centralizer is installed 19.55 m from the drill bit, a total of six places exceed the wellbore instability coefficient of 1.When the packed hole centralizer is installed 28.55 m from the drill bit, a total of eight places exceed the wellbore instability coefficient of 1, and when it is installed 10.55 m from the drill bit, there are no points exceeding the instability coefficient of 1 in the whole well.Compared with the placement position of the three types of centralizers, the centralizer is placed 10.55 m away from the drill bit, which can effectively ensure the stability and safety of the wellbore.Observing the collision frequency diagram of the three centralizers (see Figure 20), the collision frequency of the 327.2 mm centralizer in the BHA part is significantly higher than the 330.2 mm and 329.7 mm centralizers, and the collision frequency reaches a peak of 1592 times at 105 m from the drill bit.The collision frequency of the full-eye centralizer at BHA is relatively reduced for the 327.2 mm centralizer, reaching a peak of 775 times at 105 m from the drill bit.The collision frequency of the 329.2 mm centralizer is slightly smaller than that of other sizes on the whole, reaching a peak of 645 times at 160 m from the drill bit, and its collision frequency is reduced by 59% and 17% compared with the other two sizes.Observing the collision frequency diagram of the three centralizers (see Figure 20), the collision frequency of the 327.2 mm centralizer in the BHA part is significantly higher than the 330.2 mm and 329.7 mm centralizers, and the collision frequency reaches a peak of 1592 times at 105 m from the drill bit.The collision frequency of the full-eye centralizer at BHA is relatively reduced for the 327.2 mm centralizer, reaching a peak of 775 times at 105 m from the drill bit.The collision frequency of the 329.2 mm centralizer is slightly smaller than that of other sizes on the whole, reaching a peak of 645 times at 160 m from the drill bit, and its collision frequency is reduced by 59% and 17% compared with the other two sizes.According to the collision force diagram in Figure 21, it can be seen that the combination of drilling tools with centralizers with packed hole sizes of 330.2 mm, 329.2 mm, and 327.2 mm does not exceed the force required for a failure at BHA, and the combination of drilling tools with a 329.2 mm and 327.2 mm centralizer has a large collision force at the drill pipe and exceeds the force required for a single failure.The combination of drilling tools with a 329.2 mm centralizer has a maximum collision force of 82,789.5 N at 500 m from the drill bit, and the maximum collision force of 70,017 N occurs at 425 m from the drill bit with a combination of drilling tools installed with a 327.2 mm centralizer, which far exceeds the force required for a single failure.According to the collision force diagram in Figure 21, it can be seen that the combination of drilling tools with centralizers with packed hole sizes of 330.2 mm, 329.2 mm, and 327.2 mm does not exceed the force required for a failure at BHA, and the combination of drilling tools with a 329.2 mm and 327.2 mm centralizer has a large collision force at the drill pipe and exceeds the force required for a single failure.The combination of drilling tools with a 329.2 mm centralizer has a maximum collision force of 82,789.5 N at 500 m from the drill bit, and the maximum collision force of 70,017 N occurs at 425 m from the drill bit with a combination of drilling tools installed with a 327.2 mm centralizer, which far exceeds the force required for a single failure.In the schematic diagram of the wellbore instability coefficient in Figure 22, it can be seen that when the packed hole centralizer is installed, the instability coefficient of the whole wellbore is below 1.When the 329.2 mm centralizer was installed, the wellbore instability coefficient exceeded 1 in 20 places, and when the 327.2 mm centralizer was installed, the instability coefficient of the wellbore exceeded 1 in 6 places.It can be seen that the combination of drilling tools installed with 329.2 mm and 327.2 mm centralizers results in unstable wellbore phenomena at both the drill pipe and the drilling collar, so it is recommended to use a packed-hole centralizer.

Comparative Analysis of the Influence of Single and Double Centralizers on the Stability of the Wellbore
The effect of single and double centralizers on the stability of the wellbore is unknown.Thus, the combination of drilling tools, shown in Figure 23, was simulated and studied.
The combinations of single-and double-centralizer drilling tools were simulated with a drilling pressure of 16 t and speed of 50 rpm, and the collision frequency of the two was compared (see Figure 24).The collision frequency of the double-centralizer drilling tool combination reached a peak of 667 times at 80 m from the drill bit, indicating that the bottom movement of the well is complex, the collision frequency of the drill pipe part is relatively stable, and the curve fluctuation is not large.In the combination of singlecentralizer drilling tools, the collision frequency at BHA was higher, reaching a peak of 775 collision frequencies at 105 m from the bottom of the well, and the collision frequency of the double-centralizer drilling tool combination was reduced by 14% compared with the single centralizer.
In the schematic diagram of the wellbore instability coefficient in Figure 22, it can be seen that when the packed hole centralizer is installed, the instability coefficient of the whole wellbore is below 1.When the 329.2 mm centralizer was installed, the wellbore instability coefficient exceeded 1 in 20 places, and when the 327.2 mm centralizer was installed, the instability coefficient of the wellbore exceeded 1 in 6 places.It can be seen that the combination of drilling tools installed with 329.2 mm and 327.2 mm centralizers results in unstable wellbore phenomena at both the drill pipe and the drilling collar, so it is recommended to use a packed-hole centralizer.The effect of single and double centralizers on the stability of the wellbore is unknown.Thus, the combination of drilling tools, shown in Figure 23, was simulated and studied.The combinations of single-and double-centralizer drilling tools were simulated with a drilling pressure of 16 t and speed of 50 rpm, and the collision frequency of the two was compared (see Figure 24).The collision frequency of the double-centralizer drilling tool combination reached a peak of 667 times at 80 m from the drill bit, indicating that the bottom movement of the well is complex, the collision frequency of the drill pipe part is relatively stable, and the curve fluctuation is not large.In the combination of single-centralizer drilling tools, the collision frequency at BHA was higher, reaching a peak of 775 collision frequencies at 105 m from the bottom of the well, and the collision frequency of the double-centralizer drilling tool combination was reduced by 14% compared with the single centralizer.From the perspective of the collision force of single and double centralizers (see Figure 25), regarding the collision force of the combination of double-centralizer drilling tools, four places exceed the force required for one failure, and at 155 m from the drill bit, the maximum collision force is 75,730.5N. The collision force of the whole well of the single-centralizer drilling tool combination did not exceed the force required for one failure, and the maximum collision force was 62,139 N at 60 m from the drill bit, which was reduced by 18% compared to the double-centralizer drilling tool.
tool combination reached a peak of 667 times at 80 m from the drill bit, indicating that the bottom movement of the well is complex, the collision frequency of the drill pipe part is relatively stable, and the curve fluctuation is not large.In the combination of single-centralizer drilling tools, the collision frequency at BHA was higher, reaching a peak of 775 collision frequencies at 105 m from the bottom of the well, and the collision frequency of the double-centralizer drilling tool combination was reduced by 14% compared with the single centralizer.In the comparison and analysis of the wellbore instability coefficient of the combination of single-and double-centralizer drilling tools (see Figure 26), it is found that nine positions of double-centralizer drilling tools exceed the wellbore instability coefficient of 1, but the combination of single-centralizer drilling tools does not exceed the wellbore instability coefficient of 1 in the whole well.It can be seen that the stability of the wellbore of the single-centralizer drilling tool combination is significantly higher than that of the double-centralizer drilling tool combination.the comparison and analysis of the wellbore instability coefficient of the combination of single-and double-centralizer drilling tools (see Figure 26), it is found that nine positions of double-centralizer drilling tools exceed the wellbore instability coefficient of 1, but the combination of single-centralizer drilling tools does not exceed the wellbore instability coefficient of 1 in the whole well.It can be seen that the stability of the wellbore of the single-centralizer drilling tool combination is significantly higher than that of the doublecentralizer drilling tool combination.Based on the above, it is recommended to install a single centralizer with packed holes 10.55 m away from the bit, which can effectively maintain wellbore stability.

Conclusions
In this paper, the mechanism of borehole instability caused by drill string vibration is studied, a dynamic simulation model of the drill string is established, and the concept of wellbore stability coefficient is proposed.The wellbore stability is visualized and analyzed by the numerical simulation research method, and the following conclusions are drawn: (1) When the single centralizer is placed close to the drill bit and the wellbore instability coefficient is low, the stability of the wellbore is improved compared with other situations; (2) When the packed hole centralizer is installed close to the drill bit, it will further reduce the borehole instability coefficient and improve the stability of the wellbore; (3) Compared to the double-centralizer drilling tool combination, the wellbore will be more stable with the single-centralizer drilling tool combination; (4) It is verified that this method can predict wellbore stability and optimize the drilling tool combination and drilling parameters according to the wellbore stability.The minimum force required to destroy a rock, N F n The normal force of the collision of the drill string with the wellbore, N F z The friction resistance in the Z-axis direction when the drill string collides with the wellbore, N k h Wellbore stiffness, N/m L Drill string length, m l The chord length of the arc corresponding to central angle α, m l 1 The length of the BHA, m M o The frictional moment of the drill string, N • m N Predict the number of collisions between drill string and wellbore per unit length N 1 The number of collisions per unit time per unit contact area N 2 The number of collisions between the string and the wellbore per unit length in the simulation time n The instability evaluation coefficient of the wellbore S The contact area when the drill string collides with the wellbore, m 2 S 1 Surface area of the inner wall of the wellbore per unit length, m 2 T The pure drilling time for drilling through the length of the BHA, s T 1 Simulation time, s v 1 The velocity before the drill string collides with the wellbore, m/s v 2 The velocity of the drill string after collision with the wellbore, m/s The rock strength of the wellbore, Pa

Figure 1 .
Figure 1.Contact collision force diagram of drill column section.


-Radial displacement of the center of mass of the drill string, m; -Damping factor; -Radial velocity of drill string center of mass, m/s; -Wellbore stiffness, N/m; -Wellbore diameter, m; -Drill string diameter, m.

Figure 1 .
Figure 1.Contact collision force diagram of drill column section.
c r -Recovery factor; µ r -Radial displacement of the center of mass of the drill string, m; c f -Damping factor; v r -Radial velocity of drill string center of mass, m/s; k h -Wellbore stiffness, N/m; d o -Wellbore diameter, m; d i -Drill string diameter, m.

Figure 2 .
Figure 2. Simplified schematic diagram of drilling tool combination.

Figure 2 .
Figure 2. Simplified schematic diagram of drilling tool combination.

Figure 4 .
Figure 4. Structure and physical diagram of vibration measuring device.

Figure 4 .
Figure 4. Structure and physical diagram of vibration measuring device.

Figure 4 .
Figure 4. Structure and physical diagram of vibration measuring device.

Figure 5 .Table 5 .
Figure 5. Schematic diagram of drilling tool combination and installation of vibration measurement equipment.Table 5. Drilling tool assembly dimensions.

Figure 5 .
Figure 5. Schematic diagram of drilling tool combination and installation of vibration measurement equipment.

Figure 8 .
Figure 8.Comparison of real drilling and simulated radial acceleration.

Figure 9 .
Figure 9.Comparison of real drilling and simulated tangential acceleration.

Figure 8 .
Figure 8.Comparison of real drilling and simulated radial acceleration.

Figure 8 .
Figure 8.Comparison of real drilling and simulated radial acceleration.

Figure 9 .
Figure 9.Comparison of real drilling and simulated tangential acceleration.

Figure 9 .
Figure 9.Comparison of real drilling and simulated tangential acceleration.

Figure 11 .
Figure 11.Schematic diagram of collision force under different rotational speed conditions.

Figure 11 .
Figure 11.Schematic diagram of collision force under different rotational speed conditions.

Figure 12 .
Figure 12.Schematic diagram of collision force.Figure 12. Schematic diagram of collision force.

Figure 12 .
Figure 12.Schematic diagram of collision force.Figure 12. Schematic diagram of collision force.

Figure 15 .
Figure 15.Borehole instability coefficient with drilling pressure of 16 t and speed of 50 rpm.

Figure 15 .
Figure 15.Borehole instability coefficient with drilling pressure of 16 t and speed of 50 rpm.

4. 2 .
Influence of Drill String Vibration on Wellbore Stability under Different Drilling Tool Combinations 4.2.1.Influence of Single Centralizer Placement on Wellbore Stability

Figure 16 .
Figure 16.Schematic diagram showing the placement positions of packed hole centralizers.

Figure 16 .
Figure 16.Schematic diagram showing the placement positions of packed hole centralizers.

Figure 19 .
Figure 19.Evaluation diagram of the wellbore instability coefficient.

4. 2 . 2 .
Effect of Single Centralizer Size on Wellbore Stability According to the above numerical simulation results, it is recommended to install the single centralizer at a distance of 10.55 m from the drill bit.The centralizers with packed hole sizes of 330.2 mm, 329.2 mm, and 327.2 mm are placed at a distance of 10.55 m from the drill bit, and the drilling pressure of 16 t and the speed of 50 rpm are used for simulation.

Figure 19 .
Figure 19.Evaluation diagram of the wellbore instability coefficient.4.2.2.Effect of Single Centralizer Size on Wellbore Stability According to the above numerical simulation results, it is recommended to install the single centralizer at a distance of 10.55 m from the drill bit.The centralizers with packed hole sizes of 330.2 mm, 329.2 mm, and 327.2 mm are placed at a distance of 10.55 m from the drill bit, and the drilling pressure of 16 t and the speed of 50 rpm are used for simulation.Observing the collision frequency diagram of the three centralizers (see Figure20), the collision frequency of the 327.2 mm centralizer in the BHA part is significantly higher than the 330.2 mm and 329.7 mm centralizers, and the collision frequency reaches a peak of 1592 times at 105 m from the drill bit.The collision frequency of the full-eye centralizer at BHA is relatively reduced for the 327.2 mm centralizer, reaching a peak of 775 times at 105 m from the drill bit.The collision frequency of the 329.2 mm centralizer is slightly smaller than that of other sizes on the whole, reaching a peak of 645 times at 160 m from the drill bit, and its collision frequency is reduced by 59% and 17% compared with the other two sizes.

Figure 20 .Figure 20 .
Figure 20.Collision frequency diagram of different centralizer sizes.According to the collision force diagram in Figure 21, it can be seen that the combination of drilling tools with centralizers with packed hole sizes of 330.2 mm, 329.2 mm, and 327.2 mm does not exceed the force required for a failure at BHA, and the combination of drilling tools with a 329.2 mm and 327.2 mm centralizer has a large collision force at the drill pipe and exceeds the force required for a single failure.The combination of drilling Figure 20.Collision frequency diagram of different centralizer sizes.

Figure 20 .
Figure 20.Collision frequency diagram of different centralizer sizes.

Figure 21 .
Figure 21.Schematic diagram of collision force of different centralizer sizes.Figure 21.Schematic diagram of collision force of different centralizer sizes.

Figure 21 .
Figure 21.Schematic diagram of collision force of different centralizer sizes.Figure 21.Schematic diagram of collision force of different centralizer sizes.

Figure 22 .
Figure 22.Schematic diagram of borehole instability coefficient of different centralizer sizes.4.2.3.Comparative Analysis of the Influence of Single and Double Centralizers on the Stability of the Wellbore

Figure 22 .
Figure 22.Schematic diagram of borehole instability coefficient of different centralizer sizes.

Figure 23 .
Figure 23.Schematic diagram of the combinations of single-and double-centralizer drilling tools.

Figure 23 .
Figure 23.Schematic diagram of the combinations of single-and double-centralizer drilling tools.

Figure 24 .
Figure 24.Schematic diagram of collision frequency of single and double centralizers.Figure 24.Schematic diagram of collision frequency of single and double centralizers.

Figure 24 .
Figure 24.Schematic diagram of collision frequency of single and double centralizers.Figure 24.Schematic diagram of collision frequency of single and double centralizers.

From
the perspective of the collision force of single and double centralizers (see Figure25), regarding the collision force of the combination of double-centralizer drilling tools, four places exceed the force required for one failure, and at 155 m from the drill bit, the maximum collision force is 75,730.5N. The collision force of the whole well of the single-centralizer drilling tool combination did not exceed the force required for one failure, and the maximum collision force was 62,139 N at 60 m from the drill bit, which was reduced by 18% compared to the double-centralizer drilling tool.

Figure 25 .
Figure 25.Schematic diagram of collision force of single and double centralizers.

Figure 25 .
Figure 25.Schematic diagram of collision force of single and double centralizers.

Figure 26 .
Figure 26.Schematic diagram of borehole instability coefficient of single and double centralizers.

NomenclatureA
The number of units per contact area c f Damping coefficient of drill string c r Recovery factor of drill string d i Drill string diameter, m d o Wellbore diameter, m F The maximum collision force per meter, N F N V h Mechanical drilling speed, m/s v r Radial velocity of drill string center of mass, m/s v s Stick-slip speed of the drill string, m/s α Surrounding angle of drill string ∆ The width of the ring gap, m δ Clay layer thickness, m µ r Radial displacement of the center of mass of the drill string, m µ(v s ) Coefficient of friction of the drill string [σ]

Table 2 .
Constraint relationships between components.

Table 2 .
Constraint relationships between components.

Table 4 .
Contact force parameters.

Table 6 .
Vibration measuring device measures sensor installation parameters.

Table 6 .
Vibration measuring device measures sensor installation parameters.

Table 6 .
Vibration measuring device measures sensor installation parameters.

Table 8 .
Evaluation table of borehole instability.

Table 8 .
Evaluation table of borehole instability.