煤炭是我国的基础能源，也是重要的工业原料，到2030年，煤炭科学开采预测产能总量达30～35亿t[1]。所以，未来一段时间内煤炭生产仍将持续进行，煤矿的安全生产仍需加强保障。煤炭形成过程中伴生有大量的瓦斯气体，瓦斯既是瓦斯爆炸、煤与瓦斯突出灾害的诱因，又是一种高效的清洁能源[2-4]。掌握煤层中瓦斯流动规律是减少瓦斯事故、扩大瓦斯采集利用的前提，而煤体的渗透率是决定瓦斯流动规律的关键参数[5]，所以对煤体渗透率进行研究是十分必要的。我国大部分矿区的煤层渗透率较低，需通过采动卸压、水力压裂、深孔爆破等手段进行人工增透[6-8]，此外煤与瓦斯突出事故多发于掘进工作面[9-11]。因此，建立工作面前方采动应力下的煤体渗透率模型对提高瓦斯抽采和减少瓦斯事故具有重要意义。

以往的煤体渗透率模型较多以煤层气开采为背景[12-14]，同时基于单轴应变条件和上覆载荷不变的假设建立，而在煤矿采掘活动中，煤层的受力特点有别于煤层气开采。工作面推进过程中采煤工作面前方煤体经受采动应力路径的变化过程，其竖直方向的支承压力逐渐升高后下降、水平方向的应力持续降低，煤层受力呈现典型的三向应力变化特征(图1(a))。随着支承压力的升高以及人工增透活动的进行，煤体内部继发出现损伤和破裂现象，煤层内裂隙扩展贯通造成瓦斯压力下降，同时伴随瓦斯解吸扩散的发生，最终煤体渗透率急剧增长。采动应力决定了煤体在由原岩应力状态渐近为侧面临空过程中的受力状态、物理力学特征、温度梯度及渗透率变化有别于煤层气的抽采过程，以往的煤体渗透率模型并不能完全适用于采动应力条件。

笔者首先概述了煤体渗透率经典模型，然后结合采动过程中煤体内的力学变化机制及渗透率的控制因素提出了采动应力下煤体渗透模型构建中的3个关键问题，即采动煤体各向异性特征、采动煤体损伤破裂特征和煤体吸附解吸特征，并将这3个方面的国内外研究成果进行了总结和分析。在此基础上，将各向异性特征的煤体渗透率模型归纳为由有效应力变化表征和几何参数变化表征的两类，又将损伤破裂特征的煤体渗透率模型划分为本构方程中含损伤变量和渗透率表达式中含损伤变量的两类。最后，对采动应力下煤体渗透率模型的研究进展进行了介绍并对当前采动应力下煤体渗透率模型研究的不足和未来的发展方向进行了探讨。

1 煤体渗透率经典模型概述

为了对煤体的渗透率演化过程进行定量描述，国内外学者建立了大量的渗透率模型，这些研究成果主要集中在煤层气开采方面，其中一些模型具有较高的代表性，这些经典的模型成果在煤体渗透率模型的研究领域中发挥了重要作用。

GRAY[15]认为煤储层的渗透率与垂直于割理结构的有效水平应力大小有直接关系，较早建立了煤体渗透率演化模型。SEIDLE等[16]将储层煤体简化为“捆绑的火柴棍”几何结构，得到了煤体渗透率与有效静水压力的关系，从理论角度推导出渗透率变化的指数函数形式。PALMER和MANSOORI[17]假设储层为单轴应变条件，考虑煤层割理压缩性的变化，基于渗透率k与孔隙率φ的立方关系，结合非等温多孔弹性理论建立了煤层渗透率模型。SHI和DURUCAN[18]从各向同性的弹性本构方程出发，将煤储层视为连续介质，基于单轴应变的假设条件求出煤储层水平方向有效应力的变化，建立了指数函数形式的渗透率模型。CUI和BUSTIN[19]考虑整体平均有效应力对煤体孔隙度和渗透率的影响，利用多孔弹性本构方程进行推导得到解吸过程中煤体渗透率模型。ZHANG等[20]将煤体视为固体颗粒和孔隙空间组成，假设煤的孔隙空间同样存在吸附应变，建立了一个广义孔隙度的渗透率模型。LIU等[21]认为煤体渗透率由基质渗透率和裂隙渗透率2部分组成，考虑基质和裂隙的相互作用得到合成的渗透率表达式。经典的煤体渗透率模型汇总见表1。

表1 煤体渗透率经典模型
Table 1 Classical permeability models of coals

2 采动应力下煤体渗透率模型构建的关键问题

2.1 采动煤体各向异性特征

煤属于沉积岩，在沉积方向具有明显的层理结构，而在垂直层理的方向又存在空间狭窄的裂隙结构，这种结构被称为割理结构。割理可分为面割理和端割理结构，通常将两种割理结构简化为相互正交的形式。因此，煤体具有显著的各向异性特征(图1(b))，进而煤体在各个方向上存在不同的渗透率变化趋势。除此之外，采动应力下煤体渗透率模型的构建必须把握三向受力的特点，建立可表征竖向加卸载、水平卸载过程的渗透率理论模型。因此，无论是煤体内部结构特征还是外部受力条件均要求使用各向异性的渗透率模型对采动应力下煤体的渗透率演化过程进行描述。笔者根据渗透率模型的建模角度将现有的各向异性煤体渗透率模型成果归纳为有效应力变化和几何参数变化两类进行介绍。各向异性特征的煤体渗透率模型汇总见表2。

表2 各向异性渗透率模型
Table 2 Anisotropic permeability models

2.1.1 有效应力变化表征的煤体渗透率模型

有效应力变化表征的煤体渗透率模型都以煤体所受有效应力为自变量，在建模过程中较多将煤基质和裂隙视为一个整体，利用弹性本构方程对煤体在坐标轴3个方向的应力和应变关系进行推导，建立渗透率和有效应力之间的关系。此类模型[22-25]均可以归纳为以下形式

(1)

式中，ki0为每个方向的初始渗透率;为不同方向的割理压缩系数;为不同方向的有效应力;下标0为初始状态。

式(1)为归纳出的以有效应力为自变量的渗透率模型概化函数形式，不同学者以有效应力变化角度所建立的模型基本都采用这个函数形式，但相互之间的不同之处在于割理压缩系数和有效应力是否也为各向异性，即分别使用沿坐标轴3个方向的分量进行表达。CHEN等[22]建立了割理压缩系数和有效应力同时为各向异性的渗透率模型;而LIU等[23]仅考虑各方向割理压缩系数的异性特征，但各个方向的渗透率均由平均有效应力控制。除此以外，部分学者将不同方向的割理压缩系数视为恒定常数，建立了各方向有效应力控制各方向渗透率变化的模型。ZHANG等[24]认为模型表达式中的系数不是恒定常数3，根据煤体内各方向裂隙的分布特征，将系数修正为αi;PAN和CONNELL[25]基于单轴应变和上覆载荷不变的假设，并认为煤层符合横观各向同性的特征推导了水平方向的有效应力增量表达式:

(2)

式中，为i方向有效应力增量;Ei为i方向的煤体弹性模量;νij为i平面上j方向的泊松比;Δεsx为x方向的煤体吸附应变增量。

2.1.2 几何参数变化表征的煤体渗透率模型

几何参数变化表征的煤体渗透率模型大多以煤体内孔裂隙尺寸、裂隙应变和孔隙度为自变量，建模过程基本是将煤体划分为基质和裂隙两部分，通过计算应力、吸附作用下裂隙尺寸或孔隙度变化，建立煤体裂隙开度、孔隙度和渗透率的关系。为了考虑煤体各向异性的结构特征(图1(b))，这些渗透率模型[26-34]构建过程中大多将煤体抽象为立方体几何结构模型(图1(c))。

一些学者将气体在煤层裂隙中的流动简化为平板间流体流动过程，通过POISEUILLE定律得到渗透率与裂隙开度之间的立方关系(式(3))，并以此开始推导,即

(3)

式中，a，b分别为基质单元和裂隙单元的宽度。

GU和CHALATURNYK[26]将煤体简化为火柴棍模型，而马天然等[27]将煤体简化为立方体模型，分别对式(4)进行变换直接得到了渗透率和裂隙开度的3次方关系式。LIU等[28]将煤体抽象为立方体几何结构，认为煤体为各向同性的线弹性材料，提出了模量损失率来表征各方向弹性模量的变化过程，建立了不同边界条件下的各向异性渗透率模型。YANG等[29-31]认为工程分析中应忽略基质渗透率的影响，借此对LIU等[28]模型加以修正，构建出横观各向同性的渗透率模型。WANG等[32]从式(4)出发对压缩过程中裂隙孔径的变化进行推导，得出一个与LIU等[28]模型形式类似的表达式，同时定义了基质与整体单元应变比的新参数，建立了一个2维各向异性的渗透率模型。

另外，一些学者定量表示出三向有效应力变化过程中煤体内部孔隙度的变化，并将孔隙度变化结果代入渗透率和孔隙度的立方关系[17](式(4))建立了相应的各向异性渗透率模型。WU等[33]将煤体的物理力学性质视为各向同性，假设坐标轴3个方向下热膨胀和吸附作用产生的煤体应变相等而有效应力不相等，建立了双孔隙度的煤体各向异性渗透率模型。WANG等[34]基于立方体模型提出了内膨胀应变系数以对吸附过程中煤体基质的内膨胀量进行修正，建立了一个以三向有效应力和气体压力为自变量的渗透率模型。

(4)

式中，k为渗透率;φ为孔隙率。

从建模过程上来说，有效应力变化和几何参数变化的两个建模角度可认为是互逆的相反过程。有效应力的角度是将煤体所受有效应力通过孔隙度、应变、模量、泊松比、吸附常数等参数综合表示出来，而几何参数的角度是将裂隙开度或孔隙度用应力、模量、泊松比、吸附常数等参数进行表达，进而各自代入相应的框架基础公式。此外，根据各自建模所依据的框架基础，有效应力角度建模结果基本为指数型函数、几何参数角度建模的结果多数为3次方的幂函数。

渗透率模型中如果能将煤体的边界条件和结构参数均考虑为各向异性且持续变化的特征是最为理想和贴近实际的情况，然而受目前试验条件和测试手段所限，大量学者也只能将模型中的部分参数做出近似和假设处理。从上述各向异性的模型成果中可以看出，外部边界的三向应力和三向应变参数较易跟踪和获取，所以三向应力的特征在模型中容易实现[22-25]。结构异性参数大多较难准确把握，比如3个方向的割理压缩系数、弹性模量和泊松比等，因此一些模型中将结构参数视为各向同性[26-33]。

2.2 采动煤体损伤破裂特征

受开采扰动的影响，支承压力的空间分布具有分区化特征，煤体由原岩应力逐渐升高至应力峰值后迅速卸压(图1(a))，在此过程中煤体由原生初始状态渐进为损伤破裂状态，伴随煤体内微裂纹的扩展贯通，渗透率将会迅速增大，可达初始值的数十倍，乃至百倍[35-36]。因此，对于采场前方煤体必须考虑开采扰动引起的损伤破裂对渗透率的影响。

目前，国内外学者已经开始探究损伤破裂影响下渗透率的变化规律。开采扰动下煤体中瓦斯流动及渗透率的变化是一个典型的气固耦合问题。一种解决此问题的手段就是将应力场和渗流场的控制方程进行耦合求解。根据损伤力学理论可以定义多种形式的损伤变量，损伤变量对渗透率的控制作用通过耦合方程的形式进行传递。这种调控方式下，损伤变量处于固体力学的本构方程之中，而渗透率表达式内并不含损伤变量。另一种处理此问题的方式就是将定义出的损伤变量直接嵌入渗透率模型的表达式之中。笔者依据损伤变量所处的位置不同将损伤特征的实现方式归纳为两种:一种为损伤变量存在于本构方程中，另一种为损伤变量存在于渗透率表达式中。考虑损伤破裂特征的渗透率模型成果汇总见表3。

表3 考虑损伤破裂特征的渗透率模型
Table 3 Damage-based permeability models

2.2.1 本构方程中含损伤变量的煤体渗透率模型

一些学者基于宏观唯象的连续损伤力学理论建立了煤体损伤本构并将其代入固体力学控制方程，与瓦斯渗流方程联立进行求解。基于LEMAITRE[37]提出的应变等效假设(式(5))，魏明尧等[38]建立了弹塑性损伤本构;HU等[39]定义了双重孔隙结构下煤体的等效弹性损伤柔度张量，构建了双重孔隙介质损伤的本构模型;WANG等[40-41]将HU[39]定义的等效弹性损伤柔度张量进行简化后应用在单一孔隙率的渗流方程之中,即

(5)

式中，E0为无损材料的弹性模量;为受损材料的弹性模量;D为损伤变量。

另外，一部分学者通过细观损伤力学的方法建立了煤体的损伤本构，将煤体视为大量代表性体积单元(RVE)组成，通过研究每个RVE的力学破坏特征来反映煤体内损伤演化过程。ZHU和TANG[42-43]考虑岩石的非均质性，假设岩石参数符合Weibull分布(式(6))，将最大拉应力和Mohr-Coulomb准则作为判断煤体RVE破坏的准则(式(7))，建立了煤体的弹性损伤本构。沿用这种研究思路，一些学者建立了不同的气固耦合模型用以探究煤层中损伤-应力-渗流的问题，间接搭建了损伤变量和渗透率之间的关系[44-52]。除了根据渗透率与孔隙度立方关系(式(4))推导渗透率表达式的常见方法[44-45]之外，YANG等[48]利用分段函数依次建立了煤体单元在压缩和拉伸条件下的渗透率函数式;刘力源等[51-52]将煤体视为双重孔隙介质，结合煤体吸附CO2出现的损伤演化过程分别得到煤基质和裂隙的渗透率计算式。

(6)

(7)

式中，C1，C2为两个损伤破裂阈值函数;u为单元参数;u0为单元参数的平均值;m为函数的形状参数;σ1,σ3分别为单元的最大主应力和最小主应力;ft0为单轴抗拉强度;fc0为单轴抗压强度;φ′为内摩擦角。

2.2.2 渗透率表达式中含损伤变量的煤体渗透率模型

通过理论推导构建损伤变量并将其嵌入渗透率表达式的研究主要限于连续损伤力学的框架内进行，大多数学者基于应变等效假设展开研究，而损伤变量的构建方法多样。XUE等[53]通过孔隙度和有效应力定义了损伤变量;ZHANG等[54]采用煤体的体积应变、吸附应变和初始孔隙率构建了损伤变量;LI等[55]以轴向应变构建了损伤变量。另外，也有学者采用回归分析的方法，通过拟合系数来确定渗透率和不同损伤程度下应力状态之间的关系。CHEN等[56]为了描述煤体从弹性状态到峰后破坏的渗透率变化全过程，提出了修正系数L，通过回归函数得到有效偏应力与修正系数L的“S”形关系。ZHANG等[57]基于达西定律和裂隙的本构方程，通过拟牛顿算法的回归分析确定了以应力表示的渗透率函数。为了搭建渗透率和损伤变量之间的“桥梁”，另外一些学者立足于试验结果反映出的变化规律，在前人建立的渗透率模型基础之上进行改造，直接添加反映损伤影响的作用项或修正系数对渗透率变化过程进行推演，也取得了一些研究成果。ZHU等[58]在渗透率和孔隙度的立方关系(式(4))上进行改造，添加了指数形式的损伤作用项;ZHENG等[59]以CUI和BUSTIN[19]模型为基础，在该模型中添加了损伤作用项，假设了渗透率和损伤变量的指数关系;XUE等[60]通过室内试验的渗透率和声发射数据反分析提出在峰后阶段引入一个作用项来反映损伤影响，并假设峰后阶段渗透率和损伤变量呈线性关系;XIE等[61]假设损伤作用下渗透率的变化符合式(8)的关系,即

(8)

式中，ζ为与材料有关的系数。

可以看出，考虑损伤破裂特征的渗透率模型研究在近十年才逐渐展开，该领域的研究尚处于起始阶段，模型中损伤效应的研究深度尚浅，损伤变量的构建方式也较为分散。在模型适用性上，损伤变量存在本构方程中的实现方式可以适用较多边界条件情况，具有更广的适用范围，但是该方式下求解煤体渗透率需要依次计算应力场和渗流场控制方程，计算工作量较大，所以该方式更适合采用数值计算。损伤变量存在渗透率表达式的实现方式使得研究渗透率与各因素之间的数量关系更加便捷，但其大多基于已有模型成果进行改进，受原模型构建过程中的假设条件限制，诸如单轴应变条件、上覆载荷不变等。损伤变量所处不同位置的两种处理方式都可以达到在渗透率计算过程中考虑损伤破裂影响的目的。究其本质而言，损伤变量存在于本构方程的情况可视为损伤变量存在渗透率表达式情况的隐式函数形式。

2.3 煤体吸附解吸特征的表征方法

不同于其他气体在岩石中的渗透行为，瓦斯在煤层中运移时存在吸附解吸的特殊现象。煤层采掘过程中，煤体损伤破裂伴随着瓦斯渗流通道的形成。与此同时，开采扰动打破了煤层内原有的瓦斯压力平衡状态造成瓦斯压力持续下降，煤体内吸附态瓦斯释放引起解吸现象的发生。吸附解吸作用直接影响煤基质的胀缩变形，在渗透率模型构建的过程中，大多通过吸附应变这一中间变量间接反映煤体所受有效应力的大小或煤体内部渗流通道几何尺寸的变化[62]，因此关于吸附应变的计算方法就成为了研究吸附解吸影响下渗透率变化的焦点。

在早期的煤体渗透率模型中，通常将吸附应变与气体压力的增量视为线性关系进行处理。GRAY[15]通过提出等效吸附压力的概念建立了吸附应变与气体压力之间的线性关系。LEVINE[63]利用基质收缩系数确定了基质吸附应变和气体压力增量之间的线性函数。国内外学者普遍将煤对瓦斯的吸附认为是物理吸附过程且为单分子层吸附。基于此，Langmuir等温吸附模型[64]可以较好的用于煤对瓦斯气体吸附量的计算(式(9))。HARPALANI和CHEN[65-66]认为在气体吸附过程中吸附应变和气体吸附量成正比，采用Langmuir理论公式计算瓦斯气体吸附量，建立了煤体吸附应变和瓦斯气体吸附量的线性关系(式(10))。鉴于Langmuir等温吸附模型可以较好的描述气体压力改变下煤对瓦斯的吸附气体量，LEVINE[63]对Langmuir模型作出进一步改进，采用一个和Langmuir吸附模型具有相同数学函数形式的方程(式(11))来描述煤体吸附应变。

(9)

εs=αsVs(10)

(11)

式中，Vs为吸附的瓦斯气体含量;VL为Langmuir体积常数;αs为煤的吸附膨胀系数;εL为Langmuir吸附体积应变常数;pL为Langmuir吸附压力常数。

煤体对瓦斯吸附解吸的内在机理可通过热力学理论进行探究，所以除了采用上述Langmuir形式的方程对吸附应变进行表征之外，一些学者通过吸附热力学方法建立了吸附应变的计算式[67-72]，详见表4。

表4 吸附热力学角度建立的吸附应变表达式
Table 4 Expressions of adsorption strain based on adsorption thermodynamics

采用以气体压力增量为自变量的线性函数对煤吸附应变进行描述是最早的研究方法，其优点是瓦斯气体压力的测定较为简单，但是线性函数却有夸大瓦斯压力对吸附应变影响的缺点。通过对比试验结果发现Langmuir等温吸附理论可以很好的描述煤对瓦斯的吸附解吸过程，上述Langmuir形式的吸附应变表达式具有一定的理论基础且方程形式较为简洁，表达式中的VL,εL,pL等参数也较易测试，所以式(10)[73-76]和式(11)[77-80]被后来学者大量效仿使用。通过热力学方法建立的吸附应变大多以吸附热力学为基础进行推导，由此得出的吸附应变表达式具有更为严密的理论基础，然而此类方法的模型参数较多且部分参数难以准确测定，所以此类吸附应变表达式的实用性有所限制。

3 采动应力下煤体渗透率模型研究进展

现有的煤体渗透率模型中同时考虑各向异性、损伤破裂、吸附解吸等特征的渗透率模型还较少。近些年，部分学者已经开始关注这一问题，并建立了相应的采动应力下煤体渗透率模型。笔者将这些模型成果归纳为有效应力型、几何参数型和系数拟合型的3种进行介绍。

3.1 有效应力型的煤体渗透率模型

有效应力型的煤体渗透率模型主要以有效应力、气体压力、弹性模量和泊松比等力学参数为基本变量对渗透率的变化进行表征。其中，ZHANG等[81]基于立方体模型结构(图1(c))，运用损伤和断裂力学理论提出了一个可考虑采动应力引起微裂纹扩展、吸附解吸造成煤体变形的有效应力型渗透率模型(式(12))，并根据谢和平等[82]提出的3种典型开采方式下采动应力路径(图2)，对采动煤体渗透率变化规律进行了拟合分析。

(12)

式中，Di为裂纹损伤变量;ci为体现弹性常数和裂纹展布系数的综合参数;φi0为不同方向的面孔隙率;为不同方向的有效应力增量;为不同方向的煤体弹性模量;为不同方向的基质弹性模量;为与裂纹周边煤基质硬度有关的弹性模量;为沿设定方向的煤体泊松比;为沿设定方向的基质泊松比;εLi,pLi为不同方向的Langmuir吸附常数;Ti为局部拉应力;FLi为不同方向的吸附影响系数;αi为幂函数系数;下标t为j,k方向的可轮换的坐标轴的方向。

图2(a)～(c)为不同开采方式下煤体渗透率的试验和拟合结果。可以看出在3种开采方式下，随着采动卸压过程的进行，煤体内微裂隙不断萌生和扩展，损伤积累不断增大，渗透率的实测和拟合结果持续升高。无煤柱和放顶煤开采方式的渗透率在加卸载途中出现陡然上升的情况，但拟合结果却没有出现陡然升高的特征。不过整体来看，式(12)基本可以反映出煤体渗透率随着采动卸压过程的进行而持续升高的变化规律。

3.2 几何参数型的煤体渗透率模型

几何参数型的煤体渗透率模型主要以孔隙率和煤体应变为基本参数来构建渗透率模型。ZHENG等[83]将LIU等[28]提出的煤体各向异性渗透率模型和ZHU等[58]定义的损伤破裂特征修正项进行联结。通过细观损伤力学的方法构建损伤变量以对渗透率进行修正，合成出可表征采动应力下煤体损伤的几何参数型渗透率表达式:

(13)

式中，εj为j方向的应变;εt0,εc0分别为最大拉伸主应变和最大压缩主应变;ε1,ε3分别为最大主应变和最小主应变。

其后，利用式(13)分别对考虑和不考虑损伤效应的保护层底板下煤岩体渗透率进行了模拟分析。由图3可以看出距离保护层底板下方80 m处的岩体渗透率为初始状态;在与底板距离缩进过程中煤岩体渗透率逐渐升高;在距离底板12 m之外，考虑损伤与否的渗透率结果相同;而在12 m范围之内，考虑损伤效应的煤岩体渗透率比率出现急剧上升，从初始渗透率比率的180倍陡增至680倍左右，而不考虑损伤的煤岩体渗透率比率仅升高至200倍左右。

3.3 系数拟合型的煤体渗透率模型

系数拟合型的煤体渗透率模型主要通过试验数据拟合的手段来确定一系列待定系数的变化与不同区段渗透率的关系。XUE等[84]根据煤体三轴压缩试验过程中应力-应变-渗透率关系的变化特征(图4)，将全应力-应变过程划分为3阶段:弹性阶段、塑性软化阶段、残余状态阶段。同时认为弹性阶段没有微裂纹产生，塑性软化阶段萌生大量裂纹，进而结合巷道开挖后围岩损伤的3区特征，将围岩的渐进损伤区与全应力-应变过程3阶段进行了一一对应:弹性阶段对应无损区、塑形软化阶段对应开挖损伤区、残余状态阶段对应强开挖损伤区。由此，得到了巷道开挖后，掘进工作面周围煤体的系数拟合型渗透率表达式

(14)

式中，αki,βki,γki(i=1,2,3)分别为不同应力和应变条件下描述渗透率变化的系数;εip为煤进入塑性状态的应变;εir为煤进入残余状态的应变。

图5为掘进工作面推进过程中周围煤体内渗透率分布演化结果，周围煤体的渗透率分布同样符合三区特征。距离巷道中心较远处，煤体渗透率处于初始状态;随着与巷道中心距离的缩减，无损区内的煤体渗透率由于初始裂隙的闭合而渗透率逐渐减小(图5 Ⅰ)，距离至塑性区的边界时，渗透率达到最小值;开挖损伤区内由于产生大量裂纹而造成渗透率持续升高(图5 Ⅱ);强开挖损伤区内煤体破裂程度较为严重，由此出现大量贯通裂隙引起渗透率急剧升高(图5 Ⅲ)。对比不同开挖时间下的渗透率分布结果可知，随着时间的推移，煤体内渗透率的变化幅度呈上升趋势。

开采扰动下工作面前方煤体或底板下煤岩体的渗透率均呈现分区段变化的特征，上述3个渗透率模型分别从断裂力学、细观损伤力学和连续损伤力学的角度考虑了采动过程中煤体损伤破裂造成的渗透率分区特征。在此之中，有效应力型的模型运用断裂力学理论，确定了微裂纹的开裂判据，通过损伤理论构建了微裂纹和损伤变量的关系，建立的模型较好的体现了采动应力和损伤破裂特征，模型的理论推导过程较为严密，具有较好的基础理论支撑。几何参数型的模型加入了煤体参数的微观统计分布，更加贴近煤体内部孔裂隙结构导致的非均匀性特征，同时模型采用拉伸和剪切破坏准则作为单元失效判据更为全面，模型从细观损伤本构出发考虑损伤破裂的因素对煤体变形和渗透的影响，适用范围较广，但计算中将煤体划分为细观单元后，求解计算的工作量较大，所以较为适合进行数值计算。系数拟合型的模型从形式上来看较为简洁，可显著区分工作面前方煤体的损伤分布特征，可以较好的反映各损伤区煤体渗透率的变化规律，此外该模型也易于进行理论计算，尤其是便于解析计算。

4 结论与展望

(1)相比天然气和煤层气的开采过程，采动应力环境下的煤层在受力方式和结构形式上更为复杂。工作面推进过程中采动应力的3个应力主方向出现不同的变化趋势，采动应力可视为外部边界条件，而煤层由于沉积作用出现的层理和割理结构属于内部结构条件。煤炭开采过程中，其内部结构和外部边界条件的主方向并非一致，所以需要考虑外部边界条件和内部渗透率主方向的夹角问题，统一两者之间的关系。与此同时，割理结构的各向异性也造成煤体物理力学参数存在异性特征。因此，需要结合地质条件和各向异性的物理力学参数建立一个考虑采动应力、煤层产状和渗透率方向的模型，以此确定煤层瓦斯渗透的最优化方向，指导煤层增透和瓦斯抽采工作。

(2)采动卸压是一个打破原始地应力后应力持续重分布引起煤层损伤破裂的增透过程。现有的煤体渗透率模型已经开始逐渐关注到这个问题，然而在已有成果中完全基于损伤理论推导出的损伤渗透率模型还较少，部分经验模型的实用性较好，但推导过程存在较多前提假设。因此，损伤渗透率模型的理论基础仍需进一步夯实。除此之外，结合断裂力学理论对煤体内裂隙演化过程进行定量描述，探究裂隙连通度对煤体渗透性的影响，确定采场尺度下不同裂隙网络空间分布与煤层渗透性的关系[85-88]，建立一个多理论交叉、宏细观结合的渗透率模型是十分必要的。

(3)煤体吸附应变的研究中，大量的实验和理论成果都是基于等温条件获得，而采面前方煤体内的气体压力和温度由表及内却是不断升高的。此外，在瓦斯吸附和解吸过程中，煤体应变严格来讲并不具有完全可逆的时空演化特征[89-91]。因此，加强非等温条件下瓦斯吸附特性的研究，掌握瓦斯解吸和吸附过程中煤体应变的差异性对采动影响下煤层渗透率演化过程的精确描述大有裨益。

总体来说，多因素影响下的采动煤体渗透率模型的理论体系并未完全建立，能够较为系统和完善体现上述特征的采动煤体渗透率模型也还鲜有出现，但煤体渗透率模型在各向异性特征、损伤破裂特征和吸附解吸特征的方面均已取得大量研究成果，以这些成果为基础的采动煤体渗透率模型将会迅速发展。同时，以上各方面的研究深度及广度仍需进一步加强。除此之外，采动应力下煤体渗透率模型的研究还有一些问题需要面对。如今，煤炭开采深度的迅速增大给煤体渗透率的研究带来巨大挑战。一方面，深部煤体的赋存环境具有“三高”特征，即高温、高地应力和高岩溶水压[92-94]。深部开采是多场多相环境共同影响的工程活动，采掘和增透过程中“三高”因素的变化相互耦合，任一因素的变化都将严重影响煤层的透气性效果。因此，深部采动煤体渗透率模型需要解决多场多相耦合的问题。另一方面，深部煤体开采过程中具有显著的强扰动和强时效特征[95-96]。强扰动是指采掘活动对深部高应力区岩体形成的强烈卸荷效应;强时效是指深部环境下采动岩体具有与采掘活动无关的、明显的流变效应。显而易见，强扰动和强时效特征将会强烈影响深部煤体的渗透特性。所以，为了全面反映深部煤体渗透率的控制因素，深部采动煤体渗透率模型还需要考虑强扰动和强时效特征的影响。

参考文献

[1] 谢和平,钱鸣高,彭苏萍,等.煤炭科学产能及发展战略初探[J].中国工程科学,2011,13(6):44-50.

XIE Heping,QIAN Minggao,PENG Suping,et al.Sustainable capacity of coal mining and its strategic plan[J].Engineering Sciences,2011,13(6):44-50.

[2] 林柏泉.瓦斯爆炸动力学特征参数的测定及其分析[J].煤炭学报,2002,27(2):164-167.

LIN Baiquan.The measurement and analysis of dynamics feature parameter in gas explosion[J].Journal of China Coal Society,2002,27(2):164-167.

[3] 聂百胜,何学秋,王恩元,等.煤与瓦斯突出预测技术研究现状及发展趋势[J].中国安全科学学报,2003,13(6):43-46,83.

NIE Baisheng,HE Xueqiu,WANG Enyuan,et al.Present situation and progress trend of prediction technology of coal and gas outburst[J].China Safety Science Journal,2003,13(6):43-46,83.

[4] 谢和平,鞠杨,高明忠,等.煤炭深部原位流态化开采的理论与技术体系[J].煤炭学报,2018,43(5):1210-1219.

XIE Heiping,JU Yang,GAO Mingzhong,et al.Theories and technologies for in-situ fluidized mining of deep underground coal resources[J].Journal of China Coal Society,2018,43(5):1210-1219.

[5] 周世宁.瓦斯在煤层中流动的机理[J].煤炭学报,1990,15(1):15-24.

ZHOU Shining.Mechanism of gas flow in coal seams[J].Journal of China Coal Society,1990,15(1):15-24.

[6] 谢和平,周宏伟,薛东杰,等.我国煤与瓦斯共采:理论、技术与工程[J].煤炭学报,2014,39(8):1391-1397.

XIE Heping,ZHOU Hongwei,XUE Dongjie,et al.Theory,technology and engineering of simultaneous exploitation of coal and gas in China[J].Journal of China Coal Society,2014,39(8):1391-1397.

[7] 袁亮.我国深部煤与瓦斯共采战略思考[J].煤炭学报,2016,41(1):1-6.

YUAN Liang.Strategic thinking of simultaneous exploitation of coal and gas in deep mining[J].Journal of China Coal Society,2016,41(1):1-6.

[8] 谢和平,高峰,周宏伟,等.煤与瓦斯共采中煤层增透率理论与模型研究[J].煤炭学报,2013,38(7):1101-1108.

XIE Heping,GAO Feng,ZHOU Hongwei,et al.On theoretical and modeling approach to mining enhanced permeability for simultaneous exploitation of coal and gas[J].Journal of China Coal Society,2013,38(7):1101-1108.

[9] 姜福兴,尹永明,朱权洁,等.基于掘进面应力和瓦斯浓度动态变化的煤与瓦斯突出预警试验研究[J].岩石力学与工程学报,2014,33(S2):3581-3588.

JIANG Fuxing,YIN Yongming,ZHU Quanjie,et al.Experimental study of precaution technology of heading face coal and gas outburst based on dynamic changes of stress and methane concentration[J].Chinese Journal of Rock Mechanics and Engineering,2014,33(S2):3581-3588.

[10] 马念杰,赵希栋,赵志强,等.掘进巷道蝶型煤与瓦斯突出机理猜想[J].矿业科学学报,2017,2(2):137-149.

MA Nianjie,ZHAO Xidong,ZHAO Zhiqiang,et al.Conjecture about mechanism of butterfly-shape coal and gas outburst in excavation roadway[J].Journal of Mining Science and Technology,2017,2(2):137-149.

[11] 舒龙勇,王凯,齐庆新,等.煤巷掘进面应力场演化特征及突出危险性评价模型[J].采矿与安全工程学报,2017,34(2):259-267.

SHU Longyong,WANG Kai,QI Qingxin,et al.Stress field evolution characteristics and coal-gas outburst hazard evaluation model of the heading face in coal roadway[J].Journal of Mining & Safety Engineering,2017,34(2):259-267.

[12] PALMER I.Permeability changes in coal:Analytical modeling[J].International Journal of Coal Geology,2009,77(1-2):119-126.

[13] WANG G G X,ZHANG X D,WEI X R,et al.A review on transport of coal seam gas and its impact on coalbed methane recovery[J].Frontiers of Chemical Science and Engineering,2011,5(2):139-161.

[14] PAN Z J,CONNELL L D.Modelling permeability for coal reservoirs:A review of analytical models and testing data[J].International Journal of Coal Geology,2012,92:1-44.

[15] GRAY I.Reservoir engineering in coal seams:Part 1 the physical process of gas storage and movement in coal seams[J].SPE Reservoir Engineering,1987,2(1):28-34.

[16] SEIDLE J P,JEANSONNE M W,ERICKSON D J.Application of matchstick geometry to stress dependent permeability in coals[A].SPE Rocky Mountain Regional Meeting[C].Society of Petroleum Engineers,1992:433-444.

[17] PALMER I,MANSOORI J.How permeability depends on stress and pore pressure in coalbeds:A new model[J].SPE Reservoir Evaluation & Engineering,1998,1(6):539-544.

[18] 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.

[19] CUI X J,BUSTIN R M.Volumetric strain associated with methane desorption and its impact on coalbed gas production from deep coal seams[J].American Association of Petroleum Geologists Bulletin,2005,89(9):1181-1202.

[20] ZHANG H B,LIU J S,ELSWORTH D.How sorption-induced matrix deformation affects gas flow in coal seams:A new FE model[J].International Journal of Rock Mechanics and Mining Sciences,2008,45(8):1226-1236.

[21] LIU J S,CHEN Z W,ELSWORTH D,et al.Evaluation of stress-controlled coal swelling processes[J].International Journal of Coal Geology,2010,83(4):446-455.

[22] CHEN D,PAN Z J,LIU J S,et al.Characteristic of anisotropic coal permeability and its impact on optimal design of multi-lateral well for coalbed methane production[J].Journal of Petroleum Science and Engineering,2012,88-89:13-28.

[23] LIU Y B,LI M H,YIN G Z,et al.Permeability evolution of anthracite coal considering true triaxial stress conditions and structural anisotropy[J].Journal of Natural Gas Science and Engineering,2018,52:492-506.

[24] ZHANG Z T,ZHANG R,XIE H P,et al.The relationships among stress,effective porosity and permeability of coal considering the distribution of natural fractures:Theoretical and experimental analyses[J].Environmental Earth Sciences,2015,73(10):5997-6007.

[25] PAN Z J,CONNELL L D.Modelling of anisotropic coal swelling and its impact on permeability behaviour for primary and enhanced coalbed methane recovery[J].International Journal of Coal Geology,2011,85(3-4):257-267.

[26] GU F G,CHALATURNYK R.Permeability and porosity models considering anisotropy and discontinuity of coalbeds and application in coupled simulation[J].Journal of Petroleum Science and Engineering,2010,74(3-4):113-131.

[27] 马天然,刘卫群,JONNY Rutqvist,等.裂隙煤岩各向异性渗透率模型和煤层注气THM耦合行为[J].煤炭学报,2017,42(S2):407-416.

MA Tianran,LIU Weiqun,JONNY Rutqvist,et al.Anisotropy permeability model for highly fractured coal seams associated with coupled THM analysis of CO2-ECBM[J].Journal of China Coal Society,2017,42(S2):407-416.

[28] LIU J S,CHEN Z W,ELSWORTH D,et al.Linking gas-sorption induced changes in coal permeability to directional strains through a modulus reduction ratio[J].International Journal of Coal Geology,2010,83(1):21-30.

[29] YANG D S,QI X Y,CHEN W Z,et al.Numerical investigation on the coupled gas-solid behavior of coal using an improved anisotropic permeability model[J].Journal of Natural Gas Science and Engineering,2016,34:226-235.

[30] 亓宪寅,王威.基于结构异性比的含瓦斯煤渗透各向异性研究[J].岩土工程学报,2017,39(6):1030-1037.

QI Xianyin,WANG Wei.Anisotropic permeability model for coal containing methane based on anisotropic structure ratio[J].Chinese Journal of Geotechnical Engineering,2017,39(6):1030-1037.

[31] 亓宪寅,李家卓,王威.基于不同方向模量损失率的含瓦斯煤各向异性渗透模型[J].岩土力学,2018,39(2):635-643.

QI Xianyin,LI Jiazhuo,WANG Wei.An anisotropic permeability model of coal containing methane based on different directional modulus reduction ratios[J].Rock and Soil Mechanics,2018,39(2):635-643.

[32] WANG J G,LIU J S,KABIR A.Combined effects of directional compaction,non-Darcy flow and anisotropic swelling on coal seam gas extraction[J].International Journal of Coal Geology,2013,109-110:1-14.

[33] 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.

[34] WANG K,ZANG J,WANG G D,et al.Anisotropic permeability evolution of coal with effective stress variation and gas sorption:Model development and analysis[J].International Journal of Coal Geology,2014,130:53-65.

[35] ZHANG H,CHENG Y P,LIU Q Q,et al.A novel in-seam borehole hydraulic flushing gas extraction technology in the heading face:Enhanced permeability mechanism,gas flow characteristics,and application[J].Journal of Natural Gas Science and Engineering,2017,46:498-514.

[36] HUANG Q M,WU B,CHENG W M,et al.Investigation of permeability evolution in the lower slice during thick seam slicing mining and gas drainage:A case study from the Dahuangshan coalmine in China[J].Journal of Natural Gas Science and Engineering,2018,52:141-154.

[37] LEMAITRE J.A course on damage mechanics[M].Berlin:SpringerVerlag,1992.

[38] 魏明尧,王春光,崔光磊,等.损伤和剪胀效应对裂隙煤体渗透率演化规律的影响研究[J].岩土力学,2016,37(2):574-582.

WEI Mingyao,WANG Chunguang,CUI Guanglei,et al.Influences of damage and shear dilation on permeability evolution of fractured coal[J].Rock and Soil Mechanics,2016,37(2):574-582.

[39] HU S B,WANG E Y,KONG X G.Damage and deformation control equation for gas-bearing coal and its numerical calculation method[J].Journal of Natural Gas Science and Engineering,2015,25:166-179.

[40] WANG E Y,KONG X G,HU S B,et al.Multi-scale fractured coal gas-solid coupling model and its applications in engineering projects[J].Transport in Porous Media,2018,121(3):703-724.

[41] KONG X G,WANG E Y,LIU Q L,et al.Dynamic permeability and porosity evolution of coal seam rich in CBM based on the flow-solid coupling theory[J].Journal of Natural Gas Science and Engineering,2017,40:61-71.

[42] ZHU W C,TANG C A.Micromechanical model for simulating the fracture process of rock[J].Rock Mechanics and Rock Engineering,2004,37(1):25-56.

[43] ZHU W C,WEI C H,LIU J,et al.Impact of gas adsorption induced coal matrix damage on the evolution of coal permeability[J].Rock Mechanics and Rock Engineering,2013,46(6):1353-1366.

[44] 孟磊.含瓦斯煤体损伤破坏特征及瓦斯运移规律研究[D].北京:中国矿业大学(北京),2013.

MENG Lei.Research on coal damage characteristics and gas migration law of coal containing gas[D].Beijing:China University of Mining and Technology(Beijing),2013.

[45] 薛熠.采动影响下损伤破裂煤岩体渗透性演化规律研究[D].徐州:中国矿业大学,2017.

XUE Yi.Study on the permeability evolution of fractured coal under the influence of mining[D].Xuzhou:China University of Mining and Technology,2017.

[46] 徐涛,唐春安,宋力,等.含瓦斯煤岩破裂过程流固耦合数值模拟[J].岩石力学与工程学报,2005,24(10):1667-1673.

XU Tao,TANG Chunan,SONG Li,et al.Numerical simulation of coupled gas flow in failure process of gassy coal-rock[J].Chinese Journal of Rock Mechanics and Engineering,2005,24(10):1667-1673.

[47] XU T,TANG C A,YANG T H,et al.Numerical investigation of coal and gas outbursts in underground collieries[J].International Journal of Rock Mechanics and Mining Sciences,2006,43(6):905-919.

[48] YANG T H,XU T,LIU H Y,et al.Stress-damage-flow coupling model and its application to pressure relief coal bed methane in deep coal seam[J].International Journal of Coal Geology,2011,86(4):357-366.

[49] 张春会,于永江,赵全胜.非均匀煤岩渗流-应力弹塑性耦合数学模型及数值模拟[J].岩土力学,2009,30(9):2837-2842.

ZHANG Chunhui,YU Yongjiang,ZHAO Quansheng.Seepage-stress elastoplastic coupling model of heterogeneous coal and numerical simulation[J].Rock and Soil Mechanics,2009,30(9):2837-2842.

[50] LU Y L,ELSWORTH D,WANG L G.Microcrack-based coupled damage and flow modeling of fracturing evolution in permeable brittle rocks[J].Computers and Geotechnics,2013,49:226-244.

[51] 刘力源,朱万成,魏晨慧,等.气体吸附诱发煤强度劣化的力学模型与数值分析[J].岩土力学,2018,39(4):1500-1508.

LIU Liyuan,ZHU Wancheng,WEI Chenhui,et al.Mechanical model and numerical analysis of mechanical property alterations of coal induced by gas adsorption[J].Rock and Soil Mechanics,2018,39(4):1500-1508.

[52] ZHU W C,LIU L Y,LIU J S,et al.Impact of gas adsorption-induced coal damage on the evolution of coal permeability[J].International Journal of Rock Mechanics and Mining Sciences,2018,101:89-97.

[53] XUE Y,GAO F,GAO Y N,et al.Quantitative evaluation of stress-relief and permeability-increasing effects of overlying coal seams for coal mine methane drainage in Wulan coal mine[J].Journal of Natural Gas Science and Engineering,2016,32:122-137.

[54] ZHANG N,LI X R,CHENG H M,et al.A coupled damage-hydro-mechanical model for gas drainage in low-permeability coalbeds[J].Journal of Natural Gas Science and Engineering,2016,35:1032-1043.

[55] LI X C,NIE B S,LIU F B,et al.The relation of gas seepage and coal body damage under the true three dimension stress[J].Procedia Engineering,2011,26:1462-1466.

[56] CHEN D,PAN Z J,SHI J Q,et al.A novel approach for modelling coal permeability during transition from elastic to post-failure state using a modified logistic growth function[J].International Journal of Coal Geology,2016,163:132-139.

[57] ZHANG X M,ZHANG D M,LEO C J,et al.Damage evolution and post-peak gas permeability of raw coal under loading and unloading conditions[J].Transport in Porous Media,2017,117(3):465-480.

[58] ZHU W C,WEI C H,LI S,et al.Numerical modeling on destress blasting in coal seam for enhancing gas drainage[J].International Journal of Rock Mechanics and Mining Sciences,2013,59:179-190.

[59] ZHENG C S,MEHMET K,CHEN Z W,et al.Effects of coal damage on permeability and gas drainage performance[J].International Journal of Mining Science and Technology,2017,27(5):783-786.

[60] XUE Y,GAO F,LIU X G.Effect of damage evolution of coal on permeability variation and analysis of gas outburst hazard with coal mining[J].Natural Hazards,2015,79(2):999-1013.

[61] XIE J,GAO M Z,YU B,et al.Coal permeability model on the effect of gas extraction within effective influence zone[J].Geomechanics and Geophysics for Geo-Energy and Geo-Resources,2015,1(1-2):15-27.

[62] LIU J S,CHEN Z W,ELSWORTH D,et al.Interactions of multiple processes during CBM extraction:A critical review[J].International Journal of Coal Geology,2011,87(3-4):175-189.

[63] LEVINE J R.Model study of the influence of matrix shrinkage on absolute permeability of coal bed reservoirs[J].Geological Society London Special Publications,1996,109(1):197-212.

[64] Langmuir I.The adsorption of gases on plane surfaces of glass,mica and platinum.[J].Journal of Chemical Physics,2015,40(12):1361-1403.

[65] HARPALANI S,CHEN G L.Estimation of changes in fracture porosity of coal with gas emission[J].Fuel,1995,74(10):1491-1498.

[66] HARPALANI S,CHEN G L.Influence of gas production induced volumetric strain on permeability of coal[J].Geotechnical and Geological Engineering,1997,15(4):303-325.

[67] PAN Z J,CONNELL L D.A theoretical model for gas adsorption-induced coal swelling[J].International Journal of Coal Geology,2007,69(4):243-252.

[68] LIU S M,HARPALANI S.A new theoretical approach to model sorption-induced coal shrinkage or swelling[J].AAPG Bulletin,2013,97(7):1033-1049.

[69] LIU Z X,FENG Z C,ZHANG Q M,et al.Heat and deformation effects of coal during adsorption and desorption of carbon dioxide[J].Journal of Natural Gas Science and Engineering,2015,25:242-252.

[70] VERMOREL R,PIJAUDIER-CABOT G.Enhanced continuum poromechanics to account for adsorption induced swelling of saturated isotropic microporous materials[J].European Journal of Mechanics A:Solids,2014,44(1):148-156.

[71] HOL S,PEACH C J,SPIERS C J.Effect of 3-D stress state on adsorption of CO2 by coal[J].International Journal of Coal Geology,2012,93:1-15.

[72] LIU J F,PEACH C J,ZHOU H W,et al.Thermodynamic models for swelling of unconfined coal due to adsorption of mixed gases[J].Fuel,2015,157:151-161.

[73] SHI J Q,DURUCAN S.A model for changes in coalbed permeability during primary and enhanced methane recovery[J].SPE Reservoir Evaluation & Engineering,2005,8(4):291-299.

[74] CUI X J,BUSTIN R M,CHIKATAMARLA L.Adsorption-induced coal swelling and stress:Implications for methane production and acid gas sequestration into coal seams[J].Journal of Geophysical Research Solid Earth,2007,112(B10):1-8.

[75] CONNELL L D.Coupled flow and geomechanical processes during gas production from coal seams[J].International Journal of Coal Geology,2009,79(1-2):18-28.

[76] QU H Y,LIU J S,CHEN Z W,et al.Complex evolution of coal permeability during CO2 injection under variable temperatures[J].International Journal of Greenhouse Gas Control,2012,9:281-293.

[77] SHI J Q,DURUCAN S,SHIMADA S.How gas adsorption and swelling affects permeability of coal:A new modelling approach for analysing laboratory test data[J].International Journal of Coal Geology,2014,128-129:134-142.

[78] KUMAR H,ELSWORTH D,LIU J S,et al.Optimizing enhanced coalbed methane recovery for unhindered production and CO2 injectivity[J].International Journal of Greenhouse Gas Control,2012,11:86-97.

[79] WANG J G,LIU J S,KABIR A.Combined effects of directional compaction,non-Darcy flow and anisotropic swelling on coal seam gas extraction[J].International Journal of Coal Geology,2013,109-110:1-14.

[80] LIU Q Q,CHENG Y P,WANG H F,et al.Numerical assessment of the effect of equilibration time on coal permeability evolution characteristics[J].Fuel,2015,140:81-89.

[81] ZHANG Z T,ZHANG R,XIE H P,et al.An anisotropic coal permeability model that considers mining-induced stress evolution,microfracture propagation and gas sorption-desorption effects[J].Journal of Natural Gas Science and Engineering,2017,46:664-679.

[82] 谢和平,周宏伟,刘建锋,等.不同开采条件下采动力学行为研究[J].煤炭学报,2011,36(7):1067-1074.

XIE Heping,ZHOU Hongwei,LIU Jianfeng,et al.Mining-induced mechanical behavior in coal seams under different mining layouts[J].Journal of China Coal Society,2011,36(7):1067-1074.

[83] ZHENG C S,KIZIL M S,CHEN Z W,et al.Role of multi-seam interaction on gas drainage engineering design for mining safety and environmental benefits:linking coal damage to permeability variation[J].Process Safety and Environmental Protection,2018,114:310-322.

[84] XUE Y,GAO F,LIU X G,et al.Permeability and pressure distribution characteristics of the roadway surrounding rock in the damaged zone of an excavation[J].International Journal of Mining Science and Technology,2017,27(2):211-219.

[85] ZHOU H W,LIU J F,XUE D J,et al.Numerical simulation of gas flow process in mining-induced crack network[J].International Journal of Mining Science and Technology,2012,22(6):793-799.

[86] XUE D J,ZHOU H W,CHEN C F,et al.A combined method for evaluation and prediction on permeability in coal seams during enhanced methane recovery by pressure-relieved method[J].Environmental Earth Sciences,2015,73(10):5963-5974.

[87] ZHOU H W,ZHONG J C,REN W G,et al.Characterization of pore-fracture networks and their evolution at various measurement scales in coal samples using X-ray μCT and a fractal method[J].International Journal of Coal Geology,2018,189:35-49.

[88] ZHOU H W,YANG S.Fractional derivative approach to non-Darcian flow in porous media[J].Journal of Hydrology,2018,566:910-918.

[89] MAJEWSKA Z,MAJEWSKI S,ZIETEK J.Swelling of coal induced by cyclic sorption/desorption of gas:Experimental observations indicating changes in coal structure due to sorption of CO2 and CH4[J].International Journal of Coal Geology,2010,83(4):475-483.

[90] 祝捷,张敏,传李京,等.煤吸附/解吸瓦斯变形特征及孔隙性影响实验研究[J].岩石力学与工程学报,2016,35(S1):2620-2626.

ZHU Jie,ZHANG Min,CHUAN Lijing,et al.Experimental study on coal strain induced by methane sorption/desorption and effect of pore features[J].Chinese Journal of Rock Mechanics and Engineering,2016,35(S1):2620-2626.

[91] 李祥春,张良,赵建飞,等.瓦斯气体吸附解吸过程煤变形响应特征[J].矿业科学学报,2018,3(1):46-54.

LI Xiangchun,ZHANG Liang,ZHAO Jianfei,et al.Coal deformation characteristics in gas adsorption and desorption[J].Journal of Mining Science and Technology,2018,3(1):46-54.

[92] 何满潮,谢和平,彭苏萍,等.深部开采岩体力学研究[J].岩石力学与工程学报,2005,24(16):2803-2813.

HE Manchao,XIE Heping,PENG Suping,et al.Study on rock mechanics in deep mining engineering[J].Chinese Journal of Rock Mechanics and Engineering,2005,24(16):2803-2813.

[93] 谢和平,周宏伟,薛东杰,等.煤炭深部开采与极限开采深度的研究与思考[J].煤炭学报,2012,37(4):535-542.

XIE Heping,ZHOU Hongwei,XUE Dongjie,et al.Research and consideration on deep coal mining and critical mining depth[J].Journal of China Coal Society,2012,37(4):535-542.

[94] 周宏伟,谢和平,左建平.深部高地应力下岩石力学行为研究进展[J].力学进展,2005,35(1):91-99.

ZHOU Hongwei,XIE Heping,ZUO Jiangping.Developments in researches on mechanical behaviors of rocks under the condition of high ground pressure in the depths[J].Advances in Mechanics,2005,35(1):91-99.

[95] 谢和平.“深部岩体力学与开采理论”研究构想与预期成果展望[J].工程科学与技术,2017,49(2):1-16.

XIE Heping.Research framework and anticipated results of deep rock mechanics and mining theory[J].Advanced Engineering Sciences,2017,49(2):1-16.

[96] 薛东杰,周宏伟,彭瑞东,等.基于应力降的非连续支承压力强扰动特征研究[J].岩石力学与工程学报,2018,37(5):1080-1095.

XUE Dongjie,ZHOU Hongwei,PENG Ruidong,et al.Strong disturbance of discontinuous abutment pressure[J].Chinese Journal of Rock Mechanics and Engineering,2018,37(5):1080-1095.