王 欣, 赵其华 , 张朝俊,王 红

(1.地质灾害防治与地质环境保护国家重点实验室(成都理工大学)，四川 成都 610059;2.成都理工大学 环境与土木工程学院，四川 成都 610059)

地质环境是生态环境的载体，生态环境是地质环境的“屏障”，生态环境的优劣与发展演化在很多方面受地质环境的影响，而两者的稳定运行又为人类生存环境提供了坚实的保障。生态地质环境作为一个组合概念，前人给出了不同的定义与研究方向[1-4]，但始终离不开岩石圈近地表受人为影响的地带。其本底值取决于地质环境，而承载力则受地质环境、生态环境、社会环境三者的共同制约与影响。目前，许多学者针对环境承载力进行了许多研究并取得丰硕的成果[5-10]，但其评价指标主要集中在生态环境系统与人类社会环境系统的组合上，有关地质环境系统方面的评价指标却很少。随着我国生态环境的日益恶化与山区城镇建设步伐的加快，许多学者围绕不同区域开展了生态地质环境承载力的综合评价，并取得丰硕的成果[11-14]，但由于评价指标众多且难以量化，导致评价结果适用性不强。因此，许多学者利用GIS强大的空间可视化及统计分析功能开展了不同区域的生态地质环境承载力研究[15-19]，其研究成果不仅直观地展现了承载力在空间上的差异性，而且为决策提供了有力的依据。但是，由于不同区域生态地质环境系统的复杂性与差异性，目前在评价过程中，权重的确定上仍存在一定的分歧，没有形成统一的认识。因此，本文从权重计算入手，通过将权重组合，降低权重计算的人为主观性，以期对金阳县生态地质环境承载力进行客观的评价。

1 金阳县生态地质环境背景

1.1 社会环境

金阳县地理坐标为102°56′～103°30′E，27°22′～27°58′N，位于四川省凉山彝族自治州东部边缘，坐落于金沙江大峡谷北岸，东南与云南省昭通市隔金沙江相望，西与布拖县毗邻，东北与雷波县相接，西北与昭觉县相连。金阳县域总面积约1 587 km2，共辖34个乡镇，截止2016年末总人口达20.6万人，人口密度129.8人/km2，国内生产总值28.9亿元，人均GDP达14 052元，对外交通主要依赖于省道S208线，交通安全通道单一，道路网络建设不均衡。

1.2 生态环境

金阳县全县年平均气温15.7℃，年平均降水量800 mm，属中亚热带季风季候，降雨集中在5月—10月，时空分布不均，局部暴雨现象时有发生。地下水资源贫乏，共发育大小23条河流，均属金沙江水系，三条主干流近南北向纵贯全县，从西向东依次为西溪河、金阳河、金沙江，其中金沙江境内全长98.8 km，沿河发育多条冲沟，受自然和人类活动的影响，水土流失严重。据遥感影像统计，金阳县耕地面积247.2 km2，占全县总面积的15.5%，且大都处于金沙江及其支流沿岸的陡坡地带。植被归一化指数介于-0.26～0.56，以草地、灌木林为主，植被覆盖率不均衡，全年空气质量良好，人均生态承载力0.85 hm2/人，略高于全国人均0.8 hm2/人[20]，随着人口的增长和社会经济发展，会出现生态赤字的情况。另外，金阳县耕地与草地生态承载力占比达78%，生态系统结构较差，一经破坏，将很难恢复。

1.3 地质环境

金阳县总体地势西北高东南低，是山原地貌向山地地貌过渡的区域，波洛梁子与东山梁子为境内两大山脉，以波洛梁子为界，西部为凉山山原地貌，东部为中山山地地貌，最高海拔4 036 m，最低海拔447 m，落差3 589 m。金阳县属上扬子台拗凉山隆褶区，为“川滇南北构造带”与“华夏式”构造交汇处，主要构造体系金沙江断裂束，属小江断裂波及范围，构造骨架为燕山运动期产物，以北东向和南北向构造为主[21]。断裂构造十分发育，其中金阳断裂近南北向穿越县城主城区，构造裂隙发育，断裂面两侧岩体挤压破碎强烈。区内地层除白垩系、石炭系缺失外，其余均有出露，寒武系、奥陶系、志留系地层广泛分布于县域中部，岩性以碳酸盐岩夹碎屑灰岩为主，二叠系峨眉山玄武岩在北部山原地区大面积出露，震旦系、二叠系、侏罗系地层以长条状分布于金沙江、金阳河沿岸，地层岩性主要为变质砂岩、灰岩、泥岩。据已有地震资料显示，金阳县未有7级以上地震记载，地震烈度Ⅶ度。截止2018年，全县共发育地质灾害(隐患)点211处，灾害点密度达14.36个/100km2，其中滑坡137处，崩塌23处，泥石流39处，不稳定斜坡12处，主要集中分布在断裂及水系附近(见图1)。

2 数据处理与评价模型

2.1 数据来源与评价流程

本研究采用的数据包括：(1)30 m分辨率的DEM(数字高程模型)数据，用于提取地形地貌数据、河网数据；(2)30 m分辨率的Landsat 8遥感影像数据(2018年)，用于提取植被归一化指数(NDVI)与土地利用情况数据，生态承载力数据在土地利用情况数据的基础上，利用生态足迹法计算得到；

图1 金阳县地质环境示意图

(3)1：5万数字地质图，用于工程岩组分类及地质构造矢量化；(4)地震动峰值加速度数据由国家地震科学数据共享中心获取；(5)四川蜀通岩土公司提供的2018年金阳县地质灾害编录数据及野外调查资料；(6)降雨量数据由中国气象数据网下载后空间插值处理得到；(7)空气质量数据由中国空气质量在线监测分析平台获取后空间插值处理得到；(8)人口、经济数据通过金阳县统计年鉴获取；(9)道路数据通过公开发表资料进行矢量化处理得到。

通过对金阳县生态地质环境系统的调查与相关资料的收集与分析，将金阳县生态地质环境系统划分为地质环境、生态环境和社会环境三个子系统。通过层次分析法构建了一套递阶式评价指标体系。地质环境子系统从地形地貌、地质灾害、地层岩性、地质构造四个方面考虑，共选取10个二级评价指标；生态环境子系统评价指标借鉴前人的研究成果[13-19]，主要从水、土资源与环境质量三方面进行考虑，共选取6个二级指标；社会环境子系统从人口、经济、基础设施建设状况三方面进行考虑，共选取3个二级指标。具体评价流程见图2。

2.2 评价指标体系分级

在原始数据收集的基础上，使用ArcGIS软件建立各评价指标的专题图层，并进行坐标系的统一转换与投影，具体处理流程在此不再赘述。在评价过程中，首先统计各指标数据的分布特征，依据自然间断法和金阳县实际情况对其进行4级划分，得出分级标准(见表1)。其次，将各指标图层转换成30 m×30 m空间分辨率的栅格数据，最后利用重分类功能，将相应栅格按分值4、3、2、1进行赋值。

图2 评价流程图

2.3 权重的确定

2.3.1 AHP法

层次分析法(AHP)可以合理地确定指标权重的排序，一般不会出现与实际重要程度相悖的情况[22]。每个指标的权重按递阶层次结构进行分系统、分层次计算得到。在层次模型建立的基础上，使用专家打分法对不同指标按重要性等级进行相互比较，从而构建出判断矩阵，采用“和积法”算出判断矩阵的特征向量，根据式(1)求得最大特征根，用式(2)、式(3)进行判断矩阵的一致性检验，若CR<0.1，则一致性检验通过，否则，修正原判断矩阵，直到一致性检验通过。经上述计算后得到主观权重(见表2)。

(1)

CI=(λmax-n)/(n-1)

(2)

CR=CI/RI

(3)

式中:λmax为最大特征根;(AW)i为向量AW的第i个分量;n为指标个数;CR为一致性比例;CI为一致性检验指标;RI为随机一致性指标。

表1 金阳县生态地质环境承载力评价指标体系与分级标准

注：MPa:表示岩石饱和单轴抗压强度；hm2表示公顷；表中的1、2、3、4为不同等级评分。

表2 基于层次分析法确定的主观权重

2.3.2 CRITIC法

CRITIC法能反映各指标之间差异分布的特性，通过指标间的对比强度和冲突性来计算权重[23]。其中，对比强度用标准差表示，标准差越大，差异性越明显，对承载力的影响越大。冲突性由相关系数计算得到，相关性越大，冲突性越小，指标越相似，对承载力的影响越相同。计算过程如下：首先计算各指标数据的标准差STD与协方差，确定协方差矩阵，由协方差矩阵计算得到相关系数rij矩阵。在相关系数计算的基础上，冲突性由式(4)获得，根据式(5)、式(6)分别计算得到各指标的信息量与客观权重(见表3)。

(4)

Cj=σj·Rj

(5)

(6)

式中：Rj表示冲突性;rij为指标i和j之间的相关系数;n为指标数量;Cj表示第j项指标所包含的信息量;σj为标准差;Wj表示第j项指标的客观权重。

表3 基于CRITIC法确定的客观权重

2.3.3 权重组合

本文采用线性功效系数法将主观权重与客观权重进行组合[24]，经线性加权后得到综合权重。计算过程如下：首先通过式(7)建立权重之间的函数关系，然后根据式(8)求得主、客观权重的分配系数，最后依据式(9)对主、客观权重进行线性加和，求得综合权重(见表4)。

(7)

(8)

wc=αwsi+βwoi

(9)

式中:ws、wo为主、客观权重;i为指标个数;α、β为主、客观权重分配系数;wc为综合权重。

表4 评价指标综合权重

以地质环境子系统B1各指标权重为例，将主观权重、客观权重、组合权重三种不同类型的权重值进行比较(见图3)发现，组合权重值介于主、客观权重值之间，一定程度上弥补了主、客观权重的不足。其他两个子系统类似。

2.4 评价模型

在建立指标体系分级标准和权重计算的基础上，构建金阳县生态地质环境承载力评价模型如下：

(10)

(11)

式中：C为生态地质环境承载力综合值；WBi为第i个子系统的权重值;Vi为第i个子系统的承载力评价值;Wij为第i个子系统的第j个指标权重;Cij为第i个子系统的第j个方案层指标归一化后的值;n为子系统的个数;m为第i个子系统中方案层评价指标的个数。

图3 地质环境子系统指标权重值比较

2.5 评价结果与分析

综合上述评价方法与模型，通过ArcGIS的栅格计算功能将各评价指标栅格图层乘以综合权重进行加权总和后得到各子系统的承载力空间分级图。利用自然断点分级方法，将承载力评价值分为优、良、中、差四个等级。对各系统不同等级所占面积进行统计后得到图4。

图4 各系统对应不同等级面积统计

结合图4各系统生态地质环境承载力分级图发现：

承载力为优的区域，主要分布在县城西北部，马依足乡与老寨子乡部分区域。该区地形较为平坦，平均海拔在2 200 m以上，地质灾害基本不发育，地层岩性以坚硬岩为主。发育小型断裂，地震动峰值加速度0.1g，植被茂密，人口分布较少。各子系统承载力均为优，地质环境系统、生态环境系统与综合评价结果面积相近，约占全县总面积的17.3%左右，社会环境系统为优的区域较多主要原因在于金阳县人口集中分布在少数几个乡镇，交通不发达，经济落后。

承载力为良的区域，主要分布在县城中西部，坡度8°～15°，地质灾害发育较少，地层岩性以较硬岩类为主，断裂以小型断裂为主，地震动峰值加速度0.1g，垂向上主要分布在海拔2 200 m以上的中高山区域以及1 300 m以上的河谷区域，植被覆盖中等，交通较为落后。各子系统与综合评价结果面积相近，约占全县面积的37%左右。

承载力为中的区域，主要分布在南东部及北东部，垂向上分布在海拔1 300 m～2 100 m之间的半山人口聚集区域，该区域地形起伏较大，往往位于缓坡与陡坎交界地带，人口密度中等，地质灾害中等易发，中部及南东沿江断裂较发育，地震动峰值加速度0.15g，地层岩性以软岩为主。地质环境系统、生态环境系统与综合评价结果面积相近，约占全县总面积的30.8%左右，社会环境系统为中的区域有100.4 km2，主要为县城中部人类工程活动加剧所致。

承载力为差的区域，主要分布在县城北部及中东部沿江区域，该区地质环境复杂，断裂构造密集，由于地形起伏较大，且山区河流流量随季节变化较大，河谷冲刷下切，侵蚀作用强烈，导致两岸岸坡陡峻，斜坡地质灾害极为发育，沿交通干线崩塌现象也比较普遍。地层岩性以泥岩、页岩、砂岩等互层状较硬岩为主。地震动峰值加速度为0.15g。植被破坏严重，为金阳县主要人口聚集地，耕地较多，北部区域生态环境脆弱，水土流失严重。地质环境系统，生态环境系统与综合评价结果十分相近，约占全县总面积的14.7%，社会环境系统急剧下降，只有30.1 km2，体现了人类社会经济活动相对集中的情况。

3 结论与讨论

本文以山区城镇金阳县为例，构建了金阳县生态地质环境承载力评价指标体系，将组合权重应用到评价模型中去，并进行了实证评价，得到了金阳县生态地质环境承载力的现状和空间格局，评价结果与实际调查结果较吻合。结果显示，金阳县整体生态承载力中等，生态环境优良，为长江上游生态环境建设提供了保障，但是复杂的地质环境又决定了其生态环境的脆弱性，一旦遭到破坏则很难恢复。另一方面，基于组合权重建立的评价模型，可以从更客观的角度对生态地质环境承载力进行定量综合评价。

与此同时，在研究过程中，由于数据收集的难易程度，没有做时间序列上的比较。对于降水与空气质量数据，县域尺度的空间分布只能通过地统计方法进行空间插值得到。此外，针对生态地质环境承载力的量化模型和评价方法应用相对较少，在进行综合评价时，可以结合不同评价目的对不同子系统进行单独的评价。