CN111244962B - 一种统一求解的交直流混合配电网潮流计算方法 - Google Patents
一种统一求解的交直流混合配电网潮流计算方法 Download PDFInfo
- Publication number
- CN111244962B CN111244962B CN202010042263.4A CN202010042263A CN111244962B CN 111244962 B CN111244962 B CN 111244962B CN 202010042263 A CN202010042263 A CN 202010042263A CN 111244962 B CN111244962 B CN 111244962B
- Authority
- CN
- China
- Prior art keywords
- power
- node
- equation
- nodes
- branch
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 31
- 238000002940 Newton-Raphson method Methods 0.000 claims abstract description 5
- 238000000034 method Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 239000011541 reaction mixture Substances 0.000 claims description 3
- 238000011217 control strategy Methods 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910002056 binary alloy Inorganic materials 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明是一种统一求解的交直流混合配电网潮流计算方法,其特点是,通过建立交流‑直流变流器的稳态潮流模型,结合二进制描述的网络矩阵,得到了交直流混合配电网拓扑结构的统一潮流计算模型,针对不同控制策略建立对应潮流方程,并基于牛顿‑拉夫逊法对潮流方程进行求解,该方法提供的结果是分析配电网的稳态运行,设计相应控制方式和研究配备相关保护设置的前提和基础,对于不同控制策略和不同拓扑结构的交直流混合配电网均能够迅速、稳定收敛,其计算速度快,适用性强,应用效果佳。
Description
技术领域
本发明涉及交直流混合配电网领域,特别是涉及一种统一求解的交直流混合配电网潮流计算方法。
背景技术
近年来,城市规模日益扩展,用电负荷逐年攀升,重要负荷,非线性负荷的占比越来越大,对城市配电网的供电可靠性和经济性提出了更高的要求。现有的交流配电网面临供电容量不足,线路饱和及电能变换环节多等问题。而直流配电网具有可控性好、损耗低、电能质量高、输电容量大等特点,并且非常适合分布式电源的多点接入,同时,分布式电源、储能装置可从直流侧并入系统,减少了电力电子设备并网的逆变环节,降低成本和损耗,提高供电可靠性,为适应不同类型的电源和负荷,城市配电网逐渐向交直流混合配电网模式发展。
电力系统潮流计算是电力系统运行状态分析的基础,根据给定的各元件状态和网络结构计算电力系统中各节点电压及线路潮流分布,各种稳态暂态运行分析、电力系统规划、可靠性评估等都需要进行潮流计算。现有方法的交直流潮流模型对于网架拓扑结构的适应能力有限,无法完全适应城市配电网的某些特定的拓扑结构。
发明内容
本发明针对现有技术存在的问题,提供一种计算速度快,适用性强,应用效果佳的交直流混合配电网潮流计算方法,该方法提供的结果是分析配电网的稳态运行,设计相应控制方式和研究配备相关保护设置的前提和基础,对于不同控制策略和不同拓扑结构的交直流混合配电网均能够迅速、稳定收敛。
1.实现本发明的技术方案是:一种统一求解的交直流混合配电网潮流计算方法,其特征是,它包括建立二进制网络矩阵、建立交直流支路功率模型、建立功率平衡方程、求解功率平衡方程,具体内容有:
1)建立二进制网络矩阵:
①N行1列节点种类矩阵B,其中,N是节点数,N∈[2,9999];对于直流节点i∈[1,N],j∈[1,N],Bi=1,Bj=1,对于交流节点Bi=0,Bj=0;
②N行N列网络连接矩阵C,对于两不同节点i,j∈[1,N]之间有支路连接,Cij=1;对于两不同节点i,j之间无支路连接,Cij=0;
③N行N列支路种类矩阵L,对于两不同节点i,j之间有直流支路连接,Lij=1;对于两不同节点i,j之间有交流支路连接,Lij=0;
2)建立交直流支路功率模型:
①两个交流节点i,j相连的功率方程,
其中,Pij、Qij分别是i、j两节点间交流支路首端环节输入的有功功率和无功功率;θij是i、j两节点间电压相角差;Vi、Vj分别为i、j两节点电压幅值;Gij、Bij分别是i、j两节点间交流支路电导和电纳,结合网络矩阵,得到的潮流方程为:
其中,Pij (a)、Qij (a)是中间变量;
②两个直流节点i,j相连的其功率方程为,
结合网络矩阵,得到的潮流方程为:
其中,Pij (b)是中间变量;
则:
Vi=MikVk (5)
Pkj=VkIkj=Vk(Vk-Vj)Gkj (6)
其中,Vk是交流-直流换流器交流侧节点i在直流侧等效节点k节点电压幅值,Mik是交流-直流变流器交直流两侧i、k节点电压幅值间折算系数,Gkj是k、j节点间连接的等效电导参数,Ikj是k节点流向j节点的电流,结合网络矩阵,得到潮流方程,
④一个直流节点i和一个交流节点j相连,功率方程为
结合网络矩阵,得到潮流方程,
其中,Pij (d)是中间变量;
⑤两个交流节点i、j通过直流线路相连,列功率方程和潮流方程:
其中,k、k’节点分别是i、j节点交流-直流变流器直流侧对应节点,a2、b2、Pij (e)、Qij (e)是中间变量,
且
3)建立功率平衡方程:
建立交直流混合配电网功率平衡方程
其中:Pis、Qis分别是i节点注入的有功功率和无功功率,且
4)应用牛顿-拉夫逊法求解功率平衡方程:
①求出修正方程的不平衡量;
②求出修正方程的系数矩阵雅可比矩阵;
③应用高斯消去法求解修正方程,得到修正量;
④求出新值;
⑤判断修正量精度,若不满足自第一步迭代计算,若满足根据电功率公式求潮流分布。
本发明一种统一求解的交直流混合配电网潮流计算方法的有益效果体现在:通过建立交流-直流变流器的稳态潮流模型,结合二进制描述的网络矩阵,得到了交直流混合配电网拓扑结构的统一潮流计算模型,针对不同控制策略建立对应潮流方程,并基于牛顿-拉夫逊法对潮流方程进行求解,本发明对于不同控制策略和不同拓扑结构的交直流混合配电网均能够迅速、稳定收敛,其计算速度快,适用性强,应用效果佳。
附图说明
图1为本发明的一种统一求解的交直流混合配电网潮流计算方法框图;
图2为图1的实施例示意图;
图3为图2的节点电压幅值、节点电压相角、节点有功功率和节点无功功率潮流分布示意图;
图4为图2的支路首端有功功率、支路末端有功功率、支路首端无功功率和支路末端无功功率潮流分布示意图。
具体实施方式
以下结合图1至图4和具体实施例对本发明作进一步详细说明,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
2.参照图1,一种统一求解的交直流混合配电网潮流计算方法,它包括的内容有:
1)建立二进制网络矩阵:
①N行1列节点种类矩阵B,其中,N是节点数,N∈[2,9999];对于直流节点i∈[1,N],j∈[1,N],Bi=1,Bj=1,对于交流节点Bi=0,Bj=0;
②N行N列网络连接矩阵C,对于两不同节点i,j∈[1,N]之间有支路连接,Cij=1;对于两不同节点i,j之间无支路连接,Cij=0;
③N行N列支路种类矩阵L,对于两不同节点i,j之间有直流支路连接,Lij=1;对于两不同节点i,j之间有交流支路连接,Lij=0;
2)建立交直流支路功率模型:
①两个交流节点i,j相连的功率方程,
其中,Pij、Qij分别是i、j两节点间交流支路首端环节输入的有功功率和无功功率;θij是i、j两节点间电压相角差;Vi、Vj分别为i、j两节点电压幅值;Gij、Bij分别是i、j两节点间交流支路电导和电纳,结合网络矩阵,得到的潮流方程为:
其中,Pij (a)、Qij (a)是中间变量;
②两个直流节点i,j相连的其功率方程为,
结合网络矩阵,得到的潮流方程为:
其中,Pij (b)是中间变量;
则:
Vi=MikVk (5)
Pkj=VkIkj=Vk(Vk-Vj)Gkj (6)
其中,Vk是交流-直流换流器交流侧节点i在直流侧等效节点k节点电压幅值,Mik是交流-直流变流器交直流两侧i、k节点电压幅值间折算系数,Gkj是k、j节点间连接的等效电导参数,Ikj是k节点流向j节点的电流,结合网络矩阵,得到潮流方程,
④一个直流节点i和一个交流节点j相连,功率方程为
结合网络矩阵,得到潮流方程,
其中,Pij (d)是中间变量;
⑤两个交流节点i、j通过直流线路相连,列功率方程和潮流方程:
其中,k、k’节点分别是i、j节点交流-直流变流器直流侧对应节点,a2、b2、Pij (e)、Qij (e)是中间变量,
且
3)建立功率平衡方程:
建立交直流混合配电网功率平衡方程
其中:Pis、Qis分别是i节点注入的有功功率和无功功率,且
4)应用牛顿-拉夫逊法求解功率平衡方程:
①求出修正方程的不平衡量;
②求出修正方程的系数矩阵雅可比矩阵;
③应用高斯消去法求解修正方程,得到修正量;
④求出新值;
⑤判断修正量精度,若不满足自第一步迭代计算,若满足根据电功率公式求潮流分布。3.实施例:
参照图2,利用本发明的一种统一求解的交直流混合配电网潮流计算方法,在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下,计算有功功率、无功功率及电压在电力系统中的分布。
参照图3,利用本发明的一种电力系统潮流计算方法,计算出的结果,包括节点电压幅值、节点电压相角、节点有功功率和节点无功功率。
参照图4,利用本发明的一种电力系统潮流计算方法,计算出的结果,包括支路首端有功功率、支路末端有功功率、支路首端无功功率和支路末端无功功率。
以上所述仅是本发明的优选实施例,应当指出的是,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应该视为本发明的保护范围。
Claims (1)
1.一种统一求解的交直流混合配电网潮流计算方法,其特征是,它包括建立二进制网络矩阵、建立交直流支路功率模型、建立功率平衡方程、求解功率平衡方程,具体内容有:
1)建立二进制网络矩阵:
①N行1列节点种类矩阵B,其中,N是节点数,N∈[2,9999];对于直流节点i∈[1,N],j∈[1,N],Bi=1,Bj=1,对于交流节点Bi=0,Bj=0;
②N行N列网络连接矩阵C,对于两不同节点i,j∈[1,N]之间有支路连接,Cij=1;对于两不同节点i,j之间无支路连接,Cij=0;
③N行N列支路种类矩阵L,对于两不同节点i,j之间有直流支路连接,Lij=1;对于两不同节点i,j之间有交流支路连接,Lij=0;
2)建立交直流支路功率模型:
①两个交流节点i,j相连的功率方程,
其中,Pij、Qij分别是i、j两节点间交流支路首端环节输入的有功功率和无功功率;θij是i、j两节点间电压相角差;Vi、Vj分别为i、j两节点电压幅值;Gij、Bij分别是i、j两节点间交流支路电导和电纳,结合网络矩阵,得到的潮流方程为:
其中,Pij (a)、Qij (a)是中间变量;
②两个直流节点i,j相连的其功率方程为,
Pij=Gij(Vi 2-ViVj) (3)
结合网络矩阵,得到的潮流方程为:
其中,Pij (b)是中间变量;
则:
Vi=MikVk (5)
Pkj=VkIkj=Vk(Vk-Vj)Gkj (6)
其中,Vk是交流-直流换流器交流侧节点i在直流侧等效节点k节点电压幅值,Mik是交流-直流变流器交直流两侧i、k节点电压幅值间折算系数,Gkj是k、j节点间连接的等效电导参数,Ikj是k节点流向j节点的电流,结合网络矩阵,得到潮流方程,
④一个直流节点i和一个交流节点j相连,功率方程为
结合网络矩阵,得到潮流方程,
其中,Pij (d)是中间变量;
⑤两个交流节点i、j通过直流线路相连,列功率方程和潮流方程:
其中,k、k’节点分别是i、j节点交流-直流变流器直流侧对应节点,a2、b2、Pij (e)、Qij (e)是中间变量,
且
3)建立功率平衡方程:
建立交直流混合配电网功率平衡方程
其中:Pis、Qis分别是i节点注入的有功功率和无功功率,且
4)应用牛顿-拉夫逊法求解功率平衡方程:
①求出修正方程的不平衡量;
②求出修正方程的系数矩阵雅可比矩阵;
③应用高斯消去法求解修正方程,得到修正量;
④求出新值;
⑤判断修正量精度,若不满足自第一步迭代计算,若满足根据电功率公式求潮流分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010042263.4A CN111244962B (zh) | 2020-01-15 | 2020-01-15 | 一种统一求解的交直流混合配电网潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010042263.4A CN111244962B (zh) | 2020-01-15 | 2020-01-15 | 一种统一求解的交直流混合配电网潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111244962A CN111244962A (zh) | 2020-06-05 |
CN111244962B true CN111244962B (zh) | 2022-06-28 |
Family
ID=70879529
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010042263.4A Expired - Fee Related CN111244962B (zh) | 2020-01-15 | 2020-01-15 | 一种统一求解的交直流混合配电网潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111244962B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111884268B (zh) * | 2020-07-14 | 2021-10-01 | 东北电力大学 | 考虑可靠性指标的交直流混合配电网最优潮流控制方法 |
CN111799800A (zh) * | 2020-07-21 | 2020-10-20 | 国网上海市电力公司 | 一种交直流混合配电网潮流计算方法 |
CN112909952B (zh) * | 2021-03-09 | 2022-10-04 | 武汉船用电力推进装置研究所(中国船舶重工集团公司第七一二研究所) | 一种舰船直流综合电力推进系统环网交直流混合潮流计算方法 |
CN113517713B (zh) * | 2021-09-13 | 2021-12-31 | 国网江西省电力有限公司电力科学研究院 | 适用于交直流混联系统的静态电压安全域分析方法及装置 |
CN115313529B (zh) * | 2022-08-05 | 2023-09-12 | 国网安徽省电力有限公司经济技术研究院 | 一种考虑空间特性和交直流两侧耦合作用的电能质量评估方法 |
CN115036962B (zh) * | 2022-08-15 | 2022-11-04 | 山东大学 | 柔性直流输电系统潮流计算及交直流混联潮流计算方法 |
CN115995818A (zh) * | 2023-03-22 | 2023-04-21 | 国网湖北省电力有限公司经济技术研究院 | 一种多能互补交直流混合配电网潮流计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105958485A (zh) * | 2016-06-14 | 2016-09-21 | 中国电力科学研究院 | 一种柔性互联交直流混合配电网的潮流计算方法 |
CN106953329A (zh) * | 2017-04-11 | 2017-07-14 | 国网安徽省电力公司电力科学研究院 | 含直流系统的复杂拓扑有源配电网潮流计算方法 |
CN108988316A (zh) * | 2018-06-15 | 2018-12-11 | 四川大学 | 一种交直流混合配电系统网架结构优化配置方法 |
CN109494746A (zh) * | 2018-11-08 | 2019-03-19 | 国网甘肃省电力公司电力科学研究院 | 基于改进自适应下垂控制的孤岛交直流混联微电网潮流计算方法 |
CN110176760A (zh) * | 2019-04-03 | 2019-08-27 | 国网浙江省电力有限公司电力科学研究院 | 一种交直流混合网络配置规划模型 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10243836B2 (en) * | 2017-04-28 | 2019-03-26 | The Boeing Company | Communications traffic optimization and obfuscation in delay and disruption tolerant networks |
-
2020
- 2020-01-15 CN CN202010042263.4A patent/CN111244962B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105958485A (zh) * | 2016-06-14 | 2016-09-21 | 中国电力科学研究院 | 一种柔性互联交直流混合配电网的潮流计算方法 |
CN106953329A (zh) * | 2017-04-11 | 2017-07-14 | 国网安徽省电力公司电力科学研究院 | 含直流系统的复杂拓扑有源配电网潮流计算方法 |
CN108988316A (zh) * | 2018-06-15 | 2018-12-11 | 四川大学 | 一种交直流混合配电系统网架结构优化配置方法 |
CN109494746A (zh) * | 2018-11-08 | 2019-03-19 | 国网甘肃省电力公司电力科学研究院 | 基于改进自适应下垂控制的孤岛交直流混联微电网潮流计算方法 |
CN110176760A (zh) * | 2019-04-03 | 2019-08-27 | 国网浙江省电力有限公司电力科学研究院 | 一种交直流混合网络配置规划模型 |
Non-Patent Citations (2)
Title |
---|
A Planning Approach for the Network Configuration of AC–DC Hybrid Distribution Systems;Haytham M. A. Ahmed et al.;《IEEE TRANSACTIONS ON SMART GRID》;20180531;第9卷(第3期);第2203-2213页 * |
基于二进制粒子群算法的交直流混合配电网故障恢复方法;马天祥 等;《电力系统保护与控制》;20190501;第47卷(第9期);第112-119页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111244962A (zh) | 2020-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111244962B (zh) | 一种统一求解的交直流混合配电网潮流计算方法 | |
CN106410835A (zh) | 一种多端柔性直流输电系统的电网仿真实现方法 | |
Beerten et al. | Generalized steady-state VSC MTDC model for sequential AC/DC power flow algorithms | |
CN104638678B (zh) | 一种多风机并网谐振分析方法 | |
CN108134401B (zh) | 交直流混合系统多目标潮流优化及控制方法 | |
CN103123664B (zh) | 一种模块化多电平换流器动态模型建模方法 | |
US20060282239A1 (en) | Method of setting-up steady state model of VSC-based multi-terminal HVDC transmission system | |
CN106911132A (zh) | 基于mmc的分布式潮流控制器的电磁暂态数学模型、控制系统及建模方法 | |
CN102570477B (zh) | 基于线性化潮流算法模型的svc及tcsc混合布点规划方法 | |
CN105356481B (zh) | 一种基于多馈入短路比的动态无功补偿选点方法 | |
CN111245002A (zh) | 基于mmc的双极柔性直流电网短路和接地故障电流预测方法 | |
CN110489806B (zh) | 包含多电压源型变流器电网的电磁暂态建模和计算方法 | |
CN105958485A (zh) | 一种柔性互联交直流混合配电网的潮流计算方法 | |
CN107769213A (zh) | 一种多换流器并联的交直流配电网潮流计算方法 | |
CN103793755B (zh) | 基于原对偶内点法的混合直流输电系统最优潮流方法 | |
CN103558481A (zh) | 同塔并架多回直流线路耦合作用引发换相失败的判别方法 | |
CN109038641A (zh) | 高压直流输电系统快速提升功率能力的计算方法及系统 | |
CN109241678A (zh) | 模块化多电平换流器多拓扑通用仿真模型 | |
CN104899396A (zh) | 一种修正系数矩阵的快速分解法潮流计算方法 | |
CN108233408B (zh) | 一种mmc-mtdc系统自适应下垂控制方法 | |
CN104934973B (zh) | 一种考虑温度的电力线路损耗计算方法 | |
CN109861233B (zh) | 交直流混合配电网解耦潮流确定方法 | |
CN106229970A (zh) | 基于换流器控制特性的微电网状态估计方法 | |
CN111985066B (zh) | 适用于多电压等级直流配电系统的平均动态相量模型 | |
CN206004320U (zh) | 模块化多电平换流器的仿真装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220628 |
|
CF01 | Termination of patent right due to non-payment of annual fee |