[发明专利]一种暂态稳定裕度概率分布的计算方法

权利要求书

1.一种暂态稳定裕度概率分布的计算方法，其特征在于，所述方法包括如下步骤：

（1）确定初始状态下受扰系统的失稳模式和稳定裕度初始值η₀，并由用户确定采用的计算模式，包括基于纯解析法的暂态稳定裕度概率分布计算和综合时域仿真与解析法的暂态稳定裕度概率分布计算；

所述基于纯解析法的暂态稳定裕度概率分布计算包括如下步骤：

（2）计算所述受扰系统稳定裕度对临界群S驱动功率的灵敏度系数I；

（3）计算所述临界群S驱动功率的概率数字特征；

（4）计算所述受扰系统稳定裕度的各阶中心矩；

（5）计算所述受扰系统稳定裕度的累积分布函数；

所述综合时域仿真和解析法的暂态稳定裕度概率分布计算包括如下步骤：

（6）计算确定节点功率变化范围内的失稳模式变化情况；

（7）对所述各失稳模式计算稳定裕度对临界群驱动功率P_ms、临界群负荷功率P_ls、余下群驱动功率P_ma、余下群负荷功率P_la的灵敏度系数II；

（8）对所述各失稳模式计算该失稳模式下稳定裕度线性变化对应的各节点功率等于其概率分布均值时的稳定裕度初始值η′_0(i)；

（9）对所述失稳模式计算临界群驱动功率的概率数字特征、临界群负荷功率的概率数字特征、余下群驱动功率的概率数字特征和余下群负荷功率的概率数字特征；

（10）对所述各失稳模式计算稳定裕度的各阶中心矩；

（11）通过分段函数求取稳定裕度的累积分布函数。

2.如权利要求1所述的计算方法，其特征在于，所述概率数字特征即各阶半不变量γ_υ。

3.如权利要求1所述的计算方法，其特征在于，所述步骤（1）是令所述受扰系统中各节点功率等于其概率分布均值，对所述受扰系统中的设定故障，利用详细时域仿真和扩展等面积准则（EEAC），得到所述受扰系统的失稳模式和稳定裕度初始值η₀。

4.如权利要求1所述的计算方法，其特征在于，所述步骤（2）是令所述临界群S的驱动功率P_ms和余下群A的驱动功率P_ma的初始值分别等于各自群内机组驱动功率均值的总和，通过静态EEAC（SEEAC）算法求取所述受扰系统的稳定裕度对所述临界群S驱动功率P_ms的灵敏度系数I；所述灵敏度系数I包括一阶灵敏度系数和二阶灵敏度系数

5.如权利要求1所述的计算方法，其特征在于，所述步骤（3）是根据所述临界群S内每台发电机组驱动功率P_mk的概率分布，求取每台发电机驱动功率P_mk的各阶半不变量γ_υ(P_mk)，按所述临界群S中发电机组个数对所述各阶半不变量γ_υ(P_mk)求和；求取P_ms的各阶半不变量γ_υ(P_ms)，为

γυ(Pms)=Σk=1Nsγυ(Pmk)]]>

其中，k为机组编号，υ表示半不变量的阶数，Ns为临界群机组个数。