松花江白山河段冰情数值计算及分析-水科学进展.pdf

第２７卷 第６期 水 科 学 进 展 Ｖｏｌ．２７，Ｎｏ．６ ２０１６年１１月　 　 ＡＤＶＡＮＣＥＳ ＩＮ ＷＡＴＥＲ ＳＣＩＥＮＣＥ Ｎｏｖ．，２０１６ ＤＯＩ： １０１４０４２／ ｊｃｎｋｉ３２１３０９２０１６０６０１２ 松花江白山河段冰情数值计算及分析 １ １ ２ ３ 樊　 霖 ， 茅泽育 ， 吴剑疆 ， 齐文彪 （１ 清华大学水利水电工程系， 北京　 １０００８４；２ 水利部水利水电规划设计总院， 北京　 １００１２０； ３ 吉林省水利水电勘测设计研究院， 吉林 长春　 １３００２１） 摘要： 根据水动力学、 河流动力学、 热力学、 河冰水力学及固体力学等原理， 针对松花江上游白山河段具体特征， 建立了河冰数学模型， 并应用有限差分计算方法， 对白山河段冬季冰情演变进行了精确模拟。 应用该河段 １９５８— １９７３年共 １５年完整的地形、 气象、 水文、 热力、 冰情等原型观测资料， 分别进行了模型参数率定及冰情数值模 拟。 研究结果表明： 白山河段的封冻首先开始于白山坝址下游４ｋｍ 的大崴子河段， 然后封冻边缘逐渐上溯， 最终 到达松 １４断面； 计算的白山河段冰花堆积体外形与河床纵剖面呈相似趋势； 各种水力及冰情要素的数值计算结果 和实测值吻合较好， 所建立的数值模型能较好地模拟该河道的冰情。 研究结果对东北地区河流冬季冰情研究及冰 害防治具有一定的借鉴意义。 关键词： 河流冰情；河冰数值模拟；天然河道；松花江白山河段 中图分类号： ＴＶ１２２　 　 　 文献标志码： Ａ　 　 　 文章编号： １００１⁃６７９１（２０１６）０６⁃０８９０⁃０８ 河冰是冬季高纬度地区及寒冷区域河流中普遍存在的自然现象。 在寒冷的冬季， 河流中的水体经历十分 显著而且极为复杂的变化， 出现一系列的河冰现象， 例如河流中各种形式的冰情的形成、 演变、 输移、 耗散 和消融等。 这些现象通常是水动力、 机械力和热力等相互作用的结果， 同时还受到气象、 水文条件和河道地 形及河势的影响， 是一系列极其复杂的物理过程的组合。 对于寒带地区水利水电工程、 水资源工程的设计及 调度运行、 冬季航运交通、 引水发电、 给水排水以及涉及环境、 生态等的研究问题来说， 正确理解河冰演变 ［１］ 过程非常重要 。 多年来， 随着河冰水力学理论的发展， 国内外学者们提出了不少模拟河道冰情的数值模型。 例如， Ｌａｌ ［２］ ［３］ 和Ｓｈｅｎ 根据热力和两层输冰模型的理论分析框架， 建立了冰情演变河冰数值模型ＲＩＣＥ模型。 Ｓｈｅｎ等 提出了二维与一维结合的河冰动力学模型ＤｙｎａＲＩＣＥ， 主要是模拟冰坝形成的动力学过程等。 Ｚｕｆｅｌｔ 和Ｅｔｔｅ⁃ ［４］ ［５⁃６］ ｍａ 建立了一维冰水耦合的冰坝动力学模型， 模拟封河期和开河期水位变化。 茅泽育等 针对天然河道弯 曲复杂的特点， 首次建立了适体坐标下的二维河冰数值模型， 对黄河河曲段从龙口至禹庙５６ｋｍ河段进行了 数值模拟及验证， 取得了较好的效果， 并将模型应用于万家寨水库库区冰凌、 吉林中部城市供水长春干线伊 ［７］ 丹河冰期输水研究等。 Ｆｕ等 在ＲＩＣＥ和ＲＩＣＥＮ模型的基础上， 开发了用于模拟冰期洪水的ＹＲＩＤＭ（ＹＲＣＣ Ｒｉｖｅｒ Ｉｃｅ Ｄｙｎａｍｉｃ Ｍｏｄｅｌ）模型， 并用 ２００８—２００９ 年冬季黄河宁蒙河段的冰情数据对模型的精确度进行了 验证。 如上所述， 河流冰情演变受到水力、 气象、 热力、 河道特征等众多因素的影响， 这些众多影响因素的不 同组合导致