本文在二维潮波方程的基础上,采用变分方法,给出其伴随方程。并基于Arakawa C网格给出了潮波方程及其伴随方程的差分格式,建立了二维潮波模式及其伴随模式。在此基础上,采用渤海沿岸19个验潮站的调和常数资料,实现了渤海M2、S2、K1、O1四个主要分潮的伴随同化数值模拟。所得四个主要分潮调和常数的精度为:振幅绝均差分别为2.04cm、0.76cm、1.38cm、0.55cm,迟角绝均差分别为1.81°、1.29°、2.48°、1.73°,均方根偏差分别为4.42cm、1.41cm、2.81cm、1.23cm,拟合程度分别为99.70%、99.65%、98.91%、99.58%,比前人结果有所改善。利用伴随同化方法,本文对潮波模式中的水平涡动粘性系数的优化进行了探讨。设计了三组数值实验用以研究水平涡动粘性系数的优化效果:实验1只优化底摩擦系数,实验2只优化水平涡动粘性系数,实验3同时优化底摩擦系数和水平涡动粘性系数。这三组实验得到的M2分潮调和常数的模拟值和观测值之间的偏差值分别为:实验1,振幅绝均差4.83cm,迟角绝均差6.36°;实验2,振幅绝均差7.14cm,迟角绝均差7.80°;实验3,振幅绝均差2.04cm,迟角绝均差1.81°。实验结果表明,同时优化底摩擦系数和水平涡动粘性系数对提高模式结果的精度有重要作用。对于S2、K1、O1分潮,也得到了同样的结论。