在传统电力系统分析中，负荷的波动、电网运行方式的变化和发电机的停运等因素造成了电力系统一定程度上的不确定性。随着电力工业的发展，以太阳能和风能等为代表的新能源接入电网，给电网带了明显的间歇性和随机性；微网、分布式电源和电动汽车等配电网新概念的发展，大大增强了电源、负荷与电网之间的互动性川，其直接结果导致了电力系统的不确定性显著增加，用于电力系统分析的概率潮流算法的研究日益重要。

1974年，Borkowaka 提出概率潮流计算方法叫，用以解决电力系统中诸多不确定因素。在随后四十年的时间里，概率潮流理论与方法得到了发展。与其几乎同时出现的随机潮流叫和概率潮流相互补充融合，逐渐形成处理电力系统不确定因素的体系：一般认为对于电力系统短期不确定因素采用随机潮流处理，而对于长期的具备规律性的不确定因素采用概率潮流处理，后者更趋向于概率分布的计算。

概率潮流计算的提出与发展，其最显著的意义是在进行电力系统分析时，考虑了系统各种不确定因素的随机性，从而使得计算分析更加贴合实际电网的运行状态。概率潮流的研究问题，主要集中在3个层面：系统模型、计算模型和计算方法。

系统元件的不确定性是引入概率潮流的根本原因，主要体现在发电机、负荷、输电线路和变压器的随机性。近些年，随着可再生能源并网规模的日益提高和电力用户的市场行为日趋突出，发电机和负荷的功率模型越发复杂。

就概率潮流计算模型而言，以四大类模型为主。Borkowaka基于简化的直流模型团提出了概率潮流计算方法。为了提高潮流计算的精度，Allan 分别在1976年和1981年提出了线性化交流模型和分段线性化交流模型，Sokierajski在1978年提出保留非线性的交流模型 。目前概率潮流的计算方法都是基于这4种模型进行。

对概率潮流算法的研究是概率潮流分析中的热点，具备广阔的研究空间与研究意义。一种性能良好的概率潮流计算方法应满足以下指标：能够求出输出随机变量的数字特征（包括均值和方差）及概率分布；能够处理多个随机变量间的相关性；满足实用化要求，结果具有足够精度的情况下尽量减少计算时间；满足通用性要求，对输入变量的数学模型不应有太高要求。这4项指标构成概率潮流算法研究的重点和难点，专家学者们从一方面或多方面入手展开研究，形成了当前的多种概率潮流计算方法。就目前已有研究成果而言，概率潮流算法可以大致分为模拟法、近似法和解析法3类。解析法结合电力系统复杂卷积计算的简化特性，通过对卷积计算的特殊处理衍生而来，广义上也可归入近似法。纵观现有概率潮流算法，存在的突出问题是难以将计算准确度和时效性有机统一。

以模拟法、近似法和解析法为分类基准，综述各类方法中具体算法的原理步骤、优缺点评估以及发展趋势。以算法中现有的问题及不足为依托，结合当前电力系统的研究热点与新兴技术，对概率潮流的发展方向进行展望。

潮流计算是电力系统分析的基础，其典型计算潮流方程如下所示：

式中：W为系统节点功率注入向量；X为节点电压向量；Y为系统网络参数；Z为系统支路潮流向量。概率潮流的计算正是在上式的基础上，通过考虑输入变量W和Y的概率特性，获得系统状态变量X和Z的分布情况，从而全面地给出系统的运行状况和概率特征。

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

1）电力系统规划，包括电源规划、电网规划和无功规划等规划问题 。规划问题的求解均以潮流计算作为基础，在不确定性增加的电力系统中，概率潮流将成为含随机因素规划问题的求解前提。

2）静态安全分析，作为电力系统分析的基本问题，采用概率潮流的静态安全分析方法可以更加真实地反映电力系统全而信息 。

3）电力系统运行状态实时在线分析，包括机组组合、在线调度、电力市场机制下的源一网嗬互动。应用某些概率潮流算法，可以在不显著增加计算次数与时间的条件下，更为精确地分析电力系统的运行状态及变化趋势，为电力系统实时分析提供了强有力的工具 。

4）其他基于潮流运算的系统分析，例如最优化潮流计算、电力系统风险评估、互动型配电网潮流计算等。

模拟法概率潮流，是将电力系统中的不确定因素作为随机变量建立概率模型，然后抽取概率分布的样本，最后统计输出变量的分布特征。传统的模拟法概率潮流计算方法一般是指随机采样的蒙特卡洛模拟法 ，后来基于随机模拟法改进衍生出重要抽样法 、拉丁超立方采样法和拟蒙特卡洛方法等。

电力系统概率潮流算法随机采样的蒙特卡洛模拟法

蒙特卡洛模拟是二战时期美国物理学家Metropoli 在执行曼哈顿计划的过程中提出的。蒙特卡洛模拟法以随机模拟和统计实验为手段，是一种从随机变量的概率分布中，通过随机选择数字的方法产生一种符合该随机变量概率分布特性的随机数值序列，作为输入变量序列进行特定分析的求解方法。其计算关键与核心步骤如下：①对潮流方程的输入变量W构造相应的概率模型；②产生随机数序列，作为系统的抽样输入进行大量的数字模拟，每一组采样值通过潮流计算得到相应的模拟实验值；③系统计算，对模拟实验结果进行统计处理，给出所求问题的解。

蒙特卡洛模拟的优点在于样本数量足够大时，计算结果足够精确；并且计算量一般不受系统规模的影响，该方法的抽样次数与抽样精度的平方成反比。缺点在于为提高计算精度，往往需要提高系统抽样规模，从而导致计算时长过大。考虑其精度优势，随机采样的蒙特卡洛模拟法一般用来作为基准方法进行比较，是衡量其他方法准确性的重要参考。

电力系统概率潮流算法重要抽样法

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

如何选取新分布中系统的概率分布使得随机变量在期望不变的情况下减小方差是重要抽样法的关键步骤。有文献采用迭代法搜索重要分布函数，给出了若干重要分布函数的定义方法，并结合分散抽样的技巧提高重要抽样法的收敛速度。也有文献利用蒙特卡洛方法模拟出负荷样本，然后利用核密度估计方法估计出负荷模型的密度函数，将之作为重要抽样密度函数，计算出支路潮流和节点电压的概率密度函数 。

重要抽样法在电力系统的概率估计中有着广泛应用，该方法可以快速准确地计算出系统运行状态的期望值，为系统分析提供参考。但重要抽样法中仅以期望为研究对象，对于概率变量的方差、概率分布等参数分析存在天然缺陷，计算结果局限性较大。

电力系统概率潮流算法拉丁超立方采样法

为了避免随机采样的蒙特卡洛模拟法的大规模抽样，Mckay等人于1979年提出了拉丁超立方采样法 。它是一种分层采样法，通过改进输入随机变量的样本生成过程，保证其采样值能够有效地反映随机变量的整体分布，算法的出发点就是确保所有的采样区域都能够被采样点覆盖。其基本运算过程分为如下两个步骤：采样和排列。

拉丁超立方采样法的不足是对输入随机变量的处理较为复杂，一方面要求已知输入随机变量的概率分布函数或累积分布函数，另一方面对不同类型概率分布的随机变量相关性需要特殊变换处理困。但该方法作为一种非常有效的估计输出随机变量期望值的方法，由于采样值能够确保覆盖所有输入随机变量的整个分布区域，无须大规模抽样，并且可以有效处理输入变量之间的相关性和随机性，在准确性、稳健性和时效性上都有较大的优势。

电力系统概率潮流算法拟蒙特卡洛法

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

低差异序列，又称伪随机数列，是一系列数值确定的[0,1]区间中的数。在d维变量的空间中，低差异序列中己有n-1个数，生成第n个数的方法是：将这个数插入已有数列中最大的“空白”处，即避免数列在局部空间聚集，从而保证了有限数据的空间全覆盖。

目前拟蒙特卡洛法己经被应用于概率最优潮流计算和含互动式新能源的电网静态稳定分析中。由于采样过程中一次性生成所需序列，该方法具有比拉丁超立方采样法更高的计算效率。但拟蒙特卡洛法对多变量的高维度问题理论基础薄弱、计算效果差，因此目前多用于小规模电力系统分析计算。

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

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

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

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

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

近似法是利用输入随机变量的数字特征近似描述系统状态变量统计特性的方法。该方法避开了大规模的重复抽样，因而求解速度较快，又因其能够计及系统输入变量之间的互相关性，因而受到重视。目前研究应用较多的有点估计法、一次二阶矩法和状态变换法 。

电力系统概率潮流算法点估计法

点估计法是一种概率统计方法，目前所做的应用研究都是基于1998年Hong在已知输入随机变量的连续分布下提出的点估计法。该方法能够根据已知随机变量的概率分布，求得待求随机变量的各阶矩。

点估计法属于逼近技术的一种，利用输入随机变量的统计信息来逼近输出随机变量的数字特征。其主要运算过程分为以下几步。

1）用潮流方程中输入随机变量W的各个分布函数求出相应的前2M-1阶中心矩。

2）通过构造的方式，利用前2M-1阶中心矩独立求出每个输入随机变量的M个离散状态，使得这M个离散状态包含了前2M-1阶中心矩的所有信息。

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

4）用求得的潮流方程输出随机变量X和Z的M´K个离散状态逼近相应的期望值与方差等相关数字特征。

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

1）该方法中实际的输入量为输入随机变量前2M-1阶中心矩，此中心矩可以由概率分布函数直接求出，也可以由大量样本逼近拟合方程式展开得到，这样就不必受限于必须已知输入变量概率分布的条件约束。

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

3）输出随机变量有2M-1阶多项式逼近的精度，为了提高估计的精度，可以增加输入变量的高阶矩信息，即增加取点个数。但实际应用中点个数M大于3时不仅急剧增大计算量，而且往往造成解的结果非实数，因此M通常取2或3，即构成常用的两点估计法和三点估计法。

两点估计法计算简单、容易实现，但其只利用输入变量的前三阶矩信息，计算精度低；三点估计法既能得到较高精度的估计值，又保持了简易性，在点估计法中广为使用。点估计法的缺点在于计算结果中随机变量的高阶矩不够精确，无法准确获得变量的概率分布函数。同时在处理输入变量的时间和空间相关性上具有一定的计算复杂度。

电力系统概率潮流算法一次二阶矩法

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

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

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

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

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

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

电力系统概率潮流算法状态变换法

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

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

电力系统概率潮流算法近似法概率潮流展望

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

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

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

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

电力系统概率潮流算法快速傅里叶变换

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

电力系统概率潮流算法半不变量法

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

半不变量法具备的最大优势在于计算方法简单、计算效率高，虽然计算结果精确度存在一定争议，但满足工程应用要求，因此受到广泛研究。关于半不变量法的算法改进研究主要集中在3个问题上：一是如何处理输入随机变量之间的相关性；二是如何正确分析静态安全稳定问题；三是如何更为精确地描述系统运行状态。有文献提出一种基于Cholesky分解的计及输入变量相关性的半不变量法概率潮流计算方法，并提出基于蒙特卡洛抽样的方法解决一些输入变量的半不变量难以被常规数值方法求解的问题。

电力系统概率潮流算法序列运算理论

为了将原本复杂的卷积计算转换为简单的算术运算，康重庆教授创建了序列运算理论，在此基础上扩展了概率性序列运算方法，从而形成了全新的电力系统不确定性分析框架。

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

电力系统概率潮流算法解析法概率潮流展望

解析法概率潮流的研究重点在卷积运算的高效处理上，快速傅里叶运算难以处理大规模电力系统的多变量计算，序列运算理论体系架构的特殊性使之难以在短时间内大规模推广，因此半不变量法是解析法中研究的热点算法。半不变量法在处理正态分布的变量相关性上己有较多的研究，该算法的研究重点是非正态分布的变量相关性处理，改善各阶矩和半不变量的计算精度与计算效率，获得更为准确的概率密度函数。

电力系统中不确定因素的增加使得概率潮流成为研究热点之一。从概率潮流算法原理的角度出发，将电力系统概率潮流算法分为模拟法、近似法和解析法3类，对各算法的研究现状和发展进行了综述和总结，并对各类算法的发展趋势提供了分析思路。

概率潮流与确定性潮流本质上的区别在于输入变量的不确定性，计算目标是合理反映随机分布本身的概率特性对于系统潮流的影响，因此，概率潮流的研究目的可概括为随机变量的概率输入及其对潮流输出的概率贡献。模拟法分析输入随机变量的抽样方法，近似法关注输出与输入之间的转化关系，解析法重点讨论系统方程卷积的计算与简化，而3类方法中随机变量概率分布的表征和相关性分析都是关键因素。

而对新环境下电力系统的全新特征，未来概率潮流算法研究的总体发展趋势应该包含以下几个主要方向：可再生能源功率输出的概率密度函数并不一定精确满足特定的典型分布，如何生成适用于概率潮流计算的准确的可再生能源功率输出的概率分布非常重要，目前已有的核密度估计法。具有很大的研究空间；在大规模新能源集中并网的趋势下，如何处理区域内新能源功率输出的复杂相关性变得日益重要；电动汽车作为电网互动性负荷，与电网的功率交互是双向的，如何对双向功率负荷进行概率建模也是研究难点之一；电力市场机制下的电源和负荷具备开放性特点，基于市场特征和需求响应的最优概率潮流分析具备很高的研究意义。2100433B

遗传算法是80年代出现的新型优化算法，近年来迅速发展，它的机理源于自然界中生物进化的选择和遗传，通过选择、杂交、变异得核心操作，实现“优胜劣汰”。它的主要特点是:可从多初值点开始，沿多路径搜索实现全局或准全局最优;可方便地处理混合整数离散J陛问题;是一种有效的自适应优化方法。

GA应用于潮流优化问题时，一般步骤为:首先随机给出一组初始潮流解，受各种约束条件约束，然后通过目标函数评价其优劣，然对其编码，通过遗传操作—选择、杂交和变异，使其重新组合，评价值低的被抛弃，只有评价值高的有机会将其特征迭代至下一轮解，最后这码串对应的解将趋向优化 。

遗传算法优点是具有很好的全局寻优能力，优化结果普遍比传统优化方法好。缺点是计算量比较大，计算时间长。现在遗传算法的研究主要集中在以下两方面:通过改进目标函数计算方法以提高其计算速度，通过改进遗传算法的操作改进整体收敛J陛和寻优性能。

在遗传算法操作研究方面，在一个103节点系统上研究了使用不同的算子参数对迭代次数和优化结果的影响，还研究了控制变量约束的影响，建议在寻优过程中不断缩小解空间。研究了多种用于提高GA效率及精度的方法，表明同时变罚因子及变权重因子的GA应用于经济调度中最有效，它最能保证收敛精度，虽然它牺牲了一些收敛时间。

针对目标函数计算加速，也就是潮流计算加速，将潮流方程中PV节点转为ve节点作为控制变量，同时将网络按节点联系进行分层，以形成一个带状稀疏阵，然后针对网络分层的特点使用一种高效的改进高斯消去法求解线性方程组。在lEEE57, 118, 300和KT896, EvI试验网络上的对比计算表明其方法的速度比浅解耦潮流算法要快。

此外，GA还用于解决含电力电子设备的灵活交流输电系统这样的非凸性的优化。对此进行了研究，结果表明遗传算法在这种非线性、非光滑、不可微的函数优化上十分适合。

利用电子计算机进行潮流计算从20世纪50年代中期就已经开始。此后，潮流计算曾采用了各种不同的方法，这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几点：

（1）算法的可靠性或收敛性

（2）计算速度和内存占用量

（3）计算的方便性和灵活性

电力系统潮流计算属于稳态分析范畴，不涉及系统元件的动态特性和过渡过程。因此其数学模型不包含微分方程，是一组高阶非线性方程。非线性代数方程组的解法离不开迭代，因此，潮流计算方法首先要求它是能可靠的收敛，并给出正确答案。随着电力系统规模的不断扩大，潮流问题的方程式阶数越来越高，已达到几千阶甚至上万阶，对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的。这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法。

在用数字计算机求解电力系统潮流问题的开始阶段，人们普遍采用以节点导纳矩阵为基础的高斯-赛德尔迭代法（一下简称导纳法）。这个方法的原理比较简单，要求的数字计算机的内存量也比较小，适应当时的电子数字计算机制作水平和电力系统理论水平，于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法（以下简称阻抗法）。

20世纪60年代初，数字计算机已经发展到第二代，计算机的内存和计算速度发生了很大的飞跃，从而为阻抗法的采用创造了条件。阻抗矩阵是满矩阵，阻抗法要求计算机储存表征系统接线和参数的阻抗矩阵。这就需要较大的内存量。而且阻抗法每迭代一次都要求顺次取阻抗矩阵中的每一个元素进行计算，因此，每次迭代的计算量很大。

阻抗法改善了电力系统潮流计算问题的收敛性，解决了导纳法无法解决的一些系统的潮流计算，在当时获得了广泛的应用，曾为我国电力系统设计、运行和研究做出了很大的贡献。但是，阻抗法的主要缺点就是占用计算机的内存很大，每次迭代的计算量很大。当系统不断扩大时，这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺点，后来发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统，在计算机内只需存储各个地区系统的阻抗矩阵及它们之间的联络线的阻抗，这样不仅大幅度的节省了内存容量，同时也提高了计算速度。

克服阻抗法缺点的另一途径是采用牛顿-拉夫逊法（以下简称牛顿法）。牛顿法是数学中求解非线性方程式的典型方法，有较好的收敛性。解决电力系统潮流计算问题是以导纳矩阵为基础的，因此，只要在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿潮流程序的计算效率。自从20世纪60年代中期采用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、计算速度方面都超过了阻抗法，成为后来仍被广泛采用的方法。

在牛顿法的基础上，根据电力系统的特点，抓住主要矛盾，对纯数学的牛顿法进行了改造，得到了P-Q分解法。P-Q分解法在计算速度方面有显著的提高，迅速得到了推广。

牛顿法的特点是将非线性方程线性化。20世纪70年代后期，有人提出采用更精确的模型，即将泰勒级数的高阶项也包括进来，希望以此提高算法的性能，这便产生了保留非线性的潮流算法。另外，为了解决病态潮流计算，出现了将潮流计算表示为一个无约束非线性规划问题的模型，即非线性规划潮流算法。

近20多年来，潮流算法的研究仍然非常活跃，但是大多数研究都是围绕改进牛顿法和P-Q分解法进行的。此外，随着人工智能理论的发展，遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。但是，这些新的模型和算法还不能取代牛顿法和P-Q分解法的地位。由于电力系统规模的不断扩大，对计算速度的要求不断提高，计算机的并行计算技术也将在潮流计算中得到广泛的应用，成为重要的研究领域。

风电出力受自然天气条件的影响很大，当系统中含有大规模风电场时，确定性潮流无法全而反映系统稳态运行情况，由前述随机潮流的概念可知，系统随机潮流分析是解决该问题的有效方法。国内外学者对含风电场的电力系统随机潮流研究早已有涉及。从现有研究成果来看，系统随机潮流模型算法与常规系统基本一致。事实上，风电这一非常规随机因素引入后，进行系统随机潮流分析需要解决的关键问题是模型中如何计及风电场出力随机性。此外，由于风电出力随机性区别于常规系统中一般随机性因素而具有一定特殊性，对此计算中如何处理亦是含风电场随机潮流分析需要进行深入研究的问题。现主要从上述两个方而综述含风电场的随机潮流方法，具体如下。

电力系统随机潮流及其安全评估模型中风电出力随机性处理

蒙特卡罗抽样法是进行风电场出力随机性分析最直接简单的方法，如有文献均采用蒙特卡罗法分析风电场出力的随机性，并采用模拟法计算含风电场的电力系统随机潮流，比较分析了风电接入对系统节点电压和支路潮流的影响。

为采用计算速度相对较快的解析法计算含风电场随机潮流，有文献分别提出了风电出力随机性不同处理方法 。

电力系统随机潮流及其安全评估含风电场随机潮流计算几个特殊问题

(1)随机潮流模型中风电场潮流计算模型问题。现有含风电场随机潮流研究时，大多数文献假设风场与外界系统按恒功率因数发生无功功率交换，即将风电场等效为PQ节点，该处理方式与系统实际并不完全相符，为此有文献将风场风机分为恒速}巨频异步发电机、恒功率因数控制方式的变速恒频双馈机、恒压控制方式的变速恒频双馈机三种类型，分别建立对应无功概率模型，采用半不变量和Gram-Charlier级数展开法研究了风电场不同控制方式下系统的随机潮流 。

(2)风场间出力相关性问题。随机潮流计算中往往为了简化模型、加快计算速度而假设系统各节点的注入功率相互独立。事实上，这一假设与常规系统实际运行已经存在一定差异，而风电并网将进一步加大该差异。因为多个风电场地理位置相对靠近而使得其基本处于同一风速带，从而使得各风场风速具有较强的相关性，各风场出力也具有较强的相关性。因此，随机潮流分析时应考虑风速相关性以更准确地评估风场并网后系统静态运行特性。不少学者对此进行了探讨。

(3)风场内有功无功出力相关性问题。有文献提出含风电场的电力系统随机潮流分析时，应计及风电场有功、无功相关性尤其是采用异步发电机的风电场，该文将风电场有功无功同其灵敏度矩阵中对应的权重元素组合为一个随机变量，采用半不变量法进行随机潮流分析，研究表明计及风场有功无功相关性后计算结果更准确，更接近系统电压或潮流的真实分布情况。

电力系统随机潮流及其安全评估含风电场的电力系统随机潮流研究展望

由上述分析可知，目前含风电场的电力系统随机潮流基本模型算法已建立。但是，当前含风电系统随机潮流研究大多均基于风速服从Weibull分布的基本假设。事实上，含风电系统的随机潮流计算应根据分析周期的不同分为以下两种:中长期随机潮流问题;短期随机潮流问题。其中，中长期随机潮流计算时一般假定风速服从Weibull分布，考虑风电出力从零到满发随机波动，其分析结果适用于系统长期或中长期评估。而短期随机潮流计算时一般考虑风电出力值为预测结果加上随机波动的预测误差，即系统随机变量为风功率预测误差，此时的结果适用于系统短期或在线评估分析。基于上述随机潮流问题的划分，给出几个值得深入研究的相关问题，以供参考。具体如下:

(1)含大规模风电并网的中长期随机潮流解析化方法研究。现有含风电场的随机潮流解析方法，虽然其实用性较强，但对于风电并网规模较大且需考虑风电出力波动范围较大的中长期分析这一特殊场合，可能存在不适之处，具体分别表现模型与算法两个方而。首先，现有解析方法中大多采用最常用的线性化交流模型，其精度能否适应上述特殊场合有待深入研究。其次，在算法方而，现有解析算法以基于半不变量和级数展开法最为常用，其中级数大多为A型Gram-Charlier级数或Edgeworth级数。而数学界研究已表明，当随机变量的三阶或四阶矩超出一定范围时，其逼近所得概率密度函数可能出现负值，导致结果不满足基本概率公理。因此，将此方法直接用于上述场合分析值得探究。

(2)计及动态调度策略的中长期随机潮流方法研究。在现有随机潮流文献中，亦即中长期随机潮流研究，为简化分析或受模型、方法限制，负荷或间歇性电源出力波动时均没有考虑系统中各机组出力调整，而是将其视为不变的恒量，完全依靠单一的平衡节点实现功率平衡，这显然不相符合电网实际运行。事实上，电网运行时会根据负荷水平的不同及时调整系统发电计划、重新安排机组出力（尤其是承担峰荷的机组），以力求系统运行经济性。此外，一方而随着智能电网建设工作的深入开展，系统的互动性、兼容性、经济性等特征愈发明显。这在一定程度上也预示着未来电网分析工作必须而对诸多随机因素并恰当处理系统在此环境下运行状态的动态调整问题，从而使系统分析方法满足实际需求。另一方而，我国风电“建设大基地、融入大电网”的发展路线，也进一步突显了计及动态调度策略是系统分析的必然要求。因此，研究计及动态调度策略的中长期随机潮流分析方法具有重要意义，应作为后续研究的重要方向。

(3)短期随机潮流方法研究。由于短期随机潮流研究时一般考虑功率波动范围相对较小，风电出力随机变量取决于风功率预测误差，因而当前含风电场的电力系统随机潮流计算模型与求解方法基本适用。然而，难点问题是如何建立准确的风电功率预测误差分布模型，用于描述风电出力短期随机波动性。因此，研究一种考虑风电功率预测误差分布的随机潮流实用方法是短期随机潮流研究的重点。