CN107221929B - 基于市场利益的能效电厂优化配置与厂网协调规划方法 - Google Patents
基于市场利益的能效电厂优化配置与厂网协调规划方法 Download PDFInfo
- Publication number
- CN107221929B CN107221929B CN201710388796.6A CN201710388796A CN107221929B CN 107221929 B CN107221929 B CN 107221929B CN 201710388796 A CN201710388796 A CN 201710388796A CN 107221929 B CN107221929 B CN 107221929B
- Authority
- CN
- China
- Prior art keywords
- power plant
- stage
- cost
- power
- plant
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000008901 benefit Effects 0.000 title claims abstract description 19
- 230000000903 blocking effect Effects 0.000 claims abstract description 54
- 230000002146 bilateral effect Effects 0.000 claims abstract description 33
- 230000035945 sensitivity Effects 0.000 claims abstract description 18
- 238000012795 verification Methods 0.000 claims abstract description 5
- 230000005540 biological transmission Effects 0.000 claims description 34
- 239000011159 matrix material Substances 0.000 claims description 21
- 239000013598 vector Substances 0.000 claims description 21
- 238000010276 construction Methods 0.000 claims description 17
- 238000012423 maintenance Methods 0.000 claims description 17
- 230000002068 genetic effect Effects 0.000 claims description 10
- 238000007726 management method Methods 0.000 claims description 5
- 238000010248 power generation Methods 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 2
- HTIQEAQVCYTUBX-UHFFFAOYSA-N amlodipine Chemical compound CCOC(=O)C1=C(COCCN)NC(C)=C(C(=O)OC)C1C1=CC=CC=C1Cl HTIQEAQVCYTUBX-UHFFFAOYSA-N 0.000 claims description 2
- 101100353051 Neurospora crassa (strain ATCC 24698 / 74-OR23-1A / CBS 708.71 / DSM 1257 / FGSC 987) epp-1 gene Proteins 0.000 description 8
- 230000003068 static effect Effects 0.000 description 8
- 239000000243 solution Substances 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000002301 combined effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010977 unit operation 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
-
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- 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/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/10—Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier
-
- 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)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Power Engineering (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
本发明涉及一种基于市场利益的能效电厂优化配置与厂网协调规划方法,包括以下步骤:1)根据灵敏度计算双边交易分摊的阻塞成本;2)对基于市场利益的含能效电厂的厂网协调规划构建三层规划模型;3)采对三层规划模型进行求解,获取电厂优化配置与厂网协调规划最优方案,并对方案进行连通性校验,确保最优方案不存在孤岛网架。与现有技术相比,本发明具有有效地缓解阻塞、实现阻塞成本的快速分配、方法简单、实用性强等优点。
Description
技术领域
本发明涉及电力市场规划领域,尤其是涉及一种基于市场利益的基于市场利益的能效电厂优化配置与厂网协调规划方法。
背景技术
在电力系统中,通常把超过系统安全性限制的状态称为输电阻塞,这些安全性限制包括系统正常运行状态和非正常运行状态下的各种约束,如输电线路或变压器有功潮流热容量限制,节点电压限制,系统稳定限制等。为消除阻塞,应采取一系列的技术、经济措施,使系统运行于一定的安全和可靠性水平内,维持电力市场的有效运行,这一过程称为阻塞管理。
在垄断的体制下,系统调度、各发电厂和输电网同属一个经济实体,阻塞管理只需由调度员向各个电厂发出调度指令,根据需要调节机组出力,同时采取其它调度手段,如断开过负荷线路、调节变压器抽头或移相器、投入无功补偿设备、停止供应可中断负荷等,直至电力系统运行在安全性限制之内。但是在电力市场环境下,输电网和发电厂是独立的各具经济利益的实体,负责系统运行的独立系统调度员,必须平等对待各个市场成员,在社会效益最大化或购电成本最小化等目标下寻求合理的阻塞管理方法,并准确计算阻塞成本,为网络规划提供经济信号,促进系统在长期内稳定、健康地发展同时应将阻塞成本分摊至引起阻塞的市场参与者,为系统短期运行提供激励。
随着输电网络的开放,双边交易的参与者及交易电量迅速增多。在这种趋势下,如何有效地缓解阻塞并实现阻塞成本的快速分配成为电力市场研究的重要内容之一。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种有效地缓解阻塞、实现阻塞成本的快速分配、方法简单、实用性强的基于市场利益的能效电厂优化配置与厂网协调规划方法。
本发明的目的可以通过以下技术方案来实现:
一种基于市场利益的能效电厂优化配置与厂网协调规划方法,包括以下步骤:
1)根据灵敏度计算双边交易分摊的阻塞成本;
2)对基于市场利益的含能效电厂的厂网协调规划构建三层规划模型;
3)采对三层规划模型进行求解,获取电厂优化配置与厂网协调规划最优方案,并对方案进行连通性校验,确保最优方案不存在孤岛网架。
所述的步骤1)中,双边交易分摊的阻塞成本F的表达式为:
F=CT·AT·diag(B)/(ATB)·diag(C)·D·diag(E)
其中,diag(·)表示以矩阵元素为对角元的对角矩阵,A为系统阻塞成本对阻塞线路功率的灵敏度,B为阻塞线路过负荷量,C为以阻塞线路功率倒数为元素的矩阵,D为阻塞线路功率对双边交易功率的灵敏度,E为实际的双边交易功率,CT为全网阻塞成本。
所述的步骤2)具体包括以下步骤:
21)以阻塞成本和投资成本之和最小作为目标函数构建上层规划模型;
22)以投资成本最小为目标函数建立中层规划模型,投资成本包括输电线路投资成本、常规电厂投资成本、EPP投资成本以及切负荷惩罚费用;
23)以电厂的投资成本最小为目标函数建立下层规划模型,电厂投资成本包括新建常规电厂建设成本、新建能效电厂建设成本及运行维修成本。
所述的步骤21)中,上层规划模型的目标函数为:
min{F1+CT,1,F2+CT,2,...,Fh+CT,h,...,FH+CT,H}
其中,F1、F2…Fh...FH分别为第1、2…h…H种情形的各双边交易的分摊结果,CT,1、CT,2…CT,h、CT,H分别为第1、2…h…H种情形的全网阻塞成本,H为基于下层电源规划结果的新增电源节点配置情形种类总数,T为阶段总数,Y为第t阶段包含的年数,为第i个发电机增投标报价,为第i个发电机第t阶段第y年第h种情形下的增发出力,为第j个发电机减投标报价,为第i个发电机第t阶段第y年第h种情形下的减发出力,Th为年最大运行小时数;
所述的步骤21)中,上层规划模型的约束条件为:
其中,为第t阶段第y年第h种情形下的发电功率向量,为第t阶段第y年第h种情形下的负荷功率向量,为第t阶段第y年第h种情形下的系统节点导纳矩阵,为第t阶段第y年第h种情形下的电压相角向量,分别为第i个发电机增投标功率的下限和上限,为发电机减投标功率的下限和上限,分别为第t阶段发电机的功率下限向量和上限向量,为第t阶段第y年第h种情形下线路ij的传输功率,为第t阶段第h种情形下线路ij的最大传输功率。
所述的步骤22)中,中层规划模型的目标函数为:
其中,Fh为第h种情形的中层目标值,r为贴现率,1/(1+r)(t-1)Y为资金的折算系数,Ω为节点的集合,为第t阶段节点i到j的输电线路单位投资费用,Lij为节点i到j之间输电线路的长度,为第h种情形下第t阶段节点i到j之间新建输电线路回数,Gt为电厂的投资成本,为第t阶段新建常规电厂投资成本,为第t阶段新建能效电厂投资成本,为第t阶段常规电厂的运行维修费用,为第h种情形下第t阶段N、N-1网络安全下的切负荷惩罚费用之和,为第h种情形下第t阶段正常运行状态下节点i的切负荷量,为第h种情形下第t阶段第s条线路断开的N-1运行状态下节点i的切负荷量,a、b为正常状态下的切负荷惩罚系数和N-1状态下的切负荷惩罚系数,I为节点数,S为线路回路数,Y为第t阶段包含的年数。
所述的步骤22)中,中层规划模型的约束条件包括正常运行状态下潮流约束、线路有功潮流约束、N-1运行状态下潮流约束、线路有功潮流约束以及输电线路约束。
所述的步骤23)中,下层规划模型的目标函数为:
其中,Gt为电厂的投资成本,为第t阶段常规电厂的建设成本,为第t阶段能效电厂的建设成本,为第t阶段电厂的运行维修费用,M为常规电厂机组类型数,CG,m为新建常规发电机组类型为m的单位容量投资成本,PG,m为m型常规机组的容量,为第t阶段新建m型常规机组的台数,K为能效电厂机组类型数,CE,k为新建能效电厂类型为k的单位容量投资成本,PE,k为k型能效电厂的容量,为第t阶段新建k型能效电厂的台数,N为第t阶段包含所选择的待建常规电厂的常规机组数目,Y为第t阶段包含的年数,为第t阶段第n个常规电厂机组的容量,IOM&E,n为第n个常规电厂机组的单位运行维修费用,为第t阶段第n个常规电厂机组第y年的运行小时数,L为第t阶段包含所选择的待建能效电厂的能效电厂数目,为第t阶段第l个能效电厂机组的容量,IOM,l为第l个能效电厂的单位运行费用,为第t阶段第l个能效电厂机组第y年的运行小时数。
所述的下层规划模型的约束条件包括:
机组备用容量约束、电量约束、系统中能效电厂的容量约束、常规电厂运行时间约束、能效电厂运行时间约束、新建常规电厂机组台数约束、新建能效电厂机组台数约束。
所述的步骤3)中,采用自适应遗传算法对上层决策变量进行求解,并且采用原始-对偶内点法对中层和下层进行快速求解。
与现有技术相比,本发明具有以下优点:
本发明中充分考虑在随着输电网络的开放,双边交易的参与者及交易电量迅速增多的趋势下,有效地缓解阻塞并实现阻塞成本的快速分配的重要性,本发明基于市场利益以能效电厂优化配置与厂网协调规划问题为出发点,通过建立模型和有效算法计算,得出最优规划方案具有方法简单、实用性强的优点。
附图说明
图1为本发明的发明流程图。
图2为本发明混合算法求解厂网协调规划流程图。
图3为实施例的IEEE-RTS 24系统网络结构图。
图4市场利益下不考虑EPP的厂网协调规划阶段一规划方案。
图5市场利益下不考虑EPP的厂网协调规划阶段二规划方案。
图6市场利益下不考虑EPP的厂网协调规划阶段三规划方案。
图7市场利益下考虑EPP的厂网协调规划阶段一规划方案。
图8市场利益下考虑EPP的厂网协调规划阶段二规划方案。
图9市场利益下考虑EPP的厂网协调规划阶段三规划方案。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
如图1所示,本发明提供一种基于市场利益的基于市场利益的能效电厂优化配置与厂网协调规划方法,包括以下步骤:
S1根据灵敏度计算双边交易分摊的阻塞成本,步骤S1中根据灵敏度计算双边交易分摊的阻塞成本,其具体步骤为:
步骤S11:生成以矩阵形式表达的双边交易分摊的阻塞成本,其中各阻塞线路的阻塞成本Cij由下式确定:
式中,CT为全网阻塞成本;分数项表示线路i-j上单位功率所引起的阻塞成本,ΔPij为线路i-j的过负荷量。
从m点到n点的双边交易在阻塞线路i-j中分摊的阻塞成本为:
式中,Pij为线路i-j的功率,Pij,mn为阻塞线路i-j中由m点到n点的双边交易引起的功率,这一项可以通过下式来计算:
式中,表示在m点注入单位功率从n点流出单位功率时线路i-j中的功率变化量;Tmn为双边交易实际交易功率。
从m点到n点的双边交易在全部阻塞线路中所分摊的阻塞成本为:
变形后可得:
式中,表示线路i-j单位功率所引起的阻塞成本;表示单位双边交易在线路i-j中引起的功率变化量;二者之积表示单位交易功率所引起的阻塞成本。ΔPij/Pij是一个系数,表示阻塞成本应由线路中过负荷的部分承担,而不是由全部的功率承担。
假设系统在某时段存在T个双边交易,L条线路发生阻塞。各线路的首末节点分别用i、j表示,各双边交易的功率注入点和接收点分别用m、n表示。那么,对于这个时段产生的阻塞成本,各双边交易的分摊结果可用矩阵表示为:
F=AT·diag(B)·diag(C)·D·diag(E)
式中各矩阵分别为:
其中,diag(·)表示以矩阵元素为对角元的对角矩阵。以上各矩阵中,矩阵A为系统阻塞成本对阻塞线路功率的灵敏度,B为阻塞线路过负荷量;C为以阻塞线路功率倒数为元素的矩阵;D为阻塞线路功率对双边交易功率的灵敏度;E为实际的双边交易功率;F为双边交易分摊的阻塞成本。其中,矩阵E可以根据计算条件确定,B、C可以根据潮流计算获得。
步骤S12:计算灵敏度矩阵,其中阻塞线路最终分摊的成本为:
式中,Cij为各阻塞线路的阻塞成本;C′ij为线路l最终分摊的阻塞成本。用矩阵表示的各双边交易的分摊结果可修改为:
F=CT·AT·diag(B)/(ATB)·diag(C)·D·diag(E)
矩阵D为线路功率对双边交易功率的灵敏度,可以借助于线路功率对节点注入功率的灵敏度来获得。将双边交易分解成两个子交易:从功率注入点到参考节点的交易与从参考节点到功率接收点的交易。这两个交易的合成作用与双边交易的作用相当。线路功率对双边交易功率的灵敏度就等于线路功率对这两个子交易所对应的功率注入点的灵敏度之差。
线路i-j对双边交易的灵敏度Kij,mn为:
式中,Kij,m、Kij,n分别为线路i-j对m、n两点功率的灵敏度;
S2采用三层规划方法对基于市场利益的含能效电厂的厂网协调规划进行建模,具体步骤为:
步骤S21:以阻塞成本(即发电机再调度成本)和投资成本之和最小为目标函数建立上层规划模型,上层目标函数如下:
S=min{F1+CT,1,F2+CT,2,...,Fh+CT,h,...,FH+CT,H}
阻塞成本函数表示如下:
约束条件表示如下:
约束条件包括潮流约束、发电机的增、减投标容量约束、发电机出力约束、线路输电容量约束。
以上步骤S21式中,H为基于下层电源规划结果的H种新增电源节点配置情形;CT,h为第h种情形的全网阻塞成本;Th为年最大运行小时数;分别为第i个发电机增投标报价及第t阶段第y年第h种情形下的增发出力; 分别为第j个发电机减投标报价及第t阶段第y年第h种情形下的减发出力;分别为发电机增投标功率的下限和上限;分别为发电机减投标功率的下限和上限;为第t阶段第y年第h种情形下的发电功率向量;为第t阶段第y年第h种情形下的负荷功率向量;为第t阶段第y年第h种情形下的系统节点导纳矩阵;为第t阶段第y年第h种情形下的电压相角向量;分别为第t阶段发电机的功率下限向量和上限向量;为第t阶段第y年第h种情形下线路ij的传输功率;为第t阶段第h种情形下线路ij的最大传输功率。
步骤S22:以投资成本最小为目标函数建立中层规划模型,投资成本包括输电线路投资成本、电厂(常规电厂和EPP)投资成本以及切负荷惩罚费用,中层目标函数如下:
电源规划的投资运行模型如下:
N、N-1网络安全下的切负荷惩罚费用之和表示如下:
正常运行状态下潮流约束、线路有功潮流约束表示如下:
N-1运行状态下潮流约束、线路有功潮流约束表示如下:
输电线路约束条件表示如下:
以上步骤S22各式中,I为节点数;S为线路回路数;T为规划周期包含的阶段数;Fh为第h种情形的中层目标值;r为贴现率(%),1/(1+r)(t-1)Y为资金的折算系数;Y为第t阶段包含的年数;Ω为节点的集合;为第t阶段节点i到j的输电线路单位投资费用(万元/km);Lij为节点i到j之间输电线路的长度;为第h种情形下第t阶段节点i到j之间新建输电线路回数;为第t阶段新建常规电厂投资成本;为第t阶段新建能效电厂投资成本;为第t阶段常规电厂的运行维修费用;为第h种情形下第t阶段N、N-1网络安全下的切负荷惩罚费用之和;为第h种情形下第t阶段正常运行状态下节点i的切负荷量,为其列向量;为第h种情形下第t阶段第s条线路断开的N-1运行状态下节点i的切负荷量;a、b分别为正常状态下的切负荷惩罚系数和N-1状态下的切负荷惩罚系数;为第h种情形下第t阶段系统节点导纳矩阵;为第h种情形下第t阶段节点相角列向量;为第h种情形下第t阶段常规电厂节点出力列向量, 分别为第t阶段常规电厂节点出力最小值和最大值的列向量;为第h种情形下第t阶段能效电厂节点出力列向量;为第h种情形下第t阶段节点负荷列向量;为第h种情形下第t阶段节点i到节点j之间原有的线路回路数;xij为节点i到j之间单条输电线路的电抗;为第h种情形下第t阶段节点i到节点j之间输电线路上的总潮流;为第t阶段节点i到节点j之间单条回路容量上限;为第t阶段节点i到j之间新建输电线路回数上限。公式中上标有“∧”符号的均代表线路N-1条件下的网络参数及对应的潮流。
步骤S23:以电厂的投资成本最小为目标函数建立下层规划模型,电厂投资成本包括新建常规电厂建设成本、新建能效电厂建设成本及运行维修成本,下层目标函数表示如下:
第t阶段常规电厂的建设成本表示如下:
第t阶段能效电厂的建设成本表示如下:
第t阶段电厂的运行维修费用表示如下,包括常规电厂的燃料成本、环境成本和维修费用以及能效电厂的运行维修费用。
各约束条件表示如下,依次为机组备用容量约束、电量约束、系统中能效电厂的容量约束、常规电厂运行时间约束、能效电厂运行时间约束、新建常规电厂机组台数约束、新建能效电厂机组台数约束。
以上步骤S23各式中M为常规电厂机组类型数;K为能效电厂机组类型数;N为第t阶段包含所选择的待建常规电厂的常规机组数目;L为第t阶段包含所选择的待建能效电厂的能效电厂数目;CG,m为新建常规发电机组类型为m的单位容量投资成本;PG,m为m型常规机组的容量;为第t阶段新建m型常规机组的台数;CE,k为新建能效电厂类型为k的单位容量投资成本;PE,k为k型能效电厂的容量;为第t阶段新建k型能效电厂的台数;为第t阶段第n个常规电厂机组的容量;IOM&E,n为第n个常规电厂机组的单位运行维修费用,包括常规电厂的燃料成本、环境成本及维修费用;为第t阶段第n个常规电厂机组第y年的运行小时数;为第t阶段第l个能效电厂机组的容量;IOM,l为第l个能效电厂的单位运行费用;为第t阶段第l个能效电厂机组第y年的运行小时数;Th为年最大运行小时数;R为机组的备用容量系数;为第t阶段的最大负荷;α为能效电厂与负荷的最小比值;β为能效电厂与负荷的最大比值;为第t阶段第y年的负荷;为第t阶段新建m型常规机组的最大台数;为第t阶段新建k型能效电厂的最大台数。
步骤S24:对所生成的方案进行连通性校验,以确保求得的最优方案不存在孤岛网架。对存在孤岛的网架,通过随机选择孤岛与网架中节点的待架线路来消除孤岛,实现网架连通;对存在独立小网的不连通的网架,通过随机选择独立小网中节点与网架节点间的待架线路来消除独立小网,实现网架连通,最终保证了最优规划方案的网架连通性要求;
S3采用自适应遗传算法和原始-对偶内点法的混合算法对三层规划模型进行求解,具体步骤为:
步骤S31:采用自适应遗传算法对上层决策变量进行求解,并进行全局寻优;
步骤S32:采用原始-对偶内点法对中层和下层进行快速求解,并反馈至上层;
基于阻塞成本和投资成本最小的混合算法求解厂网协调规划流程图如图2所示,具体步骤如下:
(1)输入节点负荷数据、已有机组数据、待建机组数据、待建EPP数据等原始参数,令阶段计数器t=1;
(2)计算该阶段末的最大负荷和每一年末的负荷;
(3)随机生成初始种群,遗传代数计数器gen=1;
(4)判断个体是否满足备用容量约束和EPP容量约束,若满足,计算常规电厂和EPP的建设费用,并以各个机组运行小时数为变量,使用内点法计算机组的运行维修费用,将建设费用和运行维修费用的总和设为下层目标函数值,若不满足,将该个体下层目标函数值置0;
(5)计算下层初始种群的目标函数值和适应度值,根据适应度值进行选择、交叉和变异操作,产生新个体;
(6)将子代作为新的原始种群,并计算子代的目标函数值和适应度值;
(7)判断遗传代数gen是否达到最大值,若达到最大值,则输出最优解及新建机组的节点配置集合,若没有达到最大值,遗传代数计数器gen=gen+1,返回步骤(4);
(8)给节点配置集合编号,取节点配置情况h=1,进入中层电网规划阶段;
(9)随机生成初始种群,遗传代数计数器gen=1;
(10)对个体进行连通性校验及修正,计算线路建设费用,并用内点法计算切负荷惩罚费用;
(11)计算上层初始种群的目标函数值和适应度值,根据适应度值进行选择、交叉和变异操作,产生新个体;
(12)将子代作为新的原始种群,并计算子代的目标函数值和适应度值;
(13)判断遗传代数gen是否达到最大值,若达到最大值,则输出最优解并存储,若没有达到最大值,遗传代数计数器gen=gen+1,返回步骤(10);
(14)判断h是否达到最大值,若达到最大值,则按照中层规划目标函数值由小到大进行排序,输出前m个解及对应规划方案,若没有达到最大值,h=h+1,返回步骤(9);
(15)给规划方案集合编号,取规划方案n=1,进入上层规划阶段;
(16)利用内点法计算最优潮流,并计算出发电机增发出力或减发出力;
(17)根据发电机的增投标或减投标报价计算方案n的阻塞成本及上层目标函数值;
(18)判断n是否达到最大值m,若达到最大值,则输出最优解及最优规划方案,若没有达到最大值,n=n+1,返回步骤(16);
(19)判断t是否达到最大值,若达到最大值,则结束运算,若没有达到最大值,t=t+1,返回步骤(2)。
实施例网络结构图如图3所示,为IEEE-RTS系统,该系统现有24个节点,基础总负荷为5700MW,图中实线为现有输电线路,虚线为待选输电线路。本算例中,将规划周期分为三个阶段,每个阶段为5年,负荷按每年6%的速度增长。输电线路单位投资成本均取200万元/千米,年最大运行小时数取5500小时;正常运行状态下的切负荷惩罚系数α和N-1运行状态下的切负荷惩罚系数β均取100万元/MW,中心参数σ取0.1,对偶间隙取10-7,贴现率为0.08,在matlab2014a平台上进行编程实现。
分析市场利益下不考虑能效电厂的厂网协调规划和考虑能效电厂的厂网协调规划两种情况。两种情况下各阶段厂网协调规划情况如表1、表2所示。
表1市场利益下不考虑EPP的厂网协调规划数据
表2市场利益下考虑EPP的厂网协调规划数据
由表1可知,不考虑EPP的厂网协调规划第一阶段节点3新增300MW和600MW燃煤机组各1台,节点6新增300MW燃煤机组1台,节点18新增600MW燃煤机组2台,新增线路76条,静态投资成本为570.403亿元,动态投资成本为570.403亿元,SO2排放量为57.05万吨,NOX排放量为19.99万吨,CO2排放量为8078.15万吨;第二阶段节点13新增300MW和600MW燃煤机组各1台,节点15新增1000MW燃煤机组1台,节点19新增600MW燃煤机组2台,新增线路35条,静态投资成本为664.3867亿元,动态投资成本为452.1704亿元,SO2排放量为92.66万吨,NOX排放量为31.91万吨,CO2排放量为12134.97万吨;第三阶段节点6新增600MW燃煤机组1台,节点7新增390MW和220MW燃气机组各1台,节点13新增300MW燃煤机组1台,节点15新增300MW燃煤机组2台,节点19新增1000MW燃煤机组1台,节点20新增1000MW核电机组1台,新增线路15条,静态投资成本为868.7347亿元,动态投资成本为402.3923亿元,SO2排放量为129.97万吨,NOX排放量为44.74万吨,CO2排放量为17013.27万吨。三个阶段的静态投资总成本为2103.5244亿元,动态投资总成本为1424.9657亿元,SO2排放总量为279.68万吨,NOX排放总量为96.64万吨,CO2排放总量为37226.39万吨。规划结果如图4、图5、图6所示。
由表2可知,考虑EPP的厂网协调规划第一阶段节点7新增60MW的EPP1,节点9新增90MW的EPP1,节点13新增70MW的EPP1和70MW的EPP3,节点14新增50MW的EPP1,节点15新增80MW的EPP1和80MW的EPP3,节点18新增90MW的EPP1;节点6新增600MW燃煤机组各1台,节点15新增600MW燃煤机组1台,节点18新增600MW燃煤机组1台,新增线路78条,静态投资成本为507.6063亿元,动态投资成本为507.6063亿元,SO2排放量为52.15万吨,NOX排放量为18.29万吨,CO2排放量为7277.07万吨;第二阶段节点2新增45MW的EPP3,节点3新增100MW的EPP4,节点6新增70MW的EPP1,节点8新增80MW的EPP4,节点18新增90MW的EPP4;节点13新增600MW燃煤机组各1台,节点15新增300MW燃煤机组2台,节点19新增600MW和1000MW燃煤机组各1台,新增线路54条,静态投资成本为632.1017亿元,动态投资成本为430.1978亿元,SO2排放量为85.51万吨,NOX排放量为29.47万吨,CO2排放量为11213.05万吨;第三阶段节点1新增50MW的EPP2,节点5新增50MW的EPP2,节点10新增50MW的EPP1和50MW的EPP5,节点16新增50MW的EPP3;节点3新增600MW燃煤机组1台,节点6新增600MW燃煤机组1台,节点7新增220MW燃气机组各2台,节点13新增300MW燃煤机组2台,节点18新增1000MW燃煤机组1台,节点20新增1000MW核电机组1台,新增线路14条,静态投资成本为822.2234亿元,动态投资成本为380.8485亿元,SO2排放量为126.16万吨,NOX排放量为43.33万吨,CO2排放量为16376.83万吨。三个阶段的静态投资总成本为1961.9314亿元,动态投资总成本为1318.6526亿元,SO2排放总量为263.82万吨,NOX排放总量为91.09万吨,CO2排放总量为34866.95万吨。规划结果如图7、图8、图9所示。
Claims (8)
1.一种基于市场利益的能效电厂优化配置与厂网协调规划方法,其特征在于,包括以下步骤:
1)根据灵敏度计算双边交易分摊的阻塞成本,双边交易分摊的阻塞成本F的表达式为:
F=CT·AT·diag(B)/(ATB)·diag(C)·D·diag(E)
其中,diag(·)表示以矩阵元素为对角元的对角矩阵,A为系统阻塞成本对阻塞线路功率的灵敏度,B为阻塞线路过负荷量,C为以阻塞线路功率倒数为元素的矩阵,D为阻塞线路功率对双边交易功率的灵敏度,E为实际的双边交易功率,CT为全网阻塞成本;
2)对基于市场利益的含能效电厂的厂网协调规划构建三层规划模型,具体包括以下步骤:
21)以阻塞成本和投资成本之和最小作为目标函数构建上层规划模型;
22)以投资成本最小为目标函数建立中层规划模型,投资成本包括输电线路投资成本、电厂的投资成本以及切负荷惩罚费用;
23)以电厂的投资成本最小为目标函数建立下层规划模型,电厂的投资成本包括新建常规电厂建设成本、新建能效电厂建设成本及运行维修成本;
3)采对三层规划模型进行求解,获取电厂优化配置与厂网协调规划最优方案,并对方案进行连通性校验,确保最优方案不存在孤岛网架。
2.根据权利要求1所述的一种基于市场利益的能效电厂优化配置与厂网协调规划方法,其特征在于,所述的步骤21)中,上层规划模型的目标函数为:
min{F1+CT,1,F2+CT,2,...,Fh+CT,h,...,FH+CT,H}
其中,F1、F2…Fh...FH分别为第1、2…h…H种情形的各双边交易的分摊结果,CT,1、CT,2…CT,h…CT,H分别为第1、2…h…H种情形的全网阻塞成本,H为基于下层规划结果的新增电源节点配置情形种类总数,T为阶段总数,Y为第t阶段包含的年数,为第i个发电机增投标报价,为第i个发电机第t阶段第y年第h种情形下的增发出力,为第j个发电机减投标报价,为第i个发电机第t阶段第y年第h种情形下的减发出力,Th为年最大运行小时数。
3.根据权利要求2所述的一种基于市场利益的能效电厂优化配置与厂网协调规划方法,其特征在于,所述的步骤21)中,上层规划模型的约束条件为:
其中,为第t阶段第y年第h种情形下的发电功率向量,为第t阶段第y年第h种情形下的负荷功率向量,为第t阶段第y年第h种情形下的系统节点导纳矩阵,为第t阶段第y年第h种情形下的电压相角向量,分别为第i个发电机增投标功率的下限和上限,为发电机减投标功率的下限和上限,分别为第t阶段发电机的功率下限向量和上限向量,为第t阶段第y年第h种情形下线路ij的传输功率,为第t阶段第h种情形下线路ij的最大传输功率。
4.根据权利要求1所述的一种基于市场利益的能效电厂优化配置与厂网协调规划方法,其特征在于,所述的步骤22)中,中层规划模型的目标函数为:
min
其中,Fh为第h种情形的中层目标值,r为贴现率,1/(1+r)(t-1)Y为资金的折算系数,Ω为节点的集合,为第t阶段节点i到j的输电线路单位投资费用,Lij为节点i到j之间输电线路的长度,为第h种情形下第t阶段节点i到j之间新建输电线路回数,Gt为电厂的投资成本,为第t阶段新建常规电厂建设成本,为第t阶段新建能效电厂建设成本,为第t阶段常规电厂和能效电厂的运行维修费用,为第h种情形下第t阶段N、N-1网络安全下的切负荷惩罚费用之和,为第h种情形下第t阶段正常运行状态下节点i的切负荷量,为第h种情形下第t阶段第s条线路断开的N-1运行状态下节点i的切负荷量,a、b为正常状态下的切负荷惩罚系数和N-1状态下的切负荷惩罚系数,I为节点数,S为线路回路数,Y为第t阶段包含的年数。
5.根据权利要求4所述的一种基于市场利益的能效电厂优化配置与厂网协调规划方法,其特征在于,所述的步骤22)中,中层规划模型的约束条件包括正常运行状态下潮流约束、线路有功潮流约束、N-1运行状态下潮流约束、线路有功潮流约束以及输电线路约束。
6.根据权利要求1所述的一种基于市场利益的能效电厂优化配置与厂网协调规划方法,其特征在于,所述的步骤23)中,下层规划模型的目标函数为:
min
其中,Gt为电厂的投资成本,为第t阶段新建常规电厂建设成本,为第t阶段新建能效电厂建设成本,为第t阶段常规电厂和能效电厂的运行维修费用,M为常规电厂机组类型数,CG,m为新建常规发电机组类型为m的单位容量投资成本,PG,m为m型常规机组的容量,为第t阶段新建m型常规机组的台数,K为能效电厂机组类型数,CE,k为新建能效电厂类型为k的单位容量投资成本,PE,k为k型能效电厂的容量,为第t阶段新建k型能效电厂的台数,N为第t阶段包含所选择的待建常规电厂的常规机组数目,Y为第t阶段包含的年数,为第t阶段第n个常规电厂机组的容量,IOM&E,n为第n个常规电厂机组的单位运行维修费用,为第t阶段第n个常规电厂机组第y年的运行小时数,L为第t阶段包含所选择的待建能效电厂的能效电厂数目,为第t阶段第l个能效电厂机组的容量,IOM,l为第l个能效电厂的单位运行费用,为第t阶段第l个能效电厂机组第y年的运行小时数。
7.根据权利要求6所述的一种基于市场利益的能效电厂优化配置与厂网协调规划方法,其特征在于,所述的下层规划模型的约束条件包括:
机组备用容量约束、电量约束、系统中能效电厂的容量约束、常规电厂运行时间约束、能效电厂运行时间约束、新建常规电厂机组台数约束、新建能效电厂机组台数约束。
8.根据权利要求1所述的一种基于市场利益的能效电厂优化配置与厂网协调规划方法,其特征在于,所述的步骤3)中,采用自适应遗传算法对上层决策变量进行求解,并且采用原始-对偶内点法对中层和下层进行快速求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710388796.6A CN107221929B (zh) | 2017-05-27 | 2017-05-27 | 基于市场利益的能效电厂优化配置与厂网协调规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710388796.6A CN107221929B (zh) | 2017-05-27 | 2017-05-27 | 基于市场利益的能效电厂优化配置与厂网协调规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107221929A CN107221929A (zh) | 2017-09-29 |
CN107221929B true CN107221929B (zh) | 2019-12-27 |
Family
ID=59948186
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710388796.6A Active CN107221929B (zh) | 2017-05-27 | 2017-05-27 | 基于市场利益的能效电厂优化配置与厂网协调规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107221929B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108446796A (zh) * | 2018-03-01 | 2018-08-24 | 国网福建省电力有限公司 | 考虑电动汽车负荷需求响应的网-源-荷协调规划方法 |
CN108985502B (zh) * | 2018-07-02 | 2022-02-15 | 国网江西省电力有限公司经济技术研究院 | 基于两层电能替代规划模型的规划方法及算法框架 |
CN110569555B (zh) * | 2019-08-14 | 2023-04-18 | 上海电力设计院有限公司 | 基于配电网过渡方案三层规划模型的计算机辅助规划方法 |
CN111969658B (zh) * | 2020-08-26 | 2021-12-28 | 重庆大学 | 考虑风电的发输电系统防御性-常规协调规划方法 |
CN113191827B (zh) * | 2021-05-28 | 2024-08-13 | 中国电力科学研究院有限公司 | 市场化负荷的阻塞成本确定方法、系统、设备及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104158188A (zh) * | 2014-08-29 | 2014-11-19 | 上海电力学院 | 一种可中断负荷参与下的输电阻塞消除方法 |
CN105160450A (zh) * | 2015-07-08 | 2015-12-16 | 上海电力学院 | 一种考虑环境因素的含能效电厂电源规划方法 |
CN106159944A (zh) * | 2016-08-04 | 2016-11-23 | 上海电力学院 | 低碳环境下基于二层规划模型的多阶段输电网规划方法 |
CN106327007A (zh) * | 2016-08-10 | 2017-01-11 | 上海电力学院 | 一种考虑能效电厂的多阶段输电网二层规划方法 |
-
2017
- 2017-05-27 CN CN201710388796.6A patent/CN107221929B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104158188A (zh) * | 2014-08-29 | 2014-11-19 | 上海电力学院 | 一种可中断负荷参与下的输电阻塞消除方法 |
CN105160450A (zh) * | 2015-07-08 | 2015-12-16 | 上海电力学院 | 一种考虑环境因素的含能效电厂电源规划方法 |
CN106159944A (zh) * | 2016-08-04 | 2016-11-23 | 上海电力学院 | 低碳环境下基于二层规划模型的多阶段输电网规划方法 |
CN106327007A (zh) * | 2016-08-10 | 2017-01-11 | 上海电力学院 | 一种考虑能效电厂的多阶段输电网二层规划方法 |
Non-Patent Citations (2)
Title |
---|
兼顾输电利润和社会成本的输电网二层规划方法;范宏等;《电力系统保护与控制》;20111216;第39卷(第24期);第106-111页 * |
双边交易模式下阻塞成本的快速分配;肖宏飞等;《电网技术》;20050331;第29卷(第6期);第37-42页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107221929A (zh) | 2017-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107221929B (zh) | 基于市场利益的能效电厂优化配置与厂网协调规划方法 | |
CN107392395A (zh) | 一种基于价格激励机制的配电网和微电网协调优化方法 | |
CN113011670B (zh) | 配电网故障应急抢修调度-故障恢复协同方法及装置 | |
CN111079972A (zh) | 一种主动配电网可靠性规划方法、设备及介质 | |
Salkuti et al. | Multi-objective day-ahead real power market clearing with voltage dependent load models | |
CN117937623B (zh) | 配电网智能储能软开关规划方法及装置 | |
CN114358430A (zh) | 一种基于鲁棒优化的新型农村配电网规划策略 | |
CN113364043A (zh) | 一种基于条件风险价值的微电网群优化方法 | |
Lu et al. | A model for balance responsible distribution systems with energy storage to achieve coordinated load shifting and uncertainty mitigation | |
Wang et al. | Two-stage distributionally robust offering and pricing strategy for a price-maker virtual power plant | |
Hinojosa | Static generation capacity expansion planning using linear transmission distribution factors | |
Yin et al. | Distributionally robust transactive control for active distribution systems with SOP-connected multi-microgrids | |
CN111276965B (zh) | 一种基于松弛罚因子的电能量市场优化方法、系统及设备 | |
Jimenez et al. | Long term hydrothermal scheduling linear programming model for large scale power systems | |
Zhu et al. | A clearing method for cascade hydropower spot market considering coupling relations of hydraulic power, electric power and electricity price | |
CN107169607B (zh) | 一种基于成本的能效电厂优化配置与厂网协调规划方法 | |
Bhongade et al. | Genetic algorithm based PID controller design for a multi-area AGC scheme in a restructured power system | |
Rajesh et al. | Congestion management in deregulated environment using Generation Rescheduling with an intelligent approach | |
Nazari et al. | A bilevel scheduling approach for modeling energy transaction of virtual power plants in distribution networks | |
CN111404154A (zh) | 一种配电网供电能力优化方法、设备及存储介质 | |
Siahi et al. | Competitive distribution system planning model integration of DG, interruptible load and voltage regulator devices | |
Xu | Essays on Regional Power System Investment: Value of Planning Model Enhancements, Transmission Generation Storage Co-Optimization, and Border Carbon Adjustment | |
Hermoza et al. | Electricity Market Simulator to Assess the Implementation of Free Prices considering the Transmission Network | |
Boisseleau et al. | A Supply Function Equilibrium Model with Forward Contracts-An Application to Wholesale Electricity Markets | |
Mobarak et al. | Implementation effects of economics and market operations based model for traditionally integrated power systems |
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 |