CN112202177B - 一种计算配电网支路系统潮流的方法 - Google Patents

一种计算配电网支路系统潮流的方法 Download PDF

Info

Publication number
CN112202177B
CN112202177B CN202011048173.2A CN202011048173A CN112202177B CN 112202177 B CN112202177 B CN 112202177B CN 202011048173 A CN202011048173 A CN 202011048173A CN 112202177 B CN112202177 B CN 112202177B
Authority
CN
China
Prior art keywords
power
constant
load
node
phase
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.)
Active
Application number
CN202011048173.2A
Other languages
English (en)
Other versions
CN112202177A (zh
Inventor
周燕
李展高
李雅琳
皇甫武军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Luoyang Power Supply Co of State Grid Henan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Luoyang Power Supply Co of State Grid Henan Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Luoyang Power Supply Co of State Grid Henan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202011048173.2A priority Critical patent/CN112202177B/zh
Publication of CN112202177A publication Critical patent/CN112202177A/zh
Application granted granted Critical
Publication of CN112202177B publication Critical patent/CN112202177B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种计算配电网支路系统潮流的方法，通过采用收敛更快的前推回代法对支路系统的节点电压和电流进行迭代，从而计算出电网潮流，计算准确度高，计算速度块，适用于配电网支路系统，可为完善配电网的智能监测提供进一步的技术支撑。

Description

一种计算配电网支路系统潮流的方法
技术领域
本发明涉及电网监测及维护技术，具体是一种计算配电网支路系统潮流的方法。
背景技术
电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性、可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要的基础计算。
目前，PQ分解法是输电网中采用最多的潮流计算方法。该方法忽略相连的两节点间电压相角差，使雅克比矩阵变成常数阵。但是在配电网中，线路的R/X 值不能忽略，这使PQ分解法所假设的前提条件不成立。另外配电网还有很多区别于输电网的特性：结构成环状，运行时呈开环，稳态运行时拓扑结构呈辐射型、三相线路包括阻抗等参数不对称、三相负荷的不平衡，抑或称之为功率损耗的不平衡问题明显等。因此，以PQ分解法计算配电网潮流是不合适的。
发明内容
针对背景技术中的问题，本发明的目的是提出一种计算配电网支路系统潮流的方法，采用收敛更快的前推回代法来进行潮流计算，计算准确度高，计算速度块，适用于配电网支路系统，可为完善配电网的智能监测提供进一步的技术支撑。
为实现上述目的，本发明采用以下技术方案：
一种计算配电网支路系统潮流的方法，包括以下步骤：
步骤1、获取系统节点与支路数据以及各个节点的电压初值UK；节点与支路数据包括节点数量、支路数量、每个支路的连接方式、每个支路上的元件及负荷数量；
步骤2、根据节点注入功率与节点电压初值UK计算所有支路的末节点注入电流IK；
步骤3、自各支路末端节点向首端节点计算当前电压迭代值UK之下的系统各支路电流IK+1，此时仅计算各元件中的功率损耗而不计算电压降落；
步骤4、不计算功率损耗，自各支路首端节点向末端节点根据线路阻抗与线路电流推算各个支路电压压降U；
步骤5、根据电压压降U计算出本次迭代各个节点电压UK+1；
步骤6、设△U＝∣UK+1-UK|，△I＝∣IK+1-IK|，并设定一个阈值ε，当max(△ U、△I)≥ε时，设定各个节点上一次计算出的节点电压为本次迭代计算的节点电压初值，并转入步骤2；否则转入下一步；
步骤7：计算结束。
2、根据权利要求1所述的一种计算配电网支路系统潮流的方法，其特征是：所述步骤2中，注入电流IK的计算方法为：根据负荷的连接方式及负荷类型，列写相应的负荷参数方程，并对方程求解，从而得到注入电流IK；
其中，负荷的连接方式包括星形连接和三角形连接，负荷类型包括恒功率负荷、恒电流负荷和恒阻抗负荷；
对于星形连接的恒功率负荷，其恒功率参数方程为：
Figure BDA0002708675900000031
上式中，
Figure BDA0002708675900000032
为恒功率阵；
Figure BDA0002708675900000033
为A相恒功率，
Figure BDA0002708675900000034
为B相恒功率，
Figure BDA0002708675900000035
为C相恒功率；
对于星形连接的恒电流负荷，其恒电流参数方程为：
Figure BDA0002708675900000036
上式中，
Figure BDA0002708675900000037
为A相电压，
Figure BDA0002708675900000038
为B相电压，
Figure BDA0002708675900000039
为C相电压；
Figure BDA00027086759000000310
为恒电流，IK为支路末端注入电流；
对于星形连接的恒阻抗负荷，其恒阻抗参数方程为：
Figure BDA00027086759000000311
上式中，U为支路电压压降，
Figure BDA00027086759000000312
为恒阻抗的倒数，
Figure BDA00027086759000000313
为A相恒功率的标幺值，
Figure BDA00027086759000000314
为A相恒功率的标幺值，
Figure BDA00027086759000000315
为C相恒功率的标幺值；
对于三角形连接的负荷，首先设定其相间功率为：
Figure BDA00027086759000000316
上式中，
Figure BDA00027086759000000317
为AB相间功率，
Figure BDA00027086759000000318
为BC相间功率，
Figure BDA00027086759000000319
为CA相间功率；那么对于三角形连接的恒功率负荷，其恒功率参数方程为：
Figure BDA0002708675900000041
对于三角形连接的恒电流负荷，其恒电流参数方程为：
Figure BDA0002708675900000042
上式中，*号表示取标幺值；
对于三角形连接的恒阻抗负荷，其恒阻抗参数方程为：
Figure BDA0002708675900000043
上式中，*号表示取标幺值；
分别对上述方程式(1)～(7)求解，即可求出对应节点的注入电流IK。
本发明的有益效果是：本发明综合了支路电流法和前推回代法对电网潮流进行计算，算法收敛性能好，计算精度高，适用于配电网支路系统，可为完善配电网的智能监测提供进一步的技术支撑。
附图说明
图1为本发明的流程图。
图2为支路的统一模型示意图。
具体实施方式
下面结合说明书附图对本发明作进一步的详细说明。
配电网支路系统的统一模型示意图如图2所示，定义图中i为首节点，j为到根节点距离比i大1的末节点。
计算该支路系统潮流的方法，包括以下步骤：
步骤1、获取系统节点与支路数据以及各个节点的电压初值UK；节点与支路数据包括节点数量、支路数量、每个支路的连接方式、每个支路上的负荷数量；
步骤2、根据节点注入功率与节点电压初值UK计算所有支路的末节点注入电流IK；所述注入电流IK的计算方法为：根据负荷的连接方式及负荷类型，列写相应的负荷参数方程，并对方程求解，从而得到注入电流IK；
负荷的连接方式包括星形连接和三角形连接，负荷类型包括恒功率负荷、恒电流负荷和恒阻抗负荷；其中，
对于星形连接的恒功率负荷，其恒功率参数方程为：
Figure BDA0002708675900000051
上式中，
Figure BDA0002708675900000052
为恒功率阵；
Figure BDA0002708675900000053
为A相恒功率，
Figure BDA0002708675900000054
为B相恒功率，
Figure BDA0002708675900000055
为C相恒功率；
对于星形连接的恒电流负荷，其恒电流参数方程为：
Figure BDA0002708675900000056
上式中，
Figure BDA0002708675900000057
为A相电压，
Figure BDA0002708675900000058
为B相电压，
Figure BDA0002708675900000059
为C相电压；
Figure BDA00027086759000000510
为恒电流，IK为支路末端注入电流；
对于星形连接的恒阻抗负荷，其恒阻抗参数方程为：
Figure BDA0002708675900000061
上式中，U为支路电压压降，
Figure BDA0002708675900000062
为恒阻抗的倒数，
Figure BDA0002708675900000063
为A相恒功率的标幺值，
Figure BDA0002708675900000064
为A相恒功率的标幺值，
Figure BDA0002708675900000065
为C相恒功率的标幺值；
对于三角形连接的负荷，首先设定其相间功率为：
Figure BDA0002708675900000066
上式中，
Figure BDA0002708675900000067
为AB相间功率，
Figure BDA0002708675900000068
为BC相间功率，
Figure BDA0002708675900000069
为CA相间功率；
那么对于三角形连接的恒功率负荷，其恒功率参数方程为：
Figure BDA00027086759000000610
对于三角形连接的恒电流负荷，其恒电流参数方程为：
Figure BDA00027086759000000611
上式中，*号表示取标幺值；
对于三角形连接的恒阻抗负荷，其恒阻抗参数方程为：
Figure BDA0002708675900000071
上式中，*号表示取标幺值；
分别对上述方程式(1)～(7)求解，即可求出对应节点的注入电流IK；
步骤3、自各支路末端节点向首端节点计算当前电压迭代值UK之下的系统各支路电流IK+1，此时仅计算各元件中的功率损耗而不计算电压降落，计算公式为：
SK＝UK·(IK+1)*
上式中，
Figure BDA0002708675900000072
为有效功率，(IK＝1)*为IK+1的标幺值；
步骤4、不考虑功率损耗，自各支路首端节点向末端节点根据线路阻抗与线路电流推算各个支路电压压降U，计算公式为：
U＝Z·IK+1
上式中，Z为线路阻抗；
步骤5、根据电压压降U计算出本次迭代各个节点电压UK+1；
步骤6、设△U＝∣UK+1-UK|，△I＝∣IK+1-IK|，并设定一个阈值ε，当max(△ U、△I)≥ε时，设定各个节点上一次计算出的节点电压为本次迭代计算的节点电压初值，并转入步骤2；否则转入下一步；
步骤7：计算结束。
在步骤3中，由末端节点向首端节点推算时，仅计算各元件中的功率损耗而不计算电压降落；待求得始端功率后，再运用给定的始端电压和求得的始端功率由始端向末端逐段推算电压降落，但这时不再重新计算功率损耗。
本发明未详述部分为现有技术。

Claims (2)

1.一种计算配电网支路系统潮流的方法，其特征是：包括以下步骤：
步骤1、获取系统节点与支路数据以及各个节点的电压初值UK；节点与支路数据包括节点数量、支路数量、每个支路的连接方式、每个支路上的元件及负荷数量；
步骤2、根据节点注入功率与节点电压初值UK计算所有支路的末节点注入电流IK；
步骤3、自各支路末端节点向首端节点计算当前电压迭代值UK之下的系统各支路电流IK+1，此时仅计算各元件中的功率损耗而不计算电压降落；
步骤4、不计算功率损耗，自各支路首端节点向末端节点根据线路阻抗与线路电流推算各个支路电压压降U；
步骤5、根据电压压降U计算出本次迭代各个节点电压UK+1；
步骤6、设△U＝∣UK+1-UK|，△I＝∣IK+1-IK|，并设定一个阈值ε，当max(△U、△I)≥ε时，设定各个节点上一次计算出的节点电压为本次迭代计算的节点电压初值，并转入步骤2；否则转入下一步；
步骤7：计算结束。
2.根据权利要求1所述的一种计算配电网支路系统潮流的方法，其特征是：所述步骤2中，注入电流IK的计算方法为：根据负荷的连接方式及负荷类型，列写相应的负荷参数方程，并对方程求解，从而得到注入电流IK；
其中，负荷的连接方式包括星形连接和三角形连接，负荷类型包括恒功率负荷、恒电流负荷和恒阻抗负荷；
对于星形连接的恒功率负荷，其恒功率参数方程为：
Figure FDA0002708675890000021
上式中，
Figure FDA0002708675890000022
为恒功率阵；
Figure FDA0002708675890000023
为A相恒功率，
Figure FDA0002708675890000024
为B相恒功率，
Figure FDA0002708675890000025
为C相恒功率；
对于星形连接的恒电流负荷，其恒电流参数方程为：
Figure FDA0002708675890000026
上式中，
Figure FDA0002708675890000027
为A相电压，
Figure FDA0002708675890000028
为B相电压，
Figure FDA0002708675890000029
为C相电压；
Figure FDA00027086758900000210
为恒电流，IK为支路末端注入电流；
对于星形连接的恒阻抗负荷，其恒阻抗参数方程为：
Figure FDA00027086758900000211
上式中，U为支路电压压降，
Figure FDA00027086758900000212
为恒阻抗的倒数，
Figure FDA00027086758900000213
为A相恒功率的标幺值，
Figure FDA00027086758900000214
为A相恒功率的标幺值，
Figure FDA00027086758900000215
为C相恒功率的标幺值；
对于三角形连接的负荷，首先设定其相间功率为：
Figure FDA00027086758900000216
上式中，
Figure FDA00027086758900000217
为AB相间功率，
Figure FDA00027086758900000218
为BC相间功率，
Figure FDA00027086758900000219
为CA相间功率；
那么对于三角形连接的恒功率负荷，其恒功率参数方程为：
Figure FDA0002708675890000031
对于三角形连接的恒电流负荷，其恒电流参数方程为：
Figure FDA0002708675890000032
上式中，*号表示取标幺值；
对于三角形连接的恒阻抗负荷，其恒阻抗参数方程为：
Figure FDA0002708675890000033
上式中，*号表示取标幺值；
分别对上述方程式(1)～(7)求解，即可求出对应节点的注入电流IK。
CN202011048173.2A 2020-09-29 2020-09-29 一种计算配电网支路系统潮流的方法 Active CN112202177B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011048173.2A CN112202177B (zh) 2020-09-29 2020-09-29 一种计算配电网支路系统潮流的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011048173.2A CN112202177B (zh) 2020-09-29 2020-09-29 一种计算配电网支路系统潮流的方法

Publications (2)

Publication Number Publication Date
CN112202177A CN112202177A (zh) 2021-01-08
CN112202177B true CN112202177B (zh) 2022-03-04

Family

ID=74007825

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011048173.2A Active CN112202177B (zh) 2020-09-29 2020-09-29 一种计算配电网支路系统潮流的方法

Country Status (1)

Country Link
CN (1) CN112202177B (zh)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3668051B2 (ja) * 1999-06-04 2005-07-06 ティーエム・ティーアンドディー株式会社 充停電判定および潮流計算における系統認識方法
JP4556946B2 (ja) * 2006-12-19 2010-10-06 株式会社日立製作所 潮流計算方法及びその装置
CN102447254A (zh) * 2010-10-09 2012-05-09 江苏方程电力科技有限公司 基于补偿电流环状配电网潮流方法
CN103248044A (zh) * 2013-05-07 2013-08-14 广东电网公司电力科学研究院 配电网潮流计算方法与系统
CN107069739A (zh) * 2017-02-15 2017-08-18 河海大学 一种基于前推回代法的配电网潮流计算方法
CN109802392B (zh) * 2019-02-26 2020-10-13 中国农业大学 大规模配电网潮流计算方法及装置
CN111555273A (zh) * 2020-04-30 2020-08-18 华北电力大学 基于多直流弱送端电网连锁故障传播的无功电压分析方法

Also Published As

Publication number Publication date
CN112202177A (zh) 2021-01-08

Similar Documents

Publication Publication Date Title
CN108363306B (zh) 基于线性二次型优化的微电网分布式控制器参数确定方法
WO2018103317A1 (zh) 一种含upfc电力系统的通用潮流计算方法
CN111082427A (zh) 一种基于全纯函数的微电网潮流计算方法
CN109167366A (zh) 一种电力系统静态电压稳定临界点的计算方法
CN110401184A (zh) 多直流馈入受端电网紧急控制优化方法及系统
CN117973096B (zh) 基于双层分网的直驱风电场并行仿真方法、系统、终端及介质
CN112255567A (zh) 一种含光伏电源配电网的短路电流快速确定方法
CN109149584B (zh) 提高ipfc功率注入模型在潮流计算中收敛性的方法
CN112202177B (zh) 一种计算配电网支路系统潮流的方法
CN104934973A (zh) 一种考虑温度的电力线路损耗计算方法
CN108306299B (zh) 一种配电网异步迭代分布式三相潮流计算方法及系统
CN106991229B (zh) 一种针对复杂拓扑的风电场等值建模方法
CN112531715B (zh) 基于虚拟电阻的下垂控制多端直流微电网潮流计算方法
CN117458534A (zh) 一种新型液流储能调峰调频方法及装置
CN102593876A (zh) 一种大规模交直流互联电力系统的连续潮流算法
CN114629122A (zh) 一种交直流混合配电网线性潮流计算方法
CN116316852A (zh) 一种多光伏电源故障穿越协调方法
CN110456223B (zh) 一种含分布式电源和电动汽车的配电网短路电流测算方法
CN113922384A (zh) 一种风电场分散式无功电压优化协调控制方法
CN113659559A (zh) 一种用于dfig并网电力系统暂态稳定分析的方法
Kumar et al. A load flow algorithm for radial systems having renewable energy generation sources
CN115800285A (zh) 一种含环网有源配电网的快速潮流计算方法
CN114912285B (zh) 一种适用于大型风电场动态稳定性分析的降阶方法及系统
CN112198454B (zh) 一种在含新能源电源的网络中进行短路电流局部迭代的方法
CN116760039B (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