一种区域空气PM2.5浓度预测方法技术

本发明专利技术公开了一种区域空气PM2.5浓度预测方法，该方法首先通过历史数据来构造待训练的支持向量机回归模型的训练样本数据，然后通过训练样本数据得到训练的支持向量机回归模型，将该训练后的支持向量机回归模型作为PM2.5浓度预测模型；再将粒子群寻优算法与PM2.5浓度预测模型相结合，通过粒子群寻优算法不断寻优迭代，从而不断采用粒子的位置去重构PM2.5浓度预测模型的输入参数，直至迭代完成后得到粒子群最终的全局极性，采用粒子群最终的全局极值对应的粒子的位置去重构PM2.5浓度预测模型的输入参数，将该输入参数输入PM2.5浓度预测模型中得到的输出即为PM2.5浓度；优点是可降低PM2.5浓度预测模型输入参数的维数，提高PM2.5浓度预测准确率。

【技术实现步骤摘要】

本专利技术涉及一种PM2. 5浓度预测方法，尤其是涉及一种区域空气PM2. 5浓度预测 方法。
技术介绍
随着工业化的快速推进，大气雾霾现象越来越严重，PM2. 5是造成雾霾现象的主要 元凶之一，其颗粒粒径小，能长时间悬浮于空气中并传播，并可携带有毒有害物质进入呼吸 道和肺部，对人体健康造成直接威胁。为此，基于大气的历史监测数据，对PM2. 5进行浓度 预测具有重大意义。但PM2. 5浓度的影响因素复杂繁多，如：即时工业排放、大气中污染颗 粒物浓度、气象条件、季节及阳光辐射，对其进行浓度预测具有挑战性。 最早应用于空气污染物浓度预测的为多元回归模型，随着人工神经网络的发展， 不少学者研究用神经网络模型对空气污染物浓度进行了预测，所展现的效果均优于多元回 归模型。然而神经网络普遍存在由于大气污染物浓度的训练数据不充分而引起过拟合的问 题。支持向量机（Support Vector Machine，SVM)是20世纪90年代提出的一种新学习机， SVM因在解决小样本和非线性问题方面表现出良好的泛化能力，已在分类、时间预测、函数 估计等领域得到广泛应用。Sanchez等、Xue等、Nieto等提出利用支持向量机回归（Support Vector Regression, SVR)模型f (X) =〈co，X>+b，这里，c〇，b为待辨识的参数，X为模型 的输入，〈，> 表示点积运算。对空气污染物进行浓度预测的方法。无论是SVR模型还是人 工神经网络模型，对污染物浓度的预测都需考虑变量的时序特性，并在模型输入参数中考 虑时滞因子，以致模型的输入参数维数高，训练样本数据不足，最终导致模型的预测准确率 不高。为此，有学者提出一种基于遗传-神经网络的（GA-ANN)模型，GA-ANN利用神经网络 (ANN)模型对大气污染物浓度进行预测建模，用遗传算法（GA)根据预测性能指标对输入参 数进行0或1选择，这样虽然降低了输入参数的维数，一定程度上避免了模型的过拟合，但 同时也损失某些输入参数的有用信息，预测准确率仍然不高。
技术实现思路
本专利技术所要解决的技术问题是提供一种在降低预测模型输入参数维数的基础上， 可以提高PM2. 5浓度预测准确率的区域空气PM2. 5浓度预测方法。 本专利技术解决上述技术问题所采用的技术方案为：一种区域空气PM2. 5浓度预测方 法，包括以下步骤： ①构造待训练的支持向量机回归模型的训练样本数据： ①-1将当前时刻记为t时刻，t-n时刻表示当前时刻的前η个时亥I」，η = 1，2, 3, 4,…； 将t-Ι时刻测量得到的区域空气的大气能见度记SxN(t_l)，浊度记Sx z(t_l)，气 温记为Xt (t-ι)，气压记为Xp (t-ι)，风速记为Xw (t-ι); 将t_td时刻测量得到的区域空气中S02的浓度记为X SC]2 (t_td)，N02的浓度记 为xNci2(t-td)，PMlO的浓度记为xpnlQ(t-t d)，ΡΜ2· 5的浓度记为Xpn25(t-td)，03的浓度记为 x〇3 (t-td)，CO的浓度记为xco (t-td) ;td为时滞因子且t尸1，2,…，d，d为大于等于2的整 数；将t时刻预测的PM2. 5浓度记为ypni25(t); ①-2将待训练的支持向量机回归模型f (X) =〈 ω，X>+b的输入参数记为Xin⑴， 待训练的支持向量机回归模型的输出记为Yciut(T)，Xin(T)和Yciut(T)构成训练参数对，令 Xin (T) = γ_(Τ) =Xpni25 (T+l) (2) Xin(t-1) = 其中 T = t-2，t-3，t-4,…t-501，□为矩阵符号，Xin(T)为 IX (6Xd+5)的矩阵， 符号X为乘运算符号； ①-3获取t时刻之前区域空气中S02浓度、N02浓度、PMlO浓度、PM2. 5浓 度、03浓度、CO浓度、大气能见度、浊度、气温记、气压和风速的历史测量数据，得到 {Χιηα-2)，Υ_α-2)}，{Χιηα-3)，Υ_α-3)}，…，{Χιηα-501)，Υ_α-501)}共 500 组训练 参数对，将该500组训练参数对作为待训练的支持向量机回归模型的训练样本数据；②采 用训练样本数据对待训练的支持向量机回归模型f (X) =〈 ω，X>+b进行训练，得到支持向 量机回归模型的模型参数ω和b ;将得到的支持向量机回归模型的模型参数ω和b代入 支持向量机回归模型f (X) = <?，X>+b中得到训练后的支持向量机回归模型，将该训练后 的支持向量机回归模型作为PM2. 5浓度预测模型； ③初始化粒子群参数：在内随机产生种群规模为N的粒子群，N = 100,每 个粒子均含有位置和速度属性，随机初始化粒子群中各粒子的速度和位置，将第i个粒子 初始化后的当前位置记为W1(O),当前速度记为V1(O), i = 1，2，···，N，其中W1(O)和V1(O) 均为 I Xm 的矩阵，Wi (0) = ，i = 1，2,…，N，Wli (0)，W2i (0)，… ，wmi(〇) e ，m = 6Xd+5;将第i个粒子的在第k次迭代后的当前位置记为議， Wi (k) = ，第 k 次迭代后的当前速度记为 Vi (k)，Vi (k)= ，k = 1，2，…，kmax;k _为粒子群的总迭代次数； ④将粒子群中每个粒子的当前位置定义为PM2. 5浓度预测模型输入参数的权系 数矩阵；采用N个粒子的当前位置分别对PM2. 5浓度预测模型的输入参数Xin (t-Ι)进行重 构，由此得到初始状态下N个重构后的PM2. 5浓度预测模型的输入参数，将初始状态下第i 个重构后的PM2. 5浓度预测模型的输入参数记为A= (/ - U); CN 105184012 A 说明书 3/11 页 ⑤将粒子群中第i个粒子的个体极值的初始值记为，，f,,. = 6(0)，F1 (0)= Xp^t-D-y-do) |，其中I I为取绝对值符号，yciutl(t，〇)为At(/-U)作为输入参数输入 PM2. 5浓度预测模型中得到PM2. 5浓度预测模型的输出；将粒子群的全局极值的初始值记 为G。，（:广=min 丨f:(0),/ = 1,2,…翁}; ?-σι?τι ι-min 1 ⑥采用粒子群每次迭代后更新后的位置对重构后的PM2. 5浓度预测模型的输入 参数进行更新，粒子群第k次迭代后重构后的PM2. 5浓度预测模型的输入参数进行第k次 更新，将第k次更新后的第i个重构后的PM2. 5浓度预测模型的输入参数记为J^(f-l,i): 将(/-U)作为输入参数输入PM2. 5浓度预测模型中，得到PM2. 5浓度预测模型 的输出，记为;将凡代入公式Fi (k) = I Xpm25 (t-1) Iciuti (t, k) I中，计算得到第 k次更新后第i个粒子对应的适应度值F1 (k) ; I I为取绝对值符号；第k次迭代后粒子群中 第i个粒子的个体极值记为Pi, , Zy',,, =mmif:;(./)，./ = 0，l,2,·· ,将第k次迭代后第i个 粒子的个体极值对应的粒子位置记为W1 (a)，a为大于等本文档来自技高网...

【技术保护点】
一种区域空气PM2.5浓度预测方法，其特征在于包括以下步骤：①构造待训练的支持向量机回归模型的训练样本数据：①‑1将当前时刻记为t时刻，t‑n时刻表示当前时刻的前n个时刻，n＝1,2,3,4,…；将t‑1时刻测量得到的区域空气的大气能见度记为xN(t‑1)，浊度记为xZ(t‑1)，气温记为xT(t‑1)，气压记为xp(t‑1)，风速记为xw(t‑1)；将t‑td时刻测量得到的区域空气中SO2的浓度记为xso2(t‑td)，NO2的浓度记为xNo2(t‑td)，PM10的浓度记为xpm10(t‑td)，PM2.5的浓度记为xpm25(t‑td)，O3的浓度记为xo3(t‑td)，CO的浓度记为xco(t‑td)；td为时滞因子且td＝1,2,…,d，d为大于等于2的整数；将t时刻预测的PM2.5浓度记为ypm25(t)；①‑2将待训练的支持向量机模型f(X)＝<ω,X>+b的输入参数记为Xin(T)，待训练的支持向量机回归模型的输出记为Yout(T)，Xin(T)和Yout(T)构成训练参数对，令Xin(T)＝[xso2(T),xso2(T‑1),…,xso2(T‑d),xno2(T),xno2(T‑1),…,xno2(T‑d),xpm10(T),xpm10(T‑1),…,xpm10(T‑d),xpm25(T),xpm25(T‑1),…,xpm25(T‑d),  (1)xO3(T),xO3(T‑1),…,xO3(T‑d),xCO(T),xCO(T‑1),…,xCO(T‑d),xN(t‑1),xZ(t‑1),xT(t‑1),xP(t‑1),xw(t‑1)]Yout(T)＝xpm25(T+1)  (2)Xin(t‑1)＝[xso2(t‑1),xso2(t‑2),…,xso2(t‑1‑d),xno2(t‑1),xno2(t‑2),…,xno2(t‑1‑d),xpm10(t‑1),xpm10(t‑2),…,xpm10(t‑1‑d),xpm25(t‑1),xpm25(t‑2),…,xpm25(t‑1‑d),  (3)xO3(t‑1),xO3(t‑2),…,xO3(t‑1‑d),xCO(t‑1),xCO(t‑2),…,xCO(t‑1‑d),xN(t‑1),xZ(t‑1),xT(t‑1),xP(t‑1),xw(t‑1)]其中T＝t‑2，t‑3，t‑4，…t‑501，[ ]为矩阵符号，Xin(T)为1×(6×d+5)的矩阵，符号×为乘运算符号；①‑3获取t时刻之前区域空气中SO2浓度、NO2浓度、PM10浓度、PM2.5浓度、O3浓度、CO浓度、大气能见度、浊度、气温记、气压和风速的历史测量数据，得到{Xin(t‑2)，Yout(t‑2)}，{Xin(t‑3)，Yout(t‑3)}，…，{Xin(t‑501)，Yout(t‑501)}共500组训练参数对，将该500组训练参数对作为待训练的支持向量机回归模型的训练样本数据；②采用训练样本数据对待训练的支持向量机回归模型f(X)＝<ω,X>+b进行训练，得到支持向量机回归模型的模型参数ω和b；将得到的支持向量机回归模型的模型参数ω和b代入支持向量机回归模型f(X)＝<ω,X>+b中得到训练后的支持向量机回归模型，将该训练后的支持向量机回归模型作为PM2.5浓度预测模型；③初始化粒子群参数：在[0,1]内随机产生种群规模为N的粒子群，N＝100，每个粒子均含有位置和速度属性，随机初始化粒子群中各粒子的速度和位置，将第i个粒子初始化后的当前位置记为Wi(0)，当前速度记为Vi(0)，i＝1,2,…,N，其中Wi(0)和Vi(0)均为1×m的矩阵，Wi(0)＝[w1i(0),w2i(0),…,wmi(0)],i＝1,2,…,N，w1i(0),w2i(0),…,wmi(0)∈[0,1]，Vi(0)＝[v1i(0),v2i(0),…,vmi(0)],i＝1,2,…,N，v1i(0),v2i(0),…,vmi(0)∈[0,1]，m＝6×d+5；将第i个粒子的在第k次迭代后的当前位置记为Wi(k)，Wi(k)＝[w1i(k),w2i(k),…,wmi(k)]，第k次迭代后的当前速度记为Vi(k)，Vi(k)＝[v1i(k),v2i(k),…,vmi(k)]，k＝1，2，…，kmax；kmax为粒子群的总迭代次数；④将粒子群中每个粒子的当前位置定义为PM2.5浓度预测模型输入参数的权系数矩阵；采用N个粒子的当前位置分别对PM2.5浓度预测模型的输入参数Xin(t‑1)进行重构，由此得到初始状态下N个重构后的PM2.5浓度预测模型的输入参数，将初始状态下第i个重构后的PM2.5浓度预测模型的输入参数记为Xin0(t-1,i)=[xso2(t-1)×w1i(0...

【技术特征摘要】

【专利技术属性】
技术研发人员：史旭华，俞杰，童楚东，傅晓钦，汪伟峰，蓝艇，杨忠，李海琴，陈煜琛，
申请(专利权)人：宁波大学，宁波市环境监测中心，
类型：发明
国别省市：浙江;33