.

以每平方公里为取样单位的话，只要位置波动值能控制在一千以内，那就说明我们的模型的稳定性是比较好的。我们分别改变模型中常数，

?t，?l，Vd和n的

值，以模型中的值为基准上下波动10%，并计算相应的污染源，（程序代码）和模型中一样，见下表： 变动的参数 参数较小10%情况参数不变的位置参数增加10%情况下的位置变化 变化 下的位置变化 300～630 0 400～700 ?t ?l Vd200～530 400～800 0 0 180～460 350～460 200～600 0 330～650 从上表中可以看出，数值波动都在1000一下，也也就是说，当参数变化是，对于单元素污染源的位置波动比较大，但是动控制在1000以内，也就是说，没有超过其极限值。所以验证了我们都模型的稳定性。

我们通过对单元素污染源点坐标周围的点用模型公式中进行计算，得出相应的模拟浓度，并进行Excel软件进行拟合，相应的点在表中对应相应的浓度，得到如下图形：

模型检验图18000160001400012000100008000600040002000013579111315171921232527293133353739实际线模型线n

通过图可以看出，峰值相差比较大，拟合效果不是很好，但他任然展示了相似的趋势，并且有很好的相似性。证明了我们的模型正确性。

七、问题四模型优缺点的分析及优化

模型优缺点的分析

内梅罗综合指数模型

对用于污染程度评价的模型一所使用的内梅罗指数法模型进行分析

优点：数学过程简便。物理概念清晰，评价方式简单便于决策。对数据的处理考

虑到了各个散点数据间的联系，加入权重进行综合内梅罗指数排名。

缺点：其描述的环境质量是非连续的，分级标准建立在二值逻辑基础上，它的截

..

.

然性和非连续性造成了相差很小的污染指数强度间可能会出于两种不同的等级。 污染源定位模型

优点：

? 模型有很好的实践性，而所给的算法几乎没有时间的复杂性。对于所

给问题的数据规模，我们采用格点搜索法求最优解。

? 模型得出了数据与计算值鱼很好的一致性，它是快速、有效和稳定的。 ? 至于对数据的简化计算，准确性并没有降低。 缺点：

? 如果考虑的区域比较大，就会有一些误差。

? 为了降低计算的复杂性，我们间滑落影响污染扩散的纵向地下水流这

会影响到结果的精确度，还有在地形方面只考虑了污染源的水平定位没有考虑其海拔高度的定位。

八、模型的推广与改进

8.1.1模型的推广：

由于我们的模型具有一定实用性和稳定性，所以可以将该模型推广到具有稳定地形的关于水污染的模型中，我们充分考虑到在模型中遇到的一些问题，设定参数在可变的情况下具有一定的稳定性。所以也可以应用到一些地形较稳定的地形中，也能达到良好的效果。

8.1.2模型的改进：

我们所建立的模型是将两种弥散系数a近似估计为一实数25ft，这是在外界条件相同的情况下的一种假定，而实际中不同地方土壤的PH值是不同的，特别是污染区。水-土壤系统pH升高能明显地降低各元素在土壤中的吸附,促进其在土壤中的迁移,且吸附常数(Kf)与土壤有机质含量、粘土含量呈正相关,而与土壤pH呈负相关.pH值高金属元素在土壤中淋溶贡献较大,且淋溶量随雨量的增大而增大.同时在土壤中的淋溶与土壤性质密切相关,有机质含量和粘粒含量较高的土壤对其的持留能力较强，由此我们可对模型进行优化。

PH对传播的影响主要可以体现在阻滞力系数Rd以及渗透力系数D上，对此我们假设

Rd?

并且令

1,D?PH （17） PHs,D?kPH (18) PH带入模型得到优化后的模型如下：

Rd?..

.

C(x,y,z,t)???r?VdtRd12??r?VdtRd12??sqC0VV?Vdexp(d)?exp(d)erfc?()??exp()erfc?()??8?nDrPH2kPH?2kPHDt?2kPHDt???2?2 (19)

在针对模型进行取样采点进行统计计算，即可得到s，k的值，获得后就是完善的模型.

九、参考文献

【1】孙树瑜，曾爱武．王树楹。等．规整填料塔中精馏过程的三维模拟lⅡ)：

模型的验证及液相分布和混台行为对精馏过程的影响 J] 化工学报 1998，49(5)：503—565

【2】夏星辉,陈静生.土壤重金属污染治理方法研究进展[J].环境科学,1997，

18(3)72-76

【3】郑州市郊区农业区划办公室．郑州市郊区土壤普查办公室．郑州市郊区壤

R1．1986，39．

【4】魏秀琴．郑州市东北郊污水灌区重金属元素对土壤的污染【J】．河南

地质，1998，16(3)：46-48．

【5】刘玉燕，刘敏，刘浩峰．城市土壤重金属污染特征分析【J】．土壤通报，

2006，37(1)：184一188．

【6】魏秀琴．郑州市东北郊污水灌区重金属元素对土壤的污染【J】．河南

地质，1998，16(3)：46-48． 【7】刑文训，谢金星，《现代优化算法》[M].北京：清华大学出版社，1999 【8】[1] 孙铁珩.污染生态学[M].北京：科学出版社，2002：18-24.

【9】[2] 鲍桐，廉梅花，孙丽娜，等.重金属污染土壤植物修复研究进展[J].

..

.

生态环境，2008，l 7（2）：858-865.

【10】 夏家淇，骆永明.我国土壤环境质量研究几个值得探讨的问题[J].

生态与农村环境学报，2007，23（1）：1-6.

【11】 阎伍玖. 芜湖市城市郊区土壤重金属的初步研究[J]. 环境科学学

报, 1999, 19(3): 339 - 41.

【12】 张中一, 朱长会. 南京市郊菜地土壤重金属污染状况[J]. 南京农专

学报, 1995, 26(4): 6 - 11.

【13】 张庆利, 史学正, 黄标, 等. 南京城郊蔬菜基地土壤有效态铅、锌、铜

和镉的空间分异及其驱动因子研究[J]. 土壤, 2005, 37(1): 41 - 47. 【14】 李静,谢正苗,徐建明, 等.杭州市郊蔬菜地土壤重金属环境质量评价[J].

生态环境, 2003, 12(3): 277 - 280.

【15】 陈晶中, 陈杰, 谢学俭, 等. 北京城市边缘区土壤重金属污染物分布特

征[J]. 土壤学报, 2005, 42(1): 149 - 152.

【16】潘剑君主编.土壤资源调查与评价[M].中国农业出版社, 2004

【17】 丁爱芳, 潘根兴. 南京城郊零散菜地土壤与蔬菜重金属含量及健康风险

分析[J]. 生态环境, 2003, 12(4): 409 - 411

【18】 陈芳，邹修明.基于二次曲面拟合的插值算法[J].计算机应用与软件，

2006，（ 11）：105-107.

【19】〔4〕华中工学院数学，软件教研室.算法语言·计算方法[M].高等教育出

版社，1978. 【20】〔3〕崔国华.计算方法[M].华中科技大学出版社，1996.

十、附表及附录

附表一：取样点位置及其所属功能区 编号 x(m) y(m) 海拔(m) 1 74 781 5 2 1373 731 11 3 4 0 1787 4 5 1049 2127 12 6 1647 2728 6 7 2883 3617 15 8 2383 3692 7 9 2708 2295 22 1 11 4233 895 6 12 4043 1895 14 14 3526 4357 7 ..

功能区 4 4 2 4 2 4 2 4 5 1 1 4 4 .

16 17 19 2 21 22 23 24 25 26 27 28 29 3 3 32 33 34 35 36 37 38 39 4 41 42 43 44 45 46 47 48 49 5 51 52 53 54 55 56 57 ..

4777 5868 5481 2486 3299 3573 4741 5375 5635 5394 5291 4742 7004 7304 7048 8180 9328 9090 8049 8077 6869 7056 7747 8457 9460 9062 9319 1 1 11702 11730 11482 1 1 11678 11902 4897 4904 6004 5999 6018 6213 6434 8643 7965 8631 7349 7293 6226 5230 4600 4496 4311 5365 5439 6401 7286 8348 8260 8991 8311 7639 6799 4480 5532 6354 8618 7709 8 16 0 2 4 5 5 15 29 12 10 9 11 10 24 15 24 20 18 29 18 37 49 21 45 45 49 5 71 54 61 17 30 1 4 1 4 1 4 1 5 1 4 4 4 2 4 1 4 4 1 4 4 1 4 4 1 4 4 4 4 4 4 4 4 4 4 4