会员体验
专利管家(专利管理)
工作空间(专利管理)
风险监控(情报监控)
数据分析(专利分析)
侵权分析(诉讼无效)
联系我们
交流群
官方交流:
QQ群: 891211   
微信请扫码    >>>
现在联系顾问~
首页 / 专利库 / 风能 / 风力发电场 / 风力发电机组 / 双馈感应发电机 / 一种双馈感应发电机参数辨识方法

一种双馈感应发电机参数辨识方法

阅读:933发布:2020-05-13

IPRDB可以提供一种双馈感应发电机参数辨识方法专利检索,专利查询,专利分析的服务。并且本发明公开了一种双馈感应发电机(DFIG)参数辨识方法,属电力系统建模领域。首先,确定双馈感应发电机在端口三相金属性短路情况下定子轴短路电流的解析表达式,根据该解析表达式确定待辨识参数;基于Matlab/Simulink平台搭建含双馈风电机组的单机无穷大系统,在双馈感应发电机端口设置三相金属性短路,仿真获得定子轴短路电流的受扰轨线;采用轨迹灵敏度方法,分析各待辨识参数的可辨识性以及参数辨识的难易程度;最后采用粒子群优化算法辨识各待辨识参数,基于多次辨识结果确定参数辨识值。本发明基于定子轴短路电流进行DFIG参数辨识,与基于有功功率和无功功率相比,受变流器控制环节的影响小;采用多次参数辨识结果确定参数辨识值,避免了单次辨识可能导致的误差大的问题。,下面是一种双馈感应发电机参数辨识方法专利的具体信息内容。

1.一种双馈感应发电机参数辨识方法,其特征在于,包括以下具体步骤:步骤1,根据双馈感应发电机DFIG在端口三相金属性短路情况下的定子dq轴短路电流的解析表达式,确定待辨识参数;

步骤2,基于Matlab/Simulink平台搭建含双馈风电机组的单机无穷大系统,在双馈感应发电机端口设置三相金属性短路故障,仿真获得定子dq轴短路电流的受扰轨线;

步骤3,采用轨迹灵敏度方法,以定子dq轴短路电流作为观测量,计算待辨识参数的轨迹灵敏度;

步骤4,根据步骤3中计算得到的轨迹灵敏度曲线,判断待辨识参数的可辨识性;

步骤5,计算短路后不同时间窗口内各待辨识参数的轨迹灵敏度的绝对平均值,选择轨迹灵敏度的绝对平均值最大的时间窗口进行参数辨识;

步骤6,根据步骤5选择的时间窗口,采用粒子群优化算法进行参数辨识,辨识的目标函数为误差指标err达到最小,即 其中,K为辨识时间窗口内定子dq轴短路电流的总点数;isd_est(n)为辨识时间窗口内第n个基于辨识值的定子d轴短路电流;isd(n)为辨识时间窗口内第n个定子d轴短路电流实际值;isq_est(n)为辨识时间窗口内第n个基于辨识值的定子q轴短路电流;isq(n)为辨识时间窗口内第n个定子q轴短路电流实际值;

步骤7,单独进行多次参数辨识,计算多次参数辨识结果的均值,剔除与均值偏差大于

100%的辨识结果;再将剩下的辨识结果重新求取均值,作为DFIG各待辨识参数的最终辨识结果。

2.根据权利要求1所述的一种双馈感应发电机参数辨识方法,其特征在于,步骤1中定子dq轴短路电流的解析表达式为:式中,isd和isq分别为定子d轴短路电流和q轴短路电流;is∞和δsu分别为定子侧稳态电流的幅值和相位;M为暂态电感, Ls为定子自感,Ls=Lσs+Lm,Lσs为定子漏感,Lm为定转子互感,Lr为转子自感,Lr=Lσr+Lm,Lσr为转子漏感; 和 分别为定子和转子磁链幅值的初值;ωs为同步转速;t为时间;δ1和δ2分别为定子和转子磁链相位的初值;TS和Tr分别为定子和转子衰减时间常数, ω0为额定同步转速;Rs为定子电阻,Rr为转子电阻;s为滑差。

3.根据权利要求2所述的一种双馈感应发电机参数辨识方法,其特征在于,步骤1中待辨识参数为Rs、Rr、LsLr和Lm。

4.根据权利要求1所述的一种双馈感应发电机参数辨识方法,其特征在于,步骤2中双馈风电机组采用相量模型phasor model。

5.根据权利要求1所述的一种双馈感应发电机参数辨识方法,其特征在于,步骤3中轨迹灵敏度的计算方法为:a:将一个待辨识参数θ的数值增大Δθ,仿真得到定子dq轴短路电流受扰轨迹y,其中,Δθ为增量,Δθ=10%θ0,θ0为待辨识参数θ的初始值;

b:将待辨识参数θ的数值减少Δθ,仿真得到定子dq轴短路电流受扰轨迹y′;

c:根据公式 计算θ的轨迹灵敏度Sj,其中,K为辨识时间

窗口内定子dq轴短路电流的总点数;y(n)和y′(n)分别为第n个定子d轴和q轴短路电流。

6.根据权利要求1所述的一种双馈感应发电机参数辨识方法,其特征在于,步骤4中待辨识参数可辨识性的判断方法为:如果有大于或等于2个待辨识参数的轨迹灵敏度曲线同时过零点,则判定这些参数不是唯一可辨识;如果所有参数的灵敏度都不同时过零,则判定这些参数唯一可辨识。

7.根据权利要求1所述的一种双馈感应发电机参数辨识方法,其特征在于,步骤3中选择变流器控制环节动作结束后的时间窗口进行轨迹灵敏度计算。

说明书全文

一种双馈感应发电机参数辨识方法

技术领域

[0001] 本发明属于电力系统建模领域,特别涉及一种双馈感应发电机参数辨识方法。

背景技术

[0002] 建模是电力系统计算、分析以及运行控制的基础。合理而精准的风电机组模型参数对电力系统电磁暂态和机电暂态仿真分析的正确性尤为重要,如何得到更接近实际的风电机组模型参数,是电力学术界和工业界的研究热点。
[0003] 在辨识双馈感应发电机DFIG参数时,目前一般基于风电机组在端口电压跌落下有功功率和无功功率的受扰轨线进行参数辨识。但由于变流器控制环节动态与DFIG动态耦合较强,电压跌落下风电机组出口的有功功率、无功功率受扰轨线动态不仅与DFIG密切相关,受变流器控制环节的影响也较大。因此基于有功功率和无功功率辨识DFIG参数时,如果控制器参数给定值不准确,将导致发电机参数的较大辨识误差。

发明内容

[0004] 本发明所要解决的技术问题是:针对目前DFIG参数辨识精度受变流器控制环节影响的问题,提出一种基于定子dq轴短路电流作为观测量的参数辨识新方法。本发明提出的双馈感应发电机参数辨识方法,基于定子dq轴短路电流进行参数辨识;本发明还基于多次辨识结果获得参数辨识值,避免了单次辨识可能带来误差大的问题。
[0005] 本发明为解决上述技术问题采用以下技术方案:
[0006] 本发明提供一种双馈感应发电机参数辨识方法,包括以下具体步骤:
[0007] 步骤1,根据双馈感应发电机在端口三相金属性短路情况下的定子dq轴短路电流的解析表达式,确定待辨识参数;
[0008] 步骤2,基于Matlab/Simulink平台搭建含双馈风电机组的单机无穷大系统,在双馈感应发电机端口设置三相金属性短路故障,仿真获得定子dq轴短路电流的受扰轨线;
[0009] 步骤3,采用轨迹灵敏度方法,以定子dq轴短路电流作为观测量,计算待辨识参数的轨迹灵敏度;
[0010] 步骤4,根据步骤3中计算得到的轨迹灵敏度曲线,判断待辨识参数的可辨识性;
[0011] 步骤5,计算短路后不同时间窗口内各待辨识参数的轨迹灵敏度的绝对平均值,选择轨迹灵敏度的绝对平均值最大的时间窗口进行参数辨识;
[0012] 步骤6,根据步骤5选择的时间窗口,采用粒子群优化算法进行参数辨识,辨识的目标函数为误差指标err达到最小,即 其中,K为辨识时间窗口内定子dq轴短路电流的总点数;isd_est(n)为辨识时间窗口内第n个基于辨识值的定子d轴短路电流;isd(n)为辨识时间窗口内第n个定子d轴短路电流实际值;
isq_est(n)为辨识时间窗口内第n个基于辨识值的定子q轴短路电流;isq(n)为辨识时间窗口内第n个定子q轴短路电流实际值;
[0013] 步骤7,单独进行多次参数辨识,计算多次参数辨识结果的均值,剔除与均值偏差大于100%的辨识结果;再将剩下的辨识结果重新求取均值,作为DFIG各待辨识参数的最终辨识结果。
[0014] 作为本发明的进一步优化方案,步骤1中定子dq轴短路电流的解析表达式为:
[0015]
[0016]
[0017] 式中,isd和isq分别为定子d轴短路电流和q轴短路电流;is∞和δsu分别为定子侧稳态电流的幅值和相位;M为暂态电感, Ls为定子自感,Ls=Lσs+Lm,Lσs为定子漏感,Lm为定转子互感,Lr为转子自感,Lr=Lσr+Lm,Lσr为转子漏感; 和 分别为定子和转子磁链幅值的初值;ωs为同步转速;t为时间;δ1和δ2分别为定子和转子磁链相位的初值;TS和Tr分别为定子和转子衰减时间常数, ω0为额定同步转速;Rs为定子电阻,Rr为转子电阻;s为滑差。
[0018] 作为本发明的进一步优化方案,步骤1中待辨识参数为Rs、Rr、LsLr和Lm。
[0019] 作为本发明的进一步优化方案,步骤2中双馈风电机组采用相量模型phasor model。
[0020] 作为本发明的进一步优化方案,步骤3中轨迹灵敏度的计算方法为:
[0021] a:将一个待辨识参数θ的数值增大Δθ,仿真得到定子dq轴短路电流受扰轨迹y,其中,Δθ为增量,Δθ=10%θ0,θ0为待辨识参数θ的初始值;
[0022] b:将待辨识参数θ的数值减少Δθ,仿真得到定子dq轴短路电流受扰轨迹y′;
[0023] c:根据公式 计算θ的轨迹灵敏度Sj,其中,K为辨识时间窗口内定子dq轴短路电流的总点数;y(n)和y′(n)分别为第n个定子d轴和q轴短路电流。
[0024] 作为本发明的进一步优化方案,步骤4中待辨识参数可辨识性的判断方法为:如果有大于或等于2个待辨识参数的轨迹灵敏度曲线同时过零点,则判定这些参数不是唯一可辨识;如果所有参数的灵敏度都不同时过零,则判定这些参数唯一可辨识。
[0025] 作为本发明的进一步优化方案,步骤3中选择变流器控制环节动作结束后的时间窗口进行轨迹灵敏度计算。
[0026] 本发明采用以上技术方案与现有技术相比,具有以下技术效果:本发明提出采用双馈感应发电机在端口三相短路情况下定子dq轴短路电流进行参数辨识,受变流器控制环节的参数影响小,辨识精度高;采用多次参数辨识结果确定参数辨识值,避免了单次辨识可能导致的误差较大的问题,整体上提高了参数的辨识精度。

附图说明

[0027] 图1为含双馈风电机组的单机无穷大系统结构图。
[0028] 图2为双馈风电机组的结构图。
[0029] 图3为DFIG各参数的轨迹灵敏度,其中,(a)为定子电阻Rs的轨迹灵敏度,(b)转子电阻Rr的轨迹灵敏度,(c)定转子漏感Lσs、Lσr的轨迹灵敏度,(d)定转子互感Lm的轨迹灵敏度。
[0030] 图4为各待辨识参数的50次辨识结果,其中,(a)为定子电阻Rs的50次辨识结果,(b)为转子电阻Rr的50次辨识结果,(c)为定转子自感的乘积LsLr的50次辨识结果,(d)为定转子互感Lm的50次辨识结果。
[0031] 图5为基于参数辨识结果与实际短路电流的受扰轨线,其中,(a)为定子d轴电流isd的受扰轨线,(b)定子q轴电流isq的受扰轨线。
[0032] 图6为本发明的工作流程图。

具体实施方式

[0033] 下面结合附图以及具体实施例对本发明的技术方案做进一步的详细说明:
[0034] 仿真系统如图1所示,为一台基于DFIG的风电机组接入无穷大系统,该仿真系统搭建于Matlab 2013b软件平台,系统中所有元件均取自Matlab自带的风电场算例,双馈风电机组采用相量模型(phasor model),双馈感应发电机的参数见表1。
[0035] 表1 DFIG仿真参数
[0036]参数名 数值
额定功率/MW 1.5
额定电压/kV 0.575
初始转差/pu 0.2
定子电阻/pu 0.00706
转子电阻/pu 0.005
定子漏感/pu 0.171
转子漏感/pu 0.156
互感/pu 2.9
[0037] 双馈风电机组的结构如图2所示,其中,β为风力机叶片的桨距角,ωr为DFIG转子转速,ir为转子绕组中的电流(以下简称“转子电流”),vdr和vqr分别为转子励磁电压的d轴和q轴分量;vDC代表直流侧电容的电压;vg、ig分别代表网侧变流器的输出电压、电流;vdg、vqg分别代表网侧变流器输出电压的d、q轴分量。
[0038] 在如图1所示的双馈感应发电机端口B1母线处设置三相金属性短路故障,在时间t=1s时,B1母线电压降为零,故障持续0.2s,在时间t=1.2s故障消失,系统恢复至原状态。
[0039] 如图6所示,本发明提供的DFIG参数辨识的具体步骤如下:
[0040] 步骤1:根据双馈风电机组的定子dq轴短路电流解析表达式,分析参数的可辨识性。
[0041] 定子dq轴短路电流表达式如下:
[0042]
[0043]
[0044] 式中,is∞和δsu分别为定子侧稳态电流的幅值与相位;M为暂态电感,Ls为定子自感,Ls=Lσs+Lm,Lσs为定子漏感,Lm为定转子互感,Lr为转子自感,Lr=Lσr+Lm,Lσr为转子漏感; 和 分别为定子和转子磁链幅值的初值;ωs为同步转速;t为时间;δ1和δ2分别为定子和转子磁链相位的初值;TS和Tr分别为定子和转子衰减时间常数, ω0
为额定同步转速;Rs为定子电阻, Rr为转子电阻;s为滑差。
[0045] 根据上述定子dq轴短路电流的解析表达式可知:如果定子dq轴短路电流的解析表达式已知,则参数Lr/M、Lm/M、Ts以及Tr可知,且由于ω0与初始状态有关,数值已知。则根据Lr/M、Lm/M、Ts以及Tr可求得参数Rs和Rr。根据Lr/M和Lm/M可求得LsLr以及Lm。因此如果定子dq轴短路电流已知,则参数Rs、Rr、LsLr以及Lm可辨识,但参数Ls和Lr不可单独辨识。
[0046] 步骤2:以定子dq轴短路电流为观测量,计算各参数的轨迹灵敏度,如图3所示。从图3中参数Lσs和Lσr的轨迹灵敏度曲线,可以判断参数Lσs和Lσr无法区分辨识,由于Ls=Lσs+Lm,Lr=Lσr+Lm,则参数Ls与Lr也无法区分辨识。该结论与步骤1中的结论一致。
[0047] 步骤3:为定量比较各参数轨迹灵敏度的大小,进一步计算各参数的轨迹灵敏度的绝对平均值。选择受扰轨线时间窗口为变流器控制环节动作结束后的时段,以便忽略控制器的动态特性对短路电流的影响。本发明选择从1.06s开始到故障结束前的时间窗口进行灵敏度计算。不同观察时窗内各参数轨迹灵敏度的结果见表2。
[0048] 表2不同观察时窗下各参数的轨迹灵敏度
[0049]时窗/s Rs/pu Lσs/pu Rr/pu Lσr/pu Lm/pu
1.06-1.10 0.0992 0.2657 0.0462 0.2103 0.0924
1.06-1.12 0.1031 0.2582 0.0340 0.2095 0.1998
1.06-1.14 0.1022 0.2489 0.0272 0.1801 0.2064
1.06-1.16 0.1030 0.2711 0.0230 0.2191 0.2092
1.06-1.18 0.1025 0.2607 0.0210 0.1966 0.2164
1.06-1.20 0.1034 0.2609 0.0200 0.1931 0.2306
[0050] 步骤4:从表2中选择轨迹灵敏度绝对平均值最大的时间窗口的受扰轨线进行参数辨识。如选择时间窗口[1.06-1.2s]内的定子dq轴短路电流受扰轨线辨识定子电阻Rs与互感Lm,选择时间窗口[1.06-1.1s]内定子dq轴短路电流受扰轨线辨识转子电阻Rr,选择时间窗口[1.06-1.16s]内定子dq轴短路电流受扰轨线辨识定转子自感之积LsLr。
[0051] 步骤5:采用粒子群优化算法辨识各参数。单独进行多次(本发明采用50次)参数辨识,多次辨识结果见图4。将多次参数辨识结果平均,剔除与均值偏差大于100%的辨识结果;再将剩下的参数重新平均,作为DFIG各参数的最终辨识结果,见表3。
[0052] 表3DFIG参数辨识结果
[0053]参数 真值 均值 方差×10-4 均值偏差/%
Rs/pu 0.00706 0.00745 33.00 5.5241%
Rr/pu 0.005 0.00529 70.00 5.8000%
LsLr/pu 9.385 9.43672 14.43 0.5511%
Lm/pu 2.9 2.90209 1.64 0.0721%
[0054] 步骤6:比较基于参数辨识结果与实际短路电流的受扰轨线,见图5。
[0055] 以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭示的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
高效检索全球专利

IPRDB是专利检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,专利查询、专利分析

电话:13651749426

侵权分析

IPRDB的侵权分析产品是IPRDB结合多位一线专利维权律师和专利侵权分析师的智慧,开发出来的一款特色产品,也是市面上唯一一款帮助企业研发人员、科研工作者、专利律师、专利分析师快速定位侵权分析的产品,极大的减少了用户重复工作量,提升工作效率,降低无效或侵权分析的准入门槛。

立即试用