河流冰塞数值模拟进展 Recent development in simulation of river ice jam.pdfVIP

第19卷第4期 水科学进展 V01．19．No．4 2008年7月 ADVANCESINWA7rERSCIENCE JIll．，2008 河流冰塞数值模拟进展 王 军，赵慧敏 (合肥工业大学土木学院，安徽合肥230009) 摘要：冰塞是高纬度地区河流凌汛的主要诱因。为防凌减灾，围绕着河冰问题的试验研究、原型观测和数值模拟， 国内外学者展开了大量的研究。就数值模拟方面，综述了国内外学者在河冰一维、二维数值模拟方面的研究成果； 重点描述了Shen的一维模型和模拟冰塞体溃决的二维模型及其在实际中的应用；在评述和分析各模型特点的基础 上，尝试提出了今后的发展方向。 关键词：冰塞；数值模拟；耦合；热力 中图分类号：P332．8；P338．4；G353．11文献标识码：A 文章编号：1001．6791【2008)04．0597．08 冬季，北方地区河流常因冰塞而成灾，围绕着冰塞灾害问题。国内外众多学者从试验研究、原型观测、数 值模拟等不同的方面进行研究，取得了大量有益的成果。就数值模拟方面而言，从20世纪70年代开始，首先 出现了一维河冰数值模拟，迄今已发展出现了二维河冰数值模拟，并在工程实际中得以应用。本文就该方面的 研究成果进行综述，并尝试提出了今后的发展方向。 1现有模型分析 1．1 一维模型 近年来，新近研究所建立的一维河冰模型主要有zufelt…、She和Hicl【s怛J的研究成果。 Zufeh等建立了一维输冰耦合模型，该模型的基本方程分为冰水质量守恒方程和动量守恒方程。 利用该模型，Zufelt等研究了冰的动量和水的非恒定性对冰塞厚度值大小及其分布的影响，得出了由 非均匀性是造成冰塞前缘头部较大的原因，虽然该结论未考虑到水流条件与输冰量相互之间的调整关系，但已 经打破了千篇一律的冰塞厚度分布常规的描述，即较薄的上游过渡段；厚度均匀的平衡段；坚冰盖前的厚度最 大的冰趾段。 通过研究，Zufelt等还提出了考虑冰动量和水流非均匀性与否的两个主要的无量纲数。 Zufeh等耦合求解了一维冰水运动方程，模拟了冰水运动及其对冰塞厚度的影响，这种冰水偶合解法在一 维模型中是一个进步。 Zufelt等模型的不足在于未考虑当冰塞运动或坍塌时，冰塞强度以及冰塞与岸边摩擦系数的变化。 混合形式，单独列出冰质量守恒方程且非耦合求解该方程。 冰质量纵向释放扩散效果，研究了在冰塞体溃决时冰对洪水波推进和冰塞体再次形成的影响。为检验该模型， 收稿日期：2007．06．15 基金项目：国家自然科学基金资助项目 万方数据 598 水科学进展 第19卷 John 首先同Shen和Liu[4]的算例进行了对比，结果吻合较好。使用该模型对Saint 拟表明，冰的影响主要体现在波的释放的第一阶段，而对Athabasca 型中考虑冰的影响并未对模拟产生改善的效果，虽然模拟取得了成功，但结论却有点似是而非。 Tao Hung 对较为成熟的阶段。 Shen[5】基于热交换的相关原理，提出了河流冰盖消长模型，模型将河流系统视为气．冰．水．床系统，考虑了 所有交界面上的热交换和冰盖的热传导，对于热传导方程采用一维准稳态的形式，水温分布采用纵向扩散方程 表示，在冰盖的热力消长方面，Shen提出了冰盖的热交换受水内冰花堆积和雪盖的影响。该模型较为全面的 河上的44．7km河段的冰盖厚度和水温进行了模拟。 用于模拟河湖冰盖的热力消长，提出计算公式： h=AOSl尼 (1) 式中h为冰塞厚度；A。为与特定场所有关的经验系数；S为累积负气温。 Shen的研究中提出了一种改进的度一13一法，其计算公式为 h=(h3一aS)1