CN110299712B - 一种直驱风机次同步振荡的控制方法及系统 - Google Patents

一种直驱风机次同步振荡的控制方法及系统 Download PDF

Info

Publication number
CN110299712B
CN110299712B CN201910695773.9A CN201910695773A CN110299712B CN 110299712 B CN110299712 B CN 110299712B CN 201910695773 A CN201910695773 A CN 201910695773A CN 110299712 B CN110299712 B CN 110299712B
Authority
CN
China
Prior art keywords
phase
svg
locked loop
loop
grid
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
CN201910695773.9A
Other languages
English (en)
Other versions
CN110299712A (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.)
North China Electric Power University
Original Assignee
North China Electric Power University
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN201910695773.9A priority Critical patent/CN110299712B/zh
Publication of CN110299712A publication Critical patent/CN110299712A/zh
Application granted granted Critical
Publication of CN110299712B publication Critical patent/CN110299712B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明涉及一种直驱风机次同步振荡的控制方法及系统,属于风力发电技术领域。该方法包括以下步骤:S1,当直驱风机发生次同步振荡时,获取直驱风机并联SVG并网系统的运行数据;S2,根据所述运行数据计算得到网侧变流器锁相环和电流环的耗散能量、SVG锁相环和电流环的耗散能量、直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量;S3,根据所述计算结果及控制参数调整最优化模型对网侧变流器锁相环和电流环的控制参数、SVG锁相环和电流环的控制参数进行调整,利用调整后的控制参数以对所述次同步振荡进行抑制。本发明实现直驱风机系统中次同步振荡主动控制。

Description

一种直驱风机次同步振荡的控制方法及系统
技术领域
本发明涉及风力发电技术领域,尤其涉及一种直驱风机次同步振荡的控制方法及系统。
背景技术
大规模开发利用新能源已成为世界各国发展的共识,特别是风电大量代替常规电源,这对电力系统的运行特性发生了深刻变化。一方面,系统惯量和阻尼水平持续弱化;另一方面,风电机组与电网相互耦合,诱发新的宽频带稳定问题,极大威胁电力系统的安全稳定运行。近年来,电网多次发生次同步振荡事故。因此,亟需对次同步振荡抑制方法进行深入研究。
目前,针对风电并网次同步振荡问题,国内外学者提出了多种次同步振荡抑制方法,但现有控制策略未考虑实际系统运行方式变化导致的频率漂移等问题,控制策略缺乏自适应性,且未充分考虑风机自身控制环节以及风电场无功补偿装置控制环节参与平抑系统振荡的能力,难以实现振荡的主动控制。
发明内容
鉴于上述的分析,本发明旨在提供一种直驱风机次同步振荡的控制系统及方法,以解决目前风电并网次同步振荡的控制策略存在频率漂移以及难以实现振荡的主动控制问题。
本发明的目的主要是通过以下技术方案实现的:
一方面,本发明提供了一种直驱风机次同步振荡的控制方法,包括以下步骤:S1,当直驱风机发生次同步振荡时,获取直驱风机并联SVG并网系统的运行数据;S2,根据所述运行数据计算得到网侧变流器锁相环和电流环的耗散能量、SVG锁相环和电流环的耗散能量、直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量;S3,根据所述计算结果及控制参数调整最优化模型对网侧变流器锁相环和电流环的控制参数、SVG锁相环和电流环的控制参数进行调整,利用调整后的控制参数以对所述次同步振荡进行抑制。
进一步地,所述运行数据包括:网侧变流器电压、电流的d、q轴分量u2d、u2q、i2d、i2q,直流电压udc,电网电角速度ω2,网侧变流器向电网输出有功P2,q轴外环控制目标的实测值交流电压us或无功q,次同步频率分量的电压幅值U-,基频电压幅值U0和次同步分量频率ωs
进一步地,所述网侧变流器锁相环和电流环的控制参数包括网侧变流器锁相环和电流环PI控制器的比例系数k、kp1和积分系数k、ki1;所述SVG锁相环和电流环的控制参数包括SVG锁相环和电流环PI控制器的比例系数k′、k′p1和积分系数k′、k′i1
进一步地,所述控制参数调整最优化模型为:
Figure BDA0002149327650000031
其中,WhPMSG为直驱风机控制环节的总耗散能量,WhSVG为SVG控制环节的总耗散能量,λ,λ′分别为直驱风机和SVG的目标函数加权系数,W′pll,W′I分别为SVG锁相环和电流环提供的耗散能量,Wpll,WI分别为网侧变流器锁相环和电流环提供的耗散能量,gi(k)为控制参数的约束条件,K为控制参数的可调整区间范围,Ω为次同步振荡频率在dq坐标系下对应的频率范围。
进一步地,所述控制参数调整最优化模型包括参数调整的目标函数和控制参数的约束条件;
所述参数调整的目标函数f(k,kpi)为锁相环动态角Δθpll幅值的计算式:
Figure BDA0002149327650000032
其中,A表示Δθpll的幅值,U-为次同步频率分量的电压幅值,U0为基频电压幅值,ωs为次同步分量频率;
Figure BDA0002149327650000033
k为网侧变流器锁相环PI控制器的比例系数,k为网侧变流器锁相环PI控制器的积分系数,u2q为网侧变流器电压q轴分量;
所述控制参数的约束条件包括直驱风机控制参数约束条件和SVG控制参数约束条件;
所述直驱风机控制参数约束条件为:
Figure BDA0002149327650000041
所述SVG控制参数约束条件为:
D′(k,k)<D3
进一步地,所述直驱风机控制参数约束条件中D1和D2的取值通过电流环环节控制参数调整策略得到,所述SVG控制参数约束条件中D3的取值通过SVG控制参数调整策略得到;
所述电流环环节控制参数调整策略包括:
令Wkp=0,得到改变Wkp能量特性的锁相环动态角临界幅值A1,对应得到锁相环比例系数k和积分系数k取值对应的D1
当D∈(0,D1)时,Wkp<0,增大电流环比例系数kp,直驱风机输出的耗散能量减小;
Figure BDA0002149327650000043
时,Wkp>0,减小电流环比例系数kp,直驱风机输出的耗散能量减小;
令Wki=0,得到锁相环动态角幅值A2与对应的D2
当D∈(0,D2)时,Wki<0,增大电流环积分系数ki,直驱风机输出的耗散能量减小;
Figure BDA0002149327650000042
时,Wki>0,减小电流环积分系数ki,直驱风机输出的耗散能量减小;
所述SVG控制参数调整策略包括:
保持锁相环参数不变,调整电流环比例系数令Wpll=0,求得此时电流环比例系数的临界值kp1=C1
当kp1>C1时,Wpll<0,锁相环提供负的耗散能量,此时增大锁相环的比例、积分系数使SVG发出的耗散能量减小;
当kp1<C1时,Wpll>0,锁相环提供正的耗散能量,此时减小锁相环的比例、积分系数使SVG发出的耗散能量减小;
令Wkp=0,改变锁相环动态角的幅值A,求得此时的临界值A=A3,对应得到(k,k)取值的边界D3
当D∈(0,D3)时,Wkp>0,减小电流环比例系数kp1,SVG输出的耗散能量减小;
Figure BDA0002149327650000051
时,Wki<0,增大电流环比例系数kp1,SVG输出的耗散能量减小。
进一步地,所述电流环环节控制参数调整策略和所述SVG控制参数调整策略的能量特征Wkp和Wki表达式为:
Figure BDA0002149327650000052
进一步地,所述网侧变流器锁相环耗散能量的计算模型为:
Figure BDA0002149327650000053
所述网侧变流器电流环耗散能量的计算模型为:
Figure BDA0002149327650000054
所述SVG锁相环耗散能量的计算模型为:
Figure BDA0002149327650000055
所述SVG电流环耗散能量的计算模型为:
Figure BDA0002149327650000061
所述直驱风机控制环节的总耗散能量的计算模型为:
Figure BDA0002149327650000062
所述SVG控制环节的总耗散能量的计算模型为:
Figure BDA0002149327650000063
其中,i2d、i2q分别为网侧变流器电流的d、q轴分量,Δθpll为锁相环动态角,L2为进线电抗器的等效电抗值,udc为直流电压,
Figure BDA0002149327650000064
为直流电压参考值,ω2为电网电角速度,
Figure BDA0002149327650000065
为d轴电流参考值,
Figure BDA0002149327650000066
为q轴电流参考值,P2为网侧变流器向电网输出有功,x*为q轴外环控制目标的参考值,,x为q轴外环控制目标的实测值。
进一步地,所述利用调整后的控制参数以对所述次同步振荡进行抑制包括改变控制参数后按照以下公式计算锁相环的耗散能量:
Figure BDA0002149327650000067
其中:D=A/U-
Figure BDA0002149327650000068
D改变控制参数后表示为
Figure BDA0002149327650000069
Figure BDA00021493276500000610
U-为次同步频率分量的电压幅值,ωs为次同步分量频率,k,k为锁相环PI控制器的比例系数、积分系数,U0为基频电压幅值,Wpll为网侧变流器锁相环环节的耗散能量。
另一方面,一种直驱风机次同步振荡的控制系统,包括数据采集模块、耗散能量分析模块、参数调整模块、结果输出模块;
所述数据采集模块,用于采集直驱风机并联SVG并网系统的运行数据;
所述耗散能量分析模块,用于根据耗散能量分析模型及所述数据采集模块输出的数据分别计算并分析网侧变流器锁相环和电流环的耗散能量、SVG锁相环和电流环的耗散能量、直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量;
所述参数调整模块,用于根据控制参数调整最优化模型及所述耗散能量分析模块输出的各环节耗散能量调整网侧变流器锁相环和电流环以及SVG锁相环和电流环的控制参数;
所述结果输出模块,用于输出调整后的网侧变流器锁相环和电流环的控制参数、SVG锁相环和电流环的控制参数,还用于输出控制参数调整后直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量。
本发明技术方案的有益效果:本发明公开了一种直驱风机次同步振荡的控制方法及系统,具有较强自适应性,且充分考虑直驱风机自身控制环节以及风电场无功补偿装置控制环节参与平抑系统振荡的能力,实现直驱风机系统中次同步振荡主动控制。
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分的从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制,在整个附图中,相同的参考符号表示相同的部件。
图1为本发明实施例的一种直驱风机次同步振荡的控制方法流程图;
图2为本发明实施例的直驱风机并联SVG的并网系统;
图3为本发明实施例的PMSG耗散能量对比图;
图4为本发明实施例的SVG耗散能量对比图;
图5为本发明实施例的总耗散能量对比图;
图6为本发明实施例的有功功率输出对比图;
图7为本发明实施例的一种直驱风机次同步振荡的控制系统结构图。
具体实施方式
下面结合附图来具体描述本发明的优选实施例,其中,附图构成本申请一部分,并与本发明的实施例一起用于阐释本发明的原理,并非用于限定本发明的范围。
本发明的一个具体实施例,如图1所示,公开了一种直驱风机次同步振荡的控制方法,包括以下步骤:
S1,当直驱风机发生次同步振荡时,获取直驱风机并联SVG并网系统的运行数据;
S2,根据所述运行数据计算得到网侧变流器锁相环和电流环的耗散能量、SVG锁相环和电流环的耗散能量、直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量;
S3,根据所述计算结果及控制参数调整最优化模型对网侧变流器锁相环和电流环的控制参数、SVG锁相环和电流环的控制参数进行调整,利用调整后的控制参数以对所述次同步振荡进行抑制。
与现有技术相比,本发明实施例基于控制参数调整的次同步振荡最优控制模型具有较强的自适应性,且充分利用直驱风机自身控制环节以及风电场无功补偿装置控制环节平抑系统振荡的能力,实现直驱风机系统中次同步振荡主动控制。
本发明的一个具体实施例,所述运行数据包括:网侧变流器电压、电流的d、q轴分量u2d、u2q、i2d、i2q,直流电压udc,电网电角速度ω2,网侧变流器向电网输出有功P2,q轴外环控制目标的实测值交流电压us或无功q,次同步频率分量的电压幅值U-,基频电压幅值U0和次同步分量频率ωs。也即,本发明实施例充分利用直驱风机自身控制环节以及风电场无功补偿装置控制环节的运行数据,为实现直驱风机系统中次同步振荡主动控制提供了充足的数据。
本发明的一个具体实施例,所述网侧变流器锁相环和电流环的控制参数包括网侧变流器锁相环和电流环PI控制器的比例系数k、kp1和积分系数k、ki1;所述SVG锁相环和电流环的控制参数包括SVG锁相环和电流环PI控制器的比例系数k′、k′p1和积分系数k′、k′i1
本发明的一个具体实施例,所述控制参数调整最优化模型为:
Figure BDA0002149327650000101
其中,WhPMSG为直驱风机控制环节的总耗散能量,WhSVG为SVG控制环节的总耗散能量,λ,λ′分别为直驱风机和SVG的目标函数加权系数,W′pll,W′I分别为SVG锁相环和电流环提供的耗散能量,Wpll,WI分别为网侧变流器锁相环和电流环提供的耗散能量,gi(k)为控制参数的约束条件,K为控制参数的可调整区间范围,Ω为次同步振荡频率在dq坐标系下对应的频率范围。根据暂态过程中的电压支撑,以及直驱风机与SVG并联母线的无功水平等约束,得到各控制参数的可调整区间K。
本发明的一个具体实施例,所述控制参数调整最优化模型包括参数调整的目标函数和控制参数的约束条件;
所述参数调整的目标函数f(k,kpi)为锁相环动态角Δθpll幅值的计算式:
Figure BDA0002149327650000102
其中,A表示Δθpll的幅值,U-为次同步频率分量的电压幅值,U0为基频电压幅值,ωs为次同步分量频率;
Figure BDA0002149327650000103
k为网侧变流器锁相环PI控制器的比例系数,k为网侧变流器锁相环PI控制器的积分系数,u2q为网侧变流器电压q轴分量;
需要说明的是,对于锁相环动态角Δθpll的积分表达式中,积分对象为控制系统dq坐标下的q轴电流,含有Δθpll的二阶扰动量,忽略Δθpll的二阶扰动项时:
Figure BDA0002149327650000111
Δθpll与q轴电流的频率一致,可以表示为:
Figure BDA0002149327650000118
式中:A和
Figure BDA0002149327650000112
分别表示Δθpll的幅值与相位,A和
Figure BDA0002149327650000113
的表达式分别为:
Figure BDA0002149327650000114
可知,Δθpll的幅值主要受次同步频率分量电压幅值U-,次同步分量频率ωs,基频电压幅值U0以及锁相环比例积分系数影响。
所述控制参数的约束条件包括直驱风机控制参数约束条件和SVG控制参数约束条件;
所述直驱风机控制参数约束条件为:
Figure BDA0002149327650000115
所述SVG控制参数约束条件为:
D′(k,k)<D3
需要说明的是,SVG控制参数理想的约束条件为:
Figure BDA0002149327650000116
在实际计算中,式(-∫ΔθplldP2)和
Figure BDA0002149327650000117
为正,且数量级与受kp1影响的两项相当,所以Wpll<0难以满足。此时,SVG的控制参数约束条件可修改为:
D′(k,k)<D3
本发明的一个具体实施例,所述直驱风机控制参数约束条件中D1和D2的取值通过电流环环节控制参数调整策略得到,所述SVG控制参数约束条件中D3的取值通过SVG控制参数调整策略得到;
所述电流环环节控制参数调整策略包括:
令Wkp=0,得到改变Wkp能量特性的锁相环动态角临界幅值A1,对应得到锁相环比例系数k和积分系数k取值对应的D1
当D∈(0,D1)时,Wkp<0,增大电流环比例系数kp,直驱风机输出的耗散能量减小;
Figure BDA0002149327650000121
时,Wkp>0,减小电流环比例系数kp,直驱风机输出的耗散能量减小;
令Wki=0,得到锁相环动态角幅值A2与对应的D2
当D∈(0,D2)时,Wki<0,增大电流环积分系数ki,直驱风机输出的耗散能量减小;
Figure BDA0002149327650000122
时,Wki>0,减小电流环积分系数ki,直驱风机输出的耗散能量减小;
所述SVG控制参数调整策略包括:
保持锁相环参数不变,调整电流环比例系数令Wpll=0,求得此时电流环比例系数的临界值kp1=C1
当kp1>C1时,Wpll<0,锁相环提供负的耗散能量,此时增大锁相环的比例、积分系数使SVG发出的耗散能量减小;
当kp1<C1时,Wpll>0,锁相环提供正的耗散能量,此时减小锁相环的比例、积分系数使SVG发出的耗散能量减小;
令Wkp=0,改变锁相环动态角的幅值A,求得此时的临界值A=A3,对应得到(k,k)取值的边界D3
当D∈(0,D3)时,Wkp>0,减小电流环比例系数kp1,SVG输出的耗散能量减小;
Figure BDA0002149327650000131
时,Wki<0,增大电流环比例系数kp1,SVG输出的耗散能量减小。
本发明的一个具体实施例,所述电流环环节控制参数调整策略和所述SVG控制参数调整策略的能量特征Wkp和Wki表达式为:
Figure BDA0002149327650000132
需要说明的是,由于控制环节的相互耦合,电流环的控制参数不仅仅对电流环的耗散能量有所影响,同时还影响锁相环的耗散能量。电流环的控制参数不决定锁相环耗散能量的正负,但会影响其绝对值大小。单纯的增大电流环比例、积分系数,虽然电流环消耗的耗散能量增加,但同时,锁相环发出的耗散能量也有可能增加,因此,在调整电流环参数时,要综合考虑这两部能量的大小,因此将比例、积分系数的影响分开讨论:
Figure BDA0002149327650000133
Figure BDA0002149327650000134
在网侧变流器控制的时间尺度下,认为风速不发生变化,即机侧获取的风功率不变,那么此时d轴电流参考值
Figure BDA0002149327650000135
主要受直流电容电压波动的影响,变化很小,且
Figure BDA0002149327650000136
可近似认为不变,所以Wkp和Wki的正负取决于锁相环参数。在Wkp表达式中后两项能量为正,且均与Δθpll有关,dΔθpll的幅值与A成正比。
本发明的一个具体实施例,所述网侧变流器锁相环耗散能量的计算模型为:
Figure BDA0002149327650000141
所述网侧变流器电流环耗散能量的计算模型为:
Figure BDA0002149327650000142
需要说明的是,在直驱风机中,所述电流环耗散能量模型表达式中的两部分能量均为负值,电流环提供负的耗散能量,意味着暂态过程中吸收能量,对振荡起到正阻尼作用,有利于平息次同步振荡。因此,应主动增加这一部分能量的绝对值,即增大电流环的比例、积分系数,以增加负的耗散能量。
所述SVG锁相环耗散能量的计算模型为:
Figure BDA0002149327650000143
所述SVG电流环耗散能量的计算模型为:
Figure BDA0002149327650000144
所述直驱风机控制环节的总耗散能量的计算模型为:
Figure BDA0002149327650000145
需要说明的是,若直驱风机控制环节提供的总耗散能量为正,则风机在暂态过程中向系统输入的耗散能量不断增加,对振荡的影响为加剧振荡发散;若直驱风机控制环节提供的耗散能量为负,则风机在暂态过程中向系统输入的耗散能量逐渐收敛,表现出的特性为抑制振荡发展。
所述SVG控制环节的总耗散能量的计算模型为:
Figure BDA0002149327650000151
其中,i2d、i2q分别为网侧变流器电流的d、q轴分量,Δθpll为锁相环动态角,L2为进线电抗器的等效电抗值,udc为直流电压,
Figure BDA0002149327650000152
为直流电压参考值,ω2为电网电角速度,
Figure BDA0002149327650000153
为d轴电流参考值,
Figure BDA0002149327650000154
为q轴电流参考值,P2为网侧变流器向电网输出有功,x*为q轴外环控制目标的参考值,即恒交流电压us或恒无功q,x为q轴外环控制目标的实测值,即交流电压us或无功q。
需要说明的是,次同步振荡过程中,SVG工作在恒电压控制方式下耗散能量受振荡影响更大,且SVG在两种控制方式下均有可能导致振荡发散。以SVG工作在恒电压控制方式下进行分析,得到的规律同样适用于恒无功控制方式下。
由于SVG输出的q轴电流与d轴电流在同一数量级上,因此表达式∫(i2ddi2q-i2qdi2d)的正负不能直观的判断出来,当其大于零时,SVG的参数调整策略与直驱风机相同。
以下分析∫(i2ddi2q-i2qdi2d)<0时,控制参数变化对耗散能量的影响:所述锁相环耗散能量表达式中ω2L2kp1∫[Δθpll(i2ddi2q-i2qdi2d)]将变为负值,其他三项符号不变,而这一项同时受锁相环参数和电流环比例系数影响,决定了电流环和锁相环的参数调整策略,因此将两个环节的控制参数调整统一考虑。
本发明的一个具体实施例,所述利用调整后的控制参数以对所述次同步振荡进行抑制包括改变控制参数后按照以下公式计算锁相环的耗散能量:
Figure BDA0002149327650000161
其中:D=A/U-
Figure BDA0002149327650000162
D改变控制参数后表示为
Figure BDA0002149327650000163
Figure BDA0002149327650000164
U-为次同步频率分量的电压幅值,ωs为次同步分量频率,k,k为锁相环PI控制器的比例系数、积分系数,U0为基频电压幅值,Wpll为网侧变流器锁相环环节的耗散能量。
需要说明的是,由于次同步频率扰动出现时,其初相位已经确定,在不改变锁相环延时的前提下,可认为Δθpll的相位不发生变化。因此改变锁相环比例、积分系数不改变乘积项:Δθpll(i2ddi2q-i2qdi2d)、ΔθplldP2
Figure BDA0002149327650000165
的正负,而改变控制参数引起的Δθpll幅值变化ΔA,可等价于所述锁相环环节的耗散能量表达式整体乘以系数
Figure BDA0002149327650000166
记:D=A/U-,则锁相环参数改变后,新的耗散能量大小为:
Figure BDA0002149327650000167
Figure BDA0002149327650000168
时,锁相环提供的耗散能量减小,因而控制环节提供的总耗散能量减小,有利于振荡收敛,得到使耗散能量减小的锁相环的比例、积分系数变化规律为:
ΔD<0;
锁相环比例、积分系数一个保持不变另一个参数减小或两个参数都减小时,ΔD<0;而两个变化趋势不相同时,ΔD的正负取决于参数调整的相对大小。
实际应用中,如图2所示直驱风机并联SVG并网系统,利用所述控制参数调整最优化模型计算出调整后的系统控制参数,比较在两组控制参数下的次同步振荡情况,验证结果的正确性。
未进行参数调整时直驱风机锁相环和电流环控制参数为:kp1=0.034,ki1=1.68,k=3.26,k=65.2;SVG的控制参数为:k′p1=0.023,k′i1=1.15,k=3.26,k=65.2。
调整后直驱风机锁相环和电流控制参数为:kp1=0.121,ki1=6.166,k=0.55,k=25.39;SVG的控制参数为:k′p1=0.53,k′i1=10.84,k′=0.80,k′=38.3。
系统中发生次同步振荡时,分别计算两种参数控制下直驱风机与并联SVG的耗散能量,如图3与4所示,实施参数调整后,直驱风机和SVG均发出负的耗散能量,验证了参数调整策略有效性。
直驱风机和SVG的总耗散能量如图5所示,施加参数调整之后,直驱风机与SVG的总耗散能量由正变负。如图6有功功率图所示,施加参数调整后系统次同步振荡在0.3s时间内收敛,实现了振荡主动抑制,且要快于不并联SVG系统的收敛速度,因此SVG控制参数的调整增强了振荡抑制的效果。
由此仿真结果验证了本发明技术方案提出的基于直驱风机和SVG控制参数优化调整的次振荡主动抑制方法。
本发明的一个具体实施例,如图7所示,一种直驱风机次同步振荡的控制系统,包括数据采集模块、耗散能量分析模块、参数调整模块、结果输出模块;
所述数据采集模块,用于采集直驱风机并联SVG并网系统的运行数据;
所述耗散能量分析模块,用于根据耗散能量分析模型及所述数据采集模块输出的数据分别计算并分析网侧变流器锁相环和电流环的耗散能量、SVG锁相环和电流环的耗散能量、直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量;
所述参数调整模块,用于根据控制参数调整最优化模型及所述耗散能量分析模块输出的各环节耗散能量调整网侧变流器锁相环和电流环以及SVG锁相环和电流环的控制参数;
所述结果输出模块,用于输出调整后的网侧变流器锁相环和电流环的控制参数、SVG锁相环和电流环的控制参数,还用于输出控制参数调整后直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量。
需要说明的是,所述数据采集模块采集系统的运行数据,发送给具有耗散能量分析模型的耗散能量分析模块,所述耗散能量分析模块输出网侧变流器锁相环和电流环、SVG锁相环和电流环、直驱风机控制环节和SVG控制环节各环节耗散能量,发送给具有控制参数调整最优化模型的参数调整模块,所述参数调整模块输出的各环节耗散能量调整网侧变流器锁相环和电流环以及SVG锁相环和电流环的控制参数,改变后的控制参数由所述结果输出模块输出,并反馈给所述耗散能量分析模块后由所述结果输出模块输出控制参数调整后直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量。以上所述各模块之间的连接关系如图7所示。
综上所述,本发明公开了一种直驱风机次同步振荡的控制方法,包括以下步骤:S1,当直驱风机发生次同步振荡时,获取直驱风机并联SVG并网系统的运行数据;S2,根据所述运行数据计算得到网侧变流器锁相环和电流环的耗散能量、SVG锁相环和电流环的耗散能量、直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量;S3,根据所述计算结果及控制参数调整最优化模型对网侧变流器锁相环和电流环的控制参数、SVG锁相环和电流环的控制参数进行调整,利用调整后的控制参数以对所述次同步振荡进行抑制。同时公开了与上述方法构成同一发明构思的控制系统,包括数据采集模块、耗散能量分析模块、参数调整模块、结果输出模块;所述数据采集模块,用于采集直驱风机并联SVG并网系统的运行数据;所述耗散能量分析模块,用于根据耗散能量分析模型及所述数据采集模块输出的数据分别计算并分析网侧变流器锁相环和电流环的耗散能量、SVG锁相环和电流环的耗散能量、直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量;所述参数调整模块,用于根据控制参数调整最优化模型及所述耗散能量分析模块输出的各环节耗散能量调整网侧变流器锁相环和电流环以及SVG锁相环和电流环的控制参数;所述结果输出模块,用于输出调整后的网侧变流器锁相环和电流环的控制参数、SVG锁相环和电流环的控制参数,还用于输出控制参数调整后直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量。本发明技术方案建立适用于网侧变流器锁相环控制、电流环控制以及并联SVG控制环节的耗散能量模型,计算网侧变流器锁相环、电流环以及并联SVG控制环节的耗散能量,并基于控制参数调整的次同步振荡最优化模型。与现有技术相比,该系统及方法具有较强自适应性,且充分考虑直驱风机自身控制环节以及风电场无功补偿装置控制环节参与平抑系统振荡的能力,实现直驱风机系统中次同步振荡主动控制。
本领域技术人员可以理解,实现上述实施例中方法的全部或部分流程,可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于计算机可读存储介质中。其中,所述计算机可读存储介质为磁盘、光盘、只读存储记忆体或随机存储记忆体等。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。

Claims (9)

1.一种直驱风机次同步振荡的控制方法,其特征在于,包括以下步骤:
S1,当直驱风机发生次同步振荡时,获取直驱风机并联SVG并网系统的运行数据;
S2,根据所述运行数据计算得到网侧变流器锁相环和电流环的耗散能量、SVG锁相环和电流环的耗散能量、直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量;
S3,根据计算结果及控制参数调整最优化模型对网侧变流器锁相环和电流环的控制参数、SVG锁相环和电流环的控制参数进行调整,利用调整后的控制参数以对所述次同步振荡进行抑制;
所述控制参数调整最优化模型为:
Figure FDA0002813923100000011
其中,WhPMSG为直驱风机控制环节的总耗散能量,WhSVG为SVG控制环节的总耗散能量,λ,λ′分别为直驱风机和SVG的目标函数加权系数,W′pll,W′I分别为SVG锁相环和电流环提供的耗散能量,Wpll,WI分别为网侧变流器锁相环和电流环提供的耗散能量,gi(k)为控制参数的约束条件,K为控制参数的可调整区间范围,Ω为次同步振荡频率在dq坐标系下对应的频率范围,ωs为次同步分量频率,k,k分别为网侧变流器锁相环PI控制器的比例系数、积分系数,kp1,ki1分别为网侧变流器电流环PI控制器的比例系数、积分系数,k′,k′分别为SVG锁相环PI控制器的比例系数、积分系数,k′p1,k′i1分别为SVG电流环PI控制器的比例系数、积分系数,D=A/U-,A为锁相环动态角Δθpll的幅值,U-为次同步频率分量的电压幅值,D1、D2和D3分别为与锁相环动态角的临界幅值A1、A2和A3对应的值。
2.根据权利要求1所述的方法,其特征在于,所述运行数据包括:网侧变流器电压、电流的d、q轴分量u2d、u2q、i2d、i2q,直流电压udc,电网电角速度ω2,网侧变流器向电网输出有功P2,q轴外环控制目标的实测值交流电压us或无功q,次同步频率分量的电压幅值U-,基频电压幅值U0和次同步分量频率ωs
3.根据权利要求1所述的方法,其特征在于,所述网侧变流器锁相环和电流环的控制参数包括网侧变流器锁相环和电流环PI控制器的比例系数k、kp1和积分系数k、ki1;所述SVG锁相环和电流环的控制参数包括SVG锁相环和电流环PI控制器的比例系数k′、k′p1和积分系数k′、k′i1
4.根据权利要求1所述的方法,其特征在于,所述控制参数调整最优化模型包括参数调整的目标函数和控制参数的约束条件;
所述参数调整的目标函数f(k,kpi)为锁相环动态角Δθpll幅值的计算式:
Figure FDA0002813923100000021
其中,A表示Δθpll的幅值,U-为次同步频率分量的电压幅值,U0为基频电压幅值,ωs为次同步分量频率;
Figure FDA0002813923100000022
k为网侧变流器锁相环PI控制器的比例系数,k为网侧变流器锁相环PI控制器的积分系数,u2q为网侧变流器电压q轴分量;
所述控制参数的约束条件包括直驱风机控制参数约束条件和SVG控制参数约束条件;
所述直驱风机控制参数约束条件为:
Figure FDA0002813923100000031
所述SVG控制参数约束条件为:
D′(k,k)<D3
5.根据权利要求4所述的方法,其特征在于,所述直驱风机控制参数约束条件中D1和D2的取值通过电流环环节控制参数调整策略得到,所述SVG控制参数约束条件中D3的取值通过SVG控制参数调整策略得到;
所述电流环环节控制参数调整策略包括:
令Wkp=0,得到改变Wkp能量特性的锁相环动态角临界幅值A1,对应得到锁相环比例系数k和积分系数k取值对应的D1
当D∈(0,D1)时,Wkp<0,增大电流环比例系数kp,直驱风机输出的耗散能量减小;
Figure FDA0002813923100000032
时,Wkp>0,减小电流环比例系数kp,直驱风机输出的耗散能量减小;
令Wki=0,得到锁相环动态角幅值A2与对应的D2
当D∈(0,D2)时,Wki<0,增大电流环积分系数ki,直驱风机输出的耗散能量减小;
Figure FDA0002813923100000033
时,Wki>0,减小电流环积分系数ki,直驱风机输出的耗散能量减小;
所述SVG控制参数调整策略包括:
保持锁相环参数不变,调整电流环比例系数令Wpll=0,求得此时电流环比例系数的临界值kp1=C1
当kp1>C1时,Wpll<0,锁相环提供负的耗散能量,此时增大锁相环的比例、积分系数使SVG发出的耗散能量减小;
当kp1<C1时,Wpll>0,锁相环提供正的耗散能量,此时减小锁相环的比例、积分系数使SVG发出的耗散能量减小;
令Wkp=0,改变锁相环动态角的幅值A,求得此时的临界值A=A3,对应得到(k,k)取值的边界D3
当D∈(0,D3)时,Wkp>0,减小电流环比例系数kp1,SVG输出的耗散能量减小;
Figure FDA0002813923100000044
时,Wki<0,增大电流环比例系数kp1,SVG输出的耗散能量减小;
其中,Wkp和Wki分别为电流环环节控制参数调整策略和SVG控制参数调整策略的能量特征。
6.根据权利要求5所述的方法,其特征在于,所述电流环环节控制参数调整策略和所述SVG控制参数调整策略的能量特征Wkp和Wki表达式为:
Figure FDA0002813923100000041
其中,i2d、i2q分别为网侧变流器电流的d、q轴分量,L2为进线电抗器的等效电抗值,ω2为电网电角速度,
Figure FDA0002813923100000042
为d轴电流参考值,
Figure FDA0002813923100000043
为q轴电流参考值,kp1,ki1分别为网侧变流器电流环PI控制器的比例系数、积分系数。
7.根据权利要求1所述的方法,其特征在于,所述网侧变流器锁相环耗散能量的计算模型为:
Figure FDA0002813923100000051
所述网侧变流器电流环耗散能量的计算模型为:
Figure FDA0002813923100000052
所述SVG锁相环耗散能量的计算模型为:
Figure FDA0002813923100000053
所述SVG电流环耗散能量的计算模型为:
Figure FDA0002813923100000054
所述直驱风机控制环节的总耗散能量的计算模型为:
Figure FDA0002813923100000055
所述SVG控制环节的总耗散能量的计算模型为:
Figure FDA0002813923100000056
其中,i2d、i2q分别为网侧变流器电流的d、q轴分量,Δθpll为锁相环动态角,L2为进线电抗器的等效电抗值,udc为直流电压,
Figure FDA0002813923100000061
为直流电压参考值,ω2为电网电角速度,
Figure FDA0002813923100000062
为d轴电流参考值,
Figure FDA0002813923100000063
为q轴电流参考值,P2为网侧变流器向电网输出有功,x*为q轴外环控制目标的参考值,x为q轴外环控制目标的实测值,k,k分别为网侧变流器锁相环PI控制器的比例系数、积分系数,kp1,ki1分别为网侧变流器电流环PI控制器的比例系数、积分系数,k′,k′分别为SVG锁相环PI控制器的比例系数、积分系数,k′p1,k′i1分别为SVG电流环PI控制器的比例系数、积分系数,kp2,ki2分别为直流电压外环PI控制器的比例系数、积分系数,k′p2,k′i2分别为SVG电压外环PI控制器的比例系数、积分系数。
8.根据权利要求1或7所述的方法,其特征在于,所述利用调整后的控制参数以对所述次同步振荡进行抑制包括改变控制参数后按照以下公式计算锁相环的耗散能量:
Figure FDA0002813923100000064
其中:D=A/U-
Figure FDA0002813923100000065
D改变控制参数后表示为
Figure FDA0002813923100000066
Figure FDA0002813923100000067
U-为次同步频率分量的电压幅值,ωs为次同步分量频率,k,k为锁相环PI控制器的比例系数、积分系数,U0为基频电压幅值,Wpll为网侧变流器锁相环环节的耗散能量。
9.一种直驱风机次同步振荡的控制系统,其特征在于,包括数据采集模块、耗散能量分析模块、参数调整模块、结果输出模块;
所述数据采集模块,用于采集直驱风机并联SVG并网系统的运行数据;
所述耗散能量分析模块,用于根据耗散能量分析模型及所述数据采集模块输出的数据分别计算并分析网侧变流器锁相环和电流环的耗散能量、SVG锁相环和电流环的耗散能量、直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量;
所述参数调整模块,用于根据控制参数调整最优化模型及所述耗散能量分析模块输出的各环节耗散能量调整网侧变流器锁相环和电流环以及SVG锁相环和电流环的控制参数;
所述控制参数调整最优化模型为:
Figure FDA0002813923100000071
其中,WhPMSG为直驱风机控制环节的总耗散能量,WhSVG为SVG控制环节的总耗散能量,λ,λ′分别为直驱风机和SVG的目标函数加权系数,W′pll,W′I分别为SVG锁相环和电流环提供的耗散能量,Wpll,WI分别为网侧变流器锁相环和电流环提供的耗散能量,gi(k)为控制参数的约束条件,K为控制参数的可调整区间范围,Ω为次同步振荡频率在dq坐标系下对应的频率范围,ωs为次同步分量频率,k,k分别为网侧变流器锁相环PI控制器的比例系数、积分系数,kp1,ki1分别为网侧变流器电流环PI控制器的比例系数、积分系数,k′,k′分别为SVG锁相环PI控制器的比例系数、积分系数,k′p1,k′i1分别为SVG电流环PI控制器的比例系数、积分系数,D=A/U-,A为锁相环动态角Δθpll的幅值,U-为次同步频率分量的电压幅值,D1、D2和D3分别为与锁相环动态角的临界幅值A1、A2和A3对应的值;
所述结果输出模块,用于输出调整后的网侧变流器锁相环和电流环的控制参数、SVG锁相环和电流环的控制参数,还用于输出控制参数调整后直驱风机控制环节的总耗散能量以及SVG控制环节的总耗散能量。
CN201910695773.9A 2019-07-30 2019-07-30 一种直驱风机次同步振荡的控制方法及系统 Active CN110299712B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910695773.9A CN110299712B (zh) 2019-07-30 2019-07-30 一种直驱风机次同步振荡的控制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910695773.9A CN110299712B (zh) 2019-07-30 2019-07-30 一种直驱风机次同步振荡的控制方法及系统

Publications (2)

Publication Number Publication Date
CN110299712A CN110299712A (zh) 2019-10-01
CN110299712B true CN110299712B (zh) 2021-02-09

Family

ID=68032276

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910695773.9A Active CN110299712B (zh) 2019-07-30 2019-07-30 一种直驱风机次同步振荡的控制方法及系统

Country Status (1)

Country Link
CN (1) CN110299712B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113471981B (zh) * 2020-03-31 2022-08-09 新疆金风科技股份有限公司 风电场群集中并网的无功控制方法、装置及系统
CN111864767B (zh) * 2020-07-15 2022-04-15 华北电力大学 一种双馈风电机组次同步振荡主动阻尼控制方法及系统
CN112186758B (zh) * 2020-09-28 2022-12-09 西安热工研究院有限公司 一种自适应捕捉频点的直驱风电机组次同步振荡抑制方法
CN113346523B (zh) * 2021-07-05 2022-05-20 华北电力大学 一种风电场次同步振荡风险评估、抑制方法及系统
CN114759593B (zh) * 2021-08-26 2024-08-06 华北电力大学 基于能量耗散率的虚拟同步双馈风机参数优化方法及系统
CN114759575B (zh) * 2021-08-26 2024-08-06 华北电力大学 一种虚拟同步双馈风机次同步振荡抑制方法及系统
CN114744665B (zh) * 2021-09-03 2024-08-06 华北电力大学 一种风电系统经柔直外送的振荡传导路径确定方法及装置
CN114759577B (zh) * 2022-02-10 2024-08-06 华北电力大学 一种直驱风电柔直外送系统的次同步振荡抑制方法及系统
CN114744628B (zh) * 2022-02-28 2024-08-02 华北电力大学 基于耗散能量的交直流系统故障暂态特性优化方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108599236A (zh) * 2018-04-24 2018-09-28 华北电力科学研究院有限责任公司 双馈风电场次同步振荡svg抑制方法及装置
CN109217362A (zh) * 2018-08-22 2019-01-15 华北电力大学 一种双馈风机并网系统低频振荡扰动源定位系统及方法
CN109888776A (zh) * 2019-03-12 2019-06-14 深圳大学 针对直驱型风电场次同步谐振频率的预测方法及终端设备
CN109962495A (zh) * 2019-01-11 2019-07-02 四川大学 一种超低频振荡扰动源定位及抑制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108599236A (zh) * 2018-04-24 2018-09-28 华北电力科学研究院有限责任公司 双馈风电场次同步振荡svg抑制方法及装置
CN109217362A (zh) * 2018-08-22 2019-01-15 华北电力大学 一种双馈风机并网系统低频振荡扰动源定位系统及方法
CN109962495A (zh) * 2019-01-11 2019-07-02 四川大学 一种超低频振荡扰动源定位及抑制方法
CN109888776A (zh) * 2019-03-12 2019-06-14 深圳大学 针对直驱型风电场次同步谐振频率的预测方法及终端设备

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Stability analysis of power system with multiple operating conditions considering the stochastic characteristic of wind speed;Jing Ma ; Yang Qiu ; Yinan Li ; James S. Thorp;《IET Generation, Transmission & Distribution》;20161231;第10卷(第4期);1056-1066 *

Also Published As

Publication number Publication date
CN110299712A (zh) 2019-10-01

Similar Documents

Publication Publication Date Title
CN110299712B (zh) 一种直驱风机次同步振荡的控制方法及系统
WO2024021206A1 (zh) 一种基于构网型变流器的储能系统控制方法、系统、存储介质及设备
WO2020252813A1 (zh) 逆变型分布式电源双层自适应惯量控制方法及装置
US9124140B2 (en) Intelligent power converter control for grid integration of renewable energies
CN114759574B (zh) 一种并网直驱风电机组振荡主动阻尼控制方法及系统
CN108964040B (zh) 电网不平衡下虚拟同步发电机功率-电流协调控制方法
WO2024007465A1 (zh) 双馈型变速抽水蓄能机组自适应变参数调频方法及装置
CN110380432B (zh) 一种直驱风电场次同步振荡抑制方法及其系统
CN110739721A (zh) 一种电压源型风电机组控制方法及系统
CN110518631B (zh) 一种直驱风电机组的稳定性评估方法及系统
CN113113929B (zh) 柔性直流输电系统的电网构造型控制方法、装置及介质
CN110350587B (zh) 基于换流器控制的双馈风机次同步振荡抑制方法及其系统
CN110311414B (zh) 基于能量稳定域的双馈风机控制参数优化方法及系统
WO2023077740A1 (zh) 一种并网变流器直流电容同步控制系统及方法
CN109787274B (zh) 一种虚拟同步控制方法及转子侧变频器控制器
CN108964095B (zh) 一种抑制风电场低频振荡的储能控制方法及系统
CN113991705A (zh) 一种含储能风电场参与一次调频控制的风电并网系统调频控制方法
Hemeyine et al. Robust takagi sugeno fuzzy models control for a variable speed wind turbine based a DFI-generator
Xue et al. Influence and comparison of P/Q-control based VSC-HVDC system on the grid power oscillation damping
Li et al. VSG virtual inertial control strategy based on lead-lag link and fuzzy logic control
CN114243762A (zh) 一种风机并网的分析与控制方法
CN111769574B (zh) 基于能量耗散的直驱风场次/超频振荡抑制方法及装置
CN109713664B (zh) 直流孤岛频率稳定的网源协调控制策略计算方法及系统
CN109787267B (zh) 一种双馈风机经vsc-hvdc接入弱受端系统的判定方法
Boroujeni et al. Dynamic stability enhancement of a multimachine electric power system using STATCOM

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