CN104578159B - 一种含分布式电源的配电网三相潮流修正方法 - Google Patents

一种含分布式电源的配电网三相潮流修正方法 Download PDF

Info

Publication number
CN104578159B
CN104578159B CN201510015450.2A CN201510015450A CN104578159B CN 104578159 B CN104578159 B CN 104578159B CN 201510015450 A CN201510015450 A CN 201510015450A CN 104578159 B CN104578159 B CN 104578159B
Authority
CN
China
Prior art keywords
mrow
msub
mtd
mtr
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
CN201510015450.2A
Other languages
English (en)
Other versions
CN104578159A (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.)
Electric Power Research Institute of State Grid Chongqing Electric Power Co Ltd
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd
Electric Power Research Institute of State Grid Tianjin Electric Power Co Ltd
Nari Technology Co Ltd
Original Assignee
Electric Power Research Institute of State Grid Chongqing Electric Power Co Ltd
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd
Electric Power Research Institute of State Grid Tianjin Electric Power Co Ltd
Nari Technology 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 Electric Power Research Institute of State Grid Chongqing Electric Power Co Ltd, State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd, Electric Power Research Institute of State Grid Tianjin Electric Power Co Ltd, Nari Technology Co Ltd filed Critical Electric Power Research Institute of State Grid Chongqing Electric Power Co Ltd
Priority to CN201510015450.2A priority Critical patent/CN104578159B/zh
Publication of CN104578159A publication Critical patent/CN104578159A/zh
Application granted granted Critical
Publication of CN104578159B publication Critical patent/CN104578159B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种含分布式电源的配电网三相潮流修正方法,具体包括以下几个步骤:构建分布式电源三相模型;建立含分布式电源的配电网三相系统模型;对配电网三相系统模型中的节点和支路进行编号;按节点的编号顺序存储配电网各个节点的三相负荷参数,按支路的编号顺序存储配电网各个支路的三相阻抗参数;对各个节点的三相负荷参数和各个支路的三相阻抗参数进行标幺化;分析含分布式电源配电网潮流计算存在的问题,对分布式电源并网电压进行修正;选择收敛精度,对含分布式电源的配电网三相系统进行三相前推回代潮流计算。本发明考虑了分布式电源模型输出电流和并网电压对称的特性,选择网侧电压的正序分量作为分布式电源的并网电压,参与潮流计算。

Description

一种含分布式电源的配电网三相潮流修正方法
技术领域
本发明涉及一种含分布式电源的配电网三相潮流修正方法,属于电力系统配电网潮流计算技术领域。
背景技术
能源和环境问题是当前世界各国关注的焦点。在我国,为实现国家“节能减排”的既定目标,电力行业在发电侧大力发展各种可再生能源发电,朝着电能生产的“低碳化”、“零排放”方向迈进,由此导致风力发电、光伏发电等分布式电源的快速发展。
随着用户对供电可靠性及电能质量要求的提高,作为电力系统相对薄弱的配电网日益得到重视,配网环网结构、辐射状运行,联络开关经常切换,运行方式多变,潮流计算作为电力系统分析最基本的计算,不仅可以计算网损、校验各种运行方式的合理性等,也可以为暂态计算提供初值,配电网基本潮流计算重要性不言而喻。
考虑到随着分布式电源在电力系统中的大量随机接入,分布式电源渗透率不断提高,随之产生的源侧三相不对称性对配电网产生了较大影响。而分布式电源模型本身的特性,也将对配电网潮流计算产生一定的影响。
从处理三相的方式上来看,电力系统潮流计算可分为相分量法、序分量法,但与输电网不同的是配电网的运行和网络参数不对称,其中网络参数的不对称使得序分量法解耦失效,配电网潮流计算必须基于相分量法,也即考虑全耦合的abc模型。
从传统潮流算法上来看,配电网潮流计算应用最多的是前推回代法和牛顿拉夫逊法:前推回代法充分利用了配电网辐射状的特点,具有收敛性好、占用内存少、计算速度快,是公认的高效率算法,但前推回代对于配网中存在的弱环处理能力弱,需要对配网各节点及支路合理编号;牛顿拉夫逊法理论上适合一切非线性方程组的求解、不受环网影响,但存在对初值敏感、收敛性差、数据存储量大、计算速度慢的缺点,而且牛顿法是利用注入功率代替了注入电流得到非线性方程组继而迭代求解,不好改变注入电流为正序的情况,不能保证实际计算过程中负荷分布式电源的模型特征。
发明内容
针对现有技术存在的不足,考虑到牛顿拉夫逊法采用注入功率代替注入电流,不能保证注入电流为正序,本发明目的是提供一种采用前推回代算法的含分布式电源的配电网三相潮流修正方法。
为了实现上述目的,本发明是通过如下的技术方案来实现:
本发明的一种含分布式电源的配电网三相潮流修正方法,具体包括以下几个步骤:
步骤S1,构建分布式电源三相模型;
步骤S2,建立含分布式电源的配电网三相系统模型;
步骤S3,根据分布式电源三相模型和配电网三相系统模型构成含分布式电源的配电网三相系统,对所述配电网三相系统模型中的节点和支路进行编号;
步骤S4,按节点的编号顺序存储配电网各个节点的三相负荷参数,按支路的编号顺序存储配电网各个支路的三相阻抗参数;
步骤S5,对各个节点的三相负荷参数和各个支路的三相阻抗参数进行标幺化;
步骤S6,分析含分布式电源配电网潮流计算存在的问题,对分布式电源并网电压进行修正;
步骤S7,选择收敛精度,对含分布式电源的配电网三相系统进行三相前推回代潮流计算。
步骤S1中,分布式电源三相模型的构建方法如下:
将分布式电源源侧视为直流电压源,分布式电源电路拓扑中主要由直流电压源、逆变器和滤波器构成;分布式电源并网采用PQ控制策略;由于逆变器的控制是在dq轴解耦进行的,采用正交park变换得到dq轴下电路模型,变换阵相角采用和网侧电压相位一致。
步骤S2中,含分布式电源的配电网三相系统模型的建立方法如下:
将配电网的上级变电站出口母线视为无穷大电源,作为潮流计算中平衡节点,等值电压幅值和相角为恒定已知量,并假定三相电压对称;
将配电网线路等效为集中参数模型,考虑各相自阻抗和相间互阻抗;
将支路上所有负荷等效为支路末节点集中负荷,负荷模型为不对称的恒功率PQ负荷;将分布式电源看作是对称的恒功率PQ负荷。
步骤S3中,对配电网三相系统模型中的节点和支路进行编号的方法如下:
将平衡节点的编号设为0,按照广度优先对其他节点进行编号;各支路编号取支路末节点的节点编号。
步骤S4中,如节点i,应当存储三相的有功功率Pia、Pib、Pic,无功功率Qia、Qib、Qic
如支路Li,应当存储三相的阻抗参数为其中,Zaa、Zbb、Zcc表示各相自阻抗,Zab、Zab、Zbc表示相间互阻抗,各相自阻抗含实部电阻和虚部电抗,两相间耦合阻抗只有虚部电抗。
步骤S5中,选取上级变电站出口母线电压作为基准电压,选取上级变电站容量为三相基准功率,并计算基准阻抗值,对各个节点的三相负荷参数和各个支路的三相阻抗参数进行标幺化。
步骤S6中,分析含分布式电源配电网潮流计算存在的问题,对分布式电源并网电压进行修正的方法如下:
由于分布式电源采用逆变器接口模型,控制器采用PQ控制方式,而PQ控制要求电网侧电压是对称的,所以需要对潮流算法进行修正;注入电流求解时,分布式电源节点先求取出各节点不对称电压的正序分量,功率对称、电压对称情形下求取得到的分布式电源节点注入电流为正序。
步骤S7中,三相前推回代潮流计算的方法如下:
选取收敛精度,依次经过计算各节点注入电流、用并网点正序电压修正分布式电源并网电压、回代计算各支路电流、前推计算各节点电压,循环迭代至各节点有功及无功功率偏差收敛。
注入电流的计算方法如下:
初始假设所有节点各相电压幅值为1,相角依次为0、-120°、120°,计算各节点注入电流,如节点i的三相注入电流为:
其中Sia、Sib、Sic为负荷三相功率;是节点i所在所有支路并联导纳及节点i并联电容器等效导纳组合得到的等效导纳矩阵;
回代计算支路电流的方法如下:
从末端节点到首端节点回代计算各支路电流,某支路电流等于该支路所接所有下一支路电流加上本支路末节点注入电流:
ILk表示k支路的支路电流,Ik表示k节点注入电流,ΣIl(k,x)表示所有k下一支路的支路电流求和,也即所有支路首节点编号为k、末节点编号任意的支路电流求和;
前推节点电压的计算方法如下:
从首端节点到末端节点前推计算各节点电压,某支路末节点电压为本支路首节点电压减去本支路压降,也即:
i、j分别表示支路的首、末节点编号,该支路编号为j记为支路Lj,阻抗阵表示支路Lj的串联阻抗矩阵,表示支路Lj从首节点流向末节点的电流;
迭代收敛方法如下:
利用前推回代潮流算法每步计算得到的节点注入电流和节点电压计算节点注入功率,与已知的负荷功率及并联导纳、电容器充电功率计算得到的注入功率作差,若各节点功率偏差小于收敛精度,迭代结束,否者将上一步中计算的电压重新带入注入电流,进行迭代计算;功率偏差如下:
式(5)右侧大括号中两项分别是电容器和并联导纳的注入功率、初始负荷功率。
本发明采用前推回代算法,在计算过程中对分布式电源数据进行修正;本发明考虑了分布式电源模型输出电流和并网电压对称的特性,选择网侧电压的正序分量作为分布式电源的并网电压,参与潮流计算。
附图说明
图1为本发明的一种含分布式电源的配电网三相潮流修正方法的工作流程图;
图2为本发明提供的一种含分布式电源的配电网三相潮流修正方法的实施例的工作流程图;
图3为本发明提供的一种含分布式电源的配电网三相潮流修正方法的配电网三相系统测试案例模型图。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施方式,进一步阐述本发明。
参见图1,本发明提供一种含分布式电源的配电网三相潮流修正方法,包括:
步骤S1,构建分布式电源三相模型:采用逆变器接口型分布式电源模型,控制策略采用PQ控制;
步骤S2,建立含分布式电源的配电网三相系统模型;
步骤S3,忽略配电网三相系统模型中的合环支路,构成配电网三相系统,对该配电网三相系统模型中的节点和支路进行编号;
步骤S4,按节点的编号顺序存储各个节点的三相负荷参数,按支路的编号顺序存储各个支路的三相阻抗参数;
步骤S5,对各个节点的三相负荷参数和各个支路的三相阻抗参数标幺化;
步骤S6,分析含分布式电源配电网潮流计算存在的问题,提出分布式电源并网电压修正方法;
步骤S7,选择收敛精度,对含分布式电源的配电网三相系统进行潮流计算。
步骤S1中构建分布式电源三相模型包括:
将所述分布式电源源侧视为直流电压源,分布式电源电路拓扑中主要由直流电压源、逆变器和滤波器构成;分布式电源并网采用PQ控制策略;由于逆变器的控制是在dq轴解耦进行的,可采用正交park变换得到dq轴下电路模型,变换阵相角采用和网侧电压相位一致。
步骤S2中建立含分布式电源的配电网三相系统模型包括:
将所述配电网的上级变电站出口母线视为无穷大电源,作为所述潮流计算中作为平衡节点,等值为电压幅值和相角为恒定已知量并假定三相电压对称;将所述支路上所有负荷等效为支路末节点集中负荷,模型为可考虑不对称的恒功率PQ负荷;将分布式电源看作是对称的恒功率PQ负荷;所述支路采用集中参数模型,考虑三相之间的耦合电抗。
步骤S3中对所述配电网三相系统模型中的节点和支路进行编号的方法包括:
对所述配电网三相系统,将所述平衡节点的编号设为0,按照广度优先对其他节点进行编号;对所述配电网三相系统,按照各支路末端节点编号对各支路进行编号。
步骤S4中存储节点i的三相负荷参数包括a、b和c三相的有功功率Pia、Pib、Pic和无功功率Qia、Qib、Qic
存储支路Li的三相阻抗参数为其中,各相阻抗含实部电阻和虚部电抗,两相间耦合阻抗只有虚部电抗。
步骤S5中选取所述上级变电站出口母线电压作为基准电压,选取所述上级变电站容量为三相基准功率,对所述各个节点的三相负荷参数和所述各个支路的三相阻抗参数进行标幺化。
步骤S6中分析含分布式电源的配电网潮流计算存在的问题,提出分布式电源并网电压修正方法包括:
由于分布式电源采用逆变器接口模型,控制器采用PQ控制方式,而PQ控制要求电网侧电压是对称的,所以需要对潮流算法进行修正;注入电流求解时,分布式电源节点先求取出各节点不对称电压的正序分量,功率对称、电压对称情形下求取得到的分布式电源节点注入电流为正序。
步骤S7中选择收敛精度,对含分布式电源的配电网三相系统进行潮流计算包括:
为消除负序和零序功率,从潮流中注入电流入手,只要注入电流为正序无论注入电压对称或不对称,注入功率只有正序注入功率;由于牛顿法是利用注入功率代替了注入电流得到非线性方程组继而迭代求解,不好改变注入电流为正序的情况,此处选择三相前推回代法;收敛判断时,分布式电源节点利用三相功率求和得到正序功率,与给定三相功率作差得到功率偏差,最后进行功率偏差收敛校验。
参见图2和图3,该配电网三相系统网络即为纯辐射网,看作一棵树,源节点编号为0,其他节点按树的广度进行编号,各支路编号取支路末节点的节点编号,包括L1~L32。
按节点编号顺序存储配电网各节点三相负荷参数,如节点i应当存储有功功率Pia、Pib、Pic和无功功率Qia、Qib、Qic;按支路编号顺序存储配电网各支路三相阻抗参数,如支路Li应当存储阻抗参数为各相自阻抗含实部电阻和虚部电抗,两相间耦合阻抗只有虚部电抗。
选取上级变电站出口母线电压作为基准电压,选取上级变电站容量为基准功率,对各个节点的三相负荷参数和各个支路的三相阻抗参数进行标幺化,本实施例中取UB=12.67kV,SB=10MW,对上一步中存储的各个节点的三相负荷参数和各个支路的三相阻抗参数进行标幺化,考虑节点的负荷功率是单相的,其基准值应当取阻抗不存在单相和三相概念,基准值取
传统三相潮流计算注入点的各相功率都是按照式(1)计算的,由于节点
电压的不对称,按此式计算分布式电源的注入电流即便是对称的三相功率(分布式电源给定的是三相对称功率),也会得出不对称的电流,分布式电源的注入功率就会包含负序和零序功率,与实际情况明显不符合。
为消除负序和零序功率,从潮流中注入电流入手,只要注入电流为正序无论注入电压对称或不对称,注入功率Sia+Sib+Sic只有正序注入功率。由于牛顿法是利用注入功率代替了注入电流得到非线性方程组继而迭代求解,不好改变注入电流为正序的情况,此处选择三相前推回代法。
三相系统前推回代计算:选取收敛精度,依次经过计算各节点注入电流、修正分布式电源并网电压、回代计算各支路电流、前推计算各节点电压再循环迭代至各节点有功、无功功率偏差收敛,本实施例中收敛精度选取为0.000001。
注入电流:初始假设所有节点各相电压幅值为1,相角依次为0、-120°、120°,计算各节点注入电流,如节点i的三相注入电流为:
其中Sia、Sib、Sic为负荷三相功率;是节点i所在所有支路并联导纳及节点i并联电容器等效导纳组合得到的等效导纳矩阵。
回代支路电流:从末端节点到首端节点回代计算各支路电流,某支路电流等于该支路所接所有下一支路电流加上本支路末节点注入电流:
ILk表示k支路的支路电流,ΣIl(k,x)表示所有k下一支路的支路电流求和,也即所有支路首节点编号为k、末节点编号任意的支路电流求和。
前推节点电压:从首端节点到末端节点前推计算各节点电压,某支路末节点电压为本支路首节点电压减去本支路压降,也即:
i、j分别表示支路的首、末节点编号,该支路编号为j记为支路Lj,阻抗阵表示支路Lj的串联阻抗矩阵。
迭代收敛判据:利用计算得到的节点注入电流和节点电压计算节点注入功率,与负荷功率及并联导纳、电容器充电功率计算得到的注入功率作差,若各节点功率偏差小于收敛精度,迭代结束,否者将上一步中计算的电压重新带入注入电流,进行迭代计算。具体功率偏差如下:
式(5)右侧大括号中两项分别是电容器和并联导纳的注入功率、初始负荷功率。
存储各节点电压幅值、相角。
案例采用IEEE33节点数据,在节点5、16、22三处并入分布式电源,分布式电源相关数据见下表1:
表1分布式电源相关参数
DG1 DG2 DG3
节点 5 16 22
功率(kW) 290+j20 310-j20 300+j20
结合案例仿真,利用Matlab编程处理,最终结果潮流迭代5次,误差0.000000075373425。
表2为利用三相前推回代法计算得到的各节点电压幅值和相角。
表2计算结果
节点号 a相幅值 b相幅值 c相幅值 a相相角 b相相角 c相相角
0 1 1 1 0 -120 120
1 0.99771 0.997795 0.99759 0.028067 -119.953 120.03
2 0.987206 0.987767 0.986507 0.184225 -119.7 120.1957
3 0.981646 0.982552 0.980571 0.290129 -119.528 120.3093
4 0.976304 0.977127 0.974736 0.398776 -119.347 120.4262
5 0.96251 0.962404 0.95965 0.485918 -119.082 120.516
6 0.959557 0.958767 0.956624 0.35214 -119.144 120.3649
7 0.956517 0.95612 0.95334 0.420122 -119.024 120.4274
8 0.952917 0.952485 0.94931 0.460111 -118.959 120.4454
9 0.949784 0.949264 0.945723 0.512582 -118.914 120.4725
10 0.949431 0.948828 0.945267 0.52972 -118.898 120.4886
11 0.948886 0.948144 0.944534 0.560447 -118.868 120.5169
12 0.946468 0.945166 0.941371 0.65306 -118.798 120.578
13 0.945556 0.944003 0.940168 0.68881 -118.777 120.5925
14 0.945628 0.943779 0.939919 0.73482 -118.743 120.6264
15 0.945788 0.943897 0.940018 0.787543 -118.692 120.678
16 0.946141 0.944351 0.940451 0.939051 -118.543 120.8302
17 0.945546 0.943766 0.939867 0.930529 -118.553 120.821
18 0.997189 0.997267 0.997055 0.01702 -119.964 120.0194
19 0.993747 0.993657 0.993325 -0.04855 -120.032 119.9536
20 0.993381 0.992916 0.99261 -0.07475 -120.052 119.935
21 0.992795 0.992204 0.991994 -0.09418 -120.075 119.9155
22 0.984555 0.985136 0.98378 0.181548 -119.693 120.1915
23 0.977959 0.978578 0.977023 0.088343 -119.765 120.093
24 0.974428 0.975398 0.973811 0.039556 -119.796 120.0438
25 0.960753 0.960255 0.958637 0.523776 -118.983 120.5436
26 0.958395 0.957606 0.956054 0.571769 -118.924 120.6013
27 0.947901 0.945796 0.944498 0.634417 -118.834 120.6908
28 0.940364 0.937346 0.93616 0.695829 -118.75 120.7732
29 0.937096 0.933706 0.932539 0.785662 -118.642 120.8819
30 0.933178 0.929491 0.928357 0.703844 -118.725 120.8002
31 0.932285 0.928556 0.92748 0.68144 -118.747 120.7781
32 0.93201 0.928268 0.927203 0.673737 -118.754 120.7705
鉴于三相前推回代算法主要应用于辐射状电网,对于环网的处理比较困难,本发明在含有环网较多的场景比较难以实现,但这并不影响对分布式电源并网时网侧电压的修正。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (3)

1.一种含分布式电源的配电网三相潮流修正方法,其特征在于,具体包括以下几个步骤:
步骤S1,构建分布式电源三相模型;
步骤S2,建立含分布式电源的配电网三相系统模型;
步骤S3,根据分布式电源三相模型和配电网三相系统模型构成含分布式电源的配电网三相系统,对所述配电网三相系统模型中的节点和支路进行编号;
步骤S4,按节点的编号顺序存储配电网各个节点的三相负荷参数,按支路的编号顺序存储配电网各个支路的三相阻抗参数;
步骤S5,对各个节点的三相负荷参数和各个支路的三相阻抗参数进行标幺化;
步骤S6,分析含分布式电源配电网潮流计算存在的问题,对分布式电源并网电压进行修正;
步骤S7,选择收敛精度,对含分布式电源的配电网三相系统进行三相前推回代潮流计算;步骤S1中,分布式电源三相模型的构建方法如下:
将分布式电源源侧视为直流电压源,分布式电源电路拓扑中主要由直流电压源、逆变器和滤波器构成;分布式电源并网采用PQ控制策略;由于逆变器的控制是在dq轴解耦进行的,采用正交park变换得到dq轴下电路模型,变换阵相角采用和网侧电压相位一致;步骤S2中,含分布式电源的配电网三相系统模型的建立方法如下:
将配电网的上级变电站出口母线视为无穷大电源,作为潮流计算中平衡节点,等值电压幅值和相角为恒定已知量,并假定三相电压对称;
将配电网线路等效为集中参数模型,考虑各相自阻抗和相间互阻抗;
将支路上所有负荷等效为支路末节点集中负荷,负荷模型为不对称的恒功率PQ负荷;将分布式电源看作是对称的恒功率PQ负荷;步骤S3中,对配电网三相系统模型中的节点和支路进行编号的方法如下:
将平衡节点的编号设为0,按照广度优先对其他节点进行编号;各支路编号取支路末节点的节点编号;步骤S4中,如节点i,应当存储三相的有功功率Pia、Pib、Pic,无功功率Qia、Qib、Qic
如支路Li,应当存储三相的阻抗参数为其中,Zaa、Zbb、Zcc表示各相自阻抗,Zab表示a、b两相间互阻抗,Zac表示a、c两相间互阻抗,Zbc表示b、c两相间互阻抗,各相自阻抗含实部电阻和虚部电抗,两相间耦合阻抗只有虚部电抗;步骤S5中,选取上级变电站出口母线电压作为基准电压,选取上级变电站容量为三相基准功率,并计算基准阻抗值,对各个节点的三相负荷参数和各个支路的三相阻抗参数进行标幺化;步骤S6中,分析含分布式电源配电网潮流计算存在的问题,对分布式电源并网电压进行修正的方法如下:
由于分布式电源采用逆变器接口模型,控制器采用PQ控制方式,而PQ控制要求电网侧电压是对称的,所以需要对潮流算法进行修正;注入电流求解时,分布式电源节点先求取出各节点不对称电压的正序分量,功率对称、电压对称情形下求取得到的分布式电源节点注入电流为正序。
2.根据权利要求1所述的含分布式电源的配电网三相潮流修正方法,其特征在于,
步骤S7中,三相前推回代潮流计算的方法如下:
选取收敛精度,依次经过计算各节点注入电流、用并网点正序电压修正分布式电源并网电压、回代计算各支路电流、前推计算各节点电压,循环迭代至各节点有功及无功功率偏差收敛。
3.根据权利要求2所述的含分布式电源的配电网三相潮流修正方法,其特征在于,
注入电流的计算方法如下:
初始假设所有节点各相电压幅值为1,相角依次为0、-120°、120°,计算各节点注入电流,如节点i的三相注入电流为:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>i</mi> <mi>a</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>i</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mo>-</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mrow> <mi>i</mi> <mi>a</mi> </mrow> </msub> <mo>/</mo> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>a</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> </msub> <mo>/</mo> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mrow> <mi>i</mi> <mi>c</mi> </mrow> </msub> <mo>/</mo> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>c</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>*</mo> </msup> <mo>-</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>a</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>a</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>a</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>a</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>b</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>b</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>a</mi> <mi>c</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>b</mi> <mi>c</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>c</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>a</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
其中Sia、Sib、Sic为负荷三相功率;是节点i所在所有支路并联导纳及节点i并联电容器等效导纳组合得到的等效导纳矩阵;
回代计算支路电流的方法如下:
从末端节点到首端节点回代计算各支路电流,某支路电流等于该支路所接所有下一支路电流加上本支路末节点注入电流:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>L</mi> <mi>k</mi> <mi>a</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>L</mi> <mi>k</mi> <mi>b</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>L</mi> <mi>k</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mo>-</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>k</mi> <mi>a</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>k</mi> <mi>b</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>k</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>l</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <mi>i</mi> <mo>)</mo> </mrow> </mrow> </munder> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>L</mi> <mi>i</mi> <mi>a</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>L</mi> <mi>i</mi> <mi>b</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <mi>L</mi> <mi>i</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
其中,ILia、ILib、ILic表示编号为i的支路流过的电流a、b、c三相的分量,ILka、ILkb、ILkc表示k支路的三相电流,Ika、Ikb、Ikc表示k节点的三相注入电流,l(k,i)表示所有首节点为k的支路;
前推节点电压的计算方法如下:
从首端节点到末端节点前推计算各节点电压,某支路末节点电压为本支路首节点电压减去本支路压降,也即:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>U</mi> <mrow> <mi>j</mi> <mi>a</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>U</mi> <mrow> <mi>j</mi> <mi>b</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>U</mi> <mrow> <mi>j</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>a</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>Z</mi> <mrow> <mi>a</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Z</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Z</mi> <mrow> <mi>a</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Z</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Z</mi> <mrow> <mi>b</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Z</mi> <mrow> <mi>b</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Z</mi> <mrow> <mi>a</mi> <mi>c</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Z</mi> <mrow> <mi>b</mi> <mi>c</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Z</mi> <mrow> <mi>c</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <msub> <mi>L</mi> <mi>j</mi> </msub> </msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <msub> <mi>L</mi> <mi>j</mi> </msub> <mi>a</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <msub> <mi>L</mi> <mi>j</mi> </msub> <mi>b</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mrow> <msub> <mi>L</mi> <mi>j</mi> </msub> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
i、j分别表示支路的首、末节点编号,该支路编号为j记为支路Lj,阻抗阵表示支路Lj的串联阻抗矩阵,表示支路Lj从首节点流向末节点的电流;
迭代收敛方法如下:
利用前推回代潮流算法每步计算得到的节点注入电流和节点电压计算节点注入功率,与已知的负荷功率及并联导纳、电容器充电功率计算得到的注入功率作差,若各节点功率偏差小于收敛精度,迭代结束,否者将上一步中计算的电压重新带入注入电流,进行迭代计算;功率偏差如下:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;S</mi> <mrow> <mi>i</mi> <mi>a</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;S</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;S</mi> <mrow> <mi>i</mi> <mi>c</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>a</mi> </mrow> </msub> <msubsup> <mi>I</mi> <mrow> <mi>i</mi> <mi>a</mi> </mrow> <mo>*</mo> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> </msub> <msubsup> <mi>I</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> <mo>*</mo> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>c</mi> </mrow> </msub> <msubsup> <mi>I</mi> <mrow> <mi>i</mi> <mi>c</mi> </mrow> <mo>*</mo> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>*</mo> </msup> <mo>-</mo> <mo>{</mo> <mo>-</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>a</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>a</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>a</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>a</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>b</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>b</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>a</mi> <mi>c</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>b</mi> <mi>c</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mi>c</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>*</mo> </msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>U</mi> <mrow> <mi>i</mi> <mi>a</mi> </mrow> <mn>2</mn> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>U</mi> <mrow> <mi>i</mi> <mi>b</mi> </mrow> <mn>2</mn> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>U</mi> <mrow> <mi>i</mi> <mi>c</mi> </mrow> <mn>2</mn> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;S</mi> <mrow> <mi>i</mi> <mi>a</mi> <mn>0</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;S</mi> <mrow> <mi>i</mi> <mi>b</mi> <mn>0</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;S</mi> <mrow> <mi>i</mi> <mi>c</mi> <mn>0</mn> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
式(5)右侧大括号中两项分别是电容器和并联导纳的注入功率、初始负荷功率。
CN201510015450.2A 2015-01-13 2015-01-13 一种含分布式电源的配电网三相潮流修正方法 Active CN104578159B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510015450.2A CN104578159B (zh) 2015-01-13 2015-01-13 一种含分布式电源的配电网三相潮流修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510015450.2A CN104578159B (zh) 2015-01-13 2015-01-13 一种含分布式电源的配电网三相潮流修正方法

Publications (2)

Publication Number Publication Date
CN104578159A CN104578159A (zh) 2015-04-29
CN104578159B true CN104578159B (zh) 2018-02-13

Family

ID=53093626

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510015450.2A Active CN104578159B (zh) 2015-01-13 2015-01-13 一种含分布式电源的配电网三相潮流修正方法

Country Status (1)

Country Link
CN (1) CN104578159B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105022885B (zh) * 2015-07-24 2018-12-18 天津大学 基于改进的Big Bang-Big Crunch的分布式光伏并网接纳能力计算方法
CN106815243B (zh) * 2015-11-29 2020-06-19 国网江西省电力公司赣西供电分公司 基于Web-GIS的分布式电源智能分析与决策系统
CN105591404A (zh) * 2016-03-01 2016-05-18 国家电网公司 双馈型风力发电机在电力系统潮流计算中的节点处理方法
CN107134779B (zh) * 2017-03-31 2022-02-01 中国电力科学研究院 一种配电网潮流算法的自动选择方法及装置
CN109478781B (zh) * 2017-05-15 2021-08-31 深圳大学 获取直流电力网功率传输系数的均衡电导补偿型对称方法
CN109257949B (zh) * 2017-05-15 2021-08-06 深圳大学 获取直流电力网潮流的等量电导补偿型全局线性偏心方法
CN109417295B (zh) * 2017-05-15 2021-08-10 深圳大学 获取直流电力网潮流的无损耗全局线性偏心方法
CN107257130B (zh) * 2017-07-25 2020-01-07 武汉大学 基于区域量测解耦的低压配电网损耗计算方法
CN108173272B (zh) * 2017-12-08 2021-04-06 云南电网有限责任公司 一种基于在线灵敏度的分布式电压控制方法
US11169187B2 (en) 2019-06-28 2021-11-09 King Fahd University Of Petroleum And Minerals Zig zag based load flow method and system for extended radial distribution systems
CN111308270B (zh) * 2020-03-06 2021-08-17 西南交通大学 一种基于节点故障注入电流的输电线路故障检测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103094900A (zh) * 2012-12-20 2013-05-08 河海大学 考虑相序混合法的分布式发电配电网三相潮流计算方法
CN103353918A (zh) * 2013-04-24 2013-10-16 国家电网公司 一种配电网三相牛顿拉夫逊法潮流计算初值选取方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5444131B2 (ja) * 2010-06-15 2014-03-19 株式会社日立製作所 配電系統潮流シミュレーション装置、配電系統潮流シミュレーション方法およびそのプログラム
CN103066593B (zh) * 2012-12-20 2014-08-13 河海大学 含多类型分布式电源的弱环配电网三相潮流计算方法
CN103248044A (zh) * 2013-05-07 2013-08-14 广东电网公司电力科学研究院 配电网潮流计算方法与系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103094900A (zh) * 2012-12-20 2013-05-08 河海大学 考虑相序混合法的分布式发电配电网三相潮流计算方法
CN103353918A (zh) * 2013-04-24 2013-10-16 国家电网公司 一种配电网三相牛顿拉夫逊法潮流计算初值选取方法

Also Published As

Publication number Publication date
CN104578159A (zh) 2015-04-29

Similar Documents

Publication Publication Date Title
CN104578159B (zh) 一种含分布式电源的配电网三相潮流修正方法
CN103353918B (zh) 一种配电网三相牛顿拉夫逊法潮流计算初值选取方法
Beerten et al. Generalized steady-state VSC MTDC model for sequential AC/DC power flow algorithms
CN102541621B (zh) 风光储联合发电系统的仿真等值方法
CN103094900B (zh) 考虑相序混合法的分布式发电配电网三相潮流计算方法
CN106446458B (zh) 考虑分布式电源的弱环配电网潮流计算方法
CN103971026B (zh) 主动配电网潮流计算通用方法
CN106295160B (zh) 交直流互联电网戴维南等值参数在线计算方法
CN105260949A (zh) 一种基于矩阵变换的配电网络短路计算方法
CN107272645B (zh) 中性点经电阻接地的光伏电站并网故障模型与分析方法
CN102879668A (zh) 含逆变型分布式电源配电网的非对称故障分析方法
CN105162099A (zh) 一种确定分布式发电接入电网不对称短路电流的运算曲面法
CN107576886A (zh) 含逆变型分布式电源小电阻接地系统的接地故障分析方法
CN106532710A (zh) 计及电压稳定约束的微电网潮流优化方法
CN107104463B (zh) 换流站在黑启动过程中考虑无功的有功投入量优化方法
CN107196307B (zh) 一种变压器故障后电网有功潮流快速估算的方法
CN104917197B (zh) 一种并行计算主动配电网三相不平衡潮流的方法
CN112052543B (zh) 一种基于混合整数二阶锥规划的保底网架搜索建模方法
CN106992522B (zh) 电压源换流器潮流计算模型及有源配电网潮流计算方法
CN103560494B (zh) 配电网保护整定的短路电流获取方法
CN107968398B (zh) 一种节点归并简化建模方法
CN107959292A (zh) 一种基于新型vsc潮流计算模型的交直流混合系统潮流计算方法
CN110165712B (zh) 一种基于网络流约束推导的骨干网架规划建模方法
CN104092213A (zh) 一种基于优化方法的不确定潮流支路功率分析方法
Yang et al. Power flow calculation for unbalanced three-phase distribution network with DGs based on phase-sequence hybrid modeling

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant