[发明专利]一种基于稀疏技术的电力系统直角坐标PQ分解法潮流数据存贮方法

权利要求书

1.一种基于稀疏技术的电力系统直角坐标PQ分解法潮流数据存贮方法，其特征包括以下步骤：

步骤1：定义数据文件数组A(n,d)；

(1)定义与行号i对应的节点为父节点，与其相连的非零元素节点为子节点，零元素的子节点及参数不出现在存贮单元中；

(2)定义存放Y、B′、B″阵元素的三个虚拟数组为Y(n,3d₁)、B′(n-1,d₂)、B″(m,d₂)；Y(n,3d₁)数组中存放所有父节点和非零元素子节点的列号j及相应的参数，对应的最大非零元素数为S_1max；而B′(n-1,d₂)、B″(m,d₂)数组中均只存放所有父节点和非零元素子节点相应的参数，其非零元素的列号j均可由Y(n,3d₁)中得到，对应的最大非零元素数均为S_2max，且B″(m,d₂)中仅存放PQ节点对应的参数；

(3)将Y(n,3d₁)、B′(n-1,d₂)、B″(m,d₂)三个数组的数据共存于A(n,d)数组中；

(4)按各行中最大的非零元素数定义为最大存贮单元数U_max.new；

(5)按各行S_i实际统计的非零元素数定义为实际存贮单元数U_act.new；

(6)将A(n,d)数组分为“行号组”、“节点数组”、“Y阵组”、“B′阵组”、“B″阵组”，存放格式如下：

行号组i：为检验数据之用，存放与父节点对应的行号；

节点数组S_i1、S_i2：为高效读写数据之用，S_i1存放Y阵中各行父节点及非零元素子节点数之和，S_i2存放B′和/或B″阵中各行父节点及非零元素子节点数之和，其S_i1、S_i2值由程序自动累加；

Y阵组：存放与Y(n,3d₁)数组对应的Y阵数据，按列号递增顺序存贮父节点和非零元素子节点的列号j及相应的自导纳、互导纳的实部、虚部值；

B′阵组：存放与B′(n-1,d₂)数组对应的B′阵数据，按列号递增顺序存贮父节点和非零元素子节点相应的自导纳、互导纳的虚部值，其非零元素列号j由Y(n,3d₁)中得到；

B″阵组：存放与B″(m,d₂)数组对应的B″阵数据，按列号递增顺序存贮父节点和非零元素子节点相应的自导纳、互导纳的虚部值，其非零元素列号j由Y(n,3d₁)中得到；

步骤2：从数据文件读入全部支路数据；

步骤3：计算Y、B′、B″阵的所有元素并将数据存于A(n,d)数组中；

Y阵元素通常包括所有支路参数所形成的Y阵的实部和虚部元素，仅用于后续PQ分解法潮流程序中I_pi、I_qi或ΔP_i、ΔQ_i的计算；B′、B″阵元素用于后续程序中求解Δf_i、Δe_i，在计算B′阵元素的过程中去掉线路对地电容c、而在计算B″阵元素的过程中去掉线路电阻r；

步骤4：将A(n,d)数组的数据写入数据文件以备后续程序使用。