一种电力系统的概率潮流计算方法与流程

技术特征：

1.一种电力系统的概率潮流计算方法，其特征在于，包括步骤：

采集电力系统中的有功功率变量，建立所述有功功率变量的样本矩阵；

根据所述样本矩阵构造三次样条函数，并通过所述三次样条函数获得所述有功功率变量的累积分布函数F(p_k)，k＝1,2,3，……，n；

将所述累积分布函数F(p_k)的曲线纵坐标划分为N个子区间，利用拉丁超立方抽样方法分别对所述N个子区间进行抽样，形成有功功率变量的采样矩阵；

利用Cholesky分解法处理所得采样矩阵，以获得相关性最小的采样样本，将所述采样样本作为输入量进行确定性潮流计算。

2.如权利要求1所述的概率潮流计算方法，其特征在于，所述根据样本矩阵构造三次样条函数，并通过所述三次样条函数获得所述有功功率变量的累积分布函数F(p_k)具体包括步骤：

设定区间[0,1]，将所述区间划分为N个分段，且满足0＝p₁<p₂<p₃…<p_n+1＝1，有功功率变量p_k为所述区间内的随机变量，k＝1,2,……n；

设定三次样条函数s(p)的方程：根据所述区间的N个分段，所述三次样条函数s(p)有4n个未知系数待求解，需要4n个方程求解，而根据三次样条函数的性质，可建立(3n+3)个方程，三次样条函数s(p)剩余(n-3)个方程未建立；

根据所述样本矩阵的前n-3阶矩阵构建出n-3个方程；

根据建立的4n个方程求解获得三次样条函数s(p)；

利用所述三次样条函数s(p)逼近有功功率变量的概率密度函数，并在所述区间上对逼近的三次样条函数s(p)求积分，以获取所述有功功率变量的累积分布函数F(pk)。

3.如权利要求2所述的概率潮流计算方法，其特征在于，所述将累积分布函数F(p_k)的曲线纵坐标划分为N个子区间，利用拉丁超立方抽样方法分别对所述N个子区间进行抽样，形成有功功率作为随机变量的采样矩阵，具体包括步骤：

将累积分布函数F(p_k)的曲线按照纵坐标在[0,1]范围内等分为N个子区间n＝1,2,……,N，并取每一个子区间的中点值r，

根据有功功率p_k在[0,1]范围的累积分布函数F(p_k)的逆函数计算出r对应有功功率的值并将p_kr作为子区间的样本点；

分别获取每一个有功功率p_k所在子区间的采样点，并根据所述采样点形成有功功率变量的采样矩阵，所述采样矩阵大小为n×N。

4.如权利要求3所述的概率潮流计算方法，其特征在于，所述利用Cholesky分解法处理包括：对有功功率变量的相关系数矩阵进行Cholesky分解得到下三角矩阵B。

5.如权利要求4所述的概率潮流计算方法，其特征在于，所述将采样矩阵作为输入量进行概率潮流计算之前，还包括生成顺序矩阵L_S的步骤：

假定向量Y由n个服从标准正态分布相互独立的随机变量y_i(i＝1，2，…，n)构成；

根据所述下三角矩阵B和向量Y得到矩阵Z，且，所述矩阵Z中随机变量Z₁，Z₂，…，Z_n服从标准正态分布；

将矩阵Z按行生成顺序矩阵L_S，所述顺序矩阵L_S为n×N矩阵，每一行的数值从1到N排序，且其排列顺序与Z所对应行的数值大小顺序相同。

6.如权利要求5所述的概率潮流计算方法，其特征在于，所述生成顺序矩阵L_S之后，还包括步骤：根据顺序矩阵L_S重新排列所述采样矩阵，以获取新采样矩阵，所述新采样矩阵用于作为输入量进行概率潮流计算。