具体实施方式

下面结合具体实施例对本发明进一步进行描述。

实施例1

一种基于高斯模式的大气污染溯源方法，其特征在于：包括如下处理步骤：(1)污染源信息处理步骤：通过污染源信息处理模块将研究区域内各污染源进行统一标识，赋值污染源ID；本模块需要对污染源的属性进行标识，包括污染排放类型、主要污染物类型排放信息；根据研究区域的污染源实际调研结果，制作形成标准分辨率的GRID格式文件，并导入数据库；(2)气象数据预处理：利用气象数据预处理模块对气象要素进行预处理，形成标准分辨率的GRID文件；(3)污染源实时数据处理：利用污染源实时数据处理模块通过在线监测系统实时获取污染源排放信息，通过标准化的处理流程，形成标准分辨率的GRID文件；(4)二段式溯源计算：分为模糊溯源和精准溯源，能够实现不同污染场景的追踪溯源计算。

所述的污染排放类型包括点源、线源、面源、体源，所述的主要污染物类型包括VOCs、PM₁₀、SO₂、氮氧化物，所述的排放信息包括管径、排口高度、排放效率、气体温度。

所述的污染源信息处理模块中调研的所有污染源集合记为{PS_i}，PS_i表示第i个污染源，每个污染源带有类型、位置、主要污染物、排口等众多属性。

步骤(2)中通过克里金插值法将相关要素进行空间插值处理，形成标准的GRID数据文件。

所述的气象要素中的风速(u)、风向(dir)、云量(CA)、辐射(rad)、温度(T)、湿度(hum)对污染源的追溯有重要影响，本发明中通过克里金插值法将相关要素进行空间插值处理，形成标准的GRID数据文件，避免大量的实时运算，提高模型计算效率，降低溯源时间成本，能够实现污染源分钟级溯源追踪，提升对行政管理的支撑水平。

其中σ_y，σ_z与风速(u)、风向(dir)、云量(CA)、辐射(rad)、温度(T)、湿度(hum)存在关联关系，需要根据不同的条件进行选择判断，在实际应用时，需进行优化。

步骤(3)中所述的污染源排放信息包括污染源ID、排放位置、类型、排放浓度。

各个厂区都需安装污染源排放在线监测系统，需要监测的要素包括排放速率(u_ps)、排放效率(q_ps)、污染物浓度({C_ps})、排放温度(T_ps)，用于计算不同排口的抬升高度(ΔH)和单位时间污染物排放量q_ps。

ΔH＝φ(u,u_ps,T_ps)

结合线下调研成果，对当前污染场景进行污染物类型识别，作为模块输入(污染物类型和浓度值等)。本模块根据输入条件计算出所有超标污染物，需针对每一种超标污染物进行溯源分析。针对某一种污染物，模糊溯源是利用当前污染物类型、风速风向和预处理完成的污染物类型分布等信息初步判断污染来源，定性分析出哪些污染源对污染场景负责，将该污染源定义为模糊源(Vague Source，VS)。

所述模糊溯源具体包括如下步骤：①确定输入的污染物类型和浓度集合{Type_i,C_i}；②计算各个类型污染物的超标程度，对于超标的污染物提取为新的集合{Over_i}；③根据风速(u)、风向(dir)以及污染源{PS_i}属性信息(污染物类型、位置)等信息，获取对当前污染场景可能造成贡献的模糊源{VS_i}

{VS_i}＝Θ(u,dir,{Over_i},{PS_i})

Θ为判断函数，判断PS_i是否为模糊源，其中

所述的精准溯源是在模糊溯源{VS_i}基础上，利用大气扩散理论高斯模式公式计算各个模糊源单独排放造成的空间污染分布特征，然后实现多模糊源在空间上的矢量叠加计算，接着利用当前污染场景的浓度信息，配合全局优化算法逆向分析不同模糊源的贡献系数；当前场景有多种超标污染物时，单独分析各类型污染物的模糊源的贡献系数，然后按照权重系数进行综合计算，分析出不同污染源对本次污染场景的综合贡献系数。

所述的精准溯源的步骤如下：①针对模糊源集合{VS_i}，假设第i个模糊源VS_i的排放效率为Q_VS,i，根据排放类型和相关基础信息计算出在监测点位置的污染物浓度为C_VS,i，计算公式如下

C_vs,i＝F(Q_vs,i,VS_i,u)

F为污染物浓度模拟函数，计算特定环境下污染物的大气扩散过程。

②根据{VS_i}、{Q_vs,i}、{C_vs,i}以及污染物浓度C_i，将每个模糊源的不同组合参数进行矢量叠加计算，通过优化算法分析出最可靠的模糊源排放组合{Q_vs,i}_s，使得下面的公式成立

Abs，绝对值函数；ε为置信值，一般取1％。

③对于某一种污染物按照①②步骤可以求出一个最优的{C_vs,i}，可以计算出{PS_i}在第种污染物的贡献值{R_i}

其中不在模糊源{VS_i}中的污染源的R_i＝0。

④对于多种污染物{Over_i}的场景，针对不同的Over_i重复①②③步骤，可以得到第i污染物的企业贡献率{R_j,i}，R_j,i表示第j个企业污染源对第i种超标污染物的贡献值。

⑤根据不同类型污染物的重要程度，按照专家经验分配不同污染物的计算权重{ω_i}，计算得到当前污染场景下各企业污染源的综合贡献系数

R_j表示第j个企业污染源综合贡献系数。

应用例1

某中部工业城市有重工业园区，园区内有水泥厂、钢铁厂、发电厂等多种类型的企业，园区主要监测的污染物有悬浮颗粒物(TSP)、二氧化硫(SO2)、氮氧化物(NOX)、挥发性有机物(VOCS)等。

(1)园区污染物信息收集

园区内收集12家企业单位基础资料，每一个家企业都会产生不同类型的污染物，以园区左下角为坐标原点，正东方向为X轴正方向，正北方向为Y轴正方向建立平面坐标系，详细信息统计如下：

表1

(2)气象信息收集

夏季某日14时，园区天气晴，无云，10米高测得风速2.3m/s；风向正西南；1.5米测得气温为28℃，大气稳定度为较稳定。

(3)实时信息获取

某日14时，在线监测系统显示某固定监测点(1000,800)实时监测到各项指标的值摘录如下：

表2

某日14时，12家企业的排气口污染物类型统计如下：

表3

(4)二段式溯源分析

根据实测的污染物指标值，TSP、SO₂存在超标现象，针对两种污染物超标场景进行溯源分析。

1)以TSP超标溯源分析为例

①企业的排污类型刷选出可能污染源编号集合，除了QY005都包含在内

{QY001,QY002,QY003,QY004,QY006,QY007,QY008,QY009,QY0010,QY0011,QY0012}

②风速、风向以及监测点和污染源(企业排口)的相对位置进行进一步筛选，认为处于监测点风向下游的企业对本次污染场景不造成影响，概化为风向向量与某污染物到监测点的位移向量的内积小于0，即

以企业QY001为例，风向向量为位移向量计算两向量的内积故认为QY001对本次场景中TSP超标有影响，即QY001为模糊源。

根据上述原理可以得到本次污染物超标场景中，TSP超标所有的模糊源集合为：{VS_i}＝{QY001,QY002,QY003,QY004,QY006,QY008}；

根据不同污染源气体中TSP的监测数据统计，模糊源排放效率摘录如下表：

表4

③计算污染源QY001在监测点位(1000,800)处形成的污染物浓度C_QY001，假设监测点高度为Z＝1.8m；带入公式(4)计算，得到浓度值C_QY001＝1600，单位μg/m³。

同理，计算出所有模糊源根据公式(6)判断是否满足条件。计算可知：

不满足条件则需调整模型计算参数，重新计算所有模糊源的扩散浓度；经过多次计算得到所有模糊源的优化后的计算浓度值：最终获取{VS_i}＝{QY001,QY002,QY003,QY004,QY006,QY008}的合理计算浓度为计算置信值：满足条件。

④根据公式(7)计算所有污染源对TSP超标的贡献值依次为

{R_TSP,j}＝{0.216,0.278,0.33,0.172,0,0.004,0,0,0,0,0,0}

⑤按照同样的原理，计算所有污染源中对SO₂超标的贡献值：

⑥依据专家经验判断，两种指标的权重分布为ω_TSP＝0.4，计算得到所有污染源对本次污染场景的综合贡献值

上述实施案例仅为本发明中部分实施案例，但本发明的实施方式并不受上述实施案例的限制，如实施例中方案的各种形式的组合，其他的任何未背离本发明的精神实质与原理下所做的改变、修饰、替代、组合均应为等效的替换方式，都在本发明的保护范围之内。