Permeability Model of Coal Measure Gas Reservoirs Considering Dynamic Diffusion
-
摘要: 为了预测煤系气开采过程中储层渗透率的演化规律,考虑气体在基质中的动力学扩散作用,基于储层的应力-应变本构关系和渗透率-孔隙率的立方关系,提出了储层的有效应力-渗透率模型,分别建立了储层在常体积和单轴应变条件下的渗透率解析模型。利用现场和实验室测试的渗透率数据,探讨了两种模型的有效性。结果表明,与常体积条件下的渗透率模型和C-M模型相比,单轴应变条件下的渗透率模型能够较好地拟合现场和实验室的渗透率数据,建立渗透率模型须考虑气体在基质中的动力学扩散作用。研究了模型参数对渗透率的影响,结果表明模型参数对渗透率的演化规律以及反弹压力有显著的影响。Abstract: In order to predict evolution of permeability for reservoir in the process of coal measure gas exploitation, based on the stress-strain constitutive of reservoir and the cubical relation of permeability and porosity, the model of the effective stress-permeability of reservoir was presented, which considers the kinetic diffusion of gases in the matrix. The analytical models of reservoir under constant volume and uniaxial strain conditions were established respectively. Furthermore, the effectiveness of the two models are investigated using permeability data from field and laboratory tests respectively. The results show that, compared with the model under constant volume condition and C-M model, the permeability model under uniaxial strain condition can better fit the permeability data from the field and laboratory. It is very important to take the dynamic diffusion of gas in the matrix into consideration during establishing the permeability model. Besides, the effect of model parameters on permeability was studied. It is shown that the model parameters have a significant influence on the permeability evolution and rebound pressure.
-
煤系气泛指煤系中赋存的各类天然气,主要包括煤层气、页岩气及煤系砂岩气等[1]。截至2017年底,我国煤层气累计提交探明储量为7.118 × 1011 m3,2017年产气量为4.95 × 109 m3[2–3]。我国页岩气技术可采资源量约为2.18 × 1013 m3[4],煤系页岩气占页岩气总资源量的19%~37%,截至2019年,我国页岩气年产气量达到1.53 × 1010 m3[5]。目前,我国页岩气勘探开发的主体为海相页岩气和陆相页岩气。规模开发煤系气,能够优化调整我国的能源结构。
煤系气储层是典型的天然裂缝性储层,渗透率是评价其商业开发可行性的重要参数之一。目前,国内外学者对煤系气储层渗透率解析模型的研究成果颇丰。Palmer等[6]认为储层压力降低导致渗透率降低,而储层压降又导致甲烷解吸,煤基质收缩,裂隙张开度增大,渗透率随之增大,在此基础上提出了计算煤层渗透率随有效应力和煤基质收缩变化的理论模型,该模型适用于单轴应变条件。Shi等[7]假设储层水平主应力垂直于裂隙,通过直接类比热收缩和基质收缩,得到等温气体解吸煤层的本构关系,基于此,推导了煤层渗透率模型。Cui等[8]利用线性等温多孔弹性理论,结合孔隙率的定义,假设煤层处于单轴应变条件,推导了煤层的渗透率变化模型。张宏学等[9]通过引入裂隙法向刚度,建立了页岩储层的有效应力-渗透率模型,进一步分析了页岩储层在单轴应变和常体积条件下的渗透率模型。文献[10–11]均将煤层简化为双孔双渗模型(基质和裂隙),分别提出了煤基质系统和裂隙系统的渗透率模型。文献[12–14]将煤系气储层简化为三孔双渗模型,分别建立了储层的渗透率解析模型。
上述文献在建立储层的渗透率模型时,均认为在开采过程中储层的裂隙压力始终等于基质压力,但是在实际开采过程中,由于气体在基质中的动力学扩散作用,储层基质压力不等于裂隙压力,尤其在开采初期,这一现象特别明显。本研究将煤系气储层简化为双孔单渗模型,基于多孔弹性力学理论以及渗透率与孔隙率的立方关系,考虑气体在基质中的动力学扩散作用,提出与应力相关的渗透率模型,在此基础上,分别得到煤系气储层在常体积和单轴应变条件下的渗透率模型,并利用现场实测和实验室测试的渗透率数据验证模型的有效性,并分析模型中的各参数对渗透率的影响。
1. 渗透率模型
1.1 有效应力-渗透率模型
将煤系气储层简化为双孔单渗模型,如图1所示。该模型中基质渗透率近似为零,吸附在基质上的气体解吸以后向裂缝扩散,通过裂缝流入生产井。设压应力和压应变为正,考虑气体在基质中的动力学扩散作用,煤系气储层的应力-应变本构关系为
σij=λεkkδij+2Gεij+αpmδij+βpfδij+Kmεsδij (1) 式中:σij为应力分量,Pa;
λ 为拉梅常数,Pa,λ= Eμ/[(1+μ)(1−2μ)] ,其中E为岩体的弹性模量,Pa,μ 为泊松比;εkk 为体应变,εkk =εx +εy +εz ;δij 为Kronecker符号;G为剪切模量,Pa,G=E/[2(1+μ)] ;εij 为应变分量;Km和Ks分别为基质和固体颗粒的体积模量,Pa,Km=Em/[3(1−2μ)] ,Em为基质的弹性模量,Pa;α 和β 均为Boit系数,α=1−Km/Ks ,β=1−K/Km ,K为体积模量,Pa;pm和pf分别为基质和裂隙的压力,Pa;εs 为吸附/解吸引起的基质的体应变,εs=εLpm/(pm+pL) ,εL 为朗缪尔体积应变,pL 为朗缪尔压力,Pa。煤系气储层和基质系统的体积分别为
V=Vm+Vf (2) Vm=Vp+Vs (3) 式中:V、Vm、Vs、Vf和Vp分别为储层体积、基质体积、固体颗粒体积、裂隙体积和基质中孔隙的体积,m3。
储层的体应变增量为
dεb=−dVV=1K(d¯σ−αdpm−βdpf)−dεs (4) 式中:d
εb 为储层体应变增量;¯σ 为平均应力,Pa。基质孔隙和裂隙孔隙的体应变增量分别为
dεp=−dVpVp=1Kpd¯σ−(1Kp−1Ks)dpm−dεs (5) dεf=−dVfVf=1Kfd¯σ−(1Kf−1Km)dpf (6) 式中:d
εp 和dεf 分别为基质孔隙和裂隙孔隙的体应变增量;Kp和Kf分别为孔隙体积模量和裂隙体积模量,Pa,Kp =ϕ 0K,ϕ 0为初始孔隙率。煤系气储层的孔隙率为
ϕ=Vp+VfV (7) dϕ=d(Vp+VfV)=Vp+VfV[d(Vp+Vf)Vp+Vf−dVV] (8) 式中:
ϕ 为孔隙率。由于Vm >> Vf,根据基质孔隙率和裂隙孔隙率的定义,对式(7)和式(8)进行化简,得到
dϕϕ=ϕmϕdVpVp+ϕfϕdVfVf−dVV (9) 式中:
ϕm 和ϕf 分别为基质孔隙率和裂隙孔隙率。将式(4)、式(5)和式(6)代入式(9),得到
dϕϕ=(1K−ϕmϕKp−ϕfϕKf)d¯σ+[ϕmϕ(1Kp−1Ks)−αK]dpm+[ϕfϕ(1Kf−1Km)−βK]dpf+(ϕmϕ−1)dεs (10) 对于双孔单渗模型,基质中孔隙体积比裂隙体积大若干个数量级,简化式(10),得到
dϕϕ=(1K−1Kp)d¯σ+(1Kp−1Ks−αK)dpm−βKdpf (11) 假设K、Km、Ks、Kp和Kf为常数,对式(11)进行积分,得到
ϕϕ0=exp[(1K−1Kp)(¯σ−¯σ0)+(1Kp−1Ks−αK)(pm−pm0)−βK(pf−pf0)] (12) 式中:
¯σ0 为初始平均应力,Pa;pm0和pf0分别为初始基质压力和初始裂隙压力,Pa。由于Km >> K >> Kf,Ks >> Km >> Kp,K >> Kp,基于渗透率和孔隙率之间的立方关系,简化式(12),得到煤系气储层与有效应力相关的渗透率模型
kk0=exp[−3Kp(¯σ−¯σ0)+3Kp(pm−pm0)−3K(pf−pf0)] (13) 式中:k和k0分别为储层渗透率和初始渗透率。
若煤系气储层裂隙压力和基质压力相等,即pm = pf = p,则式(13)简化为与文献[8]中相同的有效应力-渗透率模型
kk0=exp{−3Kp[(¯σ−¯σ0)−(p−p0)]} (14) 1.2 常体积条件下的渗透率模型
在常体积条件下,储层沿着竖直和水平方向的应变增量均为零,即d
ε x = dε y = dε z = 0,由于Ks >> Km >> K(α=β=1 ),由式(1)得到储层平均应力的变化量为¯σ−¯σ0=(pm−pm0)+(pf−pf0)+Km(εs−εs0) (15) 式中:
εs0 为初始孔隙压力下气体吸附引起的体应变。将式(15)代入式(13),且K >> Kp,得到煤系气储层在常体积条件下的渗透率解析模型
kk0=exp{−3Kp[(pf−pf0)+Km(εLpmpm+pL−εLpm0pm0+pL)]} (16) 煤系气储层在常体积条件下的渗透率模型表明,储层渗透率与裂隙压力、基质压力和气体解吸/吸附有关。随着煤系气的开采,裂隙压力降低,作用在裂隙上的有效应力增大,裂隙张开度减小,渗透率降低;另一方面,基质压力降低,吸附气解吸导致基质收缩,裂隙张开度增大,渗透率增大。
若煤系气储层裂隙压力和基质压力相等,则煤系气储层在常体积条件下的渗透率解析模型简化为
kk0=exp{−3Kp[(p−p0)+Km(εLpp+pL−εLp0p0+pL)]} (17) 1.3 单轴应变条件下的渗透率模型
假设煤系气储层处于单轴应变条件,即
εx =εy = 0,σz 为常数,且Ks >> Km >> K(α =β = 1),由式(1)可得页岩储层水平应力为σx=σy=μ1−μσz+1−2μ1−μpm+1−2μ1−μpf+1−2μ1−μKmεs (18) 由式(18)可得平均法向应力的变化量为
¯σ−¯σ0=2(1−2μ)3(1−μ)[(pm−pm0)+(pf−pf0)+Km(εs−εs0)] (19) 将式(19)代入式(13),得到煤系气储层的渗透率模型为
kk0=exp{3Kp[1+μ3(1−μ)(pm−pm0)−[2(1−2μ)3(1−μ)+KpK](pf−pf0)−2Em9(1−μ)(εs−εs0)]} (20) 由于
KpK<<2(1−2μ)3(1−μ) ,因此式(20)简化为kk0=exp{3Kp[1+μ3(1−μ)(pm−pm0)−2(1−2μ)3(1−μ)(pf−pf0)−2Em9(1−μ)(εs−εs0)]} (21) 将吸附引起的体应变代入式(21),即可得到煤系气储层在单轴应变条件下的渗透率解析模型
kk0=exp{3Kp[1+μ3(1−μ)(pm−pm0)−2(1−2μ)3(1−μ)(pf−pf0)−2Em9(1−μ)(εLpmpm+pL−εLpm0pm0+pL)]} (22) 煤系气储层渗透率解析模型中的第1项表示基质压力变化对渗透率的影响,第2项表示裂隙压力变化对渗透率的影响,第3项表示煤系气解吸/吸附对渗透率的影响。该渗透率模型认为储层基质压力和裂隙压力不相等,因此,式(19)表示的渗透率模型比C-M模型多了一项。与常体积条件下的渗透率解析模型相比,单轴应变条件下的渗透率模型考虑了基质压力变化对渗透率的影响。
若煤系气储层裂隙压力和基质压力相等,则式(19)所示的渗透率模型简化为
kk0=exp{3Kp[5μ−13(1−μ)(p−p0)−2Em9(1−μ)(εLpp+pL−εLp0p0+pL)]} (23) 1.4 反弹压力
由于裂隙压缩和基质收缩对渗透率的影响是相反的,因此随着煤系气的开采,当储层压力下降到某一数值时,煤系气储层的渗透率会出现反弹,渗透率反弹时的储层压力称为反弹压力,当渗透率反弹到初始渗透率时的储层压力称为开采压力。由式(23)可得开采压力prc和反弹压力prb分别为
prc=2EmεLpL3(5μ−1)(p0+pL)−pL (24) prb=√2EmεLpL3(5μ−1)−pL (25) 式(25)表明,反弹压力与基质弹性模量、泊松比、朗缪尔体应变和朗缪尔压力有关。各参数取值不同时,反弹压力随弹性模量的演化规律如图2所示。
由图2可知,当其他参数不变时,反弹压力随弹性模量的增大而增大。各参数相同时,本研究的反弹压力大于C-M模型的反弹压力。当其他参数不变时,反弹压力随着泊松比的减小而增大,随着朗缪尔压力的减小而减小。
2. 模型有效性分析
2.1 模型验证
为了利用渗透率解析模型评价煤系气藏商业开发的可行性,解决关于煤系气开采过程中储层为常体积条件还是单轴应变条件的争议,必须利用实验数据验证模型的有效性和准确性。本研究分别建立了储层在常体积条件和单轴应变条件下的渗透率模型,分别利用San Juan盆地现场实测的渗透率数据以及室内实验得到的渗透率数据对两种条件下的渗透率模型进行拟合,如图3所示,其中:Model 1、Model 2分别表示单轴应变模型和常体积模型。
由图3可知,当常体积条件和单轴应变条件下的渗透率模型中的参数相同时,如表1[15-16]和表2[8, 16]所示,模型中的初始渗透率
k0=b30/12a0 ,初始孔隙率ϕ 0 =2b0/a0, 其中:a0为初始裂隙间距,m;b0为初始裂隙张开度,m。单轴应变条件下的渗透率模型能够较好地反映San Juan盆地以及室内实验岩心的渗透率演化规律,因此煤系气储层在单轴应变条件下的渗透率解析模型是有效、准确的,能够较准确地预测储层以及室内岩心的渗透率演化规律。由图3还可以看出,随着孔隙压力的增大,模型渗透率和现场渗透率的比值逐渐减小。在本模型和C-M模型(没有考虑动力学扩散作用)中的参数取值相同的条件下,本模型所体现的渗透率比值随孔隙压力的演化规律与现场渗透率、室内实验测试的渗透率演化规律基本一致。在孔隙压力增大初期,C-M模型的渗透率比值下降较平缓,而本模型和实验测试的渗透率比值急剧下降。这是由于C-M模型认为基质压力和裂隙压力相等,储层压力增大表示基质压力和裂隙压力均同时增大,作用在基质和裂隙上的有效应力相同且随压力的增大而减小,基质尺寸和裂隙张开度随压力的增大逐渐增大,裂隙张开度增大时,渗透率增大,而基质尺寸增大时,渗透率减小,因为基质弹性模量大于裂隙弹性模量,在相同的有效应力作用下,基质尺寸的增大量小于裂隙张开度的增大量,因此储层渗透率增大。另一方面,基质压力增大导致气体吸附,基质膨胀,裂隙张开度减小,渗透率随之减小。当气体吸附引起的渗透率减小量略大于基质和裂隙变形引起的渗透率增大量时,渗透率随着孔隙压力的增大呈现出平缓的下降趋势。本研究在建立储层的渗透率模型时,认为基质压力与裂隙压力不相等,储层压力增大表示裂隙压力先增大,气体因此向基质流动,基质压力随之增大,在储层压力增大过程中,基质压力小于裂隙压力,因此作用在基质上的有效应力大于作用在裂隙上的有效应力;与C-M模型相比,在同一压力作用下,本模型的裂隙压力不变,裂隙张开度不变,而基质压力降低,基质受到的有效应力增大,基质尺寸的变形量减小,因此由基质尺寸变化导致本模型的渗透率增加量小于C-M模型的渗透率增加量,所以因压力变化导致基质尺寸和裂隙张开度变化而引起的渗透率增加量减小;另一方面,基质压力增大导致气体吸附,基质膨胀,裂隙张开度减小,渗透率随之减小。综上所述,与C-M模型相比,本模型在孔隙压力增大初期,渗透率比值急剧下降。
2.2 渗透率演化规律
对于式(19)所示的煤系气储层渗透率解析模型,基质压力和裂隙压力是不相等的,基质初始压力等于裂隙初始压力。下面将分别讨论渗透率随基质压力(裂隙压力不变)和裂隙压力(基质压力不变)的演化规律。
当裂隙压力不变时,储层渗透率随基质压力的演化规律如图4所示。由图4可以看出,当裂隙压力不变时,随着基质压力的降低,渗透率先平缓地下降,当基质压力下降到一定程度时,渗透率出现反弹,然后迅速增大。裂隙压力越小,随着基质压力的下降,渗透率下降得越急剧,则随后反弹得越剧烈。渗透率之所以随基质压力降低出现反弹现象,是因为基质变形和气体解吸对渗透率的影响是相反的。
当基质压力不变时,储层渗透率随裂隙压力的演化规律如图5所示。由图5可以看出,当基质压力不变时,渗透率随裂隙压力的下降逐渐增大,渗透率没有出现反弹。基质压力越大,渗透率增大得越急剧。随着裂隙压力的降低,渗透率之所以没有出现反弹现象,是因为本研究在建立渗透率模型时不考虑气体对裂隙的吸附作用。
3. 模型参数对渗透率的影响
本研究所建立的煤系气储层渗透率解析模型表明,渗透率的主要影响因素为弹性模量、泊松比、初始孔隙压力、初始渗透率、初始裂隙间距等。下面将分别分析上述因素对储层渗透率的影响。
模型参数对储层渗透率的影响如图6所示。由图6(a)可知,泊松比越小,随着孔隙压力的下降,渗透率变化越平缓,但是渗透率反弹越剧烈。由图6(b)可知,弹性模量越大,随着孔隙压力的下降,渗透率变化得越平缓,但是渗透率反弹越剧烈。由图6(c)可知,初始孔隙压力不同时,渗透率随着孔隙压力的下降呈现出相同的变化趋势,孔隙压力下降初期,渗透率急剧下降,当压力下降到某一数值时,出现急剧的反弹现象。由图6(d)可知,当初始渗透率相同时,初始裂隙间距越大,随着孔隙压力的下降,渗透率下降越剧烈,渗透率反弹越剧烈。由图6(e)可知,当初始裂隙间距相同时,初始渗透率越小,渗透率下降越剧烈,渗透率反弹越剧烈。
综上所述,泊松比越小,储层渗透率反弹越剧烈;弹性模量越大,渗透率反弹越剧烈;初始孔隙压力对渗透率反弹趋势影响较小;当初始渗透率相同时,初始裂隙间距越大,渗透率反弹越剧烈;当初始裂隙间距相同时,初始渗透率越小,渗透率反弹越剧烈。
4. 结 论
(1)考虑气体在基质中的动力学扩散作用,基于多孔弹性力学理论,建立了煤系气储层的有效应力-渗透率解析模型,分别推导了储层在常体积条件和单轴应变条件下的渗透率模型。该模型表明,弹性模量、泊松比、初始孔隙压力、初始渗透率、初始裂隙间距是影响储层渗透率的主要因素。
(2)分别利用San Juan盆地现场实测的渗透率数据以及室内实验测试的渗透率数据验证了模型的有效性。结果表明,当考虑气体在基质中的动力学扩散作用时,单轴应变条件下的渗透率模型随孔隙压力的演化规律和实验测试的渗透率演化规律基本吻合。C-M模型(未考虑气体在基质中的动力学扩散作用)体现的渗透率演化规律与实验测试的渗透率演化规律的误差较大。因此,在建立渗透率模型时,必须考虑气体在基质中的动力学扩散作用,不能简单地认为基质压力等于裂隙压力。
(3)渗透率随着压力的降低而减小,当压力下降到一定程度时,由于基质变形和气体解吸对渗透率的影响是相反的,进而导致渗透率开始出现反弹。
(4)在孔隙压力下降初期,泊松比越大,弹性模量越小;初始裂隙间距越大(初始渗透率不变),初始渗透率越小(初始裂隙间距不变),渗透率下降越剧烈。在孔隙压力下降后期,泊松比越小,弹性模量越大;初始裂隙间距越大(初始渗透率不变),初始渗透率越小(初始裂隙间距不变),渗透率反弹越剧烈。
-
Field experiment Laboratory experiment Reservoir pressure/MPa k/k0 Reservoir pressure/MPa k/k0 0.10 7.30 0.38 12.99 0.35 6.62 0.68 10.51 0.79 5.65 1.39 3.93 1.10 4.95 2.04 2.49 1.38 4.25 3.06 1.68 2.07 3.20 4.11 1.36 2.76 2.45 5.14 1.20 3.45 1.75 6.26 1.04 4.14 1.35 4.83 1.05 5.52 1.02 -
[1] 秦勇. 中国煤系气共生成藏作用研究进展 [J]. 天然气工业, 2018, 38(4): 26–36.QIN Y. Research progress of symbiotic accumulation of coal measure gas in China [J]. Natural Gas Industry, 2018, 38(4): 26–36. [2] 李勇, 王延斌, 孟尚志, 等. 煤系非常规天然气合采地质基础理论进展及展望 [J]. 煤炭学报, 2020, 45(4): 1406–1418.LI Y, WANG Y B, MENG S Z, et al. Theoretical basis and prospect of coal measure unconventional natural gas co-production [J]. Journal of China Coal Society, 2020, 45(4): 1406–1418. [3] 曹毅民, 丁蓉, 赵启阳, 等. 煤层气可采储量计算方法的评价与应用 [J]. 天然气工业, 2020, 38(Suppl 1): 50–56.CAO Y M, DING R, ZHAO Q Y, et al. Evaluation and application of the calculation method of recoverable reserves of coalbed methane [J]. Natural Gas Industry, 2020, 38(Suppl 1): 50–56. [4] 张宏学, 刘卫群. 海陆过度相煤系页岩气的渗流特征 [J]. 高压物理学报, 2018, 32(5): 055901.ZHANG H X, LIU W Q. Seepage of marine-terrigenous facies coal measures shale [J]. Chinese Journal of High Pressure Physics, 2018, 32(5): 055901. [5] 王红岩, 周尚文, 刘德勋, 等. 页岩气地质评价关键实验技术的进展与展望 [J]. 天然气工业, 2020, 40(6): 1–17.WANG H Y, ZHOU S W, LIU D X, et al. Progress and prospect of key experimental technologies for shale gas geological evaluation [J]. Natural Gas Industry, 2020, 40(6): 1–17. [6] PALMER I, MANSOORI J. How permeability depends on stress and pore pressure in coalbeds: a new model [J]. SPE Reservoir Evaluation and Engineering, 1998, 1(6): 539–544. doi: 10.2118/52607-PA [7] SHI J Q, DURUCAN S. Drawdown induced changes in permeability of coalbeds: a new interpretation of the reservoir response to primary recovery [J]. Transport in Porous Media, 2004, 56(1): 1–16. doi: 10.1023/B:TIPM.0000018398.19928.5a [8] CUI X J, BUSTIN R M. Volumetric strain associated with methane desorption and its impact on coalbed gas production from deep coal seams [J]. AAPG Bulletin, 2005, 89(9): 1181–1202. doi: 10.1306/05110504114 [9] 张宏学, 刘卫群, 朱立. 页岩储层裂隙渗透率模型和试验研究 [J]. 岩土力学, 2015, 36(3): 719–729.ZHANG H X, LIU W Q, ZHU L. Fracture permeability model and experiments of shale gas reservoirs [J]. Rock and Soil Mechanics, 2015, 36(3): 719–729. [10] WU Y, LIU J S, ELSWORTH D, et al. Development of anisotropic permeability during coalbed methane production [J]. Journal of Natural Gas Science and Engineering, 2010, 2(4): 197–210. doi: 10.1016/j.jngse.2010.06.002 [11] LIU T, LIN B Q, YANG W, et al. Coal permeability evolution and gas migration non-equilibrium state [J]. Transport in Porous Media, 2017, 118(3): 393–416. doi: 10.1007/s11242-017-0862-8 [12] WEI Z J, ZHANG D X. Coupled fluid-flow and geomechanics for triple-porosity/dual-permeability modeling of coalbed methane recovery [J]. International Journal of Rock Mechanics and Mining Sciences, 2010, 47(8): 1242–1253. doi: 10.1016/j.ijrmms.2010.08.020 [13] SANG G J, ELSWORTH D, MIAO X X, et al. Numerical study of a stress dependent triple porosity model for shale gas reservoirs accommodating gas diffusion in kerogen [J]. Journal of Natural Gas Science and Engineering, 2016, 32: 423–438. doi: 10.1016/j.jngse.2016.04.044 [14] WANG G, WANG K, WANG S G, et al. An improved permeability evolution model and its application in fractured sorbing media [J]. Journal of Natural Gas Science and Engineering, 2018, 56: 222–232. doi: 10.1016/j.jngse.2018.05.038 [15] SHI J Q, DURUCAN S. A model for changes in coalbed permeability during primary and enhanced methane recovery [J]. SPE Reservoir Evaluation and Engineering, 2005, 8(4): 291–299. doi: 10.2118/87230-PA [16] MITRA A, HARPALANI S, LIU S M. Laboratory measurement and modeling of coal permeability with continued methane production: Part 1-laboratory results [J]. Fuel, 2012, 94(4): 110–116. -