CN115765022A - 确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法 - Google Patents
确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法 Download PDFInfo
- Publication number
- CN115765022A CN115765022A CN202211542966.9A CN202211542966A CN115765022A CN 115765022 A CN115765022 A CN 115765022A CN 202211542966 A CN202211542966 A CN 202211542966A CN 115765022 A CN115765022 A CN 115765022A
- Authority
- CN
- China
- Prior art keywords
- doubly
- fed
- wind turbine
- fed wind
- voltage
- 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.)
- Pending
Links
Images
Classifications
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Control Of Eletrric Generators (AREA)
Abstract
本发明提供的一种确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法,包括以下步骤:S1.向双馈风机的并网点注入电压扰动小信号,并确定出并网点的扰动电压和扰动电流,并同时记录每次注入电压扰动信号时的风速vw;S2.基于并网点的扰动电压和扰动电流确定出在不同风速下的并网点的序阻抗Zsc;S3.构建包含双馈风机机侧变流器控制参数以及风速的序阻抗模型Zdfig(fp,vw);S4.将步骤S2中的序阻Zsc代入到序阻抗模型Zdfig(fp,vw)中,并采用Gauss‑Newton迭代最小二乘法求解,得到双馈风机机侧变流器控制参数。
Description
技术领域
本发明涉及一种风电参数确定方法,尤其涉及一种确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法。
背景技术
在风电领域中,对风电变流器参数辨识方法主要有时域辨识法、频域辨识法两大类,分别是基于时域并网故障/阶跃信号注入的传统识别方法以及基于系统频域响应特性的线性识别方法。
传统时域辨识方法对含有不可观测量的系统模型基本参数进行辨识,需要求解复杂的微分方程组,增加了辨识难度,且会使辨识过程不平稳,影响辨识精度。
传统频域辨识方法耗时短,能反映实际运行状况,且可用于系统稳定性分析。基于dq轴阻抗的辨识方法具有阻抗难以直接测量的缺陷,且双馈风机具有风速随机波动特性,端口阻抗特性随工况时变的特征较为明显,基于单一风速模型的参数辨识可能会大范围偏离实际参数。
因此,为了解决上述技术问题,亟需提出一种新的技术手段。
发明内容
有鉴于此,本发明的目的是提供一种确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法,基于双馈风机并网点的序阻抗进行控制参数的确定,在其中考虑了实际风速的影响,解决了阻抗数据测量工况波动限制,而且有效提升了工况适应能力,提高控制参数辨识精度,为实际风电系统的稳定性分析与控制提供准确的数据支持。
本发明提供的一种确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法,包括以下步骤:
S1.向双馈风机的并网点注入电压扰动小信号,并确定出并网点的扰动电压和扰动电流,并同时记录每次注入电压扰动信号时的风速vw;
S2.基于并网点的扰动电压和扰动电流确定出在不同风速下的并网点的序阻抗Zsc;
S3.构建包含双馈风机机侧变流器控制参数以及风速的序阻抗模型Zdfig(fp,vw);
S3.将步骤S2中的序阻Zsc代入到序阻抗模型Zdfig(fp,vw)中,并采用Gauss-Newton迭代最小二乘法求解,得到双馈风机机侧变流器控制参数,其中,控制参数为双馈风机机侧变流器控制参数,包括电流控制比例系数krp、电流控制积分系数kri,锁相环比例系数kpllp和锁相环积分系数kplli。
进一步,步骤S1中,通过如下方法确定出并网点的扰动电压和扰动电流:
其中:Ip为并网点扰动电流,ia为双馈风机的a相电流,ib为双馈风机的b相电流,Vp为并网点的扰动电压,vab为双馈风机的a相和b相之间的线电压,vbc为双馈风机的b相和c相之间的线电压。
进一步,步骤S2中,通过如下方法确定不同风速下的序阻抗Zsc:
进一步,序阻抗模型Zdfig(fp,vw)具体为:
Tpll(fp-f1,vw)=Hpll(fp-f1)/(1+Vs(vw)·Hpll(fp-f1));
Hpll(fp-f1)=(kpllp+kplli/(fp-f1))/(fp-f1);
ηp(fp)=ke/(jωpLm),μp(fp)=ke(Rs+jωpLs)/(jωpLm);
其中:fp为注入的电压扰动小信号的频率,Lm为双馈风机的励磁电感,Lr为双馈风机转子侧电感,ωr(vw)为在风速vw下双馈风机转子电气角速度,ke为双馈风机的定子与转子的绕组的匝数比;Krd(vw)为转子侧变流器电流控制的解耦项;Hri()表示双馈风机中PI控制器的电流控制传递函数,Tpll()为双馈风机的锁相环闭环传递函数,Rs为双馈风机定子侧电阻,为双馈风机转子侧电阻经过绕组折算至定子侧后的值;Dr(vw)为双馈风机转子侧变流器PWM控制信号的占空比,f1为电网的基波频率,Vdc为双馈风机直流母线的电压,Ir(vw)为双馈风机转子侧电流。
进一步,通过如下方法确定在风速vw下双馈风机转子电气角速度ωr(vw):
ωr(vw)=qωmopt(vw);其中,ωmopt(vw)为双馈风机的最有机械转速,q为风机转子的极对数;其中:ωmopt(vw)=Nωt(vw),N为双馈风机的齿轮箱的变速比;
进一步,通过如下方法确定双馈风机转子侧电流Ir(vw)以及双馈风机转子侧变流器PWM控制信号的占空比Dr(vw):
式(1)中:η=ke/(jω1Lm),μ=ke(Rs+jω1Ls)/(jω1Lm);
联系式(1)-(4)式形成方程组,从而解出Dr(vw)和Ir(vw);
其中:*表示共轭运算,Re{}表示取变量的实部,Is(vw)为双馈风机定子侧的电流,Rs为双馈风机的定子侧的电阻,Ug0表示电网电压,Vs(vw)为双馈风机定子侧电压,Ls为双馈风机定子侧的电感,ω1为电网的基频角速度,为双馈风机转子侧电压经过绕组折算至定子侧后的值,Pmopt(vw)为双馈风机的最优机械功率,Peopt(vw)为双馈风机定子侧的最优功率,A为双馈风机叶片迎风面积,ρ为空气密度。
本发明的有效果:通过本发明,基于双馈风机并网点的序阻抗进行控制参数的确定,在其中考虑了实际风速的影响,解决了阻抗数据测量工况波动限制,而且有效提升了工况适应能力,提高控制参数辨识精度,为实际风电系统的稳定性分析与控制提供准确的数据支持。
附图说明
下面结合附图和实施例对本发明作进一步描述:
图1为本发明的流程图。
图2为本发明的双馈风机并网系统示意图。
图3为本发明的序阻抗测量示意图。
具体实施方式
以下进一步对本发明做出详细说明:
在双馈风机并网系统中,如图1所示:双馈风机阻抗特性由机侧变流器和发电机阻抗主导,为评估双馈风机运行特性以及接入电网的稳定性,通常需要参数辨识获取系统的机侧控制器控制参数。双馈风机并网系统如图1所示,双馈风机机侧变流器主要包括电流控制和锁相环两个控制环节,Hri(s)和Hpll(s)分别表示转子电流控制和单同步锁相环的PI控制器传递函数;Hri(s)=krp+kri/s,krp、kri为需要辨识的电流控制比例系数和积分系数。Hpll(s)=(kpllp+kplli/s)/s,kpllp、kplli即为需要辨识的锁相环比例系数和积分系数。图1中,Vdc表示直流母线电压,vsa、vsb与vsc表示三相定子电压,isa、isb与isc表示三相定子电流,vra、vrb与vrc表示三相转子绕组电压,ira、irb与irc表示三相转子绕组电流。θpll表示锁相环输出相位,θr表示转子电压相位。ωs、ωr分别表示转子侧电气旋转速度和机械旋转速度。Ir_ref、Ir分别表示转子电流相量参考值和测量值,Krd表示矢量控制解耦系数,dra、drb、drc表示变流器占空比。
本发明提供的一种确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法,包括以下步骤:
S1.向双馈风机的并网点注入电压扰动小信号,并确定出并网点的扰动电压和扰动电流,并同时记录每次注入电压扰动信号时的风速vw;其中,采用现有的频率分析仪(frequency response analyzer,FRA)向并网点注入电压扰动小信号;
S2.基于并网点的扰动电压和扰动电流确定出在不同风速下的并网点的序阻抗Zsc;
S3.构建包含双馈风机机侧变流器控制参数以及风速的序阻抗模型Zdfig(fp,vw);
S3.将步骤S2中的序阻Zsc代入到序阻抗模型Zdfig(fp,vw)中,并采用Gauss-Newton迭代最小二乘法求解,得到双馈风机机侧变流器控制参数,其中,控制参数为双馈风机机侧变流器控制参数,包括电流控制比例系数krp、电流控制积分系数kri,锁相环比例系数kpllp和锁相环积分系数kplli,通过上述方法,基于双馈风机并网点的序阻抗进行控制参数的确定,在其中考虑了实际风速的影响,解决了阻抗数据测量工况波动限制,而且有效提升了工况适应能力,提高控制参数辨识精度,为实际风电系统的稳定性分析与控制提供准确的数据支持。
本实施例中,步骤S1中,通过如下方法确定出并网点的扰动电压和扰动电流:
其中:Ip为并网点扰动电流,ia为双馈风机的a相电流,ib为双馈风机的b相电流,Vp为并网点的扰动电压,vab为双馈风机的a相和b相之间的线电压,vbc为双馈风机的b相和c相之间的线电压。
其中:步骤S2中,通过如下方法确定不同风速下的序阻抗Zsc:
本实施例中,序阻抗模型Zdfig(fp,vw)具体为:
Tpll(fp-f1,vw)=Hpll(fp-f1)/(1+Vs(vw)·Hpll(fp-f1));
Hpll(fp-f1)=(kpllp+kplli/(fp-f1))/(fp-f1);
ηp(fp)=ke/(jωpLm),μp(fp)=ke(Rs+jωpLs)/(jωpLm);
其中:fp为注入的电压扰动小信号的频率,Lm为双馈风机的励磁电感,Lr为双馈风机转子侧电感,ωr(vw)为在风速vw下双馈风机转子电气角速度,ke为双馈风机的定子与转子的绕组的匝数比;Krd(vw)为转子侧变流器电流控制的解耦项;Hri()表示双馈风机中PI控制器的电流控制传递函数,Tpll()为双馈风机的锁相环闭环传递函数,Rs为双馈风机定子侧电阻,Rr s为双馈风机转子侧电阻经过绕组折算至定子侧后的值;Dr(vw)为双馈风机转子侧变流器PWM控制信号的占空比,f1为电网的基波频率,Vdc为双馈风机直流母线的电压,Ir(vw)为双馈风机转子侧电流。
上述模型通过下属的过程建立:
构建序阻抗计算公式:
在公式(2)中加入扰动小信号后,通过线性化处理(处理过程为现有)的模型如下:
其中fs=fp-fr,表示转子电压电流扰动频率,fr为转子电气旋转频率;考虑图2中RSC(双馈风机变流器)电流内环控制的影响,可以将占空比信号Drp1(fp)表示为:
扰动电压会使锁相环输出角度中引入角度扰动,进而使得经过Park变换的变量额外增加扰动分量;其中,变流器锁相环的影响表示如下:
本实施例中,通过如下方法确定在风速vw下双馈风机转子电气角速度ωr(vw):
ωr(vw)=qωmopt(vw);其中,ωmopt(vw)为双馈风机的最有机械转速,q为风机转子的极对数;其中:ωmopt(vw)=Nωt(vw),N为双馈风机的齿轮箱的变速比;
本实施例中,通过如下方法确定双馈风机转子侧电流Ir(vw)以及双馈风机转子侧变流器PWM控制信号的占空比Dr(vw):
在建立公式(1)至(5)之前,需要确定双馈风机的稳态方程:
其中:定子侧电压电流Vs(vw)、Is(vw)和转子侧电压电流Vr(vw)、Ir(vw)均为关于风速的函数,ωs为转差角速度,Lls为双馈风机的定子绕组的漏感,为转子绕组的漏感通过绕组折算至定子侧后的值;
式(1)中:η=ke/(jω1Lm),μ=ke(Rs+jω1Ls)/(jω1Lm);
Vs(vw)=Ug0 (5);
联系式(1)-(4)式形成方程组,从而解出Dr(vw)和Ir(vw);
其中:*表示共轭运算,Re{}表示取变量的实部,Is(vw)为双馈风机定子侧的电流,Rs为双馈风机的定子侧的电阻,Ug0表示电网电压,Vs(vw)为双馈风机定子侧电压,Ls为双馈风机定子侧的电感,ω1为电网的基频角速度,为双馈风机转子侧电压经过绕组折算至定子侧后的值,Pmopt(vw)为双馈风机的最优机械功率,Peopt(vw)为双馈风机定子侧的最优功率,A为双馈风机叶片迎风面积,ρ为空气密度。
上述中将一个参数经过绕组折算至定子侧的值为将该参数的值除以定子和转子的匝数比来确定。
本实施例中,Gauss-Newton迭代最小二乘法求解过程为一个现有技术,下面简要说明下计算过程:
将序阻抗模型转换成用于参数辨识的模型函数,
f(x1,x2;β1,β2,β3,β4)=20Lg[Zdfig(fp,vw;krp,kri,kpllp,kplli)] (10)
其中,向量x1对应序阻抗模型中的扰动频率,向量β对应变流器控制器参数。值得注意的是,本专利引入了风速变量vw,对应模型函数中的变量x2。假设有m个观测点,建立优化目标函数为,
yi表示第i个观测点函数值。第i次观测点预测偏差记为,
则偏差向量可以表示为r=[r1…rm]T,则目标函数可以进一步表示为:
对目标函数求偏导,得到目标函数的梯度,
可以推导得到偏差向量函数的雅可比矩阵为
进一步推导参数优化递推公式,
式中(k)和(k+1)表示当前的参数和将优化的参数,考虑到预测误差ri趋于0的情况,可以得到求解控制器参数的最终递推公式,
式中Δβ(k+1)表示下次迭代控制参数增量。根据序阻抗模型,推导式(18)所示的Gauss-Newton回归递推公式,然后对控制参数进行递推求解,直到满足精度要求。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法,其特征在于:包括以下步骤:
S1.向双馈风机的并网点注入电压扰动小信号,并确定出并网点的扰动电压和扰动电流,并同时记录每次注入电压扰动信号时的风速vw;
S2.基于并网点的扰动电压和扰动电流确定出在不同风速下的并网点的序阻抗Zsc;
S3.构建包含双馈风机机侧变流器控制参数以及风速的序阻抗模型Zdfig(fp,vw);
S4.将步骤S2中的序阻Zsc代入到序阻抗模型Zdfig(fp,vw)中,并采用Gauss-Newton迭代最小二乘法求解,得到双馈风机机侧变流器控制参数,其中,控制参数为双馈风机机侧变流器控制参数,包括电流控制比例系数krp、电流控制积分系数kri,锁相环比例系数kpllp和锁相环积分系数kplli。
4.根据权利要求1所述确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法,其特征在于:序阻抗模型Zdfig(fp,vw)具体为:
Tpll(fp-f1,vw)=Hpll(fp-f1)/(1+Vs(vw)·Hpll(fp-f1));
Hpll(fp-f1)=(kpllp+kplli/(fp-f1))/(fp-f1);
ηp(fp)=ke/(jωpLm),μp(fp)=ke(Rs+jωpLs)/(jωpLm);
其中:fp为注入的电压扰动小信号的频率,Lm为双馈风机的励磁电感,Lr为双馈风机转子侧电感,ωr(vw)为在风速vw下双馈风机转子电气角速度,ke为双馈风机的定子与转子的绕组的匝数比;Krd(vw)为转子侧变流器电流控制的解耦项;Hri()表示双馈风机中PI控制器的电流控制传递函数,Tpll()为双馈风机的锁相环闭环传递函数,Rs为双馈风机定子侧电阻,Rr s为双馈风机转子侧电阻经过绕组折算至定子侧后的值;Dr(vw)为双馈风机转子侧变流器PWM控制信号的占空比,f1为电网的基波频率,Vdc为双馈风机直流母线的电压,Ir(vw)为双馈风机转子侧电流。
6.根据权利要求4所述确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法,其特征在于:通过如下方法确定双馈风机转子侧电流Ir(vw)以及双馈风机转子侧变流器PWM控制信号的占空比Dr(vw):
式(1)中:η=ke/(jω1Lm),μ=ke(Rs+jω1Ls)/(jω1Lm);
Vs(vw)=Ug0 (5);
联系式(1)-(4)式形成方程组,从而解出Dr(vw)和Ir(vw);
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211542966.9A CN115765022A (zh) | 2022-12-02 | 2022-12-02 | 确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211542966.9A CN115765022A (zh) | 2022-12-02 | 2022-12-02 | 确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115765022A true CN115765022A (zh) | 2023-03-07 |
Family
ID=85343383
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211542966.9A Pending CN115765022A (zh) | 2022-12-02 | 2022-12-02 | 确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115765022A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117277412A (zh) * | 2023-09-26 | 2023-12-22 | 中国电力科学研究院有限公司 | 风电机组变流器的控制参数识别方法、系统、设备和介质 |
-
2022
- 2022-12-02 CN CN202211542966.9A patent/CN115765022A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117277412A (zh) * | 2023-09-26 | 2023-12-22 | 中国电力科学研究院有限公司 | 风电机组变流器的控制参数识别方法、系统、设备和介质 |
CN117277412B (zh) * | 2023-09-26 | 2024-02-27 | 中国电力科学研究院有限公司 | 风电机组变流器的控制参数识别方法、系统、设备和介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111525611B (zh) | 计及频率耦合效应的双馈并网系统次同步振荡分析方法 | |
CN112165115B (zh) | 一种直驱风电机组灰箱模型的参数辨识方法及装置 | |
CN109327171B (zh) | 一种适用于轨道交通牵引电机参数在线辨识的方法 | |
CN112907075B (zh) | 一种电力系统综合负荷模型参数辨识方法 | |
CN110071667B (zh) | 一种双馈风力发电机参数辨识方法 | |
CN103091109B (zh) | 用于风力发电机动态特性测试的风轮机模拟器的控制方法 | |
CN115765022A (zh) | 确定基于序阻抗特性的双馈风机机侧变流器控制参数的方法 | |
CN112615393B (zh) | 基于矢量拟合的直驱风电机组控制器参数辨识方法及装置 | |
CN109873587B (zh) | 一种永磁同步电机多参数自动辨识方法 | |
CN109921421B (zh) | 双馈风电机组输出谐波电流模型建立方法 | |
CN109782173A (zh) | 异步电机励磁互感曲线测量系统及其测量方法 | |
CN110212575A (zh) | 双馈风机单机无穷大系统的小干扰稳定性模态分析方法 | |
Ruan et al. | Magnetizing curve estimation of induction motors in single-phase magnetization mode considering differential inductance effect | |
CN115392110A (zh) | 基于ppo算法的数据模型混合驱动风电场建模方法 | |
CN111092579B (zh) | 带有定子温度在线监测的异步电机自适应矢量控制系统 | |
CN106788056B (zh) | 电机定子电阻的在线辨识方法、装置和电机控制系统 | |
WO2024120487A1 (zh) | 一种永磁同步电机性能的标定方法 | |
CN110601175A (zh) | 一种传统xad基值体系下励磁电流改进标幺方法 | |
CN113378347B (zh) | 一种基于模块化多端口的风电机组频域阻抗建模方法 | |
Pan et al. | Hierarchical parameter estimation of DFIG and drive train system in a wind turbine generator | |
CN110417047B (zh) | 基于复转矩系数分析双馈风机ssci阻尼特性的方法 | |
CN102081691A (zh) | 考虑频率变化的电动机实用化建模方法 | |
CN111293693A (zh) | 基于扩展卡尔曼滤波的双馈风机变流器控制参数辨识方法 | |
CN112994113B (zh) | 一种双馈风力发电机序阻抗建模方法及系统 | |
CN110412353B (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 |