CN102891485A - 基于序分量法的弱环配电网三相解耦潮流计算方法 - Google Patents
基于序分量法的弱环配电网三相解耦潮流计算方法 Download PDFInfo
- Publication number
- CN102891485A CN102891485A CN2012103890561A CN201210389056A CN102891485A CN 102891485 A CN102891485 A CN 102891485A CN 2012103890561 A CN2012103890561 A CN 2012103890561A CN 201210389056 A CN201210389056 A CN 201210389056A CN 102891485 A CN102891485 A CN 102891485A
- Authority
- CN
- China
- Prior art keywords
- node
- sequence
- phase
- distribution network
- power distribution
- 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.)
- Granted
Links
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于序分量法的弱环配电网三相解耦潮流计算方法,包括以下步骤:(1)基于序分量法对弱环配电网进行三相潮流建模,并实现三序解耦;(2)在步骤(1)的基础上,建立弱环配电网三相解耦潮流计算方法,所述弱环配电网三相解耦潮流算法是基于道路的回路分析法。因此,该算法保留了面向支路的前推回推法的优点,而其收敛性不受回路多少的影响,具有很好的收敛性。同时,由于采用三序解耦,减少***潮流计算节点数目,降低了潮流计算量,提高了计算效率。最后,通过6母线和33母线测试算例验证了本发明方法的正确性和良好的收敛性,以及较强的回路处理能力,迭代次数随着闭合环路的增加而减少。
Description
技术领域
本发明属于电力***分析与计算领域,具体涉及一种基于序分量法的弱环配电网三相解耦潮流计算方法。
背景技术
配电网潮流计算是配电***分析的一项重要内容,它是对配电***规划设计和运行方式的合理性、可靠性及经济性进行定量分析的重要依据。国内外学者根据配电网的特点提出了各种配电网潮流算法,如前推回代法、隐式Zbus高斯法、回路阻抗法、改进牛顿法及快速解耦法等。其中前推回代法由于充分利用了配电网的结构特点,并且其具有物理概念明晰、编程简单、没有大矩阵计算、计算速度快、收敛性较好、非常适合于求解辐射状配电网潮流等优点而被广泛应用。但这些方法处理环网能力较弱。
相分量法对三相直接建模,物理概念清晰,易处理不对称元件和负荷。但是,对称元件三相不能解耦,导致潮流计算时节点数目较多,计算量较大,计算效率不高。同时,平衡节点和PV节点的处理也存在困难。
因此,需要一种新的弱环配电网三相潮流计算方法以解决上述问题。
发明内容
发明目的:本发明针对现有技术中配电网三相潮流计算方法存在的缺陷,提供一种基于序分量法的弱环配电网三相解耦潮流计算方法。
技术方案:为解决上述技术问题,本发明的基于序分量法的弱环配电网三相解耦潮流计算方法采用如下技术方案:
一种基于序分量法的弱环配电网三相解耦潮流计算方法,其特征在于,包括以下步骤:
(1)基于序分量法对弱环配电网进行三相潮流建模,并实现三序解耦;
(2)在步骤(1)的基础上,建立弱环配电网三相解耦潮流计算方法,所述弱环配电网三相解耦潮流算法是基于道路的回路分析法。
发明原理:采用序分量法时,可以对三相对称元件实现三序解耦,从而减少***潮流计算节点数目,提高计算效率。因此基于序分量法的三相潮流建模方法在处理线路空间对称、负荷三相不平衡的配电网具有很好的效果和计算优势。可以将一组不对称的“a”、“b”、“c”三相分量分解为三组三相对称的序分量,那么,三相潮流计算就变成了计算三组三相对称序分量中的一相。因此,配电网三相不平衡潮流计算的计算量会减少2/3,在保持较好的收敛性的情况下,可以给配电网三相潮流计算带来更快的计算速度。本发明针对现有配电网三相潮流计算方法存在的一些不足,提出了一种基于序分量法的弱环配电网三相解耦潮流计算方法,该方法首先基于序分量法进行配电***三相潮流建模,实现三序解耦。
有益效果:本发明针对现有配电网三相潮流计算方法存在的一些不足,提出了一种基于序分量法的弱环配电网三相解耦潮流计算方法,该方法首先基于序分量法进行配电***三相潮流建模,实现三序解耦。接着,在配电网三序解耦的情况下,基于道路的回路分析法提出一种改进的弱环配电网三序解耦潮流算法。该算法保留了面向支路的前推回推法的优点,而其收敛性不受回路多少的影响,具有很好的收敛性。同时,由于采用三序解耦,减少***潮流计算节点数目,降低了潮流计算量,提高了计算效率。最后,通过6母线和33母线测试算例验证了本发明方法的正确性和良好的收敛性,以及较强的回路处理能力,迭代次数随着闭合环路的增加而减少。
附图说明
图1为本发明流程结构示意图;
图2为6母线三相不平衡弱环配电网***;
图3为33母线三相不平衡弱环配电网***;
图4为算例2在表4中的7种运行方式下各节点C相电压幅值分布图。
具体实施方式
下面结合附图、具体实施方式和算例,进一步阐明本发明,应理解这些算例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
请参阅图1所示,本发明的基于序分量法的弱环配电网三相解耦潮流计算方法,包括以下步骤:(1)基于序分量法对弱环配电网进行三相潮流建模,并实现三序解耦;(2)在步骤(1)的基础上,建立弱环配电网三相解耦潮流计算方法,弱环配电网三相解耦潮流算法基于道路的回路分析法。
其中,步骤(1)弱环配电网三相潮流建模包括以下步骤:
A、确定弱环配电网络的树,树具有节点、树支和连支,分别给节点、树支和连支进行编号;
B、确定弱环配电网拓扑结构参数,弱环配电网拓扑结构参数包括节点数N、支路数b和独立回路数l;根节点是电源且作为参考节点,独立节点个数为n=N-1,支路数b=n+l,其中,弱环配电网络的树包括树支和连支,树支条数与独立节点个数相等,连支条数与独立回路数相等,即树支n条,连支l条,支路数等于树支条数与连支条数的和。
步骤(2)弱环配电网三相解耦潮流计算方法包括以下步骤:
E、计算三序网络道路矩阵Ts(n×b阶)和回-支关联矩阵Bs(l×b阶),另对于注入序电流为零的节点,在三序网络道路矩阵Ts中把该节点所对应行删去后形成新矩阵为Tsg(m×b阶),其中,下标s=0、1、2,分别表示三序网络模型中的零序、正序和负序网络模型,m为节点注入序电流不为零的节点个数。
J、计算k次迭代时的其中,为k次迭代时去除注入序电流为零的节点后形成的新注入序电流矩阵(m×1阶),m为节点注入序电流不为零的节点个数;为回路序阻抗矩阵的逆矩阵,即且s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型;
K、计算k次迭代时的其中,1n=[1,1,…,1]T,共n个1;s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型;
L、基于逆变换计算k次迭代时节点i三相电压相量i=1,2,…,n;
请参阅图1所示,本发明的基于序分量法的弱环配电网三相解耦潮流计算方法,具体包括如下步骤:
1)确定弱环配电网络的树,给节点编号,规定树中的根节点的编号最小,设为“0”,其余节点按其离根节点的远近来编号,离根节点越远的节点编号越大。而树支的编号则规定为取两端节点编号中的大者,连支的编号依次从树支的最大编号后加1。
2)确定弱环配电网拓扑结构参数,包括节点数,支路数,独立回路数。不妨设三相弱环配电网有N个节点、b条支路和l个独立回路,假设首节点(即根节点)是电源且作为参考节点,则独立节点个数为n=N-1,支路数b=n+l。其中,弱环配电网络的树包括树支和连支,树支条数与独立节点个数相等,连支条数与独立回路数相等,即,树支n条,连支l条;支路数等于树支条数与连支条数的和。
5)计算三序网络道路矩阵Ts(n×b阶)和回-支关联矩阵Bs(l×b阶);另对于注入序电流为零的节点,在三序网络道路矩阵Ts中把该节点所对应行删去后形成新矩阵为Tsg(m×b阶)。其中,下标s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型,m为节点注入序电流不为零的节点个数。
7)给弱环配电网各节点三相电压赋初始值其中,En=[E,E,…,E]T,共n个E,E为3×3单位矩阵。
10)计算k次迭代时的其中,为k次迭代时去除注入序电流为零的节点后形成的新注入序电流矩阵(m×1阶),m为节点注入序电流不为零的节点个数;为回路序阻抗矩阵的逆矩阵,即且s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型。
步骤10)的公式推导如下:
针对具有N个节点、b条支路和l个独立回路的三相弱环配电网,假设首节点(根节点)是电源且作为参考节点,则独立节点个数为n=N-1,支路数b=n+l。其中,树支n条,连支l条。
对于连通图中一颗选定的树,由于基本回路中仅包含一条连支,基本回路数等于连支数,基本回路-支路(下面简称“回-支”)关联矩阵B描述基本回路、树支、连支之间的联系。其中回-支关联矩阵B是一个l×b阶矩阵,假定连支支路的正方向都是从大号节点指向小号节点,基本回路的正方向与连支支路的正方向相同,如果支路k在回路j内,且二者方向相同,则B(j,k)=1,如果支路k在回路j内,且二者方向相反,则B(j,k)=-1,如果支路k不在回路j内,则B(j,k)=0。
一个节点的道路是指节点沿树到根所经过的路径上的支路集合,节点的道路强调的是路径上的支路,对于一个给定的树,节点的道路是唯一的,节点的道路只由树支支路组成,用道路矩阵T描述道路。其中道路矩阵T是一个n×b阶矩阵,假定道路的正方向都是从电源点(即根节点)指向各节点,各树支支路正方向与道路正方向相同,如果支路k在道路i上,则T(i,k)=1,反之T(i,k)=0。道路矩阵T是一个稀疏矩阵,利用稀疏技术可以降低内存需求。
在弱环配电序网络中,定义节点注入序电流向量矩阵为支路序电流向量矩阵为连支序电流(也即回路序电流)向量矩阵为在解耦的各序网模型电路中,可以获得各序网络的道路矩阵和回-支关联矩阵分别为T0、T1、T2和B0、B1、B2,并依据KCL电流定律,支路序电流与节点注入序电流回路序电流满足如下等式:
其中,s=0,1,2,分别表示序网络模型中的零序、正序和负序网络模型。上式表明支路序电流由两部分组成:一是节点注入序电流对在道路上的支路的贡献,用表示,二是回路序电流对在回路上的支路的贡献,用表示,则有
对任一弱环配电***序分量电路模型中,基于欧姆定律,支路特性约束为
依据KVL电压定律,基本回路约束为
由式(2)、(4)、(5)、(6)代入式(7)可得:
设电源节点三相电压相量矩阵为各节点三相电压相量矩阵为在配电序网络中,可以得出电源节点的三序电压矩阵为 各节点三序电压矩阵为那么,在各序网络模型中,可知任一节点与电源节点的序电压差等于从此节点开始沿着该节点的道路到达电源节点所经支路的支路序电压之和,即(设1n=[1,1,…,1]T,共n个1;s=0,1,2,分别表示序网络模型中的零序、正序和负序网络模型):
算例分析
算例1
图2为6母线三相不平衡弱环配电网,变压器Yn-yn接线方式,对***进行了一些调整,线路空间对称,即线路相分量阻抗矩阵完全对称,而三相负荷不平衡,其单线图如图2所示。
下面就开环和闭环2种运行方式进行讨论,2种运行方式下,基于本文算法的潮流收敛情况如表1所示,其中收敛精度为10-6。
表1算例1的2种运行方式及潮流收敛迭代次数
从表1可以看出,方式2的潮流计算比方式1的收敛更快,这是因为:在图1的***中,环路闭合后负荷离电源的电气距离缩短,负荷电流的变化引起负荷节点的电压波动变小,所以算法的迭代次数减少了,程序收敛更快了。
在方式1情况下,基于本发明计算方法的潮流计算结果如表2所示。
表2在方式1情况下潮流计算结果
在方式2情况下,基于本发明计算方法的潮流计算结果如表3所示。
表3在方式2情况下潮流计算结果
比较表2(方式1)和表3(方式2)中对应节点的三相电压也可以看出,测试***在闭环运行(即方式2)时,整个***的三相电压水平都得到了较大的提升,因为环路闭合后改善了整个***的功率流向,对电压有很好的提升作用,所以***的节点电压水平都得到了极大的改善,同时,也提高了算法的收敛性。这些在算例1的仿真中都已经得到了验证。
算例2
图3为33母线***,对负荷参数进行了调整,扩展为三相,但三相负荷不平衡,假定线路空间对称,即线路相分量阻抗矩阵完全对称,5条环路,其单线图如图3所示。
下面就表4所示的7种运行方式进行讨论,在7种运行方式下,各节点C相电压幅值分布如图4所示。
表4算例2的7种运行方式及潮流收敛迭代次数
由表4可以看出,Case1时,5条环路全部打开,算法的迭代次数为7次,迭代次数最多,收敛最慢。当Case7时,5条环路全部闭合,算法的迭代次数为4次,迭代次数最少,收敛最快。就算法收敛性的总体趋势来看,算法的迭代次数随着闭合环路的增加而减少。
图4为算例2在表4中的7种运行方式下各节点C相电压幅值分布图,从图4可以看出,Case1时,5条环路全部打开,节点电压分布最差,这也充分的反映了配电***负荷三相不平衡的特性,同时,也表明了表4中迭代次数最多的原因所在。随着环路的闭合,节点电压水平都得到了提高。特别是在Case7时,5条环路全部闭合,节点电压水平提高到最好。这说明:随着环路闭合的越多,由于***负荷三相不平衡给配电网带来的节点电压水平低问题都会得到更好的改善和提升,而且算法的收敛性也得到了提高。
Claims (7)
1.一种基于序分量法的弱环配电网三相解耦潮流计算方法,其特征在于,包括以下步骤:
(1)基于序分量法对弱环配电网进行三相潮流建模,并实现三序解耦;
(2)在步骤(1)的基础上,建立弱环配电网三相解耦潮流计算方法,所述弱环配电网三相解耦潮流算法是基于道路的回路分析法。
3.如权利要求2所述的基于序分量法的弱环配电网三相解耦潮流计算方法,其特征在于,所述电源节点是根节点且作为参考节点,独立节点个数为n=N-1,支路数b=n+l,其中,弱环配电网络的树包括树支和连支,树支条数与独立节点个数相等,连支条数与独立回路数相等,即树支n条,连支l条,支路数等于树支条数与连支条数的和。
6.如权利要求1所述的基于序分量法的弱环配电网三相解耦潮流计算方法,其特征在于,步骤(2)所述的弱环配电网三相解耦潮流计算方法包括以下步骤:
E、计算三序网络道路矩阵Ts(n×b阶)和回-支关联矩阵Bs(l×b阶),其中,下标s=0、1、2,分别表示三序网络模型中的零序、正序和负序网络模型;
J、计算k次迭代时的其中,为k次迭代时去除注入序电流为零的节点后形成的新注入序电流矩阵(m×1阶),m为节点注入序电流不为零的节点个数;为回路序阻抗矩阵的逆矩阵,即且s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型;
K、计算k次迭代时的其中,1n=[1,1,…,1]T,共n个1;s=0,1,2,分别表示三序网络模型中的零序、正序和负序网络模型;
7.如权利要求1所述的基于序分量法的弱环配电网三相解耦潮流计算方法,其特征在于,步骤J中,对于注入序电流为零的节点,在三序网络道路矩阵Ts中把该节点所对应行删去后形成新矩阵为Tsg(m×b阶),其中,下标s=0、1、2,分别表示三序网络模型中的零序、正序和负序网络模型,m为节点注入序电流不为零的节点个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210389056.1A CN102891485B (zh) | 2012-10-15 | 2012-10-15 | 基于序分量法的弱环配电网三相解耦潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210389056.1A CN102891485B (zh) | 2012-10-15 | 2012-10-15 | 基于序分量法的弱环配电网三相解耦潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102891485A true CN102891485A (zh) | 2013-01-23 |
CN102891485B CN102891485B (zh) | 2014-11-12 |
Family
ID=47534911
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210389056.1A Active CN102891485B (zh) | 2012-10-15 | 2012-10-15 | 基于序分量法的弱环配电网三相解耦潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102891485B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103679551A (zh) * | 2013-12-05 | 2014-03-26 | 江苏省电力公司 | 一种配电网潮流在线计算方法 |
CN103795057A (zh) * | 2014-01-23 | 2014-05-14 | 华北电力大学(保定) | 配电网潮流计算中基于搜索的配电网拓扑编号生成方法 |
CN105528530A (zh) * | 2016-02-22 | 2016-04-27 | 国家电网公司 | 配电网三相潮流计算分析方法 |
CN106159947A (zh) * | 2016-08-09 | 2016-11-23 | 河海大学 | 一种基于序分量的孤岛交直流混联微电网三相解耦潮流的计算方法 |
CN109980651A (zh) * | 2019-05-21 | 2019-07-05 | 华北水利水电大学 | 一种环网型配电***潮流计算方法 |
CN113964831A (zh) * | 2021-10-29 | 2022-01-21 | 国网北京市电力公司 | 负载调节方法、装置、电子设备及计算机可读存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102708250A (zh) * | 2012-05-10 | 2012-10-03 | 天津大学 | 基于树形分层双向迭代的电力***数字混合仿真方法 |
-
2012
- 2012-10-15 CN CN201210389056.1A patent/CN102891485B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102708250A (zh) * | 2012-05-10 | 2012-10-03 | 天津大学 | 基于树形分层双向迭代的电力***数字混合仿真方法 |
Non-Patent Citations (5)
Title |
---|
何禹清等: "基于序电流注入模型的三相潮流计算方法", 《中国电力》 * |
吴文传等: "配网潮流回路分析法", 《中国电机工程学报》 * |
李红伟等,: "含多台变压器支路的三相弱环配电网潮流计算", 《电力***保护与控制》 * |
李红伟等: "一种基于回路分析法的弱环配电网潮流直接算法", 《电力***保护与控制》 * |
王良缘等: "基于序分量的辐射网三相潮流求解", 《电力自动化设备》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103679551A (zh) * | 2013-12-05 | 2014-03-26 | 江苏省电力公司 | 一种配电网潮流在线计算方法 |
WO2015081679A1 (zh) * | 2013-12-05 | 2015-06-11 | 国家电网公司 | 一种配电网潮流在线计算方法 |
CN103795057A (zh) * | 2014-01-23 | 2014-05-14 | 华北电力大学(保定) | 配电网潮流计算中基于搜索的配电网拓扑编号生成方法 |
CN105528530A (zh) * | 2016-02-22 | 2016-04-27 | 国家电网公司 | 配电网三相潮流计算分析方法 |
CN106159947A (zh) * | 2016-08-09 | 2016-11-23 | 河海大学 | 一种基于序分量的孤岛交直流混联微电网三相解耦潮流的计算方法 |
CN106159947B (zh) * | 2016-08-09 | 2018-09-21 | 河海大学 | 一种基于序分量的孤岛交直流混联微电网三相解耦潮流的计算方法 |
CN109980651A (zh) * | 2019-05-21 | 2019-07-05 | 华北水利水电大学 | 一种环网型配电***潮流计算方法 |
CN109980651B (zh) * | 2019-05-21 | 2022-11-15 | 华北水利水电大学 | 一种环网型配电***潮流计算方法 |
CN113964831A (zh) * | 2021-10-29 | 2022-01-21 | 国网北京市电力公司 | 负载调节方法、装置、电子设备及计算机可读存储介质 |
CN113964831B (zh) * | 2021-10-29 | 2023-11-24 | 国网北京市电力公司 | 负载调节方法、装置、电子设备及计算机可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102891485B (zh) | 2014-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102891485A (zh) | 基于序分量法的弱环配电网三相解耦潮流计算方法 | |
CN103926538B (zh) | 基于aic准则的变阶数rc等效电路模型及实现方法 | |
CN102842907A (zh) | 基于道路矩阵的配电网三相解耦潮流计算方法 | |
CN103066593B (zh) | 含多类型分布式电源的弱环配电网三相潮流计算方法 | |
CN103632018B (zh) | 一种基于Simscape平台的电池建模方法 | |
CN103094900B (zh) | 考虑相序混合法的分布式发电配电网三相潮流计算方法 | |
CN103199525A (zh) | 一种基于等效节点注入电流变换的配电网潮流计算方法 | |
CN102609598B (zh) | 一种对大型电力***进行电磁暂态仿真的方法 | |
CN104181470A (zh) | 一种基于非线性预测扩展卡尔曼滤波的电池soc估计方法 | |
CN108054757B (zh) | 一种内嵌无功和电压的n-1闭环安全校核方法 | |
CN103810646B (zh) | 一种基于改进投影积分算法的有源配电***动态仿真方法 | |
CN102842908A (zh) | 含多变压器支路的配电网三相解耦潮流计算方法 | |
CN103023013A (zh) | 一种主动配电网三相潮流的计算方法 | |
CN104113061A (zh) | 一种含分布式电源的配电网三相潮流计算方法 | |
CN104009472B (zh) | 一种基于协同粒子群的配电网状态估计方法 | |
CN103955594A (zh) | 一种电力***动态等值方法 | |
CN106875055A (zh) | 基于ahp和ga的储能装置参与电网调频的参数优化方法 | |
CN104979840B (zh) | 一种主动配电网三相无功优化方法 | |
CN107846022A (zh) | 基于ilutp预处理并行迭代法的大规模配电网潮流分析方法 | |
CN105896558B (zh) | 一种基于vsc的upfc机电暂态模块化建模方法 | |
CN105896547A (zh) | 一种风电接入下的大电网分级电压控制方法 | |
Husain et al. | Load flow analysis of radial and mesh distribution system using ZIP model | |
CN105608256A (zh) | 一种吸纳不同步长接口延时的并行计算分网方法 | |
CN102157930B (zh) | 一种同塔双回直流输电线路直流侧谐波电流的计算方法 | |
CN104462815B (zh) | 一种基于潮流分布均衡度的目标网架安全性分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |