CN115276088A - 一种联合发电系统中光热电站储热容量优化配置方法 - Google Patents

一种联合发电系统中光热电站储热容量优化配置方法 Download PDF

Info

Publication number
CN115276088A
CN115276088A CN202210923499.8A CN202210923499A CN115276088A CN 115276088 A CN115276088 A CN 115276088A CN 202210923499 A CN202210923499 A CN 202210923499A CN 115276088 A CN115276088 A CN 115276088A
Authority
CN
China
Prior art keywords
heat storage
power station
photo
thermal power
optimal
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.)
Granted
Application number
CN202210923499.8A
Other languages
English (en)
Other versions
CN115276088B (zh
Inventor
徐敏
白望望
崔炎
冯智慧
王涛
贾春蓉
张洪源
张耀忠
陈宝奇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Lanzhou University of Technology
Economic and Technological Research Institute of State Grid Gansu Electric Power Co Ltd
Original Assignee
Lanzhou University of Technology
Economic and Technological Research Institute of State Grid Gansu Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Lanzhou University of Technology, Economic and Technological Research Institute of State Grid Gansu Electric Power Co Ltd filed Critical Lanzhou University of Technology
Priority to CN202210923499.8A priority Critical patent/CN115276088B/zh
Publication of CN115276088A publication Critical patent/CN115276088A/zh
Application granted granted Critical
Publication of CN115276088B publication Critical patent/CN115276088B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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
    • G06Q30/00Commerce
    • G06Q30/02Marketing; Price estimation or determination; Fundraising
    • G06Q30/0201Market modelling; Market analysis; Collecting market data
    • G06Q30/0206Price or cost determination based on market factors
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/22The renewable source being solar energy
    • H02J2300/24The renewable source being solar energy of photovoltaic origin

Abstract

本发明公开了一种联合发电系统中光热电站储热容量优化配置方法,该方法包括:针对光伏出力预测和负荷预测的不确定性,采用场景分析法获取出现概率最大的光伏出力预测值和负荷预测值;建立光热电站储热容量两阶段优化配置模型;在第一阶段中,以光热电站综合投资成本最低为优化目标,求解得到最优光热电站镜场面积和储热系统最大储热容量;在第二阶段中,以联合发电系统的经济效益最大和电能供需差波动最小作为优化目标,结合获得的光伏出力预测值和负荷预测值及最优镜场面积与最大储热容量、计算得到每个时段最优的储热量;通过本方法可以最大化地平抑电能供需差的波动,提高整个系统的运行效益。

Description

一种联合发电系统中光热电站储热容量优化配置方法
技术领域
本发明涉及光热电站储热容量配置技术领域,特别涉及一种联合发电系统中光热电站储热容量优化配置方法。
背景技术
当前能源短缺、环境污染已成为全世界公认的头号问题,制约着人类发展的速度。作为解决此类问题的关键,以太阳能为代表的新能源发电得到了迅速的发展。但由于光伏发电过多依赖于天气的状况,造成其出力的随机性与不确定性。光伏发电出力的随机性与不确定性势必限制其大规模并网发电。
相较于光伏发电,光热发电作为新兴的发电技术,因其配置有大容量的储热系统,可以灵活地调节其出力。在光伏电站附近建设光热电站形成光伏-光热联合发电将成为常见的能源开发形式。
为了利用光热电站的出力来缓解光伏出力的不确定性,需要对联合发电系统中的光热电站的储热容量进行合理的配置。
现有的关于光热电站储热容量优化配置,存在着以下几点不足:(1)在光热电站的个体优化中,常以平准化度电成本最低为优化目标,设立固定的太阳倍数来优化储热时间,未将太阳倍数作为决策变量去优化;(2)在联合发电系统中光热电站储热容量优化配置,常以综合成本最低为优化目标去优化储热容量,但这类发明中又未考虑光伏出力、负荷预测的不确定性。
发明内容
本发明的目的在于克服现有技术的缺陷,并为此提供一种联合发电系统中光热电站储热容量优化配置方法,该方法可以在光伏-光热组成的联合发电系统中,利用光热电站的出力来缓解光伏出力的不确定性造成的电能供需差距的波动,针对光伏预测出力、负荷预测的不确定性,本发明采用场景分析法来描述;并建立了光热电站储热容量两阶段优化配置模型,可以计算出每个时段最优的储热量;进而可以最大限度地提高光伏上网电量,同时提高整个系统的运行效益。
为实现上述目的,本发明采取的技术方案为:
本发明提供一种联合发电系统中光热电站储热容量优化配置方法,包括以下步骤:
S1、针对光伏出力预测和负荷预测的不确定性,采用场景分析法获得出现概率最大的光伏出力预测值和负荷预测值;
S2、建立光热电站储热容量两阶段优化配置模型;
S3、基于光热电站储热容量两阶段优化配置模型,在第一阶段中,建立光热电站综合投资成本最低的目标函数,求解得到最优太阳倍数与储热时间,结合光热电站的太阳倍数模型与储热时间模型,进而得到光热电站的最优镜场面积与储热系统的最大储热容量;
S4、在第二阶段中,建立联合发电系统的经济效益最大的目标函数和电能供需差波动最小的目标函数,结合场景分析法获得的光伏出力预测值和负荷预测值以及第一阶段求得的所述最优镜场面积与所述最大储热容量,计算得到每个时段最优的储热量。
进一步地,所述步骤S1中,所述采用场景分析法获取出现概率最大的光伏出力预测值和负荷预测值,具体为:
构建光伏出力预测误差和负荷预测误差的概率密度函数;
根据所述概率密度函数,采用拉丁超立方抽样的方法来生成光伏出力预测及负荷预测的初始场景;
采用k-means算法将生成的所述初始场景进行缩减,获得出现概率最大的光伏出力预测值和负荷预测值。
进一步地,所述步骤S3中,所述光热电站综合投资成本最低的目标函数为:
min f1=FSF+FTES+FPB
式中,f1为光热电站综合投资成本折现值;FSF为集热系统投资成本折现值;FTES为储热系统投资成本折现值;FPB为发电系统投资成本折现值;具体表达式如下:
FSF=CSFSMδCRF(r,n)
FTES=CTEShTESδCRF(r,n)
FPB=CPBPCSP
Figure BDA0003778650720000031
式中,CSF为光热电站太阳倍数投资成本;SM为光热电站太阳倍数;CTES为光热电站储热时间投资成本,hTES为光热电站储热时间;CPB为光热电站发电系统单位容量投资成本;PCSP为光热电站额定容量;δCRF(r,n)为等年值折算系数;r为贴现率;n为光热电站运行年限;
通过求解所述综合投资成本最低的目标函数,确定最优FSF与FTES
根据所述最优FSF与FTES,确定光热电站的最优太阳倍数SM与储热时间hTES
进一步地,所述步骤S3中,所述光热电站的最优镜场面积与储热系统的最大储热容量的获取过程具体为:
将所述最优太阳倍数SM与储热时间hTES分别带入太阳倍数模型与储热时间模型中,求解得到光热电站的最优镜场面积与储热系统的最大储热容量。
进一步地,所述太阳倍数模型为:
Figure BDA0003778650720000032
式中,SM为太阳倍数;SSF为镜场面积;IDNI为标准工况下太阳直接辐射强度;ηPB为发电系统热电转换效率;ηSF为集热系统光热转换效率;
所述储热时间模型为:
Figure BDA0003778650720000033
式中,hTES为储热时间;ηdis为储热系统放热效率;
Figure BDA00037786507200000410
为储热系统最大储热容量。
进一步地,所述步骤S4中,所述联合发电系统的经济效益最大的目标函数为:
max f2=F1+F2-F3-F4
式中,f2表示联合系统运行效益;F1表示联合系统市场售电效益;F2表示联合系统并网环境效益;F3表示光伏弃光惩罚成本;F4表示联合系统运行维护成本;具体表达式如下:
Figure BDA0003778650720000042
Figure BDA0003778650720000043
Figure BDA0003778650720000044
Figure BDA0003778650720000045
式中,T表示调度时段;
Figure BDA0003778650720000046
表示t时刻市场售电的分时电价;Cprice为光热电站市场售电电价;Pt pv、Pt csp分别表示t时刻光伏、光热实际出力;Kpv、Kcsp分别表示光伏、光热并网环境效益系数;
Figure BDA0003778650720000047
表示光伏弃光惩罚成本系数;Cpv、CTES和CPB分别表示光伏、光热电站储热系统和光热电站发电系统运维成本系数。
进一步地,所述步骤S4中,所述电能供需差波动最小的目标函数,由等效负荷方差最小的目标函数来表示,具体为:
Figure BDA0003778650720000048
Plt=PLt-Pt pv-Pt csp
Figure BDA0003778650720000049
式中,f3表示等效负荷方差;Plt表示t时刻等效负荷值;Plt.av表示t时刻等效负荷平均值;PLt表示t时刻负荷预测值;T表示调度时段。
与现有技术相比,本发明具有如下有益效果:
该方法可以在光伏-光热组成的联合发电系统中,便于利用光热电站的出力来缓解光伏出力的不确定性造成的电能供需差距的波动,针对光伏预测出力、负荷预测的不确定性,本发明采用场景分析法来描述;并建立了光热电站储热容量两阶段优化配置模型,第一阶段以光热电站综合投资成本最低为优化目标;可以获取光热电站的最优太阳倍数与储热时间,根据最优太阳倍数与储热时间可以计算出光热电站的最优镜场面积与最大储热容量;第二阶段以联合发电系统的经济效益最大和电能供需差波动最小为优化目标,可以计算出每个时段最优的储热容量;进而可以最大限度地提高光伏上网电量,同时提高整个系统的运行效益。
附图说明
图1为本发明提供的一种联合发电系统中光热电站储热容量优化配置方法的流程图。
图2为本发明实施例提供的光伏-光热联合发电系统的结构图。
图3为本发明实施例提供的光热电站能量流动示意图。
图4为本发明实施例提供的光伏出力预测初始场景示意图。
图5为本发明实施例提供的负荷预测初始场景示意图。
图6为本发明实施例提供的光伏出力预测典型场景示意图。
图7为本发明实施例提供的负荷预测典型场景示意图。
图8为本发明实施例提供的光热电站各时段储热容量与充、放热功率的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参见图1所示,本发明实施例提供了一种联合发电系统中光热电站储热容量优化配置方法,具体方法如下:
S1、针对光伏出力预测和负荷预测的不确定性,采用场景分析法获得出现概率最大的光伏出力预测值和负荷预测值;
S2、建立光热电站储热容量两阶段优化配置模型;
S3、基于光热电站储热容量两阶段优化配置模型,在第一阶段中,建立光热电站综合投资成本最低的目标函数,求解得到最优太阳倍数与储热时间,结合光热电站的太阳倍数模型与储热时间模型,进而得到光热电站的最优镜场面积与储热系统的最大储热容量;
S4、在第二阶段中,建立联合发电系统的经济效益最大的目标函数和电能供需差波动最小的目标函数,结合场景分析法获得的光伏出力预测值和负荷预测值以及第一阶段求得的所述最优镜场面积与所述最大储热容量,计算得到每个时段最优的储热量。
下面分别对上述各个步骤进行详细的说明,本发明实施例整体方法如下:
首先需要说明的是:本实施例中的联合发电系统为光伏-光热联合发电系统,其结构如图2所示,该系统主要由光伏发电系统、光热发电系统组成。其中,光热发电系统主要由集热子系统、储热子系统、发电子系统构成。
光热电站靠传热介质作为能量的传递媒介,进而完成整个运行过程;其能量流动过程如图3所示;
图3中:It DNI为t时刻太阳直接辐射强度;It DNI,curt为t时刻的弃光量;Pt curt为t时刻舍弃的热功率;Pt SF为t时刻从集热系统中获取到的热能;Pt SP为t时刻从集热系统供给发电系统的热能;Pt ST为t时刻从集热系统储存在储热系统中的热能;Pt TP为t时刻由储热系统传递至发电系统中的热能;Pt CSP为t时刻发电系统输出的电功率;ηSF为集热系统光热转换效率;ηPB为发电系统热电转换效率;ηcha为储热系统储热效率;ηdis为储热系统放热效率。
进一步地,本实施例中,针对光伏出力预测和负荷预测的不确定性,采用场景分析法获取出现概率最大的光伏出力预测值和负荷预测值,具体包括:
构建光伏出力预测误差和负荷预测误差的概率密度函数;
进一步地,根据概率密度函数,采用拉丁超立方抽样的方法来生成光伏出力预测及负荷预测的初始场景;采用k-means算法将生成的大量初始场景进行缩减,选取出现概率最大的光伏出力预测值和负荷预测值。
进一步地,光热电站因配有大容量的储热系统,在光伏-光热联合发电系统中,可以平抑联合系统的短期出力波动。因此,建立光热电站储热容量两阶段优化配置模型。
在本实施例中,第一阶段优化配置模型为:
建立光热电站综合投资成本最低的目标函数,表示为:
min f1=FSF+FTES+FPB (1)
式中,f1为光热电站综合投资成本折现值;FSF为集热系统投资成本折现值;FTES为储热系统投资成本折现值;FPB为发电系统投资成本折现值;
具体表达式如下:
FSF=CSFSMδCRF(r,n) (2)
FTES=CTEShTESδCRF(r,n) (3)
FPB=CPBPCSP (4)
Figure BDA0003778650720000071
式中,CSF为光热电站太阳倍数投资成本;SM为光热电站太阳倍数;CTES为光热电站储热时间投资成本,hTES为光热电站储热时间;太阳倍数和储热时间为光热电站规划建造过程中的重要指标;CPB为光热电站发电系统单位容量投资成本;PCSP为光热电站额定容量;δCRF(r,n)为等年值折算系数;r为贴现率;n为光热电站运行年限。
通过求解目标函数min f1=FSF+FTES+FPB,可以确定最优FSF和FTES;根据FSF=CSFSMδCRF(r,n),求得光热电站的最优太阳倍数SM;根据FTES=CTEShTESδCRF(r,n),确定光热电站的最优储热时间hTES
通常采用太阳倍数来描述光场面积,太阳倍数是指在标准工况下(DNI为950W/m2,风速为5m/s,温度为25℃)聚光集热系统产生的热功率与发电系统在额定出力情况下所需热功率之比,一般大于1;
本实施例中,针对集热环节建立的太阳倍数模型为:
Figure BDA0003778650720000081
式中,SM为太阳倍数;SSF为镜场面积;IDNI为标准工况下太阳直接辐射强度;ηPB为发电系统热电转换效率;ηSF为集热系统光热转换效率;PCSP为光热电站额定输出功率。
本实施例中,针对储热环节建立的储热时间模型为:
Figure BDA0003778650720000082
式中,hTES为储热时间;ηdis为储热系统放热效率;
Figure BDA0003778650720000083
、为储热系统最大储热容量。
进一步地,根据上述目标函数求得的最优太阳倍数SM和给定的光热电站额定容量PCSP,代入太阳倍数模型中进一步求得最优镜场面积SSF
进一步地,根据上述目标函数求得的储热时间hTES和给定的光热电站额定容量PCSP,代入储热时间模型求得储热系统的最大储热容量
Figure BDA0003778650720000084
进一步地,第二阶段优化配置模型为:
建立联合发电系统的经济效益最大的目标函数,表示为:
max f2=F1+F2-F3-F4 (8)
式中,f2表示联合系统运行效益;F1表示联合系统市场售电效益;F2表示联合系统并网环境效益;F3表示光伏弃光惩罚成本;F4表示联合系统运行维护成本;具体表达式如下:
Figure BDA0003778650720000091
Figure BDA0003778650720000092
Figure BDA0003778650720000093
Figure BDA0003778650720000094
式中,T表示调度时段;
Figure BDA0003778650720000095
表示t时刻市场售电的分时电价;Cprice为光热电站市场售电电价;Pt pv、Pt csp分别表示t时刻光伏、光热实际出力;Kpv、Kcsp分别表示光伏、光热并网环境效益系数;
Figure BDA0003778650720000096
表示光伏弃光惩罚成本系数;Cpv、CTES和CPB分别表示光伏、光热电站储热系统和光热电站发电系统运维成本系数。
约束条件为:
基于上述图3的光热电站能流图,将传热介质看作成一个节点,忽略在传热介质中的能量损失;
则光热电站内部瞬时热功率平衡约束为:
Pt SF+Pt TP=Pt ST+Pt SP (13)
式中,Pt SF为t时刻从集热系统中获取到的热能;Pt TP为t时刻由储热系统传递至发电系统中的热能;Pt ST为t时刻从集热系统储存在储热系统中的热能;Pt SP为t时刻从集热系统供给发电系统的热能。
光热电站集热环节约束为:
集热系统t时刻将太阳能转化成的热能:
Pt thr solar=ηSFSSFIDNI t (14)
式中:ηSF为光-热转化效率;SSF为镜场面积;IDNI t为t时刻的光照直接辐射指数DNI。
光热电站储热环节中储热系统热能平衡约束为:
Et TES=(1-γTES)Et-1 TES+(Pt cha-Pt dis)Δt (15)
式中,Et TES为t时刻储热系统储热容量;Et-1 TES为t-1时刻储热系统储热容量;Pt cha、Pt dis分别为t时刻光热电站储热与放热功率;γTES为热耗散系数;Δt为时间间隔。
储热系统t时刻储热容量约束:
Figure BDA0003778650720000101
式中,
Figure BDA0003778650720000102
分别为储热系统最大、最小储热容量。
光热电站储热系统在储热与放热的过程中都会伴随着热量的损失,可以通过引入储热效率与放热效率来描述,即:
Pt cha=ηchaPt ST (17)
Figure BDA0003778650720000103
式中,ηcha、ηdis分别为光热电站储热系统的储热效率、放热效率。
光热电站储热系统储热、放热功率不能同时进行且放热状态只有在机组启动时发生:
Figure BDA0003778650720000104
式中,
Figure BDA0003778650720000105
分别为储热系统最大储热、放热功率;
Figure BDA0003778650720000106
分别为储、放热二进制变量,1表示系统储热,0表示系统放热;
Figure BDA0003778650720000107
为光热电站的工作状态变量。
发电系统输入的热功率和输出电功率之间的函数关系约束为:
Pt SP≈f(Pt csp)=Pt cspPB+Ut PBEPB,SU (20)
式中,Pt SP为t时刻从集热系统供给发电系统的热能;Pt CSP为t时刻发电系统输出的电功率;Ut PB为t时刻光热电站的状态;EPB,SU为光热电站开始发电所需的启动能量;
光热电站出力约束为:
Figure BDA0003778650720000111
式中,
Figure BDA0003778650720000112
为光热电站的工作状态变量;
Figure BDA0003778650720000113
分别为光热电站的最大、最小出力;
光热机组最小启、停时间约束为:
Figure BDA0003778650720000114
式中,
Figure BDA0003778650720000115
分别为光热机组最小启、停时间;
光热机组的爬坡速率约束为:
Figure BDA0003778650720000116
式中,
Figure BDA0003778650720000117
分别为光热机组向上、向下爬坡速率。
进一步地,建立光伏-光热联合发电系统的电能供需差波动最小的目标函数,由等效负荷方差最小的目标函数来表示,具体为:
Figure BDA0003778650720000118
Plt=PLt-Pt pv-Pt csp (25)
Figure BDA0003778650720000119
式中,f3表示等效负荷方差;Plt表示t时刻等效负荷值;Plt.av表示t时刻等效负荷平均值;PLt表示t时刻负荷预测值;T表示调度时段。
约束条件为:
功率平衡约束:
Pt pv+Pt csp+Plt=PLt (27)
式中,Pt pv表示t时刻光伏电站实际出力;Pt csp表示t时刻光热电站实际出力;
光伏出力约束模型为:
Figure BDA0003778650720000123
式中,
Figure BDA0003778650720000124
表示光伏预测出力。
本实施例中,太阳倍数SM由步骤S3中求解第一阶段的光热电站综合投资成本最低的目标函数获得。太阳倍数SM与镜场面积SSF的关系如公式(6)所示。将最优太阳倍数代入公式(6),进而可以确定光热电站的最优镜场面积SSF。通过求解光热电站综合投资成本最低的目标函数,还可以获得光热电站的最优储热时间hTES,将最优储热时间hTES代入公式(7)中可以得到储热系统的最大储热容量
Figure BDA0003778650720000121
进一步地,在第二阶段模型中,以联合发电系统经济效益最大和电能供需差波动最小为优化目标,以场景分析法获得的光伏出力预测值和负荷预测值以及第一阶段求得的最优镜场面积与最大储热容量,作为第二阶段t时刻储热量、发电量以及运维成本的约束条件,并兼顾上述各种约束条件,可以确定各时段光热最优发电占比,则每个时段的最优储热容量也随之确定,进而可以最大化地平抑电能供需差的波动,提高整个系统的运行效益。
下面以装机容量为100MW的光伏电站以及装机容量为100MW的光热电站为例对本发明方法进行说明。
本实施例中,光热电站聚光集热系统单位面积投资成本为1000元/m2,光热电站储热系统单位容量投资成本为3×104元/MWh,光热电站发电系统投资成本8000元/kW。光伏电站市场售电电价为0.9元/KWh,光热电站市场售电电价为1.2元/KWh。光伏、光热并网环境效益系数Kpv=Kcsp=230元/MW,弃光惩罚成本系数取1000元/MW,光伏、光热电站运维成本系数分别取20元/MW、30元/MW。光热电站运行参数如表1所示。
表1光热电站运行参数
Figure BDA0003778650720000122
Figure BDA0003778650720000131
采用光伏电站某一典型日预测出力与实际出力数据,数据每1小时采集一次。通过对光伏预测误差及负荷预测误差统计分析,本发明采用拉丁超立方抽样生成960个随机预测误差,再结合预测误差与实际值之间的关系,可得到40组光伏出力预测与负荷预测值。光伏出力预测、负荷预测初始场景生成分别如图4、图5所示。
通过k-means算法,对光伏出力预测、负荷预测初始场景进行缩减;可得到典型的光伏出力预测场景、负荷预测场景;分别如图6、7所示。每个典型场景出现的概率分别如表2、表3所示。
表2光伏出力预测典型场景概率
聚类数 Z1 Z2 Z3 Z4 Z5
占比 12.5% 42.5% 20% 17.5% 7.5%
表3负荷预测典型场景概率
聚类数 Z1 Z2 Z3 Z4 Z5
占比 5% 17.5% 30% 12.5% 3.5%
本发明提出的光伏-光热联合发电系统中光热电站储热容量两阶段优化配置模型是一个混合整数规划问题,借助MATLAB软件,调用CPLEX求解器求解此模型。通过对第一阶段模型求解可得到最优的太阳倍数与储热时间。
为了分析光热电站集热系统中镜场面积对储热系统中储热容量的影响,本发明统计了不同的太阳倍数对应的储热时长,结果如表4所示。
表4不同的太阳倍数对应的储热时长
太阳倍数 储热时间/h 系统运行效益/10<sup>6</sup>元 弃光惩罚成本/10<sup>3</sup>元
1 2.86 0.91 2.66
1.33 3.21 0.94 2.66
2 4.92 1.07 3.11
2.33 3.88 1.11 3.31
3 3.86 1.27 3.53
从表4可以看出来,太阳倍数从1增加到2的过程中,弃光惩罚成本保持不变,光伏电站弃光量没有增加,此阶段光热电站的储热系统处于储热状态,发电系统没有发电,尽可能为光伏发电提供上网空间。随着太阳倍数的增大,即当太阳倍数从2增加到3时,弃光惩罚成本持续增长,在此阶段光热电站集热系统获得的热功率逐渐增大,需要光热电站发电。因此,光伏电站的弃光惩罚成本逐渐增大。从表4同样也可以看出,随着太阳倍数从1增加到2的过程中,储热时间即储热容量呈现出增大的趋势,随后随着太阳倍数增大,储热容量呈现出下降的趋势。
本实施例中,太阳倍数最优为2.33,储热最优时间为3.88h时,光热电站综合投资成本达到最低。此时,将求到的最优太阳倍数与储热时间带入公式(6)与公式(7)中可得到光热电站集热系统镜场面积为1900000m2,储热系统中最大储热容量为1000MWh。将其作为第二阶段优化模型的约束条件,可得到每个时段充、放热功率与储热容量;光热电站各时段储热容量与充放功率如图8所示。
从图8中可以看出,从凌晨1:00至6:00时,此时负荷处于低谷时段,光热电站通过在白天太阳辐射较强时储存的热来发电,此时段光热电站储热系统处于放热状态;在白天7:00至16:00时,此时段太阳辐射逐渐增加,光伏电站开始发电,光热电站此时储存热能,储热系统容量逐渐增大;在晚上18:00至24:00时,光热电站用白天储存的热能来发电。由此可知,由于光热电站具有储热系统,可以允许光热电站将发电转移到日出和日落时段。
本发明提供的一种联合发电系统中光热电站储热容量优化配置方法,针对光伏出力预测和负荷预测的不确定性,采用场景分析法来描述;为了利用光热电站的出力来缓解光伏出力的不确定性,对联合发电系统中的光热电站的储热容量进行合理的配置,在考虑光伏预测出力、负荷预测不确定性的基础上,建立了联合系统中光热电站储热容量两阶段优化配置模型;通过两阶段优化配置模型,兼顾光热电站的约束条件,可以计算出最优的镜场面积与最大储热容量,进而可以计算出每个时段最优的储热容量;最大限度地提高光伏上网电量,同时提高整个系统的运行效益。
以上对本发明所提供的一种联合发电系统中光热电站储热容量优化配置方法进行了详细介绍,本实施例中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上,本说明书内容不应理解为对本发明的限制。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本实施例中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本实施例所示的这些实施例,而是要符合与本实施例所公开的原理和新颖特点相一致的最宽的范围。

Claims (7)

1.一种联合发电系统中光热电站储热容量优化配置方法,其特征在于,包括以下步骤:
S1、针对光伏出力预测和负荷预测的不确定性,采用场景分析法获得出现概率最大的光伏出力预测值和负荷预测值;
S2、建立光热电站储热容量两阶段优化配置模型;
S3、基于光热电站储热容量两阶段优化配置模型,在第一阶段中,建立光热电站综合投资成本最低的目标函数,求解得到最优太阳倍数与储热时间,结合光热电站的太阳倍数模型与储热时间模型,进而得到光热电站的最优镜场面积与储热系统的最大储热容量;
S4、在第二阶段中,建立联合发电系统的经济效益最大的目标函数和电能供需差波动最小的目标函数,结合场景分析法获得的光伏出力预测值和负荷预测值以及第一阶段求得的所述最优镜场面积与所述最大储热容量,计算得到每个时段最优的储热量。
2.根据权利要求1所述的一种联合发电系统中光热电站储热容量优化配置方法,其特征在于,所述步骤S1中,所述采用场景分析法获取出现概率最大的光伏出力预测值和负荷预测值,具体为:
构建光伏出力预测误差和负荷预测误差的概率密度函数;
根据所述概率密度函数,采用拉丁超立方抽样的方法来生成光伏出力预测及负荷预测的初始场景;
采用k-means算法将生成的所述初始场景进行缩减,获得出现概率最大的光伏出力预测值和负荷预测值。
3.根据权利要求2所述的一种联合发电系统中光热电站储热容量优化配置方法,其特征在于,所述步骤S3中,所述光热电站综合投资成本最低的目标函数为:
minf1=FSF+FTES+FPB
式中,f1为光热电站综合投资成本折现值;FSF为集热系统投资成本折现值;FTES为储热系统投资成本折现值;FPB为发电系统投资成本折现值;具体表达式如下:
FSF=CSFSMδCRF(r,n)
FTES=CTEShTESδCRF(r,n)
FPB=CPBPCSP
Figure FDA0003778650710000021
式中,CSF为光热电站太阳倍数投资成本;SM为光热电站太阳倍数;CTES为光热电站储热时间投资成本,hTES为光热电站储热时间;CPB为光热电站发电系统单位容量投资成本;PCSP为光热电站额定容量;δCRF(r,n)为等年值折算系数;r为贴现率;n为光热电站运行年限;
通过求解所述综合投资成本最低的目标函数,确定最优FSF与FTES
根据所述最优FSF与FTES,确定光热电站的最优太阳倍数SM与储热时间hTES
4.根据权利要求3所述的一种联合发电系统中光热电站储热容量优化配置方法,其特征在于,所述步骤S3中,所述光热电站的最优镜场面积与储热系统的最大储热容量的获取过程具体为:
将所述最优太阳倍数SM与储热时间hTES分别带入太阳倍数模型与储热时间模型中,求解得到光热电站的最优镜场面积与储热系统的最大储热容量。
5.根据权利要求4所述的一种联合发电系统中光热电站储热容量优化配置方法,其特征在于,所述太阳倍数模型为:
Figure FDA0003778650710000022
式中,SM为太阳倍数;SSF为镜场面积;IDNI为标准工况下太阳直接辐射强度;ηPB为发电系统热电转换效率;ηSF为集热系统光热转换效率;PCSP为光热电站额定输出功率;
所述储热时间模型为:
Figure FDA0003778650710000031
式中,hTES为储热时间;ηdis为储热系统放热效率;
Figure FDA0003778650710000032
为储热系统最大储热容量。
6.根据权利要求5所述的一种联合发电系统中光热电站储热容量优化配置方法,其特征在于,所述步骤S4中,所述联合发电系统的经济效益最大的目标函数为:
maxf2=F1+F2-F3-F4
式中,f2表示联合系统运行效益;F1表示联合系统市场售电效益;F2表示联合系统并网环境效益;F3表示光伏弃光惩罚成本;F4表示联合系统运行维护成本;具体表达式如下:
Figure FDA0003778650710000033
Figure FDA0003778650710000034
Figure FDA0003778650710000035
Figure FDA0003778650710000036
式中,T表示调度时段;
Figure FDA0003778650710000037
表示t时刻市场售电的分时电价;Cprice为光热电站市场售电电价;Pt pv、Pt csp分别表示t时刻光伏、光热实际出力;Kpv、Kcsp分别表示光伏、光热并网环境效益系数;
Figure FDA0003778650710000038
表示光伏弃光惩罚成本系数;Cpv、CTES和CPB分别表示光伏、光热电站储热系统和光热电站发电系统运维成本系数。
7.根据权利要求6所述的一种联合发电系统中光热电站储热容量优化配置方法,其特征在于,所述步骤S4中,所述电能供需差波动最小的目标函数,由等效负荷方差最小的目标函数来表示,具体为:
Figure FDA0003778650710000041
Plt=PLt-Pt pv-Pt csp
Figure FDA0003778650710000042
式中,f3表示等效负荷方差;Plt表示t时刻等效负荷值;Plt.av表示t时刻等效负荷平均值;PLt表示t时刻负荷预测值;T表示调度时段。
CN202210923499.8A 2022-08-02 2022-08-02 一种联合发电系统中光热电站储热容量优化配置方法 Active CN115276088B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210923499.8A CN115276088B (zh) 2022-08-02 2022-08-02 一种联合发电系统中光热电站储热容量优化配置方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210923499.8A CN115276088B (zh) 2022-08-02 2022-08-02 一种联合发电系统中光热电站储热容量优化配置方法

Publications (2)

Publication Number Publication Date
CN115276088A true CN115276088A (zh) 2022-11-01
CN115276088B CN115276088B (zh) 2023-08-08

Family

ID=83746722

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210923499.8A Active CN115276088B (zh) 2022-08-02 2022-08-02 一种联合发电系统中光热电站储热容量优化配置方法

Country Status (1)

Country Link
CN (1) CN115276088B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116667346A (zh) * 2023-08-02 2023-08-29 国网山东省电力公司广饶县供电公司 分布式光伏源网荷储协同优化方法、系统、终端及介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140039699A1 (en) * 2012-07-31 2014-02-06 Joseph W. Forbes, Jr. System, method, and apparatus for electric power grid and network management of grid elements
CN105488357A (zh) * 2016-01-26 2016-04-13 清华大学 一种光热电站-风电场联合系统的有功功率滚动调度方法
CN106685315A (zh) * 2016-12-21 2017-05-17 中广核太阳能开发有限公司 光伏光热互补发电系统以及发电方法
CN110119850A (zh) * 2019-05-22 2019-08-13 长沙理工大学 基于光热发电调节的储热量双阶段优化调度方法
WO2020261960A1 (ja) * 2019-06-25 2020-12-30 住友電気工業株式会社 集光型太陽光発電装置の監視システム、追尾ずれの検出方法、及び、追尾ずれの検出プログラム
CN114548544A (zh) * 2022-02-18 2022-05-27 国网甘肃省电力公司经济技术研究院 一种光伏光热互补发电系统优化配置方法
CN114723136A (zh) * 2022-04-07 2022-07-08 重庆大学 计及光热电站镜场面积和储热容量的优化运行方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140039699A1 (en) * 2012-07-31 2014-02-06 Joseph W. Forbes, Jr. System, method, and apparatus for electric power grid and network management of grid elements
CN105488357A (zh) * 2016-01-26 2016-04-13 清华大学 一种光热电站-风电场联合系统的有功功率滚动调度方法
CN106685315A (zh) * 2016-12-21 2017-05-17 中广核太阳能开发有限公司 光伏光热互补发电系统以及发电方法
CN110119850A (zh) * 2019-05-22 2019-08-13 长沙理工大学 基于光热发电调节的储热量双阶段优化调度方法
WO2020261960A1 (ja) * 2019-06-25 2020-12-30 住友電気工業株式会社 集光型太陽光発電装置の監視システム、追尾ずれの検出方法、及び、追尾ずれの検出プログラム
CN114548544A (zh) * 2022-02-18 2022-05-27 国网甘肃省电力公司经济技术研究院 一种光伏光热互补发电系统优化配置方法
CN114723136A (zh) * 2022-04-07 2022-07-08 重庆大学 计及光热电站镜场面积和储热容量的优化运行方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CUI YANG;YANG ZHIWEN;ZHONG WUZHI: ""Capacity configuration of thermal energy storagewithin CSP plant considering scheduling economy"", 《2018 INTERNATIONAL CONFERENCE ON POWER SYSTEM TECHNOLOGY (POWERCON)》, pages 1 - 6 *
姚元玺: ""计及调度经济性的光热电站储热容量配置方法"", 《太阳能学报》, vol. 40, no. 11, pages 3039 - 3045 *
毛维宙,陆华军,单宝奇: ""塔式光热电站储热系统容量优化配置"", 《电气技术》, vol. 19, no. 11, pages 37 - 52 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116667346A (zh) * 2023-08-02 2023-08-29 国网山东省电力公司广饶县供电公司 分布式光伏源网荷储协同优化方法、系统、终端及介质
CN116667346B (zh) * 2023-08-02 2023-09-26 国网山东省电力公司广饶县供电公司 分布式光伏源网荷储协同优化方法、系统、终端及介质

Also Published As

Publication number Publication date
CN115276088B (zh) 2023-08-08

Similar Documents

Publication Publication Date Title
CN108321837B (zh) 一种风电-光热联合发电系统及其运行方法
CN111030094B (zh) 一种多能互补园区供能可靠性评估方法
Kong et al. Optimization of the hybrid solar power plants comprising photovoltaic and concentrating solar power using the butterfly algorithm
CN110994606B (zh) 一种基于复杂适应系统理论的多能源电源容量配置方法
CN111740408B (zh) 一种基于鲁棒随机模型的光热电站最优报价决策方法
CN112541609A (zh) 风光热和水蓄能联合可再生能源发电系统容量优化模型
CN112953364A (zh) 一种计及光热电站寿命的光热-风电-光伏的联合系统的运行优化模型
CN114498617A (zh) 一种考虑随机性的多能源发电系统两阶段优化调度方法
CN114069688A (zh) 一种基于时序生产模拟的多电源容量布局规划方法
CN112785065A (zh) 基于混合人工鱼群算法的综合能源系统规划方法及系统
CN110198052B (zh) 一种光热-风电联合并网发电协调控制方法
CN111585279A (zh) 一种基于新能源消纳的微电网优化调度方法
CN115276088A (zh) 一种联合发电系统中光热电站储热容量优化配置方法
CN115577929A (zh) 基于多场景分析的乡村综合能源系统随机优化调度方法
CN111049179A (zh) 一种考虑不确定性的新能源发电系统多目标优化调度方法
CN111969591B (zh) 一种光蓄能水电站互补的离网式混合能源系统优化方法
CN112653137A (zh) 计及碳交易的光热电站与风电系统、低碳调度方法及系统
CN110994698B (zh) 一种太阳能光伏-光热联合发电系统优化运行方法
CN107346474B (zh) 水冷光伏光热一体化系统发电量计算的三维梯形模糊方法
Zhao et al. Research on reliability evaluation of power generation system with solar thermal power
Xingping et al. Operation mode and economic analysis of concentrating solar power station
CN107292768B (zh) 一种光伏发电系统日发电量模糊概率计算方法及装置
CN107332272B (zh) 一种空气冷却光伏-光热发电系统输出功率计算方法
CN116542475B (zh) 一种计及光热电站集热器工作模式的两阶段优化调度方法
CN107330291B (zh) 一种光伏发电日发电量的二型点值Zadeh模糊计算方法及装置

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