一种用于大规模电力系统的振荡源定位方法技术方案

本发明专利技术公开了一种用于大规模电力系统的振荡源定位方法，属于电力系统稳定分析领域，其主要步骤包括：首先读取一段电力系统扰动后多台发电机的角速度和有功功率量测信号；然后对所述信号求取变化量；利用TLS‑ESPRIT方法对上述变化量进行分析，得到各信号在不同振荡模式下的低频振荡参数；再选择所关注的模式，在该模式下对有功功率、角速度进行重构；最后计算振荡能量，若某台发电机能量斜率为负，则可判定该机组振荡源。本发明专利技术方法利用PMU实测信息，能反应系统真实运行情况，同时，该方法能适用于大规模电力系统，具有较高的工程实用价值。

【技术实现步骤摘要】

本专利技术涉及电力系统稳定分析领域，具体涉及一种用于大规模电力系统的振荡源定位方法。
技术介绍
低频振荡是威胁电力系统安全运行的一个突出问题，快速并准确地定位出负阻尼机组对电网的低频振荡分析和控制有着极为重要的意义。目前，低频振荡产生机理主要有负阻尼机理和强迫振荡机理。负阻尼机理认为发电机阻尼转矩不足导致了小扰动功角失稳，原因在于高倍率的励磁放大系统产生的负阻尼抵消了系统固有的正阻尼，使得系统总阻尼下降甚至成为负值。负阻尼机理是目前发展最为完善的低频振荡机理，目前已成为解析低频振荡现象和设计电力系统稳定器(PowerSystemStabilizer，PSS)的主要理论依据。在负阻尼振荡中，负阻尼机组对模式阻尼削弱最严重，通常被定义为振荡源。目前，研究振荡源定位的方法主要有：混合仿真法、能量函数法和波形相似法。混合仿真法利用PMU单元的实时同步相量数据和分布模式信息将电网等效成若干子系统，通过对比仿真曲线和实测曲线来定位含扰动源的子系统，然后反复将含扰动源的子系统细分为更小子系统，并最终确定含扰动源的最小子系统，但是混合仿真法的解耦方法有待优化，以满足现场PMU的布点，并且混合仿真法不适用于在线定位。能量函数法应用较为广泛，但该方法的理论基于单机无穷大系统推导而来，应用在多机系统中会存在较大误差，因此对于实际大规模电力系统而言，该方法不够准确。波形相似法利用低频振荡时电压行波的特殊形态传播到电网中各节点所具有的相似性和时延性来实现对扰动源的定位，但是实际电网结构复杂，在电网参数时变情况下，波形相似法会产生较大误差。
技术实现思路
本专利技术所要解决的技术问题是提供一种用于大规模电力系统的振荡源定位方法，本方法便于量测、准确的高、适用于大规模电力系统振荡源定位，以及检测在不同振荡模式中各发电机是否提供了负阻尼。为解决上述技术问题，本专利技术采用的技术方案是：一种用于大规模电力系统的振荡源定位方法，包括以下步骤：步骤1：获取扰动后一段时间内各机组的量测量，所述量测量包括第i台发电机所连母线支路在k时刻的采样值，所述采样值包括有功功率Pi,k、角速度ωi,k，其中i＝1,2,…,s，s为系统中发电机个数，k＝1,2,…,L，L为设定时间段内的采样个数；步骤2：根据有功功率Pi,k、角速度ωi,k获得相应的有功功率变化量ΔPi、角速度变化量Δωi，即ΔPi＝Pi-Pi,0，Δωi＝ωi-ωi,0，Pi,0、ωi,0分别表示稳态有功功率、稳态角速度；步骤3：利用TLS-ESPRIT算法进行低频振荡模式辨识，分别辨识出有功功率变化量、无功功率变化量、电压幅值变化量以及电压相角变化量在不同振荡模式下的振荡幅值Amn、衰减因子σmn、振荡频率fmn和初始相位其中m＝1,2分别对应有功功率变化量、角速度变化量，n＝1,2,…,s-1，n表示振荡模式的个数；步骤4：选择弱阻尼振荡模式，并在该振荡模式下对有功功率变化量分量、角速度变化量分量进行重构，得到第n个振荡模式下各变化量分量ΔPi,n(t)、Δωi,n(t)，所述重构的公式为：步骤5：在第n个振荡模式下，计算振荡能量，判断有功功率变化量ΔPi,n和角速度变化量Δωi,n的相位关系，计算公式为：步骤6：根据第n个振荡模式下各发电机端口上的振荡能量判断振荡源所在机组，当振荡能量为负值时，表明有功功率变化量ΔPi,n和角速度变化量Δωi,n相位为90°～270°，该机组提供负阻尼，该机组为振荡源。与现有技术相比，本专利技术的有益效果是：a)不需要发电机内部参数，采用PMU实测数据定位振荡源，反应了系统的真实运行情况；b)适用于实际大规模电力系统，具有很高的工程实用价值。附图说明图1为本专利技术实施例中一种用于大规模电力系统的振荡源定位方法之流程图。图2为四机两区系统结构图。图3为故障后各发电机的有功功率振荡曲线。图4为发电机1有功功率在模式1下的重构曲线。图5为各发电机振荡能量计算结果。具体实施方式下面结合附图和具体实施方式对本专利技术作进一步详细的说明，本专利技术包括以下步骤：步骤1：获取扰动后一段时间内各机组的量测量，该量测量包括第i台发电机所连母线支路在k时刻的采样值，采样值包括有功功率Pi,k、角速度ωi,k，其中i＝1,2,…,s，s为系统中发电机个数，k＝1,2,…,L，L为设定时间段内的采样个数。步骤2：根据有功功率Pi,k、角速度ωi,k获得相应的有功功率变化量ΔPi，角速度变化量Δωi。步骤3：利用TLS-ESPRIT算法进行低频振荡模式辨识，分别辨识出有功功率变化量、无功功率变化量、电压幅值变化量以及电压相角变化量在不同振荡模式下的振荡幅值Amn、衰减因子σmn、振荡频率fmn和初始相位其中m＝1,2分别对应有功功率变化量、角速度变化量，n＝1,2,…,s-1，n表示振荡模式的个数。步骤4：选择弱阻尼模式，并在该振荡模式下对有功功率变化量分量、角速度变化量分量进行重构，得到第n个振荡模式下各变化量分量ΔPi,n(t)，Δωi,n(t)。步骤5：在第n个振荡模式下，按式计算振荡能量，判断有功功率变化量ΔPi,n和角速度变化量Δωi,n的相位关系。步骤6：根据上述在某一特定模式下各发电机端口上的振荡能量判断振荡源所在机组，当振荡能量为负值时，表明有功功率变化量ΔPi,n和角速度变化量Δωi,n相位为90°～270°，该机组提供负阻尼，该机组即为振荡源。步骤2中第n个振荡模式下获取变化量分量的公式为：步骤4中，对变化量进行重构的公式为：本专利技术基于模式分解的大规模电力系统振荡源定位方法的理论依据推导详述如下。多机系统的PH模型为：其中，Δδ，Δω，ΔEq′，ΔEfd分别为发电机的功角、转速、暂态电势、励磁电压列向量；M，T′d0，TA,KA,Dm分别为发电机惯性时间常数、励磁绕组时间常数、励磁调压器的时间常数、励磁调压器的放大倍数、机械阻尼系数对角阵；K1～K6为系数矩阵。令式(1)的系统矩阵为A，状态向量为X，式(1)经Laplace变换得(A-sI)X＝0(2)由式(2)可以推得ΔEq′＝(T′d0s+K3+(TAs+I)-1KAK6)-1·(-K4-(TAs+I)-1KAK5)Δδ(3)可知ΔTE＝K1Δδ+K2ΔEq′(4)将式(3)代入式(4)得ΔTE＝KE(s)Δδ(5)其中，KE(s)＝K1+K2(Td′0s+K3+(TAs+I)-1KAK6)-1·(-K4-(TAs+I)-1KAK5)其是n维方阵。由(1)可知，当不考虑阻尼系数D时考虑ω0＝1，得到若将第k台机的功角Δδk作为未知量(不失一般性，令k＝1)，上式变为由此得式中，R为一列向量。上式表明，系统中所有n-1台发电机的功角Δδi(i＝1,2,…,k-1,k+1,…,n)表示为剩下的第k台机的功角Δδk的线性函数形式，结合式(9)和式(6)，得到ΔTei＝Ai(s)Δδi(10)从式(10)看出，在多机系统中，第i台发电机的电磁转矩ΔTei是该发电机功角Δδi或角速度Δωi的函数，这种定义与单机无穷大系统中的电磁转矩系数在形式上一致。对于式(10)，令s＝jωd，得到ΔTei＝KsΔδi+KdjωdΔδi＝KsΔδi+KdΔωi(11)由线性系统理论知识可知，对于存在n个振荡模式的系统，小扰动后系统的零输本文档来自技高网...

【技术保护点】
一种用于大规模电力系统的振荡源定位方法，其特征在于，包括以下步骤：步骤1：获取扰动后一段时间内各机组的量测量，所述量测量包括第i台发电机所连母线支路在k时刻的采样值，所述采样值包括有功功率Pi,k、角速度ωi,k，其中i＝1,2,…,s，s为系统中发电机个数，k＝1,2,…,L，L为设定时间段内的采样个数；步骤2：根据有功功率Pi,k、角速度ωi,k获得相应的有功功率变化量ΔPi、角速度变化量Δωi，即ΔPi＝Pi‑Pi,0，Δωi＝ωi‑ωi,0，Pi,0、ωi,0分别表示稳态有功功率、稳态角速度；步骤3：利用TLS‑ESPRIT算法进行低频振荡模式辨识，分别辨识出有功功率变化量、无功功率变化量、电压幅值变化量以及电压相角变化量在不同振荡模式下的振荡幅值Amn、衰减因子σmn、振荡频率fmn和初始相位其中m＝1,2分别对应有功功率变化量、角速度变化量，n＝1,2,…,s‑1，n表示振荡模式的个数；步骤4：选择弱阻尼振荡模式，并在该振荡模式下对有功功率变化量分量、角速度变化量分量进行重构，得到第n个振荡模式下各变化量分量ΔPi,n(t)、Δωi,n(t)，所述重构的公式为：步骤5：在第n个振荡模式下，计算振荡能量，判断有功功率变化量ΔPi,n和角速度变化量Δωi,n的相位关系，计算公式为：Wi,n=∫t=t0t=t1ΔPi,ndΔωi,n]]>步骤6：根据第n个振荡模式下各发电机端口上的振荡能量判断振荡源所在机组，当振荡能量为负值时，表明有功功率变化量ΔPi,n和角速度变化量Δωi,n相位为90°～270°，该机组提供负阻尼，该机组为振荡源。...

【技术特征摘要】
1.一种用于大规模电力系统的振荡源定位方法，其特征在于，包括以下步骤：步骤1：获取扰动后一段时间内各机组的量测量，所述量测量包括第i台发电机所连母线支路在k时刻的采样值，所述采样值包括有功功率Pi,k、角速度ωi,k，其中i＝1,2,…,s，s为系统中发电机个数，k＝1,2,…,L，L为设定时间段内的采样个数；步骤2：根据有功功率Pi,k、角速度ωi,k获得相应的有功功率变化量ΔPi、角速度变化量Δωi，即ΔPi＝Pi-Pi,0，Δωi＝ωi-ωi,0，Pi,0、ωi,0分别表示稳态有功功率、稳态角速度；步骤3：利用TLS-ESPRIT算法进行低频振荡模式辨识，分别辨识出有功功率变化量、无功功率变化量、电压幅值变化量以及电压相角变化量在不同振荡模式下的振荡幅值Amn、衰减因子σmn、振荡频率fmn和初...

【专利技术属性】
技术研发人员：王晓茹，谢剑，
申请(专利权)人：西南交通大学，
类型：发明
国别省市：四川;51