Excel散点图趋势线回归分析在隧道工程监控量测中的应用

施工论文Excel散点图趋势线回归分析在隧道工程监控量测中的应用沈永林中铁五局集团第四工程有限责任公司第十二隧道作业队贵阳【摘要】论述了Excel2007散点图对数函数趋势线回归分析在隧道工程监控量测中的应用。并结合具体实例，对回归分析步骤进行了详细讲解与说明，总结出了用Excel散点图趋势线进行回归分析的方法。【关键词】隧道工程；监控量测；Excel散点图对数函数趋势线回归分析【引言】目前部分技术人员对隧道监控量测收敛分析心存疑虑，尤其对回归分析与趋势预测这方面的知识讨论最为频繁且争议最大。规范要求监控量测数据分析一般采用散点图和回归分析方法，但并未详细介绍使用何种软件与如何分析。就以前隧道施工经验出发，笔者曾采用过《变形数据分析与预测系统》软件与Excel自带的回归分析功能分析等方法进行回归分析与预测，相比而言Excel自带分析功能进行回归分析与预测最为方便快捷与有效。【正文】隧道监控量测分为以下几个步骤进行:阅读监控量测方案与规范―布设测点并取得初始监测值―现场监控量测及分析―信息反馈与对策。现对几个步骤的操作要求进行详细解释并重点说明第三步中数据分析时采用Excel散点图对数函数趋势线回归分析的方法。一、阅读监控量测方案与规范进行隧道监控量测前必须熟读理解《铁路隧道监控量测技术规程》与监控量测实施方案，并按已批准的方案实施监控量测。二、现场观测得到原始数据按照《TB10121-2007铁路隧道监控量测技术规程》的要求做好量测点的埋设，并按规范要求的频率做好数据采集工作。根据规范，拱顶下沉和净空变化是普通隧道内监控量测必测项目。必测项目埋点间距，埋点位置、量测频率以及量测仪器具体要求如表2-1至表2-4所示。勇于跨越3追求卓越 施工论文表2-1-监控量测断面间距表围岩级别断面间距（m）Ⅴ-Ⅵ5-10Ⅳ10-30Ⅲ30-50表2-2-净空变化量测测线数开挖方法地段一般地段特殊地段全断面法一条水平测线台阶法每台阶一条水平测线每台阶一条水平测线，两条斜测线分部开挖法每分部一条水平测线CD或CRD法上部、双侧壁导坑法左右侧部，每分部一条水平侧线，两条斜测线、其余分部一条水平测线表2-3按距开挖面距离确定的监控量测频率监控量测断面距开挖面距离（m）监控量测频率(0-1)B2次/d(1-2)B1次/d(2-5)B1次/2-3d>5B1次/7d注:B为隧道开挖宽度。表2-4按位移速度确定的监控量测频率位移速度（mm/d）监控量测频率≥52次/d1～51次/d0.5～11次/2-3d0.2～0.51次/3d<0.21次/7d必测项目的监控量测频率应根据测点距开挖面的距离及位移速度分别确定，原则上取2者中频率的较高值。出现异常情况或不良地质时，还应增大监控量测频率。勇于跨越3追求卓越 施工论文表4-监控量测必测项目仪器要求表序号监控量测项目常用量测仪器备注1拱顶下沉水准仪挂钢尺或全站仪2净空变化收敛仪、全站仪监控量测系统的测试精度应满足设计要求。拱顶下沉、净空变化测试精度为0.5-1mm，故当采用m为单位时，记录数据需至小数点后3位数。当采用收敛仪记录净空变化数据时，记录位数由测试仪器的精度而定，一般为mm后2位小数。三、现场监控量测及分析1、数据采集及录入严格根据规范和监测方案要求，采集记录好相关数据，录入Excel后备用。以拱顶下沉为例，录入后的数据具备《表3-1-拱顶下沉测量记录表》所示的形式。勇于跨越3追求卓越 施工论文表3-1-拱顶下沉量测记录表标段名称站前工程施工标施工单位中铁五局单位工程跳厂坡隧道施工方法台阶法施工部位LgDK6+447埋设日期2011-6-24测线编号量测时间　温度修正值修正后测点高程相对初次下沉值（Δμ）相对上次下沉值时间间隔下沉速度备注年-月-日时第一次第二次第三次平均值mmmmmmmmmmmdmm/d12011-6-2420:011247.10501247.10471247.10471247.104801247.105　　　　　12011-6-244:271247.10291247.10301247.10341247.103101247.1031.71.70.53.4　12011-6-258:071247.10191247.10151247.10131247.101601247.1023.21.60.53.1　12011-6-2611:591247.10021247.10031247.09991247.100101247.1004.71.411.4　12011-6-2719:531247.09891247.09851247.09871247.098701247.0996.11.411.4　12011-6-288:491247.09731247.09781247.09751247.097501247.0987.31.211.2　12011-6-293:211247.09641247.09631247.09671247.096501247.0968.31.011.0　12011-6-3020:571247.09551247.09571247.09561247.095601247.0969.20.910.9　12011-7-215:391247.09411247.09481247.09471247.094501247.09510.31.020.5　12011-7-423:521247.09371247.09361247.09371247.093601247.09411.20.920.5　12011-7-714:551247.09271247.09221247.09231247.092401247.09212.41.230.4　12011-7-1012:071247.09141247.09181247.09141247.091501247.09213.30.930.3　12011-7-1320:331247.09091247.09071247.09041247.090601247.09114.20.930.3　12011-7-1610:271247.09011247.08951247.08961247.089701247.09015.10.930.3　12011-7-235:281247.08961247.09001247.08961247.089701247.09015.10.070.0　12011-7-3023:571247.08981247.08991247.08961247.089701247.09015.10.070.0　勇于跨越4追求卓越施工论文勇于跨越4追求卓越 施工论文2、数据分析监控量测数据分析一般采用散点图和回归分析方法。利用Excel中自带的公式功能将原始数据表与分析用数据表链接成一个整体，然后从Excel中插入散点图并链接分析数据进行回归分析与趋势预测。具体步骤如下。（1）原始数据处理监控量测记录的原始数据为一系列高程值或宽度值，需要进行简单的加减运算得出变化速率值以及累计下沉值或累计收敛值。（2）插入散点图在Excel工具条上依次点击[插入]-[图表]-[散点图]-[带平滑曲线的散点图]后即可插入散点图。（3）链接数据鼠标单击散点图，使散点图处于激活状态。点击Excel工具条上的[选择数据]按钮，选择[添加]，依次添加X轴系列值与Y轴系列值。X轴系列值为量测累计时间，Y轴系列为累计下沉值收敛值或变化速率值。数据链接完成点击确定即可生成实测数据线图。生成实测数据线后根据需要点击[图表布局]选项，可在预设的11种图表样式中自定义图表样式。（4）添加趋势线图1-添加趋势线在实测数据线上单击鼠标右键即出现添加趋势线对话框，如《图1-添加趋势线》所示。勇于跨越10追求卓越 施工论文出现设置趋势线格式对话框后，可以选择指数、线性、对数、多项式、幂、移动平均等趋势预测/回归分析类型，隧道施工监控量测一般应选择采用对数函数进行趋势预测与回归分析。中铁隧道勘测设计院有限公司王胜涛硕士在《隧道监控量测的数据回归分析探讨》一文中描述使用对数函数作为回归分析函数，从理论上分析是合理的，工程实例证明也是有效的。函数选择结果如《图2-选择趋势线函数类型》所示。图2-选择趋势线函数类型勇于跨越10追求卓越 施工论文（5）趋势预测当将图2中的预测趋势设置好一个前推值后能预测围岩总变形量以及围岩稳定时间，为隧道开挖预留变形量以及二次衬砌施工时间提供指导性意见。围岩稳定时间趋势预测如《图3-围岩稳定时间预测图》所示。图中选择已实测的6个下沉速率值进行预测，通过函数预测判定开挖第83天后围岩变形速率趋于零，即已基本稳定。围岩总变形量预测方法与围岩稳定时间趋势预测方法相似。变形速率预测曲线不与横轴相交并成上升趋势时意味着隧道变形加剧，需加强支护，及早施工二次衬砌。累计下沉值预测曲线趋于平滑时的最大值为隧道围岩总变形量，若累计下沉值预测曲线不趋于平滑则预示隧道围岩将持续变形。在选择数据步骤时选择非6个实测数据进行分析时预测结果会与上述刚好选择6个实测数据分析时略有差异，具体选择几个数据进行分析由不同的各种情况自定。上述预测是基于统计学基本原理的，即预测具有一般性意义但对某个特定的隧道施工而言最终实际数据不一定恰好是预测的数据。图3-围岩稳定时间预测图勇于跨越10追求卓越 施工论文通过原始数据处理-插入散点图-链接数据-添加趋势线-趋势预测系列五个步骤形成的数据自动分析表如《表3-2-拱顶下沉量测分析表》所式。表3-2-拱顶下沉量测分析表量测编号：T0P-1测点位置：拱顶量测点里程：LgDK6+447埋设时间：2011-6-24　　序号量测日期本次下沉值（mm）累计下沉值（mm）时间间隔（d）量测累计时间（d）下沉速率（mm/d)警报等级施工工序支护参照地质情况12011-6-2400000Ⅲ首次量测按设计施工，无需调整白云岩夹灰岩22011-6-241.61.50.50.53.1Ⅲ开挖32011-6-251.53.10.513.0Ⅲ/42011-6-261.64.7121.6Ⅲ/52011-6-271.66.2131.6Ⅲ/62011-6-281.57.8141.5Ⅲ下部开挖72011-6-291.49.1151.4Ⅲ/82011-6-301.210.4161.2Ⅲ/92011-7-21.311.6280.6Ⅲ/102011-7-41.112.82100.6Ⅲ/112011-7-71.214.03130.4Ⅲ/122011-7-100.914.93160.3Ⅲ/132011-7-130.615.53190.2Ⅲ仰拱开挖142011-7-160.616.13220.2Ⅲ仰拱封闭152011-7-230.016.17290.0Ⅲ/26　量测数据处理结果：1、推算出最终位移量：Umax=16.1mm。2、推算出基本稳定时间：第36天以后，开挖面距该量测断面76m围岩收敛值达16.1mm,收敛率达94％≥90％；位移速度为0.00mm/天≤0.2mm/天，围岩基本稳定。3、评定:该段围岩稳定性适中,初期支护36天后该围岩收敛完成、形成自稳,可进行二次衬砌施工。计算:校核:日期:勇于跨越10追求卓越 施工论文四、信息反馈与对策信息反馈应以位移反馈为主，主要依据时态曲线的形态对围岩稳定性、支护结构的工作状态、对周围环境的影响程度进行判定，验证和优化设计参数，指导施工。铁路隧道监控量测规范对位移管理等级的规定如《表7-位移管理等级表》所式。表7-位移管理等级表管理等级距开挖面1B距开挖面2BⅢUU1B/3U>U2B/3注：u为位移实测值，U1B与U2B为位移控制基准，可按图纸取值或按规范查表得。监控量测结果综合评价设计施工措施，加强监控量测暂停施工继续施工位移（应力）达到Ⅰ级管理工程对策位移（应力）是否超过Ⅱ级管理位移（应力）是否超过Ⅲ级管理否否安全不安全图4-工程安全性评价流程图在实际施工时应根据监控量测结果，按《图4-工程安全性评价流程图》指导施工。根据工程安全性评价的结果，需要变更设计时，应根据有关铁路工程变更管理办法及时进行设计变更。勇于跨越10追求卓越 施工论文工程对策主要包括下列内容1、一般措施（1）稳定开挖工作面措施（2）调整开挖方法（3）调整初期支护强度和刚度并及时支护（4）降低爆破振动影响（5）围岩与支护结构间回填注浆2、辅助措施（1）地层预处理，包括注浆加固、降水、冻结等方法（2）超前支护，包括超前锚杆（管）、管棚、超前插板、水平高压旋喷法、预切槽法等。位移（应力）是未超过Ⅲ级管理时，按正常程序进行施工，并根据预测出的总变形值，对同类围岩开挖时留置优化后的预留变形量,根据预测出的围岩稳定时间指导二次衬砌施工时间。【结束语】合理快速的监控量测分析与预测能优化施工措施，提高企业施工生产率，节约企业施工成本，创造良好效益。运用Excel2007对6个以上的实测数据进行统计分析便基本可对同一性质围岩的总变形量以及围岩稳定时间进行预测。将预测的围岩总变形量值与设计预留变形量值进行对比，可进行优化设计。围岩稳定时间的预测可指导二次衬砌施做时间。采用Excel2007进行回归分析与预测主要过程是：采集数据-通过散点图链接数据-添加趋势线函数-选择对数函数进行趋势预测与回归分析。采用Excel散点图对数函数趋势线回归分析是一种快速有效的监控量测分析方法。【主要参考文献】[1]喻渝、倪光斌等.《铁路隧道监控量测技术规程》[S].北京.中国铁道出版社,2009[2]MicrosoftCorporation.《Excel帮助文件》[K].Redmond.MicrosoftCorporation,2006[3]王秀成、管德鹏等.《铁路隧道工程施工技术指南》[S].北京.中国铁道出版社,2009[4]王毅才、陶纪南等.《隧道工程》[M].北京.人民交通出版社,2005[5]王胜涛、梁小勇等.《隧道监控量测的数据回归分析探讨》[J].隧道建设,2009勇于跨越10追求卓越