CN109725254B - 一种实心转子大型调相机动态参数辨识方法 - Google Patents

一种实心转子大型调相机动态参数辨识方法 Download PDF

Info

Publication number
CN109725254B
CN109725254B CN201810383015.9A CN201810383015A CN109725254B CN 109725254 B CN109725254 B CN 109725254B CN 201810383015 A CN201810383015 A CN 201810383015A CN 109725254 B CN109725254 B CN 109725254B
Authority
CN
China
Prior art keywords
axis
reactance
operational
current
rotor
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
CN201810383015.9A
Other languages
English (en)
Other versions
CN109725254A (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.)
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 CN201810383015.9A priority Critical patent/CN109725254B/zh
Publication of CN109725254A publication Critical patent/CN109725254A/zh
Application granted granted Critical
Publication of CN109725254B publication Critical patent/CN109725254B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种实心转子大型调相机动态参数辨识方法,包括:在实心转子大型调相机处于静止状态时,在电枢绕组a、b两相加频率的交流电压,并拖动转子,将实心转子大型调相机的转子固定在励磁绕组感应电压最小的位置;将励磁绕组短接,在电枢绕组a、b两相间输入sinc信号,并由得到的电枢绕组a相电流确定在当前转子位置下的实心转子大型调相机的d轴电流;确定考虑d轴存在多个阻尼绕组时,d轴运算电抗的表达式;将时域下的d轴电流进行频域变换,利用帕德逼近,将频域下的d轴电流表示成有理分式的形式,并根据d轴下的电压平衡关系得到d轴运算电抗,进而确定d轴各阶电抗与对应的时间常数。本发明简便易行、安全性好且精确度高。

Description

一种实心转子大型调相机动态参数辨识方法
技术领域
本发明属于电机试验的参数辨识领域,更具体地,涉及一种基于sinc信号电压注入与帕德逼近的实心转子大型调相机动态参数辨识方法。
背景技术
针对我国现有特高压直流输电工程中长期存在的动态无功补偿能力不足的问题,具有高动态响应能力且无功补偿容量大的实心转子大型调相机被视为解决此类问题的有效措施,并将广泛加装于特高压直流输电工程的换流站中。
实心转子大型调相机优良的动态无功输出特性,与其主机性能参数,特别是d轴的动态参数紧密相关。因此在实心转子大型的主机制造上,对其d轴动态参数提出了更高的要求。这些要求可以概括为:直轴各阶电抗尽量小,以保证动态时的无功输出幅值足够大;直轴次暂态时间常数尽量大,增大次暂态过程的动态无功输出;直轴暂态时间常数尽量小,保证实心转子大型调相机能够迅速响应其励磁系统的调节作用。
鉴于以上情况,实心转子大型调相机出厂时必须对影响其动态无功输出特性的d轴动态参数进行试验校核。由于实心转子大型调相机采用隐极实心转子等结构,阻尼效果很强,阻尼回路复杂,因此需要提出合理的动态参数试验方法及参数辨识方法,以对实心转子大型调相机d轴动态参数进行准确辨识。
实心转子大型调相机的制造是基于传统汽轮发电机的制造经验,并对其结构进行了较大程度的优化,属于电励磁同步电机范畴,目前尚未有针对于实心转子大型调相机的d轴动态参数辨识方法提出,仅能沿用电励磁同步电机的相关方法。电励磁同步电机的d轴动态参数试验与辨识方法是基于三相突然短路试验,根据其d轴短路电流公式,对得到的机端短路电流用最小二乘法或其他方法进行拟合辨识。
基于三相突然短路试验与最小二乘法的实心转子大型调相机d轴动态参数试验与辨识方法存在以下缺陷:(1)三相突然短路试验准备过程复杂,且会对电机本体形成较大损伤,由于新型调相机容量大、造价昂贵不适宜采用此类试验方法。(2)最小二乘法辨识d轴动态参数首先需要将机端三相突然短路电流根据Park变换,变换到dq0坐标系下,得到d轴电流。由于采用最小二乘法进行拟合的前提是需要求出短路时d轴电流的表达式,而表达式的求解过程繁琐且运用到了许多假设进行简化,因此精确度不高,在对新型调相机这类d轴参数精确度要求高的电机不适用。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种实心转子大型调相机动态参数辨识方法,由此解决采用三相突然短路试验以及最小二乘法在针对实心转子大型调相机d轴动态参数辨识存在的操作过程复杂、对电机造成损伤及精确度不高的技术问题。
为实现上述目的,本发明提供了一种实心转子大型调相机动态参数辨识方法,包括:
将所述实心转子大型调相机的转子固定在励磁绕组感应电压最小的位置;
将励磁绕组短接,在所述电枢绕组a、b两相间输入sinc信号,并由输入sinc信号后电枢绕组的a相电流确定在当前转子位置下的所述实心转子大型调相机的d轴电流,其中,所述电枢绕组的a相电流与所述d轴电流成正比;
确定d轴电压、d轴电流与d轴运算电抗的关系及在考虑d轴存在多个阻尼绕组时,所述d轴运算电抗的表达形式;
将时域下的所述d轴电流进行频域变换,利用帕德逼近,将频域下的d轴电流表示成有理分式的形式,并根据d轴电压、d轴电流与d轴运算电抗的关系得到所述d轴运算电抗的有理分式的表达形式,进而确定d轴各阶电抗与对应的时间常数。
优选地,所述将所述实心转子大型调相机的转子固定在励磁绕组感应电压最小的位置,包括:
在所述实心转子大型调相机的转子处于静止状态时,通过在所述实心转子大型调相机的电枢绕组a、b两相加交流电压,并拖动转子,将转子固定在励磁绕组感应电压最小的位置。
优选地,所述确定d轴电压、d轴电流与d轴运算电抗的关系及在考虑d轴存在多个阻尼绕组时,所述d轴运算电抗的表达形式,包括:
根据所述实心转子大型调相机的d轴的电压平衡关系得到d轴电压、d轴电流与d轴运算电抗的关系,由考虑d轴存在多个阻尼绕组的实心转子大型调相机的数学模型得到d轴运算电抗的表达形式。
优选地,所述d轴运算电抗的表达形式为:
Figure BDA0001641460840000031
其中,Xd表示d轴同步电抗,Td表示短路时间常数,Td0为开路时间常数。
优选地,所述将时域下的所述d轴电流进行频域变换,利用帕德逼近,将频域下的d轴电流表示成有理分式的形式,包括:
将时域下的所述d轴电流进行频域变换,基于d轴电压在输入sinc信号后的表达形式、d轴电压、d轴电流与d轴运算电抗的关系及所述d轴运算电抗的表达形式,将频域下的d轴电流表示为多阶幂级数之和;
利用帕德逼近,将表示为多阶幂级数之和形式的d轴电流表示成有理分式的形式。
优选地,所述根据d轴电压、d轴电流与d轴运算电抗的关系得到所述d轴运算电抗的有理分式的表达形式,进而确定d轴各阶电抗与对应的时间常数,包括:
根据d轴电流的有理分式表示及d轴电压、d轴电流与d轴运算电抗的关系得到所述d轴运算电抗的有理分式的表达形式;
根据所述d轴运算电抗的有理分式的表达形式及在考虑d轴存在多个阻尼绕组时,所述d轴运算电抗的表达形式确定d轴同步电抗Xd与d轴各阶时间常数Td与Td0
根据所述d轴运算电抗的倒数表达形式及所述d轴运算电抗的有理分式的表达形式确定d轴各阶电抗的值。
所述d轴运算电抗的倒数表达形式为:
Figure BDA0001641460840000041
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
提出了一种针对实心转子大型调相机d轴动态参数辨识方法,此方法在实心转子大型调相机处于静止状态时,向电枢两相绕组注入sinc信号,根据注入信号后的电枢电流响应,结合小波变换与帕德逼近即可辨识得到实心转子大型调相机d轴的动态参数。由于注入的sinc信号的电压幅值很小,此方法极大程度上避免了采用三相突然短路法可能对电机造成的损伤,试验便捷安全,辨识精度高,不仅适用于实心转子大型调相机,也适用于其他类型的同步电机。
附图说明
图1是本发明实施例公开的一种实心转子大型调相机动态参数辨识方法的流程示意图;
图2是本发明实施例公开的一种调相机的sinc信号注入法试验的接线图;
图3是本发明实施例公开的一种sinc信号的时域波形图;
图4是本发明实施例公开的一种sinc信号的频域波形图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明提供了一种简便易行、安全性好且精确度高的基于sinc信号注入试验和帕德逼近的新型调相机的d轴动态参数辨识方法。
如图1所示是本发明实施例公开的一种实心转子大型调相机动态参数辨识方法的流程示意图,在图1所示的方法中包括以下步骤:
(1)将实心转子大型调相机的转子固定在励磁绕组感应电压最小的位置;
在本发明实施例中,将实心转子大型调相机的转子固定在励磁绕组感应电压最小的位置可以通过以下方式实现:
在实心转子大型调相机的转子处于静止状态时,通过在实心转子大型调相机的电枢绕组a、b两相加交流电压,缓慢拖动转子,将转子固定在励磁绕组感应电压最小的位置。
(2)将励磁绕组短接,在电枢绕组a、b两相间输入sinc信号,并由输入sinc信号后电枢绕组的a相电流确定在当前转子位置下的实心转子大型调相机的d轴电流,其中,电枢绕组的a相电流与d轴电流成正比;
在本发明实施例中,在硬件实现上,可以利用信号发生器在实心转子大型调相机电枢绕组a、b两相间输入sinc信号电压,用录波仪对此时的定子a相电流录波,进而由录波得到的定子a相电流确定d轴电流。
其中,a相绕组电流幅值与d轴电流幅值的关系为:
Figure BDA0001641460840000061
其中,Id表示d轴电流幅值,Is表示a相绕组电流幅值。
(3)确定d轴电压、d轴电流与d轴运算电抗的关系及在考虑d轴存在多个阻尼绕组时,d轴运算电抗的表达形式;
在本发明实施例中,确定d轴电压、d轴电流与d轴运算电抗的关系及在考虑d轴存在多个阻尼绕组时,d轴运算电抗的表达形式可以通过以下方式实现:
根据实心转子大型调相机的d轴的电压平衡关系得到d轴电压、d轴电流与d轴运算电抗的关系,由考虑d轴存在多个阻尼绕组的实心转子大型调相机的数学模型得到d轴运算电抗的表达形式。
(4)将时域下的所述d轴电流进行频域变换,利用帕德逼近,将频域下的d轴电流表示成有理分式的形式,并根据d轴电压、d轴电流与d轴运算电抗的关系得到所述d轴运算电抗的有理分式的表达形式,进而确定d轴各阶电抗与对应的时间常数。
在本发明实施例中,将时域下的d轴电流进行频域变换,利用帕德逼近,将频域下的d轴电流表示成有理分式的形式可以通过以下方式实现:
将时域下的d轴电流进行频域变换,基于d轴电压在输入sinc信号后的表达形式、d轴电压、d轴电流与d轴运算电抗的关系及d轴运算电抗的表达形式,将频域下的d轴电流表示为多阶幂级数之和;
利用帕德逼近,将表示为多阶幂级数之和形式的d轴电流表示成有理分式的形式。
下面结合附图对本发明作进一步的详细说明。
附图2是实心转子大型调相机sinc信号注入试验的接线图,试验时首先让转子静止,在电枢a、b两相加一定频率的交流电压,并缓慢拖动转子,将转子位置定于励磁绕组感应电压最小的位置。然后,短接励磁绕组,sinc信号从定子a、b相绕组输入,用录波仪对此时的定子a相电流进行录波。sinc信号的时域与频域波形如附图3与附图4所示,由附图4可知,sinc信号在频域的低频段值为常数,对高频段的衰减十分明显,因此在频率较低的范围内,可以直接将sinc信号频域中的值视为常数。
由电机理论可知,在dq0坐标系下,d轴的电压平衡关系为式(1):
Ud(s)=sXd(s)Id(s)+rsId(s) (1)
其中,Ud(s)表示频域下的d轴电压,rs表示电枢绕组电阻。
考虑d轴有多阶阻尼绕组时,可以得到d轴运算电抗表达式Xd(s),如式(2)。
Figure BDA0001641460840000071
其中,Xd表示定子d轴同步电抗;Td表示短路时间常数;Td0为开路时间常数。
由于试验测得的电流是时域下的电流Id(t),考虑到所输入的sinc信号是脉冲形式,因此需要利用离散化的小波变换将时域下的电流Id(t)转化为频域下的Id(s)。
因为Ud(s)在sinc信号的截断频率之前是一个常数,再由式(1)与式(2)可知,Id(s)可以表示为多阶幂级数之和的形式,如式(3):
Figure BDA0001641460840000072
以幂级数表达的函数可以采用帕德逼近方法转化为用有理分式表达的函数,如式(4):
Figure BDA0001641460840000073
由式(4),已经将Id(s)转化为有理分式表达的形式,根据Xd(s)与Id(s)的关系,得到用有理分式表达的Xd(s),如式(5):
Figure BDA0001641460840000081
对比所得到的式(5)与式(2),即可确定出d轴同步电抗Xd与d轴各阶时间常数Td与Td0。再由于Xd(s)的倒数可以表示为:
Figure BDA0001641460840000082
再结合式(5)与式(6)可以得到d轴各阶电抗的值,至此已辨识出d轴的各阶电抗与其对应的短路与开路时间常数。由以上过程可知,采用基于sinc信号注入与帕德逼近的新型调相机d轴动态参数辨识方法可以有效辨识出实心转子大型调相机d轴的各阶电抗与时间常数,具有试验方法易行,辨识过程简单且精度高等优点,为包括实心转子大型调相机在内的同步电机的动态参数辨识提供了很好的思路。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种实心转子大型调相机动态参数辨识方法,其特征在于,包括:
将所述实心转子大型调相机的转子固定在励磁绕组感应电压最小的位置;
将励磁绕组短接,在电枢绕组a、b两相间输入sinc信号,并由输入sinc信号后电枢绕组的a相电流确定在当前转子位置下的所述实心转子大型调相机的d轴电流,其中,所述电枢绕组的a相电流与所述d轴电流成正比;
确定d轴电压、d轴电流与d轴运算电抗的关系及在考虑d轴存在多个阻尼绕组时,所述d轴运算电抗的表达形式;
将时域下的所述d轴电流进行频域变换,利用帕德逼近,将频域下的d轴电流表示成有理分式的形式,并根据d轴电压、d轴电流与d轴运算电抗的关系得到所述d轴运算电抗的有理分式的表达形式,进而确定d轴各阶电抗与对应的时间常数。
2.根据权利要求1所述的方法,其特征在于,所述将所述实心转子大型调相机的转子固定在励磁绕组感应电压最小的位置,包括:
在所述实心转子大型调相机的转子处于静止状态时,通过在所述实心转子大型调相机的电枢绕组a、b两相加交流电压,并拖动转子,将转子固定在励磁绕组感应电压最小的位置。
3.根据权利要求1或2所述的方法,其特征在于,所述确定d轴电压、d轴电流与d轴运算电抗的关系及在考虑d轴存在多个阻尼绕组时,所述d轴运算电抗的表达形式,包括:
根据所述实心转子大型调相机的d轴的电压平衡关系得到d轴电压、d轴电流与d轴运算电抗的关系,由考虑d轴存在多个阻尼绕组的实心转子大型调相机的数学模型得到d轴运算电抗的表达形式。
4.根据权利要求3所述的方法,其特征在于,所述d轴运算电抗的表达形式为:
Figure FDA0002401529370000021
其中,Xd表示d轴同步电抗,Td表示短路时间常数,Td0为开路时间常数。
5.根据权利要求4所述的方法,其特征在于,所述将时域下的所述d轴电流进行频域变换,利用帕德逼近,将频域下的d轴电流表示成有理分式的形式,包括:
将时域下的所述d轴电流进行频域变换,基于d轴电压在输入sinc信号后的表达形式、d轴电压、d轴电流与d轴运算电抗的关系及所述d轴运算电抗的表达形式,将频域下的d轴电流表示为多阶幂级数之和;
利用帕德逼近,将表示为多阶幂级数之和形式的d轴电流表示成有理分式的形式。
6.根据权利要求5所述的方法,其特征在于,所述根据d轴电压、d轴电流与d轴运算电抗的关系得到所述d轴运算电抗的有理分式的表达形式,进而确定d轴各阶电抗与对应的时间常数,包括:
根据d轴电流的有理分式表示及d轴电压、d轴电流与d轴运算电抗的关系得到所述d轴运算电抗的有理分式的表达形式;
根据所述d轴运算电抗的有理分式的表达形式及在考虑d轴存在多个阻尼绕组时,所述d轴运算电抗的表达形式确定d轴同步电抗Xd与d轴各阶时间常数Td与Td0
根据所述d轴运算电抗的倒数表达形式及所述d轴运算电抗的有理分式的表达形式确定d轴各阶电抗的值。
7.根据权利要求6所述的方法,其特征在于,所述d轴运算电抗的倒数表达形式为:
Figure FDA0002401529370000022
CN201810383015.9A 2018-04-26 2018-04-26 一种实心转子大型调相机动态参数辨识方法 Active CN109725254B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810383015.9A CN109725254B (zh) 2018-04-26 2018-04-26 一种实心转子大型调相机动态参数辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810383015.9A CN109725254B (zh) 2018-04-26 2018-04-26 一种实心转子大型调相机动态参数辨识方法

Publications (2)

Publication Number Publication Date
CN109725254A CN109725254A (zh) 2019-05-07
CN109725254B true CN109725254B (zh) 2020-05-19

Family

ID=66293726

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810383015.9A Active CN109725254B (zh) 2018-04-26 2018-04-26 一种实心转子大型调相机动态参数辨识方法

Country Status (1)

Country Link
CN (1) CN109725254B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102008063295A1 (de) * 2007-12-28 2009-07-30 Hitachi Industrial Equipment Systems Co. Ltd. Motorregeleinrichtung und Motorregelsystem
CN102073014A (zh) * 2011-01-07 2011-05-25 华北电网有限公司 一种同步发电机参数的辨识方法
CN102520353A (zh) * 2011-12-09 2012-06-27 清华大学 同步发电机模型参数分步辨识方法
CN103178769A (zh) * 2013-04-03 2013-06-26 哈尔滨工业大学 永磁同步电机静止状态下参数离线辨识方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102008063295A1 (de) * 2007-12-28 2009-07-30 Hitachi Industrial Equipment Systems Co. Ltd. Motorregeleinrichtung und Motorregelsystem
CN102073014A (zh) * 2011-01-07 2011-05-25 华北电网有限公司 一种同步发电机参数的辨识方法
CN102520353A (zh) * 2011-12-09 2012-06-27 清华大学 同步发电机模型参数分步辨识方法
CN103178769A (zh) * 2013-04-03 2013-06-26 哈尔滨工业大学 永磁同步电机静止状态下参数离线辨识方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"基于 Park 模型的同步发电机参数辨识";孙黎霞;《中国电机工程学报》;20090705;第29卷(第19期);第50-55页 *
"基于双模态自适应小波粒子群的永磁同步电机";赵海森 等;《中国电机工程学报》;20141025;第34卷(第30期);第5386-5393页 *

Also Published As

Publication number Publication date
CN109725254A (zh) 2019-05-07

Similar Documents

Publication Publication Date Title
Prince et al. Modeling, parameter measurement, and control of PMSG-based grid-connected wind energy conversion system
CN107167695B (zh) 基于磁场分布监测的永磁同步电机匝间短路故障诊断方法
CN110518852B (zh) 基于谐波注入的永磁同步电机电流谐波抑制方法
CN110146777A (zh) 一种发电机/调相机转子匝间短路故障检测方法
Sarma et al. Modeling, analysis, and validation of controller signal interharmonic effects in DFIG drives
CN113783488B (zh) 永磁同步电机全参数辨识方法及永磁同步电机系统
CN109725254B (zh) 一种实心转子大型调相机动态参数辨识方法
CN109471026B (zh) 一种实心转子大型调相机多阻尼回路动态参数辨识方法
Refaat et al. ANN-based diagnosis of incipient stator winding turn faults for three-phase induction motors in the presence of unbalanced supply voltage
CN109932648B (zh) 一种同步电机q轴电感饱和特性试验测量方法
CN109490774B (zh) 一种实心转子大型调相机动态参数试验测量方法
CN104267243A (zh) 同步发电机内电势和电抗参数的测量方法及装置
CN2684208Y (zh) 一种实现异步电机真正空载的实验装置
CN114567222A (zh) 永磁同步电机静态初始位置估算方法及定子电感辨识方法
Zheng et al. Diagnostic strategy and modeling of PMSM stator winding fault in electric vehicles
CN108107315B (zh) 一种凸极永磁同步电机定子绕组匝间短路抗干扰故障诊断方法及系统
Liu et al. Analysis of direct power control strategies applied to doubly fed induction generator
Zhang et al. Research on output waveform of generator with rectifier load considering commutation overlap angle
Chen The study on source of vibration and acoustic noise of permanent magnet machines by inverter
CN109782172A (zh) 一种同步电机参数试验测量方法
Dai et al. Speed sensorless control of electrically excited synchronous motor based on the reduce-order state observer flux observation
Ma et al. A New Parameters Identification Method for Novel Synchronous Condenser Based on Particle Swarm Algorithm and Wavelet Transform
CN112636326B (zh) 风电及柔直系统的次同步谐振控制方法及装置
CN114244215B (zh) 基于电流注入的电机谐波转矩抑制方法
Zhang et al. Study on permanent magnet synchronous demagnetization fault performance

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