主成分分析法在水质评价中的应用.docxVIP

主成分分析法在水质评价中的应用 主成分分析法确定环境质量等级的依据 主成分分析作为一种常见的多指标统计方法，将原来的多个变量转化为多个索引指数。从数学的角度来看，这是一种减少维处理技术。主成分分析能够在最大限度地保留原始数据信息的基础上,对高维变量进行综合和简化,并且能够客观地确定各个指标的权重,避免了主观随意性,有一定 的优越性,已经被广泛地应用于环境领域。 在环境质量评价中,主成分分析法在确定环境质量等级时,其分级依据主要由以下两种方法得到:方法1,将环境标准单独进行主成分分析,以环境标准的主成分综合得分作为评价样本等级的判定依据;方法2,将全部或部分环境标准视为某年份的监测值加入评价样本中同时进行主成分分析,同样以标准的主成分得分作为分级依据。 文章先将主成分分析法用于该河流的水质评价,得出水质变化趋势,然后按方法1与方法2划分该河流水质等级,从得出的评价结果出发,对主成分分析法用于环境质量评价的合理性进行探讨。 1 主要成分分析用于评估实例 1.1 水质指标及溶解氧污染指数计算 某河流2001年-2007年的污染物浓度如表1所示,将各年份监测值与五个类别的水质标准值进行比较以确定水质级别。计算各年份的综合污染指标P?P=∑Pi=∑CiSi?CiΡ?Ρ=∑Ρi=∑CiSi?Ci为指标监测值,Si在此采用Ⅳ级标准,其中溶解氧污染指数的计算采用地面水环境影响评价技术导则(HJ/T2.3-93)推荐的溶解氧标准指数。 从表1可知影响各年份等级确定的主要因素是BOD和氨氮。 1.2 主要成分分析用于河流水质分析 1.2.1 主成分分析法计算水质变化趋势 根据主成分分析法原理,主要评价步骤有:(1)数据的标准化处理,用SPSS软件进行数据标准化之前,先将DO的监测值取倒数,使其与其它指标成为同向指标,即数值越大,污染越严重;(2)建立相关系数矩阵;(3)计算特征值、特征向量、主成分累计贡献率;(4)确定主成分。主成分个数提取原则为主成分对应的特征值大于1的前m个主成分;(5)计算综合得分。得分越大,水质越差。将主成分分析法的评价结果见表2。 根据表1中的污染指数值与表2中的主成分得分,将水质变化趋势见图1。 图1显示两条曲线变化趋势基本吻合,主成分分析法得出的水质变化趋势是2001年-2006年水质越来越差,2007年水质好转。这与综合污染指数P所反映的变化趋势基本相符。 1.2.2 主成分得分与关键因素对照 按照上述主成分分析法的步骤,对各级标准值组成的样本进行主成分分析,以各级标准的主成分得分作为分级依据,见表3。将表2主成分得分与表3进行对照确定评价等级,如表4所示2002年的主成分得分为-1.81164,其介于Ⅰ类标准得分-2.498和Ⅱ类标准得分-1.688之间,则确定其为Ⅱ类,其它依次类推。 表4显示,除2004年外其他年份的评价等级与实际水质完全不符,显然方法1的评价结果极不可靠。 1.2.3 各年份主成分分析法+标准 方法2是将监测值与各级水质标准放在一起进行主成分分析。两者的组合方法有多种,为了获取分级依据,本例将五级标准全部加入到样本监测值中,这种组合记为a0。另外,针对本例主要有三类、四类、五类和劣五类水,与之相关的标准就有二类、三类、四类和五类。于是另有四种组合:a1,监测值与二、三类标准;a2,监测值与三、四类标准;a3,监测值与四、五类标准;a4,监测值与第五类标准。评价等级的确定方法与方法1类似,同样是将各年份主成分得分与标准得分进行比较以确定等级。 见表5评价等级与水质等级完全不符。见表6显示除a1组合中2001年的评价等级与实际水质相符外,其它组合的评价结果与水质等级都不相符。显然,方法2的评价结果也不可靠。 综上,主成分分析法能准确分析水质随时间的变化趋势,若案例中的年份换成不同河流或不同地区,指标为水质指标或大气污染指标,则用主成分分析法可得出不同地区的环境质量优劣顺序。但无论是方法1还是方法2,由于不能将环境质量与评价标准相对应,都得不出符合实际的评价结果,因此要对其进行探讨。 2 监测数据的主成分分析 在主成分分析法中,数据集合的变异信息可以用全部变量方差的总和来测量。方差越大,数据中含的信息就越多,方差越小,所包含的信息量就越少。 主成分分析法的核心步骤采用数学方法求特征值,特征向量,确定主成分个数。我们知道有D(F1)D(F2)…D(Fm),Fm是前m个主成分。 所以,用数据变异大小反映数据中的信息,则第一主成分F1携带的信息量最大,F2次之,以此类推。根据确定主成分的原则:特征值大于1.0的或累计的方差贡献率85%以上的因子数定为主成分个数,可知,主成分分析法提取出的前m个主成分所携带的信息量总和为:∑h=1mD(Fh)=∑h=1mλh∑h=1mD(Fh)=∑h=1mλh。 为了便