CN104915892A - 电力系统随机时滞稳定性分析方法及装置 - Google Patents
电力系统随机时滞稳定性分析方法及装置 Download PDFInfo
- Publication number
- CN104915892A CN104915892A CN201510325561.3A CN201510325561A CN104915892A CN 104915892 A CN104915892 A CN 104915892A CN 201510325561 A CN201510325561 A CN 201510325561A CN 104915892 A CN104915892 A CN 104915892A
- Authority
- CN
- China
- Prior art keywords
- stability
- upper limit
- omega
- time
- module
- 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
Links
Landscapes
- Control Of Eletrric Generators (AREA)
Abstract
一种电力系统随机时滞稳定性分析方法及装置,所述装置包括数据采集模块、状态方程构建模块、时滞上限求解模块和结果输出模块,其中数据采集模块用于采集网络结构参数、系统潮流值、发电机功角与转速和风速;状态方程构建模块用于构建互联系统状态方程;时滞上限求解模块用于构造互联系统在随机激励下的稳定性判据,并确定系统时滞稳定上限的变化情况;结果输出模块用于输出系统时滞稳定上限结果。通过利用本发明的电力系统随机时滞稳定性分析方法及装置,能够解决现有技术中未考虑随机激励对系统时滞稳定性建模和分析带来的问题,在不同工况下保证电力系统稳定性分析的正确性与有效性。
Description
技术领域
本发明涉及电力系统稳定性分析技术领域,尤其是涉及到电力系统随机时滞稳定性分析技术。
背景技术
广域量测系统(Wide-Area Measurement System,WAMS)为电力系统的全局稳定分析、运行优化和协调控制带来了新的契机,但WAMS存在明显的时滞,对系统控制策略的制定以及控制器的效果将产生不利影响,已成为系统不稳定和控制器性能变差的重要诱因。因此,分析电力系统的时滞稳定性和确定电力系统所能承受的最大时滞对保证其稳定运行意义重大。
目前,对电力系统时滞稳定性的分析方法主要有频域法和时域法:频域法通过计算特征根的分布情况判别系统的稳定性,是分析时滞系统稳定性较早的方法;时域法主要基于李雅普诺夫-克拉索夫斯基稳定性理论以及拉实米辛稳定性理论,是当前分析时滞系统稳定性主流的方法。然而,上述方法均采用常微分方程组描述系统状态,属于确定性分析方法。近年来,随着波动性新能源装机容量的增大、接入电压等级的升高,大量随机激励使常规电力系统时滞稳定性建模和分析方法面临严峻挑战。
发明内容
有鉴于此,本发明的目的旨在解决现有技术中未考虑随机激励对系统时滞稳定性建模和分析带来的问题,提出一种考虑风电并网的电力系统随机时滞稳定性分析方法及装置。以双馈异步电机为例,首先,建立考虑随机风功率的双馈异步电机维纳模型,在此基础上,构建基于随机微分方程(Stochastic DifferentialEquation,SDE)描述的互联系统状态方程,然后,构造考虑风功率随机性的李雅普诺夫-克拉索夫斯基目标泛函,并借助伊藤微分公式求解该目标泛函的弱无穷小生成算子,并据此制定互联系统在随机激励下稳定性判据,最后,利用广义特征值(Generalized Eigenvalue Problem,GEVP)计算系统时滞稳定上限随风功率随机激励强度、发电机励磁放大系数以及阻尼系数的变化情况。系统仿真结果表明本发明在不同运行工况下的正确性及有效性。
为了实现此目的,本发明采取的技术方案为如下。
一种电力系统随机时滞稳定性分析装置,所述装置包括顺序相连的数据采集模块、状态方程构建模块、时滞上限求解模块和结果输出模块,其中,
数据采集模块用于采集网络结构参数、系统潮流值、发电机功角与转速和风速,并将采集数据发送至状态方程构建模块;
状态方程构建模块用于构建互联系统状态方程;
时滞上限求解模块用于构造互联系统在随机激励下的稳定性判据,并确定系统时滞稳定上限的变化情况;
结果输出模块用于输出系统时滞稳定上限结果。
其中所述状态方程构建模块根据数据采集模块采集到的数据,建立考虑随机风功率的双馈异步电机维纳模型,在此基础上,构建基于随机微分方程描述的互联系统状态方程。
所述时滞上限求解模块包括稳定性判据构造单元和系统时滞稳定上限单元,其中,
稳定性判据构造单元用于构造考虑风功率随机性的李雅普诺夫-克拉索夫斯基目标泛函,并借助伊藤微分公式确定该目标泛函的弱无穷小生成算子,制定互联系统在随机激励下稳定性判据;
系统时滞稳定上限单元用于利用广义特征值确定系统时滞稳定上限随随机激励强度、发电机励磁放大系数以及阻尼系数的变化情况。
一种电力系统随机时滞稳定性分析方法,所述方法包括步骤:
A、采集网络结构参数、系统潮流值、发电机功角与转速、风速;
B、构建互联系统状态方程;
C、构造互联系统在随机激励下的稳定性判据,并确定系统时滞稳定上限的变化情况;
D、输出系统时滞稳定上限结果。
步骤B中根据采集到的数据,建立考虑随机风功率的双馈异步电机维纳模型,在此基础上,构建基于随机微分方程描述的互联系统状态方程。
另外,步骤B中建立考虑随机风功率的双馈异步电机维纳模型为:
其中,
x3为跟踪定子d轴电压的误差积累,
θ为同步旋转的x-y坐标系和d-q坐标系的夹角,
s为转差,
为双馈异步电机平衡点,
ω(t)为维纳过程,
t为时间,
H为惯性系数,
isd0、isq0分别为定子交、直轴稳态电流值,
Ls为定子自感,Lr为转子自感,Lm为定子与转子间的互感,
U0为风机定子电压,
Kp3、Ki3分别为风机锁相环比例、积分系数,
为随机激励强度系数。
且步骤B中构建基于随机微分方程描述的互联系统状态方程为:
dx=Axdt+Adx(t-ht)dt+Bωdωt,
其中,
x=[x3,θ,s,δ,ω,E′q,Ef]T,
δ、ω、E′q、Ef分别为发电机功角、角速度、电抗后电势和励磁电势,ωs为系统额定转速,ωt为维纳过程,D为发电机阻尼系数,T′d、TA分别为发电机定子和励磁回路时间常数,Pe为发电机输出功率,Vt为发电机机端电压,Eq为发电机空载电动势,ht为系统时滞,M为同步发电机组惯性时间常数。
另外,步骤C进一步包括:
C1、构造考虑风功率随机性的李雅普诺夫-克拉索夫斯基目标泛函,并借助伊藤微分公式求解该目标泛函的弱无穷小生成算子,制定互联系统在随机激励下稳定性判据;
C2、利用广义特征值确定系统时滞稳定上限随随机激励强度、发电机励磁放大系数以及阻尼系数的变化情况。
所述互联系统在随机激励下稳定性判据为:
对于给定标量和μ,若存在P=PT>0,Q=QT≥0,R=RT≥0,K=KT>0,Z=ZT≥0和标量ε>0,以及适维矩阵L、M和S,使得如下线性矩阵不等式成立,则所述互联系统在随机激励下稳定,
Ω2=[L+SA SAd -L 0 -S],
通过采用本发明的电力系统随机时滞稳定性分析方法和装置,能解决现有技术中未考虑随机激励对系统时滞稳定性建模和分析带来的问题,在不同工况下能够保证电力系统稳定性分析的正确性与有效性。
附图说明
图1是本发明具体实施方式中电力系统随机时滞稳定性分析装置的结构示意图。
图2是为说明本发明具体实施方式中电力系统随机时滞稳定性分析装置及方法的单机无穷大系统的结构示意图。
图3是本发明具体实施方式中处于不同时滞时同步发电机功角的变化曲线示意图。
图4是本发明具体实施方式中系统时滞稳定上限随随机激励强度系数的变化情况示意图。
图5是本发明具体实施方式中处于不同随机激励强度系数时同步发电机功角的变化曲线示意图。
图6是本发明具体实施方式中系统时滞稳定上限发电机随励磁放大倍数KA的变化情况示意图。
图7是本发明具体实施方式中处于不同发电机励磁放大倍数KA时发电机功角的变化曲线示意图。
图8是本发明具体实施方式中系统时滞稳定上限随发电机阻尼系数D的变化情况示意图。
图9是本发明具体实施方式中处于不同发电机阻尼系数D时发电机功角的变化曲线
具体实施方式
下面结合附图,对本发明作详细说明。
以下公开详细的示范实施例。然而,此处公开的具体结构和功能细节仅仅是出于描述示范实施例的目的。
然而,应该理解,本发明不局限于公开的具体示范实施例,而是覆盖落入本公开范围内的所有修改、等同物和替换物。在对全部附图的描述中,相同的附图标记表示相同的元件。
同时应该理解,如在此所用的术语“和/或”包括一个或多个相关的列出项的任意和所有组合。另外应该理解,当部件或单元被称为“连接”或“耦接”到另一部件或单元时,它可以直接连接或耦接到其他部件或单元,或者也可以存在中间部件或单元。此外,用来描述部件或单元之间关系的其他词语应该按照相同的方式理解(例如,“之间”对“直接之间”、“相邻”对“直接相邻”等)。
介绍本发明的具体实施方式之前,首先介绍本发明电力系统随机时滞稳定性分析方法及装置的原理,并结合本发明的技术方案进行了分析和演算。
双馈风机的电磁模型如下:
1.定子电压方程:
2.转子电流方程:
3.转子电流控制方程:
4.锁相环跟踪方程:
5.d-q坐标系运动方程:
式中,usd和usq分别为双馈感应电机定子的d轴和q轴电压,ψsd和ψsq分别为定子的d轴和q轴磁链,isd和isq分别为定子的d轴和q轴电流,ωs为坐标系的旋转速度,rs为电枢的电阻,ird和irq分别为双馈异步电机转子的d轴和q轴电流,urd和urq分别为转子的d轴和q轴电压,ωr为转子转速,rr为转子电阻,Ls为定子自感,Lr为转子自感,Lm为定子与转子间的互感,和分别为转子d轴和q轴电流的给定参考值,x1和x2分别为转子d轴和q轴电流跟踪误差的积累,x3为跟踪定子d轴电压的误差积累,θ为同步旋转的x-y坐标系和d-q坐标系的夹角,ωn为电网频率,式(3)中,
电磁部分8阶详细模型的代数部分如下式(9)-式(15):
1.定子电流方程:
2.转子磁链方程:
ψrd=Lmisd+Lrird (11)
ψrq=Lmisq+Lrirq (12)
3.转子电压方程:
4.d-q坐标系旋转速度方程:
ωs=-Kp3usd+Ki3x3 (15)
其中,ψrd和ψrq分别为转子d轴和q轴的磁链,Kp1,Kp2,Kx1,Kx2为转子电流控制器的参数,Kp3和Ki3为控制器参数。
为获得一般性结论,需对双馈风力发电机8阶电磁模型进行降阶,从而获得进行稳定分析的简化模型。由于双馈感应电机中,定子磁链暂态衰减过程比转子电流暂态衰减过程慢得多,因此,双馈感应电机的定子磁链暂态过程对系统的机电暂态影响不大。同时,数值计算表明,仅保留锁相环2阶动态的简化模型与8阶电磁模型具有相近的稳定边界。另一方面,风电机组轴系将引入相关的振荡模态,同时电力系统网络故障时,可能激起风力机械系统某些频率的振荡和摇摆,因此,本发明在锁相环动态2阶简化模型的基础上,引入风机的机械模型,避免了因忽略机械部分带来的误差。风电机组机电动力学方程为:
其中,H为惯性系数;Pe为电磁功率,s为转差,Pm为风机输入机械功率。Pe和Pm分别按式(17)和式(18)计算:
Pm=P0+PΔ (17)
式中,P0为风机机械功率的确定部分;PΔ为由于风速随机变化引起的机械功率波动,可以采用维纳过程描述为:
式中,ω(t)为维纳过程;为随机激励强度系数,表示功率随机波动的强度。式(11)、式(12)和式(16)组成双馈异步风机的3阶动态模型,设双馈风机的一个平衡点为在平衡点处将3阶动态模型线性化,状态方程可表示为:
式中:
以双馈异步风机和同步发电机并联接入无穷大系统为例,分析电力系统随机时滞稳定性,本发明实施方式中双机无穷大系统结构图如图2所示。
考虑励磁系统输出电压存在时滞,且系统模型为4阶微分方程,则时滞系统状态方程可表示为:
式中:δ,ω,E′q,Ef分别为发电机功角、角速度、电抗后电势和励磁电势;Pm为原动机输出功率;ωs为系统额定转速;D为发电机阻尼系数;T′d,TA分别为发电机定子和励磁回路时间常数;KA为励磁回路放大系数;Vref为机端电压的参考值;xd,x′d分别为发电机稳态和暂态电抗;xe为线路电抗;V0为无穷大母线电压;Pe为发电机输出功率;Vt,id分别为发电机机端电压和纵轴输出电流;M为同步发电机组惯性时间常数。
在平衡点(δ0,ω0,E′q0,Ef0)处,将4阶动态模型线性化,状态方程如下所示:
其中,
由式(20)与式(23)组成的2机系统的随机状态方程为:
dx=Axdt+Adx(t-ht)dt+Bωdωt
=f(t)dt+Bωdωt (25)
其中:
x=[x3,θ,s,δ,ω,E′q,Ef]T,
由式(25)可以看出,时滞系统含有随机因素,传统分析时滞稳定性的方法已不适用,本发明借助伊藤微分理论分析随机系统(25)的时滞稳定性。
前提1:对于式(25)所示的随机时滞电力系统,给定标量和μ,若存在P=PT>0,Q=QT≥0,R=RT≥0,K=KT>0,Z=ZT≥0和标量ε>0,以及适维矩阵L、M和S,,使得如下线性矩阵不等式(26)成立,则随机时滞系统(25)均方稳定。
其中:
Ω2=[L+SA SAd -L 0 -S],
说明:定义向量城y(t)∈Rn,使其满足:
构造如下形式的李雅普诺夫-克拉索夫斯基泛函:
由微分公式,V(xt,t)沿着系统(25)的随机微分为:
dV(t,x(t))=ζVdt+2xT(t)PBωdω(t) (29)
其中,弱无穷小生成算子为:
由牛顿-莱布尼兹公式可知,对于任意的适维矩阵L和M,有下列式子成立:
其中,
同时,由f(t)=Ax(t)+Adx(t-ht)可知,对于任意适维矩阵S,以下等式恒成立:
2ζT(t)S[Ax(t)-Adx(t-ht)-f(t)]=0 (34)
进一步,利用詹森不等式,以下不等式均成立:
综合以上分析,将式(21)-(26)代入式(20),整理可得:
由于Z>0,若下式成立:
则对于任意的ξ(t)≠0,均有E{ζV(t,x(t))}<0成立,随机系统(15)在概率空间上渐近均方稳定。基于舒尔补引理,式(39)与式(26)等价,定理1得证。
式(26)表征的矩阵不等式仅能判定系统是否稳定,而无法获取系统时滞稳定上限等信息,考虑到时滞稳定上限的求解是一个具有线性不等式约束的凸优化问题,具有广义特征值(GEVP)的形式。因此,本发明提出利用GVEP方法计算系统的时滞稳定上限。由于式(26)不是标准的GVEP形式,需要进行必要的处理。
根据舒尔补引理,式(26)可变形为:
令 以及 得:
在式(41)中,将-Z1,-R1分别代替得到:
由式(40)-式(42)可知,系统时滞稳定上限问题即可转化为如下优化问题:
min τ
s.t.(41),(42) (43)
通过求解式(43)中以式(41)和式(42)为约束的最小τ,然后,利用可计算出系统时滞稳定上限。
因此,如图1所示,本发明具体实施方式中公开了一种电力系统随机时滞稳定性分析装置,所述装置包括顺序相连的数据采集模块、状态方程构建模块、时滞上限求解模块和结果输出模块,其中,
数据采集模块用于采集网络结构参数、系统潮流值、发电机功角与转速和风速,并将采集数据发送至状态方程构建模块;
状态方程构建模块用于构建互联系统状态方程;
时滞上限求解模块用于构造互联系统在随机激励下的稳定性判据,并确定系统时滞稳定上限的变化情况;
结果输出模块用于输出系统时滞稳定上限结果。
因此,从电力系统随机时滞稳定性分析装置的配置可以看出,通过采用本发明的技术方案,能解决现有技术中未考虑随机激励对系统时滞稳定性建模和分析带来的问题。
本发明具体实施方式的电力系统随机时滞稳定性分析装置中,所述状态方程构建模块根据数据采集模块采集到的数据,建立考虑随机风功率的双馈异步电机维纳模型,在此基础上,构建基于随机微分方程描述的互联系统状态方程。
另外,所述时滞上限求解模块包括稳定性判据构造单元和系统时滞稳定上限单元,其中稳定性判据构造单元用于构造考虑风功率随机性的李雅普诺夫-克拉索夫斯基目标泛函,并借助伊藤微分公式确定该目标泛函的弱无穷小生成算子,制定互联系统在随机激励下稳定性判据;系统时滞稳定上限单元用于利用广义特征值确定系统时滞稳定上限随随机激励强度、发电机励磁放大系数以及阻尼系数的变化情况。
与本发明的电力系统随机时滞稳定性分析装置相适应,本发明还公开了一种电力系统随机时滞稳定性分析方法,所述方法包括步骤:
A、采集网络结构参数、系统潮流值、发电机功角与转速、风速;
B、构建互联系统状态方程;
C、构造互联系统在随机激励下的稳定性判据,并确定系统时滞稳定上限的变化情况;
D、输出系统时滞稳定上限结果。
具体实施方式的步骤B中根据采集到的数据,建立考虑随机风功率的双馈异步电机维纳模型,在此基础上,构建基于随机微分方程描述的互联系统状态方程。
特别地,步骤B中建立考虑随机风功率的双馈异步电机维纳模型为:
其中,
步骤B中构建基于随机微分方程描述的互联系统状态方程为:
dx=Axdt+Adx(t-ht)dt+Bωdωt,
其中,
x=[x3,θ,s,δ,ω,E′q,Ef]T,
另外,步骤C进一步包括:
C1、构造考虑风功率随机性的李雅普诺夫-克拉索夫斯基目标泛函,并借助伊藤微分公式求解该目标泛函的弱无穷小生成算子,制定互联系统在随机激励下稳定性判据;
C2、利用广义特征值确定系统时滞稳定上限随随机激励强度、发电机励磁放大系数以及阻尼系数的变化情况。
所述互联系统在随机激励下稳定性判据为:
对于给定标量和μ,若存在P=PT>0,Q=QT≥0,R=RT≥0,K=KT>0,Z=ZT≥0和标量ε>0,以及适维矩阵L、M和S,使得如下线性矩阵不等式成立,则所述互联系统在随机激励下稳定,
Ω2=[L+SA SAd -L 0 -S],
以下,通过具体实施例来说明本发明的技术效果,所述具体实施例的应用场景如图2所示,为单机无穷大系统。
其中同步机系统参数如下:Pm=1.3,ωs=377,D=10,T′d=10,TA=1,KA=190,Vref=1.05,xd=1,x′d=0.4,xe=0.5,Vt=1.0。
双馈异步风机系统参数如下:rr=0.006298,H=3.025,Ls=0.004872,Lr=0.005098,Lm=0.004603,L=0.00075,Kp1=0.155,Kp2=0.155,Kx2=15,Kp3=1,Ki3=345。利用上述参数可以确定:当随机激励强度时,随机时滞系统(25)在平衡点处的线性化状态方程系数矩阵为:
为验证本方法的正确性及有效性,将式(44)-式(48)代入线性矩阵不等式工具箱中的GEVP求解器,计算随机时滞系统在时滞变化率μ=0.1时的系统时滞稳定上限为即当0≤h≤328.1ms时,随机系统时滞稳定;当h>328.1ms时,系统失去稳定,将随机系统时滞大小分别设置为50ms、328.1ms以及360ms时,系统发电机功角的变化曲线如图3所示。由图3可以看出,当时滞h满足0≤h≤328.1ms时,随机系统时滞稳定;当h>328.1ms时,系统失去稳定,这与利用本方法得到的结论一致。
1随机激励强度变化对系统时滞稳定上限的影响
首先,令随机激励强度在0.02~0.2之间变动,且固定系统的其他参数,基于GEVP计算系统相应的系统时滞稳定上限,结果如图4所示,可以看出,系统时滞稳定上限随的变化规律较为复杂:在增大的过程中,系统的时滞稳定上限先减小后增大,因此,系统随机时滞稳定域是一个非凸集,存在一个明显的凹区域。需要注意的是,在空间中,曲线对应系统的随机时滞稳定域的上边界,其下方为时滞系统的随机时滞稳定域。
为验证图4的正确性,分别取以及此时系统(25)的时滞设置为h=220ms,在上述典型运行方式下,对系统进行时域仿真,发电机功角的变化曲线如图5所示。由图5可知,当时滞h=250ms时,系统处于或的运行状态下稳定,而系统在的状态下不稳定,这表明当或时,系统的时滞稳定上限均大于250ms,而当时,时滞稳定上限小于250ms,这与图4中时滞稳定上限随的变化情况吻合,在图4中,当以及时,时滞稳定上限分别为259.81ms、202.61ms、271.64ms。同时,对比和的发电机功角曲线,可以发现,虽然系统均处于稳定状态,但是当时发电机功角达到稳定值的时间较长,这说明随机激励较大时会削弱系统的稳定性。
2发电机励磁放大系数KA变化对系统时滞稳定上限的影响
进一步,令发电机励磁放大系数KA在50~210之间变动,固定其他参数不变,计算结果见图6。由图6可以看出,系统的时滞稳定上限随发电机励磁放大倍数KA的增大而减小,这种变化趋势可以解释为:发电机励磁放大系数KA越大表明发电机励磁系统灵敏度越高,励磁电压对机端电压与电压参考值的差值越灵敏,由时滞产生的任何误差对系统稳定性影响越大,因此,发电机励磁放大倍数越大,系统的时滞稳定上限越小。
为验证图6的正确性,分别取发电机励磁放大倍数KA=180、KA=190以及KA=200,此时系统(25)的时滞设置为h=320ms,在上述三种运行方式下对系统进行了时域仿真,同步发电机功角曲线如图7所示。由图7可知,当时滞h=320ms时,发电机励磁放大倍数KA=180或KA=190的系统处于稳定运行状态,而在KA=200的状态下系统不稳定,这说明当KA=180或KA=190时,系统时滞稳定上限均大于320ms,而当KA=200时,时滞稳定上限小于320ms,符合图6中时滞稳定上限随KA的变化趋势,验证了本方法的正确性。
3发电机阻尼系数D变化对系统时滞稳定上限的影响
进一步,令发电机阻尼系数D在0~40之间变动,固定其他参数不变,计算结果如图8所示。理论分析可知,电力系统的稳定性随着阻尼的增加而增强,因此系统保持稳定性所容许的时滞也随之增大。由图8可以看出,当阻尼系数由0增大到60时,系统时滞稳定上限随D的增大而增大,这种变化趋势与理论分析结果一致。
为验证图8的合理性,分别取发电机阻尼系数D=0、D=10以及D=20,此时系统(25)的时滞设置为h=320ms,在上述状态下对系统进行了时域仿真,同步发电机功角曲线如图9所示,可知在时滞h=320ms时,当D=10或D=20时,系统均处于稳定运行状态,而当D=0时,系统处于不稳定状态,由此表明当D=10或D=20时,系统时滞稳定上限均大于320ms,而当D=0时,时滞稳定上限小于320ms,这符合图8中时滞稳定上限的变化情况。
需要说明的是,上述实施方式仅为本发明较佳的实施方案,不能将其理解为对本发明距离保护范围的限制,在未脱离本发明构思前提下,对本发明所做的任何微小变化与修饰均属于本发明的距离保护范围。
Claims (9)
1.一种电力系统随机时滞稳定性分析装置,所述装置包括顺序相连的数据采集模块、状态方程构建模块、时滞上限求解模块和结果输出模块,其中,
数据采集模块用于采集网络结构参数、系统潮流值、发电机功角与转速和风速,并将采集数据发送至状态方程构建模块;
状态方程构建模块用于构建互联系统状态方程;
时滞上限求解模块用于构造互联系统在随机激励下的稳定性判据,并确定系统时滞稳定上限的变化情况;
结果输出模块用于输出系统时滞稳定上限结果。
2.根据权利要求1中所述的电力系统随机时滞稳定性分析装置,其特征在于,所述状态方程构建模块根据数据采集模块采集到的数据,建立考虑随机风功率的双馈异步电机维纳模型,在此基础上,构建基于随机微分方程描述的互联系统状态方程。
3.根据权利要求1中所述的电力系统随机时滞稳定性分析装置,其特征在于,所述时滞上限求解模块包括稳定性判据构造单元和系统时滞稳定上限单元,其中,
稳定性判据构造单元用于构造考虑风功率随机性的李雅普诺夫-克拉索夫斯基目标泛函,并借助伊藤微分公式确定该目标泛函的弱无穷小生成算子,制定互联系统在随机激励下稳定性判据;
系统时滞稳定上限单元用于利用广义特征值确定系统时滞稳定上限随随机激励强度、发电机励磁放大系数以及阻尼系数的变化情况。
4.一种电力系统随机时滞稳定性分析方法,所述方法包括步骤:
A、采集网络结构参数、系统潮流值、发电机功角与转速、风速;
B、构建互联系统状态方程;
C、构造互联系统在随机激励下的稳定性判据,并确定系统时滞稳定上限的变化情况;
D、输出系统时滞稳定上限结果。
5.根据权利要求4中所述的电力系统随机时滞稳定性分析方法,其特征在于,步骤B中根据采集到的数据,建立考虑随机风功率的双馈异步电机维纳模型,在此基础上,构建基于随机微分方程描述的互联系统状态方程。
6.根据权利要求5中所述的电力系统随机时滞稳定性分析方法,其特征在于,步骤B中建立考虑随机风功率的双馈异步电机维纳模型为:
其中,
x3为跟踪定子d轴电压的误差积累,
θ为同步旋转的x-y坐标系和d-q坐标系的夹角,
s为转差,
为双馈异步电机平衡点,
ω(t)为维纳过程,
t为时间,
H为惯性系数,
isd0、isq0分别为定子交、直轴稳态电流值,
Ls为定子自感,Lr为转子自感,Lm为定子与转子间的互感,
U0为风机定子电压,
Kp3、Ki3分别为风机锁相环比例、积分系数,
为随机激励强度系数。
7.根据权利要求5中所述的电力系统随机时滞稳定性分析方法,其特征在于,步骤B中构建基于随机微分方程描述的互联系统状态方程为:
dx=Axdt+Adx(t-ht)dt+Bωdωt,
其中,
x=[x3,θ,s,δ,ω,E′q,Ef]T,
δ、ω、E′q、Ef分别为发电机功角、角速度、电抗后电势和励磁电势,ωs为系统额定转速,ωt为维纳过程,D为发电机阻尼系数,T′d、TA分别为发电机定子和励磁回路时间常数,Pe为发电机输出功率,Vt为发电机机端电压,Eq为发电机空载电动势,ht为系统时滞,M为同步发电机组惯性时间常数。
8.根据权利要求7中所述的电力系统随机时滞稳定性分析方法,其特征在于,步骤C进一步包括:
C1、构造考虑风功率随机性的李雅普诺夫-克拉索夫斯基目标泛函,并借助伊藤微分公式求解该目标泛函的弱无穷小生成算子,制定互联系统在随机激励下稳定性判据;
C2、利用广义特征值确定系统时滞稳定上限随随机激励强度、发电机励磁放大系数以及阻尼系数的变化情况。
9.根据权利要求8中所述的电力系统随机时滞稳定性分析方法,其特征在于,所述互联系统在随机激励下稳定性判据为:
对于给定标量和μ,若存在P=PT>0,Q=QT≥0,R=RT≥0,K=KT>0,Z=ZT≥0和标量ε>0,以及适维矩阵L、M和S,使得如下线性矩阵不等式成立,则所述互联系统在随机激励下稳定,
Ω2=[L+SA SAd -L 0 -S],
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510325561.3A CN104915892B (zh) | 2015-06-12 | 2015-06-12 | 电力系统随机时滞稳定性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510325561.3A CN104915892B (zh) | 2015-06-12 | 2015-06-12 | 电力系统随机时滞稳定性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104915892A true CN104915892A (zh) | 2015-09-16 |
CN104915892B CN104915892B (zh) | 2018-12-18 |
Family
ID=54084935
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510325561.3A Active CN104915892B (zh) | 2015-06-12 | 2015-06-12 | 电力系统随机时滞稳定性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104915892B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105335904A (zh) * | 2015-11-30 | 2016-02-17 | 广东工业大学 | 一种电力系统Lyapunov稳定性分析方法 |
CN107069774A (zh) * | 2017-04-01 | 2017-08-18 | 华北电力大学 | 电力系统能量函数优选及暂态稳定裕度评估系统及方法 |
CN107069712A (zh) * | 2017-04-01 | 2017-08-18 | 华北电力大学 | 基于李雅普诺夫的电力系统临界切除时间计算系统及方法 |
CN107832264A (zh) * | 2017-11-20 | 2018-03-23 | 宁德师范学院 | 一种电机设计的综合评价系统及评价方法 |
CN112383057A (zh) * | 2020-11-27 | 2021-02-19 | 哈尔滨工程大学 | 基于电网潮流的电力耦合系统互感和自感设计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040167731A1 (en) * | 2002-12-23 | 2004-08-26 | Abb, Inc. | Failure rate adjustment for electric power network reliability analysis |
CN104680323A (zh) * | 2015-03-02 | 2015-06-03 | 华北电力大学 | 电力系统能量稳定域构建系统及方法 |
-
2015
- 2015-06-12 CN CN201510325561.3A patent/CN104915892B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040167731A1 (en) * | 2002-12-23 | 2004-08-26 | Abb, Inc. | Failure rate adjustment for electric power network reliability analysis |
CN104680323A (zh) * | 2015-03-02 | 2015-06-03 | 华北电力大学 | 电力系统能量稳定域构建系统及方法 |
Non-Patent Citations (2)
Title |
---|
和萍等: "风电场并网对互联系统小干扰稳定及低频振荡特性的影响", 《电力系统自动化》 * |
马静等: "基于改进自由权矩阵与广义特征值的时滞稳定上限计算方法研究", 《电力系统保护与控制》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105335904A (zh) * | 2015-11-30 | 2016-02-17 | 广东工业大学 | 一种电力系统Lyapunov稳定性分析方法 |
CN107069774A (zh) * | 2017-04-01 | 2017-08-18 | 华北电力大学 | 电力系统能量函数优选及暂态稳定裕度评估系统及方法 |
CN107069712A (zh) * | 2017-04-01 | 2017-08-18 | 华北电力大学 | 基于李雅普诺夫的电力系统临界切除时间计算系统及方法 |
CN107832264A (zh) * | 2017-11-20 | 2018-03-23 | 宁德师范学院 | 一种电机设计的综合评价系统及评价方法 |
CN112383057A (zh) * | 2020-11-27 | 2021-02-19 | 哈尔滨工程大学 | 基于电网潮流的电力耦合系统互感和自感设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104915892B (zh) | 2018-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104915892A (zh) | 电力系统随机时滞稳定性分析方法及装置 | |
Ullah et al. | Adaptive fractional order terminal sliding mode control of a doubly fed induction generator-based wind energy system | |
CN103259479B (zh) | 一种永磁同步电机神经网络左逆状态观测方法 | |
Shabani et al. | Real-time transient stability detection in the power system with high penetration of DFIG-based wind farms using transient energy function | |
CN109753689A (zh) | 一种电力系统机电振荡模态特征参数在线辩识方法 | |
CN104680323B (zh) | 电力系统能量稳定域构建系统及方法 | |
CN103915839A (zh) | 一种用于分析含风电电力系统随机稳定性的方法 | |
CN108988320A (zh) | 电力系统动态元件响应特性对暂态电压稳定性影响分析方法 | |
CN104680426A (zh) | 基于伊藤微分的时滞电力系统随机稳定性分析方法及系统 | |
CN106058897A (zh) | 一种基于相量的发电机强迫扰动源定位方法 | |
CN103605858A (zh) | 直线感应电机稳态特性分析方法 | |
CN106385053A (zh) | 风电极限穿透功率优化方法和装置 | |
Azizipanah‐Abarghooee et al. | Modeling DFIG‐based system frequency response for frequency trajectory sensitivity analysis | |
Xie et al. | Blended chaos control of permanent magnet linear synchronous motor | |
Masaud et al. | Modeling and analysis of self-excited induction generator for wind energy conversion | |
Ouambo et al. | Parameters and states estimation by moving horizon estimation, high gain observer and unscented Kalman filter of a doubly-fed induction generator driven by wind turbine: A comparative study | |
CN109546650A (zh) | 输电通道暂稳极限的估算方法及装置 | |
Thomas et al. | Investigations on transient behavior of an energy conservation chopper fed DC series motor subjected to a change in duty cycle | |
Bayat et al. | Modeling and linearization of DFIG based wind turbine | |
CN103279641A (zh) | 对除发电机外复杂电力系统进行多模态动态等值的方法 | |
CN105243214A (zh) | 一种互联电网强迫功率振荡节点敏感性评估方法 | |
CN107342599B (zh) | 一种机电暂态仿真中控制系统稳定性自动诊断和参数调整方法 | |
Quang et al. | Observer-based integral sliding mode control for sensorless PMSM drives using FPGA | |
Bonaccorso et al. | EKF-based MPPT control for vertical axis wind turbines | |
Sivakumar et al. | Investigations on the impacts of uncertain wind power dispersion on power system stability and enhancement through PSO technique |
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 |