电力系统概率潮流算法综述_刘宇

ＤＯＩ：１０．７５００／ＡＥＰＳ２０１３１０１４０１７

电力系统概率潮流算法综述

刘　宇１，

２，高　山１，２

，杨胜春３，姚建国３（１．东南大学电气工程学院，江苏省南京市２１００１８；

２．江苏省智能电网技术与装备重点实验室，东南大学，江苏省南京市２１００１８；

３．中国电力科学研究院（南京），江苏省南京市２１０００３

）摘要：概率潮流是解决电力系统不确定因素的重要基础。随着间歇性能源的发展与电力系统随机

性的提升，概率潮流在近些年来得到了广泛的研究。文中以算法的原理与优缺点为立足点，对电力

系统概率潮流算法研究进行综述。首先，

对概率潮流的研究问题进行阐述，简要介绍了概率潮流理论的发展、计算模型分类以及评价指标，并简述了概率潮流在电力系统中的应用情况。然后，按照

算法的不同原理将概率潮流算法进行分类，

基于不同类别的方法对实际应用的具体算法进行详细分析，分别介绍了不同算法的原理步骤以及优劣性和适用性，并针对各类方法进行了算法总体评价和发展趋势分析。最后，结合电力系统的最新发展要求对概率潮流算法的研究方向做出展望。关键词：不确定性；概率潮流；相关性；模拟采样；近似计算

收稿日期：２０１３－１０－１４；修回日期：２０１４－０７－

０１。国家高技术研究发展计划（８６３计划）资助项目（２０１１ＡＡ０５Ａ１０５）；国家电网公司科技项目（ＤＺ７１－１３－０３６）；北京市自然科学基金资助项目（３１３２０３５

）。０　引言

在传统电力系统分析中，负荷的波动、电网运行方式的变化和发电机的停运等因素造成了电力系统一定程度上的不确定性。随着电力工业的发展，以太阳能和风能等为代表的新能源接入电网，给电网

带来了明显的间歇性和随机性；

微网、分布式电源和电动汽车等配电网新概念的发展，大大增强了电源、负荷与电网之间的互动性［１］

，其直接结果导致了电力系统的不确定性显著增加，用于电力系统分析的

概率潮流算法的研究日益重要。

１９７４年，Ｂｏｒｋｏｗａｋａ提出概率潮流计算方法［２］

，用以解决电力系统中诸多不确定因素。在随后四十年的时间里，概率潮流理论与方法得到了发展。与其几乎同时出现的随机潮流［３］

和概率潮流相互补充融合，逐渐形成处理电力系统不确定因素的体系：一般认为对于电力系统短期不确定因素采用随机潮流处理，而对于长期的具备规律性的不确定因素采用概率潮流处理，后者更趋向于概率分布的计算。

概率潮流计算的提出与发展，其最显著的意义是在进行电力系统分析时，考虑了系统各种不确定因素的随机性，从而使得计算分析更加贴合实际电

网的运行状态。概率潮流的研究问题，主要集中在

３个层面：

系统模型、计算模型和计算方法。系统元件的不确定性是引入概率潮流的根本原

因，主要体现在发电机、负荷、输电线路和变压器的随机性。近些年，随着可再生能源并网规模的日益

提高和电力用户的市场行为日趋突出，

发电机和负荷的功率模型越发复杂。文献［４

］提出Ｋ均值聚类负荷模型，对研究时段内具有相近特征的系统负荷状态进行分析归类，构成多个等值负荷水平，实现了复杂负荷模型的快速计算。

就概率潮流计算模型而言，以四大类模型为主。

Ｂｏｒｋｏｗａｋａ基于简化的直流模型［２］

提出了概率潮流计算方法。为了提高潮流计算的精度，Ａｌｌａｎ分别在１９７６年和１９８１年提出了线性化交流模型［５］和分

段线性化交流模型［６］

，Ｓｏｋｉｅｒａｊｓｋｉ在１９７８年提出保留非线性的交流模型［

７］

。目前概率潮流的计算方法都是基于这４种模型进行。

对概率潮流算法的研究是概率潮流分析中的热点，具备广阔的研究空间与研究意义。一种性能良

好的概率潮流计算方法应满足以下指标［８］

：①能够

求出输出随机变量的数字特征（

包括均值和方差）及概率分布；②能够处理多个随机变量间的相关性；

③满足实用化要求，

结果具有足够精度的情况下尽量减少计算时间；④满足通用性要求，

对输入变量的数学模型不应有太高要求。这４项指标构成概率潮流算法研究的重点和难点，专家学者们从一方面或多方面入手展开研究，形成了当前的多种概率潮流

—

７２１—第３８卷　第２３期２０１４年１２月１０

日Ｖｏｌ．３８　Ｎｏ．２３

Ｄｅｃ．１０，２０１４

计算方法。就目前已有研究成果而言，概率潮流算

法可以大致分为模拟法、近似法和解析法３类。解析法结合电力系统复杂卷积计算的简化特性，通过

对卷积计算的特殊处理衍生而来，

广义上也可归入近似法。纵观现有概率潮流算法，存在的突出问题

是难以将计算准确度和时效性有机统一。

本文以模拟法、近似法和解析法为分类基准，综

述各类方法中具体算法的原理步骤、

优缺点评估以及发展趋势。以算法中现有的问题及不足为依托，

结合当前电力系统的研究热点与新兴技术，对概率潮流的发展方向进行展望。

１　概率潮流在电力系统中的应用

潮流计算是电力系统分析的基础，其典型计算

潮流方程为：

Ｗ＝ｆ（Ｘ，Ｙ）

Ｚ＝ｇ（Ｘ，Ｙ）

｛

（１

）式中：Ｗ为系统节点功率注入向量；Ｘ为节点电压向量；Ｙ为系统网络参数；Ｚ为系统支路潮流向量。

概率潮流的计算正是在上式的基础上，通过考虑输入变量Ｗ和Ｙ的概率特性，获得系统状态变量

Ｘ和Ｚ的分布情况，

从而全面地给出系统的运行状况和概率特征。

随着电力系统随机性的显著提升，概率潮流算法在电力系统分析中获得广泛应用，其主要应用方向可以分为以下几类。

１

）电力系统规划，包括电源规划、电网规划和无功规划等规划问题［

９］

。规划问题的求解均以潮流计算作为基础，在不确定性增加的电力系统中，概率潮流将成为含随机因素规划问题的求解前提。

２）静态安全分析，作为电力系统分析的基本问题，采用概率潮流的静态安全分析方法可以更加真

实地反映电力系统全面信息［

１０］

。３

）电力系统运行状态实时在线分析，包括机组组合、在线调度、电力市场机制下的源－网－荷互动［１１］

等。应用某些概率潮流算法，可以在不显著增加计算次数与时间的条件下，更为精确地分析电力系统的运行状态及变化趋势，为电力系统实时分析提供了强有力的工具。

４

）其他基于潮流运算的系统分析，例如最优化潮流计算［１２］、电力系统风险评估［１３

］、互动型配电网

潮流计算等。

２　模拟法概率潮流

模拟法概率潮流，是将电力系统中的不确定因

素作为随机变量建立概率模型，然后抽取概率分布的样本，最后统计输出变量的分布特征。传统的模拟法概率潮流计算方法一般是指随机采样的蒙特卡

洛模拟法［１４－

１６］，后来基于随机模拟法改进衍生出重

要抽样法［１７－１９］、拉丁超立方采样法［２０－２５］和拟蒙特卡洛方法［２６－

２８］等。

２．１　随机采样的蒙特卡洛模拟法

蒙特卡洛模拟是二战时期美国物理学家Ｍｅｔｒｏｐｏｌｉｓ在执行曼哈顿计划的过程中提出的。蒙特卡洛模拟法以随机模拟和统计实验为手段，是一种从随机变量的概率分布中，通过随机选择数字的方法产生一种符合该随机变量概率分布特性的随机数值序列，作为输入变量序列进行特定分析的求解方法。其计算关键与核心步骤如下：①对潮流方程的输入变量Ｗ构造相应的概率模型；②产生随机数序列，

作为系统的抽样输入进行大量的数字模拟，每一组采样值通过潮流计算得到相应的模拟实验值；③系统计算，

对模拟实验结果进行统计处理，给出所求问题的解。

蒙特卡洛模拟的优点在于样本数量足够大时，计算结果足够精确；并且计算量一般不受系统规模

的影响，

该方法的抽样次数与抽样精度的平方成反比。缺点在于为提高计算精度，往往需要提高系统

抽样规模，从而导致计算时长过大。考虑其精度优势，随机采样的蒙特卡洛模拟法一般用来作为基准方法进行比较，是衡量其他方法准确性的重要参考。２．２　重要抽样法

重要抽样法认为期望值附近的采样值对计算结果具有更大的影响力，因此可以重点关注期望值附近的点。基于此，重要抽样法的基本思路是保持原有样本期望值不变，通过改变已知变量概率分布来减小其方差，从而达到减少运算时间的目的。

重要抽样法的经典计算［１７］

公式为：Ｅ（Ｆ）＝∑Ｘ

Ｆ（Ｘ）Ｐ（Ｘ）＝∑Ｘ

Ｆ′（Ｘ）Ｐ′（Ｘ）Ｆ′（Ｘ）＝

Ｆ（Ｘ）Ｐ（Ｘ）

Ｐ′（Ｘ）烅

烄烆

（２

）式中：Ｐ（Ｘ）和Ｆ（Ｘ）分别为原分布中系统的概率函数和状态函数；Ｐ′（Ｘ）和Ｆ′（Ｘ）分别为新分布中系统的概率函数和状态函数。

如何选取新分布中系统的概率分布Ｐ′（Ｘ）使得随机变量在期望不变的情况下减小方差是重要抽

样法的关键步骤。文献［１８］采用迭代法搜索重要分布函数，给出了若干重要分布函数的定义方法，并结

合分散抽样的技巧提高重要抽样法的收敛速度。文献［１９

］利用蒙特卡洛方法模拟出负荷样本，然后利—

８２１—２０１４，３８（２３

）　

用核密度估计方法估计出负荷模型的密度函数，将

之作为重要抽样密度函数，计算出支路潮流和节点电压的概率密度函数。

重要抽样法在电力系统的概率估计中有着广泛应用，该方法可以快速准确地计算出系统运行状态的期望值，

为系统分析提供参考。但重要抽样法中仅以期望为研究对象，对于概率变量的方差、概率分布等参数分析存在天然缺陷，

计算结果局限性较大。２．３　拉丁超立方采样法

为了避免随机采样的蒙特卡洛模拟法的大规模抽样，Ｍｃｋａｙ等人于１

９７９年提出了拉丁超立方采样法［２０］

。它是一种分层采样法，通过改进输入随机

变量的样本生成过程，

保证其采样值能够有效地反映随机变量的整体分布，算法的出发点就是确保所

有的采样区域都能够被采样点覆盖。其基本运算过

程分为如下两个步骤［

２１］

。１）采样。假设Ｘ１，Ｘ２，…，ＸＫ是一个待求概

率问题中的Ｋ个输入随机变量，Ｘｋ为Ｘ１，Ｘ２，

…，ＸＫ中任意一个随机变量，

其累计概率分布函数为：Ｙｋ＝Ｆｋ（Ｘｋ）（３

） 设Ｎ表示采样规模，

拉丁超立方采样示意图如图１所示

。

图１　拉丁超立方采样示意图

Ｆｉｇ．１　Ｓａｍｐｌｉｎｇ　ｓｔｅｐ　

ｉｎ　Ｌａｔｉｎ　ｈｙｐｅｒｃｕｂｅ　ｓａｍｐｌｉｎｇ其具体步骤为：将图１纵轴分为Ｎ个等间距不

重叠的区间，选择每个区间的中点作为Ｙｋ的采样

值，然后用式（４

）所示的反函数计算Ｘｋ的第ｎ个采样值。

ｘｋｎ＝Ｆ

－

１ｋ

（ｎ－０．５

Ｎ

）

（４） 式（

４）中，可以选取任意［０，１］区间中的随机数代替区间中点０．５，

从而实现采样的不同特性。Ｋ个输入随机变量的Ｎ次采样结果构成规模为Ｋ×Ｎ阶的初始采样矩阵Ｘｓ。

２

）排列。改变各随机变量采样值的排列顺序，使相互独立的随机变量采样值的相关性趋于最小。在排列过程中先形成一个Ｋ×Ｎ阶的顺序矩阵Ｌｓ，

该矩阵中的每一行元素值代表初始采样矩阵Ｘｓ对

应行的元素应该排列的位置，Ｘｓ的元素按照顺序矩

阵Ｌｓ排列后得到最终的采样矩阵。文献［２２］

列举了计算顺序矩阵Ｌｓ的方法，核心思想是利用如序列正交化等方法最大程度降低相关性。通过排列步骤中顺序矩阵的优化计算问题，又衍生出诸多的改进

拉丁超立方采样算法［２

３－

２５］，提高了拉丁超立方采样法的应用范围和计算精度。

拉丁超立方采样法的不足是对输入随机变量的

处理较为复杂，一方面要求已知输入随机变量的概

率分布函数或累积分布函数［

２５］

，另一方面对不同类型概率分布的随机变量相关性需要特殊变换处

理［

８］

。但该方法作为一种非常有效的估计输出随机变量期望值的方法，由于采样值能够确保覆盖所有输入随机变量的整个分布区域，无须大规模抽样，并且可以有效处理输入变量之间的相关性和随机性，在准确性、稳健性和时效性上都有较大的优势。２．４　拟蒙特卡洛法

拟蒙特卡洛法（ＱＭＣＳ

）的出发点与拉丁超立方采样法相同，希望通过有效的空间覆盖采样法来规避蒙特卡洛模拟法中的随机抽样。但与拉丁超立方采样法的处理方式不同，ＱＭＣＳ采用低差异序列实

现多维随机变量的空间采样［

２６］

。低差异序列，又称伪随机数列，是一系列数值确定的［０，１］

区间中的数。在ｄ维变量的空间中，低差异序列中已有ｎ－１个数，

生成第ｎ个数的方法是：将这个数插入已有数列中最大的“空白”处，即避

免数列在局部空间聚集，从而保证了有限数据的空间全覆盖。Ｖａｎ　ｄｅｒ　Ｃｏｒｐ

ｕｔ序列是典型的一维空间上的低差异序列，

以２为基数前８位数字的生成情况如表１所示［

２７］。表１　Ｖａｎ　ｄｅｒ　Ｃｏｒｐ

ｕｔ序列生成过程Ｔａｂｌｅ　１　Ｂｕｉｌｄｉｎｇ　

ｐｒｏｃｅｄｕｒｅ　ｏｆ　ａ　Ｖａｎ　ｄｅｒ　Ｃｏｒｐｕｔ　ｓｅｑｕｅｎｃｅ采样点２进制表示２进制反转序列值００００．００．００００．０００１　００１．０　０．１００　０．５００２　０１０．０　０．０１０　０．２５０３　０１１．０　０．１１０　０．７５０４　１００．０　０．００１　０．１２５５　１０１．０　０．１０１　０．６２５６　１１０．０　０．０１１　０．３７５７　

１１１．０　

０．１１１　

０．８７５

基于上述基本序列，为了解决实际应用中高维度变量的问题，需要进行维度的扩展。根据低差异序列的生成方法，随着维度的上升，在多维空间中寻找“空白”的难度就越大，因此，在高维度上的数据覆盖性会显著降低，这是目前已有低差异序列的主要

—

９２１—·综述·　刘　宇，等　电力系统概率潮流算法综述

不足。文献［２８

］中讨论了常用的３种多维度低差异序列，即Ｈａｌｔｏｎ序列、Ｆａｕｒｅ序列和Ｓｏｂｏｌ序列的优缺点，其中：Ｈａｌｔｏｎ序列变量维度最高６～８维；

Ｆａｕｒｅ序列可达２５维，但其计算时间急剧增加；Ｓｏｂｏｌ序列优越性较为明显，

其变量维度可达百位。目前ＱＭＣＳ已经被应用于概率最优潮流计

算［２６］和含互动式新能源的电网静态稳定分析［２８］

中。由于采样过程中一次性生成所需序列，该方法

具有比拉丁超立方采样法更高的计算效率。但

ＱＭＣＳ对多变量的高维度问题理论基础薄弱、计算效果差，

因此目前多用于小规模电力系统分析计算。２．５　模拟法概率潮流展望

模拟法的基础是蒙特卡洛模拟，因此该方法是目前电力系统中处理含不确定因素概率潮流问题中计算结果最为精确的一类方法，但与其高精度计算能力相对应的是低时效性。在蒙特卡洛模拟基础上改进的各种算法都是以减少抽样计算量为目的：重要抽样法计算速度快，期望值结果准确，但其计算结果单一，无法全面描述系统运行特征；拉丁超立方采样法的采样规模相对于蒙特卡洛模拟法大幅降低，是目前精度和时效较高的概率潮流算法，但该方法在处理变量之间的相关性时具备较高的复杂度，且相关性结果往往与实际情况存在一定的偏差；ＱＭＣＳ理论上计算速度最快，

但低差异序列的维数问题是该方法应用瓶颈。

随着电子计算机技术的发展与电力系统计算平台的壮大，模拟法必然在概率潮流计算中扮演越来越重要的角色，本文认为其发展趋势可以包含以下几个方向。

１）采用多种采样方式相结合的混合采样法，例如ＱＭＣＳ对维度较低的问题具备优越性，当问题维度较高时，采用ＱＭＣＳ与拉丁超立方采样法的混合采样法计算，从而发挥不同采样方法的优势。

２）采样过程中针对不同的研究重点分段处理，期望值附近的概率曲线在处理过程中保持良好的计算准确度，但分布曲线的尾部计算精度往往很差，在对精度要求很高时，可以对分布曲线分段，对尾部采样进行特殊处理。

３

）并行计算的应用。目前计算机的多核并行处理计算能力得到了大幅提升，同时基于图像处理单元的并行计算得到了发展，该硬件平台对重复性抽样计算的提速效果显著。

３　近似法概率潮流

近似法是利用输入随机变量的数字特征近似描

述系统状态变量统计特性的方法。该方法避开了大规模的重复抽样，因而求解速度较快，又因其能够计

及系统输入变量之间的互相关性，因而受到重视。

目前研究应用较多的有点估计法［２９－

３５］、一次二阶矩

法［

３６－３８］和状态变换法［２３，３９－

４０］。３．１　点估计法

点估计法是一种概率统计方法，目前所做的应用研究都是基于１９９８年Ｈｏｎｇ在已知输入随机变

量的连续分布下提出的点估计法［

２９］

。该方法能够根据已知随机变量的概率分布，求得待求随机变量的各阶矩。

点估计法属于逼近技术的一种，利用输入随机变量的统计信息来逼近输出随机变量的数字特

征［

３０－

３１］。其主要运算过程分为以下几步。１

）用潮流方程中输入随机变量Ｗ的各个分布函数求出相应的前２　

Ｍ－１阶中心矩。２）通过构造的方式，利用前２　Ｍ－１阶中心矩独立求出每个输入随机变量的Ｍ个离散状态，使得这

Ｍ个离散状态包含了前２　

Ｍ－１阶中心矩的所有信息。

３

）用所求得的每个输入随机变量的Ｍ个离散状态和它们的均值，构造Ｍ×Ｋ个输入随机变量的离散状态，求出对应输出随机变量的Ｍ×Ｋ个离散状态。

４）用求得的潮流方程输出随机变量Ｘ和Ｚ的Ｍ×Ｋ个离散状态逼近相应的期望值与方差等相关数字特征。

由以上步骤可以分析点估计法的特点如下。１

）该方法中实际的输入量为输入随机变量前２　

Ｍ－１阶中心矩，此中心矩可以由概率分布函数直接求出，也可以由大量样本逼近拟合方程式展开得到［３２

］，这样就不必受限于必须已知输入变量概率分布的条件约束。

２

）点估计法不需要知道输入与输出之间的具体函数关系表达式，仅要求每个输入有唯一对应的输出。

３）输出随机变量有２　Ｍ－１阶多项式逼近的精度，为了提高估计的精度，可以增加输入变量的高阶矩信息，即增加取点个数。但实际应用中点个数Ｍ大于３时不仅急剧增大计算量，而且往往造成解的

结果非实数，因此Ｍ通常取２或３［３３］

，即构成常用

的两点估计法和三点估计法。

两点估计法计算简单、容易实现，但其只利用输入变量的前三阶矩信息，计算精度低；三点估计法既能得到较高精度的估计值，又保持了简易性，在点估计法中广为使用。点估计法的缺点在于计算结果中随机变量的高阶矩不够精确，无法准确获得变量的

—

０３１—２０１４，３８（２３

）　

概率分布函数［３４］；同时在处理输入变量的时间和空间相关性上具有一定的计算复杂度［３５］。

３．２　一次二阶矩法

一次二阶矩法作为一种近似概率仿真方法，已被广泛应用于机械、结构可靠性分析中。该方法通过将状态方程泰勒展开，近似保留一次线性项，形成包含前两阶矩（即均值和方差）的计算方程式［３６－３７］。在电力系统概率潮流分析中，其具体步骤如下［３８］。

步骤１：将输入随机变量对输出状态变量的潮流方程按泰勒级数展开为一次项形式。

步骤２：计算输入变量均值方程式。

步骤３：计算输入变量协方差方程式。

步骤４：由步骤２和３联合，通过输入变量的均值和协方差计算输出状态变量的数字特征。

一次二阶矩法计算简单，效率高；但其计算能力有限，仅能处理输出与输入之间均值和方差的数值计算，算法模型误差较大，并且计算精度受到系统概率潮流模型约束很大，因此研究较少。

３．３　状态变换法

状态变换按照变换方法分为线性变换、多项式变换和无迹变换等。线性变换法基于正态变量线性变换不变性定理，假设节点注入随机变量均服从正态分布，将潮流方程线性化后可得系统状态变量为节点注入变量的线性组合并且仍服从正态分布［３９］。多项式变换多用作其他计算方法的辅助手段，用以表征电力系统随机因素的模型转换等问题［２３］。无迹变换认为：拟合一个概率分布比求解非线性变换容易得多［４０］，基于此，通过较少的样本点和相应的样本权重准确捕获状态分布参数，通过非线性函数传递后输出状态变量的期望与方差［４１］。

状态变换法的优点在于变换过程数学理论清晰，意义明确，计算规模不大，但由于该方法以高斯正态分布为变换基础，使其在新能源（不满足高斯分布）并网问题下的概率分析存在一些不足。

３．４　近似法概率潮流展望

近似在各个科学领域均有应用，在电力系统概率潮流中近似计算的使用也较为普遍。与确定性潮流相比，概率潮流的计算规模显著提升，快速的近似计算显得更为迫切。目前已有的近似法概率潮流算法中，计算结果以均值和方差为主要目标，如何获得状态变量较为准确的整体概率分布是改进的重要方向。另一方面，近似法概率潮流算法计算结果的可信度与误差分析也是研究的内容之一。

在近似法中，随机变量状态的变换是一种基础计算手段，而泛函分析领域的空间转换与之具有相似的计算逻辑，把发展成熟的泛函空间变换应用于

近似法概率潮流，是值得探索的方向。

４　解析法概率潮流

解析法概率潮流主要关注如何利用随机变量间的关系进行卷积计算得到状态量的概率分布，即核心思想在于有效处理复杂的卷积计算。在卷积计算的过程中，往往会做一些近似处理，因此解析法广义上可归为近似法。传统的卷积技术是一种获得概率潮流的基本方法［２，５，４２］，假设变量独立，在已知注入概率分布的情况下，可以通过卷积技术得到待求状态变量的概率密度函数。但对于多元线性方程的卷积计算量十分庞大，因此限制了其使用。已有解析法中卷积计算多采用快速傅里叶变换［４３］、半不变量法［４４－５１］和序列运算理论［５２－５４］。

解析法处理卷积计算的应用前提是假设输入随机变量相互独立，导致解析法在处理变量相关性上具有固有缺陷，因此，解析方法中变量相关性处理是研究的热点［４８－４９，５１］。

４．１　快速傅里叶变换

在信号处理学科领域，快速傅里叶变换是处理卷积问题的最佳方式［４３］，其具备良好的精度和效率，对于处理小规模数据系统具备很大的优越性。在电力系统方面，随着系统规模的增大，系统分析的输入变量急剧增加，使得快速傅里叶变换不再具备精度和效率上的优势，因此，在２０世纪８０年代经历了短暂的研究后，快速傅里叶变换逐渐退出电力系统概率潮流分析领域。

４．２　半不变量法

目前在电力系统广泛使用的处理卷积的方法是半不变量法。该方法的核心思想是将复杂的卷积运算转换为半不变量之间简单的算术运算，从而大大降低计算过程的复杂度。具体计算步骤可以概括如下［４４－４７］：将潮流方程中的随机变量Ｗ进行概率分布拟合，经过中心矩计算出输入随机变量的各阶半不变量；然后，通过雅可比矩阵和灵敏度矩阵进行简单的数学运算获得潮流方程输出变量Ｘ和Ｚ的各阶半不变量，结合不同级数扩展方式得到相应输出变量的概率密度扩展方程，从而获得输出随机变量的概率分布。文献［４８］深入分析了半不变量法计算随机潮流时各环节的假设条件及可能引起的误差，并提出如何处理节点功率相关性、故障列表和调度策略等问题，使得半不变量法更加实用化。

半不变量法具备的最大优势在于计算方法简单、计算效率高，虽然计算结果精确度存在一定争议，但满足工程应用要求，因此受到广泛研究。关于半不变量法的算法改进研究主要集中在３个问题

—

１

３

１

—

·综述·　刘　宇，等　电力系统概率潮流算法综述

上：一是如何处理输入随机变量之间的相关性；二是

如何正确分析静态安全稳定问题；三是如何更为精确地描述系统运行状态。文献［４９］提出一种基于Ｃｈｏｌｅｓｋｙ分解的计及输入变量相关性的半不变量法概率潮流计算方法，并提出基于蒙特卡洛抽样的方法解决一些输入变量的半不变量难以被常规数值

方法求解的问题。文献［５０］以受时域和地域不均匀分布影响的系统当前运行状态为基础构造预想事故

集，

将全概率理论和半不变量法概率潮流相结合，提出一种考虑电网静态安全风险的概率潮流算法。文献［５１］在当前结合半不变量和级数展开的概率潮流计算方法基础上，采用分段线性化手段减少潮流方

程线性化误差，

提高系统状态量的求解精度，并引入Ｃ型Ｇｒａｍ－Ｃｈａｒｌｉｅｒ级数改进概率密度函数准确

性。

４．３　序列运算理论

为了将原本复杂的卷积计算转换为简单的算术

运算，康重庆教授创建了序列运算理论［５２］

，在此基

础上扩展了概率性序列运算方法［５３］

，从而形成了全

新的电力系统不确定性分析框架。文献［５４

］提出了基于序列运算的概率直流潮流计算方法。该理论中，将输电元件（包括输电线路和变压器）定义为双向元件，将发电机组和节点负荷定义为单向元件，通

过确定性直流潮流方程（式（５））计算输电线路潮流分布。

Ｐｌ，Ｔ＝∑

ＮＢ

ｋ＝１

Ｇｌ?ｋＰｋ（５

）式中：ｌ＝１，２，…，ＮＴ，其中ＮＴ为输电元件数；Ｐｌ，Ｔ

为输电元件ｌ的有功潮流；Ｐｋ为节点ｋ的有功功率注入；Ｇｌ－ｋ为节点ｋ与输电元件ｌ之间的功率转移分布因子；ＮＢ为节点数。直流潮流中变量的概率特性序列化与自定义的序乘、卷和运算是该概率潮流算法的关键。具体计算步骤可以概括如下。

１

）通过离散化步长将节点功率Ｐｋ和分布因子Ｇｌ－ｋ作序列化，

形成相应的概率性序列。２）根据式（５

）通过序乘运算将节点ｋ的概率功率转换为输电元件ｌ的单节点潮流贡献。

３）考虑到所有节点的影响，通过卷和运算计算输电元件ｌ上的累积潮流。其中，若单个节点有多台发电机组或负荷，则需要在相应变量序乘运算后进行卷和运算，从而获得单个节点对某一输电元件的有功潮流贡献。

基于序列运算的概率直流潮流是以简化序列卷积为出发点，其自定义的卷和、序乘等运算方法计算简单，在效率上具有很大的优势。但由于序列的建立和运算的定义都要满足全新的规则和要求，从而

限制了该方法的大规模推广应用。

４．４　解析法概率潮流展望

解析法概率潮流的研究重点在卷积运算的高效处理上，快速傅里叶运算难以处理大规模电力系统的多变量计算，

序列运算理论体系架构的特殊性使之难以在短时间内大规模推广，因此半不变量法是解析法中研究的热点算法。半不变量法在处理正态分布的变量相关性上已有较多的研究，该算法的研究重点是非正态分布的变量相关性处理，

改善各阶矩和半不变量的计算精度与计算效率，获得更为准确的概率密度函数。

５　结语

电力系统中不确定因素的增加使得概率潮流成

为研究热点之一。本文从概率潮流算法原理的角度出发，将电力系统概率潮流算法分为模拟法、近似法

和解析法３类，

对各算法的研究现状和发展进行了综述和总结，并对各类算法的发展趋势提供了分析

思路。

概率潮流与确定性潮流本质上的区别在于输入变量的不确定性，计算目标是合理反映随机分布本身的概率特性对于系统潮流的影响，因此，概率潮流的研究目的可概括为随机变量的概率输入及其对潮流输出的概率贡献。模拟法分析输入随机变量的抽样方法，

近似法关注输出与输入之间的转化关系，解析法重点讨论系统方程卷积的计算与简化，而３类方法中随机变量概率分布的表征和相关性分析都是关键因素。

面对新环境下电力系统的全新特征，未来概率潮流算法研究的总体发展趋势应该包含以下几个主要方向：①可再生能源功率输出的概率密度函数并不一定精确满足特定的典型分布，

如何生成适用于概率潮流计算的准确的可再生能源功率输出的概率

分布非常重要，目前已有的核密度估计法［３４］

具有很大的研究空间；②在大规模新能源集中并网的趋势下，如何处理区域内新能源功率输出的复杂相关性变得日益重要；③电动汽车作为电网互动性负荷，与电网的功率交互是双向的，如何对双向功率负荷进行概率建模也是研究难点之一；④电力市场机制下

的电源和负荷具备开放性特点［１１］

，基于市场特征和

需求响应的最优概率潮流分析［５５－５６

］具备很高的研究意义。

参考文献

［１］ＯＫＥ　Ｏ　Ａ，ＴＨＯＭＡＳ　Ｄ　Ｗ　Ｐ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｉｎ　ｍｉｃｒｏ－

—

２３１—２０１４，３８（２３

）　

ｇｒｉｄ　ａｓｓｅｓｓｍｅｎｔ　ａｎｄ　ｐｌａｎｎｉｎｇ　ｓｔｕｄｉｅｓ［Ｃ］／／ＩＥＥＥ　ＥｌｅｃｔｒｉｃａｌＰｏｗｅｒ　ａｎｄ　Ｅｎｅｒｇｙ　Ｃｏｎｆｅｒｅｎｃｅ，Ｏｃｔｏｂｅｒ　１０－１２，２０１２，Ｌｏｎｄｏｎ，Ｂｒｉｔａｉｎ：１５１－１５６．

［２］ＢＯＲＫＯＷＳＫＡ　Ｂ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎＰｏｗｅｒ　Ａｐｐａｒａｔｕｓ　ａｎｄ　Ｓｙｓｔｅｍｓ，１９７４，９３（３）：７５２－７５９．

［３］ＤＯＰＡＺＯ　Ｊ　Ｆ，ＫＬＩＴＩＮ　Ｏ　Ａ，ＳＡＳＳＯＮ　Ａ　Ｍ．Ｓｔｏｃｈａｓｔｉｃ　ｌｏａｄｆｌｏｗｓ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ａｐｐａｒａｔｕｓ　ａｎｄ　Ｓｙｓｔｅｍｓ，１９７５，９４（２）：２９９－３０９．

［４］丁明，李生虎，黄凯．基于蒙特卡洛模拟的概率潮流计算［Ｊ］．电网技术，２００１，２５（１１）：１０－１５．

ＤＩＮＧ　Ｍｉｎｇ，ＬＩ　Ｓｈｅｎｇｈｕ，ＨＵＡＮＧ　Ｋａｉ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗａｎａｌｙｓｉｓ　ｂａｓｅｄ　ｏｎ　Ｍｏｎｔｅ－Ｃａｒｌｏ　ｓｉｍｕｌａｔｉｏｎ［Ｊ］．Ｐｏｗｅｒ　ＳｙｓｔｅｍＴｅｃｈｎｏｌｏｇｙ，２００１，２５（１１）：１０－１５．

［５］ＡＬＬＡＮ　Ｒ　Ｎ，ＡＬＳＨＡＫＡＲＣＨＡＩ　Ｍ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　Ａ．Ｃ．ｌｏａｄｆｌｏｗ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ｔｈｅ　Ｉｎｓｔｉｔｕｔｉｏｎ　ｏｆ　Ｅｌｅｃｔｒｉｃａｌ　Ｅｎｇｉｎｅｅｒｓ，１９７６，１２３（６）：５３１－５３６．

［６］ＡＬＬＡＮ　Ｒ　Ｎ，ＬＥＩＴＥ　ＤＡ　ＳＩＬＶＡ　Ａ　Ｍ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗｕｓｉｎｇ　ｍｕｌｔｉｌｉｎｅａｒｉｓａｔｉｏｎｓ［Ｊ］．ＩＥＥ　Ｐｒｏｃｅｅｄｉｎｇｓ　Ｃ：Ｇｅｎｅｒａｔｉｏｎ，Ｔｒａｎｓｍｉｓｓｉｏｎ　ａｎｄ　Ｄｉｓｔｒｉｂｕｔｉｏｎ，１９８１，１２８（５）：２８０－２８７．

［７］ＳＯＫＩＥＲＡＪＳＫＩ　Ｍ　Ａ．Ｍｅｔｈｏｄ　ｏｆ　ｓｔｏｃｈａｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｃａｌｃｕｌａｔｉｏｎ［Ｊ］．Ａｃｈｉｖ　Ｆｕｒ　Ｅｌｅｋｔｒｏｔｅｃｈｎｉｋ，１９７８，６０（２）：３７－４０．

［８］陈雁，文劲宇，程时杰．考虑输入变量相关性的概率潮流计算方法［Ｊ］．中国电机工程学报，２０１１，３１（２２）：８０－８７．

ＣＨＥＮ　Ｙａｎ，ＷＥＮ　Ｊｉｎｙｕ，ＣＨＥＮＧ　Ｓｈｉｊｉｅ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄｆｌｏｗ　ａｎａｌｙｓｉｓ　ｃｏｎｓｉｄｅｒｉｎｇ　ｄｅｐｅｎｄｅｎｃｉｅｓ　ａｍｏｎｇ　ｉｎｐｕｔ　ｒａｎｄｏｍｖａｒｉａｂｌｅｓ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ｔｈｅ　ＣＳＥＥ，２０１１，３１（２２）：８０－８７．［９］邓威，李欣然，刘志勇，等．考虑无功补偿影响因素的间歇性分布式电源综合优化配置［Ｊ］．中国电机工程学报，２０１２，３２（１０）：８０－８８．

ＤＥＮＧ　Ｗｅｉ，ＬＩ　Ｘｉｎｒａｎ，ＬＩＵ　Ｚｈｉｙｏｎｇ，ｅｔ　ａｌ．Ｃｏｍｐｒｅｈｅｎｓｉｖｅｏｐｔｉｍａｌ　ａｌｌｏｃａｔｉｏｎ　ｏｆ　ｉｎｔｅｒｍｉｔｔｅｎｔ　ｇｅｎｅｒａｔｉｏｎ　ｃｏｎｓｉｄｅｒｉｎｇ　ｒｅａｃｔｉｖｅｐｏｗｅｒ　ｃｏｍｐｅｎｓａｔｉｏｎ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ｔｈｅ　ＣＳＥＥ，２０１２，３２（１０）：８０－８８．

［１０］张传成．概率潮流在静态安全分析中的应用［Ｄ］．北京：华北电力大学，２００８．

［１１］姚建国，杨胜春，王珂，等．智能电网“源－网－荷”互动运行控制概念及研究框架［Ｊ］．电力系统自动化，２０１２，３６（２１）：１－６．

ＹＡＯ　Ｊｉａｎｇｕｏ，ＹＡＮＧ　Ｓｈｅｎｇｃｈｕｎ，ＷＡＮＧ　Ｋｅ，ｅｔ　ａｌ．Ｃｏｎｃｅｐｔａｎｄ　ｒｅｓｅａｒｃｈ　ｆｒａｍｅｗｏｒｋ　ｏｆ　ｓｍａｒｔ　ｇｒｉｄ“ｓｏｕｒｃｅ－ｇｒｉｄ－ｌｏａｄ”

ｉｎｔｅｒａｃｔｉｖｅ　ｏｐｅｒａｔｉｏｎ　ａｎｄ　ｃｏｎｔｒｏｌ［Ｊ］．Ａｕｔｏｍａｔｉｏｎ　ｏｆ　ＥｌｅｃｔｒｉｃＰｏｗｅｒ　Ｓｙｓｔｅｍｓ，２０１２，３６（２１）：１－６．

［１２］ＴＡＭＴＵＭ　Ａ，ＳＣＨＥＬＬＥＮＢＥＲＧ　Ａ，ＲＯＳＥＨＡＲＴ　Ｗ　Ｄ．Ｅｎｈａｎｃｅｍｅｎｔｓ　ｔｏ　ｔｈｅ　ｃｕｍｕｌａｎｔ　ｍｅｔｈｏｄ　ｆｏｒ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｏｐｔｉｍａｌｐｏｗｅｒ　ｆｌｏｗ　ｓｔｕｄｉｅｓ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙｓｔｅｍｓ，２００９，２４（４）：１７３９－１７４６．

［１３］ＤＯＮＤＥＲＡ　Ｄ，ＰＯＰＡ　Ｒ，ＶＥＬＩＣＥＳＣＵ　Ｃ．Ｔｈｅ　ｍｕｌｔｉ－ａｒｅａｓｙｓｔｅｍｓ　ｒｅｌｉａｂｉｌｉｔｙ　ｅｓｔｉｍａｔｉｏｎ　ｕｓｉｎｇ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｂｙＧｒａｍｍ－Ｃｈａｒｌｉｅｒ　ｅｘｐａｎｓｉｏｎ［Ｃ］／／Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｃｏｎｆｅｒｅｎｃｅ　ｏｎＣｏｍｐｕｔｅｒ　ａｓ　ａ　Ｔｏｏｌ，Ｓｅｐｔｅｍｂｅｒ　９－１２，２００７，Ｗａｒｓａｗ，Ｐｏｌａｎｄ：１４７０－１４７４．

［１４］ＬＩ　Ｊ　Ｆ，ＺＨＡＮＧ　Ｂ　Ｈ，ＬＩＵ　Ｙ　Ｆ．Ｄａｔａ　ｍｉｎｉｎｇ　ｉｎ　ｎｏｎｌｉｎｅａｒｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｂａｓｅｄ　ｏｎ　Ｍｏｎｔｅ　Ｃａｒｌｏ　ｓｉｍｕｌａｔｉｏｎ［Ｃ］／／Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　１ｓｔ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｃｏｎｆｅｒｅｎｃｅ　ｏｎ　ＩｎｆｏｒｍａｔｉｏｎＳｃｉｅｎｃｅ　ａｎｄ　Ｅｎｇｉｎｅｅｒｉｎｇ，Ｄｅｃｅｍｂｅｒ　２６－２８，２００９，Ｎａｎｊｉｎｇ，Ｃｈｉｎａ：８３３－８３６．［１５］丁明，吴兴龙，陆巍，等．含多个不对称光伏并网系统的配电网三相潮流随机计算［Ｊ］．电力系统自动化，２０１２，３６（１６）：４７－５２．

ＤＩＮＧ　Ｍｉｎｇ，ＷＵ　Ｘｉｎｇｌｏｎｇ，ＬＵ　Ｗｅｉ，ｅｔ　ａｌ．Ｔｈｒｅｅ－ｐｈａｓｅｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｃａｌｃｕｌａｔｉｏｎ　ｉｎ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ｓｙｓｔｅｍｓ　ｗｉｔｈｍｕｌｔｉｐｌｅ　ｕｎｓｙｍｍｅｔｒｉｃａｌ　ｇｒｉｄ－ｃｏｎｎｅｃｔｅｄ　ｐｈｏｔｏｖｏｌｔａｉｃ　ｓｙｓｔｅｍｓ［Ｊ］．Ａｕｔｏｍａｔｉｏｎ　ｏｆ　Ｅｌｅｃｔｒｉｃ　Ｐｏｗｅｒ　Ｓｙｓｔｅｍｓ，２０１２，３６（１６）：４７－５２．

［１６］邓威，李欣然，徐振华，等．考虑风速相关性的概率潮流计算及影响分析［Ｊ］．电网技术，２０１２，３６（４）：４５－５０．

ＤＥＮＧ　Ｗｅｉ，ＬＩ　Ｘｉｎｒａｎ，ＸＵ　Ｚｈｅｎｈｕａ，ｅｔ　ａｌ．Ｃａｌｃｕｌａｔｉｏｎ　ｏｆｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｃｏｎｓｉｄｅｒｉｎｇ　ｗｉｎｄ　ｓｐｅｅｄ　ｃｏｒｒｅｌａｔｉｏｎ　ａｎｄａｎａｌｙｓｉｓ　ｏｎ　ｉｎｆｌｕｅｎｃｅ　ｏｆ　ｗｉｎｄ　ｓｐｅｅｄ　ｃｏｒｒｅｌａｔｉｏｎ［Ｊ］．ＰｏｗｅｒＳｙｓｔｅｍ　Ｔｅｃｈｎｏｌｏｇｙ，２０１２，３６（４）：４５－５０．

［１７］ＨＵＡＮＧ　Ｊ，ＸＵＥ　Ｙ　Ｓ，ＤＯＮＧ　Ｚ　Ｙ，ｅｔ　ａｌ．Ａｎ　ａｄａｐｔｉｖｅｉｍｐｏｒｔａｎｃｅ　ｓａｍｐｌｉｎｇ　ｍｅｔｈｏｄ　ｆｏｒ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｏｐｔｉｍａｌ　ｐｏｗｅｒｆｌｏｗ［Ｃ］／／ＩＥＥＥ　Ｐｏｗｅｒ　ａｎｄ　Ｅｎｅｒｇｙ　Ｓｏｃｉｅｔｙ　Ｇｅｎｅｒａｌ　Ｍｅｅｔｉｎｇ，Ｊｕｌｙ　２４－２９，２０１１，Ｄｅｔｒｏｉｔ，ＵＳＡ：６ｐ．

［１８］宋晓通，谭震宇．改进的重要抽样法在电力系统可靠性评估中的应用［Ｊ］．电网技术，２００５，２９（１３）：５６－６０．

ＳＯＮＧ　Ｘｉａｏｔｏｎｇ，ＴＡＮ　Ｚｈｅｎｙｕ．Ａｐｐｌｉｃａｔｉｏｎ　ｏｆ　ｉｍｐｒｏｖｅｄｉｍｐｏｒｔａｎｃｅ　ｓａｍｐｌｉｎｇ　ｍｅｔｈｏｄ　ｉｎ　ｐｏｗｅｒ　ｓｙｓｔｅｍ　ｒｅｌｉａｂｉｌｉｔｙｅｖａｌｕａｔｉｏｎ［Ｊ］．Ｐｏｗｅｒ　Ｓｙｓｔｅｍ　Ｔｅｃｈｎｏｌｏｇｙ，２００５，２９（１３）：５６－６０．

［１９］周阳洋．重要抽样法在概率潮流中的应用［Ｄ］．北京：华北电力大学，２０１１．

［２０］ＭＣＫＡＹ　Ｍ　Ｄ，ＢＥＣＫＭＡＮ　Ｒ　Ｊ，ＣＯＮＯＶＥＲ　Ｗ　Ｊ．Ａｃｏｍｐａｒｉｓｏｎ　ｏｆ　ｔｈｒｅｅ　ｍｅｔｈｏｄｓ　ｆｏｒ　ｓｅｌｅｃｔｉｎｇ　ｖａｌｕｅｓ　ｏｆ　ｉｎｐｕｔｖａｒｉａｂｌｅｓ　ｉｎ　ｔｈｅ　ａｎａｌｙｓｉｓ　ｏｆ　ｏｕｔｐｕｔ　ｆｒｏｍ　ａ　ｃｏｍｐｕｔｅｒ　ｃｏｄｅ［Ｊ］．Ｔｅｃｈｎｉｍｅｔｒｉｃｓ，１９７９，２１（２）：２３９－２４５．

［２１］于晗，钟志勇，黄杰波，等．采用拉丁超立方采样的电力系统概率潮流计算方法［Ｊ］．电力系统自动化，２００９，３３（２１）：３２－３５．

ＹＵ　Ｈａｎ，ＣＨＵＮＧ　Ｃｈｉｙｏｎｇ，ＷＯＮＧ　Ｋｉｔｐｏ，ｅｔ　ａｌ．Ａｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｃａｌｃｕｌａｔｉｏｎ　ｍｅｔｈｏｄ　ｗｉｔｈ　Ｌａｔｉｎ　ｈｙｐｅｒｃｕｂｅｓａｍｐｌｉｎｇ［Ｊ］．Ａｕｔｏｍａｔｉｏｎ　ｏｆ　Ｅｌｅｃｔｒｉｃ　Ｐｏｗｅｒ　Ｓｙｓｔｅｍｓ，２００９，３３（２１）：３２－３５．

［２２］李俊芳，张步涵．基于进化算法改进拉丁超立方抽样的概率潮流计算［Ｊ］．中国电机工程学报，２０１１，３１（２５）：９０－９６．

ＬＩ　Ｊｕｎｆａｎｇ，ＺＨＡＮＧ　Ｂｕｈａｎ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｂａｓｅｄ　ｏｎｉｍｐｒｏｖｅｄ　Ｌａｔｉｎ　ｈｙｐｅｒｃｕｂｅ　ｓａｍｐｌｉｎｇ　ｗｉｔｈ　ｅｖｏｌｕｔｉｏｎａｒｙａｌｇｏｒｉｔｈｍ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ＣＳＥＥ，２０１１，３１（２５）：９０－９６．［２３］ＹＵ　Ｈ，ＣＨＵＮＧ　Ｃ　Ｙ，ＷＯＮＧ　Ｋ　Ｐ，ｅｔ　ａｌ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄｆｌｏｗ　ｅｖａｌｕａｔｉｏｎ　ｗｉｔｈ　ｈｙｂｒｉｄ　Ｌａｔｉｎ　ｈｙｐｅｒｃｕｂｅ　ｓａｍｐｌｉｎｇ　ａｎｄＣｈｏｌｅｓｋｙ　ｄｅｃｏｍｐｏｓｉｔｉｏｎ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙｓｔｅｍｓ，２００９，２４（２）：６６１－６６７．

［２４］丁明，王京景，李生虎．基于扩展拉丁超立方采样的电力系统概率潮流计算［Ｊ］．中国电机工程学报，２０１３，３３（４）：１６３－１７１．

ＤＩＮＧ　Ｍｉｎｇ，ＷＡＮＧ　Ｊｉｎｇｊｉｎｇ，ＬＩ　Ｓｈｅｎｇｈｕ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄｆｌｏｗ　ｅｖａｌｕａｔｉｏｎ　ｗｉｔｈ　ｅｘｔｅｎｄｅｄ　Ｌａｔｉｎ　ｈｙｐｅｒｃｕｂｅ　ｓａｍｐｌｉｎｇ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ｔｈｅ　ＣＳＥＥ，２０１３，３３（４）：１６３－１７１．

［２５］ＣＡＩ　Ｄ　Ｆ，ＳＨＩ　Ｄ　Ｙ，ＣＨＥＮ　Ｊ　Ｆ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗｃｏｍｐｕｔａｔｉｏｎ　ｗｉｔｈ　ｐｏｌｙｎｏｍｉａｌ　ｎｏｒｍａｌ　ｔｒａｎｓｆｏｒｍａｔｉｏｎ　ａｎｄ　Ｌａｔｉｎｈｙｐｅｒｃｕｂｅ　ｓａｍｐｌｉｎｇ［Ｊ］．ＩＥＴ　Ｇｅｎｅｒａｔｉｏｎ，Ｔｒａｎｓｍｉｓｓｉｏｎ　＆Ｄｉｓｔｒｉｂｕｔｉｏｎ，２０１３，７（５）：４７４－４８２．

［２６］ＺＨＯＵ　Ｂ，ＸＩＡＯ　Ｑ．Ｓｏｌｖｉｎｇ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｏｐｔｉｍａｌ　ｐｏｗｅｒ　ｆｌｏｗｐｒｏｂｌｅｍ　ｕｓｉｎｇ　ｑｕａｓｉ　Ｍｏｎｔｅ　Ｃａｒｌｏ　ｍｅｔｈｏｄ　ａｎｄ　ｎｉｎｔｈ－ｏｒｄｅｒ

—

３

３

１

—

·综述·　刘　宇，等　电力系统概率潮流算法综述

ｐｏｌｙ

ｎｏｍｉａｌ　ｎｏｒｍａｌ　ｔｒａｎｓｆｏｒｍａｔｉｏｎ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　ＰｏｗｅｒＳｙ

ｓｔｅｍｓ，２０１４，２９（１）：３００－３０６．［２７］ＨＡＪＩＡＮ　Ｍ，ＲＯＳＥＨＡＲＴ　Ｗ　Ｄ，ＺＡＲＥＩＰＯＵＲ　Ｈ．

Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｐｏｗｅｒ　ｆｌｏｗ　ｂｙ　

Ｍｏｎｔｅ　Ｃａｒｌｏ　ｓｉｍｕｌａｔｉｏｎ　ｗｉｔｈ　Ｌａｔｉｎｓｕｐｅｒｃｕｂｅ　ｓａｍｐｌｉｎｇ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙｓｔｅｍｓ，２０１３，２８（２）：１５５０－

１５５９．［２８］ＨＵＡＮＧ　Ｈ　Ｚ，ＣＨＵＮＧ　Ｃ　Ｙ，ＣＨＡＮ　Ｋ　Ｗ，ｅｔ　ａｌ．Ｑｕａｓｉ－

Ｍｏｎｔｅ　Ｃａｒｌｏ　ｂａｓｅｄ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｓｍａｌｌ　ｓｉｇｎａｌ　ｓｔａｂｉｌｉｔｙ　ａｎａｌｙｓｉｓｆｏｒ　ｐｏｗｅｒ　ｓｙｓｔｅｍ　ｗｉｔｈ　ｐｌｕｇ－ｉｎ　ｅｌｅｃｔｒｉｃ　ｖｅｈｉｃｌｅ　ａｎｄ　ｗｉｎｄ　ｐ

ｏｗｅｒｉｎｔｅｇｒａｔｉｏｎ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙｓｔｅｍｓ，２０１３，２８（３）：３３３５－

３３４３．［２９］ＨＯＮＧ　Ｈ　Ｐ．Ａｎ　ｅｆｆｉｃｉｅｎｔ　ｐｏｉｎｔ　ｅｓｔｉｍａｔｅ　ｍｅｔｈｏｄ　ｆ

ｏｒｐｒｏｂａｂｉｌｉｓｔｉｃ　ａｎａｌｙｓｉｓ［Ｊ］．Ｒｅｌｉａｂｉｌｉｔｙ　Ｅｎｇｉｎｅｅｒｉｎｇ　ａｎｄ　ＳｙｓｔｅｍＳａｆｅｔｙ，１９９８，５９（３）：２６１－２６７．［３０］ＳＵ　Ｃ　Ｌ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ－ｆｌｏｗ　ｃｏｍｐｕｔａｔｉｏｎ　ｕｓｉｎｇ　ｐ

ｏｉｎｔｅｓｔｉｍａｔｅ　ｍｅｔｈｏｄ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙｓｔｅｍｓ，２００５，２０（４）：１８４３－

１８５１．［３１］ＭＯＲＡＬＥＳ　Ｊ　Ｍ，ＰＥＲＥＺ－ＲＵＩＺ　Ｊ．Ｐｏｉｎｔ　ｅｓｔｉｍａｔｅ　ｓｃｈｅｍｅｓ　

ｔｏｓｏｌｖｅ　ｔｈｅ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｐｏｗｅｒ　ｆｌｏｗ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　ＰｏｗｅｒＳｙ

ｓｔｅｍｓ，２００７，２２（４）：１５９４－１６０１．［３２］潘炜，刘文颖，杨以涵．概率最优潮流的点估计算法［Ｊ］．

中国电机工程学报，２００８，２８（１６）：２８－

３３．ＰＡＮ　Ｗｅｉ，ＬＩＵ　Ｗｅｎｙｉｎｇ，ＹＡＮＧ　Ｙｉｈａｎ．Ｐｏｉｎｔ　ｅｓｔｉｍａｔｅｍｅｔｈｏｄ　ｆｏｒ　ｐｒｏｂａｂｉｌｉｓｔｉｃａｌｌｙ　

ｏｐｔｉｍａｌ　ｐｏｗｅｒ　ｆｌｏｗ　ｃｏｍｐｕｔａｔｉｏｎ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ｔｈｅ　ＣＳＥＥ，２００８，２８（１６）：２８－３３．［３３］艾小猛，文劲宇，吴桐，等．基于点估计和Ｇｒａｍ－Ｃｈａｒｌｉｅｒ展开的

含风电电力系统概率潮流实用算法［Ｊ］．中国电机工程学报，２０１３，３３（１６）：１６－

２３．ＡＩ　Ｘｉａｏｍｅｎｇ，ＷＥＮ　Ｊｉｎｙｕ，ＷＵ　Ｔｏｎｇ，ｅｔ　ａｌ．Ａ　ｐｒａｃｔｉｃａｌａｌｇｏｒｉｔｈｍ　ｂａｓｅｄ　ｏｎ　ｐｏｉｎｔ　ｅｓｔｉｍａｔｅ　ｍｅｔｈｏｄ　ａｎｄ　Ｇｒａｍ－Ｃｈａｒｌｉｅｒｅｘｐａｎｓｉｏｎ　ｆｏｒ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｃａｌｃｕｌａｔｉｏｎ　ｏｆ　ｐｏｗｅｒｓｙｓｔｅｍｓ　ｉｎｃｏｒｐｏｒａｔｉｎｇ　

ｗｉｎｄ　ｐｏｗｅｒ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ｔｈｅＣＳＥＥ，２０１３，３３（１６）：１６－２３．［３４］ＳＯＬＥＩＭＡＮＰＯＵＲ　Ｎ，ＭＯＨＡＭＭＡＤＩ　Ｍ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　

ｌｏａｄｆｌｏｗ　ｂｙ　ｕｓｉｎｇ　ｎｏｎｐａｒａｍｅｔｒｉｃ　ｄｅｎｓｉｔｙ　ｅｓｔｉｍａｔｏｒｓ［Ｊ］．ＩＥＥＥＴｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙ

ｓｔｅｍｓ，２０１３，２８（４）：３７４７－３７５５．［３５］ＳＡＵＮＤＥＲＳ　Ｃ　Ｓ．Ｐｏｉｎｔ　ｅｓｔｉｍａｔｅ　ｍｅｔｈｏｄ　ａｄｄｒｅｓｓｉｎｇ　

ｃｏｒｒｅｌａｔｅｄｗｉｎｄ　ｐｏｗｅｒ　ｆｏｒ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｏｐｔｉｍａｌ　ｐｏｗｅｒ　ｆｌｏｗ［Ｊ］．ＩＥＥＥＴｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙ

ｓｔｅｍｓ，２０１４，２９（３）：１０４５－１０５４．［３６］ＬＩ　Ｘ，ＬＩ　Ｙ　Ｚ，ＺＨＡＮＧ　Ｓ　Ｈ．Ａｎａｌｙｓｉｓ　ｏｆ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｏｐ

ｔｉｍａｌｐｏｗｅｒ　ｆｌｏｗ　ｔａｋｉｎｇ　

ａｃｃｏｕｎｔ　ｏｆ　ｔｈｅ　ｖａｒｉａｔｉｏｎ　ｏｆ　ｌｏａｄ　ｐｏｗｅｒ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙｓｔｅｍｓ，２００８，２３（３）：９９２－９９９．［３７］ＷＡＮ　Ｃ，ＸＵ　Ｚ，ＤＯＮＧ　Ｚ　Ｙ，ｅｔ　ａｌ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　

ｆｌｏｗｃｏｍｐｕｔａｔｉｏｎ　ｕｓｉｎｇ　ｆｉｒｓｔ－ｏｒｄｅｒ　ｓｅｃｏｎｄ－ｍｏｍｅｎｔ　ｍｅｔｈｏｄ［Ｃ］／／ＩＥＥＥ　Ｐｏｗｅｒ　ａｎｄ　Ｅｎｅｒｇｙ　Ｓｏｃｉｅｔｙ　Ｇｅｎｅｒａｌ　Ｍｅｅｔｉｎｇ，Ｊｕｌｙ　２２－２６，２０１２，Ｓａｎ　Ｄｉｅｇｏ，Ｃａｎａｄａ：８ｐ

．［３８］李雪，李渝曾，李海英．几种概率潮流算法的比较与分析［Ｊ］．

电力系统及其自动化学报，２００９，２１（３）：１２－

１７．ＬＩ　Ｘｕｅ，ＬＩ　Ｙｕｚｅｎｇ，ＬＩ　Ｈａｉｙｉｎｇ．Ｃｏｍｐａｒｉｓｏｎ　ａｎｄ　ａｎａｌｙｓｉｓ　ｏｆｓｅｖｅｒａｌ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｐｏｗｅｒ　ｆｌｏｗ　ａｌｇｏｒｉｔｈｍｓ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆｔｈｅ　ＣＳＵ－ＥＰＳＡ，２００９，２１（３）：１２－

１７．［３９

］刘艳丽．随机潮流中正态变量线性变换不变性定理应用的研究［Ｄ］．天津：天津大学，２００９．

［４０］ＪＵＬＩＥＲ　Ｓ　Ｊ，ＵＨＬＭＡＮＮ　Ｊ　Ｋ．Ｕｎｓｃｅｎｔｅｄ　ｆｉｌｔｅｒｉｎｇ　

ａｎｄｎｏｎｌｉｎｅａｒ　ｅｓｔｉｍａｔｉｏｎ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇ

ｓ　ｏｆ　ｔｈｅ　ＩＥＥＥ，２００４，９２（３）：４０１－

４２２．［４１］ＭＯＲＴＥＺＡ　Ａ，ＭＡＨＭＵＤ　Ｆ　

Ｆ，ＦＡＲＲＯＫＨ　Ａ．Ｐｒｏｂａｂｉｌｉｓｔｉｃｌｏａｄ　ｆｌｏｗ　ｉｎ　ｃｏｒｒｅｌａｔｅｄ　ｕｎｃｅｒｔａｉｎ　ｅｎｖｉｒｏｎｍｅｎｔ　ｕｓｉｎｇ　ｕｎｓｃｅｎｔｅｄｔｒａｎｓｆｏｒｍａｔｉｏｎ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙ

ｓｔｅｍｓ，２０１２，２７（４）：２２３３－２２４１．［４２］ＬＥＩＴＥ　ＤＡ　ＳＩＬＶＡ　Ａ　Ｍ，ＡＲＥＩＥＮＴＩ　Ｖ　Ｌ，ＡＬＬＡＮ　

Ｒ　Ｎ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｃｏｎｓｉｄｅｒｉｎｇ　ｄｅｐｅｎｄｅｎｃｅ　ｂｅｔｗｅｅｎ　ｉｎｐｕｔｎｏｄａｌ　ｐｏｗｅｒｓ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ａｐｐａｒａｔｕｓ　ａｎｄＳｙ

ｓｔｅｍｓ，１９８４，１０３（６）：１５２４－１５３０．［４３］ＫＯＬＢＡ　Ｄ，ＰＡＲＫＳ　Ｔ　Ｗ．Ａ　ｐ

ｒｉｍｅ　ｆａｃｔｏｒ　ＦＦＴ　ａｌｇｏｒｉｔｈｍ　ｕｓｉｎｇｈｉｇｈ－ｓｐｅｅｄ　ｃｏｎｖｏｌｕｔｉｏｎ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ａｃｏｕｓｔｉｃｓ，Ｓｐ

ｅｅｃｈａｎｄ　Ｓｉｇｎａｌ　Ｐｒｏｃｅｓｓｉｎｇ，１９７７，２５（４）：２８１－２９４．［４４］王锡凡，王秀丽．电力系统的随机潮流分析［Ｊ］．

西安交通大学学报，１９８８，２２（３）：８７－

９７．ＷＡＮＧ　Ｘｉｆａｎ，ＷＡＮＧ　Ｘｉｕｌｉ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ａｎａｌｙｓｉｓ　ｉｎｐｏｗｅｒ　ｓｙｓｔｅｍｓ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｘｉ’ａｎ　Ｊｉａｏｔｏｎｇ　

Ｕｎｉｖｅｒｓｉｔｙ，１９８８，２２（３）：８７－９７．［４５］ＺＨＡＮＧ　Ｐ，ＬＥＥ　Ｓ　Ｔ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｃｏｍｐ

ｕｔａｔｉｏｎｕｓｉｎｇ　ｔｈｅ　ｍｅｔｈｏｄ　ｏｆ　ｃｏｍｂｉｎｅｄ　ｃｕｍｕｌａｎｔｓ　ａｎｄ　Ｇｒａｍ－Ｃｈａｒｌｉｅｒｅｘｐａｎｓｉｏｎ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙｓｔｅｍｓ，２００４，１９（１）：６７６－

６８２．［４６］ＳＣＨＥＬＬＥＮＢＥＲＧ　Ａ，ＲＯＳＥＨＡＲＴ　Ｗ，ＡＧＵＡＤＯ　Ｊ

．Ｃｕｍｕｌａｎｔ－ｂａｓｅｄ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｏｐｔｉｍａｌ　ｐ

ｏｗｅｒ　ｆｌｏｗ（Ｐ－ＯＰＦ）ｗｉｔｈ　Ｇａｕｓｓｉａｎ　ａｎｄ　Ｇａｍｍａ　ｄｉｓｔｒｉｂｕｔｉｏｎｓ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎＰｏｗｅｒ　Ｓｙ

ｓｔｅｍｓ，２００５，２０（２）：７７３－７８１．［４７］董雷，程卫东，杨以涵．含风电场的电力系统概率潮流计算［Ｊ］．

电网技术，２００９，３３（１６）：８７－

９１．ＤＯＮＧ　Ｌｅｉ，ＣＨＥＮＧ　Ｗｅｉｄｏｎｇ

，ＹＡＮＧ　Ｙｉｈａｎ．Ｐｒｏｂａｂｉｌｉｓｔｉｃｌｏａｄ　ｆｌｏｗ　ｃａｌｃｕｌａｔｉｏｎ　ｆｏｒ　ｐｏｗｅｒ　ｇｒｉｄ　ｃｏｎｔａｉｎｉｎｇ　ｗｉｎｄ　ｆａｒｍｓ［Ｊ］．Ｐｏｗｅｒ　Ｓｙｓｔｅｍ　Ｔｅｃｈｎｏｌｏｇｙ

，２００９，３３（１６）：８７－９１．［４８］胡泽春，王锡凡．基于半不变量法的随机潮流误差分析［Ｊ］．

电网技术，２００９，３３（１８）：３２－

３７．ＨＵ　Ｚｅｃｈｕｎ，ＷＡＮＧ　Ｘｉｆａｎ．Ｅｒｒｏｒ　ａｎａｌｙｓｉｓ　ｏｆ　ｔｈｅ　ｐｒｏｂａｂｉｌｉｓｔｉｃｌｏａｄ　ｆｌｏｗ　ｂａｓｅｄ　ｏｎ　ｃｕｍｕｌａｎｔ　ｍｅｔｈｏｄ［Ｊ］．Ｐｏｗｅｒ　ＳｙｓｔｅｍＴｅｃｈｎｏｌｏｇｙ

，２００９，３３（１８）：３２－３７．［４９

］石东源，蔡德福，陈金富，等．计及输入变量相关性的半不变量法概率潮流［Ｊ］．中国电机工程学报，２０１２，３２（２８）：１０４－１１４．ＳＨＩ　Ｄｏｎｇｙｕａｎ，ＣＡＩ　Ｄｅｆｕ，ＣＨＥＮ　Ｊｉｎｆｕ，ｅｔ　ａｌ．Ｐｒｏｂａｂｉｌｉｓｔｉｃｌｏａｄ　ｆｌｏｗ　ｃａｌｃｕｌａｔｉｏｎ　ｂａｓｅｄ　ｏｎ　ｃｕｍｕｌａｎｔ　ｍｅｔｈｏｄ　ｃｏｎｓｉｄｅｒｉｎｇｃｏｒｒｅｌａｔｉｏｎ　ｂｅｔｗｅｅｎ　ｉｎｐｕｔ　ｖａｒｉａｂｌｅｓ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ｔｈｅＣＳＥＥ，２０１２，３２（２８）：１０４－

１１４．［５０

］刘怡芳，张步涵，李俊芳，等．考虑电网静态安全风险的随机潮流计算［Ｊ］．中国电机工程学报，２０１１，３１（１）：５９－

６４．ＬＩＵ　Ｙｉｆａｎｇ，ＺＨＡＮＧ　Ｂｕｈａｎ，ＬＩ　Ｊｕｎｆａｎｇ，ｅｔ　ａｌ．Ｐｒｏｂａｂｉｌｉｓｔｉｃｌｏａｄ　ｆｌｏｗ　ａｌｇｏｒｉｔｈｍ　ｃｏｎｓｉｄｅｒｉｎｇ　ｓｔａｔｉｃ　ｓｅｃｕｒｉｔｙ　ｒｉｓｋ　ｏｆ　ｔｈｅ　ｐｏｗｅｒｓｙｓｔｅｍ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇ

ｓ　ｏｆ　ＣＳＥＥ，２０１１，３１（１）：５９－６４．［５１

］朱星阳，刘文霞，张建华．考虑大规模风电并网的电力系统随机潮流［Ｊ］．中国电机工程学报，２０１３，３３（７）：７７－

８６．ＺＨＵ　Ｘｉｎｇｙａｎｇ，ＬＩＵ　Ｗｅｎｘｉａ，ＺＨＡＮＧ　Ｊｉａｎｈｕａ．Ｐｒｏｂａｂｉｌｉｓｔｉｃｌｏａｄ　ｆｌｏｗ　ｍｅｔｈｏｄ　ｃｏｎｓｉｄｅｒｉｎｇ　ｌａｒｇｅ－ｓｃａｌｅ　ｗｉｎｄ　ｐｏｗｅｒｉｎｔｅｇｒａｔｉｏｎ［Ｊ］．Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ｔｈｅ　ＣＳＥＥ，２０１３，３３（７）：７７－８６．［５２］康重庆，夏清，相年德，等．序列运算理论及其应用［Ｊ］．

电力系统自动化，２００２，２６（１７）：６－

１１．ＫＡＮＧ　Ｃｈｏｎｇｑｉｎｇ，ＸＩＡ　Ｑｉｎｇ，ＸＩＡＮＧ　Ｎｉａｎｄｅ，ｅｔ　ａｌ．Ｓｅｑｕｅｎｃｅ　ｏｐｅｒａｔｉｏｎ　ｔｈｅｏｒｙ　

ａｎｄ　ｉｔｓ　ａｐｐｌｉｃａｔｉｏｎｓ［Ｊ］．Ａｕｔｏｍａｔｉｏｎ—

４３１—２０１４，３８（２３

）　

ｏｆ　Ｅｌｅｃｔｒｉｃ　Ｐｏｗｅｒ　Ｓｙ

ｓｔｅｍｓ，２００２，２６（１７）：６－１１．［５３］徐玮，康重庆，夏清．概率性序列运算的扩展及其运算性质［Ｊ］．

清华大学学报：自然科学版，２０１０，５０（１）：３１－

３４．ＸＵ　Ｗｅｉ，ＫＡＮＧ　Ｃｈｏｎｇｑｉｎｇ，ＸＩＡ　Ｑｉｎｇ．Ｅｘｔｅｎｓｉｏｎｓ　ａｎｄｐｒｏｐｅｒｔｉｅｓ　ｏｆ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｓｅｑｕｅｎｃｅ　ｏｐ

ｅｒａｔｉｏｎｓ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆＴｓｉｎｇｈｕａ　Ｕｎｉｖｅｒｓｉｔｙ：Ｓｃｉｅｎｃｅ　＆Ｔｅｃｈｎｏｌｏｇｙ，２０１０，５０（１）：３１－

３４．［５４］康重庆，夏清，徐玮．电力系统不确定性分析［Ｍ］．

北京：科学出版社，２０１１．

［５５］ＹＵ　Ｐ，ＶＥＮＫＡＴＥＳＨ　Ｂ．Ａ　ｐ

ｒａｃｔｉｃａｌ　ｒｅａｌ－ｔｉｍｅ　ＯＰＦ　ｍｅｔｈｏｄｕｓｉｎｇ　ｎｅｗ　ｔｒｉａｎｇｕｌａｒ　ａｐｐｒｏｘｉｍａｔｅ　ｍｏｄｅｌ　ｏｆ　ｗｉｎｄ　ｅｌｅｃｔｒｉｃｇ

ｅｎｅｒａｔｏｒｓ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙｓｔｅｍｓ，２０１２，２７（４）：２０３６－

２０４６．［５６］ＳＩＡＮＯ　Ｐ，ＭＯＫＲＹＡＮＩ　Ｇ．Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ａｓｓｅｓｓｍｅｎｔ　ｏｆ　

ｔｈｅｉｍｐａｃｔ　ｏｆ　ｗｉｎｄ　ｅｎｅｒｇｙ　ｉｎｔｅｇｒａｔｉｏｎ　ｉｎｔｏ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ｎｅｔｗｏｒｋｓ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｐｏｗｅｒ　Ｓｙ

ｓｔｅｍｓ，２０１３，２８（４）：４２０９－４２１７．刘　宇（１９９０—），男，博士研究生，主要研究方向：电力系统运行规划、新能源并网。Ｅ－ｍａｉｌ：ｌｉｕｙ

ｕ＿ｓｅｕ＠１２６．ｃｏｍ高　山（１９７３—），男，通信作者，副教授，博士生导师，主要研究方向：

电力系统运行规划、新能源并网、柔性输电技术在电力系统中的应用、负荷预测。Ｅ－ｍａｉｌ：ｓｈａｎｇａｏ＠ｓｅｕ．ｅｄｕ．ｃｎ

杨胜春（１９７３—），男，研究员级高级工程师，主要研究方向：智能电网及电网调度自动化。Ｅ－ｍａｉｌ：ｙａｎｇｓｈｅｎｇｃｈｕｎ＠ｅｐｒｉ．ｓｇ

ｃｃ．ｃｏｍ．ｃｎ（编辑　万志超）

Ｒｅｖｉｅｗ　ｏｎ　Ａｌｇｏｒｉｔｈｍｓ　ｆｏｒ　Ｐｒｏｂａｂｉｌｉｓｔｉｃ　Ｌｏａｄ　Ｆｌｏｗ　ｉｎ　Ｐｏｗｅｒ　Ｓｙ

ｓｔｅｍＬＩＵ　Ｙｕ１，

２，ＧＡＯ　Ｓｈａｎ１，

２，ＹＡＮＧ　Ｓｈｅｎｇｃｈｕｎ３，ＹＡＯ　Ｊｉａｎｇｕｏ３

（１．Ｓｃｈｏｏｌ　ｏｆ　Ｅｌｅｃｔｒｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ，Ｓｏｕｔｈｅａｓｔ　Ｕｎｉｖｅｒｓｉｔｙ，Ｎａｎｊｉｎｇ　

２１００１８，Ｃｈｉｎａ；２．Ｋｅｙ　Ｌａｂｏｒａｔｏｒｙ　ｏｆ　Ｓｍａｒｔ　Ｇｒｉｄ　Ｔｅｃｈｎｏｌｏｇｙ　

ａｎｄ　Ｅｑｕｉｐｍｅｎｔ　ｉｎ　Ｊｉａｎｇｓｕ　Ｐｒｏｖｉｎｃｅ，Ｓｏｕｔｈｅａｓｔ　Ｕｎｉｖｅｒｓｉｔｙ，Ｎａｎｊｉｎｇ　

２１００１８，Ｃｈｉｎａ；３．Ｃｈｉｎａ　Ｅｌｅｃｔｒｉｃ　Ｐｏｗｅｒ　Ｒｅｓｅａｒｃｈ　Ｉｎｓｔｉｔｕｔｅ（Ｎａｎｊｉｎｇ），Ｎａｎｊｉｎｇ　

２１０００３，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｉｓ　ｔｈｅ　ｆｏｕｎｄａｔｉｏｎ　ｏｆ　ｓｏｌｖｉｎｇ　

ｐｒｏｂｌｅｍｓ　ｗｉｔｈ　ｕｎｃｅｒｔａｉｎｔｉｅｓ　ｉｎ　ｔｈｅ　ｐｏｗｅｒ　ｓｙｓｔｅｍ．Ｗｉｔｈ　ｔｈｅｄｅｖｅｌｏｐｍｅｎｔ　ｏｆ　ｉｎｔｅｒｍｉｔｔｅｎｔ　ｅｎｅｒｇｙ　ａｎｄ　ｔｈｅ　ｒａｉｓｅ　ｏｆ　ｒａｎｄｏｍｎｅｓｓ　ｏｆ　ｐｏｗｅｒ　ｓｙｓｔｅｍｓ，ｒｅｓｅａｒｃｈｅｓ　ｏｎ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｈａｖｅｅｘｔｅｎｓｉｖｅｌｙ　ｉｎｃｒｅａｓｅｄ　ｉｎ　ｒｅｃｅｎｔ　ｙｅａｒｓ．Ｔｈｉｓ　ｐａｐｅｒ　ｉｓ　ｉｎｔｅｎｄｅｄ　ｔｏ　ｍａｋｅ　ａ　ｂｒｉｅｆ　ｓｕｒｖｅｙ　

ｏｎ　ａｌｇｏｒｉｔｈｍｓ　ｆｏｒ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｂｙｐｒｏｃｅｅｄｉｎｇ　ｆｒｏｍ　ｔｈｅ　ｐｒｉｎｃｉｐｌｅ　ｏｆ　ｔｈｅ　ａｌｇｏｒｉｔｈｍｓ　ａｎｄ　ｔｈｅｉｒ　ｍｅｒｉｔｓ　ａｎｄ　ｄｅｍｅｒｉｔｓ．Ｆｉｒｓｔｌｙ，ｐｒｏｂｌｅｍｓ　ｓｔｕｄｉｅｄ　ｂｙ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄｆｌｏｗ　ａｒｅ　ｔｒｅａｔｅｄ　ｉｎｃｌｕｄｉｎｇ　ｔｈｅ　ｄｅｖｅｌｏｐｍｅｎｔ　ｏｆ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ，ｃａｌｃｕｌａｔｉｎｇ　ｍｏｄｅｌｓ　ａｎｄ　ｅｖａｌｕａｔｉｎｇ　ｉｎｄｉｃａｔｏｒｓ，ｆｏｌｌｏｗｅｄ　ｂｙａｐｐｌｉｃａｔｉｏｎｓ　ｏｆ　ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｉｎ　ｐｏｗｅｒ　ｓｙｓｔｅｍｓ．Ｔｈｅ　ｐｒｉｍａｒｙ　

ｃｏｖｅｒａｇｅ　ｉｓ　ｔｈｅ　ｒｅｖｉｅｗ　ｏｆ　ｄｉｆｆｅｒｅｎｔ　ａｌｇｏｒｉｔｈｍｓ　ｏｆｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ　ｃｌａｓｓｉｆｉｅｄ　ａｃｃｏｒｄｉｎｇ　ｔｏ　ｔｈｅ　ｄｉｆｆｅｒｅｎｃｅｓ　ｉｎ　ｔｈｅｏｒｙ．Ｔｈｅ　ｍａｉｎ　ｃｏｎｔｅｎｔｓ　ｉｎｖｏｌｖｅ　ｔｈｅ　ｐｒｉｎｃｉｐｌｅｓ，ａｄｖａｎｔａｇｅｓａｎｄ　ｄｒａｗｂａｃｋｓ　ａｎｄ　ａｐｐｌｉｃａｂｉｌｉｔｙ　ｏｆ　ｄｉｆｆｅｒｅｎｔ　ａｌｇｏｒｉｔｈｍｓ，ｆｏｌｌｏｗｅｄ　ｂｙ　ａｎ　ｏｖｅｒａｌｌ　ｅｖａｌｕａｔｉｏｎ　ａｎｄ　ｔｈｅｉｒ　ｄｅｖｅｌｏｐｍｅｎｔ　ｔｒｅｎｄ．Ｖｉｓｔａｓｏｆ　ｒｅｓｅａｒｃｈ　ａｒｅ　ｐｒｏｊ

ｅｃｔｅｄ　ｉｎｔｏ　ｔｈｅ　ｆｕｔｕｒｅ．Ｔｈｉｓ　ｗｏｒｋ　ｉｓ　ｓｕｐｐｏｒｔｅｄ　ｂｙ　Ｎａｔｉｏｎａｌ　Ｈｉｇｈ　Ｔｅｃｈｎｏｌｏｇｙ　

Ｒｅｓｅａｒｃｈ　ａｎｄ　Ｄｅｖｅｌｏｐｍｅｎｔ　Ｐｒｏｇｒａｍ　ｏｆ　Ｃｈｉｎａ（８６３Ｐｒｏｇｒａｍ）（Ｎｏ．２０１１ＡＡ０５Ａ１０５），Ｓｔａｔｅ　Ｇｒｉｄ　Ｃｏｒｐｏｒａｔｉｏｎ　ｏｆ　Ｃｈｉｎａ（Ｎｏ．ＤＺ７１－１３－０３６）ａｎｄ　Ｂｅｉｊｉｎｇ　

Ｎａｔｕｒａｌ　Ｓｃｉｅｎｃｅ　Ｆｏｕｎｄａｔｉｏｎ（Ｎｏ．３１３２０３５）．Ｋｅｙ　

ｗｏｒｄｓ：ｕｎｃｅｒｔａｉｎｔｙ；ｐｒｏｂａｂｉｌｉｓｔｉｃ　ｌｏａｄ　ｆｌｏｗ；ｃｏｒｒｅｌａｔｉｏｎ；ｓｉｍｕｌａｔｉｏｎ　ｏｆ　ｓａｍｐｌｉｎｇ；ａｐｐｒｏｘｉｍａｔｅ　ｃａｌｃｕｌａｔｉｏｎ櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧櫧（上接第７０页　ｃｏｎｔｉｎｕｅｄ　ｆｒｏｍ　ｐａｇ

ｅ　７０）Ｃａｌｃｕｌａｔｉｏｎ　ｏｆ　ＤＣ　Ｓｉｄｅ　Ｈａｒｍｏｎｉｃ　Ｃｕｒｒｅｎｔｓ　ｆｏｒ　ＬＣＣ－ＭＭＣ　Ｈｙｂｒｉｄ　ＨＶＤＣ　Ｔｒａｎｓｍｉｓｓｉｏｎ　Ｓｙ

ｓｔｅｍＺＨＡＮＧ　Ｚｈｅｒｅｎ，ＸＵ　Ｚｈｅｎｇ，ＸＵＥ　Ｙｉｎｇ

ｌｉｎ（Ｃｏｌｌｅｇｅ　ｏｆ　Ｅｌｅｃｔｒｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ，Ｚｈｅｊｉａｎｇ　

Ｕｎｉｖｅｒｓｉｔｙ，Ｈａｎｇｚｈｏｕ　３１００２７，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｆｏｒ　ａ　ｈｙｂｒｉｄ　ｈｉｇｈ－ｖｏｌｔａｇｅ　ｄｉｒｅｃｔ　ｃｕｒｒｅｎｔ（ＨＶＤＣ）ｔｒａｎｓｍｉｓｓｉｏｎ　ｓｙｓｔｅｍ　ｕｓｉｎｇ　

ｔｈｅ　ｌｉｎｅ　ｃｏｍｍｕｔａｔｅｄ　ｃｏｎｖｅｒｔｅｒ（ＬＣＣ）ａｔｔｈｅ　ｒｅｃｔｉｆｉｅｒ　ｓｉｄｅ　ａｎｄ　ｔｈｅ　ｍｏｄｕｌａｒ　ｍｕｌｔｉｌｅｖｅｌ　ｃｏｎｖｅｒｔｅｒ（ＭＭＣ）ａｔ　ｔｈｅ　ｉｎｖｅｒｔｅｒ　ｓｉｄｅ，ａ　ｃａｌｃｕｌａｔｉｏｎ　ｍｅｔｈｏｄ　ａｎｄ　ｔｈｅ　ｃｏｍｐｌｅｔｅｐｒｏｃｅｓｓ　ｏｆ　ｈａｒｍｏｎｉｃ　ｃｕｒｒｅｎｔ　ａｔ　ｔｈｅ　ＤＣ　ｓｉｄｅ　ａｒｅ　ｐｒｏｐｏｓｅｄ．Ｆｉｒｓｔｌｙ，ｔｈｅ　ＬＣＣ　ａｔ　ｔｈｅ　ｒｅｃｔｉｆｉｅｒ　ｓｉｄｅ　ｉｓ　ｒｅｐｌａｃｅｄ　ｂｙ　

ｔｈｅ　ｔｈｒｅｅ－ｐｕｌｓｅｈａｒｍｏｎｉｃ　ｖｏｌｔａｇｅ　ｓｏｕｒｃｅ，ｗｈｏｓｅ　ｈａｒｍｏｎｉｃ　ｏｕｔｐｕｔ　ｃｈａｒａｃｔｅｒｉｓｔｉｃ　ｉｓ　ｅｑｕｉｖａｌｅｎｔ　ｔｏ　ａ　ｔｗｅｌｖｅ－ｐｕｌｓｅ　ｃｏｎｖｅｒｔｅｒ．Ｓｅｃｏｎｄｌｙ，ｔｈｅ　ＭＭＣ　ｉｓｒｅｐｒｅｓｅｎｔｅｄ　ｂｙ　ａｐａｓｓｉｖｅ　ｅｌｅｃｔｒｉｃ　ｎｅｔｗｏｒｋ，ｗｉｔｈ　ａｎ　ｉｍｐｒｏｖｅｄ　ｃａｌｃｕｌａｔｉｎｇ　

ｍｅｔｈｏｄ　ｆｏｒ　ｔｈｅ　ｃｏｕｐｌｅｄ　ｌｉｎｅ　ｍｏｄｅｌ　ｉｎｔｒｏｄｕｃｅｄ　ｆｏｒｃａｌｃｕｌａｔｉｎｇ　ｔｈｅ　ａｄｍｉｔｔａｎｃｅ　ｍａｔｒｉｘ　ｏｆ　ＤＣ　ｔｒａｎｓｍｉｓｓｉｏｎ　ｌｉｎｅｓ．Ｆｉｎａｌｌｙ，ｂｙ　ｃｏｍｐａｒｉｎｇ　ｗｉｔｈ　ｈａｒｍｏｎｉｃ　ｃｕｒｒｅｎｔ　ｃａｌｃｕｌａｔｉｏｎ　ｒｅｓｕｌｔｓ　ａｔｔｈｅ　ＤＣ　ｓｉｄｅ　ｏｆ　ａ　ｄｉｇｉｔａｌ　ｓｉｍｕｌａｔｉｏｎ　ｍｏｄｅｌ　ｂａｓｅｄ　ｏｎ　ＰＳＣＡＤ／ＥＭＴＤＣ，ｔｈｅ　ａｃｃｕｒａｃｙ　

ｏｆ　ｔｈｅ　ｐｒｏｐｏｓｅｄ　ｍｅｔｈｏｄ　ｉｓ　ｖｅｒｉｆｉｅｄ．Ｔｈｉｓ　ｗｏｒｋ　ｉｓ　ｓｕｐｐｏｒｔｅｄ　ｂｙ　Ｎａｔｉｏｎａｌ　Ｈｉｇｈ　Ｔｅｃｈｎｏｌｏｇｙ　Ｒｅｓｅａｒｃｈ　ａｎｄ　Ｄｅｖｅｌｏｐｍｅｎｔ　Ｐｒｏｇｒａｍ　ｏｆ　Ｃｈｉｎａ（８６３Ｐｒｏｇｒａｍ）（Ｎｏ．２０１２ＡＡ０５１７０４）．

Ｋｅｙ　

ｗｏｒｄｓ：ｈｙｂｒｉｄ　ｈｉｇｈ－ｖｏｌｔａｇｅ　ｄｉｒｅｃｔ　ｃｕｒｒｅｎｔ　ｔｒａｎｓｍｉｓｓｉｏｎ；ＤＣ　ｌｉｎｅ；ｈａｒｍｏｎｉｃ　ｃｕｒｒｅｎｔ；ｃａｌｃｕｌａｔｉｏｎ　ｐｒｏｃｅｓｓ—

５３１—·综述·　刘　宇，等　电力系统概率潮流算法综述