Zoning Productivity Calculation Method of Fractured Horizontal Wells in High-Water-Cut Tight Sandstone Gas Reservoirs under Complex Seepage Conditions

: Tight sandstone gas reservoirs generally contain water. Studying the impact of water content on the permeability mechanism of tight gas reservoirs is of positive signiﬁcance for the rational development of gas reservoirs. Selected cores from a tight sandstone gas reservoir in the Ordos Basin were used to establish the variation in its seepage mechanism under different water saturations. The experimental results show that the gas slip factor in tight water-bearing gas reservoirs decreases as the water saturation increases. The stress sensitivity coefﬁcient and the threshold pressure gradient (TPG) increase with increasing water saturation, characterizing the relationships between stress sensitivity coefﬁcients, TPG, permeability, and water saturation. As the water saturation gradually increases, the relative gas phase permeability of tight sandstone gas reservoirs will sharply decrease. When the water saturation exceeds 80%, the gas phase permeability becomes almost zero, resulting in gas almost ceasing to ﬂow. Through the analysis of experimental results, we deﬁned high-water-cut tight sandstone gas reservoirs and analyzed the permeability characteristics of high-water-cut tight sandstone gas reservoirs in different regions. Combining stress sensitivity coefﬁcients and the TPG with permeability and water saturation relationships, we established a zoning productivity calculation method of fractured horizontal wells in high-water-cut tight sandstone gas reservoirs under complex seepage conditions and validated the practicality of the model through example calculations.


Introduction
Countries around the world possess numerous tight sandstone natural gas reservoirs, displaying significant development potential.For a considerable time to come, they will continue to play a leading role in the field of natural gas development.Non-tight sandstone gas, as an unconventional natural gas resource, is found within tight sandstone reservoirs, typically with a permeability of less than 0.1 mD [1][2][3][4].Such reservoirs generally lack natural productivity or have extremely low productivity, necessitating hydraulic fracturing [5].With the development and transformation of tight sandstone reservoirs in progress, gas wells often experience water production phenomena [6].These phenomena affect the seepage mechanisms, yet there are currently no mature solutions available for quantitative assessment.Therefore, it is imperative to establish a new model for the dual-phase hydraulic fracturing productivity of tight gas reservoirs, considering multiple factors simultaneously, based on a deeper understanding of the impact of water saturation on seepage mechanisms.
When gas flows through porous media, it is influenced by the slip flow effect [7,8].This is because the interaction forces between gas molecules and solid molecules are relatively weak.At the pipe wall, gas molecules maintain certain movement states and undergo directional motion along the wall through momentum exchange, leading to the directional flow of adjacent gas molecules [9].The water saturation in tight sandstone gas reservoirs is higher than that in conventional gas reservoirs [10].Therefore, it is essential to consider the influence of water saturation when studying the slip flow effect in tight gas reservoirs.Some scholars believe that as water saturation increases, the slip flow effect in tight sandstone reduces [11,12].However, other scholars have different conclusions, suggesting that the slip flow effect increases with water saturation [13,14].Unfortunately, due to limited research on this topic, the reasons behind this contradiction remain unknown.Hence, it is necessary to conduct a study on the slip flow effect characteristics in tight water-bearing gas reservoirs.
The study of stress sensitivity in reservoirs has a long history.Geertsma [15] defined the rock compressibility coefficient in 1957 to quantitatively describe the phenomenon of pore volume changes caused by changes in reservoir pore pressure.Stress sensitivity has a significant impact on the development of low-permeability oil and gas fields.Scholars have conducted extensive research on the stress sensitivity characteristics of low-permeability tight reservoirs and their impact on gas well productivity [16][17][18][19][20].As the water saturation of the reservoir gradually increases, the stress sensitivity of tight reservoirs becomes more significant [21][22][23][24].If there is bound water in the rock pores, it will have a certain impact on the strength of argillaceous reservoirs.Bound water forms a water film on the surface of rock particles, reducing pore volume and throat radius and leading to lower permeability.When the stress applied to the rock changes, causing deformation, the permeability of the rock changes more dramatically.Therefore, reservoirs saturated with water exhibit higher stress sensitivity.
The threshold pressure gradient (TPG) is the minimum pressure gradient required to establish a continuous flow of the non-wetting phase in the pores of a rock.The concept of the TPG was first proposed in 1951, and several foreign scholars have since proven through experimental studies the existence of the TPG during fluid seepage processes in low-permeability tight reservoirs [25,26].The TPG significantly influences gas well productivity.In previous studies, the TPG was often treated as a constant value.However, research has shown that the TPG is not fixed [27,28].Zafar (2020) conducted a sensitivity analysis of the two main influencing factors of the TPG, namely, permeability and water saturation.The results showed that the TPG is a power function of these two factors.It decreases with increasing permeability but increases with increasing water saturation [29].Through experimental research, Song (2015) concluded that the TPG in low-permeability tight reservoirs gradually decreases as permeability decreases [30].Zhu (2022), using experiments with tight sandstone core samples, found that higher water saturation gradually increases the TPG [31].Wang (2022)'s experimental results indicated that under water conditions, seepage in tight gas reservoirs exhibits nonlinear characteristics and forms the TPG.The TPG for the gas phase has a close power-law relationship with water saturation [32].For the water-producing gas wells in water-bearing-inclined gas reservoirs considering stress sensitivity, Fu (2022) has developed a new gas well production equation that accurately determines the relationship between gas well production and stress sensitivity and water production [33].
In summary, most scholars have primarily focused on the impact of single-phase gas or single-factor gas-water interactions on well productivity.However, during the production process of gas wells, the fracturing process leads to different flow patterns in various regions, which conventional productivity models have not taken into account.On the other hand, tight gas reservoirs generally contain water, but most researchers have not quantitatively characterized the impact of water saturation on the permeability mechanisms of gas wells.In this study, the influence of water saturation on permeability mechanisms was analyzed through experiments.Formulas describing the relationship between the stress sensitivity coefficient, initiation pressure gradient, permeability, and water saturation were derived.Through the analysis of experimental results, a definition for tight gas reservoirs with high water saturation was proposed.A multifactor complex flow condition model considering the influence of water saturation on permeability mechanisms was established for fractured horizontal wells in tight gas reservoirs.This model was then compared with the unobstructed flow rates obtained from productivity tests to validate its accuracy.The new model presented in this paper provides a basis for predicting the productivity of tight sandstone gas reservoirs.

Experimental Study on the Seepage Mechanism of Tight Cores
Dongsheng gas field is a tight sandstone gas reservoir located in the northeastern Ordos Basin, northern China.The main gas-bearing horizon is the Shihezi Formation and the Shanxi Formation with depths from 3000 to 3600 m.The predominant lithology of the reservoir rocks is lithic quartz sandstone and quartz sandstone, with a porosity ranging from 5.0% to 18.8% and permeability ranging from 0.20 to 3.99 mD.The liquid-gas ratio during production ranges from 5.7 to 7.4 m 3 /10 4 m 3 .This gas field is characterized as a lithologic, lithologic-structural, and structural-fracture composite gas reservoir.

Characterization of Gas Slip Effects in Tight Water-Bearing Gas Reservoirs
The gas slip effect refers to the gas in the low porosity porous media flow.When the average free range of gas molecules (gas molecules and other molecules between two successive collisions, through the average value of straight-line distance) is close to the radius of the pore, gas molecule and media pore wall collision increases and the wall of each molecule is in a state of motion.Therefore, in the wall of the tube, a slip flow is generated on the tube wall, resulting in a non-Darcy flow (non-linear flow) when measuring permeability with gas, making the measured gas permeability greater than the actual permeability of the phenomenon, as was discovered by Klinkenberg in 1941 in a test [34].
The corrected model for gas permeability measurements: where K is the gasometric permeability, K ∞ is the absolute permeability, b is the slip factor, and p m is the average pore pressure.Natural cores from the Dongsheng gas field were used in the experiments to study the gas slip effects of tight cores.Relevant basic petrophysical parameters of the samples are listed in Table 1.To conduct the experiments, the core samples were initially dried for 24 h and their dry weights were recorded.Subsequently, we saturated the samples using a combination of vacuum suction and pressure injection methods to calculate the pore volume and porosity.Finally, the samples were placed in a core holder, and high-purity nitrogen gas was used for displacement to achieve the desired experimental water saturation.The prepared rock samples with varying water saturations were individually placed into an inert gas permeability measurement system for low-permeability rocks.This allowed us to measure the permeability of the samples under different confinement pressures and injection pressures.To ensure the reliability of the data, this article conducted tests on four sets of samples for each water saturation level and obtained experimental results by measuring their slip factor.For specific data, please refer to Table 2.The data from Figure 1 (X1-X4) clearly demonstrate that as the water saturation of the rock samples increases, the slip factor gradually decreases.This indicates that with the increasing water saturation of the rock samples, the gas slip effect diminishes progressively, and at a water saturation of 40%, the gas slip effect is almost completely inhibited.

Stress-Sensitive Characterization of Tight Water-Bearing Gas Reservoirs
Rock stress sensitivity refers to the extent to which the permeability of rock changes under different effective stress conditions.There are noticeable differences in stress sensitivity among different types of reservoir rocks.Stress sensitivity is assessed by applying external stress from overlying strata and internal fluid pressure.Therefore, measuring the permeability of rocks under different confining pressure conditions to evaluate their stress sensitivity is one of the most commonly used methods in laboratory experiments.This method complies with the petroleum industry standard SY/T5358-2010, titled "Evaluation Method for Reservoir Sensitivity Flow Experiments".The conventional stress sensitivity evaluation experimental method involves altering the confining pressure of the core under atmospheric pressure conditions and subsequently measuring the permeability under the corresponding pressures.
This study utilized nine core samples from a tight gas reservoir in the Ordos Basin, and the basic rock properties of these samples are listed in Table 3.In the stress sensitivity experiments conducted under confined water conditions, the establishment of water saturation within the cores was crucial.Firstly, the cores were dried and subjected to vacuum treatment to simulate the saturation of formation water, and the mass of the saturated

Stress-Sensitive Characterization of Tight Water-Bearing Gas Reservoirs
Rock stress sensitivity refers to the extent to which the permeability of rock changes under different effective stress conditions.There are noticeable differences in stress sensitivity among different types of reservoir rocks.Stress sensitivity is assessed by applying external stress from overlying strata and internal fluid pressure.Therefore, measuring the permeability of rocks under different confining pressure conditions to evaluate their stress sensitivity is one of the most commonly used methods in laboratory experiments.This method complies with the petroleum industry standard SY/T5358-2010, titled "Evaluation Method for Reservoir Sensitivity Flow Experiments".The conventional stress sensitivity evaluation experimental method involves altering the confining pressure of the core under atmospheric pressure conditions and subsequently measuring the permeability under the corresponding pressures.
This study utilized nine core samples from a tight gas reservoir in the Ordos Basin, and the basic rock properties of these samples are listed in Table 3.In the stress sensitivity experiments conducted under confined water conditions, the establishment of water saturation within the cores was crucial.Firstly, the cores were dried and subjected to vacuum treatment to simulate the saturation of formation water, and the mass of the saturated cores was measured.Subsequently, the cores were placed in core holders, and gas displacement was conducted.During this operation, appropriate displacement pressures were selected based on the core's porosity and permeability data.Throughout the displacement process, the cores were taken out and rotated at both ends to ensure an even distribution of confined water.The mass of the cores was repeatedly measured until the desired level of confined water saturation was achieved.Subsequently, stress sensitivity experiments were conducted under these water saturation conditions.
The stress sensitivity coefficient (Ss) is calculated using the following formula: where σ i denotes the confining pressure, k i is the permeability of the core at different confining pressures, σ * is the reference stress point, and σ i is the permeability corresponding to σ * .The dimensionless net confining pressure and dimensionless permeability can be fitted into the following relationship: where K i represents the initial permeability.
Taking core X1 as an example, the experimental data indicate that it exhibits different stress sensitivities under various levels of confined water saturation.For each distinct level of confined water saturation, the stress sensitivity coefficient (S s ) of core X1 varies, and these coefficients can be normalized based on stress sensitivity experimental data.According to the results shown in Figure 2, it can be observed that the stress sensitivity coefficients (S s ) of the nine cores involved in the experiment change under different levels of confined water saturation and water saturation conditions.Further analysis reveals that, at zero water saturation, core X1 exhibits the lowest stress sensitivity, but as water saturation increases, its sensitivity gradually intensifies.As permeability stabilizes, an increase in confined water saturation leads to a significant exponential growth trend in the stress sensitivity coefficient of X1.
where Ki represents the initial permeability.
Taking core X1 as an example, the experimental data indicate that it exhibits different stress sensitivities under various levels of confined water saturation.For each distinct level of confined water saturation, the stress sensitivity coefficient (Ss) of core X1 varies, and these coefficients can be normalized based on stress sensitivity experimental data.According to the results shown in Figure 2, it can be observed that the stress sensitivity coefficients (Ss) of the nine cores involved in the experiment change under different levels of confined water saturation and water saturation conditions.Further analysis reveals that, at zero water saturation, core X1 exhibits the lowest stress sensitivity, but as water saturation increases, its sensitivity gradually intensifies.As permeability stabilizes, an increase in confined water saturation leads to a significant exponential growth trend in the stress sensitivity coefficient of X1.In the experimental process, in order to study the influence of permeability on stress sensitivity at a specific bound water saturation, the bound water saturation values established for each core sample were kept relatively consistent.The difference between the maximum and minimum values for each of the nine core samples at each bound water In the experimental process, in order to study the influence of permeability on stress sensitivity at a specific bound water saturation, the bound water saturation values established for each core sample were kept relatively consistent.The difference between the maximum and minimum values for each of the nine core samples at each bound water saturation point did not exceed 3%.Figure 3 illustrates the relationship between permeability and the stress sensitivity coefficient for the nine core samples at a specific bound water saturation.
Processes 2023, 11, x FOR PEER REVIEW 8 of 37 saturation point did not exceed 3%.Figure 3 illustrates the relationship between permeability and the stress sensitivity coefficient for the nine core samples at a specific bound water saturation.When the water saturation is at 0%, an exponential relationship is observed between the stress sensitivity coefficient and permeability.As permeability decreases, the stress sensitivity coefficient shows an increasing trend.Similarly, at water saturations of 45%, 55%, and 65%, there continues to be an exponential relationship between permeability and When the water saturation is at 0%, an exponential relationship is observed between the stress sensitivity coefficient and permeability.As permeability decreases, the stress sensitivity coefficient shows an increasing trend.Similarly, at water saturations of 45%, 55%, and 65%, there continues to be an exponential relationship between permeability and the stress sensitivity coefficient.Therefore, it can be reasonably assumed that the relationship between permeability and the stress sensitivity coefficient can be described by an exponential function: The relationships between the exponent coefficients a and b in Equation ( 2) and water saturation are arranged as shown in Figures 4 and 5.When the water saturation is at 0%, an exponential relationship is observed between the stress sensitivity coefficient and permeability.As permeability decreases, the stress sensitivity coefficient shows an increasing trend.Similarly, at water saturations of 45% 55%, and 65%, there continues to be an exponential relationship between permeability and the stress sensitivity coefficient.Therefore, it can be reasonably assumed that the relationship between permeability and the stress sensitivity coefficient can be described by an exponential function: The relationships between the exponent coefficients a and b in Equation ( 2) and water saturation are arranged as shown in Figures 4 and 5.The relationship between the multiplicative power coefficient a and the bound wate saturation Sw can be obtained from Figure 4: (5 The relationship between the multiplicative power index b and the bound water sat The relationship between the multiplicative power coefficient a and the bound water saturation S w can be obtained from Figure 4: a = 0.28421 + 0.000148694 * e (S w /0.09258) (5) The relationship between the multiplicative power index b and the bound water saturation S w can be obtained from Figure 5: b = −0.0317* e (S w /0.17024) + 0.31759 (6) By substituting Equations ( 3) and (4) into Equation (2), we obtain: S S = (0.28421 + 0.000148694 * e (S w /0.09258) ) * K −(−0.0317* e (Sw /0.17024) +0.31759) where S s is the stress sensitivity factor for bound water conditions; S w is the binding water saturation; and K is the absolute permeability, 10 −3 µm 2 .Formula (5) represents the stress sensitivity relationship of low-permeability tight gas reservoirs under confined water conditions.Since conducting stress sensitivity experiments under confined water conditions involves complex procedures and significant difficulties, the significance of Formula (5) lies in its ability to determine the reservoir's stress sensitivity under confined water saturation conditions using the permeability and water saturation values of the tight gas reservoir.In this case, this formula provides a convenient method, thereby avoiding the extensive use of manpower, resources, and financial investment required for experiments.Therefore, this approach can better assist researchers in understanding the mechanical properties of reservoir rocks and provide more reliable support for the development of tight gas reservoirs.

Characteristics of the TPG in Tight Water-Bearing Gas Reservoirs
The Jamin effect is one of the main reasons for the generation of the TPG in tight sandstone gas reservoirs.In porous media within tight reservoirs, there are pores of different sizes and narrow throats connecting these pores, forming a complex pore network.When bubbles flow from larger pores into narrower throats, they encounter resistance and need to overcome the deformation of the bubbles to continue flowing.This phenomenon is known as the Jamin effect.When the external driving pressure cannot overcome the capillary pressure, bubbles accumulate and block the throats.Over time, bubbles accumulate energy until gas and water phases can break through these constraints and start flowing.Therefore, a necessary condition for gas flow is that the pressure difference between the two sides of the bubble's surface reaches a certain level.This effect is crucial for understanding the productivity of gas reservoirs.
Five core samples from a tight sandstone gas reservoir in the Ordos Basin were selected, with the physical properties of the core samples as shown in Table 4. Experiments were conducted under different water saturation conditions to determine the TPG.The sentence describes the permeability characteristics under different water saturation conditions, as shown in Figure 6, revealing an important phenomenon: as the water saturation of the rock core samples increases, the TPG gradually increases.Especially when the permeability remains constant, with the increasing saturation of the bound water, the TPG exhibits a pronounced exponential growth trend.This discovery highlights the significant influence of water saturation on rock permeability and startup characteristics.

X3
2 The sentence describes the permeability characteristics under different water saturation conditions, as shown in Figure 6, revealing an important phenomenon: as the water saturation of the rock core samples increases, the TPG gradually increases.Especially when the permeability remains constant, with the increasing saturation of the bound water, the TPG exhibits a pronounced exponential growth trend.This discovery highlights the significant influence of water saturation on rock permeability and startup characteristics.To study the influence of permeability on stress sensitivity under a certain bound water saturation, in the experimental process, it is ensured that the bound water saturation values established for each core sample are essentially consistent.Furthermore, the difference between the maximum and minimum values at each water saturation point does not exceed 3%.As shown in Figure 7, the study investigated the relationship between permeability and initiation pressure gradient for these five core samples at different water saturation levels.To study the influence of permeability on stress sensitivity under a certain bound water saturation, in the experimental process, it is ensured that the bound water saturation values established for each core sample are essentially consistent.Furthermore, the difference between the maximum and minimum values at each water saturation point does not exceed 3%.As shown in Figure 7, the study investigated the relationship between permeability and initiation pressure gradient for these five core samples at different water saturation levels.According to the results of the pressure gradient testing experiment, it is evident that there is a linear relationship between the reciprocal of the core sample's permeability and the startup pressure gradient, while the coefficient "a" exhibits an exponential relationship with water saturation (Sw) (as detailed in Figure 8).Additionally, the coefficient "b" shows a linear relationship with the reciprocal of water saturation (as detailed in Figure 9).By fitting the test data from the five core samples, a relationship formula for the TPG in tight sandstone gas reservoirs was derived.This formula is related to both permeability and water saturation.According to the results of the pressure gradient testing experiment, it is evident that there is a linear relationship between the reciprocal of the core sample's permeability and the startup pressure gradient, while the coefficient "a" exhibits an exponential relationship with water saturation (S w ) (as detailed in Figure 8).Additionally, the coefficient "b" shows a linear relationship with the reciprocal of water saturation (as detailed in Figure 9).By fitting the test data from the five core samples, a relationship formula for the TPG in tight sandstone gas reservoirs was derived.This formula is related to both permeability and water saturation.
According to the results of the pressure gradient testing experiment, it is evident tha there is a linear relationship between the reciprocal of the core sample's permeability and the startup pressure gradient, while the coefficient "a" exhibits an exponential relationship with water saturation (Sw) (as detailed in Figure 8).Additionally, the coefficient "b" show a linear relationship with the reciprocal of water saturation (as detailed in Figure 9).By fitting the test data from the five core samples, a relationship formula for the TPG in tigh sandstone gas reservoirs was derived.This formula is related to both permeability and water saturation.

Two-Phase Seepage Characteristics in Tight Water-Bearing Gas Reservoirs
In order to study the gas-water co-permeability characteristics of low-permeability reservoirs, this article selected seven tight sandstone core samples for gas-driven wate experiments.The experimental procedure is illustrated in Figure 10 and includes the fol lowing steps: First, the core samples are placed in core holders, and appropriate displacement pres sures are determined based on the physical properties of the core and water permeability test data.Next, the gas pressure inside the intermediate container is stabilized at the se displacement pressure, and both the inlet and outlet of the core are opened to initiate th

Two-Phase Seepage Characteristics in Tight Water-Bearing Gas Reservoirs
In order to study the gas-water co-permeability characteristics of low-permeability reservoirs, this article selected seven tight sandstone core samples for gas-driven water experiments.The experimental procedure is illustrated in Figure 10 and includes the following steps: First, the core samples are placed in core holders, and appropriate displacement pressures are determined based on the physical properties of the core and water permeability test data.Next, the gas pressure inside the intermediate container is stabilized at the set displacement pressure, and both the inlet and outlet of the core are opened to initiate the gas displacement water experiment.During the experiment, the weight of water exiting the core and the gas flow rate are continuously measured using a precision balance and a gas mass flow meter, respectively.Simultaneously, a computer collects data on cumulative water and gas production at different time intervals.Finally, using pressure, water production, and gas production data at different time points, gas-water relative permeability curves under various water saturation conditions during gas displacement water experiments are calculated using software for processing.These curves reflect the characteristics of gas-water co-seepage and the seepage behavior of gas and water.Pressure gradients have a significant impact on the permeation ability of gas and water phases.With an increase in pressure gradient, the relative permeability of the gas phase exhibits a leftward trend.This implies that under the same water saturation conditions, higher displacement pressure gradients lead to a decrease in the relative permeability of the gas phase.This effect is more pronounced at lower water saturations and gradually diminishes as water saturation increases to around 75% or so.On the other hand, the relative permeability of the water phase increases with an increase in pressure gradient, especially under conditions with higher water saturation, where the impact is more significant.
Regardless of the pressure gradient conditions, as water saturation increases, the relative permeability of the gas phase in tight sandstone gas reservoirs decreases sharply.When the water saturation reaches 80% or higher, the gas phase permeability becomes almost negligible, resulting in minimal gas penetration.
These findings emphasize the significant influence of pressure gradients and water saturation on the gas-water co-permeation characteristics of tight sandstone gas reservoirs, especially in low-permeability reservoirs.This has important practical implications for the development and management of gas reservoirs.Pressure gradients have a significant impact on the permeation ability of gas and water phases.With an increase in pressure gradient, the relative permeability of the gas phase exhibits a leftward trend.This implies that under the same water saturation conditions, higher displacement pressure gradients lead to a decrease in the relative permeability of the gas phase.This effect is more pronounced at lower water saturations and gradually diminishes as water saturation increases to around 75% or so.On the other hand, the relative permeability of the water phase increases with an increase in pressure gradient, especially under conditions with higher water saturation, where the impact is more significant.
Regardless of the pressure gradient conditions, as water saturation increases, the relative permeability of the gas phase in tight sandstone gas reservoirs decreases sharply.When the water saturation reaches 80% or higher, the gas phase permeability becomes almost negligible, resulting in minimal gas penetration.
These findings emphasize the significant influence of pressure gradients and water saturation on the gas-water co-permeation characteristics of tight sandstone gas reservoirs, especially in low-permeability reservoirs.This has important practical implications for the development and management of gas reservoirs.

Derivation of Productivity Equation for Fractured Horizontal Wells in High-Water-Cut Tight Gas Reservoirs
In high-water-cut tight gas reservoirs, the gas flow patterns vary in different regions due to hydraulic fracturing, accompanied by different flow effects.This section first defines a high-water-cut tight gas reservoir through the analysis of experimental results and then derives the productivity equations for different regions.These research findings contribute to a better understanding and management of the development process in high-water-cut tight gas reservoirs, aiming to enhance reservoir productivity and efficiency.

Definition of High-Water-Cut Tight Gas Reservoir
(1) Gas slippage effect The analysis of Figure 12a leads to the following conclusion: as the water saturation increases, the slip factor gradually decreases.This is because as the average free path of the gas molecules approaches the capillary diameter, the collision frequency between the gas molecules and the pore walls of the porous media increases, leading to the occurrence of sidewall slip flow phenomena.This in turn leads to a continuous increase in gas velocity.Therefore, when gas is used as a permeability measurement medium, non-Darcy flow phenomena will inevitably occur, causing the measured permeability to be significantly higher than the true permeability.

Derivation of Productivity Equation for Fractured Horizontal Wells in High-Water-Cut Tight Gas Reservoirs
In high-water-cut tight gas reservoirs, the gas flow patterns vary in different regions due to hydraulic fracturing, accompanied by different flow effects.This section first defines a high-water-cut tight gas reservoir through the analysis of experimental results and then derives the productivity equations for different regions.These research findings contribute to a better understanding and management of the development process in highwater-cut tight gas reservoirs, aiming to enhance reservoir productivity and efficiency.

Definition of High-Water-Cut Tight Gas Reservoir
(1) Gas slippage effect The analysis of Figure 12a leads to the following conclusion: as the water saturation increases, the slip factor gradually decreases.This is because as the average free path of the gas molecules approaches the capillary diameter, the collision frequency between the gas molecules and the pore walls of the porous media increases, leading to the occurrence of sidewall slip flow phenomena.This in turn leads to a continuous increase in gas velocity.Therefore, when gas is used as a permeability measurement medium, non-Darcy flow phenomena will inevitably occur, causing the measured permeability to be significantly higher than the true permeability.In addition, when gas circulates in water-bearing rock samples, water molecules typically adhere to pore walls in the form of a water film or enter smaller pores due to the pressure exerted by the gas and capillaries.However, as gas molecules are constantly alternating within the porous rock system, they have difficulty colliding with solid pore walls and can only collide with water molecules.This limitation limits the occurrence of slip flow.As water saturation increases, more water molecules cover the pore walls, reducing the frequency of collisions between gas molecules and solid pore walls and further limiting slip flow.In particular, when water saturation reaches 40%, the gas slip effect is almost completely restricted.Therefore, when studying gas reservoirs with water saturations above 40%, the gas slip effect can be neglected.
The analysis from Figure 12b leads to the following conclusion: an increase in water saturation results in a decrease in the slip factor, with changes in pore structure being the primary cause of this phenomenon.When water enters the rock framework, interfacial tension forms a water film, reducing the radii of pores and throats, leading to a decline in the formation's permeability and increasing the stress sensitivity of the reservoir.The presence of water reduces the compressive strength of the rock framework, causing further compression of particles and pores.Under high water saturation conditions, water forms a film on the surface of quartz particles.When effective stress increases, the water film loses balance, leading to pressure dissolution and SiO 2 precipitation.Consequently, this reduces the radii of pores and throats as well as the permeability.
When the water saturation is less than 40%, the influence of water saturation on the stress sensitivity coefficient is minimal and can almost be neglected.However, when the water saturation exceeds 40%, the impact of water saturation on the stress sensitivity coefficient exhibits an exponential increase trend.Therefore, when studying gas reservoirs with water saturation exceeding 40%, it is crucial to consider the significance of water saturation's effect on stress sensitivity.
(3) Threshold pressure gradient Based on the analysis from Figure 12c, the following conclusions can be drawn: when the water saturation is less than 30%, no TPG occurs during gas penetration.However, when the water saturation exceeds 40%, the TPG significantly changes with the increase in water saturation.This phenomenon is due to the fact that at higher water saturations, gas in the reservoir cannot exist in a continuous phase but disperses into multiple small bubbles for flow.As these small bubbles pass through the throats, they induce the Jamin effect at each throat, causing capillary pressures to add up in the displacement direction.On a macroscopic scale, the higher the reservoir's water saturation, the larger the value of the TPG, as this additive effect is more likely to occur in throats with thicker water films.
In summary, once the water saturation in tight gas reservoirs exceeds 40%, there is a significant change in the seepage mechanisms.The gas slip flow effect no longer plays a role, and stress sensitivity effects and the TPG become significantly influenced by water saturation.Therefore, it is necessary to categorize tight sandstone gas reservoirs with water saturation exceeding 40% separately.These can be referred to as "High-Water-Cut Tight Gas Reservoir", aiming for a better understanding and description of their seepage characteristics and behaviors.

Zoning Productivity Physicochemical Mode of Fractured Horizontal Wells in High-Water-Cut Tight Gas Reservoirs under Complex Seepage Conditions
Based on the performance of fractured gas wells in different production stages and the characteristics of gas flow during each stage, the fluid flow region from the reservoir to the fractured horizontal well can be divided into three distinct zones (As shown in Figure 13).These three zones include Zone I: a high-speed non-Darcy flow region inside the fractures.In this zone, gas flows at high speed in the fractures with predominant features being highspeed non-Darcy effects and stress-sensitive effects within the fractures.It is noteworthy that, unlike in vertical fractured wells, the cross-sectional area of horizontal fractured wells is much larger.This leads to differences in fluid seepage compared to vertical fractures, such as radial flow near the wellbore due to the convergence effect after the fluid linearly enters the fractures.Zone II: an elliptical non-Darcy flow region caused by boundary leakage at the fracture boundaries.In this zone, leakage at the fracture boundaries results in an elliptical non-Darcy flow region, where stress-sensitive effects in the matrix and the TPG effect on gas flow dominate.Zone III: an elliptical low-speed non-Darcy seepage region formed by fluid inflow from the surrounding formations.In this zone, fluid inflow from the surrounding formations creates an elliptical low-speed non-Darcy seepage region.The main characteristics of this zone are the simultaneous presence of stress-sensitive effects and TPG effects.
to the fractured horizontal well can be divided into three distinct zones (As shown in Figure 13).These three zones include Zone I: a high-speed non-Darcy flow region inside the fractures.In this zone, gas flows at high speed in the fractures with predominant features being high-speed non-Darcy effects and stress-sensitive effects within the fractures.It is noteworthy that, unlike in vertical fractured wells, the cross-sectional area of horizontal fractured wells is much larger.This leads to differences in fluid seepage compared to vertical fractures, such as radial flow near the wellbore due to the convergence effect after the fluid linearly enters the fractures.Zone II: an elliptical non-Darcy flow region caused by boundary leakage at the fracture boundaries.In this zone, leakage at the fracture boundaries results in an elliptical non-Darcy flow region, where stress-sensitive effects in the matrix and the TPG effect on gas flow dominate.Zone III: an elliptical low-speed non-Darcy seepage region formed by fluid inflow from the surrounding formations.In this zone, fluid inflow from the surrounding formations creates an elliptical low-speed non-Darcy seepage region.The main characteristics of this zone are the simultaneous presence of stress-sensitive effects and TPG effects.
This zoning helps to better understand the gas flow characteristics in different regions of the gas wells and the impact of various seepage effects.This information is crucial for optimizing gas well production and reservoir management.When establishing seepage models, the reservoir and fractures are often regarded as two independent permeable systems that mutually influence each other.They are connected based on the principle of equal seepage rates and pressures between the reservoir and fractures.This connection condition not only effectively describes the interaction between the reservoir and fractures but also enhances the accuracy and reliability of the model.Therefore, during the model design, it is essential to consider the mutual influence between the reservoir and fractures to better simulate real-world situations.
Basic Assumptions of the Mathematical Model: 1. Perforated and fractured horizontal wells are located in a central position of an isotropic, isothermal, top and bottom-sealed gas reservoir.Only the heterogeneity of permeability is considered.2. These artificial fractures are evenly distributed vertically along the wellbore, mirroring the gas reservoir's height.They exhibit symmetrical distribution around the wellbore and possess infinite conductivity, situated at the midpoint of the gas reservoir.This zoning helps to better understand the gas flow characteristics in different regions of the gas wells and the impact of various seepage effects.This information is crucial for optimizing gas well production and reservoir management.
When establishing seepage models, the reservoir and fractures are often regarded as two independent permeable systems that mutually influence each other.They are connected based on the principle of equal seepage rates and pressures between the reservoir and fractures.This connection condition not only effectively describes the interaction between the reservoir and fractures but also enhances the accuracy and reliability of the model.Therefore, during the model design, it is essential to consider the mutual influence between the reservoir and fractures to better simulate real-world situations.
Basic Assumptions of the Mathematical Model: 1. Perforated and fractured horizontal wells are located in a central position of an isotropic, isothermal, top and bottom-sealed gas reservoir.Only the heterogeneity of permeability is considered.

2.
These artificial fractures are evenly distributed vertically along the wellbore, mirroring the gas reservoir's height.They exhibit symmetrical distribution around the wellbore and possess infinite conductivity, situated at the midpoint of the gas reservoir.

3.
Isothermal gas seepage occurs, following the high-speed non-Darcy flow law inside the artificially fractured fractures.Fluid flow between matrix and fracture systems adheres to the low-speed non-Darcy seepage law.

4.
Both the reservoir and gas are compressible.The influence of gravity and capillary forces is neglected.

5.
The steady-state seepage of homogeneous gas-water two-phase fluids is assumed.Temperature changes and special physicochemical phenomena are ignored.6.
Contamination on the fracture walls is not considered, and contamination in the near-wellbore zone and the reservoir is neglected.The effects of gravity and capillary forces are ignored, as well as the pressure drop in the wellbore.
Based on the given assumptions, the physical model of the fractured horizontal well is formulated as shown in Figure 14.This model is designed for tight sandstone gas reservoirs.In the model, the horizontal well has a length of L and a radius of r w .Additionally, the model includes N fractures with a half-length of L f , each fracture has a width of W f , the drainage radius is r e , and the gas layer thickness of h.The model is situated within the heterogeneous space O-XYZ, where the primary permeability is oriented along the X, Y, and Z directions.The main permeability aligned with the X, Y, and Z directions are denoted as K x , K y , and K z , respectively.K z points in the vertical direction, the angle between the horizontal wellbore and the primary permeability K x is θ, and the fractures are perpendicular to the wellbore with permeability of K fi .
3. Isothermal gas seepage occurs, following the high-speed non-Darcy flow law inside the artificially fractured fractures.Fluid flow between matrix and fracture systems adheres to the low-speed non-Darcy seepage law.4. Both the reservoir and gas are compressible.The influence of gravity and capillary forces is neglected.5.The steady-state seepage of homogeneous gas-water two-phase fluids is assumed.
Temperature changes and special physicochemical phenomena are ignored.6. Contamination on the fracture walls is not considered, and contamination in the nearwellbore zone and the reservoir is neglected.The effects of gravity and capillary forces are ignored, as well as the pressure drop in the wellbore.
Based on the given assumptions, the physical model of the fractured horizontal well is formulated as shown in Figure 14.This model is designed for tight sandstone gas reservoirs.In the model, the horizontal well has a length of L and a radius of rw.Additionally, the model includes N fractures with a half-length of Lf, each fracture has a width of Wf, the drainage radius is re, and the gas layer thickness of h.The model is situated within the heterogeneous space O-XYZ, where the primary permeability is oriented along the X, Y, and Z directions.The main permeability aligned with the X, Y, and Z directions are denoted as Kx, Ky, and Kz, respectively.Kz points in the vertical direction, the angle between the horizontal wellbore and the primary permeability Kx is θ, and the fractures are perpendicular to the wellbore with permeability of Kfi.

Transformation of Heterogeneity
In the heterogeneous space O-XYZ, coordinate transformation is used to equivalently transform a permeability heterogeneous gas reservoir into a homogeneous gas reservoir.

Transformation of Heterogeneity
In the heterogeneous space O-XYZ, coordinate transformation is used to equivalently transform a permeability heterogeneous gas reservoir into a homogeneous gas reservoir.
The following non-uniform coordinate transformation is performed: After coordinate transformation, the seepage problem in the heterogeneous space O-XYZ is transformed into a problem in the homogeneous space O-X 1 Y 1 Z 1 .At this point, in space O-X 1 Y 1 Z 1 , the average permeability is K 1 , the half-length of the fracture is L f1 , the width of the fracture is W f1 , the gas release radius is r el , the gas reservoir thickness is h 1 , and the wellbore radius is r wl .
After a coordinate transformation, the horizontal cross-section of the hydraulic fracturing wellbore in the O-XYZ space has transitioned from a circle with a radius of r w to an ellipse within the O-X 1 Y 1 Z 1 space, characterized by two semi-axes r w ( ε 2/ε 3 ) 1 /3 and r w ( ε 3/ε1) 1 /3 .When the fluids near the wellbore reach a state of equilibrium, the wellbore radius is equal to the arithmetic mean of the major and minor semi-axes.Therefore, the wellbore diameter of the hydraulic fracturing well after the coordinate transformation is:

Derivation of Zoning Productivity Calculation Method for Fractured Horizontal Wells
The production of tight high-water gas wells is divided into three zones: Zone I represents the high-speed non-Darcy flow in fractures, Zone II reflects the elliptical Darcy flow controlled by fracture boundaries, and Zone III represents the low-permeability non-Darcy flow in the matrix.To establish a production forecasting model for fractured gas wells, the concept of perturbed ellipses and the idea of equivalent developing rectangles are introduced.This model comprehensively considers the effects of gas slippage, TPG, and stress sensitivity on productivity, providing a theoretical basis for the rational determination of productivity in high-water-cut tight gas reservoirs.

Zone I Productivity Prediction Model
In the fracturing operation of horizontal gas wells, the fracturing stage is significant and crucial.During this crucial period, the permeability in the fractures is usually at a relatively high level, enabling gas to flow within the fractures at high speeds through non-Darcy flow.Additionally, due to the significantly larger cross-sectional area of the transverse fractures in fractured horizontal wells compared to the cross-sectional area of the wellbore, the fluid flow process differs from that in vertical wells with vertical fractures.After linear flow into the fractures, fluid flow near the wellbore transitions to radial flow due to the confluence effect.In this complex process, the high-speed non-Darcy effect and stress-sensitive effect are very significant.It is worth noting that the characteristics of this region allow us to temporarily exclude the effects of the TPG and slippage.
In the fluid movement within the fracture, from the tip of the fracture to the point where the linear and radial flows intersect, the flow conditions can be equivalent to the flow in a strip-shaped formation with a length of L f1 , a width of W f1 , and a height of h 1 .On this basis, assuming that the radius at the intersection of linear and radial flows within the fracture is r x , the corresponding pressure is p x .
The equations governing the movement of gas and water phases within the fracture are: Gas phase: where r is the linear seepage distance along the centerline of the fracture (m); K f is the absolute permeability of the fracture (mD); β f is the Forchheimer coefficient (m −1 ); v g is the gas velocity within the fracture (m/d); v w is the water phase velocity within the fracture (m/d); and K rg is the relative permeability of the gas phase under two-phase gas-water seepage conditions within the fracture.K rw is the relative permeability of the water phase under two-phase gas-water seepage conditions within the fracture.w f1 is the width of the fracture (m).h 1 is the height of the fracture (m). Where , under conditions where different proppants are used, m and n values are as shown in Table 6.Considering the stress-sensitive effect of water saturation: 000148694 * e (S w /0.09258) b = −0.0317* e (S w /0.17024) + 0.31759 (18) where K f represents the initial permeability of the fracture (mD); a is the coefficient parameter for the stress sensitivity of water saturation; and b is the exponent parameter for the stress sensitivity of water saturation.
Considering the stress sensitivity and turbulent effects during the linear flow stage of gas in the fracture, the flow equation for hydraulic fracturing horizontal wells can be derived as follows: Gas phase: The following equations can be obtained by combining Equations ( 17) and ( 18): The definition of two-phase pseudo-pressure within the fracture is as follows: Integrating the equation resulting from the combination of Equations ( 19) and ( 20): where p f is the pressure at the fracture tip (MPa); p x is the pressure at the boundary between the linear seepage field and the radial seepage field within the fracture (MPa); r = h/2; r x is the radial radius at the boundary between the linear seepage field and the radial seepage field within the fracture (m). (

2) Radial flow stage
When the fluid in the fracture starts to undergo radial flow, it is equivalent to flowing from the boundary to the wellbore with radius r w1 in a circular formation with supply radius r x , supply boundary pressure p x , and formation thickness W f1 .Based on the previous derivation, the following seepage equations of motion can be obtained for the gas and water phases in a single fracture of a fractured horizontal well during the radial flow stage in the fracture, respectively, under the consideration of the high-velocity non-Darcy and stress-sensitive effects.
Gas phase: The following equations can be obtained by combining Equations ( 22) and ( 23): (5.787 ρ w q w +ρ gsc q gsc πrw f 1 The pseudo-pressure for the two phases within the fracture is defined.Combining Equations ( 20) and ( 22), and integrating from the boundary of the radial seepage field within the fracture (radius r x ) to the wellbore, the productivity equation for any single fracture in the radial seepage stage of the fracture, considering stress sensitivity and high-speed non-Darcy effects, can be obtained under steady-state gas-water flow conditions.
Due to the principle of water and electricity similitude, the fluid flow rates in the two seepage fields are equal.Therefore, by simultaneously considering Equations ( 21) and ( 25), a productivity prediction model for a single fracture in hydraulic fracturing horizontal wells in tight sandstone gas reservoirs is derived.This model takes into account permeability heterogeneity, high-speed non-Darcy flow, gas slippage phenomenon, and stress sensitivity.

.2. Region II Productivity Prediction Model
In Region II of horizontal gas wells, the fracture boundaries induce a planar twodimensional elliptical seepage phenomenon within the formation (As shown in Figure 15).In this scenario, gas seepage is influenced by the TPG, and stress-sensitive effects become apparent.It is worth mentioning that the influence of slip flow can be neglected.Gas seepage under these conditions conforms to the characteristics of low-speed non-Darcy elliptical seepage.

Region II Productivity Prediction Model
In Region II of horizontal gas wells, the fracture boundaries induce a planar twodimensional elliptical seepage phenomenon within the formation (As shown in Figure 15).In this scenario, gas seepage is influenced by the TPG, and stress-sensitive effects become apparent.It is worth mentioning that the influence of slip flow can be neglected.Gas seepage under these conditions conforms to the characteristics of low-speed non-Darcy elliptical seepage.In the region around the horizontally fractured well, the external seepage field shows a series of elliptical curves (see Figure 14).When the formation is in the non-Darcy flow zone, the seepage form is elliptical seepage, and the elliptical seepage law can be applied to describe its seepage characteristics.There is the following relationship between elliptical coordinates and right-angle coordinates: In the region around the horizontally fractured well, the external seepage field shows a series of elliptical curves (see Figure 14).When the formation is in the non-Darcy flow zone, the seepage form is elliptical seepage, and the elliptical seepage law can be applied to describe its seepage characteristics.There is the following relationship between elliptical coordinates and right-angle coordinates: where a is the half-axis length of the elliptical seepage field, m; b is the short half-axis length of the elliptical seepage field, m.
The seepage field of elliptical seepage is usually described by a developing family of rectangles, so the average short half-axis is: The pressure gradient in the elliptical coordinate system is: In gas-water two-phase percolation, the continuity equations are as follows: Gas phase: Water phase: aK wi e bp K rw v w = 2500 ρ w q w ρ w x f cosh ξ (33) where p is the formation pressure, MPa; r is the radial radius, m; m r is gas mass flow at radius r, kg/s; m w is the gas mass flow at radius r, kg/s; µ g is the viscosity of gas, MPa•s; µ g is the viscosity of water, MPa•s; K rg is the gas phase relative permeability, 10 −3 µm 3 ; and K rw is the relative permeability of water phase, 10 −3 µm 3 .
Considering the influence of water saturation on stress sensitivity and TPG, the permeability is revised as follows: Gas phase: The following equations can be obtained by combining Equations ( 32) and ( 33): The pseudo-pressure of two-phase formation is defined as: The pseudo-startup pressure of two phases is defined as: Integrating the equation resulting from the combination of Equations ( 34)-(36): where a is the semi-major axis length of the hydraulic fracture flow ellipse (m); ξ corresponds to an ellipse with a long axis of r el and ξ f corresponds to an ellipse with an average short axis of W f1 /2; p f is the fracture tip pressure, MPa, and ξ is an elliptic coordinate; ξ e is the coordinates of the gas reservoir supply boundary in the elliptic coordinate system, dimensionless; and ξ f is the coordinates of the fracture tip in the elliptic coordinate system, dimensionless.

Region III Productivity Prediction Model
In Region III, due to the continuous depletion of reservoir energy, fractures deform, resulting in relatively low permeability.During this stage, gas production mainly relies on the release of gas from the formation matrix in the far-field area under low-pressure conditions.It is important to emphasize that this study focuses on tight water-bearing gas reservoirs.Therefore, the effect of slippage is no longer considered.Moreover, during this process, stress sensitivity gradually decreases, the effect of the TPG impacts gas seepage, and it complies with the laws of low-speed non-Darcy flow.
Referring to the seepage equations in Section 3.4.2, the gas-water biphasic productivity equation for a single fracture in hydraulic fracturing horizontal wells, considering the TPG and stress-sensitive effects under stable seepage conditions in the formation, can be derived as follows: where r e is the supply radius (m); p ξ is the pressure at the interface (MPa); and r ξ is the distance from the well axis at the interface (m).

Establishment of Hydraulic Fracturing Horizontal Well Productivity Equation (1) Single Fracture Productivity Equation Development
The production of fluids from the formation results from the series connection of the seepage field within the formation and the linear seepage field within the fractures.When all three flow regions are simultaneously seeping, according to the principle of mass conservation, the seepage rates and pressures at the interfaces between adjacent flow regions should be equal.This allows the elimination of the interface pressures.By simultaneously solving Equations ( 26), (37) and (38), a hydraulic fracturing horizontal well prediction model is obtained, considering factors such as non-Darcy flow at low and high speeds, stress sensitivity, and TPG.This model can be used to predict the productivity of tight, highly saturated gas reservoirs.
(2) Equivalent wellbore radius model Under steady-state seepage conditions, if a productivity calculation model for a certain type of gas well, influenced by complex factors, is known, it can be compared with a general vertical well productivity calculation model that does not consider those complex factors under the general Darcy conditions.When the production rates of the two models are equal, the calculated equivalent radius of the general vertical well is referred to as the equivalent wellbore radius, denoted as r equ , and is expressed as follows: (3) Establishment of productivity equation for N fractures in hydraulic fracturing horizontal wells.
Building on any single fracture productivity calculation model derived from the previous text, employing the equivalent wellbore diameter theory, the single fracture is transformed into a regular vertical well.Accounting for inter-fracture interference, according to the pressure superposition principle, it can be expressed as: By solving the equations simultaneously, the production of N fractures in hydraulic fracturing horizontal wells can be obtained.The total production of hydraulic fracturing horizontal wells is: q gsc = q gsc1 + q gsc2 + q gsc3 + . . .q gscn + . . .q gscN (44)

Example Calculation
A horizontal hydraulic fracturing well in a tight, highly saturated gas reservoir undergoes a large-scale proppant fracturing transformation.The proppant used for hydraulic fracturing is a mixture of 70/140 mesh quartz sand and 40/70 mesh coated sand.The original formation pressure is 50 MPa, and the drainage radius of the reservoir is 200 m.The bottom hole flowing pressure is 25 MPa, the formation salinity is 1.06, and the reservoir thickness is 20 m.The density of water is 1 g/cm 3 , with a reservoir water saturation of 65%.The viscosity of water is 1 mPa•s, and the formation temperature is 340 K.The initial formation pressure is 40 MPa, and the fracture width is 0.008 m.The effective half-length of the fracture is 106.2 m, and the absolute permeability of the fracture is 50 Darcy.
Using the model described in this paper and existing productivity calculation models, the productivity of the gas well is calculated separately, and the results are shown in the graph.
Figure 16 shows that the gas well productivity calculated using the zonal productivity model is significantly higher than the results obtained from existing models.This is because, compared to existing models, the model in this paper considers the seepage characteristics of the fracturing zone.After the fracturing transformation, the reservoir properties and seepage capacity in the transformed area significantly improved.As a result, the gas well productivity increases significantly compared to the traditional models that do not consider zonal divisions.

Sensitivity Analysis
Based on the validation from the practical example, a comparative analysis was conducted by adjusting different parameters in the model.This study aimed to investigate the impact of various factors on the zonal productivity calculation results.

The Impact of Differential Stress Sensitivity in Different Regions
Initially, by adjusting the stress sensitivity coefficients inside and outside the fracturing zone, the study explored the impact of differential stress sensitivity in different regions on the productivity of hydraulic fracturing gas wells in tight gas reservoirs.The results are shown in the figure .The analysis of Figure 17 reveals that, considering zonal stress sensitivity, the productivity of the gas well significantly decreases compared to the scenario where stress sensitivity is not considered.This decline is attributed to the continuous reduction in reservoir permeability as the development process progresses, significantly affecting the seepage capacity of fluids.Consequently, the productivity of the gas well experiences a significant decrease.

Sensitivity Analysis
Based on the validation from the practical example, a comparative analysis was conducted by adjusting different parameters in the model.This study aimed to investigate the impact of various factors on the zonal productivity calculation results.

The Impact of Differential Stress Sensitivity in Different Regions
Initially, by adjusting the stress sensitivity coefficients inside and outside the fracturing zone, the study explored the impact of differential stress sensitivity in different regions on the productivity of hydraulic fracturing gas wells in tight gas reservoirs.The results are shown in the figure .The analysis of Figure 17 reveals that, considering zonal stress sensitivity, the productivity of the gas well significantly decreases compared to the scenario where stress sensitivity is not considered.This decline is attributed to the continuous reduction in reservoir permeability as the development process progresses, significantly affecting the seepage capacity of fluids.Consequently, the productivity of the gas well experiences a significant decrease.The results in Figure 17 also indicate that, when considering stress sensitivity in a single region alone, the gas well productivity is lower compared to the scenario where stress sensitivity is not considered at all.Moreover, when stress sensitivity in Region I alone is considered, the gas well productivity is lower than when stress sensitivity in Regions II and III is taken into account.This is because although Region I has a higher sensitivity to stress and a more significant impact on seepage due to the presence of fractures its area is much smaller than that of Regions II and III.Therefore, when stress sensitivity in Region I alone is considered, the productivity is lower than when stress sensitivity in Regions II and III is considered.In other words, only when the stress sensitivity characteristics of all three regions are considered simultaneously, the gas well productivity can be accurately determined.

The Impact of Dynamic TPG in Tight Reservoirs
The impact of dynamic TPG on the productivity of hydraulic fracturing gas wells in tight gas reservoirs is shown in Figure 18.
From Figure 18, it can be observed that considering the TPG results in a significant decrease in gas well productivity compared to not considering the TPG.This is because when the TPG is taken into account, the additional resistance to fluid flow increases, leading to a decrease in productivity.Furthermore, when considering the dynamic TPG, gas well productivity further decreases compared to the conventional fixed TPG.This is primarily because when accounting for the dynamic TPG, the TPG in the reservoir gradually increases as the development progresses.The lower the reservoir pressure, the larger the TPG, leading to a more pronounced decrease in productivity.The results in Figure 17 also indicate that, when considering stress sensitivity in a single region alone, the gas well productivity is lower compared to the scenario where stress sensitivity is not considered at all.Moreover, when stress sensitivity in Region I alone is considered, the gas well productivity is lower than when stress sensitivity in Regions II and III is taken into account.This is because although Region I has a higher sensitivity to stress and a more significant impact on seepage due to the presence of fractures, its area is much smaller than that of Regions II and III.Therefore, when stress sensitivity in Region I alone is considered, the productivity is lower than when stress sensitivity in Regions II and III is considered.In other words, only when the stress sensitivity characteristics of all three regions are considered simultaneously, the gas well productivity can be accurately determined.

The Impact of Dynamic TPG in Tight Reservoirs
The impact of dynamic TPG on the productivity of hydraulic fracturing gas wells in tight gas reservoirs is shown in Figure 18.
From Figure 18, it can be observed that considering the TPG results in a significant decrease in gas well productivity compared to not considering the TPG.This is because when the TPG is taken into account, the additional resistance to fluid flow increases, leading to a decrease in productivity.Furthermore, when considering the dynamic TPG, gas well productivity further decreases compared to the conventional fixed TPG.This is primarily because when accounting for the dynamic TPG, the TPG in the reservoir gradually increases as the development progresses.The lower the reservoir pressure, the larger the TPG, leading to a more pronounced decrease in productivity.

The Impact of the Number of Fractures on Productivity
Analyzing Figure 19 shows that as the number of fractures increases, the IPR (Inflow Performance Relationship) curve shifts towards higher productivity.The gas well productivity increases, but the rate of increase becomes smaller.Figure 20 shows the variation in gas well productivity with the number of fractures.For the same hydraulic fracturing gas well, within a certain range, having more fractures leads to a larger effective fracture area, which is beneficial for increasing productivity.However, if the number of fractures continues to increase, the interference between fractures intensifies, leading to mutual inhibition among the fractures.As a result, the increase in gas well productivity becomes smaller, indicating that there exists an optimal number of fractures.Therefore, in practical production, it is crucial to optimize the number of fractures reasonably and avoid productivity losses caused by inter-fracture interference.

The Impact of the Number of Fractures on Productivity
Analyzing Figure 19 shows that as the number of fractures increases, the IPR (Inflow Performance Relationship) curve shifts towards higher productivity.The gas well productivity increases, but the rate of increase becomes smaller.Figure 20 shows the variation in gas well productivity with the number of fractures.For the same hydraulic fracturing gas well, within a certain range, having more fractures leads to a larger effective fracture area, which is beneficial for increasing productivity.However, if the number of fractures continues to increase, the interference between fractures intensifies, leading to mutual inhibition among the fractures.As a result, the increase in gas well productivity becomes smaller, indicating that there exists an optimal number of fractures.Therefore, in practical production, it is crucial to optimize the number of fractures reasonably and avoid productivity losses caused by inter-fracture interference.

The Impact of the Number of Fractures on Productivity
Analyzing Figure 19 shows that as the number of fractures increases, the IPR (Inflow Performance Relationship) curve shifts towards higher productivity.The gas well productivity increases, but the rate of increase becomes smaller.Figure 20 shows the variation in gas well productivity with the number of fractures.For the same hydraulic fracturing gas well, within a certain range, having more fractures leads to a larger effective fracture area, which is beneficial for increasing productivity.However, if the number of fractures continues to increase, the interference between fractures intensifies, leading to mutual inhibition among the fractures.As a result, the increase in gas well productivity becomes smaller, indicating that there exists an optimal number of fractures.Therefore, in practical production, it is crucial to optimize the number of fractures reasonably and avoid productivity losses caused by inter-fracture interference.

The Impact of Water-Gas Ratio (WRG) on Gas Well Productivity
As shown in the Figure 21, different IPR curves are plotted under various water-gas ratios (WRGs): one without considering the stress sensitivity and threshold pressure gradient (TPGASS), one considering a fixed initiation pressure gradient and stress sensitivity, one considering the influence of water saturation on the TPG and stress sensitivity, and one showing the unobstructed flow rate variation curve with changes in the water-gas ratio considering the impact of water saturation on TPG and stress sensitivity.The analysis indicates that gas well productivity decreases after water production, and as the watergas ratio increases, the curve shifts to the left.After water production, the water phase occupies the large pores initially occupied by the gas phase.More and more gas becomes discontinuous, reducing the relative permeability of the gas phase, which affects gas well productivity.With the increase in the water-gas ratio, the TPG and stress sensitivity effects increase, and the shift of the curve to the right becomes smaller.The lower the bottom hole flowing pressure, the greater the impact of TPG and stress sensitivity effects on gas well productivity.The variation curve of the production capacity with water-gas ratio is shown in Figure 22.Therefore, in productivity modeling, the influence of water saturation on the TPG and stress sensitivity effects should be considered.In practical production processes, appropriate operating procedures should be adopted to avoid early water breakthroughs or delay the time when the gas well starts producing water.

The Impact of Water-Gas Ratio (WRG) on Gas Well Productivity
As shown in the Figure 21, different IPR curves are plotted under various water-gas ratios (WRGs): one without considering the stress sensitivity and threshold pressure gradient (TPGASS), one considering a fixed initiation pressure gradient and stress sensitivity, one considering the influence of water saturation on the TPG and stress sensitivity, and one showing the unobstructed flow rate variation curve with changes in the water-gas ratio considering the impact of water saturation on TPG and stress sensitivity.The analysis indicates that gas well productivity decreases after water production, and as the watergas ratio increases, the curve shifts to the left.After water production, the water phase occupies the large pores initially occupied by the gas phase.More and more gas becomes discontinuous, reducing the relative permeability of the gas phase, which affects gas well productivity.With the increase in the water-gas ratio, the TPG and stress sensitivity effects increase, and the shift of the curve to the right becomes smaller.The lower the bottom hole flowing pressure, the greater the impact of TPG and stress sensitivity effects on gas well productivity.The variation curve of the production capacity with water-gas ratio is shown in Figure 22.Therefore, in productivity modeling, the influence of water saturation on the TPG and stress sensitivity effects should be considered.In practical production processes, appropriate operating procedures should be adopted to avoid early water breakthroughs or delay the time when the gas well starts producing water.

Summary
This study delves into the intricate permeation mechanisms of dense, highly watersaturated gas reservoirs through comprehensive experimentation.As the water saturation increases, the slip factor of gas gradually diminishes, becoming negligible at a water satu ration of 40%.This suggests that, under conditions of elevated water saturation, the gas flow tends towards Darcy flow, no longer exhibiting characteristics of non-Darcy flow.The dense gas reservoirs exhibit pronounced stress sensitivity, with the stress sensitivity coefficient escalating with rising water saturation, particularly becoming more pronounced beyond 40% water saturation.This observation underscores that, under conditions of high water saturation, the mechanical properties of rocks are more susceptible to stress influences, necessitating a deeper exploration and analysis to comprehend and effectively manage this stress sensitivity.Upon exceeding a water saturation of 30%, dense water-saturated gas reservoirs manifest a TPG, intensifying with the elevation of water saturation.This trend becomes notably significant, particularly beyond a water saturation of 40%.This result highlights the challenges in initiating gas flow within reservoirs under

Summary
This study delves into the intricate permeation mechanisms of dense, highly watersaturated gas reservoirs through comprehensive experimentation.As the water saturation increases, the slip factor of gas gradually diminishes, becoming negligible at a water satu ration of 40%.This suggests that, under conditions of elevated water saturation, the gas flow tends towards Darcy flow, no longer exhibiting characteristics of non-Darcy flow.The dense gas reservoirs exhibit pronounced stress sensitivity, with the stress sensitivity coefficient escalating with rising water saturation, particularly becoming more pronounced beyond 40% water saturation.This observation underscores that, under conditions of high water saturation, the mechanical properties of rocks are more susceptible to stress influences, necessitating a deeper exploration and analysis to comprehend and effectively manage this stress sensitivity.Upon exceeding a water saturation of 30%, dense water-saturated gas reservoirs manifest a TPG, intensifying with the elevation of water saturation.This trend becomes notably significant, particularly beyond a water saturation of 40%.This result highlights the challenges in initiating gas flow within reservoirs under high water saturation conditions, emphasizing the need to consider this factor for more accurate predictions of reservoir productivity.With the increase in water saturation, the relative permeability of gas in tight sandstone gas reservoirs sharply declines, reaching close to zero when water saturation exceeds 80%.This implies that under conditions of high water saturation, gas infiltration becomes nearly impossible, holding significant practical implications for the production and development of gas reservoirs.Addressing this challenge requires the implementation of corresponding technological measures.
Through a comprehensive analysis and summary of the experimental results, this study provided a clear definition of tight water-bearing gas reservoirs and delved into the different regional flow patterns.Given the heterogeneity and unique flow characteristics of tight sandstone gas reservoirs, various theories and methods were employed, including coordinate transformation, the elliptical coordinate theory, the perturbation ellipse theory, the pressure superposition principle, and the water-gas similarity principle, to establish a zoning productivity prediction model.This model takes into account multiple factors affecting productivity, including heterogeneity, the impact of gas well water production on productivity, and the different flow characteristics exhibited by gas wells in different regions.Using the newly developed method for predicting the production capacity of fractured wells in tight gas reservoirs, this study conducted analyses on the effects of regional differences in stress sensitivity, dynamic TPG, the number of fractures, and the water-gas ratio on well productivity.
In Region I, the high-speed non-Darcy flow within the fractures and the effects of stress sensitivity were thoroughly taken into account.This region is crucial for reservoir production since the high permeability of the fractures significantly impacts productivity.In Region II, particular attention was given to the stress sensitivity effects due to elliptical non-Darcy flow caused by fracture boundary drainage and the initiation pressure gradient effects experienced by the gas flow.The characteristic of this region lies in the significant influence of matrix stress sensitivity on productivity.Finally, in Region III, the effects of reservoir stress sensitivity and the TPG experienced by gas flow were comprehensively considered.This region requires special attention to reservoir characteristics because reservoir stress sensitivity significantly affects reservoir productivity.
This article established a mathematical model for the two-phase flow of gas and water in tight gas reservoirs.Through specific calculation examples, the practicality and accuracy of the model were verified, providing an important theoretical foundation and engineering guidance for the production capacity prediction of unconventional oil and gas resources.

Figure 1 .
Figure 1.The curves depicting the variation in slip factors in wells (X1-X4) with changing water saturation.

Figure 2 .
Figure 2. Relationship between water saturation and the stress sensitivity coefficient.

Figure 2 .
Figure 2. Relationship between water saturation and the stress sensitivity coefficient.

Figure 3 .
Figure 3.The influence of permeability on the stress sensitivity coefficient.

Figure 3 .
Figure 3.The influence of permeability on the stress sensitivity coefficient.

Figure 3 .
Figure 3.The influence of permeability on the stress sensitivity coefficient.

Figure 4 .
Figure 4. Relation between water saturation Sw and a.

Figure 4 . 3 Figure 5 .
Figure 4. Relation between water saturation S w and a. Processes 2023, 11, x FOR PEER REVIEW 9 of 3

Figure 5 .
Figure 5. Relation between water saturation S w and b.

Processes 2023 , 37 Figure 7 .
Figure 7.The impact of permeability on the TPG

Figure 7 .
Figure 7.The impact of permeability on the TPG.

8 Figure 8 .
Figure 8. Curve of water saturation versus TPG coefficient a.

Figure 8 .
Figure 8. Curve of water saturation versus TPG coefficient a.

Figure 9 .
Figure 9. Curve of water saturation versus TPG coefficient b.

Figure 9 .
Figure 9. Curve of water saturation versus TPG coefficient b.

Figure 12 .
Figure 12.The permeability law of tight water-bearing gas reservoirs.(a) The permeability law of tight water-bearing gas reservoirs; (b) the stress sensitivity characteristics of tight water-bearing gas reservoirs; (c) the TPG characteristics of the startup pressure gradient.

Figure 15 .
Figure 15.Conceptual diagram of the elliptical flow field within the formation.

Figure 15 .
Figure 15.Conceptual diagram of the elliptical flow field within the formation.

Figure 16 .
Figure 16.Comparison of results from different productivity calculation methods.

Figure 16 .
Figure 16.Comparison of results from different productivity calculation methods.

Processes 2023 , 37 Figure 17 .
Figure 17.The impact of zonal stress sensitivity on the productivity of hydraulic fracturing gas wells in tight gas reservoirs.

Figure 17 .
Figure 17.The impact of zonal stress sensitivity on the productivity of hydraulic fracturing gas wells in tight gas reservoirs.

Figure 18 .
Figure 18.Effect of dynamic threshold pressure gradient on fractured gas well productivity in tight gas reservoirs.

Figure 19 .
Figure19.The impact of the number of fractures on gas well productivity.

Figure 18 .
Figure 18.Effect of dynamic threshold pressure gradient on fractured gas well productivity in tight gas reservoirs.

Processes 2023 , 37 Figure 18 .
Figure 18.Effect of dynamic threshold pressure gradient on fractured gas well productivity in tight gas reservoirs.

Figure 19 .
Figure 19.The impact of the number of fractures on gas well productivity.Figure 19.The impact of the number of fractures on gas well productivity.

Figure 19 .
Figure 19.The impact of the number of fractures on gas well productivity.Figure 19.The impact of the number of fractures on gas well productivity.

Figure 20 .
Figure 20.The curve depicting the variation in gas well productivity with the number of fractures.

Figure 20 .
Figure 20.The curve depicting the variation in gas well productivity with the number of fractures.

Figure 22 .
Figure 22.Unobstructed flow rate variation curve with changes in WRG considering the impact of water saturation on the TPG and stress sensitivity effects.

Figure 22 .
Figure 22.Unobstructed flow rate variation curve with changes in WRG considering the impact of water saturation on the TPG and stress sensitivity effects.

Table 1 .
Petrophysical parameters of the core.

Table 2 .
Slip factor for different perimeter pressures; different water saturations for different core conditions.

Table 3 .
The results of stress sensitivity experiments at different water saturation levels.

Table 4 .
Basic rock properties for the TPG testing of the core samples.
Table 5 summarizes the experimental results data.

Table 5 .
Gas-water phase permeation experiment results.

Table 5 .
Gas-water phase permeation experiment results.

Table 6 .
Under different proppants, m and n value table.