将s级2s阶的辛Gauss方法用于电力系统暂态稳定性计算,提出了一种新的并行计算方法。该算法首先将微分一代数方程组经多级差分后转化为大规模非线性方程组,并利用牛顿法对其进行求解。在此基础上,利用矩阵分解方法将整体计算任务分解为两部分:一部分计算任务可按相应的级数或在不同的时间点上进行“解耦”,因而具有完全的时间并行性;对剩下的一部分计算任务,采用预处理GMRES方法对其进行空间并行求解,并为此提出了一种新的预处理方法。利用三个不同规模的算例系统,对所提算法的收敛性进行了测试,并在GPU上对算法进行了实际测试。测试结果表明,该算法可以获得很高的加速比,可以用于大规模电网暂态稳定性的实时分析计算。
This paper presents a novel parallel algorithm for power system transient stability computation. The proposed algorithm uses the s-stage 2s-order symplectic Gauss method to convert the differential-algebraic system simultaneously at s time points into a set of nonlinear algebraic equations, and the algebraic system is then solved using Newton's method. By the use of the matrix factorization technique, the solution of the linear equations involved in Newton's process is decomposed into two parts: the first part is decoupled according to the stage or at different time points, thus it is fully parallelizable-in-time, and the second part is solved using a preconditioned GMRES method while a new preconditioning way has been proposed. The convergence of the proposed algorithm has been tested on three example power systems. Furthermore, the proposed algorithm has been implemented on a single GPU based computer, and the results show the proposed algorithm achieves high speed-up ratio, and can be applied to the real-time computation of large-scale power system transient stability.