大规模电力系统暂态稳定并行计算研究

当前是我国电力系统发展的重要时期,跨区域的大规模互联电网正在逐步形成。为满足互联电网分析、规划和安全稳定运行对高性能暂态稳定仿真工具的迫切需要,本文对大规模电力系统暂态稳定并行仿真技术进行了研究,内容涉及暂态稳定串行算法改进、并行算法的构造、并行任务的划分及算法在高性能并行计算平台上的实现等方面。论文的主要工作如下： 1)提出了一种基于Shamanskii算法和非诚实牛顿法(Very Dishonest Newton Method, VDHN)的可变步长暂态稳定仿真组合算法。在基本的微分代数方程组联立求解框架下,根据隐式梯形积分局部截断误差理论,对步长进行控制,在保证精度的条件下,减少了积分步数；考虑牛顿类算法的收敛性,使用Shamanskii算法控制雅可比矩阵的更新,减少了不必要的更新计算；应用VDHN原理简化了迭代过程中电压向量的计算。算法可适应不同规模算例,在不同故障下都能较好提升计算速度,后文工作以此算法为基础展开。 2)提出了一种基于多核处理器的并行变步长VDHN算法用于暂态稳定仿真。引入α动态调度策略,将算法中最耗时的部分——发电机组的计算配置到多核CPU并行处理,并在仿真中动态调整各核心的计算量,以获得更好的负载平衡性能。进一步,在并行环境中考虑Newton类算法迭代时间的改变,自适应地调整雅可比矩阵更新策略,减少了并行计算过程中的串行部分和并行开销。算法复合加速比达到6.01倍,并可方便灵活地部署在多种软硬件平台上,适用性广。 3)提出了基于CPU-GPU (Graphics Processing Unit)异构平台的一种非诚实牛顿-稳定双共轭梯度(BiConjugate Gradient Stabilized Method, BiCGSTAB)暂态稳定并行算法。算法依据联立矩阵的双层对角加边结构,将整体计算分解为3部分：1.动态元件相关计算；2.子分区系统计算；3.联络系统计算。在异构平台上,第1、2部分被分配到多核CPU上进行处理。第3部分则采用可完全并行化的稳定双共轭梯度法在GPU上计算,并且为了减少迭代次数使用了稀疏近似逆预处理技术。并行任务的划分使用了超图算法,可以对电网进行更精细的描述,并用数学语言表示地理区域信息,显著提高了划分效果。算法可对万级节点电网进行超实时仿真,仿真时间仅为实际暂态过程的63.4%。 4)提出了基于CPU-GPU异构平台的一种多速率并行算法,应用于交-直流互联电网仿真。算法采用双层并行结构：第一层为“交-直并行”,考虑直流系统动态的独立性,将交流系统与直流系统解耦,分别部署在CPU和GPU上采用不同速率计算,系统间通过一定的接口时序交换数据；第二层为“直流系统时间并行”,直流系统在小步长下采用详细模型仿真,使用GPU实现了基于时间并行算法的流水线计算,可灵活设置流水线条数,对多个直流系统多积分时步并行求解。在此基础上,基于模型-视图-控制器模式构建了适应于大规模交直流互联电网暂态仿真的云计算原型系统,可方便地调用前文所述多种并行算法,并根据网络请求,分配合适的计算资源供用户使用。 本文实现的方法可用于电力系统的规划、分析及安全稳定控制等方面,进一步推动了电力系统暂态稳定并行计算技术的发展,为大规模互联电网暂态稳定并行仿真和新型的电力系统仿真工具提供新的研究思路。