[0049] 本发明实施例提供的一种考虑PI节点的HELM潮流计算方法,包括以下步骤:
[0050] 步骤一,建立包含PQ节点和PI节点的电网的数学模型。
[0051] 根节点为平衡节点,另外有PQ节点和PI节点:
[0052] 若节点为PQ节点,由节点功率方程:
[0053]
[0054] 式中m为PQ节点集合,i、k为节点号;Vi、Vk分别为节点i和节点k的电压;Yik为节点i、节点k之间的导纳;Si为节点i的注入功率;N为节点总数;上标*表示共轭关系;
[0055] 若节点为PI节点,由节点功率方程:
[0056]
[0057]
[0058] 式中p为PI节点的集合,Pi表示i节点的注入有功功率,Ii表示i节点的注入电流,|Ii|表示其绝对值, 表示i节点的电流幅值,Vi、Vk分别为节点i和节点k的电压;Yik为节点i、节点k之间的导纳,上标*表示共轭关系;N为电网节点总数;Re(·)表示复数向量的实部。
[0059] 步骤二、利用HELM方法构造全纯函数数学模型。
[0060] 由于节点电压与节点注入功率等有关系,而注入功率包括有功和无功,常用复数表示,而全纯函数是复数分析方法,因此根据内嵌全纯函数方法,可以构造一个内嵌纯虚函数:
[0061]
[0062] 式中s为内嵌参数因子;Vi(s)为i节点的s级数展开;n为s级数阶数;Vi[n]为i节点n的电压的s级数n阶项系数。s表示频域算子s的n次项。
[0063] 若节点为PQ节点,将式(4)代入式(1)并构造相应的全纯函数可得:
[0064]
[0065] 式中: Yii为节点导纳矩阵中i节点的自导纳;Yik为节点导纳矩阵中i和k节点间的互导纳;Yi,shunt为i节点的接地导纳。
[0066] 若节点为PI节点,将式(4)代入式(2)并构造相应的全纯函数可得:
[0067]
[0068] 式中Pi为i节点的注入有功功率,Qi(s)为i节点的注入无功功率Qi(s)=Qi[0]+Qi2 n
[1]s+Qi[2]s+…+Qi[n]s。
[0069] 由PI节点的约束条件:
[0070]
[0071] 假设:
[0072]
[0073] 当s=0时,由式(5)、(6)可以得到:
[0074]
[0075] 由式(9)求解所有节点的Vi[0]
[0076] Vi[0]=1 (10)
[0077] 将式(4)代入式(7)可得:
[0078]
[0079] 展开式为
[0080]
[0081] 由式(12)根据s级数的系数相等得到:
[0082]
[0083]
[0084]
[0085] 若节点为PQ节点,将式(4)带入式(5)得到:
[0086]
[0087] 再将式(8)代入式(16)得到:
[0088]
[0089] 由式(17)根据s级数的系数相等得到:
[0090]
[0091] 若节点为PI节点,将式(4)带入式(6)得到:
[0092]
[0093] 再将式(8)代入式(19)得到:
[0094]
[0095] 由式(20)根据s级数的系数相等得到:
[0096]
[0097] 由
[0098]
[0099] 根据s级数的系数相等,可以得到:
[0100]
[0101] 步骤三、根据构造的全纯函数潮流计算模型递推计算潮流。
[0102] 1)首先根据式(10)求出Vk[i](i≤n‑1),根据式(23)求出dk[i](i≤n‑1),对于PI节点根据式(13)‑(15)求出Qk[i](i≤n‑1);
[0103] 2)利用式(18)和式(21)的递推方程求出Vk[n];
[0104] 重复1)、2)以上二步,逐次求出所有的Vk[n]、Qi[n],i∈p。
[0105] 从而可以得到:
[0106]
[0107] 当s=1时,可以得到潮流的解。
[0108] HELM潮流解法不依靠初始值设定,HELM可以清楚的告知潮流解是否存在,可以预测电压崩溃点。
[0109] 以下结合附图以一具体应用实例做进一步说明。
[0110] a、根据图1原始数据建立电网的数学模型,原始数据如表1所示。该电网既包含PQ节点也包含PI节点,其中1节点是平衡节点,16节点是PI节点(加入一个P=300KW,I=50A的分布式电源),其他节点是PQ节点。
[0111] 若节点为PQ节点,由节点功率方程
[0112]
[0113] 式中:m为PQ节点集合,i、k为节点号;Vi、Vk分别为节点i和节点k的电压;Yik为节点i、节点k之间的导纳;Si为节点i的注入功率;N为节点总数;上标*表示共轭关系。
[0114] 若节点为PI节点,由节点功率方程:
[0115]
[0116]
[0117] 式中p为PI节点的集合,Pi表示注入功率的实部, 表示i节点的电流幅值;
[0118] b、利用HELM方法构造全纯函数数学模型。
[0119] c、根据构造的全纯函数潮流计算模型递推计算潮流,HELM方法和牛拉法计算包含PI节点的电网潮流结果如图2所示。HELM方法和牛拉法计算结果误差如图3所示。
[0120] 从图2和图3可以看出该配电网HELM潮流计算结果与牛拉法潮流计算结果基本一致,说明本文发明的利用HELM方法计算包含PI节点的电网潮流方法是正确可行的。
[0121] 表1 33节点电网系统标准数据
[0122]
[0123]