采集宿州矿区朱仙庄煤矿(Z)、芦岭煤矿(L)、祁东煤矿(D)主采煤层8号煤层及10号煤层的原生结构煤和碎软煤样片5块，其中4块不同变形程度的碎软煤分别为典型的初碎裂煤、碎裂煤、片状煤和碎斑煤，1块原生结构煤用于对比分析。采样过程依据煤层煤样采取方法(GB/T 482-2008)和煤岩样品采取方法(GB/T 19222-2003)进行，选取不同变形程度的碎软煤，用塑料盒装取以防样品破坏。样品制备：选择样品中典型构造部分，依次用240目、480目、600目、1000目、1500目、2000目、3600目和5000目砂纸打磨，将样品打磨成表面光滑，上下底面水平的圆柱状，其中圆柱的上、中、下直径误差小于0.1 mm，从而可以准确标定样品体积。样品基本特征见 表1 。

宿州矿区初碎裂煤Z1，煤岩条带清晰可见，以亮煤为主；层理可见，原生结构构造保存完好，内生裂隙可以分辨( 图2(a) )。显微镜下，原生结构保存较好，发育少量显微构造裂隙，且以张裂隙为主，裂隙的成组性和方向性较明显( 图2(a') )。

碎裂煤Z2较初碎裂煤Z1变形加强，发育有碎裂结构，原生结构、构造尚可识别，无定向性，节理发育较密集，且相互交切，内生裂隙尚可分辨( 图2(b) )。显微镜下，原生结构保存较好，发育多组显微构造裂隙，垂直裂隙、顺层裂隙和斜交裂隙均有发育，以张裂隙为主；构造裂隙切割煤岩条带以及内生裂隙，构造裂隙间以微裂隙相连通( 图2(b') )。

片状煤Z3，宏观特征：层理清晰可见，镜煤和亮煤中可见内生裂隙；片状结构，可捏碎为碎片和碎粒( 图2(c) )。显微特征：主要受到剪切构造应力作用，发育一组平行排列、延伸稳定的剪切裂隙，裂隙多平直稳定、平行排列，裂隙间常见其他方向的细微裂隙发育( 图2(c') )。

碎斑煤D1，煤体变形较强，肉眼可见碎斑结构，节理十分发育，内生裂隙较难分辨；煤体易破碎为碎块( 图2(d) )。显微镜下，煤体破碎较为严重，煤体裂隙发育密集，不同方向裂隙多呈小角度相交，将煤体切割成碎块状、碎斑状，裂隙规模较大，延伸不稳定( 图2(d') )。

渗透率实验是在中国矿业大学煤层气资源与成藏过程教育部重点实验室进行。样品经砂纸打磨成圆柱状，利用实验室孔隙度测试仪进行煤样孔隙度测试，然后再采用美国生产的PDP-200型渗透率仪进行渗透率测试。PDP-200型渗透率仪采用非稳态法的压力脉冲衰减测量技术。其测量范围为：0.00001~10 md。测试方法：将样品放入岩芯，拧紧上下压阀，并关闭围压加油阀，开始加围压，其具体围压加载方案如 表2 ，气体压力设定为2 MPa。

注：测试气体压力均为2 MPa，红色部分为降压；由于煤样存在裂缝，若加载较小围压，无法准确测出其相应的渗透率，故各煤样初始围压不同。

本文根据李传亮建立的多孔介质的有效应力体系 [21] ，根据其计算公式：

${\sigma }_e=P_c-\phi P_m$(1)

式中：σ_e——有效应力，MPa；

P_c——围压，MPa；

$\phi$——为孔隙度；

P_m——平均压力，即进口压力与出口压力的平均值，MPa。

由上式计算每个煤样每个压力点的有效应力，并拟合有效应力值与渗透率之间的相互关系，如 图3(a)~(e) 。

在加载过程中，宿州矿区碎软煤的有效应力与渗透率呈负指数关系(R²= 0.92~0.99)，即随着有效应力的增加，煤储层的渗透率不断减小( 图3(a)~(e) )，原生结构煤的有效应力与渗透率相关性较碎软煤差。就初始渗透率而言，碎斑煤最高，但随着有效应力的增大，其渗透率迅速下降，变化幅度较大( 图3(e) )；片状煤渗透率变化幅度相比碎斑煤、碎裂煤以及初碎裂煤最小；原生结构煤初始渗透率最低，随着有效应力的增加，渗透率变化幅度不大。

初碎裂煤与原生结构煤孔隙度相差不大( 图3(f) )，但初碎裂煤的渗透率是原生结构煤渗透率的9430倍( 图3(a) 和 图3(b) )，表明这与煤样中发育的裂隙以及孔裂隙连通性有关，初碎裂煤受剪切应力的作用，发育多组构造裂隙，为瓦斯渗流提供通道，从而使渗透率大大增加。从初碎裂煤到碎裂煤到碎斑煤再到片状煤，随着煤体变形程度的增加，煤基质的孔隙度和渗透率大体呈增大趋势，但也有反常。碎斑煤孔隙度较小( 图3(f) )，但渗透率最大( 图3(e) )，猜测随着构造应力的不断增大，煤储层发生严重变形，部分内生孔在构造应力下闭合，但是构造裂隙大大发育，且孔隙连通性较其他碎软煤好，所以碎斑煤渗透率最大。片状煤渗透率较初碎裂煤、碎裂煤以及碎斑煤低( 图3 )，表明片状煤孔裂隙连通性差，导致其渗透率低。所以宿州矿区碎斑煤及碎裂煤发育区，煤与瓦斯突出危险性高，但也是煤层气开发和利用的有利区。

1) 根据石油天然气行业标准(SY/T 5336，SY/T 5358，SY/T 6385)对储层敏感性进行评价，引入渗透率损害率以及渗透率损害系数来反应在有效应力作用下煤储层渗透率损害的百分数，应力敏感性引起的渗透率损害系数D_kp [22] ，即：

$D_{kp}=\frac{K_i-K_{i+1}}{\left|P_{i+1}-P_i\right|}$(2)

式中：D_kp为渗透率损害系数，mD/MPa；K_i、P_i、K_i+1及P_i+1分别为i、i + 1点的渗透率和有效应力值，mD、MPa。

煤层渗透率的影响因素十分复杂，地质构造、应力状态、煤层埋深、煤体结构、煤岩煤质特征以及煤阶及天然裂隙等都不同程度地影响着煤层渗透率，因此可通过定义渗透率对有效应力的敏感系数，即：

${\alpha }_k=-\frac{1}{K_0}\frac{\partial K}{\partial P}$(3)

式中：K₀为初始渗透率，mD。为更直观的描述煤储层有效应力对气相渗透率的影响，引入无因次渗透率，即K/K₀(煤样渗透率与初始渗透率的比值)。

当围压从2.07 MPa增加到6 MPa时，片状煤和原生结构煤的无因次渗透率分别为0.099、0.242，片状煤渗透率降低为初始渗透率的9.9%，原生结构煤降低为初始值的24.2%；当围压从10 MPa增加到14 MPa时，初碎裂煤的无因次渗透率为0.123，其渗透率较初始渗透率下降12.3%，碎裂煤无因次渗透率为0.521，其渗透率降低为原始值的52.1%。总体而言，碎裂煤随着有效应力的增强衰减较快( 图4(a) )，表明连通裂隙随着有效应力的增强逐渐趋于闭合，这与碎裂煤具有较好的孔隙配置有关。本文根据李传亮建立的多孔介质的有效应力体系 [21] ，根据其计算公式：

从渗透率损害系数变化特征来看，在有效应力作用下，碎斑煤的渗透率损害系数最高且随着有效应力增加减小较快，其损害系数介于0.3148~0.1287 md/MPa (平均0.1496 md/MPa)；原生结构煤损害系数较小且变化不大，其损害系数介于0.000025~0.00000063 md/MPa (平均0.000009 md/MPa)；有效应力对于初碎裂煤和碎裂煤具有相似的损伤规律，其损害系数分别介于0.2635~0.0011 md/MPa (平均0.0100 md/MPa)，0.2529~0.0392 md/MPa (平均0.1483 md/MPa)；对于片状煤(平均0.0907 md/MPa)而言，随着有效应力的增大，渗透率损害系数骤减(0.2461~0.0184 md/MPa)，随后，随着有效应力的增大，渗透率损害系数变化不大(0.0044~0.0006 md/MPa)。总体而言，从初碎裂煤–碎裂煤–碎斑煤，随着变形程度的增加，碎软煤渗透率损害系数增加，说明煤体变形程度越高，其压力敏感性越强。

2) 不同变形煤降压渗透率变化

由 图5 可知，随着有效应力的减小(绿色曲线)，碎软煤的渗透率增大。其中，碎斑煤降压前后渗透率变化最大，降压后的煤体渗透率远小于降压前的煤体渗透率( 图5(b) )，碎裂煤次之( 图5(a) )，片状煤降压前后渗透率几乎没有变化( 图5(c) )，表明不同煤体变形程度的碎软煤压力敏感性碎斑煤 > 碎裂煤 > 片状煤。其中碎斑煤发生强烈的脆性变形，围压增加，孔裂隙趋于闭合，煤体发生塑性变形；围压减小，煤体无法恢复原始状态，部分孔裂隙依然为闭合状态，渗透率较降压前小。所以碎斑煤发育区，煤储层渗透率更容易受地应力的影响，碎裂煤次之，片状煤影响最小。