(19)中华人民共和国国家知识产权局
(12)发明专利申请
(10)申请公布号 CN 112531715 A(43)申请公布日 2021.03.19
(21)申请号 202011293340.X(22)申请日 2020.11.18
(71)申请人 天津大学
地址 300072 天津市南开区卫津路92号 申请人 国网新源张家口风光储示范电站有
限公司(72)发明人 王守相 刘琪 赵倩宇 王振宇
田云峰 沈宇 董文琦 马步云 (74)专利代理机构 天津才智专利商标代理有限
公司 12108
代理人 张文华(51)Int.Cl.
H02J 3/06(2006.01)
权利要求书2页 说明书5页
(54)发明名称
基于虚拟电阻的下垂控制多端直流微电网潮流计算方法(57)摘要
本发明涉及基于虚拟电阻的下垂控制多端直流微电网潮流计算方法，包括如下步骤：获取直流微电网拓扑结构、换流阀设备参数、线路参数、负荷参数及分布式电源参数；根据换流阀下垂控制参数构建虚拟电阻模型；得到直流微电网潮流计算的雅可比矩阵；采用牛顿‑拉夫逊迭代法计算节点电压，当前后两次迭代计算的节点电压差满足收敛条件时，计算完成，输出潮流计算结果。本发明通过在直流潮流模型中引入虚拟电阻，将下垂控制模型转化为精确的潮流模型，避免了直流潮流计算中松弛母线电压的更新，提高了计算速度。
CN 112531715 ACN 112531715 A
权 利 要 求 书
1/2页
1.基于虚拟电阻的下垂控制多端直流微电网潮流计算方法，其特征在于，包括如下步骤：
获取直流微电网拓扑结构、换流阀设备参数、线路参数、负荷参数及分布式电源参数；根据换流阀下垂控制参数构建虚拟电阻模型；得到直流微电网潮流计算的雅可比矩阵；采用牛顿‑拉夫逊迭代法计算节点电压，当前后两次迭代计算的节点电压差满足收敛条件时，计算完成，输出潮流计算结果。
2.根据权利要求1所述的基于虚拟电阻的下垂控制多端直流微电网潮流计算方法，其特征在于：根据换流阀下垂控制参数构建虚拟电阻模型的具体方法为：
在下垂控制的直流微电网中，第n个换流阀的直流侧电压UVSC,n和输出功率PVSC,n满足：
式中，UVSC,n,ref和PVSC,n,ref分别是第n个换流阀的运行参考电压和运行参考功率；kVSC,n是第n个换流阀的下垂系数；
在所构造的虚拟电阻模型中，换流阀下垂控制模型相当于一个理想的电压源和一条线路，下垂节点相当于负荷节点；其中，理想电压源的输出电压即为换流阀的运行参考电压，同时，通过向下垂节点增加负的负载来补偿虚拟电阻的功率损耗，从而保证下垂节点电压和虚拟线路的传输功率满足下垂控制方程；
对于第n个换流阀，抵消等效线路损耗的等效荷载Ploss,n可表示为：
式中，RVSC,n为第n个换流阀的等效线路电阻；UVSC,n,sou为第n个换流阀的理想电压源的电压；
第n个换流阀下垂节点视为负荷节点，负荷大小为：Pload,n＝‑Ploss,n‑PVSC,n,ref (3)
第n个换流阀的虚拟电阻的潮流模型为：
其中，Pn是第n个换流阀的虚拟线路的传输功率。根据(1)，当换流阀的输出功率变化ΔP时，换流阀的输出电压变化量ΔUVSC,n可以表示为：
ΔUVSC,n＝kVSC,nΔP (5)由(4)可知，第n个换流阀由于换流阀输出功率的变化而产生的输出电压变化可以表示为：
因此，第n个换流阀下垂控制对应的虚拟电阻大小为：RVSC,n＝‑kVSC,nUVSC,n,sou (7)。
2
CN 112531715 A
权 利 要 求 书
2/2页
3.根据权利要求1所述的基于虚拟电阻的下垂控制多端直流微电网潮流计算方法，其特征在于：令负荷节点序号为1，2，…，h，下垂节点的序号为h+1，h+2，…，h+m，通过虚拟电阻的应用，在直流配电网中引入了m个理想电压源，分布式电源输出功率假定为负的负载，第i条负载母线的节点功率Pi可表示为：
Pi＝PDG,i‑PL,i i＝1,2,L,h (8)其中，PDG,i和PL,i是位于节点i的分布式电源输出功率和负载功率；设第n个换流阀的下垂节点和理想电压源节点分别为节点i和节点k，下垂节点的节点功率可以表示为：
Pi＝PVSC,n,ref+gVSC,n(UVSC,n,ref‑Un)2 i＝h+1,h+2,L,h+m (9)式中，gVSC,n是第n个换流阀的虚拟线路的导纳；直流微电网模型为：
在每次迭代中，负荷节点和下垂节点的节点功率变化可以表示为：
牛顿‑拉夫逊迭代方法的迭代形式为：
其中，ΔF＝[ΔP1,ΔP2,L,ΔPi],x＝[U1,U2,L,Ui]，k是迭代次数，J是雅可比矩阵；根据(11)，雅可比矩阵中的元素可以表示为：
4.根据权利要求1所述的基于虚拟电阻的下垂控制多端直流微电网潮流计算方法，其特征在于：采用牛顿‑拉夫逊迭代法计算节点电压，当前后两次迭代计算的节点电压差小于10‑8p.u.时，计算完成，输出潮流计算结果。
3
CN 112531715 A
说 明 书
1/5页
基于虚拟电阻的下垂控制多端直流微电网潮流计算方法
技术领域
[0001]本发明属于微电网计算技术领域，尤其涉及基于虚拟电阻的下垂控 制多端直流微电网潮流计算方法。
背景技术
[0002]随着直流型分布式电源的大规模接入和电动汽车、数据中心等直流 负荷的不断增加，配电系统的直流特征日益突出。发展直流微电网能有 效提高对光伏等直流型分布式电源的接纳能力，减少换流环节，降低线 路损耗，提高能源利用效率。[0003]目前，潮流计算是电力系统的基本工具，前推回代法和牛顿迭代法 是配电网最常用的潮流算法：前推回代法计算速度快，适用于辐射状配 电网，虽然牛顿迭代法比前推回代法需要更多的执行时间来收敛，但是 牛顿迭代法的有效性不受网络拓扑结构的影响。因此，基于牛顿迭代法 的潮流算法更适合于直流配电网的潮流计算。[0004]在直流微电网中，主从控制和下垂控制是两种常见的控制方式。在 下垂控制的直流微电网中，交/直流变换器的输出电压与交/直流变换器 的输出功率有关。因此，在直流微电网中没有松弛节点维持母线节点电 压恒定。
[0005]在已有的基于下垂控制的直流微电网潮流计算方法中，通常首先假 设一个最大容量的交/直流变换器的输出电压是恒定的。然后，根据每次 迭代的结果，得到该交/直流变换器的输出功率，并更新相应的输出电压。 经过多次迭代，满足了交/直流变换器输出电压的收敛条件，完成潮流计 算。在上述算法中，潮流计算多次迭代增加了计算负担，因此，基于这 些问题，提供一种更为快速准确的直流微电网潮流计算方法，具有重要 的现实意义。
发明内容
[0006]本发明的目的在于克服现有技术的不足，引入虚拟电阻模型来表示 直流微电网
偿虚拟电阻的功率损耗，保证下垂控制对潮流的影响，通过向下垂节点加入注入功率来补
了虚拟电阻模型的精度，在此模型基础上 提出了基于牛顿‑拉夫逊迭代法的潮流计算方法。
[0007]本发明解决其技术问题是采取以下技术方案实现的：
[0008]基于虚拟电阻的下垂控制多端直流微电网潮流计算方法，包括如下 步骤：[0009]获取直流微电网拓扑结构、换流阀设备参数、线路参数、负荷参数 及分布式电源参数；
[0010]根据换流阀下垂控制参数构建虚拟电阻模型；[0011]得到直流微电网潮流计算的雅可比矩阵；[0012]采用牛顿‑拉夫逊迭代法计算节点电压，当前后两次迭代计算的节点 电压差满足收敛条件时，计算完成，输出潮流计算结果。[0013]进一步的，根据换流阀下垂控制参数构建虚拟电阻模型的具体方法 为：
4
CN 112531715 A[0014]
说 明 书
2/5页
在下垂控制的直流微电网中，第n个换流阀的直流侧电压UVSC,n和输出 功率PVSC,n
满足：
[0015]
式中，UVSC,n,ref和PVSC,n,ref分别是第n个换流阀的运行参考电压和运行参考 功率；
kVSC,n是第n个换流阀的下垂系数；
[0017]在所构造的虚拟电阻模型中，换流阀下垂控制模型相当于一个理想 的电压源和一条线路，下垂节点相当于负荷节点；其中，理想电压源的 输出电压即为换流阀的运行参考电压，同时，通过向下垂节点增加负的 负载来补偿虚拟电阻的功率损耗，从而保证下垂节点电压和虚拟线路的 传输功率满足下垂控制方程；[0018]对于第n个换流阀，抵消等效线路损耗的等效荷载Ploss,n可表示为：
[0019][0020]
[0016]
式中，RVSC,n为第n个换流阀的等效线路电阻；UVSC,n,sou为第n个换流阀 的理想电压源的电压；
[0021]第n个换流阀下垂节点视为负荷节点，负荷大小为：[0022]P＝‑Ploss,n‑PVSC,n,ref (3)load,n
[0023]第n个换流阀的虚拟电阻的潮流模型为：
[0024]
其中，Pn是第n个换流阀的虚拟线路的传输功率。
[0026]根据(1)，换流阀的输出电压变 化量ΔUVSC,n可以当换流阀的输出功率变化ΔP时，表示为：[0027]ΔU (5)VSC,n＝kVSC,nΔP[0028]由(4)可知，第n个换流阀由于换流阀输出功率的变化而产生的输 出电压变化可以表示为：
[0029]
[0025]
因此，第n个换流阀下垂控制对应的虚拟电阻大小为：
[0031]R (7) 。VSC,n＝‑kVSC,nUVSC,n,sou[0032]进一步的，令负荷节点序号为1，2，…，h，下垂节点的序号为h+1， h+2，…，h+m，通过虚拟电阻的应用，在直流配电网中引入了m个理 想电压源，分布式电源输出功率假定为负的负载，第i条负载母线的节 点功率Pi可表示为：[0033]P＝P i＝1,2,L,h (8)iDG,i‑PL,i[0034]其中，PDG,i和PL,i是位于节点i的分布式电源输出功率和负载功率；[0035]设第n个换流阀的下垂节点和理想电压源节点分别为节点i和节点 k，下垂节点的节点功率可以表示为：
2
[0036]P＝P+g(U‑U) (9)iVSC,n,refVSC,nVSC,n,refn i＝h+1,h+2,L,h+m
5
[0030]
CN 112531715 A[0037][0038][0039][0040]
说 明 书
3/5页
式中，gVSC,n是第n个换流阀的虚拟线路的导纳；直流微电网模型为：
在每次迭代中，负荷节点和下垂节点的节点功率变化可以表示为：
[0041]
[0042][0043][0044]
牛顿‑拉夫逊迭代方法的迭代形式为：
其中，ΔF＝[ΔP1,ΔP2,L,ΔPi],x＝[U1,U2,L,Ui]，k是迭代次数，J是雅可比矩 根据(11)，雅可比矩阵中的元素可以表示为：
阵；
[0045]
[0046]
进一步的，采用牛顿‑拉夫逊迭代法计算节点电压，当前后两次迭 代计算的节点
电压差小于10‑8p.u.时，计算完成，输出潮流计算结果。[0048]本发明的优点和积极效果是：
[0049]本发明通过在直流潮流模型中引入虚拟电阻，将下垂控制模型转化 为精确的潮流模型，避免了直流潮流计算中松弛母线电压的更新，提高 了计算速度。具体实施方式[0050]首先，需要说明的是，以下将以示例方式来具体说明本发明的具体 结构、特点和优点等，然而所有的描述仅是用来进行说明的，而不应将 其理解为对本发明形成任何。此外，在本文所提及各实施例中予以 描述或隐含的任意单个技术特征，仍然可在这些技术特征(或其等同物) 之间继续进行任意组合或删减，从而获得可能未在本文中直接提及的本 发明的更多其他实施例。[0051]需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中 的特征可以相互组合。
[0052]本实施例提供的基于虚拟电阻的下垂控制多端直流微电网潮流计算 方法，包括如下步骤：
[0053]获取直流微电网拓扑结构、换流阀设备参数、线路参数、负荷参数 及分布式电源
6
[0047]
CN 112531715 A
说 明 书
4/5页
参数；
根据换流阀下垂控制参数构建虚拟电阻模型；
[0055]得到直流微电网潮流计算的雅可比矩阵；[0056]采用牛顿‑拉夫逊迭代法计算节点电压，当前后两次迭代计算的节点 电压差满足收敛条件时，计算完成，输出潮流计算结果。[0057]具体的，在本实施例中，根据换流阀下垂控制参数构建虚拟电阻模 型的具体方法为：
[0058]在下垂控制的直流微电网中，第n个换流阀的直流侧电压UVSC,n和输出 功率PVSC,n满足：
[0059][0054]
式中，UVSC,n,ref和PVSC,n,ref分别是第n个换流阀的运行参考电压和运行参考 功率；
kVSC,n是第n个换流阀的下垂系数；[0061]在虚拟电阻模型中，换流阀下垂控制模型相当于一个理想的电压源 和一条线路。其中，理想电压源的输出电压即为换流阀的运行参考电压。
[0062]通过向下垂节点增加负的负载来补偿虚拟电阻的功率损耗，从而保 证下垂节点电压和虚拟线路的传输功率满足下垂控制方程。对于第n个 换流阀，抵消等效线路损耗的等效荷载Ploss,n可表示为：
[0063][00]
[0060]
式中，RVSC,n为第n个换流阀的等效线路电阻；UVSC,n,sou为第n个换流阀 的理想电压
源的电压。
[0065]第n个换流阀下垂节点视为负荷节点，负荷大小为：[0066]P (3)load,n＝‑Ploss,n‑PVSC,n,ref
[0067]第n个换流阀的虚拟电阻的潮流模型为：
[0068]
其中，Pn是第n个换流阀的虚拟线路的传输功率。
[0070]根据(1)，当换流阀的输出功率变化ΔP时，换流阀的输出电压变 化量ΔUVSC,n可以表示为：[0071]ΔU (5)VSC,n＝kVSC,nΔP[0072]由(4)可知，第n个换流阀由于换流阀输出功率的变化而产生的输 出电压变化可以表示为：
[0073][0074][0075][0076]
[0069]
因此，第n个换流阀下垂控制对应的虚拟电阻大小为：RVSC,n＝‑kVSC,nUVSC,n,sou (7)令负荷节点序号为1，2，…，h，下垂节点的序号为h+1，h+2，…， h+m，通过虚拟电阻
7
CN 112531715 A
说 明 书
5/5页
的应用，在直流配电网中引入了m个理想电压源， 分布式电源输出功率假定为负的负载，第i条负载母线的节点功率Pi可 表示为：[0077]P＝P i＝1,2,L,h (8)iDG,i‑PL,i[0078]其中，PDG,i和PL,i是位于节点i的分布式电源输出功率和负载功率；[0079]设第n个换流阀的下垂节点和理想电压源节点分别为节点i和节点 k，下垂节点的节点功率可以表示为：
2
[0080]P＝P+g(U‑U) (9)iVSC,n,refVSC,nVSC,n,refn i＝h+1,h+2,L,h+m[0081]式中，gVSC,n是第n个换流阀的虚拟线路的导纳；[0082]直流微电网模型为：
[0083][0084]
在每次迭代中，负荷节点和下垂节点的节点功率变化可以表示为：
[0085]
[0086][0087][0088]
牛顿‑拉夫逊迭代方法的迭代形式为：
其中，ΔF＝[ΔP1,ΔP2,L,ΔPi],x＝[U1,U2,L,Ui]，k是迭代次数，J是雅可比矩 根据(11)，雅可比矩阵中的元素可以表示为：
阵；
[00]
[0090]
采用牛顿‑拉夫逊迭代法计算节点电压：当前后两次迭代计算的节点 电压差小于
10‑8p.u.时，计算完成，输出潮流计算结果。[0092]相比于传统的针对下垂控制直流配电网的潮流计算，所提出的潮流 计算方法不需要多次潮流迭代修正下垂节点的节点电压，针对多个直流 微电网的仿真对比表明，本发明所提的针对下垂控制的直流微电网的潮 流计算方法的计算速度是已有的潮流计算方法的3倍。
[0093]以上实施例对本发明进行了详细说明，但所述内容仅为本发明的较 佳实施例，不能被认为用于限定本发明的实施范围。凡依本发明申请范 围所作的均等变化与改进等，均应仍归属于本发明的专利涵盖范围之 内。
[0091]
8
