{'fraction_non_alphanumeric': 0.02197802197802198, 'fraction_numerical': 0.011545886545886546, 'mean_word_length': 5.077842960288809, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Study on the Water Invasion and Its Effect on the Production from Multilayer Unconsolidated Sandstone Gas Reservoirs

Water invasion is a common occurrence in multilayer unconsolidated gas reservoirs, which results in excessive water production and reduces the economic life of gas wells. However, due to multiple layers, active edge water, and strong heterogeneity, the mechanism of water invasion and its effect in the unconsolidated sandstone gas reservoir require understanding in order to improve efficiency and minimize economic cost. In this study, an experimental study on edge water invasion of the multilayer commingled production in unconsolidated sandstone gas reservoirs was conducted to understand the water invasion process along with different permeability layers. The results show that the edge water invasion in the commingling production is mainly affected by two major factors including reservoir permeability and gas production rate, which jointly control the encroaching water advance path and speed. The nonuniform invade of edge water may occur easily and water prefers to invade toward the gas well along with high permeability layers when the commingling production is in the condition of large permeability gradient and high production rate. The bypass flow will occur when there are high permeability channels between the layers, which causes water blocking to low-permeability layers and periphery reservoirs far away from gas wells. The findings of this study can help for a better understanding of water invasion and the effects of reservoir properties so as to optimize extraction conditions and predict gas productivity in unconsolidated sandstone gas reservoirs.


Introduction
As a clean and efficient low-carbon energy, natural gas has been attracting greatly wide concerns from governments around the world, which exhibits an increasingly broad development prospect [1]. The natural gas resources in China are rather abundant, and the gas reserves and production are increasing rapidly. According to the statistics from PetroChina Company Limited, China's natural gas geological resources have exceeded 220 trillion cubic meters in 2015, and the recoverable gas resources were about 31% of total resources, and the annual production increased from 17.4 billion cubic meters in 1995 to 135 billion cubic meters in 2015 with an average annual growth rate of 11% [2][3][4]. The Quaternary unconsolidated sandstone gas reservoir in the Qaidam Basin, Northwest China, is an important strategic successor gas source for the West-East Gas Pipeline project, which belongs to a rare biogenetic gas reservoir [5,6]. Geologically, the gas reservoir is characterized by long interval, multiple layers, and active edge water, which is based on the multilayer commingled production, and there exist lots of difficult problems to solve in the production and development [7]. Therefore, it is very important to develop the gas reservoir effectively and guarantee a stable gas supply in the sustainable development of the national economy.
Most of the gas reservoirs in China belong to water-drive gas reservoirs surrounded by aquifers, in which there are about 40 to 50 percent of gas reservoirs with active edgeand-bottom water [8,9]. Due to the existence of edge-andbottom water in gas reservoirs, the reservoir pressure will gradually decrease with the exploitation of gas reservoirs, which results in the invasion of the edge-and-bottom water into the gas pay zones. The gas-water two-phase flow occurs in gas reservoirs with the edge-and-bottom water drive, and the gas phase permeability will be decreased, which consequently interferes with gas production and impacts on the recovery factor of gas reservoirs [10][11][12]. In addition, water invasion can also block a large amount of gas through snapping, bypassing and water locking in gas reservoirs, which will greatly reduce the gas production of a single well [13][14][15]. Thus, understanding the questions such as water invasion in gas reservoirs and its effects on gas production is significant for predicting gas productivity and improving recovery rate during the development of gas reservoirs.
The problem of water invasion in gas reservoirs has always been one of the major concerns in terms of productivity, increased operating costs, and environmental effects, and there have been many studies carried out on water invasion in gas reservoirs and its effects on gas productivity in the past years. Persoff and Pruess [16] and Zhou et al. [17] used a transparent microphysical model of a water-driven gas experiment to analyze the water invasion, and they considered that water would trap the gas in the matrix block by means of circumfluence, cut off phenomenon, and water locking. Tokunaga and Wan [18,19] conducted water films flow on fracture surfaces under near-zeromatric potentials and examines the possibility of fast unsaturated flow along the macroscopic rock fracture surface. Saad et al. [20] and Bahrami et al. [21] analyzed the problem of water coning in naturally fractured reservoirs with the experiment and field data, respectively. Perez et al. [22] applied a coning radial model to analyze the occurrence of coning in naturally fractured reservoirs. Hu et al. [23] conducted the water invasion mechanism of gas reservoirs under horizontal fracture conditions and analyzed the changes of water saturation of reservoirs at different locations from the bottom of the well during gas reservoir depletion. Azim [24] developed a fully coupled poroelastic multiphase fluid flow model to evaluate the water invasion in naturally fractured reservoir under the effects of fluid properties. Shen et al. [25] and Fang et al. [26] used full-diameter core experiments to analyze the effects of different influencing factors on water invasion in pore gas reservoirs. Based on the material balance equation, Kabir et al. [27] and Patacchini [28] analyzed and predicted the water invasion intensity and speed combined with the production data. Although there were a lot of studies conducted on the water invasion of gas reservoirs, the above studies focused on naturally fractured gas reservoirs. And water invasion in the multilayer unconsolidated sandstone gas reservoir was somewhat lacking. The unconsolidated sandstone gas reservoir in the Qaidam Basin is characterized by multiple layers and strong heterogeneity, and there exist many uncertainties in the gas production. Hence, there is an extreme necessity to understand the mechanism of generating edge water nonuniform water invasion and evaluate the effects on different edge water so as to predict gas production and optimize the producing method in the unconsolidated sandstone gas reservoir.
In this study, the physical experiment of edge water invasion in the commingling production was established to understand the water invasion process along with the reservoirs with different permeability layers, based on the geological characteristics of unconsolidated sandstone gas reservoirs. Then, the edge water invasion characteristics and mechanism of the multilayer production were analyzed. Furthermore, the effects of edge water invasion on the development of gas reservoirs such as gas phase seepage capability, recovery rate, and residual gas storage were determined. The results can provide the theoretical guidance for establishing a reasonable gas recovery rate and water control development in unconsolidated sandstone gas reservoirs.

Gas Reservoir Geological Characteristics and
Experimental Method 2.1. Gas Reservoir Geological Characteristics. The Quaternary unconsolidated sandstone gas reservoir of the Sebei gas field is located in the Qaidam Basin, Northwest China, which is mainly a structural gas reservoir [6]. The internal structures of the gas reservoir are intact without any faults developed and their lithological changes are small. The gas reservoir has many gas layers in structural high positions with large thickness while there are few gas layers in structural low positions with small thickness, which is surrounded by edge water. There exists a set of lacustrine facies beach-bar sandstone with sedimentary subfacies including shore swamp, shore lake, shallow lake and semideep lake subfacies, and sedimentary microfacies including sand bank, beach sand, sand sheet, mud bank, and marsh mud microfacies in the gas reservoir.

Experimental
Method. Based on the geological multilayer characteristics of the unconsolidated sandstone gas reservoir in the Qaidam Basin, Northwest China, the characteristics of reservoir heterogeneity and edge water invasion are considered, and a set of physical simulation experimental methods on edge water invasion during the commingling production is established to understand the edge water invasion in the unconsolidated sandstone gas reservoirs. The experimental apparatus for physical simulation of water invasion in multilayer commingled production of gas reservoirs are shown in Figure 1, which is mainly composed of the confining pressure, water body, gas source, experimental model, gas rate control, and gas-water separation measurement system. The experiment features four basic characteristics as follows: (1) four sets of cores with different magnitudes of permeability are selected to assembly connect so as to reproduce the longitudinal multigas layer geological model, and the maximum permeability difference is 12.6; (2) high-precision back pressure and confining pressure control systems are used to 2 Geofluids achieve the gas well exploitation, which can be simulated under any production rates or any production pressure difference; (3) the high-pressure transparent pipelines are connected in series between the two core holders of each group to visualize observation of the water invasion progress process as illustrated in Figure 2; (4) the encroaching water advance speed can be calculated in the experiment. When the edge water can be pushed to the transparent pipeline, it can be calculated according to the length of the core and the time for making it into the transparent pipeline. When it cannot be pushed into the transparent pipeline, the core is removed at the end of the experiment, and the advanced speed can be calculated based on the length of the core with water invasion and the time.
To understand the mechanism of edge water invasion and its effects during the multilayer commingled production in unconsolidated sandstone gas reservoir, the following experimental schemes are designed in this study: (1) four sets of cores with different permeability are used to simulate the commingled production with the vertical four gas layers, as shown in Table 1; (2) the encroaching water advance process in different permeability layers can be observed by the series connection experiment flow of high-pressure transparent pipeline; (3) different production rates (20 mL/min, 50 mL/min, 80 mL/min,100 mL/min, 150 mL/min) are adopted in the experiment to simulate the gas well exploitation process and study the effects of production rate on the water invasion; (4) the core is saturated with gas until the pore pressure is at around 4 MPa, and the outside of the gas layer is connected to a water body with a constant pressure as shown in Figure 2; (5) in the experimental process, the parameters such as experimental time, reservoir pressure, water invasion path, and gas production rate are recorded to analyze the encroaching water advance speed and the effects of edge water invasion on the gas phase seepage capacity and recovery rate as well as the residual gas storage. In this study, the characteristics of the encroaching water advance path in the gas reservoir multilayer production are recorded through the visualized experimental process, which is illustrated in Figure 1. According to the analysis results, the characteristics diagram of the encroaching water advance path is shown in Figure 2. From the result of Figure 2, it can be observed that the encroaching water advance path is mainly influenced by reservoir permeability and distribution and gas production rate. The commingling production in gas reservoirs with four layers of different permeability is conducted in the experiment; due to the permeability difference in the gas layers, the edge water exhibits the nonuniform advancement characteristics. The water invasion preferentially advances into the gas well along with two layers with the permeability of 24.4 mD and 9.34 mD, while the water invasion in the layers with the permeability of 5.37 mD and 1.93 mD is weak. The results imply that the reservoir  3 Geofluids permeability and distribution are the dominant factors of the water invasion path. When the reservoir permeability exceeds 10 mD, it will be the key channel for water invasion. Besides, the gas production rate is also an important factor on the water invasion path. When the gas production rate is 20 mL/min, the edge water is uniformly pushed along three layers with the permeability of 24.4 mD, 9.34 mD, and 5.37 mD, respectively. Once the gas production rate reaches 150 mL/min, the edge water rushes along the single direction in the reservoir with the permeability of 24.4 mD, which implies that the gas production rate is a key factor affecting water invasion advance path. In addition, when the reservoir is characterized by strong heterogeneous, a larger gas production rate will result in nonuniform water invasion easily.

Characteristics of the Encroaching Water Advance
Speed. According to the encroaching water advance path and experimental time, the calculation of the encroaching water advance speed can be expressed as follows: where v is the encroaching water advance speed, cm/min; L is the encroaching water advance path, cm; and t is the experimental time, min. Based on the experimental results, the effects and comparisons of production rate on different encroaching water advance speeds are illustrated in Figure 3 and Table 2 (no water invasion is abbreviated as NWI). From the results of Figure 2 and Table 2, it can be seen that the encroaching water advance speed is also influenced by both reservoir permeability and production rate. When the production rate is low, the edge water is relatively uniformly pushed along the multiple layers, and the encroaching water advance speed is slower. For instance, when the production rate is 20 mL/min in the experiment, the edge water is uniformly pushed along three gas layers with the permeability of 24.4 mD, 9.34 mD, and 5.37 mD, and the advance speeds are 0.64 cm/min, 0.60 cm/min, and 0.34 cm/min, respectively; when the production rate is 50 mL/min, 80 mL/min, or 100 mL/min, and the edge water permeability is 24.4 mD and 9.34 mD, it is pushed nonuniformly in both gas layers and the advance speed is obviously increased between 0.92 cm/min and 1.71 cm/min. When the production rate reaches 150 mL/min, it will rush unidirectional advance along with the gas layer with the permeability of 24.4 mD, and the encroaching water advance speed is 2.19 cm/min.
As previously mentioned, the rule of edge water invasion during the commingling production in gas reservoirs is mainly affected by the reservoir permeability and distribution and gas production rate. The two influencing factors jointly control the encroaching water advance path and speed. When the commingling production development is under the conditions with large permeability contrast and high production rate, the edge water is prone to nonuniform rush. The water invasion preferentially advances along with the high permeability layers into gas wells, and the encroaching water advance speed is fast. If the gas production rate is properly reduced, the edge water is pushed relatively uniformly along multiple layers, and the advance speed will be much slower than that in a single high permeability layer.

Nonuniform Rush of Edge Water along with High
Permeability Layers. Permeability is an important property of the reservoir rock that measures the capacity of the formation to transmit fluid, and there exist a few preferred pathways that water flow through the high permeability reservoir [29][30][31]. In this study, the threshold pressure of gas-water seepage flows in an unconsolidated sandstone reservoir is illustrated in Figure 4. Threshold pressure is the minimum pressure at which the fluid begins to flow in the low permeability reservoir. From the result of Figure 4, as can be seen, there are obvious differences in the threshold pressure of gas-water seepage flow in different permeability reservoirs. For unconsolidated sandstone reservoir with the permeability greater than 10 mD, the threshold pressure of gas flooding is generally less than 1.0 MPa. When the permeability of the unconsolidated sandstone reservoirs is no more than 10 mD, the threshold pressure of gas-water seepage flow is generally greater than 1.0 MPa. The lower the reservoir permeability is, and the higher the threshold pressure of gas-water seepage flow is. Thus, when the commingling production is in multiple gas layers with different magnitudes of reservoir permeability, gas and water will be prone to choosing preferential a seepage path along with high permeability layers under the same displacement pressure. As the displacement pressure increases, i.e., the production pressure difference is increased, gas and water can be able to form seepage flows along with the relatively low permeability reservoirs. This is the primary reason why the edge water will choose to rush along high permeability layers during commingling production in gas reservoirs.

Effect of Gas Production Rate on Nonuniform
Advancement of Water Invasion. The gas production rate is 4 Geofluids a significant parameter during the gas reservoir exploitation, which affects water invasion and ultimate recovery of gas reservoirs [8,32]. In order to understand the effect of gas production rate on nonuniform water invasion, the pressure difference between water bodies and different permeability reservoirs is measured by conducting water invasion physical simulation experiments on different production rates. Figure 5 shows the pressure difference between water bodies and different permeability reservoirs on different production rates. From the result of Figure 5, as can be seen, water invasion is influenced by gas production rates, when the pressure difference between water bodies and different permeability reservoirs is large. The gas supply is mainly from the high permeability layers when the gas production rate is more than 100 mL/min, and the pressure difference between the water body and high permeability layers is much larger than that between the water body and low-permeability layers. Due to the pressure difference between the water body and the reservoir layers, when the pressure difference is 3.99 MPa (PW − PA = 4:11 MPa, PW − PD = 0:12 MPa), the water body is pushed forward along the high permeability layers. When the gas production rate is 20 mL/min, the pressure difference between the water body and the reservoir with each scale of permeability will be reduced. For example, when the pressure difference is 2.80 MPa (PW − PA = 3:0 MPa, PW − PD = 0:2 MPa) in the experiment, which implies that if the production rate is properly controlled, the pressure differences between the edge water and the reservoir with different scales of permeability reservoir are close, and the water invasion is more easily pushed uniformly along the permeability reservoirs. When the gas production rate is high, the gas supply capacity of the high permeability layer is much larger than that of the low permeability layer, and the high permeability layer will be produced faster than the low permeability layer, which results in the pressure difference between the water body and the high permeability layers greater than that between the water body and the low permeability layers. Consequently, the water body will

Effects of the Edge Water Invasion
4.1. Effect on the Gas Phase Seepage Capability. The water production from gas-producing well is a common occurrence observed in gas reservoirs, which results in the watergas two-phase filtration and reduces gas phase seepage capability [33][34][35]. In this study, the gas permeability measurement under residual water saturation is conducted on the unconsolidated sandstone cores with different permeabilities. The relationship between conventional permeability and gas permeability under the residual water saturation is shown in Figure 6. From the result of Figure 6, it can be observed that there are significant differences in the effect of residual water saturation on the gas phase seepage capacity of different permeability sandstone reservoirs. Compared with the permeability of dry cores, there is a great effect on the gas phase seepage capacity of the reservoirs with the permeability less than 5 mD, and the gas phase permeability under residual water saturation decreases by nearly 100%. For the cores with the permeability of 5~50 mD, the gas phase permeability under residual water saturation drops of more than 90%, and the gas phase permeability under residual water saturation declines by more than 80% when the permeability exceeds 50 mD. Hence, the water will have a great influence on the gas flow seepage capacity of the reservoirs with the permeability less than 5 mD, and water invasion rushing along the high permeability layers should be avoided so as not to form water blocking during the development in the low permeability gas reservoir.

Effect on the Recovery Rate.
Due to the existence of edgeand-bottom water in gas reservoirs, the edge and bottom water will invade the gas pay zones with the gas reservoir exploitation, which leads to the lower recovery rate of gas reservoirs, and some previous studies have indicated the recovery rate of water-driven gas reservoirs was significantly lower than that of dry gas reservoirs [36][37][38]. In the study, three physical models, including no water invasion, water invasion without bypass flow, and water invasion with bypass flow, are designed to analyze the effect on the recovery rate during the production. Figure 7 illustrates the recovery rate comparison of gas reservoirs with and without water invasion. From the result of Figure 7, it can be seen that the recovery rates of the two water invasion models are significantly lower compared with the recovery rate of the no water invasion model, which indicates that it will have a major influence on gas reservoir recovery rate when the water invasion occurs in the gas reservoir development. The water invasion with bypass flow has the lowest recovery rate, which is about 30% lower than that of the no water invasion. While the recovery rate of water invasion without bypass flow ranges between 80% and 90%, which is about 10% lower than that of the no water invasion model. Although this is an experimental result, it can still reflect the changing trend, which is used as a reference during the gas reservoir development evaluation. However, the results seem to be optimistic in terms of the recovery rate itself, and it is difficult to reflect the overall gas recovery, and the simulation results can be regarded as the recovery rate of the near well areas.

Effect on the Residual Gas
Storage. The gas recovery under water drive in gas reservoirs appears to depend on an important way on the residual gas saturation, and some studies have indicated that gas recovery may be as low as 45 percent of the initial gas in place due to the water invasion [39,40]. In order to understand the effect on the residual gas storage, the residual gas ratio of each group of cores is analyzed at the end of the physical simulation experiment in this study, and the residual gas ratio is defined as the ratio between residual gas and the reserves. The experimental results of the residual gas ratio in different permeability cores are shown in Table 3. From the result of Table 3, as can be seen, there is a certain amount of residual gas that can hardly be recovered during the multilayers commingled production in gas reservoir. The ratio of residual gas to the reserves is in the range between 31% and 36%, and the residual gas is mainly distributed in the low-permeability layers. For  6 Geofluids instance, the ratio of residual gas in the gas layer with the permeability of 3.6 mD is apparently higher than that in other gas layers. The gas production rate has a certain effect on the residual gas, especially in the low-permeability layers. The larger the gas production is, and the more the residual gas is. The overall ratio of residual gas in the gas layers is 31% when the gas production rate is 20 mL/min, while the overall residual gas ratio is 36% when the production rate is 150 mL/min.

Summary and Conclusions
In this study, based on the geological characteristics of an unconsolidated sandstone gas reservoir, the physical experiment of edge water invasion was conducted to understand the mechanism of water invasion and its effect on the development phase during the commingling production. Then, the edge water invasion characteristics and mechanism of the multilayer production were discussed, and the effects such as gas phase seepage capability, recovery rate, and residual gas storage on the edge water invasion were analyzed. According to the above results, the conclusions from this study are summarized as follows: (1) During the commingled production phase, the edge water invasion is mainly affected by reservoir permeability and gas production rate, and the two influencing factors jointly control the encroaching water advance path and speed. When the gas reservoir is characterized by strong heterogeneity, a larger gas production rate will cause easily nonuniform water invasion. (2) When the commingling production is in large permeability gradation and high production rate, the nonuniform rush of edge water along the high permeability layers will occur. While there exists interflow in the high-permeability channel, the bypass flow will emerge, and water block will be formed in the lowpermeability layers and the peripheral reservoirs, which results in a significant decline in productive capacity and recovery rate. (3) The residual water saturation is a great influence on the gas phase seepage capacity of the reservoir less than 5 mD. Compared with no water invasion, the recovery rates of water invasion with and without bypass flow decrease by 30% and 10%, respectively, and there exists a certain amount of residual gas that can hardly be recovered during the commingling production, which is mainly distributed in the low-permeability layer.

Data Availability
The data used to support the findings of this study are included within the article.

Conflicts of Interest
The authors declare that they have no conflicts of interest.