CN114718669A - 汽轮机调速系统执行机构参数辨识误差的分析方法及装置 - Google Patents
汽轮机调速系统执行机构参数辨识误差的分析方法及装置 Download PDFInfo
- Publication number
- CN114718669A CN114718669A CN202210514864.XA CN202210514864A CN114718669A CN 114718669 A CN114718669 A CN 114718669A CN 202210514864 A CN202210514864 A CN 202210514864A CN 114718669 A CN114718669 A CN 114718669A
- Authority
- CN
- China
- Prior art keywords
- sample
- time
- standard deviation
- actual measurement
- determining
- 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.)
- Granted
Links
- 230000001105 regulatory effect Effects 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 61
- 230000007246 mechanism Effects 0.000 title claims abstract description 55
- 238000004458 analytical method Methods 0.000 claims abstract description 24
- 238000004088 simulation Methods 0.000 claims description 146
- 238000005259 measurement Methods 0.000 claims description 129
- 230000000630 rising effect Effects 0.000 claims description 52
- 238000005070 sampling Methods 0.000 claims description 50
- 238000010606 normalization Methods 0.000 claims description 36
- 238000004590 computer program Methods 0.000 claims description 20
- 238000003860 storage Methods 0.000 claims description 11
- 230000001174 ascending effect Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 21
- 238000010586 diagram Methods 0.000 description 16
- 238000004891 communication Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 14
- 238000012545 processing Methods 0.000 description 14
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 239000000872 buffer Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000003303 reheating Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
- F01D21/003—Arrangements for testing or measuring
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/32—Circuit design at the digital level
- G06F30/33—Design verification, e.g. functional simulation or model checking
- G06F30/3308—Design verification, e.g. functional simulation or model checking using simulation
- G06F30/3312—Timing analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/32—Circuit design at the digital level
- G06F30/333—Design for testability [DFT], e.g. scan chain or built-in self-test [BIST]
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05D—INDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
- F05D2260/00—Function
- F05D2260/81—Modelling or simulation
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05D—INDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
- F05D2270/00—Control
- F05D2270/01—Purpose of the control system
- F05D2270/02—Purpose of the control system to control rotational speed (n)
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mechanical Engineering (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Control Of Turbines (AREA)
Abstract
本申请提供一种汽轮机调速系统执行机构参数辨识误差的分析方法及装置,涉及电力设备性能检测领域,包括:根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差;根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间;根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差。本申请能够根据预设的滑动时间窗口参数确定汽轮机调速系统执行机构的参数辨识误差。
Description
技术领域
本申请涉及电力设备性能检测领域,具体是一种汽轮机调速系统执行机构参数辨识误差的分析方法及装置。
背景技术
通过分析汽轮机调速系统的参数,建立汽轮机调速系统的数学模型,可系统地分析在各种扰动条件下,电网频率响应及负荷响应的变化曲线,对电网稳定性的分析具有重要的实用价值。汽轮机调速系统模型是电网最重要的数据模型之一,其真实性直接影响着电网的仿真精度。
目前,对汽轮机调速系统执行机构的辨识结果进行误差分析,通常需要分析人员在实测曲线和仿真曲线上手动进行标点,采用此方法进行数据处理的过程繁杂,处理效率低且受主观判断的影响较大。
发明内容
针对现有技术中的问题,本申请提供一种汽轮机调速系统执行机构参数辨识误差的分析方法及装置,能够根据预设的滑动时间窗口参数确定汽轮机调速系统执行机构的参数辨识误差。
为解决上述技术问题,本申请提供以下技术方案:
第一方面,本申请提供一种汽轮机调速系统执行机构参数辨识误差的分析方法,包括:
根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差;
根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间;
根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差。
进一步地,所述滑动时间窗口参数包括:所述滑动时间窗口的宽度及数量;所述参数样本包括:仿真样本及实测样本;所述样本标准差包括:仿真样本标准差及实测样本标准差;所述根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差,包括:
根据所述滑动时间窗口的宽度及数量对所述仿真样本进行分析,得到各所述滑动时间窗口对应的仿真样本标准差;
根据所述滑动时间窗口的宽度及数量对所述实测样本进行分析,得到各所述滑动时间窗口对应的实测样本标准差。
进一步地,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间,包括:
对各所述滑动时间窗口对应的仿真样本标准差进行归一化,得到各所述滑动时间窗口对应的仿真样本归一化标准差;
按照所述仿真样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,首个小于所述阶跃判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的起始阶跃时间;
按照所述仿真样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,最后一个大于所述稳态判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的稳态进入时间。
进一步地,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间,包括:
对各所述滑动时间窗口对应的实测样本标准差进行归一化,得到各所述滑动时间窗口对应的实测样本归一化标准差;
按照所述实测样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,首个小于所述阶跃判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的起始阶跃时间;
按照所述实测样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,最后一个大于所述稳态判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的稳态进入时间。
进一步地,所述参数辨识误差包括:上升参数辨识误差及调节参数辨识误差;所述根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差,包括:
根据所述起始阶跃时间及稳态进入时间确定所述参数样本的阶跃量;
根据所述参数样本的阶跃量确定所述参数样本的上升时间偏差值及调节时间偏差值;
比较预设的上升时间偏差允许值与所述上升时间偏差值,得到所述上升参数辨识误差;
比较预设的调节时间偏差允许值与所述调节时间偏差值,得到所述调节参数辨识误差。
进一步地,所述根据所述起始阶跃时间及稳态进入时间确定所述参数样本的阶跃量,包括:
根据所述仿真样本的起始阶跃时间及所述仿真样本的稳态进入时间确定所述仿真样本的阶跃量;
根据所述实测样本的起始阶跃时间及所述实测样本的稳态进入时间确定所述实测样本的阶跃量。
进一步地,所述根据所述参数样本的阶跃量确定所述参数样本的上升时间偏差值及调节时间偏差值,包括:
根据所述仿真样本的阶跃量确定所述仿真样本的第一上升时间偏差值及第一调节时间偏差值;
根据所述实测样本的阶跃量确定所述实测样本的第二上升时间偏差值及第二调节时间偏差值;
根据所述第一上升时间偏差值及所述第二上升时间偏差值确定所述上升时间偏差值;
根据所述第一调节时间偏差值及所述第二调节时间偏差值确定所述调节时间偏差值。
第二方面,本申请提供一种汽轮机调速系统执行机构参数辨识误差的分析装置,包括:
样本标准差确定单元,用于根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差;
阶跃稳态时间确定单元,用于根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间;
参数辨识误差确定单元,用于根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差。
进一步地,所述滑动时间窗口参数包括:所述滑动时间窗口的宽度及数量;所述参数样本包括:仿真样本及实测样本;所述样本标准差包括:仿真样本标准差及实测样本标准差;所述样本标准差确定单元,包括:
仿真标准差确定模块,用于根据所述滑动时间窗口的宽度及数量对所述仿真样本进行分析,得到各所述滑动时间窗口对应的仿真样本标准差;
实测标准差确定模块,用于根据所述滑动时间窗口的宽度及数量对所述实测样本进行分析,得到各所述滑动时间窗口对应的实测样本标准差。
进一步地,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述阶跃稳态时间确定单元,包括:
仿真归一化标准差确定模块,用于对各所述滑动时间窗口对应的仿真样本标准差进行归一化,得到各所述滑动时间窗口对应的仿真样本归一化标准差;
仿真起始阶跃时间确定模块,用于按照所述仿真样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,首个小于所述阶跃判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的起始阶跃时间;
仿真稳态进入时间确定模块,用于按照所述仿真样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,最后一个大于所述稳态判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的稳态进入时间。
进一步地,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述阶跃稳态时间确定单元,包括:
实测归一化标准差确定模块,用于对各所述滑动时间窗口对应的实测样本标准差进行归一化,得到各所述滑动时间窗口对应的实测样本归一化标准差;
实测起始阶跃时间确定模块,用于按照所述实测样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,首个小于所述阶跃判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的起始阶跃时间;
实测稳态进入时间确定模块,用于按照所述实测样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,最后一个大于所述稳态判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的稳态进入时间。
进一步地,所述参数辨识误差包括:上升参数辨识误差及调节参数辨识误差;所述参数辨识误差确定单元,包括:
阶跃量确定模块,用于根据所述起始阶跃时间及稳态进入时间确定所述参数样本的阶跃量;
偏差值确定模块,用于根据所述参数样本的阶跃量确定所述参数样本的上升时间偏差值及调节时间偏差值;
上升参数辨识误差确定模块,用于比较预设的上升时间偏差允许值与所述上升时间偏差值,得到所述上升参数辨识误差;
调节参数辨识误差确定模块,用于比较预设的调节时间偏差允许值与所述调节时间偏差值,得到所述调节参数辨识误差。
进一步地,所述阶跃量确定模块,包括:
仿真阶跃量确定子模块,用于根据所述仿真样本的起始阶跃时间及所述仿真样本的稳态进入时间确定所述仿真样本的阶跃量;
实测阶跃量确定子模块,用于根据所述实测样本的起始阶跃时间及所述实测样本的稳态进入时间确定所述实测样本的阶跃量。
进一步地,所述偏差值确定模块,包括:
仿真偏差值确定子模块,用于根据所述仿真样本的阶跃量确定所述仿真样本的第一上升时间偏差值及第一调节时间偏差值;
实测偏差值确定子模块,用于根据所述实测样本的阶跃量确定所述实测样本的第二上升时间偏差值及第二调节时间偏差值;
上升时间偏差值确定子模块,用于根据所述第一上升时间偏差值及所述第二上升时间偏差值确定所述上升时间偏差值;
调节时间偏差值确定子模块,用于根据所述第一调节时间偏差值及所述第二调节时间偏差值确定所述调节时间偏差值。
第三方面,本申请提供一种电子设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现所述汽轮机调速系统执行机构参数辨识误差的分析方法的步骤。
第四方面,本申请提供一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现所述汽轮机调速系统执行机构参数辨识误差的分析方法的步骤。
第五方面,本申请提供一种计算机程序产品,包括计算机程序/指令,该计算机程序/指令被处理器执行时实现所述汽轮机调速系统执行机构参数辨识误差的分析方法的步骤。
针对现有技术中的问题,本申请提供的汽轮机调速系统执行机构参数辨识误差的分析方法及装置,能够基于改进的滑动窗口法对汽轮机调速系统执行机构参数辨识的误差进行分析与计算,避免了传统手动标点法中人为判断的影响,计算结果更加准确可靠。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例中汽轮机调速系统执行机构参数辨识误差的分析方法的流程图之一;
图2为本申请实施例中得到对应的样本标准差的流程图;
图3为本申请实施例中确定参数样本的起始阶跃时间及稳态进入时间的流程图之一;
图4为本申请实施例中确定参数样本的起始阶跃时间及稳态进入时间的流程图之二;
图5为本申请实施例中确定汽轮机调速系统执行机构的参数辨识误差的流程图;
图6为本申请实施例中确定参数样本的阶跃量的流程图;
图7为本申请实施例中确定参数样本的上升时间偏差值及调节时间偏差值的流程图;
图8为本申请实施例中汽轮机调速系统执行机构参数辨识误差的分析装置的结构图;
图9为本申请实施例中样本标准差确定单元的结构图;
图10为本申请实施例中阶跃稳态时间确定单元的结构图;
图11为本申请实施例中阶跃稳态时间确定单元的结构图;
图12为本申请实施例中参数辨识误差确定单元的结构图;
图13为本申请实施例中阶跃量确定模块的结构图;
图14为本申请实施例中偏差值确定模块的结构图;
图15为本申请实施例中的电子设备的结构示意图;
图16为本申请实施例中执行机构阶跃响应特性示例曲线示意图;
图17为本申请实施例中某300MW机组GV1参数辨识结果示意图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
一实施例中,参见图1,为了能够根据预设的滑动时间窗口参数确定汽轮机调速系统执行机构的参数辨识误差,本申请提供一种汽轮机调速系统执行机构参数辨识误差的分析方法,包括:
S101:根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差;
可以理解的是,汽轮机调速系统(下文简称“系统”)的执行机构辨识结果一般包含实测曲线及仿真曲线。由于实测曲线及仿真曲线均符合图16中示例曲线的基本样式,可以认为实测曲线与仿真曲线的上升时间tup与调节时间ts满足相同的计算方法。
具体地,一实施例中,参见图2,所述滑动时间窗口参数包括:所述滑动时间窗口的宽度及数量;所述参数样本包括:仿真样本及实测样本;所述样本标准差包括:仿真样本标准差及实测样本标准差;所述根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差,包括:根据所述滑动时间窗口的宽度及数量对所述仿真样本进行分析,得到各所述滑动时间窗口对应的仿真样本标准差(S201);根据所述滑动时间窗口的宽度及数量对所述实测样本进行分析,得到各所述滑动时间窗口对应的实测样本标准差(S202)。
步骤S201至步骤S202采用滑动窗口法实现:首先,在时序上分别对仿真样本及实测样本进行划分,使仿真样本及实测样本分别形成多个连续的时间窗口。通过计算样本数据的标准差,进而获取到样本数据的归一化标准差的变化情况,从而确定出汽轮机调速系统的起始阶跃时间及稳态进入时间,标准差的计算公式如下:
其中,SSC(Steady State Criteria)为样本标准差;N为时窗宽度,t为当前样本时刻。仿真样本与实测样本分别有其自身对应的样本标准差SSC,但计算公式均为公式(2-1),具体计算时,只需根据样本类型(仿真样本或实测样本)选取对应的数据(x)代入公式进行计算即可。在对仿真样本进行处理时的时窗宽度N及当前样本时刻t的取值应与对实测样本进行处理时相同。
具体实施时,依据时窗公式(公式2-1)选择系统检测信号样本。x为选择用于判断准稳态特性的变量,需针对研究对象选择合适的变量。对于汽轮机执行机构参数辨识结果误差判断推荐采用高压缸调节门开度反馈。时窗宽度N的选取与过程的惯性延迟相关,通常,可取为10s到20s(样本取样时间间隔为0.1s),在本申请实施例中,选取N为14s,但本申请不以此为限。
S102:根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间;
可以理解的是,在计算得到时窗内样本数据的标准差后,需要对样本标准差进行归一化处理;该归一化处理同样也是针对仿真样本及实测样本分别进行,但处理方法(所涉及的公式算法)相同。
一实施例中,参见图3,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间,包括:对各所述滑动时间窗口对应的仿真样本标准差进行归一化,得到各所述滑动时间窗口对应的仿真样本归一化标准差(S301);按照所述仿真样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,首个小于所述阶跃判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的起始阶跃时间(S302);按照所述仿真样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,最后一个大于所述稳态判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的稳态进入时间(S303)。
一实施例中,参见图4,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间,包括:对各所述滑动时间窗口对应的实测样本标准差进行归一化,得到各所述滑动时间窗口对应的实测样本归一化标准差(S401);按照所述实测样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,首个小于所述阶跃判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的起始阶跃时间(S402);按照所述实测样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,最后一个大于所述稳态判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的稳态进入时间(S403)。
具体实施时,在计算得到时窗内样本数据的标准差后,可以计算得到样本数据标准差在[0,1]之间的分布情况,然后根据设定的判别阈值(也称样本阶跃阈值)来确定阶跃信号发生的时间点(也称起始阶跃时间)及汽轮机调速系统再次趋稳的时间点(也称稳态进入时间),如式(2-2)。
其中,判别阈值θ为无量纲量,反应了对系统阶跃判定的严格程度。对于阶跃发生时间的判定而言,θ越大表明阶跃越明显,因此对于起始阶跃时间而言,判别阈值并非越大越好。在一较佳地实施例中,将起始阶跃时间的判定θup设定为0.95,即通过SSCi′首次小于0.95时对应的时间判定A点(参见图16);对于阶跃发生后系统再次趋稳的时间判定而言,θ越大表明系统越稳定,但由于系统运行自身存在多重不确定性及非稳态情况,在一较佳地实施例中,将系统阶跃后趋稳的判定θlow设定为0.05,即通过SSCi′最后一次大于0.05时对应的时间判定M点(参见图16);但本申请不以此为限。
S103:根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差。
具体地,一实施例中,参见图5,所述参数辨识误差包括:上升参数辨识误差及调节参数辨识误差;所述根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差,包括:
S501:根据所述起始阶跃时间及稳态进入时间确定所述参数样本的阶跃量。
其中,参见图6,步骤S501包括:根据所述仿真样本的起始阶跃时间及所述仿真样本的稳态进入时间确定所述仿真样本的阶跃量(S601);根据所述实测样本的起始阶跃时间及所述实测样本的稳态进入时间确定所述实测样本的阶跃量(S602)。
S502:根据所述参数样本的阶跃量确定所述参数样本的上升时间偏差值及调节时间偏差值。
其中,参见图7,步骤S502包括:根据所述仿真样本的阶跃量确定所述仿真样本的第一上升时间偏差值及第一调节时间偏差值(S701);根据所述实测样本的阶跃量确定所述实测样本的第二上升时间偏差值及第二调节时间偏差值(S702);根据所述第一上升时间偏差值及所述第二上升时间偏差值确定所述上升时间偏差值(S703);根据所述第一调节时间偏差值及所述第二调节时间偏差值确定所述调节时间偏差值(S704)。
可以理解的是,上述步骤S501至步骤S502,需要针对仿真样本及实测样本分别执行,但执行方法相同。
参见图16,假设参数样本的起始阶跃点(对应于起始阶跃时间)为A(t0,U0)、最先达到稳态的点(对应于稳态进入时间)为M(t1,U1)。则阶跃量为ΔU=U1-U0。
曲线中被控量(本发明指高压缸调节汽门开度反馈)变化至90%阶跃量所对应的点为B(t2,U0+0.9ΔU),由此可得:
tup=t2-t0
通过编程,运用程序中的循环函数可以找到曲线中所有对应95%阶跃量和105%阶跃量的点C1(tc1,U0+0.95ΔU)、C2(tc2,U0+1.05ΔU)、C3(tc3,U0+0.95ΔU)、C4(tc4,U0+1.05ΔU)……Cn(tcn,U0+1.05ΔU)。
找到横坐标最大的点Cn(tcn,Un)
tcn=max(tc1,tc2,…)
由此可得:
ts=tcn-t0
分别对仿真样本及实测样本执行上述公式算法,可以分别得到仿真样本及实测样本对应的上升时间tup及调节时间ts。
设tup0、ts0分别为实测曲线的上升时间、实测曲线的调节时间;tup’、ts’分别为仿真曲线的上升时间、仿真曲线的调节时间。
则上升时间偏差值Δtup和调节时间偏差值Δts的计算公式:
Δtup=tup0-tup’
Δts=ts0-ts'
S503:比较预设的上升时间偏差允许值与所述上升时间偏差值,得到所述上升参数辨识误差;
S504:比较预设的调节时间偏差允许值与所述调节时间偏差值,得到所述调节参数辨识误差。
可以理解的是,行业标准DL/T 1235-2019《同步发电机原动机及其调节系统参数实测与建模导则》(以下简称:《导则》)规定了参数实测的试验内容包括静态试验和负载试验。
依据《导则》,汽轮机执行机构仿真与实测的偏差要求如表1-1所示。
表1-1汽轮机执行机构仿真与实测的偏差允许值
tup:阶跃试验中,从阶跃量加入开始到被控量变化至90%阶跃量所需时间,如图16所示。
ts:从起始时间开始,到被控量与最终稳态值之差的绝对值始终不超过5%阶跃量的最短时间,参见图16所示。
从上述描述可知,本申请提供的汽轮机调速系统执行机构参数辨识误差的分析方法及装置,能够基于改进的滑动窗口法对汽轮机调速系统执行机构参数辨识的误差进行分析与计算,避免了传统手动标点法中人为判断的影响,计算结果更加准确可靠。
为了更清楚地说明本申请提供的方法,下面举一实例。
其中,参见图17,实测曲线并非连续的曲线,而是采样周期为0.001s的二维数列,横坐标是时间,纵坐标是被控量(本发明为高压缸调节汽门开度反馈的阶跃量)。
以某300MW超临界、一次中间再热燃煤发电机组为例,对其执行机构(高调门GV1)进行参数辨识,得到实测数据与仿真数据如图17所示,然后分别运用传统方法和本发明算法进行误差分析计算,所得结果如表2-1所示:图17为某300MW机组GV1参数辨识结果
由表2-1可见,本发明所述算法误差计算结果与传统手动标点方法误差计算结果相近且满足《导则》要求,具有较好的实用性。
表2-1误差分析计算结果比较
综上所述,本申请至少存在如下有益效果:汽轮机及其调节系统参数辨识对于电网稳定性的分析具有重要的价值。汽轮机调速系统执行机构参数辨识是其中的重要组成部分,其辨识误差受到《导则》的严格限制。当前的误差算法受人为判断的影响较大且数据处理效率低,亟需改进。本发明基于改进滑动窗口法完成了汽轮机调速系统执行机构参数辨识的误差分析与计算。该方法避免传统手动标点法中人为判断的影响,使得计算结果更科学、可靠。并且本发明的计算过程迅速,计算结果准确,满足《导则》要求。
基于同一发明构思,本申请实施例还提供了一种汽轮机调速系统执行机构参数辨识误差的分析装置,可以用于实现上述实施例所描述的方法,如下面的实施例所述。由于汽轮机调速系统执行机构参数辨识误差的分析装置解决问题的原理与汽轮机调速系统执行机构参数辨识误差的分析方法相似,因此汽轮机调速系统执行机构参数辨识误差的分析装置的实施可以参见基于软件性能基准确定方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的系统较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
一实施例中,参见图8,为了能够根据预设的滑动时间窗口参数确定汽轮机调速系统执行机构的参数辨识误差,本申请提供一种汽轮机调速系统执行机构参数辨识误差的分析装置,包括:样本标准差确定单元801、阶跃稳态时间确定单元802及参数辨识误差确定单元803。
样本标准差确定单元801,用于根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差;
阶跃稳态时间确定单元802,用于根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间;
参数辨识误差确定单元803,用于根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差。
一实施例中,参见图9,所述滑动时间窗口参数包括:所述滑动时间窗口的宽度及数量;所述参数样本包括:仿真样本及实测样本;所述样本标准差包括:仿真样本标准差及实测样本标准差;所述样本标准差确定单元801,包括:仿真标准差确定模块901及实测标准差确定模块902。
仿真标准差确定模块901,用于根据所述滑动时间窗口的宽度及数量对所述仿真样本进行分析,得到各所述滑动时间窗口对应的仿真样本标准差;
实测标准差确定模块902,用于根据所述滑动时间窗口的宽度及数量对所述实测样本进行分析,得到各所述滑动时间窗口对应的实测样本标准差。
一实施例中,参见图10,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述阶跃稳态时间确定单元802,包括:仿真归一化标准差确定模块1001、仿真起始阶跃时间确定模块1002及仿真稳态进入时间确定模块1003。
仿真归一化标准差确定模块1001,用于对各所述滑动时间窗口对应的仿真样本标准差进行归一化,得到各所述滑动时间窗口对应的仿真样本归一化标准差;
仿真起始阶跃时间确定模块1002,用于按照所述仿真样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,首个大于所述阶跃判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的起始阶跃时间;
仿真稳态进入时间确定模块1003,用于按照所述仿真样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,首个小于所述稳态判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的稳态进入时间。
一实施例中,参见图11,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述阶跃稳态时间确定单元802,包括:实测归一化标准差确定模块1101、实测起始阶跃时间确定模块1102及实测稳态进入时间确定模块1103。
实测归一化标准差确定模块1101,用于对各所述滑动时间窗口对应的实测样本标准差进行归一化,得到各所述滑动时间窗口对应的实测样本归一化标准差;
实测起始阶跃时间确定模块1102,用于按照所述实测样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,首个小于所述阶跃判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的起始阶跃时间;
实测稳态进入时间确定模块1103,用于按照所述实测样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,最后一个大于所述稳态判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的稳态进入时间。
一实施例中,参见图12,所述参数辨识误差包括:上升参数辨识误差及调节参数辨识误差;所述参数辨识误差确定单元803,包括:
阶跃量确定模块1201,用于根据所述起始阶跃时间及稳态进入时间确定所述参数样本的阶跃量;
偏差值确定模块1202,用于根据所述参数样本的阶跃量确定所述参数样本的上升时间偏差值及调节时间偏差值;
上升参数辨识误差确定模块1203,用于比较预设的上升时间偏差允许值与所述上升时间偏差值,得到所述上升参数辨识误差;
调节参数辨识误差确定模块1204,用于比较预设的调节时间偏差允许值与所述调节时间偏差值,得到所述调节参数辨识误差。
一实施例中,参见图13,所述阶跃量确定模块1201,包括:
仿真阶跃量确定子模块1301,用于根据所述仿真样本的起始阶跃时间及所述仿真样本的稳态进入时间确定所述仿真样本的阶跃量;
实测阶跃量确定子模块1302,用于根据所述实测样本的起始阶跃时间及所述实测样本的稳态进入时间确定所述实测样本的阶跃量。
一实施例中,参见图14,所述偏差值确定模块1202,包括:
仿真偏差值确定子模块1401,用于根据所述仿真样本的阶跃量确定所述仿真样本的第一上升时间偏差值及第一调节时间偏差值;
实测偏差值确定子模块1402,用于根据所述实测样本的阶跃量确定所述实测样本的第二上升时间偏差值及第二调节时间偏差值;
上升时间偏差值确定子模块1403,用于根据所述第一上升时间偏差值及所述第二上升时间偏差值确定所述上升时间偏差值;
调节时间偏差值确定子模块1404,用于根据所述第一调节时间偏差值及所述第二调节时间偏差值确定所述调节时间偏差值。
从硬件层面来说,为了能够根据预设的滑动时间窗口参数确定汽轮机调速系统执行机构的参数辨识误差,本申请提供一种用于实现所述汽轮机调速系统执行机构参数辨识误差的分析方法中的全部或部分内容的电子设备的实施例,所述电子设备具体包含有如下内容:
处理器(Processor)、存储器(Memory)、通讯接口(Communications Interface)和总线;其中,所述处理器、存储器、通讯接口通过所述总线完成相互间的通讯;所述通讯接口用于实现所述汽轮机调速系统执行机构参数辨识误差的分析装置与核心业务系统、用户终端以及相关数据库等相关设备之间的信息传输;该逻辑控制器可以是台式计算机、平板电脑及移动终端等,本实施例不限于此。在本实施例中,该逻辑控制器可以参照实施例中的汽轮机调速系统执行机构参数辨识误差的分析方法的实施例,以及汽轮机调速系统执行机构参数辨识误差的分析装置的实施例进行实施,其内容被合并于此,重复之处不再赘述。
可以理解的是,所述用户终端可以包括智能手机、平板电子设备、网络机顶盒、便携式计算机、台式电脑、个人数字助理(PDA)、车载设备、智能穿戴设备等。其中,所述智能穿戴设备可以包括智能眼镜、智能手表、智能手环等。
在实际应用中,汽轮机调速系统执行机构参数辨识误差的分析方法的部分可以在如上述内容所述的电子设备侧执行,也可以所有的操作都在所述客户端设备中完成。具体可以根据所述客户端设备的处理能力,以及用户使用场景的限制等进行选择。本申请对此不作限定。若所有的操作都在所述客户端设备中完成,所述客户端设备还可以包括处理器。
上述的客户端设备可以具有通讯模块(即通讯单元),可以与远程的服务器进行通讯连接,实现与所述服务器的数据传输。所述服务器可以包括任务调度中心一侧的服务器,其他的实施场景中也可以包括中间平台的服务器,例如与任务调度中心服务器有通讯链接的第三方服务器平台的服务器。所述的服务器可以包括单台计算机设备,也可以包括多个服务器组成的服务器集群,或者分布式装置的服务器结构。
图15为本申请实施例的电子设备9600的系统构成的示意框图。如图15所示,该电子设备9600可以包括中央处理器9100和存储器9140;存储器9140耦合到中央处理器9100。值得注意的是,该图15是示例性的;还可以使用其他类型的结构,来补充或代替该结构,以实现电信功能或其他功能。
一实施例中,汽轮机调速系统执行机构参数辨识误差的分析方法功能可以被集成到中央处理器9100中。其中,中央处理器9100可以被配置为进行如下控制:
S101:根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差;
S102:根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间;
S103:根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差。
从上述描述可知,本申请提供的汽轮机调速系统执行机构参数辨识误差的分析方法及装置,能够基于改进的滑动窗口法对汽轮机调速系统执行机构参数辨识的误差进行分析与计算,避免了传统手动标点法中人为判断的影响,计算结果更加准确可靠。
在另一个实施方式中,汽轮机调速系统执行机构参数辨识误差的分析装置可以与中央处理器9100分开配置,例如可以将数据复合传输装置汽轮机调速系统执行机构参数辨识误差的分析装置配置为与中央处理器9100连接的芯片,通过中央处理器的控制来实现汽轮机调速系统执行机构参数辨识误差的分析方法的功能。
如图15所示,该电子设备9600还可以包括:通讯模块9110、输入单元9120、音频处理器9130、显示器9160、电源9170。值得注意的是,电子设备9600也并不是必须要包括图15中所示的所有部件;此外,电子设备9600还可以包括图15中没有示出的部件,可以参考现有技术。
如图15所示,中央处理器9100有时也称为控制器或操作控件,可以包括微处理器或其他处理器装置和/或逻辑装置,该中央处理器9100接收输入并控制电子设备9600的各个部件的操作。
其中,存储器9140,例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息,此外还可存储执行有关信息的程序。并且中央处理器9100可执行该存储器9140存储的该程序,以实现信息存储或处理等。
输入单元9120向中央处理器9100提供输入。该输入单元9120例如为按键或触摸输入装置。电源9170用于向电子设备9600提供电力。显示器9160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器,但并不限于此。
该存储器9140可以是固态存储器,例如,只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器,其即使在断电时也保存信息,可被选择性地擦除且设有更多数据,该存储器的示例有时被称为EPROM等。存储器9140还可以是某种其它类型的装置。存储器9140包括缓冲存储器9141(有时被称为缓冲器)。存储器9140可以包括应用/功能存储部9142,该应用/功能存储部9142用于存储应用程序和功能程序或用于通过中央处理器9100执行电子设备9600的操作的流程。
存储器9140还可以包括数据存储部9143,该数据存储部9143用于存储数据,例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器9140的驱动程序存储部9144可以包括电子设备的用于通讯功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通讯模块9110即为经由天线9111发送和接收信号的发送机/接收机9110。通讯模块(发送机/接收机)9110耦合到中央处理器9100,以提供输入信号和接收输出信号,这可以和常规移动通讯终端的情况相同。
基于不同的通讯技术,在同一电子设备中,可以设置有多个通讯模块9110,如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通讯模块(发送机/接收机)9110还经由音频处理器9130耦合到扬声器9131和麦克风9132,以经由扬声器9131提供音频输出,并接收来自麦克风9132的音频输入,从而实现通常的电信功能。音频处理器9130可以包括任何合适的缓冲器、解码器、放大器等。另外,音频处理器9130还耦合到中央处理器9100,从而使得可以通过麦克风9132能够在本机上录音,且使得可以通过扬声器9131来播放本机上存储的声音。
本申请的实施例还提供能够实现上述实施例中的执行主体为服务器或客户端的汽轮机调速系统执行机构参数辨识误差的分析方法中全部步骤的一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的执行主体为服务器或客户端的汽轮机调速系统执行机构参数辨识误差的分析方法的全部步骤,例如,所述处理器执行所述计算机程序时实现下述步骤:
S101:根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差;
S102:根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间;
S103:根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差。
从上述描述可知,本申请提供的汽轮机调速系统执行机构参数辨识误差的分析方法及装置,能够基于改进的滑动窗口法对汽轮机调速系统执行机构参数辨识的误差进行分析与计算,避免了传统手动标点法中人为判断的影响,计算结果更加准确可靠。
本领域内的技术人员应明白,本发明的实施例可提供为方法、装置、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(装置)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (17)
1.一种汽轮机调速系统执行机构参数辨识误差的分析方法,其特征在于,包括:
根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差;
根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间;
根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差。
2.根据权利要求1所述的汽轮机调速系统执行机构参数辨识误差的分析方法,其特征在于,所述滑动时间窗口参数包括:所述滑动时间窗口的宽度及数量;所述参数样本包括:仿真样本及实测样本;所述样本标准差包括:仿真样本标准差及实测样本标准差;所述根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差,包括:
根据所述滑动时间窗口的宽度及数量对所述仿真样本进行分析,得到各所述滑动时间窗口对应的仿真样本标准差;
根据所述滑动时间窗口的宽度及数量对所述实测样本进行分析,得到各所述滑动时间窗口对应的实测样本标准差。
3.根据权利要求2所述的汽轮机调速系统执行机构参数辨识误差的分析方法,其特征在于,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间,包括:
对各所述滑动时间窗口对应的仿真样本标准差进行归一化,得到各所述滑动时间窗口对应的仿真样本归一化标准差;
按照所述仿真样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,首个小于所述阶跃判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的起始阶跃时间;
按照所述仿真样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,最后一个大于所述稳态判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的稳态进入时间。
4.根据权利要求2所述的汽轮机调速系统执行机构参数辨识误差的分析方法,其特征在于,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间,包括:
对各所述滑动时间窗口对应的实测样本标准差进行归一化,得到各所述滑动时间窗口对应的实测样本归一化标准差;
按照所述实测样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,首个小于所述阶跃判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的起始阶跃时间;
按照所述实测样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,最后一个大于所述稳态判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的稳态进入时间。
5.根据权利要求3或4所述的汽轮机调速系统执行机构参数辨识误差的分析方法,其特征在于,所述参数辨识误差包括:上升参数辨识误差及调节参数辨识误差;所述根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差,包括:
根据所述起始阶跃时间及稳态进入时间确定所述参数样本的阶跃量;
根据所述参数样本的阶跃量确定所述参数样本的上升时间偏差值及调节时间偏差值;
比较预设的上升时间偏差允许值与所述上升时间偏差值,得到所述上升参数辨识误差;
比较预设的调节时间偏差允许值与所述调节时间偏差值,得到所述调节参数辨识误差。
6.根据权利要求5所述的汽轮机调速系统执行机构参数辨识误差的分析方法,其特征在于,所述根据所述起始阶跃时间及稳态进入时间确定所述参数样本的阶跃量,包括:
根据所述仿真样本的起始阶跃时间及所述仿真样本的稳态进入时间确定所述仿真样本的阶跃量;
根据所述实测样本的起始阶跃时间及所述实测样本的稳态进入时间确定所述实测样本的阶跃量。
7.根据权利要求5所述的汽轮机调速系统执行机构参数辨识误差的分析方法,其特征在于,所述根据所述参数样本的阶跃量确定所述参数样本的上升时间偏差值及调节时间偏差值,包括:
根据所述仿真样本的阶跃量确定所述仿真样本的第一上升时间偏差值及第一调节时间偏差值;
根据所述实测样本的阶跃量确定所述实测样本的第二上升时间偏差值及第二调节时间偏差值;
根据所述第一上升时间偏差值及所述第二上升时间偏差值确定所述上升时间偏差值;
根据所述第一调节时间偏差值及所述第二调节时间偏差值确定所述调节时间偏差值。
8.一种汽轮机调速系统执行机构参数辨识误差的分析装置,其特征在于,包括:
样本标准差确定单元,用于根据预设的滑动时间窗口参数分析汽轮机调速系统执行机构的参数样本,得到对应的样本标准差;
阶跃稳态时间确定单元,用于根据预设的样本阶跃阈值及所述样本标准差确定所述参数样本的起始阶跃时间及稳态进入时间;
参数辨识误差确定单元,用于根据所述起始阶跃时间及稳态进入时间确定所述汽轮机调速系统执行机构的参数辨识误差。
9.根据权利要求8所述的汽轮机调速系统执行机构参数辨识误差的分析装置,其特征在于,所述滑动时间窗口参数包括:所述滑动时间窗口的宽度及数量;所述参数样本包括:仿真样本及实测样本;所述样本标准差包括:仿真样本标准差及实测样本标准差;所述样本标准差确定单元,包括:
仿真标准差确定模块,用于根据所述滑动时间窗口的宽度及数量对所述仿真样本进行分析,得到各所述滑动时间窗口对应的仿真样本标准差;
实测标准差确定模块,用于根据所述滑动时间窗口的宽度及数量对所述实测样本进行分析,得到各所述滑动时间窗口对应的实测样本标准差。
10.根据权利要求9所述的汽轮机调速系统执行机构参数辨识误差的分析装置,其特征在于,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述阶跃稳态时间确定单元,包括:
仿真归一化标准差确定模块,用于对各所述滑动时间窗口对应的仿真样本标准差进行归一化,得到各所述滑动时间窗口对应的仿真样本归一化标准差;
仿真起始阶跃时间确定模块,用于按照所述仿真样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,首个小于所述阶跃判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的起始阶跃时间;
仿真稳态进入时间确定模块,用于按照所述仿真样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的仿真样本归一化标准差,并确定在所述仿真样本的采样时序上,最后一个大于所述稳态判断阈值的仿真样本归一化标准差所对应的滑动时间窗口,得到所述仿真样本的稳态进入时间。
11.根据权利要求9所述的汽轮机调速系统执行机构参数辨识误差的分析装置,其特征在于,所述样本阶跃阈值包括:阶跃判断阈值及稳态判断阈值;所述阶跃稳态时间确定单元,包括:
实测归一化标准差确定模块,用于对各所述滑动时间窗口对应的实测样本标准差进行归一化,得到各所述滑动时间窗口对应的实测样本归一化标准差;
实测起始阶跃时间确定模块,用于按照所述实测样本的采样时序,比较所述阶跃判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,首个小于所述阶跃判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的起始阶跃时间;
实测稳态进入时间确定模块,用于按照所述实测样本的采样时序,比较所述稳态判断阈值与各所述滑动时间窗口对应的实测样本归一化标准差,并确定在所述实测样本的采样时序上,最后一个大于所述稳态判断阈值的实测样本归一化标准差所对应的滑动时间窗口,以确定所述实测样本的稳态进入时间。
12.根据权利要求10或11所述的汽轮机调速系统执行机构参数辨识误差的分析装置,其特征在于,所述参数辨识误差包括:上升参数辨识误差及调节参数辨识误差;所述参数辨识误差确定单元,包括:
阶跃量确定模块,用于根据所述起始阶跃时间及稳态进入时间确定所述参数样本的阶跃量;
偏差值确定模块,用于根据所述参数样本的阶跃量确定所述参数样本的上升时间偏差值及调节时间偏差值;
上升参数辨识误差确定模块,用于比较预设的上升时间偏差允许值与所述上升时间偏差值,得到所述上升参数辨识误差;
调节参数辨识误差确定模块,用于比较预设的调节时间偏差允许值与所述调节时间偏差值,得到所述调节参数辨识误差。
13.根据权利要求12所述的汽轮机调速系统执行机构参数辨识误差的分析装置,其特征在于,所述阶跃量确定模块,包括:
仿真阶跃量确定子模块,用于根据所述仿真样本的起始阶跃时间及所述仿真样本的稳态进入时间确定所述仿真样本的阶跃量;
实测阶跃量确定子模块,用于根据所述实测样本的起始阶跃时间及所述实测样本的稳态进入时间确定所述实测样本的阶跃量。
14.根据权利要求13所述的汽轮机调速系统执行机构参数辨识误差的分析装置,其特征在于,所述偏差值确定模块,包括:
仿真偏差值确定子模块,用于根据所述仿真样本的阶跃量确定所述仿真样本的第一上升时间偏差值及第一调节时间偏差值;
实测偏差值确定子模块,用于根据所述实测样本的阶跃量确定所述实测样本的第二上升时间偏差值及第二调节时间偏差值;
上升时间偏差值确定子模块,用于根据所述第一上升时间偏差值及所述第二上升时间偏差值确定所述上升时间偏差值;
调节时间偏差值确定子模块,用于根据所述第一调节时间偏差值及所述第二调节时间偏差值确定所述调节时间偏差值。
15.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求1至7任一项所述的汽轮机调速系统执行机构参数辨识误差的分析方法的步骤。
16.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至7任一项所述的汽轮机调速系统执行机构参数辨识误差的分析方法的步骤。
17.一种计算机程序产品,包括计算机程序/指令,其特征在于,该计算机程序/指令被处理器执行时实现权利要求1至7任一项所述的汽轮机调速系统执行机构参数辨识误差的分析方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210514864.XA CN114718669B (zh) | 2022-05-11 | 2022-05-11 | 汽轮机调速系统执行机构参数辨识误差的分析方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210514864.XA CN114718669B (zh) | 2022-05-11 | 2022-05-11 | 汽轮机调速系统执行机构参数辨识误差的分析方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114718669A true CN114718669A (zh) | 2022-07-08 |
CN114718669B CN114718669B (zh) | 2024-01-26 |
Family
ID=82230791
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210514864.XA Active CN114718669B (zh) | 2022-05-11 | 2022-05-11 | 汽轮机调速系统执行机构参数辨识误差的分析方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114718669B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN85205476U (zh) * | 1985-12-16 | 1986-11-26 | 华中工学院 | 模型辨识仪 |
US6549861B1 (en) * | 2000-08-10 | 2003-04-15 | Euro-Celtique, S.A. | Automated system and method for spectroscopic analysis |
CN1487268A (zh) * | 2003-08-12 | 2004-04-07 | 上海交通大学 | 基于阶跃响应测试的多变量系统结构化闭环辨识方法 |
US20050200347A1 (en) * | 2004-03-15 | 2005-09-15 | Maria Agoston | Real-time time drift adjustment for a TDR step stimulus |
CN101924533A (zh) * | 2010-07-19 | 2010-12-22 | 浙江工业大学 | 基于fir模型辨识的多变量时滞参数估计方法 |
CN103577647A (zh) * | 2013-11-11 | 2014-02-12 | 国家电网公司 | 一种汽轮机及其调速系统模型的评估分析方法及装置 |
CN105068423A (zh) * | 2015-06-03 | 2015-11-18 | 贵州电力试验研究院 | 一键式实现汽轮机及其调速系统参数智能辨识的方法 |
CN105912779A (zh) * | 2016-04-12 | 2016-08-31 | 国网上海市电力公司 | 汽轮机调速系统模型参数辨识方法 |
US20160252562A1 (en) * | 2013-10-21 | 2016-09-01 | Washington State University | Performance analysis of power grid monitors |
CN108564205A (zh) * | 2018-03-27 | 2018-09-21 | 昆明理工大学 | 一种基于实测数据的负荷模型及参数辨识优化方法 |
US10690722B1 (en) * | 2019-02-08 | 2020-06-23 | Real Intent, Inc. | Methods and systems for efficient identification of glitch failures in integrated circuits |
CN112286043A (zh) * | 2020-10-13 | 2021-01-29 | 国网浙江省电力有限公司电力科学研究院 | 基于被控对象阶跃响应特性数据的pid参数整定方法 |
-
2022
- 2022-05-11 CN CN202210514864.XA patent/CN114718669B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN85205476U (zh) * | 1985-12-16 | 1986-11-26 | 华中工学院 | 模型辨识仪 |
US6549861B1 (en) * | 2000-08-10 | 2003-04-15 | Euro-Celtique, S.A. | Automated system and method for spectroscopic analysis |
CN1487268A (zh) * | 2003-08-12 | 2004-04-07 | 上海交通大学 | 基于阶跃响应测试的多变量系统结构化闭环辨识方法 |
US20050200347A1 (en) * | 2004-03-15 | 2005-09-15 | Maria Agoston | Real-time time drift adjustment for a TDR step stimulus |
CN101924533A (zh) * | 2010-07-19 | 2010-12-22 | 浙江工业大学 | 基于fir模型辨识的多变量时滞参数估计方法 |
US20160252562A1 (en) * | 2013-10-21 | 2016-09-01 | Washington State University | Performance analysis of power grid monitors |
CN103577647A (zh) * | 2013-11-11 | 2014-02-12 | 国家电网公司 | 一种汽轮机及其调速系统模型的评估分析方法及装置 |
CN105068423A (zh) * | 2015-06-03 | 2015-11-18 | 贵州电力试验研究院 | 一键式实现汽轮机及其调速系统参数智能辨识的方法 |
CN105912779A (zh) * | 2016-04-12 | 2016-08-31 | 国网上海市电力公司 | 汽轮机调速系统模型参数辨识方法 |
CN108564205A (zh) * | 2018-03-27 | 2018-09-21 | 昆明理工大学 | 一种基于实测数据的负荷模型及参数辨识优化方法 |
US10690722B1 (en) * | 2019-02-08 | 2020-06-23 | Real Intent, Inc. | Methods and systems for efficient identification of glitch failures in integrated circuits |
CN112286043A (zh) * | 2020-10-13 | 2021-01-29 | 国网浙江省电力有限公司电力科学研究院 | 基于被控对象阶跃响应特性数据的pid参数整定方法 |
Non-Patent Citations (2)
Title |
---|
孙海顺等: "辨识理论在发电机励磁系统参数测试中的应用", 《高电压技术》, vol. 32, no. 10, pages 102 - 105 * |
王修中等: "二阶加滞后连续模型的直接辨识", 《自动化学报》, vol. 27, no. 5, pages 728 - 731 * |
Also Published As
Publication number | Publication date |
---|---|
CN114718669B (zh) | 2024-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108833458B (zh) | 一种应用推荐方法、装置、介质及设备 | |
CN108388509B (zh) | 一种软件测试方法、计算机可读存储介质及终端设备 | |
CN110291540A (zh) | 批再归一化层 | |
KR20210017342A (ko) | 과거 데이터에 기초한 시계열 예측 방법 및 장치 | |
CN109460943A (zh) | 关键绩效指标管理方法、设备、存储介质及装置 | |
CN113537510A (zh) | 基于不均衡数据集的机器学习模型数据处理方法及装置 | |
CN115392037A (zh) | 设备故障预测方法、装置、设备及存储介质 | |
CN105022686A (zh) | 安卓平台上进程的cpu占用率确定方法和装置 | |
CN111027643B (zh) | 深度神经网络模型的训练方法、人机交互识别方法、装置、电子设备及存储介质 | |
CN116684330A (zh) | 基于人工智能的流量预测方法、装置、设备及存储介质 | |
US8327196B2 (en) | Identifying an optimized test bit pattern for analyzing electrical communications channel topologies | |
CN111177500A (zh) | 数据对象分类方法、装置、计算机设备和存储介质 | |
CN113077815B (zh) | 一种音频评估方法及组件 | |
US20190042681A1 (en) | Method and apparatus for selecting integrated circuit device neural network modeling sample | |
CN114718669B (zh) | 汽轮机调速系统执行机构参数辨识误差的分析方法及装置 | |
CN111882423A (zh) | 存款利率信息推送方法及装置 | |
US8700686B1 (en) | Robust estimation of time varying parameters | |
CN113190429B (zh) | 服务器性能预测方法、装置及终端设备 | |
CN111310794B (zh) | 目标对象的分类方法、装置和电子设备 | |
US20220215271A1 (en) | Detection device, detection method and detection program | |
Pop et al. | The estimation of the lifetime variation for power devices | |
CN115374555A (zh) | 电功率阶跃响应仿真辨识误差的分析方法及装置 | |
CN113435058B (zh) | 配网自愈测试模型的数据降维方法、系统、终端及介质 | |
CN116302661B (zh) | 一种异常预测方法、装置、电子设备及存储介质 | |
CN114676656B (zh) | 多响应cfd模型的一致性度量方法、装置、设备及存储介质 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |