气象资料变分同化高效优化算法研究及并行实现

【摘要】：数值天气预报模式对初始条件十分敏感，因此初值的优劣是决定数值预报效果的关键之一。利用气象观测资料、通过资料同化技术形成数值预报模式的初值，是数值天气预报研究领域的一个重要方向。如何利用卫星、雷达等现代观测手段所产生的非常规观测资料是资料同化技术面临的一个重大挑战。 变分资料同化技术是解决上述问题的有效方法，它将资料同化问题描述为一个无约束(三维变分同化)或一个以动力学模式为约束(四维变分同化)的极小化问题。 变分同化问题的核心算法是最优化算法。本文适应“中国气象数值预报系统技术创新研究”的需要，实现了多种针对大规模问题的无约束优化算法。详细分析了每个算法的适用范围，并通过大量数值实验数据对最速下降算法、非线性和线性共轭梯度法、L-BFGS方法和离散截断牛顿法的收敛速度、计算效率和所需存储空间进行了详细比较，指出了各种算法的优劣。结合业务化变分同化系统的实际应用特点，提供了最优化算法的选择方案。研究成果在有我国自主知识产权的GRAPES三维变分同化系统中得到了应用。 本文还针对线性共轭梯度法和L-BFGS法的并行计算进行了研究，设计了这两个算法的并行算法及实现方案，分析了它们的并行效率，并给出了测试结果。