CN113036820B - 一种光热电站参与的双边调峰辅助服务市场模拟运行方法 - Google Patents
一种光热电站参与的双边调峰辅助服务市场模拟运行方法 Download PDFInfo
- Publication number
- CN113036820B CN113036820B CN202110280732.0A CN202110280732A CN113036820B CN 113036820 B CN113036820 B CN 113036820B CN 202110280732 A CN202110280732 A CN 202110280732A CN 113036820 B CN113036820 B CN 113036820B
- Authority
- CN
- China
- Prior art keywords
- thermal power
- photo
- power station
- auxiliary service
- peak shaving
- 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
- 230000002146 bilateral effect Effects 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004088 simulation Methods 0.000 title claims abstract description 16
- 230000033228 biological regulation Effects 0.000 claims abstract description 68
- 238000005338 heat storage Methods 0.000 claims abstract description 48
- 230000008901 benefit Effects 0.000 claims abstract description 44
- 230000005611 electricity Effects 0.000 claims abstract description 26
- 238000005457 optimization Methods 0.000 claims abstract description 13
- 238000011156 evaluation Methods 0.000 claims abstract description 6
- 238000010248 power generation Methods 0.000 claims description 26
- 230000009467 reduction Effects 0.000 claims description 15
- 230000007246 mechanism Effects 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 230000009194 climbing Effects 0.000 claims description 9
- 239000003245 coal Substances 0.000 claims description 9
- 238000005265 energy consumption Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 5
- 238000010977 unit operation Methods 0.000 claims description 4
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 3
- 239000001569 carbon dioxide Substances 0.000 claims description 3
- 230000017525 heat dissipation Effects 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 6
- 238000010521 absorption reaction Methods 0.000 abstract description 4
- 239000002699 waste material Substances 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000029087 digestion Effects 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- 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/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
- H02J3/0075—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
-
- 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]
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/40—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation wherein a plurality of decentralised, dispersed or local energy generation technologies are operated simultaneously
-
- 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/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种光热电站参与的双边调峰辅助服务市场模拟运行方法,根据光热电站运行状态,建立含储热装置的光热电站的简化运行模型;根据光热电站参与调峰辅助服务市场的报价模型,建立双边调峰辅助服务市场模型,并求解得到出清电量;通过按报价支付的方式进行结算;根据机组申报电价,得到机组实时出清结算电价;根据光热电站、火电、风电和光伏在调峰辅助服务市场协同优化运行时的出清电量和机组实时出清结算电价,进行调峰效益评估,实现光热电站和火电共同参与的双边调峰辅助服务市场的模拟运行。本发明能全面反映光热电站在提供调峰辅助服务中的经济性和发展潜能,有效利用光热电站的调峰价值,进一步提高风电和光伏的消纳率。
Description
技术领域
本发明属于电力系统领域,涉及一种光热电站参与的双边调峰辅助服务市场模拟运行方法。
背景技术
随着新能源装机和发电量的快速增长,弃电问题日益突出,新能源消纳压力日益增加。在传统调峰辅助服务市场下,火电机组被要求提供调峰服务,出让发电空间。然而,火电机组在进行深度调峰时,被迫运行在投油深度调峰状态,利益受损,调峰积极性不高。在这种情况下,火电机组能够带来的调峰资源相比日益增长的新能源带来的调峰需求十分有限,这为提高系统消纳新能源带来了挑战。
随着太阳能光热发电技术(Concentrating Solar Power,CSP)的迅猛发展,光热发电的成本有望进一步降低。CSP通过配置储热装置能够抑制太阳能随机波动对发电功率的影响,保持出力稳定。同时,光热发电机组可对自身出力进行快速调节,最快可达到每分钟20%的装机容量,远高于普通火电机组,从而可以为系统调峰提供一定的爬坡支撑,在调峰辅助服务市场具有巨大潜能。此外,与传统火电机组相比,CSP作为一种利用太阳能的可再生能源发电技术,具有很大的节能减排效益。因此,在光热发电技术发展初期,探索CSP在调峰辅助服务市场的获益方式,对于提高新能源消纳,引导光热发电技术健康发展均具有重大意义。
目前对CSP的研究大多集中在其参与需求侧响应和电网电度运行方面,缺少对CSP参与辅助服务市场的相关研究和模型建立。现有系统内含有光热电站的省份其调峰辅助服务市场没有考虑到光热电站参与调峰辅助服务市场的运作流程和出清机制,光热电站的调峰价值没有得到充分的利用,因而系统的弃风率和弃光率没有得到进一步地降低,且火电作为调峰资源的唯一提供方,经常处于深度调峰状态,经济效益较差。该专利提出一种光热电站和火电厂共同参与的调峰辅助服务市场的模型,可以有效解决光热电站与火电同时参与调峰辅助服务市场的市场运作问题,提高风电和光伏的消纳率,并实现对光热电站调峰价值的有效利用。
发明内容
本发明的目的在于提供一种光热电站参与的双边调峰辅助服务市场模拟运行方法,以填补现有光热电站参与调峰市场机制建模的空白:充分考虑光热电站的运行特性和其他电网技术约束,首先建立含储热的光热电站的运行模型,再建立光热电站在双边调峰辅助服务市场的报价模型,最终完成在综合效益最大的目标函数下光热电站和火电厂共同参与双边调峰辅助服务市场的模型构建。
为了实现上述目的,本发明采用如下技术方案:
一种光热电站参与的双边调峰辅助服务市场模拟运行方法,包括以下步骤:
步骤一:根据光热电站运行状态,建立含储热装置的光热电站的简化运行模型;
步骤二:基于调峰辅助服务市场机制下火电机组的报价模型,建立光热电站参与调峰辅助服务市场的报价模型;
步骤三:根据步骤二的光热电站参与调峰辅助服务市场的报价模型和步骤一光热电站的简化运行模型,建立光热电站、火电、风电、光伏共同参与的双边调峰辅助服务市场模型;
步骤四:通过求解步骤三中的光热电站、火电、风电、光伏共同参与的双边调峰辅助服务市场模型,得到光热电站、火电、风电和光伏在调峰辅助服务市场协同优化运行时的出清电量;
步骤五:通过按报价支付的方式进行结算;根据机组申报电价,得到机组实时出清结算电价;
步骤六:根据步骤四得到的光热电站、火电、风电和光伏在调峰辅助服务市场协同优化运行时的出清电量和步骤五得到的机组实时出清结算电价,进行调峰效益评估,评估市场状况,实现光热电站和火电共同参与的双边调峰辅助服务市场的模拟运行。
本发明进一步的改进在于,步骤一中,光热电站运行状态为充热模式运行、边充热边释放热能发电模式运行或释放热能发电模式运行。
本发明进一步的改进在于,步骤一中,含储热装置的光热电站的简化运行模型包括光热电站的能量流动约束和储热罐系统运行约束;其中,光热电站的能量流动约束如下:
其中,为光热电站c在时刻t收集的热量,ηsf为集热场的光热转换效率,Ssf为集热场的面积,Rt表示时刻t太阳能直接辐射强度,Δt表示时间间隔;其中表示时刻t由集热场传输到导热介质的热量,表示时刻t集热场c弃置的太阳能热量;表示时刻t导热介质传输到汽轮发电机系统用于发电的热量,则分别表示时刻t储热罐储存/释放的热量;Pc,t表示时刻t光热电站的输出功率,ηc表示光热电站的汽轮发电系统热转电的能量转换效率,则分别表示光热电站c的最大、最小出力限制,分别表示光热机组c的上下爬坡速率;
储热罐系统运行约束如下:
|Qc,T-Qc,0|≤δQc,0 (7)
其中,Qc,t/Qc,t+1分别为时刻t/t+1储热罐内的储热量,γc为时刻t的热量耗散比例,ηch、ηdis分别表示储热罐储存/释放热量的过程中的损耗系数;分别表示单位时间储热罐储存/释放热量的最大能量限制,0-1变量xc,t、yc,t分别表示储热罐储存/释放热量的启停状态;表示储热罐为保持正常运行需保留的最小热量,ρc为光热电站在最大出力下能够连续运行的时长;Qc,0、Qc,T分别为一个周期起始和结束时储热罐内热量;δ为波动比例。
本发明进一步的改进在于,步骤二中,调峰辅助服务市场机制下火电机组的报价模型如下式所示:
其中,其中,分别表示火电机组义务调峰、一档调峰、二档调峰的报价,和对应表示一档调峰和二档调峰的报价上限;分别表示火电机组义务调峰、一档调峰和二档调峰申报电量,ΔPi,t表示火电机组i在时刻t各档次总申报电量,Pi f和Pi min分别表示火电机组i在时刻t的日前计划出力、有偿调峰基准、深度调峰基准和最小技术出力。
本发明进一步的改进在于,光热电站参与调峰辅助服务市场的报价模型如下:
本发明进一步的改进在于,步骤三的具体过程为:
首先根据步骤二的光热电站参与调峰辅助服务市场的报价模型,得到报价变量火电一档调峰报价火电二档调峰报价火电一档调峰报量火电二档调峰报量光热电站一档调峰减出力报价光热电站一档调峰减出力报量以及双边调峰辅助服务市场下风电申报的调峰需求电量ΔPw,t和光伏申报的调峰需求电量ΔPs,t和申报的对应价格qw,t、qs,t,再建立光热电站、火电、风电、光伏共同参与的双边调峰辅助服务市场模型,双边调峰辅助服务市场模型的目标函数为调峰辅助服务市场的综合效益。
本发明进一步的改进在于,调峰辅助服务市场的综合效益如下:
Rw,t=qw,tΔPw,t (11)
Rs,t=qs,tΔPs,t (12)
Re,t=(qcoal+qc)ΔPi,t (13)
其中,Rw,t为出力空间出售给风电的收益,Rs,t为出力空间出售给光伏的收益,Re,t为新能源消纳的节能减排效益,qcoal、qc分别表示火电机组单位减发电量带来的煤耗减少成本和二氧化碳减排收益;分别为调用火电机组参与一档调峰和二档调峰的总成本,为调用光热电站参与调峰的成本;
其次,根据步骤一建立的含储热装置的光热电站的简化运行模型和火电机组运行模型建立光热电站、火电、风电、光伏共同参与的双边调峰辅助服务市场模型的约束条件。
其中,表示光热电站c在时刻t的增加的出力,Pt Dw代表风电的总调峰需求,Pt Ds对应表示光伏总调峰需求,αc,t、βc,t分别表示光热电站增加出力和减小出力的0-1变量,则分别表示火电机组i上下爬坡速率,分别表示风电站w和光伏电站s时刻t的日前计划出力,分别表示风电站w和光伏电站s在时刻t的预测最大出力,rw、rs分别表示周期T内弃风率和弃光率要求,表示光热电站c在时刻t的最大技术出力。
本发明进一步的改进在于,光热电站、火电、风电和光伏在调峰辅助服务市场协同优化运行时的出清电量包括光热电站和火电减小的出力空间出清电量ΔPi,t,以及光热电站增加的电量风电和光伏在某时刻t增加的出力空间ΔPw,t、ΔPs,t。
本发明进一步的改进在于,步骤六的具体过程包括以下步骤:
获得风电、光伏、火电和光热电站的上网电价pwind、psolar、pthermal、pCSP以及火电的第一档、第二档调峰成本γI、γII和光热电站的发电成本cCSP,然后再根据步骤四得到的光热电站、火电、风电和光伏在调峰辅助服务市场协同优化运行时的出清电量和机组实时出清结算电价,按照式(25)-(28)分别计算光热电站、火电、风电或光伏的调峰效益;
其中,NCSP、Nthermal、Nw、Ns分别为参与双边调峰辅助服务市场的光热电站、火电、风电和光伏机组数量,RCSP、Rthermal、Rw、Rs分别表示光热电站、火电、风电和光伏参与双边调峰辅助服务市场所获收益,CCSP、Cthermal、Cw、Cs分别表示光热电站、火电、风电和光伏参与双边调峰辅助服务市场的成本,BCSP、Bthermal、Bw、Bs分别表示光热电站、火电、风电和光伏参与双边调峰辅助服务市场的效益。
相对于现有调峰辅助服务模型相比,本发明具有以下有益效果:
本发明从风电和光伏消纳造成的巨大调峰压力出发,考虑火电机组调峰经济性、参与调峰动力不足而光热电站具有巨大调峰潜能的电力系统实际问题,对光热电站提供调峰资源进行建模,基于双边调峰辅助服务市场机制,考虑光热电站和火电同时进行调峰服务报价,建立光热电站参与的双边调峰模型,能全面反映光热电站在提供调峰辅助服务中的经济性和发展潜能,有效利用光热电站的调峰价值,进一步提高风电和光伏的消纳率。相比已有方法,本发明更加全面完整地刻画了光热电站和火电厂共同参与双边调峰辅助服务市场的模型构建问题,并提供了一种对光热电站参与调峰市场的效益评估方法,优势显著。在风电和光伏消纳压力持续增加的背景下,本发明考虑光热电站参与的双边调峰辅助服务市场、节能减排效益的综合效益最大化的目标函数以及光热电站在市场机制下的运行约束下,研究光热-火电-风电-光伏在双边辅助服务市场中协同优化,充分扩展系统调峰资源,提高新能源消纳率,这对促进光热发电技术进一步发展有重大基础研究促进意义,并能填补已有研究尚未关注的空白。
附图说明
图1为光热电站简化运行机理说明图。
图2为改进的IEEE 30节点算例系统说明图。
图3为本发明的整体流程。
图4为算例系统中新能源日前预测出力和DNI预测值。
图5为各类型电源出力功率堆叠图。
图6为光热电站和火电机组和出力情况图。其中,(a)为光热电站,(b)为火电厂光热电站。
具体实施方式
下面结合附图对本发明进行详细描述。
本发明的一种光热电站和火电厂共同参与的双边调峰辅助服务市场模拟运行方法,具体包括以下步骤:
步骤一:总结归纳光热电站参与调峰辅助服务市场时的三种运行状态,建立含储热装置的光热电站的简化运行模型。主要包括:在负荷低谷时刻或者系统发生弃风弃光率超过政策要求的时刻仅充热模式运行,在有光照且系统弃风弃光程度较小的时刻(即弃风弃光率达到政策要求的、但仍存在一定弃风弃光的时刻),边充热边释放热能发电,在其他时刻仅释放热能发电运行。建立含储热装置的光热电站的简化运行模型,如图1所示,该模型主要包括光热电站的能量流动约束和储热罐系统运行约束。其中,光热电站的能量流动约束如下:
其中,约束(1)表示光热电站的能量输入,为光热电站c在时刻t收集的热量,ηsf为集热场的光热转换效率,Ssf为集热场的面积,Rt表示时刻t太阳能直接辐射强度(DirectNormal Irradiance,DNI),Δt表示时间间隔。约束(2)和(3)表示光热电站系统的内部能量流动约束。其中表示时刻t由集热场传输到导热介质的热量,表示时刻t集热场c弃置的太阳能热量。表示时刻t导热介质传输到汽轮发电机系统用于发电的热量,则分别表示时刻t导入/导出储热罐的热量。约束(4)表示光热电站的输出功率限制,其中,Pc,t表示时刻t光热电站的输出功率,ηc表示光热电站的汽轮发电系统热转电的能量转换效率,则分别表示光热电站c的最大、最小出力限制,分别表示光热机组c的上下爬坡速率。
约束(5)-(7)表示储热罐系统的运行约束。
|Qc,T-Qc,0|≤δQc,0 (7)
约束(5)表示储热罐热量满足守恒定律。其中,Qc,t/Qc,t+1分别为时刻t/t+1储热罐内的储热量,γc为时刻t的热量耗散比例,ηch、ηdis分别表示热量导入/导出储热罐的过程中的损耗系数。分别表示单位时间热量导入/导出储热罐的最大能量限制,0-1变量xc,t、yc,t分别表示热量导入/导出储热罐的启停状态。xc,t取值为0时代表无热量导入储热罐,取值为1是表示有热量导入到储热罐,同理,yc,t取值则对应表示热量导出的启停状态。
用约束(6)表示储热罐的储热量Qc,t需在一定范围不能越限。其中,表示储热罐为保持正常运行需保留的最小热量,ρc为光热电站在最大出力下能够连续运行的时长。为了保证光热电站系统的连续稳定运行,一般要求在一个周期结束后储热罐内热量与该周期起始储热罐热量相差不大,用约束(7)表示。其中,Qc,0、Qc,T分别为一个周期起始和结束时储热罐内热量。δ为波动比例,一般取10%-30%。
步骤二:基于现有调峰辅助服务市场机制下火电机组的报价模型,建立光热电站参与调峰辅助服务市场的报价模型。现有调峰辅助服务市场机制下火电机组的报价模型如下式(8)所示,为两档报价模型,还包括义务调峰。
其中,其中,分别表示火电机组义务调峰、一档调峰、二档调峰的报价,和对应表示一档调峰和二档调峰的报价上限。分别表示火电机组义务调峰、一档调峰和二档调峰申报电量,ΔPi,t表示火电机组i在时刻t各档次总申报电量,Pi base、Pi f和Pi min为火电机组运行参数,分别表示为火电机组i在时刻t的日前计划出力(通过日前运行模拟得到)、有偿调峰基准(一般取为火电机组额定出力的50%)、深度调峰基准(一般取为火电机组额定出力的40%)和最小技术出力。
光热电站参与调峰辅助服务市场的报价模型如下:
式(9)分别表示光热电站的报价模型,光热电站按一档报价建立模型,不包括义务调峰。其中,为光热电站在调峰辅助服务市场的减出力报价,表示光热电站在调峰辅助服务市场的报价上限,为光热电站在调峰市场申报减出力电量,为光热电站的日前计划出力。
步骤三:根据步骤二的光热电站参与调峰辅助服务市场的报价模型和步骤一光热电站的简化运行模型,建立光热电站、火电、风电、光伏共同参与的双边调峰辅助服务市场模型。
其中,光热电站和火电为调峰辅助服务的提供方,通过在弃风弃光超过政策要求的时刻减出力出让出力空间;风电、光伏为调峰辅助服务的需求方,通过参与调峰辅助服务市场获得出力空间,增加出力,从而充分利用该时刻的风力、太阳能资源,提高新能源消纳率。
步骤三的具体过程如下:首先根据步骤二火电机组和光热电站报价模型得到的报价变量以及双边调峰辅助服务市场下风电和光伏申报的调峰需求电量ΔPw,t、ΔPs,t和价格qw,t、qs,t,建立该双边调峰辅助服务市场的目标函数。目标函数如式(10)所示,目标函数为调峰辅助服务市场的综合效益,由收益减去成本得到。式(11)-(15)分别表示综合效益的各子项计算公式。
Rw,t=qw,tΔPw,t (11)
Rs,t=qs,tΔPs,t (12)
Re,t=(qcoal+qc)ΔPi,t (13)
其中,收益包括:将出力空间出售给风电和光伏的收益Rw,t、Rs,t,提高新能源消纳的节能减排效益Re,t。qcoal、qc分别表示火电机组单位减发电量带来的煤耗减少成本和二氧化碳减排收益;成本包括调用火电机组参与调峰的成本和调用光热电站参与调峰的成本
其次,根据步骤一建立的含储热装置的光热电站的简化运行模型和已有的火电机组运行模型建立该双边调峰辅助服务市场模型的约束条件。约束(16)-(23)分别描述了系统级电力平衡约束、系统级调峰需求约束、光热电站增减出力约束、光热电站汽轮发电机出力约束、光热电站出力爬坡约束、火电机组出力爬坡约束、风电和光伏参与调峰辅助服务市场增加出力约束和新能源弃风弃光约束。具体地,式(19)和式(20)由式(4)得到:
其中,表示光热电站c在时刻t的增加的出力,Pt Dw代表风电的总调峰需求,计算方法是风电最大预测出力减去日前计划出力之和,Pt Ds对应表示光伏总调峰需求。约束(18)表示增加出力和减小出力不能同时进行。αc,t、βc,t分别表示光热电站增加出力和减小出力的0-1变量。αc,t为1时,表示光热电站增加出力,为0时则表示不增加出力,同理,βc,t为1时,表示光热电站减出力,为0时则表示光热电站不减出力。则分别表示火电机组i上下爬坡速率。分别表示风电站w和光伏电站s时刻t的日前计划出力,分别表示风电站w和光伏电站s在时刻t的预测最大出力。rw、rs分别表示周期T内弃风率和弃光率要求。表示光热电站c在时刻t的最大技术出力。
步骤四:通过求解步骤三中的双边调峰辅助服务市场模型,得到光热电站、火电、风电和光伏在调峰辅助服务市场协同优化运行时的出清电量。
步骤五:根据步骤二机组报价模型得到的光热电站、火电、风电和光伏在调峰辅助服务市场的申报电价,通过按报价支付(Pay as Bid,PCB)的方式进行结算,得到机组实时出清结算电价。
步骤六:根据步骤四得到的光热电站、火电、风电和光伏在调峰辅助服务市场协同优化运行时的出清电量和步骤五得到的机组实时出清结算电价,进行调峰效益评估,评估市场状况,实现光热电站和火电共同参与的双边调峰辅助服务市场的模拟运行。
根据步骤四得到的光热电站、火电、风电和光伏在调峰辅助服务市场协同优化运行时的出清电量ΔPi,t、ΔPw,t、ΔPs,t和步骤五得到的机组实时出清结算电价βc,t、βw,t、βs,t对各类型机组进行效益评估。首先获得风电、光伏、火电和光热电站的上网电价pwind、psolar、pthermal、pCSP以及火电的第一档、第二档调峰成本γI、γII和光热电站的发电成本cCSP。
然后按照式(25)-(28)分别计算光热电站、火电、风电或光伏的调峰效益。
其中,NCSP、Nthermal、Nw、Ns分别为参与双边调峰辅助服务市场的光热电站、火电、风电和光伏机组数量,RCSP、Rthermal、Rw、Rs分别表示光热电站、火电、风电和光伏参与双边调峰辅助服务市场所获收益,CCSP、Cthermal、Cw、Cs分别表示光热电站、火电、风电和光伏参与双边调峰辅助服务市场的成本,BCSP、Bthermal、Bw、Bs分别表示光热电站、火电、风电和光伏参与双边调峰辅助服务市场的效益,由相应收益减去成本得到。
下面以一个修改的IEEE 30算例系统说明本方法实施流程。
该实例以IEEE 30节点系统为背景,用CSP电站替换系统中的2号火电机组,在节点1处接入光伏电站,在节点8处接入风电场。考虑到大规模风电场和光伏电站在节点8的接入,本算例将IEEE-30节点系统中支路10的线路传输容量扩建为原来的5倍,支路40的线路传输容量扩建为原来的2倍,如图2所示。光热电站的运行参数取值如表1所示,调峰辅助服务市场的参数取值如表2所示。
表1光热电站运行参数
表2调峰市场参数
表3新能源消纳情况
表4交易结果
表5效益分析
利用本发明的方法,依据各步骤具体实施,可以得到风电和光伏的消纳情况如表3所示,得到市场交易结果如表4及图5所示,得到火电机组和光热电站参与市场情况如图6中(a)和(b)所示,得到市场主体光热电站、火电、风电和光伏的效益分析如表5所示。
从表3中可以看到,该方法下风电和光伏的消纳程度得到进一步提高。从表5中可以看到,该方法下光热电站、火电、风电和光伏均能获得正的市场效益,实现协作共赢。从图6中可以看到,该方法下光热电站在风电和光伏出力较大的时刻8、12、21-24减小出力,在风电出力较小负荷较大时刻14-15增加出力,有效参与调峰。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案,本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都属于本发明的保护范围。
Claims (10)
1.一种光热电站参与的双边调峰辅助服务市场模拟运行方法,其特征在于,包括以下步骤:
步骤一:根据光热电站运行状态,建立含储热装置的光热电站的简化运行模型;
步骤二:基于调峰辅助服务市场机制下火电机组的报价模型,建立光热电站参与调峰辅助服务市场的报价模型;
步骤三:根据步骤二的光热电站参与调峰辅助服务市场的报价模型和步骤一光热电站的简化运行模型,建立光热电站、火电、风电、光伏共同参与的双边调峰辅助服务市场模型;
步骤四:通过求解步骤三中的光热电站、火电、风电、光伏共同参与的双边调峰辅助服务市场模型,得到光热电站、火电、风电和光伏在调峰辅助服务市场协同优化运行时的出清电量;
步骤五:通过按报价支付的方式进行结算;根据机组申报电价,得到机组实时出清结算电价;
步骤六:根据步骤四得到的光热电站、火电、风电和光伏在调峰辅助服务市场协同优化运行时的出清电量和步骤五得到的机组实时出清结算电价,进行调峰效益评估,评估市场状况,实现光热电站和火电共同参与的双边调峰辅助服务市场的模拟运行。
2.根据权利要求1所述的一种光热电站参与的双边调峰辅助服务市场模拟运行方法,其特征在于,步骤一中,光热电站运行状态为充热模式运行、边充热边释放热能发电模式运行或释放热能发电模式运行。
3.根据权利要求1所述的一种光热电站参与的双边调峰辅助服务市场模拟运行方法,其特征在于,步骤一中,含储热装置的光热电站的简化运行模型包括光热电站的能量流动约束和储热罐系统运行约束;其中,光热电站的能量流动约束如下:
其中,为光热电站c在时刻t收集的热量,ηsf为集热场的光热转换效率,Ssf为集热场的面积,Rt表示时刻t太阳能直接辐射强度,Δt表示时间间隔;其中表示时刻t由集热场传输到导热介质的热量,表示时刻t集热场c弃置的太阳能热量;表示时刻t导热介质传输到汽轮发电机系统用于发电的热量,则分别表示时刻t储热罐储存/释放的热量;Pc,t表示时刻t光热电站的输出功率,ηc表示光热电站的汽轮发电系统热转电的能量转换效率,Pc max、Pc min则分别表示光热电站c的最大、最小出力限制,分别表示光热机组c的上下爬坡速率;
储热罐系统运行约束如下:
|Qc,T-Qc,0|≤δQc,0 (7)
7.根据权利要求6所述的一种光热电站参与的双边调峰辅助服务市场模拟运行方法,其特征在于,调峰辅助服务市场的综合效益如下:
Rw,t=qw,tΔPw,t (11)
Rs,t=qs,tΔPs,t (12)
Re,t=(qcoal+qc)ΔPi,t (13)
其中,Rw,t为出力空间出售给风电的收益,Rs,t为出力空间出售给光伏的收益,Re,t为新能源消纳的节能减排效益,qcoal、qc分别表示火电机组单位减发电量带来的煤耗减少成本和二氧化碳减排收益;分别为调用火电机组参与一档调峰和二档调峰的总成本,为调用光热电站参与调峰的成本;
其次,根据步骤一建立的含储热装置的光热电站的简化运行模型和火电机组运行模型建立光热电站、火电、风电、光伏共同参与的双边调峰辅助服务市场模型的约束条件。
10.根据权利要求1所述的一种光热电站参与的双边调峰辅助服务市场模拟运行方法,其特征在于,步骤六的具体过程包括以下步骤:
获得风电、光伏、火电和光热电站的上网电价pwind、psolar、pthermal、pCSP以及火电的第一档、第二档调峰成本γI、γII和光热电站的发电成本cCSP,然后再根据步骤四得到的光热电站、火电、风电和光伏在调峰辅助服务市场协同优化运行时的出清电量和机组实时出清结算电价,按照式(25)-(28)分别计算光热电站、火电、风电或光伏的调峰效益;
其中,NCSP、Nthermal、Nw、Ns分别为参与双边调峰辅助服务市场的光热电站、火电、风电和光伏机组数量,RCSP、Rthermal、Rw、Rs分别表示光热电站、火电、风电和光伏参与双边调峰辅助服务市场所获收益,CCSP、Cthermal、Cw、Cs分别表示光热电站、火电、风电和光伏参与双边调峰辅助服务市场的成本,BCSP、Bthermal、Bw、Bs分别表示光热电站、火电、风电和光伏参与双边调峰辅助服务市场的效益。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110280732.0A CN113036820B (zh) | 2021-03-16 | 2021-03-16 | 一种光热电站参与的双边调峰辅助服务市场模拟运行方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110280732.0A CN113036820B (zh) | 2021-03-16 | 2021-03-16 | 一种光热电站参与的双边调峰辅助服务市场模拟运行方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113036820A CN113036820A (zh) | 2021-06-25 |
CN113036820B true CN113036820B (zh) | 2023-04-28 |
Family
ID=76470803
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110280732.0A Active CN113036820B (zh) | 2021-03-16 | 2021-03-16 | 一种光热电站参与的双边调峰辅助服务市场模拟运行方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113036820B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114665514A (zh) * | 2021-12-09 | 2022-06-24 | 湖北省电力勘测设计院有限公司 | 储能的光伏消纳与辅助调峰协调控制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107565613A (zh) * | 2017-09-18 | 2018-01-09 | 中国电力工程顾问集团西北电力设计院有限公司 | 一种考虑电力辅助服务市场的光热电站日前优化调度方法 |
CN108256670A (zh) * | 2017-12-22 | 2018-07-06 | 甘肃省电力公司风电技术中心 | 基于热电联产的光热发电及热电机组联合调峰优化模型 |
CN110021932A (zh) * | 2019-04-29 | 2019-07-16 | 西安交通大学 | 双边参与的调峰辅助服务容量选取和经济模型构建方法 |
CN110535185A (zh) * | 2019-08-09 | 2019-12-03 | 东北电力大学 | 一种考虑需求响应的含光热电站多源优化运行方法 |
CN110707756A (zh) * | 2019-10-12 | 2020-01-17 | 华北电力大学 | 一种高比例风电接入电网的光热电站日前调峰优化控制方法 |
CN110930073A (zh) * | 2019-12-14 | 2020-03-27 | 东北电力大学 | 计及价格型需求响应的风-光-光热联合发电系统日前调度方法 |
CN112184000A (zh) * | 2020-09-25 | 2021-01-05 | 华北电力大学 | 含光热电站的新能源发电系统调峰效能评价方法及系统 |
-
2021
- 2021-03-16 CN CN202110280732.0A patent/CN113036820B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107565613A (zh) * | 2017-09-18 | 2018-01-09 | 中国电力工程顾问集团西北电力设计院有限公司 | 一种考虑电力辅助服务市场的光热电站日前优化调度方法 |
CN108256670A (zh) * | 2017-12-22 | 2018-07-06 | 甘肃省电力公司风电技术中心 | 基于热电联产的光热发电及热电机组联合调峰优化模型 |
CN110021932A (zh) * | 2019-04-29 | 2019-07-16 | 西安交通大学 | 双边参与的调峰辅助服务容量选取和经济模型构建方法 |
CN110535185A (zh) * | 2019-08-09 | 2019-12-03 | 东北电力大学 | 一种考虑需求响应的含光热电站多源优化运行方法 |
CN110707756A (zh) * | 2019-10-12 | 2020-01-17 | 华北电力大学 | 一种高比例风电接入电网的光热电站日前调峰优化控制方法 |
CN110930073A (zh) * | 2019-12-14 | 2020-03-27 | 东北电力大学 | 计及价格型需求响应的风-光-光热联合发电系统日前调度方法 |
CN112184000A (zh) * | 2020-09-25 | 2021-01-05 | 华北电力大学 | 含光热电站的新能源发电系统调峰效能评价方法及系统 |
Non-Patent Citations (2)
Title |
---|
批发竞争电力市场双边交易结算模型与报价研究;韩正华;《电气应用》;20080905(第17期);27-32 * |
考虑需求响应与光热电站参与的多源系统优化经济调度;贠韫韵等;《电力系统保护与控制》;20200716(第14期);146-155 * |
Also Published As
Publication number | Publication date |
---|---|
CN113036820A (zh) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115018260B (zh) | 一种梯级水光互补发电系统的调峰优化调度方法 | |
CN110661246A (zh) | 一种城市轨道交通光伏储能系统的容量优化配置方法 | |
CN111799772B (zh) | 一种考虑机组深度调峰的电热系统优化调度方法 | |
CN110163415B (zh) | 一种变工况特性下的多能流系统多目标模糊协同优化方法 | |
CN109993345B (zh) | 一种面向园区的孤岛运行多能互补系统动态经济调度方法 | |
CN110930073B (zh) | 计及价格型需求响应的风-光-光热联合发电系统日前调度方法 | |
CN111030104B (zh) | 一种含储氢装置的多能源系统提升能源利用率的方法 | |
CN105186583B (zh) | 基于多智能体建模的能量路由器及其能量调度方法 | |
CN110994606B (zh) | 一种基于复杂适应系统理论的多能源电源容量配置方法 | |
CN103955862A (zh) | 一种区域电网风电发展的规划方法 | |
CN111030101B (zh) | 一种基于多元化大数据清洁能源消纳联动调控方法及系统 | |
CN112531786A (zh) | 区域电网电力外送方法 | |
CN114021911A (zh) | 含碳捕捉装置的综合能源系统低碳优化调度方法 | |
CN115577929A (zh) | 基于多场景分析的乡村综合能源系统随机优化调度方法 | |
CN110957722B (zh) | 一种含电转气设备的微型能源网日前优化调度方法 | |
CN113036820B (zh) | 一种光热电站参与的双边调峰辅助服务市场模拟运行方法 | |
CN115882460A (zh) | 一种考虑需求侧管理的两阶段新能源微电网优化调度方法 | |
CN116128553B (zh) | 基于绿证和碳交易交互的综合能源调度方法、系统 | |
CN110571868B (zh) | 微电网的优化配置方法 | |
CN112653137A (zh) | 计及碳交易的光热电站与风电系统、低碳调度方法及系统 | |
CN110994698A (zh) | 一种太阳能光伏-光热联合发电系统优化运行方法 | |
CN113988392B (zh) | 一种考虑可靠性需求响应的微网优化规划方法 | |
CN114362152B (zh) | 综合能源系统多时间尺度调度方法 | |
CN112685879B (zh) | 一种区域电热互联能源系统的多目标优化方法 | |
CN114825466A (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 |