CN105956370B - 一种用于燃煤电厂实时计算煤耗等经济指标的方法 - Google Patents

一种用于燃煤电厂实时计算煤耗等经济指标的方法 Download PDF

Info

Publication number
CN105956370B
CN105956370B CN201610254860.7A CN201610254860A CN105956370B CN 105956370 B CN105956370 B CN 105956370B CN 201610254860 A CN201610254860 A CN 201610254860A CN 105956370 B CN105956370 B CN 105956370B
Authority
CN
China
Prior art keywords
coal
power plant
burning power
consumption
real time
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
Application number
CN201610254860.7A
Other languages
English (en)
Other versions
CN105956370A (zh
Inventor
何翔
马达夫
孟永杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Power Equipment Research Institute Co Ltd
Original Assignee
Shanghai Power Equipment Research Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Power Equipment Research Institute Co Ltd filed Critical Shanghai Power Equipment Research Institute Co Ltd
Priority to CN201610254860.7A priority Critical patent/CN105956370B/zh
Publication of CN105956370A publication Critical patent/CN105956370A/zh
Application granted granted Critical
Publication of CN105956370B publication Critical patent/CN105956370B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Control Of Steam Boilers And Waste-Gas Boilers (AREA)
  • Regulation And Control Of Combustion (AREA)

Abstract

本发明提供了一种用于燃煤电厂实时计算煤耗等经济指标的方法,对采集的数据形成标准数据库;再对燃煤电厂汽轮机热耗率、管道效率和锅炉效率进行同时计算,流量参数只用到除氧器入口给水流量,用热平衡计算的方式避开众多汽水参数的流量测量值,尤其在锅炉效率的计算过程中,建立了核心方程组,仅飞灰、大渣含碳量是手动输入数据,其余全部是采集的运行数据和热平衡计算得出的中间参数,从而得出锅炉效率的实时计算值,结合汽轮机热耗率和管道效率,最终得出燃煤电厂实时煤耗等指标。本发明方法使得燃煤电厂经济性指标可以实时计算得到,为后续的耗差分析、不确定度分析等也提供了实时计算的可能,使得燃煤电厂深度节能降耗得到了极大的便利。

Description

一种用于燃煤电厂实时计算煤耗等经济指标的方法
技术领域
本发明涉及一种燃煤电厂煤耗等经济指标的计算方法,用于日常机组运行实时的经济性指标计算与统计,采用该方法后,更将为燃煤电厂设备运行性能监测、节能减排等工作的开展带来极大的便利性和实时性,本发明提供的方法属于火力发电厂燃煤机组经济指标领域。
背景技术
我国目前火电占据所有发电规模的70%左右,根据电力工业中长期发展新的预测,到2020年,我国装机总容量将达到13.4亿千瓦,其中煤电装机总容量为9.1亿千瓦。以2005年3.8亿千瓦煤电装机量为基准,到2020年煤电装机新增容量为5.3亿千瓦。预计到2050年,我国火电机组发电方式仍高居所有发电量的55%~60%。
燃煤电厂实时经济性指标与性能监测在当前遇到两个方面的难题:
第一方面(真实性):
燃煤电厂对日常煤耗、设备效率等一些重要的经济性指标计算方法掌握不够全面,对于煤耗涉及到的锅炉效率主要采用正平衡或者简化后的反平衡进行计算。由于正平衡计算需要精确测量较多重要参数的流量,如主蒸汽流量、给煤流量、给水流量等,而关于流量测量的误差对最终结果的影响较大,使得煤耗计算等经济性指标实际上并不准确,就导致了数据参考性不强的结果。为了知道设备性能,燃煤电厂一般会在机组大、小修后,委托电力科学研究院或其他具有相关资质的单位进行机组设备的效率、性能考核试验,从而得出试验条件下的机组煤耗、锅炉效率、汽轮机热耗率等重要指标数据。这些试验均根据相关标准进行,如:美国机械工程师ASME标准、我国电站锅炉、汽轮机性能试验的国标等。而且根据不同的标准,在性能试验过程要求系统一些参数隔离,如:锅炉试验过程中不吹灰、汽机试验过程中除氧器不补水等。但是在机组日常运行过程中,这些参数的隔离是不可能的,也即,机组大、小修后的设备性能试验没有给出机组的实际运行经济性,得到的机组指标性数据的参考性也不强,本质就是试验条件下的机组性能与实际日常运行下的是不同的。
第二方面(实时性):
火电厂希望知道机组真实、实时的能耗和排放水平,煤炭供应市场变化频繁,电厂每个月甚至每周燃用的煤来源并不相同,而且入炉煤煤样的元素分析化验是当日取样、送检,结果出来后已经是次日后了,而次日后有可能入炉煤煤样已经改变。目前大多是电厂日常能耗统计是通过简易的反平衡计算(入炉煤元素分析大多是假设、修正得出),再结合正平衡对比得出一个报表结果,由于性能试验不可能每天进行,故更为准确的、快速的日常运行中的经济性指标的获得就显得日益重要。
对于上述困难,已经有一些学者、工程师开展了研究。他们利用燃煤电厂自身SIS系统数据进行经济性指标计算的案例已有很多,也有成熟的商业软件,但是这些案例中或多或少存在一些缺陷,比如:湿蒸汽湿度的测量,如汽轮机末端1~2级抽汽、低压缸排汽。湿蒸汽湿度的准确测量现在较难,就算利用先进设备测量准确,而将其设置在汽机设备上进行长时间测量的成本会很高。又比如:实时煤耗计算涉及的锅炉效率计算,需要实时的知道入炉煤元素分析,燃煤电厂通常只进行入炉煤工业分析化验,故这方面的难题难以克服。
发明内容
本发明的目的是:使得燃煤电厂获得实时煤耗等经济性指标,为电厂节能降耗带来更大的空间。
为了达到上述目的,本发明的技术方案是提供了一种用于燃煤电厂实时计算煤耗等经济指标的方法,其特征在于,包括以下步骤:
第一步、利用燃煤电厂自身的SIS系统或DCS系统将数据实时采集进入就地服务器后,整理成标准数据库;
第二步、在标准数据库基础上,对汽轮机回热系统进行热平衡计算,同时将锅炉侧减温水流量的热平衡纳入回热系统一起进行热平衡计算,从而计算出包括全水分Mar在内的涉及的汽水流量,全硫分Sar根据烟气进出口SO2含量的测量值估算得到;
第三步、在第二步已知各汽水流量参数的基础上,计算锅炉输出热量;
第四步、至少计算全厂管道效率、厂用电率、汽轮机热耗率;
第五步、脱硫单元运行参数估算入炉煤全硫分;制粉系统单元运行参数计算入炉煤全水分;
第六步、将第三步及第四步整理完善后,建立主方程组,主方程组包括:理论空气容积计算公式;理论烟气容积计算公式;Car、Har、Oar、Nar、Sar、Mar、Aar组份和等式;组分元素拟合发热量公式;锅炉效率正、反平均计算模块;组分元素计算烟气焓模块;Car、Har、Oar、Nar、Aar分别为收到基的全碳分、全氢分、全氧分、全氮分、全灰分;
第七步、通过主方程组计算得到锅炉效率,再结合上述第四步的结果,根据标准计算汽机热耗率、供电煤耗;
第八步、至少将供电煤耗、中间结果锅炉效率和重要流量参数、汽轮机热耗率通过就地服务器返回SIS系统或DCS系统。
本发明将机组性能试验标准、燃煤电厂SIS数据系统这两部分作为基础,提出燃煤机组实际运行条件下的实时煤耗等经济性指标计算方法。其中,仅依靠凝汽器排汽后的凝结水流量作为唯一流量测量基准值,依次建立汽轮机热耗率、锅炉效率与管道效率的实时计算模型,使得燃煤电厂获得实时煤耗等经济性指标,为电厂节能降耗带来更大的空间。
附图说明
图1为本发明的流程图,图中:Err为精度,定义为:
Abs[Q′net,ar (k+1)-Q′net,ar (k)]<0.00001&Abs[η′net,ar (k+1)-η′net,ar (k)]<0.00001,式中,Abs是求本次计算和上一次计算之间的差的绝对值,k为迭代次数。
具体实施方式
为使本发明更明显易懂,兹以优选实施例,并配合附图作详细说明如下。
本发明利用燃煤电厂自身已有的SIS数据系统中锅炉、汽轮机等相关测点的数据,进行煤耗、设备效率等经济性指标参数的计算。汽机热耗率、管道效率的计算根据美国机械工程师标准(ASME标准)所载公式进行。
锅炉效率的计算是利用正平衡、反平衡建立两个方程模块,再结合上述汽轮机热耗率计算过程中涉及到的参数和公式,建立方程组,寻找其余方程,一起求解入炉煤元素分析和锅炉效率。除了飞灰、大渣含碳量必须需要手动输入以外,锅炉效率就能通过方程组求解。
首先确定所求的未知数及其个数,Car、Har、Oar、Nar、Sar、Mar、Aar、Qnet,ar、η,未知数为9个,分别是收到基C、H、O、N、S、全水、灰、低位发热量和锅炉效率;其次简化未知数。对于全水分Mar,根据权威教科书上通过制粉系统热平衡计算,即可求出原煤全水分。对于全硫分Sar,目前300MW及以上的电站含脱硫设备投入运行,而且烟气进出口SO2的含量都有测量值,尽管Sar包含有机硫、无机硫,不是全部的S元素转化为SO2,但是由于Sar在整个煤元素分析中所占分量较小(一般1%以下),所以可以根据SO2含量进行估算;最后确定所需求解的未知数是:Car、Har、Oar、Nar、Aar、Qnet,ar、η,共计7个。
本发明提供了一种用于燃煤电厂实时计算煤耗等经济指标的方法,包括以下步骤:
第一步、利用燃煤电厂运行数据库(SIS系统)进行数据(上百个压力、温度等)实时采集进入就地服务器(可以1分钟1个点甚至30秒1个点的采集频率),整理成标准数据库。
第二步、根据标准,在标准数据库基础上,对汽轮机回热系统进行热平衡计算,同时将锅炉侧减温水流量的热平衡纳入回热系统一起进行热平衡计算,也即,锅炉减温水流量是计算值而不是测量值。这一步就将涉及的汽水流量都计算出来,比如抽汽流量、主蒸汽流量、减温水流量、给水流量等。作为所需的已知汽水流量参数只有一个就是除氧器入口给水流量,这就提升了整体参数的准确性。
第三步、在第二步已知各汽水参数的基础上,计算锅炉输出热量,为主方程组计算做准备。
第四步、计算全厂管道效率、厂用电率、汽轮机热耗率等参数。
第五步、脱硫单元运行参数估算入炉煤全硫分;制粉系统单元运行参数计算入炉煤全水分(结合美国机械工程师标准——ASME标准)。
第六步、将第三、四步整理完善后,进行主方程组的计算,主方程组见下表。
表1构建主方程组
表1中,Iy是烟焓;α是过量空气系数;εa是各级受热面处的过量空气系数不同而对上述公式的修正系数;Cy rθy是通用烟气焓,即1m3通用烟气焓由0℃-θ℃的平均定压容积热容与烟气温度的乘积。上述7个方程7个未知数构成的方程线性无关,该线性方程组可求得唯一解。
第七步、通过主方程组计算得到锅炉效率,再结合第上述第四步的结果,根据标准计算供电煤耗。
第八步、将最终结果供电煤耗、中间结果锅炉效率、汽轮机热耗率等重要指标参数通过就地服务器返回SIS数据库系统,接入DCS,实现实时对接。

Claims (1)

1.一种用于燃煤电厂实时计算煤耗等经济指标的方法,其特征在于,包括以下步骤:
第一步、利用燃煤电厂自身的SIS系统或DCS系统将数据实时采集进入就地服务器后,整理成标准数据库;
第二步、在标准数据库基础上,对汽轮机回热系统进行热平衡计算,同时将锅炉侧减温水流量的热平衡纳入回热系统一起进行热平衡计算,从而计算出包括全水分Mar在内的涉及的汽水流量,全硫分Sar根据烟气进出口SO2含量的测量值估算得到;
第三步、在第二步已知各汽水流量参数的基础上,计算锅炉输出热量;
第四步、至少计算全厂管道效率、厂用电率、汽轮机热耗率;
第五步、脱硫单元运行参数估算入炉煤全硫分;制粉系统单元运行参数计算入炉煤全水分;
第六步、将第三步及第四步整理完善后,建立主方程组,主方程组包括:理论空气容积计算公式;理论烟气容积计算公式;Car、Har、Oar、Nar、Sar、Mar、Aar组份和等式;组分元素拟合发热量公式;锅炉效率正、反平均计算模块;组分元素计算烟气焓模块;Car、Har、Oar、Nar、Aar分别为收到基的全碳分、全氢分、全氧分、全氮分、全灰分;
第七步、通过主方程组计算得到锅炉效率,再结合上述第四步的结果,根据标准计算汽机热耗率、供电煤耗;
第八步、至少将供电煤耗、中间结果锅炉效率和重要流量参数、汽轮机热耗率通过就地服务器返回SIS系统或DCS系统。
CN201610254860.7A 2016-04-22 2016-04-22 一种用于燃煤电厂实时计算煤耗等经济指标的方法 Active CN105956370B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610254860.7A CN105956370B (zh) 2016-04-22 2016-04-22 一种用于燃煤电厂实时计算煤耗等经济指标的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610254860.7A CN105956370B (zh) 2016-04-22 2016-04-22 一种用于燃煤电厂实时计算煤耗等经济指标的方法

Publications (2)

Publication Number Publication Date
CN105956370A CN105956370A (zh) 2016-09-21
CN105956370B true CN105956370B (zh) 2019-02-05

Family

ID=56916559

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610254860.7A Active CN105956370B (zh) 2016-04-22 2016-04-22 一种用于燃煤电厂实时计算煤耗等经济指标的方法

Country Status (1)

Country Link
CN (1) CN105956370B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109460885B (zh) * 2018-09-11 2023-01-20 广东粤电靖海发电有限公司 一种热力发电厂能量平衡单元的发电能耗评价方法
CN109726460B (zh) * 2018-12-21 2022-11-29 国网浙江省电力有限公司电力科学研究院 一种等效煤质的获取方法
CN112488435A (zh) * 2019-09-11 2021-03-12 范文 火力发电企业大数据实时正平衡供电煤耗计算与控制方法
CN110594780B (zh) * 2019-09-29 2021-02-26 上海发电设备成套设计研究院有限责任公司 一种燃煤电厂锅炉在线实时燃烧优化技术方法
CN111400875B (zh) * 2020-03-02 2023-09-19 神华神东电力有限责任公司 一种汽轮机组运行经济性评价方法及系统
CN112417685A (zh) * 2020-11-20 2021-02-26 西安热工研究院有限公司 一种计算汽轮机通流改造后最终考核热耗率的方法
CN113253018A (zh) * 2021-04-20 2021-08-13 苏州西热节能环保技术有限公司 一种燃煤电厂电能平衡监测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1369712A (zh) * 2002-03-15 2002-09-18 山东中实股份有限公司 电站锅炉入炉煤质实时监测方法及装置
CN102494714A (zh) * 2011-11-11 2012-06-13 东南大学 电站锅炉效率与燃煤热值、灰分和水分的同步测算方法
CN104698149A (zh) * 2015-03-20 2015-06-10 常州英集动力科技有限公司 一种燃煤锅炉煤质在线软测量方法及测量系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1369712A (zh) * 2002-03-15 2002-09-18 山东中实股份有限公司 电站锅炉入炉煤质实时监测方法及装置
CN102494714A (zh) * 2011-11-11 2012-06-13 东南大学 电站锅炉效率与燃煤热值、灰分和水分的同步测算方法
CN104698149A (zh) * 2015-03-20 2015-06-10 常州英集动力科技有限公司 一种燃煤锅炉煤质在线软测量方法及测量系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
入炉煤元素分析和飞灰含碳量的软测量实时监测;常建平;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑 》;20070815;C042-387

Also Published As

Publication number Publication date
CN105956370A (zh) 2016-09-21

Similar Documents

Publication Publication Date Title
CN105956370B (zh) 一种用于燃煤电厂实时计算煤耗等经济指标的方法
CN102494714B (zh) 电站锅炉效率与燃煤热值、灰分和水分的同步测算方法
CN103697946B (zh) 一种火电厂燃煤锅炉烟气流量的计算方法及污染物排放量的控制方法
CN104698149B (zh) 一种燃煤锅炉煤质在线软测量方法
AU762836B2 (en) Input/loss method for determining fuel flow, chemistry, heating value and performance of a fossil-fired system
CN101697179A (zh) 基于正反热平衡关系的电站锅炉燃煤热值趋势测算方法
CN102734782B (zh) 一种燃煤锅炉能效监测方法
CN101825502B (zh) 汽机带疏水冷却器的加热器出水及疏水温度测算方法
CN101865867B (zh) 一种实时计算煤质元素和工业成分的方法
CN110375787A (zh) 一种计量仪表运行状态在线评价方法
CN106768068A (zh) 火电机组燃煤热值与发电标煤耗指标在线监测系统与方法
CN105091944A (zh) 一种火力发电厂机组燃煤发热值及煤耗率指标在线监测的方法
WO2005083447A1 (fr) Procede et dispositif de controle en ligne pour un appareil de conversion de combustible fossile
CN109785187B (zh) 一种发电机组供电煤耗检测数据校正方法
CN111753389B (zh) 一种燃煤锅炉入炉原煤热值在线软测量方法和装置
CN108197723B (zh) 煤电机组供电煤耗与污染物排放的优化节能调度方法
CN112131517A (zh) 一种垃圾焚烧电厂入炉垃圾低位热值的测算方法
CN105223036A (zh) MgO烟气脱硫性能现场验证方法与系统
CN108073763B (zh) 一种电站锅炉飞灰含碳量的测量方法
CN103267539B (zh) 一种卧式三段式给水加热器上端差及下端差的测算方法
CN106096763B (zh) 一种预测运行机组的燃煤发热值的方法
CN103697958A (zh) 燃煤机组汽包出口饱和蒸汽质量流量的实时计量方法
CN106596156B (zh) 一种除尘器性能试验误差测评方法
Shi et al. Real-time online monitoring of thermal efficiency in coal-fired power plant boiler based on coal heating value identification
CN113836794A (zh) 一种软硬结合的飞灰含碳量在线监测方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant