CN108304610B - 一种大气高污染过程动态追踪方法 - Google Patents
一种大气高污染过程动态追踪方法 Download PDFInfo
- Publication number
- CN108304610B CN108304610B CN201711404530.2A CN201711404530A CN108304610B CN 108304610 B CN108304610 B CN 108304610B CN 201711404530 A CN201711404530 A CN 201711404530A CN 108304610 B CN108304610 B CN 108304610B
- Authority
- CN
- China
- Prior art keywords
- pollution
- state
- atmospheric pollution
- atmospheric
- days
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种大气高污染过程动态追踪方法。该方法首先利用连续多日大气污染状态估计,结合高污染发生可能性诊断,实现当前日期下对未来连续多日的高污染发生的追踪诊断;并跟随当前日期的更新,实现追踪时期的动态更新,动态循环未来连续多日高污染发生追踪诊断,最终实现对大气高污染过程的动态追踪。
Description
技术领域
本发明涉及环境质量预报预警技术领域,具体涉及一种大气高污染过程动态追踪方法。
背景技术
目前在中国,冬季大气高污染带来的灰霾天气频发,高温低湿常天气下也常现光化学烟雾高污染事件,严重危害居民健康和生活,对其进行有效预防与控制至关重要。目前大气污染预测技术往往对高污染过程预测不精确,急需面向高污染过程的污染变化预测新方法,提高大气污染物高污染的预测能力。
目前常用的大气污染或者空气质量预测方法主要有基于污染扩散过程模拟的数值方法和基于数据驱动的统计预测方法两大类。数值预报由于计算十分复杂,所需输入的排放清单和边界气象场不确定性较大,对高污染事件的预测精度不高。随着大气污染监测网络发展带来更多的研究数据,数据驱动方法的力量逐渐超于数值方法,主要有回归拟合、神经网络、小波分析、支持向量机、贝叶斯网络等。此类方法,虽使用方便,但是由于其对大气污染与气象关联关系表征不清晰,往往预测结果难以跟踪到高污染事件的增长趋势。因此,本研究面向大气高污染预测需求,建立了一种大气高污染过程追踪方法,用于实现大气高污染过程的饱和预测。
发明内容
本发明的目的是克服上述现有技术的缺点,提供一种大气高污染过程动态追踪方法,该方法基于大气污染与气象条件的影响关系,利用连续多日大气污染状态估计,结合高污染发生可能性诊断,实现了对大气高污染过程的动态追踪。
本发明是通过以下技术方案来实现的:
一种大气高污染过程动态追踪方法,具体步骤如下:
1S.确定追踪开始日期的大气污染状态以及未来连续多日的气象条件。
2S.基于气象对大气污染状态的影响关系,建立连续多日大气污染状态估计方法,滚动估计未来连续多日的大气污染状态。
3S.通过估计结果进行未来高污染发生可能性的诊断。
4S.当日期有更新时,更新开始日期,重复1S-3S。
作为上述方案的改进,所述的1S的具体步骤如下:
1.1S.根据国家大气污染的相关标准,确定大气污染状态的描述标准;根据研究区域气象条件对大气污染影响的分析,确定与大气污染影响相关的气象因子;大气污染通常指PM2.5等颗粒物污染、臭氧污染、氮氧化物污染、一氧化碳污染等;气象因子包括风向、风速、气温、相对湿度、降水量、大气压、天气形势。
2.2S.确定要追踪的时期范围,将追踪时期的前一日确定为追踪开始日期;获取开始日期的大气污染监测结果,确定其大气污染状态;获取未来连续多日的气象因子的预报信息,确定为未来连续多日的的气象条件;未来连续多日的长度由气象预报时长以及目前预报水平决定,通常为3~7日。
作为上述方案的改进,所述的2S中基于气象对大气污染状态的影响关系,建立连续多日大气污染状态估计方法如下:污染状态时间演变可看作是一个马尔可夫过程,污染状态演变过程模拟成一个n步状态转变链模型,以时间为节点,公式如下:
{XT,XT+1,…,XT+n-1,XT+n}={f(PT,XT-1),f(PT+1,XT),…,f(PT+n,XT+n-1)},n∈Z+ (1)
其中,函数f的含义:f(PT+n,XT+n-1)=PT+n·XT+n-1;
XT+n的形式:X=[x1,x2,…,xi,…,xs]’;
XT+n指的T+n时刻的污染状态向量;XT+n-1是T+n-1时刻的污染状态向量;
XT+n依赖于XT+n-1,PT+n是污染状态从T+n-1时刻到T+n时刻的转变关系矩阵;其中,X是大气污染状态向量形式,x1,x2,…,xi,…,xs可取值0~1;当大气污染状态为i时,xi为1,其余元素为0;s为大气污染状态的总数目;pij指状态xi转变为xj的发生概率。
因此,假设开始日期的大气污染状态向量为X0,则未来连续多日大气污染状态向量的估计方法公式如下:
X1=P1·X0
X2=P2·X1
X3=P3·X2
......
Xn=Pn·Xn-1 (2)
与一般马尔可夫过程中的状态转移矩阵不同,在此,转变关系矩阵Pn是一个动态变化量,受前一时刻污染状态以及当前气象条件影响,计算方法如下公式所示:
Pn=P(Xn=xi|Xn-1=xj,Yn=mk1,mf2,…mgk) (3)
在这里,Yn指的是n时刻气象条件的状态,xi,xj指大气污染的某个状态,mk1,mf2,…mgk分别指k个气象因子的某个状态。每一步Pn的计算基于历史污染事件计数实现。
基于连续多日大气污染状态向量估计结果,滚动估计未来连续多日的最大可能性发生的大气污染状态,公式如下:
X1最大可能状态→MAX(X1)=max(P1·X0)
X2最大可能状态→MAX(X2)=max(P2·X1)
X3最大可能状态→MAX(X3)=max(P3·X2)
......
Xn最大可能状态→MAX(Xn)=max(Pn·Xn-1) (4)
基于未来连续多日的大气污染状态估计结果,逐日判断大气污染状态是否达到或超过国家空气质量轻度污染水平:若存在某一日为是,则判定为追踪日期内有高污染过程发生,提取高污染过程的发生日期以及可能性概率;若所有日期均为否,则判定为追踪日期内无高污染过程发生。
本发明具有以下有益效果:有效利用连续多日大气污染状态估计,对大气高污染的过程实现准确性高、快速高效的动态追踪。
附图说明
图1是本发明的技术方案流程图。
图2是2017年1月13日至20日北京PM2.5高污染过程追踪结果图。
具体实施方式
实施例1
一种大气高污染过程动态追踪方法,具体步骤如下:
1S.确定追踪开始日期的大气污染状态以及未来连续多日的气象条件。
2S.基于气象对大气污染状态的影响关系,建立连续多日大气污染状态估计方法,滚动估计未来连续多日的大气污染状态。
3S.通过估计结果进行未来高污染发生可能性的诊断。
4S.当日期有更新时,更新开始日期,重复1S-3S。
作为上述方案的改进,所述的1S的具体步骤如下:
1.1S.根据国家大气污染的相关标准,确定大气污染状态的描述标准;根据研究区域气象条件对大气污染影响的分析,确定与大气污染影响相关的气象因子;大气污染通常指PM2.5等颗粒物污染、臭氧污染、氮氧化物污染、一氧化碳污染等;气象因子包括风向、风速、气温、相对湿度、降水量、大气压、天气形势。
2.2S.确定要追踪的时期范围,将追踪时期的前一日确定为追踪开始日期;获取开始日期的大气污染监测结果,确定其大气污染状态;获取未来连续多日的气象因子的预报信息,确定为未来连续多日的的气象条件;未来连续多日的长度由气象预报时长以及目前预报水平决定,通常为3~7日。
作为上述方案的改进,所述的2S中基于气象对大气污染状态的影响关系,建立连续多日大气污染状态估计方法如下:污染状态时间演变可看作是一个马尔可夫过程,污染状态演变过程模拟成一个n步状态转变链模型,以时间为节点,公式如下:
{XT,XT+1,…,XT+n-1,XT+n}={f(PT,XT-1),f(PT+1,XT),…,f(PT+n,XT+n-1)},n∈Z+ (1)
其中,函数f的含义:f(PT+n,XT+n-1)=PT+n·XT+n-1;
XT+n的形式:X=[x1,x2,…,xi,…,xs]’;
XT+n指的T+n时刻的污染状态向量;XT+n-1是T+n-1时刻的污染状态向量;
XT+n依赖于XT+n-1,PT+n是污染状态从T+n-1时刻到T+n时刻的转变关系矩阵;其中,X是大气污染状态向量形式,x1,x2,…,xi,…,xs可取值0~1;当大气污染状态为i时,xi为1,其余元素为0;s为大气污染状态的总数目;pij指状态xi转变为xj的发生概率。
因此,假设开始日期的大气污染状态向量为X0,则未来连续多日大气污染状态向量的估计方法公式如下:
X1=P1·X0
X2=P2·X1
X3=P3·X2
......
Xn=Pn·Xn-1 (2)
与一般马尔可夫过程中的状态转移矩阵不同,在此,转变关系矩阵Pn是一个动态变化量,受前一时刻污染状态以及当前气象条件影响,计算方法如下公式所示:
Pn=P(Xn=xi|Xn-1=xj,Yn=mk1,mf2,…mgk) (3)
在这里,Yn指的是n时刻气象条件的状态,xi,xj指大气污染的某个状态,mk1,mf2,…mgk分别指k个气象因子的某个状态。每一步Pn的计算基于历史污染事件计数实现。
基于连续多日大气污染状态向量估计结果,滚动估计未来连续多日的最大可能性发生的大气污染状态,公式如下:
X1最大可能状态→MAX(X1)=max(P1·X0)
X2最大可能状态→MAX(X2)=max(P2·X1)
X3最大可能状态→MAX(X3)=max(P3·X2)
......
Xn最大可能状态→MAX(Xn)=max(Pn·Xn-1) (4)
实施例2
2017年1月13日至1月20日北京PM2.5高污染过程追踪
步骤1、确定大气污染状态的描述标准,以及表征气象条件的因素。
根据国家空气质量标准,确定PM2.5污染状态描述如表1所示,PM2.5浓度监测日均数据根据这一对应规则均可转化为六大等级水平。
表1 PM2.5污染的污染状态描述
参考北京PM2.5受气象条件影响的大量研究结果,影响北京PM2.5的气象条件主要有地面气象条件和宏观天气形势两大类,对以上两大类别中的气象因子进行与北京PM2.5污染的相关性分析,分析结果见表2,得到风速、相对湿度、环流形势与PM2.5污染的相关性最为显著。另外,虽然降水量的相关关系并不显著,但是众所周知,降水水平较大时对PM2.5污染冲刷作用显著,其作用并未与风速、相对湿度重叠。因此,最终确定影响北京PM2.5污染的主要气象因子为降水、风速、相对湿度、环流形势四大因子。
表2七大基础气象因子与PM2.5污染的相关关系
获取开始日期2017年1月12日的大气污染监测结果为61μg/m3,确定其大气污染状态为2级;从气象数据发布网站获取2017年1月13日至1月20日期间气象因子的预报信息。
步骤2、北京PM2.5高污染的高发时期为秋冬季节,因此收集2013至2016年秋冬季节的北京PM2.5污染数据,以及相匹配日期的气象因子数据,作为每一步Pn计算所需的历史污染事件集。
未来连续五日的大气污染状态向量估计结果为:
日期 | 估计的大气污染状态向量 |
2017年1月13日 | [83% 0% 17% 0% 0% 0%]' |
2017年1月14日 | [25% 50% 25% 0% 0% 0%]' |
2017年1月15日 | [0% 67% 0% 33% % 0%]' |
2017年1月16日 | [0% 25% 25% 50% 0% 0%]' |
2017年1月17日 | [0% 0% 0% 0% 10% 90%]' |
基于连续五日大气污染状态向量估计结果,滚动估计未来连续五日的最大可能性发生的大气污染状态,结果为:
步骤3、基于连续五日的大气污染状态估计结果,逐日判断大气污染状态是否达到或超过国家空气质量轻度污染水平(4级):
判别结果为:1月13日至17日将有高污染过程发生,高污染发生日期为1月16日至17日,发生可能性为50%~90%。
步骤4、时间更新至1月2日时,则更新追踪的开始日期为1月2日,重复步骤2-3,滚动估计14日至18日的PM2.5污染状态。
时间更新至1月3日时,则更新追踪的开始日期为1月3日,重复步骤2-3,滚动估计15日至19日的PM2.5污染状态。
时间更新至1月4日时,则更新追踪的开始日期为1月4日,重复步骤2-3,滚动估计16日至20日的PM2.5污染状态。
循环至开始日期更新为1月4日结束,完成了1月13日至1月20日的PM2.5高污染过程动态追踪。最终得到的追踪结果如图2所示。从图中可看出,经过追踪可预测到:2017年1月13日至20日期间将可能有PM2.5高污染发生,发生日期可能为1月16日至1月17日,发生可能性概率为50%~90%。追踪估计结果与实际监测结果基本一致,分析得出两个结果的相关系数达到0.97。
上列详细说明是针对本发明可行实施例的具体说明,该实施例并非用以限制本发明的专利范围,凡未脱离本发明所为的等效实施或变更,均应包含于本案的专利范围中。
Claims (2)
1.一种大气高污染过程动态追踪方法,其特征在于,具体步骤如下:
1S.确定追踪开始日期的大气污染状态以及未来连续多日的气象条件;
2S.基于气象对大气污染状态的影响关系,建立连续多日大气污染状态估计方法,滚动估计未来连续多日的大气污染状态;
3S.通过估计结果进行未来高污染发生可能性的诊断;
4S.当日期有更新时,更新开始日期,重复1S-3S;
所述的1S的具体步骤如下:
1.1S.根据国家大气污染的相关标准,确定大气污染状态的描述标准;根据研究区域气象条件对大气污染影响的分析,确定与大气污染影响相关的气象因子;大气污染通常指PM2.5颗粒物污染、臭氧污染、氮氧化物污染、一氧化碳污染;气象因子包括风向、风速、气温、相对湿度、降水量、大气压、天气形势;
2.2S.确定要追踪的时期范围,将追踪时期的前一日确定为追踪开始日期;获取开始日期的大气污染监测结果,确定其大气污染状态;获取未来连续多日的气象因子的预报信息,确定为未来连续多日的的气象条件;未来连续多日的长度由气象预报时长以及目前预报水平决定,通常为3~7日,
所述的2S中基于气象对大气污染状态的影响关系,建立连续多日大气污染状态估计方法如下:污染状态时间演变可看作是一个马尔可夫过程,污染状态演变过程模拟成一个n步状态转变链模型,以时间为节点,公式如下:
{XT,XT+1,…,XT+n-1,XT+n}={f(PT,XT-1),f(PT+1,XT),…,f(PT+n,XT+n-1)},n∈Z+ (1)
其中,函数f的含义:f(PT+n,XT+n-1)=PT+n XT+n-1;
XT+n的形式:X=[x1,x2,…,xi,…,xs]’;
XT+n指的T+n时刻的污染状态向量;XT+n-1是T+n-1时刻的污染状态向量;XT+n依赖于XT+n-1,PT+n是污染状态从T+n-1时刻到T+n时刻的转变关系矩阵;其中,X是大气污染状态向量形式,x1,x2,…,xi,…,xs可取值0~1;当大气污染状态为i时,xi为1,其余元素为0;s为大气污染状态的总数目;pij指状态xi转变为xj的发生概率;
因此,假设开始日期的大气污染状态向量为X0,则未来连续多日大气污染状态向量的估计方法公式如下:
与一般马尔可夫过程中的状态转移矩阵不同,在此,转变关系矩阵Pn是一个动态变化量,受前一时刻污染状态以及当前气象条件影响,计算方法如下公式所示:
Pn=P(Xn=xi|Xn-1=xj,Yn=mk1,mf2,…mgk) (3)
在这里,Yn指的是n时刻气象条件的状态,xi,xj指大气污染的某个状态,mk1,mf2,…mgk分别指k个气象因子的某个状态,每一步Pn的计算基于历史污染事件计数实现。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711404530.2A CN108304610B (zh) | 2017-12-22 | 2017-12-22 | 一种大气高污染过程动态追踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711404530.2A CN108304610B (zh) | 2017-12-22 | 2017-12-22 | 一种大气高污染过程动态追踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108304610A CN108304610A (zh) | 2018-07-20 |
CN108304610B true CN108304610B (zh) | 2021-05-18 |
Family
ID=62870333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711404530.2A Active CN108304610B (zh) | 2017-12-22 | 2017-12-22 | 一种大气高污染过程动态追踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108304610B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110472354B (zh) * | 2019-08-20 | 2023-07-18 | 东南大学 | 一种新能源汽车渗透的环境影响评估方法 |
CN116070923B (zh) * | 2023-02-15 | 2024-05-07 | 中科三清科技有限公司 | 大气污染情景模拟方法、装置和电子设备 |
CN117933466A (zh) * | 2024-01-23 | 2024-04-26 | 中科三清科技有限公司 | 污染预测方法、装置、存储介质及电子设备 |
CN118070110B (zh) * | 2024-04-22 | 2024-06-28 | 无锡九方科技有限公司 | 基于时域数据分析的大气污染物监测方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105493109A (zh) * | 2013-06-05 | 2016-04-13 | 微软技术许可有限责任公司 | 使用多个数据源的空气质量推断 |
CN105738974A (zh) * | 2016-02-05 | 2016-07-06 | 北京市环境保护监测中心 | 空气重污染天气的预报方法和*** |
WO2017051411A1 (en) * | 2015-09-24 | 2017-03-30 | Agt International Gmbh | Near real-time modeling of pollution dispersion |
CN106650017A (zh) * | 2016-11-23 | 2017-05-10 | 大连理工大学 | 一种利用固定位置污染物探测器辨识城市空间多污染物源的方法 |
CN106682381A (zh) * | 2015-11-10 | 2017-05-17 | 中国科学院沈阳计算技术研究所有限公司 | 一种面向环境空气质量的动态数据模拟预报方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106650825B (zh) * | 2016-12-31 | 2020-05-12 | 中国科学技术大学 | 一种机动车尾气排放数据融合*** |
-
2017
- 2017-12-22 CN CN201711404530.2A patent/CN108304610B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105493109A (zh) * | 2013-06-05 | 2016-04-13 | 微软技术许可有限责任公司 | 使用多个数据源的空气质量推断 |
WO2017051411A1 (en) * | 2015-09-24 | 2017-03-30 | Agt International Gmbh | Near real-time modeling of pollution dispersion |
CN106682381A (zh) * | 2015-11-10 | 2017-05-17 | 中国科学院沈阳计算技术研究所有限公司 | 一种面向环境空气质量的动态数据模拟预报方法 |
CN105738974A (zh) * | 2016-02-05 | 2016-07-06 | 北京市环境保护监测中心 | 空气重污染天气的预报方法和*** |
CN106650017A (zh) * | 2016-11-23 | 2017-05-10 | 大连理工大学 | 一种利用固定位置污染物探测器辨识城市空间多污染物源的方法 |
Non-Patent Citations (1)
Title |
---|
"城市空气污染指数预报方法研究";张雅斌 等;《陕西气象》;20030322;第24-27页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108304610A (zh) | 2018-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108304610B (zh) | 一种大气高污染过程动态追踪方法 | |
JP6129028B2 (ja) | 建物動力設備のエネルギー消費量予測方法 | |
CN106779129A (zh) | 一种考虑气象因素的短期电力负荷预测方法 | |
CN113919231A (zh) | 基于时空图神经网络的pm2.5浓度时空变化预测方法及*** | |
CN106600037B (zh) | 一种基于主成分分析的多参量辅助负荷预测方法 | |
CN111191193A (zh) | 一种基于自回归滑动平均模型的长期土壤温湿度高精度预测方法 | |
CN112116149B (zh) | 一种考虑预报不确定性关联演化特征的多站中长期径流滚动概率预测方法 | |
CN116151621A (zh) | 一种基于数据分析的大气污染治理风险检测*** | |
CN113836808A (zh) | 一种基于重污染特征约束的pm2.5深度学习预测方法 | |
CN117892165B (zh) | 一种基于灾情分析的低温灾害农业影响预测方法 | |
CN114295778B (zh) | 一种基于火电厂烟气排放质量评价方法及评价*** | |
CN111898093A (zh) | 适用于气象集合预报的多维后处理方法、装置以及*** | |
CN113033910B (zh) | 一种光伏发电功率的预测方法、存储介质及终端设备 | |
CN112307675A (zh) | 一种基于神经网络的温度敏感型负荷分离辨识方法及*** | |
CN117893057A (zh) | 一种基于深度学习的光伏组串积灰程度评价方法及*** | |
CN117200223A (zh) | 日前电力负荷预测方法和装置 | |
Libasin et al. | Evaluation of Single Missing Value Imputation Techniques for Incomplete Air Particulates Matter (PM 10) Data in Malaysia. | |
CN115511230B (zh) | 一种电能替代潜力分析预测方法 | |
CN111680888A (zh) | 基于rfid数据确定路网容量的方法 | |
CN113225391B (zh) | 基于滑动窗口异常检测的大气环境监测质量监控方法及计算设备 | |
CN110108609B (zh) | 一种基于卫星多通道的pm2.5浓度值分布情况获取方法及*** | |
CN115291302B (zh) | 基于频率与百分位匹配的逐小时降水订正预报方法及设备 | |
CN117192063B (zh) | 基于耦合卡尔曼滤波数据同化的水质预测方法及*** | |
CN116522166B (zh) | 一种基于大数据的空气污染监测*** | |
CN114971072B (zh) | 一种基于多因子相似分析的水文预报*** |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |