预测区域不同方位、观测点PM_2.5背景浓度计算方法

预测区域不同方位、观测点PM_2.5背景浓度计算方法
【专利摘要】一种预测区域不同方位PM2.5背景浓度的计算方法，包括：在预测区域的不同方位选取多个气象观测站点及相邻的PM2.5观测点；对气象观测数据及PM2.5值观测数据进行处理，剔除极值数据；采用裁剪平均的计算方法排除高、低监测极值；用20％裁剪平均求得数理统计法的背景浓度；基于情景分析法得到数值模拟的背景浓度；依据数理统计法的背景浓度与数值模拟的背景浓度，得到不同方位PM2.5背景浓度的空间分布。以及一种预测区域不同观测点PM2.5背景浓度的计算方法。本方法综合采用数理统计、物理识别、数值模拟方法，能较准确地得出模拟背景值浓度，能大幅降低背景浓度分析的不确定性，为管理部门提供多角度决策支持信息。
【专利说明】
预测区域不同方位、观测点PM2.5背景浓度计算方法
技术领域
[0001] 本发明设及大气环境PM2.5浓度分布技术领域，更具体地设及一种预测区域不同方 位、观测点PM2.5背景浓度的计算方法。
【背景技术】
[0002] W往的研究关于某一地区环境空气中污染物背景值的确立一般依靠设立的环境 背景站获得相应的监测数据，然而随着国家城市化进程的发展，连绵成片的城市群排放的 污染物在一定气象条件下被远距离输送，必然导致背景区域污染物浓度的抬升。
[0003] 运些研究PM2.5背景值方法多为数理统计、物理识别的方法，如何合理的剔除本地 污染源贡献，同时对区域浓度水平给予客观评价，反映出相对的背景浓度水平，是关键点同 时也是难点所在。

【发明内容】

[0004] 有鉴于此，本发明的目的在于提供一种预测区域不同方位PM2.5背景浓度的计算方 法。
[0005] 为了达到上述目的，作为本发明的一个方面，本发明提供了一种预测区域不同观 测点PM2.5背景浓度的计算方法，包括W下步骤：
[0006] 步骤1，在预测区域选取多个气象观测站点及相邻的PM2.5观测点；
[0007] 步骤2，对所述气象观测数据及对应的PM2.5值观测数据进行处理，剔除极值数据； [000引步骤3,采用裁剪平均的计算方法排除高、低监测极值，所述裁剪平均是样本数据 排序后去掉最低和最高两个等比例部分之后剩下的样本平均；
[0009] 步骤4,用20%裁剪平均代表统计学的结果，并对其20%裁剪平均范围内样本依据 风向时次加权平均，加权后的PM2.5值代表依据数理统计方法得出的背景浓度；
[0010] 步骤5,基于情景分析法模拟预测区域不同网格的背景浓度，分别得到所述预测区 域内及周边的所有污染源均不关闭和均关闭两种情况下的PM2.5污染分布;根据预测结果计 算所述预测区域的所有不同网格内均关闭与均不关闭两种情况下PM2.5污染分布模拟结果 的比值，即不同网格的周边源的贡献率;利用得到的贡献率乘W不同观测点位的实测值，得 到依据数值模拟的点位的背景浓度；
[0011] 步骤6,将所述依据数理统计方法得出的背景浓度与所述对应模拟网格的背景浓 度进行加权平均，得到任意所述气象观测站点或PM2.5观测点的最终PM2.5背景浓度。
[0012] 作为本发明的另一个方面，本发明还提供了一种预测区域不同方位PM2.5背景浓度 的计算方法，包括W下步骤：
[OOU]步骤1，在预测区域的不同方位选取多个气象观测站点及相邻的PM2.5观测点；
[0014] 步骤2,对所述不同方位的气象观测数据及对应的PM2.5值观测数据进行处理，剔除 极值数据；
[0015] 步骤3,采用裁剪平均的计算方法排除高、低监测极值，所述裁剪平均是样本数据 排序后去掉最低和最高两个等比例部分之后剩下的样本平均；
[0016] 步骤4,用20%裁剪平均代表统计学的结果，并对其20%裁剪平均范围内样本依据 风向时次加权平均，加权后的PM2.5值代表依据数理统计方法得出的背景浓度；
[0017] 步骤5,基于情景分析法模拟预测区域不同网格的背景浓度，分别得到所述预测区 域内及周边的所有污染源均不关闭和均关闭两种情况下的PM2.5污染分布;根据预测结果计 算所述预测区域的所有不同网格内均关闭与均不关闭两种情况下PM2.5污染分布模拟结果 的比值，即不同网格的周边源的贡献率;利用得到的贡献率乘W不同观测点位的实测值，得 到依据数值模拟的点位的背景浓度；
[0018] 步骤6,依据所述依据数理统计方法得出的背景浓度与所述对应模拟网格的背景 浓度，得到不同方位的加权系数，将步骤5中模拟的预测区域人为源关闭后的PM2.5空间浓度 分布乘W所述不同方位的加权系数，最终得到所述预测区域不同方位的PM2.5背景浓度的空 间分布。
[0019] 其中，步骤1中所述的不同方位包括刚k、北、东北、东、东南、南、西南、西八个方位 中的至少=个W上方位。
[0020] 其中，步骤2中所述剔除极值数据的步骤包括：
[0021] 步骤21，将风向设为4个象限，单小时风向跨越两个象限的数据予W剔除。
[0022] 其中，步骤2中所述剔除极值数据的步骤进一步包括：
[0023] 步骤22,在步骤21之后，结合每日08:00和20:00的探空数据资料及实况天气图做 进一步的判别，保留中低空天气形势配置中尺度过程清楚、850hPa至低空风向较为一致、风 向差《60°的数据，由此得到整体风向较为统一的气象数据时间序列。
[0024] 其中，步骤2中所述剔除极值数据的步骤进一步包括：
[0025] 步骤23,在步骤22之后，结合遥感、实测天气图资料对特殊污染时段的数据进行进 一步剔除；
[0026] 其中，所述特殊污染时段包括沙尘过程、烟花爆竹燃放、賴杆燃烧、污染回流时段。
[0027] 其中，步骤4中所述的20%裁剪平均步骤是通过对样本数据的20%裁剪平均范围 内进行样本算术平均或依据风向时次加权平均来得到最终结果，其中优选采用依据风向时 次加权平均计算。
[0028] 其中，步骤5中所述情景分析法是通过空气质量数值模式CMAQ来实现的。
[0029] 其中，步骤6中依据所述依据数理统计方法得出的背景浓度与所述对应模拟网格 的背景浓度，得到不同方位的加权系数的步骤中，所述不同方位的加权系数根据历史数据 及经验综合给出。
[0030] 其中，步骤6中依据所述依据数理统计方法得出的背景浓度与所述对应模拟网格 的背景浓度，得到不同方位的加权系数的步骤中，同一个方位不同模拟网格的加权系数相 同。
[0031 ]基于上述技术方案可知，本发明的计算方法具有如下有益效果:本发明的计算方 法采用观测数据及数值模拟技术相结合的方式，综合采用数理统计、物理识别、数值模拟的 方法，能够相对比较准确地得出模拟背景值浓度，为管理部口提供多角度的决策支持信息； 此外，采用模式情景分析方法，集成模型模拟相对值较为准确、观测点位实测浓度较为准确 的各自优点，能够大幅降低背景浓度分析的不确定性并能得出一定区域内PM2.5背景浓度空 间分布；当运用于北京市的预测分析时，研究结果及方法可为北京市大气污染防治及十= 五规划提供决策及管理支持。
【附图说明】
[0032] 图1是北京地区的监测站分布和分类情况图；
[0033] 图2是本发明方法的完整技术路线图；
[0034] 图3是北京市不同风向不同站点PM2.5背景浓度箱式分布图；
[0035] 图4是北京市PM2.5背景浓度的空间分布图。
【具体实施方式】
[0036] 为使本发明的目的、技术方案和优点更加清楚明白，W下结合具体实施例，并参照 附图，对本发明作进一步的详细说明。
[0037] 本发明中的PM2.5背景值是指未叠加待预测区域局地的人为排放污染时，大气环境 中PM2.5的浓度水平;既包括待预测区域自然环境背景，又包括一定范围内长距离输送的贡 献。为了结合实际情况阐述本发明的技术方案，下面W北京市2013年为情景年进行示例性 地阐述和说明。
[003引（1)数理统计
[0039] 首先，在北京市西北、东北、东、东南、南、西南6个不同方位选取了密云水库、八达 岭、东高村、琉璃河、永乐店、愉堡6个站点，挑选出北京市不同方位的6个气象观测站点，如 图1所示。
[0040] 其次，进行背景时段筛选，各风向设为4个象限，同一象限下，风向角度持续2个时 次W上，单小时风向跨越两个象限将予W剔除;在依据上述标准进行风向筛选后，结合每日 08:00和20:00的探空数据资料及实况天气图做进一步的判别，中低空天气形势配置中尺度 过程清楚，且85化化至低空风向较为一致，风向差《60%运样得到了整体风向较为统一的 气象数据时间序列;再结合遥感、实测天气图等资料对特殊污染(沙尘过程、烟花爆竹、賴杆 燃烧、污染回流)时段的数据也进一步剔除，最终获得全年8760个时段。
[0041] 最后，采用裁剪平均的计算方法排除较高和较低的监测极值后的统计结果，裁剪 平均是样本数据排序后去掉最低和最高两个等比例部分之后剩下的样本平均，去掉极端值 后剩下的结果更能接近稳定的值。各站各风向20%-50%裁剪平均的结果差异全部在IOii g-Hf3W内，区间比较窄，说明W80%(20%裁剪平均）的样本，已经可W代表背景值的整体 水平，可W通过进一步计算获得背景浓度的参考值。
[0042] 最终用20%裁剪平均代表统计学的结果，并对其20%裁剪平均范围内样本算术平 均和依据风向时次加权平均两种计算方法结果进行比较，加权结果略高于前者且计算方式 相对更为合理，最后采用加权结果。
[0043] (2)数值模式
[0044] 研究采用情景分析的方法模拟北京市PM2.5的背景浓度，步骤如下：（1)利用美国环 保局推荐的CMAQ(Community Multiscale Air Quality Modeling System,https:// www.cmascenter.org/)模拟一个完整年北京市及周边所有源都不关闭时的PM2.日污染分布； (2)利用CMAQ模拟一个完整年北京人为源关闭时的PM2.5污染分布，其结果代表了北京市的 背景浓度或天然源与周边源贡献相对分布；（3)计算北京市所有网格第(2)步骤的模拟结果 与第(1)步骤的模拟结果的比值，得出不同网格的周边源的贡献率；（4)利用第(3)步骤计算 的贡献率乘W不同观测点位的实测值，得到模拟的背景浓度绝对值的点位分布。
[0045] 数理统计得出的点位背景浓度与数值模拟得出的点位背景浓度进行加权平均，可 W得到最终的不同点位的最终背景浓度;而将最终的不同点位的最终背景浓度除W数值模 拟得出的点位背景浓度，可W得到不同方位的加权系数的初值，再结合历史数据和经验对 该值进行微调即得到最终的不同方位的加权系数。对于没有监测站的方位的加权系数可W 根据相邻方位加权系数及经验进行拟合得出。该系数默认同一个方位不同模拟网格的加权 系数相同，且默认在一定较长时期内该值数值不变。用模拟的北京人为源关闭时的PM2.5浓 度空间分布乘W不同方位的加权系数，得到最终的PM2.5背景浓度的空间分布。运也表明，即 使某一网格区域没有观测点位，甚至某一方位没有观测点位也能通过该方法得到其背景浓 度。
[0046] 上述整个方法的技术路线图如图2所示。
[0047] 根据上述方法对北京市2013年不同方位的PM2.5背景浓度值计算结果如下所示： [004引1.统计学结果：
[0049] 6个监测站点PM2.5北风时段背景值最低，西风时段次之，南风、东风时段明显偏高， 不同风向下背景浓度平均值分别在6.5~27.9、22.4~73.4、67.2~91.7、40.7~116.化邑? Hf3之间，表现出北京东、南方向PM2.5背景浓度较高分布特点。
[(K)加]2.数值模式结果：
[0051 ]表1不同站点PM2.5背景浓度统计
[0化2]
[0053] 模拟的北京市PM2.5年均背景浓度的空间分布与实测年均浓度空间分布相比，背景 浓度整体低于实测浓度30~40%左右;模拟的2013年北京市PM2.5背景浓度空间分布呈现出 南高北低的特征，周边区域对北京市PM2.5背景浓度空间分布影响显著。
[0054] 经过多次实际使用发现，本发明的方法极大地消除了不确定性，可W获得相对较 为准确的预测值。此外，通过上述说明和介绍，本领域技术人员清楚本发明的方法不仅适用 于北京市2013年的背景浓度的研究分析，也适用于北京市其他年份的分析，W及除了北京 市的其他任意区域的背景浓度的计算分析。
[0055] W上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详 细说明，应理解的是，W上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在 本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护 范围之内。
【主权项】
1. 一种预测区域不同观测点PM2.5背景浓度的计算方法，其特征在于，包括以下步骤： 步骤1，在预测区域选取多个气象观测站点及相邻的PM 2.5观测点； 步骤2，对所述气象观测数据及对应的PM2.5值观测数据进行处理，剔除极值数据； 步骤3,采用裁剪平均的计算方法排除高、低监测极值，所述裁剪平均是样本数据排序 后去掉最低和最高两个等比例部分之后剩下的样本平均； 步骤4,用20%裁剪平均代表统计学的结果，并对其20%裁剪平均范围内样本依据风向 时次加权平均，加权后的PM2.5值代表依据数理统计方法得出的背景浓度； 步骤5,基于情景分析法模拟预测区域不同网格的背景浓度，分别得到所述预测区域内 及周边的所有污染源均不关闭和均关闭两种情况下的PM2.5污染分布;根据预测结果计算所 述预测区域的所有不同网格内均关闭与均不关闭两种情况下PM 2.5污染分布模拟结果的比 值，即不同网格的周边源的贡献率;利用得到的贡献率乘以不同观测点位的实测值，得到依 据数值模拟的点位的背景浓度； 步骤6,将所述依据数理统计方法得出的背景浓度与所述对应模拟网格的背景浓度进 行加权平均，得到任意所述气象观测站点或PM2.5观测点的最终PM2. 5背景浓度。2. -种预测区域不同方位PM2.5背景浓度的计算方法，其特征在于，包括以下步骤： 步骤1，在预测区域的不同方位选取多个气象观测站点及相邻的PM2.5观测点； 步骤2,对所述不同方位的气象观测数据及对应的PM2.5值观测数据进行处理，剔除极值 数据； 步骤3,采用裁剪平均的计算方法排除高、低监测极值，所述裁剪平均是样本数据排序 后去掉最低和最高两个等比例部分之后剩下的样本平均； 步骤4,用20%裁剪平均代表统计学的结果，并对其20%裁剪平均范围内样本依据风向 时次加权平均，加权后的PM2.5值代表依据数理统计方法得出的背景浓度； 步骤5,基于情景分析法模拟预测区域不同网格的背景浓度，分别得到所述预测区域内 及周边的所有污染源均不关闭和均关闭两种情况下的PM2.5污染分布;根据预测结果计算所 述预测区域的所有不同网格内均关闭与均不关闭两种情况下PM 2.5污染分布模拟结果的比 值，即不同网格的周边源的贡献率;利用得到的贡献率乘以不同观测点位的实测值，得到依 据数值模拟的点位的背景浓度； 步骤6,依据所述依据数理统计方法得出的背景浓度与所述对应模拟网格的背景浓度， 得到不同方位的加权系数，将步骤5中模拟的预测区域人为源关闭后的PM2.5空间浓度分布 乘以所述不同方位的加权系数，最终得到所述预测区域不同方位的PM 2.5背景浓度的空间分 布。3. 如权利要求2所述的PM2.5背景浓度的计算方法，其特征在于，步骤1中所述的不同方 位包括西北、北、东北、东、东南、南、西南、西八个方位中的至少三个以上方位。4. 如权利要求1或2所述的PM2.5背景浓度的计算方法，其特征在于，步骤2中所述剔除极 值数据的步骤包括： 步骤21，将风向设为4个象限，单小时风向跨越两个象限的数据予以剔除。5. 如权利要求1或2所述的PM2.5背景浓度的计算方法，其特征在于，步骤2中所述剔除极 值数据的步骤进一步包括： 步骤22，在步骤21之后，结合每日08:00和20:00的探空数据资料及实况天气图做进一 步的判别，保留中低空天气形势配置中尺度过程清楚、850hPa至低空风向较为一致、风向差 <60°的数据，由此得到整体风向较为统一的气象数据时间序列。6. 如权利要求1或2所述的PM2.5背景浓度的计算方法，其特征在于，步骤2中所述剔除极 值数据的步骤进一步包括： 步骤23,在步骤22之后，结合遥感、实测天气图资料对特殊污染时段的数据进行进一步 剔除； 其中，所述特殊污染时段包括沙尘过程、烟花爆竹燃放、秸杆燃烧、污染回流时段。7. 如权利要求1或2所述的PM2.5背景浓度的计算方法，其特征在于，步骤4中所述的20% 裁剪平均步骤是通过对样本数据的20%裁剪平均范围内进行样本算术平均或依据风向时 次加权平均来得到最终结果，其中优选采用依据风向时次加权平均计算。8. 如权利要求1或2所述的PM2.5背景浓度的计算方法，其特征在于，步骤5中所述情景分 析法是通过空气质量数值模式CMAQ来实现的。9. 如权利要求2所述的PM2.5背景浓度的计算方法，其特征在于，步骤6中依据所述依据 数理统计方法得出的背景浓度与所述对应模拟网格的背景浓度，得到不同方位的加权系数 的步骤中，所述不同方位的加权系数根据历史数据及经验综合给出。10. 如权利要求9所述的PM2.5背景浓度的计算方法，其特征在于，步骤6中依据所述依据 数理统计方法得出的背景浓度与所述对应模拟网格的背景浓度，得到不同方位的加权系数 的步骤中，同一个方位不同模拟网格的加权系数相同。
【文档编号】G06Q10/04GK105956691SQ201610262180
【公开日】2016年9月21日
【申请日】2016年4月25日
【发明人】李云婷, 程念亮, 张大伟, 孙瑞雯, 董欣
【申请人】北京市环境保护监测中心