| 基于Lp范数正则化的V字型密度界面重力反演 |
| |
| 引用本文： | 冯旭亮,王万银,宋立军,袁炳强.基于Lp范数正则化的V字型密度界面重力反演[J].地球物理学报,2019,62(3):1022-1036. |
| |
| 作者姓名： | 冯旭亮 王万银 宋立军 袁炳强 |
| |
| 作者单位： | 1. 西安石油大学地球科学与工程学院, 西安 710065;2. 长安大学重磁方法技术研究所, 长安大学地质工程与测绘学院, 西安 710054 |
| |
| 基金项目： | 陕西省自然科学基础研究计划项目（2018JQ4034）、陕西省教育厅科研计划项目（17JK0617）以及中国地质调查局项目（DD20160227，DD20160227-05）联合资助. |
| |
| 摘 要： | V字型密度界面是一类常见的密度界面，如海沟、半地堑以及俯冲带之下的莫霍面，利用重力数据刻画此类密度界面形态对于区域构造研究、油气勘探以及物理海洋学等都具有重要意义.本文首先建立了Lp-范数形式的模型约束函数，并利用正则化原理将其与重力数据误差函数和已知深度约束函数结合形成V字型密度界面反演的目标函数，推导了目标函数的梯度表达式，并以非线性共轭梯度法为核心给出了反演流程.二维简单模型试算结果表明p=5时该方法能准确地刻画V字型密度界面起伏特征，且亦能准确地应用于二维复杂密度界面和三维界面的反演.最后将反演方法应用于挑战者深渊及邻区的实际资料处理之中，利用研究区海底地形数据和沉积层厚度数据对自由空间重力异常逐层剥离而得到莫霍面引起的重力异常，用本文方法对此重力异常进行反演，结果呈现了板块俯冲作用引起的V字型莫霍面起伏特征.
|
| 关 键 词： | V字型密度界面 重力异常 Lp-范数 非线性反演 正则化 |
| 收稿时间： | 2017-12-28 |
Gravity inversion for V-shaped density interface based on Lp-norm regularization |
| |
| FENG XuLiang,WANG WanYin,SONG LiJun,YUAN BingQiang.Gravity inversion for V-shaped density interface based on Lp-norm regularization[J].Chinese Journal of Geophysics,2019,62(3):1022-1036. |
| |
| Authors: | FENG XuLiang WANG WanYin SONG LiJun YUAN BingQiang |
| |
| Institution: | 1. School of Earth Sciences and Engineering, Xi'an Shiyou University, Xi'an 710065, China;2. Institute of Gravity and Magnetic Technology, College of Geology Engineering and Geomatics, Chang'an University, Xi'an 710054, China |
| |
| Abstract: | The V-shaped density interface is commonly seen such as the ocean trench, half-graben, and the Moho beneath a subduction zone. Delineating the relief of these V-shaped density interfaces is of significance in regional tectonic research, physical oceanography and oil and gas exploration.
Among inversion methods of the density interface, the regularization approach, which applies prior information, can guarantee the inversion result to fit the known depth of the interface to be estimated. It can also control the shape of the interface via the model constraint function in the regularization item, making the inversion result agree with the real geologic feature. We have designed a model constraint function in the form of Lp-norm and then integrated it with the gravity data misfit function and the known depth constraint function to establish an objective function for V-shaped density interface inversion. Subsequently, we have derived the gradient of the objective function and built the inversion process taking the nonlinear conjugate gradient algorithm as the core.#br#The effect of our proposed inversion method is related to the value of p. The test results conducted with a two-dimensional simple model show that the method can be used to accurately estimate the V-shaped density interface relief when p=5, and the correctness of the proposed method is further confirmed by a two-dimensional V-shaped density interface with complex shape. The test results of a three-dimensional V-shaped density interface show that our proposed method is also appropriate for three-dimensional inversion. At the end of this article, the inversion method is tested with real gravity data of the Challenger Deep and adjacent regions. The gravity anomalies caused by the Moho were calculated by removing extra gravity effects layer by layer from the free-air gravity anomalies under the constraint of the submarine topography and sedimentary thickness data in this area. Inversion of these gravity anomalies using our proposed method shows that the method is able to delineate the relief of a V-shaped Moho under a subduction zone. The estimated depth of the Moho beneath the Challenger Deep ranges 18 km to 20 km, while the depth at the edge of the western Pacific Ocean is 8 km to 12 km, which, to some extent, suggests crustal complexity in the study area. Furthermore, a large inclination of the Pacific subduction slab is inferred from the locations of the maximum of the submarine topography and Moho depth in this area, indicating a relatively slow convergence between the Pacific and Philippine plates. |
| |
| Keywords: | V-shaped density interface Gravity anomaly Lp-norm Nonlinear inversion Regularization |
| 本文献已被 CNKI 等数据库收录！ |
| 点击此处可从《地球物理学报》浏览原始摘要信息 |
| 点击此处可从《地球物理学报》下载免费的PDF全文 |
|
