[发明专利]一种电力系统潮流计算的支路功率矩阵计算方法

权利要求书

1.一种电力系统潮流计算的支路功率矩阵计算方法，其特征在于：包括以下步骤：

E1、读支路首节点号数组I、末节点号数组J、电阻数组R、电抗数组X、对地电纳数组B、变压器变比数组K；

所述的首节点号数组I、末节点号数组J、电阻数组R、电抗数组X、对地电纳数组B、变压器变比数组K分别按顺序存放所有支路的首节点号i_m、末节点号j_m、电阻r_m、电抗x_m、对地电纳b_m、变压器变比k_m，其中变压器支路对地电纳b_m为0，变压器非标准变比k_m侧的节点号加个负号，输电线路变比k_m设为1，其中下标m为支路序号；

E2、形成支路首节点电压相量数组UI和支路末节点电压相量数组UJ；

根据首节点号数组I和已计算出的节点电压相量数组U形成支路首节点电压相量数组UI，根据末节点号数组J和已计算出的节点电压相量数组U形成支路末节点电压相量数组UJ；

式中，abs为Matlab求绝对值函数；

E3、令Z＝R+j X；

E4、根据支路首节点号i_m为负的标志形成变比k_m位于首节点侧的变压器支路号数组Ik，根据支路末节点号j_m为负的标志形成变比k_m位于末节点侧的变压器支路号数组Jk；

形成数组Ik和数组Jk使用Matlab的find函数实现：

Ik＝find(I<0) (2)

Jk＝find(J<0) (3)

E5、按输电线路类型计算支路功率；

由于支路中输电线路比例较大，先按输电线路类型计算所有支路的功率；输电线路支路功率计算公式为：

式中，SIJ为从支路首节点i流向末节点j的复功率数组，SJI为从支路末节点j流向首节点i的复功率数组，节点号i和节点号j为带负号标志的节点号i_m和节点号j_m去掉负号标志的节点号，“.*”表示两数组对应元素相乘，“./”表示两数组对应元素相除，conj为Matlab求复数共轭函数；

E6、计算变比k_m位于首节点i侧的变压器支路的支路功率；

式中，SIJ(Ik)表示以数组Ik所有元素为下标提取数组SIJ的元素，SJI(Ik)表示以数组Ik所有元素为下标提取数组SJI的元素，UI(Ik)表示以数组Ik所有元素为下标提取数组UI的元素，UJ(Ik)表示以数组Ik所有元素为下标提取数组UJ的元素，K(Ik)表示以数组Ik所有元素为下标提取数组K的元素，Z(Ik)表示以数组Ik所有元素为下标提取数组Z的元素；

E7、计算变比k_m位于末节点j侧的变压器支路的支路功率；

式中，SIJ(Jk)表示以数组Jk所有元素为下标提取数组SIJ的元素，SJI(Jk)表示以数组Jk所有元素为下标提取数组SJI的元素，UI(Jk)表示以数组Jk所有元素为下标提取数组UI的元素，UJ(Jk)表示以数组Jk所有元素为下标提取数组UJ的元素，K(Jk)表示以数组Jk所有元素为下标提取数组K的元素，Z(Jk)表示以数组Jk所有元素为下标提取数组Z的元素。