基于相轨迹与相对动能变化率的电力系统暂态稳定性判别方法技术方案

本发明专利技术涉及基于相轨迹与相对动能变化率的电力系统暂态稳定性判别方法。本发明专利技术采集电力系统某一时刻中发电机组的信息；再判断电力系统是否发生故障，如果发生故障时把整个机群等值为单机；计算等值单机的相轨迹凹凸性指标并通过此指标初步判断系统稳定性；计算系统相对动能变化率，对其求导作为暂态失稳的辅助判据并依据此判据判断系统是否暂态失稳；如果不满足暂态失稳的辅助判据，则再通过判断系统是否到达不稳定平衡点来判断系统的暂态稳定性，如果没有到达不稳定平衡点则再采集下一时刻的发电机组信息循环上述步骤判断系统的暂态稳定性。本发明专利技术综合应用相轨迹凹凸性和相对动能变化率判据，可以准确、快速、实时判断多机系统的暂态稳定性。

【技术实现步骤摘要】
基于相轨迹与相对动能变化率的电力系统暂态稳定性判别方法
本专利技术涉及基于相轨迹与相对动能变化率的电力系统暂态稳定性判别方法，属于电力系统运行

技术介绍
现代电力系统发展为大规模的区域互联电网，这种发展可以带来巨大的经济性，同时也会使暂态稳定性问题更加复杂。暂态不稳定性仍然是现代电力系统面临的最大威胁之一。有效的实时暂态稳定性预测和应急控制至关重要。传统的暂态稳定控制策略主要采用“离线制定策略表，实时匹配”的控制方法。而系统模型和参数存在一些偏差，从而影响计算结果的精度。目前，常用的暂态稳定分析方法有时域仿真、基于李雅普诺夫能量函数的直接法、数值仿真法与直接法相结合的混合法、以扩展等面积法则为代表的系统等值法以及以人工智能算法为指导的暂稳分析方法等。随着广域测量系统(WAMS)的广泛应用，基于PMU/WAMS的实时电力系统暂态不稳定性识别已成为热点研究课题之一。出现了一系列利用相轨迹特征判断电力系统暂态稳定性多的方法。基于广域测量系统的轨迹几何特征的方法，没有网络动态等值，容易实现，只需要最新两个时间间隔的信息就可以判断暂态稳定性。但是判断标准是基于二维自制系统，对于复杂的非哈密顿系统会导致误判。因此，旨在提高相轨迹凹凸性方法的鲁棒性，以应对复杂的电力系统，从而引入暂态能量函数的方法，把相轨迹凹凸性与暂态不平衡能量变化率相结合，应用到暂稳识别中。
技术实现思路
本专利技术提供了基于相轨迹与相对动能变化率的电力系统暂态稳定性判别方法，以针对现存的利用相轨迹特征方法判断电力系统暂态稳定性存在精度低的问题。本专利技术的技术方案是：基于相轨迹与相对动能变化率的电力系统暂态稳定性判别方法，首先采集电力系统某一时刻中发电机组的信息；再判断电力系统是否发生故障，如果发生故障时把整个机群等值为单机；计算等值单机的相轨迹凹凸性指标并通过此指标初步判断系统稳定性；计算系统相对动能变化率，对其求导作为暂态失稳的辅助判据并依据此判据判断系统是否暂态失稳；如果不满足暂态失稳的辅助判据，则再通过判断系统是否到达不稳定平衡点来判断系统的暂态稳定性，如果没有到达不稳定平衡点则再采集下一时刻的发电机组信息循环上述步骤判断系统的暂态稳定性。所述判别方法的具体步骤如下：步骤S1、在t时刻，通过WAMS系统实时采集电力系统第i台发电机的功角δi、角速度Δωi、电磁功率Pei、机械功率Pmi，i＝1、2、3……n，n为电力系统中所有发电机的台数；步骤S2、根据步骤S1提取到的发电机信息，判断电力系统是否发生故障，如果没有发生故障则返回步骤S1，如果发生故障对所有发电机按功角从大到小进行排序，计算相邻两台机组的功角差，选出最大的功角差对应的两台发电机组作为分割线，分割线上面的机组为S群，下面的机组为A群，其中，最大的功角差对应的两台发电机组中功角大的机组分到S群，功角小的分到A群，并运用公式(1)和(2)进行两机群等值为单机：其中，δs、δa分别为S群和A群等值功角；δi、Mi、Pmi、Pei分别为第i台发电机的功角、转动惯量、机械功率及电磁功率；Ms和Ma分别为S群和A群等值转动惯量；Pms和Pma分别为S群和A群等值机械输入功率；Pes和Pea分别为S群和A群等值电气输出功率；步骤S3，计算等值单机的相轨迹凹凸性指标并通过此指标初步判断系统稳定性；当轨迹位于凹区域，能够判断系统暂态稳定，进入步骤S5；轨迹位于凸区域，系统有暂态失稳迹象，进入步骤S4；轨迹位于凹凸性分界线上，系统临界稳定，进入步骤S5；其中，Δω是对所有发电机等值为单机后的角速度；δ是对所有发电机等值为单机后的功角，δ＝δs‐δa；步骤S4，计算系统相对动能变化率，对其求导作为暂态失稳的辅助判据并依据此判据判断系统是否暂态失稳；当等值后的单机满足公式(3)时，判断系统暂态失稳，否则，不能判断系统失稳，进入步骤S5；VKE'＝(1/2)MΔωΔω'＞0(3)其中，VKE'是对所有发电机等值为单机后暂态动能一阶倒数，即系统暂态失稳的辅助判据；M是对所有发电机等值为单机后的转动惯量；步骤S5，再通过判断系统是否到达不稳定平衡点来判断系统的暂态稳定性；即判断所有发电机等值为单机后的功角δ是否大于δUEP；当满足公式(4)时，判断系统暂态失稳，假如不满足公式(4)则返回步骤S1，等待采集下一个时刻系统的信息；δ＞δUEP(4)其中，δUEP为单机等值后的系统不稳定平衡点。所述步骤S3，具体步骤如下：(1)建立单机无穷大系统的暂态过程；当线路发生故障时，在哈密顿单机无穷大系统中，忽略阻尼，不计调节器和调速器作用，描述发电机暂态过程的数学表达式如下：其中，δi为发电机功角；Δωi为发电机角速度；ω0为同步电角速度，电力系统中ω0＝100π；Mi为发电机转动惯量；Pmi为发电机机械功率；Pmax为发电机电磁功率的极值；(2)建立故障时角速度与功角的关系；对(5)式移相、积分，并且根据故障发生前后相轨迹的连续性可得：其中，Δω1是故障后的角速度；P1max是故障发生后电磁功率的极值；C1是根据故障前后相平面轨迹的连续性求出来的值；下标i＝1、2、3.......n，n代表电力系统所有发电机；(3)建立故障切除后角速度与功角的关系；故障切除时，Δω1，δi为已知，根据故障切除前后相平面轨迹的连续性，故障切除后相平面轨迹式为：其中，Δω2为故障切除后的角速度；P2max为故障切除后电磁功率的极值；C2是根据故障切除前后相平面轨迹的连续性求出来的值；(4)根据单机相轨迹是否到达初步判断系统的稳定性；根据(7)式，令当满足当若(7)式有解，则相轨迹在处取到极小值，Δω2将会单调增大，系统将失稳，若(7)式无解，说明相轨迹不会到达处，系统稳定；(5)根据单机相轨迹凹凸性初步判断系统稳定性；当相轨迹满足相轨迹在凹区域内，不会到达能够判断系统暂态稳定，进入步骤S5；当相轨迹进入凸区域，一定到达系统有暂态失稳迹象，进入步骤S4；相轨迹位于凹凸性分界线上，不能判断是否到达系统临界稳定，进入步骤S5。采用所述步骤S4的具体原因如下：(1)定义发电机i的暂态动能如下：VKEi＝(1/2)Mi(Δωi)2(8)(2)暂态动能的一阶导数如下：VKEi'＝MiΔωi'Δωi(9)(3)当相轨迹达到凸区域时，计算暂态动能二阶导数：(4)建立系统暂态失稳的辅助判据。当满足凹凸性指标μ＝Δω·l＞0，也就是-McosδP2maxω0Δω2＞(Pm-sinδP2max)2，此时暂态相对动能的二阶导数也恒大于零。根据不稳定平衡点的定义，暂态不稳定时，相轨迹将保持在凸区域，暂态动能的一阶导数将保持单调增长。因此，一定存在一个时间t0符合以下式子：VKEi'(t0)＝(1/2)MiΔωiΔωi'＞0(11)其中，VKEi是发电机的暂态动能；δi为发电机功角；Δωi为发电机角速度；Mi为发电机转动惯量。本专利技术的有益效果是：1)本专利技术可以实时在线判断电力系统任意摆暂态稳定性问题。2)本专利技术提出的相对动能变化率的辅助判据可以避免仅仅依靠相轨迹特征在系统临界稳定造成的误判。3)综合应用相轨迹凹凸性和相对动能变化率判据，可以准确判断所有机组的暂态稳定性，这对系统故障时准确切除不稳定机组具有重要实际意义。附图说明图1本专利技术的流程图；图2新本文档来自技高网...

【技术保护点】
基于相轨迹与相对动能变化率的电力系统暂态稳定性判别方法，其特征在于：首先采集电力系统某一时刻中发电机组的信息；再判断电力系统是否发生故障，如果发生故障时把整个机群等值为单机；计算等值单机的相轨迹凹凸性指标并通过此指标初步判断系统稳定性；计算系统相对动能变化率，对其求导作为暂态失稳的辅助判据并依据此判据判断系统是否暂态失稳；如果不满足暂态失稳的辅助判据，则再通过判断系统是否到达不稳定平衡点来判断系统的暂态稳定性，如果没有到达不稳定平衡点则再采集下一时刻的发电机组信息循环上述步骤判断系统的暂态稳定性。

【技术特征摘要】
1.基于相轨迹与相对动能变化率的电力系统暂态稳定性判别方法，其特征在于：首先采集电力系统某一时刻中发电机组的信息；再判断电力系统是否发生故障，如果发生故障时把整个机群等值为单机；计算等值单机的相轨迹凹凸性指标并通过此指标初步判断系统稳定性；计算系统相对动能变化率，对其求导作为暂态失稳的辅助判据并依据此判据判断系统是否暂态失稳；如果不满足暂态失稳的辅助判据，则再通过判断系统是否到达不稳定平衡点来判断系统的暂态稳定性，如果没有到达不稳定平衡点则再采集下一时刻的发电机组信息循环上述步骤判断系统的暂态稳定性。2.根据权利要求1所述的基于相轨迹与相对动能变化率的电力系统暂态稳定性判别方法，其特征在于：所述判别方法的具体步骤如下：步骤S1、在t时刻，通过WAMS系统实时采集电力系统第i台发电机的功角δi、角速度Δωi、电磁功率Pei、机械功率Pmi，i＝1、2、3……n，n为电力系统中所有发电机的台数；步骤S2、根据步骤S1提取到的发电机信息，判断电力系统是否发生故障，如果没有发生故障则返回步骤S1，如果发生故障对所有发电机按功角从大到小进行排序，计算相邻两台机组的功角差，选出最大的功角差对应的两台发电机组作为分割线，分割线上面的机组为S群，下面的机组为A群，其中，最大的功角差对应的两台发电机组中功角大的机组分到S群，功角小的分到A群，并运用公式(1)和(2)进行两机群等值为单机：其中，δs、δa分别为S群和A群等值功角；δi、Mi、Pmi、Pei分别为第i台发电机的功角、转动惯量、机械功率及电磁功率；Ms和Ma分别为S群和A群等值转动惯量；Pms和Pma分别为S群和A群等值机械输入功率；Pes和Pea分别为S群和A群等值电气输出功率；步骤S3，计算等值单机的相轨迹凹凸性指标并通过此指标初步判断系统稳定性；当轨迹位于凹区域，能够判断系统暂态稳定，进入步骤S5；轨迹位于凸区域，系统有暂态失稳迹象，进入步骤S4；轨迹位于凹凸性分界线上，系统临界稳定，进入步骤S5；其中，Δω是对所有发电机等值为单机后的角速度；δ是对所有发电机等值为单机后的功角，δ＝δs‐δa；步骤S4，计算系统相对动能变化率，对其求导作为暂态失稳的辅助...

【专利技术属性】
技术研发人员：刘志坚，孙应毕，宋琪，王洪亮，谢静，刘亚锦，罗琳玲，尹冬晖，海世芳，周洁，
申请(专利权)人：昆明理工大学，
类型：发明
国别省市：云南,53