异构并行算法快速构建全球扰动重力梯度全张量图

摘要

扰动重力梯度是扰动重力位的二阶导数,相对于其他重力场元素能更多地反映变化的不规则地球产生的高频信息。在使用高阶次球谐系数模型获取大范围高分辨率的扰动重力梯度数据时,存在重复运算多、计算效率低下、耗时较长的问题。针对该问题,推导了简化计算公式,将中间变量提取出来作为全局参数和局部参数单独进行计算、存储,从而有效减少重复运算;并在简化公式的基础上,提出了扰动重力梯度张量快速异构并行算法,利用CUDA(compute unified device architecture)实现了梯度全张量在GPU端的并行计算。根据T_(xx)、T_(yy)、T_(zz)三个分量满足Laplace条件验证了算法可靠性,并与传统串行算法进行了计算效率对比,实验结果表明,相较于串行算法,所提算法可减少90%以上计算耗时,可将计算效率提高60倍以上。最后利用该算法基于2 190阶EIGEN-6C4模型快速构建了5′×5′分辨率的全球扰动重力梯度全张量图,计算结果显示了扰动重力梯度同地形、地球质量分布变化的相关性及其在全球范围内的数值特征。