利用高阶径向导数带限模型进行重力向下延拓计算

邓凯亮; 黄谟涛; 吴太旗; 王伟平; 欧阳永忠; 陈欣; 熊雄; 刘敏; 王许

doi:10.13203/j.whugis20210630

利用高阶径向导数带限模型进行重力向下延拓计算

doi: 10.13203/j.whugis20210630

1 海军研究院, 天津 300061;
2 自然资源部海洋环境探测技术与应用重点实验室, 广州 510300;
3 国家海洋局南海调查技术中心, 广州 5103004 91001 部队, 北京 100830

基金项目:

国家自然科学基金项目（42174013，41804011，41774021）；国家重点研发计划（2016YFB0501704，2016YFC0303007）；军科委173计划（2019-JCJQ-ZD-017；2020-JCJQ-ZD-139）

详细信息

作者简介:
邓凯亮,博士,高级工程师,主要从事海洋重力场理论方法研究。dengkailiang036@163.com

Downward Continuation of Gravity Using the Band-Limited Models for High-Order Radial Derivatives of Gravity Anomaly

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;

Funds:

The National Natural Science Foundation of China, Nos.42174013, 41804011, 41774021

摘要: 泰勒（Taylor）级数展开是实施位场向下延拓解算的主要方法之一。该方法的有效性主要取决于位场延拓参量各阶垂向（或径向）偏导数的求取精度及其可靠性。为了避免使用封闭解析核函数在球边界面出现奇异性带来的不确定性问题，本文依据各类重力观测经滤波处理后均表现为一类有限频谱带宽信号的特点，提出将地球外部重力异常Poisson积分式的核函数表示为球谐级数展开式，并将其截断为与重力观测值频谱范围相一致的带限求和式，进而通过直接求导方法，推导得到一组与带限核函数相对应的重力异常高阶径向导数带限计算公式，同时对该组公式进行了实用性改化，最后将其应用于重力异常向下延拓泰勒级数展开式计算。使用超高阶地球位模型EGM2008并设计两个阶段的数值计算检验方案，分别对重力异常高阶径向偏导数带限模型及其Taylor级数展开向下延拓模型的计算精度进行了检核评估，表明新模型具有良好的可靠性和有效性，在解算稳定性和计算精度两个方面都优于其他同类模型。
- 重力向下延拓 /
- 重力异常径向导数 /
- 带限模型 /
- 封闭解析核函数 /
- 球谐级数展开截断式 /
- 泰勒级数展开
Abstract: Taylor series expansion is often used in the downward continuation of potential field, and its performance depends on the accuracy and reliability of the vertical or radial partial derivatives (VPDs or RPDs) of potential field parameters. In order to avoid the singularity on spherical boundary and the uncertainty to the computational results by using the closed analytic kernel function to solve the partial derivative, considering the fact that all kinds of gravity observations behave as a type of limited spectrum bandwidth signal after being filtered, this research proposes to express the kernel function of the Poisson integral for the external gravity anomaly by a spherical harmonic series expansion, which is then truncated into a band-limited summation that has the same spectrum range as the gravity observation. After that, we derive a set of band-limited formulas to calculate the high-order RPDs, which are modified and applied to the downward continuation of the gravity anomaly by Taylor series expansion. The formulas are validated using the ultra-high-degree geopotential model EGM2008 by a two-stage procedure. The numerical tests of the band-limited formulas and the Taylor series expansion downward continuation model show that the proposed band-limited formulas have good reliability and validity, and are superior to other models in terms of stability and accuracy.

[1]	Moritz H. Advanced Physical Geodesy[M]. Karlsruhe:Herbert Wichmann Verlag, 1980
[2]	Sebera J, Šprlák M, Novák P, et al. Iterative Spherical Downward Continuation Applied to Magnetic and Gravitational Data from Satellite[J]. Surveys in Geophysics, 2014, 35(4):941-958
[3]	Pitoňák M, Novák P, Eshagh M, et al. Downward Continuation of Gravitational Field Quantities to an Irregular Surface by Spectral Weighting[J]. Journal of Geodesy, 2020, 94(7):1-26
[4]	Leão J W D, Silva J B C. Discrete Linear Transformations of Potential Field Data[J]. GEOPHYSICS, 1989, 54(4):497-507
[5]	Pawlowski R S. Preferential Continuation for Potential-Field Anomaly Enhancement[J]. GEOPHYSICS, 1995, 60(2):390-398
[6]	Fedi M, Florio G. A Stable Downward Continuation by Using the ISVD Method[J]. Geophysical Journal International, 2002, 151(1):146-156
[7]	Trompat H, Boschetti F, Hornby P. Improved Downward Continuation of Potential Field Data[J]. Exploration Geophysics, 2003, 34(4):249-256
[8]	Heiskanen W A, Moritz H. Physical Geodesy[M]. San Francisco:Freeman and Company, 1967
[9]	Martinec Z. Stability Investigations of a Discrete Downward Continuation Problem for Geoid Determination in the Canadian Rocky Mountains[J]. Journal of Geodesy, 1996, 70(11):805-828
[10]	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
[11]	SansòF, Sideris M. Geoid Determination:Theory and Methods[M]. Berlin:Springer, 2013
[12]	Xu P L. Truncated SVD Methods for Discrete Linear Ill-Posed Problems[J]. Geophysical Journal International, 1998, 135(2):505-514
[13]	Hansen P C, O'Leary D P. The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems[J]. SIAM Journal on Scientific Computing, 1993, 14(6):1487-1503
[14]	Kern M. An Analysis of the Combination and Downward Continuation of Satellite, Airborne and Terrestrial Gravity Data[D]. Calgary:University of Calgary, 2003
[15]	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
[16]	Hwang C, Hsiao Y S, Shih H C, et al. Geodetic and Geophysical Results from a Taiwan Airborne Gravity Survey:Data Reduction and Accuracy Assessment[J]. Journal of Geophysical Research:Solid Earth, 2007, 112(B4):93-101
[17]	Bjerhammar A. A New Theory of Geodetic Gravity[M]. Stockholm:Tekniska Hogskolan, 1964
[18]	Dampney C N G. The Equivalent Source Technique[J]. GEOPHYSICS, 1969, 34(1):39-53
[19]	Sünkel H. The Generation of a Mass Point Model from Surface Gravity Data[R]. Ohio:Ohio State University, 1983
[20]	Boschetti F, Therond V, Hornby P. Feature Removal and Isolation in Potential Field Data[J]. Geophysical Journal International, 2004, 159(3):833-841
[21]	Novák P, Kern M, Schwarz K. Numerical Studies on the Harmonic Downward Continuation of Band-Limited Airborne Gravity[J]. Studia Geophysica et Geodaetica, 2001, 45:327-345
[22]	Novák P, Kern M, Schwarz K P, et al. On Geoid Determination from Airborne Gravity[J]. Journal of Geodesy, 2003, 76(9/10):510-522
[23]	Kern M, Schwarz K K P P, Sneeuw N. A Study on the Combination of Satellite, Airborne, and Terrestrial Gravity Data[J]. Journal of Geodesy, 2003, 77(3/4):217-225
[24]	Mansi A H, Capponi M, Sampietro D. Downward Continuation of Airborne Gravity Data by Means of the Change of Boundary Approach[J]. Pure and Applied Geophysics, 2018, 175(3):977-988
[25]	Ma G Q, Liu C, Huang D N, et al. A Stable Iterative Downward Continuation of Potential Field Data[J]. Journal of Applied Geophysics, 2013, 98:205-211
[26]	Zhang C, Lü Q, Yan J Y, et al. Numerical Solutions of the Mean-Value Theorem:New Methods for Downward Continuation of Potential Fields[J]. Geophysical Research Letters, 2018, 45(8):3461-3470
[27]	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, 2019, 220(2):1316-1329
[28]	Wei Ziqing. High-Order Radial Derivatives of Harmonic Function and Gravity Anomaly[J]. Journal of Physical Science and Application, 2014, 4(7):454-467
[29]	Wong L, Gore R. Accuracy of Geoid Heights from Modified Stokes Kernels[J]. Geophysical Journal International, 1969, 18(1):81-91
[30]	Vaníček P, Featherstone W E. Performance of Three Types of Stokes's Kernel in the Combined Solution for the Geoid[J]. Journal of Geodesy, 1998, 72(12):684-697
[31]	Wang Y M, Saleh J, Li X, et al. The US Gravimetric Geoid of 2009(USGG2009):Model Development and Evaluation[J]. Journal of Geodesy, 2012, 86(3):165-180
[32]	Pavlis N K, Holmes S A, Kenyon S C, et al. The Development and Evaluation of the Earth Gravitational Model 2008(EGM2008)[J]. Journal of Geophysical Research:Solid Earth, 2012, 117(B4)

[1]	杜劲松, 梁青, 张祎, 赵艳彬, 廖鹤. 基于径向泰勒级数展开的月球表面磁场快速球谐解算方法 . 武汉大学学报 ● 信息科学版, doi: 10.13203/j.whugis20190164
[2]	徐新强, 赵俊. 航空重力向下延拓的多参数正则化方法 . 武汉大学学报 ● 信息科学版, doi: 10.13203/j.whugis20180335
[3]	伍丰丰, 黄海军, 任青阳, 樊文有, 陈洁, 潘雄. 综合半参数核估计与正则化方法的航空重力向下延拓模型分析 . 武汉大学学报 ● 信息科学版, doi: 10.13203/j.whugis20180491
[4]	刘敏, 黄谟涛, 邓凯亮, 欧阳永忠, 翟国君, 吴太旗. 顾及地形效应的地面重力向上延拓模型分析与检验 . 武汉大学学报 ● 信息科学版, doi: 10.13203/j.whugis20150519
[5]	陈欣, 翟国君, 暴景阳, 欧阳永忠, 陆秀平, 邓凯亮. 航空重力向下延拓的最小二乘配置Tikhonov正则化法 . 武汉大学学报 ● 信息科学版, doi: 10.13203/j.whugis20150728
[6]	李建成, 徐新禹, 赵永奇, 万晓云. 由GOCE引力梯度张量不变量确定卫星重力模型的半解析法 . 武汉大学学报 ● 信息科学版, doi: 10.13203/j.whugis20150554
[7]	肖雄武, 郭丙轩, 潘飞, 张春森, 薛万唱. 利用泰勒展开的点特征子像素定位方法 . 武汉大学学报 ● 信息科学版,
[8]	黄谟涛, 欧阳永忠, 刘敏, 翟国君, 邓凯亮. 海域航空重力测量数据向下延拓的实用方法 . 武汉大学学报 ● 信息科学版,
[9]	刘鹏程, 艾廷华, 毕旭. 傅立叶级数支持下的等高线多尺度表达模型 . 武汉大学学报 ● 信息科学版,
[10]	顾勇为, 归庆明, 韩松辉, 王靳辉. 航空重力向下延拓分组修正的正则化解法 . 武汉大学学报 ● 信息科学版,
[11]	翟振和, 孙中苗, 王兴涛, 刘晓刚. 利用级数展开法确定区域大地水准面的重力位 . 武汉大学学报 ● 信息科学版,
[12]	杜劲松, 陈超, 梁青, 张毅. 基于球冠滑动平均的球面重力异常分离方法 . 武汉大学学报 ● 信息科学版,
[13]	朱广彬, 李建成, 文汉江, 徐新禹. 极空白对三类重力梯度边值问题球谐分析解的影响分析 . 武汉大学学报 ● 信息科学版,
[14]	王伟, 宁津生, 罗志才, 汪海洪. 利用小波变换分离俯冲带重力异常 . 武汉大学学报 ● 信息科学版,
[15]	李宗春, 邓勇, 张冠宇, 杨晓晖. 变形测量异常数据处理中小波变换最佳级数的确定 . 武汉大学学报 ● 信息科学版,
[16]	黄贤源, 隋立芬, 范澎湃. 有色观测噪声的随机模型级数表示及其补偿法 . 武汉大学学报 ● 信息科学版,
[17]	申文斌, 鄢建国, 晁定波. 重力场的局部虚拟向下延拓以及利用EGM96模型的模拟实验检验 . 武汉大学学报 ● 信息科学版,
[18]	罗志才, 晁定波, 宁津生. 重力梯度张量的球谐分析 . 武汉大学学报 ● 信息科学版,
[19]	罗志才, 宁津生, 晁定波. 卫星重力梯度分量的球谐综合 . 武汉大学学报 ● 信息科学版,
[20]	白亿同, 罗玉芳, 胡永旭. 分数微分积分及其级数展开式 . 武汉大学学报 ● 信息科学版,

点击查看大图

计量

文章访问数: 49
HTML全文浏览量: 9
PDF下载量: 3
被引次数: 0

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

利用高阶径向导数带限模型进行重力向下延拓计算

doi: 10.13203/j.whugis20210630

作者简介:
邓凯亮,博士,高级工程师,主要从事海洋重力场理论方法研究。dengkailiang036@163.com

Downward Continuation of Gravity Using the Band-Limited Models for High-Order Radial Derivatives of Gravity Anomaly

计量

利用高阶径向导数带限模型进行重力向下延拓计算

doi: 10.13203/j.whugis20210630

1 海军研究院, 天津 300061;

2 自然资源部海洋环境探测技术与应用重点实验室, 广州 510300;

3 国家海洋局南海调查技术中心, 广州 5103004 91001 部队, 北京 100830

作者简介:
邓凯亮,博士,高级工程师,主要从事海洋重力场理论方法研究。dengkailiang036@163.com

English Abstract

Downward Continuation of Gravity Using the Band-Limited Models for High-Order Radial Derivatives of Gravity Anomaly

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;

全文HTML

目录

留言板

利用高阶径向导数带限模型进行重力向下延拓计算

doi: 10.13203/j.whugis20210630

作者简介: 邓凯亮,博士,高级工程师,主要从事海洋重力场理论方法研究。dengkailiang036@163.com

Downward Continuation of Gravity Using the Band-Limited Models for High-Order Radial Derivatives of Gravity Anomaly

计量

出版历程

利用高阶径向导数带限模型进行重力向下延拓计算

doi: 10.13203/j.whugis20210630

1 海军研究院, 天津 300061; 2 自然资源部海洋环境探测技术与应用重点实验室, 广州 510300; 3 国家海洋局南海调查技术中心, 广州 5103004 91001 部队, 北京 100830

作者简介: 邓凯亮,博士,高级工程师,主要从事海洋重力场理论方法研究。dengkailiang036@163.com

English Abstract

Downward Continuation of Gravity Using the Band-Limited Models for High-Order Radial Derivatives of Gravity Anomaly

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;

全文HTML

目录

作者简介:
邓凯亮,博士,高级工程师,主要从事海洋重力场理论方法研究。dengkailiang036@163.com

1 海军研究院, 天津 300061;

2 自然资源部海洋环境探测技术与应用重点实验室, 广州 510300;

3 国家海洋局南海调查技术中心, 广州 5103004 91001 部队, 北京 100830

作者简介:
邓凯亮,博士,高级工程师,主要从事海洋重力场理论方法研究。dengkailiang036@163.com

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;