重力异常向上延拓严密改化模型及向下延拓应用

doi:10.11947/j.AGCS.2022.20200547

测绘学报 ›› 2022, Vol. 51 ›› Issue (1): 41-52.doi: 10.11947/j.AGCS.2022.20200547

重力异常向上延拓严密改化模型及向下延拓应用

黄谟涛^1,2, 邓凯亮¹, 吴太旗¹, 欧阳永忠², 陈欣¹, 刘敏³, 王许¹

1. 海军研究院, 天津 300061;
2. 自然资源部海洋环境探测技术与应用重点实验室, 广东广州 510300;
3. 91001部队, 北京 100830

收稿日期:2020-12-11 修回日期:2021-09-24 发布日期:2022-02-15
通讯作者: 邓凯亮 E-mail:dengkailiang036@163.com
作者简介:黄谟涛(1961-),男,博士生导师,研究方向为海洋重力场测定理论方法及应用。E-mail:dengkailiang036@163.com
基金资助:
国家自然科学基金（41804011；41706111；41774021）；国家重点研发计划（2016YFB0501704，2016YFC0303007）

Rigorous modification model of upward continuation and its applications on the downward continuation of gravity anomaly

HUANG Motao^1,2, DENG Kailiang¹, WU Taiqi¹, OUYANG Yongzhong², CHEN Xin¹, LIU Min³, WANG Xu¹

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. Troops 91001, Beijing 100830, China

Received:2020-12-11 Revised:2021-09-24 Published:2022-02-15
Supported by:
The National Natural Science Foundation of China(Nos. 41804011;41706111;41774021);The National Key Research and Development Program of China(Nos. 2016YFB0501704;2016YFC0303007)

摘要/Abstract

摘要： 重力异常向上延拓全球积分模型在航空重力测量数据质量评估和向下延拓迭代计算等领域具有广泛的应用。为了消除积分核函数奇异性影响，需要对该模型进行基于积分恒等式的移去-恢复转换及全球积分域的分区改化处理。在此过程中，传统改化处理方法往往忽略了全球积分过渡到局域积分引起的积分恒等式偏差影响，从而导致不必要的计算模型误差，最终影响向上延拓计算结果的可靠性，甚至影响向下延拓迭代解算结果的稳定性。针对此问题，本文开展了重力异常向上延拓积分模型改化及向下延拓应用分析研究，依据实测数据保障条件和积分恒等式适用条件要求，导出了重力异常向上延拓积分模型的分步改化公式，提出了补偿传统改化模型缺陷的修正公式，并将最终的严密改化模型应用于重力异常向下延拓迭代解算。使用超高阶地球位模型EGM2008作为标准位场开展数值计算检验，分别对重力异常向上延拓分步改化模型的计算精度及在向下延拓迭代解算中的应用效果进行了检核评估，验证了采用严密改化模型的必要性和有效性。

关键词: 重力向上向下延拓, 模型改化, 奇异积分, 积分恒等式, 精度评估

Abstract: The global integral model of upward continuation for gravity anomaly is often used in the data quality assessment for airborne gravity measurement and the iterative solution of downward continuation. In order to eliminate the singularity of the integral kernel function, we have to make a remove-restore transformation to the integral model based on an equation of integral identity. And at the same time, the global integral domain has to be divided into two parts of near and far zones. During this processing, the deviation of identical formula caused by the transformation of integral domain from global to local integration is often ignored in the traditional modified method, which leads to unnecessary model errors and ultimately affects the reliability of upward continuation solution, and even affects the stability of iterative solution of downward continuation. In response to this problem, some analysis and research have been made on the integral model modification of upward continuation and its application of downward continuation for gravity anomaly in this paper. Based on the measured data situation and the applied condition of the integral identity, a set of step-by-step modified formulas of the upward continuation integral model are derived, and a correction formula is proposed to compensate the omission of the traditional modified model. Finally, the rigorous modified model is applied to the iterative solution of downward continuation for gravity anomaly. The ultra-high-degree geopotential model EGM2008 is used as a standard potential field to carry out a set of numerical calculation test, and the calculated accuracy of step-by-step modified models of the upward continuation and the applied effect in the iterative solution of downward continuation for gravity anomaly are estimated and evaluated separately. It proves the necessity and effectiveness of adopting the rigorous modified model.

Key words: upward and downward continuation of gravity anomaly, model modification, singular integration, integral identity, accuracy evaluation

中图分类号:

P223

黄谟涛, 邓凯亮, 吴太旗, 欧阳永忠, 陈欣, 刘敏, 王许. 重力异常向上延拓严密改化模型及向下延拓应用[J]. 测绘学报, 2022, 51(1): 41-52.

HUANG Motao, DENG Kailiang, WU Taiqi, OUYANG Yongzhong, CHEN Xin, LIU Min, WANG Xu. Rigorous modification model of upward continuation and its applications on the downward continuation of gravity anomaly[J]. Acta Geodaetica et Cartographica Sinica, 2022, 51(1): 41-52.

参考文献

[1] HEISKANEN W A, MORITZ H. Physical geodesy[M]. San Francisco:W. H. Freeman and Company,1967.
[2] MORITZ H. Advanced physical geodesy[M]. England:Abacus Press, 1980.
[3] 李建成,陈俊勇,宁津生,等. 地球重力场逼近理论与中国2000似大地水准面的确定[M]. 武汉:武汉大学出版社,2003. LI Jiancheng, CHEN Junyong, NING Jinsheng, et al. Theory of the Earth's gravity field approximation and determination of China quasi-geoid 2000[M]. Wuhan:Wuhan University Press, 2003.
[4] 黄谟涛,翟国君,管铮,等. 海洋重力场测定及其应用[M]. 北京:测绘出版社, 2005. HUANG Motao, ZHAI Guojun, GUAN Zheng, et al. Determination and application of marine gravity field[M]. Beijing:Surveying and Mapping Press, 2005.
[5] SANSÒ F, SIDERIS M. Geoid determination:theory and methods[M]. Berlin:Springer,2013.
[6] HIRVONEN R A, MORITZ H. Practical computation of gravity at high altitudes[M]. Ohio:Ohio State University, 1963.
[7] CRUZ J Y, LASKOWSKI P. Upward continuation of surface gravity anomalies[M]. Ohio:Ohio State University, 1984.
[8] 翟振和,孙中苗,李迎春,等. 航空重力测量在近海区域的精度评估与分析[J]. 测绘学报, 2015, 44(1):1-5. ZHAI Zhenhe, SUN Zhongmiao, LI Yingchun. The accuracy evaluation and analysis of airborne gravimetry in coastal area[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(1):1-5.
[9] HUANG Motao,DENG Kailiang,WU Taiqi,et al.Research and evaluation on key technological indicators for airborne and shipborne gravimetry[J].Journal of Geodesy and Geoinformation Science,2019,2(3):44-54.
[10] MORITZ H. The computation of the external gravity field and the geodetic boundary-value problem[M]//Orlin H ed. Gravity anomalies:unsurveyed areas. Washington D C:American Geophysical Union, 1966:127-136.
[11] 王兴涛,石磐,朱非洲. 航空重力测量数据向下延拓的正则化算法及其谱分解[J].测绘学报,2004,33(1):33-38. WANG Xingtao, SHI Pan, ZHU Feizhou. Regularization methods and spectral decomposition for the downward continuation of airborne gravity data[J]. Acta Geodaetica et Cartographica Sinica, 2004, 33(1):33-38.
[12] 徐世浙. 位场延拓的积分-迭代法[J]. 地球物理学报, 2006, 49(4):1176-1182. XU Shizhe. The integral-iteration method for continuation of potential fields[J]. Chinese Journal of Geophysics, 2006, 49(4):1176-1182.
[13] 黄谟涛,刘敏,邓凯亮, 等. 基于向上延拓的航空重力向下解析延拓解[J]. 地球物理学报, 2018, 61(12):4746-4757. DOI:10.6038/cjg2018L0594. HUANG Motao, LIU Min, DENG Kailiang, et al. Analytical solution of downward continuation for airborne gravimetry based on upward continuation method[J]. Chinese Journal of Geophysics, 2018, 61(12), 4746-4757. DOI:10.6038/cjg2018L0594.
[14] TRAN K V, NGUYEN T N. A novel method for computing the vertical gradients of the potential field:application to downward continuation[J]. Geophysical Journal International, 2020, 220, 1316-1329. DOI:10.1093/gji/ggz524.
[15] 马在田,曹景忠,王家林, 等. 计算地球物理学[M]. 上海:同济大学出版社, 1997. MA Zaitian, CAO Jingzhong, WANG Jialin, et al. Computational geophysics[M]. Shanghai:Tongji University Press, 1997.
[16] SANSÒ F, SIDERIS M G. Geodetic boundary value problem:the equivalence between Molodensky's and Helmert's solutions[M]. New York:Springer, 2017.
[17] 刘敏,黄谟涛,邓凯亮, 等. 顾及地形效应的地面重力向上延拓模型分析与检验[J]. 武汉大学学报(信息科学版), 2018, 43(1):112-119. LIU Min, HUANG Motao, DENG Kailiang, et al. Test and analysis of upward continuation models for Earth surface gravity with regard to the effect of topographic height[J]. Geomatics and Information Science of Wuhan University, 2018, 43(1):112-119.
[18] 翟振和,王兴涛,李迎春. 解析延拓高阶解的推导方法与比较分析[J]. 武汉大学学报(信息科学版), 2015, 40(1):134-138. ZHAI Zhenhe, WANG Xingtao, LI Yingchun. Solution and comparison of high order term of analytical continuation[J]. Geomatics and Information Science of Wuhan University, 2015, 40(1):134-138.
[19] 刘长弘. 改进的直接积分方法计算低空扰动引力[D]. 郑州:信息工程大学, 2016. LIU Changhong. Improved direct integral methods calculating exterior disturbing gravity in low-altitude airspace[D]. Zhengzhou:Information Engineering University, 2016.
[20] 黄谟涛,刘敏,邓凯亮, 等. 海域流动点外部扰动引力无奇异计算模型[J]. 地球物理学报, 2019, 62(7):2394-2404. DOI:10.6038/cjg2019L0739. HUANG Motao, LIU Min, DENG Kailiang, et al. A nonsingular model for computing the external gravity field at a mobile point in a sea area[J]. Chinese Journal of Geophysics, 2019, 62(7):2394-2404. DOI:10.6038/cjg2019L0739.
[21] NOVÁK P, HECK B. Downward continuation and geoid determination based on band-limited airborne gravity data[J]. Journal of Geodesy, 2002, 76(5):269-278.
[22] 刘敏,黄谟涛,欧阳永忠, 等. 顾及地形效应的重力向下延拓模型分析与检验[J]. 测绘学报,2016, 45(5):521-530. LIU Min, HUANG Motao, OUYANG Yongzhong, et al. Test and analysis of downward continuation models for airborne gravity data with regard to the effect of topographic height[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(5):521-530.
[23] WONG L, GORE R. Accuracy of geodesy heights from modified Stokes kernels[J]. Geophysical Journal International, 1969, 18(1), 81-91.
[24] PAUL M. A method of evaluating the truncation error coefficients for geoidal heights[J]. Bulletin Geodesique, 1973, 110:413-425.
[25] NOVÁK P, KERN M, SCHWARZ K P. Numerical studies on the harmonic downward continuation of band-limited airborne gravity[J]. Studia Geophysica et Geodaetica, 2001, 45:327-345.
[26] ALBERTS B, KLEES R. A comparison of methods for the inversion of airborne gravity data[J]. Journal of Geodesy, 2004, 78(1-2):55-65.
[27] BJERHAMMAR A. A new theory of geodetic gravity[M]. Stockholm:Tekniska Hogskolan, 1964.
[28] MORITZ H. The figure of the Earth:theoretical geodesy and the Earth's interior[M]. Karlsruhe:Wichmann, Germany, 1990.
[29] TROMPAT H, BOSCHETTI F, HORNBY P. Improved downward continuation of potential field data[J]. Exploration Geophysics, 2003, 34(4):249-256.
[30] MARTINEC Z. Stability investigations of a discrete downward continuation problem for geoid determination in the Canadian rocky mountains[J]. Journal of Geodesy, 1996(70):805-828.
[31] PAVLIS N K, HOLMES S A, KENYON S C, et al. The development and evaluation of the earth gravitational model 2008(EGM2008)[J]. Journal Geophysical Research, 2012, 117(B4):B04406. DOI:10.1029/2011JB008916.

重力异常向上延拓严密改化模型及向下延拓应用

Rigorous modification model of upward continuation and its applications on the downward continuation of gravity anomaly

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 8

编辑推荐

Metrics

本文评价

[1]	袁运斌, 李敏, 霍星亮, 李子申, 王宁波. 北斗三号全球导航卫星系统全球广播电离层延迟修正模型(BDGIM)应用性能评估[J]. 测绘学报, 2021, 50(4): 436-447.
[2]	巩秀强, 袁俊军, 胡小工, 王彬, 陈俊平, 周善石. 北斗广播电文钟差模型精度评估及改善策略[J]. 测绘学报, 2021, 50(2): 181-188.
[3]	章传银, 马旭, 章磊, 丁剑. 基于GNSS水准和重力场误差特性的大地水准面精度评估方法[J]. 测绘学报, 2021, 50(1): 12-17.
[4]	刘伟平, 郝金明, 吕志伟, 谢建涛, 刘婧, 焦博. 北斗三号空间信号测距误差评估与对比分析[J]. 测绘学报, 2020, 49(9): 1213-1221.
[5]	邓晨, 游雄, 张威巍, 智梅霞. 基于2D地图的城市户外ARGIS视觉辅助地理配准技术[J]. 测绘学报, 2019, 48(10): 1305-1319.
[6]	陈明, 郭际明, 武军郦, 宋伟伟, 张龙平. 北斗精密定轨精度评估及站点分布影响分析[J]. 测绘学报, 2018, 47(S0): 93-100.
[7]	翟振和, 孙中苗, 李迎春, 肖云. 航空重力测量在近海区域的精度评估与分析[J]. 测绘学报, 2015, 44(1): 1-5.
[8]	李姗姗,吴晓平,张传定,欧阳永忠. 顾及地形与完全球面布格异常梯度项改正的区域似大地水准面精化[J]. 测绘学报, 2012, 41(4): 0-590.