CN105633947B - 一种upfc阻尼控制的定位方法 - Google Patents
一种upfc阻尼控制的定位方法 Download PDFInfo
- Publication number
- CN105633947B CN105633947B CN201510698346.8A CN201510698346A CN105633947B CN 105633947 B CN105633947 B CN 105633947B CN 201510698346 A CN201510698346 A CN 201510698346A CN 105633947 B CN105633947 B CN 105633947B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- msup
- 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.)
- Expired - Fee Related
Links
Landscapes
- Control Of Eletrric Generators (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种UPFC阻尼控制的定位方法,依据全系统线性化模型,应用分解阻尼矩阵分析法DTA进行UPFC稳定控制器的定位,使得DTA计算的结果与控制器向电力系统中发电机提供的阻尼转矩和发电机对指定机电振荡模态的参与性联系起来,物理意义更加清晰。本发明的有益效果为:具体的给出了控制器影响指定的电力系统机电振荡模态的过程,从而为选择控制器稳定信号提供物理意义清晰的指导。
Description
技术领域
本发明涉及控制技术领域,尤其是一种UPFC阻尼控制的定位方法。
背景技术
统一潮流控制器(UPFC)是迄今为止功能最为强大的FACTS装置,其主要由一台激励变压器(Excitation Transformer,ET)、一台增压变压器(Boosting Transformer,BT)、2个三相电压源逆变器(Voltage Source Converter,VSC)和一个直流连接电容器组成,其结构如图1所示。mE,mB和δE,δB分别为并联电压源逆变器和串联电压源逆变器的幅值调值比率和相位角,作为UPFC的控制信号输入,通过一定的控制策略分别完成潮流控制、电压调节等不同功能。UPFC每个控制回路都可装设附加阻尼控制器来抑制弱阻尼模态,这一特性为选择附加阻尼控制器回路提供了极大的灵活性,但同时阻尼信号也会极大影响稳定控制的效果。因此,稳定控制器理想阻尼信号的选择,是一个十分值得研究的问题。
现有技术中,对于选择和分析UPFC稳定控制器反馈信号通常使用的是源于现代控制理论的分析方法(MA:model analysis),因为它适合于大型系统计算和分析,编程简单,且利于软件的开发。其中用于分析UPFC阻尼控制定位对系统小干扰下振荡稳定性的影响主要是通过计算模态可控性指标、可观性指标和它们的乘积(残差)进行的。然后由残差Ri计算的模态分析给出的只是UPFC稳定控制器与指定机电振荡模态的数学关系,它实际上是一种黑箱方法。这种模态分析不能表达指标的来源、分配和传递的物理过程,即不能清晰地表示电力系统的物理过程,所以这种模态分析对于分析探索UPFC阻尼控制信号的选择,常常较为困难。
发明内容
本发明所要解决的技术问题在于,提供一种UPFC阻尼控制的定位方法,依据全系统线性化模型,应用分解阻尼转矩分析法(DTA)进行UPFC稳定控制器的定位,使得DTA计算的结果与控制器向电力系统中发电机提供的阻尼转矩和发电机对指定机电振荡模态的参与性(灵敏度)联系起来,物理意义更清晰。
为解决上述技术问题,本发明提供一种UPFC阻尼控制的定位方法,包括如下步骤:
(1)收集数据;收集发电机内部电抗数据、励磁系统数据,通过收据采集和监控系统SCADA系统、能量管理系统EMS获得电力系统稳态数据和静态数据;(2)计算包含UPFC的开环系统线性化矩阵;
式中,ΔmE,ΔδE,ΔmB,ΔδB为线性化后的UPFC输入控制信号,δ为发电机功角状态变量向量,ω为发电机转速状态变量向量,Eq为励磁电流空载电动势,Efd为发电机励磁电压,VDC为UPFC直流电容电压;Δ为线性化算子,变量加点为微分算子;M为发电机惯性常数对角矩阵,KA与TA分别是自动电压调节器的增益和时间常数,Td'0为励磁绕组时间常数;K1~K9与Kpd、Kqd、Kvd、Kpe、Kpde、Kpb、Kpdb、Kve、Kvde、Kvde、Kqb、Kqdb、Kvb、Kvdb为线性化系数;
将(1)式中的有关系数写成如下向量
将(1)用分块矩阵表示,则进一步写为:
式中,ω0Ι为对角阵,ω0为额定角速度;δ为发电机功角状态变量向量,ω为发电机转速状态变量,u为控制器输出信号,j为发电机个数;A21,A22,A23,A31,A32,A33,B2,B3为上述公式中的分块矩阵;Z为除了功角和转速之外的发电机状态变量,还包括UPFC自身的状态变量(不包括附加阻尼控制器的状态变量);
(3)计算全系统线性化模型;假设阻尼控制器传递函数为G(s),则
式中,y为输出变量,C为状态变量到反馈量y的传递函数;
联立式(3)和(4)即可得到全系统线性化方程:
其中,X为系统状态变量,A为系统线性化矩阵,B为控制矩阵,C为输出矩阵;
(4)计算前向通道函数Bj(s);根据状态方程,控制信号Δu到发电机机电振荡环节的前向通道函数为:
Bj(s)=A23(sI-A33)-1B3+B2 (6)
(5)计算阻尼转矩系数Dij;UPFC阻尼控制器针对第i个振荡模态向系统中第j台发电机提供的阻尼转矩系数Dij为:
Dij=M[Bj(λi)γj(λi)G(λi)],j=1,2,…N (7)
若定义,则式(7)可表示为:
其中,M为发电机惯性常数对角矩阵;γj(λi)为控制器输入反馈信号用各台发电机转速状态变量重构的重构系数;
(6)计算重构系数γj(λi);根据线性控制理论,输出信号y是状态变量的组合,即有:
y=γj(s)Δωj,j=1,2,…N (10)
(7)计算发电机灵敏度Sij;定义模态λi对第j台发电机转矩TDij的灵敏度参数来衡量影响转矩对模态的影响程度为Sij:
(8)计算DTA指标系统第i个振荡模态可表示为:
根据式(16),上式进一步化简为:
这表明UPFC控制器通过两种渠道向第i个振荡模态提供阻尼,首先通过各个机组振荡模态的参与性Sij对模态的阻尼起作用;还通过对各个机组的机电振荡提供阻尼转矩,因此DTA指标可定义为:
(9)选择UPFC阻尼控制器反馈通道;通过输出步骤(7)中得到的DTA指标并比较其大小即可选择出UPFC阻尼控制器反馈通道。
优选的,步骤(6)中计算重构系数的方法如下:根据线性控制理论得,
其中,m为状态变量的总个数,ai为状态量初值,则输出变量:
而状态变量:
因此可以得到:
假设只针对第i个振荡模态,则可以选取初值使得aj=0,(j=1,2...i-1,i+1,...m),则上式又可写为:
其中Vi为相应于λi的右特征向量,为Vi中对应于Δωj的分量,j为第j台发电机。
本发明的有益效果为:本发明的UPFC阻尼控制的定位方法,能够具体的给出控制器影响指定的电力系统机电振荡模态的过程,从而为选择控制器稳定信号提供物理意义方面清晰的指导。
附图说明
图1是本发明的统一潮流控制器UPFC的结构示意图。
图2是本发明的装有UPFC多机系统线性化Phillips-Heffron模型示意图。
图3是本发明的UPFC阻尼控制器传递函数框图。
图4是本发明的模态分析分解式的物理意义示意图。
图5是本发明的装有UPFC及其阻尼控制器的简单四机二区域电力系统结构示意图。
图6(a)、图6(b)、图6(c)、图6(d)是本发明的DTA指标计算结果物理意义示意图;图6(a)在调制信号Δmc上附加阻尼控制信号;图6(b)在调制信号Δmb上附加阻尼控制信号;图6(c)在调制信号Δδc上附加阻尼控制信号;图6(d)在调制信号Δδb上附加阻尼控制信号。
图7是本发明的阻尼稳定器对模态的传递单向通道示意图。
图8是本发明的调制通道选择仿真的结果。
具体实施方式
如图所示,一种UPFC阻尼控制的定位方法,包括如下步骤:
(1)收集数据;收集发电机内部电抗数据、励磁系统数据,通过收据采集和监控系统SCADA系统、能量管理系统EMS获得电力系统稳态数据和静态数据;
(2)计算包含UPFC的开环系统线性化矩阵;
式中,ΔmE,ΔδE,ΔmB,ΔδB为线性化后的UPFC输入控制信号,δ为发电机功角状态变量向量,ω为发电机转速状态变量向量,Eq为励磁电流空载电动势,Efd为发电机励磁电压,VDC为UPFC直流电容电压;Δ为线性化算子,变量加点为微分算子;M为发电机惯性常数对角矩阵,KA与TA分别是自动电压调节器的增益和时间常数,Td'0为励磁绕组时间常数;K1~K9与Kpd、Kqd、Kvd、Kpe、Kpde、Kpb、Kpdb、Kve、Kvde、Kvde、Kqb、Kqdb、Kvb、Kvdb为线性化系数;
将(1)式中的有关系数写成如下向量
如图2所示,为式(1)所示线性化模型的传递函数框图;
将(1)用分块矩阵表示,则进一步写为:
式中,ω0Ι为对角阵,ω0为额定角速度;δ为发电机功角状态变量向量,ω为发电机转速状态变量向,u为控制器输出信号,j为发电机个数;A21,A22,A23,A31,A32,A33,B2,B3为上述公式中的分块矩阵;Z为除了功角和转速之外的发电机状态变量,还包括UPFC自身的状态变量(不包括附加阻尼控制器的状态变量);
(3)计算全系统线性化模型;假设阻尼控制器传递函数为G(s),则
式中,y为输出变量,C为状态变量到反馈量y的传递函数;
联立式(3)和(4)即可得到全系统线性化方程:
其中,X为系统状态变量,A为系统线性化矩阵,B为控制矩阵,C为输出矩阵;
(4)计算前向通道函数Bj(s);根据状态方程,可得传递函数框图如图3所示,控制信号Δu到发电机机电振荡环节的前向通道函数为:
Bj(s)=A23(sI-A33)-1B3+B2 (6)
(5)计算阻尼转矩系数Dij;UPFC阻尼控制器针对第i个振荡模态向系统中第j台发电机提供的阻尼转矩系数Dij为:
Dij=M[Bj(λi)γj(λi)G(λi)],j=1,2,…N (7)
若定义,则式(7)可表示为:
其中,M为发电机惯性常数对角矩阵;γj(λi)为控制器输入反馈信号用各台发电机转速状态变量重构的重构系数;
(6)计算重构系数γj(λi);根据线性控制理论,输出信号y是状态变量的组合,即有:
y=γj(s)Δωj,j=1,2,…N (10)
根据线性控制理论,可得:
其中,m为状态变量的总个数,ai为状态量初值,则输出变量:
而状态变量:
因此可以得到:
假设只针对第i个振荡模态,则可以选取初值使得aj=0,(j=1,2...i-1,i+1,...m),则上式又可写为:
其中Vi为相应于λi的右特征向量,为Vi中对应于Δωj的分量,j为第j台发电机;
(7)计算发电机灵敏度Sij;定义模态λi对第j台发电机转矩TDij的灵敏度参数来衡量影响转矩对模态的影响程度为Sij:
(8)计算DTA指标系统第i个振荡模态可表示为:
根据式(16),上式进一步化简为:
这表明UPFC控制器通过两种渠道向第i个振荡模态提供阻尼,首先通过各个机组振荡模态的参与性Sij对模态的阻尼起作用;还通过对各个机组的机电振荡提供阻尼转矩,这一关系的物理意义如图4所示,因此DTA指标可定义为:
(9)选择UPFC阻尼控制器反馈通道;通过输出步骤(7)中得到的DTA指标并比较其大小即可选择出UPFC阻尼控制器反馈通道。
如图4所示,为装有UPFC及其阻尼控制器的简单四机二区域电力系统,DTA分析的内容是比较在UPFC四个控制信号me,mB,δe,δb上加设阻尼控制器的阻尼效果。
四机二区域电力系统的参数为:
稳态时,发电机1~4的输出功率均为900MW;母线7和9的负荷为967MW和1767MW;UPFC的参数为:Cdc=1,vdc0=1,xE=0.02,xB=0.005;发电机、输电线路详细参数分别见表1和表2:
表1四机两区系统的发电机参数
发电机 | TJ | Xd | Xq | X'd | T′d0 | D | KA | TA |
G1 | 13 | 1.8 | 1.7 | 0.3 | 8 | 0 | 20 | 0.5 |
G2 | 13 | 1.8 | 1.7 | 0.3 | 8 | 0 | 20 | 0.5 |
G3 | 12.35 | 1.8 | 1.7 | 0.3 | 8 | 0 | 20 | 0.5 |
G4 | 12.35 | 1.8 | 1.7 | 0.3 | 8 | 0 | 20 | 0.5 |
表2四机两区系统的线路参数
根据式(5)求得全系统线性化模型的模态为:
其中第一个模态为局部振荡模态(G1-G2);第二个模态为局部振荡模态(G3-G4);第三个模态为区域振荡模态(G1G2-G3G4),以下针对区域间模态λ3对区域模态的振荡进行分析;针对区域间模态λ3,根据式(6)计算控制信号到发电机机电振荡环节的前向通道函数Bi(s),计算结果列于表3;根据式(12)计算UPFC控制器的输入反馈信号用各台发电机转速状态变量重构的重构系数γj(λ3),j=1,2,3,4,计算结果列于表3;
表3 Bj(s)与γj(λ3)计算结果
根据式(8)计算计算结果列于表4;根据式(16)计算系统中各台发电机相应于指定振荡模态的灵敏系数Sij,计算结果列于表4;由步骤(6)和步骤(7)中得到的结果,由式(19)即可计算DTA指标计算结果列于表4,其结果的物理意义示意图如图4所示。
表4 DTA指标分析结果及验证
通过比较DTA指标大小可以发现,应该选择调制信号me以构成UPFC阻尼控制器信号的最优反馈通道。另外,以抑制四机两区域系统区域间振荡的效果为根据,四个调制信号作为反馈信号的优先级应依次为me、δe、mb、δb。
UPFC阻尼控制器对模态的阻尼传递通道如图7所示,可以利用相位补偿法对UPFC阻尼控制器参数进行整定,使得稳定器传递函数
参数的整定结果列于表5,系统仿真结果如图8所示。
表5 UPFC阻尼控制器参数结果
调制信号 | me | δe | mb | δb |
Kw | 3.6462 | 713.6518 | 34.3731 | 20.2430 |
T1 | 1.4147 | 0.09 | 0.4523 | 5.2422 |
T2 | 0.06 | 0.9910 | 0.09 | 0.04 |
T3 | 1.4147 | 0.09 | 0.4523 | 5.2422 |
T4 | 0.06 | 0.9910 | 0.09 | 0.04 |
显然,由图8所示,系统仿真结果与DTA指标计算结果得出的结论想一致,即UPFC阻尼控制器最优反馈信号是me。
尽管本发明就优选实施方式进行了示意和描述,但本领域的技术人员应当理解,只要不超出本发明的权利要求所限定的范围,可以对本发明进行各种变化和修改。
Claims (2)
1.一种UPFC阻尼控制的定位方法,其特征在于,包括如下步骤:
(1)收集数据;收集发电机内部电抗数据、励磁系统数据,通过收据采集和监控系统SCADA系统、能量管理系统EMS获得电力系统稳态数据和静态数据;
(2)计算包含UPFC的开环系统线性化矩阵;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mover>
<mi>&delta;</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mover>
<mi>&omega;</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>E</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>q</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>E</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>f</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>V</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<mi>I</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>D</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>K</mi>
<mrow>
<mi>p</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<msup>
<mi>T</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>d</mi>
<mn>0</mn>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mn>4</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<msup>
<mi>T</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>d</mi>
<mn>0</mn>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<msup>
<mi>T</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>d</mi>
<mn>0</mn>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<msup>
<mi>T</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>d</mi>
<mn>0</mn>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mrow>
<mi>q</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mi>A</mi>
</msub>
<msub>
<mi>K</mi>
<mn>5</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mi>A</mi>
</msub>
<msub>
<mi>K</mi>
<mn>6</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mi>A</mi>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>v</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>K</mi>
<mn>7</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msub>
<mi>K</mi>
<mn>8</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>9</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>&delta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>&omega;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;E</mi>
<mi>q</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;E</mi>
<mrow>
<mi>f</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;V</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>K</mi>
<mrow>
<mi>p</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>K</mi>
<mrow>
<mi>p</mi>
<mi>d</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>K</mi>
<mrow>
<mi>p</mi>
<mi>b</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>K</mi>
<mrow>
<mi>p</mi>
<mi>d</mi>
<mi>b</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<msup>
<mi>T</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>d</mi>
<mn>0</mn>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mrow>
<mi>v</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<msup>
<mi>T</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>d</mi>
<mn>0</mn>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mrow>
<mi>v</mi>
<mi>d</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<msup>
<mi>T</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>d</mi>
<mn>0</mn>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mrow>
<mi>q</mi>
<mi>b</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<msup>
<mi>T</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>d</mi>
<mn>0</mn>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mrow>
<mi>q</mi>
<mi>d</mi>
<mi>b</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mi>A</mi>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>v</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mi>A</mi>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>v</mi>
<mi>d</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mi>A</mi>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>v</mi>
<mi>b</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>K</mi>
<mi>A</mi>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>v</mi>
<mi>d</mi>
<mi>b</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>K</mi>
<mrow>
<mi>c</mi>
<mi>e</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>K</mi>
<mrow>
<mi>c</mi>
<mi>&delta;</mi>
<mi>e</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>K</mi>
<mrow>
<mi>c</mi>
<mi>b</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>K</mi>
<mrow>
<mi>c</mi>
<mi>&delta;</mi>
<mi>b</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;m</mi>
<mi>E</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;&delta;</mi>
<mi>E</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;m</mi>
<mi>B</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;&delta;</mi>
<mi>B</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,ΔmE,ΔδE,ΔmB,ΔδB为线性化后的UPFC输入控制信号,δ为发电机功角状态变量向量,ω为发电机转速状态变量向量,Eq为励磁电流空载电动势,Efd为发电机励磁电压,VDC为UPFC直流电容电压;Δ为线性化算子,变量加点为微分算子;M为发电机惯性常数对角矩阵,KA与TA分别是自动电压调节器的增益和时间常数,Td'0为励磁绕组时间常数;K1~K9与Kpd、Kqd、Kvd、Kpe、Kpde、Kpb、Kpdb、Kve、Kvde、Kqb、Kqdb、Kvb、Kvdb为线性化系数;
将(1)式中的有关系数写成如下向量
将(1)用分块矩阵表示,则进一步写为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mover>
<mi>&delta;</mi>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mover>
<mi>&omega;</mi>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mover>
<mi>Z</mi>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<mi>I</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>21</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>22</mn>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>A</mi>
<mn>23</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>31</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>A</mi>
<mn>32</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>A</mi>
<mn>33</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>&delta;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>&omega;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>B</mi>
<mn>3</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>&Delta;</mi>
<mi>u</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,ω0Ι为对角阵,ω0为额定角速度;δ为发电机功角状态变量向量,ω为发电机转速状态变量向,u为控制器输出信号,j为发电机个数;A21,A22,A23,A31,A32,A33,B2,B3为上述公式中的分块矩阵;Z为除了功角和转速之外的发电机状态变量,还包括UPFC自身的状态变量,不包括附加阻尼控制器的状态变量;
(3)计算全系统线性化模型;假设阻尼控制器传递函数为G(s),则
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>u</mi>
<mo>=</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>y</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>=</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;&delta;</mi>
<mi>j</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;&omega;</mi>
<mi>j</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,y为输出变量,C为状态变量到反馈量y的传递函数;
联立式(3)和(4)即可得到全系统线性化方程:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>A</mi>
<mi>&Delta;</mi>
<mi>X</mi>
<mo>+</mo>
<mi>B</mi>
<mi>&Delta;</mi>
<mi>u</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>=</mo>
<mi>C</mi>
<mi>&Delta;</mi>
<mi>X</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,X为系统状态变量,A为系统线性化矩阵,B为控制矩阵,C为输出矩阵;
(4)计算前向通道函数Bj(s);根据状态方程,控制信号Δu到发电机机电振荡环节的前向通道函数为:
Bj(s)=A23(sI-A33)-1B3+B2 (6)
(5)计算阻尼转矩系数Dij;UPFC阻尼控制器针对第i个振荡模态向系统中第j台发电机提供的阻尼转矩系数Dij为:
Dij=M[Bj(λi)γj(λi)G(λi)],j=1,2,…N (7)
若定义,则式(7)可表示为:
其中,M为发电机惯性常数对角矩阵;γj(λi)为控制器输入反馈信号用各台发电机转速状态变量重构的重构系数;
(6)计算重构系数γj(λi);根据线性控制理论,输出信号y是状态变量的组合,即有:
y=γj(s)Δωj,j=1,2,…N (10)
(7)计算发电机灵敏度Sij;定义模态λi对第j台发电机转矩TDij的灵敏度参数来衡量影响转矩对模态的影响程度为Sij:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
(8)计算DTA指标系统第i个振荡模态可表示为:
<mrow>
<msub>
<mi>&Delta;&lambda;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&Delta;D</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
根据式(16),上式进一步化简为:
UPFC控制器通过两种渠道向第i个振荡模态提供阻尼,首先通过各个机组振荡模态的参与性Sij对模态的阻尼起作用;还通过对各个机组的机电振荡提供阻尼转矩,因此DTA指标可定义为:
(9)选择UPFC阻尼控制器反馈通道;通过输出步骤(7)中得到的DTA指标并比较其大小即可选择出UPFC阻尼控制器反馈通道。
2.如权利要求1所述的UPFC阻尼控制的定位方法,其特征在于,步骤(6)中计算重构系数的方法如下:根据线性控制理论得,
其中,m为状态变量的总个数,ai为状态量初值,则输出变量:
<mrow>
<mi>y</mi>
<mo>=</mo>
<mi>C</mi>
<mi>X</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</msup>
<msub>
<mi>CV</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</msup>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
而状态变量:
<mrow>
<msub>
<mi>&Delta;&omega;</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</msup>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<msub>
<mn>2</mn>
<mi>j</mi>
</msub>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
因此可以得到:
<mrow>
<msub>
<mi>&gamma;</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>y</mi>
<mrow>
<msub>
<mi>&Delta;&omega;</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</msup>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</msup>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<msub>
<mn>2</mn>
<mi>j</mi>
</msub>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
假设只针对第i个振荡模态,则可以选取初值使得aj=0,(j=1,2...i-1,i+1,...m),则上式又可写为:
<mrow>
<msub>
<mi>&gamma;</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</msup>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</msup>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<msub>
<mn>2</mn>
<mi>j</mi>
</msub>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</msup>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</msup>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<msub>
<mn>2</mn>
<mi>j</mi>
</msub>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<msub>
<mi>V</mi>
<mrow>
<mi>i</mi>
<msub>
<mn>2</mn>
<mi>j</mi>
</msub>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Vi为相应于λi的右特征向量,为Vi中对应于Δωj的分量,j为第j台发电机。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510698346.8A CN105633947B (zh) | 2015-10-23 | 2015-10-23 | 一种upfc阻尼控制的定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510698346.8A CN105633947B (zh) | 2015-10-23 | 2015-10-23 | 一种upfc阻尼控制的定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105633947A CN105633947A (zh) | 2016-06-01 |
CN105633947B true CN105633947B (zh) | 2018-04-06 |
Family
ID=56048615
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510698346.8A Expired - Fee Related CN105633947B (zh) | 2015-10-23 | 2015-10-23 | 一种upfc阻尼控制的定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105633947B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108964540A (zh) * | 2018-07-27 | 2018-12-07 | 中国电力科学研究院有限公司 | 一种计算励磁附加调差对同步发电机同步转矩影响的方法及系统 |
CN110571831B (zh) * | 2019-09-16 | 2021-01-22 | 广东电网有限责任公司 | 一种计及新能源接入的主动配电网互联系统稳定控制方法 |
CN110829462B (zh) * | 2019-11-15 | 2023-07-14 | 国家电网公司西北分部 | 基于mfo算法的光伏高占比电网upfc应用优化方法及系统 |
CN116131284A (zh) * | 2023-01-29 | 2023-05-16 | 南京理工大学 | 一种电力系统超低频振荡分析的阻尼转矩分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007325359A (ja) * | 2006-05-30 | 2007-12-13 | Central Res Inst Of Electric Power Ind | 電力系統の制御系定数設定方法、装置及びプログラム |
CN101447679A (zh) * | 2008-09-17 | 2009-06-03 | 中国电力科学研究院 | 一种并行电力系统稳定器的实现方法 |
CN103956755A (zh) * | 2014-04-23 | 2014-07-30 | 国家电网公司 | 一种具有抑制超低频振荡能力的电力系统稳定器设计方法 |
CN104300555A (zh) * | 2014-10-13 | 2015-01-21 | 国家电网公司 | 一种抑制电力系统共振机理低频振荡的方法 |
-
2015
- 2015-10-23 CN CN201510698346.8A patent/CN105633947B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007325359A (ja) * | 2006-05-30 | 2007-12-13 | Central Res Inst Of Electric Power Ind | 電力系統の制御系定数設定方法、装置及びプログラム |
CN101447679A (zh) * | 2008-09-17 | 2009-06-03 | 中国电力科学研究院 | 一种并行电力系统稳定器的实现方法 |
CN103956755A (zh) * | 2014-04-23 | 2014-07-30 | 国家电网公司 | 一种具有抑制超低频振荡能力的电力系统稳定器设计方法 |
CN104300555A (zh) * | 2014-10-13 | 2015-01-21 | 国家电网公司 | 一种抑制电力系统共振机理低频振荡的方法 |
Non-Patent Citations (1)
Title |
---|
《利用多端统一潮流控制器阻尼多模态振荡》;胡晓波;《电力系统自动化》;20070910;第31卷(第17期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105633947A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105633947B (zh) | 一种upfc阻尼控制的定位方法 | |
US9300142B2 (en) | Method for emulation of synchronous machine | |
CN102150356B (zh) | 具有分量分离的直接功率控制 | |
Zhu et al. | Dual-loop control strategy for DFIG-based wind turbines under grid voltage disturbances | |
CN104716886B (zh) | 基于虚拟阻抗电压型变流器虚拟同步电机控制方法 | |
CN104102788B (zh) | 一种基于联合循环燃机系统模型的调速方法 | |
CN103197559A (zh) | 一种改善双馈机组接入后系统小干扰稳定性的方法 | |
EP2070174A2 (de) | Konditionierungseinrichtung für energieversorgungsnetze | |
CN101917003B (zh) | 一种电力系统小干扰振荡稳定性分解式模态分析方法 | |
CN101127445A (zh) | 基于非线性鲁棒电力系统稳定器的励磁控制方法 | |
CN115313524A (zh) | 一种基于构网型变流器的光伏发电并网控制方法及系统 | |
CN104808147B (zh) | 一种三相短路下双馈发电机仿真方法及其系统 | |
CN104993494B (zh) | 一种基于四象限电力电子变流器的电机模拟装置及方法 | |
CN103050992A (zh) | 风电并网对电力系统小干扰稳定影响的阻尼转矩分析方法 | |
CN108880300B (zh) | 一种基于双闭环控制的双馈风机整流器阻抗计算方法 | |
CN103972912B (zh) | 一种含风电电力系统频率响应的频域分析方法 | |
CN109921421A (zh) | 基于谐波电流传递函数的双馈风电机组输出谐波电流模型建立方法 | |
CN106786715A (zh) | 一种多端vsc‑hvdc系统下垂控制系数确定方法 | |
CN107124126A (zh) | 一种双馈感应发电机的无锁相环电流控制方法及装置 | |
Zhou et al. | A power hardware-in-loop based testing bed for auxiliary active power control of wind power plants | |
CN109217365A (zh) | 一种无刷双馈电机虚拟同步控制方法 | |
CN105140957B (zh) | 基于风电场和光伏电站聚合模型的机电振荡模式估算方法 | |
Almas et al. | RT-HIL testing of an excitation control system for oscillation damping using external stabilizing signals | |
Nayeripour et al. | An advanced analytical calculation and modeling of the electrical and mechanical harmonics behavior of Doubly Fed Induction Generator in wind turbine | |
CN107994565A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180406 Termination date: 20191023 |
|
CF01 | Termination of patent right due to non-payment of annual fee |