一种预测黄河河道冲刷的方法
【技术领域】：
[0001] 本发明涉及一种水利工程汛期险情预测方法，特别涉及一种洪水期间对黄河河 道是否发生冲刷进行预测，为防汛抢险提供预警手段的方法，其IPC国际专利分类号E02B 1/00,属于水利领域。
【背景技术】：
[0002] 黄河洪水含沙量大，而且洪水的流速越大，挟沙力也越大，对河床的冲刷力更强， 往往会造成黄河主河道的冲刷下切。黄河下游的险工工程靠近主河道，河道的冲刷下切可 能会导致险工的坝基被主流严重淘刷，导致根石流失、坝体墩蛰入水的重大险情。除了每年 的小浪底调水调沙外，黄河没有发生过较大洪水，险工坝前的河道缺少被冲刷的机会，细颗 粒沉积增加，河道表面逐步形成块状的固结沉积物，在河床上构成了一层沉积块保护层。但 是，实际观测及室内试验证明，这层保护层并不是不可冲刷的，在某些来水条件下这层沉积 块会被揭掉冲走，进而导致河床被冲刷下切。研究沉积块被揭起冲刷的影响因素，并确定沉 积块揭起的临界条件，在各级洪水时发生险情的可能性进行预测，具有重要的现实意义。
【发明内容】
：
[0003] 本发明的目的在于提供一种能够确定影响沉积块揭起的影响因素的数学模型，并 并结合实测数据，预测黄河河道冲刷的方法。
[0004] 为实现上述目的，本发明采用如下技术方案。
[0005] -种预测黄河河道冲刷的方法，包括以下步骤：
[0006] (一）建立水流作用下河道沉积块起动流速的数学模型：
[0007] F_+Ff-G-Fs= 0 (1)
[0012] 式中：F_为水流作用下沉积块起动时的最大脉动上举力；Ff为沉积块受到的浮 力；Fs为沉积块受到的剪断力；G为沉积块重力；Y"为浑水容重；d、h分别为沉积块的直径 与高；ys为沉积块的湿容重；J为河道比降；K为河道沉积块起动临界系数；vs为沉积块起 动时的水流速度；g为重力加速度；C为沉积块土体粘聚力；巾为沉积块土体内摩擦角。
[0013] (二）选取需要预测的险工工程所在河段作为试验河段，在当年的黄河调水调沙 前，在试验河段的河道布置视频监控设备和流速测定仪；
[0014] (三）通过试验或现场实测获得如下参数：河道沉积块的形状参数、沉积块的湿容 重ys、沉积块的强度参数、浑水容重^和河道比降J;方法如下：
[0015] 1.河道沉积块的形状参数：沉积块的高度h和直径d
[0016] 根据河道沉积块的高度和面积，将每个河段的沉积块分为多个类别，选取其中数 量最多的一类沉积块作为试验沉积块；测量从现场获得的试验沉积块，得到所有块体的高 度，采用取平均值或中值的方法确定试验沉积块的高度h;测量每个试验沉积块的周长L， 并采用等效周长法确定每个沉积块的直径，然后采用取平均值或中值的方法确定试验沉积 块直径d;
[0017] 2.进行室内土工试验确定沉积块的湿容重Ys;
[0018] 3.进行室内土工试验确定河道沉积块强度参数：土体粘聚力c和土体内摩擦角 伞；
[0019] 4.根据水情预报获得浑水容重Y"，或者提取发生沉积块起动时的水样，测量水样 的浑水容重Yn;
[0020] 5.通过现场实测或参考河流地形资料确定河道比降J;
[0021] (四）实测河道内试验沉积块的实际起动流速vs。
[0022] 利用视频监控设备和流速测定仪，测定所述试验沉积块实际起动流速'。，测定 3~5次，计算其平均值巧〇 ;
[0023](五）确定河道沉积块起动临界系数K
[0024] 将步骤（三）得到的各项参数和步骤（四）得到的实际起动流速平均值代入 步骤（一）的数学模型中，联立求解，得到河道沉积块起动临界系数K;
[0025] (六）将步骤（五）得到河道沉积块起动临界系数K代入步骤（一）的数学模型 中，联立求解，即得到该试验河段或类似河段的河道沉积块起动临界流速测算公式：
[0026]
(15)
[0027] (七）每次洪水来临时，根据水情预报得到洪水的浑水容重y将步骤（三）获 得的沉积块的高度h和直径d、湿容重ys、沉积块土体内摩擦角巾及土体粘聚力c和河道 比降J代入公式（15)，即可计算得到本次洪水时沉积块起动临界流速vs;根据水情预报，若 本次洪水的预计流速大于该临界流速vs，则河道会发生冲刷，应当做好相应的抢险准备；若 本次洪水的预计流速小于该临界流速VS，则河道不会发生冲刷。
[0028] 本发明提供的预测黄河河道冲刷的方法，在建立数学模型时，引入了代表沉积块 本身物理特性的参数，并结合待预测河段的试验及实测数据计算河道沉积块起动临界流 速，能够更真实地反映沉积块起动时的受力情况，准确预测每次洪水时河道是否会发生冲 刷，并做好相应当抢险准备。
【具体实施方式】：
[0029] 下面对本发明做详细说明。
[0030] 本发明的预测黄河河道冲刷的方法，包括以下步骤：
[0031]( -）首先，建立本发明所采用的数学模型，同时也阐明本发明的原理。黄河河 道天然形成的沉积块，多为黏土形成的胶泥块，为了便于理论分析，假设沉积块为一薄圆柱 形，根据静力平衡原理可知，河道沉积块起动时的临界受力情况如下式所示：
[0032]F_+Ff-G-Fs= 0 (1)
[0033] 其中：G为沉积块重力；FS为沉积块受到的剪断力；Ff为沉积块受到的浮力为 高含沙水流作用下沉积块起动时的上举力。其中：
[0034]Ff=ynV (2)
[0039] 式中：yn为浑水容重，与水的含沙量有关；d、h分别为沉积块的直径与高；Y3为 沉积块的湿容重；v为沉积块体积；j为河道比降；a为沉积块表面积，;K为河道 沉积块起动临界系数，由试验率定；vs为沉积块起动时的水流速度；g为重力加速度。
[0040] 发明人根据大量的试验数据和实际观测发现，由于各段河道的沉积块的化学成 分、沉积时间、受力历史等不同，导致沉积块的土体粘聚力、土体内摩擦角等物理性质的不 同，这些物理性质，也是影响沉积块起动流速的关键因素。因此，本发明建立了公式（7)、 (8)：
[0041] tAxz= 〇nAxy (7)
[0042] X=c+ 〇ntan<i> (8)
[0043] 式中：t为抗剪强度；o。为沉积块所受正应力；Axz为沉积块x、z轴向面积，Axz = Jrdh;Axy为沉积块x、y轴向面积，S卩表面积，
c为沉积块土体粘聚力；小为 沉积块土体内摩擦角。
[0044] 贝丨J
[0053] 如通过多次实测得到某段河道沉积块的平均起动流速为&〇,则可得到该段河道 沉积块起动临界系数K，从而可以率定该段河道沉积块起动临界系数K。
[0054]
[0055] 由此可见，沉积块起动流速vs，即河道发生冲刷的临界条件，与河道沉积块的形状 (高度h和直径d)、湿容重ys、沉积块的抗剪强度（内摩擦角巾和粘聚力c)、含沙量（由 浑水容重丫"体现出来）、河道比降J等有关。这些参数与土的力学特性指标及河床形态密 切相关，均可根据现场实测和室内试验取得。
[0056] 利用上述数学模型，本发明进一步采取以下步骤：
[0057] (二）自2002年以来，每年的6月份，在黄河水利委员会的同一指挥、调度下，进 行一次调水调沙。黄河调水调沙，就是利用小浪底水库，或者对小浪底水库及上游的多个水 库进行联合