隧道监控位移量测末段数据回归分析法
【技术领域】
[0001] 本发明属于隧道监控位移量测数据处理方法，特别是指隧道监控位移量测末段数 据回归分析法。
【背景技术】
[0002] 铁路隧道施工主要工序包括开挖、初期支护、二次衬砌等主要工序。根据围岩级别 的不同开挖方法可分为全断面法、台阶法等分部开挖等方法;初期支护主要包括钢支撑、挂 网、喷射混凝土等工序；二次衬砌是指模筑混凝土。《铁路隧道监控量测技术规程（Q/CR 9218-2015)》（以下简称规程）的技术要求明确了上述主要工序中二次衬砌的施作时机，需 要通过隧道监控量测来确定。
[0003] 规程中明确了隧道监控量测的目的，一是可确保施工安全及结构的长期稳定性； 二是验证支护结构效果，确认支护参数和施工方法的合理性，为调整支护参数和施工方法 提供依据;三是确定二次衬砌施作时间；四是积累量测数据，为信息化设计及施工提供依据 等。其中隧道位移相对值已达到总相对位移量的90%以上是二次衬砌施作的必要条件之 〇
[0004] 规程中指出监控量测数据分析可采用散点图法和回归分析法。散点图法是以时间 t为横坐标，累计实测位移值（即规程中所述的隧道位移相对值，下同）u为纵坐标形成散点 对，对散点图中的各点用光滑曲线连接即为实测曲线，然后利用回归分析法对实测曲线的 进一步变化趋势进行预测，预测其位移值(即规程中的中相对位移量，下同），以判断隧道围 岩变化是否稳定，由此可见，趋势预测对下一道工序施工具有重要的作用。
[0005] 规程中还对监控量测回归分析中趋势预测函数类型作了规定，常用的回归函数有 指数函数、对数函数、双曲线函数三种类型。随开挖方法的区别对围岩扰动规律也存在不 同，实测曲线形状也各异，有的是一条完整光滑曲线，有的则呈由几段平滑曲线组成且中间 有突变的台阶状曲线。实际工作中通常有如下两种处理方法:一是整体回归分析法，即用三 种函数分别对该断面所有数据进行回归分析;二是分段回归分析法，即根据实测曲线特性 分段进行，每段都用三种函数进行回归分析。实际应用表明，不加分析的应用整体回归分析 法会出现回归曲线与实测曲线发生偏离的情形，达不到拟合效果，因此预测值与实测值的 偏差量较大而失去对后续施工环节预测指导的意义;而分段回归法因为要对每段曲线都要 用三种函数进行拟合，拟合函数个数会是分段数的三倍，大大增加了工作量，尤其在断面密 集且监测时间长的隧道施工数据分析中，大量数据的处理以及无疑会使工作量成倍增大。
【发明内容】
[0006] 本发明的目的在于提供一种方法简单、数据处理效率高的隧道监控位移量测末段 数据回归分析法。
[0007] 本发明的整体技术构思是：
[0008] 隧道监控位移量测末段数据回归分析法，包括如下工艺步骤：
[0009] A、根据量测位移数据制作散点图，绘制实测曲线；
[0010] B、对实测曲线进行回归分析；
[0011] C、计算预测位移值；
[0012] D、依据《铁路隧道监控量测技术规程(Q/CR 9218-2015)》要求，将累计实测位移值 与预测位移值相比较，判定围岩是否稳定；
[0013] 所述的步骤B是选择实测曲线中靠近最终累计实测位移值不少于5天的末段数据， 对所选择的末段数据进行回归分析。
[0014] 本发明的具体技术构思还有：
[0015] 所述的步骤B是将所选择的末段数据采用对数函数、指数函数、双曲函数进行回归 分析。
[0016] 为保证回归分析的精度并满足其数据量的基本要求，优选的技术方案是，所述的 步骤B是选择实测曲线中靠近最终累计实测位移值不少于7天的末段数据。
[0017] 所述的步骤C是将步骤B中所获得的回归函数的相关系数进行对比，选择相关系数 最大的回归函数，计算预测位移值。
[0018] 所述的步骤D是将累计实测位移值与预测位移值的比值在90%以上作为判定围岩 稳定的依据之一。
[0019] 本发明所取得的实质性特点和显著的技术进步在于：
[0020] 1、本发明提供的方法可针对不同的隧道开挖方法所形成的位移量实测数据进行 回归分析，除可适用于全断面开挖方法形成的整条光滑曲线，尤其适用于台阶法等其他开 挖方法所形成的分段曲线。
[0021] 2、本发明的方法使参与回归分析的数据段或数据量大大减少，可快速的对实测曲 线进行回归，并预测其终值，极大提高了数据处理效率，尤其对隧道中因软弱围岩、密集断 面等情况造成的量测数据量大的情形效果更为明显。
[0022] 3、本发明巧妙地利用了相关性原理对现有分段回归分析法进行了改进，因为越是 末段实测数据对预测值影响越大，经申请人试验证实，该方法不仅具有与分段回归法同等 的预测准确性，而且改进了原整体回归分析法容易出现偏差、分段回归工作量大等技术缺 陷，仅需选择合理可行的末段少量数据便可分析预测围岩的稳定终值，从而大大提高了数 据处理效率。
【附图说明】
[0023] 图1是本发明实施例中全部数据用三种函数整体回归的效果示意图。
[0024]图1中说明实测位移曲线与三条回归曲线拟合程度不高，出现明显偏差。
[0025]图2是本发明实施例中全部数据用三种函数分段回归的效果示意图。
[0026] 图2中说明分段回归中每一段采用三种函数回归，工作量大。
[0027] 图3是本发明实施例中实测数据散点图及实测曲线。
[0028] 图4是对本发明实施例实测曲线中末段数据进行回归的效果示意图。
[0029] 图4中说明仅用末段少量数据即可完成回归，且拟合程度很高。
【具体实施方式】
[0030] 以下结合实施例对本发明做进一步描述，但不应理解为对本发明的限定，本发明 的保护范围以权利要求记载的内容为准，任何依据说明书所作出的等效技术手段替换，均 不脱离本发明的保护范围。
[0031] 实施例1
[0032]以某铁路隧道DK512+116断面为例对比两种现有方法的工作量及曲线回归程度。 图1中的对应函数见表1;图2中对应的函数见表2.a-表2.C。表中t为时间，u为累计实测位移 值，其中t趋于无穷大时的u值为预测值(下同）。
[0038]表2.b分段回归分析第二段结果表
[0033]
[0034]
[0035]
[0036]
[0037]
[0039]
[0040]表2.c分段回归分析第三段结果表 [0041 ] 实施例2
[0042]因台阶法开挖数据特性呈现分段式曲线特征，因此具有代表性，而整体回归可看 成分段回归的特例。故以某铁路隧道DK512+116断面台阶法开挖分段数据回归分析为例，主 要步骤如下：
[0043] A、根据量测位移数据制作散点图，绘制实测曲线。根据表3中的累计位移值原始数 据，以时间t为横坐标，累计实测位移值u为纵坐标，形成散点图并绘制实测曲线见图3。 Γηη/Μ ?
[0046] Β、选择实测曲线中靠近最终累计实测位移值不少于5天的末段数据，对所选择的 末段数据进行回归分析。
[0047]从曲线发展规律看，实测曲线大致分为三段，l-22d为第一段;23-30d为第二段;余 下的为第三段，根据相关性原理，越是末段数据对后期的预测值影响越大，故仅选取第三段 数据(靠近最终累计实测位移值第31d-42d的数据)参与回归分析；回归分析时可不考虑第 一段、第二段对末段数据的影响，直接在第二段数据的基础上，对31d以后的第三段数据用 前述的三种函数进行回归拟合，三种函数的拟合效果均与实测曲线重合度极高，求得回归 函数见表4。
[0048]
[0049] 表4 DK512+116断面净空收敛末段数据回归分析结果表 [0050] C、计算预测位移值；
[0051 ]将表4中的回归函数的相关系数进行对比，选择相关系数最大的回归函数，计算预 测位移值。
[0052]三种函数经对比，双曲函数的相关系数最大，故选择双曲函数作为回归曲线，则其 预测位移值为1/0.0313 = 31.95mm。
[0053] D、将累计实测位移值与预测位移值的比值在90%以上作为判定围岩稳定的依据 之一。
[0054] 从表3中可知，实测最大累计值为31.6mm，则累计实测位移值/预测位移值= 31.6/ 31.95 = 98.9% >90%，按照规程要求还应结合位移速度是否明显下降等因素判定围岩及 支护结构的稳定性，作为是否施作隧道二次衬砌的判定依据。
【主权项】
1. 隧道监控位移量测末段数据回归分析法，包括如下工艺步骤： A、 根据量测位移数据制作散点图，绘制实测曲线； B、 对实测曲线进行回归分析； C、 计算预测位移值； D、 依据《铁路隧道监控量测技术规程(Q/CR 9218-2015)》要求，将累计实测位移值与预 测位移值相比较，判定围岩是否稳定； 其特征在于所述的步骤B是选择实测曲线中靠近最终累计实测位移值不少于5天的末 段数据，对所选择的末段数据进行回归分析。2. 根据权利要求1所述的隧道监控位移量测末段数据回归分析法，其特征在于所述的 步骤B是将所选择的末段数据采用对数函数、指数函数、双曲函数进行回归分析。3. 根据权利要求1所述的隧道监控位移量测末段数据回归分析法，其特征在于所述的 步骤B是选择实测曲线中靠近最终累计实测位移值不少于7天的末段数据。4. 根据权利要求1-3中任一项所述的隧道监控位移量测末段数据回归分析法，其特征 在于步骤C是将步骤B中所获得的回归函数的相关系数进行对比，选择相关系数最大的回归 函数，计算预测位移值。5. 根据权利要求1所述的隧道监控位移量测末段数据回归分析法，其特征在于所述的 步骤D是将累计实测位移值与预测位移值的比值在90%以上作为判定围岩稳定的依据之 〇
【专利摘要】本发明属于隧道监控位移量测数据处理方法，特别是指隧道监控量测末段数据回归分析法。包括根据量测位移数据制作散点图、绘制实测曲线，选择实测曲线中靠近最终累计实测位移值不少于5天的末段数据,对所选择的末段数据进行回归分析，计算预测位移值，依据《铁路隧道监控量测技术规程(Q/CR？9218-2015)》要求，将累计实测位移值与预测位移值相比较，判定围岩是否稳定。本发明克服了现有技术存在原整体回归分析法容易出现偏差、分段回归工作量大等技术缺陷。具有可针对不同的隧道开挖方法所形成的实测数据进行回归分析、数据处理效率高、具有与分段回归法同等的预测准确性等优点。
【IPC分类】G06F17/18
【公开号】CN105608061
【申请号】CN201510979402
【发明人】郝士华, 王艳璐, 娄国充, 牛耀文, 贺雷宁, 周海强, 耿俊杰, 魏晓楠, 李慧勇, 柴波
【申请人】郝士华
【公开日】2016年5月25日
【申请日】2015年12月23日