高含沙洪水“揭河底”冲刷临界条件综合判别方法技术

本发明专利技术涉及一种水文预报方法，特别涉及一种高含沙洪水“揭河底”冲刷临界水流条件综合判别方法，属于泥沙动力学技术领域，IPC国际专利分类号为E02B1/00。该方法根据“揭河底”发生机理，利用胶泥块抗折强度试验及抗剪强度试验数据，构建了四种胶泥块可能揭掀的力学分析模型，并考虑胶泥块揭掀综合影响因素及最不利荷载组合条件，分别计算胶泥块揭掀时的临界水流流速，最后取计算结果的最小值作为“揭河底”现象发生时的临界水流条件。本发明专利技术的判别方法，能够较准确预报某次洪水是否能导致“揭河底”现象发生，为防汛抢险的调度、准备工作提供了决策依据，同时，也为“揭河底”现象发时的实时数据采集提供了可能。

【技术实现步骤摘要】

【技术保护点】
一种高含沙洪水“揭河底”冲刷临界条件综合判别方法，包括以下步骤：(一)、“揭河底”冲刷易发河段胶泥块抗折强度试验：试验土样为取自“揭河底”冲刷易发河段的胶泥块，首先对试验土样进行重塑，具体方法是：将从“揭河底”冲刷易发河段原型河床取回的胶泥块碾碎，放入容器中加水搅拌，使泥沙絮凝沉积，通过风干形成胶泥块试样，并利用风干时间不同控制胶泥块试样的含水率；对不同组次的胶泥块试样，利用抗折试验机测定抗折强度；试验方法如下：将胶泥块试样横放在抗折试验机的两根支撑圆柱上，试样长轴垂直于抗折试验机支撑圆柱，以50N/s±10N/s的速率通过抗折试验机的加荷圆柱将垂直荷载均匀地加在棱柱体的水平面上，直至胶泥块试样折断；利用式(1)计算出胶泥块的抗折强度M，M=1.5FfLB3---(1)式中，Ff?试件折断时垂直加在其顶部平面中部的荷载，单位为N；L?支撑圆柱之间的距离，单位为mm；B?棱柱体正方形截面的边长，单位为mm；以一组三个胶泥块试样抗折强度的平均值作为试验结果，利用实验结果建立胶泥块抗折强度与含水率的关系方程，由该关系方程计算出“揭河底”冲刷易发河段的胶泥块饱和含水率下的抗折强度Ms；(二)、“揭河底”冲刷易发河段胶泥块抗剪强度试验：采用与步骤(一)相同的胶泥块试样，对试样进行密封浸泡，使胶泥块试样达到饱和含水率，利用应变控制式直剪仪进行饱和快剪试验，得到胶泥块的抗剪强度τ’；试验方法：胶泥块试样分别在100kPa、200kPa、300kPa和400kPa四级垂直压力下，以0.8mm/min的剪切速率快速施加水平剪力使试样剪切破坏，试样加载时间为3～5min，当测力计读数不再增加或开始倒退时，记下破坏值；(三)、构建“揭河底”冲刷时胶泥块因力矩作用折断揭掀的受力数学模型：FD水流正面推力：FD=γmCD2Vbi2bcg---(2)式中，Vbi为胶泥块底部水流流速，i＝1，2，3，4，CD为阻力系数，γm为浑水容重，b为胶泥块宽度，c为胶泥块厚度，F’D为折断面上FD的反作用力：FD,=-γmCD2Vbi2bvg---(3)τ为胶泥块上表面的水流拖曳力：τ=γmCD2Vbi2lbg---(4)Cs胶泥块重量为：Gs＝γslbcg??????(5)式中，γs为胶泥块的容重，l为冲蚀悬臂长度，g为重力加速度，P1为胶泥块所受水压力：P1＝γmlbhg???????(6)式中，h为胶泥块表面处水深，P2为胶泥块的浮力：P2＝γmlb(h+c)g???????(7)最大水流脉动上举力：Fdmax=KγmJVbi2lb---(8)式中，K为线性系数，单位kg·m·s，取值范围：3～4.2，J为水面坡降，构建力矩平衡平衡方程：Fdmaxl+P2l2+τc-Gsl2-P1l2=Msbc---(9)将式(2)?(8)代入式(9)中并求解，得到式(10)：Vb12=2Msc+l2cg(γs-γm)2KJl2γm+lcgγm---(10)胶泥块在Δt时间段内被淘刷的冲蚀悬臂长度l采用公式(11)计算：l=C1Δt[0.16D2ρmVb12+AD(γc-γm)(1+ξ)]---(11)式中，C1为河床冲刷系数，γc为河床淤积物的容重，D为胶泥块周围可动淤积物的平均粒径，ρm为高含沙洪水的密度，A为综合影响系数，其取值范围为0.015～0.034，中值为0.0245，系数ξ的取值为0～1，将式(10)和(11)联立求解，即可求出Vb1；(四)、构建“揭河底”冲刷时作用在胶泥块水深方向的受力数学模型：作用在胶泥块水深方向的力平衡方程：Fdmax+P2?P1?Gs＝τ“bc?????????（12)式中，τ’为胶泥块的抗剪强度，由步骤(二)的抗剪强度试验给出，将式（5)?(8)代入式(12)，得到式(13)：Vb22=τ′c+lcg(γs-γm)KJlγm---(13)联立式(11)和式(13)，并将式(11)的水流流速用式(13)替换，即可求出Vb2；(五)、构建“揭河底”冲刷时因力矩作用胶泥整片揭掀受力数学模型：胶泥块瞬时揭掀力矩方程：KJVb32gc=ga2l2(γsγm-1)+πq0γm-1δ03δ12(al-1)2D1(D1+2t)2[δ12t2-1](b+2fc2b...

【技术特征摘要】

【专利技术属性】
技术研发人员：董文胜，江恩惠，刘晓文，王永军，翟丽平，刘雪梅，王卉，李军华，翟莹莹，
申请(专利权)人：华北水利水电大学，
类型：发明
国别省市：