可变网格下的水库滑坡涌浪数值模拟研究 Horizontal 2-D flow model with variable grid for simulating surges due to landslide in reservoirs.pdfVIP

第19卷第4期 水科学进展 V01．19。No．4 2008年7月 ADVANCESINWATElt JuJ．，2008 SCIENCE 可变网格下的水库滑坡涌浪数值模拟研究 袁 晶1，张小峰2，张 为2 (1．长江水利委员会水文局。湖北武汉430010； 2．武汉大学水利水电学院水资源与水电工程科学国家重点实验室，湖北武汉430072) 摘要：采用动网格技术，建立了可变网格下的平面二维水库滑坡涌浪数学模型，与传统固定网格下的数学模型相 比，克服了难以精确计算水库滑坡体对水体作用的困难。以一模型试验为研究对象，采用数学模型对其进行了复 演，模型将滑坡体对水体的作用采用可变网格的思想实现，使滑坡体所占区域的计算网格随滑体下滑。模型计算值 与实测值比较表明，两者吻合较好，说明基于可变网格思想的数学模型能够较准确地模拟滑坡涌浪过程，为滑坡涌 浪的计算、预测提供了一种新的研究方法。 关键词：滑坡；涌浪；水库；可变网格；数学模型 中图分类号：TVl39．232文献标识码：A 文章编号：1001．6791(2008)04．0546—06 水库的修建将改变库区边岸的自然地质条件，其边岸的岩体在水库波浪及水位变化等因素的共同作用下， 将发生塌陷，甚至滑坡。其中尤其以近坝的高速滑坡对水库的危害最大。如1963年意大利瓦依昂水库滑坡…， XlOs 超过2．4 m3的滑坡体滑入水库激起约100m的涌浪漫过坝顶，摧毁坝下游一个村庄，致使2000多人丧生。 1961年我国的柘溪水库塘岩光滑坡【2J，2000年西藏易贡巨型滑坡13J，均给人民生命财产造成了不同程度的损 失。三峡大坝下闸蓄水至135m以来，三峡库区两岸已有几处滑坡体失稳下滑，随着水库蓄水位的抬高，将有 更多的潜在滑坡体可能因目前的稳定条件被破坏而失稳下滑【4J。研究滑坡涌浪的传播过程，对滑坡涌浪的预 报。边坡加固设计与防范工程措施的实施，水库运行决策等，均有十分重要的现实意义【5J。 数学模型以其计算速度快，精度高的优点被广泛地应用于滑坡涌浪的计算中。金峰和姜治兵L41建立了非规 则网格下的数学模型，对一概化地形进行了滑坡涌浪计算，同时考虑了体积守恒与水平方向动量守恒来计算滑 体对水体的作用。严骏龙[6】应用二维有限元分析方法，引入质量集中的概念对一概化矩形水槽进行了滑坡涌浪 模拟，将滑体对水体的作用分解成水平分力和垂直分力进行计算，这一处理方法与Yong等【7】的思想类似，汪 洋【8’9]从地面运动和水下运动两个阶段对水库库岸滑坡的受力特征进行了分析，用牛顿定律和运动学基本原理 求解滑坡地面运动的加速度和速度，并运用动量定理对其初始涌浪高度进行了计算。由于问题的复杂性，以上 方法在计算滑体对水体水平作用的过程中，计算参数的选取具有较大的随机性和主观性。 采用动网格技术模拟了挑流冲刷时因河床边界变化引起的流场变化过程。本文采用动网格技术，建立了可变网格 下的二维水库滑坡涌浪数学模型，以一模型试验为研究对象，将滑体对水体的作用采用可变网格的思想实现，随 时间变动计算网格，使计算网格随滑体下滑，达到真实反映物理运动过程，减小模型计算误差的目的。 收稿日期：2007．08．08 源与水电工程科学国家重点实验室开放基金(2007C016) 作者简介：衰晶(1979一)，女。湖北潜江人，工程师，博士，主要从事河流水沙数值模拟及河道演变等研究。 F_,-mail：删@cjh．COrn 万方数据 第4期 袁 晶等：可变网格下的水库滑坡涌浪数值模拟研究 547 1模型的建立 1．1可变网格计算的基本原理 数值计算中，在物理空间里，点(茗，Y，f)在计算空间里可表示为(搴，叩，r)，即 r=t ∈=e(鼻，Y，t)叩=17(石，，，，t) 式(1)的3个方程分别对t，茗和Y求偏导得：