Study on the Impact Pressure of Swirling-Round Supercritical CO2 Jet Flow and Its Influencing Factors

Supercritical carbon dioxide (SC-CO2) jet is capable of decreasing the threshold pressure of rock breakage and mitigating formation damage, owing to its low viscosity, high diffusivity, and extremely-low surface tension. The swirling-round jet holds the advantages of both a swirling jet and a round jet. Therefore, the comprehensive technique, swirling-round SC-CO2 (SR-SC-CO2) jet, is expected to substantially enhance rock-breaking efficiency. However, theoretical analysis of the flow field characteristics of SR-SC-CO2 has not been reported yet. This work aims to lay a theoretical foundation for employing SR-SC-CO2 in drilling and fracturing. The flow field is simulated using Naiver-Stokes equations and the RNG k-ε turbulence model. Sensitivity analysis, regarding pressure drop of the nozzle, confining pressure, fluid temperature, jetting distance, the diameter of the nozzle’s central hole, and grooving area, are performed. We show that the combined swirling-round SCCO2 jet flow could maintain a relatively larger axial as well as tangential velocity compared to a single approach of swirling jet or round jet, enabling one to acquire a deeper oillet and expand the perforation area effectively. The simulation results substantiate the enormous potential of SR-SC-CO2 in improving rock-breaking efficiency and clarify the influence of relevant parameters on the impact pressure of the jet flow.


Introduction
With the development of petroleum exploration and exploitation toward a deeper formation, the quantity of deep and ultra-deep wells is ever-increasing [1][2][3]. One critical issue of successful deep well drilling is to enhance the rate of penetration (ROP) and lower drilling costs [4]. High-pressure water (HPW) jet has played a crucial role in deep-well exploitation. However, several problems of HPW, including high water consumption, clay swelling, groundwater contamination, treatment of flow-back fluid, and high threshold pressure of rock breaking, remain to be solved [5][6][7][8][9][10][11][12]. By improving the jet pressure, modifying nozzle structures, or introducing new fluid, a series of novel jet technologies have been proposed.
Li and Shen [13] put forward a dual-jet concept that combines a swirling jet and a round jet. Buckman et al. [14] designed the swirling-round jet bit to effectively eliminate the central low-speed zone and the convex plate of the swirling jet. The swirling-round jet is a mixed jet flow of the swirling jet passing through the swirling groove of the impeller and the round jet flowing through the impeller's central channel. The rock-breaking effect and flow field characteristics of the swirling-round water (SRW) jet have been effect and flow field characteristics of the swirling-round water (SRW) jet have been extensively investigated [15,16]. It has been concluded that SRW synthesizes the greater rock-breaking depth of the round jet and larger rock-breaking area of the swirling jet, enabling one to decrease the threshold pressure of rock breakage and sharply increase rockbreaking efficiency.
Supercritical CO2 has been recognized as a successful candidate for drilling and fracturing fluid, given its low viscosity, high diffusivity, and extremely-low surface tension [17][18][19][20][21][22][23][24]. Those unique physicochemical properties provide significant advantages to the SC-CO2 jet. Firstly, the threshold pressure of rock breaking of the SC-CO2 jet is far below that of the water jet [21,22]. Secondly, considering the stronger adsorption capability of CO2 onto the rock surface [25][26][27][28], CO2 is expected to replace initially adsorbed shale-gas and coal bed methane in a matrix, enhancing hydrocarbon recovery [29][30][31][32][33][34][35]. Furthermore, clay swelling and formation damage caused by water injection are mitigated to the maximum extent [17,18]. Experimental results suggested that comparing with the water jet, SC-CO2 jet is more efficient in rock-breaking due to its higher rock-breaking speed and lower threshold pressure. Kolle [17] reported that the ROP of SC-CO2 jet in the Mancos shale is 3.3 times as high as the water jet, while the specific energy of rock breakage is only 20% of the water jet. Kolle and Marvin [36] substantiated the feasibility of using SC-CO2 as drilling fluid by conducting field trials of SC-CO2 assisted drilling. Tian et al. [16] investigated the rock-breaking effect of the swirling SC-CO2 jet. He et al. [37] numerically studied the flow field of the SC-CO2 abrasive jet and clarified the influence of confining pressure, fluid temperature, and jetting distance on rock breakage. Cheng et al. [38] concluded that the SC-CO2 jet exhibits higher flow speed, longer core zone, and stronger capability of enhancing pressurizing effect if compared with the water jet.
Tian et al. [39] experimentally studied the rock-breaking efficiency of swirling-round SC-CO2 jet. Their results suggested that the radial rock-breaking capability of swirling jet can be perfectly combined with the axial rock-breaking capacity of round jet. Figure 1 shows the rock-breaking effects of the swirling jet, the round jet, and their combination. The stronger oillet-enlarging ability and greater perforation depth of the swirling-round jet can be observed. The swirling-round SC-CO2 jet allows us to find a balance between perforating depth and area, and thus the rock-breaking efficiency is enhanced. In this paper, we present numerical simulation on the rock-breaking efficiency and impact pressure of the swirling-round SC-CO2 jet. Navier-Stokes equations and the RNG k-ε turbulence model are used to simulate the three-dimensional flow field.

Numerical Simulation
In this section, we demonstrate the geometric model and meshing (Section 2.1), governing equations (Section 2.2), and structure parameter setting (Section 2.3). In this paper, we present numerical simulation on the rock-breaking efficiency and impact pressure of the swirling-round SC-CO 2 jet. Navier-Stokes equations and the RNG k-ε turbulence model are used to simulate the three-dimensional flow field.

Numerical Simulation
In this section, we demonstrate the geometric model and meshing (Section 2.1), governing equations (Section 2.2), and structure parameter setting (Section 2.3). nozzle are 4 mm 2 , and the central hole's diameter of a spiral impeller is 2 mm. A cylinder of 30 mm in diameter and 12 mm in height is designed as the impactor at the nozzle outlet. The established model incorporates a swirling-round jet (including an inlet section, rotary section, mixing section, convergent section, outlet section, and lengthened section) and impact zone of jet flow. When a low-speed fluid enters the nozzle, a part of fluid flows through the spiral impeller, forming a swirling jet, while the other part of fluid flows through the central hole, generating a round-straight jet. The swirling and round-straight jet flows are thoroughly blended in the mixing chamber, giving rise to a mixed flow of a low speed. The mixed flow passes through the nozzle outlet and the lengthened section, then the created high-speed swirling-round jet flow impacts the downstream surface and eventually flows out the field. The pressure at the nozzle inlet is set as the inlet boundary value of pressure, and the ambient pressure determines the outlet-boundary pressure at the nozzle outlet. The central channel, convergent section, outlet section, and lengthened section of the nozzle and the impact zone of dual-jet flow are partitioned with a structured grid. In contrast, an unstructured grid is used for the nozzle entrance, helical channel, and mixing section. Figure 2b shows the amplification of meshing at the rotating impeller and central channel.

Geometric Model and Meshing
Energies 2020, 13, x FOR PEER REVIEW 3 of 1 Figure 2a demonstrates the three-dimensional geometric model for the simulation o the swirling-round SC-CO2 jet flow field. The cross-section areas of the swirling-round je nozzle are 4 mm 2 , and the central hole's diameter of a spiral impeller is 2 mm. A cylinde of 30 mm in diameter and 12 mm in height is designed as the impactor at the nozzle outle The established model incorporates a swirling-round jet (including an inlet section, rotar section, mixing section, convergent section, outlet section, and lengthened section) an impact zone of jet flow. When a low-speed fluid enters the nozzle, a part of fluid flow through the spiral impeller, forming a swirling jet, while the other part of fluid flow through the central hole, generating a round-straight jet. The swirling and round-straigh jet flows are thoroughly blended in the mixing chamber, giving rise to a mixed flow of low speed. The mixed flow passes through the nozzle outlet and the lengthened section then the created high-speed swirling-round jet flow impacts the downstream surface an eventually flows out the field. The pressure at the nozzle inlet is set as the inlet boundar value of pressure, and the ambient pressure determines the outlet-boundary pressure a the nozzle outlet. The central channel, convergent section, outlet section, and lengthene section of the nozzle and the impact zone of dual-jet flow are partitioned with a structure grid. In contrast, an unstructured grid is used for the nozzle entrance, helical channel, an mixing section. Figure 2b shows the amplification of meshing at the rotating impeller an central channel.

Governing Equations
The swirling-round SC-CO2 jet flow obeys the conservation laws of mass, momen tum, and energy. The RNG k-ε model is used for the turbulence closure, given that the SC CO2 jet is a turbulent flow. The complete governing equations are given as follows.

Continuity Equation
The equation for continuity equation, or conservation of mass, can be described b [40]: (1 where ρ is density, kg/m 3 ; ui denotes the ith component of time-averaged velocity, m/s; is time, s; xi represents the ith coordinates, m. i = 1, 2, 3.

Momentum Conservation Equation
Conservation of momentum in each direction is written as [40]:

Governing Equations
The swirling-round SC-CO 2 jet flow obeys the conservation laws of mass, momentum, and energy. The RNG k-ε model is used for the turbulence closure, given that the SC-CO 2 jet is a turbulent flow. The complete governing equations are given as follows.

Continuity Equation
The equation for continuity equation, or conservation of mass, can be described by [40]: ∂ρ ∂t where ρ is density, kg/m 3 ; u i denotes the ith component of time-averaged velocity, m/s; t is time, s; x i represents the ith coordinates, m. i = 1, 2, 3.

Momentum Conservation Equation
Conservation of momentum in each direction is written as [40]: where µ is molecular viscosity, Pa·s; p is static pressure, Pa; σ ij is the Kronecker delta function; Reynolds stresses, −ρu i u j representing the effect of turbulence, is given by: where k is the turbulent kinetic energy per unit mass, J/kg; µ t is turbulent viscosity, m 2 /s, calculated by µ t = ρc µ k 2 /ε; ε is dissipation rate, m 2 /s 3 ; c µ is a constant with a default value of 0.09.

Energy Conservation Equation
Conservation of energy is given as [41,42]: where c p is specific heat, J/kg·K; T is temperature, K; E is the total energy, J; k h is thermal conductivity, W/m·K; S h denotes the heat of chemical reaction or any other volumetric heat sources; Pr t is the turbulent Prandtl number for energy; (τ ij ) eff is the deviatoric stress tensor, defined as: The effective viscosity µ eff is given by:

RNG k-ε Equations
The RNG k-ε equations are given by [42,43]: ∂(ρε) ∂t where α k and α ε are the inverse effective Prandtl numbers for k and ε, respectively, with both given by 1.39; C 3ε is the buoyancy coefficient, equal to one if the flow direction is perpendicular to the direction of gravity and zero if the flow direction is parallel to the direction of gravity; C 1ε and C 2ε are constants, given by 1.42 and 1.68, respectively. G k is the mean velocity gradient induced turbulence kinetic energy, determined by: G b is the buoyancy induced turbulence kinetic energy, determined by: where g i is the component of the gravitational vector in the ith direction, m/s 2 ; β is the coefficient of thermal expansion, K −1 ; Y M denotes the contribution of fluctuating dilatation to the overall dissipation rate, determined by Energies 2021, 14, 106 where M t is the turbulent Mach number, defined as where a is the sound speed, m/s. R ε is given by: where η = Sk/ε, η 0 = 4.38, β = 0.012.

Boundary Conditions Inlet Boundary
The nozzle inlet serves as the inlet of the geometrical model with the pressure inlet condition (Figure 2a). The inlet pressure and fluid temperature are given by 30 MPa and 400 K, respectively.

Outlet Boundary
The outlet of the geometrical model with the pressure outlet condition is defined at the jet nozzle exit (Figure 2a). The outlet pressure and fluid temperature are given by 40 MPa and 400 K, respectively.

Wall Boundary Conditions (i) Wall boundary of u
The law-of-the-wall for mean velocity yields [42]: where κ is the von Karman constant, set as 0.4187; E m is an empirical constant, given by 9.793; y * is the dimensionless distance from the wall, defined by , and u * is dimensionless velocity, denoted by k 1/2 p c 1/4 µ u p τ w /ρ . k p is turbulence kinetic energy at the nearwall node P, J/kg; y p is the distance from point P to the wall, m; u p is the mean velocity of the fluid at node P, m/s; c µ is a constant with a default value of 0.09; τ w is the wall-shear stress, Pa, determined by: y * varies from 11.225 to 300 to guarantee the validity of the logarithmic distribution of velocity. When y * is less than 11.225 at the wall-adjacent cells, the laminar stress-strain relationship is employed, i.e., u * = y * .
(ii) wall boundary of T The wall boundary of T is determined by: where Pr is the Prandtl number; u c is the average velocity magnitude at y * = y * T , m/s; y * T represents the non-dimensional thermal sublayer thickness. q is the density of heat flux at the wall, W/m 2 ; P e is given by Jayatilleke (1969) for a smooth wall: Energies 2021, 14, 106 6 of 15 (iii) wall boundary of k and ε In the k-ε model, the k equation is solved in the whole domain incorporating the wall-adjacent cells. The boundary condition for k imposed on the wall is: where n is the local coordinate normal to the wall, m. The wall boundary of the generation term of turbulent kinetic energy G k and dissipation rate ε are determined from the local equilibrium assumption:

Structure Parameter Setting
Structure parameters of the nozzle used in the three-dimensional flow field simulation of the swirling-round SC-CO 2 jet are summarized in Table 1.  [44]. The velocity field corresponding to the initial guess of the pressure field is obtained by solving the momentum equation. The corrected value of pressure is acquired from the continuity equation, resulting in an updated pressure and velocity field. The updated pressure field is set as a new guess for the pressure field if the calculated results are not convergent, and the calculation process is repeated until convergence.

Results and Discussion
In this section, we discuss the obtained simulation results. Section 3.1 presents the flow field characteristics of the swirling-round SC-CO 2 jet. Section 3.2 performs sensitivity analysis of critical parameters, including the pressure drop of the nozzle, confining pressure, jet flow temperature, jetting distance, the diameter of the nozzle's central hole, and the grooving area of the nozzle.

Flow Field Simulation
The supercritical CO 2 enters the nozzle and is split into two jets. One of the jets accesses the straight channel, forming a round jet, while the other one enters the rotating channel, generating a swirling jet. The round jet and swirling jet are fully mixed in the mixing section, giving rise to a swirling-round SC-CO 2 jet flow, which hits the wall after flowing through the nozzle outlet and lengthened section. Figures 3 and 4 reveal the flow Energies 2021, 14, 106 7 of 15 field simulation results. When the SC-CO 2 flow through the convergent section of the nozzle, fluid pressure decreases, and flow velocity sharply rises, resulting in the swirlinground SC-CO 2 jet flow that hits the wall on the right side. Not surprisingly, the hitting speed reached its maximum on the jet axis, while getting lower to both sides.
1 Figure 3. The velocity contour of swirling-round SC-CO 2 jet flow field.
The axial and tangential velocities of the mixed jet are investigated to unveil the jet energy's attenuation law. Figure 5 exhibits the distribution and evolution of axial velocity. The axial radius r is given as the abscissa. Sensitivity analysis is performed in regard to the dimensionless jetting distance, defined as the ratio of jetting distance L and the nozzle's outlet diameter d. It is shown that the flow speed reaches its maximum at the central line, demonstrating more concentrated energy if L/d is zero. In contrast, energy drops sharply in the radial direction. This leads to a more substantial jet impact and a certain depth in perforating. When ejecting fluid of high speed, the entrainment effect of a swirling jet favors driving the flow of surrounding fluid and exchanging energy, which attenuates the axial velocity of jet flow. Figure 6 shows the distribution and evolution of tangential velocity. Given a particular jetting distance, the tangential velocity first increases and then decreases with an increasing jetting radius, reaching a maximum at a specific position away from the central line. Meanwhile, the action range of tangential velocity rises with the jetting distance, and the maximum value of tangential velocity declines with the jetting distance.
According to the distribution profiles shown in Figure 6, the tangential velocity at the central line reaches a local minimum far less than the values of the two maxima. However, in comparison, the axial velocity reaches its peak at the central line ( Figure 5). Compared to a single approach, the combined swirling-round SC-CO 2 jet flow could maintain a relatively larger axial as well as tangential velocity, enabling one to obtain a deeper oillet and expand the perforation area effectively, thus acquiring satisfactory rock breakage. The axial and tangential velocities of the mixed jet are investigated to unveil the je energy's attenuation law. Figure 5 exhibits the distribution and evolution of axial velocity. The axial radius r i given as the abscissa. Sensitivity analysis is performed in regard to the dimensionless jet ting distance, defined as the ratio of jetting distance L and the nozzle's outlet diameter d It is shown that the flow speed reaches its maximum at the central line, demonstrating more concentrated energy if L/d is zero. In contrast, energy drops sharply in the radia direction. This leads to a more substantial jet impact and a certain depth in perforating When ejecting fluid of high speed, the entrainment effect of a swirling jet favors driving the flow of surrounding fluid and exchanging energy, which attenuates the axial velocity of jet flow. Figure 6 shows the distribution and evolution of tangential velocity. Given a partic ular jetting distance, the tangential velocity first increases and then decreases with an in creasing jetting radius, reaching a maximum at a specific position away from the centra line. Meanwhile, the action range of tangential velocity rises with the jetting distance, and the maximum value of tangential velocity declines with the jetting distance.     According to the distribution profiles shown in Figures 6, the tangential velocity a the central line reaches a local minimum far less than the values of the two maxima. How ever, in comparison, the axial velocity reaches its peak at the central line ( Figure 5). Com pared to a single approach, the combined swirling-round SC-CO2 jet flow could maintai a relatively larger axial as well as tangential velocity, enabling one to obtain a deeper oille and expand the perforation area effectively, thus acquiring satisfactory rock breakage.  Figure 7 shows the axial velocity profiles of swirling-round SC-CO 2 jet and water jet to compare their flow fields under the same parameters. Temperature and confining pressure are set as 450 K and 30 MPa, respectively. The wall surface is located at −20 mm of the axial coordinate, and 0 mm at the axial coordinate represents the nozzle's outlet. As depicted in Figure 7, the maximum speeds of swirling-round SC-CO 2 jet flow exhibit higher values compared with that of water jet under all pressure drops. For example, at the pressure drop of 30 MPa, the maximum speed of swirling-round SC-CO 2 jet reaches 358 m/s, significantly higher than 263 m/s of the water jet. According to the Bernoulli equation, if the nozzle's pressure drop maintains constant, then the SC-CO 2 jet speed is higher than the water jet speed, due to the lower density of SC-CO 2 .
Energies 2020, 13, x FOR PEER REVIEW 9 of 1 Figure 7 shows the axial velocity profiles of swirling-round SC-CO2 jet and water je to compare their flow fields under the same parameters. Temperature and confining pres sure are set as 450 K and 30 MPa, respectively. The wall surface is located at −20 mm o the axial coordinate, and 0 mm at the axial coordinate represents the nozzle's outlet. A depicted in Figure 7, the maximum speeds of swirling-round SC-CO2 jet flow exhibi higher values compared with that of water jet under all pressure drops. For example, a the pressure drop of 30 MPa, the maximum speed of swirling-round SC-CO2 jet reache 358 m/s, significantly higher than 263 m/s of the water jet. According to the Bernoull equation, if the nozzle's pressure drop maintains constant, then the SC-CO2 jet speed i higher than the water jet speed, due to the lower density of SC-CO2.

Pressure Drop of the Nozzle
The nozzle's pressure drop plays a significant role in determining the kinetic energy of jet flow, which directly influences the characteristics of the jet flow field. As the nozzle' pressure drop increases, jet energy rises, improving the effect of rock breakage. Figure 8 depicts impact pressure under variant pressure drops of the nozzle. It is suggested that a higher pressure drop yields higher impact pressure. As a result, if the confining pressure remains invariant, the impact pressure increases with an increasing pressure drop. The pressure drop induced kinetic energy of jet flow rises as the nozzle's pressure drop in creases, and thus the impact pressure acting on the wall surface arisen from the jet flow

Pressure Drop of the Nozzle
The nozzle's pressure drop plays a significant role in determining the kinetic energy of jet flow, which directly influences the characteristics of the jet flow field. As the nozzle's pressure drop increases, jet energy rises, improving the effect of rock breakage. Figure 8 depicts impact pressure under variant pressure drops of the nozzle. It is suggested that a higher pressure drop yields higher impact pressure. As a result, if the confining pressure remains invariant, the impact pressure increases with an increasing pressure drop. The pressure drop induced kinetic energy of jet flow rises as the nozzle's pressure drop increases, and thus the impact pressure acting on the wall surface arisen from the jet flow. However, it should be noted that the impact range does not show a noticeable enlargement at a higher pressure drop.
The nozzle's pressure drop plays a significant role in determining the kinetic energy of jet flow, which directly influences the characteristics of the jet flow field. As the nozzle's pressure drop increases, jet energy rises, improving the effect of rock breakage. Figure 8 depicts impact pressure under variant pressure drops of the nozzle. It is suggested that a higher pressure drop yields higher impact pressure. As a result, if the confining pressure remains invariant, the impact pressure increases with an increasing pressure drop. The pressure drop induced kinetic energy of jet flow rises as the nozzle's pressure drop increases, and thus the impact pressure acting on the wall surface arisen from the jet flow. However, it should be noted that the impact range does not show a noticeable enlargement at a higher pressure drop.

Confining Pressure
The jet flow field can be affected by the downhole confining pressure. As shown in Figure 9, with increasing confining pressure, the pressure level of the wall exhibits an increase of the same magnitude. The strength and action range of impact pressure, however, is not influenced. Although the confining pressure allows determining the pressure level of the whole flow field, it shows no direct influence on jet kinetic energy and, thus, on the strength and action range of impact pressure.

Confining Pressure
The jet flow field can be affected by the downhole confining pressure. As shown in Figure 9, with increasing confining pressure, the pressure level of the wall exhibits an increase of the same magnitude. The strength and action range of impact pressure, however, is not influenced. Although the confining pressure allows determining the pressure level of the whole flow field, it shows no direct influence on jet kinetic energy and, thus, on the strength and action range of impact pressure.  Figure 10 demonstrates fluid temperature effect on impact pressure. It is anticipated that temperature might influence the flow field by affecting the physical properties of a fluid. However, as shown in the figure, impact pressure and the action range of impact pressure are hardly influenced by temperature changes varying from 430 K to 550 K. It is worth noting that further variation of temperature demonstrates similar results (data not shown). Therefore, we conclude that fluid temperature can be ignored for engineering applications. noting that further variation of temperature demonstrates similar results (data not shown). Therefore, we conclude that fluid temperature can be ignored for engineering applications. Figure 9. Influence of the confining pressure on impact pressure. Figure 10 demonstrates fluid temperature effect on impact pressure. It is anticipated that temperature might influence the flow field by affecting the physical properties of a fluid. However, as shown in the figure, impact pressure and the action range of impact pressure are hardly influenced by temperature changes varying from 430 K to 550 K. It is worth noting that further variation of temperature demonstrates similar results (data not shown). Therefore, we conclude that fluid temperature can be ignored for engineering applications. Figure 10. Influence of fluid temperature on impact pressure. Figure 10. Influence of fluid temperature on impact pressure.

Jetting Distance
Jetting distance has a significant influence on rock breakage. Proper adjustment is essential to use jetting energy efficiently. As shown in Figure 11, impact pressure decreases with an increasing jetting distance. The impact range, however, expands slightly. When the jetting distance is small, the development of jet flow is inadequate and, therefore, the impact range is narrowed. Nevertheless, if the jetting distance is oversized, the high dissipation of jet energy might occur and lead to a slump of impact pressure. Accordingly, only if the jetting distance is appropriately selected can we achieve perfect impact pressure and action range of the SC-CO 2 jet flow, thus acquiring the desired rock-breaking effect.

Jetting Distance
Jetting distance has a significant influence on rock breakage. Proper adjustment i essential to use jetting energy efficiently. As shown in Figure 11, impact pressure de creases with an increasing jetting distance. The impact range, however, expands slightly When the jetting distance is small, the development of jet flow is inadequate and, therefor the impact range is narrowed. Nevertheless, if the jetting distance is oversized, the hig dissipation of jet energy might occur and lead to a slump of impact pressure. Accordingly only if the jetting distance is appropriately selected can we achieve perfect impact pres sure and action range of the SC-CO2 jet flow, thus acquiring the desired rock-breakin effect. Figure 11. Influence of jetting distance on impact pressure. Figure 12 shows the influence of the diameter of the nozzle's central hole on impac pressure. 1 mm, 1.5 mm, 2 mm, and 2.5 mm are selected as typical values of the centra hole's diameter. It is shown that impact pressure and its action range increase with in creasing hole size. If the nozzle's pressure drop remains unchanged, larger hole diamete gives rise to the higher kinetic energy of round jet and thus the axial velocity and impac pressure acting on the wall. Therefore, increasing the nozzle's pressure drop and using Figure 11. Influence of jetting distance on impact pressure. Figure 12 shows the influence of the diameter of the nozzle's central hole on impact pressure. 1 mm, 1.5 mm, 2 mm, and 2.5 mm are selected as typical values of the central hole's diameter. It is shown that impact pressure and its action range increase with increasing hole size. If the nozzle's pressure drop remains unchanged, larger hole diameter gives rise to the higher kinetic energy of round jet and thus the axial velocity and impact pressure acting on the wall. Therefore, increasing the nozzle's pressure drop and using a nozzle of larger central hole size favor enhancing the impact pressure and action range of swirling-round SC-CO 2 and improving the efficiency of rock breakage. Figure 11. Influence of jetting distance on impact pressure. 3.2.5. Diameter of the Nozzle's Central Hole Figure 12 shows the influence of the diameter of the nozzle's central hole on impact pressure. 1 mm, 1.5 mm, 2 mm, and 2.5 mm are selected as typical values of the central hole's diameter. It is shown that impact pressure and its action range increase with increasing hole size. If the nozzle's pressure drop remains unchanged, larger hole diameter gives rise to the higher kinetic energy of round jet and thus the axial velocity and impact pressure acting on the wall. Therefore, increasing the nozzle's pressure drop and using a nozzle of larger central hole size favor enhancing the impact pressure and action range of swirling-round SC-CO2 and improving the efficiency of rock breakage.

Grooving Area of the Nozzle
One may expect that the flow channel's area can have a considerable effect on the generation of swirling jet. Figure 13 presents the influence of the grooving area (3 mm 2 , 4 mm 2 , 5 mm 2 , and 6 mm 2 ) on swirling-round jet flow. As depicted, the larger grooving area of the impeller yields a slightly higher impact pressure. This could be explained by the fact that, rather than swirling jet, the round jet is the primary control of impact pressure. As a result, the grooving area could be neglected in practice due to its weak influence on the strength and impact range of the jet flow.

Grooving Area of the Nozzle
One may expect that the flow channel's area can have a considerable effect on th generation of swirling jet. Figure 13 presents the influence of the grooving area (3 mm 2 , mm 2 , 5 mm 2 , and 6 mm 2 ) on swirling-round jet flow. As depicted, the larger grooving are of the impeller yields a slightly higher impact pressure. This could be explained by th fact that, rather than swirling jet, the round jet is the primary control of impact pressure As a result, the grooving area could be neglected in practice due to its weak influence on the strength and impact range of the jet flow.

Conclusions
We present a synthetic simulation on the impact pressure of swirling-round super critical CO2 jet flow. Sensitivity analysis of relevant parameters are performed. Our re search allows drawing the following conclusions: The swirling-round SC-CO2 jet exhibits a higher flow speed and a more substantia Figure 13. Influence of grooving area on impact pressure.

Conclusions
We present a synthetic simulation on the impact pressure of swirling-round supercritical CO 2 jet flow. Sensitivity analysis of relevant parameters are performed. Our research allows drawing the following conclusions: The swirling-round SC-CO 2 jet exhibits a higher flow speed and a more substantial rock-breaking effect compared with the water jet flow, showing combined advantages of the round jet and swirling jet.
Compared to a single approach of swirling jet or round jet, the combined swirlinground SC-CO 2 jet flow could maintain a relatively larger axial as well as tangential velocity, enabling one to obtain a deeper oillet and expand the perforation area effectively.
Higher pressure drop of the nozzle yields higher impact pressure. Moreover, the strength and action range of impact pressure also increase with the increasing diameter of the nozzle's central hole. As a result, rock breakage could be improved by enlarging the nozzle's pressure drop and central hole size.
As the jetting distance increases, the strength of impact pressure declines, while its action range expands. An appropriate jetting distance is essential to achieving perfect impact pressure and action range for swirling-round SC-CO 2 jet flow, thus acquiring the desired rock-breaking effect.
Confining pressure, fluid temperature, and grooving area have a negligible influence on the performance of swirling-round SC-CO 2 jet flow.