Next Article in Journal
Development on Surrogate Models for Predicting Plume Evolution Features of Groundwater Contamination with Natural Attenuation
Previous Article in Journal
Geochemical Characteristics and Genesis of Brine Chemical Composition in Cambrian Carbonate-Dominated Succession in the Northeastern Region of Chongqing, Southwestern China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of Gas–Water Two-Phase Flow Patterns in Fracture: Implication for Enhancing Natural Gas Production

1
State Key Laboratory of Intelligent Construction and Healthy Operation & Maintenance of Deep Underground Engineering, Xuzhou 221116, China
2
College of Mining Engineering and Geology, Xinjiang Institute of Engineering, Urumqi 830091, China
*
Author to whom correspondence should be addressed.
Water 2024, 16(19), 2860; https://doi.org/10.3390/w16192860
Submission received: 3 September 2024 / Revised: 6 October 2024 / Accepted: 7 October 2024 / Published: 9 October 2024
(This article belongs to the Section Hydraulics and Hydrodynamics)

Abstract

:
The main objective of this paper is to investigate the evolution of rock fracture slug structures and decongestion strategies for natural gas extraction processes. For this purpose, the level set method was used to simulate the evolution of the slug structure under the effect of different flow ratios, fracture surface wettability, and fracture tortuosity. The results show that an increase in the water-to-gas flow ratio and fracture tortuosity leads to a significant increase in the proportion of slug structures in the fracture, while an increase in the surface contact angle leads to a decrease in the proportion of slug structures in the fracture. Based on the above slug structure evolution law, a quantitative characterization method for the slug structure of two-phase fluids considering the combined effects of the water–gas flow ratio, average wall contact angle, and flow channel tortuosity was developed. Subsequently, we engage in further discussion on the optimization of the extraction and decongestion process in natural gas extraction.

1. Introduction

Reducing carbon emissions has become a global consensus and the demand for clean energy is increasing worldwide [1,2]. As an important component of clean fossil energy, coalbed methane has gradually become the focus of global exploration and development [3], and abandoned mine gas extraction is gradually becoming a new way of coalbed methane extraction. Since the permeability of fractures in abandoned mines is much higher than that of the matrix, fractured gas reservoirs occupy an important position in the coalbed methane extraction process [4]. The complexity of two-phase fluid flow arises mainly from the complex flow interface between the wetting and non-wetting phases [5]. Different flow structures imply different flow characteristics. Currently, common flow structure classifications include bubble flow, slug flow, stratified flow, etc. [6,7]. Among them, the main feature of slug flow is a series of long bubbles separated by slugs [8]. It is precisely for this reason that the pressure in the flow direction exhibits significant fluctuations. During the coalbed methane extraction process in abandoned mines, the presence of slug flow is a serious limitation to stable extraction [9]. Therefore, the investigation of the gas–liquid slug flow in rock fractures and its controlling factors is of great importance to instruct reservoir applications.
In a study of two-phase fluid flow, the flow structure map consisting of the flow structure under the action of different two-phase fluid apparent flow velocities is usually used to reflect the evolution law of the two-phase fluid flow structure [10,11]. However, the flow structure map only reflects the flow interface where the flow structure is significantly transformed and fails to reflect the variability of the same flow structure under the action of different two-phase fluid flow rates. For the slug structure, there are also significant differences in the two-phase fluid flow characteristics represented by different slug lengths and the number of slugs in a rock fracture. Therefore, it is necessary to study the evolution of the slug flow structure under the action of different flow ratios to improve the understanding of the two-phase fluid slug flow structure.
Fracture tortuosity variation affects two-phase fluid flow behavior [12,13]. The two-phase fluid flow structure in rough fractures is more complex and tortuous [14,15].Through visualization tests, it was found that the two-phase fluid has an intermittent flow structure in the rough fracture and the fluid channels in the fracture become more concentrated as the fracture tortuosity and inhomogeneity increase [16]. Meanwhile, similar investigation results were also obtained using a numerical simulation [17]. However, the evolution pattern of the slug structure under the influence of fracture tortuosity has not been adequately studied.
The wetting angle is the angle between the liquid–solid interface and the tangent line of the liquid surface at the contact point between the liquid phase and the solid phase [18]. The variation in it changes the flow structure of the two-phase fluid. The research shows significant differences in flow structures in fractures with different wettability [19]. For droplet and bubble wettability, related reports have appeared in recent years [20]. However, little has been reported on the evolution of the slug flow structure in fractures under the action of surface wettability. Meanwhile, due to the non-homogeneity of rock mineral components, the wetting angle of the fracture surface shows obvious non-homogeneity [21]. However, most of the studies on the wettability of fracture surfaces are based on the average wettability angle of the fracture surface, so the current studies on the two-phase fluid flow structure in rough fractures differ greatly from the two-phase fluid distribution in real rock fractures [19]. At present, there are still very few research reports on the evolution of two-phase fluid flow structures in rough fractures under the influence of surface wetting angle inhomogeneity.
To study the evolution of slug flow structures in rough rock fractures, the transformation of slug length and slug number with variation in the two-phase fluid flow ratio, fracture roughness and the surface contact angle are investigated by performing computational fluid dynamics modeling. A quantitative characterization method for the evolution of the fracture slug flow structure is proposed. Finally, a simple prediction method for the two-phase fluid flow structure in rough fractures is proposed by comprehensively analyzing the range of values of the characteristic parameters of the slug flow structure.

2. Numerical Model of Two-Phase Flow

2.1. Governing Equations of Two-Phase Fluid Flow

The diffusion effect at the interface is equal for an incompressible and insoluble two-phase Newtonian fluid flow via rock fractures, and the convection velocity complies with the requirement that the divergence be zero [22,23,24]. The Navier–Stokes equations with time dependence and the continuity equation are used as follows:
u = 0
ρ u t + ρ ( u ) u = [ p I + μ ( u + ( u ) T ] + F g + F s t + F
where u is the velocity vector, ρ is the fluid density, p is the pressure, I is the identity matrix, Fg is the gravitational force, F commprises the other external forces in the system, and Fst is the surface tension force which can be expressed as [25]
F s t = [ ( τ ( I n n T ) ) δ ]
where n is the unit normal to the interface, τ is the surface tension, and δ is the Dirac delta function.
The following are definitions for fluid density and dynamic viscosity:
ρ = ρ w + ( ρ g ρ w ) φ
μ = μ w + ( μ g μ w ) φ
where φ is a smooth continuous level set function, ρg is the gas density, ρw is the water density, μg is the dynamic viscosity of the gas, and μw is the dynamic viscosity of water.
Due to the current lack of experimental and theoretical studies on gas–water two-phase flow under high-temperature and high-pressure conditions, numerical simulations based on these conditions face challenges in being fully validated through relevant experiments or theoretical models. As a result, many researchers have used the gas–water physical and mechanical parameters under ideal conditions to ensure the feasibility of experiments and improve the reliability of simulation results. For computational fluid dynamics modeling, the ambient temperature is 20 °C, the pressure differential is simulated to be low, and the gas dissolution in water is neglected. Moreover, the atmospheric pressure is much higher than the gas pressure drop, so the compressibility of the gas is not considered. Table 1 lists the physical properties of the two-phase fluids.
The two-phase fluid flow interface in the level set method is described by the following advection equation:
φ t + v φ = 0
A constant interface thickness is used for numerical modeling and smooth parameter variation is allowed. A phase reinjection scheme is added to the two-phase fluid flow interface equation to enhance the stability of the solution process. The transfer equation of the fluid interface is expressed as
φ t + v φ = γ ( ε φ φ ( 1 φ ) φ φ )
where ε and γ are model parameters. ε is related to the mesh size and controls the thickness of the interface and γ is related to the velocity field and controls the re-initialization and stabilization amount. In this study, γ is 1 m/s.

2.2. Flow Channel

As this chapter aims to investigate the effect of the inhomogeneity of the fracture sur-face wettability on the flow behavior of the gas–liquid two-phase flow, a flow channel consisting of elements with different tortuosity is selected for this chapter [26]. Figure 1 shows a schematic diagram of the flow channel. By dividing the two-phase fluid flow channel into ten regions from the middle and setting the wettability parameters in each of them, the two-phase flow behavior in a non-homogeneous surface wettability rock fracture can be simulated. Meanwhile, to consider the effect of fracture roughness on the slug structure, the distance between the center of the U-shaped channel and the outermost edge (D) is changed to form different flow channel tortuosity to reflect the change in fracture surface roughness.

2.3. Boundary Condition

In the whole simulation process, the gas and liquid phases are injected from the left side and flow out from the right side of the horizontal fracture. Velocity control is used at the inlet side, where the liquid phase is fixed at 0.5 m/s and the gas phase is set according to the water–gas fluid flow ratio (n) of 0.5, 1, 1.5, and 2, respectively, to meet the needs of different gas–liquid flow ratio simulations. The two-phase fluid water–gas flow ratio (n) can be expressed as
n = v w v g
The outlet side is uniformly set to a zero-pressure outlet boundary (P = 0) by using pressure control. To analyze the effect caused by a different surface wettability on the two-phase fluid transport, the wall contact angles were set according to a hydrophilic fracture surface (30°), a no significant wettability fracture surface (90°), a hydrophobic fracture surface (120°), and a superhydrophobic fracture surface (150°). The contact angle formed at each stage of the simulation conforms to Young’s equation. The Young–Laplace equation is expressed as
P c = P g P w = σ κ 0
where Pc is capillary pressure, σ is the interfacial tension, and κ0 is the interface curvature. The slip boundary conditions in numerical modeling are directly related to the wettability of the fracture surface through the two-phase fluid flow. In this study, wall slip is set to the Navier slip condition and the minimum length factor of slip length is set to one.
For the non-homogeneous boundary condition setting of the fracture surface, we assign the surface contact angle parameters to the two fracture wall sections by generating the same mean value but with deviation values (d) of 0°, 5°, 10°, 15°, and 20° to make each section exhibit significant non-homogeneity.
For different fracture tortuosity boundary conditions, the distance (D) between the center of the U-shaped channel and the outermost edge was set to 0.1, 0.2, 0.3, and 0.4 mm to make the flow channels show different tortuosity. In this study, the non-homogeneity of the fracture pore size distribution is disregarded, and all fracture pore sizes are set to 0.1 mm.

2.4. Grid Division and Solver Settings

The generated rough fractures are modeled using Comsol Multiphysics 6.0. In the numerical simulation, the flow parameters of each phase of the two-phase fluid were calculated and the control equations were solved using a separation solver. A generalized minimal residual solver (GMRES) was used to calculate the model parameters with an error of 10−3. To optimize convergence, the GMRES was solved in combination with a preprocessing method. The parallel sparse direct linear solver (PARDISO) was used as a hyperfine mesh preprocessor to investigate the phase initialization and time dependence problems. Since the calculations are time-dependent, the simulation time range is 0–0.02 s with a time step of 0.001 s. To avoid the effect of different rough fractures due to different mesh sizes, a mesh size of 0.3 μm~5 μm is used for the whole simulation. The maximum unit growth rate was 1.13 and the curvature factor was 0.3.

2.5. Model Verification

The validity and scientific validity of using the level set method to study the structural evolution of two-phase fluid flow has been reported in our previous studies. Meanwhile, the structural evolution of two-phase fluid flow in a trickle bed reactor by the level set method was also reported by [19]. To verify the effectiveness of the level set method in characterizing the slug structure of two-phase fluids, we used the simulation results to compare with different slug structures reported by [26]. The obtained slug structures are in good agreement with the slug structures found in the experiments (Figure 2). In this paper, the red color represents water, and the blue color represents gas.

3. Evolutionary Law of Slug Structure

3.1. Effect of Water–Gas Flow Ratio

Figure 3 shows the evolution of the slug structure with the increase in the water–gas flow ratio when the mean value of the contact angle (θ) of the fracture surface is 90° and the deviation value (d) is 0°. It can be seen that the two-phase fluid flow structure in the fracture is slug flow under several two-phase fluid flow ratio conditions, but the slug section structure in the fracture also shows variability as the two-phase flow ratio increases. The average length of the gas slug section in the fracture is shortened by nearly two times when the water–gas flow ratio rises from 0.5 to 2. The number of slug sections in the fracture gradually increases and the length of the single liquid slug section gradually decreases as the two-phase fluid flow ratio rises. In addition, when n = 0.5, the longest gas slug section present in the fracture is about double the length of the shortest gas slug section, while when n = 2, the lengths between individual gas slug sections in the fracture are basically the same. This indicates that the structural length of gas slug segments in the fracture is more stable when the water–gas flow rate is relatively high. Detailed data for the simulations in this paper are presented in Appendix A.
For the slug flow structure, the length of a slug section is certainly a quantitative description index and likewise, the number of slug sections is also an important parameter to evaluate the slug flow structure in the fracture. It can be seen that the number of slug segments increases significantly when the water–gas flow ratio increases from 0.5 to 1.5, and the number of slug segments in the fracture begins to stabilize from region to region when the water–gas flow ratio increases from 1.5 to 2. This is because the change in the slug flow structure in the fracture when the water–gas flow ratio increases from 0.5 to 1.5 is mainly reflected in the decrease in the length of the gas slug segments and the increase in the number of segments, while the liquid slug segments are much less affected than the gas slug segments at this stage. At this stage, the slug flow structure shows a significant increase in the number of slug sections. But when the water–gas flow ratio is higher than 1.5 and continues to rise, the fracture is gradually dominated by water when the liquid slug structure is stable, and the main change trend in the slug flow structure in the fracture changes in the direction of increasing the length of the gas slug section and decreasing the number of segments.

3.2. Effect of Mean Contact Angle

Figure 4 shows the two-phase fluid flow structure distribution in the fracture under the effect of the mean surface contact angle when the water–gas flow (n) ratio is one and the deviation value (d) is 0°. It can be seen that the two-phase fluid flow structure in the fracture is also slug flow under several average surface contact angle conditions. However, unlike the slug flow structure under the action of the water-to-gas flow ratio, the structural instability of the slug section in the fracture increases significantly with the increase in the average wall contact angle, and the slug section in the fracture starts to become different lengths.
Taking the average contact angle of 90° as the benchmark, it can be found that when the contact angle decreases, the fracture presents as hydrophilic, and the transformation of the slug flow structure in the fracture is mainly reflected in the elevated length of the gas slug section and the transient stratified flow structure of the two-phase fluid near the inlet. Meanwhile, when the contact angle of the fracture is higher than 90°, many short slug segments and even the presence of some slug segments in the fracture begin to convert to a droplet structure.
The evolution of the flow structure for two-phase fluids is analyzed on the basis of the apparent flow velocities of different two-phase fluids. Normally, the two-phase fluid flow structure in the fracture develops gradually along the bubble flow–slug flow–annular flow sequence as the apparent gas flow rate increases. However, the results in Figure 4 show that the flow structure in the fracture similarly develops gradually along the bubble flow–slug flow–annular flow sequence with the decrease in the average contact angle.

3.3. Effect of Non-Homogeneous Wettability

Table 2 lists the value of the surface contact angle at each stage of fracture. Figure 5 shows a schematic diagram of the change in the slug structure due to surface wetting heterogeneity. Figure 6 shows the two-phase fluid flow structure for different surface wettability heterogeneities. There is some variability in the slug flow structure in different regions, and the structural distribution characteristics exhibited are generally consistent with the evolution of the two-phase fluid flow structure under different mean contact angle effects found in the previous section.
The flow structure of the two-phase fluid in the fracture exhibited locally depends mainly on the two-phase fluid flow ratio and the surface contact angle at the local location, and when the fracture surface wettability cannot be reflected by a single contact angle size, the two-phase fluid flow structure in the fracture also shows obvious inhomogeneity. However, unlike the water–gas flow ratio and the average surface contact angle on the number of slug sections in the channel, there is a clear pattern of evolution as the contact angle variance in the fracture surface increases, causing the slug structure to change, but there is no clear trend in this change.

3.4. Effect of Fracture Tortuosity

Figure 7 shows the two-phase fluid flow structure distribution in the fracture under the influence of fracture tortuosity when the water–gas flow ratio is one and the mean value of the wall contact angle is 90° with zero deviation values. It can be seen that under several fracture tortuosity conditions, the two-phase fluid flow structure in the fracture is also slug flow; with an increase in fracture tortuosity, the structural instability of the slug section in the fracture increases significantly, and the slug section in the fracture starts to become different lengths. Meanwhile, when the fracture tortuosity decreases, the transformation of the slug flow structure in the fracture is mainly reflected in the decrease in the number of gas slug segments and the decrease in the gas content in the fracture.

4. Quantitative Characterization Method for Two-Phase Fluid Slug Structure

To investigate the evolution of the slug structure in the fracture, the ratio of the total length of the slug section in the fracture to the total length of the flow channel (LD) is introduced in this section to describe the slug structure. Figure 8 shows a schematic diagram of gas plunger segments within a flow channel. LD in the figure is equal to the ratio of the total length of the two slugs to the total length of the flow channel.

4.1. Effect of Water–Gas Flow Ratio on Slug Structure

To investigate the evolution of the slug structure in the fracture under the effect of the water–gas flow ratio and to better study the flow characteristics of two-phase fluid flow in the fracture, the evolution of the slug structure in the fracture under the effect of different flow ratios is studied. Figure 9 shows the relationship between the water–gas flow ratio and the total length of the slug section to the total length of the flow channel (LD) in the fracture. LD increases linearly as the water–gas flow ratio increases from 0.5 to 2. Table 3 lists the values of the correlation coefficient (R2).
These correlation coefficients are all higher than 0.9, indicating a significant linear relationship between the water–gas flow ratio and LD. The above linear relationship can be expressed as the following equation:
L D = a 1 n + b 1
where a1 and b1 are fitting coefficients.

4.2. Effect of Mean Contact Angle on Slug Structure

Figure 10 shows the relationship between the mean contact angle and the total length of the slug section to the total length of the flow channel (LD) in the fracture. LD decreases linearly as the mean contact angle increases from 30° to 150°. Table 4 lists the values of the correlation coefficient (R2). These correlation coefficients are all higher than 0.9, indicating a significant linear relationship between the mean contact angle and LD. The above linear relationship can be expressed as the following equation:
L D = a 2 θ + b 2
where a2 and b2 are fitting coefficients.

4.3. Effect of Non-Homogeneous Wettability on Slug Structure

Figure 11 shows the evolution of the ratio of the total length of the slug section to the total length of the flow channel in the fracture under the action of different heterogeneities of surface wettability. It can be seen that LD is of a high level when d is zero and decreases with an increase in d. However, it is difficult to establish a uniform functional relationship between these two parameters.
Therefore, the effect of the non-homogeneity of surface wetting on the two-phase fluid slug flow structure can be summarized qualitatively. The non-homogeneity of the contact angle of the fracture wall can lead to the difficulty of stable formation of the slug structure and LD decreases; however, LD is difficult to be characterized quantitatively by constructing a function of the non-homogeneous wettability on the fracture surface.

4.4. Effect of Fracture Tortuosity on Slug Structure

Figure 12 shows the relationship between the fracture tortuosity and the total length of the slug section to the total length of the flow channel (LD) in the fracture. LD increases linearly as the fracture tortuosity increases from 0.1 to 0.4 mm. Table 5 lists the values of the correlation coefficient (R2). These correlation coefficients are all higher than 0.9, indicating a significant linear relationship between the mean contact angle and LD. The above linear relationship can be expressed as the following equation:
L D = a 3 D + b 3
where a3 and b3 are fitting coefficients.

4.5. Quantitative Characterization Method for the Evolution of Slug Structure

Equations (10)–(12) show that the relationship between the LD and the water–gas flow ratio, the mean contact angle, and the flow channel tortuosity satisfies the linear functional form. The LD of the slug structure can be expressed by Equation (13) [26].
L D = A 1 n + A 2 θ + A 3 D + B 1
where A1, A2, A3, and B1 are the best-fitting coefficients. Based on all simulation results, the best-fit equation is
L D = 0.046 × n 9.4 × 10 4 × θ + 0.133 × D + 0.507
The correlation coefficient (R2) is 0.96, showing that Equation (14) can describe the simulation data well. It is worth mentioning that in combining the value of LD with the flow structure in the corresponding state, it can be found that LD ranges from 0.4 to 0.6 when there is a slug structure in the fracture. When LD is lower than 0.4, sparse bubbles or droplet structures appear in the fracture (Figure 4a), while a large number of gas–liquid layered flow structures exist in the fracture when it is higher than 0.6 and less than one (Figure 4d).
Based on the above studies, it is known that the water–gas flow ratio, surface wettability, and tortuosity change the slug structure in the fracture by altering the gas-to-water ratio, spatial distribution, and gas-to-water interfacial forces within the fracture.

5. Optimization of Natural Gas Extraction Decongestion Technology

The experiment results show that the formation of slug flow is closely related to the two-phase fluid flow ratio, fracture tortuosity, and fracture surface wettability. To ensure efficient gas production, the gas–water two-phase fluid in the fracture should be maintained in an annular flow structure as far as possible (0.6 < LD < 1). Figure 13 shows the relationship between the flow structure and LD. There are several issues that should be considered for the efficient and stable extraction of natural gas.
The slug structure is more likely to be found in tortuous fractures than in straight fractures. Optimizing the pumping process to reduce the curvature of fractures in complex fracture networks can therefore help to prevent the formation of slug flow. In addition, optimization of the proppant composition is also effective in reducing the surface roughness of the fracture and thus preventing the formation of a slug structure.
The elevated surface wettability drives the formation of the slug structure. Therefore, the effect of surface wettability on the slug structure is reduced by selecting a more suitable proppant to ensure that the fracture opening is increased.

6. Conclusions

This paper explores the effects of the two-phase fluid flow ratio, fracture tortuosity, and surface wettability on the slug structure of two-phase fluids in fractures during natural gas extraction using the horizontal set method. The study aims to identify key strategies to prevent the emergence of slug structures during natural gas extraction. The key findings are summarized as follows:
  • Controlling the flow ratio of the water–gas phase at the inlet can have a significant effect on the slug structure in the fracture. An increase in the water–gas flow ratio resulted in an increase in the number of slug segments and a decrease in the length of individual liquid slug segments. When the water–gas flow ratio increased from 0.5 to 2, the average length of the gas slug segment in the fracture nearly doubled. Additionally, when n = 0.5, the longest gas slug section present in the fracture was approximately double the length of the shortest gas slug section, while when n = 2, the length of each gas slug section in the fracture was essentially the same. This suggests that the structure of the gas slug segments in the fracture is more stable when the water–gas flow rate is relatively high.
  • As the average wall contact angle increases, the structural instability of the slug section in the fracture increases. With an average contact angle of 90°, the transformation of the slug flow structure in the fracture is mainly reflected in the increase in the length of the gas slug section and the transient layered flow structure of the two-phase fluid near the inlet. When the contact angle of the fracture wall is higher than 90°, many short liquid slug segments appear in the fracture and some liquid slug segments start to transform into a small droplet structure.
  • The non-homogeneity of the contact angle at the fracture wall can lead to a difficult formation of the slug structure and a decrease in the ratio of the total length of the slug section in the fracture to the total length of the flow channel. However, it is difficult to characterize quantitatively by constructing a function of the variance in the contact angle at the fracture surface.
  • A quantitative characterization method for the evolution of slug structures in cracks, considering a two-phase fluid slug structure in terms of the water–gas flow ratio, mean wall contact angle, and flow channel torque, has been constructed. The results show that the method can accurately reflect the evolution of gas–liquid two-phase fluid slug structures in fractures.
  • A comprehensive analysis of the evolution of the flow structure in the fracture and the correlation between the corresponding parameter LD reveals that the corresponding LD values range from 0.4 to 0.6 when a slug flow structure is present in the fracture.

Author Contributions

Methodology, D.L. and H.P.; investigation, K.X.; resources, D.L.; data curation, K.X. and D.L.; writing—original draft preparation, D.L.; writing—review and editing, D.L. and K.X.; visualization, H.N. and D.L.; supervision, H.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 2023YFC3804204, No. 52374147, No. 2023A01002, and No. 41931284).

Data Availability Statement

Data are available upon request due to privacy restrictions.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. Summary table of numerical simulation parameters related to Figure 9.
Table A1. Summary table of numerical simulation parameters related to Figure 9.
nCase 1Case 2Case 3Case 4
10.528930.481430.436630.41884
1.50.55060.512310.480910.44302
20.567040.53840.503890.47891
2.50.585680.564450.521290.49436
Table A2. Summary table of numerical simulation parameters related to Figure 10.
Table A2. Summary table of numerical simulation parameters related to Figure 10.
θ (°)Case 1Case 2Case 3Case 4
300.528930.55060.567040.59568
900.481430.512310.53840.56445
1200.436630.489910.503890.52129
1500.418840.443020.478910.49436
Table A3. Summary table of numerical simulation parameters related to Figure 12.
Table A3. Summary table of numerical simulation parameters related to Figure 12.
D (mm)Case 1Case 2Case 3Case 4
0.10.528930.481430.436630.41884
0.20.562640.512110.464460.44553
0.30.573510.5220.473430.45414
0.40.584840.533220.484530.46519

References

  1. Luo, P.; Zhang, Z.; Geng, X.; Xue, K.; Guang, W. Evaluation of ScCO2-water performance on bituminous coal: Insights from experiments and 3D CT image reconstruction. Geomech. Geophys. Geo-Energy Geo-Resour. 2022, 8, 118. [Google Scholar] [CrossRef]
  2. Huang, X.; Zhang, L.; Zhang, R.; Chen, X.; Zhao, Y.; Yuan, S. Numerical simulation of gas-liquid two-phase flow in the micro-fracture networks in fractured reservoirs. J. Nat. Gas Sci. Eng. 2021, 94, 104101. [Google Scholar] [CrossRef]
  3. Ren, L.; Lin, R.; Zhao, J. Stimulated Reservoir Volume Estimation and Analysis of Hydraulic Fracturing in Shale Gas Reservoir. Arab. J. Sci. Eng. 2018, 43, 6429–6444. [Google Scholar] [CrossRef]
  4. Bagherzadeh, P.; Goshtasbi, K.; Kazemzadeh, E.; Kashef, M.; Aloki Bakhtiari, H. Stress-dependence of the permeability, porosity, and compressibility in fractured porous media regarding fracturing condition. Bull. Eng. Geol. Environ. 2021, 80, 5091–5110. [Google Scholar] [CrossRef]
  5. Zhu, Q.; Wang, G.; Schlegel, J.P.; Yan, Y.; Yang, X.; Liu, Y.; Ishii, M.; Buchanan, J.R. Experimental study of two-phase flow structure in churn-turbulent to annular flows. Exp. Therm. Fluid Sci. 2021, 129, 110397. [Google Scholar] [CrossRef]
  6. Hu, Z.; Yang, Y.; Liu, L.; Zhou, F. Local flow regime transition criteria of gas-liquid two-phase flow in vertical upward tube with a horizontal rod. Chin. J. Chem. Eng. 2006, 14, 442–449. [Google Scholar] [CrossRef]
  7. Razali, M.A.B.; Xie, C.G.; Loh, W.L. Experimental investigation of gas-liquid flow in a vertical venturi installed downstream of a horizontal blind tee flow conditioner and the flow regime transition. Flow Meas. Instrum. 2021, 80, 101961. [Google Scholar] [CrossRef]
  8. Kim, H.G.; Kim, S.M. Slug flow characteristics of air-liquid two-phase flow in horizontal pipes over a wide range of liquid viscosities. Int. J. Heat Mass Transf. 2023, 208, 124005. [Google Scholar] [CrossRef]
  9. Su, X.; Yao, S.; Song, J.; Wang, Q.; Li, F. The discovery and control of slug flow in coalbed methane reservoir. J. Pet. Sci. Eng. 2019, 172, 115–123. [Google Scholar] [CrossRef]
  10. Cao, Z.; Wu, Z.; Sundén, B. Dimensionless analysis on liquid-liquid flow patterns and scaling law on slug hydrodynamics in cross-junction microchannels. Chem. Eng. J. 2018, 344, 604–615. [Google Scholar] [CrossRef]
  11. Sudarja; Haq, A.; Deendarlianto; Indarto; Widyaparaga, A. Experimental study on the flow pattern and pressure gradient of air-water two-phase flow in a horizontal circular mini-channel. J. Hydrodyn. 2019, 31, 102–116. [Google Scholar] [CrossRef]
  12. Fourar, M.; Bories, S.; Lenormand, R.; Persoff, P. Two-phase flow in smooth and rough fractures: Measurement and correlation by porous-medium and pipe flow models. Water Resour. Res. 1993, 29, 3699–3708. [Google Scholar] [CrossRef]
  13. Wang, Y.; Zhang, Z.; Ranjith, P.G.; Luo, Y. Water-gas flow in rough rock fractures: Insights from coupled triaxial compression experiments. Hydrogeol. J. 2022, 30, 1569–1581. [Google Scholar] [CrossRef]
  14. Chen, C.Y.; Horne, R.N.; Fourar, M. Experimental study of liquid-gas flow structure effects on relative permeabilities in a fracture. Water Resour. Res. 2004, 40, W083011–W0830115. [Google Scholar] [CrossRef]
  15. Suri, Y.; Islam, S.Z.; Hossain, M. Effect of fracture roughness on the hydrodynamics of proppant transport in hydraulic fractures. J. Nat. Gas Sci. Eng. 2020, 80, 103401. [Google Scholar] [CrossRef]
  16. Chen, C.Y.; Horne, R.N. Two-phase flow in rough-walled fractures: Experiments and a flow structure model. Water Resour. Res. 2006, 42, W03430. [Google Scholar] [CrossRef]
  17. Ye, Z.; Liu, H.H.; Jiang, Q.; Zhou, C. Two-phase flow properties of a horizontal fracture: The effect of aperture distribution. Adv. Water Resour. 2015, 76, 43–54. [Google Scholar] [CrossRef]
  18. Cheng, L.; Li, X.; Rong, G.; Zhou, C. The Effect of Surface Wettability and Wall Roughness on the Residual Saturation for the Drainage Process in Sinusoidal Channels. Transp. Porous Media 2019, 129, 203–229. [Google Scholar] [CrossRef]
  19. Zaidani, M.; Hasan, A.; Al-musharfy, M.; Sassi, M. Numerical investigation of surface wettability on gas—Liquid flow hydrodynamics in microchannel: Application to trickle bed reactors. J. Pet. Sci. Eng. 2020, 184, 106576. [Google Scholar] [CrossRef]
  20. Xiao, Y.; Zheng, J.; He, Y.; Wang, L. Droplet and bubble wetting behaviors: The roles of surface wettability and roughness. Colloids Surf. A Physicochem. Eng. Asp. 2022, 653, 130008. [Google Scholar] [CrossRef]
  21. Tan, X.; Konietzky, H.; Chen, W. Numerical Simulation of Heterogeneous Rock Using Discrete Element Model Based on Digital Image Processing. Rock Mech. Rock Eng. 2016, 49, 4957–4964. [Google Scholar] [CrossRef]
  22. Zou, L.; Håkansson, U.; Cvetkovic, V. Two-phase cement grout propagation in homogeneous water-saturated rock fractures. Int. J. Rock Mech. Min. Sci. 2018, 106, 243–249. [Google Scholar] [CrossRef]
  23. Wood, T.; Huang, H. Experimental and Modeling Studies of Episodic Air-Water Two-Phase Flow in Fractures and Fracture Networks. In Fluid Dynamics in Complex Fractured-Porous Systems; John Wiley & Sons: Hoboken, NJ, USA, 2015; pp. 209–228. [Google Scholar]
  24. Xing, K.; Shi, X.; Kang, X.; Qian, J.; Wu, J. Influence of aperture field anisotropy on the immiscible displacement regime transition in rough fractures. J. Hydrol. 2024, 632, 130908. [Google Scholar] [CrossRef]
  25. Wang, Y.; Zhang, Z.; Ranjith, P.G.; Han, X. Flow structure transition and identification of two-phase fluid flow through rough rock fractures. Eur. Phys. J. Plus 2023, 138, 376. [Google Scholar] [CrossRef]
  26. Wang, Q.; Su, X.; Feng, Y.; Wang, H.; Song, J. Experimental study of gas-water two-phase flow patterns in fracture: Implication for enhancing coalbed methane production. J. Pet. Sci. Eng. 2021, 207, 109168. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the geometric model.
Figure 1. Schematic diagram of the geometric model.
Water 16 02860 g001
Figure 2. Comparison of simulation results and experimental results.
Figure 2. Comparison of simulation results and experimental results.
Water 16 02860 g002
Figure 3. Two-phase fluid flow structures for different flow ratios: (a) n = 0.5; (b) n = 1; (c) n = 1.5; and (d) n = 2 (red is water and gas is blue).
Figure 3. Two-phase fluid flow structures for different flow ratios: (a) n = 0.5; (b) n = 1; (c) n = 1.5; and (d) n = 2 (red is water and gas is blue).
Water 16 02860 g003
Figure 4. Two-phase fluid flow structures for different mean contact angles: (a) θ = 30°; (b) θ = 90°; (c) θ = 120°; and (d) θ = 150° (red is water and gas is blue).
Figure 4. Two-phase fluid flow structures for different mean contact angles: (a) θ = 30°; (b) θ = 90°; (c) θ = 120°; and (d) θ = 150° (red is water and gas is blue).
Water 16 02860 g004
Figure 5. Diagram of slug structure variations due to surface wetting heterogeneity.
Figure 5. Diagram of slug structure variations due to surface wetting heterogeneity.
Water 16 02860 g005
Figure 6. Two-phase fluid flow structures for different surface wettability heterogeneities: (a) d = 0°; (b) d = 5°; (c) d = 10°; (d) d = 15°; and (e) d = 20° (red is water and gas is blue).
Figure 6. Two-phase fluid flow structures for different surface wettability heterogeneities: (a) d = 0°; (b) d = 5°; (c) d = 10°; (d) d = 15°; and (e) d = 20° (red is water and gas is blue).
Water 16 02860 g006
Figure 7. Two-phase fluid flow structures for different surface wettability heterogeneities: (a) D = 0.1 mm; (b) D = 0.2 mm; (c) D = 0.3 mm; and (d) D = 0.4 mm (red is water and gas is blue).
Figure 7. Two-phase fluid flow structures for different surface wettability heterogeneities: (a) D = 0.1 mm; (b) D = 0.2 mm; (c) D = 0.3 mm; and (d) D = 0.4 mm (red is water and gas is blue).
Water 16 02860 g007
Figure 8. Schematic diagram of slug flow in a fracture.
Figure 8. Schematic diagram of slug flow in a fracture.
Water 16 02860 g008
Figure 9. Evolution of the ratio of the total length of the slug section to the total length of the flow channel in the fracture under the action of different flow ratios (n): case 1: θ = 30°, D = 0.1 mm; case 2: θ = 90°, D = 0.1 mm; case 3: θ = 120°, D = 0.1 mm; and case 4: θ = 150°, D = 0.1 mm (black for case 1, red for case 2, blue for case 3, and green for case 4).
Figure 9. Evolution of the ratio of the total length of the slug section to the total length of the flow channel in the fracture under the action of different flow ratios (n): case 1: θ = 30°, D = 0.1 mm; case 2: θ = 90°, D = 0.1 mm; case 3: θ = 120°, D = 0.1 mm; and case 4: θ = 150°, D = 0.1 mm (black for case 1, red for case 2, blue for case 3, and green for case 4).
Water 16 02860 g009
Figure 10. The evolution of the relative permeability under different mean contact angles (θ): case 1: n = 0.5, D = 0.1 mm; case 2: n = 1, D = 0.1 mm; case 3: n = 1.5, D = 0.1 mm; and case 4: n = 2, D = 0.1 mm (black for case 1, red for case 2, blue for case 3, and green for case 4).
Figure 10. The evolution of the relative permeability under different mean contact angles (θ): case 1: n = 0.5, D = 0.1 mm; case 2: n = 1, D = 0.1 mm; case 3: n = 1.5, D = 0.1 mm; and case 4: n = 2, D = 0.1 mm (black for case 1, red for case 2, blue for case 3, and green for case 4).
Water 16 02860 g010
Figure 11. The evolution of the ratio of the total length of the slug section to the total length of the flow channel in the fracture under the action of different heterogeneities of surface wettability.
Figure 11. The evolution of the ratio of the total length of the slug section to the total length of the flow channel in the fracture under the action of different heterogeneities of surface wettability.
Water 16 02860 g011
Figure 12. Evolution of the ratio of the total length of the slug section to the total length of the flow channel in the fracture under the action of a different fracture tortuosity: case 1: θ = 30°, σ2 = 0, n = 1; case 2: θ = 90°, σ2 = 0, n = 1.5; case 3: θ = 120°, σ2 = 0, n = 2; and case 4: θ = 150°, σ2 = 0, n = 2.5 (black for case 1, red for case 2, blue for case 3, and green for case 4).
Figure 12. Evolution of the ratio of the total length of the slug section to the total length of the flow channel in the fracture under the action of a different fracture tortuosity: case 1: θ = 30°, σ2 = 0, n = 1; case 2: θ = 90°, σ2 = 0, n = 1.5; case 3: θ = 120°, σ2 = 0, n = 2; and case 4: θ = 150°, σ2 = 0, n = 2.5 (black for case 1, red for case 2, blue for case 3, and green for case 4).
Water 16 02860 g012
Figure 13. Relationship between flow structure and LD.
Figure 13. Relationship between flow structure and LD.
Water 16 02860 g013
Table 1. Physical properties of two-phase fluids.
Table 1. Physical properties of two-phase fluids.
PhaseTemperature (K)Density (kg/m3)Dynamic Viscosity (Pa·s)
Gas293.151.1841.81 × 10−5
Water293.1510001.25
Table 2. The value of the surface contact angle at each stage of fracture.
Table 2. The value of the surface contact angle at each stage of fracture.
No.d = 0 (°)d = 5 (°)d = 10 (°)d = 15 (°)d = 20 (°)
19085807570
29095100105110
Table 3. The values of the correlation coefficient (R2) in Figure 9.
Table 3. The values of the correlation coefficient (R2) in Figure 9.
CaseCase 1Case 2Case 3Case 4
R20.990.970.900.96
Table 4. The values of the correlation coefficient (R2) in Figure 10.
Table 4. The values of the correlation coefficient (R2) in Figure 10.
CaseCase 1Case 2Case 3Case 4
R20.980.950.960.96
Table 5. The values of the correlation coefficient (R2) in Figure 12.
Table 5. The values of the correlation coefficient (R2) in Figure 12.
CaseCase 1Case 2Case 3Case 4
R20.910.920.930.93
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, D.; Pu, H.; Xue, K.; Ni, H. Numerical Simulation of Gas–Water Two-Phase Flow Patterns in Fracture: Implication for Enhancing Natural Gas Production. Water 2024, 16, 2860. https://doi.org/10.3390/w16192860

AMA Style

Liu D, Pu H, Xue K, Ni H. Numerical Simulation of Gas–Water Two-Phase Flow Patterns in Fracture: Implication for Enhancing Natural Gas Production. Water. 2024; 16(19):2860. https://doi.org/10.3390/w16192860

Chicago/Turabian Style

Liu, Dejun, Hai Pu, Kangsheng Xue, and Hongyang Ni. 2024. "Numerical Simulation of Gas–Water Two-Phase Flow Patterns in Fracture: Implication for Enhancing Natural Gas Production" Water 16, no. 19: 2860. https://doi.org/10.3390/w16192860

APA Style

Liu, D., Pu, H., Xue, K., & Ni, H. (2024). Numerical Simulation of Gas–Water Two-Phase Flow Patterns in Fracture: Implication for Enhancing Natural Gas Production. Water, 16(19), 2860. https://doi.org/10.3390/w16192860

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop