测绘学报 ›› 2024, Vol. 53 ›› Issue (8): 1552-1563.doi: 10.11947/j.AGCS.2024.20230161

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

基于非负约束的重力向下延拓BTTB-MRNSD方法

刘天佑1(), 曾小牛1,2(), 李夕海1   

  1. 1.火箭军工程大学核工程学院,陕西 西安 710025
    2.海洋油气勘探国家工程研究中心,北京 100028
  • 收稿日期:2023-05-24 发布日期:2024-09-25
  • 通讯作者: 曾小牛 E-mail:tianyou_liu1@163.com;tianyou_liu1@163.com;xiaoniuzeng@163.com
  • 作者简介:刘天佑(2000—),男,博士生,主要研究方向为重磁数据分析与处理。E-mail:tianyou_liu1@163.com
  • 基金资助:
    国家自然科学基金(41804136)

BTTB-MRNSD method for downward continuation of gravity field based on nonnegative constraints

Tianyou LIU1(), Xiaoniu ZENG1,2(), Xihai LI1   

  1. 1.School of Nuclear Engineering, Rocket Force University of Engineering, Xi'an 710025, China
    2.National Engineering Research Center of Offshore Oil and Gas Exploration, Beijing 100028, China
  • Received:2023-05-24 Published:2024-09-25
  • Contact: Xiaoniu ZENG E-mail:tianyou_liu1@163.com;tianyou_liu1@163.com;xiaoniuzeng@163.com
  • About author:LIU Tianyou (2000—), male, PhD candidate, majors in gravity and magnetic data analysis and processing. E-mail: tianyou_liu1@163.com
  • Supported by:
    The National Natural Science Foundation of China(41804136)

摘要:

向下延拓作为处理与解释重力数据的一项重要技术,因其计算的不适定性而备受研究人员的关注。空间域延拓方法一般延拓精度较高,但通常计算复杂度较大。本文首先根据重力异常数据的特点,借鉴图像复原问题的处理思路,提出了一种重力向下延拓等效数学模型。然后,基于系数矩阵的块-托普利茨-托普利茨-块(BTTB)结构特点,本文提出了一种具有非负约束的空间-波数混合域下延迭代方法。该方法克服了空间域延拓计算复杂度大的问题,计算效率较高。基于理论重力模型和实测异常数据的对比试验表明,本文提出的重力下延方法具有相对较高的下延精度和稳定性,且收敛性较好。

关键词: 重力, 非负约束, 向下延拓, 最速下降法, BTTB矩阵

Abstract:

As an important technology for processing and interpreting gravity data, downward continuation has attracted the attention of researchers because of its ill-posedness. The space domain continuation method generally has high continuation accuracy, but the computational complexity is usually large. In this paper, firstly, according to the characteristics of gravity anomaly data and the idea of image restoration, an equivalent mathematical model of gravity downward continuation is proposed. Then, based on the block-Toeplitz Toeplitz-block (BTTB) structure of the coefficient matrix, we propose a space-wavenumber mixed domain downward continuation iterative method with nonnegative constraints. The method overcomes the disadvantage of large computational complexity in space domain continuation and has high computational efficiency. The comparison experiments of theoretical gravity model and real anomaly data show that the gravity field downward continuation method proposed in this paper has high downward continuation accuracy and stability, and has good convergence.

Key words: gravity, nonnegative constraints, downward continuation, steepest descent method, BTTB matrix

中图分类号: