Abstract：Three-dimensional inversion of gravity gradient data can obtain a subsurface density structure model,which can be used for geological resource exploration and other fields.Gravity gradient data from airborne,ground,and borehole observations contain information at different frequencies.Joint inversion of these data can reduce the non-uniqueness of inversion and improve imaging resolution.However,for underground anomalous bodies with complex morphology,the vertical resolution of current joint inversion of multi-scale data,especially the imaging resolution of the anomalous body bottom still needs improvement.To address this issue,the joint inversion of airborne-ground-borehole vertical gravity gradient data was studied in this paper.Firstly,a fuzzy C-means clustering algorithm was introduced into the regularization inversion,and the clustering constraint was added during the iteration process to reduce non-uniqueness.Secondly,a joint inversion method was proposed by combining airborne,ground,and borehole vertical gravity gradient data,and GPU acceleration was utilized for computation.Thirdly,the inversion was applied to theoretical model data and gravity gradient data from the Vinton Dome,USA,so as to verify the effectiveness of the method and discuss the influence of borehole locations on the results.Finally,the performance of GPU acceleration-based parallel inversion was analyzed.The data test demonstrates that the fuzzy C-means clustering can reduce the non-uniqueness of inversion,and the joint inversion can obtain accurate density distributions with certain anti-noise ability.Using the data from anomaly-side and anomaly-through boreholes can improve imaging resolution.The computed density distribution of the Vinton Dome area is close to the conclusions of other researchers,demonstrating the effectiveness and feasibility of the proposed method.The test results also show that the GPU parallel method has a high acceleration ratio,and the proposed method can provide technical support for geological explorations and other studies.
ZHANG Henglei,GENG Meixia,HU Xiangyun.Three-dimensional inversion of gravity/magnetic anomalies based on curvelet compression and its applications[J].Oil Geophysical Prospecting,2023,58(4):993-1001.张恒磊,耿美霞,胡祥云.基于曲波压缩的重磁异常三维反演及其应用[J].石油地球物理勘探,2023,58(4):993-1001.
[2]
XIANG Peng,TAN Shaoquan,CHEN Xueguo,et al.Gravity inversion method based on quasi-neural network featuring Gaussian radial basis function[J].Oil Geophysical Prospecting,2021,56(6):1409-1418.相鹏,谭绍泉,陈学国,等.利用高斯径向基函数的拟神经网络重力反演方法[J].石油地球物理勘探,2021,56(6):1409-1418.
[3]
QIN Pengbo,HUANG Danian.Integrated gravity and gravity gradient data focusing inversion[J].Chinese Journal of Geophysics,2016,59(6):2203-2224.秦朋波,黄大年.重力和重力梯度数据联合聚焦反演方法[J].地球物理学报,2016,59(6):2203-2224.
[4]
PILKINGTON M.Analysis of gravity gradiometer inverse problems using optimal design measures[J].Geophysics,2012,77(2):G25-G31.
[5]
PILKINGTON M.Evaluating the utility of gravity gradient tensor components[J].Geophysics,2014,79(1):G1-G14.
[6]
LI Y,OLDENBURG D W.Joint inversion of surface and three-component borehole magnetic data[J].Geophysics,2000,65(2):540-552.
[7]
ZHANG Dalian,LIU Tianyou,CHEN Shixian,et al.The application of PSO to joint inversion of surface and borehole magnetic data[J].Geophysical and Geochemical Exploration,2009,33(5):571-575,591.张大莲,刘天佑,陈石羡,等.粒子群算法在磁测资料井地联合反演中的应用[J].物探与化探,2009,33(5):571-575,591.
[8]
GENG M,YANG Q,HUANG D.3D joint inversion of gravity-gradient and borehole gravity data[J].Exploration Geophysics,2017,48(2):151-165.
[9]
WANG Taihan,MA Guoqing,XIONG Shengqing,et al.Joint regularized density inversion method of airborne,surface and borehole gravity anomaly data[J].Chinese Journal of Geophysics,2020,63(7):2737-2750.王泰涵,马国庆,熊盛青,等.空-地-井重力异常正则化协同密度反演方法[J].地球物理学报,2020,63(7):2737-2750.
[10]
GAO Xiuhe,XIONG Shengqing,YU Changchun,et al.A study of joint inversion of gravity data from multi-planes and boreholes[J].Geophysical and Geochemical Exploration,2020,44(6):1361-1367.高秀鹤,熊盛青,于长春,等.不同高度重力数据和井中重力数据融合反演研究[J].物探与化探,2020,44(6):1361-1367.
[11]
PENG G M,SUN Z Y,LIU Z.3D joint inversion of surface and borehole gravity data using zeroth-order minimum entropy regularization[J].Applied Geophysics,2021,18(2):131-144.
[12]
YANG Peijie,YIN Xingyao,ZHANG Guangzhi.Cluster analysis of seismic attributes by fuzzy C-mean algorithm[J].Oil Geophysical Prospecting,2007,42(3):322-324.杨培杰,印兴耀,张广智.模糊C均值地震属性聚类分析[J].石油地球物理勘探,2007,42(3):322-324.
[13]
SUN J,LI Y.Geophysical inversion using petrophysical constraints with application to lithology differentiation[C].SEG Technical Program Expanded Abstracts,2011,30:2644-2648.
[14]
MAAG E,LI Y.Discrete-valued gravity inversion using the guided fuzzy c-means clustering technique[J].Geophysics,2018,83(4):G59-G77.
[15]
XU Kaijun,LI Meng,JI Chunhui,et al.Data space joint inversion of ground and borehole magnetic data with fuzzy C-means clustering constraints[J].Journal of China University of Petroleum (Edition of Natural Science),2021,45(3):55-64.徐凯军,李猛,季春晖,等.基于模糊C均值聚类约束的井-地磁法联合数据空间反演[J].中国石油大学学报(自然科学版),2021,45(3):55-64.
[16]
YI Ke,ZHANG Zhiyong,LI Man,et al.Joint inversion of resistivity and permittivity for two dimensional RMT data based on FCM clustering[J].Chinese Journal of Geophysics,2022,65(6):2340-2350.易柯,张志勇,李曼,等.基于FCM聚类算法的二维RMT电阻率与介电常数联合反演[J].地球物理学报,2022,65(6):2340-2350.
[17]
SUN J,LI Y.Multidomain petrophysically constrained inversion and geology differentiation using guided fuzzy C-means clustering[J].Geophysics,2015,80(4):ID1-ID18.
[18]
LI Y,OLDENBURG D W.3-D inversion of magnetic data[J].Geophysics,1996,61(2):394-408.
[19]
TIKHONOV A N,ARSENINV Y.Solutions of Ill-posed Problems[M].Halsted Press,New York,1977.
[20]
PORTNIAGUINE O,ZHDANOV M S.3-D magnetic inversion with data compression and image focusing[J].Geophysics,2000,67(5):1532-1541.
[21]
CHEN Zhaoxi,MENG Xiaohong,GUO Lianghui,et al.Three-dimensional fast forward modeling and the inversion strategy for large scale gravity and gravimetry data based on GPU[J].Chinese Journal of Geophysics,2012,55(12):4069-4077.陈召曦,孟小红,郭良辉,等.基于GPU并行的重力、重力梯度三维正演快速计算及反演策略[J].地球物理学报,2012,55(12):4069-4077.
[22]
HOU Z,SUN B,QIN P B,et al.Joint nonlinear inversion of full tensor gravity gradiometry data and its parallel algorithm[J].IEEE Transactions on Geoscience and Remote Sensing,2022,60:1-12.
[23]
COKER M O,BHATTACHARYA J P,MARFURT K J.Fracture patterns within mudstones on the flanks of a salt dome:syneresis or slumping?[J].Gulf Coast Association of Geological Societies,2007,57:125-137.
[24]
OLIVEIRA V C,BARBOSA V C F.3-D radial gravity gradient inversion[J].Geophysical Journal International,2013,195(2):883-902.
