CN102608286B - 一种实时监测煤矿瓦斯监测值异常的方法 - Google Patents

一种实时监测煤矿瓦斯监测值异常的方法 Download PDF

Info

Publication number
CN102608286B
CN102608286B CN201210052314.7A CN201210052314A CN102608286B CN 102608286 B CN102608286 B CN 102608286B CN 201210052314 A CN201210052314 A CN 201210052314A CN 102608286 B CN102608286 B CN 102608286B
Authority
CN
China
Prior art keywords
gas concentration
gas
mrow
time
monitoring
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201210052314.7A
Other languages
English (en)
Other versions
CN102608286A (zh
Inventor
肖成
邸鹏浩
殷大发
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Coal Research Institute CCRI
Original Assignee
China Coal Research Institute CCRI
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Coal Research Institute CCRI filed Critical China Coal Research Institute CCRI
Priority to CN201210052314.7A priority Critical patent/CN102608286B/zh
Publication of CN102608286A publication Critical patent/CN102608286A/zh
Application granted granted Critical
Publication of CN102608286B publication Critical patent/CN102608286B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Sampling And Sample Adjustment (AREA)

Abstract

本发明公开一种实时监测煤矿瓦斯监测值异常的方法,包括:S1:在一段时间内按照预定周期重复采集回风巷固定区域瓦斯浓度值,获得历史数据样本;S2:获得用于反映在一个预定周期内的瓦斯浓度的分布情况的重排样本;S3:将所述重排样本按照时间划分成对应于连续的时间段间隔的多个子重排样本;S4:计算对应时间段间隔内的瓦斯浓度置信区间,利用计算得到所有时间段间隔的瓦斯浓度置信区间获得所述预定周期内的瓦斯浓度置信区间分布;S5:判断所述回风巷固定区域瓦斯浓度的实时监测值是否在监测时间所属的时间段间隔对应的瓦斯浓度置信区间内,若是,则判断煤矿瓦斯监测值正常,若否,则判断煤矿瓦斯监测值异常。

Description

一种实时监测煤矿瓦斯监测值异常的方法
技术领域
本发明属于煤矿瓦斯安全实时监测技术领域,尤其涉及一种实时监测煤矿瓦斯监测值异常的方法。
背景技术
煤矿安全监控系统是煤矿中广泛使用的一种监控设备,主要用于监测煤矿井下的各种环境和设备状态,当瓦斯超限时进行断电控制,发现其他异常时进行报警或断电控制。很多用于煤矿瓦斯的安全监测监控系统对于预防和降低煤矿安全生产事故起了很大的作用。但随着煤矿数量的增多、产量的增加和不规范的操作,安全事故还是时有发生。传统的安全监控系统是按照煤矿井下不同的区域来设定报警上下限的,并且按照一定的等级划分报警级别。一般规定在低瓦斯和高瓦斯矿井的采煤工作面瓦斯的报警浓度为>=1.0%CH4、断电浓度为>=1.5%CH4、复电浓度为<1.0%CH4,对于煤(岩)瓦斯突出矿井的巷道、工作面和其它工作区域报警上下限的值另有规定,但这些规定是对所有矿井的统一规定。
目前煤矿安全监控系统对(CH4、CO、风速等)模拟信号异常的判断方法主要是量值超过或低于一定值(也称为超限报警)。这些方法都是基于单个传感器以测量值的大小为判别依据,不同煤矿不同区域其瓦斯涌出规律不尽相同,但是却使用同样的报警规则,对于煤矿监控系统使用中的一些异常(如传感器零点漂移或催化元件老化或探头被异物堵塞等等)不能发现。
公开号为CN 101649751A中国专利公开一种煤矿安全监控系统的检测方法,该专利公开的“监测数据正常表现模型”虽然能发现一些不按规定使用的异常,但该模型不能全面反映煤矿的实际情况,建模型比较复杂,其实用性需要进一步提高。加上矿井的生产系统、巷道布置、通风系统、采煤工艺、采煤工序、班次制度、煤层特点等等都有可能成为影响瓦斯涌出波动的因素,因此想针对每种影响因素来制定不同的预警规则几乎不可能。
瓦斯涌出异常预报煤与瓦斯突出,秦汝祥,煤炭学报Vol.31,NO.5,599-602页,公开了一种瓦斯涌出异常预报煤与瓦斯突出,并具体公开了以下技术特征:通过对掘进工作面突出前夕煤壁瓦斯浓度变化规律的研究,提出了利用瓦斯浓度变化预报煤与瓦斯突出的方法,并建立了突出预报数学模型。研究表明:对瓦斯浓度序列的移动平均线、振幅、频次和方差等几个方面的综合评判,可预报掘进工作面突出的发生,并以潘一矿突出实例进行验证,表明在突出前215h瓦斯浓度变化出现异常,根据这一异常可实现突出预报。
发明内容
为了克服现有技术中的缺陷和不足,本发明提出一种实用性强、准确性高,可灵活扩展的实时监测煤矿瓦斯监测值异常的方法。该方法根据煤矿的实际情况,对特定区域内的瓦斯历史数据进行分析,得出这一区域内的瓦斯波动规律,并根据这一规律制定出符合本区域的自身报警预警方法。
本发明公开一种实时监测煤矿瓦斯监测值异常的方法,包括如下步骤:
S1:在一段时间内按照预定采样周期重复采集相对当前最近2个月内的回风巷固定区域瓦斯浓度值,获得历史数据样本;
S2:根据步骤S1中的历史数据样本,按照每天0到24点中相对时间的同一时刻的数据进行重新排列,得到一个在每天0到24点内相对时间的瓦斯浓度分布情况重排样本;
S3:将所述重排样本按照时间划分成对应于连续的时间段间隔的多个子重排样本;
S4:针对每个子重排样本按正态分布计算对应时间段间隔内的瓦斯浓度置信区间,利用计算得到所有时间段间隔的瓦斯浓度置信区间获得0到24点内的瓦斯浓度置信区间分布,所述瓦斯浓度置信区间表示在对应时间段间隔内测量位置的正常的瓦斯浓度范围;
S5:判断所述回风巷固定区域瓦斯浓度的实时监测值是否在监测时间所属的时间段间隔对应的瓦斯浓度置信区间内,若是,则判断煤矿瓦斯监测值正常,若否,则判断煤矿瓦斯监测值异常。
进一步地,所述步骤S3中多个子重排样本以5分钟作为时间段间隔划分得到。
进一步地,所述步骤S4中瓦斯浓度置信区间分布区间值能够根据置信水平P进行调整,其中,置信水平P为小于1的正数。
进一步地,所述步骤S4中计算瓦斯浓度置信区间的步骤包括:
S401:根据下述公式(1)计算瓦斯浓度置信区间所对应时间段间隔内数据样本的数学期望
x &OverBar; = 1 n &Sigma; i = 1 n x i    (式1)
其中,n为瓦斯浓度置信区间所对应时间段间隔内的瓦斯监测值个数;xi为所在对应时间段间隔内第i个瓦斯浓度检测值;
S402:根据下述公式(2)计算瓦斯浓度置信区间所对应时间段间隔内数据样本的标准差s,
s = 1 ( n - 1 ) &Sigma; i = 1 n ( x i - x &OverBar; ) 2    (式2);
S403:根据下述公式(3)计算瓦斯浓度置信区间的上下限值,得到瓦斯浓度置信区间,
( &theta; &OverBar; , &theta; &OverBar; ) = x &OverBar; &PlusMinus; z &alpha; / 2 &CenterDot; s n    (式3)
其中,为双侧置信区间,θ分别为双侧置信下限与双侧置信上限,zα/2为双侧置信区间的置信因子,根据显著性水平α,zα/2值通过查正态分布表得到。
进一步地,步骤S401中还包括通过选取数据样本中位数的方法得到对应时间段间隔内数据样本的数学期望。
进一步地,该方法能够根据回风巷瓦斯监测值异常判断回风巷上游工作面和上隅角的瓦斯信息异常。
本发明实时监测煤矿瓦斯监测值异常的方法克服了现有煤矿瓦斯监控系统的阈值报警方法单一、无法判断瓦斯异常等不足,根据同一时刻煤矿瓦斯浓度变化符合正态分布规律,选择回风巷固定区域相对当前最近一段时间的历史数据建立瓦斯浓度置信区间分布,相对当前最近一段时间历史数据与当前瓦斯传感器浓度监测值具有较强的相关性,能够较准确的反映瓦斯监测值情况,异常判断准确性高;同时,该方法可根据实际情况进行调整,可灵活扩展,具有较强的实用性。
附图说明
图1为本发明实时监测煤矿瓦斯监测值异常时传感器分布示意图;
图2为本发明实施例的一种实时监测煤矿瓦斯监测值异常的方法流程图;
图3为本发明方法中瓦斯浓度置信区间分布示意图;
图4为本发明方法中瓦斯浓度置信区间分布组成示意图。
具体实施方式
为详细说明本发明的技术内容、所实现目的及效果,以下结合实施方式并配合附图予以详细说明。
如图1所示,为采用本发明实时监测煤矿瓦斯监测值异常时传感器分布示意图,其中T0为上隅角、T1为工作面、T2为回风巷、T3为进风巷。工作面T1和上隅角T0的瓦斯浓度据随时间的推移分布极不规律,波动起伏很大,主要原因是工作面受采动影响,瓦斯释放不均匀。上隅角T0瓦斯随风流呈现涡流状态,浓度极不稳定。回风巷T2中的瓦斯已与空气充分混合,并且风流相对未定,浓度变化较均匀。所以回风巷T2数据相对平滑稳定,并一定程度上呈现周期性变化。
回风巷T2中的瓦斯大都来自工作面(自身也有少量释放,但涌出较均匀,波动较小)。在正常的工作面瓦斯涌出情况下(不包括瓦斯突出),工作面T1和上隅角T0发生瓦斯超限报警时,回风巷T2中的瓦斯浓度也会升高。据此可以得出,回风巷T2与上游的工作面T1和上隅角T0存在一定的波动关系,虽然有些滞后。当上游瓦斯异常涌出之前,位处下游不远的回风巷中瓦斯已经有所上升,据此根据回风巷的瓦斯波动情况可以反推出工作面和上隅角的瓦斯异动。本方法选择回风巷T2作为研究对象进行分析。在回风巷T2的固定区域放置传感器进行瓦斯浓度异常监测,根据回风巷T2的瓦斯浓度情况,在工作面T1和上隅角T0用传感器对瓦斯浓度进行辅助监测。
图2为本发明实施例的一种实时监测煤矿瓦斯监测值异常的方法流程图。该方法包括:
S1:在一段时间内按照预定周期重复采集回风巷固定区域瓦斯浓度值,获得历史数据样本。
预定周期可根据实际情况进行设定。煤矿的生产活动安排通常是以天(24小时)为周期进行班次轮循,本技术方案以天(24小时)为预定周期进行数据采集,选择相对当前最近2个月的历史数据作为数据样本进行分析,原因在于,时间太久的数据可能已经不符合当前的工作面环境,时间太短又不能充分表达当前的瓦斯涌出规律。相对当前最近2个月的历史数据具有可参考性。历史数据样本可以根据工作面的实际生产变化情况来选取,可以是1个月或3个月等,重点在于选取的数据样本能够更好的体现出当前情况下的瓦斯涌出情况。其中,历史数据样本需要通过过滤方法去除伪数据,如异常大数、传感器调校数据等。
S2:以数据样本采集在所述预定周期内的相对时间为参考重新排列所述历史数据样本,获得用于反映在一个预定周期内的瓦斯浓度的分布情况的重排样本;
将过滤筛选后的数据样本分布到一个二维坐标上,纵轴为瓦斯浓度,横轴为时间(预定周期)。根据以预定周期(24小时)为周期的历史数据样本,按照预定周期(24小时)中相对时间的同一时刻的数据进行重新排列,得到一个在预定周期内相对时间的瓦斯浓度分布情况重排样本。如回风巷T2监测点在相对当前近两个月的历史数据样本,按照每天24小时的相对时间的同一时刻瓦斯浓度值进行排列,形成以24小时为周期的从0点到24点的瓦斯浓度分布数据样本。
S3:将所述重排样本按照时间划分成对应于连续的时间段间隔的多个子重排样本。
所述时间段间隔可以根据实际情况进行选择,可以按照分钟进行划分或者按照工作班次等进行划分,如5分钟或10分钟等。
S4:针对每个子重排样本按正态分布计算对应时间段间隔内的瓦斯浓度置信区间,利用计算得到所有时间段间隔的瓦斯浓度置信区间获得所述预定周期内的瓦斯浓度置信区间分布,所述瓦斯浓度置信区间表示在对应时间段间隔内测量位置的正常的瓦斯浓度范围。
瓦斯涌出量和生产活动有很大关系,产量大时,瓦斯涌出量也大,产量小,瓦斯涌出量也随之减少。由于瓦斯涌出量不规律,对整个数据样本而言,不可能通过数学公式得到如图3所示的两条平滑的曲线作为瓦斯浓度置信区间分布的边界,因此,本发明通过对预定周期按照时间划分成多个时间段间隔来对数据进行分割,每个对应时间段间隔内的数据样本符合正态分布规律,对每个时间段间隔进行计算得到一个瓦斯浓度置信区间,再将所有瓦斯浓度置信区间组成瓦斯浓度置信区间分布。参见图4,为所述瓦斯浓度置信区间分布组成示意图。设以24小时为预定周期,将以24小时为周期的瓦斯浓度重排数据样本按照5分钟的时间间隔分成288个连续的子重排样本。按照此种分割方式,相对当前最近2个月的重排历史数据样本中,每5分钟内的样本数据量大该在300~~600个之间,根据每5分钟时间段间隔内的数据样本符合正态分布规律计算得到对应5分钟时间段间隔的瓦斯浓度置信区间。在24小时预定周期中,需要计算得到共288(24*60/5)个瓦斯浓度置信区间。
选取所述相对当前最近2个月历史数据样本重排后得到的反映在一个预定周期内的瓦斯浓度的分布情况的重排样本中第一个5分钟段(00:00~~00:05)内的重排数据样本计算得到一个瓦斯浓度置信区间。
设x1,x2,…,xn是取自总体X的一个样本,n即第一个5分钟内的全部监测值个数。
设P=1-α(0<α<1)为置信水平,在本发明的实施例中取P=0.99,则显著性水平α=1-P=0.01。P=0.99的置信水平只是一个常用的值,这个值可以根据实际使用情况和现场情况进行调整。例如在实际应用中有很多异常点落在了置信水平P=0.99的置信区间之内,则说明P的取值过大;如果很多正常点落在了置信水平P=0.99的置信区间之外,则说明P的取值过小。因此要根据具体的应用情况对置信水平P的取值做出适当调整。
将所述5分钟内数据样本的数学期望(平均值)作为较优估计量μ,既计算第一个5分钟的瓦斯浓度置信区间,步骤包括:
S401:根据如下公式(1)并利用样本数据计算得出瓦斯浓度置信区间对应的第一个5分钟内的数据样本的数学期望
x &OverBar; = 1 n &Sigma; i = 1 n x i    (式1)
其中,n为瓦斯浓度置信区间所对应时间段间隔内的瓦斯监测值个数;xi为所在对应时间段间隔内第i个瓦斯浓度检测值;
本发明实施例中,所述数据样本属于离散型的随机的瓦斯浓度值,对于离散型数据样本求数学期望也就是求数据样本的平均值(较优估计量)。此外,本实施例中的数据样本的数学期望也可以根据情况通过选取数据样本中位数的方法将得到。
S402:根据如下公式(2)计算得到第一个5分钟内的样本标准差s
s = 1 ( n - 1 ) &Sigma; i = 1 n ( x i - x &OverBar; ) 2    (式2);
S403:根据如下公式(3)计算得到第一个5分钟的瓦斯浓度置信区间,
( &theta; &OverBar; , &theta; &OverBar; ) = x &OverBar; &PlusMinus; z &alpha; / 2 &CenterDot; s n    (式3),
其中,设为双侧置信区间,分别称θ为双侧置信下限与双侧置信上限;设zα/2为双侧置信区间的置信因子,根据显著性水平α和置信水平P的关系(P=1-α),zα/2值通过查正态分布表得到。当取置信水平P=0.99时,显著性水平α=1-P=0.01,通过查表得到zα/2=2.58,从而求得第一个瓦斯浓度置信区间。
依照上述方法分别求得共计288个瓦斯浓度置信区间,此288个瓦斯浓度置信区间作为一个整体得到瓦斯浓度置信区间分布。此瓦斯浓度置信区间分布作为判断回风巷瓦斯监测值异常的依据。
上述方案除以5分钟作为时间段间隔进行分割外,还可根据具体情况选择其他时间段间隔或改为按班次分段取值计算瓦斯浓度置信区间。
根据具体的应用情况,瓦斯浓度置信区间分布可通过调整置信水平P的取值适当调整区间分布范围,以达到最优的判断结果。
S5:判断所述回风巷固定区域瓦斯浓度的实时监测值是否在监测时间所属的时间段间隔对应的瓦斯浓度置信区间内,如是,则判断煤矿瓦斯监测值正常,若否,则判断煤矿瓦斯监测值异常。
回风巷固定区域的传感器对瓦斯浓度值进行实时监测,将所述监测值与在监测时间所属的时间段间隔对应的瓦斯浓度置信区间进行匹配,如监测值在瓦斯浓度置信区间内,则判断煤矿瓦斯监测值正常,如瓦斯浓度置信区间外(越界),则判断煤矿瓦斯监测值异常。
参见图3,为本发明检测煤矿瓦斯信息量异常的方法的瓦斯浓度置信区间分布示意图。如图3所示,例如在6:02监测值为0.5,把这一点与瓦斯浓度置信区间分布进行匹配,得到一个A1点,这一点正好落在第73个瓦斯浓度置信区间中,由此判断,这是一个正常的监测值。在9:11分监测值为0.62,把这一值与瓦斯浓度置信区间分布进行匹配,得到点A2,这一点落在第110个置信区间之外,并且超越上界。虽然监测值没有达到国家规定的报警值1.0,但按照长时间的数据累计来看,此时的瓦斯涌出明显高于正常情况,属于小概率事件,由此判断,回风巷中的瓦斯非常规涌出,同时关注上隅角T0和工作面T1点,此时上隅角T0和工作面T1可能已经超限报警或即将发生超限报警。在17:11分监测值为0.05,把这一值与瓦斯浓度置信区间分布进行匹配,得到点A3,这一点落在第206个置信区间之外,并且超越该置信区间的下界。按照长时间的数据累计来看,此时的瓦斯涌出明显低于正常情况,由此判断,回风巷中的瓦斯量非常规涌出,有可能是停工停产导致工作面瓦斯涌量确实较低,也有可能是传感器出现异常(传感器零点漂移或催化元件老化或探头被异物堵塞等等)。如果是传感器发生异常需及时排查原因,采取措施,保证传感器能正常工作。
按照上述方法,选取某矿一个回风巷的监测数据(测点编号102A)进行具体分析。选取最近2个月(2011-10-1~2011-11-30)的数据量作为历史数据样本,数据样本如表1所示:
序号 测点编号 监测值 监测时间
1 102A 0.2 2011-10-01 00:00:01
2 102A 0.21 2011-10-01 00:01:41
3 102A 0.15 2011-10-01 00:01:56
4 102A 0.17 2011-10-01 00:02:11
...... ...... ...... ......
95159 102A 0.16 2011-11-30 23:58:01
95160 102A 0.20 2011-11-30 23:58:31
95161 102A 0.21 2011-11-30 23:59:21
95162 102A 0.20 2011-11-30 23:59:51
表1
从表1中选取每天的第一个5分钟时间段间隔内的数据作为第1个瓦斯浓度置信区间的计算样本。如表2所示:
序号 测点编号 监测值 监测时间
1 102A 0.2 2011-10-01 00:00:01
2 102A 0.21 2011-10-01 00:01:41
3 102A 0.35 2011-10-01 00:01:56
4 102A 0.37 2011-10-01 00:02:11
...... ...... ...... ......
312 102A 0.18 2011-11-30 00:03:01
313 102A 0.21 2011-11-30 00:03:31
314 102A 0.24 2011-11-30 00:04:21
315 102A 0.25 2011-11-30 00:04:51
表2
该瓦斯浓度置信区间,根据上述公式(1)计算得到第一个5分钟内的样本平均值其中n=315(每个置信区间的样本数可能不同)。通过上述公式(2)计算得到第一个5分钟内的样本标准差s=0.51。当置信水平P=0.99时,根据显著性水平α=1-P=0.01,查询正态分布表得到zα/2=2.58。根据上述公式(3)可求得第一个5分钟的置信区间
以此类推,计算得到全部288个瓦斯浓度置信区间如表3所示:
表3
取测点102A一条近期监测数据(见表4第1条数据),监测值为0.20,监测时间为2011-12-0110:20:32,计算可知这条数据对应表3中的第10×12+20/5+1=125个瓦斯浓度置信区间;已知则0.24-0.05<0.2<0.24+0.05,对比得知,此条数据的监测值0.2落在置信区间内,由此判断,此条数据无异常,此时瓦斯监测值正常。
取测点102A一条近期监测数据(见表4第128条数据),监测值为0.39,时间为2011-12-01 16:20:32,计算可知这条数据对应表3中第16×12+27/5+1=198个瓦斯浓度置信区间;已知则0.23-0.07<0.23+0.07<0.39,对比得知,此条数据的监测值0.39落在置信区间之外并高于上限,由此判断,此条数据有异常(虽然没有超限,但超过了平时正常的浓度)。通过查找102A对应的工作面瓦斯测点105A(表5)相应时间内的数据,发现在表5中相应的时间点(与表4中第128条数据时间相近)第130条数据的监测值已经有所上升,并在1分钟后105A测点发生超限报警。
取测点102A一条近期监测数据(见表4第3893条数据),监测值为0.01,时间为2011-12-29 23:55:12,计算可知这条数据对应表3中第23×12+55/5+1=288个瓦斯浓度置信区间,已知则0.01<0.24-0.1<0.24+0.1,对比得知,此条数据的监测值0.01落在置信区间之外并低于下限,由此判断,此条数据有异常。经过检查,发现102A测点传感器被粉尘堵住了探头气孔,清理后传感器正常。
序号 测点编号 监测值 监测时间
1 102A 0.20 2011-12-01 10:20:32
...... ...... ...... ......
128 102A 0.39 2011-12-01 16:20:32
129 102A 0.40 2011-12-01 16:21:01
...... ...... ...... ......
3893 102A 0.01 2011-12-29 23:55:12
...... ...... ...... ......
表4
序号 测点编号 监测值 监测时间
1 105A 0.40 2011-12-01 10:20:11
...... ...... ...... ......
130 105A 0.65 2011-12-01 16:19:50
131 105A 0.97 2011-12-01 16:20:28
132 105A 1.17 2011-12-01 16:21:31
...... ...... ...... ......
表5
本发明实时监测煤矿瓦斯监测值异常的方法根据煤矿瓦斯浓度随时间周期性变化,不同周期的同一时刻监测值符合正态分布规律,选择回风巷固定区域的相对当前最近一段时间的历史数据得到瓦斯浓度置信区间分布,通过判断当前瓦斯传感器监测值与该监测时间所属的时间段间隔对应的瓦斯浓度置信区间匹配情况,判断煤矿瓦斯监测值异常情况,该方法异常判断准确性高,同时,该方法可根据实际情况进行调整,可灵活扩展,具有较强的实用性。
上述仅为本发明的较佳实施例及所运用技术原理,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围内。

Claims (6)

1.一种实时监测煤矿瓦斯监测值异常的方法,包括如下步骤:
S1:按照预定采样周期重复采集相对当前最近2个月内的回风巷固定区域瓦斯浓度值,获得历史数据样本;
S2:根据步骤S1中的历史数据样本,按照每天0到24点中相对时间的同一时刻的数据进行重新排列,得到一个在每天0到24点内相对时间的瓦斯浓度分布情况重排样本;
S3:将所述重排样本按照时间划分成对应于连续的时间段间隔的多个子重排样本;
S4:针对每个子重排样本按正态分布计算对应时间段间隔内的瓦斯浓度置信区间,利用计算得到所有时间段间隔的瓦斯浓度置信区间获得0到24点内的瓦斯浓度置信区间分布,所述瓦斯浓度置信区间表示在对应时间段间隔内测量位置的正常的瓦斯浓度范围;
S5:判断所述回风巷固定区域瓦斯浓度的实时监测值是否在监测时间所属的时间段间隔对应的瓦斯浓度置信区间内,若是,则判断煤矿瓦斯监测值正常,若否,则判断煤矿瓦斯监测值异常。
2.根据权利要求1所述的实时监测煤矿瓦斯监测值异常的方法,其特征在于,所述步骤S3中多个子重排样本以5分钟作为时间段间隔划分得到。
3.根据权利要求1所述的实时监测煤矿瓦斯监测值异常的方法,其特征在于,所述步骤S4中瓦斯浓度置信区间分布区间值能够根据置信水平P进行调整,其中,置信水平P为小于1的正数。
4.根据权利要求1所述的实时监测煤矿瓦斯监测值异常的方法,其特征在于,所述步骤S4中计算瓦斯浓度置信区间的步骤包括:
S401:根据下述公式(1)计算瓦斯浓度置信区间所对应时间段间隔内数据样本的数学期望
x &OverBar; = 1 n &Sigma; i = 1 n x i     (式1)
其中,n为瓦斯浓度置信区间所对应时间段间隔内的瓦斯监测值个数;xi为所在对应时间段间隔内第i个瓦斯浓度监测值;
S402:根据下述公式(2)计算瓦斯浓度置信区间所对应时间段间隔内数据样本的标准差s,
s = 1 ( n - 1 ) &Sigma; i = 1 n ( x i - x &OverBar; ) 2     (式2);
S403:根据下述公式(3)计算瓦斯浓度置信区间的上下限值,得到瓦斯浓度置信区间,
( &theta; &OverBar; , &theta; &OverBar; ) = x &OverBar; &PlusMinus; z &alpha; / 2 &CenterDot; s n     (式3)
其中,(θ,)为双侧置信区间,θ分别为双侧置信下限与双侧置信上限,zα/2为双侧置信区间的置信因子,根据显著性水平α,zα/2值通过查正态分布表得到。
5.根据权利要求4所述的实时监测煤矿瓦斯监测值异常的方法,其特征在于,所述步骤S401中还包括通过选取数据样本中位数的方法得到对应时间段间隔内数据样本的数学期望。
6.根据权利要求1所述的实时监测煤矿瓦斯监测值异常的方法,其特征在于,该方法能够根据回风巷瓦斯监测值异常判断回风巷上游工作面和上隅角的瓦斯信息异常。
CN201210052314.7A 2012-03-01 2012-03-01 一种实时监测煤矿瓦斯监测值异常的方法 Active CN102608286B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210052314.7A CN102608286B (zh) 2012-03-01 2012-03-01 一种实时监测煤矿瓦斯监测值异常的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210052314.7A CN102608286B (zh) 2012-03-01 2012-03-01 一种实时监测煤矿瓦斯监测值异常的方法

Publications (2)

Publication Number Publication Date
CN102608286A CN102608286A (zh) 2012-07-25
CN102608286B true CN102608286B (zh) 2015-07-08

Family

ID=46525826

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210052314.7A Active CN102608286B (zh) 2012-03-01 2012-03-01 一种实时监测煤矿瓦斯监测值异常的方法

Country Status (1)

Country Link
CN (1) CN102608286B (zh)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103984938B (zh) * 2014-05-30 2017-05-24 中国科学院遥感与数字地球研究所 一种遥感时间序列异常检测方法
CN104537199B (zh) * 2014-10-20 2017-09-05 天地(常州)自动化股份有限公司 利用瓦斯浓度时间序列识别矿井工作面生产工序的方法
CN104950285B (zh) * 2015-06-02 2017-08-25 西安理工大学 一种基于近邻标签信号差值变化的rfid室内定位方法
CN105181017A (zh) * 2015-09-11 2015-12-23 中国矿业大学 一种煤与瓦斯突出检测方法
CN105158610B (zh) * 2015-09-14 2018-12-11 广西电网有限责任公司电力科学研究院 一种变压器状态预警数据可疑值的筛查处理方法
CN106952106A (zh) * 2016-01-06 2017-07-14 阿里巴巴集团控股有限公司 一种数据的监控方法及装置
CN105718748B (zh) * 2016-01-28 2018-08-03 许继集团有限公司 一种在线监测数据处理方法
CN105653879B (zh) * 2016-02-03 2018-03-30 中煤科工集团重庆研究院有限公司 煤矿井下安全监控数据综合滤噪方法
CN106096633A (zh) * 2016-06-05 2016-11-09 丁旭秋 煤矿瓦斯浓度测量方法
CN106295214B (zh) * 2016-08-18 2017-11-10 西安科技大学 一种矿井瓦斯浓度预警方法
CN106894841B (zh) * 2017-04-28 2019-01-15 华北科技学院 基于正态分布检验辨识瓦斯涌出异常的瓦斯灾害预警方法
CN107843811B (zh) * 2017-11-02 2019-12-06 广东电网有限责任公司中山供电局 一种电网设备在线监测数据的分析方法及系统
CN109752504B (zh) * 2019-01-25 2021-11-30 西安科技大学 一种工作面瓦斯传感器调校辅助决策方法
CN110593954B (zh) * 2019-08-28 2020-06-19 精英数智科技股份有限公司 一种煤矿采煤工作面瓦斯传感器位置异常识别方法及装置
CN110631624B (zh) * 2019-09-04 2020-12-15 精英数智科技股份有限公司 一种识别矿井传感器运行数据异常的方法、装置及系统
CN110763809B (zh) * 2019-11-15 2022-03-29 中国石油大学(华东) 一种气体探测器优化布置方案的实验验证方法
CN112649034B (zh) * 2020-12-28 2023-09-19 精英数智科技股份有限公司 一种煤矿井下传感器调校状态的识别方法及装置
CN112763662B (zh) * 2020-12-29 2023-08-22 精英数智科技股份有限公司 气体传感器数据异常的识别方法及相关装置
CN116877203B (zh) * 2023-08-23 2024-03-12 河南理工大学 一种煤与瓦斯突出监测预警装置

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU626405A1 (ru) * 1977-04-14 1978-09-30 Государственный Макеевский Ордена Октябрьской Революции Научно-Исследовательский Институт По Безопастности Работ В Горной Промышленности Министерства Угольной Промышленности Ссср Способ контрол горючих газов
DE3327154C2 (zh) * 1983-07-28 1988-06-09 Hans-Joerg Dipl.-Kfm. 4400 Muenster De Huebner
JP4066921B2 (ja) * 2003-09-17 2008-03-26 日立電線株式会社 多点光路切替式ガス濃度検知方法及びその装置
CN101532397A (zh) * 2009-04-10 2009-09-16 重庆梅安森科技发展有限责任公司 煤与瓦斯突出实时诊断方法
CN101776596A (zh) * 2010-02-03 2010-07-14 中北大学 瓦斯浓度智能测试系统及方法
CN101975092A (zh) * 2010-11-05 2011-02-16 中北大学 基于径向基函数神经网络集成的矿井瓦斯浓度中短期实时预测方法
CN102109842A (zh) * 2010-12-15 2011-06-29 辽宁工程技术大学 矿用智能瓦斯安全信息监控装置及方法
CN102155264A (zh) * 2011-03-29 2011-08-17 天地(常州)自动化股份有限公司 一种检测煤矿监控系统中的瓦斯量监测信息异常的方法
CN102192894A (zh) * 2011-03-24 2011-09-21 周孟然 光纤传感煤矿瓦斯的激光光谱吸收检测方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU626405A1 (ru) * 1977-04-14 1978-09-30 Государственный Макеевский Ордена Октябрьской Революции Научно-Исследовательский Институт По Безопастности Работ В Горной Промышленности Министерства Угольной Промышленности Ссср Способ контрол горючих газов
DE3327154C2 (zh) * 1983-07-28 1988-06-09 Hans-Joerg Dipl.-Kfm. 4400 Muenster De Huebner
JP4066921B2 (ja) * 2003-09-17 2008-03-26 日立電線株式会社 多点光路切替式ガス濃度検知方法及びその装置
CN101532397A (zh) * 2009-04-10 2009-09-16 重庆梅安森科技发展有限责任公司 煤与瓦斯突出实时诊断方法
CN101776596A (zh) * 2010-02-03 2010-07-14 中北大学 瓦斯浓度智能测试系统及方法
CN101975092A (zh) * 2010-11-05 2011-02-16 中北大学 基于径向基函数神经网络集成的矿井瓦斯浓度中短期实时预测方法
CN102109842A (zh) * 2010-12-15 2011-06-29 辽宁工程技术大学 矿用智能瓦斯安全信息监控装置及方法
CN102192894A (zh) * 2011-03-24 2011-09-21 周孟然 光纤传感煤矿瓦斯的激光光谱吸收检测方法
CN102155264A (zh) * 2011-03-29 2011-08-17 天地(常州)自动化股份有限公司 一种检测煤矿监控系统中的瓦斯量监测信息异常的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于数据挖掘法的矿井瓦斯联动监测;牛立东;《中国安全科学学报》;20110731;第21卷(第7期);62-68 *
瓦斯涌出异常预报煤与瓦斯突出;秦汝祥;《煤炭学报》;20061031;第31卷(第5期);599-602 *

Also Published As

Publication number Publication date
CN102608286A (zh) 2012-07-25

Similar Documents

Publication Publication Date Title
CN102608286B (zh) 一种实时监测煤矿瓦斯监测值异常的方法
CN112161648B (zh) 异常传感器的识别方法、装置、可读存储介质及电子设备
AU2012200296B2 (en) System and method for identifying likely geographical locations of anomalies in a water utility network
US10554077B2 (en) Automated monitoring for changes in energy consumption patterns
CN110593954B (zh) 一种煤矿采煤工作面瓦斯传感器位置异常识别方法及装置
CN102155264B (zh) 一种检测煤矿监控系统中的瓦斯量监测信息异常的方法
CN103147799B (zh) 一种基于瓦斯涌出特征的突出危险性预警方法
CN114417697A (zh) 一种基于神经网络的tbm滚刀磨损实时预测方法及系统
CN104535735B (zh) 水质远程在线监测方法
CN117268455B (zh) 一种工程建筑施工质量检测设备的监测系统
CN104931989A (zh) 运动轨迹中异常点的检测方法与装置
CN106324179A (zh) 一种催化燃烧传感器的自动校准方法
EP3289370B1 (en) Methods and systems for alerting a user to the presence of a fault in an electromechanical system in a railway infrastructure
CN116105802A (zh) 一种基于物联网的地下设施安全监测预警方法
CN103277137A (zh) 掘进煤巷瓦斯突出预测方法
US20220137022A1 (en) Estimation Method
CN102155266B (zh) 检测煤矿监控系统中的瓦斯量监测信息异常的方法
CN109808641A (zh) 防盗告警方法及装置、电子设备和可读存储介质
CN111709606A (zh) 基于大数据分析的采掘工作面瓦斯涌出异常预警方法
CN111830208A (zh) 一种甲烷传感器位置状态监测方法及系统
CN103487275A (zh) 基于二维投射的煤矿设备状态识别和预警方法
CN116227781A (zh) 一种区域性电力能源结构合理性分析方法
CN101923605A (zh) 铁路防灾风预警方法
CN107809764A (zh) 一种基于Markov链的多事件探测方法
CN118038627B (zh) 基于车辆大数据的隧道路段安全监控方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: CHINA COAL RESEARCH INSTITUTE CO., LTD.

Free format text: FORMER OWNER: COAL GENERAL ACADEMY

Effective date: 20140801

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20140801

Address after: 100013 Beijing, Hepingli Youth ditch East Road, No. 5, No.

Applicant after: Co., Ltd of coal science and technology academy

Address before: 100013 Beijing, Hepingli Youth ditch East Road, No. 5, No.

Applicant before: General Research Institute of Coal Science

C14 Grant of patent or utility model
GR01 Patent grant