基于图像处理技术的岩体裂隙及飞石运动分析系统及方法

本发明属于岩土工程领域，具体涉及一种基于图像处理技术的岩体裂隙及飞石运动分析系统及方法。

背景技术：

近年来，随着建筑工程、水利水电工程、地下工程等国家基础设施的不断发展，产生了大量的边坡、隧道、矿山工程项目；众多深部开采的金属矿山、大型水利工程的建设都涉及到岩土领域的安全问题，开展建设项目时，常出现岩体垮塌、滑坡、飞石等地质事故，对安全生产及社会效益造成了巨大的不利影响。自然界中存在的岩石，具有较多的天然损伤，是一种典型的非连续介质，岩体内结构相当复杂，体现在晶面、空隙、片理、层理、节理、断层等非连续性形式。在工程建设过程中，岩体接受多种异向载荷作用，内部的微裂隙开始发育，起裂、扩展、贯通，最后发生破坏，例如山体滑坡、地层涌水、突泥、顶板垮塌、瓦斯突出等，都与岩体裂隙发育有着直接关系。当进行深部开采时，岩体具有非常高的构造应力，甚至会发生突发式岩体破坏，即岩爆。当岩层发生岩爆现象时，伴随产生的巨大声波及气浪，将破碎的岩块抛出，可能造成飞石伤人事故，同时砸坏机器设备。因此，监测岩体工程项目中岩石裂纹萌发、扩展及破坏具有重大的应用意义。

为监测在工程开发、资源开采过程中岩体裂隙发育岩体破坏后产生的飞石运动规律，减少岩体地质破坏事故，国内外开展了诸多研究：常见的岩体微细观分析方法有断口形貌扫描仪、热红外辐射、声发射、电磁辐射法、ct扫描等。断口形貌扫描仪利用各种显微镜如光学显微镜、扫描电子显微镜，将其应用于岩石断面、裂纹扩展分析，清华大学、中国矿业大学与日本岛津制作社研发的sem加载试验系统，可实现在各种温度、载荷下对试样表面的损伤演化情况；热红外辐射利用岩石破裂时产生的红外辐射，研究其能量耗散过程中辐射的各种波段的遥感前兆；电磁辐射法主要利用压电效应，对岩体破裂失稳过程中电磁辐射产生机理、特征、规律以及传播特性进行分析，常用于突出、冲击地亚事故的预报中；声发射是常用的岩石损伤分析技术，其利用岩石受载时，岩体内部微裂纹产生、扩展以及裂纹面间的摩擦产生的次声波、声波、超声波，进而分析声波的振幅和频率等参数，反映出岩石裂纹的数量、大小、位置等；ct扫描技术采用x射线对岩石进行断面扫描，根据扫描分辨单元，结合三维重构法，可以判断与重现出岩体内部结构图像，从而无损地检测出岩石内部裂隙变化情况。上述技术可在一定程度可反映和解决岩石裂纹发育规律问题，但仍有待进一步的研究发展，例如断口形貌扫描仪、电测辐射法、ct扫描等技术多应用于实验室研究，对于工程建设现场岩体检测缺乏指导性，同时对岩石破碎后产生飞石运动规律的研究较少；声发射及热红外辐射对仪器要求高，在较为复杂地下空间或环境条件使用时精度往往不高，在实际监测时也非常依赖操作者对仪器的操控能力及熟练程度。

现有技术中的岩体裂隙发育及破碎飞石运动规律分析方法如：

cn111879858a公开了一种基于声发射震源主频唯一性的岩石破坏预测方法，包括标准岩石试件的制备，声发射震源定位，声发射震源类型辨识，不同类型震源对应的声发射主频的计算与岩石破坏预测方法。本发明解决了目前基于声发射主频特征的岩石破坏预测方法，存在声发射主频带宽，声发射主频特征不显著的缺陷。

cn111767631a公开了一种基于多相数字岩心模拟岩石裂纹扩展的方法及系统，包括获取岩心样本的三维灰度图像；根据多个灰度阈值范围对三维灰度图像中每个像素点进行分类、编号，得到每个编号对应的像素点集合；根据编号对应的像素点集合建立初始数字岩心数据库；初始数字岩心数据库包括每个像素点对应的像素点编号和像素点位置；根据初始数字岩心数据库，采用近场动力学模型进行岩心裂纹扩展模拟，建立破裂后的数字岩心数据库；对破裂后的数字岩心数据库进行可视化，得到岩心裂纹扩展过程的模拟结果。

上述现有技术中的岩体裂隙发育及飞石运动规律分析方法在一定程度上均可以对岩体裂隙产生、扩展和贯通进行分析或预测，但其还存在一定局限性，例如方法主体主要面向岩体破坏过程的数值模拟分析，以及实验室分析，缺乏现场监测的实用性，且实验过程较为繁琐，对仪器精度有极大的依赖性。

技术实现要素：

本发明的目的之一在于提供一种基于图像处理技术的岩体裂隙及飞石运动分析系统，该系统通过对岩体裂隙面进行实时的监测，并采用图像处理技术技术对岩石裂隙及飞石运动的检测数据进行准确分析和处理。

本发明的目的之二在于提供一种基于图像处理技术的岩体裂隙及飞石运动分析系统的方法。

本发明提供的这种基于图像处理技术的岩体裂隙及飞石运动分析系统，包括裂隙及周边环境监测装置1、第一数据传输线9、计算机处理系统11和数据反应系统13，裂隙及周边环境监测装置1通过第一数据传输线9连接计算机处理系统11，计算机处理系统11连接数据反应系统13；裂隙及周边环境监测装置1用于监测岩体裂隙与周围环境，并将监测数据发送到计算机处理系统11中；计算机处理系统11用于对裂隙及周边环境监测装置1输出的数据进行处理，获得岩体裂隙发育规律，以及岩体裂隙发育完全发生破碎后飞石运动规律，并将数据发送到数据反应系统13；数据反应系统13用于对计算机处理系统11处理后的数据产生动作。

所述的裂隙及周围环境监测装置1包括摄像头2和第一信号接收装置7，摄像头2连接第一信号接收装置7，第一信号接收装置7连接计算机处理系统11；摄像头2对外部的岩体裂隙面3、第一岩体裂隙4、第二岩体裂隙5和破碎飞石6进行实时的图像监测，将得到的图像数据通过第一信号接收装置7和第一数据传输线9传输到计算机处理系统11。

所述的计算机处理系统11包括计算机12，计算机12分别连接第一数据传输线9和数据反应系统13，计算机12通过图像处理技术，对裂隙及周围环境监测装置1输入的图像数据进行处理，并将处理后的参数输出到数据反应系统13。

所述的数据反应系统13包括操作台14、数据显示台15、第二信号接收装置8、第二数据传输线10和报警装置17；数据显示台15用于接收计算机12输出的数据并显示，同时当输出数据超过设定阈值之后，采用操作台14将控制信号通过第二数据传输线10传输到第二信号接收装置8，从而控制报警装置17进行紧急报警。

本发明还提供了一种基于图像处理技术的岩体裂隙及飞石运动分析系统的方法，包括如下步骤：

s1.通过摄像头2采集原色图，采集内容为岩体裂隙面3岩体裂隙发育，以及裂隙发育完全后岩体发生破坏产生的破碎飞石6的运动状态数据；

s2.将采集到的原色图转换成灰度图，并生成灰度值；

s3.对步骤s2中转换成的灰度图，进行二值化处理；

s4.定义结构元素，对步骤s3二值化处理后的灰度图像进行图像开运算去噪；

s5.对图像中的裂隙边界进行检测，绘制裂隙轮廓，并将得到的裂隙轮廓大小表示为裂隙的大小，将得到的裂隙轮廓数量表示为裂隙的数量；通过飞石运动规律分析方法处理飞石运动状态数据；

s6.对步骤s5中裂隙的数量及大小和飞石运动进行评价，并将其与设定阈值标准进行对比，将对比数据输出至数据反应系统13，并进行报警。

步骤s2，具体包括将采集到的原色图进行加权平均灰度化：

gray＝0.2989*r+0.5870*g+0.1140*b

其中，r代表图像红色通道的像素值；g代表绿色通道的像素值；b代表蓝色通道的像素值，gray为图像经灰度化后像素点的灰度值。

步骤s3具体为采用局部阈值的大津阈值分割算法，对灰度化处理后的图像进行二值化处理；假设灰度化图像像素数为n，灰度范围为[0,k-1]，灰度级i的像素数为ni，灰度级i出现的概率pi为：

根据灰度图灰度阈值t，将图像像素分为a、b两部分，a的像素阈值范围为：[0,t]，同时b的像素阈值范围为[t+1,k-1]，根据像素灰度分布概率，求得：

像素区域a的灰度均值为：

像素区域b的灰度均值为：

其中，a的像素阈值范围为：[0,t]，同时b的像素阈值范围为[t+1,k-1]，pi为灰度级i出现的概率；n为灰度化图像像素数，灰度范围为[0,k-1]；

已知，灰度图像的灰度均值ut、像素区域a的灰度均值ua和像素区域b的灰度均值ub，可得灰度图像的灰度均值ut：

ut＝ω0ua+ω1ub

其中，a的像素阈值范围为：[0,t]，同时b的像素阈值范围为[t+1,k-1]，pi为灰度级i出现的概率；n为灰度化图像像素数，灰度范围为[0,k-1]；

像素区域a和b的方差为：

σ²＝ω0(ua-ut)²+ω1(ub-ut)²

其中，a的像素阈值范围为：[0,t]，同时b的像素阈值范围为[t+1,k-1]，pi为灰度级i出现的概率；n为灰度化图像像素数，灰度范围为[0,k-1]；ut为灰度图像的灰度均值；ua为像素区域a的灰度均值；ub为像素区域b的灰度均值；

按照最大类间方差，灰度阈值t在[0,k-1]范围中依次取值，当阈值t取值使σ²取得最大值时，则阈值t即为算法的最佳阈值。

步骤s4，采用的开运算去噪具体体现为对灰度图像二值化后，先腐蚀，再膨胀的运算方法；

假设待处理的二值化图像像素集合为i，结构元素为p：

p对i的腐蚀运算具体定义为：

其中，θ为腐蚀算子，x∈i

p对i的膨胀运算定义为：

其中，为膨胀算子；

结构元素集合p在二值化处理后的灰度图的集合i内进行平移并集，开运算具体定义为：

其中，θ为腐蚀算子，为膨胀算子；

具体表现为对s4输出的图像，以从左到右、从上到下的顺序，进行像素值扫描，若扫描到像素值为255时，对目标轮廓进行检测框绘制，并针对检测框对角线方向，只保留检测方向的终点坐标，并对检测框进行数量标定；否则继续扫描；

将实际裂隙面积大小设为：

y1＝k1x1+b1

其中，k1和b1为常数；y1为裂隙实际面积大小；x1为图像中裂隙面积大小；实际裂隙及飞石面积大小方程的确定方法为：在实际情况下，对裂隙绘制轮廓，测算轮廓面积大小；将实际测算的裂隙轮廓面积大小与在图像上得出的面积数据进行多次拟合，得出k1，b1值，最后取平均值，得出对应方程；

将飞石面积大小方程设为：

y'1＝k'1x'1+b'1

其中，k’1和b’1为常数；y’1为飞石实际面积大小；x’1为图像中飞石面积大小；实际飞石面积大小方程的确定方法为：在实际情况下，对飞石绘制轮廓，测算轮廓面积大小；将实际测算的飞石轮廓面积大小与在图像上得出的面积数据进行多次拟合，得出k’1，b’1值，最后取平均值，得出对应方程；

将实际裂隙密度大小方程设为：

y2＝k2x2+b2

其中，k2和b2为常数；y2为裂隙实际密度大小；x2为图像中裂隙密度大小；实际裂隙密度大小方程的确定方法为：在实际情况下，对裂隙绘制轮廓，测算轮廓面积大小，在摄像头检测范围内，除以裂隙面面积，将实际裂隙密度与在图像上得出的数据进行多次拟合，得出k2，b2值，最后取平均值，得出密度确定方程；

将飞石密度大小方程设为：

y'2＝k'2x'2+b'2

其中，k’2和b’2为常数；y’2为飞石实际密度大小；x2为图像中飞石密度大小；实际飞石密度大小方程的确定方法为：在实际情况下，对飞石绘制轮廓，测算轮廓面积大小，在摄像头检测范围内，除以飞石分布空间体积，将实际的飞石密度与在图像上得出的数据进行多次拟合，得出k’2，b’2值，最后取平均值，得出密度确定方程。

步骤s5中，飞石运动规律分析方法，包括：飞石的体积和质量分析、飞石的初始动能和势能分析；

其中，v为飞石体积；y1为裂隙或飞石实际面积大小；ρ为岩石已知密度；

监测图像中飞石初速度

δx代表图像中最初两帧时间差内，飞石运动坐标方向上的位移之差；t代表图像两帧的时间差；

实际飞石速度大小方程设为：

y3＝k3x3+b3

其中，k3和b3为常数；y3为飞石实际速度大小，x3为监测图像中飞石初速度v；该方程的确定方法为：在实际情况下，图像采集帧间时间差中，测得飞石运动实际距离，从而得出对应实际速度，采集若干组数据，将采集到的若干组数据与在图像上得出的数据进行若干次拟合，得出k3，b3值，最后取平均值，得出飞石实际速度确定方程；

已知监测岩体的高度h，动能、势能参数具体表现为：

ep为飞石的初始动能、ek为飞石的初始势能，m代表飞石的质量，g代表重力加速度。

本发明提供的这种基于图像处理技术的岩体裂隙及飞石运动分析系统及方法，通过裂隙检测装置对岩体裂隙面进行实时的监测，将监测数据以图像数据的形式传输给高性能计算机，计算机利用图像处理技术，对图像进行灰度化处理、二值化、去噪、轮廓扫描、检测框绘制等，得到裂隙的数量、大小参数，以及裂隙发育完全后产生破碎飞石数量、大小、密度的基本参数，初速度、初始动能、势能等宏观层面上的运动学参数。将得到的参数与设定阈值标准进行比较，从而在裂隙发育未完全，有岩爆、地表塌陷等事故危险情况下，起到一定的预警作用，有利于保障仪器设备、人员安全，促进工程建设的顺利正常进行。

附图说明

图1为本发明系统的结构示意图。

图2为本发明实施例的结构示意图。

图3为本发明方法的流程示意图。

图4为本发明实施例的观测结果示意图。

图5为本发明实施例的opencv处理示意图。

图6为本发明实施例的飞石数量和大小示意图。

图7为本发明实施例的标记图片示意图。

具体实施方式

如图1为本发明系统的结构示意图；本发明提供的这种基于图像处理技术的岩体裂隙及飞石运动分析系统，包括裂隙及周边环境监测装置1、第一数据传输线9、计算机处理系统11和数据反应系统13，裂隙及周边环境监测装置1通过第一数据传输线9连接计算机处理系统11，计算机处理系统11连接数据反应系统13；裂隙及周边环境监测装置1用于监测岩体裂隙与周围环境，并将监测数据发送到计算机处理系统11中；计算机处理系统11用于对裂隙及周边环境监测装置1输出的数据进行处理，获得岩体裂隙发育规律，以及岩体裂隙发育完全发生破碎后飞石运动规律，并将数据发送到数据反应系统13；数据反应系统13用于对计算机处理系统11处理后的数据产生动作。

如图2为本发明实施例的结构示意图；如图2所示1为裂隙及周围环境监测装置；2为摄像头(在本实施例中为高速运动摄像头)；3为岩体裂隙面；4为第一岩体裂隙；5为第二岩体裂隙；6为破碎飞石；7为第一信号接收装置；8为第二信号接收装置；9为第一数据传输线；10为第二数据传输线；11为计算机处理系统；12为计算机；13为数据反应系统；14为操作台；15为数据显示台；16为地平面；17为报警装置。

裂隙及周围环境监测装置1包括摄像头2和第一信号接收装置7，摄像头2连接第一信号接收装置7，第一信号接收装置7连接计算机处理系统11；摄像头2能够自由转向，对外部的岩体裂隙面3、第一岩体裂隙4、第二岩体裂隙5和破碎飞石6进行实时的图像监测，以每秒为单位时间间隔，将得到的图像数据通过第一信号接收装置7和第一数据传输线9传输到计算机处理系统11；计算机处理系统11包括计算机12，计算机12分别连接第一数据传输线9和数据反应系统13，计算机12通过图像处理技术，对裂隙及周围环境监测装置1输入的图像数据进行处理，并将处理后的参数输出到数据反应系统13；数据反应系统13包括操作台14、数据显示台15、第二数据传输线10、第二信号接收装置8和报警装置17；数据显示台15用于接收计算机12输出的数据并显示，同时当输出数据超过设定阈值之后，采用操作台14将控制信号通过第二数据传输线10传输到第二信号接收装置8，从而控制报警装置17进行紧急报警。在本实施例中，警铃工作，同时进行灯光和钟鸣警示，指导安全撤离工作。

如图3为本发明方法的流程示意图：本发明提供的这种基于图像处理技术的岩体裂隙及飞石运动分析系统的方法，包括如下步骤：

s1.通过摄像头2采集原色图，采集内容为岩体裂隙面3岩体裂隙发育，以及可能产生的，裂隙发育完全后岩体发生破坏产生的破碎飞石6的运动状态数据；

s2.将采集到的原色图转换成灰度图，并生成灰度值；

s3.对步骤s2中对原图进行灰度化处理后的灰度图，进行二值化处理，将整个图像呈现出明显的黑白效果，达到凸显出目标的轮廓的效果；

s4.定义结构元素，对步骤s3二值化处理后的灰度图像进行图像开运算去噪；

s5.对图像中的裂隙边界进行检测，绘制裂隙轮廓，并将得到的裂隙轮廓大小表示为裂隙的大小，将得到的裂隙轮廓数量表示为裂隙的数量；通过飞石运动规律分析方法处理飞石运动状态数据；

当裂隙发育完全，岩体发生破碎时，对图像中产生的飞石边界进行检测，同时绘制飞石轮廓，轮廓数量表示为飞石数量；建立线性方程，裂隙和飞石的轮廓大小与实际对应目标的面积关系，同时，利用飞石数量及所得图像面积，通过建立另一线性方程解算得出图像密度与实际浓度之间的联系。通过飞石运动轨迹得到飞石运动初速度、动能、势能等运动学参数。

s6.对步骤s5中裂隙的数量及大小进行评价，并将其与设定阈值标准进行对比，将对比数据输出至数据反应系统13，若任一指标超出设定阈值，则进行报警，若裂隙的发育状态和飞石运动状态均在在设定阈值内，则继续监测，直到裂隙变大或飞石过多等情况。

步骤s2具体包括将采集到的原色图进行加权平均灰度化：

gray＝0.2989*r+0.5870*g+0.1140*b

其中，r代表图像红色通道的像素值；g代表绿色通道的像素值；b代表蓝色通道的像素值，gray为图像经灰度化后像素点的灰度值。

步骤s3具体为采用局部阈值的大津阈值分割算法，对灰度化处理后的图像进行二值化处理；假设灰度化图像像素数为n，灰度范围为[0,k-1]，灰度级i的像素数为ni，灰度级i出现的概率pi为：

根据灰度图灰度阈值t，将图像像素分为a、b两部分，a的像素阈值范围为：[0,t]，同时b的像素阈值范围为[t+1,k-1]，根据像素灰度分布概率，得到：

像素区域a的灰度均值为：

像素区域b的灰度均值为

其中，a的像素阈值范围为：[0,t]，同时b的像素阈值范围为[t+1,k-1]，pi为灰度级i出现的概率；n为灰度化图像像素数，灰度范围为[0,k-1]；

已知，灰度图像的灰度均值ut、像素区域a的灰度均值ua和像素区域b的灰度均值ub，可得灰度图像的灰度均值ut为：

ut＝ω0ua+ω1ub

其中，a的像素阈值范围为：[0,t]，同时b的像素阈值范围为[t+1,k-1]，pi为灰度级i出现的概率；n为灰度化图像像素数，灰度范围为[0,k-1]；

像素区域a和b的方差为：

σ²＝ω0(ua-ut)²+ω1(ub-ut)²

其中，a的像素阈值范围为：[0,t]，同时b的像素阈值范围为[t+1,k-1]，pi为灰度级i出现的概率；n为灰度化图像像素数，灰度范围为[0,k-1]；ut为灰度图像的灰度均值；ua为像素区域a的灰度均值；ub为像素区域b的灰度均值。

按照最大类间方差，灰度阈值t在[0,k-1]范围中依次取值，当阈值t取值使σ²取得最大值时，则阈值t即为算法的最佳阈值。

步骤s4，采用的开运算去噪具体体现为对灰度图像二值化后，先腐蚀，再膨胀的运算方法；

假设待处理的二值化图像像素集合为i，结构元素为p：

p对i的腐蚀运算具体定义为：

其中，θ为腐蚀算子，x∈i

p对i的膨胀运算定义为：

其中，为膨胀算子。

结构元素集合p在二值化处理后的灰度图的集合i内进行平移并集，开运算具体定义为：

其中，θ为腐蚀算子，为膨胀算子；

具体表现为对s4输出的图像，以从左到右、从上到下的顺序，进行像素值扫描，若扫描到像素值为255时，对目标轮廓进行检测框绘制，并针对检测框对角线方向，只保留检测方向的终点坐标，并对检测框进行数量标定；否则继续扫描。

将实际裂隙面积大小设为：

y1＝k1x1+b1

其中，k1和b1为常数；y1为裂隙实际面积大小；x1为图像中裂隙面积大小；实际裂隙及飞石面积大小方程的确定方法为：在实际情况下，对裂隙绘制轮廓，测算轮廓面积大小；将实际测算的裂隙轮廓面积大小与在图像上得出的面积数据进行多次拟合，得出k1，b1值，最后取平均值，得出对应方程；

将飞石面积大小方程设为：

y'1＝k'1x'1+b'1

其中，k’1和b’1为常数；y’1为飞石实际面积大小；x’1为图像中飞石面积大小；实际飞石面积大小方程的确定方法为：在实际情况下，对飞石绘制轮廓，测算轮廓面积大小；将实际测算的飞石轮廓面积大小与在图像上得出的面积数据进行多次拟合，得出k’1，b’1值，最后取平均值，得出对应方程；

将实际裂隙密度大小方程设为：

y2＝k2x2+b2

其中，k2和b2为常数；y2为裂隙实际密度大小；x2为图像中裂隙密度大小；实际裂隙密度大小方程的确定方法为：在实际情况下，对裂隙绘制轮廓，测算轮廓面积大小，在摄像头检测范围内，除以裂隙面面积，将实际裂隙密度与在图像上得出的数据进行多次拟合，得出k2，b2值，最后取平均值，得出密度确定方程。

将飞石密度大小方程设为：

y'2＝k'2x'2+b'2

其中，k’2和b’2为常数；y’2为飞石实际密度大小；x2为图像中飞石密度大小；实际飞石密度大小方程的确定方法为：在实际情况下，对飞石绘制轮廓，测算轮廓面积大小，在摄像头检测范围内，除以飞石分布空间体积，将实际的飞石密度与在图像上得出的数据进行多次拟合，得出k’2，b’2值，最后取平均值，得出密度确定方程。

步骤s5中，飞石运动规律分析方法，包括：飞石的体积和质量分析；飞石的初始动能和势能分析。

其中，v为飞石体积；y1为裂隙或飞石实际面积大小；ρ为岩石已知密度；

监测图像中飞石初速度v：

δx代表图像中最初两帧时间差内，飞石运动坐标方向上的位移之差；t代表图像两帧的时间差，在本实施例中为1/12s。

实际飞石速度大小方程设为：

y3＝k3x3+b3

其中，k3和b3为常数；y3为飞石实际速度大小，x3为监测图像中飞石初速度v。该方程的确定方法为：在实际情况下，图像采集帧间时间差中，测得飞石运动实际距离，从而得出对应实际速度，采集多组数据，将其与在图像上得出的数据进行多次拟合，得出k3，b3值，最后取平均值，得出飞石实际速度确定方程。

已知监测岩体的高度h，动能、势能参数具体表现为：

ep为飞石的初始动能、ek为飞石的初始势能，m代表飞石的质量，g代表重力加速度，在本实施例中为9.80m/s²。

在具体实施方式中，以某铁路隧道作为试验基地。该隧道是由某单位承建的国外铁路隧道项目，共发生过约3000次岩爆崩落事故，最严重的一次崩落长度甚至达24米，崩落方量约为1500m³，频繁、难以预料的岩爆事故不仅严重影响工程建设正常进度(影响工期长达3个月以上)，大大增加了施工成本，还带来了巨大的人员、设备的安全风险。鉴于此，在该隧道3号在建隧道段进行本次岩体裂隙发育及飞石运动规律监测分析。试验步骤如下：

步骤1、准备硬件设施，包括高速运动摄像机2、信号接收装置7，8、高性能计算机12、数据显示台15、报警装置17。

步骤2、高速运动摄像头2输入采集到的原色图，采集内容为裂隙面3岩体裂隙发育，以及可能产生的，裂隙发育完全后岩体发生破坏产生的飞石6运动；

步骤3、将采集到的原色图转换成灰度图，并生成灰度值；

步骤4、在步骤3中对原图进行灰度化处理后，进行二值化处理，将整个图像呈现出明显的黑白效果，达到凸显出目标的轮廓的效果；

步骤5、定义结构元素，对步骤4二值化处理后的灰度图像进行图像开运算去噪；

步骤6、对图像中的裂隙边界进行检测，绘制裂隙轮廓，并将得到的裂隙轮廓大小表示为裂隙大小，将得到的裂隙轮廓数量表示为岩体裂隙的数量；

当裂隙发育完全，岩体发生破碎时，对图像中产生的飞石边界进行检测，同时绘制飞石轮廓，轮廓数量表示为飞石数量；建立线性方程，裂隙和飞石的轮廓大小与实际对应目标的面积关系，同时，利用飞石数量及所得图像面积，通过建立另一线性方程解算得出图像密度与实际浓度之间的联系。通过飞石运动轨迹得到飞石运动初速度、动能、势能等运动学参数。

步骤7、对步骤6中裂隙发育数量及大小和飞石运动情况进行评价，并将其与设定阈值标准进行对比，根据对比结果将数据输出至数据反应系统13中的数据显示台15。

步骤8、当裂隙面裂隙发育的数量、大小任一参数超出设定阈值时，报警装置17发生动作，进行实时紧急报警，班组长立即要求所有施工人员撤离施工现场。

对该隧道3号建设段运用本方法监测裂隙面裂隙发育情况，观测结果如图4所示，如图4为本发明实施例的观测结果示意图；经计算机通过opencv计算对图4进行处理后，结果如图5为本发明实施例的opencv处理示意图。当图5中裂隙数量或大小超出设定的阈值标准时，进行报警，并有序撤离了施工人员，短时间后，隧道发生岩体崩落，如图6为本发明实施例的飞石数量和大小示意图。，经计算机处理后，输出图6崩落的飞石进行数量、大小。并可在后期进行飞石密度和动力学规律的详细分析。

如图7为本发明实施例的标记图片示意图，从图5、图7可知，在该隧道3号建设段运用一种基于图像处理技术的岩体裂隙发育及飞石运动规律分析方法，计算机成功对岩体裂隙面进行裂隙、破碎飞石的实时监测，并通过opencv技术，显现出了裂隙、飞石的数量及其标记，同时在图片上记录了相应的大小和出现的坐标。该裂隙及飞石监测过程符合工作岩体实际情况。