可变网格下的水库滑坡涌浪数值模拟研究-水科学进展.PDF

第 19 卷 第 4 期 水 科 学 进 展 Vol19 ,No4 2008 年 7 月 　　 ADVANCES IN WATER SCIENCE 　 Jul . , 2008 　 可变网格下的水库滑坡涌浪数值模拟研究 1 2 2 袁 　晶 , 张小峰 , 张 　为 ( 1 长江水利委员会水文局 , 湖北 武汉 　430010 ; 2 武汉大学水利水电学院水资源与水电工程科学国家重点实验室 , 湖北 武汉 　430072) 摘要 : 采用动网格技术 , 建立了可变网格下的平面二维水库滑坡涌浪数学模型 , 与传统固定网格下的数学模型相 比 , 克服了难以精确计算水库滑坡体对水体作用的困难 。以一模型试验为研究对象 , 采用数学模型对其进行了复 演 , 模型将滑坡体对水体的作用采用可变网格的思想实现 , 使滑坡体所占区域的计算网格随滑体下滑 。模型计算值 与实测值比较表明 , 两者吻合较好 , 说明基于可变网格思想的数学模型能够较准确地模拟滑坡涌浪过程 , 为滑坡涌 浪的计算 、预测提供了一种新的研究方法 。 关 　键 　词 : 滑坡 ; 涌浪 ; 水库 ; 可变网格 ; 数学模型 ( ) 中图分类号 : TV139232 　　　文献标识码 : A 　　　文章编号 : 2008 水库的修建将改变库区边岸的自然地质条件 , 其边岸的岩体在水库波浪及水位变化等因素的共同作用下 , 将发生塌陷 , 甚至滑坡 。其中尤其以近坝的高速滑坡对水库的危害最大 。如 1963 年意大利瓦依昂水库滑坡[ 1] , 超过 24 ×105 m3 的滑坡体滑入水库激起约 100 m 的涌浪漫过坝顶 , 摧毁坝下游一个村庄 , 致使 2 000 多人丧生 。 1961 年我国的柘溪水库塘岩光滑坡[2 ] , 2000 年西藏易贡巨型滑坡[3 ] , 均给人民生命财产造成了不同程度的损 失 。三峡大坝下闸蓄水至 135 m 以来 , 三峡库区两岸已有几处滑坡体失稳下滑 , 随着水库蓄水位的抬高 , 将有 更多的潜在滑坡体可能因目前的稳定条件被破坏而失稳下滑[4 ] 。研究滑坡涌浪的传播过程 , 对滑坡涌浪的预 报 , 边坡加固设计与防范工程措施的实施 , 水库运行决策等 , 均有十分重要的现实意义[5 ] 。 [4 ] 数学模型以其计算速度快 , 精度高的优点被广泛地应用于滑坡涌浪的计算中。金峰和姜治兵 建立了非规 则网格下的数学模型 , 对一概化地形进行了滑坡涌浪计算 , 同时考虑了体积守恒与水平方向动量守恒来计算滑