首页 > 专利 > 杭州电子科技大学 > 考虑PI节点的HELM潮流计算方法专利详情

考虑PI节点的HELM潮流计算方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2021-06-22
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2021-09-24
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2022-04-22
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2041-06-22
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN202110691838.X 申请日 2021-06-22
公开/公告号 CN113363990B 公开/公告日 2022-04-22
授权日 2022-04-22 预估到期日 2041-06-22
申请年 2021年 公开/公告年 2022年
缴费截止日
分类号 H02J3/06 主分类号 H02J3/06
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 0
权利要求数量 1 非专利引证数量 1
引用专利数量 1 被引证专利数量 0
非专利引证 1、2014.06.05赵奕岚.基于 HELM 网损灵敏度的配电网无功优化研究《.浙江电力》.2020,李思儒.基于全纯函数嵌入的暂态稳定网络方程算法《.能源互联网》.2021,Muthu Kumar Subramanian.PV busmodeling in a holomorphically embeddedpower-flow formulation《.2013 NorthAmerican Power Symposium (NAPS)》.2013,;
引用专利 US2014156094A 被引证专利
专利权维持 1 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 高慧敏、钟毅、赵嘉敏、黄卓、罗平 第一发明人 高慧敏
地址 浙江省杭州市下沙高教园区二号路 邮编 310018
申请人数量 1 发明人数量 5
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
浙江永鼎律师事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
陆永强
摘要
本发明公开了考虑PI节点的HELM潮流计算方法。该方法首先建立包含PQ和PI节点的配电网数学模型,然后根据HELM方法构造PI节点对应的全纯函数,最后根据递推方法推导出PI节点全纯函数递推公式来计算包含PI节点的配电网潮流。本发明利用HELM法考虑了PI节点数学模型,推导出PI节点全纯函数递推公式,可以计算包含PI节点的配电网潮流,可用于包含PI分布式电源的配电网潮流计算和电压分析控制。
  • 摘要附图
    考虑PI节点的HELM潮流计算方法
  • 说明书附图:图1
    考虑PI节点的HELM潮流计算方法
  • 说明书附图:图2
    考虑PI节点的HELM潮流计算方法
  • 说明书附图:图3
    考虑PI节点的HELM潮流计算方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2022-04-22 授权
2 2021-09-24 实质审查的生效 IPC(主分类): H02J 3/06 专利申请号: 202110691838.X 申请日: 2021.06.22
3 2021-09-07 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.考虑PI节点的HELM潮流计算方法,其特征在于:具体包括以下步骤:
步骤一,建立电网的数学模型,建立包含PI节点的电网的数学模型:
若节点为PQ节点,由节点功率方程:
式中m为PQ节点集合,i、k为节点号;Vi、Vk分别为节点i和节点k的电压;Yik为节点i、节点k之间的导纳;Si为节点i的注入功率;N为节点总数;上标*表示共轭关系;
若节点为PI节点,由节点功率方程:
式中p为PI节点的集合,Pi表示i节点的注入有功功率,Ii表示i节点的注入电流,|Ii|表示其绝对值, 表示i节点的电流幅值,Vi、Vk分别为节点i和节点k的电压;Yik为节点i、节点k之间的导纳,上标*表示共轭关系;N为电网节点总数;Re(·)表示复数向量的实部;
步骤二,利用HELM方法构造PI节点的内嵌纯虚函数:
其中Vi[n]表示HELM潮流计算中的i节点电压的第n个电压分项,sn表示频域算子s的n次项;
若节点为PQ节点,将式(4)代入式(1)并构造相应的全纯函数得到:
式中: Yii为节点导纳矩阵中i节点的自导纳;Yik为节点导纳
矩阵中i和k节点间的互导纳;Yi,shunt为i节点的接地导纳;
若节点为PI节点,将式(4)代入式(2)并构造相应的全纯函数得到:
式中Pi为i节点的注入有功功率,Qi(s)为i节点的注入无功功率Qi(s)=Qi[0]+Qi[1]s+
2 n
Qi[2]s+…Qi[n]s;
由PI节点的约束条件:
假设:
当s=0时,由式(5)、(6)得到:
由式(9)求解所有节点的Vi[0]
Vi[0]=1              (10)
将式(4)代入式(7)得到:
展开式为
由式(12)根据s级数的系数相等得到:
若节点为PI节点,将式(4)代入式(6)得到:
再将式(8)代入式(19)得到:
由式(20)根据s级数的系数相等得到:

根据s级数的系数相等,得到:
步骤三,根据构造的全纯函数潮流计算模型递推计算潮流;
1)首先根据式(10)求出Vk[i](i≤n‑1),根据式(23)求出dk[i](i≤n‑1),对于PI节点根据式(13)‑(15)求出Qk[i](i≤n‑1);
2)利用式(21)的递推方程求出Vk[n];
重复1)、2)以上二步,逐次求出所有的Vk[n]、Qi[n],i∈p,
从而得到:
当s=1时,得到潮流的解。
说明书

技术领域

[0001] 本发明属于电力信息技术领域,具体涉及一种考虑PI节点的HELM潮流计算方法。

背景技术

[0002] 随着新能源技术在配电网领域的发展,分布式电源(Distributed Generation,DG)逐渐大量接入。分布式电源是指利用可再生能源并且发电功率在50MW以下小型发电机组,放置在用户附近,开展短距离供电工作,并将其就近接入配电网当中以满足用户的电力需求。分布式电源一般包括光伏发电系统、微型燃气轮机、燃料电池、风力机组等。传统的PQ分解法和前推回代法求解电网潮流中可以求解包含PQ、PV节点的潮流,但包含DG的配电网潮流计算与普通潮流计算的区别之一是DG的潮流计算模型与传统发电机组计算模型不一致。传统发电机节点在潮流计算中一般取为PQ节点,PV节点、或平衡节点,而DG有特殊性,其中双馈感应机和同步风机能够通过频率转换控制系统的转换信号对功率做出相应调整与控制,通常认为是PQ节点;微型燃气轮机、燃料电池、电压控制的光伏发电系统通常认为是PV节点;光伏发电系统如果采用电流控制逆变器策略,则认为是PI节点,输出有功功率及电流注入到配电网。包含DG的电网潮流也可以用PQ分解法和前推回代法求解包含PQ、PV节点的潮流,将PI节点在迭代中转化为PQ节点来处理。
[0003] 最近两年提出了一种基于全纯函数嵌入(Holomorphic embedding method,HELM)的潮流计算方法,这种潮流方法完全颠覆了传统的牛拉法,不需要依靠节点初始值,就可以明确潮流解是否存在。这种方法可以完全改变传统的电力系统潮流计算,无功优化,电压稳定分析等问题。目前HELM可以求解包含PQ、PV节点的潮流,但如何利用HELM计算包含PI节点的电网潮流,目前还没有相关文献研究。由于HELM方法是递归算法,并不是迭代算法,无法用常规的将PI节点在迭代中转化为PQ节点的方法来处理。

发明内容

[0004] 针对HELM潮流计算现有技术的不足,本发明提出了一种考虑PI节点的HELM潮流计算方法,可以用于包含PI节点的分布式电源的潮流计算,以及包含分布式电源的配电网电压稳定预防控制,计算速度快,计算简单,有较高的理论意义和应用价值。
[0005] 本发明技术方案具体是:
[0006] 一种考虑PI节点的电网HELM潮流计算方法,包括以下步骤:
[0007] 步骤一,建立电网的数学模型,建立包含PI节点的电网的数学模型:
[0008] 若节点为PI节点,由节点功率方程:
[0009]
[0010]
[0011] 式中p为PI节点的集合,Pi表示i节点的注入有功功率,Ii表示i节点的注入电流,|Ii|表示其绝对值, 表示i节点的电流幅值,Vi、Vk分别为节点i和节点k的电压;Yik为节点i、节点k之间的导纳,上标*表示共轭关系;N为电网节点总数;Re(·)表示复数向量的实部。
[0012] 步骤二,利用HELM方法构造PI节点的内嵌纯虚函数:
[0013]
[0014] 其中Vi[n]表示HELM潮流计算中的i节点电压的第n个电压分项,sn表示频域算子s的n次项;
[0015] 若节点为PI节点,将式(4)代入式(2)并构造相应的全纯函数可得:
[0016]
[0017] 式中Pi为i节点的注入有功功率,Qi(s)为i节点的注入无功功率Qi(s)=Qi[0]+Qi2 n
[1]s+Qi[2]s+…+Qi[n]s , Yii为节点导纳矩阵中i节点的自导
纳;Yik为节点导纳矩阵中i和k节点间的互导纳;Yi,shunt为i节点的接地导纳。
[0018] 由PI节点的约束条件:
[0019]
[0020] 将式(4)代入式(7)可得:
[0021]
[0022] 展开式为
[0023]
[0024] 由式(12)根据s级数的系数相等得到:
[0025]
[0026]
[0027]
[0028] 若节点为PI节点,将式(4)带入式(6)得到:
[0029]
[0030] 再将式(8)代入式(19)得到:
[0031]
[0032] 由式(20)根据s级数的系数相等得到:
[0033]
[0034] 由
[0035]
[0036] 根据s级数的系数相等,可以得到:
[0037]
[0038] 步骤三,根据构造的全纯函数潮流计算模型递推计算潮流;
[0039] 1)首先根据式(10)求出Vk[i](i≤n‑1),根据式(23)求出dk[i](i≤n‑1),对于PI节点根据式(13)‑(15)求出Qk[i](i≤n‑1);
[0040] 2)利用式(21)的递推方程求出Vk[n];
[0041] 重复1)、2)以上二步,逐次求出所有的Vk[n]、Qi[n],i∈p。
[0042] 从而可以得到:
[0043]
[0044] 当s=1时,可以得到潮流的解。
[0045] 本发明具有以下有益效果:可以用于包含PI节点分布式电源的配电网的潮流计算、电压灵敏度计算、电压稳定分析及电压控制,包含分布式电源的配电网的电压分析控制问题,新能源的接入和电压稳定分析及控制等问题,使无功补偿点的选取更为合理,补偿方案经济性更好,进一步降低网损,提高电压水平和功率因数,对电网的无功优化及电压稳定分析控制有较高的理论意义和应用价值。

实施方案

[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]

附图说明

[0046] 图1为具体应用实例中为原始数据建立33节点电网的数学模型。
[0047] 图2为一具体应用实例中潮流结果示意图;
[0048] 图3为包含HELM方法与牛拉法潮流结果示意图。