少沙河流推移质计算方法.docVIP

少沙河流推移质计算方法 刘晓华 吴建军 (山东工业大学 济南 250061) 摘 要 将柘溪水库认为一天然实验水槽，利用柘溪水库近20年悬移质和推移质的实测资料与渣洋滩水文站入库的实测水力要素根据河流动力学原理，建立了少沙河流起动流速公式。从而推导出少沙河流推移质输沙模数，经12个水文站验算，成果十分合理，可推广到少沙河流上应用。另外，利用新建的起动流速公式与水力学连续方程、动力方程三式联解求得推移质淤积比降Jk。再用图解法经试算后，即可求得推移质与悬移质的分界点，大大节省野外勘测及实验室工作，能普遍应用到少沙河流的大、中、小型水库。 关键词 少沙河流 泥沙淤积 推移质 悬移质 1 引言 我国的少沙河流，水库淤积造成的损失是较大的。例如：山东省46座大中型水库，每年淤积库容0.61 亿 m3 ，加上小型水库的淤积，每年约淤积库容 1亿m3, 相当于每年损失一座大型水库，又如：丹江口水电站从 1960年到 1965 年之中就淤积泥沙 9?641.5万 m3 。在少沙河流上，水库淤积问题反映在通航、防洪和推移质颗粒对水轮机的磨损等方面的严重性和迫切性已经引起了重视。为了研究和解决水库泥沙淤积问题，从柘溪电站的实测资料和实验分析提出了少沙河流推移质的简便计算方法，曾在南方地区和山东省得到了验证，可在国内外推广应用。 2 柘溪水库的自然概况及运行情况 柘溪水库位于长江流域资水干流，控制流域面积22?640km2，年径流总量196亿m3。 原设计正常高水位167.5m (后提高至169.5m)，坝高104m，死水位 144.0m， 极限死水位140.m，总库容35.7亿m3，有效库容22.58亿m3，属河道型带状水库。 资水流域地处亚热带季风气候区，降雨集中在4~6月，占年总雨量的50%~60%，丰枯流量变化大，4~8月来水占全年水量的70%。流域地质为泥盆系， 石灰系、二叠系的石英砂岩、灰岩、硅质岩、页岩组成， 以红壤、黄壤、灰棕壤为主，覆盖较好，水土流失面积为1270km2，占总流域面积的5.6%，近坝库区为砂岩、页岩和板岩, 库尾为灰岩和红砂跞岩。多年平均悬移质输沙量为307万吨。 柘溪电站1974年12月当水库水位降落到死水位以下，籍近100km的库区悬移质和推移质采用篩分成果淤积暴露于水外的最佳时机，利用原设横断面对库内淤积体进行了测量，求得总淤积体3。在24个断面上进行淤积物取样，水上部分在两岸挖槽，分层取样。水下部分以戽斗捞取，对各断面的沙样分为7个粒径组进行现场篩分（水篩），其粒径为0.5、0.5～1.0、1.0～5.0、5.0～10.0、10.0～20.0、20.0～40.0mm，求其各组百分比，绘制各断面的级配曲线。同时，还取样回室内作滴定分析。当d0.5mm，采用篩分成果；当d≤0.5mm，采用滴定分析成果，结合现场勘测测量，明显地区别了悬移质和推移质大致的分界点是在CS14断面。干幺重γ分别为1.4t/m3（d0.01mm）和0.9t/m3(d0.01mm)。在采样桶内取水样1∽2cm3量积，浓缩、过滤、烘干、称重，求得γ=1.4t/m3占总量的59.5%，γ=0.9t/m3占总量的40.5%，即可求得11年总淤量为5426万吨。由各断面级配曲线求得推移质淤积量为726万吨，占总淤量14%，悬移质淤积量达4664万吨。 3 入库水文站测验情况 渣洋滩水文站于1956年12月设立，为柘溪专用站，1963年改为基本站，站址距坝158.01km，正处于回水末端。控制流域面积为16236km2，多年平均流量为424m3/s，多年平均悬移质243万吨，平均输沙率为77 kg/s。历年单沙与断沙关系良好。该站用横式采样器进行输沙率测量，全年测次不少于25次。用采样桶取水作悬移质泥沙颗粒分析和含沙量分析，全年均不少于50次。 4 长时段推移质输沙量的估算 将柘溪水库认定为一天然试验水槽，利用柘溪水库泥沙淤积成果，与入库站（渣洋滩）实测的水力要素建立数学模型，以寻求适用于各种不同特性的少沙河流推移质输沙量与入库站水力要素的一般规律，影响推移质运动的水力要素主要是流量、流速、水深、河宽泥沙粒径等。利用柘溪水库实测的河床质颗粒级配曲线（见图1）及入库站平均流速、水深及粒径（见表1）建立推移质起动流速公式。通常实测流速分布公式以指数型表示，应用河流动力学原理，经推导可得起动流速公式 (1) 由表1所列数据可得(由计算机算出，其程序及过程从略)。 (2) 式中: dα — 给定流量相应的颗粒平均粒径(