CN109474014B - 一种对双馈风场接入电网友好性的量化评估方法 - Google Patents
一种对双馈风场接入电网友好性的量化评估方法 Download PDFInfo
- Publication number
- CN109474014B CN109474014B CN201811526140.7A CN201811526140A CN109474014B CN 109474014 B CN109474014 B CN 109474014B CN 201811526140 A CN201811526140 A CN 201811526140A CN 109474014 B CN109474014 B CN 109474014B
- Authority
- CN
- China
- Prior art keywords
- input
- local
- output
- gain
- integral
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000011158 quantitative evaluation Methods 0.000 title claims abstract description 22
- 230000009977 dual effect Effects 0.000 title description 2
- 238000004088 simulation Methods 0.000 claims abstract description 20
- 238000004458 analytical method Methods 0.000 claims description 22
- 230000004044 response Effects 0.000 claims description 14
- 238000011217 control strategy Methods 0.000 claims description 12
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 230000001360 synchronised effect Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000013178 mathematical model Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 3
- 230000006641 stabilisation Effects 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 230000006698 induction Effects 0.000 claims description 2
- 238000011160 research Methods 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 4
- 238000012795 verification Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Images
Classifications
-
- H02J3/386—
-
- 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]
-
- 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
- Supply And Distribution Of Alternating Current (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
Description
技术领域
本发明属于电力系统控制技术领域,具体涉及一种对双馈风场接入电网友好性的量化评估方法。
背景技术
随着全球能源短缺和环境污染问题日益严重,大力发展可再生能源已成为全世界的共识,其中风能是最重要的可再生能源之一,近年来已经得到了人们的极大关注。与此同时,由于风电具有间歇性,为了提高风力发电对电网的友好性,需要有效评估风力发电并网对电力系统的影响。
国内外现有对于DFIG并网稳定性的分析研究大多集中在功角稳定问题。由于DFIG对电网电压变化敏感,在系统故障时可能会使系统的动态响应性能恶化,甚至引起电压失稳。因此,研究DFIG并网对互联电力系统的短期电压稳定性的影响具有重要意义。
时域仿真法广泛用于风电场的稳定性分析,虽然时域仿真可以对系统进行精确的稳定性评估,但是计算量较大,并且通过仿真无法获得系统稳定裕度的量化分析结果。目前,Lyapunov理论已广泛应用于电力系统的稳定性分析,通过计算初始状态到给定稳定区域边界的距离可以得到量化的稳定性分析结果。然而,缺乏用于构造Lyapunov函数的系统算法使得该理论难以应用于一般非线性系统,并且Lyapunov理论研究无外部扰动的系统的稳定性,无法量化分析外部扰动对系统的影响。
研究非线性系统稳定性的另一种途径是使用输入-状态稳定性(ISS)理论,ISS理论将Lyapunov方法和输入-输出稳定理论相结合,将互联系统分解为若干子系统,然后利用子系统的ISS特性及其连接关系研究互联系统的稳定性。ISS理论可以研究系统在外部干扰下的动态响应特性,并量化系统可承受的外部干扰的大小,因此,ISS理论适用于研究含风力发电的电力系统的稳定性问题。
应用ISS理论的关键问题之一是估计子系统的稳定属性,已有研究主要是通过人工观察法或基于系统的数学模型进行分析。然而,实际系统难以精确建模,简化又可能导致不准确的分析结果。
发明内容
为克服现有技术中的问题,本发明的目的在于提出一种对双馈风场接入电网友好性的量化评估方法。
为实现上述目的,本发明采用如下的技术方案:
一种对双馈风场接入电网友好性的量化评估方法,包括以下步骤:
S1,基于输入-状态稳定的等价定义形式积分-积分估计,采用仿真的方法预先估计子系统的局部输入-状态稳定/局部输入-输出稳定属性,其中局部输入-状态稳定/局部输入-输出稳定属性包括最大输入积分能量、局部稳定域、和渐近增益,渐近增益采用分段线性形式进行估计;
S2,根据互联系统稳定性量化评估判据,计算获取和双馈风场互联后电力系统的小增益条件和局部稳定极限,量化评估双馈风场接入对电网的友好性。
本发明进一步的改进在于,
(1)子系统层面:基于积分-积分估计,根据不同控制策略对双馈风场输出特性性的量化评估;
(2)互联系统层面:基于小增益定理和积分-积分估计,对双馈风场接入电网友好性的量化评估。
本发明进一步的改进在于,具有外部输入的非线性系统:
y=h(x,u)
其中x∈Rn,y∈Rp,u∈Rm,f:Rn×m→Rn,h:Rn×m→Rp,f和h在x和u中是连续的局部Lipschitz,f(0,0)=0,h(0,0)=0;
LISS的原始定义如下所示:
则上述非线性系统是ISS的,欧几里德范数由|·|表示,γISS称为输入到状态渐近增益;||u||∞表示外部输入u的上确界,它是使得|u(t)|≤a对所有t成立的最小的a;比较函数的定义如下:
定义2.对于任何x0∈ΩIOS∈Rn,u∈UIOS∈Rm;如果存在函数
x0∈ΩISS∈Rn,u∈UISS∈Rm,下式成立:
则上述非线性系统称为局部输入输出稳定系统;
定理1.在非线性系统的LISS/LIOS分析中,ΩISS,UISS,和ΩIOS,UIOS分别表示上述非线性系统的LISS域和LIOS域;在实际系统中,上述非线性系统的稳定域用Ω和U来表示,其中Ω=ΩISS∩ΩIOS,U=UISS∩UIOS;
定义3.对于任意x0∈ΩISS∈Rn,u∈UISS∈Rm,如果下列不等式成立
同样,局部输入-输出积分-积分估计的定义为:
由上述积分-积分估计的定义看出,系统状态、外部输入和系统输出均由积分表示,包含外部扰动的幅值信息以及随时间波动的信息;
基于定义1-3及定理1,采用电力系统系统LISS/LIOS分析的实用算法,利用系统初始状态和外部输入的局部定义域和渐近增益对非线性系统的动态响应进行量化评估。
本发明进一步的改进在于,利用系统初始状态和外部输入的局部定义域和渐近增益对非线性系统的动态响应进行量化评估的具体过程为:
设系统初始状态位于平衡点,当考虑线性增益函数时,渐近增益表示如下:
其中,T为步长;
通过改变输入信号u的幅值,得到一系列的渐近增益,选择整个输入域上的最大γ作为系统线性形式的渐近增益将导致估计结果非常保守,为了降低LISS/LIOS分析的保守性,采用分段线性渐近增益。
本发明进一步的改进在于,分段线性渐近增益的估计过程如下:
第一步,根据式(1)估计区间[0,a1]上的线性渐近增益函数w=l1υ;
第二步,令输入u∈[a1,a2],根据第一步计算得到的线性渐近增益函数的端点[a1,b1]进行估计;
d2=b1-a1l2
其中b1=l1a1·
第三步,剩余输入区间的渐近增益估计如下:
di=bi-1-ai-1li
其中bi-1=li-1ai-1+di-1。
minz=g(m1)+w(m2)-m3
其中g是预先给定的函数,以近似估计w是渐近增益函数,C是仿真得到的点(m1,m2,m3)的集合;选择外部输入的局部区域和初始状态的局部区域作为和|x0|的最大值,使系统状态和输出保持在最大容许范围内,最大容许范围分别为xlim和ylim。
本发明进一步的改进在于,对于一个由n个子系统组成的互联系统,其中第i个子系统数学模型表达式如下:
基于积分-积分估计的互联系统量化稳定性判据如下:
如果满足下列四个条件,则上述n个子系统组成的互联系统是LISS和LIOS的,说明对双馈风场接入电网友好:
II.每个子系统的输入/输出满足以下条件:
|u(s)|≤Z|y(s)|
其中Zij≥0,当选择电压和电流作为输入输出时,由系统的导纳矩阵导出Z;
III.以下小增益条件成立:
G(υ*)<<υ*
IV.系统输出满足预先给定如下条件:
与现有技术相比,本发明具有的有益效果:
本发明首先基于离线仿真,提出了一种分析子系统LISS特性的实用算法,实现了不同控制策略对双馈风场输出特性影响的量化评估。其次,基于小增益定理,提出了双馈风场并网的量化稳定性评估指标。该方法的优点是:首先,基于所提子系统层面的量化评估指标,可以得到不同控制策略对双馈风场输出特性影响的量化评估。其次,通过提出的双馈风场并网的量化稳定性指标可以量化评估双馈风场并网后系统的动态响应,进而为系统的控制提供依据,本发明通过时域仿真验证了所提出的方法正确性。通过本发明的对双馈风场接入电网友好性的量化评估方法,能够研究双馈风电场并网对电力系统短期电压定量稳定性的影响。
附图说明
图1为不同DFIG控制策略的γ拟合。
图2为各控制策略的详细仿真结果和相应的γ拟合图,其中,(a)为PI控制,(b)为磁通阻尼控制,(c)为改进的前馈控制,(d)为基于SDRE的非线性控制。
图3为电力系统结构图。
图4为DFIG渐近增益的估计。
图5为初始状态、外部输入和系统输出的关系。
图6为小增益条件的仿真验证。
图7为风电场端电压的动态变化。
图9为n=16时风电场端电压的动态变化。
图10为本发明的流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了分析双馈风机并网后电力系统的稳定性,需要研究互联系统在外部干扰下系统的动态响应特性,并量化互联系统可承受的外部干扰的大小,量化评估双馈风场接入电网的友好性,通过采用分段线性形式的渐近增益,可以减少子系统LISS分析的保守性。
参见图10,本发明的一种双馈风电场并网友好性量化评估方法包括:
S1,采用输入-状态稳定(Input-to-StateStability,ISS)的等价定义形式积分-积分估计,采用仿真的方法预先估计子系统的局部输入-状态稳定/局部输入-输出稳定(LISS/LIOS)属性,其中LISS/LIOS属性包括最大输入积分能量、局部稳定域、和渐近增益,渐近增益采用分段线性形式进行估计。
S2,提出互联系统稳定性量化评估判据,计算获取和双馈风场互联后电力系统的小增益条件和局部稳定极限,量化评估双馈风场接入对电网的友好性;
与DFIG互联的电力系统包括双馈风机、同步发电机、感应电动机、恒功率负载和恒阻抗负载。
在进行S1之前,本专利先给出以下定义以及推论:
考虑以下具有外部输入的非线性系统:
y=h(x,u)
其中x∈Rn,y∈Rp,u∈Rm,f:Rn×m→Rn,h:Rn×m→Rp,f和h在x和u中是连续的局部Lipschitz,f(0,0)=0,h(0,0)=0。
LISS的原始定义如下所示:
则上述非线性系统是ISS的,欧几里德范数由|·|表示,γISS称为输入到状态渐近增益。||u||∞表示外部输入u的(本质)上确界,它是使得|u(t)|≤a对所有t成立的最小的a,比较函数的定义如下:
定义2.对于任何x0∈ΩIOS∈Rn,u∈UIOS∈Rm。如果存在函数
x0∈ΩISS∈Rn,u∈UISS∈Rm,以下成立:
则上述非线性系统称为局部输入输出稳定(LIOS)系统。
定理1.在非线性系统的LISS/LIOS分析中,初始状态的局部定义域和外部输入的局部定义域是两个关键因素,如上述定义所示,ΩISS,UISS,和ΩIOS,UIOS分别表示上述非线性系统的LISS域和LIOS域。在实际系统中,上述非线性系统的稳定域可以用Ω和U来表示,其中Ω=ΩISS∩ΩIOS,U=UISS∩UIOS。
渐近增益反映了外部扰动对系统状态和输出的影响,并可量化系统对外部扰动的承受能力。基于上述定义的外部输入由其上确界表示,仅能反映外部扰动的幅值信息,不能反映其随时间的波动信息。为克服这一问题,本发明采用ISS等价定义形式——积分-积分估计。
定义3.对于任意x0∈ΩISS∈Rn,u∈UISS∈Rm。如果下列不等式成立
同样,本发明给出局部输入-输出积分-积分估计的定义:
由上述积分-积分估计的定义可以看出,系统状态、外部输入和系统输出均由积分表示的,包含外部扰动的幅值信息以及随时间波动的信息。
基于上述定义及定理,本发明提出一种电力系统系统LISS/LIOS分析的实用算法,利用系统初始状态和外部输入的局部定义域和渐近增益对非线性系统的动态响应进行量化评估。
设系统初始状态位于平衡点,当考虑线性增益函数时,渐近增益可以表示如下:
其中,T为步长,T是由实际需要决定的。
通过改变输入信号u的幅值,可以得到一系列的渐近增益,选择整个输入域上的最大γ作为系统线性形式的渐近增益将导致估计结果非常保守。为了降低LISS/LIOS分析的保守性,本发明采用分段线性形式的渐近增益。
取LIOS以分析为例,选择如下分段线性渐近增益:
分段线性渐近增益的详细的估计过程如下:
第一步,根据式(1)估计区间[0,a1]上的渐近增益函数w=l1υ。
第二步,令输入u∈[a1,a2],根据第一步计算得到的线性渐近增益函数的端点[a1,b1]进行估计。
d2=b1-a1l2
其中b1=l1a1.
第三步,与第二步类似,剩余输入区间的渐近增益可以估计如下。
di=bi-1-ai-1li
其中bi-1=li-1ai-1+di-1.
为了进一步估计外部输入的局部区域和初始状态的局部区域,分别选择系统状态和输出的上限为xlim和ylim。从定义3可以看出,函数与非线性系统的初始状态有关,估计需要考虑|x0|的变化。设m1=|x0|,改变m1和m2的值,可以得到一系列的点(m1,m2,m3)。函数的估计可以表示为下面的优化问题:
minz=g(m1)+w(m2)-m3
其中g是预先给定的函数,以近似估计w是渐近增益函数,C是仿真得到的点(m1,m2,m3)的集合。从C中仿真结果可以看出,选择外部输入的局部区域和初始状态的局区部域作为和|x0|的最大值,使系统状态和输出保持在最大容许范围内,最大容许范围分别为xlim和ylim。
对于子系统层面:基于积分-积分估计,根据不同控制策略对双馈风场输出特性性的量化评估;
对于互联系统层面:基于小增益定理和积分-积分估计,对双馈风场接入电网友好性的量化评估。
本发明基于子系统的LISS/LI0S特性及其连接,提出互联后系统的量化稳定性判据,本发明还提出了一种适用于积分-积分估计的实用性算法以校验所提稳定性判据。
考虑一个由n个子系统组成的互联系统,其中第i个子系统数学模型表达式如下:
本发明给出基于积分-积分估计的互联系统量化稳定性判据如下:
如果满足下列四个条件,则上述n个子系统组成的互联系统是LISS和LIOS的,则说明对双馈风场接入电网友好:
II.每个子系统的输入/输出满足以下条件:
|u(s)|≤Z|y(s)|
其中Zii≥0,当选择电压和电流作为输入输出时,可以由系统的导纳矩阵导出Z。
III.以下小增益条件成立:
G(υ*)<<υ*
IV.系统输出满足预先给定如下条件:
上述稳定性判据中的第四个条件可通过计算复合函数的值来进行校验的。因此,还需验证条件III中的局部小增益条件。
本发明给出了局部小增益条件的等价形式,详细的描述如下:
引理1:设不动点满足G(υ*)<<υ*,考虑一系列离散序列{υ(k)}满足υ(k+1)=G(υ(k)),其中k=0,1,2...,υ(0)=υ*。当且仅当υ(k)→0为k→∞时局部小增益条件III成立。
从引理1可以知道,需要对序列{υ(k)}的收敛性进行验证。由于本发明使用积分-积分估计,输入能量不能完全收敛到零。从实际应用的角度看,如果系统的状态和输出收敛到平衡点附近区域,则该系统被认为是稳定的。设系统存在保持稳定的最小输入能量υ0,则局部小增益的验证条件可以转换为当k→∞时候υ(k)→υ0。
引理1中不动点υ*可由各子系统的输入所组成的向量表示。一般地,系统的运行约束范围总是小于系统的稳定范围,因此可以有系统的运行约束表示不动点υ*。由于研究短期电压稳定问题,考虑电压、电流的安全运行约束为Ulim和Ilim,不动点υ*可记为[Ilim,Ulim]T。
综上,仅需验证序列{υ(k)}的收敛性,其中υ(k+1)=G(υ(k)),υ(0)=υ*。具体地,给定最大迭代次数记为kmax,校验序列{υ(k)}是否能够在kmax次迭代内收敛至υ0的某一邻域内,如果收敛,则系统满足小增益条件,反之系统不稳定。
为了更好的验证小增益条件,增益函数应该是平滑的。因此,利用一个满足跟随条件的光滑非线性增益函数来逼近分段线性渐近增益函数:
下面通过两个实例对本发明实施例进行进一步说明,以下仅为本发明实施例的两个实例,本发明实施例并不以此为限。
实施例1
首先,本发明通过Matlab/Simulink对一台DFIG进行量化稳定性评估,研究不同控制对双馈风场输入输出特性的影响,DFIG的详细参数如下:
SB=1.67MVA,VB=575V,H=0.685sRs=0.023p.u.,Ls=0.18p.u.,Lm=2.9p.u.,Rr=0.016p.u.,Lr=0.16p.u.
选择端电压降作为输入信号,由于转子电流能反映风力发电机的功率控制能力,因此选择转子电流作为输出信号,转子电流极限设为Ir≤2p.u.,本发明基于所提子系统量化评估算法,分析了DFIG在常规PI控制、基于PI的改进控制和非线性控制下的输入输出特性不同策略的渐近增益估计如下:
不同控制策略的渐近增益估计如图1所示。各控制策略的详细仿真结果和相应的γ拟合如图2所示。
点v1、v2、v3和v4表示每个控制策略的最大输入。在一定的范围内,转子电流越大,在过渡过程中就会对系统产生越多的支撑,这对风力发电系统的是有利的。利用图1得出DFIG的在不同控制策略下的量化稳定性评估结果如下:
1、单纯的PI控制与基于SDRE控制相比,PI控制在暂态过程中提供了较少的控制支撑,但是转子绕组的过电流相对较低。
2、基于SDRE的非线性控制策略在暂态过程中充分考虑了DFIG的非线性特性,并且转子侧换流器具有较强的输出控制能力,从而提高了对系统无功的支撑和端电压的动态响应性能。然而,在初始故障阶段可能发生转子过电流,应考虑保护电路的动作。
实施例2:
对双馈风场并网后电力系统的稳定性进行量化评估,测试系统的总体结构如图3所示。
为了研究上述电力系统的短期电压稳定性,应合理选择每个子系统的输入和输出。同步发电机的输入和输出分别选择为电流和电压,DFIG的输入和输出分别选择为端电压和定子电流,负载的输入和输出分别选择为电压和电流。设定电力系统的安全运行约束为:在故障清除后系统内各节点电压应不小于0.8p.u,DFIG的转子电流最大限制值设为2p.u。
(1)DFIG的LISS/LIOS属性估计
DFIG的LIOS特性估算结果,如图4所示。本发明考虑一个含有20台双馈风机的风电场,其非线性渐近增益近似为其中s2=0.000683,DFIG的其他LIOS特性通过求解上述分段线性渐近增益进行估计,其中g(m1):=kDm1,优化结果为kD=0.00075。图5为当最大输出能量为σDFIG=0.0048时初始状态、外部输入和系统输出之间的总体关系。
(2)同步发电机的LISS/LIOS属性估计
基于所提方法估计同步发电机的LIOS特性,其中g(m1):=kGm1,kG>0,w(m2):=lGm2。增益函数的估计结果:kG≈10-6,lG≈0.0432,最大输入能量值为τG=0.83,相应的最大输出能量为σG=0.3972。同步发电机参数:额定电压UG=13.8kV.额定功率PG=200MVA,Xd=1.305,X′d=0.296,X″d=0.252,Xq=0.474,X″q=0.243.
(3)恒定阻抗负载的LISS/LIOS属性估计
为了研究互联的电力系统,这里的增益函数根据测试系统的基值进行归一化。基于系统安全运行约束,最大输入能量为τL=0.5,并且相应的最大输出能量为σL=0.186。恒定阻抗负载参数:ZL1=ZL2=ZL3=2.4+j1.2
(4)恒定阻抗负载的输入-输出连接关系
由于每个子系统使用电流和电压作为输入和输出,因此可以网络方程推导系统的输入-输出连接关系。详细过程描述如下。
其中,IP和UP表示公共母线的电流和电压,同步发电机的电流和电压用IG和UG表示,IL和UL表示负载和DFIG的电流和电压。消去公共节点,上式可以转化为:
其中,
综上,利用所提的稳定性量化评估方法对系统的稳定性进行分析:
考虑短路故障情况,研究故障清除后系统稳定性的量化结果,根据上面所得结果,可以得到最大输出能量为:通过在给定的输入局部域内验证局部小增益条件,图6描绘了序列υ(k)的最大值,从图6中可以看出该结果满足小增益条件,并且因此,还需验证以保证系统输出在安全约束内。
系统状态将保持在安全约束内,但是,Ψ中第二个元素接近安全约束边界。
因此,稳定性判据的第四个条件不成立。通过仿真验证上述理论分析结果,风电场由20个双馈风机组成,在t=3s时,B4发生三相接地短路故障,短路故障持续时间为0.27s。故障清除后DFIG的初始值为0.255,风电场端电压的动态响应如图7所示。从图7可以看出,在故障清除后,端电压降到0.8p.u以下,电压低于安全约束,与理论分析结果一致。
为了进一步验证所提出的定量稳定性分析方法的有效性,本发明分析了风电场中风力发电机个数对系统稳定性的影响。从上面的分析中,Ψ中第二个元素超出了极限值,而其他元素则在给定的范围中。令Ψ和中的第二个元素分别由ΨDFIG和表示,图8描述了风电场中风力发电机数n对ΨDFIG和的影响。
分析结果表明,当n=16时,系统在操作极限值范围内运行。仿真设置相同的故障条件,并给出系统的动态响应,风电场的端电压仿真结果如图9所示。从图9可以看出,将风力发电机的数量减少到n=16,风电场的端电压保持在运行状态内,验证了理论分析的正确性。
本发明通过数值仿真,提出LISS/LIOS属性的实用估计算法,降低LISS/LIOS分析的保守性。以小增益条件和局部输出条件作为双馈风场接入系统友好性的量化评估指标;通过验证输入能量收敛到平衡点附近的某一邻域内以校验友好性指标中的小增益条件;为了能够更好地验证小增益条件,渐近增益函数应该是平滑的,本发明使用平滑非线性增益函数来对分段线性渐近增益函数进行拟合。
Claims (4)
1.一种对双馈风场接入电网友好性的量化评估方法,其特征在于,包括以下步骤:
S1,基于输入-状态稳定的等价定义形式积分-积分估计,采用仿真的方法预先估计子系统的局部输入-状态稳定/局部输入-输出稳定属性,其中局部输入-状态稳定/局部输入-输出稳定属性包括最大输入积分能量、局部稳定域、和渐近增益,渐近增益采用分段线性形式进行估计;
S2,根据互联系统稳定性量化评估判据,计算获取和双馈风场互联后电力系统的小增益条件和局部稳定极限,量化评估双馈风场接入对电网的友好性;
和DFIG互联的电力系统包括双馈风机、同步发电机、感应电动机、恒功率负载和恒阻抗负载;具有外部输入的非线性系统:
y=h(x,u)
其中x∈Rn,y∈Rp,u∈Rm,f:Rn×m→Rn,h:Rn×m→Rp,f和h在x和u中是连续的局部Lipschitz,f(0,0)=0,h(0,0)=0;
局部输入-状态稳定LISS的原始定义如下所示:
则上述非线性系统是ISS的,欧几里德范数由|·|表示,γISS称为输入到状态渐近增益;||u||∞表示外部输入u的上确界,它是使得|u(t)|≤a对所有t成立的最小的a;比较函数的定义如下:
定义2.对于任何x0∈ΩIOS∈Rn,u∈UIOS∈Rm;如果存在函数
x0∈ΩISS∈Rn,u∈UISS∈Rm,下式成立:
则上述非线性系统称为局部输入输出稳定系统;
定理1.在非线性系统的局部输入-状态稳定LISS/局部输入-输出稳定LIOS分析中,ΩISS,UISS,和ΩIOS,UIOS分别表示上述非线性系统的LISS域和LIOS域;在实际系统中,上述非线性系统的稳定域用Ω和U来表示,其中Ω=ΩISS∩ΩIOS,U=UISS∩UIOS;
定义3.对于任意x0∈ΩISS∈Rn,u∈UISS∈Rm,如果下列不等式成立
同样,局部输入-输出积分-积分估计的定义为:
由上述积分-积分估计的定义看出,系统状态、外部输入和系统输出均由积分表示,包含外部扰动的幅值信息以及随时间波动的信息;
基于定义1-3及定理1,采用电力系统LISS/LIOS分析的实用算法,利用系统初始状态和外部输入的局部定义域和渐近增益对非线性系统的动态响应进行量化评估;
利用系统初始状态和外部输入的局部定义域和渐近增益对非线性系统的动态响应进行量化评估的具体过程为:
设系统初始状态位于平衡点,当考虑线性增益函数时,渐近增益表示如下:
其中,T为步长;
通过改变输入信号u的幅值,得到一系列的渐近增益,选择整个输入域上的最大γ作为系统线性形式的渐近增益将导致估计结果非常保守,为了降低LISS/LIOS分析的保守性,采用分段线性渐近增益;
考虑一个由n个子系统组成的互联系统,其中第i个子系统数学模型表达式如下:
基于积分-积分估计的互联系统量化稳定性判据如下:
如果满足下列四个条件,则上述n个子系统组成的互联系统是LISS和LIOS的,说明对双馈风场接入电网友好:
II.每个子系统的输入/输出满足以下条件:
|u(s)|≤Z|y(s)|
其中Zij≥0,当选择电压和电流作为输入输出时,由系统的导纳矩阵导出Z;
III.以下小增益条件成立:
G(υ*)<<υ*
IV.系统输出满足预先给定如下条件:
2.根据权利要求1所述的一种对双馈风场接入电网友好性的量化评估方法,其特征在于,
(1)子系统层面:基于积分-积分估计,根据不同控制策略对双馈风场输出特性性的量化评估;
(2)互联系统层面:基于小增益定理和积分-积分估计,对双馈风场接入电网友好性的量化评估。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811526140.7A CN109474014B (zh) | 2018-12-13 | 2018-12-13 | 一种对双馈风场接入电网友好性的量化评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811526140.7A CN109474014B (zh) | 2018-12-13 | 2018-12-13 | 一种对双馈风场接入电网友好性的量化评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109474014A CN109474014A (zh) | 2019-03-15 |
CN109474014B true CN109474014B (zh) | 2021-03-23 |
Family
ID=65675316
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811526140.7A Active CN109474014B (zh) | 2018-12-13 | 2018-12-13 | 一种对双馈风场接入电网友好性的量化评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109474014B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111555311B (zh) * | 2020-05-20 | 2021-10-15 | 国网陕西省电力公司电力科学研究院 | 一种电力系统即插即用的稳定性分析与控制方法 |
CN114362167A (zh) * | 2022-01-12 | 2022-04-15 | 西安交通大学 | 一种电力系统暂态电压稳定性的评估方法 |
CN115360719B (zh) * | 2022-08-30 | 2024-04-12 | 东北大学秦皇岛分校 | 基于plnn的电力系统短期电压稳定性评估方法 |
CN116577995B (zh) * | 2023-07-05 | 2023-11-21 | 广州惠特节能科技有限公司 | 热泵烘干机组参数智能控制方法及系统 |
CN118013862A (zh) * | 2024-04-09 | 2024-05-10 | 湖南大学 | 一种风电场短期电压稳定性仿真模型构建、评估及监测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107123997A (zh) * | 2017-04-18 | 2017-09-01 | 燕山大学 | 一种大型多机电力系统的分散协调分层鲁棒控制方法 |
CN107579519A (zh) * | 2017-09-25 | 2018-01-12 | 华中科技大学 | 一种考虑多状态切换的交直流混联电网鲁棒最优控制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102611381A (zh) * | 2012-03-12 | 2012-07-25 | 浙江工业大学 | 永磁同步电机直接转矩控制系统 |
CN108964127B (zh) * | 2018-07-27 | 2020-03-13 | 西安交通大学 | 一种双馈风力发电系统故障穿越的控制方法 |
-
2018
- 2018-12-13 CN CN201811526140.7A patent/CN109474014B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107123997A (zh) * | 2017-04-18 | 2017-09-01 | 燕山大学 | 一种大型多机电力系统的分散协调分层鲁棒控制方法 |
CN107579519A (zh) * | 2017-09-25 | 2018-01-12 | 华中科技大学 | 一种考虑多状态切换的交直流混联电网鲁棒最优控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109474014A (zh) | 2019-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109474014B (zh) | 一种对双馈风场接入电网友好性的量化评估方法 | |
Yang et al. | Perturbation observer based fractional-order sliding-mode controller for MPPT of grid-connected PV inverters: Design and real-time implementation | |
Divya et al. | Models for wind turbine generating systems and their application in load flow studies | |
Hossain et al. | Control strategies for augmenting LVRT capability of DFIGs in interconnected power systems | |
Liu et al. | Coordinated fault‐ride‐through strategy for doubly‐fed induction generators with enhanced reactive and active power support | |
Koessler et al. | Integration of large wind farms into utility grids pt. I-Modeling of DFIG | |
CN108964127B (zh) | 一种双馈风力发电系统故障穿越的控制方法 | |
Jena et al. | Setting a fostered energy network by decarbonizing the grid: H ybridization, control, and future solutions upon storage | |
Okedu et al. | Comparative study of the effects of machine parameters on DFIG and PMSG variable speed wind turbines during grid fault | |
Yang et al. | Control of superconducting magnetic energy storage systems in grid‐connected microgrids via memetic salp swarm algorithm: An optimal passive fractional‐order PID approach | |
Sun et al. | MPPT adaptive controller of DC-based DFIG in resistances uncertainty | |
Zheng et al. | Transient stable region of power systems incorporating stochasticity of variable renewable energies and system contingencies | |
CN113852120A (zh) | 一种新能源发电最大并网容量的确定方法及系统 | |
CN114465280A (zh) | 一种新能源并网系统动态等效建模方法 | |
Liu et al. | Temporary overvoltage assessment and suppression in heterogeneous renewable energy power systems | |
Acevedo et al. | Dynamic voltage stability comparison of thermal and wind power generation with different static and dynamic load models | |
Olanite et al. | Grid Integration of Wind Power System for Voltage Stability Enhancement | |
Liao et al. | A study on transient overvoltage caused by HVDC tie blocking in multi-rectifier systems | |
Belabbas et al. | Hierarchical energy management and control to improve the reliability and efficiency of wind farms connected to the grid | |
Datta et al. | Three-phase steady state model for unbalanced operation of grid-connected wind generation unit | |
Tan et al. | A cut-out strategy for wind turbines that ensures low-voltage ride-through capability | |
Li et al. | Wind farms reactive power optimization using genetic/tabu hybrid algorithm | |
Jervekani et al. | LVRT capability based on P‐V curve fitting under partial shading conditions | |
Ma et al. | Cyber-physical Modeling Technique based Dynamic Aggregation of Wind Farm Considering LVRT Characteristic | |
CN111555354B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |