张振义，邹 磊

(1.国华徐州发电有限公司，江苏徐州 221166；2.江苏方天技术有限公司，江苏南京 211102)

火力发电厂烟气排放流量是环保监测的重要数据[1]，烟气在线监测设备大多安装在烟囱前的烟道上。由于烟道内烟气流场复杂，随意布置的流量测点无法合理反映实际流量，烟气流量的测量精度无法达到规范要求。燃煤机组改建脱硫工程后，由于场地空间限制，脱硫系统后净烟气通往烟囱的烟道较短，且存在多处拐角，这些均给烟气流量测点的安装带来了困难。由于测点安装位置不当，在线监测设备监测到的流量值往往波动较大、精度较低，有的甚至出现较长时间无流量值的情况，而且在不同工况下，由于烟气量的不同，相同位置的流场分布也可能会发生改变。若测点安装在涡流区还会引起监测设备的强烈震动[2]。目前，火电厂普遍存在由于烟气流量测点安装位置不当造成测量值准确性差的问题。数值模拟作为一种成熟的研究技术在电站锅炉流场分析中得到了广泛的应用[3，4]。

1 计算模型及边界条件

1.1 计算模型

数值计算的物理模型采用实际的净烟气烟道的结构数据，计算区域为从脱硫系统的烟气换热器（GGH）出口连接处到烟囱烟道连接处的净烟气烟道的连接面，计算模型如图1 所示。从图1 中可看出，烟气自GGH 流出后至下行烟道，经转弯至下部水平烟道，再转弯至上行烟道，进一步转弯至上部水平弯曲烟道，最后至烟囱烟道连接处。

计算模型采用烟道的实际结构尺寸，其中的坐标轴为：X 向基本为烟气的主流方向，零点设定在GGH出口烟道的轴线处；Y 向为高度方向，坐标值采用工程实际的标高数值；Z 向为烟道的深度方向，零点设定在下部水平烟道的轴线位置。本文中所指的下部烟道包括了GGH 后下行烟道、转弯、下部水平烟道及进一步转弯后的上行烟道；上部水平弯曲烟道则主要是指下部上行烟道转水平后的、为适应烟囱位置的水平弯曲烟道。

图1 净烟气烟道流场计算区域示意图

采用k-ε 双方程湍流模型来模拟烟道内气相湍流流动，其控制方程包括连续性方程、动量方程、能量方程、k 方程、ε 方程。若不考虑热交换的单纯流场问题，则不需要包含能量方程[5]。上述控制方程可表示为如下通用形式：

在三维直角坐标系下，k-ε 双方程模型的控制方程如下：

式（1，2）中：Φ为待求解的物理量；Γ为扩散系数；div（ρVΦ）为对流项；div（ΓgradΦ）为扩散项；S为源项。对于特定意义的Φ，具有特定的量Γ 和S。在通用形式方程中，当Φ 分别为1，U（u，v，w），k，ε时，分别对应连续性方程，动量方程，k 方程和ε 方程。

1.2 边界条件

按照设计烟气量1 012 604 Nm3/h 和50%设计烟气量共进行了2个工况条件的数值计算，烟气温度按照设计值80℃计算。入口采用流量边界条件，出口采用压力边界条件，压力按照设计值给定，设定烟气流动与壁面无滑移且壁面按绝热处理。收敛残差标准设定为10-4。为简化计算，视烟气为单一气体，不考虑不同烟气成分对计算结果的影响。

2 数值计算结果

2.1 净烟气烟道内流场特性分析

净烟气烟道在100%设计烟气量工况下的流场分布如图2、图3 所示。其中，图2为烟道流场速度分布立体图，以下部水平烟道轴线垂直面和上部水平弯曲烟道的3个不同标高水平面上的速度分布来反映整个烟道内的速度分布情况；图3为净烟气烟道内流场速度矢量分布图。

图2 设计烟气流量工况下烟道流场速度分布立体图

图3 烟道流场速度矢量分布图

从图2 和图3 可以看出，烟道内的速度场整体情况：GGH 后烟道下行转弯进入下部水平烟道后，烟气速度趋于均匀；烟气继续流动经水平转垂直转弯后，在上行过程中，除弯头局部出现低速区域外，速度也较均匀；但在烟气再转水平流动时，由于垂直转水平转弯和为适应烟囱布置的烟道水平弯曲的存在，导致这水平转弯烟道开始前的上部转弯区（图2、图3 中的左上角）出现一个低速回流涡流区；此后由于烟道的2个连续转弯和上部水平弯曲烟道的尺寸较短，导致整个上部水平弯曲烟道内速度分布严重不均匀。

净烟气烟道局部截面上速度分布情况如图4—6所示。其中，图4为下部烟道轴线垂直面（Z=0）上的速度分布，图5为上部水平弯曲烟道轴线水平面（标高Y=12 000 mm）上的速度分布，对应下部烟道，图6 则为下部水平烟道中部（X=7200 mm 处）垂直X 轴的截面、上行烟道（X=11 526 mm 处）垂直X 轴的截面上和烟道出口垂直面上的速度分布，以反映烟道在Z 向上速度分布的均匀性。

图4 下部烟道轴线垂直面（Z=0）上的速度分布

图5 上部水平转弯烟道轴线水平面（Y=12 000 mm）上的速度分布

图6 下部水平烟道（X=7200 mm）和上行烟道（X=11 526 mm）处垂直截面上和烟道出口的速度分布

从图4 和图5 可以看出，在上部水平弯曲烟道内侧、上部出现高速区，而外侧、下部出现低速区。从图6可以看出，为了适应烟囱的布置，通往烟囱的弯曲烟道尺寸较短，从而导致这种不均匀分布一直延续到出口面。对于下部烟道及其后转弯上行烟道，烟气速度分布基本是均匀的。分析100%设计烟气量工况下的流场分布可知，净烟气烟道下部速度分布较均匀，速度最均匀的区域出现在该段的中部水平烟道区域，在该区域内烟气速度变化较小，在15～18 m/s 范围内；上部水平弯曲烟道内烟气经历3 次转向后，在烟道垂直截面上的速度分布严重不均；烟道中各弯头的圆弧过渡处都不同程度的出现了局部低速或高速区。

2.2 不同烟气量对烟道内速度分布的影响

为了反映不同烟气量对净烟气烟道内流场分布的影响，本文对50%设计烟气量工况下烟道内流场进行数值模拟。该工况下烟道流场速度分布立体图如图7所示。为与图2 进行对照比较，图7 同样给出下部水平烟道轴线垂直面和上部水平弯曲烟道的3个不同标高水平面上的速度分布情况。

对于50%设计烟气流量的计算工况，计算结果显示，除速度的数值外，计算得到的烟道内速度分布总体上与设计烟气流量工况（图2）几乎一致，可以看出，烟气流量的变化对烟道内的速度分布影响很小。

3 结束语

（1）下部水平烟道速度分布相对较均匀，速度最均匀的区域出现在该段的中部水平烟道区域；上部水平弯曲烟道内烟气经历3 次转向，在烟道垂直和水平方向的速度分布都不均匀。

图7 50%设计烟气流量工况下烟道流场速度分布立体图

（2）烟气流量的变化对烟道内的速度分布影响不大，烟气流量从设计流量减少50%，除速度数值外，烟道内的速度分布几乎一致。

（3）通过对某电厂脱硫系统后净烟气烟道内流场的数值计算分析发现，在该烟道下部水平烟道段中部安装烟气流量测点相对较好。

（4）通过数值模拟计算结果可以合理选择烟气流量测点，从而提高烟气流量测量精度。

[1]王万林，齐小娟.火电厂烟气排放流量测量方法研究[J].浙江电力，2009(6)：17-19.

[2]白江文，魏 威，周 强，等.烟气排放连续监测系统及其常见故障分析处理[J].江苏电机工程，2011，30(3)：78-80.

[3]邹 磊，高小涛，黄 磊.基于CFX的炉内空气动力场的数值模拟研究[J].江苏电机工程，2009，28(4)：81-84.

[4]陈 雯，于向军.基于FLUENT的风粉管道风速流场模拟[J].江苏电机工程，2008，27(6)：72-75.

[5]王福军.计算流体动力学分析CFD 软件原理与应用[M].北京：清华大学出版社，2004：7-11.