CN110411572B - 承载煤岩破裂的红外辐射监测预警方法 - Google Patents
承载煤岩破裂的红外辐射监测预警方法 Download PDFInfo
- Publication number
- CN110411572B CN110411572B CN201910618496.1A CN201910618496A CN110411572B CN 110411572 B CN110411572 B CN 110411572B CN 201910618496 A CN201910618496 A CN 201910618496A CN 110411572 B CN110411572 B CN 110411572B
- Authority
- CN
- China
- Prior art keywords
- infrared radiation
- coal rock
- temperature
- frame
- bearing
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 142
- 239000003245 coal Substances 0.000 title claims abstract description 111
- 239000011435 rock Substances 0.000 title claims abstract description 111
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000012544 monitoring process Methods 0.000 title claims abstract description 20
- 230000035772 mutation Effects 0.000 claims abstract description 35
- 239000011159 matrix material Substances 0.000 claims description 35
- 238000012545 processing Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000006911 nucleation Effects 0.000 description 1
- 238000010899 nucleation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/72—Investigating presence of flaws
Abstract
本发明公开了承载煤岩破裂的红外辐射监测预警方法,该方法基于承载煤岩内部破裂会导致其表面红外辐射温度突变的特点,采集承载煤岩损伤破裂过程中的红外辐射信息,确定承载煤岩因微元破裂而产生红外辐射突变的温度幅值区间,统计承载煤岩红外辐射突变点的温度‑频数分布关系,计算红外辐射特征值,并对红外辐射特征值取倒数,基于红外辐射特征值倒数对承载煤岩损伤破裂进行监测预警。本发明方法实现了对承载煤岩损伤破裂尺度的评价,对提升矿井安全生产水平有重要意义。
Description
技术领域
本发明涉及承载煤岩破裂的红外辐射监测预警方法,属于矿山煤岩体破裂失稳的监测预警技术领域。
背景技术
承载煤岩损伤破裂会导致矿井冲击矿压、遗留煤柱坍塌、顶板垮落和突水等灾害,造成严重的人员伤亡与经济损失。煤岩损伤破裂过程伴随着内部微裂纹、微孔洞等缺陷的扩展、成核,最终形成宏观裂纹,导致煤岩失去承载能力。这种微缺陷由无序向有序发展的过程是煤岩表面红外辐射变化的根本原因。因此,红外辐射观测作为一种无损遥感探测技术,是对煤岩损伤破裂进行监测预警的便捷、有效手段。目前采用红外辐射监测煤岩破裂的方法是通过红外辐射特征的异常现象判断其是否破裂,无法有效反应煤岩的破裂尺度,不能准确判断煤岩的破裂状态。因此,研究一种能有效反应煤岩破裂尺度的监测预警指标是目前亟需解决的问题。
发明内容
本发明所要解决的技术问题是:提供承载煤岩破裂的红外辐射监测预警方法,通过观测煤岩表面红外辐射信息,判断煤岩破裂尺度,实现煤岩损伤破裂监测预警。
本发明为解决上述技术问题采用以下技术方案:
承载煤岩破裂的红外辐射监测预警方法,包括如下步骤:
步骤1,采集承载煤岩损伤破裂过程中的红外辐射信息以及参照煤岩的红外辐射信息,并对红外辐射信息进行中值和高斯滤波处理,去除噪声干扰,得到参照煤岩和承载煤岩各自对应的原始红外辐射温度矩阵;
步骤2,对参照煤岩和承载煤岩各自对应的原始红外辐射温度矩阵分别进行前向一阶差分处理,并取绝对值,得到参照煤岩和承载煤岩各自对应的红外辐射变温矩阵;
步骤3,确定承载煤岩产生红外辐射突变的温度幅值区间,计算参照煤岩红外辐射变温矩阵中所有元素的平均值,将该平均值作为温度幅值区间的最小值,将承载煤岩红外辐射变温矩阵中所有元素的最大值作为温度幅值区间的最大值;
步骤4,在第p帧对应的温度幅值区间[Min(p),Max(p)]中,以Min(p)为起始温度,以0.01℃为温度间隔,统计红外辐射突变点Tp的频数Np(Tp),得到承载煤岩红外辐射变温矩阵中红外辐射突变点的温度-频数分布,Tp∈[Min(p),Max(p)],Tp=Min(p)+0.01×s,s∈[0,a],且s为整数,Min(p)为第p帧温度幅值区间的最小值,Max(p)为第p帧温度幅值区间的最大值;
步骤5,对红外辐射突变点的频数Np(Tp)进行对数运算,得到红外辐射突变点的温度-对数频数,根据红外辐射突变点的温度-对数频数求得第p帧对应的红外辐射特征值;
步骤6,对第p帧对应的红外辐射特征值取倒数,得到所有帧红外辐射特征值倒数,求取所有帧红外辐射特征值倒数的平均值和最小值,并用平均值的2倍减去最小值,将所得值作为判定承载煤岩红外辐射特征值倒数发生突变的阈值,将红外辐射特征值倒数超过该阈值的时刻,作为承载煤岩发生破裂的时刻。
作为本发明的一种优选方案,步骤2所述红外辐射变温矩阵的计算公式为:
作为本发明的一种优选方案,步骤3所述确定承载煤岩产生红外辐射突变的温度幅值区间的计算公式为:
Max(p)=Max(φjz(i,j,p))
其中,φcz(i,j,p)表示参照煤岩第p帧的红外辐射变温矩阵的第i行第j列元素,φjz(i,j,p)表示承载煤岩第p帧的红外辐射变温矩阵的第i行第j列元素,m、n分别为矩阵的行、列数,Min(p)为第p帧温度幅值区间的最小值,Max(p)为第p帧温度幅值区间的最大值。
作为本发明的一种优选方案,步骤5所述红外辐射特征值的计算方法为:
对红外辐射突变点的频数Np(Tp)进行对数运算,得到红外辐射突变点的对数频数Mp(Tp):
Mp(Tp)=log10(Np(Tp))
其中,为第p帧对应的红外辐射特征值, Min(p)为第p帧温度幅值区间的最小值,Max(p)为第p帧温度幅值区间的最大值,h、k均为参数,Tp(s)为第p帧承载煤岩的红外辐射温度幅值区间的某温度,Mp(s)为第p帧承载煤岩的红外辐射温度幅值区间的某温度对应的对数频数。
作为本发明的一种优选方案,步骤6所述对第p帧对应的红外辐射特征值取倒数,表示为:
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
本发明确定了承载煤岩因微元破裂而产生的红外辐射突变的温度幅值区间,统计了承载煤岩红外辐射突变点的温度-频数分布关系,计算了红外辐射特征值,并以红外辐射特征值倒数为指标,对承载煤岩损伤破裂进行监测预警,实现了煤岩破裂尺度的评价,提高了矿山的安全生产水平。
附图说明
图1是本发明承载煤岩破裂的红外辐射监测预警方法的流程图。
图2是本发明中某两帧红外辐射变温矩阵的红外辐射突变点温度-频数(Np(Tp)-Tp)分布,其中,(a)为第7帧,(b)为第7447帧。
图3是本发明中某两帧红外辐射变温矩阵的红外辐射突变点温度-对数频数(Mp(Tp)-Tp)分布,其中,(a)为第7帧,(b)为第7447帧。
图4是本发明中某两帧红外辐射变温矩阵的红外辐射突变点温度-对数频数的回归方程曲线,其中,(a)为第7帧,(b)为第7447帧。
图5是本发明中某岩样的红外辐射特征值倒数(RIRE)和应力随时间变化的曲线。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
如图1所示为本发明承载煤岩破裂的红外监测预警方法的流程图,包括以下步骤:
a、采集承载煤岩损伤破裂过程中的红外辐射信息,并对红外辐射数据进行中值和高斯滤波处理,去除噪声干扰。
b、建立参照和承载煤岩的红外辐射变温矩阵:对去噪后的参照煤岩和加载煤岩原始红外辐射温度矩阵分别进行前向一阶差分处理,并取绝对值,得到红外辐射变温矩阵序列:
c、确定承载煤岩因微元破裂而产生的红外辐射突变的温度幅值区间:计算每一帧参照煤岩红外辐射变温矩阵的平均值,为承载煤岩因微元破裂而产生的红外辐射突变阈值,以红外辐射突变阈值为温度幅值区间的最小值;计算加载煤岩每一帧红外辐射变温矩阵的最大值,为温度幅值区间的最大值:
Max(p)=Max(φjz(i,j,p))
Tp∈[Min(p),Max(p)]
其中:φcz(i,j,p)为参照煤岩第p帧的红外辐射变温矩阵,φjz(i,j,p)为加载煤岩第p帧的红外辐射变温矩阵,Min(p)为第p帧温度幅值最小值,Max(p)为第p帧温度幅值最大值,Tp为第p帧承载煤岩因微元破裂而产生的红外辐射突变温度。
d、统计每帧加载煤岩红外辐射变温矩阵中红外辐射突变点的温度-频数分布:以Min(p)为起始温度,以0.01℃为温度间隔,统计温度为Tp的红外辐射突变点频数Np(Tp),得到加载煤岩红外辐射变温矩阵中红外辐射突变点的温度-频数分布,具体示例如图2的(a)、(b)所示。
e、计算红外辐射特征值:对红外辐射突变频数Np进行对数运算,得到红外辐射突变点的对数频数,具体示例如图3的(a)、(b)所示:
Mp=log10(Np)
其中:Mp为红外辐射突变点的对数频数。
f、对承载煤岩损伤破裂进行监测预警:对步骤e得到的红外辐射特征值取倒数:
其中:RIREp为第p帧红外辐射特征值倒数。
求取所有帧RIRE的平均值和最小值,用平均值的2倍减去最小值,并将所得值作为判定承载煤岩RIRE发生突变的阈值,当RIRE超出阈值时视为发生突变。找出红外辐射特征值倒数(RIRE)的突变时刻,作为煤岩发生大尺度破裂时刻,且RIRE的突变幅度越大,承载煤岩的损伤破裂尺度越大。
本发明确定了承载煤岩因微元破裂而产生的红外辐射突变的温度幅值区间,统计了承载煤岩红外辐射突变点的温度-频数分布关系,计算了红外辐射特征值,并以红外辐射特征值倒数为指标,对承载煤岩损伤破裂进行监测预警,实现了煤岩破裂尺度的评价,提高了矿山的安全生产水平。
如图4的(a)、(b)所示,红外辐射突变点的温度-对数频数分布符合一次线性关系曲线。图5所示承载煤岩在发生较大破裂前,RIRE会发生突变,且破裂尺度越大,RIRE的突变幅度越大。因此,利用RIRE能够很好地对承载煤岩损伤破裂进行监测预警。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (5)
1.承载煤岩破裂的红外辐射监测预警方法,其特征在于,包括如下步骤:
步骤1,采集承载煤岩损伤破裂过程中的红外辐射信息以及参照煤岩的红外辐射信息,并对红外辐射信息进行中值和高斯滤波处理,去除噪声干扰,得到参照煤岩和承载煤岩各自对应的原始红外辐射温度矩阵;
步骤2,对参照煤岩和承载煤岩各自对应的原始红外辐射温度矩阵分别进行前向一阶差分处理,并取绝对值,得到参照煤岩和承载煤岩各自对应的红外辐射变温矩阵;
步骤3,确定承载煤岩产生红外辐射突变的温度幅值区间,计算参照煤岩红外辐射变温矩阵中所有元素的平均值,将该平均值作为温度幅值区间的最小值,将承载煤岩红外辐射变温矩阵中所有元素的最大值作为温度幅值区间的最大值;
步骤4,在第p帧对应的温度幅值区间[Min(p),Max(p)]中,以Min(p)为起始温度,以0.01℃为温度间隔,统计第p帧承载煤岩的红外辐射温度幅值区间内温度Tp的频数Np(Tp),得到承载煤岩红外辐射变温矩阵中红外辐射突变点的温度-频数分布,Tp∈[Min(p),Max(p)],Tp=Min(p)+0.01×s,s∈[0,a],且s为整数,Min(p)为第p帧温度幅值区间的最小值,Max(p)为第p帧温度幅值区间的最大值;
步骤5,对第p帧承载煤岩的红外辐射温度幅值区间内温度Tp的频数Np(Tp)进行对数运算,得到红外辐射突变点的温度-对数频数,根据红外辐射突变点的温度-对数频数求得第p帧对应的红外辐射特征值;
步骤6,对第p帧对应的红外辐射特征值取倒数,得到所有帧红外辐射特征值倒数,求取所有帧红外辐射特征值倒数的平均值和最小值,并用平均值的2倍减去最小值,将所得值作为判定承载煤岩红外辐射特征值倒数发生突变的阈值,将红外辐射特征值倒数超过该阈值的时刻,作为承载煤岩发生破裂的时刻。
4.根据权利要求1所述承载煤岩破裂的红外辐射监测预警方法,其特征在于,步骤5所述红外辐射特征值的计算方法为:
对第p帧承载煤岩的红外辐射温度幅值区间内温度Tp的频数Np(Tp)进行对数运算,得到第p帧承载煤岩的红外辐射温度幅值区间内温度Tp对应的对数频数Mp(Tp):
Mp(Tp)=log10(Np(Tp))
Min(p)为第p帧温度幅值区间的最小值,Max(p)为第p帧温度幅值区间的最大值,h、k均为参数,Tp(s)为第p帧承载煤岩的红外辐射温度幅值区间的某温度,Mp(s)为第p帧承载煤岩的红外辐射温度幅值区间内某温度对应的对数频数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910618496.1A CN110411572B (zh) | 2019-07-10 | 2019-07-10 | 承载煤岩破裂的红外辐射监测预警方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910618496.1A CN110411572B (zh) | 2019-07-10 | 2019-07-10 | 承载煤岩破裂的红外辐射监测预警方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110411572A CN110411572A (zh) | 2019-11-05 |
CN110411572B true CN110411572B (zh) | 2020-12-15 |
Family
ID=68360849
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910618496.1A Active CN110411572B (zh) | 2019-07-10 | 2019-07-10 | 承载煤岩破裂的红外辐射监测预警方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110411572B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111811924B (zh) * | 2020-07-06 | 2023-12-12 | 安徽理工大学 | 一种判别岩石扩容起始点的红外试验方法 |
CN111811933B (zh) * | 2020-07-31 | 2022-03-11 | 中国矿业大学 | 一种承载煤岩损伤破裂过程中的红外辐射信息去噪方法 |
CN113188909B (zh) * | 2021-04-28 | 2022-03-01 | 中国矿业大学 | 承载煤岩裂纹萌生和稳定扩展起点的红外辐射识别方法 |
CN113484145B (zh) * | 2021-07-06 | 2022-07-05 | 中国矿业大学 | 煤岩变形破裂过程的红外辐射信息去噪与监测预警方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104764528A (zh) * | 2015-04-03 | 2015-07-08 | 中国矿业大学 | 一种煤岩裂隙发育过程中的热红外信息去噪方法 |
CN106018096A (zh) * | 2016-07-20 | 2016-10-12 | 中国矿业大学 | 煤岩破裂过程中裂隙发育区的红外辐射监测定位方法 |
CN109443543A (zh) * | 2018-11-23 | 2019-03-08 | 中国矿业大学 | 承载煤岩损伤演化的红外辐射量化表征方法 |
-
2019
- 2019-07-10 CN CN201910618496.1A patent/CN110411572B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104764528A (zh) * | 2015-04-03 | 2015-07-08 | 中国矿业大学 | 一种煤岩裂隙发育过程中的热红外信息去噪方法 |
CN106018096A (zh) * | 2016-07-20 | 2016-10-12 | 中国矿业大学 | 煤岩破裂过程中裂隙发育区的红外辐射监测定位方法 |
CN109443543A (zh) * | 2018-11-23 | 2019-03-08 | 中国矿业大学 | 承载煤岩损伤演化的红外辐射量化表征方法 |
Non-Patent Citations (3)
Title |
---|
"Characteristics of Infrared Radiation of Coal Specimens Under Uniaxial Loading";Liqiang Ma et.al;《Rock Mech Rock Eng》;20161231;1567-1572页 * |
"复合煤岩受载破裂内部红外辐射温度变化规律";杨帧 等;《煤炭学报》;20160331;618-624页 * |
"煤岩受压过程中内部红外辐射温度变化特征研究";马立强 等;《中国矿业大学学报》;20130531;331-336页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110411572A (zh) | 2019-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110411572B (zh) | 承载煤岩破裂的红外辐射监测预警方法 | |
CN110648480B (zh) | 一种基于变化速率的单变量报警系统及方法 | |
CN112348237B (zh) | 一种动态钻井数据异常趋势检测方法 | |
CN110703214B (zh) | 一种气象雷达状态评估和故障监测方法 | |
CN116881745B (zh) | 基于大数据的压力变送器异常监测方法 | |
CN110765983B (zh) | 一种基于微震频谱特征的煤柱型冲击地压危险性评价方法 | |
CA3142064A1 (en) | Infrared radiation quantitative characterization method for damage evolution of coal-bearing rocks | |
CN111859732A (zh) | 船闸闸门及其支承运转件损伤程度自动监测系统及监测方法 | |
CN114135477A (zh) | 一种机泵设备状态监测动态阈值预警方法 | |
CN117072460B (zh) | 一种基于振动数据和专家经验的离心泵状态监测方法 | |
US8416086B2 (en) | Methods for improved forewarning of condition changes in monitoring physical processes | |
CN108508860B (zh) | 一种基于耦合关系的流程工业生产系统数据监测方法 | |
CN110609324B (zh) | 一种用于深埋隧道岩爆预警微震事件筛选的方法 | |
CN115950666B (zh) | 一种淀粉葡萄糖生产设备故障智能检测方法 | |
CN114113217A (zh) | 一种煤岩体损伤程度的红外辐射量化评价方法 | |
CN115041544A (zh) | 冲压件的异常检测方法、装置及存储介质 | |
CN115578841A (zh) | 基于Z-score模型的盾构异常数据检测方法及系统 | |
CN115163195A (zh) | 煤矿采煤工作面瓦斯传感器数值失真识别方法与系统 | |
CN115083123A (zh) | 一种以实测数据为驱动的矿井煤自燃智能分级预警方法 | |
CN109555561B (zh) | 一种矿压预测预报方法及其系统 | |
CN114233397A (zh) | 一种基于空陆两栖机器人的隧道施工岩爆预警系统 | |
CN113685166A (zh) | 一种钻井事故预警方法及系统 | |
JP6033563B2 (ja) | 地殻破壊予測方法 | |
CN106291344B (zh) | 基于分闸压力曲线的断路器刚分点检测方法及系统 | |
CN116992389B (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 |