龙门山前河流迁移的历史研究

《Water》：Historical Study on the Migration of Rivers in Front of Longmen Mountains

　　

摘要：冲积平原上的河流迁移反映了地形、水文过程和泥沙输运之间的复杂相互作用。然而，在复杂地貌背景下，这些因素在控制相互作用的河流系统长期迁移中的相对作用仍未被充分理解。本研究利用基于弗劳德相似律（Froude similarity）的几何比尺为1:1200的动床物理模型，研究了龙门山前河流的长期迁移特征。研究人员再现了岷江、湔江和石亭江在不同水沙条件下的演变过程，并应用三维（3D）激光扫描技术量化了岸线迁移和侵蚀-淤积格局。结果表明，不同河流系统表现出不同的迁移响应，在72小时时不同水流条件下受影响岸线面积为0.93至3.11平方米。河流迁移受地形、河道曲率和泥沙供给的综合控制。具体而言，湔江表现出与丰富泥沙供给、地形坡度骤减和不对称岸滩地形相关的显著迁移变异性。通过将物理模型观测与历史和考古证据相结合，本研究为河流迁移、景观演化与三星堆遗址相关环境条件之间的相互作用提供了见解。

论文解读：龙门山前河流迁移的物理模型研究

一、研究背景与科学问题

冲积平原上的河流迁移是地形、水文过程和泥沙输运共同作用的复杂地表过程。河流通过侵蚀和淤积不断重塑河谷地貌，形成阶地、泛滥平原等典型冲积地貌单元。然而，在复杂地貌背景下，特别是多个河流系统相互作用的区域，地形梯度、泥沙积累和河流间相互作用如何共同调节长期河道演变，以往研究仍缺乏系统认识。龙门山前冲积平原——成都平原，由源自龙门山的多条河流长期泥沙堆积和地貌相互作用形成，具有边缘高、中心低的扇形冲积平原地貌特征。其中，湔江在出山进入成都平原后流向发生显著改变，历史记录显示其经历了多次改道和分汊。理解这些变化的控制机制，对于揭示地形、河道曲率和泥沙输运如何共同控制山前冲积平原河流迁移具有重要意义。已有成都平原地质和地貌研究主要基于沉积记录和地貌重建分析，侧重于结果地貌演化，而对相互作用河流动态、泥沙再分配与河道迁移的过程机制缺乏约束。因此，本研究通过构建包含岷江、湔江和石亭江水系的大比尺动床物理模型，探讨三个核心问题：控制龙门山河流迁移的主导因素是什么？长期河流迁移和泥沙再分配如何塑造成都平原？河流迁移及水文变化如何改变三星堆遗址环境条件？

二、研究技术与方法

研究人员基于弗劳德相似律设计了几何比尺为1:1200的动床物理模型，模拟了石亭江、湔江、岷江出山口河段及其周边的典型地貌特征，模型中还纳入了湔江支流鸭子河。模型平台长8米、宽4米，采用中值粒径d₅₀ = 1毫米的均匀细沙作为床沙。实验以2020年长江流域第5号洪水为参考水文条件，根据三河多年流量比例确定了各河初始流量：湔江300 mL/s（输沙率80 mg/s）、岷江800 mL/s（输沙率270 mg/s）、石亭江300 mL/s（输沙率30 mg/s）。研究人员还开展了敏感性试验，仅将湔江流量和输沙率同步增加20%。实验过程中采用RIEGL VZ-400i地面三维激光扫描仪周期性获取高精度点云数据，利用Riscan Pro和Geomagic Qualify软件进行点云处理、三维模型重建和侵蚀-淤积定量分析。

三、主要研究结果

3.1 试验过程与河道演变

在初始水流条件下，三条河流均表现出不同程度的侧向迁移和河岸调整，高曲率河段河岸变化尤为显著。在增加水流条件下，虽然岷江自身流量和泥沙供给未改变，但其河道演变与初始条件相比出现差异，表明一个河流系统的变化可能影响相连河网中相邻河道的形态响应。

3.2 湔江迁移特征

初始条件下，湔江出口河段最大局部位移在24小时达36.4厘米，48小时入口附近侧向位移达61.0厘米，72小时下游段位移达92.4厘米。增加流量条件下，24小时入口左岸最大位移达67.6厘米，48小时达89.2厘米，96小时水流发生显著偏转，主流逐渐从湔江主河道转向鸭子河支流，最大侧向位移增至156.4厘米。

3.3 岷江迁移特征

初始条件下，岷江在24小时形成弧形河道和扇状三角洲，最大位移118.5厘米；48小时主流转向左岸，位移增至131.3厘米；72小时左岸侵蚀加剧，位移达168.9厘米。增加流量条件下，24小时位移达153.9厘米，48小时入口河道呈“沙漏”形，96小时位移增至225.9厘米。

3.4 石亭江迁移特征

初始条件下，石亭江24小时最大位移87.7厘米，48小时增至144.3厘米，72小时达159.5厘米。增加流量条件下，24小时位移仅42.7厘米，48小时为44.5厘米，96小时受鸭子河影响下游段形成新弯道，位移增至77.9厘米。

3.5 岸线迁移定量分析

72小时时，初始条件下岷江、湔江、石亭江高曲率侧受影响岸线面积分别为3.11平方米、1.95平方米和2.31平方米；增加流量条件下分别为2.96平方米、1.42平方米和0.93平方米。按几何比尺换算后对应初始条件下4.48平方公里、2.81平方公里和3.33平方公里，增加流量条件下4.26平方公里、2.04平方公里和1.34平方公里。相对变化分别为-4.8%、-27.2%和-59.7%，表明增加水力强迫并未导致河岸调整的均匀响应。河道演变经历三个阶段：0-3小时河道展宽阶段、3-72小时不稳定调整阶段、72-96小时趋势稳定阶段。历史文献记载的湔江多次分汊和改道与模型观察到的往复岸线调整和汊道发育具有一致性。

3.6 侵蚀-淤积格局演变

72小时时，初始条件下岷江、湔江、石亭江侵蚀面积分别为4.14平方米、3.26平方米和3.65平方米，淤积面积分别为0.09平方米、0.36平方米和0平方米；增加流量条件下侵蚀面积分别为3.49平方米、2.44平方米和1.58平方米，淤积面积分别为1.66平方米、0.75平方米和0平方米。与初始条件相比，增加流量条件下侵蚀面积分别减少了约15.7%、25.2%和56.7%，但淤积面积增加，尤其是岷江和湔江。

四、讨论与结论

讨论部分重点分析了地形因素、河道形态和泥沙输运对河流迁移的控制作用。地形方面，湔江两岸较高而河道间区域较低，形成“高-低-高”地形格局，主流向低海拔区域迁移，横向上地形梯度影响了湔江迁移的空间格局。河道曲率方面，岷江和石亭江在高曲率河段表现出更明显的迁移，符合曲率驱动河岸调整的既有认识。泥沙供给方面，湔江虽然流量最小，但其较高单位泥沙供给导致显著河道调整，表明泥沙可用性强烈影响迁移强度。对于成都平原的形成，模型试验表明不同河流间泥沙输运、侵蚀-淤积差异可产生地貌调整的空间差异，为理解边缘高、中心低的扇形平原发育提供了过程机制解释。对于三星堆遗址，物理模型显示湔江流量增加会显著增强向鸭子河的分流，为极端洪水事件可能加剧遗址区域洪涝风险提供了水力学机制。结合遗址内厚约20-50厘米的洪水淤积层和考古证据，频繁的河道变迁和洪水事件对三星堆文明产生了深远影响，但外部势力入侵和政治中心转移也是其衰亡的重要因素。研究结论可归纳为：河流迁移受地形、河道曲率和泥沙输运综合控制；湔江偏转迁移受高泥沙供给、坡度骤减和不对称岸滩地形共同作用，增加流量主要改变泥沙再分配格局而非简单加强侵蚀；成都平原通过长期河流沉积和多河流相互作用形成；河流迁移和水文变化可能影响了三星堆文明的环境条件和洪涝风险，但其衰落受自然和社会多重因素共同影响。