熵权法并联组合模型在大坝变形监测中的应用

水力发电

熵权法并联组合模型在大坝变形监测中的应用

1，21，21，21，2

郑旭东，陈天伟，邓捷利，段青达，甘

1，2

若，王

雷

1，2

(1.桂林理工大学测绘地理信息学院，广西桂林541004;

2.广西空间信息与测绘重点实验室，广西桂林541004)

摘

“小样本、贫信息”要:针对大坝沉降变形观测存在的相关问题，结合灰色理论适用于的特点，马尔可夫链能够

准确拟合波动性较大数据的特点，ARIMA模型能够挖掘非平稳时间序列的内部信息、建模简单的特点，基于最大熵的原理对模型进行了合理组合，建立了基于最大熵原理的混合预测模型。将此组合模型应用于大坝的实际水平径向位移预测中，分别进行单一模型预测和定权后的组合模型预测，并进行对比，结果显示:熵权法定权后的组合预测模型精度优于ARIMA模型和灰色—加权马尔可夫链模型，且适用性更高，稳定性更强。关键词:ARIMA;灰色—加权马尔可夫链;熵权法;变形预测;大坝

ApplicationofParallelCombinationModelofEntropyEightMethodinDamDeformationMonitoring

ZHENGXudong1，2，CHENTianwei1，2，DENGJieli1，2，DUANQingda1，2，GANRuo1，2，WANGLei1，2

(1．CollegeofGeomaticsandGeoinformation，GuilinUniversityofTechnology，Guilin541004，Guangxi，China;

2．GuangxiKeyLaboratoryofSpatialInformationandGeomatics，Guilin541004，Guangxi，China)

Abstract:Inviewoftheissuesrelatedtothedeformationobservationofdamsettlement，acombinedmodelbasedonentropyweightmethodisestablishedtopredictdamdeformation，inwhich，thegreytheorythatissuitableforsmallsampleandpoorinformation，theMarkovchainthatcanaccuratelyfitthedatawithlargevolatility，andtheARIMAmodelthatcanminetheinternalinformationofnon-stationarytimeseriesandeasytomodelingarerationallycombinedbasedontheprincipleofmaximumentropy．Thecombinedmodelisappliedtothehorizontalradialdisplacementpredictionofadam，andtheresultiscomparedwiththepredictionsofothersingle-factormodels．ThecomparisonshowsthatthecombinedmodelhashigheraccuracythantheARIMAmodelandgreyweightedMarkovchainmodel，andismoreapplicableandstable．KeyWords:ARIMA;gray-weightedmarkovchain;entropyweightmethod;deformationprediction;dam中图分类号:P258;TV698.1

文献标识码:A

文章编号:0559-9342(2019)03-0115-05

0引言

相对误差为0.71%的预测结果［3］

;覃邵峰使用改进的ARIMA模型在大坝变形预测中平均相对误差为0.48%［4］。为了能在实际预测中尽可能地发挥单一模型的预测优势，又能约束各自的缺陷，本文基于信息熵原理提出熵权法并联组合模型。使用灰色—目前，大坝变形预测的方法主要有BP神经网络法、支持向量机法、相关向量机法、时间序列法等［1］

。由于大坝工程的复杂性，以及外界条件的不确定性，决定了使用单一预测模型进行变形预测时，难以得到预测精度和预测期数都较为优秀的预测结果。针对大坝变形预测中存在的种种问题，部分学者进行了如下研究:王利等使用改进灰色模型在大坝沉降预测中取得相对误差小于2%的预测结果［2］

28收稿日期:2018-06-基金项目:广西自然科学基金项目(2017GXNSFAA198308);广140-07-09)西空间信息与测绘重点实验室主任基金项目(15-作者简介:郑旭东(1994—)，男，河南郑州人，硕士研究生，主要研究方向为变形数据处理与应用;陈天伟(通讯作者)．

WaterPowerVol．45No．3

;何启等使用改进的加权马尔可夫链模型得到了平均115

水力发电

2019年3月

加权马尔可夫链预测模型用于解决变形预测中“少数据，贫信息”以及灰色模型预测周期短的问题，同时使用ARIMA模型来提高预测结果的精度，利用最大熵原理对两模型合理组合，得到一个精度更高、更加合理的预测模型。1

模型介绍

1.1

灰色—加权马尔可夫链

1.1.1灰色预测模型GM(1.1)

灰色理论主要研究“小样本，贫信息”的不确定系统的问题［5］

。建立灰色模型时，为了淡化原始数据随机性误差的影响，首先对原始数据进行累加处理，再用微分方程进行建模，最后对模型值进行还原，得到预测值。由于该模型在其他文献中出现较多，本文不再加以详述。1.1.2加权马尔可夫链预测

对于一组随机变量，在验证其满足“马氏性”后，通过各步长的状态转移矩阵，以及其相对应的状态做出预测后，利用其各阶自相关系数能够描述出各种滞时指标值相关关系的强弱这一特性［6］

，对各步长的预测结果按照相依关系的强弱进行加权平均，这就是加权马尔可夫链预测法的基本思想。加权马尔可夫链预测的步骤:(1)计算指标值序列的均值x，均方差s，确定分级标准，根据分级标准确定序列各指标值所属状态。(2)“马氏性”检验。(3)计算各阶自相关系数rk，k∈Ern－kn

k=∑

2

l=1

(xl－x)(xl+k－x)/∑l=1(xl－x)(1)式中，rk为第k阶自相关系数;xi为第i时段的指标值;x为各个指标值的均值;n为指标值序列长度。对各阶自相关系数规范化，即wk=rk/∑m

k=1rk

(2)式中，wk为各个滞时(即步长)的马尔可夫链的权重;m为预测的最大阶数。(4)根据步骤(1)所划分的状态，统计出各步长的状态转移概率矩阵。(5)分别把所需指标值所对应的状态作为初始状态，结合相对应的状态转移矩阵即可得出指标值所对应的状态的概率P(k

)

i，其中，i∈E;k为滞时(步长)，k=1，2，…，m。(6)将同一状态的各预测概率加权求和作为该指标值处于该状态的预测概率，即Pi=∑

mk)k=1

wkP(i，i∈E(3)则最终该指标值的预测状态，即max{Pi，i∈E}。116

WaterPowerVol．45No．3

(7)求出级别特征值对预测值进行估算［7］。待该指标值预测完成后，将预测值加入到原始序列中，重复步骤(1)～(7)即可完成对下一指标值的滚动预测。1.1.3灰色—加权马尔可夫链

灰色预测模型的优点是需要的数据信息较少，但其预测时间较短，而且对于波动性较大的数据，预测效果也不尽人意。而马尔可夫链预测模型是通过统计原始序列得到状态转移概率矩阵进行预测，能够在对变化较大的随机序列的预测中取得不错的效果，但其需要足够多的原始数据来统计出原始序列的内在规律。故将两模型进行串联组合，用马尔可夫链预测模型对灰色模型的结果进行改正，以提高预测精度。其具体步骤如下:(1)将灰色模型的拟合值和实际观测值进行对比，求出相对误差。(2)用求出的相对误差进行马尔可夫链建模，由马尔可夫链预测模型预测出某时段的相对误差。(3)用马尔可夫链预测模型预测出的相对误差对灰色模型的预测值进行改正。1.2ARIMA模型

ARIMA模型也写作为ARIMA(p，d，q)模型，它的实质是将ARIMA(p，q)进行d阶差分得到的。通过这样的操作，把一个非平稳的时间序列转化成一个平稳的时间序列，再通过观察相关函数结尾和拖尾特征得到自回归阶数p和移动平均阶数q，对时间序列进行ARIMA建模，达到预测的效果［8］

。ARIMA模型为φ(B)(1－B)d

Yt=θ(B)et或φ(B)!d

Y(4)t=θ(B)et

式中，φ(B)为自回归算子;d为差分阶数;Yt为时间序列{Yt}在t时刻的值;θ(B)为滑动平均算子;et

为均值为零、方差为δ2

e的噪声序列{et}在t时刻的值;!为差分符号，且(1－B)d

=!d。1.3基于熵权法的组合模型

并联组合模型的关键在于对单一模型的确定及单一模型权重的确定。本文采用熵权法确定各个单一模型权重，建立组合模型来预测大坝变形的数据。熵权法是根据评价对象的指标值构成的判断矩阵来确定权重的一种方法。如果把信息熵理解为某种特定信息xi出现的概率p(xi)，那么信息熵H(x)的定义为H(x)=－∑

mi=1

p(xi)lnp(xi)(5)根据信息熵的定义可知，我们将特定信息认为是单项预测模型的误差，如果单项模型的误差的信第45卷第3期郑旭东，等:熵权法并联组合模型在大坝变形监测中的应用

期的观测值作为原始起算数据，后4期数据和预测数据进行对比。观测原始数据如图2所示。息熵越小，其变异程度也就越大，则该单项预测模型的权重也越小。反之，模型误差的信息越大，那么其变异程度越小，则模型的权重也就越大权重。具体步骤如下:(1)计算单个模型对应时刻的相对误差权重pij=eij

［9］

。因此，可以利用信息熵原理来确定单一模型的∑

nj=1

(6)eij

图2

原始数据

式中，eij为第i个模型第j时刻的相对误差。(2)计算第i个模型相对误差的熵值Hi=－k∑j=1pijlnpij

n

(7)2.1单项模型预测

将灰色模型预测得到的前16期的拟合值与实际式中，k=1/lnn，其中，n为选择的单一模型数量。(3)计算第i个模型相对误差序列的变异程度系数di=1－Hi

(4)计算第i个模型的权重系数ωi=n1(1－di/∑i=1di)n－1观测值对比得到的相对误差，对相对误差进行马氏性检验，用马尔可夫链模型进行建模，预测得到后4期相对误差，用此结果对灰色模型的预测值进行改进，即可得到灰色—加权马尔可夫链预测模型的预测值;将前16期数据进行一阶差分处理，对处理的数据进行自回归系数与偏回归系数分析，确定ARIMA(p，d，q)模型的参数。经分析，确定p=3，d=1，q=3。利用SPSS软件进行ARIMA建模，获得预测值。其结果见表1。(8)(9)(5)将单一模型的预测值和求得的权重加权求和，获得组合模型的预测值Y=∑

ni=1

ωiyi

(10)yi为单一模型的预测值。式中，组合预测模型建模基本流程如图1所示。2.2单项模型熵权计算

根据前文介绍的熵权计算公式以及两个单一模型的变形值的预测结果，可得到灰色—加权马尔可夫链预测模型，ARIMA预测模型在并联组合模型中所占的权重分别为:ω1=0.75，ω2=0.25;再由式(10)求得最终大坝水平径向位移的预测值Y={87.837，87.555，87.196，86.862}。2.3结果分析

为了更加直观的表现出3种模型的预测效果，把3种模型的预测结果与实测值的对比用折线图的形式表现，如图3所示。表中，模型一是指灰色—加权马尔可夫链模型，图1

组合模型建模流程

模型二是指ARIMA模型，模型三是指基于熵权法的组合模型。从表1、图3分析相对误差可得，组合模型预测效果最好，平均相对误差约为－0.26%，其次为ARIMA预测模型，平均相对误差约为0.27%;加权2具体实例

本文以某大坝水平变形监测点水平径向位移观2，10］。这里取前16测值为例，数据来源于文献［表1

期数17181920

实测值/

mm87.9687.8587.4987.07

模型一

预测值/mm87.7487.4287.0386.66

相对误差/%

－0.25－0.49－0.53－0.48

各模型预测值对比

模型二

预测值/mm88.1387.9887.7087.49

相对误差/%

0.190.150.240.48

预测值/mm87.8487.5587.2086.86

模型三

相对误差/%

－0.14－0.34－0.34－0.24

WaterPowerVol．45No．3117

水力发电

表3

期数实测值/mm模型一

预测值/mm相对误差/%

预测值/mm预测结果对比

模型二

相对误差/%

预测值/mm模型三

2019年3月

相对误差/%

17.127.272.0027.097.566.6937.137.353.0

7.03

7.20

2.47

图3

各模型预测值与真实值对比

马尔可夫链预测模型相对较差，平均相对误差为0.46%。单就平均误差而言，组合模型和ARIMA模型差别不大。为了进一步分析3种模型的平稳性以及精确度，按照整体评价预测方法的原则，引入以下评价误差平方和=∑

n(yi－^i=1yi)2

(10)[∑n

^1

22U=1(yi－yi)/n](∑

n^y21

1

(11)22i=1

i)+(∑

n

i=1y2i

)式中，U为Theil不等系数是衡量预测模型预测能力一种指标［11］

，0≤U≤1。U越小说明精度越高，越大说明精度越低。yi为预测值;^yi为真实值;n为数据个数。根据以上两种评价指标，得出预测评价结果见表2。表2

模型的评价指标对比预测模型误差平方和Theil不等系数模型一0.61930.0023模型二0.26630.0015模型三

0.2326

0.0014

实例分析表明，组合模型的预测精度最好，与实际观测数据相比较，Theil不等系数和误差平方和均优于单一的预测模型，显示出它在大坝变形预测精度上的优越性;其次，组合预测模型可以有效的弥补单一预测模型实际预测中存在的不足。对于现有的大坝变形预测模型而言，并联加权模型可以合理的结合各单一模型的优点，提高预测结果精度。为避免其偶然性，笔者使用文献［3］中的数据再次进行实验，求得的模型一和模型二的权重为ω1=118

WaterPowerVol．45No．3

7.221.367.221.437.150.877.201.517.04－1.327.07－0.846.93－1.35

6.96

－0.93

0.14，ω2=0.86。由于篇幅，在此只对结果进行呈现。其结果见表3、4。表4

预测结果精度指标对比

预测模型平均相对误差/%

误差平方和Theil不等系数模型一3.550.31780.0399模型二

1.230.03110.0126模型三

1.18

0.0297

0.0123

由表3、4结果可知，对文献中的大坝水平位移分别进行灰色—加权马尔可夫链模型，ARIMA模型和组合模型的预测，得到的预测结果从平均相对误差，误差平方和，以及Theil不等系数来看组合模型均为最优，由此可见基于熵权法的并联组合模型对于大坝变形预测来说是一种行之有效的预测方法。3结论

本文针对有限的观测数据，结合灰色理论适用于“小样本、贫信息”的特点;马尔可夫链能够准确拟合波动性较大的数据的特点;ARIMA模型能够挖掘非平稳时间序列的内部信息，拟合出变化趋势，从而实现预测目的，具有建模简单的特点;并且结合最大熵的原理对模型进行合理的组合，建立了基于最大熵原理的混合预测模型，通过实例验证了该模型对于小样本序列预测精度较高。具体结论如下:(1)基于信息熵原理，建立并联组合模型来预测大坝的形变量，该方法在精度上较单一模型有所提高，并且平衡了单一模型的实测数据和样本数据之间在质量和数量上的需求，同时能够平衡各单一模型的拟合及预测精度以及普适性等方面的优点及缺点。(2)基于信息熵原理，可对在真实值上下浮动的单一模型的预测值进行合理组合，减少了人为因素，使得组合模型在预测结果上更加接近真实值，从而提高了预测模型的预测精度以及合理性，为数据处理以及大坝变形预测提供了一条新的途经。参考文献:

［1］于鹏，顾冲时．大坝安全监测的组合预测模型［J］．人民黄河，

2006(1):67-68．

第45卷第3期郑旭东，等:熵权法并联组合模型在大坝变形监测中的应用

2013．庆:重庆交通大学，

［8］孙轶轩，邵春福，计寻，等．基于ARIMA与信息粒化SVR组合

模型的交通事故时序预测［J］．清华大学学报:自然科学版，2014，54(3):348-353．

［9］杨锡运，刘欢，张彬，等．基于熵权法的光伏输出功率组合预测

2014，35(5):744-749．模型［J］．太阳能学报，

［10］王利，李亚红，刘万林．卡尔曼滤波在大坝动态变形监测数据

2006(3):353-357．处理中的应用［J］．西安科技大学学报，

［11］王圃，陈荣艳，孙晓楠，等．加权组合模型在城市用水量预测中

2010，18(3):428-434．的应用［J］．应用基础与工程科学学报，

［2］王利，张双成，李亚红．动态灰色预测模型在大坝变形监测及预

2005(3):328-332．报中的应用研究［J］．西安科技大学学报，

［3］何启，戴波．基于灰色神经网络-加权马尔可夫链的大坝变形监

控模型及预报研究［J］．中国农村水利水电，2016(10):146-150，155．

［4］覃劭峰．基于GM(1，1)-ARIMA最优组合的大坝预测模型［J］．

2014(S2):66-69．测绘通报，

［5］邓聚龙．灰色系统理论简介［J］．内蒙古电力，1993(3):51-52．［6］夏乐天．马尔可夫链预测方法及其在水文序列中的应用研究

［D］．南京:河海大学，2005．

［7］朱惠群．马尔可夫链理论在滑坡变形预测中的应用研究［D］．重

(责任编辑焦雪梅)

櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓(上接第53页)作用，地震作用对白鹤滩地下厂房结构不会带来严重危害。［3］刘依松，陈灯红．基于无质量地基模型的重力坝地震响应分析

［J］．三峡大学学报:自然科学版，2014，36(4):6-10．

［4］李静，陈健云．动力相互作用分析中无质量地基的应用研究［J］．

2007，23(2):58-62．世界地震工程，

［5］GB/T13441．1—2007机械振动与冲击人体暴露于全身振动的

评价第1部分:一般要求［S］．

［6］NB350471—2015水电工程水工建筑物抗震设计规范［S］．

参考文献:

［1］樊启祥，汪志林，吴关叶．金沙江白鹤滩水电站工程建设的重大

2018，44(6):5-10．作用［J］．水力发电，

［2］NB/T35011—2016水电站厂房设计规范［S］．

(责任编辑焦雪梅)

櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓櫓(上接第103页)表1

开度/mm2575100150250

充水阀水损

水力损失/m

(出口105°弯管)

1.520.960.850.680.45

度对阀门过流能力影响不大;开度为75～250mm时，105°弯管相比90°弯管能有效提高阀门的过流能力，整个开度范围内的平均过流量增加14.1%，最大过流量增加16.39%。水力损失/m(出口90°弯管)

1.470.960.0.720.58

参考文献:

［1］夏维洪，贾春英．一种充水阀的空穴特性［J］．河海大学学报:自

1980(2):74-87．然科学版，

［2］刘素英．谈广蓄电站上库导流泄洪闸门充水阀设计的体会［J］．

1997(2):12-13．水利水电，

［3］肖段龙，周育桢，董笑波．清远抽水蓄能电站上水库事故检修闸

2014，40(7):79-81．门充水阀设计探讨［J］．水力发电，

［4］龚洪波，万红雁．某电站快速闸门充水方式的优化改造设计［J］．

2015(10):128，221．江西建材，

［5］邓亚新，浦睿诣．新安江水电站水轮发电机组充水方式改进［J］．

2016，37(4):4-491．电力与能源，

［6］陶云冬．抽水蓄能电站进/出水口闸门充水阀的设计探讨［J］．中

2016，16(8):156-157，160．国水运(下半月)，

［7］王福军，唐学林，陈鑫，等．泵站内部流动分析方法研究进展

［J］．水利学报，2018，49(1):47-61．

注:压力用相应的水柱高度表示，进口质量平均总压均为2.04m。

(1)流场随阀门开度变化规律为:0～100mm开度下，弯管在内、外侧壁形成漩涡，减小过流面积，增大能量损耗，开度增大能快速降低漩涡影响范围，增大过流能力;100～250mm开度下，流体顺畅通过阀门，弯管出口流速较均匀，增大开度对流场的改善作用有限。此外，当出口未出现回流现象时，105°弯管出口流场分布更加均匀、流畅;而在小开度下，两种结构紊乱程度均非常高。(2)阀门开度增大能够提高阀门过流能力、降低阻力系数。开度小于100mm时，过流能力提升效果显著，每毫米开度的平均流量增长率可达5.82kg/s;开度大于150mm时候，过流能力提升效果降低，每毫米开度的平均流量增长率为2.0kg/s。(3)开度为0～75mm时，改变阀门出口弯管角(责任编辑高瑜)

WaterPowerVol．45No．3119