CN107069739A - 一种基于前推回代法的配电网潮流计算方法 - Google Patents
一种基于前推回代法的配电网潮流计算方法 Download PDFInfo
- Publication number
- CN107069739A CN107069739A CN201710081025.2A CN201710081025A CN107069739A CN 107069739 A CN107069739 A CN 107069739A CN 201710081025 A CN201710081025 A CN 201710081025A CN 107069739 A CN107069739 A CN 107069739A
- Authority
- CN
- China
- Prior art keywords
- node
- voltage
- power
- search
- 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.)
- Withdrawn
Links
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
技术领域:
本发明涉及配电网分析、运行与控制技术领域,更具体的说是涉及一种配电网潮流计算方法。
背景技术:
配电网潮流计算是电网经济运行、***分析的重要基础。配电网不仅呈辐射状运行结构,而且分支多,各馈线之间基本没有联系,与输电网络结构有明显差异,正常运行的配电网具有辐射状网络结构、负荷节点数量很多、线路R/X较大等特点,所以传统的潮流计算方法如:牛顿法、PQ分解法等在配电网潮流计算中不再适用。配电网潮流计算问题在数学上是一组多元非线性方程式求解问题,其解法都离不开迭代。因此,对潮流计算方法,首先要求它能可靠地收敛,并给出正确答案。由于配电网***结构及参数的一些特点,并且随着配电网***不断扩大,潮流计算的方程式阶数也越来越高,对这样的方程式并不是任何数学方法都能保证给出正确答案的。这种情况成为促使配电网***计算人员不断寻求新的更可靠方法的重要因素。
近年来,许多学者对配电网潮流计算展开大量的研究,并出现了许多计算配电网潮流的算法,主要有:回路阻抗法,改进牛顿法,快速解耦法等。虽然有些学者为使快速解偶法能在配电网得以继续应用而做了一些有益的尝试,如应用补偿技术处理R/X较大的线路,但这些方法都使算法复杂化,丧失了快速解偶算法原有的计算量小,收敛可靠的特点。目前还没有一种很好的方法能够满足配电网潮流计算的各项要求,因此急需开发出一种适合配电网的潮流计算方法,因此,鉴于配电网潮流计算算法的现状予以改进。
发明内容:
本发明的目的就是针对现有技术之不足,而提供一种基于前推回代法的配电网潮流计算方法,其具有编程简单、计算速度快、收敛性好、数值稳定性好的特点,此计算方法是配电网潮流计算的有效算法,具有很强的实用性。
本发明的技术解决措施如下:
一种基于前推回代法的配电网潮流计算方法,包括以下步骤:
1)首先要搜索节点关系,确定拓扑结构表;
2)给定配电馈线根节点电压并为其它节点电压赋初值k=0;
3)考虑负荷电压静特性,计算负荷数据;
4)从各负荷节点出发,先子节点后父节点,通过前推计算,由节点电压分布求支路功率分布;
5)从节点出发,先父节点后子节点,通过回推计算,由支路功率分布求节点电压分布
6)判断相邻两次迭代电压差的模分量的最大值是否小于给定的收敛指标ε,若是,则停止计算,否则,k=k+1,转步2)。
作为上述技术方案的优选,所述的步骤a)中采用原始数据输入结构,自动搜索节点关系,确定网络的拓扑结构,具体步骤如下,
a)搜索末梢节点作为第一层节点;
b)搜索末梢节点的父节点作为第二层节点;
c)继续搜索第二层节点的父节点作为第三层节点,这样反的搜索下去,直到搜索到某层节点的父节点全部是根节点时停止搜索;
d)删除在后面层次中有重复的前面层次中的节点,形成真正的层次关系,确定潮流计算的节点顺序。
作为上述技术方案的优选,所述的原始数据输入结构为:
节点结构体:{节点号 节点有功 节点无功}
支路结构体:{支路首端节点号 支路末端节点号 支路电阻 支路电抗}。
作为上述技术方案的优选,所述的前推后代的计算方法包括如下步骤,
a)首先假定各节点的电压幅值为1,幅角为0;
b)从第一层节点开始,根据基尔霍夫电流定律,求支路上的电流:
式中,是节点j的功率,是节点j的电压;
c)从第二层开始逐层计算非末梢节点的注入电流:
d)由步骤b)和c)可求出所有支路的支电流,再利用已知的根节点电压,从根节点向后顺次求得各个负荷节点的电压
其中i为父节点,j为子节点,Z为i、j间支路的阻抗;
e)计算各节点的电压幅值修正量:
计算节点电压修正量的最大值;
f)判别收敛条件:
其中k为迭代次数,若最大电压修正量小于阈值,则跳出循环,输出电压计算结果;否则重复步骤b)~f),直到满足条件为止;
g)在得到各个节点的电压电流后,就可以计算线路潮流S:
本发明的有益效果在于:本发明以馈线为基本计算单位,假设电压为额定电压,根据负荷功率由末端向始端逐段推算,仅计算各元件中的功率损耗而不计算节点电压,求得各支路上的电流和功率损耗,并据此获得始端功率,再根据给定的始端电压和求得的始端功率,由始端向末端逐段推算电压降落,求得各节点电压,如此重复上述过程,直至各个节点的功率偏差满足允许条件为止。本发明具有编程简单、计算速度快、收敛性好、数值稳定性好的特点,此计算方法是配电网潮流计算的有效算法,具有很强的实用性。
附图说明:
图1为本发明配电网潮流计算方法的流程图。
具体实施方式:
实施例:以下由特定的具体实施例说明本发明的实施方式,本领域技术人员可由本说明书所揭示的内容轻易地了解本发明的其他优点及功效。
见图1所示,一种基于前推回代法的配电网潮流计算方法,包括以下步骤:
1)首先要搜索节点关系,确定拓扑结构表;
2)给定配电馈线根节点电压并为其它节点电压赋初值k=0;
3)考虑负荷电压静特性,计算负荷数据;
4)从各负荷节点出发,先子节点后父节点,通过前推计算,由节点电压分布求支路功率分布;
5)从节点出发,先父节点后子节点,通过回推计算,由支路功率分布求节点电压分布
6)判断相邻两次迭代电压差的模分量的最大值是否小于给定的收敛指标ε,若是,则停止计算,否则,k=k+1,转步骤2)。
首先搜索节点关系,确定拓扑结构表。为了配合算法和避免复杂的网络编号,采用以下原始数据输入结构,不用形成节点导纳矩阵,就可以自动搜索节点关系,确定网络的拓扑结构。
节点结构体:
{节点号 节点有功 节点无功}
支路结构体:
{支路首端节点号 支路末端节点号 支路电阻 支路电抗}
根据线路首末节点,就可以确定每个节点连接的节点及其关系,从而可以形成整体的树状的关系结构。为了形成层次关系,确定节点计算顺序,要利用网络拓扑结构,经过多次按层遍历的广度优先搜索,形成层次关系,确定前推后代潮流算法的节点计算顺序。具体方法如下:
a)搜索末梢节点作为第一层节点;
b)搜索末梢节点的父节点作为第二层节点;
c)继续搜索第二层节点的父节点作为第三层节点,这样反的搜索下去,直到搜索到某层节点的父节点全部是根节点时停止搜索;
d)删除在后面层次中有重复的前面层次中的节点,形成真正的层次关系,确定潮流计算的节点顺序。
前推回代法基于支路电流进行,具体计算步骤为:
a)首先假定各节点的电压幅值为1,幅角为0,
b)从第一层节点开始,根据基尔霍夫电流定律,求支路上的电流:
式中,是节点j的功率,是节点j的电压。
c)从第二层开始逐层计算非末梢节点的注入电流:
d)由步骤b)和c)可求出所有支路的支电流,再利用已知的根节点电压,从根节点向后顺次求得各个负荷节点的电压
其中i为父节点,j为子节点,Z为i、j间支路的阻抗。
e)计算各节点的电压幅值修正量:
计算节点电压修正量的最大值
f)判别收敛条件:
其中k为迭代次数,若最大电压修正量小于阈值,则跳出循环,输出电压计算结果;否则重复步骤b)~f),直到满足条件为止。
g)在得到各个节点的电压电流后,就可以计算线路潮流S:
从树的末梢节点出发,利用已知的负荷功率,逐一计算,即可求得根节点处的电流,再从根节点出发,对树进行后向遍历,可求各节点电压。这样就完成一次前推后代的计算,迭代重复进行,直至满足收敛标准为止。
所述实施例用以例示性说明本发明,而非用于限制本发明。任何本领域技术人员均可在不违背本发明的精神及范畴下,对所述实施例进行修改,因此本发明的权利保护范围,应如本发明的权利要求所列。
Claims (4)
1.一种基于前推回代法的配电网潮流计算方法,其特征在于,包括以下步骤:
1)首先要搜索节点关系,确定拓扑结构表;
2)给定配电馈线根节点电压并为其它节点电压赋初值k=0;
3)考虑负荷电压静特性,计算负荷数据;
4)从各负荷节点出发,先子节点后父节点,通过前推计算,由节点电压分布求支路功率分布;
5)从节点出发,先父节点后子节点,通过回推计算,由支路功率分布求节点电压分布
6)判断相邻两次迭代电压差的模分量的最大值是否小于给定的收敛指标ε,若是,则停止计算,否则,k=k+1,转步骤2)。
2.根据权利要求1所述一种基于前推回代法的配电网潮流计算方法,其特征在于:所述步骤a)中采用原始数据输入结构,自动搜索节点关系,确定网络的拓扑结构,具体步骤如下,
a)搜索末梢节点作为第一层节点;
b)搜索末梢节点的父节点作为第二层节点;
c)继续搜索第二层节点的父节点作为第三层节点,这样反的搜索下去,直到搜索到某层节点的父节点全部是根节点时停止搜索;
d)删除在后面层次中有重复的前面层次中的节点,形成真正的层次关系,确定潮流计算的节点顺序。
3.根据权利要求2所述一种基于前推回代法的配电网潮流计算方法,其特征在于:所述原始数据输入结构为:
节点结构体:{节点号节点有功节点无功}
支路结构体:{支路首端节点号支路末端节点号支路电阻支路电抗}。
4.根据权利要求1所述一种基于前推回代法的配电网潮流计算方法,其特征在于:所述前推后代的计算方法包括如下步骤,
a)首先假定各节点的电压幅值为1,幅角为0;
b)从第一层节点开始,根据基尔霍夫电流定律,求支路上的电流:
式中,是节点j的功率,是节点j的电压;
c)从第二层开始逐层计算非末梢节点的注入电流:
d)由步骤b)和c)可求出所有支路的支电流,再利用已知的根节点电压,从根节点向后顺次求得各个负荷节点的电压
其中i为父节点,j为子节点,Z为i、j间支路的阻抗;
e)计算各节点的电压幅值修正量:
计算节点电压修正量的最大值;
f)判别收敛条件:
其中k为迭代次数,若最大电压修正量小于阈值,则跳出循环,输出电压计算结果;否则重复步骤b)~f),直到满足条件为止;
g)在得到各个节点的电压电流后,就可以计算线路潮流S:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710081025.2A CN107069739A (zh) | 2017-02-15 | 2017-02-15 | 一种基于前推回代法的配电网潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710081025.2A CN107069739A (zh) | 2017-02-15 | 2017-02-15 | 一种基于前推回代法的配电网潮流计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107069739A true CN107069739A (zh) | 2017-08-18 |
Family
ID=59598506
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710081025.2A Withdrawn CN107069739A (zh) | 2017-02-15 | 2017-02-15 | 一种基于前推回代法的配电网潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107069739A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110988594A (zh) * | 2019-12-06 | 2020-04-10 | 国网上海市电力公司 | 一种配电网故障位置确定方法及确定装置 |
CN111555285A (zh) * | 2020-04-03 | 2020-08-18 | 浙江工业大学 | 分布式冷热电联供综合能源***能量流解耦分析计算方法 |
CN111835011A (zh) * | 2020-06-19 | 2020-10-27 | 国网河南省电力公司鹤壁供电公司 | 一种配电网潮流计算方法 |
CN112202177A (zh) * | 2020-09-29 | 2021-01-08 | 国网河南省电力公司洛阳供电公司 | 一种计算配电网支路***潮流的方法 |
CN112350318A (zh) * | 2020-11-16 | 2021-02-09 | 北京交通大学 | 基于广度优先搜索算法的交流配电网拓扑辨识方法 |
CN112510687A (zh) * | 2020-11-20 | 2021-03-16 | 广西电网有限责任公司电力科学研究院 | 一种基于最大最小负荷潮流的线损异常校核方法及*** |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103746388A (zh) * | 2014-01-16 | 2014-04-23 | 南京软核科技有限公司 | 配电网无功电压三级协调控制方法 |
-
2017
- 2017-02-15 CN CN201710081025.2A patent/CN107069739A/zh not_active Withdrawn
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103746388A (zh) * | 2014-01-16 | 2014-04-23 | 南京软核科技有限公司 | 配电网无功电压三级协调控制方法 |
Non-Patent Citations (2)
Title |
---|
廖强: "基于实时数据的配电网线损分析及管理***研究", 《中国优秀博士学位论文全文数据库工程科技Ⅱ辑》 * |
祖凯涛等: "配电网精细化线损计算及降损分析", 《陕西电力》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110988594A (zh) * | 2019-12-06 | 2020-04-10 | 国网上海市电力公司 | 一种配电网故障位置确定方法及确定装置 |
CN110988594B (zh) * | 2019-12-06 | 2022-04-08 | 国网上海市电力公司 | 一种配电网故障位置确定方法及确定装置 |
CN111555285A (zh) * | 2020-04-03 | 2020-08-18 | 浙江工业大学 | 分布式冷热电联供综合能源***能量流解耦分析计算方法 |
CN111555285B (zh) * | 2020-04-03 | 2022-06-17 | 浙江工业大学 | 分布式冷热电联供综合能源***能量流解耦分析计算方法 |
CN111835011A (zh) * | 2020-06-19 | 2020-10-27 | 国网河南省电力公司鹤壁供电公司 | 一种配电网潮流计算方法 |
CN112202177A (zh) * | 2020-09-29 | 2021-01-08 | 国网河南省电力公司洛阳供电公司 | 一种计算配电网支路***潮流的方法 |
CN112350318A (zh) * | 2020-11-16 | 2021-02-09 | 北京交通大学 | 基于广度优先搜索算法的交流配电网拓扑辨识方法 |
CN112510687A (zh) * | 2020-11-20 | 2021-03-16 | 广西电网有限责任公司电力科学研究院 | 一种基于最大最小负荷潮流的线损异常校核方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107069739A (zh) | 一种基于前推回代法的配电网潮流计算方法 | |
CN105656031B (zh) | 基于高斯混合分布特征的含风电电力***安全风险评估方法 | |
CN110212521B (zh) | 配电网线损计算方法及*** | |
CN103236691B (zh) | 基于复仿射数学理论的三相不平衡潮流计算方法 | |
Arboleya et al. | BFS algorithm for voltage-constrained meshed DC traction networks with nonsmooth voltage-dependent loads and generators | |
CN108599154A (zh) | 一种考虑不确定性预算的三相不平衡配电网鲁棒动态重构方法 | |
CN107591807B (zh) | 一种新能源接入下的输电网规划的优化方法 | |
CN106911140A (zh) | 一种储能规划方法 | |
CN106532711A (zh) | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 | |
CN104037764A (zh) | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 | |
CN105514971A (zh) | 一种适用于各种运行模式微电网的潮流计算方法 | |
CN104751006B (zh) | 一种计及变量相关性的概率潮流计算方法 | |
CN103795057A (zh) | 配电网潮流计算中基于搜索的配电网拓扑编号生成方法 | |
CN108092284A (zh) | 一种基于线性模型的三相不平衡智能配电网网络重构方法 | |
CN104113061A (zh) | 一种含分布式电源的配电网三相潮流计算方法 | |
CN105656084A (zh) | 一种计及新能源发电预测误差的改进随机潮流算法 | |
CN104037759B (zh) | 电力***注入转移分布因子的抗差最小二乘估计方法 | |
Liu et al. | Remote voltage control using the holomorphic embedding load flow method | |
CN107808200B (zh) | 配电网短时检修计划优化方法和装置 | |
CN107257130A (zh) | 基于区域量测解耦的低压配电网损耗计算方法 | |
CN106056466A (zh) | 基于FP‑growth算法的大电网关键线路识别方法 | |
CN106709243A (zh) | 含小阻抗支路电网的补偿法极坐标牛顿法潮流计算方法 | |
CN106532712B (zh) | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 | |
CN107968399A (zh) | 一种快速搜索静态电压稳定域边界的方法 | |
CN102780220A (zh) | 含有pv恒定型分布式电源的配电网潮流计算方法 |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20170818 |