TA的每日心情 | 奋斗 2024-7-1 22:21 |
|---|
签到天数: 2014 天 [LV.Master]伴坛终老
- 自我介绍
- 数学中国站长
 群组: 数学建模培训课堂1 群组: 数学中国美赛辅助报名 群组: Matlab讨论组 群组: 2013认证赛A题讨论群组 群组: 2013认证赛C题讨论群组 |
摘 要
0 l$ Q* k# [* V2 l3 k6 Y本文主要探讨的是PM2.5 扩散、衰减模式的问题,根据该模式分析探究9 s+ p* ]7 @5 p# K" ~
PM2.5 的危机治理与后5 年的治理问题。建立了PM2.5 与其它污染物之间的多4 b y8 G- t! Y0 m% k
元非线性对数模型;在静态、风力和湿度等因素下,分别探究污染物颗粒的运动
: R9 `6 U! C) Y8 @模式,并建立了PM2.5 扩散演变模型;在污染物浓度突变的情况下,依据该模
0 w' _) r1 W8 C4 q/ m8 K. K4 n型得出不同区域污染物的浓度,最后确定安全区域的范围。建立综合费用和专项
& K2 |( T9 W3 H费用的多目标优化模型,并运用系统动力学理论对目标值进一步优化。* P$ r8 Z6 z+ d# S, M
针对问题一,首先,运用主成分分析法,按照方差贡献率的大小剔除臭氧这9 y; m, {' s( M m+ |/ K1 y
个指标;其次,运用SPSS 软件分析剩余指标之间的相关性及独立性,并建立了5 x9 G/ ~( h: k( F) b" `8 {9 G" b
PM2.5 与其它污染物之间的多元非线性对数模型,得出西安市的拟合优度
2 t3 R# s6 f9 e' y2 R 0.90;最后,搜集西安市的相对湿度数据,运用该指标对模型进行再度优' o3 k- h4 g# |( Z$ U
化,优化后的拟合优度2 R 0.94,因而相对湿度是PM2.5 影响因素。
! _ X3 b! a8 o0 b7 F针对问题二:(1)运用统计学原理分析13 个监测点PM2.5 的浓度,描绘( f- I2 f% c1 d: `0 ^
了西安市PM2.5 的空间分布云图。同时添加时间因素,探究PM2.5 颗粒四维空
' T, q2 Y1 M# T0 z2 D间分布情况,得出采集点之间的PM2.5 具有较高的协同性。' {$ H* W, `% O* c) ]
(2)分析了静态下PM2.5 粒子受力,漂移模式,通过结合风速、湿度、大4 ~3 c# \! y- X
气稳定度等季节性因素从点源、面源两方面分析了PM2.5 扩散模式;建立了
' Q& K4 ] }8 M2 ^6 }' iPM2.5 点源和面源扩散的偏微分方程模型;通过利用P-G 曲线近似法,布里吉
, M; N, Y$ }8 Y! a* r斯扩散参数以及现有数据对季节参数进行求解,得出PM2.5 扩散衰减模型。计
; y% B" Y. R5 [算结果与西安市地理位置和提供数据相吻合,说明模型所刻画传播衰减模式与事( _$ m' Z# c) N J7 H% b! O% Z
实相符。$ Z( b/ f o5 k+ s1 `9 M- F h5 T) _
(3)通过第2 小问所得的PM2.5 点源扩散模型与PM2.5 面源扩散模型,以
6 e O+ y& w3 a高压开关厂为参考点,在3 级北风状态下,运用matlab 软件仿真模拟出的点源
. K7 V8 [& z5 m0 H& r4 u2
5 u; e4 D0 ~5 {/ _' W8 W与面源扩散情况,其结果展示如下表:5 O0 _& U0 C/ [4 v
扩散方向/ a# l8 G: e! ^. I) U% N
向东
3 o' h1 X5 j7 a+ u向西" h5 k" m. i" r; S! }5 T
向南" x% o! P' J. B: m, M, d
向北
2 B. @$ ` \/ d/ w( D扩散
. A3 G- o7 \9 z- J+ g7 M距离(m), ^2 X1 O$ [9 O2 l& P
点源模型
( L. V6 Y7 Z5 ~0 K9 r: f S50+ X3 `$ \( d7 b9 [, z/ c
50% K% i: C4 H1 D3 ~% T( V
200
! j2 V% C# Y7 o1 ]& l' x% A108 n" f) c; H' J7 \0 _4 t) w$ c
面源模型
# z0 t! n# j3 E9 C0 |# Y500
( \% _7 n# ]6 I$ G+ ~4 H' E2 B5 \500
& C; s/ u& T9 {3000
8 `, U# r* s8 d2009 o4 `1 D* s9 ^ W0 r
(4)结合西安市各个地区的地理位置和天气、气候等条件,建立了各个区之间的PM2.5扩散分析体系,利用西安市2013年1月8日—2月8日的数据,通过模型求解出各个区之间的PM2.5的相互扩散量,然后计算仿真出各个区的PM2.5的浓度,通过与原始值进行对比,发现模型所得结果与实际相差在10%之内,说明模型可信。仿真与原始对比如下:+ A- Q) J, p B6 G; M
高压开关厂
% b/ }( x- I2 d f9 S兴庆小区
& I( y2 h/ @& h$ z日期
3 J9 r- T8 D- m真实值6 m) D. f+ k) i: Z6 a. X- j0 g
计算结果
! p. M2 r, z2 c$ d# \8 A7 U R( d误差率9 s1 v7 I. N/ W% K
真实值
0 B( H& n' a% {+ C7 S计算结果
1 F2 |8 y* _! n& H; l1 ]误差率 2013-1-8! j' v1 w, _8 r2 a( E
383& f' } d% [; B$ Z: q5 b6 \
356.7054
2 R+ w" m, _/ P0 m! }-6.87%
* `/ o& [& L5 p! x+ z373& i2 e) w, Q. n7 W" o O8 `
381.17837 @, b- m" c+ A+ ?; O% o" ] s
2.19% 2013-1-9
+ g. Q( f" W9 s: c216
5 X7 B) u8 a2 Y# t7 E+ Z2 {211.39
+ n' m U. x! ~9 V8 o-2.13%
3 B( J5 f. u4 O2 w5 ?$ q% c236
3 A" a1 Z7 D8 w217.71479 Y n* ~6 H0 t! Y& |
-7.75%, z0 _2 P7 j& _8 [' y* w. \$ a* h
针对问题三:基于系统动力学理论,考虑治理效果,建立了系统动力学多目标复合治理的最优化模型。利用贝叶斯支持向量机方法对武汉市基本面数据进行宏观预测,对PM2.5进行系统性预测,并且仿真求解出PM由280单位到35单位的五年治理方法,结果表明将综合治理与专项治理结合时治理效果最好。其最优相结合治理计划为: x, H( N0 q! m
年份5 U* T9 H2 f3 J8 k; Q" x
2013
# W" ^. _# s1 s5 u0 H+ R' C2014, {4 [6 ? E x! u
2015& I* W7 L/ l3 u3 H; N
2016 [0 b$ B5 i! j3 T3 M
2017
# T6 c5 U7 z( s* y! U f综合6 y% u5 M8 z& i! o
治理& Q! a9 E# m, G6 s6 b' H
投入费用(百万)
" r- J. Y6 V) c8 r51
) ]& E6 b+ Y/ S! O& Z" S42$ i% b' Z: J" p: ]- u+ U% k% K9 {
32( Z; z N. y5 n/ F# A
222 r* m1 {: J2 T D1 o% L
12: e* o7 w9 T4 H# H7 h
PM2.5减少浓度
; h0 x$ e. `( m5 N. E+ f! `4.5! @. N/ n3 a+ _2 `7 x& J$ E+ f M
19.3
: \9 K- ]4 H/ e2 l+ k8 Z346 V5 b% e" G$ J; u3 s- z! z K
48.76 C! E# g& \1 i' t# B0 ]
63.5% [! K9 A7 Q0 K0 h# V6 i
专项
R& q, O+ w6 A6 E* J& O治理
, \4 r* I1 |( W* j0 @5 |2 A投入费用(百万)
5 a% v: [# b/ c7 C0 y20* {! h" b2 i1 f4 N
214 M* v- X5 }" w7 o
19" l# @1 z( b: P9 B! w+ ~2 `' P3 V
20
6 u1 Q$ W- h. X; A7 Q1 k18& h9 t$ v; _9 H' X2 h, c5 ^
PM2.5减少浓度( O7 D; Z2 f& m+ G
285 h& n. \( e4 x0 g! g' [+ J3 g
216 p/ ]# u9 l" l3 \3 }% {, `- g
15
8 u. m- H+ ~+ K K; ~! `& S8.5! {. B1 L0 o0 ~8 `7 @
2
3 N! U; ]. X5 \7 W- X9 b最后结合本文研究结果,对研究实施进行总结撰写了一份研究试验报告。
3 c) x" u! a+ E1 I' L本文创新点在于,建立了基于贝叶斯理论的支持向量机方法和基于系统动力学的多目标治理模型。/ { `8 t! ]% ]3 \9 j+ U
关键词:多元非线性对数模型,扩散衰减微分模型,系统动力学理论4 F) t3 G. ]5 Q8 H5 ^. g
7 a8 a& a* n$ G5 Q2 X |
zan
|