深部煤体采动应力下双曲函数型渗透率模型

2019-04-11王路军周宏伟荣腾龙任伟光

煤炭学报 2019年3期

王路军，周宏伟，荣腾龙，任伟光

(1.中国矿业大学(北京) 力学与建筑工程学院，北京 100083; 2.中国矿业大学(北京) 能源与矿业学院，北京 100083; 3.中国矿业大学(北京) 煤炭资源与安全开采国家重点实验室，北京 100083)

随着煤炭资源开采深度的增加，深部煤岩体处于高地应力、高地温、高渗透压环境中，在开采扰动的条件下，岩爆、煤与瓦斯突出、巷道难以支护等现象发生概率随之增大。目前，深部矿产资源开采趋于常态化，深部强扰动和强时效下多场多相渗流理论研究成为深部岩体力学理论的关键问题之一[1-4]，而煤体作为典型的多孔介质，其渗透特性成为煤炭开采技术的研究重点之一。当前应用广泛的煤体渗透率模型主要以应力为自变量[5-7]，仅能描述渗透率随应力增加而降低的过程，难以完整地表征煤体在采动应力下渗透率先减小后增大的演化规律。

针对煤体开采中煤岩体变形与流体耦合的研究，为了能更好地反映多场的相互作用，以Darcy定律为基础增加了渗透率随应力变化的耦合项，从而得到加载应力路径下煤体应力-渗流的耦合规律。李世平[8]通过岩石全应力-应变渗流实验，采用拟合手段获得了3个不同阶段的渗透率与应变的关系。李晓泉和尹光志[9]在考虑三轴压缩状态下损伤的影响条件下，得出了含瓦斯煤样有效体积应力和渗透率存在指数关系的规律，在峰值前后指数系数会发生变化。祝捷等[10]进行了不同气体压力下煤样应力-应变的瓦斯渗流实验，建立了加载煤样变形与渗透率的相关性模型。

为了切实反映地下煤炭开采过程中，煤体具有的特殊应力路径，谢和平等[11]通过分析煤矿典型的3种开采条件下工作面超前支承压力分布规律，凝练出表征煤岩体采动力学的应力路径，为地下工程煤岩体室内实验提供了理论基础。谢和平等[12]通过考虑开采扰动的影响，研究了不同开采方式下煤体应力场-裂隙场-渗流场耦合行为，推导出渗透率增大率与体应变的关系式。为了从机理阐述煤体渗透率的变化，ZHOU等[13-15]采用μCT实验手段、分形理论对煤体内裂隙演化进行了定量描述，并提出采场尺度下不同裂隙网络空间分布与渗透性关系。周宏伟等[16]基于捆绑的火柴棍模型，建立了三向应力条件下煤体渗透率动态演化模型，并进行了实验验证。随着岩体力学理论研究的深入，一些学者研究深部煤体多场耦合行为时，进行了考虑蠕变影响的渗流实验。尹光志等[17]通过含瓦斯煤卸围压蠕变-渗流实验，得出卸围压可使煤体加速破坏从而引起瓦斯流动速度显著增加的结论，并改进Chaboche黏塑性本构模型以描述煤体卸压短期蠕变破坏的规律。许江等[18]利用型煤试样研究了不同温度及不同有效应力条件下三轴蠕变对渗透率的影响规律，得出随着有效应力增大和温度升高，蠕变前后试样渗透率降低的规律。蒋长宝等[19]进行了不同初始围压和瓦斯压力条件下，煤体多级卸围压变形破坏及渗透率演化实验，得出随着围压的降低，煤体蠕变速度和渗透率均加速增大的规律。YANG和ZOBACK[20]为了研究气体吸附对次烟煤的力学及流动特征的影响，通过4种不同气体(He,N2,CH4和 CO2)的吸附和渗流实验，得出CO2，CH4，He对煤体黏塑性的影响依次减弱的规律。CHEN等[21]将Nishihara流变模型引入孔隙弹性的应力-应变本构方程中，推导出包含黏弹性特征的煤体渗透率模型，并通过数值模拟得出蠕变对煤体渗透率和瓦斯抽采具有重要的影响。

深部煤体采动应力条件下必然表现出异于常规三轴岩石实验的特征，同时其特有的渗透率变化规律需要给出理论公式，以便于数值分析和现场应用。笔者在总结常规三轴应力-渗流实验和采动应力下三轴实验规律的基础上，得出深部煤体在采动条件下渗透率随体积应变的变化规律，并对渗透网络在扩展演化中进行合理假设，推导出能反映采动应力条件的渗透率表达式，同时采用已有的实验数据进行了拟合分析。

1 采动应力下渗透率随体积应变变化规律

深部煤体在开采过程中其应力变化路径是典型的轴压增加、围压卸载的过程。文献[11]通过对煤体超前支承压力变化过程分析得出可操作于室内实验的采动应力路径，以模拟煤体在原位开采中应力变化，其实验过程异于传统岩石力学的三轴实验，如图1所示。文献[22]进行了常规三轴压缩实验和采动应力实验下煤体渗透率研究。根据文献[22]得到两种实验条件下煤体应力-应变、渗透率-应变变化规律，如图2所示。常规三轴压缩实验中，煤体经历初始压缩阶段、线弹性阶段、损伤累积阶段、峰值破坏及峰后阶段，在体积应变-轴向应变曲线中表现为损伤起始点对应于体积应变压缩最大值(扩容点)，峰值点对应于体积应变从压缩转变为膨胀。采动应力三轴实验中，煤体按照原位条件被加载至静水压力，然后再根据采动应力路径进行加轴压、卸围压过程，因此应力-应变曲线的起始点为静水压力状态。此时，应力-应变曲线表现为初始裂隙压缩阶段、线弹性阶段、损伤累计阶段、峰值及峰后阶段，其中初始裂隙压缩阶段很短或没有。在采动应力条件下，压缩体积应变很快达到最大值，一般发生在应力-应变曲线初期，而在差应力未达到峰值时，体积应变由压缩转变为膨胀，应力峰值对应于体积应变膨胀段，如图2(b)所示。

图2 煤体应力应变曲线Fig.2 Stress-strain curves of deep coal

图3 常规三轴压缩实验下煤体渗透率演化及裂纹扩展机制Fig.3 Permeability evolution and crack propagation mechanism of coal under conventional triaxial compression test

图4 采动应力下煤体渗透率演化及裂纹扩展机制Fig.4 Permeability evolution and crack propagation mechanism of coal under mining-induced stress test

在常规岩石三轴实验中，渗透率随体积应变的变化趋势大致可分为2个阶段，在岩石损伤起始点之前，渗透率随着体积应变的压缩而逐渐减小，当岩石进入损伤积累阶段时，渗透率随体积应变从压缩变为膨胀呈逐渐增加的趋势，如图3(a)所示。由于煤体内孔隙、裂隙分布错综复杂，为了清晰地阐述煤体裂隙扩展演化过程，将煤体内无规则分布的裂隙通过水平及垂直投影简化为水平裂隙和垂直裂隙。当煤体试样在三向应力下被压缩时，煤基质和孔隙、裂隙被逐渐压缩直至体积应变达到压缩最大值，此间渗透网络被压密导致渗透率下降;煤体试件轴向和环向继续加载，由于内部孔隙、裂隙尖端产生应力集中导致孔隙融合贯通，裂隙迅速扩展形成宏观裂纹，导致渗透率逐渐增加，甚至出现加速增长现象;在应力峰后阶段，尽管轴向和环向应力卸载，但煤基质内弹性能释放导致渗透网络进一步扩展，造成渗透率持续增加，如图3(b)所示。

采动应力条件下，渗透率演化更复杂，以扩容点为界分为两种演化规律。渗透率-体应变曲线中，从起始点至扩容点，渗透率随体应变呈现向下凹的曲线，表现为先增加后降低;扩容点至煤体破坏阶段，渗透率随体应变曲线呈向上凹的形态，表现为先降低再增加的规律，如图4(a)所示。以图4(b)为例，采动应力条件下煤体渗透率演化是煤基质体积变化和孔隙裂隙体积变化共同作用的结果。针对煤体下凹型渗透率形态，原位煤体在未采动条件下处于压缩状态，煤基质和孔隙、裂隙均被压缩(图4(b)中A点)，随着轴向应力增加环向卸载，基质体积被压缩，而孔隙、裂隙因环向卸载而膨胀导致渗透率增加，进一步加轴压卸围压，煤基质发生局部膨胀引起孔隙、裂隙空间被压缩，而孔隙、裂隙体积变化量整体处于减小状态致使渗透率发生下降。在采动应力下，起始状态至扩容点阶段占据渗透率演化的很小部分，因此本文不作为研究重点。针对扩容点至煤体破坏阶段，渗透率呈现出异于常规三轴压缩实验的演化形式，可以分为3个阶段：第1阶段，渗透率随压缩体积应变的减小而减小。该阶段内，煤基质膨胀，原有孔隙、裂隙因煤基质膨胀而被压缩，同时新裂纹扩展，在两种介质变化的博弈状态下，导致压缩体积应变减小，渗透率减小。第2阶段，渗透率处于临界稳定状态，煤基质变化和孔隙、裂隙出现平衡状态。煤基质逐渐从局部膨胀变为整体膨胀而新的裂隙逐渐继续扩展，局部区域仍被压缩，致使渗透率表现为平稳波动状态，但仍高于静水压力下的渗透水平。第3阶段，渗透率随体积应变的膨胀而增加。该阶段内，煤基质整体发生膨胀，裂隙扩展贯通导致煤体渗透网络畅通，当煤体处于峰后状态时，轴向、环向应力均处于卸载状态，导致煤体渗透网络整体处于扩张状态，渗透率持续增高直至由渗流变为流动形态。

2 采动应力下双曲函数型渗透率模型

2.1 基于体积应变的渗透率表达式

为了更好地描述渗透率随体积应变的变化规律，取渗透率比率作为因变量，体积应变作为自变量。由图4可知，渗透率比率呈现上凹的曲线状。本文提出假设:在体积应变空间内，真实渗透网络所呈现的形式在所有可能的渗透网络中使流体渗透耗能最小。

(1)

式中，V为渗透率泛函;a为煤体常量;kr=k/k0，kr为渗透率比率，k为渗透率，k0为初始渗透率;εV为体积应变;εV0为初始体积应变;εVp为终值体积应变。

这是泛函分析中的等周问题，引入拉格朗日乘子，可得拉格朗日量:

(2)

式中，f为拉格朗日函数;λ为拉格朗日乘子。

欧拉-拉格朗日方程可表示为

⟹

(3)

(4)

结合式(3)，(4)，欧拉-拉格朗日方程首次积分得到式(5)，即

(5)

式中，α为常量。

(6)

令akr+λ=αz，微分方程(6)可化解为(7)。

(7)

对式(7)积分可得

(8)

式中，b1为积分常数。

(9)

所以，渗透率比率可表示为体积应变的函数，即

(10)

通过选择合理的坐标系及适当的常数代换，渗透率关于体积应变的函数可表示为

(11)

(12)

文献[23]认为煤岩渗透率与变形(体积应变、轴向应变和环向应变)的关系可采用同一形式的方程式表达，因此，渗透率比率以轴向应变为自变量的方程可以表示为

(13)

式中，c1，c2，d为待定常数;ε1为轴向应变。

2.2 考虑体积蠕变的渗透率表达式

周宏伟等[24]基于分数阶导数理论推导出能反映蠕变加速阶段的流变模型，而文献[25]在此基础上得到以轴向应变为自变量的体积蠕变表达式，通过适当的参数代换可表示为式(14)。将式(14)代入式(12)则可得到体积蠕变条件下以轴向应变为自变量的渗透率表达式(15)。

(14)

(15)

式中，ν为泊松比;h，m，n，p为待定常数。

3 深部煤体渗透率模型验证

3.1 采动应力下渗透率演化模型验证

深部含瓦斯煤体在开采过程中常出现两种扰动环境，其一为保护层开采引起被保护层煤体应力卸压，其二为本煤层开采引起工作面前方煤体卸压破坏。文献[26]利用物理相似试验模型研究了下保护层开采过程中被保护层体积应变变化规律。下保护层开采时，被保护层经历了加载、卸载过程，随之渗透率发生先减小后增大的现象。利用本文渗透率公式(11)进行拟合(图5)，以获得双曲渗透率模型的相关参数，见表1。选取文献[22]中卸围压加轴压渗流实验数据，以体积应变扩容点为起始点，绘制渗透率比率-体积应变图，选择式(12)进行参数拟合，渗透率拟合曲线如图6所示，参数见表1。图5,6中，体积应变负值表示压缩应变，正值表示膨胀应变。

表1 双曲函数型渗透率模型参数拟合结果Table 1 Parameter determination of developed permeability model

图5 保护层开采下煤体渗透率比率-体积应变曲线Fig.5 Permeability ratio with volumetric strain of coal under protective seam mining

图6 卸围压煤体渗透率比率-体积应变曲线Fig.6 Permeability ratio with volumetric strain of coal under unloading confining pressure

从图5，6可以看出，体积应变从压缩转变为膨胀，即体积应变零值，并未对应着渗透率的最低值，而是煤体还处于压缩阶段，在煤基质变形和裂纹扩展共同作用下导致渗透率增加，且这种增长趋势为不可逆行为。

3.2 蠕变条件下渗透率演化模型验证

深部煤体蠕变渗流实验中测得渗透率随轴向应变呈现降低、稳态、增加的规律。采用渗透率比率表达式(13)对文献[27]蠕变条件下渗透率随轴向应变的数据进行拟合，如图7中Ⅰ型渗透率模型所示。利用体积蠕变表达式获得的由轴向应变表示的渗透率方程(15)同样可以对蠕变条件下渗透率数据进行拟合分析，如图7中Ⅱ型渗透率模型曲线所示。两种渗透率模型均能很好地模拟渗透率先减小、稳定、后增加的趋势，相关参数见表2。