一维泥石流的静动力阻力特征研究及数值模拟-工程地质学报.PDF

ＪｏｕｒｎａｌｏｆＥｎｇｉｎｅｅｒｉｎｇＧｅｏｌｏｇｙ　　工程地质学报　　１００４－９６６５／２０１０／１８（６）０８５７０５ 一维泥石流的静动力阻力特征研究及数值模拟 ① ①② ① 樊 　王思敬 　王恩志 （ 清华大学　水沙科学与水利水电工程国家重点实验室　北京　１０００８６） ① （ 中国科学院　地质与地球物理研究所　北京　１０００２９） ② 摘　要　泥石流是山区多发的一种地质灾害，它的发生和发展威胁着人们的生命和财产安全，影响着人们的正常生活，因而 需要加强对其发生和发展过程的研究。结合泥石流的动力模型方程采用数值模拟方法再现泥石流发生和发展的过程，是研 究和预测模拟泥石流灾害的有效手段。目前的动力模型方程大多只关注动力过程，却忽视了静动力过程的统一，这将导致在 一些情况下产生错误的结果。本文研究了一维泥石流的静动力阻力特征，通过修正泥石流动力学方程的阻力项，得到了具有 静动力统一特征的模型方程。并以Ｒｏｅ格式的近似Ｒｉｅｍａｎｎ解为基础，采用ＭＵＳＣＬ线性重构方法建立了具有较高精度和分 辨率的有限体积数值求解。具体算例的数值验证表明，方程阻力项的修正是合理的，所建立的数值求解也是稳定和有效的。 关键词　泥石流　阻力特征　Ｒｏｅ格式 中图分类号：Ｐ６４２２，Ｘ１４１　　文献标识码：Ａ ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＴＡＴＩＣＡＮＤＤＹＮＡＭＩＣＲＥＳＩＳＴＡＮＣＥＯＦＯＮＥ ＤＩＭＥＮＳＩＯＮＡＬＤＥＢＲＩＳＦＬＯＷ ＡＮＤＩＴＳＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴＩＯＮ ① ①② ① ＦＡＮＹｕｎｙｕｎ　ＷＡＮＧＳｉｊｉｎｇ 　ＷＡＮＧＥｎｚｈｉ （ ＳｔａｔｅＫｅｙＬａｂｏｒａｔｏｒｙｏｆＨｙｄｒｏｓｃｉｅｎｃｅａｎｄＥｎｇｉｎｅｅｒｉｎｇ，ＴｓｉｎｇｈｕａＵｎｉｖｅｒｓｉｔｙ，Ｂｅｉｊｉｎｇ　１０００８４） ① （②ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｌｏｇｙａｎｄＧｅｏｐｈｙｓｉｃｓ，ＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ，Ｂｅｉｊｉｎｇ　１０００２９） Ａｂｓｔｒａｃｔ　Ｄｅｂｒｉｓｆｌｏｗｉｓｏｎｅｏｆｔｈｅｃｏｍｍｏｎｍｏｕｎｔａｉｎｈａｚａｒｄｓａｎｄｓｅｒｉｏｕｓｌｙｔｈｒｅａｔｅｎｓｈｕｍａｎｌｉｖｅｓａｎｄｂｅｌｏｎｇｉｎｇｓ ａｎｄａｆｆｅｃｔｓｓｏｃｉｅｔｙｄｅｖｅｌｏｐｍｅｎｔ．Ｔｈｅｒｅｆｏｒｅ，ｉｔｉｓｎｅｃｅｓｓａｒｙｔｏｓｔｕｄｙｍｏｒｅａｂｏｕｔｉｔｓｈａｐｐｅｎｉｎｇａｎｄｄｅｖｅｌｏｐｍｅｎｔ．Ｔｈｅ ａｐｐｒｏａｃｈｏｆｃｏｍｂｉｎｉｎｇｔｈｅｄｙｎａｍｉｃｍｏｄｅｌｅｑｕａｔｉｏｎｓｗｉｔｈｔｈｅｕｓｅｏｆｔｈｅｎｕｍｅｒｉｃａｌｍｅｔｈｏｄｓｃａｎｌｅａｄｔｏｉｎｄｅｐｔｈｕｎ ｄｅｒｓｔａｎｄｉｎｇｏｆｔｈｅｄｅｖｅｌｏｐｍｅｎｔｏｆｄｅｂｒｉｓｆｌｏｗ．Ｔｈｉｓｉｓａｎｅｆｆｅｃｔｉｖｅａｐｐｒｏａｃｈｔｏｓｔｕｄｙａｎｄｆｏｒｅｃａｓｔｔｈｅｄｅｂｒｉｓｆｌｏｗ ｄｉｓａｓｔｅｒｓ．Ｈｏｗｅｖｅｒ，ａｍａｊｏｒｉｔｙｏｆｔｈｅｐｒｅｓｅｎｔｄｙｎａｍｉｃｍｏｄｅｌｅｑｕａｔｉｏｎｓｆｏｃｕｓｏｎｔｈｅｄｙｎａｍｉｃｐｒｏｃｅｓｓ，ａｎｄｉｇｎｏｒｅｔｈｅ ｃｏｍｂｉｎａｔｉｏｎｏｆｔｈｅｓｔａｔｉｃａｎｄｄｙｎａｍｉｃｐｒｏｃｅｓｓ．Ｉｎｆａｃｔ，ｔｈｉｓｎｅｇｌｅｃｔｍａｙｒｅｓｕｌｔｉｎｓｏｍｅｍｉｓｌｅａｄｉｎｇｒｅｓｕｌｔｓｕｎｄｅｒ ｓｏｍｅｃｏｎｄｉｔｉｏｎ．Ｉｎｔｈｉｓｐａｐｅｒ，ｔｈｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｓｔａｔｉｃａｎｄｄｙｎａｍｉｃｒｅｓｉｓｔａｎｃｅｏｆｏｎｅｄｉｍｅｎｓｉｏｎａｌｄｅｂｒｉｓｆｌｏｗａｒｅ ｓｔｕｄｉｅｄ．Ｂｙｍｏｄｉｆｙｉｎｇｔｈｅｆｒｉｃｔｉｏｎｒｅｓｉｓｔａｎｃｅｔｅｒｍｏｆｔｈｅｋｉｎｅｔｉｃｅｑｕａｔｉｏｎ，ｔｈｅｍｏｄｅｌｅｑｕａｔｉｏｎｉｓｏｂｔａｉｎｅｄｆｏｒｔｈｅｕｎｉ ｆｏｒｍｃｈａｒａｃｔｅｒｉｓｔｉｃｏｆｓｔａｔｉｃａｎｄｄｙｎａｍｉｃ．ＢａｓｅｄｏｎｔｈｅａｐｐｒｏｘｉｍａｔｅＲｉｅｍａｎｎｓｏｌｖｅｒｏｆＲｏｅｓｃｈｅｍｅ，ｔｈｅｎｕｍｅｒｉｃａｌ ｓｏｌｕｔｉｏｎｗｉｔｈｈｉｇｈｅｒｐｒｅｃｉｓｉｏｎａｎｄｒｅｓｏｌｕｔｉｏｎｉｓａｃｈｉｅｖｅｄｂｙｕｓｅｏｆＭＵＳＣＬｌｉｎｅａｒｒｅｃｏｎｓｔｒｕｃｔｉｏｎ．Ｔｈｅｎｕｍｅｒｉｃａｌｖａｌｉ ｄａｔｉｏｎｏｆａｎａｃｔｕａｌｅｘａｍｐｌｅｐｒｏｖｅｓｔｈｅｒａｔｉｏｎａｌｉｔｙｏｆｔｈｅｍｏｄｉｆｉｃ