CN106650068B - 一种预测燃煤电厂碳排放量的计算方法 - Google Patents
一种预测燃煤电厂碳排放量的计算方法 Download PDFInfo
- Publication number
- CN106650068B CN106650068B CN201611137401.7A CN201611137401A CN106650068B CN 106650068 B CN106650068 B CN 106650068B CN 201611137401 A CN201611137401 A CN 201611137401A CN 106650068 B CN106650068 B CN 106650068B
- Authority
- CN
- China
- Prior art keywords
- coal
- power plant
- flue gas
- calculating
- 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.)
- Active
Links
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 40
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 34
- 238000004364 calculation method Methods 0.000 title claims abstract description 31
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims abstract description 80
- 239000003245 coal Substances 0.000 claims abstract description 48
- 239000000446 fuel Substances 0.000 claims abstract description 47
- 229910002092 carbon dioxide Inorganic materials 0.000 claims abstract description 44
- 239000003546 flue gas Substances 0.000 claims abstract description 43
- 239000001569 carbon dioxide Substances 0.000 claims abstract description 40
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims abstract description 32
- 238000013178 mathematical model Methods 0.000 claims abstract description 7
- 238000004458 analytical method Methods 0.000 claims description 26
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 claims description 22
- 238000012417 linear regression Methods 0.000 claims description 18
- 238000000034 method Methods 0.000 claims description 17
- 229910052760 oxygen Inorganic materials 0.000 claims description 16
- RHZUVFJBSILHOK-UHFFFAOYSA-N anthracen-1-ylmethanolate Chemical compound C1=CC=C2C=C3C(C[O-])=CC=CC3=CC2=C1 RHZUVFJBSILHOK-UHFFFAOYSA-N 0.000 claims description 14
- 239000003830 anthracite Substances 0.000 claims description 14
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 14
- 239000001301 oxygen Substances 0.000 claims description 14
- 239000003077 lignite Substances 0.000 claims description 12
- 238000002485 combustion reaction Methods 0.000 claims description 8
- 239000007789 gas Substances 0.000 claims description 8
- 230000001419 dependent effect Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 239000000779 smoke Substances 0.000 claims description 5
- 238000000611 regression analysis Methods 0.000 claims description 3
- 238000007605 air drying Methods 0.000 claims description 2
- 238000000491 multivariate analysis Methods 0.000 claims 1
- 239000003039 volatile agent Substances 0.000 claims 1
- 238000001514 detection method Methods 0.000 description 8
- 229910002091 carbon monoxide Inorganic materials 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000003915 air pollution Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 239000002893 slag Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/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
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/84—Greenhouse gas [GHG] management systems
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Computer Hardware Design (AREA)
- General Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Regulation And Control Of Combustion (AREA)
Abstract
本发明公开了一种预测燃煤电厂碳排放量的计算方法,该方法包括下列步骤:S1、建立数学模型,得出燃料特性系数β的关系式,S2、计算电厂燃烧的煤的燃料特性系数β值,S3、计算烟气中二氧化碳的百分比含量CO2,S4、计算烟气的排放总量Vy,S5、计算二氧化碳的排放总量本发明利用Matlab软件进行多元逐步线性回归拟合出的燃料特性系数β的线性方程和根据检测到的尾气中氧气和二氧化硫的浓度来计算二氧化碳排放量的方法计算结果精确可靠,更有利于碳排放的监测和控制。该计算方法所采用数据均是常规燃煤电厂的基本检测数据,不需要价格昂贵的计量设备或检测手段来检测二氧化碳排放量,降低了检测成本,实用性较强。
Description
技术领域
本发明涉及碳排放量计算技术领域,具体涉及一种预测燃煤电厂碳排放量的计算方法。
背景技术
2010年10月中央政府发布的“十二五”规划中已明确提出了实现低碳发展的核心内容是向低碳经济转型,调整优化以低碳经济为核心的产业结构和能源结构,建立资源环境友好型社会,坚持大力推广节能技术,制定围绕低碳经济的消费模式和方针政策。而电力行业在能源转换过程中排放的二氧化碳约占全国排放总量的50%,是二氧化碳减排的关键部门之一。截止到2014年底,中国全国单位国内生产总值二氧化碳排放同比下降了6.1%,比2010年累计下降了15.8%,“十二五”规划要求下降17%的目标已经完成。
在“十三五”规划中,能源局明确提出我国要完成到2020年单位GDP碳排放要比2005年下降40%~45%的国际承诺低碳目标,并且要为完成中美气候变化联合声明中提出的我国在2030年左右要达到碳排放的峰值的中长期低碳发展目标奠定基础,同时要在大气污染防治等环境指标方面取得明显成效。要完成上述目标,需要各省市结合经济发展实际,结合本地特点,制定出科学的碳减排规划和对策,有步骤、有计划的完成碳减排指标。由于我国的水资源和耕地资源的缺乏,我国以煤为主的能源消耗结构将会持续相当长的一段时间,所以为了更好的实现碳排放减排,需要比较精确的碳排放量。
目前,我国燃煤电厂主要是通过《2006IPCC指南》中提供的方法一来计算CO2排放量,即根据燃料特性直接获得碳的缺省排放因子,根据电厂实际原煤耗量,乘上各煤种的排放因子,从而估算二氧化碳排放量。实际上,宽范围的煤种排放因子对CO2预测误差较大,同时影响CO2排放因子的主要影响因素除燃料类型外,还有机组装机容量、机组使用年限与维护质量以及燃料燃烧方式等。不同类型机组由于发电热效率不同CO2排放绩效不同;相同容量机组由于使用年限和维护质量的不同,CO2排放因子也会有差异甚至产生较大差异。此外,由于燃烧方式不同,炉渣中未完全燃烧的碳含量不同,CO2排放因子也会产生一定差异。因此,利用缺省因子只能得到相当粗糙的结果,从而引入很大的不确定性。
目前为获得较精确的二氧化碳排放,各研究者大多采用了一些基于碳平衡的精确数学模型,通过对原煤的煤质分析进行预测。而这种基于碳平衡的二氧化碳排放预测结果虽然较精确,但过程却较复杂,成本相对较高。由于通过烟气侧进行计算也是一种途径,考虑到电厂一般具有尾气氧气浓度和二氧化硫浓度的检测,以及燃煤的工业分析和热值这些基本数据,本发明拟采用烟气排放成分的方法,尽量利用电厂现有数据,在不增加检测的情况下,较精确的预测二氧化碳排放量。
发明内容
本发明的目的是为了解决现有技术中的上述缺陷,提供一种预测燃煤电厂碳排放量的计算方法,该计算方法根据电厂烟气排放量和燃煤的燃料特性系数来预测二氧化碳排放量,是一种能够适应我国燃煤电厂实际情况的碳排放量的计算方法。
本发明的目的可以通过采取如下技术方案达到:
一种预测燃煤电厂碳排放量的计算方法,所述计算方法包括下列步骤:
S1、建立数学模型,得出燃料特性系数β的关系式:
根据现有的煤质分析数据,利用多元线性回归模型拟合出燃煤的工业分析数据Aad、Vad、FCad、Qgr,ad与燃料特性系数β之间的关系式,其中Aad、Vad、FCad、Qgr,ad分别表示为收到基下燃煤中的灰分含量(单位:%)、挥发分含量(单位:%)、固定碳含量(单位:%)和燃煤的高位发热量(单位:MJ/kg),得到褐煤和无烟煤关于燃料特性系数和工业分析数据之间的线性回归方程;
S2、计算电厂燃烧的煤的燃料特性系数β值:
将电厂测得的燃煤工业分析数据Aad、Vad、FCad、Qgr,ad根据煤种的不同带入所述的褐煤或无烟煤关于燃料特性系数和工业分析数据之间的线性回归方程中,计算出电厂燃煤的燃料特性系数β值;
S3、计算烟气中二氧化碳的百分比含量CO2:
将电厂检测的尾气中氧气和二氧化硫的百分比含量和所述的计算出电厂燃煤的燃料特性系数β值带入二氧化碳含量计算式计算出烟气中二氧化碳的百分比含量;
S4、计算烟气的排放总量Vy;
S5、计算二氧化碳的排放总量
将所述的烟气中二氧化碳的百分比含量CO2和所述的烟气的排放总量Vy带入公式即可计算得到电厂中CO2的实际排放量。
进一步地,所述的步骤S1、建立数学模型,得出燃料特性系数β的关系式具体如下:
S11、将无烟煤的工业分析数据和煤的燃料特性系数作为样本数据,并将其带入仿真软件的线性回归分析工具中,则X和Y分别表示为:
式中,Ai、Vi、FCi、Qi分别表示为第i个样本中的空气干燥基灰分(单位:%)、挥发分(单位:%)、固定碳(单位:%)和高位发热量(单位:MJ/kg),βi为第i个样本的燃料特性系数;
S12、根据仿真软件分别对褐煤和无烟煤样本数据进行逐步回归分析,得到褐煤和无烟煤关于燃料特性系数和工业分析数据之间的线性回归方程分别如下所示,
进一步地,所述仿真软件为MATLAB,所述线性回归分析工具为stepwise命令,其中,stepwise命令的格式为:stepwise(X,Y),其中X表示自变量数据,为n×m阶矩阵,m表示自变量种类数,n表示每种自变量的样本数,Y表示因变量数据,为n×1阶矩阵,n表示因变量样本数。
进一步地,所述的二氧化碳含量计算式如下:
式中,CO2表示烟气中二氧化碳的体积分数,O2表示烟气中氧气的体积分数,SO2表示烟气中二氧化硫的体积分数。
进一步地,所述的步骤S4、计算烟气的排放总量Vy具体如下:
若电厂已经测得烟气排放的总量Vy,则直接使用该数据;
若电厂未测得烟气排放的总量Vy,则可根据燃煤量和过量空气系数进行计算,如下所示:
其中,M表示燃烧所用的燃料量(单位:t),Qnet,ar表示燃料的低位发热量(单位:MJ/kg),α表示空气过量系数,V0表示每燃烧1吨煤所需的理论空气量,其计算式为:
V0=(1.05×Qnet,ar+0.278)×1000
其中,O2表示已检测得到的烟气中氧气的体积分数。
本发明相对于现有技术具有如下的优点及效果:
(1)相比我国目前较为通用的利用IPCC指南中的缺省因子进行碳排放量计算的方法,本发明利用Matlab软件进行多元逐步线性回归拟合出的燃料特性系数β的线性方程和根据检测到的尾气中氧气和二氧化硫的浓度来计算二氧化碳排放量的方法更加精确,结果更为可靠,相对而言更有利于碳排放的监测和控制。
(2)本发明公开的碳排放量的计算方法所采用数据均是常规燃煤电厂的基本检测数据,计算简单方便,因此并不需要价格昂贵的计量设备或检测手段来检测二氧化碳排放量,与其他一些基于碳平衡的计算方法相比,大大降低了检测成本,实用性较强,具备推广条件。
附图说明
图1是本发明中公开的预测燃煤电厂碳排放量的计算方法的流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例
本实施例利用电厂烟气成分来计算二氧化碳排放量,其计算原理为:
式中,CO2表示烟气中二氧化碳的体积分数,Vy为烟气排放总量。
根据燃烧方程,可以得到燃料特性系数和烟气成分的关系如下所示:
21=RO2+O2+0.605CO+β(RO2+CO) (2)
其中,RO2—烟气中三原子气体的体积分数,主要有二氧化硫和二氧化碳;O2—烟气中氧气的体积分数;CO—烟气中一氧化碳的体积分数;β—燃料的特性系数。
由于实际燃煤电厂烟气中一氧化碳浓度十分低,并且在管道中高温烟气的作用下,部分一氧化碳会被继续氧化成二氧化碳,因此可忽略烟气中的一氧化碳含量,即此时的二氧化碳含量计算式为:
电厂一般具有尾气氧气和二氧化硫浓度的检测,因此由原理(3)可知,在已知烟气中氧气和二氧化硫含量的条件下,只需要得到燃料特性系数β就可以得到烟气中二氧化碳的体积分数,其具体实施步骤包括:
S1、建立数学模型,得出燃料特性系数β的关系式:
根据大量的现有的煤质分析数据,利用多元线性回归模型拟合出燃煤的工业分析数据Aad、Vad、FCad、Qgr,ad与燃料特性系数β之间的关系式,其中Aad、Vad、FCad、Qgr,ad分别表示为收到基下燃煤中的灰分含量(单位:%)、挥发分含量(单位:%)、固定碳含量(单位:%)和燃煤的高位发热量(单位:MJ/kg)。具体实施方法如下所示:
在MATLAB中,stepwise命令的格式为:stepwise(X,Y),其中X表示自变量数据,为n×m阶矩阵(m表示自变量种类数,这里为4;n表示每种自变量的样本数),Y表示因变量数据,为n×1阶矩阵(n表示因变量样本数)。在《中国煤炭性质、分类与利用》上分别选取51组褐煤和111组无烟煤的工业分析数据和煤的燃料特性系数作为样本数据,并将其带入stepwise命令中,则X和Y分别表示为:
式中,Ai、Vi、FCi、Qi分别表示为第i个样本中的空气干燥基灰分(单位:%)、挥发分(单位:%)、固定碳(单位:%)和高位发热量(单位:MJ/kg);βi为第i个样本的燃料特性系数。再根据MATLAB分别对褐煤和无烟煤样本数据进行逐步回归分析,并经过进一步调试,最终可得到褐煤和无烟煤关于燃料特性系数和工业分析数据之间的线性回归方程分别如下所示,
其误差范围均在-2.5%—+2.5%之间,都在可接受的合理范围内。
S2、计算电厂燃烧的煤的燃料特性系数β值:
将电厂测得的燃煤工业分析数据Aad、Vad、FCad、Qgr,ad根据煤种的不同带入上述步骤1中利用线性回归方法拟合得到的相应的褐煤或无烟煤关于燃料特性系数和工业分析数据之间的线性回归方程中,计算出电厂燃煤的燃料特性系数β值。
S3、计算烟气中二氧化碳的百分比含量CO2:
将电厂检测的尾气中氧气和二氧化硫的百分比含量和上述步骤2所求的燃煤的燃料特性系数β值一起带入式(3)即可计算出烟气中二氧化碳的百分比含量。
S4、计算烟气的排放总量Vy:
若电厂已经测得烟气排放的总量Vy,则可直接使用该数据;
若电厂未测得烟气排放总量,则可根据燃煤量和过量空气系数进行计算,如下所示:
其中,M—表示燃烧所用的燃料量,t;Qnet,ar—表示燃料的低位发热量,MJ/kg;α—表示空气过量系数;V0—表示每燃烧1吨煤所需的理论空气量,其计算式为:
V0=(1.05×Qnet,ar+0.278)×1000 (8)
其中,O2—表示已检测得到的烟气中氧气的体积分数。
S5、计算二氧化碳的排放总量
将上述步骤3中所计算得到的烟气中的CO2百分比含量和步骤4中计算得到的烟气排放总量Vy带入式(1)即可计算得到电厂中CO2的实际排放量。
如上所述,便可较好地实现本发明。
综上所述,本发明利用Matlab软件进行多元逐步线性回归拟合出的燃料特性系数β的线性方程和根据检测到的尾气中氧气和二氧化硫的浓度来计算二氧化碳排放量的方法计算结果精确可靠,更有利于碳排放的监测和控制。该计算方法所采用数据均是常规燃煤电厂的基本检测数据,不需要价格昂贵的计量设备或检测手段来检测二氧化碳排放量,降低了检测成本,实用性较强。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (5)
1.一种预测燃煤电厂碳排放量的计算方法,其特征在于,所述计算方法包括下列步骤:
S1、建立数学模型,得出燃料特性系数β的关系式:
根据现有的煤质分析数据,利用多元线性回归模型拟合出燃煤的工业分析数据Aad、Vad、FCad、Qgr,ad与燃料特性系数β之间的关系式,其中Aad、Vad、FCad、Qgr,ad分别表示为收到基下燃煤中的灰分含量、挥发分含量、固定碳含量和燃煤的高位发热量,得到褐煤和无烟煤关于燃料特性系数和工业分析数据之间的线性回归方程,Aad、Vad、FCad、Qgr,ad单位分别为%、%、%和MJ/kg;
S2、计算电厂燃烧的煤的燃料特性系数β值:
将电厂测得的燃煤工业分析数据Aad、Vad、FCad、Qgr,ad根据煤种的不同带入所述的褐煤或无烟煤关于燃料特性系数和工业分析数据之间的线性回归方程中,计算出电厂燃煤的燃料特性系数β值;
S3、计算烟气中二氧化碳的百分比含量CO2:
将电厂检测的尾气中氧气和二氧化硫的百分比含量和所述的计算出电厂燃煤的燃料特性系数β值带入二氧化碳含量计算式计算出烟气中二氧化碳的百分比含量;
S4、计算烟气的排放总量Vy;
S5、计算二氧化碳的排放总量
将所述的烟气中二氧化碳的百分比含量CO2和所述的烟气的排放总量Vy带入公式即可计算得到电厂中CO2的实际排放量。
2.根据权利要求1所述的一种预测燃煤电厂碳排放量的计算方法,其特征在于,所述的步骤S1、建立数学模型,得出燃料特性系数β的关系式具体如下:
S11、将无烟煤的工业分析数据和煤的燃料特性系数作为样本数据,并将其带入仿真软件的线性回归分析工具中,则X和Y分别表示为:
式中,Ai、Vi、FCi、Qi分别表示为第i个样本中的空气干燥基灰分、挥发分、固定碳和高位发热量,βi为第i个样本的燃料特性系数,Ai、Vi、FCi、Qi单位分别为%、%、%和MJ/kg;
S12、根据仿真软件分别对褐煤和无烟煤样本数据进行逐步回归分析,得到褐煤和无烟煤关于燃料特性系数和工业分析数据之间的线性回归方程分别如下所示,
3.根据权利要求2所述的一种预测燃煤电厂碳排放量的计算方法,其特征在于,
所述仿真软件为MATLAB,所述线性回归分析工具为stepwise命令,其中,stepwise命令的格式为:stepwise(X,Y),其中X表示自变量数据,为n×m阶矩阵,m表示自变量种类数,n表示每种自变量的样本数,Y表示因变量数据,为n×1阶矩阵,n表示因变量样本数。
4.根据权利要求1所述的一种预测燃煤电厂碳排放量的计算方法,其特征在于,
所述的二氧化碳含量计算式如下:
式中,CO2表示烟气中二氧化碳的体积分数,O2表示烟气中氧气的体积分数,SO2表示烟气中二氧化硫的体积分数。
5.根据权利要求1所述的一种预测燃煤电厂碳排放量的计算方法,其特征在于,
所述的步骤S4、计算烟气的排放总量Vy具体如下:
若电厂已经测得烟气排放的总量Vy,则直接使用该数据;
若电厂未测得烟气排放的总量Vy,则根据燃煤量和过量空气系数进行计算,如下所示:
其中,M表示燃烧所用的燃料量,单位为t,Qnet,ar表示燃料的低位发热量,单位为MJ/kg,α表示空气过量系数,V0表示每燃烧1吨煤所需的理论空气量,其计算式为:
V0=(1.05×Qnet,ar+0.278)×1000
其中,O2表示已检测得到的烟气中氧气的体积分数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611137401.7A CN106650068B (zh) | 2016-12-12 | 2016-12-12 | 一种预测燃煤电厂碳排放量的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611137401.7A CN106650068B (zh) | 2016-12-12 | 2016-12-12 | 一种预测燃煤电厂碳排放量的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106650068A CN106650068A (zh) | 2017-05-10 |
CN106650068B true CN106650068B (zh) | 2019-12-10 |
Family
ID=58824471
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611137401.7A Active CN106650068B (zh) | 2016-12-12 | 2016-12-12 | 一种预测燃煤电厂碳排放量的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106650068B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110681254B (zh) * | 2019-08-26 | 2022-07-08 | 光大环保技术研究院(南京)有限公司 | 一种基于模型的半干法烟气处理控制系统 |
CN113282868A (zh) * | 2020-02-20 | 2021-08-20 | 赫普能源环境科技股份有限公司 | 一种火电厂度电碳排放强度在线监测系统及计算分析方法 |
CN111477286A (zh) * | 2020-03-19 | 2020-07-31 | 北京国电龙源环保工程有限公司 | 一种脱硫废水零排放系统中粉煤灰氯含量的测算方法 |
CN111707782A (zh) * | 2020-03-30 | 2020-09-25 | 江苏方天电力技术有限公司 | 基于氧量的火电机组二氧化碳排放浓度异常检测方法 |
CN112966216A (zh) * | 2021-02-08 | 2021-06-15 | 龙源(北京)碳资产管理技术有限公司 | 一种基于碳氢元素的燃煤热值的计算方法及其用途 |
CN113281466A (zh) * | 2021-05-17 | 2021-08-20 | 吴迅海 | 碳排放量检测的校准方法、装置及计算机存储介质 |
CN113219932B (zh) * | 2021-06-02 | 2023-09-05 | 内蒙古自治区计量测试研究院 | 一种火力发电行业碳排放数字化分析系统 |
CN114140294B (zh) * | 2021-11-26 | 2022-07-12 | 碳阻迹(北京)科技有限公司 | 一种企业碳排放信息实时分析展示的服务方法及系统 |
CN114897226A (zh) * | 2022-04-25 | 2022-08-12 | 广东邦普循环科技有限公司 | 一种碳排放量预测方法、装置、设备及存储介质 |
CN115034454A (zh) * | 2022-05-25 | 2022-09-09 | 北京电链科技有限公司 | 一种基于电量数据的碳排放量预测的方法及装置 |
CN117788021B (zh) * | 2024-02-28 | 2024-05-31 | 广东省特种设备检测研究院顺德检测院 | 一种控排企业入炉煤碳排放数据质量快速甄别方法及系统 |
CN117825640B (zh) * | 2024-03-05 | 2024-06-11 | 国网湖北省电力有限公司电力科学研究院 | 实现电力行业源-网-荷-储全方位碳计量监测的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103294928A (zh) * | 2013-06-28 | 2013-09-11 | 华北电力大学(保定) | 一种碳排放组合预测方法 |
CN104484749A (zh) * | 2014-12-04 | 2015-04-01 | 广东电网有限责任公司电力科学研究院 | 用于燃煤电厂碳排放量预测方法和系统 |
CN105242000A (zh) * | 2015-10-29 | 2016-01-13 | 广东电网有限责任公司电力科学研究院 | 一种燃煤电厂碳排放量的精确测量方法 |
-
2016
- 2016-12-12 CN CN201611137401.7A patent/CN106650068B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103294928A (zh) * | 2013-06-28 | 2013-09-11 | 华北电力大学(保定) | 一种碳排放组合预测方法 |
CN104484749A (zh) * | 2014-12-04 | 2015-04-01 | 广东电网有限责任公司电力科学研究院 | 用于燃煤电厂碳排放量预测方法和系统 |
CN105242000A (zh) * | 2015-10-29 | 2016-01-13 | 广东电网有限责任公司电力科学研究院 | 一种燃煤电厂碳排放量的精确测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106650068A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106650068B (zh) | 一种预测燃煤电厂碳排放量的计算方法 | |
Cao et al. | Toward a better practice for estimating the CO2 emission factors of cement production: An experience from China | |
CN107844682B (zh) | 基于煤气热值和烟气成分的转炉煤气成分软测量方法 | |
CN104008297A (zh) | 煤粉与高炉煤气混烧锅炉的热效率计算方法 | |
CN102734782B (zh) | 一种燃煤锅炉能效监测方法 | |
CN110991756B (zh) | 基于ts模糊神经网络的mswi炉膛温度预测方法 | |
CN102880795A (zh) | 一种锅炉优化运行经济参数目标值获取方法 | |
CN103528846A (zh) | 煤气锅炉热效率与煤气特性的实时监测方法 | |
CN107103421A (zh) | 工业锅炉节能方案的综合评价体系及其构建方法 | |
CN103995987A (zh) | 掺烧高炉煤气的煤粉锅炉的热效率测算方法 | |
Libao et al. | Prediction of CO2 emissions based on multiple linear regression analysis | |
CN105242000A (zh) | 一种燃煤电厂碳排放量的精确测量方法 | |
CN117132123A (zh) | 一种基于原料侧数据的燃煤二氧化碳排放量确定方法及确定系统 | |
CN104598761A (zh) | 混烧锅炉运行参数变化对机组发电煤耗影响的分析方法 | |
CN116757500A (zh) | 一种工业互联网平台监测系统及方法 | |
CN107476903A (zh) | 一种增压柴油机egr性能评估及最佳egr率确定方法 | |
CN114166998A (zh) | 一种水泥生产企业碳排放计量方法和系统 | |
CN104008307A (zh) | 煤粉与高炉煤气混烧锅炉的入炉燃煤量的测算方法 | |
CN114118726A (zh) | 一种燃煤电厂基于煤炭全价值链管理的科学配煤掺烧系统 | |
CN114427688A (zh) | 一种基于气温对锅炉运行影响构造的闭环燃烧优化控制方法 | |
CN109655489B (zh) | 一种燃煤机组标准煤耗的计量方法及装置 | |
CN1322955A (zh) | 基于烟气成分分析的锅炉效率在线监测方法 | |
CN104200285A (zh) | 一种电厂混合燃煤配煤优化方法 | |
CN109635464A (zh) | 一种钢厂煤气含硫量软测量方法 | |
CN115293453A (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 |