空气污染问题研究-数学建模论文

PAGE2015年第十二届五一数学建模联赛承诺书我们仔细阅读了五一数学建模联赛的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与本队以外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其它公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们愿意承担由此引起的一切后果。我们授权五一数学建模联赛赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。 我们参赛选择的题号为（从A/B/C中选择一项填写）：B我们的参赛报名号为：1238参赛组别（研究生或本科或专科）：本科所属学校（请填写完整的全名）徐州工程学院参赛队员(打印并签名)：1.唐惠2.李俊3.王媛媛日期：2015年5月3日获奖证书邮寄地址：徐州工程学院中心校区敬知楼514室邮政编码：221018收件人姓名：数学建模协会联系电话：83105225

2015年第十二届五一数学建模联赛编号专用页竞赛评阅编号（由竞赛评委会评阅前进行编号）：评阅记录评阅人评分备注裁剪线裁剪线裁剪线竞赛评阅编号（由竞赛评委会评阅前进行编号）：参赛队伍的参赛号码：（请各参赛队提前填写好）：2015第十二届五一数学建模联赛题目空气污染问题研究摘要本文通过建立高斯模型对空气污染物及其扩散进行研究，求出空气污染的浓度变化和对应的空气质量等级。针对问题一，为了解决衡量空气质量优劣程度等级问题，参考美标和国标下的污染物浓度限值，引用空气质量分指数公式，建立空气质量指数模型。求出不同标准下的空气质量指数，再由美标和国标下的空气质量等级，判断出同一污染物浓度在不同标准下的空气质量优劣程度。针对问题二，运用因子分析法，根据变量共同度，求出京津冀地区的主要污染参数为、、和。由主要污染参数得出主要污染源为煤炭和石油的燃烧源、排放污染气体的化学工业、重工业。通过应用分析，按其化学性质可将主要污染源划分为：还原型污染源和氧化型污染源；按其能源性质可将主要污染源划分为四个种类：煤炭型、石油型、混合型、特殊型。针对问题三的单污染源空气污染扩散问题，建立高斯烟羽模型：利用编程求出公里范围内的沿地面轴线的污染物浓度（具体数据见正文表），并结合污染物限值表中的限值，求出在这范围内不同时刻的空气质量等级（具体数据见正文表）。针对问题四中多污染源空气污染扩散的问题，建立高斯线源模型：利用编程得出三个时间点的各处污染物浓度梯度变化情况（具体数据见正文表）。分析可知：二环和四环之间的平均距离为4.9km<5km，四环和六环平均距离达16.22>6km公里，因此将二环和四环的汽车尾气叠加，六环线路作独立线源分析。将处理后的三个环路段的浓度与美标和国标下空气质量等级对应，得出各个环路在各时刻的空气质量等级（具体数据见正文表）。关键词空气质量指数因子分析法高斯烟羽模型高斯线源模型PAGE26一、问题背景与重述1.1问题背景空气污染，又称大气污染，按照国际标准化组织（ISO）的定义，空气污染通常是指由于人类活动或自然过程引起某些物质进入大气中，呈现出足够的浓度，达到足够的时间，并因此危害了人类的舒适、健康和福利或环境的现象。随着我国改革开放带来的经济迅速发展，城市人口规模扩大、数量增加，也带来了一些空气污染，为了强化科技支撑大气污染防治工作，科技部与环境保护部编制了《大气污染防治先进技术汇编》[1]，其中提到要对电站锅炉烟气排放、工业锅炉及煤窑锅炉排放、典型有毒有害工业废气净化、机动车尾气排放、居室及公共场所典型空气污染物净化、无组织排放源等进行控制，给予大气复合污染检测模拟与决策支持、扩大清洁生产等八个领域的关键技术。近十年来，我国持续快速增长，但经济增长模式相对传统落后，对生态平衡和自然环境造成一定的破坏，空气污染的弊病日益突出，特别是日益加重的雾霾天气已经干扰到社会的出行秩序和生活质量。国家能源委员会《新能源产业振兴和发展规划》等“国家新能源发展战略”政策的出台，说明国家已经把能源环境问题上升到国家安全级别，经济发展转型、节能减排、能源利用新途径和发展新能源等方面的问题亟待解决。环保部中国环境监测总站2013年逐月发布了《城市空气质量状况报告》[2]，报告空气质量综合指数排名“空气质量综合指数最差前十名”中，包含了河北省的六个市，其中有五个市并列第一，且每月上榜，另一个市排第六，上榜十次；天津市排名第九，上榜六次，而北京市作为中国首都，排名第二十，由此可见，研究并改善京津冀地区的空气质量状况具有十分重大的意义。1.2问题重述一般认为影响空气质量的主要影响因素有、、、、、、硫化氢、碳氢化合物和烟尘等，以京津冀地区为研究对象解决以下问题：1、参考现有国标和美标，建立衡量空气质量优劣程度等级的数学模型。2、查找分析数据，得出京津冀地区主要污染物，推测主要污染参数及列出其污染源，分析影响空气质量的主要污染源的性质和种类。3、现知河北境内某一工厂废气排放烟囱高50m，主要排放物为氮氧化物。早上9点至下午3点期间的排放浓度为406.92，排放速度为1200；晚上10点-凌晨4点期间的排放浓度为1160，排放速度为5700；通过你的扩散模型求解该工厂方圆51公里分别在早上8点、中午12点、晚上9点空气污染浓度分布和空气质量等级。建立单污染源空气污染扩散模型，描述其对周围空气污染的动态影响规律。4、以汽车尾气污染源为例求解分析：北京在2015年1月15日已经连续三天发生重污染，假设从16日开始北京启动汽车单双号限行交通管制措施。由单污染源空气污染扩散模型推广到多污染源空气污染扩散模型，求解北京市二环、四环、六环路在16日早上8点、中午12点、晚上9点时空气污染浓度梯度变化及空气质量等级。5、根据建立的模型和求解结果，分析总结影响空气质量的关键参数，为京津冀地区环保部门撰写一份建议报告，给出实现“”蓝天的可行性措施和建议。二、问题分析2.1对问题一的分析李新在《建筑材料对室内空气质量的影响及其评价》[3]中，通过对室内环境质量的研究，根据质量平衡的原理，建立了室内空气质量的基本数学模型，求出各污染物的浓度并依据室内空气质量标准，结合室内空气污染特点，判断出室内空气质量的优劣与等级。此建立模型的思路便于解决本问题，因此我们采用其思路解决本问题。在本问题中中也要通过分析主要污染物、依据空气质量分指数公式原理，参考美标和国标建立衡量空气质量优劣程度的数学模型。因为空气质量指数是在各空气质量分指数中的最大值，而空气质量分指数是根据、、、、、等各项污染物的实测浓度值分别计算出的，所以只要查阅相关数据，结合各项污染物的分级浓度限值和对应的浓度数据求出各污染物的空气质量分指数，就可以确定总的空气质量指数，进而判断空气质量优劣程度。2.2对问题二的分析针对主要污染源的确定及分析其性质和种类的问题，张新在论文[3]中也提到基于室内空气质量的基本数学模型的理论研究，通过应用分析，确定主要污染物及污染源并分析出主要污染源的性质与种类。其解决的问题与本问题相近，因此我们可以采用其应用分析的思路解决主要污染源的性质与种类问题。为了求出京津冀地区主要污染源和污染参数，并分析影响空气质量的主要污染源的性质和种类，对影响空气质量的相关数据进行分析处理。京津冀地区包括北京市、天津市以及河北省的保定、廊坊、沧州、秦皇岛、唐山、承德、张家口、衡水、邢台、邯郸和石家庄。在这里，我们选取河北省的省会石家庄和北京市、天津市进行具体分析。对各项空气污染物进行分析，找出对空气质量影响较大的因素，即污染参数，并由主要污染物和具体影响因素的大小确定主要污染源。通过应用分析，再根据各污染物的性质推出主要污染源的性质和所属种类。2.3对问题三的分析针对大气污染扩散的研究，孙志宽在《高斯烟羽扩散模型再研究》论文中对以往的高斯烟羽模型进行修正，通过边界层污染气象观测、烟羽对比观测和对比计算证明后，新的高斯烟羽模型更接近实际[4]，因此本文采用此类模型。问题中要求建立单污染源空气扩散模型，描述其对周围空气污染的动态影响规律。可对单污染源中的一种情形进行具体分析，在这里，我们通过定性分析河北省境内某一工厂废气排放模型来研究污染物在空气中的扩散问题。采用修正后的高斯烟羽模型来分析工厂废气向周围扩散的过程，进而得出在不同时间段空气中污染物浓度分布情况以及对应的空气质量等级。2.4对问题四的分析刘永红等人在《一种线性扩散模型的建立及算法的实现》中，根据线源扩散的特点，在高斯扩散理论的基础上，建立了一种可用于城市开阔路、高架桥上层路的机动车尾气扩散模型[5]，此模型可用于风向与线源成任意角度时线源在大气中扩散计算上，同时计算较简便，结果偏差较小。因此在模型的建立中，我们采用高斯线源扩散模型。问题中要建立多污染源扩散模型，描述在污染源头不单一的情况下，污染源对周围空气污染的规律，可对多污染源中的一种情形具体分析，以北京市汽车尾气污染源为研究对象分析多污染源向周围环境扩散的过程。汽车在行驶时尾气的扩散现象严格地说是随机流动点源，但是在研究公路两侧空间污染物分布时，通常将车辆排放源当作车道上连续线源处理，因此我们建立高斯线源模型。考虑到相邻环路距离不一，所以我们要考虑相邻两环路上汽车尾气叠加的情形。三、模型假设结合本题的实际，为了确保模型求解的准确性和合理性，我们排除了一些未知因素的干扰，提出以下几点假设：1、实时浓度按24小时平均浓度计算出值；2、污染物的浓度在、轴上的分布是高斯分布（正态分布）的；3、污染源的源强是连续且均匀的，初始时刻云团内部的浓度、温度呈均匀分布；4、污染气体是理想气体，遵守理想气体状态方程；5、地面对污染气体起全反射作用，不发生吸收、吸附、沉降、分解，不发生任何化学反应等；6、在建立高斯扩散模型时，道路两侧不存在高层建筑。四、符号说明为了便于问题的求解，我们给出以下符号说明（未在此说明的符号已在文中作出说明）：各项污染物浓度各项污染物的空气质量分指数对应于的指数限值，常量对应于的指数限值，常量小于或等于的浓度限值，常量大于或等于的浓度限值，常量五、模型的建立与求解经过以上的分析和准备，我们将逐步建立以下数学模型，进一步阐述模型的实际建立过程。5.1问题一模型的建立与求解为了衡量空气质量优劣程度等级，应了解对应的空气质量指数。参考国标和美标，引用空气质量分指数公式，建立空气质量指数模型。在已知各分污染物浓度值的前提下，要计算，应确定各种污染物在不同空气质量水平下的浓度限值。在此，我们主要研究各项污染物24小时平均浓度限值，其中臭氧浓度限值为8小时平均浓度限值。通过查阅相关资料，美国环保局对空气污染物的限值定义如下：表1美国空气污染物限值表污染物项目浓度限值\o"可吸入颗粒物"\o"细颗粒物"单位0~500~790~500~500~1551~10080~36750~1506~10/15~40101~150368~576150~25011~14/40~65151~200577~786250~35015~17351~41565~150201~300786~1572350~42018~34416~865150~250301~5001572~2620420~60035~57866~1284250~500对照美国各项污染物的分级浓度限值和各污染物的浓度，计算出各项空气质量分指数：已知空气质量指数（）是取各项污染物的空气质量分指数（）中的最大值，因此为：同理，研究各污染物在中国空气污染物限值表(见附录)的求得的空气质量分指数和空气质量指数。根据两表空气污染物浓度限值，将空气中污染物的浓度带入公式可求得在两国制定的空气质量标准下各空气污染物的空气质量分指数，最后取它们中的最大值即为我们所要研究的在国标和美标空气质量标准下的空气质量指数。经过查阅资料可知如下表所示的美国空气质量标准和中国空气质量标准：表2美国和中国空气质量标准美国空气质量标准0-5051-100101-199200-299300-350350-500级别一级二级三级四级五级六级等级良好中等不良极不良有害危险中国空气质量标准0-5051-100101-150151-200201-300300-500级别一级二级三级四级五级六级等级优良轻度污染中度污染重度污染严重污染将所求得的空气质量指数和上表对照，即可得到在美标和国标两种不同标准下，我们所研究的空气的质量优劣程度。应当注意的是，当污染物浓度超出500时，因为在这之上的指数是不存在的，所以这种情况下，称为“爆表”或“”[6]，也就是用空气质量指数已经无法描述这极坏的空气质量。5.2问题二模型的建立与求解5.2.1模型建立为了研究京津冀地区主要污染源及污染参数，运用因子分析法，以各污染物为变量，求出变量共同度，若变量共同度大于0.4，则为主要污染物。先应对影响这些地区空气质量的污染物数据进行预处理，通过画出三市原数据（具体北京、天津和石家庄三市空气污染数据见附录）的直方图与正态分布曲线，推测原数据是否呈正态分布，具体如下：图1北京市图2天津市图3石家庄由图可知，三市的直方图的变化趋势与正态分布曲线的走势大致相同，可推测出三市的空气污染数据呈正态分布。接下来研究三市相关数据偏、峰度数，具体说明原数据呈正态分布。运用分析的出三市相关数据偏、峰度数值如下：表3三市相关数据偏、峰度北京天津石家庄偏斜度0.4270.4270.427峰度0.8330.8330.833由表可知：北京、天津、石家庄的偏斜度都为0.427，峰度都为0.833，因为0.427<1.96，0.833<1.96，所以与推测一样原数据呈正态分布，说明数据可进行标准化处理（具体数据见附录），进而可以进行因子分析。5.2.2模型求解经过数据的预处理后，以污染物、、、、、为变量，对其进行因子分析。变量共同度就是变量方差，是每个原始变量在每个共同因子的负荷量的平方和，也就是指原始变量方差中由共同因子所决定的比率。共同度表明了原始变量方差中能被共同因子解释的部分，共同度越大，变量能被因子说明的程度越高，原始变量的信息被保留的程度越大。一般认为，变量共同度大于0.4即可接受，大于0.5则比较好。下表为三市空气污染物对空气质量指数的变量共同度：表4三市空气污染物对空气质量指数的变量共同度北京天津石家庄0.8440.7520.6110.7190.1490.5720.8540.8510.6410.7660.6880.8430.8410.7480.8590.0180.0230.159根据上表将各污染物对空气质量指数的变量共同度的大小排序：对于北京市来说：>>>>>0.5>0.4>；对于天津市来说：>>>>0.5>0.4>>；对于石家庄来说：>>>>>0.5>0.4>。总的来说：>0.5、>0.5、>0.5、>0.5、<0.4，所以京津冀地区的主要污染物即污染参数为、、和。通过查阅资料和结合污染参数的变量共同度，分析出各污染物的主要污染源如下：1、的主要污染源：燃煤、机动车尾气、扬尘、二次粒子、生物质燃烧、有机物、海盐粒子等，且各来源的贡献率因所在功能区不同而异[7]。北京市中燃煤、扬尘、有机物及二次硫酸盐和硝酸盐贡献率较大[8]；天津市的热力供应、工业燃煤及居民生活源是主要贡献来源[9]；河北省的以煤为主的得能源结构是主要贡献来源[10]。2、的主要污染源：含硫能源（如煤炭、石油）的燃烧、金属的熔炼和其他含硫原材料的生产加工等过程[11]。其中，北京市的采暖源和一般工业锅炉源行业为主要贡献行业[12]；天津市从季节来看二氧化硫季节性污染特征明显，采暖期均高于非采暖期，说明煤燃烧是主要贡献来源[13]；河北省石家庄在2013年每月连续排于空气质量最差前十，且大多数原因是重工业尤其是钢铁企业[2]，由此看来，重工业是的主要贡献来源。3、的主要来源：燃料的燃烧、城市汽车尾气的排放以及工业生产过程中也会产生，此外，闪电时由于空气中电厂极强，空气分子被撕裂而导致雷电电流通过时产生大量的热，使已经呈游离状态的空气成分和结合成；北京、天津、石家庄，经济发达，城市来往车辆多，因而城市汽车尾气的排放为主要贡献来源，若为采暖期，则燃料的燃烧也为主要贡献来源；因为石家庄受重工业的影响，所以工业生产过程中产生也是主要贡献来源。4、的来源：含碳物质的不完全燃烧和汽车尾气的排放，同，北京、天津、石家庄的城市汽车尾气的排放为主要贡献来源，若为采暖期，则燃料的燃烧也为主要贡献来源。总结以上各污染物的主要污染源分析可得：大多数主要污染源都是煤炭、石油的燃烧源、排放污染气体的化学工业、重工业企业、工厂。因此可概括出主要污染源的性质和种类：按化学性质性质和大气污染物的组成和反应可划为两种性质[14]：还原型空气污染源、氧化型空气污染源。1、还原型空气污染源：指以煤、石油等为燃料的采暖源、一般企业、工业锅炉源和排放一些化学物质或具有污染性的气体的企业污染源。2、氧化型空气污染：指以石油为燃料的石油化工厂和汽车的尾气排放源。按能源性质和大气污染物的组成和反应可划分为四种种类[14]：煤炭型空气污染源、石油型空气污染源、混合型空气污染源、特殊性空气污染源。1、煤炭型空气污染源：指一次污染物是烟气、粉尘和二氧化硫，二次污染物是硫酸及其盐类所构成的气溶胶的采暖源、一般企业锅炉源。2、石油型空气污染源：指一次污染物是烯烃、二氧化氮以及烷、醇、羰基化合物等，二次污染物主要是臭氧、氢氧基、过氧化氢基等自由基以及醛、酮和的石油化工厂和汽车的尾气排放源。3、混合型空气污染源：指以煤炭和石油为主的排放一些化学物质的工业、企业污染源。4、特殊性空气污染源：指排放一些特殊气体，如氯气、金属蒸汽或硫化氢、氟化氢等气体所造成的污染的某些工矿企业排放源。5.3问题三模型的建立与求解5.3.1模型的建立为了研究单污染源空气污染扩散的问题，建立专门解决物质扩散问题的高斯烟羽模型。具体以河北省境内某一工厂废气排放为例，以该厂烟囱正下方的地面为坐标原点，风向为轴，且指定下风向为的正轴，铅直方向为轴，且向上的方向为正方向，水平垂直于风向轴（轴）的为轴，建立空间坐标系（如下图示），设空气污染物排放点距规定水平地面的高度为，则易得排放点位置坐标为。图4烟囱排放污染物浓度的空间坐标系示意图设空间任意一点的空气污染物浓度为，它代表下风向米、横向米、地面上方米处的扩散的气体浓度，单位为；在假设中已提到污染物的浓度在轴上的分布是高斯分布即正态分布，因此可以导出下风向上任意一点处泄漏气体的浓度函数为：(1)我们可以根据概率统计原理列出如下方差的表达式：(2)假设可知如下所示的源强公式：(3)将（1）式代入（2）式，经过积分得：（4）再将（1）（4）两式代入（3）式，经过积分得:（5）最后将公式（4）（5）代入（1）式，经过积分得出任意一点的污染物浓度公式：（6）上式仅考虑在无界空间连续点源扩散的高斯烟羽模型，在实际中，不仅烟羽的扩散是有界的，地面也会对排放到空气中的污染物造成一定的反射。因此接下来我们采用像源法来处理地面对污染物的反射问题，把任一点的浓度看做是两部分的和：第一部分是假设不存在地面的情况下污染物在空气中的扩散浓度；第二部分是由于地面反射作用增加的污染物浓度。简单形象地来说，空间中一点的污染物浓度是由若不存在地面时污染实源和存在地面反射时污染像源在该点叠加的浓度之和。图5像源法原理示意图结合上述公式和原理示意图我们可以求出实源和像源分别对某一点污染物浓度的贡献，具体公式如下：实源的贡献：（7）像源的贡献：（8）综合上述两种源头，空间中任一点的污染物浓度为：（9）令y=z=0,即可得到地面轴线气体浓度计算公式：（10）在上式中的指的是平均风速，单位为；指的是空气污染物的源强（即源释放速率），单位为，已知空气中污染物的源强计算公式为:（11）其中为在排放口处污染物被排出的浓度，即污染物在排放口被排出的速度。即指污染物传出口的有效原高，考虑到在实际情况中，这个高度不仅指废气排放口距离地面的高度，还应包含受污染物初始速度、初始温度、传出口直径、环境风速、环境温度及大气稳定度影响被抬升的高度，则有效源高：（12）其中，为传出口的几何架高，为污染物被抬升高度。实验表明，传出口污染物被抬升高度可以用下面公式进行近似计算：（13）这样就可以计算出污染物在传出口处的有效高度。对于(2)式中的方差表达式即（9）式中的污染物的扩散系数，为水平扩散参数，单位为；为垂直扩散参数，单位为；它们的选取与大气的稳定度有关，具体如下表所示：表5水平扩散参数、垂直扩散参数与大气的稳定度的关系大气稳定度因此要确定污染物在两方向的标准差，就要确定对应的大气稳定度。按照的分类方法，大气的稳定度可以分为6类：，其中三类表示气象条件不稳定，类表示中性气象条件，表示气象条件稳定。下表是大气稳定度的具体分类（强：10：00-14：00；中：8：00-10：00/14：00-16：00；弱：6：00-8：00/16：00-18：00）：表6大气稳定度的确定风速白天光照（6:00-18:00）夜间条件（18:00-6:00）强中弱多云无云<22-33-44-6>6由上表可知，只要已知风速，就可以知道对应的光照强度，从而确定大气稳定度的类型，得出水平扩散参数、垂直扩散参数即对应的的值，在根据已知数据计算出、，将各值带入公式（10）即可求出某一点的沿地面轴线的污染物浓度。5.3.2模型的求解在模型的建立部分我们得到如下污染物浓度随空间、时间变化方程式：所研究的河北境内的工厂工作的时间段为早上9点到下午3点和晚上10点到凌晨4点，现要计算早上8点、中午12点一级晚上9点空气污染浓度分布和空气质量等级。由于早上8点距离凌晨4点排放污染物结束时间较长，因此可认为此时空气污染浓度为0。通过查阅相关资料，取河北省石家庄2015年4月的天气状况（具体数据见附录），分析数据得出其平均速度为4.4m/s，大部分天气为晴，因此根据风速和时间点可推测出该工厂中午12点和晚上9点的相关情况，如下表所示：表7中午12点和晚上9点的风速及稳定度等级风速稳定度等级午12点4.4晚9点4.4再根据排放口处污染物被排出的浓度，污染物在排放口被排出的速度算出不同时间空气污染物的源强。将表格中的数据和代入上述公式（9），利用编程（代码见附录）得到该工厂具体在不同时刻的51公里范围内空气污染浓度分布(详细数据表格见附录)：表851公里范围内空气污染浓度分布浓度单位：白天12时晚上9时154.043737.64233.255453.9319.858271.0450005100由于该工厂的主要排放物为氮氧化物，因此在这里我们可以近似将如上求出来的空气污染物浓度数值看作是氮氧化物的相关浓度数值，根据问题一中判断空气质量等级的空气质量指数原理：并结合问题一中中美污染物限值表中对的限值，即可求出在这51公里范围内的空气质量等级,如下表：表9随着范围变化的空气质量等级统计表美标国标空气质量等级中午12点晚上9点中午12点晚上9点0-1公里中等危险良严重污染1-2公里良好危险优严重污染2-3公里良好极不良优重度污染3-4公里良好不良优中度污染4-5公里良好中等优良5-51公里良好良好优优5.4问题四模型的建立与求解5.4.1模型的建立研究多污染源空气污染扩散问题，建立高斯线源模型。具体以北京市二、四、六环路上汽车尾气的排放为例。通过查阅资料了解到这三环的线路长和平均环间距，示意图如下：图6北京市二环、四环、六环示意图由图可知：北京市二环线路长约33.3公里，四环线路长约65.1公里，六环线路长约186.3公里，二环和四环间平均距离约为4.9公里，四环和六环间平均距离为16.22公里。由第三问的分析我们了解到随着距离的增加，污染源对远距离点的污染强度逐渐变弱，且当距离达到一定值时，所研究的点就几乎不受污染。由于二环和四环之间的平均距离为4.9km<5km，因此这两个环路汽车尾气叠加的情况不可忽略，而了解到四环和六环平均距离达16.22>6km公里，所以四环对六环可看作没有影响，将六环线路看作独立线源进行分析。在较长街道或公路上行驶的车辆密度，足以在道路两侧形成连续稳定的浓度场线源，称其为无限长线源，而在街道上行驶的车辆只能在街道两侧形成断续稳定浓度的线源，称其为有限长线源。对于汽车排放的污染，我们结合具体情况进行分析，示意图如下：图7汽车排放示意图由图我们可以直观地发现：由于线源排放路径相对固定，具有一定的方向性，若取平均风向为轴，则线源与平均风向未必相同，所以线源的情况较为复杂，应当考虑想源于风向夹角以及线源长度问题。1、风向和线源夹角时，结合第三问结论易得无限长连续线源下风向地面浓度分布为：（14）2、风向和线源夹角时，以上模式不能应用。当风向和线源的夹角垂直时，，此时：（15）考虑北京市二环、四环和六环线路长，我们将这三条环路看作有限线源来考虑。对于有限长的线源，线源末端引起的“边缘效应”将对污染物的浓度分布有很大的影响，随着污染物接收点和线源之间距离的增加，“边缘效应”将在横风向距离的更远处起作用，因此我们在估算有限长污染源形成的浓度分布时，将“边缘效应”考虑进去。对于横风向的有限长线源，以污染物接收点的平均风向为轴。若线源的范围在范围内变化，且，则有限长线源地面浓度分布为：（16）其中的人为取定汽车重心高度，值为0.5米，可以通过2015年1月16日6m/s的平均风速和问题三模型建立中的光强参照表5算出，为当日的风速。通过查阅资料，我们发现汽车尾气中主要成分仍然是氮氧化物，具体汽车尾气成分比例如下所示：表10汽车尾气成分比例统计表公路类型平均车速（km/h）kg/(1000车.km)碳氢化合物一氧化碳氮氧化物大马路160.0060.0993.6马路290.0040.059公路390.0040.048高速公路720.0020.028因此我们可以参照问题一中美国和中国空气质量标准来考虑二环、四环和六环路段上汽车尾气扩散导致的空气污染浓度梯度变化及空气质量等级的问题。要了解汽车尾气排放氮氧化物的情况，就需了解汽车尾气的源强，根据上表，左边一列为车速，记为，右边一列为每1000列车行驶1千米排放的污染物的量，记为，则每1000列车的源强为：（17）下表为北京市二环、四环和六环每日不同时刻平均车流量统计（具体数据见附录）：表11北京市二环、四环和六环每日不同时刻平均车流量统计表北京市二环四环六环7:30-8:3013684.2117789.4720526.328:30-9:305730.9947450.2928596.4919:30-10:307017.5449122.80710526.3218:30-19:3014035.0918245.6121052.6319:30-20:3011653.5214653.7117098.23根据上表我们可以知道北京市二环、四环和六环路在16日早上8点、中午12点和晚上9点的平均车流量，进而可以计算出这三条环路上的汽车在这三个时间段的总源强。接下来可以考虑二环路和四环路的污染物叠加情况：对于二环路来说，它的污染物浓度由其本身路上汽车排放的污染物和四环路在二环路上的分污染叠加而来；对于四环路来说，它的污染物由本身路上汽车排放的污染物和二环路在四环路上的分污染叠加而来。5.4.2模型的求解我们通过分析求解距离二环路、四环路和六环路各0-5公里范围内在白天8时、白天12时和晚上9时三个时间点的各处的污染物浓度，利用高斯线源模型：其中的值为0.5米，可以通过2015年1月16日6m/s的平均风速和问题三模型建立中的光强参照表5算出；为当日的风速；由已知的车速和每1000列车行驶1千米排放的污染物的量，可知每1000列车的源强。利用编程（代码见附录）得出三个时间点的各处的污染物浓度梯度的变化：表11三个时间点的各处的污染物浓度梯度的变化白天8时白天12时晚上9时距离2环4环6环2环4环6环2环4环6环0-16728731008620806930637939251-273094910955827568725917438682-371893410784175426254245336223-44525886791061381591081361594-5831081由上表我们很清楚地了解到在不考虑二环和四环污染物叠加的情况下，各环路1-5范围内污染物浓度的相关情况，在模型建立部分已了解到在实际问题中对于二环路来说，它的污染物浓度由其本身路上汽车排放的污染物和四环路在二环路上的分污染叠加而来；对于四环路来说，它的污染物由本身路上汽车排放的污染物和二环路在四环路上的分污染叠加而来。叠加后的三个环路段的浓度以及对应的在美标和国标下空气质量等级如下表所示：表12三个环路段的浓度以及对应的在美标、国标下空气质量等级2环空气质量等级4环空气质量等级6环空气质量等级美标国标美标国标美标国标白8时683危险严重污染882危险严重污染1008危险严重污染12时620危险严重污染806危险严重污染930危险严重污染晚9时630危险严重污染793危险严重污染925危险严重污染5.5关于京津冀地区实现“”蓝天的建议报告“”蓝天，是2014年一个新生的网络词汇，它是用来形容会议期间北京蓝蓝的天空。北京这几年深受雾霾天气的困扰，在会议期间，北京人民难得享受了一回“特供”蓝天，舆论将其调侃为“”蓝，国家主席习近平也呼吁，要努力把“”蓝天保持下去。自那之后，“”蓝慢慢升级为环境友好、人类与自然和谐的代名词。本文通过对空气污染问题的研究（包括单污染源和多污染源），发现影响空气质量的关键参数是、、和。因此我们针对这些参数的来源以及目前相关机构在环境保护方面的不足给出相应的建议。影响空气质量的主要污染源为燃料的燃烧、机动车尾气的排放以及扬尘等。对于机动车尾气排放及扬尘产生的部分，环保部门可以根据机动车尾气污染状况和功能区达标要求在扬尘污染，其中包括工厂烟尘、沙尘暴、建筑修建、道路修建方面进行重点治理，在汽车尾气方面进行治理和清尾，具体采取以下措施：1、鼓励公共交通工具，使用清洁能源；2、限制机动车车型；3、限制机动车行驶时间、区域、路段；4、禁止生产、销售、使用超过国家规定标准的高污染燃料；5、禁止新建、扩建燃煤供热装置。对于燃料燃烧的部分，环保部门可以通过鼓励工厂提高技术水平，提高燃料利用率以及适当更改某些工厂运作的方式等措施减轻燃料燃烧造成的污染，具体措施如下：1、限制工业燃料的燃烧废物排放量，根据具体情况，可以关闭或者整改相应的工厂；2、环保部门投资革新汽车技术和燃油技术；针对目前环境保护机构在政策制定方面的不足，提出如下建议：1、环境保护行政主管部门对行驶中的机动车污染物排放状况，应当会同公安交通管理部门进行监督抽测。在机动车停放地，对在用机动车的污染物排放状况进行监督抽测，经检测不符合污染物排放标准的，应当立即书面通知公安,且交通管理部门责令其停止运行；2、应在交通路线两旁栽植吸收各种有害气体的植被；3、应在可行条件下适当减少城市建筑和生物质燃烧；4、应增加城市森林覆盖面积。六、模型的检验6.1因子分析模型的检验针对问题二，运用因子分析法，利用得出标准数据的变量共同度，求出主要污染物。首先要对数据进行预处理，判断原数据是否呈正态分布，若呈正态分布，则可以进行因子分析；否则，不能进行因子分析。我们是先用得出原数据的直方图和正态分布曲线，推测出它可能呈正态分布，再求出原数据的偏斜度与峰度，得出其与推测相同即原数据呈正态分布，可进行因子分析。为了确定所建模型的正确性，我们对因子分析模型进行检验。主要利用求出检定值来检验原数据是呈正态分布，可进行因子分析。对三市进行检定的值如下表所示：表13三市检定的值显著性北京市0.6840天津市0.5940石家庄0.6600由上表可知：0.684>0.5、0.594>0.5、0.660>0.5即三市的的值都大于0.5，说明原数据可进行因子分析；又因为显著性值都为0，说明原数据中没有太大的数据值出现，因此做因子分析是合理的。从而也验证了因子分析模型建立的准确性与合理性。七、模型的评价与改进7.1模型的评价7.1.1模型的优点1、因子分析把多个指标中的大部分信息浓缩到几个比较少的指标中，便于问题的求解；2、高斯烟羽模型和高斯线源扩散模型能较好、较具体解释物质在单污染源和多污染源两种类型下的污染物扩散问题；3、使用高斯修正模型，极大提高了大气污染物扩散精度和大气环境控制水平。7.1.2模型的缺点1、因子分析中浓缩后的因子意义不能完全确定，且会有一些信息没有被获取；2、高斯扩散模型只适用于道路周围平坦的高速公路以及建筑物少、低矮的道路。7.2模型的改进1、实际生活中污染物的种类不同，其扩散方式及速度就会不同，应对污染物进行分类讨论；2、实际生活中污染物的扩散过程与周围阻挡物、大气的水平运动、垂直运动、水平面存在的气压梯度力、地球自身的自转以及地面摩擦力均有关，应对不同时刻的天气情况和大气运动详细讨论。八、模型的推广因子分析模型的是主成分分析的推广。它利用了降维的思想，把一些具有错综复杂关系的变量归结为少数几个综合因子的一种多变量的统计分析方法，其目的是用有限个不可观测的隐变量来解释原始变量之间的相关关系。与主成分分析相比，因子分析更倾向于描述原始变量之间的相关关系。早期的因子分析主要应用于心理学与教育学方面的问题。近年来，随着电子计算机的高速发展，因子分析理论广泛的应用与医学、气象、地质、经济、社会、考古、体育科学、生物等各个领域。高斯模型又分为高斯烟团模型和高斯烟羽模型。高斯模型主要解决大气污染物泄露问题。大气污染物泄露可分为瞬时泄露和连续泄露，瞬时泄露指的是污染物泄露的时间相对于污染物扩散的时间较短如图发泄泄露等情形，采用高斯烟团模型进行模拟；连续泄露则是指污染物泄放的时间较长的情形，采用高斯烟羽模型进行模拟。高斯模型一般适用于非重云气体（包括轻气云和中性气云气体），且风俗均匀稳定。高斯扩散模式可分为四种：连续点源的扩散、连续线源的扩散、连续面源的扩散、连续体源的扩散。连续点源一般指排放大量污染物的烟囱、放散管、通风口等。若排放口安置在地面则称为地面点源，处于高空位置的称为高架点源。连续线源是污染物沿一水平方向连续排放，如汽车行驶在开阔的公路上。连续面源的扩散是众多的污染源在一地区内排放，如城市中小工厂、生活锅炉、居民炉灶等数量众多、分布面广、排放高度低的污染源。连续体源扩散是污染物自低矮烟道或房屋的排放系统排放，常会受到附近建筑物的空气动力学效应的影响，或者在一些工厂的无组织排放和泄露排放的情形。九、参考文献[1]中国环境科学学会，大气污染防治先进技术汇编[J]，中国环境科学，34（4）：875-876，2014；[2]庄庆鸿陈晓孙悦杨雪，空气最差城市全国人大代表12%来自高污染行业[J]，中国青年报，14（4）：7-8，2014；[3]李新，建筑材料对室内空气质量的影响及其评价[D]，中国优秀硕士学位论文全文数据库，2004；[4]孙志宽，高斯烟羽扩散模型再研究[J]，环境与可持续发展，2013（5）；[5]刘永红余志田丰蔡铭，一种线性扩散模型的建立及算法的实现[J]，环境科学与技术，29（10）：13-21，2006；[6]王宗爽武婷车飞等，中外环境空气质量标准比较[J]，环境科学研究，23（3）：256-257，2013；[7]陈璐赵东风薛建良等，环境空气中PM2.5的分布特征及来源解析[J]，现代化工，33（9）：5-6，2013；[8]朱先磊张远航曾立民王玮，北京市大气细颗粒物PM2.5的来源研究[J]，环境科学研究，18（5）：1-2，2005；[9]陈璐周阳王伟等，天津市PM2.5排放源构成分析及防治措施评估[J]，城市环境与城市生态，27（3）：26-27，2014；[10]张荣李波米慧欣文丹王磊，河北省PM2.5污染特征[J]，统计分析及对策研究，2014（19）；[11]刘玉香，SO2的危害及其流行病学与毒理学研究[J]，生态毒理学报，2（2）：225-231，2007；[12]黄青程水源陈东升等，北京市一次SO2污染过程来源分析[J]，环境科学与技术，33（1）：90-92，2010；[13]刘际芳李洁，天津市环境空气中污染物的种类分析[J]，天津科技，33（3）：60-61，2006；[14]唐晓慧，大气污染的主要类型及防治技术探讨[J]，科技与企业，2014(21)。附录1问题三模型求解中求不同时刻方圆51公里的污染物浓度的代码：1、中午12点：forx=1:1:51y=(134.28/(3.14*4.4*(0.22*x/((1+0.0004*x)^0.5)))*0.2*x)*exp(-2500/(2*0.04*x^2));disp(y)end2、晚上9点：forx=1:1:51y=(1832.8/(3.14*4.4*(0.22*x/((1+0.0004*x)^0.5)))*0.2*x)*exp(-2500/(2*0.04*x^2));disp(y)end问题四模型求解中2、4、6环在不同时段的汽车排放污染气体浓度的代码：1、白天8时，2环forx=1:1:5f1=(13684*39)/(15*0.2*x);f2=exp(-0.25/(0.04*x^2));f3=0.4*exp(-0.5)-(1/85.75)*exp(-588.245);y=f1*f2*f3;disp(y)end2、白天8时4环forx=1:1:5f1=(17789*39)/(15*0.2*x);f2=exp(-0.25/(0.04*x^2));f3=0.4*exp(-0.5)-(1/85.75)*exp(-588.245);y=f1*f2*f3;disp(y)end3、白天8时6环forx=1:1:5f1=(20526*39)/(15*0.2*x);f2=exp(-0.25/(0.04*x^2));f3=0.4*exp(-0.5)-(1/85.75)*exp(-588.245);y=f1*f2*f3;disp(y)end4、白天12时2环forx=1:1:5f1=(11461*39)/(15*0.14*x*(1+0.0003*x)^(-0.5));f2=exp(-0.25/(0.0196*(x^2)*(1+0.0003*x)^(-1)));f3=0.4*exp(-0.5)-(1/85.75)*exp(-2184.605);y=f1*f2*f3;disp(y)end5、白天12时4环forx=1:1:5f1=(14900*39)/(15*0.14*x*(1+0.0003*x)^(-0.5));f2=exp(-0.25/(0.0196*(x^2)*(1+0.0003*x)^(-1)));f3=0.4*exp(-0.5)-(1/85.75)*exp(-2184.605);y=f1*f2*f3;disp(y)end6、白天12时6环forx=1:1:5f1=(17192*39)/(15*0.14*x*(1+0.0003*x)^(-0.5));f2=exp(-0.25/(0.0196*(x^2)*(1+0.0003*x)^(-1)));f3=0.4*exp(-0.5)-(1/85.75)*exp(-2184.605);y=f1*f2*f3;disp(y)end7、晚上9时2环forx=1:1:5f1=(11653*39)/(15*0.14*x*(1+0.0003*x)^(-0.5));f2=exp(-0.25/(0.0196*(x^2)*(1+0.0003*x)^(-1)));f3=0.4*exp(-0.5)-(1/85.75)*exp(-2184.605);y=f1*f2*f3;disp(y)end8、晚上9时4环forx=1:1:5f1=(14653*39)/(15*0.14*x*(1+0.0003*x)^(-0.5));f2=exp(-0.25/(0.0196*(x^2)*(1+0.0003*x)^(-1)));f3=0.4*exp(-0.5)-(1/85.75)*exp(-2184.605);y=f1*f2*f3;disp(y)end9、晚上9时6环forx=1:1:5f1=(17098*39)/(15*0.14*x*(1+0.0003*x)^(-0.5));f2=exp(-0.25/(0.0196*(x^2)*(1+0.0003*x)^(-1)));f3=0.4*exp(-0.5)-(1/85.75)*exp(-2184.605);y=f1*f2*f3;disp(y)end附录2问题一的模型建立与求解中中国空气污染物限值表表1中国空气污染物限值表污染物项目浓度限制\o"可吸入颗粒物"\o"细颗粒物"单位0~500~500~400~500~20~1000~3551~10050~15040~8050~1502~4100~16035~75101~150150~47580~180150~2504~14160~21575~115151~200475~800180~280250~35014~24215~265115~150201~300800~1600280~565350~42024~36265~800150~250301~5001600~2100565~750420~50036~48800~250~350401~5002100~2620750~940500~60048~60-3350~500问题二的模型建立中北京、天津和石家庄三市空气污染物浓度的原数据：表2北京市2015年4月份的空气质量状况日期2015/4/158.62148.8362015/4/0.69925.865.52015/4/30.56826.988.52015/4/462.33636.190.92015/4/5536229.5662015/4/620.58913.570.12015/4/7234822.7692015/4/897.2123.712.30.84848.551.92015/4/9161.7189.916.61.19258.957.32015/4/10150.1185.591.07161.9602015/4/11141.476.516.31.74746.473.82015/4/1211.926.120.2641276.42015/4/1326.47422.396.12015/4/140.70554.466.12015/4/15108.728614.10.9397054.22015/4/1625.713230.23513.486.52015/4/1767.3173.1160.6754.854.32015/4/18137.7272.526.91.49788.748.42015/4/1959.80829.976.82015/4/2066.68143.262.32015/4/2162.7129.180.64243.868.32015/4/2243.2148.8702015/4/2354.2107.710.60.61849.160.82015/4/2451.92647.475.62015/4/2571.6326499.12015/4/2698.3181.819.80.96952.8122.62015/4/2781.2133.812.20.89945.9103.72015/4/2843.5685.10.71138.885.52015/4/29115135.610.81.32448100.22015/4/30114.26743.4123.4表3石家庄2015年4月份的空气质量状况日期2015/4/146.374.616.21.09932.721.62015/4/248.5575.50.95626.524.42015/4/342.23327.233.32015/4/468.45142.356.62015/4/558.380.715.20.97533.957.22015/4/6254411.60.79319.560.22015/4/745.17118.50.97332.852.82015/4/867.8638.747.32015/4/976.794.620.70.98138.550.62015/4/10105.36343.457.32015/4/11109.411624.41.30136.551.52015/4/1227.538.680.83617.260.52015/4/1342.4618.30.92729.578.22015/4/14105.7134.427.81.49850.662.52015/4/1591188.818.10.87528.835.92015/4/1628.85610.245.42015/4/1754.49826.232.52015/4/1881.4192.818.71.01345.916.72015/4/1944.810712.10.8942818.62015/4/2067.8126.312.30.85826292015/4/2163.712613.10.71119.6372015/4/2254.29931.5302015/4/2370.5136.2160.80825.3382015/4/2450.5146.4100.59829272015/4/2575.8166.920.20.85929.347.82015/4/2659.5106.912.50.7215.662.72015/4/2778.1512293.22015/4/2854.48110.30.73823.874.62015/4/2976.9103110.80330.358.82015/4/30111.2144.716.80.85128.7103.7表4石家庄2015年4月的天气质量状况日期2015/4/169.990.314.70.80231.339.92015/4/252.371.519.30.73528.160.12015/4/3110.8151.926.61.09650.560.12015/4/492.3118.523.70.5238.570.42015/4/570.7135.8240.51439.751.82015/4/629.955.313.20.3061956.72015/4/777.6129.538.30.61643.8522015/4/8100.70960.632.82015/4/9148.33657.844.42015/4/1082.3142.836.50.96258.249.22015/4/1196.74656.4522015/4/1233.40312.876.32015/4/1330.450.628.60.529.172.52015/4/1443.90257.166.62015/4/15101.926263.554.42015/4/1653.4821982.42015/4/1779.73154.955.72015/4/18130.4304.7561.28571.824.82015/4/1982.5138.426.30.99439.256.62015/4/2079.1145.537.21.195148.42015/4/2163.3134.454.31.0495351.32015/4/2234.3100.443.60.78750.559.12015/4/2347.46956702015/4/2448.4127.837.70.77538.483.22015/4/2583.2186.863.61.29579.763.42015/4/2687.28867.3612015/4/2769.9139.361.90.9953.866.62015/4/2895153.647.30.9842.570.92015/4/29113.119541.61.0595690.12015/4/30105.4173.1450.99745.3103.6问题二中因为原数据呈正态分布，所以原数据可进行因子分析。对三市的标准化处理数据为：表5北京市各污染物的标准化数据-0.29649-0.59528-0.73343-0.124380.33062-1.87996-0.86587-1.49363-0.883-0.18647-0.97819-0.46915-1.03306-1.16655-0.883-0.55621-0.915590.6308-0.20684-0.64502-0.70019-0.36429-0.392070.74558-0.43217-0.2893-0.99933-0.85539-0.76764-0.44524-1.21962-1.28562-1.11566-1.34368-1.67811-0.24916-1.15905-1.18011-0.883-1.17715-1.15459-0.301770.638770.090540.562810.234080.31354-1.119562.201561.088371.277411.2050.90535-0.861311.92051.022050.01440.863481.07606-0.732191.7097-0.620911.227562.771460.19404-0.07221-1.42799-1.38058-1.1489-1.41424-1.763470.05213-1.07667-1.18916-1.01595-1.10377-1.177350.99426-0.40552-0.51087-0.60048-0.169540.64928-0.440460.91742.536880.861950.490921.53699-1.00956-1.093630.21564-0.98271-1.49609-1.68380.53515-0.085690.835141.1777-0.268320.67204-1.004781.620052.33342.989122.065852.6011-1.28694-0.267410.22318-0.55063-0.44331-0.744880.07126-0.10265-0.07979-0.367820.044970.01195-0.62219-0.197150.17193-0.15178-0.347350.04609-0.33525-0.66962-0.26368-0.11855-0.437670.33062-0.25395-0.40309-0.150630.2803-0.415090.34769-0.69393-0.458820.56534-0.26811-0.392510.250950.013870.01850.610560.579430.188921.195561.137730.665420.966281.809210.575590.558232.26160.25110.242780.54620.378020.165591.35772-0.66235-0.74903-0.63372-0.1526-0.238430.487331.070050.269910.313541.577570.285091.190341.050660.428170.396630.852190.023332.29986表6处理后的天津市污染物的标准化数据-0.76182-0.837680.176170.939080.33577-1.28934-0.66915-1.26986-1.636350.25533-0.34678-1.15676-0.93451-1.2895-0.80632-0.33279-0.26972-0.735350.16904-0.579840.92151.665871.392630.36791-0.25637-0.687890.006780.346180.467880.39632-1.65898-1.58908-0.60304-0.52405-1.117410.53837-0.81236-0.926080.565780.336620.346780.187980.14377-0.481621.531320.274460.99631-0.072450.51864-0.346560.938440.374870.974290.083811.723290.282061.243351.24511.513730.401061.895980.178931.56521.904940.754110.12642-1.55368-1.72168-1.21286-0.31845-1.370620.55258-0.92609-1.17163-1.162040.11667-0.016511.390671.740140.630752.141142.84692.306380.647281.120971.966580.49802-0.13197-0.09358-0.61224-1.498930.75108-1.92432-2.13541-2.14124-0.16241-0.420641.19062-0.09486-0.97829-0.37981-0.773230.716612.064810.599650.527881.78896-1.52136-0.825-0.04207-0.51835-0.04112-0.18165-1.431390.143770.43185-0.48447-0.21325-0.40183-0.93895-0.028920.42449-0.34895-0.91613-1.1064-0.56015-0.42907-0.52828-0.78938-0.495360.20367-0.89160.25750.674950.14229-0.45233-0.47889-0.5128-0.584910.92542-0.87407-1.45644-0.07156-1.033650.480731.428810.85375-0.20847-0.03853-0.04877-0.20583-0.04453-0.45059-0.8731-1.546760.656750.577610.503060.97232-0.24672-0.842192.10092-0.42064-0.68052-0.82326-0.78703-0.644021.220210.52707-0.14029-0.70468-0.476240.071560.472081.97180.883680.27781-0.24672-0.104592.5981表7处理后的石家庄污染物的标准化数据-0.24125-0.9857-1.50502-0.32625-1.02195-1.27141-0.82987-1.29577-1.22567-0.54338-1.22389-0.047071.126630.03029-0.782370.626570.18975-0.047070.50791-0.52059-0.95847-1.24017-0.567560.57722-0.21449-0.23525-0.94026-1.25961-0.49183-0.55014-1.57902-1.56297-1.59611-1.93371-1.79819-0.253150.01628-0.33916-0.07186-0.92905-0.23308-0.538020.788840.002250.596140.992780.82715-1.701752.380790.66199-0.260111.728460.65044-0.998660.17346-0.1198-0.181170.192290.67569-0.707730.655060.429430.796541.112690.56209-0.53802-1.46196-0.90323-1.80865-1.61935-2.189460.93482-1.5623-1.64048-0.66091-1.30499-1.160790.7045-1.1108-0.937870.37752-0.326250.606270.34690.828972.309671.203410.840461.01017-0.39256-0.793081.17328-1.47465-2.01149-1.798191.304550.086511.616950.899770.091820.46743-0.313761.782142.550481.003011.239091.53397-2.186640.18015-0.19237-0.800580.296-0.52338-0.259210.06644-0.07526-0.138660.931210.2213-0.75622-0.46198-0.258340.899770.474250.34752-0.58045-1.43186-0.819120.24999-0.374860.18975-0.10769-0.99374-0.6080.984790.214980.536850.55297-0.9603-0.3672-0.1083-0.41375-0.573871.353040.203560.605911.464531.27152.032530.152940.337340.478911.73780.924721.249980.00748-0.24125-0.177521.36130.283030.398010.34690.598210.058330.474680.25063-0.315120.607521.203550.741160.128540.506650.536851.771250.946030.379950.335010.30572-0.138422.5895问题二模型求解中：将表格中的数据和代入上述公式（9），即可得到该工厂具体在不同时刻的51公里范围内空气污染浓度分布：表851公里范围内空气污染浓度分布浓度单位：白天12时晚上9时154.043737.64233.255453.9319.858271.04411.478156.6756.404187.4163.438346.9371.770324.16288.71E-0111.88494.07E-015.5593101.80E-012.4618117.52E-021.0263122.93E-024.00E-01131.06E-021.45E-01143.55E-034.85E-02151.09E-031.48E-02163.00E-044.10E-03177.42E-051.01E-03181.62E-052.21E-04193.07E-064.18E-05204.96E-076.77E-06216.71E-089.15E-07227.40E-091.01E-07236.47E-108.84E-09244.34E-115.93E-10252.15E-122.93E-11267.45E-141.02E-12271.71E-152.34E-14282.43E-173.32E-16291.96E-192.68E-18308.08E-221.10E-20311.49E-242.03E-23321.04E-271.43E-26332.25E-313.08E-30341.15E-351.57E-34359.70E-411.32E-393