基于主成分分析的郑州市重金属污染评价

基于主成分分析的郑州市重金属污染评价

0土壤地质环境随着城市经济的快速发展和城市人口的增加，人类活动对城市环境的影响日益突出。根据不同城市土壤和地质环境的评估，以及通过评估方法获得的大量数据，对城市环境进行评估。人类活动的影响下城市地质环境的发展模式日益成为人们的中心。按照功能划分,城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等,分别记为1类区、2类区……5类区,不同的区域环境受人类活动影响的程度不同.本文对郑州市城区土壤地质环境进行调查.将所考察的城区划分为间距1km左右的网格子区域,按照每平方千米1个采样点对表层土(0~10cm深度)进行取样、编号,并用GPS记录采样点的位置.应用仪器测试分析获得了每个样本所含的多种化学元素的浓度数据.另外,按照2km的间距在远离人群及工业活动的自然区取样,将其作为该城区表层土壤中元素的背景值.1测量精度和污染程度.假设每类区表层土壤中每种重金属含量连续变化;在对采样点分析时,忽略土壤pH值、温度、气候等因素对测量带来的影响;在对每类区域的污染程度进行分析时,忽略考察的8种元素以外的重金属污染;忽略重金属由于土壤特性和氧化还原等一系列的物理和化学作用带来的对污染程度评价的误差.由于不易判断5类区的地理分布情况,我们将问题理想化,考虑每类区的8种重金属污染程度;假设附件给出的数据对于建立模型已经足够充分;假设每个采样点的采集时间相同,重金属的含量不会因为采集时间的差异出现变化;忽略污染化合物之间的相关影响;忽略重金属污染物之间的交叉污染;将土壤理想化成成分均匀的介质.1.1不同地区重金属污染本文使用单因子指数法和综合因子指数法分析该城区内不同区域重金属的污染程度.1.1.1综合因子指数法首先采用单因子污染指数法对土壤重金属污染进行评价,计算公式为Ρi=1mm∑i=1Ci/Si.Pi=1m∑i=1mCi/Si.然后采用综合因子指数法对土壤重金属污染进行评价,计算公式为P综=√(1nn∑i=1Ρi)2+(maxΡi)22=(1n∑i=1nPi)2+(maxPi)22−−−−−−−−−−−−−−−ue001⎷ue000ue000,式中,Pi为土壤污染物i的环境质量指数,Ci为土壤污染物i的实际测量值,Si为土壤污染物i的背景值,m是土壤污染物i采样点的个数,n为重金属种类总数.1.1.2评估标准为标准评价土壤环境的污染程度,依据国家环保局颁布的土壤环境质量标准,土壤中各元素环境质量评价分级含量范围采用表1分级值.1.1.3重金属污染程度使用EXCEL进行数据分析与处理,得到如下评价结果:生活区中Cu和Zn属于重度污染,工业区中Hg和Cu的污染程度极其严重,山区各重金属污染程度均较低,交通区中也是Hg的污染相当严重,公园绿地区各重金属污染程度相对较好.1.2评价指标权重的确定显然,用此方法计算的结果会突出最严重污染物产生的影响,但没有考虑污染因子的危害性差异.会引入很大的误差,因此需对此式进行修改.将污染指数平均值改为加权平均,根据评价指标对环境及人体的危害性来确定各评价指标的权重.首先将综合评价指标按由小到大的顺序排列为P综1,P综2,…,P综n.将最大的一个记为Pmax,令γi=Pmax/Pi表示第i种评价指标的相对重要性比值,则ωi=γi/5∑i=1γiωi=γi/∑i=15γi为各评价指标的权重.修正后得到不同城区的综合污染程度:工业区污染指数为8.374,属于重度污染;其次是交通区,污染指数为5.784;生活区和公园绿地区属于中度污染,污染指数分别为2.569和2.227;山区属于轻度污染,污染指数为1.214.2重金属污染的主要原因2.1初步分析通过对重金属污染的数据分别进行横向和纵向分析,找出重金属污染的主要原因.2.1.1工业污染的污染程度.城市水系城市的污染程度比人民生活城市更横向比较5类综合污染指数,可以清晰地看出,工业区所占比例明显高出其他四类区域,交通区次之,生活区和公园绿化区持平居后,山区的污染程度最轻.这说明工业污染导致重金属污染的成分最重,交通区主要以排放的污染气体为主,被污染的大气长期滞留,这样就会污染到土壤表层.对于生活区和公园绿化区,人们的日常生活基本一致,而公园绿化区的污染程度要比生活区稍轻一些,这跟实际情况是相吻合的.横向分析可以看出,工业区的综合指数最高,这与通常情况相符,说明工业区周围的治理污染工作还不到位,横向主要原因为工业污染.2.1.2hg污染原因分析纵向比较污染程度最严重的区域(工业区),分析8种重金属的单因子指数,这样可以明显找出导致污染的主要因素.可知Hg的污染指数最高,工业区土壤中Hg明显富集,这与工业区的Hg元素及含Hg的化合物渗透土壤的污染物比较多有关,下面采用主成分分析法确定污染物产生的主要原因.2.2---先利用上述处理结果,建立5个区域和8种元素的5行8列的单因子矩阵,X=(1.742.232.233.742.661.492.233.432.013.021.729.6618.351.613.004.031.121.171.261.311.171.261.181.061.592.771.874.7112.771.432.053.521.742.161.412.293.291.241.962.24).X=⎛⎝⎜⎜⎜⎜⎜⎜1.742.011.121.591.742.233.021.172.772.162.231.721.261.871.413.749.661.314.712.292.6618.351.1712.773.291.491.611.261.431.242.233.001.182.051.963.434.031.063.522.24⎞⎠⎟⎟⎟⎟⎟⎟.然后利用MATLAB计算X的相关系数方阵R=(rij),其中rij为R中i行j列的元素,其计算公式为rij=n∑k=1(xki-ˉxi)(xkj-ˉxj)√n∑k=1(xki-ˉxi)2n∑k=1(xkj-ˉxj)2.rij=∑k=1n(xki−x¯i)(xkj−x¯j)∑k=1n(xki−x¯i)2∑k=1n(xkj−x¯j)2√.由此公式可以知道rij=rji,然后计算上述方阵的特征值λ,令|λΙ-R|=0|λI−R|=0,由MATLAB计算出结果然后进行排序λ1≥λ2≥…≥λp≥0,得到特征向量ei‚p∑j=1eij2=1ei‚∑j=1peij2=1,得到贡献率λi/p∑k=1λkλi/∑k=1pλk,累计贡献率λ′i=i∑k=1λk/p∑k=1λkλ′i=∑k=1iλk/∑k=1pλk,利用MATLAB计算结果见表2.取累计贡献率达85%~95%的特征值λ1,λ2,…,λm所对应的第一、第二……第m(m≤p)m(m≤p)个主成分,则得到3个主成分,分别为第一主成分F1、第二主成分F2和第三主成分F3,F1=0.3851x1+0.2011x2+0.00538x3+0.25894x4-0.7180x5-0.1248x6-0.4541x7-0.0829x8,F2=0.3839x1-0.0680x2-0.2628x3-0.3274x4-0.2553x5-0.5994x6-0.0429x7-0.4955x8,F3=0.3635x1+0.1182x2+0.4779x3-0.4283x4-0.1244x5-0.2610x6+0.2952x7+0.5209x8.由此确定,工业区的Hg、Cu,交通区的As、Hg,生活区的Cu分别为重金属污染的主要元素.3确定震源位置3.1污染位置的确定对于重金属污染物的空间分布,结合数据的变化情况,首先从几个侧面分析污染物的传播特征,然后根据最速下降法进行方向追踪污染源的位置.整合数据分析可得,平面角度:以污染物浓度较高的地理位置为圆心扩散,区域划分角度:工业区、交通区、生活区、公园绿化区、山区.然而值得注意的是,从元素的种类上看,某些重金属不会扩散,也就是说,其在土壤的稳定性很强,在一定时期内不会传播.所以,在确定传染源的时候,需要进行模糊聚类,将一些传播明显的元素进行梯度分析,进而得到最速下降最快的方向,也就得到了污染源的位置.3.2土壤重金属污染的成因土壤中重金属元素主要有自然来源和人为干扰输入两种途径.在自然因素中,成土母质和成土过程对土壤重金属含量的影响很大.在各种人为因素中,则主要包括工业和交通等来源引起的土壤重金属污染.以下主要对受人为作用影响的土壤重金属污染来源进行介绍.3.2.1企业固体废弃物排放废水工业过程中广泛使用重金属元素,工矿企业将未经严格处理的废水直接排放,使得周围的土壤容易富集高含量的有毒重金属.企业排放的烟尘、废气中也含有重金属,并最终通过自然沉降和雨淋沉降进入土壤.矿业和工业固体废弃物在堆放或处理过程中,由于日晒、雨淋、水洗等原因,重金属极易移动,以辐射状、漏斗状向周围土壤扩散,固体废弃物也可以通过风的传播而使污染范围扩大.3.2.2重金属道路一般会在道路,一般会在道路,一般分为5105和5.道路两侧土壤中的污染物主要来自汽车尾气排放及汽车轮胎磨损产生的大量含重金属的有害气体和粉尘的沉降,而污染元素则主要为Pb、Cu、Zn等元素.它们一般以道路为中心呈条带状分布,强度因距公路、铁路、城市的距离以及交通量的大小有明显的差异.3.3确定震源位置3.3.1污染源的估计为简化模型,从Cu的污染分布着手讨论污染源位置确定的模型,其他元素依次代入模型即可.污染源(x0,y0,z0)=(0,0,0)使用INGLADA方法进行初值计算,设观测点坐标为(xi,yi,zi),初步估计污染源多于2个,污染源到观测点的距离为R,对于污染源i有(x0-xi)2+(y0-yi)2+(z0-zi)2=R2i(x0−xi)2+(y0−yi)2+(z0−zi)2=R2i.(1)对于污染源j有(x0-xj)2+(y0-yj)2+(z0-zj)2=R2j(x0−xj)2+(y0−yj)2+(z0−zj)2=R2j.(2)将(1)式、(2)式展开后相减得到线性方程(xi-xj)x0+(yi-yj)y0+(zi-zj)z0=12(r2i-r2j+R2j-R2i)(xi−xj)x0+(yi−yj)y0+(zi−zj)z0=12(r2i−r2j+R2j−R2i),(3)其中r2i2i=x2i2i+y2i2i+z2i2i.由于土壤为均匀介质,传播速度恒定,得(xi-x0)2+(yi-y0)2+(zi-z0)2(Ci-C0)=(xj-x0)2+(yj-y0)2+(zj-z0)2(Cj-C0)(xi−x0)2+(yi−y0)2+(zi−z0)2(Ci−C0)=(xj−x0)2+(yj−y0)2+(zj−z0)2(Cj−C0).(4)将观测点坐标的数据每3个一组代入(3)式,得到一系列(x0,y0,z0)(x0,y0,z0),从而反演得源位置.数据每4个一组共分成80组,代入(4)式中,得到80组(x0,y0,z0)(x0,y0,z0)污染源位置.统计得x0的波动范围为0~2454,y0的波动范围为987~1570,从而确定平面污染源的区域范围为{0≤x0≤2454987≤y0≤1570{0≤x0≤2454987≤y0≤1570.由于采样点是在间距1km左右的网格区域内,|zi-zj||zi−zj|远小于|xi-xj||xi−xj|和|yi-yj|‚(3)|yi−yj|‚(3)式中的系数矩阵会有奇异性,在z0解发散,充分说明模型适宜此问题.3.3.2其他接触问题的补偿基于GEIGER法的基本思想,建立空间直角坐标系ΔC+∂Ck∂xΔx+∂Ck∂yΔy+∂Ck∂zΔz=rkΔC+∂Ck∂xΔx+∂Ck∂yΔy+∂Ck∂zΔz=rk,(5)式中Δx,Δy,Δz分别为污染源坐标校正量,ΔC为浓度校正量,∂Ck∂x,∂Ck∂y,∂Ck∂z∂Ck∂x,∂Ck∂y,∂Ck∂z为空间偏导数,rk为测定浓度的偏差,可写成矢量形式AΔ→x=→r‚其中Δ→x=(ΔC,Δx,Δy,Δz)‚→r=(r1,r2⋯‚rm).(6)(6)式中系数矩阵A为m×4矩阵A=(1∂C1∂x∂C1∂y∂C1∂z1∂C2∂x∂C2∂y∂C2∂z⋮⋮⋮⋮1∂Cm∂x∂Cm∂y∂Cm∂z).(7)将INGLADA方法计算的初值代入(5)式求解,再用求解结果校正r,如此反