一种长期服役抽水蓄能机组性能异常预警方法与流程

1.本发明属于水电机组运行维护技术领域，具体涉及一种长期服役抽水蓄能机组性能异常预警方法。


背景技术：

2.抽水蓄能机组各设备随着累积运行时间的增加，不断劣化。如何有效地确定机组的真实状态，以便及时发现机组异常，合理安排检修，是抽水蓄能机组由计划检修向状态检修转变的重要课题。
3.目前抽水蓄能机组故障诊断技术还不能满足现场需求，其主要原因就是抽水蓄能机组故障的复杂性和故障样本的匮乏。随着电站状态监测系统不断完善，机组运行正常时的状态数据越来越多。如何充分利用现有状态监测数据，挖掘出机组可用的信息至关重要。为及早发现机组异常，需对机组劣化趋势进行提前预警，尽可能减少事故发生的机率。


技术实现要素：

4.针对现有技术中的上述不足，本发明的目的在于提供一种长期服役抽水蓄能机组性能异常预警方法。该方法可以有效地应用于抽水蓄能机组异常状态预警，为机组状态分析提供了一种新的方法。
5.为了达到上述发明目的，本发明采用的技术方案为：
6.提供一种长期服役抽水蓄能机组性能异常预警方法，该方法包括以下步骤：
7.s1、获取并分析抽水蓄能机组运行初期不同工况下的状态监测数据，确定机组的标准健康状态；
8.s2、根据机组的标准健康状态，选取有功功率和工作水头作为反映机组运行状态的敏感特征参数，建立基于分箱最近邻法的回归模型，获得机组标准健康模型；
9.s3、获取机组实时在线监测数据中的运行工况数据，将其输入机组标准健康模型，计算当前工况下的健康值；
10.s4、将健康值和实测值比较，获得偏差，并进行阈值校验，若发生超限，则生成对应报警事件。
11.进一步地，步骤s2中，建立分箱最近邻法的回归模型包括以下步骤：
12.s2-1、根据公式：
[0013][0014]
获取两个对象之间的相似性s
st
；其中：d
st
是对象s和t之间的马氏距离，xs和x
t
是对象s和t的特征向量，s-1
是训练集协方差矩阵的逆矩阵；
[0015]
s2-2、根据公式：
[0016][0017]
获取预测响应ys；其中：y
t
和w
t
分别是第t个邻居的响应和权重，总和超过k个邻居；s
st
是对象s和t之间的相似性，并且总和再次运行在k个最近的邻居上；
[0018]
s2-3、根据公式：
[0019][0020]
获取bin宽度其中：m＝1，2，
…
，10，α为bin阈值；用于定义bin阈值的最佳指数α在[0.1，2]范围内采用圆形搜索算法进行搜索。
[0021]
进一步地，圆形搜索算法具体包括以下步骤：
[0022]
s1、将搜索代理在搜索空间的上限值(ub)和下限值(lb)之间初始化，如下式所示：
[0023]
x
t
＝lb+r
×
(ub-lb)
[0024]
式中：x
t
为搜索个体，另有n个；b和lb分别为变量的上、下限值；r为0到1之间均匀分布的随机数；
[0025]
s2、计算所有搜索个体的适应度函数f(x
t
)：
[0026]
将分别输入到分箱最近邻法模型，根据公式：
[0027][0028]
获取每个位置的代价函数值ji；其中：m为采集数据的个数；yi为实际输出的第i个真实值；为分箱最近邻法模型输出的第i个回归值；
[0029]
s3、令f(xc)＝minf(x
t
)，xc＝x
t
；
[0030]
s4、输入常数值c(c∈[0，1])，令iter＝0，maxiter＝500；
[0031]
s5、如果iter小于maxiter，进入s6；否则，进入s9；
[0032]
s6、更新搜索位置；
[0033]
根据以下公式分别计算w、a、p：
[0034]
w＝w
×
rand-w
[0035][0036][0037]
式中：rand是介于0和1之间的随机数，iter代表迭代计数器，maxiter代表最大迭代次数，随着迭代次数的增加，变量w从-π变为0；变量a从π变为0；变量p从1变为0；
[0038]
根据以下公式计算θ：
[0039][0040]
式中：角度θ从-π到0之间变化；c是介于0和1之间的常数，代表最大迭代次数的百
分比；
[0041]
根据评估的最佳位置xc更新搜索x
t
的位置，如下式所示：
[0042]
x
t
＝xc+(x
c-x
t
)
×
tan(θ)
[0043]
如果更新的搜索位置超出边界，则设置搜索位置等于边界；
[0044]
s7、计算f(x
t
)：
[0045]
如果f(x
t
)《f(xc)，令f(xc)＝f(x
t
)，xc＝x
t
，进入s8；否则，直接进入s8；
[0046]
s8、iter＝iter+1，进入s5；
[0047]
s9、输出f(xc)和xc。
[0048]
本发明的有益效果为：
[0049]
本发明提出了基于分箱最近邻法和圆形搜索法的长期服役抽水蓄能机组性能异常预警方法，引入圆形搜索法的寻优功能对分箱最近邻法的参数进行优化，使模型具有一定的有效性和泛化性，并以某抽蓄机组运行数据为例，对不同负荷、不同水头下机组振动进行了分析，并通过振动实测值与模型预测目标值对比，实现对机组运行状态的评估，为机组运行优化及故障预警提供了参考依据。实际预测结果表明，该方法能获得很好的计算精度，具有很好的应用前景。
附图说明
[0050]
图1为本发明实施例提供的机组运行功率数据示意图；
[0051]
图2为本发明实施例提供的机组运行工作水头数据示意图；
[0052]
图3为本发明实施例提供的机组上机架x向水平振动示意图；
[0053]
图4为本发明实施例提供的抽水蓄能机组振动-功率-水头三维曲面图；
[0054]
图5为本发明实施例提供的测试样本的计算结果示意图；
[0055]
图6为本发明实施例提供的测试样本的绝对误差示意图；
[0056]
图7为本发明实施例提供的测试样本的相对误差示意图；
[0057]
图8为本发明实施例提供的机组健康值与实测值对比示意图；
[0058]
图9为本发明实施例提供的机组偏离度与预警示意图。
具体实施方式
[0059]
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
[0060]
实施例
[0061]
1分箱最近邻法
[0062]
分箱最近邻法根据多数投票标准通过可变数量的k个邻居来预测目标响应。主要思想是考虑与目标具有最大和可比相似性的所有那些邻居进行预测。为了选择最相似的邻居，预先定义了相似区间(即bin)，并根据它们与目标的相似性将邻居分布到这些区间中。所有落入具有最大相似性的bin的邻居都被考虑用于预测。
[0063]
在本发明中，两个对象之间的相似性计算如下：
[0064][0065]
式中：d
st
是对象s和t之间的马氏距离，xs和x
t
是对象s和t的特征向量，s-1
是训练集协方差矩阵的逆矩阵。计算k个邻居的加权平均值，获得预测响应ys，其中权重作为相似性的函数计算，如下：
[0066][0067]
式中：y
t
和w
t
分别是第t个邻居的响应和权重，总和超过k个邻居。s
st
是对象s和t之间的相似性，并且总和再次运行在k个最近的邻居上。
[0068]
分箱最近邻法首先如下定义的相似性度量的分箱过程。
[0069]
相似度箱是通过优化一个参数α，它定义了bin宽度：
[0070][0071]
式中：m＝1，2，
…
，10，α为bin阈值。
[0072]
在本发明中，用于定义bin阈值的最佳指数α在[0.1,2]范围内采用圆形搜索算法方法进行搜索。
[0073]
2圆形搜索算法
[0074]
该算法的具体过程是：
[0075]
(1)初始化。
[0076]
这一步在圆形搜索法中很重要，其中每个搜索位置的整个维度应该是同样随机化的。然后，搜索代理在搜索空间的上限值(ub)和下限值(lb)之间初始化，如下式所示：
[0077]
x
t
＝lb+r
×
(ub-lb)
[0078]
式中：x
t
为搜索个体，另有n个。ub和lb分别为变量的上、下限值；r为0到1之间均匀分布的随机数。
[0079]
(2)计算所有搜索个体的适应度函数f(x
t
)。
[0080]
将分别输入到分箱最近邻法模型，利用下式计算每个位置的代价函数值ji。
[0081][0082]
式中：m为采集数据的个数；yi为实际输出的第i个真实值；为分箱最近邻法模型输出的第i个回归值。
[0083]
(3)令f(xc)＝minf(x
t
)，xc＝x
t
。
[0084]
(4)输入常数值c(c∈[0,1])，令iter＝0，maxiter＝500。
[0085]
(5)如果iter小于maxiter，进入第(6)步。否则，进入第(9)步。
[0086]
(6)更新搜索位置。
[0087]
利用以下公式分别计算w、a、p。
[0088]
w＝w
×
rand-w
[0089][0090][0091]
式中：rand是介于0和1之间的随机数，iter代表迭代计数器，maxiter代表最大迭代次数，随着迭代次数的增加，变量w从-π变为0。变量a从π变为0。变量p从1变为0。
[0092]
利用以下公式计算θ。
[0093][0094]
式中：角度θ在圆形搜索法中起着重要作用，角度θ从-π到0之间变化。c是介于0和1之间的常数，代表最大迭代次数的百分比。
[0095]
根据评估的最佳位置xc更新搜索x
t
的位置，如下式所示：
[0096]
x
t
＝xc+(x
c-x
t
)
×
tan(θ)
[0097]
如果更新的搜索位置超出边界，则设置搜索位置等于边界。
[0098]
(7)计算f(x
t
)。
[0099]
如果f(x
t
)《f(xc)，令f(xc)＝f(x
t
)，xc＝x
t
。进入第(8)步
[0100]
否则，直接进入第(8)步。
[0101]
(8)iter＝iter+1，进入第(5)步。
[0102]
(9)输出f(xc)和xc。
[0103]
以某抽水蓄能电站一台机组在2008年9月22日～2011年12月15日的实测状态监测数据为样本进行研究，验证本技术提供的基于分箱最近邻法和圆形搜索法的长期服役抽水蓄能机组异常状态预警方法的有效性。选取该机组的上机架x向水平振动数据作为研究对象。图1-3给出了2008年9月22日20:49～2008年9月28日16:38机组上机架x向水平振动及对应运行工况的实测数据，从图中可以看出，由于抽水蓄能机组结构复杂、工况转换频繁，使得其振动参数时间序列非常复杂，难以从图中准确地分析其实际运行状态。
[0104]
通过分析可以发现，抽水蓄能机组抽水/发电工况、有功功率、工作水头对机组振动参数有重要影响，由于机组工作水头及运行工况的不断转换，使得机组振动变化及其复杂，不能从振动数据中直接获取机组的真实状态。为实时获取抽水蓄能机组的真实运行状态，需要对其建立健康状态模型。考虑到抽水蓄能机组功率、工作水头对机组振动特性的重要影响，以及分箱最近邻法回归模型对散乱数据良好的拟合性能，本发明基于分箱最近邻法回归模型，建立如图4的抽水蓄能机组振动-功率-工作水头三维曲面模型，获得机组健康状态下功率(p)、水头(h)和振动参数(v)之间的映射关系v＝f(p,h)。引入圆形搜索法的寻优功能对分箱最近邻法的参数进行优化，使模型具有一定的有效性和泛化性。
[0105]
本发明采用机组运行初期，运行状态良好无故障数据，建立机组健康状态下振动标准模型。在2008年9月22日～2009年10月22日1000组数据中，抽取800组能覆盖机组可能的工作水头和有功功率变化区间的健康标准数据建立分箱最近邻法健康模型，将剩下的200组数据作为测试样本进行模型验证。经分析可以看出，基于分箱最近邻法回归模型的机组上机架x向水平振动健康标准模型计算值和实测值几乎完全吻合，计算平均相对误差为0.91％。图5给出了基于分箱最近邻法回归模型的200个测试样本的计算结果。图6-7给出了
测试样本的绝对误差和相对误差，可以清楚地看出，绝大部分样本误差在0附近，健康模型具有很高精度。
[0106]
将机组运行2年后的机组状态监测中的功率、水头等实时在线数据(2011年5月27日～2011年12月12日)代入机组健康模型v(t)＝f(p(t),h(t))，计算当前工况下的状态参数健康标准值v(t)，并和实测值比较，如图8所示。
[0107]
从图中可以看出，机组经过2年运行后，机组振动状态参数开始退化，实测值和健康值之间的偏差越来越大。从图9可以进一步看出，振动偏离度持续波动上升，从2011年11月20日开始，出现显著上升，已超过30％，如果预警线设为30％，即产生相关报警，第一时间定位设备上机架x向水平振动异常，提醒现场工程师关注相关事件。
[0108]
本发明提出了基于分箱最近邻法和圆形搜索法的长期服役抽水蓄能机组性能异常预警方法，引入圆形搜索法的寻优功能对分箱最近邻法的参数进行优化，使模型具有一定的有效性和泛化性，并以某抽蓄机组运行数据为例，对不同负荷、不同水头下机组振动进行了分析，并通过振动实测值与模型预测目标值对比，实现对机组运行状态的评估，为机组运行优化及故障预警提供了参考依据。实际预测结果表明，该方法能获得很好的计算精度，具有很好的应用前景。
[0109]
于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
[0110]
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。