基于多态贝叶斯网络的喜马拉雅山冰碛湖溃决灾害评价与应用研究

基于多态贝叶斯网络的喜马拉雅山冰碛湖溃决灾害评价与应用研究一、引言1.1研究背景与意义1.1.1喜马拉雅山冰碛湖溃决灾害现状喜马拉雅山脉作为世界屋脊的重要组成部分，是中低纬度地区最大的现代冰川分布区，发育着世界上罕见的海洋性山谷冰川。在漫长的地质历史时期与复杂的气候环境作用下，这里形成了数量众多的冰碛湖。据相关研究统计，中国喜马拉雅山地区共有1680个冰碛湖，总面积达215.28平方千米。这些冰碛湖广泛分布在不同的海拔高度与地形条件下，从高耸的山谷到相对平缓的冰川槽谷，都能发现它们的踪迹。它们的存在不仅是该地区独特自然景观的一部分，更是区域水文循环系统的重要组成部分，对维持当地生态平衡、水资源补给等方面发挥着重要作用。近年来，随着全球气候变暖趋势的加剧，喜马拉雅山地区的气温呈现出明显的上升态势。相关数据显示，从1961年至2020年，该地区年平均气温每10年上升0.35℃，超过全球同期增温速率的2倍。气温的持续升高导致冰川加速消融，为冰碛湖提供了更多的水源补给，使得冰碛湖的面积不断扩张，数量也有所增加。同时，冰川的消融还导致冰体结构不稳定，容易引发冰崩、雪崩等灾害，这些灾害一旦发生，大量的冰雪物质涌入冰碛湖，会造成湖水位急剧上升，增加了冰碛湖溃决的风险。在过去的几十年里，喜马拉雅山地区发生了多起冰碛湖溃决事件，给当地的生态环境、居民生命财产安全以及基础设施建设带来了巨大的破坏。1981年，喜马拉雅中部波曲河流域的次仁玛错冰碛湖发生溃决，强大的洪水瞬间冲毁了超过55km的中尼公路路段，导致交通中断，数百人在这场灾难中失去生命。又如，当地时间2023年10月4日午夜，位于喜马拉雅山南麓的印度东北部锡金邦北部，受暴雨影响，冰川湖洛纳克湖发生溃决，大量冰冷的高山湖水倾泻而下，涌向提斯塔河谷，形成山洪。此次洪水导致超过50人遇难，150人失踪，受洪水影响的人口达到22000人，村庄乡镇被淹没，锡金邦多条重要的公路交通线中断，下游的提斯塔三号水坝被冲毁，15座桥梁也悉数被冲毁。这些溃决事件不仅造成了人员伤亡和财产损失，还对当地的生态环境产生了长期的负面影响，如破坏了河流生态系统、导致水土流失加剧等。此外，研究表明，1980年以来，冰碛湖溃决洪水的发生频率呈较弱的增长趋势，特别是2010-2018年间，喜马拉雅山脉中段发生8起冰湖溃决洪水事件，远高于天山山脉、喜马拉雅山脉东段和念青唐古拉山等地区，成为新的高发区。冰碛湖溃决的诱因复杂多样，其中冰-雪崩或冰川滑塌是导致冰碛湖溃决的主要原因，占比达到50.1%；埋藏冰融化或管涌、强降水或上游来水、滑坡-岩崩以及地震等因素也不容忽视，其占比分别为23.1%、18.5%、7.4%和0.9%。1.1.2研究意义喜马拉雅山冰碛湖溃决灾害的研究具有极其重要的现实意义，对灾害防控和保障居民安全至关重要。通过深入研究冰碛湖溃决灾害，能够准确识别潜在危险性冰碛湖，提前对其进行监测和预警，为当地居民争取更多的疏散时间，从而有效减少人员伤亡。通过对溃决风险的评估，可以制定针对性的防控措施，如修建防护堤坝、疏导引流等，降低灾害发生时的损失。在基础设施建设方面，对于穿越冰碛湖分布区的铁路、公路等交通线路以及水电设施等的规划和建设，冰碛湖溃决灾害的研究成果能够提供科学依据，避免在高风险区域进行建设，或者采取相应的工程措施提高基础设施的抗灾能力，保障其安全稳定运行。从理论意义来看，该研究对冰川学和灾害学的发展具有重要推动作用。在冰川学领域，研究冰碛湖的形成、演化与冰川消融之间的关系，有助于深入理解冰川的动态变化过程，揭示冰川物质平衡、能量交换等方面的内在机制，丰富和完善冰川学的理论体系。在灾害学方面，冰碛湖溃决灾害作为一种特殊的地质灾害类型，对其溃决机制、风险评估方法等的研究，可以拓展灾害学的研究范畴，为其他类似灾害的研究提供借鉴和参考，推动灾害学在多灾种、多尺度研究方面的发展。此外，冰碛湖溃决灾害的研究还涉及到气象学、水文学、地质学等多个学科领域，通过跨学科的研究方法，能够促进不同学科之间的交叉融合，为解决复杂的地球科学问题提供新的思路和方法。1.2国内外研究现状冰碛湖溃决灾害的研究在国内外都受到了广泛关注，学者们从不同角度进行了深入探索，取得了一系列有价值的成果。在冰碛湖溃决灾害的形成机制方面，国内外研究较为深入。国外学者如Huggel等通过对阿尔卑斯山地区冰碛湖的研究，指出冰碛湖溃决主要由冰崩、雪崩、滑坡等因素触发，这些因素会导致冰碛坝的稳定性降低，进而引发溃决。例如，在阿尔卑斯山的某些地区，由于气温升高，冰川融化加速，冰体的稳定性变差，容易发生冰崩，大量冰体落入冰碛湖，引发涌浪，当涌浪高度超过冰碛坝时，就会导致冰碛湖溃决。在国内，徐道明等对我国西藏地区溃决冰碛湖进行考察后，提出坝体特征（如坝体高度、宽度、物质组成等）和气候因素（如气温、降水的变化）是判别危险冰碛湖的重要指标。在西藏的一些地区，冰碛坝内存在埋藏冰，随着气温升高，埋藏冰融化，坝体结构被破坏，增加了冰碛湖溃决的风险。刘时银等研究发现，冰碛湖溃决与冰川的退缩、冰碛坝的稳定性以及冰湖的水文条件密切相关。当冰川退缩速度加快，冰碛湖的补给水源发生变化，同时冰碛坝在长期的水流侵蚀和冻融作用下，稳定性下降，容易引发溃决。对于冰碛湖溃决灾害的评价方法，国内外也有诸多研究成果。国外常用的评价方法包括基于地理信息系统（GIS）和遥感（RS）技术的冰碛湖特征提取与分析，以及利用数值模型对溃决洪水进行模拟和风险评估。如Richardson和Reynolds利用遥感影像对喜马拉雅地区的冰碛湖进行监测，分析其面积、水位等变化情况，评估溃决风险。他们通过对不同时期遥感影像的对比，能够及时发现冰碛湖的扩张或收缩，以及冰碛坝的形态变化，从而判断溃决风险的高低。国内学者则提出了多种评价指标和模型。吕儒仁等将冰碛湖与母冰川变坡点（或裂隙面）下游冰体的体积比定义为危险性指数，用于评估由冰崩涌浪造成坝体溃决的风险；陈晓清等从喜马拉雅地区冰碛湖相关因子中选取7个参数作为直接判别的依据，结合冰碛湖溃决危险指数建立起一种定性-半定量的个体冰碛湖危险性评价方法；舒有锋等通过粗糙集理论确定各指标权重，建立了冰碛湖溃决危险性评价方法。王欣等将冰碛湖溃决机制的定性描述与溃决概率进行转换，建立了冰碛湖溃决的事故树模型，对我国青藏高原以南冰碛湖溃决危险性进行分析，通过对溃决机制的详细分析，确定各个因素对溃决概率的影响，从而更准确地评估冰碛湖的溃决危险性。在冰碛湖溃决灾害的应用研究方面，国外在灾害预警和防治措施方面取得了一定进展。在一些冰碛湖分布集中的地区，建立了完善的监测体系，通过实时监测冰碛湖的水位、水温、冰碛坝的变形等参数，及时发出溃决预警。同时，采取了一系列防治措施，如加固冰碛坝、修建排水通道等，以降低溃决风险。国内则针对具体区域的冰碛湖溃决灾害进行了深入研究和应用实践。以西藏定结县龙巴萨巴湖为例，科研人员对其开展了溃决风险评价、溃决洪水模拟和减灾减缓的应用研究，通过建立数学模型，模拟溃决洪水的演进过程，为制定合理的减灾措施提供科学依据。还结合当地的实际情况，提出了如生态修复、移民搬迁等减灾减缓措施，以减少冰碛湖溃决灾害对当地居民和生态环境的影响。尽管国内外在冰碛湖溃决灾害研究方面取得了显著成果，但仍存在一些不足。在形成机制研究方面，虽然对主要的触发因素有了一定认识，但对于多种因素相互作用下的溃决过程，以及冰碛坝内部结构和力学特性对溃决的影响，还缺乏深入研究。在评价方法上，现有的评价指标和模型大多基于有限的数据和特定的区域，普适性和准确性有待进一步提高，不同评价方法之间的对比和验证也相对较少。在应用研究中，灾害预警的时效性和准确性仍需提升，防治措施的有效性和可持续性也需要进一步评估和改进，且跨区域的冰碛湖溃决灾害协同研究相对薄弱。1.3研究内容与方法1.3.1研究内容冰碛湖溃决灾害评价方法构建：深入分析冰碛湖溃决的形成机制，全面梳理导致冰碛湖溃决的各类因素，包括冰碛坝的稳定性、冰川的消融状况、气象条件以及地质构造等。基于此，选取具有代表性的评价指标，如冰碛坝的高度、宽度、物质组成，冰碛湖的面积、水深、蓄水量，以及周边冰川的厚度、流速等，构建科学合理的冰碛湖溃决灾害评价指标体系。同时，对现有的冰碛湖溃决灾害评价模型进行深入研究和对比分析，结合喜马拉雅山地区的实际情况，对模型进行改进和优化，提高模型的准确性和适用性。例如，考虑到喜马拉雅山地区复杂的地形地貌和多变的气候条件，在模型中增加地形坡度、降水强度等变量，以更准确地模拟冰碛湖溃决的过程和风险。冰碛湖溃决灾害评价模型验证：收集喜马拉雅山地区历史冰碛湖溃决事件的数据，包括溃决发生的时间、地点、原因、溃决规模以及造成的损失等信息。运用构建的评价模型对这些历史事件进行模拟和分析，将模型预测结果与实际情况进行对比验证，评估模型的可靠性和精度。通过对模型的验证，发现模型在某些情况下存在一定的误差，进一步分析误差产生的原因，如数据的准确性、模型参数的选取等，并对模型进行调整和改进，以提高模型的预测能力。同时，利用实地观测数据对模型进行验证，在喜马拉雅山地区选取多个具有代表性的冰碛湖，设置监测站点，实时监测冰碛湖的水位、水温、冰碛坝的变形等参数，将监测数据代入模型进行验证，确保模型能够准确反映冰碛湖的实际情况。喜马拉雅山冰碛湖溃决灾害案例分析：选取喜马拉雅山地区具有代表性的冰碛湖，如西藏定结县的龙巴萨巴湖等，对其进行详细的溃决风险评价。综合考虑冰碛湖的自然条件、周边环境以及人类活动等因素，运用构建的评价方法和模型，对冰碛湖的溃决风险进行量化评估，确定冰碛湖的危险等级。以龙巴萨巴湖为例，该湖位于喜马拉雅山北麓，周边冰川众多，气候复杂多变。通过对龙巴萨巴湖的地形、地质、气象等数据的分析，结合冰碛湖溃决灾害评价模型，评估出该湖的溃决风险较高，属于高危险冰碛湖。针对评估结果，制定相应的防灾减灾措施，如加强监测预警、修建防护堤坝、制定应急预案等，以降低冰碛湖溃决灾害的损失。同时，对防灾减灾措施的实施效果进行评估，根据评估结果对措施进行调整和完善，确保措施的有效性和可行性。1.3.2研究方法文献研究法：全面收集国内外关于冰碛湖溃决灾害的研究文献，包括学术论文、研究报告、专著等，了解冰碛湖溃决灾害的研究现状、形成机制、评价方法以及防治措施等方面的研究成果。对这些文献进行系统的梳理和分析，总结前人研究的优点和不足，为本文的研究提供理论基础和研究思路。通过对文献的研究，发现目前冰碛湖溃决灾害的研究在评价方法的普适性和准确性方面还存在一定的问题，需要进一步改进和完善。案例分析法：深入研究喜马拉雅山地区及其他地区的冰碛湖溃决灾害案例，详细分析这些案例中冰碛湖的形成过程、溃决原因、溃决模式以及造成的损失等。从案例中总结冰碛湖溃决灾害的发生规律和特点，为构建评价方法和模型提供实际案例支持。以1981年喜马拉雅中部波曲河流域的次仁玛错冰碛湖溃决事件为例，通过对该案例的分析，发现冰崩涌浪漫顶是导致此次溃决的主要原因，同时冰碛坝的稳定性较差也是一个重要因素。这些发现为冰碛湖溃决灾害评价指标的选取提供了重要参考。模型构建法：依据冰碛湖溃决灾害的形成机制和影响因素，运用数学、物理学等相关知识，构建冰碛湖溃决灾害评价模型。在模型构建过程中，充分考虑各种因素之间的相互关系和作用，确保模型能够准确地模拟冰碛湖溃决的过程和风险。例如，运用流体力学原理构建冰碛湖溃决洪水模拟模型，考虑冰碛湖的水位、流量、流速等因素，模拟溃决洪水的演进过程，评估洪水对下游地区的影响。实地调研法：深入喜马拉雅山地区，对冰碛湖进行实地考察和调研。通过实地观测，获取冰碛湖的地理位置、形态特征、周边环境等第一手资料，了解冰碛湖的实际情况。同时，与当地居民和相关部门进行交流，了解冰碛湖溃决灾害的历史情况和当地的防灾减灾措施。在实地调研过程中，发现一些冰碛湖周边存在人类活动，如放牧、修路等，这些活动可能会对冰碛湖的稳定性产生影响，需要在评价过程中加以考虑。此外，还可以在实地设置监测站点，安装相关监测设备，对冰碛湖的水位、水温、冰碛坝的变形等参数进行长期监测，为模型验证和灾害预警提供数据支持。二、喜马拉雅山冰碛湖溃决灾害形成机制2.1冰碛湖的形成喜马拉雅山冰碛湖的形成与冰川的演化密切相关，是一个长期且复杂的地质过程，主要源于冰川消退以及冰碛物的堆积阻塞。在漫长的地质历史时期，喜马拉雅地区经历了多次冰期与间冰期的交替。在冰期，气候寒冷，大量积雪在重力作用下逐渐形成冰川。冰川在自身重力和地形的影响下，沿着山谷缓慢向下滑动。在滑动过程中，冰川犹如一台巨大的刨蚀机器，对其经过的山体和谷底进行强烈的侵蚀，将大量的岩石碎屑、泥沙等物质裹挟其中。这些被冰川携带的物质统称为冰碛物，它们大小混杂，从细小的黏土颗粒到巨大的石块都有，且磨圆度差，具有明显的冰川搬运特征。随着气候逐渐变暖，进入间冰期，冰川开始消退。冰川消退时，其携带的冰碛物会在不同的位置堆积下来。当冰碛物在山谷中堆积形成天然的堤坝，或者在冰川退缩后形成的低洼区域周围堆积，拦截了冰川融水以及降水等水源时，冰碛湖便得以形成。例如，在一些山谷的狭窄处，冰川消退后留下的冰碛物堆积成坝，上游的冰川融水不断汇聚，水位逐渐上升，最终形成了冰碛湖。这些冰碛湖的形状往往不规则，受冰碛物堆积形态和地形的影响，可能呈现出长条状、圆形或不规则多边形等。冰碛湖的形成还受到地形地貌的制约。喜马拉雅山地区地势起伏大，高山峡谷众多，这种复杂的地形为冰碛湖的形成提供了有利条件。在一些地势低洼的区域，更容易汇聚冰川融水和降水，从而形成冰碛湖。地质构造活动也对冰碛湖的形成有一定影响。地震等地质活动可能导致山体滑坡，增加冰碛物的数量和堆积规模，进而促进冰碛湖的形成。同时，断裂带等地质构造可能影响地下水的分布，为冰碛湖提供额外的水源补给。2.2溃决灾害的诱发因素2.2.1自然因素冰崩与雪崩：喜马拉雅山地区的冰川在重力、温度变化以及地震等因素的影响下，冰体结构会变得不稳定，容易发生冰崩。当大量冰体突然崩落并坠入冰碛湖时，会产生巨大的冲击力和涌浪。如果涌浪高度超过冰碛坝的高度，就会导致湖水漫顶，进而引发冰碛湖溃决。雪崩也可能带来大量的积雪和碎屑物质，冲入冰碛湖后同样会引起湖水位急剧上升，增加溃决风险。在喜马拉雅山的某些区域，由于冰川表面的裂缝不断扩大，冰体在重力作用下突然崩塌，大量冰块落入冰碛湖，瞬间激起数米高的涌浪，对冰碛坝造成严重威胁。相关研究表明，冰崩引发的冰碛湖溃决事件在该地区的溃决案例中占比较高，是导致溃决的重要自然因素之一。地震：喜马拉雅山处于板块碰撞的强烈构造活动带上，地震频繁发生。地震产生的强烈震动会破坏冰碛坝的结构，使其稳定性降低。地震引发的山体滑坡、崩塌等次生灾害，还会导致大量土石等物质进入冰碛湖，使湖水位迅速上升。当地震的震级较高且震中靠近冰碛湖时，冰碛坝可能会出现裂缝、塌陷等情况，无法承受湖水的压力，从而引发溃决。例如，在历史上的某些地震事件后，冰碛湖周边地区出现了明显的山体滑坡，大量的滑坡物质涌入冰碛湖，短时间内湖水位急剧上升，最终导致冰碛湖溃决，给下游地区带来了严重的灾害。降雨：强降雨是引发冰碛湖溃决的重要因素之一。当喜马拉雅山地区遭遇持续的强降雨时，大量的雨水会迅速汇入冰碛湖，使湖水位在短时间内大幅上升。强降雨还可能导致冰碛坝的土体饱和，降低坝体的抗剪强度，增加坝体滑坡和坍塌的风险。如果此时冰碛坝的排水能力不足，无法及时排出过多的湖水，就容易引发溃决。在雨季，一些冰碛湖会因为强降雨导致水位迅速超过冰碛坝的安全高度，湖水漫溢，冲毁冰碛坝，引发溃决灾害。相关统计数据显示，在部分年份的雨季，冰碛湖溃决事件的发生频率明显增加，与强降雨的关联性显著。冰川融化：随着全球气候变暖，喜马拉雅山地区的气温持续升高，导致冰川加速融化。冰川融水是冰碛湖的主要补给水源之一，冰川融化速度的加快使得冰碛湖的水量不断增加，湖水位持续上升。长期的冰川融化还可能导致冰碛坝内的埋藏冰融化，坝体结构被削弱，稳定性降低。当冰碛湖的水位超过冰碛坝的承受能力，或者冰碛坝因结构破坏无法阻挡湖水时，就会发生溃决。据观测，近年来喜马拉雅山地区的一些冰碛湖面积不断扩大，水位持续上升，这与冰川加速融化密切相关，也预示着冰碛湖溃决的风险在不断增加。2.2.2人为因素工程建设：在喜马拉雅山地区进行的公路、铁路、水电等工程建设活动，可能会对冰碛湖的稳定性产生影响。在冰碛湖周边进行的工程开挖、爆破等作业，可能会破坏冰碛坝的结构，使其抗滑和抗渗能力下降。工程建设过程中产生的弃渣如果随意堆放在冰碛湖附近，遇到暴雨等情况时，这些弃渣可能会被冲入冰碛湖，导致湖水位上升，增加溃决风险。例如，在一些公路建设项目中，为了开辟道路，对冰碛湖周边的山体进行开挖，破坏了冰碛坝的部分结构，使得冰碛坝在后续的降雨和湖水压力作用下，出现了裂缝和滑坡，增加了冰碛湖溃决的隐患。水资源开发：对喜马拉雅山地区水资源的不合理开发利用，也可能影响冰碛湖的稳定性。过度抽取冰碛湖的湖水用于灌溉、工业用水等，会导致湖水位下降，使冰碛坝的浸润线发生变化，坝体的稳定性受到影响。在冰碛湖上游修建水库等水利设施，改变了冰碛湖的来水条件，如果调度不当，可能会导致冰碛湖的水量突然增加，引发溃决。一些地区为了发展农业，大量抽取冰碛湖的湖水进行灌溉，使得冰碛湖的水位明显下降，冰碛坝的部分区域暴露在空气中，在后续的气候和地质作用下，坝体出现了干裂和松动，增加了冰碛湖溃决的风险。2.3溃决过程与危害当冰碛湖溃决时，其湖水倾泻与洪水演进过程极为迅猛且复杂。一旦冰碛湖的溃决条件满足，如冰碛坝因各种诱发因素失去稳定性，湖水就会瞬间突破冰碛坝的束缚，以极高的速度向外倾泻。在最初阶段，湖水会形成强大的水头，如同决堤的洪流，对冰碛坝剩余部分以及周边的地形产生巨大的冲击力，进一步破坏冰碛坝结构，使溃口迅速扩大。随着溃口的扩大，湖水的流量急剧增加，形成汹涌的洪水。洪水在演进过程中，受地形地貌的影响显著。在狭窄的山谷中，洪水会加速下泄，水位迅速抬升，形成高耸的洪峰，流速可达到数米每秒甚至更高，其携带的巨大能量能够轻易地冲毁沿途的一切障碍物。当洪水进入相对开阔的区域时，流速虽会有所减缓，但洪水会迅速扩散，淹没大面积的土地。在这个过程中，洪水会裹挟着大量的泥沙、石块等物质，这些物质随着洪水的流动不断撞击沿途的物体，增加了洪水的破坏力。相关研究表明，冰碛湖溃决引发的洪水，其洪峰流量往往是正常河流流量的数倍甚至数十倍，如1981年次仁玛错冰碛湖溃决时，洪峰流量高达1000立方米每秒以上，远远超过了下游河流的承受能力。冰碛湖溃决对下游生命财产安全构成了巨大威胁。在人口密集的地区，溃决引发的洪水可能在短时间内淹没村庄、城镇，导致大量居民被困，造成人员伤亡。许多居民来不及撤离，被洪水冲走，生命安全受到严重威胁。洪水还会冲毁房屋、桥梁、道路等基础设施，使得交通中断，救援物资和人员难以进入受灾地区，进一步加剧了灾害的影响。在一些山区，由于道路被冲毁，救援工作往往会延迟数小时甚至数天，给被困人员的生存带来了极大的挑战。冰碛湖溃决对生态环境也会产生深远的负面影响。洪水携带的大量泥沙和杂物会导致下游河流的水质恶化，水中的含沙量急剧增加，影响水生生物的生存环境，导致许多水生生物死亡或迁移。溃决还可能引发山体滑坡、泥石流等次生地质灾害，进一步破坏生态系统的稳定性。大量的泥沙和石块堆积在下游地区，改变了原有的地形地貌，破坏了植被，导致水土流失加剧，土壤肥力下降，影响了当地的农业生产和生态平衡。三、冰碛湖溃决灾害评价方法3.1传统评价方法概述在冰碛湖溃决灾害评价领域，传统评价方法主要包括定性评价方法和定量评价方法，它们在冰碛湖溃决灾害研究的历史进程中发挥了重要作用，为后续研究奠定了基础。定性评价方法侧重于从冰碛湖的相关特征和影响因素进行分析判断。徐道明等通过对我国西藏地区溃决冰碛湖的考察，提出从坝体特征（如坝体高度、宽度、物质组成等）、气候因素（如气温、降水变化等）判别危险冰碛湖的几类指标。吕儒仁等分析了我国西藏冰碛湖溃决特点，统计了有利于溃决事件发生的7类指标，涵盖了冰碛湖与母冰川的关系、冰碛坝的结构等方面。Huggel等利用遥感影像从冰碛湖、坝、激发条件、下游洪水路径等方面提出冰碛湖潜在危险性评判的三级准则和相应指标，建立了较为系统的定性评估体系。这些定性评价方法的优点在于能够直观地反映冰碛湖溃决的一些关键因素，操作相对简单，不需要复杂的数学计算和大量的数据支持。通过对坝体物质组成的观察，就可以初步判断坝体的稳定性。但定性评价方法也存在明显的局限性，其评价结果往往依赖于研究者的经验和主观判断，缺乏量化的标准，不同研究者可能会得出不同的结论，难以进行精确的比较和分析。定量评价方法则试图通过数学模型和具体的数据指标来量化冰碛湖溃决的风险。吕儒仁等将冰碛湖与母冰川变坡点（或裂隙面）下游冰体的体积比定义为危险性指数，提出了评估由冰崩涌浪造成坝体溃决的定量化指标，为冰崩涌浪引发的溃决风险评估提供了具体的量化依据。陈晓清等从喜马拉雅地区冰碛湖相关因子中选取7个参数作为直接判别的依据，结合冰碛湖溃决危险指数建立起一种定性-半定量的个体冰碛湖危险性评价方法，在一定程度上实现了对冰碛湖溃决风险的量化评估。舒有锋等从冰川、湖盆、终碛堤、气候等特征中提取冰碛湖溃决危险性评价指标，通过粗糙集理论确定各指标权重，建立了冰碛湖溃决危险性评价方法，使评价结果更具客观性和科学性。定量评价方法的优势在于能够通过具体的数据和模型，较为精确地评估冰碛湖溃决的风险程度，便于不同冰碛湖之间的风险比较和排序。但这种方法对数据的依赖性很强，需要大量准确的数据来支撑模型的建立和运算。而在实际情况中，喜马拉雅山地区的冰碛湖往往地处偏远，获取数据的难度较大，数据的准确性和完整性也难以保证。定量评价方法中的模型往往具有一定的假设条件和局限性，难以完全准确地反映冰碛湖溃决这一复杂的自然过程。3.2多态贝叶斯网络（B-Ns）评价方法原理3.2.1B-Ns基本原理多态贝叶斯网络（B-Ns）是一种基于概率推理的图形模型，它巧妙地结合了概率论与图论的知识，能够有效处理不确定性问题，在复杂系统的分析和预测中具有独特优势。从结构上看，B-Ns由节点和有向边组成，节点代表随机变量，有向边则表示变量之间的因果关系。每个节点都有一个条件概率表（CPT），用于描述该节点在其所有父节点不同状态组合下的概率分布。在实际应用中，B-Ns的推理过程基于贝叶斯定理。贝叶斯定理是概率论中的一个重要定理，它描述了在已知一些先验信息的情况下，如何更新对某个事件的概率估计。对于B-Ns来说，当给定某些节点的观测值（证据）时，通过贝叶斯推理可以计算出其他节点的后验概率，从而实现对系统状态的推断和预测。假设我们有一个简单的B-Ns模型，其中节点A是节点B的父节点，已知节点A处于某种状态的先验概率，以及在节点A不同状态下节点B的条件概率。当我们观测到节点B的实际状态时，就可以利用贝叶斯定理计算出节点A在这个证据下的后验概率，进而更准确地了解整个系统的状态。与其他模型相比，B-Ns在处理不确定性问题上具有显著优势。传统的确定性模型往往假设系统中的变量是完全确定的，然而在现实世界中，尤其是在自然灾害研究领域，如冰碛湖溃决灾害，存在大量的不确定性因素。冰碛湖溃决的发生受到多种因素的综合影响，这些因素之间的关系复杂且具有不确定性，很难用传统的确定性模型进行准确描述。B-Ns能够充分考虑这些不确定性，通过概率分布来表示变量的不确定性，以及变量之间关系的不确定性。它可以融合多源信息，包括历史数据、专家知识等，从而更全面地描述系统的状态和行为，提高预测和分析的准确性。在确定冰碛湖溃决风险时，B-Ns可以将冰碛坝的稳定性、冰川的变化、气象条件等多种因素作为节点纳入网络，并根据历史数据和专家经验确定这些节点之间的概率关系，从而更准确地评估冰碛湖溃决的可能性。3.2.2构建针对冰碛湖溃决的B-Ns模型构建针对冰碛湖溃决的B-Ns模型，首要任务是明确节点变量。通过对青藏高原南缘74例冰碛湖溃决灾害实例的深入统计分析可知，冰碛坝失效的主要模式可归纳为5种，将其作为一级指标，分别为冰崩涌浪漫顶、管涌、洪水漫顶、冰碛坝坍塌和外界动力激发。在此基础上，结合冰碛湖溃决灾害形成的基本条件，进一步为每个一级指标确定各自的二级分析指标。对于冰崩涌浪漫顶模式（y1），其影响因子众多。危险冰体坡度（x1）对冰体崩落起动的难易程度起着关键作用，统计数据表明，当冰川末端的危险冰体坡度超过8°时，冰崩灾害发生的可能性显著增加。危险冰体指数（x2），即危险冰体体积与冰碛湖湖水体积的比值，该指数越大，意味着冰体滑入冰碛湖后产生的涌浪和水位骤涨综合形成的水头高度越大，坝体溃口处起动泥沙颗粒越容易，相应的溃决可能性也就越高，参照冰碛湖溃决案例，溃决前危险冰体指数一般位于0.11-0.52。冰舌冰碛湖距离（x3）也是重要影响因素，冰崩涌浪的首要形成条件是冰川与冰碛湖处于邻近状态，已溃决的冰碛湖此距离大多小于1km。冰川裂隙发育程度（x4）与冰崩灾害的易发性呈正相关关系，在气候反常年份，冰川容易在裂隙发育面发生断裂。年际气温变异系数（x5）可代表气温突变程度，系数值越大，气温年际变化越剧烈，越有利于冰崩灾害的形成，冰碛湖溃决区的气温变异系数值通常位于0.23-0.27。管涌模式（y2）主要受坝体厚度和坝体物质组成的影响。目前对已溃决冰碛湖的统计显示，大部分冰碛坝顶宽（x6）不超过60m。当冰碛坝体内部含有埋藏冰时，随着气温升高，埋藏冰融化容易引发冰碛坝的管涌现象，虽然当前缺乏有效判别冰碛坝内是否存在埋藏冰的方法，但研究认为在小冰期末时期所形成的冰碛堰塞湖的坝体中，埋藏冰普遍存在且受气温变化影响明显，因此可以采用冰碛坝顶宽、坝体物质状态（x7）与年际气温变异系数对管涌的危险性进行分析。洪水漫顶模式（y3）通常在气候异常年份发生，当冰川融水与降雨导致的超标洪水对冰碛湖的补给超过溢流、蒸发的损耗时，冰碛湖容易发生漫流型溃决。采用年际降雨变异系数表征降雨的不均匀度，冰碛湖漫流溃决影响因子可以采用冰碛湖满溢程度（x8）、年际降雨变差系数（x9）与年际气温变异系数表示。冰碛坝坍塌模式（y4）中，冰碛坝多为松散岩屑或积雪（冰）夹带岩屑的稳定性较差的自然堆积体。坝体的坡度越大，稳定性越差，在内、外部因素作用下容易坍塌，对青藏高原南缘冰碛湖溃决案例的调查表明，溃决冰碛堰塞坝背水面的坡度多位于20°-33°。坝顶若存在漫流下切槽道，则坝体容易在溯源侵蚀的作用下坍塌，此外，冰碛坝内埋藏冰消融也可能导致坝体坍塌，因此，可以采用背水坡度（x10）与坝顶漫流槽道分布状态（x11）以及年际气温变异系数进行分析。外界动力激发模式（y5）主要由地震作用引起。地震力一方面促使冰碛湖附近岩体、冰雪崩滑激发涌浪，另一方面地震直接激发水体涌浪导致漫顶现象的发生。研究发现，地震峰值加速度（PGA，x12）与涌浪波高正相关，当x12>0.25g（超越概率等于10%）时，涌浪波高可达到冰碛坝颗粒起动临界条件，可能引发冰碛湖溃决灾害。在确定了上述5种溃坝机制作为中间节点变量后，依据节点变量之间的因果关系建立B-Ns拓扑结构。其中，叶节点事件T表示冰碛堰塞湖溃决，yj为中间节点事件（j=1,2,・・・,m），xi为根节点事件（i=1,2,・・・,n）。为了便于分析和计算，依据有利于形成灾害的节点状态分界值，对B-Ns根节点变量进行离散化处理，得到各根节点变量状态及取值范围，状态值0、1分别表示安全、危险状态。传统的事故树模型以及二态B-Ns系统仅将节点状态简化为“正常”和“失效”两种，但在实际中中间节点与叶节点往往存在着中间状态。以冰碛湖溃坝机制中占比最高的冰崩涌浪漫顶为例，漫顶涌浪往往需要达到一定的水头高度才能使得堰塞坝溃口处粗颗粒泥沙连续起动形成溃决连锁反应。因此，根据事故严重程度，将中间节点y1的状态分为低、中、高危险3类，状态值分别用0、1、2表示；其余溃决机制（y2-y5）由于发生几率较小，为降低贝叶斯网络复杂度，仍将其简化为两个状态（以0、1表示）。对于叶节点T，从堰塞湖溃决物理过程角度，其溃决方式可以分为全溃与部分溃，其中全溃形成的相对洪峰流量最大，危险性也最高，根据溃决危险性严重程度，将叶节点T分为低、中、高3种危险状态。3.3节点变量与概率分布确定3.3.1节点变量选择与定义在构建针对冰碛湖溃决的多态贝叶斯网络（B-Ns）模型中，节点变量的选择与定义至关重要，它们直接影响着模型对冰碛湖溃决风险评估的准确性和可靠性。基于对青藏高原南缘74例冰碛湖溃决灾害实例的深入剖析，确定了5种主要的冰碛坝失效模式作为中间节点变量，这些模式涵盖了冰碛湖溃决的主要机制，反映了冰碛湖溃决过程中复杂的因果关系。冰崩涌浪漫顶模式（y1）中，危险冰体坡度（x1）是一个关键的根节点变量。它决定着冰体崩落起动的难易程度，从统计数据来看，当冰川末端的危险冰体坡度超过8°时，冰崩灾害发生的可能性显著增加。这是因为较大的坡度会使冰体在重力作用下更容易失去平衡，从而引发冰崩。危险冰体指数（x2）同样重要，它是危险冰体体积与冰碛湖湖水体积的比值。当危险冰体滑入冰碛湖后，会产生涌浪并导致冰碛湖水位骤涨，二者综合形成的水头高度越大，坝体溃口处起动泥沙颗粒就越容易，相应的溃决可能性也就越高。参照冰碛湖溃决案例，溃决前危险冰体指数一般位于0.11-0.52，这一范围为评估冰崩涌浪漫顶模式下的溃决风险提供了重要的量化依据。冰舌冰碛湖距离（x3）也是影响冰崩涌浪漫顶模式的重要因素，冰崩涌浪的首要形成条件是冰川与冰碛湖处于邻近状态，已溃决的冰碛湖此距离大多小于1km，较短的距离使得冰体崩落能够更直接地影响冰碛湖，增加涌浪和溃决的风险。冰川裂隙发育程度（x4）与冰崩灾害的易发性呈正相关关系，在气候反常年份，冰川容易在裂隙发育面发生断裂，从而引发冰崩。年际气温变异系数（x5）可代表气温突变程度，系数值越大，气温年际变化越剧烈，越有利于冰崩灾害的形成，冰碛湖溃决区的气温变异系数值通常位于0.23-0.27，反映了气温变化对冰崩的影响。管涌模式（y2）主要受坝体厚度和坝体物质组成的影响。目前对已溃决冰碛湖的统计显示，大部分冰碛坝顶宽（x6）不超过60m，较窄的坝顶宽度会影响坝体的稳定性，增加管涌发生的可能性。当冰碛坝体内部含有埋藏冰时，随着气温升高，埋藏冰融化容易引发冰碛坝的管涌现象。虽然当前缺乏有效判别冰碛坝内是否存在埋藏冰的方法，但研究认为在小冰期末时期所形成的冰碛堰塞湖的坝体中，埋藏冰普遍存在且受气温变化影响明显，因此可以采用冰碛坝顶宽、坝体物质状态（x7）与年际气温变异系数对管涌的危险性进行分析。坝体物质状态（x7）可以反映坝体的组成成分和结构特征，不同的物质状态对管涌的敏感性不同，例如，含有较多细颗粒物质的坝体更容易发生管涌。洪水漫顶模式（y3）通常在气候异常年份发生，当冰川融水与降雨导致的超标洪水对冰碛湖的补给超过溢流、蒸发的损耗时，冰碛湖容易发生漫流型溃决。采用年际降雨变异系数表征降雨的不均匀度，冰碛湖漫流溃决影响因子可以采用冰碛湖满溢程度（x8）、年际降雨变差系数（x9）与年际气温变异系数表示。冰碛湖满溢程度（x8）直接反映了冰碛湖的水位与坝顶的相对关系，当满溢程度较高时，洪水漫顶的风险就会增加。年际降雨变差系数（x9）体现了降雨在不同年份之间的变化程度，较大的变差系数意味着降雨的不稳定性增加，更容易引发超标洪水。冰碛坝坍塌模式（y4）中，冰碛坝多为松散岩屑或积雪（冰）夹带岩屑的稳定性较差的自然堆积体。坝体的坡度越大，稳定性越差，在内、外部因素作用下容易坍塌，对青藏高原南缘冰碛湖溃决案例的调查表明，溃决冰碛堰塞坝背水面的坡度多位于20°-33°。坝顶若存在漫流下切槽道，则坝体容易在溯源侵蚀的作用下坍塌，此外，冰碛坝内埋藏冰消融也可能导致坝体坍塌，因此，可以采用背水坡度（x10）与坝顶漫流槽道分布状态（x11）以及年际气温变异系数进行分析。背水坡度（x10）影响着坝体在水压力作用下的稳定性，较大的背水坡度会使坝体更容易发生滑动和坍塌。坝顶漫流槽道分布状态（x11）反映了坝顶水流的冲刷情况，槽道越多、越深，坝体在溯源侵蚀作用下坍塌的风险就越高。外界动力激发模式（y5）主要由地震作用引起。地震力一方面促使冰碛湖附近岩体、冰雪崩滑激发涌浪，另一方面地震直接激发水体涌浪导致漫顶现象的发生。研究发现，地震峰值加速度（PGA，x12）与涌浪波高正相关，当x12>0.25g（超越概率等于10%）时，涌浪波高可达到冰碛坝颗粒起动临界条件，可能引发冰碛湖溃决灾害。地震峰值加速度（x12）是衡量地震强度的重要指标，较大的峰值加速度意味着更强的地震作用力，对冰碛湖和冰碛坝的稳定性产生更大的影响。3.3.2概率分布获取方法在构建多态贝叶斯网络（B-Ns）模型时，确定节点变量的概率分布是关键步骤之一，它直接影响模型的准确性和可靠性。节点概率的获取方法主要有资料统计法和专家知识法两种，这两种方法各有优势，将它们结合使用可以更全面、准确地确定节点概率分布。根节点先验概率可依据资料统计法获得。通过对所研究区域的冰川、冰碛湖、遥感影像与数字高程模型（DEM）地形等资料的分析，对根节点因子落入各状态区间频数进行统计，从而确定根节点各状态的先验概率。例如，在确定危险冰体坡度（x1）的先验概率时，利用研究区的遥感影像和DEM地形数据，测量大量冰川末端危险冰体的坡度，并统计不同坡度区间内的冰碛湖数量，以此计算出不同坡度状态下的概率。在获取冰碛坝顶宽（x6）的先验概率时，对遥感影像（2014年-2019年）与DEM地形（2009年）进行解译，测量已溃决冰碛湖的坝顶宽度，统计不同宽度区间的频数，进而确定各状态的先验概率。本文所使用统计资料中，冰川、冰湖数据来自国际山地综合发展中心（ICIMOD，2015年）与中科院寒旱所（2014年）提供的第二次冰川编目，这些权威的数据来源为准确统计根节点先验概率提供了有力支持。气候数据来自中国气象数据网提供的研究区气象站近30年气象资料，利用这些气象数据，可以统计年际气温变异系数（x5）、年际降雨变差系数（x9）等气候相关根节点变量在不同状态下的概率。在获取根节点先验概率之后，还需要确定其他节点（中间节点和叶节点）的条件概率表（CPT）。节点处CPT可利用模糊数学法将专家的定性描述转化为定量的发生概率。例如，对于冰崩涌浪漫顶模式（y1）这一中间节点，邀请多位在冰川灾害研究领域具有丰富经验的专家，让他们根据自己的专业知识和实际经验，对在不同根节点状态组合下y1处于低、中、高危险状态的可能性进行定性判断。然后，运用模糊数学法将这些定性判断转化为具体的概率值，构建y1的条件概率表。针对专家知识法获得的CPT存在主观性的问题，可以将获取的条件概率与实际案例互相验证，以降低主观性。以1981年7月11日波曲河流域次仁玛错冰碛湖溃决案例为例，根据溃决前该冰碛湖的实际状态，更新贝叶斯网络中各根节点状态量，结合利用专家知识法构建的中间节点条件概率进行推理计算，得到叶节点T（冰碛堰塞湖溃决）处于不同危险状态的概率。将计算结果与实际溃决情况进行对比，如果计算结果与实际情况相符，则说明构建的条件概率表具有一定的可靠性；如果存在偏差，则进一步分析原因，调整条件概率表，直到计算结果与实际情况较为吻合。通过这种方式，可以不断优化条件概率表，提高模型的准确性。3.4冰碛湖溃决概率计算与危险性分级在多态贝叶斯网络（B-Ns）模型中，冰碛湖溃决概率的计算是基于根节点的先验概率分布以及各节点之间的条件概率关系。根据B-Ns的运算规则，从根节点的先验概率出发，可以正向推算叶节点事件（即冰碛堰塞湖溃决）的发生概率。在多态B-Ns中，设根节点为x_i，中间节点为y_j，叶节点为T，a_i、b_j和q分别代表根节点、中间节点、叶节点的风险状态，其中a_i=0,1,\cdots,u_i−1；b_j=0,1,\cdots,v_j−1；q=0,1,\cdots,r−1，u_i、v_j、r分别为根节点、中间节点、叶节点的状态数。当根节点各状态的概率分别为P(x_{i}=a_{i})时，叶节点T处于风险状态T_q的概率计算式为：P(T=T_q)=\sum_{b_1=0}^{v_1-1}\cdots\sum_{b_m=0}^{v_m-1}\left[P(T=T_q|\pi(T))\prod_{j=1}^{m}P(y_j=b_j|\pi(y_j))\prod_{i=1}^{n}P(x_i=a_i)\right]式中，\pi(y_j)为中间节点y_j的父节点集合；\pi(T)为叶节点T的父节点集合。这个公式综合考虑了所有根节点和中间节点对叶节点状态的影响，通过对各个节点不同状态组合下的概率进行求和，得到叶节点处于特定风险状态的概率。依据待评估冰碛湖的证据信息，利用B-Ns的正向推理，将根节点的先验概率和中间节点的条件概率代入上述公式，即可计算叶节点T各状态的发生概率。为了确定冰碛湖的危险性等级，需要进一步确定各风险状态的隶属度函数。这里采用最大隶属度原则，即通过比较叶节点处于不同风险状态的概率大小，确定叶节点隶属的状态。例如，若计算得到P(T=T_0)=0.1，P(T=T_1)=0.3，P(T=T_2)=0.6，由于P(T=T_2)最大，则该冰碛湖的溃决危险性评价结果为高危险状态。根据溃决危险性严重程度，将叶节点T分为低、中、高3种危险状态。具体的分级标准可以根据实际情况和研究目的进行设定。一般来说，当叶节点处于低危险状态时，冰碛湖溃决的概率相对较低，可能在未来很长一段时间内都不会发生溃决，对下游地区的威胁较小；处于中危险状态时，冰碛湖溃决的可能性有所增加，需要密切关注其变化，加强监测和预警；而处于高危险状态时，冰碛湖溃决的概率较高，随时可能发生溃决，对下游地区的生命财产安全构成严重威胁，需要立即采取有效的防治措施，如疏散下游居民、修建防护堤坝等。四、评价方法的验证与应用4.1模型验证4.1.1选择验证案例选择1981年波曲河流域次仁玛错冰碛湖溃决事件作为验证案例，具有多方面的考量因素。该事件在冰碛湖溃决灾害研究领域具有较高的知名度和典型性，是研究冰碛湖溃决机制和影响的重要案例之一，相关的研究资料和数据较为丰富，为模型验证提供了坚实的数据基础。次仁玛错冰碛湖位于喜马拉雅中部波曲河流域，该地区是冰碛湖溃决灾害的高发区域，其地理环境和气候条件具有显著的代表性。喜马拉雅山地区作为中低纬度地区最大的现代冰川分布区，冰川活动频繁，冰碛湖众多，而波曲河流域在这一区域中，无论是冰川的规模、冰碛湖的形成条件还是气候的变化特征，都与其他地区具有相似性。选择该地区的冰碛湖溃决事件进行验证，能够更好地反映出模型在喜马拉雅山地区的适用性和有效性，为其他冰碛湖的风险评估提供可靠的参考。从溃决机制来看，1981年次仁玛错冰碛湖溃决是多种因素共同作用的结果。根据吕儒仁等的调查，溃决当年喜马拉雅南坡的气候背景为暖干气候，降水因量小而作用很小，突然升温导致冰崩，叠加管涌效应导致冰碛湖溃决。这种复杂的溃决机制涵盖了多态贝叶斯网络（B-Ns）模型中所考虑的多个关键因素，如冰崩涌浪漫顶模式中的危险冰体坡度、危险冰体指数、冰舌冰碛湖距离等因素，以及管涌模式中的坝体厚度、坝体物质组成等因素。通过对这一案例的验证，可以全面检验B-Ns模型在处理多种溃决机制共同作用下冰碛湖溃决概率计算的准确性和可靠性。4.1.2验证过程与结果分析按照B-Ns模型的评价流程，首先需要确定节点概率分布。根节点各风险状态区间的先验概率依据资料统计法获得。根据所研究区域的冰川、冰碛湖、遥感影像与数字高程模型（DEM）地形对根节点因子落入各状态区间频数进行统计，从而确定根节点各状态的先验概率。本文所使用统计资料中，冰川、冰湖数据来自国际山地综合发展中心（ICIMOD，2015年）与中科院寒旱所（2014年）提供的第二次冰川编目；冰碛坝等几何与结构参数源自对遥感影像（2014年-2019年）与DEM地形（2009年）的解译；气候数据来自中国气象数据网提供的研究区气象站近30年气象资料。通过这些权威的数据来源，统计得到次仁玛错冰碛湖溃决事件中根节点各状态的先验概率。依据各风险因素间的因果关系，结合案例资料与专家知识构造中间节点y1-y5及叶节点T各状态区间的条件概率分布。以中间节点y1（冰崩涌浪漫顶模式）为例，其条件概率表（CPT）的构建过程如下：邀请多位在冰川灾害研究领域具有丰富经验的专家，让他们根据自己的专业知识和实际经验，对在不同根节点状态组合下y1处于低、中、高危险状态的可能性进行定性判断。然后，运用模糊数学法将这些定性判断转化为具体的概率值，构建y1的条件概率表。针对专家知识法获得的CPT存在主观性的问题，将获取的条件概率与实际案例互相验证，以降低主观性。在确定节点概率分布后，依据溃决前次仁玛错冰碛湖的状态，更新贝叶斯网络中各根节点状态量。溃决前夕冰舌到达海拔4700m位置，距离冰碛湖非常近，这表明冰舌冰碛湖距离（x3）这一根节点处于危险状态；冰川前缘地形为陡崖状，且在陡崖处呈不连续状态，裂纹发育，说明危险冰体坡度（x1）和冰川裂隙发育程度（x4）也处于危险状态；溃决之前，冰碛湖处于满溢状态，坝顶槽道发育，已有湖水从坝体中连续渗出，这对冰崩涌浪漫顶和管涌模式都有影响，相应的根节点状态也需更新。结合中间节点的条件概率进行推理计算，可得叶节点T（冰碛堰塞湖溃决）各状态的概率。假设计算结果为P(T=0)=0.0670，P(T=1)=0.2381，P(T=2)=0.6949。根据最大隶属度原则，由于P(T=2)最大，所以该冰碛湖的溃决危险性评价结果为高危险。这与1981年次仁玛错冰碛湖实际发生溃决的情况相一致，从而对本文B-Ns模型的准确性进行了进一步验证。从计算结果可以看出，B-Ns模型能够较为准确地评估次仁玛错冰碛湖的溃决危险性，反映出模型在处理复杂的冰碛湖溃决问题上具有较高的可靠性和有效性。通过对这一案例的验证，也为B-Ns模型在其他冰碛湖溃决风险评估中的应用提供了有力的支持，证明了该模型在喜马拉雅山地区冰碛湖溃决灾害评价中具有重要的应用价值。4.2应用案例分析-中尼铁路线路方案比选4.2.1研究区概况中尼铁路是拟建穿越喜马拉雅山脉的重要铁路项目，线路从日喀则出发，先在高原面上行走，而后下高原面至加德满都。在其规划过程中，面临经吉隆廊道与经樟木廊道这两个局部走向方案的比选问题，这两个廊道均位于喜马拉雅中部希夏邦马峰冰川群分布区，冰碛湖数量众多，冰碛湖溃决洪水已成为对线路方案起控制作用的重大灾害类型。吉隆流域位于喜马拉雅中部，东西走向的喜马拉雅山脉将其分割成南北两部分。北部海拔较高，地势相对平缓，多为高原地貌，这种地形条件使得冰碛湖的形成和分布具有一定的特点，由于地势开阔，冰碛物的堆积和冰川融水的汇聚相对分散。南部海拔低，但山体陡峻，高差较大，属于高山峡谷地貌。在这种地形中，冰川活动频繁，冰碛湖更容易在山谷的低洼处形成，且由于地形狭窄，冰碛湖的形态往往较为狭长。远古时期吉隆地区分布有吉隆、沃玛等湖泊，由于青藏高原的不均匀抬升，导致吉隆藏布溯源侵蚀切穿原有湖泊，湖水下泄，这样原古湖盆与现代河流形成了吉隆藏布上游段宽窄相间的河谷地貌形态。此外，由于海洋性冰川塑造地貌效应，中游段也形成了一些宽谷地形，这些宽谷和峡谷段的分布对冰碛湖溃决洪水的演进以及铁路选线都有着重要影响。在宽谷段，溃决洪水的流速相对较慢，对铁路的冲击力相对较小，但淹没范围可能较大；而在峡谷段，洪水的流速会加快，冲击力增强，对铁路的威胁更大。樟木廊道位于波曲流域，波曲河发源于西藏聂拉木县波绒乡，经聂拉木县城到樟木镇。流域内地形高差较大，最高峰为海拔超过8km的希夏邦马峰。波曲的河谷地貌以聂拉木县城为分界线，其北基本位于高原面上，主要以宽谷地貌为主，地势较为平坦，冰碛湖的分布相对稀疏，且规模较大。其南属于溯源侵蚀河段，以峡谷地貌为主，河谷狭窄，两岸山势陡峭，冰碛湖数量较多，且多分布在支沟内。这种地形条件使得冰碛湖溃决洪水在该区域的演进过程更加复杂，一旦冰碛湖溃决，洪水在峡谷中容易形成高水头，流速极快，对下游的铁路线路安全构成严重威胁。4.2.2冰碛湖溃决危险性评价在对中尼铁路沿线的冰碛湖溃决危险性进行评价时，首先利用筛选判据来确定需要进行危险性评价的冰碛湖。该筛选判据综合考虑了冰碛湖的体积以及溃决洪水到线路工程的路径长度等因素，通过借用Washakh等提出的冰碛湖溃决洪水淹没长度计算式（L=0.52V+36.13），并定义溃决洪水到线路工程的路径长度为Ld，以此筛选出可能对铁路线路造成影响的冰碛湖。然后，运用前文建立的多态贝叶斯网络（B-Ns）评价模型，对吉隆、樟木廊道流域内冰碛湖溃决危险性进行评价。在构建B-Ns模型时，明确了5种主要的冰碛坝失效模式作为中间节点变量，分别为冰崩涌浪漫顶、管涌、洪水漫顶、冰碛坝坍塌和外界动力激发，并确定了各自对应的根节点变量，如危险冰体坡度、冰碛坝顶宽、冰碛湖满溢程度等。通过资料统计法获得根节点各状态的先验概率，利用模糊数学法将专家的定性描述转化为中间节点和叶节点的条件概率，从而构建出完整的B-Ns模型。经过评价可知，吉隆流域危险性较大的冰湖总数为7个，其中高危险冰碛湖4个，中危险冰碛湖3个。在吉隆流域的高危险冰碛湖中，部分冰碛湖由于冰川末端的危险冰体坡度较大，超过8°，且危险冰体指数位于0.11-0.52之间，冰舌冰碛湖距离小于1km，同时年际气温变异系数处于0.23-0.27的范围，这些因素综合作用，使得这些冰碛湖在冰崩涌浪漫顶模式下具有较高的溃决风险。樟木流域危险性较大的冰碛湖总数为13个，其中高危险冰碛湖6个，中危险冰碛湖7个。樟木流域的一些高危险冰碛湖，除了存在与吉隆流域类似的冰崩涌浪漫顶风险因素外，还由于其所处的河谷地形狭窄，在洪水漫顶和冰碛坝坍塌模式下，溃决风险也相对较高。当遇到强降雨或冰川融水过多时，冰碛湖容易发生漫流型溃决，而狭窄的河谷会加剧洪水的冲击力，对下游铁路构成更大威胁。樟木流域的冰碛坝由于多为松散岩屑或积雪（冰）夹带岩屑的自然堆积体，稳定性较差，在背水坡度较大以及坝顶漫流槽道发育的情况下，更容易发生坍塌，引发冰碛湖溃决。4.2.3铁路选线方案评价中尼铁路跨喜马拉雅段的樟木、吉隆局部走向方案各具特点。从线路长度来看，樟木方案在冰碛湖分布区的线路长度小于吉隆方案，从选线的角度，樟木方案属于线路短直方案，在常规的铁路选线考量中，短直的线路方案通常具有工程成本相对较低、线路运营效率较高等优势。在考虑冰碛湖溃决风险这一关键因素时，情况发生了变化。依据前文建立的铁路方案冰碛湖溃决危险性评价指标对两方案进行比选，结果显示樟木方案的高、次高、中风险区段线路长度都远超过吉隆方案。樟木方案的高风险区段线路长度较长，主要是因为该方案经过的区域冰碛湖数量较多，且部分冰碛湖的溃决危险性较高。如前文所述，樟木流域危险性较大的冰碛湖总数为13个，其中高危险冰碛湖6个，中危险冰碛湖7个，这些危险冰碛湖分布较为集中，一旦发生溃决，洪水可能会对较长路段的铁路造成影响。而吉隆方案虽然线路在冰碛湖分布区相对较长，但由于其经过的区域冰碛湖分布相对分散，且高危险冰碛湖数量相对较少，仅为4个，所以在冰碛湖溃决风险方面相对较低。从冰碛湖溃决风险角度综合分析，吉隆方案优于樟木方案。这表明在铁路选线过程中，不能仅仅考虑线路的长度等常规因素，对于像冰碛湖溃决这样的重大灾害风险因素必须给予充分的重视。在冰碛湖溃决危险区进行铁路选线时，应优先选择经过冰碛湖溃决风险较低区域的线路方案，即使该方案的线路长度相对较长。因为一旦冰碛湖溃决，对铁路造成的破坏可能是灾难性的，不仅会导致铁路设施的损毁，还可能影响铁路的正常运营，造成巨大的经济损失和社会影响。在后续的铁路选线和建设过程中，需要进一步结合其他因素，如工程地质条件、生态环境影响等，对线路方案进行综合评估和优化，以确保铁路建设的安全性和可持续性。五、结论与展望5.1研究成果总结本研究成功建立了基于多态贝叶斯网络（B-Ns）的冰碛湖溃决灾害评价方法。通过对青藏高原南缘74例冰碛湖溃决灾害实例的深入统计分析，明确了5种主要的冰碛坝失效模式作为中间节点变量，分别为冰崩涌浪漫顶、管涌、洪水漫顶、冰碛坝坍塌和外界动力激发，并确定了各自对应的根节点变量，如危险冰体坡度、冰碛坝顶宽、冰碛湖满溢程度等。这些节点变量全面涵盖了冰碛湖溃决过程中的关键影响因素，通过对这些因素