电力系统暂态稳定分析的时域仿真法.docx
- 文档编号:9142580
- 上传时间:2023-05-17
- 格式:DOCX
- 页数:17
- 大小:190.84KB
电力系统暂态稳定分析的时域仿真法.docx
《电力系统暂态稳定分析的时域仿真法.docx》由会员分享,可在线阅读,更多相关《电力系统暂态稳定分析的时域仿真法.docx(17页珍藏版)》请在冰点文库上搜索。
电力系统暂态稳定分析的时域仿真法电力系统暂态稳定分析的时域仿真法引言暂态能量函数法是基于一个古典的力学概念发展而来的,该概念中指出:
“对于一个自由的(无外力作用的)动态系统,若系统的总能量V(V(X)0,X为系统状态量)随时间变化率恒为负,则系统总能量不断减少,直至最终达到一个最小值,即平衡状态,则此系统是稳定的”。
图9-1所示的滚球系统在无扰动时,球位于稳定平衡点(stableequilibriumpoint,SEP);受扰后,小球在扰动结束时位于高度h处(以SEP为参考点),并具有速度V。
该质量为m的小球,总能量V由动能jmv2及势能mgh(g为重力加速度)的和组成,即12Vmvmgh0若小球与壁有摩擦力,则受扰后能量在摩擦力作用下逐步减少;设小球所在容器的壁高为H(以SEP为参考点),当小球位于壁沿上,且速度为零时(即处于不稳定平衡状态),相应的势能为mgH,称此位置为不稳定平衡点(unstableequilibriumpoint,UEP),相应的势能为系统临界能量Vcr,即mgH根据运动原理,我们知道,若忽略容器壁摩擦,在扰动结束时小球总能量V大于临界能量乂时,则小球最终将滚出容器,而失去稳定性;反之V:
V“,则小球将在摩擦力作用下,能量逐步减少,最终静止于SEP。
而V二V“为临界状态,显然可根据(Vcr-V)判别稳定裕度。
对于一个实际系统要解决两个关键问题:
一是对于一个实际系统如何构造(定义)一个合理的暂态能量函数,它的大小应能正确地反映系统失去稳定的严重性;二是如何确定和系统临界稳定相对应的函数值,即临界能量,从而可通过对扰动结束时暂态能量函数值(即上例中的1mv2mgh)和临界值(即上例中的mgH)的比较来判别稳定性或确定稳定域。
这种判别稳定的方法统称为暂态能量函数法(transientenergyfunctionTEF法)。
它的特点是从能量的观点来判别稳定性,而不是根据系统运动的轨迹(物理量随时间变化的曲线)来判别稳定性,从而计算量少,速度快。
直接法的最大优点:
(1)能计及非线性,适应较大系统。
(2)计算速度快,不必逐步积分求摇摆曲线,而是通过能量判据来判别稳定。
(3)能给出稳定度。
直接法的两个较大缺点:
(1)模型较简单。
目前真正实用的软件采用发电机二阶经典模型,恒定阻抗负荷,尚不能计及励磁系统对稳定的作用。
(2)分析结果容易偏于保守。
这是因为李雅普诺夫直接法相应的稳定准则是充分条件,而不是必要条件。
此外在系统很大,或受到一系列扰动(如重合闸过程)时,直接法的速度、精度较差,故目前仅用于判别第一摇摆稳定性。
单机无穷大系统的直接法暂态稳定分析本节介绍在最简单的单机无穷大系统中如何构造暂态能量函数、确定系统临界能量,并进行暂态稳定分析。
对于图9-2中系统,若发电机采用经典二阶模型,忽略原动机及调速器动态,忽略励磁系统动态,则系统完整的标幺值数学模型为fd豹豹MPm-Pedtd(9-1)=codt式中,为转子角速度和同步速的偏差;:
为转子角;Pm=const.,为机械功率;Pe二EUsin:
为电磁功率,X、为发电机内电势E,及无穷XZ大系统电压U0间的系统总电抗(设电阻为零);其中E和U为常数;M为发电机惯性时间常数(即)。
设图9-2中系统在稳态时一.0,功角特性为Pe;在t=0时,线路上受到三相故障扰动,功角特性变为Pe
(2),此时发电机加速,转子角增加,直到-c时,切除故障线路,功角特性变为Pe(3)。
图9-2单机无穷大系统直接法分析要求研究的问题是:
如何用直接法判别故障切除后系统的第一摇摆稳定性。
显然对故障后的系统,稳定平衡点为S(对应s),不稳定平衡点为U(对应),在这二点上,均有电磁功率平衡,即Pe二Pm。
下面定义系统的暂态能量函数,设系统动能Vk为(注意,为与同步速之偏差,故稳态时vk=o)12Vk二2M(9-2)显然可以将式(9-1)的加速方程的二边对积分而求得故障切除时的动能,即VkcM;=;M牛d:
(Pm-Pe)d=加速面积A(9-3)若定义系统的势能Vp为以故障切除后系统稳定平衡点S为参考点的减速面积(反映系统吸收动能的性能),则故障切除时的系统势能为Vp厂J:
(P(3)-Pm)少二面积B(9-4)从而系统在扰动结束时总暂态能量V为Vc=VkcVp2m2:
(PJ一Pm)d面积(AB)(9-5)若将系统处于不稳定平衡点U(转子角:
u)时,系统以S点为参考点的势能作为临界能量Vcr(此值相当于滚球系统的Vcr=mgH),则Vcr二:
(Pe(3)-Pm)d二面积(BC)(9-6)和滚球系统相似,可作稳定判别如下。
当Vc:
Vcr,即图9-2中面积(A+B)v面积(B+C),或者说面积AV面积C时,则系统第一摆稳定;反之若VcVcr,则系统不稳定;Vc二Vcr时系统为临界状态。
显然这和等面积准则完全一致,是一个准确的稳定判据。
这里假定系统有足够的阻尼,若第一摆稳定,则以后作衰减振荡,趋于S点。
下面对此作一个简单讨论。
(1)从上面分析可知,最关键问题是如何构造(定义)一个反映系统稳定性的暂态能量函数,以及如何正确确定系统的临界能量,并以此作为稳定判别的标准。
(2)直接法分析稳定时,不必逐步积分求(t),而只要求出c和c,计算Vc,并设法确定Vcr,通过比较Vc与Vcr来判别稳定。
(3)对单机无穷大系统,UEP点不仅功率平衡(Pe二Pm),且系统巳二Pm来求解二及算Vcr,也可搜索Vp-max点,并取Vcr=Vp,max在多机系统中前一种途径又称UEP法。
后一种途径又称势能边界面法,即PEBS法(potentialenergyboundarysurface)。
在多机系统中也可把系统作单机一一无穷大等值,再用上述等面积准则判稳,又称EEAC法(extendedequalareacriteria。
(4)用本方法只能解决第一摇摆稳定问题。
(5)在分析中一般把转子阻尼忽略,使结果更保守些。
(6)可以用Ver-Vc作为系统稳定度的定量描述,从而对事故严重性排队,以便作动态安全分析,实际应用中使用的是规格化的稳定度V,通常定义文献建议2安全二12预警Vn=0.51警告二00.5严重警告0潜在危机(7)暂态能量函数同元件模型紧密相关,当采用复杂的元件模型和计及调节器动态时,相应的暂态能量函数将十分复杂,直接法暂态稳定分析过程也将复杂化。
对于图9-2系统,若在状态空间,即一相平面上作故障切除后系统的定常能量曲线族,即由126EUVC,)二VkVp二cM2(-Sin、-Pm)d2也12EU.M(cos-coss)-巳(-s)=C2x(9-8)作曲线,式中X、为故障切除后系统的视在电抗;C为参变量。
该曲线族见图9-3所示对系统能量取微分可知dV二Md(Pe-Pm)d=M-(Pm-P,3)ddt而由运动方程可知,在故障切除后系统运动轨迹上M2Pm-Pe(3),dt故其运动轨迹上dV=0,V=const.,即在故障切除后系统运动轨迹必为上述定常能量曲线族中一支。
当系统稳定时,由图9-2可知,发电机转子将围绕s点摇摆,可以证明在相平面上其轨迹为一围绕s点的封闭曲线。
设系统临界失稳,贝做障切除时,系统(,)到达临界轨迹,并沿临界轨迹运动而失稳;当系统转子角达;u时,、0(图9-3中T点)。
显然当故障切除时系统相应的forc)位于图9-3阴影域内任一点,系统均为稳定的;若位于此域外,贝S系统不稳定。
而临界轨迹所对应的系统总能量即为临界能量(Vc=Vcr),相应故障切除时间为临界切除时间。
在(j,0)这一点,系统动能为0=0),全部转化为势能,该势能与最大减速面积对应反映了临界能量。
图9-3相平面上V=const.曲线族以上介绍了单机无穷大系统的暂态能量函数的定义及临界能量的确定,以及系统暂态稳定判别准则(它和等面积准则完全一致),并在相平面上作了相应讨论。
扩展等面积法(EEAC)暂态稳定分析直接暂态稳定分析的RUEP法和PEBS法都是在多机系统条件下进行稳定分析的,而EEAC法则在单机无穷大等值条件下进行稳定分析,其特点是速度特别快,缺点是在一些特殊情况下,稳定分析的精度问题有时较突出。
下面对它作概要介绍。
元件模型假定,即发电机为经典二阶模型,负荷线性,网络线性,忽略原动机、调速系统及励磁系统动态。
设系统导纳阵收缩到只剩发电机内节点,则系统同步坐标下的数学模型如式(9-9)所示。
”WiMiPmiPei即:
即:
d(i=1,2,n)(9-9)-dti下面用EEAC法分析系统发生简单故障,在故障切除后,系统的暂态稳定性。
EEAC法分析假定系统失稳为双机模式,设系统主导UEP或失稳模式已知,把受扰严重的机群称为S,其余机群称为A,在同步坐标基础上定乂S及A机群的等值角度及速度为fFs=(送M浮J/MsHES佗弋M)/MsMS=送MiiH=S及AMjj)/MAIjA(9-10)+CMjJ/MajAMA八MjIP=A设S机群中各机组的转子角间无相对摆动,A机群类同,即或十SiS,0jS)匕=%+,0(jA)(9-12)在此假定上,对全系统作双机等值,并最终化为单机无穷大系统,用等面积准则判别稳定。
显然双机等值时,惯量中心S和A的运动方程为Mss八(Pmi-Pei)USMaAY(PmjPej)(9-13)-jA由式(9-12)假定,若再进一步简化,设rsrj,可知2巳iiS二EjGji*EjEkGikEi二EiBiisin(s-a)Giicos(s-a)IWA(9-14)同理Pej二Ej2Gjj*EjElGjllA,l韵nEjEkBjksin(A-S)GjkCOS(A-S)“A(9-15)式中,GjjBj二Yj为Y阵中元素。
进一步作单机对无穷大系统等值,取单机转子角为(无穷大系统转子角恒为零,作参考点)由式(9-13)可知1八(Pmi-Pei)(PmjPej)SiSMAjA(9-16)若定义单机惯性时间常数则式(9-16)改写为(9-17)式中式中EjEjGjj;D二二EjEjBjjjAiS从而系统的单机无穷大等值数学模型为MPm-Pe二Pm一FC一PmaxSin(-)(9-19)mAMSMtS;S一A;Mt二MaMMs,Ma,sA的定义见式(9-11);1Pm(MaPmiMsPmj)二COnst.;MTiSjAPc,Pmax,计算式见式(9-18)。
(a)正常运行时P-曲线;(b)第一摇摆稳定分析示意图所示。
根据式(9-19),可在功角平面上作出功角曲线如图9-11(a)Rd,若设故障前系统的功角特性为Peo(见图9-11(b),故障时为故障后为PeP,相应的功角特性表达式为(9-20)各时段Pc,Pmax及可由式(9-18)及相应节点导纳阵参数确定。
可推导得图9-11(b)中故障前稳定平衡点为0二siP叮巴。
0max0(9-21a)故障后稳定平衡点为1Pm_PCPp=SInpp(9-21b)pmaxP故障后不稳定平衡点转子角为-p2p。
从而图9-11(b)中加速面积Acc及最大减速面积Alec分别为Aacc=Pm-氐-PmaxDs|n(-d)d=(Pm_%)(-0)PmaxDCOS(-d)-C0S(o-D)7:
-Jp2pAdec=PCDPmaxPSin(-P)-PmdI=(Pcp-Pm)5-J-J+2Yp)+PmaxP【COS(J-Yp)+COSp-%)(9-22)式中,1为故障切除时等值转子角。
则可据式(9-22)用等面积准则判别第一摇摆稳定性,并可相似地定义稳定度Aacc为计算Adec及Aacc,都需要计算故障切除时单机无穷大等值系统的发电机转子角1,由于实际分析中已知的是故障切除时间t,因此要根据单机无穷大等值转子运动方程(9-19)计算1,其中的Pe与系统故障时功角特性即式(9-20)的第二式相对应。
贝肉于I5=co.(9-24)MPm-PcDPmaxDSIn(-d)(丿可将区间0,t分为若干时步,当t小时可取作一步,用高阶泰勒级数法快速计算t时刻对应的1。
也可在多机系统下,用高阶泰勒级数法计算t对应的发电机转速和转子角,再作单机无穷大等值得则更为准确。
设计算步长为t,而计算区间为tn,切1-J二tn-tn_i讥。
若记(m)表示:
对时间的m阶导数,并设tn时刻的n已知,则时rn2)en3)en4)的计算式由式(9-24)导出为与同步速偏差)。
(1)由梯形积分法则,当tn和时J=ndPm-FCd-PmaxDsin(、n-D)Pm-PCD-PmaxDsin(、n-D)2M
(1)At八育2Pm-2Pcd-PmaxDSinD)Si门(-D)2M(9-25a)当tn=0时,鳥1)=*n八0=0,不必计算。
(2)ndt=PmPCD-PmaxDSin(-D)/M(9-25b)tn耗(3)6ndt2tn1和PmaxDC。
心D)川(9-25c)当tn=0时,鳥4):
.
(1)n.3d0=0,贝ln3)=0,不必计算。
dt3tn普Cn口)(們2COS(-D)/0,则当tn=0时,.(4)nd3dt3PPmaxDcosCd)V2)tn/讥丄n-t2-n4)-=t423!
4!
(9-26a)=0,故上式可简化为根据式(9-25),可用四阶泰勒级数近似计算n1如下(也可导出更高阶泰勒级数计算si):
ndn当仁=0时,因J=打3)6=%+26(2t24d4t4(9-26b)若取,t=t,即一时步计算到故障切除,则(9-27)t2丄*240对式(9-27)引入修正因子和2,以便在t较大时修正用四阶泰勒级数计算引起的误差,则式(9-27)可改为J=%+订;602(2七)2+加绐打4(9-28)式中,:
i和2为经验性修正因子,通常取-i=0.30.6,另外有;=1。
1和2引入的机理可参阅有关文献。
当系统发生简单故障,故障切除时间不长时,用式(9-28)计算一是十分简捷的。
若要计算ter,则只要据式(9-22)求解1使Adec=Aacc,然后据式(9-28)求解.相应的t即为ter,可用迭代法或插值法求解。
在简单模型及简单故障条件下用EEAC法进行暂态稳定分析的步骤如下
(1)输入潮流计算结果,计算初始值
(2)形成系统故障前、故障时、故障切除后的节点导纳阵,分别追加负荷阻抗及发电机Xd支路并收缩到发电机内节点。
(3)据扰动计算t二0时各台机的加速功率和惯性时间常数之比Pacc,i/Mi,据此排队,决定可能的失稳模式。
(4)对每一种可能的失稳模式作如下计算:
(a)据式(9-18)及式(9-19)分别计算式(9-20)中的参数;(b)据式(9-21)计算故障前稳定平衡点及故障后稳定平衡点P和故障后不稳定平衡点-p2p);(c)据式(9-25)和式(9-26)用多步泰勒级数计算或据式(9-25)和式(9-27)或式(9-28)用单步泰勒级数计算故障切除时的等值转子角:
.,由式(9-22)和式(9-23)计算厶Vn进行稳定分析。
若计算tcr,可由Aacc=Adec求解cr再据高阶泰勒级数展开式由:
cr求解相应的G。
(5)对上述各种可能的失稳模式中以Vn最小值或tcr最小值相应的失稳模式为最终的失稳模式,相应的Vn及tcr为该故障下,系统的稳定裕度及临界切除时间。
大量工程实例计算表明,EEAC法的计算结果在大多数情况下的工程精度良好。
进一步研究主要是研究方法失效的判据,以鉴别出有大误差时的计算结果并加以改进;改进其对元件模型及扰动的适应性;扩展EEAC的应用领域,并使之实用化;以及改进失稳模式的判别方法和改进分析精度等。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电力系统 稳定 分析 时域 仿真