合空气污染确定出个监测点的综合污染等级。用层次分析法求出空气污染、水污染和土壤污染的权重，经一致性检验后，对各等级进行评分，根据得分标准得出个监测点的污染程度，根据污染程度对居民进行相应的经济赔偿。
5.2 问题一具体分析及数据的处理
5.2.1 垃圾焚烧厂附近监测点的确定
根据附件中所给的谷歌地图经纬度坐标具体确定到该垃圾焚烧厂在深圳市宝安区白鸽湖路67号，如图一，由图可知，居住点多分布在北侧和东侧，其南侧为平湖生态园对污染气体有净化作用，其中有一些居民点离该焚烧厂较远，因此我们选取离污染源较近的居民点作为监测点，对其污染程度进行分析，所以我们选取如下点作为监测点：
点A：新田社区 （西北/1000m）
点B：平湖日航塑胶机械商行 （东北/653m）
点C：白虎头小区 （东/1500m）
点D：华侨山庄 （北/2000m）
括号中分别表示各监测点相对于垃圾焚烧厂所处的方位和离垃圾焚烧厂的距离。
图一 监测点分布图
4
5.2.2 基于高斯（烟羽）气体扩散模型的污染评价模型
对于本问题，影响气体污染的因素较多，因此我们考虑高斯气体扩散模型
?1?，该
模型能有效的将焚烧炉污染物排放量、居住点离开垃圾焚烧厂的距离、风力和建筑物的遮挡程度等因素综合的考虑进来，又因为降雨对烟羽中的颗粒物及气溶胶具有清洗作用，因此我们结合高斯烟羽气体扩散模型
?2?对各监测点污染浓度进行最终的求解。
在高架点源的实际扩散中，污染物可能受到地面障碍物的阻挡，因此应当考虑地面对扩散的影响。处理的方法是，或者假定污染物在扩散过程中的质量不变，到达地面时不发生沉降或化学反应而全部反射；或者污染物在没有反射而被全部吸收，实际情况应在这两者之间。
点源在地面上的投影点O作为坐标原点，有效源位于z轴上某点， z＝H。高架有效源的高度由两部分组成，即H＝h＋Δh，其中h为排放口的有效高度，Δh是热烟流的浮升力和烟气以一定速度竖直离开排放口的冲力使烟流抬升的一个附加高度，如图二。
图二 地面全反射的高架连续点源扩散 其中
5
vsD?Ts?Ta???h?1.5?2.7D? ?? (5-1) k?Ts?式中
Vs—烟流出口速度，m/s； D—烟囱出口内径，m；
k—烟囱出口的环境平均风速，m/s；
Ta—环境平均气温度，K。
Ts—烟气出口温度，K；
其中K为开氏温度单位。
在以排放口为原点的坐标系中, 点源下风向任一点( x,y, z )的浓度分布函数为:
（5－2）
当污染物到达地面后被全部反射时，可以按照全反射原理，用“像源法”来求解空间某点k的浓度。图二中k点的浓度显然比大空间点源扩散公式（5－2）计算值大，它是位于(0，0，H)的实源在k点扩散的浓度和反射回来的浓度的叠加。反射浓度可视为由一与实源对称的位于(0，0，－H)的像源（假想源）扩散到k点的浓度。
?1?y2z2??qC?x,y,z??exp???2?2?????2?u?y?z2?yz?????? 由图可见，k点在以实源为原点的坐标系中的垂直坐标为(z-H)，则实源在k点扩散
的浓度为式（5－2）的坐标沿z轴向下平移距离H：
22????z?H??q?1y?Cs?exp???2???22?u?y?z2??z????y???? （5－3）
k点在以像源为原点的坐标系中的垂直坐标为(z＋H)，则像源在k点扩散的浓度为
式（5－2）的坐标沿z轴向上平移距离H：
22????z?H??q1y??Cx?exp???2???22?u?y?z2??z????y???? （5－4）
由此，实源Cs与像源Cx之和即为k点的实际污染物浓度：
222????????z?H???z?H???q?y??C?x,y,z,H??exp?2??exp???exp???22??2?u?y?z?????2?y???2?z??2?z????
6
式中
C— 空间点( x, y, z )的污染物浓度, m g /m3;
q— 源强, 单位时间污染物排放量, m g / s;
u— 各监测点的平均风速, m / s.
?y、?z为烟气的扩散系数, 与大气稳定度和水平距离x有关, 并随x 的增大而增加。
我国GB3840－91《制定地方大气污染物排放标准的技术方法》定扩散参数?y、?z：
?3?采用如下经验公式确
?y??1x?及?z??2x?12 （5－5）
式中，γ1、α1、γ2及α2称为扩散系数。这些系数由实验确定，在一个相当长的x距离内为常数，可从GB3840－91的表中查取，由GB3840－91中可知深圳市的大气稳定度为B-C级，由此可得扩散系数的值如表1。
表1 稳定度 ?1 ?1 下风距离/m 0-1000 >1000 下风距离/m 0-500 >500 B-C 0.919325 0.875086 ? 20.229500 0.314238 ?2 0.114682 0.0757182 0.941015 1.0077
降雨对烟羽中的颗粒物及气溶胶具有清洗作用，可溶性气体与蒸汽亦可溶于降雨中，降雨过程造成的这类沉积是导致放射性气溶胶和气体向地面沉积的另一重要机制，因此我们结合高斯烟羽气体扩散模型对污染物浓度进行求解。通常以冲洗系数?描述降雨对烟羽中污染物清洗作用的大小。?与雨强的关系表达式为： ?其中：I为雨强（mm/h); a,b为经验系数。
式中，按放射性物质含碘，不含碘情况取值，对含碘的放射性物质，取a?8?10?5,b?0.6;对于不含碘情况，取a?1.2?10?5,b?0.5；对于湿沉积导致的烟羽消耗，采用湿沉积消耗因子对源强q进行修正：
?aIb
7
?x) q1(x)?qexp(?k式中各参数的取值见表2、表3.
表2 项目 Vs（m/s） D（m） k（m/s） Ts（K） Ta（K) I（mm/h)
A 表3 B C D 数值 15 2 2.7 365.5 296.15 119 ?y 131.4491221 88.83881363 189.0652431 243.1889495 79.854663 51.97425753 120.156548 160.5640107 2.6 1.8 1.5 1.9 ?z u?m/s? ?x,y,z? (1000,0,15) (653,0,15) (1500,0,15) (2000,0,18) 表2中，k根据附件4中所给一年风速平均值求得；表3中z为居民楼按每层3米计算
求得的高度。
将各参数值代入式（5-1）中解得?h?28。
将各数据代入入高斯（烟羽）公式用Matlab（代码见附件一）得到各监测点的气体污染浓度值，见下表4。
表4 SO2 NOx Hcl A 0.045567038 0.159484631 0.027340223 B 0.033455628 0.117094700 0.020073377 0.010036688 C 0.034430693 0.120507424 0.020658415 0.010329208 D 0.023900456 0.083651595 0.014340273 0.007170136 颗粒物 0.013670111
用Excel表格将各监测点的污染物浓度进行绘图，绘制的柱状图如下：
8