CN112994111A - 考虑功率控制的永磁同步风机接入弱电网稳定性分析方法 - Google Patents

考虑功率控制的永磁同步风机接入弱电网稳定性分析方法 Download PDF

Info

Publication number
CN112994111A
CN112994111A CN202110451859.4A CN202110451859A CN112994111A CN 112994111 A CN112994111 A CN 112994111A CN 202110451859 A CN202110451859 A CN 202110451859A CN 112994111 A CN112994111 A CN 112994111A
Authority
CN
China
Prior art keywords
side converter
grid
generator
converter controller
power
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.)
Granted
Application number
CN202110451859.4A
Other languages
English (en)
Other versions
CN112994111B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202110451859.4A priority Critical patent/CN112994111B/zh
Publication of CN112994111A publication Critical patent/CN112994111A/zh
Application granted granted Critical
Publication of CN112994111B publication Critical patent/CN112994111B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/241—The oscillation concerning frequency
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Power Engineering (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明涉及一种考虑功率控制影响的永磁同步风机接入弱电网稳定性分析方法。新能源发电多采用扰动观察法进行最大功率跟踪，其中的非线性不连续环节导致稳定性分析困难。本发明基于描述函数法分析了永磁同步风机接入弱电网***的稳定性，充分考虑了功率环中的非线性不连续环节，使分析结果更加精确。同时，描述函数法是一种可以量化地计算振荡功率和幅值的方法，本发明所述分析方法可为振荡抑制及控制器设计提供有力良好参照。

Description

考虑功率控制的永磁同步风机接入弱电网稳定性分析方法
技术领域
本发明涉及一种考虑功率控制影响的永磁同步风机接入弱电网稳定性分析方法，属于电力***中新能源发电***稳定性分析领域。
背景技术
可再生能源发电占比不断提升，风力发电装机容量持续增长，风机并网***的稳定运行成为重要问题。永磁同步风电机组(PMSG)具有能量转换效率高、运行可靠性好、可控性强的优势，装机占比持续提高。为了更好地利用风能，必须采用合适的最大功率跟踪算法，其中扰动观察法无需测量风速，具有较好的实用价值。然而，扰动观察法隐含非线性不连续环节，导致传统的小信号稳定性分析方法难以适用，常规的稳定性分析方法一般忽视基于扰动观察法的功率控制环。描述函数法可以很好地完成非线性环节建模和***的稳定性分析，其基本思想是，当***满足一定假设条件时，在正弦输入作用下，***非线性环节的输出可用一次谐波分量来近似表示，从而获得非线性环节的近似等效频率特性。因此，本发明基于描述函数法提出一种考虑功率控制影响的永磁同步风机接入弱电网稳定性分析方法，充分考虑功率环中的非线性环节，提高稳定性分析的精确性。
发明内容
本发明旨在提出一种考虑功率控制影响的永磁同步风机接入弱电网稳定性分析方法。所述永磁同步风机包括风力机、发电机、机侧变流器、直流电容、网侧变流器、滤波器、机侧变流器控制器和网侧变流器控制器。机侧变流器控制器中，包含功率环，转速环，机侧电流环；网侧变流器控制器中，包含电压环和网侧电流环。
本发明采用如下技术方案：
步骤一：获取永磁同步风机的主要参数，分别建立风力机、发电机及机侧变流器、机侧变流器控制器、直流电容、网侧变流器及滤波器、网侧变流器控制器的数学模型，并在dq坐标系下进行线性化，计算稳态运行参数，得到风力机、发电机及机侧变流器、机侧变流器控制器、直流电容、网侧变流器及滤波器、网侧变流器控制器的小信号模型。
步骤二：基于描述函数法，对机侧变流器控制器中的功率环进行建模，其数学表达式为
Figure BDA0003039022330000011
其中，ε为功率环扰动步长，Tp为功率环控制周期，Pref为永磁同步风机输出功率参考值，P表示永磁同步风机输出功率，Pn为当前采样时刻n的永磁同步风机输出功率，Pn-1为上一采样时刻永磁同步风机输出功率，ωg表示发电机转速，
Figure BDA0003039022330000025
表示发电机转速的参考值，ωg，n为当前采样时刻发电机转速，ωg，n-1为上一采样时刻发电机转速。sgn(x)为符号函数，当x≥0时，sgn(x)＝1，当x＜0时，sgn(x)＝-1。考虑实际永磁同步风机的功率-转速曲线，可以认为
Figure BDA0003039022330000021
其中，ωmpp表示最大功率点处的发电机转速。进而功率环模型可简化为
Figure BDA0003039022330000022
式中的符号函数可以采用描述函数进行建模，其描述函数为
Figure BDA0003039022330000023
式中A表示输入信号的幅值。
步骤三：考虑交流弱电网影响，将弱电网、功率环模小信号模型的线性部分与步骤一中建立的永磁同步风机小信号模型联立，推导***线性部分的传递函数G(s)。
步骤四：在复平面中绘制G(s)与-1/N(A)曲线，基于描述函数法分析***稳定性，具体方法为，若G(s)包含右半平面极点，则***必不稳定；若G(s)不包含右半平面极点，则通过G(s)轨迹与-1/N(A)轨迹的关系判断***稳定性：
a、若G(s)曲线不包围-1/N(A)曲线，则***是稳定的，不发生振荡；
b、若G(s)曲线与-1/N(A)曲线相交，则***是临界稳定的，此时***发生恒幅恒频振荡，可以通过下式计算振荡的频率和幅值
Figure BDA0003039022330000024
其中，G(jω)＝GRe(ω)+jGIm(ω)，GRe表示求取复数实部，GIm表示求取复数虚部，ω0为振荡角频率，A0为振荡幅值；
c、若G(s)曲线包围-1/N(A)曲线，则***是不稳定的。
进一步地，步骤一中，风力机、发电机及机侧变流器、机侧变流器控制器的小信号模型建立如下：
建立风力机的数学模型为
sJωg＝Tm-Te-Bωg
式中，J为风力机和发电机的等效集中质量块转动惯量，Tm为发电机机械转矩，Te为发电机电磁转矩，B为自阻尼系数，此处认为B＝0，s为拉普拉斯变化引入的参变量。对该模型进行线性化可得
sJΔωg＝ΔTm-ΔTe
发电机的电磁转矩为
Figure BDA0003039022330000031
np为发电机极对数，iqr为q轴发电机定子电流，ψf为发电机永磁体磁链。对该式进行线性化可得
Figure BDA0003039022330000032
发电机的机械转矩为
Tm＝Btωg
式中，
Figure BDA0003039022330000033
为风力机机械转矩线性化常数，ωg *为发电机转速稳态值，在平衡点处为特定常数。对该式线性化可得
ΔTm＝BtΔωg
由此可得风力机小信号模型，
Figure BDA0003039022330000034
令
Figure BDA0003039022330000035
则Δωg＝Giqω·Δidqr，Δωe＝npGiqω·Δidqr。
发电机及机侧变流器的数学模型为
Figure BDA0003039022330000036
式中，Rs、Ls分别为发电机转子电阻和电枢电感，ωe为转子电角速度，ωe＝npωg。idr、iqr为dq坐标系下发电机定子电流，ddr、dqr为dq坐标系下机侧变流器控制器输出占空比，udc为直流电压。对该模型进行线性化，可得发电机及机侧变流器的小信号模型为
Figure BDA0003039022330000037
其中，
Figure BDA0003039022330000038
Figure BDA0003039022330000039
大写字母及上标*表示对应小写变量的稳态分量，Δ表示对应变量的小信号分量。
建立机侧变流器控制器的数学模型为
Figure BDA0003039022330000041
其中，
Figure BDA0003039022330000042
Kcpr、Kcir分别为机侧电流环PI控制的比例参数和积分参数，
Figure BDA0003039022330000043
Kωp、Kωi分别为转速环PI控制的比例参数和积分参数。上标c表示机侧变流器控制器dq坐标系。
Figure BDA0003039022330000044
为发电机转速的参考值。对其进行线性化，可得机侧变流器控制器的小信号模型为
Figure BDA0003039022330000045
其中，
Figure BDA0003039022330000046
Figure BDA00030390223300000412
为转子电角速度稳态值。受转速扰动影响，动态过程中机侧变流器控制器的dq坐标系将与发电机dq坐标系存在相角差。发电机电角度为
θe＝npωg/s
对其进行线性化，可得
Figure BDA0003039022330000047
式中，
Figure BDA0003039022330000048
由此，机侧变流器控制器dq坐标系与发电机dq坐标系之间变量的转换关系为
Figure BDA0003039022330000049
其中，
Figure BDA00030390223300000410
则机侧变流器控制器的小信号模型为
Figure BDA00030390223300000411
进一步地，步骤一中，小信号模型建立过程如下：
建立直流电容的数学模型为
sCdcudc＝idc2-idc1＝1.5(ddridr+dqriqr)-1.5(ddgidg+dqgiqg)
式中，Cdc为直流电容，idc1为网侧直流电流，idc2为机侧直流电流，idg、iqg分别为网侧变流器交流端口d、q轴电流，ddg、dqg为dq坐标系下网侧变流器控制器输出占空比。udc为直流电压，idr、iqr分别为d、q坐标系发电机定子电流，ddr、dqr为dq坐标系下机侧变流器控制器输出占空比，对其进行线性化，可得直流电容小信号模型为
Figure BDA0003039022330000051
式中，
Figure BDA0003039022330000052
大写字母及上标*表示对应小写变量的稳态分量，Δ表示对应变量的小信号分量。
进一步地，步骤一中，所述网侧变流器、网侧变流器控制器的小信号模型建立过程如下：
建立网侧变流器及滤波器的数学模型为
Figure BDA0003039022330000053
式中，Lf为滤波电感，ω为工频角频率，ω＝100πrad/s，idg、iqg分别为网侧变流器交流端口d、q轴电流，ddg、dqg为dq坐标系下网侧变流器控制器输出占空比，udc为直流电压，udg、uqg分别为并网点d、q轴电压。对该模型进行线性化，可得网侧变流器及滤波器的小信号模型为
Figure BDA0003039022330000054
其中，
Figure BDA0003039022330000055
大写字母及上标*表示对应小写变量的稳态分量，Δ表示对应变量的小信号分量。
建立网侧变流器控制器的数学模型为
Figure BDA0003039022330000056
其中，
Figure BDA0003039022330000057
Kcpg、Kcig分别为网侧电流环PI控制的比例参数和积分参数，
Figure BDA0003039022330000058
Kvp、Kvi分别为网侧电压环PI控制的比例参数和积分参数，Udcref为直流电压参考值。在网侧变流器控制器中，采用锁相环使风机与电网保持同步。上标c表示网侧变流器控制器dq坐标系。对该模型进行线性化，可得网侧变流器控制器的小信号模型为
Figure BDA0003039022330000059
其中，
Figure BDA0003039022330000061
此外，网侧变流器中还应考虑锁相环动态，其数学模型为
Figure BDA0003039022330000062
其中，
Figure BDA0003039022330000063
Kppll、Kipll分别为锁相环PI控制器的比例参数和积分参数，
Figure BDA0003039022330000064
为网侧变流器控制器dq坐标系下的并网点q轴电压。对其进行线性化，可得
Figure BDA0003039022330000065
其中，***dq坐标系与控制器dq坐标系存在一定偏差，二者之间可以通过如下方程相互转换
Figure BDA0003039022330000066
式中变量Δxd、Δxq可以表示网侧变流器输出电流Δidg、Δiqg、并网点电压Δudg、Δuqg或网侧控制器输出占空比Δddg、Δdqg，
Figure BDA0003039022330000067
表示对应变量的稳态分量
Figure BDA0003039022330000068
Figure BDA0003039022330000069
由此可推导出锁相环小信号模型，即
Δθ＝Gpll·Δuqg
其中，
Figure BDA00030390223300000610
从而可以得出控制器dq坐标系与***dq坐标系之间的关系为
Figure BDA00030390223300000611
其中，
Figure BDA00030390223300000612
则网侧变流器控制器的小信号模型转变为
Figure BDA00030390223300000613
进一步地，步骤三具体为：
交流弱电网采用理想电压源串联等效电感表示，建立其数学模型为
Figure BDA00030390223300000614
式中，Lg为弱电网等效电感，uds、uqs分别为d、q轴理想电压源电压，idg、iqg分别为网侧变流器交流端口d、q轴电流。将该式线性化，可得
Zg·Δidqg＝Δudqg
式中，
Figure BDA0003039022330000071
永磁同步风机的输出功率为
P＝1.5(idgudg+iqguqg)
将该式线性化可得
Figure BDA0003039022330000072
式中，
Figure BDA0003039022330000073
结合前述永磁同步风机并网***各部分模型的线型部分，可得***线性部分传递函数G(s)为
Figure BDA0003039022330000074
Figure BDA0003039022330000075
Figure BDA0003039022330000076
Figure BDA0003039022330000077
Figure BDA0003039022330000078
Figure BDA0003039022330000079
Figure BDA00030390223300000710
Figure BDA00030390223300000711
Figure BDA00030390223300000712
式中，Tf表示功率采样滤波器周期，1/(1+Tfs)为功率采样滤波器延时，1/(1+1.5Tps)为控制器及PWM延时。
与现有技术相比，本发明的优点在于：
(1)本发明将描述函数法应用于永磁同步风机并网***稳定性分析中，充分考虑基于扰动观察法的功率环中的非线性环节，克服了传统小信号分析法无法适用于不连续、非线性环节的问题。
(2)当***处于临界稳定状态时，发生恒幅恒频振荡。描述函数法可以量化地计算振荡频率和幅值，为振荡预防和抑制提供了重要依据。
附图说明
图1为永磁同步风机拓扑结构(A)及其控制器结构(B、D为坐标变换，C为锁相环，E为机侧变流器控制器，F为网侧变流器控制器)
图2为G(s)极点分布图
图3为改变Lg值时G(s)与-1/N(A)曲线
图4为Lg＝0.1mH时的仿真验证波形(a)与FFT频谱分析(b)
图5为Lg＝0.4mH时的仿真验证波形(a)与FFT频谱分析(b)
具体实施方式
本发明所述永磁同步风机拓扑结构及其控制器如图1所示，包含风力机、发电机、机侧变流器、直流电容、网侧变流器、滤波器。风力机捕捉风能并转化为机械能，机侧变流器和网侧变流器均为两电平电压源型变换器，其中，机侧变流器将永磁同步发电机输出的交流电转换成直流电，网侧变流器将直流电逆变成工频交流电并入电网。机侧变流器和网侧变流器均采用dq坐标系下的矢量控制方法，通过dq变换将abc三相电压电流转换为d轴和q轴下的电压电流。机侧变流器控制器中，包含功率环(P&O)，转速环(Hω)，机侧电流环(Hcr)；网侧变流器控制器中，包含电压环(Hv)和网侧电流环(Hcg)。此外，还包含锁相环(PLL)和坐标变换环节(abc/dq和dq/abc)。下面，结合具体实施例对本发明作进一步说明：
本发明的一个实施例中，***的主要参数如表1所示。
表1***主要参数
Figure BDA0003039022330000081
Figure BDA0003039022330000099
在本发明的实施例中，第一步，获取永磁同步风机的主要参数如表1所示，分别建立风力机、发电机及机侧变流器、机侧变流器控制器、直流电容、网侧变流器及滤波器、网侧变流器控制器的数学模型，并在dq坐标系下进行线性化，计算稳态运行参数，得到风力机、发电机及机侧变流器、机侧变流器控制器、直流电容、网侧变流器及滤波器、网侧变流器控制器的小信号模型：
sJΔωg＝ΔTm-ΔTe
Figure BDA0003039022330000091
ΔTm＝BtΔωg
Figure BDA0003039022330000092
Figure BDA0003039022330000093
Figure BDA0003039022330000094
Figure BDA0003039022330000095
Figure BDA0003039022330000096
Figure BDA0003039022330000097
Figure BDA0003039022330000098
s为拉普拉斯变化引入的参变量，J为风力机和发电机的等效集中质量块转动惯量，ωg表示发电机转速，Tm为发电机机械转矩，Te为发电机电磁转矩。np为发电机极对数，ψf为发电机永磁体磁链，Bt为风力机机械转矩线性化常数。
Figure BDA0003039022330000101
Figure BDA0003039022330000102
Rs、Ls分别为发电机转子电阻和电枢电感，ωe为转子电角速度，ωe＝npωg。idr、iqr为dq坐标系下发电机定子电流，ddr、dqr为dq坐标系下机侧变流器控制器输出占空比，udc为直流电压。
Figure BDA0003039022330000103
Figure BDA0003039022330000104
Kcpr、Kcir分别为机侧电流环PI控制的比例参数和积分参数，
Figure BDA0003039022330000105
Kωp、Kωi分别为转速环PI控制的比例参数和积分参数，
Figure BDA0003039022330000106
为发电机转速的参考值，
Figure BDA0003039022330000107
为转子电角速度稳态值，上标c表示变流器控制器dq坐标系。
Figure BDA0003039022330000108
Figure BDA0003039022330000109
Figure BDA00030390223300001010
Cdc为直流电容，idg、iqg分别为网侧变流器交流端口d、q轴电流，ddg、dqg为dq坐标系下网侧变流器控制器输出占空比。
Figure BDA00030390223300001011
Figure BDA00030390223300001012
Lf为滤波电感，ω为工频角频率，ω＝100πrad/s，idg、iqg分别为网侧变流器交流端口d、q轴电流，ddg、dqg为dq坐标系下网侧变流器控制器输出占空比，udg、uqg分别为并网点d、q轴电压。
Figure BDA00030390223300001013
Figure BDA00030390223300001014
Kcpg、Kcig分别为网侧电流环PI控制的比例参数和积分参数，
Figure BDA00030390223300001015
Figure BDA00030390223300001016
Kvp、Kvi分别为网侧电压环PI控制的比例参数和积分参数。
Figure BDA00030390223300001017
Figure BDA00030390223300001018
Kppll、Kipll分别为锁相环PI控制器的比例参数和积分参数。
大写字母及上标*表示对应变量的稳态分量，具体的计算方式如下。
Figure BDA0003039022330000111
第二步，基于描述函数法，对机侧变流器控制器中的功率环进行建模，其数学模型为
Figure BDA0003039022330000112
式中的符号函数可以采用描述函数进行建模，其描述函数为
Figure BDA0003039022330000113
第三步，考虑交流弱电网影响，弱电网模型为
Zg·Δidqg＝Δudqg
式中，
Figure BDA0003039022330000114
永磁同步风机的输出功率小信号模型为
Figure BDA0003039022330000115
将弱电网模型、功率环模型的线性部分与步骤一中的永磁同步风机小信号模型联立，推导***线性部分的传递函数G(s)为
Figure BDA0003039022330000116
其中，
Figure BDA0003039022330000117
Figure BDA0003039022330000118
Figure BDA0003039022330000119
Figure BDA00030390223300001110
Figure BDA00030390223300001111
Figure BDA0003039022330000121
Figure BDA0003039022330000122
Figure BDA0003039022330000123
第四步，分析***稳定性。首先绘G(s)极点图，如图2所示，可以看出G(s)不包含右半平面(实部大于0)极点，因此满足***稳定的第一个条件。在复平面中绘制G(s)与-1/N(A)图像，如图3所示，G(s)与-1/N(A)相交，说明此时***处于临界稳定状态。通过计算可知，当Lg＝0.1mH时，***的振荡频率约为129rad/s(20.5Hz)，振荡幅值约为61kW；当Lg＝0.4mH时，***的振荡频率约为131rad/s(20.9Hz)，振荡幅值约为56kW。当电网强度降低(Lg增大)时，***振荡幅值减小，说明在特定条件下，电网等效阻抗的增加有利于***保持稳定。
图4、图5分别为Lg＝0.1mH、0.4mH时***仿真及FFT频谱分析结果。图4中，直流分量为0.6889MW，对应振荡频率为19.5Hz，幅值为58kW，与理论分析结果基本一致；图5中，直流分量为0.6889MW，对应振荡频率为19.5Hz，幅值为53kW，与理论分析结果基本一致。仿真结果验证了所述分析方法的有效性和精确性。

Claims (5)

1.一种考虑功率控制影响的永磁同步风机接入弱电网稳定性分析方法，其特征在于，所述永磁同步风机包括风力机、发电机、机侧变流器、直流电容、网侧变流器、滤波器、机侧变流器控制器和网侧变流器控制器。机侧变流器控制器中，包含功率环，转速环，机侧电流环；网侧变流器控制器中，包含电压环和网侧电流环。该方法包括以下步骤：
步骤一：获取永磁同步风机的主要参数，分别建立风力机、发电机及机侧变流器、机侧变流器控制器、直流电容、网侧变流器及滤波器、网侧变流器控制器的数学模型，并在dq坐标系下进行线性化，计算稳态运行参数，得到风力机、发电机及机侧变流器、机侧变流器控制器、直流电容、网侧变流器及滤波器、网侧变流器控制器的小信号模型。
步骤二：基于描述函数法，对机侧变流器控制器中的功率环进行建模，其数学表达式为
Figure FDA0003039022320000011
其中，ε为功率环扰动步长，Tp为功率环控制周期，Pref为永磁同步风机输出功率参考值，P表示永磁同步风机输出功率，Pn为当前采样时刻n的永磁同步风机输出功率，Pn-1为上一采样时刻永磁同步风机输出功率，ωg表示发电机转速，
Figure FDA0003039022320000012
表示发电机转速的参考值，ωg，n为当前采样时刻发电机转速，ωg，n-1为上一采样时刻发电机转速。sgn(x)为符号函数，当x≥0时，sgn(x)＝1，当x＜0时，sgn(x)＝-1。考虑实际永磁同步风机的功率-转速曲线，则
Figure FDA0003039022320000013
其中，ωmpp表示最大功率点处的发电机转速。进而功率环模型可简化为
Figure FDA0003039022320000014
式中的符号函数可以采用描述函数进行建模，其描述函数为
Figure FDA0003039022320000015
式中A表示输入信号的幅值。
步骤三：考虑交流弱电网影响，将弱电网、功率环小信号模型的线性部分与步骤一中建立的小信号模型联立，推导***线性部分的传递函数G(s)。
步骤四：在复平面中绘制G(s)与-1/N(A)曲线，基于描述函数法分析***稳定性，具体方法为，若G(s)包含右半平面极点，则***必不稳定；若G(s)不包含右半平面极点，则通过G(s)轨迹与-1/N(A)轨迹的关系判断***稳定性：
a、若G(s)曲线不包围-1/N(A)曲线，则***是稳定的，不发生振荡；
b、若G(s)曲线与-1/N(A)曲线相交，则***是临界稳定的，此时***发生恒幅恒频振荡，可以通过下式计算振荡的频率和幅值
Figure FDA0003039022320000021
其中，G(jω)＝GRe(ω)+jGIm(ω)，GRe表示求取复数实部，GIm表示求取复数虚部，ω0为振荡角频率，A0为振荡幅值；
c、若G(s)曲线包围-1/N(A)曲线，则***是不稳定的。
2.如权利要求1所述方法，其特征在于，步骤一中，风力机、发电机及机侧变流器、机侧变流器控制器的小信号模型建立如下：
建立风力机的数学模型为
sJωg＝Tm-Te-Bωg
式中，J为风力机和发电机的等效集中质量块转动惯量，Tm为发电机机械转矩，Te为发电机电磁转矩，B为自阻尼系数，此处认为B＝0，s为拉普拉斯变化引入的参变量。对该模型进行线性化可得
sJΔωg＝ΔTm-ΔTe
发电机的电磁转矩为
Figure FDA0003039022320000022
np为发电机极对数，iqr为q轴发电机定子电流，ψf为发电机永磁体磁链。对该式进行线性化可得
Figure FDA0003039022320000023
发电机的机械转矩为
Tm＝Btωg
式中，
Figure FDA0003039022320000024
为风力机机械转矩线性化常数，ωg *为发电机转速稳态值，在平衡点处为特定常数。对该式线性化可得
ΔTm＝BtΔωg
由此可得风力机小信号模型，
Figure FDA0003039022320000031
令
Figure FDA0003039022320000032
则Δωg＝Giqω·Δidqr，Δωe＝npGiqω·Δidqr。
建立发电机及机侧变流器的数学模型为
Figure FDA0003039022320000033
式中，Rs、Ls分别为发电机转子电阻和电枢电感，ωe为转子电角速度，ωe＝npωg。idr、iqr为dq坐标系下发电机定子电流，ddr、dqr为dq坐标系下机侧变流器控制器输出占空比，udc为直流电压。对该模型进行线性化，可得发电机及机侧变流器的小信号模型为
Figure FDA0003039022320000034
其中，
Figure FDA0003039022320000035
Figure FDA0003039022320000036
大写字母及上标*表示对应小写变量的稳态分量，Δ表示对应变量的小信号分量。
建立机侧变流器控制器的数学模型为
Figure FDA0003039022320000037
其中，
Figure FDA0003039022320000038
Kcpr、Kcir分别为机侧电流环PI控制的比例参数和积分参数，
Figure FDA0003039022320000039
Kωp、Kωi分别为转速环PI控制的比例参数和积分参数。上标c表示机侧变流器控制器dq坐标系。
Figure FDA00030390223200000310
为发电机转速的参考值。对其进行线性化，可得机侧变流器控制器的小信号模型为
Figure FDA00030390223200000311
其中，
Figure FDA00030390223200000312
Figure FDA00030390223200000313
为转子电角速度稳态值。受转速扰动影响，动态过程中机侧变流器控制器的dq坐标系将与发电机dq坐标系存在相角差。发电机电角度为
θe＝npωg/s
对其进行线性化，可得
Figure FDA0003039022320000041
式中，
Figure FDA0003039022320000042
由此，机侧变流器控制器dq坐标系与发电机dq坐标系之间变量的转换关系为
Figure FDA0003039022320000043
其中，
Figure FDA0003039022320000044
则机侧变流器控制器的小信号模型为
Figure FDA0003039022320000045
3.如权利要求2所述方法，其特征在于，步骤一所述直流电容的小信号模型建立过程如下：
建立直流电容的数学模型为
sCdcudc＝idc2-idc1＝1.5(ddridr+dqriqr)-1.5(ddgidg+dqgiqg)
式中，Cdc为直流电容，idc1为网侧直流电流，idc2为机侧直流电流，idg、iqg分别为网侧变流器交流端口d、q轴电流，ddg、dqg为dq坐标系下网侧变流器控制器输出占空比。udc为直流电压，idr、iqr分别为d、q坐标系发电机定子电流，ddr、dqr为dq坐标系下机侧变流器控制器输出占空比，对其进行线性化，可得直流电容小信号模型为
Figure FDA0003039022320000046
式中，
Figure FDA0003039022320000047
大写字母及上标*表示对应小写变量的稳态分量，Δ表示对应变量的小信号分量。
4.如权利要求3所述方法，其特征在于，步骤一所述网侧变流器及滤波器、网侧变流器控制器的小信号模型建立过程如下：
建立网侧变流器及滤波器的数学模型为
Figure FDA0003039022320000048
式中，Lf为滤波电感，ω为工频角频率，ω＝100πrad/s，idg、iqg分别为网侧变流器交流端口d、q轴电流，ddg、dqg为dq坐标系下网侧变流器控制器输出占空比，udc为直流电压，udg、uqg分别为并网点d、q轴电压。对该模型进行线性化，可得网侧变流器及滤波器的小信号模型为
Figure FDA0003039022320000051
其中，
Figure FDA0003039022320000052
大写字母及上标*表示对应小写变量的稳态分量，Δ表示对应变量的小信号分量。
建立网侧变流器控制器的数学模型为
Figure FDA0003039022320000053
其中，
Figure FDA0003039022320000054
Kcpg、Kcig分别为网侧电流环PI控制的比例参数和积分参数，
Figure FDA0003039022320000055
Kvp、Kvi分别为网侧电压环PI控制的比例参数和积分参数，Udcref为直流电压参考值。在网侧变流器控制器中，采用锁相环使风机与电网保持同步。上标c表示网侧变流器控制器dq坐标系。对该模型进行线性化，可得网侧变流器控制器的小信号模型为
Figure FDA0003039022320000056
其中，
Figure FDA0003039022320000057
此外，网侧变流器中还应考虑锁相环动态，其数学模型为
Figure FDA0003039022320000058
其中，
Figure FDA0003039022320000059
Kppll、Kipll分别为锁相环PI控制器的比例参数和积分参数，
Figure FDA00030390223200000510
为网侧变流器控制器dq坐标系下的并网点q轴电压。对其进行线性化，可得
Figure FDA00030390223200000511
其中，***dq坐标系与控制器dq坐标系存在一定偏差，二者之间可以通过如下方程相互转换
Figure FDA00030390223200000512
式中变量Δxd、Δxq可以表示网侧变流器输出电流Δidg、Δiqg、并网点电压Δudg、Δuqg或网侧控制器输出占空比Δddg、Δdqg，
Figure FDA00030390223200000513
表示对应变量的稳态分量
Figure FDA00030390223200000514
Figure FDA00030390223200000515
由此可推导出锁相环小信号模型，即
Δθ＝Gpll·Δuqg
其中，
Figure FDA0003039022320000061
从而可以得出控制器dq坐标系与***dq坐标系之间的关系为
Figure FDA0003039022320000062
其中，
Figure FDA0003039022320000063
则网侧变流器控制器的小信号模型转变为
Figure FDA0003039022320000064
5.如权利要求4所述方法，其特征在于，所述步骤三具体为：
交流弱电网采用理想电压源串联等效电感表示，建立其数学模型为
Figure FDA0003039022320000065
式中，Lg为弱电网等效电感，uds、uqs分别为d、q轴理想电压源电压，idg、iqg分别为网侧变流器交流端口d、q轴电流。将该式线性化，可得
Zg·Δidqg＝Δudqg
式中，
Figure FDA0003039022320000066
永磁同步风机的输出功率为
P＝1.5(idgudg+iqguqg)
将该式线性化可得
Figure FDA0003039022320000067
式中，
Figure FDA0003039022320000068
结合风力机、发电机及机侧变流器、机侧变流器控制器、直流电容、网侧变流器及滤波器、网侧变流器控制器小信号模型的线型部分，可得***线性部分传递函数G(s)为
Figure FDA0003039022320000069
Figure FDA00030390223200000610
Figure FDA0003039022320000071
Figure FDA0003039022320000072
Figure FDA0003039022320000073
Figure FDA0003039022320000074
Figure FDA0003039022320000075
Figure FDA0003039022320000076
Figure FDA0003039022320000077
式中，Tf表示功率采样滤波器周期，1/(1+Tfs)为功率采样滤波器延时，1/(1+1.5Tps)为控制器及PWM延时。
CN202110451859.4A 2021-04-26 2021-04-26 考虑功率控制的永磁同步风机接入弱电网稳定性分析方法 Active CN112994111B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110451859.4A CN112994111B (zh) 2021-04-26 2021-04-26 考虑功率控制的永磁同步风机接入弱电网稳定性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110451859.4A CN112994111B (zh) 2021-04-26 2021-04-26 考虑功率控制的永磁同步风机接入弱电网稳定性分析方法

Publications (2)

Publication Number Publication Date
CN112994111A true CN112994111A (zh) 2021-06-18
CN112994111B CN112994111B (zh) 2022-02-22

Family

ID=76340219

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110451859.4A Active CN112994111B (zh) 2021-04-26 2021-04-26 考虑功率控制的永磁同步风机接入弱电网稳定性分析方法

Country Status (1)

Country Link
CN (1) CN112994111B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115133578A (zh) * 2022-08-04 2022-09-30 四川大学 无网压传感器的电压源并网换流器功率自同步控制方法
CN117713088A (zh) * 2024-02-05 2024-03-15 云南电网有限责任公司 高铁电力机车并网对电网谐波影响的分析方法及***

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017016617A1 (de) * 2015-07-30 2017-02-02 Siemens Aktiengesellschaft Verfahren zur rechnergestützten parametrierung eines umrichters in einem stromnetz
CN109755964A (zh) * 2019-03-18 2019-05-14 哈尔滨工业大学 一种提高弱电网条件下双馈风电机组稳定性的控制方法
CN109861265A (zh) * 2018-12-28 2019-06-07 四川大学 一种风电场经mmc-hvdc接入弱电网的虚拟惯性控制方法
CN111342484A (zh) * 2019-11-26 2020-06-26 湖南工业大学 一种常规控制策略下直驱式风力发电***的动态特性分析方法
CN112260290A (zh) * 2020-10-13 2021-01-22 合肥工业大学 电压源型的永磁同步风电机组在弱网下的并网控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017016617A1 (de) * 2015-07-30 2017-02-02 Siemens Aktiengesellschaft Verfahren zur rechnergestützten parametrierung eines umrichters in einem stromnetz
CN109861265A (zh) * 2018-12-28 2019-06-07 四川大学 一种风电场经mmc-hvdc接入弱电网的虚拟惯性控制方法
CN109755964A (zh) * 2019-03-18 2019-05-14 哈尔滨工业大学 一种提高弱电网条件下双馈风电机组稳定性的控制方法
CN111342484A (zh) * 2019-11-26 2020-06-26 湖南工业大学 一种常规控制策略下直驱式风力发电***的动态特性分析方法
CN112260290A (zh) * 2020-10-13 2021-01-22 合肥工业大学 电压源型的永磁同步风电机组在弱网下的并网控制方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115133578A (zh) * 2022-08-04 2022-09-30 四川大学 无网压传感器的电压源并网换流器功率自同步控制方法
CN115133578B (zh) * 2022-08-04 2024-05-28 四川大学 无网压传感器的电压源并网换流器功率自同步控制方法
CN117713088A (zh) * 2024-02-05 2024-03-15 云南电网有限责任公司 高铁电力机车并网对电网谐波影响的分析方法及***
CN117713088B (zh) * 2024-02-05 2024-05-10 云南电网有限责任公司 高铁电力机车并网对电网谐波影响的分析方法及***

Also Published As

Publication number Publication date
CN112994111B (zh) 2022-02-22

Similar Documents

Publication Publication Date Title
CN111277001B (zh) 基于虚拟同步发电机参数自适应控制的风机并网控制方法
WO2022226709A1 (zh) 考虑功率控制的永磁同步风机接入弱电网稳定性分析方法
Zhang et al. Model-free predictive current control of a DFIG using an ultra-local model for grid synchronization and power regulation
CN107968591B (zh) 基于预测公共点电压的并网逆变器虚拟惯性功率解耦控制方法
Yang et al. A speed-adaptive reduced-order observer for sensorless vector control of doubly fed induction generator-based variable-speed wind turbines
BOUDJEMA et al. Second order sliding mode control of a dual-rotor wind turbine system by employing a matrix converter
CN108988391B (zh) 基于转速控制的双馈风机转子侧变换器的稳定性分析方法
CN111525611A (zh) 计及频率耦合效应的双馈并网***次同步振荡分析方法
CN112994111B (zh) 考虑功率控制的永磁同步风机接入弱电网稳定性分析方法
CN110676874B (zh) 计及频率耦合效应的直驱式风机次同步振荡电气量分析方法
Djoudi et al. Sliding mode control of DFIG powers in the case of unknown flux and rotor currents with reduced switching frequency
Pang et al. Stator harmonic current suppression for DFIG system considering integer harmonics and interharmonics
CN109039180B (zh) 双馈感应发电机并网过程的分数阶控制方法
CN112952901B (zh) 一种针对多风机并网***的分散式稳定性分析方法
Phan et al. An improved control strategy using a PI-resonant controller for an unbalanced stand-alone doubly-fed induction generator
CN110460113B (zh) 一种逆变器电源的机械惯性模拟方法
CN114928076B (zh) 无交流电压传感器的虚拟同步机双闭环控制方法
Kenne et al. Adaptive PI control strategy for a self-excited induction generator driven by a variable speed wind turbine
CN113783183B (zh) 弱电网下双馈风机在故障穿越期间的暂态稳定性评估方法
Phan et al. Enhanced proportional-resonant current controller for unbalanced stand-alone DFIG-based wind turbines
Mossa Field orientation control of a wind driven dfig connected to the grid
CN113328457A (zh) 一种提高并网友好性的风电机组虚拟控制方法
Hong et al. Grid voltage modulated direct power control for output smoothing in battery energy storage systems
Nasim et al. Field control grid connected dfig turbine System
Unais et al. Optimum Torque-Zero d-axis Current Control of Direct Driven PMSG Based Wind Energy Conversion System

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