建立
在海拔最高与最低相差308米，且坡度高低起伏的总体成坡型分布的城区，我们不应该武断的忽略海拔，特别是与海拔相关的因素坡度对此问题引起的微妙的影响
为此本实验引出梯度场分析法：在已知函数某元素的空间分布散点图，和地势图后，可模拟出该二维浓度图的梯度场如下图所示：
图7同理可模拟出该城区总的地势分布的梯度：
f图8据此不难求得与之差，即为平面扩散梯度，但随着实验的进行，新的问题出现了，此时只知道的方向，并不知道两个矢量的数值关系，因此若以E为一给定矢量，设E2为E的k倍，至此问题并没有结束，k是变数还是常熟此时还不得而知，本实验的一大特点就是用计算机进行了半定量的模拟实验，将k值设定为一系列可供选择的参考值，逐一带入矢量差值表达式EE1观察比对做差后所形成的新的梯度场，在此仅将该过程简记如下：
图9
f图10
图11随着k的增大，矢量做差后的梯度场图像逐渐逼近一个成中心辐射状分布的梯度图，易推知该梯度场的等势线接近一族同心圆，最后用matlab就此时的数据反向
f画其等势线，证实猜想，由此可知，原扩散场为一回转场，且场源就是污染源将上图存为jpg格式图像后用matlab中imgread函数定出场源的X，Y坐标。以上基于matlab的计算机模拟实验的意义在于，将一个受空间的多因素影响的不同高度差下的扩散模型简化为一同一平面的空间扩散模型。
因本问题属于扩散问题，可借助于菲克定律建立数学模型
进一步分析问题可发现，固体的扩散缓慢近于稳态模型
5模型的求解：
pλ，ρ
式中λ为重金属污染扩散系数，α为实系数，p为某元素的密度，ρ为单位面积上的质
量。
通过MATLAB数据拟合可以得出λ，分别为：λ30291，00001。模型既为：
p30291
，ρ
传统指数评价模型
对土壤环境质量进行评价时所采取的评价方法多种多样主要有污染指数法、污染程
度法T值分级法、基准分级法、密切值法和模糊数学综合评判法、灰色聚类法等。文
献9给出了采用污染指数法和GIS空间分析技术对德兴地区土壤的环境质量进行现状
评价。由于地域差异等各种因素的影响污染指数法没有统一的评价标准。通常进行土壤
环境质量评价时，参照中国土壤环境质量标准（GB156181995）来确定研究的评价标准
（表1）。具体的单因子污染指数法如下：
Pi

CiSi
（1）
式中：Ci为土壤中污染中污染元素i的实测值；Si为土壤中污染元素i的评价标准。采
用Xa、Xc、Xp分别代表土壤污染积累起始值、中度污染起始值和重度污染起始值（如
r