CN111489083A - 考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法 - Google Patents
考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法 Download PDFInfo
- Publication number
- CN111489083A CN111489083A CN202010281250.2A CN202010281250A CN111489083A CN 111489083 A CN111489083 A CN 111489083A CN 202010281250 A CN202010281250 A CN 202010281250A CN 111489083 A CN111489083 A CN 111489083A
- Authority
- CN
- China
- Prior art keywords
- oxygen
- gas
- power plant
- heat
- carbon
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 title claims abstract description 389
- 229910052760 oxygen Inorganic materials 0.000 title claims abstract description 389
- 239000001301 oxygen Substances 0.000 title claims abstract description 389
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 236
- 238000002485 combustion reaction Methods 0.000 title claims abstract description 165
- 238000005516 engineering process Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 19
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 203
- 238000012423 maintenance Methods 0.000 claims abstract description 12
- 239000007789 gas Substances 0.000 claims description 259
- 238000003860 storage Methods 0.000 claims description 81
- 238000000926 separation method Methods 0.000 claims description 75
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 62
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 40
- 238000005338 heat storage Methods 0.000 claims description 40
- 238000000746 purification Methods 0.000 claims description 40
- 238000007906 compression Methods 0.000 claims description 39
- 230000006835 compression Effects 0.000 claims description 39
- 239000003345 natural gas Substances 0.000 claims description 31
- 238000006243 chemical reaction Methods 0.000 claims description 25
- 230000009194 climbing Effects 0.000 claims description 22
- 238000012946 outsourcing Methods 0.000 claims description 14
- 238000010521 absorption reaction Methods 0.000 claims description 12
- 230000008929 regeneration Effects 0.000 claims description 12
- 238000011069 regeneration method Methods 0.000 claims description 12
- 238000004519 manufacturing process Methods 0.000 claims description 11
- 239000000446 fuel Substances 0.000 claims description 10
- 238000007599 discharging Methods 0.000 claims description 9
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 8
- 230000005611 electricity Effects 0.000 claims description 8
- 238000006477 desulfuration reaction Methods 0.000 claims description 6
- 230000023556 desulfurization Effects 0.000 claims description 6
- 238000010438 heat treatment Methods 0.000 claims description 6
- 238000006213 oxygenation reaction Methods 0.000 claims description 6
- 238000002407 reforming Methods 0.000 claims description 6
- 230000009919 sequestration Effects 0.000 claims description 6
- 239000001569 carbon dioxide Substances 0.000 claims description 4
- 238000007789 sealing Methods 0.000 claims description 4
- 239000002699 waste material Substances 0.000 claims description 4
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 3
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 claims description 3
- 230000033228 biological regulation Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 239000003546 flue gas Substances 0.000 claims description 3
- 239000002994 raw material Substances 0.000 claims description 3
- 239000002918 waste heat Substances 0.000 claims description 3
- 238000009413 insulation Methods 0.000 claims description 2
- 230000002035 prolonged effect Effects 0.000 claims description 2
- 230000007246 mechanism Effects 0.000 abstract description 12
- 230000000694 effects Effects 0.000 abstract description 7
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000011049 filling Methods 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000029305 taxis Effects 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Images
Classifications
-
- 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/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
- G06Q10/06375—Prediction of business process outcome or impact based on a proposed change
-
- 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
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems 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)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Separation By Low-Temperature Treatments (AREA)
Abstract
本发明提出一种考虑富氧燃烧技术的电‑气‑热综合能源系统低碳经济调度方法,其特点是,包括:富氧燃烧电厂电能时移特性分析、富氧燃烧电厂净输出功率模型构建、富氧燃烧电厂的电碳特性分析、电‑气‑热综合能源系统低碳经济调度模型构建等步骤,通过挖掘富氧燃烧电厂运行机理,研究其内部电能时移特性,构建净输出功模型并分析其电碳特性;其次,考虑电‑气‑热综合能源系统对富氧燃烧电厂的协调作用,基于其运行机理设定约束条件,建立以运行成本、维护成本、折旧成本最低为目标函数的电‑气‑热综合能源系统低碳经济调度模型。具有科学合理,适用性强,效果佳等优点。
Description
技术领域
本发明涉及一种考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法。
背景技术
随着经济的不断发展,化石燃料被大量利用,由此带来的全球变暖问题日趋加剧。我国电力行业是主要的能源消耗者,具有较大的碳排放,但通过能源转型能实现碳减排,构建综合能源系统是实现能源转型的重要举措。综合能源系统是一种内含电、气、热等多种能源,通过能源转换设施、储存设施的协调优化所形成的产供销一体化系统,可实现能源的高效利用,而综合能源系统低碳运行的研究主要集中在两方面,其一是政策类,如建立碳交易机制等;其二是技术类,如碳捕集技术等。通过二者协调配合能够实现整体系统的低碳经济运行。
当前,将碳交易引入综合能源系统是实现低碳运行的主要研究方法,碳交易是一种通过建立合法的碳排放权,并可利用市场对碳排放权进行交易的市场机制。文献[1]将碳交易引入电-气综合能源系统,证明碳交易有利于综合系统的低碳性和经济性,推动低碳机组的发展。文献[2]在电-气-热区域综合能源系统中引入碳交易,证明碳交易的设立有利于清洁机组上网,推动系统低碳经济运行。
上述研究中,火电作为二氧化碳的主要产生源仍为综合能源系统供应基荷,因此仅引入碳交易机制并不能实现碳排量的骤降,需要形成技术手段和政策机制协调配合,如文献[3-5]。
文献[3-5]通过将碳捕集技术与碳税、碳交易机制等市场机制相结合,实现综合能源系统的低碳经济运行。碳捕集技术指对燃煤电厂进行低碳化改造的技术,主要分为燃烧前捕集、燃烧后捕集、富氧燃烧技术,其中燃烧前捕集技术只能应用于新建电厂,而富氧燃烧技术和燃烧后捕集技术能对传统燃煤电厂进行低碳化改造,文献[3-5]所应用的碳捕集技术为燃烧后捕集技术,存在改造投资成本高、最高碳捕集水平较低等缺陷。而近年来碳捕集技术中的富氧燃烧技术对上述问题有所改善,但鲜有在综合能源系统中引入富氧燃烧技术并配合碳交易机制的研究。
富氧燃烧技术是碳捕集技术的一种,其对传统燃煤电厂进行低碳化改造,工作原理是将氧气与循环CO2流的混合气代替空气输入燃煤电厂的燃烧室,产生高浓度CO2流,更易于CO2压缩纯化装置进行捕获。富氧燃烧电厂具有相对成本低、易规模化、可改造存量机组等诸多优势,被认为是最可能大规模推广和商业化的碳捕集技术之一。
发明内容
本发明所要解决的技术问题是,提出一种考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法,其科学合理,适用性强,效果佳,成本低,能够计及综合能源系统的协调配合,提高经济性和降低碳排放量。
解决其技术问题采用的技术方案是:一种考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法,其特征是,它通过考虑富氧燃烧电厂的运行特点及电-气-热综合能源系统的协调配合,构建了电-气-热综合能源系统协调调度模型,其具体包括以下步骤:
步骤一,富氧燃烧电厂电能时移特性分析:
富氧燃烧电厂的电能损耗主要集中于空分制氧装置和碳捕集压缩纯化装置;为应对快速变负荷要求,空分制氧装置配有储存液态氧的隔热储氧罐;富氧燃烧电厂的电能主要流向电负荷、空分制氧装置和碳捕集压缩纯化装置,正常运行情况下,电负荷用电能最多,空分制氧装置其次,碳捕集压缩纯化装置最少;当电负荷过多时通过降低空分制氧装置和碳捕集压缩纯化装置的功率为电负荷提供支撑,而富氧燃烧过程所需氧气量和燃料量成正比,此时空分制氧装置不能产生足够的氧气,但缺少的氧气能从储氧罐中获取;因此,富氧燃烧电厂能够利用储氧罐在多时段进行电能调节,在电负荷低谷时让空分制氧装置满功率运行,在电负荷高峰时减少空分制氧装置的出力,其中的氧气差值利用储氧罐进行调整,实现电能的时移;
步骤二,富氧燃烧电厂净输出功率模型构建:
富氧燃烧电厂产生的电能主要流向空分制氧装置、碳捕集压缩纯化装置与电负荷,用(1)式表示:
PZ=PN+PCPU+PASU (1)
式中,PZ为火电厂的输出总功率;PN为火电厂的净输出功率;PCPU为碳捕集压缩纯化装置所需功率;PASU为空分制氧装置所需功率;
设常规火电厂的碳排放强度为EG,其燃烧化石能源所产生的CO2总量EZ为:
EZ=EGPZ (2)
若富氧燃烧电厂的碳捕集水平为α,则所能捕获的CO2量EB为:
EB=αEZ (3)
此时,空分制氧装置所需功率PASU和碳捕集压缩纯化装置所需功率PCPU为(4)式:
式中,OC为空分制氧装置产生氧气量;OCQ为储氧罐供氧量;OZ为进行富氧燃烧所需氧气量;β为空分制氧装置单位时间制造单位氧气所需电能;δ为碳捕集压缩纯化装置单位时间捕捉单位CO2所需电能;χ为富氧燃烧电厂进行富氧燃烧产生单位功率所需氧气量;
从而推出净输出功率为(5)式:
PN=(1-αδEG-βχ)PZ+βOCQ (5)
从中得知富氧燃烧电厂净输出功率和碳捕集水平、空分制氧装置的制氧量、储气罐供氧量有关,改变储氧罐的充放氧能够改变某一时段的净输出功率,因此富氧燃烧能够利用储氧罐做到电能时移;
富氧燃烧电厂存在两种运行模式,一是富氧运行,二是空气运行,这两种模式能够相互切换,另外在空气燃烧的过程中,碳捕集压缩纯化装置不进行捕碳,但这并不影响空分制氧装置的制氧,其能将产生的氧气装入储氧罐中;
步骤三,富氧燃烧电厂的电碳特性分析:
考虑空分制氧装置是否投运对富氧燃烧电厂运行的影响,分别对空分制氧装置投运时和空分制氧装置停运时的电碳特性进行分析,电碳特性指电厂净输出功率和二氧化碳净排放量的关系;
1)空分制氧装置投运
在空分制氧装置运行时,由于存在储氧罐,空分制氧装置能够在限定内的功率区间运行,即在火电厂最小运行功率时,能够使空分制氧装置以最大功率运行,在火电厂最大运行功率时,使空分制氧装置以最小功率运行,因此火电厂具有更深的调峰深度;而碳捕集压缩纯化装置能够根据当前火电厂输出总功率调节碳捕集电能;净碳排放量为火电厂实际排放量减去捕碳量,则空分制氧装置运行时电碳特性关系为(6)式:
式中,EJP为净排碳量;αmax为最大碳捕集水平;
2)空分制氧装置停运
在空分制氧装置停运时,由于存在储氧罐,富氧燃烧电厂依旧能工作在碳捕集状态,则空分制氧装置停运时电碳特性关系为(7)式:
当储氧罐的储氧量不足时,富氧燃烧电厂需退出富氧运行状态,转为空气运行状态;
步骤四,电-气-热综合能源系统低碳经济调度模型构建:
电-气-热综合能源系统,其中,电负荷、电转气设备的用电由风电场、火电厂、燃气轮机进行电能的供给,热负荷的用热由燃气轮机、燃气锅炉、储热罐进行热能供给,气负荷、燃气轮机、燃气锅炉的用气由电转气设备、气源、储气罐进行气能供给;
利用储热罐协调燃气轮机出力支撑空分制氧装置和碳捕集压缩纯化装置的运行;燃气轮机是以天然气为原材料,通过燃烧天然气使空气膨胀推动汽轮机转动从而获取电能,其高温烟气通过余热锅炉能够获取热能,因此会出现以热定电的情况;通过加入储热罐能在一定程度打破以热定电,从而使燃气轮机能在电负荷峰值时提供更多的电能,减小火电厂净出力,增加对空分制氧装置和碳捕集压缩纯化装置的电能输送,提高碳捕集水平;
在电-气-热综合能源系统中,电转气设备能够将弃风接纳转换成天然气,同时能接收一部分富氧燃烧电厂捕捉到的CO2,从而节约封存成本和购买高纯度CO2成本,使经济性得到改善;
以电-气-热综合能源系统运行成本、碳捕集设备折旧成本以及维护成本最小为目标,考虑各机组、元件的约束条件构建低碳经济调度模型;其中电-气-热综合能源系统运行成本包括碳交易成本、火电燃料成本、外购天然气成本、运输封存CO2成本、外购CO2成本,而碳捕集设备包括燃烧后捕集电厂中的吸收塔、再生塔、压缩机以及富氧燃烧电厂中的碳捕集压缩纯化装置、空分制氧装置,目标函数如下:
min F=min(CT+CH+CQ+CW+CYF+CMCO2+CTB)
式中,F为电-气-热综合能源系统总成本;CT为电-气-热综合能源系统碳交易成本;CH为电-气-热综合能源系统火电燃料成本;CQ为外购天然气成本;CW为火电厂维护成本;CYF为运输封存CO2成本;CMCO2为外购CO2成本,指购买高纯度CO2供给电转气设备的成本;CTB为碳捕集设备折旧成本;σ为单位碳排放权交易价格;Ec为电-气-热综合能源系统的实际碳排放量;N1为火电厂的个数;N2为燃气机组的个数;Dk1为第k1个火电厂的碳排量配额;Dk2为第k2个燃气机组的碳排放配额;Pk1,t为第k1个火电厂t时段的输出电功率;ak1、bk1、ck1为火电厂的成本系数;μgas为购买天然气价格;GRZ为天然气热值;Zer为单位电能所能转化的热能;Rrg,k2,t为第k2个燃气机组为燃气锅炉时t时段输出的热功率;Pk2,t为第k2个燃气机组为燃气轮机时t时段输出的电功率;ηrg为燃气锅炉的热效率;ηrq为燃气轮机的电效率;QP2G,t为t时段电转气设备产生的天然气量;CCG为火电厂成本;ρCG为火电厂维护系数;CTL为脱硫装置成本;ρTL为脱硫装置维护系数;CTX为脱硝装置成本;ρTX为脱硝装置维护系数;CGZ为富氧燃烧改造锅炉成本;CASU为空分制氧装置成本;ρASU为空分制氧装置维护系数;CCPU为碳捕集压缩纯化装置成本;ρCPU为碳捕集压缩纯化装置维护系数;CRHGZ为燃烧后捕集电厂锅炉改造成本;CTBJ为燃烧后捕集电厂中吸收塔、再生塔、压缩机的成本;ρTB为燃烧后捕集电厂中吸收塔、再生塔、压缩机的维护系数;μYF为运输封存单位CO2的成本;BCO2为捕捉到的CO2总量;BP2G为捕捉到的CO2供电转气设备用掉的量;μCO2为外购单位CO2所需价格;θCO2为输入电转气设备单位电功率所需CO2量;PP2G,t为t时段输入电转气设备的电功率;NT为富氧燃烧电厂中空分制氧装置、碳捕集压缩纯化装置和燃烧后捕集电厂中吸收塔、再生塔、压缩机的折旧年限;ω为折旧率;PCY为单位体积储氧罐价格;VCY为储氧罐体积;NC为储氧罐折旧年限;
约束条件包括功率平衡约束、富氧燃烧电厂约束、电转气设备约束、风电场约束、燃气轮机约束、燃气锅炉约束、储氧\热\气罐约束;
1)功率平衡约束
对于电负荷、热负荷、气负荷分别满足功率平衡约束,如下:
式中,Pel,t为t时段的电负荷;Pk1J,t为t时段k1火电厂净输出功率;Pk2,t为t时段k2燃气机组为燃气轮机时所供电功率;Pwind,t为风电在t时段的输出电功率;PP2G,t为t时段输入电转气设备的电功率;Rrl,t为t时段的热负荷;Rrg,k2,t为t时段k2燃气机组为燃气锅炉时所供热功率;Rrq,k2,t为t时段k2燃气机组为燃气轮机时所供热功率;σrx为储热罐充放热效率;RXr,t为t时段储热罐充热量;RFr,t为t时段储热罐放热量;Qql,t为t时段的气负荷;QQY,t为t时段的气源出力;QXQ,t为t时段的储气罐充气量;QFQ,t为t时段的储气罐放气量;QP2G,t为t时段的电转气设备产气量;Qrg,k2,t为t时段k2燃气机组为燃气锅炉时所需天然气量;Qrq,k2,t为t时段k2燃气机组为燃气轮机时所需天然气量;
2)富氧燃烧电厂约束
富氧燃烧电厂需要满足火电平衡约束、空分制氧装置制氧约束、空分制氧装置出力约束、碳捕集压缩纯化装置工作约束、捕碳量约束、碳捕集水平约束、净排碳量约束、氧气平衡约束、富氧燃烧所需氧气量约束、火电出力限值约束、火电出力爬坡约束为(10)式:
式中,Pkl,t为t时段k1火电厂的总电功率;PASUk1,t为k1火电厂在t时段的空分制氧装置所需功率;PCPUk1,t为k1火电厂在t时段碳捕集压缩纯化装置的所需功率;β为空分制氧装置单位时间制造单位氧气所需电能;OCk1,t为k1火电厂的空分制氧装置在t时段产生氧气量;λASU为空分制氧装置最小运行系数;PASUk1max为火电厂k1的空分制氧装置最大功率;δ为碳捕集压缩纯化装置单位时间捕捉单位CO2所需电能;EBk1,t为k1火电厂在t时段的捕碳量;ζ为0-1变量,ζ=1表示富氧燃烧电厂处于富氧燃烧状态,ζ=0表示富氧燃烧电厂不处于富氧燃烧状态;α为碳捕集水平;EGk1为k1火电厂的碳排放强度;αmax为最大碳捕集水平;EJPk1,t为k1火电厂在t时段的净排碳量;OZk1,t为k1火电厂在t时段进行富氧燃烧所需氧气量;OFQk1,t为火电厂k1的储氧罐t时段放氧量;OXQk1,t为火电厂k1的储氧罐t时段充氧量;OQYLk1,t为t时段k1火电厂允许弃氧量;χ为富氧燃烧电厂进行富氧燃烧产生单位功率所需氧气量;Pk1min为火电厂k1最小出力值;Pk1max为火电厂k1最大出力值;Pk1Xmin为火电厂k1最大下爬坡率;Pk1Smax为火电厂k1最大上爬坡率;
3)电转气设备约束
电转气设备主要满足其最大功率约束为:
0≤PP2G,t≤PP2Gmax (11)
式中,PP2Gmax为电转气设备运行的最大功率;
4)风电场约束
风电场主要满足其最大功率约束为:
0≤Pwind,t≤PW,t (12)
式中,Pw,t为风电预测功率;
5)燃气轮机约束
燃气轮机主要满足出力上下限约束、爬坡率约束和热电比约束为:
Pk2min≤Pk2,t≤Pk2max (13)
Pk2Xmin≤Pk2,t-Pk2,t-1≤Pk2Smax (14)
式中,Pk2min为第k2个燃气机组为燃气轮机时的最小电功率;Pk2max为第k2个燃气机组为燃气轮机时的最大电功率;Pk2Xmin为第k2个燃气机组为燃气轮机时的最大下爬坡率;Pk2Smax为第k2个燃气机组为燃气轮机时的最大上爬坡率;ΨREC为燃气轮机的热电比;
6)燃气锅炉约束
燃气锅炉主要满足出力上下限约束和爬坡率约束为:
Rrg,k2min≤Rrg,k2,t≤Rrg,k2max (16)
Rrg,k2Xmin≤Rrg,k2,t-Rrg,k2,t-1≤Rrg,k2Smax (17)
式中,Rrg,k2min为第k2个燃气机组为燃气锅炉时最小热功率;Rrg,k2max为第k2个燃气机组为燃气锅炉时最大热功率;Rrg,k2Xmin为第k2个燃气机组为燃气锅炉时的最大下爬坡率;Rrg,k2Smax为第k2个燃气机组为燃气锅炉时的最大上爬坡率;
7)储氧/气/热罐约束
约束包括储氧/气/热容量上下限约束,储氧/气/热周期始末相等约束,充放氧/气/热限值约束,储氧/气/热平衡约束为:
0≤Echu,t≤Echumax (18)
Echu,1=Echu,24 (19)
0≤EX,t≤ωEXmax (20)
0≤EF,t≤(1-ω)EFmax (21)
Echu,t=Echu,t-1+EX,t-EF,t (22)
式中,Echu,t为储氧/气/热罐在t时段的储氧/气/热量;Echumax为储氧/气/热罐容量上限;Echu,1为储氧/气/热罐在周期内初始储氧/气/热量;Echu,24为储氧/气/热罐在周期结束时储氧/气/热量;EX,t为t时段的充氧/气/热量;EF,t为t时段的放氧/气/热量;EXmax和EFmax分别为单位时段充氧/气/热和放氧/气/热量上限;ω为二进制变量,当ω=1时进行充氧/气/热,当ω=0时进行放氧/气/热。
本发明提出一种考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法,其特点是,将富氧燃烧技术引入电-气-热综合能源系统,考虑富氧燃烧技术的投资成本和低碳效益,同时计及综合能源系统的协调配合,提高经济性的同时降低碳排量。首先,通过挖掘富氧燃烧电厂运行机理,研究其内部电能时移特性,构建净输出功模型并分析其电碳特性;其次,考虑电-气-热综合能源系统对富氧燃烧电厂的协调作用,基于其运行机理设定约束条件,建立以运行成本、维护成本、折旧成本最低为目标函数的电-气-热综合能源系统低碳经济调度模型。具有科学合理,适用性强,效果佳等优点。
附图说明
图1本发明的考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法方框图;
图2是燃煤火电厂分别为常规电厂、燃烧后捕集电厂、富氧燃烧电厂下火电厂出力图;
图3是燃煤火电厂分别燃烧后捕集电厂、富氧燃烧电厂下系统碳捕集总量图;
图4是燃煤火电厂分别为常规电厂、燃烧后捕集电厂、富氧燃烧电厂下系统净碳排量图。
具体实施方式
下面利用附图和实施例对本发明一种考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法作进一步说明。
参照图1,本发明的一种考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法,其包括富氧燃烧电厂的运行特点及与电-气-热综合能源系统的协调配合,构建了电-气-热综合能源系统协调调度模型,其具体包括以下步骤:
步骤一,富氧燃烧电厂电能时移特性分析:
与常规火电厂相比,富氧燃烧电厂的电能损耗主要集中于空分制氧装置(AirSeparation Unit,ASU)和碳捕集压缩纯化装置(Compression Purification Unit,CPU);为应对快速变负荷要求,ASU一般配有隔热储氧罐,能接收ASU制造的液态氧;富氧燃烧电厂的电能主要流向三个部分,分别为电负荷、ASU、碳捕集CPU,正常运行情况下,电负荷用电能最多,ASU其次,碳捕集CPU最少;当电负荷过多时通过降低ASU和碳捕集CPU的功率为电负荷提供支撑,而富氧燃烧过程所需氧气量和燃料量成正比,此时ASU不能产生足够的氧气,但缺少的氧气能从储氧罐中获取;因此,富氧燃烧电厂能够利用储氧罐在多时段进行电能调节,在电负荷低谷时让ASU满发,在电负荷高峰时减少ASU的出力,其中的的氧气差值利用储氧罐进行调整,实现电能的时移;
步骤二,富氧燃烧电厂净输出功率模型构建:
富氧燃烧电厂产生的电能主要流向ASU、碳捕集CPU与电负荷,用(1)式表示:
PZ=PN+PCPU+PASU (1)
式中,PZ为火电厂的输出总功率;PN为火电厂的净输出功率;PCPU为碳捕集CPU所需功率;PASU为ASU所需功率;
设常规火电厂的碳排放强度为EG,其燃烧化石能源所产生的CO2总量EZ为:
EZ=EGPZ (2)
若富氧燃烧电厂的碳捕集水平为α,则所能捕获的CO2量EB为:
EB=αEZ (3)
此时,ASU所需功率PASU和碳捕集CPU所需功率PCPU为(4)式:
式中,OC为ASU产生氧气量;OCQ为储氧罐供氧量;OZ为进行富氧燃烧所需氧气量;β为ASU单位时间制造单位氧气所需电能;δ为碳捕集CPU单位时间捕捉单位CO2所需电能;χ为富氧燃烧电厂进行富氧燃烧产生单位功率所需氧气量;
从而推出净输出功率为(5)式:
PN=(1-αδEG-βχ)PZ+βOCQ (5)
从中得知富氧燃烧火电厂净输出功率和碳捕集水平、ASU制氧量、储气罐供氧量有关,改变储氧罐的充放氧能够改变某一时段的净输出功率,因此富氧燃烧能够利用储氧罐做到电能时移;
富氧燃烧电厂存在两种运行模式,一是富氧运行,二是空气运行,这两种模式能够相互切换,另外在空气燃烧的过程中,碳捕集CPU不进行捕碳,但这并不影响ASU制氧,其能将产生的氧气装入储氧罐中;
步骤三,富氧燃烧电厂的电碳特性分析:
考虑ASU是否投运对富氧燃烧电厂运行的影响,分别对ASU投运时和ASU停运时的电碳特性进行分析,电碳特性指电厂净输出功率和二氧化碳净排放量的关系;
1)ASU投运
在ASU运行时,由于存在储氧罐,ASU能够在限定内的功率区间运行,即在火电厂最小运行功率时能使ASU以最大功率运行,在火电厂最大运行功率时使ASU以最小功率运行,因此火电厂具有更深的调峰深度;而碳捕集CPU能够根据当前火电厂输出总功率调节碳捕集电能;净碳排放量为火电厂实际排放量减去捕碳量,则ASU投运时电碳特性关系式为(6)式:
式中,EJP为净排碳量;αmax为最大碳捕集水平;
2)ASU停运
在ASU停运时,由于存在储氧罐,富氧燃烧电厂依旧能工作在碳捕集状态,则ASU停运时电碳特性关系式为(7)式:
当储氧罐的储氧量不足时,富氧燃烧电厂需退出富氧运行状态,转为空气运行状态;
步骤四,电-气-热综合能源系统低碳经济调度模型构建:
电-气-热综合能源系统,其中,电负荷、电转气(power to gas,P2G)设备的用电由风电场、火电厂、燃气轮机进行电能的供给,热负荷的用热由燃气轮机、燃气锅炉、储热罐进行热能供给,气负荷、燃气轮机、燃气锅炉的用气由P2G设备、气源、储气罐进行气能供给;
利用储热罐协调燃气轮机出力支撑ASU和碳捕集CPU的运行;燃气轮机是以天然气为原材料,通过燃烧天然气使空气膨胀推动汽轮机转动从而获取电能,其高温烟气通过余热锅炉能够获取热能,因此会出现以热定电的情况;通过加入储热罐能在一定程度打破以热定电,从而使燃气轮机能在电负荷峰值时提供更多的电能,减小火电厂净出力,增加对ASU和碳捕集CPU的电能输送,提高碳捕集水平;
在电-气-热综合能源系统中,P2G设备能够将弃风接纳转换成天然气,同时能接收一部分富氧燃烧电厂捕捉到的CO2,从而节约一定的封存成本和购买高纯度CO2成本,使经济性得到改善;
以电-气-热综合能源系统运行成本、碳捕集设备折旧成本以及维护成本最小为目标,考虑各机组、元件的约束条件构建低碳经济调度模型;其中,电-气-热综合能源系统运行成本包括碳交易成本、火电燃料成本、外购天然气成本、运输封存CO2成本、外购CO2成本,而碳捕集设备包括燃烧后捕集电厂中的吸收塔、再生塔、压缩机以及富氧燃烧电厂中的碳捕集CPU、ASU,目标函数如下:
min F=min(CT+CH+CQ+CW+CYF+CMCO2+CTB)
式中,F为电-气-热综合能源系统总成本;CT为电-气-热综合能源系统碳交易成本;CH为电-气-热综合能源系统火电燃料成本;CQ为外购天然气成本;CW为火电厂维护成本;CYF为运输封存CO2成本;CMCO2为外购CO2成本,指购买高纯度CO2供给P2G设备的成本;CTB为碳捕集设备折旧成本;σ为单位碳排放权交易价格;Ec为电-气-热综合能源系统的实际碳排放量;N1为火电厂的个数;N2为燃气机组的个数;Dk1为第k1个火电厂的碳排量配额;Dk2为第k2个燃气机组的碳排放配额;Pk1,t为第k1个火电厂t时段的输出电功率;ak1、bk1、ck1为火电厂的成本系数;μgas为购买天然气价格;GRZ为天然气热值;Zer为单位电能所能转化的热能;Rrg,k2,t为第k2个燃气机组为燃气锅炉时t时段输出的热功率;Pk2,t为第k2个燃气机组为燃气轮机时t时段输出的电功率;ηrg为燃气锅炉的热效率;ηrq为燃气轮机的电效率;QP2G,t为t时段P2G设备产生的天然气量;CCG为火电厂成本;ρCG为火电厂维护系数;CTL为脱硫装置成本;ρTL为脱硫装置维护系数;CTX为脱硝装置成本;ρTX为脱硝装置维护系数;CGZ为富氧燃烧改造锅炉成本;CASU为ASU成本;ρASU为ASU维护系数;CCPU为碳捕集CPU成本;ρCPU为碳捕集CPU维护系数;CRHGZ为燃烧后捕集电厂锅炉改造成本;CTBJ为燃烧后捕集电厂中吸收塔、再生塔、压缩机的成本;ρTB为燃烧后捕集电厂中吸收塔、再生塔、压缩机的维护系数;μYF为运输封存单位CO2的成本;BCO2为捕捉到的CO2总量;BP2G为捕捉到的CO2供P2G设备用掉的量;μCO2为外购单位CO2所需价格;θCO2为输入P2G设备单位电功率所需CO2量;PP2G,t为t时段输入P2G设备的电功率;NT为富氧燃烧电厂中ASU、碳捕集CPU和燃烧后捕集电厂中吸收塔、再生塔、压缩机的折旧年限;ω为折旧率;PCY为单位体积储氧罐价格;VCY为储氧罐体积;NC为储氧罐折旧年限;
约束条件包括功率平衡约束、富氧燃烧电厂约束、P2G设备约束、风电场约束、燃气轮机约束、燃气锅炉约束、储氧\热\气罐约束;
1)功率平衡约束
对于电负荷、热负荷、气负荷分别满足功率平衡约束,如下:
式中,Pel,t为t时段的电负荷;Pk1J,t为t时段k1火电厂净输出功率;Pk2,t为t时段k2燃气机组为燃气轮机时所供电功率;Pwind,t为风电在t时段的输出电功率;PP2G,t为t时段输入P2G设备的电功率;Rrl,t为t时段的热负荷;Rrg,k2,t为t时段k2燃气机组为燃气锅炉时所供热功率;Rrq,k2,t为t时段k2燃气机组为燃气轮机时所供热功率;σrx为储热罐充放热效率;RXr,t为t时段储热罐充热量;RFr,t为t时段储热罐放热量;Qql,t为t时段的气负荷;QQY,t为t时段的气源出力;QXQ,t为t时段的储气罐充气量;QFQ,t为t时段的储气罐放气量;QP2G,t为t时段的P2G设备产气量;Qrg,k2,t为t时段k2燃气机组为燃气锅炉时所需天然气量;Qrq,k2,t为t时段k2燃气机组为燃气轮机时所需天然气量;
2)富氧燃烧电厂约束
富氧燃烧电厂需要满足火电平衡约束、ASU制氧约束、ASU出力约束、碳捕集CPU工作约束、捕碳量约束、碳捕集水平约束、净排碳量约束、氧气平衡约束、富氧燃烧所需氧气量约束、火电出力限值约束、火电出力爬坡约束为(10)式:
式中,Pkl,t为t时段k1火电厂的总电功率;PASUk1,t为k1火电厂在t时段的ASU所需功率;PCPUk1,t为k1火电厂在t时段碳捕集CPU的所需功率;β为ASU单位时间制造单位氧气所需电能;OCk1,t为k1火电厂的ASU在t时段产生氧气量;λASU为ASU最小运行系数;PASUk1max为火电厂k1的ASU最大功率;δ为碳捕集CPU单位时间捕捉单位CO2所需电能;EBk1,t为k1火电厂在t时段的捕碳量;ζ为0-1变量,ζ=1表示富氧燃烧电厂处于富氧燃烧状态,ζ=0表示富氧燃烧电厂不处于富氧燃烧状态;α为碳捕集水平;EGk1为k1火电厂的碳排放强度;αmax为最大碳捕集水平;EJPk1,t为k1火电厂在t时段的净排碳量;OZk1,t为k1火电厂在t时段进行富氧燃烧所需氧气量;OFQk1,t为火电厂k1的储氧罐t时段放氧量;OXQk1,t为火电厂k1的储氧罐t时段充氧量;OQYLk1,t为t时段k1火电厂允许弃氧量;χ为富氧燃烧电厂进行富氧燃烧产生单位功率所需氧气量;Pk1min为火电厂k1最小出力值;Pk1max为火电厂k1最大出力值;Pk1Xmin为火电厂k1最大下爬坡率;Pk1Smax为火电厂k1最大上爬坡率;
3)P2G设备约束
P2G设备主要满足其最大功率约束为:
0≤PP2G,t≤PP2Gmax (11)
式中,PP2Gmax为P2G设备运行的最大功率;
4)风电场约束
风电场主要满足其最大功率约束为:
0≤Pwind,t≤PW,t (12)
式中,Pw,t为风电预测功率;
5)燃气轮机约束
燃气轮机主要满足出力上下限约束、爬坡率约束和热电比约束为:
Pk2min≤Pk2,t≤Pk2max (13)
Pk2Xmin≤Pk2,t-Pk2,t-1≤Pk2Smax (14)
式中,Pk2min为第k2个燃气机组为燃气轮机时的最小电功率;Pk2max为第k2个燃气机组为燃气轮机时的最大电功率;Pk2Xmin为第k2个燃气机组为燃气轮机时的最大下爬坡率;Pk2Smax为第k2个燃气机组为燃气轮机时的最大上爬坡率;ΨREC为燃气轮机的热电比;
6)燃气锅炉约束
燃气锅炉主要满足出力上下限约束和爬坡率约束为:
Rrg,k2min≤Rrg,k2,t≤Rrg,k2max (16)
Rrg,k2Xmin≤Rrg,k2,t-Rrg,k2,t-1≤Rrg,k2Smax (17)
式中,Rrg,k2min为第k2个燃气机组为燃气锅炉时最小热功率;Rrg,k2max为第k2个燃气机组为燃气锅炉时最大热功率;Rrg,k2Xmin为第k2个燃气机组为燃气锅炉时的最大下爬坡率;Rrg,k2Smax为第k2个燃气机组为燃气锅炉时的最大上爬坡率;
7)储氧/气/热罐约束
约束包括储氧/气/热容量上下限约束,储氧/气/热周期始末相等约束,充放氧/气/热限值约束,储氧/气/热平衡约束为:
0≤Echu,t≤Echumax (18)
Echu,1=Echu,24 (19)
0≤EX,t≤ωEXmax (20)
0≤EF,t≤(1-ω)EFmax (21)
Echu,t=Echu,t-1+EX,t-EF,t (22)
式中,Echu,t为储氧/气/热罐在t时段的储氧/气/热量;Echumax为储氧/气/热罐容量上限;Echu,1为储氧/气/热罐在周期内初始储氧/气/热量;Echu,24为储氧/气/热罐在周期结束时储氧/气/热量;EX,t为t时段的充氧/气/热量;EF,t为t时段的放氧/气/热量;EXmax和EFmax分别为单位时段充氧/气/热和放氧/气/热量上限;ω为二进制变量,当ω=1时进行充氧/气/热,当ω=0时进行放氧/气/热。
本实施例是以某电网实际负荷、风电出力为依据,在IEEE-30节点电网、7节点气网和6节点热网系统进行算例分析,并设定三种不同的场景进行对比,以此验证本方法所构建模型的有效性。
实施例计算条件说明如下:
1)IEEE-30节点电网、7节点气网和6节点热网系统中包括4个火电厂,容量分别为200MW、50MW、35MW、30MW,一个80MW燃气轮机,一个200MW风电场,一个150MW的P2G设备,以及一个80MW燃气锅炉,1个100MWh的储热罐,1个10000m3的储气罐,2个最大速率为50000m3/h的气源。
2)所设定的三种场景分别是:
场景一:燃煤火电厂为常规火电厂的电-气-热综合能源系统低碳调度模型;
场景二:将燃煤火电厂改造为燃烧后捕集电厂的电-气-热综合能源系统低碳调度模型;
场景三:将燃煤火电厂改造为富氧燃烧电厂的电-气-热综合能源系统低碳调度模型。
在上述场景下,应用本发明方法的优化调度结果如下:
1、富氧燃烧技术在减少系统碳排量和电能时移方面的作用与效果。
图2是燃煤火电厂分别为常规电厂、燃烧后捕集电厂、富氧燃烧电厂下火电厂出力图。图3是燃煤火电厂分别燃烧后捕集电厂、富氧燃烧电厂下系统碳捕集总量图。图4是燃煤火电厂分别为常规电厂、燃烧后捕集电厂、富氧燃烧电厂下系统净碳排量图。
从图2得知,富氧燃烧火电厂相较于燃烧后捕集电厂在“削峰填谷”方面更具优势,在1-9点、19-24点非火电高发时期大多数时间场景三火电出力高于场景二,在10-18点火电高发时期,场景三火电出力小于等于场景二,原因在于ASU和储氧罐的协调配合,实现了电能时移。
从图2和图3得知,场景二在电负荷峰值时(10-18点)不能达到正常的碳捕集水平(和火电厂输出总功率同比例变化),而场景三碳捕集量随其电厂输出功率同比例变化,原因在于富氧燃烧电厂的电能时移特性;其在电负荷低谷时利用ASU产生更多的氧气储存在氧气罐中,消耗更多的火电功率;在电负荷峰值时利用储存的氧气进行碳捕集,使ASU无需耗费过多的电能,同时不降低碳捕集水平。
从图4得知,相较于场景二,场景三的系统净碳排放量在24h内都较小,原因在于富氧燃烧电厂具有更高的碳捕集水平,同时通过储氧罐进行电能时移,使峰时有足够的电能捕碳,保持正常碳捕集水平,最终实现碳排量的降低。
2、验证本方法所构建综合能源系统对低碳经济调度的有效性
表1为三种场景下综合能源系统的火电燃料成本、外购天然气成本、运输封存CO2成本、外购CO2成本、碳捕集设备折旧成本、维护成本、总成本与排碳量的具体数值。
表1系统成本表
由表1能看出,场景二相对于场景一总成本上升了10424元,即0.437%,而碳排量减少了3659.7t,即67.43%;场景三相对于场景一总成本下降了51726元,即2.177%,且碳排量下降了4263.6t,即78.55%;场景三相对于场景二在总成本上下降了2.614%,在碳排量上下降了11.12%,能够证明富氧燃烧电厂在低碳经济调度方面的有效性。
参考文献:
[1]卫志农,张思德,孙国强,等.基于碳交易机制的电-气互联综合能源系统低碳经济运行[J].电力系统自动化,2016,40(15):9-16.
Wei Zhinong,Zhang Side,Sun Guoqiang,et al.Carbon Trading Based Low-carbon Economic Operation for Integrated Electricity and Natural Gas EnergySystem[J].Automation of Electric Power Systems,2016,40(15):9-16.
[2]秦婷,刘怀东,王锦桥,等.基于碳交易的电-热-气综合能源系统低碳经济调度[J].电力系统自动化,2018,42(14):8-13+22.
Qin Ting,Liu Huaidong,Wang Jinqiao,et al.Carbon Trading Based Low-carbon Economic Dispatch for Integrated Electricity-Heat-Gas Energy System[J].Automation of Electric Power Systems,2018,42(14):8-13+22.
[3]卢志刚,杨宇,耿丽君,等.基于Benders分解法的电热综合能源系统低碳经济调度[J].中国电机工程学报,2018,38(7):1922-1934.
Lu Zhigang,Yang Yu,Geng Lijun,et al.Low-carbon Economic Dispatch ofthe Integrated Electrical and Heating Systems Based on Benders Decomposition[J].Proceedings of the CSEE,2018,38(7):1922-1934.
[4]周任军,肖钧文,唐夏菲,等.电转气消纳新能源与碳捕集电厂碳利用的协调优化[J].电力自动化设备,2018,38(07):61-67.
Zhou Renjun,Xiao Junwen,Tang Xiafei,et al.Coordinated Optimization ofCarbon Utilization between Power-to-gas Renewable Energy Accommodation andCarbon Capture Power Plant[J].Electric Power Automation Equipment,2018,38(07):61-67.
[5]He L,Lu Z,Zhang J,et al.Low-carbon economic dispatch forelectricity and natural gas systems considering carbon capture systems andpower-to-gas[J].Applied energy,2018,224:357-370.
本发明实施例中的计算条件、图例等仅用于对本发明作进一步的说明,并非穷举,并不构成对权利要求保护范围的限定,本领域技术人员根据本发明实施例获得的启示,不经过创造性劳动就能够想到其它实质上等同的替代,均在本发明保护范围内。
Claims (1)
1.一种考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法,其特征是,它通过考虑富氧燃烧电厂的运行特点及电-气-热综合能源系统的协调配合,构建了电-气-热综合能源系统协调调度模型,其具体包括以下步骤:
步骤一,富氧燃烧电厂电能时移特性分析:
富氧燃烧电厂的电能损耗主要集中于空分制氧装置和碳捕集压缩纯化装置;为应对快速变负荷要求,空分制氧装置配有储存液态氧的隔热储氧罐;富氧燃烧电厂的电能主要流向电负荷、空分制氧装置和碳捕集压缩纯化装置,正常运行情况下,电负荷用电能最多,空分制氧装置其次,碳捕集压缩纯化装置最少;当电负荷过多时通过降低空分制氧装置和碳捕集压缩纯化装置的功率为电负荷提供支撑,而富氧燃烧过程所需氧气量和燃料量成正比,此时空分制氧装置不能产生足够的氧气,但缺少的氧气能从储氧罐中获取;因此,富氧燃烧电厂能够利用储氧罐在多时段进行电能调节,在电负荷低谷时让空分制氧装置满功率运行,在电负荷高峰时减少空分制氧装置的出力,其中的氧气差值利用储氧罐进行调整,实现电能的时移;
步骤二,富氧燃烧电厂净输出功率模型构建:
富氧燃烧电厂产生的电能主要流向空分制氧装置、碳捕集压缩纯化装置与电负荷,用(1)式表示:
PZ=PN+PCPU+PASU (1)
式中,PZ为火电厂的输出总功率;PN为火电厂的净输出功率;PCPU为碳捕集压缩纯化装置所需功率;PASU为空分制氧装置所需功率;
设常规火电厂的碳排放强度为EG,其燃烧化石能源所产生的CO2总量EZ为:
EZ=EGPZ (2)
若富氧燃烧电厂的碳捕集水平为α,则所能捕获的CO2量EB为:
EB=αEZ (3)
此时,空分制氧装置所需功率PASU和碳捕集压缩纯化装置所需功率PCPU为(4)式:
式中,OC为空分制氧装置产生氧气量;OCQ为储氧罐供氧量;OZ为进行富氧燃烧所需氧气量;β为空分制氧装置单位时间制造单位氧气所需电能;δ为碳捕集压缩纯化装置单位时间捕捉单位CO2所需电能;χ为富氧燃烧电厂进行富氧燃烧产生单位功率所需氧气量;
从而推出净输出功率为(5)式:
PN=(1-αδEG-βχ)PZ+βOCQ (5)
从中得知富氧燃烧电厂净输出功率和碳捕集水平、空分制氧装置的制氧量、储气罐供氧量有关,改变储氧罐的充放氧能够改变某一时段的净输出功率,因此富氧燃烧能够利用储氧罐做到电能时移;
富氧燃烧电厂存在两种运行模式,一是富氧运行,二是空气运行,这两种模式能够相互切换,另外在空气燃烧的过程中,碳捕集压缩纯化装置不进行捕碳,但这并不影响空分制氧装置的制氧,其能将产生的氧气装入储氧罐中;
步骤三,富氧燃烧电厂的电碳特性分析:
考虑空分制氧装置是否投运对富氧燃烧电厂运行的影响,分别对空分制氧装置投运时和空分制氧装置停运时的电碳特性进行分析,电碳特性指电厂净输出功率和二氧化碳净排放量的关系;
1)空分制氧装置投运
在空分制氧装置运行时,由于存在储氧罐,空分制氧装置能够在限定内的功率区间运行,即在火电厂最小运行功率时,能够使空分制氧装置以最大功率运行,在火电厂最大运行功率时,使空分制氧装置以最小功率运行,因此火电厂具有更深的调峰深度;而碳捕集压缩纯化装置能够根据当前火电厂输出总功率调节碳捕集电能;净碳排放量为火电厂实际排放量减去捕碳量,则空分制氧装置运行时电碳特性关系为(6)式:
式中,EJP为净排碳量;αmax为最大碳捕集水平;
2)空分制氧装置停运
在空分制氧装置停运时,由于存在储氧罐,富氧燃烧电厂依旧能工作在碳捕集状态,则空分制氧装置停运时电碳特性关系为(7)式:
当储氧罐的储氧量不足时,富氧燃烧电厂需退出富氧运行状态,转为空气运行状态;
步骤四,电-气-热综合能源系统低碳经济调度模型构建:
电-气-热综合能源系统,其中,电负荷、电转气设备的用电由风电场、火电厂、燃气轮机进行电能的供给,热负荷的用热由燃气轮机、燃气锅炉、储热罐进行热能供给,气负荷、燃气轮机、燃气锅炉的用气由电转气设备、气源、储气罐进行气能供给;
利用储热罐协调燃气轮机出力支撑空分制氧装置和碳捕集压缩纯化装置的运行;燃气轮机是以天然气为原材料,通过燃烧天然气使空气膨胀推动汽轮机转动从而获取电能,其高温烟气通过余热锅炉能够获取热能,因此会出现以热定电的情况;通过加入储热罐能在一定程度打破以热定电,从而使燃气轮机能在电负荷峰值时提供更多的电能,减小火电厂净出力,增加对空分制氧装置和碳捕集压缩纯化装置的电能输送,提高碳捕集水平;
在电-气-热综合能源系统中,电转气设备能够将弃风接纳转换成天然气,同时能接收一部分富氧燃烧电厂捕捉到的CO2,从而节约封存成本和购买高纯度CO2成本,使经济性得到改善;
以电-气-热综合能源系统运行成本、碳捕集设备折旧成本以及维护成本最小为目标,考虑各机组、元件的约束条件构建低碳经济调度模型;其中电-气-热综合能源系统运行成本包括碳交易成本、火电燃料成本、外购天然气成本、运输封存CO2成本、外购CO2成本,而碳捕集设备包括燃烧后捕集电厂中的吸收塔、再生塔、压缩机以及富氧燃烧电厂中的碳捕集压缩纯化装置、空分制氧装置,目标函数如下:
minF=min(CT+CH+CQ+CW+CYF+CMCO2+CTB)
式中,F为电-气-热综合能源系统总成本;CT为电-气-热综合能源系统碳交易成本;CH为电-气-热综合能源系统火电燃料成本;CQ为外购天然气成本;CW为火电厂维护成本;CYF为运输封存CO2成本;CMCO2为外购CO2成本,指购买高纯度CO2供给电转气设备的成本;CTB为碳捕集设备折旧成本;σ为单位碳排放权交易价格;Ec为电-气-热综合能源系统的实际碳排放量;N1为火电厂的个数;N2为燃气机组的个数;Dk1为第k1个火电厂的碳排量配额;Dk2为第k2个燃气机组的碳排放配额;Pk1,t为第k1个火电厂t时段的输出电功率;ak1、bk1、ck1为火电厂的成本系数;μgas为购买天然气价格;GRZ为天然气热值;Zer为单位电能所能转化的热能;Rrg,k2,t为第k2个燃气机组为燃气锅炉时t时段输出的热功率;Pk2,t为第k2个燃气机组为燃气轮机时t时段输出的电功率;ηrg为燃气锅炉的热效率;ηrq为燃气轮机的电效率;QP2G,t为t时段电转气设备产生的天然气量;CCG为火电厂成本;ρCG为火电厂维护系数;CTL为脱硫装置成本;ρTL为脱硫装置维护系数;CTX为脱硝装置成本;ρTX为脱硝装置维护系数;CGZ为富氧燃烧改造锅炉成本;CASU为空分制氧装置成本;ρASU为空分制氧装置维护系数;CCPU为碳捕集压缩纯化装置成本;ρCPU为碳捕集压缩纯化装置维护系数;CRHGZ为燃烧后捕集电厂锅炉改造成本;CTBJ为燃烧后捕集电厂中吸收塔、再生塔、压缩机的成本;ρTB为燃烧后捕集电厂中吸收塔、再生塔、压缩机的维护系数;μYF为运输封存单位CO2的成本;BCO2为捕捉到的CO2总量;BP2G为捕捉到的CO2供电转气设备用掉的量;μCO2为外购单位CO2所需价格;θCO2为输入电转气设备单位电功率所需CO2量;PP2G,t为t时段输入电转气设备的电功率;NT为富氧燃烧电厂中空分制氧装置、碳捕集压缩纯化装置和燃烧后捕集电厂中吸收塔、再生塔、压缩机的折旧年限;ω为折旧率;PCY为单位体积储氧罐价格;VCY为储氧罐体积;NC为储氧罐折旧年限;
约束条件包括功率平衡约束、富氧燃烧电厂约束、电转气设备约束、风电场约束、燃气轮机约束、燃气锅炉约束、储氧\热\气罐约束;
1)功率平衡约束
对于电负荷、热负荷、气负荷分别满足功率平衡约束,如下:
式中,Pel,t为t时段的电负荷;Pk1J,t为t时段k1火电厂净输出功率;Pk2,t为t时段k2燃气机组为燃气轮机时所供电功率;Pwind,t为风电在t时段的输出电功率;PP2G,t为t时段输入电转气设备的电功率;Rrl,t为t时段的热负荷;Rrg,k2,t为t时段k2燃气机组为燃气锅炉时所供热功率;Rrq,k2,t为t时段k2燃气机组为燃气轮机时所供热功率;σrx为储热罐充放热效率;RXr,t为t时段储热罐充热量;RFr,t为t时段储热罐放热量;Qql,t为t时段的气负荷;QQY,t为t时段的气源出力;QXQ,t为t时段的储气罐充气量;QFQ,t为t时段的储气罐放气量;QP2G,t为t时段的电转气设备产气量;Qrg,k2,t为t时段k2燃气机组为燃气锅炉时所需天然气量;Qrq,k2,t为t时段k2燃气机组为燃气轮机时所需天然气量;
2)富氧燃烧电厂约束
富氧燃烧电厂需要满足火电平衡约束、空分制氧装置制氧约束、空分制氧装置出力约束、碳捕集压缩纯化装置工作约束、捕碳量约束、碳捕集水平约束、净排碳量约束、氧气平衡约束、富氧燃烧所需氧气量约束、火电出力限值约束、火电出力爬坡约束为(10)式:
式中,Pkl,t为t时段k1火电厂的总电功率;PASUk1,t为k1火电厂在t时段的空分制氧装置所需功率;PCPUk1,t为k1火电厂在t时段碳捕集压缩纯化装置的所需功率;β为空分制氧装置单位时间制造单位氧气所需电能;OCk1,t为k1火电厂的空分制氧装置在t时段产生氧气量;λASU为空分制氧装置最小运行系数;PASUk1max为火电厂k1的空分制氧装置最大功率;δ为碳捕集压缩纯化装置单位时间捕捉单位CO2所需电能;EBk1,t为k1火电厂在t时段的捕碳量;ζ为0-1变量,ζ=1表示富氧燃烧电厂处于富氧燃烧状态,ζ=0表示富氧燃烧电厂不处于富氧燃烧状态;α为碳捕集水平;EGk1为k1火电厂的碳排放强度;αmax为最大碳捕集水平;EJPk1,t为k1火电厂在t时段的净排碳量;OZk1,t为k1火电厂在t时段进行富氧燃烧所需氧气量;OFQk1,t为火电厂k1的储氧罐t时段放氧量;OXQk1,t为火电厂k1的储氧罐t时段充氧量;OQYLk1,t为t时段k1火电厂允许弃氧量;χ为富氧燃烧电厂进行富氧燃烧产生单位功率所需氧气量;Pk1min为火电厂k1最小出力值;Pk1max为火电厂k1最大出力值;Pk1Xmin为火电厂k1最大下爬坡率;Pk1Smax为火电厂k1最大上爬坡率;
3)电转气设备约束
电转气设备主要满足其最大功率约束为:
0≤PP2G,t≤PP2Gmax (11)
式中,PP2Gmax为电转气设备运行的最大功率;
4)风电场约束
风电场主要满足其最大功率约束为:
0≤Pwind,t≤PW,t (12)
式中,Pw,t为风电预测功率;
5)燃气轮机约束
燃气轮机主要满足出力上下限约束、爬坡率约束和热电比约束为:
Pk2min≤Pk2,t≤Pk2max (13)
Pk2Xmin≤Pk2,t-Pk2,t-1≤Pk2Smax (14)
式中,Pk2min为第k2个燃气机组为燃气轮机时的最小电功率;Pk2max为第k2个燃气机组为燃气轮机时的最大电功率;Pk2Xmin为第k2个燃气机组为燃气轮机时的最大下爬坡率;Pk2Smax为第k2个燃气机组为燃气轮机时的最大上爬坡率;ΨREC为燃气轮机的热电比;
6)燃气锅炉约束
燃气锅炉主要满足出力上下限约束和爬坡率约束为:
Rrg,k2min≤Rrg,k2,t≤Rrg,k2max (16)
Rrg,k2Xmin≤Rrg,k2,t-Rrg,k2,t-1≤Rrg,k2Smax (17)
式中,Rrg,k2min为第k2个燃气机组为燃气锅炉时最小热功率;Rrg,k2max为第k2个燃气机组为燃气锅炉时最大热功率;Rrg,k2Xmin为第k2个燃气机组为燃气锅炉时的最大下爬坡率;Rrg,k2Smax为第k2个燃气机组为燃气锅炉时的最大上爬坡率;
7)储氧/气/热罐约束
约束包括储氧/气/热容量上下限约束,储氧/气/热周期始末相等约束,充放氧/气/热限值约束,储氧/气/热平衡约束为:
0≤Echu,t≤Echumax (18)
Echu,1=Echu,24 (19)
0≤EX,t≤ωEXmax (20)
0≤EF,t≤(1-ω)EFmax (21)
Echu,t=Echu,t-1+EX,t-EF,t (22)
式中,Echu,t为储氧/气/热罐在t时段的储氧/气/热量;Echumax为储氧/气/热罐容量上限;Echu,1为储氧/气/热罐在周期内初始储氧/气/热量;Echu,24为储氧/气/热罐在周期结束时储氧/气/热量;EX,t为t时段的充氧/气/热量;EF,t为t时段的放氧/气/热量;EXmax和EFmax分别为单位时段充氧/气/热和放氧/气/热量上限;ω为二进制变量,当ω=1时进行充氧/气/热,当ω=0时进行放氧/气/热。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010281250.2A CN111489083B (zh) | 2020-04-11 | 2020-04-11 | 考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010281250.2A CN111489083B (zh) | 2020-04-11 | 2020-04-11 | 考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111489083A true CN111489083A (zh) | 2020-08-04 |
CN111489083B CN111489083B (zh) | 2022-03-18 |
Family
ID=71811739
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010281250.2A Active CN111489083B (zh) | 2020-04-11 | 2020-04-11 | 考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111489083B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112417652A (zh) * | 2020-10-30 | 2021-02-26 | 东北电力大学 | 一种电-气-热综合能源系统优化调度方法及系统 |
CN113178865A (zh) * | 2021-04-23 | 2021-07-27 | 东北电力大学 | 基于碳氧环流的能源集线器及其优化调度方法 |
CN113327180A (zh) * | 2021-07-05 | 2021-08-31 | 华北电力大学 | 一种考虑氢能应用的电力系统低碳经济调度方法及系统 |
CN113627011A (zh) * | 2021-08-06 | 2021-11-09 | 深圳供电局有限公司 | 一种基于多目标优化的综合能源运行仿真方法 |
CN113705991A (zh) * | 2021-08-17 | 2021-11-26 | 国网四川省电力公司技能培训中心 | 一种多能源园区的建立及低碳调度方法 |
CN114033638A (zh) * | 2021-11-30 | 2022-02-11 | 东北电力大学 | 一种电-热-气-肥联产系统及其协调控制方法 |
CN116499272A (zh) * | 2023-06-19 | 2023-07-28 | 广东中鹏热能科技有限公司 | 一种多能源燃烧碳排放追踪的智能窑炉控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103106621A (zh) * | 2013-01-18 | 2013-05-15 | 长沙理工大学 | 碳排放权交易下碳捕集机组的厂内优化运行方法 |
CN104268668A (zh) * | 2014-09-09 | 2015-01-07 | 国家电网公司 | 一种闭环低碳优化调度方法 |
CN104595909A (zh) * | 2015-01-29 | 2015-05-06 | 威海鲁源科技环保设备有限公司 | 一种医疗垃圾焚烧净化装置 |
CN106650219A (zh) * | 2016-10-20 | 2017-05-10 | 广东电网有限责任公司电力科学研究院 | 一种富氧燃烧煤粉锅炉热平衡计算方法 |
CN108229865A (zh) * | 2018-03-30 | 2018-06-29 | 天津大学 | 一种基于碳交易的电热气综合能源系统低碳经济调度方法 |
CN108625990A (zh) * | 2018-04-13 | 2018-10-09 | 华中科技大学 | 一种天然气富氧燃烧与跨临界co2循环耦合的发电系统 |
-
2020
- 2020-04-11 CN CN202010281250.2A patent/CN111489083B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103106621A (zh) * | 2013-01-18 | 2013-05-15 | 长沙理工大学 | 碳排放权交易下碳捕集机组的厂内优化运行方法 |
CN104268668A (zh) * | 2014-09-09 | 2015-01-07 | 国家电网公司 | 一种闭环低碳优化调度方法 |
CN104595909A (zh) * | 2015-01-29 | 2015-05-06 | 威海鲁源科技环保设备有限公司 | 一种医疗垃圾焚烧净化装置 |
CN106650219A (zh) * | 2016-10-20 | 2017-05-10 | 广东电网有限责任公司电力科学研究院 | 一种富氧燃烧煤粉锅炉热平衡计算方法 |
CN108229865A (zh) * | 2018-03-30 | 2018-06-29 | 天津大学 | 一种基于碳交易的电热气综合能源系统低碳经济调度方法 |
CN108625990A (zh) * | 2018-04-13 | 2018-10-09 | 华中科技大学 | 一种天然气富氧燃烧与跨临界co2循环耦合的发电系统 |
Non-Patent Citations (2)
Title |
---|
崔杨等: "计及综合成本的风电-光伏-光热联合出力调度策略", 《高电压技术》 * |
曲艳伟等: "基于富氧燃烧CCS技术对我国低碳发展的研究", 《锅炉制造》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112417652A (zh) * | 2020-10-30 | 2021-02-26 | 东北电力大学 | 一种电-气-热综合能源系统优化调度方法及系统 |
CN113178865A (zh) * | 2021-04-23 | 2021-07-27 | 东北电力大学 | 基于碳氧环流的能源集线器及其优化调度方法 |
CN113178865B (zh) * | 2021-04-23 | 2022-05-24 | 东北电力大学 | 基于碳氧环流的能源集线器及其优化调度方法 |
CN113327180A (zh) * | 2021-07-05 | 2021-08-31 | 华北电力大学 | 一种考虑氢能应用的电力系统低碳经济调度方法及系统 |
CN113327180B (zh) * | 2021-07-05 | 2023-09-29 | 华北电力大学 | 一种考虑氢能应用的电力系统低碳经济调度方法及系统 |
CN113627011A (zh) * | 2021-08-06 | 2021-11-09 | 深圳供电局有限公司 | 一种基于多目标优化的综合能源运行仿真方法 |
CN113705991A (zh) * | 2021-08-17 | 2021-11-26 | 国网四川省电力公司技能培训中心 | 一种多能源园区的建立及低碳调度方法 |
CN113705991B (zh) * | 2021-08-17 | 2023-08-29 | 国网四川省电力公司技能培训中心 | 一种多能源园区的建立及低碳调度方法 |
CN114033638A (zh) * | 2021-11-30 | 2022-02-11 | 东北电力大学 | 一种电-热-气-肥联产系统及其协调控制方法 |
CN114033638B (zh) * | 2021-11-30 | 2023-11-24 | 东北电力大学 | 一种电-热-气-肥联产系统及其协调控制方法 |
CN116499272A (zh) * | 2023-06-19 | 2023-07-28 | 广东中鹏热能科技有限公司 | 一种多能源燃烧碳排放追踪的智能窑炉控制方法 |
CN116499272B (zh) * | 2023-06-19 | 2024-01-16 | 广东中鹏热能科技有限公司 | 一种多能源燃烧碳排放追踪的智能窑炉控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111489083B (zh) | 2022-03-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111489083B (zh) | 考虑富氧燃烧技术的电-气-热综合能源系统低碳经济调度方法 | |
CN110854930A (zh) | 一种考虑多能协同的园区综合能源系统优化调度方法 | |
EP3155237A1 (en) | System and method for supplying an energy grid with energy from an intermittent renewable energy source | |
CN110543157A (zh) | 一种多能互补智慧供应热电氢的系统及方法 | |
EP3155238A1 (en) | System and method for supplying an energy grid with energy from an intermittent renewable energy source | |
CN113794227B (zh) | 一种实现碳循环的综合能源系统 | |
CN113451612A (zh) | 一种绿色高效的电力-氨-电力能源系统 | |
CN110957722B (zh) | 一种含电转气设备的微型能源网日前优化调度方法 | |
CN114462889A (zh) | 一种氢-电耦合多能源跨区域优化配置方法及系统 | |
CN115796323A (zh) | 含ccus-p2g-chp的综合能源系统低碳经济调度方法 | |
CN210199571U (zh) | 一种多能互补智慧供应热电氢的系统 | |
CN116777054A (zh) | 考虑碳捕集与垃圾焚烧参与的虚拟电厂及其优化运行方法 | |
CN217362587U (zh) | 一种新能源驱动的富碳可再生燃烧循环系统 | |
CN113972662B (zh) | 综合能源生产单元及方法 | |
CN116316888A (zh) | 氢电耦合系统的优化调度方法、系统及装置 | |
CN115660142A (zh) | 一种园区综合能源系统源荷储协调优化调度方法 | |
CN113178865B (zh) | 基于碳氧环流的能源集线器及其优化调度方法 | |
Li et al. | Low-carbon economic dispatch of power system based on mobile hydrogen storage | |
CN219792508U (zh) | 一种基于源网荷储氢一体化的工业硅低碳制造系统 | |
CN114362152B (zh) | 综合能源系统多时间尺度调度方法 | |
Mu et al. | An Operation Scheduling Model for Carbon Neutrality in Industrial Integrated Energy System | |
CN115619006B (zh) | 一种计及辅助服务的电-气-氢混联综合能源系统优化调度方法 | |
CN117216963B (zh) | 一种含碳-氢利用的综合能源系统运行方法 | |
CN117910722A (zh) | 基于源荷不确定性的区域多能互补能源系统优化调度方法 | |
CN117638862A (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 |