2 x* c# j7 O: \- G" o) x2.PM2.5的分布与演变及应急处理 / ?% O3 J, _( @- G# N % w6 R& u+ l+ k. X9 u2.1第一小问的模型建立与求解) p+ M" c. c# A4 x, [" C
3 Z% |& j- V' B$ ^" K& z+ e' X
该问题要求描述该地区内PM2.5的时空分布及其规律,并结合环境保护部新修订的《环境空气质量标准》分区进行污染评估。 " K% E- |) H$ I1 ~/ S: u5 Z, a1 I5 v& k( }7 y) `
原文思路如下: ) C# ^" C/ Z6 ~$ j. H3 U* N ( r& x" Y7 q# J0 z5 ~# @(1)绘制PM2.5浓度与时间的散点图,得出分布趋势; & ~. X& c1 k1 H1 O3 R% _6 e" T, {% T 7 A- [* N9 J7 Q$ A(2)以全市平均浓度代替各个监测点浓度刻画PM2.5时间分布;6 V% X$ E. A( V; Y R8 ~; T
. J/ q7 J$ W1 m i6 M& h(3)使用Shepard反距离插值法刻画PM2.5空间分布;2 G1 O. j2 h1 X# @4 B
# o- @) f& R+ Q5 h! w
(4)污染评估。 9 ]9 ]9 r' }9 j& \0 C C% ?5 {* p# Y, C* c7 Bshepard反距离插值结果 \1 M: n* B# n7 \# d4 X
& a2 I6 Q; F. X) c/ {加入人工神经网络训练的结果/ w2 @. R! C% _) H4 y
原文思路清晰,以PM2.5全市平均浓度随时间的变化刻画西安市PM2.5浓度的时间分布;通过Shepard插值求得西安市PM2.5浓度的空间分布。但是原文程序在计算PM2.5浓度空间分布时用的却是立方插值,并没有使用上文提到的Shepard反距离插值。$ y5 j6 B/ v7 ~9 Y
4 R' p& u. f* m4 U g) k" d
n8 ~* {1 S, L: ?* F' M2 x6 T2.2第二小问的模型建立与求解6 y, w( H, i c
问题要求建立能够刻画该地区PM2.5的发生和演变(扩散与衰减等)规律的数学模型,合理考虑风力、湿度等天气和季节因素的影响,并利用该地区的数据进行定量与定性分析。' d, m7 A- V% r
" }2 F/ R. n+ T原文思路如下: 2 N' d; u/ t) }0 M, C. N! u, T0 W, a4 |4 J2 m, q7 \3 {
(1) 利用该地区的数据进行PM2.5与风力、湿度等天气和季节因素间的多元回归分析,即建立PM2.5浓度的估算模型; * R o$ r- n* ^) i 9 i b" ~8 H: ?( T4 W. {0 U(2) 利用大气污染物的扩散原理建立PM2.5扩散的偏微分方程模型,主要参考《城市大气污染扩散监测模型的理论与试验研究》这篇论文。! s5 K8 M$ g4 ]
0 q8 G; O/ V8 X3 e8 Z+ t$ z9 z5 c
本人能力有限,反正我看不懂,对于该问题的求解,我主要参考了国防科技大学同学的一篇优秀论文,利用经典的对流扩散方程进行建模,并使用ComsolMultiphysics对问题进行仿真实验。2 O& `9 G) n# \& A
% Y7 ^" W2 g1 f3 j. D最终的对流扩散方程如下,等号左边后两项为对流项,等号右边前两项为扩散项,第三项为源项。 2 M% s. M9 G" {+ F X 9 c5 c, Z) z1 X# k% Q! N& R原文的方法对我而言晦涩难懂,对环境工程、污染物传播的相关知识要有良好的基础,超出了我专业能力范畴,所以我选择以国防科技大学的论文作为研究该问题的主要方法,国防科技大学的同学在论文中对所使用的方法进行了详细的阐述,通俗易懂,使我获益匪浅,并且该方法有良好的可扩展性强。南京航空航天大学的论文在进行PM2.5与风力、湿度等天气和季节因素间的多元回归分析时,提到“对于风速,根据风速与风级的对照表,由于基本都为小于等于三级风,若为平均值或等级值将为定常数不具有统计作用,对此我们另外在相关气息官网上查阅到的该地区风速情况的具体数据以及历史数据,对此为风速数据”,然而找遍全文我也没有看到风速历史数据的出处,对于一篇超过50页的科技论文而言,用一句话说明你论文中用到数据的来源应该不过分吧。 % I) a% g6 I- F2.3第三小问的模型建立与求解 + I! g+ K) x& o$ s( ^ 7 Y- R p4 C F% X& F! V% |题目要求以该地区PM2.5 监测数据最高的一天为例,在全地区PM2.5浓度最高点处的浓度增至2倍,持续2小时,利用建立的模型进行预测和评估。 $ t1 G* t }* u( n8 E使用Comsol Multiphysic进行无风情况下仿真,无风情况对流为零,以扩散为主,而PM2.5的扩散系数很小,所以呈现堆积现象。" ]4 b! L6 E. Q& R
使用Comsol Multiphysic进行有风情况下仿真! I+ L& _9 A9 Z
0 u& S0 g2 o& F) A
3.总结 - Q: k5 ?" s; x* b. z" F. b % N% L( @4 ?9 `$ E; [* ^自己在重现的过程中发现自己用具体方法实现可能还可以,但具体的数学逻辑推导可能就不能表达得很清楚,在数学逻辑写作方面需要向两篇优秀论文学习。- E8 e0 s; a) m, A9 N9 H: D
! M* h! e- r& Z5 I1 j
南京航空航天大学同学论文的问题: # E6 ~; }4 F5 s/ E+ j4 e& u: N4 r' o/ p& {4 Q* M- ^7 o
(1) 对于风速数据的出处没有说明; 6 ^ R( R' _+ n, j! g1 C" _ ! H$ F, u' s, v! p9 z, X' w4 i4 q(2) 论文中求PM2.5的空间分布使用的是Shepard插值法,但是程序中使用的仅仅是立方插值; * n4 E1 V) w9 K @8 R9 N6 [' K( {% v. g2 u1 z* q
(3) Lingo程序与原文逻辑不一致; ; r4 {- I8 |! j N 9 [- Y* t2 N9 _( W& {考虑到建模时间的紧迫性,在论文写作的过程中难免会有瑕疵,总的来说南航同学的论文还是非常优秀的,很值得我们吸收和借鉴。3 \/ T' {% N& g
8 U. S( T4 U4 M6 n) @/ G6 f4 c {* Q; c5 f0 p' D! { : H/ `3 V0 X8 v% @7 ?& M6 P5 W) t" y, Q