APP下载

层次分析法在水资源承载力评估中的应用

2017-03-27王弛

东北水利水电 2017年3期
关键词:浑河特征向量分析法

王弛

(辽宁省汤河水库管理局,辽宁辽阳111000)

层次分析法在水资源承载力评估中的应用

王弛

(辽宁省汤河水库管理局,辽宁辽阳111000)

文中利用层次分析法从流域的水资源、社会、经济及生态环境的实际出发,评价该流域水资源承载能力的现状及未来发展水平,摸清流域的水资源承载状况,并且对规划年的水资源承载能力进行定量计算。

水资源;承载力;评估;层次分析法;浑河流域

1 层次分析法基本原理

层次分析法是把复杂问题中的各因素划分成相关联的有序层次,使之条理化的多目标、多准则的决策方法,是一种定量分析与定性分析结合的有效方法。用层次分析法做决策分析,首先要把问题层次化。根据问题的性质和预计达到的总目标,将问题分解为不同的组成因素,并按照因素间的相互影响以及隶属关系将因素按不同层次聚集组合,形成一个多层次的分析结构模型。最终把系统分析归结为目标层、准则层、指标层的相对重要性权值的确定或相对优劣次序的排序问题,从而为决策方案的选择提供依据。

2 确定多目标权重系数

2.1 建立递阶层次结构

这是层次分析法中最重要的一步。将问题分为若干组成部分或元素,按照属性的不同把这些元素分成若干组,每一组构成一个层次,层次之间互不相交。同一层次的元素作为准则,对下层次的全部或部分元素起支配作用,同时它又受上一层次元素的支配,因而形成了自上而下的逐层支配关系的递阶层次结构形式。

层次数与问题的复杂程度和所需要分析的详细程度有关。每一层次中的元素一般不超过9个,因同一层次中包含数目过多的元素会给两两比较判断带来困难,一个好的层次结构对于解决问题是极为重要的。层次结构是建立在决策者对所面临的问题具有全面深入的认识基础之上的,如果在层次的划分和确定层次之间的支配关系上举棋不定,最好重新分析问题,弄清问题各部分相互之间的关系。

2.2 构造比较判断矩阵

1)在建立递阶层次结构以后,上下层之间的隶属关系就被确定了。假定上一层的元素Ak作为准则,对下一层次的元素Cl,C2,…,Cn有支配关系,目的是在准则Ak之下按它们的相对重要性赋予Ci相应的权重Wi(i=1,2,…,n)。对于大多数社会经济问题,特别是对于人的判断起重要作用的问题,直接得到这些元素的权重并不容易,往往需要通过适当的方法来导出它们的权重。层次分析法采用T. L.Saaty提出的1~9标度法对不同评价指标进行两两比较,构造判断矩阵,此过程将思维数量化。

2)浑河流域水资源承载能力评价指标的AHP权重。根据图1所示的浑河流域水资源承载能力评价指标体系,容易得出浑河流域水资源承载能力自上而下的逐层支配关系的递阶层次结构,按照AHP模型,构建如图1。

构造判断矩阵,计算一致性指标(CI)和一致性比率(CR)。

图1 浑河流域水资源承载能力评价的层次结构图

3 层次单排序及一致性检验

求解判断矩阵的最大特征值max及其对应的特征向量W,即得到同一层次各指标相对于上一层指标的相对重要性的权重排序。在准则Ak下,n个元素Cl,C2,…,Cn排序权重计算,可以求解判断矩阵C的最大特征根max,并满足关系式:

其中:W为对应λmax特征向量,其分量Wi,i= 1,2,…,n为对应元素Cl,C2,…,Cn在准则Ak下单排序的权值。对于n阶判断矩阵,其最大特征根为单根,且λmax≥n。当λmax=n,其余特征根均为0时,则A有完全一致性。如果λmax稍大于n,而其余特征根接近于0时,则A具有满意的一致性。一般来说,由于客观事物的复杂性和人们认识上的多样性以及主观上的片面性和不稳定性,判断者只能给出估计判断或大概判断,即给出的值与客观存在的是有偏差的,这样就不能保证判断矩阵具有一致性。为了保证应用层次分析法得到的结论合理、正确,还需要对所构成的判断矩阵进行一致性检验。检验判断矩阵的一致性,需计算一致性指标:

式中:CR——判断矩阵的随机一致性比率;n——判断矩阵的阶数;CI——判断矩阵的一般一致性指标;RI——判断矩阵的平均随机一致性指标;RI的具体值见表1。

当随机一致性比率CR=CI/RI<0.10时,认为层次单排序结果具有满意的一致性。否则,就认为初步建立的判断矩阵是不能令人满意的,需要重新调整元素标度值,直到具有满意的一致性为止。采用MATLAB软件求出各个矩阵的最大特征值以及对应的特征向量。

表1 1-9阶矩阵RI值

对于Z-C比较判断矩阵一致性检验采用MATLAB软件求得的最大特征值λmax=4.198 0,对应的特征向量W=[0.118 6,0.398 9,0.208应用公式(2)求出

查表1得到RI=0.9;

其他计算出结果后依照一致性检验公式计算结果如下。

1)C1-A比较判断矩阵一致性检验

采用MATLAB软件求出最大特征值λmax= 8.933 2,对应的特征向量W=[-0.065 7,-0.087 7,-0.103 3,-0.385 8,-0.290 3,-0.210 7,-0.159 3,-0.821 2],

查表1得到RI=1.41;

2)C2-A比较判断矩阵一致性检验

采用MATLAB软件求出最大特征值λmax= 8.791 1,对应的特征向量W=[-0.6267,-0.506 7,-0.237 0,-0.106 9,-0.090 5,-0.290 1,-0.181 1,-0.397 2];

查表1得到RI=1.41;

3)C3-A比较判断矩阵一致性检验×

采用MATLAB软件求出最大特征值λmax= 8.935 6,对应的特征向量W=[-0.398 8,-0.131 5,-0.251 1,-0.213 2,-0.169 7,-0.084 3,-0.056 5,-0.822 2];

查表1得到RI=1.41;

4)C4-A比较判断矩阵一致性检验

采用MATLAB软件求出最大特征值λmax= 8.767 8,对应的特征向量W=[-0.122 2,-0.220 5,-0.600 0,-0.507 9,-0.320 0,-0.223 6,-0.402 5,-0.063 8];

查表1得到RI=1.41;

有上述结果可知,各个判断矩阵均满足一致性检验。见表2层次单排序表。

表2 层次单排序表

4 层次总排序

计算同一层次所有因素对整个总目标相对重要性的排序权值,称层次总排序。层次总排序是层次分析法的最终目的。为实现此目的,还需在层次单排序的基础上进行层次总排序。它是用下一层次各个因素的权值和上一层次因素的组合权值,得到最下层因素相对于整个总目标的相对重要性。

由表2可知各判断矩阵的特征向量,将特征向量归一化处理,得到二级指对于一级指标的权重,以及三级指标对于二级指标的权重,只需将三级指标的相对权重乘以二级指标的相对权重即可得到三级指标的最终的权重值。

首先将各矩阵的特征向量进行归一化处理,就是将对应的特征向量值除以各特征向量值之和,即计算某一指标的相对权重值的计算公式如下:

由公式(4)可以进行表3中的计算。

表3 层次总排序权值

从表3中可以看出,在二级指标中,水资源子系统和生态子系统的权重较大,说明水资源本身和生态环境对浑河流域水资源可持续发展有着重要的影响。P层次权重的结果显示需水模数和人口密度所占权重最大,说明了需水模数和人口密度这两个评价指标对浑河流域综合评价体系有着很大的影响力。

5 结语

利用层次分析法对浑河流域水资源评价指标体系中的8个指标进行权重的确定,水资源子系统和生态子系统的权重较大,说明水资源本身和生态环境对浑河流域水资源可持续发展有着重要的影响。需水模数和人口密度所占权重最大,说明了需水模数和人口密度这两个评价指标对浑河流域综合评价体系有着很大的影响力。

[1]许树伯.层次分析法原理[M].天津:天津出版社,1988.

[2]辽宁省水利厅.辽宁省水资源公报[R].2005.

[3]温淑瑶,马占青,周之豪,等.层次分析法在区域水资源可持续展评价中的应用[J].长江流域资源与环境,2000(2).

[4]刘哲,张江山.层次分析法在水库水质评价中的模型设计及应用[J].江苏环境科技,2005,18(2).

[5]郑哲敏.基于直觉模糊层次分析法的邯郸市水资源承载力评价研究[D].河北工程大学,2015.

[6]王海亮,杨永红,张颖.基于模糊层次分析法的水资源承载力研究[J].水科学与工程技术,2014(03).

[7]刘晓峰,许静.石羊河流域水环境承载力分析研究[J].水利规划与设计,2011(03).

[8]魏景忠,贾志国.水资源可持续利用初探[J].水利规划与设计,2009(05).

TV213.9

A

1002-0624(2016)12-0047-04

2016-08-16

猜你喜欢

浑河特征向量分析法
二年制职教本科线性代数课程的几何化教学设计——以特征值和特征向量为例
异步机传统分析法之困难及其克服
克罗内克积的特征向量
浑河各河段水质随季节变化趋势及影响因素研究
一类特殊矩阵特征向量的求法
流域尺度的浑河水污染防治环境经济政策
基于随机森林的浑河沈阳段水质评价
EXCEL表格计算判断矩阵近似特征向量在AHP法检验上的应用
基于时间重叠分析法的同车倒卡逃费探析
层次分析法在SWOT分析法中的应用