1. Introduction
Pneumatic down-the-hole (DTH) hammer drilling technology has been widely used in the fields of thermal storage mining and geological exploration [
1,
2]. Its principle is that the piston inside the DTH hammer is mainly driven by compressed air to carry out high-speed reciprocating motion to break the rock [
3], which has the advantages of high drilling efficiency and small pollution to the formation compared with traditional rotary drilling [
4,
5]. In the process of pneumatic DTH hammer drilling, the near-bit measurement while drilling (MWD) is crucial. It is a device that is used to monitor various important drilling parameters relating to the near drill bit in real time. By accurately measuring data, the real-time drilling status can be monitored, and the drilling efficiency can be optimized. However, the vibration at the drill bit during pneumatic DTH hammer drilling is as high as 40 g, which poses a serious challenge to the performance of the near-bit MWD. The existing MWD tools are mainly designed for oil drilling and geological drilling. Its vibration level is only 20 g. Adapting to the high-intensity and low-frequency vibrations in pneumatic DTH hammer drilling is difficult. Therefore, for the special vibration environment of pneumatic DTH hammer drilling, it is important to research and develop a new type of anti-vibration MWD to improve drilling efficiency and ensure operational safety.
Since the vibration response of a circuit unit is caused by the vibration of the drill string near the drill bit, studying the vibration characteristics of the drill string is crucial for establishing a vibration model. Many scholars have carried out research on the vibration mechanics and vibration control of drill strings. For example, Rey O F [
6] used the finite difference method to analyze the effect of damping, mass, and other factors on vibration by simplifying the drilling system. Aarrestand [
7] and Wei et al. [
8] focused on the vibration dissipation process of a drill string by using the linear elasticity theory and experimental methods. They pointed out that damping plays an important role in the vibration of a drill string. Zhang et al. [
9,
10] investigated the effect of the damper’s installation position on the longitudinal vibration of a drill string by establishing a mathematical model and using a computer simulation. They determined the optimal installation position of the damper. Hakimi H [
11] calculated the intrinsic frequencies of a drill string’s longitudinal, transverse, and torsional vibration by using the differential quadrature method (DQM). Bu [
12] investigated the longitudinal forced vibration of a drill string under the action of cyclic impacts of the drilling cycle of a DTH hammer in hard rock. Dong [
13] explored vibration and impact-suppression techniques for drill strings. They proposed three control methods: passive control, active control, and semi-active control. Further, they comparatively analyzed the characteristics of the various methods, equipment, and application scopes. Tian [
14] analyzed the effects of stiffness and damping of the damper and the installation position on the vibrations of downhole drilling tools. Riane et al. [
15] eliminated the severe stick–slip vibrations that appear along the drill string of a rotary drilling system according to the LQG observer-based controller approach. The above studies revealed the vibration mechanisms of a drill pipe and quantified its vibration characteristics. However, there is still a gap in the research related to the forced vibration response of the internal unit of the near-bit MWD during DTH hammer drilling and its damping technology.
In other fields, the research and development of vibration-damping methods is more mature. Du [
16] proposed the principle of elastic vibration isolation and the frequency design criterion. She established a dynamic model of an engine’s elastic system and determined the dynamic characteristic parameters. Wang [
17] addressed the problem of the vibration resistance of MEMS IMUs in high-g shock environments. He embedded gyroscopes and accelerometers into a symmetric hexahedral frame and utilized a viscoelastic-damping damper to achieve internal vibration damping. Zhou et al. [
18,
19,
20,
21] addressed the problem of damper–damping matching in a vehicle’s steel plate spring suspension system. They established mathematical models of the suspension system for the optimal damping ratio, and the root mean square values of dynamic deflection and vibration velocity. The optimal damping ratio of the damper in the suspension system was determined. Li [
22] proposed four new passive vibration control damping forms based on inertia–spring–damping (ISD) structures. Through the establishment of dynamic equations and the acquisition of a steady-state amplitude amplification factor, the optimal damping ratio and stiffness ratio of the system were obtained. Shatskyi and Velychkovych [
23] presented a new dry friction shell shock absorber design. Such shell shock absorbers are projected to be used in the mining, oil and gas industries. The development of the above damping methods is more mature, and they have good damping effects in related fields. However, these methods are mostly applied to vibration suppression in precision instruments and have not yet been utilized in the more severe environment of low-frequency, high-intensity impact vibrations.
Accordingly, this study proposes a vibration-damping system, which is based on non-dominated sorting genetic algorithm (NSGA II) parameter optimization. This method is designed to solve the problems relating to the significant longitudinal vibrations in pneumatic DTH hammer drilling processes. First, the forced vibration model of a short-section circuit unit under complex drilling conditions was constructed and analyzed. Then, the vibration-damping model of the circuit unit was established via model simplification, the separation of variables, the superposition of vibration shapes, etc. For the problem relating to the complexity of the model parameter solution, the forced vibration full response function of the circuit unit was used for analysis. This study proposes two key parameters, the logarithmic attenuation rate of the amplitude and the displacement amplification factor, with both affecting the performance of the vibration-damping system. NSGA II is used to find the optimal stiffness-damping parameters of the system. The optimal stiffness and damping parameters of the system are derived using NSGA II(Surrogates and Evolutionary Algorithms Laboratory, Indian Institute of Technology, Surat, India). Finally, a vibration-damping scheme adapted to the drilling conditions of a DTH hammer was designed. It effectively suppress the vibration amplitude of the circuit unit in a violent vibration environment. This enhances the efficiency and safety of the DTH hammer drilling operations.
2. Overall Design for the Near-Bit MWD
The MWD is installed near the drill bit for the real-time measurement of downhole pressure, temperature, and other parameters. The main structure includes the measuring sub-shell, sensor, circuit, battery, and connecting structure between the upper and lower drill pipes. Among them, the sensors are mounted on the measuring sub-shell and the circuit unit, which are used to measure the real-time downhole data; the circuit and battery are mounted on the circuit unit, which are used for signal acquisition, signal processing, and power supply. The vibration-damping device support frame is mounted between the vibration-damping device and the circuit unit, which is used to provide space for the sensor to be mounted. The vibration-damping devices at both ends are mounted tightly on the inside of the measuring sub-shell by the reducer union. The reducer union connects the MWD to the upper and lower drill pipes. The pressure detection hole is opened on the wall of the shell, which is used to measure the pressure of the annulus. According to the actual project background, the preliminary design of the MWD is 127 mm in diameter and 750 mm in length. The internal installation structure of the MWD is shown in
Figure 1.
In pneumatic DTH hammer drilling operation processes, a drill string will produce very complex coupled vibrations, of which longitudinal vibrations are the most intense. Although the focus of this study is on longitudinal vibrations, it is worth noting that transverse and torsional vibrations in the drilling string system can similarly affect the performance of a system. Although these forms of vibrations were not explored in depth in this study, their presence cannot be ignored. Future research efforts will be directed towards the development of more comprehensive vibration-damping solutions to cover all types of vibration. This will further improve the reliability and measurement accuracy of the system.
To protect the circuit unit, this study adopts spring-damper elements at both ends of the circuit unit. These elements are flexibly connected to the upper and lower reducer unions. The elastic deformation of the spring element and the energy dissipation of the damping element can isolate the vibration. However, the structure and parameters of the stiffness damper directly affect the vibration-damping performance. Therefore, it is necessary to optimize the parameters of the damper.
3. The Vibration Model of Circuit Unit
To solve and design the structure and parameters of the stiffness damper, and to investigate the damping performance of the spring-damping system, it is necessary to study the vibration response characteristics of the circuit unit.
3.1. Vibration Response Model of Circuit Unit
A DTH hammer piston-bit-rock-drill string system is a complex and highly nonlinear coupled vibration system. During rock-drilling processes, there is not only the impact vibration generated by the impactor but also the reverse force generated by the rock on the drill pipe. These forces jointly provoke the complex vibration response of the drill string system, which, in turn, will directly act on the near-bit MWD [
24,
25]. To facilitate analysis and research, the following assumptions are made [
26,
27]:
(i) The drill string is in the state of linear elastic deformation, within the range of elastic deformation; (ii) the effect of shear force on the deformation of the drill string is not considered; (iii) the effect of temperature is omitted; (iv) both ends of the drill string are free, and only the horizontal displacement is limited; (v) transverse and torsional vibrations of the drill string are not considered, and only the longitudinal vibration characteristics of the drill string are investigated; and (vi) the vibration-damper is regarded as a stiff and damped damper, and the mass of the damper is on the circuit unit.
For the drill string, the impact force generated by the drill bit is an external excitation, and the drill string system is a forced vibration. For the circuit unit, the vibration response of the drill string provides excitation, and the circuit unit is forced vibration. Therefore, it is necessary to establish a dynamic model of the drill string system, including the DTH hammer impactor, drill string, circuit unit impactor, etc. The analytical model is shown in
Figure 2.
3.2. Dynamic Equation of Circuit Unit
Since only the longitudinal vibration characteristics of the circuit unit are investigated, the system can be regarded as a damped single-degree-of-freedom system under external excitation [
28]. Based on the mechanical structure of the circuit unit mounting described above, the circuit unit is simplified to a mass, and the mass is represented by
m. The vibration-damping unit is simplified as a set of spring-damper systems, where stiffness and damping are
k and
c.
Due to the influence of gravity, the static equilibrium position of the damping system is not the center point. The static equilibrium position is taken as the origin of the coordinate system to establish a coordinate system. The downward direction is specified as the positive direction. Given the role of damping in energy consumption, vibration analysis is often used in the principle of energy equivalence. Moreover, the circuit unit quality does not change with movement. In this system, the elastic force of the system can be simplified to linear, and the damping can be simplified to linear viscous damping. As the deformations of the two groups of springs are consistent, they can be regarded as a set of springs in parallel with equivalent stiffness:
. Damping is the same, with equivalent damping
. Therefore, the model is further simplified, as shown in
Figure 3.
The force analysis of the circuit unit is carried out, and the mass
m is taken as the separator. According to Newton’s second law and D’Alembert principle, the vibration equation can be obtained.
This is expressed as follows:
The above equation is the forced vibration response of the circuit unit under external excitation, where m is the mass of the circuit unit in kg, k is the stiffness of the damping mechanism, c is the damping coefficient of the damping mechanism, is the acceleration of the measuring mechanism in m/s2, is the velocity of the measuring mechanism in m/s, is the speed of the drill string in m/s, is the displacement of the measuring mechanism in m, and is the displacement of the drill string in m.
4. Solving the Vibration Response of Circuit Unit
4.1. Decomposition of the Vibration Response
The forced total vibration response under simple harmonic force is similar to the superposition of a general solution of the damped free vibration response and a special solution of simple harmonic vibration. Equation (2) is a second-order linear non-homogeneous ordinary differential equation, whose solution should be a superposition of the general solution of the homogeneous equation and a special solution of the non-homogeneous equation. The general solution of the homogeneous equation is the solution of the damped free vibration response of the system under no external excitation. Over time, the amplitude of the part of the general solution gradually decays. After a while, it decays to zero. This part is the transient response of the circuit unit. The special solution part of the response amplitude does not change with time; the frequency and excitation frequency are consistent, whereas the amplitude and phase depend on the excitation amplitude and system parameters, independent of the initial conditions. This part is the steady-state response of the circuit unit. The forced vibration of the whole system is the superposition of transient vibration and steady-state vibration. The full response of the vibration system is the superposition of the steady-state response and transient response of the circuit unit.
4.1.1. Transient Response
The homogeneous part of
is written as follows:
The damping coefficient and stiffness are known to be
,
, which are obtained by substituting it into Equation (3):
Let the general solution be
, and substitute it into Equation (4) to obtain the general solution of the equation.
When , it is an over-damped case in which the circuit unit will gradually return to the equilibrium position after only one time at most, with no periodic change in the factor and no oscillation characteristics. In this case, the energy input to the system by the initial excitation is quickly consumed by the damping. The system is too late to generate reciprocating vibration and is not conducive to the damping of continuous excitation.
When
, it is an under-damped case. The transient response of the circuit unit is expanded and collated, using Euler’s formula to obtain:
Denoted by
where
is the vibration amplitude of the transient response and
is the phase of the transient response.
This is an under-damped condition in which the circuit unit vibrates back and forth near the system’s equilibrium position. Its amplitude decreases with the exponential law
. The logarithmic decay rate of the amplitude is introduced here to describe the speed of amplitude decay.
Equation (8) is transformed into an expression of the stiffness and damping coefficients. Denoted as
:
For this system, the larger the damping rate, the faster the amplitude decay of the transient response and the shorter the time to reach the steady state. However, too large a damping rate can also result in the excitation of the system by the drill string being applied directly to the circuit unit. This will cause the circuit unit to be subjected to excessive transient forces. In engineering applications, it is generally taken that
. That is
4.1.2. Steady-State Response
The periodic excitation
can be decomposed into the sum of infinitely many harmonic components using a Fourier expansion. That is
where
.
For each harmonic component, a non-homogeneous special solution is obtained. According to the principle of the superposition of linear systems, the steady-state characteristics of the system can be obtained. The second-order non-homogeneous differential equation of arbitrary harmonics is written as follows:
Let the special solution of the above equation be
, and substituting it into the differential equation:
The above equation can hold for any time, so the coefficients of the terms
and
on both sides of the equal sign must be equal. That is
The solution is
,
, where
is the amplification factor:
Equation (14) is transformed into an expression of the stiffness and damping coefficients. Denoted as
:
indicates the magnification of the static displacement amplitude under dynamic vibration. The smaller the amplification factor, the better the damping effect. Furthermore, only when the amplification factor is greater than 0 but less than 1 will the amplitude of the vibration will be attenuated. Therefore, in the parameter design of the vibration-suppression system for the MWD, the amplification factor should be ensured to be small enough. That is
At the same time, the vibration isolation effect of the damping system is achieved when the frequency ratio is greater than
. However, to ensure the stability and reliability of the system, it is usually no greater than 5. That is
4.2. Solving for Optimal Stiffness and Damping Coefficient of the Damper
Through the decomposition of the vibration response of the above circuit unit, two main parameters, the logarithmic attenuation rate of the amplitude and the displacement amplification factor, are obtained. They are both functions of damping and stiffness, and the two parameters will affect each other when seeking their optimum values. Therefore, the two functions are used as the objective functions, which are transformed into a multi-objective optimization problem. The optimal stiffness and damping parameters of the vibration system can be obtained by seeking the optimal solution through a suitable algorithm.
The multi-objective optimization problem is more than a single-objective optimization problem, multiple objective functions may conflict. Therefore, the result is an optimal solution set. In this study, NSGA-II was chosen to solve the optimal stiffness damping parameters of the vibration system, which has the advantages of fast convergence, high computational efficiency, and good stability compared with other multi-objective optimization algorithms.
Based on the previous section, Equations (9) and (16) are used as the objective functions of this multi-objective optimization problem, and Equations (10), (17) and (18) are used as the boundary conditions of this function. This is sorted out as follows:
A non-dominated sorting genetic algorithm was developed based on the above objective function and boundary conditions. To optimize programming and calculation, take the reciprocal of the displacement amplification factor and find the maximum of both objective functions simultaneously using the algorithm.
An initial set of populations is first generated randomly when the algorithm is executed. Each spring-damping parameter represents a set of solutions and ensures that these initial values satisfy the boundary conditions in Equation (19). Then, a hierarchy is performed based on the dominance of the spring-damped parameter over the set of multi-objective functions. This can ensure that each individual parameter solved maintains the property of being a non-inferior solution on at least one objective. Next, a crowding calculation is performed for each spring-damped parameter to quantify the intensity of the solution. Thus, this promotes the diversity and equilibrium of the solution set. After that, the spring-damped parameter with a higher grade is preferred. According to the congestion index, the spring-damped parameter with a sparse distribution is selected in each layer. Then, it is added to the population of the subsequent evolution. Then, new populations are obtained through the use of two genetic operations, crossover, and mutation, to increase the diversity of solutions within the population. After the newly generated individuals are merged with the existing population through non-dominated sorting and crowding evaluation, the individuals with higher priority are again filtered to form the new generation of the population. This is carried out until the maximum number of iterations is reached. Eventually, the algorithm will converge to a set of high-quality optimal solution sets, resulting in the optimal stiffness and damping coefficients that are applicable to the model. The optimization process is shown in
Figure 4.
The known circuit unit prototype is 8 kg, and the DTH hammer (China University of Geosciences, Beijing, China) used in the mining of geothermal resources has an operating frequency range of 5–15 Hz. Let the impact frequency be 10 Hz, and then substitute it into Equation (19). A set of optimal stiffness and damping coefficients are obtained using MATLAB R2023a (MathWorks, Natick, MA, USA), which are rounded to , , respectively. The optimal displacement amplification factor and amplitude attenuation rate of and are obtained, respectively. The parameter that measures the damping effect of the system is the absolute motion conductivity——which damps up to 90.78%.