主成分分析法在清水河健康评估体系指标筛选中应用.docVIP

主成分分析法在清水河健康评估体系指标筛选中应用 　　中图分类号：TU973+.19 　　1.引言 　　河湖水系是地表水资源的主要载体，是维系生态系统健康的重要因子，也是哺育人类历史文明的摇篮。在大规模经济开发和全球气候变化双重作用下，一些河流、湖泊出现不同程度地恶化，形态、结构、水文条件变化，生境退化及重要或敏感水生生物消失等问题 。主成分分析也称主分量分析法，旨在利用降维的思想，把多指标转化为少数几个综合指标，是系统工程处理定性与定量相结合的简单易行且行之有效的一种系统分析方法。通过对现有指标的统计分析，提取出可以最大程度上代表和说明主体状况的指标，从而达到通过用较少指标得到全面结果的目的。 　　本文以云南省丘北县清水河为研究对象，运用主成分分析法对2011年-2012年进行的三次调查数据结果进行统计分析，结合实际的调查和评估工作，筛选出清水河健康评估体系指标。 　　2.研究区地理背景及数据源分析 　　2.1地理背景 　　清水河位于云南省文山州丘北县，东经103°50′～104°19′，北纬23°59′～24°13′。河长61.5km，落差727m，平均比降7.7‰，集水面积1533.5km2。《云南省水功能区划》将清水河从河源依次划为清水河丘北开发利用区和清水河丘北保留区，其中开发利用区兼有饮用、农业和景观用水功能，因此又划分为红旗水库饮用、农业用水区和普者黑景观、农业用水区两个二级区，规划水平年水质目标除红旗水库为Ⅱ类，其他两个区均为Ⅲ类。红旗水库为中型水库，总库容5370万m3，是丘北县的主要供水水源，现状水质为Ⅱ类。 　　2.2 数据来源及应用分析 　　研究中的数据主要来源于2011年-2012年的清水河健康评估调查，内容主要包括：水文水资源、物理结构、水质、水生生物、社会服务功能等方面。 　　水文水资源是以流量过程变异程度和生态流量满足程度作为判断河流是否健康的指标，河流流量的变化直接影响到周围社会经济活动和水生生物对水的需求。其评价的依据是最低生态水位，湖体运行水位是维持湖体各项功能正常运行的必要条件。 　　物理结构数据包括：湖滨带岸坡的基质、高度、倾角、植被覆盖、冲刷程度，以及河岸带区域人工干扰程度，数据获取有现场测量和卫星遥感影像数据。 　　水质数据包括由溶解氧代表的DO健康状况，由高锰酸盐指数、氨氮含量共同构成的耗氧有机物污染状况，考虑河流两岸可能会有加工工厂，在评价中添加了重金属指标。 　　生物多样性数据是一项综合指标，由于其具有复杂性和多尺度性，涉及多种因素影响，考虑生物多样性监测的适用性和可操作性，用简化的生物或环境特征参数来表征人类活动压力对生物多样性的影响，选择生物多样性指数、均匀度指数、硅藻生物指数、生物完整性指数等影响因子，作为生物多样性这一综合指标评估依据。 　　社会服务功能数据，指河流在当地的社会经济活动中的作用，包括判断水体是否能满足其生态与环境保护和资源利用（饮用水、工业用水、农业用水）的水功能区达标率；提示流域经济社会活动对水量影响，反应流域开发程度、社会经济发展与生态环境保护之间协调性的水资源开发利用率；以及防洪指标和公众满意度调查。 　　3. 清水河健康评估体系指标的筛选 　　3.1 方法说明 　　主成分分析数学模型 　　F2=a12ZX1+a22ZX2……+ap2ZXp 　　…… 　　Fp=a1mZX1+a2mZX2+……+apmZXp 　　其中a1i， a2i， ……，api（i=1，……，m）为X的协方差阵Σ的特征值多对应的特征向量，ZX1， ZX2， ……，ZXp是原始变量经过标准化处理的值，因为在实际应用中，往往存在指标的量纲不同，所以在计算之前须先消除量纲的影响，而将原始数据标准化。 　　A=（aij）p×m=（a1，a2，…am，），Rai=λiai，R为相关系数矩阵，λi、ai是相应的特征值和单位特征向量，λ1≥λ2≥…≥λp≥0。 　　进行主成分分析主要步骤如下： 　　① . 指标数据标准化（SPSS软件自动执行）； 　　② . 指标之间的相关性判定； 　　③ . 确定主成分个数m； 　　④ . 主成分Fi表达式，找到各指标的系数； 　　3.2 指标筛选 　　对指标的筛选过程所涉及的计算由SPSS（17.0）分析软件进行，其详细步骤为： 　　① 将调查数据导入SPSS（17.0）。 　　② 选择分析―降维―因子分析，将各数值变量进行因子选择，描述―原始分析结果―系数―显著性水平―继续―确定。得到所要分析数据的相关性矩阵、主成分分析的初始特征方差、各成份矩阵。 　　③ 将主成分矩阵的结果继续导入SPSS（17.0），进行各成分指标系数的计算，找出关键性指标。 　　用主成分分析法对清水河2011年10月，2012年4月、10月的三次水质监测数据中所用的11个指标进