分区变权顶板突水危险性评价方法改进及实例

2024-03-31刘传泽谢成龙李宏泽

合肥工业大学学报（自然科学版） 2024年3期

刘传泽, 谢成龙, 詹润, 李宏泽, 孙贵, 韩锋

(1.合肥工业大学资源与环境工程学院,安徽合肥 230009; 2.安徽省煤田地质局勘查研究院,安徽合肥 230088; 3.中煤新集能源股份有限公司,安徽淮南 232001)

0 引言

煤层顶板突水作为煤矿开采的主要灾害之一,是在多种地质因素综合影响下形成的一种非线性动力现象[1-6]。针对煤层顶板复杂、变化、量多的指标控制体系,基于地理信息系统(geographic information system,GIS)的多源信息集成技术与评价方法,目前已得到广泛的应用[7-13]。但以往评价方法在指标权重的确定上相对固定,大多数未考虑各类控制因素在不同组合状态下的相对重要性与均衡程度。为此,近年来一些学者提出了分区变权理论模型,旨在解决常权模型下突水评价、预测的难题[14]。这种理论模型不仅能综合考虑多种指标之间的复杂关系,而且通过对主控因素设定相应的“惩罚”和“激励”机制,能够达到调整各类指标在不同单元权重的目的,是一种更为先进合理的突水评价方法。然而,前人在进行分区变权评价过程中,其常权向量的确定通常采用层次分析法,主观判断性仍然较强,且计算过程相对繁琐。而且,其所构建的状态变权向量对“极大”或“极小”值过于凸显,一定程度上弱化了其余高值或低值的权重。再者,以往分区变权权重与脆弱性指数的计算,仅通过钻孔数据插值获取,常忽略构造因素的影响,没有刻画出因构造指标变化造成的突水情况。

基于上述研究背景,本文以口孜东矿13-1煤层顶板突水危险性为例,引入一种简单、客观的常权赋权方法——熵权法,将该方法与分区变权理论相结合,采用保守型函数改进状态变权向量;同时运用ArcGIS强大的数据处理与空间分析功能,建立地质构造、含水层与隔水层等多源因素耦合的分区变权改进判识模型。将该方法评价结果与仅通过钻孔构建的风险型分区变权模型及常权模型进行比较,确定本次改进模型的合理性与准确性,为煤矿顶板水害防治提供技术方法参考。

1 研究区地质背景

口孜东矿位于阜阳市颍东区与颍上县交界处,地处淮河冲积平原,年均降雨893.74 mm,面积43.447 5 km2。研究区东边界为F12断层,西边界为DF3断层,北边界为陈桥背斜露头区,南边界为阜凤逆冲推覆断裂。其基本构造特征表现为一个不完整的向斜构造。

区内以发育NNE的张性断层为主,地质构造复杂程度中等,13-1煤层构造如图1所示。

图1 研究区13-1煤层构造

13-1煤层是口孜东矿的矿井首采煤层,其直接充水水源为顶板砂岩裂隙含水层,顶板突水危险性较高。但以往对该煤层的突水危险性评价工作较为有限,评价结果与地质事实吻合度不高。该煤层顶板岩性主要为灰色、灰白色细砂岩,具有垂向裂隙。由于区内裂隙发育不均一,导致砂岩富水性差异较大,总体以静态存储量为主。

2 主控因素选择及数据归一化

利用收集到的钻孔抽水试验数据,对13-1煤层顶板各主控因素进行量化处理,借助ArcGIS空间插值功能可直接建立导水裂隙带高度、关键层厚度、含水层厚度、隔水层厚度、脆塑比5个主控因素栅格分布情况。断层强度、断层交叉点与尖灭点密度采用ArcGIS线密度分析与点密度分析工具直接将线文件与点文件量化处理为栅格图。

不同主控因素数值分布的栅格图如图2、图3所示。

图2 不同主控因素数值分布栅格图

图3 断层交叉与尖灭点密度分布栅格图

由于各主控因素分属不同类别,需对各因素量值进行归一化处理。各指标按照对煤层突水控制作用区别,可分为正向与逆向指标。

正向指标与突水呈正相关的关系,即量化值越大,越容易突水;反之,逆向指标主要对煤层突水起抑制作用。

正向指标归一化公式为:

yj=(xj-xj,min)/(xj,max-xj,min)

(1)

逆向指标归一化公式为:

yj=(xj,max-xj)/(xj,max-xj,min)

(2)

其中:yj为归一化处理后数据;xj为归一化前各主控因素的量化值;xj,max和xj,min分别为各主控因素量化值的最大值和最小值;j为样本编号。

利用ArcGIS中模糊分类工具中的线性函数分类功能对主控因素专题图进行归一化处理,得到主控因素归一化专题过程图。该过程图是建立分区变权模型的基础。

3 分区变权模型构建

3.1 熵权法确定常权向量

常规多通过层次分析法确定分区变权模型中的常权向量,但层次分析法操作过程繁琐,且需反复校验调整评价参数。熵权法作为一种客观赋值方法,可依据实际信息量的大小确定影响因素的权重[16-17]。信息量越大的因素,在评价中所起到的作用越大,其权重也越大;相反,信息量越小的因素,在评价中所起到的作用越小,其权重也越小。为排除人为因素影响,本文引入熵权法对分区变权模型进行改进。

在运用熵权法确定各主控因素的权重值时,首先要计算第i项主控因素中第j个数据的比重,其公式为:

(3)

其中:yij为各主控因素归一化后数值;fij为第j个数据占第i项主控因素的比重。

第i项主控因素的信息熵为:

(4)

其中:Hi为信息熵,可反映信息的无序程度;n为主控因素中数据的个数。

第i项主控因素的常权权重为:

(5)

其中:wi0为第i项主控因素的权重;m为主控因素个数。运用熵权法计算的各常权主控因素常权权重见表1所列。

表1 主控因素常权权重

3.2 构建分区状态变权向量

构建变权模型的目的是为凸显主控因素突变时对煤层顶板突水的影响程度。设定合理的“惩罚”与“激励”机制可对主控因素原始权重进行调整,使主控因素指标高值或低值获得更大权重,从而提高评价结果的准确性[18-19]。“惩罚”区间的因素权重随量化值的增大而减小,“激励”区间的因素权重则随量化值的增大而增大。值得注意的是,单一指标值的激增会明显提高顶板突水危险性等级,但单一指标值的显著降低不一定会使顶板突水危险性等级降低。因此,在构建变权模型时应保证“激励”幅度比“惩罚”幅度大,以满足实际评价需求。

不同函数变权向量模型如图4所示。

图4 不同函数的状态变权向量模型

以往的状态变权向量多为风险型函数[20],仅凸显了“极大”或“极小”指标值权重,而“较大”或“较小”指标值的权重并未得到明显加强。由于煤层顶板突水的发生是由量变到质变的过程,且受多种地质因素综合影响,单指标异常也不一定引发突水。因此,仅凸显“极大”或“极小”指标值的权重并不能体现各影响因素对煤层突水的均衡控制作用。本文采用保守型函数对状态变权向量进行改进,使“较大”或“较小”指标值的权重同样获得较大加强,从而使主控因素权重在高、低异常值处变化更为均衡。构建的状态变权向量函数为:

Si(x)=(Si1,Si2,Si3,Si4)

(6)

(7)

其中:Si为第i项主控因素的状态变权向量;d1、d2、d3为各主控因素变权区间阈值;a1、a2、a3、c为状态变权向量的调权参数。经计算与调整,确定c=0.18、a1=0.533、a2=1.592、a3=0.865。

3.3 基于GIS的聚类分析及变权区间确定

在应用分区变权模型对各主控因素的常权权重进行调整时,需对主控因素的指标值进行分区处理。根据指标值的相似性和差异性作出分类,确定“惩罚”与“激励”区间的范围。利用GIS的isodata聚类功能,即迁移平均值法对各主控因素指标值进行迭代运算[21-22]。相较于K-均值聚类法,该算法更加灵活合理,通过设置特定的分裂与合并条件,可在聚类过程中自动增减类别数目。

在聚类过程中,首先将不同参数指标值输入分配至各聚类中心,然后计算输入值在各聚类区间的距离指标函数,再依据指定条件对聚类区间进行分裂或合并处理,从而获得新的聚类区间。通过上述步骤反复迭代运算,直到聚类结果收敛,最终获得各主控因素聚类分区图与分区阈值。

通过20次迭代运算,将栅格图中所有像元分为4类,对分区数据进行提取与归一化处理,得到研究区各主控因素变权区间见表2所列。

表2 主控因素变权区间

3.4 建立变权权重栅格图

基于确定的主控因素变权区间阈值和调权参数,结合式(6),建立研究区13-1煤顶板突水主控因素的变权模型公式[23]为:

(8)

其中:wi为第i项主控因素的变权权重向量;wi0为第i项主控因素的常权权重;Si为第i项主控因素的状态变权向量。

依据分区变权模型,借助ArcGIS软件,确定不同主控因素变权权重值的空间变化。首先,根据表2,利用GIS栅格计算器con函数提取各主控因素不同区间的像元信息,然后使用式(6)、式(8)并结合常权权重wi0对各区间栅格图进行运算、组合、叠加,最终得到的各主控因素变权权重栅格分区如图5所示。

图5 不同主控因素的变权权重栅格图

图5中每个栅格图像元的信息值代表各主控因素专题图对应位置的变权权重。

对比图5与图2、图3可看出,主控因素专题图中正向指标红色区域的权重最高,橙色区与绿色区中等,黄色区与浅绿色区的权重最低。反向指标则为绿色区权重最高,浅绿色与红色区中等,黄色区与橙色区最低。

以上特点说明,各主控因素指标值在高值或低值时权重得到加强,这与变权模型的基本原理相符,侧面反映本次所构建的分区变权模型与计算方法准确性较高。

4 顶板突水脆弱性评价

采用脆弱性指数法[10],将分区变权模型确定的各主控因素变权权重与构造条件、岩性因素相结合,建立研究区顶板突水脆弱性评价模型:

(9)

其中:IV为脆弱性指数;wi为第i各主控因素的变权权重向量;fi(x,y)为第i项主控因素归一化后数据;x,y为地理坐标。

在ArcGIS中使用栅格计算器将变权权重栅格图与主控因素归一化专题图相乘后叠加,可计算出13-1煤层顶板突水脆弱性指数。采用自然间断点分级法对脆弱性指数预测分区,最终得到保守型变权模型顶板突水脆弱性评价结果如图6所示。

图6 保守型变权模型顶板突水脆弱性评价结果

5 模型对比与验证

5.1 模型对比

为进一步证实本次建立的分区变权模型的准确性,将其与常权模型和仅通过钻孔构建的风险型分区变权模型进行对比,常权模型及仅通过钻孔构建的风险型变权模型如图7所示。由对比结果可知,3种模型中各分区空间分布总体趋势相似。其中脆弱区与较脆弱区均呈NNE向展布,与矿井主体断裂走向保持一致,显示断裂构造是造成本区突水的主要因素。矿井南部阜凤断裂与NNE向断层交汇处,脆弱性指数多呈现高值,说明构造交点处应力集中是造成顶板突水的高发地段。

图7 传统的风险型变权模型与常权模型评价结果

在分区结果上,3种模型存在明显的差异。相较于图7b常权模型,图6变权模型缩小了突水较安全区,扩大了较脆弱区与脆弱区的范围。造成这种差异的原因在于构建状态变权向量的“激励”程度大于“惩罚”程度。与处在“惩罚”区间的指标相比,“激励”区间的指标被赋予了更高的权重,因此其安全区范围变小。同时,“强激励”区间的指标被赋予的权重最高,对突水的影响最大,因此脆弱区与较脆弱区范围扩大。此外,图6中过渡区范围较图7a明显减少,较脆弱区与脆弱区的范围增加,其原因是改进的保守型变权模型不仅提高了指标值在“极大”或“极小”值的权重,也提高了“较大”或“较小”指标值的权重,从而增加了脆弱区与较脆弱区的范围,减少了过渡区的范围。对比图6与图7a也可明显看出,利用栅格图构建的评价模型受构造影响更为显著,在断层发育处突水危险性程度普遍提高,这更符合实际地质情况。

从局部地区来看,图6中变权模型A区DF3断层处为较脆弱区,而在常权模型图7b中相应位置却为过渡区。究其原因是该地区关键层厚度比周围薄许多,仅为2 m左右,且该地区导水裂隙带高度远高于周围,均大于55 m。因此,相对于其他顶板突水主控因素,该区域关键层厚度和导水裂隙带高度的权重得到加强,指标突变对评价分区结果的影响得到加强。

对比图7a与图6中的B区可以看出,通过栅格图构建的分区变权模型中,B区DF14断层附近为较脆弱区,但在仅通过钻孔构建的变权模型中却为过渡区。造成这种差异的原因是该区域断层密度较高,通过钻孔数据生成的模型仅能反映岩性场分布特点,无法精确地展示断层的影响和分布。而栅格图数据对断层密集程度的表达具有更好的效果,从而使该区域脆弱性指数相对变高。

5.2 模型验证

利用GIS将口孜东矿13-1煤实际突水点信息投影至本次构建的脆弱性评价模型图6中,共7个点,其中3个点位于较脆弱区,4个点位于过渡区,表明预测模型与实际地质情况一致。位于较脆弱区内的3个点均离断层较近,受构造影响明显,且导水裂隙带高度或含水层厚度较大,隔水层厚度较低,脆弱性指数偏大。位于过渡区内的4个突水点与断层的距离相对较远,构造影响减弱,但其导水裂隙带厚度较大或关键层厚度较小,因而处于过渡区。

尽管通过钻孔数据构建的风险型变权模型及常权模型针对上述7个突水点的评价结果与本次工作基本一致,但两者与保守型变权模型相比,其总体的有效性较低。例如,突水点3本次评价在脆弱区,但在常权模型中则处于较安全区,与地质事实不符;2个较严重的突水点5、6本次评价中位于较脆弱区,但在风险型变权模型中位于过渡区,与地质事实也不符合。综上验证,证明本文构建的分区变权模型相较另2种传统的模型,与实际情况更吻合,具有较高的精度。