简化计算方法在泥沙分析中的应用

《简化计算方法在泥沙分析中的应用》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、简化计算方法在泥沙分析中的应用陈光洪李学通（四川省水利水电勘测设计研究院规划设计分院，德阳618000）【摘要】大中型水利水电工程泥沙分析计算通常采用二维或三维数学模型，对资料要求较高，计算较为复杂。根据工程特点和资料条件，可通过简化计算方法进行泥沙分析。本文以撒鱼沱电站水库的泥沙计算分析为例，结合与数学模型计算成果的对比，分析了简化计算方法的可行性。【关键词】电站；水库；泥沙分析；简化计算；应用列，得到1957～2000年横江水沙系列。豆沙关水文站缺测1978年资料，采用本站水沙相关分析而得。白水江牛

2、街水文站有1966～2000年实测资料。根据三站1966～2000年同期悬沙资料统计各代表站水沙系列平均特征值和多年平均悬移质输沙量年内分配，各代表站多年平均悬移质输沙量年内分配见图1。1概况撒鱼沱电站是横江干流六级河床式开发的第一级，上距盐津新县城3.5km。坝址处控制集雨面积13339km2，多年平均流量240m3/s。电站正常蓄水位414.0m，死水位3411.6m，日调节库容267万m，装机容量60.0MW，多年平均发电量2.84亿kW·h。电站闸坝最大壅水高约28.0m，布置3孔8m宽的冲沙闸

3、和6孔12m宽的泄洪闸，泄洪总宽度96m，采取汛期来流量大于分界流量900m3/s时电站停机拉沙的运行方式。电站水库为河道型水库，库区天然河道比降约为2.27‰，回水长度约10.6km，平均水面宽约96m，库区河道断面形态以“V”型为主。水库多年平均入库悬移质输沙量1125万t（合803.6万m3），推移质112.5万t（合51.1万m3），库沙比1.22，工程泥沙问题较严重。根据横江、豆沙关及牛街水文站1966～2000年同期悬沙资料，求得区间悬移质输沙模数，用区间面积求出区间产沙量，从而求出撒鱼沱梯

4、级入库悬沙量为1125万t。2悬移质特征值1）悬移质泥沙颗粒级配根据横江站悬移质泥沙观测资料，横江站悬沙D50=0.0052mm，dcp=0.017mm，dmax=0.8mm，横江站悬移质泥沙颗粒级配成果见图2。2泥沙分析简化计算2.1基本资料2.1.1悬移质1悬移质输沙量工程河段上下游及支流具有较长泥沙资料的水文站有干流的横江、豆沙关站，支流白水江的牛街水文站。横江水文站控制集水面积14781km，控制了横江流域面积的98.7%，从1957年至今，除1961年～1964年缺测泥沙外，其余均有实测水沙资

5、料，其实测含沙量在0～34.3kg/m3之间。通过横江站实测水沙系列作相关分析，得水沙相关方程（此略），用相关方程插补缺测泥沙资料后再加实测系22）分界粒径dk据造床流量（敞泄时坝前水位400.00m）坝前流速、水深，用式（2）计算悬移质泥沙分界粒径，按计算分界粒径重新对悬沙级配进行统计，求平均粒径、中值粒径，并采用式（3）计算悬沙中床沙质沉降速度ω。计算成果见表1。dk?0.052(v0.85h)（2）??d2g(rs?r?0.674d218rv)v（3）式中：ν为浑水粘滞系数（可查表）；γs为沙容重

6、；γ为水容重；d为悬沙中床沙质的中值粒径。表1建库后悬沙特征值造床流量流速V（m/s）2.72水深h（m）7.38分界粒径dk（mm）0.0223悬沙中平均粒径Dcp（mm）0.0659床沙质中值粒径d50（mm）0.0435部分浑水粘滞系数ν（10-6）1.0324沉降速度ω0.00124冲泻质（%）：床沙质（%）81.8：18.22.1.2推移质１库区床沙特征根据撒鱼沱库区床沙取样分析，库区床沙dcp=133mm，d50=77mm，dmax=400mm，床沙级配曲线见图3。２推移质沙量采用沙莫夫公式

7、计算：V3H?3.83dh6（4）gb?KD2(V?VV3H)(dV)()4Hh（5）式中：gb为单位宽度的推移质输沙率；D为非均匀沙最粗一组的平均粒径；K=2.5。经计算，撒鱼沱梯级库区推移质输沙量为71.9万t（合32.7万m3），小于悬移质输沙量的10%。参考专家意见偏安全考虑，取用悬移质输沙量的10%，即推移质输沙量112.5万t（合51.1万m3）。2.2泥沙分析简化计算2.2.1悬移质1基本参数1）造床流量Q造的选择根据经验，造床流量一般相当于多年平均出现2天～3天的流量级。若按多年平均出现

8、3天以上的日平均流量选择，为1730m3/s；若按多年平均出现2天以上的日平均流量，为1840m3/s。即撒鱼沱梯级造床流量应在1730m3/s～1840m3/s左右。若用黄委水科所公式：Q0.85造?7.7Q汛?90Q汛3（6）计算，则造床流量为1868m3/s（汛期平均流量384.7m3/s），与按经验选择的结果较为接近。因此，撒鱼沱梯级造床流量采用黄委水科所公式计算结果。2）挟沙能力S*的确定采用张瑞瑾公式计算：3S*?k（VmgR?