CN103217896B - 基于自由权矩阵方法的多facts抗时滞协调控制方法 - Google Patents
基于自由权矩阵方法的多facts抗时滞协调控制方法 Download PDFInfo
- Publication number
- CN103217896B CN103217896B CN201310106531.4A CN201310106531A CN103217896B CN 103217896 B CN103217896 B CN 103217896B CN 201310106531 A CN201310106531 A CN 201310106531A CN 103217896 B CN103217896 B CN 103217896B
- Authority
- CN
- China
- Prior art keywords
- tau
- damping ratio
- individual
- facts
- phi
- 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 36
- 238000013016 damping Methods 0.000 claims abstract description 35
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 13
- 230000002068 genetic effect Effects 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 31
- 230000008054 signal transmission Effects 0.000 abstract 1
- 241000036569 Carp sprivivirus Species 0.000 description 11
- 230000000875 corresponding effect Effects 0.000 description 7
- 230000001052 transient effect Effects 0.000 description 7
- 230000005284 excitation Effects 0.000 description 5
- 230000006872 improvement Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000003190 augmentative effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Feedback Control In General (AREA)
Abstract
本发明提供了一种基于自由权矩阵方法的多FACTS抗时滞协调控制方法,包括:(1).得到含多FACTS协调控制的电力系统单时滞模型的状态方程;(2).应用量子遗传算法生成个体;(3).设定约束条件,求解满足约束条件的个体的特征值和最小阻尼比;(4).判断迭代次数是否已到达限值,如果是则输出最小阻尼比对应的个体,否则进行步骤6;(5).判断个体的解是否满足约束条件,满足的,输出最小阻尼比对应的个体;否则返回步骤2;(6).根据得到的个体设计抗时滞协调控制器。通过本发明,可协调多台FACTS设备联合运行提高系统的动态稳定性;并且本发明对广域信号传输时滞具有较好的忍耐性,在较大的时滞范围内仍能发挥良好的协调作用。
Description
技术领域
本发明属于电力系统领域,具体涉及一种基于自由权矩阵方法的多FACTS抗时滞协调控制方法。
背景技术
当前广域测量系统(WAMS)的PMU装置以GPS为采样基准,能全网同步采集机组和线路的电压、电流以及重要的开关保护信号;并能计算得到电压和电流相量、频率和频率变化率、机组和线路功率、发电机内电势(功角)以及根据机组键相信号实测机组功角;同时还能提供扰动触发的暂态记录。
随着广域测量系统(wide‐area measurement system,WAMS)的逐步推广与完善,获取系统的部分状态变量和输出变量不再困难,而电力系统稳定本质上是一个全局问题,引入全局信息也有助于进一步提高全网稳定性。因此,结合WAMS技术进行协调控制研究也逐渐兴起,主要集中在系统阻尼控制和直流系统控制方面,与此相比采用WAMS技术的多FACTS协调控制研究较少。
发明内容
为克服上述缺陷,本发明提供了一种基于自由权矩阵方法的多FACTS抗时滞协调控制方法,利用量子遗传算法寻优获得多FACTS协调控制器的最佳增益;而且在多台SVC的多机算例中进行时域仿真证明了所提算法的有效性。
为实现上述目的,本发明提供一种基于自由权矩阵方法的多FACTS抗时滞协调控制方法,其改进之处在于,所述方法包括如下步骤:
(1).得到含多FACTS协调控制的电力系统单时滞模型的状态方程;
(2).应用量子遗传算法生成个体;
(3).设定约束条件,求解满足约束条件的个体的特征值和最小阻尼比最大化;
(4).判断迭代次数是否已到达限值,如果是则输出最小阻尼比最大化对应的个体,否则进行步骤(5);
(5).判断个体的解是否满足约束条件,满足的,输出最小阻尼比最大化对应的个体;否则返回步骤(2);
(6).根据得到的个体设计抗时滞协调控制器。
本发明提供的优选技术方案中,在所述步骤(1)中,含多FACTS协调控制的电力系统单时滞模型的状态方程:
其中,Ad=BKC;A、B、C是系数矩阵;K是待设计的多FACTS协调控制器系数矩阵。
本发明提供的第二优选技术方案中,在所述步骤(3)中,设定的约束条件为:使线性矩阵不等式(10)和(11)成立;
其中, τ>0;
本发明提供的第三优选技术方案中,在所述步骤(3)中,求解满足约束条件的个体的特征值和最小阻尼比包括:
设定一个最小阻尼比的阈值,当各个体求得对应的系统阻尼比超出此阈值时停止求解;要求解系统的阻尼比必先求解其特征值λ;时滞系统(5)的特征方程是:
det(λI-A-Ade-τλ)=0 (12)
方程(12)是超越方程,采用PDE离散化方法来近似求解。
本发明提供的第四优选技术方案中,在所述步骤(5)、(6)中,求得的最小阻尼比最大化对应的个体就是所需的控制器系数矩阵K。
与现有技术比,本发明提供的一种基于自由权矩阵方法的多FACTS抗时滞协调控制方法,采用WAMS的输出反馈信号,针对信号的传输时滞,使用自由权矩阵时滞稳定定理作为判据,同时保证时滞系统的最小阻尼比在一定阈值之上,利用量子遗传算法寻优获得多FACTS协调控制器的最佳增益。最后在多台SVC的多机算例中进行时域仿真证明了所提算法的有效性。
附图说明
图1为基于自由权矩阵方法的多FACTS抗时滞协调控制方法的流程图。
图2为SVC数学模型的控制框图。
图3为基于WAMS的集中式协调控制流程图。
图4为含协调控制量的SVC动态模型。
具体实施方式
一种基于自由权矩阵方法的多FACTS抗时滞协调控制方法,其改进之处在于,所述方法包括如下步骤:
(1).得到含多FACTS协调控制的电力系统单时滞模型的状态方程;
(2).应用量子遗传算法生成个体;
(3).设定约束条件,求解满足约束条件的个体的特征值和最小阻尼比最大化;
(4).判断迭代次数是否已到达限值,如果是则输出最小阻尼比最大化对应的个体,否则进行步骤6;
(5).判断个体的解是否满足约束条件,满足的,输出最小阻尼比最大化对应的个体;否则返回步骤2;
(6).根据得到的个体设计抗时滞协调控制器。
在所述步骤1中,含多FACTS协调控制的电力系统单时滞模型的状态方程:
其中,Ad=BKC;A、B、C是系数矩阵;K是控制器系数矩阵。
在所述步骤3中,设定的约束条件为:使线性矩阵不等式(10)和(11)成立;
其中, τ>0;
在所述步骤3中,求解满足约束条件的个体的特征值和最小阻尼比包括:
设定一个最小阻尼比的阈值,当各个体求得对应的系统阻尼比超出此阈值时停止求解;要求解系统的阻尼比必先求解其特征值λ;时滞系统(5)的特征方程是:
det(λI-A-Ade-τλ)=0 (12)
方程(12)是超越方程,采用PDE离散化方法来近似求解。
在所述步骤4、5中,求得的最小阻尼比最大化对应的个体就是所需的控制器系数矩阵K。
通过以下实施例对基于自由权矩阵方法的多FACTS抗时滞协调控制方法做进一步描述。
在基于自由权矩阵方法的多FACTS抗时滞协调控制方法中:
1.确定系统模型和状态空间方程
1)本方法所涉及动态元件主要包括发电机和FACTS元件。其中发电机采用4阶实用模型,并安装一阶快速励磁。第i台发电机可以描述为:
其中:δ表示发电机功角;ω表示发电机转子角速度;Hi表示发电机惯性时间常数的一半;Pm表示发电机机械功率;E′d表示发电机交轴暂态电动势;E′q表示发电机直轴暂态电动势;Id表示发电机定子电流直轴分量;Iq表示发电机定子电流交轴分量;D表示发电机转子阻尼系数;T′d0表示直轴开路暂态时间常数;
T′q0表示交轴开路暂态时间常数;xd表示直轴同步电抗;xq表示交轴同步电抗;x′d表示直轴暂态电抗;x′q表示交轴暂态电抗;Efd表示励磁电动势;TA表示励磁环节时间常数;KA表示励磁环节放大倍数;Vtref表示励磁参考电压;Vt表示机端电压。
2)本方法拟采用多台SVC为例阐述所提协调控制算法。SVC数学模型的控制框图如图2所示。
SVC的动态表达式可以考虑为:
其中,BSVC为SVC的等效输出电纳值;B1为中间变量;K为SVC控制器测量环节的增益;T1和T0表示时间常数;Vref为参考电压;Vt为SVC控制节点的测量电压。
3)本发明以多FACTS协调控制机制集中式协调为例,提出如图3所示的协调控制示意图。协调控制器利用WAMS数据平台同步测量系统多点状态量和输出量,进行计算后得到各台FACTS设备的控制指令并下发,实现FACTS间的协调控制。
对于每台FACTS设备来说,接收到的协调控制变量u将作为控制器输入信号的一部分,与本地控制变量一起送入控制器。以图1的SVC模型为例,叠加了协调控制变量u的控制器模型如图4所示。
由此,全系统的状态空间模型可写成:
其中,x为系统状态变量;u为输入FACTS的协调控制变量;y为系统输出变量;A、B、C是系数矩阵。
4)本文设计的多FACTS协调控制器采用系统输出变量y作为输入。获取远端输出信号y经过WAMS传输需要一定时间;而且控制变量u由协调控制器发送至FACTS设备也存在时滞。假设最终总的时滞大小为τ,则多FACTS的静态输出反馈协调控制器可设计为
u(t)=Ky(t-τ) (4)
式中,K待设计的多FACTS协调为控制器的系数矩阵。
5)最终,由式(3)和式(4)可以得到含多FACTS协调控制的电力系统单时滞模型的状态方程:
其中,Ad=BKC。
因此,设计抗时滞协调控制器的主要任务便是寻找合适的K,使得系统状态方程(5)保持稳定。
2、系统稳定性的判断依据:
本发明使用Lyapunov‐Krasovskii稳定性定理,它的主要思想是通过构造一个合适的Lyapunov‐Krasovskii泛函来获得使系统(5)稳定的充分条件。目前常采用的是含有二次型双积分项的泛函:
其中,P=PT>0,Q=QT≥0和Z=ZT>0是待定矩阵。对于任意合适维数的矩阵 有
其中,η1(t)=[xT(t)xT(t-d(t))]T。同时,根据牛顿—莱布尼茨公式,对于时滞项有
那么,根据自由权矩阵理论,对于任意合适维数的矩阵N1和N2,有
此时,计算V(t,xt)沿系统(5)的导数,并结合式(8)和式(9),可得如下定理:
定理1给定时滞τ>0,如果存在满足上述要求的矩阵P,Z,X,N1和N2,使得如下两个线性矩阵不等式(10)和(11)成立,则时滞系统(5)是渐进稳定的。
其中
定理1给出了判断时滞系统稳定的充分条件。在许多电力系统实际工程中,WAMS系统架设好后时滞τ通常是可以估计到的,此时再结合系统的状态方程,便可利用该定理作为系统稳定的判据。
3、利用寻优算法来获取合适的K
1)对系统的状态方程降节,本文采用Hankel降阶方法,它通过保证降阶前后系统的Hankel奇异值误差在较小范围内来实现降阶。
2)本文采用量子遗传算法寻求合适的控制器系数矩阵K。它的基本算法是:生成个体→测量→求解目标函数→与上一步最优值比较→生成量子旋转策略→转到第一步。
3)约束条件可设置为满足定理1中的式(10)和式(11)的解是否存在,如此来首先保证求出的系统是渐进稳定的。但可能会存在多组可行解,因此结合电力系统的实际需求,目标函数设定为使系统最小阻尼比最大化,即最小系统阻尼比越大,系统稳定水平越高,如此来保证系统能运行在强阻尼模式。同时为了确保控制的合理性,应设定一个最小阻尼比的阈值(如7%),当超出此阈值时可停止寻优。
要求解系统的阻尼比必先求解其特征值λ。时滞系统(5)的特征方程是
det(λI-A-Ade-τλ)=0 (12)
这一方程是超越方程,不能直接求解。因此本文采用PDE离散化方法来近似求解。由于系统(5)的时滞微分方程在相互关联的边界条件下可转换为一组双曲型偏微分方程(Hyperbolic Partial Differential Equation,H‐PDE),因此通过对PDE进行精细离散化得到的增广矩阵,其特征值可近似于式(12)的解。
4)抗时滞的协调控制算法流程如图1所示。
需要声明的是,本发明内容及具体实施方式意在证明本发明所提供技术方案的实际应用,不应解释为对本发明保护范围的限定。本领域技术人员在本发明的精神和原理启发下,可作各种修改、等同替换、或改进。但这些变更或修改均在申请待批的保护范围内。
Claims (4)
1.一种基于自由权矩阵方法的多FACTS抗时滞协调控制方法,其特征在于,所述方法包括如下步骤:
(1).得到含多FACTS协调控制的电力系统单时滞模型的状态方程;
(2).应用量子遗传算法生成个体;
(3).设定约束条件,求解满足约束条件的个体的特征值和最小阻尼比最大化;
(4).判断迭代次数是否已到达限值,如果是则输出最小阻尼比最大化对应的个体,否则进行步骤(5);
(5).判断个体的解是否满足约束条件,满足的,输出最小阻尼比最大化对应的个体;否则返回步骤(2);
(6).根据得到的个体设计抗时滞协调控制器;
在所述步骤(1)中,含多FACTS协调控制的电力系统单时滞模型的状态方程:
其中,Ad=BKC;A、B、C是系数矩阵;K是待设计的多FACTS协调控制器系数矩阵。
2.根据权利要求1所述的方法,其特征在于,在所述步骤(3)中,设定的约束条件为:使线性矩阵不等式(10)和(11)成立;
其中,
3.根据权利要求1所述的方法,其特征在于,在所述步骤(3)中,求解满足约束条件的个体的特征值和最小阻尼比包括:
设定一个最小阻尼比的阈值,当各个体求得对应的系统阻尼比超出此阈值时停止求解;要求解系统的阻尼比必先求解其特征值λ;时滞系统(5)的特征方程是:
det(λI-A-Ade-τλ)=0 (12)
方程(12)是超越方程,采用PDE离散化方法来近似求解。
4.根据权利要求1所述的方法,其特征在于,在所述步骤(4)、(5)中,求得的最小阻尼比对应的个体就是所需的控制器系数矩阵K。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310106531.4A CN103217896B (zh) | 2013-03-29 | 2013-03-29 | 基于自由权矩阵方法的多facts抗时滞协调控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310106531.4A CN103217896B (zh) | 2013-03-29 | 2013-03-29 | 基于自由权矩阵方法的多facts抗时滞协调控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103217896A CN103217896A (zh) | 2013-07-24 |
CN103217896B true CN103217896B (zh) | 2016-08-31 |
Family
ID=48815803
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310106531.4A Active CN103217896B (zh) | 2013-03-29 | 2013-03-29 | 基于自由权矩阵方法的多facts抗时滞协调控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103217896B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103984239B (zh) * | 2014-02-13 | 2018-12-04 | 国家电网公司 | 一种基于wams的多facts协调控制数模混合仿真平台 |
CN103838965B (zh) * | 2014-02-26 | 2017-01-04 | 华北电力大学 | 基于广义特征值的时滞稳定上限计算系统及其计算方法 |
CN103972881B (zh) * | 2014-03-20 | 2016-05-18 | 国网新疆电力公司 | 基于多facts结构的链式结构高压电网安全控制方法 |
CN104199347B (zh) * | 2014-08-18 | 2017-02-15 | 国家电网公司 | 一种基于wams多facts协调控制的实时控制系统 |
CN105977969B (zh) * | 2016-06-08 | 2017-10-03 | 山东大学 | 基于sod‑lms的大规模多时滞电力系统稳定性判别方法 |
CN106529013B (zh) * | 2016-11-03 | 2019-03-08 | 国网宁夏电力公司电力科学研究院 | 一种高压直流附加频率控制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5517422A (en) * | 1993-10-12 | 1996-05-14 | Massachusetts Institute Of Technology | Method and apparatus for direct control of the inter-area dynamics in large electric power systems |
CN102570477A (zh) * | 2012-01-04 | 2012-07-11 | 河南省电力公司 | 基于线性化潮流算法模型的svc及tcsc混合布点规划方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1737098A1 (en) * | 2005-06-24 | 2006-12-27 | Abb Research Ltd. | Damping electromagnetic oscillations in power system |
WO2013037846A1 (en) * | 2011-09-12 | 2013-03-21 | Alstom Technology Ltd | Sub-synchronous oscillation damping by shunt facts apparatus |
-
2013
- 2013-03-29 CN CN201310106531.4A patent/CN103217896B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5517422A (en) * | 1993-10-12 | 1996-05-14 | Massachusetts Institute Of Technology | Method and apparatus for direct control of the inter-area dynamics in large electric power systems |
CN102570477A (zh) * | 2012-01-04 | 2012-07-11 | 河南省电力公司 | 基于线性化潮流算法模型的svc及tcsc混合布点规划方法 |
Non-Patent Citations (3)
Title |
---|
FACTS协调控制研究进展及展望;黄柳强等;《电力系统保护与控制》;20120301;第40卷(第5期);第141-142页 * |
考虑广域反馈信号时滞影响的互连电网区间阻尼控制;王伟岸;《中国优秀硕士学位论文全文数据库》;20111215(第S2期);摘要,第31、33、39页 * |
计及时滞影响的FACTS广域阻尼控制器的设计;白碧蓉;《中国优秀硕士学位论文全文数据库》;20050915(第5期);摘要,第14、35、37页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103217896A (zh) | 2013-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103217896B (zh) | 基于自由权矩阵方法的多facts抗时滞协调控制方法 | |
Visconti et al. | Measurement-based load modeling using transfer functions for dynamic simulations | |
Nuqui et al. | Hybrid linear state estimation utilizing synchronized phasor measurements | |
Yu et al. | An unscented particle filtering approach to decentralized dynamic state estimation for DFIG wind turbines in multi-area power systems | |
CN106294993B (zh) | 一种考虑逆变器电流饱和的暂态能量函数分析方法 | |
US20120150497A1 (en) | Power Flow Analysis for Balanced Power Distribution Systems | |
CN101964525B (zh) | 一种支持大规模电流量测的配电网状态估计方法 | |
CN103401238B (zh) | 一种基于总体测辨法的电力负荷建模方法 | |
Zenelis et al. | Wide-area damping control for interarea oscillations in power grids based on PMU measurements | |
CN101957872A (zh) | 一种交直流电力系统的混合实时仿真方法 | |
CN102163844B (zh) | 基于相量测量装置的电力系统状态检测方法 | |
Feng et al. | A practical multi-phase distribution state estimation solution incorporating smart meter and sensor data | |
Resende et al. | Application of dynamic equivalence techniques to derive aggregated models of active distribution network cells and microgrids | |
CN104981954A (zh) | 电力系统控制系统以及其中所使用的分散控制器 | |
CN103972912B (zh) | 一种含风电电力系统频率响应的频域分析方法 | |
Li et al. | Application of ensemble Kalman filter in power system state tracking and sensitivity analysis | |
CN105403834A (zh) | 一种发电机动态状态估计方法 | |
Luan et al. | Estimation of state of charge for hybrid unmanned aerial vehicle Li-ion power battery for considering rapid temperature change | |
Ram et al. | Analysis of ZIP load modeling in power transmission system | |
Chouhan et al. | A literature review on optimal placement of PMU and voltage stability | |
CN103647284B (zh) | 一种解决单时间断面问题的电压稳定预测方法 | |
Ashrafidehkordi et al. | Improved accuracy of the power hardware-in-the-loop modeling using multirate discrete domain | |
CN102280884A (zh) | 一种电网等值方法 | |
CN103413031B (zh) | 一种基于线路电压稳定指标的连续潮流计算方法 | |
CN104300536B (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 |