王晓东,陈予恕

(哈尔滨工业大学 航天学院，哈尔滨 150001)

电力系统振荡是关系电力系统安全稳定运行的重要问题，除了正常运行的周期振荡外，电力系统还时有发生混沌振荡现象，其表现为非周期的、似乎无规则的、突发性或阵发性的机电振荡，在振荡严重的情况下会导致互联系统的解列[1]。为此，越来越多的科学家和工程技术人员认识到电力系统稳定性问题的重要性，并将其作为系统安全运行的重要方面在理论和工程技术应用领域进行研究[2-6]。

对于单机无穷大系统，Duan等[7]分析了带有一系列电容补偿器的单机无穷大系统的次同步谐振，应用分岔理论分析了该系统出现的复杂的周期轨道分岔，环面折叠分岔；Zhu等[8]研究了单机无穷大系统中的次同步谐振，用Hopf分岔定理分析了Hopf分岔现象和稳定的周期解；Alberto等[9]应用Melnikov方法计算单机无穷大系统中异宿轨道横截相交的条件；Wei等[10]研究了高斯白噪声对电力系统安全盆侵蚀的影响； Chen等[11]用数值方法研究了单机无穷大系统中混沌现象并对其进行控制；Zhang等[12]分析了一个非线性三参数电力系统的异宿分支，用Melnikov函数确定了发生混沌现象的阈值，获得了发生混沌现象的参数区域；Chen等[13]研究了一类电力系统的混沌和次谐振荡，给出了Melnikov函数的计算方法，并获得电力系统发生混沌振荡和次谐振荡的锥形参数区域；Nayfeh等[14]用多尺度方法分析了一种电力系统中的亚谐共振和混沌现象，并用数值方法得到了发生混沌现象的途径；Abed等[15]采用的修正单机无穷大模型考虑了系统状态对机组阻尼系数的影响，分析了Hopf分岔现象。

现有大多数文献[10-14]中同步发电机采用的是常数阻尼，即令阻尼系数D为常数。实际上，阻尼在系统正常运行的情况下是与系统的状态变量相关的，而且可能为正，也可能为负[15]。因此若D取用常值可能会改变系统的非线性动力学特性，不能真实反映系统的稳定性情况，影响了分析的准确性。

为此本文在文献[15]提出的修正单机无穷大模型基础上，建立了单机无穷大系统的数学模型，研究具有非线性阻尼时单机无穷大系统的动力学行为，运用多尺度法分析了周期负荷扰动下单机无穷大系统的主共振响应，求得了二阶近似周期解的解析表达式，得到了主共振响应的幅频曲线，应用C-L方法分析了不同参数区域内的局部分岔现象及其稳定性，通过数值仿真所得的结论，解释了同步发电机在不同系统参数下对应的实际运行工况，也为同步发电机的稳定同步运行提供定性理论指导。

1 单机无穷大电力系统的数学模型

图1 单机无穷大系统

单机无穷大系统模型如图1所示。 图1中“1”和“2”分别为同步发电机和无穷大系统，“3”和“4”为主变压器，“5”为负荷，“6”为断路器，“7”为系统联络线。设无限大母线电压U∠0°为电压参考相量，同步发电机的内电动势为E，δ为发电机转子q轴与同步参考轴间的相对角度，俗称“功角”，ω为相对角速度，发电机的惯性时间常数为TJ，ω0为同步参考轴线的同步角速度，记M=TJ/ω0，发电机的机械输入功率PM为恒定值，发电机到无限大母线之间的等效电抗为X，PE为发电机的电磁功率，则系统的数学模型为：

(1)

当发电机为隐极式时有：

考虑系统受到某周期性负荷扰动并考虑同步发电机的阻尼作用，则得：

(2)

式中:D为等值阻尼系数，PD为扰动功率幅值，Ω为扰动频率。

做无量纲化处理，令:

并记:

则方程(2)化为：

(3)

为理解单机无穷大电力系统的动力学行为，首先对于无扰系统即c=0,f=0时的系统进行分析，可得其相图如图2所示：

图2 0<ρ<1时无扰系统的相图

图2是不同初始条件的无扰系统的相轨迹，“1”代表小幅振荡轨道，“2”为联接两个平衡点的同宿轨道，为“1”与“3”的分界线，“3”代表旋转轨道。

记无扰系统的平衡点为E,P0，则其坐标分别为：

(x1,0)=(arcsinρ,0), (x2,0)=(π-arcsinρ,0)

易知E(x1,0)为中心，P0(x2,0)为鞍点。

随着参数ρ的增大，中心和鞍点的位置逐渐靠近，当参数ρ=1时，中心和鞍点重合在一起退化成为中心-鞍点(如图3所示)，此时对应同步发电机稳定与不稳定的临界点，系统处于静态稳定极限，同时表明静态稳定极限所对应的功角正好与最大电磁功率的功角一致。

图3 ρ=1时，无扰系统相图

考虑系统正常运行的情况下，同步发电机中阻尼与系统的状态变量相关性作用，文献[15]给出的修正单机无穷大模型考虑了系统状态对机组阻尼系数的影响，结合文献[16]给出的非线性阻尼模型(经验公式)，本文用vander Pol阻尼代替常系数阻尼，即令式(3)中c=ξ+γx2，则得:

(4)

2 小幅振荡的主共振分析

下面讨论周期性负荷扰动时系统的主共振响应，将式(4)中的非线性项sinx在中心点(x1,0)=(arcsinρ,0)处Taylor展开，由奇异性理论[17]中的识别条件，保留至二次非线性项得：

(5)

重新标度方程令:

ξ→ε2ξ,γ→ε2γ,α→εα,f→ε2f, 0<ε≪1为小参数，则方程化为：

(6)

现用多尺度法[18]寻求方程的渐近解：

x(t,ε)=x0(T0,T1,T2)+εx1(T0,T1,T2)+ε2x2(T0,T1,T2)

(7)

式中:T0=t,T1=εt,T2=ε2t，T0为快变时间尺度，T1,T2为慢变时间尺度。

研究主共振情况，因而设:

(8)

由多尺度法，比较ε同次幂的系数，得下列微分方程组：

(9)

(10)

(11)

方程(9)的复数形式的解为:

(12)

(13)

式中:cc表示等式右端函数的共轭复数部分，为使解中不出现久期项，须有:

D1A=0

由上式得：A=A(T2)，则方程(13)的特解为：

(14)

于是系统周期振荡响应的二阶近似解为:

将式(12),式(14)代入式(11)得：

(15)

式中:NST表示不产生久期项的项，消除式(15)中的久期项得：

(16)

设:

(17)

式中:a,φ为T2的实函数，将式(17)代入式(16)分离实部与虚部得：

(18)

图4 幅频特性曲线

(19)

从方程组(19)消去φ得频率响应方程：

(20)

讨论定常运动的稳定性，设:

(21)

其中:a0,φ0为式(20)的一组解，a1,φ1为小的扰动量，将式(21)代入式(18)并对a1,φ1展开，注意到a0,φ0满足式(19)并略去a1,φ1二阶及以上微量，得：

(22)

因而定常运动的稳定性依赖于式(22)的右端的系数矩阵的特征值。

利用式(19)，可得特征方程为：

(23)

由特征方程的特征根可知，在幅频曲线图像(图4)中，对于每一幅频曲线，具有铅垂切线的曲线之间的部分为不稳定，其余部分稳定，即当有多值解时，中间的一个解为不稳定的，其余两个为稳定的，其在该区间内有可能发生振幅的跳跃现象。

图5 响应幅值-激励幅值曲线

3 主共振的奇异性分析

根据C-L方法[17]对方程(20)作奇异性分析，令:

z=a2

(24)

得:

c3z3+c2z2+c1z+c0=0

(25)

其中:

记:

c3(z3+ε2z2+λz+ε1)=c3G(z,λ,ε1,ε2)

G(z,λ,ε1,ε2)=z3+ε2z2+λz+ε1

(26)

取λ为分岔参数，ε1,ε2为开折参数，显然G(z,λ,ε1,ε2)为芽g=(z,λ)=z3+λz的普适开折，且余维数为2。由奇异性理论[19]通过计算得：

图6 转迁集

图7 不同保持域内的分岔图

转迁集如图6所示，可以看出，转迁集将平面ε1-ε2分成四个不同的保持域，在不同区域内的分岔模式如图7所示，依之可调整系统参数达到选择工程中允许的动态行为模式的目的，因而可为同步发电机不同运行参数和结构参数下的动力学分析、分岔控制和工程优化设计提供理论指导。

4 数值分析

为更好的理解单机无穷大电力系统随着周期性负荷扰动的变化所呈现的动态特性，以下用数值仿真的方法讨论系统(4)的动力学行为变化情况。

图8 f不同时的分岔图

由图8可以看出，随着激励幅值f的增大，系统响应出现了倍周期分岔现象，轨道经历了周期1→周期2→混沌解的变化。当激励的幅值较小时，系统出现了周期1解，其相图如图9(a)所示，此时同步发电机可以保持同步运行，电力系统处于暂态稳定状态，随着外激励幅值f增大，当f=0.72时出现了周期2解，其相图如图9(b)所示；当f=0.725 9时出现了混沌解，此时发电机功角出现了不规则也不衰减的类似随机的混沌振荡，功角不再稳定可靠，系统失稳。图9(c)和图9(d)分别为混沌状态的相图和对应的混沌吸引子，由全局分岔图(图8)可见，混沌状态仅存在于扰动幅值f的一个很窄的数值范围内，当f=0.79时出现了旋转周期轨道，其相图如图9(e)所示，此时随着功角δ的增大，加速转矩不断增大，使发电机转子不断加速，发电机就不可能再恢复到同步运行，系统由于增幅振荡最终失去稳定，发电机失去同步后将在系统间产生功率和电压的强烈振荡，结果使发电机组和负荷被迫切除，甚至导致系统解列或瓦解。

图9 f不同时的相图及混沌吸引子

5 结 论

本文应用非线性动力学的解析方法，分析了在周期性负荷扰动下电力系统振荡的失稳现象，考虑了同步发电机非线性阻尼转矩的作用，比传统文献更贴近工程实际，由于保留了电力系统的非线性特征，因此所得结果能够更好地反应振荡失稳的机理，可以为电力系统的稳定运行提供依据，具体有：

(1)应用多尺度方法和C-L方法，对受周期负荷扰动的单机无穷大电力系统的主共振进行了研究，得到了系统小幅振荡周期解的二次解析表达式和幅频响应方程，分析了定常解的稳定性并得到了系统稳态响应分岔方程的转迁集和对应的分岔图，研究表明该系统的不同分岔模式与其运行参数和结构参数有密切联系，可为同步发电机的动力学分析、分岔控制和工程优化设计提供理论指导。

(2)电力系统在周期性激励下，当扰动激励幅值较小时，同步发电机功角具有小幅振荡周期解，数值结果与理论结果一致，同步发电机能够保持同步运行，扰动幅值达到某个临界值时，发电机功角出现不规则的振荡现象，系统处于混沌振荡状态，影响了系统的正常运行，同步发电机失步，电力系统失稳，扰动幅值若继续增大，发电机功角持续增大，同步发电机失步，造成系统解列。

参 考 文 献

[1]卢 强，梅生伟，孙元章.电力系统非线性控制[M].北京：科学出版社，1993.

[2]贾宏杰，余贻鑫，王成山.电力系统混沌现象及相关研究[J].中国电机工程学报，2001，21(7)：26-30.

JIA Hong-jie，YU Yi-xin，WANG Cheng-shan.Chaotic phenomena in power systems and its studies[J].Proceedings of the CSEE，2001,21(7):26-30.

[3]Gerardo E P,Paul M O,Arnau D C,et al.Output-feedback IDA stabilisation of an SMIB system using a TCSC[J].International Journal of Control,2010,83(12): 2471-2482.

[4]Shi J,Tang Y,Dai T,et al.Determination of SMES capacity to enhance the dynamic stability of power system[J].Physica C,2010,470 (20):1707-1710.

[5]Manjarekar N S,Banavar R N,Ortega R.Application of interconnection and damping assignment to the stabilization of a synchronous generator with a controllable series capacitor[J].International Journal of Electrical Power and Energy Systems,2010,32(1): 63-70.

[6]Kumaran R C,Venkatesh T G,Swaru K S.Voltage stability-Case study of saddle node bifurcation with stochastic load dynamics[J].International Journal of Electrical Power & Energy Systems,2011,33(8): 1384-1388.

[7]Duan X Z,Wen J Y,Cheng S J.Bifurcation analysis for an SMIB power system with series capacitor compensation associated with sub-synchronous resonance[J].Science in China.Series E: Technological Sciences,2009,52(2): 436-441.

[8]Zhu W,Mohler R R.Hopf bifurcations in a SMIB power system with SSR[J].IEEE Transactions on Power Systems,1996,11(3): 1579-1584.

[9]Alberto L F C,Bretas N G.Application of Melnikovs method for computing heteroclinic orbits in a classical SMIB power system model[J].IEEE Transactions on Automatic Control,2000,47(7): 1085-1089.

[10]Wei D Q,Zhang B,Qiu D Y,et al.Effect of noise on erosion of safe basin in power system[J].Nonlinear Dynamics,2010,61: 477-482.

[11]Chen H K,Lin T N,Chen J H.Dynamic analysis,controlling chaos and chaotification of a SMIB power system[J].Chaos,Solitons and Fractals,2005,24: 1307-1315.

[12]Zhang W N,Zhang W D.Chaotic Ocillation of a Nonlinear Power System[J].Applied Mathematics and Mechanics,1999,20(10): 1175-1183.

[13]Chen X W,Zhang W N,Zhang W D.Chaotic and subharmonic oscillations of a nonlinear power system[J].IEEE Transactions on Circuits and Systems Ⅱ ,2005,52(12): 811-815.

[14]Nayfeh M A,Hamdan A M A,Nayfeh A H.Chaos and stability in a power system: subharmonic-resonant case[J].Nonlinear Dynamics,1991,2: 53-72.

[15]Abed E H,Varalya P P.Nonlinear Oscillations in Power systems[J].International Journal of Electrical Power & Energy Systems,1984,6(1): 37-43.

[16]Liwschitz M M.Positive and negative damping in synchronous machines[J],AIEE Trans,1941,60: 210-213.

[17]陈予恕.非线性振动系统的分岔和混沌理论[M].北京：高等教育出版社，1993.

[18]陈予恕.非线性振动[M].北京：高等教育出版社，2002.