随机森林融合SVM预测盾构引起近接结构变形方法及装置与流程

本发明属于隧道盾构施工诱发临近建筑物变形沉降的智能算法预测领域，具体涉及一种随机森林融合svm预测盾构引起近接结构变形方法及装置。

背景技术：

近些年来，随着我国社会经济快速发展，城市现代化建设程度不断提高，城市地下空间开发与利用成为重点发展方向，相应地各大城市地铁隧道的施工规模在飞速扩大。地铁隧道的开挖,使地层产生了变形和应力重分布,不可避免会导致地表的不均匀沉降,尤其对既有建筑物而言,往往很小的变形沉降,会带来较大的危害,影响结构的安全稳定性和使用功能。因此,预测城市地铁盾构施工对地表建筑物的沉降具有十分重要的意义。

目前，国内外学者研究的隧道盾构施工引起地表建筑物沉降方法主要有理论解析、实测分析法、有限单元法、智能算法等。其中用于预测地表建筑物沉降的智能算法一般有bp神经网络、支持向量机、遗传算法、灰色预测模型等。这些算法在对建筑物沉降的预测上有其可行性，但在实际运用过程中存在重大限制。例如bp神经网络存在易陷入局部极值和收敛速度慢等缺点，支持向量机在输入变量过多时，训练时间较长且对缺失值较敏感；遗传算法涉及到大数据处理时计算时间长，得到的结果稳定性差；灰色预测模型的精度很大程度依赖于假设的模型服从函数的正确性。

技术实现要素：

针对现有方法的以上缺陷或改进需求，本发明提供一种随机森林融合svm预测盾构引起近接结构变形方法及装置，用以解决在隧道开挖过程中诱发的地表建筑物变形沉降的快速预测问题。该方法利用随机森林特征选择，能在过多影响因素中，剔除与预测值相关性小的变量，筛选出用于建模的关键变量得到最优变量组合，降低支持向量机训练模型的维度，提高预测精度，得到更贴近实际的预测结果，实现了建筑物变形沉降更稳定、更准确的预测。

为了实现上述目的，按照本发明的一个方面，提供一种随机森林融合svm预测盾构引起近接结构变形方法，包括：

根据影响建筑物变形沉降的主要因素，并收集所述主要因素对应的数据；

建立随机森林模型，根据所述随机森林模型训练所述数据并进行重要性度量，得到最优特征集；

将所述最优特征集降维处理后输入支持向量机模型进行训练，得到建筑物沉降量预测结果。

进一步的，对所述建筑物沉降量预测结果进行误差分析，并进行全局敏感性评价。

进一步地，选择高斯核函数建立所述支持向量机模型，并采用交叉验证进行参数寻优。

进一步地，确定所述支持向量机模型的最优参数组合，基于所述最优特征集对样本训练集和测试集进行训练预测，输出所述训练集和所述测试集的实际值与预测值的拟合曲线，根据所述拟合曲线判断预测效果。

进一步地，根据下列公式计算一阶敏感度和总敏感度，

其中，si为参数的一阶敏感度值；vi为某个参数的方差；v为系统的总方差；k＝1，2，…ns；sti为参数的总敏感度值；ns为参数样本数，ns＝2m(wmax+1)；wi为参数特征频率；wmax为wi的最大值；m为4或6。

进一步的，所述影响建筑物变形沉降的主要因素包括，隧道相关因素、水文地质因素、盾构施工因素、建筑物因素。

进一步的，所述隧道相关因素包括隧道埋深和覆跨比。

进一步的，所述水文地质因素包括内摩擦角、弹性模量和粘聚力.

进一步的，所述盾构施工因素包括推进速度、刀盘扭矩、推进力、刀盘转速、上部土仓压力和注浆量。

进一步的，所述建筑物因素包括相对水平位置、相对垂直位置、相对纵向位置、建筑物自身结构和建筑物完好程度。

第二方面，本发明提供一种预测隧道盾构施工引起建筑物变形沉降的装置，其特征在于，所述装置包括：

获取模块，根据影响建筑物变形沉降的主要因素，并收集所述主要因素对应的数据；

处理模块，建立随机森林模型，根据所述随机森林模型训练所述数据并进行重要性度量，得到最优特征集；

确定模块，将所述最优特征集降维处理后输入支持向量机模型进行训练，得到建筑物沉降量预测结果。

第三方面，本发明提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现本发明第一方面中所述一种随机森林融合svm预测盾构引起近接结构变形方法。

第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现本发明第一方面中所述一种随机森林融合svm预测盾构引起近接结构变形方法。

总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：

1.本发明中提供的随机森林融合svm预测盾构引起近接结构变形方法充分利用了随机森林能够有效处理有多个指标、噪声干扰的样本数据，筛选有用特征变量的特点，创新性地提高了支持向量机模型的预测精度，使得建筑物沉降变形预测结果更加准确、可靠。

2.本发明提供的随机森林融合svm预测盾构引起近接结构变形方法将随机森林与支持向量机结合建模，用以预测盾构施工引起地表建筑物变形沉降，规避了现有方法的缺陷，降低模型维度，加快训练速度且解决了预测结果不稳定的问题，为实现建筑物变形预测提供了一种快捷有效的工具。

附图说明

图1为本发明实施例提供的一种预测隧道盾构施工引起建筑物变形沉降方法的流程图。

图2为本发明实施例提供的一种预测隧道盾构施工引起建筑物变形沉降方法中影响因素指标体系图。

图3是本发明实施例中各变量重要性排序图。

图4是本发明实施例中支持向量机模型训练集拟合结果。

图5是本发明实施例中支持向量机模型测试集预测结果。

图6是本发明实施例中建筑物沉降各特征变量敏感性分析图。

图7是本发明实施例提供的预测隧道盾构施工引起建筑物变形沉降装置结构示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

本发明实施例提出一种随机森林融合svm预测盾构引起近接结构变形方法，包括：

根据影响建筑物变形沉降的主要因素，并收集主要因素对应的数据；

建立随机森林模型，根据随机森林模型训练数据并进行重要性度量，得到最优特征集；

将最优特征集降维处理后输入支持向量机模型进行训练，得到建筑物沉降量预测结果。

对所述建筑物沉降量预测结果进行误差分析，获得一阶敏感度和总敏感度，并进行全局敏感性评价，根据下列公式计算一阶敏感度和总敏感度，

其中，si为参数的一阶敏感度值；vi为某个参数的方差；v为系统的总方差；k＝1，2，…ns；sti为参数的总敏感度值；ns为参数样本数，ns＝2m(wmax+1)；wi为参数特征频率；wmax为wi的最大值；m为4或6。

如图1，本发明实施例提出的一种随机森林融合svm预测盾构引起近接结构变形方法，主要包括以下步骤：

步骤1：影响因素指标体系样本数据采集

进一步的，以隧道相关因素包括隧道埋深(x1)和覆跨比(x2)，水文地质因素包括内摩擦角(x3)、弹性模量(x4)和粘聚力(x5)，盾构施工因素包括推进速度(x6)、刀盘扭矩(x7)、推进力(x8)、刀盘转速(x9)、上部土仓压力(x10)和注浆量(x11)，建筑物因素包括相对水平位置(x12)、相对垂直位置(x13)、相对纵向位置(x14)、建筑物自身结构(x15)和建筑物完好程度(x16)等16个指标作为输入变量，以地表建筑物沉降作为输出变量，选取监测的400组数据作为原始训练集，部分数据如表1：

表1样本数据

步骤2：随机森林特征选择

进一步的，将全部数据样本划分为容量为320的训练数据集和容量为80的测试数据集两部分，输入mtry＝5(16/3)，ntree＝800，进行建模。

通过r软件中randomforest程序包来实现importance函数的重要性评价指标可视化绘图，并将不同变量的重要性降序排列，其排列分布如图3所示。由图3可知，隧道埋深，建筑物完好程度、覆跨比、弹性模量等变量重要性度量值比较大，说明这些变量对沉降有较大的影响。

经过5折交叉验证得到不同变量组合的rmse和r²的数值，如表2。

表2不同变量组合时rmse和r²变化表

当变量组合的变量数达到6个时，均方根误差值达到最小，拟合优度也达到了最小，此时模型精度最高。确定最优指标集为隧道埋深，建筑物完好程度、覆跨比、弹性模量，粘聚力，推进力将用于支持向量机模型的构建。

本发明实施例提供的随机森林融合svm预测盾构引起近接结构变形方法充分利用了随机森林能够有效处理有多个指标、噪声干扰的样本数据，筛选有用特征变量的特点，创新性地提高了支持向量机模型的预测精度，使得建筑物沉降变形预测结果更加准确、可靠。

步骤3：支持向量机建模和评估

基于隧道埋深(x1)、建筑物完好程度(x16)、覆跨比(x2)、弹性模量(x4)，粘聚力(x5)，推进力(x8)共6特征，加载r语言中e10711.6-7程序，采用10-cv方法进行参数寻优，输出寻优结果如表3所示。

表3参数优选表

进一步的，参数优化的结果为：bestc＝0.1，bestg＝10，cvmse＝0.6410707(最小)，输入该参数建立支持向量机模型，对训练样本训练拟合结果如图4所示，对测试样本的预测结果如图5所示。从图4可以看出，模拟值和实际值很接近，模拟效果较好。利用训练好的随机森林模型对测试集进行预测，如图5，可以看出随机森林模型测试集上的预测值曲线较贴近真实值。

步骤4：预测结果评价

为了检验随机森林融合支持向量机模型(rf-svm)的优越性，选择未进行特征选择的支持向量机以及人工神经网络进行建模并做对比分析，选用公式(2)均方根误差rmse和公式(3)确定性系数r²和来衡量模型的预测精度，得到误差对比结果如表4。

表4误差比较

结果显示：rf-svm模型、支持向量机模型、人工神经网络模型的均方根误差分别为0.126、0.847、2.05，确定性系数分别为0.997、0.937、0.822。可以看出rf-svm模型预测结果的均方根误差最小且确定性系数最为接近1，说明此模型预测结果最为贴近实际值，精度最高，效果更好。

进一步的，用sobol指数法对监测的400组数据六个指标之间进行一阶和全局敏感性进行分析，得到各指标的敏感度，如图6所示，以建筑物沉降量为目标函数，一阶敏感度和全局总敏感度最高的安全指标是隧道埋深，分别为0.338和0.347，隧道埋深的敏感度明显高于其他参数。建筑物完好程度和推进力的一阶敏感度和全局总敏感度分别为0.238、0.247和0.125、0.135次之，弹性模量，粘聚力，覆跨比的一阶敏感度和全局总敏感度分别为0.128、0.12、0.112和0.134、0.089、0.080，依次减小，

进一步的，说明这3个参数对建筑沉降的影响相对隧道埋深、建筑物完好程度、推进力较小。各参数的一阶敏感性和总敏感性比较接近，说明各指标对建筑沉降的影响规律比较相似。

本发明提供的随机森林融合svm预测盾构引起近接结构变形方法将随机森林与支持向量机结合建模，用以预测盾构施工引起地表建筑物变形沉降，规避了现有方法的缺陷，降低模型维度，加快训练速度且解决了预测结果不稳定的问题，为实现建筑物变形预测提供了一种快捷有效的工具。

如图7所示，本发明实施例提供一种预测隧道盾构施工引起建筑物变形沉降的装置，包括：

获取模块，根据影响建筑物变形沉降的主要因素，并收集所述主要因素对应的数据；

处理模块，建立随机森林模型，根据随机森林模型训练所述数据并进行重要性度量，得到最优特征集；

确定模块，将最优特征集降维处理后输入支持向量机模型进行训练，得到建筑物沉降量预测结果。

本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现本发明实施例中一种随机森林融合svm预测盾构引起近接结构变形方法。

本发明还提供一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现本发明实施例中一种随机森林融合svm预测盾构引起近接结构变形方法。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。

显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。