CN109443543B - 承载煤岩损伤演化的红外辐射量化表征方法 - Google Patents
承载煤岩损伤演化的红外辐射量化表征方法 Download PDFInfo
- Publication number
- CN109443543B CN109443543B CN201811412926.6A CN201811412926A CN109443543B CN 109443543 B CN109443543 B CN 109443543B CN 201811412926 A CN201811412926 A CN 201811412926A CN 109443543 B CN109443543 B CN 109443543B
- Authority
- CN
- China
- Prior art keywords
- infrared radiation
- coal rock
- bearing
- damage
- coal
- 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 124
- 239000011435 rock Substances 0.000 title claims abstract description 108
- 239000003245 coal Substances 0.000 title claims abstract description 105
- 238000012512 characterization method Methods 0.000 title claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 238000000034 method Methods 0.000 claims abstract description 15
- 230000001186 cumulative effect Effects 0.000 claims description 19
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 238000005065 mining Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 2
- 206010016807 Fluid retention Diseases 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 230000036962 time dependent Effects 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
- G01J5/0003—Radiation pyrometry, e.g. infrared or optical thermometry for sensing the radiant heat transfer of samples, e.g. emittance meter
-
- 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
- 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
- G01J5/48—Thermography; Techniques using wholly visual means
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种承载煤岩损伤演化的红外辐射量化表征方法，此方法以参照煤岩红外辐射变温矩阵中的最大值，作为承载煤岩因微元破坏导致红外辐射变化的门槛值，筛选出承载煤岩损伤演化过程中的红外辐射累计数，基于连续损伤力学理论，构建承载煤岩损伤演化的红外辐射量化表征方法，通过实时监测红外辐射数据，可以动态精准评价煤岩体损伤演化的实时过程、提升矿山安全生产水平以及提高现代岩石力学与工程领域煤岩体监测的稳定性和精准度。
Description
技术领域
本发明涉及一种承载煤岩损伤演化的红外辐射量化表征方法，属于保水采煤与岩层控制领域。
背景技术
煤岩承载受力、变形及破裂的过程，是一个从内部弥散的微缺陷无序分布向有序发展，最终形成宏观裂纹的过程，这也是导致煤岩表面红外辐射发生变化的根源。煤岩每一个体元的破裂都会对红外辐射的变化有一份贡献，煤岩损伤演化程度与其表面红外辐射响应信息之间存在着规律性的量化关系。因此，通过观测煤岩表面红外辐射，能够及时掌握煤岩的变形情况，预测煤岩体的损伤演化程度。但是，现有技术中只是对煤岩表面红外辐射进行定性分析，尚未找到合适红外辐射表征参数去实现损伤参数的精细标定，更没有形成量化表征方法。
发明内容
本发明目的在于提供一种承载煤岩损伤演化的红外辐射量化表征方法，能够通过实时监测采掘面红外辐射数据，动态精准评价煤岩体损伤演化。
为实现上述目的，本发明一种承载煤岩损伤演化的红外辐射量化表征方法，包括以下步骤：
第一步、采集煤岩的原始红外辐射矩阵序列；
第二步、建立红外辐射变温矩阵序列：
对第一步中的原始红外辐射矩阵序列进行一阶前向差分处理，得到煤岩的差分红外辐射矩阵序列，对差分红外辐射矩阵序列中每个元素取绝对值，得到红外辐射变温矩阵序列；
F(i,j,p)＝|f(i,j,p+1)-f(i,j,p)|
其中：F(i,j,p)为煤岩第p帧的红外辐射变温矩阵，
f(i,j,p)为煤岩第p帧的原始红外辐射矩阵，
i为红外辐射矩阵的行号；j为红外辐射矩阵的列号；
第三步、取参照煤岩的红外辐射变温矩阵序列中每一帧的最大值，作为承载煤岩红外辐射变温矩阵序列相应时刻因微元破坏导致红外辐射变化的门槛值：
M(p)＝Max(F(i,j,p))
其中：M(p)为承载煤岩第p帧红外辐射变温矩阵因微元破坏导致红外辐射变化的门槛值；
第四步、统计承载煤岩因微元破坏的红外辐射计数：
对承载煤岩红外辐射变温矩阵序列中大于相应时刻门槛值的红外辐射点进行计数，记为Q，承载煤岩第p帧的红外辐射计数为Q(p)；
第五步、计算承载煤岩因微元破坏的红外辐射累积计数：
对第四步中获得的承载煤岩红外辐射计数Q进行求和，得到承载煤岩的红外辐射累积计数，记为N，承载煤岩第p帧的红外辐射累积计数为N(p)：
其中：k为红外辐射矩阵序列的帧号；
承载煤岩全破坏时刻红外辐射累积计数为N(m)：
其中N为承载煤岩全破坏时刻红外辐射累积计数；m为承载煤岩全破坏时刻红外辐射变温矩阵序列的帧数；k为红外辐射矩阵序列的帧号；
第五步、承载煤岩损伤参量的红外辐射量化表征：
D(p)＝N(p)/N(m)
其中：D(p)为承载煤岩第p帧的损伤参量。
本发明确定了承载煤岩因微元破坏导致红外辐射变化的门槛值，筛选出承载煤岩损伤演化过程中的红外辐射累计数，构建承载煤岩损伤演化的红外辐射量化表征方法，实现了利用煤岩损伤过程中表面红外辐射响应信息来表征损伤演化，通过实时监测红外辐射数据，可以动态精准评价煤岩体损伤演化的实时过程、提升矿山安全生产水平以及提高现代岩石力学与工程领域煤岩体监测的稳定性和精准度。
附图说明
图1为本发明中某一承载煤岩损伤参量和应力随时间变化曲线图。
具体实施方式
下面结合附图对本发明进一步说明：
一种承载煤岩损伤演化的红外辐射量化表征方法，包括以下步骤：
第一步、采集煤岩的原始红外辐射矩阵序列：
在采掘面布置参照煤岩，利用红外热像仪探测并存储采掘面煤岩体和参照煤岩的表面红外辐射信息，得到承载煤岩和参照煤岩的原始红外辐射矩阵序列。
第二步、建立红外辐射变温矩阵序列：
对第一步中的原始红外辐射矩阵序列进行一阶前向差分处理，得到煤岩的差分红外辐射矩阵序列，对差分红外辐射矩阵序列中每个元素取绝对值，得到红外辐射变温矩阵序列；
F(i,j,p)＝|f(i,j,p+1)-f(i,j,p)|
其中：F(i,j,p)为煤岩第p帧的红外辐射变温矩阵，
f(i,j,p)为煤岩第p帧的原始红外辐射矩阵，
i为红外辐射矩阵的行号；j为红外辐射矩阵的列号；
第三步、取参照煤岩的红外辐射变温矩阵序列中每一帧的最大值，作为承载煤岩红外辐射变温矩阵序列相应时刻因微元破坏导致红外辐射变化的门槛值：
M(p)＝Max(F(i,j,p))
其中：M(p)为承载煤岩第p帧红外辐射变温矩阵因微元破坏导致红外辐射变化的门槛值；
第四步、统计承载煤岩因微元破坏的红外辐射计数：
对承载煤岩红外辐射变温矩阵序列中大于相应时刻门槛值的红外辐射点进行计数，记为Q，承载煤岩第p帧的红外辐射计数为Q(p)；
第五步、计算承载煤岩因微元破坏的红外辐射累积计数：
对第四步中获得的承载煤岩红外辐射计数Q进行求和，得到承载煤岩的红外辐射累积计数，记为N，承载煤岩第p帧的红外辐射累积计数为N(p)：
其中：k为红外辐射矩阵序列的帧号；
承载煤岩全破坏时刻红外辐射累积计数为N(m)：
其中N为承载煤岩全破坏时刻红外辐射累积计数；m为承载煤岩全破坏时刻红外辐射变温矩阵序列的帧数；k为红外辐射矩阵序列的帧号；
第五步、承载煤岩损伤参量的红外辐射量化表征：
D(p)＝N(p)/N(m)
其中：D(p)为承载煤岩第p帧的损伤参量。
本发明确定了承载煤岩因微元破坏导致红外辐射变化的门槛值，筛选出承载煤岩损伤演化过程中的红外辐射累计数，构建承载煤岩损伤演化的红外辐射量化表征方法，实现了利用煤岩损伤过程中表面红外辐射响应信息来表征损伤演化，通过实时监测红外辐射数据，可以动态精准评价煤岩体损伤演化的实时过程、提升矿山安全生产水平以及提高现代岩石力学与工程领域煤岩体监测的稳定性和精准度。
从图1可知，以红外辐射累积计数为表征参数的煤岩损伤演化曲线，有明显的阶段性变化特征，能够很好地反映承载煤岩裂隙缺陷产生、发展和破坏的演化过程。
Claims (1)
1.一种承载煤岩损伤演化的红外辐射量化表征方法，包括以下步骤：
第一步、采集煤岩的原始红外辐射信息，得到承载煤岩和参照煤岩的原始红外辐射矩阵序列；
第二步、建立红外辐射变温矩阵序列：
对第一步中的原始红外辐射矩阵序列进行一阶前向差分处理，得到煤岩的差分红外辐射矩阵序列，对差分红外辐射矩阵序列中每个元素取绝对值，得到红外辐射变温矩阵序列；
F(i，j，p)＝|f(i，j，p+1)-f(i，j，p)|
其中：F(i，j，p)为煤岩第p帧的红外辐射变温矩阵，
f(i，j，p)为煤岩第p帧的原始红外辐射矩阵，
i为红外辐射矩阵的行号；j为红外辐射矩阵的列号；
第三步、取参照煤岩的红外辐射变温矩阵序列中每一帧的最大值，作为承载煤岩红外辐射变温矩阵序列相应时刻因微元破坏导致红外辐射变化的门槛值：
M(p)＝Max(F(i,j,p))
其中：M(p)为承载煤岩第p帧红外辐射变温矩阵因微元破坏导致红外辐射变化的门槛值；
第四步、统计承载煤岩因微元破坏的红外辐射计数：
对承载煤岩红外辐射变温矩阵序列中大于相应时刻门槛值的红外辐射点进行计数，记为Q，承载煤岩第p帧的红外辐射计数为Q(p)；
第五步、计算承载煤岩因微元破坏的红外辐射累积计数：
对第四步中获得的承载煤岩红外辐射计数Q进行求和，得到承载煤岩的红外辐射累积计数，记为N，承载煤岩第p帧的红外辐射累积计数为N(p)：
其中：k为红外辐射矩阵序列的帧号；
承载煤岩全破坏时刻红外辐射累积计数为N(m)：
其中N(m)为承载煤岩全破坏时刻红外辐射累积计数；m为承载煤岩全破坏时刻红外辐射变温矩阵序列的帧数；k为红外辐射矩阵序列的帧号；
第六步、承载煤岩损伤参量的红外辐射量化表征：
D(p)＝N(p)/N(m)
其中：D(p)为承载煤岩第p帧的损伤参量。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811412926.6A CN109443543B (zh) | 2018-11-23 | 2018-11-23 | 承载煤岩损伤演化的红外辐射量化表征方法 |
| PCT/CN2019/102488 WO2020103508A1 (zh) | 2018-11-23 | 2019-08-26 | 承载煤岩损伤演化的红外辐射量化表征方法 |
| CA3142064A CA3142064A1 (en) | 2018-11-23 | 2019-08-26 | Infrared radiation quantitative characterization method for damage evolution of coal-bearing rocks |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811412926.6A CN109443543B (zh) | 2018-11-23 | 2018-11-23 | 承载煤岩损伤演化的红外辐射量化表征方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109443543A CN109443543A (zh) | 2019-03-08 |
| CN109443543B true CN109443543B (zh) | 2021-05-04 |
Family
ID=65554250
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811412926.6A Active CN109443543B (zh) | 2018-11-23 | 2018-11-23 | 承载煤岩损伤演化的红外辐射量化表征方法 |
Country Status (3)
| Country | Link |
|---|---|
| CN (1) | CN109443543B (zh) |
| CA (1) | CA3142064A1 (zh) |
| WO (1) | WO2020103508A1 (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109443543B (zh) * | 2018-11-23 | 2021-05-04 | 中国矿业大学 | 承载煤岩损伤演化的红外辐射量化表征方法 |
| CN110411572B (zh) * | 2019-07-10 | 2020-12-15 | 中国矿业大学 | 承载煤岩破裂的红外辐射监测预警方法 |
| CN112665731A (zh) * | 2020-12-31 | 2021-04-16 | 中国矿业大学 | 红外辐射技术监测煤岩破坏失稳的分级预警方法 |
| CN112816072B (zh) * | 2021-01-12 | 2024-05-03 | 江苏师范大学 | 水岩作用下煤岩压缩热辐射温度时空分布及预测的方法 |
| CN113029995B (zh) * | 2021-03-10 | 2022-09-27 | 太原理工大学 | 一种线性调频的煤岩辐射检测装置及方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0707206A1 (en) * | 1994-10-11 | 1996-04-17 | United Technologies Corporation | High temperature crack monitoring apparatus |
| CN102095755A (zh) * | 2010-12-09 | 2011-06-15 | 重庆建工市政交通工程有限责任公司 | 一种混凝土结构的无损检测方法 |
| CN102853916A (zh) * | 2012-09-26 | 2013-01-02 | 中国神华能源股份有限公司 | 一种用于对煤垛表面进行远距离红外测温的方法及系统 |
| CN106018096A (zh) * | 2016-07-20 | 2016-10-12 | 中国矿业大学 | 煤岩破裂过程中裂隙发育区的红外辐射监测定位方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103983514B (zh) * | 2014-05-22 | 2016-06-01 | 中国矿业大学 | 一种煤岩裂隙发育红外辐射监测试验方法 |
| CN104764528B (zh) * | 2015-04-03 | 2018-01-12 | 中国矿业大学 | 一种煤岩裂隙发育过程中的热红外信息去噪方法 |
| CN109443543B (zh) * | 2018-11-23 | 2021-05-04 | 中国矿业大学 | 承载煤岩损伤演化的红外辐射量化表征方法 |
-
2018
- 2018-11-23 CN CN201811412926.6A patent/CN109443543B/zh active Active
-
2019
- 2019-08-26 WO PCT/CN2019/102488 patent/WO2020103508A1/zh active Application Filing
- 2019-08-26 CA CA3142064A patent/CA3142064A1/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0707206A1 (en) * | 1994-10-11 | 1996-04-17 | United Technologies Corporation | High temperature crack monitoring apparatus |
| CN102095755A (zh) * | 2010-12-09 | 2011-06-15 | 重庆建工市政交通工程有限责任公司 | 一种混凝土结构的无损检测方法 |
| CN102853916A (zh) * | 2012-09-26 | 2013-01-02 | 中国神华能源股份有限公司 | 一种用于对煤垛表面进行远距离红外测温的方法及系统 |
| CN106018096A (zh) * | 2016-07-20 | 2016-10-12 | 中国矿业大学 | 煤岩破裂过程中裂隙发育区的红外辐射监测定位方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CA3142064A1 (en) | 2020-05-28 |
| CN109443543A (zh) | 2019-03-08 |
| WO2020103508A1 (zh) | 2020-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109443543B (zh) | 承载煤岩损伤演化的红外辐射量化表征方法 | |
| Huang et al. | Machine learning-based automatic control of tunneling posture of shield machine | |
| JP2020500420A (ja) | マシンラーニング基盤の半導体製造の収率予測システム及び方法 | |
| CN110442666B (zh) | 一种基于神经网络模型的矿产资源预测方法及系统 | |
| CN110298503B (zh) | 基于微震信息和深度卷积神经网络的隧洞岩爆预警方法 | |
| CN105260599A (zh) | 一种基于bp神经网络建模的冲击地压动态预测方法 | |
| CN107194478A (zh) | 融合寿命数据和性能退化数据的单机剩余寿命预测方法 | |
| CN107368463B (zh) | 基于光纤光栅传感网络数据的巷道非线性变形预测方法 | |
| CN112699572A (zh) | 一种基于时滞相关性分析预测滑坡变形的方法 | |
| CN105239608A (zh) | 一种基于wt-rs-svr相结合的滑坡位移预测方法 | |
| CN112948932A (zh) | 一种基于TSP预报数据与XGBoost算法的围岩等级预测方法 | |
| CN109252855B (zh) | 确定气井最终累积产量的方法及装置 | |
| CN114444258A (zh) | 一种基于切线角与多种滑坡模型的智能滑坡预测方法 | |
| CN106997334A (zh) | 一种基于时间加权的矿山压力数据处理系统及方法 | |
| CN114548482A (zh) | 一种基于蠕变型滑坡动能变化率临滑预警方法 | |
| CN105205002A (zh) | 一种基于测试工作量的软件安全缺陷发现模型的建模方法 | |
| CN110222825B (zh) | 一种水泥成品比表面积预测方法及系统 | |
| CN115759409A (zh) | 多时间模式注意力机制优化lstm模型的水闸形变预测方法 | |
| CN114662793A (zh) | 基于可解释分层模型的业务流程剩余时间预测方法与系统 | |
| CN107290800A (zh) | 实用的钻前测井曲线预测方法 | |
| Morgenroth et al. | A novel long-short term memory network approach for stress model updating for excavations in high stress environments | |
| CN117113644A (zh) | 一种基于深部监测变形序列的边坡临滑预报方法、系统及介质 | |
| CN109272007B (zh) | 基于深度神经网络的初撑力、末阻力识别方法、存储介质 | |
| CN110780342A (zh) | 一种岩石边坡变形预警方法 | |
| CN105069214A (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 |