CN115310259B - 基于软行动者-批评者的园区综合能源系统梯级优化方法 - Google Patents

基于软行动者-批评者的园区综合能源系统梯级优化方法 Download PDF

Info

Publication number
CN115310259B
CN115310259B CN202210762777.6A CN202210762777A CN115310259B CN 115310259 B CN115310259 B CN 115310259B CN 202210762777 A CN202210762777 A CN 202210762777A CN 115310259 B CN115310259 B CN 115310259B
Authority
CN
China
Prior art keywords
power
energy
energy storage
gas
heat
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
Application number
CN202210762777.6A
Other languages
English (en)
Other versions
CN115310259A (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.)
North China Electric Power University
Original Assignee
North China Electric Power University
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN202210762777.6A priority Critical patent/CN115310259B/zh
Publication of CN115310259A publication Critical patent/CN115310259A/zh
Application granted granted Critical
Publication of CN115310259B publication Critical patent/CN115310259B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • 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/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/06Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • General Business, Economics & Management (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Engine Equipment That Uses Special Cycles (AREA)

Abstract

一种基于软行动者‑批评者的园区综合能源系统梯级优化方法:对园区综合能源系统内的负荷按能源品位进行划分:热负荷划分为高品位的蒸汽负荷、高温热负荷和低品位的中温热负荷;构建园区综合能源系统的燃气轮机与余热锅炉的热电联产模型、燃气锅炉的制热模型、变压器和换热器的转换模型、电热泵模型、尖峰加热器的转换模型和储能设备模型;建立园区综合能源系统的能源供需平衡模型;建立园区综合能源系统梯级优化模型;设计智能体的状态、动作和奖励函数,运用软行动者‑批评者对园区综合能源系统梯级优化模型进行求解,得到电力、天然气以及热力的优化运行结果。本发明实现了各品位能源的合理梯级利用,提升了能源利用效率。

Description

基于软行动者-批评者的园区综合能源系统梯级优化方法
技术领域
本发明涉及一种园区综合能源系统梯级优化方法。特别是涉及一种基于软行动者-批评者的园区综合能源系统梯级优化方法。
背景技术
低碳经济发展是全球为应对气候变化而提出的全新经济发展模式,能源利用效率的提高将助力全球能源系统的低碳化转型。综合能源系统耦合了电,热,气等多类能源形式,在运行过程中需要通过优化多类型能源设备出力实现综合能源系统经济运行,大大提高了能源利用效率。
现有研究从时空维度针对综合能源系统的协调优化运行提供了可行方案,但未从供能结构角度对提高能源转化效率深入探讨。从求解方法角度分析,现有研究所采用求解方法往往基于历史数据,依赖各类可再生能源出力以及用户负荷的精确预测,且对于大规模综合能源系统高维优化问题的求解速度较慢。
随着人工智能和电力系统及综合能源系统的融合日益紧密,将人工智能方法应用到综合能源系统优化运行领域,能够有效提升复杂系统应对多元不确定性的能力,从而提高综合能源系统运营的安全性和经济性。
发明内容
本发明所要解决的技术问题是,为了克服现有技术的不足,提供一种考虑多种能源品位的基于软行动者-批评者的园区综合能源系统梯级优化方法。
本发明所采用的技术方案是:一种基于软行动者-批评者的园区综合能源系统梯级优化方法,包括如下步骤:
1)对园区综合能源系统内的负荷按能源品位进行划分:将园区综合能源系统内的负荷分为电负荷、气负荷和热负荷,其中热负荷划分为高品位的蒸汽负荷、高温热负荷和低品位的中温热负荷;
2)构建园区综合能源系统的能源转换设备模型,包括有;燃气轮机与余热锅炉的热电联产模型、燃气锅炉的制热模型、变压器和换热器的转换模型、电热泵模型、尖峰加热器的转换模型和储能设备模型;
3)建立园区综合能源系统的能源供需平衡模型;
4)建立园区综合能源系统梯级优化模型,包括,以经济成本F最小为园区综合能源系统梯级优化模型的目标函数,以网络购能约束和设备运行约束为约束条件;
5)针对园区综合能源系统梯级优化模型设计智能体的状态、动作和奖励函数,运用软行动者-批评者对园区综合能源系统梯级优化模型进行求解,得到电力、天然气以及热力的优化运行结果。
步骤2)中:
(1)所述的燃气轮机与余热锅炉的热电联产模型为:
PGT=ηGT,eHngGGT
HGT,m=ηGT,reηHR(1-ηGT,e)HngGGT
HGT,st=ηGT,st(1-ηGT,e)HngGGT
其中,GGT为燃气轮机的进气量;PGT为燃气轮机发出的电功率;HGT,m为余热锅炉通过回收燃气轮机的废气产生的中温热功率;HGT,st为燃气轮机抽气蒸汽功率;ηGT,e,ηGT,re,ηGT,st分别为燃气轮机的发电效率、余热比例系数和抽气比例系数;ηHR为余热锅炉余热回收效率;Hng为天然气的低位燃烧热值;
(2)燃气锅炉的制热模型为:
HGB=ηGBHngGGB
其中,GGB为燃气锅炉的进气量;HGB为燃气锅炉产生的高温蒸汽所含热能;ηGB为燃气锅炉产热效率;
(3)变压器和换热器的转换模型为:
PT=ηTPgrid
HHE=ηHEHgrid
其中,Pgrid和Hgrid分别为购电功率和购热功率;PT和HHE分别为变压器的输出电功率和换热器的输出热功率;ηT和ηHE分别为变压器和换热器的转换效率;
(4)电热泵模型为:
HHP=CHPPHP
其中,PHP,HHP分别为电热泵耗电功率和制热功率;CHP为电热泵能效比;
(5)尖峰加热器的转换模型为:
HPLC,out=ηPLCHPLC,in
其中,HPLC,in,HPLC,out分别为尖峰加热器消耗的蒸汽功率和产生的高温热功率;ηPLC为尖峰加热器转换效率;
(6)储能设备模型为:
EES(t)=EES(t-1)(1-δES)+ηESPES(t)Δt
EGS(t)=EGS(t-1)(1-δGS)+ηGSGGS(t)Δt
EHS(t)=EHS(t-1)(1-δHS)+ηHSHHS(t)Δt
其中,EES(t),EGS(t),EHS(t)分别为当前时刻的电储能能量状态、气储能能量状态和热储能能量状态;EES(t-1),EGS(t-1),EHS(t-1)分别为上一时刻的电储能能量状态、气储能能量状态和热储能能量状态;PES(t),GGS(t),HHS(t)分别为当前时刻的电储能充/放能功率、气储能充/放能功率和热储能充/放能功率;δES,δGS,δHS分别为电储能自损率、气储能自损率和热储能自损率;ηESGSHS分别为电储能的充/放能系数、气储能的充/放能系数和热储能的充/放能系数;Δt为时间间隔。
步骤3)中所述的园区综合能源系统的能源供需平衡模型表示为如下公式:
Figure GDA0003877384300000031
其中,Le,Lg,Lm,Lh,Lst,Ll分别为电负荷、气负荷、中温热负荷、高温热负荷、蒸汽负荷和低温热负荷;Pgrid,Ggrid,Hgrid分别为购电功率、购气功率和购热功率;Pwt和Ppv分别为风力发电功率和光伏发电功率;GGT,HHR,PHP,GGB,HPLC,in分别为燃气轮机的进气量、余热锅炉收集的热功率、电热泵耗电功率、燃气锅炉的进气量和尖峰加热器消耗的蒸汽功率;PES,GGS,HHS分别为电储能充/放能功率、气储能充/放能功率和热储能充/放能功率;ηT,ηHE分别为变压器和换热器的转换效率;ηGT,e,ηGT,re,ηGT,st,ηHR,ηGB,ηPLC,CHP分别为燃气轮机的发电效率、余热比例系数、抽气比例系数、余热锅炉余热回收效率、燃气锅炉产热效率、尖峰加热器转换效率和电热泵能效比;ηESGSHS分别为电储能的充/放能系数、气储能的充/放能系数和热储能的充/放能系数;Hng为天然气的低位燃烧热值。
步骤4)中所述的目标函数表示为:
Figure GDA0003877384300000032
其中,
Figure GDA0003877384300000033
分别为购能成本以及储能的充放电折旧成本;
购能成本
Figure GDA0003877384300000034
具体形式为:
Figure GDA0003877384300000035
其中,Pgrid(t),Ggrid(t),Hgrid(t)分别为当前时刻的购电功率、购气功率和购热功率;εe(t),εs(t),εh(t)分别为当前时刻的购电单价、购气单价和购热单价;Δt为时间间隔;
储能充放电折旧成本
Figure GDA0003877384300000036
具体形式为:
Figure GDA0003877384300000037
Figure GDA0003877384300000038
其中,
Figure GDA0003877384300000039
分别为当前时刻的电储能充放能折旧成本、气储能充放能折旧成本和热储能充放能折旧成本;ρES,ρGS,ρHS分别表示电储能的折旧系数、气储能的折旧系数和热储能的折旧系数;PES(t),GGS(t),HHS(t)分别为当前时刻的电储能充/放能功率、气储能充/放能功率和热储能充/放能功率。
步骤4)中所述的
(1)网络购能约束为:
Figure GDA0003877384300000041
其中,Pgrid(t),Ggrid(t),Hgrid(t)分别为当前时刻的购电功率、购气功率和购热功率;
Figure GDA0003877384300000042
和/>
Figure GDA0003877384300000043
为购电功率上、下限;/>
Figure GDA0003877384300000044
和/>
Figure GDA0003877384300000045
为购气功率上、下限;/>
Figure GDA0003877384300000046
和/>
Figure GDA0003877384300000047
为购热功率上、下限;
(2)设备运行约束为:
Figure GDA0003877384300000048
Figure GDA0003877384300000049
Figure GDA00038773843000000410
Figure GDA00038773843000000411
Figure GDA00038773843000000412
其中,PGT(t),HGB(t),PHP(t),HPLC,out(t)分别为当前时刻燃气轮机发出的电功率、燃气锅炉产生的高温蒸汽所含热能、电热泵耗电功率和尖峰加热器产生的高温热功率;EES(t),EGS(t),EHS(t)分别为当前时刻的电储能能量状态、气储能能量状态和热储能能量状态;
Figure GDA00038773843000000413
和/>
Figure GDA00038773843000000414
为燃气轮机发出电功率上、下限;/>
Figure GDA00038773843000000415
和/>
Figure GDA00038773843000000416
为燃气锅炉发出热功率上、下限;/>
Figure GDA00038773843000000417
和/>
Figure GDA00038773843000000418
为电热泵耗电功率上、下限;/>
Figure GDA00038773843000000419
和/>
Figure GDA00038773843000000420
为尖峰加热器产生的高温热功率上、下限;/>
Figure GDA00038773843000000421
和/>
Figure GDA00038773843000000422
为电储能的能量状态上、下限;/>
Figure GDA00038773843000000423
和/>
Figure GDA00038773843000000424
为气储能的能量状态上、下限;/>
Figure GDA00038773843000000425
和/>
Figure GDA00038773843000000426
为热储能的能量状态上、下限。
步骤5)中所述的针对园区综合能源系统梯级优化模型设计智能体的状态、动作和奖励函数,包括:
(1)智能体的状态st设计为:
st={Le(t),Lg(t),Lm(t),Lh(t),Lst(t),Pwt(t),Ppv(t),EES(t),EGS(t),EHS(t),t}
其中,Le(t),Lg(t),Lm(t),Lh(t),Lst(t)分别为当前时刻的电负荷、气负荷、中温热负荷、高温热负荷和蒸汽负荷;Pwt(t)和Ppv(t)分别为当前时刻的风力发电功率和光伏发电功率;EES(t),EGS(t),EHS(t)分别为当前时刻的电储能能量状态、气储能能量状态和热储能能量状态;t表示当前时刻;
(2)智能体的动作at设计为:
at={PGT(t),PHP(t),HGB(t),HPLC,out(t),PES(t),GGS(t),HHS(t),Pgrid(t),Ggrid(t),Hgrid(t)}
其中,PGT(t),PHP(t),HGB(t),HPLC,out(t)分别为当前时刻燃气轮机发出的电功率、电热泵耗电功率、燃气锅炉产生的高温蒸汽所含热能和尖峰加热器产生的高温热功率;PES(t),GGS(t),HHS(t)分别为当前时刻的电储能充/放能功率、气储能充/放能功率和热储能充/放能功率;Pgrid(t),Ggrid(t),Hgrid(t)分别为当前时刻的购电功率、购气功率和购热功率;
(3)智能体的奖励函数rt设计为:
Figure GDA0003877384300000051
其中,
Figure GDA0003877384300000052
分别为购能成本以及储能的充放电折旧成本;α1,α2分别为购能成本的权重和储能的充放电折旧成本的权重。
步骤5)中所述的运用软行动者-批评者对园区综合能源系统梯级优化模型进行求解,得到电力、天然气以及热力的优化运行结果,表示为在训练过程中寻找最优策略π*
Figure GDA0003877384300000053
其中,π为智能体策略;st为智能体的状态;at为智能体的动作;rt为智能体的奖励函数;ρπ为智能体依照智能体策略π生成的状态-动作组集合;H(π(·|st))为智能体依照智能体策略π输出的动作熵值;α为熵值对应参数。
本发明的基于软行动者-批评者的园区综合能源系统梯级优化方法,具有如下优点:
1、本发明所提出的园区综合能源系统梯级优化模型结合园区综合能源系统负荷多样化需求,从供能结构角度进行考虑,将负荷按能源品位进行划分,实现了各品位能源的合理梯级利用,提升了能源利用效率。
2、本发明所设计方法充分发挥了多能互补优势,提高了综合能源系统的经济性。
3、本发明所采用的软行动者-批评者方法能够对设备的连续动作进行优化控制,有效避免了动作空间的“维数灾”问题。
4、本发明所设计方法无需对源荷进行精准预测,能够有效提升复杂系统应对多元不确定性的能力,从而提高了综合能源系统运营的安全性。
5、本发明所提出的园区综合能源系统梯级优化运行方法可有效实现园区综合能源系统经济运行,对工程应用具有实际意义。
附图说明
图1是本发明基于软行动者-批评者的园区综合能源系统梯级优化方法的流程图;
图2是本发明实例中新能源出力曲线、负荷曲线和电价曲线示意图;
图3是本发明实例中电力优化运行结果示意图;
图4是本发明实例中天然气优化运行结果示意图;
图5是本发明实例中热力优化运行结果示意图。
具体实施方式
下面结合实施例和附图对本发明的基于软行动者-批评者的园区综合能源系统梯级优化方法做出详细说明。
如图1所示,本发明的基于软行动者-批评者的园区综合能源系统梯级优化方法,包括如下步骤:
1)对园区综合能源系统内的负荷按能源品位进行划分:将园区综合能源系统内的负荷分为电负荷、气负荷和热负荷,其中热负荷划分为高品位的蒸汽负荷、高温热负荷和低品位的中温热负荷;
2)构建园区综合能源系统的能源转换设备模型,包括有;燃气轮机与余热锅炉的热电联产模型、燃气锅炉的制热模型、变压器和换热器的转换模型、电热泵模型、尖峰加热器的转换模型和储能设备模型;其中:
(1)所述的燃气轮机与余热锅炉的热电联产模型为:
PGT=ηGT,eHngGGT
HGT,m=ηGT,reηHR(1-ηGT,e)HngGGT
HGT,st=ηGT,st(1-ηGT,e)HngGGT
其中,GGT为燃气轮机的进气量;PGT为燃气轮机发出的电功率;HGT,m为余热锅炉通过回收燃气轮机的废气产生的中温热功率;HGT,st为燃气轮机抽气蒸汽功率;ηGT,e,ηGT,re,ηGT,st分别为燃气轮机的发电效率、余热比例系数和抽气比例系数;ηHR为余热锅炉余热回收效率;Hng为天然气的低位燃烧热值;
(2)燃气锅炉的制热模型为:
HGB=ηGBHngGGB
其中,GGB为燃气锅炉的进气量;HGB为燃气锅炉产生的高温蒸汽所含热能;ηGB为燃气锅炉产热效率;
(3)变压器和换热器的转换模型为:
PT=ηTPgrid
HHE=ηHEHgrid
其中,Pgrid和Hgrid分别为购电功率和购热功率;PT和HHE分别为变压器的输出电功率和换热器的输出热功率;ηT和ηHE分别为变压器和换热器的转换效率;
(4)电热泵模型为:
HHP=CHPPHP
其中,PHP,HHP分别为电热泵耗电功率和制热功率;CHP为电热泵能效比;
(5)尖峰加热器的转换模型为:
HPLC,out=ηPLCHPLC,in
其中,HPLC,in,HPLC,out分别为尖峰加热器消耗的蒸汽功率和产生的高温热功率;ηPLC为尖峰加热器转换效率;
(6)储能设备模型为:
EES(t)=EES(t-1)(1-δES)+ηESPES(t)Δt
EGS(t)=EGS(t-1)(1-δGS)+ηGSGGS(t)Δt
EHS(t)=EHS(t-1)(1-δHS)+ηHSHHS(t)Δt
其中,EES(t),EGS(t),EHS(t)分别为当前时刻的电储能能量状态、气储能能量状态和热储能能量状态;EES(t-1),EGS(t-1),EHS(t-1)分别为上一时刻的电储能能量状态、气储能能量状态和热储能能量状态;PES(t),GGS(t),HHS(t)分别为当前时刻的电储能充/放能功率、气储能充/放能功率和热储能充/放能功率;δES,δGS,δHS分别为电储能自损率、气储能自损率和热储能自损率;ηESGSHS分别为电储能的充/放能系数、气储能的充/放能系数和热储能的充/放能系数;Δt为时间间隔。
3)建立园区综合能源系统的能源供需平衡模型;所述的园区综合能源系统的能源供需平衡模型表示为如下公式:
Figure GDA0003877384300000071
其中,Le,Lg,Lm,Lh,Lst,Ll分别为电负荷、气负荷、中温热负荷、高温热负荷、蒸汽负荷和低温热负荷;Pgrid,Ggrid,Hgrid分别为购电功率、购气功率和购热功率;Pwt和Ppv分别为风力发电功率和光伏发电功率;GGT,HHR,PHP,GGB,HPLC,in分别为燃气轮机的进气量、余热锅炉收集的热功率、电热泵耗电功率、燃气锅炉的进气量和尖峰加热器消耗的蒸汽功率;PES,GGS,HHS分别为电储能充/放能功率、气储能充/放能功率和热储能充/放能功率;ηT,ηHE分别为变压器和换热器的转换效率;ηGT,e,ηGT,re,ηGT,st,ηHR,ηGB,ηPLC,CHP分别为燃气轮机的发电效率、余热比例系数、抽气比例系数、余热锅炉余热回收效率、燃气锅炉产热效率、尖峰加热器转换效率和电热泵能效比;ηESGSHS分别为电储能的充/放能系数、气储能的充/放能系数和热储能的充/放能系数;Hng为天然气的低位燃烧热值。
4)建立园区综合能源系统梯级优化模型,包括,以经济成本F最小为园区综合能源系统梯级优化模型的目标函数,以网络购能约束和设备运行约束为约束条件;其中,
(1)所述的目标函数表示为:
Figure GDA0003877384300000072
其中,
Figure GDA0003877384300000081
分别为购能成本以及储能的充放电折旧成本;
购能成本
Figure GDA0003877384300000082
具体形式为:
Figure GDA0003877384300000083
其中,Pgrid(t),Ggrid(t),Hgrid(t)分别为当前时刻的购电功率、购气功率和购热功率;εe(t),εs(t),εh(t)分别为当前时刻的购电单价、购气单价和购热单价;Δt为时间间隔;
储能充放电折旧成本
Figure GDA0003877384300000084
具体形式为:
Figure GDA0003877384300000085
Figure GDA0003877384300000086
其中,
Figure GDA0003877384300000087
分别为当前时刻的电储能充放能折旧成本、气储能充放能折旧成本和热储能充放能折旧成本;ρES,ρGS,ρHS分别表示电储能的折旧系数、气储能的折旧系数和热储能的折旧系数;PES(t),GGS(t),HHS(t)分别为当前时刻的电储能充/放能功率、气储能充/放能功率和热储能充/放能功率。
(2)所述的网络购能约束为:
Figure GDA0003877384300000088
其中,Pgrid(t),Ggrid(t),Hgrid(t)分别为当前时刻的购电功率、购气功率和购热功率;
Figure GDA0003877384300000089
和/>
Figure GDA00038773843000000810
为购电功率上、下限;/>
Figure GDA00038773843000000811
和/>
Figure GDA00038773843000000812
为购气功率上、下限;/>
Figure GDA00038773843000000813
和/>
Figure GDA00038773843000000814
为购热功率上、下限;
(3)所述的设备运行约束为:
Figure GDA00038773843000000815
Figure GDA00038773843000000816
Figure GDA00038773843000000817
Figure GDA00038773843000000818
Figure GDA00038773843000000819
其中,PGT(t),HGB(t),PHP(t),HPLC,out(t)分别为当前时刻燃气轮机发出的电功率、燃气锅炉产生的高温蒸汽所含热能、电热泵耗电功率和尖峰加热器产生的高温热功率;EES(t),EGS(t),EHS(t)分别为当前时刻的电储能能量状态、气储能能量状态和热储能能量状态;
Figure GDA00038773843000000820
和/>
Figure GDA00038773843000000821
为燃气轮机发出电功率上、下限;/>
Figure GDA00038773843000000822
和/>
Figure GDA00038773843000000823
为燃气锅炉发出热功率上、下限;/>
Figure GDA00038773843000000824
和/>
Figure GDA00038773843000000825
为电热泵耗电功率上、下限;/>
Figure GDA0003877384300000091
和/>
Figure GDA0003877384300000092
为尖峰加热器产生的高温热功率上、下限;/>
Figure GDA0003877384300000093
和/>
Figure GDA0003877384300000094
为电储能的能量状态上、下限;/>
Figure GDA0003877384300000095
和/>
Figure GDA0003877384300000096
为气储能的能量状态上、下限;/>
Figure GDA0003877384300000097
和/>
Figure GDA0003877384300000098
为热储能的能量状态上、下限。
5)针对园区综合能源系统梯级优化模型设计智能体的状态、动作和奖励函数,运用软行动者-批评者对园区综合能源系统梯级优化模型进行求解,得到电力、天然气以及热力的优化运行结果;其中,
(1)所述的针对园区综合能源系统梯级优化模型设计智能体的状态、动作和奖励函数,包括:
(1.1)智能体的状态st设计为:
st={Le(t),Lg(t),Lm(t),Lh(t),Lst(t),Pwt(t),Ppv(t),EES(t),EGS(t),EHS(t),t}
其中,Le(t),Lg(t),Lm(t),Lh(t),Lst(t)分别为当前时刻的电负荷、气负荷、中温热负荷、高温热负荷和蒸汽负荷;Pwt(t)和Ppv(t)分别为当前时刻的风力发电功率和光伏发电功率;EES(t),EGS(t),EHS(t)分别为当前时刻的电储能能量状态、气储能能量状态和热储能能量状态;t表示当前时刻;
(1.2)智能体的动作at设计为:
at={PGT(t),PHP(t),HGB(t),HPLC,out(t),PES(t),GGS(t),HHS(t),Pgrid(t),Ggrid(t),Hgrid(t)}
其中,PGT(t),PHP(t),HGB(t),HPLC,out(t)分别为当前时刻燃气轮机发出的电功率、电热泵耗电功率、燃气锅炉产生的高温蒸汽所含热能和尖峰加热器产生的高温热功率;PES(t),GGS(t),HHS(t)分别为当前时刻的电储能充/放能功率、气储能充/放能功率和热储能充/放能功率;Pgrid(t),Ggrid(t),Hgrid(t)分别为当前时刻的购电功率、购气功率和购热功率;
(1.3)智能体的奖励函数rt设计为:
Figure GDA0003877384300000099
其中,
Figure GDA00038773843000000910
分别为购能成本以及储能的充放电折旧成本;α1,α2分别为购能成本的权重和储能的充放电折旧成本的权重。
(2)所述的运用软行动者-批评者对园区综合能源系统梯级优化模型进行求解,得到电力、天然气以及热力的优化运行结果,表示为在训练过程中寻找最优策略π*
Figure GDA00038773843000000911
其中,π为智能体策略;st为智能体的状态;at为智能体的动作;rt为智能体的奖励函数;ρπ为智能体依照智能体策略π生成的状态-动作组集合;H(π(·|st))为智能体依照智能体策略π输出的动作熵值;α为熵值对应参数。
下面给出实例:
采用本发明基于软行动者-批评者的园区综合能源系统梯级优化方法,对图2所示的新能源出力曲线、负荷曲线和电价曲线进行梯级优化运行,得到如图3、图4、图5所示的优化运行结果。其中,
图3给出了电力优化运行结果,电价谷段时,电力负荷基本由风电及外网购电满足,燃气轮机发出少量电能,电热泵消耗较多电能用于发热,电储能充电以备后续电价升高。而在电价的平段以及峰段时,光伏发出部分电能,燃气轮机由于其热电联产带来的经济性在此时段几乎满出力运行,电热泵消耗少量电能,电储能放电赚取峰谷电价差。
图4给出了天然气优化运行结果,谷电价时燃气轮机由于较低负荷运行,消耗少量天然气,燃气锅炉消耗较多天然气,平电价以及峰电价时,燃气轮机消耗大量天然气而燃气锅炉消耗少量天然气。
图5给出了热力优化运行结果,高温热负荷完全由尖峰加热器满足,燃气轮机在谷电价时发出少量蒸汽热能并结合余热锅炉发出少量中温热能,在平电价以及峰电价时发出大量蒸汽热能以及中温热能,燃气锅炉在谷电价时发出大量蒸汽热能,在平电价和峰电价时发出少量蒸汽热能,部分蒸汽热能输入尖峰加热器用于生产高温热能,电热泵在谷电价时生产大量中温热能,在平电价和峰电价时生产少量中温热能,其余中温热能缺额从热力网购能补足。

Claims (5)

1.一种基于软行动者-批评者的园区综合能源系统梯级优化方法,其特征在于,包括如下步骤:
1)对园区综合能源系统内的负荷按能源品位进行划分:将园区综合能源系统内的负荷分为电负荷、气负荷和热负荷,其中热负荷划分为高品位的蒸汽负荷、高温热负荷和低品位的中温热负荷;
2)构建园区综合能源系统的能源转换设备模型,包括有;燃气轮机与余热锅炉的热电联产模型、燃气锅炉的制热模型、变压器和换热器的转换模型、电热泵模型、尖峰加热器的转换模型和储能设备模型;其中:
(1)所述的燃气轮机与余热锅炉的热电联产模型为:
PGT=ηGT,eHngGGT
HGT,m=ηGT,reηHR(1-ηGT,e)HngGGT
HGT,st=ηGT,st(1-ηGT,e)HngGGT
其中,GGT为燃气轮机的进气量;PGT为燃气轮机发出的电功率;HGT,m为余热锅炉通过回收燃气轮机的废气产生的中温热功率;HGT,st为燃气轮机抽气蒸汽功率;ηGT,e,ηGT,re,ηGT,st分别为燃气轮机的发电效率、余热比例系数和抽气比例系数;ηHR为余热锅炉余热回收效率;Hng为天然气的低位燃烧热值;
(2)燃气锅炉的制热模型为:
HGB=ηGBHngGGB
其中,GGB为燃气锅炉的进气量;HGB为燃气锅炉产生的高温蒸汽所含热能;ηGB为燃气锅炉产热效率;
(3)变压器和换热器的转换模型为:
PT=ηTPgrid
HHE=ηHEHgrid
其中,Pgrid和Hgrid分别为购电功率和购热功率;PT和HHE分别为变压器的输出电功率和换热器的输出热功率;ηT和ηHE分别为变压器和换热器的转换效率;
(4)电热泵模型为:
HHP=CHPPHP
其中,PHP,HHP分别为电热泵耗电功率和制热功率;CHP为电热泵能效比;
(5)尖峰加热器的转换模型为:
HPLC,out=ηPLCHPLC,in
其中,HPLC,in,HPLC,out分别为尖峰加热器消耗的蒸汽功率和产生的高温热功率;ηPLC为尖峰加热器转换效率;
(6)储能设备模型为:
EES(t)=EES(t-1)(1-δES)+ηESPES(t)Δt
EGS(t)=EGS(t-1)(1-δGS)+ηGSGGS(t)Δt
EHS(t)=EHS(t-1)(1-δHS)+ηHSHHS(t)Δt
其中,EES(t),EGS(t),EHS(t)分别为当前时刻的电储能能量状态、气储能能量状态和热储能能量状态;EES(t-1),EGS(t-1),EHS(t-1)分别为上一时刻的电储能能量状态、气储能能量状态和热储能能量状态;PES(t),GGS(t),HHS(t)分别为当前时刻的电储能充/放能功率、气储能充/放能功率和热储能充/放能功率;δES,δGS,δHS分别为电储能自损率、气储能自损率和热储能自损率;ηESGSHS分别为电储能的充/放能系数、气储能的充/放能系数和热储能的充/放能系数;Δt为时间间隔;t表示当前时刻;
3)建立园区综合能源系统的能源供需平衡模型;所述的园区综合能源系统的能源供需平衡模型表示为如下公式:
Figure FDA0004238572470000021
其中,Le,Lg,Lm,Lh,Lst,Ll分别为电负荷、气负荷、中温热负荷、高温热负荷、蒸汽负荷和低温热负荷;Pgrid,Ggrid,Hgrid分别为购电功率、购气功率和购热功率;Pwt和Ppv分别为风力发电功率和光伏发电功率;GGT,HHR,PHP,GGB,HPLC,in分别为燃气轮机的进气量、余热锅炉收集的热功率、电热泵耗电功率、燃气锅炉的进气量和尖峰加热器消耗的蒸汽功率;PES,GGS,HHS分别为电储能充/放能功率、气储能充/放能功率和热储能充/放能功率;ηT,ηHE分别为变压器和换热器的转换效率;ηGT,e,ηGT,re,ηGT,st,ηHR,ηGB,ηPLC,CHP分别为燃气轮机的发电效率、余热比例系数、抽气比例系数、余热锅炉余热回收效率、燃气锅炉产热效率、尖峰加热器转换效率和电热泵能效比;ηESGSHS分别为电储能的充/放能系数、气储能的充/放能系数和热储能的充/放能系数;Hng为天然气的低位燃烧热值;
4)建立园区综合能源系统梯级优化模型,包括,以经济成本F最小为园区综合能源系统梯级优化模型的目标函数,以网络购能约束和设备运行约束为约束条件;
5)针对园区综合能源系统梯级优化模型设计智能体的状态、动作和奖励函数,运用软行动者-批评者对园区综合能源系统梯级优化模型进行求解,得到电力、天然气以及热力的优化运行结果。
2.根据权利要求1所述的基于软行动者-批评者的园区综合能源系统梯级优化方法,其特征在于,步骤4)中所述的目标函数表示为:
Figure FDA0004238572470000031
其中,
Figure FDA0004238572470000032
分别为购能成本以及储能的充放电折旧成本;t表示当前时刻;
购能成本
Figure FDA0004238572470000033
具体形式为:
Figure FDA0004238572470000034
其中,Pgrid(t),Ggrid(t),Hgrid(t)分别为当前时刻的购电功率、购气功率和购热功率;εe(t),εs(t),εh(t)分别为当前时刻的购电单价、购气单价和购热单价;Δt为时间间隔;
储能充放电折旧成本
Figure FDA0004238572470000035
具体形式为:
Figure FDA0004238572470000036
Figure FDA0004238572470000037
其中,
Figure FDA0004238572470000038
分别为当前时刻的电储能充放能折旧成本、气储能充放能折旧成本和热储能充放能折旧成本;ρES,ρGS,ρHS分别表示电储能的折旧系数、气储能的折旧系数和热储能的折旧系数;PES(t),GGS(t),HHS(t)分别为当前时刻的电储能充/放能功率、气储能充/放能功率和热储能充/放能功率。
3.根据权利要求1所述的基于软行动者-批评者的园区综合能源系统梯级优化方法,其特征在于,步骤4)中所述的
(1)网络购能约束为:
Figure FDA0004238572470000039
其中,Pgrid(t),Ggrid(t),Hgrid(t)分别为当前时刻的购电功率、购气功率和购热功率;
Figure FDA00042385724700000310
和/>
Figure FDA00042385724700000311
为购电功率上、下限;/>
Figure FDA00042385724700000312
和/>
Figure FDA00042385724700000313
为购气功率上、下限;/>
Figure FDA00042385724700000314
和/>
Figure FDA00042385724700000315
为购热功率上、下限;t表示当前时刻;
(2)设备运行约束为:
Figure FDA00042385724700000316
Figure FDA00042385724700000317
Figure FDA00042385724700000318
Figure FDA00042385724700000319
Figure FDA0004238572470000041
其中,PGT(t),HGB(t),PHP(t),HPLC,out(t)分别为当前时刻燃气轮机发出的电功率、燃气锅炉产生的高温蒸汽所含热能、电热泵耗电功率和尖峰加热器产生的高温热功率;EES(t),EGS(t),EHS(t)分别为当前时刻的电储能能量状态、气储能能量状态和热储能能量状态;
Figure FDA0004238572470000042
Figure FDA0004238572470000043
为燃气轮机发出电功率上、下限;/>
Figure FDA0004238572470000044
和/>
Figure FDA0004238572470000045
为燃气锅炉发出热功率上、下限;/>
Figure FDA0004238572470000046
Figure FDA0004238572470000047
为电热泵耗电功率上、下限;/>
Figure FDA0004238572470000048
和/>
Figure FDA0004238572470000049
为尖峰加热器产生的高温热功率上、下限;
Figure FDA00042385724700000410
和/>
Figure FDA00042385724700000411
为电储能的能量状态上、下限;/>
Figure FDA00042385724700000412
和/>
Figure FDA00042385724700000413
为气储能的能量状态上、下限;/>
Figure FDA00042385724700000414
Figure FDA00042385724700000415
为热储能的能量状态上、下限。
4.根据权利要求1所述的基于软行动者-批评者的园区综合能源系统梯级优化方法,其特征在于,步骤5)中所述的针对园区综合能源系统梯级优化模型设计智能体的状态、动作和奖励函数,包括:
(1)智能体的状态st设计为:
st={Le(t),Lg(t),Lm(t),Lh(t),Lst(t),Pwt(t),Ppv(t),EES(t),EGS(t),EHS(t),t}
其中,Le(t),Lg(t),Lm(t),Lh(t),Lst(t)分别为当前时刻的电负荷、气负荷、中温热负荷、高温热负荷和蒸汽负荷;Pwt(t)和Ppv(t)分别为当前时刻的风力发电功率和光伏发电功率;EES(t),EGS(t),EHS(t)分别为当前时刻的电储能能量状态、气储能能量状态和热储能能量状态;t表示当前时刻;
(2)智能体的动作at设计为:
at={PGT(t),PHP(t),HGB(t),HPLC,out(t),PES(t),GGS(t),HHS(t),Pgrid(t),Ggrid(t),Hgrid(t)}
其中,PGT(t),PHP(t),HGB(t),HPLC,out(t)分别为当前时刻燃气轮机发出的电功率、电热泵耗电功率、燃气锅炉产生的高温蒸汽所含热能和尖峰加热器产生的高温热功率;PES(t),GGS(t),HHS(t)分别为当前时刻的电储能充/放能功率、气储能充/放能功率和热储能充/放能功率;Pgrid(t),Ggrid(t),Hgrid(t)分别为当前时刻的购电功率、购气功率和购热功率;
(3)智能体的奖励函数rt设计为:
Figure FDA00042385724700000416
其中,
Figure FDA00042385724700000417
分别为购能成本以及储能的充放电折旧成本;α1,α2分别为购能成本的权重和储能的充放电折旧成本的权重。
5.根据权利要求1所述的基于软行动者-批评者的园区综合能源系统梯级优化方法,其特征在于,步骤5)中所述的运用软行动者-批评者对园区综合能源系统梯级优化模型进行求解,得到电力、天然气以及热力的优化运行结果,表示为在训练过程中寻找最优策略π*
Figure FDA0004238572470000051
其中,π为智能体策略;st为智能体的状态;at为智能体的动作;rt(st,at)为智能体的奖励函数;ρπ为智能体依照智能体策略π生成的状态-动作组集合;H(π(·|st))为智能体依照智能体策略π输出的动作熵值;α为熵值对应参数;t表示当前时刻。
CN202210762777.6A 2022-06-30 2022-06-30 基于软行动者-批评者的园区综合能源系统梯级优化方法 Active CN115310259B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210762777.6A CN115310259B (zh) 2022-06-30 2022-06-30 基于软行动者-批评者的园区综合能源系统梯级优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210762777.6A CN115310259B (zh) 2022-06-30 2022-06-30 基于软行动者-批评者的园区综合能源系统梯级优化方法

Publications (2)

Publication Number Publication Date
CN115310259A CN115310259A (zh) 2022-11-08
CN115310259B true CN115310259B (zh) 2023-06-30

Family

ID=83854822

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210762777.6A Active CN115310259B (zh) 2022-06-30 2022-06-30 基于软行动者-批评者的园区综合能源系统梯级优化方法

Country Status (1)

Country Link
CN (1) CN115310259B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022048127A1 (zh) * 2020-09-04 2022-03-10 华北电力大学 一种电热泵-热电联合系统的优化调控方法及系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104181900B (zh) * 2014-09-04 2017-01-18 徐雪松 一种多能源介质分层动态调控方法
CN106845701B (zh) * 2017-01-11 2019-11-08 东南大学 一种基于热网和房屋热惯性的综合能源系统优化方法
WO2019211528A2 (fr) * 2018-05-02 2019-11-07 Buendia Jose La conduction du champ electrique par la physique quantique
US20210133376A1 (en) * 2019-11-04 2021-05-06 Global Energy Interconnection Research Institute Co. Ltd Systems and methods of parameter calibration for dynamic models of electric power systems

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022048127A1 (zh) * 2020-09-04 2022-03-10 华北电力大学 一种电热泵-热电联合系统的优化调控方法及系统

Also Published As

Publication number Publication date
CN115310259A (zh) 2022-11-08

Similar Documents

Publication Publication Date Title
Yan et al. An integrated design for hybrid combined cooling, heating and power system with compressed air energy storage
CN108629462B (zh) 含储能的综合能源微网优化规划方法及综合能源微网系统
CN108154309B (zh) 计及冷热电多负荷动态响应的能源互联网经济调度方法
CN109474025B (zh) 一种园区级综合能源系统优化调度模型
CN109242350A (zh) 计及可平移负荷的冷热电联供系统容量优化配置方法
CN111027846B (zh) 一种考虑热氢联产的电氢综合能源系统及其容量配置方法
CN108521132B (zh) 多能互补支撑电网频率调节的多时间尺度优化控制方法
CN113159407B (zh) 基于区域综合能源系统的多能储存模块容量优化配置方法
CN110084465B (zh) 基于储能的风力发电系统成本/供电可靠性评估方法
CN110286590A (zh) 考虑多类型能源存储的冷热电联供系统日前优化调度方法
CN117081143A (zh) 促进分布式光伏就地消纳的园区综合能源系统协调优化运行方法
CN113141005B (zh) 一种面向新能源消纳的综合能源系统多时间尺度调度方法
CN109217373B (zh) 商业园区混合能源系统的源储荷协调优化方法
CN115310259B (zh) 基于软行动者-批评者的园区综合能源系统梯级优化方法
Wang et al. Thermodynamic and economic analysis of a multi-energy complementary distributed cchp system coupled with solar thermochemistry and active energy storage regulation process
CN116502921A (zh) 一种园区综合能源系统优化管理系统及其协调调度方法
CN116316888A (zh) 氢电耦合系统的优化调度方法、系统及装置
Zhu et al. Joint optimal scheduling model of clean energy considering CSP regulation
CN115660142A (zh) 一种园区综合能源系统源荷储协调优化调度方法
CN114386256A (zh) 一种考虑电热设备灵活性约束及热网特性的区域电热系统优化调度方法
CN110363397B (zh) 一种基于可转换自由度的综合能源系统规划方法
Wei et al. The day-ahead optimization scheduling of CCHP energy system with multi-mode energy storages
CN113722895A (zh) 一种基于多站融合的综合能源系统优化配置方法
Wang et al. Research on multi-objective planning and optimization of integrated energy system based on economy and environmental protection
CN114362152B (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