CN110503250B - 考虑电-热转移负荷转移量不确定性的综合能源规划方法 - Google Patents

考虑电-热转移负荷转移量不确定性的综合能源规划方法 Download PDF

Info

Publication number
CN110503250B
CN110503250B CN201910729085.XA CN201910729085A CN110503250B CN 110503250 B CN110503250 B CN 110503250B CN 201910729085 A CN201910729085 A CN 201910729085A CN 110503250 B CN110503250 B CN 110503250B
Authority
CN
China
Prior art keywords
year
electric
constraint conditions
node
formula
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
CN201910729085.XA
Other languages
English (en)
Other versions
CN110503250A (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.)
Guangzhou Dayu Chuangfu Technology Co ltd
Original Assignee
Yanshan 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 Yanshan University filed Critical Yanshan University
Priority to CN201910729085.XA priority Critical patent/CN110503250B/zh
Publication of CN110503250A publication Critical patent/CN110503250A/zh
Application granted granted Critical
Publication of CN110503250B publication Critical patent/CN110503250B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems 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)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提出了一种考虑电‑热转移负荷转移量不确定性的综合能源规划方法,其具体步骤如下:步骤1、以能源中心模型为基础,以电能和天然气为输入,以电能和热能供应负荷需求为输出,建立综合能源系统规划模型;步骤2、考虑电‑热可转移负荷转移量对综合能源系统规划的影响,采用两阶段鲁棒优化方法处理电‑热可转移负荷转移量的不确定性;步骤3、求解综合能源系统规划模型,得到规划方案。该方法,可以减少系统碳排放量,提升系统的可靠性和多能互补的效益。另一方面,采用两阶段鲁棒优化理论处理电‑热可转移负荷转移量的不确定性,可以更加精确地考虑其不确定性,灵活调整综合能源系统规划方案的保守性。

Description

考虑电-热转移负荷转移量不确定性的综合能源规划方法
技术领域
本发明涉及综合能源系统规划领域和电力市场领域,具体涉及一种考虑电-热转移负荷转移量不确定性的综合能源规划方法。
背景技术
随着能源需求和环境压力的不断增加,天然气开采技术、电转气(power to gas,P2G)技术和热电联产(combined heat and power,CHP)等技术得到了迅速发展。电力系统和天然气系统的耦合日益加深,电力需求和天然气需求可以相互转化,例如热负荷需求可以由供热机组直接供热,也可以由电加热设备供热,使得这类负荷具有电-热可转移的性质,某些电负荷被热负荷替代可以减少系统CO2排放量。传统的电力系统与天然气系统独立规划的模式已经不能满足电气综合能源系统的需要。因此,为了提升综合能源系统供能的可靠性,在进行综合能源系统规划时就有必要考虑电-热可转移负荷转移量的不确定性。
发明内容
本发明的目的在于提供一种考虑电-热转移负荷转移量不确定性的综合能源规划方法,提高系统的可靠性和整体效益。
为解决上述技术问题,本发明提出一种考虑电-热转移负荷转移量不确定性的综合能源规划方法,所述方法的具体步骤如下:
步骤1、以能源中心模型为基础,以电能和天然气为输入,以电能和热能供应负荷需求为输出,建立综合能源系统规划模型;
步骤2、考虑电-热可转移负荷转移量对综合能源系统规划的影响,采用两阶段鲁棒优化方法处理电-热可转移负荷转移量的不确定性;
步骤3、求解综合能源系统规划模型,得到规划方案。
优选的,所述综合能源系统规划模型包括目标函数和约束条件;
所述目标函数以规划年限内总成本净现值最小为优化目标,所述总成本包括投资成本和运行成本;
所述目标函数表达式为:
Figure BDA0002159919770000021
Figure BDA0002159919770000022
Figure BDA0002159919770000023
式中:F为规划总成本,N为规划年限,r为贴现率,
Figure BDA0002159919770000024
为第t年机组的投资成本,
Figure BDA0002159919770000025
为第t年机组的运行成本,G为候选设备的种类,M为运行设备的种类;ΩG为候选设备集合,ΩY为运行设备集合,
Figure BDA0002159919770000026
为第G类候选设备的额定容量,
Figure BDA0002159919770000027
为第G类候选设备的单位容量投建成本;xi,t为设备i第t年的状态变量,xi,t-1为设备i第t-1年的状态变量,所述设备包括运营设备和候选设备;所述运行设备为已经投入生产运营的设备,包括气源点、常规机组、燃气锅炉、CHP机组和P2G厂站,所述运营设备的状态变量设为1;所述候选设备为尚未投入运营的设备,包括常规机组、燃气锅炉、CHP机组、P2G厂站、输电线路和天然气管道,所述候选设备投建之前的状态变量为0,投建之后的状态变量为1;
Figure BDA0002159919770000028
为第M类运行设备的出力,
Figure BDA0002159919770000029
为第M类运行设备单位容量的运行费用。
优选的,所述约束条件包括系统装机类型及数目约束条件、电力系统运行约束条件、天然气系统运行约束条件和能源耦合单元约束条件。
优选的,所述系统装机类型及数目约束条件为:
dtk=max{0,ωt,kt-1,k} (4)
Figure BDA00021599197700000210
式中:dtk表示第t年k型机组的新增装机数量,ωt,k为k型机组在第t年的数量,ωt-1,k为k型机组在第t-1年的数量,
Figure BDA00021599197700000211
表示k型机组的最大数目。
优选的,所述电力系统运行约束条件包括常规发电机组出力约束条件、风电机组运行约束条件、输电线路有功功率传输极限约束条件、支路潮流约束条件、节点电功率平衡约束条件和备用容量约束条件;
所述常规发电机组出力约束条件为:
Figure BDA0002159919770000031
式中,
Figure BDA0002159919770000032
为常规发电机组i的出力上限,
Figure BDA0002159919770000033
为常规发电机组i的出力下限,Pi gen(t,τ)为第t年时段τ火电机组i的出力;
所述风电机组运行约束条件为:
Figure BDA0002159919770000034
式中,
Figure BDA0002159919770000035
为风电机组i的出力下限,
Figure BDA0002159919770000036
为风电机组i的出力上限,Pi wind(t,τ)为第t年时段τ风电机组i的出力;
所述输电线路有功功率传输极限约束条件为:
Figure BDA0002159919770000037
式中,
Figure BDA0002159919770000038
为输电线路(i,j)传输有功功率的最大值,
Figure BDA0002159919770000039
为输电线路(i,j)传输的有功功率;
所述支路潮流约束条件为:
Figure BDA00021599197700000310
式中,θi(t,τ)为节点i处在第t年时段τ的电压相角;Xij为输电线路(i,j)的电抗;
所述节点电功率平衡约束条件为:
Figure BDA00021599197700000311
式中,Pi CHP(t,τ)为节点i处CHP机组在第t年时段τ输出的有功功率;Pi P2G(t,τ)为节点i处P2G设备在第t年时段τ消耗的有功功率;Pi′(t,τ)为节点i处所需电功率;
Figure BDA00021599197700000312
为输电线路(i,j)在第t年时段τ传输的有功功率;
Figure BDA00021599197700000313
为输电线路(j,i)在第t年时段τ传输的有功功率;
所述备用容量约束条件为:
Figure BDA0002159919770000041
式中,Re(t,τ)为第t年时段τ电力系统所需的备用容量。
优选的,所述天然气系统运行约束条件包括气源点的出气限额约束条件、管道传输流量约束条件、节点天然气平衡约束条件和节点热功率平衡约束条件;
所述出气限额约束条件为:
Figure BDA0002159919770000042
式中,Fi gs(t,τ)表示节点i处气源点在第t年时段τ的进气量或者出气量,所述Fi gs(t,τ)为正值表示进气量,所述Fi gs(t,τ)为负值表示出气量;
Figure BDA0002159919770000043
表示节点i处气源点在第t年时段τ的最小出气量;
Figure BDA0002159919770000044
表示节点i处气源点在第t年时段τ的最大出气量;
所述管道传输流量约束条件为:
Figure BDA0002159919770000045
式中,
Figure BDA0002159919770000046
表示天然气管道(i,j)在第t年时段τ的输送流量;
Figure BDA0002159919770000047
表示天然气管道(i,j)在第t年时段τ的最大输送流量;
所述节点天然气平衡约束条件为:
Figure BDA0002159919770000048
所述节点天然气平衡约束条件表达式中:
Figure BDA0002159919770000049
表示天然气管道(j,i)在第t年时段τ的输送流量;
所述节点热功率平衡约束条件为:
Figure BDA00021599197700000410
式中:
Figure BDA00021599197700000411
表示点i处在第t年时段τCHP机组的供热负荷;
Figure BDA00021599197700000412
表示点i处在第t年时段τ燃气锅炉的供热负荷;H′i(t,τ)表示点i处在第t年时段τ的热负荷。
优选的,所述能源耦合单元约束包括CHP设备、P2G设备和燃气锅炉在运行时需要满足的容量约束和能量转换时的能量守恒约束:
所述CHP设备的能量转换约束条件为:
Figure BDA0002159919770000051
Figure BDA0002159919770000052
式中,
Figure BDA0002159919770000053
表示CHP机组的发电效率;
Figure BDA0002159919770000054
表示CHP机组的供热效率;
所述CHP设备的容量约束条件为:
Figure BDA0002159919770000055
式中,Fi chp(t,τ)表示节点i处CHP机组在第t年时段τ的耗气功率;
Figure BDA0002159919770000056
表示节点i处CHP机组在第t年时段τ的最小耗气功率;
Figure BDA0002159919770000057
表示节点i处CHP机组在第t年时段τ的最大耗气功率;
所述P2G设备的能量转换约束条件为:
Fi p2g(t,τ)=θpPi p2g(t,τ) (19)
式中,θp表示P2G厂站的电转气效率;Fi p2g(t,τ)表示节点i处P2G厂站在第t年时段τ的电转气功率;
所述P2G设备的容量约束条件为:
Figure BDA0002159919770000058
式中,
Figure BDA0002159919770000059
为节点i处P2G设备在第t年时段τ消耗的最小有功功率;
Figure BDA00021599197700000510
为节点i处P2G设备在第t年时段τ消耗的最大有功功率;
所述燃气锅炉设备的能量转换约束条件为:
Figure BDA00021599197700000511
式中,
Figure BDA00021599197700000512
燃气锅炉的供热效率;
Figure BDA00021599197700000513
表示节点i处燃气锅炉在第t年时段τ的供热功率;Fi f(t,τ)表示节点i处燃气锅炉在第t年时段τ的耗气功率;
所述燃气锅炉设备的容量约束条件为:
Figure BDA00021599197700000514
式中,
Figure BDA00021599197700000515
表示节点i处燃气锅炉在第t年时段τ的最小耗气功率;
Figure BDA00021599197700000516
表示节点i处燃气锅炉在第t年时段τ的最大耗气功率。
优选的,所述综合能源系统规划的确定性模型为:
Figure BDA0002159919770000061
所述确定性模型中,x为投资决策变量,y为运行优化变量;所述确定性模型需要满足的约束条件包括所述约束条件(4)、(5)以及所述约束条件(6)-(22)。
优选的,所述步骤2中,采用两阶段鲁棒优化方法处理后的电-热可转移负荷转移量的不确定性模型,包括电负荷折算公式和电-热可转移负荷转移量的波动范围:
所述电负荷折算公式为:
Figure BDA0002159919770000062
式中:P为供电量;
Figure BDA0002159919770000063
为发热折算成供电量的折算系数,即发1kWh电所需的热量,取6MJ/kWh;H为供热量;
所述电-热可转移负荷转移量的波动范围为:
Figure BDA0002159919770000064
Figure BDA0002159919770000065
式中,
Figure BDA0002159919770000066
为第t年电负荷需求可转移到热负荷需求的量的潜在值,
Figure BDA0002159919770000067
为第t年电负荷需求可转移到热负荷需求的量的预测值,
Figure BDA0002159919770000068
为第t年电负荷需求可转移到热负荷需求的量的波动偏差,在不考虑所述电-热可转移负荷转移量的不确定性时,
Figure BDA0002159919770000069
的初始值为0;Γz为规划期内电-热负荷可转移量的潜在值的不确定性调节参数。
优选的,所述步骤2中,采用两阶段鲁棒规划模型方法处理处理电-热可转移负荷转移量的不确定性,得到min-max-min结构的两阶段鲁棒优化形式的所述目标函数表达式为:
Figure BDA00021599197700000610
所述目标函数表达式(27)需要满足约束条件(4)-(22)和(25),且需要采用Benders算法对所述目标函数表达式(27)进行转化和求解。
与现有的技术相比,本发明的技术方案所带来的有益效果是:
本发明提供的考虑电-热转移负荷转移量不确定性的综合能源规划方法,通过考虑电-热可转移负荷,可以减少系统碳排放量,提升系统的可靠性和多能互补的效益。另一方面,采用两阶段鲁棒优化理论处理电-热可转移负荷转移量的不确定性,可以更加精确地考虑其不确定性,灵活调整综合能源系统规划方案的保守性,提升多能互补的效益。
附图说明
图1是本发明实施例中综合能源规划方法的流程图;
图2是能源中心结构示意图;以及
图3是本发明实施例中9节点综合能源测试系统结构图。
具体实施方式
下面结合附图对本发明做进一步说明。
本实施例的考虑电-热转移负荷转移量不确定性的综合能源规划方法,具体步骤如图1所示:
步骤一,以图2所示的能源中心结构为基础,输入能源为电能和天然气,输出电能和热能供应负荷需求。
步骤二,建立综合能源系统规划模型。
步骤三,考虑电-热可转移负荷转移量对综合能源系统规划的影响,采用两阶段鲁棒优化方法处理电-热可转移负荷转移量的不确定性。
步骤四,求解综合能源系统规划模型,得到规划方案。
作为上述技术方案的补充,综合能源系统规划模型包括目标函数和约束条件:
以规划年限内总成本净现值最小为优化目标,其中总成本包含投资成本、运行成本两部分:
Figure BDA0002159919770000071
Figure BDA0002159919770000081
Figure BDA0002159919770000082
式中:F为规划总成本;N为规划年限;r为贴现率;
Figure BDA0002159919770000083
为第t年机组的投资成本;
Figure BDA0002159919770000084
为第t年机组的运行成本;
Figure BDA0002159919770000085
Figure BDA0002159919770000086
分别为第t年的投资成本和运行成本;G和M分别为候选设备与运行设备的种类;ΩG和ΩY分别为候选设备和运行设备集合(候选设备包含常规机组、燃气锅炉、CHP机组、P2G厂站、输电线路和天然气管道,运行设备包含气源点、常规机组、燃气锅炉、CHP机组和P2G厂站);
Figure BDA0002159919770000087
分别表示第G类候选设备的额定容量、单位容量投建成本;xi,t、xi,t-1为不同设备的状态变量,已有设备的状态变量设为1,候选设备第t年投运后,其状态变量由0变为1;
Figure BDA0002159919770000088
Figure BDA0002159919770000089
分别为第M类运行设备的出力和单位容量的运行费用。
约束条件主要包括系统装机类型及数目约束、电力系统运行约束、天然气系统运行约束、能源耦合单元约束四个方面。
(1)系统装机类型及数目约束:
dtk=max{0,ωt,kt-1,k} (4)
Figure BDA00021599197700000810
式中:dtk表示第t年k型机组的新增装机数量;ωt,k和ωt-1,k分别为k型机组在第t年和第t-1年的数量;
Figure BDA00021599197700000811
表示k型机组的最大数目。
(2)电力系统运行约束包括常规发电机组出力约束、风电机组运行约束、输电线路有功功率传输极限约束、支路潮流约束、节点电功率平衡约束、备用容量约束:
Figure BDA00021599197700000812
Figure BDA00021599197700000813
Figure BDA0002159919770000091
Figure BDA0002159919770000092
Figure BDA0002159919770000093
Figure BDA0002159919770000094
式中,
Figure BDA0002159919770000095
Figure BDA0002159919770000096
分别表示常规发电机组i的出力上下限;Pi gen(t,τ)为第t年时段τ火电机组i的出力;
Figure BDA0002159919770000097
Figure BDA0002159919770000098
分别表示风电机组i的出力上下限;Pi wind(t,τ)为第t年时段τ风电机组i的出力;
Figure BDA0002159919770000099
表示输电线路(i,j)传输有功功率的最大值;
Figure BDA00021599197700000910
表示输电线路(i,j)传输的有功功率;θi(t,τ)为节点i处在第t年时段τ的电压相角;Xij为输电线路(i,j)的电抗;Pi CHP(t,τ)为节点i处CHP机组在第t年时段τ输出的有功功率;Pi P2G(t,τ)为节点i处P2G设备在第t年时段τ消耗的有功功率;Re(t,τ)为第t年时段τ电力系统所需的备用容量。
(3)天然气系统运行约束包括:气源点的出气限额约束见式、管道传输流量约束见式、节点天然气平衡约束见式、节点热功率平衡约束。
Figure BDA00021599197700000911
Figure BDA00021599197700000912
Figure BDA00021599197700000913
Figure BDA00021599197700000914
式中:Fi gs(t,τ)表示节点i处气源点在第t年时段τ的进/出气量(正值表示进气,负值表示出气);
Figure BDA00021599197700000915
表示天然气管道(i,j)在第t年时段τ的输送流量;H′i(t,τ)表示点i处在第t年时段τ的热负荷。
(4)能源耦合单元约束包括CHP、P2G设备和燃气锅炉在运行时需要满足的容量约束和能量转换时的能量守恒约束:
Figure BDA0002159919770000101
Figure BDA0002159919770000102
Fi p2g(t,τ)=θpPi p2g(t,τ) (18)
Figure BDA0002159919770000103
Figure BDA0002159919770000104
Figure BDA0002159919770000105
Figure BDA0002159919770000106
式中:
Figure BDA0002159919770000107
θp分别表示CHP机组的发电和供热效率、P2G厂站的电转气效率;Fi chp(t,τ)、Fi p2g(t,τ)分别表示节点i处CHP机组在第t年时段τ的耗气功率、节点i处P2G厂站在第t年时段τ的电转气功率;
Figure BDA0002159919770000108
代表CHP机组在第t年时段τ的供热功率;
Figure BDA0002159919770000109
Fi f(t,τ)、
Figure BDA00021599197700001010
分别表示节点i处燃气锅炉在第t年时段τ的供热功率、节点i处燃气锅炉在第t年时段τ的耗气功率、燃气锅炉的供热效率。
作为上述技术方案的补充,综合能源系统规划的确定性模型为:
综合能源系统规划模型的目标函数为式(1);投资决策变量为x,运行优化变量为y,表达式为:
Figure BDA00021599197700001011
需要满足的约束条件包括:与规划层面相关的约束有式(4)-(5),与运行层面相关的约束有式(6)-(22)。
作为上述技术方案的补充,电-热负荷可转移量不确定性建模,电负荷折算公式和电-热可转移负荷转移量的波动范围如下:
Figure BDA0002159919770000111
Figure BDA0002159919770000112
Figure BDA0002159919770000113
式中:P为供电量;
Figure BDA0002159919770000114
为发热折算成供电量的折算系数,即发1kWh电所需的热量,取6MJ/kWh;H为供热量;
Figure BDA0002159919770000115
Figure BDA0002159919770000116
分别为第t年电负荷需求可转移到热负荷需求的量的潜在值、预测值和波动偏差,假设
Figure BDA0002159919770000117
Γz为规划期内电-热负荷可转移量的潜在值的不确定性调节参数。
作为上述技术方案的补充,综合能源系统规划的两阶段鲁棒规划模型的建立,将目标函数式(1)转变为min-max-min结构的两阶段鲁棒优化形式:
Figure BDA0002159919770000118
并且需要满足约束式(4)-(22)、式(25)。采用Benders算法对模型进行转化和求解。
第五步:采用Benders算法对模型进行转化求解,调用CPLEX获得综合能源系统规划最优方案。
本实施例以如图3所示的一个9节点能源中心测试系统为例,来解释本发明的实际应用。
9节点能源中心测试系统结构图如附图1所示,该能源中心测试系统包含3个火电机组、2个气源点、1个风电场、4个燃气锅炉、9条输电线路和9条天然气管道。给定规划周期为10年,年均折现率为6%,规划期内电负荷和热负荷的年均增长率分别为5%和2%。
电-热负荷可转移量不确定性方面,设第t年电负荷需求可转移到热负荷需求的量的预测值为0,即
Figure BDA0002159919770000119
不确定域的波动范围
Figure BDA00021599197700001110
参数Γz在1~10之间浮动。
表1 综合能源系统参数
Figure BDA00021599197700001111
根据本发明所提出的方法,求解综合能源系统的规划结果,为说明考虑电-热转移负荷转移量不确定性的综合能源规划方法的有效性与优越性,本发明分析了以下两种场景:
场景1:不考虑电-热转移负荷转移量,即电-热转移负荷转移量为0。
场景2:考虑电-热转移负荷转移量,及其不确定性。
针对上述2种场景,求解得规划结果如表2所示。
场景 总成本(106$) 投资成(106$) 运行成本(106$) 碳排放量(104t)
场景1 421.73 35.97 385.76 163.25
场景2 426.45 38.91 387.54 126.83
比较上述2种场景的结果,在综合能源系统的规划中考虑电-热转移负荷转移量不确定性后,虽然投资成本、运行成本略有升高,但是碳排放量大幅减少。
本发明考虑了电-热转移负荷转移量不确定性对综合能源规划的影响,充分挖掘了多种能源形式之间互补互济、灵活规划的潜力,可以减少碳排放量,提高综合能源系统的低碳性和灵活性。
尽管上面结合附图对本发明的功能及工作过程进行了描述,但本发明并不局限于上述的具体功能和过程,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可以做出很多形式,这些均属于本发明的保护之内。

Claims (10)

1.一种考虑电-热转移负荷转移量不确定性的综合能源规划方法,其特征在于:所述方法的具体步骤如下:
步骤1、以能源中心模型为基础,以电能和天然气为输入,以电能和热能供应负荷需求为输出,建立综合能源系统规划模型;
步骤2、考虑电-热可转移负荷转移量对综合能源系统规划的影响,采用两阶段鲁棒优化方法处理电-热可转移负荷转移量的不确定性;
步骤3、求解综合能源系统规划模型,得到规划方案。
2.根据权利要求1所述的考虑电-热转移负荷转移量不确定性的综合能源规划方法,其特征在于:所述综合能源系统规划模型包括目标函数和约束条件;
所述目标函数以规划年限内总成本净现值最小为优化目标,所述总成本包括投资成本和运行成本;
所述目标函数表达式为:
Figure FDA0004046573430000011
Figure FDA0004046573430000012
Figure FDA0004046573430000013
式中:F为规划总成本,N为规划年限,r为贴现率,
Figure FDA0004046573430000014
为第t年机组的投资成本,
Figure FDA0004046573430000015
为第t年机组的运行成本,G为候选设备的种类,M为运行设备的种类;ΩG为候选设备集合,ΩM为运行设备集合,
Figure FDA0004046573430000016
为第G类候选设备的额定容量,
Figure FDA0004046573430000017
为第G类候选设备的单位容量投建成本;xi,t为设备i第t年的状态变量,xi,t-1为设备i第t-1年的状态变量,所述设备包括运营设备和候选设备;所述运行设备为已经投入生产运营的设备,包括气源点、常规机组、燃气锅炉、CHP机组和P2G厂站,所述运营设备的状态变量设为1;所述候选设备为尚未投入运营的设备,包括常规机组、燃气锅炉、CHP机组、P2G厂站、输电线路和天然气管道,所述候选设备投建之前的状态变量为0,投建之后的状态变量为1;
Figure FDA0004046573430000018
为第M类运行设备的出力,
Figure FDA0004046573430000019
为第M类运行设备单位容量的运行费用。
3.根据权利要求2所述的考虑电-热转移负荷转移量不确定性的综合能源规划方法,其特征在于:所述约束条件包括系统装机类型及数目约束条件、电力系统运行约束条件、天然气系统运行约束条件和能源耦合单元约束条件。
4.根据权利要求3所述的考虑电-热转移负荷转移量不确定性的综合能源规划方法,其特征在于:所述系统装机类型及数目约束条件为:
dtk=max{0,ωt,kt-l,k} (4)
Figure FDA0004046573430000021
式中:dtk表示第t年k型机组的新增装机数量,ωt,k为k型机组在第t年的数量,ωt-1,k为k型机组在第t-1年的数量,
Figure FDA0004046573430000022
表示k型机组的最大数目。
5.根据权利要求4所述的考虑电-热转移负荷转移量不确定性的综合能源规划方法,其特征在于:所述电力系统运行约束条件包括常规发电机组出力约束条件、风电机组运行约束条件、输电线路有功功率传输极限约束条件、支路潮流约束条件、节点电功率平衡约束条件和备用容量约束条件;
所述常规发电机组出力约束条件为:
Figure FDA0004046573430000023
式中,
Figure FDA0004046573430000024
为常规发电机组i的出力下限,
Figure FDA0004046573430000025
为常规发电机组i的出力上限,Pi gen(t,τ)为第t年时段τ火电机组i的出力;
所述风电机组运行约束条件为:
Figure FDA0004046573430000026
式中,
Figure FDA0004046573430000027
为风电机组i的出力下限,
Figure FDA0004046573430000028
为风电机组i的出力上限,Pi wind(t,τ)为第t年时段τ风电机组i的出力;
所述输电线路有功功率传输极限约束条件为:
Figure FDA0004046573430000029
式中,
Figure FDA00040465734300000210
为输电线路(i,j)传输有功功率的最大值,
Figure FDA00040465734300000211
为输电线路(i,j)传输的有功功率;
所述支路潮流约束条件为:
Figure FDA0004046573430000031
式中,θi(t,τ)为节点i处在第t年时段τ的电压相角;Xij为输电线路(i,j)的电抗;
所述节点电功率平衡约束条件为:
Figure FDA0004046573430000032
式中,Pi CHP(t,τ)为节点i处CHP机组在第t年时段τ输出的有功功率;Pi P2G(t,τ)为节点i处P2G设备在第t年时段τ消耗的有功功率;Pi′(t,τ)为节点i处所需电功率;
Figure FDA0004046573430000033
为输电线路(i,j)在第t年时段τ传输的有功功率;
Figure FDA0004046573430000034
为输电线路(j,i)在第t年时段τ传输的有功功率;
所述备用容量约束条件为:
Figure FDA0004046573430000035
式中,Re(t,τ)为第t年时段τ电力系统所需的备用容量。
6.根据权利要求5所述的考虑电-热转移负荷转移量不确定性的综合能源规划方法,其特征在于:所述天然气系统运行约束条件包括气源点的出气限额约束条件、管道传输流量约束条件、节点天然气平衡约束条件和节点热功率平衡约束条件;
所述出气限额约束条件为:
Figure FDA0004046573430000036
式中,Fi gs(t,τ)表示节点i处气源点在第t年时段τ的进气量或者出气量,所述Fi gs(t,τ)为正值表示进气量,所述Fi gs(t,τ)为负值表示出气量;
Figure FDA0004046573430000037
表示节点i处气源点在第t年时段τ的最小出气量;
Figure FDA0004046573430000038
表示节点i处气源点在第t年时段τ的最大出气量;
所述管道传输流量约束条件为:
Figure FDA0004046573430000039
式中,
Figure FDA00040465734300000310
表示天然气管道(i,j)在第t年时段τ的输送流量;
Figure FDA00040465734300000311
表示天然气管道(i,j)在第t年时段τ的最大输送流量;
所述节点天然气平衡约束条件为:
Figure FDA0004046573430000041
所述节点天然气平衡约束条件表达式中:
Figure FDA0004046573430000042
表示天然气管道(j,i)在第t年时段τ的输送流量;
所述节点热功率平衡约束条件为:
Figure FDA0004046573430000043
式中:
Figure FDA0004046573430000044
表示点i处在第t年时段τCHP机组的供热负荷;
Figure FDA0004046573430000045
表示点i处在第t年时段τ燃气锅炉的供热负荷;Hi'(t,τ)表示点i处在第t年时段τ的热负荷。
7.根据权利要求6所述的考虑电-热转移负荷转移量不确定性的综合能源规划方法,其特征在于:所述能源耦合单元约束包括CHP设备、P2G设备和燃气锅炉在运行时需要满足的容量约束和能量转换时的能量守恒约束:
所述CHP设备的能量转换约束条件为:
Figure FDA0004046573430000046
Figure FDA0004046573430000047
式中,
Figure FDA0004046573430000048
表示CHP机组的发电效率;
Figure FDA0004046573430000049
表示CHP机组的供热效率;
所述CHP设备的容量约束条件为:
Figure FDA00040465734300000410
式中,Fi chp(t,τ)表示节点i处CHP机组在第t年时段τ的耗气功率;
Figure FDA00040465734300000411
表示节点i处CHP机组在第t年时段τ的最小耗气功率;
Figure FDA00040465734300000412
表示节点i处CHP机组在第t年时段τ的最大耗气功率;
所述P2G设备的能量转换约束条件为:
Fi p2g(t,τ)=θpPi p2g(t,τ)(19)
式中,θp表示P2G厂站的电转气效率;Fi p2g(t,τ)表示节点i处P2G厂站在第t年时段τ的电转气功率;
所述P2G设备的容量约束条件为:
Figure FDA0004046573430000051
式中,
Figure FDA0004046573430000052
为节点i处P2G设备在第t年时段τ消耗的最小有功功率;
Figure FDA0004046573430000053
为节点i处P2G设备在第t年时段τ消耗的最大有功功率;
所述燃气锅炉设备的能量转换约束条件为:
Figure FDA0004046573430000054
式中,
Figure FDA0004046573430000055
燃气锅炉的供热效率;
Figure FDA0004046573430000056
表示节点i处燃气锅炉在第t年时段τ的供热功率;Fi f(t,τ)表示节点i处燃气锅炉在第t年时段τ的耗气功率;
所述燃气锅炉设备的容量约束条件为:
Figure FDA0004046573430000057
式中,
Figure FDA0004046573430000058
表示节点i处燃气锅炉在第t年时段τ的最小耗气功率;
Figure FDA0004046573430000059
表示节点i处燃气锅炉在第t年时段τ的最大耗气功率。
8.根据权利要求7所述的考虑电-热转移负荷转移量不确定性的综合能源规划方法,其特征在于:所述综合能源系统规划的确定性模型为:
Figure FDA00040465734300000510
所述确定性模型中,x为投资决策变量,y为运行优化变量;所述确定性模型需要满足的约束条件包括所述约束条件(4)、(5)以及所述约束条件(6)-(22)。
9.根据权利要求8所述的考虑电-热转移负荷转移量不确定性的综合能源规划方法,其特征在于:所述步骤2中,采用两阶段鲁棒优化方法处理后的电-热可转移负荷转移量的不确定性模型,包括电负荷折算公式和电-热可转移负荷转移量的波动范围:
所述电负荷折算公式为:
Figure FDA00040465734300000511
式中:P为供电量;
Figure FDA00040465734300000512
为发热折算成供电量的折算系数,即发1kWh电所需的热量,取6MJ/kWh;H为供热量;
所述电-热可转移负荷转移量的波动范围为:
Figure FDA0004046573430000061
Figure FDA0004046573430000062
式中,
Figure FDA0004046573430000063
为第t年电负荷需求可转移到热负荷需求的量的潜在值,
Figure FDA0004046573430000064
为第t年电负荷需求可转移到热负荷需求的量的预测值,
Figure FDA0004046573430000065
为第t年电负荷需求可转移到热负荷需求的量的波动偏差,在不考虑所述电-热可转移负荷转移量的不确定性时,
Figure FDA0004046573430000066
的初始值为0;Γz为规划期内电-热负荷可转移量的潜在值的不确定性调节参数。
10.根据权利要求9所述的考虑电-热转移负荷转移量不确定性的综合能源规划方法,其特征在于:所述步骤2中,采用两阶段鲁棒规划模型方法处理电-热可转移负荷转移量的不确定性,得到min-max-min结构的两阶段鲁棒优化形式的所述目标函数表达式为:
Figure FDA0004046573430000067
所述目标函数表达式(27)需要满足约束条件(4)-(22)和(25),且需要采用Benders算法对所述目标函数表达式(27)进行转化和求解。
CN201910729085.XA 2019-08-08 2019-08-08 考虑电-热转移负荷转移量不确定性的综合能源规划方法 Active CN110503250B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910729085.XA CN110503250B (zh) 2019-08-08 2019-08-08 考虑电-热转移负荷转移量不确定性的综合能源规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910729085.XA CN110503250B (zh) 2019-08-08 2019-08-08 考虑电-热转移负荷转移量不确定性的综合能源规划方法

Publications (2)

Publication Number Publication Date
CN110503250A CN110503250A (zh) 2019-11-26
CN110503250B true CN110503250B (zh) 2023-04-07

Family

ID=68587121

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910729085.XA Active CN110503250B (zh) 2019-08-08 2019-08-08 考虑电-热转移负荷转移量不确定性的综合能源规划方法

Country Status (1)

Country Link
CN (1) CN110503250B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111030110B (zh) * 2019-12-28 2021-05-18 福州大学 一种考虑电转气消纳风电的电力-天然气耦合系统鲁棒协同调度方法
CN111049135B (zh) * 2019-12-30 2023-02-03 国网吉林省电力有限公司 一种多区域电气耦合系统的分布式两阶段合作运行方法
CN111222694A (zh) * 2019-12-31 2020-06-02 新奥数能科技有限公司 考虑负荷预测不确定性的综合能源系统优化方法和装置
CN111738498B (zh) * 2020-06-08 2024-05-24 国网安徽省电力有限公司经济技术研究院 一种综合能源系统鲁棒规划方法及系统
CN111799777A (zh) * 2020-06-11 2020-10-20 国网山东省电力公司经济技术研究院 一种考虑天然气与电力耦合的综合能源规划方法
CN112016195B (zh) * 2020-08-10 2022-11-29 浙江大学 一种考虑供能可靠性的电气耦合系统柔性规划方法
CN112907030B (zh) * 2021-01-20 2023-12-19 国网山东省电力公司寿光市供电公司 一种考虑需求侧响应的能源中心配置方法及系统
CN114362168B (zh) * 2022-01-11 2024-02-13 云南电网有限责任公司电力科学研究院 一种能源互联系统的设备选型方法
CN115907155B (zh) * 2022-11-22 2023-07-18 南方电网能源发展研究院有限责任公司 引入碳排放成本的发电机组规划方法以及装置
CN115906488B (zh) * 2022-11-25 2024-01-23 国网江苏省电力有限公司电力科学研究院 一种用于楼宇综合能源系统的优化调度方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109149656A (zh) * 2018-09-19 2019-01-04 南京师范大学 一种气-电互联综合能源系统机组组合方法
CN109255471A (zh) * 2018-08-17 2019-01-22 国网山东省电力公司电力科学研究院 一种含风电的电-气-热综合能源系统扩展规划优化方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105633948B (zh) * 2015-12-09 2018-05-29 长沙理工大学 一种分布式能源接入电力系统随机模糊潮流算法
CN106097154A (zh) * 2016-07-12 2016-11-09 广东电网有限责任公司电力科学研究院 一种工业园区能源中心的调度方法及系统
CN108062022A (zh) * 2016-11-05 2018-05-22 南京理工大学 一种热电协同系统的优化控制方法
CN107730129B (zh) * 2017-10-24 2021-06-22 重庆大学 考虑光热热电联产与电锅炉的电-气-热互联系统风险评估方法
CN108805328A (zh) * 2018-04-30 2018-11-13 国网浙江省电力有限公司经济技术研究院 光热电站热电联供微网系统的优化运行方法
CN110009152B (zh) * 2019-04-03 2021-12-07 东南大学 一种考虑电转气和不确定性的区域综合能源系统运行鲁棒优化方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109255471A (zh) * 2018-08-17 2019-01-22 国网山东省电力公司电力科学研究院 一种含风电的电-气-热综合能源系统扩展规划优化方法
CN109149656A (zh) * 2018-09-19 2019-01-04 南京师范大学 一种气-电互联综合能源系统机组组合方法

Also Published As

Publication number Publication date
CN110503250A (zh) 2019-11-26

Similar Documents

Publication Publication Date Title
CN110503250B (zh) 考虑电-热转移负荷转移量不确定性的综合能源规划方法
CN110245878B (zh) 智能楼宇群体的分布式综合能源需求响应协同优化方法
CN109978625B (zh) 一种计及电热气网络的综合能源系统多目标运行优化方法
CN110571789B (zh) 数据驱动下基于风电不确定性的电热气网三阶段调度方法
CN110210747B (zh) 一种电热气互联能源系统灵活性调度方法
CN111652441B (zh) 考虑气电联合需求响应的气电综合能源系统配网优化方法
CN103426032B (zh) 一种热电联产机组的经济优化调度方法
CN112701687B (zh) 考虑价格型联合需求响应的气电配网系统鲁棒优化运行方法
CN111222257B (zh) 基于凸优化的电热水多能流协同调度方法
CN106055773B (zh) 结合热网模型的多区域综合能源系统配置模型的建立方法
CN107862418A (zh) 基于储能能源枢纽的热电耦合系统优化调度方法及装置
CN114266165B (zh) 考虑碳排放的蒸汽动力系统中蒸汽透平布局优化方法
CN115455719A (zh) 一种基于全生命周期碳成本的综合能源系统优化仿真方法
CN110244568B (zh) 工业企业微网的能源枢纽模型及其多能互补优化控制方法
CN113806952A (zh) 一种考虑源-荷-储的冷热电综合能源系统及其优化运行方法
CN111342453B (zh) 计及多类型备用资源的电气综合能源系统备用决策方法
Varganova et al. Energy-efficient distribution of heat between the boiler units of industrial power palnts with the use of computers
CN112564182A (zh) 一种基于分布鲁棒机会约束的多能联合系统统一调度方法
CN115222102B (zh) 输配一体化综合能源系统的分布式调度方法
CN113191083B (zh) 考虑全工况外部参数变化的烟气余热回收系统优化设计方法
Hamedi et al. Operation and flexibility contribution of natural gas tanks in multi-carrier energy systems
Makwana Energy Retrofit Debottlenecking of Total Sites
Khokhlov et al. The economic effect determination for combined cycle power plant with pressurized heat recovery steam generator
CN112968479B (zh) 考虑p2x变工况特性的电力系统调度运行方法
CN109800480A (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
TR01 Transfer of patent right

Effective date of registration: 20230822

Address after: Room 801, 85 Kefeng Road, Huangpu District, Guangzhou City, Guangdong Province

Patentee after: Guangzhou Dayu Chuangfu Technology Co.,Ltd.

Address before: 066004 No. 438 west section of Hebei Avenue, seaport District, Hebei, Qinhuangdao

Patentee before: Yanshan University

TR01 Transfer of patent right