西宁市人工湖工程水流泥沙二维数模计算研究.pdfVIP

十、河流模拟技术 。1465。 西宁市人工湖工程水流泥沙二维数模计算 李国斌 陈金杭 王丽芝 高亚军 f南京水利科学研究院 南京 210024) 摘要：本文运用正交曲线拟合坐标下的平面二维水沙运动数学模型对西宁市人工湖工程水流泥沙进行 了计算，计算中考虑了高含沙量条件下泥沙沉速、分组挟沙能力等，计算得到百年一遇洪峰流量时人工 湖沿程水位及遭遇1981年型水沙过程考虑上游减沙与否两种情况下沉沙池、人工湖的泥沙淤积量和淤 积分布，研究成果为工程设计提供了依据。 关键词：正交曲线拟合坐标 平面二维 水沙运动 数学模型 泥沙淤积 1前言 西宁市为青海省的省会城市，是多民族共同聚居的高原名城。西宁市政府根据国家 西部开发的宏观战略，结合自身的社会和经济发展需要，计划在市区东郊附近(技术开发 区北侧)，利用湟水河沙塘川入口一小峡121河段兴建景IX(人工湖)，由人工湖面和湖岸景 物形成山水相映、塞北江南的优美环境，并依此作为西宁市东大门的标志性景观，从而 提高西宁市的人居环境质量，改善区域投资环境，促进当地经济发展。为论证人工湖工 程建成后水流泥沙运动状况，我们建立了水流泥沙二维数学模型，计算研究了遭遇百年 一遇洪峰流量时人工湖沿程水面线及流场、遭遇1981年型水沙过程考虑上游减沙与否两 种情况下人工湖本身及上游沉沙池的泥沙淤积量及淤积分布，其成果可供工程决策和运 行参考。 2二维水沙运动数学模型 2．1基本方程组”1 针对人工湖河道蜿蜒曲折的特点，本研究选用正交贴体曲线坐标拟合边界，以便较 好地模拟河道边界，又使废网格大大减少。为此，首先依据势流的流线与等势线正交原 理生成正交曲线网格。 2．1．1水流运动方程 连续方程： … az+上f虫型+型]：o 、’ Ot c≥c0 a善 a刁 作者简介：李国斌，男，1963年出生，教授级高级工程师。 ·1466· 第六届全国泥沙基本理论研究学术讨论会论文集 动量方程 掣+壶降咖曲+南睡酬m髻一肼割=一呼一器+ 壶砖心‰)+南候％小№。等一％。割 c2， 掣+壶降咖_I)+南睡酬+协髫一肼雳]=一呼一等考+ 壶睦(钟嘞)+南睡锄o+‰等一鲰筹l c3， 式中：甜、V分别为善、刁方向流速分量；Z为水位；H为水深；C为谢才系数；眙、 仃仰、吒口、盯描为紊动切应力，h为紊动黏性系数，¨=G竺，七为紊动动能，占为紊 动动能耗散率，可应用k一占紊流模型计算确定；一般情况下，u：Ku．H，K：0．5～1．0。 玑为摩阻流速，或者取u为某一常数。 2．1．2悬移质不平衡输沙方程 根据上游来沙情况，将非均匀悬沙按其粒径大小分成NS组，用S表示第f组粒径 的含沙量，用只表示此组粒径悬沙在悬沙总含沙量S中所占的百分数，则有： Ns Si=只S s=∑St (4) f=l 对于非均匀悬沙中第i组粒径的含沙量，可以写出窦国仁、赵士清[21提出的水平二 维分组悬沙不平衡输沙方程式如下 考字+上qq『L旦a；b协刚+南心啪)] @’ =去[毒(詈昌嚣)+南(毒苦器]]+伪劬届研一墨) 式中：日，、Oi、筇分别为分组粒径恢复饱和系数、沉速和挟沙力；鬈、勖分别