CN106446458A - 考虑分布式电源的弱环配电网潮流计算方法 - Google Patents
考虑分布式电源的弱环配电网潮流计算方法 Download PDFInfo
- Publication number
- CN106446458A CN106446458A CN201610942684.6A CN201610942684A CN106446458A CN 106446458 A CN106446458 A CN 106446458A CN 201610942684 A CN201610942684 A CN 201610942684A CN 106446458 A CN106446458 A CN 106446458A
- Authority
- CN
- China
- Prior art keywords
- node
- looped
- current
- matrix
- network
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 claims abstract description 47
- 230000008569 process Effects 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000010248 power generation Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 63
- 238000002347 injection Methods 0.000 claims description 43
- 239000007924 injection Substances 0.000 claims description 43
- 239000003990 capacitor Substances 0.000 claims description 4
- 238000000205 computational method Methods 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims description 4
- 239000000243 solution Substances 0.000 claims description 4
- 230000033228 biological regulation Effects 0.000 claims description 3
- 238000010937 topological data analysis Methods 0.000 claims description 3
- 239000000446 fuel Substances 0.000 abstract description 10
- 238000004458 analytical method Methods 0.000 description 10
- 239000007789 gas Substances 0.000 description 9
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000007363 ring formation reaction Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000004804 winding Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000567 combustion gas Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003032 molecular docking Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Water Supply & Treatment (AREA)
- Power Engineering (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑分布式电源的弱环配电网潮流计算方法,包括:通过利用叠加原理,将环网回路进行处理,得到辐射状网络以及端口补偿电路,解决了弱环配电网的潮流计算问题。给出了风力发电、燃料电池、光伏以及微型燃气轮机四种分布式电源在潮流计算中的处理方式。最终得出了含分布式电源的弱环配电网的潮流计算方法。本发明有益效果:该潮流算法简单清晰,使计算过程大幅简化,而且不需要对电网进行复杂的搜索。
Description
技术领域
本发明涉及弱环配电网潮流计算技术领域,尤其涉及考虑分布式电源的弱环配电网潮流计算方法。
背景技术
配电网是电力系统的重要组成部分,对配电网的分析必不可少。配电网潮流计算是配电网分析的一项基础内容。目前,已经提出了多种配电网潮流算法,有隐式Zbus高斯法、改进牛顿法、前推回代法和改进快速解耦法。而配电网有时会由于分段开关或联络开关闭合而形成弱环网,这些方法对弱环网的处理能力不强。而且随着分布式电源(DG)的大量加入,使得配电网的结构发生很大改变,需要对配电网潮流算法做进一步的研究。
现有技术中基于回路分析法,利用道路矩阵和回路矩阵,对各节点电压以及注入电流的关系矩阵进行了推导,具有较强的处理回路的能力,但是没考虑分布式电源的处理。
现有技术中提出的配电网潮流计算方法将弱环配电网等效为环网和辐射形网络两个部分,实现了含光伏电源的配电网潮流计算,然而对其它类型的分布式电源并未加以考虑。
现有技术中对多种分布式电源的运行和控制性能进行详细分析,采用将PV节点和环网打开的方式,提出了适用于含多种分布式电源的弱环配电网的前推回代三相潮流算法,对传统前推回代法本身不能处理环网的不足进行了改进。然而前推回代法需要对支路首尾节点进行繁琐的描述,将PV节点和环网打开又进一步增加了迭代的复杂性。
发明内容
本发明的目的就是为了解决上述问题,提出了一种考虑分布式电源的弱环配电网潮流计算方法,该方法相比于传统的前推回代法,该潮流算法简单清晰,使计算过程大幅简化,而且不需要对电网进行复杂的搜索;同时,该方法考虑了多种分布式电源在潮流计算中的处理。
为实现上述目的,本发明的具体方案如下:
一种考虑分布式电源的弱环配电网潮流计算方法,包括以下步骤:
(1)读入弱环配电网络参数,搜寻环网回路的连接点,然后将其断开,得到纯辐射状网络;在回路的解环点两侧各形成一个大小相等、极性相反的注入电流来模拟回路;
(2)进行拓扑分析,分别得到纯辐射状网络道路矩阵T、环网端口阻抗矩阵ZL、节点-端口关联矩阵ML以及PV节点的节点电抗矩阵X;
(3)初始化系统参数及各类分布式电源的参数,计算各类分布式电源的注入电流;
(4)计算各节点负荷的注入电流,在分布式电源的接入点处叠加分布式电源的注入电流;
(5)计算各支路电流和节点电压;
(6)计算环网回路端口的电压差,并对环网回路端口的注入电流进行修正;
(7)计算纯辐射状网络中各个PI节点的无功功率;
若计算得到的某个PI节点的无功功率发生越限,则将该PI节点处理成PQ节点,然后对该PI节点的注入电流进行修正;
(8)计算纯辐射状网络中各个PV节点的无功功率注入;
若计算得到的某个PV节点的无功功率注入量发生越限,则将该PV节点按照PQ节点处理,并对该PV节点的注入电流进行修正;
(9)判断系统中各节点相邻两次迭代的电压的幅值、PV节点的电压差以及环网回路端口电压差的变化量是否满足收敛规则,如果所有条件同时满足,转到下一步,否则跳转到步骤(4);
(10)对异步发电机的转差率进行计算,求解风力发电机所吸收的无功功率;
判断P恒定、Q=f(V)类型的分布式电源是否满足功率因数要求;若满足,潮流计算结束,并将结果输出;否则需调整接入的并联电容器的组数,然后跳转到步骤(4)。
进一步地,所述步骤(2)中,纯辐射状网络道路矩阵T的计算方法为:
规定道路的方向为电源点指向各个节点,道路的方向与各条支路的方向相同;如果支路j在道路i上,则T(i,j)=1,否则T(i,j)=0。
进一步地,所述步骤(4)中,各节点负荷的注入电流为:
其中,Ij为注入电流矩阵,在三相潮流计算中其阶数为解环后节点个数的3倍;ML为纯辐射状网络的节点-端口关联矩阵;IL为端口补偿电流,其计算方法为:
ZL为环网端口阻抗矩阵;为将环网解开得到的两个节点之间的电压差;k表示第k次潮流迭代计算。
进一步地,所述步骤(5)中,计算支路电流的方法为:
设支路电流为Ib,依据KCL定律,支路电流Ib与注入电流Ij有如下关系:
Ib=TTIj;
其中,T为纯辐射状网络的道路矩阵。
进一步地,所述步骤(5)中,计算节点电压的方法为:
设根节点电压为U0,各节点电压为U(n×1阶);
从任一节点开始沿着该节点所在的道路到达根节点所经过的支路电压之和,即为该节点与根节点之间的电压差:ΔU=U0E-U=TZ bIb;
各节点电压可以表示为:U=U0E-ΔU;
其中,E为n×1阶全1矩阵;Zb为n×n阶支路阻抗矩阵;Ib为支路电流。
进一步地,所述步骤(6)中,环网回路端口的电压差根据解环后两节点第k次潮流计算的电压值求得;
根据求得端口补偿电流IL,即环网回路端口的注入电流;ZL为环网端口阻抗矩阵。
进一步地,所述步骤(7)中,当计算得到的无功功率越限时,无功功率Qmax保持最大值且恒定;采用Qmax的值为输出功率,重新计算得到节点的注入电流;并将修正后的量代入下一次迭代过程。
本发明的有益效果:
本发明方法在前推回代法的基础上,提出了一种基于道路矩阵的潮流算法。对于弱环配电网,将其在合环点处打开,按照辐射型网络处理。然后根据端口补偿法,在回路解环点两侧各形成一个大小相等、极性相反的注入电流来模拟回路。由于道路矩阵为稀疏下三角阵,利用稀疏技术可以降低内存需求。相比于传统的前推回代法,该潮流算法简单清晰,使计算过程大幅简化,而且不需要对电网进行复杂的搜索。然后研究了多种分布式电源在潮流计算中的处理方法,得出了含分布式电源的弱环配电网潮流计算方法。利用IEEE 33节点模型,验证了所提出的算法的有效性。
附图说明
图1为简单网络示意图;
图2为多端口电路示意图;
图3为IEEE 33节点弱环网络图;
图4为不同接入方式下各节点的电压示意图。
具体实施方式:
下面结合附图对本发明进行详细说明:
本发明公开了一种考虑分布式电源的利用道路矩阵来计算弱环配电网潮流计算方法,通过利用叠加原理,将环网回路进行处理,得到辐射状网络以及端口补偿电路,解决了弱环配电网的潮流计算问题。给出了风力发电、燃料电池、光伏以及微型燃气轮机四种分布式电源在潮流计算中的处理方式。最终得出了含分布式电源的弱环配电网的潮流计算方法。通过算例分析对所提算法进行了验证,具有良好的收敛性。
下面对本发明的实现方法具体介绍如下:
1、基于道路矩阵的潮流算法
道路—支路关联矩阵,又称为道路矩阵。节点沿树到根所经过的路径上的所有支路集合称为道路。道路的特点为:节点的道路完全由树支支路组成;对于某一给定的树,节点的道路是唯一的。
道路可以用n×n阶的道路矩阵T来描述。规定道路的方向为电源点指向各个节点,道路的方向与各条支路的方向相同。如果支路j在道路i上,则T(i,j)=1,否则T(i,j)=0。
为了说明道路矩阵的形成过程,下面以图1所示的简单网络图来描述道路矩阵的建立过程。首先,以根节点为搜索起点,采用深度优先搜索算法进行搜索,并对各个节点和各条支路进行编号,如图1所示。Ij表示节点j负荷的注入电流(n×1阶)。该网络所建立的道路矩阵如式(1)所示。
式中,矩阵左侧纵向的数字表示的是各个节点的编号,上方数字则表示支路编号。可以看出,道路矩阵T是一个下三角阵。设支路电流为Ib(n×1阶),依据KCL定律,支路电流Ib与注入电流Ij的有如下关系:
Ib=TTIj (2)
对于辐射状配电网络,根据欧姆定律可以得到:
Ub=ZbIb (3)
式中,Zb为n×n阶支路阻抗矩阵。设根节点电压为U0,各节点电压为U(n×1阶)。可知从任一节点开始沿着该节点所在的道路到达电源节点所经过的支路电压之和,即为该节点与根节点之间的电压差:
ΔU=U0E-U=TZ bIb (4)
其中,E为n×1阶全1矩阵。所以,各节点电压可以表示为:
U=U0E-ΔU (5)
由于道路矩阵T为稀疏下三角阵,利用稀疏技术可以降低内存需求。而且,只需将矩阵中的元素更新为三相系统参数,就可实现三相潮流的计算。相比于传统的前推回代法,基于道路矩阵的潮流算法简单清晰,大大简化了计算过程。
另外,由于辐射状配电网节点互阻抗具有如下规律:节点i和节点k之间的互阻抗为两节点到根节点之间所经过的公共支路的阻抗和。因此,当电网结构为辐射状树支形网络时,其节点阻抗矩阵可以等效为将各条树支的阻抗逐条列写出来,利用道路矩阵可以得到,即:
Z=TZbTT (6)
T为道路矩阵,Zb为支路阻抗矩阵。这样就可以将节点阻抗矩阵直接求出,不需要对电网进行复杂的搜索,充分体现了道路矩阵的优越性。
2、环网回路的处理方法
在对有环网结构的配电网进行潮流计算时,首先需要在合环点将环网回路打开,得到辐射状配电网。然后根据端口补偿法,在回路的解环点两侧各形成一个大小相等、极性相反的注入电流来模拟回路。首先计算辐射状电网的潮流,然后对端口回路进行处理。该端口补偿电流IL可由下式计算得到:
其中为将环网解开得到的两个节点之间的电压差。该值可以由解环后两节点第k次潮流计算电压值求得。ZL为环网端口阻抗矩阵。当进行三相潮流计算时,其阶数为变为环网回路数的3倍。环网端口阻抗的参数可根据以下方法求得。
对于如图2所示的多端口电路,当从端口ij注入单位电流时,该端口的电压差即为端口ij的自阻抗,而从端口ij注入单位电流时,在端口pq产生的电压差即为端口pq的互阻抗,可以表示成:
其中,Zij,ij为端口ij的自阻抗,Zij,pq为端口ij和pq的互阻抗,Mij和Mpq为关联矢量;Zii,Zjj为节点i,j的自阻抗,Zij,Zip,Zjq,Ziq,Zjp为节点i,j,p,q之间的互阻抗。通过式(6)可知,节点阻抗矩阵Z可以由支路阻抗矩阵Zb和道路矩阵T计算得出,因此利用式(8)和(9)可以求解出多端口阻抗矩阵。
然后将所有解环端口的关联矢量按列排在一起,组成环网的节点-端口关联矩阵ML。
以图2所示电路为例,有:
ML=[M ij,M pq] (12)
则计及环网回路补偿后的系统注入电流为:
Ij为注入电流矩阵,在三相潮流计算中其阶数为解环后节点个数的3倍。
3、DG在潮流计算中的数学模型
3.1风力电源模型
风能作为一种清洁能源被广泛应用,其中双馈风力发电机组逐渐成为主流机型。双馈机一般采用恒功率因数运行方式,即在运行时保持定子侧的输出功率具有恒定的功率因数。根据这一特性,一般可将其处理成PQ节点。然而,现在一般要求双馈风力机组除了发出有功功率,还要对具有一定的维持系统电压水平的能力。这就要求在运行中根据系统的电压水平对双馈机所发出的无功功率进行调节。因而,处理双馈机的潮流计算模型时,需要分别考虑恒功率因数和恒电压两种控制方式。
首先,当双馈机采用恒功率因数控制方式时,通过对转子绕组外接电压的幅值和相角进行调节,以维持风电机组定子侧的功率因数保持不变,因此可以将风机视为PQ模型。弱功率因数为则有而由变流器发出或者吸收的无功则很小,因此将风机的无功功率近似等效为定子绕组的无功功率,即:
而在恒电压运行方式下,为了维持机端电压恒定,双馈发电机可以发出或者吸收无功功率。在风电机组无功调节的范围之内,风电场可以视为PV节点,其潮流计算模型为:
当所需无功越限时,则可看作无功功率为限值的PQ节点。
3.2燃料电池模型
燃料电池一般通过换流器与电力系统连接。换流器既可看成相角可调的电压源又能等效为电流源,可以对接口母线的电压和逆变器的输出功率进行单独调节。因此将燃料电池处理成PV结点。当逆变器无功输出达到上限时,可以将其转化为PQ结点。同时,当燃料电池的输出电流增加时,其终端电压也会下降。燃料电池输出的功率可以表示为:
式中,X=XT+Xs,XT和Xs分别为变压器和系统等效电抗,E为电池电动势,RFC为燃料电池内阻,IFC为燃料电池输出的直流,UFC为电池的直流输出电压,δFC为功率调节单元的超前角,m为换流器调制度,IT为变压器输入电流,UT为接口母线电压,Us为系统电压。
3.3光伏电池模型
光伏电池将太阳能直接转换成直流电能,通过逆变器并网。光伏发电系统逆变器的控制方式为外环电压控制、内环电流控制的双环控制结构。通过电压外环来稳定光伏阵列的直流输出电压,并给出电流内环的电流参考值信号。电流内环实现对并网电流的跟踪控制,并保证电流跟踪的快速性和准确性。
通常可以将光伏发电系统看作可维持输出电流恒定的PI节点,其潮流模型可描述为:
3.4微型燃气轮机模型
微型燃气轮机发电系统主要有两种结构:单轴结构和分轴结构。其中,单轴结构的微型燃气轮机,通过电压型控制逆变器并网的微型燃气轮机可以处理为PV节点,而通过电流型控制逆变器并网的微型燃气轮机可以处理成PI节点。分轴结构微型燃气轮的动力涡轮和燃气涡轮采用不同的转轴,动力涡轮通过同步发电机直接并网。一般要求同步发电机具有励磁调节能力。采用电压控制的微型燃气轮机在潮流计算中可处理成PV节点,采用功率因数控制的微型燃气轮机可处理成PQ节点。
4、含DG弱环配电网的潮流计算方法
含DG的弱环配电网潮流计算的具体步骤可以表示为:
(1)首先,读入弱环配电网络参数,搜寻环网回路的连接点,然后将其断开,得到纯辐射状网络;
(2)进行拓扑分析,得到道路矩阵T、环网端口阻抗矩阵ZL、节点-端口关联矩阵ML以及PV节点的节点电抗矩阵X;
(3)初始化系统参数及各类DG的参数,计算各类DG的注入电流;
(4)计算各节点负荷的注入电流,在DG的接入点处叠加DG的注入电流;
(5)计算各支路电流和节点电压;
(6)计算环网回路端口的电压差,并对环网回路端口的注入电流进行修正;
(7)计算PI节点的无功功率。若计算得到的无功功率发生越限,则将PI节点处理成PQ节点,然后将PI节点的注入电流进行修正。
当计算得到的无功功率越限时,无功功率Qmax保持最大值且恒定。电流的修正是指的,采用Qmax的值为输出功率,此时重新计算得到节点的注入电流。由于采用迭代计算的方法,将修正后的量代入下一次迭代即可。
(8)计算PV节点的无功功率注入。若计算得到的无功功率注入量发生越限,则将按照PQ节点处理,并对PV节点的注入电流进行修正。
(9)判断系统中各节点相邻两次迭代的电压的幅值、PV节点的电压差以及环网回路端口电压差的变化量是否满足收敛规则,如果所有条件同时满足,转到下一步,否则跳转到步骤(4);
收敛规则是指迭代计算停止迭代的条件,即两次计算所得结果的差值小于某一定值,则表示计算的精度符合要求,迭代结束。比如,对于电压量,判断相邻两次迭代的电压幅值差的最大值是否满足预先给定的精度要求,即:max|U(k+1)-U(k)|<ε,若满足,则停止计算,否则跳转到步骤(4)。
(10)对异步发电机的转差率进行计算,求解风力发电机所吸收的无功功率。判断P恒定、Q=f(V)类型的DG是否满足功率因数要求。若满足,潮流计算结束,并将结果输出。否则需调整接入的并联电容器的组数,然后跳转到步骤(4)。
5、算例分析
以图3所示IEEE 33节点网络为例,对风力发电、燃料电池、光伏以及微型燃气轮机四种分布式电源接入配网进行测试分析。该系统的三相参数见文献:一种少环配电网三相潮流计算新方法。变压器的变比为12.66/0.38,为Y0/Y0接线,归算到低压侧的阻抗为0.01+j0.04Ω。用于无功补偿的电容器QN-Unit=40Kvar,UN=12.66kV。
本发明对不含分布式电源以及两种分布式电源的接入方式进行潮流计算分析,结果如表1所示。
表1分布式电源多种接入方式下的分析测试结果
不同接入方式下各节点的电压如图4所示。为了简便,图中只展示了A相电压。
由表1可以看出,迭代次数并没有因为分布式电源的接入而增多,说明本文提出的潮流算法能够有效应用于含有多种DG的配电网潮流计算。另外,根据分析,当不含分布式电源时,采用传统的前推回代算法,迭代次数也为7次。这是因为本质上本文所提出的基于道路矩阵的方法仍是面向支路的前推回代法,因此本算法迭代次数与传统前推回代法无异,都具有良好的收敛性。然而由于采用了道路矩阵,省去了繁琐的支路首尾节点描述,计算程序清晰简单,大大简化了迭代过程,提高了算法计算的效率。
由图4可以看出,接入DG的节点处,电压升高明显。由此可以看出,DG的接入对于支撑电网电压具有一定作用。但是不同类型的DG对电压的支撑能力是不一样的,一般与DG的类型、运行和控制方式有关。因此,需要综合考虑多种因素,合理的选择DG的容量和接入点。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (7)
1.一种考虑分布式电源的弱环配电网潮流计算方法,其特征是,包括以下步骤:
(1)读入弱环配电网络参数,搜寻环网回路的连接点,然后将其断开,得到纯辐射状网络;在回路的解环点两侧各形成一个大小相等、极性相反的注入电流来模拟回路;
(2)进行拓扑分析,分别得到纯辐射状网络道路矩阵T、环网端口阻抗矩阵ZL、节点-端口关联矩阵ML以及PV节点的节点电抗矩阵X;
(3)初始化系统参数及各类分布式电源的参数,计算各类分布式电源的注入电流;
(4)计算各节点负荷的注入电流,在分布式电源的接入点处叠加分布式电源的注入电流;
(5)计算各支路电流和节点电压;
(6)计算环网回路端口的电压差,并对环网回路端口的注入电流进行修正;
(7)计算纯辐射状网络中各个PI节点的无功功率;
若计算得到的某个PI节点的无功功率发生越限,则将该PI节点处理成PQ节点,然后对该PI节点的注入电流进行修正;
(8)计算纯辐射状网络中各个PV节点的无功功率注入;
若计算得到的某个PV节点的无功功率注入量发生越限,则将该PV节点按照PQ节点处理,并对该PV节点的注入电流进行修正;
(9)判断系统中各节点相邻两次迭代的电压的幅值、PV节点的电压差以及环网回路端口电压差的变化量是否满足收敛规则,如果所有条件同时满足,转到下一步,否则跳转到步骤(4);
(10)对异步发电机的转差率进行计算,求解风力发电机所吸收的无功功率;
判断P恒定、Q=f(V)类型的分布式电源是否满足功率因数要求;若满足,潮流计算结束,并将结果输出;否则需调整接入的并联电容器的组数,然后跳转到步骤(4)。
2.如权利要求1所述的一种考虑分布式电源的弱环配电网潮流计算方法,其特征是,所述步骤(2)中,纯辐射状网络道路矩阵T的计算方法为:
规定道路的方向为电源点指向各个节点,道路的方向与各条支路的方向相同;如果支路j在道路i上,则T(i,j)=1,否则T(i,j)=0。
3.如权利要求1所述的一种考虑分布式电源的弱环配电网潮流计算方法,其特征是,所述步骤(4)中,各节点负荷的注入电流为:
其中,Ij为注入电流矩阵,在三相潮流计算中其阶数为解环后节点个数的3倍;ML为纯辐射状网络的节点-端口关联矩阵;IL为端口补偿电流,其计算方法为:
ZL为环网端口阻抗矩阵;为将环网解开得到的两个节点之间的电压差;k表示第k次潮流迭代计算。
4.如权利要求1所述的一种考虑分布式电源的弱环配电网潮流计算方法,其特征是,所述步骤(5)中,计算支路电流的方法为:
设支路电流为Ib,依据KCL定律,支路电流Ib与注入电流Ij有如下关系:
Ib=TTIj;
其中,T为纯辐射状网络的道路矩阵。
5.如权利要求1所述的一种考虑分布式电源的弱环配电网潮流计算方法,其特征是,所述步骤(5)中,计算节点电压的方法为:
设根节点电压为U0,各节点电压为U(n×1阶);
从任一节点开始沿着该节点所在的道路到达根节点所经过的支路电压之和,即为该节点与根节点之间的电压差:ΔU=U0E-U=TZbIb;
各节点电压可以表示为:U=U0E-ΔU;
其中,E为n×1阶全1矩阵;Zb为n×n阶支路阻抗矩阵;Ib为支路电流。
6.如权利要求1所述的一种考虑分布式电源的弱环配电网潮流计算方法,其特征是,所述步骤(6)中,环网回路端口的电压差根据解环后两节点第k次潮流计算的电压值求得;
根据求得端口补偿电流IL,即环网回路端口的注入电流;ZL为环网端口阻抗矩阵。
7.如权利要求1所述的一种考虑分布式电源的弱环配电网潮流计算方法,其特征是,所述步骤(7)中,当计算得到的无功功率越限时,无功功率Qmax保持最大值且恒定;采用Qmax的值为输出功率,重新计算得到节点的注入电流;并将修正后的量代入下一次迭代过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610942684.6A CN106446458B (zh) | 2016-11-01 | 2016-11-01 | 考虑分布式电源的弱环配电网潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610942684.6A CN106446458B (zh) | 2016-11-01 | 2016-11-01 | 考虑分布式电源的弱环配电网潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106446458A true CN106446458A (zh) | 2017-02-22 |
CN106446458B CN106446458B (zh) | 2017-11-10 |
Family
ID=58179024
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610942684.6A Expired - Fee Related CN106446458B (zh) | 2016-11-01 | 2016-11-01 | 考虑分布式电源的弱环配电网潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106446458B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107171329A (zh) * | 2017-07-11 | 2017-09-15 | 国网江苏省电力公司宿迁供电公司 | 一种含新能源并网的配电网潮流计算方法 |
CN108599167A (zh) * | 2018-01-15 | 2018-09-28 | 国网吉林省电力有限公司电力科学研究院 | 一种辐射状配电网的线性化潮流计算方法 |
CN109066696A (zh) * | 2018-03-26 | 2018-12-21 | 吉林大学 | 一种含有双馈感应式风电机组电网的稳态潮流计算方法 |
CN109361219A (zh) * | 2018-11-08 | 2019-02-19 | 国网新疆电力有限公司经济技术研究院 | 一种分布式电源配电网潮流计算方法及装置 |
CN109980651A (zh) * | 2019-05-21 | 2019-07-05 | 华北水利水电大学 | 一种环网型配电系统潮流计算方法 |
CN110445156A (zh) * | 2019-08-22 | 2019-11-12 | 浙江大学 | 含环网的三相不平衡的中低压主动配电网潮流计算方法 |
CN110676837A (zh) * | 2019-08-16 | 2020-01-10 | 广东电网有限责任公司 | 考虑分布式计算框架的配电系统电压功率灵敏度计算方法 |
CN110728032A (zh) * | 2019-09-20 | 2020-01-24 | 河海大学 | 一种计及环网的电-热互联综合能源系统快速潮流计算方法 |
CN111725812A (zh) * | 2020-06-04 | 2020-09-29 | 海仿(上海)科技有限公司 | 大规模配电系统的潮流计算方法、装置、设备和存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103208818A (zh) * | 2013-04-16 | 2013-07-17 | 山东科技大学 | 一种计及分布式电源的配网合环电流计算方法 |
CN105391086A (zh) * | 2015-11-30 | 2016-03-09 | 国家电网公司 | 基于分布式发电的环网系统潮流计算方法 |
-
2016
- 2016-11-01 CN CN201610942684.6A patent/CN106446458B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103208818A (zh) * | 2013-04-16 | 2013-07-17 | 山东科技大学 | 一种计及分布式电源的配网合环电流计算方法 |
CN105391086A (zh) * | 2015-11-30 | 2016-03-09 | 国家电网公司 | 基于分布式发电的环网系统潮流计算方法 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107171329A (zh) * | 2017-07-11 | 2017-09-15 | 国网江苏省电力公司宿迁供电公司 | 一种含新能源并网的配电网潮流计算方法 |
CN108599167B (zh) * | 2018-01-15 | 2020-10-20 | 国网吉林省电力有限公司电力科学研究院 | 一种辐射状配电网的线性化潮流计算方法 |
CN108599167A (zh) * | 2018-01-15 | 2018-09-28 | 国网吉林省电力有限公司电力科学研究院 | 一种辐射状配电网的线性化潮流计算方法 |
CN109066696A (zh) * | 2018-03-26 | 2018-12-21 | 吉林大学 | 一种含有双馈感应式风电机组电网的稳态潮流计算方法 |
CN109361219A (zh) * | 2018-11-08 | 2019-02-19 | 国网新疆电力有限公司经济技术研究院 | 一种分布式电源配电网潮流计算方法及装置 |
CN109980651A (zh) * | 2019-05-21 | 2019-07-05 | 华北水利水电大学 | 一种环网型配电系统潮流计算方法 |
CN109980651B (zh) * | 2019-05-21 | 2022-11-15 | 华北水利水电大学 | 一种环网型配电系统潮流计算方法 |
CN110676837A (zh) * | 2019-08-16 | 2020-01-10 | 广东电网有限责任公司 | 考虑分布式计算框架的配电系统电压功率灵敏度计算方法 |
CN110676837B (zh) * | 2019-08-16 | 2021-07-13 | 广东电网有限责任公司 | 考虑分布式计算框架的配电系统电压功率灵敏度计算方法 |
CN110445156A (zh) * | 2019-08-22 | 2019-11-12 | 浙江大学 | 含环网的三相不平衡的中低压主动配电网潮流计算方法 |
CN110728032A (zh) * | 2019-09-20 | 2020-01-24 | 河海大学 | 一种计及环网的电-热互联综合能源系统快速潮流计算方法 |
CN110728032B (zh) * | 2019-09-20 | 2021-04-27 | 河海大学 | 一种计及环网的电-热互联综合能源系统快速潮流计算方法 |
CN111725812A (zh) * | 2020-06-04 | 2020-09-29 | 海仿(上海)科技有限公司 | 大规模配电系统的潮流计算方法、装置、设备和存储介质 |
CN111725812B (zh) * | 2020-06-04 | 2023-04-11 | 海仿(上海)科技有限公司 | 大规模配电系统的潮流计算方法、装置、设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106446458B (zh) | 2017-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106446458B (zh) | 考虑分布式电源的弱环配电网潮流计算方法 | |
CN103094900B (zh) | 考虑相序混合法的分布式发电配电网三相潮流计算方法 | |
CN104331619B (zh) | 一种含逆变器接口电源电力系统的短路电流分析方法 | |
CN102148510A (zh) | 含分布式电源的配电网潮流算法 | |
CN104485684A (zh) | 一种含双馈风电机组的电力系统故障电流计算方法 | |
CN104578159A (zh) | 一种含分布式电源的配电网三相潮流修正方法 | |
CN109742763A (zh) | 一种考虑含多类型能源结构的配电网双向潮流计算方法 | |
CN103971026A (zh) | 主动配电网潮流计算通用方法 | |
CN106300427A (zh) | 一种统一潮流控制器在风电机组并网运行中应用的方法 | |
CN111900718A (zh) | 基于多级优化追赶变分迭代法的有源配电网动态仿真方法 | |
CN114511418A (zh) | 含逆变型分布式电源配电网短路电流计算方法 | |
CN109149583A (zh) | 有源配电网简洁前推回代潮流计算方法 | |
CN107171329A (zh) | 一种含新能源并网的配电网潮流计算方法 | |
Lei et al. | Modeling and analysis of MW-level grid-connected PV plant | |
CN111799812A (zh) | 一种基于无功电压灵敏度的海上风电场无功协调控制方法 | |
Yang et al. | Power flow calculation for unbalanced three-phase distribution network with DGs based on phase-sequence hybrid modeling | |
Xuefeng et al. | A brief analysis of distributed generation connected to distribution network | |
Kang et al. | Research on grid connection of wind farm based on VSC-HVDC | |
CN107994607A (zh) | 一种风电场无功补偿装置的容量优化配置方法及系统 | |
CN109149645B (zh) | 一种含有双馈感应式风电机组电网的暂态稳定计算方法 | |
Zhou et al. | Dynamic Model Validation and Harmonic Stability Analysis of Offshore Wind Power Plants | |
Liu et al. | Research on power flow algorithm for power system including wind farm | |
Xinggang et al. | Total transfer capability analysis in power system including large-scale wind farms | |
Abumeteir et al. | Impact of High Penetration Renewable Energy Systems on Low-Frequency Oscillations | |
Wang et al. | Analysis and Evaluation of High Voltage Distribution Network with Access to Distributed Energy Based on Power Flow Calculation |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171110 Termination date: 20181101 |