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

基于GM （1,1）和DGM （1,1）模型的变形监测数据处理与分析

郭亚荣

（山东科技大学资源与土木工程系，山东 泰安，271000）

摘 要：本文以济南路上国奥城项目为例，采用了灰色系统理论中的GM （1,1）和DGM （1,1）模型以及Excel 的相关函数功能，对累积沉降量的时间变化曲线进行了预测，证明了模型在变形监测中的有效性和预测的精确性，为提高信息预测的精度和自动化有着积极的意义。 关键词：灰色系统，DGM （1,1），GM （1,1），变形监测，数据处理

0 引入

近些年来，人类对变形监测的重要性逐渐有了深刻地认识，在工程施工与运营期间在建筑物周围布设了监测网并进行了重复的周期性观测，积累了变形观测数据。但由于条件的限制与人为的疏忽，观测资料的保存、分析及利用还不很完善，又因为观测数据自身所隐含的形变信息不能直接获得，所以必须对观测资料做仔细地分析与挖掘，才能更好地对变形做出正确的预测。因此对观测数据做出正确分析和处理、建立合理的预报模型是十分必要的。本文基于GM （1,1）和DGM （1,1）两个模型，以济南鲁商国奥城项目为例进行了相关的数据处理与分析。

1 理论介绍

1.1 GM （1,1）模型

GM （1，1）模型的动态模型为：

b k az k x =+)()()1()0(

1）含义为一阶的一个变量的模型；

2）a 称为发展系数，因为a 的大小及符号，反映了)

0(x

(及)

1(x )的发展态势。

3） b 为灰作用量，因为b 的内涵为系统的作用量，然而b 不是可以直接观测的，是通过计算得到的，是等效的作用量，是具有灰的信息覆盖的作用量，故称灰作用量。 4）)()

1(k z

的序列

)()1(k z =（)2()1(z ，)3()1(z ，...，)()1(n z ）

)()1(k z =0.5)()1(k x +0.5)1()1(-k x

称为白化背景值序列。基于每个白化背景值)()1(k z 都是)()1(k x 与)1()

1(-k x 的平均值，故记)1(z 为

MEAN )1(x 。

1.2 DGM （1,1）模型

设非负模型：

{}

)(,),2(),1()()0()0()0()0(n x x x k X ⋯=

其一次累加生成的序列为：

{}

)(,),2(),1()()1()1()1()1(n x x x k X ⋯=

称(1)

(1)

12(1)()X

k x k ββ∧∧+=+为DGM （1,1）模型，或称为GM （1,1）模型的离散形式。

(1)

()x k ∧是原始序列数据的拟合值，1β、2β为待定参数，(1)

(1)x ∧为迭代基值。

1.3 DGM （1,1）模型与GM （1,1）模型的关系

1）DGM(1，1)模型全面符合灰色预测模型的建模机理.是一种新的灰色预测模型.或者说是灰色预测模型的一种新形式 ；

2）原GM(1，1)模型存在的缺陷在DGM(1，1)中得到了解决,DGM(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.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.03 0.18 5.13

2014.04.10 0.41 5.84

2014.04.15 0.49 5.91

2014.04.20 0.55 5.79

2014.04.25 0.72 5.87

2014.04.30 0.83 5.81

2014.05.04 0.96 5.96

2014.05.09 0.96 5.56

2014.05.14 1.05 5.49

2014.05.19 1.15 5.41

2.1.2 基于GM （1,1）模型的数据处理

使用EXCEL 进行数据建模计算，该数据处理过程以L8点为例。

GM(1，1)模型原始数据序列：

()

)()0()0(i X X == （0.18 0.41 0.49 0.55 0.72 0.83 0.96 0.96 1.05 1.15）

据式∑==n

i i X i X

1

)0()

1()()( 求出相应的1-AGO 序列：

{}

)()1()1(i X X ==（0.1800 0.5900 1.0800 1.6300 2.3500 3.1800 4.1400 5.1000 6.1500

7.3000）

根据)1(5.0)(5.0)()1()1()

1(-+=i X i X i z

得出)1(X 的紧邻均值生成序列)1(Z =（0.3850

0.8350 1.3550 1.9900 2.7650 3.6600 4.6200 5.6250 6.7250）

若T

b a a ),(=∧

为参数列，且

⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡⋯⋯=)()3()2()0()0()0(N X X X Y ，⎢⎢⎢⎢⎢⎣⎡-⋯⋯--=)

()3()2(z )1()1()1(n z z B ⎥⎥⎥⎥⎦⎤

⋯111 则微分方程系数向量

Y

B B B b

T T 1)(-∧

=⎥⎦

⎤

⎢⎣⎡=αα

a= -0.115533134 b= 0.432188173

得微分方程动态模型及时间响应函数：

b X dt

dx

=+)1(α αααb

e b X k X k +⎥⎦⎤⎢⎣

⎡-=+-∧∧

)1()1()0()

1( (k=0，1，…，N-1)

求出生成数的回代计算值)()

1(i X

然后计算原始数据的还原值

)1()()()

1()

1()

0(--=∧∧∧i X

i X

i X

=（0.1800 0.4802 0.5390 0.6050 0.6791 0.7623 0.8556 0.9604

1.0781 1.2101）.

计算残差：