CN109325470B - 基于瓦斯浓度参数的井下工作面作业类型智能识别方法 - Google Patents
基于瓦斯浓度参数的井下工作面作业类型智能识别方法 Download PDFInfo
- Publication number
- CN109325470B CN109325470B CN201811243460.1A CN201811243460A CN109325470B CN 109325470 B CN109325470 B CN 109325470B CN 201811243460 A CN201811243460 A CN 201811243460A CN 109325470 B CN109325470 B CN 109325470B
- Authority
- CN
- China
- Prior art keywords
- gas concentration
- signal sequence
- local
- sequence
- working face
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000012545 processing Methods 0.000 claims abstract description 23
- 238000013528 artificial neural network Methods 0.000 claims abstract description 17
- 230000009467 reduction Effects 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 230000003044 adaptive effect Effects 0.000 claims abstract description 5
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 47
- 230000006870 function Effects 0.000 claims description 30
- 239000003245 coal Substances 0.000 claims description 22
- 238000009499 grossing Methods 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 10
- 238000012549 training Methods 0.000 claims description 7
- 238000005553 drilling Methods 0.000 claims description 6
- 238000012804 iterative process Methods 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000005422 blasting Methods 0.000 claims description 3
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 238000005520 cutting process Methods 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 3
- 230000000739 chaotic effect Effects 0.000 claims description 2
- 238000012423 maintenance Methods 0.000 claims description 2
- 230000000717 retained effect Effects 0.000 claims description 2
- 210000002569 neuron Anatomy 0.000 claims 3
- 238000012544 monitoring process Methods 0.000 abstract description 18
- 238000005065 mining Methods 0.000 abstract description 12
- 238000004458 analytical method Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 4
- 238000004519 manufacturing process Methods 0.000 description 7
- 230000002265 prevention Effects 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000011282 treatment Methods 0.000 description 5
- 238000007726 management method Methods 0.000 description 4
- 238000000926 separation method Methods 0.000 description 4
- 230000005641 tunneling Effects 0.000 description 4
- 238000009412 basement excavation Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000001537 neural effect Effects 0.000 description 3
- 238000003062 neural network model Methods 0.000 description 3
- 206010063385 Intellectualisation Diseases 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 239000002817 coal dust Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000009131 signaling function Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Signal Processing (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
本发明公开了一种基于瓦斯浓度参数的井下工作面作业类型智能识别方法。方法首先从矿井安全监测系统中获取工作面瓦斯浓度数据，并存入瓦斯参数数据库，对数据库中的瓦斯浓度数据进行自适应降噪滤波处理，而后在作业类型划分和敏感特征参数提取的基础上建立历史瓦斯浓度数据样本集，构建自适应模块化粗神经网络预测模型，利用构建的预测模型对其对应的作业类型进行辨识。本发明方法降噪效果显著，在有效去除噪声的同时能保留获得瓦斯浓度时间序列中的有效信息，能实现工作面作业类型的分析和识别，能够较好的实现瓦斯浓度数据的深度挖掘和二次利用。
Description
技术领域
本发明涉及矿井瓦斯灾害监测防治领域，尤其涉及一种基于瓦斯浓度数据特征的矿井采掘工作面作业类型智能分析与识别方法。
背景技术
中国是世界第一煤炭大国，根据我国《能源中长期发展规划纲要(2004-2020年)》明确指出在未来一段时期内我国仍将坚持以煤炭为主体能源。目前，我国煤炭开采业已进入深部开采时期，且开采深度正逐年增加。在深部开采时期，矿井瓦斯压力增大、地温增高、地应力增大，矿井灾害日益严重。煤矿事故灾害严重妨害着我国煤炭安全生产形势，阻碍了煤炭工业的健康、持续、稳定发展。健全煤矿事故灾害防治安全管理策略，深化煤矿事故灾害防控意识理念，提升煤矿事故灾害监管防控技术，对于确保煤矿安全生产、促进国民经济全面、健康发展具有重要意义。
煤矿生产中的瓦斯事故、突出、冲击地压、煤尘事故、水灾、火灾、顶板事故等主要灾害绝大多数发生在采掘工作面，与采掘、打钻等作业关系密切，且以掘进作业期间事故居多。矿井安全监测是矿井事故灾害防控的重要手段，及时获取和有效利用井下工作面的有效信息，有助于提升事故预防和灾害控制的工作质量。当前，国内外对工作面作业环境监测监控技术研究的主要侧重于瓦斯浓度、风速、温度等作业环境参数的监测。高瓦斯煤巷掘进工作面，工作面的瓦斯状况是与作业状况密切相关的环境参数，监测监控系统的长期运行使得煤矿现场积累了大量的瓦斯监测数据。国内外研究人员对瓦斯监测数据的利用早已有所关注，并用数据挖掘、人工智能等先进技术手段进行了大量的研究，如利用瓦斯参数序列进行瓦斯浓度预测、利用瓦斯涌出量来判读突出危险性等，这些均对矿井安全监测工作起到了巨大的促进作用。
然而，前人对工作面前方作业情况的直接关注较少，井下的作业施工情况主要来源于工作调度安排计划或者人为的统计记录，难免存在错误和纰漏，并且在时效性上受到了极大的限制。目前，作为事故灾害防控的常态化手段，矿井安全监测现阶段已经由单纯的安全监测面向灾害诊断、作业调度和生产管理等功能多元化、综合化、智能化方向发展并且呈现出非接触、动态化、多参量综合化、模块集成化的特点。采动作业的影响往往是很多事故的直接原因，通过井下的相关信息主动、精准、快速的辨识工作面的作业工序，并纳入系统数据库，对于事故机理分析、灾害诊断、应急处置和事故调查等工作均具有重大的实际意义，将为煤矿安全监测预警技术的发展提供基础信息支持。
发明内容
针对现有技术应用的缺失，本发明提供了一种基于矿井瓦斯浓度参数的井下工作面作业类型智能识别方法。该方法特殊针对瓦斯监测的浓度数据进行预处理，有效降低原始时间序列的噪声信号，保留和凸显有效瓦斯浓度序列特征，在次基础上通过构建瓦斯参数敏感参数体系和作业类型样本训练库，采用误差逆向传播算法训练的多层前馈神经网络，对井下工作面的作业工序类型进行智能辨识。
本发明的技术方案具体包括如下主要步骤：
步骤1、通过井下瓦斯浓度传感器获取瓦斯浓度的原始信号序列x(t)，作为历史数据存入瓦斯浓度数据库；
步骤2、根据矿井生产的实际情况，将各个原始信号序列x(t)根据工作面的作业类型进行划分；
步骤3、对瓦斯浓度数据库中的瓦斯浓度数据进行基于局域均值分解的特殊自适应滤波降噪处理，获得降噪后的瓦斯浓度序列；
步骤4、将降噪后的瓦斯浓度序列作为混沌时间序列，根据作业类型的划分对各类作业类型对应的数据分别提取敏感特征参数，构建作业类型划分敏感参数体系；
步骤5、构建井下工作面的作业类型神经网络智能识别模型，根据已建立的瓦斯浓度数据库和敏感特征参数输入到模型中进行训练；
步骤6、利用训练后的模型对新获取的待测瓦斯浓度数据进行作业类型的智能识别，获得待测瓦斯浓度数据对应的作业类型结果。
本发明创造性地利用瓦斯浓度参数来判断井下工作面作业类型，发挥了数据的作用，其中自适应滤波降噪处理作出了特殊的降噪处理使得看似与作业类型不相关的瓦斯浓度参数能体现表征井下工作面作业类型。
所述步骤2中，工作面的作业类型划分为如放炮、割煤、落煤、交接班、打钻作业、临时支护、安装锚杆、检修和其他无关作业。也可以添加其他作业形式或者有规律的组合其中若干个作业类型，其主要依据根据安全监测预警实际需求和作业类型辨识的可行性而确定的。
所述步骤4中，所述的敏感特征参数包括瓦斯浓度数据的均值、均方差、方差、极差、差分特征值等。具体也可以选用特定的乘积分量作为敏感特征参数。具体实施中，或者此类参数基础上的改良参数或其他二次挖掘参数；也可以使用在数据降噪预处理过程中产生的包括乘积分量序列等在内的其他敏感特征参数。
所述步骤3中，自适应滤波降噪处理具体为：
3.1，获得原始信号序列x(t)的全部极值点ni，i＝1,2,3…,I，I为极值点数，然后根据极值点，求得所有相邻局域极值点的平均值作为局域均值mi＝(ni+ni+1)/2和局域包络估计值作为局域包络值ai＝(ni-ni+1)/2；
3.2，对所有局域均值mi和所有局域包络值ai采用滑动平均法进行多次平滑处理分别得到局域均值函数m′1(t)、m′2(t)、…和局域包络函数a′1(t)、a′2(t)、…，下标1、2表示第一次、第二次平滑处理所得出的函数结果；
3.3，步骤3.2的每一次平滑处理后均处理迭代直至获得纯调频信号sk(t)：
在步骤3.2的第一次平滑处理后，从原始信号序列x(t)中分离出第一次平滑处理获得的局域均值函数m′1(t)而保留得到第一中间分离信号序列h1(t)＝x(t)-m′1(t)，然后再根据第一次平滑处理获得的局域包络函数a′1(t)进行解调得第一中间解调信号序列s1(t)＝h1(t)/a′1(t)；
然后将第一中间解调信号序列s1(t)视为原始信号序列x(t)重复步骤，结合第二次平滑处理获得的局域均值函数m′1(t)和局域包络函数a′1(t)计算获得第二中间分离信号序列和第二中间解调信号序列；
以此迭代处理直到第K次的中间解调信号序列满足-1≤sK(t)≤1且其下一次的局域包络函数a′K+1(t)＝1，K为进行迭代的次数，则将第K次的中间解调信号序列sK(t)作为纯调频信号，停止迭代处理；具体迭代过程表示为：
式中，K为进行迭代的总次数；
3.4，将以上纯调频信号sK(t)在迭代过程产生的所有局域包络函数a′k(t)相乘获得第一乘积分量(Production function，PF)的包络信号A1(t)＝a′1(t)a′2(t)…a′K(t)，k＝1,2,3…，K，接着将第一乘积分量的包络信号A1(t)和其对应的纯调频信号sK(t)相乘得到原始信号序列的第一乘积分量PF1(t)＝A1(t)sk(t)；
3.5，将第一乘积分量PF1(t)从原始信号序列x(t)中分离出来，保留得到第一残余信号u1(t)，并以第一残余信号u1(t)为原始信号序列x(t)重复上述步骤3.1-3.4迭代处理，直至通过J次循环分离出所有乘积分量，且满足第J残余信号uJ(t)为一个单调函数，由此构建以下原始信号序列x(t)的分解表达式：
其中，j表示步骤3.1-3.4迭代处理的次序数，j＝1,2,3…，J，J表示步骤3.1-3.4迭代处理的总次数；
原始瓦斯浓度信号变换为一组乘积分量和一个残余分量的和。
3.6，对获得的各个乘积分量和原始信号序列x(t)一起进行相关性分析，取所有乘积分量的能量值构成序列的全局极小值作为高低频成分的调频分界点，将调频分界点之前的乘积分量的高频分量全部作为噪声乘积分量，将剔除噪声乘积分量后对信号进行重构处理，得到降噪后的瓦斯浓度序列dx(t)。
对信号进行重构处理是指将各个乘积分量进行加合。
所述的作业类型神经网络智能识别模型采用粗糙集理论和人工神经网络相结合的粗神经网络(rough membership neural network，RMNN)，以实现对井下采掘工作面工作类型的智能识别，利用神经网络输入层的粗神经元分别计算各作业类型时瓦斯浓度数据敏感特征参数，并采用模糊神经元代替普通神经元，采用部分历史瓦斯浓度数据的敏感特征参数作为训练样本，完成对新近瓦斯数据对应的作业类型进行智能识别。
本发明针对当前矿井瓦斯浓度数据中通常含有大量噪声影响预测模型预测精度的问题，通过局域均值分解将历史瓦斯浓度序列数据分解成多个纯调频信号分量，再将分解后的调频信号函数经过低通滤波器自适应阈值降噪滤波，在有效去除噪声的同时，能够最大限度的保留瓦斯浓度时间序列中的有效信息，除噪效果理想。
本发明的有益效果是：
本发明针对人工记录作业工序及类型的时效延迟性和不准确性等缺陷，构建一个自适应模块化的粗神经网络预测模型，模型可以基于多种瓦斯浓度参数的统计参数或其构成的综合参数根据瓦斯浓度序列，用于实现工作面作业类型的分析和识别，能够较好的实现瓦斯浓度数据的深度挖掘和二次利用。
本发明方法降噪效果显著，在有效去除噪声的同时能保留获得瓦斯浓度时间序列中的有效信息。
本发明方法可实现对井下实际工作类型的精准判识，符合当前矿井安全监测系统的发展要求。
附图说明
图1为本发明实施例方法流程图；
图2为本发明实施例提供的原始瓦斯浓度数据时间序列图；
图3为本发明实施例提供的降噪后的瓦斯浓度数据时间序列图；
图4为本发明实施例提供的基于综合参数的粗神经网络模型作业类型分类图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
本发明实施例具体如下：
通过瓦斯浓度传感器采集某煤矿采掘工作面的瓦斯浓度数据，作为实施例研究对象，同时人工进行周期段内的作业类型跟踪和记录，作为样本库，具体实施流程如图1所示，具体方法如下。
步骤1、通过矿井安全监控系统自动连续获取掘进面瓦斯涌出动态时间序列数据，本实施例中矿井采用KJ90(分站-F16)安全监控系统获取的瓦斯数据，其瓦斯浓度传感器为KG9701A甲烷传感器，瓦斯浓度数据采集频率为每5分钟一个数据。记井下瓦斯浓度传感器获取原始瓦斯浓度原始序列x(t)，作为历史数据存入数据库，任意截取某时段的瓦斯浓度信号如图2所示；
步骤2、根据矿井生产的实际情况，对工作面的作业类型进行合理划分。本实施例中矿井作业方式采用“四六制”，整个作业周期放炮、割煤、落煤、交接班、打钻作业、临时支护、安装锚杆、检修和其他无关作业等多种作业形式，为便于阐述，本实施例中仅根据采动作业的强烈程度将作业类型划分为“无作业”、“打钻作业”和“掘进作业”三种。
步骤3、根据对作业类型的人工记录，在相应时段内分别筛选n组瓦斯浓度序列样本，并对数据进行基于局域均值分解的自适应滤波降噪处理，其具体实施过程如下；
3.1，获得原始信号序列x(t)的全部极值点ni，i＝1,2,3…,I，I为极值点数，然后根据极值点，求得所有相邻局域极值点的平均值作为局域均值mi＝(ni+ni+1)/2和局域包络估计值作为局域包络值ai＝(ni-ni+1)/2；
3.2，对所有局域均值mi和所有局域包络值ai采用滑动平均法进行多次平滑处理分别得到局域均值函数m′1(t)、m′2(t)、…和局域包络函数a′1(t)、a′2(t)、…；
3.3，步骤3.2的每一次平滑处理后均处理迭代直至获得纯调频信号sk(t)：
在步骤3.2的第一次平滑处理后，从原始信号序列x(t)中分离出第一次平滑处理获得的局域均值函数m′1(t)而保留得到第一中间分离信号序列h1(t)＝x(t)-m′1(t)，然后再根据第一次平滑处理获得的局域包络函数a′1(t)进行解调得第一中间解调信号序列s1(t)＝h1(t)/a′1(t)；
然后将第一中间解调信号序列s1(t)视为原始信号序列x(t)重复步骤，结合第二次平滑处理获得的局域均值函数m′1(t)和局域包络函数a′1(t)计算获得第二中间分离信号序列和第二中间解调信号序列；
以此迭代处理直到第K次的中间解调信号序列满足-1≤sK(t)≤1且其下一次的局域包络函数a′K+1(t)＝1，K为进行迭代的次数，则将第K次的中间解调信号序列sK(t)作为纯调频信号，停止迭代处理；具体迭代过程表示为：
式中，K为进行迭代的总次数；
3.4，将以上纯调频信号sK(t)在迭代过程产生的所有局域包络函数a′k(t)相乘获得第一乘积分量的包络信号A1(t)＝a′1(t)a′2(t)…a′K(t)，k＝1,2,3…，K，接着将第一乘积分量的包络信号A1(t)和其对应的纯调频信号sK(t)相乘得到原始信号序列的第一乘积分量PF1(t)＝A1(t)sk(t)；
3.5，将第一乘积分量PF1(t)从原始信号序列x(t)中分离出来，保留得到第一残余信号u1(t)，并以第一残余信号u1(t)为原始信号序列x(t)重复上述步骤3.1-3.4迭代处理，直至通过J次循环分离出所有乘积分量，且满足第J残余信号uJ(t)为一个单调函数，从而将原始瓦斯浓度信号变换为一组乘积分量和一个残余分量的和。
3.6，对获得的各个乘积分量和原始信号序列x(t)一起进行相关性分析，取乘积分量的全局极小值作为高低频成分的调频分界点，将调频分界点之前的乘积分量的高频分量全部作为噪声乘积分量，将剔除噪声乘积分量后对信号进行重构处理，将各个乘积分量进行加合，得到降噪后的瓦斯浓度序列dx(t)，如图3所示。
步骤4、取降噪后时间步长为T的瓦斯浓度数据，构建作业类型划分敏感参数体系。本实施例选用的敏感参数主要包括有瓦斯浓度均值、均方差、方差、极差、差分特征值等统计参数，分别代表不同作业类型条件下的瓦斯涌出特征，各统计参数如表1所示。
表1瓦斯浓度序列敏感特征参数集
其中，T表示用以计算各敏感参数瓦斯数据的时间步长，Xmax表示所取时段内瓦斯浓度的最大值和Xmin分别表示所取时段内瓦斯浓度的最小值。
步骤5、构建井下工作面的作业类型神经网络智能识别模型，根据已建立的瓦斯浓度数据库和敏感特征参数输入到模型中进行训练。
本实施例选用了三个作业类型，即掘进作业、钻孔作业和无作业三种情况，分别记为w1、w2和w3，根据不同作业类型的样本分别计算各自的统计参数特征，汇总如表2所示。采用粗糙集理论和人工神经网络相结合的粗神经网络(rough membership neural network，RMNN)，提取瓦斯浓度均值、均方差、方差、极差、差分特征值等作为分类的特征量，根据各个特征量的特点进行作业类型分类。实施例中，分别根据不同的参数参数进行样本训练，并根据各参数的精度和贡献赋予权重i，而后确定综合参数Y＝d1i1+d2i2+…+d8i8。其中i1～i8分别表示前述8个敏感指标的贡献权重。
使用综合Y值和前述构建的神经网络模型，重新进行样本训练，直到能够较好的实现作业类型的精准识别，该实施例根据综合Y值的识别结果如图4所示。
表2不同作业类型样本统计参数特征值汇总
步骤6、利用训练后的模型对新获取的待测瓦斯浓度数据进行作业类型的智能识别，并自适应地进行改进。
由上述提供的算法方案可以看出，本发明所述的基于瓦斯浓度特征的作业工序智能识别方法，主要以瓦斯浓度参数为研究对象，以实现工作面作业类型智能识别为目的，进行特殊降噪处理，采用粗神经网络模型实现对井下采掘工作面工作类型的智能识别，算法的实现有助于快速有效的获取井下的作业情况，能够扩宽矿井安全监测系统在灾害诊断、作业调度和生产管理等方面的功能多样化和智能化，同时，辨识结果对于事故机理分析、应急处置和事故调查等工作也具有参考价值。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
Claims (5)
1.一种基于瓦斯浓度参数的井下工作面作业类型智能识别方法，其特征在于方法主要包括如下主要步骤：
步骤1、通过井下瓦斯浓度传感器获取瓦斯浓度的原始信号序列x(t)，存入瓦斯浓度数据库；
步骤2、将各个原始信号序列x(t)根据工作面的作业类型进行划分；
步骤3、对瓦斯浓度数据库中的瓦斯浓度数据进行自适应滤波降噪处理，获得降噪后的瓦斯浓度序列；
步骤4、将降噪后的瓦斯浓度序列作为混沌时间序列，根据作业类型的划分对各类作业类型对应的数据分别提取敏感特征参数；
步骤5、构建井下工作面的作业类型神经网络智能识别模型，根据已建立的瓦斯浓度数据库和敏感特征参数输入到模型中进行训练；
步骤6、利用训练后的模型对新获取的待测瓦斯浓度数据进行作业类型的智能识别。
2.根据权利要求1所述的一种基于瓦斯浓度参数的井下工作面作业类型智能识别方法，其特征在于：所述步骤2中，工作面的作业类型划分为放炮、割煤、落煤、交接班、打钻作业、临时支护、安装锚杆、检修和其他无关作业。
3.根据权利要求1所述的一种基于瓦斯浓度参数的井下工作面作业类型智能识别方法，其特征在于：所述步骤4中，所述的敏感特征参数包括瓦斯浓度数据的均值、均方差、方差、极差、差分特征值。
4.根据权利要求1所述的一种基于瓦斯浓度参数的井下工作面作业类型智能识别方法，其特征在于：所述步骤3中，自适应滤波降噪处理具体为：
3.1，获得原始信号序列x(t)的全部极值点ni，i＝1,2,3…,I，I为极值点数，然后根据极值点，求得所有相邻局域极值点的平均值作为局域均值mi＝(ni+ni+1)/2和局域包络估计值作为局域包络值ai＝(ni-ni+1)/2；
3.2，对所有局域均值mi和所有局域包络值ai采用滑动平均法进行多次平滑处理分别得到局域均值函数m′1(t)、m′2(t)、…和局域包络函数a′1(t)、a′2(t)、…；
3.3，步骤3.2的每一次平滑处理后均处理迭代直至获得纯调频信号sk(t)：
在步骤3.2的第一次平滑处理后，从原始信号序列x(t)中分离出第一次平滑处理获得的局域均值函数m′1(t)而保留得到第一中间分离信号序列h1(t)＝x(t)-m′1(t)，然后再根据第一次平滑处理获得的局域包络函数a′1(t)进行解调得第一中间解调信号序列s1(t)＝h1(t)/a′1(t)；
然后将第一中间解调信号序列s1(t)视为原始信号序列x(t)重复步骤，结合第二次平滑处理获得的局域均值函数m′1(t)和局域包络函数a′1(t)计算获得第二中间分离信号序列和第二中间解调信号序列；
以此迭代处理直到第K次的中间解调信号序列满足-1≤sK(t)≤1且其下一次的局域包络函数a′K+1(t)＝1，K为进行迭代的次数，则将第K次的中间解调信号序列sK(t)作为纯调频信号，停止迭代处理；具体迭代过程表示为：
式中，K为进行迭代的总次数；
3.4，将以上纯调频信号sK(t)在迭代过程产生的所有局域包络函数a′k(t)相乘获得第一乘积分量的包络信号A1(t)＝a′1(t)a′2(t)…a′K(t)，k＝1,2,3…，K，接着将第一乘积分量的包络信号A1(t)和其对应的纯调频信号sK(t)相乘得到原始信号序列的第一乘积分量PF1(t)＝A1(t)sk(t)；
3.5，将第一乘积分量PF1(t)从原始信号序列x(t)中分离出来，保留得到第一残余信号u1(t)，并以第一残余信号u1(t)为原始信号序列x(t)重复上述步骤3.1-3.4迭代处理，直至通过J次循环分离出所有乘积分量，且满足第J残余信号uJ(t)为一个单调函数，由此构建以下原始信号序列x(t)的分解表达式：
其中，j表示步骤3.1-3.4迭代处理的次序数，j＝1,2,3…，J，J表示步骤3.1-3.4迭代处理的总次数；
3.6，对获得的各个乘积分量和原始信号序列x(t)一起进行相关性分析，取所有乘积分量的能量值构成序列的全局极小值作为高低频成分的调频分界点，将调频分界点之前的乘积分量的高频分量全部作为噪声乘积分量，将剔除噪声乘积分量后对信号进行重构处理，得到降噪后的瓦斯浓度序列dx(t)，其中，乘积分量的能量值按照分量所有序列数据的平方和计算得出。
5.根据权利要求1所述的一种基于瓦斯浓度参数的井下工作面作业类型智能识别方法，其特征在于：
所述的作业类型神经网络智能识别模型采用粗糙集理论和人工神经网络相结合的粗神经网络(rough membership neural network，RMNN)，利用神经网络输入层的粗神经元分别计算各作业类型时瓦斯浓度数据敏感特征参数，并采用模糊神经元代替普通神经元。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811243460.1A CN109325470B (zh) | 2018-10-24 | 2018-10-24 | 基于瓦斯浓度参数的井下工作面作业类型智能识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811243460.1A CN109325470B (zh) | 2018-10-24 | 2018-10-24 | 基于瓦斯浓度参数的井下工作面作业类型智能识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109325470A CN109325470A (zh) | 2019-02-12 |
| CN109325470B true CN109325470B (zh) | 2021-09-03 |
Family
ID=65263225
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811243460.1A Active CN109325470B (zh) | 2018-10-24 | 2018-10-24 | 基于瓦斯浓度参数的井下工作面作业类型智能识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109325470B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112308306A (zh) * | 2020-10-27 | 2021-02-02 | 贵州工程应用技术学院 | 一种多模态输入的煤与瓦斯突出危险预测方法 |
| CN112686513B (zh) * | 2020-12-23 | 2024-07-09 | 精英数智科技股份有限公司 | 井下工作面作业状态的识别方法、装置及生产决策系统 |
| CN114046179B (zh) * | 2021-09-15 | 2023-09-22 | 山东省计算中心（国家超级计算济南中心） | 一种基于co监测数据智能识别和预测井下安全事故的方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6289328B2 (en) * | 1998-04-17 | 2001-09-11 | The United States Of America As Represented By The Secretary Of The Navy | Chemical sensor pattern recognition system and method using a self-training neural network classifier with automated outlier detection |
| CN104537199A (zh) * | 2014-10-20 | 2015-04-22 | 天地（常州）自动化股份有限公司 | 利用瓦斯浓度时间序列识别矿井工作面生产工序的方法 |
| CN106296435A (zh) * | 2016-08-18 | 2017-01-04 | 西安科技大学 | 一种矿井瓦斯监测异常数据识别方法 |
| EP3126787A2 (en) * | 2014-04-01 | 2017-02-08 | Saudi Arabian Oil Company | Tomographic imaging of multiphase flows |
| CN106503798A (zh) * | 2015-09-08 | 2017-03-15 | 中国科学院沈阳自动化研究所 | 基于粗糙集与bp神经网络的泵的故障诊断方法 |
| US9652684B2 (en) * | 2012-11-08 | 2017-05-16 | Total Sa | Image processing for classification and segmentation of rock samples |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8374974B2 (en) * | 2003-01-06 | 2013-02-12 | Halliburton Energy Services, Inc. | Neural network training data selection using memory reduced cluster analysis for field model development |
-
2018
- 2018-10-24 CN CN201811243460.1A patent/CN109325470B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6289328B2 (en) * | 1998-04-17 | 2001-09-11 | The United States Of America As Represented By The Secretary Of The Navy | Chemical sensor pattern recognition system and method using a self-training neural network classifier with automated outlier detection |
| US9652684B2 (en) * | 2012-11-08 | 2017-05-16 | Total Sa | Image processing for classification and segmentation of rock samples |
| EP3126787A2 (en) * | 2014-04-01 | 2017-02-08 | Saudi Arabian Oil Company | Tomographic imaging of multiphase flows |
| CN104537199A (zh) * | 2014-10-20 | 2015-04-22 | 天地（常州）自动化股份有限公司 | 利用瓦斯浓度时间序列识别矿井工作面生产工序的方法 |
| CN106503798A (zh) * | 2015-09-08 | 2017-03-15 | 中国科学院沈阳自动化研究所 | 基于粗糙集与bp神经网络的泵的故障诊断方法 |
| CN106296435A (zh) * | 2016-08-18 | 2017-01-04 | 西安科技大学 | 一种矿井瓦斯监测异常数据识别方法 |
Non-Patent Citations (3)
| Title |
|---|
| 基于多因素模式识别的煤与瓦斯突出预测研究;毕慧杰等;《中国安全生产科学技术》;20170630;98-103 * |
| 基于瓦斯涌出特征的机掘面作业类型判识技术;王川等;《煤矿安全》;20131231;71-74 * |
| 瓦斯浓度的分形分析与混沌预测模型研究;李刚;《中国博士学位论文全文数据库工程科技Ⅰ辑》;20110315;B021-11 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109325470A (zh) | 2019-02-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109325470B (zh) | 基于瓦斯浓度参数的井下工作面作业类型智能识别方法 | |
| AU2019214077B2 (en) | Method for dividing ecological and geological environment types based on coal resource development | |
| CN111832506B (zh) | 一种基于长时序植被指数的重建植被遥感判别方法 | |
| Nourani et al. | Applications of hybrid wavelet–artificial intelligence models in hydrology: a review | |
| CN106503792B (zh) | 一种基于自适应模块化神经网络的瓦斯浓度预测方法 | |
| CN113570000A (zh) | 一种基于多模型融合的海洋单要素观测质量控制方法 | |
| CN103646167A (zh) | 一种基于遥测数据的卫星异常状态检测系统 | |
| US11060899B2 (en) | Method for determining a maximum allowable volume of water that can be removed over time from an underground water source | |
| CN109636104A (zh) | 矿山地质环境质量演变规律的综合评价方法 | |
| Liu et al. | Developments in quantitative assessment and modeling of mineral resource potential: an overview | |
| CN110348568B (zh) | 一种适用于强电磁干扰地区的深部采空区探测方法 | |
| CN116699612A (zh) | 融合InSAR形变量与降雨量的滑坡危险性评价方法 | |
| CN104373153A (zh) | 一种煤矿井下综放工作面的煤岩性状识别方法与系统 | |
| CN110956155B (zh) | 基于co数据的综采工作面作业工序模糊聚类识别方法 | |
| CN118035890A (zh) | 基于灾情自动获取的自适应地震烈度评定模型及构建方法 | |
| CN116658246A (zh) | 一种煤矿水害监测预警系统及方法 | |
| Xie et al. | A novel quality control method of time-series ocean wave observation data combining deep-learning prediction and statistical analysis | |
| CN110348510B (zh) | 一种基于深水油气钻井过程阶段性特征的数据预处理方法 | |
| CN116011811A (zh) | 基于多源异构数据的岩溶隧道地表塌陷施工风险评价方法 | |
| CN110991363B (zh) | 不同采煤工序下煤矿安全监控系统co涌出特征提取方法 | |
| Hosseini Duki et al. | Evaluation of base flow separation methods for determining water extraction (case study: Gorganroud River Basin) | |
| CN113887049A (zh) | 一种基于机器学习的石油钻井的钻速预测方法及系统 | |
| CN113834547A (zh) | 一种河流虚拟站水位时序重建方法及系统 | |
| CN109857783B (zh) | 基于深度学习的声纳渗流运动规律发现方法 | |
| CN114818493B (zh) | 一种隧道岩体完整程度定量评价的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |