CN103838965B - 基于广义特征值的时滞稳定上限计算***及其计算方法 - Google Patents
基于广义特征值的时滞稳定上限计算***及其计算方法 Download PDFInfo
- Publication number
- CN103838965B CN103838965B CN201410067208.5A CN201410067208A CN103838965B CN 103838965 B CN103838965 B CN 103838965B CN 201410067208 A CN201410067208 A CN 201410067208A CN 103838965 B CN103838965 B CN 103838965B
- Authority
- CN
- China
- Prior art keywords
- time
- lag
- time lag
- upper limit
- matrix
- 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
- 238000000205 computational method Methods 0.000 title abstract 2
- 239000011159 matrix material Substances 0.000 claims abstract description 87
- 238000000034 method Methods 0.000 claims abstract description 39
- 230000009467 reduction Effects 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 15
- 239000013598 vector Substances 0.000 claims description 12
- 230000006641 stabilisation Effects 0.000 claims description 11
- 238000011105 stabilization Methods 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 5
- 238000013016 damping Methods 0.000 description 16
- 230000004044 response Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 230000010355 oscillation Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 230000005284 excitation Effects 0.000 description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 3
- 238000005457 optimization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了电力***稳定分析技术领域中的一种基于广义特征值的时滞稳定上限计算***及其计算方法。***包括顺序相连的数据采集模块、时滞***处理模块、时滞上限求解模块和结果输出模块;方法包括:采集建立时滞***状态方程所需的网络结构参数、发电机频率和发电机功角;建立时滞***状态方程,并对时滞***状态方程中的参数矩阵进行降阶处理,得到降阶后的时滞***状态方程;生成基于改进自由权矩阵的时滞稳定判据;利用时滞稳定判据求解时滞稳定上限。本发明能够有效降低时滞稳定上限计算过程中的保守性,具有很好的正确性和有效性。
Description
技术领域
本发明属于电力***稳定分析技术领域,尤其涉及一种基于广义特征值的
时滞稳定上限计算***及其计算方法。
背景技术
基于广域信息的反馈控制器的设计存在时滞问题,而这必然会造成控制器的控制效果降低,甚至出现负阻尼的情况。因此,迫切需要对***的时滞稳定上限进行研究。
目前国内外学者在***时滞稳定上限的研究方面取得了大量有益成果,主要可以分为3类:时域法、频域法和直接法。时域法可判定***在特定场景下是否稳定,但在稳定程度、时滞稳定上限等信息的获取方面还需要进一步研究。频域法通过在实数空间中搜索时滞***的关键特征根,能够在一定程度上揭示时滞***变化规律,但计算量较大,求解速度有待提升。直接法借助Lyapunov理论和线性矩阵不等式(Linear Matrix Inequality,LMI)技术,可同时考虑时滞的随机变动、存在切换环节等情况,适用范围更为广泛,但该方法具有一定的保守性。
针对以上问题,本发明提出一种基于广义特征值的时滞稳定上限计算***及其计算方法。本发明从电力***实际出发,能够有效解决直接法求解时滞上限过程中保守性较高的问题。首先利用读入的广域信号数据建立时滞***状态方程,并在有效保留***中低频振荡成分的基础上降阶***。其次,该发明形成基于改进自由权矩阵的时滞稳定判据,在此基础上,将时滞稳定判据等价变换,利用广义特征值法求解***的时滞稳定上限。基于IEEE4机11节点***和IEEE16机68节点***仿真表明,本发明可以较好的降低传统方法的保守性,具有很好的有效性和正确性。
发明内容
本发明的目的在于,提供一种基于广义特征值的时滞稳定上限计算***及其计算方法,用于解决直接法求解时滞上限过程中保守性较高的问题。
为了实现上述目的,本发明提出的技术方案是,一种基于广义特征值的时滞稳定上限计算***,其特征是所述***包括顺序相连的数据采集模块、时滞***处理模块、时滞上限求解模块和结果输出模块;
所述数据采集模块用于采集建立时滞***状态方程所需的网络结构参数、发电机频率和发电机功角,并将采集的数据发送至时滞***处理模块;
所述时滞***处理模块用于建立时滞***状态方程,并对时滞***状态方程中的参数矩阵进行降阶处理;
所述时滞上限求解模块用于生成基于改进自由权矩阵的时滞稳定判据,并利用时滞稳定判据求解时滞稳定上限;
所述结果输出模块用于输出时滞稳定上限结果。
一种基于广义特征值的时滞稳定上限计算方法,其特征是所述方法包括:
步骤1:采集建立时滞***状态方程所需的网络结构参数、发电机频率和发电机功角;
步骤2:建立时滞***状态方程,并对时滞***状态方程中的参数矩阵进行降阶处理,得到降阶后的时滞***状态方程;
步骤3:生成基于改进自由权矩阵的时滞稳定判据;
步骤4:利用时滞稳定判据求解时滞稳定上限。
所述降阶后的时滞***状态方程为其中,x(t)为降阶后的电力***状态向量;
Α为降阶后的电力***状态矩阵;
Αd为降阶后的电力***时滞矩阵;
为降阶后的电力***状态量对应的状态值;
d(t)为时滞,0≤d(t)≤h且
h为时滞稳定上限且h>0;
μ为时滞最大变化率。
所述基于改进自由权矩阵的时滞稳定判据为:
其中,
Φ2=[N+M-N+S-M-S];
Ac=[AAd0];
N、M和S为改进自由权矩阵;
改进自由权矩阵N满足
改进自由权矩阵S满足
改进自由权矩阵M满足
?1(t)=[xT(t)xT(t-d(t))xT(t-h)]T;
P、Q、R、Z1和Z2为待定矩阵。
所述利用时滞稳定判据求解时滞稳定上限包括:
子步骤A1:将时滞稳定判据等价变换为
其中,Y1和Y2为附加矩阵,并且Y1=Y1 T≥0,Y2=Y2 T≥0,
子步骤A2:以v最小为目标,以时滞稳定判据等价变换后的矩阵不等式 以及 为约束条件,计算待定矩阵P、Q、R、Z1、Z2、Y1和Y2,进而得到时滞稳定上限h。
本发明利用改进自由权矩阵建立时滞***稳定判据,借助广义特征值法对***时滞上限进行求解,其能够有效降低时滞稳定上限计算过程中的保守性,具有很好的正确性和有效性。
附图说明
图1是基于广义特征值的时滞稳定上限计算***结构图;
图2是IEEE4机11节点***结构图;
图3是IEEE4机11节点SMA降阶前后特征根对比图;其中,(a)是全阶开环状态矩阵和降阶开环状态矩阵的特征根对比图,(b)是全阶闭环状态矩阵和降阶闭环状态矩阵的特征根对比图;
图4是IEEE4机11节点不同延时情况下发电机1-4相对功角动态响应曲线图;
图5是IEEE4机11节点不同延时情况下发电机2-3相对功角动态响应曲线图;
图6是IEEE4机11节点***G1与G4功角差各时滞时间下的阻尼比结果表;
图7是IEEE4机11节点***G2与G3功角差各时滞时间下的阻尼比结果表;
图8是IEEE16机68节点***结构图;
图9是IEEE16机68节点Schur降阶前后频率响应对比图;其中,(a)是全阶开环状态矩阵和降阶开环状态矩阵的特征根对比图,(b)是全阶闭环状态矩阵和降阶闭环状态矩阵的特征根对比图;
图10是IEEE16机68节点不同延时情况下发电机1-16相对功角动态响应曲线图;
图11是IEEE16机68节点不同延时情况下发电机3-14相对功角动态响应曲线图;
图12是IEEE16机68节点不同延时情况下发电机10-15相对功角动态响应曲线图;
图13是16机***G1与G16功角差各时滞时间下的阻尼比表;
图14是16机***G3与G14功角差各时滞时间下的阻尼比表;
图15是16机***G10与G15功角差各时滞时间下的阻尼比表。
具体实施方式
下面结合附图,对优选实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
实施例1
图1是本发明提供的基于广义特征值的时滞稳定上限计算***结构图。如图1所示,本发明提供的基于广义特征值的时滞稳定上限计算***结构图包括顺序相连的数据采集模块、时滞***处理模块、时滞上限求解模块和结果输出模块。
数据采集模块用于采集网络结构参数、发电机频率和发电机功角,并将采集的数据发送至时滞***处理模块。其中,网络结构参数是用于建立时滞***状态方程所需的参数,而发电机频率和发电机功角则是时滞***状态方程中的状态向量。
时滞***处理模块用于建立时滞***状态方程,并对时滞***状态方程中的参数矩阵进行降阶处理,得到降阶后的时滞***状态方程。
时滞上限求解模块用于生成基于改进自由权矩阵的时滞稳定判据;并对所述时滞稳定判据进行等价变换,利用广义特征值法求解***的时滞稳定上限。
结果输出模块用于输出时滞稳定上限结果。
本发明提供的基于广义特征值的时滞稳定上限计算方法包括:
步骤1:采集建立时滞***状态方程所需的网络结构参数、发电机频率和发电机功角。
网络结构参数包括时滞***中线路的阻抗、导纳、发电机的内阻抗和负荷的等值阻抗。发电机频率用于计算时滞***转速,发电机功角和时滞***转速的变化量为时滞***状态方程的状态向量。
步骤2:建立时滞***状态方程,并对时滞***状态方程中的参数矩阵进行降阶处理,得到降阶后的时滞***状态方程。
多输入多输出电力***的状态方程可表示为:
其中,x′(t)∈Rn为电力***状态向量,u′(t)∈Rm为电力***控制输入向量,A′∈Rn×n为电力***状态矩阵,B′∈Rn×m为电力***控制矩阵。
通过状态反馈后得到相应的闭环***为:
其中,C′为闭环状态矩阵,当***经过状态反馈时,闭环状态矩阵为C′=A′+B′K′1,其中,K′1∈Rm×n为各附加控制器的综合状态反馈矩阵。
在实际电力***中,控制输入向量通过SCADA/WAMS***向各控制器传达,信号传递过程必然存在一定时滞,则相应闭环***可描述为:
由式(3)可知,电力***时滞矩阵A′d=B′K′1。对于含时滞环节的***,其状态方程有如下形式:
其中,矩阵A′和A′d分别为电力***状态矩阵和电力***时滞矩阵,h为时滞稳定上限。公式(4)中时滞d(t)满足条件:
0≤d(t)≤h (5)
公式(4)-式(6)即为时滞***状态方程,μ为时滞最大变化率。
实际***分析中,仅对低频振荡特征根相应的特征向量中与***状态量Δω(发电机转速变化量)和Δδ(发电机功角变化量)相对应的元素感兴趣,以便了解在Δωi中所含的该振荡模式分量的相对幅值及相位。常用的***降阶方法包括SMA降阶方法和Schur平衡降阶方法。由于***降阶方法已经是本领域常用的方法,因此本发明只以选择模式分析法SMA为例,对***降阶做简单介绍。
对于***状态方程将其按下式划分
其中,X1=[ΔωT,ΔδT]为保留变量,X2为其他变量,待消去。
由式(7)可消去X2,得:
其中,I为单位阵,p为微分算子。
将上式改写为
其中,Xr=X1为保留变量,Ar(p)为运算形式的“降阶”***系数阵。
由式(7)-(9)可以得到两个重要性质:
(1)如果p=λ1(i=1,2,…,N)为式(7)相应***特征根,即|λiI-A|=0,则p=λi也为式(8)或(9)形式上降阶的***特征根,即亦有|λiI-Ar(λi)|=0,特征根不发生变化,***模式不变。
(2)对于原***,λi的特征向量ui,有Aui=λiui。设降阶***λi相应的特征向量为uri,即Ar(λi)uri=λiuri,则uri和ui中保留变量Xr相对应的元素相等,即特征向量的相应元素不变。因此,在Xr保留变量处去观察同一模式λi的振荡时,相对幅值即相位不变,或者说模态不变。这样,所关心频带输入输出特性被完整的保留下来。
通过降阶处理,可以得到降阶后的时滞***状态方程为:
公式(10)中,x(t)是降阶后的电力***状态向量,Α为降阶后的电力***状态矩阵,Αd为降阶后的电力***时滞矩阵,为降阶后的电力***状态量对应的状态值,d(t)、h和μ的含义同公式(4)且满足式(5)和式(6)。
步骤3:生成基于改进自由权矩阵的时滞稳定判据。
构造如下形式的Lyapunov-Krasovskii泛函(李雅普诺夫-克拉索夫斯基泛函):
公式(11)中,P=PT>0,Q=QT≥0,R=RT≥0和Zi=Zi T≥0(i=1,2)都是待定矩阵。
由Newton-Leibniz公式(牛顿-莱布尼兹公式)可知,对于改进自由权矩阵N、S和M,式(12)-式(14)成立:
其中,ζ1(t)=[xT(t)xT(t-d(t))xT(t-h)]T。
计算式(11)中V(x)关于t的导数为:
将式(10)中第一个式子及式(12)-式(14)代入式(15),并加入必要的松散项可得:
上述公式中,
Ac=[AAd0]。
N1、N2和N3为矩阵N的分块矩阵且N1的维度与PA+ATP+Q+R的维度相等,N2的维度与的维度相等,N3的维度与R的维度相等。
S1、S2和S3为矩阵S的分块矩阵且S1的维度与PA+ATP+Q+R的维度相等,S2的维度与的维度相等,S3的维度与R的维度相等。
M1、M2和M3为矩阵M的分块矩阵且M1的维度与PA+ATP+Q+R的维度相等,M2的维度与的维度相等,M3的维度与R的维度相等。
考虑到式(11)中Zi=Zi T≥0,(i=1,2),因此公式(18)-式(20)中θi=θi T≥0,(i=1,2,3)。若Φ+Φs≤0,则式(16)中对于Φ+Φs,利用Schur补后可得式(10)表征的时滞***稳定判据如下:
对于给定标量h>0和μ,若存在P=PT>0,Q=QT≥0,R=RT≥0和Zi=Zi T≥0(i=1,2), 和 使得如下线性矩阵不等式成立:
则对于同时满足条件式(5)和式(6)的时滞***(10)渐进稳定。
步骤4:利用时滞稳定判据求解时滞稳定上限。
式(21)表征的线性矩阵不等式仅能判定***是否稳定,而无法获取***时滞稳定上限等信息。考虑到广义特征值法能够求解优化问题的全局最小值,因此,本发明提出利用广义特征值法计算***的时滞稳定上限。由于式(21)不是标准的广义特征值形式,无法直接利用广义特征值法进行求解,因此,做如下变换,将式(21)同时左乘与右乘式(22),如下:
可得:
其中,v=1h,I为单位阵。
为了求解最大时滞稳定上限h,即最小的v,本发明引入附加矩阵Yi=Yi T≥0(i=1,2),该矩阵需满足式(24)成立:
再将式(24)代入式(23)可得:
由此,时滞稳定上限问题已经转化为以下优化问题:
以v最小为目标,以 和 为约束条件,计算待定矩阵P、Q、R、Z1、Z2、Y1和Y2,进而得到时滞稳定上限h。公式表示如下:
min v
P,Q,R,Zi,Yi
s.t.
通过求解式(26),可以计算出以矩阵P、Q、R、Z1、Z2、Y1和Y2为变量,以式(25)和式(26)为约束的最小v。最终,利用h=1/v可以求出时滞稳定上限。
实施例2
基于MATLAB仿真软件搭建的图2所示的IEEE4机11节点***,发电机采用6阶详细模型,励磁***采用快速励磁,基准模型下的负荷采用50%恒阻抗和50%恒电流模型。首先,通过模态分析法得到四机***的状态矩阵,并利用SMA方法分别对开、闭环状态矩阵进行降阶,如图3所示。
图3(a)是全阶开环状态矩阵和降阶开环状态矩阵的特征根对比图,(b)是全阶闭环状态矩阵和降阶闭环状态矩阵的特征根对比图。由图3可以看出,开、闭环状态矩阵在利用SMA进行降阶后,均保留了***中低频振荡成分。因此,利用降阶后***矩阵能够求解***允许的最大时滞。将降阶后状态矩阵Α和时滞矩阵Αd代入式(26),求得最大时滞边界h=281.88ms。其中,降阶后状态矩阵Α如下:
降阶后时滞矩阵Αd如下:
时滞分别设置为h1=50ms,h2=100ms,h3=288.88ms和h4=350ms。所观察的物理量为发电机G1和G4之间的功角差,以及G2和G3之间的功角差,分别如图4和图5所示。
由图4和图5均可看出,在未加广域阻尼控制时,发电机之间的功角发生了严重的低频振荡。加入广域阻尼控制后,在不考虑时滞的情况下,低频振荡可得到有效抑制;然而随着时滞的增加,阻尼效果随之减弱,在最大时滞288ms的情况下,虽然仍存在一定的阻尼,但其阻尼比已经降到10%以下,说明此时控制器已经不满足控制要求。
利用prony算法对各时滞下的功角差曲线进行阻尼比分析,结果如图6和图7所示。
实施例3
基于MATLAB仿真软件搭建的图8所示的IEEE16机68节点***进一步考查基于广义特征值法的最大时滞求解方法的有效性和通用性。该***重要联络线为区域4和区域5的联络线1-2,1-27和8-9。发电机采用6阶详细模型,励磁采用IEEE-DC1型励磁,负荷模型15%恒有功功率,25%的恒有功电流和15%的恒无功功率,25%的恒无功功率和60%恒阻抗。首先利用Schur平衡降阶方法,在保证所关心频带的输入输出特性保持不变的前提下,对***进行降阶,如图9所示。
图9(a)中实线为开环全阶***的频率响应,虚线为开环降阶***的频率响应,可以看出,降阶***和全阶***的输入输出特性相同。图9(b)中实线为闭环全阶***的频率响应,虚线为闭环降阶***的频率响应,可以看出,闭环降阶***同样保留了全阶***的特性。因此,利用降阶后的***矩阵求解全阶***的时滞稳定上限具有有效性和可行性。
降阶后的状态矩阵Αre与降阶后的时滞矩阵Αdre代入式(26),可得最大的时滞上限为h=93ms。时滞分别设置为h1=50ms和h2=93ms。以发电机G1和G16之间的功角差,G3和G14之间的功角差与G10和G15之间的功角差,为所观察的物理量,分别如图10,图11和图12所示。其中,降阶后的状态矩阵Αre如下:
降阶后的时滞矩阵Αdre如下:
随着时滞的增加,阻尼效果随之减弱,在最大时滞93ms的情况下,虽然存在一定的阻尼效果,但是其阻尼比已经降到10%以下,此时控制器同样也不能满足控制的要求。利用prony方法对各时滞下的功角差曲线进行阻尼比分析,结果如图13、图14和图15所示。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (3)
1.一种基于广义特征值的时滞稳定上限计算方法,其特征是所述方法包括:
步骤1:采集建立时滞***状态方程所需的网络结构参数、发电机频率和发电机功角;
步骤2:建立时滞***状态方程,并对时滞***状态方程中的参数矩阵进行降阶处理,得到降阶后的时滞***状态方程;
步骤3:生成基于改进自由权矩阵的时滞稳定判据;
步骤4:利用时滞稳定判据求解时滞稳定上限;
所述降阶后的时滞***状态方程为
其中,x(t)为降阶后的电力***状态向量;
Α为降阶后的电力***状态矩阵;
Αd为降阶后的电力***时滞矩阵;
为降阶后的电力***状态量对应的状态值;
d(t)为时滞,0≤d(t)≤h且
h为时滞稳定上限且h>0;
μ为时滞最大变化率;
所述基于改进自由权矩阵的时滞稳定判据为:
其中,
Φ2=[N+M -N+S -M-S];
Ac=[A Ad 0];
N、M和S为改进自由权矩阵;
改进自由权矩阵N满足
改进自由权矩阵S满足
改进自由权矩阵M满足
P、Q、R、Z1和Z2为待定矩阵。
2.根据权利要求1所述的计算方法,其特征是所述利用时滞稳定判据求解时滞稳定上限包括:
子步骤A1:将时滞稳定判据等价变换为
其中,Y1和Y2为附加矩阵,并且Y1=Y1 T≥0,Y2=Y2 T≥0,v=1/h;
子步骤A2:以v最小为目标,以时滞稳定判据等价变换后的矩阵不等式以及为约束条件,计算待定矩阵P、Q、R、Z1、Z2、Y1和Y2,进而得到时滞稳定上限h。
3.一种用于执行权利要求1所述方法的基于广义特征值的时滞稳定上限计算***,其特征是所述***包括顺序相连的数据采集模块、时滞***处理模块、时滞上限求解模块和结果输出模块;
所述数据采集模块用于采集建立时滞***状态方程所需的网络结构参数、发电机频率和发电机功角,并将采集的数据发送至时滞***处理模块;
所述时滞***处理模块用于建立时滞***状态方程,并对时滞***状态方程中的参数矩阵进行降阶处理;
所述时滞上限求解模块用于生成基于改进自由权矩阵的时滞稳定判据,并利用时滞稳定判据求解时滞稳定上限;
所述结果输出模块用于输出时滞稳定上限结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410067208.5A CN103838965B (zh) | 2014-02-26 | 2014-02-26 | 基于广义特征值的时滞稳定上限计算***及其计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410067208.5A CN103838965B (zh) | 2014-02-26 | 2014-02-26 | 基于广义特征值的时滞稳定上限计算***及其计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103838965A CN103838965A (zh) | 2014-06-04 |
CN103838965B true CN103838965B (zh) | 2017-01-04 |
Family
ID=50802453
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410067208.5A Active CN103838965B (zh) | 2014-02-26 | 2014-02-26 | 基于广义特征值的时滞稳定上限计算***及其计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103838965B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104505830B (zh) * | 2015-01-14 | 2017-08-04 | 华北电力大学 | 时滞电力***稳定性分析方法和装置 |
CN104615882B (zh) * | 2015-02-03 | 2017-06-27 | 山东大学 | 基于eigd的大规模时滞电力***特征值计算方法 |
CN104680426B (zh) * | 2015-03-03 | 2018-06-22 | 华北电力大学 | 基于伊藤微分的时滞电力***随机稳定性分析方法及*** |
CN104932256B (zh) * | 2015-05-15 | 2018-04-17 | 河南理工大学 | 基于优化迭代算法的时滞广域电力***控制器 |
CN108808705B (zh) * | 2018-07-13 | 2020-05-22 | 山东大学 | 基于低阶sod-ps-ii-r算法的时滞电力***机电振荡模式计算方法 |
CN108808702A (zh) * | 2018-07-13 | 2018-11-13 | 山东大学 | 基于低阶igd-lms算法的时滞电力***机电振荡模式计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101350523A (zh) * | 2008-09-02 | 2009-01-21 | 天津大学 | 一种多时滞电力***稳定的判别方法 |
CN103217896A (zh) * | 2013-03-29 | 2013-07-24 | 国家电网公司 | 基于自由权矩阵方法的多facts抗时滞协调控制方法 |
CN103279035A (zh) * | 2013-05-20 | 2013-09-04 | 武汉大学 | 考虑wams信号时延的电力***广域输出反馈控制方法 |
-
2014
- 2014-02-26 CN CN201410067208.5A patent/CN103838965B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101350523A (zh) * | 2008-09-02 | 2009-01-21 | 天津大学 | 一种多时滞电力***稳定的判别方法 |
CN103217896A (zh) * | 2013-03-29 | 2013-07-24 | 国家电网公司 | 基于自由权矩阵方法的多facts抗时滞协调控制方法 |
CN103279035A (zh) * | 2013-05-20 | 2013-09-04 | 武汉大学 | 考虑wams信号时延的电力***广域输出反馈控制方法 |
Non-Patent Citations (5)
Title |
---|
Delay-dependent criteria for robust stability of time-varying delay systems;Min Wu et al.;《Automatica》;20041231;第40卷(第8期);第1435-1439页 * |
Improved Delay-Dependent Stability Criteria for Time-Delay Systems;Shengyuan Xu et al.;《IEEE Transactions On Automatic Control》;20050331;第50卷(第3期);第384-387页 * |
具有时变滞后的随机***的时滞依赖鲁棒稳定性与H∞分析;孙继涛 等;《应用数学和力学》;20100215;第31卷(第2期);第236-244页 * |
时滞电力***稳定性分析与网络预测控制研究;姚伟;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20101115(第11期);参见第3页,第21页第3段至第23页第2段,第36页第1段,第38页倒数第1段至第41页第3段 * |
时滞***的若干控制问题研究;曾启杰;《中国博士学位论文全文数据库 信息科技辑》;20130815(第8期);参见第35页第3段至第41页第1段 * |
Also Published As
Publication number | Publication date |
---|---|
CN103838965A (zh) | 2014-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103838965B (zh) | 基于广义特征值的时滞稳定上限计算***及其计算方法 | |
Borsche et al. | Effects of rotational inertia on power system damping and frequency transients | |
CN107800146B (zh) | 兼顾一次调频和超低频振荡抑制的调速器参数优化方法 | |
CN103886209B (zh) | 基于马尔科夫的跳变电力***时滞稳定性分析***及方法 | |
CN109962495B (zh) | 一种超低频振荡扰动源定位及抑制方法 | |
Yang et al. | Perturbation estimation based coordinated adaptive passive control for multimachine power systems | |
CN107240921A (zh) | 基于积分自适应逆推的svc滑模控制方法 | |
CN105048917A (zh) | 基于eso的双馈风力发电***积分滑模控制器的控制方法 | |
Gurung et al. | Optimal oscillation damping controller design for large-scale wind integrated power grid | |
WO2018145498A1 (zh) | 基于强化学习算法的双馈感应风力发电机自校正控制方法 | |
Swain et al. | Design of static synchronous series compensator based damping controller employing real coded genetic algorithm | |
CN106712055A (zh) | 一种与低励限制功能相协调的电力***稳定器配置方法 | |
CN108107720A (zh) | 基于状态空间分析的水轮机调速器参数整定方法及*** | |
CN108429501B (zh) | 一种永磁同步电机负载扰动的观测方法 | |
Bi et al. | Mode-based Damping Torque Analysis in Power System Low-frequency Oscillations | |
Patel et al. | A study on wide-area controller design for inter-area oscillation damping | |
Mandour et al. | Damping of power systems oscillations using FACTS power oscillation damper–design and performance analysis | |
Roy et al. | A nonlinear adaptive backstepping approach for coordinated excitation and steam-valving control of synchronous generators | |
Sadhana et al. | Revamped Sine Cosine Algorithm Centered Optimization of System Stabilizers and Oscillation Dampers for Wind Penetrated Power System | |
Orchi et al. | Feedback linearizing model predictive excitation controller design for two-axis models of synchronous generators in single machine power systems | |
Razali et al. | Power system stabilizer placement and tuning methods for inter-area oscillation damping | |
CN110556873B (zh) | 一种基于罚函数的vsg自适应转动惯量控制方法 | |
Ibrahim et al. | Modified particle swarm optimization based proportional-derivative power system stabilizer | |
Akkawi et al. | Comparative study between various controllers for power system stabilizer using particle swarm optimization | |
CN107248742B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |