基于运行参数相关性的凝汽器真空低跳机征兆捕捉方法

1.本发明涉及大型火力发电机组凝汽器真空低跳机前异常检测技术领域，更确切地说涉及一种基于运行参数相关性的凝汽器真空低跳机征兆捕捉方法。


背景技术：

2.凝汽器系统是火电机组的重要系统之一，它位于汽轮机之后，从汽轮机中排出的蒸汽在凝汽器中遇冷液化成凝结水，从而促成了凝汽器中真空的形成，凝汽器通过凝结水泵将凝结水重新送回锅炉进行循环，循环水泵将循环水注入到表面式凝汽器的铜管中，吸收蒸汽中的热量。一旦其发生异常跳机，不仅会造成安全事故，而且对电厂将造成巨大的经济损失。然而由于其特殊的工作条件，凝汽器系统极易产生异常甚至发生故障，因此，如果能够提前检测出可能发生的跳机故障，就可以做好准备，甚至能够阻止故障的发生。目前，越来越多的研究将目光投向火电机组的异常检测，方法主要分为基于模型和基于数据两大类。基于模型的方法需要我们能够充分理解火电机组的机理，并建立准确的模型，这在复杂的工况下是十分困难的。基于数据的方法有如基于机器学习的方法，通过大量的案例数据，训练出一个合适的模型用于异常的检测。但是大量的用于训练模型的异常案例数据一般是难以获得的。


技术实现要素：

3.本发明要解决的技术问题是，提出一种基于运行参数相关性的凝汽器真空低跳机征兆捕捉方法，根据凝汽器真空及与其相关的运行参数历史数据，在线捕捉到可能发生的跳机异常。
4.本发明的技术解决方案如下：
5.一种基于运行参数相关性的凝汽器真空低跳机征兆捕捉方法，其包括如下步骤：
6.s1、获取待监测汽轮机在未出现跳机异常的正常时刻从凝汽器系统中采集的特征参数的历史采样数据，其中所述特征参数包括由凝汽器甲侧循环水出口压力和凝汽器乙侧循环水出口压力组成的第一特征参数组合、由第一凝汽器甲侧循环水出口温度和第一凝汽器乙侧循环水出口温度组成的第二特征参数组合、由第二凝汽器甲侧循环水出口温度和第二凝汽器乙侧循环水出口温度组成的第三特征参数组合，第一凝汽器甲侧循环水出口温度和第二凝汽器甲侧循环水出口温度由布置在凝汽器甲侧循环水出口的不同温度传感器获取，第一凝汽器乙侧循环水出口温度第二凝汽器乙侧循环水出口温度由布置在凝汽器乙侧循环水出口的不同温度传感器获取；所述历史采样数据所涵盖的时段包含多个时间窗；
7.s2、基于所述历史采样数据，针对每一个特征参数组合，分别计算每个时间窗内两个特征参数的相关性系数并作为线性相关特征，从而得到每一个特征参数组合在正常时刻的线性相关特征范围；
8.s3、对待监测凝汽器系统中的所述特征参数进行实时线上采样，并基于最近一个时间窗内的线上采样数据，计算每一个特征参数组合中两个特征参数的相关性系数并作为
实时线性相关特征；针对第一特征参数组合、第二特征参数组合和第三特征参数组合分别判断其实时线性相关特征是否超出各自在正常时刻的线性相关特征范围，若第一特征参数组合被判断超出，或者第二特征参数组合和第三特征参数组合同时被判断超出，则判定为会有凝汽器真空低跳机异常发生。
9.作为优选，所述相关性系数为皮尔逊相关性系数。
10.作为优选，所述的历史采样数据预先经过数据清洗，剔除凝汽器系统非正常运行状态下的采样数据和其他异常值。
11.作为优选，所述第一特征参数组合在正常时刻的线性相关特征范围为[0.8,1]。
[0012]
作为优选，所述第二特征参数组合在正常时刻的线性相关特征范围为[0.8,1]。
[0013]
作为优选，所述第三特征参数组合在正常时刻的线性相关特征范围为[0.8,1]。
[0014]
作为优选，所述历史采样数据的采样间隔为1秒。
[0015]
作为优选，所述历史采样数据中相邻时间窗的间隔步长为1。
[0016]
相对于现有技术而言，本发明的有益效果如下：
[0017]
由于本发明的基于运行参数相关性的凝汽器真空低跳机征兆捕捉方法是根据运行参数的相关关系是否发生异常变化来进行跳机前异常征兆捕捉的，计算简单，容易实现，泛化能力强。基于相关性系数的方法不需要增加额外的检测设备，也不需要针对汽轮机凝汽器系统建立复杂的物理模型，只需要凝汽器真空和与凝汽器真空相关的机组运行参数的历史数据就能完成跳机前异常征兆捕捉模型的建立，方便推广，且所提出的方法能够较早的检测出可能发生的跳机异常。
附图说明
[0018]
图1为本发明实施例中的凝汽器真空原始曲线图。
[0019]
图2为本发明实施例中的凝汽器甲侧循环水出口压力曲线图。
[0020]
图3为本发明实施例中的凝汽器乙侧循环水出口压力曲线图。
[0021]
图4为本发明实施例中的凝汽器甲侧循环水出口压力与凝汽器乙侧循环水出口压力的相关性系数曲线图。
[0022]
图5为本发明实施例中的凝汽器甲侧循环水出口温度1和凝汽器乙侧循环水出口温度1的相关性系数曲线图。
[0023]
图6为本发明实施例中的凝汽器甲侧循环水出口温度2和凝汽器乙侧循环水出口温度2的相关性系数曲线图。
[0024]
图7为本发明实施例中的归一化后的凝汽器真空与异常指示时间段对照图。
具体实施方式
[0025]
为了更好得理解本技术，将参考附图对本技术的各个方面做出更详细的说明。应理解，这些详细说明只是对本技术的示例性实施方式的描述，而非以任何方式限制本技术的范围。
[0026]
下面本具体实施例用一个火力发电厂凝汽器真空低的真实案例来阐述具体操作步骤以及验证所提出方法的有效性。
[0027]
在该实施例中，提供了一种基于运行参数相关性的凝汽器真空低跳机征兆捕捉方
法，，其具体实现步骤如下：
[0028]
s1、采集凝汽器系统运行参数的历史采样数据
[0029]
本步骤主要是获取待监测机组在未出现跳机异常的正常时刻的多种系统运行参数，然后分析参数间的相关性关系，进而构建历史采样数据集，具体步骤如下：
[0030]
s101、挑选与凝汽器真空相关的机组运行参数变量作为特征参数，特征参数一共6个，包括：凝汽器甲侧循环水出口压力、凝汽器乙侧循环水出口压力、凝汽器甲侧循环水出口温度1、凝汽器乙侧循环水出口温度1、凝汽器甲侧循环水出口温度2、凝汽器乙侧循环水出口温度2。其中，凝汽器甲侧循环水出口温度 1和凝汽器甲侧循环水出口温度2由布置在凝汽器甲侧循环水出口的不同温度传感器获取，凝汽器乙侧循环水出口温度1凝汽器乙侧循环水出口温度2由布置在凝汽器乙侧循环水出口的不同温度传感器获取。
[0031]
s102、不断从凝汽器系统中采样上述6个特征参数的数据，采样的频率为1 秒钟，建立历史采样数据库。
[0032]
s103、识别凝汽器系统在上述采样时段内的运行状态，从历史采样数据中剔除凝汽器系统非正常运行状态下的采样数据和其他异常值。
[0033]
s104、对历史采样数据中各特征参数之间的相关关系进行分析，最终形成三组具有线性相关关系的特征参数组合，具体为：由凝汽器甲侧循环水出口压力和凝汽器乙侧循环水出口压力组成的第一特征参数组合、由凝汽器甲侧循环水出口温度1和凝汽器乙侧循环水出口温度1组成的第二特征参数组合、由凝汽器甲侧循环水出口温度2和凝汽器乙侧循环水出口温度2组成的第三特征参数组合。
[0034]
本实施例中，根据步骤s1，历史采样数据库一共包含了3组特征参数组合的。其中对应的凝汽器真空如图1所示，图中灰色阴影部分表示凝汽器的异常时段。通过图1可以看出，仅通过凝汽器真空曲线并不能及时有效的发现可能的异常，即使后面出现凝汽器真空低的现象，但此时已经太晚了。因此还需选取与凝汽器真空相关的上述三组特征参数组合进行辅助判断。凝汽器甲乙侧循环水出口压力曲线如图2、图3所示。
[0035]
上述第二特征参数组合和第三特征参数组合实际上都是凝汽器甲乙侧的循环水出口温度。设置两组同指标参数的目的在于，可以通过不同的温度传感器来获取凝汽器甲侧循环水出口温度和凝汽器乙侧循环水出口温度消除凝汽器内温度检测过程中不可避免的偶然性，降低因温度误检导致的跳机误报率。不同温度传感器在相应循环水出口侧的具体布置位点不限，保持一定空间间距即可。
[0036]
另外，为了满足后续的统计分析的要求，历史采样数据所涵盖的时段应当足够长，这些时段的历史采样数据按照设定步长用时间窗进行数据提取，每一个时间窗内的数据作为单个时段样本进行统计分析，用于计算一组线性相关特征。由于这些历史采样数据中已经剔除了凝汽器系统非正常时刻运行时刻的数据，因此剩余的数据均为正常运行时刻的数据。
[0037]
s2、计算每一个特征参数组合在正常时刻的线性相关特征范围
[0038]
本步骤中主要对上述三组特征参数组合，基于前述历史采样数据，计算每个时间窗内两个特征参数之间的相关性系数并作为线性相关特征，第一特征参数组合、第二特征参数组合、第三特征参数组合的线性相关特征可分别命名为特征1、特征2、特征3。由此，根据大量的正常时刻历史数据得出正常情况下线性相关系数的取值范围，即得到每一个特征
参数组合在正常时刻的线性相关特征范围。具体包括以下步骤：
[0039]
s201、针对历史采样数据中的每一个时间窗，计算凝汽器甲侧循环水出口压力与凝汽器乙侧循环水出口压力之间的相关性系数，相关性系数采用皮尔逊相关性系数，计算公式为该式中x,y分别为当前时间窗内凝汽器甲侧循环水出口压力与凝汽器乙侧循环水出口压力运行参数时间序列，e表示期望，ρ
1x,y
为x,y之间的相关性系数即特征1，范围为
‑
1至1，负数表示负相关，0表示不相关，相关性系数绝对值越大，相关性越强。
[0040]
s202、针对历史采样数据中的每一个时间窗，计算凝汽器甲侧循环水出口温度1和凝汽器乙侧循环水出口温度1之间的相关性系数，相关性系数采用皮尔逊相关性系数，计算公式为该式中 x,y分别为当前时间窗内凝汽器甲侧循环水出口温度1和凝汽器乙侧循环水出口温度1运行参数时间序列，e表示期望，ρ
2x,y
为x,y之间的相关性系数即特征2，范围为
‑
1至1，负数表示负相关，0表示不相关，相关性系数绝对值越大，相关性越强。
[0041]
s203、针对历史采样数据中的每一个时间窗，计算凝汽器甲侧循环水出口温度2和凝汽器乙侧循环水出口温度2之间的相关性系数，相关性系数采用皮尔逊相关性系数，计算公式为该式中 x,y分别为当前时间窗内凝汽器甲侧循环水出口温度2和凝汽器乙侧循环水出口温度2运行参数时间序列，e表示期望，ρ
3x,y
为x,y之间的相关性系数即特征3，范围为
‑
1至1，负数表示负相关，0表示不相关，相关性系数绝对值越大，相关性越强。
[0042]
s204、历史采样数据中每一个时间窗均获得一个正常时刻的特征1、特征2、特征3后，即可根据这些正常时刻的特征值确定在没有跳机异常发生时的正常时刻的三个线性相关特征范围。
[0043]
在本实施例中，如图4所示凝汽器甲侧循环水出口压力与凝汽器乙侧循环水出口压力的相关性系数曲线，可以看出，在正常时段，凝汽器甲侧循环水出口压力与凝汽器乙侧循环水出口压力之间的相关性系数约为1，在发生异常前，二者的相关性系数突然出现一个向下的波峰，下降到0.8以下后，一段时间后回升。在这期间，相关性系数发生了突变，而在正常的时段内，相关性系数变化平缓且在一个较小的范围内，因此满足s204提出的作为指示变量的要求。然后再根据大量的历史数据，取正常时段凝汽器甲侧循环水出口压力与凝汽器乙侧循环水出口压力的相关性系数取值范围为1至0.8。根据s202、s203，分别计算凝汽器甲侧循环水出口温度1和凝汽器乙侧循环水出口温度1、凝汽器甲侧循环水出口温度2和凝汽器乙侧循环水出口温度2的相关性系数曲线，如图5、图6所示。从图5、图6中可以看出，在正常时段，凝汽器甲侧循环水出口温度1和凝汽器乙侧循环水出口温度1、凝汽器甲侧循环水出口温度2和凝汽器乙侧循环水出口温度2的相关性几乎为1，出现小的波动也在0.9以上，而在异常前，二者的相关性突然减弱，表现在相关性系数上为从1突变成0.7和0.6，均出
现一个向下的波峰，并在较短的时间内出现回升。异常前后，相关性系数有明显的变化，符合 s204中指示变量的要求。然后再根据大量的历史数据，取正常时段凝汽器甲侧循环水出口温度1和凝汽器乙侧循环水出口温度1、凝汽器甲侧循环水出口温度 2和凝汽器乙侧循环水出口温度2相关性系数取值范围为1至0.8。
[0044]
因此，最终确定的三个正常时刻的线性相关特征范围如下：特征1为1至 0.8，即当相关性系数小于0.8的时候判定为凝汽器甲侧循环水出口压力与凝汽器乙侧循环水出口压力相关关系发生了异常，特征2、特征3的阈值范围也为1至 0.8，即当相关性系数小于0.8的时候判定为凝汽器甲侧循环水出口温度1和凝汽器乙侧循环水出口温度1、凝汽器甲侧循环水出口温度2和凝汽器乙侧循环水出口温度2相关关系发生了异常。
[0045]
s3、基于线上采样数据的跳机异常征兆捕捉
[0046]
本步骤主要对对待监测凝汽器系统中的上述6个特征参数进行实时线上采样，并基于最近一个时间窗内的线上采样数据，计算每一个特征参数组合中两个特征参数的相关性系数并作为实时线性相关特征。针对第一特征参数组合、第二特征参数组合和第三特征参数组合分别判断其实时线性相关特征是否超出各自在正常时刻的线性相关特征范围，若第一特征参数组合被判断超出，或者第二特征参数组合和第三特征参数组合同时被判断超出，则判定为会有凝汽器真空低跳机异常发生。具体包括以下步骤：
[0047]
s301、对待监测凝汽器系统中的上述6个特征参数进行实时线上采样，然后基于最近一个时间窗内的线上采样数据(所谓最近一个时间窗是指以当前时刻为基准的最近一段采样时间，该时段内的数据可以代表待监测凝汽器系统最新的运行状态，时间窗大小可以根据实际进行调整)，分别按照s201、s202、s203中的公式计算凝汽器甲侧循环水出口压力与凝汽器乙侧循环水出口压力之间、凝汽器甲侧循环水出口温度1和凝汽器乙侧循环水出口温度1之间、凝汽器甲侧循环水出口温度2和凝汽器乙侧循环水出口温度2之间的皮尔逊相关性系数，依次记为实时特征1、实时特征2、实时特征3，分别判断实时特征1、实时特征2、实时特征3是否小于阈值0.8，如果小于记为1，不小于记为0；
[0048]
s302、根据s301的结果，如果实时特征1的判断结果为1，或者实时特征 2和实时特征3的判断结果同时为1，则最终输出的判断结果为1，表示有异常即会有凝汽器真空低跳机异常发生，否则为0，表示没有异常及没有凝汽器真空低跳机异常发生。
[0049]
本实施例中，根据s302，对之前的三种相关性特征进行综合，即在实时特征1超过阈值范围，或者实时特征2与实时特征3二者同时超过阈值范围(防止温度波动导致误检)的时候才判定为有异常发生。为了便于实现逻辑运算，用1 表示此时为异常状态，0表示为正常状态，检测出异常的曲线如图7所示，为了便于观察，本实施例将归一化后的凝汽器真空曲线也绘制在图中，归一化的公式为其中y
new
为归一化后的值，y
max
为原始一段凝汽器真空最大值，y
min
为原始一段凝汽器真空最小值，y为需要归一化的凝汽器真空。通过图 4可以看出，凝汽器甲侧循环水出口压力与凝汽器乙侧循环水出口压力相关性特征在第1378秒至第1480秒之间超出了阈值，通过图5可以看出，凝汽器甲侧循环水出口温度1和凝汽器乙侧循环水出口温度1相关性特征在第1416秒至第 1476秒之间超出了阈值，通过图6可以看出，凝汽器甲侧循环水出口温度2和凝汽器乙侧循环水出口温度2相关性特征在第1475秒至第1556秒之间超出了阈值。最终结论可以通过图7读出，首次检测到异常的时间
约在第1378秒的时候。而通过电厂的现场记录报警的时间在第1593秒，提前了215秒，相较于最终的跳机时间第1623秒提前了245秒。
[0050]
由此可见，本发明基于皮尔逊相关性系数的方法，利用凝汽器真空和与凝汽器真空相关的参数历史运行数据。在数据预处理阶段，本发明将历史运行参数中的异常数据进行剔除。根据相关性关系，将参数之间的相关性分为线性相关和非线性相关，对于线性相关的参数，使用皮尔逊相关性系数计算公式求得线性相关特征，最终选取的特征为凝汽器甲侧循环水出口压力与凝汽器乙侧循环水出口压力相关性特征、凝汽器甲侧循环水出口温度1和凝汽器乙侧循环水出口温度1 相关性特征、凝汽器甲侧循环水出口温度2和凝汽器乙侧循环水出口温度2相关性特征，记为特征1、特征2、特征3。最后结合三种特征作为最终的异常判断结果，即只有实时特征1超过阈值或者实时特征2、实时特征3同时超过阈值发生异常的时候，才可以判断为有跳机异常发生。
[0051]
需注意的是，上述实施例中两个特征参数之间的相关性系数均采用皮尔逊相关性系数，但事实上也可以采用其他的相关性系数代替。
[0052]
以上所述的实施例只是本发明的一种较佳的方案，然其并非用以限制本发明。有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型。因此凡采取等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。