一种隧道爆破施工超欠挖优化方法

文献发布时间：2024-04-18 19:58:30

技术领域

本发明涉及隧道施工领域，特别是一种隧道爆破施工超欠挖优化方法。

背景技术

隧道工程中，钻爆法是一种常用的施工方法，但在实际施工中往往出现爆破效果不良和超欠挖现象的问题。其中的主要原因是隧道不同部位的围岩结构面分布情况不同，导致传统爆破方案无法满足各部位的要求，爆破施工常常受到围岩结构面分布不均的影响，使得爆破效果不一致，进而引发超欠挖现象。如何优化隧道爆破施工过程中的超欠挖问题成了钻爆施工设计方案中亟待解决的重要问题。为了解决这一问题，现有技术通常无法准确选择适应不同结构面分布的围岩的爆破参数，从而无法最小化超欠挖的程度。

发明内容

本发明的目的在于克服现有技术的缺点，提供一种隧道爆破施工超欠挖优化方法。本发明的目的通过以下技术方案来实现：一种隧道爆破施工超欠挖优化方法，包括以下步骤：S1、选择初始爆破参数集；

S2、划分岩体区域，获取各岩体区域的岩体质量评估参数；

S3、搭建输入参数与输出参数的映射关系模型，输入参数为岩体质量评估参数和初始爆破参数，输出参数包括超欠挖面积参数、炮孔利用率、拱顶下沉参数和块石尺寸参数；

S4、基于输入参数对所述映射关系模型进行训练直至稳定；

S5、基于光面爆破参数范围表并利用寻优算法获取最优光面爆破参数并输出。

具体的，S1中所述的初始爆破参数包括炮孔数量、孔深、孔径、孔距、总药量和炸药密度。具体的，S2中划分岩体区域时将挖面划分成六个小块的区域，包括上中下三层，每层两块区域，分别对应拱顶至拱肩、拱肩至拱腰以及拱腰至拱底区域。

具体的，S2中岩体质量评估参数包括结构面信息参数和围岩质量指标，所述结构面信息参数包括倾向、倾角、走向、迹长和体积节理数；围岩质量指标采用YZ岩体质量评价体系，按下式计算YZ值：

YZ＝90+3δ

式中，δ

式中，T

具体的，S3中所述映射关系模型包括输入层、隐含层和输出层，并通过均方误差损失函数对所述映射关系模型进行优化，其中，所述输入层、隐含层和输出层之间的激活函数如下公式：

h＝max(0,a

其中，h为模型基本处理单元的输出，a

其中，s

具体的，S4中对于采集的输入参数，90％数据用于训练，10％的数据用于预测，映射关系模型预测的决定系数、均方误差和均方根误差作为评价指标，当满足映射关系模型要求误差条件时认为映射关系模型已经训练完成，然后绘制出真实值与预测值的图像以及真实值与预测值的误差图，其中，所述映射关系模型的适应度函数按下式进行计算：

其中，D为测试样本数量；b

具体的，S5中先设置超欠挖优化算法，超欠挖优化算法搜索过程始于一个候选解，其中每个爆破参数被赋予随机的初始位置，爆破参数为超欠挖面积参数、炮孔利用率、拱顶下沉参数和块石尺寸参数，在评估候选解的适应度后，超欠挖优化算法开始迭代直到满足终止条件，将超欠挖优化算法代入训练好的映射关系模型，计算出超欠挖面积参数、炮孔利用率、拱顶下沉参数和块石尺寸参数，通过寻优算法寻找出满足适应度值最小的光面爆破参数，输出最优光面爆破参数，其中寻优过程中每个候选解的适应度值由下式计算：

其中，N为光面爆破参数初始群体值数量，e

具体的，寻优算法优化的步骤为：

S10、初始化阶段；

给定搜索空间，超欠挖优化算法通过在d维搜索空间中随机初始化若干个候选解的位置

来启动，如下式：

其中，W是所有候选解的位置，d代表维度，

候选解的初始位置按下式随机生成：

其中，x

然后插入一个用户定义的适应度评估函数，利用适应度评估函数对每个候选解的位置进

行评估，相应的适应度值被存储在数组中如下式所示：

将所述适应度值进行初始化并排序，接着初始化全局最佳位置和最佳位置，将排序后的位置中的第一个值赋给全局最佳位置，将排序后的适应度中的第一个值赋给最佳位置；S20、超欠挖优化算法初始化参数：

两个超欠挖优化算法中的基本参数按下式给出：

Gg＝0.1×log(2.75×(t/T)

式中，t代表当前迭代次数，T代表迭代总次数。

S30、超欠挖优化算法的迭代过程：

循环迭代T次，进行优化过程，同时生成随机候选解记为a(i)，按下列三种情况计算候选

解新位置：

Case1:

Case2:

Case3:

其中，

S40、在每次更新候选解位置后，检查当前候选解的适应度是否比全局最佳候选解位置的适应度更好，若更好，则更新全局最佳候选解位置为当前候选解的位置，经过所有候选解的迭代后，全局最佳候选解即为所求最优爆破参数。

本发明具有以下优点：

本发明的映射关系模型，用于生成爆破参数的范围并优选最佳参数，以有效减少隧道的超欠挖，利用现代隧道技术和电子设备获取实时的隧道围岩结构面信息，并根据这些信息生成爆破参数的范围。通过分析围岩结构面在隧道拱顶、拱肩、拱腰、拱底等不同部位的分布情况，能够确定不同部位所需的爆破参数，并进行优选，以使超欠挖的程度最小化。这种非均匀布控的爆破参数选择方法能够解决围岩不均匀分布导致爆破效果不良的问题，并显著减少超欠挖现象的发生。

通过本发明，隧道施工方能更加准确地选择合适的爆破参数，从而有效减少超欠挖现象的发生。这不仅降低了隧道施工过程中的排岩量和初期支护工程量，减少了建设成本，还能避免超欠挖带来的安全隐患，提高了隧道的长期使用安全性和稳定性。

附图说明

图1为本发明的隧道岩体区域分块示意图；

图2为本发明的真实值与预测值训练的拟合示意图；

图3为本发明的真实值与预测值验证的拟合示意图；

图4为本发明的真实值与预测值测试的拟合示意图；

图5为本发明的真实值与预测值的拟合示意图；

图6为本发明的的真实值与预测值的误差示意图。

具体实施方式

下面结合附图对本发明做进一步的描述，但本发明的保护范围不局限于以下所述。如图1-图6所示，一种隧道爆破施工超欠挖优化方法，包括以下步骤：

S1、选择初始爆破参数集，所述的初始爆破参数包括炮孔数量、孔深、孔径、孔距、总药量和炸药密度；

S2、划分岩体区域，将挖面划分成六个小块的区域，包括上中下三层，每层两块区域，分别对应拱顶至拱肩、拱肩至拱腰以及拱腰至拱底区域，获取各岩体区域的岩体质量评估参数；岩体质量评估参数包括结构面信息参数和围岩质量指标，所述结构面信息参数包括倾向、倾角、走向、迹长和体积节理数，本发明采用数字摄影测量技术获取结构面信息参数，其从岩体不同角度对指定区域进行数字成像，利用三维重构技术建立岩体表面三维模型，经后期处理获取岩体结构面信息；围岩质量指标采用YZ岩体质量评价体系，按下式计算YZ值：

YZ＝90+3δ

式中，δ

式中，T

S3、搭建输入参数与输出参数的映射关系模型，输入参数为岩体质量评估参数和初始爆破参数，输出参数包括超欠挖面积参数、炮孔利用率、拱顶下沉参数和块石尺寸参数；所述映射关系模型包括输入层、隐含层和输出层，并通过均方误差损失函数对所述映射关系模型进行优化，其中，所述输入层、隐含层和输出层之间的激活函数如下公式：

h＝max(0,a

其中，h为模型基本处理单元的输出，a

其中，s

S4、基于输入参数对所述映射关系模型进行训练直至稳定；对于采集的输入参数，90％数据用于训练，10％的数据用于预测，映射关系模型预测的决定系数、均方误差和均方根误差作为评价指标，当满足映射关系模型要求误差条件时认为映射关系模型已经训练完成，然后绘制出真实值与预测值的图像以及真实值与预测值的误差图，如图2和图3所示，其中，所述映射关系模型的适应度函数按下式进行计算：

其中，D为测试样本数量；b

S5、基于光面爆破参数范围表并利用寻优算法获取最优光面爆破参数并输出；在隧道施工前需对围岩情况进行地质勘探，获取围岩类型以及围岩的坚硬程度，从而根据已有规范获取当前挖面围岩所对应的光面爆破参数范围表，如表1所示：

表1、光面爆破参数范围表

在计算机系统中设置超欠挖优化算法，超欠挖优化算法搜索过程始于一个候选解，其中每个爆破参数被赋予随机的初始位置，爆破参数为超欠挖面积参数、炮孔利用率、拱顶下沉参数和块石尺寸参数，在评估候选解的适应度后，超欠挖优化算法开始迭代直到满足终止条件，在每次迭代中，超欠挖优化算法会更新每个爆破参数的位置，这样，超欠挖优化算法在每次迭代中通过更新位置来寻求更好的解决方案，通过迭代的过程，超欠挖优化算法不断优化每个爆破参数的位置，以期最终找到全局最优解或接近最优解的解决方案。

将超欠挖优化算法代入上步训练好的映射关系模型，计算出超欠挖面积参数、炮孔利用率、拱顶下沉参数和块石尺寸参数，通过寻优算法寻找出满足适应度值最小的光面爆破参数，输出最优光面爆破参数，其中寻优过程中每个候选解的适应度值由下式计算：

其中，N为光面爆破参数初始群体值数量。e

需要说明的是，每一种岩石都有相对应的拱顶下沉参数范围，在寻优过程中首先判断拱顶下沉参数是否满足要求，然后计算候选解的适应度值，若拱顶下沉参数不符合要求则在适应度值的计算过程中适当的添加一惩罚条件使其值偏大。

利用寻优算法优化的具体步骤如下：

S10、初始化阶段：

给定搜索空间，超欠挖优化算法是通过在d维搜索空间中随机初始化若干个候选解的位置来启动的，如下式：

其中，W是所有候选解的位置，d代表维度，

候选解的初始位置按下式随机生成：

其中，l

插入一个用户定义的适应度评估函数，利用适应度评估函数对每个候选解的位置进行评估；相应的适应度值被存储在数组中如下所示：

将所述适应度值进行初始化并排序，接着初始化全局最佳位置和最佳位置。将排序后的位置中的第一个值赋给全局最佳位置，将排序后的适应度中的第一个值赋给最佳位置。

S20、优化算法初始化参数：