CN108280047A - 一种基于现场监测数据的火电机组碳排放核算方法 - Google Patents
一种基于现场监测数据的火电机组碳排放核算方法 Download PDFInfo
- Publication number
- CN108280047A CN108280047A CN201810068380.0A CN201810068380A CN108280047A CN 108280047 A CN108280047 A CN 108280047A CN 201810068380 A CN201810068380 A CN 201810068380A CN 108280047 A CN108280047 A CN 108280047A
- Authority
- CN
- China
- Prior art keywords
- unit
- hourly
- emission
- data
- hour
- 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
- 238000012544 monitoring process Methods 0.000 title claims abstract description 93
- 238000000034 method Methods 0.000 title claims abstract description 61
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 47
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 47
- 238000004364 calculation method Methods 0.000 claims abstract description 72
- 238000010248 power generation Methods 0.000 claims abstract description 37
- 238000000738 capillary electrophoresis-mass spectrometry Methods 0.000 claims abstract description 16
- 230000005611 electricity Effects 0.000 claims abstract description 6
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 134
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 67
- 239000007789 gas Substances 0.000 claims description 38
- 239000001569 carbon dioxide Substances 0.000 claims description 34
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 29
- 239000003546 flue gas Substances 0.000 claims description 29
- 230000008569 process Effects 0.000 claims description 24
- 238000007781 pre-processing Methods 0.000 claims description 19
- 238000006477 desulfuration reaction Methods 0.000 claims description 12
- 230000023556 desulfurization Effects 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 230000008030 elimination Effects 0.000 claims description 6
- 238000003379 elimination reaction Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 description 12
- 239000000463 material Substances 0.000 description 5
- 238000000691 measurement method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 239000003245 coal Substances 0.000 description 4
- 239000005431 greenhouse gas Substances 0.000 description 4
- 235000019738 Limestone Nutrition 0.000 description 3
- 238000013480 data collection Methods 0.000 description 3
- 239000006028 limestone Substances 0.000 description 3
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 2
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 230000003009 desulfurizing effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- FYYHWMGAXLPEAU-UHFFFAOYSA-N Magnesium Chemical compound [Mg] FYYHWMGAXLPEAU-UHFFFAOYSA-N 0.000 description 1
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 1
- 239000003513 alkali Substances 0.000 description 1
- 229910021529 ammonia Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010170 biological method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 229910000019 calcium carbonate Inorganic materials 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000011777 magnesium Substances 0.000 description 1
- 229910052749 magnesium Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000011593 sulfur Substances 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Combustion & Propulsion (AREA)
- Probability & Statistics with Applications (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于现场监测数据的火电机组碳排放核算方法，包括以下步骤：步骤一，监测点选取：将CO2测量仪安装在火电机组烟囱排口处；步骤二，监测数据采集：设置在烟囱排口处的CEMS和CO2测量仪监测到的监测数据通过数据采集模块采集后传输给主站，由主站进行存储；步骤三，小时CO2排放量计算：以1小时作为计算周期，以1小时内采集到的监测数据为基础，剔除其中的停机数据和超限数据，然后计算小时CO2排放量；步骤四，CO2排放量统计，以机组小时CO2排放量为基础，通过求和运算得到CO2排放量的汇总值；步骤五，CO2排放强度计算，机组CO2排放强度是机组每发一度电排放的CO2量，即该时段内CO2排放量与该时段内发电量的比值。
Description
技术领域
本发明涉及一种基于现场监测数据的火电机组碳排放核算方法。
背景技术
在所有温室气体中，二氧化碳是引起全球气候变化的主要气体之一。据统计，二氧化碳对全球升温的影响约为55％。电力行业由于存在动力煤燃烧过程，一直是二氧化碳的主要排放源，据电监会发布的《2010年及“十一五”电力行业节能减排情况通报》显示，2010年全国电力行业二氧化碳排放量约占全国的50％。
江苏作为全国经济大省，也是温室气体排放大省，纳入碳排放权交易市场的重点排放单位达400家之多，位居全国前列。为了更好地贯彻落实国家“全面启动全国碳排放交易市场”的政策部署，省政府办公厅印发了《江苏省碳排放权交易市场建设实施方案》，要求从2017年到2020年建立健全与全国碳排放权交易市场机制相配套的碳排放权交易市场体系。碳排放交易市场建立的前提是碳排放企业及时准确地上报温室气体排放数据。
目前，企业二氧化碳排放量计算方法有三种：
1)物料平衡法：根据质量守恒定律，用输入物料中的含碳量减去输出物料中的含碳量进行平衡计算得到二氧化碳排放量。由于入炉煤固定碳含量、入炉煤灰分含量及灰分中含碳量、脱硫剂石灰石消耗量及石灰石中碳酸钙含量、入炉煤含硫率等参数无法在线测量，只能依靠人工定时统计数据进行计算，因此物料平衡法误差相对较大。而且，物料平衡法仅适用于以石灰石为脱硫剂的脱硫系统，诸如双碱法、氨水法、生物法、镁法等使用其他脱硫剂的系统则不能使用这种算法。
2)排放因子法：主要通过排放因子计算CO2排放量，核心算法可以用公式“活动量×排放因子”来表示。排放因子法由于其计算过程相对简单，目前世界上官方碳排放量计算、碳交易核算等基本都采用这种方法。但是，排放因子的数据主要来源于IPCC(Intergovernmental Panel on Climate Change，政府间气候变化专业委员会)、美国石油学会等专业机构提供的技术资料，是在国外相关燃料、设备、技术的基础上研究出来的。然而在不同技术水平、生产状况、能源使用情况、工艺过程等因素的影响下排放因子存在很大差异。因此，使用排放因子法存在的不确定性比较大，是否适用于我国燃料类型复杂的火电机组也有待商榷。
3)实测法：这是本专利采用的方法。这种方法通过安装监测仪表、设备，如烟气排放连续监测系统(Continuous Emissions Monitoring System，CEMS)，并按照相关技术规范的要求测量和记录排放烟气中CO2的浓度、烟气量等主要参数数据，从而计算单个电力企业的CO2排放量。实测法计算结果准确性较高，而且江苏省内135MW及以上燃煤机组都在烟囱上装有CEMS系统，测量参数包括烟气流速、压力、温度、湿度等，这些测点可用于计算机组二氧化碳排放量。因此，现场也具备了采用实测法核算企业二氧化碳排放量的测量条件。而且，通过连续监测系统实时采集碳排放数据，可以实现火电机组二氧化碳排放量的在线量化计算，有助于提高企业碳排放数据的时效性。
发明内容
本发明为了解决现有技术中存在的上述缺陷和不足，提供了一种基于现场监测数据的火电机组碳排放核算方法，通过连续排放监测的方式进行自动化的温室气体排放数据收集；然后采用实测法，利用现场监测数据对火电机组二氧化碳排放量和排放强度进行量化计算和汇总统计。
为解决上述技术问题，本发明提供一种基于现场监测数据的火电机组碳排放核算方法，其特征在于包括以下步骤：
步骤一，监测点选取：将CO2测量仪安装在火电机组烟囱排口处；
步骤二，监测参数数据采集：设置在烟囱排口处的CEMS和CO2测量仪监测到的监测数据通过数据采集模块采集后传输给主站，由主站进行存储；
步骤三，小时CO2排放量计算：以1小时作为计算周期，1小时内采集到的监测参数数据称为机组小时监测数据集，计算过程如下：
3.1机组停机判断：将机组小时监测数据集中满足停机条件的数据进行剔除，得到机组小时运行数据集；
3.2数据预处理：对机组小时运行数据集进行统一数据单位、剔除异常数据处理，得到机组小时预处理数据集；
3.3机组小时均值计算：由机组小时预处理数据集计算监测点监测参数的算术平均值，计算出的算术平均值即为监测参数的机组小时均值，计算公式如下：
上式中：为监测参数的小时均值；xi是机组小时预处理数据集中第i条记录中监测参数的取值；n是机组小时预处理数据集中的记录个数；
3.4小时CO2排放量的计算，具体过程如下：
3.4.1机组小时运行时间的确定：机组小时运行数据集的数据更新频率为1分钟，即每条记录对应1分钟，因此，机组小时运行时间即为机组小时运行数据集中数据记录个数，单位为分钟，用RunMins表示；
3.4.2机组小时烟气量计算：通过监测参数的机组小时均值和小时运行时间计算机组小时标干烟气量；
3.4.3机组小时CO2排放量计算：通过机组CO2排放浓度小时均值和机组小时标干烟气量计算机组小时CO2排放量；
步骤四，CO2排放量统计：以机组小时CO2排放量为基础，通过求和运算得到CO2排放量的汇总值；
步骤五，CO2排放强度计算，具体过程如下：
5.1机组小时发电量计算：通过对机组负荷进行积分运算得到小时发电量；
5.2某时段内CO2排放强度计算：机组CO2排放强度是机组每发一度电排放的CO2量，即该时段内CO2排放量与该时段内发电量的比值。
其中，步骤二中，数据采集模块的采集方式有两种：
其一，工控机方式：CO2测量仪和CEMS监测到的监测参数数据经由数采仪进行数据采集后，通过部署在工控机上的数据终端发送给主站；
其二，DCS方式：通过CEMS小室与脱硫DCS之间的备用电缆，将CO2测量仪接入脱硫DCS中，经由脱硫DCS将数据写入电厂SIS系统，最后经由电厂SIS系统传输给主站；同时，CEMS监测到的监测参数数据经由数采仪进行数据采集后，通过部署在工控机上的数据终端发送给主站。
步骤3.1中，所述停机条件为：
燃煤机组：FHsteam≤0
燃气机组：FHgas≤0
上式中，FHsteam为汽机负荷，FHgas为燃机负荷。
步骤3.2中，数据预处理的过程如下：
3.2.1标定故障剔除：当CO2测量仪表处于标定(BDvent＝1)或故障(GZvent＝1)状态时，监测参数数据无效，需剔除此段时间的数据记录；
3.2.2单位转换：将监测参数数据的单位均转换为统一单位制后再进行后续计算；
3.2.3超限数据剔除：若监测参数的监测数据超出其正常取值范围，则剔除该条数据记录。
步骤3.4.2中，通过机组小时均值和小时运行时间计算机组小时标干烟气量，计算公式如下：
上式中：Qgas,vent,N是标干状态下机组小时烟气排放量；RunMins是机组小时运行时间；A是烟囱排口的截面积；Vgas.vent、Pgas.vent、Tgas.vent、Xgas.vent分别是烟囱排口处烟气的流速、压力、温度和湿度。
步骤3.4.3中，通过机组CO2排放浓度小时均值和机组小时标干烟气量计算机组小时CO2排放量，计算过程如下：
上式中，是机组小时二氧化碳排放量；Qgas,vent,N是标干状态下机组小时烟气排放量；是二氧化碳排放浓度小时均值，是将二氧化碳排放浓度从百分比浓度转换成百万分比浓度ppm的转换系数；44/22.4是二氧化碳体积浓度和质量浓度的转换系数；10-9是毫克mg到吨t的转换系数。
步骤四中，以机组小时CO2排放量为基础，通过求和运算得到某时段内CO2排放量的汇总统计值，计算公式如下：
上式中，ECO2是该时段的机组CO2排放量；是该时段中某个小时的机组CO2排放量；k是该时段内的小时数，i是小时下标。
步骤5.1中，小时发电量是负荷在时间上的积分值，计算公式如下：
上式中：FDLh是机组小时发电量；FH是机组负荷时刻值；Δt是时间间隔，除以60是将其从分钟转换为小时；t1、t2分别是积分的开始和结束时刻；
对于燃煤机组，利用上述公式即可计算出机组的小时发电量；
对于燃气机组，利用上述公式分别计算出燃机和汽机的小时发电量，然后将两者相加得到机组的小时发电量，计算公式如下：
式中：FDLh、分别是机组、燃机和汽机的小时发电量。
步骤5.2中，某时段内CO2排放强度为机组每发一度电排放的CO2量，计算公式如下：
式中：是该时段的机组CO2排放强度；是该时段的机组CO2排放量，106是将CO2排放量从吨t转换成克g的转换系数；FDL是该时段的机组发电量，计算公式为其中是该时段中某个小时的机组发电量，k是该时段内的小时数，i是小时下标，103是将发电量从兆瓦时MWh转换成千瓦时kWh的转换系数。
本发明所达到的有益技术效果：本发明提供的一种基于现场监测数据的火电机组碳排放核算方法，通过连续监测系统实时采集碳排放数据，然后采用实测法在线计算火电机组二氧化碳排放量和排放强度。这种方法可以有效缩短数据流路线、减少人工干预，提高了企业碳排放数据的准确性。且该方法逻辑简单，易于实现在线量化计算，提高了企业碳排放数据的时效性，为政府部门及时掌握火电企业二氧化碳排放情况、开展火电企业碳排放交易、制定碳减排目标、分配碳排放配额提供坚实可靠的数据支撑。
附图说明
图1本发明计算流程示意图；
图2本发明实施例中测试机组于2017年10月31日二氧化碳小时排放强度趋势图。
具体实施方式
下面结合具体实施例对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
下面结合附图和实施例对本发明专利进一步说明。
如图1所示，本发明提供一种基于现场监测数据的火电机组碳排放核算方法，包括以下步骤：
步骤一，监测点选取：由于火电机组在脱硫过程中也会产生二氧化碳排放，而且二氧化碳在脱硫系统中的变化量无法通过计算实时得出。如果将CO2测量仪表安装在锅炉出口，无法准确实时计算出机组总的二氧化碳排放量。同时火电机组烟囱排口处的CEMS都设有流速、压力、温度、湿度测点，也便于计算烟气排放量，因此将CO2测量仪表安装在烟囱排口位置；
步骤二，监测参数数据采集：设置在烟囱排口处的CEMS和CO2测量仪监测到的监测数据通过数据采集模块采集后传输给主站，由主站进行存储；监测参数数据采集的方式有两种：
其一，工控机方式：CO2测量仪和CEMS监测到的监测参数数据经由数采仪进行数据采集后，通过部署在工控机上的数据终端发送给主站；
其二，DCS方式：通过CEMS小室与脱硫DCS之间的备用电缆，将CO2测量仪接入脱硫DCS中，经由脱硫DCS将数据写入电厂SIS系统，最后经由电厂SIS系统传输给主站；同时，CEMS监测到的监测参数数据经由数采仪进行数据采集后，通过部署在工控机上的数据终端发送给主站。
其中监测参数如表1所示：
表1机组的监测参数
步骤三，小时CO2排放量计算：机组正常运行过程中，只要设备没有出现故障，在1小时内系统运行状态变化不大；其次，考虑到现场运行参数数据自身的波动性和数据采集过程中尖峰脉冲数据处理等因素，对碳排放进行实时计算容易放大数据波动的影响而忽视数据固有的内在特征。同时，在国家关于固定污染源烟气排放连续监测的技术规范和要求中，对监测数据的统计、审核、失控处理、缺失处理等工作都是以小时为周期。因此，本方法以1小时作为计算周期。目前，现场碳排放监测参数的数据采集频率最快可以到秒级，最慢也在1分钟左右。因此，1小时计算周期内会有60～3600条记录，形成了机组小时监测数据集。机组小时CO2排放量，以1小时作为计算周期，1小时内采集到的数据称为机组小时监测数据集，计算过程如下：
3.1机组停机判断：机组停机判断，一是为了保证机组处于运行状态，二是因为测量仪表的读数在机组停机时并不总是会归零，因此需要通过机组停机判断来剔除这部分数据。将机组小时监测数据集中满足停机条件的数据进行剔除，得到机组小时运行数据集；所述停机条件为：
燃煤机组：FHsteam≤0
燃气机组：FHgas≤0
上式中，FHsteam为汽机负荷，FHgas为燃机负荷。
3.2数据预处理：对机组小时运行数据集进行统一数据单位、剔除异常数据处理，得到机组小时预处理数据集；数据预处理的过程如下：
3.2.1标定故障剔除：CO2测量仪表处于标定(BDvent＝1)或故障(GZvent＝1)状态时，监测参数数据无效，需剔除此段时间的数据记录；
3.2.2单位转换：将监测参数的单位均转换为统一单位制再进行后续计算；
3.2.3超限数据剔除：若监测参数的监测数据超出其正常取值范围，则剔除该条数据记录，机组的监测参数数据正常取值范围如表2所示：
表2机组的监测参数数据正常取值范围
3.3机组小时均值计算：由机组小时预处理数据集计算监测点监测参数的算术平均值，计算出的算术平均值即为机组监测参数的小时均值，计算公式如下：
上式中：为监测参数的小时均值；xi是机组小时预处理数据集中第i条记录中监测参数的取值；n是机组小时预处理数据集中的记录个数；
3.4小时CO2排放量的计算，具体过程如下：
3.4.1机组小时运行时间的确定：机组小时运行数据集的数据更新频率为1分钟，即每条记录对应1分钟，因此，机组小时运行时间即为机组小时运行数据集中数据记录个数，单位为分钟，用RunMins表示；
3.4.2机组小时烟气量计算：通过机组小时均值和小时运行时间计算机组小时标干烟气量，计算公式如下：
上式中：Qgas，vent，N是标干状态下机组小时烟气排放量；RunMins是机组小时运行时间；A是烟囱排口的截面积；Vgas.vent、Pgas.vent、Tgas.vent、Xgas.vent分别是烟囱排口处烟气的流速、压力、温度和湿度。3.4.3机组小时CO2排放量计算：通过机组CO2排放浓度小时均值和机组小时标干烟气量计算机组小时CO2排放量，计算过程如下：
上式中，是机组小时二氧化碳排放量；Qgas,vent,N是标干状态下机组小时烟气排放量；是二氧化碳排放浓度小时均值，是将二氧化碳排放浓度从百分比浓度转换成百万分比浓度ppm的转换系数；44/22.4是二氧化碳体积浓度和质量浓度的转换系数；10-9是毫克mg到吨t的转换系数。
步骤四，CO2排放量统计：以机组小时CO2排放量为基础，通过求和运算得到某时段内CO2排放量的汇总统计值，计算公式如下：
上式中，是该时段的机组CO2排放量；是该时段中某个小时的机组CO2排放量；k是该时段内的小时数，i是小时下标。
步骤五，CO2排放强度计算，具体过程如下：
5.1机组小时发电量计算：小时发电量是负荷在时间上的积分值，计算公式如下：
上式中：FDLh是机组小时发电量；FH是机组负荷时刻值；Δt是时间间隔，除以60是将其从分钟转换为小时；t1、t2分别是积分的开始和结束时刻；
对于燃煤机组，利用上述公式即可计算出机组的小时发电量；
对于燃气机组，利用上述公式分别计算出燃机和汽机的小时发电量，然后将两者相加得到机组的小时发电量，计算公式如下：
式中：FDLh、分别是机组、燃机和汽机的小时发电量。
5.2.某时段内CO2排放强度计算：某时段内CO2排放强度为机组每发一度电排放的CO2量，计算公式如下：
式中：是该时段的机组CO2排放强度；是该时段的机组CO2排放量，106是将CO2排放量从吨t转换成克g的转换系数；FDL是该时段的机组发电量，计算公式为其中是该时段中某个小时的机组发电量，k是该时段内的小时数，i是小时下标，103是将发电量从兆瓦时MWh转换成千瓦时kWh的转换系数。
实施例
下面以具体实施例对本发明做进一步的阐述：
以一台9F燃机作为测试机组为例，选取该机组2017年10月31日全天1分钟间隔的发电负荷和烟囱排口烟气参数监测数据进行说明。测试机组燃机负荷、汽机负荷分别为282.3、156.51MW，烟囱排口处截面积为40.69m2，采用工控机方式实现数据传输。该机组在2017年10月31日6点31分启动，7点30分启动过程基本结束，之后一直处于运行状态。
1、小时碳排放量计算
(1)正常运行过程
表3是2017年10月31日8点到9点这1小时的监测数据列表，时间间隔为1分钟。
表3测试机组2017年10月31日8点到9点1分钟间隔的监测数据
根据机组停机判断条件，测试机组在2017年10月31日8点这一个小时内都处于运行状态，因此小时运行时间为60分钟。经过数据预处理剔除标定故障数据和超限数据记录后，得到排口各烟气参数的小时均值和小时二氧化碳排放数据如表4所示。
表4测试机组2017年10月31日8点小时碳排放数据
(2)机组启动过程
表5是2017年10月31日7点至8点这1小时的监测数据列表。从表5中可以看出，在这1小时中虽然汽机有一段时间处于停机状态，但是燃机一直在运行，因此这小时的运行时间仍然为60分钟。
表5测试机组2017年10月31日7点到8点1分钟间隔的监测数据
测试机组2017年10月31日7点的小时碳排放计算结果如表6所示。
表6测试机组2017年10月31日7点小时碳排放数据
比较表4和表6的计算结果，可以看出测试机组启动过程中二氧化碳排放强度明显大于正常运行过程，这主要是因为启动过程中发电量较小，而二氧化碳排放量却没有成比例下降，使得两者的比值有较大的增长。
2、全天碳排放量统计
表7列出了2017年10月31日从机组启动开始每个小时的碳排放数据，图2展示了该日二氧化碳小时排放强度变化趋势。从中可以看出，6～7点这2个小时由于机组启动，负荷处于爬升过程没有达到稳定状态，因此发电量偏小使得二氧化碳排放强度明显偏大。在剩余时间里，机组一直处于运行状态，排放强度稳定在319.43g/kWh左右。而全天的二氧化碳排放强度为325.89g/kWh，比运行状态的数值略有上升，反映了测试机组全天二氧化碳排放平均水平。
表7测试机组2017年10月31日小时碳排放数据
以上已以较佳实施例公布了本发明，然其并非用以限制本发明，凡采取等同替换或等效变换的方案所获得的技术方案，均落在本发明的保护范围内。
Claims (9)
1.一种基于现场监测数据的火电机组碳排放核算方法，其特征在于包括以下步骤：
步骤一，监测点选取：将CO2测量仪安装在火电机组烟囱排口处；
步骤二，监测参数数据采集：设置在烟囱排口处的CEMS和CO2测量仪监测到的监测数据通过数据采集模块采集后传输给主站，由主站进行存储；
步骤三，小时CO2排放量计算：以1小时作为计算周期，1小时内采集到的监测参数数据称为机组小时监测数据集，计算过程如下：
3.1机组停机判断：将机组小时监测数据集中满足停机条件的数据进行剔除，得到机组小时运行数据集；
3.2数据预处理：对机组小时运行数据集进行统一数据单位、剔除异常数据处理，得到机组小时预处理数据集；
3.3机组小时均值计算：由机组小时预处理数据集计算监测点监测参数的算术平均值，计算出的算术平均值即为监测参数的机组小时均值，计算公式如下：
上式中：为监测参数的小时均值；xi是机组小时预处理数据集中第i条记录中监测参数的取值；n是机组小时预处理数据集中的记录个数；
3.4小时CO2排放量的计算，具体过程如下：
3.4.1机组小时运行时间的确定：机组小时运行数据集的数据更新频率为1分钟，即每条记录对应1分钟，因此，机组小时运行时间即为机组小时运行数据集中数据记录个数，单位为分钟，用RunMins表示；
3.4.2机组小时烟气量计算：通过监测参数的机组小时均值和小时运行时间计算机组小时标干烟气量；
3.4.3机组小时CO2排放量计算：通过机组CO2排放浓度小时均值和机组小时标干烟气量计算机组小时CO2排放量；
步骤四，CO2排放量统计：以机组小时CO2排放量为基础，通过求和运算得到CO2排放量的汇总值；
步骤五，CO2排放强度计算，具体过程如下：
5.1机组小时发电量计算：通过对机组负荷进行积分运算得到小时发电量；
5.2某时段内CO2排放强度计算：机组CO2排放强度是机组每发一度电排放的CO2量，即该时段内CO2排放量与该时段内发电量的比值。
2.根据权利要求1所述的基于现场监测数据的火电机组碳排放核算方法，其特征在于：步骤二中，数据采集模块的采集方式有两种：
其一，工控机方式：CO2测量仪和CEMS监测到的监测参数数据经由数采仪进行数据采集后，通过部署在工控机上的数据终端发送给主站；
其二，DCS方式：通过CEMS小室与脱硫DCS之间的备用电缆，将CO2测量仪接入脱硫DCS中，经由脱硫DCS将数据写入电厂SIS系统，最后经由电厂SIS系统传输给主站；同时，CEMS监测到的监测参数数据经由数采仪进行数据采集后，通过部署在工控机上的数据终端发送给主站。
3.根据权利要求1所述的基于现场监测数据的火电机组碳排放核算方法，其特征在于：步骤3.1中，所述停机条件为：
燃煤机组：FHsteam≤0
燃气机组：FHgas≤0
上式中，FHsteam为汽机负荷，FHgas为燃机负荷。
4.根据权利要求1所述的基于现场监测数据的火电机组碳排放核算方法，其特征在于：步骤3.2中数据预处理的过程如下：
3.2.1标定故障剔除：当CO2测量仪表处于标定(BDvent＝1)或故障(GZvent＝1)状态时，监测参数数据无效，需剔除此段时间的数据记录；
3.2.2单位转换：将监测参数数据的单位均转换为统一单位制后再进行后续计算；
3.2.3超限数据剔除：若监测参数的监测数据超出其正常取值范围，则剔除该条数据记录。
5.根据权利要求1所述的基于现场监测数据的火电机组碳排放核算方法，其特征在于：步骤3.4.2中，通过机组小时均值和小时运行时间计算机组小时标干烟气量，计算公式如下：
上式中：Qgas,vent,N是标干状态下机组小时烟气排放量；RunMins是机组小时运行时间；A是烟囱排口的截面积；Vgas.vent、Pgas.vent、Tgas.vent、Xgas.vent分别是烟囱排口处烟气的流速、压力、温度和湿度。
6.根据权利要求1所述的基于现场监测数据的火电机组碳排放核算方法，其特征在于：步骤3.4.3中，通过机组CO2排放浓度小时均值和机组小时标干烟气量计算机组小时CO2排放量，计算过程如下：
上式中，是机组小时二氧化碳排放量；Qgas,vent,N是标干状态下机组小时烟气排放量；是二氧化碳排放浓度小时均值，是将二氧化碳排放浓度从百分比浓度转换成百万分比浓度ppm的转换系数；44/22.4是二氧化碳体积浓度和质量浓度的转换系数；10-9是毫克mg到吨t的转换系数。
7.根据权利要求1所述的基于现场监测数据的火电机组碳排放核算方法，其特征在于：步骤四中，以机组小时CO2排放量为基础，通过求和运算得到某时段内CO2排放量的汇总统计值，计算公式如下：
上式中，是该时段的机组CO2排放量；是该时段中某个小时的机组CO2排放量；k是该时段内的小时数，i是小时下标。
8.根据权利要求1所述的基于现场监测数据的火电机组碳排放核算方法，其特征在于：步骤5.1中，小时发电量是负荷在时间上的积分值，计算公式如下：
上式中：FDLh是机组小时发电量；FH是机组负荷时刻值；Δt是时间间隔，除以60是将其从分钟转换为小时；t1、t2分别是积分的开始和结束时刻；
对于燃煤机组，利用上述公式即可计算出机组的小时发电量；
对于燃气机组，利用上述公式分别计算出燃机和汽机的小时发电量，然后将两者相加得到机组的小时发电量，计算公式如下：
式中：FDLh、分别是机组、燃机和汽机的小时发电量。
9.根据权利要求1所述的基于现场监测数据的火电机组碳排放核算方法，其特征在于：步骤5.2中，某时段内CO2排放强度为机组每发一度电排放的CO2量，计算公式如下：
式中：是该时段的机组CO2排放强度；是该时段的机组CO2排放量，106是将CO2排放量从吨t转换成克g的转换系数；FDL是该时段的机组发电量，计算公式为其中是该时段中某个小时的机组发电量，k是该时段内的小时数，i是小时下标，103是将发电量从兆瓦时MWh转换成千瓦时kWh的转换系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810068380.0A CN108280047A (zh) | 2018-01-24 | 2018-01-24 | 一种基于现场监测数据的火电机组碳排放核算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810068380.0A CN108280047A (zh) | 2018-01-24 | 2018-01-24 | 一种基于现场监测数据的火电机组碳排放核算方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108280047A true CN108280047A (zh) | 2018-07-13 |
Family
ID=62804836
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810068380.0A Pending CN108280047A (zh) | 2018-01-24 | 2018-01-24 | 一种基于现场监测数据的火电机组碳排放核算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108280047A (zh) |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109190967A (zh) * | 2018-08-29 | 2019-01-11 | 华电电力科学研究院有限公司 | 一种用于火力发电机组的碳排放量核算方法和系统 |
| CN109443423A (zh) * | 2018-10-10 | 2019-03-08 | 江苏方天电力技术有限公司 | 一种基于积差分析的燃气机组碳排放数据关联度检验方法 |
| CN109557869A (zh) * | 2018-11-27 | 2019-04-02 | 江苏方天电力技术有限公司 | 一种火电机组碳排放在线监测管理系统 |
| CN109740301A (zh) * | 2019-03-14 | 2019-05-10 | 华北电力大学 | 一种基于bp神经网络的燃气机组碳排放量的核算方法 |
| CN110044852A (zh) * | 2019-04-23 | 2019-07-23 | 华南理工大学 | 基于激光诱导击穿光谱法的燃煤电厂碳排放在线监测方法 |
| CN111191881A (zh) * | 2019-12-13 | 2020-05-22 | 大唐东北电力试验研究院有限公司 | 一种基于大数据的火电机组工业设备状态监测方法 |
| CN111241465A (zh) * | 2020-03-06 | 2020-06-05 | 赫普能源环境科技股份有限公司 | 一种用户电表耦合碳排放量监测系统和方法 |
| CN111399483A (zh) * | 2020-03-31 | 2020-07-10 | 大唐贵州发耳发电有限公司 | 一种带有故障剔除功能的和利时dcs模拟量均值计算方法 |
| CN112669061A (zh) * | 2020-12-23 | 2021-04-16 | 华电电力科学研究院有限公司 | 火电厂碳配额盈缺分析方法、系统、装置及可读存储介质 |
| CN113282868A (zh) * | 2020-02-20 | 2021-08-20 | 赫普能源环境科技股份有限公司 | 一种火电厂度电碳排放强度在线监测系统及计算分析方法 |
| CN113344462A (zh) * | 2021-07-09 | 2021-09-03 | 国家电网有限公司华东分部 | 电力现货市场的碳排放水平量化方法、装置及电子设备 |
| CN113640466A (zh) * | 2021-08-03 | 2021-11-12 | 广东电网有限责任公司 | 一种碳排放强度测量方法、设备及介质 |
| CN113778316A (zh) * | 2021-08-18 | 2021-12-10 | 神华神东电力有限责任公司 | 一种模拟量小时均值的计算方法、装置及电子设备 |
| CN113837600A (zh) * | 2021-09-18 | 2021-12-24 | 西安热工研究院有限公司 | 一种燃煤电厂碳排放量指标实时监测系统及方法 |
| CN114062759A (zh) * | 2021-10-28 | 2022-02-18 | 阿凡提物联网科技(沈阳)有限公司 | 一种碳排放监测核对系统及方法 |
| CN114382547A (zh) * | 2022-01-17 | 2022-04-22 | 天地（常州）自动化股份有限公司 | 煤炭开采领域温室气体实时监测计量方法以及监测系统 |
| CN114384203A (zh) * | 2021-12-03 | 2022-04-22 | 国网青海省电力公司 | 一种钢铁企业碳排放量在线监测方法及装置 |
| CN114418791A (zh) * | 2021-11-29 | 2022-04-29 | 华能集团技术创新中心有限公司 | 一种火电厂小时级碳排放在线监测方法 |
| CN114579818A (zh) * | 2022-03-11 | 2022-06-03 | 张巧铃 | 一种可视化碳排放检测管理系统及方法 |
| CN114894981A (zh) * | 2022-05-26 | 2022-08-12 | 广东信拓网络科技有限公司 | 一种高时效的碳排放环境监测计量设备及其计量系统 |
| WO2023038579A3 (en) * | 2021-09-10 | 2023-05-19 | Envision Digital International Pte. Ltd. | Method and apparatus for calculating carbon intensities, terminal and storage medium |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103323273A (zh) * | 2013-05-24 | 2013-09-25 | 国家电网公司 | 一种检测电站锅炉空气预热器性能的方法 |
| CN103399536A (zh) * | 2013-07-15 | 2013-11-20 | 冶金自动化研究设计院 | 一种长流程钢铁企业co2排放量的监测系统及方法 |
| CN103399535A (zh) * | 2013-07-15 | 2013-11-20 | 冶金自动化研究设计院 | 一种短流程钢铁企业co2排放量的监测系统及方法 |
| CN106339585A (zh) * | 2016-08-23 | 2017-01-18 | 华电电力科学研究院 | 一种燃煤发电机组碳排放强度评估方法 |
| CN106651171A (zh) * | 2016-12-19 | 2017-05-10 | 江苏方天电力技术有限公司 | 一种超低排放机组在线监控管理系统 |
| CN206919930U (zh) * | 2017-06-22 | 2018-01-23 | 广西壮族自治区环境保护科学研究院 | 一种碳排放监测系统 |
-
2018
- 2018-01-24 CN CN201810068380.0A patent/CN108280047A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103323273A (zh) * | 2013-05-24 | 2013-09-25 | 国家电网公司 | 一种检测电站锅炉空气预热器性能的方法 |
| CN103399536A (zh) * | 2013-07-15 | 2013-11-20 | 冶金自动化研究设计院 | 一种长流程钢铁企业co2排放量的监测系统及方法 |
| CN103399535A (zh) * | 2013-07-15 | 2013-11-20 | 冶金自动化研究设计院 | 一种短流程钢铁企业co2排放量的监测系统及方法 |
| CN106339585A (zh) * | 2016-08-23 | 2017-01-18 | 华电电力科学研究院 | 一种燃煤发电机组碳排放强度评估方法 |
| CN106651171A (zh) * | 2016-12-19 | 2017-05-10 | 江苏方天电力技术有限公司 | 一种超低排放机组在线监控管理系统 |
| CN206919930U (zh) * | 2017-06-22 | 2018-01-23 | 广西壮族自治区环境保护科学研究院 | 一种碳排放监测系统 |
Non-Patent Citations (1)
| Title |
|---|
| 国家环境监测总站等: "《环境监测方法标准实用手册 第5册 监测技术规范》", 30 April 2013 * |
Cited By (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109190967A (zh) * | 2018-08-29 | 2019-01-11 | 华电电力科学研究院有限公司 | 一种用于火力发电机组的碳排放量核算方法和系统 |
| CN109190967B (zh) * | 2018-08-29 | 2020-09-01 | 华电电力科学研究院有限公司 | 一种用于火力发电机组的碳排放量核算方法和系统 |
| CN109443423A (zh) * | 2018-10-10 | 2019-03-08 | 江苏方天电力技术有限公司 | 一种基于积差分析的燃气机组碳排放数据关联度检验方法 |
| CN109443423B (zh) * | 2018-10-10 | 2021-09-07 | 江苏方天电力技术有限公司 | 一种基于积差分析的燃气机组碳排放数据关联度检验方法 |
| CN109557869A (zh) * | 2018-11-27 | 2019-04-02 | 江苏方天电力技术有限公司 | 一种火电机组碳排放在线监测管理系统 |
| CN109740301A (zh) * | 2019-03-14 | 2019-05-10 | 华北电力大学 | 一种基于bp神经网络的燃气机组碳排放量的核算方法 |
| CN110044852A (zh) * | 2019-04-23 | 2019-07-23 | 华南理工大学 | 基于激光诱导击穿光谱法的燃煤电厂碳排放在线监测方法 |
| CN111191881A (zh) * | 2019-12-13 | 2020-05-22 | 大唐东北电力试验研究院有限公司 | 一种基于大数据的火电机组工业设备状态监测方法 |
| CN111191881B (zh) * | 2019-12-13 | 2024-05-14 | 大唐东北电力试验研究院有限公司 | 一种基于大数据的火电机组工业设备状态监测方法 |
| CN113282868A (zh) * | 2020-02-20 | 2021-08-20 | 赫普能源环境科技股份有限公司 | 一种火电厂度电碳排放强度在线监测系统及计算分析方法 |
| CN111241465A (zh) * | 2020-03-06 | 2020-06-05 | 赫普能源环境科技股份有限公司 | 一种用户电表耦合碳排放量监测系统和方法 |
| CN111399483A (zh) * | 2020-03-31 | 2020-07-10 | 大唐贵州发耳发电有限公司 | 一种带有故障剔除功能的和利时dcs模拟量均值计算方法 |
| CN112669061A (zh) * | 2020-12-23 | 2021-04-16 | 华电电力科学研究院有限公司 | 火电厂碳配额盈缺分析方法、系统、装置及可读存储介质 |
| CN113344462A (zh) * | 2021-07-09 | 2021-09-03 | 国家电网有限公司华东分部 | 电力现货市场的碳排放水平量化方法、装置及电子设备 |
| CN113640466A (zh) * | 2021-08-03 | 2021-11-12 | 广东电网有限责任公司 | 一种碳排放强度测量方法、设备及介质 |
| CN113778316A (zh) * | 2021-08-18 | 2021-12-10 | 神华神东电力有限责任公司 | 一种模拟量小时均值的计算方法、装置及电子设备 |
| WO2023038579A3 (en) * | 2021-09-10 | 2023-05-19 | Envision Digital International Pte. Ltd. | Method and apparatus for calculating carbon intensities, terminal and storage medium |
| CN113837600A (zh) * | 2021-09-18 | 2021-12-24 | 西安热工研究院有限公司 | 一种燃煤电厂碳排放量指标实时监测系统及方法 |
| CN114062759A (zh) * | 2021-10-28 | 2022-02-18 | 阿凡提物联网科技(沈阳)有限公司 | 一种碳排放监测核对系统及方法 |
| CN114418791A (zh) * | 2021-11-29 | 2022-04-29 | 华能集团技术创新中心有限公司 | 一种火电厂小时级碳排放在线监测方法 |
| CN114418791B (zh) * | 2021-11-29 | 2024-05-31 | 华能集团技术创新中心有限公司 | 一种火电厂小时级碳排放在线监测方法 |
| CN114384203B (zh) * | 2021-12-03 | 2023-10-20 | 国网青海省电力公司 | 一种钢铁企业碳排放量在线监测方法及装置 |
| CN114384203A (zh) * | 2021-12-03 | 2022-04-22 | 国网青海省电力公司 | 一种钢铁企业碳排放量在线监测方法及装置 |
| CN114382547A (zh) * | 2022-01-17 | 2022-04-22 | 天地（常州）自动化股份有限公司 | 煤炭开采领域温室气体实时监测计量方法以及监测系统 |
| CN114579818A (zh) * | 2022-03-11 | 2022-06-03 | 张巧铃 | 一种可视化碳排放检测管理系统及方法 |
| CN114894981A (zh) * | 2022-05-26 | 2022-08-12 | 广东信拓网络科技有限公司 | 一种高时效的碳排放环境监测计量设备及其计量系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108280047A (zh) | 一种基于现场监测数据的火电机组碳排放核算方法 | |
| Zhang et al. | Energy-exergy analysis and energy efficiency improvement of coal-fired industrial boilers based on thermal test data | |
| CN106650068B (zh) | 一种预测燃煤电厂碳排放量的计算方法 | |
| CN103697946A (zh) | 一种火电厂燃煤锅炉烟气流量的计算方法及污染物排放量的控制方法 | |
| CN103955202B (zh) | 一种基于燃煤电厂脱硫系统数据自动诊断甄别方法 | |
| CN101865867A (zh) | 一种实时计算煤质元素和工业成分的方法 | |
| CN100526868C (zh) | 大型锅炉排烟含氧量软测量方法 | |
| CN111522323A (zh) | 一种基于物联网技术的锅炉能效在线诊断及智能控制方法 | |
| CN116664161B (zh) | 一种基于燃煤火电厂的二氧化碳排放核算技术选择方法 | |
| LU503892B1 (en) | Online carbon emission monitoring system and method for thermal power units | |
| CN113744079B (zh) | 一种综合能源系统的低碳优化运行方法 | |
| CN114663266A (zh) | 燃煤电厂碳排放量在线监测装置及其核算方法、存储介质 | |
| CN117557111A (zh) | 服务于实时碳交易的燃煤电厂碳排放因子计算方法及系统 | |
| CN111339479B (zh) | 一种火电机组污染物达标排放实时浓度限值的计算方法 | |
| CN103728071B (zh) | 一种火电机组最大出力测量方法 | |
| CN117350458A (zh) | 一种基于在线监测的火电机组碳排放分析方法 | |
| CN115187012A (zh) | 一种基于日数据验证燃煤电厂燃烧产生的碳排放方法 | |
| CN114118726A (zh) | 一种燃煤电厂基于煤炭全价值链管理的科学配煤掺烧系统 | |
| CN116068128A (zh) | 一种生物质气化耦合机组co2在线监测计量系统及方法 | |
| CN106527382A (zh) | 一种焦炉热效率在线监测方法 | |
| CN117314706A (zh) | 一种燃煤机组二氧化碳排放量计算方法 | |
| Liu et al. | Comparative Analysis of CO₂ Emission Accounting Methods for Coal-fired Power Plants | |
| CN106557078A (zh) | 一种1000mw机组配煤优化系统 | |
| CN114627976A (zh) | 一种电站水系统运行过程碳排放计算方法 | |
| CN109489766B (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: 20180713 |
|
| RJ01 | Rejection of invention patent application after publication |