首页 / 专利库 / 电脑编程 / 算法 / 一种煤岩动力灾害多系统多参量集成综合预警方法及系统

一种岩动灾害多系统多参量集成综合预警方法及系统

阅读:26发布:2022-10-01

专利汇可以提供一种岩动灾害多系统多参量集成综合预警方法及系统专利检索,专利查询,专利分析的服务。并且本 发明 提供一种 煤 岩动 力 灾害多系统多参量集成综合预警方法及系统,该方法包括:获取多个煤岩动力灾害监测系统的监测数据,基于监测 数据挖掘 出各监测系统中能反映煤岩动力灾害前兆信息的多元特征参量;基于多元特征参量,筛选出每一监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值;基于各监测系统的预警效能最高的特征参量组合及其中各特征参量的最优临界值,计算出每一监测系统的综合预警指数和预警效能;基于各监测系统的综合预警指数和预警效能,计算出煤岩动力灾害多系统综合预警结果。本发明将多个煤岩动力灾害监测系统优势互补,提高了预警准确率并剔除了主观因素的影响,可为矿山灾害防治工作的开展起到一定指导作用。,下面是一种岩动灾害多系统多参量集成综合预警方法及系统专利的具体信息内容。

1.一种岩动灾害多系统多参量集成综合预警方法,其特征在于,所述煤岩动力灾害多系统多参量集成综合预警方法包括:
获取多个煤岩动力灾害监测系统的监测数据,基于所述监测数据挖掘出各监测系统中能反映煤岩动力灾害前兆信息的多元特征参量;
基于所述多元特征参量,筛选出每一监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值;
基于各监测系统的预警效能最高的特征参量组合及其中各特征参量的最优临界值,计算出每一监测系统的综合预警指数和预警效能;
基于各监测系统的综合预警指数和预警效能,计算出煤岩动力灾害多系统综合预警结果。
2.如权利要求1所述的煤岩动力灾害多系统多参量集成综合预警方法,其特征在于,在计算出煤岩动力灾害多系统综合预警结果后,所述方法还包括:
将所述煤岩动力灾害多系统综合预警结果与煤岩动力灾害分级预警评判准则相对比,确定动力灾害危险级别。
3.如权利要求1所述的煤岩动力灾害多系统多参量集成综合预警方法,其特征在于,所述监测系统包括在线监测系统和便携监测系统;其中,所述在线监测系统包括微震监测系统、地音监测系统和液压支架监测系统。
4.如权利要求3所述的煤岩动力灾害多系统多参量集成综合预警方法,其特征在于,所述监测数据包括所述在线监测系统实时采集得到的连续监测数据和所述便携监测系统采集得到的离散监测数据。
5.如权利要求4所述的煤岩动力灾害多系统多参量集成综合预警方法,其特征在于,所述多元特征参量的维度包括时间维度、空间维度和强度维度,所述多元特征参量包括:
所述微震监测系统中反映时序信息的频次、频次比和频次偏差值;反映空间信息的离散度;反映强度信息的微震能量、微震能量偏差值和离散度;
所述地音监测系统中反映强度信息的地音能量、地音能量偏差值、地音能量平均值和脉冲因子;反映时序信息的脉冲;
所述液压支架监测系统中反映时序信息的闪烁间隔危险度;反映强度信息的记录频次和支架压力。
6.如权利要求5所述的煤岩动力灾害多系统多参量集成综合预警方法,其特征在于,所述筛选出每一监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值,包括:
将同一监测系统中属于不同维度的特征参量进行两两组合;
利用遗传算法训练并优选出每一监测系统中适应度值最高的特征参量组合并计算出其中各特征参量的最优临界值;
将各监测系统中适应度值最高的特征参量组合作为该监测系统对应的预警效能最高的特征参量组合。
7.如权利要求6所述的煤岩动力灾害多系统多参量集成综合预警方法,其特征在于,在利用遗传算法训练并优选出每一监测系统中适应度值最高的特征参量组合时,对于所述微震监测系统优选出适应度值最高的前两组特征参量组合;对于所述地音监测系统和液压支架监测系统,优选出每一传感器对应的适应度值最高的一组特征参量组合。
8.如权利要求7所述的煤岩动力灾害多系统多参量集成综合预警方法,其特征在于,所述基于各监测系统的预警效能最高的特征参量组合及其中各特征参量的最优临界值,计算出每一监测系统的综合预警指数和预警效能,包括:
基于各监测系统的适应度值最高的特征参量组合及其中各特征参量的最优临界值,根据下式计算出各监测系统的单系统预警程度:
其中,对于所述微震监测系统,单系统预警程度为适应度值最高的每一特征参量组合的预警程度,此时Wi表示第i组适应度值最高的特征参量组合的预警程度,ai和bi分别表示第i组适应度值最高的特征参量组合中两个特征参量的实时值,V1和V2分别表示第i组适应度值最高的特征参量组合中各特征参量对应的最优临界值;
对于所述地音监测系统和液压支架监测系统,单系统预警程度为各传感器的预警程度,此时Wi表示第i个传感器的预警程度,ai和bi分别表示第i个传感器对应的适应度值最高的特征参量组合中两个特征参量的实时值,V1和V2分别表示第i个传感器对应的适应度值最高的特征参量组合中各特征参量对应的最优临界值;
基于各监测系统的单系统预警程度,根据下式计算出各监测系统的综合预警指数W综:
其中,对于所述微震监测系统,Ri为第i组适应度值最高的特征参量组合的适应度值;对于所述地音监测系统和液压支架监测系统,Ri为第i个传感器对应的适应度值最高的特征参量组合的适应度值;
基于各监测系统的综合预警指数,根据下式计算出每一监测系统的预警效能RI,RI表示第I个监测系统的预警效能:
其中, 为在监测时间段内产生预警并且报准的次数;N1为监测时间段内大能量或冲击事件的总次数;t0为产生预警占用时长;T0为监测总时长;当监测系统对应的综合预警指数在冲击或大能量事件发生前预设天数内超过预设阈值则预警为报准,否则为漏报或误报。
9.如权利要求8所述的煤岩动力灾害多系统多参量集成综合预警方法,其特征在于,所述基于各监测系统的综合预警指数和预警效能,计算出煤岩动力灾害多系统综合预警结果,包括:
基于各监测系统的综合预警指数和预警效能,根据下式计算出煤岩动力灾害多系统综合预警结果I综:
其中,W综I表示第I个监测系统的综合预警指数,RI表示第I个监测系统的预警效能。
10.一种煤岩动力灾害多系统多参量集成综合预警系统,其特征在于,所述煤岩动力灾害多系统多参量集成综合预警系统包括:
多元特征参量挖掘模,用于获取多个煤岩动力灾害监测系统的监测数据,并基于所述监测数据挖掘出各监测系统中能反映煤岩动力灾害前兆信息的多元特征参量;
特征参量组合筛选模块,用于基于所述多元特征参量,筛选出每一监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值;
单系统预警指数及预警效能计算模块,用于基于各监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值,计算出每一监测系统的综合预警指数和预警效能;
综合预警结果计算模块,用于基于各监测系统的综合预警指数和预警效能,计算出煤岩动力灾害多系统综合预警结果。

说明书全文

一种岩动灾害多系统多参量集成综合预警方法及系统

技术领域

[0001] 本发明涉及煤岩动力灾害预警技术领域,特别是指一种煤岩动力灾害多系统多参量集成综合预警方法及系统。

背景技术

[0002] 煤岩动力灾害是煤矿主要的灾害类型之一,其破坏性大,严重时导致井巷闭合,大量煤岩体喷出,造成人员伤亡和财产损失。由于煤炭长期开采,现阶段大多数矿井已进入深部开采,伴随而来的是煤岩动力灾害频次、强度和破坏程度均呈上升趋势。煤岩动力灾害成因复杂多变,给监测预警和防治带来挑战。因此,准确高效的监测预警方法成为煤矿提前了解灾害危险,采取防治措施,确保安全生产的主要方法。
[0003] 目前,矿山中对于煤岩动力灾害的实时监测手段主要包括微震监测、地音监测、支架应力监测和电磁辐射监测等,其中微震监测系统主要针对低频高能事件,可实现对高尺度煤岩断裂的时序监测和对破裂源的空间定位,用于大范围区域监测;地音监测系统主要针对高频低能事件,可实现对低尺度煤岩体微破裂的局部监测,传感器监测范围达30m左右;支架应力监测系统主要针对采场顶板来压情况进行监测;电磁辐射监测系统可同时对煤岩体内部应力集中状况和破裂进行监测,监测范围达20m。然而,由于各种装备在监测原理、监测方式和监测范围等方面存在差异,上述多种监测装备一般各自为政、多源监测数据的预警结果可能不一致,导致仪器装备使用效率不高,甚至产生漏报、误报。同时,由于缺乏对各监测系统监测数据的深入挖掘,预警准确率较低,严重影响井下安全生产。

发明内容

[0004] 针对上述现有技术所存在的缺陷,本发明提供一种煤岩动力灾害多系统多参量集成综合预警方法及系统,对各监测系统的监测数据进行深入挖掘分析,并将得到的监测预警信息进行有机融合,构建多系统多参量综合集成监测预警模型,形成预警准则统一、预警指标统一、预警临界值统一的综合集成高可靠性预警方法,提高对煤岩动力灾害危险的预警精度,减少漏报、误报次数,为矿井灾害防治工作的开展起到一定指导作用。
[0005] 为解决上述技术问题,本发明提供一种煤岩动力灾害多系统多参量集成综合预警方法,所述煤岩动力灾害多系统多参量集成综合预警方法包括:
[0006] 获取多个煤岩动力灾害监测系统的监测数据,基于所述监测数据挖掘出各监测系统中能反映煤岩动力灾害前兆信息的多元特征参量;
[0007] 基于所述多元特征参量,筛选出每一监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值;
[0008] 基于各监测系统的预警效能最高的特征参量组合及其中各特征参量的最优临界值,计算出每一监测系统的综合预警指数和预警效能;
[0009] 基于各监测系统的综合预警指数和预警效能,计算出煤岩动力灾害多系统综合预警结果。
[0010] 进一步地,在计算出煤岩动力灾害多系统综合预警结果后,所述方法还包括:将所述煤岩动力灾害多系统综合预警结果与煤岩动力灾害分级预警评判准则相对比,确定动力灾害危险级别。
[0011] 可选地,所述监测系统包括在线监测系统和便携监测系统;其中,所述在线监测系统包括微震监测系统、地音监测系统和液压支架监测系统。
[0012] 可选地,所述监测数据包括所述在线监测系统实时采集得到的连续监测数据和所述便携监测系统采集得到的离散监测数据。
[0013] 可选地,所述多元特征参量的维度包括时间维度、空间维度和强度维度,所述多元特征参量包括:
[0014] 所述微震监测系统中反映时序信息的频次、频次比和频次偏差值;反映空间信息的离散度;反映强度信息的微震能量、微震能量偏差值和离散度;
[0015] 所述地音监测系统中反映强度信息的地音能量、地音能量偏差值、地音能量平均值和脉冲因子;反映时序信息的脉冲;
[0016] 所述液压支架监测系统中反映时序信息的闪烁间隔危险度;反映强度信息的记录频次和支架压力。
[0017] 进一步地,所述筛选出每一监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值,包括:
[0018] 将同一监测系统中属于不同维度的特征参量进行两两组合;
[0019] 利用遗传算法训练并优选出每一监测系统中适应度值最高的特征参量组合并计算出其中各特征参量的最优临界值;
[0020] 将各监测系统中适应度值最高的特征参量组合作为该监测系统对应的预警效能最高的特征参量组合。
[0021] 进一步地,在利用遗传算法训练并优选出每一监测系统中适应度值最高的特征参量组合时,对于所述微震监测系统优选出适应度值最高的前两组特征参量组合;对于所述地音监测系统和液压支架监测系统,优选出每一传感器对应的适应度值最高的一组特征参量组合。
[0022] 进一步地,基于各监测系统的预警效能最高的特征参量组合及其中各特征参量的最优临界值,计算出每一监测系统的综合预警指数和预警效能,包括:
[0023] 基于各监测系统的适应度值最高的特征参量组合及其中各特征参量的最优临界值,根据下式计算出各监测系统的单系统预警程度:
[0024]
[0025] 其中,对于所述微震监测系统,单系统预警程度为适应度值最高的每一特征参量组合的预警程度,此时Wi表示第i组适应度值最高的特征参量组合的预警程度,ai和bi分别表示第i组适应度值最高的特征参量组合中两个特征参量的实时值,V1和V2分别表示第i组适应度值最高的特征参量组合中各特征参量对应的最优临界值;
[0026] 对于所述地音监测系统和液压支架监测系统,单系统预警程度为各传感器的预警程度,此时Wi表示第i个传感器的预警程度,ai和bi分别表示第i个传感器对应的适应度值最高的特征参量组合中两个特征参量的实时值,V1和V2分别表示第i个传感器对应的适应度值最高的特征参量组合中各特征参量对应的最优临界值;
[0027] 基于各监测系统的单系统预警程度,根据下式计算出各监测系统的综合预警指数W综:
[0028]
[0029] 其中,对于所述微震监测系统,Ri为第i组适应度值最高的特征参量组合的适应度值;对于所述地音监测系统和液压支架监测系统,Ri为第i个传感器对应的适应度值最高的特征参量组合的适应度值;
[0030] 基于各监测系统的综合预警指数,根据下式计算出每一监测系统的预警效能RI,RI表示第I个监测系统的预警效能:
[0031]
[0032] 其中,为在监测时间段内产生预警并且报准的次数;N1为监测时间段内大能量或冲击事件的总次数;t0为产生预警占用时长;T0为监测总时长;当监测系统对应的综合预警指数在冲击或大能量事件发生前预设天数内超过预设阈值则预警为报准,否则为漏报或误报。
[0033] 进一步地,所述基于各监测系统的综合预警指数和预警效能,计算出煤岩动力灾害多系统综合预警结果,包括:
[0034] 基于各监测系统的综合预警指数和预警效能,根据下式计算出煤岩动力灾害多系统综合预警结果I综:
[0035]
[0036] 其中,W综I表示第I个监测系统的综合预警指数,RI表示第I个监测系统的预警效能。
[0037] 相应地,为解决上述问题,本发明还提供一种煤岩动力灾害多系统多参量集成综合预警系统,所述煤岩动力灾害多系统多参量集成综合预警系统包括:
[0038] 多元特征参量挖掘模,用于获取多个煤岩动力灾害监测系统的监测数据,并基于所述监测数据挖掘出各监测系统中能反映煤岩动力灾害前兆信息的多元特征参量;
[0039] 特征参量组合筛选模块,用于基于所述多元特征参量,筛选出每一监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值;
[0040] 单系统预警指数及预警效能计算模块,用于基于各监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值,计算出每一监测系统的综合预警指数和预警效能;
[0041] 综合预警结果计算模块,用于基于各监测系统的综合预警指数和预警效能,计算出煤岩动力灾害多系统综合预警结果。
[0042] 本发明的上述技术方案的有益效果如下:
[0043] 本发明通过对微震、地音、应力和电磁辐射等设备监测到的数据进行深入挖掘分析,得出各系统的预警特征参量,将各系统特征参量两两配对后输入遗传算法中得到该特征参量组合的最优临界值V1、V2和相应的适应度值,优选出适应度值最高的组合作为该监测系统的特征参量组合,之后计算出各单系统的预警程度Wi和综合预警指数W综,并利用R值评分法得出各单系统预警效能RI,最后计算出煤岩动力灾害多系统综合预警结果I综,确定当前煤岩动力灾害危险程度以及采取的措施。该方法将多监测系统优势互补,提高了预警准确率并剔除了主观因素的影响,可为矿山灾害防治工作的开展起到一定指导作用。附图说明
[0044] 图1为离散度时序变化曲线;
[0045] 图2为频次比时序变化曲线;
[0046] 图3为能量偏差值时序变化曲线;
[0047] 图4为频次偏差值时序变化曲线;
[0048] 图5为能量平均值(取对数)时序变化曲线;
[0049] 图6为脉冲因子时序变化曲线;
[0050] 图7为记录频次时序变化曲线;
[0051] 图8为闪烁间隔危险度时序变化曲线;
[0052] 图9为多系统综合预警结果;
[0053] 图10为本发明的煤岩动力灾害多系统多参量集成综合预警方法流程图
[0054] 图11为本发明的遗传算法计算最优临界值的流程示意图;
[0055] 图12为本发明采用的遗传算法在MATLAB 2014b环境中的运行过程图。

具体实施方式

[0056] 为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
[0057] 第一实施例
[0058] 请参阅图1至图12,本实施例提供一种煤岩动力灾害多系统多参量集成综合预警方法,如图10所示(图中“……”表示各监测系统中其它可实时在线读取的连续监测数据或离散监测数据),所述方法包括:
[0059] S101,获取多个煤岩动力灾害监测系统的监测数据,基于所述监测数据挖掘出各监测系统中能反映煤岩动力灾害前兆信息的多元特征参量;
[0060] 需要说明的是,上述监测系统为矿井已安装的监测系统,其包括在线监测系统和便携监测系统;其中,在线监测系统包括但不限于微震监测系统、地音监测系统和液压支架监测系统。上述监测数据包括在线监测系统实时采集得到的连续监测数据和便携监测系统采集得到的离散监测数据,进一步地,监测数据不仅包括能够直接读取的数据,而且还包括经数据挖掘后得到的深层次时间、空间、强度多元特征参量。
[0061] 多元特征参量为时间、空间、强度多元特征参量,其包括但不限于:
[0062] 微震监测系统中反映时序信息的频次、频次比和频次偏差值;反映空间信息的离散度;反映强度信息的微震能量、微震能量偏差值和离散度;
[0063] 地音监测系统中反映强度信息的地音能量、地音能量偏差值、地音能量平均值和脉冲因子;反映时序信息的脉冲;
[0064] 液压支架监测系统中反映时序信息的闪烁间隔危险度;反映强度信息的记录频次和支架压力。
[0065] 具体地,上述多元特征参量的计算方式如下:
[0066] 频次比P(t)的取值由下式计算:
[0067]
[0068] 其中,a(t)表示大能量事件发生前第n小时(可以取8小时)内的平均频次,b(T)表示大能量事件前n天(可以取2天)内的平均频次;
[0069] 频次偏差值DP(t)的取值由下式计算:
[0070]
[0071] 其中,N(t)表示该天发生微震事件的个数, 表示两次大能量矿震事件(可以取能量>106J)之间平均每日发生微震事件的频次;
[0072] 离散度Q(t)的取值由下式计算:
[0073]
[0074] 其中,TE(t)为一段时间内最大能量和中位数能量之和;LE(t)为微震临界特征能量,具有计算步骤如下:
[0075] 1)计算近一段时间内(可以取前8小时)发生的所有微震事件位置的加权平均距离d:
[0076]
[0077] 其中,前一段时间内共发生n次微震事件,x、y、z的平均值是指这n次微震事件的均值,对于第i次微震事件,其坐标为xi、yi、zi。wi为第i个微震事件能量的开方值;
[0078] 2)将历史微震事件计算的加权平均距离从小到大排序;
[0079] 3)找出其中位于a1、a2、a3、a4…ax位置所对应的加权平均距离的序号,然后向下取整,共得到x个加权平均距离对应的序号(其中x可以取7,对应的ai各取值可以分别为1.2%,7.04%,27.5%,80.84%,88.51%,94.84%,98.34%);
[0080] 4)将上述得到的每个加权平均距离左右半径一定范围内(可以取一米内)的所有特征能量点分别按从小到大的顺序排列;
[0081] 5)将每个加权平均距离周围一定范围内的特征能量点取出(该范围个数与3)中的x相同,对应范围可以取(40%-45%),(45%-50%),(98%-99%),(98%-99%),(98%-99%),(90%-98%),(85%-90%)内的所有值,将取出的特征能量点计算平均值后得到临界特征能量点,从而得到临界特征能量LE(t);
[0082] 6)将临界特征能量LE(t)对应的加权平均距离找出,即在加权平均距离和特征能量对应的坐标系中找出x个点,利用插值法得到临界特征能量曲线;
[0083] 7)将数据代入 中即可得到离散度;
[0084] 微震能量偏差值DE(t)的取值由下式计算:
[0085]
[0086] 其中,E(t)表示该t时刻发生微震事件的能量, 表示两次大能量矿震事件(可6
以取能量>10J)之间平均每次微震事件的能量;
[0087] 地音能量偏差值D(t)的取值由下式计算:
[0088]
[0089] 其中, P(t)为一段时间(可以取10分钟)内研究参量的累积和除以一段时间的时间窗长度,T(t)表示时间间隔,这一时间间隔与t时刻相关,是t时刻之前一段时间(可以取24小时)到t时刻之间的时间间隔,n为T(t)时间间隔内的地音加权能量的个数;
[0090] 地音能量平均值AE为给定时间(可以取6小时)内的平均能量;
[0091] 脉冲因子P得取值由下式计算:
[0092]
[0093] 其中,n是近一段时间(可以取6h)某一时间窗(可以取10min)的个数, 的计算方法为:将近一段时间(可以取6h)地音某一时间窗(可以取10min)能量累计值按从大到小排序,其是前n/10个数(向下取整)的均值, 是近一段时间(可以取6h)地音某一时间窗(可以取10min)能量累计值的均值;
[0094] 记录频次H为近一段时间(可以取6h)内的支架压力采集频次;
[0095] 闪烁间隔危险度的计算方法为: 其中t1表示前一次闪烁的时间,t2表示随后一次闪烁的时间。闪烁为将一段时间(可以取2天)内的平均压力值记为均值,后续压力数据若从高于均值,下降到均值的一定范围以下(可以为均值的1/10),并再次上升到均值以上时,记为一次闪烁。
[0096] S102,基于所述多元特征参量,筛选出每一监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值;
[0097] 需要说明的是,本实施例中,上述步骤具体是利用遗传算法训练并优选出每一监测系统中预警效能最高的特征参量组合及对应的最优临界值;
[0098] 其中,遗传算法是一种模仿自然界生物进化机制的在目标空间内进行有导向的过程搜索最优解算法,先天具有可进行多参数、多组合同时优化的特征,能直接对结构对象进行操作,不存在求导和函数连续性的限定。其实质是利用自然界生物“优胜劣汰”的机制,由选择、交叉、变异等遗传操作通过迭代的方式最终筛选出对环境适应度高的样本,并且它采用概率化的寻优方法,利用适应度函数的评价信息而不是传统导数的目标函数,形式对目标函数没有要求,因而有良好的适应性和可规模化特征。
[0099] 具体地,本实施例中利用遗传算法训练并优选出每一监测系统中预警效能最高的特征参量组合及对应的最优临界值的过程包括:
[0100] 1)将同一监测系统中属于不同维度的特征参量进行两两组合;也即组合时需将各监测系统中属于时间、空间、强度三种不同维度的特征参量进行组合;
[0101] 2)利用遗传算法训练并优选出每一监测系统中适应度值最高的特征参量组合并计算出其中各特征参量的最优临界值;
[0102] 3)将各监测系统中适应度值最高的特征参量组合作为该监测系统对应的预警效能最高的特征参量组合。
[0103] 进一步地,上述利用遗传算法训练并优选出每一监测系统中适应度值最高的特征参量组合并计算出最优临界值的过程如图11所示,包括:
[0104] (i)输入不同监测系统的特征参量组件;
[0105] (ii)根据步骤(i)输入的数据产生初始样本;
[0106] (iii)利用遗传算法的适应度函数计算适应度;
[0107] (iv)根据步骤(iii)中计算的适应度判断是否满足遗传算法中设定的终止条件;
[0108] (v)如果步骤(iv)的判断满足终止条件,则输出特征参数的最优临界值;如果步骤(iv)的判断不满足终止条件,则进入步骤(vi);
[0109] (vi)重新依次进行选择、交叉与变异步骤(ii)产生的初始样本,保留适应度值最高的个体同时淘汰适应度低的个体,产生新的样本,并返回步骤(iii)进行计算,直至满足所述终止条件。
[0110] 其中,上述步骤(iii)中,遗传算法的适应度函数采用R值评分法,R值评分法是20世纪70年代地球物理所提出的用于评价地震预报效能的方法,其思路清晰、便于操作;既不鼓励虚报,又不鼓励漏报,几十年来被广泛应用于地震预测预报中,取得了良好的效果。R值评分法公式如下式所示:
[0111]
[0112] 其中, 表示产生预警并且报对的次数;N1表示有大能量或冲击事件的总次数;表示产生预警但误报的次数;N0表示无大能量或冲击事件的总次数,其中,大能量指能量大于105J的微震事件,冲击指有冲击显现的事件;定义特征参量值在冲击或大能量事件发生前两天内超过最优临界值则为报准,否则为漏报或误报;R值越接近于1表示预警准确率越高,选取适应度值R最大的特征参量组合a和b及其对应的最优临界值V1和V2。
[0113] 遗传算法在MATLAB 2014b环境中的运行过程如图12所示,其中a为适应度函数在每次迭代过程中的惩罚值(为1-适应度值);b为迭代次数。
[0114] 此外,在利用遗传算法训练并优选出每一监测系统中适应度值最高的特征参量组合时,对于微震监测系统,由于其监测全采区范围的微震事件,所有探头综合采集一组监测结果,故需选取适应度值最高的前两组特征参量组合;
[0115] 对于地音监测系统和液压支架监测系统,由于其分别采集不同范围的监测信息,每个传感器都有一组监测结果,故选取每一传感器对应的适应度值最高的一组特征参量组合。
[0116] S103,基于各监测系统的预警效能最高的特征参量组合及其中各特征参量的最优临界值,计算出每一监测系统的综合预警指数和预警效能;
[0117] 需要说明的是,本实施例中,上述步骤具体包括:
[0118] 1)基于各监测系统的适应度值最高的特征参量组合及其中各特征参量的最优临界值,根据下式计算出各监测系统的单系统预警程度:
[0119]
[0120] 其中,对于微震监测系统,单系统预警程度为适应度值最高的每一特征参量组合的预警程度,此时Wi表示第i组适应度值最高的特征参量组合的预警程度,ai和bi分别表示第i组适应度值最高的特征参量组合中两个特征参量的实时值,V1和V2分别表示第i组适应度值最高的特征参量组合中各特征参量对应的最优临界值;
[0121] 对于地音监测系统和液压支架监测系统,单系统预警程度为各传感器的预警程度,此时Wi表示第i个传感器的预警程度,ai和bi分别表示第i个传感器对应的适应度值最高的特征参量组合中两个特征参量的实时值,V1和V2分别表示第i个传感器对应的适应度值最高的特征参量组合中各特征参量对应的最优临界值;
[0122] 2)基于各监测系统的单系统预警程度,根据下式计算出各监测系统的综合预警指数W综:
[0123]
[0124] 其中,对于微震监测系统,Ri为第i组适应度值最高的特征参量组合的适应度值;对于地音监测系统和液压支架监测系统,Ri为第i个传感器对应的适应度值最高的特征参量组合的适应度值;
[0125] 3)基于各监测系统的综合预警指数,利用R值评分法,根据下式计算出每一监测系统的预警效能RI,RI表示第I个监测系统的预警效能:
[0126]
[0127] 其中, 为在监测时间段内产生预警并且报准的次数;N1为监测时间段内大能量或冲击事件的总次数;t0为产生预警占用时长;T0为监测总时长;当监测系统对应的综合预警指数在冲击或大能量事件发生前预设天数内超过预设阈值则预警为报准,否则为漏报或误报。
[0128] 具体地,本实施例定义W综值在冲击或大能量事件发生前2天内超过0.5则为报准,否则为漏报或误报;每次产生预警的时长为2天。
[0129] S104,基于各监测系统的综合预警指数和预警效能,计算出煤岩动力灾害多系统综合预警结果。
[0130] 需要说明的是,在本实施例中,上述步骤具体为:
[0131] 根据下式计算出煤岩动力灾害多系统综合预警结果I综:
[0132]
[0133] 其中,W综I表示第I个监测系统的综合预警指数,RI表示第I个监测系统的预警效能。
[0134] 进一步地,在S104之后,本实施例的方法还包括:
[0135] S105,将煤岩动力灾害多系统综合预警结果与煤岩动力灾害分级预警评判准则相对比,确定动力灾害危险级别。
[0136] 需要说明的是,在本实施例中,上述步骤根据对比结果确定的危险级别为:当0≤I综≤0.25时,为无危险状态;当0.25<I综≤0.5时,为弱危险状态;当0.5<I综≤0.75时,为中危险状态;当0.75<I综≤1时,为强危险状态。
[0137] 根据本发明的预警结果,可以根据表1采取相应的措施。
[0138] 表1 煤岩动力灾害危险状态分级表
[0139]
[0140]
[0141] 下面,结合具体的应用场景,对本实施例的方法作进一步说明:
[0142] 选取某矿井的某一工作面,在工作面附近分别布置微震、地音和液压支架监测系统,预警步骤如下:
[0143] (1)收集并分析该矿井各煤岩动力灾害监测系统的历史监测数据,挖掘出各系统中可反映煤岩动力灾害前兆信息的时间、空间、强度多元特征参量,包括微震监测系统中离散度、频次比、能量偏差值和频次偏差值,地音监测系统中能量平均值和脉冲因子,支架监测系统中闪烁间隔危险度和记录频次。其中离散度为基于空间维度的特征参量,一段时期内微震事件震源位置离散程度越小,则之后诱发冲击事件的可能性越大;频次比和频次偏差值为时间维度特征参量,冲击事件发生前微震事件发生频次较之前降低,出现“平静期”现象,为后期更大破坏积蓄能量;能量平均值和能量偏差值为强度维度特征参量,地音或微震监测到的能量值越大,则煤岩体内部破裂现象越严重;脉冲因子为时间维度特征参量,表示地音能量产生相对峰值的大小;闪烁间隔危险度为时间维度特征参量,支架压力的快速变化称为闪烁,闪烁能排除移架造成的影响且其发生频率越大则煤岩体内部能量积蓄越严重;记录频次为强度维度特征参量,煤岩体内部应力变化时液压支架传感器的记录频次会随之升高。以两个月为时间窗口绘制如下各特征参量的时序变化曲线,结果如图1-图8所示。
[0144] (2)将各系统不同维度的特征参量两两组合输入遗传算法中计算最优临界值V1和V2及其对应的适应度值R,优选出微震监测系统的适应度值最高的两组特征参量组合,分别是离散度和频次比,能量偏差值和频次偏差值;地音监测系统的能量平均值(取对数)和脉冲因子;支架监测系统的记录频次和闪烁间隔危险度。计算结果如表2所示:
[0145] 表2 遗传算法计算结果
[0146]
[0147] (3)计算各单系统预警程度Wi和综合预警指数W综,并利用R值评分法得出各单系统预警效能RI。在历史监测时间段内该工作面共发生9次冲击或大能量事件,为便于表述将其编号为①-⑨;并将预警成败的依据定义为:冲击或大能量事件发生前3天内综合预警指数W综超过0.5则为预警成功,反之则失败。计算结果如下表3所示。
[0148] 表3 单系统综合预警指数W综和预警效能RI
[0149]
[0150] 在冲击或大能量事件发生前各系统综合预警指数W综基本都达到了预警标准,预警结果较准确,具有一定参考性。
[0151] (4)根据步骤(3)中表3的计算结果,计算出该矿煤岩动力灾害多系统综合预警结果I综,以两个月为绘图时间窗口,结果如图7所示:
[0152] (5)将多系统综合预警结果I综与煤岩动力灾害分级预警评判准则相对比以确定危险级别:当0≤I综≤0.25时,为无危险状态;当0.25<I综≤0.5时,为弱危险状态;当0.5<I综≤0.75时,为中危险状态;当0.75<I综≤1时,为强危险状态。
[0153] 本实施例的方法采用多系统综合预警结果作为预警指标对监测大范围区域进行煤岩动力灾害危险预警,实现了多系统多参量集成时空综合预警,且准确性高、操作简单、可实现编程操作,为煤岩动力灾害的高效预警奠定了基础
[0154] 第二实施例
[0155] 本实施例提供一种煤岩动力灾害多系统多参量集成综合预警系统,所述煤岩动力灾害多系统多参量集成综合预警系统包括:
[0156] 多元特征参量挖掘模块,用于获取多个煤岩动力灾害监测系统的监测数据,并基于所述监测数据挖掘出各监测系统中能反映煤岩动力灾害前兆信息的多元特征参量;
[0157] 特征参量组合筛选模块,用于基于所述多元特征参量,筛选出每一监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值;
[0158] 单系统预警指数及预警效能计算模块,用于基于各监测系统中预警效能最高的特征参量组合及其中各特征参量的最优临界值,计算出每一监测系统的综合预警指数和预警效能;
[0159] 综合预警结果计算模块,用于基于各监测系统的综合预警指数和预警效能,计算出煤岩动力灾害多系统综合预警结果。
[0160] 进一步地,本实施例的煤岩动力灾害多系统多参量集成综合预警系统还包括:动力灾害危险级别确定模块,用于将煤岩动力灾害多系统综合预警结果与煤岩动力灾害分级预警评判准则相对比,确定动力灾害危险级别。
[0161] 本实施例的煤岩动力灾害多系统多参量集成综合预警系统与上述第一实施例的煤岩动力灾害多系统多参量集成综合预警方法相对应;其中,该煤岩动力灾害多系统多参量集成综合预警系统中各模块的作用与上述煤岩动力灾害多系统多参量集成综合预警方法中的流程步骤相对应;故,在此不再赘述。
[0162] 此外,需要说明的是,本领域内的技术人员应明白,本发明实施例的实施例可提供为方法、装置、或计算机程序产品。因此,本发明实施例可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
[0163] 本发明实施例是参照根据本发明实施例的方法、终端设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、嵌入式处理机或其他可编程数据处理终端设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理终端设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0164] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理终端设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理终端设备上,使得在计算机或其他可编程终端设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程终端设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0165] 还需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者终端设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者终端设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者终端设备中还存在另外的相同要素。
[0166] 以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,一旦得知了本发明的基本创造性概念,在不脱离本发明所述原理的前提下,还可以对本发明的这些实施例做出若干变更和修改,这些变更和修改也应视为本发明的保护范围。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明实施例范围的所有变更和修改。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈