CN107273995A - 空气质量预报方法 - Google Patents
空气质量预报方法 Download PDFInfo
- Publication number
- CN107273995A CN107273995A CN201610218043.6A CN201610218043A CN107273995A CN 107273995 A CN107273995 A CN 107273995A CN 201610218043 A CN201610218043 A CN 201610218043A CN 107273995 A CN107273995 A CN 107273995A
- Authority
- CN
- China
- Prior art keywords
- data
- air quality
- historical
- current
- prediction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000003915 air pollution Methods 0.000 title abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 40
- 238000013277 forecasting method Methods 0.000 claims description 17
- 238000012952 Resampling Methods 0.000 claims description 14
- 230000007613 environmental effect Effects 0.000 claims description 13
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 claims description 12
- JCXJVPUVTGWSNB-UHFFFAOYSA-N nitrogen dioxide Inorganic materials O=[N]=O JCXJVPUVTGWSNB-UHFFFAOYSA-N 0.000 claims description 8
- MGWGWNFMUOTEHG-UHFFFAOYSA-N 4-(3,5-dimethylphenyl)-1,3-thiazol-2-amine Chemical compound CC1=CC(C)=CC(C=2N=C(N)SC=2)=C1 MGWGWNFMUOTEHG-UHFFFAOYSA-N 0.000 claims description 4
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 4
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 claims description 4
- 229910002091 carbon monoxide Inorganic materials 0.000 claims description 4
- 239000007789 gas Substances 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 7
- 238000005070 sampling Methods 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 18
- 239000003344 environmental pollutant Substances 0.000 description 17
- 231100000719 pollutant Toxicity 0.000 description 17
- 238000011160 research Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 239000000443 aerosol Substances 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000002253 acid Substances 0.000 description 1
- 238000003916 acid precipitation Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004002 angle-resolved photoelectron spectroscopy Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 239000000356 contaminant Substances 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000013618 particulate matter Substances 0.000 description 1
- 238000001303 quality assessment method Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 239000000779 smoke Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/84—Greenhouse gas [GHG] management systems
- Y02P90/845—Inventory and reporting systems for greenhouse gases [GHG]
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- Development Economics (AREA)
- General Physics & Mathematics (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种空气质量预报方法,包括以下步骤:步骤一：采集待预报地区的历史环境数据，分别得到两个或两个以上历史空气质量预测数据；步骤二：将所述两个或两个以上历史空气质量预测数据以实际观测数据为参考数据进行修正，计算出各所述空气质量预报模型在最终集合计算结果中所占权重值；步骤三：采集所述待预测地区的当前环境数据，分别得到两个或两个以上初始预测数据；步骤四：根据所述权重值对所述初始预测数据进行修正，得到最终预报数据。本发明所涉及的空气预报方法简单实用，可针对较小范围的区域进行预报，预报精度高，且计算过程简单，方便采样。
Description
技术领域
本发明涉及数值分析技术领域，适用于大气环境分析领域,特别涉及一种空气质量预报方法。
背景技术
随着工业化进程的发展，空气污染成为越来越严重的一个问题。城市区域是空气污染最常发生的地区，分布于珠江三角洲、四川盆地、长江三角洲、京津冀、东北平原地区的都市化城市和重工业城市空气污染的状况比其他地区更为严重。目前城市空气质量问题已受到普遍关注，其形成机理、影响、控制途径及预报预测系统的建立与发展已成为亟需解决的问题。
目前，城市空气质量多通过统计学预报法和数值预报法进行预报。统计学预报法是分析空气污染出现时的大尺度天气系统活动规律，分析中低空大气温度、相对湿度和风速等气象要素和污染物的统计学关系。该方法简单方便，但由于气象条件、气溶胶浓度、谱分布、化学成分的关系非常复杂，建立准确度较高的统计预报方程难度较大，且这种方法需要同期城市污染物和气象监测资料，预报成本较高。
数值预报法则是通过空气质量数值预报模式计算污染物浓度分布和变化规律。该方法虽具有完善的理论基础，但预报过程难度大且计算量。以美国环保局(EPA)为代表,从1970年至今已开发了三代空气质量模型。第三代空气质量模式,通称为Models-3,该类模式拟将所有的大气问题均考虑进模式之中,可以有效地进行较为全面的空气质量控制策略的评估。目前类似的比较常用的空气质量模式还有：EURAD模式、ADOM模式、RADM模式、UAM模式、STEM模式等，这些模式各有特色,各有其优缺点。前三种主要用于酸雨问题的研究，而后几种用于大气光化学和气溶胶二次污染的模拟。虽然各模式的动力框架结构不一,但所采取的化学反应机理主要是CBM-IV、RADM II、以及SAPRC等。
国内空气质量预报根据其应用于环境决策和环境研究的目标、污染物的排放特征和污染物在大气中的相互作用、形成和转化的特点，对模型的参数和功能作修订和改进,也发展形成了具有适应区域特点的各种尺度多模型体系。第一代模型主要是高斯烟流模型,国内第二代模式有城市尺度的空气质量预报模式(如中国科学院大气物理研究所HRDM)及区域尺度污染物欧拉输送模式(如中国科学院大气物理研究所RAQM、南京大学区域酸沉降模式RegADM等)。目前国内模式发展也进入第三代，以一个大气的概念，建构了从全球尺度、区域尺度及套网格的大气环境模式系列，有南京大学的区域大气环境模式系统(RegAEMS)、中国科学院大气物理研究所的嵌套网格空气质量预报系统(NAQPMS[10])、中国科学院大气物理研究所全球环境大气输送模式等。
上述空气质量预报模型依赖各种排放源清单能够从全球或区域尺度进行大范围的空气质量模拟预报。每个模型根据不同的大气物理化学机理从不同角度反应了大气成分的变化，但是没有一个模型能全面包括所有的大气物化反应机理。在现有技术中，虽然采用了多模式集合预报的方法，但该方法只可针对较大范围尺度的区域(如全球)或是中等范围尺度的区域(如亚洲、中国东北地区等)进行预报，预报精度较低，并且预报过程中计算量大，耗时多。
发明内容
为了解决上述问题，本发明的目的在于提供一种空气质量预报方法,对空气中一种或多种气体或颗粒物的含量进行预报，包括以下步骤、:
步骤一：采集待预报地区的历史环境数据，并将历史环境数据依次输入两个或两个以上空气质量预报模型中，分别得到两个或两个以上历史空气质量预测数据；
步骤二：将两个或两个以上历史空气质量预测数据以实际观测数据为参考数据进行修正，计算出各空气质量预报模型在最终集合计算结果中所占权重值；
步骤三：采集待预测地区的当前环境数据，并将当前环境数据依次输入两个或两个以上空气质量预报模型中，分别得到两个或两个以上初始预测数据；
步骤四：根据权重值对初始预测数据进行修正，得到最终预报数据。
本发明所涉及的空气预报方法简单实用，可针对较小范围的区域(如城市、区县等)进行预报，预报精度高，且计算过程简单，方便采样。
进一步地，历史环境数据包括第一历史排放数据、历史第一地形数据以及历史气象数据。
优选地,当前环境数据包括第一当前排放数据、当前第一地形数据以及当前气象数据。
进一步地,在步骤一中，增加采集历史第二地形数据，并将历史气象数据与历史第二地形数据分别与两个或两个以上空气质量预报模型得到的基于历史第一地形数据的历史空气质量预测数据进行计算，分别得到两个或两个以上基于历史第二地形数据的初始预测数据。
优选地,步骤三中，增加采集当前第二地形数据，并将当前气象数据与当前第二地形数据分别与两个或两个以上空气质量预报模型得到的基于当前第一地形数据的当前空气质量预测数据进行计算，分别得到两个或两个以上基于当前第二地形数据的初始预测数据。
进一步地，在步骤一中，将历史第二地形数据，历史气象数据以及第二历史排放数据一起依次输入到两个或两个以上空气质量预报模型中，以嵌套计算的方式进行运算，分别得到两个或两个以上历史空气质量预测数据。
优选地，在步骤一中，以历史第二地形数据、历史人口分布数据和/或历史工厂分布数据对历史空气质量预测数据进行重采样运算，得到基于历史第一地形数据的历史空气质量预测数据。
进一步地，在步骤三中，将当前第二地形数据、当前气象数据以及第二当前排放数据一起依次输入两个或两个以上空气质量预报模型中，以嵌套计算的方式进行运算，分别得到两个或两个以上基于当前第二地形数据的初始预测数据。
优选地，在步骤三中，以当前第二地形数据、人口分布数据和/或工厂分布数据对历史空气质量预测数据进行重采样运算，得到基于当前第二地形数据的初始预测数据。
进一步地，空气质量预报模型为三个。
优选地，空气质量预报方法主要对空气中PM2.5、PM10、二氧化硫、二氧化氮、一氧化碳、臭氧中一种或多种的含量进行预报。
如上，本发明所涉及的空气质量预报方法预报过程简单，采样方便，且预报精确度高，还可改善因排放源数据缺乏造成的城市尺度空气质量预报精度较差的问题。
附图说明
图1为本发明空气质量预报方法的流程图。
图2为空气预报分析区域示例图；
图3(a)为PM10含量预报结果示例图；
图3(b)为二氧化硫含量预报结果示例图；
图3(c)为二氧化氮含量预报结果示例图；
图3(d)为臭氧含量预报结果示例图；
图3(e)为一氧化碳含量预报结果示例图；
图3(f)为PM2.5含量预报结果示例图。
具体实施方式
具体来说，本发明所涉及的空气质量预报方法是一种基于数据同化方法和待预测地区环境数据的多模式空气质量预报方法。
以下由特定的具体实施例说明本发明的实施方式，本领域技术人员可由本说明书所揭示的内容轻易地了解本发明的其他优点及功效。虽然本发明的描述将结合较佳实施例一起介绍，但这并不代表此发明的特征仅限于该实施方式。恰恰相反，结合实施方式作发明介绍的目的是为了覆盖基于本发明的权利要求而有可能延伸出的其它选择或改造。为了提供对本发明的深度了解，以下描述中将包含许多具体的细节。本发明也可以不使用这些细节实施。此外，为了避免混乱或模糊本发明的重点，有些具体细节将在描述中被省略。
另外，在以下的说明中所使用的“上”、“下”、“左”、“右”、“顶”、“底”，不应理解为对本发明的限制。
如图1所示，本发明提供一种空气质量预报方法,对空气中一种或多种污染物含量进行预报，在本发明的空气质量预报方法主要对空气中PM2.5、PM10、二氧化硫、二氧化氮、一氧化碳、臭氧中一种或多种的含量进行预报，包括以下步骤:
步骤一：采集待预报地区的历史环境数据，并将历史环境数据依次输入两个或两个以上空气质量预报模型中，分别得到两个或两个以上历史空气质量预测数据；
步骤二：将两个或两个以上历史空气质量预测数据以实际观测数据为参考数据进行修正，计算出各空气质量预报模型在最终集合计算结果中所占权重值；
步骤三：采集待预测地区的当前环境数据，并将当前环境数据依次输入两个或两个以上空气质量预报模型中，分别得到两个或两个以上初始预测数据；
步骤四：根据权重值对初始预测数据进行修正，得到最终预报数据。
在本发明中,历史环境数据优选包括第一历史排放数据、历史第一地形数据以及历史气象数据；当前环境数据优选包括当前排放数据、当前第一地形数据以及当前气象数据。
在步骤一中，首先采集第一历史排放数据、历史第一地形数据以及历史气象数据，利用两个或两个以上空气质量预报模型对待预报地区进行历史时段的空气质量预报，得到历史空气质量预测数据。
进一步地，待预测地区优选为城市，为了使计算所得结果更为精确，在步骤一中，增加采集历史第二地形数据，并将历史气象数据与历史第二地形数据分别与两个或两个以上空气质量预报模型得到的基于历史第一地形数据的历史空气质量预测数据进行计算，分别得到两个或两个以上基于第二历史地形数据的历史空气质量预测数据，具体来说，即增加采集历史第二地形数据，并将历史环境数据与历史第二地形数据与基于第一地形数据的历史空气质量预测数据进行如下方式之一的操作：
方式一：采集历史第二地形数据以及待预报城市区域所对应的第二历史排放数据，并将历史第二地形数据、历史气象数据以及第二历史排放数据一起依次输入两个或两个以上空气质量预报模型中，以嵌套计算的方式分别得到两个或两个以上基于第二历史地形数据的空气质量预测数据；在该方式中，对第二历史排放数据的测量精度要求较高，以获得较为精确的计算结果。
方式二：由于在实际操作中，高精度的第二历史排放数据往往较难获得，因此，在本发明中，针对实际上缺乏高精度的第二历史排放源数据的大多数城市区域的情况，利用方式一中基于第二历史排放数据的多层嵌套计算可用基于人口分布数据和/或工厂分布数据的重采样计算替代，用以提高城市区域预报的精度。即利用历史第二地形数据，以及待预报地区的人口分布数据和/或工厂分布数据对基于历史第一地形数据计算得到的历史空气质量预测数据进行重采样运算，得到基于第二地形数据的空气质量预测数据。
在步骤二中，将步骤一中所得到的历史空气质量预测数据与已采集的实际观测数据进行比对，以实际观测数据为参考对历史空气质量预测数据进行数据同化，对历史空气质量预测数据进行修正，得到修正后的历史空气质量预测数据和误差值；再分别计算各个空气质量预报模型所得到的数据与实际观测数据的相关性，并根据相关系数计算出各空气质量预报模型在最终集合计算结果中所占权重值。
在步骤三中，采集待预测地区的当前环境数据，即第一当前排放数据、当前第一地形数据以及当前气象数据，利用多个空气质量预报模式进行未来时段的预报，并利用前述计算得到的权重进行加权集成计算得到最终的集合预报结果。
进一步地，在步骤三中，待预测地区优选为城市，为了使计算所得结果更为精确，增加采集当前第二地形数据，并将当前气象数据与当前第二地形数据分别与两个或两个以上空气质量预报模型得到的基于当前第一地形数据的当前空气质量预测数据进行计算，分别得到两个或两个以上基于当前第二地形数据的初始预测数据，具体来说，即增加采集当前第二地形数据，并将当前气象数据与当前第二地形数据与基于当前第一地形数据的当前空气质量预测数据进行如下方式之一的操作：
方式一：采集当前第二地形数据以及待预报城市区域所对应的第二当前排放数据，并将当前第二地形数据、当前气象数据以及第二当前排放数据一起依次输入两个或两个以上空气质量预报模型中，以嵌套计算的方式分别得到两个或两个以上基于第二地形数据的空气质量预测数据；在该方式中，对第二当前排放数据的测量精度要求较高，以获得较为精确的计算结果。
方式二：进一步地，由于在实际操作中，高精度的第二当前排放数据往往较难获得，因此，在本发明中，针对实际上缺乏高精度的第二当前排放源数据的大多数城市区域的情况，利用方式一中基于第二当前排放数据的多层嵌套计算可用基于人口分布数据和/或工厂分布数据的重采样计算替代，用以提高城市区域预报的精度。即利用当前第二地形数据，以及待预报地区的人口分布数据和/或工厂分布数据对基于历史第一地形数据计算得到的历史空气质量预测数据进行重采样运算，得到基于当前第二地形数据的空气质量预测数据。
其中，当前/历史第二地形数据和第二当前/历史排放数据的分辨率高于当前/历史第一地形数据和第一当前/历史排放数据的分辨率，也就是说，当前/历史第一地形数据和第一当前/历史排放数据的每一个值覆盖的空间范围较大，当前/历史第二地形数据和第二当前/历史排放数据的每一个值覆盖的空间范围较小。
进一步地，在本发明中，基于当前第二地形数据、人口分布数据和/或工厂分布数据的重采样计算方式如下：
计算待预测城市区域范围内的当前第二地形数据平均值A、人口分布数据平均值P、工厂分布数据平均值F，按下式计算城市范围内的当前第二地形数据、人口分布数据，工厂分布数据每个网格的权重w：
其中a,p,f为各类数据每个网格的实际值。
将基于历史第一地形数据的历史空气质量预测数据中的每个值按下式重采样得到基于历史第二地形数据的预测数据中的多个值：
其中，v1是基于历史第一地形数据的预测数据的值，v2是基于历史第二地形数据的预测数据的值，n为每个v1的网格对应的v2的网格数量，wi为n个网格中第i个的权重。
因此，本发明所涉及的空气质量预报方法能利用多种空气质量预报模型的侧重于不同的空气物理化学变化机理的优点，改善传统算术平均集成方法所忽略的各模型精度不同的问题，同时对于数据缺乏的地区能够较大程度的提高区域模拟的精度。更进一步地，在本发明中，空气质量预报模型优选为三个：气象研究预报-区域建模分析系统(WeatherResearch&Forecast-Community Modeling and Analysis System，下文简称WRF-CMAQ)、化学物质含量预报模型(Weather Research Forecast coupled with Chemistry，下文简称WRF-Chem)、综合空气质量扩展模型(Comprehensive Air quality Model with Extensions，下文简称CAMx)。在本发明中，可以根据实际所需预报的污染物种类，选择上述模型中的一个或多个进行计算。
进一步地，在本发明中，还可以采用其他各种模型替代上述模型，如WRF-Chem模型可用MM5/RAMS/ARPS等模型进行替代；WRF-Chem/CAMx模型可用AERMOD/ADMS/UAM等模型进行替代。其中，MM5(Mesoscale Model 5)为中尺度预报模型、RAMS(Regional Atmosphere Modeling System)为区域大气建模系统、ARPS(AdvancedRegional Prediction System)为区域预测系统、AERMOD(AMS/EPA Regulatory Model)为大气模拟预测管理模型、ADMS(Atmospheric Diffusion Model System)为大气扩散模型系统、UAM(Urban Atmosphere Model)为城市大气模型。
但本发明并不局限于此，还可以用其他各种模型对上述模型进行替代，以获得所需的空气污染物含量。
在本发明中，在获得将对两个或两个以上历史空气质量预测数据后，以实际观测数据为参考数据进行修正所采用的方法为数据同化的方法。数据同化是分析处理随空间和时间分布的观测资料为数值预报提供初始场的一个过程，具体包括两层基本涵义：合理利用各种不同精度的非常规观测资料，把它们与常规观测资料融合为有机的整体，为数值预报提供更好的初始场，并综合利用不同时间段所采集的观测资料，将这些资料中所包含的时间演变信息转化为要素场的空间分布状况。数据同化从最初的逐步订正法，最优插值法，到现在的变分同化和集合卡尔曼滤波已经有多种同化方法。具体来说：
初步预报值定义为背景场xb，最终修正后的目标值为分析场xa，观测值构成观测场yo1)三维变分同化中，
xa＝xb+(B-1＋HTR-1H)-1HTR-1(yo-H(x))
其中，x，y变量为矩阵，H(x)为前向算子，形如y＝H(x)，是分析变量到观测变量的一个计算函数，最简单的为插值函数，即：(d为观测值与背景值的空间距离)；
T为矩阵转置运算；R为观测误差协方差矩阵；B为背景误差协方差矩阵。
建立目标函数如下：
Jb和Jo表示对背景场和观测场的拟合程度。J(x)对x的梯度函数为：：
通过迭代法求最小的x则对应最终的xa。其中，是J(x)对x的梯度函数，也就是J(x)对x的偏导数。
2)集合最优插值同化中，
xa＝xb+W(yo-H(xb))
其中，
W＝BHT(HTBH＋R)-1
更进一步，在本发明中，在计算各空气质量预报模型所得出结果的权重时，相关或相关系数是用于反映变量之间相关关系密切程度的统计指标，相关及相关系数可通过如下方法确定：
其中，x，y分别为观测值集合，和与观测值位置时刻一致的同化后的分析值集合。
对应的基于相关系数的权重确定方案按下式确定，其中，c1、c2、c3分别代表三个模型所计算出的数值，w1、w2、w3分别代表三个数值的权重：
w1＝c1/(c1+c2+c3)
w2＝c2/(c1+c2+c3)
w3＝c3/(c1+c2+c3)
基于上述计算方法，本发明的空气质量预报方法可以细化为如下步骤：
1、利用历史时段的数据进行预报和同化，求取多模式集合所需权重：
a)将历史气象数据，历史排放源数据和历史第一地形数据分别输入所选三个空气质量预报模式WRF-CMAQ，WRF-Chem，WRF-CAMx中，计算得到初始大尺度低分辨率的历史空气质量预测数据。
b)利用较为精确的排放源数据，历史第二地形数据和历史空气质量预测数据进行嵌套计算，得到三个模式的初始城市尺度高分辨率空气质量预报结果。
c)如果缺乏较为精确的排放源数据，则采用历史第二地形数据，人口分布数据和/或工厂分布数据对初始低分辨率的历史空气质量预测数据进行重采样计算，首先计算权重，再按下式重采样得到城市尺度高分辨率空气质量预报结果。
d)采用三维变分同化方法，对城市尺度预报结果和历史空气质量预测数据进行同化，得到修正后的城市空气质量预报结果和误差项。
e)计算三个修正后城市空气质量预报结果和观测数据的相关性，并计算对应三个模式的权重。
2、利用当前时段的数据对未来时段的空气质量进行预报：
a)将当前气象分析数据，当前排放源数据和当前第一地形数据分别输入所选三个空气质量预报模式WRF-CMAQ，WRF-Chem，WRF-CAMx中，计算得到初始大尺度低分辨率初始预测数据。
b)利用较为精细的当前排放源数据，当前第二地形数据和当前空气质量预测数据进行嵌套计算，得到三个模式的初始城市尺度高分辨率空气质量预报结果。
c)如果缺乏较为精细的当前排放源数据，则可采用当前第二地形数据，当前人口分布数据和/或工厂分布数据对初始低分辨率预测数据进行重采样计算，在计算权重后，得到当前城市尺度高分辨率空气质量预报结果。
3、将得到的三个当前城市尺度高分辨率空气质量预测数据按照历史空气质量预测数据计算得到的集成权重值进行集合加权计算，得到最终的空气质量预报结果。
下面通过两具体实施例对本发明进行详细说明：
第一实施例：
数据源采用美国国家环境预报中心(下文简称NCEP)提供的全球地理数据库，全球基础排放源数据，历史和即时的全球预报系统(下文简称GFS)气象场数据。实验所面向的区域为华北某城市范围，在3km分辨率尺度下由50*50的格网覆盖，图2显示了在嵌套计算方式下各层所对应的范围。分析时间为2015年12月5日，预报时间为12月6日00时至9日00时，总计72小时。历史数据采用12月1日至12月5日共120小时的数据。
步骤1：输入地理数据，36km全球排放源数据和1-5日的GFS气象数据至WRF-CMAQ，WRF-Chem，WRF-CAMx三种模式中，计算初始无嵌套预报结果，得到36km格网分析结果。统计PM2.5，PM10，NO2，SO2，O3，CO六种主要污染物。
步骤2：以嵌套计算方式处理1-5日36km分析结果，即修改三个模型中的嵌套计算方式，并输入嵌套计算后两层所需9km和3km北京区域详细排放源数据，处理得到3km分析结果。
步骤3：将1-5日3km污染物预报结果和地面观测站数据输入三维变分数据同化计算模型，计算得到同化后的三种模式的3km污染物数据。
步骤4：计算同化后三种模式3km污染物数据与地面观测数据的相关系数，分别为：0.782,0.634,0.699。并以三个相关系数，确定三个模型对应的权重值分别为：0.370，0.300，0.330.
步骤5：输入地理数据，36km全球排放源数据和6-8日的GFS气象分析数据计算三个模式的未来时段预报结果，格网大小36km。
步骤6：以嵌套方式计算6-8日3km分析结果。
步骤7：以步骤4得到的权重为系数，计算集合预报的最终结果。
步骤8，输出空气中各污染物/颗粒物浓度结果(如图3(a)至图3(f)所示)。
第二实施例：
步骤1：输入地理数据，36km全球排放源数据和1-5日的GFS气象数据至WRF-CMAQ，WRF-Chem，WRF-CAMx三种模式中，计算初始无嵌套预报结果，得到36km格网分析结果。统计PM2.5，PM10，NO2，SO2，O3，CO六种主要污染物。
步骤2：以重采样方式处理1-5日36km分析结果。
步骤2-1：以3km分辨率的北京城区地形高度数据，人口分布数据，工厂分布数据计算重采样权重值。
步骤2-2：以3km重采样权重值为依据对36km分析结果对应区域进行重采样。
步骤3：将1-5日3km污染物预报结果和地面观测站数据输入三维变分数据同化计算模型，计算得到同化后的三种模式的3km污染物数据。
步骤4：计算同化后三种模式3km污染物数据与地面观测数据的相关系数，分别为：0.754,0.579,0.649。并以三个相关系数，确定三个模型对应的权重值：0.382，0.291，0.327
步骤5：输入地理数据，3km全球排放源数据和6-8日的GFS气象数据计算三个模式的预报结果。
步骤6：以重采样方式计算6-8日3km分析结果。
步骤7：以步骤4得到的权重为系数，计算集合预报的最终结果。
步骤8，输出颗粒物浓度结果。
综上，本发明所涉及的空气质量预报方法通过数据同化后的精度和相关性分析来评估多个空气质量预报模式在最终集合预报结果中的贡献度，并且利用城市地形数据，人口经济数据与雾霾产生的相关性来提高初始预报的分辨率，并从而解决精细排放源数据难以获取的问题。此外，本发明经过数据同化辅助提供多模式集合预报的结果精度得到提高，还可在在显著降低数据获取成本的同时提高城市级别空气预报的精度，值得进行广泛推广应用。
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。
Claims (11)
1.一种空气质量预报方法,对空气中一种或多种气体或颗粒物的含量进行预报，其特征在于,包括以下步骤:
步骤一：采集待预报地区的历史环境数据，并将所述历史环境数据依次输入两个或两个以上空气质量预报模型中，分别得到两个或两个以上历史空气质量预测数据；
步骤二：将所述两个或两个以上历史空气质量预测数据以实际观测数据为参考数据进行修正，计算出各所述空气质量预报模型在最终集合计算结果中所占权重值；
步骤三：采集所述待预测地区的当前环境数据，并将所述当前环境数据依次输入两个或两个以上空气质量预报模型中，分别得到两个或两个以上初始预测数据；
步骤四：根据所述权重值对所述初始预测数据进行修正，得到最终预报数据。
2.如权利要求1所述的空气质量预报方法，其特征在于，所述历史环境数据包括第一历史排放数据、历史第一地形数据以及历史气象数据。
3.如权利要求1所述的空气质量预报方法，其特征在于，所述当前环境数据包括第一当前排放数据、当前第一地形数据以及当前气象数据。
4.如权利要求2所述的空气质量预报方法，其特征在于，在所述步骤一中，增加采集历史第二地形数据，并将所述历史气象数据与所述历史第二地形数据分别与两个或两个以上空气质量预报模型得到的基于所述历史第一地形数据的所述历史空气质量预测数据进行计算，分别得到两个或两个以上基于所述第二历史地形数据的历史空气质量预测数据。
5.如权利要求3所述的空气质量预报方法，其特征在于，所述步骤三中，增加采集当前第二地形数据，并将所述当前气象数据与所述当前第二地形数据分别与两个或两个以上空气质量预报模型得到的基于所述当前第一地形数据的所述当前空气质量预测数据进行计算，分别得到两个或两个以上基于所述当前第二地形数据的初始预测数据。
6.如权利要求4所述的空气质量预报方法，其特征在于，在步骤一中，将所述历史第二地形数据，所述历史气象数据以及第二历史排放数据一起依次输入到所述两个或两个以上空气质量预报模型中，以嵌套计算的方式进行运算，分别得到两个或两个以上历史空气质量预测数据。
7.如权利要求4所述的空气质量预报方法，其特征在于，在步骤一中，以所述历史第二地形数据、历史人口分布数据和/或历史工厂分布数据对所述历史空气质量预测数据进行重采样运算，得到所述基于所述历史第一地形数据的历史空气质量预测数据。
8.如权利要求5所述的空气质量预报方法，其特征在于，在步骤三中，将所述当前第二地形数据、所述当前气象数据以及第二当前排放数据输入到所述两个或两个以上空气质量预报模型中，以嵌套计算的方式进行运算，得到基于所述当前第二地形数据的初始预测数据。
9.如权利要求5所述的空气质量预报方法，其特征在于，在步骤三中，以所述当前第二地形数据、人口分布数据和/或工厂分布数据对所述历史空气质量预测数据进行重采样运算，得到基于所述当前第二地形数据的初始预测数据。
10.如权利要求1所述的空气质量预报方法，其特征在于，空气质量预报模型为三个。
11.如权利要求1所述的空气质量预报方法，其特征在于，所述空气质量预报方法主要对空气中PM2.5、PM10、二氧化硫、二氧化氮、一氧化碳、臭氧中一种或多种的含量进行预报。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610218043.6A CN107273995A (zh) | 2016-04-08 | 2016-04-08 | 空气质量预报方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610218043.6A CN107273995A (zh) | 2016-04-08 | 2016-04-08 | 空气质量预报方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107273995A true CN107273995A (zh) | 2017-10-20 |
Family
ID=60052519
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610218043.6A Pending CN107273995A (zh) | 2016-04-08 | 2016-04-08 | 空气质量预报方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107273995A (zh) |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109492907A (zh) * | 2018-11-08 | 2019-03-19 | 成都市环境保护科学研究院 | 基于cmaq模型的空气质量措施评估方法、系统、存储介质和终端 |
| CN109543906A (zh) * | 2018-11-23 | 2019-03-29 | 长三角环境气象预报预警中心（上海市环境气象中心） | 一种大气能见度预测的方法及设备 |
| CN109709577A (zh) * | 2018-12-28 | 2019-05-03 | 南京雨后地软环境技术有限公司 | 一种基于WRF-Chem模式的气溶胶激光雷达反演PM2.5的三维变分同化方法 |
| CN109858686A (zh) * | 2019-01-14 | 2019-06-07 | 南京大学 | 一种基于EnKF的地面排放清单反演优化方法 |
| CN109948840A (zh) * | 2019-03-08 | 2019-06-28 | 宁波市气象台 | 一种空气质量预报方法 |
| CN109991379A (zh) * | 2019-05-10 | 2019-07-09 | 苏州华能检测技术有限公司 | 一种空气中硫化氢含量检测方法 |
| CN110009245A (zh) * | 2019-04-12 | 2019-07-12 | 阳江职业技术学院 | 基于神经网络的室内空气质量预测方法、装置及电子设备 |
| CN110045440A (zh) * | 2019-04-26 | 2019-07-23 | 国网湖南省电力有限公司 | 一种大气雾霾数值预报的集合预报方法及系统 |
| CN110261547A (zh) * | 2019-07-04 | 2019-09-20 | 北京思路创新科技有限公司 | 一种空气质量预报方法和设备 |
| CN110346518A (zh) * | 2019-07-25 | 2019-10-18 | 中南大学 | 一种交通排放污染可视化预警方法及其系统 |
| CN110346517A (zh) * | 2019-07-25 | 2019-10-18 | 中南大学 | 一种智慧城市工业大气污染可视化预警方法及其系统 |
| CN110968926A (zh) * | 2018-09-29 | 2020-04-07 | 中国科学院国家空间科学中心 | 一种基于改进背景误差协方差矩阵预测大气参数的方法 |
| CN110989043A (zh) * | 2019-12-25 | 2020-04-10 | 中科三清科技有限公司 | 空气质量指数级别概率预报方法、装置及存储介质 |
| CN111487216A (zh) * | 2020-05-06 | 2020-08-04 | 北京中科锐景科技有限公司 | 一种二氧化碳通量反演方法、系统 |
| CN111860996A (zh) * | 2020-07-15 | 2020-10-30 | 中科三清科技有限公司 | 一种空气质量模型的排放源数据生成方法及装置、介质 |
| CN113642223A (zh) * | 2021-10-13 | 2021-11-12 | 中科三清科技有限公司 | 空气质量预报方法、装置、设备和计算机可读存储介质 |
| CN114565136A (zh) * | 2022-01-27 | 2022-05-31 | 内蒙古工业大学 | 一种基于生成对抗网络的空气质量预测优化方法 |
| CN116663786A (zh) * | 2023-08-02 | 2023-08-29 | 北京东润环能科技股份有限公司 | 一种风光资源评估方法、装置、电子设备和存储介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102799772A (zh) * | 2012-07-03 | 2012-11-28 | 中山大学 | 面向空气质量预报的样本优化方法 |
| CN103294928A (zh) * | 2013-06-28 | 2013-09-11 | 华北电力大学（保定） | 一种碳排放组合预测方法 |
| WO2014020778A1 (ja) * | 2012-07-31 | 2014-02-06 | 日揮株式会社 | 気象再現方法、気象再現装置、及び空気利用装置 |
| CN104881582A (zh) * | 2015-05-28 | 2015-09-02 | 象辑知源（武汉）科技有限公司 | 空气质量的预测方法和装置 |
| CN105069537A (zh) * | 2015-08-25 | 2015-11-18 | 中山大学 | 一种组合式空气质量预报模型的构建方法 |
-
2016
- 2016-04-08 CN CN201610218043.6A patent/CN107273995A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102799772A (zh) * | 2012-07-03 | 2012-11-28 | 中山大学 | 面向空气质量预报的样本优化方法 |
| WO2014020778A1 (ja) * | 2012-07-31 | 2014-02-06 | 日揮株式会社 | 気象再現方法、気象再現装置、及び空気利用装置 |
| CN103294928A (zh) * | 2013-06-28 | 2013-09-11 | 华北电力大学（保定） | 一种碳排放组合预测方法 |
| CN104881582A (zh) * | 2015-05-28 | 2015-09-02 | 象辑知源（武汉）科技有限公司 | 空气质量的预测方法和装置 |
| CN105069537A (zh) * | 2015-08-25 | 2015-11-18 | 中山大学 | 一种组合式空气质量预报模型的构建方法 |
Non-Patent Citations (3)
| Title |
|---|
| 曹飞等: ""中国碳排放的组合预测与减排缺口分析"", 《工业技术经济》 * |
| 王自发等: ""北京空气质量多模式集成预报系统的建立及初步应用"", 《南京信息工程大学学报(自然科学版)》 * |
| 黄思等: ""利用多模式集合和多元线性回归改进北京PM10预报"", 《环境科学学报》 * |
Cited By (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110968926A (zh) * | 2018-09-29 | 2020-04-07 | 中国科学院国家空间科学中心 | 一种基于改进背景误差协方差矩阵预测大气参数的方法 |
| CN109492907A (zh) * | 2018-11-08 | 2019-03-19 | 成都市环境保护科学研究院 | 基于cmaq模型的空气质量措施评估方法、系统、存储介质和终端 |
| CN109543906A (zh) * | 2018-11-23 | 2019-03-29 | 长三角环境气象预报预警中心（上海市环境气象中心） | 一种大气能见度预测的方法及设备 |
| CN109543906B (zh) * | 2018-11-23 | 2024-04-16 | 长三角环境气象预报预警中心（上海市环境气象中心） | 一种大气能见度预测的方法及设备 |
| CN109709577A (zh) * | 2018-12-28 | 2019-05-03 | 南京雨后地软环境技术有限公司 | 一种基于WRF-Chem模式的气溶胶激光雷达反演PM2.5的三维变分同化方法 |
| CN109858686A (zh) * | 2019-01-14 | 2019-06-07 | 南京大学 | 一种基于EnKF的地面排放清单反演优化方法 |
| CN109948840B (zh) * | 2019-03-08 | 2020-02-21 | 宁波市气象台 | 一种空气质量预报方法 |
| CN109948840A (zh) * | 2019-03-08 | 2019-06-28 | 宁波市气象台 | 一种空气质量预报方法 |
| CN110009245A (zh) * | 2019-04-12 | 2019-07-12 | 阳江职业技术学院 | 基于神经网络的室内空气质量预测方法、装置及电子设备 |
| CN110045440A (zh) * | 2019-04-26 | 2019-07-23 | 国网湖南省电力有限公司 | 一种大气雾霾数值预报的集合预报方法及系统 |
| CN110045440B (zh) * | 2019-04-26 | 2021-03-09 | 国网湖南省电力有限公司 | 一种大气雾霾数值预报的集合预报方法及系统 |
| CN109991379A (zh) * | 2019-05-10 | 2019-07-09 | 苏州华能检测技术有限公司 | 一种空气中硫化氢含量检测方法 |
| CN110261547A (zh) * | 2019-07-04 | 2019-09-20 | 北京思路创新科技有限公司 | 一种空气质量预报方法和设备 |
| CN110346517A (zh) * | 2019-07-25 | 2019-10-18 | 中南大学 | 一种智慧城市工业大气污染可视化预警方法及其系统 |
| CN110346518A (zh) * | 2019-07-25 | 2019-10-18 | 中南大学 | 一种交通排放污染可视化预警方法及其系统 |
| CN110989043A (zh) * | 2019-12-25 | 2020-04-10 | 中科三清科技有限公司 | 空气质量指数级别概率预报方法、装置及存储介质 |
| CN111487216B (zh) * | 2020-05-06 | 2023-01-24 | 北京中科锐景科技有限公司 | 一种二氧化碳通量反演方法、系统 |
| CN111487216A (zh) * | 2020-05-06 | 2020-08-04 | 北京中科锐景科技有限公司 | 一种二氧化碳通量反演方法、系统 |
| CN111860996A (zh) * | 2020-07-15 | 2020-10-30 | 中科三清科技有限公司 | 一种空气质量模型的排放源数据生成方法及装置、介质 |
| CN113642223A (zh) * | 2021-10-13 | 2021-11-12 | 中科三清科技有限公司 | 空气质量预报方法、装置、设备和计算机可读存储介质 |
| CN113642223B (zh) * | 2021-10-13 | 2022-05-06 | 中科三清科技有限公司 | 空气质量预报方法、装置、设备和计算机可读存储介质 |
| CN114565136A (zh) * | 2022-01-27 | 2022-05-31 | 内蒙古工业大学 | 一种基于生成对抗网络的空气质量预测优化方法 |
| CN114565136B (zh) * | 2022-01-27 | 2024-01-26 | 内蒙古工业大学 | 一种基于生成对抗网络的空气质量预测优化方法 |
| CN116663786A (zh) * | 2023-08-02 | 2023-08-29 | 北京东润环能科技股份有限公司 | 一种风光资源评估方法、装置、电子设备和存储介质 |
| CN116663786B (zh) * | 2023-08-02 | 2023-11-21 | 北京东润环能科技股份有限公司 | 一种风光资源评估方法、装置、电子设备和存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107273995A (zh) | 空气质量预报方法 | |
| CN112905560B (zh) | 一种多源时空大数据深度融合的空气污染预测方法 | |
| Vignati et al. | Global scale emission and distribution of sea-spray aerosol: Sea-salt and organic enrichment | |
| CN111508569A (zh) | 一种基于土壤传递函数的目标土壤性质含量预测方法 | |
| CN111368401A (zh) | 污染源的溯源方法、装置和存储介质 | |
| CN111896680B (zh) | 基于卫星遥感数据的温室气体排放分析方法及系统 | |
| CN109146161A (zh) | 融合栈式自编码和支持向量回归的pm2.5浓度预测方法 | |
| CN113348471A (zh) | 一种大气污染预测中对区域边界优化的方法 | |
| CN111651873B (zh) | 城市大气VOCs排放清单反演方法、系统、终端以及储存介质 | |
| CN111353634B (zh) | 一种基于模式初始场同化的二氧化硫排放源优化方法 | |
| Deng et al. | The MR-CA models for analysis of pollution sources and prediction of PM 2.5 | |
| CN110766222A (zh) | 基于粒子群参数优化和随机森林的pm2.5浓度预测方法 | |
| Metia et al. | Urban air pollution estimation using unscented Kalman filtered inverse modeling with scaled monitoring data | |
| CN108764527B (zh) | 一种土壤有机碳库时空动态预测最优环境变量筛选方法 | |
| Gao et al. | Sensitivity analysis of surface ozone to emission controls in Beijing and its neighboring area during the 2008 Olympic Games | |
| CN113011455B (zh) | 一种空气质量预测svm模型构建方法 | |
| CN113095499A (zh) | 一种绝缘子等值附盐密度的预测方法 | |
| CN117669270B (zh) | 一种对微站组网数据进行时空一致性订正的方法 | |
| Ren et al. | Analysis of the spatial characteristics of inhalable particulate matter concentrations under the influence of a three-dimensional landscape pattern in Xi'an, China | |
| CN116611547A (zh) | 大气污染物来源预报方法、装置、设备及存储介质 | |
| CN116223395A (zh) | 近地表痕量气体浓度反演模型及反演方法 | |
| Zhu et al. | The improvement and comparison of diffuse radiation models in different climatic zones of China | |
| CN115881239A (zh) | 基于扩展wrf和cmaq模型的大气污染弱扩散区域的划分方法 | |
| CN113834902B (zh) | 一种基于四维变分同化的二氧化硫排放源反演方法 | |
| CN114529035A (zh) | 一种基于cart的多模式集成模型的风速预报方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171020 |