基于主成分分析法南四湖流域水污染物排放水平空间布局分析.docVIP

基于主成分分析法南四湖流域水污染物排放水平空间布局分析 　　摘要：本文基于2007年第一次全国污染源普查数据，运用主成分分析法，对南四湖流域山东省境内的28个县市区行政单元的水污染物排放水平进行了综合评价，结果显示：滕州市的污染物排放水平明显高于流域内的其他区域，济宁市中区、邹城市、菏泽市区、曹县的污染物排放水平也较高。因此这五个县市是今后南四湖流域治污的重点。 　　关键词：主成分分析；南四湖流域；污染物排放；空间布局 　　中图分类号： X5 文献标识码：A DOI： 10.3969/j.issn1003-8256.2013.06.013 　　1 引言 　　南四湖是南水北调东线工程的重要调蓄水库[1]。国家要求流域内的水质达到国家Ⅲ类水质要求[2]。这对南四湖流域来说是个严峻的挑战，如何在尽量不影响经济发展的同时保证调水的水质，是急需解决的问题。本文根据2007年的污染普查数据，对当时南四湖流域的污染物空间分布进行分析，提出了相应的对策建议。 　　南四湖流域北起大汶河左岸，南至废黄河南堤，东以泰沂山脉的尼山为分水岭，西至黄河右堤，属于淮河流域泗河水系。总流域面积约3.2万平方公里，包括山东、江苏、河南、安徽四省38个县（市、区）。本文研究范围为南四湖流域山东境内所有区域，流域面积约2.9 万平方公里，包括济宁市、枣庄市、菏泽市和宁阳县三市一县[3-6]。 　　全国污染源普查是为贯彻落实科学发展观，加强环境监督管理，了解各类企事业单位与环境有关的基本信息，建立健全各类重点污染源档案和各级污染源信息数据库，为制定经济社会政策提供依据，国务院决定开展第一次全国污染源普查所得数据。全国污染源普查对象是我国境内排放污染物的工业污染源（简称“工业源”）、农业污染源（简称“农业源”）、生活污染源（简称“生活源”）和集中式污染治理设施。普查内容包括各类污染源的基本情况、主要污染物的产生和排放数量、污染治理情况等[7]。 　　2 分析方法介绍 　　2.1 主成分分析法 　　主成分分析法最早是由美国心理学家 Charles Spearman在1904年提出。他将实测的多个指标所表达的信息用少数几个潜在的独立的主成分指标来表达，而这些相???独立的主成分指标是由实测指标的线性组合来表示，来反映原来多个实测指标所要反映的主要信息[8]，即主成分分析法是把原来多个变量划为了少数几个综合指标的一种统计方法[9]。 　　2. 2 3σ法 　　3σ法通常用于测量误差的处理，它的基本思想是由于随机误差是服从于正态分布，误差的绝对值主要集中在均值 3σ附近，如有大于3σ的误差，是粗大误差，应予以剔除。从概率来看，有： 　　P（|d|3σ）=0.0027 （4） 　　所测的数据大于3σ的概率为2.7‰，这一概率是非常小的。因此如果所测数据大于3σ，可以认为其是极端的水平。基于3σ法原理，通过计算各个指标数据系列的中心值和标准差，本文确定以偏离中心值0.5倍标准差的区间即[ E-0.5σ，E+0.5σ ]属于中；偏离中心值1.5倍标准差和0.5倍标准差的区间即[E-1.5σ，E-0.5σ]属于中低水平和[E+0.5σ，E+1.5σ]属于中高水平区间；偏离中心值1.5倍标准差以上的区间即[-∞，E-1.5σ]属于低水平和[E+1.5σ，+∞ ]属于高水平区间。这样就得到了水污染物排放水平的5个区间。 　　2.3 数据处理 　　由于南四湖是南水北调东线工程重要的调蓄水库，其水质的好坏直接影响着调水的成败。本文以2007年全国污染物普查的数据，选取四个与水质相关的生活源污染物排放量，即COD排放量X1（吨）、氨氮排放量X2（吨）、总氮排放量X3（吨）、总磷排放量X4（吨），和四个与水质密切相关的工业源污染物排放量：化学需氧量排放量X5（吨），氨氮排放量X6（吨），石油类排放量X7（吨）和生化需氧量排放量X8（吨）等。运用SPSS18.0统计软件对污染物普查指标数据进行主成分分析，利用SPSS求得特征值和方差贡献率及累计贡献率，见表1。 　　表1 相关主成分累计贡献率 　　成份 初始特征值 　　特征根 方差的 % 累积 % 　　1 4.308 53.847 53.847 　　2 1.445 18.068 71.914 　　3 1.019 12.735 84.650 　　4 .893 11.164 95.814 　　5 .201 2.510 98.324 　　6 .128 1.604 99.928 　　从表中得出前4个主成分的累积百分比即贡献率已达到95.814%大于85%，故取前4个样本主成分。污染指标载荷值见表4，每个载荷量表示主成分与对应变量的相关系数。 　　表2 污染指标载荷值 　　成份 　　1 2 3 4