基于网络分割电力系统潮流并行算法研究.docVIP

基于网络分割电力系统潮流并行算法研究 　　【摘要】本文介绍了一种基于网络分割的电力系统潮流并行算法，按照支路切割有效地将大规模电力系统的潮流计算问题转化为若干个子网络的潮流计算问题，通过分配于各个PC机的子网络之间的分解计算获取加速比。该算法在IEEE 14节点、IEEE 30节点和IEEE 118节点网络上的验算结果表明：此算法具有较高的收敛速度和计算精度，应用该算法获得的潮流计算结果正确，适合在网络计算环境中实现。 　　【关键词】PC机群；潮流计算；并行算法；支路切割；网络分块 　　 　　 1.引言 　　 随着我国电力工业的发展，我国的电网已经形成大规模的超高压、大机组、大电网的全国性互联电网，根据互联电网安全稳定运行的需要，对全网一体化仿真计算已经成为必要。有着巨大的存储容量和较快计算速度的PC计算机群在电网计算中得到广泛的应用。 　　 2.网络分块的必要性和方法[1] 　　 2.1 网络分块的必要性 　　 潮流计算的计算量与电力系统的节点数目有关，一个具有N个独立节点的网络，其节点导纳矩阵和节点阻抗矩阵元素的数目为N2N，且矩阵各元素为复数，在应用牛顿一拉夫逊法进行潮流的迭代计算时，其雅可比矩阵的阶数为2N-2。对于大规模电力系统而言，其迭代计算量非常庞大，即使是PQ分解法，在忽略了电压幅值对有功功率的影响以及电压相位对无功功率的影响并将雅可比矩阵常数化后，其迭代计算依然要消耗可观的机时。 　　 2.2 网络分块的方法 　　 上世纪50年代初，Kron提出了网络分块算法理论，网络分块方法之一是支路切割法，即在子网络之间的联络线处加入电流源，以代替各子网络之间状态变量的耦合影响，电流源矢量作为该方法的协调变量。因为把电流源矢量作为协调变量，所以应把电力系统中大部分的量测转换为电流量测。网络的物理分块方法和并行算法的基本意义都在于降低状态估计中雅可比矩阵的阶数以缩短计算时间。本文采用的是基于支路的网络数学分块方法。 　　 3.潮流计算数学模型[2][3] 　　 电力系统潮流计算问题经过线性化处理后，通常表示为求解如下一组线性代数方程：，其中为非奇异对称稀疏矩阵，向量为待求???变量，向量为给定的独立变量。该式的求解具有多种解法，目前在电力系统中使用得最有效的串行解法是三角分解法，通过前代和回代过程求解，即： 　　 ；； 　　 要进行并行求解，就要使式三个方程能够同时计算。 　　 稀疏矢处理量直接法和逆矩阵方法是潮流稀疏线性代数方程组并行计算的主流。它们各有其优势和不足：稀疏矢量法，能够保持因子表的稀疏性，但操作间依赖关系较强；逆矩阵法并行性较好，但逆矩阵的使用大大增加了矩阵中的非零量，增加了计算量。因此，要想提高潮流问题求解的并行效率还在十把稀疏矢量直接法和逆矩阵方法结合起来。 　　 4.网络协调与分解[4][5] 　　 基于BBDF模型的潮流网络分解与协调计算分解如下图1所示。 　　 　　图1 BBDF模型 　　Figure 1 Model of BBDF 　　 BBDF算法是电力系统潮流并行算法中比较成功的一种实现方法。其任务划分如图1所示，其中A1，A2...An之间相互独立，它们只与An十1有交互关系，而且图中每一分块还可进一步进行BBDF划分。 　　 潮流计算基本矩阵可分解为：，、为上三角和下三角矩阵。 　　 其中： 　　 对A的并行分解如下： 　　 　　 L矩阵的并行分解如下： 　　 　　 　　 最后计算： 　　 U矩阵的并行求解如下： 　　 　　 　　 　　 先计算， 　　 BBDF划分方法结合了电力系统物理意义，是一种非常典型的划分方法，将BBDF应用到在网络计算平台上执行，也是一种非常好的划分方案。 　　 5.网络计算环境下数学模型[6][7][8] 　　 网络计算也是一种广域分布式计算，象并行计算一样，也要进行任务的划分，任务的划分也一件比较困难的事情，到目前为止，仍然没有一个比较容易实现的非常成熟的方法来进行任务划分。并行计算领域内有一些比较成熟的区域划分算法，如BBDF算法等，但这些方法要么是针对应特殊应用背景，要么是根据特殊的数学变换提出来的，其应用都有一定局限性。任务的划分是进行网络计算的第一步，任务划分如果很复杂，将会直接影响网络计算的应用。本文对潮流计算采取了基于电流注入法，下面看一个简单的例子。 　　 　　图2 互联电力系统 　　Figuer 2 Interconnected power system 　　 　　图3 互联电力系统的分解 　　Figure3 Decomposition of interconnected power system 　　 通常潮流方程可表示如下： 　　 　　 其中：