从灾害到资源：喜马拉雅中段波曲河流域冰湖水电潜力和开发适宜性的流域尺度筛选评估

《Water Resources Research》：From Hazard to Resource: A Screening-Level Assessment of Glacial Lake Hydropower Potential and Development Suitability in the Poiqu River Basin, Central Himalaya

　　

气候变暖加速冰川退缩和冰湖扩张，在增加洪水灾害的同时也创造了水电开发机遇。尽管针对冰湖灾害的研究已十分广泛，但关于冰湖水电潜力和开发适宜性的评估仍然有限。本研究通过整合多情景气候预估与年度湖泊水量收支重建，估算了喜马拉雅中段波曲河流域37个冰湖的理论水电潜力，并据此计算2025–2100年四个运行阶段中与水电相关的入流量；同时，通过分级河网与地形梯度分析推求有效水头。研究人员构建了一个多准则框架，用以评估环境、技术和经济适宜性并对湖泊进行优先级排序。在四种气候情景下，流域尺度年平均发电量在第一阶段（风险减缓期，2025–2030）约为212.8 GWhr yr^?1，随后在第二至第四阶段（适应性发电期，2031–2100）分别约为200.6、197.0和194.3 GWhr yr^?1。帕库普科（Paqupuco）被确定为最高优先级场址，年发电潜力约为14.3 GWhr yr^?1，适宜性评分为80.6。在区域筛选尺度上，模拟运行可将高风险湖泊（如岗西错Gangxico和嘎龙错Galongco）的最大库容削减至其初始库容的约71%。敏感性分析表明，在所测试的水量收支与权重假设下，流域尺度估算值和主要候选湖泊的结果总体稳定。通过将水电资源评估与气候适应性灾害治理相结合，本研究提供了一个区域筛选框架，用于在冰川作用山区联合评估清洁能源机遇与冰湖溃决洪水缓解措施。

研究背景与问题：全球气候变暖加速冰川退缩与冰湖扩张，这一方面增加了冰湖溃决洪水（GLOF）风险，另一方面也使冰湖成为潜在的水电资源。然而，现有研究多聚焦于冰湖灾害评估，缺乏对冰湖水电潜力与开发适宜性的系统性、湖盆尺度和流域尺度联合评估。波曲河流域位于喜马拉雅中段中尼边境，面积约2,150 km²，拥有106个面积大于0.01 km²的冰湖，其中多数为终碛湖，GLOF风险较高，同时该区域地形高差大、冰川径流丰富，具备水电开发的潜在条件。为此，研究人员开展了一项流域尺度筛选评估，旨在量化37个面积大于0.01 km²的冰湖的理论水电潜力、开发适宜性及其在GLOF风险减缓中的潜在作用。

研究内容与结论：研究人员整合多源地理空间与气候数据，采用多情景气候预估和年度湖泊水量收支重建方法，估算了2025–2100年四个运行阶段的冰湖水电相关入流量；通过分级河网与地形梯度分析推求有效水头，并基于经典水电公式计算理论发电量。研究进一步构建了包含环境、技术和经济三个维度的多准则适宜性评价框架，并结合K-means聚类建立优先度矩阵。结果表明：四种共享社会经济路径（SSP126、SSP245、SSP370、SSP585）情景下，流域尺度年平均发电量在第一阶段（2025–2030，风险减缓）约为212.8 GWhr yr^?1，随后在第二阶段（2031–2050）、第三阶段（2051–2075）和第四阶段（2076–2100）分别降至200.6、197.0和194.3 GWhr yr^?1。冰川融水径流占净水电相关入流量的83.2%，非冰川区径流占12.8%，湖面降水占5.0%，湖面蒸发为1.0%的损失项。7个湖泊的年均发电量超过10 GWhr yr^?1，其中查玛曲丹错（Chamaqudanco）最高（17.9 GWhr yr^?1），帕库普科（Paqupuco）以14.3 GWhr yr^?1的年发电量和80.6的适宜性评分被确定为最高优先级场址。5个湖泊属于高潜力中等适宜性类别，2个湖泊属于中等潜力高适宜性类别。此外，模拟运行表明，高风险湖泊如岗西错（Gangxico）和嘎龙错（Galongco）的最大库容可降至初始库容的约71%，显示出水电运行与GLOF风险减缓协同的潜力。敏感性分析显示，流域尺度估算和前几位候选湖在水量收支参数及权重假设变化下总体稳定，但GCM（全球气候模式）离散度是数值最大的不确定性来源（变异系数6.16%）。

主要技术方法：研究使用Copernicus 30 m数字高程模型（COPDEM）进行地形分析和河网提取，采用基于Shreve河网分级的分水岭提取方法确定冰湖集水区；使用Open Global Glacier Model（OGGM v1.6.2）模拟冰川径流，并利用6个CMIP6全球气候模式（BCC-CSM2-MR、CAMS-CSM1-0、FGOALS-f3-L、INM-CM4-8、INM-CM5-0、MRI-ESM2-0）在4种SSP情景（SSP126、SSP245、SSP370、SSP585）下的输出作为气候驱动；采用经验径流系数法和固定湖面蒸发深度（300 mm yr^?1）重建年度湖泊水量收支；通过自动化的四步河网交汇点识别方法确定候选水电场址和水头；使用层次分析法式多准则评分框架评估环境（保护区、受威胁哺乳动物丰富度、土地覆盖）、技术（GLOF风险等级、水力梯度、道路可达性、径流调节能力）和经济（平准化度电成本LCOE、居民点距离、交通基础设施）适宜性，并采用K-means聚类进行优先级分类。样本冰湖清单来自Zheng等（2021）的波曲河流域冰湖数据库，冰川边界来自Randolph冰川目录（RGI v6），保护区数据来自WDPA，道路数据来自OpenStreetMap。

主要结果：按章节总结如下。（1）湖泊入流量与径流组成：模型重建表明，流域尺度上冰川径流在2025–2100年间始终是主导水源，SSP585情景下因增强的消融抵消了冰川退缩导致的径流减少，后期冰川径流仍维持较高水平；计入非冰川径流、湖面降水、蒸发和第一阶段排水后，水电可用总水量高于单纯冰川径流。（2）冰湖水电潜力：流域尺度理论年均发电量在第一至第四阶段分别为212.8、200.6、197.0和194.3 GWhr yr^?1；高发电量湖泊集中在具有较大水头（平均245 m，最高776 m）和较大库容的湖泊，水头而非流量是控制潜力的主要地形因素。（3）开发适宜性：环境维度评分普遍较高，无湖泊与保护区重叠，67.6%的场址为裸地，受威胁哺乳动物平均2.9种/km²；技术维度障碍显著，97.3%的湖泊为冰碛湖，仅4个湖泊邻近可通车道路，45.9%的湖泊库容比>150%表明径流调节能力不足；经济维度上45.9%的湖泊LCOE低于0.1 USD kWhr^?1，但可达性差会推高成本。（4）综合分析：K-means聚类将37个湖泊分为9类中的8类，帕库普科是唯一高潜力高适宜性湖泊；5个高潜力中等适宜性湖泊（查玛曲丹错、嘎龙错、查热错、L12等）因可达性差或高风险评分受限；岗西错和龙决错属于中等潜力高适宜性类别。（5）不确定性与敏感性：GCM离散度（6.16%）、SSP差异（2.23%）、时间变率（2.11%）和DEM水头误差（2.02%）分别报告；水电效率扰动产生最大响应（±6.25%），蒸发假设变化（200–1000 mm yr^?1）使发电量变化-0.41%至-1.94%；权重替代方案下Spearman秩相关为0.95–0.98，前10名中9–10个保持稳定。

讨论与结论：讨论部分指出，有控的湖泊水位调节已在秘鲁布兰卡山脉等地区用于GLOF风险减缓，第一阶段20%–30%初始库容的排水分配对应37个湖泊的中位降深5.70 m（均值7.16 m），与尼泊尔伊姆贾错（Imja Tsho）、措尔罗帕（Tsho Rolpa）、不丹托托米湖（Thorthormi）及秘鲁513号冰川湖等已实施干预的降深量级相当。但针对具体湖泊的安全降深需考虑冰碛物组成、湖盆形态、出口几何和降深速率。替代开发模式如抽水蓄能（PSH）和多水库梯级在特定湖群（如岗西错-嘎龙错-贡错组合）可能提升运行灵活性，但需进一步研究。不确定性方面，非冰川径流系数和蒸发深度未独立率定且设为常数，未考虑未来植被演替和湖温变化；年尺度水量平衡适用于区域筛选，但56.8%的湖泊储流入库比<0.5年，需月尺度或更高分辨率模拟。环境方面，冰湖生态敏感性虽低但湖水位调节可能影响底栖无脊椎动物、硅藻、微生物和滨岸栖息地，下游河道冷脉冲和浊度影响需纳入未来评估。下游水文容量方面，第一阶段排水总量约321.1 Mm³，均匀释放相当于额外平均流量1.70 m³ s^?1，约占波曲干流年均流量（31.7 m³ s^?1）的5.4%，但需避免在季风洪水期排放。研究结论：该流域尺度筛选评估确定了波曲河流域37个冰湖的理论水电潜力与开发适宜性，帕库普科为最佳场址；四个阶段框架将风险减缓与适应性发电相结合，可在适当工程措施配合下支持当地清洁能源供应并降低下游洪水风险，为脆弱山地地区联合评估清洁能源潜力与冰湖灾害管理提供了途径。