一种采用模糊数学理论圈定采空区电阻率影像异常的方法

摘要

本发明涉及一种采用模糊数学理论圈定采空区电阻率影像异常的方法，包括通过电阻率法采集的采空区电阻率影像数据，其特征在于，对于数据的处理包括如下步骤：1)数据评估，构建剖面数据量计算公式，数据直方图统计分析评价；2)数据增强，确定采空区电阻率影像异常是否存在的量化评价指标DR、DT和DH；3)DR的提取和赋值运算；4)DT的提取和赋值运算；5)DH的提取和赋值运算；6)DR、DT和DH的综合评价；7)采空区电阻率影像异常区域的定量评价及采空区异常范围的定量圈定，供潜孔钻钻探实时验证。本发明的优点是：利用模糊数学理论，使采空区异常圈定定量化和智能化，实现了海量地球物理异常数据的快速定量评价。

说明书

技术领域

本发明属于金属矿山安全生产领域，具体涉及一种采用模糊数学理论圈定采空区电阻率影像异常的方法。

背景技术

在金属矿山安全生产领域，特别是露天金属矿山，经常会出现一些由于不明采空区的变形、沉降和塌陷引发的生产安全事故，轻则造成设备损失，重则造成人员伤亡，严重影响矿山的安全生产作业，成为时刻威胁露天金属矿山生产安全的主要难题。为了解决这些不明采空区存在诱发的生产安全技术难题，科研院所和生产单位采用多种技术方法和手段探测采空区，且主要是大量采用高效快捷的地球物理探测方法，尤其是其中的高密度电阻率法。电阻率法以其成本低、方便快捷高效和无损，是非接触的地质CT探测技术方法，属于超前的“绿色探测”，适应于对面积性露天采场区域内的不明采空区异常快速扫面筛查和定位。但电阻率法成果显示采空区异常具有间接性、多解性和圈定指标不确定性，这往往也成为物探成果难以及时服务生产的一大限制因素和缺陷，使得采空区的成果解译变成了“有”或“无”的难以判定，这使得探测成果的解译和应用必须要有丰富经验的地质与地球物理专家级的人才联合参与，从而严重制约地球物理这种快捷高效技术方法的推广与应用。与此同时，面临的现实情况是在露天矿山快速持续不断的平台作业生产中，每天采集的大量电阻率法影像数据却难以得到具有丰富经验的专家级人才及时处理，致使采矿作业平台采空区异常边界得不到及时准确圈定，造成钻探验证和采空区安全处理严重滞后，在制约矿山正常生产的同时威胁矿山安全，使得大量不明采空区分布区的铁矿资源被禁锢，无法安全开采与释放。这一切问题的关键在于，如何使物探解译成果的模型化和智能化，最大限度的摆脱专家人力解译缓慢和滞后性，从而实现计算机自动处理。

为此，面对上述应用地球物理问题仅依靠专家经验解决的局限性，本发明基于模糊数学理论、采空区异常模型和计算机科学技术的进步发展，特别是目前深凹露天金属矿山面临充水型采空区的安全威胁，提出了一种基于模糊数学理论的充水型采空区电阻率法异常定量圈定方法。该方法主要基于多年的BIF型铁矿山采空区探测技术成果支撑，特别是不明充水型采空区的地球物理特征研究成果的总结，应用模糊数学理论研究和模型评价的方法，基于GIS平台构建的采空区空间地球物理异常影像快速处理与定量评价模型，实现金属矿山不明采空区地球物理异常圈定的模型化和智能化。

发明内容

本发明的目的在于提供一种采用模糊数学理论圈定采空区电阻率影像异常的方法，该方法可以快速高效的评价采空区电阻率影像异常，及时准确高效圈出生产区域内的采空区异常，供潜孔钻钻探验证，而后安全处理，全力保障矿山生产安全。

本发明的的目的是通过以下技术方案实现的：

本发明的一种采用模糊数学理论圈定采空区电阻率影像异常的方法，包括通过高密度电阻率法测量仪采集的采空区电阻率影像数据，其特征在于，对于采空区电阻率影像数据的处理包括如下步骤：

步骤1、电阻率影像数据评估

对利用高密度电阻率法测量仪采集到的电阻率剖面影像数据，利用测量电极布设数和数据采集层数来构建电阻率影像数据可采总量公式，对单个剖面采集到的R个电阻率数据进行直方图统计分析，在其影像数据分布符合正态分布或对数正态分布时方可进入下一步骤进行影像数据处理；

步骤2、电阻率影像数据增强，确定采空区电阻率影像异常是否存在的三个量化评价指标

对电阻率影像数据进行矩阵中值滤波降噪增强和数据反演计算，根据计算结果提出确定采空区电阻率影像异常是否存在的三个量化综合评价指标，分别为采空区电阻率影像异常高低评价指标D

步骤3、电阻率影像异常高低评价指标D

求出单个剖面电阻率法R个影像数据的均值

步骤4、采空区电阻率影像异常空间变化梯度评价指标D

对单个剖面的电阻率法R个影像数据采用离散傅里叶变换垂向求导获得负向梯度值圈闭区域，并使其全面正值化；利用Jenks自然间断点分级法进行归一化10级空间分类，利用模糊数学理论中的模糊高斯模型

步骤5、采空区电阻率影像异常埋深条件指标D

对单个剖面R个影像数据点所在的深度由浅到深进行Jenks自然间断点分级法的归一化10级空间分类；D

步骤6、采空区电阻率影像异常的三个量化指标综合评价

采用模糊数学中的FuzzyAnd模型FuzzyAnd＝min(arg1,...,argn)对D

步骤7、采空区电阻率影像异常区域的定量评价及采空区异常范围的定量圈定

依据上述采空区电阻率影像异常的三个量化评价指标的选取原则，为了无遗漏采空区异常，最大限度保证安全生产，在采空区电阻率影像异常的概率分布影像图D

另外，在步骤1中，获取的高密度电阻率法数据必须进行直方图统计分析和采空区异常降噪增强；在步骤2～步骤6中，对于提取充水型采空区电阻率影像异常的三个定量评价指标D

与现有技术相比，本发明的优点是：

(1)提取了充水型采空区探测过程中电阻率值大小、空间分布和变化所隐含的采空区异常信息，实现了采空区地球物理异常信息的挖掘。

(2)把采空区定性解译过程中的专家经验和依据的采空区地球物理特征信息量化，实现了采空区数据的重构和归一化处理。

(3)利用模糊数学理论，使采空区异常圈定定量化和智能化，可以实现海量地球物理异常数据的快速定量评价。

附图说明

图1、电阻率法测线G-1原始数据图。

图2、电阻率法测线G-1的矩阵降噪滤波和数据反演成果图。

图3、电阻率法测线G-1电阻率异常高低评价指标D

图4、电阻率法测线G-1电阻率异常空间变化梯度评价指标D

图5电阻率法测线G-1电阻率异常埋深条件指标D

图6电阻率法测线G-1中充水型采空区电阻率异常的概率分布影像图D

具体实施方式

下面结合附图对本发明作进一步说明。

具体实施方式为对某一露天铁矿的充水型采空区电阻率影像异常定量圈定。

如图1～图6所示，本发明的一种采用模糊数学理论圈定采空区电阻率影像异常的方法，包括通过高密度电阻率法测量仪采集某一露天铁矿的充水型采空区电阻率影像数据，其特征在于，对于采空区电阻率影像数据的处理包括如下步骤：

步骤1、电阻率影像数据评估

对利用高密度电阻率法测量仪采集到的电阻率剖面影像数据，依据单条测线的测量电极布设数DJn和影像数据采集层数Cn，可知获取了R个影像数据(R＝((DJn-3)+(DJn-3Cn))*Cn/2)，如图1显示的电阻率法测线G-1为60个测量电极和19层数据采集，则其采集影像数据量为570个。然后对单条测线剖面中的R个电阻率影像数据(x

步骤2、电阻率影像数据增强，确定采空区电阻率影像异常是否存在的三个量化评价指标

对上述可靠的电阻率影像数据进行3×3的矩阵中值滤波降噪增强和数据反演，通过电阻率法测线G-1的处理(见图2)，明显突出了潜在的充水型采空区异常及其边界，然后根据充水型采空区存在的地球物理异常特征，提出决定采空区电阻率影像异常是否存在的三个量化指标如下：分别是采空区电阻率异常高低评价指标D

步骤3、电阻率异常高低评价指标D

首先求出一个测线剖面R个电阻率影像数据的均值

步骤4、采空区电阻率影像异常空间变化梯度评价指标D

充水型采空区地球物理特征上为典型负向梯度值圈闭，即从充水型采空区异常的边界到核心，电阻率逐渐减少，对测线剖面的电阻率数据采用离散傅里叶变换垂向求导获得负向梯度值圈闭区域，并使其全面正值化。利用Jenks自然间断点分级法对负向梯度数据进行归一化10级空间分类，指定x

步骤5、采空区电阻率影像异常埋深条件指标D

该条件指标可以有效地将浅部接触不良导致的异常影像数据和可信度低的深边部探测影像数据有效筛除，以降低采空区异常误判概率，对R个影像数据点所在的深度由浅到深进行Jenks自然间断点分级法的归一化10级空间分类。D

步骤6、采空区电阻率影像异常的三个量化指标综合评价

采用模糊数学中的FuzzyAnd模型FuzzyAnd＝min(arg1,...,argn)对D

步骤7、采空区电阻率影像异常区域的定量评价及采空区异常范围的定量圈定

依据上述采空区电阻率影像异常的三个量化评价指标指标的选取原则，为了无遗漏采空区异常，最大限度保证安全生产，可以依据电阻率异常高低评价指标D

表1充水型采空区电阻率影像异常判定表

由图6可确定具体采空区异常范围的大小。图6一共圈出两处采空区电阻率影像异常，左侧影像异常明显小于右侧影像异常，且左侧采空区异常概率在0.5-0.7之间，集中于0.6，属于二级采空区异常，需要钻孔验证；右侧采空区异常概率在0.5-0.95，集中于0.8，属于一级采空区异常，需要钻孔重点验证和树立警示牌。

另外，在步骤1中，获取的高密度电阻率法数据必须进行直方图统计分析和采空区异常降噪增强；在步骤2～步骤6中，对于提取充水型采空区电阻率影像异常的三个定量评价指标D

本发明提取了充水型采空区探测过程中电阻率值大小、空间分布和变化所隐含的采空区异常信息，实现了采空区地球物理异常信息的深度挖掘；而且把采空区定性解译过程中的专家经验和依据的采空区地球物理特征信息量化，实现了采空区数据的重构和归一化处理。另外，本发明利用模糊数学理论，使采空区异常圈定定量化、自动化和智能化，可以实现海量地球物理异常数据的快速定量评价。