Numerical Investigation on Wellbore Temperature Prediction during the CO2 Fracturing in Horizontal Wells

A novel model is established to predict the temperature field in the horizontal wellbore during CO2 fracturing. The pressure work and viscous dissipation are considered, and the transient energy, mass and momentum equations as well as the CO2 physical properties are solved fully coupled. The model passes the convergence test and is verified through a comparison using the COMSOL software. Then, a sensitivity analysis is performed to study the effects of the treating parameters. Results illustrate that the relationship between the injection rate and the stable bottom-hole temperature (hereinafter referred to as BHT) is non-monotonic, which is different from the hydraulic fracturing. The existence of the horizontal section will increase the BHT at 2 m3/min condition but reduce the BHT at 10 m3/min condition. The problem of high wellbore friction can be alleviated through tube size enhancement, and the ultimate injection rate allowed increased from 2.7 m3/min to 29.6 m3/min when the tube diameter increased from 50.3 mm to 100.3 mm. Additionally, the open-hole completion method of the horizontal section can increase the BHT to 2.7 ◦C but reduce the near formation temperature to 24.5 ◦C compared with the casing completion method.


Introduction
Shale gas is one of the most important unconventional oil and gas resources [1,2]. Slickwater fracturing and horizontal well drilling are the basic technologies for the development of the shale gas reservoirs [3]. However, the usage of the water-based fracturing fluid can generate formation damage and consume a large amount of water resources. CO 2 fracturing can simultaneously avoid the above disadvantages and achieve the effects of carbon sequestration (CCS) [4,5] and gas recovery enhancement (EGR) [6][7][8][9], which is a potential alternative technology to the slick-water fracturing. Unlike the water-based fracturing fluid, the physical properties of CO 2 vary significantly with temperature and pressure [10,11], which can affect the sand-carrying capacity and the performance of the stimulation. Additionally, it is impractical to modify the CO 2 temperature during the fracturing process based on the temperature monitor data because of the expense and time consumed. Therefore, accurately predicting the temperature and pressure of the horizontal well CO 2 fracturing process is necessary in terms of both the environmental and economic aspects.
The study of the wellbore temperature field can be traced back to 1937 [12]. After a long-term development, the wellbore temperature prediction methods can be divided into three categories. The first category is the quasi-steady numerical method represented by the Ramey's model [13], which assumes that the fluid and wellbore zones are steady-state heat transfer, and the formation zone is transient heat transfer. The quasi-steady numerical wells. Moreover, an efficient wellbore temperature prediction model remains lacking for horizontal wells, which are extensively used in unconventional oil and gas reservoirs.
In this article, a transient fully coupled model of wellbore temperature for CO 2 fracturing in horizontal wells is established, which combines the 1D non-isothermal flow in the axial direction with the 1D transient heat transfer in the radial direction. The effects of the pressure work and viscous dissipation are considered in the energy equation, and the transient mass and momentum equations as well as the CO 2 physical property models are solved simultaneously. The 1D + 1D model is verified through a convergence test and a comparison using the COMSOL software. Finally, the parameter analysis is performed to provide guidance for the treatment of the CO 2 fracturing in horizontal wells. Figure 1 illustrates the physical model. The wellbore is divided into vertical, curved, and horizontal sections. In the vertical and curved sections, the wellbore structure includes the tubing, annulus, casing, and the cement sheath. In the horizontal section, as the openhole completion method is considered here, the wellbore structure only includes tubing and annulus.

Basic Assumption
The horizontal well CO 2 fracturing has three characteristics. First, the aspect ratio of the wellbore of oil and gas wells is usually very large (1 × 10 3 to 1 × 10 4 ) given the enormous gap between the length of the axial direction and the width of the radial direction. Second, the injection rate is typically above 3 m 3 /min, which results in a high Reynolds number (Re > 1 × 10 6 ) and a turbulent flow during fracturing. Last but not least, the L-shaped structure of the horizontal wellbore cannot be simplified through the method of axisymmetric modeling as the vertical well in 2D or 3D space. All of the above necessitate the meshes in the 2D, or the 3D Computational Fluid Dynamics (hereinafter referred to as CFD) model becomes a huge number to obtain stability and convergence. An example is the k-ε model. The dimensionless wall distance function y+ of the first boundary layer grids should be 30-200 [58,60], and the dimensionless wall distance function can be calculated as y+ = v∆y/µ, where y is the thickness of the first boundary layer (please notice that parameters that are not further explained in the main text hereinafter are listed at the Nomenclature section). For CO 2 fracturing, this process requires the thickness of the first boundary layer to be less than 0.001 m (as shown in Figure 2). To solve the above problems, considering the following two points: (1) the radial dimension is pretty small compared with the axial dimension, (2) the effect of the axial heat conduction on the wellbore temperature can be ignored [53,54], the following assumptions are used:

1.
The CO 2 zone is considered a one-dimensional axial non-isothermal flow; 2.
The wellbore and formation zones are considered a one-dimensional radial transient heat transfer; 3.
The heterogeneity and anisotropy of the formation are ignored; 4.
At the initial moment, the wellbore is full of fluid, and the fluid is stationary.

Mathematical Model
The calculation zone includes the CO 2 , wellbore, and formation zones. The CO 2 zone is considered a non-isothermal flow, and the continuity equation, momentum equation, as well as the energy equation must be solved simultaneously. The wellbore zone is considered a fluid natural convective heat transfer, and the formation zone is considered a solid heat transfer process, for which the energy equation must be solved.
Continuity equation of the CO 2 zone: Momentum equation of the CO 2 zone: Energy equation of the CO 2 zone: where f D is the friction coefficient. The fitting formula (R 2 > 0.999) [61] of the CO 2 fracturing wellbore friction plate [62] is used here to obtain the expression of f D .
Q wall , Q fric and Q pres are the heat flows of the unit axial length generated by the heat transfer, viscous dissipation, and pressure work, respectively, W/(m 2 ·K). They are defined as follows: where h t is the heat transfer coefficient of the forced heat convection. The Gnielinski correlation, which is valid for 0.5 ≤ Pr ≤ 2000 and 3000 ≤ Re ≤ 5 × 10 6 , is applied here to calculate h t [63]: Pr is the Prandtl number, defined as Pr = C p µ/λ. Unless otherwise specified, the heat flow mentioned in the following refers to the heat flow of the unit axial length.

Energy Equation of the Wellbore Zone
The wellbore zone includes the tubing, annulus, casing, and the cement sheath. Except for the annulus, the remaining regions are all solid, and only the heat conduction requires consideration. For the annulus, the internal fluid is in a state of natural convection, which is characterized here by the natural convection heat transfer coefficient. The overall governing equation of the wellbore zone is as follows: In the formula, ρ, C p and λ take the values of different materials according to different spatial positions, and the temperature at the interface of different zones is continuous. The specific form of the governing equation of each part is different, and the wellbore structure in the horizontal section is also different from the vertical and curved sections. The governing equations for each part of the wellbore are given below.
Vertical Section (Curved Section) Figure 3 presents the radial grids of the vertical section (curved section), in which the white point with blue edge indicates the location of the various radius of the tubing, casing, cement sheath, and formation, while the red point with black edge indicates the location of the calculation temperature. The energy equation of the tubing: The energy equation of the annulus: The energy equation of the casing: The energy equation of the cement sheath: where h an is the natural convection heat transfer coefficient in annulus, W/(m 2 · • C), defined as h an = 0.049(GrPr) 1/3 Pr 0.074 λ an r to ln(r ci /r to ) .
Gr is the Grashof number of the annulus fluid, defined as where λ c,h and λ h,r are the composite thermal conductivity, which can have the following definition based on Fourier law: Horizontal Section Figure 4 illustrates the radial grids of the horizontal section. The implications of the two kinds of points are the same as in Figure 3. The wellbore in the horizontal section only includes the tubing and the annulus. The energy governing equation of the tubing is the same as the vertical (curved) section. The annulus energy governing equation is as follows:

Energy Equation of the Formation Zone
No flow occurs in the formation. Accordingly, the energy governing equation in the formation can be simplified to the heat conduction equation: Owing to the difference between the radial grids of the vertical (curved) and horizontal sections, the equations at the boundary of the formation and the wellbore are different.
The first layer grids of the vertical (curved) section: The first layer grids of the horizontal section:

Initial Conditions
The initial wellbore temperature is calculated from the original geothermal gradient. As for the pressure, considering that the wellbore is assumed to be full of CO 2 , the initial wellbore hydrostatic pressure can be obtained through an iterative method combined with the CO 2 density model. The governing equations for the initial conditions are as follows:

Boundary Conditions
For the flow field in the CO 2 zone, a given injection rate at the wellhead is used as the velocity boundary, and a pressure boundary is required at the bottom of the well. The current study uses the fracture closure pressure. For the temperature field of the CO 2 zone, the injection temperature at the wellhead is given as the inlet boundary, and the bottom-hole of the well is insulated. However, the thermal convection is considered at the bottom-hole. The governing equations for the boundary conditions are as follows: (24)

Solving Method
In the 1D + 1D transient fully coupled model, the continuity, momentum, and energy equations are coupled together through the interaction among pressure, temperature, and CO 2 physical properties to form a transient nonlinear system, which requires an iteration solving method. The Span-Wagner EoS [44], Fenghour model [45] and Vesovic model [46] are used here as the physical properties models, which are proven to be highly accurate [10,11]. Two iteration loops are set up during the solution process: the pressurespeed and temperature iterations. Figure 5 shows the model solution flowchart.

Test of the Convergence
This paper used the iterative method to solve the velocity, pressure, and temperature. To test the convergence of the numerical algorithm, the changes in the calculation results must be analyzed during the iteration. Given that the velocity-pressure iteration is nested in the temperature iteration, the convergence test was performed on the temperature calculation results. Table 1 details the parameters of the benchmark case. The convergence condition of the temperature calculation is defined as where X 2 represents the 2-norm of vector X. Figure 6a presents the T k f − T k−1 f 2 with various iteration steps and the number of iterations used to reach the convergence condition at different time steps. Figure 6b indicates that as the number of iteration steps increased, the error continued to decrease; as the time step increased, the number of iteration steps used to achieve convergence condition also continued to decrease. Thus, the numerical solution method used in this paper met the requirements of convergence.

Comparison of the Non-Isothermal Flow Simulation in a Horizontal Tube with COMSOL
Currently, measured data of CO 2 fracturing are lacking in horizontal wells. Accordingly, a comparison method with commercial software was used to verify the calculation results of the 1D + 1D transient fully coupled model in the horizontal section. COMSOL is a multi-physics coupled finite element software. The non-isothermal tube flow module can simulate thermal convection and heat conduction in the tube. However, in COMSOL, the wellbore wall layer (tubing, annulus, casing, cement sheath, etc.) notably adopts the steady-state heat transfer assumption and does not consider the natural convection in the annulus [64].
The calculation parameters were based on Table 1. The vertical section was ignored, and the horizontal section was extended to 1000 m. Simultaneously, the initial values were set as the parameters of the depth of 1600 m, implying that the initial temperature was 68 • C, and the initial pressure was 32 MPa, based on Table 1. As the horizontal wellbore was completed with the open-hole completion method, the casing and cement sheath were removed. The natural convection in the annulus was also ignored in the 1D + 1D fully coupled model for comparative analysis.  Figure 7a shows the variation of the BHT with time. Apparently, after ignoring the natural convection in the annulus, the results of the 1D + 1D transient fully coupled model were close to the results of COMSOL. However, at the early stage, the change rate of BHT calculated by COMSOL was faster than that of the 1D + 1D transient fully coupled model because of the steady-state wall heat transfer assumption used in COMSOL. The latter ignored the temperature stabilization process of each wall surface of the wellbore and underestimated the CO 2 temperature. However, due to the large injection rate and the neglection of the natural convection, the effect of the wall heat transfer on the CO 2 temperature was small; the temperature difference between the two was only 0.2 • C after 120 min. Considering that the treatment of wall heat transfer was the main difference between COMSOL and the 1D + 1D transient fully coupled model, another comparison was performed in which the wall heat transfer term was removed both from COMSOL and the 1D + 1D transient fully coupled model. As shown in Figure 7b, the calculation results of the two were highly consistent after neglecting the effect of wall heat transfer, and the deviation did not exceed 0.03 • C. This finding validated the model accuracy in the horizontal section.

Results and Discussions
This section presents an analysis of the effects of the treating parameters using the 1D + 1D fully coupled model. Table 1 details the parameters of the benchmark case.

Analysis of the Injection Rate
The cases were considered when the injection rates were 2 m 3 /min, 6 m 3 /min, and 10 m 3 /min, based on the benchmark case. 4.1.1. Analysis of the Time Domain Figure 8 shows the variation of the BHT with time under various injection rates. When the injection time was less than 50 min, an increase in the injection rate speeded up the decrease in the temperature. This finding was due to the change in the convective heat transfer efficiency with the CO 2 velocity. Figure 9 shows the average CO 2 velocity under different injection rates. As the size of the tubing was unchanged, the velocity was basically linear with the injection rate. However, faster convective heat transfer will also hasten the stability of the CO 2 temperature. Using 0.03 • C/min as the criterion of the temperature stability, Figure 10 shows that the temperatures of the cases of 2 m 3 /min, 4 m 3 /min, and 6 m 3 /min did not reach the stable state after injection for 120 min. The stable temperature of CO 2 was related to the heat transfer mechanism determined by the combination of the pressure work, viscous dissipation, and wall heat transfer. The analysis of space domain will illustrate this process. Figure 11 shows the variation of the wellhead pressure (hereinafter referred to as WHP) with time under various injection rates. Given that no flow occurred at the initial time, the initial WHP was a constant with the value of 17.5 MPa. When fracturing started, CO 2 flowed into the wellbore and generated friction. To keep the BHP constant, the WHP increased. The enhancement was determined by the friction, which increased significantly as the injection rate increased. Figure 12 shows that the wellbore friction reached 67.7 MPa at 10 m 3 /min. In addition, Figure 11 shows that the pressure increased to a maximum within a short period of time, and then the pressure decreased slightly. This phenomenon was caused by the variation of the CO 2 density, which was affected by the CO 2 temperature. As shown in Figure 13, after 120 min, the CO 2 density at the bottom-hole increased by 185 kg/m 3 compared with that after 1 min.       Figure 14 shows the temperature profiles of CO 2 with various injection rates after 120 min. With the increase in the measured depth, the axial direction and the wellbore structure changed. This was divided into three sections: vertical, curving, and horizontal sections. In the vertical section, the temperature of CO 2 decreased first and then increased with the increase in the injection rate. This phenomenon was caused by the heating effect of the viscous dissipation. As shown in Figure 15a, the heat flow generated by the viscous dissipation was only 120 W/m at 2 m 3 /min but reached 4796 W/m at 10 m 3 /min. In the curved section, the well inclination angle slowly changed from 0 • to 90 • , and the gradient of the hydrostatic pressure gradually disappeared. Thus, the total pressure drop in the curved section was significantly increased, which resulted in an increase in the cooling effect of pressure work. As shown in Figure 15b    Therefore, in the curved section, the cooling effect of pressure work was obviously enhanced, leading to a decrease in the temperature gradients of all three cases. In the horizontal section, the absolute value of the Q pres under the injection rate of 10 m 3 /min continued to increase, whereas the Q pres under the injection rate of 2 m 3 /min did not change much. Simultaneously, as shown in Figure 15c, compared with the vertical and curved sections, the Q wall in the horizontal section increased significantly, which was due to the open-hole completion method. Section 4.4 discusses the effects of horizontal completions. Generally, with the combined effects of the pressure work and wall heat transfer, the temperature gradient in the horizontal section was higher than that of the vertical section under the 2 m 3 /min case but lower than that of the vertical section under the 10 m 3 /min case. The 6 m 3 /min case fell in between the other two cases. Therefore, compared with the vertical section, the existence of the horizontal section increased the BHT under a low injection rate but reduced the BHT under a high injection rate. Figure 16 shows the pressure profiles of CO 2 with various injection rates after 120 min. CO 2 pressure was determined by the hydrostatic pressure and wellbore friction. In the vertical section, when the injection rate was 2 m 3 /min, the hydrostatic column pressure was higher than the wellbore friction, resulting in a positive pressure gradient. Additionally, the cases of 6 m 3 /min and 10 m 3 /min were contrary to the 2 m 3 /min case. In the horizontal section, as the hydrostatic pressure no longer increased, the CO 2 pressure decreased along the axial direction under any injection rates. The curved section was the transition between the vertical and horizontal sections, during which the hydrostatic column pressure gradient slowly decreased to zero.

Analysis of the Injection Temperature
The cases where the injection temperatures were −20 • C, 0 • C and 20 • C were considered on the bases of the benchmark case. Figure 17 shows the variation of the BHT with time under various injection temperatures. Evidently, the increase in the injection temperature effectively enhanced the BHT. When the injection temperature increased to 20 • C, the BHT maintained above the supercritical temperature during the entire fracturing process. Figure 18 shows the variation of the WHP with time under various injection temperatures. As the lower temperatures increased the CO 2 density, which enhanced the hydrostatic pressure, the WHP of the −20 • C case was 1.8 MPa lower than that of the 20 • C case. Figure 19 shows the temperature profiles of CO 2 with various injection temperatures after 120 min. When the injection temperature increased from −20 • C to 20 • C, the temperature of CO 2 increased significantly. However, the temperature gradient of CO 2 decreased continuously, especially in the horizontal section. Figure 20 indicates that when the temperature increased from −20 • C to 20 • C, the temperature difference between the wellhead and the bottom-hole decreased by 5.2 • C. This phenomenon was caused by the effects of the injection temperature on Q pres and Q wall . As shown in Figure 21a,b, respectively, the heating effect of Q wall reduced, whereas the cooling effect of Q pres increased with the enhancement of the injection temperature. Therefore, with the increase in the injecting temperature, the CO 2 temperature gradient decreased. Figure 22 illustrates the pressure profiles of CO 2 with various injection temperatures after 120 min. As stated above, the hydrostatic pressure increased with the reduction of the temperature. Therefore, the lower the temperature, the higher the pressure in the vertical section. However, given that the hydrostatic pressure gradient disappeared, the pressure profiles of the three cases in the horizontal section were slightly different.

Analysis of the Tube Inner Diameters
The cases where the tube inner diameters were 62, 76, and 100.3 mm were considered on the bases of the benchmark case. Table 2 shows the specific size of the tube.    Figure 23 shows the variation of the BHT with time under various tube inner diameters. The smaller the inner diameter of the tube, the faster the cooling rate. However, the CO 2 temperature of the 62-mm case was 11.8 • C higher than that of the 100.3-mm case after 120 min. This phenomenon was caused by the variation of the convective heat transfer efficiency and Q fric . As shown in Figure 24, the flow velocity enhanced obviously with the decrease in the tube inner diameter, which resulted in a higher convective heat transfer effectiveness. However, the extra friction from the higher flow velocity Q fric increased by 3489 W/m, which made the lowest temperature higher in the 62-mm case. Figure 25 shows the variation of the WHP with time under various tube inner diameters. Its effects on WHP were similar to that of the injection rate. When the inner diameter decreased from 100.3 mm to 62 mm, the wellhead pressure increased by 79.6 MPa.    Figure 26 shows the temperature profiles of CO 2 with various tube inner diameters after 120 min. In the vertical section, when the inner diameter increased, the effect of viscous dissipation decreased and the CO 2 temperature decreased. However, in the horizontal section, the temperature gradient of the 62-mm case decreased significantly. Although the viscous dissipation increased with the increase in the wellbore friction, the cooling effect of the pressure work also became increasingly obvious. The combined effect of the two factors determined the temperature gradient of CO 2 . Figure 27a provides the sum of Q fric and Q pres under various tube inner diameters after 120 min, indicating that the smaller the inner diameter, the greater the heat flow gradient, especially in the horizontal section. The heat flow reduction rate was very fast in the 62-mm case: at a depth of 2300 m, the cooling effect of the pressure work exceeded the warming effect of the viscous dissipation, which created a negative sum of the two heat flows. Although the total heat flow remained positive throughout the wellbore due to the effect of the wall heat transfer (as shown in Figure 27b), the temperature gradient was anticipated to be negative along the wellbore if the horizontal section continued to extend in the 62-mm case.   Figure 28 indicates the pressure profiles of CO 2 with various tube inner diameters after 120 min. In the vertical section, the hydrostatic pressure of the 100.3-mm case was greater than the wellbore friction, and the pressure gradient was positive; the wellbore friction of the 76-and 62-mm cases was greater than the hydrostatic pressure, and the pressure gradient was negative. In the horizontal section, as the hydrostatic pressure gradient was 0, the pressure gradients of the three cases were controlled by the wellbore friction, which made them all negative. The simulation results showed that the wellbore friction could be reduced significantly by increasing the inner diameter. Therefore, by increasing the inner diameter, the maximum of the CO 2 injection rate could be increased, which is beneficial for improving the sand-carrying capacity of the CO 2 fracturing fluid. At present, the newly developed 3000 type fracturing pump can stably provide 140 MPa pumping pressure [65]. Under the parameters of the benchmark case, using the wellhead pressure of 140 MPa as the ultimate pressure, the ultimate injection rate corresponding to different inner diameters was obtained as shown in Figure 29. Apparently, when the inner diameter was enlarged from 50.3 mm to 100.3 mm, the ultimate injection rate was increased from 2.7 m 3 /min to 29.6 m 3 /min.

Analysis of the Completion Method of the Horizontal Section
The cases where the completion methods of the horizontal section were casing completion and open-hole completion were considered on the bases of the benchmark case. The main difference between the two methods was that casing and cement sheath existed between the annulus and the formation in the casing completion, which were otherwise absent in the open-hole completion. Figure 30 shows the variation of the BHT with time under various completion methods. The BHT of the open-hole completion in the horizontal section was 2.7 • C higher than that of the casing completion after 120 min. The variation of the annulus natural convection explained this phenomenon. In the open-hole completion case, the radius of the annulus was much larger than that of the casing case, which enhanced the efficiency of the natural convection. As shown in Figure 31a, compared with the casing completion case, the average annulus natural convection heat transfer coefficient of open-hole completion was 30 W/(m 2 · • C) higher. The enhancement of the natural convection in the annulus improved the wall heat transfer efficiency. Figure 31b indicates that the average Q wall with open-hole completion method was 313 W/m higher than that of the casing completion method. Figure 32 shows the variation of the WHP with time under various completion methods. Apparently, the completion method had little effect on the WHP. This finding is further discussed in the analysis of the space domain.     Figure 33 shows the temperature profiles of CO 2 with various completion methods after 120 min. In the vertical and curved sections, the temperature profiles coincided, and the difference between the two cases only appeared after entering the horizontal section. As discussed above, the difference in the Q wall was the main cause of the temperature variation. Figure 34 provides the Q wall profiles with different horizontal section completion methods, which indicated that Q wall was the same in the vertical and curved sections, whereas a jump (∆Q wall = 817 W/m) occurred at the boundary between the curved and horizontal sections, which increased the CO 2 temperature gradient in the horizontal section in the open-hole completion method. Therefore, open-hole completion for the horizontal section was beneficial to increase the BHT.   Figure 35 shows the variation of the CO 2 pressure profiles with time under various completion methods after 120 min. Firstly, as the completion method did not change the flow characteristics, the wellbore friction was the same at different wellbore sections. Secondly, the variation of CO 2 density occurred at the horizontal section, in which the hydrostatic pressure gradient was 0. Therefore, the completion method had little effect on CO 2 pressure in both the time and space domains.

Analysis of the Space Domain
In addition, the usage of the open-hole completion method also significantly increased the cooling effect of CO 2 on the formation. As shown in Figure 36a,b, after injecting 120 min, the formation temperature near the wellbore of the open-hole completion case was lower by 24.5 • C than in the casing completion case.

Conclusions
A 1D + 1D numerical model was established to investigate the horizontal wellbore temperature field in CO 2 fracturing, which can provide an inexpensive and efficient way to investigate and guide the CO 2 fracturing design. The convergence test and the comparison with the COMSOL software were performed to verify the model. The following conclusions were drawn on the bases of the assumptions and conditions of this work: 1.
The higher the injection rate, the faster the reduction of temperature and the shorter the time for the BHT to reach stability. The relationship between the stable BHT and the injection rates is determined by the heat transfer mechanisms, which are non-monotonic. Under the parameters of the benchmark case, after 120 min, the corresponding injection rates were 2 m 3 /min, 10 m 3 /min, and 6 m 3 /min, according to the sequence of BHT from large to small. 2.
In the horizontal section, the cooling effect of the pressure work and the heating effect of the wall heat transfer are enhanced. Under the combined effect of the two, the temperature gradient in the horizontal section was higher than that in the vertical section when the injection rate was 2 m 3 /min. The opposite was the case where the injection rate was 10 m 3 /min. Therefore, compared with the vertical section, the existence of the horizontal section will increase the BHT when the injection rate is low and reduce the BHT when the injection rate is high.

3.
Wellbore friction increases exponentially with the increase in the injection rate. Enhancing the tube diameter can effectively reduce the wellbore friction, thus extra attention to casing fracturing is necessary. Based on the benchmark case and taking WHP at 140 MPa as the ultimate pressure, the ultimate injection rate increased from 2.7 m 3 /min to 29.6 m 3 /min when the tube diameter increased from 50.3 mm to 100.3 mm.

4.
The open-hole completion method in the horizontal section enlarges the annulus space more than the casing completion method, which improves the natural convection heat transfer efficiency of the annulus. Based on the benchmark case, the BHT of the open-hole completion in the horizontal section was 2.7 • C higher than that of the casing completion after 120 min, whereas its effects on pressure were negligible. In addition, the open-hole completion in the horizontal section will significantly enhance the cooling effect of CO 2 on the formation near the wellbore. After 120 min, the formation temperature of the bottom-hole with the open-hole method was 24.5 • C lower than that in the casing completion method. Informed Consent Statement: Not applicable.

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest:
The authors declare no conflict of interest.