基于塔式扩散模型的火电厂污染物排放监测显示方法
基于塔式扩散模型的火电厂污染物排放监测显示方法
【技术领域】
[0001] 本发明涉及监控领域，特别涉及污染物在大气中的排放监测，用于在电子地图上 显示火电厂排放的烟气污染物在各区域的浓度。
【背景技术】
[0002] 环境污染问题日益严重，雾霾天气不断威胁人类的健康。近年来，国家电网公司 积极落实相关政策要求，大力开展了污染物排放系统功能平台的建设工作，实现对固定污 染源排放的污染物进行连续地、实时地跟踪测定，连续测定颗粒物、气态污染物浓度和排放 率。
[0003] 现在的技术，只是通过表格或曲线、棒图显示当前污染物排放量，不能直观地反映 污染物排放对周边环境的影响。
[0004] 大气污染物扩散模式是模拟大气污染物的扩散、迀移过程，预测在不同污染源条 件、气象条件及下垫面条件下某污染物浓度时空分布的数学模型，是低层大气中污染物迀 移和扩散规律的、简单化的数学描述。根据不同的建模理论体系、污染物迀移、扩散过程以 及不同的描述对象，模式的形式也各不相同。
[0005] 受地形、气象、大气污染物的物理化学特征、污染源特征等多种因素的制约，不同 的扩散模式都有各自不同的考虑因素和适用范围，选择恰当的扩散模式能较为准确地模拟 污染物的扩散及分布。目前，已经有许多发展成熟的污染物扩散模式应用，国内的相关研宄 中尤其以高斯类模式的应用最为广泛。
[0006] 高斯扩散模型一种经验模型，大量实验数据表明高斯分布是污染物在分布的较好 的近似形式，高斯模型数学形式简单，通过大规模的实验获得模型计算的主要参数，可以在 各种气象条件下选用，高斯扩散模型有两个特点：浓度计算在垂直方向和水平方向都采用 高斯分布假设，湍流分类和扩散参数采用离散化的分类方法，它的优点是仅用常规气象观 测资料就可以评估扩散参数。
[0007] 高斯扩散模型基于以下条件： 坐标系：原点为地面源的排放点或高架源排放点在地面上的投影，右手坐标系，x轴为 主风向，y轴为横风向，z轴为垂直向。
[0008] 几点假设：1、污染物浓度在y、z方向为正态分布，2、全部高度风速均匀稳定，3、源 强连续均匀，4、扩散中不考虑转化，污染物是守恒的，5、x向上，平流作用远大于扩散作用， 6、地面足够平坦。
[0009] 无界空间连续点源扩散模式为： c(x,y,z) =q/ (2 31y〇y 〇 z)exp[- (y2/2 0 y2+z2/2 0 z2)] 高架连续点源扩散模式为： c(x,y,z,H) =q/(2 31y〇y 〇 z)exp(-y2/2 〇 y2) *{exp[-(z-H) 2/2 0 y2] +exp[-(z+H)2/20z2)]} 其中，y为平均风速，q为源强，〇y、〇z为侧向和竖向扩散参数，H为有效烟囱高度。
[0010] 具体到火电厂烟气污染物排放，由于火电厂的特殊性，如污染物在出口的温度较 高、污染物的组成比较复杂以及气压等因素，上述扩散模式预测此类污染物的扩散会产生 较大的偏差。
【发明内容】
[0011] 为了更加准确地预测火电厂烟气污染物的排放，实时在电子地图上显示污染物浓 度，提出了本发明。
[0012] 为了达到上述目的，本发明采用的技术方案是：基于塔式扩散模型的火电厂污染 物排放监测显示方法，基于火电厂污染物排放检测系统实现，系统中包括主机，GIS服务器 和数据采集系统，数据采集系统连接CEMS系统，关键在于：主机中存储当地全年的气象资 料、每个烟囡的位置和高度信息；在每个烟囡内部出口位置设置温度传感器，在室外设置温 度传感器，在每个烟囡顶端外侧设置风向风速传感器，所有传感器编号并连接到数据采集 系统，主机采用塔式扩散模型计算污染物的扩散，塔式扩散模型的数学表达式为：
其中：x、y、z为空间坐标，H是烟囱口距地面的高度，巾为扩散因数，y为风速，t为扩 散时间，A为源强，a为颗粒属性系数，0为高度系数，空间坐标系是右手坐标系，原点是 烟囱在地面上的投影，x轴为主风向，y轴为横风向，z轴为垂直向； 主机计算各时间点不同位置的污染物浓度，并根据污染物浓度等级，以不同颜色显示 在电子地图上。
[0013] 颗粒属性系数a如下表：
根据高度得到的高度系数9如下表：
扩散因数巾根据温度、气压、风速等多方面因素影响来确定，本算法按照高度系数0 来划分不同高度时的扩散因数：
1、塔式扩散模型的推导过程： 假定由污染源排放出的大气污染物全都由单位面元dydz中通过，根据质量守恒定律 和连续性原理，单位时间内通过单位面源的污染物就等于污染源单位时间内排放出的污染 物质量，故可用公式表示如下：
式中，污染源排放的污染物浓度，Q:污染源单位时间内排放的污染物质量。
[0014] 将高斯曲线的参数代入得到：
将其进行转换可得：
根据普哇松积分公式得到
将该值代回得
在考虑到高度以及风速的影响，最终可得到任一点的污染物浓度公式为：
小为扩散因数，U为风速，t为扩散时间，A为源强，a为颗粒属性系数，0为高度系 数。
[0015] 2、上升高度AH应用国家标准，定义如下： 烟囱的高度为实际测量值， 当QH彡 2100kW和（TS-Ta)彡 35K时， Qh=0. 35PaQvAT/TS,AT=TS-Ta AH=n0QHnl.Hsn2.y-1 QH为烟囱的热排放率，这是一个固定参数，存储在主机中，Qv是实际排烟率，由CEMS系 统中的烟气流量计量仪器获得，单位是m3/S。
[0016] y为烟囱出口处的平均风速，1；和Ta分别为烟气出口温度和环境大气温度，n^为 烟气热机地表状况系数，^为烟气热释放率指数，n2为烟囡高度指数，Pa为大气压力，大气 压力根据烟囱出口处的高度确定。下表为系数的选取：
CEMS是英文ContinuousEmissionMonitoringSystem的缩写，是指对大气污染源排 放的气态污染物和颗粒物进行浓度和排放总量连续监测并将信息实时传输到主管部门的 装置，被称为"烟气自动监控系统"，亦称"烟气排放连续监测系统"或"烟气在线监测系统"。 从该系统中，可以实时得到排放的污染物中的成分和源强。
[0017] 采用本发明，使用适合火电厂污染物排放的塔式扩散模型，能够更加准确地计算 火电厂烟气污染物的排放，根据计算结果，在电子地图上实时显示烟囱口周围的污染物浓 度，同时，还可以预测污染物的扩散，模拟显示污染物扩散过程。
【附图说明】
[0018] 图1为本发明的系统组成框图， 图2为柱状图显示污染物浓度的示意图， 图3为曲线图显示污染物浓度的示意图。
[0019] 其中，1为测量室外温度的温度传感器，2为烟囱内部出口位置的温度传感器，3为 风向风速传感器，5为烟囡。
【具体实施方式】
[0020] 基于塔式扩散模型的火电厂污染物排放监测显示方法，基于火电厂污染物排放检 测系统实现，系统中包括主机，GIS服务器和数据采集系统，数据采集系统连接CEMS系统。 GIS服务器用来存储和显示电子地图，数据采集系统采集算法需要的各种参数。
[0021] 塔式扩散模型的数学表达式为：
其中：x、y、z为空间坐标，H是烟囱口距地面的高度，巾为扩散因数，y为风速，t为扩 散时间，A为源强，a为颗粒属性系数，0为高度系数，空间坐标系是右手坐标系，原点是 烟囱在地面上的投影，X轴为主风向，y轴为横风向，z轴为垂直向。
[0022] 算法中，y的取值是平均风速。一般的，是按照统计规律，从气象资料中取得当前 日期的平均风速，主机中存储当地全年的气象资料就是为了取得上述数值；另外，主机还存 储每个烟囱的位置和高度信息；在每个烟囱内部出口位置设置温度传感器测量烟囱出口温 度，在室外设置温度传感器测量室外温度，在每个烟囱顶端外侧设置风向风速传感器测量 实时风向风速。
[0023] 火电厂烟囱排放的污染物是气体和颗粒的混合物，如气体的二氧化硫、氮氧化物， 固体颗粒的烟尘等。CEMS系统可以获得排放物的种类和源强。
[0024] CEMS系统主要测量S02、NOx、02、粉尘的排放强度，也就是源强。这里的粉尘大小 是取颗粒的平均尺寸。
[0025] 主机获得实时参数，根据塔式扩散模型，计算各时间点不同位置的污染物浓度，并 根据污染物浓度等级，以不同颜色显示在电子地图上。
[0026] 为了区分不同颗粒大小的污染物，算法中使用了颗粒属性系数a。颗粒属性系数 a与污染物颗粒的直径相关：直径为0-14Mm颗粒属性系数为0,直径为15_3〇Mm颗粒属性 系数为0. 3,直径为31-47Mm颗粒属性系数为0. 5,直径为48-75Mm颗粒属性系数为0. 8,直 径大于76Mm颗粒属性系数为1 ;气体污染物按小于15um的颗粒计算。在现有技术基础上， 只计算平均尺寸的颗粒污染物的排放，也就是说按照从CEMS系统中得到的数据进行计算。
[0027] 本发明中，还充分考虑了高度对污染物扩散的影响，增加了高度系数0和与 高度相关的扩散因数:高度系数9与高度z的关系如下：H+AH>z多H，0为1;H >z彡 0? 75H，9 为 0? 75 ;0? 75H>z彡 0? 5H，9 为 0? 5 ;0? 5H>z彡 0? 25H，9 为 0? 25 ;0? 25H >z彡 0，9 为 0? 1〇
[0028] 扩散因数巾与高度z的关系如下： H+AH>z>H，<i)S(3 0/2)yto-87; H>z彡 0? 75H，<i> 为（0 /2)yt0.84; 0? 75H>z彡 0? 5H，<i> 为（9 /2)yt0.82; 0? 5H>z彡 0? 25H，<i> 为（9 /5)yt0.8; 0? 25H>z彡 0, <i> 为（9 /6)yt0.78。
[0029] 上述公式中，H是烟囱口距地面的高度，AH是排放物从烟囱口出来后上升的高度， y为风速，t为扩散时间。
[0030] 风速y为平均风速，是污染物在扩散过程中的平均值。风向和平均风速一般从气 象资料中获得。为了使结果更精确，在具体计算时，数据采集系统通过风向风速传感器获得 风速y和主风向，开始计算时使用获得的主风向作为x轴，气象资料中的平均风速作为y。 连续取得实时风速并计算平均值，经过一段时间后，使用实时取得的平均风速作为y，并在 其后的计算过程中不断修正。
[0031] 火电厂烟囱排放的污染物是气体和颗粒的混合物，如气体的二氧化硫、氮氧化物， 固体颗粒的烟尘等。在计算各时空点的污染物浓度时，不能只计算单一的污染物浓度，要按 照不同尺寸的污染物分别计算。由数据采集系统连 接CEMS系统获得排放污染物中不同尺 寸的污染物和各自的源强A，将探测到的参数代入扩散模型分别计算。目前只是计算S02、 NOx和平均粉尘尺寸的污染物浓度，当检测技术可以实时得到排放的污染物中其他气体和 粉尘中各种尺寸污染物的源强时，本算法会得到更精确的结果。
[0032] 将污染物浓度显示到有电子地图的屏幕上时，使用离散点的浓度值来代表某一区 域的浓度。为了在计算量和显示效果之间取得平衡，将显示的区域等分成N*M个网格，每个 网格的污染物浓度按照其中心点的污染物浓度显示，或者在每个网格内选取3~9个点，分 别计算选取点的污染物浓度，取其平均值，每个网格的污染物浓度按照平均值显示，或者每 个网格内选取5个点，位置分别为网格的中心及四角。N的取值范围是10~50,M的取值范 围是10~30。根据实际情况，N和M的取值可以更大。
[0033] 由于污染物由多种成分组成，因此，由数据采集系统从CEMS获得排放污染物中的 成分和各自的源强A，在计算某一位置的污染物浓度时，按照污染物中不同成分分别计算。 [0034] 显示时，可以分别显示单独的污染物浓度，也可以将不同成分的污染物浓度加权 求和为综合浓度，显示综合浓度。权重的比例分配可以按需要设定，如302分配权重50%、N0x30%，粉尘20%。加权求和的过程为各成分的污染物浓度乘以权重后相加。
[0035] 为了精确显示污染物浓度，显示时要求有高度参数。设定高度参数，按照高度参数 进行计算，并显示设定高度的计算结果。这里实现了分时分层显示，如显示地表的污染物分 布情况、显示与人体等高的1. 7米处污染物分布情况、显示离地高度5米的污染物分布情况 等。
[0036] 另外，还可以不设定高度参数，系统根据风速风向等因素，自动计算各位置污染物 浓度最大的值并显示。如在烟囱出口附近，与烟囱等高处污染物浓度最大，x轴方向20米 位置，低于烟囱出口 2米处污染物浓度最大，x轴方向200米位置，低于烟囱出口 4米处污 染物浓度最大，选取上述最大值并显示。
[0037] 在显示一定范围内污染物分布的同时，可以在地图上选取一个点，计算该点不同 高度的各种污染物浓度，使用柱状图或曲线图显示，高度间隔及高度范围人为设定。
[0038] 显示时，使用二维坐标系，X轴为高度，Y轴为污染物浓度。
[0039] 污染物从烟囱口排出后，上升一定高度后开始下降，因此，这里只考虑从地面开始 到H+AH，再高的地方不会有该烟囱排放的污染物出现。
[0040] 曲线图方式显示时，在上述高度范围内，如果没有设定高度间隔及范围，每隔 0. 5~3米取一点进行计算。
[0041] 根据计算机的计算速度和选取点的密集度，污染物浓度的显示更新周期为1-10 秒钟。
[0042] 国家电网对烟囱出口的污染物排放有明确的规定，如S02，脱硫前其浓度在2000 mg.nT3左右，排放前需要脱硫处理，排放浓度应小于200mg.nT3。
[0043] 污染物浓度显示颜色的选取可以按照一般的习惯，如污染物浓度〈30mg.nT3时显 不绿色，30mg.m3-50mg.m3时显不蓝色，50mg.m3-100mg.m3时显不黄色，100mg.m3-200 mg.nT3时显示橙色，大于200mg.nT3时显示红色，某种颜色的显示可以根据浓度从淡到浓；或 者可以按照灰度显示，污染物浓度0-500mg.nT3对应无色-黑色；污染物浓度越大，颜色越 深。
[0044] 当改变显示比例时，按照上述方式，重新进行计算并显示。
[0045] 将上述计算结果存储，可以在显示器上回放某一时间段污染物扩散的情况。
[0046] 根据该扩散模型，可以预测污染物的扩散情况。
[0047] 下面结合附图对本发明做进一步说明。
[0048] 参看图1，主机与显示器、GIS服务器、数据采集系统连接，数据采集系统连接CEMS系统，数据采集系统通过网络连接烟囱5内部设置的烟囱内部出口位置的温度传感器2、测 量室外温度的温度传感器1、设置在烟囡5顶端外侧的风向风速传感器3。所有传感器都有 编号。
[0049] 主机采用塔式扩散模型计算污染物的扩散，塔式扩散模型的数学表达式为： 塔式扩散模型的数学表达式为：
其中：x、y、z为空间坐标，H是烟囱口距地面的高度，巾为扩散因数，y为风速，t为 扩散时间，A为源强，a为颗粒属性系数，0为高度系数，空间坐标系是右手坐标系，原点 是烟囱在地面上的投影，x轴为主风向，y轴为横风向，z轴为垂直向。
[0050] 数据采集系统通过风向风速传感器3取得当前风向，待其稳定后选取主风向作为x轴，同时定义y和z轴。风速y为平均风速，初始值从气象数据中取得，经过一定时间后 通过实测的平均风速进行修正。
[0051] 烟囱内部出口位置的温度传感器2和测量室外温度的温度传感器1分别测量烟气 出口温度Ts和环境大气温度Ta，这两个参数在计算上升高度AH时使用。
[0052] 主机根据上述参数实时计算各成分的污染物浓度，并以不同颜色显示在电子地图 上。当在电子地图上选取一点时，使用柱状图或曲线图显示该点不同高度的污染物浓度，如 图2和图3所示。
[0053] 下面的表格是实验参数：
上表中，距离为主风向方向距烟囱口的距离，具体到公式中，x=距离，y=〇,z=2m;风速 为平均风速。 可以看出，采用塔式扩散模型得出的结果比采用普通高斯算法得出的结果更接近实测 值。作为显示的依据，最终在显示器上模拟显示的污染物扩散过程更接近实际情况。
【主权项】
1. 基于塔式扩散模型的火电厂污染物排放监测显示方法，基于火电厂污染物排放检测 系统实现，系统中包括主机，GIS服务器和数据采集系统，数据采集系统连接CEMS系统，其 特征在于：主机中存储当地全年的气象资料、每个烟囱的位置和高度信息；在每个烟囱内 部出口位置设置温度传感器，在室外设置温度传感器，在每个烟囡顶端外侧设置风向风速 传感器，所有传感器编号并连接到数据采集系统，主机采用塔式扩散模型计算污染物的扩 散，塔式扩散模型的数学表达式为：其中：x、y、z为空间坐标，H是烟囱口距地面的高度，Φ为扩散因数，μ为风速，t为扩 散时间，λ为源强，α为颗粒属性系数，Θ为高度系数，空间坐标系是右手坐标系，原点是 烟囱在地面上的投影，X轴为主风向，y轴为横风向，Z轴为垂直向； 主机计算各时间点不同位置的污染物浓度，并根据污染物浓度等级，以不同颜色显示 在电子地图上。2. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其 特征在于：所述颗粒属性系数α与污染物颗粒的直径相关：直径为O-HMffl颗粒属性系数 为0,直径为15-3〇Mm颗粒属性系数为0. 3,直径为31-47Mm颗粒属性系数为0. 5,直径为 48-75Mm颗粒属性系数为0. 8,直径大于76Mm颗粒属性系数为1 ;气体污染物按小于15um的 颗粒计算。3. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其特 征在于：所述高度系数Θ与高度z的关系如下：Η+ΔΗ >z彡Η，Θ为I ;H >z彡0. 75H，Θ 为 0· 75 ;0· 75H >z 彡 0· 5H，Θ 为 〇· 5 ;0· 5H >z 彡 0· 25H，Θ 为 〇· 25 ;0· 25H >z 彡 0，Θ 为 〇. I ;H是烟囱口距地面的高度，是排放物从烟囱口出来后上升的高度。4. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其特 征在于：所述扩散因数Φ与高度z的关系如下： Η+ΔΗ>ζ$Η，Φ*(3θ/2)μL0·87; H >ζ 彡 0· 75Η，Φ 为（Θ /2) yt0.84; 0· 75H >z 彡 0· 5H，Φ 为（Θ /2) μ t0.82; 0· 5H >z 彡 0· 25H，Φ 为（Θ /5) μ t0.8; 0· 25H >z 彡 0, Φ 为（Θ /6) yt0.78; H是烟囱口距地面的高度，是排放物从烟囱口出来后上升的高度，μ为风速，t为扩散 时间。5. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其特 征在于：数据采集系统通过风向风速传感器获得风速μ和主风向。6. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其特 征在于：由数据采集系统从CEMS获得排放污染物中的成分和各自的源强λ，在计算某一位 置的污染物浓度时，按照污染物中不同成分分别计算，分别显示。7. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其特 征在于：由数据采集系统从CEMS获得排放污染物中的成分和各自的源强λ，在计算某一位 置的污染物浓度时，按照污染物中不同成分分别计算，不同成分的污染物浓度加权求和为 综合浓度，显示综合浓度。8. 根据权利要求6或7所述的基于塔式扩散模型的火电厂污染物排放监测显示方法， 其特征在于：设定高度参数，显示的内容为设定高度的计算结果。9. 根据权利要求6或7所述的基于塔式扩散模型的火电厂污染物排放监测显示方法， 其特征在于：选取地图上一个点，使用柱状图或曲线图显示该点不同高度的污染物浓度。10. 根据权利要求6或7所述的基于塔式扩散模型的火电厂污染物排放监测显示方法， 其特征在于：将显示的区域等分成N*M个网格，每个网格的污染物浓度按照其中心点的污 染物浓度显示；或者在每个网格内选取3~9个点，分别计算选取点的污染物浓度，取其平均 值，每个网格的污染物浓度按照平均值显示。11. 根据权利要求10所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其 特征在于：N的取值范围是10~50, M的取值范围是10~30。12. 根据权利要求10所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其 特征在于：每个网格内选取5个点，位置分别为网格的中心及四角。
【专利摘要】基于塔式扩散模型的火电厂污染物排放监测显示方法，特别涉及污染物在大气中的排放监测，用于在电子地图上显示火电厂排放的烟气污染物在各区域的浓度。基于火电厂污染物排放检测系统实现，系统中包括主机，GIS服务器和数据采集系统，数据采集系统连接CEMS系统，采用塔式扩散模型计算污染物的扩散，塔式扩散模型的数学表达式为：主机计算各时间点不同位置的污染物浓度，并根据污染物浓度等级，以不同颜色显示在电子地图上。
【IPC分类】G01N33/00, G01N15/06
【公开号】CN104897853
【申请号】CN201510276062
【发明人】李春海, 陈洪雨, 陈贺, 芦斌, 袁琳琳, 刘晓龙, 崔振伟, 支智勇, 梁岩
【申请人】石家庄科林电气股份有限公司
【公开日】2015年9月9日
【申请日】2015年5月27日
【技术领域】
[0001] 本发明涉及监控领域，特别涉及污染物在大气中的排放监测，用于在电子地图上 显示火电厂排放的烟气污染物在各区域的浓度。
【背景技术】
[0002] 环境污染问题日益严重，雾霾天气不断威胁人类的健康。近年来，国家电网公司 积极落实相关政策要求，大力开展了污染物排放系统功能平台的建设工作，实现对固定污 染源排放的污染物进行连续地、实时地跟踪测定，连续测定颗粒物、气态污染物浓度和排放 率。
[0003] 现在的技术，只是通过表格或曲线、棒图显示当前污染物排放量，不能直观地反映 污染物排放对周边环境的影响。
[0004] 大气污染物扩散模式是模拟大气污染物的扩散、迀移过程，预测在不同污染源条 件、气象条件及下垫面条件下某污染物浓度时空分布的数学模型，是低层大气中污染物迀 移和扩散规律的、简单化的数学描述。根据不同的建模理论体系、污染物迀移、扩散过程以 及不同的描述对象，模式的形式也各不相同。
[0005] 受地形、气象、大气污染物的物理化学特征、污染源特征等多种因素的制约，不同 的扩散模式都有各自不同的考虑因素和适用范围，选择恰当的扩散模式能较为准确地模拟 污染物的扩散及分布。目前，已经有许多发展成熟的污染物扩散模式应用，国内的相关研宄 中尤其以高斯类模式的应用最为广泛。
[0006] 高斯扩散模型一种经验模型，大量实验数据表明高斯分布是污染物在分布的较好 的近似形式，高斯模型数学形式简单，通过大规模的实验获得模型计算的主要参数，可以在 各种气象条件下选用，高斯扩散模型有两个特点：浓度计算在垂直方向和水平方向都采用 高斯分布假设，湍流分类和扩散参数采用离散化的分类方法，它的优点是仅用常规气象观 测资料就可以评估扩散参数。
[0007] 高斯扩散模型基于以下条件： 坐标系：原点为地面源的排放点或高架源排放点在地面上的投影，右手坐标系，x轴为 主风向，y轴为横风向，z轴为垂直向。
[0008] 几点假设：1、污染物浓度在y、z方向为正态分布，2、全部高度风速均匀稳定，3、源 强连续均匀，4、扩散中不考虑转化，污染物是守恒的，5、x向上，平流作用远大于扩散作用， 6、地面足够平坦。
[0009] 无界空间连续点源扩散模式为： c(x,y,z) =q/ (2 31y〇y 〇 z)exp[- (y2/2 0 y2+z2/2 0 z2)] 高架连续点源扩散模式为： c(x,y,z,H) =q/(2 31y〇y 〇 z)exp(-y2/2 〇 y2) *{exp[-(z-H) 2/2 0 y2] +exp[-(z+H)2/20z2)]} 其中，y为平均风速，q为源强，〇y、〇z为侧向和竖向扩散参数，H为有效烟囱高度。
[0010] 具体到火电厂烟气污染物排放，由于火电厂的特殊性，如污染物在出口的温度较 高、污染物的组成比较复杂以及气压等因素，上述扩散模式预测此类污染物的扩散会产生 较大的偏差。
【发明内容】
[0011] 为了更加准确地预测火电厂烟气污染物的排放，实时在电子地图上显示污染物浓 度，提出了本发明。
[0012] 为了达到上述目的，本发明采用的技术方案是：基于塔式扩散模型的火电厂污染 物排放监测显示方法，基于火电厂污染物排放检测系统实现，系统中包括主机，GIS服务器 和数据采集系统，数据采集系统连接CEMS系统，关键在于：主机中存储当地全年的气象资 料、每个烟囡的位置和高度信息；在每个烟囡内部出口位置设置温度传感器，在室外设置温 度传感器，在每个烟囡顶端外侧设置风向风速传感器，所有传感器编号并连接到数据采集 系统，主机采用塔式扩散模型计算污染物的扩散，塔式扩散模型的数学表达式为：
其中：x、y、z为空间坐标，H是烟囱口距地面的高度，巾为扩散因数，y为风速，t为扩 散时间，A为源强，a为颗粒属性系数，0为高度系数，空间坐标系是右手坐标系，原点是 烟囱在地面上的投影，x轴为主风向，y轴为横风向，z轴为垂直向； 主机计算各时间点不同位置的污染物浓度，并根据污染物浓度等级，以不同颜色显示 在电子地图上。
[0013] 颗粒属性系数a如下表：
根据高度得到的高度系数9如下表：
扩散因数巾根据温度、气压、风速等多方面因素影响来确定，本算法按照高度系数0 来划分不同高度时的扩散因数：
1、塔式扩散模型的推导过程： 假定由污染源排放出的大气污染物全都由单位面元dydz中通过，根据质量守恒定律 和连续性原理，单位时间内通过单位面源的污染物就等于污染源单位时间内排放出的污染 物质量，故可用公式表示如下：
式中，污染源排放的污染物浓度，Q:污染源单位时间内排放的污染物质量。
[0014] 将高斯曲线的参数代入得到：
将其进行转换可得：
根据普哇松积分公式得到
将该值代回得
在考虑到高度以及风速的影响，最终可得到任一点的污染物浓度公式为：
小为扩散因数，U为风速，t为扩散时间，A为源强，a为颗粒属性系数，0为高度系 数。
[0015] 2、上升高度AH应用国家标准，定义如下： 烟囱的高度为实际测量值， 当QH彡 2100kW和（TS-Ta)彡 35K时， Qh=0. 35PaQvAT/TS,AT=TS-Ta AH=n0QHnl.Hsn2.y-1 QH为烟囱的热排放率，这是一个固定参数，存储在主机中，Qv是实际排烟率，由CEMS系 统中的烟气流量计量仪器获得，单位是m3/S。
[0016] y为烟囱出口处的平均风速，1；和Ta分别为烟气出口温度和环境大气温度，n^为 烟气热机地表状况系数，^为烟气热释放率指数，n2为烟囡高度指数，Pa为大气压力，大气 压力根据烟囱出口处的高度确定。下表为系数的选取：
CEMS是英文ContinuousEmissionMonitoringSystem的缩写，是指对大气污染源排 放的气态污染物和颗粒物进行浓度和排放总量连续监测并将信息实时传输到主管部门的 装置，被称为"烟气自动监控系统"，亦称"烟气排放连续监测系统"或"烟气在线监测系统"。 从该系统中，可以实时得到排放的污染物中的成分和源强。
[0017] 采用本发明，使用适合火电厂污染物排放的塔式扩散模型，能够更加准确地计算 火电厂烟气污染物的排放，根据计算结果，在电子地图上实时显示烟囱口周围的污染物浓 度，同时，还可以预测污染物的扩散，模拟显示污染物扩散过程。
【附图说明】
[0018] 图1为本发明的系统组成框图， 图2为柱状图显示污染物浓度的示意图， 图3为曲线图显示污染物浓度的示意图。
[0019] 其中，1为测量室外温度的温度传感器，2为烟囱内部出口位置的温度传感器，3为 风向风速传感器，5为烟囡。
【具体实施方式】
[0020] 基于塔式扩散模型的火电厂污染物排放监测显示方法，基于火电厂污染物排放检 测系统实现，系统中包括主机，GIS服务器和数据采集系统，数据采集系统连接CEMS系统。 GIS服务器用来存储和显示电子地图，数据采集系统采集算法需要的各种参数。
[0021] 塔式扩散模型的数学表达式为：
其中：x、y、z为空间坐标，H是烟囱口距地面的高度，巾为扩散因数，y为风速，t为扩 散时间，A为源强，a为颗粒属性系数，0为高度系数，空间坐标系是右手坐标系，原点是 烟囱在地面上的投影，X轴为主风向，y轴为横风向，z轴为垂直向。
[0022] 算法中，y的取值是平均风速。一般的，是按照统计规律，从气象资料中取得当前 日期的平均风速，主机中存储当地全年的气象资料就是为了取得上述数值；另外，主机还存 储每个烟囱的位置和高度信息；在每个烟囱内部出口位置设置温度传感器测量烟囱出口温 度，在室外设置温度传感器测量室外温度，在每个烟囱顶端外侧设置风向风速传感器测量 实时风向风速。
[0023] 火电厂烟囱排放的污染物是气体和颗粒的混合物，如气体的二氧化硫、氮氧化物， 固体颗粒的烟尘等。CEMS系统可以获得排放物的种类和源强。
[0024] CEMS系统主要测量S02、NOx、02、粉尘的排放强度，也就是源强。这里的粉尘大小 是取颗粒的平均尺寸。
[0025] 主机获得实时参数，根据塔式扩散模型，计算各时间点不同位置的污染物浓度，并 根据污染物浓度等级，以不同颜色显示在电子地图上。
[0026] 为了区分不同颗粒大小的污染物，算法中使用了颗粒属性系数a。颗粒属性系数 a与污染物颗粒的直径相关：直径为0-14Mm颗粒属性系数为0,直径为15_3〇Mm颗粒属性 系数为0. 3,直径为31-47Mm颗粒属性系数为0. 5,直径为48-75Mm颗粒属性系数为0. 8,直 径大于76Mm颗粒属性系数为1 ;气体污染物按小于15um的颗粒计算。在现有技术基础上， 只计算平均尺寸的颗粒污染物的排放，也就是说按照从CEMS系统中得到的数据进行计算。
[0027] 本发明中，还充分考虑了高度对污染物扩散的影响，增加了高度系数0和与 高度相关的扩散因数:高度系数9与高度z的关系如下：H+AH>z多H，0为1;H >z彡 0? 75H，9 为 0? 75 ;0? 75H>z彡 0? 5H，9 为 0? 5 ;0? 5H>z彡 0? 25H，9 为 0? 25 ;0? 25H >z彡 0，9 为 0? 1〇
[0028] 扩散因数巾与高度z的关系如下： H+AH>z>H，<i)S(3 0/2)yto-87; H>z彡 0? 75H，<i> 为（0 /2)yt0.84; 0? 75H>z彡 0? 5H，<i> 为（9 /2)yt0.82; 0? 5H>z彡 0? 25H，<i> 为（9 /5)yt0.8; 0? 25H>z彡 0, <i> 为（9 /6)yt0.78。
[0029] 上述公式中，H是烟囱口距地面的高度，AH是排放物从烟囱口出来后上升的高度， y为风速，t为扩散时间。
[0030] 风速y为平均风速，是污染物在扩散过程中的平均值。风向和平均风速一般从气 象资料中获得。为了使结果更精确，在具体计算时，数据采集系统通过风向风速传感器获得 风速y和主风向，开始计算时使用获得的主风向作为x轴，气象资料中的平均风速作为y。 连续取得实时风速并计算平均值，经过一段时间后，使用实时取得的平均风速作为y，并在 其后的计算过程中不断修正。
[0031] 火电厂烟囱排放的污染物是气体和颗粒的混合物，如气体的二氧化硫、氮氧化物， 固体颗粒的烟尘等。在计算各时空点的污染物浓度时，不能只计算单一的污染物浓度，要按 照不同尺寸的污染物分别计算。由数据采集系统连 接CEMS系统获得排放污染物中不同尺 寸的污染物和各自的源强A，将探测到的参数代入扩散模型分别计算。目前只是计算S02、 NOx和平均粉尘尺寸的污染物浓度，当检测技术可以实时得到排放的污染物中其他气体和 粉尘中各种尺寸污染物的源强时，本算法会得到更精确的结果。
[0032] 将污染物浓度显示到有电子地图的屏幕上时，使用离散点的浓度值来代表某一区 域的浓度。为了在计算量和显示效果之间取得平衡，将显示的区域等分成N*M个网格，每个 网格的污染物浓度按照其中心点的污染物浓度显示，或者在每个网格内选取3~9个点，分 别计算选取点的污染物浓度，取其平均值，每个网格的污染物浓度按照平均值显示，或者每 个网格内选取5个点，位置分别为网格的中心及四角。N的取值范围是10~50,M的取值范 围是10~30。根据实际情况，N和M的取值可以更大。
[0033] 由于污染物由多种成分组成，因此，由数据采集系统从CEMS获得排放污染物中的 成分和各自的源强A，在计算某一位置的污染物浓度时，按照污染物中不同成分分别计算。 [0034] 显示时，可以分别显示单独的污染物浓度，也可以将不同成分的污染物浓度加权 求和为综合浓度，显示综合浓度。权重的比例分配可以按需要设定，如302分配权重50%、N0x30%，粉尘20%。加权求和的过程为各成分的污染物浓度乘以权重后相加。
[0035] 为了精确显示污染物浓度，显示时要求有高度参数。设定高度参数，按照高度参数 进行计算，并显示设定高度的计算结果。这里实现了分时分层显示，如显示地表的污染物分 布情况、显示与人体等高的1. 7米处污染物分布情况、显示离地高度5米的污染物分布情况 等。
[0036] 另外，还可以不设定高度参数，系统根据风速风向等因素，自动计算各位置污染物 浓度最大的值并显示。如在烟囱出口附近，与烟囱等高处污染物浓度最大，x轴方向20米 位置，低于烟囱出口 2米处污染物浓度最大，x轴方向200米位置，低于烟囱出口 4米处污 染物浓度最大，选取上述最大值并显示。
[0037] 在显示一定范围内污染物分布的同时，可以在地图上选取一个点，计算该点不同 高度的各种污染物浓度，使用柱状图或曲线图显示，高度间隔及高度范围人为设定。
[0038] 显示时，使用二维坐标系，X轴为高度，Y轴为污染物浓度。
[0039] 污染物从烟囱口排出后，上升一定高度后开始下降，因此，这里只考虑从地面开始 到H+AH，再高的地方不会有该烟囱排放的污染物出现。
[0040] 曲线图方式显示时，在上述高度范围内，如果没有设定高度间隔及范围，每隔 0. 5~3米取一点进行计算。
[0041] 根据计算机的计算速度和选取点的密集度，污染物浓度的显示更新周期为1-10 秒钟。
[0042] 国家电网对烟囱出口的污染物排放有明确的规定，如S02，脱硫前其浓度在2000 mg.nT3左右，排放前需要脱硫处理，排放浓度应小于200mg.nT3。
[0043] 污染物浓度显示颜色的选取可以按照一般的习惯，如污染物浓度〈30mg.nT3时显 不绿色，30mg.m3-50mg.m3时显不蓝色，50mg.m3-100mg.m3时显不黄色，100mg.m3-200 mg.nT3时显示橙色，大于200mg.nT3时显示红色，某种颜色的显示可以根据浓度从淡到浓；或 者可以按照灰度显示，污染物浓度0-500mg.nT3对应无色-黑色；污染物浓度越大，颜色越 深。
[0044] 当改变显示比例时，按照上述方式，重新进行计算并显示。
[0045] 将上述计算结果存储，可以在显示器上回放某一时间段污染物扩散的情况。
[0046] 根据该扩散模型，可以预测污染物的扩散情况。
[0047] 下面结合附图对本发明做进一步说明。
[0048] 参看图1，主机与显示器、GIS服务器、数据采集系统连接，数据采集系统连接CEMS系统，数据采集系统通过网络连接烟囱5内部设置的烟囱内部出口位置的温度传感器2、测 量室外温度的温度传感器1、设置在烟囡5顶端外侧的风向风速传感器3。所有传感器都有 编号。
[0049] 主机采用塔式扩散模型计算污染物的扩散，塔式扩散模型的数学表达式为： 塔式扩散模型的数学表达式为：
其中：x、y、z为空间坐标，H是烟囱口距地面的高度，巾为扩散因数，y为风速，t为 扩散时间，A为源强，a为颗粒属性系数，0为高度系数，空间坐标系是右手坐标系，原点 是烟囱在地面上的投影，x轴为主风向，y轴为横风向，z轴为垂直向。
[0050] 数据采集系统通过风向风速传感器3取得当前风向，待其稳定后选取主风向作为x轴，同时定义y和z轴。风速y为平均风速，初始值从气象数据中取得，经过一定时间后 通过实测的平均风速进行修正。
[0051] 烟囱内部出口位置的温度传感器2和测量室外温度的温度传感器1分别测量烟气 出口温度Ts和环境大气温度Ta，这两个参数在计算上升高度AH时使用。
[0052] 主机根据上述参数实时计算各成分的污染物浓度，并以不同颜色显示在电子地图 上。当在电子地图上选取一点时，使用柱状图或曲线图显示该点不同高度的污染物浓度，如 图2和图3所示。
[0053] 下面的表格是实验参数：
上表中，距离为主风向方向距烟囱口的距离，具体到公式中，x=距离，y=〇,z=2m;风速 为平均风速。 可以看出，采用塔式扩散模型得出的结果比采用普通高斯算法得出的结果更接近实测 值。作为显示的依据，最终在显示器上模拟显示的污染物扩散过程更接近实际情况。
【主权项】
1. 基于塔式扩散模型的火电厂污染物排放监测显示方法，基于火电厂污染物排放检测 系统实现，系统中包括主机，GIS服务器和数据采集系统，数据采集系统连接CEMS系统，其 特征在于：主机中存储当地全年的气象资料、每个烟囱的位置和高度信息；在每个烟囱内 部出口位置设置温度传感器，在室外设置温度传感器，在每个烟囡顶端外侧设置风向风速 传感器，所有传感器编号并连接到数据采集系统，主机采用塔式扩散模型计算污染物的扩 散，塔式扩散模型的数学表达式为：其中：x、y、z为空间坐标，H是烟囱口距地面的高度，Φ为扩散因数，μ为风速，t为扩 散时间，λ为源强，α为颗粒属性系数，Θ为高度系数，空间坐标系是右手坐标系，原点是 烟囱在地面上的投影，X轴为主风向，y轴为横风向，Z轴为垂直向； 主机计算各时间点不同位置的污染物浓度，并根据污染物浓度等级，以不同颜色显示 在电子地图上。2. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其 特征在于：所述颗粒属性系数α与污染物颗粒的直径相关：直径为O-HMffl颗粒属性系数 为0,直径为15-3〇Mm颗粒属性系数为0. 3,直径为31-47Mm颗粒属性系数为0. 5,直径为 48-75Mm颗粒属性系数为0. 8,直径大于76Mm颗粒属性系数为1 ;气体污染物按小于15um的 颗粒计算。3. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其特 征在于：所述高度系数Θ与高度z的关系如下：Η+ΔΗ >z彡Η，Θ为I ;H >z彡0. 75H，Θ 为 0· 75 ;0· 75H >z 彡 0· 5H，Θ 为 〇· 5 ;0· 5H >z 彡 0· 25H，Θ 为 〇· 25 ;0· 25H >z 彡 0，Θ 为 〇. I ;H是烟囱口距地面的高度，是排放物从烟囱口出来后上升的高度。4. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其特 征在于：所述扩散因数Φ与高度z的关系如下： Η+ΔΗ>ζ$Η，Φ*(3θ/2)μL0·87; H >ζ 彡 0· 75Η，Φ 为（Θ /2) yt0.84; 0· 75H >z 彡 0· 5H，Φ 为（Θ /2) μ t0.82; 0· 5H >z 彡 0· 25H，Φ 为（Θ /5) μ t0.8; 0· 25H >z 彡 0, Φ 为（Θ /6) yt0.78; H是烟囱口距地面的高度，是排放物从烟囱口出来后上升的高度，μ为风速，t为扩散 时间。5. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其特 征在于：数据采集系统通过风向风速传感器获得风速μ和主风向。6. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其特 征在于：由数据采集系统从CEMS获得排放污染物中的成分和各自的源强λ，在计算某一位 置的污染物浓度时，按照污染物中不同成分分别计算，分别显示。7. 根据权利要求1所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其特 征在于：由数据采集系统从CEMS获得排放污染物中的成分和各自的源强λ，在计算某一位 置的污染物浓度时，按照污染物中不同成分分别计算，不同成分的污染物浓度加权求和为 综合浓度，显示综合浓度。8. 根据权利要求6或7所述的基于塔式扩散模型的火电厂污染物排放监测显示方法， 其特征在于：设定高度参数，显示的内容为设定高度的计算结果。9. 根据权利要求6或7所述的基于塔式扩散模型的火电厂污染物排放监测显示方法， 其特征在于：选取地图上一个点，使用柱状图或曲线图显示该点不同高度的污染物浓度。10. 根据权利要求6或7所述的基于塔式扩散模型的火电厂污染物排放监测显示方法， 其特征在于：将显示的区域等分成N*M个网格，每个网格的污染物浓度按照其中心点的污 染物浓度显示；或者在每个网格内选取3~9个点，分别计算选取点的污染物浓度，取其平均 值，每个网格的污染物浓度按照平均值显示。11. 根据权利要求10所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其 特征在于：N的取值范围是10~50, M的取值范围是10~30。12. 根据权利要求10所述的基于塔式扩散模型的火电厂污染物排放监测显示方法，其 特征在于：每个网格内选取5个点，位置分别为网格的中心及四角。
【专利摘要】基于塔式扩散模型的火电厂污染物排放监测显示方法，特别涉及污染物在大气中的排放监测，用于在电子地图上显示火电厂排放的烟气污染物在各区域的浓度。基于火电厂污染物排放检测系统实现，系统中包括主机，GIS服务器和数据采集系统，数据采集系统连接CEMS系统，采用塔式扩散模型计算污染物的扩散，塔式扩散模型的数学表达式为：主机计算各时间点不同位置的污染物浓度，并根据污染物浓度等级，以不同颜色显示在电子地图上。
【IPC分类】G01N33/00, G01N15/06
【公开号】CN104897853
【申请号】CN201510276062
【发明人】李春海, 陈洪雨, 陈贺, 芦斌, 袁琳琳, 刘晓龙, 崔振伟, 支智勇, 梁岩
【申请人】石家庄科林电气股份有限公司
【公开日】2015年9月9日
【申请日】2015年5月27日
最新回复(0)