CFD-DEM Simulation of Reverse Circulation Pneumatic Cuttings Removal during Coal Seam Drilling

To solve the problems that the borehole depth is shallow and the borehole formation rate is low during the gas drainage drilling in soft coal seam with current cuttings removal method, a new technology of reverse circulation pneumatic cuttings removal is put forward. (e CFD-DEM coupling method is used to establish the simulation model of cuttings-air two-phase flow in drill pipe. (e effects of the air velocity for cuttings removal and the mass flow rate of cuttings on the flow characteristics, cuttings removal effect and pressure drop of cuttings-gas two-phase flow are analysed. (e results show that the drag force of drilling cuttings becomes larger with the increase of air velocity and the stratified flow characteristic is obvious.(e drill cuttings migration ratio is positively correlated with the air velocity for cuttings removal and negatively correlated with the mass flow rate of cuttings. When the mass flow rate of cuttings is constant, the increase of air velocity for cuttings removal leads to the increase of pressure drop in the inner hole of drill pipe. When the air velocity of cuttings removal is constant, the mass flow rate of cuttings and the pressure drop in the inner hole of drill pipe increases. (erefore, the appropriate air velocity should be selected considering the energy consumption during cuttings removal.


Introduction
Gas extraction is the main means to prevent and control gas outburst in coal seam. And, drill cuttings conveying during the borehole drilling in coal seam is the key to ensure the drilling depth and the drilling efficiency. e mechanical removal and the hydraulic removal of cuttings have a great impact on the stability of borehole wall, which is easy to cause the borehole collapse and blockage, sticking of the drill pipe, and other accidents [1]. e current cuttings removal method with pressure air removes the cuttings through the annular hole between drill pipe and borehole wall by the pressure air through the inner hole of drill pipe. It is difficult to collect the drill cuttings at the outlet of boreholes, which is easy to cause the dust pollution [2][3][4]. So, a new reverse circulation pneumatic cuttings removal technology is proposed that the annular hole between drill pipe and borehole wall is used as the air inlet and the inner hole of drill pipe is used as the cuttings removal channel [5,6].
Numerical methods based on computational fluid dynamics (CFD) have mainly been used. Niu and Zhang [7] studied the critical wind speed under different cuttings removal rate and solid gas ratio as well as the pressure loss of air flow in drill pipe inner hole and annular hole. It is found that the critical wind speed and pressure loss are positively correlated with the solid gas ratio. Guo et al. [8] performed three dimensional and time-dependent calculations by using the finite volume CFD of sudden pipe expansion flows. Riaz and Sadiq [9] set up the governing equations with the help of similarity transformations and handled the solution of boundary value problems by perturbation procedure. So, the analytical solutions for fluid and particulate phase velocities, mean flow rates, and pressure gradient profile have been presented, while a numerical treatment has been carried out for pressure rise. And, the exact solutions of non-Newtonian multiphase fluid through peristaltic pumping characteristics in an annulus having complaint walls and applied magnetic field was studied [10]. Abdelsalam et al. [11] and Ijaz et al. [12] studied the two-phase flows characteristics in the biologically inspired pumping systems.
e Computational Fluid Dynamics and Discrete Element Method Coupling Algorithm (CFD-DEM) is a well-developed method that has been proven effective in the study of liquidsolid and gas-solid two-phase flow system. e bilateral coupled CFD-DEM method was used to numerically simulate the coal particles flow pattern in the swirling pneumatic conveying process and the interaction between coal particles and pipe wall to reveal the mechanism of coal particle swirling pneumatic conveying [13,14]. Based on the CFD-DEM coupling algorithm, the dust cyclone separator in negative pressure fixedpoint sampling is simulated to study the effect of air velocity and coal mass flow rate on the coal-air two-phase separation process [15]. Akhshik et al. [16,17] applied the CFD-DEM coupling method to study the effects of drill pipe rotation and drill cuttings shape on the cuttings flow characteristics in different inclined angle drilling holes. And, the numerical simulation is consistent with the experimental results. Based on the CFD-DEM coupling method, Shao et al. [18] simulated the migration of large-diameter unconventional drill cuttings particles to study the drill cuttings removal rule in borehole and reveal the influence of cuttings shape and content on the removal effect of cuttings. Sun et al. [19] used the CFD-DEM coupling method to simulate the trajectory of drill cuttings. e simulation results of the cuttings concentration coincided with the experimental data by Kim et al. [20], which proved the accuracy and feasibility of the CFD-DEM coupling method in the drill cuttings removal. erefore, the CFD-DEM coupling method is also used to simulate the cuttings-gas two-phase flow in the inner hole of drill pipe during the horizontal drilling borehole in coal seam. e research can reveal the effects of the air velocity for cuttings removal and the mass flow rate of cuttings on the flow characteristics, cuttings removal effect, and pressure drop of cuttings-gas two-phase flow to improve the drilling efficiency in coal seam. Figure 1, the reverse circulation pneumatic cuttings removal system mainly includes drilling rig, drill pipe, drilling bit, cyclone separator, and Roots vacuum pump. e Roots vacuum pump provides the air force for the pneumatic conveying of drill cuttings. Under the suction effect of vacuum pump, air enters the bottom of the drill hole along the annular hole between drill pipe and borehole wall. e drill cuttings generated by the operating drilling bit are wrapped and carried into the inner hole of drill pipe through the discharge hole. And, the inner hole of drill pipe is used as the cuttings removal channel. e cyclone separator acts as the dust removal device to separate the drill cuttings from the cuttings-air two-phase flow. e drill cuttings are collected in the separator under the gravity while the filtered air is discharged through the Roots vacuum pump.

Introduction of the Reverse Circulation Pneumatic Cuttings Removal System. As shown in
e system has high removal efficiency of drill cuttings and plays a role of dust control and dust suppression by collecting the discharged cuttings.

Drill Cuttings-Air Two-Phase Flow Control Equation.
e drill cuttings generated by the drilling bit breaking coal are smaller in size. Considering the continuity and fluidity of a large number of drill cuttings along the inner hole of drill pipe are similar to the movement of fluid, the dual-Eulerian model is used to describe the two-phase flow of drill cuttings and air. e air phase control equation includes the mass conservation and momentum conservation equations [16,17,21]: where, ε a is the porosity, ρ a is the air density, v a is the air velocity, ∇ is the divergence operator, p is the pressure on a fluid mesh element, η is the aerodynamic viscosity, and S is the momentum sink which is the ratio of the sum of the forces between drill cuttings and air in the fluid mesh unit to the volume of the fluid mesh unit: where N is the number of drill cuttings in the fluid mesh unit, F c−a,i is the interaction force between drill cuttings and air, and V cell is the fluid grid element volume.
where C D is the resistance coefficient, d c is the drill cuttings diameter, and v c is the drilling speed. e drill cuttings motion follows Newton's second law. So, where m i is the quality of drill cuttings, v i is the drilling speed, F c is the contact force of drill cuttings, I i is the rotary inertia of drill cuttings, ω i is the angular velocity of drill cuttings, and k j�1 T ij is the resultant moment of drill cuttings: where, r i is the vector from the center of mass of cutting i to the contact point and T DT,i is the moment produced by the slip-rotation: where ρ is the cutting density, d i is the mean diameter of ith cutting, C DR is the rotational drag coefficient, and Ω is the relative angular velocity of the cutting to the air.

Drill Cuttings Contact Model.
During the movement of drill cuttings along the inner hole of drill pipe, the Hertz-Mindilin contact model is used to describe the cuttings-cuttings contact and the contact between the cuttings and the inner wall of drill pipe. e contact forces of drill cuttings are shown in Figure 2. e contact forces between the cuttings or between the cuttings and the inner wall of drill pipe are divided into normal force F c n , normal damping force F d n , tangential force F c t , and tangential damping force F d t . e normal force F c n is the function of normal overlap δ n , equivalent Young's modulus Y * , and equivalent radius R * [13]: 1 where υ i , Y i , and R i and υ j , Y j , and R j are Poisson's ratio, Young's modulus, and radius of cuttings i and cuttings j, respectively. e normal damping force F d n is c � ln e �������� ln 2 e + π 2 , where S n is the normal stiffness, m * is the equivalent mass, v rel n is the normal component of relative velocity between drill cuttings, e is the recovery coefficient, And m i andm j are the quality of drill cuttings i and drill cuttings j, respectively. e tangential force F c t depends on the tangential overlap δ t and the tangential stiffness S t : where G * is the equivalent shear modulus and G i and G j are the shear modulus of cuttings i and cuttings j, respectively. e tangential damping force F d t is where v rel t is the normal component of relative velocity between drill cuttings.
When the tangential force between the drill cuttings is greater than the static friction force, the slip between the drill cuttings will occur. e tangential force is where φ s is the sliding friction coefficient.

Simulation Model.
According to the working principle of the reverse circulation pneumatic cuttings removal system, the geometric model including borehole, drilling bit, and drill pipe is built, as shown in Figure 3. According to the working principle of the reverse circulation pneumatic cuttings removal system, the diameter of drilling bit is 94 mm and the diameter of discharge hole in drilling bit is 40 mm. e outer diameter of drill pipe is 73 mm, and the inner hole diameter of drill pipe is 40 mm. e drilling bit is connected with the drill pipe. e diameter of the borehole in coal seam is about 100 mm. In order to analyze the flow characteristics of the drill cuttings along the inner hole of drill pipe, the length of drill pipe is set as 4.9 m and the drilling bit length is set as 0.1 m. So, the geometry model length of the drill pipe-bit reverse circulation pneumatic cuttings removal simulation is 5 m. e hexahedral mesh is used to mesh the geometric model. A total of 53209 mesh elements and 47428 mesh nodes are generated. e mesh generation is shown in Figure 4. e annular hole between drill pipe and borehole wall is set as the inlet and the end of the inner hole of drill pipe is set as the outlet.
To solve the air phase, the unsteady pressure solver with a gravitational acceleration of 9.81 m/s 2 along the negative Y direction is set. e turbulence model is the realizable k-ε turbulence model. e inlet is set as speed inlet with the air velocity of 5.13 m/s, 6.85 m/s、8.56 m/s、10.27 m/s, and 11.98 m/s, respectively. Correspondingly, the air velocity in the inner hole of drill pipe is 15 m/s, 20 m/s, 25 m/s, 30 m/s, and 35 m/s, respectively. e turbulence intensity is 5%, and the turbulent viscosity ratio is 0.5. e outlet is set as pressure outlet with the pressure of −3000 Pa. e finite volume method and the QUICK algorithm are used to discretize the air phase control equation and the momentum equation.
e second-order upwind scheme is used to discretize the turbulent kinetic energy and the turbulent diffusion equation. e SIMPLE algorithm is used for the phase-to-phase coupling of pressure-velocity phase. e simulation parameters are set as Table 1. e drill cuttings generated by the drilling bit in coal seam were collected, and the particle size analysis was carried out. It was found that the particle size of the cuttings was mainly less than 1 mm. erefore, the spherical particles with the particle diameter of 1 mm were used to simulate the drill cuttings. e drill cuttings are generated at the entrance of the discharge hole of drilling bit. And, the cuttings production is set as 0.06 kg/s, 0.08 kg/s, 0.10 kg/s, 0.12 kg/s, and 0.14 kg/s, respectively. In order to avoid the drill cuttings accumulating at the discharge hole when the cuttings are generated, the initial velocity of cuttings is set as 2 m/s in the X direction. e drill cuttings phase solution time step is set to be 2 × 10 −6 s, and the air phase time step is set to be 1 × 10 −4 s.

Model Verification.
To validate the simulation predictions, the reverse circulation pneumatic cuttings removal test device is built which is shown in Figure 5.
e Roots vacuum pump typed of QZSR125A is used as the power source of the experimental system. e boost range is −50∼0 kPa, and the maximum flow is 6.38 m 3 /min. During the test, the flow of Roots vacuum pump can be changed by adjusting the speed of variable frequency motor of Roots vacuum pump to adjust the air velocity. A feeding device is used to simulate the generation of drill cuttings. e feeding device is composed of bunker and screw conveyor. e mass flow rate of cuttings is controlled by adjusting the speed of screw conveyor. e transparent rigid plexiglass tube is used to simulate the drill pipe and borehole. e two-phase separation of drill cuttings and air is completed in the cyclone separator. In order to obtain the pressure drop of drill pipe under the steady flow of cuttings and air, the pressure measuring point 1 is set at the distance of 3 m from the drill pipe inlet, and the pressure measuring point 2 is set at the downstream 2 m from the pressure measuring point 1. Figure 6 is the actual experimental equipment.
In the experiment, the mass flow rate of cuttings was selected as 0.06 kg/s, 0.10 kg/s, and 0.14 kg/s, and the air velocity was selected as 15 m/s, 20 m/s, 25 m/s, 30 m/s, and 35 m/s. e pressure drop per unit length between measuring point 1 and measuring point 2 was tested and is shown in Figure 7. From the comparison of pressure drop per unit length, it shows that the results of CFD-DEM simulation reasonably agree with those of the experiment.

Drilling Cuttings Flow Pattern in Inner Hole of Drill
Pipe.
e cuttings flow patterns of 0∼0.5 m and 4.5∼5 m segments in the inner hole of drill pipe under different air velocities are shown in Figure 8 when the cuttings mass flow rate is 0.14 kg/s. From Figure 8, it can be seen that the drilling cuttings are evenly generated at the entrance of the discharge hole of drilling bit. en, the cuttings enter the inner hole of drill pipe at the initial speed of 2 m/s. Because the cuttings speed is less than the air velocity, the drilling cuttings accelerate and form obvious cuttings flow under the drag force of air flow, and then under the gravity action, the drilling cuttings flow gathers at the bottom of the inner hole of drill pipe and moves along the direction of the air flow until they are discharged from the inner hole of drill pipe.

Mathematical Problems in Engineering
When the air velocity is 15 m/s, the cuttings flow pattern is shown in Figure 8(a). Because of the low air velocity, the characteristic length of the cuttings flow is relatively short and deposits at the bottom of the inner hole of drill pipe. e typical characteristics of pipe bottom flow are presented along the flow direction and a thicker cuttings layer is formed. e moving speed of cuttings in the 4.5∼5 m segment is obviously lower than that in the 0∼0.5 m segment, which indicates that the cuttings have decelerated during the movement. When the air velocity is 20 m/s, the cuttings flow pattern is shown in Figure 8 e stratified flow characteristic is obvious. e stratified flow can be roughly divided into two layers: most of the cuttings move along the bottom of the inner hole to form the cuttings conveying layer, and some cuttings suspension movement forms the mixed conveying layer. Figure 8(e) shows that the amount of cuttings floating in the flow direction increases and the stratified flow characteristics is obvious when the air velocity is 35 m/s.

Cuttings Distribution along Flow Direction.
e distribution of drill cuttings along the flow direction can reflect the accumulation of drill cuttings in the inner hole of drill pipe.
e drill pipe-bit model is equally divided into 25 segments along the flow direction, and the amount of drill cuttings at different positions along the flow direction is counted at the time of 2 s, respectively. And, the ratio of the amount of cuttings in different positions to the total number of cuttings in the inner hole of drill pipe is defined as the amount proportion of cuttings along the flow direction which is used to characterize the distribution of cuttings along the flow direction, as shown in Figure 9.  erefore, the amount proportion of cuttings along the flow direction decreases gradually and the maximum proportion increases with the increase of the air velocity.

Cuttings Velocity Variation along Flow Direction.
e variation of drill cuttings velocity along the flow direction can reflect the acceleration and motion characteristics of drill cuttings along the flow direction. e average velocity of drill cuttings at different positions along the flow direction is calculated at time of 2 s separately. e variation of drill cuttings velocity along the flow direction is shown in Figure 10.  As can be seen from Figure 10, when the air velocity is 15 m/s, the drill cuttings enter the inner hole of drill pipe at the initial speed of 2 m/s and accelerate under the drag force of air flow. Drill cuttings decelerate along the flow direction after short-distance acceleration movement. With the increase of the mass flow rate of drill cuttings, the acceleration distance of drill cuttings becomes shorter, the maximum velocity decreases, and the deceleration phenomenon is obvious. When the mass flow rate of cuttings is 0.12 kg/s and 0.14 kg/s, the cuttings decelerate along the flow direction first and accelerate at the end of the drill pipe. As the mass flow rate of drill cuttings increases, the accumulated cuttings increase, the cuttings layer becomes thicker, the area of air flow passage decreases and the air flow velocity increases at this position, which makes the cuttings located at the upper part of the cuttings layer accelerate again. When the air velocity are 20 m/s, 25 m/s, 30 m/s, and 35 m/s, respectively, the cuttings velocity increases along the flow direction. Under the same air velocity, the cuttings velocity at the same position decreases with the increase of the mass flow rate of cuttings. e drag force on the cuttings is proportional to the velocity difference between the cuttings and air flow. When the cuttings enter the inner hole of drill pipe at a lower initial velocity, the drag force on the cuttings is larger because of the large velocity difference between the cuttings and the air flow. erefore, the cuttings move in the inner hole of drill pipe with a higher acceleration. With the increase of cuttings velocity, the velocity difference between cuttings and air flow decreases, the drag force decreases, and the acceleration of the cuttings decreases. So, the velocity curve of drill cuttings in the 1∼5 m segment is relatively flat compared with that in the 0∼1 m segment. Figures 9 and 10 show the distribution and velocity of the drill cuttings at different positions along the flow direction under different working conditions. In summary, when the air velocity is 15 m/s, the velocity of drill cuttings along the flow direction decreases and accumulates. e cuttings cannot effectively discharge from the inner hole of drill pipe in time and even block the inner hole of drill pipe, which makes the cuttings removal difficult and affects the drilling depth and drilling rate of borehole. When the air velocity is greater than 20 m/s, the velocity of drill cuttings increases along the flow direction, and no accumulation phenomenon occurs, so the cuttings removal can be achieved smoothly.

Mathematical Problems in Engineering
When the air velocity is constant, the cuttings mass flow rate increases from 0.06 kg/s to 0.14 kg/s, the amount of drill cuttings increases, and the individual cuttings energy obtained from the air flow decreases. erefore, the cuttings velocity decreases along the flow direction with the increase of the mass flow rate of cuttings. Similarly, when the mass flow rate of drill cuttings is constant, the individual cuttings energy obtained from the air flow increases with the increase of air velocity. e drill cuttings accelerate faster and move in the inner hole of drill pipe at a higher speed. erefore, in the application of the new technology of reverse circulation pneumatic cuttings removal, whether the drill cuttings can be discharged from the drilling hole directly depends on the air velocity and the mass flow rate of drill cuttings. In order to ensure that the cuttings are smoothly discharged from the drilling hole, the air velocity cannot be lower than 20 m/s.

Flow Characteristics Analysis of Air
Phase. In order to analyze the influence of operating parameters on the air flow characteristics, the axial air flow velocity distribution at the position of 4.75 m under different cuttings mass flow rate and different air velocity is compared and analysed, as shown in Figure 11.

Forces between Cuttings and Air under Different Working
Conditions. In the process of reverse circulation pneumatic cuttings removal, there is a relationship between the action force and the reaction force besides the momentum exchange between drill cuttings and air. e interaction forces between drill cuttings and air in 0∼2 s under different working conditions are extracted. As shown in Figure 12, the variation of interaction forces between drill cuttings and air under different working conditions with time is analysed. Figure 12(a) shows that the interaction forces between drill cuttings and air increase continuously with time when the air velocity is 15 m/s. When the mass flow rate of cuttings is 0.06 kg/s, the interaction forces between two phases increase slowly with time. With the increase of the mass flow rate of cuttings, the slope of the force curves between two phases increase and the forces between two phases increase. When the air velocity is 20 m/s, 25 m/s, 30 m/s, and 35 m/s, respectively, the interaction forces between the cuttings and air increase first. And then, the forces stabilize with time. e time-varying curves of the interaction forces between two phases are divided into the "growth section" and the "stable section" for analysis. It can be seen that, under the same air velocity, with the increase of the mass flow rate of cuttings, the slope of the "growth section" increases, the corresponding time length increases, and the interaction forces between the two phases of the "stable section" increase. When the mass flow rate of cuttings is constant, the slope of the "growth section" curve increases with the increases of the air velocity, and the corresponding time length becomes shorter, while the interaction forces between the two phases in the "stable section" become larger. When the air velocity is greater than 20 m/s, there is no accumulation of drill cuttings in the inner hole of drill pipe. e amount of drill cuttings in the drill pipe reaches the dynamic balance state and the interaction forces between the two phases remain relatively stable. When the mass flow rate of cuttings keeps constant, with the increase of the air velocity, the movement speed of drill cuttings becomes faster, and the time of reaching the "stable section" of the interaction forces between two phases becomes shorter. Although the increase of air velocity decreases the number of cuttings in the drill pipe, the drag force of a single cuttings by the flow field increases significantly at this time. So, the interaction forces between the two phases tend to increase. At the same air velocity, with the increase of the mass flow rate of cuttings, the amount of cuttings in the drill pipe increases and the interaction forces between the two phases increase. At the same time, the increase of the number of drill cuttings slows down the cuttings velocity. So, the time when the interaction forces between two phases reach the "stable section" becomes longer. It can be seen that the variation of the interaction forces between drill cuttings and air with time can be used to judge whether the cuttings-air two-phase flow reaches a stable state in the simulation process.

Cuttings Removal Effect under Different Working
Conditions. In order to analyze the influence of different working conditions on the effect of reverse circulation pneumatic cuttings removal, the cuttings migration ratio is used to evaluate the cuttings removal effect [22]. e cuttings migration ratio is defined as the ratio of the average cuttings velocity in the drill pipe to the air velocity under the corresponding working conditions. Figure 13 shows the cuttings migration ratio under different working conditions. Figure 13 shows that, when the mass flow rate of cuttings is constant, the cuttings migration ratio is positively correlated with the air velocity, while when the air velocity is constant, the cuttings migration ratio is negatively correlated with the mass flow rate of cuttings. In addition, when the air velocity is low, the difference between the cuttings migration ratios under different mass flow rate of cuttings is large. With the increase of the air velocity, the difference between the cuttings migration ratios under different mass flow rate of cuttings decreases gradually. Drill cuttings move along the flow direction with the characteristics of pipe bottom flow under the action of airflow drag force and self-gravity after entering the inner hole of drill pipe. When the air velocity is low, the pipe bottom flow characteristics are significant. At this time, the collision and friction between cuttings-cuttings, drill pipecuttings have a greater influence on the cuttings movement, which result in the average velocity of drill cuttings far lower than the air velocity. e effect of cuttings removal is poor. With the increase of air velocity, the characteristics of pipe bottom flow weaken, and some drill cuttings move in suspension in the inner hole of drill pipe. e collision and friction between cuttings-cuttings, drill pipe-cuttings decrease. So, the cuttings migration ratio increases. When the air velocity is constant, the increase of the mass flow rate of cuttings increases the amount of cuttings in the drill pipe. e influence of collision and friction between cuttings- cuttings, drill pipe-cuttings on the cuttings movement is intensified. So, the cuttings removal effect becomes worse. Compared with the low air velocity, the higher air velocity can make the drill cuttings be discharged from the drill pipe in time and effectively. e amount of drill cuttings in the drill pipe is relatively small. e increase of the mass flow rate of cuttings has little effect on the cuttings movement. erefore, under the condition of low air velocity, the cuttings migration ratio varies greatly with different mass flow rate of cuttings. Under the condition of high air velocity, the increase of mass flow rate of cuttings has little effect on the cuttings migration ratio.

Pressure Drop of Cuttings-Air Two-Phase Flow.
e drill cuttings generated by the drilling bit breaking coal enter the inner hole of drill pipe under the action of air flow and discharge along the flow direction. In this process, all kinds of energy consumed by air flow and drill cuttings movement are compensated by the pressure energy of the air flow. So, the pressure drop of air flow can be used to characterize the energy consumption of cuttings-air twophase flow [23]. e pressure drop of cuttings-air twophase flow in the inner hole of drill pipe under different working conditions is shown in Figure 14.
As can be seen from Figure 14, the pressure drop of cuttings-air two-phase flow is significantly higher than that in the inner hole of drill pipe under the pure air flow, which indicates that the drill cuttings cause the greater energy consumption to the air flow. e pressure drop shows the same trend with the increase of air velocity under different mass flow rate of cuttings. At the same mass flow rate of cuttings, the pressure drop increases with the increase of air velocity. And, at the same air velocity, the pressure drop also increases with the increase of the mass flow rate of cuttings. e pressure drop in the inner hole of drill pipe is mainly composed of two parts. One part is the friction pressure drop of air flow. e other part is the pressure drop caused by the drill cuttings which includes the energy consumption caused by the collision and friction among cuttings-cuttings and drill pipe-cuttings during the cuttings removal and the kinetic energy obtained by the drill cuttings. When the mass flow rate of cuttings is constant, the mixing ratio of drill cuttings to air decreases with the increase of the air velocity, but the friction pressure drop of air flow and the kinetic energy of drill cuttings increases. So, the pressure drop increases with the increase of the air velocity. At the same air velocity, with the increase of the mass flow rate of cuttings, the mixing ratio of cuttings and air increases and the porosity in the inner hole of drill pipe decreases. At this time, the energy consumption increases due to the collision and friction among cuttings-cuttings and cuttings-drill pipe, so the pressure drop also increases with the increase of the mass flow rate of cuttings.
In conclusion, the higher air velocity increases the pressure drop of air flow and the air consumption. erefore, in order to avoid unnecessary energy consumption, on the premise of ensuring the smooth discharge of drill cuttings from the drill hole, the appropriate air velocity should be selected to determine the air volume of Roots vacuum pump. Taking the simulation model in this paper as an example, when the air velocity is 15 m/s, the amount of cuttings increases along the flow direction, and the velocity of cuttings decreases. So, the drill cuttings will accumulate in the inner hole of drill pipe. If the air velocity is greater than 20 m/s, the amount of cuttings will decrease along the flow direction. e speed and pressure drop of drill cuttings will increase. So, 20 m/s is selected as the appropriate air velocity of cuttings removal.

Conclusions
In this paper, the CFD-DEM coupling method is used to simulate the two-phase flow of drill cuttings and air in the inner hole of drill pipe with the reverse circulation pneumatic cuttings removal technology during the horizontal drilling in coal seam. e effects of the air velocity and the mass flow rate of cuttings on the twophase flow characteristics, cuttings removal effect, and pressure drop of the flow field are analysed. e results show that the following: (1) With the increase of air velocity, the drag force of drilling cuttings becomes larger. e characteristic length of drilling cuttings flow becomes longer. e thickness of the cuttings layer decreases further along the flow direction. A small amount of drilling cuttings distributes in the upper part of the inner hole of drill pipe and the cuttings movement speed is higher than that of the bottom cuttings. e stratified flow characteristic is obvious. (2) When the air velocity of cuttings removal is 15 m/s, the proportion of drill cuttings along the flow direction increases, the forces between the cuttings and air increase continuously with time, and the cuttings velocity decreases. So, the drill cuttings are easy to accumulate in the inner hole of drill pipe. If the air velocity is greater than 20 m/s, the proportion of cuttings along the flow direction decreases and the cuttings velocity increases. (3) When the mass flow rate of cuttings is constant, the cuttings migration ratio is positively correlated with the air velocity. When the air velocity is constant, the cuttings migration ratio is negatively correlated with the mass flow rate of cuttings. e interaction forces between drill cuttings and air increase first. And then, the forces stabilize with time. (4) At the same mass flow rate of cuttings, the pressure drop of the cuttings-air two-phase flow increases with the increase of air velocity. At the same air velocity, the pressure drop of the cuttings-air twophase flow increases with the increase of the mass flow rate of cuttings. So, 20 m/s is selected as the appropriate air velocity for cuttings removal considering the pressure drop of cuttings-air twophase flow and the flow characteristics of drill cuttings.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.