论文部分内容阅读
潮汐潮流运动是海洋中的基本运动之一,它是动力海洋学研究的重要组成部分,对它的研究直接影响着波浪、风暴潮、环流、水团等其他海洋现象的研究,在大陆架浅海海洋中,对潮汐潮流的研究更具重要性.伴随同化方法是建立在严格的数学基础之上的一种有效的四维变分同化技术,它将变分原理与最优控制理论相结合,以海洋动力模型作为约束条件,通过同化可观测到的海洋要素优化控制变量,达到反演无法观测到的海洋要素的目的.潮汐潮流数值模拟中最大的难点莫过于开边界条件和底摩擦系数的确定,在数学物理反问题的框架下,伴随同化方法可以把开边界条件和底摩擦系数的确定变成依赖于区域内部观测资料的数值迭代过程,实现了调试控制参数的自动化和最优化.本文对伴随同化方法在潮汐潮流数值模拟中的应用进行了深入研究.本文建立了一个球坐标系下的二维非线性潮波伴随同化模型,并基于渤、黄、东海M2分潮的数值模拟,对空间分布底摩擦系数处理方法进行了伴随法反演研究.孪生实验结果表明,空间分布底摩擦系数与伴随同化方法相结合,可以成功反演出给定的空间结构非常复杂的底摩擦系数分布,这是传统的处理方法所不能做到的.底摩擦效应由海底地形决定,本文尝试了根据海底地形的空间分布特征选取独立底摩擦系数的做法;数值实验结果表明,与均匀选取相比,根据底地形的空间分布特征选取独立参数,能在减少独立变量个数的情况下提高模拟精度.通过同化T/P高度计资料,本文数值模拟了渤、黄、东海的M2潮波;实验结果表明,与传统的底摩擦系数处理方法相比,空间分布的底摩擦系数能有效提高模拟精度,实验结果准确体现了渤、黄、东海M2潮波的分布特征.使用上述二维非线性潮波伴随同化模型,考虑分潮间的非线性相互作用,本文通过同化T/P高度计资料,优化模式开边界条件及空间分布的底摩擦系数,数值模拟了渤、黄、东海的M2、S2、K1、O1四个主要分潮,模拟结果与验潮站资料(独立资料)之间的绝均差分别为(5 .7cm ,5.2)、(3 .2cm ,6.9)、( 2.6cm ,6.3)和( 2.1cm ,8.5).实验所得同潮图准确体现了渤、黄、东海四个主要分潮的潮位分布特征,半日分潮共有三个无潮点和一个退化型无潮点,全日分潮有三个无潮点,分别位于渤海海峡、南黄海和对马海峡.基于内外模态分离技术,本文建立了一个三维正压非线性潮汐潮流模型,外模态采用ADI方法离散,时间步长不受CFL条件的限制;内模态的时间步长则可以远大于外模态时间步长,从而显著提高计算效率.将上述模型作为正向模型,本文建立了其伴随模型.此三维伴随模式建立之后,运算过程中开边界条件、底摩擦系数和垂向涡动粘性系数等控制参数的选取、数值模式和观测结果的有机结合等问题就可迎刃而解.与前人的有关三维潮汐潮流伴随同化模型的工作不同,本模式将底应力项设置为底层流速的函数,从而在物理上更加合理.在数学物理反问题的框架下,本文设计孪生实验研究了三维正压潮汐潮流模型中开边界条件、空间分布底摩擦系数和垂向混合涡动粘性系数的伴随法反演问题,以验证伴随模型的有效性和正确性.同时,本文分别讨论了参数的不同类型分布、参数初值、观测误差以及观测数量对参数估计反问题的影响,对反问题适定性进行了初步研究.