数学建模获奖论文

word文档精品文档分享****大学〔第七届〕数学建模竞赛承诺书本次竞赛的章程与?全国大学生数学建模竞赛章程?〔以下简称为“竞赛章程和参赛规那么〞，可从全国大学生数学建模竞赛下载〕一样，我们已认真阅读。我们明白，在竞赛开场后参赛队员不能以任何方式〔包括、电子、网上咨询等〕与队外的任何人研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反参赛规那么的。如果引用别人的成果或其他公开的资料〔包括网上查到的资料〕，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们X重承诺，严格遵守竞赛章程和参赛规那么，以保证竞赛的公正、公平性。如有违反参赛规那么的行为，我们的比赛成绩将被取消。我们参赛选择的题号是〔从A/B中选择一项填写〕：B参赛队员(打印并签名)：12.3.〔论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。〕为汇总参赛信息，请认真填写以下电子表格〔不是手写〕，如表格空间不够，可缩小字体：XX学院专业学号手机X健胡静娴刁梦雯word文档精品文档分享B题摘要本文主要分析津京冀地区主要污染源扩散〔分为点源扩散模型和线源扩散模型〕及其对大气空气质量的影响，在现有的国家最新空气污染物监测标准〔HJ633-2021环境空气质量指数(AQI)技术规定〕的根底上利用层次分析原理对现有空气质量计算模型进展改良。利用改良后的评价模型对空气质量进展评估，并根据污染源的相关特征提出可行性建议与措施。对于问题一，联系实际说明现有模型的局限性，然后引入层次分析法根据污染物的危害性对不同的污染物赋予相应的权重得到优化后的空气质量评价模型。对于问题二，通过相关查的津京冀地区的主要污染源以及污染源的主要成分，分析污染源的种类、成因、性质、危害。对于问题三，以高斯烟羽扩散模型作为工厂污染气体排放的根本模型，并作相应修改使之适用于高架点源扩散，核心问题在于解决无污染气体排放阶段到有污染气体排放阶段中间过渡段污染气体浓度的变化，最终得出了不同时刻工厂周围污染源的浓度值。对于问题四，以高斯线源扩散为根本模型，综合考虑二环、四环、六环所处位置周围地形地貌的差异，风向的影响以及城区与郊区不同扩散因素，车流量情况等求出浓度分布与梯度分布，并以图表的方式形象的展示了污染物浓度随时间地点的梯度变化，从理论上证明了政府限车令对于缓解环境压力的作用。在前四问的根底上，我们综合统筹模型的优缺点，结合津京冀地区开展实际，提出了一系列有助于蓝天方案的详实措施。AQI高斯点源扩散高斯线源扩散污染物津京冀层次分析法1问题重述近年来，我国GDP持续快速增长，对生态平衡和自然环境造成一定的破坏，空气污染的加重，雾霾天气严重影响人们的生活和XX。一般认为影响空气质量的主要因素有PM2.5、PM10、二氧化氮、二氧化硫、一氧化碳、臭氧、硫化氢、碳氢化合物和烟尘等，以京津冀地区为研究对象解决以下问题：1word文档精品文档分享〔〕参考现有国标和美标，建立衡量空气质量优劣程度等级的数学模型。〔〕分析京津冀地区主要污染源及其污染参数，分析影响空气质量的主要污染源的性质和种类。〔描述其对周围空气污染的动态影响规律。现有XX境内某一工厂废气排放烟囱高50m，主要排放物为氮氧化物。早上9点至下午3点期间的排放浓度为406.92mg/m31200m3/h10点-凌晨4点期间的排放浓度为1160mg/m35700m3/h型求解该工厂方圆51公里分别在早上8129点空气污染浓度分布和空气质量等级。〔〕建立多污染源空气污染扩散模型，并以汽车尾气污染源为例求解分析以下问题：在2021年1月15日已经连续三天发生重污染，假设从16日开场北京启动汽车单双号限行交通管制措施，求解市二环、四环、六环路在16日早上8点、中午12点、晚上9点时空气污染浓度梯度变化及空气质量等级。〔〕根据建立的模型和求解结果，分析总结影响空气质量的关键参数，为京津冀地区环保部门撰写一份建议报告，给出实现“APEC2根本假设①风的平均流场稳定，风速均匀，风向平直；②污染物的浓度在y、z轴方向符合正态分布；③污染物在输送扩散中质量守恒；④污染源的源强均匀、连续；⑤假设地面对放射性气体起全反射作用3符号说明I空气质量指数，即AQI，输出值σz垂直方向的标准差，即z方向的扩散参数C污染物浓度，输入值u平均风速，m/slow小于或等于C的浓度限值，常量Δh是热烟流的浮升力和烟气以一定速度竖直离开排放口的冲力使烟流抬升的一个附加高度high大于或等于C的浓度限值，常量qμg/s2word文档精品文档分享Ilow对应于Clow的指数限值，常量Ihigh对应于Chigh的指数限值，常量σy分别为水平方向的标准差，即y方向的扩散参数4问题分析问题一的分析我国现有空气质量计算模型主要依据?HJ633-2021环境空气质量指数(AQI)技术规定?里一维插值分析法。一维差值公式，可计算出每一种污染物对应的AQI:注：中美两国浓度限值见附录一我国的空气质量评价表与美国的空气质量评价表有一点区别，主要在于评价空气质量标准中各污染物参考上限值不一样，美国的空气质量优良的标准偏高，中国针对津京冀地区而言，我们应以国标为重，美标为辅建立空气质量评价标准。根据美标和国标AQIAQI的最大值来评价整体的空气质量。公式如下：,IAQI2,IAQI3,IAQI......IAQIn}这就出现了一个问题：以津京冀地区某一时刻的各主要污染气体的IAQ为例〔表3PM2.5、PM10在不同的时刻均保持较高值，且3的AQI分指数也相对高，如果仅仅按照最大值来评估，该最大值只是其中的一种污染物，其他污染物对空气造成的影响没有考虑在内。试想想，如果有一个地区各项AQI分指标为〔75，，，，，10AQI=75。而另一个地区的各项AQI分指数为〔75，75，75，75，75，75〕那么按国标算法公布的AQI也应为75这两个地区的环境质量差距是可观的，按照现有算法很可能会给群众带来误解，认为这两个地区环境质量相当。为了解决这个问题，我们利用了层次分析法，赋予每一个污染项一定的权值，最3word文档精品文档分享后以加权的形式给出最终的AQI值。表1：某一时刻各主要污染气体的AQIi污染物分指数〔AQI〕7:009:0012:0015:0020:00PM2.57328222275SO202101NO276661432527303340CO33336PM1015172504573问题二的分析通过网络、图书馆等途径查找我们所需要的具体数据，针对京津冀地区的地理位置和经济状况分析其主要污染源，找到具体有力的相关参数，分析出影响空气质量的污染源的性质及其种类。同时也为后面问题的解答提供根底知识，从而更好地验证第三第四小题解答的正确性。问题三的分析对于问题3，为了处理单污染源扩散模型，我们建立了高斯烟羽扩散模型，考虑到题目中要解决的是工厂烟囱排放废气问题，有效源高〔烟囱的架高与烟云抬升高度之和〕将在问题中起到一定影响，于是我们将高斯模型进展修改，综合考虑有效源高和地面反射的作用，最终确定了高架点源高斯扩散模型来对工厂污染物持续释放时间内周围污染物浓度进展求解。至于题目中要求对不同时间点浓度值的计算，此时持续排放的时刻与未进展排放时刻，工厂周围污染物浓度必然会有所不同，另外，不同的时刻大气稳定度也不一样，这会导致扩散参数的不同，为了区分这些不同，建立了高斯扩散下排放源强为0mg/s时的污染物扩散模型。问题四的分析对于建立多污染源空气污染扩散模型，我们只需在问题三的根底上，将高斯点源扩散模型转化为高斯线源扩散模型即可，线源模型下的污染物浓度C即为高斯点源模型下污染物浓度沿线源方向的定积分或广义积分。为了解决题目中政策实施4word文档精品文档分享后，不同地点，不同时间污染物浓度C的变化，对于不同环路早上8点，中午12点，晚上9点这三个特殊时间点我们综合考虑了气象〔大气稳定度〕和社会方面〔上下班顶峰期等〕的因素，从而使算出的数据更具实际意义。最后，通过对空气质量等级的评估，来阐述政府调控在保护环境当中的重要作用。问题五的分析在问题一中我们通过层次分析法建立了新的空气质量评价标准，根据表3中的数据，我们肉眼可以直观的看出PM2.5、PM10、3等污染气体浓度较高。问题二中对津京冀地区分别进展分析得出：表2对津京冀地区主要污染源的统计地区主要污染源机动车尾气排放XX工业废气污染XX燃煤、工业废气污染津京冀地区的主要污染原因不尽一样，问题三、四针对和XX举出典型的两个例子分析浓度分布变化，我们可以直观的看出不同时间段浓度变化的程度并可以在此题中提出针对各个时间段的改善措施。模型的建立与求解问题一的建模和求解1国标与美标浓度限度的比拟见附录一2在研究各种指标集对某物影响的过程中，不仅指标集中的最大值具有最重要的作用，次大值等的作用也不容无视，甚至具有与最大值类似的影响。因此，为了在评价体系中引入次大值等的影响。我们采用层次分析加权法得到各个污染物IAQI所占的权重。AQI1*IAQI1a2*IAQI23*IAQI3a4*IAQI45*IAQI5a6*IAQI66ai1i1ai即为各个污染物IAQI的权值。具体方法如下：1构造判断矩阵通过对指标之间两两重要程度进展比拟和分析判断，构造判断5word文档精品文档分享矩阵。层次分析法在对指标的相对重要程度进展测量时，我们引入了九分位的相对重要的比例标度。令A为判断矩阵，用以表示同一层次各个指标的相对重要性的判断值。根据心理学家提出的“人区分信息等级能力为〞的研究结论，有如下评分规那么：表3权重的评分规那么不甲指标与乙指标极端强烈明显比拟重较不很不极不重比拟重要重要重要重要要重要重要重要要甲指标评分975311\31\51\71\9注：取8,6,4,2,1/2,1/4,1/6,1/8为上述评价值的中间值。在综合考虑了京津冀地区的污染情况后，我们建立了比拟矩阵A如下：表4空气质量指数空气评价判断矩阵AO3SO2NO2COPM10PM2.5O314331/21/2SO21/413/43/41/81/8NO21/34/3111/61/6CO1/34/3111/61/6PM10286611PM2.52866113对判断矩阵进展一致性检验：3.1用matlab求出判断矩阵A的特征向量和最大特征值，最大特征值对应的特征向量归一化后即为各分指标IAQI的权重。3.2计算判断矩阵的一致性指标。具体公式为：CI=(λ-n)/(λ-1)6word文档精品文档分享CI为0时，A一致；CI越大，A的不一致程度就越高。表5随机一致性指标RI：N1234567891011RI000.580.901.121.241.321.411.451.491.51CR=CI/RI,当CR<0.1时，A的不一致性程度在允许X围内，此时可用A的特征向量作为权向量。表6各指标权重和一次性指标值O3SO2NO2COPM10PM2.5CICRλ0.16900.04230.05630.05630.33800.3380-3.552-2.86567e-161e-16然后把各污染物的IAQI与其相应的权重相乘得到综合指数AQI。注：中国国家AQI评价标准见附录二表7各污染物AQI分指标SO2NO2O3COPM2.5PM10地区一5830185073地区二32523073地区一与地区二按国标算法两者AQI环境优于地区一，因此国标算法在这个问题上存在缺陷。而依据我们改良的算法，地区一和地区二的AQI分别为48.32和36.00地的环境差异。当然基于此算法环境优良等级也应做相应调整。问题二的建模和求解：有关数据显示，京津冀地区主要污染源分别为机动车、燃煤、扬尘、流动源以及工业生产，其中燃煤和机动车占主要因素，见以下表格及饼图。近几年，主要污染源为机动车，见以下图及表格。2021年，京津冀机动车氮氧化物排放量68.2万吨，占氮氧化物排放总量的30%，其中机动车氮氧化物排放量占本地区氮氧化物的比重达45%28.8个和XX13.9个百分点。7word文档精品文档分享XX工业污染影响最大,数据显示，2021年，京津冀工业二氧化硫排放量占二氧化硫排放总量的91.2%68.4%；工业烟〔粉〕尘排放量占烟〔粉〕尘排放总量的82.6%。分省市看，XX工业污染影响最大。XX燃煤污染尤其严重。XX二氧化硫排放占京津冀八成,数据显示，2021年，京津冀燃煤消费总量38927万吨，XX煤炭消费量占其能源消费总量的88.8%，远远高于的25.4%和XX的59.6%。煤炭消费排放出大量二氧化硫，对大气环境造成很大影响，2021年XX二氧化硫排放量占京津冀的80.8%。表8主要污所含主要污染成分地点危害染源固体悬浮微粒、一氧化碳、碳氢汽车尾造成雾霾，颗粒物对人体危害化合物、氮氧化合物、铅及硫氧气大化合物等工业生SO2、粉尘、NOx对工农业建立和人民XX危害XX产极大XX燃煤SO2、CO、粉尘、NOx酸雨，破坏生态环境表9京津冀区域ＰＭ2.5来源所占比例燃煤34%机动车16%工业、外来输送、扬尘、餐饮和其他50%34%16%50%表10雾霾天气的成因汽车尾气燃煤污染建筑扬尘工业排放22.20%16.70%16.30%15.70%8word文档精品文档分享图1京津冀区域ＰＭ2.5来源所占比例图2雾霾成因由以上查找数据可以总结：主要污染源机动车尾气的所含污染成分为CO和氮氧化物，燃煤所含污染物为二氧化硫和一氧化碳，工业生产含有的污染成分为二氧化硫和氮氧化物及粉尘。问题三的建模和求解：对于单污染源的扩散，我们以高斯烟羽扩散为根本模型：此题为高架点源的高斯扩散模式。有效源位于坐标原点ox轴平行，并与x轴正向同向。假设点源在没有任何障碍物的自由空间扩散，不考虑下垫面的存在。大气中的扩散是具有y与z两个坐标方向的二维正态分布，当两坐标方向的随机变量独立时，分布密度为每个坐标方向的一维正态分布密度函数的乘积。由正态分布的假设条件②，参照正态分布函数的根本形式，取μ＝0，那么在点源下风向任一点的浓度分布函数为：〔3-1〕式中C—空间点〔x，y，z〕的污染物的浓度，mg/m3；A〔x〕—待定函数；σyz—分别为水平、垂直方向的标准差，即yx方向的扩散参数，。由守恒和连续假设条件③和④，在任一垂直于x轴的烟流截面上有：〔3-2〕式中q—源强，即单位时间内排放的污染物，μg/s；u—平均风速，m/s。将式〔3-1〕代入式〔3-2〕,由风速稳定假设条件①，A与y、z无关，考虑到③和④，积分可得待定函数〔x〕：9word文档精品文档分享〔3-3〕将式〔3-3〕代入式〔3-1〕，得大空间连续点源的高斯扩散模式〔3-4〕式中，扩散系数σy、σz与大气稳定度和水平距离x有关，并随x的增大而增加。为了解决题中烟囱排放污染物问题，我们对高斯点扩散模型做了一些修改，得到适用于高架源的高斯高架源扩散模型见以下图。点源在地面上的投影点o作为坐标原点，有效源位于z轴上某点，z＝。高架有效源的高度由两局部组成，即＝h＋Δh，其中h为排放口的有效高度，Δh是热烟流的浮升力和烟气以一定速度竖直离开排放口的冲力使烟流抬升的一个附加高度，如图3所示。当污染物到达地面后被全部反射时，可以按照全反射图3地面全反射高架点源模型原理，用“像源法〞来求解空间某点k的浓度。图3中k点的浓度显然比大空间点源扩散公式〔3-4〕计算值大，它是位于(0，，H)的实源在k点扩散的浓度和反射回来的浓度的叠加。反射浓度可视为由一与实源对称的位于(0，，－H)的像源〔假想源〕扩散到k点的浓度。由图可见，k点在以实源为原点的坐标系中的垂直坐标为(z-H)k点扩散的浓度为式〔3-4z轴向下平移距离：〔3-5〕k点在以像源为原点的坐标系中的垂直坐标为(z＋H)，那么像源在k点扩散的浓度为式〔3-4〕的坐标沿z轴向上平移距离：10word文档精品文档分享〔3-6〕由此，实源s与像源x之和即为k点的实际污染物浓度：〔3-7〕假设污染物到达地面后被完全吸收，那么Cx＝0，污染物浓度〔x，y，z，〕＝Cs，即式〔3-5〕。〔2实际中，高架点源扩散问题中最关心的是地面浓度的分布状况，尤其是地面最大浓度值和它离源头的距离。在式〔3-7〕中，令z＝，可得高架点源的地面浓度公式：〔3-8〕上式中进一步令y＝0那么可得到沿x轴线上的浓度分布：〔3-9〕考虑到题目中对于不同时间点污染物浓度的计算，我们对高斯扩散模型进展了一些修改，得到了与时刻t有关的高斯扩散模型〔此模型更加适用于源强为0时的高斯扩散〕:〔3-10〕题中问题求解：1.参数确实定。高斯扩散公式的应用效果依赖于公式中的各个参数的准确程度，尤其是扩散参数σ、σz及烟流抬升高度Δh的估算。其中，平均风速u取多年观测的常规气象11word文档精品文档分享数据；源强q可以计算或测定，σy、σz及Δh与气象条件和地面状况密切相关。Δh确实定：Δh可用如下公式近似计算：在查阅了XX省历年平均风速以及污染物排放指数后，大计算出Δh=12.8mH=h+Δh=50+12.8=62.8m扩散参数σy、σz确实定：σy、σz确实定依赖于当地的风速，太阳照射角，云量，对于它们确实定一般都是在确定了稳定度〔依据大气稳定度的等级表〕的情况下，根据扩散参数函数表〔添加表名称〕来确定σy、σz关于X的具体表达式。表11为Pasquill稳定度分级法表3.3Pasquill稳定度分级法地面日间日射强度夜间天空状况风速〔m/s〕强中等弱薄云遮天云量或云量≤4/8≤3/8<2AA-BB2-3A-BBCEF3-5BB-CCDE5-6CC-DDDD>6CDDDD不管何种天空状况，夜晚及最后稳定度为中。不管何种天空状况，日落前一个小时和日出后一个小时稳定度为中，且仲夏晴〔日落前至日出后1h为夜晚〕。仲夏晴天中午为强日照，寒冬晴天中午为弱日照。天中午为强日照，寒冬晴天中午为弱日照。表12横向扩散参数幂函数表达式数据12word文档精品文档分享表3.11-a横向扩散参数幂函数表达式数据扩散参数稳定度等级11下风距离A0.9010740.8509340.4258090.6020520~1000>1000B0.9143700.8650140.2818460.3963530~1000>1000B~C0.9193250.2295000~10000.8750860.314238>1000C0.9242790.8851570.1771540.2321230~1000>1000y=x1C~D0.9268490.8869400.1439400.1893960~1000>1000D0.9294180.1107260~10000.8887230.14669>1000D~E0.9251180.8927940.09856310.1243080~1000>1000E0.9208180.08640010~10000.8968640.101947>1000F0.9294180.05536340~10000.8887230.0733348>1000表13Briggs扩散参数表Briggs扩散参数(单位：m)(a)开阔乡间条件稳定度级别yzA0.22x/(1+0.0001x)-1/20.20xB0.16x/(1+0.0001x)-1/20.12xC0.11x/(1+0.0001x)-1/20.08x/(1+0.0002x)-1/2D0.08x/(1+0.0001x)-1/20.06x/(1+0.0015x)-1/2E0.06x/(1+0.0001x)-1/20.03x/(1+0.0003x)-1F0.04x/(1+0.0001x)-1/20.016x/(1+0.0003x)-1(b)城市条件yzA-B0.32x/(1+0.0004x)-1/20.24x/(1+0.001x)-1/2C0.22x/(1+0.0004x)-1/20.20xD0.16x/(1+0.0004x)-1/20.14x/(1+0.0003x)-1/2E-F0.11x/(1+0.0004x)-1/20.08x/(1+0.0015x)-1/2根据XX省气象站1971-2021年统计数据，我们大致确定出早上8点大气稳定度等级为，中午12点大气稳定度等级为，晚上9点大气稳定度等级为，那么不同时刻不同地点σy、σz可以相继算出。确定参数后，我们用matlab对题中问题进展了简答并对不同时刻工厂周围的污染物浓度作出图形，以便形象的观测。13word文档精品文档分享-3早上8点污染物浓度分布x109876C度5浓物染污4321000.511.522.5距工厂距离X4x10图4早上8点污染物浓度分布3.5x10-3中午12点污染物浓度分布32.5C度浓物2染1.5污10.5000.511.522.5距工厂距离X4x10图5中午12点污染物浓度分布14word文档精品文档分享-3晚上9点污染物浓度分布x1032.52C度浓物染污1.510.5000.511.522.5距工厂距离X4x10图6晚上9点污染物浓度分布从图中我们大致可以得到结论：污染物浓度整体变化趋势是随着距工厂的距离的增加递减的，另外，从中午１２点与早上８点同一地点的污染物浓度比拟来看，早上８点污染物的浓度比起中午１２点要低，这说明正在排污时污染物对周围环境的破坏是最大的，因此从保护环境的角度来看，限制污染源的排放是一个很有效的措施。对于3个点质量等级的判断，我们不妨取3个时刻污染物浓度最大的地点进展计算，因为污染物浓度最大的地点是最容易形成污染物聚集而对环境造成破坏的。根据matlab所计算数据，IAQI计算公式以及第一问建立的AQI模型可以对题中3个时间点进展环境质量等级评估。AQIaIAQIaIAQIaIAQIaIAQIaIAQIaIAQI1*12*23*34*45*56*66ai1i1早上８点浓度最大值为：0.0082/m3其对周围氮氧化物IAQI的奉献值为4.1；15word文档精品文档分享中午１２点浓度最大值为：0.035mg/m3其对周围氮氧化物IAQI的奉献值为17.5；晚上９点浓度最大值为：0.0028mg/m3其对周围氮氧化物IAQI的奉献值为1.4。根据XX环境监测局数据XX市郊区各个污染源平均IAQI分布如下：表14XX省郊区各污染物AQI分指标PM2.5SO2NO2O3COPM107317253072根据matlab计算后可知污染物排放后各个污染源AQI分布如下：表15XX省郊区各污染物AQI分指标PM2.5SO2NO2O3COPM108:0073111.125307212:0073124.525307221:007318.4253072表16根据题一所建模型，各时刻AQI值如下AQI排污前54.96638:0055.22312:0055.977521:0055.071注：对应代码见附录三问题四的建模求解在问题三的根底上，我们将高斯点源模型扩展为高斯线源模型，并对不同情况下线源模型的适用进展分类说明〔为了简化计算，以下模型均建立在风向与线源垂1.有风时(u>0.4m/s)采用高斯烟流模式设x轴与风向一致，坐标原点在线源的中点，线源在y轴上的长度为2y0。由它所产生的在任一点〔x,y,0〕的浓度为：16word文档精品文档分享C(x,y,0；h)=错误！未找到引用源。*exp(-错误！未找到引用源。)*错误！未找到引用源。(4-1)式中；u为平均风速〔m/sy、σz为水平和垂直方向的扩散参数；H为有效源高；C为污染物浓度；q为线源强度。当为无限长线源时，地面浓度公式为：(4-2)2.准静风时(u<0.4m/s)采用高斯烟团模式通常平均风速小于1m/s时，适用于连续点源的准定常烟流模式不再适用。此时用高斯烟团模式代替烟流模式。取y轴与线源一致，坐标远点与线源重点重合时，对于线源排放高度为h的地面浓度计算公式为：(4-3)当为无限长线源时，准静止风条件下的地面浓度计算公式为：(4-4)参数确实定：1.有风时，参数σy、σz确实定与问题三中确定方法一样，在这不作重复说明。2.准静风时〔u<0.4m/s〕扩散参数确实定可参照Turner分类法进展的Pasquill稳定度分类表。2.线源源强的计算车辆平均排放因子*车流量排放因子：指车辆行驶单位距离平均排出污染物的量车流量：每小时通过的车辆数目实例解答：在查阅了历年气象资料后，发现的情况适合用有风时的高斯线源模型来17word文档精品文档分享求解。通过查阅相关资料，我们得出了二，四，六环在题中3个时间点的车流量如下表所示：表17政策实施前车流量单位〔辆/小时〕环数时8:0012:0021:00间二环1200060008600四环600030004000六环300020002000表18政策实施后车流量单位〔辆/小时〕环数时8:0012:0021:00间二环800040005600四环400020003000六环200015001500表19车辆平均排放因子为：汽车尾气中不同污染物排放因子单位〔毫克/米.辆〕污染物CONOx排放因子40.052.03二环：综合考虑地理位置和时间的影响，二环早上8点大气稳定度为12点为，晚上9点为E根据不同时刻大气稳定度可确定出不同地点σy、σz。用matlab将各数据带入18word文档精品文档分享高斯线源扩散模型即可求出各个时刻的汽车尾气浓度变化。二环CO污染物浓度曲线12蓝色8点红色12点黄色21点108C度浓OC642001002003004005006007008009001000距线源距离x图7二环浓度曲线二环NOx浓度曲线0.35蓝线8点0.3红线12点黄线21点0.25C0.2度浓xON0.150.10.05001002003004005006007008009001000距线源距离x图8二环NOx浓度曲线19word文档精品文档分享四环：综合考虑地理位置和时间的影响，四环早上8点大气稳定度为12点为，晚上9点为D根据不同时刻大气稳定度可确定出不同地点σy、σz。用matlab将各数据带入高斯线源扩散模型即可求出各个时刻的汽车尾气污染物浓度分布情况。四环CO浓度曲线3蓝线8点红线12点黄线21点2.52C度浓OC1.510.5001002003004005006007008009001000距线源距离x图9四环浓度曲线四环NOx浓度曲线0.16蓝线8点红线12点0.14黄线21点0.120.1C度浓xON0.080.060.040.02001002003004005006007008009001000距线源距离x20word文档精品文档分享图10四环浓度曲线六环：综合考虑地理位置和时间的影响，六环平均风速2.1m/s，六环早上8点大气稳定度为，中午12点为，晚上9点为D根据不同时刻大气稳定度可确定出不同地点σy、σz。用matlab将各数据带入高斯线源扩散模型即可求出各个时刻的汽车尾气污染物浓度分布情况。六环CO浓度曲线1蓝线8点0.9红线12点黄线21点0.80.70.6C度浓OC0.50.40.30.20.1001002003004005006007008009001000距线源距离x图11六环浓度曲线21word文档精品文档分享六环NOx浓度曲线0.05蓝线8点0.045红线12点黄线21点0.040.0350.03C度浓xON0.0250.020.0150.010.005001002003004005006007008009001000距线源距离x图12六环浓度曲线分析：经过上面几幅图的比照，我们发现，就同一地点而言，比方二环，早晨8点时刻比其它两个时刻CO和NOx的浓度都要高，这正符合实际情况，因为早晨8点处于上班顶峰期，车流量的增多必然会引起污染物浓度相较其他时间偏高。而从不同地点即不同环度来比拟，同一时刻污染物浓度也不尽一样，这也正说明了气象条件〔风速〕对污染物扩散的影响。但从整体来看，车流量是影响污染物扩散的主要因素，因此单双号限行措施会在很大程度上缓解的空气环境压力。题中问题解答：污染物在同一地点不同时间的梯度变化可以上面图形中看出。至于空气质量等级的评估，我们更愿意去考虑线源对于周围环境影响的大小。故列表如下：表20政策实施政策实施距污染源距离X前CO浓度单项API值后浓度单项API值mg/smg/s1016.28112.5610.85101.7503.5435.42.3523.51001.7717.71.1811.822word文档精品文档分享2000.898.90.595.9从图表中可以看出，车流量的减少在一定程度上释放了周围空气环境的压力，同一时间点，政策实施前后AQI值有所改变，这在一定程度说明了限行政策的可行性。另外从空间来看，当离线源距离超过200时，线源对周围环境的影响便不再明显，而影响较大的区域为以内，这也就指导人们房屋建筑应与道路拉开一定距离，不能靠的太近。尽管车流数据不大准确，但我们还是可以看出，政府政策调控对于环境保护的有力作用。注：对应代码见附录三问题五的求解由第二问的分析以及第三问的高斯烟羽扩散模型和第四问的高斯线源模型所获得结果和数据可以直观的比照出影响津京冀地区空气质量的主要参数。表21津京冀地区主要污染源成分参数分析地区主要污染源主要成分机动车尾气排放PM2.5XX工业废气污染SO2NOxXX燃煤、工业废气污染SO2CONOx对津京冀地区环保局部的建议：地区建议根据第四问可知在晚上九点钟的时候解除限行指令；推广使用天然气；改用优质低硫低灰煤XX发电厂燃煤机组的清洁化利用、开发新型能源、严格执行相关法律法规；定期清洗马路控制粉尘XX大力开发使用新型能源，例如风能、太阳能、利用燃煤企业余热发电等，严格监视工厂做到合格排污，根据第三问分析可知，中午要注重调控工厂废气的释23word文档精品文档分享放对于APEC蓝的可行性措施和建议：车辆限行，提高车辆质量，加大对电力汽车的研发和使用力度。在华北地区建立起联动控制工期质量的工作机制，加大环境执法的执行力度。模型的优缺点及推广优点：在建模过程中，我们针对不同情况，对模型进展了相应修改，到达了根底模型的通用性和特定模型的针对性的结合，并且由模型求出的解在一定程度上与实际吻合较好，对实际生活有一定的指导意义。缺点：以上模型大多为了求解方便都是建立在理想情况下的，未能细节深入到实际细小的方面，比方污染源周围建筑物，污染物的沉降效应等都未作考虑，这就使得模型在一些特定情况下不适用。推广：以上模型还可以推广到重金属的污染扩散求解等问题，另假设将模型在做一维扩展，可以用来求解面源扩散问题。参考文献[1]姜启源2003年8月[2]谢绍东等?汽车排放污染物浓度的预测?，中国环境科学，1995年10月[3]卓金武，魏永生，秦建，李必文．Matlab在数学建模中的应用[M]航空航天大学，2021[4]X斌才，赵军．大气污染扩散的高斯烟羽模型及其GIS集成研究[J]监测管理与技术，2021，20〔517-19[5]环境空气质量指数(AQI)技术规定〔试行〕HJ633-2021[Z].：中国环境科学[6]环境空气质量标准GB3095-2021[Z].：中国环境科学附录附录一24word文档精品文档分享附录二中国空气质量等级分类空气质量指数值中国级别中国分类25word文档精品文档分享0-50一级优51-100二级良101-150三级轻度污染151-200四级中度污染201-250五级重度污染251-300301-350351-400401-450六级严重污染451-500500+附录三问题一对应MATLAB代码利用层次分析加权法计算各污染物的权重A=[1,6,2,5,1/2,1/21/6,1,1/3,5/6,1/12,1/121/2,3,1,2/5,1/4,1/41/5,6/5,2/5,1,1/10,1/102,12,4,10,1,12,12,4,10,1,1];[D,X]=eig(A);ZD=max(max(X));[m,m]=size(D);p=max(X);fori=1:mifp(i)==ZDq=iendend26word文档精品文档分享z=sum(D);fori=1:mforj=1:mD(i,j)=D(i,j)/z(j);endendB=(D(:,q))'RI=[0,0,0.58,0.9,1.12,1.24,1.32,1.41,1.45,1.49,1.51];CI=(ZD-m)/(m-1)CR=CI/RI(m)Q=zeros(11,1);fori=1:11Q(i,1)=0;forj=1:6Q(i,1)=Q(i,1)+sum(i,j)*B(1,j);endEnd问题三对应matlab码中午12点图形x=0:100:20000;oy=0.433