重建青藏高原东南部澜沧江（湄公河上游）大峡谷全新世高幅度溃决洪水：对灾害评估的启示

《Journal of Hydrology: Regional Studies》：Reconstructing Holocene high-magnitude outburst floods in the Lancang River (Upper Mekong) Grand Cayon, Southeastern Tibetan Plateau: Implication for hazard assessment

　　

研究区：澜沧江（湄公河上游）大峡谷，位于青藏高原东南部，是三江并流区的一个高山峡谷。研究焦点：该研究整合了澜沧江大峡谷溃决洪水事件的沉积学、地貌学、年代学和水力学（hydraulic）证据。研究人员的调查识别出三个主要的古洪水记录：细粒滞水沉积物（fine-grained slackwater deposits）、砾石坝沉积物（gravel bar deposits，GBDs）和漂砾沉积物（boulder deposits）。这些具有相似地球化学组成的洪水沉积物表明同一洪水来源。上游两个滑坡坝引发的溃决洪水被定年为11.9–6.4 ka（千年）。区域新水文认知：全新世溃决洪水规模被估算为39,500–45,900 m³ s^-1（立方米每秒），超过最大观测洪水的7.5–10倍。这样的流量是由部分滑坡坝溃决产生的，侵蚀了坝高的50%以上。水动力变化和河道形态导致洪水期间的侧向加积和回水-涡流效应。基岩河道侵蚀优先发生在狭窄河曲弯道内的局部地形特征处。高幅度溃决洪水产生高剪切应力（0.4–1.1 kPa，千帕）和单位水流功率（4.5–10.5 kW m^-2，千瓦每平方米），并在澜沧江上游的地貌演化中起关键作用。基于古洪水最坏情景的深度-速度（depth-velocity，DV）指标表明，桥梁面临最高危害（DV > 240 m² s^-1），而道路/村庄则面临回水淹没（DV < 65 m² s^-1）。该研究为评估青藏高原东南部类似的级联洪水灾害提供了宝贵的见解。

论文解读文章

研究背景：青藏高原东南部的澜沧江（湄公河上游）大峡谷位于三江并流区，是一个典型的高山峡谷。高幅度溃决洪水由滑坡坝或冰碛坝溃决触发，其流量远超正常季风洪水，对河流地貌演化和下游基础设施构成重大威胁。然而，当前主要存在两个问题：一是短期水文观测记录（如年洪峰值）难以捕捉百年至千年尺度的极端洪水事件，导致工程防洪设计可能低估实际风险；二是该区域滑坡坝溃决洪水的水文地貌效应（如侵蚀、搬运和沉积过程）尚未得到系统研究。因此，研究人员开展此项研究，旨在利用古洪水沉积记录重建全新世高幅度溃决洪水的规模、时序和水力特征，填补澜沧江大峡谷古洪水知识的空白，并为类似级联洪水灾害评估提供科学依据。该论文发表在《Journal of Hydrology: Regional Studies》。

研究人员开展了什么研究：研究人员通过对澜沧江大峡谷（Quzika和Rengda两个地点）的古洪水沉积物进行系统调查，整合沉积学、地貌学、年代学和水力学证据。首先，他们识别并描述了三种主要古洪水沉积类型：滞水沉积物（SWDs）、砾石坝沉积物（GBDs）和漂砾沉积物。其次，利用光释光（OSL）和加速器质谱¹⁴C（AMS ¹⁴C）测年建立洪水事件年代序列。第三，通过二维水力学模型（HEC-RAS）和经验公式（如古流速-粒径关系）重建洪水的水位、流速、剪切应力和单位水流功率等参数。最后，基于深度-速度（DV）指标评估潜在灾害。主要结论：全新世（11.9–6.4 ka）期间发生了两次由上游滑坡坝部分溃决（侵蚀坝高超过50%）引发的溃决洪水，峰值流量达39,500–45,900 m³ s^-1，为最大观测季风洪水（2017年20年一遇洪水）的7.5–10倍。这些高能量洪水在狭窄弯曲段产生高剪切应力（0.4–1.1 kPa）和单位水流功率（4.5–10.5 kW m^-2），主导了澜沧江上游的地貌演化（如基岩侵蚀和砾石坝沉积）。基于DV指标的危害评估表明，位于河道中心线的桥梁面临最高危害（DV > 240 m² s^-1），而道路和村庄主要受回水淹没影响（DV < 65 m² s^-1）。该研究的重要意义在于为青藏高原东南部高余地形区的级联洪水灾害评估提供了长期基准和定量阈值。

关键的技术方法（不超过250字）：主要关键技术方法包括：（1）古洪水地貌调查与沉积学分析，利用无人机（UAV）影像和激光粒度分析仪（Camsizer P4及Malvern Mastersizer-S）测定GBDs和SWDs的粒度和形状参数；（2）地球化学分析，通过波长色散X射线荧光光谱仪（XRF）测量主要和微量元素含量，以上陆壳（UCC）标准化判断沉积物源；（3）年代学，采用光释光（OSL）单颗粒再生剂量法（SAR）和AMS ¹⁴C测年（由Beta实验室完成）建立洪水事件时序；（4）水力学重建，使用二维HEC-RAS模型（基于12.5 m分辨率Tan DEM-X地形数据）模拟洪水过程，并结合古流速-粒径经验公式（Costa 1983和O'Connor 1993）验证；（5）灾害评估，通过深度-速度（DV）指标（DV = 水深×流速）量化危害等级。样本来源为澜沧江大峡谷Quzika（QZK）和Rengda（RD）两个地点。

研究结果：
4.1 古代溃决洪水事件的地貌证据
4.1.1 Quzika（QZK）地点的古洪水GBDs和漂砾：通过无人机影像和野外调查，在QZK地点右岸识别出厚4–6 m的涡流砾石坝沉积（GBDs），含粗砂、砾石及直径1.2–3.5 m的漂砾。这些沉积物高于2017年洪水高水位标志（HWMs）8–11 m，表明古洪水规模远超现代季风洪水。沉积特征显示块状结构、向上变细序列以及漂砾上的光滑侵蚀表面，指示高能量溃决洪水搬运。
4.1.2 Rengda（RD）地点的古洪水GBDs和SWDs：在RD地点发现两种关键证据：滞水沉积物（SWDs）覆盖在侧向砾石坝沉积（GBDs）之上。侧向GBDs显示韵律层理和倾斜/水平层理，反映水动力变化；SWDs分为两类：SWD-1（松散粉砂质砂）和SWD-2（黏土质粉砂，含根孔和棱块状结构）。这些沉积物位置高于现代河床30–50 m，表明洪水水深极大。
4.2 古洪水GBDs和SWDs的粒度分布与地球化学：粒度分析显示：QZK的涡流GBDs砾石含量>80%，RD的侧向GBDs砾石含量>95%，且侧向GBDs粒度更粗。颗粒形状分析表明70%的GBDs长宽比>0.6，80%球度>0.7，指示高能量水动力搬运。SWDs粒度以粉砂为主（峰值约100 μm），细于GBDs。地球化学元素（经上陆壳UCC标准化）显示QZK和RD两地点GBDs与SWDs具有相似的配分模式，但QZK的GBDs因长期风化在Na₂O、CaO等元素上略有差异，证实两处沉积物源自上游同一溃决洪水事件。
4.3 古洪水事件的年代学：通过OSL和AMS ¹⁴C测年建立时序。RD地点底部侧向GBDs的OSL年龄为11.9 ± 2.1 ka，QZK地点涡流GBDs的OSL年龄为8.1 ± 1.5 ka至6.4 ± 1.1 ka。RD地点SWD-2的OSL年龄为8.6–7.3 ka，AMS ¹⁴C年龄为7.7–7.5 cal ka BP。综合结果确定全新世溃决洪水事件发生在11.9–6.4 ka（含年龄误差）。
4.4 全新世古洪水的水文重建：基于SWDs和GBDs顶面高程，利用HEC-RAS二维模型模拟得出峰值流量为39,500–45,900 m³ s^-1，最大水深51 m，最大流速16–17 m s^-1。漂砾古流速估算为5.9–6.3 m s^-1。高剪切应力（>0.25 kPa）和单位水流功率（>4.7 kW m^-2）出现在狭窄河道（QZK地点），SWDs沉积于低流速区（<3 m s^-1）。通过经验公式，推断更大滑坡坝（高130 m）的部分溃决（溃口深度约65 m，超坝高50%）可产生约32,850 m³ s^-1的峰值，而实际重建流量（39,500–45,900 m³ s^-1）更大，表明溃口深度超坝高50%。
4.5 高幅度溃决洪水的潜在灾害评估：基于古洪水最坏情景的DV指标（水深×流速）量化危害。模拟表明，所有位于河道中心线的桥梁面临最高危害（DV > 240 m² s^-1），道路和村庄受回水淹没（DV < 65 m² s^-1），部分低洼建筑被淹没，河岸易被高DV（130–200 m² s^-1）冲毁。古洪水重建为类似低频率高幅度洪水事件的减灾策略提供了关键依据。

总结讨论部分：讨论指出，高幅度洪水的水力学特性（高流速、高剪切应力、高单位水流功率）是形成不同高程洪水沉积物（SWDs、GBDs、漂砾）的主控因素。季风洪水无法搬运这些沉积物。在窄深峡谷向宽谷过渡段，流速骤降导致巨量粗颗粒沉积；在河曲弯道内，局部地形引发高剪切应力（最大约13 kPa）和单位水流功率（约350 kW m^-2），大于2018年白格溃决洪水（6.6 kPa, 173 kW m^-2）。与2017年20年一遇季风洪水对比，溃决洪水的水深为前者的2–9倍，流速为1.5–2.2倍，剪切应力为2–2.2倍，单位水流功率为3–5倍，表明高幅度溃决洪水主导了澜沧江大峡谷的地貌演化（如形成年轻阶地T1，约9–5 ka）。完全或部分滑坡坝溃决的相似情景在高原高余地形区可能重现，古洪水提供的长期基准（DV阈值）可指导桥梁（>240 m² s^-1）、河岸（>130 m² s^-1）和低洼定居点（<65 m² s^-1）的防灾规划。

研究结论翻译：对澜沧江（湄公河上游）大峡谷沿线的全新世高幅度洪水沉积物和现代洪水高水位标志（2017年洪水事件）进行了详细调查。从宏观到微观尺度，识别出不同高程的涡流/侧向GBDs、SWDs和漂砾的沉积特征。涡流GBDs含>80%砾石，侧向GBDs砾石含量超95%。两种SWD单元与相关GBDs的相似地球化学组成表明它们由同一溃决洪水事件沉积。年代约束的地层序列揭示全新世溃决洪水阶段发生在11.9–6.4 ka。以2018年白格滑坡溃决情景为基准的古水文重建表明，三江并流区内最大洪水水深达51 m，流速达17 m s^-1。全新世溃决洪水流量估计为39,500–45,900 m³ s^-1，超过当前最大观测季风洪水7.5–10倍。该规模由上游两个滑坡坝溃决推断，回溯模拟显示溃口深度超过坝高约50%。古洪水年龄与低级阶地年代同步表明高幅度洪水在全新世河流地貌演化中的重要作用。基于模拟古洪水范围通过DV指标量化，所有位于河道中心线的桥梁面临最高危害（>240 m² s^-1）。这些发现利用高幅度洪水的水力特性和地貌影响，为青藏高原东南部高余地形区提供了长期灾害评估依据。