CN110729720A - 一种蓄热罐容量规划方法和电转气参数敏感性分析方法 - Google Patents
一种蓄热罐容量规划方法和电转气参数敏感性分析方法 Download PDFInfo
- Publication number
- CN110729720A CN110729720A CN201910991901.4A CN201910991901A CN110729720A CN 110729720 A CN110729720 A CN 110729720A CN 201910991901 A CN201910991901 A CN 201910991901A CN 110729720 A CN110729720 A CN 110729720A
- Authority
- CN
- China
- Prior art keywords
- storage tank
- gas
- heat storage
- cost
- 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.)
- Pending
Links
- 238000005338 heat storage Methods 0.000 title claims abstract description 111
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000010206 sensitivity analysis Methods 0.000 title claims description 8
- 238000010276 construction Methods 0.000 claims abstract description 23
- 238000003860 storage Methods 0.000 claims abstract description 17
- 230000005611 electricity Effects 0.000 claims abstract description 16
- 239000002699 waste material Substances 0.000 claims abstract description 7
- 239000007789 gas Substances 0.000 claims description 98
- 238000006243 chemical reaction Methods 0.000 claims description 40
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 31
- 238000010977 unit operation Methods 0.000 claims description 19
- 238000005457 optimization Methods 0.000 claims description 17
- 239000003345 natural gas Substances 0.000 claims description 15
- 238000013486 operation strategy Methods 0.000 claims description 15
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 238000010586 diagram Methods 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- RZZGGSCKACYXBQ-UHFFFAOYSA-N 3,5,6,7,8,9-hexahydro-2h-cyclohepta[c]pyridazin-3-ylhydrazine Chemical compound C1CCCCC2=CC(NN)NN=C21 RZZGGSCKACYXBQ-UHFFFAOYSA-N 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000011084 recovery Methods 0.000 abstract description 11
- 230000008878 coupling Effects 0.000 abstract description 4
- 238000010168 coupling process Methods 0.000 abstract description 4
- 238000005859 coupling reaction Methods 0.000 abstract description 4
- 230000035945 sensitivity Effects 0.000 abstract description 2
- 239000000243 solution Substances 0.000 description 9
- 238000010248 power generation Methods 0.000 description 7
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 5
- 229910002092 carbon dioxide Inorganic materials 0.000 description 5
- 229910052739 hydrogen Inorganic materials 0.000 description 5
- 239000001257 hydrogen Substances 0.000 description 5
- 239000000126 substance Substances 0.000 description 4
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000009194 climbing Effects 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 230000029087 digestion Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- XLYOFNOQVPJJNP-ZSJDYOACSA-N heavy water Substances [2H]O[2H] XLYOFNOQVPJJNP-ZSJDYOACSA-N 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- AYEKOFBPNLCAJY-UHFFFAOYSA-O thiamine pyrophosphate Chemical compound CC1=C(CCOP(O)(=O)OP(O)(O)=O)SC=[N+]1CC1=CN=C(C)N=C1N AYEKOFBPNLCAJY-UHFFFAOYSA-O 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- 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/381—Dispersed generators
-
- 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
-
- 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)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Power Engineering (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请公开了一种蓄热罐容量规划方法和电转气参数敏感性分析方法，包括：初始化系统，输入网络拓扑的节点信息、支路的电阻电抗值和风机节点的风电出力数据；确定蓄热罐投资决策经济性目标函数，所述蓄热罐投资决策经济性目标函数包括：蓄热罐建设成本，火电、热电联产、蓄热罐、电转气、燃气锅炉、气源运行成本和弃风成本；确定弃风率目标函数和约束条件；根据所述蓄热罐投资决策经济性目标函数和所述弃风率目标函数，构建蓄热罐投资决策经济性和弃风率多目标模型；通过法线边界交叉法对所述多目标模型求解。考虑电转气热回收的综合能源系统蓄热罐容量规划方法，可以加深电、气、热三网耦合，减少蓄热罐投资与弃风率。
Description
技术领域
本申请涉及电力系统规划技术领域，尤其涉及一种蓄热罐容量规划方法和电转气参数敏 感性分析方法。
背景技术
世界能源正处于大调整、大变革时期，我国也已走上从“扩能保供”向“提质增效”的能源 变革之路。面对化石能源的日益短缺，能源消费结构必须做出改变。从提高各类能源的利用 效率出发，综合能源系统(integrated energy system，IES)成为新的能源发展方向与研究热点。
电转气(power to gas,P2G)实现了能量从电网向天然气网的流动，从而在燃气轮机等耦 合元件促成天然气-电力能量流流向的基础上实现了电-气耦合系统能量双向流动，P2G的弃 风消纳作用与削峰填谷作用是其中关注的热点。鉴于P2G的运行成本在未来短期内无法大幅 降低，对于P2G的研究开始侧重于考虑其运行成本条件下的容量规划与经济调度。然而，在 对电转气的研究中，其强放热的反应特性尚未得到充分关注与利用，这也使得电转气与热网 的互动研究尚未开始
另一方面，中国弃风问题依然严峻，据统计，2017年“三北”地区仍有五省弃风率超10％， 其中甘肃省弃风率高达33％。面对不断增长的风电装机容量和成为隐患的弃风困境，如何最 大化的消纳风电出力，成为运营商和风电企业都非常关心的问题。为了解决热电联产 (combined heat and power,CHP)最小出力难以压缩因而难以给风电提供消纳空间的难题， 蓄热罐(heat accumulator,HA)作为重要的热电解耦元件，其弃风消纳价值受到了广泛研究。
在传统研究分析中，往往在分析蓄热罐对系统的风电消纳及运行成本的影响时，大都将 容量设为给定值或使其足够大，容量选取方式不够合理，且未能考虑电转气设备热量回收对 综合能源系统中蓄热罐容量的影响，导致在计算系统总投资运行成本时，不符合实际应用要 求。
发明内容
本申请提供了一种蓄热罐容量规划方法和电转气参数敏感性分析方法，以提高风电的利 用率以及综合能源系统蓄热罐投资和运行的经济性。
为了解决上述技术问题，本申请实施例公开了如下技术方案：
第一方面，本申请实施例公开了一种蓄热罐容量规划方法，包括：初始化系统，输入网 络拓扑的节点信息、支路的电阻电抗值和风机节点的风电出力数据；
确定蓄热罐投资决策经济性目标函数，所述蓄热罐投资决策经济性目标函数包括：蓄热 罐建设成本，火电、热电联产、蓄热罐、电转气、燃气锅炉、气源运行成本和弃风成本；
确定弃风率目标函数和约束条件；
根据所述蓄热罐投资决策经济性目标函数和所述弃风率目标函数，构建蓄热罐投资决策 经济性和弃风率多目标模型；
通过法线边界交叉法对所述多目标模型求解。
可选的，所述蓄热罐投资决策经济性目标函数为：
cop＝cop,tpp+cop,CHP+cop,HA+cop,p2g+cop,gfb+cop,gs
式中：下标wind、tpp、CHP、p2g、HA、gb、gs分别表示风机机组、火电机组、CHP 机组、P2G机组、HA、燃气锅炉、气源；Fcost为投资规划总成本；cinv、cop、cwind分别为投资成 本、运行成本和弃风成本；n为计算周期数；k为年份；τ为夏/冬季典型日负荷场景数；t为 时段数；κ为折现率；cinv、cop、cwind分别为HA建设、运行、弃风成本；σHA为HA建设的0-1 变量(1表示建设)；CHA、ΔT、ω、KHA分别为HA规划容量、高温供水和低温回水的温差、蓄 热容积率、热损失因数；UHA为HA建设成本；Ufix、Uvar分别为机组运行的固定、可变成本；cper为机组单位运行成本；CE、ε、CM分别为P2G用电电价、生成单位天然气所需CO2系数、CO2价格系数；ζ为单位弃风成本；
分别为节点i风机在场景τ下第t时段的理论 出力、实际出力。
可选的，所述约束条件包括：电力系统约束、天然气网络约束、热力网络约束。
可选的，所述热力网络约束包括：供热约束和蓄热罐运行约束，其中，所述供热约束的 表达式为：
HCHP,a+δhhHA,a+Hgb,a+Hp2g,a≥HL,a
其中：HCHP,a、hHA,a、Hgb,a、Hp2g,a、HL,a分别表示a节点的热电联产、蓄热罐、锅炉、电转气 机组的热出力和热负荷；δh表示蓄热罐的蓄/放热状态；h HA表示蓄热罐的热出力的下限，表示蓄热罐的热出力的上限。
可选的，所述通过法线边界交叉法对所述多目标模型求解包括：
将经济性模型、弃风率模型归一化；
通过改变权重系数，将原目标优化问题转化为一系列单目标优化问题；
通过对一系列单目标问题的寻优，解得Pareto前沿。
可选的，还包括：在MATLAB中采用Yalmip和Gurobi进行建模求解，利用判断是否满足 收敛条件，未满足则返回重新计算，若满足则输出最优解。
第二方面，本申请实施例公开了一种电转气参数敏感性分析方法，包括：
采取投资规划经济性最优为目标函数进行分析；
分析电转气单位运行成本对系统运行策略的影响；
分析电转气放热效率对系统运行策略的影响。
可选的，所述分析电转气单位运行成本对系统运行策略的影响包括：
定义电转气的运行成本系数，确定运行成本系数的上限与仿真步距；
仿真得到蓄热罐规划容量随电转气单位成本变化的曲线图；
在曲线图中标记出多个表现出不同规划策略的运行点；
得到所述运行点处的电转气出力策略、蓄/放热策略，并与不接入电转气的情况对比，得 出接入电转气以及电转气单位运行成本对系统运行策略的影响；
输入电转气单位运行成本并得到对应的各机组出力成本，得到了“蓄热罐规划容量-电转 气单位运行成本图”。
可选的，所述分析电转气放热效率对系统运行策略的影响包括：
定义电转气的放热效率并设置电转气放热效率的区间0～1.2；
按照预定步距输入放热效率的值得到对应的不同P2G放热效率下的蓄热罐出力时刻值， 绘制成曲线，讨论得到电转气放热效率对蓄热罐容量规划的影响。
与现有技术相比，本申请的有益效果为：
本申请提供了一种蓄热罐容量规划方法，包括：初始化系统，输入网络拓扑的节点信息、 支路的电阻电抗值和风机节点的风电出力数据；确定蓄热罐投资决策经济性目标函数，所述 蓄热罐投资决策经济性目标函数包括：蓄热罐建设成本，火电、热电联产、蓄热罐、电转气、 燃气锅炉、气源运行成本和弃风成本；确定弃风率目标函数和约束条件；根据所述蓄热罐投 资决策经济性目标函数和所述弃风率目标函数，构建蓄热罐投资决策经济性和弃风率多目标 模型；通过法线边界交叉法对所述多目标模型求解。建立了同时考虑弃风率与规划经济性的 多目标优化模型，使电转气与蓄热罐协调配合，考虑电转气热回收的综合能源系统蓄热罐容 量规划方法，可以加深电、气、热三网耦合，减少蓄热罐投资与弃风率。同时创新性的考虑 了计及电转气热回收对蓄热罐容量规划与系统运行策略的影响。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制 本申请。
附图说明
为了更清楚地说明本申请的技术方案，下面将对实施例中所需要使用的附图作简单地介 绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动的前提下，还可以根 据这些附图获得其他的附图。
图1为本申请实施例建立的综合能源系统能量流示意图；
图2为本申请实施例提供的蓄热罐容量规划方法流程示意图；
图3为本申请实际算例分析中的9节点综合能源系统结构示意图；
图4为多目标求解后的Pareto前沿；
图5为不同P2G单位成本下的蓄热罐规划容量图；
图6为不同P2G放热效率下的弃风率及蓄热罐规划容量图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中 的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅 是本申请一部分实施例，而不是全部实施例。基于本申请中的实施例，本领域普通技术人员 在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
结合图1所示，为本申请实施例建立的综合能源系统能量流示意图，箭头方向为综合能 源系统能量流动方向。
如图2所示，为本申请实施例提供的蓄热罐容量规划方法流程示意图，本申请提供了一 种综合能源系统蓄热罐容量规划方法，包括：
初始化系统，输入网络拓扑的节点信息、支路的电阻电抗值和风机节点的风电出力数据。
确定蓄热罐投资决策经济性目标函数，所述蓄热罐投资决策经济性目标函数包括：蓄热 罐建设成本，火电、热电联产、蓄热罐、电转气、燃气锅炉、气源运行成本和弃风成本。
所述蓄热罐投资决策经济性目标函数为：
cop＝cop,tpp+cop,CHP+cop,HA+cop,p2g+cop,gfb+cop,gs
式中：下标wind、tpp、CHP、p2g、HA、gb、gs分别表示风机机组、火电机组、CHP 机组、P2G机组、HA、燃气锅炉、气源；Fcost为投资规划总成本；cinv、cop、cwind分别为投资成 本、运行成本和弃风成本；n为计算周期数；k为年份；τ为夏/冬季典型日负荷场景数；t为 时段数；κ为折现率；cinv、cop、cwind分别为HA建设、运行、弃风成本；σHA为HA建设的0-1 变量(1表示建设)；CHA、ΔT、ω、KHA分别为HA规划容量、高温供水和低温回水的温差、蓄 热容积率、热损失因数；UHA为HA建设成本；Ufix、Uvar分别为机组运行的固定、可变成本；cper为机组单位运行成本；CE、ε、CM分别为P2G用电电价、生成单位天然气所需CO2系数、CO2价格系数；ζ为单位弃风成本；分别为节点i风机在场景τ下第t时段的理论 出力、实际出力。
以下各公式中相同的参数代表意义相同，不再一一赘述。
其中：
1)火电成本函数：
式中，tpp表示火电机组；Ufix、Uvar分别表示运行的固定成本和可变成本；P表示火电机 组出力；cper表示机组单位运行成本；τ表示夏/冬季典型日负荷场景。
2)热电联产成本函数：
式中，CHP表示热电联产机组。
3)蓄热罐成本函数：
式中，HA表示蓄热罐；CE表示电价；ρHA表示电泵蓄/放单位热量时的耗电量；hHA表示蓄热罐热出力。
4)电转气成本函数：
式中，α、CM分别表示生成单位天然气所需CO2系数、CO2价格系数。
5)燃气锅炉成本函数：
6)气源成本函数：
7)弃风成本成本函数：
式中，ζ表示单位弃风成本；
分别表示风机在τ场景下的理论发电量和实际 发电量。
确立投资规划经济性目标函数，为各机组运行成本、弃风成本和蓄热罐建设成本之和。 采用等年金法，将总成本分摊到蓄热罐寿命期内的各年年末。
蓄热罐建设成本：
式中，UHA表示蓄热罐建设费用，σHA表示HA建设的0-1变量，CHA、KHA、ΔT、ω分别表示蓄热罐规划容量、热损失因数、高温供水与低温回水的温差、蓄热容积率。
蓄热罐投资规划经济性目标函数：
式中，Fcost表示投资规划总成本，n为计算期，κ为折现率，cinv、cop、cwind分别表示蓄热罐 建设、运行和弃风成本。
确定弃风率目标函数和约束条件。
其中：确立弃风率目标，令某一时刻的弃风量等于该时刻的理论风电发电量与实际发电 量之差。因此规划期每一年的最小化弃风率目标可以表示为：
式中，表示弃风率；分别表示风机在τ场景下的理论发电量和实际发 电量。
确定约束条件，约束条件主要表现为针对系统的约束，包括电力系统约束、天然气网络 约束、热力网络约束。
1.电力系统约束
(a)电功率平衡约束：
其中，σHA为HA的建设状态；ρHA为电泵蓄、放单位热量时的耗电量；Ω(p)为所有与电 力节点p相连的电力节点集合；Pl,pq为节点p和节点q之间的线路潮流。
(b)电力线路潮流约束：
θref＝0
-Plmax≤Pl≤Plmax
其中，θ为节点相角，θref为平衡节点电压相角参考值；x为线路电抗；Plmax为线路潮流 上限。
(c)热电联产出力约束：
其中，
为CHP出力上限，P CHP,p为CHP出力下限。
(d)电转气出力约束：
其中，为P2G出力上限，P p2g,p为P2G出力下限
(e)机组爬坡约束：
其中，PG,up、PG,down分别为CHP和火电机组上行、下行爬坡功率。
2.天然气网络约束：
(a)天然气功率平衡约束。
其中，G为各机组天然气出力或天然气消耗；pl(m)为与燃气节点m连接的所有节点集合； Fpl,mn为节点m和节点n之间的天然气潮流。
(b)管道流量约束。
其中，F为管道流量；sgn(fprem,fpren)为节点m、节点n间的气流方向，fpre为节点气压大小； K为管道传输系数；F、
分别为管道流量下、上限；fpre 、
分别为节点气压下、上限。
(c)气源出力约束：
其中，G gs、
分别为气源出力的下、上限。
3.热力网络约束：
(a)供热约束：
HCHP,a+δhhHA,a+Hgb,a+Hp2g,a≥HL,a
其中：HCHP,a、hHA,a、Hgb,a、Hp2g,a、HL,a分别表示a节点的热电联产、蓄热罐、锅炉、电转气 机组的热出力和热负荷；δh表示蓄热罐的蓄/放热状态；h HA、分别表示蓄热罐的热出力上 下限。
其中：将电转气的总反应过程分为电解制氢和甲烷化两大部分考虑，建立了单位时间内 的电转气甲烷产率与反应放热模型，具体包括步骤：
假设电转气的甲烷化过程中不生成支路产物，因此电转气的总反应可以分解为以下两步 反应：
4H2(g)+CO2(g)——→CH4(g)+2H2O(g)
建立单位时间内电转气的甲烷生产模型(考虑电转气出力的功率单位与化学反应中的体 量单位之间的换算)：
建立单位时间内电转气的反应放热模型：
参数的参考值如下：
甲烷热值
甲烷、氢气密度为
甲烷、 氢气摩尔质量为电转氢气速率
即 
反应热注入热网比重
(b)蓄热罐运行约束。
hHA,t＝(1-KHA)SHA,t-1-SHA,t
SHA,t≤SHA,max,SHA(0)＝SHA(T)
其中，SHA,t表示HA在t时刻的容量状态，单位为MW·h；hHA表示蓄热罐出力； SHA,t≤SHA,max,SHA(0)＝SHA(T)表示周期T(h)结束时蓄热罐的储热量等于其初始的储热量；KHA表 示蓄热罐热损失因数。
将目标函数和约束条件中的非线性模型线性化，比如天然气管道模型、供热约束等。
根据所述蓄热罐投资决策经济性目标函数和所述弃风率目标函数，构建蓄热罐投资决策 经济性和弃风率多目标模型。通过法线边界交叉法对所述多目标模型求解。具体步骤包括：
1)将蓄热罐投资决策经济性模型、弃风率模型归一化；
应用法线边界交叉法(normal boundary intersection，NBI)对投资规划经济性、弃风率目标 归一化后，有：
式中，
分别表示归一化后的规划成本、弃风率目标；
分别表示规 划成本、弃风率目标函数单独优化时的最优解。
2)通过改变权重系数，将原目标优化问题转化为一系列单目标优化问题；
将权重系数表示为β1，则转化后的单目标优化问题可表示为：
maxλ
式中，g(x)＝0、c(x)≤0分别表示模型中的等式约束和不等式约束。
3)通过对一系列单目标问题的寻优，解得Pareto前沿。
在MATLAB中采用Yalmip和Gurobi进行建模求解，利用判断是否满足收敛条件，未满 足则返回重新计算，若满足则输出最优解。
本申请还提供了一种电转气参数敏感性分析方法，具体如下：
1.采取投资规划经济性最优为目标函数进行分析。
2.分析电转气单位运行成本对系统运行策略的影响：
定义电转气的运行成本系数，确定运行成本系数的上限与仿真步距。
取P2G单位运行成本与成本基值的比值为P2G运行成本系数，并以P2G运行成本系数 表示影响运行策略的P2G成本参数，取成本基准值为1000元/MW·h。
仿真得到蓄热罐规划容量随电转气单位成本变化的曲线图。
在曲线图中标记出多个表现出不同规划策略的运行点，本实施例中在曲线图中标记出6 个表现出不同规划策略的运行点。
比较典型运行点处的电转气出力策略、蓄/放热策略，并与不接入电转气的情况对比，得 出是否接入电转气以及电转气单位运行成本对系统运行策略的影响。
输入电转气单位运行成本并得到对应的各机组出力成本，得到了“蓄热罐规划容量-电转 气单位运行成本图”，取P2G单位运行成本与成本基值的比值为P2G运行成本系数，并以P2G 运行成本系数表示影响运行策略的P2G成本参数，取成本基准值为1000元/MW·h。从“蓄 热罐规划容量-电转气单位运行成本图”可分析得到各扰动点形成原因。
3.分析电转气放热效率对系统运行策略的影响：
定义电转气的放热效率并设置电转气放热效率的区间0～1.2；
考虑到P2G的放热效率与P2G的运行效率、反应热注入热网的比例成正相关。因此将 P2G放热效率定义为：
其中，
ηp2g分别表示电转氢效率、甲烷化效率、反应热注入热网比 例和P2G运行效率。
取
为基准值，放热效率上限取
为了使P2G出力最大化，取P2G成本系数cper,p2g＝0。
以0.1为步距输入放热效率的值，得到对应的不同P2G放热效率下的蓄热罐出力时刻值， 绘制成曲线，讨论得到电转气放热效率对蓄热罐容量规划的影响。
以下通过一个实际算例对本申请作进一步详细描述。
算例采用如图3所示的9节点综合能源系统，图3中：黑色线条代表电力系统；灰色线 条代表天然气网络，箭头代表热负荷，TPP代表火电机组，WT代表风机机组，P2G代表 电转气机组，GS代表气源，CHP代表冷热电联产，HA代表蓄热罐，GFB代表燃气锅炉。
其中节点1接入常规火力发电机，节点3接入CHP机组和HA，节点4建设有风电场和P2G，节点6接入燃气锅炉设备，用户电/热负荷分布于5,7,9节点，气源建设于1,2节点。年等金法中计算期n取20年，HA在第一年建设并投入使用。在求解多目标优化结果时，P2G 成本仅取二氧化碳成本，为90美元/t，放热效率取基值0.68。
设置优化目标：1)得到多目标优化的Pareto前沿和折中解；2)研究P2G单位运行成本 对系统运行策略的影响；3)研究P2G放热效率对HA容量优化的影响。
多目标优化结果：解得的Pareto前沿如图4所示。求得能最大程度满足各目标函数的折 中解权重系数为β1＝0.5555，此时优化结果为：系统运行成本为1.8906×108元，弃风率ηWT＝2.12％， 规划容量为93.8864MW。
根据图5，蓄热罐规划容量呈现随cper,p2g的降低而减小的趋势。但是变化曲线并非一条平 滑曲线，这说明蓄热罐的运行策略会随P2G单位运行成本发生改变。
根据图6不同P2G放热效率下的弃风率及蓄热罐规划容量，随着P2G放热效率的提高， 首先，蓄热罐规划容量从不计热回收转变为计及热回收时，弃风率发生跃升，蓄热罐规划容 量发生突降。原因为：(1)计及P2G热回收后，P2G在出力同时可以向蓄热罐蓄热，但由于 蓄热罐存在蓄/放热不能同时进行的限制，因此蓄热罐在放热时P2G将不出力(7～9,16～20时 段)，这使得风电消纳空间暂时性压缩，从而造成了弃风率的跃升；(2)计及P2G热回收后， 蓄热罐将采取负荷低谷时段利用P2G平稳蓄热的策略，从而替代不计P2G热回收时的负荷高 峰时段利用CHP机组增发出力的“临时”蓄热策略，这造成了蓄热罐热出力峰值的削减，使 得蓄热罐容量突降。
为了在电-气-热协同调度系统中规划容量经济且合理的HA，本申请建立了同时考虑弃风 率与规划经济性的多目标优化模型，并创新性的考虑了计及电转气热回收对蓄热罐容量规划 与系统运行策略的影响。首先，通过场景对比研究，初步描述了接入HA和P2G对风电消纳 的促进作用，并着重分析了HA的运行方式。在对多目标优化模型求解后，针对得到的不均 匀Pareto前沿，在规划经济性最佳的单目标决策场景下研究了不同P2G运行参数对系统运行 策略的影响。基于较为典型的算例研究，得出了以下结论：(1)随着P2G单位运行成本降低， 系统将增发P2G出力，从而增大风电消纳空间，此时HA会一改未接入P2G(或不计P2G热 回收)时在负荷高峰时段利用增发CHP机组和燃气锅炉出力的“临时”蓄热策略，转变为在 P2G出力时段利用P2G放热的“稳定”蓄热策略。(2)随着P2G放热效率的提高，HA的蓄热峰值被逐渐削减，放热峰值则逐渐增大，当放热峰值增大至代替蓄热峰值而成为HA容量规划的主导因素时，HA容量将需要提高。
由于以上实施方式均是在其他方式之上引用结合进行说明，不同实施例之间均具有相同 的部分，本说明书中各个实施例之间相同、相似的部分互相参见即可。在此不再详细阐述。
需要说明的是，在本说明书中，术语“包括”、“包含”或者其任何其他变体意在涵盖非 排他性的包含，从而使得包括一系列要素的电路结构、物品或者设备不仅包括那些要素，而 且还包括没有明确列出的其他要素，或者是还包括为这种电路结构、物品或者设备所固有的 要素。
本领域技术人员在考虑说明书及实践这里发明的公开后，将容易想到本申请的其他实施 方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化，这些变型、用途或者适应 性变化遵循本申请的一般性原理并包括本申请未公开的本技术领域中的公知常识或惯用技术 手段。说明书和实施例仅被视为示例性的，本申请的真正范围和精神由权利要求的内容指出。
以上所述的本申请实施方式并不构成对本申请保护范围的限定。
Claims (9)
1.一种蓄热罐容量规划方法，其特征在于，包括：
初始化系统，输入网络拓扑的节点信息、支路的电阻电抗值和风机节点的风电出力数据；
确定蓄热罐投资决策经济性目标函数，所述蓄热罐投资决策经济性目标函数包括：蓄热罐建设成本，火电、热电联产、蓄热罐、电转气、燃气锅炉、气源运行成本和弃风成本；
确定弃风率目标函数和约束条件；
根据所述蓄热罐投资决策经济性目标函数和所述弃风率目标函数，构建蓄热罐投资决策经济性和弃风率多目标模型；
通过法线边界交叉法对所述多目标模型求解。
2.根据权利要求1所述的蓄热罐容量规划方法，其特征在于，所述蓄热罐投资决策经济性目标函数为：
cop＝cop,tpp+cop,CHP+cop,HA+cop,p2g+cop,gfb+cop,gs
式中：下标wind、tpp、CHP、p2g、HA、gb、gs分别表示风机机组、火电机组、CHP机组、P2G机组、HA、燃气锅炉、气源；Fcost为投资规划总成本；cinv、cop、cwind分别为投资成本、运行成本和弃风成本；n为计算周期数；k为年份；τ为夏/冬季典型日负荷场景数；t为时段数；κ为折现率；cinv、cop、cwind分别为HA建设、运行、弃风成本；σHA为HA建设的0-1变量(1表示建设)；CHA、ΔT、ω、KHA分别为HA规划容量、高温供水和低温回水的温差、蓄热容积率、热损失因数；UHA为HA建设成本；Ufix、Uvar分别为机组运行的固定、可变成本；cper为机组单位运行成本；CE、ε、CM分别为P2G用电电价、生成单位天然气所需CO2系数、CO2价格系数；ζ为单位弃风成本；
分别为节点i风机在场景τ下第t时段的理论出力、实际出力。
3.根据权利要求1所述的蓄热罐容量规划方法，其特征在于，所述约束条件包括：电力系统约束、天然气网络约束、热力网络约束。
4.根据权利要求3所述的蓄热罐容量规划方法，其特征在于，所述热力网络约束包括：供热约束和蓄热罐运行约束，其中，所述供热约束的表达式为：
HCHP,a+δhhHA,a+Hgb,a+Hp2g,a≥HL,a
其中：HCHP,a、hHA,a、Hgb,a、Hp2g,a、HL,a分别表示a节点的热电联产、蓄热罐、锅炉、电转气机组的热出力和热负荷；δh表示蓄热罐的蓄/放热状态；h HA表示蓄热罐的热出力的下限，
表示蓄热罐的热出力的上限。
5.根据权利要求1所述的蓄热罐容量规划方法，其特征在于，所述通过法线边界交叉法对所述多目标模型求解包括：
将经济性模型、弃风率模型归一化；
通过改变权重系数，将原目标优化问题转化为一系列单目标优化问题；
通过对一系列单目标问题的寻优，解得Pareto前沿。
6.根据权利要求1所述的蓄热罐容量规划方法，其特征在于，还包括：在MATLAB中采用Yalmip和Gurobi进行建模求解，利用判断是否满足收敛条件，未满足则返回重新计算，若满足则输出最优解。
7.一种电转气参数敏感性分析方法，其特征在于，包括：
采取投资规划经济性最优为目标函数进行分析；
分析电转气单位运行成本对系统运行策略的影响；
分析电转气放热效率对系统运行策略的影响。
8.根据权利要求7所述的电转气参数敏感性分析方法，其特征在于，所述分析电转气单位运行成本对系统运行策略的影响包括：
定义电转气的运行成本系数，确定运行成本系数的上限与仿真步距；
仿真得到蓄热罐规划容量随电转气单位成本变化的曲线图；
在曲线图中标记出多个表现出不同规划策略的运行点；
得到所述运行点处的电转气出力策略、蓄/放热策略，并与不接入电转气的情况对比，得出是否接入电转气以及电转气单位运行成本对系统运行策略的影响；
输入电转气单位运行成本并得到对应的各机组出力成本，得到了“蓄热罐规划容量-电转气单位运行成本图”。
9.根据权利要求7所述的电转气参数敏感性分析方法，其特征在于，所述分析电转气放热效率对系统运行策略的影响包括：
定义电转气的放热效率并设置电转气放热效率的区间0～1.2；
按照预定步距输入放热效率的值得到对应的不同P2G放热效率下的蓄热罐出力时刻值，绘制成曲线，讨论得到电转气放热效率对蓄热罐容量规划的影响。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910991901.4A CN110729720A (zh) | 2019-10-18 | 2019-10-18 | 一种蓄热罐容量规划方法和电转气参数敏感性分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910991901.4A CN110729720A (zh) | 2019-10-18 | 2019-10-18 | 一种蓄热罐容量规划方法和电转气参数敏感性分析方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110729720A true CN110729720A (zh) | 2020-01-24 |
Family
ID=69220313
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910991901.4A Pending CN110729720A (zh) | 2019-10-18 | 2019-10-18 | 一种蓄热罐容量规划方法和电转气参数敏感性分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110729720A (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111445350A (zh) * | 2020-03-30 | 2020-07-24 | 华能国际电力股份有限公司上安电厂 | 应用于供热管网的蓄热罐选型计算方法 |
| CN112035984A (zh) * | 2020-08-10 | 2020-12-04 | 浙江华云电力工程设计咨询有限公司 | 一种电-气-储区域综合能源系统协同规划方法 |
| CN112163780A (zh) * | 2020-10-14 | 2021-01-01 | 国网能源研究院有限公司 | 一种风光互补制-储氢容量规划方法及系统 |
| CN112288216A (zh) * | 2020-09-04 | 2021-01-29 | 广西电网有限责任公司电力科学研究院 | 一种基于合作博弈的电转气装置容量规划方法及系统 |
| CN114417625A (zh) * | 2022-01-24 | 2022-04-29 | 太原理工大学 | 一种考虑风气互补特性的季节储能解决方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110097294A (zh) * | 2019-05-16 | 2019-08-06 | 长沙理工大学 | 一种考虑虚拟电厂经济运行的蓄热罐定容配置决策方法 |
-
2019
- 2019-10-18 CN CN201910991901.4A patent/CN110729720A/zh active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110097294A (zh) * | 2019-05-16 | 2019-08-06 | 长沙理工大学 | 一种考虑虚拟电厂经济运行的蓄热罐定容配置决策方法 |
Non-Patent Citations (1)
| Title |
|---|
| 李东森等: "计及电转气热回收的综合能源系统蓄热罐容量规划与运行策略" * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111445350A (zh) * | 2020-03-30 | 2020-07-24 | 华能国际电力股份有限公司上安电厂 | 应用于供热管网的蓄热罐选型计算方法 |
| CN112035984A (zh) * | 2020-08-10 | 2020-12-04 | 浙江华云电力工程设计咨询有限公司 | 一种电-气-储区域综合能源系统协同规划方法 |
| CN112035984B (zh) * | 2020-08-10 | 2024-04-09 | 浙江华云电力工程设计咨询有限公司 | 一种电-气-储区域综合能源系统协同规划方法 |
| CN112288216A (zh) * | 2020-09-04 | 2021-01-29 | 广西电网有限责任公司电力科学研究院 | 一种基于合作博弈的电转气装置容量规划方法及系统 |
| CN112163780A (zh) * | 2020-10-14 | 2021-01-01 | 国网能源研究院有限公司 | 一种风光互补制-储氢容量规划方法及系统 |
| CN114417625A (zh) * | 2022-01-24 | 2022-04-29 | 太原理工大学 | 一种考虑风气互补特性的季节储能解决方法 |
| CN114417625B (zh) * | 2022-01-24 | 2024-03-26 | 太原理工大学 | 一种考虑风气互补特性的季节储能解决方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110729720A (zh) | 一种蓄热罐容量规划方法和电转气参数敏感性分析方法 | |
| CN106208157B (zh) | 基于电转气的电-气互联综合能源系统削峰填谷方法 | |
| JP7261507B2 (ja) | 電気ヒートポンプ－熱電併給システムを最適化する調整方法及びシステム | |
| CN108832656B (zh) | 基于电转气和可再生能源利用微能源网多目标规划方法 | |
| Zhang et al. | Optimization analysis of a novel combined heating and power system based on biomass partial gasification and ground source heat pump | |
| CN106058942B (zh) | 考虑风电不确定性的含电转气和cchp的能量枢纽优化模型 | |
| CN110544025B (zh) | 电转气-储气罐组合的气电综合能源系统优化调度方法 | |
| CN111859683B (zh) | 基于动态能量集线器的园区综合能源系统优化配置方法 | |
| CN112347607A (zh) | 一种基于凸松弛的热电联合调度方法 | |
| CN113850474A (zh) | 一种热电氢多能流综合能源系统及其优化调度方法 | |
| CN105244870A (zh) | 一种快速计算电网风电场弃风率及机组发电量的方法 | |
| CN112531687A (zh) | 含热电联合机组的综合能源系统周前机组组合优化方法 | |
| CN116468215A (zh) | 计及源荷不确定性的综合能源系统调度方法及装置 | |
| CN115659651A (zh) | 一种考虑多种灵活性资源的综合能源协同优化调度方法 | |
| CN116681228A (zh) | 一种考虑场景不确定性和碳排放的综合能源系统设备容量优化配置方法 | |
| CN114861546A (zh) | 一种基于燃氢燃气轮机的电热氢多能互补调度方法及系统 | |
| CN113011655B (zh) | 一种基于两阶段随机规划的社区多能源系统的混合储能系统规划方法 | |
| CN106100002A (zh) | 一种交直流混合微网的优化运行方法 | |
| CN111724026B (zh) | 一种多能源网络与配水网耦合运行的优化方法 | |
| CN110457745B (zh) | 一种具有紧密耦合关系的水-能联系网络模型的建立方法 | |
| CN116596123A (zh) | 一种园区综合能源系统低碳优化调度方法 | |
| CN109241676B (zh) | 综合能源系统中热网及建筑物时间分辨率的选择方法 | |
| CN112767189B (zh) | 一种计及风电并网下电网灵活性评估方法 | |
| CN111428940B (zh) | 一种基于储热装置降低多能源系统运行成本的优化方法 | |
| CN113794242A (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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200124 |
|
| RJ01 | Rejection of invention patent application after publication |