张玉军,申晨辉,张志巍,李友伟
(1.天地科技股份有限公司 开采设计事业部,北京 100013;2.中煤科工开采研究院有限公司,北京 100013;3.煤炭资源高效开采与洁净利用国家重点实验室,北京 100013;4.煤炭科学研究总院 开采研究分院,北京 100013)
目前煤炭仍是我国的主体能源,在碳达峰、碳中和目标下,对煤炭绿色开采提出了更高的要求。我国厚煤层可采储量大,其安全、高效、绿色开采对我国煤炭能源的兜底保障作用具有重要意义。随着煤炭科学技术和煤炭开采装备的不断进步,针对厚及特厚煤层的大采高、快速推进及综放等高强度开采技术近年来飞速发展,加之深部矿井、西部特殊地质和生态条件等使得煤矿水害发生机理日趋复杂,防控难度显著增大,水害防控与水资源协同保护成为今后煤炭资源开发面临的重大科学和技术难题。而煤矿开采引起的导水裂缝带高度是水害防控和保水采煤的关键参数,对于煤矿水害防控、地表建构筑物防护、水资源保护都具有极其重要的意义。
国内各矿区应用不同的方法进行了大量针对不同覆岩和开采条件的导水裂缝带高度实测,并分析了导水裂缝带高度的发育特征,总结出了单层采厚在1~3 m,厚煤层分层开采的“两带”高度预计公式,写入到《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规范》中(以下简称《规范》)[1]。许延春等[2-3]利用华北矿区多个矿井中硬和软弱覆岩综放开采导水裂缝带高度数据,回归了综放开采经验公式。但由于厚煤层综放开采或者一次采全高导水裂缝带实测成果相对较少,且不同地质条件下高强度开采导水裂缝带发育高度存在较大的差异性,预计公式在实际应用中还不具有普适性。郭文兵等[3]从厚煤层(大采高)、推进速度快、工作面尺寸大等角度给出了高强度开采的定义;文献[4-6]针对黄陇区域深埋特厚煤层综放开采覆岩导水裂缝带的发育高度开展了一系列的现场实测和理论研究,基于实测数据得出了白垩纪特厚煤层综放开采导水裂缝带高度与采高的关系;文献[7-9]采用钻孔冲洗液、分布式光纤和并行电法等手段实测了软弱覆岩地层大采高工作面覆岩破坏高度,得出了上海庙、新疆矿区软弱覆岩厚煤层综放开采工作面的垮采比和裂采比;李江华等[10]采用钻孔冲洗液消耗量和钻孔电视观测相结合的方法,测得了坚硬和软弱地层高强度开采条件下覆岩破坏裂隙发育特征,得出坚硬顶板高强度开采综合垮采比为9.39~9.62,裂采比大于17.80,约为软弱覆岩的2倍。近年来,随着陕北矿区浅埋特厚煤层的大规模开采,文献[11]研究了榆神府矿区不同隔水岩组厚度与采高比条件下的采动顶板垮落和导水裂缝带发育高度。张玉军[1]从水体下采煤视角,总结了采动覆岩破坏规律和控制技术的发展现状,重点划分了浅埋薄基岩高强度开采导水裂缝带异常发育特征模式。曹祖宝等[12]研究了4 种典型覆岩结构对导水裂缝带发育高度的影响特征;王振荣等[13]针对布尔台煤矿多煤层重复采动条件,实测了导水裂缝带高度。
虽然以上研究针对厚及特厚煤层各矿区进行了导水裂缝带高度实测,但是获得的数据还相对比较零散,大多是针对局部煤矿或矿区的测试与研究成果,缺乏全国范围的高强度开采导水裂缝带分布特征与规律的深入总结研究。为解决这一突出问题,以《建筑物及水体、铁路及主要井巷煤柱留设与压煤开采指南》及笔者单位近年的实测数据以及查阅文献获取的数据为基础,分析了我国厚及特厚煤层高强度开采导水裂缝带发育高度区域分布规律和影响因素,并以此数据为基础回归了导水裂缝带高度预计公式。
系统收集整理了179组国内近年来厚及特厚煤层综放和一次采全高条件下导水裂缝带高度的实测数据,坚硬、中硬和软弱覆岩类型数据分别为51组、100组和28组。数据主要来源于《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采指南》以及查阅文献获取的数据等。通过对实测数据的统计、归类,最终以我国14个大型煤炭基地为区域,分别列出了各大基地煤矿区不同覆岩类型条件下的导水裂缝带高度和裂采比(导水裂缝带高度与采厚的比值)的范围以及平均值(图1和表1)。
从区域上划分,蒙东基地的鹤岗矿区覆岩类型以坚硬为主,裂采比为15.43~22.93,平均18.51;蒙东基地(辽宁矿区)中硬覆岩类型裂采比为9.46~11.23,平均10.23;蒙东基地(扎赉诺尔、敏东、大雁等矿区)软弱覆岩类型裂采比为5.6~11.67,平均9.25。神东基地覆岩类型以中硬为主,基岩厚度足够大的情况下裂采比为9.03~22.06,平均15.52;对于局部埋深浅基岩薄的区域,裂采比集中在19~28,最大达到35。新疆基地覆岩类型以中硬偏软弱为主,裂采比为11.38~13.59,平均为12.24。宁东基地坚硬覆岩类型裂采比为23.05,中硬覆岩类型为12.46,软弱覆岩类型为7.75。陕北基地覆岩类型以中硬为主,裂采比为10.47~19.41,平均14.98。黄陇基地坚硬覆岩类型裂采比为13.33~23.33,平均为18.58,中硬覆岩类型为9.95~13.67,平均为11.76。晋北基地覆岩类型以软弱为主,裂采比为5.91~7.82,平均6.87。晋东基地覆岩类型以坚硬为主,裂采比为15.85~20.15,平均17.94。冀中基地中硬覆岩类型裂采比为11.08~16.1,平均13.97,软弱覆岩类型裂采比为4.45~8.25,平均6.35。鲁西基地中硬覆岩类型裂采比为7.52~14.77,平均10.61,软弱覆岩类型裂采比为6.8~9.69,平均8.67。河南基地中硬覆岩类型裂采比平均17.58,软弱覆岩类型裂采比平均8.29。两淮基地坚硬覆岩类型裂采比为15.74~20.8,平均17.67,中硬覆岩类型裂采比为8.85~14.66,平均为11.45,软弱覆岩类型裂采比为6.83~10.1,平均8.08。
总体来看,目前现有的实测数据基本涵盖了我国的主要产煤基地,能够为全国各大矿区的安全生产提供一定的依据。
依据实测数据统计发现,不同覆岩条件下厚及特厚煤层高强度开采的导水裂缝带高度分布规律有较大差异,见表2和图2。通过对比分析可知,实测平均裂采比随着覆岩由坚硬、中硬到软弱依次降低,坚硬覆岩条件下的平均裂采比值是中硬覆岩的1.52倍,是软弱覆岩的2.29倍,而中硬覆岩平均裂采比值则是软弱覆岩的1.5倍。
依据实测数据绘制了不同覆岩类型条件下导水裂缝带高度以及裂采比分别与采厚的关系图,如图3和图4所示。从图3中可以看出,坚硬、中硬和软弱覆岩条件下的导水裂缝带高度整体上均随着采厚的增加而呈现逐步增大的趋势,尤其是采厚达到6 m时,增加的速率明显加快,且坚硬覆岩条件下增加的速率明显要高于中硬和软弱覆岩。由图4可知,坚硬覆岩的裂采比明显高于中硬和软弱覆岩,但裂采比均随着采厚的增加而呈现逐渐降低的趋势,且减小速率随着采厚达到一定厚度后有趋于稳定的趋势。
图3 不同覆岩类型导水裂缝带高度与采厚关系Fig.3 Relationship between coal seam thickness and water-conducting crack zone height
图4 不同覆岩类型裂采比与采厚关系Fig.4 Relationship between fracture-mining ratio and mining thickness of different ovelying rock types
14个大型煤炭基地生产能力达到全国的95%以上,其中神东、陕北、黄陇、蒙东等矿区煤炭赋存主要以厚及特厚煤层为主。因此,重点分析以下3个典型矿区的覆岩破坏特征。
黄陇煤炭基地包括彬长(含永陇)、黄陵、旬耀、蒲白、澄合、韩城、华亭矿区。其中彬长矿区是我国新兴的西部煤炭生产基地,地处鄂尔多斯盆地南缘,主采煤层为侏罗系延安组含煤地层,厚度大,是我国大规模的综放式开采矿区,覆岩类型为中硬。从实测数据中整理了彬长矿区7个煤矿25组数据,同时选取全国东部矿区同类覆岩(采厚大于5 m)24个煤矿33组数据进行统计对比分析,得出该矿区导水裂缝带发育高度具有以下3个方面特征:
1)导水裂缝带发育高度、裂采比与采厚、工作面长度的关系,总体上具有一定的规律性,随着煤层采厚的增加,导水裂缝带发育高度随之增加,但增加趋势逐步减缓;随着采厚的增加,裂采比则呈降低趋势(图5、图6)。
图5 彬长及东部矿区导水裂缝带高度与采厚关系Fig.5 Relationship between height of water-flowing fracture zone and mining thickness in Binchang and eastern mining area of China
图6 彬长矿区导水裂缝带高度、裂采比与工作面长度关系Fig.6 Relationship between height of water-flowing fractured zone,fracture-mining ratio and working face length in Binchang Mining Area
2)中硬覆岩条件下,国内东部矿区在采厚5.0~13.4 m,平均7.1 m的条件下实测平均裂采比12.4,而彬长矿区在采厚6.0~13.5 m,平均9.4 m的条件下实测平均裂采比为17.3,同类覆岩条件下采厚相同,导水裂缝带高度和裂采比明显高于东部矿区,其原因在于彬长矿区的开采强度(工作面尺寸、推进速度)明显高于东部矿区。
3)从图6可以看到,工作面长度以170 m为分界点,成为影响彬长矿区导水裂缝带高度发育的主要影响因素。工作面长度大于170 m时导水裂缝带发育高度受工作面长度影响较大,其高度及裂采比随着工作面长度增加,呈明显增高增大趋势,且增加速率加快。
蒙东(东北)基地主要包括内蒙古东部的呼伦贝尔市、伊敏、通辽、赤峰平庄、兴安盟和锡林郭勒盟,同时也包括辽宁阜新、铁法矿区,以及黑龙江鸡西、鹤岗、双鸭山等矿区。覆岩分为内蒙古东部的软弱-极软弱类型和鹤岗矿区的中硬-坚硬类型。
2.2.1 蒙东基地极软弱覆岩矿区
扎赉诺尔矿区顶板岩石单向抗压强度多小于20 MPa,且固结性差,属于极软弱类覆岩。对于该类覆岩条件,由于采动裂缝被压实和闭合的时间很短,且软弱覆岩对导水裂缝具有很好的抑制作用,因此很难实测获得其最大高度,因此收集到的实测数据主要是基于现场开采实践并根据采动影响范围内含水层的水位变化而确定的,如该矿区的铁北煤矿在采厚3.2 m时推断导水裂缝带高度约29 m,裂采比约为9.1;灵东矿在采厚不超过11.60 m时,第三含水层距煤层超过62 m的情况下,导水裂缝带没有直接波及含水层,推断裂采比小于5.6;灵东煤矿采厚按照2.5~16.0 m、导水裂缝带高度与采厚的比值约为5~10进行预计,矿井涌水量无明显增加,实现了首采工作面含水层下的安全回采[14-15]。由此,提出采用趋势分析方法来预计极软弱覆岩条件下的导水裂缝带发育高度,图7为绘制的趋势图,由图7可知,该类覆岩条件下裂采比明显小于国内软弱覆岩条件下的裂采比,最小裂采比仅为5,采厚达到16 m时仅为全国软弱覆岩平均裂采比的68%。
图7 极软弱覆岩条件下的导水裂缝带发育高度及趋势Fig.7 Development height and trend of water-flowing fracture zone under extremely weak overburden
2.2.2 蒙东基地鹤岗矿区坚硬覆岩
鹤岗矿区煤层覆岩平均抗压强度为58.3~89.8 MPa,属于中硬-坚硬覆岩类型,整个矿区开采具有断层发育、开采煤层多、累计开采厚度大、地层受反复扰动等特征。收集整理了该矿区8个煤矿不同条件下的导水裂缝带高度实测数据,见表1和图8。由实测数据分析可知,在采厚为3.5~12.0 m,裂采比绝大多数介于16~23,但是受断层影响相同采厚条件下导水裂缝带发育高度差异较大。无断层存在的条件下平均裂采比为15.26,断层存在条件下平均裂采比为17.69,断层构造条件下综放开采覆岩破坏高度增大了16%。
图8 鹤岗矿区导水裂缝带高度及裂采比分布Fig.8 Height and ratio distribution of water-flowing fracture zone in Hegang Mining Area
神东基地位于陕西以北榆林地区、内蒙古东部东胜地区,以神府煤田、东胜煤田为主。该矿区具有“埋深浅、基岩薄、土层厚”高强度开采的特点,因而大部分矿井仅存在 “两带”,无弯曲下沉带,甚至常发生切冒、抽冒、台阶下沉等非均衡性覆岩破坏形式。受以上条件因素的影响,实测获得的导水裂缝带高度差异性极大(图9、图10),在采厚3.00~9.06 m,中硬覆岩条件下导水裂缝带高度实测值整体离散性较大,不具有规律性。在采厚3~6 m时,裂采比集中在19~28,最大达到35,明显高于全国尤其是东部矿区中硬覆岩条件的裂采比。但是部分矿井由于厚土层的存在,尤其是厚层红土层对导水裂缝带发育具有明显抑制作用,实测裂采比7.16~7.78,较类似软弱覆岩条件下开采裂采比仍降低了15%~53%。因此,该矿区导水裂缝带高度特征有以下3个方面:
图9 神东矿区实测导水裂缝带高度、裂采比与采厚关系Fig.9 Relationship between measured water-flowing fractured zone height,fracture-mining ratio and mining thickness in Shendong Mining Area
图10 神东矿区实测裂采比柱状Fig.10 Column chart of measured fracture-mining ratio in Shendong Mining Area
1)基岩厚度足够大,“两带”发育完整、充分,其中的导水裂缝带和垮落带高度、分布形态和破坏程度均发育正常,但是由于采深较小,受基岩风化带和土层的影响导水裂缝带高度将随着采深的减小而降低[16-17]。
2)基岩厚度足够小且开采高度较大时,覆岩(土)层被垂直裂缝完全贯穿,导水裂缝带甚至直接发育至地表,此时实测获得的裂采比通常较大,接近甚至超过坚硬覆岩条件下实测获得的裂采比。
3)基岩厚度较小时,但存在较厚的隔水性较好的土层,导水裂缝带高度将受到黏性土层的抑制,“两带”高度会大幅度降低,接近软弱覆岩条件实测获得的裂采比。
从上述各个典型矿区及全国导水裂缝带实测数据发育趋势来看,总体具有一定的规律性,但也同时发现不仅在同一矿区即使同一矿井实测值有时差别较大,分析其原因主要有以下5个方面:
1)地层及覆岩结构的差异性影响。同一矿区甚至是同一矿井覆岩结构存在着差异,比如神东矿区实测导水裂缝带高度受地层结构、埋深等因素影响显著。
2)导水裂缝带高度实测方法确定标准存在差异性。目前,导水裂缝带高度实测的主要方法有地面钻孔冲洗液漏失量、井下仰上钻孔注水试验、钻孔电视、物探、微震以及现场水位观测等。各方法确定裂缝带顶点的标准不同,导致实测结果存在较大的差异。据初步统计,以上收集的实测数据中采用传统的钻孔冲洗液漏失量观测方法占到了65%以上,该方法钻孔从地面施工,自上而下从完整岩层逐步过渡到采动影响的破碎岩层寻找冲洗液漏失量的突变值,依据裂隙的导水性确定裂缝顶点,容易把握,取得的数据可靠且真实。
3)实测位置的差异性。对于地面钻孔观测方法而言,钻孔位置设计决定了导水裂缝带高度实测值是否为最大值。
4)实测时间的差异性。实测时间视地层与工作面推进情况而定,一般滞后于工作面对应位置回采后1~2个月可获得导水裂缝带最大高度。
5)综放开采放煤厚度的不均匀性影响。对于综放开采,放煤高度是否均匀决定了采厚的波动,从而会导致导水裂缝带高度的变化,甚至会出现台阶状分布规律。
根据收集的实测数据对导水裂缝带高度与采厚的关系分别按照线性和分式函数形式进行拟合,并求得公式的中误差。其中分式函数拟合过程如下:将原始实测数据按照式(1)做相对应的变换,规范经验公式形式为
(1)
(2)
(3)
得到3种覆岩条件导水裂缝带高度预计公式,见表3。由表3可知,无论是线性拟合还是分式函数拟合,导水裂缝带高度均符合线性和分式函数增长规律,尤其是坚硬覆岩拟合决定系数R2达到0.90~0.94,其次为中硬和软弱覆岩,R2也分别达到了0.81~0.84。
表3 厚及特厚煤层高强度开采导水裂缝带预计公式Table 3 Prediction formula of water-flowing fractured zone in high strength mining of thick and extra thick coal seam
现将回归得到的线性和分式函数预计公式、以及《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采指南》(简称《指南》)中的经验公式一[18]计算结果与不同采厚和覆岩条件下的实测值进行对比,见表4。
表4 不同覆岩导水裂缝带计算结果(部分)Table 4 Calculation results of different overburden water-flowing fractured zone (part)
续表
由表5和图11—图13可知,坚硬覆岩条件下用回归的线性公式计算得到的导水裂缝带高度与实测值差值率绝对值平均8.5%;按照分式函数关系回归的计算公式计算得到的导水裂缝带高度与实测值差值率绝对值平均差值率5.54%。中硬覆岩条件下分别为12.20%、6.62%和5.73%。软弱覆岩条件下分别为26.99%、11.94%和9.93%。综上分析,回归得出的厚及特厚煤层高强度开采导水裂缝带公式明显优于《规范》中的经验公式,具有良好的适用性,且分式函数公式优于线性函数公式。
图11 坚硬覆岩导水裂缝带高度预计公式计算结果与实测值对比Fig.11 Comparison between predicted calculation results and measured values of height of water-flowing fractured zone in hard overburden
表5 计算导水裂缝带高度与实测值绝对差值率Table 5 Absolute difference rate between calculated fracture zone height and measured value
图12 中硬覆岩导水裂缝带高度预计公式计算结果与实测值对比Fig.12 Comparison between predicted calculation results and measured values of height of water-flowing fractured zone in medium hard overburden
图13 软弱覆岩导水裂缝带高度预计公式计算结果与实测值对比Fig.13 Comparison between predicted formula and measured value of water-flowing fractured zone height of weak overburden
1)基于179组国内厚及特厚煤层综放或者大采高综采条件下导水裂缝带高度实测数据,从我国14大煤炭基地区域角度分析了不同覆岩类型条件下的导水裂缝带高度和裂采比的范围和分布特征,初步绘制了中国煤矿导水裂缝带高度分布图,对我国的主要产煤基地高强度开采条件下导水裂缝带发育高度区域分布规律有了新的认识。
2)覆岩类型和采厚是影响导水裂缝带高度的重要因素。实测平均裂采比随着岩性由软弱、中硬到坚硬依次递增,坚硬覆岩的平均裂高分别是中硬覆岩和软弱覆岩的1.52倍和2.29倍,中硬覆岩平均裂高是软弱覆岩的1.5倍;相同覆岩类型的导水裂缝带高度都随着采厚的增加而增加,裂采比则随着采厚的增加而降低,且减小速率随着采厚达到一定厚度有趋于稳定的趋势。
3)从矿区地层及覆岩结构的赋存特征、导水裂缝带高度实测方法确定标准、实测时间和位置等方面,分析了实测导水裂缝带高度存在较大差异性的原因。黄陇基地永陇-彬长矿区导水裂缝带高度和裂采比,明显高于同类覆岩和采厚条件下的东部矿区,其原因在于彬长矿区的开采强度(工作面尺寸、推进速度)明显高于东部矿区;蒙东基地扎赉诺尔矿区覆岩条件下裂采比明显小于国内软弱覆岩条件下的裂采比,最小裂采比仅为5;神东基地实测获得的导水裂缝带高度差异性极大,实测结果与工作面的埋深、基岩厚度以及上覆黏性土层的采动隔水性能密切相关。
4)基于统计的实测数据,按照线性和分式函数2种形式回归得到了厚及特厚煤层高强度开采导水裂缝带高度预计公式,预计结果明显优于《规范》中的经验公式,具有良好的适用性,且分式函数公式优于线性函数公式,可为全国不同覆岩类型条件下厚及特厚煤层高强度开采导水裂缝带高度预计提供依据。