CN105275730B - 一种水轮机调速器微分环节系数及时间常数辨识方法 - Google Patents

一种水轮机调速器微分环节系数及时间常数辨识方法 Download PDF

Info

Publication number
CN105275730B
CN105275730B CN201510670936.XA CN201510670936A CN105275730B CN 105275730 B CN105275730 B CN 105275730B CN 201510670936 A CN201510670936 A CN 201510670936A CN 105275730 B CN105275730 B CN 105275730B
Authority
CN
China
Prior art keywords
mrow
msub
gate control
control signal
differentiation element
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.)
Active
Application number
CN201510670936.XA
Other languages
English (en)
Other versions
CN105275730A (zh
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.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
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 State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd, State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201510670936.XA priority Critical patent/CN105275730B/zh
Publication of CN105275730A publication Critical patent/CN105275730A/zh
Application granted granted Critical
Publication of CN105275730B publication Critical patent/CN105275730B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/20Hydro energy

Landscapes

  • Complex Calculations (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了一种水轮机调速器微分环节系数及时间常数辨识方法,包括:实际微分环节数据获取、数据预处理、构建实际的导叶控制信号的响应函数及实际微分环节参数辨识。所述实际微分环节数据获取使用信号发生器或调速器综合测试仪模拟机组频率扰动信号,并使用调速器综合测试仪录取导叶控制信号及模拟机频信号;对工程实际中激励信号与理想激励信号分析,确定出理论分析激励信号传递函数表达式;所述实际微分环节参数辨识基于实测数据实现对实际微分环节参数辨识。本发明通过实现了调速器实际微分环节的参数辨识,实现了多次频率扰动的数据自动分段,适用于现场实测数据及仿真数据,具有较强的抗干扰性。

Description

一种水轮机调速器微分环节系数及时间常数辨识方法
技术领域
本发明涉及一种辨识方法,特别涉及一种水轮机调速器微分环节系数及时间常数辨识方法。
背景技术
获取准确的水轮机调节系统(包括:调速器、液压执行机构、水轮机及引水系统)参数对于保证电力系统仿真的准确性,预测电力系统的动态响应及电力系统故障重现具有重要意义。目前电力系统中励磁调节系统、发电机等环节的参数辨识较为成熟,但水轮机调节调节系统参数辨识相对滞后,严重制约了电力系统仿真的精度,故水轮机调节系统参数辨识理论与实践亟需完善。在开展调速器参数辨识中,需要辨识调速器实际微分环节的参数。通常激励信号采用阶跃扰动,系统响应通常采用导叶控制信号yPID。然而工程实际中由于实际阶跃扰动信号并非理想的阶跃扰动,扰动过程持续时间较长,导致阶跃扰动信号中难以激发实际微分环节所有模态;另一方面,部分调速器仿真测试仪使用阶跃扰动,部分调速器使用斜坡信号扰动;此外,工程实测信号含有大量的干扰,信号具有“毛刺”特征,影响了辨识结果的准确性,限制了一些辨识方法的应用。
因此,急需提出一种适合工程实测数据的水轮机调速器微分环节系数及时间常数辨识方法,要求该方法能够适用于导叶接力器行程信号,适用于不同的激励信号,且对工程实测信号具有一定的抗干扰性。
发明内容
本发明的目的在于提供一种水轮机调速器微分环节系数及时间常数辨识方法,提高微分环节参数辨识精度、确保调速器建模的准确度。
一种水轮机调速器微分环节系数及时间常数辨识方法,包括以下步骤:
步骤1:获取水轮机实际微分环节数据;
选择斜坡函数作为水轮机调速器的激励信号,通过信号发生器或调速器综合测试仪发出模拟机组频率扰动,采集导叶控制信号YPID和模拟机组频率信号;
步骤2:对步骤1获取的导叶控制信号进行去均值处理;
步骤3:构建实际的导叶控制信号的响应函数;
其中,Kd为待辨识的微分系数,T1v为待辨识的微分环节时间常数,t为时间,k为斜坡信号斜率,yPID为去均值后的实测导叶控制信号,Tmax为实测的导叶控制信号yPID取最大值yPIDmax所对应的时刻;
步骤4:依据t=Tmax+T1v时,yPID从最大值衰减1个T1v时间后,yPID值将下降36.8%,并结合以下公式辨识微分环节时间常数T1v
步骤5:利用当t=Tmax时,获得待辨识的微分系数Kd
所述步骤2中对导叶控制信号进行去均值处理公式如下:
式中,t为采样时刻,YPID(i·Δt)为t=i·Δt时刻去均值前的导叶控制信号,yPID(i·Δt)为t=i·Δt时刻去均值后的导叶控制信号,Δt为时域采样间隔,N为阶跃扰动前稳态导叶控制信号数值点数,N=S/Δt,S表示导叶控制信号扰动前的稳态数据时长。
对导叶控制信号进行去均值处理前,先进行数据分段处理,采用机组频率微分法确定扰动发生时刻,自动获得导叶控制信号扰动前稳态数据时长;
所述机组频率微分法是指对机组频率信号进行差分操作,根据突变值可以确定扰动的发生时刻。
有益效果
本发明提供了一种水轮机调速器微分环节系数及时间常数辨识方法,包括:实际微分环节数据获取、数据预处理、构建实际的导叶控制信号的响应函数及实际微分环节参数辨识。所述实际微分环节数据获取使用信号发生器或调速器综合测试仪模拟机组频率扰动信号,并使用调速器综合测试仪录取导叶控制信号及模拟机频信号;所述数据预处理包括扰动数据分段、数据去均值处理;对工程实际中激励信号与理想激励信号分析,确定出理论分析激励信号传递函数表达式;通过理论推导建立了微分环节参数与微分环节响应输出特征值间关系;所述实际微分环节参数辨识基于实测数据实现对实际微分环节参数辨识。本发明通过实现了调速器实际微分环节的参数辨识,实现了多次频率扰动的数据自动分段,提出的实际微分环节时间参数T1v和Kd系数的辨识公式,适用于现场实测数据及仿真数据,具有较强的抗干扰性。
附图说明
图1为调速器实际微分环节参数测试流程图;
图2为实际微分环节数据获取现场试验接线示意图;
图3为机组频率微分法示意图,其中,(a)为机组频率实测波形,(b)机组频率微分波形;
图4为理想阶跃信号与实际阶跃信号对比图其中,(a)为理想阶跃扰动,(b)为实际阶跃扰动信号放大图,(c)为理想上升斜坡激励信号,(d)为理想下降斜坡激励信号;
图5为实际阶跃信号与实际导叶PID响应图,其中,(a)为实际频率阶跃激励信号,(b)为实际导叶PID响应曲线。
具体实施方式
下面将结合附图和实施例对本发明做进一步的说明。
如图1所示,一种水轮机调速器微分环节系数及时间常数辨识方法,包括以下步骤:
步骤1:获取水轮机实际微分环节数据;
选择斜坡函数作为水轮机调速器的激励信号,通过信号发生器或调速器综合测试仪发出模拟机组频率扰动,采集导叶控制信号YPID和模拟机组频率信号;
实际微分环节数据获取包括现场试验接线、模拟机组频率扰动。如图2所示,现场试验接线主要为调速器电器柜引取导叶控制信号YPID、调速器综合测试仪或信号发生器发送给调速器电器柜的模拟机组频率信号。模拟机组频率扰动前将调速器人工频率死区设置为零、比例项系数、积分项系数、永态差值系数设置为零,断路器信号短接,操作调速器综合测试仪或信号发生器按照50.0Hz—>50.1Hz—>50.0Hz—>49.9Hz—>50.0Hz、50.0Hz—>50.2Hz—>50.0Hz—>49.8Hz—>50.0Hz、50.0Hz—>50.3Hz—>50.0Hz—>49.7Hz—>50.0Hz的顺序进行频率扰动,使用调速器综合测试仪录取导叶控制信号YPID、模拟机组频率信号。
步骤2:对步骤1获取的导叶控制信号进行去均值处理;
对导叶控制信号进行去均值处理前,先进行数据分段处理,采用机组频率微分法确定扰动发生时刻,自动获得导叶控制信号扰动前稳态数据时长;
所述机组频率微分法是指对机组频率信号进行差分操作,根据突变值可以确定扰动的发生时刻。
所述数据分段处理的关键在于找到频率扰动的发生时刻,由于频率扰动为近似阶跃扰动,对机组频率信号进行差分操作,可以发现在扰动时刻机组频率信号发生显著突变,其示意图见图3。根据突变值可以确定扰动的发生时刻,根据扰动时刻可以选择指定的扰动。
实测的导叶控制信号往往在50%开度附近对调速器施加激励信号,因此,在辨识微分环节前需要去除导叶控制信号稳态均值。
所述步骤2中对导叶控制信号进行去均值处理公式如下:
式中,t为采样时刻,YPID(i·Δt)为t=i·Δt时刻去均值前的导叶控制信号,yPID(i·Δt)为t=i·Δt时刻去均值后的导叶控制信号,Δt为时域采样间隔,N为阶跃扰动前稳态导叶控制信号数值点数,N=S/Δt,S表示导叶控制信号扰动前的稳态数据时长。
如采样间隔为0.01s,实测数据稳态较长,可取N为500,即取导叶控制发生扰动前的5s内的导叶控制信号平均值为导叶控制信号的稳态值。
步骤3:构建实际的导叶控制信号的响应函数;
其中,Kd为待辨识的微分系数,T1v为待辨识的微分环节时间常数,t为时间,k为斜坡信号斜率,yPID为去均值后的实测导叶控制信号,Tmax为实测的导叶控制信号yPID取最大值yPIDmax所对应的时刻;
所述激励信号分析主要涉及实际阶跃扰动信号上升时间、扰动幅度分析。图4(a)给出了理想阶跃扰动激励信号,但是实际中的测试仪器所发出的阶跃扰动均存在一定的时间,表现出图4(b)所示的斜坡特征,通常斜坡时间较短,但要求实际微分环节参数辨识应按照斜坡函数来分析。
水轮机调速器常采用PID控制器,微分控制器常采用实际微分环节,表示为:
式中,Kd为待辨识的微分系数,T1v为待辨识的微分环节时间常数,s为拉普拉斯算子yPID为去均值后的导叶控制信号,x为机组频率信号
当输入激励信号为上升斜率为k的斜坡信号时(见图4(c)),实际微分环节输出响应的传递函数表示为:
式中,Kd为待辨识的微分系数,T1v为待辨识的微分环节时间常数,s为拉普拉斯算子,yPID为去均值后的导叶控制信号,x为机组频率信号
使用拉普拉斯逆变换可将式(3)从s域转换为时域。
式中,Kd为待辨识的微分系数,T1v为待辨识的微分环节时间常数,t为时间,k为斜坡信号斜率,yPID为去均值后的导叶控制信号。
当输入激励信号为下降斜率为-k的斜坡信号时(见图4(d)),实际微分环节输出响应的传递函数表示为:
式中,Kd为待辨识的微分系数,T1v为待辨识的微分环节时间常数,s为拉普拉斯算子,k为斜坡信号斜率,yPID为去均值后的导叶控制信号。
使用拉普拉斯逆变换可将式(5)从s域转换为时域。
式中,Kd为待辨识的微分系数,T1v为待辨识的微分环节时间常数,t为时间,k为斜坡信号斜率,yPID为去均值后的导叶控制信号。
式(4)、(6)给出的yPID响应仅适用于斜坡激励图4(c)-图4(d),对于图4(b)所示的斜坡激励,可以根据图4(c)和图4(d)叠加而得,由于实际微分环节为线性环节,故其输出为斜坡响应激励的分段叠加。
步骤4:依据t=Tmax+T1v时,yPID从最大值衰减1个T1v时间后,yPID值将下降36.8%,并结合以下公式辨识微分环节时间常数T1v
步骤5:利用当t=Tmax时,获得待辨识的微分系数Kd
为了验证该方法的辨识效果,设置调速器Kd=1.0,T1v=0.4s,对调速器施加机组频率扰动,实测数据见图5。从图5(b)中可以看出,扰动前导叶控制信号为50.24%,扰动时刻为20.65s,导叶控制信号YPID最大值为51.19%,最大值时刻为20.76s,故Tmax=20.76s-20.65s=0.11s,yPIDmax=51.19%-50.24%=0.95%;最大值衰减36.8%为0.95%*0.368=0.3496%,故衰减36.8%后的YPID值=50.24%+0.3496%=50.58%,对应的时刻为21.17s。根据最大值时刻和衰减36.8%后的时刻计算出T1v=21.17s-20.76s=0.4s。计算得到Kd如下:

Claims (3)

1.一种水轮机调速器微分环节系数及时间常数辨识方法,其特征在于,包括以下步骤:
步骤1:获取水轮机实际微分环节数据;
选择斜坡函数作为水轮机调速器的激励信号,通过信号发生器或调速器综合测试仪发出模拟机组频率扰动,采集导叶控制信号YPID和模拟机组频率信号;
步骤2:对步骤1获取的导叶控制信号进行去均值处理;
步骤3:构建实际的导叶控制信号的响应函数;
<mrow> <msub> <mi>y</mi> <mrow> <mi>P</mi> <mi>I</mi> <mi>D</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>kK</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mi>t</mi> <mrow> <mi>T</mi> <mn>1</mn> <mi>v</mi> </mrow> </mfrac> </mrow> </msup> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>t</mi> <mo>&amp;le;</mo> <msub> <mi>T</mi> <mi>max</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>K</mi> <mi>d</mi> </msub> <mi>k</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mi>T</mi> <mi>max</mi> </mrow> <mrow> <mi>T</mi> <mn>1</mn> <mi>v</mi> </mrow> </mfrac> </mrow> </msup> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mi>t</mi> <mo>-</mo> <mi>T</mi> <mi>max</mi> </mrow> <msub> <mi>T</mi> <mrow> <mn>1</mn> <mi>v</mi> </mrow> </msub> </mfrac> </mrow> </msup> </mrow> </mtd> <mtd> <mrow> <mi>t</mi> <mo>&gt;</mo> <msub> <mi>T</mi> <mi>max</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,Kd为待辨识的微分系数,T1v为待辨识的微分环节时间常数,t为时间,k为斜坡信号斜率,yPID为去均值后的实测导叶控制信号,Tmax为实测的导叶控制信号yPID取最大值yPIDmax所对应的时刻;
步骤4:依据t=Tmax+T1v时,yPID从最大值衰减1个T1v时间后,yPID值将下降36.8%,并结合以下公式辨识微分环节时间常数T1v
<mrow> <msub> <mi>y</mi> <mrow> <mi>P</mi> <mi>I</mi> <mi>D</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>kK</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <msub> <mi>T</mi> <mrow> <mn>1</mn> <mi>v</mi> </mrow> </msub> </mfrac> </mrow> </msup> <mo>)</mo> </mrow> </mrow>
步骤5:利用当t=Tmax时,获得待辨识的微分系数Kd
2.根据权利要求1所述的方法,其特征在于,所述步骤2中对导叶控制信号进行去均值处理公式如下:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mi>P</mi> <mi>I</mi> <mi>D</mi> </mrow> </msub> <mo>(</mo> <mi>i</mi> <mo>&amp;CenterDot;</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>)</mo> <mo>=</mo> <msub> <mi>Y</mi> <mrow> <mi>P</mi> <mi>I</mi> <mi>D</mi> </mrow> </msub> <mo>(</mo> <mi>i</mi> <mo>&amp;CenterDot;</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>)</mo> <mo>-</mo> <msub> <mi>y</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>y</mi> <mn>0</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>y</mi> <mrow> <mi>P</mi> <mi>I</mi> <mi>D</mi> <mi>o</mi> <mi>r</mi> <mi>g</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>&amp;CenterDot;</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
式中,t为采样时刻,YPID(i·Δt)为t=i·Δt时刻去均值前的导叶控制信号,yPID(i·Δt)为t=i·Δt时刻去均值后的导叶控制信号,Δt为时域采样间隔,N为阶跃扰动前稳态导叶控制信号数值点数,N=S/Δt,S表示导叶控制信号扰动前的稳态数据时长。
3.根据权利要求2所述的方法,其特征在于,对导叶控制信号进行去均值处理前,先进行数据分段处理,采用机组频率微分法确定扰动发生时刻,自动获得导叶控制信号扰动前稳态数据时长;
所述机组频率微分法是指对机组频率信号进行差分操作,根据突变值可以确定扰动的发生时刻。
CN201510670936.XA 2015-10-15 2015-10-15 一种水轮机调速器微分环节系数及时间常数辨识方法 Active CN105275730B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510670936.XA CN105275730B (zh) 2015-10-15 2015-10-15 一种水轮机调速器微分环节系数及时间常数辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510670936.XA CN105275730B (zh) 2015-10-15 2015-10-15 一种水轮机调速器微分环节系数及时间常数辨识方法

Publications (2)

Publication Number Publication Date
CN105275730A CN105275730A (zh) 2016-01-27
CN105275730B true CN105275730B (zh) 2017-09-12

Family

ID=55145358

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510670936.XA Active CN105275730B (zh) 2015-10-15 2015-10-15 一种水轮机调速器微分环节系数及时间常数辨识方法

Country Status (1)

Country Link
CN (1) CN105275730B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107834610B (zh) * 2017-11-29 2019-08-20 西南交通大学 一种考虑水轮机水锤效应的电网频率动态分析方法
CN111130414B (zh) * 2020-01-03 2021-09-14 沈机(上海)智能系统研发设计有限公司 电机平均电流平滑处理方法及系统、电机电流采样设备
CN113464354B (zh) * 2021-08-13 2023-06-23 西华大学 一种应用于具有长有压引水道水电站的水轮机控制方法
CN114856901A (zh) * 2021-12-22 2022-08-05 国家能源集团新疆开都河流域水电开发有限公司 一种水轮机调速器导叶开度调节方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101382967A (zh) * 2008-09-02 2009-03-11 陕西电力科学研究院 水轮机调速器主接力器的仿真方法
CN101404043A (zh) * 2008-09-10 2009-04-08 陕西电力科学研究院 水轮机及其有压过水系统仿真方法
CN103577647A (zh) * 2013-11-11 2014-02-12 国家电网公司 一种汽轮机及其调速系统模型的评估分析方法及装置
CN104063584A (zh) * 2014-06-10 2014-09-24 国家电网公司 一种汽轮机调速系统控制参数整定方法
CN104573197A (zh) * 2014-12-22 2015-04-29 国家电网公司 一种水轮机模型参数的测试方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101382967A (zh) * 2008-09-02 2009-03-11 陕西电力科学研究院 水轮机调速器主接力器的仿真方法
CN101404043A (zh) * 2008-09-10 2009-04-08 陕西电力科学研究院 水轮机及其有压过水系统仿真方法
CN103577647A (zh) * 2013-11-11 2014-02-12 国家电网公司 一种汽轮机及其调速系统模型的评估分析方法及装置
CN104063584A (zh) * 2014-06-10 2014-09-24 国家电网公司 一种汽轮机调速系统控制参数整定方法
CN104573197A (zh) * 2014-12-22 2015-04-29 国家电网公司 一种水轮机模型参数的测试方法及系统

Also Published As

Publication number Publication date
CN105275730A (zh) 2016-01-27

Similar Documents

Publication Publication Date Title
CN105275730B (zh) 一种水轮机调速器微分环节系数及时间常数辨识方法
CN107590317B (zh) 一种计及模型参数不确定性的发电机动态估计方法
CN106707760B (zh) 一种用于压电驱动器动态迟滞补偿的非线性逆控制方法
Crow et al. The matrix pencil for power system modal extraction
CN102156045B (zh) 水轮机调速器主接力器反应时间常数测试方法
CN108896857B (zh) 一种基于深度学习的变压器复杂工况识别方法
CN105353789B (zh) 连续振动信号时间历程复现控制方法
CN107066712B (zh) 基于导叶开度—功率特性的水轮机模型参数辨识方法及装置
CN111507637B (zh) 电力系统稳定计算用水轮机及引水系统参数建模与测试方法
CN102176653A (zh) 带指数渐消因子的卡尔曼滤波器的感应电机转速观测方法
CN104500150B (zh) 一种汽轮机调速器执行机构分段线性模型及参数获取方法
CN104391497A (zh) 一种关节伺服系统参数辨识和控制参数在线优化方法
CN103942355B (zh) 一种模拟水轮发电机组电液随动系统的建模方法
CN110098610A (zh) 故障扰动下电力系统振荡主导模式的实时辨识方法及系统
CN108152629B (zh) 一种用于火电厂绘制电压耐受力曲线的方法
CN103743955B (zh) 一种电机电感的测量方法
CN109190302A (zh) 一种接触器电磁系统机电关系的神经网络分段拟合方法及系统
CN104239614A (zh) 压气机气动失稳信号的模拟方法
CN107167733B (zh) 一种同步发电机励磁系统性能的基础评估数据的获取方法
CN115576189A (zh) 基于自适应齐次微分器的进气环境模拟系统pid控制方法
CN105391365B (zh) 一种同时检测转子电阻与励磁电感的静态辨识方法
CN108362940A (zh) 一种半周波畸变信号的动态谐波频率提取方法
Arastou et al. Dynamic state and parameter estimation of the improved Heffron-Phillips model using a fast UKF-based algorithm and a novel rotor angle measurement approach
CN107563047A (zh) 一种电力系统并列虚拟实验实现方法和装置
CN107342599B (zh) 一种机电暂态仿真中控制系统稳定性自动诊断和参数调整方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant