rj — 某地区第j个重金属污染物的权重值,且 ?Cj?1，?rj?1

j?1j?1nn将各金属因子的实测浓度值、毒性系数和选定的评价标准分别代入上式,可得到各金属因子的权重值,由此组成某个区域各金属因子的权重向量:

R = [ r1、r2 , r3 , ?,rn ] (9)

将权重向量R和隶属度矩阵D带入模型（4）可得区域金属对评价等级的隶属度N

C、隶属度函数的确定 ：

本文所选用的八种元素，选用国家《土壤环境质量标准》[ 18 ] ( GB15618 -

1995)进行土壤环境评价,土壤环境质量标准见表（4）。根据表（1）和表（5）的数据，利用公式（5）计算各重金属元素对应于各土壤重金属环境质量等级的隶属函数,得到 关系模糊矩阵。

表（4） 国家土壤环境质量标准

一级 二级 三级 自然景值 PH﹤6.5 PH=6.5～7.5 PH ＞7.5 PH＞6.5 As (mg/kg) 15 30 25 20 30 Cd (mg/kg) 0.20 0.30 0.30 0.60 1.00 Cr (mg/kg) 90 150 200 250 300 Cu (mg/kg) 35 50 100 100 400 Hg (mg/kg) 0.15 0.3 0.5 1.0 1.5 Ni (mg/kg) 40 40 50 60 200 Pb (mg/kg) 35 250 300 350 500 Zn (mg/kg) 100 200 250 300 500 根据Hakanson制定的标准化重金属毒性响应系数,分别对各重金属对生物的毒性级别指数赋值：Zn=l， Cr=2，Cu=5，Pb=5，Ni=5， As=10，Cd=30，Hg=40，指数越大，代表毒性越大。按照（7）、（8）两式，将表（1）的数据和毒性指数代入计算，得各区域各重金属权重值，其中

R = [ r1、r2 , r3 , ?,r8 ]= [0. 0962 、0. 2838、 0. 0680、 0. 0585、 0. 3723、 0. 0211、0.0563、0.0538] （10） 结合公式（4）、（6）、（10）及表（1）可得： N1=2.5158 N2=8.9537

N3= 1.0546 N4=6.4033 N5=2.5056

9

所以：N2＞N4＞N1＞N5＞N3，即，工业区>交通区>生活区>公园绿地区>山区

结合1）、2）两个模型可以得到该城区各区域重金属的污染程度为工业区>交通区>生活区>公园绿地区>山区

二．重金属污染成因：

根据上述两个模型的结合与对比，我们可以得出下面的关系式： 工业区>交通区>生活区>公园绿地区>山区

1.从整体的角度分析：

工业区污染最大，占整体区域的36%。分析其污染的主要来源有金属矿山的开采、冶炼、重金属尾矿、冶炼废渣和矿渣堆放等。而且其扩散和传播污染的途径也是非常广泛的。

交通区重金属的污染程度相对也较大，其污染主要体现为大气污染，来自工业生产排放的废气、汽车尾气排放生的大量含重金属的有害气体和粉尘。 生活区重金属污染居中，重金属的主要来源有生活垃圾、污水、医疗垃圾和细菌等。 公园绿地区重金属污染程度较轻，其主要来源有农业农药和化肥。但是，植物有吸收金属矿物的作用，相对减轻了重金属的危害。

山区重金属污染最轻，只占到整体区域的8%，表明不仅重金属的来源相对较少，而且，重金属污染与重力和海拔也有一定关系，山区中植物和一些微生物都有净化的作用。

2.从局部的角度分析：

从表（3）可以看出八种重金属的主要分布区域，再根据这些区域的区别可以分析出各金属污染物的成因。其主要来源如下表所示： 污染物 来源 砷As 矿山开采、冶炼、加工排放的三废、飘尘、农药等 镉Cd 电镀、塑料电池电子、染料化工工业废水 铬Cr 电镀、染料化工工业废水、、飘尘、矿山开采、冶炼、主要来源是汽车燃料和轮胎磨损。 铜Cu 电镀工业废水、农药 汞Hg 塑料电池电子工业废水、矿山开采、冶炼、加工排放的三废 镍Ni 电镀工业废水、矿山开采、冶炼、加工排放的三废 铅Pb 电镀工业废水、矿山开采、冶炼、加工排放的三废、汽车尾气 锌Zn 电镀工业废水 三、重金属污染物传播特性及污染源位置确定模型：

影响重金属污染物传播的因素有很多，如：地质条件（如海拔，地壳运动，土质情况等）、气候条件（如降雨量，气温，季风等）、工业发展情况、人口密度、河流流向等。在这里，我们由于条件有限我们仅考虑地势、河流等。

重金属污染物的传播有许多种，但主要有三种方式,即：水体传播、土壤传播、大

10

气沉降。在此，我们将模型理想化，假设只有这三种方式，分别对其建立模型，之后再将三种模型进行时间和空间上的几何叠加。

1.水体传播模型的建立及求解：

在此模型中，我们将问题理想化，认为水体传播仅仅是靠河流传播的。假设河流没有多大的改变。

在河流中，重金属污染物部分以可溶状态存在，随水流运动，部分以沙粒为载体，随沙运动，两者不断地发生化学和生物化学反应，一般情况下，泥沙吸附—解吸是重金属污染物从不饱和液相转移到固相的主要途径。因此，在考虑重金属在水环境中的物理迁移过程下，我们用水力学和泥沙运动学的原理处理。

假设，泥沙浓度和水相金属离子浓度反应分布均匀。在天然河流中, 水环境化学条件对某一河段、某一时段基本不变,变化较大并难处理的是水流泥沙因素。

泥沙颗粒在水溶液中对重金属污染物的吸附遵循兰格茂尔吸附运动学方程

dN?K1c(b?N)?K2N，将其对时间的偏微分 ，即对均匀沙有： dt?N?k1c(b?N)?k2N （11） ?t水流泥沙运动数学模拟研究已有较长的历史，在水流泥沙运动确定后, 根据建立的重金属迁移转化方程, 吸附动力学方程及边界、初始条件就可确定河流重金属污染物迁移转化过程, 构成河流一维重金属污染物迁移转化数学模型的基本方程如下：

??y水流连续方程 (Bhu)?B?0 （12）

?x?t?u?u?yu2水流运动方程 ?u?g?g2?0 （13）

?t?x?tc1k??y(Bhu)?B(hs)???B?(s?s0) （14） ?x?t?y??y河床变形方程 (Bhu)?B(hs)???B0?0 （15）

?x?t?t泥沙连续方程

有A?Bh,A3?By0及L=2h?B，则

?A1?N?c?c1??c1?s?N1???N?u?(AE1')???N13??N3L???s1?us1?E1'?（16） ?t?xA?x?xA?tA?x?x?x???t 假定悬浮泥沙的吸附量N1?N1?N1?N1'其中，N=N1'=0，对（11）式全断面积分得：

?N1?k1sc(bs?N1)?k2sN1 （17） ?t 解决重金属迁移转换模型可将（12）至（17）式简化成：

2?y0?N?Nv3m?c?c'?ch=常数，u?常数，?0,s??k(),?u?E12??s1?us1,?tgkw?t?x?x?t?x ?N1?k1c(b?N1)?k2N1?x 11

根据静态吸附的结果，取b=0.543g/kg,k1=0.0076 1/(ppm·s)，

k2=0.00084 1/s ,重金属污染物的初始条件为：

cx?0?c0?(x)式中：c0?1ppm?1?0x?0,Ntx?0?(x)=?t?0?0

边界条件为：

cx?0,?1?c0?(t),?(t)=??0t?0?2c,t?0?x2?0,N1x?Lx?0?0

在相同条件下，污染物的运动遵循以下方程：

2?c?c'?c?u?E12?0,ct?0?c0?(x),cx?0?c0?(t),?2c?x2?t?x?xx?L?0。用向后显式差分求解，结

果如下图2:

由图2可知：在不同时刻, 水相重金属浓度由于泥沙吸附, 峰值明显减少, 而且由于泥沙吸附需要时间, 因此河段上泥沙吸附量的峰值总是滞后于水相浓度的峰值。

从3中可以看出：污染物不受泥沙运动影响, 在输移过程中, 由于纵向综合紊动扩散作用, 峰型变化相同条件下重金属污染物输移, 尽管图形相似, 但重金属污染物的峰值因泥沙吸附而下降很多, 特别是后半部减小更为明显, 且随时间增加明显加大。

所以，重金属在水体中的传播特征是随着离污染源的距离增加而减少。速度也随之减小。

2.土壤传播模型的建立及求解：

重金属在土壤中的传播的影响因素也很多，不能单纯的看作是线性变化，在此，我们基于重金属在土壤中的质量分数，根据在水中传播的方式类似的定义了一个非线性沉淀模型。并推导出非线性沉淀情况下重金属速度的控制方程。

对于饱和土壤空隙内充填的溶液质量分数为ｗ，而土粒表面会吸附溶质，因此，

12