CN106921164B - 配网电压无功协同优化的混合整数二阶锥规划方法与*** - Google Patents
配网电压无功协同优化的混合整数二阶锥规划方法与*** Download PDFInfo
- Publication number
- CN106921164B CN106921164B CN201710217565.9A CN201710217565A CN106921164B CN 106921164 B CN106921164 B CN 106921164B CN 201710217565 A CN201710217565 A CN 201710217565A CN 106921164 B CN106921164 B CN 106921164B
- Authority
- CN
- China
- Prior art keywords
- power
- distribution network
- voltage
- photo
- power supply
- 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
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000006243 chemical reaction Methods 0.000 claims abstract description 21
- 238000010248 power generation Methods 0.000 claims description 62
- 230000007613 environmental effect Effects 0.000 claims description 27
- 239000003990 capacitor Substances 0.000 claims description 11
- 238000004146 energy storage Methods 0.000 claims description 10
- 230000005611 electricity Effects 0.000 claims description 9
- 238000010276 construction Methods 0.000 claims description 6
- 230000009977 dual effect Effects 0.000 claims description 5
- 230000007246 mechanism Effects 0.000 abstract description 5
- 230000008569 process Effects 0.000 abstract description 4
- 238000013486 operation strategy Methods 0.000 abstract 1
- 238000012544 monitoring process Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000005855 radiation Effects 0.000 description 3
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 2
- 102000016550 Complement Factor H Human genes 0.000 description 2
- 108010053085 Complement Factor H Proteins 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007596 consolidation process Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
-
- 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]
-
- H02J3/383—
-
- 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
本发明提供一种配网电压无功协同优化的混合整数二阶锥规划方法与***,获取主动配电网中光伏电源输出功率,获取主动配电网的电压无功多目标协同优化的目标函数,确定主动配电网安全运行需满足的约束条件,构建主动配电网的电压无功多目标协同优化模型,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,求解混合整数二阶锥规划模型,获得决策变量的运行策略,并进行协同优化。整个过程中,基于光伏电源输出功率,考虑光伏电源出力的不确定性对配电网电压分布和无功潮流的影响机理,并且基于混合整数二阶锥规划的请求结果对配网电压无功协同优化的混合整数二阶锥规划,能够确保优化的最优性和准确度,确保电力电网的安全运行。
Description
技术领域
本发明涉及电力电网技术领域,特别是涉及配网电压无功协同优化的混合整数二阶锥规划方法与***。
背景技术
电力需求的持续增长、传统能源的短缺以及电力市场的开放正驱动电网朝高效、灵活、智能和可持续方式发展,以适应未来的技术需求。可持续性是未来电网的基础特征,其本质表现为分布式电源尤其是可再生能源的规模化接入与应用。
分布式发电近年来发展迅速,特别是光伏发电发展迅速,分布式光伏电源对配电网电压的影响将从两方面考虑。首先接入分布式电源后,在稳态情况下,会使得沿馈线的各负荷节点处的电压有所提高,甚至某些节点的电压会超过上限,对用户产生影响,但是分布式电源的退出,使得依靠分布式电源支撑的馈线电压下跌,带来电能质量问题,所以分布式电源接入配电网要进行调压。其次分布式电源接入后能在一定程度上改善配电网的无功不足,但是有些分布式电源不但不能提供无功功率,反而会消耗无功功率,所以分布式电源接入配电网也需要进行一定的无功补偿。
在实际运行中,配网电压还会受到负荷不断变动的影响,为了避免运行电压过高或过低,无功补偿不足或过剩,就必须在不同负荷水平下相应调整电压无功控制策略,以使配电网中的电压和无功功率得到合理控制和优化。
发明内容
基于此,有必要针对目前尚无针对包含光伏电源的配网电压无功协同优化的混合整数二阶锥规划方法的问题,提供一种配网电压无功协同优化的混合整数二阶锥规划方法与***。
一种配网电压无功协同优化的混合整数二阶锥规划方法,应用于包含光伏电源的主动配电网,包括步骤:
获取所述主动配电网中光伏电源输出功率;
获取所述主动配电网的电压无功多目标协同优化的目标函数;
获取主动配电网运行数据,确定主动配电网安全运行需满足的约束条件;
根据所述光伏电源输出功率、所述目标函数以及所述约束条件,构建主动配电网的电压无功多目标协同优化模型;
引入中间变量,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,获得混合整数二阶锥规划模型;
求解所述混合整数二阶锥规划模型,获得决策变量的运行策略;
根据所述决策变量的运行策略,对主动配电网电压无功协同进行优化。
一种配网电压无功协同优化的混合整数二阶锥规划***,应用于包含光伏电源的主动配电网,包括:
输出功率获取模块,用于获取所述主动配电网中光伏电源输出功率;
目标函数获取模块,用于获取所述主动配电网的电压无功多目标协同优化的目标函数;
约束条件确定模块,用于获取主动配电网运行数据,确定主动配电网安全运行需满足的约束条件;
模型构建模块,用于根据所述光伏电源输出功率、所述目标函数以及所述约束条件,构建主动配电网的电压无功多目标协同优化模型;
模型转换模块,用于引入中间变量,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,获得混合整数二阶锥规划模型;
求解模块,用于求解所述混合整数二阶锥规划模型,获得决策变量的运行策略;
优化模块,用于根据所述决策变量的运行策略,对主动配电网电压无功协同进行优化。
本发明配网电压无功协同优化的混合整数二阶锥规划方法与***,获取所述主动配电网中光伏电源输出功率,获取所述主动配电网的电压无功多目标协同优化的目标函数,确定主动配电网安全运行需满足的约束条件,构建主动配电网的电压无功多目标协同优化模型,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,求解所述混合整数二阶锥规划模型,获得决策变量的运行策略,并进行协同优化。整个过程中,基于光伏电源输出功率,考虑光伏电源出力的不确定性对配电网电压分布和无功潮流的影响机理,并且基于混合整数二阶锥规划的请求结果对配网电压无功协同优化的混合整数二阶锥规划,能够确保优化的最优性和准确度,确保电力电网的安全运行。
附图说明
图1为本发明配网电压无功协同优化的混合整数二阶锥规划方法其中一个实施例的流程示意图;
图2为主动配电网***结构示意图;
图3为本发明配网电压无功协同优化的混合整数二阶锥规划方法其中一个应用实施例的流程示意图;
图4为本发明配网电压无功协同优化的混合整数二阶锥规划***其中一个实施例的结构示意图。
具体实施方式
太阳能是一种巨量的可再生能源,在化石能源日益枯竭的世纪,光伏发电必将从太空能源、应急能源、军用能源、便携式能源、扶贫能源,最终发展成为最重要的、可再生的替代能源之一。随着光伏并网技术的发展,近年来越来越多的光伏电站并网运行,在向电网馈送电能的基础上,也具备了一定的无功调节和辅助服务能力。但是,太阳能光伏电源发电量与太阳照射强度有关,由于到达某一地面的太阳直接福射受气候、季节等因素的影响,所以太阳能光伏发电是不稳定的。光伏电源并网后,电网的稳定,安全,经济运行和供电质量等将会受到一定的影响。配电网中除了通过投切电容器调节电外,一般很少具有其他的动态无功调节设备,如果光伏电站所发电量比例较大,则由于其易变性将会使配电网线路中的潮流变化较大且极易波动,从而使得电网正常运行时的电压调整难度加大了,调节不好可能会使电正质量不合格。如果将大量的光伏电站接入到配电网的终端,配电网中将会产生反向潮流,这样由于反向潮流所产生的压降将可能使得负荷侧的电压越限。
在数学上,电压无功优化是典型的非线性规划问题,具有非线性、小连续、不确定因素较多等特点。对于电压无功优化的研究方法,传统的数学规划方法主要有非线性规划法和线性规划法等。采用常规算法求解电压无功优化问题时遇到的主要困难就是离散变量的归整问题和多峰多极值问题。二阶锥规划是在有限个二阶锥的笛卡尔乘积的仿射子空间的交集上最小化或最大化一个线性函数,即在非空尖凸锥导入偏序下,线性等式、线性不等式约束条件下的线性目标函数的问题。通过将复杂的优化模型转化成锥模型,可以把变量间复杂的关系以特殊结构的锥集表示,大大简化了原模型的求解,加快收敛速度。
基于上述实际应用情况以及科学理论依据,本发明提供一种配网电压无功协同优化的混合整数二阶锥规划方法,应用于包含光伏电源的主动配电网,如图1所示,其具体包括步骤:
S100:获取所述主动配电网中光伏电源输出功率。
主动配电网中光伏电源输出功率可以从历史数据分析获得,或者从光伏发电监控中心直接获取,又或者,基于当前采集的到参数进行处理与计算获得。具体来说,光伏电源输出功率需考虑到各节点位于不同地理区域,由于各区域的环境温度、日照强度、日照时间、日照阴影、日照偏角等因素造成光伏电源出力存在不确定性于随机性。
S200:获取所述主动配电网的电压无功多目标协同优化的目标函数。
具体来说,主动配电网的电压无功多目标协同优化的目标函数F为:
F=min(αf1+βf2)
式中,f1、f2分别为子目标1配电网有功网损和子目标2电压偏移最优值程度,α、β分别为两个子目标的权重系数,且α+β=1。更进一步来说,目标函数的子目标f1、f2的计算步骤分别为:
1、子目标f1为配电网有功网损,计算公式为:
式中,Vi为节点i的电压幅值,Vj为节点j的电压幅值,θij为支路ij间的电压相位角度。
2、子目标f2为电压偏移最优值程度,计算公式为:
式中,VN为节点N电压幅值。
S300:获取主动配电网运行数据,确定主动配电网安全运行需满足的约束条件。
主动配电网运行数据具体可以包括节点上的负荷功率最大值和最小值S Di,分布式电源接入功率的最大值和最小值S DGi,无功补偿电容器组投切容量的最大值和最小值Q Ci S DGi,储能设备吸放功率的最大值和最小值S DSi,外用用电设备(电动汽车)充放电功率的最大值和最小值S EVi,电压幅值的最大值和最小值V i。更具体来说,电网安全运行需满足的约束条件为:
SDi+SDSi+SEVi+LLi=SDGi+QCi+SBi
S400:根据所述光伏电源输出功率、所述目标函数以及所述约束条件,构建主动配电网的电压无功多目标协同优化模型。
基于步骤S100获得的光伏电源输出功率,步骤S200获得的目标函数以及步骤S300确定的约束条件,构建主动配电网的电压无功多目标协同优化模型,考虑采用数学模型的方式来对主动配电网的电压无功多目标进行协同优化。
S500:引入中间变量,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,获得混合整数二阶锥规划模型。
引入中间变量,将主动配电网的电压无功多目标协同优化模型求解的问题转换为混合整数二阶锥规划问题,由于二阶锥规划是在有限个二阶锥的笛卡尔乘积的仿射子空间的交集上最小化或最大化一个线性函数,即在非空尖凸锥导入偏序下,线性等式、线性不等式约束条件下的线性目标函数的问题。通过将复杂的优化模型转化成锥模型,可以把变量间复杂的关系以特殊结构的锥集表示,大大简化了原模型的求解。具体来说,可以引入中间变量,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,将搜索空间限制在有限的凸锥范围内,获得混合整数二阶锥规划模型。更进一步来说,引入的中间变量为:
uij=ViVj sinθij
vij=ViVjcosθij。
式中,Vi为节点i的电压幅值,Vj为节点j的电压幅值,θij为支路ij间的电压相位角度。在这里引入中间变量之后,过二阶锥转换后主动配电网的电压无功多目标协同优化的目标函数F'为:
F'=min(αf1'+βf2')
式中,f1'和f2'分别为经过二阶锥转换后的子目标1配电网有功网损和子目标2电压偏移最优值程度。
S600:求解所述混合整数二阶锥规划模型,获得决策变量的运行策略。
对混合整数二阶锥规划模型求解,获得决策变量的运行策略。具体来说,决策变量的运行策略包括随光伏电源出力变化的无功补偿电容器组的投切容量QCi、储能设备的吸放功率SDSi以及外部用电设备(电动汽车)的充放电功率SEVi等。更具体来说,可以采用MOSEK7.1求解器利用分支定界-原对偶内点法来求解计算优化后的决策变量。
S700:根据所述决策变量的运行策略,对主动配电网电压无功协同进行优化。
本发明配网电压无功协同优化的混合整数二阶锥规划方法,获取所述主动配电网中光伏电源输出功率,获取所述主动配电网的电压无功多目标协同优化的目标函数,确定主动配电网安全运行需满足的约束条件,构建主动配电网的电压无功多目标协同优化模型,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,求解所述混合整数二阶锥规划模型,获得决策变量的运行策略,并进行协同优化。整个过程中,基于光伏电源输出功率,考虑光伏电源出力的不确定性对配电网电压分布和无功潮流的影响机理,并且基于混合整数二阶锥规划的请求结果对配网电压无功协同优化的混合整数二阶锥规划,能够确保优化的最优性和准确度,确保电力电网的安全运行。
在其中一个实施例中,所述获取所述主动配电网中光伏电源输出功率的步骤包括:
步骤一:获取光伏电源所在地预设周期内各时段的相关环境数据信息,计算光伏电源中光伏发电板的槽温。
从公共气象数据平台获取光伏电源所在地一定周期内各时段的相关环境数据信息,即环境温度TAi、平均光照辐射强度sai,计算光伏发电板的槽温TCi为:
式中,TOT为光伏发电板正常工作温度,单位为℃。
步骤二:获取光伏电源接入点的短路电流值,并根据所述光伏发电板的槽温和所述短路电流值,计算光伏发电板在当前环境条件下的工作电流。
从光伏发电监控中心获取光伏电源接入点的短路电流值ISCi,计算光伏发电板在该环境条件下的工作电流IPVi为:
IPVi=sai[ISCi+Ki(TCi-25)]
式中,KIi为电流温度系数,单位为A/℃。
步骤三:获取光伏电源接入点的开路电压值,并根据所述光伏发电板的槽温和所述开路电压值,计算光伏发电板在当前环境条件下的工作电压。
从光伏发电监控中心获取光伏电源接入点的开路电压值VOCi,计算光伏发电板在该环境条件下的工作电压VPVi为:
VPVi=VOCi-KVi·TCi
式中,KVi为电压温度系数,单位为V/℃。
步骤四:获取光伏发电板的最大功率点电流和最大功率点电压,根据光伏电源接入点的短路电流值、光伏电源接入点的开路电压值以及光伏发电板的最大功率点电流和最大功率点电压,计算光伏电源填充因数。
从光伏发电监控中心获取光伏发电板最大功率点电流IMPPi和光伏发电板最大功率点电压VMPPi,计算光伏电源填充因数H为:
步骤五:根据光伏发电板在当前环境条件下的工作电流、光伏发电板在当前环境条件下的工作电压以及光伏电源填充因数,获取光伏电源输出功率。
根据光伏发电***输出功率与环境因素之间的特征关系,采用光伏发电模型模拟的方法计算确定光伏发电***输出功率PPV为:
PPV=N·H·VPVi·IPVi
式中,N为光伏发电板的总光照面积。
为更进一步详细解释本发明配网电压无功协同优化的混合整数二阶锥规划方法的技术方案及其实现的效果,下面将采用一具体应用实例,并集合附图2以及附图3进行详细的说明。
如图2为主动配电网结构示意图,在图2中1为110kV变电站主变压器高压侧母线;2为主变压器高压侧功率;3为主变压器阻抗;4为主变压器低压侧功率;5为主变压器损耗;6为110kV变电站主变压器低压侧母线;7为主变压器低压侧母线的10kV侧有1条配电支线,其注入功率为;8为主变压器低压侧母线上接入的无功补偿电容器组;9为主变压器低压侧母线上接入的储能设备;10为主变压器低压侧母线上的负荷;11为主变压器低压侧母线上接入的分布式电源;12为主变压器低压侧母线上接入的电动汽车;13为10kV配电线路上的第i个节点;14为节点i上有1条配电支线,其注入功率为;15为节点i上接入的无功补偿电容器组;16为节点i上接入的储能设备;17为节点i上接入的负荷;18为节点i上接入的分布式电源;19为节点i上接入的电动汽车;20为配电线路ij的首端功率;21为配电线路ij上的线路损耗;22为配电线路ij的末端功率;23为配电线路末端。其中i=1,2,…,n。
在具体应用实例中,主动配电网为10kV电压等级的配电***,且由n个负荷节点组成,通过一条馈线与上一级变压器低压侧连接,假设第i个节点上的负荷功率、分布式电源接入功率、无功补偿电容器组的投切容量、储能设备的吸放功率、外部用电设备(电动汽车)的充放电功率、电压幅值分别为SDi、SDGi、QCi、SDSi、SEVi、Vi。第i个节点上的负荷功率最大值和最小值分别为 S Di,分布式电源接入功率的最大值和最小值分别为 S DGi,无功补偿电容器组投切容量的最大值和最小值分别为 Q Ci,储能设备吸放功率的最大值和最小值分别为 S DSi,电动汽车充放电功率的最大值和最小值别为 S EVi,电压幅值的最大值和最小值分别为 V i。
如图3所示,在该具体应用实例中,本发明配网电压无功协同优化的混合整数二阶锥规划方法包括步骤:
S1、从公共气象数据平台获取光伏电源所在地一定周期内各时段的相关环境数据信息,即环境温度TAi、平均光照辐射强度sai,计算光伏发电板的槽温TCi为:
式中,TOT为光伏发电板正常工作温度,单位为℃。
S2、从光伏发电监控中心获取光伏电源接入点的短路电流值ISCi,计算光伏发电板在该环境条件下的工作电流IPVi为:
IPVi=sai[ISCi+Ki(TCi-25)]
式中,KIi为电流温度系数,单位为A/℃。
S3、从光伏发电监控中心获取光伏电源接入点的开路电压值VOCi,计算光伏发电板在该环境条件下的工作电压VPVi为:
VPVi=VOCi-KVi·TCi
式中,KVi为电压温度系数,单位为V/℃。
S4、从光伏发电监控中心获取光伏发电板最大功率点电流IMPPi和光伏发电板最大功率点电压VMPPi,计算光伏电源填充因数H为:
S5、根据光伏发电***输出功率与环境因素之间的特征关系,采用光伏发电模型模拟的方法计算确定光伏发电***输出功率PPV为:
PPV=N·H·VPVi·IPVi
式中,N为光伏发电板的总光照面积。
S6、建立主动配电网电压无功多目标协同优化的目标函数。主动配电网的电压无功多目标协同优化的目标函数F为:
F=min(αf1+βf2)
式中,f1、f2分别为子目标1配电网有功网损和子目标2电压偏移最优值程度,α、β分别为两个子目标的权重系数,且α+β=1。
S7、考虑电网安全运行需要满足的约束条件。电网安全运行需满足的约束条件为:
SDi+SDSi+SEVi+LLi=SDGi+QCi+SBi
S8、引入中间变量的计算。
引入的中间变量为:
uij=ViVjsinθij
vij=ViVjcosθij。
式中,Vi为节点i的电压幅值,Vj为节点j的电压幅值,θij为支路ij间的电压相位角度。
S9、引入中间变量之后,过二阶锥转换后主动配电网的电压无功多目标协同优化的目标函数F'为:
F'=min(αf1'+βf2')
式中,f1'和f2'分别为经过二阶锥转换后的子目标1配电网有功网损和子目标2电压偏移最优值程度。
S10、可以采用MOSEK7.1求解器利用分支定界-原对偶内点法来求解计算优化后的决策变量。包括随光伏电源出力变化的无功补偿电容器组的投切容量QCi、储能设备的吸放功率SDSi以及外部用电设备(电动汽车)的充放电功率SEVi。
具体来说,在其中一个应用实例中,本发明配网电压无功协同优化的混合整数二阶锥规划方法是一种基于混合整数二阶锥规划的配网电压无功协同优化的混合整数二阶锥规划的方法,可以计算出在调度周期内,分布式电源环境下配电网中无功补偿电容器组的投切容量、储能设备的吸放功率、电动汽车的充放电功率。反映了环境温度、日照强度、日照时间、日照阴影、日照偏角等因素对光伏发电***输出功率的影响机理,反映了分布式光伏发电***出力的不确定性对配电网电压分布和无功潮流的影响机理,反映了混合整数二阶锥规划问题利用严格的数学理论进行建模,将搜索空间限制在有限的凸锥范围内,使其不仅能在有效的时间内对原问题的求解,还能保证所求解的最优性和精度,为电网安全运行和提高电能质量提供了一种可靠的运行策略。
如图4所示,一种配网电压无功协同优化的混合整数二阶锥规划***,应用于包含光伏电源的主动配电网,包括:
输出功率获取模块100,用于获取所述主动配电网中光伏电源输出功率。
目标函数获取模块200,用于获取所述主动配电网的电压无功多目标协同优化的目标函数。
约束条件确定模块300,用于获取主动配电网运行数据,确定主动配电网安全运行需满足的约束条件。
模型构建模块400,用于根据所述光伏电源输出功率、所述目标函数以及所述约束条件,构建主动配电网的电压无功多目标协同优化模型。
模型转换模块500,用于引入中间变量,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,获得混合整数二阶锥规划模型。
求解模块600,用于求解所述混合整数二阶锥规划模型,获得决策变量的运行策略。
优化模块700,用于根据所述决策变量的运行策略,对主动配电网电压无功协同进行优化。
本发明配网电压无功协同优化的混合整数二阶锥规划***,输出功率获取模块100获取所述主动配电网中光伏电源输出功率,目标函数获取模块200获取所述主动配电网的电压无功多目标协同优化的目标函数,约束条件确定模块300确定主动配电网安全运行需满足的约束条件,模型构建模块400构建主动配电网的电压无功多目标协同优化模型,模型转换模块500将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,求解模块600求解所述混合整数二阶锥规划模型,获得决策变量的运行策略,优化模块700进行协同优化。整个过程中,基于光伏电源输出功率,考虑光伏电源出力的不确定性对配电网电压分布和无功潮流的影响机理,并且基于混合整数二阶锥规划的请求结果对配网电压无功协同优化的混合整数二阶锥规划,能够确保优化的最优性和准确度,确保电力电网的安全运行。
在其中一个实施例中,所述输出功率获取模块100包括:
槽温获取单元,用于获取光伏电源所在地预设周期内各时段的相关环境数据信息,计算光伏电源中光伏发电板的槽温。
工作电流计算单元,用于获取光伏电源接入点的短路电流值,并根据所述光伏发电板的槽温和所述短路电流值,计算光伏发电板在当前环境条件下的工作电流。
工作电压计算单元,用于获取光伏电源接入点的开路电压值,并根据所述光伏发电板的槽温和所述开路电压值,计算光伏发电板在当前环境条件下的工作电压。
填充因素计算单元,用于获取光伏发电板的最大功率点电流和最大功率点电压,根据光伏电源接入点的短路电流值、光伏电源接入点的开路电压值以及光伏发电板的最大功率点电流和最大功率点电压,计算光伏电源填充因数。
输出功率计算单元,用于根据光伏发电板在当前环境条件下的工作电流、光伏发电板在当前环境条件下的工作电压以及光伏电源填充因数,获取光伏电源输出功率。
在其中一个实施例中,所述模型转换模块500还用于引入中间变量,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,将搜索空间限制在有限的凸锥范围内,获得混合整数二阶锥规划模型。
在其中一个实施例中,所述求解模块600还用于通过分支定界-原对偶内点法,求解所述混合整数二阶锥规划模型,获得决策变量的运行策略。
在其中一个实施例中,所述主动配电网的电压无功多目标协同优化的目标函数具体为:
F=min(αf1+βf2)
式中,F为所述主动配电网的电压无功多目标协同优化的目标函数,f1为所述主动配电网的有功网损,f2为所述主动配电网中电压偏移最优值程度,α和β分别为预设权重系数,且α+β=1。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种配网电压无功协同优化的混合整数二阶锥规划方法,应用于包含光伏电源的主动配电网,其特征在于,包括步骤:
获取所述主动配电网中光伏电源输出功率;
获取所述主动配电网的电压无功多目标协同优化的目标函数;
获取主动配电网运行数据,确定主动配电网安全运行需满足的约束条件;
根据所述光伏电源输出功率、所述目标函数以及所述约束条件,构建主动配电网的电压无功多目标协同优化模型;
引入中间变量,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,获得混合整数二阶锥规划模型;
求解所述混合整数二阶锥规划模型,获得决策变量的运行策略;
根据所述决策变量的运行策略,对主动配电网电压无功协同进行优化;
所述主动配电网运行数据包括节点上的负荷功率最大值和最小值S Di,分布式电源接入功率的最大值和最小值S DGi,无功补偿电容器组投切容量的最大值和最小值Q Ci,储能设备吸放功率的最大值和最小值S DSi,外用用电设备充放电功率的最大值和最小值S EVi,电压幅值的最大值和最小值V i,所述约束条件为SDi+SDSi+SEVi+LLi=SDGi+QCi+SBi,
2.根据权利要求1所述的配网电压无功协同优化的混合整数二阶锥规划方法,其特征在于,所述获取所述主动配电网中光伏电源输出功率的步骤包括:
获取光伏电源所在地预设周期内各时段的相关环境数据信息,计算光伏电源中光伏发电板的槽温;
获取光伏电源接入点的短路电流值,并根据所述光伏发电板的槽温和所述短路电流值,计算光伏发电板在当前环境条件下的工作电流;
获取光伏电源接入点的开路电压值,并根据所述光伏发电板的槽温和所述开路电压值,计算光伏发电板在当前环境条件下的工作电压;
获取光伏发电板的最大功率点电流和最大功率点电压,根据光伏电源接入点的短路电流值、光伏电源接入点的开路电压值以及光伏发电板的最大功率点电流和最大功率点电压,计算光伏电源填充因数;
根据光伏发电板在当前环境条件下的工作电流、光伏发电板在当前环境条件下的工作电压以及光伏电源填充因数,获取光伏电源输出功率。
3.根据权利要求1或2所述的配网电压无功协同优化的混合整数二阶锥规划方法,其特征在于,所述引入中间变量,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,获得混合整数二阶锥规划模型的步骤包括:
引入中间变量,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,将搜索空间限制在有限的凸锥范围内,获得混合整数二阶锥规划模型。
4.根据权利要求1或2所述的配网电压无功协同优化的混合整数二阶锥规划方法,其特征在于,所述求解所述混合整数二阶锥规划模型,获得决策变量的运行策略的步骤包括:
通过分支定界-原对偶内点法,求解所述混合整数二阶锥规划模型,获得决策变量的运行策略。
5.根据权利要求1或2所述的配网电压无功协同优化的混合整数二阶锥规划方法,其特征在于,所述主动配电网的电压无功多目标协同优化的目标函数具体为:
F=min(αf1+βf2)
式中,F为所述主动配电网的电压无功多目标协同优化的目标函数,f1为所述主动配电网的有功网损,f2为所述主动配电网中电压偏移最优值程度,α和β分别为预设权重系数,且α+β=1。
6.一种配网电压无功协同优化的混合整数二阶锥规划***,应用于包含光伏电源的主动配电网,其特征在于,包括:
输出功率获取模块,用于获取所述主动配电网中光伏电源输出功率;
目标函数获取模块,用于获取所述主动配电网的电压无功多目标协同优化的目标函数;
约束条件确定模块,用于获取主动配电网运行数据,确定主动配电网安全运行需满足的约束条件,所述主动配电网运行数据包括节点上的负荷功率最大值和最小值S Di,分布式电源接入功率的最大值和最小值S DGi,无功补偿电容器组投切容量的最大值和最小值Q Ci,储能设备吸放功率的最大值和最小值S DSi,外用用电设备充放电功率的最大值和最小值S EVi,电压幅值的最大值和最小值V i,所述约束条件为SDi+SDSi+SEVi+LLi=SDGi+QCi+SBi,
模型构建模块,用于根据所述光伏电源输出功率、所述目标函数以及所述约束条件,构建主动配电网的电压无功多目标协同优化模型;
模型转换模块,用于引入中间变量,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,获得混合整数二阶锥规划模型;
求解模块,用于求解所述混合整数二阶锥规划模型,获得决策变量的运行策略;
优化模块,用于根据所述决策变量的运行策略,对主动配电网电压无功协同进行优化。
7.根据权利要求6所述的配网电压无功协同优化的混合整数二阶锥规划***,其特征在于,所述输出功率获取模块包括:
槽温获取单元,用于获取光伏电源所在地预设周期内各时段的相关环境数据信息,计算光伏电源中光伏发电板的槽温;
工作电流计算单元,用于获取光伏电源接入点的短路电流值,并根据所述光伏发电板的槽温和所述短路电流值,计算光伏发电板在当前环境条件下的工作电流;
工作电压计算单元,用于获取光伏电源接入点的开路电压值,并根据所述光伏发电板的槽温和所述开路电压值,计算光伏发电板在当前环境条件下的工作电压;
填充因素计算单元,用于获取光伏发电板的最大功率点电流和最大功率点电压,根据光伏电源接入点的短路电流值、光伏电源接入点的开路电压值以及光伏发电板的最大功率点电流和最大功率点电压,计算光伏电源填充因数;
输出功率计算单元,用于根据光伏发电板在当前环境条件下的工作电流、光伏发电板在当前环境条件下的工作电压以及光伏电源填充因数,获取光伏电源输出功率。
8.根据权利要求6或7所述的配网电压无功协同优化的混合整数二阶锥规划***,其特征在于,所述模型转换模块还用于引入中间变量,将构建的主动配电网的电压无功多目标协同优化模型进行锥转换,将搜索空间限制在有限的凸锥范围内,获得混合整数二阶锥规划模型。
9.根据权利要求6或7所述的配网电压无功协同优化的混合整数二阶锥规划***,其特征在于,所述求解模块还用于通过分支定界-原对偶内点法,求解所述混合整数二阶锥规划模型,获得决策变量的运行策略。
10.根据权利要求6或7所述的配网电压无功协同优化的混合整数二阶锥规划***,其特征在于,所述主动配电网的电压无功多目标协同优化的目标函数具体为:
F=min(αf1+βf2)
式中,F为所述主动配电网的电压无功多目标协同优化的目标函数,f1为所述主动配电网的有功网损,f2为所述主动配电网中电压偏移最优值程度,α和β分别为预设权重系数,且α+β=1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710217565.9A CN106921164B (zh) | 2017-04-05 | 2017-04-05 | 配网电压无功协同优化的混合整数二阶锥规划方法与*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710217565.9A CN106921164B (zh) | 2017-04-05 | 2017-04-05 | 配网电压无功协同优化的混合整数二阶锥规划方法与*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106921164A CN106921164A (zh) | 2017-07-04 |
CN106921164B true CN106921164B (zh) | 2019-10-22 |
Family
ID=59567144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710217565.9A Active CN106921164B (zh) | 2017-04-05 | 2017-04-05 | 配网电压无功协同优化的混合整数二阶锥规划方法与*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106921164B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108023364B (zh) * | 2017-11-24 | 2019-07-26 | 天津大学 | 基于凸差规划的配电网分布式电源最大接入能力计算方法 |
CN108039711B (zh) * | 2017-12-22 | 2020-06-12 | 国网冀北电力有限公司秦皇岛供电公司 | 一种电压无功协调优化控制方法及装置 |
CN108306342B (zh) * | 2018-01-19 | 2021-04-02 | 广东电网有限责任公司佛山供电局 | 含dg离散功率因数和设备动作限制的配网运行控制方法 |
CN108599142B (zh) * | 2018-03-21 | 2021-03-12 | 广东电网有限责任公司电力科学研究院 | 一种基于权重的混合整数二次规划多目标调度方法 |
CN108985561B (zh) * | 2018-06-08 | 2021-09-07 | 天津大学 | 一种基于机会约束的有源配电网孤岛划分方法 |
CN109245113A (zh) * | 2018-10-19 | 2019-01-18 | 东北大学 | 分布式电源和电动汽车大规模接入的配电网无功电压优化方法 |
CN109638874A (zh) * | 2018-10-24 | 2019-04-16 | 中国电力科学研究院有限公司 | 一种分布式光伏集群控制方法和装置 |
CN111130121A (zh) * | 2019-11-14 | 2020-05-08 | 广东电网有限责任公司 | 一种dg和ev环境下配电网无功补偿***的模糊协调控制计算方法 |
CN111900732B (zh) * | 2020-08-04 | 2022-08-26 | 四川大学 | 考虑风电不确定性及无功潮流与电压约束的风电集群接入方式规划方法 |
CN112052543B (zh) * | 2020-08-11 | 2022-12-23 | 广西大学 | 一种基于混合整数二阶锥规划的保底网架搜索建模方法 |
CN113346504B (zh) * | 2021-06-17 | 2022-06-28 | 国家电网有限公司 | 基于数据知识驱动的主动配电网电压控制方法 |
CN115051392A (zh) * | 2021-12-24 | 2022-09-13 | 国网浙江省电力有限公司衢州供电公司 | 一种含新能源、储能及柔性负荷的无功协调优化控制方法 |
CN116632837B (zh) * | 2023-07-19 | 2023-12-26 | 国网江西省电力有限公司电力科学研究院 | 一种有源配电网的电压调节方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537258A (zh) * | 2015-01-12 | 2015-04-22 | 天津大学 | 分布式储能参与有源配电网运行调节的锥优化建模方法 |
CN104701867A (zh) * | 2015-03-27 | 2015-06-10 | 河海大学 | 基于分支定界法以及原对偶内点法的日前无功优化方法 |
CN105740973A (zh) * | 2016-01-25 | 2016-07-06 | 天津大学 | 基于混合整数锥规划的智能配电网综合电压无功优化方法 |
CN106292743A (zh) * | 2016-10-13 | 2017-01-04 | 哈尔滨工业大学 | 太阳能双轴自动跟踪***及跟踪方法 |
-
2017
- 2017-04-05 CN CN201710217565.9A patent/CN106921164B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537258A (zh) * | 2015-01-12 | 2015-04-22 | 天津大学 | 分布式储能参与有源配电网运行调节的锥优化建模方法 |
CN104701867A (zh) * | 2015-03-27 | 2015-06-10 | 河海大学 | 基于分支定界法以及原对偶内点法的日前无功优化方法 |
CN105740973A (zh) * | 2016-01-25 | 2016-07-06 | 天津大学 | 基于混合整数锥规划的智能配电网综合电压无功优化方法 |
CN106292743A (zh) * | 2016-10-13 | 2017-01-04 | 哈尔滨工业大学 | 太阳能双轴自动跟踪***及跟踪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106921164A (zh) | 2017-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106921164B (zh) | 配网电压无功协同优化的混合整数二阶锥规划方法与*** | |
CN107069814B (zh) | 配网分布式电源容量布点的模糊机会约束规划方法与*** | |
CN105811407B (zh) | 一种基于分布式牛顿法的微电网一次调频控制方法 | |
CN104361403B (zh) | 一种分布式电源与微电网的优化分组配置方法 | |
CN106549392B (zh) | 一种配电网协调控制方法 | |
CN104734196B (zh) | 一种可并网运行的风光储一体微电网的监控装置 | |
Zhang et al. | Research on bi-level optimized operation strategy of microgrid cluster based on IABC algorithm | |
CN104734194B (zh) | 一种可并网运行的风光储一体微电网 | |
CN107025520A (zh) | 配电网新能源消纳能力确定的双层二阶锥规划方法与*** | |
CN106058914B (zh) | 基于Elman算法的配网发电预测技术的电压优化方法 | |
CN110912177A (zh) | 一种多端柔性直流输电***多目标优化设计方法 | |
CN115545290A (zh) | 一种含光伏发电配电网中的分布式储能经济优化配置方法 | |
CN108629499A (zh) | 一种基于二阶锥理论的配电网光伏电站选址定容方法 | |
CN111884263A (zh) | 多目标配电网优化方法 | |
Halder | A maximum power point tracker (MPPT) using the incremental conductance (INC) technique | |
CN109120018A (zh) | 基于一致性迭代算法的混合配电网分布式控制方法及*** | |
Li et al. | Optimal dispatching method for smoothing power fluctuations of the wind-photovoltaic-battery hybrid generation system | |
CN114243692A (zh) | 一种源-网-荷协调优化调度方法 | |
CN110994632A (zh) | 一种基于机会约束规划考虑电压和环保指标的分布式电源布点定容优化的计算方法 | |
CN107910891B (zh) | 一种分布式光伏集群电压双层优化下垂控制方法 | |
CN106780119A (zh) | 基于多主动管理的配电网可再生能源消纳方法 | |
Jadhav | Optimal power flow in wind farm microgrid using dynamic programming | |
CN110380419A (zh) | 一种混合随机和区间变量的不确定潮流计算方法 | |
CN115549216A (zh) | 一种风光储场站有功-无功协调控制方法及*** | |
CN108718093A (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 |