1.一种煤岩变形破裂过程的红外辐射信息去噪与监测预警方法，其特征在于：首先进行煤岩单轴压缩红外辐射监测试验，在压力机上设置一个加载煤岩试样，并在加载煤岩试样的左右设置两个对比煤岩试样，然后布置防辐射区域并设置均匀的无辐射主动光源激励环境；采集加载的煤岩试样从开始加载到压裂后全过程的红外辐射数据和两个对比煤岩试样同时间的红外辐射数据，对加载的煤岩试样进行初始去噪得到初始红外辐射温度序列，之后计算两个对比煤岩试样的对比试样噪声因子序列，过滤掉初始红外辐射温度序列中的噪声因子得到加载煤岩试样真实红外辐射温度序列，通过差值点运算计算得到加载煤岩试样在压力过程中的真实红外辐射温度差分序列，将真实红外辐射温度差分序列线性映射到灰度空间，得到真实红外辐射灰度差分图序列最后计算真实红外辐射灰度差分图序列的红外辐射信号值，之后利用红外辐射信号值曲线参考应力时间曲线，即可分析获得测试区域被测煤岩所在压力下红外辐射突变状态，从而实现红外辐射突变预警；

具体包括以下步骤：

步骤1，在压力机加载平台上待加载的煤岩试样左右两侧同一水平高度放置的与煤岩试样尺寸相同的对比煤岩试样Ⅰ和对比煤岩试样Ⅱ作对比用以消除噪声因子对试验去噪效果的影响；然后在加载煤岩试样和对比煤岩试样Ⅰ和对比煤岩试样Ⅱ的后边、左边、右边分别设置电磁屏蔽板将三个煤岩试样包围以减少辐射干扰，在没有设置电磁屏蔽板的前面使用红外热像仪录制监测；

步骤2，将压力机加载平台与红外热像仪区间内1米范围内的空间区域用防辐射黑色幕布包围，在包围区域顶部，给予无辐射主动光源激励，保证区域内光照相同，且试样表面没有由于光照不均匀产生的阴影遮掩，使煤岩试样、压力机加载平台、红外热像仪放置在同样光照环境下使光照环境保持不变，以减少误差；

步骤3，利用红外热像仪采集被煤岩试样从开始加载到压力机停止加载过程中的红外辐射变化数据，同时采集相同时间内的对比煤岩试样Ⅰ和对比煤岩试样Ⅱ的红外辐射数据，分别以加载煤岩试样和对比试样的边界为采集轮廓，定义相同大小、同一水平高度的三个监测区域，采集监测区域的红外辐射数据，得到一个加载煤岩试样和两个对比煤岩试样的红外热像温度矩阵序列，每个序列的第p帧温度矩阵表示如下：式中，IRTEGSp、IRTRGS1p、IRTRGS2p分别表示加载煤岩试样、对比煤岩试样Ⅰ、对比煤岩p p p p试样Ⅱ第p帧的红外辐射温度矩阵，xm,n xm,n、ym,n、zm,n分别为加载煤岩试样、对比煤岩试样Ⅰ、对比煤岩试样Ⅱ第p帧的红外辐射温度矩阵的第m行、n列的矩阵元素；其中，m，n分别为矩阵的最大行号和列号；

步骤4，利用对比煤岩试样Ⅰ对加载煤岩试样进行初始去噪，得到初始红外辐射温度序列；

步骤5，利用对比煤岩试样Ⅰ的红外辐射温度矩阵与对比煤岩试样Ⅱ的红外辐射温度矩阵进行差值点运算得到对比试样噪声因子序列；

步骤6，将初始红外辐射温度序列中的噪声因子过滤，得到加载煤岩试样的真实红外辐射温度序列；

步骤7，将加载煤岩试样的真实红外辐射温度序列进行间隔差值点运算，得到加载煤岩试样的真实红外辐射温度差分序列；

步骤8，将加载煤岩试样的真实红外辐射温度差分序列线性映射到灰度颜色空间，得到真实红外辐射灰度差分序列图；

步骤9，计算加载煤岩试样真实红外辐射灰度差分图序列的信号值，该信号值即为红外辐射的前兆预警指标，利用前兆预警指标即可实现煤岩破裂预警。

2.根据权利要求1所述的煤岩变形破裂过程的红外辐射信息去噪与监测预警方法，其特征在于，步骤 4 所述的利用对比煤岩试样 I 对加载煤岩试样进行初始去噪的具体计算方法为：将加载煤岩试样的第p帧IRTEGSp与对比煤岩试样Ⅰ的第p帧IRTRGS1p进行差值点运算，得到加载煤岩试样序列与对比煤岩试样Ⅰ序列的差分序列，即为初始红外辐射温度序列，记为：式中，IIRTSp为第p帧加载煤岩试样与对比煤岩试样Ⅰ红外辐射温度序列一次差分序列，pwm,n为初始红外辐射温度序列第p帧第m行、n列的矩阵元素2其中p为帧序列号，m，n分别为矩阵的行号和列号。

3.根据权利要求1所述的煤岩变形破裂过程的红外辐射信息去噪与监测预警方法，其特征在于，步骤 5 所述的运算得到对比试样噪声因子序列方法为：将对比煤岩试样Ⅰ的第p帧IRTRGS1p与对比煤岩试样Ⅱ的第p帧IRTRGS2p进行差值点运算，得到对比试样噪声矩阵，记作NFCRSP，表达式如下：式中，NFCRSP为第p帧对比煤岩试样Ⅰ与对比煤岩试样Ⅱ的红外辐射温度序列进行差值运算得到的对比试样噪声因子2，其中p为噪声因子序列号，m,n分别为矩阵的行号和列号。

4.根据权利要求3所述的煤岩变形破裂过程的红外辐射信息去噪与监测预警方法，其特征在于，步骤6所述的计算加载煤岩试样真实红外辐射温度序列方法为：将步骤4所得的初始红外辐射温度序列中的第p帧温度矩阵IIRTSp与步骤5所得的对比试样噪声因子序列中的第p个噪声因子矩阵NFCRSP进行差值点运算，得到加载煤岩试样真实红外辐射温度序列，其第p帧记作NIRTp，表达式如下所示：p

式中，NIRTp为加载煤岩试样第p帧真实红外辐射温度矩阵，tm,n表示真实红外辐射温度图序列第p帧的温度矩阵第m行、n列的元素值。

5.根据权利要求4所述的煤岩变形破裂过程的红外辐射信息去噪与监测预警方法，其特征在于，步骤 7 所述得到加载煤岩试样的真实红外辐射温度差分序列具体为：DRIRTSp(m,n)＝NIRTp(m,n)‑NIRTp‑1(m,n)

式中，DRIRISp(m,n)表示第p帧真实红外辐射温度差分矩阵，NIRTp加载煤岩试样第p帧真实红外辐射温度矩阵，NIRTp‑1加载煤岩试样第p‑1帧真实红外辐射温度矩阵.

6.根据权利要求1所述的煤岩变形破裂过程的红外辐射信息去噪与监测预警方法，其特征在于，步骤 8 所述的将加载煤岩试样的真实红外辐射温度差分序列线性映射到灰度颜色空间进行灰度空间线性映射计算方法为：p

将步骤7得到的加载煤岩试样真实红外辐射温度差分序列中的每一帧中的元素tm,n线性映射到灰度颜色空间，得到对应的加载煤岩试样二次去噪温度矩阵的灰度像素矩阵序列，称为真实红外辐射灰度差分图序列，记作NIRGISp，其表达式为：p

式中：NIRGISp(m,n)为加载煤岩试样第p帧真实红外辐射灰度差分图，tm,n表示真实红p p外辐射温度图序列第p帧的温度矩阵第m行、n列的元素值，min(tm,n)、max(tm,n)分别表示第p帧真实红外辐射温度图中像素值的最小、最大值。

7.根据权利要求6所述的煤岩变形破裂过程的红外辐射信息去噪与监测预警方法，其特征在于，步骤9所述的计算真实红外辐射灰度差分图序列的信号值的方法为：将第q帧真实红外辐射灰度差分图中像素点(i,j)处的灰度值定义为g(i,j)，该灰度值的像素点在本帧图像所有像素点所占比例定义为Pg(i,j)，则第q帧真实红外辐射灰度差分图的红外辐射信号值为：式中，IRSVq为第q帧真实红外辐射灰度差分图的红外辐射信号值。