河口几何特征对河口区盐淡水混合的影响的研究.pdfVIP

中国环境与生态水力学2008 河口几何特征对河口区盐淡水 混合的影响研究 袁丽蓉程香菊 赵文谦 (中山大学工学院海洋工程研究中心，广东广州 510275) 摘 要 为了研究河口几何特征对河口区盐淡水混合的影响，本文建立非结构化 的三维数值模型，针对断面宽度接指数规律变化，横断面为抛物线，具有不同扩 展率和不同坡度变化的河口，研究河口岸线变化和地形变化对河口流场和盐度场 的影响。该模型的特点是采用非结构化的网格，可以很好的模拟岸线的曲折变 化；采用欧拉一拉格朗日方法计算动量方程和物质输运方程的对流项，在速度很 大、方向不断改变的流场中仍然能够以较大的时间步长进行计算。本文采用理想 化的岸线和地形，得到一些普遍的规律，可以更好的理解河口的动力过程，认识 河口环流和成水入侵现象，为河口整治、航道改造等改变河口岸线和地形的工程 提供设计依据，减少因人为原因加剧成潮入侵。 关键字 河口几何特征 河口环流 盐淡水混合 非结构化三维数值模型 1 前言 河口是淡水和盐水交汇混合的区域，也是河流动力和海洋动力相互作用的场所，河口 的流场和盐度场由特定的岸线和地形以及径流、潮汐、风力、科氏力等因素决定，形成一 个复杂的系统。发生在河口的咸水人侵既是一个社会问题，也是一个科学问题。很多人海 的河流承担着农业灌溉和城市饮水的功能，但是近些年人类活动加剧了咸水入侵，甚至造 成有的取水口现在已经无淡水可取。加强河口盐淡水混合问题的研究，揭示盐水入侵机理 具有重要的科学价值和现实意义。 Bowden等通过实验、现场观测和概念模型对密度环流和盐淡水混合做了比较全面的研 基金项目：国家自然科学基金，环境保护公益性行业科研专项(200709049) 作者简介：袁丽蓉(1976·)，女，湖南邵阳人。讲师，主要研究方向：环境水力学，近岸环境动力 学，E-mail：yuanlr@mail．sysu．edu．cn · 61· 口垂向断面上的各种紊动量，经过分析认为该模型用于河口的潮流计算还有很多需要改进 的地方。龚政[53和罗小峰‘61对长江口的盐度场做了数值模拟和分析工作。朱建荣等口]对河 口环流和盐水入侵做了理论分析和数值模拟。包芸和任杰[8]运用三维斜压模型对珠江口伶 仃洋进行了数值模拟，模拟结果反映出盐度高度层化的整体分布。 通过分析认为，对河口盐淡水混合问题进行数值模拟的难点在于河口介于河流与海 洋之间，空间尺度变化很大，模型需要在网格划分上具有很大的灵活性，才能同时适应 河流和海洋的空间尺度特征；盐度不是被动输运物质，盐度场本身也是速度场的驱动 力，必须与速度场耦合求解，而盐度驱动与其它驱动力比较其量级小，变化速度慢，需 要对长时间序列做计算才能反映出河口盐淡水混合的实际情况，模型必须具有稳定性和 高效性。通过对现有模型的综合评价，本文选择ELCIRC模型进行数值试验研究。EL— CIRC模型是非结构化的三维斜压模型，适用于研究空间上跨河流和海洋尺度的河口问 题。它采用了控制体积法、差分法和拉格朗日方法求解浅水方程，其中包含了河流、海 洋和大气驱动等物理过程。该数值模型是低阶的，但是保证了体积守恒和计算的稳定性 高效性，并且加入了潮滩干湿边界条件。模型最初是为模拟Columbia河口而开发的， 现在已经对近岸河口问题作了广泛的测试，也正在被越来越多的应用于世界各地的河口 和近岸研究‘。．10]。 2数值模型 2．1 模型介缁 度、盐度守恒方程)来得到自由面、三维速度场和温度、盐度。下面列出具体的控制 方程： 差+寰+装=o (1) 争毫腮池+南牌地一。 ㈣ c3， 器一扣一毫{咖一匆，+导卜盖r如爱dz+毫(K。爰)