电力系统暂态稳定性

1  电力系统暂态稳定性概述

1.1 电力系统暂态稳定及其意义

对于某一特定的稳定运行状态,以及对于某一特定的扰动,如果在扰动后系统可以达到一个可以接受的稳定运行状态,则对此初始状态及此扰动而言,称之为暂态稳定。

电力系统是一个复杂的动态系统,一方面它必须时刻保证必要的电能质量及数量;另一方面它又处于不断的扰动之中,扰动发生的时间、地点、类型、严重性均有随机性,扰动发生后的系统动态过程中一旦发生稳定性问题,系统可能在几秒内发生严重后果,造成极大的经济损失和社会影响。

电力系统暂态分析的主要目的是检查系统在大的扰动下(如故障、切机、切负荷、重合闸操作等情况),各发电机组间能否保持同步运行,如果能同步运行,并具有可接受的频率和电压水平,则称此电力系统在这一大扰动下是暂态稳定的。在电力系统规划、设计、运行等工作中都要进行大量的暂态分析。通过暂态分析还可以考察和研究各种稳定措施的效果以及稳定控制的性能,因此通过仿真来验证所求结果是否正确,即电力系统在某一状态时是否是稳定的具有重要意义。

电力线系统稳定的破坏,往往会导致系统的解列和崩溃,造成大面积停电,所以保证电力系统稳定是电力系统安全运行的必要条件。判定电力系统暂态稳定性主要是在大扰动下检查系统中各发电机组间能否保持同步运行水平,并具有可以接受的电压和频率水平。对这项工作已深入多年,并在离线计算中有成熟的算法取得了良好的成果。然而,随着电力市场化和区域联网的不断推进,电网运行状态越发复杂多变且接近其极限水平,同时,近年来国内外时有发生由于系统暂态稳定性不足且对策迟缓引起的大停电事故,造成了经济上的巨大损失,这些都提示着我们提高对电网稳定性分析的重视。只有更快速更准确地对系统处于的稳定水平进行判断,找到稳定性遭到威胁的环节才能进而有效地采取措施改善系统运行状态。为此,进一步深入研究改进原有暂态稳定分析方法,开拓新的创造性的方法,解决其在实际系统中的应用的难题,仍然是我们面对的重要课题。 

1.2 国内发展研究的现状

对于我国电网来说,其覆盖面积大,结构薄弱,负荷密度极不均匀,而电源有往往远离负荷中心。由于长期以来输电线路总长度年增长率比总装机容量增长率小得多,故进一步恶化了系统的安全稳定性。

电力系统的互联,可以带来显著地经济效益,但是长期以来,“分省平衡”的策略成为我国电力发展的重要弊端,严重地制约着我国电力资源的优化配置,全国联网的进程明显滞后。同时,电网的互联使得电力系统的规模变大,从而引起事故的可能性也越大。如果电网不够强壮,自动安全装置不够健全管理不得当,都有可能破坏系统的稳定,导致大面积停电,甚至全网崩溃。

以厂网分开为主要内容的电力体制改革实施后,我国电网建设的步伐明显加快,并且根据我国电网的特点和发展趋势,制定了“西电东送、南北互供、全国联网”的电网发展战略,大力推进跨区输电、跨区联网,其目标就是为了促进电力资源在更大范围内的优化配置。

截至到2006年,以三峡工程为核心,以华中电网为依托,向东南西北四个方向辐射联网的输电线路已基本建成。以北、中、南三大西电东送通道为主体南北网间多点互联、纵向通道联系较为紧密的全国电网互联的格局已基本形成预计到2010年,西电东送的规模将达到5500万kW;2020年将再增加到1亿kW以上。

2电力系统暂态稳定研究的内容和方法

2.1电力系统暂态稳定研究的内容

1.电力系统暂态(电磁和机电暂态)分析,电力系统暂态稳定性及提高措施的研究

2.掌握MATLAB的仿真工具simulink及power system 工具箱

3.运用MATLAB电力系统仿真程序psb构建模型

4.利用搭建的模型对电力系统暂态及提高稳定性的措施进行仿真研究

2.2研究方法

    研究暂态稳定的方法有许多种,下面简单介绍其中几种。

(1)时域仿真判定法

时域法是将电力系统各元件模型根据元件拓扑关系形成全系统模型,这是一组联立的微分方程组和代数方程组,然后以稳态工况或潮流解为初值,求扰动下的数值解,即逐步求得系统状态量和代数量随时间的变化曲线,并根据发电机功角值大于某一特定阀值来判别系统能否在大扰动后维持暂态稳定运行。

(2)能量函数判定法     

暂态能量函数法的理论基础是李亚普洛夫稳定性定理,因此也称为拟李亚普洛夫直接法(简称直接法)。该方法是通过比较扰动结束时暂态能量函数值与系统临界稳定时的暂态能量值来判别系统稳定性的。

(3)扩展等面积(EEAC)判定法

等面积法则从理论上较为完美地解决了单机-无穷大系统的暂态稳定评估问题,因此人们投入了大量的研究将其应用于多机系统中,其中以我国薛禹胜院士提出的扩展等面积法最为著名。扩展等面积法的基本思想是:在给定故障扰动下系统的机组分为临界机组和非临界机组两群,基于轨迹聚合和保稳变换,分别对这两群机组用等值机进行动态等效,进而又等值为单机-无穷大系统,从而可利用等面积法则评估系统的暂态稳定性。

(4)人工智能判定法

人工智能方法可进行非模型的电力系统暂态稳定判别,具有在线计算速度快、容易生成决策用的启发规则等优点,因此与传统暂态稳定分析方法构成了良好的互补。目前,人工智能方法主要有:模式识别、模糊理论、神经网络、支持向量机法。

2.3提高电力系统暂态稳定的方法

由于大扰动后,发电机机械功率和电磁功率的差额,即加速功率是导致系统暂态稳定破坏的主要原因。因此,提高暂态稳定,应从减小发电机转轴上的不平衡功率、减小转子相对加速以及减少转子相对动能变化量等方面着手。根据这种原则,提高电力系统暂态稳定的措施主要有:                                               

(1)自动重合闸装置:重合闸成功就可以增加减速面积,从而提高暂态稳定。

图2-1 自动重合闸对暂态稳定的影响

(2)快速切除故障:

减轻电气设备因短路电流产生的热效应等不良影响,而且加快切除故障可以减小切除角,这样既减小了加速面积,有增大了减速面积,对于提高暂态稳定起着决定作用。应该指出的是,切除故障时间是继电保护装置动作时间和开关接到跳闸脉冲到触头分开后电弧熄灭为止的时间总和。

                     

图2-2快速切除故障对暂态稳定性的影响

(3)变压器中性点经小阻抗接地:

电压器中性点经小电阻接地只对接地短路起作用。原因是短路时零序电流通过接地电阻时消耗有功功率,其中的一部分由发电机来负担,因而使发电机输出的电磁功率增加,从而减小了加速功率,提高暂态稳定。

(4)发电机采用电气制动:

在系统发生短路故障后,有控制地在发电机端投入电阻负荷,则可以增加发电机的电磁功率,产生制动作用从而提高暂态稳定。

2.4研究设计的内容

应用MATLAB 仿真进行电力系统仿真的主要步骤是:

(1)建立系统模型:建立的是单机无穷大系统(SIMB)。

单机—无穷大系统认为功率无穷大,频率恒定,电压恒定,是工程上最常用的手段,也是电力系统模拟仿真最简单、最基本的的运行方式,即对现实进行近似处理,以简化模型,更有利于得出结论,简化计算过程。

图2-5 无输电线的单机—无穷大系统原理图

假定联络阻抗为纯电感,则由发电机向无穷大系统送出去的有功功率的P为:

                                          

式中 分别为:包括发电机阻抗在内的发电机电动势到无穷大系统母线的总阻抗; 功角; 发电机电势; 系统母线电压。

由图1可得MATLAB单机无穷大仿真模型可能用到的模块:

PSB电力系统工具箱:

1.Electrical Sources中的Three-Phase Source(三相电源)模块

2.Elements 中的Three-Phase Parallel RLC Load(三相负载RLC并联)模块用来仿真输电线路,和Ground(交流接地)模块 ,Three-Phase Fault (三相故障整流器)模块用来仿真短路故障,Three-Phase Transformer(Two Windings)(三相变压器绕组)模块用来仿真变压器,

3.Machines里Synchronous Machine pu Standad(标么标准同步电机)模块用来仿真同步发电机

(2)设置仿真参数

(3)进行动态仿真

增加提高暂态稳定的模型,分别进行快速切除故障,变压器中性点经小电阻接地,改变强励倍数,记录结果。

(4)结果分析

   仿真之后,分析所得的图形。从电压图、转速图着手分析是否能够提高暂态稳定。

3电力系统常用仿真简介

3.1常用的电力系统仿真软件

目前常用的电力系统仿真软件有:

(1)邦纳维尔电力局开发的BPA程序和EMTP程序;

(2)曼尼托巴高压直流输电研究中心考法的PSCAD/EMTDC程序;

(3)德国西门子公司研制的电力系统仿真软件NETOMAC;

(4)中国电力科学研究院开发的电力系统分析综合程序PSASP;

(5)Math Works公司开发的科学与工程计算软件MATLAB。

以上各个电力系统仿真软件的结构和功能不同,它们各自的应用领域也有所侧重。EMTP主要用来进行电磁暂态过程数字仿真,PSCAD/EMTDC、NETOMAC主要用来进行电磁暂态和控制环节的仿真,BPA、PSASP主要用来进行潮流和机电暂态数字仿真。

3.2 MATLAB简介

MATLAB是美国MATH WORK公司1984年开始推出的一种简洁的工程计算语言,它是以矩阵计算为基础,把计算、可视化、程序设计融合到一个交互的工作环境中。

Simulink 是在MATLAB 环境下用于动态建模和仿真应用最广泛的软件包之一,意思是仿真链接。该软件有两个特别明显的功能:仿真与链接。也就是说,可以直接利用鼠标在模型窗口中划出所需要的控制系统模型,然后再利用该软件提供的功能来对控制系统直接进行模拟,很明显,这种所发是的一个原本很复杂的系统变得相当容易输入。

SimPowerSystems库是SIMULINK下面的一个专用模块库,是在SIMULINK环境下进行电力、电子系统建模和仿真的先进工具。SimPowerSystems库提供了一种类似电路建模的方式进行模型绘制,在仿真前自动将仿真系统图变化为状态方程描述的系统形式,然后在SIMULINK下进行仿真分析。它为电路、电力电子系统、电机系统、发电、输变电系统和配电计算提供了强有力的解决方法,尤其是当设计开发内容涉及控制系统设计时,优势更为突出。PSB(Power System Blocket)是一个图编辑器,在Simulink环境下能建立系统原理并进行仿真计算。PSB库提供了电力系统仿真通用的原件和装置,包括RLC支路和负载、变压器、传输线、避雷器、电机、电力电子装置等。只需通过点击和拖放PSB库内的模型即可建立用户所需要的电力系统仿真原理图,并利用模型元件的对话框来设置相关参数,使用Simulink提供的示波器模型,可现实仿真结果及其波形

4基于SIMULINK的单机无穷大系统的暂态稳定性仿真

4.1单机-无穷大系统的建模

 单机无穷大系统原理图如图4-1.

4-1单机-无穷大系统原理图

根据图4-1建立的单机无穷大系统仿真原理图如图4-2.

图4-2单机无穷大仿真原理图

图中使用到的基本模块有:

示波器scope,逻辑加,增益gain,电机测量模块Machines Measurement Demux,同步电机(标幺值单位)模块Synchronous Machine pu Standard,励磁系统Excitation System,通用电力系统稳定器Generic Power System Stabilizer,三相并联RLC负载Three  Phase  Parallel RLC Load,双绕组三相变压器Three-Phase Transformer(Two  Windings),电压测量模块Voltage Measurement,分布参数线路模块Distributed Parameters Line,三相电源Three-Phase Source,三相故障Three-Phase Fault,POWERGUI

运行经验表明,在110KV及以上的大接地电流系统的高压架空线路上,短路故障中70%以上短路故障都是接地短路,尤其在220-500KV的架空线路上,由于线间距离大,其中单相接地短路故障的几率高达90%以上,所以仿真时的故障采取的是单相接地短路。仿真模型中,发电机采用的是标幺制的同步电机,设置的有功是900MW;变压器采用的是双绕组三相变压器,有功也是900MW,频率均设置为我国电力的工频50Hz。

4.2采用的模块及其参数设置

  1. 同步发电机:

图4-3

凸极机,200MVA,13.8kV,60Hz,xd=1.305,xd′=0.296,xd″=0.252, xq=0.474 , xq′=0.243 ,xq″=0.18 , Rs=0.0028544, Td0′=1.01s , Td0″=0.053s ,Tq0″=0.1s,H=3.2s,p=32。其中各电抗电阻均为标幺值(p.u.);

  1. 变压器:

图4-4

210MVA, 60Hz , Delta(D1)/Yg ,13.8/230kV, R2=0.0027p.u.,L2=0.08p.u.,Rm=Lm=500p.u;

  1. 无穷大系统:

图4-5

采用三相电源模块代替,参数设置为:230kVA,A 相相角为0,X/R=10;

(4)励磁系统:

图4-6

Tr(s)=20e-3,励磁倍数Ke=1,Efmin=-11.5 ,Efmax=11.5 ;

(5)短路模块:

图4-7

设置为0.1s开始短路,10秒结束。

(6)断路器

图4-8

设置为0.2s开始切除故障线路,10秒结束。

1.快速切除故障

快速切除短路故障,除了能减轻电气设备因短路流产生的热效应等不良影响外,对于提高电力系统暂态稳定性,还有着决定性的意义。加快切除速度,可以减小切除角.这样既减小了加速面积,有增大了可能的减速面积,从而提高了暂态稳定性

由计算可得极限切除时间为t_c=0.4s.

当断路器切除时间设置为0.2s时:

图4-8 电机转速

图4-9 线电压

当设置的时间为0.5s时:

图4-10 电机转速

图4-11线电压

5 结论分析

当故障切除时间t_𝑐<t_(𝑐.𝑙𝑖𝑚) 时,系统保持暂态稳定。

当故障切除时间t_𝑐>t_(𝑐.𝑙𝑖𝑚) 时,系统不稳定,发电机失去同步,电压不稳定。

虽然系统有各种调节器,但快速切除短路故障对于提高电力系统的暂态稳定有着决定性的作用。

参考文献

[1]张晓春,郭小进,邹登海等.基于MATLAB语言的电力系统暂态稳定仿真[J]. 湖北电力,2000,03:9-10+14.

[2]李玉梅,马伟明,张波涛,等.PSB——基于SIMULINK的一种新型电力系统仿真工具[J],,电工技术杂志,2001,08:5-7.

[3]谢小荣,唐义良,崔文进. MATLAB在电力系统仿真中的应用[J].电工技术,2000,09:4-5.

[4]都伟杰,张俊芳,刘鹏,等. 基于MATLAB的电力系统暂态稳定性仿真研究[J]. 电网与清洁能源,2009,01:17-20.

[5]时宇琳,王宝华.基于Matlab的电力系统暂态稳定仿真实验与分析[J].实验室研究与探索,2010,04:40-43.

[6]华梁,史志平. 基于MATLAB的电力系统暂态稳定性仿真[J]. 硅谷,2010,11:63+55.

[7]糜作维,周遥,王林. 电力系统仿真工具综述[J].电气开关,2010,04:8-10+14.

[8]陶海英.基于MATLAB的电力系统稳态仿真分析[J].江西电力职业技术学院学报,2010,03:80-82.

[9]李广凯,李庚银.电力系统仿真软件综述[J].电气电子教学学报,2005,03:61-65.

[10]彭建飞,任岷,王树锦.MATLAB在电力系统仿真研究中的应用[J].计算机仿   真,2005,06:193-196.

[11]全宇文,李宝兴,李慧敏.电力系统暂态稳定分析方法的探讨[J].西北电力技术,2005,06:28-30+33.

[12]盛义发,唐耀庚,苏泽光,等.基于Matlab的电力系统故障的仿真分析[J].南华大学学报(理工版),2003,04:45-47+51.

[13]张颖,金维香.基于MATLAB的电力系统暂态稳定仿真研究[J].大电机技术,2004,01:46-48.

[14]樊艳芳,蔺红.MATLAB/SIMULINK在电力系统仿真中的应用[J].新疆大学学报(自然科学版),2004,02:205-207.

[15]盛义发,洪镇南.MATLAB在电力系统仿真中的应用[J].计算机仿真,2004,11:197-199.

[16]周兆庆,陈星莺.Matlab电力系统工具箱在电力系统机电暂态仿真中的应用[J]. 电力自动化设备,2005,04:51-55.

[17]蔡克红,李升.MATLAB PSB在电力系统暂态稳定中的仿真应用[J].江西电力,2006,03:4-7.

[18]曾江华,陈晓明,金伟,等.MATLAB在电力系统仿真中的运用[J].人民长江,2006,11:41-42.

[19]王辉,衡卫华.MATLAB在电力系统仿真中的应用和研究[J].安阳工学院学报,2007,01:61-66.

[20]陶海英,杜琴.基于MATLAB/SimPowerSystems的电力系统暂态仿真分析[J].江西电力职业技术学院学报,2011,03:46-48+60.

[21]孙联喜,牛文东.基于Matlab工具的电力系统故障分析[J].山东电力高等专科学校学报,2012,02:18-22.

[22]李晨,蒋德珑,程生安.电力系统暂态稳定分析方法的现状与发展[J].现代电子技术,2012,16:89-94.

[23]许允之,刘昊,冯宇,等.Matlab在电力系统仿真实验中的应用[J].实验技术与管理,2007,01:103-105+125.

  • 6
    点赞
  • 39
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 4
    评论
clear; clc; f=50; %系统额定频率 Tj=8.47; %归算后的发电机惯性时间常数 PT=1; %正常运行时发电机向无穷大系统传输的有功功率 P2M=0.48; %故障存在时发电机的最大功率 P3M=1.38; %故障切除后发电机的最大功率 %下面是利用改进欧拉法进行逐段计算 %需要注意故障切除前后电磁功率有跃变 h=0.05; %设置步长0.05s Duration=2; %设置计算时段长度2s Delta(1)=33.92; %初始功角 Omega(1)=1; %初始转速 t(1)=0; Delta_h=pi-asin(1/1.38); Delta_cm=acos((PT*(Delta_h-Delta(1)*pi/180)+P3M*cos(Delta_h)-P2M*cos(Delta(1)*pi/180))/(P3M-P2M))*180/pi; d=Delta(1); for i=1:round(Duration/h) if d<Delta_cm d_Delta(i)=(Omega(i)-1)*360*f; d_Omega(i)=(PT-P2M*sin(Delta(i)*pi/180))/Tj; Delta0(i+1)=Delta(i)+d_Delta(i)*h; Omega0(i+1)=Omega(i)+d_Omega(i)*h; d_Delta0(i+1)=(Omega0(i+1)-1)*f*360; d_Omega0(i+1)=(PT-P2M*sin(Delta0(i+1)*pi/180))/Tj; d_Deltaa(i+1)=(d_Delta(i)+d_Delta0(i+1))/2; d_Omegaa(i+1)=(d_Omega(i)+d_Omega0(i+1))/2; Delta(i+1)=Delta(i)+d_Deltaa(i+1)*h; Omega(i+1)=Omega(i)+d_Omegaa(i+1)*h; d=Delta(i+1); t(i+1)=i*h; T=t(i); end; end; s=sprintf(' 最大摇摆角 Delta_h=%f\n 最大切除角 Delta_cm=%f\n 最大稳定切除时间 Tmax=%f \n',Delta_h*180/pi,Delta_cm,T); disp(s); CutTime=input('输入故障切除时间\n'); %故障发生时的功角变化过程 Delta(1)=33.92; %初始功角 Omega(1)=1; %初始转速 t(1)=0; for i=1:round(CutTime/h) d_Delta(i)=(Omega(i)-1)*360*f; d_Omega(i)=(PT-P2M*sin(Delta(i)*pi/180))/Tj; Delta0(i+1)=Delta(i)+d_Delta(i)*h; Omega0(i+1)=Omega(i)+d_Omega(i)*h; d_Delta0(i+1)=(Omega0(i+1)-1)*f*360; d_Omega0(i+1)=(PT-P2M*sin(Delta0(i+1)*pi/180))/Tj; d_Deltaa(i+1)=(d_Delta(i)+d_Delta0(i+1))/2; d_Omegaa(i+1)=(d_Omega(i)+d_Omega0(i+1))/2; Delta(i+1)=Delta(i)+d_Deltaa(i+1)*h; Omega(i+1)=Omega(i)+d_Omegaa(i+1)*h; t(i+1)=i*h; end; %故障切除后的功角变化过程 for i=round(CutTime/h)+1:round(Duration/h) d_Delta(i)=(Omega(i)-1)*360*f; d_Omega(i)=(PT-P3M*sin(Delta(i)*pi/180))/Tj; Delta0(i+1)=Delta(i)+d_Delta(i)*h; Omega0(i+1)=Om

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

扬辰鑫

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值