非等间距多点变形预测模型及其应用

尹 晖，周晓庆，张晓鸣

1.武汉大学测绘学院，湖北 武汉 430079；2.地球空间信息技术协同创新中心，湖北 武汉 430079



非等间距多点变形预测模型及其应用

尹 晖1,2，周晓庆1，张晓鸣1,2

1.武汉大学测绘学院，湖北 武汉 430079；2.地球空间信息技术协同创新中心，湖北 武汉 430079

将彼此关联的多个变形监测点纳入整体建模，将单点的变形分析扩展到空间多点的整体变形分析，采用非等间距等距化处理的改进方法，建立了基于非等间距的多点变形预测模型。本文通过实例分析与对比，验证了非等间距多点变形预测模型的可行性和有效性，是一种非线性时空整体变形分析与预测的新方法。

非等间距序列；多点；整体变形分析；预测模型

采用空间多点的立体布网模式监测工程建筑物的安全状态，用变形观测序列的潜在信息进行定量预测是当前变形监测分析与预报的有效方法。由于受工程现状、观测条件及突发因素等影响，变形观测所获得的数据不一定是大样本的长期观测序列，而是小样本、非等间距的观测数据，因此，针对小样本数据处理的灰色系统理论就具有优越性[1-2]。由于传统GM(1,1)建模要求数据序列是等间距的，不少学者就非等间距GM(1,1)建模方法进行了研究[3-11]，但仍局限在单点或单方向时间序列的研究上，事实上监测工程建筑物的变形状态都布设有大量的变形监测点，监测点间是相互影响、彼此关联的，单点的处理没有利用监测点间相互关系的信息，不足以反映变形体的整体变形趋势和变形规律[12-14]。而事实上，同一变形块体上的变形点，其各向变形规律应具有相似性，通过变形观测序列进行关联分析与描述，在结合工程地质综合考虑的基础上，将彼此关联的多个变形点视为同一块体上的点，从而建立空间多点的整体变形预测模型。

针对变形体的整体变形分析与模型研究，文献[12]首次提出了将单点GM(1,1)模型扩展为多点灰色预测模型的新方法，建立了多点变形预测模型。文献[13]利用多点动态预测模型描述了建筑物整体变形的趋势和规律。近年来，不少学者也进行了相关研究[15-21]，并在系统工程、岩土工程等领域得到了应用，但上述研究均是基于等间距序列所建立的预测模型，对于实际工程可能存在非等间距观测或某几期缺测的实际状况，如何建立非等间距观测的多点变形模型却鲜有报道。本文基于非等间距的变形观测序列，首先提出了一种非等间距等距化处理的改进方法，在此基础上，将单点的GM(1,1)模型扩展为空间多点的变形预测模型，最后再处理为非等间距多点变形预测模型。为验证非等间距多点变形预测模型的可行性和有效性，首先以实例给出了建模的具体步骤及模型拟合和预测结果，并与单点GM(1,1)模型进行了对比分析；再利用文献[12]和文献[13]的实例与等间距多点预测模型进行了对比分析；最后利用文献[14]的实例，采用人为删去部分观测数据，使原等间距序列变为缺测某几期数据的不等间距序列进行了建模与对比，结果表明非等间距多点变形预测模型既适用于观测时间超前或滞后、以及因缺测某几期观测数据导致的非等间距序列的建模，也适用于等间距序列的建模，且具有较高的模型拟合精度与预测精度，是一种非线性时空整体变形分析与预测的新方法。相比单点GM(1,1)建模，利用彼此关联的变形监测点所建立的空间多点整体变形预测模型的拟合及预测精度更高。

1 非等间距多点变形预测模型

1.1 模型建立

与现有非等间距序列等距化处理方法不同[5-7]，考虑到非等间距序列在时间轴上相对等间距而言，可能出现超前或滞后的现象，本文提出一种非等间距序列等距化处理的改进方法，具体步骤如下：

(1) 求平均时间间隔

(1)

(2) 求等距处理后的时间序列

tk′=t1+(k-1)Δt0k=1，2，…,n

(2)

(3) 求等距处理前后各期的时间差

Δtk=tk-tk'=tk-t1-(k-1)Δt0k=1，2，…,n

(3)

(4) 计算各序列值的修正系数

(4)

(5) 计算各序列值的修正值

(5)

式中，Δtk≥0时，μi取μi(tk)；Δtk<0时，μi取μi(tk+1)。

(6) 求得处理后的等间距序列值

(6)

(7)

其一阶累加生成序列为

考虑m个点相互关联，将GM(1,1)模型扩展到空间多点的变形预测模型，其一阶白化微分方程为

(8)

用矩阵表示为

(9)

式中

k=1,2,…,n

(10)

1.2 模型参数求解

(11)

由式(11)可得

(12)

(13)

对式(13)作转置处理，并将k=2,3,…,n代入，得

(14)

按最小二乘准则求得

(15)

根据式(15)求得模型参数A的最小二乘解

(16)

(17)

式(11)整理为

(18)

将k=1,2,…,n代入式(17)，整理得

(19)

其相应的残差方程为

(20)

同理，可根据最小二乘准则求得

(21)

式中

（3）引进大数据处理的相关技术人员。在大数据时代下建设新型电网规划体系需要多方面技术人员的共同合作和努力，大数据不仅仅与计算机科学技术相关，它还涉及到社会经济、电子信息、社会调查等多个学科领域，所以除了电力人员以外，引入大数据相关的专业性人才是非常有必要的，这对新型电网规划体系的建设有着十分重要的意义。

1.3 非等间距多点预测模型

由式(18)并顾及均值化处理式(7)，可得等间距序列的预测模型为

(22)

式(22)中的时序k=1,2,…,n，其时间间隔为1；考虑实际观测序列为非等间距序列tk-tk-1≠const(k=2,3,…,n)，则可以得到非等间距多点预测模型在任一时刻tk的预测式为

(23)

1.4 非等间距多点变形预测模型的检验

为了检验非等间距多点变形灰色模型预测结果的可靠性，需要对其进行检验。一般采用模型的拟合精度和模型预测的平均相对误差来检验[12,19]。

模型的拟合精度

(24)

记ρtk为tk时预测值的平均相对误差，则

(25)

2 实例分析与验证

2.1 非等间距多点预测模型

通过关联性分析，选取位于某基坑北侧边坡上相互关联的4个沉降观测点的监测数据进行整体变形分析与建模，表1给出了A点非等间距的累积沉降观测数据及其等距化处理计算结果。本文取前9期序列用于建模，后5期数据作为模型预测结果的对比分析。下面给出建模的具体步骤及模型拟合预测结果。

2.1.1 原始序列等距化处理

根据上面的计算公式，首先对4个监测点前9期非等间距序列进行等间距化处理，表1给出了A点的等距化处理计算结果。

2.1.2 计算模型参数

由预测式(23)知，预测式只需要模型参数A、C，求解结果为

2.1.3 计算模型拟合值与预测值

经过计算，所建立的多点模型的拟合精度为σ2=0.80，模型的拟合与预测结果如表2所示。

表1 A点等距化处理计算结果

表2 多点模型的拟合与预测结果

2.2 与单点预测模型的比较

为了比较非等间距多点模型与单点模型的模拟和预测结果，对4个监测点分别建立了非等间距单点GM(1,1)模型，各监测点的拟合和预测结果如表3所示。由表2和表3可知，无论是残差均值还是模型的模拟和预测精度，多点整体变形模型的拟合和预测结果均优于单点模型，说明了空间多点模型由于利用了具有相似变形性质、彼此关联和相互影响的多个变形监测点建模，可以获得更好的模型拟合和预测效果。图1给出了非等间距的单点与多点模型的拟合与预测结果的对比。

表3 单点模型的拟合与预测结果

续表

图1 单点与多点模型的对比图Fig.1 Comparison between single point and multi-points modeling

2.3 与等间距多点预测模型比较

为验证非等间距多点变形预测模型的有效性

和适用性，本文对等间距多点变形预测模型也进行了对比分析，数据来源分别为文献[12]和文献[13]的等间距多点变形监测实例。由解算结果可知:①本文计算的模型参数与文献[12]基本一致，但由于文献[12]给出的是图形结果，因此可以量化比较的差异在于本文模型的模拟精度为0.22，文献[12]为0.28，其差异在于预测模型式(23)的细微差别;②本文提出的建模思路因与文献[13]在模型参数求解方式及模型预测公式上完全不同，两者对比结果略有差异，但都具有很好的预测效果。表4给出了取自文献[13]的实测值和残差(实测值减拟合或预测值)以及本文计算的对比结果，前8期为模拟结果，后两期为预测结果，两种模型的相对误差和残差绝对值的均值由本文计算给出。

表4 两种模型的拟合和预测结果对照

2.4 缺测数据的建模与比较

对工程建筑物进行周期性的重复观测是变形监测的基本要求，但由于多种原因可能会出现漏测或缺测的情况，导致变形观测数据为非等间距序列，对此，本文利用文献[14]所给的5个监测点共计24期等间距沉降观测资料，人为假设第9、10期缺测及第12、13、14期缺测的两种情况，如图2、图3实线所示为两组不同的非等间距观测序列，其中断链为数据缺测处。

利用上述非等间距多点变形模型分别对两组数据进行建模，并对后6期进行预测，图2、图3给出了模型的拟合和预测结果。表5给出了取自文献[14]与利用图2、图3缺测数据序列计算得到的预测精度的对比结果，再次验证了非等间距多点变形预测模型的有效性和适用性。

图2 模型1的拟合与预测结果Fig.2 Simulation and prediction results of model 1

图3 模型2的拟合与预测结果Fig.3 Simulation and prediction results of model 2

3 结 论

由以上理论分析与实例计算可以得到以下结论。

(1) 借鉴单变量非等间距等距化的数据处理方法并加以改进，可以将等间距的多点变形预测模型拓展到非等间距的多点变形预测模型，这是一种解决非等间距观测序列多点整体建模与预测问题的有效方法。

表5 不同模型的预测结果对照

(2) 通过实例计算与对比分析，验证了本文所提出的利用彼此关联和相互影响的多个变形监测点，所建立的非等间距多点变形预测模型不仅适用于非等间距序列的建模与预测，也适用于等间距序列的建模与预测，且优于单个点的GM(1,1)建模，是一种通用的多点整体变形预测模型。

(3) 非等间距多点变形预测模型既可针对观测时间与原定观测周期出现超前或滞后而产生的非等间距序列建模，也可直接针对缺测或漏测部分观测周期数据的非等间距序列建模，不需要对缺测数据作内插等数据预处理，具有工程实用性。

(4) 将单点单变量的变形分析拓展到空间多点的整体分析与建模，是时空变形监测分析与预测更为科学和合理的模式和方法，但需要顾及变形监测中多点观测序列的时空关联性，具有时空关联的变形监测点才能纳入整体模型建模，这将是本文后续研究探讨的内容之一。

[1] DENG Julong.Control Problems of Grey Systems[J].Systems & Control Letters,1982,1(5): 288-294.

[2] 尹晖.时空变形分析与预报的理论和方法[M].北京: 测绘出版社,2002.

YIN Hui.Theory and Methodology for Temporal-spatial Deformation Analysis and Forecasting[M].Beijing: Surveying and Mapping Press,2002.

[3] SHEN Chunguang.Grey Models for Non-equidistant Sequence Based on Trapezoid Formula and Central Difference Quotient[J].Journal of Grey System,2010,22(1): 17-24.

[4] WANG Qin,WEI Yong.A Kind of New Strengthening Buffer Operator Suitable for Non-equigap GM(1,1) Model[J].Journal of Grey System,2009,21(1): 105-112.

[5] LI Qiufeng,DANG Yaoguo,WANG Zhengxin,et al.An Extended GM(1,1) Power Model for Non-equidistant Series[J].Journal of Grey System,2012,24(3): 269-274.

[6] 孙虎元,魏绪钧.非等间隔灰色模型及应用[J].应用基础与工程科学学报,1996,4(4): 407-411.

SUN Huyuan,WEI Xujun.Non-equal Time-interval Grey Model and Its Application[J].Journal of Basic Science and Engineering,1996,4(4): 407-411.

[7] 黄声享,李志成.工程建筑沉降预测的非等间距灰色建模[J].地理空间信息,2004,2(1): 41-43.

HUANG Shengxiang,LI Zhicheng.Grey Modeling of Non-equidistant Data Sepuent for Forecasting Subsidence of the Engineering Buildings[J].Geospatial Information,2004,2(1): 41-43.

[8] 李斌,朱健.非等间隔灰色GM(1,1)模型在沉降数据分析中的应用[J].测绘科学,2007,32(4): 52-55.

LI Bin,ZHU Jian.Application of Unequal Interval Grey Model in Analysis of Settlement Data[J].Science of Surveying and Mapping,2007,32(4): 52-55.

[9] 戴文战,李俊峰.非等间距GM(1,1)模型建模研究[J].系统工程理论与实践,2005,26(9): 89-93.

DAI Wenzhan,LI Junfeng.Modeling Research on Non-equidistance GM(1,1) Model[J].Systems Engineering-Theory & Practice,2005,26(9): 89-93.

[10] 史雪荣,王作雷,张正娣.变参数非等间距GM(1,1)模型及应用[J].数学的实践与认识,2006,36(6): 216-220.

SHI Xuerong,WANG Zuolei,ZHANG Zhengdi.GM(1,1) Model with Variational Parameter for Non-equidistant Sequence and Its Application[J].Mathematics in Practice and Theory,2006,36(6): 216-220.

[11] 姜爱平,张启敏.非等间距近似非齐次指数序列的灰色建模方法及其优化[J].系统工程理论与实践,2014,34(12): 3199-3203.JIANG Aiping,ZHANG Qimin.Methods and Optimum of Grey Modeling for Approximation Non-homogenous and Non-equidistant Series[J].Systems Engineering-Theory & Practice,2014,34(12): 3199-3203.

[12] 尹晖,陈永奇,张琰.贫信息条件下的多点变形预测模型及其应用[J].测绘学报,1997,26(4): 365-372.YIN Hui,CHEN Yongqi,ZHANG Yan.Multi-point Deformation Prediction Model with Poor Data Information and Its Application[J].Acta Geodaetica et Cartographica Sinica,1997,26(4): 365-372.

[13] 潘国荣,王穂辉.多点变形动态灰色模型辨识及预测[J].测绘学报,2002,31(S1): 66-69.DOI: 10.3321/j.issn:1001-1595.2002.z1.014.

PAN Guorong,WANG Suihui.Dynamic Grey Modelling Identification and Predication of Multi-point Deformation[J].Acta Geodaetica et Cartographica Sinica,2002,31(S1): 66-69.DOI: 10.3321/j.issn:1001-1595.2002.z1.014.

[14] 尹晖,王尚庆.基于灰关联聚类分析的多点时空非线性模型及其应用[J].武汉测绘科技大学学报,1998,23(2): 100-104.YIN Hui,WANG Shangqing.Multi-point Spatial Nonlinear Model Based on Grey Correlation Clustering and Its Application[J].Journal of Wuhan Technical University of Surveying and Mapping,1998,23(2): 100-104.

[15] 崔立志,刘思峰,吴正朋.基于向量连分式理论的MGM(1,n)模型[J].系统工程,2008,26(10): 47-51.

CUI Lizhi,LIU Sifeng,WU Zhengpeng.MGM(1,n) Based on Vector Continued Fractions Theory[J].Systems Engineering,2008,26(10): 47-51.

[16] 熊萍萍,党耀国,王正新.MGM(1,m)模型背景值的优化[J].控制与决策,2011,26(6): 806-810,815.

XIONG Pingping,DANG Yaoguo,WANG Zhengxin.Optimization of Background Value in MGM(1,m) Model[J].Control and Decision,2011,26(6): 806-810,815.

[17] 刘寒冰,向一鸣,阮有兴.背景值优化的多变量灰色模型在路基沉降预测中的应用[J].岩土力学,2013,34(1): 173-181.LIU Hanbing,XIANG Yiming,RUAN Youxing.A Multivariable Grey Model Based on Background Value Optimization and Its Application to Subgrade Settlement Prediction[J].Rock and Soil Mechanics,2013,34(1): 173-181.

[18] 夏卫国,米传民,刘思峰,等.基于初值改进的多变量MGM(1,m)模型研究[J].中国管理科学,2013,21(S1): 81-85.

XIA Weiguo,MI Chuanmin,LIU Sifeng,et al.The Study on the Improved Multiple Variable MGM(1,m) Model Based on Improving the Initial Value[J].Chinese Journal of Management Science,2013,21(S1): 81-85.

[19] WANG Qijie,WANG Changcheng,XIE Rong’an,et al.An Improved SCGM(1,m) Model for Multi-point Deformation Analysis[J].Geosciences Journal,2014,18(4): 477-484.

[20] GUO Xiaojun,LIU Sifeng,WU Lifeng,et al.A Multi-variable Grey Model with A Self-memory Component and Its Application on Engineering Prediction[J].Engineering Applications of Artificial Intelligence,2015,42: 82-93.

[21] 李晓蕾,刘睿,田永瑞,等.基于灰色预测的空间多点残差修正模型研究[J].大地测量与地球动力学,2010,30(5): 125-128.LI Xiaolei,LIU Rui,TIAN Yongrui,et al.Study on Spatial Multi-point Residual Model Based on Grey Prediction[J].Journal of Geodesy and Geodynamics,2010,30(5): 125-128.

(责任编辑：陈品馨)

Non-equidistant Multi-point Deformation Prediction Model and Its Application

YIN Hui1,2，ZHOU Xiaoqing1，ZHANG Xiaoming1,2

1.School of Geodesy and Geomatics,Wuhan University,Wuhan 430079,China；2.Collaborative Innovation Center for Geospatial Technology,Wuhan 430079,China

Taking consideration of the integrity modeling with correlative multiple deformation points,this paper extends the single point deformation analysis into spatial multi-point integrity analysis and presents a non-equidistant multi-point modeling by improving equispacing processing for the non-equidistant series.The real practical analysis and comparative results indicate that the non-equidistant multi-point prediction model is feasible and effective,which is a new nonlinear approach to the integrated deformation analysis and prediction in time and space domain.

non-equidistant series; multi-point; integrated deformation analysis; prediction model

The National Natural Science Foundation of China (No.51077105)；The Open Research Fund Program of the State Key Laboratory of Geodesy and Earth’s Dynamics (No.SKLGED2013-3-6-E)；The Research Project of State Grid Corporation of China (No.SGSX0000YJJS(2014)457)

YIN Hui(1962—)，female,PhD,professor,majors in deformation analysis and prediction,the theory and methods of spatial information processing.

尹晖，周晓庆，张晓鸣.非等间距多点变形预测模型及其应用[J].测绘学报，2016,45(10)：1140-1147.

10.11947/j.AGCS.2016.20160005.

YIN Hui，ZHOU Xiaoqing，ZHANG Xiaoming.Non-equidistant Multi-point Deformation Prediction Model and Its Application[J].Acta Geodaetica et Cartographica Sinica,2016,45(10):1140-1147.DOI:10.11947/j.AGCS.2016.20160005.

P258

A

1001-1595(2016)10-1140-08

国家自然科学基金(51077105)；大地测量与地球动力学国家重点实验室开放基金(SKLGED2013-3-6-E)；国家电网公司总部科技项目(SGSX0000YJJS(2014)457)

2016-01-04

修回日期：2016-09-22

尹晖(1962—)，女，博士，教授，研究方向为变形分析与预测、空间信息处理理论和方法研究。

E-mail：hyin@sgg.whu.edu.cn