本文用流体动力学、热力学和κ—ε两方程湍流闭合模式相结合的方法对渤海温度混合和层化的演变过程进行三维数值模拟。 首先在流体动力学和热动力学的基础上，应用雷诺时均方程法，详细推导得到海水在湍流状态下的三维原始方程。然后引入湍流闭合模式来使原始方程封闭。最开始用零方程湍流闭合模式，其计算结果和多年统计结果有较大差距。后来应用了κ—ε两方程湍流闭合模式，该模式合理地表达了由湍流产生的混合效应。用这种流体动力学、热动力学和κ—ε两方程湍流闭合模式相结合的方法，较好地再现了渤海温度场从混合至层化到再混合的季节性变化过程，以及温跃层的季节性变化。三月中旬，渤海层化首先出现在莱州湾和渤海湾，四月和五月层化不断扩大范围和增加强度，到七月中旬层化强度达到最大，随后层化开始缓慢减弱，到九月层化减弱加剧，范围也减小，到十月底，层化完全消失，渤海全区水温垂直分布呈均匀状态。模拟结果与多年统计结果符合较好。且完整、全面地反映了整个渤海温跃层空间分布特征，及其随时间演变的细致过程。 这种方法具有适普性，可用于其它陆架海域混合和层化的研究。