铁路隧道现场监辅助坑道交叉口施工技术分析.doc

铁路隧道现场监控量测回归分析蒲荣宇摘 要：介绍了包西铁路通道大保当至张桥段铁路中新九燕山隧道现场监控量测记录和数据的回归分析方法，对施工中的围岩稳定性、二次衬砌施工时间进行了分析，证明应用铁路隧道现场监控量测结果，能够有利于指导铁路隧道现场施工，科学合理确定开挖方法，及时施工二次衬砌，从而保证铁路隧道施工顺利进行。关键词：铁路隧道 监控量测 围岩稳定 回归分析1 工程概况 包西铁路通道大保当至张桥段铁路新九燕山隧道位于陕西省延安市南川河与劳川河上游分水岭处的劳山川右岸黄土梁峁区，隧道于三十里一沟左侧进洞，下穿既有线西延铁路洪市沟二号隧道，再穿过九燕山分水岭从前黄土沟出洞，地面高程为11581135米。隧道起讫里程为DK514+049DK523+402，全长9353米，为双线隧道，隧道除进口312.77米位于R-2804.53m的曲线上，出口1725.09位于R-4004.53m的曲线上外，其余地段均位于直线上。洞身依次为长5901米的3、长2150米的7的连续上坡，出口为1302米3的下坡。隧道穿越了III、IV、V级围岩 隧道经过区出露主要地层为第四系全新统坡积砂质黄土，中更新统风积砂质黄土，中更新统风积黏质黄土，上第三系红粘土及侏罗纪页岩夹砂岩。隧道通过主要不良地质为滑坡，局部围岩稳定性差，尤其DK521+177DK523+402段洞身位于第三系红粘土层，红色粘土岩为中等膨胀土，含较多疆石结核层富水，受地下水浸泡，对隧道工程影响极大，工程性质较差。地下水类型主要为第四系孔隙水和基岩裂隙水，是造成黄土山体变形的重要条件之一。 新九燕山隧道共设3个辅助坑道，均位于线路右侧，1号斜井与右线线路中线交于DK517+500，平面夹角35，相交处夹角90，最大坡度值11.6%，全长513.06米；2号斜井与右线线路中线交于DK519+720，平面夹角40，最大坡度值11.6%，全长791.7米；3号斜井与右线线路中线交于DK521+400，平面夹角553428，相交处夹角80，最大坡度值11.8%，全长367.15米。 本工点地震动峰值加速度值采用0.05g（相当于地震基本烈度六度），地震动反应谱特征周期为0.35S。2 施工方法 本隧道施工严格按照新奥法施工，充分利用围岩的承载自稳能力，初期支护和二次衬砌对围岩起约束作用，使围岩的变形在有限的空间里，以发挥其承载能力，又阻止围岩过度变形而产生失稳。IV、V级围岩采用环形开挖预留核心土法，III级围岩采用台阶法施工。本隧道施工工期极度紧张，我们在开挖施工过程中，根据围岩监控量测的数据以及回归曲线分析，及时调整开挖方法，变形较小地段采用三台阶法施工，台阶之间的距离为5米，确保仰拱和掌子面的距离40米，二衬和掌子面的距离90米这安全红线。3 监控量测 现场监控量测不仅检测各施工阶段围岩和支护动态，确保施工安全，而且是调整初期支护设计参数，确定二次衬砌和仰拱的施做时间提供合理的科学依据。 洞内外观察、水平相对净空变化值的量测、拱顶下沉量测、红粘土地段的隧底上鼓量测为必须进行的量测项目。量测断面间距V级围岩取5米，IV级围岩取20米，III级围岩取40米，当检测结果有异常时，缩小量测断面间距，加快量测频率。保证一个台阶一条量测线。本隧道在施工过程中，把围岩监控量测作为一项工序纳入到施工中。 新九燕山隧道开挖采用上下台阶法开挖，上台阶开挖采用环形开挖预留核心土法，下台阶开挖左右两侧跳槽开挖。测点布设为：拱顶位置布设一个测点，上台阶布设一对测点，下台阶布设一对测点，共计5个测点，对于地质变化处，位移收敛较明显位置，可以加强测点的布设，不同断面的测点应布置在相同的部位，测点应尽量对称布置，以便数据的相互验证。4 数据处理 现场监控量测数据取得后，及时进行校对、整理、打印相关监控量测报表（以拱顶下沉原始记录表为例），并根据数据绘制位移时间（ut）散点图（图1），以便分析时间效应和空间效应的影响。当ut关系曲线图趋于平稳时，根据散点图的数据分布情况，选择合适的函数进行回归分析，对最大值（最终值）进行预测，并于控制基准值进行比较，综合分析围岩和支护结构的工作状态，对施工工法、工序进行评价，及时反馈评价结论，并提出相应工程对策建议。5 回归分析 从图1可以看出，由于量测的偶然误差所造成的离散型，绘制的ut散点图是上下波动和不规则的，应此，必须进行数据处理才能获得合理的曲线，并以相应的数学公式进行描述。回归分析就是处理测读数据，最终绘制曲线的一种方法。 回归分析是对一系列具有内在规律性的测试数据进行处理，通过处理和计算得到两个变量之间的函数式关系。用这个函数式做出的曲线能代表测试数据的散点分布，用以推算围岩的最终位移和掌握位移变化的规律。项目名称包西铁路新九燕山隧道出口合同段BXS-2标桩号DK523+300埋设日期2008-4-3施工单位中铁十二局包西铁路延安项目部第四分部测量编号量测时间（dt)测点高程修正后测点高程（m)相对初次下沉值（mm)相对上次下沉值（mm)间隔时间下沉速率年月日第一次第二次平均12008431134.16200 1134.16200 1134.16200 1134.16200 22008441134.16056 1134.16057 1134.16057 1134.16057 1.43 1.43 11.4 32008451134.15970 1134.15971 1134.15971 1134.15971 2.30 0.86 10.9 42008461134.15939 1134.15938 1134.15939 1134.15939 2.62 0.32 10.3 52008471134.15920 1134.15921 1134.15921 1134.15921 2.80 0.18 10.2 62008481134.15904 1134.15905 1134.15905 1134.15905 2.96 0.16 10.2 72008491134.15899 1134.15899 1134.15899 1134.15899 3.01 0.05 10.1 820084101134.15897 1134.15898 1134.15898 1134.15898 3.03 0.02 10.0 920084111134.15890 1134.15890 1134.15890 1134.15890 3.10 0.07 10.1 1020084121134.15887 1134.15886 1134.15887 1134.15887 3.14 0.04 10.0 1120084131134.15886 1134.15886 1134.15886 1134.15886 3.14 0.00 10.0 1220084141134.15882 1134.15882 1134.15882 1134.15882 3.18 0.04 10.0 1320084151134.15881 1134.15882 1134.15882 1134.15882 3.19 0.01 10.0 1420084161134.15878 1134.15877 1134.15878 1134.15878 3.23 0.04 10.0 1520084171134.15877 1134.15878 1134.15878 1134.15878 3.23 0.00 10.0 1620084181134.15878 1134.15877 1134.15878 1134.15878 3.23 0.00 10.0 1720084191134.15877 1134.15878 1134.15878 1134.15878 3.23 0.00 10.0 1820084201134.15877 1134.15877 1134.15877 1134.15877 3.23 0.00 10.0 图1 位移时间(ut)散点图5.1 确定变量位移（u）和时间（t）之间的函数类型 确定变量u和t之间的函数类型时，需要根据散点图的数据特点进行分析，从图1可以看出，位移u随时间t的增加而增加，且最初几天增加较快，以后逐步趋于稳定。根据这个特点，采用不同的函数来计算相关系数R，采取相关系数R最接近1的函数来进行回归分析。5.1.1 双曲线函数计算 u=t/(a+bt) 式中，a,b为回归常数，u为位移变形值（mm），t为测点的观测时间（d）。5.1.1.1 将非线性回归方程转化为线性回归方程 从图一中可以看出，双曲线函数的两个变量u和t不是线性关系，需要通过数学处理，替代换元，才能把非线性回归化为线性回归，从而确定回归常数A和B。 将u=t/(a+bt)数学处理、换元如下： 首先将双曲线函数u=t/(a+bt)做以下变形; u=t/(a+bt)两边取倒数，得到： 1/u=(a+bt)/t=b+a*1/t 换元令Y=1/u，X=1/t，b=A，a=B可以得到线性方程为： Y=A+BX5.1.1.2 求回归常数A和B 对图一中的数据进行回归计算，过程如下: X=Xi/n=0.194 Y=Yi/n=0.351 LXX=Xi2-n(X)2=0.912 LYY=Yi2-n(Y)2=0.146 LXY=XiYI- Xi Yi /n=0.341 B=LXY/LXX=0.374 A=Y-bX=0.278 相关系数r=LXY/LXXLYY=0.933 将回归常数A、B代入双曲线方程得： u=t/（0.374+0.278t）5.1.1.3 回归方程线性检验 根据回归方程，判断X和Y是否存在线性关系，对X和Y之间回归直线方程得有效性，用相关系数检验得： r=0.9335.1.2 指数函数计算 u=a*e-b/t 式中，a,b为回归常数，u为位移变形值（mm），t为测点的观测时间（d）5.1.2.1 将非线性回归方程转化为线性回归方程 从图一中可以看出，指数函数的两个变量u和t不是线性关系，需要通过数学处理，替代换元，才能把非线性回归化为线性回归，从而确定回归常数A和B。 将指数函数u=a* e-b/t数学处理、换元如下： 首先将指数函数u=a* e-b/t做以下变形; u=a* e-b/t两边取对数： （u/a）=e-b/t u=a+b（-1/t） 换元令Y= u，X=-1/t,A=a,B=b，可以得到线性方程为： Y=A+BX5.1.2.2 求回归常数A和B 对图一中的数据进行回归计算，过程如下: X=Xi/n=-0.194 Y=Yi/n=1.069 LXX=Xi2-n(X)2=0.912 LYY=Yi2-n(Y)2=0.673 LXY=XiYI- Xi Yi /n=0.781 B=LXY/LXX=0.857 A=exp（Y-BX）=3，441 相关系数r=LXY/LXXLYY=0.998 将回归常数A、B代入指数方程得： u=3.441*e-0.857/t5.1.2.3 回归方程线性检验 根据回归方程，判断X和Y是否存在线性关系，对X和Y之间回归直线方程得有效性，用相关系数检验得： r=0.9985.1.3 对数函数计算 u=a-b/(lg(1+t)式中，a,b为回归常数，u为位移变形值（mm），t为测点的观测时间（d）5.1.3.1 将非线性回归方程转化为线性回归方程 从图一中可以看出，对数函数的两个变量u和t不是线性关系，需要通过数学处理，替代换元，才能把非线性回归化为线性回归，从而确定回归常数A和B。 将对数函数u=a-b/(lg(1+t)数学处理、换元如下： 将u=a-b/(lg(1+t)数学处理如下： 换元令Y=u，X=-/(lg(1+t),A=a,B=b，可以得到线性方程为： Y=A+BX5.1.3.2 求回归常数A和B对图一中的数据进行回归计算，过程如下: X=Xi/n=-1.215 Y=Yi/n=2.96 LXX=Xi2-n(X)2=6.723 LYY=Yi2-n(Y)2=3.548 LXY=XiYI- Xi Yi /n=4.878 B=LXY/LXX=0.726 A= Y-BX=3，841 相关系数r=LXY/LXXLYY=0.999 将回归常数A、B代入指数方程得： u= 3.841-0.726/(lg(1+t)5.1.3.3 回归方程线性检验 根据回归方程，判断X和Y是否存在线性关系，对X和Y之间回归直线方程得有效性，用相关系数检验得： r=0.999 通过以上三种回归方程，指数函数相关系数r得绝对值最靠近1，回归精度较高，所以选用该回归方程来预测最终收敛情况。5.2 求变形稳定时间t（d）及相应得变形量u（mm） 对数函数u= 3.841-0.726/(lg(1+t)最大值 Umax=3.841 当时间为t=40d时，u=3.391，位移速率为0.003mm/d0.2 mm/d 因此，当开挖40d后，可以进行二次衬砌。 图2量测数据处理与回归分析表6 结论 本例仅对拱顶下沉量测数据为依据进行了判定，在实际操作中，要结合监控量测要求得围岩支护状态、周边位移、拱顶下沉、地表沉降等项目，综合考虑，才能得出合理结果，指导施工。新九燕山隧道地质条件复杂，尤其出口1997米为红粘土，遇水膨胀，监控量测工作指导施工尤其重要，本隧道在施工前特制定了监控量测实施细则，专人负责数据的收集、校对、整理、分析、反馈，主要依据时态曲线的形态对围岩的稳定性、支护结构的工作状态、对周围环境的影响程度进行判定，验证和优化隧道施工参数，指导施工。新九燕山隧道在施工中把监控量测工作作为一道工序来做，初期支护未出现掉拱和塌方事故。根据监控量测数据指导二次衬砌的施做时间，很好的保证了隧道施工的安全，成为全线的样板示范工程。一、 建立下沉量(u)随时间(t)发展的时态函数,根据施工规范选取下列函数:指数函数u=ae-(bt);对数函数u=a+b/log(1+t);双曲函数u=t/(a+bt)时间t（天）实测下沉值u（mm)计算回归值(mm)u-里程桩号DK523+380围岩类别 1 指数函数：u=3.442e-0.857/t r=0.998第1天1.431.44-0.01第2天2.302.33-0.032 对数函数：u=3.841-0.723/lg(1+t) r=0.999第3天2.622.64-0.02第4天2.802.81-0.01第5天2.962.910.053 双曲函数u=t/(