****大学（第七届）数学建模竞赛
承 诺 书
本次竞赛的章程与《全国大学生数学建模竞赛章程》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）相同，我们已认真阅读。
我们明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反参赛规则的。如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反参赛规则的行为，我们的比赛成绩将被取消。
我们参赛选择的题号是（从A/B中选择一项填写）： B 参赛队员 (打印并签名) ：1 2. 3. （论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。）
为汇总参赛信息，请认真填写以下电子表格（不是手写），如表格空间不够，可缩小字体： 姓名 学院 专业 学号 手机 邮箱 张健 胡静娴 刁梦雯
B题
摘要
本文主要分析津京冀地区主要污染源扩散（分为点源扩散模型和线源扩散模型）及其对大气空气质量的影响，在现有的国家最新空气污染物监测标准（HJ633-2012环境空气质量指数(AQI)技术规定）的基础上利用层次分析原理对现有空气质量计算模型进行改进。利用改进后的评价模型对空气质量进行评估，并根据污染源的相关特征提出可行性建议与措施。
对于问题一，联系实际说明现有模型的局限性，然后引入层次分析法根据污染物的危害性对不同的污染物赋予相应的权重得到优化后的空气质量评价模型。 对于问题二，通过相关网站查的津京冀地区的主要污染源以及污染源的主要成分，分析污染源的种类、成因、性质、危害。
对于问题三，以高斯烟羽扩散模型作为工厂污染气体排放的基本模型，并作相应修改使之适用于高架点源扩散，核心问题在于解决无污染气体排放阶段到有污染气体排放阶段中间过渡段污染气体浓度的变化，最终得出了不同时刻工厂周围污染源的浓度值。
对于问题四，以高斯线源扩散为基本模型，综合考虑二环、四环、六环所处位置周围地形地貌的差别，风向的影响以及城区与郊区不同扩散因素，车流量情况等求出浓度分布与梯度分布，并以图表的方式形象的展示了污染物浓度随时间地点的梯度变化，从理论上证明了政府限车令对于缓解环境压力的作用。
在前四问的基础上，我们综合统筹模型的优缺点，结合津京冀地区发展实际，提出了一系列有助于蓝天计划的详实措施。
【关键词】AQI 高斯点源扩散 高斯线源扩散 污染物 津京冀 层次分析法 1问题重述
近年来，我国GDP持续快速增长，对生态平衡和自然环境造成一定的破坏，空气污染的加重，雾霾天气严重影响人们的生活和健康。一般认为影响空气质量的主要因素有PM2.5、PM10、二氧化氮、二氧化硫、一氧化碳、臭氧、硫化氢、碳氢化合物和烟尘等，以京津冀地区为研究对象解决以下问题：
1
（1）参考现有国标和美标，建立衡量空气质量优劣程度等级的数学模型。 （2）分析京津冀地区主要污染源及其污染参数，分析影响空气质量的主要污染源的性质和种类。
（3）建立单污染源空气污染扩散模型，描述其对周围空气污染的动态影响规律。现有河北境内某一工厂废气排放烟囱高50m，主要排放物为氮氧化物。早上9点至下午3点期间的排放浓度为406.92mg/m3，排放速度为1200m3/h；晚上10点-凌晨4点期间的排放浓度为1160mg/m3，排放速度为5700m3/h；通过你的扩散模型求解该工厂方圆51公里分别在早上8点、中午12点、晚上9点空气污染浓度分布和空气质量等级。
（4）建立多污染源空气污染扩散模型，并以汽车尾气污染源为例求解分析以下问题：北京在2015年1月15日已经连续三天发生重污染，假设从16日开始北京启动汽车单双号限行交通管制措施，求解北京市二环、四环、六环路在16日早上8点、中午12点、晚上9点时空气污染浓度梯度变化及空气质量等级。 （5）根据建立的模型和求解结果，分析总结影响空气质量的关键参数，为京津冀地区环保部门撰写一份建议报告，给出实现“APEC”蓝天的可行性措施和建议。
2基本假设
①风的平均流场稳定，风速均匀，风向平直；②污染物的浓度在y、z轴方向符合正态分布；③污染物在输送扩散中质量守恒；④污染源的源强均匀、连续；⑤假设地面对放射性气体起全反射作用
3符号说明
I 空气质量指数，即AQI，输出值 σz 垂直方向的标准差，即z方向的扩散参数 C Clow 污染物浓度，输入值 小于或等于C的浓度限值，常量 u Δh 平均风速，m/s 是热烟流的浮升力和烟气以一定速度竖直离开排放口的冲力使烟流抬升的一个附加高度 Chigh 大于或等于C的浓度限值，常量 q 源强，即单位时间内排放的污染物，μg/s 2
Ilow Ihigh σy 对应于Clow的指数限值，常量 对应于Chigh的指数限值，常量 分别为水平方向的标准差，即y方向的扩散参数
4问题分析
问题一的分析
我国现有空气质量计算模型主要依据《HJ633-2012环境空气质量指数(AQI)技术规定》里一维插值分析法。一维差值公式，可计算出每一种污染物对应的AQI:
注：中美两国浓度限值见附录一
我国的空气质量评价表与美国的空气质量评价表有一点区别，主要在于评价空气质量标准中各污染物参考上限值不一样，美国的空气质量优良的标准偏高，中国偏低，这是由我国的基本国情决定的。针对津京冀地区而言，我们应以国标为重，美标为辅建立空气质量评价标准。
根据美标和国标AQI算法，观测点取各种污染物所对应AQI的最大值来评价整体的空气质量。公式如下：
AQI=max{IAQI1, IAQI2, IAQI3, IAQI4...... IAQIn }
这就出现了一个问题：以津京冀地区北京某一时刻的各主要污染气体的IAQ为例（表3），其中主要污染气体PM2.5、PM10在不同的时刻均保持较高值，且O3的AQI分指数也相对高，如果仅仅按照最大值来评估，该最大值只是其中的一种污染物，其他污染物对空气造成的影响没有考虑在内。试想想，如果有一个地区各项AQI分指标为（75，1，1，2，3，10），，那么最终公布的AQI=75。而另一个地区的各项AQI分指数为（75，75，75，75，75，75）那么按国标算法公布的AQI也应为75，很显然，这两个地区的环境质量差距是可观的，按照现有算法很可能会给群众带来误解，认为这两个地区环境质量相当。
为了解决这个问题，我们利用了层次分析法，赋予每一个污染项一定的权值，最
3
后以加权的形式给出最终的AQI值。 表1：北京某一时刻各主要污染气体的AQIi
污染物 分指数（AQI） 7:00 PM2.5 SO2 NO2 O3 CO PM10 73 0 7 25 3 151 9:00 28 1 6 27 3 72 12:00 22 1 6 30 3 50 15:00 22 0 6 33 3 45 20:00 75 1 14 40 6 73 问题二的分析
通过网络、图书馆等途径查找我们所需要的具体数据，针对京津冀地区的地理位置和经济状况分析其主要污染源，找到具体有力的相关参数，分析出影响空气质量的污染源的性质及其种类。同时也为后面问题的解答提供基础知识，从而更好地验证第三第四小题解答的正确性。 问题三的分析
对于问题3，为了处理单污染源扩散模型，我们建立了高斯烟羽扩散模型，考虑到题目中要解决的是工厂烟囱排放废气问题，有效源高（烟囱的架高与烟云抬升高度之和）将在问题中起到一定影响，于是我们将高斯模型进行修改，综合考虑有效源高和地面反射的作用，最终确定了高架点源高斯扩散模型来对工厂污染物持续释放时间内周围污染物浓度进行求解。至于题目中要求对不同时间点浓度值的计算，此时持续排放的时刻与未进行排放时刻，工厂周围污染物浓度必然会有所不同，另外，不同的时刻大气稳定度也不一样，这会导致扩散参数的不同，为了区分这些不同，建立了高斯扩散下排放源强为0mg/s时的污染物扩散模型。
问题四的分析
对于建立多污染源空气污染扩散模型，我们只需在问题三的基础上，将高斯点源扩散模型转化为高斯线源扩散模型即可，线源模型下的污染物浓度C即为高斯点源模型下污染物浓度沿线源方向的定积分或广义积分。为了解决题目中政策实施
4