一、概述
随着滇藏铁路工程的分段实施，丽江-香格里拉段的规划设计已纳入日程。但是，由于该段地形地貌和地质条件非常复杂，虽然经过多轮论证，线路仍难最后确定。按照初期规划（图13-1），滇藏铁路丽江-香格里拉段共有3个走向方案可以比选：①丽江-长松坪-虎跳峡上峡口-香格里拉方案（西线方案）；②丽江-大具-白水台-小中甸-香格里拉方案（组合方案）；③丽江-大具-白水台-天生桥-香格里拉方案（东线方案）。初步分析认为，西线方案工程地质条件相对较好，可以作为推荐方案，该方案需要新建铁路隧道34座，总长87130 m，占该段线路总长的54.4%，最长的隧道是位于丽江西北的玉峰寺隧道，全长10970 m；需要新建铁路大桥39座（10253 m），涵洞182座（4547 m），桥涵占线路总长的9.2%。复杂的工程地质条件使得该方案仍存在许多问题，且工程建设难度大。
为了更好地指导该段铁路选线，我们在区域地壳稳定性评价的基础上，将基于GIS技术的层次分析法引入到丽江-香格里拉段铁路规划区的工程地质稳定性评价（工程地质条件评价）。在评价过程中，综合考虑地形坡度、工程地质岩组、斜坡结构、地质灾害发育现状、地壳稳定性、微地貌类型（地形与铁路设计高程高差）、人类工程活动、降水量、距离沟谷距离等因素，充分利用GIS技术处理海量数据信息的优势，采用层次分析法模型，进行丽江-香格里拉段铁路规划区的工程地质稳定性评价。基于评价结果，可以很好的指导该段线路比选和优化。
二、基于GIS的层次分析法原理
层次分析法（Analytical Hierarchy Process，简称AHP）是美国数学家SattyT.L.在20世纪70年代提出的一种将定性分析和定量分析相结合的系统分析方法。它适用于多准则、多目标的复杂问题的决策分析，可以将决策者对复杂系统的决策思维过程实行数量化，为选出最优决策提供依据（图13-2）。经过多年的应用实践，不少研究者开始将GIS技术与AHP方法相结合，大大提高了传统的AHP方法在地学研究中的应用效果（Harris et al.，2000；刘振军，2001；彭省临等，2005）。基于GIS的层次分析法充分利用GIS技术的空间分类和空间分析功能，在评价指标数据采集、处理和自动成图方面具有明显的优势，不仅可以对工程地质稳定性的相关影响因素进行更细致的逐次分析，而且在计算过程中不受计算单元数量的限制，因而评价结果更直观、更便于应用。
图13-1 滇藏铁路丽江-香格里拉段线路方案示意图
图13-2 基于GIS的层次分析法技术路线图
基于GIS层次分析法的工程地质稳定性分区评价过程大致可分为以下步骤：
（1）确定研究区、研究对象及研究目标，并进行数据分析，确定进行工程地质稳定性分区所需要的数据，包括数据来源、数据质量指标等。
（2）将收集的各种资料进行数据处理，包括在MapGIS 6.7软件平台上进行数字化、格式转换、投影转换、分层及属性编码等，建立研究区、研究对象的空间数据库。
（3）根据研究目标的特征，分析影响目标的因素，建立目标的层次指标模型和层次结构，构造判断矩阵，由专家对影响因素进行综合评分，并进行层次单排序、求解权向量和一致性检验，从而获得各指标因素值，并运用GIS空间分析功能提取分析因子。
（4）采用ArcGIS 9.2软件平台，对评价区域进行栅格化，每一个栅格作为模型评价的一个运算单元，并将数据库中的数据按照规则进行栅格化处理。再采用图形叠加的模型评价方式，将参与评价的各个因素权值分配到不同的栅格上。将各个因素进行图形叠加，对属性值进行代数运算，再将叠加后的栅格数据化，生成新的图形，并形成最终评价结果。
（5）工程地质稳定性分区评价的数学模型：
滇藏铁路沿线地壳稳定性及重大工程地质问题
式中：B——工程地质稳定性指数，aj——权重，Nj——指数。
（6）通过分析计算获得的工程地质稳定性指数值的分布范围，结合野外实际调查结果验证，对不同区域的铁路工程建设适宜性进行综合分区评价。