CN115455709A - 计及碳捕集设备加装的低碳综合能源系统仿真及配置方法 - Google Patents
计及碳捕集设备加装的低碳综合能源系统仿真及配置方法 Download PDFInfo
- Publication number
- CN115455709A CN115455709A CN202211140985.9A CN202211140985A CN115455709A CN 115455709 A CN115455709 A CN 115455709A CN 202211140985 A CN202211140985 A CN 202211140985A CN 115455709 A CN115455709 A CN 115455709A
- Authority
- CN
- China
- Prior art keywords
- carbon
- model
- power
- heat
- equipment
- 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
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 194
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 178
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004088 simulation Methods 0.000 title claims abstract description 16
- 239000007789 gas Substances 0.000 claims abstract description 67
- 239000002918 waste heat Substances 0.000 claims abstract description 27
- 238000009434 installation Methods 0.000 claims abstract description 17
- 238000004146 energy storage Methods 0.000 claims abstract description 16
- 230000008878 coupling Effects 0.000 claims abstract description 8
- 238000010168 coupling process Methods 0.000 claims abstract description 8
- 238000005859 coupling reaction Methods 0.000 claims abstract description 8
- 238000009826 distribution Methods 0.000 claims abstract description 8
- 238000005485 electric heating Methods 0.000 claims abstract description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 55
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 38
- 239000003345 natural gas Substances 0.000 claims description 27
- 238000003860 storage Methods 0.000 claims description 21
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 16
- 238000010521 absorption reaction Methods 0.000 claims description 10
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 8
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 8
- 239000001569 carbon dioxide Substances 0.000 claims description 8
- 229910052739 hydrogen Inorganic materials 0.000 claims description 8
- 239000001257 hydrogen Substances 0.000 claims description 8
- 230000008929 regeneration Effects 0.000 claims description 8
- 238000011069 regeneration method Methods 0.000 claims description 8
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 7
- 239000007788 liquid Substances 0.000 claims description 7
- 239000003546 flue gas Substances 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 6
- 238000010248 power generation Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000009194 climbing Effects 0.000 claims description 4
- 238000007906 compression Methods 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 4
- 238000005094 computer simulation Methods 0.000 claims description 4
- 238000011084 recovery Methods 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 229910002091 carbon monoxide Inorganic materials 0.000 claims description 2
- 238000005868 electrolysis reaction Methods 0.000 claims description 2
- 239000002803 fossil fuel Substances 0.000 claims description 2
- 230000005251 gamma ray Effects 0.000 claims description 2
- 125000004435 hydrogen atom Chemical group [H]* 0.000 claims description 2
- 238000005984 hydrogenation reaction Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 230000000295 complement effect Effects 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 abstract description 2
- 230000009102 absorption Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 238000002485 combustion reaction Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 230000009919 sequestration Effects 0.000 description 2
- 239000003245 coal Substances 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012983 electrochemical energy storage Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000005338 heat storage Methods 0.000 description 1
- 150000002431 hydrogen Chemical class 0.000 description 1
- 238000012946 outsourcing Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Primary Health Care (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种计及碳捕集设备加装的低碳综合能源系统仿真及配置方法,通过构建包含碳捕集电厂模型、P2G设备模型、光伏模型、燃气轮机模型、余热锅炉模型、燃气锅炉模型、电锅炉模型、地源热泵模型、储能设备模型、电力网络模型、热力网络模型的考虑碳捕集加装的综合能源系统耦合模型以及考虑碳捕集加装的城市区域综合能源系统经济调度模型,然后选择修正后的IEEE33节点配电网与六节点热力网构成的区域电热综合能源系统,基于matlab平台,通过yalmip和gurobi求解器进行仿真求解,最后根据求解结果分析综合能源系统的配置阶段并得到最优配置方案,实现电‑热多能互补综合能源系统的运行和配置的优化。
Description
技术领域
本发明涉及的是一种综合能源领域的技术,具体是一种计及碳捕集设备加装的低碳综合能源系统仿真及配置方法。
背景技术
城市综合能源系统因其能源利用效率高和可再生能源消纳比例高的优势,成为低碳化发展的重要方向,不同于传统的电力调度,城市能源系统是城市供热、供冷、供电和满足其他需求的能源系统的集合,现阶段继续一种能够应对综合能源系统参与碳排放的优化配置技术。
发明内容
本发明针对现有技术存在的上述不足,提出一种计及碳捕集设备加装的低碳综合能源系统仿真及配置方法,综合考虑火电机组的碳捕集结构和碳交易成本,实现电-热多能互补综合能源系统的运行和配置的优化。
本发明是通过以下技术方案实现的:
本发明涉及一种计及碳捕集设备加装的低碳综合能源系统仿真及配置方法,通过构建包含碳捕集电厂模型、P2G设备模型、光伏模型、燃气轮机模型、余热锅炉模型、燃气锅炉模型、电锅炉模型、地源热泵模型、储能设备模型、电力网络模型、热力网络模型的考虑碳捕集加装的综合能源系统耦合模型以及考虑碳捕集加装的城市区域综合能源系统经济调度模型,然后选择修正后的IEEE33节点配电网与六节点热力网构成的区域电热综合能源系统,基于matlab平台,通过yalmip和gurobi求解器进行仿真求解,最后根据求解结果分析综合能源系统的配置阶段并得到最优配置方案。
所述的考虑碳捕集加装的综合能源系统经济调度模型包括:火电机组运行成本、燃气机组运行成本、二氧化碳相关成本和新能源弃置成本目标函数以及电网运行约热网运行约束和设备约束。
本发明涉及一种实现上述方法的系统,包括:综合能源系统建模单元、综合能源系统调度单元以及综合能源系统配置单元,其中:综合能源系统建模单元根据所选的设备和网络信息,进行数学建模,得到设备和网络模型,综合能源系统调度单元根据综合能源系统的运行目标和约束条件,进行数学建模,得到调度模型,综合能源系统配置单元根据设备和网络模型和调度模型,基于matlab平台进行处理,通过Yalmip工具箱调用Gurobi求解器进行求解,得到综合能源系统配置方案。
技术效果
本发明通过分析碳捕集设备物理模型综合能源系统配置方法之间的相关性;通过热力网络的精确建模,使综合能源系统仿真更加精确;通过构建综合能源系统的边际碳价与碳捕集加装的边际碳价的拟合函数,自适应的根据碳市场的发展水平判断系统的配置方案,从而指导工程人员提升综合能源系统的运行和规划水平。
附图说明
图1为本发明流程图;
图2为本发明综合能源系统耦合模型结构示意图;
图3为本发明碳捕集系统结构示意图。
图4为综合能源系统热网结构示意图。
图5为综合能源系统热网建模对调度的影响对比示意图。
图6为综合能源系统多阶段配置方案对调度的影响对比示意图。
具体实施方式
如图1所示,为本实施例涉及一种计及碳捕集设备加装的低碳综合能源系统与配仿真置方法,包括以下步骤:
步骤1:构建包含碳捕集电厂模型、P2G设备模型、光伏模型、燃气轮机模型、余热锅炉模型、燃气锅炉模型、电锅炉模型、地源热泵模型、储能设备模型、电力网络模型、热力网络模型的考虑碳捕集加装的综合能源系统耦合模型,具体包括:
步骤1.1:建立碳捕集电厂模型,碳捕集电厂是在传统煤电厂的基础上加装碳捕集设备构成以降低系统碳排放。如图3所示,碳捕集电厂模型配置有烟气旁路系统、贫液富液存储单元、再生塔和压缩机,碳捕集电厂模型为Phd,r,t=Phd,t-Ploss,t,其中:Phd,r,t为火电机组加装碳捕集设备后实际输出,Phd,t为火电机组t时刻实际发电功率,Ploss,t为t时刻碳捕集系统电能损耗功率,其运行环节包括:烟气吸收环节和再生和压缩环节,因此碳捕集系统电能损耗功率其中:Mr,t=Mhd,t-Min,t,a1、a2、a3为耗电系数,Mr,t为火电机组实际碳排放量,火电机组在烟气吸收环节t时刻碳排放强度Mhd,t=γhdPhd,t,Min,t=αyqMhd,t,Min,t为t时刻进入富液吸收塔的量;γhd为火电机组碳排放比例,αyq为碳捕集设备吸收比例;贫富液存储器再生和压缩环节t时刻中的碳含量ML,t=ML,t-1+Min,t-Mout,t,Mout,t为t时刻从贫富液存储器进入再生塔的碳量,为t时刻系统压缩的碳量;该碳捕集电厂模型的输入为火电机组发电功率Phd,t,输出为碳捕集电厂实际电功率输出Phd,r,t,碳捕集量和实际碳排放量Mr,t。
步骤1.2:建立P2G设备运行模型,该P2G设备通过消耗二氧化碳降低系统碳排放,包括电解水制氢气和甲烷化二氧化碳加氢气制天然气,具体模型为:电解水制氢气体积 甲烷化二氧化碳加氢气制天然气体积其中:PP2H为电解水耗电量,ηp2h为能源转换效率,为氢气的热值,为氢气产量,为天然气产量,ηh2g为甲烷化的效率取0.75,为二氧化碳消耗量;该P2G设备运行模型的输入为电解水耗电量PP2H,输出为氢气产量
步骤1.3:建立光伏模型,系统中引入分布式光伏代替原有电能供给,可以降低系统碳排放。光伏出力水平由预测得到,光伏模型为Ppv,t=ηpv,tPpv,p,t,其中:Ppv,t为光伏实际上网功率,ηpv,t为光伏消纳率,Ppv,p,t为光伏实际发电量;该光伏模型的输入为光伏实际发电量Ppv,p,t,输出为光伏实际上网功率Ppv,t。
步骤1.4:建立燃气轮机模型,燃气轮机通过消耗天然气生产电能,同时会产生余热,可以通过余热锅炉输出热功率。这个过程中会因为天然气的燃烧产生碳排放,燃气轮机模型为PGT,t=ηGTLNGVGT,t,其中:ηGT为燃气轮机发电效率,LNG为天然气热值,VGT,t为燃气轮机t时刻每小时消耗的天然气量,PGT,t为燃气轮机t时刻的发电功率,ηr为余热回收效率,QGB,t为t时刻燃气轮机排气余热回收量;该燃气轮机模型的输入为燃气轮机消耗的天然气量VGT,t,输出为燃气轮机时刻的发电功率PGT,t。
步骤1.5:建立余热锅炉模型,余热锅炉收集燃气轮机产生的余热,输出热功率,余热锅炉模型为QWH,t=ηWHQGB,t,其中:ηWH为余热锅炉集热效率,QGB,t为吸收的余热,QWH,t为余热锅炉输出功率;该余热锅炉模型的输入为吸收的余热QGB,t,输出为余热锅炉功率QWH,t。
步骤1.6:建立燃气锅炉模型,燃气锅炉通过天然气产热,产生的热量与锅炉效率、燃料量有关,在天然气燃烧的过程中产生碳排放。燃气锅炉模型为QGB,t=VGB,tηGBLNG,其中:QGB,t为t时刻燃气锅炉产生的热量,ηGB为燃气轮机发电效率,LNG为天然气热值,VGB,t为燃气轮机t时刻内每小时消耗的天然气量;该燃气锅炉模型的输入为燃气轮机消耗的天然气量VGB,t,输出为燃气锅炉产生的热量QGB,t。
步骤1.7:建立电锅炉模型,该电锅炉将电能转化为热能,为系统提供可调节的热负荷,调节范围较大,但能源转化效率较低。电锅炉模型为Qgd,t=ηgdPgd,t,其中:Qgd,t为t时刻电锅炉产生的热量,ηgd为电锅炉产热效率,Pgd,t为t时刻电锅炉耗电量;该电锅炉模型的输入为电锅炉耗电量Pgd,t,输出为电锅炉产生的热量Qgd,t。
步骤1.8:建立地源热泵模型,地源热泵以岩土体、地下水或地表水为低温热源,可以为用户提供热能服务,产热过程需要消耗少量电能,地源热泵一般只能提供固定热源,作为系统供热的基载。地源热泵模型为Qrb,t=ηrbPrb,t,其中:Qrb,t为t时刻地源热泵产生的热量,ηrb为地源热泵产热效率,Prb,t为t时刻地源热泵耗电量;该地源热泵模型的输入为地源热泵耗电量Prb,t,输出为地源热泵产生的热量Qrb,t。
步骤1.9:建立储能设备模型,综合能源系统中包含电储能和热储能,分别采用电化学储能模型和储热罐模型。同时存在碳存储和气存储,损耗较低可忽略,其中:碳存储主要服务于P2G设备和碳捕集后碳封存,封存不设上限。气存储主要服务于天然气机组和P2G设备,气存储主要用于保证每日的天然气供应平衡。
步骤1.10:建立电力网络模型,针对配网中闭环设计、开环运行的现状,选用辐射状配电网线性化潮流模型描述综合能源系统中的电力网络,同时忽略支路损耗。电力网络模型由简化后的支路潮流等式构成,包括Phj=∑i→hPih-Ph,Qhj=∑i→hQih-Qh, 其中:Phj和Qhj分别为从节点h到节点j流过的有功和无功功率;Ph和Qh分别为流向节点h负荷的有功和无功功率;Rih和Xih分别为支路ih的电阻和电抗;yh和Vi分别为节点h和节点i的电压幅值。
步骤1.11:根据热水网络设计热力网络,如图4所示,该热力网络包括一次热网和二次热网,存在动态延时特性和储能特性,采用质调节的方法建模,即不改变热网中热水的质量流率,仅调节水温,同时计及回流管道热损。一次热网由具有供热设备的换热站提供热量,热用户通过二次热网获得热量。
步骤1.12:联立步骤1.1-1.9中设备模型方程和步骤1.10-1.11的网络模型方程,共同构成综合能源系统耦合模型,其中:燃气轮机、燃气锅炉、电锅炉和热泵提供热负荷;火电机组、光伏、燃气轮机和外部电网提供电负荷;储能和P2G设备用于提高能源消纳效率,改变负荷曲线;火电机组产生的碳排放经过碳捕集设备回收利用,从而降低碳排放量,燃气机组产生的碳排放难以捕集,需要和火电机组捕集后残存的碳排放一同参与碳市场交易。
步骤2:以步骤1.12中综合能源系统耦合模型为基础,构建考虑碳捕集加装的城市区域综合能源系统经济调度模型,该模型包括目标函数和约束条件。
所述的目标函数包括:火电机组运行成本、燃气机组运行成本、碳相关成本、新能源弃置成本的目标函数,具体包括:
1)火电机组运行成本目标函数,其中:T为调度周期,Ce0,,t为启动成本,Cet,t为碳捕集设备加装成本,折算到每日运行,f(Pe,t)为发电成本函数,Pe,t是发电功率,ae、be和ce为火电机组的成本系数。
2)燃气机组运行成本目标函数,Cg=∑t∈T pg,tVg,t,其中:pg,t为天然气价格;Vg,t为燃气机组天然气耗气量之和。
3)二氧化碳相关成本目标函数,包括系统化石燃料机组的碳交易成本、P2G设备的购碳成本和储碳设备的储碳成本Ccur=∑t∈Tpcur,tδtPre,t,其中:pc1,t为机组碳交易价格,μe为火电机组碳排放强度,μg为燃气机组碳排放强度,pc2,t为P2G设备购碳价格,Cp2g,t为P2G购碳量,ps,t碳存储价格,为碳储存量。
4)新能源弃置成本目标函数,Ccur=∑t∈Tpcur,tδtPre,t,其中:pcur,t为弃置惩罚系数,δt为弃置率,Pre,t为可用有功出力。
所述的约束条件包括:电网约束和热网约束,具体包括:
i)电网功率平衡约束电网线路传输约束和电网电压相角约束其中:Pi,t为j节点的所连机组出力,AG为节点所连机组集合,Phj,t为支路hj 的潮流,AF和AE为节点j为起点和终点的线路集合,和分别为潮流的上下限,Dj,t为节点j的电力负荷需求,θj,t为节点j的电压相角。
ii)热网约束:针对电力网络潮流和热力网络潮流惯性不同的问题,以电力网络潮流计算为基础,在热力网络中考虑延时特性,即通过的热能输运准动态过程温度约束 描述不同时段下进行热力网络的热水,此外,热网运行约束包含还热交换约束 汇流节点温度方程热源温度约束 c1、c2代表t时刻输出的热水质块中的权重系数,和为为质块的初始温度,Tam为环境温度,c3、c4为热水在管道中温度损失的系数,与管道长度和管道流速有关;和分别为节点ν处的热源设备的热出力和用户热负荷;mH,v,t和mL,v,t为热源和热负荷处的热水质量;和为热源和热负荷的供水温度和回水温度;S(v)+和S(v)-分别为出水管道和入水管道与节点ν相连的管道集合,为管道κ的出水温度,为节点v的流出温度。qb和qκ分别为管道b和κ的热水的质量流率。
iii)设备约束,由机组约束和存储设备约束构成,其中:
机组约束包括机组出力约束机组爬坡约束Pn,i,d≤Pn,i,t-Pn,i,t-1≤Pn,i,u,其中:Pn,i,t为t时刻第i个节点第n个设备的出力,和为该设备出力上下限;Pn,i,d和Pn,i,u为该设备爬坡上下限。
存储设备约束包括设备状态方程Si,t=Si,t-1+Si,in,t-Si,out,t总量约束Si,min≤Si,t≤Si,max,能量存储约束0≤Si,in,t≤Si,in,max(1-λi,t),能量释放约束0≤Si,out,t≤Si,out,maxλi,t,其中:Si,t为t时刻第i个设备能量储存总量;Si,in,t为该设备t时刻存储功率;Si,out,t为该设备t时刻释放功率;λi,t为01变量为设备状态。
步骤3:选择修正后的IEEE33节点配电网与六节点热力网构成的区域电热综合能源系统,基于matlab平台,通过yalmip和gurobi求解器仿真求解步骤2得到的城市区域综合能源系统经济调度模型,具体包括:
步骤3.1:设计综合能源系统运行场景,其中:设备包括大型火力机组、电转气设备、燃气轮机、余热锅炉、燃气锅炉、光伏、电锅炉、地源热泵、电热储能,同时包含电热负荷用户,该运行场景的具体设备参数包括:
表1
步骤3.2:分析热力网络建模影响:由于热力网络模型包括只保证系统热负荷与供给平衡的线性化模型和考虑热水网络结构的热网模型,因此不同建模下综合能源系统运行的结果不完全相同,其中碳排放量的区别不容忽略,具体包括:两种方案下系统运行结果如图5所示,通过观察热力相关设备运行曲线不难发现,普通线性化模型中燃气锅炉产热较少,热网模型中整体供热水平更高,说明考虑实际管道中热水流动后系统热负荷不完全等于系统用户热负荷,管道中的热水会存在时延等现象,降低系统的热负荷需求,从而降低系统碳排放,因此两种建模方法存在差异,热网模型更符合实际情况,应在综合能源系统建模中考虑。
表2
步骤4:根据求解结果分析三个综合能源系统的配置阶段并得到最优配置方案,其中:阶段1:系统不需要考虑碳交易成本,尚未加装碳捕集设备,综合能源系统正常运行;阶段2:系统内存在碳交易成本,尚未加装碳捕集设备,系统目标考虑碳交易成本;阶段3:系统内存在碳交易成本,在原有基础上加装碳捕集,系统运行考虑加装成本,具体包括:
步骤4.1:配置方案分析,三个阶段的运行结果如表3和图6所示,其中阶段1和阶段2系统设备运行没有明显区别,可见碳交易机制引入后,因为碳成本在系统中的决策占比不高,对系统运行影响很小,从表中的成本也能看出,各项发电机组的成本变化很小,碳交易引入后系统日碳排放量也只降低36.129t,说明仅仅在综合能源系统的运行中考虑碳交易机制不够充分,不能实现低碳目标。而从阶段2过度到阶段3碳排放水平明显降低,说明碳捕集机组的引入有足够的价值。从成本上看,外购电上升,天然气基本不变,火电成本上升但实际发电量下降,这是因为碳捕集设备引入后会增加火电系统的自用电,从而导致火电出力出现负值。阶段3设定的碳价为0.2千元/t,此时加装碳捕集后可以降低用能成本。
表3
步骤4.2:碳交易价格分析:对比碳价变化时步骤3.3中的阶段2和阶段3系统的运行成本,从而说明碳捕集设备的合适加装时机;然后分析碳捕集设备加装成本变动时,对应的可提高系统收益的边际碳价变化。
步骤4.3:碳价变化的影响分析。随着碳交易成本的上升,通过观察阶段2和阶段3差价的变化,当碳价为0.198千元/t时,差价约等于0,可以认为碳价超过0.198千元/t后,综合能源系统加装碳捕集设备可以提高系统收益,在此之后阶段3带来的碳交易成本下降能够覆盖碳捕集设备加装的成本,即此时边际碳价为0.198千元/t。这说明综合能源系统参与碳交易市场后,需要做好价格的预估,根据自身的碳捕集成本确定配置方案,同时碳交易市场也要需要合理调控碳价,才能助力降碳目标的实现。
表4
步骤4.4:分析边际碳价和碳捕集加装的关系。碳捕集设备加装成本带来的影响如表5所示,在以火电为核心的综合能源系统中,边际碳价同碳捕集设备的成本关系密切,由表中数据可以判断二者关系为线性曲线。对于算例中的综合能源系统而言,边际碳价pc和碳捕集加装wccs的拟合函数为wccs=1933.431(pc-0.172),因此在实际综合能源系统建设中,可以根据火电厂碳捕集设备加装成本和碳交易市场的碳价做出决策,当实际碳价确定后,带入拟合函数计算理想碳捕集成本,实际的碳捕集成本更低时即可加装。同时考虑到碳捕集设备的加装对环境具有正面作用,可将环境收益折算后,在加装成本中减去对应数值后判断是否加装。
表5
通过步骤4中的三个阶段,完成配置方案的分析,帮助综合能源系统碳捕集设备安装提高效率。
与现有技术相比,本方法通过热力网络模型,综合能源系统运行仿真更加精准;计算边际碳价和碳捕集装置函数关系后,碳捕集的工业安装应对碳市场的发展时效率更高。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。
Claims (9)
1.一种计及碳捕集设备加装的低碳综合能源系统仿真及配置方法,其特征在于,通过构建包含碳捕集电厂模型、P2G设备模型、光伏模型、燃气轮机模型、余热锅炉模型、燃气锅炉模型、电锅炉模型、地源热泵模型、储能设备模型、电力网络模型、热力网络模型的考虑碳捕集加装的综合能源系统耦合模型以及考虑碳捕集加装的城市区域综合能源系统经济调度模型,然后选择修正后的IEEE33节点配电网与六节点热力网构成的区域电热综合能源系统,基于matlab平台,通过yalmip和gurobi求解器进行仿真求解,最后根据求解结果分析综合能源系统的配置阶段并得到最优配置方案;
所述的考虑碳捕集加装的综合能源系统经济调度模型包括:火电机组运行成本、燃气机组运行成本、二氧化碳相关成本和新能源弃置成本目标函数以及电网运行约热网运行约束和设备约束。
2.根据权利要求1所述的计及碳捕集设备加装的低碳综合能源系统仿真及配置方法,其特征是,所述的碳捕集电厂模型配置有烟气旁路系统、贫液富液存储单元、再生塔和压缩机,碳捕集电厂模型为Phd,r,t=Phd,t-Ploss,t,其中:Phd,r,t为火电机组加装碳捕集设备后实际输出,Phd,t为火电机组t时刻实际发电功率,Ploss,t为t时刻碳捕集系统电能损耗功率,其运行环节包括:烟气吸收环节和再生和压缩环节,因此碳捕集系统电能损耗功率 其中:Mr,t=Mhd,t-Min,t,a1、a2、a3为耗电系数,Mr,t为火电机组实际碳排放量,火电机组在烟气吸收环节t时刻碳排放强度Mhd,t=γhdPhd,t,Min,t=αyqMhd,t,Min,t为t时刻进入富液吸收塔的量;γhd为火电机组碳排放比例,αyq为碳捕集设备吸收比例;贫富液存储器再生和压缩环节t时刻中的碳含量Ml,t=ML,t-1+Min,t-Mout,t,Mout,t为t时刻从贫富液存储器进入再生塔的碳量,为t时刻系统压缩的碳量;该碳捕集电厂模型的输入为火电机组发电功率Phd,t,输出为碳捕集电厂实际电功率输出Phd,r,t,碳捕集量和实际碳排放量Mr,t;
所述的P2G设备运行模型具体为:电解水制氢气体积甲烷化二氧化碳加氢气制天然气体积其中:PP2H为电解水耗电量,ηp2h为能源转换效率,为氢气的热值,为氢气产量,为天然气产量,ηh2g为甲烷化的效率取0.75,为二氧化碳消耗量;该P2G设备运行模型的输入为电解水耗电量PP2H,输出为氢气产量
所述的光伏模型为Ppv,t=ηpv,tPpv,p,t,其中:Ppv,t为光伏实际上网功率,ηpv,t为光伏消纳率,Ppv,p,t为光伏实际发电量;该光伏模型的输入为光伏实际发电量Ppv,p,t,输出为光伏实际上网功率Ppv,t。
3.根据权利要求1所述的计及碳捕集设备加装的低碳综合能源系统仿真及配置方法,其特征是,所述的燃气轮机模型为PGT,t=ηGTLNGVGT,t,其中:ηGT为燃气轮机发电效率,LNG为天然气热值,VGT,t为燃气轮机t时刻每小时消耗的天然气量,PGT,t为燃气轮机t时刻的发电功率,ηr为余热回收效率,QGB,t为t时刻燃气轮机排气余热回收量;该燃气轮机模型的输入为燃气轮机消耗的天然气量VGT,t,输出为燃气轮机时刻的发电功率PGT,t;
所述的余热锅炉模型为QWH,t=ηWHQGB,t,其中:ηWH为余热锅炉集热效率,QGB,t为吸收的余热,QWH,t为余热锅炉输出功率;该余热锅炉模型的输入为吸收的余热QGB,t,输出为余热锅炉功率QWH,t;
所述的燃气锅炉模型为QGB,t=VGB,tηGBLNG,其中:QGB,t为t时刻燃气锅炉产生的热量,ηGB为燃气轮机发电效率,LNG为天然气热值,VGB,t为燃气轮机t时刻内每小时消耗的天然气量;该燃气锅炉模型的输入为燃气轮机消耗的天然气量VGB,t,输出为燃气锅炉产生的热量QGB,t;
所述的电锅炉模型为Qgd,t=ηgdPgd,t,其中:Qgd,t为t时刻电锅炉产生的热量,ηgd为电锅炉产热效率,Pgd,t为t时刻电锅炉耗电量;该电锅炉模型的输入为电锅炉耗电量Pgd,t,输出为电锅炉产生的热量Qgd,t。
4.根据权利要求1所述的计及碳捕集设备加装的低碳综合能源系统仿真及配置方法,其特征是,所述的地源热泵模型为Qrb,t=ηrbPrb,t,其中:Qrb,t为t时刻地源热泵产生的热量,ηrb为地源热泵产热效率,Prb,t为t时刻地源热泵耗电量;该地源热泵模型的输入为地源热泵耗电量Prb,t,输出为地源热泵产生的热量Qrb,t;
所述的电力网络模型,针对配网中闭环设计、开环运行的现状,选用辐射状配电网线性化潮流模型描述综合能源系统中的电力网络,同时忽略支路损耗;简化后的支路潮流等式为Phj=∑i→hPih-Ph,Qhj=∑i→hQih-Qh,其中:Phj和Qhj分别为从节点h到节点j流过的有功和无功功率;Ph和Qh分别为流向节点h负荷的有功和无功功率;Rih和Xih分别为支路ih的电阻和电抗;Vh和Vi分别为节点h和节点i的电压幅值;
所述的热水网络设计热力网络包括一次热网和二次热网,存在动态延时特性和储能特性,采用质调节的方法建模,即不改变热网中热水的质量流率,仅调节水温,同时计及回流管道热损;一次热网由具有供热设备的换热站提供热量,热用户通过二次热网获得热量。
5.根据权利要求1所述的计及碳捕集设备加装的低碳综合能源系统仿真及配置方法,其特征是,所述的目标函数包括:火电机组运行成本、燃气机组运行成本、碳相关成本、新能源弃置成本的目标函数,具体包括:
1)火电机组运行成本目标函数,其中:T为调度周期,Ce0,,t为启动成本,Cet,t为碳捕集设备加装成本,折算到每日运行,f(Pe,t)为发电成本函数,Pe,t是发电功率,ae、be和ce为火电机组的成本系数;
2)燃气机组运行成本目标函数,Cg=∑t∈Tpg,tVg,t,其中:pg,t为天然气价格;Vg,t为燃气机组天然气耗气量之和;
3)二氧化碳相关成本目标函数,包括系统化石燃料机组的碳交易成本、P2G设备的购碳成本和储碳设备的储碳成本Ccur=∑t∈Tpcur,tδtPre,t,其中:pc1,t为机组碳交易价格,μe为火电机组碳排放强度,μg为燃气机组碳排放强度,pc2,t为P2G设备购碳价格,Cp2g,t为P2G购碳量,ps,t碳存储价格,为碳储存量;
4)新能源弃置成本目标函数,Ccur=∑t∈Tpcur,tδtPre,t,其中:pcur,t为弃置惩罚系数,δt为弃置率,Pre,t为可用有功出力。
6.根据权利要求1所述的计及碳捕集设备加装的低碳综合能源系统仿真及配置方法,其特征是,所述的约束条件包括:电网约束和热网约束,具体包括:
i)电网功率平衡约束电网线路传输约束和电网电压相角约束其中:Pi,t为j节点的所连机组出力,AG为节点所连机组集合,Phj,t为支路hj的潮流,AF和AE为节点j为起点和终点的线路集合,和分别为潮流的上下限,Dj,t为节点j的电力负荷需求,θj,t为节点j的电压相角;
ii)热网约束:针对电力网络潮流和热力网络潮流惯性不同的问题,以电力网络潮流计算为基础,在热力网络中考虑延时特性,即通过的热能输运准动态过程温度约束 描述不同时段下进行热力网络的热水,此外,热网运行约束包含还热交换约束 汇流节点温度方程热源温度约束 c1、c2代表t时刻输出的热水质块中的权重系数,和为为质块的初始温度,Tam为环境温度,c3、c4为热水在管道中温度损失的系数,与管道长度和管道流速有关;和分别为节点v处的热源设备的热出力和用户热负荷;mH,v,t和mL,v,t为热源和热负荷处的热水质量;和为热源和热负荷的供水温度和回水温度;S(v)+和S(v)-分别为出水管道和入水管道与节点v相连的管道集合,为管道κ的出水温度,为节点v的流出温度;qb和qκ分别为管道b和κ的热水的质量流率;
iii)设备约束,由机组约束和存储设备约束构成,其中:
机组约束包括机组出力约束机组爬坡约束Pn,i,d≤Pn,i,t-Pn,i,t-1≤Pn,i,u,其中:Pn,i,t为t时刻第i个节点第n个设备的出力,和为该设备出力上下限;Pn,i,d和Pn,i,u为该设备爬坡上下限;
存储设备约束包括设备状态方程Si,t=Si,t-1+Si,in,t-Si,out,t总量约束Si,min≤Si,t≤Si,max,能量存储约束0≤Si,in,t≤Si,in,max(1-λi,t),能量释放约束0≤Si,out,t≤Si,out,maxλi,t,其中:Si,t为t时刻第i个设备能量储存总量;Si,in,t为该设备t时刻存储功率;Si,out,t为该设备t时刻释放功率;λi,t为01变量为设备状态。
7.根据权利要求1所述的计及碳捕集设备加装的低碳综合能源系统仿真及配置方法,其特征是,所述的仿真求解,具体包括:
步骤3.1:设计综合能源系统运行场景,其中:设备包括大型火力机组、电转气设备、燃气轮机、余热锅炉、燃气锅炉、光伏、电锅炉、地源热泵、电热储能,同时包含电热负荷用户;
步骤3.2:分析热力网络建模影响。
8.根据权利要求1所述的计及碳捕集设备加装的低碳综合能源系统仿真及配置方法,其特征是,所述的分析综合能源系统的最优配置方案,具体包括:
步骤4.1:配置方案分析;
步骤4.2:碳交易价格分析;
步骤4.3:碳价变化的影响分析;
步骤4.4:分析边际碳价和碳捕集加装的关系。
9.一种实现权利要求1~8中任一所述计及碳捕集设备加装的低碳综合能源系统仿真及配置方法的系统,其特征在于,包括:综合能源系统建模单元、综合能源系统调度单元以及综合能源系统配置单元,其中:综合能源系统建模单元根据所选的设备和网络信息,进行数学建模,得到设备和网络模型,综合能源系统调度单元根据综合能源系统的运行目标和约束条件,进行数学建模,得到调度模型,综合能源系统配置单元根据设备和网络模型和调度模型,基于matlab平台进行处理,通过Yalmip工具箱调用Gurobi求解器进行求解,得到综合能源系统配置方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211140985.9A CN115455709A (zh) | 2022-09-20 | 2022-09-20 | 计及碳捕集设备加装的低碳综合能源系统仿真及配置方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211140985.9A CN115455709A (zh) | 2022-09-20 | 2022-09-20 | 计及碳捕集设备加装的低碳综合能源系统仿真及配置方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115455709A true CN115455709A (zh) | 2022-12-09 |
Family
ID=84304602
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211140985.9A Pending CN115455709A (zh) | 2022-09-20 | 2022-09-20 | 计及碳捕集设备加装的低碳综合能源系统仿真及配置方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115455709A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116305799A (zh) * | 2023-02-07 | 2023-06-23 | 新疆敦华绿碳技术股份有限公司 | 一种碳捕集方法和系统 |
CN117114262A (zh) * | 2023-04-07 | 2023-11-24 | 天津大学 | 一种碳利用点供给源的选择与分配方法 |
-
2022
- 2022-09-20 CN CN202211140985.9A patent/CN115455709A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116305799A (zh) * | 2023-02-07 | 2023-06-23 | 新疆敦华绿碳技术股份有限公司 | 一种碳捕集方法和系统 |
CN116305799B (zh) * | 2023-02-07 | 2024-05-24 | 新疆敦华绿碳技术股份有限公司 | 一种碳捕集方法和系统 |
CN117114262A (zh) * | 2023-04-07 | 2023-11-24 | 天津大学 | 一种碳利用点供给源的选择与分配方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107482688B (zh) | 一种碳捕集虚拟电厂调度优化方法 | |
CN109494809B (zh) | 包含电转气装置的电-气互联系统及运行方法 | |
CN115455709A (zh) | 计及碳捕集设备加装的低碳综合能源系统仿真及配置方法 | |
CN111681130A (zh) | 一种考虑条件风险价值的综合能源系统优化调度方法 | |
CN114169727A (zh) | 一种考虑碳捕集和电转气协同的多能流综合能源低碳调度方法 | |
CN114330827B (zh) | 多能流虚拟电厂分布式鲁棒自调度优化方法及其应用 | |
CN116341881B (zh) | 计及热网灵活性的电-热系统鲁棒超前调度方法及系统 | |
Zhang et al. | Bi-level optimization dispatch of integrated-energy systems with P2G and carbon capture | |
Qin et al. | Multi-timescale hierarchical scheduling of an integrated energy system considering system inertia | |
CN114066204A (zh) | 一种综合能源系统的一体化优化规划与运行方法及装置 | |
CN114936720A (zh) | 一种源-荷-储协调的综合能源系统低碳经济调度方法 | |
CN115659651A (zh) | 一种考虑多种灵活性资源的综合能源协同优化调度方法 | |
CN113128868B (zh) | 一种区域综合能源系统调度优化方法及装置 | |
Qi et al. | A Distributionally Robust Optimization Scheduling Model for Regional Integrated Energy Systems Considering Hot Dry Rock Co-Generation. | |
CN116562062B (zh) | 考虑电、氢耦合的生产模拟方法、系统、装置及存储介质 | |
Rahbari et al. | Real-time modeling and optimization of molten salt storage with supercritical steam cycle for sustainable power generation and grid support | |
Tabar et al. | Impact of bi-facial PV panels’ presence as the novel option on the energy management and scheduling of the interconnected grids: Comprehensive outlook | |
CN110544033B (zh) | 热电厂灵活性改造后的电力系统风电消纳评估方法 | |
CN110472364B (zh) | 一种考虑可再生能源的离网式热电气联产系统的优化方法 | |
CN108062022A (zh) | 一种热电协同系统的优化控制方法 | |
CN115204705A (zh) | 考虑电转气存储的区域综合能源系统运行优化方法及应用 | |
CN113077173B (zh) | 一种基于综合能源系统节点能价策略的综合需求响应方法 | |
CN115495906A (zh) | 高比例可再生能源接入的配网多能互补协同优化方法 | |
CN114386256A (zh) | 一种考虑电热设备灵活性约束及热网特性的区域电热系统优化调度方法 | |
CN114066056A (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 |