测绘学报 ›› 2022, Vol. 51 ›› Issue (11): 2245-2254.doi: 10.11947/j.AGCS.2022.20210222

• 大地测量学与导航 •    下一篇

广义带限航空矢量重力确定大地水准面的两步积分法

黄谟涛1,2, 邓凯亮1, 吴太旗1, 王伟平2,3, 欧阳永忠2, 陈欣1, 王许1   

  1. 1. 海军研究院, 天津 300061;
    2. 自然资源部海洋环境探测技术与应用重点实验室, 广东 广州 510300;
    3. 国家海洋局南海调查技术中心, 广东 广州 510300
  • 收稿日期:2021-05-25 修回日期:2021-10-18 发布日期:2022-11-30
  • 通讯作者: 邓凯亮 E-mail:dengkailiang036@163.com
  • 作者简介:黄谟涛(1961—),男,教授,博士生导师,研究方向为海洋重力场测定理论方法及应用。 E-mail: dkl_hych@163.com
  • 基金资助:
    国家自然科学基金(42274015;41804011;42174013;41706111;41774021);国家重点研发计划(2016YFB0501704;2016YFC0303007);军科委173项目(2019-JCJQ-ZD-017)

A two-step integral method for geoid determination using generalized band-limited airborne vector gravity data

HUANG Motao1,2, DENG Kailiang1, WU Taiqi1, WANG Weiping2,3, OUYANG Yongzhong2, CHEN Xin1, WANG Xu1   

  1. 1. Naval Research Institute, Tianjin 300061, China;
    2. Key Laboratory of Marine Environmental Survey Technology and Application, Ministry of Natural Resources, Guangzhou 510300, China;
    3. South China Sea Marine Survey and Technology Center, State Oceanic Administration, Guangzhou 510300, China
  • Received:2021-05-25 Revised:2021-10-18 Published:2022-11-30
  • Supported by:
    The National Natural Science Foundation of China(Nos. 42274015;41804011;42174013;41706111;41774021);The National Key Research and Development Program of China(Nos. 2016YFB0501704;2016YFC0303007);The PLA Basic Research Program of China(173 Program)(No. 2019-JCJQ-ZD-017)

摘要: 依据物理大地测量学的广义带限水平边值问题理论,本文提出了基于带限航空矢量重力水平分量确定大地水准面的两步积分法。首先,依据广义带限水平边值问题解将带限航空矢量重力水平分量转化为飞行高度面上的带限扰动位;然后,利用带限Dirichlet边值问题解将飞行高度面上的带限扰动位向下延拓到海平面(海拔高程起算面,也称平均海面)上的带限扰动位;最后,依据Bruns公式将海平面上的扰动位转化为大地水准面高度。利用EGM2008重力位模型开展数值仿真计算试验,结果表明,广义带限水平边值问题解具有较好的低通滤波特性,能有效抑制观测高频噪声的影响;当带限航空矢量重力水平分量观测误差取3×10-5m/s2、测量飞行高度取6 km时,基于带限Dirichlet边值问题解的带限大地水准面计算精度优于3 cm,初步验证了两步积分法的计算稳定性和有效性。

关键词: 带限航空矢量重力水平分量, 大地水准面, 两步积分法, 广义带限水平边值问题, 带限Dirichlet边值问题

Abstract: According to the generalized band-limited horizon boundary value theory of physical geodesy survey, a two-step integral method is proposed to determine the geoid based on the horizontal component of band-limited airborne vector gravity data. Firstly,based on the solution of generalized band-limited horizon boundary value problem, we transform the horizontal component of band-limited airborne vector gravity data into a band-limited disturbing gravity potential at flight level. Secondly, the band-limited disturbing gravity potential is continued downward to the sea level by using the solution of the band-limited Dirichlet boundary value problem. Finally, according to the Bruns formula, the band-limited disturbing gravity potential at sea level is converted to the geoid. Simulation experiments using EGM2008 gravity field model show that the solution of generalized band-limited horizontal boundary value problem has a good low-pass filtering characteristic, which can weaken the effect of high-frequency noises in the observed value. The accuracy of the band-limited geoid obtained by the solution to the band-limited Dirichlet boundary problem is better than 3 cm when the error of the horizontal component of band-limited airborne vector gravity data is taken as 3 mGal and the flying height as 6 km. It indicates that the two-step integral method has good stability and effectiveness.

Key words: horizontal component of band-limited airborne vector gravity data, geoid, two-step integral method, generalized band-limited horizon boundary value problem, band-limited Dirichlet boundary value problem

中图分类号: