XXXX动态电力系统(第4章).pptx
- 文档编号:14732436
- 上传时间:2023-06-26
- 格式:PPTX
- 页数:70
- 大小:2.04MB
XXXX动态电力系统(第4章).pptx
《XXXX动态电力系统(第4章).pptx》由会员分享,可在线阅读,更多相关《XXXX动态电力系统(第4章).pptx(70页珍藏版)》请在冰点文库上搜索。
动态电力系统动态电力系统20122012年秋季年秋季研究生课研究生课程程第四章直接法暂态稳定分析主要内容直接法基本原理直接法基本原理能量函数构造能量函数构造RUEPRUEP法法PEBSPEBS法法EEACEEAC法法参考文献电力系统暂态稳定性的能量函数分析刘笙上海交通大学出版社1995直接法稳定分析付书裼、倪以信、薛禹胜中国电力出版社1999一、一、概论概论从能量的角度分析稳定问题从能量的角度分析稳定问题,从而快速判别稳定性从而快速判别稳定性。
优点:
优点:
可计及非线性大系统可计及非线性大系统计算速度快,不必逐步仿真受扰运动轨迹计算速度快,不必逐步仿真受扰运动轨迹能给出稳定裕度指标能给出稳定裕度指标缺点:
缺点:
模型较简单模型较简单结果偏于保守结果偏于保守稳定性分析历史回顾稳定性分析历史回顾1818世纪末,由于瓦特发明的离心式调速器有世纪末,由于瓦特发明的离心式调速器有时会造成系统的不稳定使蒸汽机产生剧烈的振时会造成系统的不稳定使蒸汽机产生剧烈的振荡。
到了荡。
到了1919世纪又发现船舶上自动操舵机的稳定世纪又发现船舶上自动操舵机的稳定性问题。
这就迫使一些数学家用微分方程来描述和性问题。
这就迫使一些数学家用微分方程来描述和分析系统的稳定性问题。
分析系统的稳定性问题。
1867年英国物理学家年英国物理学家J.C.麦克斯韦发表论调速器的文麦克斯韦发表论调速器的文章总结了无静差调速器的理论。
章总结了无静差调速器的理论。
1876年俄国机械学家维年俄国机械学家维什涅格拉茨基在法国科学院院报上发表论调节器的一般什涅格拉茨基在法国科学院院报上发表论调节器的一般理论的文章进一步总结了调节器的理论。
理论的文章进一步总结了调节器的理论。
1877年英国年英国数学家数学家E.J.劳思提出代数稳定判据即著名的劳思稳定判劳思提出代数稳定判据即著名的劳思稳定判据。
据。
1895年德国数学家年德国数学家A.胡尔维茨提出代数稳定判据的胡尔维茨提出代数稳定判据的另一种形式即著名的胡尔维茨稳定判据。
另一种形式即著名的胡尔维茨稳定判据。
Maxwell,J.C.“OnGovernors”.ProceedingsoftheRoyalSocietyofLondon,Vol.16(1867-1868):
270283.Lyapunov第一法通过求解微分方程的解来分析运动稳第一法通过求解微分方程的解来分析运动稳定性,即通过分析非线性系统线性化方程特征值分布定性,即通过分析非线性系统线性化方程特征值分布来判别原非线性系统的稳定性;来判别原非线性系统的稳定性;Lyapunov第二法是一种定性方法,它无需求解困难的第二法是一种定性方法,它无需求解困难的非线性微分方程,而转而构造一个非线性微分方程,而转而构造一个Lyapunov函数,函数,研究它的正定性及其对时间的沿系统方程解的全导数研究它的正定性及其对时间的沿系统方程解的全导数的负定或半负定,来得到稳定性的结论。
的负定或半负定,来得到稳定性的结论。
直接法发展历程直接法发展历程18921892年,年,LyapunovLyapunov经典论文经典论文关于运动稳定性的一般理论关于运动稳定性的一般理论发表发表李亚普诺夫李亚普诺夫(18571918)Lyapunov,AleksandrMikhailovich俄罗斯数学家,物理学家。
1857年生于雅罗斯拉夫尔,1918年11月3日卒于敖德萨。
1876年入圣彼得堡大学,1892年获博士学位,1893年起任哈尔科夫大学教授。
曾先后在圣彼得堡大学、哈尔科夫大学和喀山大学执教。
李亚普诺夫最初从事流体静力学理论研究,李亚普诺夫最初从事流体静力学理论研究,1892年年开创性地提出求解非线性常微分方程的李亚普诺夫函数法,开创性地提出求解非线性常微分方程的李亚普诺夫函数法,亦称直接法,由于这个方法的明显的几何直观和简明的分析亦称直接法,由于这个方法的明显的几何直观和简明的分析技巧,从而在科学技术的许多领域中得到广泛的应用和发展技巧,从而在科学技术的许多领域中得到广泛的应用和发展,奠定了常微分方程稳定性理论的基础。
成为研究常微分方,奠定了常微分方程稳定性理论的基础。
成为研究常微分方程定性理论的重要手段。
程定性理论的重要手段。
LyapunovLyapunov第二法第二法由力学经典理论可知,对于一个振动系统,当系统由力学经典理论可知,对于一个振动系统,当系统总能量(正定函数)连续减小(这意味着总能量对时间的总能量(正定函数)连续减小(这意味着总能量对时间的导数为负定),直到平衡状态时为止,则此振动系统是稳导数为负定),直到平衡状态时为止,则此振动系统是稳定的。
定的。
LyapunovLyapunov第二法是建立在更为普遍意义的基础上第二法是建立在更为普遍意义的基础上的,即如果系统有一个渐近稳定的平衡状态,则当其运动的,即如果系统有一个渐近稳定的平衡状态,则当其运动到平衡状态的吸引域内时,系统存储的能量随着时间的增到平衡状态的吸引域内时,系统存储的能量随着时间的增长而衰减,直到在平稳状态达到极小值为止。
然而对于一长而衰减,直到在平稳状态达到极小值为止。
然而对于一些纯数学系统,毕竟还没有一个定义“能量函数”的简便些纯数学系统,毕竟还没有一个定义“能量函数”的简便方法。
为了克服这个困难,方法。
为了克服这个困难,LyapunovLyapunov定义了一个虚构的能定义了一个虚构的能量函数,称为量函数,称为LyapunovLyapunov函数。
当然,这个函数无疑比能量函数。
当然,这个函数无疑比能量更为一般,且其应用也更广泛。
实际上,任一纯量函数只更为一般,且其应用也更广泛。
实际上,任一纯量函数只要满足要满足LyapunovLyapunov稳定性定理的假设条件,都可作为稳定性定理的假设条件,都可作为LyapunovLyapunov函数(其构造可能十分困难)。
函数(其构造可能十分困难)。
基本定义基本定义系统:
系统:
设(,)xfxtu稳定性是系统本身的一种动态属性,与外部输入无关。
0u,则(,)xfxt,()xt为n维向量,(,)fxt也是n维向量,12(,)iinxfxxxt,初始状态00()xtx。
平衡状态:
平衡状态:
定义定义:
系统(,)xfxt中对所有t,必存在一些状态点ex,使(,)0exfxt,该类状态点ex称为系统的平衡状态。
意义意义:
当系统运动到ex点时,系统状态各分量将维持平衡,不再随时间变化,即:
平衡点是由系统状态在状态空间中所确定的点。
求法求法:
(1)线性定常系统()()()0extAxtAxt若A非奇异,()0ext,唯一一个平衡点,坐标原点。
若A奇异,()ext有多个。
(2)非线性系统(,)0exfxt,ex,不只一个,可能有多个。
稳定性的定义稳定性的定义李亚普诺夫稳定性:
设(,)xfxt,若任意给定一个实数0,总存在另一个实数,使当0exx时,从任意初态0x出发的解00()(,)xttxt满足0(),exxtt,则称ex在李亚普诺夫意义下稳定。
几何意义几何意义:
从()S出发的轨迹,在0tt的任何时刻总不会超出()S。
)(S)(SeX若(,t0)与初始时刻t0无关,则称平衡态xe是李雅普诺夫意义下一致渐近稳定的。
对于线性定常系统来说,上述定义中的实数(,t0)可与初始时刻t0无关,故其渐近稳定性与一致渐近稳定性等价。
但对于时变系统来说,则这两者的意义很可能不同。
渐进稳定性ex在李亚普诺夫意义下稳定,且当t时,exx,lim0etxx。
几何意义几何意义:
从(S)出发的任何一个解,当t时,最终收敛于ex。
实际上是渐近稳定。
区别区别:
稳定和渐近稳定,两者有很大的不同。
工程上常常要求渐近稳定。
)(S)(SeX必要条件:
只有一个平衡点必要条件:
只有一个平衡点。
对于非线性系统,渐近稳定性是一个局部性的概念,而非全局性的概念。
大范围渐近稳定又称全局稳定。
是大范围渐近稳定。
则称都有若对任意eetxxtxx,)(lim0线性只要渐近稳定(只有一个ex),一定是整个状态空间的渐近稳定。
非线性系统,ex不只一个。
)(S)(SeX不稳定若当0exx时,总存在一个初态0x,使00()exxtt,称平衡状态ex是不稳定的。
基本思想基本思想虚构能量函数虚构能量函数()Vx李亚普诺夫函数李亚普诺夫函数既可以描述物理系统,又可描述社会系统,满足既可以描述物理系统,又可描述社会系统,满足3个个条件:
条件:
()Vx为任一标量函数,为任一标量函数,x为系统状态变量,是为系统状态变量,是t的函数。
的函数。
()Vx是正数是正数(正定的正定的)反映能量大小。
反映能量大小。
()()dVxVxdt,连续一阶偏导,反映能量变化速,连续一阶偏导,反映能量变化速度的大小,负值为能量减少。
度的大小,负值为能量减少。
李亚普诺夫直接法:
利用李亚普诺夫直接法:
利用()Vx和和()Vx的符号性质来直接的符号性质来直接判断系统在平衡处是否稳定。
判断系统在平衡处是否稳定。
李亚普诺夫判稳定理李亚普诺夫判稳定理定理定理1设系统设系统(,)xfxt。
平衡状态为。
平衡状态为(0,)0ft,如果存,如果存在一个具有连续的一阶偏导的标量函数在一个具有连续的一阶偏导的标量函数(,)Vxt,并且满,并且满足条件:
足条件:
(1)(,)Vxt是正定的;是正定的;
(2)(,)Vxt是负定的是负定的那么,在原点处的平衡状态是渐进稳定的。
那么,在原点处的平衡状态是渐进稳定的。
若若x,有,有(,)Vxt,当满足上述条件时,在原点处,当满足上述条件时,在原点处的平衡状态为大范围的平衡状态为大范围(全局全局)渐进稳定。
渐进稳定。
2x0x)(tx1x0)(002221xVCxx=+12221Cxx=+22221Cxx=+说明说明:
(1)物理意义:
物理意义:
()Vx代表物理系统存储的能量代表物理系统存储的能量()Vx为系统运动使能量变化的速度,为系统运动使能量变化的速度,()0Vx,能量,能量。
(2)几何意义:
几何意义:
设设2212()Vxxx,取,取221212(),VxxxCcc,越逼近圆越逼近圆心,半径心,半径c越小,代表能量越小。
越小,代表能量越小。
()0,Vxt,收敛于原点。
,收敛于原点。
所以,所以,()Vx代表代表x到原点的距离,到原点的距离,()Vx代表在代表在x点趋向原点趋向原点的速度。
点的速度。
说明说明:
(3)充分条件并不是必要条件,故若找不到这样的()Vx,并不一定是不稳定的。
(4)李亚普诺夫函数不唯一,只要满足2个条件,最简形式是二次型()TVxxAx。
定理定理2设系统(,)xfxt,平衡状态为(0,)0ft,如果存在一个具有一阶连续偏导数的标量函数(,)Vxt,且满足:
(1)(,)Vxt是正定的,
(2)()Vx是负半定的,(3)(,)Vxt在0x时不恒等于零则系统在原点处的平衡状态为大范围渐近稳定。
说明说明:
(1)(,)0Vxt,则(,)Vxtc,能量保持不变,不趋向原点,不是渐近稳定,是李亚普诺夫稳定。
(2)(,)Vxt只是在某个时刻暂时为0,其他的为负,能量衰减不会终止。
1x0x)(tx0ex例例系统系统12212xxxxx试确定系统在平衡点的稳定性。
试确定系统在平衡点的稳定性。
解解:
12221210000xxxxxxx,00ex设:
设:
2212()Vxxx
(1)()0Vx
(2)21222122122()2222()20Vxxxxxxxxxxx(3)()Vx在在0x时:
时:
分析:
分析:
10x,20x,0V20x,()0Vx,()Vx不恒为零,渐近稳定。
不恒为零,渐近稳定。
定理定理3对于(,)xfxt,(0,)0ft,且有连续偏导数的()Vx,满足条件:
(1)(,)Vxt是正定的;
(2)(,)Vxt是负半定的;则系统的平衡状态在李亚普诺夫意义下稳定。
1x2x0x0)(tx定理定理4对(,)xfxt,(0,)0ft,具有(,)Vxt满足:
(1)(,)Vxt是正定的
(2)(,)Vxt是正定的或(,)Vxt是半正定的)(3)(,)Vxt在0x时不恒为零则系统在原点处的平衡状态是不稳定的。
例例系统112212xxxxxx试确定系统在平衡状态的稳定性。
解解:
(1)由0x,121200xxxx得1200eexx
(2)选2212()0Vxxx则:
11222211221212()2222220Vxxxxxxxxxxxxx()()满足定理4,该系统为不稳定系统。
线性系统的稳定性与非线性系统的稳定性比较线性系统的稳定性与非线性系统的稳定性比较在线性定常系统中,若平衡状态是局部渐近稳定的在线性定常系统中,若平衡状态是局部渐近稳定的,则它是大范围渐近稳定的。
然而在非线性系统中,不是大,则它是大范围渐近稳定的。
然而在非线性系统中,不是大范围渐近稳定的平衡状态可能是局部渐近稳定的。
因此,线范围渐近稳定的平衡状态可能是局部渐近稳定的。
因此,线性定常系统平衡状态的渐近稳定性的含义和非线性系统的含性定常系统平衡状态的渐近稳定性的含义和非线性系统的含义完全不同。
义完全不同。
如果要具体检验一个实际非线性系统平衡状态的如果要具体检验一个实际非线性系统平衡状态的渐近稳定性,则仅用前述非线性系统的线性化模型之稳定性渐近稳定性,则仅用前述非线性系统的线性化模型之稳定性分析,即分析,即Lyapunov第一法是远远不够的,必须研究没有线第一法是远远不够的,必须研究没有线性化的非线性系统。
为此有如下几种基于性化的非线性系统。
为此有如下几种基于Lyapunov第二法第二法的方法可达成这一目的。
如克拉索夫斯基方法、的方法可达成这一目的。
如克拉索夫斯基方法、Schultz-Gibson变量梯度法、鲁里叶变量梯度法、鲁里叶(Lure)法,以及波波夫方法等法,以及波波夫方法等。
小结小结:
V(x)的构造方法是关键,但李亚普诺夫方法未给出构造V(x)的一般方法。
原理简单,实用困难。
重点理解:
()()VxVx代表着能量大小,总要求正定的。
正定不稳定代表能量变化的速度和方向。
负定渐近稳定在电力系统中的应用:
在电力系统中的应用:
2020年代初,等面积定则年代初,等面积定则19581958年,年,AylettAylett,能量积分判据,能量积分判据19661966年,年,GlessGless等,明确提出等,明确提出7070年代末,取得重大进展。
三大主导流派:
年代末,取得重大进展。
三大主导流派:
加速度法:
加速度法:
Ribbens-pavellaRUEP法:
法:
Kakimoto,Athyetc.相关不稳定平衡点法相关不稳定平衡点法RUEPRelatedUnstableEquilibriumPointPEBS法:
势能界面法法:
势能界面法PEBSPotentialEnergyBoundarySurface8080年代以来:
年代以来:
单机能量函数法:
单机能量函数法:
IMEFIMEF,19831983结构保留能量函数法:
结构保留能量函数法:
SPEF,1987EEAC法:
薛禹胜法:
薛禹胜,扩展等面积法,扩展等面积法,1988EEACExtendedEqualAreaCriteriaBCU法:
江晓东,法:
江晓东,1991BoundaryofStabilitybasedControllingUnstableEquilibriumPointMethod基于稳定域边界的主导不稳定基于稳定域边界的主导不稳定平衡点法平衡点法jF直接法的简单类比:
直接法的简单类比:
,对于一个滚球系统,当滚球系统无扰动时,滚球位于稳定平衡点SEP。
受扰后,滚球位于高度h处时,其速度为v。
设球的质量为m,则滚球具有的总能量为vhHUEP0mghmv21V2势能动能当滚球到达UEP点速度v=0,这时滚球的总能量为Vcr称为临界能量根据运动学原理,若mgHVcr处,稳定到滚球在摩擦力作用下回处,临界稳定滚球停在稳定性滚球将滚出容器,失去SEPVVUEPVVVVcrcrcr关键:
关键:
如何对一个实际系统构造一个合适的能量函数。
如何确定与系统临界稳定状态相对应的临界能量,从而通过比较判别稳定性。
VVVcr二、单机无穷大系统直接法暂态稳定分析二、单机无穷大系统直接法暂态稳定分析XLXT2XT1若发电机采用经典数学模型,忽略原动机和调速器动态,忽略励磁系统的动态,则系统的数学模型可写为:
式中,为转子角,为转子角速度与的偏差M为发电机惯性时间常数为机械输入功率为电磁功率。
ddtMddtPPmemPsinXVEePn受扰条件假设为:
在t=0时线路II发生某种短路故障,时切除故障线路。
故障前、故障期间及故障切除后的电磁功率与功角的关系如图PuSSoCPmPIPIIPIIIABCttcP,和PPsu0C故障切除瞬间的动能可表示为:
故障切除瞬间的动能可表示为:
故障切除瞬间,系统的势能定义为以故障切除后系统的稳定平衡点s为参考点到故障切除后的点c的减速面积,即:
AdPPddtdMMVccIImcck加速面积00)(212BdPPVcSmIIIcp减速面积)(从而,系统在扰动结束时的总暂态能量为:
将系统处于不稳定平衡点u时,系统以s点为参考点的势能作为临界能量,则:
若,即面积(A+B)面积(B+C)或面积A面积C则系统第一摇摆稳定等面积定则等面积定则)()(212BAdPPMVVVcSmIIIccpckc面积Vcr)()(CBdPPVuSmIIIcr面积crcVV三、三、多机系统直接法暂态稳定分析多机系统直接法暂态稳定分析目标预想事故条件下的临界切除时间目标预想事故条件下的临界切除时间tcrtcr步骤:
步骤:
构造故障后系统的构造故障后系统的LyapunovLyapunov函数或能量函数函数或能量函数V(X)V(X)对给定的故障,寻找对给定的故障,寻找V(X)V(X)的临界值的临界值VcrVcr对故障后系统的暂态方程进行仿真对故障后系统的暂态方程进行仿真,直至直至V(X)=Vcr,V(X)=Vcr,对应时刻为对应时刻为tcrtcr关键关键Vcr1.多机电力系统暂态能量函数多机电力系统暂态能量函数经典模型:
经典模型:
发电机:
发电机:
n+m个节点个节点n个节点个节点负荷:
恒定阻抗负荷:
恒定阻抗网络:
收缩至发电机内节点网络:
收缩至发电机内节点n,1,=i)(eimiiiiiPPdtdMdtdGredIEYIYV电磁功率可写为:
其中为发电机的内电势为仅保留发电机内电势节点后导纳短阵中的对角元和非对角元对应的自导纳和互导纳。
Pein,1,=i)cos()sin(,12*1*jinijijjiijiiijnjijieiieeiDCGEEYERIERPiiiEEijjiijBEECijjiijGEEDiiiiiijBGYijijijjBGY暂态稳定研究的是发电机相对运动,需确定参考暂态稳定研究的是发电机相对运动,需确定参考坐标的相位角。
坐标的相位角。
MARMAR坐标(坐标(machineanglereferencemachineanglereference)以第以第nn台发电机为参考机台发电机为参考机1-n,1,=i-)()()()(nenmnieimiMPPMPPnininidtddtdCOICOI坐标(坐标(centerofinertiacenterofinertia)以系统惯性中心为参考以系统惯性中心为参考定义:
这里,并满足COIiiniMM11COIiiniMM11niiMM1ddtCOICOICOI坐标下,定义新的发电机转子角度和角速度坐标下,定义新的发电机转子角度和角速度显然,显然,n,1,=iCOIiiCOIii011iniiiniiMMCOI坐标下发电机转子运动方程为:
其中推导利用了关系式:
n,1,=iCOIieimiiiiiPMMPPdtdMdtdnininijjiijijiijjnijiijniiiimieinimiCOIDPDCGEPPPP1111,1121)cos
(2)cos()sin()(CCijjiDDijjisin()sin()ijjicos()cos()ijji2iiimiiGEPP能量函数推导能量函数推导:
为简化起见,将由变量所组成的矢量表示为(),也表示为状态空间的点。
对于故障切除后系统的稳态平衡点和不稳定平衡点,可由下式中的右端项等于零来确定,即其中稳定平衡点用表示不稳定平衡点用表示1212,nnn,1,=i0)(0COIieimiiiPMMPPf)0,(s)0,(u,多机系统的暂态能量函数实际上是单机无穷大系统中能量函数的一种推广,即定义为:
21121iCOIinieimiiniiPKdPMMPPMVVVisi)()cos()cos()cos()(21111111121ninijjijiijninijsjsijiijnisiiiiniijisjsidDCPMV推导利用了关系式:
)()sin()sin(1111,1jininijjiijijininijijdCdCjisjsiisi)()cos()(cos1111,1jininijjiijijininijijdDdDjisjsiisi01niiCOIiisidPMM注记能量函数由三部分组成:
动能项势能项(转子位能与电磁储能)耗散项与积分路径有关,近似以线性路径计及或忽略函数正定性无法确定)()cos()cos()cos()(21111111121ninijjijiijninijsjsijiijnisiiiiniijisjsidDCPMV与单机系统类似,系统的临界能量将由不稳定平衡点()处的暂态势能决定,即:
多机电力系统在故障后,系统暂态稳定的判据为:
u,0)cos()cos()()0,(1111sjsiujuininijijsiuniiucrCPVViVVcccr(,)2.稳定域估计稳定域估计经典法:
又称经典法:
又称ClosestUEPClosestUEP法。
求出全部法。
求出全部UEPUEP,共有(共有(2n-1-1)个。
每一个失稳模式对)个。
每一个失稳模式对应一个临界势能,取最小者作为应一个临界势能,取最小者作为Vcr,十分,十分保守,耗时。
保守,耗时。
不稳定平衡点不稳定平衡点UEP)0,(u1nn2n1n1nC1nCC12n台机失稳二台机失稳一台机失稳)个(机系统失稳模式)0,(min(ucrVV经典法没有考虑故障地点对稳定的影响,经典法没有考虑故障地点对稳定的影响,因而十分保守因而十分保守,无法实用。
,无法实用。
研究表明,对于一特定故障,在所有研究表明,对于一特定故障,在所有种失稳模式中,必有一种是真正合理种失稳模式中,必有一种是真正合理的,即系统以这种模式趋于失稳。
由于这种模的,即系统以这种模式趋于失稳。
由于这种模式所对应的与故障地点、故障类型等有着紧密式所对应的与故障地点、故障类型等有着紧密的关系,因此,在以后的研究中提出了的关系,因此,在以后的研究中提出了RUEP法等方法为暂态能量函数法稳定分析向实用法法等方法为暂态能量函数法稳定分析向实用法迈进了一大步。
迈进了一大步。
()211n加速度法加速度法在t=0时刻求各发电机的加速度设第i台发电机加速度最大,选取近似临界能量仿真受扰动态方程至,重新计算加速度并排列。
仿真受扰动态方程至),(其中020,)0,(siuucrVVcrVV)0,(ucruVVcrcrtVVRUEPRUEP法法又称又称CUEP法。
法。
RUEPRelatedUnstableEquilibriumPointCUEPControllingUns
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- XXXX 动态 电力系统