CN106953359B - 一种含分布式光伏配电网有功无功协调优化控制方法 - Google Patents
一种含分布式光伏配电网有功无功协调优化控制方法 Download PDFInfo
- Publication number
- CN106953359B CN106953359B CN201710265431.4A CN201710265431A CN106953359B CN 106953359 B CN106953359 B CN 106953359B CN 201710265431 A CN201710265431 A CN 201710265431A CN 106953359 B CN106953359 B CN 106953359B
- Authority
- CN
- China
- Prior art keywords
- time scale
- short
- optimizing control
- long time
- node
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000009826 distribution Methods 0.000 title claims abstract description 46
- 238000005457 optimization Methods 0.000 claims abstract description 64
- 238000005096 rolling process Methods 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000004146 energy storage Methods 0.000 claims description 30
- 230000003068 static effect Effects 0.000 claims description 24
- 239000003990 capacitor Substances 0.000 claims description 16
- 239000000243 solution Substances 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000002347 injection Methods 0.000 claims description 6
- 239000007924 injection Substances 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 239000004744 fabric Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000026676 system process Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- H02J3/383—
-
- 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
-
- 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
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
-
- 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
本发明提供一种含分布式光伏配电网有功无功协调优化控制方法,所述方法包括:S1、基于模型预测控制方法,根据不同时间尺度的控制目标和控制变量,将配电网系统控制过程分为长时间尺度优化控制、短时间尺度优化控制,并建立长时间尺度优化控制模型和短时间尺度优化控制模型;S2、将长时间尺度优化控制和短时间尺度优化控制的求解问题转化为二阶锥规划问题;S3、长时间尺度优化控制模型基于光伏出力及负荷需求预测信息,采用多步滚动优化求解各可控装置的有功无功出力,短时间尺度优化控制模型以长时间尺度的优化控制模型的求解结果为基准值,滚动求解各可控装置的有功无功出力增量。
Description
技术领域
本发明涉及电力系统控制领域,更具体地,涉及一种含分布式光伏配电网有功无功协调优化控制方法。
背景技术
近年来,随着分布式发电技术的不断发展,分布式光伏在配电网中的渗透率不断提高,对配电网的安全稳定运行提出了新的挑战。分布式光伏出力及负荷需求具有随机性与波动性,预测误差随预测时间的提前而增大,增加了系统优化控制的难度。
模型预测控制(Model Predictive Control,MPC)是解决这一问题的有效途径,它能够克服系统的不确定性,其在电力系统优化控制中的应用受到广泛重视。模型预测控制是一种基于模型的有限时域闭环优化控制算法,具有易于建模、控制效果良好、鲁棒性强等优点,能够有效应对系统的非线性、时变性和不确定性,广泛应用于炼油、化工、冶金、和电力等复杂工业过程控制中。模型预测控制由预测模型、滚动优化和反馈校正构成。与传统的开环优化一次下发全部优化指令的优化方式不同,模型预测控制采用滚动优化策略,其机理为:在每一采样时刻,根据当前系统状态及测量信息,基于预测模型对未来状态的预测结果,在线求解一个有限时长的最优控制问题,得到当前时刻和未来时段内的控制行为,并只执行当前时刻的控制行为,在下一个采样时刻,根据前一时刻控制后的系统状态及新的测量信息,重复上述过程。
但现有基于模型预测控制的电力系统优化控制研究都是对系统有功优化调度进行研究,并未考虑无功优化控制。配电网线路的R/X比值较大,有功无功的耦合性较强,二者对线路损耗、电压质量均有极大的影响,尤其是对于光伏渗透率较高的配电网,基于有功无功解耦的优化控制方法已不适用,需考虑有功无功协调优化控制。
发明内容
本发明提供一种克服上述问题或者至少部分地解决上述问题的含分布式光伏配电网有功无功协调优化控制方法,解决了现有控制方法无法适用于R/X比值较大、有功无功的耦合性较强的配电线路中,不能够对有功无功进行协调优化控制的问题,并通过多时间尺度协调控制结构应对光伏出力及负荷需求的波动性和随机性,减小网络损耗,保证系统的安全与经济运行。
根据本发明的一个方面,提供一种分布式光伏配电网有功无功协调优化控制方法,包括:
S1、基于模型预测控制方法,根据不同时间尺度的控制目标和控制变量,将系统控制过程分为长时间尺度优化控制、短时间尺度优化控制,并建立长时间尺度优化控制模型和短时间尺度优化控制模型;
S2、将长时间尺度优化控制和短时间尺度优化控制中的非凸、非线性约束条件采用二阶锥方法进行松弛处理,将长时间尺度优化控制模型和短时间尺度优化控制模型的求解问题转化为二阶锥规划求解问题;
S3、长时间尺度优化控制模型基于光伏出力及负荷需求预测信息,采用多步滚动优化求解各可控装置的有功无功出力,短时间尺度优化控制模型以长时间尺度的优化控制模型的求解结果为基准值,滚动求解各可控装置的有功无功出力增量。
作为优选的,所述步骤S1具体包括:
S11、基于模型预测控制方法,将控制过程划分为长时间尺度优化控制和短时间尺度优化控制,根据长时间尺度优化控制保证系统运行的经济性的优化目的,建立长时间尺度优化控制模型;根据短时间尺度优化控制保证系统运行的安全性的优化目的,建立短时间尺度优化控制模型;
S12、根据长时间尺度优化控制模型和短时间尺度优化控制模型的优化目的,选择长时间尺度优化控制模型和短时间尺度优化控制模型的优化目标;
S13、根据优化目标分别设定长时间尺度优化控制模型和短时间尺度优化控制模型的约束条件。
作为优选的,所述步骤S11具体包括:
基于模型预测控制方法,基于光伏出力及负荷需求的预测数据,以ΔT为时间间隔,建立长时间尺度优化控制模型,滚动求解未来M·ΔT时间段内系统中各可控装置的有功无功出力;
基于模型预测控制方法,根据当前系统运行状态及更小时间尺度的光伏出力及负荷需求预测数据,以Δt为时间间隔,建立短时间尺度优化控制模型,滚动求解未来N·Δt时间段内系统中各可控装置的有功无功出力;
所述M和N分别为长时间尺度优化控制模型和短时间尺度优化控制模型的控制步长,Δt<ΔT。
作为优选的,所述步骤S12包括:
所述长时间尺度优化控制模型中,以配电网系统网损最小为优化目标,并通过计算t0至t0+M·ΔT时间段内各支路线路的网损的总和计算配电网系统网损,t0为长时间尺度优化控制模型滚动优化起始时间;
所述短时间尺度优化控制模型中,以配电网系统网损最小为优化目标,并通过计算t0至t0+M·Δt时间段内各支路线路的网损的总和计算配电网系统网损。
作为优选的,所述长时间尺度优化控制模型中,各支路线路的网损的总和计算公式为:
iij,t=f(Kij,t,Hc,i,t,Pch,i,t,Pdis,i,t,QDG,i,t,QSVC,i,t)
所述短时间尺度优化控制模型中,各支路线路的网损的总和计算公式为:
式中,minF为系统网损最小值;t0为长时间尺度滚动优化起始时间,n为节点数,rij为线路ij的电阻,iij,t为t时刻线路ij的电流的平方,表示为各可控变量的函数;c(i)表示所有以i为线路首节点的线路末节点的集合;Kij,t有载调压变压器的可调档位;Hc,i,t为补偿电容器组调节档位;Pch,i,t、Pdis,i,t分别为储能装置的充放电功率;QDG,i,t、QSVC,i,t分别表示t时刻节点i处分布式光伏以及静止无功补偿装置的无功功率;表示t时刻储能装置充放电功率;分别表示分布式光伏及静止无功补偿装置无功出力长时间尺度优化控制的参考值;ΔPch,j,t、ΔPdis,j,t、ΔQDG,j,t、ΔQSVC,j,t分别表示t时刻短时间尺度优化控制储能装置充放电功率的调整值、分布式光伏及静止无功补偿装置无功出力的调整值。
作为优选的,所述步骤S13中包括潮流约束,所述潮流约束通过Distflow法进行潮流计算,对各支路线路节点电压幅值、电流幅值、注入的有功功率和无功功率、消耗的有功功率和无功功率进行约束计算。
作为优选的,所述长时间尺度优化控制模型中,潮流方程为:
所述短时间尺度优化控制模型中,潮流方程为:
式中,α(j)表示以节点i为首节点的支路末节点的集合;β(i)表示以节点j为首节点的支路末节点的集合;Pij,t、Qij,t分别为t时刻支路ij首端的有功功率和无功功率;Pj,t、Qj,t分别为t时刻节点j有功功率和无功功率的注入值;PDG,j,t、QDG,j,t分别为t时刻节点j处分布式光伏注入有功功率和无功功率;Pload,j,t、Qload,j,t分别为t时刻节点j处负荷消耗有功功率和无功功率;Pch,i,t、Pdis,i,t分别为储能装置的充放电功率;QSVC,j,t、Qc,j,t分别为t时刻节点j处静止无功补偿装置以及补偿电容器组的无功功率;rij+jxij为支路ij与该线路上有载调压变压器阻抗之和;ΔPch,j,t、ΔPdis,j,t、ΔQDG,j,t、ΔQSVC,j,t分别表示t时刻短时间尺度优化控制储能装置充放电功率的调整值、分布式光伏及静止无功补偿装置无功出力的调整值;Iij,t为t时刻流过支路ij的电流幅值;Vi,t为t时刻节点i的电压幅值;kij,t为t时刻支路ij处有载调压变压器的变比;vi,t、iij,t分别为t时刻节点i处电压幅值的平方以及支路ij电流幅值的平方。
作为优选的,所述步骤S13中还包括:
电压水平约束,用于限定各支路线路节点电压幅值的上下限;
支路容量约束,用于限定各支路线路电流幅值上限;
分布式光伏运行约束,用于根据支路线路节点处光伏有功无功出力的预测值及支路线路节点处的分布式光伏容量,限定该节点处分布式光伏无功出力的上下限;
有载调压变压器运行约束,用于根据有载调压变压器的标准变化和调节步长,限定有载调压变压器的可调档位及其上下限;
静止无功补偿装置运行约束,用于限定静止无功补偿装置的可调无功出力上下限;
补偿电容器组运行约束,用于根据补偿电容器组调节步长,控制补偿电容器组的运行;
储能装置运行约束,用于约束储能装置的充电限值,及充放电不同时发生。
作为优选的,所述步骤S2中,二阶锥规划的形式为:
式中,变量x∈RN;常量b∈RM,c∈RN,A∈RM×N;K为二阶锥或旋转二阶锥。
本申请提出一种含分布式光伏配电网有功无功协调优化控制方法,将系统的控制过程划分为长时间尺度优化控制及短时间尺度优化控制,不同时间尺度优化控制针对各自的控制目标及控制变量分别执行单独的模型预测控制,长时间尺度基于光伏出力及负荷需求预测信息,采用多步滚动优化求解各可控装置的有功无功出力,短时间尺度以长时间尺度的计算结果为基准值,滚动求解有功无功出力增量。能够有效应对光伏出力及负荷需求的波动性,综合考虑配电系统中的多种更可控装置,有效应对光伏接入造成的配电系统过电压现象,减小网络损耗,保证系统的安全经济运行;同时通过多时间尺度的协调,保证慢动作设备的调整次数不会过于频繁,进而保证了其的使用寿命,减少了配电网中设备购置开支;通过充分利用分布式光伏的无功出力,不对分布式光伏的有功进行削减,大大提高了配电网中分布式光伏的消纳,减少了系统从大电网中吸收有功及无功功率。长时间尺度及短时间尺度滚动优化控制中既包括连续变量又包括整数变量,其数学本质为混合整数非凸、非线性优化问题,难以求得最优解,本申请中将优化模型转化为可高效求解的二阶锥规划问题,使得计算更加简便高效。
附图说明
图1为根据本发明实施例的含分布式光伏配电网有功无功协调优化控制方法流程图;
图2为根据本发明实施例的基于MPC的控制结构示意图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
如图1和图2所示,图中示出了一种含分布式光伏配电网有功无功协调优化控制方法,包括:
S1、基于模型预测控制方法,根据不同时间尺度的控制目标和控制变量,将系统控制过程分为长时间尺度优化控制、短时间尺度优化控制,并建立长时间尺度优化控制模型和短时间尺度优化控制模型;
S2、将长时间尺度优化控制和短时间尺度优化控制中的非凸、非线性约束条件采用二阶锥方法进行松弛处理,将长时间尺度优化控制模型和短时间尺度优化控制模型的求解问题转化为二阶锥规划求解问题;
S3、长时间尺度优化控制模型基于光伏出力及负荷需求预测信息,采用多步滚动优化求解各可控装置的有功无功出力,短时间尺度优化控制模型以长时间尺度的优化控制模型的求解结果为基准值,滚动求解各可控装置的有功无功出力增量。
作为优选的,所述步骤S1具体包括:
S11、基于模型预测控制方法,将控制过程划分为长时间尺度优化控制和短时间尺度优化控制,根据长时间尺度优化控制保证系统运行的经济性的优化目的,建立长时间尺度优化控制模型;根据短时间尺度优化控制保证系统运行的安全性的优化目的,建立短时间尺度优化控制模型;
S12、根据长时间尺度优化控制模型和短时间尺度优化控制模型的优化目的,选择长时间尺度优化控制模型和短时间尺度优化控制模型的优化目标;
S13、根据优化目标分别设定长时间尺度优化控制模型和短时间尺度优化控制模型的约束条件。
作为优选的,所述步骤S11具体包括:
基于模型预测控制方法,基于光伏出力及负荷需求的预测数据,以ΔT为时间间隔,建立长时间尺度优化控制模型,滚动求解未来M·ΔT时间段内系统中各可控装置的有功无功出力;
基于模型预测控制方法,根据当前系统运行状态及更小时间尺度的光伏出力及负荷需求预测数据,以Δt为时间间隔,建立短时间尺度优化控制模型,滚动求解未来N·Δt时间段内系统中各可控装置的有功无功出力;
所述M和N分别为长时间尺度优化控制模型和短时间尺度优化控制模型的控制步长,Δt<ΔT。
作为优选的,所述步骤S12包括:
所述长时间尺度优化控制模型中,以配电网系统网损最小为优化目标,并通过计算t0至t0+M·ΔT时间段内各支路线路的网损的总和计算配电网系统网损,t0为长时间尺度优化控制模型滚动优化起始时间;
所述短时间尺度优化控制模型中,以配电网系统网损最小为优化目标,并通过计算t0至t0+M·Δt时间段内各支路线路的网损的总和计算配电网系统网损。
在本实施例中,为保证系统运行的经济性,减少网络损耗,长时间尺度优化控制的优化目标为系统网损最小,其计算公式为:
iij,t=f(Kij,t,Hc,i,t,Pch,i,t,Pdis,i,t,QDG,i,t,QSVC,i,t) (1)
式中,minF为系统网损最小值;t0为长时间尺度滚动优化起始时间,ΔT为长时间尺度优化控制的时间间隔,M为控制步长,n为节点数,rij为线路ij的电阻,iij,t为t时刻线路ij的电流的平方,表示为各可控变量的函数;c(i)表示所有以i为线路首节点的线路末节点的集合;Kij,t有载调压变压器的可调档位;Hc,i,t为补偿电容器组调节档位;Pch,i,t、Pdis,i,t分别为储能装置的充放电功率;QDG,i,t、QSVC,i,t分别表示t时刻节点i处分布式光伏以及静止无功补偿装置的无功功率。
S13、设定约束条件;长时间尺度优化控制模型的约束条件包括潮流约束、电压水平约束、支路容量约束、分布式光伏运行约束、有载调压变压器运行约束、静止无功补偿装置运行约束、补偿电容器组运行约束及储能装置运行约束;
具体的,针对辐射型配电网,所述潮流约束采用Distflow形式表示潮流方程:
α(j)表示以节点i为首节点的支路末节点的集合;β(i)表示以节点j为首节点的支路末节点的集合;Pij,t、Qij,t分别为t时刻支路ij首端的有功功率和无功功率;Pj,t、Qj,t分别为t时刻节点j有功功率和无功功率的注入值;PDG,j,t、QDG,j,t分别为t时刻节点j处分布式光伏注入有功功率和无功功率;Pload,j,t、Qload,j,t分别为t时刻节点j处负荷消耗有功功率和无功功率;Pch,i,t、Pdis,i,t分别为储能装置的充放电功率;QSVC,j,t、Qc,j,t分别为t时刻节点j处静止无功补偿装置以及补偿电容器组的无功功率;rij+jxij为支路ij与该线路上有载调压变压器阻抗之和;Iij,t为t时刻流过支路ij的电流幅值;Vi,t为t时刻节点i的电压幅值;kij,t为t时刻支路ij处有载调压变压器的变比;vi,t、iij,t分别为t时刻节点i处电压幅值的平方以及支路ij电流幅值的平方;
电压水平约束:
Vi min≤Vi,t≤Vi max (7)
式中,分别为节点i电压幅值上下限。
支路容量约束:
式中,为支路ij电流幅值上限。
分布式光伏运行约束:
式中,为t时刻节点i处分布式光伏有功出力的预测值, 分别为t时刻节点i处分布式光伏无功出力的上下限,SDG,i为节点i处分布式光伏的容量。
有载调压变压器(OLTC)运行约束:
kij,t=k0+Kij,tΔkij (12)
式中,k0为OLTC标准变比,Δkij为OLTC调节步长,Kij,t、 分别为t时刻OLTC可调档位及其上下限。
静止无功补偿装置(SVC)运行约束:
式中,分别为SVC可调无功出力上下限。
补偿电容器组运行约束:
式中,Hi,t为补偿电容器组可调档位,ΔQc,i,t为调节步长,n为最大可调档位。
储能装置运行约束:
式中,Ei,t表示储能装置的电量;Pch,i,t、Pdis,i,t分别表示储能装置的充放电功率;ηch、ηdis分别表示储能装置的充放电效率;Dch,i,t、Ddis,i,t为0-1变量,约束(18)保证储能装置的充放电不会同时发生;为储能装置的充电限值。
短时间尺度优化控制模型的约束条件包括潮流约束、电压水平约束、支路容量约束、分布式光伏运行约束、静止无功补偿装置运行约束及储能装置运行约束。由于光伏出力和负荷需求的随机波动大,长时间尺度的预测数据误差较大,无法满足精度要求,因此增加短时间尺度优化控制,基于系统当前运行状态及更短时间尺度的预测数据对长时间尺度的优化控制结果进行调整,以Δt为时间间隔,滚动求解未来N·Δt时间段内系统中各可控装置的有功无功出力,所述M和N分别为长时间尺度优化控制和短时间尺度优化控制的控制步长,Δt<ΔT;有载调压变压器和补偿电容器组响应速度较慢,且调整不宜过于频繁,因此在短时间尺度优化控制中不对二者进行调整,控制变量为储能装置充放电功率增量、静止无功补偿装置及分布式光伏无功出力增量。
由于光伏出力及负荷需求预测数据的误差随着预测时间的提前而增大,为应对二者的随机波动,防止出现电压越限,保证系统运行的安全性,同时为保证与长时间尺度整体优化方向及分析计算的一致性,短时间尺度优化控制仍以系统网损最小为优化目标,其计算公式为:
式中,minF为系统网损最小值;Δt为短时间尺度优化控制的时间间隔,N为短时间尺度优化控制的控制步长,表示t时刻储能装置充放电功率;分别表示分布式光伏及静止无功补偿装置无功出力长时间尺度优化控制的参考值;ΔPch,j,t、ΔPdis,j,t、ΔQDG,j,t、ΔQSVC,j,t分别表示t时刻短时间尺度优化控制储能装置充放电功率的调整值、分布式光伏及静止无功补偿装置无功出力的调整值。
所述短时间尺度优化控制中,约束变量包括潮流约束、电压水平约束、支路容量约束、分布式光伏运行约束、静止无功补偿装置运行约束及储能装置运行约束;具体的为:
潮流约束:
其他约束,具体包括电压水平约束、支路容量约束、分布式光伏运行约束、静止无功补偿装置运行约束及储能装置运行约束;
长时间尺度优化控制及短时间尺度优化控制均为基于预测信息的超前预测控制,由于光伏出力和负荷需求的不确定性,系统的实测值与预测值往往会出现偏差,因此通过反馈校正,将实测输出值以反馈信息反馈给长时间尺度优化控制及短时间尺度优化控制,作为下一轮滚动优化的初始值,实现闭环控制,使预测结果更接近实际值。
综合考虑系统中的有载调压变压器、补偿电容器、静止无功补偿装置、分布式光伏、储能装置等多种可控装置,将控制过程划分为长时间尺度优化控制和短时间尺度优化控制。长时间尺度优化控制保证系统运行的经济性,短时间尺度优化控制保证系统运行的安全性,二者基于各自的优化目标分别执行单独的模型预测控制。长时间尺度优化控制以系统网损最小为优化目标,基于光伏出力及负荷需求的预测数据,以ΔT为时间间隔,对未来M·ΔT时间内系统中各可控装置的有功无功出力进行优化求解,作为短时间尺度优化控制层的调节基点。短时间尺度优化控制根据当前系统运行状态及更小时间尺度的光伏出力及负荷需求预测数据,以Δt(Δt<ΔT)为时间间隔,滚动求解未来N·Δt时间内系统中各可控装置的有功无功出力增量,对上一级优化结果进行修正。控制结构如图2所示。
在本实施例中,优化问题中既包括连续变量又包括整数变量,其数学本质为混合整数非凸、非线性优化问题,难以求得最优解,本文考虑将优化模型转化为可高效求解的二阶锥规划(second-order cone programming,SOCP)问题。
所述步骤S2中,二阶锥规划的形式为:
式中,变量x∈RN;常量b∈RM,c∈RN,A∈RM×N;K为二阶锥或旋转二阶锥。
作为优选的,所述二阶锥形式为:
所述旋转二阶锥形式为:
以长时间尺度优化控制为例,潮流约束中的式(3)与(4)为非凸、非线性方程,优化目标函数及其他约束条件为线性方程,需对式(3)与(4)其进行处理。
式(3)与式(13)为对OLTC的约束,采用基于分段线性化的精确线性化建模方法对OLTC建模,将该约束转化为线性约束。
对式(4)采用二阶锥方法进行松弛,松弛之后式(4)可改写为:
将其改写为标准二阶锥形式,即:
||[2Pij,t 2Qij,t iij,t-vi,t]T||2≤iij,t+vi,t (26)
长时间尺度优化控制模型中包含离散变量及连续变量,优化模型最终转化为混合整数二阶锥规划模型;短时间尺度优化控制模型中只包含连续变量,最终转换为二阶锥规划模型。
本实施例中还示出了一种分布式光伏配电网有功无功协调优化系统,包括长时间尺度优化控制模型、短时间尺度优化控制模型;
所述长时间尺度优化控制模型用于以系统网损最小为优化目标,基于光伏出力及负荷需求的预测数据,以ΔT为时间间隔,滚动求解未来M·ΔT时间内各可控装置的有功无功出力;
所述短时间尺度优化控制模型用于基于系统当前运行状态及更短时间尺度的预测数据对长时间尺度的优化控制结果进行调整,以Δt(Δt<ΔT)为时间间隔,滚动求解未来N·Δt时间内系统中各可控装置的有功无功出力增量。
综上所述,本申请提出一种含分布式光伏配电网有功无功协调优化控制方法,将系统的控制过程划分为长时间尺度优化控制及短时间尺度优化控制,不同时间尺度优化控制针对各自的控制目标及控制变量分别执行单独的模型预测控制,长时间尺度基于光伏出力及负荷需求预测信息,采用多步滚动优化求解各可控装置的有功无功出力,短时间尺度以长时间尺度的计算结果为基准值,滚动求解有功无功出力增量。能够有效应对光伏出力及负荷需求的波动性,综合考虑配电系统中的多种更可控装置,有效应对光伏接入造成的配电系统过电压现象,减小网络损耗,保证系统的安全经济运行;同时通过多时间尺度的协调,保证慢动作设备的调整次数不会过于频繁,进而保证了其的使用寿命,减少了配电网中设备购置开支;通过充分利用分布式光伏的无功出力,不对分布式光伏的有功进行削减,大大提高了配电网中分布式光伏的消纳,减少了系统从大电网中吸收有功及无功功率。长时间尺度及短时间尺度滚动优化控制中既包括连续变量又包括整数变量,其数学本质为混合整数非凸、非线性优化问题,难以求得最优解,本申请中将优化模型转化为可高效求解的二阶锥规划问题,使得计算更加简便高效。
最后,本申请的方法仅为较佳的实施方案,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种含分布式光伏配电网有功无功协调优化控制方法,其特征在于,包括:
S1、基于模型预测控制方法,根据不同时间尺度的控制目标和控制变量,将配电网系统控制过程分为长时间尺度优化控制、短时间尺度优化控制,并建立长时间尺度优化控制模型和短时间尺度优化控制模型;所述步骤S1具体包括:
S11、基于模型预测控制方法,将控制过程划分为长时间尺度优化控制和短时间尺度优化控制;根据长时间尺度优化控制保证系统运行的经济性的优化目的,建立长时间尺度优化控制模型;根据短时间尺度优化控制保证系统运行的安全性的优化目的,建立短时间尺度优化控制模型;
S12、根据长时间尺度优化控制模型和短时间尺度优化控制模型的优化目的,选择长时间尺度优化控制模型和短时间尺度优化控制模型的优化目标;
S13、根据优化目标分别设定长时间尺度优化控制模型和短时间尺度优化控制模型的约束条件;
S2、将长时间尺度优化控制和短时间尺度优化控制中的非凸、非线性约束条件采用二阶锥方法进行松弛处理,将长时间尺度优化控制模型和短时间尺度优化控制模型的求解问题转化为二阶锥规划求解问题;
S3、长时间尺度优化控制模型基于光伏出力及负荷需求预测信息,采用多步滚动优化求解各可控装置的有功无功出力,短时间尺度优化控制模型以长时间尺度的优化控制模型的求解结果为基准值,滚动求解各可控装置的有功无功出力增量。
2.根据权利要求1所述的含分布式光伏配电网有功无功协调优化控制方法,其特征在于,所述步骤S11具体包括:
基于模型预测控制方法,建立长时间尺度优化控制模型,基于光伏出力及负荷需求的预测数据,以ΔT为时间间隔,滚动求解未来M·ΔT时间段内系统中各可控装置的有功无功出力;
基于模型预测控制方法,建立短时间尺度优化控制模型,根据当前系统运行状态及更小时间尺度的光伏出力及负荷需求预测数据,以Δt为时间间隔,滚动求解未来N·Δt时间段内系统中各可控装置的有功无功出力;
所述M和N分别为长时间尺度优化控制模型和短时间尺度优化控制模型的控制步长,Δt<ΔT。
3.根据权利要求2所述的含分布式光伏配电网有功无功协调优化控制方法,其特征在于,所述步骤S12包括:
所述长时间尺度优化控制模型中,以配电网系统网损最小为优化目标,并通过计算t0至t0+M·ΔT时间段内各支路线路的网损的总和计算配电网系统网损,t0为长时间尺度优化控制模型滚动优化起始时间;
所述短时间尺度优化控制模型中,以配电网系统网损最小为优化目标,并通过计算t0至t0+M·Δt时间段内各支路线路的网损的总和计算配电网系统网损。
4.根据权利要求3所述的分布式光伏配电网有功无功协调优化控制方法,其特征在于:
所述长时间尺度优化控制模型中,各支路线路的网损的总和计算公式为:
iij,t=f(Kij,t,Hc,i,t,Pch,i,t,Pdis,i,t,QDG,i,t,QSVC,i,t)
所述短时间尺度优化控制模型中,各支路线路的网损的总和计算公式为:
式中,minF为系统网损最小值;t0为长时间尺度滚动优化起始时间,n为节点数,rij为线路ij的电阻,iij,t为t时刻线路ij的电流的平方,表示为各可控变量的函数;c(i)表示所有以i为线路首节点的线路末节点的集合;Kij,t有载调压变压器的可调档位;Hc,i,t为补偿电容器组调节档位;Pch,i,t、Pdis,i,t分别为储能装置的充放电功率;QDG,i,t、QSVC,i,t分别表示t时刻节点i处分布式光伏以及静止无功补偿装置的无功功率;表示t时刻储能装置充放电功率;分别表示分布式光伏及静止无功补偿装置无功出力长时间尺度优化控制的参考值;ΔPch,j,t、ΔPdis,j,t、ΔQDG,j,t、ΔQSVC,j,t分别表示t时刻短时间尺度优化控制储能装置充放电功率的调整值、分布式光伏及静止无功补偿装置无功出力的调整值。
5.根据权利要求4所述的分布式光伏配电网有功无功协调优化控制方法,其特征在于,所述步骤S13中包括潮流约束,所述潮流约束通过Distflow法进行潮流计算,对各支路线路节点电压幅值、电流幅值、注入的有功功率和无功功率、消耗的有功功率和无功功率进行约束计算。
6.根据权利要求5所述的分布式光伏配电网有功无功协调优化控制方法,其特征在于:
所述长时间尺度优化控制模型中,潮流方程为:
所述短时间尺度优化控制模型中,潮流方程为:
式中,α(j)表示以节点i为首节点的支路末节点的集合;β(i)表示以节点j为首节点的支路末节点的集合;Pij,t、Qij,t分别为t时刻支路ij首端的有功功率和无功功率;Pj,t、Qj,t分别为t时刻节点j有功功率和无功功率的注入值;PDG,j,t、QDG,j,t分别为t时刻节点j处分布式光伏注入有功功率和无功功率;Pload,j,t、Qload,j,t分别为t时刻节点j处负荷消耗有功功率和无功功率;Pch,i,t、Pdis,i,t分别为储能装置的充放电功率;QSVC,j,t、Qc,j,t分别为t时刻节点j处静止无功补偿装置以及补偿电容器组的无功功率;rij+jxij为支路ij与该线路上有载调压变压器阻抗之和;ΔPch,j,t、ΔPdis,j,t、ΔQDG,j,t、ΔQSVC,j,t分别表示t时刻短时间尺度优化控制储能装置充放电功率的调整值、分布式光伏及静止无功补偿装置无功出力的调整值;Iij,t为t时刻流过支路ij的电流幅值;Vi,t为t时刻节点i的电压幅值;kij,t为t时刻支路ij处有载调压变压器的变比;vi,t、iij,t分别为t时刻节点i处电压幅值的平方以及支路ij电流幅值的平方。
7.根据权利要求4所述的分布式光伏配电网有功无功协调优化控制方法,其特征在于,所述步骤S13中还包括:
电压水平约束,用于限定各支路线路节点电压幅值的上下限;
支路容量约束,用于限定各支路线路电流幅值上限;
分布式光伏运行约束,用于根据支路线路节点处光伏有功无功出力的预测值及支路线路节点处的分布式光伏容量,限定该节点处分布式光伏无功出力的上下限;
有载调压变压器运行约束,用于根据有载调压变压器的标准变化和调节步长,限定有载调压变压器的可调档位及其上下限;
静止无功补偿装置运行约束,用于限定静止无功补偿装置的可调无功出力上下限;
补偿电容器组运行约束,用于根据补偿电容器组调节步长,控制补偿电容器组的运行;
储能装置运行约束,用于约束储能装置的充电限值,及充放电不同时发生。
8.根据权利要求1所述的分布式光伏配电网有功无功协调优化控制方法,其特征在于,所述步骤S2中,二阶锥规划的形式为:
式中,变量x∈RN;常量b∈RM,c∈RN,A∈RM×N;K为二阶锥或旋转二阶锥。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710265431.4A CN106953359B (zh) | 2017-04-21 | 2017-04-21 | 一种含分布式光伏配电网有功无功协调优化控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710265431.4A CN106953359B (zh) | 2017-04-21 | 2017-04-21 | 一种含分布式光伏配电网有功无功协调优化控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106953359A CN106953359A (zh) | 2017-07-14 |
CN106953359B true CN106953359B (zh) | 2019-08-27 |
Family
ID=59476474
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710265431.4A Expired - Fee Related CN106953359B (zh) | 2017-04-21 | 2017-04-21 | 一种含分布式光伏配电网有功无功协调优化控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106953359B (zh) |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107732982B (zh) * | 2017-10-20 | 2021-02-26 | 河海大学 | 考虑模型预测控制的综合能源系统多时间尺度调度方法 |
CN107968445B (zh) * | 2017-12-07 | 2019-09-06 | 华中科技大学 | 一种含分布式电源的配网无功优化方法 |
CN108711866B (zh) * | 2018-04-10 | 2023-06-13 | 国网安徽省电力有限公司芜湖供电公司 | 一种用于新能源电站无功电压的控制系统 |
CN109347137A (zh) * | 2018-09-21 | 2019-02-15 | 中国农业大学 | 一种分布式光伏集群协调优化控制方法及装置 |
CN109638874A (zh) * | 2018-10-24 | 2019-04-16 | 中国电力科学研究院有限公司 | 一种分布式光伏集群控制方法和装置 |
CN109638873A (zh) * | 2018-10-24 | 2019-04-16 | 中国电力科学研究院有限公司 | 一种分布式光伏集群优化调度方法和系统 |
CN109193715A (zh) * | 2018-10-24 | 2019-01-11 | 云南电网有限责任公司电力科学研究院 | 一种主动配电网的功率控制方法及系统 |
CN109474015B (zh) * | 2018-12-17 | 2020-11-17 | 国家电网有限公司 | 基于两层优化的主动配电网储能系统优化配置方法及系统 |
CN110071512B (zh) * | 2019-03-13 | 2022-07-05 | 国网浙江平阳县供电有限责任公司 | 储能装置平抑光伏输出功率波动的控制方法 |
CN110059897B (zh) * | 2019-05-23 | 2021-03-09 | 合肥工业大学 | 基于混合整数pso算法的有源配电网日内滚动优化方法 |
CN110504709B (zh) * | 2019-08-27 | 2021-05-18 | 国网河北省电力有限公司邢台供电分公司 | 光伏集群无功电压调控方法、终端设备及存储介质 |
CN110581571A (zh) * | 2019-08-29 | 2019-12-17 | 昆明理工大学 | 一种主动配电网动态优化调度方法 |
CN112787322B (zh) * | 2019-10-23 | 2024-01-23 | 滕欣元 | 基于scada系统和多时间尺度的电网动态管理方法 |
CN111106631B (zh) * | 2020-01-14 | 2021-03-12 | 清华大学 | 配电网的分布式无功调度方法、系统、设备及存储介质 |
CN111864800B (zh) * | 2020-08-10 | 2022-05-13 | 华中科技大学 | 基于变换器并网装备的多尺度惯量控制方法及附加控制器 |
CN112152226B (zh) * | 2020-08-28 | 2022-08-02 | 华北电力科学研究院有限责任公司 | 基于分布式光伏节点的调压方法及装置 |
CN112736913B (zh) * | 2020-12-29 | 2023-03-03 | 郑州轻工业大学 | 含分布式电源的配电网功率优化模式影响因素分析方法 |
CN113224769B (zh) * | 2021-02-20 | 2024-04-12 | 北京交通大学 | 考虑光伏多状态调节的多时间尺度配电网电压优化方法 |
CN113300412A (zh) * | 2021-05-11 | 2021-08-24 | 广东电网有限责任公司中山供电局 | 一种配电网分布式光伏优化调度方法及装置 |
CN113783226B (zh) * | 2021-08-27 | 2023-07-25 | 中国科学院广州能源研究所 | 一种海上风电制氢并网发电系统的分层预测能量管理方法 |
CN113890016B (zh) * | 2021-09-26 | 2022-09-23 | 天津大学 | 数据驱动的配电网多时间尺度电压协调控制方法 |
CN113807029B (zh) * | 2021-10-19 | 2022-07-29 | 华北电力大学(保定) | 基于深度强化学习的双时间尺度新能源电网电压优化方法 |
CN114006379B (zh) * | 2021-11-25 | 2023-08-15 | 华北电力大学(保定) | 基于动-静无功置换机制的双时间尺度电压控制方法 |
CN114865617A (zh) * | 2022-04-12 | 2022-08-05 | 中国电建集团华东勘测设计研究院有限公司 | 基于大网拆小网的分布式潮流优化方法 |
CN117293807A (zh) * | 2023-09-25 | 2023-12-26 | 上海能优网电力科技有限公司 | 一种配电网信息侧模型多时间尺度优化方法及系统 |
CN117117973B (zh) * | 2023-10-24 | 2024-01-12 | 国网浙江省电力有限公司宁波供电公司 | 基于时间尺度的分布式电源调度方法、装置及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103730900A (zh) * | 2013-11-13 | 2014-04-16 | 广西电网公司河池供电局 | 电力系统多时间尺度的省地县一体化无功优化方法 |
CN106487042A (zh) * | 2016-11-22 | 2017-03-08 | 合肥工业大学 | 一种多时间尺度微电网电压无功优化控制方法 |
-
2017
- 2017-04-21 CN CN201710265431.4A patent/CN106953359B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103730900A (zh) * | 2013-11-13 | 2014-04-16 | 广西电网公司河池供电局 | 电力系统多时间尺度的省地县一体化无功优化方法 |
CN106487042A (zh) * | 2016-11-22 | 2017-03-08 | 合肥工业大学 | 一种多时间尺度微电网电压无功优化控制方法 |
Non-Patent Citations (1)
Title |
---|
《基于模型预测控制的主动配电网多时间尺度动态优化调度》;董雷等;《中国电机工程学报》;20160905;第36卷(第17期);第4610页右栏第2段至第4613页左栏第1段 * |
Also Published As
Publication number | Publication date |
---|---|
CN106953359A (zh) | 2017-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106953359B (zh) | 一种含分布式光伏配电网有功无功协调优化控制方法 | |
CN108683179B (zh) | 基于混合整数线性规划的主动配电网优化调度方法及系统 | |
Mohamed et al. | Real-time energy management scheme for hybrid renewable energy systems in smart grid applications | |
CN106505579B (zh) | 一种电动汽车参与配电网电压调节的调度控制方法 | |
CN105914753B (zh) | 主动配电网电压实时模糊控制方法 | |
CN110518575A (zh) | 基于区域划分的多时间尺度主动配电网电压优化控制方法 | |
CN110472824A (zh) | 一种考虑调峰需求的梯级水电站短期多目标优化调度方法 | |
CN110581571A (zh) | 一种主动配电网动态优化调度方法 | |
Maknouninejad et al. | Cooperative control for self-organizing microgrids and game strategies for optimal dispatch of distributed renewable generations | |
CN109873447B (zh) | 一种多源协同主动配电网多时间级有功-无功调控方法 | |
CN109638873A (zh) | 一种分布式光伏集群优化调度方法和系统 | |
CN103441506A (zh) | 不同时间尺度下多目标协调分散风电场无功优化控制方法 | |
Zou | Design of reactive power optimization control for electromechanical system based on fuzzy particle swarm optimization algorithm | |
CN106253352A (zh) | 计及风电概率分布特征的鲁棒实时调度方法 | |
CN105896596B (zh) | 一种考虑需求侧响应的风电功率分层平滑系统及其方法 | |
CN104701858B (zh) | 计及分区动态无功储备的无功电压控制方法 | |
CN108599259A (zh) | 一种基于灵敏度分析的微电网主动运行决策方法 | |
CN117060470B (zh) | 一种基于灵活性资源的配电网电压优化控制方法 | |
CN105939017B (zh) | 考虑时段间耦合的无功优化的工程实用化求解方法 | |
CN111667136A (zh) | 一种区域电力市场的出清方法、装置及存储介质 | |
CN110048438A (zh) | 一种基于模型预测控制的配电网馈线级负荷功率控制方法 | |
CN105680444A (zh) | 一种输配电线路功率因数控制分析方法 | |
CN114069635A (zh) | 一种考虑灵活资源经济性的电网调控优化方法 | |
CN115036931A (zh) | 一种有源电网无功电压仿射可调鲁棒优化方法和装置 | |
CN103730900A (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 | ||
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: 20190827 |