大气污染预报问题

大气污染预报问题

摘要大气污染给人类生产生活带来了极大的危害，因此大气污染的监测和预报为人们控制污染、减少排放、保护环境提供了重要的指导意义。本文通过模糊评估和灰色预测模型，建立起对大气污染的评价和预测体系，通过构建多元线性回归方程找出空气质量与气象参数之间关系，并由此对大气污染控制特提出建议。问题一中，找出SO2、NO2、PM10之间关系，即是对三者的污染程度排序，因此我们参考了大气日污染指数API的计算方法，确定每日的首要污染物，进行统计，发现对于每个城市，PM10为最主要的污染物，SO2次之，NO2再次之；为了对这六个城市排名，我们建立模糊综合评价模型，通过确定评价权重、构造模糊综合评价矩阵、选取模糊评价算子，得出各地空气质量的综合评价结果，排序后得到各地空气质量排名，从好到差依次为：E、C、A、B、D、F。问题二中，考虑到年份以及月份的影响，我们选取离预测点较近的一段时间作为预测的原始数据，同时考虑到数据真实性与合理性，建立非等间距序列的灰色预测模型，利用MATLAB软件编程求解辨识参数a，u，得到预测方程，对选取的预测原始数据与利用预测方程得到的预测值进行误差检验，认为除个别点外，模型可以接受，从而利用模型预测了2010年9月15日到21日的各地各污染物浓度以及相应的气象参数。问题三中，由于数据过多，为简化问题，我们取2003年至2010年8年各污染物的每月月平均值作为因变量，建立多元线性回归方程，找出与气压、温度、湿度、风速等自变量之间关系，并进行了显着性和复相关性的检验，最终得到六地空气质量与气象参数之间关系。问题四中，我们根据上述三题得到结果，对控制空气质量提出了针对性建议。关键字：模糊综合评价非等间距序列的灰色预测多元线性回归大气污染预测大气污染预报问题一、问题背景随着地球上人口的急剧增加，人类经济的快速增长，地球上的大气污染日趋严重，其影响也日趋深刻。颗粒状污染物和有害气体的大量排放，不仅时刻威胁着人们的身体健康，还给工业生产带来严重的危害，影响着经济的可持续性发展。同时，局部地区大气的污染会对全球性的气候变化产生影响。因此，加强大气质量的监测和预报对保护人类健康、改善生态环境是非常必要。二、问题重述大气污染给人类生产生活带来了极大的危害，因此大气污染的监测和预报为人们控制污染、减少排放、保护环境提供了重要的指导意义。如何通过监测数据来评判城市的空气质量、预测未来几天内大气污染情况、对空气质量的控制提出建议，是亟待解决的问题。具体如下：1）找出各个城市SO2、NO2、PM10之间的特点，并将几个城市的空气质量进行排序；2）对未来一周即2010年9月15日至9月21日各个城市的SO2、NO2、PM10以及各气象参数作出预测；3）分析空气质量与气象参数之间的关系；4）对控制空气质量提出相关建议。三、基本假设1）认为数据都是真实可靠，反应实际情况的，都是在同一条件下测量的，具有可比性；2）认为A、B、C、D、E、F、G六地间的污染情况不会互相影响；

3）认为SO2、NO2、PM10三者之间不会发生反应，互不影响。四、符号说明符号意义API空气污染指数R模糊综合评价矩阵H级别变量的特征值a，u辨识参数P大气压强的数值T温度的数值f空气湿度的数值v风速的数值FF检验统计量R2预测模型的复相关系数五、模型的建立与求解问题一：.问题分析找出SO2、NO2、PM10之间关系，即是对三者的污染程度排序，找出最主要的污染物。每日的首要污染物可由大气日污染指数API来判断，对每日首要污染物进行统计可以找出每个城市的主要污染物。SO2、NO2、PM10这三个因素均影响着某地的空气质量，必须将这三个方面都纳入评价体系，通过一定方法得到对空气质量的总体评价。.各个城市SO2、NO2、PM10之间特点1）空气污染指数API的计算及确定污染物的分指数Ii，可由实测的浓度值Ci按照分段线性方程计算。对于第i种污染物质第j个转折点的分指数值和相应的浓度值（Iij，C），可由表1确定。表1空气污染指数对应的污染物浓度限值污染指数API污染物浓度（mg/L）SO2NO2PM10500.0500.0800.0501000.1500.1200.1502000.8000.2800.3503001.6000.5650.4204002.1000.7500.5005002.6200.9400.600当第i种污染物浓度CijSCiSCi,j+1时，其分指数Ii：Ii—第i种污染物的污染分指数；Ci—第i种污染物的浓度监测值；I—第i种污染物第j个转折点的污染分项指数值；Iij+1—第i种污染物第j+1个转折点的污染分项指数值；C—第j个转折点上第i种污染物（对应于Iij）浓度限值；Cij—第j+1个转折点上第i种污染物（对应于Iij+1）浓度限值。污染指数的计算结果只保留整数,小数点后的数值全部进位。多种污染物的污染分指数都计算出来后，取最大者为该城市的空气日污染指数API，则该项污染物即为该城市当天空气中的首要污染物。式中:Ii为第i种污染物的分指数，n为污染物的项目数。2）各城市的空气污染指数及首要污染物采用各2003—2010年各地污染物数据，根据空气污染指数的计算方法，应用MATLAB软件算出每个城市每天的空气污染指数，并确定每日的首要污染物。对各地的每日首要污染物进行统计分析，得表2表2各城市三种首要污染物所占天数及百分比城市监测天数各首要污染物的天数及百分比SO2百分比%NO2百分比%PM10百分比%A108418717.2530.2889482.47B108419918.3600.0088581.64C1084595.4440.37102194.19D97213914.3010.1083285.60E608203.2910.1658796.55F1002727.0000.007373.003）判断SO2、NO2、PM10之间关系由表2可以看出，对A、B、C、D、E、F六个城市来说，首要污染物是PM10即可吸入颗粒物的天数占总监测天数的百分比是最大的，其次是是二氧化硫SO2,再次是氮氧化合物NO2。因此，对各城市而言，PM10是主要的污染物，SO2其次，NO2最少。.城市的空气质量排名1）模糊评价模型的建立a）确定评价因素集根据各地的空气质量监测的观测资料，选取二氧化硫SO2，二氧化氮NO2,总悬浮颗粒PM10三项污染物作为评价因素，即有评价因素集合：U={u1（SO2）,u2（NO2）,u3（PM10）} （1）b）确定评语集根据国家制定的空气环境质量评价标准,结合各地的空气监测特点和评价目的,取评语集为：V={v1(一级),丫2(二级)』(三级),v4(四级)} (2)根据国家(GB162972-1996)制定的大气污染物综合排放标准，评价集四级各污染物的评价标准值见表3。表3污染物各级标准值(mg/L)污染物名称污染物浓度一级二级三级四级二氧化硫SO2<0.02<0.06<0.10>0.10二氧化氮no2<0.05<0.05<0.10>0.10总悬浮颗粒PM10<0.08<0.20<0.30>0.30c)单因素评价我们采用统计分析法进行单因素评价。设对某一因素i的测量取样总次数为n，经分析其中属于一、二、三、四级的次数分别为n1,n2,n3,n4次，并且有口=、+%+%+%，则单因素的评价结果为:TOC\o"1-5"\h\zr1=(ri1，ri2，ri3，ri4)=(n1/n，n2/n，n3/n，n4/n) (3)单因素的综合分析评价结果,采4用下式计算：H=^kg (4)i ik由⑶、⑷式，根据最大隶属原则，可以判定单因素的评价等级。Hi为第i个单因素评价结果的级别变量特征值,该值越小评价结果越优。d)空气环境质量的综合评价由单因素评价结果向量，组成模糊综合评价矩阵RR-(r) (5)由于各项评价指标对空气环境质室综合评价结果的作用程度不同，因此应给予不同的权重。设评价因素的权重向量为：A-(a,a,a) (6)权重确定的方法很多，此处采用1层2析3分析法中，通过成对比较法、用1-9及其倒数作为标度构造成对比较矩阵，经过一致性比率检验通过，得到各评价因素权重。由(5)、(6)式，得空气环境质量模糊综合评价矩阵由(7)式，根据最大隶属原1则2，3可判4定空气环境质量的评价等级。式中“o”为模糊综合评价算子，选取M(A,V),计算公式如下：为克服最大隶属原则的局限性，可4计算级别变量的特征值：H-^kgB (8)k由H可以判定空气环境质量的综合评价结果，比较各地H值，由前叙述知，H值越小，评价结果越优，即空气质量越好。2)模糊评价模型的求解a)确定各评价因素权重空气环境质量由SO2、NO2、PM10三个因素共同影响。通过查阅资料，我们发现这三种大气污染物对人体的危害程度不同，对人影响最大的是SO2，其次是NO2,PM10。因此我们认为SO2最为重要，NO2次之，PM10再次之。得到成对比较矩阵A：由成对矩阵A计算得出权向量ai,最大特征根入，一致性指标CI,如下表

a1a2a3入CI0.64830.22970.12203.00370.0018查表知n=3时随机一致性指标RI=0.58,一致性比率为:一致性比率检验通过，A矩阵的特征向量A=(0.6438,0.2297,0.1220)即为SO2、NO2、PM10的权重。b)求解模型针对A、B、C、D、E、F六地的空气环境质量综合评价，采用2003年至2010年8年测量的污染物数据，利用上述模糊综合评价方法，计算得各地各项污染物的单因素评价结果(以A地为例)见表4,以及考虑三项污染物权重的综合评价结果，具体见表5。表4A地空气环境质量评价结果污染物名数据个数各级别的隶属度一级二级三级四级HiSO210840.2050.4200.1460.2302.404NO210840.8390.0000.1500.0101.331PM1010840.4800.4460.0640.0111.606综合评价结果32520.3190.4200.1500.2302.529表5各地空气质量综合评价地区综合评价指标HSO2NO2PM10三项总评A地2.4041.3311.6062.529B地2.5911.2351.7672.691C地2.1491.4101.8182.225D地2.4571.5151.8792.845E地1.7431.1081.5301.731F地2.7601.1201.8003.0494.结果分析由各地空气质量综合评价表可以得出一下结论a)A、B、C、D、E、F六个城市中，空气质量由好到差排序，依次为E、C、A、B、D、F，其中E、C地达到二级标准，A、B、D、F地为三级标准；F地由于只有2004年的数据，而通过其他城市发现，随年份增长，污染物浓度有下降趋势，因此在评价时F地不占优势。

b）模糊综合评价方法评价的结果，不仅能够得到单项污染物的评价结果，而且考虑到各项污染物的权重，得到某一监测站的综合评价结果。本项目采用隶属度的级别变量特征值作为综合评价结果的指标，克服了最大隶属原则的不足，能够包含所有评价结果的信息，具有代表性。问题二1．问题分析通过分析数据，我们发现:1）污染物浓度随年份整体上呈现递减趋势,对月平均值分析如下图1-3，鉴于污染物与年份之间存在一定关系，用前几年数据对2010年进行预测不尽合理;图1某地SO月平均值图3图2 某地NO2月平均值某地PM1d月平均值2）污染物随月份、季节变化波动很大，进行预测时只有选取相同月份或者相同季节的数据进行计算才能保证预测的准确性；每一年内各月份的数据缺失比较严重，如数据表中只有04图1某地SO月平均值图3图2 某地NO2月平均值某地PM1d月平均值2）污染物随月份、季节变化波动很大，进行预测时只有选取相同月份或者相同季节的数据进行计算才能保证预测的准确性；每一年内各月份的数据缺失比较严重，如数据表中只有04、07、09年有关于9月份的数据；3）某些天三项污染物浓度同时为零，如A地2010年9月11、12日，因大气成分中固有一定量的硫、氮浓度，我们认为在实际中不太可能出现这种情况；同时在零点处数据突变很大，对预测带来极大的不确定因素，认为各零点为奇点，预测时应剔除以保证一定的数据稳定性。综上所述，在无奇点的条件下，我们选取各地各污染物在2010年9月内的浓度数据进行预测，即所提供的2010年9月1日至9月14日的浓度数据。对于这段时间内有奇点的日子，我们予以剔除，并由9月1日往前逆推，直至找到由9月14日为起点往上逆推的、没有奇点的、14个数据，作为预测的原始数据。同时，对于F地，由于只有04年的数据，数据缺失非常严重，我们认为通过2004年的数据预测2010年的污染浓度而得到的结果不具有实际意义，因而没有对F地进行预测。2.非等间距序列的灰色模型的建立考虑到进行过数据处理的预测数据不是按时间均匀分布的，我们应用非等间距序列的灰色模型原理，建立非等间距序列的GM（1,1）模型，具体步骤如下：⑴确定非等间距序列的间隔Aki（i=1,2,L,n）,n=14；设原始数据X（0）为非负序列：X（o）={x（o）（k）,X（o）（k）,L,X（o）（k）}12若间距Ak不为常数，即该序列为非等间距序列，其中：14其中X（0）（k）＞0；Ak=k-k1,i=2,3,L,14累加生成(1-AGO)序列X(i)=((i)(k),x(i)(k),K,x(i)(k)｝；（1）（2）求一次对X(0)作一次累加生成X(0)的1-AGO(accumulatinggenerationX(i)=((i)(k),x(1)(k),K,x(1)(k)｝，其中：operator）序歹列：x(1)(k)=£x(0)Q)Ak,i=1,2,L,ni jjj=1(2)（3）由一次累加生成序列x（1）（左）构造白化微分方程;idx()(t)(3) +ax(1)(t)=u,tG[0,+8dx()(t)(3)式(3)中：t为时间；a，u为待估参数，分别为发展灰数和内生控制灰数。(4)利用最小二乘法求辨识参数a，u；而x(1)(k)=£x(o)(k)Ak，i iii=1所以』1dx(i)(t)=x(i)(k)-x(而x(1)(k)=£x(o)(k)Ak，i iii=1所以』1dx(i)(t)=x(i)(k)-x(i)(k)=x(0)(k)Ak。k.ii+1i+1i+1设z(1)(k)是x(1)(t)在区间［k,k］上的背景值，则有:i+1i i+1所以：x(0)(k)=az(1)(k)+u,i=1,2,L。i+1i+1式(5)中z(1)(k)是x(1)(k)i+1x(1)(k)两点的平均值:i+1(5)z(1)(k1)x(1)(k)+x(1)(k),i=1,2,L。(6)可利用最小二乘法求参数a，u。设:3LL-z⑴(k)1nx(0)(2),x(0)(3),L,x(0)(n)]T则有参数辨识a，u：=(BtBI1BTy(7)(5)将参数a，u代入还原模型式，计算x(0)(k)的估计值x⑹(k)；将(7)式代入(3)式并令t=勺时，x(1)(kj=x(0)(k1)，可得:x"(k)=fx(0)(k)-u又因为：x(1)(k)-x(1)(k)=x(0)(k)Ak，所以:i+1i+1 i+13.非等间距序列的灰色模型求解以MATLAB为工具，对建立的大气污染预测模型求解，以A地为例NO2，PM10辨识参数SO2浓度预测模型：从而建立A地大气污染预报模型，如下:-0.000420.0274(8)分别得到A地SO2,x①(k):65.87554X。阿川)-65.8465iNO2浓度预测模型PM10浓度预测模型A地大气污染趋势预测模型拟合误差£(k)=x(0)(k)-x(0)(k),拟合相对£(k)3=一l^x100%，如表6x(0)(k)i表6A地拟合值与实际值的相对误差一览日期SO2浓度NO2浓度PM10浓度实际值预测值相对误差％实际值预测值相对误差％实际值预测值相对误差％2010/8/290.0290.02900.0330.03300.0730.07302010/8/300.0380.027-280.0320.030-60.0840.068-192010/9/10.0250.027100.0310.030-30.0680.06802010/9/20.0230.027190.0210.030430.0640.06862010/9/30.0290.027-50.0330.030-90.0610.068122010/9/40.0250.027100.0350.030-140.0630.06882010/9/50.0260.02760.0340.030-120.0660.06842010/9/60.0340.027-190.0450.030-330.0900.068-242010/9/70.0260.02760.0240.030250.0580.069182010/9/80.0190.028450.0160.030880.1050.069-352010/9/90.0230.028200.0280.03080.0500.069382010/9/100.0300.028-80.0240.030260.0590.069172010/9/130.0370.028-260.0470.030-360.0900.069-232010/9/140.0280.028-20.0320.030-50.0870.069-20利用建立的模型，可以得到未来一周即2010年9月15日至9月21日各个城市的SO2、NO2、PM10的预测数据，如表7-11。表7 A地各污染物预测日期城市ASO2NO2PM102010/9/150.0275860.0302830.0694702010/9/160.0275980.0303040.0695802010/9/170.0276090.0303260.0696892010/9/180.0276210.0303470.0697982010/9/190.0276320.0303690.0699082010/9/200.0276440.0303900.0700182010/9/210.0276550.0304120.070128表8 B地各污染物预测日期城市BSO2NO2PM102010/9/150.0475080.0258130.0505252010/9/160.0473200.0259580.0504742010/9/170.0471320.0261040.050423

2010/9/180.0469450.0262510.0503732010/9/190.0467590.0263990.0503222010/9/200.0465730.0265470.0502712010/9/210.0463880.0266960.050220表9 C地各污染物预测日期城市CSO2NO2PM102010/9/150.0202980.0416480.0293702010/9/160.0204900.0416710.0294682010/9/170.0206840.0416940.0295662010/9/180.0208790.0417160.0296652010/9/190.0210760.0417390.0297642010/9/200.0212760.0417610.0298642010/9/210.0214770.0417840.029963表10D地各污染物预测日期城市DSO2NO2PM102010/9/150.0120020.0122620.0698212010/9/160.0120410.0122750.0700422010/9/170.0120800.0122880.0702642010/9/180.0121190.0123010.0704872010/9/190.0121580.0123140.0707102010/9/200.0121970.0123270.0709342010/9/210.0122370.0123400.071159表11 E地各污染物预测日期城市ESO2NO2PM102010/9/150.0052390.0270920.0549802010/9/160.0053980.0271060.0551242010/9/170.0055610.0271210.0552692010/9/180.0057300.0271350.0554142010/9/190.0059030.0271500.0555592010/9/200.0060820.0271640.0557052010/9/210.0062660.0271790.055851同样，利用相同模型可以对各个气象参数进行预测，如表12表12各气象参数未来一周内预测日期气压温度湿度风速2010/9/15668.92818.81567.8721.2062010/9/16669.05818.82167.6951.2042010/9/17669.18818.82667.5181.2032010/9/18669.31818.83267.3421.2012010/9/19669.44818.83867.1661.2002010/9/20669.57818.84366.9901.1992010/9/21669.70818.84966.8151.1974.模型评价1）通过对污染物浓度及气象参数预测的误差分析，可以看出，变化平稳、波动较小的一组数据得到的预测值与实际值的相对误差较小，比如，对气压的预测，其14组数据的预测值与实际值之间的相对误差均在5%以下，拟合精度相当高，对其预报应相当准确；而变化剧烈，波动很大的数据得到的误差比较大，对预测带来一定影响;2）产生上述问题的原因之一是应用了灰色预测模型，灰色预测方法由于其模型特点,比较适合于具有指数增长趋势的实际问题,对于其它变化趋势则有时拟合灰度较大,导致精度不够，同时，对于突变点的值难以准确预测；3）我们还考虑了时间序列预报以及非线性回归的方法。但由于大气污染物浓度的变化呈现很强的非线性，传统的时间序列预报是以线性自回归滑动平均模型为基础，对线性系统有较好的结果，但不适用于非线性系统的时间序列；非线性回归可以完成非线性映射，但需要在固定基函数条件下完成逼近过程，针对大气污染预报这一复杂情况，无法给出最佳基函数，也不能给出满意的解答。因而，我们最终选择非等间距序列的灰色模型。考虑到灰色模型固有的一些问题，我们通过查找资料，发现利用小波分析建立的模型可以进一步提高预测精度；4）同时，针对灰色预测模型来说，在任何一个灰色系统的发展过程中，随着时间的推移，将会不断的有一些随机扰动或驱动因素进入系统，使系统的发展相继地受其影响。一般说来，越往未来发展，越是远离时间原点，模型的预测意义就越弱。在实际应用中，随着系统的发展，老数据的信息意义将逐步降低，在不断补充新信息的同时，及时地去掉老信息，建模序列更能反映系统在目前的特征。尤其是系统随着量变的积累，发生质的飞跃或突变时，与过去的系统相比，已是面目全非。去掉已根本不可能反映系统目前特征的老数据，显然是合理的。这样逐个滚动预测，依次递补，对大气污染的预测就可以得到更为精确的结果，也才能更好地指导环境保护。问题三.问题分析分析空气质量与气象参数之间的关系，即是分析各地SO2、NO2、PM10与气象参数之间关系，从而得出气象参数对空气质量的总体影响。我们采取多元线性回归的方法，将气压、温度、湿度、风速作为自变量，通过构建多元线性方程，找出它们与各地各项污染物浓度之间关系。观察数据发现，气压、温度、湿度、风速这四项气象参数随季节、月份变化较大，间接导致污染物浓度也会随季节、月份变化。因此，我们取各地各污染物浓度2003至2010年的8年月平均值作为因变量，同理处理气象参数，作为自变量，利用MATLAB求解，并自动进行检验。.多元线性回归模型的建立多元线性回归方程：:N（0,02）,P,pLP,Q2是未知数。设x「x2、x3、X4分别为8年月平均气压、温度、湿度、风速，y为8年月平均的某污染物浓度。

设(X1,X2,x.3,X4,y)(i=1,2,L,12)是(q,x2,x3,x4,y)的12个观测值，则满足:其中各，.相互独立，且£,:N(0,。2)。令0叫,外P2,P3,04>,e=(1£2,L,%>,则方程组用矩阵表达式为假定矩阵X的秩等于m+1=5,及列满秩，则解得 6=(XtX)-1XtY，。2的无偏估计3.线性回归的假设检验1)回归方程的显着性检验H0：01=02=L=0m=0H1：至少有一个0丰0(j=1,2,L,m)。当原假设H成立时，说明回归方程不显着，采用线性回归是不合适的。当备选假破H成立时，说明回归方程显着，采用线性回归有意义。1令y=1圮y，考虑总离差平方和nii=1其中S为剩余残差平方和，SR为回归平方和。在eH成立的条件下，可以证明且S与S相互独立，则s_ s_ S^mS~(n-m-1)e:F(m,n-m-1)。对给定显着水平。，可查表知F,m,n-m-1)，计算统计量F的数值f。若f>F,m,n-m-1),则拒绝H。，即认为各系数不为零，线性回归方程式显着地。否则接受H,即人认为线性回归方程不显着。2)复0相关系数的检验复相关系数离差平方和S越小，则复相关系数越大。0<R2W1,R2越接近于1,因变量y于各自变量x之间线性相关程度越强。4i模型求解利用MATLAB软件中regress函数直接求解，得到各地各污染物与气象条件之间的回归系数及残差、统计量F的值、负相关系数。以D地为例：SO2浓度与气象参数关系：复相关系数R2=0.7032,方程显着性检验F=4.1454,概率P{F>4.1454}=0.0493<a=0.05，其中a=0.05为检验水平。故回归方程显着。no2浓度与气象参数关系：复相关系数R2=0.5664，方程显着性检验F=2.2862，概率P{F>2.2862}=0.160064>x=0.05,其中a=0.05为检验水平。故回归方程不太显着。PM10浓度与气象参数关系：复相关系数 R2= 0.7943,方程显着性检验F=6.7P{F>6.7557}=0.01493<a=0.05淇中a=0.05为检验水平。故回归方程显着。我们将所求回归方程与真实值之间作图进行比较,如下图：图1D地SO2浓度的回归方程与真实值之间比较图2D地NO2浓度的回归方程与真实值之间比较0.080050.08005图3D地PM10浓度的回归方程与真实值之间比较由图可以清楚看出，对SO2、PM10的回归方程建立跟实际值吻合的很好，而对NO2的回归方程跟实际值之间还有显着性不高,模型需要改进。由D地的SO2、NO2、PM10的多元线性回归方程中各个气象参数的回归系数可知，气压都与D地的SO2、NO2、PM10成正相关，温度、湿度都与D地的SO2、NO2、PM10成负相关，风速与SO2、NO2成正相关，与PM10成负相关。对SO2、PM10气象因素影响程度由高到低依次为：风速、气压、温度、湿度；对no2,气象因素影响程度由高到低依次为：风速、气压、湿度、温度。同一气象因素对不同污染物的影响中，气压、湿度对SO2影响最大，温度、风速对PM10影响最大。同理得到其它各地的污染物与气象参数之间的关系：A地：SO2浓度与气象参数关系：no2浓度与气象参数关系：PM10浓度与气象参数关系：B地：SO2浓度与气象参数关系：no2浓度与气象参数关系：PM10浓度与气象参数关系：C地：SO2浓度与气象参数关系：no2浓度与气象参数关系：PM10浓度与气象参数关系：E地：SO2浓度与气象参数关系：NO2浓度与气象参数关系：PM10浓度与气象参数关系：F地：SO2浓度与气象参数关系：NO2浓度与气象参数关系：PM10浓度与气象参数关系：第四问目前对环境质量有较大影响的有粉尘、硫氧化如、氮氧化物、碳氢化合物和光化学烟雾等，我们就本文中提到的三种污染物谈谈其来源及危害及相关对策粉尘。1.相关危害可吸入颗粒物^忖10)：源自烟雾、尘埃、煤灰或冷凝气化物的固体或液态水珠，能长时间悬浮于空气中，直径为10微米或以下，能直达并沉积于肺部，而引发不良的健康反应，导致呼吸不适及呼吸系统症状、加重已有的呼吸系统疾病及损害肺部组织。颗粒物上容易附着多种有害物质，当沉积在绿色植物叶面时，干扰植物吸收阳光和二氧化碳与放出氧气和水分的过程，从而影响植物的健康生长。二氧化硫是(SO)：一种无色的反应性气体，主要来自燃烧含硫燃料。二氧化硫主要危害是形成工业烟雾，高2浓度时使人呼吸困难，当其进入大气层后，氧化为硫酸，在云中形成酸雨，对建筑、森林、湖泊、土壤危害大，还会形成悬浮颗粒物，随着人的呼吸进入肺部，对肺有直接损伤。氮氧化物(NO)：主要包括一氧化氮、二氧化氮和硝酸雾，以二氧化氮为主。氮氧化物主要来自于车辆废气、火力发电站和其他工业的燃料燃烧以及硝酸、氮肥、炸药的工业生产过程。二氧化氮是一种刺激性气体，可直接进入肺部，削弱肺功能，损害肺组织，引起肺水肿和持续性、阻塞性支气管炎，降低机体对传染性细菌的抵抗能力。二氧化氮被吸收后变为硝酸与血红蛋白结合变性血红蛋白，可降低血液输送氧气的能力，同时对心、肝、肾和造血器官也有影响。2.气象参数对大气污染的影响1)气压、气温、相对湿度和风速对SO的影响风速的大小对近地层污染物在大气中2扩散与稀释有很大的影响,风速越大,污染物浓度的扩散越快,空气的稀释能力越强,反之越弱。风向,决定污染物的迁移方向,污染源的不同,方位上所受的污染程度大不相同,在污染源的下风向,往往受害较重。由于风速的大小决定污染的程度,风向决定污染的去向。因此,对于城市建设与工业区的建设有指导作用。气温,一般夏季较高,空气平流层上下温差较大,利于空气上下浮动,利于热空气流动,便于污染物的扩散与抬升,不易形成逆温层。在冬季气温较低,气压较高,不利于SO的扩散,对环境空气有明显的污染,其空气中污染物的空气污染指数有时超过5级,出现重度2污染,对人类身体健康产生一定的影响。降水和相对湿度对SO2的排放存在一定的影响。在降水过程中，雨水可冲刷、吸附大气中的污染物SO2，降低污染物浓度。2)气压、气温、相对湿度和风速对NO的影响二氧化氮的来源主要是汽车尾气和部2分厂矿企业排放的废气等。二氧化氮扩散和转化也受到各种气象因素的影响。城市空气中二氧化氮的变化也呈规律性变化,冬季明显增高,夏季明显降低,其规律呈谷峰状。3）气压、气温、相对湿度和风速对PM10的影响PM10的来源:春、夏季气温较高,除部分来源于污染源之外,主要来源于裸露在地面的扬尘的和客源的污染。如春季沙尘暴中的扬尘来源于沙漠地带。在冬季北方地区,气温较低,大部分颗粒物来源于取暖锅炉和工业污染源3.相关建议1）对SO的控制①限制高硫劣质煤品的使用；②积极推广固硫型煤，要求企业采取除尘脱硫措施，大幅度减少SO2的排放量；③尽快提高城市气化率，减少居民生活用煤；④禁止三产业,特别是餐饮业的小炉、小灶使用原煤，鼓励其使用清洁能源；⑤提高锅炉脱硫除尘效率,降低二氧化硫的排放量2）对NO的控制空看中NO2的污染主要来自机动车辆尾气的排放，随着人民生活水平的提高各种机动车辆在不断增加。加强机动车辆监管;尾气排放超标的机动车严禁上路；加强机动车辆检测,严格要求尾气超标的机动车进行治理；积极开发清洁能源,减少二氧化氮的排放。3）对PM10的控制①加强锅炉、窑炉的消烟除尘工作，使之燃煤废气达标排放;，改善落后的燃烧方式，降低烟尘的排放量；②改善能源结构，提高清洁能源比例；③加快热电联网、集中供热的步伐；④加强餐饮业油烟的整治工作，安装油烟净化装置；⑤加强港务局等企业的除尘设备的运行和管理，强化进港煤车从卸煤的翻车机房、皮带转接塔、皮带机、取料机、到装船机全部过程的规范操作和现场管理；⑥督促各类建筑施工工地强化日常管理，做到施工标准化、围墙景观化、工地管理规范化；⑦要求高扬尘污染工地实施围档作业,设置防尘设施,防止粉尘飘散,风力大于四级时应停止施工；⑧在道路清扫前进行洒水或喷雾，遇连续晴天增加喷洒次数；⑨加强对清运车辆固体散装料的密闭化的管理，对沿路泄露、遗撒的车辆进行处罚；加强城市生态绿化建设,提高市区绿化覆盖率。六、模型优缺点分析模糊综合评价方法评价的结果,不仅能够得到单项污染物的评价结果，而且考虑到各项污染物的权重，得到某一污染物的综合评价结果。采用隶属度的级别变量特征值作为综合评价结果的指标，克服了最大隶属原则的不足，能够包含所有评价结果的信息，具有代表性。灰色预测方法由于其模型特点,比较适合于具有指数增长趋势的实际问题,对于其它变化趋势则有时拟合灰度较大,导致精度不够，同时，对于突变点的值难以准确预测多元线性回归能找出多个变量对同一因素影响，建立回归方程，研究因变量与自变量之间相关性，利用MATLAB还能自动检验相关性与显着性。附录一MATLAB代码问题一.API值API=[50100200300400500];SO2=[0.050.150.81.62.12.62];NO2=[0.080.120.280.5650.750.94];PM10=[0.050.150.350.420.50.6];WuRanwu=[SO2;NO2;PM10];a=xlsread('数据附件.xls','sheet1','C3:Kl086');fori=1:3fornum=1:length(a)forj=(i-1)*3+1:(i-1)*3+3fortemp=1:5if(a(num,j)>=WuRanwu(j-(i-1)*3,temp)&a(num,j)<=WuRanwu(j-(i-1)*3,temp+1))==1a(num,j)=(API(temp+1)-API(temp))/(WuRanwu(j-(i-1)*3,temp+1)-WuRanwu(j-(i-1)*3,temp))*(a(num,j)-WuRanwu(j-(i-1)*3,temp))+API(temp);endifa(num,j)<WuRanwu(j-(i-1)*3,1)a(num,j)=50;endendendendendAPIRE=zeros(length(a),6);temp1=1;fori=1:2:5fornum=1:length(a)j=(temp1-1)*3+1;ifa(num,j)<=a(num,j+1)APIRE(num,i)=a(num,j+1);APIRE(num,i+1)=2;elseAPIRE(num,i)=a(num,j);APIRE(num,i+1)=1;endifAPIRE(num,i)<=a(num,j+2)APIRE(num,i)=a(num,j+2);APIRE(num,i+1)=3;endendtemp1=temp1+1;endxlswrite('数据附件.xls',APIRE,'sheet4','Al');.计算模糊评价参数HA=[]；%此处为SO2、NO2、MP10在大气质量中max1=[0000];max2=[0000];max3=[0000];max4=[0000];max5=[0000];max6=[0000];Rl=xlsread('数据附件，'sheet3','C4:F6');R2=xlsread('数据附件'，'sheet3','K4:N6');R3=xlsread('数据附件'，'sheet3','C12:Fl4');R4=xlsread('数据附件','sheet3','K12:N14');R5=xlsread('数据附件'，'sheet3','C20:F22');R6=xlsread('数据附件','sheet3','K20:N22');fori=1:4forj=1:3ifA(j)<R1(j,i)R1(j,i)=A(j);endendmax1(i)=max(R1(:,i));endfori=1:4forj=1:3ifA(j)<R2(j,i)R2(j,i)=A(j);endendmax2(i)=max(R2(:,i));endfori=1:4forj=1:3ifA(j)<R3(j,i)R3(j,i)=A(j);endendmax3(i)=max(R3(:,i));endfori=1:4forj=1:3ifA(j)<R4(j,i)R4(j,i)=A(j);endendmax4(i)=max(R4(:,i));endfori=1:4forj=1:3ifA(j)<R5(j,i)R5(j,i)=A(j);endendmax5(i)=max(R5(:,i));endifA(j)<R6(j,i)R6(j,i)=A(j);endendmax6(i)=max(R6(:,i));endh1=max1(1)+max1(2)*2+max1(3)*3