CN112952853A - 一种含分布式电源的配电网无功优化调度方法 - Google Patents
一种含分布式电源的配电网无功优化调度方法 Download PDFInfo
- Publication number
- CN112952853A CN112952853A CN202110119826.XA CN202110119826A CN112952853A CN 112952853 A CN112952853 A CN 112952853A CN 202110119826 A CN202110119826 A CN 202110119826A CN 112952853 A CN112952853 A CN 112952853A
- Authority
- CN
- China
- Prior art keywords
- power
- distribution network
- reactive
- power distribution
- capacitor bank
- 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.)
- Pending
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 185
- 238000005457 optimization Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000003990 capacitor Substances 0.000 claims abstract description 190
- 238000004364 calculation method Methods 0.000 claims abstract description 46
- 230000003068 static effect Effects 0.000 claims abstract description 39
- 230000009471 action Effects 0.000 claims abstract description 28
- 238000000638 solvent extraction Methods 0.000 claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims description 87
- 238000005192 partition Methods 0.000 claims description 31
- 238000013178 mathematical model Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 8
- KDYAPQVYJXUQNY-OPHDRXFHSA-N 1,2-di-(alpha-linolenoyl)-3-[alpha-D-galactosyl-(1->6)-beta-D-galactosyl]-sn-glycerol Chemical compound O[C@@H]1[C@H](O)[C@@H](O)[C@H](OC[C@@H](COC(=O)CCCCCCC\C=C/C\C=C/C\C=C/CC)OC(=O)CCCCCCC\C=C/C\C=C/C\C=C/CC)O[C@@H]1CO[C@@H]1[C@H](O)[C@@H](O)[C@@H](O)[C@@H](CO)O1 KDYAPQVYJXUQNY-OPHDRXFHSA-N 0.000 claims description 7
- 230000010354 integration Effects 0.000 claims description 6
- 238000010248 power generation Methods 0.000 claims description 6
- 150000001875 compounds Chemical class 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 description 10
- 238000002474 experimental method Methods 0.000 description 7
- 238000012360 testing method Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 230000005611 electricity Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000001276 controlling effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012271 agricultural production Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
- H02J3/1821—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/50—Controlling the sharing of the out-of-phase component
-
- 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/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种含分布式电源的配电网无功优化调度方法,属于电力技术领域,解决了现有技术无功功率分布不合理、电网中电压质量低、电网损耗高的问题。该方法包括如下步骤:在配电网中,设置一定数量的电容器组和分布式电源;根据配电网能量损耗与控制参数的关系,建立配电网无功静态优化模型,并进行潮流计算,获得各电容器组的投切容量增量,进而获得各电容器组的预动作时刻表;对配电网拓扑结构进行分区,根据上述预动作时刻表,以每一时段中配电网能量损耗均最小为目标,建立配电网无功动态优化模型,并进行潮流计算,确定各电容器组的日前调度方案以及分布式电源出力方案;执行所述日前调度方案、分布式电源出力方案。
Description
技术领域
本发明涉及电力技术领域,尤其涉及一种含分布式电源的配电网无功优化调度方法。
背景技术
在配电网中,长期存在无功功率分布不合理、无功补偿设备投入少、投入率低等现象,导致电网中电压质量低、电网损耗高等问题日益凸显。其中,电压质量直接影响工农业生产、电网运行,电压质量不能满足要求将导致产品质量下降,电网传输能力减弱,电网网络损耗增加。电网损耗过高会给电力部门带来极大的经济损失,导致电价上涨。电能质量、用电经济性问题得不到解决将严重影响广大用电客户的用电积极性。
分布式发电是解决上述问题的有效手段。分布式发电是指发电功率在几千瓦至数百兆瓦的小型模块化、分散式、布置在用户附近的高效、可靠的发电单元。但是,大规模分布式电源(DG)的接入改变了当前配电网的结构,其将由单电源、辐射型结构变为遍布电源和负荷的复杂网络,同时,配电网馈线中传输有功、无功的数量和方向也发生了改变,随之将引起配电网节点电压分布及网络损耗的变化,这种变化与分布式电源的渗透率和分布式电源的出力时序特性等有着密切的关系。传统的配电网规划、运行、保护等方面将面临重大变化。
发明内容
鉴于上述的分析,本发明实施例旨在提供一种含分布式电源的配电网无功优化调度方法,用以解决现有无功功率分布不合理、电网中电压质量低、电网损耗高的问题。
一方面,本发明实施例提供了一种含分布式电源的配电网无功优化调度方法,包括如下步骤:
在配电网中,设置预设数量的电容器组和分布式电源;
根据配电网能量损耗与控制参数的关系,建立配电网无功静态优化模型;
根据配电网能量损耗与控制参数的关系,建立配电网无功静态优化模型,并进行潮流计算,获得各电容器组的投切容量增量,进而获得各电容器组的预动作时刻表;
对配电网拓扑结构进行分区,根据上述预动作时刻表,以每一时段中配电网能量损耗均最小为目标,建立配电网无功动态优化模型,并进行潮流计算,确定各电容器组的日前调度方案以及分布式电源出力方案;
执行所述日前调度方案、分布式电源出力方案。
上述技术方案的有益效果如下:配电网无功优化是改善电压质量和降低电网损耗的有效手段。上述方法限定了在已知配电网结构参数以及负荷的情况下,对电容器组的投切容量、分布式电源出力的规划及运行控制,使得配电网系统能够在满足电压、潮流、功率等的约束条件下,达到电网损耗最小。
基于上述方法的进一步改进,所述根据配电网能量损耗与控制参数的关系,建立配电网无功静态优化模型的步骤,进一步包括:
选取影响配电网能量损耗的控制参数,建立所述控制参数与一天内配电网能量损耗的关系的数学模型;
以一天内配电网能量损耗最小为目标函数,以电容器投切次数限制以及分布式电源出力限制为约束条件,建立配电网无功静态优化模型。
上述进一步改进方案的有益效果是:对建立配电网无功静态优化模型的方法进行了限定,通过对一天内配电网能量损耗最小、电容器投切次数限制以及分布式电源出力限制,使得优化后获得的配电网能量损耗最小、电压质量高。
进一步,控制参数包括各每一时段开始时刻电容器组的投切容量Ci,t、平衡节点电压Usw,t、电压控制型DG的机端电压出力UDG,t、无功补偿型DG的无功出力QDG,t;其中,t表示t时刻,i表示第i个电容器;
所述控制参数与一天内配电网能量损耗的关系的数学模型为
式中,f为一天内配电网能量损耗,Tt为一天内的第t小时,Ploss,t为t时刻的配电网能量损耗函数。
上述进一步方案的有益效果是:对配电网的控制参数,以及控制参数与一天内配电网能量损耗的关系分别进行了限定。控制参数涉及投切容量Ci,t、平衡节点电压Usw,t、电压控制型DG的机端电压出力UDG,t、无功补偿型DG的无功出力,为获得最佳的日前调度方案、分布式电源出力方案奠定了基础。
进一步,所述配电网无功静态优化模型中,目标函数为
约束条件包括等式约束和不等式约束;所述等式约束为
式中,PGi,t、PDG,t、PLi,t分别为t时刻配电网中的节点i的发电机注入有功功率、分布式电源注入有功功率、负荷消耗有功功率;QGi,t、QDG,t、QCi,t、QLi,t分别为t时刻配电网中的节点i的发电机注入无功功率、分布式电源注入无功功率、电容器补偿容量、负荷消耗无功功率;Gij,t、Bij,t、δij,t分别为t时刻配电网中的节点i和j之间的电导、电纳和电压相角差;Vi,t为t时刻配电网中的节点i的节点电压;Ni为与节点i有关联的所有节点集合,包括节点i本身,j为Ni中节点j;
所述不等式约束为
Vi,min≤Vi,t≤Vi,max i∈Ni
Usw,min≤Usw,t≤Usw,max i=SW
0≤Ci,t≤Ci,max i∈Nc
UDG,min≤UDG,t≤UDG,max
QDG,min≤QDG,t≤QDG,max
式中,Vi,min和Vi,max分别为节点i的电压上下限;SW为配电网的平衡节点;Ci,max为节点i所接电容器组投切容量上限;Nc为所有电容器补偿节点集合;nc_max为一日内电容器组允许最大投切次数;UDG,min和UDG,max分别为电压控制型分布式电源输出电压上、下限;QDG,min和QDG,max分别为无功补偿型分布式电源的无功出力上、下限;Npv为电压控制型分布式电源集合;NPQ为无功补偿型分布式电源集合;Ck(t)和Ck(t-1)分别为第k个电容器组在t时刻和t-1时刻的投切容量;为异或运算符。
上述进一步方案的有益效果是:对配电网无功静态优化模型进行了限定。目标是使得一天内配电网的能量损耗最低,对无功补偿位置、无功补偿容量、发电机机端电压、无功补偿设备投切策略等进行了规划及运行控制,使配电网系统能够在满足电压、潮流、功率等约束条件下,达到电网损耗最小。
进一步,所述QDG,max、QDG,min通过下面公式获得
式中,SN为逆变器的额定视在功率;PDG为分布式电源的有功功率出力值。
上述进一步方案的有益效果是:对QDG,max、QDG,min的获得方式进行了限定,融入了对配电网中逆变器和分布式电源的有功功率出力值的考量。
进一步,通过下面步骤获得各电容器组的预动作时刻表:
获取配电网所有母线各节点当前时刻之前24h内每一小时内的有功功率、无功功率、负荷功率以及配电网各节点的节点电压,并通过下面公式分别建立有功功率矩阵PG、无功负荷矩阵PD、负荷功率矩阵PL
将所述有功功率矩阵PG、无功负荷矩阵PD、负荷功率矩阵PL中的数据带入配电网无功静态优化模型,获得静态优化时各电容器的投切容量矩阵C
通过下面公式获得第j个电容器组在t时刻的投切容量增量ΔCj,t
ΔCj,t=Cj,t+1-Cj,t
j=1 … Nc
将获得的ΔCj,t按照从大到小顺序排列,得到每个电容器组的投切容量增量表;
根据预设的一日内电容器组允许最大投切次数nc_max,识别上述投切容量增量表前nc_max个值对应的各时刻,作为该电容器组的日前调度预动作时刻,建立每个电容器组的预动作时刻表Tj
式中,“0”表示该电容器挡位为初始挡状态。
上述进一步方案的有益效果是:通过上述配电网无功静态优化模型的潮流解算方法依次获得每个电容器组的预动作时刻表Tj。这样处理得到的电容器组的预动作时刻表能够在每个电容器最需要动作的时刻进行投切,使每个电容器组根据自身情况选择动作时刻,解决了现有技术中所有电容器在同一时刻投切的问题。既满足了离散控制设备的动作次数约束,延长其使用寿命,增大经济效益;又很好地解决了动态无功优化的时间解耦问题,增强了配电网运行的稳定性。
进一步,所述对配电网拓扑结构进行分区的步骤,进一步包括:
根据每一分区有且只有一个电容器组或平衡节点为条件,对配电网电网模型进行分区。
上述进一步方案的有益效果是:获得Nc+NSW个分区。其中,Nc、NSW分别为配电网的电容器组总数和平衡节点个数。平衡节点一般只有一个。配电网结构确定后,平衡节点位置确定。
进一步,所述根据上述预动作时刻表,以每一时段中配电网能量损耗均最小为目标,建立配电网无功动态优化模型,并进行潮流计算,确定各电容器组的日前调度方案的步骤,细化为:
根据分区结果,将有功功率矩阵PG、无功负荷矩阵PD、负荷功率矩阵PL中同一分区内的功率数据分别排列在一起,获得PGc、PDc、PLc
根据每个电容器组的预动作时刻表,对上述PGc、PDc、PLc中的每个数据分别进行下面公式中的数据整合处理
ti,e≤t≤ti,e+1
将整合处理后的PGc、PDc、PLc的数据带入下面公式中的配电网无功动态优化模型一,并采用配电网无功静态优化模型中的约束条件约束,进行潮流计算,获得动态优化时各电容器组的投切容量矩阵C′
minft=min{Ploss,t(Usw,t,C1,t,C2,t…Ci,t…CNc,t)·Tt}
式中,ft为t时刻的配电网能量损耗,Usw,t为平衡节点在t时刻的电压,Ci,t为第i个电容器组在t时刻的投切容量;i=1 … Nc,Nc为电容器组总数;
根据各电容器组的投切容量矩阵C′,结合各电容器组的预动作时刻表,确定各电容器组的日前调度方案。
上述进一步方案的有益效果是:实现将静态计算用负荷数据按分区结果重新排列,将同一个分区内的所有节点负荷数据排列在一起。根据各电容器预动作时刻表对重排后的静态计算用负荷数据进行处理,得到动态计算用负荷数据C′,用于动态优化计算。实际上,处理后得到的动态计算用有功负荷数据C′的矩阵维数仍为24×N。
进一步,所述根据各电容器组的投切容量矩阵C′,结合各电容器组的预动作时刻表,确定各电容器组的日前调度方案的步骤,进一步包括:
若依据第i个电容器组的预动作时刻表,t时刻该电容器组该动作,则依据投切容量矩阵C′将其投切容量调至Ci,t′;
若依据第i个电容器组的预动作时刻表,t时刻该电容器组不动作,则其投切容量与上一时刻投切容量相同;
依次获得每个电容器组24h的投切方案,作为各电容器组考虑其动作次数约束的日前调度方案。
上述进一步方案的有益效果是:任意电容器组初始档位的选择基于时刻1的“动态计算用负荷数据”算出。数学模型中的变量类型由各电容器预动作时刻表动态确定,若依据第i个电容器的预动作时刻表,t时刻电容器组动作,则Ci,t处理为控制变量;若t时刻电容器组不动作,则Ci,t处理为常量,取为Ci,t=Ci,t-1。计算后可得到各电容器组与其预动作时刻表相对应的电容器预动作表,即为电容器考虑其动作次数约束的日前调度方案(即第二天的分布式电源启停状态和出力方案)。完成优化计算后,能保证调度方案中各电容器组优化结果为其在各个动作时段内最优,从而降低系统能量损耗;同时,通过配电网分区能够实现分区域控制,较整体优化计算速度更快,更容易收敛。
进一步,所述确定分布式电源出力方案的步骤,包括:
控制各电容器组执行所述日前调度方案;
再次获取配电网所有母线各节点当前时刻之前24h内每一小时内的有功功率、无功功率以及负荷功率数据,分别建立有功功率矩阵、无功负荷矩阵、负荷功率矩阵;
根据分区结果,对有功功率矩阵、无功负荷矩阵、负荷功率矩阵中同一分区内的功率数据分别排列在一起,对排列后的有功功率矩阵、无功负荷矩阵、负荷功率矩阵进行整合处理;
将整合处理后的有功功率矩阵、无功负荷矩阵、负荷功率矩阵的数据带入下面公式中的配电网无功动态优化模型二,并采用配电网无功静态优化模型中的约束条件约束,进行潮流计算,获得动态优化时分布式电源初值矩阵,作为分布式电源出力方案
上述进一步方案的有益效果是:DG出力不必考虑动作次数限制,所以优化计算时将其作为连续变量处理。计算后可得到DG出力值表,与电容器预动作表共同构成了含DG的配电网动态无功优化调度方案。由于电网运行中,电容器的投切将对后续时刻配电网状态产生影响,并且控制设备过多不利于集中控制,同时,为解决现有方法中各电容器动作量仅为其在当前时刻点的优化值的问题,上述方案对电网进行分区域控制,基于配电网分区提出一种新的负荷处理方法使优化计算后系统能量损耗在整个时段内最小。
本发明中,上述各技术方案之间还可以相互组合,以实现更多的优选组合方案。本发明的其他特征和优点将在随后的说明书中阐述,并且,部分优点可从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过说明书以及附图中所特别指出的内容中来实现和获得。
附图说明
附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制,在整个附图中,相同的参考符号表示相同的部件。
图1为本发明实施例1含分布式电源的配电网无功优化调度方法步骤示意图;
图2为本发明实施例2分区后配电网结构示意图;
图3为本发明实施例2中3组对比试验节点23所接电容器组动作情况。
具体实施方式
下面结合附图来具体描述本发明的优选实施例,其中,附图构成本申请一部分,并与本发明的实施例一起用于阐释本发明的原理,并非用于限定本发明的范围。
配电网由架空线路、杆塔、电缆、配电变压器、开关设备、无功补偿电容等配电设备及附属设施组成,它在电力网中的主要作用是分配电能。平衡节点是无功功率最低的节点。
实施例1
本发明的一个具体实施例,公开了一种含分布式电源的配电网无功优化调度方法,如图1所示,包括如下步骤:
S1.在配电网中,设置预设数量的电容器组和分布式电源。
其中,在配电网未设置电容器组及分布式电源之前,采用现有配电网无功优化模型(例如专利CN104600714B)进行潮流计算;电容器组设置在节点电压及无功功率较低点(节点电压或无功功率低于预设值的节点);分布式电源设置在有功功率相对较低点(有功功率低于预设值的节点)。
S2.根据配电网能量损耗与控制参数的关系,建立配电网无功静态优化模型,并进行潮流计算,获得各电容器组的投切容量增量,进而获得各电容器组的预动作时刻表。
S3.对配电网拓扑结构进行分区,根据上述预动作时刻表,以每一时段中配电网能量损耗均最小为目标,建立配电网无功动态优化模型,并进行潮流计算,确定各电容器组的日前调度方案以及分布式电源出力方案。
S4.执行所述日前调度方案、分布式电源出力方案。
与现有技术相比,本实施例提供的配电网无功优化方法是改善电压质量和降低电网损耗的有效手段。该方法限定了在已知配电网结构参数以及负荷的情况下,对电容器组的投切容量、分布式电源出力的规划及运行控制,使得配电网系统能够在满足电压、潮流、功率等的约束条件下,达到电网损耗最小。
实施例2
在实施例1的基础上进行优化,步骤S2进一步包括:
S21.选取影响配电网能量损耗的控制参数。所述控制参数包括各每一时段开始时刻电容器组的投切容量Ci,t、平衡节点电压Usw,t、电压控制型DG的机端电压出力UDG,t、无功补偿型DG的无功出力QDG,t。其中,t表示t时刻,i表示第i个电容器。
S22.建立所述控制参数与一天内配电网能量损耗的关系的数学模型
式中,f为一天内配电网能量损耗,Tt为一天内的第t小时,Ploss,t为t时刻的配电网能量损耗函数。
示例性地,Tt=1h,Ploss,1为1点的配电网能量损耗,用功率表示。
优选地,上述数学模型可采用
即
S23.以一天内配电网能量损耗最小为目标函数,以电容器投切次数限制以及分布式电源出力限制为约束条件,建立配电网无功静态优化模型。
具体地,步骤S23中的配电网无功静态优化模型中,目标函数为
约束条件包括等式约束和不等式约束。等式约束为
式中,PGi,t、PDG,t、PLi,t分别为t时刻配电网中的节点i的发电机注入有功功率、分布式电源注入有功功率、负荷消耗有功功率;QGi,t、QDG,t、QCi,t、QLi,t分别为t时刻配电网中的节点i的发电机注入无功功率、分布式电源注入无功功率、电容器补偿容量、负荷消耗无功功率;Gij,t、Bij,t、δij,t分别为t时刻配电网中的节点i和j之间的电导、电纳和电压相角差;Vi,t为t时刻配电网中的节点i的节点电压;Ni为与节点i有关联的所有节点集合,包括节点i本身,j为Ni中节点j。
不等式约束为
Vi,min≤Vi,t≤Vi,max i∈Ni
Usw,min≤Usw,t≤Usw,max i=SW
0≤Ci,t≤Ci,max i∈Nc
UDG,min≤UDG,t≤UDG,max QDG,min≤QDG,t≤QDG,max (6)
式中,Vi,min和Vi,max分别为节点i的电压上下限;SW为配电网的平衡节点;Ci,max为节点i所接电容器组投切容量上限;Nc为所有电容器补偿节点集合;nc_max为一日内电容器组允许最大投切次数;UDG,min和UDG,max分别为电压控制型分布式电源输出电压上、下限;QDG,min和QDG,max分别为无功补偿型分布式电源的无功出力上、下限;Npv为电压控制型分布式电源集合;NPQ为无功补偿型分布式电源集合;Ck(t)和Ck(t-1)分别为第k个电容器组在t时刻和t-1时刻的投切容量;为异或运算符。
由于用电负荷时刻在变化,约束条件增加了电容器投切次数限制以及DG出力限制。
优选地,所述QDG,max、QDG,min通过下面公式获得
式中,SN为逆变器的额定视在功率;PDG为分布式电源的有功功率出力值。
优选地,步骤S2中建立配电网无功静态优化模型后,通过下面步骤获得各电容器组的预动作时刻表:
S24.获取配电网所有母线各节点当前时刻之前24h内每一小时内的有功功率、无功功率、负荷功率以及配电网各节点的节点电压,并通过下面公式分别建立有功功率矩阵PG、无功负荷矩阵PD、负荷功率矩阵PL
按小时为单位将一天划分为24个时段,每个时段内的用电负荷数据(有功功率、无功功率、负荷功率)按照积分中值定理处理为一个常量,本领域技术人员能够理解。处理后的用电负荷数据将用于静态潮流计算,也称之为“静态计算用负荷数据”。
S25.将所述有功功率矩阵PG、无功负荷矩阵PD、负荷功率矩阵PL中的数据带入配电网无功静态优化模型,获得静态优化时各电容器的投切容量矩阵C
S26.通过下面公式获得第j个电容器组在t时刻的投切容量增量ΔCj,t
ΔCj,t=Cj,t+1-Cj,t
j=1 … Nc (10)
S27.将获得的ΔCj,t按照从大到小顺序排列,得到每个电容器组的投切容量增量表。
S28.根据预设的一日内电容器组允许最大投切次数nc_max,识别上述投切容量增量表前nc_max个值对应的各时刻,作为该电容器组的日前调度预动作时刻,建立每个电容器组的预动作时刻表Tj
式中,“0”表示该电容器挡位为初始挡状态,可以是任意任何时刻点,不影响动作次数。
优选地,步骤S3中的对配电网拓扑结构进行分区,包括:
S31.根据每一分区有且只有一个电容器组或平衡节点为条件,对配电网电网模型进行分区。
优选地,步骤S3中的根据上述预动作时刻表,以每一时段中配电网能量损耗均最小为目标,建立配电网无功动态优化模型,并进行潮流计算,确定各电容器组的日前调度方案,进一步细化为如下步骤:
S32.根据分区结果,将有功功率矩阵PG、无功负荷矩阵PD、负荷功率矩阵PL中同一分区内的功率数据分别排列在一起,获得PGc、PDc、PLc
S33.根据每个电容器组的预动作时刻表,对上述PGc、PDc、PLc中的每个数据分别进行下面公式中的数据整合处理
ti,e≤t≤ti,e+1 (13)
S34.将整合处理后的PGc、PDc、PLc的数据带入下面公式中的配电网无功动态优化模型一,并采用配电网无功静态优化模型中的约束条件约束,进行潮流计算,获得动态优化时各电容器组的投切容量矩阵C′
上述配电网无功动态优化模型一minf仍采用控制参数与一天内配电网能量损耗的关系的数学模型进行动态优化。相比配电网无功静态优化模型,区别只在于UDG,t、QDG,t为0。即
Ploss,t(Ci,t,Usw,t)=Ploss,t(Ci,t,Usw,t,0,0) (16)
S35.根据各电容器组的投切容量矩阵C′,结合各电容器组的预动作时刻表,确定各电容器组的日前调度方案。
优选地,步骤S35进一步细化为:
S351.执行如下判断:若依据第i个电容器组的预动作时刻表,t时刻该电容器组该动作,则依据投切容量矩阵C′将其投切容量调至Ci,t′;若依据第i个电容器组的预动作时刻表,t时刻该电容器组不动作,则其投切容量与上一时刻投切容量相同,不再调至C1,t′;i=1… Nc。
S352.依次获得每个电容器组24h的投切方案,作为各电容器组考虑其动作次数约束的日前调度方案。
任意电容器组初始档位的选择基于时刻1的“动态计算用负荷数据”算出。数学模型中的变量类型由各电容器预动作时刻表动态确定,若依据第i个电容器的预动作时刻表,t时刻电容器组动作,则Ci,t处理为控制变量;若t时刻电容器组不动作,则Ci,t处理为常量,取为Ci,t=Ci,t-1。计算后可得到各电容器组与其预动作时刻表相对应的电容器预动作表,即为电容器考虑其动作次数约束的日前调度方案(即第二天的分布式电源启停状态和出力方案)。
当电容器调度方案确定后,接下来对电网中的分布式电源出力情况进行优化处理。根据分布式电源出力情况的控制特点,将分布式电源出力情况分为电压控制型(P已知、V可调)、无功补偿型(P已知、Q可调)、无功负荷型(P和Q均已知)3类,其中,P为分布式电源的有功功率,Q为分布式电源的无功功率,V为分布式电源的机端电压。
优选地,步骤S3中的确定分布式电源出力方案,进一步包括:
S36.控制各电容器组严格执行所述日前调度方案。即所有电容器动作量均按照预动作表中对应时刻的对应量值输入配电网系统。
S37.再次获取配电网所有母线各节点当前时刻之前24h内每一小时内的有功功率、无功功率以及负荷功率数据,分别建立有功功率矩阵、无功负荷矩阵、负荷功率矩阵。具体地,采用步骤S24中的建立方法。
S38.根据分区结果,对有功功率矩阵、无功负荷矩阵、负荷功率矩阵中同一分区内的功率数据分别排列在一起,对排列后的有功功率矩阵、无功负荷矩阵、负荷功率矩阵进行整合处理。具体地,采用步骤S31~S33的分区整合方法。
S39.将整合处理后的有功功率矩阵、无功负荷矩阵、负荷功率矩阵的数据带入下面公式中的配电网无功动态优化模型二,并采用配电网无功静态优化模型中的约束条件约束,进行潮流计算,获得动态优化时分布式电源初值矩阵,作为分布式电源出力方案
上述配电网无功动态优化模型二实质上仍是采用控制参数与一天内配电网能量损耗的关系的数学模型进行动态优化。相比配电网无功静态优化模型,区别只在于Ci,t已经通过日前调度方案确定,Usw,t在执行日前调度方案时候后通过测试可确定,只需求解UDG,t、QDG,t。
由于分布式电源出力不必考虑动作次数限制,所以优化计算时将其作为连续变量处理。计算后可得到分布式电源出力值表,与电容器预动作表共同构成了含分布式电源的配电网动态无功优化调度方案。
根据步骤S31~S39的方法进行配电网分区以及负荷数据处理,完成优化计算后能保证调度方案中各电容器组优化结果为其在各个动作时段内最优,从而降低系统能量损耗;同时,通过配电网分区能够实现分区域控制,较整体优化计算速度更快,更容易收敛。
为了更好地说明本实施例发明的效果,下面提供一个具体实例。节点配电网可采用IEEE 33节点配电网,在此基础上进行改造。改造后,该节点系统共38个节点,其中包含平衡节点1个、电容器补偿节点3个、分布式电源(DG)接入节点5个,如图2所示。系统额定电压设置为10.5kV,基准容量为100MVA,PV节点和平衡节点电压上下限设定为0.9~1.1(标准值),PQ节点电压上下限设定为0.95~1.05(标准值)。其中,电容器组安装地点为节点13、节点23、节点29,其补偿容量上限值由其所在支路配变容量总和的20%确定;PV可控型光伏发电通过电压型逆变器并网,PQ可控型光伏发电通过电流型逆变器并网。
将分布式电源、逆变器及升压变压器整体作为控制设备;分布式电源的渗透率控制在20%到25%之间,PQ可控型分布式电源无功出力上下限由下式确定,PV可控型DG控制变量为其机端电压。表1为该系统接入DG的具体类型、位置以及容量。
式中,SN为逆变器的额定视在功率;PDG为分布式电源的有功功率出力值。
表1算例中DG接入位置与参数
利用配电网分区方法将改造后的配电网分为了4个区,每个分区内包含一个电容器补偿点或者平衡节点。分区后配电网结构如图2所示。
在优化计算求解约束问题时采用MATLAB优化工具箱中的fmincon求解器进行求解。该求解器对初始值的选择无任何要求,收敛速度快,使用简便,结果准确。
表2列出了5组试验得到的优化结果。其中,试验1为优化前(电容器、DG均不参与调度,出力均为0);试验2为DG单独参与调度(电容器不动作,出力为0);试验3为按现有方法优化计算(未接入DG);试验4为按实施例2方法优化电容器(未接入DG);试验5为按实施例2调度方案优化计算(电容器、DG均参与调度)。
表2 5组试验结果对比
可见,试验2表明DG按本实施例方法参与无功调度降损效果明显,能充分发挥其无功补偿能力。试验4较试验3配电网系统能量损耗降低1.08%,表明按本实施例电容器预动作表得到的调度方案能保证配电网系统运行的经济性和稳定性,较现有方法更优。试验5表明本实施例含DG的动态无功优化调度方案能大幅降低系统损耗(较试验1降低56.03%,较现有方法降低1.370%),对DG和电容器组的协调控制保证了该调度方案的合理性和有效性。
图3为静态优化与本实施例动态优化计算得出的算例中节点23所接电容器组的投切容量对比。可以看出,利用本实施例方法计算获得的结果相比现有方法更为契合静态计算结果。由于静态结果为不考虑电容器组动作次数约束时每个时刻的最优值,那么在满足约束的情况下,优化结果越接近静态结果,说明该方法的优化效果越好。
结合算例结果,得出本发明方法相比现有技术具有如下优点:
1)DG参与配电网无功调度能大大减小系统的能量损耗。这是因为DG具有较强的电压无功调节能力,且出力是可连续调节的。因此,充分发挥DG的无功调节能力,有利于提高系统的安全性和经济性。
2)提出的电容器预动作表能严格控制电容器动作次数,保证其分散投切,计算后得到的电容器控制策略为各个动作时段内最优。实例结果表明本发明方法能有效地降低系统的损耗,提高系统运行的经济性和稳定性。
3)根据DG和电容器的出力特点,所提出的DG、电容器协调调度的控制方式能充分发挥二者的无功调节能力,算例结果验证了含DG的动态无功优化调度方法的合理性和有效性。
4)由于分布式电源的出力是连续可调的,所以电容器在系统无功调节中承担粗调任务,分布式电源则发挥细调作用。以每个时段内系统能量损耗最小为目标,以电容器投切容量为控制变量进行优化计算,得到各电容器的日前调度方案。
本领域技术人员可以理解,实现上述实施例方法的全部或部分流程,可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于计算机可读存储介质中。其中,所述计算机可读存储介质为磁盘、光盘、只读存储记忆体或随机存储记忆体等。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种含分布式电源的配电网无功优化调度方法,其特征在于,包括如下步骤:
在配电网中,设置一定数量的电容器组和分布式电源;
根据配电网能量损耗与控制参数的关系,建立配电网无功静态优化模型,并进行潮流计算,获得各电容器组的投切容量增量,进而获得各电容器组的预动作时刻表;
对配电网拓扑结构进行分区,根据上述预动作时刻表,以每一时段中配电网能量损耗均最小为目标,建立配电网无功动态优化模型,并进行潮流计算,确定各电容器组的日前调度方案以及分布式电源出力方案;
执行所述日前调度方案、分布式电源出力方案。
2.根据权利要求1所述的配电网无功优化调度方法,其特征在于,所述根据配电网能量损耗与控制参数的关系,建立配电网无功静态优化模型,进一步包括:
选取影响配电网能量损耗的控制参数,建立所述控制参数与一天内配电网能量损耗的关系的数学模型;
以一天内配电网能量损耗最小为目标函数,以电容器投切次数限制以及分布式电源出力限制为约束条件,建立配电网无功静态优化模型。
4.根据权利要求3所述的配电网无功优化调度方法,其特征在于,所述配电网无功静态优化模型中,目标函数为
约束条件包括等式约束和不等式约束;所述等式约束为
式中,PGi,t、PDG,t、PLi,t分别为t时刻配电网中的节点i的发电机注入有功功率、分布式电源注入有功功率、负荷消耗有功功率;QGi,t、QDG,t、QCi,t、QLi,t分别为t时刻配电网中的节点i的发电机注入无功功率、分布式电源注入无功功率、电容器补偿容量、负荷消耗无功功率;Gij,t、Bij,t、δij,t分别为t时刻配电网中的节点i和j之间的电导、电纳和电压相角差;Vi,t为t时刻配电网中的节点i的节点电压;Ni为与节点i有关联的所有节点集合,包括节点i本身;j为Ni中的节点j;
所述不等式约束为
Vi,min≤Vi,t≤Vi,max i∈Ni
Usw,min≤Usw,t≤Usw,max i=SW
0≤Ci,t≤Ci,max i∈Nc
UDG,min≤UDG,t≤UDG,max
QDG,min≤QDG,t≤QDG,max
6.根据权利要求5所述的配电网无功优化调度方法,其特征在于,通过下面步骤获得各电容器组的预动作时刻表:
获取配电网所有母线各节点当前时刻之前24h内每一小时内的有功功率、无功功率、负荷功率以及配电网各节点的节点电压,并通过下面公式分别建立有功功率矩阵PG、无功负荷矩阵PD、负荷功率矩阵PL
式中,N是配电网节点总数;
将所述有功功率矩阵PG、无功负荷矩阵PD、负荷功率矩阵PL中的数据带入配电网无功静态优化模型,获得静态优化时各电容器的投切容量矩阵C
通过下面公式获得第j个电容器组在t时刻的投切容量增量ΔCj,t
ΔCj,t=Cj,t+1-Cj,t
j=1 … Nc
将获得的ΔCj,t按照从大到小顺序排列,得到每个电容器组的投切容量增量表;
根据预设的一日内电容器组允许最大投切次数nc_max,识别上述投切容量增量表前nc_max个值对应的各时刻,作为该电容器组的日前调度预动作时刻,建立每个电容器组的预动作时刻表Tj
式中,0表示该电容器挡位为初始挡状态。
7.根据权利要求6所述的配电网无功优化调度方法,其特征在于,所述对配电网拓扑结构进行分区的步骤进一步包括:
根据每一分区有且只有一个电容器组或平衡节点为条件,对配电网电网模型进行分区。
8.根据权利要求7所述的配电网无功优化调度方法,其特征在于,所述根据上述预动作时刻表,以每一时段中配电网能量损耗均最小为目标,建立配电网无功动态优化模型,并进行潮流计算,确定各电容器组的日前调度方案的步骤进一步包括:
根据分区结果,将有功功率矩阵PG、无功负荷矩阵PD、负荷功率矩阵PL中同一分区内的功率数据分别排列在一起,获得PGc、PDc、PLc
根据每个电容器组的预动作时刻表,对上述PGc、PDc、PLc中的每个数据分别进行下面公式中的数据整合处理
将整合处理后的PGc、PDc、PLc的数据带入下面公式中的配电网无功动态优化模型一,并采用配电网无功静态优化模型中的约束条件约束,进行潮流计算,获得动态优化时各电容器组的投切容量矩阵C′
根据各电容器组的投切容量矩阵C′,结合各电容器组的预动作时刻表,确定各电容器组的日前调度方案。
9.根据权利要求8所述的配电网无功优化调度方法,其特征在于,所述根据各电容器组的投切容量矩阵C′,结合各电容器组的预动作时刻表,确定各电容器组的日前调度方案的步骤,进一步包括:
若依据第i个电容器组的预动作时刻表,t时刻该电容器组该动作,则依据投切容量矩阵C′将其投切容量调至Ci,t′;
若依据第i个电容器组的预动作时刻表,t时刻该电容器组不动作,则其投切容量与上一时刻投切容量相同;
依次获得每个电容器组24h的投切方案,作为各电容器组考虑其动作次数约束的日前调度方案。
10.根据权利要求9所述的配电网无功优化调度方法,其特征在于,所述确定分布式电源出力方案,进一步包括:
控制各电容器组执行所述日前调度方案;
再次获取配电网所有母线各节点当前时刻之前24h内每一小时内的有功功率、无功功率以及负荷功率数据,分别建立有功功率矩阵、无功负荷矩阵、负荷功率矩阵;
根据分区结果,对有功功率矩阵、无功负荷矩阵、负荷功率矩阵中同一分区内的功率数据分别排列在一起,对排列后的有功功率矩阵、无功负荷矩阵、负荷功率矩阵进行整合处理;
将整合处理后的有功功率矩阵、无功负荷矩阵、负荷功率矩阵的数据带入下面公式中的配电网无功动态优化模型二,并采用配电网无功静态优化模型中的约束条件约束,进行潮流计算,获得动态优化时分布式电源初值矩阵,作为分布式电源出力方案
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110119826.XA CN112952853A (zh) | 2021-01-28 | 2021-01-28 | 一种含分布式电源的配电网无功优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110119826.XA CN112952853A (zh) | 2021-01-28 | 2021-01-28 | 一种含分布式电源的配电网无功优化调度方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112952853A true CN112952853A (zh) | 2021-06-11 |
Family
ID=76238803
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110119826.XA Pending CN112952853A (zh) | 2021-01-28 | 2021-01-28 | 一种含分布式电源的配电网无功优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112952853A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114817865A (zh) * | 2022-06-28 | 2022-07-29 | 深圳市傲立电子有限公司 | 一种基于大数据的射频功率放大器用运行故障预测系统 |
CN114899825A (zh) * | 2022-04-20 | 2022-08-12 | 国家电网有限公司技术学院分公司 | 一种基于网络等值变换的配电网综合优化方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102882218A (zh) * | 2012-10-22 | 2013-01-16 | 上海交通大学 | 配电网无功/电压的静态集中与动态分布协调控制系统 |
CN104868481A (zh) * | 2015-05-06 | 2015-08-26 | 国家电网公司 | 一种使配电网络全天有功损耗最小的电容补偿方法 |
CN106849097A (zh) * | 2017-04-13 | 2017-06-13 | 云南电网有限责任公司电力科学研究院 | 一种主动配电网潮流计算方法 |
-
2021
- 2021-01-28 CN CN202110119826.XA patent/CN112952853A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102882218A (zh) * | 2012-10-22 | 2013-01-16 | 上海交通大学 | 配电网无功/电压的静态集中与动态分布协调控制系统 |
CN104868481A (zh) * | 2015-05-06 | 2015-08-26 | 国家电网公司 | 一种使配电网络全天有功损耗最小的电容补偿方法 |
CN106849097A (zh) * | 2017-04-13 | 2017-06-13 | 云南电网有限责任公司电力科学研究院 | 一种主动配电网潮流计算方法 |
Non-Patent Citations (1)
Title |
---|
刘公博等: "含分布式电源的配电网动态无功优化调度方法", 《电力系统自动化》, vol. 39, no. 15, 10 August 2015 (2015-08-10), pages 49 - 54 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114899825A (zh) * | 2022-04-20 | 2022-08-12 | 国家电网有限公司技术学院分公司 | 一种基于网络等值变换的配电网综合优化方法 |
CN114899825B (zh) * | 2022-04-20 | 2024-04-16 | 国家电网有限公司技术学院分公司 | 一种基于网络等值变换的配电网综合优化方法 |
CN114817865A (zh) * | 2022-06-28 | 2022-07-29 | 深圳市傲立电子有限公司 | 一种基于大数据的射频功率放大器用运行故障预测系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109103893B (zh) | 一种集群温控负荷参与电网agc的辅助调频方法 | |
CN109120005B (zh) | 一种多端柔性直流输电系统功率协调控制方法 | |
CN107104446B (zh) | 一种以全局最优为目标的电网自动电压控制方法及系统 | |
CN112952853A (zh) | 一种含分布式电源的配电网无功优化调度方法 | |
CN112769160B (zh) | 一种考虑并网和孤岛运行模式的微电网群自趋优控制方法 | |
CN112018790B (zh) | 基于分层分布式储能参与需求响应调节控制方法 | |
CN110311379A (zh) | 一种功能性模块化微电网组网及灵活调控方法 | |
CN112085327B (zh) | 一种配电商参与的主动配电网多层分区调控方法及系统 | |
CN109120018B (zh) | 基于一致性迭代算法的混合配电网分布式控制方法及系统 | |
CN112072672A (zh) | 含智能负荷的主动配电网优化调度方法 | |
CN114759620A (zh) | 一种风光储场站群无功协同优化调控方法、装置及系统 | |
CN110311403A (zh) | 基于分布式光伏发电的低压配电网电压功率综合控制方法 | |
CN113224769A (zh) | 考虑光伏多状态调节的多时间尺度配电网电压优化方法 | |
CN104659812B (zh) | 一种基于预测控制的多微网协调控制方法 | |
CN116565922A (zh) | 一种基于多微网互联运行结构的混合储能控制调度方法 | |
CN110350538B (zh) | 一种基于主动需求侧响应的微电网协调控制方法 | |
CN110137937B (zh) | 一种带储能设备的直流微网群系统及其分层控制方法 | |
CN112152268B (zh) | 一种交、直流子微网控制方法及子微网群间控制方法 | |
CN103094910A (zh) | 多级自动电压无功控制系统avc协调控制方法 | |
CN106410819B (zh) | 含电压不可行节点的受端电网动态无功优化方法及装置 | |
CN115967078A (zh) | 一种含分布式光伏的农村配电网电压分布式优化控制方法 | |
Tian et al. | An improved secondary voltage control method for half-wavelength transmission system | |
CN113988478A (zh) | 基于等微增率的直流微电网互联系统分散式经济优化方法 | |
Li et al. | A hierarchical multi-agent evaluation scheme for integrating distributed energy resources in distribution systems | |
CN108400583A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210611 |
|
RJ01 | Rejection of invention patent application after publication |