# Analytical Solutions for Gas-Water Two-Phase Flow in Multiseam Coalbed Methane Production

Authors: Xiaoji Shang; Zhizhen Zhang; Yixin Niu; Xiao Yang; Feng Gao
Journal: Geofluids
 (2021)

---

## Abstract

Multiseam coalbed methane (CBM) exploitation can not only reduce single-well investment but also increase the length of service of the well and significantly enhance the CBM economic recovery of the entire basin. To compare with and further to guide the actual project of CBM production, this study proposed a conceptual gas-water two-phase separate flow model for single coal seam considering the solubility of gas. This mathematical model was solved analytically by separation of variables and verified through history matching of the production data from the No. 3 seam of 1# test well of Jincheng and then applied to investigate the effect of gas solubility on the gas pressure. Furthermore, based on the coupled two-phase separate flow model of single seam, another two-phase separate flow model for the development of multicoal seam development was established. Similarly, the analytical solution of this model for multicoal seam layers was matched with the in situ data of TS-1 well of Liupanshui coal mine. It is found that the height difference and pressure difference between the two seams play key roles in the multiseam CBM development comprehensively.


---

## Body

## 1. Introduction

To improve vertical well coalbed methane (CBM) production, a great amount of well drainage control methods have been researched, such as vertical drainage of single seam layer, multiseam fracturing drainage of multicoal seam layers, and single-seam fracturing drainage of multicoal seam layers [1]. Through investigation of the multiseam CBM exploitation in Powder River Basin, the U.S. Department of Energy [2] has found that multiseam CBM exploitation can not only reduce single-well investment but also increase the length of service of the well and significantly enhance the CBM economic recovery of the entire basin [3]. Moreover, it is found that some wells adopting multiseam CBM exploitation have been working for more than 10 years in the earlier well group in the Panzhuang block of the southern Qinshui Basin, with the gas production remaining around 1000 m3/d [4]. Thus, the effective development of multicoal seam layers will be helpful to obtain high yield in vertical CBM well production [4, 5].Based on multicoal seam development, a variety of factors which can affect mining such as burial depth, coal thickness, gas content, porosity, permeability, structural setting, and damage mechanism of upper exposed producing layers have been studied [6–10]. Besides, in the technical perspective, several kinds of technology such as well test analysis, numerical simulation, and physical simulation have been applied to the development of multicoal seam [1, 11]. For instance, after the comparison between the theory and the field data from coal reservoir in southern Qinshui Basin, China, Li and Li [1] put forward that the main controlling factors influencing the multicoal seam development consist of the height of the critical production gas, the grades of reservoir pressure, the difference of ability of supply for liquid, and the permeability; in the aspect of numerical simulation, the laws of the geological mining factors to the upper and lower pillar’s stability in multicoal seam strip mining were studied by FLAC3D and the function relation between the stress increasing coefficient of upper and lower pillars with the mining depth, mining widths, interlayer spacing, mining thickness, property of interstratified rock, and the spatial relationship was obtained [11].In order to closely simulate and further to guide the practical production of CBM wells, the modeling of multicoal seam development has become the research emphasis naturally. To select favorable regions for multiseam CBM exploitation, Zhang et al. [3] have conducted a multiobjective fuzzy matter-element model to evaluate the potential of multiseam CBM exploitation in the study area based on the correlation analysis. On the other hand, in general, pore fractures of low-rank coal reservoirs are filled with water, especially in a shallow burial depth of 1000 m [12]; thus, the study of CBM exploitation process is actually the process of two-phase flow, on which many scholars have already conducted large amount of researches [13–15]. Ibrahim and Nasr-El-Din [14] have developed an analytical model for saturated CBM reservoirs where two-phase flow occurs from the start of production, and the model can be used as a prediction or a history matching tool for CBM performance. Clarkson and McGovern [16] developed a model for simulating single-phase flow of water above desorption pressure and two-phase flow of gas and water below desorption pressure, but effects of transient flow were ignored. As a supplement, using a revised DDA model, the model of transition from single-phase to multiphase flow during transient flow for CBM reservoirs has been conducted [13]. Due to the difficulty of obtaining the analytical solution, numerical simulation is always adopted to solve the model founded in the study of two-phase flow; for example, Ziarani et al. [17] have investigated the effect of sorption time on CBM recovery through numerical simulation. However, there is almost no mathematical model focused on the relationship of the fluid pressure, production rate and mining time, mining distance on the development of multiseam layers at two-phase flow stage, not to mention the analytical solution of the fluid pressure and velocity, and so on.It is generally accepted that CBM mainly exists in coal seam in adsorbed state, but more and more studies reveal that dissolved gas should not be ignored [18, 19]. Based on this conclusion, lots of experiments for methane solubility are conducted [12, 20]. Samples taken from Hailar Basin of Xinjiang, China, were physically simulated of dissolved methane content [12]. The experimental results indicate that the solubility of methane is influenced by temperature, pressure, coal reservoir water salinity, etc. But few coupling seepage models considering the solubility factor are established and discussed in theoretical research actually.This paper proposes a two-phase separate flow model for the mining of multicoal seam layers. Firstly, considering the solubility of gas, a two-phase separate flow model of single seam is proposed for the multiphysical process. Secondly, this coupled mathematical model is analytically solved by separation of variables and applied to compare with the field data of the No. 3 seam of 1# test well of Jincheng, Shanxi, for the verification of model accuracy. Thirdly, the impacts of gas solubility on the gas pressure are analyzed. Furthermore, based on the coupled two-phase separate flow model of single seam, another coupled mathematical model of two-phase separate flow for the development of multicoal seam layers is established. Similarly, the analytical solutions of this model for multicoal seam layers is matched with the field data of TS-1 well of Liupanshui coal mine, China [10]. Finally, the impact of interlayer difference on the development of multicoal seam layers is demonstrated.


## 2. Gas-Water Two-Phase Separate Flow Model for Single-Seam CBM Development

Before the establishment of the two-phase flow model, assume that the gas and water two-phase flows are two completely separate fluids with different velocities and properties for the convenience of obtaining the analytical solutions. However, the gas solubility in water is considered in the following models.


### 2.1. Gas-Water Flow Continuity Equations

In the following fluid continuity equations at the two-phase flow stage, subscripts g and w are used to represent the gas phase and water phase, respectively.The continuity equation for gas is(1)∂∂tϕRswswρw+sgρg−∇⋅ρgVg+RswρwVw=ρgqg.The continuity equation for water is(2)∂∂tϕρwsw−∇⋅ρwVw=ρwqw,where ϕ is the coal porosity and Rsw is the gas solubility in water. ρw is the density of water and ρg is the density of gas under formation conditions. Vg is the velocity of gas and Vw is the velocity of water. qg and qw denote the fluid strength of gas and water, respectively. ∇ is the Hamiltonian operator. sg and sw are the saturations of gas and water, respectively, and the pores are fully saturated by gas and water:
(3)sg+sw=1.


### 2.2. Gas-Water Flow Motion Equations

The gas-water coupled seepage equations of motion are [21]
(4)ϕswVw−Vs=−kkrwμw∇pw−ρwg∇H,(5)ϕsgVg−Vs=−kkrgμg∇pg−ρgg∇H,where Vw and Vg are the absolute velocities of water and gas, respectively. Vs is the velocity of the coal solid. k is the absolute permeability of porous media. krw and krg are the relative permeabilities of water and gas, respectively. μw and μg are the viscosities of water and gas, respectively. pw and pg are the pore pressures of water and gas within porous media, respectively. g is the gravity acceleration, and H refers to the vertical elevation.According to Equations (1)–(4), ignoring the second-order phase, the governing equations of the gas-water two-phase flow are obtained.For the water phase,(6)kkrwμwϕsw∇2pw=divVs+∇⋅ρwg∇H+qw−∂∂tϕsw.For the gas phase,(7)1kdivVs+∂∂tϕRswswρw+sgρg−∇⋅krgρgμg+krwRswρwμw∇pg=qgρg,where
(8)divVs=∂θ∂t.θ is the volumetric strain of the coal solid.


### 2.3. Deformation Field Equations of Coalbed

According to the theory of elasticity, the coalbed maintains the force balance during CBM mining, and the differential equation of equilibrium is(9)σij,j+fi=0,where σij is the total stress tensor and fi is the body force tensor.Through the principle of effective stress, the relationship between the total stress tensorσij and the effective stress tensor σij′ can be written as
(10)σij=σij′−pδij,where σij′is the effective stress tensor. p denotes the pore pressure of coalbed. δij is the Kronecker function, and it is given by
(11)δij=1,i=j,0,i≠j.As described above, the coalbed experiences elastic deformation in the CBM mining process; thus, the generalized Hooke law can be used to characterize the stress-strain relationship of the coal.(12)σij′=2Gεij+λθδij+1−αpδij,where εij is the strain tensor, α is the Biot coefficient. G is the shear modulus of elasticity, and λ is the lame constant, and they can be written as
(13)G=E21+v,λ=Ev1+v1−2v,where E is the modulus of elasticity and v is the Poisson ratio.Suppose the coal solid produced only small deformation in the CBM mining process, the following equation can be obtained through the continuity in the deformation process as(14)εij=12ui,j+uj,i,where ui and uj are the displacement components of coal solid.From Equation (14), we can obtain
(15)θ=∇⋅u,where u is the displacement vector.Combined with Equations (10)–(19), the governing equations are simplified into
(16)2G+λ∇θ+f=α∇p.


### 2.4. Fluid Pressure Model at Two-Phase Flow Stage

Combined with Equations (7), (6), and (16), the fluid pressure model in gas-water two-phase flow becomes
(17)Ikkrgρgμg+kkrwRswρwμw∇2pg=∂∂tϕRswswρw+sgρg+αλ+2Gpg+divVs+qgρg,kkrwμwϕsw∇2pw=divVs+∇⋅ρwg∇H+qw−∂∂tϕsw,2G+λ∇θ+f=α∇pThis model can be further simplified into(18)IIkkrgρgμg+kkrwRswρwμw∇2pg=∂∂tϕRswswρw+sgρg+αλ+2G∂∂tpg+qgρg,kkrwμwϕsw∇2pw=αλ+2G∂∂tpw+qw


### 2.5. Analytical Solutions

Before the analytical process, the relationship between the gas density and gas viscosity is derived.In the theoretical analysis of coal mining, the volume coefficientBα is usually introduced due to the fact that the density of ground standard state is easier to measure, (ρα=ραsc/Bα,α=g,w, ραsc is the ground standard state density).Firstly, the compressibility coefficient is introduced. Under isothermal conditions, the compressibility coefficient of the fluidcf is expressed as
(19)cf=−1VdVdp=1ρdρdp,T=constant,whereV is the volume of fluid and p is the pressure. Since the dimension of pressure is Pa or MPa, the unit of compressibility is Pa−1 or MPa−1.Integrating both sides of Equation (21) yields
(20)ρ=ρscecfp−po,where the density ρsc corresponds to the reference pressure pscunder standard atmospheric pressure.The expression between viscosity and pressure is as follows:(21)μp=μpscebp−p0,where μp and μpsc are the viscosities under pressure p and atmospheric pressure psc, respectively. b is the pressure coefficient, which is proportional to void volume and inversely proportional to absolute temperature.In this paper, for the convenience of distinction, the viscosities ofμα and μαsc are corresponding to the pressures of p and psc, respectively, with α=g,w.Combining Equations (20) and (21) yields
(22)ραμα=ραscμαsc⋅ecfp−pscebp−psc.Taylor series expansion is carried out tofx=ex at x=0, which yields
(23)ex=1+11!x+12!x2+⋯+1n!xn+⋯,x∈−∞,+∞and
(24)iρα=ραsc⋅ecfp−p0=ραsc⋅1+cfp−pα0+Rnx→ραsc⋅1+cfp−pα0(25)iiραμα=ραscμαsc⋅ecfp−pα0ebp−pα0=ραscμαsc⋅ecf−bp−pα0=ραscμαsc⋅1+cf−bp−pα0+Rnx→ραscμαsc⋅1,where Rnx=eξx/n+1!xn+10≤ξ≤1, which is the Lagrange residual term of ex, withlimn→∞Rnx=0.


#### 2.5.1. Solving Process

Figure1 shows the schematic diagram of a single-coalbed in CBM mining process. In the single-coalbed CBM development case, both the values of qw and qg from model (II) are 0. Suppose the coal seam is infinite cylinder, the solution of model can be converted into the following definite problems (III) and (IV).
(26)III∂pg∂t=η12∂2pg∂r2+1r∂pg∂r0<r<R,t>0,pgr,0=p0,pgR,t=p0,where
(27)η12=kkrwρw/μw+Rswkkrgρg/μgϕRswswρw/Bw+sgρg/Bg+α/λ+2G.p0 is the initial pore pressure in the coal seam, r is the distance to the center of the well, and t is the time.
(28)IV∂pw∂t=η22∂2pw∂r2+1r∂pw∂r0<r<R,t>0,pwr,0=p00≤r≤R,pwR,t=p0R→∞,where η22=kkrw/μwϕswλ+2G/α.Figure 1
Schematic of development of single coal seam.To simplify the calculation of definite problems (III) and (IV), a general form of definite problem (V) is constructed as(29)V∂p∂t=η2∂2p∂r2+1r∂p∂r0<r<R,t>0,pr,0=p0,pR,t=p0Method of variable separation is used in solving problem (V) as follows.Firstly, assuming that(30)pr,t=V1rT1t.Then, separation of variables is made after this expression generated into problem (V); after that, order(31)rV″+V′rV=T′η2T=−γ,(32)T′+γη2T=0,(33)V″+1rV′+γV=0,Equation (33) represents the zero-order Bessel equation. Then, eigenvalue problems will be solved as
(34)V″+1rV′+γV=0,VR=0R→∞,V0<+∞.It should be noted thatγ≤0 is not an intrinsic value. Multiply Equation (34) with rV, and we have
(35)rV″V+VV′+γrV2=0,namely,
(36)γrV2=−VrV′′.Integrate to both ends of Equation (36) from 0 to R, and we have
(37)γ∫0RrV2dr=∫0RrdVdr2dr.Ifγ≤0, then
(38)∫0RrdVdr2dr≤0,that is to say,
(39)dVdr=0.Hence(40)Vr=C,where C is constant.Nevertheless, as(41)VR=0,there must be
(42)Vr≡0.γ≤0 is not an intrinsic value.Making(43)ρ=γr.ConsideringVρ/γ as Vρ, let us discuss the condition of γ>0 in Equation (34), by
(44)Vr=Vρρr=γVρ,Vrr=γVρρ,1rVr=γργVρ=γρVρ.We can obtain(45)V″+1ρV′+V=0.The solution of Equation (45) is
(46)Vr=c1J0ρ+c2N0ρ=c1J0γr+c2N0γr,where J0 represents zero-order Bessel function and N0represents zero-order Neumann functions, which have the approximate expression, respectively.
(47)J0x=1−14x2,N0x=2πlnx2+c,where c is Euler’s constant, with a value of 0.577216. According to the natural boundary conditions, the value of c2 is 0. And the boundary condition is
(48)c1J0γR=0.Ifc1≠0, then
(49)J0γR=0.Yet, the number of zero point of zero-order Bessel function can be several. Set the zero point toμm0m=1,2,⋯; thus,
(50)γmR=μm0,γm=μm0R2m=1,2,⋯.Intrinsic function is(51)Vm=cmJ0μm0Rrm=1,2,⋯.Substituting the value ofγm into Equation (32), we have
(52)T′+μm0ηR2T=0,Tm=Dme−μm0η/R2t,pmr,t=Ame−μm0η/R2tJ0μm0Rr.To overlaypm, the constant Am is determined as follows:
(53)pr,t=∑m=1∞pmr,t=∑m=1∞Ame−μm0η/R2tJ0μm0Rr=∑m=1∞Ame−μm0η/R2tJ0μm0Rr,where
(54)pr,0=∑m=1∞AmJ0μm0Rr=p0.Owing to the orthogonal system with weightr constituted by J0μm0/Rrm=1,2,⋯,
(55)∫0RrJ0μm0RrJ0μm0Rrdr=0m≠n,Mm2m=n.Am is the generalized Fourier coefficient of pr,0=p0, namely,
(56)Am=1Mm2∫0Rrpr,0J0μm0Rrdrm=1,2,⋯,SubstitutingAm into Equation (53), the solution of the definite problem (V) is
(57)pr,t=∑m=1∞1Mm2e−μm0η/R2tJ0μm0Rr∫0Rpr,0J0μm0Rrdr,where
(58)pr,0=∑m=1∞1Mm2J0μm0Rr∫0Rrpr,0J0μm0Rrdr.Substitute Equations (51) and (62) into Equation (57); hence,
(59)Mm2=∫0RrJ0μm0RrJ0μm0Rrdr=∫0Rr1−14μm0Rr22dr.Givenm=2, by the definitionμm0=±2, the above equation becomes
(60)Mm2=∫0Rr1−14μm0Rr22dr=R26.Substituting Equation (60) into Equation (57), the analytical solutions of definite problem (V) is
(61)pr,t=p0e−2η/R2t.


#### 2.5.2. Final Solution

Combining Equation (61) with definite problems (III) and (IV), the solution of gas pressure at the two-phase flow stage is
(62)pgr,t=p0e−2η1/R2t,where
(63)η12=KKrwρw/μw+RswKKrgρg/μgϕRswswρw/Bw+sgρg/Bg+α/λ+2G.The solution of water pressure at the two-phase flow stage is(64)pwr,t=p0e−2η2/R2t,where η22=kkrw/μwϕsw⋅λ+2G/α.


### 2.6. Model Verification

#### 2.6.1. Field Data

The engineering data is from Jincheng 1# test well of Shanxi province, China. The No. 3 seam experiences single coal seam development, whose mineral depth is from 521.6 m to 527.4 m, the minimum horizontal stress on the main seam tested by hydraulic fracturing is 7.9 MPa~9.4 MPa, and the maximum horizontal stress determined from long spacing sonic and coal reservoir fracture pressure is 13.7 MPa~15.7 MPa. Meanwhile, the actual measured gas content is 19.6 m3/t, and the free gas is in dynamic equilibrium with the adsorbed gas, with gas saturation 73.0%. The single-well drainage area is within S=0.09 km2 .Table1 describes the daily actual measured gas production and water production of different time and exploitation stages, among of which I-V represents water flow stage, free gas desorption stage, fluctuation desorption stage, stability desorption stage, and attenuation desorption stage, respectively.Table 1
Production of methane and water in No. 3 seam of 1# test well, Shanxi.
StageMining time (t, d)Mean daily water production (v, m3/d)Cumulative daily water production (q, m3)Mean daily gas production (v, m3/d)Cumulative daily gas production (q, m3)(I) Water flow stage1203.95482.100(II) Free gas desorption444.8211.43398.4146011(III) Fluctuation desorption1292.96381.83844.2495907(IV) Stable desorption15000.575025003750000(V) Decreasing desorption12000.2240500600000Table2 shows physical parameters of coal bed; various physical parameters will be validated with the theoretical model.Table 2
Physical parameters of No. 3 seam of 1# test well, Shanxi.
ParameterValueYoung’s modulus of coal (E, GPa)2.735Poisson’s ratio of coalυ0.25Coal density (ρ, kg/m3)1500Seam depth (h, m)521.6Seam thickness (H, m)5.8Porosityϕ4.29%Permeability (K, 10-3 μm 2)3.801Seam initial pressure (p, MPa)3.99Water viscosity in the ground condition (μw, Pa·s)5.8×10−4Gas viscosity under standard atmospheric pressure (μgsc, Pa·s)1.47×10−5Gas density under standard atmospheric pressure (ρgsc, kg/m3)0.425Water density under standard condition (ρw, kg/m3)998Formation volume factor of gasBg0.00896Formation volume factor of waterBw1.0The initial gas saturationsg0.73Water saturationsw0.27SolubilityRsw1.3Relative permeability of gasKrg1Relative permeability of waterKrw0.0001Seam radius (R, m)300


#### 2.6.2. Verification of Theoretical Model on Fluid Pressure

As the analytical solution of the gas and water pressure has been given by Equations (62) and (64), substituting each parameter value into the theoretical model yields the gas pressure near the wellhead:
(65)pwr,t=2.409e−3.7113×10−3t,pgr,t=2.69e−5.2857×10−4t,t>120.As shown in Figure2, the gas pressure curve near the wellbore calculated by the theoretical model was compared with the field pressure data from No. 3 seam of 1# test well. The figure shows that water will be drained from the beginning in the mining process, while gas production occurred after the first 120 days. Firstly, both the pressure of gas and water are decreasing; secondly, the speed of the fluid decrease eventually tends to 0, and the pressure change approaches flattening. This is because as gas desorption occurs with time increasing, correspondingly, the increasing gas migration velocity results in the gas pressure drop. With the effect of desorption, the amount of desorbed gas gradually reduces, thereby leading to the drop of changing rate of gas pressure.Figure 2
Variation of fluid pressure with time.It is given that the original reservoir pressure is 3.99 MPa, gas pressure is 2.69 MPa, critical desorption pressure is 2.41 MPa, and the depletion pressure 0.3 MPa. The pressure obtained from this two-phase separate flow model can well fit the measured field data.


#### 2.6.3. Verification of Theoretical Model on Gas Flow Rate

The shale gas production rate is defined as [22]
(66)dGptdt=−∫ϕpadpgdtdv,where pa=101.325kPa.Based on Darcy’s law, the relationship between shale gas production and water production is obtained by Xue et al. [23].
(67)GpgGpw=krgkrwΔpgΔpwμwμg.Equation (67) can be transformed into the final form of water production rate as
(68)dGpwtdt=krwkrg⋅μgμw⋅ΔpwΔpg⋅dGptdt.Substituting each parameter of Table2 into Equations (66) and (68) can yield the results of the water and gas production rate. Figures 3 and 4 represent the performance of the water production rate and gas production rate with time, respectively. The figure shows that in the interval of 0~120 days, the coal seam is at the single-phase fluid flow stage from the in situ data, and no gas is output at this stage; in the interval of 120~295 days, the fluid flow is at a nonsaturated stage, and since the speed at this stage does not seek a mathematical expression, the curve is expressed by the boundary value of the single-phase fluid flow stage and two-phase flow stage; furthermore, the seam enters into a stable desorption stage after 295 days. In general, the overall change of the theoretical curve and measured value is of the same trend, such as the gas production and average water production reduce with time going on, and the rate of change slow down. It is necessary to note that the theoretical curves show the trend of the ideal state, while some factors in engineering will influence the trend of gas production rate, such as the amount of water pumped in hydraulic fracturing will affect the stability of daily gas production.Figure 3
Variation of water production rate with time.Figure 4
Variation of gas production rate with time.


### 2.7. The Influence of Gas Solubility on Changes of Gas Pressure

For a more comprehensive consideration, the solubility factor is intentionally added in the model for two-phase flow. Figure5 describes the effect of the solubility on gas pressure. As is shown in Figure 5, as the value of gas solubility in water becomes larger, the decreasing magnitude of gas pressure increases, and the scope of mining becomes wider. This is due to the fact that the greater the solubility of gas in water in depressurization process, the easier the desorption of gas, and the sooner the permeation rate of gas, thus the greater the pressure drop.Figure 5
Effect of gas solubility on pressure change.


## 2.1. Gas-Water Flow Continuity Equations

In the following fluid continuity equations at the two-phase flow stage, subscripts g and w are used to represent the gas phase and water phase, respectively.The continuity equation for gas is(1)∂∂tϕRswswρw+sgρg−∇⋅ρgVg+RswρwVw=ρgqg.The continuity equation for water is(2)∂∂tϕρwsw−∇⋅ρwVw=ρwqw,where ϕ is the coal porosity and Rsw is the gas solubility in water. ρw is the density of water and ρg is the density of gas under formation conditions. Vg is the velocity of gas and Vw is the velocity of water. qg and qw denote the fluid strength of gas and water, respectively. ∇ is the Hamiltonian operator. sg and sw are the saturations of gas and water, respectively, and the pores are fully saturated by gas and water:
(3)sg+sw=1.


## 2.2. Gas-Water Flow Motion Equations

The gas-water coupled seepage equations of motion are [21]
(4)ϕswVw−Vs=−kkrwμw∇pw−ρwg∇H,(5)ϕsgVg−Vs=−kkrgμg∇pg−ρgg∇H,where Vw and Vg are the absolute velocities of water and gas, respectively. Vs is the velocity of the coal solid. k is the absolute permeability of porous media. krw and krg are the relative permeabilities of water and gas, respectively. μw and μg are the viscosities of water and gas, respectively. pw and pg are the pore pressures of water and gas within porous media, respectively. g is the gravity acceleration, and H refers to the vertical elevation.According to Equations (1)–(4), ignoring the second-order phase, the governing equations of the gas-water two-phase flow are obtained.For the water phase,(6)kkrwμwϕsw∇2pw=divVs+∇⋅ρwg∇H+qw−∂∂tϕsw.For the gas phase,(7)1kdivVs+∂∂tϕRswswρw+sgρg−∇⋅krgρgμg+krwRswρwμw∇pg=qgρg,where
(8)divVs=∂θ∂t.θ is the volumetric strain of the coal solid.


## 2.3. Deformation Field Equations of Coalbed

According to the theory of elasticity, the coalbed maintains the force balance during CBM mining, and the differential equation of equilibrium is(9)σij,j+fi=0,where σij is the total stress tensor and fi is the body force tensor.Through the principle of effective stress, the relationship between the total stress tensorσij and the effective stress tensor σij′ can be written as
(10)σij=σij′−pδij,where σij′is the effective stress tensor. p denotes the pore pressure of coalbed. δij is the Kronecker function, and it is given by
(11)δij=1,i=j,0,i≠j.As described above, the coalbed experiences elastic deformation in the CBM mining process; thus, the generalized Hooke law can be used to characterize the stress-strain relationship of the coal.(12)σij′=2Gεij+λθδij+1−αpδij,where εij is the strain tensor, α is the Biot coefficient. G is the shear modulus of elasticity, and λ is the lame constant, and they can be written as
(13)G=E21+v,λ=Ev1+v1−2v,where E is the modulus of elasticity and v is the Poisson ratio.Suppose the coal solid produced only small deformation in the CBM mining process, the following equation can be obtained through the continuity in the deformation process as(14)εij=12ui,j+uj,i,where ui and uj are the displacement components of coal solid.From Equation (14), we can obtain
(15)θ=∇⋅u,where u is the displacement vector.Combined with Equations (10)–(19), the governing equations are simplified into
(16)2G+λ∇θ+f=α∇p.


## 2.4. Fluid Pressure Model at Two-Phase Flow Stage

Combined with Equations (7), (6), and (16), the fluid pressure model in gas-water two-phase flow becomes
(17)Ikkrgρgμg+kkrwRswρwμw∇2pg=∂∂tϕRswswρw+sgρg+αλ+2Gpg+divVs+qgρg,kkrwμwϕsw∇2pw=divVs+∇⋅ρwg∇H+qw−∂∂tϕsw,2G+λ∇θ+f=α∇pThis model can be further simplified into(18)IIkkrgρgμg+kkrwRswρwμw∇2pg=∂∂tϕRswswρw+sgρg+αλ+2G∂∂tpg+qgρg,kkrwμwϕsw∇2pw=αλ+2G∂∂tpw+qw


## 2.5. Analytical Solutions

Before the analytical process, the relationship between the gas density and gas viscosity is derived.In the theoretical analysis of coal mining, the volume coefficientBα is usually introduced due to the fact that the density of ground standard state is easier to measure, (ρα=ραsc/Bα,α=g,w, ραsc is the ground standard state density).Firstly, the compressibility coefficient is introduced. Under isothermal conditions, the compressibility coefficient of the fluidcf is expressed as
(19)cf=−1VdVdp=1ρdρdp,T=constant,whereV is the volume of fluid and p is the pressure. Since the dimension of pressure is Pa or MPa, the unit of compressibility is Pa−1 or MPa−1.Integrating both sides of Equation (21) yields
(20)ρ=ρscecfp−po,where the density ρsc corresponds to the reference pressure pscunder standard atmospheric pressure.The expression between viscosity and pressure is as follows:(21)μp=μpscebp−p0,where μp and μpsc are the viscosities under pressure p and atmospheric pressure psc, respectively. b is the pressure coefficient, which is proportional to void volume and inversely proportional to absolute temperature.In this paper, for the convenience of distinction, the viscosities ofμα and μαsc are corresponding to the pressures of p and psc, respectively, with α=g,w.Combining Equations (20) and (21) yields
(22)ραμα=ραscμαsc⋅ecfp−pscebp−psc.Taylor series expansion is carried out tofx=ex at x=0, which yields
(23)ex=1+11!x+12!x2+⋯+1n!xn+⋯,x∈−∞,+∞and
(24)iρα=ραsc⋅ecfp−p0=ραsc⋅1+cfp−pα0+Rnx→ραsc⋅1+cfp−pα0(25)iiραμα=ραscμαsc⋅ecfp−pα0ebp−pα0=ραscμαsc⋅ecf−bp−pα0=ραscμαsc⋅1+cf−bp−pα0+Rnx→ραscμαsc⋅1,where Rnx=eξx/n+1!xn+10≤ξ≤1, which is the Lagrange residual term of ex, withlimn→∞Rnx=0.


### 2.5.1. Solving Process

Figure1 shows the schematic diagram of a single-coalbed in CBM mining process. In the single-coalbed CBM development case, both the values of qw and qg from model (II) are 0. Suppose the coal seam is infinite cylinder, the solution of model can be converted into the following definite problems (III) and (IV).
(26)III∂pg∂t=η12∂2pg∂r2+1r∂pg∂r0<r<R,t>0,pgr,0=p0,pgR,t=p0,where
(27)η12=kkrwρw/μw+Rswkkrgρg/μgϕRswswρw/Bw+sgρg/Bg+α/λ+2G.p0 is the initial pore pressure in the coal seam, r is the distance to the center of the well, and t is the time.
(28)IV∂pw∂t=η22∂2pw∂r2+1r∂pw∂r0<r<R,t>0,pwr,0=p00≤r≤R,pwR,t=p0R→∞,where η22=kkrw/μwϕswλ+2G/α.Figure 1
Schematic of development of single coal seam.To simplify the calculation of definite problems (III) and (IV), a general form of definite problem (V) is constructed as(29)V∂p∂t=η2∂2p∂r2+1r∂p∂r0<r<R,t>0,pr,0=p0,pR,t=p0Method of variable separation is used in solving problem (V) as follows.Firstly, assuming that(30)pr,t=V1rT1t.Then, separation of variables is made after this expression generated into problem (V); after that, order(31)rV″+V′rV=T′η2T=−γ,(32)T′+γη2T=0,(33)V″+1rV′+γV=0,Equation (33) represents the zero-order Bessel equation. Then, eigenvalue problems will be solved as
(34)V″+1rV′+γV=0,VR=0R→∞,V0<+∞.It should be noted thatγ≤0 is not an intrinsic value. Multiply Equation (34) with rV, and we have
(35)rV″V+VV′+γrV2=0,namely,
(36)γrV2=−VrV′′.Integrate to both ends of Equation (36) from 0 to R, and we have
(37)γ∫0RrV2dr=∫0RrdVdr2dr.Ifγ≤0, then
(38)∫0RrdVdr2dr≤0,that is to say,
(39)dVdr=0.Hence(40)Vr=C,where C is constant.Nevertheless, as(41)VR=0,there must be
(42)Vr≡0.γ≤0 is not an intrinsic value.Making(43)ρ=γr.ConsideringVρ/γ as Vρ, let us discuss the condition of γ>0 in Equation (34), by
(44)Vr=Vρρr=γVρ,Vrr=γVρρ,1rVr=γργVρ=γρVρ.We can obtain(45)V″+1ρV′+V=0.The solution of Equation (45) is
(46)Vr=c1J0ρ+c2N0ρ=c1J0γr+c2N0γr,where J0 represents zero-order Bessel function and N0represents zero-order Neumann functions, which have the approximate expression, respectively.
(47)J0x=1−14x2,N0x=2πlnx2+c,where c is Euler’s constant, with a value of 0.577216. According to the natural boundary conditions, the value of c2 is 0. And the boundary condition is
(48)c1J0γR=0.Ifc1≠0, then
(49)J0γR=0.Yet, the number of zero point of zero-order Bessel function can be several. Set the zero point toμm0m=1,2,⋯; thus,
(50)γmR=μm0,γm=μm0R2m=1,2,⋯.Intrinsic function is(51)Vm=cmJ0μm0Rrm=1,2,⋯.Substituting the value ofγm into Equation (32), we have
(52)T′+μm0ηR2T=0,Tm=Dme−μm0η/R2t,pmr,t=Ame−μm0η/R2tJ0μm0Rr.To overlaypm, the constant Am is determined as follows:
(53)pr,t=∑m=1∞pmr,t=∑m=1∞Ame−μm0η/R2tJ0μm0Rr=∑m=1∞Ame−μm0η/R2tJ0μm0Rr,where
(54)pr,0=∑m=1∞AmJ0μm0Rr=p0.Owing to the orthogonal system with weightr constituted by J0μm0/Rrm=1,2,⋯,
(55)∫0RrJ0μm0RrJ0μm0Rrdr=0m≠n,Mm2m=n.Am is the generalized Fourier coefficient of pr,0=p0, namely,
(56)Am=1Mm2∫0Rrpr,0J0μm0Rrdrm=1,2,⋯,SubstitutingAm into Equation (53), the solution of the definite problem (V) is
(57)pr,t=∑m=1∞1Mm2e−μm0η/R2tJ0μm0Rr∫0Rpr,0J0μm0Rrdr,where
(58)pr,0=∑m=1∞1Mm2J0μm0Rr∫0Rrpr,0J0μm0Rrdr.Substitute Equations (51) and (62) into Equation (57); hence,
(59)Mm2=∫0RrJ0μm0RrJ0μm0Rrdr=∫0Rr1−14μm0Rr22dr.Givenm=2, by the definitionμm0=±2, the above equation becomes
(60)Mm2=∫0Rr1−14μm0Rr22dr=R26.Substituting Equation (60) into Equation (57), the analytical solutions of definite problem (V) is
(61)pr,t=p0e−2η/R2t.


### 2.5.2. Final Solution

Combining Equation (61) with definite problems (III) and (IV), the solution of gas pressure at the two-phase flow stage is
(62)pgr,t=p0e−2η1/R2t,where
(63)η12=KKrwρw/μw+RswKKrgρg/μgϕRswswρw/Bw+sgρg/Bg+α/λ+2G.The solution of water pressure at the two-phase flow stage is(64)pwr,t=p0e−2η2/R2t,where η22=kkrw/μwϕsw⋅λ+2G/α.


## 2.5.1. Solving Process

Figure1 shows the schematic diagram of a single-coalbed in CBM mining process. In the single-coalbed CBM development case, both the values of qw and qg from model (II) are 0. Suppose the coal seam is infinite cylinder, the solution of model can be converted into the following definite problems (III) and (IV).
(26)III∂pg∂t=η12∂2pg∂r2+1r∂pg∂r0<r<R,t>0,pgr,0=p0,pgR,t=p0,where
(27)η12=kkrwρw/μw+Rswkkrgρg/μgϕRswswρw/Bw+sgρg/Bg+α/λ+2G.p0 is the initial pore pressure in the coal seam, r is the distance to the center of the well, and t is the time.
(28)IV∂pw∂t=η22∂2pw∂r2+1r∂pw∂r0<r<R,t>0,pwr,0=p00≤r≤R,pwR,t=p0R→∞,where η22=kkrw/μwϕswλ+2G/α.Figure 1
Schematic of development of single coal seam.To simplify the calculation of definite problems (III) and (IV), a general form of definite problem (V) is constructed as(29)V∂p∂t=η2∂2p∂r2+1r∂p∂r0<r<R,t>0,pr,0=p0,pR,t=p0Method of variable separation is used in solving problem (V) as follows.Firstly, assuming that(30)pr,t=V1rT1t.Then, separation of variables is made after this expression generated into problem (V); after that, order(31)rV″+V′rV=T′η2T=−γ,(32)T′+γη2T=0,(33)V″+1rV′+γV=0,Equation (33) represents the zero-order Bessel equation. Then, eigenvalue problems will be solved as
(34)V″+1rV′+γV=0,VR=0R→∞,V0<+∞.It should be noted thatγ≤0 is not an intrinsic value. Multiply Equation (34) with rV, and we have
(35)rV″V+VV′+γrV2=0,namely,
(36)γrV2=−VrV′′.Integrate to both ends of Equation (36) from 0 to R, and we have
(37)γ∫0RrV2dr=∫0RrdVdr2dr.Ifγ≤0, then
(38)∫0RrdVdr2dr≤0,that is to say,
(39)dVdr=0.Hence(40)Vr=C,where C is constant.Nevertheless, as(41)VR=0,there must be
(42)Vr≡0.γ≤0 is not an intrinsic value.Making(43)ρ=γr.ConsideringVρ/γ as Vρ, let us discuss the condition of γ>0 in Equation (34), by
(44)Vr=Vρρr=γVρ,Vrr=γVρρ,1rVr=γργVρ=γρVρ.We can obtain(45)V″+1ρV′+V=0.The solution of Equation (45) is
(46)Vr=c1J0ρ+c2N0ρ=c1J0γr+c2N0γr,where J0 represents zero-order Bessel function and N0represents zero-order Neumann functions, which have the approximate expression, respectively.
(47)J0x=1−14x2,N0x=2πlnx2+c,where c is Euler’s constant, with a value of 0.577216. According to the natural boundary conditions, the value of c2 is 0. And the boundary condition is
(48)c1J0γR=0.Ifc1≠0, then
(49)J0γR=0.Yet, the number of zero point of zero-order Bessel function can be several. Set the zero point toμm0m=1,2,⋯; thus,
(50)γmR=μm0,γm=μm0R2m=1,2,⋯.Intrinsic function is(51)Vm=cmJ0μm0Rrm=1,2,⋯.Substituting the value ofγm into Equation (32), we have
(52)T′+μm0ηR2T=0,Tm=Dme−μm0η/R2t,pmr,t=Ame−μm0η/R2tJ0μm0Rr.To overlaypm, the constant Am is determined as follows:
(53)pr,t=∑m=1∞pmr,t=∑m=1∞Ame−μm0η/R2tJ0μm0Rr=∑m=1∞Ame−μm0η/R2tJ0μm0Rr,where
(54)pr,0=∑m=1∞AmJ0μm0Rr=p0.Owing to the orthogonal system with weightr constituted by J0μm0/Rrm=1,2,⋯,
(55)∫0RrJ0μm0RrJ0μm0Rrdr=0m≠n,Mm2m=n.Am is the generalized Fourier coefficient of pr,0=p0, namely,
(56)Am=1Mm2∫0Rrpr,0J0μm0Rrdrm=1,2,⋯,SubstitutingAm into Equation (53), the solution of the definite problem (V) is
(57)pr,t=∑m=1∞1Mm2e−μm0η/R2tJ0μm0Rr∫0Rpr,0J0μm0Rrdr,where
(58)pr,0=∑m=1∞1Mm2J0μm0Rr∫0Rrpr,0J0μm0Rrdr.Substitute Equations (51) and (62) into Equation (57); hence,
(59)Mm2=∫0RrJ0μm0RrJ0μm0Rrdr=∫0Rr1−14μm0Rr22dr.Givenm=2, by the definitionμm0=±2, the above equation becomes
(60)Mm2=∫0Rr1−14μm0Rr22dr=R26.Substituting Equation (60) into Equation (57), the analytical solutions of definite problem (V) is
(61)pr,t=p0e−2η/R2t.


## 2.5.2. Final Solution

Combining Equation (61) with definite problems (III) and (IV), the solution of gas pressure at the two-phase flow stage is
(62)pgr,t=p0e−2η1/R2t,where
(63)η12=KKrwρw/μw+RswKKrgρg/μgϕRswswρw/Bw+sgρg/Bg+α/λ+2G.The solution of water pressure at the two-phase flow stage is(64)pwr,t=p0e−2η2/R2t,where η22=kkrw/μwϕsw⋅λ+2G/α.


## 2.6. Model Verification

### 2.6.1. Field Data

The engineering data is from Jincheng 1# test well of Shanxi province, China. The No. 3 seam experiences single coal seam development, whose mineral depth is from 521.6 m to 527.4 m, the minimum horizontal stress on the main seam tested by hydraulic fracturing is 7.9 MPa~9.4 MPa, and the maximum horizontal stress determined from long spacing sonic and coal reservoir fracture pressure is 13.7 MPa~15.7 MPa. Meanwhile, the actual measured gas content is 19.6 m3/t, and the free gas is in dynamic equilibrium with the adsorbed gas, with gas saturation 73.0%. The single-well drainage area is within S=0.09 km2 .Table1 describes the daily actual measured gas production and water production of different time and exploitation stages, among of which I-V represents water flow stage, free gas desorption stage, fluctuation desorption stage, stability desorption stage, and attenuation desorption stage, respectively.Table 1
Production of methane and water in No. 3 seam of 1# test well, Shanxi.
StageMining time (t, d)Mean daily water production (v, m3/d)Cumulative daily water production (q, m3)Mean daily gas production (v, m3/d)Cumulative daily gas production (q, m3)(I) Water flow stage1203.95482.100(II) Free gas desorption444.8211.43398.4146011(III) Fluctuation desorption1292.96381.83844.2495907(IV) Stable desorption15000.575025003750000(V) Decreasing desorption12000.2240500600000Table2 shows physical parameters of coal bed; various physical parameters will be validated with the theoretical model.Table 2
Physical parameters of No. 3 seam of 1# test well, Shanxi.
ParameterValueYoung’s modulus of coal (E, GPa)2.735Poisson’s ratio of coalυ0.25Coal density (ρ, kg/m3)1500Seam depth (h, m)521.6Seam thickness (H, m)5.8Porosityϕ4.29%Permeability (K, 10-3 μm 2)3.801Seam initial pressure (p, MPa)3.99Water viscosity in the ground condition (μw, Pa·s)5.8×10−4Gas viscosity under standard atmospheric pressure (μgsc, Pa·s)1.47×10−5Gas density under standard atmospheric pressure (ρgsc, kg/m3)0.425Water density under standard condition (ρw, kg/m3)998Formation volume factor of gasBg0.00896Formation volume factor of waterBw1.0The initial gas saturationsg0.73Water saturationsw0.27SolubilityRsw1.3Relative permeability of gasKrg1Relative permeability of waterKrw0.0001Seam radius (R, m)300


### 2.6.2. Verification of Theoretical Model on Fluid Pressure

As the analytical solution of the gas and water pressure has been given by Equations (62) and (64), substituting each parameter value into the theoretical model yields the gas pressure near the wellhead:
(65)pwr,t=2.409e−3.7113×10−3t,pgr,t=2.69e−5.2857×10−4t,t>120.As shown in Figure2, the gas pressure curve near the wellbore calculated by the theoretical model was compared with the field pressure data from No. 3 seam of 1# test well. The figure shows that water will be drained from the beginning in the mining process, while gas production occurred after the first 120 days. Firstly, both the pressure of gas and water are decreasing; secondly, the speed of the fluid decrease eventually tends to 0, and the pressure change approaches flattening. This is because as gas desorption occurs with time increasing, correspondingly, the increasing gas migration velocity results in the gas pressure drop. With the effect of desorption, the amount of desorbed gas gradually reduces, thereby leading to the drop of changing rate of gas pressure.Figure 2
Variation of fluid pressure with time.It is given that the original reservoir pressure is 3.99 MPa, gas pressure is 2.69 MPa, critical desorption pressure is 2.41 MPa, and the depletion pressure 0.3 MPa. The pressure obtained from this two-phase separate flow model can well fit the measured field data.


### 2.6.3. Verification of Theoretical Model on Gas Flow Rate

The shale gas production rate is defined as [22]
(66)dGptdt=−∫ϕpadpgdtdv,where pa=101.325kPa.Based on Darcy’s law, the relationship between shale gas production and water production is obtained by Xue et al. [23].
(67)GpgGpw=krgkrwΔpgΔpwμwμg.Equation (67) can be transformed into the final form of water production rate as
(68)dGpwtdt=krwkrg⋅μgμw⋅ΔpwΔpg⋅dGptdt.Substituting each parameter of Table2 into Equations (66) and (68) can yield the results of the water and gas production rate. Figures 3 and 4 represent the performance of the water production rate and gas production rate with time, respectively. The figure shows that in the interval of 0~120 days, the coal seam is at the single-phase fluid flow stage from the in situ data, and no gas is output at this stage; in the interval of 120~295 days, the fluid flow is at a nonsaturated stage, and since the speed at this stage does not seek a mathematical expression, the curve is expressed by the boundary value of the single-phase fluid flow stage and two-phase flow stage; furthermore, the seam enters into a stable desorption stage after 295 days. In general, the overall change of the theoretical curve and measured value is of the same trend, such as the gas production and average water production reduce with time going on, and the rate of change slow down. It is necessary to note that the theoretical curves show the trend of the ideal state, while some factors in engineering will influence the trend of gas production rate, such as the amount of water pumped in hydraulic fracturing will affect the stability of daily gas production.Figure 3
Variation of water production rate with time.Figure 4
Variation of gas production rate with time.


## 2.6.1. Field Data

The engineering data is from Jincheng 1# test well of Shanxi province, China. The No. 3 seam experiences single coal seam development, whose mineral depth is from 521.6 m to 527.4 m, the minimum horizontal stress on the main seam tested by hydraulic fracturing is 7.9 MPa~9.4 MPa, and the maximum horizontal stress determined from long spacing sonic and coal reservoir fracture pressure is 13.7 MPa~15.7 MPa. Meanwhile, the actual measured gas content is 19.6 m3/t, and the free gas is in dynamic equilibrium with the adsorbed gas, with gas saturation 73.0%. The single-well drainage area is within S=0.09 km2 .Table1 describes the daily actual measured gas production and water production of different time and exploitation stages, among of which I-V represents water flow stage, free gas desorption stage, fluctuation desorption stage, stability desorption stage, and attenuation desorption stage, respectively.Table 1
Production of methane and water in No. 3 seam of 1# test well, Shanxi.
StageMining time (t, d)Mean daily water production (v, m3/d)Cumulative daily water production (q, m3)Mean daily gas production (v, m3/d)Cumulative daily gas production (q, m3)(I) Water flow stage1203.95482.100(II) Free gas desorption444.8211.43398.4146011(III) Fluctuation desorption1292.96381.83844.2495907(IV) Stable desorption15000.575025003750000(V) Decreasing desorption12000.2240500600000Table2 shows physical parameters of coal bed; various physical parameters will be validated with the theoretical model.Table 2
Physical parameters of No. 3 seam of 1# test well, Shanxi.
ParameterValueYoung’s modulus of coal (E, GPa)2.735Poisson’s ratio of coalυ0.25Coal density (ρ, kg/m3)1500Seam depth (h, m)521.6Seam thickness (H, m)5.8Porosityϕ4.29%Permeability (K, 10-3 μm 2)3.801Seam initial pressure (p, MPa)3.99Water viscosity in the ground condition (μw, Pa·s)5.8×10−4Gas viscosity under standard atmospheric pressure (μgsc, Pa·s)1.47×10−5Gas density under standard atmospheric pressure (ρgsc, kg/m3)0.425Water density under standard condition (ρw, kg/m3)998Formation volume factor of gasBg0.00896Formation volume factor of waterBw1.0The initial gas saturationsg0.73Water saturationsw0.27SolubilityRsw1.3Relative permeability of gasKrg1Relative permeability of waterKrw0.0001Seam radius (R, m)300


## 2.6.2. Verification of Theoretical Model on Fluid Pressure

As the analytical solution of the gas and water pressure has been given by Equations (62) and (64), substituting each parameter value into the theoretical model yields the gas pressure near the wellhead:
(65)pwr,t=2.409e−3.7113×10−3t,pgr,t=2.69e−5.2857×10−4t,t>120.As shown in Figure2, the gas pressure curve near the wellbore calculated by the theoretical model was compared with the field pressure data from No. 3 seam of 1# test well. The figure shows that water will be drained from the beginning in the mining process, while gas production occurred after the first 120 days. Firstly, both the pressure of gas and water are decreasing; secondly, the speed of the fluid decrease eventually tends to 0, and the pressure change approaches flattening. This is because as gas desorption occurs with time increasing, correspondingly, the increasing gas migration velocity results in the gas pressure drop. With the effect of desorption, the amount of desorbed gas gradually reduces, thereby leading to the drop of changing rate of gas pressure.Figure 2
Variation of fluid pressure with time.It is given that the original reservoir pressure is 3.99 MPa, gas pressure is 2.69 MPa, critical desorption pressure is 2.41 MPa, and the depletion pressure 0.3 MPa. The pressure obtained from this two-phase separate flow model can well fit the measured field data.


## 2.6.3. Verification of Theoretical Model on Gas Flow Rate

The shale gas production rate is defined as [22]
(66)dGptdt=−∫ϕpadpgdtdv,where pa=101.325kPa.Based on Darcy’s law, the relationship between shale gas production and water production is obtained by Xue et al. [23].
(67)GpgGpw=krgkrwΔpgΔpwμwμg.Equation (67) can be transformed into the final form of water production rate as
(68)dGpwtdt=krwkrg⋅μgμw⋅ΔpwΔpg⋅dGptdt.Substituting each parameter of Table2 into Equations (66) and (68) can yield the results of the water and gas production rate. Figures 3 and 4 represent the performance of the water production rate and gas production rate with time, respectively. The figure shows that in the interval of 0~120 days, the coal seam is at the single-phase fluid flow stage from the in situ data, and no gas is output at this stage; in the interval of 120~295 days, the fluid flow is at a nonsaturated stage, and since the speed at this stage does not seek a mathematical expression, the curve is expressed by the boundary value of the single-phase fluid flow stage and two-phase flow stage; furthermore, the seam enters into a stable desorption stage after 295 days. In general, the overall change of the theoretical curve and measured value is of the same trend, such as the gas production and average water production reduce with time going on, and the rate of change slow down. It is necessary to note that the theoretical curves show the trend of the ideal state, while some factors in engineering will influence the trend of gas production rate, such as the amount of water pumped in hydraulic fracturing will affect the stability of daily gas production.Figure 3
Variation of water production rate with time.Figure 4
Variation of gas production rate with time.


## 2.7. The Influence of Gas Solubility on Changes of Gas Pressure

For a more comprehensive consideration, the solubility factor is intentionally added in the model for two-phase flow. Figure5 describes the effect of the solubility on gas pressure. As is shown in Figure 5, as the value of gas solubility in water becomes larger, the decreasing magnitude of gas pressure increases, and the scope of mining becomes wider. This is due to the fact that the greater the solubility of gas in water in depressurization process, the easier the desorption of gas, and the sooner the permeation rate of gas, thus the greater the pressure drop.Figure 5
Effect of gas solubility on pressure change.


## 3. Gas-Water Two-Phase Separate Flow Model for Multiseam CBM Development

The case of multicoal seam development is discussed in this section. According to the pressure balance principle, due to the pressure gradient in multicoal seam development, if the pressure of the wellbore is higher than that of the coal seam, intrusion of the fluid would happen; namely, the fluid will flow into the coal seam from the wellbore; otherwise, the coal seam produces gas.Reflected into the model proposed in this paper, as shown in Figure6, the two-phase flow for multicoal seam development was studied. Assume that the fluid pressure of coal seam 1 is greater than that of coal seam 2, namely, p1>p2, in that way, the fluid output from coal seam 1 will flow into coal seam 2. For the modeling of fluid pressure from coal seam 1, it meets the solution of definite problems (III) and (IV); but for coal seam 2, it is equivalent for fluid flowing into the coal seam, with the flow rate
(69)Δqw=k1μ1w⋅πd124p1w+ρg∇H−p2w,Δqg=k1μ1g⋅πd124p1g+ρg∇H−p2g,where Δqw and Δqg are the strength difference of water and gas between the two coal seams, respectively. k1 is the permeability of coal seam 1. d1 is the diameter of the wellbore. μ1g and μ1w are the viscosities of gas and water in the coal seam 1, respectively. p1w and p1g are the pressures of water and gas in the coal seam 1, respectively. p2w and p2g are the pressures of water and gas in the coal seam 2, respectively.Figure 6
Schematic for multicoal seam development of TS-1 well.


### 3.1. Solving Process

Definite problems (VI) and (VII) of the two-phase flow in coal seam 1 are constructed as follows:(70)VIpgt−η12pgrr=0pgt=0=p0pgR,t=p0,pgrr=0=p2g−ρgg∇H−p1g0<r<R,t>0,0≤r≤R,R→∞,where
(71)η12=kkrwρw/μw+Rswkkrgρg/μgϕRswswρw/Bw+sgρg/Bg+α/λ+2G,Δqg=k1/μ1⋅πd12/4p2g−ρgg∇H−p1g..
(72)VIIpwt−η12pwrr=00<r<R,t>0,pwt=0=p00≤r≤R,pwR,t=p0,pwrr=0=p2w−ρwg∇H−p1wR→∞,where η22=kkrw/μwϕsw⋅λ+2G/α and Δqw=k1/μ1⋅πd12/4p2w′−p1w.To simplify the calculation, the similar definite problem (VIII) with the above problems is built.(73)VIIIpt−η2prr=00<r<R,t>0,pt=0=p00≤r≤R,pR,t=p0,prr=0=p2−ρg∇H−pFirstly, aiming to construct a function both meeting equation and boundary conditions, the boundary conditions will experience partial homogeneous treatment.(74)wr,t=p0+p2−ρg∇Hr.Substituting Equation (61) into the original problem, the conversion of function is made as
(75)pr,t=vr,t+wr,t.Thus,vr,t is converted into a mixed problem as
(76)IXvt−η2vrr=00<r<R,t>0,vt=0=−p2−ρg∇Hr0≤r≤R,pR,t=0,prr=0+p=0.Separation of variables is used to solve the above problem (IX), whose detailed solution procedure is similar to problem (VI). Eventually, the solution is(77)vr,t=∑k=1∞1Mk∫0r−p2−ρg∇Hξsinλkξe−η2λktsinλkr,where λ is the solution of the triangulation equation of tgυ=−υ/r, and eigenvalues can be obtained as
(78)λk=υkr2k=1,2,⋯.Then, the solution for the definite problem (IX) is(79)pr,t=p0+p2−ρg∇Hr+∑k=1∞1Mk∫0r−p2−ρg∇Hξsinλkξe−η2λktsinλkr.Corresponding to the solution of Equation (79), the expressions of gas pressure and water pressure can be obtained, respectively.For gas flow(80)pgr,t=p0+p2g−ρg∇Hr+∑k=1∞1Mk∫0r−p2g−ρg∇Hξsinλkξe−η12λktsinλkr,where
(81)η12=KKrwρw/μw+RswKKrgρg/μgϕRswswρw/Bw+sgρg/Bg+α/λ+2G.For water flow(82)pwr,t=p0+p2w−ρg∇Hr+∑k=1∞1Mk∫0r−p2w−ρg∇Hξsinλkξe−η22λktsinλkr,where η22=kkrw/μwϕsw⋅λ+2G/α.


### 3.2. Model Verification

#### 3.2.1. Field Data

As Table3 shows, the field data is from TS-1 well of Liupanshui coal mine within Tianpanjiang production area in Western Guizhou, China [10].Table 3
Field data of TS-1 well.
Fracturing segmentSeam numberVertical depth (m)Thickness (m)Reservoir pressure (MPa)Gas content (m3·t-1)Gas saturationWell testing permeability (mD)Critical desorption pressure (MPa)352531.82.686.9910.6394%0.095.759559.61.397.0212.0086%4.50Table4 shows physical parameters of coal bed, which is from the field data of No. 9 coal seam; various physical parameters will be validated with the theoretical model.Table 4
Physical parameters of coal bed.
ParameterValueYoung’s modulus of coal (E, GPa)2.735Poisson’s ratio of coalυ0.25Coal density (ρ, kg/m3)1500Seam depth (h, m)521.6Seam thickness (H, m)5.8Porosityϕ3.51%Permeability (K, mD)0.09Seam initial pressure (p, MPa)6.99Water viscosity in the ground condition (μw, Pa·s)5.8×10−4Gas viscosity under standard atmospheric pressure (μgsc, Pa·s)1.47×10−5Gas density under standard atmospheric pressure (ρgsc, kg/m3)0.425Water density under standard condition (ρw, kg/m3)0.998Formation volume factor of gasBg0.00896Formation volume factor of waterBw1.0The initial gas saturationsg0.73Water saturationsw0.27SolubilityRsw1.3Relative permeability of gasKrg0.3466Relative permeability of waterKrw0.0284Seam radius (R, m)300


#### 3.2.2. The Analysis of Noninterfering Coal Seam Mining

The case of noninterfering coal seam mining is discussed in this section. Firstly, each parameter value of both coal seams is given in Table4. Take r=0, that is, to research the area near wellhead. Besides, it is given that the gas desorption begins from the 80th day. Without considering the interaction between the two coal seams, combined with Equations (62) and (64), the pressure of the two-phase flow can be simplified into
(83)p1g=6.348e−9.6733×10−6tt>80,n!r!n−r!p1w=6.348e−3.0915×10−4tt>80,p2g=4.968e−1.0545×10−5tt>80,p2w=4.968e−7.2133×10−4tt>80,,where the subscript 1g represents gas phase of No. 52 seam, the subscript 1w represents water phase of No. 52 seam, the subscript 2g represents gas phase of No. 9 seam, and the subscript 2w represents water phase of No. 9 seam.Figure7 shows gas and water pressure curves of both the No. 52 seam and No. 9 seam in single coal seam mining. The four curves are given by the analytical solution above. As is shown in the figure, firstly, whatever the coal seam is, the gas pressure is always decreasing from the start of gas production, and the gas pressure of No. 52 seam is always greater than that of No. 9 seam; secondly, the solution of water pressure of both coal seams in the fluctuation desorption stage is not studied. For No. 52 seam, we can find that the water pressure is not always making the same trend during different water flow stages. And finally, in the 500th day, the water pressure of No. 52 seam drops to 4.75 MPa. In addition, for No. 9 seam, the water pressure maintains the downward trend, and the water pressure of No. 9 seam drops to 4.39 MPa in the 500th day.Figure 7
Fluid pressure changes of two coal seams in single coal seam mining.


#### 3.2.3. The Effects of Factors on the Development of Multicoal Seam

The effects of factors on the development of multicoal seam are discussed in this section. The impact of altitude difference between the two coal seams is considered, and the remaining factors were not considered temporarily.Due to the two-phase flow phenomenon in the CBM process, it should be emphasized that the two coal seams interact each other and intrusion of the fluid occurs between the two coal seams in multicoal seam development in the case that(1)
If(84)p2−ρg∇H>p1,the fluid of No. 9 coal seam flows to the wellhead; meanwhile, a portion of the fluid flows backward into No. 52 coal seam.
(2)
If(85)p1+ρg∇H>p2,a portion of the fluid of No. 52 coal seam will not flow to the wellhead but flows vertically downward and backward into the No. 9 coal seam.Figure8 demonstrates the relationship of pressure difference and height difference of the two seams, both of which are the key factors to influence the development of multicoal seam. As is shown, for water flow, the red solid line represents that there are no influence of water flow between the two coal seams; namely, both the No. 9 seam and No. 52 seam produce water separately and have no effect on the other coal seam. In region ①, water flows from No. 9 seam into No. 52 seam; however, in region ②+③, gas and water from No. 52 seam would flow into No. 9 seam in the mining process. For gas flow, the blue dotted line means that there is no influence of gas flow between the two coal seams. In region ①+②, gas flows from No. 9 seam into No. 52 seam; however, gas would flow from No. 52 seam into No. 9 seam in region ③.Figure 8
Relationship of pressure difference and height difference between the two coal seams.As Figures7 and 8 show, the fluid pressure of No. 52 seam is always greater than that of No. 9 seam. By calculation of Equations (78)–(81), it can be obtained that the relationship of the pressure satisfies Equation (85), which means a part of the gas and water would flow backward into the No. 9 coal seam in the mining process.


#### 3.2.4. The Analysis of Multicoal Seam Mining

As Figure6 shows, the field data is from TS-1 well of Liupanshui coal mine; a portion of the fluid (gas and water) of No. 52 coal seam will not flow to the wellhead but flows vertically downward and backward into the No. 9 coal seam.Substituting the parameter value of Tables3and 4 into Equation (80), the expressions of the gas flow can be given as
(86)pt=1.104p0+p2−ρg∇H−p1e−2.4048η/R2t.Plotting according to expression from Equation (86), the contrast curves between the development of single coal seam and multicoal seam can be obtained as follows.As Figure9 shows, the gas pressure decreases with time regardless of whether there is interaction between the two coal seams. In addition, the gas pressure for the development of multicoal seam decreases faster than that of single coal seam. What is more, the gas pressure for the development of multicoal seam keeps higher than that of single coal seam until the 200th day, and it is lower after the 200th day. It is because when the intrusion happens, the pressure gradient increases, which can accelerate the gas production, followed by the gas pressure dropping faster.Figure 9
Effect of multicoal seam development on gas pressure.For water phase, the classification discussion will be used as follows: whent≤632 d, the pressure of No. 52 is greater than that of No. 9 coal seam, and the water intrusion occurs from No. 52 to No. 9 seam; and when t>632 d, water intrusion will happen from the No. 9 to No. 52 seam, which does not match the time scope of t∈0d,500d in this paper.Thus, substituting the parameter values of Tables3 and 4 into Equation (82), the expressions of the water flow can be given as follows:For water flow of No. 52 seam
(87)p1w=6.348e−3.0915×10−4t80<t<500.For water flow of No. 9 seam(88)p2w′=6.348e−1.0305×10−3t+0.865467e−7.2133×10−4t80<t<500.According to Equations (87) and (88), the contrast curves of water pressure between the development of single coal seam and multicoal seam can be obtained as follows.As Figure10 shows, the water pressure decreases with time regardless of whether there is interaction between the two coal seams. In addition, at the two-phase flow stage, the water pressure for the development of multicoal seam decreases faster. What is more, the gas pressure for the development of multicoal seam keeps lower than that of single coal seam in the whole mining of 500 days. It is because a large amount of drainage occurs in single water flow stage, and plenty of water intrusion will affect the drainage of No. 9 coal seam while considering the impact of multicoal seam; however, the amount of drainage is relatively small when the intrusion happens, and the water pressure gradient will increase in No. 9 coal seam, which can accelerate the drainage, followed by the water pressure dropping faster.Figure 10
Effect of multicoal seam development on water pressure.


## 3.1. Solving Process

Definite problems (VI) and (VII) of the two-phase flow in coal seam 1 are constructed as follows:(70)VIpgt−η12pgrr=0pgt=0=p0pgR,t=p0,pgrr=0=p2g−ρgg∇H−p1g0<r<R,t>0,0≤r≤R,R→∞,where
(71)η12=kkrwρw/μw+Rswkkrgρg/μgϕRswswρw/Bw+sgρg/Bg+α/λ+2G,Δqg=k1/μ1⋅πd12/4p2g−ρgg∇H−p1g..
(72)VIIpwt−η12pwrr=00<r<R,t>0,pwt=0=p00≤r≤R,pwR,t=p0,pwrr=0=p2w−ρwg∇H−p1wR→∞,where η22=kkrw/μwϕsw⋅λ+2G/α and Δqw=k1/μ1⋅πd12/4p2w′−p1w.To simplify the calculation, the similar definite problem (VIII) with the above problems is built.(73)VIIIpt−η2prr=00<r<R,t>0,pt=0=p00≤r≤R,pR,t=p0,prr=0=p2−ρg∇H−pFirstly, aiming to construct a function both meeting equation and boundary conditions, the boundary conditions will experience partial homogeneous treatment.(74)wr,t=p0+p2−ρg∇Hr.Substituting Equation (61) into the original problem, the conversion of function is made as
(75)pr,t=vr,t+wr,t.Thus,vr,t is converted into a mixed problem as
(76)IXvt−η2vrr=00<r<R,t>0,vt=0=−p2−ρg∇Hr0≤r≤R,pR,t=0,prr=0+p=0.Separation of variables is used to solve the above problem (IX), whose detailed solution procedure is similar to problem (VI). Eventually, the solution is(77)vr,t=∑k=1∞1Mk∫0r−p2−ρg∇Hξsinλkξe−η2λktsinλkr,where λ is the solution of the triangulation equation of tgυ=−υ/r, and eigenvalues can be obtained as
(78)λk=υkr2k=1,2,⋯.Then, the solution for the definite problem (IX) is(79)pr,t=p0+p2−ρg∇Hr+∑k=1∞1Mk∫0r−p2−ρg∇Hξsinλkξe−η2λktsinλkr.Corresponding to the solution of Equation (79), the expressions of gas pressure and water pressure can be obtained, respectively.For gas flow(80)pgr,t=p0+p2g−ρg∇Hr+∑k=1∞1Mk∫0r−p2g−ρg∇Hξsinλkξe−η12λktsinλkr,where
(81)η12=KKrwρw/μw+RswKKrgρg/μgϕRswswρw/Bw+sgρg/Bg+α/λ+2G.For water flow(82)pwr,t=p0+p2w−ρg∇Hr+∑k=1∞1Mk∫0r−p2w−ρg∇Hξsinλkξe−η22λktsinλkr,where η22=kkrw/μwϕsw⋅λ+2G/α.


## 3.2. Model Verification

### 3.2.1. Field Data

As Table3 shows, the field data is from TS-1 well of Liupanshui coal mine within Tianpanjiang production area in Western Guizhou, China [10].Table 3
Field data of TS-1 well.
Fracturing segmentSeam numberVertical depth (m)Thickness (m)Reservoir pressure (MPa)Gas content (m3·t-1)Gas saturationWell testing permeability (mD)Critical desorption pressure (MPa)352531.82.686.9910.6394%0.095.759559.61.397.0212.0086%4.50Table4 shows physical parameters of coal bed, which is from the field data of No. 9 coal seam; various physical parameters will be validated with the theoretical model.Table 4
Physical parameters of coal bed.
ParameterValueYoung’s modulus of coal (E, GPa)2.735Poisson’s ratio of coalυ0.25Coal density (ρ, kg/m3)1500Seam depth (h, m)521.6Seam thickness (H, m)5.8Porosityϕ3.51%Permeability (K, mD)0.09Seam initial pressure (p, MPa)6.99Water viscosity in the ground condition (μw, Pa·s)5.8×10−4Gas viscosity under standard atmospheric pressure (μgsc, Pa·s)1.47×10−5Gas density under standard atmospheric pressure (ρgsc, kg/m3)0.425Water density under standard condition (ρw, kg/m3)0.998Formation volume factor of gasBg0.00896Formation volume factor of waterBw1.0The initial gas saturationsg0.73Water saturationsw0.27SolubilityRsw1.3Relative permeability of gasKrg0.3466Relative permeability of waterKrw0.0284Seam radius (R, m)300


### 3.2.2. The Analysis of Noninterfering Coal Seam Mining

The case of noninterfering coal seam mining is discussed in this section. Firstly, each parameter value of both coal seams is given in Table4. Take r=0, that is, to research the area near wellhead. Besides, it is given that the gas desorption begins from the 80th day. Without considering the interaction between the two coal seams, combined with Equations (62) and (64), the pressure of the two-phase flow can be simplified into
(83)p1g=6.348e−9.6733×10−6tt>80,n!r!n−r!p1w=6.348e−3.0915×10−4tt>80,p2g=4.968e−1.0545×10−5tt>80,p2w=4.968e−7.2133×10−4tt>80,,where the subscript 1g represents gas phase of No. 52 seam, the subscript 1w represents water phase of No. 52 seam, the subscript 2g represents gas phase of No. 9 seam, and the subscript 2w represents water phase of No. 9 seam.Figure7 shows gas and water pressure curves of both the No. 52 seam and No. 9 seam in single coal seam mining. The four curves are given by the analytical solution above. As is shown in the figure, firstly, whatever the coal seam is, the gas pressure is always decreasing from the start of gas production, and the gas pressure of No. 52 seam is always greater than that of No. 9 seam; secondly, the solution of water pressure of both coal seams in the fluctuation desorption stage is not studied. For No. 52 seam, we can find that the water pressure is not always making the same trend during different water flow stages. And finally, in the 500th day, the water pressure of No. 52 seam drops to 4.75 MPa. In addition, for No. 9 seam, the water pressure maintains the downward trend, and the water pressure of No. 9 seam drops to 4.39 MPa in the 500th day.Figure 7
Fluid pressure changes of two coal seams in single coal seam mining.


### 3.2.3. The Effects of Factors on the Development of Multicoal Seam

The effects of factors on the development of multicoal seam are discussed in this section. The impact of altitude difference between the two coal seams is considered, and the remaining factors were not considered temporarily.Due to the two-phase flow phenomenon in the CBM process, it should be emphasized that the two coal seams interact each other and intrusion of the fluid occurs between the two coal seams in multicoal seam development in the case that(1)
If(84)p2−ρg∇H>p1,the fluid of No. 9 coal seam flows to the wellhead; meanwhile, a portion of the fluid flows backward into No. 52 coal seam.
(2)
If(85)p1+ρg∇H>p2,a portion of the fluid of No. 52 coal seam will not flow to the wellhead but flows vertically downward and backward into the No. 9 coal seam.Figure8 demonstrates the relationship of pressure difference and height difference of the two seams, both of which are the key factors to influence the development of multicoal seam. As is shown, for water flow, the red solid line represents that there are no influence of water flow between the two coal seams; namely, both the No. 9 seam and No. 52 seam produce water separately and have no effect on the other coal seam. In region ①, water flows from No. 9 seam into No. 52 seam; however, in region ②+③, gas and water from No. 52 seam would flow into No. 9 seam in the mining process. For gas flow, the blue dotted line means that there is no influence of gas flow between the two coal seams. In region ①+②, gas flows from No. 9 seam into No. 52 seam; however, gas would flow from No. 52 seam into No. 9 seam in region ③.Figure 8
Relationship of pressure difference and height difference between the two coal seams.As Figures7 and 8 show, the fluid pressure of No. 52 seam is always greater than that of No. 9 seam. By calculation of Equations (78)–(81), it can be obtained that the relationship of the pressure satisfies Equation (85), which means a part of the gas and water would flow backward into the No. 9 coal seam in the mining process.


### 3.2.4. The Analysis of Multicoal Seam Mining

As Figure6 shows, the field data is from TS-1 well of Liupanshui coal mine; a portion of the fluid (gas and water) of No. 52 coal seam will not flow to the wellhead but flows vertically downward and backward into the No. 9 coal seam.Substituting the parameter value of Tables3and 4 into Equation (80), the expressions of the gas flow can be given as
(86)pt=1.104p0+p2−ρg∇H−p1e−2.4048η/R2t.Plotting according to expression from Equation (86), the contrast curves between the development of single coal seam and multicoal seam can be obtained as follows.As Figure9 shows, the gas pressure decreases with time regardless of whether there is interaction between the two coal seams. In addition, the gas pressure for the development of multicoal seam decreases faster than that of single coal seam. What is more, the gas pressure for the development of multicoal seam keeps higher than that of single coal seam until the 200th day, and it is lower after the 200th day. It is because when the intrusion happens, the pressure gradient increases, which can accelerate the gas production, followed by the gas pressure dropping faster.Figure 9
Effect of multicoal seam development on gas pressure.For water phase, the classification discussion will be used as follows: whent≤632 d, the pressure of No. 52 is greater than that of No. 9 coal seam, and the water intrusion occurs from No. 52 to No. 9 seam; and when t>632 d, water intrusion will happen from the No. 9 to No. 52 seam, which does not match the time scope of t∈0d,500d in this paper.Thus, substituting the parameter values of Tables3 and 4 into Equation (82), the expressions of the water flow can be given as follows:For water flow of No. 52 seam
(87)p1w=6.348e−3.0915×10−4t80<t<500.For water flow of No. 9 seam(88)p2w′=6.348e−1.0305×10−3t+0.865467e−7.2133×10−4t80<t<500.According to Equations (87) and (88), the contrast curves of water pressure between the development of single coal seam and multicoal seam can be obtained as follows.As Figure10 shows, the water pressure decreases with time regardless of whether there is interaction between the two coal seams. In addition, at the two-phase flow stage, the water pressure for the development of multicoal seam decreases faster. What is more, the gas pressure for the development of multicoal seam keeps lower than that of single coal seam in the whole mining of 500 days. It is because a large amount of drainage occurs in single water flow stage, and plenty of water intrusion will affect the drainage of No. 9 coal seam while considering the impact of multicoal seam; however, the amount of drainage is relatively small when the intrusion happens, and the water pressure gradient will increase in No. 9 coal seam, which can accelerate the drainage, followed by the water pressure dropping faster.Figure 10
Effect of multicoal seam development on water pressure.


## 3.2.1. Field Data

As Table3 shows, the field data is from TS-1 well of Liupanshui coal mine within Tianpanjiang production area in Western Guizhou, China [10].Table 3
Field data of TS-1 well.
Fracturing segmentSeam numberVertical depth (m)Thickness (m)Reservoir pressure (MPa)Gas content (m3·t-1)Gas saturationWell testing permeability (mD)Critical desorption pressure (MPa)352531.82.686.9910.6394%0.095.759559.61.397.0212.0086%4.50Table4 shows physical parameters of coal bed, which is from the field data of No. 9 coal seam; various physical parameters will be validated with the theoretical model.Table 4
Physical parameters of coal bed.
ParameterValueYoung’s modulus of coal (E, GPa)2.735Poisson’s ratio of coalυ0.25Coal density (ρ, kg/m3)1500Seam depth (h, m)521.6Seam thickness (H, m)5.8Porosityϕ3.51%Permeability (K, mD)0.09Seam initial pressure (p, MPa)6.99Water viscosity in the ground condition (μw, Pa·s)5.8×10−4Gas viscosity under standard atmospheric pressure (μgsc, Pa·s)1.47×10−5Gas density under standard atmospheric pressure (ρgsc, kg/m3)0.425Water density under standard condition (ρw, kg/m3)0.998Formation volume factor of gasBg0.00896Formation volume factor of waterBw1.0The initial gas saturationsg0.73Water saturationsw0.27SolubilityRsw1.3Relative permeability of gasKrg0.3466Relative permeability of waterKrw0.0284Seam radius (R, m)300


## 3.2.2. The Analysis of Noninterfering Coal Seam Mining

The case of noninterfering coal seam mining is discussed in this section. Firstly, each parameter value of both coal seams is given in Table4. Take r=0, that is, to research the area near wellhead. Besides, it is given that the gas desorption begins from the 80th day. Without considering the interaction between the two coal seams, combined with Equations (62) and (64), the pressure of the two-phase flow can be simplified into
(83)p1g=6.348e−9.6733×10−6tt>80,n!r!n−r!p1w=6.348e−3.0915×10−4tt>80,p2g=4.968e−1.0545×10−5tt>80,p2w=4.968e−7.2133×10−4tt>80,,where the subscript 1g represents gas phase of No. 52 seam, the subscript 1w represents water phase of No. 52 seam, the subscript 2g represents gas phase of No. 9 seam, and the subscript 2w represents water phase of No. 9 seam.Figure7 shows gas and water pressure curves of both the No. 52 seam and No. 9 seam in single coal seam mining. The four curves are given by the analytical solution above. As is shown in the figure, firstly, whatever the coal seam is, the gas pressure is always decreasing from the start of gas production, and the gas pressure of No. 52 seam is always greater than that of No. 9 seam; secondly, the solution of water pressure of both coal seams in the fluctuation desorption stage is not studied. For No. 52 seam, we can find that the water pressure is not always making the same trend during different water flow stages. And finally, in the 500th day, the water pressure of No. 52 seam drops to 4.75 MPa. In addition, for No. 9 seam, the water pressure maintains the downward trend, and the water pressure of No. 9 seam drops to 4.39 MPa in the 500th day.Figure 7
Fluid pressure changes of two coal seams in single coal seam mining.


## 3.2.3. The Effects of Factors on the Development of Multicoal Seam

The effects of factors on the development of multicoal seam are discussed in this section. The impact of altitude difference between the two coal seams is considered, and the remaining factors were not considered temporarily.Due to the two-phase flow phenomenon in the CBM process, it should be emphasized that the two coal seams interact each other and intrusion of the fluid occurs between the two coal seams in multicoal seam development in the case that(1)
If(84)p2−ρg∇H>p1,the fluid of No. 9 coal seam flows to the wellhead; meanwhile, a portion of the fluid flows backward into No. 52 coal seam.
(2)
If(85)p1+ρg∇H>p2,a portion of the fluid of No. 52 coal seam will not flow to the wellhead but flows vertically downward and backward into the No. 9 coal seam.Figure8 demonstrates the relationship of pressure difference and height difference of the two seams, both of which are the key factors to influence the development of multicoal seam. As is shown, for water flow, the red solid line represents that there are no influence of water flow between the two coal seams; namely, both the No. 9 seam and No. 52 seam produce water separately and have no effect on the other coal seam. In region ①, water flows from No. 9 seam into No. 52 seam; however, in region ②+③, gas and water from No. 52 seam would flow into No. 9 seam in the mining process. For gas flow, the blue dotted line means that there is no influence of gas flow between the two coal seams. In region ①+②, gas flows from No. 9 seam into No. 52 seam; however, gas would flow from No. 52 seam into No. 9 seam in region ③.Figure 8
Relationship of pressure difference and height difference between the two coal seams.As Figures7 and 8 show, the fluid pressure of No. 52 seam is always greater than that of No. 9 seam. By calculation of Equations (78)–(81), it can be obtained that the relationship of the pressure satisfies Equation (85), which means a part of the gas and water would flow backward into the No. 9 coal seam in the mining process.


## 3.2.4. The Analysis of Multicoal Seam Mining

As Figure6 shows, the field data is from TS-1 well of Liupanshui coal mine; a portion of the fluid (gas and water) of No. 52 coal seam will not flow to the wellhead but flows vertically downward and backward into the No. 9 coal seam.Substituting the parameter value of Tables3and 4 into Equation (80), the expressions of the gas flow can be given as
(86)pt=1.104p0+p2−ρg∇H−p1e−2.4048η/R2t.Plotting according to expression from Equation (86), the contrast curves between the development of single coal seam and multicoal seam can be obtained as follows.As Figure9 shows, the gas pressure decreases with time regardless of whether there is interaction between the two coal seams. In addition, the gas pressure for the development of multicoal seam decreases faster than that of single coal seam. What is more, the gas pressure for the development of multicoal seam keeps higher than that of single coal seam until the 200th day, and it is lower after the 200th day. It is because when the intrusion happens, the pressure gradient increases, which can accelerate the gas production, followed by the gas pressure dropping faster.Figure 9
Effect of multicoal seam development on gas pressure.For water phase, the classification discussion will be used as follows: whent≤632 d, the pressure of No. 52 is greater than that of No. 9 coal seam, and the water intrusion occurs from No. 52 to No. 9 seam; and when t>632 d, water intrusion will happen from the No. 9 to No. 52 seam, which does not match the time scope of t∈0d,500d in this paper.Thus, substituting the parameter values of Tables3 and 4 into Equation (82), the expressions of the water flow can be given as follows:For water flow of No. 52 seam
(87)p1w=6.348e−3.0915×10−4t80<t<500.For water flow of No. 9 seam(88)p2w′=6.348e−1.0305×10−3t+0.865467e−7.2133×10−4t80<t<500.According to Equations (87) and (88), the contrast curves of water pressure between the development of single coal seam and multicoal seam can be obtained as follows.As Figure10 shows, the water pressure decreases with time regardless of whether there is interaction between the two coal seams. In addition, at the two-phase flow stage, the water pressure for the development of multicoal seam decreases faster. What is more, the gas pressure for the development of multicoal seam keeps lower than that of single coal seam in the whole mining of 500 days. It is because a large amount of drainage occurs in single water flow stage, and plenty of water intrusion will affect the drainage of No. 9 coal seam while considering the impact of multicoal seam; however, the amount of drainage is relatively small when the intrusion happens, and the water pressure gradient will increase in No. 9 coal seam, which can accelerate the drainage, followed by the water pressure dropping faster.Figure 10
Effect of multicoal seam development on water pressure.


## 4. Conclusions

Two models of two-phase flow for the development of single coal seam layer and multicoal seam layers were developed, respectively, to simulate the fluid flow in coal seam gas reservoirs in the mining process. These models were then analytically solved by separation of variables and applied to analyze the field data from different CBM formations (the 1# test well of Jincheng and the TS-1 well of Liupanshui coal mine, China). The effects of gas solubility on the gas pressure were investigated. Besides, the impact of interlayer difference on the development of multicoal seam layers was explored. Based on these studies, the following conclusions can be made.Firstly, the two-phase models for the development of single coal seam and multicoal seam, respectively, are proved to be correct after matching the field data well, respectively. Moreover, the gas pressure is influenced by the gas solubility; if the value of gas solubility increases, the gas pressure will drop increasingly.Secondly, the key factors incorporating the height difference and the pressure difference of the two seams are proved to influence the multicoal seam development comprehensively. For the field from the TS-1 well of Liupanshui coal mine, during the seal layer mining process, the pressure of the upper No. 52 seam is greater than that of the lower No. 9 seam; thus, the gas and water will intrude into the No. 9 coal seam and stimulate the production and overall shorten the mining time.



---

*Source: 6690218-2021-02-16.xml*