一种含下垂特性电源的船舶电网潮流计算方法与流程

本发明涉及电力系统技术领域，尤其涉及船舶微电网潮流计算领域，涉及一种船舶电网潮流计算方法，特别是一种含下垂特性电源的船舶电网潮流计算方法。

背景技术：

在大型船舶电网中往往含有几个分布式电源，这些分布式电源接在不同节点，且容量有限，不能像传统的路上电力系统那样，某些电源的容量非常大，可以看成平衡节点，即能够补偿整个系统的功率不平衡。在船舶电网中，分布式电源往往按照下垂控制特性进行有功无功的分配，因此是不存在真正的平衡节点的。另外在船舶电网是低压电网，馈线电阻往往大于电感，其分布式电源的有功主要与电压有关(p-v)，而无功主要与系统频率有关(q-f)。

从当前的微电网潮流计算研究结果来看，有文献对牛顿拉夫逊法加以改进运用，这种方法不适用于船舶电网中线路电阻大于电感的情况。有部分研究采用优化算法求解，该类算法均存在参数过多，调参复杂的问题，很难获得精确解。另有文献采用了前推回代法这种适用于船舶网络特征的方法，但讨论的是电压等级较高网络中，输电线路电感远远大于电阻的情形，下垂特性主要考虑的是p-f，q-v之间的关系，不适用于低压的船舶电网。

技术实现要素：

针对上述现有技术，本发明要解决的技术问题是针对船舶电网的特征，即辐射状网络，馈线电阻大于电感的情况，提供一种继承了前推回代法易于收敛的特点、计算精度高的含下垂特性电源的船舶电网潮流计算方法。

为解决上述技术问题，本发明的一种含下垂特性电源的船舶电网潮流计算方法,包括以下步骤：

步骤1：选择网络中任意一个节点作为准平衡节点，设电网节点数为n，发电机节点数为g，电网各个节点电压ui的初始标么值为1∠0°，确定分布式电源初始有功功率和无功功率的给定值pgi,qgi，i∈g；

步骤2：利用前推回代法进行潮流计算，获得网络有功损耗ploss、无功损耗qloss，并进一步计算总负荷及网络损耗的无功功率与所有电源总无功功率的差额δq；

步骤3：计算系统频率修正量δf，并进一步计算各个电源节点的无功功率修正量δqgi；

步骤4：若满足δqgi≤ε，ε为给定误差，则进入下一步，否则更新各个电源无功功率qgi，然后转至步骤2；

步骤5：计算总负荷及网络损耗的有功功率与电源总有功出力的差额δp，计算准平衡节点电压修正量δu1；

步骤6：更新准平衡节点电压u1，然后进行前推回代潮流计算，求取各个节点新的电压值ui'，进而获得各个节点电压修正量δui，利用电压修正量计算各个电源节点的有功功率修正量δpgi；

步骤7：若满足δpgi≤ε，则终止计算，输出各个节点电压，电源节点有功功率和无功功率；若不满足δpgi≤ε，则更新各个电源有功功率pgi，转到步骤2。

本发明还包括：

1.步骤2中网络有功损耗ploss、无功损耗qloss满足：

ploss＝∑iij²·rij

qloss＝∑iij²·xij

其中iij为支路电流，rij+xij为支路阻抗，i,j∈n；

步骤2中总负荷及网络损耗的无功功率与所有电源总无功功率的差额δq满足：

其中，qgi为电源节点的无功功率；qlj为节点的无功负荷；qloss为整个网络无功损耗。

2.步骤3中系统频率修正量δf满足：

其中mqi为电源无功-频率的下垂系数；

各个电源节点的无功功率修正量δqgi满足：

3.步骤4中更新各个电源无功功率qgi具体为将qgi更新为qgi+δqgi，i∈g。

4.步骤5中总负荷及网络损耗的有功功率与电源总有功出力的差额δp满足：

其中plj为各网络节点负荷有功功率。pgi为各发电机有功功率；

步骤5中准平衡节点电压修正量δu1满足：

其中mpi为电源有功-电压的下垂系数。

5.步骤6中更新准平衡节点电压u1具体为将u1替换为u1+δu1；

步骤6中各个节点电压变化量满足：

δui＝ui'-ui

步骤6中的各个电源节点的有功功率修正量δpgi满足：

6.步骤7中更新各个电源有功功率pgi具体为：将pgi更新为pgi+δpgi。

本发明有益效果：本发明针对船舶电网的辐射状网络结构特征及馈线电阻大于电抗的特征，提供了一种含下垂控制特性电源的船舶电网潮流计算方法。算法精确性高，通过大量的算例仿真验证，该方法与pscad仿真软件的计算结果(通常被认为是精确解)相比较，最大电压幅值绝对误差为0.0003，最大相角误差为0.005。本方法由于是在传统前推回代法基础上的改进方法，易于在传统程序中加以修改，易于理解，且继承了前推回代法易于收敛的特点。

附图说明

图1为计算仿射潮流的算法流程图。

图2为6节点系统接线图。

具体实施方式

下面结合附图对本发明具体实施方式做进一步说明。

结合图1，本发明的目的是这样实现的：基于辐射状网络的经典前推回代法，分别利用内环修正各个电源无功功率，利用外环修正各个电源有功功率。

步骤1：首先选择网络中任意一个节点作为准平衡节点，设电网节点数为n，发电机节点数为g，各个节点电压初始标么值为1∠0°，确定分布式电源初始功率给定值pgi,qgi。

步骤2：利用传统前推回代法进行潮流计算，获得整个网络有功损耗ploss、无功损耗qloss。并进一步计算总负荷及网络损耗的无功功率与所有电源总无功功率的差额。

步骤3：计算系统频率修正量δf，并进一步计算各个电源节点的无功修正量。

步骤4：若满足δqgi≤ε，则进入下一步，否则更新各个电源无功功率，然后转到第2步。

步骤5：计算总负荷及网络损耗的有功功率与电源总有功功率的差额，计算准平衡节点电压修正量。

步骤6：更新准平衡节点电压，然后进行前推回代潮流计算，求取各个节点新的电压值ui'，进而获得各个节点电压修正量δui＝ui'-ui，利用电压变化量计算各个电源节点的有功修正量。

步骤7：若满足δpgi≤ε，则终止计算，输出各个节点电压，功率等变量值。若不满足δpgi≤ε，则更新各个电源有功功率，转到第2步。

本发明的实质是对传统前推回代方法进行修正，以反映船舶电网有功-电压下垂控制特性，无功-频率下垂控制特性。在求取频率变化量和准平衡节点电压变化量时，运用下垂特性公式如下：

其中，n为电网节点数，g为发电机节点，pgi、qgi分别为第i个电源节点的有功、无功功率；plj、qlj分别为第j个节点的有功、无功负荷；ploss，qloss为整个网络有功、无功损耗，mqi为电源无功-频率的下垂系数，mpi为电源有功-电压的下垂系数。

在更新各个节点电压时，通过潮流计算求得各个节点实际电压，而不是用统一的准平衡节点电压变化量，提高计算的精确性。

迭代过程中，设计单独的电压外环，防止了节点电压收敛过程中的震荡，收敛快速。

本发明的工作原理是：

1.选择网络中任意一个节点作为准平衡节点，电网节点数为n，发电机节点数为g，设定各个节点电压初始标么值为1∠0°，确定分布式电源初始功率给定值pgi,qgi。

2.已知电源节点发电功率、节点电压、各节点负荷功率及网络阻抗条件下，利用传统前推回代法进行潮流计算，获得整个网络有功损耗ploss、无功损耗qloss

ploss＝∑iij²·rij(1)

qloss＝σiij²·xij(2)

其中iij为支路电流，rij+xij为支路阻抗，i,j∈n

3.计算总负荷及网络损耗的无功功率与所有电源总无功功率的差额，如下式

其中qgi为电源节点的无功功率；qlj为节点的无功负荷；qloss为整个网络无功损耗。

4.利用无功-频率下垂特性公式，计算系统频率修正量δf

其中mqi为电源无功-频率的下垂系数

5.利用频率修正量计算各个电源节点的无功修正量。

6.若满足δqgi≤ε，则进入步骤7，否则更新各个电源无功功率

qgi＝qgi+δqgii∈g(6)

然后转到第2步。

7.计算总负荷及网络损耗的有功功率与电源总有功功率的差额，如下式

其中plj为负荷有功功率。pgi为发电机有功功率。

利用有功-电压的下垂特性公式，计算准平衡节点电压修正量

其中mpi为电源有功-电压的下垂系数。

8.利用准平衡节点电压修正量更新准平衡节点电压：

u1＝u1+δu1(9)

然后进行前推回代潮流计算，求取各个节点新的电压值ui'，进而获得各个节点电压修正量δui＝ui'-ui，利用电压修正量计算各个电源节点的有功修正量。

9.若满足δpgi≤ε，则终止计算，输出各个节点电压，功率等变量值。若不满足δpgi≤ε，则更新各个电源有功功率

pgi＝pgi+δpgi(11)

再转到第2步。

在6节点的交流系统中，该系统的接线图如图2所示。具体步骤如下：

步骤1：交流网络中节点数为6、支路数为5、支路阻抗值z12＝r12+jx12＝0.43+j0.02ω，z14＝r14+jx14＝0.43+j0.02ω，z23＝r23+jx23＝0.43+j0.02ω，z25＝r25+jx25＝0.44+j0.01ω，z36＝r36+jx36＝0.44+j0.01ω，节点负荷功率值sl1＝pl1+jql1＝3+j5kw，sl3＝pl3+jql3＝8+j3kw，sl2＝sl4＝sl5＝sl6＝0，下垂特性电源的有功-电压下垂系数mp1＝mp2＝mp3＝-0.0012、无功-频率下垂系数mq1＝mq2＝mq3＝-0.0023。基准值选取：sb＝10kva，电压基准值选取：ub＝220v。设定各个电源初始有功pgi＝0，无功qgi＝0，设定迭代次数k＝1，收敛值为0.0001。

步骤2：选定网络中某一节点为准平衡节点，这里选取节点4为准平衡节点，其初始电压标么值为u4＝1，相角为零。确定下垂特性电源的连接节点分别为4,5,6，给定各节点初始电压幅值为1.系统初始频率标么值为1。

步骤3：进行传统前推回代法潮流计算，并计算网络损耗qloss，ploss。

步骤4：内环计算：计算各个电源无功出力与负荷及网络损耗之间无功功率差值

利用下垂特性求系统频率变化量

利用i∈4,5,6，求下垂特性分布式电源节点无功变化量δqgii∈4,5,6，更新电源节点发出无功qgi^(k)＝qgi^(k-1)+δqgii∈4,5,6，返回步骤3，直至满足δqgi≤0.0001，i∈4,5,6为止。

步骤5：外环计算：计算各个电源有功出力与负荷及网损之间有功功率差值

利用有功-电压的下垂特性公式，计算准平衡节点电压变化量

利用准平衡节点电压变化量更新准平衡节点电压：

u4^(k)＝u4^(k-1)+δu4

然后进行前推回代潮流计算，求取各个节点新的电压值ui，进而获得各个节点电压变化量δui＝ui-ui-1，利用电压变化量计算各个电源节点的有功变化量。

更新各个电源有功功率

pgi＝pgi+δpgi，i∈4,5,6

再返回步骤3，直至各个电源有功功率变化量δpgi≤0.0001，i∈4,5,6为止。

表1为本发明的结果与pscad仿真结果的对比，可看出最大电压幅值误差为0.0003，最大电压相角误差为0.005°，表明算法精确度高。

表1本发明的结果与pscad仿真结果对比

本发明具体实施方式还包括：

本发明实施方式包括以下步骤：

步骤1：首先选择网络中任意一个节点作为准平衡节点，设电网节点数为n，发电机节点数为g，各个节点电压初始标么值为1∠0°，确定分布式电源初始功率给定值pgi,qgi。

步骤2：利用传统前推回代法进行潮流计算，获得整个网络有功损耗ploss、无功损耗qloss。并进一步计算总负荷及网络损耗的无功功率与所有电源总无功功率的差额δq。

步骤3：计算系统频率修正量δf，并进一步计算各个电源节点的无功修正量

步骤4：若满足δqgi≤ε，则进入下一步，否则更新各个电源无功功率，然后转到第2步。

步骤5：计算总负荷及网络损耗的有功功率与电源总有功出力的差额，计算准平衡节点电压修正量

步骤6：更新准平衡节点电压，然后进行前推回代潮流计算，求取各个节点新的电压值ui'，进而获得各个节点电压修正量，利用电压修正量计算各个电源节点的有功修正量。

步骤7：若满足δpgi≤ε，则终止计算，输出各个节点电压，功率等变量值。若不满足δpgi≤ε，则更新各个电源有功功率，转到第2步。

步骤2的整个网络有功损耗ploss、无功损耗qloss具体为：

ploss＝∑iij²·rij

qloss＝σiij²·xij

其中iij为支路电流，rij+xij为支路阻抗，i,j∈n

步骤2所述总负荷及网络损耗的无功功率与所有电源总无功功率的差额，具体为：

其中qgi为电源节点的无功出力；qlj为节点的无功负荷；qloss为整个网络无功损耗。

步骤3中系统频率修正量δf具体为：

其中mqi为电源无功-频率的下垂系数。

步骤3中各个电源节点的无功修正量，具体为：

步骤4中ε为允许误差值，ε＝0.0001。

步骤4中更新各个电源无功功率，具体为

qgi＝qgi+δqgii∈g

步骤5中总负荷及网络损耗的有功功率与电源总有功出力的差额，具体为

其中plj为负荷有功功率。pgi为发电机有功功率。

步骤5中准平衡节点电压修正量，具体为

其中mpi为电源有功-电压的下垂系数。

步骤6中更新准平衡节点电压，具体为：

u1＝u1+δu1

步骤6中各个节点电压变化量，具体为：

δui＝ui'-ui

步骤6中各个电源节点的有功修正量，具体为：

步骤7中ε为允许误差值，ε＝0.0001。

步骤7中更新各个电源有功功率，具体为：

pgi＝pgi+δpgi。