Research Article | Open Access

Lifeng Zeng, Longtan Shao, Xiaoxia Guo, "Seepage Force and Its Direct Mechanical Effects in Hydrate-Bearing Porous Media", *Geofluids*, vol. 2021, Article ID 5569665, 15 pages, 2021. https://doi.org/10.1155/2021/5569665

# Seepage Force and Its Direct Mechanical Effects in Hydrate-Bearing Porous Media

**Academic Editor:**Timothy S. Collett

#### Abstract

The direct mechanical effects of seepage force on the behavior of hydrate-bearing porous media (HBPM) have often been neglected in previous studies, which may lead to inaccurate predictions of the mechanical behavior of HBPM under seepage conditions. Here, we propose an extended three-phase physical model for unsaturated HBPM, including gas, capillary water, and generalized solid skeleton (GSS). Based on the model, the force balance equations for the three phases are formulated. Performing a force analysis of generalized solid particles under seepage conditions, we find that the tangential force acting on the generalized solid particle is the seepage force in HBPM. By combining this finding with the formulated balance equations, we derive the expression for the seepage force, which can distinguish the mechanical effects between the tangential force and normal force. The stresses, induced by the tangential force (i.e., seepage force), are divided into two types: one acts on the cross sections of generalized solid particles and the other on the contacts between particles. Neither of them is transmitted through the GSS. The former mainly causes the particles themselves to compress, whereas the latter primarily influences the sliding of the particles at contacts. Based on the mechanisms of these two stresses, their effects on the mechanical behavior of HBPM are quantified, which provides a new insight to evaluate the direct mechanical effects of seepage force.

#### 1. Introduction

Seepage is an important flow phenomenon in HBPM, and it refers to the flow of water in saturated HBPM or of gas and water in unsaturated HBPM under the gradient of total head [1]. In saturated or unsaturated HBPM, the characteristics of seepage largely determine the magnitude of fluid flow velocity [2], the fluid pressure and saturation distribution [3, 4], and the stress and strain distribution of HBPM [5], which further affects the stability of hydrate-bearing sediments [6, 7]. As a result, the analysis of the seepage characteristics of HBPM has become an active field in the efficient and safe exploitation of gas hydrate.

Some researchers have performed a large number of experiments by macroscopic and microscopic techniques in the laboratory to investigate the seepage characteristics of HBPM, such as the variation of gas/water relative permeability with different hydrate habits and saturations [8–10], the impact of various pressure gradients on gas production rate and water relative permeability [11, 12], and the evolution mechanism of relative permeability associated with capillarity and drying-wetting cycles [13, 14]. On the other hand, many numerical simulations have been conducted to evaluate the influence of different factors on the gas/water production rate during hydrate exploitation, such as the permeability of host porous media [15, 16], the methods used to cause hydrate to dissociate [17–19], and the heterogeneity of hydrate saturation [20]. It should be noted that whether the studies were performed by experimental tests or by numerical simulations, they mainly focused on the flow behavior of fluids (i.e., gas and water) in HBPM. The flow behavior, to some extent, reflects the action of the solid phase (i.e., mineral particle and hydrate) in HBPM on the fluids. However, the corresponding reaction is less considered during hydrate exploitation.

The reaction mentioned above represents the force applied by the flowing fluids to the solid phase, usually referred to as seepage force. The mechanical effects of seepage force are generally considered in such a way that the reaction of seepage force results in the variation of pore-fluid pressure, which causes the variation in the effective stress and mechanical response of HBPM [5, 21–23]. In fact, the way of considering seepage force only involves its indirect mechanical effects but ignores its direct mechanical effects.

The objective of this work is to identify the mechanism of seepage force in HBPM and to quantify the seepage force and its direct effects on the mechanical behavior of HBPM. For this objective, we first analyze the interactions between hydrates, solid particles, gas, and water and build an extended three-phase physical model including the gas phase, capillary water phase, and GSS phase. Based on this model, the independent force balance equations for these three component phases are formulated. Then, a force analysis considering the fluid flow at particle scale is performed, the mechanism of the seepage force is identified, and the expression for the seepage force is derived. Finally, the stresses, induced by the seepage force, acting on the GSS are analyzed, and their direct effects on the mechanical behavior of HBPM are quantified.

#### 2. An Extended Three-Phase Physical Model for HBPM

HBPM is generally composed of hydrate, gas, capillary water, and solid particles, particularly during the process of hydrate dissociation through depressurization, thermal stimulation, or CO_{2} injection [24–28]. The location of hydrate formation depends largely on environmental conditions, such as gas source, water saturation, and the type of host porous media [29, 30]. Different environmental conditions lead to different pore-scale habits of hydrate in HBPM. The habits of hydrate are commonly grouped into four categories: cementation, particle coating, pore filling, and load bearing [29, 31–33]. The cementing hydrates form at interparticle contacts and bond neighboring particles together. In this case, even a small number of hydrates can remarkably enhance the strength and stiffness of HBPM [34–36]; owing to the hydrophilic nature of host porous media, the particle-coating hydrates crystallize on the whole surface of particles. The hydrates of this category can, to a certain extent, bond some particles and thus give rise to a restriction of the movement of these particles [37, 38]; the pore-filling hydrates crystallize on the partial surface of particles and enlarge freely into the middle of pores without bonding adjacent particles together. Compared with the cementing hydrates, the pore-filling hydrates are believed to have a relatively small influence on the shear strength and stiffness of HBPM, but they can still restrict the movement of surrounding particles [37, 39]; the load-bearing hydrates evolving from the pore-filling hydrates bond particles together into a cluster structure when hydrate saturation is beyond 40% [40]. In particular, the hydrate saturation can be higher than 80% in some local regions where the hydrates may directly bear loads [41, 42]. Although hydrates with different pore-scale habits have different contributions to the shear strength and stiffness of HBPM, all of them can be regarded as a structure to bear and transfer loads. Additionally, the interactions between gas and capillary water create gas-water interfaces. The interfaces have the similar mechanical behavior to an elastic contractile skin and, especially those occurring at particle contacts, can pull solid particles together and thus help to increase the strength and stiffness of porous media [43, 44]. Hence, the gas-water interfaces can also be regarded as another structure to bear and transfer loads. In addition to the hydrates and gas-water interfaces, the solid particle skeleton in HBPM that is formed by the frictional and interlocking mechanisms among particles is the key structure to bear and transfer loads [45, 46]. Considering that the hydrates, gas-water interfaces, and solid particle skeleton all can bear and transfer loads, they could be combined into a generalized structure, called GSS. In such a case, the GSS, gas, and capillary water constitute an extended three-phase physical model (Figure 1) that can be employed to describe the mechanical behavior of HBPM. For simplicity, the case of particle-coating hydrates is only illustrated.

Based on the extended three-phase physical model presented above, a series of definitions can be given to describe the volume-saturation properties of HBPM. For an RVE of HBPM, the following equation must be satisfied (Figure 2): where , , and are the volumes of gas, capillary water, and GSS, respectively, and is the total volume of RVE.

The volume fraction of each phase is defined as the ratio of the volume of the corresponding phase to the total volume of the RVE: where , , and are the volume fractions of gas, capillary water, and GSS, respectively, as illustrated in Figure 2. Combining equations (1), (2), (3), and (4), we have The pore spaces occupied by the gas and by the capillary water constitute the effective pore space of HBPM, and the following equation can be given: where is the effective pore volume.

Based on the above definition, the volume fraction of the effective pore space is then The proportion of the effective pore space occupied by the capillary water is expressed as the effective saturation of the capillary water (Figure 2):

Similarly, the area fraction for each phase is defined as the ratio of the area of the corresponding phase to the total area of the RVE. In general, the area fraction is postulated to be equal to the volume fraction in the case of homogeneous porous media [47].

Based on the analysis of the pore-scale habits of hydrate and the interactions between the gas and capillary water, an extended three-phase physical model for HBPM is well established, as shown in Figure 1. The model provides a physical basis for the derivation of the force balance equations for three individual phases in HBPM.

#### 3. Force Balance Equations for Three Individual Phases in HBPM

The force balance equations for three individual phases in HBPM are the prerequisites for determining the seepage force caused by the flowing fluids and its direct effects on the mechanical behavior of HBPM. The force balance equations for individual phases in unsaturated soils have been formulated by the differential element-based force balance method [48–50]. This method is based on the widely accepted concepts of force balance and can explicitly take the interaction forces between different phases into account. Therefore, this method will be adopted to formulate the force balance equations for the gas, capillary water, and GSS contained in HBPM.

##### 3.1. Force Balance Equation for the Gas Phase

Formulating the force balance equation for the gas phase needs to calculate the surface forces, gravitational force, and interaction force that act on the differential element of the gas phase. The surface forces are calculated as the gas pressure times the area occupied by the gas. The gravitational force is equal to the product of the unit weight of the gas and the corresponding volume. The interaction force is expressed as the interaction force per unit volume times the corresponding volume. For simplicity, only the forces in the -direction are shown in Figure 3.

The equilibrium condition for the differential element of the gas phase requires that the resultant force should vanish. Summing the force components in the -direction yields the balance equation for the gas phase: where is the gas pressure, is the interaction force per unit volume between the gas and the GSS; is the density of gas, and is the gravitational acceleration.

##### 3.2. Force Balance Equation for the Capillary Water Phase

In analogy with the force analysis of the gas phase, the force components acting on the differential element of the capillary water phase in the -direction are shown in Figure 4.

The summation of the force components in the -direction gives the balance equation for the capillary water phase: where is the capillary water pressure, is the interaction force per unit volume between the capillary water and the GSS, and is the density of water.

##### 3.3. Force Balance Equation for the GSS Phase

Compared with the differential elements of the gas phase and capillary water phase that are only subjected to one type of surface force, the differential element of the GSS phase is subjected to two types of surface force. One represents the forces exerted by the adjacent differential elements of the GSS phase, and the other refers to the forces induced by the gas pressure and capillary water pressure. As the gas and capillary water together occupy the effective pore space, a difficult problem encountered is how to determine the magnitudes and corresponding areas of the stresses acting on the differential element of the GSS induced by the second type of surface force.

To solve the above problem, we shall first consider a saturated system composed of two idealized spherical solid particles coated with hydrates in equilibrium under hydrostatic pressure, as illustrated in Figure 5.

According to the equilibrium condition for section 1-1 (Figure 5(b)), we can obtain where is the stress acting on the cross section of the solid particle induced by the hydrostatic pressure and is the cross-sectional area of the solid particle.

Similarly, using the equilibrium condition for section 2-2 (Figure 5(c)), we can get where is the stress acting on the contact between particles induced by the hydrostatic pressure and is the contact area.

In analogy with the saturated system, an unsaturated system is analyzed that includes two idealized spherical generalized solid particles (composed of solid particles, particle-coating hydrates, and gas-water interfaces) in equilibrium under the gas pressure and capillary water pressure, as shown in Figure 6. In order to apply the approach developed in the saturated system to the unsaturated system, a basic assumption analogous to that adopted by the theory of mixtures [51] is made here. It is assumed that the gas and capillary water independently fill the total pore space of the unsaturated system according to their respective volume fractions. In this case, this system can be regarded as two subsystems. One refers to a system where the pore space is fully filled with the gas under a homogenized gas pressure of (Figure 6(b)); the other corresponds to a system where the capillary water under a homogenized water pressure of completely occupies the pore space (Figure 6(c)). From Figure 6(b), on the basis of the equilibrium condition of section 3-3, the following equation can be obtained: where denotes the stress acting on the cross section of the generalized solid particle induced by the homogenized gas pressure.

Likewise, according to the equilibrium condition of section 4-4, the following equation can be obtained: where denotes the stress acting on the contact between generalized solid particles induced by the homogenized gas pressure.

From Figure 6(c), using the equilibrium conditions of sections 5-5 and 6-6, equations (15) and (16) can be given, respectively: where denotes the stress acting on the cross section of the generalized solid particle induced by the homogenized capillary water pressure. where denotes the stress acting on the contact between generalized solid particles induced by the homogenized capillary water pressure.

Extending the analysis results obtained in the unsaturated system to a differential element of the GSS in HBPM, we can easily compute the magnitudes and corresponding areas of the stresses acting on the differential element induced by the homogenized gas pressure and by the homogenized capillary water pressure. These stresses are illustrated in Figure 7. Besides, Figure 7 also illustrates the gravitational force, interaction forces per unit volume, and the stresses exerted by adjacent differential elements.

Summing the force components in the -direction yields the balance equation for the GSS phase: where is the shear stress acting on the -plane in the -direction, is the normal stress acting on the -plane, is the shear stress acting on the -plane in the -direction, is the interaction force per unit volume between the GSS and the gas, is the interaction force per unit volume between the GSS and the capillary water, and is the density of the solid particle. The detailed derivation of equation (17) is presented in Appendix A.

#### 4. Identification and Quantification of Seepage Force in HBPM

When the gas and capillary water flow in the pore space of HBPM, they are subjected to the seepage resistance exerted by the GSS. On the contrary, the GSS are subjected to the seepage force imposed by the flowing gas and capillary water. To identify the seepage resistance and seepage force in unsaturated HBPM, it is necessary to make an identification of the seepage resistance and seepage force in saturated HBPM.

##### 4.1. Seepage Force in Saturated HBPM

A saturated system, which is composed of two idealized spherical solid particles coated with hydrates, is used to identify the seepage resistance and seepage forces under saturated seepage conditions, as illustrated in Figure 8. It can be seen from Figure 8(a) that the solid particles are subjected to the normal and tangential forces exerted by the flowing water. Conversely, the flowing water is also subjected to the normal and tangential forces exerted by the solid particles (Figure 8(b)). The normal forces acting on the flowing water generate pore-water pressure and thus increase the pressure potential of pore water, whereas the tangential forces acting on the water hinder the flow of water and thus reduce the pressure potential of pore water. As such, the tangential forces acting on the flowing water are the seepage resistance exerted by the solid particles. On the other hand, the normal forces acting on the solid particles cause the volumetric strain of particles, but the tangential forces acting on the solid particles promote the movement of particles. Hence, the tangential forces acting on the solid particles are the seepage forces exerted by the flowing water. Since the tangential forces acting on the flowing water and on the solid particles have the same value but opposite directions, the seepage resistance and the seepage force are numerically equal but opposite in direction.

**(a)**

**(b)**

Extending the results from a particle-level analysis to a differential element-level analysis allows us to calculate the seepage resistance and seepage force in saturated HBPM. In this case, according to equation (10), the force balance equation for the water phase in saturated HBPM can be given by where is the interaction force per unit volume between the water and the GSS.

When the water flows in the pore space of the differential element of saturated HBPM, it is subjected to the normal and tangential forces exerted by the GSS. The effects of the normal and tangential forces acting on the water phase correspond to the first and second terms on the left side of equation (18), respectively. Therefore, the interaction force per unit volume acting on the water phase is the seepage resistance in saturated HBPM. By substituting the relation between pore-water pressure and pressure head, , equation (18) can be rewritten as where is the pressure head of water and is the elevation head.

Then, the seepage resistance in saturated HBPM in the -direction can be expressed as where is the gradient of the total head of water in the -direction.

Similarly, based on equation (17), the force balance equation for the GSS phase in saturated HBPM can be obtained: where is the interaction force per unit volume between the GSS and the water and .When the water flows, the GSS is subjected to the normal and tangential forces exerted by the flowing water. The effects of the normal and tangential forces acting on the GSS associate with the fourth and fifth terms on the left side of equation (21), respectively. Hence, the interaction force per unit volume acting on the GSS is the seepage force in saturated HBPM. As demonstrated earlier, the seepage force and the seepage resistance are numerically equal, and the seepage force in saturated HBPM can therefore be expressed as where is the seepage force in the -direction acting on the GSS.

Generalizing the result obtained in the -direction to the - and -directions, we have where is the seepage force vector acting on the GSS exerted by the flowing water and is the gradient of the total head of water.

##### 4.2. Seepage Force in Unsaturated HBPM

Applying the analytical approach developed in saturated HBPM to unsaturated HBPM, we can compute the seepage resistance and seepage force in unsaturated HBPM. According to equations (9) and (10), the seepage resistance in unsaturated HBPM can be given by equations (24) and (25), respectively: where is the seepage resistance in the -direction acting on the gas phase, is the pressure head of gas and , and is the gradient of the total head of gas in the -direction. where is the seepage resistance in the -direction acting on the capillary water phase, is the pressure head of capillary water and , and is the gradient of the total head of capillary water in the -direction.

When both the gas and the capillary water flow, the GSS is subjected to the normal and tangential forces exerted by the flowing gas and capillary water. The effects of the normal and tangential forces acting on the GSS exerted by the gas correspond to the fourth and sixth terms on the left side of equation (17), respectively. In equation (17), the fifth and seventh terms on the left side are attributed to the effects of the normal and tangential forces acting on the GSS exerted by the capillary water, respectively. Therefore, the interaction forces per unit volume, and , acting on the GSS are the seepage forces exerted by the flowing gas and by the flowing capillary water, respectively. Performing similar analyses in the - and -directions and then using the conclusion drawn above that the seepage force is numerically equal to the seepage resistance, we can give the seepage forces in unsaturated HBPM as follows: where is the seepage force vector acting on the GSS exerted by the flowing gas and is the gradient of the total head of gas. where is the seepage force vector acting on the GSS exerted by the flowing capillary water and is the gradient of the total head of capillary water.

#### 5. Mechanical Effects of the Seepage Force in HBPM

The mechanical effects of the seepage force in HBPM depend on the transfer mechanism of the stresses acting on the GSS induced by the seepage force and the magnitude of these stresses. In order to identify the transfer mechanism, we need to turn to the effective stress equation for HBPM.

##### 5.1. Effective Stress Equation for Unsaturated HBPM

The effective stress equation for unsaturated HBPM can be derived based on four force balance equations for the gas phase, the capillary water phase, the GSS phase, and the whole HBPM. As the first three force balance equations have been formulated, i.e., equations (9), (10), and (17), the total force balance equation for the whole HBPM will be formulated next. Taking a differential element from HBPM and then analyzing the forces acting on it in the -direction (Figure 9), we can obtain the total balance equation for the whole HBPM: where is the total shear stress acting on the -plane in the -direction, is the total normal stress acting on the -plane, is the total shear stress acting on the -plane in the -direction, and is the density of HBPM and .

The superposition of equations (9), (10), and (17) for the gas phase, capillary water phase, and GSS phase, respectively, leads to the total balance equation for the whole HBPM: Considering that both the gas and the capillary water are unable to carry shear stress, the identities and hold. In this case, equation (29) can be rewritten as Comparing equation (30) with equation (28) and then generalizing the result obtained in the -direction to the - and -directions, we get By substituting equations (2), (3), (6), and (8) into equation (31), the effective stress equation for unsaturated HBPM can be derived: where is the effective stress, is the total stress, is the neutral stress, and is the effective saturation of the capillary water, which varies from 0 to 1.

Many studies [49, 50, 52–54] have shown that the effective stress is transmitted through the solid skeleton and largely controls the shear strength and volumetric strain behavior of soils. This statement is conducive to the identification of the transfer mechanism of the seepage force.

##### 5.2. Mechanical Effects of the Seepage Force in Saturated HBPM

According to the balance equation (equation (21)) for the GSS in saturated HBPM, it is found that there are three stresses acting on the GSS, i.e., the effective stress, the stress induced by normal forces, and the stress due to tangential forces (or seepage forces). As the effective stress is transmitted through the GSS and governs largely the mechanical behavior of HBPM, the stresses induced by normal forces and by seepage forces are not transmitted through the GSS. After identifying the transfer mechanism of the stress induced by seepage forces, the magnitude of this stress needs to be determined. We now consider a saturated system consisting of two idealized spherical solid particles coated with hydrates under saturated seepage conditions, as illustrated in Figure 10. In this figure, only the seepage forces are illustrated for simplicity because only the mechanical effects of the seepage forces will be analyzed.

For a particle-scale system (Figure 10), the seepage forces acting on the solid particle can be combined into a resultant force. By considering equations (18) and (23), this resultant force can be computed as where is the resultant force acting on the solid particle induced by the seepage forces and is the total area of the system (Figure 10(a)). In this system, , and is the cross-sectional area of the solid particle.

According to the equilibrium condition for section 1-1 (Figure 10(b)), we get where is the stress acting on the cross section of the solid particle induced by the seepage forces.

As is not transmitted through the GSS, it only produces the volumetric strain of the solid particle. Then, the contribution of to the volumetric strain of saturated HBPM can be calculated as where is the volumetric strain contributed by induced by the seepage forces under saturated seepage conditions and is the compressibility of the solid particle.

Likewise, using the equilibrium condition for section 2-2 (Figure 10(c)), we have where is the stress acting on the contact between solid particles induced by the seepage forces and is the contact area.

Owing to not being transmitted through the GSS, it only affects the sliding of solid particles. Then, the contribution of to the shear strength of saturated HBPM can be computed as where is the shear strength contributed by induced by the seepage forces under saturated seepage conditions, is the contact area ratio and , and is the intrinsic friction angle of the solid particle.

##### 5.3. Mechanical Effects of the Seepage Force in Unsaturated HBPM

The mechanical effect of seepage force in unsaturated HBPM can be quantified by employing the approach developed in saturated HBPM. We now consider an unsaturated system consisting of two idealized spherical generalized solid particles under unsaturated seepage conditions, as illustrated in Figure 11. In order to employ the approach developed in saturated HBPM, an assumption analogous to that made in Figure 6 needs to be introduced here. In this case, the unsaturated system can be considered two subsystems. One represents a system saturated with gas under a homogenized gas pressure of (Figure 11(b)); the other refers to a system saturated with capillary water under a homogenized water pressure of (Figure 11(c)).

In the light of the results obtained in the saturated system, we can easily compute the seepage forces acting on the generalized solid particle in the unsaturated system. For a subsystem fully saturated with gas (Figure 11(b)), the seepage forces acting on the particle can be composed into a resultant force. This resultant force can be expressed as where is the resultant force acting on the particle induced by the seepage forces in a subsystem fully saturated with gas.

From Figure 11(b), according to the equilibrium conditions of sections 3-3 and 4-4, the following two equations can be obtained: where and are the stresses acting on the cross section of the generalized solid particle and on the contact between particles, respectively, induced by the seepage forces in a subsystem fully saturated with gas.

Similarly, for a subsystem completely saturated with capillary water (see Figure 11(c)), the seepage forces acting on the generalized solid particle can be combined into a resultant force. This resultant force can be expressed as where is the resultant force acting on the particle induced by the seepage forces in a subsystem completely saturated with capillary water.

From Figure 11(c), using the equilibrium conditions of sections 5-5 and 6-6, equations (41) and (42) can be given, respectively: where and denote the stresses acting on the cross section of the generalized solid particle and on the contact between particles, respectively, induced by the seepage forces in a subsystem completely saturated with capillary water.

As and are not transmitted through the GSS, they only produce the volumetric strain of the solid particle. Then, the contribution of and induced by the seepage forces to the volumetric strain of unsaturated HBPM can be calculated as where is the volumetric strain contributed by that is induced by the seepage forces under unsaturated seepage conditions.

Because of and not being transmitted through the GSS, they only influence the sliding of generalized solid particles. Then, the contribution of and induced by the seepage forces to the shear strength of unsaturated HBPM can be computed as where is the shear strength contributed by that is induced by the seepage forces under unsaturated seepage conditions.

#### 6. Conclusions

By analyzing the interactions between hydrates, solid particles, gas, and water in HBPM, an extended three-phase physical model is presented, including the gas phase, capillary water phase, and GSS phase. Based on the model, the independent force balance equations for these three phases are formulated, which provide a sound theoretical basis for the derivation of seepage force.

Under saturated or unsaturated seepage conditions, the forces acting on the fluids and on the solid particles are composed of normal and tangential forces. For the fluids, the normal forces generate pore-fluid pressure and thus increase the pressure potential of fluids, whereas the tangential forces restrict the flow of fluids and thus reduce the pressure potential of fluids. Hence, the tangential forces acting on the fluids are the seepage resistance. On the contrary, for the generalized solid particles, the normal forces produce the volumetric strain of particles, while the tangential forces drive the particles to move. Therefore, the tangential forces that act on the generalized solid particles are the seepage forces. According to this finding, for saturated HBPM, the seepage force is expressed as equation (23); for unsaturated HBPM, the seepage forces induced by the gas and capillary water are calculated by equations (26) and (27), respectively. Compared with the equations for seepage force derived by previous researchers [55, 56], the equations obtained in this study distinguish the mechanical influence of the tangential force from that of the normal force, which ensures that the seepage force and seepage resistance are equal in magnitude but opposite in direction.

The seepage forces cause the stresses acting on the cross section of the generalized solid particle and on the contact between particles. Based on the effective stress equation (equation (32)) and the balance equations (equations (21) and (17)) for the GSS, it is found that the stresses acting on the cross section of the generalized solid particle and on the contact between particles are not transmitted through the GSS. The former mainly produces the volumetric strain of particles, whereas the latter primarily influences the sliding of the particle at contacts. Therefore, the mechanical effects of the stresses caused by the seepage forces can be evaluated by equations (35) and (37) for saturated HBPM and by equations (43) and (44) for unsaturated HBPM. These equations allow us to make a direct evaluation of the mechanical effects of seepage force.

#### Appendix

#### Force Balance Equation for the GSS Phase

The equilibrium condition for the differential element of the GSS phase demands that the resultant force should vanish, as shown in Figure 7. Summing the force components in the -direction yields where and denote the stresses acting on the cross sections of generalized solid particles and on the contacts between generalized solid particles of the upper surface of the differential element, respectively, induced by the homogenized gas pressure (see Figure 7), and denote the stresses acting on the cross sections of generalized solid particles and on the contacts between generalized solid particles of the lower surface, respectively, induced by the homogenized gas pressure, is the corresponding area of and , is the corresponding area of and , and denote the stresses acting on the cross sections of generalized solid particles and on the contacts between generalized solid particles of the upper surface of the differential element, respectively, induced by the homogenized capillary water pressure, and denote the stresses acting on the cross sections of generalized solid particles and on the contacts between generalized solid particles of the lower surface, respectively, induced by the homogenized capillary water pressure, is the corresponding area of and , and is the corresponding area of and .

Based on the analysis results of Figure 6, we can obtain

By expanding and simplifying equation (A.1), we have Substituting equation (5) into equation (A.3), we can obtain equation (17).

#### Data Availability

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

#### Conflicts of Interest

The authors state that they do not have any financial or nonfinancial conflict of interests.

#### Acknowledgments

This study was supported by the National Natural Science Foundation of China (Grant No. 52079018).

#### References

- W. Liu, Z. Wu, J. Li, J. Zheng, and Y. Li, “The seepage characteristics of methane hydrate-bearing clayey sediments under various pressure gradients,”
*Energy*, vol. 191, article 116507, 2020. View at: Publisher Site | Google Scholar - S. Wang, P. Wang, B. Chen, M. Yang, and Y. Li, “Velocity mapping of steady water flow through methane hydrate bearing samples,”
*Journal of Natural Gas Science and Engineering*, vol. 53, pp. 385–393, 2018. View at: Publisher Site | Google Scholar - X. Sun, Y. Li, Y. Liu, and Y. Song, “The effects of compressibility of natural gas hydrate-bearing sediments on gas production using depressurization,”
*Energy*, vol. 185, pp. 837–846, 2019. View at: Publisher Site | Google Scholar - J. Rutqvist, G. J. Moridis, T. Grover, S. Silpngarmlert, T. S. Collett, and S. A. Holdich, “Coupled multiphase fluid flow and wellbore stability analysis associated with gas production from oceanic hydrate-bearing sediments,”
*Journal of Petroleum Science and Engineering*, vol. 92-93, pp. 65–81, 2012. View at: Publisher Site | Google Scholar - X. Sun, L. Wang, H. Luo, Y. Song, and Y. Li, “Numerical modeling for the mechanical behavior of marine gas hydrate-bearing sediments during hydrate production by depressurization,”
*Journal of Petroleum Science and Engineering*, vol. 177, pp. 971–982, 2019. View at: Publisher Site | Google Scholar - X. Wang, B. Sun, H. Gao, Z. Wang, H. Li, and Y. Chen, “Numerical simulation of the stability of hydrate layer during well cementing in deep-water region,”
*Journal of Petroleum Science and Engineering*, vol. 176, pp. 893–905, 2019. View at: Publisher Site | Google Scholar - J. Rutqvist and G. J. Moridis, “Numerical studies on the geomechanical stability of hydrate-bearing sediments,”
*SPE Journal*, vol. 14, no. 2, pp. 267–282, 2009. View at: Publisher Site | Google Scholar - S. Dai and Y. Seol, “Water permeability in hydrate-bearing sediments: a pore-scale study,”
*Geophysical Research Letters*, vol. 41, no. 12, pp. 4176–4184, 2014. View at: Publisher Site | Google Scholar - P. Wu, Y. Li, W. Liu, Y. Liu, D. Wang, and Y. Song, “Microstructure evolution of hydrate-bearing sands during thermal dissociation and ensued impacts on the mechanical and seepage characteristics,”
*Journal of Geophysical Research: Solid Earth*, vol. 125, no. 5, 2020. View at: Publisher Site | Google Scholar - M. L. Delli and J. L. H. Grozic, “Experimental determination of permeability of porous media in the presence of gas hydrates,”
*Journal of Petroleum Science and Engineering*, vol. 120, pp. 1–9, 2014. View at: Publisher Site | Google Scholar - J. Ding, S. Yang, X. Nie, and Z. Wang, “Dynamic threshold pressure gradient in tight gas reservoir,”
*Journal of Natural Gas Science and Engineering*, vol. 20, pp. 155–160, 2014. View at: Publisher Site | Google Scholar - C. Lu, Y. Xia, X. Sun et al., “Permeability evolution at various pressure gradients in natural gas hydrate reservoir at the Shenhu area in the South China Sea,”
*Energies*, vol. 12, no. 19, article 3688, 2019. View at: Publisher Site | Google Scholar - G. Li, L. Zhan, T. Yun, and S. Dai, “Pore-scale controls on the gas and water transport in hydrate-bearing sediments,”
*Geophysical Research Letters*, vol. 47, no. 12, 2020. View at: Publisher Site | Google Scholar - R. L. Kleinberg, C. Flaum, D. D. Griffin et al., “Deep sea NMR: methane hydrate growth habit in porous media and its relationship to hydraulic permeability, deposit accumulation, and submarine slope stability,”
*Journal of Geophysical Research: Solid Earth*, vol. 108, no. B10, 2003. View at: Publisher Site | Google Scholar - X. Sun, T. Luo, L. Wang, H. Wang, Y. Song, and Y. Li, “Numerical simulation of gas recovery from a low-permeability hydrate reservoir by depressurization,”
*Applied Energy*, vol. 250, pp. 7–18, 2019. View at: Publisher Site | Google Scholar - X. Sun, H. Luo, T. Luo, Y. Song, and Y. Li, “Numerical study of gas production from marine hydrate formations considering soil compression and hydrate dissociation due to depressurization,”
*Marine and Petroleum Geology*, vol. 102, pp. 759–774, 2019. View at: Publisher Site | Google Scholar - G. Ahmadi, C. Ji, and D. H. Smith, “Numerical solution for natural gas production from methane hydrate dissociation,”
*Journal of Petroleum Science and Engineering*, vol. 41, no. 4, pp. 269–285, 2004. View at: Publisher Site | Google Scholar - B. Wang, Z. Fan, J. Zhao, X. Lv, W. Pang, and Q. Li, “Influence of intrinsic permeability of reservoir rocks on gas recovery from hydrate deposits via a combined depressurization and thermal stimulation approach,”
*Applied Energy*, vol. 229, pp. 858–871, 2018. View at: Publisher Site | Google Scholar - W. Sung, H. Lee, H. Lee, and C. Lee, “Numerical study for production performances of a methane hydrate reservoir stimulated by inhibitor injection,”
*Energy Sources*, vol. 24, no. 6, pp. 499–512, 2002. View at: Publisher Site | Google Scholar - J. Sun, L. Zhang, F. Ning et al., “Production potential and stability of hydrate-bearing sediments at the site GMGS3-W19 in the South China Sea: a preliminary feasibility study,”
*Marine and Petroleum Geology*, vol. 86, pp. 447–473, 2017. View at: Publisher Site | Google Scholar - X. Sun, H. Luo, and K. Soga, “A coupled thermal–hydraulic–mechanical–chemical (THMC) model for methane hydrate bearing sediments using COMSOL Multiphysics,”
*Journal of Zhejiang University-Science A*, vol. 19, no. 8, pp. 600–623, 2018. View at: Publisher Site | Google Scholar - L. Li, X. Li, Y. Wang et al., “Investigating the interaction effects between reservoir deformation and hydrate dissociation in hydrate-bearing sediment by depressurization method,”
*Energies*, vol. 14, no. 3, p. 548, 2021. View at: Publisher Site | Google Scholar - M. De La Fuente, J. Vaunat, and H. Marín-Moreno, “Thermo-hydro-mechanical coupled modeling of methane hydrate-bearing sediments: formulation and application,”
*Energies*, vol. 12, no. 11, article 2178, 2019. View at: Publisher Site | Google Scholar - T. Luo, Y. Li, B. N. Madhusudhan, X. Sun, and Y. Song, “Deformation behaviors of hydrate-bearing silty sediment induced by depressurization and thermal recovery,”
*Applied Energy*, vol. 276, article 115468, 2020. View at: Publisher Site | Google Scholar - Y. Song, Y. Kuang, Z. Fan, Y. Zhao, and J. Zhao, “Influence of core scale permeability on gas production from methane hydrate by thermal stimulation,”
*International Journal of Heat and Mass Transfer*, vol. 121, pp. 207–214, 2018. View at: Publisher Site | Google Scholar - J. W. Jung, D. N. Espinoza, and J. C. Santamarina, “Properties and phenomena relevant to CH4-CO2 replacement in hydrate-bearing sediments,”
*Journal of Geophysical Research*, vol. 115, no. B10, 2010. View at: Publisher Site | Google Scholar - Y. Song, C. Cheng, J. Zhao et al., “Evaluation of gas production from methane hydrates using depressurization, thermal stimulation and combined methods,”
*Applied Energy*, vol. 145, pp. 265–277, 2015. View at: Publisher Site | Google Scholar - Y. Song, X. Wang, M. Yang et al., “Study of selected factors affecting hydrate-based carbon dioxide separation from simulated fuel gas in porous media,”
*Energy & Fuels*, vol. 27, no. 6, pp. 3341–3348, 2013. View at: Publisher Site | Google Scholar - L. Lei, Y. Seol, J.-H. Choi, and T. J. Kneafsey, “Pore habit of methane hydrate and its evolution in sediment matrix – laboratory visualization with phase-contrast micro-CT,”
*Marine and Petroleum Geology*, vol. 104, pp. 451–467, 2019. View at: Publisher Site | Google Scholar - J. A. Priest, E. V. L. Rees, and C. R. I. Clayton, “Influence of gas hydrate morphology on the seismic velocities of sands,”
*Journal of Geophysical Research: Solid Earth*, vol. 114, no. B11, 2009. View at: Publisher Site | Google Scholar - W. F. Waite, J. C. Santamarina, D. D. Cortes et al., “Physical properties of hydrate-bearing sediments,”
*Reviews of Geophysics*, vol. 47, no. 4, 2009. View at: Publisher Site | Google Scholar - S. Dai, J. C. Santamarina, W. F. Waite, and T. J. Kneafsey, “Hydrate morphology: physical properties of sands with patchy hydrate saturation,”
*Journal of Geophysical Research: Solid Earth*, vol. 117, no. B11, 2012. View at: Publisher Site | Google Scholar - P. Wu, Y. Li, X. Sun, W. Liu, and Y. Song, “Mechanical characteristics of hydrate-bearing sediment: a review,”
*Energy & Fuels*, vol. 35, no. 2, pp. 1041–1057, 2021. View at: Publisher Site | Google Scholar - J. A. Priest, A. I. Best, and C. R. I. Clayton, “A laboratory investigation into the seismic velocities of methane gas hydrate-bearing sand,”
*Journal of Geophysical Research: Solid Earth*, vol. 110, no. B4, 2005. View at: Publisher Site | Google Scholar - W. F. Waite, W. J. Winters, and D. H. Mason, “Methane hydrate formation in partially water-saturated Ottawa sand,”
*American Mineralogist*, vol. 89, no. 8-9, pp. 1202–1207, 2004. View at: Publisher Site | Google Scholar - M. Hyodo, Y. Li, J. Yoneda et al., “Mechanical behavior of gas-saturated methane hydrate-bearing sediments,”
*Journal of Geophysical Research: Solid Earth*, vol. 118, no. 10, pp. 5185–5194, 2013. View at: Publisher Site | Google Scholar - P. Wu, Y. Li, W. Liu, X. Sun, X. Kong, and Y. Song, “Cementation failure behavior of consolidated gas hydrate-bearing sand,”
*Journal of Geophysical Research: Solid Earth*, vol. 125, no. 1, 2020. View at: Publisher Site | Google Scholar - L. Wang, Y. Li, S. Shen et al., “Mechanical behaviours of gas-hydrate-bearing clayey sediments of the South China Sea,”
*Environmental Geotechnics*, vol. 6, no. 1, pp. 1–13, 2019. View at: Publisher Site | Google Scholar - B. N. Madhusudhan, C. R. I. Clayton, and J. A. Priest, “The effects of hydrate on the strength and stiffness of some sands,”
*Journal of Geophysical Research: Solid Earth*, vol. 124, no. 1, pp. 65–75, 2019. View at: Publisher Site | Google Scholar - T. S. Yun, J. C. Santamarina, and C. Ruppel, “Mechanical properties of sand, silt, and clay containing tetrahydrofuran hydrate,”
*Journal of Geophysical Research: Solid Earth*, vol. 112, no. B4, 2007. View at: Publisher Site | Google Scholar - L. Lei and Y. Seol, “Pore-scale investigation of methane hydrate-bearing sediments under triaxial condition,”
*Geophysical Research Letters*, vol. 47, no. 5, 2020. View at: Publisher Site | Google Scholar - L. Wang, X. Sun, S. Shen et al., “Undrained triaxial tests on water-saturated methane hydrate–bearing clayey-silty sediments of the South China Sea,”
*Canadian Geotechnical Journal*, vol. 58, no. 3, pp. 351–366, 2021. View at: Publisher Site | Google Scholar - W. J. Likos, “Effective stress in unsaturated soil: accounting for surface tension and interfacial area,”
*Vadose Zone Journal*, vol. 13, no. 5, 2014. View at: Publisher Site | Google Scholar - Y. Gao, Z. Li, D. Sun, and H. Yu, “A simple method for predicting the hydraulic properties of unsaturated soils with different void ratios,”
*Soil and Tillage Research*, vol. 209, article 104913, 2021. View at: Publisher Site | Google Scholar - Y. Song, T. Luo, B. N. Madhusudhan et al., “Strength behaviors of CH
_{4}hydrate-bearing silty sediments during thermal decomposition,”*Journal of Natural Gas Science and Engineering*, vol. 72, article 103031, 2019. View at: Publisher Site | Google Scholar - Y. Song, F. Yu, Y. Li, W. Liu, and J. Zhao, “Mechanical property of artificial methane hydrate under triaxial compression,”
*Journal of Natural Gas Chemistry*, vol. 19, no. 3, pp. 246–250, 2010. View at: Publisher Site | Google Scholar - M. A. Biot, “Theory of elasticity and consolidation for a porous anisotropic solid,”
*Journal of Applied Physics*, vol. 26, no. 2, pp. 182–185, 1955. View at: Publisher Site | Google Scholar - D. G. Fredlund and N. R. Morgenstern, “Stress state variables for unsaturated soils,”
*Journal of the Geotechnical Engineering Division*, vol. 103, no. 5, pp. 447–466, 1977. View at: Publisher Site | Google Scholar - L. Shao, G. Zheng, X. Guo, and G. Liu, “Principle of effective stress for unsaturated soils,” in
*The Sixth International Conference on Unsaturated Soils*, N. Khalili, A. R. Russell, and A. Khoshghalb, Eds., pp. 239–245, CRC Press, Sydney, 2014. View at: Publisher Site | Google Scholar - C. Zhao, Z. Liu, P. Shi, J. Li, G. Cai, and C. Wei, “Average soil skeleton stress for unsaturated soils and discussion on effective stress,”
*International Journal of Geomechanics*, vol. 16, no. 6, 2016. View at: Publisher Site | Google Scholar - R. de Boer,
*Theory of Porous Media*, Springer, New York, NY, USA, 2000. View at: Publisher Site - S. J. Wheeler, R. S. Sharma, and M. S. R. Buisson, “Coupling of hydraulic hysteresis and stress–strain behaviour in unsaturated soils,”
*Géotechnique*, vol. 53, no. 1, pp. 41–54, 2003. View at: Publisher Site | Google Scholar - L. Shao, X. Guo, S. Liu, and G. Zheng,
*Effective Stress and Equilibrium Equation for Soil Mechanics*, CRC Press, Leiden, 2018. - C. Jommi, “Remarks on the constitutive modelling of unsaturated soils,” in
*Experimental Evidence and Theoretical Approaches in Unsaturated Soils*, A. Tarantio and C. Mancuso, Eds., pp. 139–153, CRC Press, Balkema, Rotterdam, 2000. View at: Publisher Site | Google Scholar - R. F. Craig,
*Craig's Soil Mechanics*, Taylor & Francis Group, London, 7th edition, 2004. View at: Publisher Site - T. W. Lambe and R. V. Whitman,
*Soil Mechanics*, John Wiley & Sons, Inc., New York, NY, USA, 1969.

#### Copyright

Copyright © 2021 Lifeng Zeng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.