[发明专利]一种基于GPU的快速求解配电网潮流的计算方法

说明书

技术领域

本发明涉及一种电网潮流计算与分析的方法，具体来说涉及一种基于GPU的快速求解配电网潮流的计算方法。

背景技术

在实际的电力系统维修、检修或运行人员执行调度方案过程中对相关设备进行开关合闸等操作时必然引起联网设备监测仪表读数的变化，相应地在虚拟仿真系统中对此等操作事件的仿真也必须能够给予正确的模拟响应，从而反映出此等联网设备仪表读数值的正确变化，否则仿真意义将大打折扣。

“潮流计算”是实现此目标的关键核心，它的任务就是根据给定的运行条件和网络结构来确定整个系统各部分的运行状态，如各母线上的电压(幅值和相角)、网络中的功率分布以及功率损耗等。其计算结果是电力网络运行分析的基础，可为电力系统运行人员进行电力系统网络重构、故障处理、无功优化和状态估计提供参考依据，在电网调度、运行分析、操作模拟和设计规划中都发挥着重要的作用，是电力系统中最基本的、应用最为广泛的一种电气计算。

基于潮流计算的重大应用价值，众多国内外科研人员对其开展了深入研究。早期主要面向环状结构的高压输电网进行潮流计算分析。主要经历了高斯-赛德尔迭代算法、牛顿拉夫逊算法和国内外广泛使用的PQ分解算法三个研究发展阶段。

90年代开始面向结构更为复杂的配电网进行潮流计算研究。因配电网在结构上具有闭环结构、开环运行的特性，稳态运行时多呈辐射状，在发生故障或倒换负荷时，则呈短时环网运行结构；且在线路参数上R/X比值较大，网络的PQ节点多、PV节点较少等特点，使得传统的牛顿拉夫逊算法、快速解耦法出现病态收敛、不再有效。

为此，许多学者针对配电网结构特性一方面对牛顿法、快速解耦法进行了改进，但使之求解复杂化，牺牲了原有算法的收敛性和稳定性；另一方面陆续提出了前推回代类算法、回路阻抗法等新配电网潮流的计算方法。

其中前推回代类算法具有线性收敛、鲁棒性好、求解速度快、内存占用少等优异特性，适用于辐射状配电网的潮流计算。但在应用于环状或含有PV节点的配电网潮流计算分析时，需进行特殊处理从而增加了求解难度，通用性、收敛性不强。

S.K.Goswani等人提出的回路阻抗法，与上述算法相比，具有网孔处理能力强、收敛性好、稳定性高，适用范围广等优点。但在计算机求解过程中存在下述问题：(1)为利用稀疏矩阵技术简化回路阻抗矩阵的LU分解，王守相等人将配电网树状结构转化成标准二叉树形式进行节点编号优化，加大了算法复杂度，限定了中间节点状态的求解顺序，改变了配电网原有的网络结构不便于网络重构；王丹提出邻接表与二叉树结合的方式进行节点编号优化便于网络重构但也无法进行有环配电网阻抗矩阵的自动化生成。(2)CPU对阻抗矩阵进行LU分解时，刘耀年等人采用一定的技巧对对角线元素和相邻元素进行了优化，但因阻抗矩阵为对称满秩矩阵在计算过程中矩阵元素仍将占用大量内存，计算速度会随着节点数目的大量增加而大幅下降。

发明内容

本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种运算速度高、适用范围广的基于GPU的快速求解配电网潮流的计算方法。

本发明的目的可以通过以下技术方案来实现：

一种基于GPU的快速求解配电网潮流的计算方法，该方法包括以下步骤：

A、对配电网上的发电站节点以及各个负荷节点进行编号，形成配电网原始接图；

B、根据节点数目和节点间的连接关系构建十字链表；

C、通过电压互感器检测配电网中各个节点间的电压，电压互感器的二次侧经交流插件和A/D转换器将电压信号输出至GPU计算处理，根据配电网的原始接图以及配电网各个节点之间的电流流向关系，构建关联矩阵；

D、根据十字链表判断配电网中是否含有环网，若是，则破坏该环网处的闭环负荷节点，增加一个虚拟节点，相应地修改十字链表，并在关联矩阵中增加闭环负荷节点破坏后的虚拟节点编号的行与列，同时调整关联矩阵中的元素值，若否，则直接执行步骤E；

E、通过十字链表和关联矩阵，构建除了发电站之外的配电网中所有负荷节点的初始阻抗矩阵Z，通过GPU迭代求解配电网潮流，并将计算结果在相应仪表中显示。

所述的十字链表中包括顶点结点和边结点。