CN111880091A - 一种双轴励磁同步调相机励磁系统参数辨识方法 - Google Patents

一种双轴励磁同步调相机励磁系统参数辨识方法 Download PDF

Info

Publication number
CN111880091A
CN111880091A CN202010003559.5A CN202010003559A CN111880091A CN 111880091 A CN111880091 A CN 111880091A CN 202010003559 A CN202010003559 A CN 202010003559A CN 111880091 A CN111880091 A CN 111880091A
Authority
CN
China
Prior art keywords
active power
error
reactive power
true
excitation
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.)
Withdrawn
Application number
CN202010003559.5A
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.)
Harbin University of Science and Technology
Original Assignee
Harbin 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 Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN202010003559.5A priority Critical patent/CN111880091A/zh
Publication of CN111880091A publication Critical patent/CN111880091A/zh
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/34Testing dynamo-electric machines

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种双轴励磁同步调相机励磁系统参数辨识方法,用于输出双轴励磁同步调相机励磁系统参数的θ(kP,kI,kD)的最优解,包括步骤:构建双轴励磁系统仿真模型M1;确定等价标准J1和J2;对双轴励磁同步调相机励磁系统以信号X1和X2进行扰动试验,获取各真实响应指标Y2i;还包括步骤:对应参数θ(kP,kI,kD)的值进行选取,将参数θ(kP,kI,kD)的值代入M1,并对其以信号X1和X2进行扰动试验,获取各仿真响应指标Y1i;计算Y2i与Y1i之间的误差ei,确定ei对应的权重Wi,其中∑Wi=1;若ei满足等价标准J1和J2,计算误差ei的加权值ew;若该ew为最小值ewmin,将参数θ(kP,kI,kD)的值作为最优解输出;其中Y1i与Y2i的数量均为2;本发明能够在双轴励磁同步调相机励磁系统的参数解集簇中找到最优解,并且在可接受的精度范围内尽量提高辨识速度。

Description

一种双轴励磁同步调相机励磁系统参数辨识方法
技术领域
本发明涉及一种系统参数辨识方法,尤其涉及一种调相机励磁系统参数辨识方法。
背景技术
为了在运行时能获取电机参数的准确信息,对电机施加特定的励磁信号,采集系统反馈量进行相应的数据处理,得到电机参数后利用程序将电机内的参数更新。这种在电机处于非正常运行状态下测定电机参数的方法,称为参数辨识。
发明内容
本发明的目的在于提供一种双轴励磁同步调相机励磁系统参数辨识方法,其能够在双轴励磁调相机励磁系统的参数解簇中找到最优解,并且在一定精度范围内提高辨识速度。
为了实现上述目的,本发明提出了一种双轴励磁同步调相机励磁系统参数辨识方法,该方法用于输出双轴励磁同步调相机励磁系统的参数θ(kP,kI,kD)的最优解,包括步骤:构建包含参数θ(kP,kI,kD)的双轴励磁同步调相机励磁系统的励磁系统仿真模型M1;参考《同步发电机励磁系统建模导则》确定有功功率等价标准J1: Pst∈[-5%PN,5%PN],PN为额定有功功率;确定无功功率等价标准J2: Qst∈[-2%QN,2%QN],QN为额定无功功率;对所述双轴励磁同步调相励磁系统以信号X1和X2作为扰动信号进行扰动试验;激励信号常选用阶跃信号、白噪声信号或伪随机二进制序列PRBS及三角波信号等;获取各真实响应指标Y2i,包括真实无功功率Q2,真实有功功率P2;其特征在于,还包括步骤:确定真实有功功率P2与仿真有功功率P1之间的有功功率误差e1,真实无功功率Q2与仿真无功功率Q1之间的误差e2分别对应的权重Wi,其中∑Wi=1;以及
S1、对所述参数θ(kP,kI,kD)的值进行选取,将该参数θ(kP,kI,kD)的值代入所述励磁系统仿真模型M1,并对该励磁系统仿真模型M1以所述信号X1和X2作为扰动信号进行扰动试验,获取各仿真响应指标Y1i,包括仿真无功功率Q1,仿真有功功率P1
S2、计算各仿真响应指标Y1i与各真实响应指标Y2i之间的误差ei,包括真实有功功率P2与仿真有功功率P1之间的有功功率误差e1,真实无功功率Q2与仿真无功功率Q1之间的误差e2;若该误差ei满足等价标准J1和J2,则进入S3,否则返回步骤S1;
S3、基于所述权重Wi计算所述误差ei的加权值ew,其中ew=∑ei·Wi;若该ew为最小值ewmin,则进入步骤S4,否则返回步骤S1;
S4、将所述参数θ(kP,kI,kD)的值作为最优解输出;
其中,各真实响应指标Y2i与各仿真响应指标Y1i的数量均为2。
本发明所述的双轴励磁同步调相机励磁系统参数辨识方法的构思是,利用导则约束和最优解约束两步判定来寻找参数θ的最优解;其中导则约束用于使误差ei保持在一定精度范围内从而平衡了精度和速度,优化了算法;,最优解约束用于提高辨识速度,用加权法将2个误差ei统一用单个加权值ew进行衡量,从而简化了算法,在提高辨识速度的同时也方便了根据对各误差ei重要程度的要求的不同而对其赋予相应的不同的权重Wi,满足了实际需求。
关于本发明的数学模型
设系统模型参数为θ时,各项指标的仿真值与试验值的误差为ei,ei是θ的函数 ei=F(θ);
根据需求,各项指标的权重分别为Wi,则指标误差加权值为:
ew=∑ei·Wi
其中权重满足其中∑Wi=1。本发明以上述方程为判定函数,旨在寻找判定函数达到最小值时相对应的系统模型参数θopt
最优解判定原则的数学模型为
Figure BDA0002354343270000021
其须满足两个约束条件:一为导则约束,二为最优解约束,
导则约束为:
等价标准J1
Pst∈[-5%PN,5%PN]
所述等价标准J2
Qst∈[-2%QN,2%QN]
其中,Pst为有功功率等价标准,ewmin为无功功率等价标准;
最优解约束为:
Figure BDA0002354343270000031
下面描述本发明进一步的技术方案
进一步地,在本发明所述的双轴励磁同步调相机励磁系统参数辨识方法中,所述信号X1和X2可为阶跃信号,白噪声信号或伪随机二进制序列PRBS及三角波信号。
通过对比仿真模型与实际模型的响应指标来进行励磁系统的参数辨识。
进一步地,在上述双轴励磁同步调相机励磁系统参数辨识方法中的一个实例中,所述各真实响应指标Y2i与各仿真响应指标Y1i的数量均为2,所述各真实响应指标Y2i包括真实有功功率P2、真实无功功率Q2
所述各仿真响应指标Y1i包括仿真有功功率P1、仿真无功功率Q1
所述误差ei包括真实有功功率P2与仿真有功功率P1之间的有功功率误差e1、真实无功功率Q2与仿真无功功率Q1之间的无功功率误差e2
所述权重Wi包括有功功率权重W1、无功功率权重W2
等价标准J1为:
Pst∈[-5%PN,5%PN]
所述等价标准J2为:
Qst∈[-2%QN,2%QN]
其中,Pst为有功功率等价标准,Qst为无功功率等价标准。
所述误差ei满足等价标准J1和J2是指有功功率误差e1、无功功率误差e2同时分别对应满足有功功率等价标准Pst、无功功率等价标准Qst
进一步地,在本发明所述的双轴励磁同步调相机励磁系统参数辨识方法中,所述参数θ包括比例参数Kp、积分参数Ki以及微分参数KD
所述参数θ为控制环节PID参数
本发明具有以下优点:
(1)引入了最优解约束,使得本发明方法的各相指标误差ei比传统人工辨识方法更小,辨识结果更优
(2)本发明方法以最小值eWmin确定最优解,具有更好的准确性和鲁棒性。
(3)本发明综合使用导则约束和最优解约束。导则约束使误差ei保持在可接受的精度范围内,从而平衡了精度和速度,优化了算法;最优解约束利用加权法简化算法,提高了辨识速度,同时方便分别对各误差ei重要程度的要求的不同而对其赋予相应不用的权重,满足实际需求。
附图说明
图1为系统辨识原理图;
图2为本发明所述的双轴励磁同步调相机励磁系统参数辨识方法在一种实施方式下的流程图。
具体实施方式
下面将结合说明书和具体的实例对本发明所述的双轴励磁同步调相机励磁参数系统辨识方法做出进一步的解释说明。
本实例中,待辨识的双轴励磁同步调相机励磁系统采用直流励磁方式,待辨识参数θ为控制环节PID参数,包括比例参数KP积分参数KI和微分参数KD,扰动试验的扰动信号X1和X2常选用阶跃信号、白噪声信号或伪随机二进制序列PRBS及三角波信号等。本发明采用经典的阶跃响应辨识方法,具有直观性强,操作简单,试验时间短等优点。相应指标考察有功功率、无功功率两个,即相应响应指标数量N取2,因此,本实例中i=1,2,。
如图2所述,本实施例的方法对上述待辨识的双轴励磁同步调相机励磁系统的参数θ进行辨识,输出其最优解,包括步骤:
A1、对待辨识的双轴励磁同步调相机励磁系统以阶跃信号X1和X2作为扰动信号进行空载扰动试验,获取各真实响应指标Y2i,包括真实有功功率P2、真实无功功率 Q2
A2、规则制定,包括:
构建待辨识的双轴励磁同步调相机励磁系统的励磁系统仿真模型M1
等价标准J1为:Pst∈[-5%PN,5%PN];所述等价标准J2为:Qst∈[-2%QN,2%QN]
其中,Pst为有功功率等价标准,Qst为无功功率等价标准;
确定各真实响应指标YRi对应的权重Wi,包括有功功率权重W1、无功功率权重 W2,本实施例中有功功率权重W1设为0.2,无功功率权重W2设为0.8,符合∑Wi=1。
A3、参数调整,包括:
对参数θ的值进行选取,将该参数θ的值代入所述励磁系统仿真模型M1
A4、标准判定,包括:
对励磁系统仿真模型M1真以阶跃信号X1和X2以作为扰动信号进行扰动试验,验获取各仿真响应指标Y1i,包括仿真有功功率P1仿真无功功率Q1
计算各仿真响应指标Y1i与各真实响应指标Y2i之间的误差ei,包括真实有功功率P2与仿真有功功率P1之间的有功功率误差e1,真实无功功率Q2与仿真无功功率Q1之间的无功功率误差e2;若该误差ei满足等价标准J1和J2,即有功功率误差e1,无功功率误差e2分别对应满足有功功率等价标准Pst,无功功率等价标准Qst,则进入步骤A5,否则,返回步骤A3。
A5、最优解判定,包括:
基于权重Wi计算误差ei的加权值ew,其中ew=∑ei·Wi
若该ew为最小值ewmin,则进入步骤A6,否则返回步骤A3;
A6、输出参数θ数包括优化比例参数KP积分参数KI和微分参数KD
本发明提出的双轴励磁同步调相机励磁系统参数辨识方法具有更好的准确性和鲁棒性,在工程中有着一定的优越性和实际应用价值。要注意的是,以上列举仅为本发明的具体实施例,显然本发明不限于以上实施例,最优有着许多类似变化。本领域的技术人员如果从本发明公开的内容直接导出或联想到的所有变形,均应属于本发明的保护范围。

Claims (5)

1.一种双轴励磁同步调相机励磁系统参数辨识方法,该方法用于输出双轴励磁同步调相机励磁系统的参数θ(kP,kI,kD)的最优解,包括步骤:构建包含参数θ(kP,kI,kD)的双轴励磁同步调相机励磁系统的励磁系统仿真模型M1;确定有功功率等价标准J1:Pst∈[-5%PN,5%PN],PN为额定有功功率;确定无功功率等价标准J2:Qst∈[-2%QN,2%QN],QN为额定无功功率;对双轴励磁同步调相励磁系统以信号X1和X2作为扰动信号进行扰动试验,获取各真实响应指标Y2i,包括真实无功功率Q2,真实有功功率P2;其特征在于,还包括步骤:确定真实有功功率P2与仿真有功功率P1之间的有功功率误差e1,真实无功功率Q2与仿真无功功率Q1之间的误差e2分别对应的权重Wi,其中∑Wi=1;以及
S1、对所述参数θ(kP,kI,kD)的值进行选取,将该参数θ(kP,kI,kD)的值代入所述励磁系统仿真模型M1,并对该励磁系统仿真模型M1以所述信号X1和X2作为扰动信号进行扰动试验,获取各仿真响应指标Y1i,包括仿真无功功率Q1,仿真有功功率P1
S2、计算各仿真响应指标Y1i与各真实响应指标Y2i之间的误差ei,包括真实有功功率P2与仿真有功功率P1之间的有功功率误差e1,真实无功功率Q2与仿真无功功率Q1之间的误差e2;若该误差ei满足等价标准J1和J2,则进入S3,否则返回步骤S1;
S3、基于所述权重Wi计算所述误差ei的加权值ew,其中ew=∑ei·Wi;若该ew为最小值ewmin,则进入步骤S4,否则返回步骤S1;
S4、将所述参数θ(kP,kI,kD)的值作为最优解输出;
其中,各真实响应指标Y2i与各仿真响应指标Y1i的数量均为2。
2.如权利要求1所述的双轴励磁同步调相机励磁系统参数识别方法,其特征在于,所述信号X1和X2可以为阶跃信号,白噪声信号或伪随机二进制序列PRBS及三角波信号。
3.如权利要求1所述的双轴励磁同步调相机励磁系统参数辨识方法,其特征在于,所述各真实响应指标Y2i与各仿真响应指标Y1i的数量均为2,所述各真实响应指标Y2i包括真实无功功率Q2,真实有功功率P2;各仿真响应指标Y1i包括仿真无功功率Q1,仿真有功功率P1
所述误差ei包括真实有功功率P2与仿真有功功率P1之间的有功功率误差e1,真实无功功率Q2与仿真无功功率Q1之间的误差e2;所述权重Wi包括有功功率权重W1,无功功率权重W2
4.如权利要求1所述的双轴励磁同步调相机励磁系统参数辨识方法,其特征在于,所述等价标准J1为:
Pst∈[-5%PN,5%PN],PN为额定有功功率
所述等价标准J2为:
Qst∈[-2%QN,2%QN],QN为额定无功功率
其中,Pst为有功功率等价标准,Qst为无功功率等价标准;所述误差ei满足等价标准J1和J2是指有功功率误差e1,无功功率误差e2同时分别对应满足有功功率等价标准Pst,无功功率等价标准Qst
5.如权利要求1所述的双轴励磁同步调相机励磁系统参数辨识方法,其特征在于,所述参数θ(kP,kI,kD)包括比例参数KP,积分参数KI和微分参数KD
CN202010003559.5A 2020-01-03 2020-01-03 一种双轴励磁同步调相机励磁系统参数辨识方法 Withdrawn CN111880091A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010003559.5A CN111880091A (zh) 2020-01-03 2020-01-03 一种双轴励磁同步调相机励磁系统参数辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010003559.5A CN111880091A (zh) 2020-01-03 2020-01-03 一种双轴励磁同步调相机励磁系统参数辨识方法

Publications (1)

Publication Number Publication Date
CN111880091A true CN111880091A (zh) 2020-11-03

Family

ID=73153930

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010003559.5A Withdrawn CN111880091A (zh) 2020-01-03 2020-01-03 一种双轴励磁同步调相机励磁系统参数辨识方法

Country Status (1)

Country Link
CN (1) CN111880091A (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005275588A (ja) * 2004-03-23 2005-10-06 Yaskawa Electric Corp 電動機制御装置制御パラメータ感度解析装置
WO2014138896A1 (en) * 2013-03-13 2014-09-18 University Of Saskatchewan Adaptive control technique for transmission line control devices to handle large disturbance scenarios
CN104378038A (zh) * 2014-11-18 2015-02-25 东华大学 基于人工神经网络的永磁同步电机参数辨识方法
US20150100282A1 (en) * 2012-10-03 2015-04-09 Operation Technology, Inc. Generator dynamic model parameter estimation and tuning using online data and subspace state space model
CN105610360A (zh) * 2016-01-22 2016-05-25 国网上海市电力公司 一种同步发电机励磁系统参数辨识方法
CN106097157A (zh) * 2016-07-22 2016-11-09 同济大学 一种同步发电机励磁系统模型参数辨识优化方法
JP2017034739A (ja) * 2015-07-28 2017-02-09 株式会社ダイヘン 系統インピーダンス推定装置、インバータ装置、および、系統インピーダンス推定方法
WO2018006499A1 (zh) * 2016-07-06 2018-01-11 南方电网科学研究院有限责任公司 电力系统的主导失稳模式识别方法和系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005275588A (ja) * 2004-03-23 2005-10-06 Yaskawa Electric Corp 電動機制御装置制御パラメータ感度解析装置
US20150100282A1 (en) * 2012-10-03 2015-04-09 Operation Technology, Inc. Generator dynamic model parameter estimation and tuning using online data and subspace state space model
WO2014138896A1 (en) * 2013-03-13 2014-09-18 University Of Saskatchewan Adaptive control technique for transmission line control devices to handle large disturbance scenarios
CN104378038A (zh) * 2014-11-18 2015-02-25 东华大学 基于人工神经网络的永磁同步电机参数辨识方法
JP2017034739A (ja) * 2015-07-28 2017-02-09 株式会社ダイヘン 系統インピーダンス推定装置、インバータ装置、および、系統インピーダンス推定方法
CN105610360A (zh) * 2016-01-22 2016-05-25 国网上海市电力公司 一种同步发电机励磁系统参数辨识方法
WO2018006499A1 (zh) * 2016-07-06 2018-01-11 南方电网科学研究院有限责任公司 电力系统的主导失稳模式识别方法和系统
CN106097157A (zh) * 2016-07-22 2016-11-09 同济大学 一种同步发电机励磁系统模型参数辨识优化方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
ANTONIO RAFFO等: "Electron Device Model Parameter Identification Through Large-Signal-Predictive Small-Signal-Based Error Functions", 《IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES 》 *
刘元尊等: "基于参数辨识的波浪发电场等效建模研究", 《电力工程技术》 *
刘立群等: "基于MMPSO算法的风电场多机等值模型参数辨识", 《可再生能源》 *
徐强等: "采样数据对综合负荷模型参数辨识的影响", 《江苏电机工程》 *
李慧磊等: "水轮发电机组励磁系统参数辨识", 《全国大中型水电厂技术协作网技术交流文集》 *
石建飞等: "永磁同步电机在线参数辨识方法研究", 《电机与控制学报》 *

Similar Documents

Publication Publication Date Title
Wulkow The simulation of molecular weight distributions in polyreaction kinetics by discrete Galerkin methods
CN107783758B (zh) 一种智能合约工程方法
Choi et al. Multiple solutions and plateau phenomenon in measurement-based load model development: Issues and suggestions
Gomes et al. HintCO: Hint-based Configuration of Co-simulations
CN113505525B (zh) 基于微分神经网络的电力系统动态元件建模方法及装置
CN104461906A (zh) 一种基于代码转换的mc/dc覆盖准则的软件测试数据生成方法
CN109325281A (zh) 一种舰船管路三维设计方法和装置
CN111880091A (zh) 一种双轴励磁同步调相机励磁系统参数辨识方法
CN111258222A (zh) 自回归滑动平均系统自适应状态估计方法及闭环控制系统
Dierkes Formal analysis of a triplex sensor voter in an industrial context
CN115993138A (zh) 一种多表冗余激光惯组故障诊断方法及系统
CN115712564A (zh) 一种scade软件模型复杂度计算方法
CN114861984A (zh) 用于预测高含co2凝析气藏油环体积的方法及处理器
CN113254435B (zh) 一种数据增强方法及系统
Nazir et al. Testability estimation model (TEM OOD)
CN109933948B (zh) 一种形式验证方法、装置、形式验证平台及可读存储介质
CN113987261A (zh) 一种基于动态信任感知的视频推荐方法及系统
Donzé BreachFlows: simulation-based design with formal requirements for industrial CPS
Saal et al. 3D Model-Based Product Definition and Production–a Mind Change with Technical Hurdles
Olga Criterion of identification of test steps in training and educational systems.
Kaur et al. Comparative analysis of the software effort estimation models
Rask et al. A COMPARISON OF ALBRECHT'S FUNCTION POINT AND SYMONS'MARK II METRICS
CN113282884B (zh) 通用根因分析方法
CN114779625B (zh) 基于vrft的pd控制器设计方法及装置、电子设备
CN114218819B (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
WW01 Invention patent application withdrawn after publication

Application publication date: 20201103

WW01 Invention patent application withdrawn after publication