CN110929913B - 一种直流跨区互联电网多目标发电计划分解协调计算方法 - Google Patents
一种直流跨区互联电网多目标发电计划分解协调计算方法 Download PDFInfo
- Publication number
- CN110929913B CN110929913B CN201910925232.0A CN201910925232A CN110929913B CN 110929913 B CN110929913 B CN 110929913B CN 201910925232 A CN201910925232 A CN 201910925232A CN 110929913 B CN110929913 B CN 110929913B
- Authority
- CN
- China
- Prior art keywords
- power
- power grid
- target
- constraint
- model
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 40
- 238000010248 power generation Methods 0.000 title claims abstract description 40
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 30
- 230000005540 biological transmission Effects 0.000 claims abstract description 50
- 238000000034 method Methods 0.000 claims abstract description 50
- 238000005457 optimization Methods 0.000 claims abstract description 38
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 14
- 231100000719 pollutant Toxicity 0.000 claims abstract description 14
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 8
- 229910052799 carbon Inorganic materials 0.000 claims abstract description 8
- 230000008878 coupling Effects 0.000 claims description 21
- 238000010168 coupling process Methods 0.000 claims description 21
- 238000005859 coupling reaction Methods 0.000 claims description 21
- 238000012545 processing Methods 0.000 claims description 9
- 230000009977 dual effect Effects 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 6
- 230000003190 augmentative effect Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000012887 quadratic function Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 230000009194 climbing Effects 0.000 claims description 3
- 239000000446 fuel Substances 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 241000238097 Callinectes sapidus Species 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009472 formulation 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
- RJIWZDNTCBHXAL-UHFFFAOYSA-N nitroxoline Chemical compound C1=CN=C2C(O)=CC=C([N+]([O-])=O)C2=C1 RJIWZDNTCBHXAL-UHFFFAOYSA-N 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
Images
Classifications
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- 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
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种直流跨区互联电网多目标发电计划分解协调计算方法,首先以发电成本最小、碳排放最少和污染物排放最少为目标构建了含直流运行工程化约束的跨区互联大电网集中式多目标日前发输电计划模型;然后利用法线边界交叉法将上述多目标优化模型转化为一系列偏好不同的集中式单目标优化子问题;最后针对每个单目标优化子问题构建基于交替方向乘子法的送端和受端电网分解协调模型,实现多目标日前发输电计划的分布式求解。本方法可以实现直流跨区互联大电网多目标日前发输电计划的分布式求解,既能兼顾多个调度目标,又能适用于我国统一调度、分层管理的调度模式。
Description
技术领域
本发明涉及电网优化运行与调度计划领域,尤其涉及一种直流跨区互联电网多目标发电计划分解协调计算方法。
背景技术
资源与需求的反向分布使得我国电网呈现明显的送受端特征,跨区直流线路的建设解决了电能远距离传输的困难,但也使得送受端电网耦合程度加深,调控运行一体化特征明显加强。考虑到送、受端电网之间存在地域、经济、政策等各类因素的影响,通常难以采用一体化方法计算全网发输电计划。常规的做法是由上层电网调制中心制定跨区直流输电计划,各分调在此基础上制定各区域的发电计划,这种模式难以充分挖掘全网资源,调度经济性有待提高。
分解协调计算是解决这一类问题的有效方法之一。通过对耦合约束的处理,可以将原问题分解为多个相互独立的子问题,既保证了各区域发电计划制定的独立性、维护了不同利益主体信息隐私,又保证了调度策略的最优性。
直流跨区互联大电网多目标日前发输电计划的分解协调计算方法的思路是在不丢失解的多样性的前提下,将多目标优化问题转化为一系列单目标优化子问题,继而针对各优化子问题采用拉格朗日松弛处理耦合约束,将原问题分解为送端和受端优化问题,通过设置收敛条件和参数更新策略实现全网发输电计划的分解协调计算。目前针对发输电计划的求解方法大致可分为三类:
第一类方法考虑了直流线路的工程化约束,构建了全网集中式发输电计划模型,并采用直接求解的方法,如文献《促进跨区新能源消纳的直流联络线功率优化模型及分析》(电力系统自动化,2017年41卷18期第152页)所述。但直接求解规模太大、未考虑信息私密性问题,且不适用于我国分层分级的调度现状。
第二类方法将送、受端电网等效为发电机与常规负荷的组合,对等效后的互联电网构建数学模型并直接求解,如文献《以直流联络线运行方式优化提升新能源消纳能力的新模式》(电力系统自动化,2015年39卷3期第36页)所述。但该方法的效果受到等效精度的影响,难以获得精确的最优解。
第三类方法采用分解协调的方法求解集中式日前发输电计划优化问题,如文献《考虑跨区联络线交易计划的多区域互联系统分散调度方法》(电力系统自动化,2018年42卷16期第32页)所述。但该方法可以在保证各分区信息隐私的同时获得全局最优解,但仅能应用于单目标的情况下,无法处理多目标问题。
发明内容
发明目的:为提升直流跨区互联大电网发输电计划的精细化程度,提出一种适用于我国分层分级调度模式的直流跨区互联大电网多目标日前发输电计划的分解协调计算方法。
技术方案:本发明提供一种直流跨区互联电网多目标发电计划分解协调计算方法,包括以下步骤:
(1)计及发电成本最小、碳排放最少和污染物排放最少三个优化目标和直流运行工程化约束条件,构建直流跨区互联大电网一体化多目标日前发输电计划模型;
(2)构建目标函数空间下的乌托邦面并规格化目标函数,采用权重等分的方法均匀求取乌托邦面上的偏好点,进而采用法线边界交叉法将多目标优化问题转化为一系列偏好不同的单目标优化子问题;
(3)对于各单目标优化子问题利用增广拉格朗日乘子法处理耦合约束条件,然后基于交替方向乘子法构建送端和受端电网分解协调模型并独立求解各自的发输电计划,通过上层调度中心协调边界变量和拉格朗日乘子的更新,最终实现全网发输电计划的分布式求解。
进一步地,所述步骤(1)包括以下步骤:
(11)利用二次函数模型描述各优化目标,表达式如下:
发电成本最小:
其中,T为时段数;NS和NWS为送端区域火电机组和风电场的数量;NR和NWR为受端区域火电机组和风电场的数量;PW,i,t为第i个风电场在时段t的计划有功出力;PG,i,t为火电机组i在时段t的有功出力;NIDR为受端电网激励型需求响应(Incentive-based DemandResponse,IDR)的数量;φi为第i个IDR的报价;PIDR,i,t为第i个IDR在时段t的调用量;
碳排放量最少:
其中,bi,2、bi,1和bi,0为机组i的CO2排放系数;
污染物排放最少:
其中,ci,2、ci,1和ci,0为机组i污染物排放系数,上述污染物包括SO2、NOX等火电机组运行过程中排放的污染气体;
(12)直流跨区互联大电网一体化多目标日前发输电计划模型约束条件为:
其中,和为送端和受端电网等效直流线路出力;PDC,i,max和PDC,i,min为直流线路i的最大、最小传输功率;ΔPDC,i,up和ΔPDC,i,down为直流线路i的功率向上调整、向下调整限值;和分别为联络线d在时段t是否爬坡、滑坡的标志,是0/1变量;M为一个大数,用以限制和的取值;σd,t表示联络线d在时段t是否调整功率,为0/1变量;和表示联络线d在时段t功率是否开始变化、是否结束变化,也为0/1变量;Kd表示联络线d每日功率调节次数上限;TN为联络线功率维持不变的最小时段数;为减少联络线计划的波动,在联络线功率结束变化之后需要在一定时间内维持功率不变;QD为跨区交易合同电量;
(13)考虑功率平衡约束、机组出力限值约束、常规机组爬坡约束、常规机组启停约束、旋转备用约束和线路传输功率约束构建了跨区互联大电网集中式多目标日前发输电计划模型,简化表达如下:
其中,xS和xR表示送端电网和受端电网的控制变量;g(x)和h(x)表示模型中的非耦合等式约束和不等式约束;I(x)表示模型中的耦合等式约束,即功率平衡约束中的等效直流功率约束。
进一步地,所述步骤(2)包括以下步骤:
(21)计算乌托邦面端点并形成乌托邦面:
分别以f1(xS,xR)、f2(xS,xR)和f3(xS,xR)为目标函数进行单目标优化可以得到三组对应的最优解,记为x1*、x2*和x3*;这三组最优解对应的三个目标函数值记为fi*(f1(xi*),f2(xi*),f3(xi*)),i=1,2,3;f1*、f2*和f3*即为目标函数空间下的乌托邦面端点,以这三点为顶点确定的三角平面记为乌托邦面;
(22)规格化目标函数:
对于第i个目标函数,规格化公式如下:
(23)乌托邦面内取点:
其中,δij为乌托邦面上的第j个点对第i个目标函数的偏好系数,满足δ1j+δ2j+δ3j=1;
(24)多目标模型转化为单目标模型:
求取乌托邦面的准法线向量,公式为:
过kj作准法线向量的平行线,该平行线与可行域的交点即为该偏好系数下多目标优化问题的解,即满足下式:
其中,d为距离变量,当d取最大时满足过kj的平行线与可行域相交;
通过增加约束式(14)可将式(7)转化为单目标优化问题,转化后的模型如下:
其中,fi,S表示送端电网第i个目标函数,fi,R表示受端电网第i个目标函数,满足fi,S+fi,R=fi,i=1,2,3。
进一步地,所述步骤(3)包括以下步骤:
(31)耦合约束处理:
模型(15)无法直接采用分解协调计算的原因在于后四组约束条件为耦合约束,即后四组约束中既含有送端电网变量,又含有受端电网变量,采用增广拉格朗日松弛将后四组约束条件松弛至目标函数中,新的目标函数为:
其中,λ1、λ2、λ3、λ4为耦合约束对应的拉格朗日乘子;ρ为惩罚系数;I1、I2、I3、I4为后四组耦合约束的简写,计算公式为:
(32)构建分解协调模型:
为使得送端和受端电网的调控中心能够并行计算各自的发输电计划,采用同步型交替方向乘子法将模型分解为送端和受端两个部分,送端电网需要计算的模型为:
受端电网需要计算的模型为:
分解协调计算的变量迭代更新公式为:
(33)分布式求解:
送端电网优化问题和受端电网优化问题都是一个简单的混合整数规划问题,可采用任意方法求解;每次迭代计算后,需要向上级电网调度中心发送边界变量及拉格朗日乘子信息,并判断是否满足如下收敛条件:
其中,r和s分别为原始残差与对偶残差;ε为迭代收敛精度;若第k次迭代满足式(21)收敛条件,则计算完成,输出最终的全网发输电计划结果;反之,转入(32),利用式(20)更新边界变量与拉格朗日乘子值,并进行下一轮迭代。
有益效果:与现有技术相比,本发明的有益效果:1、该方法计及了发电成本最小、碳排放最少和污染物排放最少三个方面的综合利益,实现了电网调度运行的综合最优;2、本发明利用分解协调技术将全局问题的解耦,降低了计算的规模、保护了分区电网隐私。
附图说明
图1为本发明流程图;
图2为直流跨区互联电网的分解协调框架示意图;
图3为仿真系统示意图;
图4为图3系统中的A区域风电预测值图;
图5为图3系统中的B区域风电预测值图;
图6为图3系统中的区域A、B的日负荷曲线图;
图7为图3系统的发输电计划在目标空间的展示图;
图8为本发明迭代原始残差的收敛曲线图;
图9为本发明迭代对偶残差的收敛曲线图;
图10为本发明与直接法计算得到的总发电成本结果对比图;
图11为本发明与加权法计算得到的发输电计划结果对比图。
具体实施方式
下面结合附图进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
图1是直流跨区互联大电网多目标日前发输电计划的分解协调计算方法流程图。如图1所示,本发明提供一种跨区互联大电网多目标日前发输电计划的分解协调计算方法,首先是计及多个调度目标构建集中式多目标日前发输电计划模型,并规格化目标函数,通过在乌托邦面内均匀取点将多目标优化问题转化为一系列单目标问题;然后采用拉格朗日松弛法处理耦合约束,将原问题分解为送端发输电计划问题和受端发输电计划问题,输入送受端电网风电预测值、负荷预测值等基本数据,并初始化迭代次数;最后迭代求解送、受端电网子问题,并更新边界变量与拉格朗日乘子,直到满足收敛条件。
图2展示了直流跨区互联电网的分解协调框架,以全网分析决策中心作为上层协调层,将直流线路的功率等效为两台等效发电机组,送受端电网通过与协调层交换边界信息与协调量实现分布式求解。参考图1和图2,该方法具有以下步骤:
步骤1、计及发电成本、碳排放和污染物排放三个目标和跨区直流工程化运行约束,构建跨区互联大电网多目标日前发输电计划模型。
(1)构建各优化目标的二次函数模型:
1)发电成本最小
式中:式中:T为时段数;NS和NWS为送端区域火电机组和风电场的数量;NR和NWR为受端区域火电机组和风电场的数量;PW,i,t为第i个风电场在时段t的计划有功出力;PG,i,t为火电机组i在时段t的有功出力;NIDR为受端电网激励型需求响应(Incentive-based DemandResponse,IDR)的数量;φi为第i个IDR的报价;PIDR,i,t为第i个IDR在时段t的调用量。Ci(·)为火电机组i在时段t的燃料成本函数,用二次函数表示:
式中:ai,2、ai,1和ai,0为机组i的燃料成本系数;ui,t为0/1变量,ui,t=0和ui,t=1分别表示机组i在时段t处于停机、启动状态;Si为机组i启动成本。
Wi(·)为第i个风电场的弃风惩罚函数:
Wi(PW,i,t)=μi(PWpre,i,t-PW,i,t) (3)
式中:μi为风电场i的弃风惩罚系数,单位为元/(MW·h);PWpre,i,t为风电场i在时段t的预测出力。
2)碳排放量最少
式中:bi,2、bi,1和bi,0为机组i的CO2排放系数。
3)污染物排放最少
式中:ci,2、ci,1和ci,0为机组i污染物排放系数,上述污染物包括SO2、NOX等火电机组运行过程中排放的污染气体。
(2)构建跨区直流运行的工程化约束条件:
式中:和为送端和受端电网等效直流线路出力;PDC,i,max和PDC,i,min为直流线路i的最大、最小传输功率;ΔPDC,i,up和ΔPDC,i,down为直流线路i的功率向上调整、向下调整限值;和分别为联络线d在时段t是否爬坡、滑坡的标志,是0/1变量;M为一个大数,用以限制和的取值;σd,t表示联络线d在时段t是否调整功率,为0/1变量;和表示联络线d在时段t功率是否开始变化、是否结束变化,也为0/1变量;Kd表示联络线d每日功率调节次数上限;TN为联络线功率维持不变的最小时段数;为减少联络线计划的波动,在联络线功率结束变化之后需要在一定时间内维持功率不变;QD为跨区交易合同电量。
(3)构建跨区互联大电网集中式多目标日前发输电计划模型:
在(1)和(2)的基础上,考虑功率平衡约束、机组出力限值约束、常规机组爬坡约束、常规机组启停约束、旋转备用约束和线路传输功率约束构建了跨区互联大电网集中式多目标日前发输电计划模型,简化表达如下:
式中:xS和xR表示送端电网和受端电网的控制变量;g(x)和h(x)表示模型中的非耦合等式约束和不等式约束;I(x)表示模型中的耦合等式约束,即上述述功率平衡约束中的等效直流功率约束。
步骤2、将多目标日前发输电计划模型转化为多个单目标日前发输电计划模型。
(1)计算乌托邦面端点并形成乌托邦面:
分别以f1(xS,xR)、f2(xS,xR)和f3(xS,xR)为目标函数进行单目标优化可以得到三组对应的最优解,记为x1*、x2*和x3*。这三组最优解对应的三个目标函数值记为fi*(f1(xi*),f2(xi*),f3(xi*)),i=1,2,3。f1*、f2*和f3*即为目标函数空间下的乌托邦面端点,以这三点为顶点确定的三角平面记为乌托邦面。
(2)规格化目标函数:
目标函数之间因存在单位、量级等差异而无法直接比较,因此需要对上述三个目标函数进行规格化处理。对于第i个目标函数,规格化公式如下:
(3)乌托邦面内取点:
式中:δij为乌托邦面上的第j个点对第i个目标函数的偏好系数,满足δ1j+δ2j+δ3j=1。采用权重等分的方法计算偏好系数,公式如下:
(4)多目标模型转化为单目标模型:
求取乌托邦面的准法线向量,公式为:
过kj作准法线向量的平行线,该平行线与可行域的交点即为该偏好系数下多目标优化问题的解,即满足下式:
式中:d为距离变量,当d取最大时满足过kj的平行线与可行域相交。
通过增加约束式(14)可将式(7)转化为单目标优化问题,公式如下:
式中:fi,S表示送端电网第i个目标函数,fi,R表示受端电网第i个目标函数,满足fi,S+fi,R=fi,i=1,2,3。
步骤3、构建分解协调计算模型并实现分布式迭代求解。
(1)耦合约束处理:
模型(15)无法直接采用分解协调计算的原因在于后四组约束条件为耦合约束,即后四组约束中既含有送端电网变量,又含有受端电网变量。采用增广拉格朗日松弛将后四组约束条件增广至目标函数中,新的目标函数为
式中:λ1、λ2、λ3、λ4为耦合约束对应的拉格朗日乘子;ρ为惩罚系数;I2、I3、I4为后三组耦合约束的简写,计算公式为:
(2)构建分解协调模型:
为使得送端和受端电网的调控中心能够并行计算各自的发输电计划,采用同步型交替方向乘子法将模型分解为送端和受端两个部分。送端电网需要计算的模型为:
受端电网需要计算的模型为:
分解协调计算的变量迭代更新公式为:
(3)分布式求解:
送端电网优化问题和受端电网优化问题都是一个简单的MIQP问题,可以直接求解。每次迭代计算后,需要向上级电网调度中心发送边界变量及拉格朗日乘子信息,并判断是否满足如下收敛条件:
式中:r和s分别为原始残差与对偶残差;ε为迭代收敛精度。若第k次迭代满足式(21)收敛条件,则计算完成,输出最终的全网发输电计划结果;反之,转入(4.2),利用式(20)更新边界变量与拉格朗日乘子值,并进行下一轮迭代。
图3为仿真系统示意图。如图3所示,作为本发明的示例,仿真系统选择两个IEEE39节点系统。为构造出送受端特征,在A区域(送端)接入4个风电场,风电场的预测曲线如图4所示;在B区域(受端)接入2个风电场,风电场的预测曲线如图5所示。A、B区域的负荷预测曲线如图6所示。利用本发明进行仿真计算,仿真结果如图7至11。
图7展示了本发明计算得到的日前发输电计划在目标函数空间下的展示图。图中粉色圆点表示离散的偏好系数,以0.1为步长,共66个点;蓝色星点表示不同调度偏好对应的多目标日前发输电计划解。粉点与蓝点一一对应,所有蓝色星点构成的三维曲面就是目标空间下日前发输电计划解集,可以看出解集分布均匀,可以较全面的表示出不同偏好系数下的调度计划。
图8、图9展示了偏好系数δ=(0.3,0.5,0.2)时本发明的收敛情况。原始残差与对偶残差均随着的迭代而变下,经过39次迭代后原始残差与对偶残差均下降于收敛精度以下,共耗时269s。
图10展示了本发明与集中式计算结果的对比。可以看出,本文方法可以实现快速收敛。在迭代13次后由本文方法计算得到的总发电成本已经趋近于集中式方法。在经过39次迭代后本文方法计算得到的总发电成本与集中式方法相比,误差仅为0.4459%。表1进一步对比了两种计算方式下三个目标函数计算结果的误差和计算时间的差异,本文方法得到的发电成本、CO2排放量和SO2排放量与集中式方法求得的结果基本一致,最大误差不到0.5%。
表1各目标计算结果的对比
图11展示了本发明与加权法计算结果的对比。可以看出由于权重系数与多目标最优解之间是非线性关系,利用线性加权法得到的帕累托解集大多分布在边缘位置,最优解较为集中且趋于一个平面,难以准确描绘不同偏好下的调度计划。而本发明得到的帕累托解集则更加均匀,可满足不同调度偏好下的要求。
Claims (3)
1.一种直流跨区互联电网多目标发电计划分解协调计算方法,其特征在于,包括以下步骤:
1、计及发电成本最小、碳排放最少和污染物排放最少三个优化目标和直流运行工程化约束条件,构建直流跨区互联大电网一体化多目标日前发输电计划模型;
2、构建目标函数空间下的乌托邦面并规格化目标函数,采用权重等分的方法均匀求取乌托邦面上的偏好点,进而采用法线边界交叉法将多目标优化问题转化为一系列偏好不同的单目标优化子问题;
3、对于各单目标优化子问题利用增广拉格朗日乘子法处理耦合约束条件,然后基于交替方向乘子法构建送端和受端电网分解协调模型并独立求解各自的发输电计划,通过上层调度中心协调边界变量和拉格朗日乘子的更新,最终实现全网发输电计划的分布式求解;
所述步骤1包括以下步骤:
1.1、利用二次函数模型描述各优化目标,表达式如下:
发电成本最小:
其中,T为时段数;Ci(·)为火电机组i在时段t的燃料成本函数;ui,t为0/1变量,ui,t=0和ui,t=1分别表示机组i在时段t处于停机、启动状态;NS和NWS为送端区域火电机组和风电场的数量;NR和NWR为受端区域火电机组和风电场的数量;PW,i,t为第i个风电场在时段t的计划有功出力;PG,i,t为火电机组i在时段t的有功出力;NIDR为受端电网激励型需求响应的数量;φi为第i个IDR的报价;PIDR,i,t为第i个IDR在时段t的调用量;
碳排放量最少:
其中,bi,2、bi,1和bi,0为机组i的CO2排放系数;
污染物排放最少:
其中,ci,2、ci,1和ci,0为机组i污染物排放系数,上述污染物包括SO2、NOX等火电机组运行过程中排放的污染气体;
1.2、直流跨区互联大电网一体化多目标日前发输电计划模型约束条件为:
其中,和为送端和受端电网等效直流线路出力;PDC,i,max和PDC,i,min为直流线路i的最大、最小传输功率;△PDC,i,up和△PDC,i,down为直流线路i的功率向上调整、向下调整限值;和分别为联络线d在时段t是否爬坡、滑坡的标志,是0/1变量;M为一个大数,用以限制和的取值;σd,t表示联络线d在时段t是否调整功率,为0/1变量;和表示联络线d在时段t功率是否开始变化、是否结束变化,也为0/1变量;Kd表示联络线d每日功率调节次数上限;TN为联络线功率维持不变的最小时段数;为减少联络线计划的波动,在联络线功率结束变化之后需要在一定时间内维持功率不变;QD为跨区交易合同电量;
1.3、考虑功率平衡约束、机组出力限值约束、常规机组爬坡约束、常规机组启停约束、旋转备用约束和线路传输功率约束构建了跨区互联大电网集中式多目标日前发输电计划模型,简化表达如下:
其中,xS和xR表示送端电网和受端电网的控制变量;g(x)和h(x)表示模型中的非耦合等式约束和不等式约束;I(x)表示模型中的耦合等式约束,即功率平衡约束中的等效直流功率约束。
2.根据权利要求1所述的一种直流跨区互联电网多目标发电计划分解协调计算方法,其特征在于,所述步骤2包括以下步骤:
2.1、计算乌托邦面端点并形成乌托邦面:
分别以f1(xS,xR)、f2(xS,xR)和f3(xS,xR)为目标函数进行单目标优化可以得到三组对应的最优解,记为x1*、x2*和x3*;这三组最优解对应的三个目标函数值记为fi*(f1(xi*),f2(xi*),f3(xi*)),i=1,2,3;f1*、f2*和f3*即为目标函数空间下的乌托邦面端点,以这三点为顶点确定的三角平面记为乌托邦面;
2.2、规格化目标函数:
对于第i个目标函数,规格化公式如下:
其中,fi min和fi max为第i个目标函数的最小值和最大值,计算公式为:
fi min=fi(xi*)i=1,2,3 式9
fi max=max{fi(x1*),fi(x2*),fi(x3*)}i=1,2,3 式10
2.3、乌托邦面内取点:
其中,δij为乌托邦面上的第j个点对第i个目标函数的偏好系数,满足δ1j+δ2j+δ3j=1;
2.4、多目标模型转化为单目标模型:
求取乌托邦面的准法线向量,公式为:
过kj作准法线向量的平行线,该平行线与可行域的交点即为该偏好系数下多目标优化问题的解,即满足下式:
其中,d为距离变量,当d取最大时满足过kj的平行线与可行域相交;
通过增加约束式14可将式7转化为单目标优化问题,转化后的模型如下:
其中,fi,S表示送端电网第i个目标函数,fi,R表示受端电网第i个目标函数,满足fi,S+fi,R=fi,i=1,2,3。
3.根据权利要求1所述的一种直流跨区互联电网多目标发电计划分解协调计算方法,其特征在于,所述步骤3包括以下步骤:
3.1、耦合约束处理:
式15无法直接采用分解协调计算的原因在于后四组约束条件为耦合约束,即后四组约束中既含有送端电网变量,又含有受端电网变量,采用增广拉格朗日松弛将后四组约束条件松弛至目标函数中,新的目标函数为:
其中,λ1、λ2、λ3、λ4为耦合约束对应的拉格朗日乘子;ρ为惩罚系数;I1、I2、I3、I4为后四组耦合约束的简写,计算公式为:
3.2、构建分解协调模型:
为使得送端和受端电网的调控中心能够并行计算各自的发输电计划,采用同步型交替方向乘子法将模型分解为送端和受端两个部分,送端电网需要计算的模型为:
受端电网需要计算的模型为:
分解协调计算的变量迭代更新公式为:
3.3、分布式求解:
送端电网优化问题和受端电网优化问题都是一个简单的混合整数规划问题,可采用任意方法求解;每次迭代计算后,需要向上级电网调度中心发送边界变量及拉格朗日乘子信息,并判断是否满足如下收敛条件:
其中,r和s分别为原始残差与对偶残差;ε为迭代收敛精度;若第k次迭代满足式21收敛条件,则计算完成,输出最终的全网发输电计划结果;反之,转入步骤3.2,利用式20更新边界变量与拉格朗日乘子值,并进行下一轮迭代。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910925232.0A CN110929913B (zh) | 2019-09-27 | 2019-09-27 | 一种直流跨区互联电网多目标发电计划分解协调计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910925232.0A CN110929913B (zh) | 2019-09-27 | 2019-09-27 | 一种直流跨区互联电网多目标发电计划分解协调计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110929913A CN110929913A (zh) | 2020-03-27 |
CN110929913B true CN110929913B (zh) | 2022-04-19 |
Family
ID=69848877
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910925232.0A Active CN110929913B (zh) | 2019-09-27 | 2019-09-27 | 一种直流跨区互联电网多目标发电计划分解协调计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110929913B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112101730A (zh) * | 2020-08-18 | 2020-12-18 | 华南理工大学 | 考虑通讯传输失误的全分布式电网-区域热网联合调度方法 |
CN112766559B (zh) * | 2021-01-14 | 2023-04-28 | 华南理工大学 | 基于改进nbi法的电气热互联系统多目标优化调度方法 |
CN113162025B (zh) * | 2021-01-27 | 2023-03-24 | 四川大学 | 一种含需求响应的电气互联网络分散式低碳经济调度方法 |
CN113505468B (zh) * | 2021-06-08 | 2022-12-13 | 西安交通大学 | 区间线性交流最大传输容量模型的构建和解析方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015090090A1 (zh) * | 2013-12-18 | 2015-06-25 | 国家电网公司 | 一种特高压交流跨区电能交易中综合网损率的确定方法 |
WO2016170912A1 (ja) * | 2015-04-22 | 2016-10-27 | 株式会社日立製作所 | 電力系統電圧無効電力監視制御装置及び方法 |
CN106786799A (zh) * | 2017-01-03 | 2017-05-31 | 国电南瑞科技股份有限公司 | 一种直流联络线功率阶梯化发电计划优化方法 |
CN107302231A (zh) * | 2017-08-23 | 2017-10-27 | 国网山东省电力公司经济技术研究院 | 一种小水电群接入电网的两目标随机经济调度方法 |
CN107529678A (zh) * | 2017-10-11 | 2018-01-02 | 中国南方电网有限责任公司 | 可扩展安全稳定约束的交直流电网多目标优化控制方法 |
CN107784373A (zh) * | 2016-08-26 | 2018-03-09 | 华北电力大学(保定) | 一种考虑节能减排的跨省购电计划安排方法 |
-
2019
- 2019-09-27 CN CN201910925232.0A patent/CN110929913B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015090090A1 (zh) * | 2013-12-18 | 2015-06-25 | 国家电网公司 | 一种特高压交流跨区电能交易中综合网损率的确定方法 |
WO2016170912A1 (ja) * | 2015-04-22 | 2016-10-27 | 株式会社日立製作所 | 電力系統電圧無効電力監視制御装置及び方法 |
CN107784373A (zh) * | 2016-08-26 | 2018-03-09 | 华北电力大学(保定) | 一种考虑节能减排的跨省购电计划安排方法 |
CN106786799A (zh) * | 2017-01-03 | 2017-05-31 | 国电南瑞科技股份有限公司 | 一种直流联络线功率阶梯化发电计划优化方法 |
CN107302231A (zh) * | 2017-08-23 | 2017-10-27 | 国网山东省电力公司经济技术研究院 | 一种小水电群接入电网的两目标随机经济调度方法 |
CN107529678A (zh) * | 2017-10-11 | 2018-01-02 | 中国南方电网有限责任公司 | 可扩展安全稳定约束的交直流电网多目标优化控制方法 |
Non-Patent Citations (2)
Title |
---|
一种求解交直流互联电网分布式最优潮流的同步ADMM 方法;梁梓均;《电力系统保护与控制》;20181201;第46卷(第23期);第29-36页 * |
大电网多目标动态优化调度的;杨柳青等;《电工技术学报》;20160331;第31卷(第6期);第177-186页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110929913A (zh) | 2020-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110929913B (zh) | 一种直流跨区互联电网多目标发电计划分解协调计算方法 | |
Lin et al. | Decentralized solution for combined heat and power dispatch through benders decomposition | |
Reynolds et al. | Operational supply and demand optimisation of a multi-vector district energy system using artificial neural networks and a genetic algorithm | |
Hemmati | Technical and economic analysis of home energy management system incorporating small-scale wind turbine and battery energy storage system | |
Ruiz-Cortes et al. | Optimal charge/discharge scheduling of batteries in microgrids of prosumers | |
Zhao et al. | Distributed model predictive control strategy for islands multimicrogrids based on noncooperative game | |
Ouammi et al. | Optimal control of power flows and energy local storages in a network of microgrids modeled as a system of systems | |
Mohamed et al. | System modelling and online optimal management of microgrid using multiobjective optimization | |
CN111030188A (zh) | 一种含分布式和储能的分层分级控制策略 | |
CN102509175B (zh) | 分布式供电系统可靠性优化方法 | |
CN111082451B (zh) | 基于场景法的增量配电网多目标优化调度模型 | |
Fan et al. | Optimally coordinated expansion planning of coupled electricity, heat and natural gas infrastructure for multi-energy system | |
CN109390973A (zh) | 一种考虑通道约束的送端电网电源结构优化方法 | |
CN110689286B (zh) | 中长期电力市场下风火打捆电厂的最优合同电量决策方法 | |
CN110535132A (zh) | 一种基于鲁棒优化的电力系统建设规划方法 | |
CN113890023A (zh) | 一种综合能源微网分布式经济调度优化方法及系统 | |
Shotorbani et al. | Enhanced real-time scheduling algorithm for energy management in a renewable-integrated microgrid | |
Saha | Adaptive model-based receding horizon control of interconnected renewable-based power micro-grids for effective control and optimal power exchanges | |
CN114742314A (zh) | 一种考虑相关性的电热综合能源系统分布鲁棒优化方法 | |
Dan et al. | A scenario-based branch-and-bound approach for MES scheduling in urban buildings | |
CN113298407B (zh) | 一种工业园区电-气综合能源系统优化调度模型建立方法 | |
Jabari et al. | Robust self-scheduling of pvs-wind-diesel power generation units in a standalone microgrid under uncertain electricity prices | |
Bhuvaneswari et al. | Online economic environmental optimization of a microgrid using an improved fast evolutionary programming technique | |
CN113869742A (zh) | 基于行动家和评论家网络的综合供需侧的电力调度系统 | |
CN114240203A (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 |