基于随机游动模型铁山港水交换数值模拟.docVIP

基于随机游动模型的铁山港水交换的数值模拟1，潘伟然,2*，骆智斌1,2，张国荣1,2，陶小琴1，郑艺妃1 (1.厦门大学海洋与地球学院， 2. 福建省海陆界面生态环境重点实验室(厦门大学)，福建 厦门361102) 摘要：采用三维浅海水动力模型粒子游动模型模拟了铁山港，分析了湾内水体时间的时空特征和物质特。研究结果表明铁山港海水交换主要受潮汐控制，层海水与湾外水交换率大，平均时间，层为25水体交换速度北暮-营盘附近海域关键词：铁山港；三维数值；随机游动；水体交换；铁山港地处北部湾东北部顶端，毗邻广东省英罗港。港湾鹿角状，湾口朝南呈喇叭状敞开与北部湾相口门宽约32。水域南北长约40，海湾面积3402，其中滩涂面积约占海湾总面积的50%。湾内掩护良好，无大的河流注入，两岸多开阔平坦的台地港资源。随着北部湾经济的快速发展，铁山港被列为重点港口开发区临海工业发展基地。港口的开建设对港域的水动力条件和水域环境容量产生一定的影响，每年排放入湾的生活污水和工业污水有增无减，水质受到了不同程度的。铁山港湾水交换物质输送扩散规律合理开发利用海湾资源。 海湾水交换能力的研究方法有箱式模型[1-2]保守物质盐度特征分布[3]放射性示踪物质浓度分布对流-扩散模型[4-7]标志质点的朗日粒子追踪模型。追踪模型流体质点，能够更直观更准确描述水交换特性，因此广泛用于水交换的研究中。早期的追踪模型忽略了扩散过程，在近岸弱流区往往低估海水的水交换能力。，，比踪更能真实反映水体的交换情况。Dimou [8]建立了水深平均的二维随机游动粒子踪模型，模拟充分混合的河口海湾地区的输运过程，并对浓度和驻留时间验证，结果符合良好。Bilgili等[9]应用随机游动模型模拟了Great 河口湾及其近岸海域各分区水体输运，阐释了各分区间水交换差异和特征。 目前国内外对于铁山港水交换的研究少，研究主要生物水体自净能力[10,11]。蒋昌波[12]二维对流扩散方程铁山港海水的半交换期，从物理角度建港前后水交换能力的差异。采用基于拉格朗日粒子追踪法的随机游动模型细致铁山港水体交换分布 模型及方法介绍 三维水动力模型和随机游动模型计算铁山湾水交换率。 .1 水动力数学模型 浅海三维动力模型数值模拟。模型在垂直方向上以固定水深分层，潮间带处理采用干/湿网格技术。在计算中采用Mellor-Yamada紊流模型[]和半隐式差分格式。模型以平均海平面作为x-y平面，z取垂直向上为正，通过Boussinesq 近似和静水压强假定得到不可压缩流体的控制方程如下 连续方程： （1） 水平方向动量方程： （2） Z方向动量方程： （3） 水位： （4） 压力计算： （5） 温度计算： （6） 盐度计算 , （7） 密度计算： （）式中，分别为方向的流速，为水位，为标准大气压，为温度，为盐度，为柯氏力参数，为压力，为海水密度标准值，为海水实际密度。分别为正压模式水平和垂直涡动粘性系数，分别为计算的水平与垂直方向的弥散系数。.2 随机游动模型 随机游动模型由Bugliarello[15]于1964年提出1971年Yang等[16]将表征物质的量粒子数转化成物质浓度该方法在流场通过释放大量粒子来物质的迁移扩散跟踪这些粒子的空间分布再反算成浓度场，进而了解物质迁移过程与特征。随机游动模型中粒子的对流输运过程由水动力学模型中求得的流确定扩散过程通过湍动扩散项粒子的位移表示如下 （） 其中，是某粒子在时刻的位置向量，则是在时刻的位置向量为质点在瞬时位置时的对流速度向量为随机游动时间步长为符合正态分布的随机数向量为水平方向上的湍动扩散系数张量。采用4阶龙格库塔法[1]求解中性粒子的 （） .2.1 水交换 本文采用两种指标来衡量铁山港的水体交换能力：一是稳态时间，示踪粒子数至原来10%所的时间，二是根据Takeoka[1]的定义表达平均驻留时间(mean residence time)流出海域的粒子数与初始粒子数之比为该海域水交换率。 （） 其中为平均驻留时间，为剩余函数指时刻域中的物质总量与初始时刻域内物质总量之比。 铁山港分湾顶湾中和湾口3个子域图1，分别标为，和，湾外记为。在随机游动方法中[1]，每个粒子代表恒定量的扩散质，假设初始时刻在研究域i中投放的粒子总数为，经过时刻从