宗敬文，李厚朴，纪 兵，欧阳永忠

1. 海军工程大学导航工程教研室，湖北 武汉 430033； 2. 自然资源部海洋环境探测技术与应用重点实验室，广东 广州 510300

海洋重力场数据的获取手段有卫星测高反演、船载及航空重力测量，随着人造地球轨道卫星技术的发展，卫星测高技术促进了大地测量学、地球物理学和海洋学等学科的交叉发展[1]，同时与传统的海洋观测手段相比，能够周期性获取除极地以外的高分辨率、全天候、长时间序列的全球海洋观测数据[2-5]。截至目前，世界上各个国家共计发射测高卫星18颗，其中包括美国海军的Geosat系列、美国宇航局和法国空间研究中心联合发射的T/P系列[6]、欧空局的ERS系列[7]及中国的HY-2A[8]。卫星测高数据遍布全球85°S—85°N的海洋区域，分辨率达1′×1′，精度达1～2 cm，利用高精度、高分辨率卫星测高数据反演重力异常精度已达到10 mGal(1 Gal=10-2m/s2)以内。文献[9]使用逆Stokes法反演中国近海30′×30′海洋重力异常，精度达到3.5 mGal。文献[10]使用Laplace方程的垂线偏差法反演全球海域的重力异常，反演的全球海域1′×1′重力场模型与NGDC船测数据的检核精度达到4～8 mGal。文献[11—14]在卫星测高数据融合和精细处理方面也作出了大量的贡献。

文献[15]利用Stokes公式将卫星测高数据转换为重力异常。文献[16]基于垂线偏差对地形信息敏感的特点，根据边值理论由重力与地形数据确定格网垂线偏差模型。文献[17]提出了利用垂线偏差来反演重力异常的方法，并系统地给出了垂线偏差、扰动重力和重力异常之间的关系表达式，但在利用大地水准面高和垂线偏差反演重力异常实际计算中，计算点及其附近区域到计算点的理论距离近似于零，会导致逆Stokes公式和逆Vening-Meinesz公式中的积分函数产生奇异。……