本发明涉及一种烟气流量测量方法，具体涉及一种燃煤电站锅炉尾部双烟道烟气流量在线软测量方法，属于电力和信息处理
技术领域：
。
背景技术：
：进入二十一世纪以来，节能和环保愈发重要。燃煤火力电站作为关系着我国国计民生的能源消耗大户，如何在保证安全生产的基础上尽可能的减少能源消耗、减少排放已经成为其面临的一个十分迫切的问题。锅炉是燃煤电站的关键设备之一，也是整个电站节能减排的重点对象。无论采取何种手段对锅炉进行节能减排优化，对锅炉整体工况的准确监测是重要前提，这其中烟气流量是其中需要监测的变量之一。我国现有燃煤电站其中相当一部分锅炉是采用尾部双烟道设计，两侧烟道的烟气挡板对烟气流量进行调节，其主要目的是通过对烟气量的调节达到对再热器蒸汽温度的调节作用。尾部双烟道的烟气比例是事先设计好的，但由于实际工况的复杂性，诸如烟气成分、烟气流速等因素的影响，双烟道的实际烟气流量并不是一直遵循比例的。要对两侧烟道的受热面工作状况进行实时监测则必须获得两侧烟道的实时烟气流量。测量烟气流量一般方法是采用烟气流量计，在锅炉设计安装阶段事先安装在相应的部位，或是在试验阶段采用测量仪器在预留的测量点进行测量。但在我国很多采用尾部双烟道设计的锅炉尾部烟道没有预先设计或者预留测点，这就导致了无法直接测量获取尾部两侧烟道实时的烟气流量，对两侧受热面的实时热力学计算造成困难。如若采用实现设计的比例来对烟气流量进行分配则不够准确，对于一些对精度要求较高的控制和优化方法则无法达到要求。常规的方法是根据两侧烟气挡板的开度信息来判断烟气流量的分配，然而由于压损系数难以获取，而且在长时间生产后，烟气挡板由于烟气的剥蚀变得密闭性不好且开度信息不准确，因而这种方法在实际现场中无法使用。经对现有文献检索发现，中国专利zl201210046212.4提出了一种“基于检测氧含量的插入式烟气流量测量方法及测量装置”，该方法从数据处理装置输入端输入燃料的各组分含量和燃料流量，计算得到该燃料的理论空气比和理论烟气比以及理论空气量和理论烟气量；通过氧含量检测装置测量探头探测得到烟气中的氧含量，最终计算出实际烟气流量。该方法仍然需要氧含量检测探头，而我国很多尾部双烟道锅炉的氧含量探头在双烟道的尾部，并不能完全准确反映烟道内的整体状况，对于本专利申请的所针对的尾部双烟道烟气流量的检测而言具有一定的局限性。王建国等人提出的“锅炉尾部烟道烟气质量流量的测量模型”以热力学平衡计算为基础，采用工控机和采集装置相结合的方式，实现了对管式空气预热器烟气质量流量的测量。该方法主要针对烟气质量流量的计算，对于燃煤波动较大的机组该方法计算出的锅炉效率与实际值存在着一定的偏差，不能满足更准确的更高一步的要求。技术实现要素：为解决现有技术的不足，本发明的目的在于提供一种基于动态能量和质量平衡的燃煤电站锅炉尾部双烟道烟气流量在线软测量方法，在机组实际运行数据的基础上，从进入尾部双烟道开始，尾部烟气挡板出来到空气预热器为止的这段烟道的一系列受热面采用烟气侧到工质侧的动态能量和质量平衡，通过软测量的方法对两侧烟道的烟气流量进行计算，同时采用每天煤质的工业数据对基准烟气物性参数进行校正，最终实现对尾部双烟道烟气流量的实时软测量。为了实现上述目标，本发明采用如下的技术方案：一种燃煤电站锅炉尾部双烟道烟气流量在线软测量方法，其特征在于，包括以下步骤：第一步：根据电站实际情况确定基准的燃煤成分及特性；选取固定的采样时间，通过opc通讯协议实时获取锅炉运行相关数据；第二步：基于烟气物理特性，根据第一步中确定的基准的燃煤成分，计算1kg燃煤燃烧后的产物的各项特性，并建立锅炉的基准烟气物性数据库；将水和水蒸气热力性质iapws-if97公式向量化用于工质焓值的实时计算；第三步：对第一步获取到的锅炉运行相关数据进行坏点处理和数据平滑处理，从而得到处理后的锅炉尾烟道的相关运行参数的实时数据库；第四步：建立基于质量和能量平衡的锅炉尾部两侧省煤器前后的烟气侧和工质侧的平衡方程，用如下公式表示：mlδhl+mrδhr＝qsm(2)ml+mr＝mtotal式中，ml、mr和mtotal分别为尾部竖井烟道左侧、右侧和总的烟气质量流量，kg/s；δhl和δhr分别为左右两侧烟道中单位质量烟气的焓值变化，kj/kg；qsm为锅炉尾烟道两侧省煤器中工质吸收的热量，kj/s；第五步：根据温度变化下工质和金属受热面的蓄热变化修正上述平衡方程，修正后的省煤器能量和质量平衡方程用如下公式表示：q′sm＝qsm+δqq+δqj(3)式中，q′sm为修正后的锅炉尾烟道两侧省煤器中工质吸收的热量，kj/s；δqq为工质蓄热变化量，kj/kg；δqj为受热面金属蓄热变化量，kj/kg；第六步：根据燃煤的工业分析数据校正燃煤的基准化学成分；第七步：以第六步得到的校正后的燃煤化学分析值为基准，从测点得到两侧烟道烟气的温度、压力、流量、工质温度以及压力流量的相关数据，结合公式(2)和公式(3)联立求得两侧烟道的烟气流量：前述的燃煤电站锅炉尾部双烟道烟气流量在线软测量方法，其特征在于，在第一步中，采样时间为5秒。前述的燃煤电站锅炉尾部双烟道烟气流量在线软测量方法，其特征在于，在第三步中，通过多项式滑动拟合的方法来判断坏点，并对其进行剔除。前述的燃煤电站锅炉尾部双烟道烟气流量在线软测量方法，其特征在于，在第三步中，通过数据加权滤波的方法来对数据进行平滑处理，具体是：式中，为滤波后的数据；yi为i时刻的实际测量值，i＝m-3,m-2,m-1,m,m+1,m+2,m+3。前述的燃煤电站锅炉尾部双烟道烟气流量在线软测量方法，其特征在于，在第六步中，校正燃煤基准化学成分采用如下的校正公式：式中，xi,ar为计算用碳、氢、氧、氮、硫元素成分分析，％；为季度化验的碳、氢、氧、氮、硫元素成分分析基准值，％；war为工业分析值，基水分，％；aar为工业分析值，基灰分，％；为季度化验的收到基水分，％；为季度化验的收到基灰分，％。本发明的有益之处在于：(1)本发明提供的在线软测量方法，可以较为准确的实时测量燃煤电站锅炉尾部双烟道烟气流量，具有良好的动态特性，同时不需要安装额外的测量仪器，具有较好的通用性和鲁棒性，适用于燃煤火力发电锅炉尾部烟道特征参数的实时监测以及在此基础上的优化调整；(2)本发明提供的在线软测量方法，对能量平衡方程进行了实时修正，考虑了工况波动的情形下受热面金属和工质蓄热变化，提高了方法的动态精度；(3)本发明提供的在线软测量方法，对入路燃煤化学成分进行校正，不需要昂贵的在线分析仪的投入的同时能够保证一定的精度；(4)本发明提供的在线软测量方法的测量结果有助于对锅炉尾部烟道工况进行更好的评估为总体的节能减排优化控制提供支持。附图说明图1是燃煤电站锅炉尾部双烟道结构示意图；图2是左侧烟道烟气流量份额图。具体实施方式以下结合附图和具体实施例对本发明作具体的介绍。第一部分：锅炉本发明所涉及的锅炉，需要满足以下条件：(1)负荷300mw，采用亚临界、自然循环、一次中间再热、双拱单炉膛、“w”火焰燃烧方式、尾部双烟道、烟气挡板调温、平衡通风、固态排渣、露天布置、全钢架悬吊式汽包炉；(2)换热器外部热流沿管长及圆周均匀分布；(3)换热器管壁金属只考虑径向换热；(4)管内工质不可压缩，且流动截面上参数均匀一致。作为一种优选的方案，我们使用hg-1025/17.3-wm18型锅炉，该锅炉的尾部双烟道结构示意图如图1所示。第二部分：双烟道烟气流量在线软测量方法本发明的燃煤电站锅炉尾部双烟道烟气流量在线软测量方法，其包括以下步骤：第一步：首先，设置锅炉尾烟道包括省煤器在内的换热面尺寸和烟道特性，具体如表1所示。表1锅炉物理特性参数名称单位数值省煤器管外径mm51省煤器管内径mm44.5横向排数(并联管数)——356纵向排数——8横向节距mm115烟气通流面积m2135.25蒸汽流通面积m20.448表面受热面积m22539壁管材质——sa210c然后，根据电站实际情况确定基准的燃煤成分及特性，具体如表2所示。表2电站燃煤代表性基准成分及特性最后，选取固定的采样时间(采样时间为5秒)，该锅炉dcs控制系统自带opc服务端，通过opc通讯协议实时获取锅炉运行相关数据。第二步：基于烟气物理特性，根据第一步中确定的现场的基准的燃煤成分，计算1kg燃煤燃烧后的产物的各项特性，利用高次多项式拟合各成分的焓值特性，建立锅炉的基准烟气物性数据库，用于后续的烟气焓值计算。对于烟气中的成分：co2、n2、空气、水蒸气以及灰的比热c，用正交多项式拟合成五次多项式计算：c＝a0+a1t+a2t2+a3t3+a4t4+a5t5式中，a0、a1、a2、a3、a4、a5为拟合系数；t为工质温度，℃。各个部分的拟合结果见表3。表3烟气及空气物性参数计算公式多项式系数注：表中气体体积均指标准状态。工质水和水蒸气的焓值根据国际标准公式iapws-if97求得。第三步：对第一步获取到的锅炉运行相关数据进行坏点处理和数据平滑处理，将处理后的锅炉尾烟道的相关运行参数数据存入锅炉运行参数实时数据库。1、坏点处理通过多项式滑动拟合的方法来判断坏点，并对其进行剔除。2、数据平滑处理通过数据加权滤波的方法来对数据进行平滑处理，具体是：式中，为滤波后的数据；yi为i时刻的实际测量值，i＝m-3,m-2,m-1,m,m+1,m+2,m+3。第四步：基于质量和能量平衡，建立锅炉尾部两侧省煤器前后的烟气侧和工质侧的质量和能量平衡方程，尾部烟道左右两侧的烟气流量之和等于尾部烟道出口总的烟气流量，省煤器烟气释放的热量等于工质吸收的热量，具体用如下公式表示：mlδhl+mrδhr＝qsm(2)ml+mr＝mtotal式中，ml、mr和mtotal分别为尾部竖井烟道左侧、右侧和总的烟气质量流量，kg/s；δhl和δhr分别为左右两侧烟道中单位质量烟气的焓值变化，kj/kg；qsm为锅炉尾烟道两侧省煤器中工质吸收的热量，kj/s。第五步：考虑工况变化情况下，金属受热面和工质在温度、压力变化下蓄热量的变化，对第四步的能量平衡方程进行修正，得到动态能量平衡方程，修正后的省煤器能量和质量平衡方程用如下公式表示：q′sm＝qsm+δqq+δqj(3)式中，q′sm为修正后的锅炉尾烟道两侧省煤器中工质吸收的热量，kj/s；δqq为工质蓄热变化量，kj/kg；δqj为受热面金属蓄热变化量，kj/kg。其中，受热面金属蓄热变化量δqj为：工质(蒸汽)蓄热变化量δqq为：而蒸汽侧吸热则为：qsm＝d(hout-hin)式中，cj、cq分别为金属和工质的平均比热容，kj/(kgg℃)；mj、mq为受热面管壁金属质量和其内工质的质量，kg；θq、θj为金属管壁温度和蒸汽温度，℃；d为流经两侧省煤器的工质质量流量，kg/s；hout、hin分别为省煤器进、出口工质焓值，kj/kg；τ为时间。第六步：根据燃煤的工业分析数据校正燃煤的基准化学成分。本发明所涉及的锅炉机组，有每天4班的燃煤工业分析数据，可以根据该数值对燃煤的基准化学成分进行校正，校正燃煤基准化学成分采用如下的校正公式：式中，xi,ar为计算用元素成分分析(碳、氢、氧、氮、硫)，％；为季度化验的元素成分分析基准值，(碳、氢、氧、氮、硫)，％；war为工业分析值，基水分，％；aar为工业分析值，基灰分，％；为季度化验的收到基水分，％；为季度化验的收到基灰分，％。第七步：以第六步得到的校正后的燃煤化学分析值为基准，从测点得到两侧烟道烟气的温度、压力、流量、工质温度以及压力流量等相关数据，结合公式(2)和公式(3)联立求得两侧烟道的烟气流量：根据锅炉的设计参数，额定工况下左右两侧烟道的烟气份额约为3：7，采用现场实际运行数据由上述方法计算所得的左侧烟气份额如图2所示。由图2可以看出，在额定工况阶段的计算结果与设计值较为吻合。需要说明的是，上述实施例不以任何形式限制本发明，凡采用等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。当前第1页12