Abstract：Gravity exploration blocks are often irregular,which causes the vacancy in the acquired gravity data.Before the processing and transformation of gravity data in wavenumber domain,the original data must be interpolated and processed with edge padding.High-frequency noise in gravity data is the main factor which causes the instability in later processing.Conventional gravity data processing generally performs interpolation,denoising and edge padding independently. These three issues were considered integrally,and an iterative method for simultaneous interpolation,denoising and edge padding of gravity data based on the projection onto convex sets was proposed.Firstly,the final cutoff wavenumber of the iterative method was determined through calculating and fitting the radial average power spectrum of the gravity data.Secondly,interpolation,edge padding and denoising were applied to gravidy data using the spectrum with the wavenumber lower than the cutoff wavenumber,until the preset iteration number was reached.Theoretical gravity model test and the application in real isostatic gravity data acquired in Afghanistan showed that the method proposed is simple theoretically and convenient to be applied.The splicing of the interpolation and edge padding results is smooth without distortion,and the method achieved good interpolation and denoising effect.The results of the method proposed in this paper are better,compared with those of the conventional joint processing method based on minimum curvature,Kriging interpolation,wavelet denosing and cosine edge padding.
曾小牛, 李夕海, 侯维君, 刘继昊. 基于凸集投影的重力同时填充扩边和去噪方法[J]. 石油地球物理勘探, 2020, 55(1): 197-205.
ZENG Xiaoniu, LI Xihai, HOU Weijun, LIU Jihao. Simultaneous interpolation,edge padding and denoising method for gravity data based on the projection onto convex sets. Oil Geophysical Prospecting, 2020, 55(1): 197-205.
王万银,邱之云,刘金兰,等.位场数据处理中的最小曲率扩边和补空方法研究[J].地球物理学进展,2009,24(4):1327-1338.WANG Wangyin,QIU Zhiyun,LIU Jinlan,et al.The research to the extending edge and interpolation based on the minimum curvature method in potential field data processing[J].Progress in Geophysics,2009,24(4):1327-1338.
[2]
骆遥,吴美平.位场向下延拓的最小曲率方法[J].地球物理学报,2016,59(1):240-251.LUO Yao,WU Meiping.Minimum curvature method for downward continuation of potential field data[J].Chinese Journal of Geophysics,2016,59(1):240-251.
[3]
Pawlowski R S.Preferential continuation for potential-field anomaly enhancement[J].Geophysics,1995,60(2):390-398.
[4]
Martinez C,Li Y.Denoising of gravity gradient data using an equivalent source technique[J].Geophysics,2016,81(4):G67-G79.
[5]
张双喜,陈超,王林松,等.二维经验模态分解及其在位场去噪和分离中的应用[J].地球物理学进展,2015,30(6):2855-2862.ZHANG Shuangxi,CHEN Chao,WANG Linsong,et al.The bidimensional empirical mode decomposition and its applications to denoising and separation of potential field[J].Progress in Geophysics,2015,30(6):2855-2862.
[6]
张旭东,詹毅,马永琴.不同信号的小波变换去噪方法[J].石油地球物理勘探,2007,42(增刊1):118-123.ZHANG Xudong,ZHAN Yi,MA Yongqin.Approaches of denoise by wavelet transform of different signals[J].Oil Geophysical Prospecting,2007,42(S1):118-123.
[7]
闫浩飞,刘国峰,薛典军,等.基于凸集投影方法的重磁数据规则缺失重建[J].地球物理学进展,2016,31(5):2192-2197.YAN Haofei,LIU Guofeng,XUE Dianjun,et al.Reconstruction of gravity/magnetic data with the projection-onto-convex-sets methods[J].Progress in Geophysics,2016,31(5):2192-2197.
[8]
Li X,Götze H J.Comparison of some gridding methods[J].The Leading Edge,1999,18(8):898-900.
[9]
Wang S Q,Zhang J H.Fast image inpainting using exponential-threshold POCS plus conjugate gradient[J].The Imaging Science Journal,2014,62(3):161-170.
[10]
Abma R,Kabir N.3D interpolation of irregular data with a POCS algorithm[J].Geophysics,2006,71(6):E91-E97.
[11]
Gao J J,Stanton A,Naghizadeh M,et al.Convergence improvement and noise attenuation considerations for beyond alias projection onto convex sets reconstruction[J].Geophysical Prospecting,2013,61(S1):138-151.
[12]
曹静杰,王本锋.基于一种改进凸集投影方法的地震数据同时插值和去噪[J].地球物理学报,2015,58(8):2935-2947.CAO Jingjie,WANG Benfeng.An improved projection onto convex sets method for simultaneous interpolation and denoising[J].Chinese Journal of Geophysics,2015,58(8):2935-2947.
[13]
刘国昌,陈小宏,郭志峰,等.基于Curvelet变换的缺失地震数据插值方法[J].石油地球物理勘探,2011,46(2):237-246.LIU Guochang,CHEN Xiaohong,GUO Zhifeng,et al.Missing seismic data rebuilding by interpolation based on Curvelet transform[J].Oil Geophysical Prospecting,2011,46(2):237-246.
[14]
王本锋,陈小宏,李景叶,等.POCS联合改进的Jitter采样理论曲波域地震数据重建[J].石油地球物理勘探,2015,50(1):20-28.WANG Benfeng,CHEN Xiaohong,LI Jingye,et al.Seismic data reconstruction based on POCS and improved Jittered sampling in the curvelet domain[J].Oil Geophysical Prospecting,2015,50(1):20-28.
[15]
陈小春,陈辉,喻勤,等.反假频POCS数据规则化及其在偏移成像中的应用[J].石油地球物理勘探,2017,52(1):13-19.CHEN Xiaochun,CHEN Hui,YU Qin,et al.Seismic data interpolation with anti-aliasing POCS method and its application in seismic migration imaging[J].Oil Geophysical Prospecting,2017,52(1):13-19.
[16]
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.
[17]
曾小牛,李夕海,刘继昊,等.基于凸集投影的重力数据扩充下延一体化方法[J].石油地球物理勘探,2019,54(5):1166-1173.ZENG Xiaoniu,LI Xihai,LIU Jihao,et al.An integration of interpolation,edge padding,and downward continuation for gravity data based on the projection onto convex sets[J].Oil Geophysical Prospecting,2019,54(5):1166-1173.
[18]
Ji S X,Wang Y F,Zou A Q.Regularizing inversion of susceptibility with projection onto convex set using full tensor magnetic gradient data[J].Inverse Problems in Science and Engineering,2017,25(2):202-217.
[19]
Spector A,Grant F S.Statistical models for interpreting aeromagnetic data[J].Geophysics,1970,35(2):293-302.
[20]
Ravat D,Pignatelli A,Nicolosi I,et al.A study of spectral methods of estimating the depth to the bottom of magnetic sources from near-surface magnetic anomaly data[J].Geophysical Journal International,2007,169(2):421-434.
[21]
曾小牛,李夕海,陈鼎新,等.基于径向谱的位场向下延拓正则参数选取方法[J].石油地球物理勘探,2015,50(4):749-754.ZENG Xiaoniu,LI Xihai,CHEN Dingxin,et al.New regularization parameter selection for potential field downward continuation based on the radial spectrum[J].Oil Geophysical Prospecting,2015,50(4):749-754.
[22]
Chen G,Cheng Q,Zhang H.Matched filtering method for separating magnetic anomaly using fractal model[J].Computers & Geosciences,2016,90(1):179-188.
[23]
Cooper G R J,Cowan D R.Enhancing potential field data using filters based on the local phase[J].Computers & Geosciences,2006,32(3):1585-1591.
