[发明专利]一种基于破碎岩体模型的地铁车站稳定性分析方法

权利要求书

1.一种基于破碎岩体模型的地铁车站稳定性分析方法，其特征在于，包括：

S1：安装拍摄设备，对掌子面进行拍摄，获取破碎岩体区域的掌子面图像；

S2：对所述掌子面图像进行多边形节理识别，以统计多边形节理的边长的长度值，并计算出所述多边形节理的边长的长度值的平均值；

S3：根据所述掌子面图像，布置种子点并进行筛选，生成节理分布平面图；

S4：处理所述节理分布平面图，以将所述节理分布平面图中多边形节理的边的端点坐标输出为命令流格式的文件；

S5：根据所述节理分布平面图中的多边形节理的边的端点坐标，构建包含多边形节理的边的端点坐标的地铁车站施工模型；

S6：根据所述地铁车站施工模型，采用对节理和岩块的抗剪强度和抗拉强度同时进行强度折减的方法，计算地铁车站施工的安全系数。

2.根据权利要求1所述的一种基于破碎岩体模型的地铁车站稳定性分析方法，其特征在于，所述S1中包括如下步骤：

S11：在掌子面前方地面上，安装拍摄底座(100)，使所述拍摄底座(100)平稳安置于地面上；

S12：将支撑杆(200)垂直固定于所述拍摄底座(100)顶部，拍摄设备(300)固定安装在支撑杆(200)顶部；

S13：将所述拍摄设备(300)与配电箱(400)进行电连接；

S14：调节所述拍摄设备(300)，对隧道开挖形成的掌子面进行拍摄，获取破碎岩体区域的掌子面图像。

3.根据权利要求1所述的一种基于破碎岩体模型的地铁车站稳定性分析方法，其特征在于，所述S2中包括如下步骤：

S2：对所述掌子面图像进行多边形节理识别，以统计多边形节理的边长的长度值，并计算出所述多边形节理的边长的长度值的平均值；

S21：对所述掌子面图片进行比例转换，以获取结构面图像；

所述比例转换系数公式如下：

其中s为图像的比例转换系数，(x₁,y₁)、(x₂,y₂)分别为所拍摄的掌子面图像上的两点的像素坐标，H为两点间的实际距离；

S22：对所述结构面图像进行预处理，以获取处理后的结构面图像；

S23：对所述处理后的结构面图像的边缘信息进行识别，以获取边缘检测图像；

S231:所述处理后的结构面图像表示为：

式中：N代表处理后的结构面图像中的像素点的数量；f(x,y)代表处理后的结构面图像；x代表处理后的结构面图像的横坐标；y代表处理后的结构面图像的纵坐标；

S232:将所述处理后的结构面图像与二维高斯滤波函数进行卷积，获取滤波后的图像；所述二维高斯滤波函数如下：

式中：σ²代表方差；

将所述高斯滤波器分解为X方向的滤波器和Y方向的滤波器，分别如下：

式中：k是计算的中间参数；

S233：将所述滤波后的图像分别与所述X方向的滤波器和Y方向的滤波器进行卷积，获取X方向的梯度和Y方向的梯度：：

式中：B_x表示X方向的梯度；B_y表示Y方向的梯度；

S234：令：

式中：A(i,j)表示处理后的结构面图像f(x,y)上的点M(i,j)处的边缘强度；α(i,j)表示处理后的结构面图像f(x,y)上的点M(i,j)处的法向矢量，M(i,j)为处理后的结构面图像上的一个像素点；i表示点M(i,j)的横坐标；j表示点M(i,j)的纵坐标；

S235：对A(i,j)进行非极大值抑制处理，以获取检测边缘点；

得到处理后的结构面图像f(x,y)上的点M(i,j)处的边缘强度即A(i,j)后，进行非极大值抑制处理，即将点M(i,j)的边缘强度与梯度方向的相邻像素点的边缘强度进行比较，判断点M(i,j)的边缘强度是否是局部极大值，若是则M＝M(i,j)记为检测边缘点，若不是将其边缘强度设为0；

S236：利用双阀值算法处理检测边缘点并进行边缘连接，获取边缘检测图像；

S24：对所述边缘检测图像中多边形节理的边长进行检测，以统计所述边缘检测图像中所有的多边形节理的边长的长度值和所述边缘检测图像中各类多边形所占比例，并计算出所述边缘检测图像中所有的多边形节理的边长的长度值的平均值：

式中，L为节理长度，为节理平均长度，n为直线段总数。