CN115496278A - 一种水电机组定子线圈温度预警方法及系统 - Google Patents

一种水电机组定子线圈温度预警方法及系统 Download PDF

Info

Publication number
CN115496278A
CN115496278A CN202211143698.3A CN202211143698A CN115496278A CN 115496278 A CN115496278 A CN 115496278A CN 202211143698 A CN202211143698 A CN 202211143698A CN 115496278 A CN115496278 A CN 115496278A
Authority
CN
China
Prior art keywords
temperature
historical
value
loss function
error
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211143698.3A
Other languages
English (en)
Inventor
李超顺
吴一凡
陈鹏
余甜
肖湘曲
黄洁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN202211143698.3A priority Critical patent/CN115496278A/zh
Publication of CN115496278A publication Critical patent/CN115496278A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/20—Hydro energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Data Mining & Analysis (AREA)
  • Marketing (AREA)
  • Mathematical Analysis (AREA)
  • Operations Research (AREA)
  • Health & Medical Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Tourism & Hospitality (AREA)
  • Computational Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Quality & Reliability (AREA)
  • Development Economics (AREA)
  • General Engineering & Computer Science (AREA)
  • Game Theory and Decision Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)
  • Probability & Statistics with Applications (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明提供了一种水电机组定子线圈温度预警方法及系统，属于水电机组状态评估与预警领域，方法包括：获取整体温度趋势数据和周期温度趋势数据；基于运行状态损失函数的GRU温度预测模型，获取机组当前运行阶段未来时刻的预测温度值；计算历史误差的标准差、均值和置信度；当历史误差的标准差小于预设标准差，结合机组当前运行阶段未来时刻的预测温度值，判断是否进行温度预警；当历史误差的标准差大于预设标准差，则结合置信度和当前运行阶段未来时刻的预测温度值，判断是否进行温度预警。本发明提高了模型训练的效果，同时也降低了环境温度对温度误报警的影响。

Description

一种水电机组定子线圈温度预警方法及系统
技术领域
本发明属于水电机组状态评估与预警领域，更具体地，涉及一种水电机组定子线圈温度预警方法及系统。
背景技术
水电机组相关部件的运行温度是水电站运行监测中重要的非电气量监控参数，运行温度与机组的运行状态和部件的工作寿命密切相关，能够直观反映机组是否正常稳定运转。研究针对水电机组运行温度的预测及预警方法，对于实现设备预测性维护，减少不必要损失，具有十分重要的实际意义。然而，受环境温度变化影响，机组运行温度在一年中会呈现出截然不同的特征，这些特征差异会影响到预测模型的效果。此外，机组温度与机组运行状态直接相关，随着机组频繁的开停机过程，机组温度也会呈现出相应地周期性变化，只有使模型深度挖掘出温度随机运行状态的变化，才能提高模型的预测精度。
发明内容
针对现有技术的缺陷，本发明的目的在于提供一种水电机组定子线圈温度预警方法及系统，旨在解决现有的水电机组定子线圈温度预警方法中存在预测精度不足以及受环境温度影响，导致误报警次数多的问题。
为实现上述目的，本发明提供了一种水电机组定子线圈温度预警方法，包括以下步骤：
S1：根据水电机组定子绕组温度变化的特征，将原始温度序列周期化处理，并使用STL分解方法将周期化后的原始温度序列进行分解，获取整体温度序列和周期温度序列；
S2：对整体温度序列滤波处理，将滤波后的残差以及STL分解方法分解后的残差加入至周期温度序列中，获取周期温度趋势数据；
S3：将机组当前运行阶段的周期温度趋势数据中若干历史温度实际值输入至基于运行状态损失函数的GRU温度预测模型中，输出机组当前运行阶段的未来若干时刻的预测温度值；
S4：将机组历史运行阶段的温度变化划分为温度上升阶段和温度稳定阶段，确定在温度上升阶段当前运行阶段温度预测点在当前运行阶段和历史运行阶段中的位置，计算历史温度预测值与历史温度实际值的误差，以及计算历史误差的标准差、均值和置信度；
S5：当历史误差的标准差小于预设标准差，则对系统误差进行校正，结合机组当前运行阶段未来时刻的预测温度值，判断是否进行温度预警；当历史误差的标准差大于预设标准差，则结合置信度和当前运行阶段未来时刻的预测温度值，判断是否进行温度预警。
进一步优选地，将原始温度序列周期化处理的方法为：
根据机组运行次数和每次运行的实际时长计算机组平均运行时间；
若机组的运行时间多于机组平均运行时间，则将机组运行时间多于平均运行时间的时间段对应原始温度序列除去；若机组的运行时间少于机组平均运行时间，则将机组运行时间少于平均运行时间的时间段进行温度补齐，补齐的温度数值为机组运行阶段温度达到稳定后的均值温度。
进一步优选地，GRU温度预测模型的损失函数包括均方误差损失函数和运行状态损失函数；运行状态损失函数包括平均温度损失函数、温度上升占比损失函数、最高温度损失函数和中位温度损失函数；
平均温度损失函数用于计算单次机组运行过程中实际温度平均值与预测温度平均值之间误差；温度上升占比损失函数用于计算机组当前运行阶段中温度上升所用时间占所有时间的比值；最高温度损失函数用于计算当前运行阶段中实际温度最高值和预测温度最高值的误差；中位数温度损失函数用于计算当前运行阶段中实际温度的中位数值和预测温度中位数值的误差。
进一步优选地，温度上升阶段的结束条件为连续的相邻五个温度预测点平均温度之间的变化不超过0.1℃。
进一步优选地，计算历史误差的标准差、均值和置信度的方法，包括以下步骤：
计算若干历史运行阶段和当前运行阶段中每个温度预测点在周期温度趋势数据中的位置；
在历史运行阶段中寻找与当前运行阶段中温度预测点相同位置的温度预测点，根据历史运行阶段中温度预测点处的历史温度预测值与历史温度实际值，计算预测误差以及置信度；
计算历史运行阶段中温度预测点的误差均值和标准差。
在历史运行阶段中任意一个温度预测点所处位置为：
Figure BDA0003854727690000031
其中，历史运行阶段中温度上升阶段为t～t+Ti r；t′表示历史运行阶段中温度预测点的时间；
在当前运行阶段中温度预测点所处位置为：
Figure BDA0003854727690000032
其中，t″表示当前运行阶段中温度预测点时间；Tr表示当前运行阶段中的温度上升时间；
Figure BDA0003854727690000033
分别为五个历史运行阶段中的温度上升时间。
进一步优选地，温度预警规则为：
Figure BDA0003854727690000034
其中，
Figure BDA0003854727690000035
为历史运行阶段中温度预测点的标准差；C为温度预测点的置信度；
Figure BDA0003854727690000036
表示当前运行阶段中未来时刻的预测温度值；TH表示阈值；
Figure BDA0003854727690000037
为温度预测点的误差均值。
另一方面，本发明提供了一种水电机组定子线圈温度预警系统，包括：
温度序列周期化处理模块，用于根据水电机组定子绕组温度变化的特征，将原始温度序列周期化处理；
STL分解模块，用于使用STL分解方法将周期化后的原始温度序列进行分解，获取整体温度序列和周期温度序列；
滤波模块，用于对整体温度序列滤波处理；
加法模块，用于将滤波后的残差以及STL分解方法分解后的残差加入至周期温度序列中，获取周期温度趋势数据；
GRU温度预测模块，其内设置有基于运行状态损失函数的GRU温度预测模型，用于以将机组当前运行阶段的周期温度趋势数据中若干历史温度实际值为输入，输出机组当前运行阶段未来时刻的预测温度值；
误差参数计算模块，用于将机组历史运行阶段的温度变化划分为温度上升阶段和温度稳定阶段，在温度上升阶段根据温度预测点在历史运行阶段中的位置计算历史温度预测值与历史温度实际值的误差，以及计算历史误差的标准差、均值和置信度；
温度预警模块，用于当历史误差的标准差小于预设标准差，则对系统误差进行校正，结合机组当前运行阶段未来时刻的预测温度值，判断是否进行温度预警；当历史误差的标准差大于预设标准差，则结合置信度和当前运行阶段未来时刻的预测温度值，判断是否进行温度预警。
进一步优选地，GRU温度预测模型的损失函数包括均方误差损失函数和运行状态损失函数；运行状态损失函数包括平均温度损失函数、温度上升占比损失函数、最高温度损失函数和中位温度损失函数；
所述平均温度损失函数用于计算单次机组运行过程中实际温度平均值与预测温度平均值之间误差；温度上升占比损失函数用于计算机组当前运行阶段中温度上升所用时间占所有时间的比值；最高温度损失函数用于计算当前运行阶段中实际温度最高值和预测温度最高值的误差；中位数温度损失函数用于计算当前运行阶段中实际温度的中位数值和预测温度中位数值的误差。
进一步优选地，温度预警规则为：
Figure BDA0003854727690000051
其中，
Figure BDA0003854727690000052
为历史运行阶段中温度预测点的标准差；C为温度预测点的置信度；
Figure BDA0003854727690000053
表示当前运行阶段中未来时刻的预测温度值；TH表示阈值；
Figure BDA0003854727690000054
为温度预测点的误差均值。
进一步优选地，计算历史误差的标准差、均值和置信度的方法，包括以下步骤：
计算若干历史运行阶段和当前运行阶段中每个温度预测点在周期温度趋势数据中的位置；
在历史运行阶段中寻找与当前运行阶段中温度预测点相同位置的温度预测点，根据历史运行阶段中温度预测点处的历史温度预测值与历史温度实际值，计算预测误差以及置信度；
计算历史运行阶段中温度预测点的误差均值和标准差。
进一步优选地，历史运行阶段中任意一个温度预测点所处位置为：
Figure BDA0003854727690000055
其中，历史运行阶段中温度上升阶段为t～t+Ti r；t′表示历史运行阶段中温度预测点的时间；
在当前运行阶段中温度预测点所处位置为：
Figure BDA0003854727690000056
其中，t″表示当前运行阶段中温度预测点时间；Tr表示当前运行阶段中的温度上升时间；
Figure BDA0003854727690000057
分别为五个历史运行阶段中的温度上升时间。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，具有以下有益效果：
本发明提供了一种水电机组定子线圈温度预警方法及系统，考虑到气候温度对机组部件温度的影响，利用STL分解方法将周期化后的原始温度序列进行分解，获取整体温度序列和周期温度序列，整体温度序列表征气候温度，在原始温度序列中剔除整体温度序列，提取出受机组运行影响的周期温度趋势数据进行预测和预警，提高了模型训练的效果，同时也降低了环境温度对温度误报警的影响。
本发明提供了一种水电机组定子线圈温度预警方法及系统，其中构建了基于运行状态损失函数的GRU温度预测模型，运行状态损失函数通过计算各运行阶段中温度上升占比、温度变化均值、温度最大值等，反映温度整体特征的量以捕捉温度趋势在机组运行阶段中的关键性特征，从而提高模型的预测效果。
本发明提供了一种水电机组定子线圈温度预警方法及系统，当历史误差的标准差小于预设标准差，则对系统误差进行校正，结合机组当前运行阶段未来时刻的预测温度值，判断是否进行温度预警；当历史误差的标准差大于预设标准差，则结合置信度和当前运行阶段未来时刻的预测温度值，判断是否进行温度预警；温度预警方法充分考虑了GRU温度预测模型在不同温度变化阶段的预测误差对预警效果的影响，提出了置信度和误差校正方法以消除模型预测误差导致的误报警和报警不及时问题，提高了预警的准确性。
附图说明
图1是本发明实施例提供的水电机组定子线圈温度预警流程图；
图2是本发明实施例提供的水电机组定子线圈原始温度序列变化趋势图；
图3是本发明实施例提供的STL分解得到的整体温度变化趋势图；
图4是本发明实施例提供的STL分解得到的周期温度变化趋势图；
图5是本发明实施例提供的运行状态损失函数不同权重α下的预测效果图；
图6是本发明实施例提供的损失函数和mse在不同预测步长下的效果对比；
图7是本发明实施例提供的水电机组定子线圈温度预测与预警效果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明提供了一种水电机组定子线圈温度预警方法及系统，整体技术方案为：根据机组温度变化的特征，对非周期性的初始温度序列进行周期化处理，使用STL分解方法对周期化后的温度序列进行分解，提取出其中的整体温度序列以及与机组运行相关的局部周期温度序列；基于运行状态损失函数的GRU温度预测模型对周期温度趋势数据进行预测，深度挖掘出温度与机组运行状态的内在联系，提出模型预测精度；采用温度预警方法，结合置信度计算方法和误差校正方法以消除不同运行阶段预测误差造成的误报警；其中，整体温度序列表征环境温度对定子线圈温度的影响；周期温度序列表征机组运行阶段对定子线圈的影响。
更为具体地，如图1所示，一方面，本发明提供了一种水电机组定子线圈温度预警方法，包括以下步骤：
S1：根据水电机组定子绕组温度变化的特征，将原始温度序列周期化处理，并使用STL分解方法将周期化后的原始温度序列进行分解，获取整体温度序列和周期温度序列；将整体温度序列进一步滤波处理，获取更加平滑的整体温度趋势数据，同时经滤波后的残差加入到周期温度序列中，获取周期温度趋势数据；
更为具体地，S1中原始温度序列周期化处理及STL分解具体包括以下步骤：
计算机组平均运行时间，具体计算公式如下：
Figure BDA0003854727690000071
其中，Ti表示每次运行的实际时长；
Figure BDA0003854727690000072
表示机组运行的平均时长；N表示机组运行次数；
对每个运行阶段进行处理，将运行时间多于平均运行时间的时间段对应原始温度序列除去，将运行时间段不足的进行温度补齐，补齐的数值即为该运行阶段机组温度达到稳定后的均值温度，则机组运行时长变为
Figure BDA0003854727690000081
将周期化的原始温度序列进行STL分解，分解后得到整体温度序列、周期温度序列和残差值；对整体温度序列进行滤波处理，得到更加平滑的整体温度趋势数据，将周期温度序列和上述分解、滤波两步中得到的残差值相加作为周期温度趋势数据；
S2：建立基于运行状态损失函数的GRU温度预测模型；该模型利用运行状态损失函数捕捉温度在每个开停机运行过程中的关键特征，提高模型的预测效果；
更为具体地，S2中基于运行状态损失函数的GRU温度预测模型，主要结构如下：
GRU温度预测模型的损失函数包括均方误差损失函数mse和运行状态损失函数；其中，运行状态损失函数包括：平均温度损失、温度上升占比损失、最高温度损失和中位温度损失；
平均温度损失用于计算单次机组运行过程中实际温度平均值与预测温度平均值之间误差，其计算公式如下：
Figure BDA0003854727690000082
其中，t表示机组开始运行的时间，y表示实际温度值，
Figure BDA0003854727690000083
表示预测温度值；Ti表示每次运行的实际时长；
温度上升占比损失函数是计算该次运行过程中温度上升所用时间占所有时间的比值，其计算公式如下：
Figure BDA0003854727690000084
其中，tup表示实际温度上升所用时间，
Figure BDA0003854727690000085
表示预测温度上升所用时间；
最高温度损失函数计算该次运行过程中实际温度最高值和预测温度最高值的误差，其计算公式如下：
Figure BDA0003854727690000091
中位数温度损失函数计算该次运行过程中实际温度的中位数值和预测温度中位数值的误差，其计算公式如下：
Figure BDA0003854727690000092
均方误差损失函数mse计算公式如下：
Figure BDA0003854727690000093
由上可知，GRU温度预测模型的总损失函数计算如下：
loss＝mse+α·(loss1+loss2+loss3+loss4)
其中，α表示运行状态损失函数的权重；
S3：将机组历史运行阶段的温度变化划分为两个阶段：温度上升阶段和温度稳定阶段；根据当前运行阶段温度预测点所处的位置计算过去五个运行阶段对应点的历史预测值与历史实际值的误差，计算获取五个历史误差的标准差和均值，同时根据误差计算得到预测值的置信度；根据误差的标准差情况对预测值进行处理，并与报警阈值进行比较，实现机组温度预警；
更为具体地，S3中温度预警的具体步骤如下：
对过去五个历史运行阶段的温度变化进行划分，过程i可划分为温度上升阶段t～t+Ti r和温度稳定阶段t+Ti r～t+Ti r+Ti s，其中，i＝1,2,...,5；温度上升阶段的结束条件为连续的相邻五个点平均温度之间的变化不超过0.1℃，公式表示如下：
Figure BDA0003854727690000094
根据上述划分规则，可以计算历史任意一个温度点所处位置p，历史运行阶段中位置p的计算过程如下：
Figure BDA0003854727690000101
其中，exp表示以e为底的指数函数，t′表示对应预测点的时间；
当前运行阶段的p具体计算方式如下：
Figure BDA0003854727690000102
其中，t″表示当前时刻的预测点时间，Tr表示当前运行阶段中的温度上升时间；由于在当前运行阶段温度上升使其并不能获取到Tr，因此选择相应的历史值进行替代，同时需要不断利用上述阶段划分规则对新的预测点进行阶段划分；
按照上述规则找到过去五个历史运行阶段中与当前运行阶段中对应阶段相同p值的点，计算其置信度公式如下：
Figure BDA0003854727690000103
C＝σ(-2ln(|ep|))
其中，|ep|为预测误差，
Figure BDA0003854727690000104
表示第n个历史运行阶段与预测点相同p值的历史温度实际值和历史温度预测值，σ表示Sigmoid函数，C即为该预测点的置信度；由公式可知，当预测误差均值为0.1时，置信度水平达到0.99；当预测误差均值为0.5时，置信度水平达到0.8；当预测误差均值为1时，置信度水平降低到0.5；
针对历史p值点的预测误差
Figure BDA0003854727690000105
计算其n＝1,2,...,5时的均值为
Figure BDA0003854727690000106
标准差为
Figure BDA0003854727690000107
当标准差较小时，表明预测误差以系统误差居多，需进行误差校正，对于温度序列而言，造成系统性误差的主要原因是温度序列发生变化，通常表现为异常增大趋势，所以当出现系统性误差时，应该适当加大预警力度；当标准差较大时，表明预测误差以随机误差居多，可利用置信度反映预测可靠性，预警规则即为：
Figure BDA0003854727690000111
其中，
Figure BDA0003854727690000112
表示当前的预测温度值，TH表示阈值大小。
另一方面，本发明提供了一种水电机组定子线圈温度预警系统，包括：
温度序列周期化处理模块，用于根据水电机组定子绕组温度变化的特征，将原始温度序列周期化处理；
STL分解模块，用于使用STL分解方法将周期化后的原始温度序列进行分解，获取整体温度序列和周期温度序列；
滤波模块，用于对整体温度序列滤波处理；
加法模块，用于将滤波后的残差以及STL分解方法分解后的残差加入至周期温度序列中，获取周期温度趋势数据；
GRU温度预测模块，其内设置有基于运行状态损失函数的GRU温度预测模型，用于以将机组当前运行阶段的周期温度趋势数据中若干历史温度实际值为输入，输出机组当前运行阶段未来时刻的预测温度值；
误差参数计算模块，用于将机组历史运行阶段的温度变化划分为温度上升阶段和温度稳定阶段，在温度上升阶段根据温度预测点在历史运行阶段中的位置计算历史温度预测值与历史温度实际值的误差，计算历史误差的标准差、均值和置信度；
温度预警模块，用于当历史误差的标准差小于预设标准差，则对系统误差进行校正，结合机组当前运行阶段未来时刻的预测温度值，判断是否进行温度预警；当历史误差的标准差大于预设标准差，则结合置信度和当前运行阶段未来时刻的预测温度值，判断是否进行温度预警。
进一步优选地，GRU温度预测模型的损失函数包括均方误差损失函数和运行状态损失函数；运行状态损失函数包括平均温度损失函数、温度上升占比损失函数、最高温度损失函数和中位温度损失函数；
所述平均温度损失函数用于计算单次机组运行过程中实际温度平均值与预测温度平均值之间误差；温度上升占比损失函数用于计算机组当前运行阶段中温度上升所用时间占所有时间的比值；最高温度损失函数用于计算当前运行阶段中实际温度最高值和预测温度最高值的误差；中位数温度损失函数用于计算当前运行阶段中实际温度的中位数值和预测温度中位数值的误差。
进一步优选地，温度预警规则为：
Figure BDA0003854727690000121
其中，
Figure BDA0003854727690000122
为历史运行阶段中温度预测点的标准差；C为温度预测点的置信度；
Figure BDA0003854727690000123
表示当前运行阶段中未来时刻的预测温度值；TH表示阈值；
Figure BDA0003854727690000124
为温度预测点的误差均值。
进一步优选地，计算历史误差的标准差、均值和置信度的方法，包括以下步骤：
计算若干历史运行阶段和当前运行阶段中每个温度预测点在周期温度趋势数据中的位置；
在历史运行阶段中寻找与当前运行阶段中温度预测点相同位置的温度预测点，根据历史运行阶段中温度预测点处的历史温度预测值与历史温度实际值，计算预测误差以及置信度；
计算历史运行阶段中温度预测点的误差均值和标准差。
进一步优选地，历史运行阶段中任意一个温度预测点所处位置为：
Figure BDA0003854727690000125
其中，历史运行阶段中温度上升阶段为t～t+Ti r；t′表示历史运行阶段中温度预测点的时间；
在当前运行阶段中温度预测点所处位置为：
Figure BDA0003854727690000131
其中，t″表示当前运行阶段中温度预测点时间；Tr表示当前运行阶段中的温度上升时间；
Figure BDA0003854727690000132
分别为五个历史运行阶段中的温度上升时间。
实施例
Step1：该温度数据集时间跨度为2018年1月2日至2018年12月29日，机组运行过程中的定子线圈温度变化原始序列如图2所示，由图可知，在一年之中，定子铁心温度随着环境温度变化有明显的变化趋势，冬季整体温度偏低，夏季整体温度偏高；将原始温度序列进行周期化处理后进行STL分解，将分解后的整体温度序列进行滤波处理得到整体温度趋势数据，如图3所示，提取出的整体温度趋势与原始温度序列中的趋势相符；将STL分解和滤波得到的残差加入周期温度序列中，得到周期温度趋势如图4所示，此时的温度变化仅与机组运行过程相关，周期温度趋势数据整体分布平稳；
Step2：利用基于运行状态损失函数的GRU温度预测模型对温度周期变化趋势进行预测，GRU模型输入序列长度为9，输出长度为1，即使用9个点预测1个点；模型的学习率设置为0.005，训练集与测试集的比例为3:1；运行状态损失函数的权重α对单步预测结果的影响如图5所示，RMSE指标为均方根误差，反映了温度预测值与温度实际值之间的误差，RMSE越小，表明预测精度越高；由图5可知，当权重取值为0.1时，各损失函数均能使预测效果达到最优，且要优于未使用状态损失函数的预测模型；更进一步地，增加模型的预测步长，实现模型的多步预测；单步预测代表使用历史数据预测最新的下一个数据点，本实施例中代表预测一分钟后的点，多步预测中预测步长与预测时长的关系以此类推；由图6可知，步长从1到10时，总损失函数相比于使用单一的mse损失函数均具有优势，能够提高模型的预测精度；
Step3：在多步预测的精度上实现温度预警功能，为了确保预警的精度，本发明的预测步长设置为3，即预测机组3分钟后的温度；如图7所示提供了多步温度预测结果和预警示意图；在达到温度预警值后，由图7可知，在实际报警点约前五分钟的时刻，模型利用历史数据对三分钟后的温度进行预测，得到的预测值满足误差校正的要求，通过误差校正后达到了报警阈值，实现了提前报警；由此可知，本发明的误差预警方法能够避免由于温度数据发生变化而造成的报警不及时问题。
本发明提供了一种水电机组定子线圈温度预警方法及系统，考虑到气候温度对机组部件温度的影响，利用STL分解方法将周期化后的原始温度序列进行分解，获取整体温度序列和周期温度序列，整体温度序列表征气候温度，在原始温度序列中剔除整体温度序列，提取出受机组运行影响的周期温度趋势数据进行预测和预警，提高了模型训练的效果，同时也降低了环境温度对温度误报警的影响。
本发明提供了一种水电机组定子线圈温度预警方法及系统，其中构建了基于运行状态损失函数的GRU温度预测模型，运行状态损失函数通过计算各运行阶段中温度上升占比、温度变化均值、温度最大值等，反映温度整体特征的量以捕捉温度趋势在机组运行阶段中的关键性特征，从而提高模型的预测效果。
本发明提供了一种水电机组定子线圈温度预警方法及系统，当历史误差的标准差小于预设标准差，则对系统误差进行校正，结合机组当前运行阶段未来时刻的预测温度值，判断是否进行温度预警；当历史误差的标准差大于预设标准差，则结合置信度和当前运行阶段未来时刻的预测温度值，判断是否进行温度预警；温度预警方法充分考虑了GRU温度预测模型在不同温度变化阶段的预测误差对预警效果的影响，提出了置信度和误差校正方法以消除模型预测误差导致的误报警和报警不及时问题，提高了预警的准确性。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种水电机组定子线圈温度预警方法，其特征在于，包括以下步骤：
S1：根据水电机组定子绕组温度变化的特征，将原始温度序列周期化处理，并使用STL分解方法将周期化后的原始温度序列进行分解，获取整体温度序列和周期温度序列；
S2：对整体温度序列滤波处理，将滤波后的残差以及STL分解方法分解后的残差加入至周期温度序列中，获取周期温度趋势数据；
S3：将机组当前运行阶段的周期温度趋势数据中若干历史温度实际值输入至基于运行状态损失函数的GRU温度预测模型中，输出机组当前运行阶段未来时刻的预测温度值；
S4：将机组历史运行阶段的温度变化划分为温度上升阶段和温度稳定阶段，确定在温度上升阶段当前运行阶段温度预测点在历史运行阶段中的位置，计算历史温度预测值与历史温度实际值的误差，以及计算历史误差的标准差、均值和置信度；
S5：当历史误差的标准差小于预设标准差，则对系统误差进行校正，结合机组当前运行阶段未来时刻的预测温度值，判断是否进行温度预警；当历史误差的标准差大于预设标准差，则结合置信度和当前运行阶段未来时刻的预测温度值，判断是否进行温度预警。
2.根据权利要求1所述的水电机组定子线圈温度预警方法，其特征在于，将原始温度序列周期化处理的方法为：
根据机组运行次数和每次运行的实际时长计算机组平均运行时间；
若机组的运行时间多于机组平均运行时间，则将机组运行时间多于平均运行时间的时间段对应原始温度序列除去；若机组的运行时间少于机组平均运行时间，则将机组运行时间少于平均运行时间的时间段进行温度补齐，补齐的温度数值为机组运行阶段温度达到稳定后的均值温度。
3.根据权利要求2所述的水电机组定子线圈温度预警方法，其特征在于，所述GRU温度预测模型的损失函数包括均方误差损失函数和运行状态损失函数；运行状态损失函数包括平均温度损失函数、温度上升占比损失函数、最高温度损失函数和中位温度损失函数；
所述平均温度损失函数用于计算单次机组运行过程中实际温度平均值与预测温度平均值之间误差；温度上升占比损失函数用于计算机组当前运行阶段中温度上升所用时间占所有时间的比值；最高温度损失函数用于计算当前运行阶段中实际温度最高值和预测温度最高值的误差；中位数温度损失函数用于计算当前运行阶段中实际温度的中位数值和预测温度中位数值的误差。
4.根据权利要求1至3任一所述的水电机组定子线圈温度预警方法，其特征在于，温度上升阶段的结束条件为连续的相邻五个温度预测点平均温度之间的变化不超过0.1℃。
5.根据权利要求3所述的水电机组定子线圈温度预警方法，其特征在于，计算历史误差的标准差、均值和置信度的方法，包括以下步骤：
计算若干历史运行阶段和当前运行阶段中每个温度预测点在周期温度趋势数据中的位置；
在历史运行阶段中寻找与当前运行阶段中温度预测点相同位置的温度预测点，根据历史运行阶段中温度预测点处的历史温度预测值与历史温度实际值，计算预测误差以及置信度；
计算历史运行阶段中温度预测点的误差均值和标准差。
6.根据权利要求5所述的水电机组定子线圈温度预警方法，其特征在于，历史运行阶段中任意一个温度预测点所处位置为：
Figure FDA0003854727680000021
其中，历史运行阶段中温度上升阶段为t～t+Ti r；t′表示历史运行阶段中温度预测点的时间；
在当前运行阶段中温度预测点所处位置为：
Figure FDA0003854727680000031
其中，t″表示当前运行阶段中温度预测点时间；Tr表示当前运行阶段中的温度上升时间；T1 r,T2 r,...,T5 r分别为五个历史运行阶段中的温度上升时间。
7.根据权利要求5所述的水电机组定子线圈温度预警方法，其特征在于，温度预警规则为：
Figure FDA0003854727680000032
其中，
Figure FDA0003854727680000033
为历史运行阶段中温度预测点的标准差；C为温度预测点的置信度；
Figure FDA0003854727680000034
表示当前运行阶段中未来时刻的预测温度值；TH表示阈值；
Figure FDA0003854727680000035
为温度预测点的误差均值。
8.一种水电机组定子线圈温度预警系统，其特征在于，包括：
温度序列周期化处理模块，用于根据水电机组定子绕组温度变化的特征，将原始温度序列周期化处理；
STL分解模块，用于使用STL分解方法将周期化后的原始温度序列进行分解，获取整体温度序列和周期温度序列；
滤波模块，用于对整体温度序列滤波处理；
加法模块，用于将滤波后的残差以及STL分解方法分解后的残差加入至周期温度序列中，获取周期温度趋势数据；
GRU温度预测模块，其内设置有基于运行状态损失函数的GRU温度预测模型，用于以将机组当前运行阶段的周期温度趋势数据中若干历史温度实际值为输入，输出机组当前运行阶段未来时刻的预测温度值；
误差参数计算模块，用于将机组历史运行阶段的温度变化划分为温度上升阶段和温度稳定阶段，在温度上升阶段根据温度预测点在历史运行阶段中的位置计算历史温度预测值与历史温度实际值的误差，以及计算历史误差的标准差、均值和置信度；
温度预警模块，用于当历史误差的标准差小于预设标准差，则对系统误差进行校正，结合机组当前运行阶段未来时刻的预测温度值，判断是否进行温度预警；当历史误差的标准差大于预设标准差，则结合置信度和当前运行阶段未来时刻的预测温度值，判断是否进行温度预警。
9.根据权利要求8所述的水电机组定子线圈温度预警系统，其特征在于，所述GRU温度预测模型的损失函数包括均方误差损失函数和运行状态损失函数；运行状态损失函数包括平均温度损失函数、温度上升占比损失函数、最高温度损失函数和中位温度损失函数；
所述平均温度损失函数用于计算单次机组运行过程中实际温度平均值与预测温度平均值之间误差；温度上升占比损失函数用于计算机组当前运行阶段中温度上升所用时间占所有时间的比值；最高温度损失函数用于计算当前运行阶段中实际温度最高值和预测温度最高值的误差；中位数温度损失函数用于计算当前运行阶段中实际温度的中位数值和预测温度中位数值的误差。
10.根据权利要求8或9所述的水电机组定子线圈温度预警系统，其特征在于，温度预警规则为：
Figure FDA0003854727680000041
其中，
Figure FDA0003854727680000042
为历史运行阶段中温度预测点的标准差；C为温度预测点的置信度；
Figure FDA0003854727680000043
表示当前运行阶段中未来时刻的预测温度值；TH表示阈值；
Figure FDA0003854727680000044
为温度预测点的误差均值。
CN202211143698.3A 2022-09-20 2022-09-20 一种水电机组定子线圈温度预警方法及系统 Pending CN115496278A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211143698.3A CN115496278A (zh) 2022-09-20 2022-09-20 一种水电机组定子线圈温度预警方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211143698.3A CN115496278A (zh) 2022-09-20 2022-09-20 一种水电机组定子线圈温度预警方法及系统

Publications (1)

Publication Number Publication Date
CN115496278A true CN115496278A (zh) 2022-12-20

Family

ID=84470647

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211143698.3A Pending CN115496278A (zh) 2022-09-20 2022-09-20 一种水电机组定子线圈温度预警方法及系统

Country Status (1)

Country Link
CN (1) CN115496278A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116432542A (zh) * 2023-06-12 2023-07-14 国网江西省电力有限公司电力科学研究院 一种基于误差序列修正的开关柜母排温升预警方法及系统
CN116501581A (zh) * 2023-06-26 2023-07-28 宜宾邦华智慧科技有限公司 一种手机温度监测预警方法
CN117542169A (zh) * 2023-11-07 2024-02-09 国网江苏省电力有限公司镇江供电分公司 基于大数据分析的自动化设备温度异常预警方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116432542A (zh) * 2023-06-12 2023-07-14 国网江西省电力有限公司电力科学研究院 一种基于误差序列修正的开关柜母排温升预警方法及系统
CN116432542B (zh) * 2023-06-12 2023-10-20 国网江西省电力有限公司电力科学研究院 一种基于误差序列修正的开关柜母排温升预警方法及系统
CN116501581A (zh) * 2023-06-26 2023-07-28 宜宾邦华智慧科技有限公司 一种手机温度监测预警方法
CN116501581B (zh) * 2023-06-26 2023-08-25 宜宾邦华智慧科技有限公司 一种手机温度监测预警方法
CN117542169A (zh) * 2023-11-07 2024-02-09 国网江苏省电力有限公司镇江供电分公司 基于大数据分析的自动化设备温度异常预警方法

Similar Documents

Publication Publication Date Title
CN115496278A (zh) 一种水电机组定子线圈温度预警方法及系统
CN110570012B (zh) 一种基于Storm的电厂生产设备故障预警方法及系统
CN110266031B (zh) 储能并网充放电量控制方法、装置、服务器及存储介质
CN108334997B (zh) 基于支撑故障事件约束机组组合的备用优化方法和装置
US7050935B1 (en) Method for assessing the reliability of technical systems
CN113657668A (zh) 基于lstm网络的电力负荷预测方法及系统
CN115453356B (zh) 一种动力设备运行状态监测分析方法、系统、终端及介质
CN111669123A (zh) 用于对光伏组串进行故障诊断的方法和装置
CN113868953B (zh) 工业系统中多机组运行优化方法、装置、系统及存储介质
CN109783873B (zh) 轴温异常状态的预测方法和装置
CN111008727A (zh) 一种配电台区负荷预测方法及装置
CN110909958A (zh) 一种计及光伏并网功率的短期负荷预测方法
CN111860611A (zh) 一种基于马尔科夫决策的配电系统弹性策略构建的方法
CN118226366A (zh) 电压互感器的在线监测方法、装置、设备及存储介质
CN111799798B (zh) 一种提高未来态潮流计算结果准确性的方法及系统
CN112801533A (zh) 一种考虑决策依赖不确定的电力系统运行可靠性评估方法
CN116502761A (zh) 一种基于大数据节约能耗的监控方法及系统
CN116662841A (zh) 基于插值法与时序智能法的多源数据修正方法及系统
CN113221248B (zh) 一种基于pf-garch模型的船舶系统设备状态参数预测方法
CN114186756A (zh) 一种配电网终端蓄电池储能容量预测方法
CN103094927A (zh) 考虑风电出力不确定性与半绝对离差风险的机组组合方法及系统
CN112257997A (zh) 一种基于神经网络的pcs转化率计算方法及系统
CN118332478B (zh) 基于改进eemd算法的光伏功率短期分段预测方法及系统
CN114118499B (zh) 一种短期风电功率预测方法及装置
CN117332901B (zh) 采用分层时间聚合策略的新能源小时间尺度功率预测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination