[发明专利]冲积河道崩岸过程的数值计算方法

说明书

技术领域

本发明涉及水利工程技术领域，具体的说是一种冲积河道崩岸过程的数值计算方法。

背景技术

崩岸是冲积河流河床变形过程的一个重要方面。大规模崩岸不仅会毁坏已建的河道整治建筑物、增加沿岸城市的防洪压力、影响岸边涉水设施的正常运行，而且可能造成主流及深槽摆动，破坏局部河段的河势稳定。然而，现有的崩岸模拟及监测技术尚不完善，不能及时预测或实时观测到崩岸的发生及发展过程。因此通过数学方法对崩岸过程进行模拟，成为研究河道崩岸规律及开展崩岸预测的一种重要手段。

现有崩岸模拟方法通常基于河流动力学或土力学中边坡稳定理论。

前者一般采用泥沙运动力学理论，结合实测数据分析，确定影响崩岸的关键水沙因子，并提出河岸崩退宽度、速率等与这些水沙因子之间的经验关系，如平滩河宽与来流流量之间的幂函数关系，河岸崩退速率与剩余水流切应力之间的幂函数关系等。这种半经验半理论方法通常能反映长时段的河道崩岸，但由于崩岸机理的复杂性及河床形态与河岸土体特性等方面的差异，它不适用于预测短时段内河道的崩岸过程，也无法确定崩岸发生的时刻及规模等。

从土力学角度出发的研究，则基于河岸崩塌的不同力学机理，提出了各类崩岸模式及相应的河岸稳定性计算方法，并部分考虑不同的冲积及非冲积因素对河岸稳定性的影响，如水流冲刷、河道水位升降及河岸土体特性变化等因素。这类方法能计算出一段时间内河岸的崩退过程，包括河岸形态的变化、崩退宽度及次数等，但由于崩岸机理涉及因素较多，要考虑各因素的综合效应目前仍具有一定难度。

发明内容

本发明的目的在于提供一种冲积河道崩岸过程的数值计算方法，能及时预测某个特定断面的二元结构河岸在一定水沙过程作用下稳定性及岸坡形态的变化过程，如崩岸发生时刻、崩塌次数及崩退宽度等，为崩岸预测及预警控制提供科学依据和技术支撑。

本发明的原理为：河岸崩退与近岸水流冲刷、河岸土体特性、河道水位变化及其引起的潜水位变化等因素密切相关，有必要综合考虑这些因素的影响。近岸水流冲刷会导致河岸变陡及其高度增加，这通常是天然冲积河道崩岸发生的主要控制因素，如长江中下游80％的崩岸均发生于迎流顶冲或水流剧烈冲刷的凹岸弯顶位置；而河道水位升降及其引起的潜水位变化会改变土体的受力条件及力学特性，尤其对退水期的崩岸发生具有明显影响。因此，本发明结合坡脚冲刷过程、潜水位变化过程与河岸稳定性三个方面的计算，提出一种二元结构河岸崩退过程的数值计算方法。

本发明冲积河道崩岸过程的数值计算方法，包括如下步骤：

步骤1，根据上年汛后的实测断面地形及当年的实测或计算水位过程，确定河岸坡顶及坡脚，计算河岸高度、初始坡度等河岸形态参数；

步骤2，计算河岸坡脚冲刷宽度及床面冲淤厚度，具体步骤如下：

步骤2.1，计算一定时段△t内的床面冲淤厚度；

步骤2.2，计算水流对坡脚的横向冲刷宽度；

步骤2.3，基于步骤2.1和2.2，修改河岸形态；

步骤3，计算河岸土体内潜水位变化过程，具体步骤如下：

步骤3.1，在河岸土体内划分计算节点，确定初始起潜水位；

步骤3.2，根据一维渗流控制方程，计算各节点在一定时段后的潜水位；

步骤3.3，基于步骤3.1-3.2，确定岸坡一定范围内潜水位的平均值；

步骤4，计算河岸稳定安全系数，判断其是否发生崩塌，具体步骤如下：

步骤4.1，基于步骤2，确定土体内部破坏面位置；

步骤4.2，计算潜在的滑动土体的稳定性，即河岸稳定安全系数F_s；

步骤4.3，对比F_s与临界值的大小，判断河岸是否发生崩塌；若是，则根据破坏面修改河岸形态，并进入步骤4.4；若否，进入步骤4.5；

步骤4.4，河岸发生崩塌后，确定坡脚堆积体的形态；

步骤4.5，当前时段的河岸稳定计算完成后，进入下一时段的计算，返回步骤2；

步骤5，划分崩岸强度，具体步骤如下：

步骤5.1，记录整个计算时段内河岸的崩塌次数N_T和崩退宽度△W，确定年均崩退宽度△w；

步骤5.2，依据年均崩退△w值，通过设定的崩岸严重程度等级划分体系，确定崩岸强度等级。

本发明一种冲积河道崩岸过程的数值计算方法的有益效果在于：

(1)考虑了多种崩岸影响因素的共同作用，能同时模拟河岸崩退过程、水流对河岸坡脚的冲刷过程及河岸土体内潜水位的变化过程；