[发明专利]一种包含随机功率单元相关性问题的随机潮流计算方法

权利要求书

1.一种包含随机功率单元相关性问题的随机潮流计算方法，其特征在于包括如下步骤：

步骤S1；对系统包含的所有随机功率单元按照相关性和地理位置进行分组，每个随机功率单元组视作一个随机变量，记录所述随机变量的个数为n；

步骤S2；根据随机功率单元的分组情况，确定简化后的系统网络拓扑结构和蒙特卡洛法的仿真次数m；

步骤S3；获取步骤S1得到的各随机变量在仿真周期内的历史时序数据，统计得到其累积分布函数；

步骤S4；利用改进联合正态变换法，即JNT法，得到各随机变量的采样值，并形成一个n维随机向量，重复m次；

步骤S5；利用蒙特卡洛法，利用m个n维随机向量作为系统输入，分别求解电力系统潮流方程，获得经典潮流问题的解的集合；

所述步骤S4具体包括以下步骤：

步骤S41：计算n个随机变量x₁至x_n之间的等级相关系数矩阵R_r，所述的R_r为n维矩阵，并根据公式将R_r转化为积差相关系数矩阵R，

其中r_r为矩阵R_r中的任一元素，r是矩阵R中与r_r处于对应位置的元素；

步骤S42：由于所述步骤S41得到的R是正定矩阵，将R矩阵分解为R＝AA^T的形式；

步骤S43：对服从标准正态分布的任意一维随机变量进行n次独立的采样，构成各分量相互独立的总体采样向量并将所述通过进行正交变换得到所需的正态边际分布域中的n维随机向量

步骤S44：利用正态分布的分布函数Φ将转换到均匀边际分布域：

U→=Φ(N→)=(u1,u2,...,un)]]>

Φ函数采用离散数据对或离散数据点来保存，两个保存的相邻离散点之间的数据点采用线性插值的方法来近似获取；

步骤S45：利用x₁至x_n各自的累积分布函数的反函数对所述中的各个对应分量进行变换，得到实际边际分布域中的随机变量x₁至x_n的采样值。

2.根据权利要求1所述的一种包含随机功率单元相关性问题的随机潮流计算方法，其特征在于：所述步骤S1具体为：所述的随机功率单元分为随机处理单元和负荷单元，将地理位置接近并且具有强相关性的随机出力单元和负荷单元分别分组，将分好的每一个随机功率单元组视作一个随机变量，并假设各组内的功率单元均为完全正相关，根据分组情况确定对应的随机出力单元变量和负荷变量。

3.根据权利要求1所述的一种包含随机功率单元相关性问题的随机潮流计算方法，其特征在于：所述步骤S3包括如下具体步骤：

步骤S31：将各随机变量记为x₁至x_n，其历史时序数据按时间顺序编号，形成形式为(编号,功率数值)的数据对；

步骤S32：通过对步骤S31中所述数据对中的“编号”进行均匀分布采样，获取数据对的采样点；

步骤S33：将步骤S32中采样获得的数据对按照“功率数值”从小到大排列，统计不同的“功率数值”的采样值个数，并进行累加和归一化，用以形成各随机变量x₁至x_n的累积分布函数。