Modelling and Simulation of Coupled Thermal-Hydraulic-Mechanical Processes in Geothermal ExtractionView this Special Issue
Two-Phase Flow Model for Numerical Investigation of Impact of Water Retention on Shale Gas Production
In this work, a triple-porosity, two-phase flow model was established to fill the knowledge gap of previous models focusing on gas production characteristics while ignoring the impacts of water injection. The proposed model considers the water flow in the fracture systems and clay minerals and the gas flow in the organic matter, inorganic pore, and fracture systems. The proposed model is solved using a finite element approach with COMSOL Multiphysics (Version 5.6) and verified with field data. Then, the evolutions of the intrinsic and relative permeabilities during water injection and gas production are examined. Furthermore, the impacts of water injection time and pressure are investigated. Good verification results are obtained; the goodness-of-fit value is 0.92, indicating that the proposed model can replicate both the water stimulation and the gas production stages. The relative gas permeability declines during water injection but recovers in the gas depletion stage. Furthermore, the intrinsic permeability increases during the water injection stage but decreases during the gas production stage. A higher water injection pressure and longer injection time would enlarge the intrinsic permeability, thus improving flow capacity. However, it would reduce gas relative permeability, thereby hindering gas flow. The shale gas production characteristic is controlled by the two abovementioned competing mechanisms. There exists a perfect combination of water injection pressure and injection time for achieving the maximum profitability of a shale gas well. This work can give a better understanding of the two-phase flow process in shale reservoirs and shed light on the field application of hydraulic fracturing.
Compared with other fossil fuels, shale gas plays an important role in the global energy mix with the exhaustion of conventional reservoirs [1, 2]. Owing to the extremely low permeability (nanodarcies, 10−21 m2) and porosity (less than 5%) of shale gas , highly effective reservoir stimulation technologies are vital for the successful exploitation of shale gas reservoirs . Hydraulic fracturing (HF), wherein large volumes of water are injected into shale reservoirs to induce gas production, is a mainstream method of stimulation (Figure 1(a)). Injected water flows into the hydraulic fractures and enters the shale matrix, thus occupying the gas flow channels and impairing the gas flowback process  (Figure 1(b)). Besides the injected water, abundant reserves of connate water are present in shale gas reservoirs [6, 7]; this is adsorbed on the surface of inorganic matter (clays) [8, 9] or trapped by capillary pressure in small inorganic pores [5, 10]. In general, connate water is initially immobile. However, during the water injection process, part of the penetrating water becomes residual water in the inorganic matrix, since it is filled to normal irreducible water saturation [11, 12]. In conclusion, both injected water and connate water can significantly impact gas production characteristics , yet most gas production models only consider the influence of gas pressure and effective stress or connate water ; there is no comprehensive investigation of the impacts of injected water.
Injected hydraulic water dramatically influences gas diffusivity and permeability [15, 16], thereby hindering gas transport in shale. According to experimental results for the Svalbard Shale, moisture content strongly decreases Ar and CO2 permeability . Increasing water saturation leads to an exponential decline in apparent effective He permeability; at a water saturation of approximately 30%, He permeability decreased strongly . Further, the permeability of crushed shale samples exhibited exponential decline when saturated with water . In addition, retained water can occupy void spaces and adsorption sites that were initially originally available for gas storage through two mechanisms : (i) a part of void space would be occupied by injected water, directly decreasing the free gas content in shale [21, 22]. (ii) Moreover, pore surface sites originally intended for methane adsorption are partially occupied by water molecules through competitive adsorption, since condensed water is blocking the pore throat, resulting in a decline in gas adsorption capacity [21, 23–25].
The two-phase flow process during shale gas production has been widely investigated analytically and semianalytically with the influence of complexity and half-length of hydraulic fracture [26, 27], increased permeability, and fracture conductivity in stimulated reservoir volume (SRV) [28, 29] considered. In Clarkson and Williams-Kovacs’ work, the short-term flowback data was adopted to characterize water and gas production characteristics and flowing pressure . In the later work, Williams-Kovacs and Clarkson  proposed a data-driven pseudoanalytical model in which the two-phase flowback process was investigated and the impacts of key frac properties were analyzed. However, the complex coupling mechanisms between reservoir deformation and gas-water two-phase flow are typically ignored in these previous approaches . By countering this problem, numerical reservoir models are widely used to predict production [2, 14], improve profitability , and understand potential hazards . Cao et al.  proposed a multiscale, multiphase model for gas-water two-phase flow within shale reservoirs. Based on an evaluation of shale gas recovery, the model can quantitatively differentiate these properties of fractured water. Later, a multidomain unstimulated reservoir volume (USRV) and stimulated reservoir volume (SRV), multiphysics model was proposed by Cui et al.  to represent two-phase flow within a heterogeneous shale gas reservoir. In Li et al.’s work [34, 35], a three-component solid system including organic kerogen, inorganic matrix, and hydraulic fracture (HF) was specified to characterize two-phase flow during the gas production stage. The dual-porosity, dual-permeability media was adapted to the first two components with a 1D cracked medium referred to as the last for simplification. Different properties in the so-called stimulated reservoir domain (SRD) and unstimulated reservoir domain (SRD) were assumed to simulate the influence of HF in the near-HF matrix.
As mentioned above, current gas production models are mostly focused on gas flow characteristics and ignore the impacts of the injection of hydraulic water. In other words, only the gas production stage is replicated in the simulation work, and the water injection process is seldom studied. In such approaches, the impacts of water injection pressure and injection time are not fully understood. To address these shortcomings, this study proposes a coupled two-phase flow model to comprehensively examine the effect of injected water on reservoir production. A group of partial differential equations governing the multiscale, two-phase flow in a shale gas reservoir is proposed, and both intrinsic and relative permeability evolutions are considered. The proposed model is elucidated via numerical simulations and sensitivity analyses and verified with field data.
This paper is organized as follows: the conceptual model and mathematical model are presented in Sections 2 and 3, respectively. Then, model verification is conducted in Section 4, and the results and discussion are shown in Section 5. Finally, the conclusions are explained in Section 6.
2. Conceptual Model
2.1. Triple-Porosity Model
A field emission scanning electron microscopy (FESEM) image (Figure 1(c)) [36–40] is selected to represent highly heterogeneous properties of shale matrix, such as mineral constituents, pore structure, and gas storage term. Based on these observations, a triple-porosity model was established in this study to probe the effect of hydraulic water and connate water on shale gas production characteristics. In this triple-porosity model, the shale reservoir is divided into the fracture system, inorganic matrix, and organic matter . Free gas exists in the fracture system and inorganic matrix, whereas absorbed gas is present in the organic matter . The hydraulic fluid enters the fracture system and inorganic matrix system, while connate water exists in the inorganic matrix and organic matter, impacting the gas storage amount.
2.2. Water Injection Process
(a)Initial equilibrium state and pressure rise stage
As shown in Figure 2(a), before HF, the water saturation is relatively low. As the injection process continues, the water injection volume increases linearly with time, and the water pressure in the wellbore rises continuously with continuous water injection. Water percolation in shale reservoir fractures can be described by Darcy’s law . When the water pressure reaches the fracture pressure, the shale reservoir is in a critical state, which manifests externally as hydraulic fractures. (b)HF stage
Under the combined action of in situ stress and high-pressure water, the radial stress on the inner surface of the fracture hole is expressed as compressive stress and the tangential stress is expressed as tensile stress . Because the compressive strength is much greater than the tensile strength, the newly generated fractures in the shale reservoir are mostly initiated by tensile failure . The fracture initiation and propagation mechanisms of HF are generally described by fracture mechanics failure criteria: the tensile strength criterion and the Mohr–Coulomb criterion . As shown in Figure 2(b), after HF is initiated, the injection water flows into the main fracture network, and generally the Darcy equation is applied. (c)Water imbibition stage
The water environment is shaped under the water pressure after the fracture penetrates, and the water imbibition process can be divided into three stages. First, water enters the fracture system (Darcy’s law). Then, it diffuses into the inorganic matrix. In this work, only clay minerals are assumed to be highly hygroscopic. In the last stage, the water diffuses into the organic matter through a competitive adsorption behavior. When the injection water diffuses into the shale matrix from the hydraulic or natural fractures, as shown in Figure 2(c), the gas relative permeability of the shale reservoir changes impacting the gas flow characteristics. (d)Final equilibrium state
During the fracturing process, a large amount of slick water is pumped into the reservoir for several days. A technique called soak-back was recently developed to boost production after HF; it sets a well shut-in period (soaking time) to soak the fracture fluids and thus increase productivity. During these two processes, the hydraulic fluid continues to diffuse until a new equilibrium state is reached. However, a new state may never be achieved, as thousands of years are needed .
2.3. Gas Production Process
The water flowback and gas production processes can be roughly split into four stages as shown in Figure 3. In the first stage (Figure 3(a)), water in the horizontal well and hydraulic fracture fluid flow back together with a small amount of free gas. This part is known as bubble flow. (b)Water-dominant two-phase flow
In the later stage (Figure 3(b)), the water in the inorganic matrix is depleted, since the water saturation between the hydraulic fracture and inorganic matrix is different. In this stage, some free gas in the inorganic matrix flows back because of the pressure difference. In this stage, the water flow plays the dominant role. The water saturation is higher than the gas saturation, leading to a significant larger water relative permeability. (c)Gas-dominant two-phase flow
In the middle and later stages of production (Figure 3(c)), the residual water on the fracture surfaces is extracted through moisture transport. In the meantime, the fracture permeability is improved, since the gas flow channels in the fractures are widened by the evaporation of the water film, which also increases the gas output. Additionally, the adsorbed gas in the organic matter begins to diffuse because of the gas concentration gradient. (d)Final equilibrium state
In the final state (Figure 3(d)), all the water and gas molecules are depleted, and a new equilibrium state is achieved. However, this state is difficult to achieve, as more than 3,000 years are needed to attain full equilibrium, as shown in a previous work .
2.4. Competitive Adsorption between Gas and Water
Shale gas adsorption under dry conditions plays a vital role in gas in place (GIP) estimation [47, 48]. Although a few studies have been conducted on methane adsorption in moisture-equilibrated shale [23, 49], the impact of water content has not been fully studied.
Figure 4 illustrates the microscopic differences between saturated shale (low/high moisture content) and dry shale. The adsorption sites in the shale matrix—the organic structure surfaces and the clay mineral surfaces—are fully occupied under dry conditions. However, these characteristics are quite different under moist conditions: (i) functional groups containing oxygen and nitrogen atoms are of stronger water absorption capacity and (ii) water molecules are easily adsorbed on clay minerals. Both conditions greatly reduce the shale gas adsorption capacity. Hence, a shale gas reservoir can achieve maximum adsorption capacity under dry conditions.
3. Mathematical Model
3.1. Gas-Water Two-Phase Flow in the Fracture System
The mass conservation equation of gas-water two-phase flow in a fracture system is given as  where subscript refers to either water () or gas (), subscript denotes the fracture system, and refers to the flow source or sink. With the consideration of gravity, the mass flux of the fluid can be calculated as  where and are the intrinsic and relative permeabilities, respectively, of the fracture system, represents the pore pressure, represents the viscosity, and represents the acceleration due to gravity.
The water mass () only contains free phase and is calculated as  where refers to the water saturation, represents the water density, refers to the porosity, and refers to the water source term.
The gas mass () is composed of free gas and mass resources provided by the inorganic pore system . where and represent the gas saturation and gas density, respectively.
The gas mass transfer term is proportional to the gas content difference between the fracture and inorganic matrix systems: where is the shape factor.
Therefore, the two-phase flow in a fracture system can be written as
3.2. Water Imbibition in the Inorganic Matrix
From microscopic images in some studies [54, 55], clay minerals in the shale matrix are usually embedded into nonclay minerals. Generally, a maximum water mass value exists, above which water mass remains constant, and the ratio of increase generally declines with the water mass content. Therefore, the rise in the water mass of clay minerals is determined by the disparity between the presented water mass () and the maximum value () : where refers to the adsorption rate (1/s).
3.3. Gas Slippage Flow in the Inorganic Matrix
The mass conservation law of gas slippage flow in an inorganic pore can be represented as  where subscript in represents inorganic minerals and refers to the gas mass flux.
Both the free-phase gas in inorganic pores and the gas mass source applied by the organic system contribute to the gas mass content : where refers to the inorganic porosity and refers to the gas density. is described as the modified Darcy equation: where refers to the apparent permeability and is written as 
3.4. Diffusion of Gas in Organic Matter
The mass transfer rate between the inorganic and organic systems depends on the diffusion duration in the organic matter . where represents the gas concentration at equipoise with the interface pressure . The diffusion time is represented as where subscript or represents the organic matter, (m2/s) is the surface diffusion coefficient, and is the shape factor.
3.5. Competitive Adsorption in Organic Matter
We assumed a pseudosteady state in this work  and found that the interface pressure was equivalent to the inorganic matrix pressure . Because of the presence of clay minerals, a water film exists on the surface of minerals in inorganic shale [21, 57] influencing the adsorption and desorption behavior of shale gas. The following equation is proposed : where is the gas adsorption volume in organic matter, refers to the solid-gas (Langmuir) saturated adsorption volume, is the Langmuir pressure constant, refers to the liquid-gas saturated adsorption capacity, and is the water molecular coverage ratio, defined as the rate of the wetting area of water molecules to the whole pore area. Moreover, the covering rate of the water molecules is relevant to irreducible water saturation, as shown by the following : where refers to the capillary water saturation (dimensionless), refers to the pore radius, and refers to the water molecular layer thickness.
3.6. Permeability Variation with Effective Stress and Water Saturation
3.6.1. Permeability Correction with Water Saturation
The following supplementary equations for saturation and capillary pressure are necessary for solving the two-phase flow in a shale reservoir : where is the capillary pressure and calculated with the Brooks and Corey formulation [61, 62]: where represents the nonwetting phase entry pressure, represents the effective saturation for the water phase, and represents a pore size distribution parameter. In our research, we assumed that the parameter was 2. Then, the effective water saturation can be expressed as where and represent the irreducible water and gas saturations in a fracture zone, respectively.
In this study, the relative permeability was governed by  where and refer to the relative permeabilities of gas and water, respectively.
Water saturation can be determined directly for the two-phase flow in a fracture system. The water saturation in clay minerals is related to the water mass and represented as where denotes the maximum water saturation in clay minerals.
3.6.2. Permeability Correction with Effective Stress
The fracture permeability correlation for gas shale with effective stress can be written as where refers to the effective stress and refers to the fracture compressibility. The effective stress, which represents the distinction between the formation stress and liquid pressure, can be formulated as
The liquid pressure is the combination of the gas and water pressures in the fracture system, which is defined as follows: and the inorganic system’s gas pressure, which is defined as follows:
4. Model Verification
4.1. Model Building
A testing well drilled into the lower sections of the Marcellus Shale  executed seven fracturing stages (each of which contained five perforation clusters) measuring over 630 m in the lateral section. Its completion design included 450,000 gallons of slick water per stage, and sand particles were used as proppants. Only 200 days of the production data were used for history matching. The distance between the different fracture stages far outweighed that between diverse perforation clusters in a special stage [29, 57]. Therefore, a computation domain including one fracture stage (five perforation clusters) and its adjacent area was selected.
Similar to Cui et al. , we connected the SRV areas induced per perforation cluster and formed a spheroid-shaped volume. The nonimpacted area is the USRV. In the process of simulation model design, the research area was configured as a cuboid with a size of , while the SRV exhibits spheroid shape with the semiaxes defined as 35, 100, and 50 m (in the , , and directions, respectively). The reservoir geometry and SRV are depicted in Figure 5(a). A highly refined mesh is required for the SRV region, while a relatively coarser mesh is suitable for the USRV. After meshing, there were a total of 113,895 elements with an average quantity of 0.85 as depicted in Figure 5(b).
Both the water injection, which lasted two days, and the gas production processes were replicated to fully verify the proposed model. For the water injection process, a constant water saturation value of 0.8 was used as the water flow boundary. For the gas production process, the bottom-hole pressure was obtained from Meyer . For the gas flow model in the organic and inorganic matrixes, due to the lack of immediate engagement with hydraulically formed fractures, flow boundary conditions were not taken into consideration. The water saturation value applied to the well was 0.2. Information about the Marcellus Shale reservoir was gathered from the previous work [66–68] and is detailed in Table 1.
4.2. Verification Results
The matched gas production ratios are shown in Figure 6. The model perfectly matched the field data, with a goodness-of-fit value of 0.92. The gas sources were also examined. As the production time increased, the gas production rate increased and then decreased. The increase stage was mainly observed in the first two days. This is because a large amount of water was discharged during this stage, which increased the relative gas permeability.
5. Results and Discussion
5.1. Water and Gas Transport Processes
5.1.1. Water and Gas Saturation Evolution
Changes in water and gas saturation directly affect gas and water relative permeabilities. In the water injection stage, the water pressure gradient forces the high-pressure water injected at the wellhead to migrate to the reservoir. Figure 7 shows the evolution of the and of the SRV and USRV regions. The average value of each domain is applied here. The results show the following: (i) the values of in the SRV and USRV regions show significant discrepancies, while the values of exhibit little differences. (ii) During the water flowback and gas production stage, the water in the fracture system first flows back, presenting a dramatic decrease in ; by contrast, shows little variation. (iii) As a result, the gas saturation () in the fracture system spikes, but the value in the inorganic system shows only a small change.
5.1.2. Water and Gas Pressure Evolution
The water and gas pressures are further investigated to fully illustrate the water and gas transport processes in the shale reservoir. The average water and gas pressures of the fracture system in the SRV and USRV regions are shown in Figure 8. During the water injection stage, the water pressures of both systems significantly increase in not only the SRV region but also the USRV region; the values in the SRV system are larger than those in the USRV system. In the gas production process, the water and gas pressures drop sharply to the extraction pressure at the wellhead; the value in the SRV area is lower than that in the USRV area, leading to a pressure gradient. This pressure gradient causes the gas in the fracture to continuously migrate to the wellhead.
5.2. Permeability Evolution
5.2.1. Intrinsic Permeability Evolution
Permeability determines water and gas flow rates and therefore controls gas production characteristics. It is thus a significant parameter indicating the flow capacity in a reservoir. In this research, permeability is assumed as the combination of relative and intrinsic permeabilities, and their variations are examined.
Figure 9 shows the evolutions of intrinsic and relative permeabilities during the water injection and gas production stages. In the water injection stage, the effective stress increases, leading to an increment in intrinsic permeability. The value in the SRV area is larger than that in the USRV area due to the higher water pressure in the former region. During the water flowback and gas production stage, both intrinsic permeabilities in the SRV and USRV regions decline because of the decreased effective stress. On the contrary, the value in the SRV area is lower than that in the USRV area, but the difference is small because of the supporting effect of proppants.
5.2.2. Gas/Water Relative Permeability Evolution
Besides intrinsic permeability, the evolutions of relative permeabilities are illustrated. Gas and water relative permeabilities are directly related to water saturation, which is directly proportional to water relative permeability but inversely proportional to gas relative permeability. During the water injection process, water saturation continues to increase, leading to an increment in water relative permeability and a decline in gas relative permeability. In the water flowback and gas production stage, in the SRV region, the gas relative permeability rises as the water relative permeability is reduced. In the USRV area, the gas and water relative permeabilities do not change significantly. The gas relative permeability in the gas production stage is far below the original value, thereby hindering gas flow.
5.3. Impact of Water Pressure
As discussed above, gas production behavior is closely related to the evolution of intrinsic permeability and water injection pressure. This impact is further studied through three scenarios: (i) scenario A (benchmark model), where intrinsic permeability has an exponential relationship with water injection pressure; (ii) scenario B, where intrinsic permeability is linearly proportional to water injection pressure; and (iii) scenario C, where intrinsic permeability has a Langmuir function relationship with water injection pressure. The equations of the abovementioned three scenarios could be found in Appendix A with the relationship illustrated in Figure 10.
Figure 11 depicts the evolution of the intrinsic permeability and water/gas relative permeabilities of scenarios B and C. Comparison of Figures 9 and 11(a) shows that during water injection, the intrinsic permeability in scenario B is greater than that in scenario A; more water is injected into the fracture, thereby increasing the water saturation. Consequently, the water relative permeability rises, and the gas relative permeability declines. As seen in Figure 11(b), the intrinsic permeability obtained using the Langmuir function (scenario C) is greater than that acquired using the exponential function (scenario B), leading to a larger gas relative permeability and less water relative permeability. Comparing Figures 11(a) and 11(b), we find that the intrinsic permeability shows little difference between two scenarios, but its difference between the SRV and USRV regions is smaller in scenario C. In addition, the gas permeability of scenario C is significantly larger than that of scenario B.
Figure 12 illustrates the gas production rates of three scenarios. As observed in the figure, scenario C exhibits the highest production rate in the early production stage and then dramatically declines. This is because the intrinsic permeability obtained using the Langmuir function is significantly higher than that of the two other scenarios at the same downward pressure. Scenario B behaves most poorly during the whole production process. The main reason is that the linear growth of intrinsic permeability, which is greater than the exponential rise of permeability at the same water pressure, causes more water to be injected into the fracture, which reduces the gas relative permeability. Under their combined effects, gas production is reduced.
5.4. Effect of Water Injection Time
The following scenarios are designed to sufficiently examine the effect of injection time on the gas production ratio: (i) scenario A (benchmark model), where intrinsic permeability is not affected by the water injection time; (ii) scenario B, where intrinsic permeability is linearly proportional to the water injection time; and (iii) scenario C, where intrinsic permeability has a Langmuir function relationship with the water injection time. The equations of the abovementioned three scenarios could be found in Appendix B with the relationship illustrated in Figure 13.
Figure 14 illustrates the evolutions of the intrinsic permeabilities and water/gas relative permeabilities of scenarios B and C. Comparison of Figures 9 and 14 indicate that during water injection, the intrinsic permeability of scenario B is almost twice that of scenario A. Therefore, after the water injection is completed, the water relative permeabilities in scenarios B and C are higher than that in scenario A.
Figure 15 illustrates the gas production rates of the three scenarios. From the figure, we can conclude that the production of scenarios B and C is significantly higher than that of scenario A. This is because the inherent intrinsic permeabilities of scenarios B and C are almost double that of scenario A after the water injection is completed. The gas production of scenario B is slightly higher than that of scenario C. The reason is that in scenario C, more water is injected, and the relative gas permeability is reduced.
5.5. Comprehensive Effect of Water Injection Pressure and Time
The comprehensive influences of water injection time and pressure on the gas depletion characteristic are examined by designing and solving nine computation scenarios. The water injection times are specified as 24, 48, and 72 h with the injection pressures defined as 20, 25, and 30 MPa. These two influencing factors are combined, and nine scenarios are obtained. In all scenarios, the intrinsic permeability has an exponential relationship with the water injection pressure and a linear relationship with the water injection time. The gas production ratios are shown in Figure 16. The gas production rate exhibits large differences between the scenarios, but no general variation laws can be identified. The scenario with 25 MPa and 48 h and the scenario with 30 MPa and 24 h show the largest gas production rates. The lowest values are observed in the scenario with 30 MPa and 48 h and the scenario with 30 MPa and 72 h. The performance of the rest of the scenarios is ranked randomly.
The nine scenarios are randomly arranged for the following reasons. As discussed above, both intrinsic and relative permeabilities affect gas production. In these nine scenarios, a larger injection pressure and longer injection time will result in a larger intrinsic permeability, thus enhancing gas flow. On the contrary, a higher intrinsic permeability will allow more water to flow into the reservoir, thereby decreasing gas relative permeability and hindering gas flow. Therefore, the gas production behavior is dominated by two competing mechanisms: increased intrinsic permeability and reduced gas relative permeability. Therefore, the impact of water injection pressure and water injection time on gas production should be considered systematically when designing HF and gas production processes. A higher water injection pressure and longer water injection time could improve reservoir intrinsic permeability to a greater extent; however, it would simultaneously cause more water to enter the reservoir and reduce the relative gas permeability, thereby dampening gas production.
5.6.1. Assumed Permeability Relationships
In this work, the relationships between reservoir intrinsic permeability, water injection pressure, and water injection time are artificially assumed and two different forms are proposed. Generally, the maximum tensile stress criterion or the Mohr–Coulomb criterion is applied to describe the permeability evolutions during the hydraulic fracturing process, while the main job of this work is to provide a feasible approach to investigate the impact of water injection or gas production. To achieve this goal, the complicated damage process is not considered.
5.6.2. Intrinsic Permeability Model
In this work, a simple function (Eq. (22)) was adopted to describe the intrinsic permeability evolution in not only the water injection process but also the gas production process. This approach ignores the impacts of gas adsorption, interactions between matrix and fracture systems, and in situ geostress. This would certainly affect the determined gas production characteristics. Nonetheless, this impact is negligible for the following reasons: (i) compared with coalbed methane, the gas adsorption behavior is much smaller and has little impact on the gas permeability evolution. (ii) The interactions between matrix and fracture systems change the way permeability evolves but not its trend. (iii) In situ geostress changes due to water injection and gas production are not considered in this work. Only the change in effective stress created by pore pressure variations is studied. This work focuses on the effect of water injection pressure and time, instead of permeability laws, on gas production characteristics. Therefore, a convenient and widely used permeability model was applied.
6. Summary and Conclusions
A triple-porosity, two-phase flow model was established and applied in this paper to examine the impact of hydraulic water on shale gas production characteristics. Different from the previous works which only investigate the water injection process, we paid our attention to the whole shale gas exploration process—both water injection and gas depletion processes. In particular, the evolutions of both intrinsic and relative permeabilities were investigated. According to the results of the studied case and numeric simulations, the following conclusions were drawn: (1)A shale reservoir can be divided into a fracture system, inorganic matrix, and organic matrix. HF fluid enters the fracture and inorganic matrix, with connate water existing in the inorganic matrix and organic matter. History matching with field gas production data verified that the proposed model can be applied to describe the evolution of physical quantities in not only the water injection stage but also the gas production stage in a shale reservoir(2)The evolutions of both relative and intrinsic permeabilities in the hydraulic fracturing and gas depletion processes were studied. The water relative permeability increases in the reservoir stimulation process due to water imbibition and then declines in the gas production stage. On the contrary, the gas relative permeability is below the original value because of the water imbibition. The intrinsic permeability increases during water injection but decreases during gas production(3)Gas production behavior is dominated by two competing mechanisms. Increasing the water injection pressure and prolonging the water injection time would produce more cracks, thereby increasing intrinsic permeability; more water will enter the cracks and inorganic matter, thus restricting gas production. For maximum efficiency of gas extraction, the perfect combination of an appropriate injection time and injection pressure was devised in the current simulation(4)The relationships of intrinsic permeability with water injection pressure and water injection time were particularly specified. However, the details of these relationships remain debatable, and more work is required. Nonetheless, this work can shed light on the field application of HF
A. Relationship of Permeability with Water Injection Pressure
The specific expressions and parameter values of scenarios B and C were as follows:
Scenario B: linear relationship where and .
Scenario C: Langmuir function relationship where .
The intrinsic permeability variations with respect to the water injection pressure of the three scenarios are shown in Figure 10. As illustrated in the figure, the value of scenario C is the largest and that of scenario A is the smallest.
B. Relationship of Permeability with Water Injection Time
The specific expressions and parameter values were as follows:
Scenario B: linear relationship where .
Scenario C: Langmuir function relationship where .
|:||Solid-gas saturated adsorption volume|
|:||Liquid-gas saturated adsorption capacity|
|:||Pore size distribution parameter|
|:||Flow source or sink|
|:||Adsorption pressure constant|
|:||Water molecular coverage ratio|
|:||Water molecular layer thickness|
|in:||Inorganic matrix system|
|or:||Organic matter system|
|:||Nonwetting phase entry pressure|
All data included in this study are available upon request by contact with the corresponding author.
Conflicts of Interest
We declare that there are no commercial (or related) interests and any other conflict of interest regarding the submitted work.
This research was advocated by the China Postdoctoral Science Foundation (Grant No. 2019M661118) and the “111” Project (Grant No. B17009).
A. Vafaie, B. Habibnia, and S. A. Moallemi, “Experimental investigation of the pore structure characteristics of the Garau gas shale formation in the Lurestan Basin, Iran,” Journal of Natural Gas Science and Engineering, vol. 27, pp. 432–442, 2015.View at: Publisher Site | Google Scholar
G. Cui, J. Liu, M. Wei, X. Feng, and D. Elsworth, “Evolution of permeability during the process of shale gas extraction,” Journal of Natural Gas Science and Engineering, vol. 49, pp. 94–109, 2018.View at: Publisher Site | Google Scholar
G. Cui, J. Liu, M. Wei, R. Shi, and D. Elsworth, “Why shale permeability changes under variable effective stresses: new insights,” Fuel, vol. 213, pp. 55–71, 2018.View at: Publisher Site | Google Scholar
S. Wang, C. Qin, Q. Feng, F. Javadpour, and Z. Rui, “A framework for predicting the production performance of unconventional resources using deep learning,” Applied Energy, vol. 295, article 117016, 2021.View at: Publisher Site | Google Scholar
Y. Cheng, “Impact of water dynamics in fractures on the performance of hydraulically fractured wells in gas-shale reservoirs,” Journal of Canadian Petroleum Technology, vol. 51, no. 2, pp. 143–151, 2012.View at: Publisher Site | Google Scholar
J. Xu, L. Xu, and Y. Qin, “Two effective methods for calculating water saturations in shale-gas reservoirs,” Geophysics, vol. 82, no. 3, pp. D187–D197, 2017.View at: Publisher Site | Google Scholar
J. Li, Z. X. Chen, X. F. Li et al., “A quantitative research of water distribution characteristics inside shale and clay nanopores,” Scientia Sinica Technologica, vol. 48, no. 11, pp. 1219–1233, 2018.View at: Publisher Site | Google Scholar
H. Xiong, D. Devegowda, and L. Huang, “Water bridges in clay nanopores: mechanisms of formation and impact on hydrocarbon transport,” Langmuir, vol. 36, no. 3, pp. 723–733, 2020.View at: Publisher Site | Google Scholar
J. Li, X. Li, K. Wu et al., “Water sorption and distribution characteristics in clay and shale: effect of surface force,” Energy & Fuels, vol. 30, no. 11, pp. 8863–8874, 2016.View at: Publisher Site | Google Scholar
H. Yanagihara, K. Yamashita, A. Endo, and H. Daiguji, “Adsorption–desorption and transport of water in two-dimensional hexagonal mesoporous silica,” The Journal of Physical Chemistry C, vol. 117, no. 42, pp. 21795–21802, 2013.View at: Publisher Site | Google Scholar
F. P. Wang and R. M. Reed, “Pore networks and fluid flow in gas shales,” in SPE Annual Technical Conference and Exhibition, p. 8, New Orleans, Louisiana, 2009.View at: Google Scholar
D. A. Handwerger, D. M. Willberg, M. Pagels, B. Rowland, and J. F. Keller, “Reconciling retort versus Dean Stark measurements on tight shales,” in SPE Annual Technical Conference and Exhibition, p. 14, San Antonio, Texas, USA, 2012.View at: Google Scholar
L. Xinjing, H. Suyun, and C. Keming, “Suggestions from the development of fractured shale gas in North America,” Petroleum Exploration and Development, vol. 34, p. 2, 2007.View at: Google Scholar
G. Cui, Y. Tan, T. Chen et al., “Multidomain two-phase flow model to study the impacts of hydraulic fracturing on shale gas production,” Energy & Fuels, vol. 34, no. 4, pp. 4273–4288, 2020.View at: Publisher Site | Google Scholar
W. Yuan, X. Li, Z. Pan, L. D. Connell, S. Li, and J. He, “Experimental investigation of interactions between water and a lower Silurian Chinese shale,” Energy & Fuels, vol. 28, no. 8, pp. 4925–4933, 2014.View at: Publisher Site | Google Scholar
T. Zhang, X. Li, Z. Sun et al., “An analytical model for relative permeability in water-wet nanoporous media,” Chemical Engineering Science, vol. 174, pp. 1–12, 2017.View at: Publisher Site | Google Scholar
R. van Noort and V. Yarushina, “Water, CO2 and argon permeabilities of intact and fractured shale cores under stress,” Rock Mechanics and Rock Engineering, vol. 52, no. 2, pp. 299–319, 2019.View at: Publisher Site | Google Scholar
J. Gao and Z. Li, “Water saturation-driven evolution of helium permeability in Carboniferous shale from Qaidam Basin, China: an experimental study,” Marine and Petroleum Geology, vol. 96, pp. 371–390, 2018.View at: Publisher Site | Google Scholar
M. Achang, J. C. Pashin, and E. A. Atekwana, “The influence of moisture on the permeability of crushed shale samples,” Petroleum Science, vol. 16, no. 3, pp. 492–501, 2019.View at: Publisher Site | Google Scholar
J. Zhang, Q. Feng, X. Zhang, Q. Hu, J. Yang, and N. Wang, “A novel data-driven method to estimate methane adsorption isotherm on coals using the gradient boosting decision tree: a case study in the Qinshui Basin, China,” Energies, vol. 13, no. 20, p. 5369, 2020.View at: Publisher Site | Google Scholar
Z. Jin and A. Firoozabadi, “Effect of water on methane and carbon dioxide sorption in clay minerals by Monte Carlo simulations,” Fluid Phase Equilibria, vol. 382, pp. 10–20, 2014.View at: Publisher Site | Google Scholar
J. Li, S. Wang, S. Lu et al., “Microdistribution and mobility of water in gas shale: a theoretical and experimental study,” Marine and Petroleum Geology, vol. 102, pp. 496–507, 2019.View at: Publisher Site | Google Scholar
M. Gasparik, P. Bertier, Y. Gensterblum, A. Ghanizadeh, B. M. Krooss, and R. Littke, “Geological controls on the methane storage capacity in organic-rich shales,” International Journal of Coal Geology, vol. 123, pp. 34–51, 2014.View at: Publisher Site | Google Scholar
Z. Liu, D. Liu, Y. Cai, Y. Yao, Z. Pan, and Y. Zhou, “Application of nuclear magnetic resonance (NMR) in coalbed methane and shale reservoirs: A review,” International Journal of Coal Geology, vol. 218, article 103261, 2020.View at: Publisher Site | Google Scholar
L. Wang and Q. Yu, “The effect of moisture on the methane adsorption capacity of shales: a study case in the eastern Qaidam Basin in China,” Journal of Hydrology, vol. 542, pp. 487–505, 2016.View at: Publisher Site | Google Scholar
C. R. Clarkson and J. Williams-Kovacs, “Modeling two-phase flowback of multifractured horizontal wells completed in shale,” SPE Journal, vol. 18, no. 4, pp. 795–812, 2013.View at: Publisher Site | Google Scholar
J. Zhang, Q. Feng, X. Zhang et al., “Multi-fractured horizontal well for improved coalbed methane production in eastern Ordos basin, China: field observations and numerical simulations,” Journal of Petroleum Science and Engineering, vol. 194, article 107488, 2020.View at: Publisher Site | Google Scholar
J. Williams-Kovacs and C. Clarkson, “A modified approach for modelling 2-phase flowback from multi-fractured horizontal shale gas wells. Unconventional Resources Technology Conference,” in Unconventional Resources Technology Conference, San Antonio, Texas, USA, 2015.View at: Google Scholar
R. Yang, Z. Huang, G. Li et al., “An innovative approach to model two-phase flowback of shale gas wells with complex fracture networks,” in SPE Annual Technical Conference and Exhibition, Dubai, 2016.View at: Google Scholar
X. Shang, J. Wang, and Z. Zhang, “Iterative analytical solutions for nonlinear two-phase flow with gas solubility in shale gas reservoirs,” Geofluids, vol. 2019, Article ID 4943582, 15 pages, 2019.View at: Publisher Site | Google Scholar
G. Cui, F. Xia-Ting, Z. Pan et al., “Impact of shale matrix mechanical interactions on gas transport during production,” Journal of Petroleum Science and Engineering, vol. 184, article 106524, 2020.View at: Publisher Site | Google Scholar
Y. Luo and J.-P. Ampuero, “Stability of faults with heterogeneous friction properties and effective normal stress,” Tectonophysics, vol. 733, pp. 257–272, 2018.View at: Publisher Site | Google Scholar
P. Cao, J. Liu, and Y.-K. Leong, “A multiscale-multiphase simulation model for the evaluation of shale gas recovery coupled the effect of water flowback,” Fuel, vol. 199, pp. 191–205, 2017.View at: Publisher Site | Google Scholar
W. Li, J. Liu, J. Zeng, Y. K. Leong, and D. Elsworth, “A fully coupled multidomain and multiphysics model for shale gas production,” in 5th ISRM Young Scholars' Symposium on Rock Mechanics and International Symposium on Rock Engineering for Innovative Future, Okinawa, Japan, 2019.View at: Google Scholar
W. Li, J. Liu, J. Zeng et al., “A fully coupled multidomain and multiphysics model for evaluation of shale gas extraction,” Fuel, vol. 278, article 118214, 2020.View at: Publisher Site | Google Scholar
C. Chen, D. Hu, D. Westacott, and D. Loveless, “Nanometer-scale characterization of microscopic pores in shale kerogen by image analysis and pore-scale modeling,” Geochemistry Geophysics Geosystems, vol. 14, no. 10, pp. 4066–4075, 2013.View at: Publisher Site | Google Scholar
A. Mehmani, M. Prodanović, and F. Javadpour, “Multiscale, multiphysics network modeling of shale matrix gas flows,” Transport in Porous Media, vol. 99, no. 2, pp. 377–390, 2013.View at: Publisher Site | Google Scholar
R. J. Ambrose, R. C. Hartman, M. D. Campos, I. Y. Akkutlu, and C. Sondergeld, “New pore-scale considerations for shale gas in place calculations,” in SPE Unconventional Gas Conference, Pittsburgh, Pennsylvania, USA, 2010.View at: Google Scholar
I. Y. Akkutlu and E. Fathi, “Multiscale gas transport in shales with local kerogen heterogeneities,” SPE Journal, vol. 17, no. 4, pp. 1002–1011, 2012.View at: Publisher Site | Google Scholar
A. Wasaki and I. Y. Akkutlu, “Permeability of organic-rich shale,” in SPE Annual Technical Conference and Exhibition, Amsterdam, The Netherlands, 2014.View at: Google Scholar
S. Wang, Q. Feng, F. Javadpour, M. Zha, and R. Cui, “Multiscale modeling of gas transport in shale matrix: an integrated study of molecular dynamics and rigid-pore-network model,” SPE Journal, vol. 25, no. 3, pp. 1416–1442, 2020.View at: Publisher Site | Google Scholar
S. Wang, X. Yao, Q. Feng et al., “Molecular insights into carbon dioxide enhanced multi-component shale gas recovery and its sequestration in realistic kerogen,” Chemical Engineering Journal, vol. 425, article 130292, 2021.View at: Publisher Site | Google Scholar
J. G. Wang, Y. Ju, F. Gao, and J. Liu, “A simple approach for the estimation of CO2 penetration depth into a caprock layer,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 8, no. 1, pp. 75–86, 2016.View at: Publisher Site | Google Scholar
A. Millard, A. Bond, S. Nakama, C. Zhang, J.-D. Barnichon, and B. Garitte, “Accounting for anisotropic effects in the prediction of the hydro-mechanical response of a ventilated tunnel in an argillaceous rock,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 5, no. 2, pp. 97–109, 2013.View at: Publisher Site | Google Scholar
P. Tan, Y. Jin, K. Han et al., “Analysis of hydraulic fracture initiation and vertical propagation behavior in laminated shale formation,” Fuel, vol. 206, pp. 482–493, 2017.View at: Publisher Site | Google Scholar
G. Cui, J. Wei, X.-T. Feng et al., “Preliminary study on the feasibility of co-exploitation of coal and uranium,” International Journal of Rock Mechanics and Mining Sciences, vol. 123, article 104098, 2019.View at: Publisher Site | Google Scholar
T. Y. Chen, X. T. Feng, and Z. J. Pan, “Experimental study of swelling of organic rich shale in methane,” International Journal of Coal Geology, vol. 150-151, pp. 64–73, 2015.View at: Publisher Site | Google Scholar
T. Chen, X.-T. Feng, and Z. Pan, “Experimental study on kinetic swelling of organic-rich shale in CO2, CH4 and N2,” Journal of Natural Gas Science and Engineering, vol. 55, pp. 406–417, 2018.View at: Publisher Site | Google Scholar
F. Yang, C. Xie, Z. Ning, and B. M. Krooss, “High-pressure methane sorption on dry and moisture-equilibrated shales,” Energy & Fuels, vol. 31, no. 1, pp. 482–492, 2017.View at: Publisher Site | Google Scholar
T. Ma, K. Zhang, W. Shen, C. Guo, and H. Xu, “Discontinuous and continuous Galerkin methods for compressible single-phase and two-phase flow in fractured porous media,” Advances in Water Resources, vol. 156, article 104039, 2021.View at: Publisher Site | Google Scholar
T. Ma, J. Rutqvist, C. M. Oldenburg, W. Liu, and J. Chen, “Fully coupled two-phase flow and poromechanics modeling of coalbed methane recovery: impact of geomechanics on production rate,” Journal of Natural Gas Science and Engineering, vol. 45, pp. 474–486, 2017.View at: Publisher Site | Google Scholar
Z. Wei and D. Zhang, “A fully coupled multiphase multicomponent flow and geomechanics model for enhanced coalbed-methane recovery and CO2 storage,” SPE Journal, vol. 18, no. 3, pp. 448–467, 2013.View at: Publisher Site | Google Scholar
J. G. Wang and Y. Peng, “Numerical modeling for the combined effects of two-phase flow, deformation, gas diffusion and CO2 sorption on caprock sealing efficiency,” Journal of Geochemical Exploration, vol. 144, pp. 154–167, 2014.View at: Publisher Site | Google Scholar
W. Yuan, Z. Pan, X. Li et al., “Experimental study and modelling of methane adsorption and diffusion in shale,” Fuel, vol. 117, pp. 509–519, 2014.View at: Publisher Site | Google Scholar
L. Wang, M. Bornert, E. Héripré, S. Chanchole, A. Pouya, and B. Halphen, “Microscale insight into the influence of humidity on the mechanical behavior of mudstones,” Journal of Geophysical Research: Solid Earth, vol. 120, no. 5, pp. 3173–3186, 2015.View at: Publisher Site | Google Scholar
J. G. Wang, A. Kabir, J. Liu, and Z. Chen, “Effects of non-Darcy flow on the performance of coal seam gas wells,” International Journal of Coal Geology, vol. 93, pp. 62–74, 2012.View at: Publisher Site | Google Scholar
R. Aringhieri, “Nanoporosity characteristics of some natural clay minerals and soils,” Clays and Clay Minerals, vol. 52, no. 6, pp. 700–704, 2004.View at: Publisher Site | Google Scholar
W. Shen, X. Li, T. Ma, J. Cai, X. Lu, and S. Zhou, “High-pressure methane adsorption behavior on deep shales: experiments and modeling,” Physics of Fluids, vol. 33, no. 6, article 063103, 2021.View at: Publisher Site | Google Scholar
J. Li, X. Li, X. Wang et al., “Water distribution characteristic and effect on methane adsorption capacity in shale clay,” International Journal of Coal Geology, vol. 159, pp. 135–154, 2016.View at: Publisher Site | Google Scholar
J. Zhang, B. Zhang, S. Xu, Q. Feng, X. Zhang, and D. Elsworth, “Interpretation of gas/water relative permeability of coal using the hybrid Bayesian-assisted history matching: new insights,” Energies, vol. 14, no. 3, p. 626, 2021.View at: Publisher Site | Google Scholar
R. H. Brooks and A. T. Corey, “Hydraulic properties of porous media and their relation to drainage design,” Transactions of the ASAE, vol. 7, no. 1, pp. 0026–0028, 1964.View at: Publisher Site | Google Scholar
M. T. Van Genuchten, “A closed-form equation for predicting the hydraulic conductivity of unsaturated soils,” Soil Science Society of America Journal, vol. 44, no. 5, pp. 892–898, 1980.View at: Publisher Site | Google Scholar
M. Leverett, “Capillary behavior in porous solids,” Transactions of the AIME, vol. 142, no. 1, pp. 152–169, 1941.View at: Publisher Site | Google Scholar
D. S. Lee, J. D. Herman, D. Elsworth, H. T. Kim, and H. S. Lee, “A critical evaluation of unconventional gas recovery from the Marcellus Shale, northeastern United States,” KSCE Journal of Civil Engineering, vol. 15, no. 4, pp. 679–687, 2011.View at: Publisher Site | Google Scholar
B. R. Meyer, L. W. Bazan, R. H. Jacot, and M. G. Lattibeaudiere, “Optimization of multiple transverse hydraulic fractures in horizontal wellbores,” in SPE Unconventional Gas Conference, Pennsylvania, USA, 2010.View at: Google Scholar
W. Yu and K. Sepehrnoori, “Simulation of gas desorption and geomechanics effects for unconventional gas reservoirs,” Fuel, vol. 116, pp. 455–464, 2014.View at: Publisher Site | Google Scholar
W. Yu, T. T. Zhang, S. du, and K. Sepehrnoori, “Numerical study of the effect of uneven proppant distribution between multiple fractures on shale gas well performance,” Fuel, vol. 142, pp. 189–198, 2015.View at: Publisher Site | Google Scholar
Y. Tan, Z. Pan, J. Liu, X.-T. Feng, and L. D. Connell, “Laboratory study of proppant on shale fracture permeability and compressibility,” Fuel, vol. 222, pp. 83–97, 2018.View at: Publisher Site | Google Scholar