摘要： 水库泥沙淤积在闸门前，会增大闸门启门力。如何通过模型试验确定有泥沙淤积时闸门启门力，本文对此进行了探讨。认为淤积在坝前的细颗粒泥沙多数情况下可当做宾汉体泥浆，根据坝前淤泥受力平衡条件，导出淤泥相似准则，并论述了模拟试验方法 。
关键词： 水库泥沙 闸门淤堵 启门力 泥浆 相似准则 模型试验
1 问题 的提出
水库建成之后，随着泥沙淤积的不断发展 ，最终将会淤积到坝前。为了排泄坝前淤泥及维持水库长期使用库容，往往需要设置底孔。在底孔关闭期间，门前会有泥沙淤积[1] 。淤积在门前的泥沙不仅增加了对闸门的水平压力，而且还增加了对闸门的附着力(粘结力)，导致闸门启门力增大。由于启门力是选择闸门启闭机容量的重要依据，所以泥沙淤积对闸门启门力的影响 是设计人员所关注的重要问题之一。
在闸门前有泥沙淤积情况下，如何确定启门力，目前 还几乎没有方法可循。工程设计中常常是在确定了清水启门力之后，然后考虑到门前泥沙淤积，再乘以一个大于1的系数，作为有泥沙淤积时的启门力，但该系数确定往往凭经验而定，带有相当的盲目性。文献 [2]针对三门峡水库具体情况，曾给出一个有泥沙淤积时启门力经验公式，但在应用 中也有一定的局限性。那么，解决该问题的另一途径就是模型试验。但由于模型试验中的一些关键技术问题至今未解决，如门前淤泥相似等，所以这方面的模型试验工作还未见诸报道，可以说还是一个空白。而工程设计又迫切需要模型试验研究 人员提供有泥沙淤积时闸门启门力，因而尽早开展这方面的研究工作非常必要。本文对门前淤泥相似准则及模拟试验方法进行了探讨，以便能够更好地开展这方面研究工作，满足工程设计需要。
2 门前淤泥相似准则
粘性细颗粒淤泥，随着其固结程度，含水量不同，其所处的物理状态也不同。当淤泥未固结，含水量较多时，属于宾汉体泥浆。随着淤泥固结，含水量变小，其力学性质发生根本变化，不再具有泥浆性质，而属于土力学所研究的范围。淹没在水下的坝前粘性细颗粒新鲜淤积物多数情况下为宾汉体泥浆。文献[3]认为，如果淤泥的干容重在810kg/m3 ～1080kg/m3 左右，那么就属于泥浆。表1统计了国内几座大型水库实测的坝前淤泥干容重[4] 。可见，这些水库坝前淤泥大部分为泥浆。
对于处在静止状态下的泥浆，作用的外力同时有重力、粘结力及压力，而粘结力起主导作用。由于影响粘结力的因素十分复杂，目前还难以用理论 关系式表达，但粘结力与切力之间存在如下关系[5]
式中：τ为泥浆切应力；P为压强；为泥沙内摩擦角；C为单位面积的粘结力。
表1 国内几座大型水库坝前淤泥干容重统计
坝前淤泥
工程名称
中值粒径/mm
干容重/(kg/m3)
刘家峡
0.020
900～1000
盐锅峡
0.0435～0.325
1030～1460
青铜峡
0.00607～0.094
1390
龚嘴
0.03～0.04
800～1000
丹江口
0.01～0.039
600～800
泥浆的内摩擦角主要与相对干容重γ′/γs 有关(其中γ′为泥浆干容重，γs 为沙粒容重).根据γ′/γs 可计算 出泥浆的内摩擦角[6，7] 。通过调整模型泥浆稠度，改变其干容重，可以使模型与原型泥浆内摩擦角相同[8] 。如果模型与原型泥浆内摩擦角相同，则由式(1)可推导出如下相似关系
式(2)、式(3)表明，在泥浆内摩擦角相同情况下，若泥浆切应力比尺能够满足相似要求，就能保证泥浆粘结力和压力相似。
为了导出泥浆切应力相似比尺关系式，列出泥浆静态平衡方程式如下[9]
(4)
式中：p、τ分别为法向应力(压强)和切应力，其下标第1个坐标表示应力分量所垂直的面，第2个坐标表示应力的方向，且有τxy =τyx ，τxz =τzx ，τyz =τzy ；ρ为泥浆密度；X、Y、Z分别为单位质量力在x、y、z方向的分量，当质量力仅是重力时，则X=Y=0,Z=-g.
若模型与原型泥浆相似，则模型与原型对应的物理量必成同一比例，即
(pxx )p /(pxx )m =(pyy )p /(pyy )m =(pzz )p /(pzz )m =λp
(τxy )p (τxy )m =(τxz )p /(τxz )m =(τyz )p /(τyz )m =λτ
xp /xm =yp /ym =zp /zm =λL
ρp /ρm =λρ ；gp /gm =λg
式中：脚标p表示原型；m表示模型。
将上述比例常数代入式(4)中，进行相似变换，得到
式中：λp 为泥浆压强比尺；λτ 为泥浆切应力比尺；λγm =λρ λg 为泥浆容重比尺；λi 为几何比尺。
宾汉体泥浆的切应力方程为
式中：τ为切应力；τb 为极限剪应力；η为刚度系数；dv/dn为流速梯度。τb 、η可通过流变试验确定。对于静态泥浆dv=0，则有
将式(9)代入式(6)中，得
满足式(10)就能保证泥浆粘结力和重力同时相似，而泥浆压力相似是粘结力和重力相似的必然结果。由式(5)可确定浆泥压强或压力比尺。
文献[8]在模拟坝前粘性淤积物局部冲刷漏斗时，正是按式(10)相似条件配制模型泥浆，并使模型与原型泥沙内摩擦角相同，进行冲刷漏斗试验，得到的试验成果较为合理。这表明若模型泥浆配制满足上述相似条件，就能保证模型与原型泥浆作用力相似。
3 模拟试验方法
满足式(10)相似条件的模型泥浆有多种选择。在几何比尺λL 已定的情况下，影响 模型泥浆选择的因素有极限剪应力τb 和泥浆容重γm .而影响τb 大小的因素又有泥浆浓度、细颗粒含量及絮凝结构 等。可见，τb 和γm 又相互影响。所以实际配制泥浆时需要反复调整泥浆体浓度和细颗粒含量，测定τb 和γm 值，使其比尺组合符合式(10)要求，并使模型与原型泥沙内摩擦角相同。
关于原型淤泥极限剪应力τb 、浆体容重γm 和内摩擦角的确定，由于在水库设计阶段无法实测淤泥参数，主要是根据坝前淤积物颗粒级配，参照类似工程实测资料并进行一些分析 计算 ，首先确定坝前淤泥干容重γ′，然后根据原型淤泥颗粒级配及干容重，选择合适的公式通过计算得到τb 和φ，或选择相同性质的天然沙样在试验室试配少量试样实际测定τb 和φ。泥浆容重γm 则根据下式计算
式中：γ为水容重。
将配制好的泥浆按门前设计淤积高程充填到模型中，再在其上部注清水至设计库水位，然后就可以利用闸门启闭装置进行启门力试验。在门前有泥沙淤积情况下，启门力的最大值通常发生在开始启门的瞬间。另外需要指出，由于泥浆的极限剪应力τb 与其配制时间的长短有关，所以配制好的泥浆应及时充填到模型中进行试验。否则τb 将发生变化，泥浆的相似性将遭到破坏。
模型上的闸门行走支承及止水型式等很难做到与原型相似，因而作用在行走支承及止水处的模型摩阻力与原型摩阻力是不相似的。为了消除摩阻力不相似因素的影响，在同一模型相同边界条件下，分别测定清水启门力和有泥沙淤积时启门力，并令
式中：Fsm 为有泥沙淤积时模型启门力；Ｆwm 为清水模型启门力，k为比例系数，大于1.
由动力相似原理知，模型的k值等于原型的k值。一旦通过模型试验确定出k值后，就可以在已知清水原型启门力Fwp 情况下，得出有泥沙淤积时原型启门力Fsp ，即
关于清水原型启门力，既可以用启门力公式计算，也可以在同一模型上通过试验确定。清水启门力的模型试验，可参见有关文献 [10、11]。
4 结语
本文在推导门前淤泥相似条件时，其前提条件是设淤泥为宾汉体泥浆。大多数水库坝前淤泥一般较细且未固结，可当作宾汉体泥浆。如果淤泥已固结或为散粒体粗沙，那么就属于土力学研究 范围，不能应用 该文所给出的淤泥相似条件进行模拟。
参 考 文 献：
[1] Xu Guobin, Bai Shilu. Modeling study on sedimentation configuration in front of the dam and in the pressure conduit of bottom outlets by flume experiment [J]。International Journal of Sediment Research, 1997,12(3).
[2] 夏毓常。关于泥沙淤积对闸门启闭力影响的估算方法[J]。人民黄河。1980,(3).
[3] 沙玉清。泥沙运动学引论[M]。北京：中国 工业 出版社，1965.
[4] 水利 电力部水利 水电规划设计院。水电站泥沙问题 总结 汇编[M]。1988.
[5] 钱宁，万兆惠。泥沙运动力学[M]。北京：科学 出版社，1983.
[6] 金腊华，石秀清。试论模型沙水下休止角[J]。泥沙研究，1990，(3).
[7] 王延贵，等。模型沙物理特性的试验研究及相似分析[J]。泥沙研究，1992,(3).
[8] 徐国宾，白世录。坝前粘性淤积物局部冲刷漏斗模拟相似性研究[J]。水利 学报，1997，(11).
[9] 希辛柯 РИ。泥浆水力学[M]。袁恩熙，陈家琅译。北京：石油工业出版社，1957.
[10] 南京水利 科学研究院，水利 水电科学研究院。水工 模型试验(第2版)[Ｍ]。北京：水利 电力出版社，1985.
[11] 哈焕文，郑大琼。快速闸门动水下降持住力的试验研究[C]。水利 水电科学研究院科学研究论文 集，第33集(水力学、泥沙、冷却水).北京：水利 电力出版社，1990.
减小字体
增大字体