大气污染源排放清单的偏差修正方法、装置及存储介质.pdf

《大气污染源排放清单的偏差修正方法、装置及存储介质.pdf》由会员分享，可在线阅读，更多相关《大气污染源排放清单的偏差修正方法、装置及存储介质.pdf（14页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201911311924.2 (22)申请日 2019.12.18 (71)申请人 中科三清科技有限公司 地址 100193 北京市海淀区东北旺西路8号 院36号楼5层523室 (72)发明人 吴剑斌王文丁肖林鸿张稳定 钟方潜陈亚飞陈焕盛 (74)专利代理机构 北京辰权知识产权代理有限 公司 11619 代理人 刘广达 (51)Int.Cl. G01N 33/00(2006.01) (54)发明名称 大气污染源排放清单的偏差修正方法、 装置 及存储介质 (57)摘要 本申请公开了。

2、一种大气污染源排放清单的 偏差修正方法、 装置及存储介质， 包括： 将目标区 域划分为若干个地区； 基于各所述地区的预设历 史时段的大气污染物数据， 计算各所述地区的大 气污染物浓度的平均相对偏差； 获取各所述地区 之间的大气污染物浓度的相互贡献率； 基于所述 平均相对偏差和所述相互贡献率， 获取各所述地 区的大气污染源排放清单系统偏差； 利用所述大 气污染源排放清单系统偏差， 对各所述地区的原 始大气污染源排放清单进行修正。 本申请的方 法， 基于平均相对偏差和相互贡献率， 计算各地 区的大气污染源排放清单系统偏差， 对各地区的 原始大气污染源排放清单进行修正， 兼顾了相互 贡献率的影响， 。

3、修正结果更加准确， 从而提高了 模式预报准确率。 权利要求书2页 说明书9页 附图2页 CN 110988269 A 2020.04.10 CN 110988269 A 1.一种大气污染源排放清单的偏差修正方法， 其特征在于， 包括： 将目标区域划分为若干个地区； 基于各所述地区的预设历史时段的大气污染物数据， 计算各所述地区的大气污染物浓 度的平均相对偏差； 获取各所述地区之间的大气污染物浓度的相互贡献率； 基于所述平均相对偏差和所述相互贡献率， 获取各所述地区的大气污染源排放清单系 统偏差； 利用所述大气污染源排放清单系统偏差， 对各所述地区的原始大气污染源排放清单进 行修正。 2.根据权。

4、利要求1所述的方法， 其特征在于， 所述将目标区域划分为若干个地区， 包括： 根据行政区划边界将目标区域划分为若干个地区； 或者 根据经度与纬度将目标区域划分为若干个地区。 3.根据权利要求2所述的方法， 其特征在于， 所述将目标区域划分为若干个地区， 还包 括： 为每一所述地区设置一独有的标识。 4.根据权利要求1所述的方法， 其特征在于， 所述大气污染物数据包括空气质量预报数 据和监测站的大气污染物监测数据； 所述监测站位于所述地区内。 5.根据权利要求4所述的方法， 其特征在于， 所述计算各所述地区的大气污染物浓度的 平均相对偏差， 包括： 计算各监测站的大气污染物浓度的平均偏差和平均值。

5、， 根据平均偏差 和平均值计算各监测站的大气污染物浓度平均相对偏差， 计算所述各监测站的大气污染物 浓度平均相对偏差的平均值， 以所述平均值作为所述各监测站所在地区的大气污染物浓度 平均相对偏差。 6.根据权利要求1所述的方法， 其特征在于， 所述获取各所述地区之间的大气污染物浓 度的相互贡献率， 包括： 收集各所述地区在所述预设历史时段内的通过数值源解析模式模 拟的源解析数据； 所述源解析数据代表所述地区的大气污染物对所述目标区域内的其他所 述地区的大气污染物浓度变化的贡献率。 7.根据权利要求1所述的方法， 其特征在于， 所述基于所述平均相对偏差和所述相互贡 献率， 获取各所述地区的大气污。

6、染源排放清单系统偏差， 包括： 利用所有所述地区的大气污染物浓度的平均相对偏差构建平均相对偏差矩阵； 利用所有所述地区的大气污染物浓度的相互贡献率构成相互贡献率矩阵； 所述相互贡献率矩阵的逆矩阵乘以所述平均相对偏差矩阵， 得到大气污染源排放清单 系统偏差矩阵； 其中， 所述大气污染源排放清单系统偏差矩阵的各元素为各所述地区的大气污染源排 放清单系统偏差。 8.根据权利要求7所述的方法， 其特征在于， 所述利用所述大气污染源排放清单系统偏 差， 对各所述地区的原始大气污染源排放清单进行修正， 包括： 计算单位列向量与所述平均相对偏差矩阵的差； 原始大气污染源排放清单矩阵乘以所述差得到修正后的大气。

7、污染源排放清单矩阵； 所 述修正后的大气污染源排放清单矩阵中的各元素为修正后的各所述地区的原始大气污染 源排放清单。 权利要求书 1/2 页 2 CN 110988269 A 2 9.一种大气污染源排放清单的偏差修正装置， 其特征在于， 包括： 划分模块， 用于将目标区域划分为若干个地区； 计算模块， 用于基于各所述地区的预设历史时段的大气污染物数据， 计算各所述地区 的大气污染物浓度的平均相对偏差； 第一获取模块， 用于获取各所述地区之间的大气污染物浓度的相互贡献率； 第二获取模块， 用于基于所述平均相对偏差和所述相互贡献率， 获取各所述地区的大 气污染源排放清单系统偏差； 修正模块， 用于。

8、利用所述大气污染源排放清单系统偏差， 对各所述地区的原始大气污 染源排放清单进行修正。 10.一种计算机可读存储介质， 其上存储有计算机程序， 其特征在于， 该程序被处理器 执行， 以实现如权利要求1-8中任一所述的大气污染源排放清单的偏差修正方法。 权利要求书 2/2 页 3 CN 110988269 A 3 大气污染源排放清单的偏差修正方法、 装置及存储介质 技术领域 0001 本申请涉及空气质量监测技术领域， 具体涉及一种大气污染源排放清单的偏差修 正方法、 装置及存储介质。 背景技术 0002 大气污染物排放清单是空气质量预测预报和污染控制措施和政策模拟的基础， 排 放清单的确定性及其。

9、在时间、 空间以及化学成分的分配精细化程度上直接影响着预报和模 拟的效果。 目前排放清单具有很大的不确定性， 需要针对模式预报效果开展排放清单的修 正， 减少不确定性， 从而提高模式预报准确率。 0003 目前排放清单修正常用的方法采用对比验证每个地区的预报与观测数据的效果， 在不考虑气象与空气质量模式的偏差情况下， 该地区的预报与观测的偏差可以认为是该地 区排放清单的偏差所造成。 但是， 这种方法没有考虑大气是一个大气的整体性， 污染物排放 出来， 会受到气象条件影响， 输送到外地， 因此， 仅仅考虑本地的预报效果来修正本地排放 清单的偏差， 会造成本地排放清单矫枉过正， 造成更大的偏差。 。

10、0004 因此， 需要考虑大气的物理传输以及化学生成过程， 从污染物贡献率的角度来量 化本地与外地的排放清单的实际偏差。 发明内容 0005 本申请的目的是提供一种大气污染源排放清单的偏差修正方法、 装置及存储介 质。 为了对披露的实施例的一些方面有一个基本的理解， 下面给出了简单的概括。 该概括部 分不是泛泛评述， 也不是要确定关键/重要组成元素或描绘这些实施例的保护范围。 其唯一 目的是用简单的形式呈现一些概念， 以此作为后面的详细说明的序言。 0006 根据本申请实施例的一个方面， 提供一种大气污染源排放清单的偏差修正方法， 包括： 0007 将目标区域划分为若干个地区； 0008 基于。

11、各所述地区的预设历史时段的大气污染物数据， 计算各所述地区的大气污染 物浓度的平均相对偏差； 0009 获取各所述地区之间的大气污染物浓度的相互贡献率； 0010 基于所述平均相对偏差和所述相互贡献率， 获取各所述地区的大气污染源排放清 单系统偏差； 0011 利用所述大气污染源排放清单系统偏差， 对各所述地区的原始大气污染源排放清 单进行修正。 0012 进一步地， 所述将目标区域划分为若干个地区， 包括： 0013 根据行政区划边界将目标区域划分为若干个地区； 或者 0014 根据经度与纬度将目标区域划分为若干个地区。 0015 进一步地， 所述将目标区域划分为若干个地区， 还包括： 为每。

12、一所述地区设置一独 说明书 1/9 页 4 CN 110988269 A 4 有的标识。 0016 进一步地， 所述大气污染物数据包括空气质量预报数据和监测站的大气污染物监 测数据； 所述监测站位于所述地区内。 0017 进一步地， 所述计算各所述地区的大气污染物浓度的平均相对偏差， 包括： 所述计 算各所述地区的污染物浓度平均相对偏差， 包括： 计算各监测站的大气污染物浓度的平均 偏差和平均值， 根据平均偏差和平均值计算各监测站的大气污染物浓度平均相对偏差， 计 算所述各监测站的大气污染物浓度平均相对偏差的平均值， 以所述平均值作为所述各监测 站所在地区的大气污染物浓度平均相对偏差。 001。

13、8 进一步地， 所述获取各所述地区之间的大气污染物浓度的相互贡献率， 包括： 收集 各所述地区在所述预设历史时段内的通过数值源解析模式模拟的源解析数据； 所述源解析 数据代表所述地区的大气污染物对所述目标区域内的其他所述地区的大气污染物浓度变 化的贡献率。 0019 进一步地， 所述基于所述平均相对偏差和所述相互贡献率， 获取各所述地区的大 气污染源排放清单系统偏差， 包括： 0020 利用所有所述地区的大气污染物浓度的平均相对偏差构建平均相对偏差矩阵； 0021 利用所有所述地区的大气污染物浓度的相互贡献率构成相互贡献率矩阵； 0022 所述相互贡献率矩阵的逆矩阵乘以所述平均相对偏差矩阵， 。

14、得到大气污染源排放 清单系统偏差矩阵； 0023 其中， 所述大气污染源排放清单系统偏差矩阵的各元素为各所述地区的大气污染 源排放清单系统偏差。 0024 进一步地， 所述利用所述大气污染源排放清单系统偏差， 对各所述地区的原始大 气污染源排放清单进行修正， 包括： 0025 计算单位列向量与所述平均相对偏差矩阵的差； 0026 原始大气污染源排放清单矩阵乘以所述差得到修正后的大气污染源排放清单矩 阵； 所述修正后的大气污染源排放清单矩阵中的各元素为修正后的各所述地区的原始大气 污染源排放清单。 0027 根据本申请实施例的另一个方面， 提供一种大气污染源排放清单的偏差修正装 置， 包括： 0。

15、028 划分模块， 用于将目标区域划分为若干个地区； 0029 计算模块， 用于基于各所述地区的预设历史时段的大气污染物数据， 计算各所述 地区的大气污染物浓度的平均相对偏差； 0030 第一获取模块， 用于获取各所述地区之间的大气污染物浓度的相互贡献率； 0031 第二获取模块， 用于基于所述平均相对偏差和所述相互贡献率， 获取各所述地区 的大气污染源排放清单系统偏差； 0032 修正模块， 用于利用所述大气污染源排放清单系统偏差， 对各所述地区的原始大 气污染源排放清单进行修正。 0033 根据本申请实施例的另一个方面， 提供一种计算机可读存储介质， 其上存储有计 算机程序， 其特征在于，。

16、 该程序被处理器执行， 以实现所述的大气污染源排放清单的偏差修 正方法。 说明书 2/9 页 5 CN 110988269 A 5 0034 本申请实施例的其中一个方面提供的技术方案可以包括以下有益效果： 0035 本申请实施例提供的大气污染源排放清单的偏差修正方法， 基于平均相对偏差和 相互贡献率， 计算各地区的大气污染源排放清单系统偏差， 利用大气污染源排放清单系统 偏差对各地区的原始大气污染源排放清单进行修正， 兼顾了相互贡献率的影响， 修正结果 更加准确， 降低了偏差， 减少了不确定性， 从而提高了模式预报准确率。 0036 本申请的其他特征和优点将在随后的说明书中阐述， 并且， 部分。

17、地从说明书中变 得显而易见， 或者， 部分特征和优点可以从说明书中推知或毫无疑义地确定， 或者通过实施 本申请实施例了解。 本申请的目的和其他优点可通过在所写的说明书、 权利要求书、 以及附 图中所特别指出的结构来实现和获得。 附图说明 0037 为了更清楚地说明本申请实施例或现有技术中的技术方案， 下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍， 显而易见地， 下面描述中的附图仅仅是本 申请中记载的一些实施例， 对于本领域普通技术人员来讲， 在不付出创造性劳动的前提下， 还可以根据这些附图获得其他的附图。 0038 图1示出了本申请的一个实施例的大气污染源排放清单的偏差修正方法。

18、流程图； 0039 图2示出了本申请的一个实施例的大气污染源排放清单的偏差修正装置的结构框 图； 0040 图3示出了本申请的另一实施例中根据行政区划边界将目标区域划分为若干个地 区的目标区域划分示意图。 具体实施方式 0041 为了使本申请的目的、 技术方案及优点更加清楚明白， 下面结合附图和具体实施 例对本申请做进一步说明。 应当理解， 此处所描述的具体实施例仅用以解释本申请， 并不用 于限定本申请。 基于本申请中的实施例， 本领域普通技术人员在没有做出创造性劳动前提 下所获得的所有其他实施例， 都属于本申请保护的范围。 0042 本技术领域技术人员可以理解， 除非另外定义， 这里使用的所。

19、有术语(包括技术术 语和科学术语)， 具有与本申请所属领域中的普通技术人员的一般理解相同的意义。 还应该 理解的是， 诸如通用字典中定义的那些术语， 应该被理解为具有与现有技术的上下文中的 意义一致的意义， 并且除非像这里一样被特定定义， 否则不会用理想化或过于正式的含义 来解释。 0043 如图1所示， 本申请的一个实施例提供了一种大气污染源排放清单的偏差修正方 法， 包括： 0044 S1、 将目标区域划分为若干个地区。 0045 S2、 基于各所述地区的预设历史时段的大气污染物数据， 计算各所述地区的大气 污染物浓度的平均相对偏差。 0046 S3、 获取各所述地区之间的大气污染物浓度的。

20、相互贡献率。 0047 S4、 基于所述平均相对偏差和所述相互贡献率， 获取各所述地区的大气污染源排 放清单系统偏差。 说明书 3/9 页 6 CN 110988269 A 6 0048 S5、 利用所述大气污染源排放清单系统偏差， 对各所述地区的原始大气污染源排 放清单进行修正。 0049 在某些实施方式中， 步骤S1、 将目标区域划分为若干个地区， 包括： 0050 根据行政区划边界将目标区域划分为若干个地区； 或者 0051 根据经度与纬度将目标区域划分为若干个地区。 0052 在某些实施方式中， 步骤S1、 将目标区域划分为若干个地区， 还包括： 为每一所述 地区设置一独有的标识。 0。

21、053 例如， 为每一所述地区分配一独有的ID。 每一地区的标识作为数值源解析模式的 输入数据， 同时也作为大气污染源排放清单偏差空间修正的最小单位。 0054 在某些实施方式中， 所述大气污染物数据包括空气质量预报数据和监测站的大气 污染物监测数据； 所述监测站位于所述地区内。 0055 例如， 大气污染物监测数据可以为PM2.5、 PM10、 SO2、 NO2、 CO、 O3的监测站检测浓度。 空气质量预报数据可以为PM2.5、 PM10、 SO2、 NO2、 CO、 O3的预报浓度。 0056 在某些实施方式中， 所述计算各所述地区的大气污染物浓度的平均相对偏差， 包 括： 计算各监测站。

22、的大气污染物浓度的平均偏差和平均值， 根据平均偏差和平均值计算各 监测站的大气污染物浓度平均相对偏差， 计算所述各监测站的大气污染物浓度平均相对偏 差的平均值， 以所述平均值作为所述各监测站所在地区的大气污染物浓度平均相对偏差。 0057 例如， 基于收集到的空气质量数值模式的空气质量预报数据以及监测站的大气污 染物监测数据， 计算各监测站的大气污染物浓度的平均偏差和平均值， 大气污染物浓度平 均相对偏差为平均偏差和平均值的比值。 计算所述各监测站的大气污染物浓度平均相对偏 差的平均值， 所述各监测站位于同一地区内， 以该平均值作为各监测站所在地区的大气污 染物浓度平均相对偏差yi。 yi代表。

23、ID为i的地区的大气污染物浓度的平均相对偏差。 0058 在某些实施方式中， 步骤S3、 获取各所述地区之间的大气污染物浓度的相互贡献 率， 包括： 收集各所述地区在所述预设历史时段内的通过数值源解析模式模拟的源解析数 据； 所述源解析数据代表所述地区的大气污染物对所述目标区域内的其他所述地区的大气 污染物浓度变化的贡献率。 0059 具体地， 源解析为各个地区的污染源对环境污染的贡献值(分担率)， 其代表了某 个地区污染源所排放的污染物， 经过在大气中的输送扩散等物理过程以及化学反应过程， 对受体地区造成的污染浓度变化的贡献率。 0060 在某些实施方式中， 采用二次污染物在线源解析模式， 。

24、其优势在于可追溯二次污 染物的前体物污染来源， 有效获得地区污染源的贡献率。 0061 在某些实施方式中， 步骤S4、 基于所述平均相对偏差和所述相互贡献率， 获取各所 述地区的大气污染源排放清单系统偏差， 包括： 0062 S401、 利用所有所述地区的大气污染物浓度的平均相对偏差构建平均相对偏差矩 阵。 0063 例如， 假设目标区域划分为n个地区， 这样就形成了n*1的平均相对偏差矩阵y： 说明书 4/9 页 7 CN 110988269 A 7 0064 0065 其中， yi代表ID为i的地区的大气污染物浓度的平均相对偏差。 0066 S402、 利用所有所述地区的大气污染物浓度的相。

25、互贡献率构成相互贡献率矩阵。 0067 例如， 假设目标区域划分为n个地区， 由于大气的整体性， 源解析数据包括了n个源 体地区和n个受体地区， 形成了n*n的二维矩阵。 例如， 某个大气污染物源解析数据Xij为j地 区排放的大气污染物对i地区的污染贡献率： 0068其中， 1in， 1jn。 0069 S403、 所述相互贡献率矩阵的逆矩阵乘以所述平均相对偏差矩阵， 得到大气污染 源排放清单系统偏差矩阵； 0070 其中， 所述大气污染源排放清单系统偏差矩阵的各元素为各所述地区的大气污染 源排放清单系统偏差。 0071 具体地， 大气污染源排放清单系统偏差矩阵 的表达式为 0072 0073。

26、 其中， i代表ID为i的地区的大气污染源排放清单的系统相对偏差。 消减比例new old*(1- i)。 0074 计算大气污染源排放清单的系统相对偏差的最简单直接的方法是 y。 但是， 由 于大气污染物在空气中受到风场的作用， 通过平流、 扩散等物理过程输送至别的地区， 也会 发生一系列化学反应， 因此， 需要考虑各个地区的相互贡献率， 才能有效准确地量化每个地 区的大气污染源排放清单的系统相对偏差。 0075 具体计算大气污染源排放清单的系统相对偏差的公式如下： 0076 X y； 0077 X-1y； 0078 0079 其中， X为数值源解析模式计算而得的每个地区的相互贡献率； y为。

27、每个地区大气 污染物浓度的平均相对偏差， 为未知数， 是需要求解的大气污染源排放清单的系统相对偏 说明书 5/9 页 8 CN 110988269 A 8 差矩阵； 由于X为n阶方阵， 且为满秩矩阵， 可求解其可逆矩阵， 因而， 具有唯一解。 0080 在某些实施方式中， 步骤S5、 利用所述大气污染源排放清单系统偏差， 对各所述地 区的原始大气污染源排放清单进行修正， 包括： 0081 S501、 计算单位列向量与所述平均相对偏差矩阵的差； 0082 计算公式为DI- ； D代表差， I代表单位列向量。 0083 S502、 原始大气污染源排放清单矩阵乘以所述差得到修正后的大气污染源排放清 。

28、单矩阵； 所述修正后的大气污染源排放清单矩阵中的各元素为修正后的各所述地区的原始 大气污染源排放清单； 0084 计算公式为e eDe(I- )， 其中， e 代表修正后的大气污染源排放清单矩 阵， e代表原始大气污染源排放清单矩阵； 0085 ei ei(1- i)， 其中， ei 代表修正后的大气污染源排放清单矩阵的元素， ei代表 原始大气污染源排放清单矩阵的元素， i代表大气污染源排放清单的系统相对偏差矩阵的 元素， i代表地区的ID。 0086 本实施例还提供了一种大气污染源排放清单的偏差修正装置， 包括： 0087 划分模块100， 用于将目标区域划分为若干个地区； 0088 计算。

29、模块200， 用于基于各所述地区的预设历史时段的大气污染物数据， 计算各所 述地区的大气污染物浓度的平均相对偏差； 0089 第一获取模块300， 用于获取各所述地区之间的大气污染物浓度的相互贡献率； 0090 第二获取模块400， 用于基于所述平均相对偏差和所述相互贡献率， 获取各所述地 区的大气污染源排放清单系统偏差； 0091 修正模块500， 用于利用所述大气污染源排放清单系统偏差， 对各所述地区的原始 大气污染源排放清单进行修正。 0092 本实施例还提供了一种计算机可读存储介质， 其上存储有计算机程序， 该程序被 处理器执行， 以实现上述的大气污染源排放清单的偏差修正方法。 009。

30、3 本申请的另一个实施例提供了一种大气污染源排放清单的偏差修正方法， 包括以 下步骤： 0094 1)将模式区域划分为若干个地区， 源解析模式以及排放清单优化均使用相同的地 区划分； 0095 2)通过空气质量数值模式系统的预报数据， 结合站点的监测数据， 计算各个地区 的污染物浓度的平均偏差y； 0096 3)通过数值源解析模式， 获得每个地区污染物浓度的相互贡献率， 形成二维矩阵 组X； 0097 4)基于各个地区的污染物浓度的平均偏差y以及相互贡献率X， 计算每个地区的大 气污染源排放清单系统偏差 。 0098 本申请的另一个实施例提供了一种大气污染源排放清单的偏差修正方法， 包括以 下。

31、步骤： 0099 (1)将目标区域划分为若干个地区， 为每个地区设置一单独的ID。 0100 划分地区的方法可采用行政边界和经纬度两种方法。 该区域划分数据(即地区ID) 作为数值源解析模式的输入数据， 也作为大气污染源排放清单偏差空间修正的最小单位， 说明书 6/9 页 9 CN 110988269 A 9 两者保持一致。 如图3所示， 可以采用行政边界划分地区。 0101 (2)收集历史时段的空气质量数值模式预报系统业务预报结果以及监测站的大气 污染物监测数据， 收集内容PM2.5、 PM10、 SO2、 NO2、 CO、 O3的站点质量浓度。 0102 (3)基于收集到的空气质量数值模式。

32、的空气质量预报数据以及监测站的大气污染 物监测数据， 计算各站点的污染物浓度的平均偏差和平均值， 再基于各站点的地理位置计 算各ID地区的平均相对偏差yi， 假设模式区域划分为n个地区， 这样就形成了n*1的平均相 对偏差矩阵y： 0103 0104 其中， yi代表ID为i的地区的污染物浓度的平均相对偏差。 0105 (4)收集历史时段的数值源解析模式模拟的源解析数据。 0106 源解析为各个污染源对环境污染的贡献值(分担率)， 其代表了某个源体(如图1所 示的某个ID地区)所排放的污染物， 经过在大气中的输送扩散等物理过程以及化学反应过 程， 对受体地区(图一所示的某个ID地区)造成的污染。

33、浓度变化的贡献率。 0107 本实施例采用二次污染物在线源解析模式,其优势在于可追溯二次污染物的前体 物污染来源， 有效获得源区的贡献率。 假设模式区域划分为n个地区， 由于大气的整体性， 源 解析数据包括了n个源体和n个受体， 形成了n*n的二维矩阵。 例如， 某个污染物源解析数据 Xij为j地区的排放对i地区的污染贡献： 0108 0109 (5)各个地区的大气污染源排放清单的系统相对偏差 是本实施例最终要获得的 一维矩阵： 0110 0111 其中， i代表ID为i的地区的大气污染源排放清单的系统相对偏差(消减比例， new old*(1- t)。 0112 计算排放清单系统相对偏差的最。

34、简单直接的方法是 y。 但是， 由于大气污染物 在空气中受到风场的作用， 通过平流、 扩散等物理过程输送至别的地区， 也会发生一系列化 学反应， 因此， 需要考虑各个地区的相互贡献率， 才可有效的量化每个地区的大气污染源排 放清单的系统相对偏差。 0113 具体计算排放清单系统相对偏差的函数如下： 0114 X y； 说明书 7/9 页 10 CN 110988269 A 10 0115 X-1y； 0116 0117 其中， X为数值源解析模式计算而得的每个地区的相互贡献率； y为每个地区大气 污染物浓度的平均相对偏差， 为未知数， 是需要求解的大气污染源排放清单的系统相对偏 差， 由于X为。

35、n阶方阵， 且基本为满秩矩阵， 可求解其可逆矩阵， 因而， 具有唯一解。 0118 (6)得到排放清单的系统相对偏差后， 就可以对每个区域的原始排放清单进行修 正， 其公式如下： 0119 ei ei(1- i) 0120 其中ei 为第i个区域校正后的排放， ei为第i个区域原始的排放， i为第i个区域排 放清单的系统相对偏差。 0121 本申请实施例提供的大气污染源排放清单的偏差修正方法， 基于平均相对偏差和 相互贡献率， 计算各地区的大气污染源排放清单系统偏差， 利用大气污染源排放清单系统 偏差对各地区的原始大气污染源排放清单进行修正， 兼顾了相互贡献率的影响， 修正结果 更加准确， 降。

36、低了偏差， 减少了不确定性， 从而提高了模式预报准确率。 0122 需要说明的是： 0123 术语 “模块” 并非意图受限于特定物理形式。 取决于具体应用， 模块可以实现为硬 件、 固件、 软件和/或其组合。 此外， 不同的模块可以共享公共组件或甚至由相同组件实现。 不同模块之间可以存在或不存在清楚的界限。 0124 在此提供的算法和显示不与任何特定计算机、 虚拟装置或者其它设备固有相关。 各种通用装置也可以与基于在此的示教一起使用。 根据上面的描述， 构造这类装置所要求 的结构是显而易见的。 此外， 本申请也不针对任何特定编程语言。 应当明白， 可以利用各种 编程语言实现在此描述的本申请的内。

37、容， 并且上面对特定语言所做的描述是为了披露本申 请的最佳实施方式。 0125 在此处所提供的说明书中， 说明了大量具体细节。 然而， 能够理解， 本申请的实施 例可以在没有这些具体细节的情况下实践。 在一些实例中， 并未详细示出公知的方法、 结构 和技术， 以便不模糊对本说明书的理解。 0126 类似地， 应当理解， 为了精简本公开并帮助理解各个发明方面中的一个或多个， 在 上面对本申请的示例性实施例的描述中， 本申请的各个特征有时被一起分组到单个实施 例、 图、 或者对其的描述中。 然而， 并不应将该公开的方法解释成反映如下意图： 即所要求保 护的本申请要求比在每个权利要求中所明确记载的特。

38、征更多的特征。 更确切地说， 如下面 的权利要求书所反映的那样， 发明方面在于少于前面公开的单个实施例的所有特征。 因此， 遵循具体实施方式的权利要求书由此明确地并入该具体实施方式， 其中每个权利要求本身 都作为本申请的单独实施例。 0127 本领域那些技术人员可以理解， 可以对实施例中的设备中的模块进行自适应性地 改变并且把它们设置在与该实施例不同的一个或多个设备中。 可以把实施例中的模块或单 元或组件组合成一个模块或单元或组件， 以及此外可以把它们分成多个子模块或子单元或 说明书 8/9 页 11 CN 110988269 A 11 子组件。 除了这样的特征和/或过程或者单元中的至少一些是。

39、相互排斥之外， 可以采用任何 组合对本说明书(包括伴随的权利要求、 摘要和附图)中公开的所有特征以及如此公开的任 何方法或者设备的所有过程或单元进行组合。 除非另外明确陈述， 本说明书(包括伴随的权 利要求、 摘要和附图)中公开的每个特征可以由提供相同、 等同或相似目的的替代特征来代 替。 0128 此外， 本领域的技术人员能够理解， 尽管在此所述的一些实施例包括其它实施例 中所包括的某些特征而不是其它特征， 但是不同实施例的特征的组合意味着处于本申请的 范围之内并且形成不同的实施例。 例如， 在下面的权利要求书中， 所要求保护的实施例的任 意之一都可以以任意的组合方式来使用。 0129 本申。

40、请的各个部件实施例可以以硬件实现， 或者以在一个或者多个处理器上运行 的软件模块实现， 或者以它们的组合实现。 本领域的技术人员应当理解， 可以在实践中使用 微处理器或者数字信号处理器(DSP)来实现根据本申请实施例的虚拟机的创建装置中的一 些或者全部部件的一些或者全部功能。 本申请还可以实现为用于执行这里所描述的方法的 一部分或者全部的设备或者装置程序(例如， 计算机程序和计算机程序产品)。 这样的实现 本申请的程序可以存储在计算机可读介质上， 或者可以具有一个或者多个信号的形式。 这 样的信号可以从因特网网站上下载得到， 或者在载体信号上提供， 或者以任何其他形式提 供。 0130 应该注。

41、意的是上述实施例对本申请进行说明而不是对本申请进行限制， 并且本领 域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。 在权利要求中， 不应将位于括号之间的任何参考符号构造成对权利要求的限制。 单词 “包含” 不排除存在未 列在权利要求中的元件或步骤。 位于元件之前的单词 “一” 或 “一个” 不排除存在多个这样的 元件。 本申请可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实 现。 在列举了若干装置的单元权利要求中， 这些装置中的若干个可以是通过同一个硬件项 来具体体现。 单词第一、 第二、 以及第三等的使用不表示任何顺序。 可将这些单词解释为名 称。 0131。

42、 应该理解的是， 虽然附图的流程图中的各个步骤按照箭头的指示依次显示， 但是 这些步骤并不是必然按照箭头指示的顺序依次执行。 除非本文中有明确的说明， 这些步骤 的执行并没有严格的顺序限制， 其可以以其他的顺序执行。 而且， 附图的流程图中的至少一 部分步骤可以包括多个子步骤或者多个阶段， 这些子步骤或者阶段并不必然是在同一时刻 执行完成， 而是可以在不同的时刻执行， 其执行顺序也不必然是依次进行， 而是可以与其他 步骤或者其他步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。 0132 以上所述实施例仅表达了本申请的实施方式， 其描述较为具体和详细， 但并不能 因此而理解为对本申请专利范围的限制。 应当指出的是， 对于本领域的普通技术人员来说， 在不脱离本申请构思的前提下， 还可以做出若干变形和改进， 这些都属于本申请的保护范 围。 因此， 本申请的保护范围应以所附权利要求为准。 说明书 9/9 页 12 CN 110988269 A 12 图1 图2 说明书附图 1/2 页 13 CN 110988269 A 13 图3 说明书附图 2/2 页 14 CN 110988269 A 14 。