CN103217896A - 基于自由权矩阵方法的多facts抗时滞协调控制方法 - Google Patents
基于自由权矩阵方法的多facts抗时滞协调控制方法 Download PDFInfo
- Publication number
- CN103217896A CN103217896A CN2013101065314A CN201310106531A CN103217896A CN 103217896 A CN103217896 A CN 103217896A CN 2013101065314 A CN2013101065314 A CN 2013101065314A CN 201310106531 A CN201310106531 A CN 201310106531A CN 103217896 A CN103217896 A CN 103217896A
- Authority
- CN
- China
- Prior art keywords
- damping ratio
- tau
- individuality
- 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.)
- Granted
Links
Images
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).判断迭代次数是否已到达限值,如果是则输出最小阻尼比对应的个体,否则进行步骤6;
(5).判断个体的解是否满足约束条件,满足的,输出最小阻尼比对应的个体;否则返回步骤2;
(6).根据得到的个体设计抗时滞协调控制器。
本发明提供的优选技术方案中,在所述步骤1中,含多FACTS协调控制的电力***单时滞模型的状态方程:
其中,Ad=BKC;A、B、C是系数矩阵;K是待设计的多FACTS协调控制器系数矩阵。
本发明提供的第二优选技术方案中,在所述步骤3中,设定的约束条件为:使线性矩阵不等式(10)和(11)成立;
其中,
本发明提供的第三优选技术方案中,在所述步骤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)成立;
其中,
在所述步骤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
2[xT(t)N1+xT(t-d(t))N2]×
和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 (5)
1.一种基于自由权矩阵方法的多FACTS抗时滞协调控制方法,其特征在于,所述方法包括如下步骤:
(1).得到含多FACTS协调控制的电力***单时滞模型的状态方程;
(2).应用量子遗传算法生成个体;
(3).设定约束条件,求解满足约束条件的个体的特征值和最小阻尼比;
(4).判断迭代次数是否已到达限值,如果是则输出最小阻尼比对应的个体,否则进行步骤6;
(5).判断个体的解是否满足约束条件,满足的,输出最小阻尼比对应的个体;否则返回步骤2;
(6).根据得到的个体设计抗时滞协调控制器。
2.根据权利要求1所述的方法,其特征在于,在所述步骤1中,含多FACTS协调控制的电力***单时滞模型的状态方程:
其中,Ad=BKC;A、B、C是系数矩阵;K是待设计的多FACTS协调控制器系数矩阵。
3.根据权利要求1所述的方法,其特征在于,在所述步骤3中,设定的约束条件为:使线性矩阵不等式(10)和(11)成立;
其中,
4.根据权利要求1所述的方法,其特征在于,在所述步骤3中,求解满足约束条件的个体的特征值和最小阻尼比包括:
设定一个最小阻尼比的阈值,当各个体求得对应的***阻尼比超出此阈值时停止求解;要求解***的阻尼比必先求解其特征值λ;时滞***(5)的特征方程是:
det(λI-A-Ade-τλ)=0 (12)
方程(12)是超越方程,采用PDE离散化方法来近似求解。
5.根据权利要求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 true CN103217896A (zh) | 2013-07-24 |
CN103217896B 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) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103972881A (zh) * | 2014-03-20 | 2014-08-06 | 国网新疆电力公司 | 基于多facts结构的链式结构高压电网安全控制方法 |
CN103984239A (zh) * | 2014-02-13 | 2014-08-13 | 国家电网公司 | 一种基于wams的多facts协调控制数模混合仿真平台 |
WO2016026243A1 (zh) * | 2014-08-18 | 2016-02-25 | 国家电网公司 | 一种基于wams多facts协调控制的实时控制*** |
CN105977969A (zh) * | 2016-06-08 | 2016-09-28 | 山东大学 | 基于sod-lms的大规模多时滞电力***稳定性判别方法 |
CN103838965B (zh) * | 2014-02-26 | 2017-01-04 | 华北电力大学 | 基于广义特征值的时滞稳定上限计算***及其计算方法 |
CN106529013A (zh) * | 2016-11-03 | 2017-03-22 | 国网宁夏电力公司电力科学研究院 | 一种高压直流附加频率控制方法 |
Citations (4)
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 |
WO2006131017A2 (en) * | 2005-06-08 | 2006-12-14 | Abb Research Ltd | Damping electromechanical oscillations in power systems |
CN102570477A (zh) * | 2012-01-04 | 2012-07-11 | 河南省电力公司 | 基于线性化潮流算法模型的svc及tcsc混合布点规划方法 |
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 (4)
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 |
WO2006131017A2 (en) * | 2005-06-08 | 2006-12-14 | Abb Research Ltd | Damping electromechanical oscillations in power systems |
WO2013037846A1 (en) * | 2011-09-12 | 2013-03-21 | Alstom Technology Ltd | Sub-synchronous oscillation damping by shunt facts apparatus |
CN102570477A (zh) * | 2012-01-04 | 2012-07-11 | 河南省电力公司 | 基于线性化潮流算法模型的svc及tcsc混合布点规划方法 |
Non-Patent Citations (3)
Title |
---|
王伟岸: "考虑广域反馈信号时滞影响的互连电网区间阻尼控制", 《中国优秀硕士学位论文全文数据库》 * |
白碧蓉: "计及时滞影响的FACTS广域阻尼控制器的设计", 《中国优秀硕士学位论文全文数据库》 * |
黄柳强等: "FACTS协调控制研究进展及展望", 《电力***保护与控制》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103984239A (zh) * | 2014-02-13 | 2014-08-13 | 国家电网公司 | 一种基于wams的多facts协调控制数模混合仿真平台 |
CN103984239B (zh) * | 2014-02-13 | 2018-12-04 | 国家电网公司 | 一种基于wams的多facts协调控制数模混合仿真平台 |
CN103838965B (zh) * | 2014-02-26 | 2017-01-04 | 华北电力大学 | 基于广义特征值的时滞稳定上限计算***及其计算方法 |
CN103972881A (zh) * | 2014-03-20 | 2014-08-06 | 国网新疆电力公司 | 基于多facts结构的链式结构高压电网安全控制方法 |
CN103972881B (zh) * | 2014-03-20 | 2016-05-18 | 国网新疆电力公司 | 基于多facts结构的链式结构高压电网安全控制方法 |
WO2016026243A1 (zh) * | 2014-08-18 | 2016-02-25 | 国家电网公司 | 一种基于wams多facts协调控制的实时控制*** |
CN105977969A (zh) * | 2016-06-08 | 2016-09-28 | 山东大学 | 基于sod-lms的大规模多时滞电力***稳定性判别方法 |
CN105977969B (zh) * | 2016-06-08 | 2017-10-03 | 山东大学 | 基于sod‑lms的大规模多时滞电力***稳定性判别方法 |
CN106529013A (zh) * | 2016-11-03 | 2017-03-22 | 国网宁夏电力公司电力科学研究院 | 一种高压直流附加频率控制方法 |
CN106529013B (zh) * | 2016-11-03 | 2019-03-08 | 国网宁夏电力公司电力科学研究院 | 一种高压直流附加频率控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103217896B (zh) | 2016-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103217896A (zh) | 基于自由权矩阵方法的多facts抗时滞协调控制方法 | |
CN102636728B (zh) | 基于力矩分解法识别负阻尼低频振荡的方法 | |
Yu et al. | An unscented particle filtering approach to decentralized dynamic state estimation for DFIG wind turbines in multi-area power systems | |
Resende et al. | Application of dynamic equivalence techniques to derive aggregated models of active distribution network cells and microgrids | |
Liu et al. | Protection and control of microgrids using dynamic state estimation | |
CN102624012B (zh) | 基于力矩分解的积分算法识别负阻尼低频振荡的方法 | |
Wang et al. | Finite-time function projective synchronization control method for chaotic wind power systems | |
CN102386832A (zh) | 基于等价补偿自抗扰的发电机转子电流控制方法及*** | |
CN104485670A (zh) | 孤立电网中电压敏感性工业负荷时变阻尼特性的控制方法 | |
Hu et al. | Inertia estimation of renewable-energy-dominated power system | |
CN102570457B (zh) | 一种基于内网量测数据的两端口外网静态等值方法 | |
CN105787604B (zh) | 一种输配协同的配网态势快速感知方法 | |
Zhang et al. | Dynamic state estimation of a grid-connected converter of a renewable generation system using adaptive cubature Kalman filtering | |
CN103217623A (zh) | 与过渡电阻和负荷电流无关的线路双端故障测距方法 | |
CN108075478A (zh) | 一种暂态稳定紧急控制方法 | |
CN103199524B (zh) | 一种适应多种分布式电源接入的潮流计算方法 | |
CN104135017B (zh) | 一种基于振荡能量的失步解列判别方法 | |
Sajadi et al. | Plane Wave Dynamic Model of Electric Power Networks with High Shares of Inverter-Based Resources | |
Srivastava et al. | Decentralized dynamic state estimation of transmission lines using local measurements | |
Gopakumar et al. | Stability control of smart power grids with artificial intelligence and wide-area synchrophasor measurements | |
Abur | Dynamic-State Estimation | |
Emami et al. | Power system dynamic state estimation using particle filter | |
Liu et al. | A new linearization method of unbalanced electrical distribution networks | |
Abur | Observability and dynamic state estimation | |
CN102945299A (zh) | 使用受控电源等效发电机节点的emtdc仿真初始化方法 |
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 |