黄河下游河道冰动力行为二维数值模拟研究（Two dimensional numerical simulation of ice dynamic behavior in the lower reaches of the Yellow River River）.pdfVIP

黄河下游河道冰动力行为二维数值模拟研究 摘 要 水在地球表面的分布范围非常广，江河、湖泊和海洋等水体约占地球表面积 的 3/4，而高纬度地区的表面水体在冬季温度下降时经常会发生结冰现象。全球 冰问题总体来说可分为海冰问题和河冰问题，由于河流沿岸区域居住着大量居 民，使河冰问题显得尤为突出。 冬季河流封冻后，水体表面出现固体冰盖，使得过水断面的湿周增大，水力 半径减小，从而导致了水流形态的变化，带来的影响包括水力发电产量的损失， 航运线路阻塞、航期缩短，河流中的水利结构破坏，严重时形成冰塞、冰坝并导 致凌洪灾害，淹没农田房屋，危害人民生命和财产安全。我国北方多数河流都存 在不同程度的冰情问题，尤其是黄河下游河道是举世闻名的地上悬河，河道上宽 下窄，流向由低纬度向高纬度，南北纬度相差 3°，冰情现象十分复杂，河冰灾 害发生频繁且破坏力巨大。因此，对黄河下游河冰的形成、演变和输移过程进行 系统、深入的研究对提前制定防灾措施具有重大意义。 本文分析了冬季河冰结冰期、封冻期和解冻期的演变过程及河冰类型。进一 步探究黄河下游冰塞和冰坝现象的特点、形成条件和变化过程，并总结了黄河下 游开河时容易形成冰坝的原因。推导了河流无冰盖、部分被冰盖覆盖和完全被冰 盖覆盖情况下的水深公式。通过对冰盖上游漂浮冰块和潜入冰盖下的冰块进行受 力分析，得到了冰盖前冰块下潜的临界流速和下潜冰块在冰盖下继续向下游输移 的临界流速。 开展物理模型试验对不同流速、水深、冰盖类型和冰盖糙率等条件下冰盖流 垂线流速进行观测和分析，得到的主要结论包括： （1）冰盖下的垂线流速呈现两头小中间大的 U 字形态，光滑冰盖模型的水 深大于同流量明流水深。床面区光滑冰盖流的流速比同位置的明流流速稍大，但 随着水深增加两者越来越接近。流量相等时光滑冰盖流的最大流速大于相同位置 的明流流速。 （2 ）粗糙冰盖流垂线最大流速值随冰盖糙率的增大而增大，最大流速点更 靠近光滑的床面一侧。在不同水深情况下，最大流速点与床面之间的距离随着冰 盖与床面糙率比的增大而减小，基本呈现出线性变化规律。 I （3 ）岸冰河流中心开敞区水深大于明流，垂线流速分布形式与明流相同， 但相同深度上的流速比明流大，并且流速随着岸冰冰盖宽度的增加而增大。岸冰 冰盖流在两岸冰盖处的水深小于相同情况的明流水深，垂线流速分布形态呈现出 中间小两端大的形式。 对用于模拟复杂河道地形中冰运动和阻塞过程的二维冰动力学数值方法进 行改进，使用欧拉有限元法进行水动力模拟，冰运动模拟使用基于光滑粒子流体 动力学（SPH ）的拉格朗日离散元法（DPM ），建立计算黄河下游河冰输移过程 的二维数值模型。 通过开展物理模型试验以验证二维冰动力学改进模型的准确性。也就是控制 入口处的输入条件与数学模型相同，对比河道内冰密集度和典型断面水深的数模 结果和实测结果，分析表明：模型能较好地模拟出了河道内的冰输移过程，且两 者水深值比较吻合。 首先将收集的地形、水深和气象等数据输入模型，模拟了黄河下游济南市徐 庄处 1970 年 1 月发生的冰坝现象。通过对比模拟河段内冰坝长度、厚度和水深 的实测值与模拟值，表明该模型较准确地模拟了冰坝形成和发展的过程。其次， 以河道狭窄弯曲、易发生冰塞、冰坝的黄河下游北店子水文站至泺口水文站河段 为模拟区域，通过冰岛的设定模拟了 2010 年冬季岸冰存在情况下该河段冰厚度 变化。结果表明水流量较小时较强的风力作用对冰的拖曳力起到主导作用，且岸 冰的存在缩小了开敞河道宽度，阻碍大块表面冰顺利通过；而在岸冰后平铺形成 冰盖并向上游延伸。模拟了 2013 年 1 月该河段的流凌情况，可以发现流冰下泄 比较顺畅，未发生冰阻塞，且泺口水文站冰密集度的数模值与实测值比较一致。 由于模拟河段S 型弯的地形易发生冰堆积，为促进河道内稳定冰盖的生成，并延 长泺口浮桥的使用时间，提出在泺口水文站前河道内设置拦冰栅，使用 2013 年 1 月的输入条件，模拟该河段增加拦冰栅后的冰运动情况，结果显示表面冰在拦 冰栅后平铺并逐渐形成稳定冰盖，并不断向上游延伸。 为便于操作，本文建立了黄河下游冰动力模拟系统，使用户在界面上可得到 清晰直观的输出结果。 关键词：黄河下游；河冰；冰动力；试验研究；数值模拟