1. 引言
在我国,关于水资源承载力的研究经历了水资源研究理论的延伸与水资源承载力概念的形成、水资源承载力的开拓性与探索性研究这两个阶段后,进入到第三个阶段,即水资源承载力的专题性和系统性研究阶段 [1] 。到目前,水资源承载力评价融合,吸收了其他领域的精髓和思想,并应用于自身的发展。我国该领域常用来综合评价水资源承载力的方法有:
常规趋势法,考虑了单承载因子发展趋势,忽略了各承载因子间的相互关联。张琳等利用常规趋势法预测2010年江苏受水区的水资源承载力 [2] ;
系统动力学法,可以较好把握系统的各种反馈关系,但是该模型的建立受主观影响较大,变量较多,不易控制。王俭等应用系统动力学方法建立了水环境承载力模型,预测了2000~2050年辽宁省水环境在不同发展方案下的承载能力 [3] ;
多目标模型分析法,综合考虑各因素及其相互作用,但其决策中的各决策因子权重的确定也是通过主观性判断。杜立新等运用多目标模型分析法对不同方案下秦皇岛市水资源承载力的变化进行预测和解释 [4] ;
模糊综合评价法,克服了常规趋势法中承载因子间相互独立的局限性,但在其过程中选取部分影响较大的因子,导致部分有用信息遗漏。段新光等应用模糊综合评判模型对新疆水资源承载力现状进行评价 [5] ;
主成分分析法,通过降维处理技术将多个变量简化为少数几个综合指标,不仅保留原有指标体系的信息,而且还减少了信息的重叠。当然这样也存在部分有效信息的遗漏,最终的评价结果相对间接。李高伟等利用主成分分析法对郑州市水资源承载力进行了综合评价 [6] 。
本文应用投影寻踪聚类模型来综合评价水资源承载力,借助加速遗传算法寻找投影方向,从而最大限度避免了评判中权重取值的主观任意性。将指标体系投影得到一组投影值。可以直接根据投影值的大小来评价水资源承载力。
2. 水资源承载力衡量
水资源承载力是指在某一具体的发展阶段下,以可以预见的技术、经济和社会发展水平为依据,以可持续发展为原则,以维护生态环境良性发展为前提,在水资源合理配置和高效利用的条件下,区域社会经济发展的最大人口容量 [7] 。
建立投影寻踪聚类模型
本设研究样本集合为:
,n为样本个数,m为评价指标个数。
1) 样本指标数据标准化
数据标准化处理主要是为了解决两个方面的问题:其一是同趋化处理;其二是,数据无量纲化。将数据转化为趋于同向和无量纲的纯数值,便于不同性质或不同单位或量级的指标能够进行正确比较。为了解决这一问题,本文采用的是max-min标准化方法:
对于评价指标越大越好的评价值:
(1)
对于评价指标越小越好的评价值:
(2)
表示处理后的新数据,
表示第j个指标的最大值,
表示第j指标的最小值。
2) 线性投影
在各个决策变量的取值变化区间随机选取投影方向
若干组,每一组投影方向得到一组投影值,即为一维向量。
(3)
3) 寻找投影指标函数
根据投影值的散布图进行方案优选。最后得到的投影值集合
的散布特征要求:局部投影点越密集表示结果越好;但是在整体上投影点团之间越疏远越好。以此为限制条件,构造投影指标函数为:
(4)
其中:
(5)
(6)
(7)
表示为类间距离,
表示所有投影值的均值,
越大表示投影团半径越大。
表示类内密度,
越大表示分类的层次越清晰。
表示一阶单位阶跃函数,R为估计局部散点密度的窗宽参数,一般取
,距离
,![](//html.hanspub.org/file/5-2570114x30_hanspub.png)
4) 优化投影方向
当方案集选定时,投影指标函数
只与投影方向有关。由于选取的投影方向是随机的,误差比较大的,而且数量较多时想要找出最优的投影方向是不容易的。为了减小误差,本文借助加速遗传算法去寻找最优投影方向。在使用加速遗传算法过程中需要人为设定交叉概率
和变异概率
。然后再选取优秀数目若干个,设迭代3次进入加速循环,然后加速次数为20次。最终可以得出最优的投影方向
。该方向应满足:
(8)
(9)
式(8)为目标函数,式(9)为约束条件。
5) 最优投影值
用最优投影方向
与各指标的优度值相乘,得出最优的投影值。
(10)
表示
的第j个分量。最后得出的最佳投影值即是衡量水资源承载力的大小。投影值越大说明水资源承载力越大;相反地,投影值越小说明水资源承载力越小 [8] 。
3. 青岛市水资源承载力评价
3.1. 评价指标选取
评价指标体系是承载力研究的核心内容之一,其选择直接影响到最终评价的准确性。根据文献 [1] [7] [9] 选取以下七个因素作为评价青岛市各区域的水资源承载力的影响因素:
水资源开发利用率u1:区域供水量/水资源可利用量(%)。衡量地区对水资源的开发程度。
降水量u2:年平均降水量(mm)。对自然状态的描述。
人均占有水量u3:区域水资源总量/总人口(m3/人)。人均自然水资源的丰沛或紧缺的相对程度。
单位GDP用水量u4:E/GDP(m3/万元)。表示总用水量,GDP是对应的生产总值。是对经济发展状况的描述。
灌溉率u5:灌溉面积/土地面积(%)。农业灌溉用水量需求程度。
工业用水率u6:工业用水量/供水量(%)。工业用水量需求程度。
生态环境用水率u7:生态环境用水量/供水量(%)。表示生态环境或农村生态补水量 [9] 。
评价指标等级 [9] 如表1。
说明:1级属于情况很差,表示水资源承载力趋于饱和,会出现短缺,阻碍经济发展,该地区综合发展最好时,不可再过多增大人口数量、过度实施经济发展与开发;3级属于情况较好,表示该地区有很大的承载空间;2级介于两者之间,表示水资源已经开发到一定程度,若是合理采取方案,在一定程度上依然可以满足当地需求。
3.2. 模型求解
青岛市管辖有市区(包含市北区、市南区、李沧区)、崂山区、城阳区、黄岛区、墨市、莱西市、平度市和胶州市八个区域。以下数据都是按照该区域分类来整理的。
将表2中的数据进行标准化处理。用遗传算法寻找最佳投影方向时,取
,取
和
。选取优秀数目为20个,设迭代3次进入加速循环。采用MATLAB编程处理,加速次数为20次。得到最优投影方向
再用投影方向与各个指标优度值相乘得出最佳投影值 [8] 。
由表3可知,青岛市各地区的水资源承载力顺序为:3级 > 崂山区 > 市区 > 1级 > 即墨市 > 莱西市 > 平度市 > 全市平均 > 黄岛区 > 胶州市 > 城阳区。
4. 结论
1) 青岛市除了崂山区和市内三区的水资源承载力处在2级,其他各地区都处在1级。并且崂山区,
![](Images/Table_Tmp.jpg)
Table 1. System resulting data of standard experiment
表1. 标准试验系统结果数据
![](Images/Table_Tmp.jpg)
Table 2. The water resources carrying capacity evaluation index statistics of Qingdao
表2. 青岛市区域水资源承载力评价指标统计
![](Images/Table_Tmp.jpg)
Table 3. The standard values and index values of each region
表3. 各区域评价指标优度值
市区,即墨市,莱西市和平度市处在全市平均承载力之上,其余地区处在全市平均承载力之下。
2) 青岛市的大部分区域都处于水资源承载力以及青岛市的平均水资源承载力均趋于饱和状态.如果不采取相应的措施,水资源将会出现短缺,阻碍经济发展,不宜过多增加人口量。
3) 根据青岛市水资源总量和社会需求量显示,该市水资源的利用对社会的支撑能力较强,但其保障程度较差;各区域水资源承载力差异不显著,大都处于饱和状态,应结合区域资源特点,制定科学合理的水资源开发利用方案和社会经济发展战略。