《中国地质灾害与防治学报》
由于矿山的地质变化情况越来越难以掌握,矿山地质灾害发生的特征呈现出多样性以及复杂化的趋势[1,2]。针对矿山地质灾害预测是有效预防地质灾害的重要手段,而矿山地质灾害预测系统成为了有关部门的调研重点。在我国,针对矿山地质灾害预测系统的研究中,尽管研究起步较早,在实际应用中存在预测精度低的问题,导致研究普遍存在局限性,很难在真正意义上实现矿山地质灾害精准预测[3]。物探磁法的出现恰好解决了出现的问题,物探磁法具备一定的空间数据基础表现能力,主要在计算机设备与良好运营环境的支撑下,采集相关磁性数据,并根据磁异常的表达特征对其实施有序的数据管理、运算、分类等功能。目前,物探磁法已经广泛应用到各个领域中,尤其是在地质勘探中占据首要位置。因此,本文进行基于物探磁法的矿山地质灾害预测研究,并通过分析其在实际当中的应用情况,证明应用的有效性。通过将物探磁法应用在矿山地质灾害预测方法设计部分,致力于为矿山地质灾害预测方法的优化设计提供更加广阔的发展空间。
1 物探磁法
物探磁法指的是以矿石中的磁性特征为基础,分析磁性特征之间的具体差异。物探磁法作为一种地球物理勘探方法,能够通过地磁场变化特征判别矿石的特征。根据物探磁法的特性,只要是针对地下含有磁性矿物的各种岩石以及其他磁性物体,都可以通过物探磁法测量出具有不同的剩余磁性以及感应磁性。这些剩余磁性以及感应磁性很容易造成程度不同的磁场异常,叠加在正常的地磁场上。通过物探磁法使用仪器对剩余磁性以及感应磁性进行测量,能够有效获得地面磁异常的特征情况,从而勘查矿石地质信息,致力于实现地质勘查的目的。根据物探磁法对勘查相对均方误差的要求,勘查相对均方误差不能超过0.8,否则就会被视为无效勘查[2]。物探磁法主要用于对弱磁性矿物的勘查以及隐伏磁性体在地表产生的弱磁异常研究等工作,适用于勘查、找矿、预测等多领域中。
2 基于物探磁法的矿山地质灾害预测方法
基于物探磁法的矿山地质灾害预测方法整体流程,如图1所示。
图1 矿山地质灾害预测方法流程图
结合图1所示,下文针对图中3步主要流程进行详细研究。
2.1 物探磁法划分矿山地质灾害易发分区单元网格
本文利用物探磁法,划分地质灾害易发分区单元网格。根据物探磁法的基本应用原理,将每个单元面积为设置为1km×2 km×3km×3km。使用仪器对网格中的磁性物体进行物探,基于物探磁法查明矿山地磁场变化特征,将地磁场变化特征以矩形阵势进行排列。设排列后的地磁场变化特征矩形阵记为(X,Y),则其计算公式,如公式(1)所示。
在公式(1)中,n指的是基于物探磁法测量得到的地磁场变化信息特征值;i指的是地磁场变化信息维数;h指的是预测时刻;a指的是实际地质横坐标;b指的是实际地质纵坐标。通过公式(1)得到排列后的地磁场变化矩形阵,进行单元格重组划分,设置变量系数,提取矿山地质灾害预测数据,地磁场变量系数越高证明矿山地质灾害发生的几率也就越高。本次划分矿山地质灾害易发分区单元网格必须遵循以下三点原则:第一点是平面网格系统划分采用国际上认可度高的角点网格系统;第二点为平面网格系统划方向与主模拟层断层方向保持一致;第三点为表层模型中储层平面非均值性差异,保证两口钻井相隔之间至少有三个空点网格[3]。以此为标准,映射矿山地质灾害预测数据到矿山地质灾害易发分区单元网格中,为下文计算矿山地质灾害预测映射误差提供依据。
2.2 计算矿山地质灾害预测映射误差
由于矿山地质灾害不确定因素大的特点,在映射矿山地质灾害预测数据的过程中,必然会产生映射误差[4]。为确保矿山地质灾害预测映射误差能够通过修正系数进行改进,且过程中不会出现权重比例失调的情况出现。将矿山地质灾害预测映射误差设为MSE,可得MSE的计算公式,如公式(2)所示。
在公式(2)中,r指的是矿山地质灾害预测映射态势值;D指的是矿山地质灾害预测特征维数。剔除计算得出的映射误差,提高矿山地质灾害预测映射精度。
2.3 预测矿山地质灾害
将剔除映射误差后的矿山地质灾害预测数据按变量权重大小与矿山地质图像信息进行集合,实现矿山地质灾害预测的数据叠加[5]。以数据为依托,构建矿山地质灾害预测方程式为,则其计算公式,如公式(3)所示。