河口混合过程的研究3.pdfVIP

( ) 中　国　科　学 A 　辑 　 第 29 卷 　第 9 期 SCIENCE IN CHINA ( Series A) 　 1999 年 9 月 河口混合过程的研究 周济福　刘青泉　李家春 ( 中国科学院力学研究所,北京 100080) 摘要 　　根据河口水流运动特点 ,应用平板振荡边界层理论及波流分解方法 ,导出了 往复运动水流的流速垂向结构 ,据此建立了河口准二维盐度数学模型 ,并得到实测资 料验证. 应用该模型研究河口混合过程 ,得到了盐度分布、盐度锋强度随径流和潮差 定量变化的规律. 模型具有物理机理清晰 ,所需 CPU 时间短的优点. 关键词 　　河口混合 　振荡边界层 　径流潮流分解 　盐度 河口为盐淡水交会处 ,河口混合通过絮凝影响泥沙输运[1 ] ,是形成拦门沙的原因之一. 盐 水入侵也是河口地区农田水利、土壤改良、城市供水等问题中不可忽视的因素. 因此 ,河口混 合的研究越来越受到人们的关注. 河口混合过程因径流和潮流的相对强弱而异 ,大体上可分为 3 种类型 , 即弱混合型、缓混 合型和强混合型. 在弱混合型河口 ,潮差较小 ,径流作用强 ,表底层盐度差别很大 ,往往出现盐 水楔而形成分层流动 ,现有研究主要是建立盐水楔的几何尺寸和位置与径流、潮流强度之间的 ( [2 ,3 ] [4 ]) 关系 , 以及探讨盐淡水交界面掺混过程的机理 Grubert 和 Kurup ;而在强混合型河口 ,潮 差较大 ,潮流作用强 ,盐度沿垂线几乎均匀分布 ,可用较简单的模式或一维数学模型进行求解 ; 韩曾萃等人[5 ] 曾用一维水量、动量和盐量平衡方程模拟钱塘江河口的盐水入侵;缓混合型河 口,盐度的分布特征则介于两者之间 ,三维特性较明显 ,流动的垂向结构和盐度的垂向分布变 得重要 ,应通过动力学方法求解. 长江口即属于这一类型的河口. ( ) 严格地说 ,河口水流及标量 如泥沙、盐度等 输运都应采用完整的三维数学模型 ,如 :匡翠 1) 萍 曾运用三维数学模型对长江口水流泥沙盐度进行模拟. 但因其涉及物理因素多 ,不易分 析机理 ,为避免计算量过大 ,垂直方向网格划分不细亦会带来误差. 因此 , 目前实际应用中仍 普遍采用二维数学模型[6 ]或准三维模型. 如Jin 等人[7 ] 建立准三维耦合模式模拟浅水环流 ; [8 ] Irene 等人 在对数速度剖面的假定下 ,用渐近方法模拟波流共同作用下的三维泥沙输运. 这 些研究主要从计算方法上进行简化 ,减少了计算工作量 ,而对物理机理的考虑不够. 如 ,对河 口往复运动的水流,对数速度剖面近似在转流时刻附近显然是不合理的. Collins 等人[9 ] 分析 Loughor 河口及 Swansea 海湾的潮间带实测流速剖面认为 ,符合对数流速分布的数据所占比例 不足 40 % , 由此可见,对河口水流完全采用对数律来近似过于粗糙. 　　收稿 收修改稿 ( ) 　　 国家自然科学基金 批准号 、空气动力学流体力学国家科技部攀登项目和中国科学院力学研究所所长择优 基金资助项目 　　1) 匡翠萍. 长江口拦门沙冲淤及悬沙沉降规律研究和水流盐度泥沙数学模型. 南京水利科学研究院博士学位论文. 19