[发明专利]一种河流险工冲刷安全分析与岸坡稳定性判定方法

摘要

本发明公开了一种河流险工冲刷安全分析与岸坡稳定性判定方法，步骤一：生成作为构建二维河道水/沙模型的数据集；步骤二、根据研究区域中河段参数，设置模型计算参数与定解条件，构建研究区域二维水/沙数值模拟模型；步骤三、根据不同洪水条件下的来沙情况设定多种工况，利用二维水/沙数值模型对研究区域河道水/沙输移及河床冲淤进行模拟，基于模拟结果对洪水前后险工河床进行分析，分别计算河流险工岸坡河床冲深距离和河侧向侵蚀距离；步骤四：构建险工岸坡渗流～稳定分析有限元模型，定量化研究各因素与岸坡稳定安全系数间关系。本发明研究了水流冲刷作用下的河道险工岸坡稳定性，为险工岸坡稳定性评估提供了定量化手段。

主权项

1.一种河流险工冲刷安全分析与岸坡稳定性判定方法，其特征在于，该方法包括以下流程：步骤一：依据河道河势变化和堤防历史出险情况，明确研究区域地形数据与险工工段几何尺寸，将水位、流量与含沙量数据作为边界条件进行该研究区域的水文数据与含沙量数据分析，生成作为构建二维河道水/沙数值模拟模型的数据集；步骤二：根据研究区域中河段的河流走向、河道断面、数字地形高程，险工几何尺寸，设置模型计算参数与定解条件，采用MIKE21中的水动力模块与泥沙模块构建研究区域二维水/沙数值模拟模型：水动力模块模型：连续性方程表达式如下：动量方程表达式如下：式中，t为时间；u为流速在x方向上的分量，η为河床底高程，d为静水深，h为水深，h＝d+η，g为重力加速度，f为科氏力；ρ为水的密度；s_xx,s_xy为辐射应力分量；p_a为大气压强，ρ₀为水的相对密度，S为源项，u_s,v_s为源项水流流速，为沿水深方向的流速平均值，T_ii为侧向应力；泥沙模块模型：悬移质非恒定不平衡输沙基本方程表达式如下：式中，h为水深，S为悬移质含沙量，t为时间，p_mod、q_mod分别为通量修正值在x,y方向上的分量，e_x，e_y分别为x,y方向的紊动扩散系数，ɑ为恢复饱和系数，与Rouse系数有关，ω_*为泥沙沉速，S_*为悬移质挟沙力。推移质输沙率表达式如下：式中，T为无量纲泥沙输移系数；D_*为无量纲泥沙粒径参数；s为泥沙相对密度；d₅₀为泥沙中值粒径。步骤三：根据不同洪水条件下的来沙情况设定多种工况，利用二维水/沙数值模型对研究区域河道水/沙输移及河床冲淤进行模拟，基于模拟结果对洪水前后险工河床进行分析，分别计算河流险工岸坡河床冲深距离和河侧向侵蚀距离；步骤四：构建险工岸坡渗流～稳定分析有限元模型进行模拟，利用模拟结果分析河道水位、河道涨退水速率、坡脚河床冲深距离和侧向侵蚀距离对岸坡稳定性的影响，分别选定不同的水位、涨退水速率、冲深距离和河侧向侵蚀距离输入渗流稳定模型，计算岸坡稳定安全系数K，得到各因素对岸坡稳定性的影响，对单因子相关关系进行分析，采用非线性拟合方法对各自变量与因变量K间关系进行拟合，提出适用研究区域的岸坡稳定安全系数K的预测经验公式如下：K＝αe^{(βΔZ+σΔB)}+λH^(ηV)    (1)其中，ΔZ为坡脚河床冲深、ΔB为侧蚀距离、H为相对水位、V为水位涨落速率，α、β、σ、λ、η为常数；利用模拟结果论证公式的可靠性；定量化研究各因素与岸坡稳定安全系数间关系，如下：当K＜1时，则岸坡稳定性较低，极有可能出现失稳破坏现象；当K≥1时，则岸坡稳定性较高，不易发生失稳现象。