海域潮汐模型开边界反演优化方法初探——以K1分潮振幅为例 Optimization inversion method for tidal model under open boundary condition A case study of K1 constituent amplitude.pdfVIP

第20卷第1期 水科学进展 V01．20。No．1 2009年1月 INWATERSClENCE Jan．．姗 ADVANCES 海域潮汐模型开边界反演优化方法初探 ——以K1分潮振幅为例 李明昌1,2，梁书秀1，孙昭晨1 (1．大连理工大学海岸和近海工程国家重点实验室，辽宁大连116024； 2．交通部天津水运工程科学研究院水路交通环境保护技术实验室，天津3410456) 摘要t将数据驱动模型理论和海域潮汐模型有机地结合在一起．建立海域潮汐模型开边界反演新方法，并纳入潮流 资料作为反演工作的依据。该方法通过海域潮汐模型设计工况的计算构建内部观测点的解集空间，并归纳建立内 部观测点潮汐(或潮流)同开边界潮汐之间的关系，将观测点潮汐和潮流实测资料带人关系之中反演开边界条件。以 辽东湾海域潮汐模型KJ分潮振幅开边界条件为例验证反演方法，表明该开边界反演方法是有效的。 关键词t开边界条件；反演；数据驱动模型；潮流；潮汐模型 中图分类号iP731．23 文献标识码：A 文章编号：1001-6791(2009)01-0129-06 开边界条件的确定是海域潮汐模型验证的主要工作内容。传统的模型验证过程中，开边界的初始值可由边 界附近潮汐测站资料插值方法给出，或由更大范围数学模型提供。在此基础上采用试错法，最终在一定精度范 围内得到满足计算域内实测值的边界条件。在实际工程中，潮汐资料可能是相当匮乏的，或者存在数据缺失和 时间上的不对应性；而大范围数值模型提供的开边界条件对于浅水海域精度较低。因此，发展系统的、优化的 开边界反演方法极具现实意义。 针对海域潮汐模型的开边界反演，国内外学者进行了大量的数值研究【1_8I，其中最普遍采用的方法是伴随 法。尽管伴随法在优化模型初始场、参数估计、开边界条件等方面取得了良好的研究成果，但是由于伴随方法 需要迭代计算正向的潮汐方程和逆向的伴随方程，其计算过程复杂并且耗费大量的机时[9|，虽然可以通过良好 的猜测初值来减少迭代步数【5J，但是实际工程中的资料并不是很充分，因此这就给实际工程研究的模型验证工 作带来了很大的不便。目前伴随方法反演所依赖的海域内部观测资料均是潮位资料，尚未应用潮流资料[10]。 而潮流作为潮波的水平运动形式，反映了海域流场的基本形态，是实际工程模型验证工作的重点。因此在开边 界条件的反演工作中，潮流资料的引入非常有意义。 本文以数据驱动模型IIIj理论的基本思想为基础，将其同海域潮汐模型有机地结合，提出开边界条件反演 的新方法，并将海域内部潮流观测资料纳入开边界反演的系统方法中。 1潮汐模型及计算海域 1．1潮汐模型 本文海域潮汐模型采用二维POM模型121(PrincetonOcean Model)，是由普林斯顿大气和海洋科学实验室研 制发展起来的一类海洋数值计算模型，目前POM模型在世界范围内的海洋数值研究中得到广泛的应用。 收稿日期：2007．12．13 基金项目：国家自然科学基金资助项目 作者简介：李明昌(1977一)，男。辽宁大连人，工程师，博士，主要从事海洋环境水动力学研究。 E-mail：Imcsql997@163．corn 通讯作者：孙昭晨。E-mail：sunzc@dlut．edu．cn 万方数据 130 水科学进展 第20卷 1．2模型计算海区及模型参数设置 辽东湾是本文的研究海域，位于渤海北部，面积 18300 km2，海底地形平坦，向渤海中部缓倾，平均坡度在 l／2000．1／2500，顶部属平原淤泥岸。海域受渤海黄海的 潮流影响，流态比较复杂，大面积滩涂的存在对潮流场有 着较大的影响。 图1是计算海区及沿潮站位置图，图中秦皇岛(西侧) 和长兴岛(东侧)连线是海域开边界