监控量测数据的分析方法.docVIP

PAGE PAGE 7 监控量测数据的分析方法 1、监控量测的意义和主要内容 新奥法施工是应用岩体力学的原理，以维护和利用围岩的自稳能力为基点，将锚杆和喷射混凝土集合在一起作为主要支护手段，及时进行支护以控制围岩的变形及松弛，使围岩成为支护体系的组成部分，形成以锚杆、喷射混凝土和隧道围岩为一体的三位一体的支护结构。通过现场监控量测，及时反馈围岩及支护体系的力学动态，以判断支护结构的合理性，指导隧道的设计与施工。 监控量测的所得的数据有一定的离散性，它包含着偶然误差的影响，不经过数据处理是难以直接利用的，监控量测对围岩变形的最终值和围岩稳定的时间的预测一般通过回归分析的方法得到。 监控量测的数据处理包括以下主要内容：根据量测值绘制时态曲线，即观测数据的散点图。确定回归函数，预测最终变形值与围岩最终稳定的时间。 下面通过实例说明监控量测的数据处理方法和如何使用计算机得到处理结果。提到的计算机软件有EXCEL和MTALAB。 2、回归分析模型的求解 回归分析模型一般采用对数模型、指数模型、双曲函数模型等，目前比较容易的求解方法是把回归模型转化为一元函数，对一元函数进行回归分析。 2.1、一元函数回归分析的处理方法 若变量初始值不为零时，即。用一元函数做回归分析 ，采用最大似然对参数进行参数估计。 参数： 参数：　　　　　　　　　 　　　　　　　　 回归精度为： 若初始值变量为零时，即，则可以选作为回归函数，采用最大似然法进行对参数进行参数估计。 参数： 2.2、可以转化为一元函数的回归模型 2.2.1指数模型 令：，， 转化成的一元函数模型为： 2.2.2 令：，，， 转化成的一元函数模型为： 2.2.3 令：，，，， 转化成的一元函数模型为： 3、围岩稳定所需时间的求解 围岩稳定的判断为隧道二次衬砌施工的时间提供准确的信息，隧道二次衬砌施做时间必须是围岩变形稳定后，当围岩变形量过大或初期支护的变形不收敛，又难以及时补强时，可以提前施做二次衬砌，此时二次衬砌必须予以加强。同类围岩支护参数必须予以加强。所以要判断围岩的变形速率和最终稳定时间。 围岩的变形速率为： 引入临界变形速率作为围岩基本稳定的判断依据 对数模型的求解： 引入临界变量后： 因此有： 4、举例说明回归模型的使用和用软件的解决方法 4.1 初始数据分析 某隧道的拱顶沉降监控量测数据如表1： 表1 观察时间 累计沉降值 观测时间 累计沉降值 2007.11.08 1.5 2007.11.30 21.2 2007.11.09 2.9 2007.12.2 22.9 2007.11.10 4.9 2007.12.4 23.5 2007.11.11 5.2 2007.12.6 25.5 2007.11.12 6.1 2007.12.8 25.3 2007.11.13 8 2007.12.10 25.9 2007.11.14 8.2 2007.12.13 26.5 2007.11.15 9.7 2007.12.16 27.3 2007.11.16 10.3 2007.12.19 27.4 2007.11.17 12.5 2007.12.26 27.9 2007.11.18 12.9 2008.1.2 28.3 2007.11.19 13.2 2008.1.9 28.4 2007.11.21 14.1 2007.11.23 16 2007.11.25 16.9 2007.11.27 18 2007.11.28 19.2 2007.11.29 20.8 以观测时间时为间序号绘制该组数据的散点图，如图1： 图1 沉降观测散点图 4.2 线性插值法的使用目的与方法 由于沉降观测频率不同，观测时间间隔不相等，无法直接将沉降观测的数据在Excel中进行分析，用线性插值的方法求出间隔1天对应的累计沉降量。 例如某天的沉降累计值为26.3，7天以后测得的沉降值为0.6，则这7中每天的累计沉降值经过线形插值计算后如表2： 表2 时间 平均沉降值 每天沉降值 　 　 26.3 第1天 0.085714286 26第2天 0.085714286 26第3天 0.085714286 26第4天 0.085714286 26第5天 0.085714286 26第6天 0.085714286 26第7天 0.085714286 26.9 为了求解方便，线性插值的过程在MATLAB程序中实现。 MTALAB中线性插值法程序的编码如下： k=0; for i=1：n for j