《变形监测数据的分析与预报研究.doc》由会员分享,可在线阅读,更多相关《变形监测数据的分析与预报研究.doc(10页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、【精品文档】如有侵权,请联系网站删除,仅供学习与交流变形监测数据的分析与预报研究.精品文档.变形监测数据的分析与预报研究王伟彪(南京工业大学 土木工程学院 江苏 南京 210009)摘要:变形数据分析与预报是变形监测数据处理的重要内容。利用有限的监测数据,选取合理的预测模型准确预报未来的变形是变形监测数据处理研究的热点问题之一。在此主要是研究灰色系统理论GM(l, 1)模型在预测模型中的应用,利用长期积累的观测资料掌握变化规律,对建筑物的未来性态做出及时有效的分析和预报。关键词:变形监测; 系灰色系统理论;GM(l,1)模型1、.引言系统的沉降观测资料是验证地基基础是否正确和分析地基事故及判别
2、施工质量的重要资料。沉降观测值是诸多因素在某一空间和特定时间上综合作用结果的反映。沉降观测应具备监测和预报双重效能。运用灰色系统理论处理观测数据11,建立灰色加权模型是作好预报的有效途径。在灰色模型预测中,采用灰色加权模型22,以减少量化的盲目性和随意性,提高预报精度。2、 GM(1,1)预测模型2-1、 数据的生成以及模型的建立33灰色理论的微分方程模型称为GM模型,而GM(1,1)模型表示一阶的一个变量的微分方程模型。设对于时间序列,由已测到的变形数据作为原始数据序列,记为:。对以上数据序列作累加生成运算可得到一组新的数据序列一一累加生成序列,即 (1)对于单变量序列若采用GM(1,1)模
3、型,其方程可写为: (2)式中为待辨识参数(称为发展系数,它反映了及发展的态势;称为灰作用量)。将待辨识参数列记为,则参数列可用最小二乘法求解,即 -(3) 式中 -(4) -(5)对于时间序列 有初始条件是的初值,且。若与中各元素是一一对应的。则(2)式的解为,或 -(6)因此,可在确定之后一次求得累加生成数的回代生成进而可通过累减生成运算得到原始数据序列的还原值: -(7)预2-2、 3 预测精度的检验 预测模型建立后,其预测效果能否满足实际要求,需要对其进行检验,再介绍这里使用较多的后差检验法。因计算值与实测值之间的残差: -(8)其相应的均值: - (9)残差之方差: -(10)又原始
4、数据序列均值: -(11)原始数据序列方差: -(12)于是,后验差检验指标为: 、后验差比值: 、小误差概率:根据值划分的精度等级4*,如表1示。数据经检验达到“合格”或以上指标,方可按式(7)式随后期进行预测。 表1 精度等级划分表指标优合格勉强合格不合格0.950.800.700.350.500.65可得:指标C越小越好,C越小,表明尽管原始数据很离散而模型所得计算值与实际值之差离散程度小。指标P越大越好,P越大,表明残差与残差平均值小于给定值的点较多,预报精度高。 3、4 GM(1,1)加权预测模型5 数据生成* 建立等时间隔的离散函数:表示观测序数,经累加生成AGO:;其中:模型的建
5、立 由原始数据列和生成数列组成矩阵 - (13) -(14)对于在某一时间内可靠性随时间成正比变化的数据序列可按下式定权。式中为精度递增因子,取值范围。 -(15)求灰参数向量: -(16)建立GM(1,1)加权模型: -(17)4、 5 精度检验应用灰色加权模型预报建筑物沉降是否可信,必须按一定的途径进行检验。在残差预测模式中,检验数不是一次算出来的,而是根据前面的数据推算出后一个数据,依次递推地检验。每一检验值对模型来说都是后验值,所以称为后验差检验。因计算值与实测值之间的残差: -(18)其相应的均值: -(19)残差之方差: -(20)又原始数据序列均值: -(21)原始数据序列方差:
6、 -(22)于是,后验差检验指标为: 、后验差比值 、小误差概率 根据值划分的精度等级,如表1所示。数据经检验达到“合格”或以上指标,方可按(17)式随后期进行预测。 5、6 算例分析为了便于应用灰色预测模型,根据上述灰色GM模型,利用基于MATLAB灰色系统沉降预测程序6*。其主要功能是可根据输入的原始数据列,通过处理得到模型计算值和预测值,并对模型进行精度等级评定。为了完成某住宅楼的安全建设,在施工过程中对01栋建筑物进行监测。图1为沉降观测点的布置情况。 图 1 沉降观测点的布置图 1 沉降观测点的布置为了验证预测模型分析的可行性以及估计预报的效果,我们取第1沉降点前15期数据进行计算及
7、预报比较。 表 2 为 沉降点1的 15期沉降观测值期号12345678本期高程(m)9.79129.79109.79099.79079.79069.79049.79039.7903沉降量(mm)00.20.10.20.10.20.10.0累积沉降量(mm)0.00.20.30.50.60.80.90.9期号9101112131415本期高程(m)9.79019.78979.78949.78909.78869.78819.7879沉降量(mm)0.20.40.30.40.40.20.2累计沉降量(mm)1.11.51.82.22.63.13.3利用GM(1,1)模型对沉降点1进行分析期号本期高
8、程(m)沉降量(mm)累积沉降量(mm) dfdfffffgddgfhgfdf()90(0909090dsfsdffd ()()19.79120029.79100.20.239.79090.10.349.79070.20.559.79060.10.669.79040.20.879.79030.10.989.79030.00.999.79010.21.1109.78970.41.5119.78940.31.8129.78900.42.2139.78860.42.6149.78810.23.1159.78790.23.3利用GM(1,1)模型对沉降点1进行分析现以1点的沉降量为例说明模型的运用方法
9、和预测精度的检验过程n=15,则原始数据序列为:=0.0 0.2 0.3 0.5 0.6 0.8 0.9 0.9 1.1 1.5 1.8 2.2 2.6 3.1 3.3,对原始数据序列作一次累加生成,得到数据序列= 0.0 0.20 0.50 1.00 1.60 2.40 3.30 4.20 5.30 6.80 8.60 10.80 13.40 16.50 19.80,则由(3)-(6)式分别计算得到:B=-0.1,1;-0.35,1;-0.75,1;-1.3,1;-2.0,1;-2.85,1;-3.75,1;-4.75,1;-6.05,1;-7.7,1;-9.7,1;-12.1,1;-14.
10、9,1;-18.15,1X=0.2;0.3;0.5;0.6;0.8;0.9;0.9;1.1;1.5;1.8;2.2;2.6;3.1;3.3= 0 0.3736 0.4463 0.5331 0.6368 0.7607 0.9087 1.0855 1.2966 1.5489 1.8502 2.2101 2.6401 3.1536 3.7671,=0 -0.1736 -0.1463 -0.0331 -0.0368 0.0393 -0.0087 -0.1855 -0.1966 -0.0489 -0.0502 -0.0101 -0.0401 -0.0536 -0.4671,由(9)-(12)式分别算得:
11、=0.1264 =1.0578进一步算得后验差比值:0.1195又均满足小误差概率:,所以因为,所以以1点的沉降量的预测精度为优,同样可以得到其他观测点的沉降量的预测模型及预测精度。则GM(1,1)非加权模型为:利用GM(1,1)加权模型对沉降点1进行分析,利用GM(1,1)加权模型取R=0.9,则按式(15)有权序列为:,P=1.5 2.25 3.375 5.0625 7.5938 11.3906 17.0859 25.6289 38.4434 57.6650 86.4976 129.7463 194.6195 291.9293 ;求灰参数向量 =则GM(1,1)加权模型为:两种预测模型拟合
12、残差对比分析见表3,加权模型非加权模型曲线见图2。表3 两种预测模型拟合残差对比分析期号本期高程(m)沉降量(mm)累积沉降量(mm) GM(1,1)非加权预测值(mm)拟合残差(mm)GM(1,1)加权预测值(mm)拟合残差(mm)19.7912000.0000.00000.00000.000029.79100.20.20.3736-0.17360.64540.1408439.79090.10.30.4463-0.14630.7544-0.10899949.79070.20.50.5331-0.03310.88190.009459.79060.10.60.6368-0.03681.03090
13、.011369.79040.20.80.7607-0.03931.20500.093779.79030.10.90.9087-0.00871.40860.052589.79030.00.91.0855-0.18551.6465-0.116899.79010.21.11.2966-0.19661.9247-0.1200109.78970.41.51.5489-0.04892.24990.0362119.78940.31.81.8502-0.05022.63000.0436129.78900.42.22.2101-0.01013.07430.0926139.78860.42.62.6401-0.0
14、4013.59360.0715149.78810.23.13.1536-0.05364.20070.0662159.78790.23.33.7671-0.46714.9104-0.3400图2 第1号点加权模型与非加权模型预测沉降量 图3 第2号点加权模型与非加权模型预测沉降量同理可以根据观测点2的沉降观测值如表4,得到加权模型非加权模型曲线见图3。表 4 2号沉降点2的 15期沉降观测值期号12345678本期高程(m)9.68839.68819.68799.68779.68769.68749.68749.6873沉降量(mm)0.30.20.20.20.10.20.00.1累积沉降量(mm
15、)0.30.50.70.91.01.21.21.3期号9101112131415本期高程(m)9.68689.68669.68639.68599.68569.68549.6850沉降量(mm)0.50.20.30.40.10.20.4累计沉降量(mm)1.82.02.32.72.83.03.4根据图2,图3 加权模型与非加权模型预测沉降量的比较可知:GM(1,1)加权预测模型预报结果精度高,置信度高。期号123456789101112131415本期高程(m)9.68839.68819.68799.68779.68769.68749.68749.68739.68689.68669.68639.
16、68599.68569.68549.6850沉降量(mm)0.30.20.20.20.10.20.00.10.50.20.30.40.10.20.4累计沉降量(mm)0.30.50.70.91.01.21.21.31.82.02.32.72.83.03.4根据图2,图3 加权模型与非加权模型预测沉降量的比较可知: GM(1,1)加权预测模型预报结果精度高,置信度高。67、 小结:灰色预测方法主要用于小样本和要求不高的情形,其优点是样本需求量少,计算简单,精度可靠,高层建筑工程开挖及施工中具有很多动态变化,因而监测数据在不断的充实,预测模型也可以随之不断更新、优化。数据的可信程度及数据间的相对精
17、度以权重表示。数据的可靠性越高,则赋予的权重越大,表示数据建模中的可信度越大。在灰色模型预测中,采用灰色加权模型,以减少量化的盲目性和随意性,提高预报精度。参考文献 邓聚龙灰色系统基本方法M武汉:华中理工大学出版社,1987. 李希灿动态平差灰色预测优化模型J测绘工程,1999,8(1):3435 陶本藻测量数据统计分析M北京:测绘出版社1992 吴翊,等应用数理统计M长沙:国防科技大学出版社,1995 中国科学院数学研究所数理统计组回归统计方法M北京:科学出版社,1974 苏金明,阮沈勇MATLAB6.1实用指南M北京:电子工业出版社,2002A study of the analysis
18、and forecast of deformation monitoring dataWangWeiBiaoCivil Engineering College NanJing of Technology 210009 JiangSu NanJingAbs tract: The analysis and forecast of data are major contents in the processing of deformation monitoring data. How to use finite deformation monitoring data and choose ratio
19、nal model to accurately forecasts deformation is one of the most important-task of the study of the processing of deformation monitoring data. This thesis research the application on the prediction model of Grey System. Theory GM (1.1) model of Grey System Theory is used in prediction . to make out the analysis and predict about the future behavior of the constructional work in time by holding the variety rule using long-term accumulated surveying datum.Key words:Deformation monitoring; Gray system theory; GM(1,1) model
限制150内