周 浩 钟 波 罗志才 张 坤

(1)武汉大学测绘学院，武汉 430079

2)地球空间环境与大地测量教育部重点实验室，武汉430079)

OpenMP并行算法在卫星重力场模型反演中的应用*

周 浩1)钟 波1，2)罗志才1，2)张 坤1)

(1)武汉大学测绘学院，武汉 430079

2)地球空间环境与大地测量教育部重点实验室，武汉430079)

利用卫星重力数据反演地球重力场需要解决重力场模型的高效计算问题。分析了最小二乘直接法求解重力场模型涉及的密集型计算任务，基于OpenMP实现了卫星重力场模型直接求解的并行算法。利用30天、5秒采样间隔的沿轨扰动位T和径向扰动重力梯度Trr数据，分别反演了60阶次的卫星重力场模型，计算结果表明，OpenMP并行算法能够有效提高直接法求解卫星重力场模型的计算效率，并具有很好的稳定性。

并行算法;重力场模型;卫星重力;重力梯度;OpenMP

1 引言

卫星重力测量技术为解决全球高覆盖率、高空间分辨率和高时间重复率的重力测量开辟了新途径，特别是新一代卫星重力计划的成功实施，已使得地球重力场模型的中长波部分精度提高了约2个量级［1-3］。

新一代卫星重力测量数据具有高采样率、全球观测的特点。利用卫星重力观测值反演地球重力场将面临海量数据的处理问题，特别是对于高阶重力场模型的反演，其计算量巨大，计算过程极为耗时。为了利用海量卫星重力观测值快速求解高精度的地球重力场模型，许多学者一方面是对现有重力场恢复方法进行改进，提出了一些快速解算法［4-6］;另一方面就是采用并行算法，即基于并行计算机，编写合适的并行算法来提高计算速度［7-9］。……