一种大坝水平位移预测方法
【专利摘要】本发明公开了一种大坝水平位移预测方法，首先基于典型成因分析，达到对高维因子的降维及去互相关性的效果，并得到对大坝上下游和左右岸两种水平位移进行整体预测的线性变换向量，其次基于聚类分析，实现对样本按性态相似性的划分，及筛选出奇异点，在此划分基础上进行分类回归和预测，可以避免因为样本质量和采样的不均衡，而造成的噪音干扰和信息损失，更为合理、成分的挖掘样本信息，得到更高可信度的预测结果，具有更实际的决策辅助功能。
【专利说明】一种大坝水平位移预测方法
【技术领域】
[0001] 本发明涉及一种预测方法，具体涉及一种大坝水平位移预测方法。
【背景技术】
[0002] 当前我国大坝安全监测总体处于特殊的阶段，很多老坝即将到达设计基准期，而 一大批新坝正投运，因此对大坝安全监测指标进行预测从而提供安全监测决策辅助，就显 得必要了，大坝安全监测和分析评估的主要目的是掌握当前大坝的运行状态，判断大坝是 否安全。目前，大坝安全监测分析的工作模式，主要是是结合日常巡查等方法，分析历史观 测资料，建立数学监控模型，对大坝的工作性态进行分析、评价和监控，并对重要观测指标 量值展开预测，实现大坝安全决策辅助的目的，而大坝水平位移对于大坝安全具有重要的 表征意义。大坝水平位移是一门涉及众多学科和环境因素的应用，难以用确定性的模型给 出精确描述，因此基于统计模型的回归拟合实现大坝水平位移预测，就成为了一种目前广 泛采用的方案。水平位移在环境因素和时效因子的作用下，既呈现出一定的周期性波动特 征，又存在一定的趋势变动，目前行业内主流的分析回归算法，虽然具有原理简单，应用广 泛等优点，但是普遍在高维变量、高噪音、采样不均衡等情况下容易噪音干扰和信息损失， 不一定能保证得到最优拟合，预测可信度不高，甚至极端情况下严重失真。
【发明内容】
[0003] 为了解决上述技术问题，本发明提供了一种大坝水平位移预测方法，解决了现有 方法预测可信度不高的问题。
[0004] 为了达到上述目的，本发明所采用的技术方案是：
[0005] -种大坝水平位移预测方法，包括以下步骤，
[0006] 步骤一，根据大坝水平形变基础成因分析理论，选择影响因子；
[0007] 步骤二，根据影响因子获取其对应的大坝历史水平位移数据，作为历史样本数据， 并归一化历史样本数据；
[0008] 步骤三，对历史样本数据进行典型成因分析，获得效应量与影响因子的线性变换 向量L;
[0009] 步骤四，基于聚类分析算法对历史样本数据进行分类；
[0010] 步骤五，对每种分类下的历史样本数据分布进行最小二乘回归分析，得到对应的 回归模型；
[0011] 步骤六，根据步骤三获得的线性变换向量，对待测的因子向量进行线性变换；
[0012] 步骤七，确定线性变换后的待测的因子向量所属的分类；
[0013] 步骤八，根据分类对应的回归模型，获得预测值。
[0014] 所述影响因子包括多时段下的水位H、气温T、时效t、水位的平方H2、水位的立方 H3、水位的四次方H4和时效的对数Int15
[0015] 所述归一化历史样本数据公式为，
【权利要求】
1. 一种大坝水平位移预测方法，其特征在于：包括以下步骤， 步骤一，根据大坝水平形变基础成因分析理论，选择影响因子； 步骤二，根据影响因子获取其对应的大坝历史水平位移数据，作为历史样本数据，并归 一化历史样本数据； 步骤三，对历史样本数据进行典型成因分析，获得效应量与影响因子的线性变换向量 L ； 步骤四，基于聚类分析算法对历史样本数据进行分类； 步骤五，对每种分类下的历史样本数据分布进行最小二乘回归分析，得到对应的回归 模型； 步骤六，根据步骤三获得的线性变换向量，对待测的因子向量进行线性变换； 步骤七，确定线性变换后的待测的因子向量所属的分类； 步骤八，根据分类对应的回归模型，获得预测值。
2. 根据权利要求1所述的一种大坝水平位移预测方法，其特征在于：所述影响因子包 括多时段下的水位H、气温T、时效t、水位的平方H 2、水位的立方H3、水位的四次方H4和时效 的对数1η\
3. 根据权利要求1所述的一种大坝水平位移预测方法，其特征在于：所述归一化历史 样本数据公式为， F/= (Fi -Vi)/ δ,.式⑴ 其中，Vi'为第i维因子归一化后的值,Vi为第i维因子的原始值，f为第i维因子的 算术平均值，丨为第i维因子的标准差。
4. 根据权利要求1所述的一种大坝水平位移预测方法，其特征在于：获得效应量与影 响因子的线性变换向量的过程为， Al)假设m个影响因子（X1, x2. . . xm)和η个效应量G1, y2. . . yn)，则 u = aTx 式（2) V = bTy 式（3) 其中，a和b分别为影响因子和效应量的线性组合因子，u和V分别为影响因子和效应 量变换后的单维投影向量； A2)求解u和V的皮尔逊相关系数Corr (u, V)，其公式和约束条件为；
3- Σ (Xj χ) 3, - I bT2 (y, y)b = I 为了使的u和v有最大关联，则C〇rr(u，v)最大，根据式（4)和约束条件可得，在式 (4)中分母不变的情况下，分子最大时Corr (u, V)最大，分子为u和V的协方差cov (u, V)= aTcov (x, y) b, 根据拉格朗日乘数法原理得公式如下，
其中，η为拉格朗日乘数法获得的函数，D(x)为X自身的协方差矩阵，D(y)为y自身的 协方差矩阵，λ和Θ表示约束条件等式算子的系数； A3)在式（5)两端分别对a和b求偏导，得如下方程组： Σ (X，y) b_ λ Σ (X，X) a = O 式（9) Σ (y, x)a- θ Σ (y, y)b = O 解方程组得λ = Θ = aTcov(x，y)b，即λ和Θ相等，都等于u和v的协方差； Α4)对式（9)进行移项得： Σ (X，X) Σ (X，y) b = λ a 式（10) Συ Σ (y，X)a =入 b 将式（10)写出矩阵与向量之积的形式，得到式（6);
A5)根据式（6)求解矩阵的本征值和本征向量，其最大的本征值对应的本征向量，即为 效应量与影响因子的第一主成分的线性组合系数，次大本征值对应的本征向量为第二主成 分的组.合系数向量，依次类推可得所有的变量主成分向量，根据设定的本征损失阈值，确 定所需的组合系数向量数目，获得了线性变换矩阵L。
5.根据权利要求1所述的一种大坝水平位移预测方法，其特征在于：基于聚类分析算 法对历史样本数据进行分类的过程为， BI)假设每个类簇的中心由若干局部密度小于P C1的点围绕，并且这些点距离其他有高 局部密度的点的距离大于d。，其中P0由人为定义，d。为截断距离，d。的值是使得平均每个 点的邻居数为样本点总数的1%?2% ; B2)计算样本点I局部密度P I
其中，X(dI;r-d。）为一种算子，当dI;r-d。为负时，值为1，否则为0, dIT为样本点I的样本 点J的距离，样本点J的局部密度高于样本点I，即样本点J是样本点I高局部密度的点； B3)计算样本点I到样本点J的距离δ Σ ;
其中Pj为样本点J的局部密度； Β4)根据计算获得的P Σ和δ i对样本点进行分类； 分类规则为：样本点的P :小于P cl，但是δ :大于d。，则认为该样本点为奇异点；样本 点的P I大于P C1,并且δ Σ大于d。，则认为该样本点为类簇中心；为每个类簇中心定义一类 另IJ，排除奇异点和类簇中心的其余样本点的类别为领域内距离样本点最近并且局部密度高 于该样本点的样本点所属的类别。
6. 根据权利要求1所述的一种大坝水平位移预测方法，其特征在于：步骤六中，对待测 的因子向量进行线性变换的公式为， Y' = VL 其中，V'为线性变换后的待测的因子向量，V为线性变换前的待测的因子向量。
7. 根据权利要求1所述的一种大坝水平位移预测方法，其特征在于：确定线性变换后 的待测的因子向量所属的分类的过程为，根据各类簇的中心，分别计算待测因子向量与各 类簇的中心的距离，距离最短的类簇即为该待测因子向量所属的分类。
【文档编号】G06Q10/04GK104376370SQ201410588297
【公开日】2015年2月25日 申请日期:2014年10月28日 优先权日:2014年10月28日
【发明者】花胜强, 蔡杰, 夏智娟, 陈意, 周锡琅, 高磊, 向南, 郑健兵 申请人:南京南瑞集团公司, 国网电力科学研究院