#### Abstract

High-content H_{2}S gas reservoirs are important for natural gas extraction. However, the precipitation and deposition of elemental sulfur in high-sulfur-content gas reservoirs eventually lead to porosity and permeability damage, resulting in the low well productivity. Therefore, it is worth developing an accurate production prediction model considering sulfur deposition for fractured horizontal wells. In this study, based on the partition model and transient percolation theory, a novel numerical model considering the damage of sulfur deposition with pressure change on reservoir porosity and permeability was first developed to predict the production from fractured horizontal wells in high-sulfur-content gas reservoirs. Then, it was validated by actual field data from a high-sulfur-content gas reservoir. After that, the influence of sulfur deposition on the production of fractured horizontal wells was revealed through theoretical calculations, and the effects of hydraulic fracture parameters on production were analyzed. The results show that elemental sulfur gradually deposits in the reservoir pores as the reservoir pressure decreases and the production time increases, which eventually leads to permeability damage and reduces reservoir productivity; this negative impact gradually increases over time. It is also shown that the production of fractured horizontal wells increases with an increase in the half-length, fracture conductivity, and fracture number. Compared with the fracture half-length, the fracture conductivity and fracture number have a greater influence on the production of a single well. The model can handle the influence of nonlinear parameters caused by sulfur deposition, which allows accurate calculations and provides guidance for the development of fractured horizontal wells in gas reservoirs with high sulfur content.

#### 1. Introduction

During the development of high-sulfur-content gas reservoirs, elemental sulfur gradually precipitates from natural gas and deposits in reservoir pores with decreasing reservoir temperature and pressure, thus affecting gas well production [1]. Many scholars have studied the sulfur solubility model [2–5]. Fractured horizontal wells have a positive impact for delaying the pressure drop of a reservoir and alleviating sulfur deposition and precipitation; therefore, they are widely used in the development of high-sulfur gas reservoirs. According to the production data of fractured horizontal wells, there is a long-term linear flow in the production process in the log-log plot of production and time. Based on this feature, Cinco et al. [6] first proposed the concept of “bilinear flow,” and Cinco and Samaniego [7] introduced this concept to the research of the production of infinite conductivity vertical fracture wells. Lee and Brockenbrough [8] developed a trilinear flow model based on bilinear flow to analyze the influence of hydraulic fracture parameters on production. Ozkan et al. [9] applied the trilinear flow model to predict the production of fractured horizontal wells in shale gas reservoirs. The results show that the trilinear flow model can meet the seepage law of shale gas flowing into fractured horizontal wells. On this basis, researchers applied the trilinear flow model to dual media and triple media conditions and divided the model into fracturing and unreformed zones considering the influence of volume fracturing [10–13]. Stalgorova and Mattar [14] proposed a five-linear flow model for hydraulic fractures and induced fractures after hydraulic fracturing and divided the reservoir into five zones from the further division of the trilinear flow model. Deng et al. [15] further considered the impact of reservoir heterogeneity on the basis of Stalgorova’s model. Guo et al. [16] extended the five-linear flow model to a seven-linear flow model and divided the reservoir into seven regions. In addition to the partition model, the point source function method was applied to divide fractures into multiple nodes and then perform superposition to obtain the production performance of fractured horizontal wells [15, 17–19].

As the reservoir porosity and permeability in high-sulfur-content gas reservoirs exhibit nonlinear changes with decreases in pressure and the gas physical parameters also change, many scholars have proposed solutions to this nonlinear problem in order to simplify the mathematical equation. Agarwal [20] proposed the pseudo time method and combined the gas viscosity and compression coefficient with the time term. Kale and Mattar [21] and Kabir and Hasan [22] linearized the nonlinear term by the perturbation method based on the Boltzmann transformation. The results show that the second-order perturbation solution is sufficient for meeting engineering requirements. Ozkan et al. [9] and Apaydin [23] incorporated all nonlinear parameters by defining a new notion of pseudo pressure for nonlinear permeability changes. Eshkalak et al. [24] used the finite difference method to deal with multiple nonlinear factors and obtained good results. Hu et al. [25] developed a production model of horizontal wells in sulfur-bearing gas reservoirs based on the transient flow theory and superposition principle and described the influence of sulfur deposition on permeability, porosity, and production performance. Hu et al. [26] simplified the nonlinear term based on perturbation change and then analyzed the influence of sulfur deposition on the gas well production performance using the Douglas–Jones correction prediction, providing the dynamic characteristic curve of sulfur deposition prediction.

Although there are many studies on the production performance of fractured horizontal wells, further research on the production of fractured horizontal wells in high-sulfur-content gas reservoirs is still required. Hu et al. [26] used the stress sensitivity theory to deal with the change in porosity and permeability caused by sulfur deposition, but the sulfur solubility model is different under different temperature and pressure conditions, and the stress sensitivity model can thus be applied only for a specified reservoir. Therefore, there remain some problems in this method to address the nonlinear change of permeability. Additionally, the existing horizontal well model is mostly applied to describe conventional gas reservoirs, which cannot explain the production performance of fractured horizontal wells in sour gas reservoirs. In this study, based on the partition model and transient flow theory, the control equations of hydraulic fracture and reservoir zone are developed, and then, the control equation is solved by explicit processing of nonlinear parameters. A production prediction model for fractured horizontal wells in sour gas reservoirs is established, which can better deal with the nonlinear problems caused by sulfur deposition and avoid the complicated mathematical derivation, and the influence of sulfur deposition and engineering parameters on production is analyzed.

#### 2. Mathematical Model of Zonal Flow in a Fractured Horizontal Well

##### 2.1. Physical Model

In this study, the physical model of fractured horizontal well flow is shown in Figure 1. A quarter of any fracture in a fractured horizontal well is taken as the research object. According to the traditional trilinear flow hypothesis, the fluid flow zone can be divided into three parts ([9]), as shown in Figure 2: the flow in hydraulic fractures, the fracturing zone (part I), and the unreformed zone (part II). In order to fully consider the nonlinear mechanism of fluid flow in gas reservoirs with high sulfur content and address the boundary problem, the partition model needs to be further improved, and the following assumptions are made: (1)The reservoir is homogeneous, with a uniform thickness and infinite levels, and the upper and lower boundaries are closed and saturated by elemental sulfur(2)Constant pressure production(3)Single-phase transient two-dimensional plane flow in the reservoir and one-dimensional flow in hydraulic fractures(4)Natural gas physical parameters under the high pressure and the change in porosity and permeability caused by sulfur deposition are considered

**(a)**

**(b)**

##### 2.2. Mathematical Model

Based on the above-mentioned partition model, the differential discrete mathematical simulation method is used to address nonlinear coupling problems. The finite difference model is shown in Figure 3.

For part I, the right boundary and the lower boundary are closed, and the left boundary intersects with the hydraulic fracture. For part II, the upper boundary, right boundary, and left boundary are all closed. The lower boundary of the hydraulic fracture is a constant pressure boundary.

The control equation and boundary conditions for the finite conductivity fracture are expressed as

The difference equation can be obtained using the numerical difference method:

Equation (2) can be simplified and written as where , , , and . The lower boundary of the fracture is a constant pressure boundary:

In the matrix system, the governing equations and boundary conditions of fluid flow in part I and part II are

The difference equation can be obtained using the numerical difference method:

Equation (6) can be simplified and written as where , , , , and

For closed boundaries,

At the boundary between part I and part II, considering the continuity of fluid flow and pressure wave, both sides of the boundary can be considered to be of equal pressure and equal rate. The nonlinear parameters in the model are explicitly processed, and the pressure is implicitly processed. The coefficients of the fracture zone and matrix zone are expressed in the form of a diagonal matrix. The pressure and rate at each time step are solved by a mathematical algorithm, and then, the permeability and porosity of the grid are updated considering the influence of sulfur deposition.

The production of a single fracture is calculated by the difference between the bottom grid of the fracture and the bottom hole pressure:

#### 3. Sulfur Solubility Prediction Model

Based on the study of Chrastil [2], Roberts [3] proposed an empirical formula to predict the solubility of sulfur using a large amount of experimental data.

Hu et al. [4] analyzed the variation of sulfur solubility in natural gas with pressure and temperature: where the parameters , , and are shown in Table 1.

With a decrease in temperature and pressure, sulfur precipitates from natural gas, and the reservoir permeability and porosity decrease due to sulfur deposition. The variation in the sulfur solubility can be expressed as

The precipitated elemental sulfur exists in the pores in the form of solid deposition, and the porosity is

According to Roberts [3], the relationship between permeability and porosity in gas reservoirs with high sulfur content is

By introducing equations (13) and (14) into equation (15), the change in the reservoir permeability with pressure considering the effect of sulfur deposition can be obtained by

The coefficients , , and in equation (16) can be selected according to different sulfur concentrations. Specific parameters can be seen in Table 1. Compared with the matrix, hydraulic fractures have a higher flow capacity, and the influence of sulfur deposition on permeability in hydraulic fractures is not considered.

If the stress sensitivity effect is taken into account, substitute it into the following equation after calculating equation (16):

#### 4. Results and Discussion

In order to validate the proposed model, the production data of a sulfur-bearing gas reservoir in Sichuan, China, were selected. The natural gas components are H_{2}S (11%), CH_{4} (77%), CO_{2} (8%), and N_{2} (2%). The natural fractures in the reservoir do not develop, and the reservoir permeability between the hydraulic fractures does not increase significantly, so part I and part II have the same reservoir physical properties. The specific parameters used in the model are shown in Table 2, and the calculation results are displayed in Figure 4. It can be seen from the figure that the calculation results of the model in this study are roughly consistent with the actual production data, validating the proposed model.

##### 4.1. Effects of Sulfur Deposition on Well Production

Figure 5 shows the effect of sulfur deposition on fractured horizontal well production at different H_{2}S concentrations and bottom hole pressure. As shown in Figure 5, the influence of sulfur deposition on production can be divided into two categories: one is the production pressure difference and the other is the concentration of H_{2}S. When the bottom hole pressure is 19 MPa, that is, the production pressure difference is 20 MPa; although the concentration of H_{2}S is only 5%, the effect of sulfur deposition on production is still obvious. And the decline gradually increases with the extension of production time. This indicates that the influence of sulfur deposition on production should be considered in the development of sulfur-bearing gas reservoirs. The effect of sulfur deposition on production was minimal when the bottom hole pressure was 30 MPa (the pressure difference was 9 MPa), and the H_{2}S concentration was 5%. Therefore, it can be considered that in this example, when the concentration of H_{2}S is less than 6% and the production pressure difference is less than 9 MPa, the effect of sulfur deposition on production can be ignored.

Figure 6 shows the reservoir permeability in the 1/4 fractured zone after ten years of production. The direction is along the wellbore direction, and the direction is along the fracture propagation orientation. As displayed in Figure 6, after ten years of production, due to the decrease in reservoir pressure, sulfur components dissolved in natural gas gradually precipitate and deposit, resulting in a decrease in the reservoir permeability, and the magnitude of the decrease in permeability near the junction of the fracture and wellbore is the largest, which eventually reduces the production of fractured horizontal wells. Therefore, in the production process of fractured horizontal wells in high-sulfur-content gas reservoirs, it is necessary to control the production pressure difference to reduce sulfur deposition.

Figure 7 shows the reservoir porosity in the 1/4 fractured zone after ten years of production. As displayed in Figure 7, after ten years of production, the porosity of the reservoir gradually decreases, and the closer it is to the fracture, the porosity decreases more. This trend is also consistent with Figure 6.

##### 4.2. Effect of Hydraulic Fracture Parameters

This proposed model is used to simulate and calculate the production change of fractured horizontal wells after one year of production. The effects of fracture half-length, fracture conductivity, and fracture number on production are analyzed as shown in Figures 8, 9, and 10, respectively.

It can be observed from Figure 8 that with a longer fracture half-length, the well productivity is higher, but the magnitude of increase in initial production gradually decreases. A larger fracture half-length leads to a larger volume of the transformation area, which benefits the production of fluid from the reservoir. But the longer the fracture is, the greater the pressure loss in the fracture will be, compared with the fracture conductivity and fracture number; the influence of fracture half-length on well production is less obvious.

Figure 9 shows the influence of fracture conductivity on production. As the fracture conductivity increases, the production of fractured horizontal wells increases. Furthermore, higher production occurs not only during early production but also in the stable production period.

The fracture number has an influence on the production of fractured horizontal wells: the production increases as the fracture number increases. A greater number of fractures lead to a larger volume of the transformation zone. This implies that the fracture number is beneficial for production.

In order to study the effect of different fracture parameters on sulfur deposition, the permeability of reservoir is used to reflect sulfur deposition in this paper. Figures 11, 12, and 13 have been studied to illustrate the effects of different fracture properties on reservoir permeability. Figure 11 shows the effect of different fracture half-lengths on reservoir permeability. From top to bottom, layers show fracture half-lengths of 25 m, 50 m, and 75 m. It can be seen from the figure that the reservoir permeability decreases gradually as the fracture half-length increases, because the larger the fracture half-length, the higher the production. In the case of constant pressure production, this means that the production pressure difference is larger, so the sulfur deposition in the reservoir also increases correspondingly and the reservoir permeability decreases.

Figure 12 shows the effect of different fracture conductivity on reservoir permeability. From top to bottom, the layers show that the fracture conductivity is 100D mm, 150D mm, and 200D mm. As can be seen from the figure, the reservoir permeability decreases gradually with the increase in fracture conductivity, especially near the fractures, because with the increase in fracture conductivity, the greater the pressure drop near the fracture, the more the sulfur deposition.

Figure 13 shows the influence of different fracture numbers on reservoir permeability. From top to bottom, the layers show that fracture numbers are 3, 5, and 7. It can be seen from the figure that the reservoir permeability gradually decreases as the number of fractures increases. Because when the length of horizontal well remains unchanged, the more fractures there are, the smaller the spacing between fractures and the greater the pressure drop between fractures. Therefore, the more obvious the difference of reservoir permeability in the direction will be.

#### 5. Conclusions

In this study, a novel numerical model considering the damage of sulfur deposition was developed, and the influence of sulfur deposition on the production of fractured horizontal wells was analyzed. Sensitivity analysis of different fracture parameters was conducted. The following conclusions can be drawn: (1)Considering the nonlinear variation characteristics of permeability in high-sulfur gas reservoirs, the “implicit pressure explicit nonlinear parameter” method is used to predict the production of fractured horizontal wells in high-sulfur-content gas reservoirs. The new model can reflect the production performance and development characteristics of fractured horizontal wells in gas reservoirs with high sulfur content(2)The accuracy of the model was verified by an example of a high-sulfur-content gas reservoir, and the influence of sulfur deposition on the production of fractured horizontal wells was revealed. As the production time increases, elemental sulfur is gradually deposited in the reservoir pores due to the decrease in reservoir pressure, resulting in a decrease in reservoir porosity and permeability, eventually leading to a decrease in production(3)The effects of hydraulic fracture half-length, fracture conductivity, and fracture number on the production of fractured horizontal wells were analyzed. With an increase in fracture half-length, fracture conductivity, and fracture number, the production of an individual well increases. Compared with the fracture half-length, the fracture conductivity and fracture number have a greater influence on well production

#### Nomenclature

: | Sulfur solubility (g/m^{3}) |

: | Initial sulfur solubility |

: | Rock compressibility (Pa^{-1}) |

: | Reservoir thickness (m) |

: | Fracture permeability (m^{2}) |

: | Matrix permeability (m^{2}) |

: | Pseudo pressure |

: | Molar mass (kg/mol) |

: | Standard pressure (Pa) |

: | Initial pressure (Pa) |

: | Wellbore pressure (Pa) |

: | Gas rate (m^{3}/s) |

: | Universal gas constant (8.315 Pa·m^{3}/mol/K) |

: | Standard temperature |

: | Production time (s) |

: | Fracture width (m) |

: | Boundary length in the direction (m) |

: | Fracture half-length (m) |

: | Boundary length in the direction (m) |

: | Gas compressibility factor |

: | Matrix porosity |

: | Initial matrix porosity |

: | Space step in the direction (m) |

: | Space step in the direction (m) |

: | Time step (s) |

: | Stress effect coefficient (Pa^{-1}) |

: | Gas density (kg/m^{3}) |

: | Solid sulfur density (kg/m^{3}) |

: | Gas viscosity (Pa·s) |

: | Gas-specific gravity. |

#### Data Availability

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

#### Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

#### Acknowledgments

The authors would like to acknowledge the support provided by the National Natural Science Foundation of China (52074248, 51804282) and Fundamental Research Funds for the Central Universities (35832019033, 35832020095).