表1），Xa、Xc、Xp与Pi的计算关系如公式（2）所示。
fCi

Xa
Pi

12
Ci
XcCi
Xa
XaXc
XpXc
3
Ci

X
p
XpXc
CiXaXaCiXcXcCiXp
CiXp
（2）
根据式（2）计算得到污染指数值来对重金属的污染程度进行评价，其中Pi1为非
污染状态；1Pi2为轻污染状态；2Pi3为中度污染状态；Pi3为重污染状态。为了全面、综合地反映多种污染物的整体污染水平，因此需要一种同时考虑多种污
染物综合污染水平的多因子评价方法即将单因子污染指数按一定方法综合。常用的方法有内梅罗指数法它兼顾了单因子污染指数的平均值和最高值能较全面地反映环境质量而且可以突出污染较重的污染物的作用其计算公式如下：
P
1



i

CiSi
2


CiSi
2max
2
（3）
式中：P为土壤污染元素综合污染指数；Ci为土壤中污染元素i的实测值；Si为土壤中污染元素i的评价标准。同单因子污染评价程度一样，P1为非污染状态；1P2为
轻污染状态；2P3为中度污染状态；P3为重污染状态。
表1土壤污染评价标准的建议范围与数值
元素
污染积累起始值Xa106
中度污染起始值Xc106
重度污染起始值Xp106
Hg
015
030
15
Cr
90
250
400
Cd
02
03
10
Cu
35
50
400
Pb
35
250
500
As
15
30
60
f问题三：问题要求分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。针对这一问题，首先我们做出假设，每种重金属都是可以在外界因素下可以自由传播的，不存在重金属元素在一个地区堆积，只有传播能力的不同。由于海拔因素对其影响不大，故忽略海拔因素对其传播的影响。首先，利用matlab作图，画出以坐标x、y及重金属浓度大小的差值等高线图。这样，我们就可以看到在整个城市中某一重金属元素的大致分布。其次，由此图及第一问计算出各点的污染程度可以看出重金属的传播特征。其传播大致符合自然界规律，高浓度向低浓度传播，但由于其重金属元素的特殊性，其传播速度较慢，且不同重金属的传播速度不同。由此就可以得到污染源的特征就是其污染元素的含量与周围含量相比相对较大，即区域极值。
六结果分析
山区工业能源大都以煤、石油类为主，它们是环境中汞、铅、镉、铬、砷等重金属污染的主要来源。在采矿、选矿、冶炼、锻造、加工、运输等工业生产过程中会产生大量的重金属污染。在农业生产中，污水灌溉、农药、劣质化肥等的不合理使用是重金属污染的重要途径。以磷肥为例，生产磷肥的磷矿石成分复杂，含有较多的重金属如锌、铬、镍、铜、镉、铅等，因此如不r