CN112487605A - 一种水轮发电机定子铁心状态确定方法 - Google Patents
一种水轮发电机定子铁心状态确定方法 Download PDFInfo
- Publication number
- CN112487605A CN112487605A CN202011168182.5A CN202011168182A CN112487605A CN 112487605 A CN112487605 A CN 112487605A CN 202011168182 A CN202011168182 A CN 202011168182A CN 112487605 A CN112487605 A CN 112487605A
- Authority
- CN
- China
- Prior art keywords
- temperature
- state
- residual error
- xobs
- iron core
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Manufacture Of Motors, Generators (AREA)
Abstract
本发明涉及一种水轮发电机定子铁心状态确定方法。本发明结合发电机监测传感器安装、数据采集传输,水轮发电机铁心故障特征,以多变量状态评估算法为基础开发水轮发电机铁心状态评估。通过特征量反映状态,能在定子铁心状态特征量温度到达阀值告警前及时发现动态过程中的定子异常,能比较同一机组不同阶段状态和同型号不同机组间状态,能定量评价机组检修前后的状态优劣、同型机组间状态差异。
Description
技术领域
本发明涉及水轮水轮发电机领域,尤其涉及一种水轮发电机定子铁心状态确定方法。
背景技术
目前流域大型水轮发电机定子铁心主要是通过在背部上中下各部位埋设的测温电阻监测温度,在端部安装垂直、水平振动传感器监测垂直、水平振动,在监控系统中设置告警值,在温度振动等越限时发出告警。但在中低负荷运行时温度一般较低离告警阈值较远,通常不容易发现潜在的缺陷,然而机组运行大多数时候都不是满载。另外目前没有较客观的方法评价机组检修前后的阶段性状态变化、同型机组间状态差异,不能为状态检修提供数据支撑。如何综合发电机监测传感器安装、数据采集传输,水轮发电机定子铁心故障特点,以多变量状态评估算法为基础开发水轮发电机铁心状态评估,是值得研究的。
发明内容
为了解决上述问题,本发明提出一种水轮发电机定子铁心状态确定方法,通过特征量反映状态,能在定子铁心状态特征量温度到达阀值告警前及时发现动态过程中的定子异常,能比较同一机组不同阶段状态和同型号不同机组间状态差异。
本发明的技术方案具体如下:
一种水轮发电机定子铁心状态确定方法,包括如下步骤:
步骤(1)、构建水轮发电机正常状态样本集矩阵
将铁心温度作为铁心的状态特征目标变量,其余发电机运行工况变量有功、无功、绕组各测点温度、冷却器冷风各测点温度、冷却器热风各测点温度作为自变量,所有变量1,2,3,...m编号即选择铁心状态相关的变量为m 维向量X;
以时间轴T为坐标在时间轴上选择1,2,3,...n个点截取n个样本构成 n行m列的矩阵D;
在选择样本1到n时,发电机正常且n足够大能涵盖所有发电机正常状态;
X=[x1,x2,x3,...,xm-1,xm];
x1到xm依次分别是有功、无功、铁芯温度1号、铁芯温度2号....绕组温度1号、绕组温度2号、....冷却器冷风温度1号、冷却器冷风温度2号、.... 冷却器热风温度1号、冷却器热风温度2号;
D=[X(1),X(2),X(3),....X(n-1),X(n)];
D由来自不同时间不同工况的n个样本向量组成n*m矩阵;构成D的 X向量足够丰富能涵盖正常状态下的各种工况;
步骤(2)、构建相似度权值矩阵
从监控系统导出与样本集相同的具有m个特征维度的观测样本集Xobs,通过下式计算观测样本集Xobs与正常状态样本集矩阵D的相似度权值矩阵 W:
步骤(3)、计算预测值
通过下式计算基于正常状态的观测样本集Xobs的预测值集Xest:
其中·表示两向量的内积(点乘/数量积);
步骤(4)、计算差值
通过下式计算观测样本集Xobs与其预测值集Xest的差值集C:
C=Xobs-Xest;
步骤(5)、绘制趋势曲线
抽取差值C中状态特征目标变量铁心温度值在时间轴上绘制曲线,即得到铁心状态的趋势曲线。
进一步地,步骤(2)中,观测样本集Xobs是1个或多个向量;
Xobs=[Xo(1),Xo(2),...,Xo(k)];
W(i)=[wi1,wi2,wi3,...,win-1,win];
i=1,2,3,...,k则
W=[W(1),W(2),W(3),.....W(k)];
K大于等于1即可,即观测样本是1个或多个。
进一步地,步骤(4)中,差值C由k个有m维度特征的向量组成,与观测样本集Xobs形状一致,均为k行m列。
进一步地,k维度为时间轴。
进一步地,步骤(4)中,差值C(i)=[ci1,ci2,ci3,...,cim-1,cim],i=1,2,3...k 则C=[C(1),C(2),C(3),.....C(k)];
ci1,ci2,ci3,...,cim-1,cim分别对应有功残差、无功残差、铁芯温度1 号残差、铁芯温度2号残差....绕组温度1号残差、绕组温度2号残差、....冷却器冷风温度1号残差、冷却器冷风温度2号残差、....冷却器热风温度1号残差、冷却器热风温度2号残差。
与现有技术相比,本发明的有益效果如下:
1.本发明将发电机正常数据建立样本集,运行中不断采集实时数据计算,通过真实值与基于正常状态的预测值之差进行状态监视,实现实时趋势预警。
2.本发明以发电机修前数据建立样本集,将修后数据作为观测数据输入,通过真实值与基于修前状态的预测值之差进行状态评估,实现检修前后状态比较。
3.本发明以一台机的数据建立样本集,将另一台机数据作为观测数据输入,实现不同机组间状态比较。
附图说明
图1为评估时段2019年5月至2020年4月运行工况调整变化图,“+”标记线为有功调整、“o”标记线为无功调整散点图;
图2为2019年5月至2020年4月定子铁心3号测点温度随运行工况调整的变化情况;图上部分“+”标记线为实际温度、“o”标记线为基于2-4月状态预测温度,图上部分“o”标记线为偏差即该测点温度在同等工况下与2-4 月相比的变化趋势;
图3为2019年5月至2020年4月定子铁心平均温度随运行工况调整的变化情况。图上部分“+”标记线为实际平均温度、“o”标记线为基于2-4月状态预测平均温度,图上部分“o”标记线为偏差即铁心整体温度在同等工况下与2-4月相比的变化趋势。
具体实施方式
下面将结合本申请实施例中的附图,对实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
除非另外定义,本申请实施例中使用的技术术语或者科学术语应当为所属领域内具有一般技能的人士所理解的通常意义。本实施例中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性,而只是用来区分不同的组成部分。“包括”或者“包含”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同,而不排除其他元件或者物件。“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。“上”、“下”、“左”、“右”、“横”以及“竖”等仅用于相对于附图中的部件的方位而言的,这些方向性术语是相对的概念,它们用于相对于的描述和澄清,其可以根据附图中的部件所放置的方位的变化而相应地发生变化。
本实施例以某水电厂650MW水轮发电机定子铁心阶段评估为例。
以该机2019年2月至4月工况数据为基础选择样本建立正常状态样本集,评估该机2019年5月至2020年4月在2019年2-4月基础上的状态的变化趋势。
本实施例的水轮发电机定子铁心状态确定方法,按以下步骤进行:
步骤(1)、建立水轮发电机铁心正常状态样本集D。
从该机组2019年2月至2019年4月监控数据128159个样本中,选出具有代表性的4603个样本,构成D,兼顾计算速度和模型准确度,D中样本即不重复又能充分体现3个月的各种正常运行工况。样本维度245个包含机组有功功率P、机组无功功率Q、定子压指温度16个、定子铁心温度54个、定子绕组温度144个、空冷器热风温度14个、空冷器冷风温度14个。
X=[x1,x2,x3,...,x244,x245];
D=[X(1),X(2),X(3),....,X(4602),X(4603)];
X为245维度样本向量,D为4603个样本向量构成的样本矩阵,矩阵大小4603*245。
步骤(2)、构建相似度权值矩阵
从监控系统导出该机2019年5月至2020年4月数据每间隔40分选择1 个样本共12924个245维度的实时样本Xobs(12924*45),通过公式(1) (2)计算观测样集Xobs与正常状态样本集矩阵D的相似度权值矩阵W。Xobs 中的12924(k)个样本依次与D的4603(n)个向量分别计算相似度权值构成12924*4603的矩阵W。
Xobs=[Xo(1),Xo(2),...,Xo(12924)];
W的矩阵计算公式如下:
n=1,2,3,...4603;则:
W(i)=[wi1,wi2,wi3,...,wi4602,wi4603];
i=1,2,3,...,12924;则:
W=[W(1),W(2),W(3),.....W(12924)]。
步骤(3)、计算预测值
通过下式计算基于正常状态的观测样Xobs预测值Xest。
Xest=D·W;
步骤(4)、计算差值
4.通过公式(4)计算观测样本Xobs与其预测值Xest的差值C,C由12924(k) 个有245(m)维度特征的向量组成,与Xobs形状一致均为12924行245列。
步骤(5)、绘制趋势曲线
抽取C中3号铁心温度列的12924个温度差值绘制曲线,对应日期时间 2019年5月1日至2020年4月23日作为横坐标,就可得到3号测点位置铁心状态的趋势曲线图2下半部分图中“o”标记线。将全部铁心温度测点求平均后就可画出铁心整体状态的趋势曲线图3下半部分图中“o”标记线。
实施效果:
从图1、图2可看出该评估方法可排除运行工况变化引起的铁心温度正常波动,实现动态过程的趋势评估。
图3中图上半部分“+”标记线为发电机不同运行工况下观测的铁心温度均值,“o”标记线为在相应工况下基于正常状态预测的铁心温度均值,图3中图下半部分“o”标记线为整体温度变化趋势线,该实例中为一条水平线即定子铁心从2019年5月至2020年4月近1年,趋势无变化、状态稳定。
尽管已经对本发明的技术方案做了较为详细的阐述和列举,应当理解,对于本领域技术人员来说,对上述实施例做出修改或者采用等同的替代方案,这对本领域的技术人员而言是显而易见,在不偏离本发明精神的基础上所做的这些修改或改进,均属于本发明要求保护的范围。
Claims (5)
1.一种水轮发电机定子铁心状态确定方法,其特征在于:包括如下步骤:步骤(1)、构建水轮发电机正常状态样本集矩阵
将铁心温度作为铁心的状态特征目标变量,其余发电机运行工况变量有功、无功、绕组各测点温度、冷却器冷风各测点温度、冷却器热风各测点温度作为自变量,所有变量1,2,3,...m编号即选择铁心状态相关的变量为m维向量X;
以时间轴T为坐标在时间轴上选择1,2,3,...n个点截取n个样本构成n行m列的矩阵D;
在选择样本1到n时,发电机正常且n足够大能涵盖所有发电机正常状态;
X=[x1,x2,x3,...,xm-1,xm];
x1到xm依次分别是有功、无功、铁芯温度1号、铁芯温度2号....绕组温度1号、绕组温度2号、....冷却器冷风温度1号、冷却器冷风温度2号、....冷却器热风温度1号、冷却器热风温度2号;
D=[X(1),X(2),X(3),....X(n-1),X(n)];
D由来自不同时间不同工况的n个样本向量组成n*m矩阵;构成D的X向量足够丰富能涵盖正常状态下的各种工况;
步骤(2)、构建相似度权值矩阵
从监控系统导出与样本集相同的具有m个特征维度的观测样本集Xobs,通过下式计算观测样本集Xobs与正常状态样本集矩阵D的相似度权值矩阵W:
步骤(3)、计算预测值
通过下式计算基于正常状态的观测样本集Xobs的预测值集Xest:
其中·表示两向量的内积(点乘/数量积);
步骤(4)、计算差值
通过下式计算观测样本集Xobs与其预测值集Xest的差值集C:
C=Xobs-Xest;
步骤(5)、绘制趋势曲线
抽取差值C中状态特征目标变量铁心温度值在时间轴上绘制曲线,即得到铁心状态的趋势曲线。
2.根据权利要求1所述的水轮发电机定子铁心状态确定方法,其特征在于:步骤(2)中,观测样本集Xobs是1个或多个向量;
Xobs=[Xo(1),Xo(2),...,Xo(k)];
W(i)=[wi1,wi2,wi3,...,win-1,win];
i=1,2,3,...,k则
W=[W(1),W(2),W(3),.....W(k)];
K大于等于1即可,即观测样本是1个或多个。
3.根据权利要求1所述的水轮发电机定子铁心状态确定方法,其特征在于:步骤(4)中,差值C由k个有m维度特征的向量组成,与观测样本集Xobs形状一致,均为k行m列。
4.根据权利要求3所述的水轮发电机定子铁心状态确定方法,其特征在于:k维度为时间轴。
5.根据权利要求1所述的水轮发电机定子铁心状态确定方法,其特征在于:
步骤(4)中,差值C(i)=[ci1,ci2,ci3,...,cim-1,cim],i=1,2,3...k则
C=[C(1),C(2),C(3),.....C(k)];
ci1,ci2,ci3,...,cim-1,cim分别对应有功残差、无功残差、铁芯温度1号残差、铁芯温度2号残差....绕组温度1号残差、绕组温度2号残差、....冷却器冷风温度1号残差、冷却器冷风温度2号残差、....冷却器热风温度1号残差、冷却器热风温度2号残差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011168182.5A CN112487605A (zh) | 2020-10-27 | 2020-10-27 | 一种水轮发电机定子铁心状态确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011168182.5A CN112487605A (zh) | 2020-10-27 | 2020-10-27 | 一种水轮发电机定子铁心状态确定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112487605A true CN112487605A (zh) | 2021-03-12 |
Family
ID=74927209
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011168182.5A Pending CN112487605A (zh) | 2020-10-27 | 2020-10-27 | 一种水轮发电机定子铁心状态确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112487605A (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103383433A (zh) * | 2013-07-03 | 2013-11-06 | 中国人民解放军海军工程大学 | 船舶发电机定子铁心的状态监测和故障预警方法 |
CN105827067A (zh) * | 2016-06-14 | 2016-08-03 | 湖南德益伟节能科技有限公司 | 一种磁矩阵无铁芯电动机的内腔通风散热系统装置 |
-
2020
- 2020-10-27 CN CN202011168182.5A patent/CN112487605A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103383433A (zh) * | 2013-07-03 | 2013-11-06 | 中国人民解放军海军工程大学 | 船舶发电机定子铁心的状态监测和故障预警方法 |
CN105827067A (zh) * | 2016-06-14 | 2016-08-03 | 湖南德益伟节能科技有限公司 | 一种磁矩阵无铁芯电动机的内腔通风散热系统装置 |
Non-Patent Citations (1)
Title |
---|
杨增杰等: "基于非线性状态估计的水轮发电机定子铁心状态评估", 《2019年云、贵、川、湘、桂、粤、青七省(区)水电站运行检修技术交流研讨会论文集》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106682814B (zh) | 一种基于故障知识库的风电机组故障智能诊断方法 | |
Zhang et al. | Machine tool thermal error modeling and prediction by grey neural network | |
CN110738360A (zh) | 一种设备剩余寿命预测方法及系统 | |
CN109976147B (zh) | 一种基于智能学习的大体积混凝土温控方法 | |
CN110006552B (zh) | 一种机组设备温度异常检测方法 | |
CN103614498A (zh) | 一种高炉炉壁三维温度场重建方法及计算机监控系统 | |
CN106597574A (zh) | 一种基于时变云模型的天气温度预测方法及装置 | |
CN116593157A (zh) | 少样本下基于匹配元学习的复杂工况齿轮故障诊断方法 | |
CN115841004A (zh) | 基于多维数据的带钢热轧过程力学性能软测量方法及装置 | |
CN104865944B (zh) | 基于pca‑lssvm的气分装置控制系统性能评估方法 | |
CN116227357A (zh) | 一种风电场尾流预测方法、模型及其系统 | |
CN113156822B (zh) | 基于Mist-边缘-雾-云计算的热误差预测系统及热误差补偿系统 | |
CN112487605A (zh) | 一种水轮发电机定子铁心状态确定方法 | |
Ozsoy et al. | Optimum scheduling of a hot rolling process by nonlinear programming | |
CN117473439A (zh) | 一种无缝钢管连轧生产过程监测与异常追溯方法及系统 | |
CN104564180B (zh) | 汽轮机转子应力在线实时监测系统 | |
CN114239708A (zh) | 一种基于质量控制图理论的燃机异常检测方法 | |
CN114326682A (zh) | 基于比例积分观测器的注塑机冷凝控制系统故障监测方法 | |
Yang et al. | Economic design of the integrated multivariate EPC and multivariate SPC charts | |
CN107654264B (zh) | 滑压曲线获取方法、系统以及压力调节方法和系统 | |
CN111695300A (zh) | 一种汽轮机阀门流量拟合方法及系统 | |
Chen et al. | An improved JITL method for soft sensing of multimodal industrial processes for search efficiency | |
Zhang et al. | Attention-based stacked supervised poisson autoencoders for defects prediction in casting-rolling process | |
CN117407773B (zh) | 基于数字孪生的风机叶片结冰状态预测方法、系统及设备 | |
JP6904473B1 (ja) | モデル作成支援装置、モデル作成支援方法及びプログラム |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210312 |