概率小扰动稳定性的蒙特卡罗模拟.docVIP

精品论文 参考文献 概率小扰动稳定性的蒙特卡罗模拟 罗 嘉 （四川电力公司遂宁公司 四川 遂宁 629000） 【摘　要】 本文基于蒙特卡罗概率模拟和特征分析方法，编制电力系统概率小扰动分析的MATLAB程序。程序中以系统特征值、阻尼比及参与因子等特征参数的样本均值、方差以及系统小扰动失稳概率为指标。用WECC3机9节点系统验证了程序的有效性。 【关键词】 电力系统，小扰动稳定，蒙特卡罗模拟，特征分析 中图分类号： Ｇ71 文献标识码： Ａ　　文章编号：ISSN1004-1621（2013）03-018-02 0 引言 电力系统几乎时时刻刻都受到小的扰动，如系统中负荷的小量变化；又如架空输电线因风吹摆动引起线间距离(影响线路电抗)的微小变化等等[1]。随着电力系统的发展，系统的小扰动稳定问题日益突出。 概率稳定分析弥补了确定性方法仅针对某些给定的运行条件进行分析的局??性，考虑了各种随机因素对系统稳定性的影响，能够更深刻地反映电力系统运行的特征和稳定性的实质，得到更有参考价值的信息，是对传统分析方法的重要突破和补充，具有重要的理论价值和工程应用意义。 1 小扰动稳定分析 描述电力系统动态特性的微分-代数方程组为 (1) 式中： x表示微分方程组中描述系统动态特性的状态向量； x为输入向量， y为描述代数方程组中系统的运行参量。 在稳定运行点( )线性化得 (2) 为维状态向量； 为维输出向量； 为维输入向量；A为 阶状态矩阵；B为 阶控制矩阵或输入矩阵；C为阶输入矩阵；D为阶前馈矩阵。 根据李雅普诺夫稳定性原理，非线性系统的小范围稳定性是由系统线性化特征方程的根即A的特征值所确定的：特征值实部全为负时，原始系统渐进稳定；至少存在一个正实部的特征值时，则系统是不稳定的。对A进行模态分析 ,则可获取更多关于系统模式、模式与状态变量的相关程度以及特征值灵敏度等方面的信息。 不同的时间、季节及气候会导致负荷组成和负荷水平的波动以及发电机出力分配的变化，同时输电元件的故障、计划检修等会导致网络拓扑发生变化。在系统的稳态运行点( )以及状态矩阵A、系统的小扰动稳定性也具有不确定性。电力系统扰动概率稳定分析的任务就是根据主要随机因素的统计特征来确定系统小扰动稳定性的概率特性。 2 基于蒙特卡罗方法的小扰动稳定仿真 2. 1 系统状态的蒙特卡罗模拟 考虑开机状态、负荷水平和网络拓扑的不确定性，蒙特卡罗状态抽样的基本方法如下： ① 发输电元件状态的抽取。假设元件失效是相互独立的，对元件 产生一个在[0，1]区间均匀分布的随机数 ，则两状态元件的状态 由式(3)确定。 (3) 式中， 为元件 的停运概率(不可用率)。对于输电线路，第 条线路的停运概率可通过下式计算： (4) 式中，分别为线路的故障率和修复率，单位均为：次/年。 当需要计及发电机的降额运行状态，设第 台发电机的降额运行概率为，则由式(5)确定发电机 的开机状态。显然，该抽样概念可方便地推广至具有多个降额状态的机组。 (5) ②负荷水平的抽取。假设节点的有功负荷 服从均值 、标准差 (以均值的百分数表示)的正态分布。节点的有功负荷水平由式(6)确定。 (6) 当需要考虑元件参数(例如线路电抗参数等)的不确定性时，一般也可应用正态分布随机变量和类似式(6)的公式来确定元件参数。 2. 2 小扰动稳定概率指标 将式(1)形成标准的线性化状态方程 (7) 则由常微分方程的稳定性理论可知，系统相应的特征方程式为 (8) 式中I为与A维数相同的单位矩阵。满足式(2.6)的 个解 便是A的特征值。若所有的特征值均有负实部，即为负实数或实部为负的复数，则式(8)所表述的系统在相应的稳定工作点上受小扰动时是稳定的；反之，若有一个或多个特征值有正实部，则系统是不稳定的。若存在特征值为零或纯虚数，而其他特征值均为负实部，则系统处于临界状态。 复数特征值总是以共轭对的形式出现，物理上把一对共轭复数特征值称为系统的一个振荡模式，把它