CN105762816B - 一种抑制次同步相互作用的方法 - Google Patents

一种抑制次同步相互作用的方法 Download PDF

Info

Publication number
CN105762816B
CN105762816B CN201610187738.2A CN201610187738A CN105762816B CN 105762816 B CN105762816 B CN 105762816B CN 201610187738 A CN201610187738 A CN 201610187738A CN 105762816 B CN105762816 B CN 105762816B
Authority
CN
China
Prior art keywords
damping controller
double
characteristic root
probability
time constant
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
CN201610187738.2A
Other languages
English (en)
Other versions
CN105762816A (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.)
Shanghai University of Electric Power
Original Assignee
Shanghai University of Electric Power
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 Shanghai University of Electric Power filed Critical Shanghai University of Electric Power
Priority to CN201610187738.2A priority Critical patent/CN105762816B/zh
Publication of CN105762816A publication Critical patent/CN105762816A/zh
Application granted granted Critical
Publication of CN105762816B publication Critical patent/CN105762816B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • H02J3/386
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明涉及一种抑制次同步相互作用的方法,考虑双馈风电场串联补偿系统的多运行状态,该方法包括以下步骤:(1)将节点电压、节点注入以及系统状态矩阵特征根视作随机变量,并采集随机变量的样本数据,获得概率特征;(2)搭建双馈风电场串联补偿系统的小干扰模型,得到系统线性化后的状态矩阵;(3)选取阻尼控制器合适的初值,所述的初值包括增益、时间常数、超前时间常数和滞后时间常数;(4)求得步骤(2)中状态矩阵的特征根;(5)对阻尼控制器的接入位置进行初步定位;(6)对阻尼控制器进行精确选址。与现有技术相比,本发明具有节约成本、效果好等优点。

Description

一种抑制次同步相互作用的方法
技术领域
本发明属于电力系统次同步相互作用技术领域,尤其是涉及一种抑制次同步相互作用的方法。
背景技术
近些年,风电作为清洁、可再生能源得到迅猛发展,风力发电机组的装机容量不断增加,对风电大规模并网及消纳提出了更高的要求。因风力资源更为丰富的地方一般远离负荷中心,需经过远距离输送,串联补偿电容技术可减小线路损耗,提高电能传输容量,有效并且经济地实现远距离输电,是目前实现大规模风电场功率外送的主要措施。与此同时也出现了不少的问题,其中第一点是串联电容补偿线路会与双馈风电机组变频器产生交互作用,引起一种抑制次同步相互作用的方法新的次同步振荡问题——次同步控制相互作用。随着双馈感应风力发电机越来越普遍的使用,由串补线路引发的次同步问题威胁到了大规模风电场的电力外送问题。第二点:当发生次同步控制相互作用时,如何抑制也是一个重要方面。目前大多数学者研究的都是风电场经串补线路并网的次同步谐振问题,对次同步相互作用问题的探讨比较少。第三点:目前对次同步控制相互作用的抑制方法主要是加装FACTs装置,但因为设备的成本昂贵,会增加投资费用,也有文献提出加装阻尼控制器抑制次同步控制相互作用,并验证了此方法的有效性,但目前所作的研究都是针对某种或是比较几种特定的运行方式,而忽略了风电并网系统中风电出力的随机性。因此,有必要研究一种抑制次同步相互作用的方法新的设计方法,使得所设计的阻尼控制器能够满足系统的多运行方式。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种节约成本、效果好、适用范围广的抑制次同步相互作用的方法。
本发明的目的可以通过以下技术方案来实现:一种抑制次同步相互作用的方法,在双馈风电场串联补偿系统中接入阻尼控制器,抑制由风机变频器和串联电容补偿线路引起的次同步相互作用,考虑双馈风电场串联补偿系统的多运行状态,包括负荷、常规发电机组出力和风力发电机组出力的随机变化,该方法包括以下步骤:
(1)将节点电压、节点注入以及系统状态矩阵特征根视作随机变量,并采集随机变量的样本数据,获得概率特征,并利用累积概率密度曲线表示;
(2)搭建双馈风电场串联补偿系统的小干扰模型,得到整个电力系统线性化后的状态矩阵;
(3)选取阻尼控制器合适的初值,所述的初值包括增益、时间常数、超前时间常数和滞后时间常数;
(4)求得步骤(2)中状态矩阵的特征根和阻尼比概率属性;
(5)对阻尼控制器的接入位置进行初步定位;
(6)对阻尼控制器进行精确选址。
所述的步骤(2)中采用插入式建模技术PMT搭建双馈风电场串联补偿系统的小干扰模型。
所述的步骤(3)中利用概率灵敏度的方法选取阻尼控制器合适的初值。
所述的步骤(4)具体为:利用概率法结合李雅普诺夫第一法求得状态矩阵的特征根和阻尼比概率属性,采用混合中心矩与半不变量结合Gram-Charlier级数展开的方法得到次同步特征根实部和阻尼比的概率密度曲线;所述的次同步特征根为特征根虚部大于(20×π)的特征根,其中,×表示乘法运算;
所述的步骤(5)具体为:计算次同步振荡模式的参与因子并将大于0.2的参与因子对应的状态变量作为与次同步振荡模式强相关的状态变量,并将所述的状态变量在变频器中的对应位置作为阻尼控制器的初步定位位置,所述的参与因子plk表示第l个状态变量对第k个模式的参与程度,其中,u,v分别表示左特征向量和右特征向量,n表示状态变量的总数。
所述的步骤(6)具体为:利用二阶概率灵敏度对阻尼控制器精确选址,并比较功率反馈型、电压反馈型、电流反馈型和速度反馈型四种输入信号对阻尼控制器阻尼效果的影响,使阻尼控制器具有概率环境下的鲁棒性,满足多运行方式下稳定次同步振荡模式的要求。
由阻尼系数和阻尼比共同反映多运行方式下系统的稳定程度。当同时满足式(1)和(2)时,表征系统符合鲁棒稳定性要求,保证系统的动态品质。
其中,α'k和ξ'k分别表示扩展的阻尼系数和阻尼比,αc和ξc分别表示阻尼系数和阻尼比的边界值。
与现有技术相比,本发明具有以下优点:
(1)阻尼控制器的加装提高了风电系统的次同步振荡稳定性,在保证次同步控制相互作用抑制效果的前提下,减少了设备投资成本;
(2)采用概率特征根法得出的结论更具普适性,并能够为次同步控制相互作用阻尼控制器方案的设计提供有价值的信息。
附图说明
图1为利用参与因子和概率灵敏度指标,在转子侧变频器选择阻尼控制器的接入位置;
图2为利用参与因子和概率灵敏度指标,在网侧变频器选择阻尼控制器的接入位置;
图3为负荷的累积概率密度曲线;
图4为风电场G5和同步发电机G1、G2、G4有功功率累积概率密度曲线;
图5为同步发电机G1~G4电压(G3作为参考电机)和风电场G5电压的累积概率密度曲线;
图6为加装阻尼控制器前后特征根实部小于0的概率密度曲线。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
一种抑制次同步相互作用的方法,在双馈风电场串联补偿系统中接入阻尼控制器,抑制由风机变频器和串联电容补偿线路引起的次同步相互作用,考虑双馈风电场串联补偿系统的多运行状态,包括负荷、常规发电机组出力和风力发电机组出力的随机变化,该方法包括以下步骤:
(1)将节点电压、节点注入以及系统状态矩阵特征根视作随机变量,并采集随机变量的样本数据,获得概率特征,并利用累积概率密度曲线表示;
(2)采用插入式建模技术PMT搭建双馈风电场串联补偿系统的小干扰模型,得到整个电力系统线性化后的状态矩阵;
(3)利用概率灵敏度的方法选取阻尼控制器合适的初值,初值包括增益、时间常数、超前时间常数和滞后时间常数;
(4)利用概率法结合李雅普诺夫第一法求得状态矩阵的特征根和阻尼比概率属性,采用混合中心矩与半不变量结合Gram-Charlier级数展开的方法得到次同步特征根实部和阻尼比的概率密度曲线;次同步特征根为特征根虚部大于(20×π)的特征根,其中,×表示乘法运算;
(5)计算次同步振荡模式的参与因子并将大于0.2的参与因子对应的状态变量作为与次同步振荡模式强相关的状态变量,并将所述的状态变量在变频器中的对应位置作为阻尼控制器的初步定位位置,所述的参与因子plk表示第l个状态变量对第k个模式的参与程度,其中,u,v分别表示左特征向量和右特征向量,n表示状态变量的总数。
(6)利用二阶概率灵敏度对阻尼控制器精确选址,并比较功率反馈型、电压反馈型、电流反馈型和速度反馈型四种输入信号对阻尼控制器阻尼效果的影响,使阻尼控制器具有概率环境下的鲁棒性,满足多运行方式下稳定次同步振荡模式的要求。
本实施例首先分析双馈风电场并网系统次同步控制相互作用:
利用PMT对含双馈风电机组的五机两区域系统搭建小扰动模型,图3、图4及图5中利用累积概率密度函数表示风电场出力、同步发电机出力、负荷等随机波动因素,求得概率特征根后分析次同步控制相互作用振荡模式的概率特征,根据表1中列出双馈风机并网新增的4个振荡模式分析可得:为特征根的期望,特征根实部反映阻尼的大小,特征根实部越远离虚轴,表明系统稳定性越强。特征根虚部表征振荡频率,由计算公式f=ω/2π得:模式1的振荡频率为14.89,属于次同步振荡范围。为阻尼比的期望,其值越大,表征动态衰减特性越好。分析系统多运行方式下的次同步振荡模式的小扰动稳定性,通常主要关注的是特征根实部的上界和阻尼比的下界。次同步模式的分别为45.62%和47.11%,稳定性不够。
根据参与因子表2可知:次同步控制相互作用与四台同步发电机无关,是由双馈风机经串补线路并网引起的。次同步振荡模式中ωt和θt的参与因子非常小,轴系状态变量不参与,说明次同步振荡是由电气量的振荡,并且是由风电机组的变频器引起的,可以判断为次同步控制相互作用。
由表3中的概率灵敏度指标可知:转子侧变频器中的参数Kp2,以及网侧变频器中的参数Kp7对次同步振荡模式影响较大,两者的符号都为负,当减小Kp2和Kp7时,次同步特征根向远离虚轴的方向移动。次同步振荡模式对参数Kp5有很高的灵敏度,并且符号相反,减小Kp5会使特征根向虚轴靠近。
由阻尼系数和阻尼比共同反映多运行方式下系统的稳定程度。当同时满足式(1)和(2)时,表征系统符合鲁棒稳定性要求,保证系统的动态品质。
其中,α'k和ξ'k分别表示扩展的阻尼系数和阻尼比,αc和ξc分别表示阻尼系数和阻尼比的边界值,为保证特征根实部的概率分布位于左半平面,本实施例中,αc取0,为保证系统的动态特性,阻尼比的概率密度曲线位于ξ=ξc的右侧,ξc取0.1。
由ak和特征根的虚部ωk的各阶半不变量计算得出阻尼比ξk的各阶半不变量γj(ξ),经过Gram-Charlier展开级数得出ak和ξk的概率密度f(αk)和f(ξk),通过式(3)和式(4)确定ak和ξk的分布概率。
其中,表示概率灵敏度指标,相对于任一参数К的灵敏度表明该参数对α'k和ξ'k调节的有效程度,也表明正确的调节方向。当参数К选择变频器中的PI参数时,式(5)和(6)中的反映了变频器的增益和时间常数在多系统运行方式下对特征根λk的综合影响,从而可用于选择阻尼控制器的接入位置。在此基础上,比较功率反馈型、电压反馈型、电流反馈型和速度反馈型四种输入信号,并进行选择,从而为系统提供合适的阻尼,改善次同步控制相互作用的抑制效果。
本实施例搭建了包含双馈风电场和阻尼控制器的5机2区域系统(系统中有5台同步风力发电机G1~G5),并利用Fortran语言编写程序进行仿真,得到风电场经串联电容补偿线路并网后引起的次同步控制相互作用振荡模式的特征概率,验证了本发明所提出的利用概率理论设计阻尼控制器抑制次同步控制相互作用的有效性和准确性。
分析加装阻尼控制器抑制次同步控制相互作用:
根据参与因子对阻尼控制器的接入位置进行初步定位,二阶概率灵敏度对阻尼控制器精确选址,在风机变频器的3个位置加装阻尼控制器,定子电磁功率作为统一的输入信号,阻尼控制器的初值通过概率灵敏度方法确定,分别取值为:KDC=0.001,Tw=10,T1=10.3,T2=0.5。将阻尼控制其加装在转子侧变频器有功电流控制环、网侧变频器有功功率控制环、网侧变频器无功电流控制环,如图1和图2所示,图1中idr和iqr分别表示转子电流的dq轴分量,id_ref和iq_ref分别表示转子电流控制参考值,vs和vref分别表示发电机出口处的电压及其参考值,Ps和Pref分别表示风机有功功率及其参考值,(x1)、(x2)、(x3)及(x4)均表示PI控制器中的状态变量;图2中vDC和vDC_ref分别表示直流环节并联电容器两端电压及其控制参考值,idg和iqg分别表示电网侧变频器电流的d轴和q轴分量,idg_ref和iqg_ref分别表示电网侧变频器电流的d轴和q轴分量控制参考值,(x5)、(x6)及(x7)均表示PI控制器中的状态变量。表1为风电场接入电力系统后,五机两区域系统新增的特征根概率;表2为次同步振荡模式的参与因子(列出部分有具有代表性的值);表3为次同步振荡模式对变频器PI控制参数和状态变量的灵敏度指标;表4为在转子侧变频器有功功率外环加装阻尼控制器后的特征根概率(定子电磁功率Ps设为输入信号);表5为在网侧变频器有功功率控制环加装阻尼控制器后的特征根概率(定子电磁功率Ps设为输入信号);表6为在网侧变频器无功电流控制环加装阻尼控制器后的特征根概率(定子电磁功率Ps设为输入信号);表7为联络线电流Il作为输入信号,阻尼控制器接入无功电流控制环时的特征根概率;表8为双馈感应风力发电机转速ω作为输入信号,阻尼控制器接入无功电流控制环时的特征根概率;表9为直流电容电压Vdc作为输入信号,阻尼控制器接入无功电流控制环时的特征根概率。
表4、表5及表6分别为得到的特征根概率。对比所得结果,选择网侧变频器无功电流控制环为合适的接入位置。
在此基础上,除了定子电磁功率,另将联络线电流、双馈感应风力发电机的转速、直流电容电压作为输入信号,得到的特征根概率如表7、表8及表9所示。对比分析可知:直流电容电压Vdc最适合作为输入信号。
因此,双馈风电场在加入本发明所设计方法的阻尼控制器之后,概率特征根向远离虚轴的方向移动,概率阻尼比从表1中0.141增加到0.718,复特征根实部小于0的概率即系统稳定的概率增加到了90.38%,阻尼比大于0.1的概率达到了100%,加装阻尼控制器前后特征根实部小于0的概率密度函数对比如图6所示,可明显看出阻尼控制器的加装提高了风电系统的次同步振荡稳定性,在保证次同步控制相互作用抑制效果的前提下,减少了设备投资成本。
表1
表2
表3
表4
表5
表6
表7
表8
表9

Claims (5)

1.一种抑制次同步相互作用的方法,在双馈风电场串联补偿系统中接入阻尼控制器,抑制由风机变频器和串联电容补偿线路引起的次同步相互作用,其特征在于,考虑双馈风电场串联补偿系统的多运行状态,该方法包括以下步骤:
(1)将节点电压、节点注入以及系统状态矩阵特征根视作随机变量,并采集随机变量的样本数据,获得概率特征;
(2)搭建双馈风电场串联补偿系统的小干扰模型,得到系统线性化后的状态矩阵;
(3)选取阻尼控制器合适的初值,所述的初值包括增益、时间常数、超前时间常数和滞后时间常数;
(4)求得步骤(2)中状态矩阵的特征根;
(5)对阻尼控制器的接入位置进行初步定位;
(6)对阻尼控制器进行精确选址;
所述的步骤(5)具体为:计算次同步振荡模式的参与因子并将大于0.2的参与因子对应的状态变量作为与次同步振荡模式强相关的状态变量,并将所述的状态变量在变频器中的对应位置作为阻尼控制器的初步定位位置,所述的参与因子plk表示第l个状态变量对第k个模式的参与程度,其中,u,v分别表示左特征向量和右特征向量,n表示状态变量的总数。
2.根据权利要求1所述的一种抑制次同步相互作用的方法,其特征在于,所述的步骤(2)中采用插入式建模技术PMT搭建双馈风电场串联补偿系统的小干扰模型。
3.根据权利要求1所述的一种抑制次同步相互作用的方法,其特征在于,所述的步骤(3)中利用概率灵敏度的方法选取阻尼控制器合适的初值。
4.根据权利要求1所述的一种抑制次同步相互作用的方法,其特征在于,所述的步骤(4)具体为:利用概率法结合李雅普诺夫第一法求得状态矩阵的特征根。
5.根据权利要求1所述的一种抑制次同步相互作用的方法,其特征在于,所述的步骤(6)具体为:利用二阶概率灵敏度对阻尼控制器精确选址。
CN201610187738.2A 2016-03-29 2016-03-29 一种抑制次同步相互作用的方法 Active CN105762816B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610187738.2A CN105762816B (zh) 2016-03-29 2016-03-29 一种抑制次同步相互作用的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610187738.2A CN105762816B (zh) 2016-03-29 2016-03-29 一种抑制次同步相互作用的方法

Publications (2)

Publication Number Publication Date
CN105762816A CN105762816A (zh) 2016-07-13
CN105762816B true CN105762816B (zh) 2018-03-30

Family

ID=56346629

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610187738.2A Active CN105762816B (zh) 2016-03-29 2016-03-29 一种抑制次同步相互作用的方法

Country Status (1)

Country Link
CN (1) CN105762816B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106972511B (zh) * 2017-04-25 2020-12-01 国网宁夏电力公司电力科学研究院 基于网侧控制回路优化的直驱风电次同步振荡抑制方法
CN108631331B (zh) * 2018-04-24 2020-07-17 华北电力科学研究院有限责任公司 一种双馈风电场次同步振荡抑制方法及装置
CN108599236B (zh) * 2018-04-24 2020-08-14 华北电力科学研究院有限责任公司 双馈风电场次同步振荡svg抑制方法及装置
CN108631338B (zh) * 2018-06-21 2021-10-22 国网宁夏电力有限公司电力科学研究院 一种用于抑制双馈风电场并网次同步振荡的方法
CN109361226B (zh) * 2018-11-05 2020-08-28 上海交通大学 次同步阻尼控制器反馈信号和安装位置的选取方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006090538A1 (ja) * 2005-02-22 2006-08-31 Kyushu Institute Of Technology 電力系統の系統安定度制御方法及びシステム
CN103729570B (zh) * 2014-01-21 2017-05-31 山东大学 基于矩阵摄动理论的电力系统振荡模式的匹配方法

Also Published As

Publication number Publication date
CN105762816A (zh) 2016-07-13

Similar Documents

Publication Publication Date Title
CN105762816B (zh) 一种抑制次同步相互作用的方法
CN102012956B (zh) 一种考虑风电场输入风速和风向随机波动的风电场等值方法
CN102856917B (zh) 一种配电网无功优化方法
CN108011364B (zh) 一种分析dfig动力学特性与电力系统动态交互影响的方法
Sun et al. A review on analysis and control of small signal stability of power systems with large scale integration of wind power
CN109347097B (zh) 一种基于改进粒子群优化算法的双馈风电系统次同步振荡抑制策略
CN106059422B (zh) 一种用于双馈风电场次同步振荡抑制的模糊控制方法
CN103558768B (zh) 一种基于风电场内风速分布特性的等值建模方法
CN112448399B (zh) 基于模拟电感的双馈风电场次同步振荡的抑制方法
CN109522607A (zh) 一种双馈风电场机电暂态等值建模方法
CN110206686A (zh) 一种用于风力发电机组的自适应最大功率跟踪控制方法
CN103094920A (zh) 一种直驱式风电机组风电场等值方法
CN107204634A (zh) 分散式风电的无功电压控制方法和系统
CN111009921A (zh) 基于奈奎斯特稳定判据的双馈风机并网系统振荡分析方法
CN104617578B (zh) 一种含风电场电力系统的可用输电能力的获取方法
CN107069799A (zh) 双馈式风力发电机组的控制方法和系统
CN106712032B (zh) 一种考虑风电机组有功电压调节能力的最优潮流模型构建方法
Wu et al. Development of an equivalent wind farm model for frequency regulation
CN114665506A (zh) 一种大规模风电并网情况下电力系统惯性常数的测量方法
Xu et al. Sub-synchronous frequency domain-equivalent modeling for wind farms based on rotor equivalent resistance characteristics
Yan et al. Analysis of subsynchronous oscillation of direct drive PMSG based wind farm under low operating condition
CN103956767B (zh) 一种考虑尾流效应的风电场并网稳定性分析方法
CN110571794B (zh) 一种适用于双馈风机风电场的暂态模型等值计算方法
CN108493988A (zh) 一种包含风电机组的电网惯性时间常数计算方法及系统
CN112271724A (zh) 一种基于电压调控的虚拟电厂分区构建模型及构建方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant