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

摘要

本发明涉及一种包含随机功率单元相关性问题的随机潮流计算方法，包括S1；对系统包含的所有随机功率单元进行分组；S2；根据随机功率单元的分组情况，确定简化后的系统网络拓扑结构和蒙特卡洛法的仿真次数m；S3；获取步骤S1得到的各随机变量在仿真周期内的历史时序数据，统计得到其累积分布函数；S4；利用改进JNT法，得到各随机变量的采样值，并形成一个n维随机向量，重复m次；S5；利用蒙特卡洛法，分别求解电力系统潮流方程，获得经典潮流问题的解的集合。本发明充分利用了JNT法基于多维正态分布函数这一特性，提出一种改进的JNT采样法对随机功率单元间的相关性结构进行建模，可以使得随机潮流计算方法中的相关性建模工作更为简便。

主权项

一种包含随机功率单元相关性问题的随机潮流计算方法，其特征在于包括如下步骤：步骤S1；对系统包含的所有随机功率单元按照相关性和地理位置进行分组，每个随机功率单元组视作一个随机变量，记录所述随机变量的个数为n；步骤S2；根据随机功率单元的分组情况，确定简化后的系统网络拓扑结构和蒙特卡洛法的仿真次数m；步骤S3；获取步骤S1得到的各随机变量在仿真周期内的历史时序数据，统计得到其累积分布函数；步骤S4；利用改进联合正态变换法，即JNT法，得到各随机变量的采样值，并形成一个n维随机向量，重复m次；步骤S5；利用蒙特卡洛法，利用m个n维随机向量作为系统输入，分别求解电力系统潮流方程，获得经典潮流问题的解的集合；所述步骤S4具体包括以下步骤：步骤S41：计算n个随机变量x1至xn之间的等级相关系数矩阵Rr，所述的Rr为n维矩阵，并根据公式将Rr转化为积差相关系数矩阵R，其中rr为矩阵Rr中的任一元素，r是矩阵R中与rr处于对应位置的元素；步骤S42：由于所述步骤S41得到的R是正定矩阵，将R矩阵分解为R＝AAT的形式；步骤S43：对服从标准正态分布的任意一维随机变量进行n次独立的采样，构成各分量相互独立的总体采样向量并将所述通过进行正交变换得到所需的正态边际分布域中的n维随机向量步骤S44：利用正态分布的分布函数Φ将转换到均匀边际分布域：U→=Φ(N→)=(u1,u2,...,un)]]>Φ函数采用离散数据对或离散数据点来保存，两个保存的相邻离散点之间的数据点采用线性插值的方法来近似获取；步骤S45：利用x1至xn各自的累积分布函数的反函数对所述中的各个对应分量进行变换，得到实际边际分布域中的随机变量x1至xn的采样值。