TA的每日心情 | 奋斗 2024-7-1 22:21 |
|---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
摘 要:
: R* K+ b0 D! h+ Z6 @9 r. g本文针对日益严重的PM2.5 大气污染问题,基于现有数据与相关研究,采* P6 y8 q2 g" Y# z. d' f% W6 ^
用相关分析、逐步回归和有限元等方法,对AQI 指标之间的相关性进行了定量
+ o( _8 [; [+ T- g q& i- f2 `分析,针对PM2.5 成因、时空分布、演变规律和应急处理建立了数学模型,结
& B0 ?, O) D9 y合气象理论知识,进行了计算结果的定性与定量分析,并以武汉市为例给出了合- p U7 P: N( I
理的PM2.5 五年治理方案。
' Q# r) J3 E- {: g8 ~; B& i问题一:基于武汉与西安市现有的空气质量监测AQI 指标数据,应用相关6 N' m4 }: {* t
分析,分别求得两城市的6 个AQI 主要指标的相关矩阵,对两组数据中PM2.59 ]% y2 p7 n) u D- ?
与另5 个AQI指标的相关性进行了定性与定量分析,得出了一致的结论,即PM2.5$ S% U9 k# o. r
与另5 个主要指标具有相关性,并与SO2、NO2、PM10 及CO 显著线性相关(相5 R! l( c; w# Q( h
关系数均明显大于0.7)。
; W' f1 s j w9 \. h问题二:首先应用回归分析方法,针对2013 年仅有的AQI 完整数据,完成, ]% J3 s( |3 p
了PM2.5 关于API 中前3 个指标数据的二次模型拟合,拟合统计量表明该拟合
! c6 K7 o# _/ h: m% p是高度显著的。并基于大量API 历史数据,计算每组API 数据对应的PM2.5 估8 V/ g/ ?" E; X: f1 }$ \2 I' R" O
算值。而后,利用PM2.5 实测数据分析PM2.5 浓度的时空分布特征,进而结合
! e8 ? k! T, E, B. a) ]- R环境保护部新修订《环境空气质量标准》,对该地区分区进行了污染评估。综合
" K& s+ m- c9 E多种气象因素,基于多元回归建立关于PM2.5 演变规律的数学模型。结果表明,; _( I7 ?0 J7 V1 M7 k" J5 P7 u+ d
模型估计值与实测数据差值较小,模型参数与实际情况基本一致,验证了模型的. u# T' J' Y( e* f: P. X7 Y
合理性和有效性。针对浓度扩散问题,本文在适当简化实际问题的基础上,建立8 B3 H* O1 p' M8 ^+ @
了描述PM2.5 浓度扩散现象的偏微分方程,通过采用商用有限元软件Comsol
* D1 U2 |, ~9 ]- IMultiphysics,建立了针对该问题的扩散预测与评估计算平台,并对PM2.5 扩散
5 c% [% A* O" a# o( R- u现象进行了分析研究,得到了浓度突增等突发事件对PM2.5 分布影响的几个重
+ _$ ?7 h/ c# z0 _26 d- Z/ s) R! R3 o4 Z% X
要规律。主要包括速度和扩散系数的取值对扩散趋势的影响。并在此基础上,基于有限元仿真平台,运用自行推导的近似公式得到了扩散系数的估计值,对测点浓度突增现象进行了仿真计算,得到了PM2.5污染预测和评估结果。最后,从模型假设条件和模型主要结论两方面,分析了模型的合理性,并在前文研究的基础上对PM2.5的成因和演变规律进行了较系统的分析和探究。4 t% L q5 c% I, u) S6 S( X
问题三:首先应用逐步回归分析方法获得了PM2.5与3个既与综合治理又与专项治理有关的因素SO2、NO2、PM10的最优二次回归方程,在将治理基准年与五年后目标的PM2.5浓度折算为对应的AQI指标之后,结合数据分析计算出了基准年的PM2.5、SO2、NO2、PM10这4个基准指标值。然后按照每年这4个指标变化最小但仍能满足治理结果约束的原则,建立了非线性规划模型,并通过将非线性约束条件简化为线性约束,将模型转化为标准的非线性规划问题,进而求出了可行的每年全年年终平均治理指标。最后,结合PM2.5综合治理与专项治理的费用经验公式,按照五年治理投入总经费最小的原则建立了一个标准的非线性规划模型,并求解出了最少的五年总投入经费及每年经费投入。对于指导PM2.5的治理工作具有一定的实际意义。
% J$ P2 U+ {* F关键词:PM2.5,相关分析,逐步回归,非线性规划,有限元
- e+ N2 ]+ j3 g( Q' i7 Z8 T: e$ k" m2 c, E5 H0 v7 @, ^* @3 S
|
zan
|