库拜煤田中低阶煤中孔-大孔发育控制因素

2022-06-06秦全灵田继军

煤炭科学技术 2022年4期

秦全灵,李鑫，吴斌，田继军

(1.新疆大学地质与矿业工程学院，新疆乌鲁木齐 830047；2.新疆维吾尔自治区煤田地质局一六一煤田地质勘探队，新疆乌鲁木齐 830009)

0 引言

孔隙结构影响着煤储层吸附、解吸、扩散及渗透率，是煤储层研究的核心问题之一[1-2]。我国煤层受多期次构造运动影响，构造煤发育，煤孔隙大小分布在毫米至纳米级，孔隙结构较复杂，不同尺度孔隙控制着煤层气的解吸、扩散与渗流等过程[3]。库拜煤田是我国南疆缺煤地区重要的煤炭资源产地与煤层气资源开发潜力区，但针对库拜煤田煤储层孔隙特征方面的研究相对较少，难以指导当地煤层气储层评价。

研究煤中孔隙的控制地质因素可间接促进煤层气的商业化开采。SAKUROV等[4]研究了澳大利亚煤的孔隙分布，发现孔隙受煤阶影响较大，并且在某些孔径范围内，镜质组含量会对孔隙产生影响。陈向军等[5]选取了不同矿区、不同变质程度的煤样进行了压汞测试，结果表明煤的变质程度对大孔体积产生影响较大。张洲等[6]对库拜煤田主要煤储层的孔隙系统进行了研究，发现主要煤储层的中孔和大孔较发育，并且孔隙连通性较好，有利于煤层气的运移。安庆等[7]对库拜煤田中部储层孔隙特征进行了研究，发现不同类型裂隙对煤层气运移有不同程度的影响，大裂隙是渗透率的主要贡献者，也是决定煤层气是否具有开采价值的关键因素。

库拜煤田煤层气资源潜力巨大，其商业开采可为当地经济发展提供长期的能源保障。通过开展高压压汞试验、工业分析、煤岩组分测定、显微组分定量分析、煤体宏观描述等实验和分析，围绕库拜煤田中低阶煤中孔-大孔发育及其控制地质因素展开研究，讨论了镜质组含量、成熟度、煤体结构等对中孔及大孔孔隙体积及结构复杂程度的影响。研究成果有助于揭示库拜煤田高渗储层发育特征并促进当地煤层气开发。

1 地质背景

库拜煤田位于塔里木盆地的北缘，面积约1 585 km2，属于塔里木盆地北部的一个次级构造单元-库车山前中新生代坳陷(图1)，区内含煤地层为向南倾斜的单斜构造，倾角为30°～89°，库车县境内煤层倾角较缓，拜城县境内地层倾角较大。煤田以北紧邻深大断裂形成的F1逆断层，以南为2条深大断裂，F2和F3正断层，以及大型含气褶皱带构造。区内的岩性主要由石英砂岩、粗砂岩、砂砾岩、细砂岩、粉砂岩和煤层组成，中生界侏罗系的塔里奇克组厚度为177～256 m，阳霞组厚度为200～571 m，侏罗系的克孜努尔组厚度为294～838 m，是库拜煤田的主要含煤地层(图2)。由于不同程度的区域变质作用和动力变质作用，使得煤样煤化程度跨度较大，由东向西煤变质程度逐渐增加[8]。

图1 库拜煤田构造纲要[9](据文献[9]修改)

图2 库拜煤田煤系柱状图[10](据文献[10]修改)

2 试验

2.1 取样及样品基本信息

样品取自于库拜煤田的铁西、阿艾、铁东3个矿区(图1)9口煤层气参数井，根据煤体结构的不同(图3)开展了45次高压压汞试验。样品包含原生结构煤15块、碎裂结构煤19块、碎粒结构煤3块和糜棱结构煤8块。铁西矿区宏观煤岩类型主要为暗淡煤、半暗煤和半亮煤，属低变质到中等变质阶段，埋深介于582.00～969.25 m，分布于A1、A3、A5～A10、B、B1、C和C2煤层。阿艾矿区宏观煤岩类型主要为半暗煤、半亮煤和亮煤，属中等变质阶段，埋深介于951.60～1 103.75 m，分布于下5和下10煤层。铁东矿区宏观煤岩类型主要为暗淡煤、半暗煤和半亮煤，属低变质到中等变质阶段，埋深介于533.95～952.96 m，分布于A3、A5、A7、A8、A9-10和A10煤层。

图3 不同煤体结构煤样照片

对库拜煤田不同煤体结构的煤岩样品进行煤质基础信息测试，其中水分(Mad)、挥发分(Vdaf)和灰分(Ad)测试遵照GB/T 212—2008国家标准进行；煤的镜质组最大反射率(Ro，max)利用显微镜在油浸物镜下测定，遵照GB/T 6948—2008；煤岩显微组分定量在偏光显微镜下进行，遵照GB/T 8899—1998。煤质基础测试结果显示样品Mad变化范围较小，为0.34%～1.68%，平均1.10%；Vdaf变化范围较大，为6.22%～38.84%，平均24.14%；Ad变化范围较大，为1.38%～44.5%，大部分在10%以下，个别样品达到44.5%，平均8.64%；镜质组最大反射率Ro，max为0.30%～1.96%，平均1.04%；镜质组含量为20.6%～93.6%，变化范围较大，平均74.9%；惰质组含量为5.0%～79.4%，平均24.8%；壳质组含量为0～1.4%，平均0.3%。

2.2 高压压汞试验

高压压汞法是在一定压力下将液态汞压入煤岩中，压入汞的量能代表孔隙的体积，通过压力与压入汞的体积的关系曲线可以得到一系列表征样品孔隙特征的参数[11]。压汞法操作简单有效，试验时间较短，一直被广泛应用于多孔材料的孔隙特征和分布等方面的研究[12]。孔隙分类方法采用苏联学者霍多特[13]孔径分类划分方案，即微孔孔径<10 nm、过渡孔孔径10～100 nm、中孔孔径100～1 000 nm、大孔孔径>1 000 nm。本次试验遵照GB/T 21650.1—2008《压汞法和气体吸附法测定固体材料孔径分布和孔隙度》执行。在进行高压压汞试验时，随着注入压力的增大，样品的原生孔隙结构会被破坏，微孔、过渡孔测定会出现较大误差，但是对中孔和大孔等较大孔隙的分析较准确[14-15]。压汞误差来源主要有2个方面：①高压阶段煤基质压缩和未进汞孔隙压缩，导致压汞测得微孔、过渡孔体积偏大[16]。但本次研究未分析微孔和过渡孔，主要分析中孔和大孔；②低压阶段的麻皮效应，使得汞在样品表面的粗糙处消耗部分汞，导致低压段进汞量增大，几十微米以上的孔隙异常发育[17]。针对此，采取粒径3～6 mm表面平直煤颗粒作为样品，并设置最低汞注入压力为0.01 MPa，尽可能弱化麻皮效应导致的误差。此外，由于高压压汞可实现快速、低成本地对大量样品中孔、过渡孔进行分析，因此本次工作采用了高压压汞法，并且采取上述措施后尽可能地减小了压汞试验误差。

3 中孔和大孔发育程度地质控制

通过高压压汞试验获得了样品的孔隙体积、退汞效率等参数，结果显示：测试样品的大孔体积为0.000 2～0.400 0 cm3，平均0.040 1 cm3，中孔体积为0.001 6～0.786 0 cm3，平均0.087 5 cm3，退汞效率为15.43%～80.08%，平均47.98%。中孔体积约是大孔体积的两倍，孔体积变化范围较大，相差2～3个数量级，退汞效率变化范围较大，总体较低。

3.1 煤岩组分对中孔、大孔体积的控制作用

由图4可知，随镜质组含量的增加，中孔体积、大孔体积以及中孔、大孔体积和呈现微弱增大趋势，煤岩组分对中孔、大孔孔隙体积具有一定控制作用。吕志发等[18]早在1991年就对块煤的孔隙特征及其影响因素做了研究，在分析煤的孔隙特征与有机显微组分时，发现镜质组孔隙主要是微孔和小孔，中孔和大孔主要发育在惰质组中。刘娜等[19]在研究煤岩孔隙度主控地质因素时，发现中大孔比例与孔隙度呈良好正相关关系，而镜质组含量又与孔隙度呈负相关关系，由此推断出镜质组含量与中孔、大孔体积呈负相关关系。图4所示的结果与前人研究成果不一致，可能原因是当镜质组含量增多时，镜质组生烃留下大量原生孔隙，在成岩演化过程或外力的作用下，原生微孔和过渡孔破坏并部分转换为中孔和大孔，从而导致高镜质组含量煤中孔、大孔体积较高[20]。此外，图4中出现了异常高点，均为糜棱煤，姜家钰[21]通过压汞法得出不同构造程度的煤体压汞曲线特征，发现随构造作用的增强，各阶段的孔容都有所增加，大孔和中孔增幅分别为537.5%和1 300%，构造作用对煤中大孔和中孔有明显改造。因此，经过构造作用形成的糜棱煤出现了较高的中孔和大孔体积。

图4 中孔及大孔体积随镜质组含量变化

3.2 成熟度对中孔、大孔体积的控制作用

图5分别表示了成熟度对中孔体积、大孔体积、中孔与大孔体积和的控制作用，结果显示：成熟度对中孔、大孔、中孔大孔体积和的控制作用并不明显；成熟度与中孔大孔体积呈现出弱相关性。成熟度增高对有机质孔的发育有促进作用[1]，刘娜等[19]认为成熟度对孔隙影响较为复杂，从低成熟度阶段到高成熟度阶段平均中大孔所占比例呈现先减小再增大再减小的趋势。

大数据理念指的是在数据关联分析的基础上进行预测，针对电网规划体系中存在的问题，本文结合大数据技术提出了新型的电网规划体系，这一体系主要由三个部分组成，即数据获取、处理以及应用。大数据时代下的新型电网规划体系，能够根据电网规划体系的具体需求，来获取电力系统内外的全面数据，并运用相应的数据处理技术，从而实现电网方案规划、风险评估防范、工程实时监测以及用户交互等项目的智能化、信息化。新型电网规划体系架构如图1所示。

杜玉娥[22]在研究煤孔隙特征对煤层气解吸的影响时发现，当Ro，max大于1.3%时，煤芳香化程度显著提高，这种分子结构形成了一系列微孔和小孔，同时也导致大、中孔减少，随成熟度增加大孔体积出现先减小再增大的结果与图5b中所示基本符合。图5中出现异常高点均为糜棱煤，Ro大致介于张吉振等[23]划分的成熟生烃溶蚀阶段(0.7%

图5 中孔及大孔体积随成熟度变化

3.3 煤体结构对中孔、大孔体积的控制作用

由图6可知，原生结构煤、碎裂煤、碎粒煤中孔、大孔体积接近，糜棱煤中孔、大孔体积明显高出。秦修培等[25]在研究构造变形对煤孔隙发育特征影响时发现随着构造变形程度的增加，煤体孔径结构分布发生变化，开放孔逐渐增多，半封闭孔逐渐减少，中孔和大孔发育逐渐增多。李明等[26]发现随着构造变形程度的加强，从原生结构到碎粒结构，各孔径段的孔容分布相对均匀，而到糜棱煤，出现中孔异常发育且大孔发育较好现象，这与研究结果有较好的一致性。

图6 煤体结构与中孔及大孔体积关系箱型图

4 孔隙结构复杂程度地质控制

分形理论能够准确地定量描述孔隙的内部复杂结构，多孔固体介质中分形维数越大，孔隙结构越复杂，非均质性越强[27]。对于用压汞法测得的数据可以采用Friesen等引入的Menger海绵模型[28]，并结合Washburn方程处理可以得公式

(1)

其中：V为孔隙体积，cm3；P为压力，MPa；D为分形维数，可通过下式计算：

D=k+4

(2)

其中：k为y与x函数的斜率，y为孔隙体积对压力的积分取对数lg(dV/dP)，x为对压力取对数lgP。利用公式分析压汞数据并构建lg(dV/dP)与lgP的分布散点图。由于压汞数据众多，限于篇幅在图7展示了部分样品的分形维数计算拟合过程，计算出中孔的分形维数D3和大孔的分形维数D4(表1)，据此分别表征中孔和大孔孔隙结构复杂程度。经计算可知中孔的分形维数D3为2.30～3.93，平均3.19；大孔的分形维数D4为2.10～3.77，平均2.98。

表1 中孔、大孔分形维数分布

图7 不同样品lg(dV/dP)与lg P的分布

4.1 镜质组含量对孔隙结构复杂程度的控制

从褐煤到长焰煤的成煤过程中，褐煤中的腐植体逐渐转化为镜质组，按照我国煤炭行业标准MT/T 1160—2011煤的镜质组含量分级，试验样品属于中高镜质组煤。样品中孔分形维数大部分大于3，大孔分形维数也有相当一部分大于3，这可能是由于压缩效应导致煤的孔隙结构遭到破坏引起的[29]。由图8a和图8b可以看出，样品镜质组含量占比表现为集中分布，大都在70%以上，并且这部分镜质组含量与中孔、大孔的复杂程度呈正相关关系；镜质组含量在70%以下的数据较少，本文不做讨论。镜质组含量大于70%时，中孔和大孔的分形维数D3、D4随镜质组含量的增加呈增大的趋势。这一现象可能由以下原因导致：当镜质组含量较低时，其自身较均质化孔隙所占比例较高，孔隙复杂程度相应较低；当镜质组含量较高时，大量脆性的镜质组受构造影响可能产生别的类型的孔隙，煤岩中孔和大孔的随机分布性和复杂性相对于镜质组较低时有所加强[20,30]。

图8 孔隙复杂程度随镜质组含量变化

4.2 成熟度对孔隙结构复杂程度的控制

图9 孔隙复杂程度随成熟度变化

4.3 煤体结构对孔隙复杂程度的控制

根据煤岩破坏的程度，可将煤样分成4种类型，分别是原生结构煤、碎裂结构煤、碎粒结构煤和糜棱结构煤[32]。图10展示了库拜煤田不同煤体结构煤岩中孔分形维数和大孔分形维数，可知不同煤体结构对应的中孔分形维数变化不大，表明库拜煤田煤岩中孔结构复杂程度并不受煤体结构的控制。由图10b可知，原生结构煤和碎裂煤的大孔结构复杂程度接近，糜棱煤与碎粒煤的大孔结构复杂程度较高，表明煤体受破坏程度影响着大孔结构的复杂程度，结合4.2节可知，成熟度不是影响大孔复杂程度的主要因素，而煤体结构对大孔的复杂程度影响较大。

图10 孔隙复杂程度与煤体结构关系箱型图

原生结构煤(图11a、11b、11c)和碎裂结构煤(图11d、11e、11f)压汞曲线划分为Ⅰ类，进汞和退汞曲线形态大致相似，基本没有分离段，退汞效率为39.60%～76.96%。

碎粒结构煤压汞曲线(图11g、图11h、图11i)可被划分为Ⅱ类，进退汞曲线形态差异较大，出现了明显的分离段，退汞效率为20.31%～39.65%。糜棱结构煤压汞曲线(图11j、图11k、图11l)可被划分为Ⅲ类，进退汞曲线形态相似程度不高，但是出现了阶段性的突变特征，退汞效率为47.11%～58.36%。随着构造作用的加强，煤岩孔容的比例也在发生转化，中孔和大孔所占比例有不同程度的增加，孔隙形态和类型也在发生变化，从而增加了中孔大孔的复杂程度，进而使得糜棱煤的退汞效率相对于原生结构煤较低[25]。

图11 不同煤体结构煤压汞曲线特征

图12展示了原生结构煤、碎裂结构煤、碎粒结构煤和糜棱结构煤的孔径分布特征，相比之下，构造煤大孔更发育；结合不同煤体结构压汞曲线，随煤体结构复杂程度增大，压汞曲线滞后回线逐渐明显，表明煤岩中孔隙形态发生了转变，从管状或平板状形态为主转向以墨水瓶形态为主[33]。构造煤促进了大孔的发育，也增加了孔隙形态的多样性，因此，与大孔孔隙的复杂程度表现出了正相关关系。