说明书基于多FACTS结构的链式结构高压电网安全控制方法
技术领域
本发明涉及电力系统中的发电机调度控制技术领域，具体地，涉及基于多FACTS结构的链式结构高压电网安全控制方法。
背景技术
电力系统是一个非线性、高维数、多目标、分层分布的动态大系统，保持其安全稳定运行对国民经济、电力企业和用户具有重大意义。随着电力系统互联规模的不断扩大，以及电力市场的快速发展和用户侧电能质量要求的显著提高，电力系统的运行环境和动态特性愈加复杂，对电网的安全稳定运行提出了更高的要求。
为了提高电力系统暂态稳定性、动态稳定性及电压稳定性，避免大停电事故的发生，保证电力可靠安全、优质经济地供应，安全稳定控制措施的作用愈加重要。作为安全稳定控制措施的重要手段，灵活交流输电技术（FACTS）在电力系统的应用越来越广泛，发挥的作用越来越显著。FACTS元件基于大功率电力电子技术，实时控制传输系统的电压、电流、电抗和相角，能够提高传输通道的输电能力和系统的电压支撑能力，增强系统的安全性和灵活性，改善系统的稳定性。因此，在新疆与西北主网联网750千伏第二通道工程中，FACTS元件得到了广泛应用，为提高联网后系统稳定起了重要作用。
新疆与西北主网联网750千伏第二通道工程，是国家电网公司实施能源集约化开发，建设能源电力大通道的重要工程，是继新疆、西藏与西北主网联网之后，进一步加强西北750千伏主网架建设的又一重大项目，也是“疆电外送”的重点项目，建成后将成为酒泉、哈密、柴达木地区风能和太阳能发电的重要外送通道，形成典型的链式结构高压网络，有力促进新疆电力开发外送。同时，将有效缓解青海电网缺电局面，提高柴达木循环经济试验区的供电能力和向西藏直流输电的可靠性。为满足链式结构高压大电网及新能源外送的动态无功补偿，该工程配备了多FACTS系统，该系统包括沙州站66kV侧360Mvar SVC、鱼卡站750kV侧330MVar磁控式母线高抗（magnetically controlled reactor，简称MCR）、沙鱼线路4*390Mvar阀控式高抗（valve-controlled reactors，简称VCR）。
多FACTS系统为新疆与西北联网系统的安全稳定运行提供了良好的支撑作用， 然而，如此大规模、多类型的FACTS设备的安装也给系统安全带来了一系列问题，其中最为显著的是多个FACTS元件间分散独立控制，未考虑多个FACTS元件的交互作用可能引发的控制不协调，导致运行效果弱化甚至系统失稳。
在实现本发明的过程中，发明人发现现有技术中至少存在可靠性低、协调性差和稳定性差等缺陷。
发明内容
本发明的目的在于，针对上述问题，提出基于多FACTS结构的链式结构高压电网安全控制方法，以实现可靠性高、协调性好和稳定性好的优点。
为实现上述目的，本发明采用的技术方案是：基于多FACTS结构的链式结构高压电网安全控制方法，主要包括：
a、建立基于暂态稳定的电力系统模型和动态元件模型；
b、基于电力系统模型和动态元件模型，针对链式结构高压电网，建立考虑多FACTS之间协调配合的协调优化模型；
c、利用遗传算法GA，求解协调优化模型的最优解。
进一步地，所述步骤b，具体包括：
通过对暂态稳定进行分析，建立电力系统暂态稳定的数学描述形式，结合动态元件的模型，建立协调优化模型和其约束条件。
进一步地，所述步骤c，具体包括：
在发生不同扰动的情况下，通过求解优化模型，确定如何根据扰动的严重程度与控制需求，优化不同FACTS设备的控制量，以满足系统稳态和暂态电压要求，同时确保系统功角稳定和电压稳定，在此前提下，使得付出的控制代价最小。
进一步地，所述求解优化模型的操作，具体包括：
目标函数的确立；各种约束条件的建立；通过遗传算法GA进行求解。
进一步地，所述目标函数的确立的操作，具体为：
min F(x,y,u)
s.t.Gi(x,y,u)=0i=1,2,3,...,NeqHj(x,y,u)<0j=1,2,3,...Nuequmin<uk<umaxk=1,2,3,...,Nuη(u1,u2,...,uNu)>ϵ]]>
其中F（x,y,u）为目标函数，Gi和Hi为等式和静态不等式约束，η为功角、电压稳定裕度约束；x是状态变量；y表示代数变量；u为控制变量；umax、umin分别是控制变量的上、下限；Neq、Nueq和Nu分别为等式约束、不等式约束以及控制变量的总数。
进一步地，所述通过遗传算法GA进行求解的操作，具体包括：
1）编码采用具有在大变异率下容易收敛，处理函数优化问题较为有效的实数编码；
2）初始化种群：采用小区间生成法，先把各待优化参数的取值范围分成群体总数个小区间，再在各小区间中分别随机地生成一个初始个体；
3）目标函数的选?。阂灾甘湍勘旰魑问媸兜哪勘旰?；
4）选择策略：采用随机均匀的选择方法；
5）交叉和变异：交叉方法选择Heuristic，变异方法选择Adaptive feasible。
进一步地，在步骤a中，所述基于暂态稳定的电力系统模型和动态元件模型，具体为：
电力系统暂态稳定性是指电力系统在遭受到大扰动后，能否过渡到新的或恢复到原来的稳定运行状态的能力；电力系统在大扰动下的稳定性，一般把系统过渡过程分为三个阶段，它们分别是：
故障前系统：
x.=f0(x,y,α)0=g0(x,y,α)t≤0-;]]>
故障中系统：
x.=fF,1(x,y,α)0=gF,1(x,y,α)0+≤t≤tcl,1-]]>
x.=fF,2(x,y,α)0=gF,2(x,y,α)tcl,1+≤t≤tcl,2-]]>
…
x.=fF,k(x,y,α)0=gF,k(x,y,α)tcl,(k-1)+≤t≤tF-;]]>
故障后系统：
x.=fPF(x,y,α)0=gPF(x,y,α)tF+≤t;]]>
x∈Rn是系统对时间连续的状态变量；y∈Rl是可发生突变的代数变量；α∈Rm是调度员可以调控的参数称为可调参数；矩阵满秩的假设能够保证电力系统可描述为常微分方程的形式。
本发明各实施例的基于多FACTS结构的链式结构高压电网安全控制方法，由于主要包括：建立基于暂态稳定的电力系统模型和动态元件模型；基于电力系统模型和动态元件模型，针对链式结构高压电网，建立考虑多FACTS之间协调配合的协调优化模型；利用遗传算法GA，求解协调优化模型的最优解；可以用于指导大规模FACTS设备的协调应用，确保大规?；チ缤陌踩榷?、高效经济运行；从而可以克服现有技术中可靠性低、协调性差和稳定性差的缺陷，以实现可靠性高、协调性好和稳定性好的优点。
本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。
下面通过实施例，对本发明的技术方案做进一步的详细描述。
具体实施方式
以下结合本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
针对链式结构高压电网，根据本发明实施例，提供了基于多FACTS结构的链式结构高压电网安全控制方法，尤其涉及一种针对链式结构高压电网，考虑多FACTS之间协调配合的安全控制方法。该基于多FACTS结构的链式结构高压电网安全控制方法，考虑多FACTS之间协调配合的安全控制方法不仅对大规模FACTS设备的协调应用具有重要的指导意义，而且对大规?；チ缤陌踩榷?、高效经济运行具有重要价值，用于指导大规模FACTS设备的协调应用，确保大规?；チ缤陌踩榷?、高效经济运行。
上述实施例的基于多FACTS结构的链式结构高压电网安全控制方法，主要包括：
步骤1：建立基于元件和控制装置详细动态模型的电力系统暂态稳定模型；
在步骤1中，基于元件和控制装置详细动态模型的电力系统暂态稳定模型为：
电力系统暂态稳定性是指电力系统在遭受到大扰动后，能否过渡到新的或恢复到原来的稳定运行状态的能力。电力系统对于这类扰动的表现为发电机潮流、相对功角、节点电压和其它状态变量相对于正常工况的大规模偏移。大扰动引发的电力系统暂态过程非常复杂，特别是在多重故障下，系统的结构和参数在扰动期间均发生变化。如果扰动导致的状态变量的运动轨线在故障后的电力系统稳定域内，那么系统是暂态稳定的。研究电力系统在大扰动下的稳定性，一般把系统过渡过程分为三个阶段，它们分别是：
故障前系统：
x.=f0(x,y,α)0=g0(x,y,α)t≤0-;]]>
故障中系统：
x.=fF,1(x,y,α)0=gF,1(x,y,α)0+≤t≤tcl,1-]]>
x.=fF,2(x,y,α)0=gF,2(x,y,α)tcl,1+≤t≤tcl,2-]]>
…
x.=fF,2(x,y,α)0=gF,2(x,y,α)tcl,1+≤t≤tcl,2-]]>
故障后系统：
x.=fPF(x,y,α)0=gPF(x,y,α)tF+≤t;]]>
x∈Rn是系统对时间连续的状态变量；y∈Rl是可发生突变的代数变量；α∈Rm是调度员可以调控的参数称为可调参数（一般地，该参数可唯一决定电力系统的暂态稳定性）。应该指出，矩阵满秩的假设可保证电力系统可描述为常微分方程的形式。
（1）发电机模型和机网接口：
在暂态稳定中，发电机的动态为Newton第二定律描述的转子运动方程和由Faraday定律等描述的电磁转换方程等，其电磁转换方程数学描述如下：
定子电压方程（由Park变换获得）
ud=ψ.d-ψqω-Raid]]>
uq=ψ.d+ψdω-Raiq]]>
u0=ψ.0-Rai0;]]>
上式中ψd（id）、ψq（iq）、ψ0（i0），分别为d轴和q轴磁链和电流；ω是转速，ud（uq）是d（q）轴的电势，Ra是电枢每相电阻。
（ii）转子各绕组电压方程：
efd=ψ.fd+Rfdifd]]>
0=ψ.1d+R1di1d]]>
0=ψ.1q+R1qi1q]]>
0=ψ.2q+R2qi2q]]>
上式中ψfd、ψ1d、ψ1q、ψ2q分别为励磁绕组f，d轴等值阻尼绕组D、q轴等值阻尼绕组q1和q2的磁链，Rfd、R1d、R1q、R2q以及ifd、i1d、i1q、i2q为上述四个绕组的等值电阻和电流，efd为励磁绕组f的电势。
（iii）定子磁链方程：
ψd=-(Lad+Ll)id+Ladifd+Ladi1d
ψq=-(Laq+Ll)iq+Laqi1q+Laqi2q
ψ0=-L0i0；
上式中Lad（Laq）、Ll是分别定子和转子d轴（q轴）的互感和漏感。（iv）气隙转矩：
Te=ψdiq-ψqid；
其转子运动方程为：
δ.=ω0Δωr;]]>
Δω.r=12H(Tm-Te-KDΔωr);]]>
其中，δ和Δωr分别是发电机在同步旋转坐标系下的转子角和转子角速度，Tm是机械力矩，Te是定义的电磁力矩。
（2）励磁系统模型：
专门为同步发电机提供励磁电流的设备，即与同步发电机转子电压的建立，调整以及必要时使其消失的有关的设备，成为励磁系统。励磁系统主要由励磁调节器和励磁功率单元组成。励磁功率单元向同步发电机的励磁绕组提供直流励磁电流；励磁调节装置AER，根据极端电压的电压变化控制励磁功率单元的输出，从而达到调节励磁电流的目的。
（i）励磁功率单元
在忽略限幅环节（饱和环节）作用时，典型地，励磁机（励磁功率单元）相应的励磁系统可以用基本的微分方程描述为：
TeE·fd=-Efd+KeUR,]]>
（ii）励磁调节器模型
发电机极端电压Ut经过测量环节与给定参考电压Uref作比较，其偏差进入电压调节器进行放大后，输出电压UR作为励磁机的励磁电压，以控制励磁机的输出电压，机发电机励磁电压Efd。为了励磁系统的稳定运行何改善其动态品质，引入励磁系统负反馈环节，即励磁系统稳定器，其一般为软反馈环节。US为励磁器的附加控制信号，往往是电力系统稳定器PSS的输出。
测量环节用一个时间常数为TR的惯性环节表示，由于TR极小，（一般为0～0.02s）故常被忽略。电压调节器通常有一个超前滞后环节和一个惯性放大环节表示。超前滞后环节反映了调节器的相位特性，由于TB和TC一般很小，可予以忽略。惯性放大环节放大倍速为KA，时间常数为TA。
其三阶系统模型是：
TeE·fd=-(KL+SE)Efd+UrTfU·F=-UF+KF/TL(UR-(KL+SE)Efd)TAU·R=-UR+KA(Uref-Ut+Us-UF);]]>
（3）原动机和调速系统模型：
电力系统中向发电机提供机械功率和机械能的装置，如汽轮机，水轮机等统称为原动机。为了控制原动机向发电机输出的机械功率，并保持电网的关于正常运行频率，以及在各自并列的发电机之间合理分配负荷。完成上述任务的系统，被称为调速系统。调速系统一般通过控制汽轮机的汽门开度或者水轮机的导水叶开度来实现功率和频率调节。通过改变调竖起的参数及给定值（一般给定速度或者给定功率）获得需要的发电机功率－频率特性。
（i）汽轮机
汽轮机是以一定温度和压力的水蒸气作为工质的叶轮式发动机，在数学上，用一阶惯性环节表示，写成微分方程形式为：
TmP·m=-Pm+μ;]]>
式中，μ为汽门开度，Pm为汽轮机的机械功率，Tm为反映蒸汽容积效应的时间常数。
（ii）水轮机
水轮机是以一定压力的水为工质的叶轮式发动机，水轮机的模型描写的是水轮开度μ和输出机械功率Pm之间的动态关系。若忽略引水管道的弹性，则刚性引水管道水锤效应可以简化如下：
ΔPmΔμ=1-Tws1+0.5Tws;]]>
其中分子的负号，表明系统是非最小相位系统，体现了水锤效应。
（iii）调速系统模型
简化的水轮机调速器在不考虑限幅等环节时，可以写成：
Tsμ=Kδ·(ωref-ω)-Ki(μ-μ0)-ξ1ξ·1=-ξ1/Ti+Kβμ·=-(Kβ/Ts+/Ti)ξ1-Kβ[Kδ(ω-ωref)+Ki(μ-μ0)]/Ts;]]>
（4）负荷模型：
（i）静态负荷模型
静态负荷模型反应了负荷的有功与无功在系统频率和电压缓慢变化而变化的规律，可用代数方程或者曲线表示。依赖于电压的静态负荷常表示为下式所示的指数函数形式：
P=P0(V‾)a]]>
Q=Q0(V‾)b;]]>
其中，P和Q是负荷的有功和无功；V是负荷所连母线电压的电压幅值； P0,Q0和V0分别为基准点运行时负荷的有功功率，无功功率和负荷母线电压幅值。
（ii）动态负荷模型
在系统电压和频率快速变化时，应考虑负荷的动态特性，并用微分方程表示，称之为动态负荷模型。一般电力系统动态负荷，约有60-70%均为负荷是电动机。本发明采用考虑机电暂态过程的感应电动机模型。具体如下：
同步坐标下电机定子绕组方程为：
Ux=E′x+X′Iy-raIx
Uy=E′y-X′Iy-raIy；
其中X′=Xs+XmXr/(Xm+Xr)为异步电动机暂态电抗。
上式可写成：
U=E′+(rs+jX′)I
转子的电压动态方程：
Td0′E·x′=-Ex′+(X-X′)Iy+sTd0′2πf0Ey′Td0′E·y′=-Ey′-(X-X′)Ix+sTd0′2πf0Ex′;]]>
上式可写成：
Td0′E·′=-E′-j(X-X′)I-jsTd0′2πf0E′;]]>
转子运动方程为：
Js·=Tm-Te;]]>
Te=Re(E′I*) ；
（5）FACTS模型：
柔性交流输电技术（Flexible AC Transmission System，简称FACTS）主要利用大功率电力电子元件构成的装置来控制或者调节交流电力系统的运行参数/网络参数从而优化电力系统的运行状态，提高电力系统的输电技术的能力。
（i）静止无功发生器（Static Var Compensator，简称SVC）
①TCR部分
②TSC部分
③总体基波分量分析：
QSVC=(ωC-2β-sin2βπωL)V2;]]>
（ii）静止同步补偿器STATCOM
STATCOM也称静止无功发生器（Advanced Static Var Generator，简称ASVG）。
交流侧输出的基波电压幅值为：
VASVG=KVCsinθ2;]]>
假设逆变器电压VASVG滞后于系统电压角度为δ，且电抗为y∠α，则逆变器从系统吸收的有功功率为： P=VsVASVGysin(δ+α)-VASVG2ysinα;]]>
STATCOM送入系统的无功功率为：
QASVG=VsVASVGycos(δ+α)-VASVG2ycosα;]]>
在稳态情况下，逆变器既不吸收也不发出有功功率。因此令P为零可知：
VASVG=Vssin(δ+α)sinα;]]>
因此有：
QASVG=Vs22rsin2δ,VC=Vssin(δ+α)Ksinαsin(θ/2);]]>
上式表明，若保持脉宽θ不变只调整相位角δ则可改变STATCOM向系统注入的无功功率，同时电容电压也随之改变。若同时调整脉宽θ和整相位角δ则可以使得电容电压保持常数而只调整无功功率。
（iii）晶闸管控制的串联电容TCSC：
XTCSC=KβXC；
XC=1/ωC ；
Kβ=1+2λ2π(λ2-1)[2cos2βλ2-1(λtanλβ-tanβ)-β-sin2β2];]]>
通常，考虑到经济因素，在TCSC中取ωL<1/ωC，且使得λ2为7左右。当β∈[0,π/2λ]时，TCSC为容性的，当β∈(π/2λ,π/2)时，TCSC是感性的。
（6）网络模型——潮流方程：
对于一个N个节点的交流电力系统，其极坐标系下的潮流方程如下：
PGi-PLi-Σj=1n(ViVjBijsinθij+Vi+VjGijcosθij)=0]]>
QGi-QLi-Σj=1n(ViVjGijsinθij+Vi+VjBijcosθij)=0;]]>
其中，PGi、QGi、PLi、QLi分别是节点i上的发电机有功、发电无功、负荷有功和负荷无功，Gij,Bij是节点i和节点j之间的电导、电纳；θij=θi-θj为节点i和节点j电压相角之差；Vi代表节点i的电压幅值。
步骤2：建立以控制代价最小、考虑系统暂态稳定约束和多FACTS之间协调配合的链式结构高压电网安全控制优化模型；
在步骤2中，以控制代价最小、考虑系统暂态稳定约束和多FACTS之间协调配合的链式结构高压电网安全控制优化模型：
对多类型、大容量FACTS协调控制措施的风险性和经济性进行综合评估和优化研究，以取得最优的混合协调控制措施。敏感严重扰动下的控制范围和控制量更大，可选择性更大，控制规律更加复杂，对于同一扰动，可能有更多的可行控制措施，需要对这些控制措施和控制对象进行综合优化研究，提出风险性和经济性综合优化方法，以取得更可行和控制代价更小的协调控制措施。
多FACTS协调控制优化问题可描述为：在发生不同扰动的情况下，如何根据扰动的严重程度与控制需求，优化不同FACTS设备的控制量，以满足系统稳态和暂态电压要求，同时确保系统功角稳定和电压稳定，在此前提下，使得付出的控制代价最小。其模型可描述为：
min F(x,y,u)
s.t.Gi(x,y,u)=0i=1,2,3,...,NeqHj(x,y,u)<0j=1,2,3,...Nuequmin<uk<umaxk=1,2,3,...,Nuη(u1,u2,...,uNu)>ϵ]]>
其中F（x,y,u）为目标函数，Gi和Hi为等式和静态不等式约束，η为功角、电压稳定裕度约束；x是状态变量，如发电机及励磁系统的快动态变量、潮流变量等；y表示代数变量，如节点电压、相角等；u为控制变量，这里选取常规串补、无功补偿和FACTS设备的控制量，例如并联无功补偿设备的电抗，可控串补的容性电抗，可控高抗的电抗，SVC的电抗，磁控式母线可控高抗的电抗等；umax、umin分别是控制变量的上、下限；Neq、Nueq和Nu分别为等式约束、不等式约束以及控制变量的总数。下面逐一详细介绍目标函数和约束条件。
1）目标函数
多FACTS协调控制优化问题的目标函数选取各控制变量的动作量成本最低，即：
min L(u)=uTRu ；
式中L（u）为成本函数，u为控制向量，包括并联无功补偿设备的电抗，可控串补的容性电抗，可控高抗的电抗，SVC的电抗，磁控式母线可控高抗的电抗等；控制变量u既包括连续变量，如串补的容性电抗以及SVC的电抗，也包括了离散变量，如并联电抗器的感性电抗。R为对角成本矩阵，该对角矩阵可以优化控制变量实施的先后顺序，比如投切并联电容器组，调节串补电容容量，调节SVC无功补偿量的成本系数分别设为5、4和3，就说明在整个控制变量实施过程当中，调节SVC无功补偿量是最低级别的，优先考虑实施，而投切电容器组是最高级别的在满足特定条件时才允许实施。
综上所述，多FACTS协调控制优化的目标函数是一个包含离散控制变量和连续控制变量的混成协调优化问题。
2）等式约束
多FACTS协调控制优化问题的等式约束除了包括常规的系统潮流约束，还包括短期动态平衡约束和网络约束，其表达式为：
f(x,y,u)=0
g(x,y,u)=0
h(x,y,u)=0；
其中x是微分状态变量，如发电机励磁系统动态变量以及负荷动态变量；y表示代数变量，如节点电压、相角等；u为控制变量，这里选取常规串补、无功补偿和FACTS 设备的控制量；f表示短期动态相关的方程，包括发电机动态以及励磁系统的动态；g表示系统潮流方程，具体包括：
①发电机节点
ΔPGi=PGi-UiΣj∈iUj(Gijcosθij+Bijsinθij)ΔQGi=QGi-UiΣj∈iUj(Gijsinθij-Bijcosθij);]]>
其中PGi、QGi分别表示发电机注入节点i的有功功率和无功功率，二者可通过联立求解发电机节点代数方程予以消去。
②负荷节点：
ΔPLi=PLi-UiΣj∈iUj(Gijcosθij+Bijsinθij)ΔQLi=QLi-UiΣj∈iUj(Gijsinθij-Bijcosθij);]]>
其中PLi、QLi分别表示负荷注入节点i的有功功率和无功功率，通过联立求解负荷代数方程可以消去这两个量。
③无能量注入节点：
ΔPi=UiΣj∈iUj(Gijcosθij+Bijsinθij)ΔQi=UiΣj∈iUj(Gijsinθij-Bijcosθij);]]>
3）静态不等式约束
静态不等式约束主要包括可调发电机有功出力的上、下限约束，可调发电机无功出力的上、下限约束，可调发电机机端电压幅值的上、下限约束，OLTC分接头的上、 下限约束，节点频率的上、下限约束以及支路传输功率的上限：
PGimin<PGi<PGimaxi∈SGQGimin<QGiQGimaxi∈SGVimin<Vi<Vimaxi∈SNTimin<TiTimaxi∈SOfimin<fi<fimaxi∈SFSij<Sijmaxi∈SL;]]>
其中为可调发电机i的有功出力上、下限；为可调发电机i的无功出力上、下限；为可调发电机i机端电压幅值的上、下限；为第i各OLTC分接头的上、下限；为节点i频率的上、下限；为支路ij传输功率的上限。
4）暂态功角不等式约束
取系统的惯量中心COI作为参考，各发电机转子相对摇摆角极限为：
|δi(tf)-δCOI(tf)|≤δmax;]]>
相应的暂态功角稳定裕度为：
η=1-|δi(tf)-δCOI(tf)|δmax≥0]]>
式中，δi（tf）为同步发电机i在tf时刻的功角，δCOI（tf）为发电机群在tf时刻的惯性中心角度，δmax为发电机转子相对摇摆角上限。
δCOI(tf)=Σj=SGδi(tf)Mj/Σj=SGMj;]]>
式中SG为系统发电机数，Mj为第j台发电机的转动惯量。
5）暂态电压不等式约束
系统电压稳定与负荷稳定直接相关，因而可采用负荷节点电压偏差的积分形式作为电压稳定的性能指标，其表达式为：
J(u)=∫t0tfψ(y,u,t)dt≤M.]]>
式中ψ(y,u,t)=[y(t)-yr(t)]TQ[y(t)-yr(t)]。Q是对角权重矩阵，y（t）是负荷节点t时刻的电压向量，yr（t）是负荷节点t时刻的参考向量，M是性能指标的目标值，具体表达式为其中VD为时间间隔tf-t0里y（t）的目标平均负荷节点电压偏差，Ny是y（t）里的负荷节点数目。该不等式保证了被控制的系统能够拥有期望的动态行为特征。由于紧急控制仿真系统轨迹的时间间隔是10s左右，因而通常采用10s，即tf-t0=10s。
步骤3：建立利用遗传算法GA求解协调优化模型。
在步骤3中，建立协调优化模型，求解目标函数使用GA遗传算法，具体步骤如下：
1）编码采用具有在大变异率下容易收敛，处理函数优化问题较为有效的实数编码。
2）初始化种群：采用小区间生成法，先把各待优化参数的取值范围分成群体总数个小区间，再在各小区间中分别随机地生成一个初始个体，这在一定程度上保证了初始群体在可行解域分布的广度，同时计算量也不大。
3）目标函数的选?。阂灾甘湍勘旰魑问媸兜哪勘旰?。
4）选择策略：采用随机均匀的选择方法。
5）交叉和变异：交叉方法选择Heuristic，变异方法选择Adaptive feasible。
本发明针对链式结构高压电网，考虑多FACTS之间协调配合提出了合理的安全控制方法，不仅对大规模FACTS设备的协调应用具有重要的指导意义，而且对大规?；チ缤陌踩榷?、高效经济运行具有重要价值。
下面结合西北新疆互联电网，对本发明的技术方案进行详细说明。应该强调的是，下 述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。
考虑多FACTS协调配合的链式结构高压电网控制方法实施步骤如下：
遗传算法步骤：
1）编码采用具有在大变异率下容易收敛，处理函数优化问题较为有效的实数编码。
2）初始化种群：采用小区间生成法，先把各待优化参数的取值范围分成群体总数个小区间，再在各小区间中分别随机地生成一个初始个体，这在一定程度上保证了初始群体在可行解域分布的广度，同时计算量也不大。
3）目标函数的选?。阂酝鹱钚〉淖钣懦绷髂Ｐ妥魑勘旰?。
4）选择策略：采用随机均匀的选择方法。
5）交叉和变异：交叉方法选择Heuristic，变异方法选择Adaptive feasible。
GA的最佳个体数设置为120，遗传代数设置为25。
针对具有链式结构的新疆与西北主网联网750千伏第二通道工程，分别在沙州站内和鱼卡站内设置协调控制器，以控制本站内部FACTS设备的协调动作。经优化后的控制策略及验证情况如下。
（1）稳态出力波动控制策略验证
利用2013年西北新疆联网夏大方式进行控制策略校核。优化前后，联网通道各变电站的无功补偿配置如下表1所示。感性为“+”，容性为“-”。
表1：联网通道各变电站的无功补偿配置
根据优化前后的无功补偿投入量可以看出，采用前述FACTS系统协调控制策略可大大减少所需无功补偿设备投入量。为进一步说明采用优化后的无功补偿配置能够满足扰动情况下系统安全稳定运行的要求，考虑风电从2400MW均匀波动至4000MW，每300MW一级，考虑SVC、磁控式母线可控高抗的跟随调压作用，沙州可控高抗、鱼卡可控高抗单独动作一级引起的电压变化如下表2、3所示。变电站中 压侧的电压变化在2.5%之内，电压变化合理。
表2：2013年夏大方式沙州可控高抗单独动作一级电压变化情况
表3：2013年夏大方式鱼卡可控高抗单独动作一级电压变化情况
根据以上计算结论可以看出，基于前述FACTS系统协调控制策略，风电出力大幅波动，近区母线电压不越限，同时沙州可控高抗、鱼卡可控高抗单独动作一级引起的电压变化在合理范围内，控制策略有效可行。
（2）新能源大规模脱网情况下控制策略校核
研究基于2013年冬大运行方式，即哈密直流投运4000MW，配套电源未建成的薄弱方式，模拟新能源大规模脱网的情况?？悸鞘芟抻贜-2暂稳的极限运行方式，此时敦煌和酒泉风电共接入4200MW，其中敦煌风电接入3200MW，酒泉风电计入1000MW。
2013年冬大极限方式下，初始方式SVC容量为0，沙州～鱼卡线路两侧四组可控 高抗均投入容量156Mvar（固定容量39Mvar+1级可控容量117Mvar）。设置酒泉风电场汇集侧玉门～嘉峪关三永N-1故障，模拟故障后甘肃4200MW风电全部脱网。风电脱网后，考虑SVC跟随调压作用，沙州和鱼卡最高电压仍达到830kV之上，启动FACTS系统协调控制策略，沙州站和鱼卡站两组线路可控高抗均投至最大容量390Mvar?？煽馗呖苟骱蟮缪鼓芄换指粗梁侠矸段?，沙州电压恢复至790kV，鱼卡电压恢复至790kV。动态无功补偿设备动作后引起的脱网故障后各站电压变化如下表4所示。
表4：2013年冬大极限方式下甘肃风电脱网后动态无功补偿设备动作引起各站电压变化表
由以上计算结果可以看出，联网通道近区风电大发方式，一旦发生风电连锁脱网事故，近区750kV母线稳态电压最高可达840kV以上，严重威胁设备安全与电网稳定运行，即使考虑SVC自动动作，鱼卡、柴达木站750kV母线电压仍达到835kV，电压严重越限；采用本项目紧急控制策略，最高电压控制在800kV，满足电网安全稳定运行要求，所提控制策略有效可行。
综上所述，本发明上述各实施例的基于多FACTS结构的链式结构高压电网安全控制方法，针对链式结构高压电网，考虑多FACTS之间协调配合的安全控制方法；该方法包括通过分析电力系统的暂态稳定模型和各动态元件模型，建立协调优化模型；并通过GA遗传算法对协调优化模型进行求解确定如何根据扰动的严重程度与控制需求，优化不同FACTS设备的控制量，以满足系统稳态和暂态电压要求，同时确保系统功角稳定和电压稳定，并在此前提下，使得付出的控制代价最小。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依 然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的?；し段е?。