[发明专利]一种基于多项式逼近的概率型暂态稳定约束的最优潮流求解方法

权利要求书

1.一种基于多项式逼近的概率型暂态稳定约束的最优潮流求解方法，其特征在于，包括以下步骤：

S1.建立概率型暂态稳定约束的最优潮流的数学模型；

S11.概率型暂态稳定约束的最优潮流的一般数学模型；

S12.基于单机等值方法的对概率型暂态稳定约束的构造；

S2.基于步骤S1所建立的数学模型，基于多项式逼近方法求解概率型暂态稳定约束的的最优潮流；

S21.求解系统暂态变量的多项式逼近式；

S22.将概率型不等式约束转化为确定性不等式约束，得到便于求解的优化数学模型；

S23.基于原对偶内点法求解概率型暂态稳定约束的最优潮流的解。

2.根据权利要求1所述的一种基于多项式逼近的概率型暂态稳定约束的最优潮流求解方法，其特征在于，所述S11中的一般数学模型包括：

min E{c(x₀,y₀,u,p)}   (1.1)

概率型静态稳定约束：

0＝g(x₀,y₀,u,p)   (1.2)

0＝f(x₀,y₀,u,p)   (1.3)

P{H₀(x₀,y₀,u,p)≤H_0.max}≥β₀   (1.4)

暂态约束：

0＝g(x(t),y(t),u,p),t∈(t₀,t_end]   (1.5)

概率型暂态稳定约束：

P{H(x(t),y(t))≤H_max}≥β,t∈(t₀,t_end]    (1.7)

其中，x为与系统暂态动态有关的状态变量，y为包含节点电压及注入电流的代数变量，u为具有上下限的控制变量，p为随机参数，x(t),y(t)为系统状态变量及代数变量在暂态仿真时域t∈(t₀,t_end]内的值，x₀,y₀为它们在t＝t₀时刻的初始值，E{·}为期望运算符，c(x₀,y₀,u,p)为目标函数，g(x₀,y₀,u,p)为表示系统稳态数学模型的函数向量，f(x₀,y₀,u,p)为表示系统暂态数学模型的函数向量，P{·}为概率运算符，H₀(x₀,y₀,u,p)和H(x(t),y(t))分别为表示系统静态稳定约束和暂态稳定约束的函数向量，H_0.max和H_max为其上限值，β₀和β为相应概率稳定约束需满足的下限制。

3.根据权利要求1所述的一种基于多项式逼近的概率型暂态稳定约束的最优潮流求解方法，其特征在于，所述S12中的基于单机等值方法的对概率型暂态稳定约束(1.7)的构造包括：

-暂态稳定裕度为：

η＝A_dec-A_acc＝-M(ω(T_u))²/2   (1.12)

其中，η为暂态稳定裕度，A_dec为减速面积，A_acc为加速面积，M为机械惯量，ω为发电机的相对转子角，T_u为等值机组的电磁功率在暂态过程中刚好等于其机械功率的时刻；

-基于以上，概率型暂态稳定约束构造为：

P{η＞0}≥β   (1.13)。