卫星测高重力与重力反演地形

卫星测高的原理

卫星测高是以装载雷达高度计的卫星为平台，主动向海面发射微波信号，接受经海表面反射的微波信号以获取海面高度、有效波高和后向散射等信息。

测高卫星在预先设定的轨道上围绕地球飞行可提供全天候、近实时、海量的卫星观测数据，其星下点轨迹具有往复性和高密度的特点。通过对卫星轨迹升弧段和降弧段的交叉点不符值进行处理，并采用相应的格网化内插方法可获得格网化大地水准面数据和垂线偏差数据，进而利甩逆Stokes公式、最小二乘配置以及逆Vening-Meinesz公式等方法可恢复海洋重力场。

本文采用的 1'×1'卫星测高重力异常 GravAlti_WHU 是联合 Geosat、ERS-1、ERS-2、Envisat、T/P、Jason-1、CryoSat-2 和 SARAL/AltiKa 等多颗卫星波形资料获取的垂线偏差信息反演的,如图1所示。数据处理过程中改进了波形重跟踪算法，选定了较优的沿轨重采样频率,在北太平洋海域、墨西哥湾海域与 NGDC(National Geophysics Data Center)船测数据之差的均方根约为 4 mGal(1 mGal=10-5 m/s),模型精度与 SIO V23.1 重力异常模型相当,优于 DTU10 重力异常模型。

图 1 全球海洋卫星测高重力异常模型 Grav_Alti_WHU

海底地形反演

根据海洋重力场信．息反演海底地形的基本原理可以描述为：依据地壳均衡理论，海底地形起伏及其均衡补偿物质的密度分布异常将会引起海面重力异常和重力异常垂直梯度等重力信息发生相应变化，因此可以说重力数据与海底地形在一定程度上具有一一对应的相关关系，基于Parker（1972）异常扰动位计算公式和Watts（1978）三个板块模型可以获得它们之间的地球物理关系。同时也可分析大量的重力和海深数据的统计特征，采用一定的数学方法构建它们之间的数学模型。利用重力异常和重力异常垂直梯度数据反演海底地形作为研究对象，围绕其理论、方法及其应用展开研究，以期为提高海洋开发能力，扩大海洋开发领域提供理论与技术支持。

（1）重力地质方法
重力地质方法（Gravity-Geologic Method，简称GGM）是一种常用的海底地形反演技术，最初被应用于冰河时期漂移沉积物下的基岩高度测量，由于陆地上地质层的密度变化大，使得该方法在陆地上的应用效果不佳。而在海底地形的探索方面，洋壳密度和海水密度的差异变化较小，通常认为海水内部密度不随深度变化。从而重力地质法（GGM）在反演海洋深度方面具有较高的可行性。
重力地质方法的实施关键在于密度差异常数的确定，一种方法为采用移去恢复技术建立海深模型，并通过数学方法构建模型海深值与船舶测量检核点之间的相关系数和均方差（ＲＭＳ）关系，利用密度差异迭代的方式获取最优的密度差异常数；另一种方法为采用频率域向下延拓方法，以参考海深面重力异常与海面重力异常之比作为密度差异常数。如今，重力地质方法依靠少数船测水深和卫星测高重力异常，便可间接获取相关海域海洋深度，与其他经典方法相比，该方法具有模型简单，易于计算的优点。
（2）导纳函数方法
重力导纳表征将海底地形转换为重力异常的能力。基于对Parker公式、Watts板块模型以及地壳均衡理论的研究，利用最新的全球地壳模型Crust1.0、“理论导纳”函数、“观测导纳”函数等相关知识，可以获取有关的海洋地球物理信息。在重力数据和海深数据频谱特征的统计分析基础上，采用相应的滤波技术对信号向下延拓过程中出现噪声放大的病态问题进行处理。可以以重力异常和重力异常垂直梯度为输入数据，分别建立不顾及均衡情况、考虑Ａｉｒｙ均衡、或者顾及地壳挠曲均衡三种不同模式的导纳函数模型。

导纳函数方法从海底地形起伏在海面对应区域产生扰动位的经典位势理论出发，经过傅里叶变换、极坐标积分变换、泰勒级数等一系列推导得到频率域内海底地形起伏计算海面重力异常的级数展开式，然后利用弹性板挠曲理论以及均衡模型最终建立了反演海底地形的导纳函数模型。该方法充分考虑了海底地形的地球物理环境，算法理论严密，但是解算相对复杂。
（3）Ｓ＆Ｓ方法
Ｓ＆Ｓ方法主要由Smith和Sandwell在１９９４年提出，相对于理论完备的导纳函数方法，在对大量的重力和海深数据的统计特征分析发现，重力异常或者重力异常垂直梯度经滤波并向下延拓后的结果与海底地形存在良好的线性关系。基于此，利用数学建模相关理论可以构建相关波段重力异常或者重力异常垂直梯度与海底地形的数学函数关系，获取重力数据与海底地形的比例系数，进而反演出相应波段的海底地形结果。
相对于重力地质方法和导纳函数方法而言，Ｓ＆Ｓ方法基于对大量的重力和海深数据统计分析，进而对数据进行数学建模。Ｓ＆Ｓ方法虽然不如导纳函数具有严密的理论推导，但是方法简单，具有较强的可操作性，同时比例系数的结果也在一定程度上反映了海底地形的某些地球物理特征。
（4）最小二乘配置方法
最小二乘配置是大地测量学中的经典方法，利用最小二乘配置方法推估海底地形要求事先知道海深与重力数据之间的自协方差函数和互协方差函数。根据重力场理论，可以采用Inverse模型、Gauss模型、Moritz模型拟合重力数据的协方差函数，但是海深与重力的协方差函数计算困难。此时，可以通过频率域功率谱密度函数（Power Spectral Density，简称PSD）模拟自协和互协方差函数，进而得到重力－海深在频率域内的转换函数，然后依据十次多项式拟合转换函数获得平滑的转换函数，最终解算相应波段的海深值。

国内利用卫星测高反演海底地形的成果

武汉大学李建成院士课题组是国内较早研究卫星测高技术应用研究的课题组之一。2014 年，课题组使用 1'×1'的 SIO V20. 1 重力异常垂直梯度数据，采用回归分析方法，联合
NGDC 发布的船测水深数据，构建了 75°S～70°N 范围 1'×1'的海底地形模型 BAT_VGG，其中反演区域以外的海底、陆地以及水 深 浅 于 100 m 海 域 的 地 形 数 据 来 源 于 SIOV15. 1。在印度洋南部海域和西北太平洋海域进行模型质 量 评估 考 察 表 明，BAT _VGG 模 型 精度 优 于ETOPO1 模型、GEBCO_08 模型和 DTU10 模型，与 SIOV15. 1
相当。2020 年，李建成院士课题组基于新构建的全球卫星测高重力异常模型 Grav_Alti_WHU，使用回归分析方法，联合水深测量资料，建立了75°S ～70°N 范围 '×1'的海底地形模型 BAT_WHU2020。其中反演区域以外的海底、陆地以及水深浅于 100 m海域的地形数据来源于 SIO V19. 1。最终模型分析评估结果显示，BAT_WHU2020 模型精度优于 ETOPO1模型、GEBCO_08 模型和 DTU10 模型; 较之前发布的BAT_VGG 模型，BAT_WHU2020 模型精度提升了约30%。

图 2 全球 １′×１′ 海底地形模型 BAT_WHU2020

钱学森实验室雷达团队和中国地质大学 ( 北京) 万晓云副教授团队联合构建了全球深海海底地形模型，该海底地形模型将 GGM 方法和经典频域法相结合，海深浅于 1 000 m 和陆地区域地形数据用ETOPO1 模型补充。

信息工程大学李姗姗教授团队长期致力于利用卫星测高技术构建海底地形模型研究。2020 年，团队利用 SIO V29. 1重力异常和SIO V29. 1 重力异常垂直梯度，融合多源水深测量结果，构建了 STO_IEU2020 全球海底地形模型，分别在大西洋海域和印度洋海域等海域，选择了共计 5 块海域作为海底地形模型质量评估区，评估结果表明，±100 m检核差值范围内，STO_IEU2020 模型对应的检核点数量分别占 78. 57%、87. 87%、72. 27%、92. 47%和 93. 25%; 检核结果平均相对精度优于 6%。STO_IEU2020 模型与 SIO V20. 1 模型贴切与吻合程度高。
用必应搜索国内的反演海底地形模型，BAT_WHU2020模型和STO_IEU2020 模型应该都没有提供下载连接。希望相关模型能够像ETOPO1和GEBCO那样作为公共产品开放给全球海洋工作者使用。

胡敏章,张胜军,金涛勇,等.新一代全球海底地形模型BAT_WHU2020[J].测绘学报, 2020, 49(8):16.DOI:10.11947/j.AGCS.2020.20190526.