基于DGM(1,1)模型的变形监测数据分析

1、基于GM（1,1）和DGM（1,1）模型的变形监测数据处理与分析郭亚荣（山东科技大学资源与土木工程系，山东 泰安，271000）摘 要：本文以济南路上国奥城项目为例，采用了灰色系统理论中的GM（1,1）和DGM（1,1）模型以及Excel的相关函数功能，对累积沉降量的时间变化曲线进行了预测，证明了模型在变形监测中的有效性和预测的精确性，为提高信息预测的精度和自动化有着积极的意义。关键词：灰色系统，DGM（1,1），GM（1,1），变形监测，数据处理0 引入近些年来，人类对变形监测的重要性逐渐有了深刻地认识，在工程施工与运营期间在建筑物周围布设了监测网并进行了重复的周期性观测，积累了变形观测数据

2、。但由于条件的限制与人为的疏忽，观测资料的保存、分析及利用还不很完善，又因为观测数据自身所隐含的形变信息不能直接获得，所以必须对观测资料做仔细地分析与挖掘，才能更好地对变形做出正确的预测。因此对观测数据做出正确分析和处理、建立合理的预报模型是十分必要的。本文基于GM（1,1）和DGM（1,1）两个模型，以济南鲁商国奥城项目为例进行了相关的数据处理与分析。1 理论介绍1.1 GM（1,1）模型GM（1，1）模型的动态模型为：1）含义为一阶的一个变量的模型；2）a称为发展系数，因为a的大小及符号，反映了(及)的发展态势。3） b为灰作用量，因为b的内涵为系统的作用量，然而b不是可以直接观测的，是通

3、过计算得到的，是等效的作用量，是具有灰的信息覆盖的作用量，故称灰作用量。4）的序列=（，.，）=0.5+0.5称为白化背景值序列。基于每个白化背景值都是与的平均值，故记为MEAN。1.2 DGM（1,1）模型设非负模型：其一次累加生成的序列为： 称为DGM（1,1）模型，或称为GM（1,1）模型的离散形式。是原始序列数据的拟合值，、为待定参数，为迭代基值。1.3 DGM（1,1）模型与GM（1,1）模型的关系1）DGM(1，1)模型全面符合灰色预测模型的建模机理.是一种新的灰色预测模型.或者说是灰色预测模型的一种新形式 ；2）原GM(1，1)模型存在的缺陷在DGM(1，1)中得到了解决,DGM

4、(1，1) 模型可以全面解释原GM(1，1)模型从离散形式到连续形式转变问题 ；3）DGM模型可以看做是GM模型的精确形式，当GM模型中的a取值很小时，二者可替换 。2 应用实例2.1 工程概况鲁商国奥城位于济南市东部，南临经十东路，北临解放东路，东临贤文路，西临规划道路。拟建区内4号楼地上42层，地下4层.按照规范和设计要求需要进行沉降观测。本工程共设水准点3个，K1、K3位于奥体中路西侧，K2位于解放东路南侧,假定K1高程为100m，水准基点采用假定高程。监测网使用的仪器为Trimble DiN03电子水准仪，水准尺为铟钢水准尺。水准点布设图如图2-1所示：图2.1 4#楼水准点布设图2.

5、2 基于GM（1,1）模型和DGM（1,1）模型的数据处理选择2014年4月3日至2014年5月19日的累积沉降量作为原始数据，每五天为一期，共取10期。下面以L8、L9号点为例，其具体方法如下：2.1.1 原始数据两点的原始数据如表2-1所示，表2.1 原始数据观测时间L8点累计沉降值/(mm)L9点累计沉降值/(mm)2014.04.030.185.132014.04.100.415.842014.04.150.495.912014.04.200.555.792014.04.250.725.872014.04.300.835.812014.05.040.965.962014.05.090.

6、965.562014.05.141.055.492014.05.191.155.412.1.2 基于GM（1,1）模型的数据处理使用EXCEL进行数据建模计算，该数据处理过程以L8点为例。 GM(1，1)模型原始数据序列：= （0.18 0.41 0.49 0.55 0.72 0.83 0.96 0.96 1.05 1.15）据式求出相应的1-AGO序列：=（0.1800 0.5900 1.0800 1.6300 2.3500 3.1800 4.1400 5.1000 6.1500 7.3000）根据得出的紧邻均值生成序列 =（0.3850 0.8350 1.3550 1.9900 2.765

7、0 3.6600 4.6200 5.6250 6.7250）若为参数列，且， 则微分方程系数向量a= -0.115533134b= 0.432188173得微分方程动态模型及时间响应函数： (k=0，1，N-1)求出生成数的回代计算值然后计算原始数据的还原值=（0.1800 0.4802 0.5390 0.6050 0.6791 0.7623 0.8556 0.9604 1.0781 1.2101）.计算残差：=（0 -0.0702 -0.049 -0.055 0.0409 0.0677 0.1044 -0.0004 -0.0281 -0.0601）.计算截图：图2.2 GM（1,1）模型计算

8、图1图2.3 GM（1,1）模型计算图2图2.4 GM（1,1）模型计算图32.1.3 基于DGM（1,1）模型的数据处理DGM(1，1)模型原始数据序列： =（0.18 0.41 0.49 0.55 0.72 0.83 0.96 0.96 1.05 1.15）据式求出相应的1-AGO序列：=（0.1800 0.5900 1.0800 1.6300 2.3500 3.1800 4.1400 5.1000 6.1500 7.3000）若为参数数列，且 则DGM（1，1）模型中最小二乘估计为根据上述两矩阵求解参数:= 0.122118491 = 0.460034315计算原始数据的还原值，取，其预

9、测模型为：其还原值为：=（0.1800 0.4820 0.5409 0.6069 0.6810 0.7642 0.8575 0.9623 1.0798 1.2116） 计算残差：=（0 -0.072 -0.0509 -0.0569 0.039 0.0658 0.1025 -0.0023 -0.0298 -0.0616） 计算截图：图2.5 DGM（1,1）模型计算图1图2.6 DGM（1,1）模型计算图2图2.7 DGM（1,1）模型计算图32.1.4 数据的统计分析计算结果统计如下表所示:表2.2 L8沉降数据处理统计表观测时间累计沉降观测值/(mm)GM（1,1）模型模拟值/(mm)DGM

10、（1,1）模型模拟值/(mm)2014.04.030.180.18000.18002014.04.100.410.48020.48202014.04.150.490.53900.54092014.04.200.550.60500.60692014.04.250.720.67910.68102014.04.300.830.76230.76422014.05.040.960.85560.85752014.05.090.960.96040.96232014.05.141.051.07811.07982014.05.191.151.21011.2116残差0.03190.0320平均相对误差7.5017%7.5647%后验差比值0.1988654240.198726197图2.8 L8原始值与模拟值折线图同理，L9原始值与模拟值折线图图2.9 L9原始值与模拟值折线图3 结论本文通过两种模型的对比分析得出，DGM(1，1)预测模型与GM（1，1）预测模型所得出的数值结果大致相同，从实例中证明了DGM（1,1）预测模型可应用于短期沉降观测的数据处理。参考文献1 苗元欣 .基于一元线性回归的变形监测数据处理与分析J.山西建筑，