Edge recognition using joint correlation coefficient of gravity gradient
HOU Zhenlong1,2, YUAN Yuan3, WANG Ende1,2, FU Jianfei1,2, ZHENG Yujun1
1. Key Laboratory of Ministry of Education on Safe Mining of Deep Metal Mines, Northeastern University, Shenyang, Liaoning 110819, China; 2. College of Resources and Civil Engineering, Northeastern University, Shenyang, Liaoning 110819, China; 3. Key Laboratory of Submarine Geosciences, Se-cond Institute of Oceanography, Ministry of Natural Resources, Hangzhou, Zhejiang 310012, China
Abstract:To improve the quality of gravity gradient ima-ging and resolution of identifying subsurface structure boundaries,we use the joint correlation coefficient of total horizontal derivative in x and y direction with gravity gradient for edge recognition.Correlation coefficients calculated at multi-directional windows are integrated to enhance image.Binarization and thresholding are also utilized to highlight edges.Model tests show that model boundaries,especially the corners,could be clearly imaged.We apply this method to the aerial data observed at the Bay St.George's,Canada and identify major subsurface structure boundaries.The tests demonstrate the feasibility of the method.We also discuss the impacts of window size,observation point spacing,and line spacing on imaging.Compared with other methods,the method is theoretically simple and can yield better image.
商宇航,邰振华,秦涛.基于双曲线密度模型的频率域界面反演[J].石油地球物理勘探,2018,53(4):858-864.SHANG Yuhang,TAI Zhenhua,QIN Tao.Interface inversion in the frequency domain based on the hyperbolic density model[J].Oil Geophysical Prospecting,2018,53(4):858-864.
[2]
林宝泽,肖锋,王明常.二维重力数据径向反演及应用[J].石油地球物理勘探,2018,53(2):403-409.LIN Baoze,XIAO Feng,WANG Mingchang.2D gra-vity data radial inversion[J].Oil Geophysical Prospecting,2018,53(2):403-409.
[3]
Cordell,L,Grauch V J S.Mapping basement magnetization zones from aeromagnetic data in the San Juan Basin,New Mexico[C].SEG Technical Program Expanded Abstracts,1982,1:246-247.
[4]
Roest W R,Verhoef J,Pilkington M.Magnetic interpretation using the 3-D analytic signal[J].Geophy-sics,1992,57(1):116-125.
[5]
Miller H G,Singh V.Potential field tilt:a new concept for location of potential field sources[J].Journal of Applied Geophysics,1994,32(2-3):213-217.
[6]
Verduzco B,Fairhead J D,Green C M,et al.New insights into magnetic derivatives for structural mapping[J].The Leading Edge,2004,23(2):116-119.
[7]
Wijns C,Perez C,Kowalczyk P.Theta map:edge detection in magnetic data[J].Geophysics,2005,70(4):L39-L43.
[8]
Sertcelik I,Kafadar O.Application of edge detection to potential field data using eigenvalue analysis of structure tensor[J].Journal of Applied Geophysics,2012,84:86-94.
[9]
朱自强,曾思红,鲁光银.重力张量数据的目标体边缘检测方法探讨[J].石油地球物理勘探,2011,46(3):482-488.ZHU Ziqiang,ZENG Sihong,LU Guangyin.Discussion on the edge detection technique for gravity gradient tensor data[J].Oil Geophysical Prospecting,2011,46(3):482-488.
[10]
Yuan Y,Geng M.Directional total horizontal derivatives of gravity gradient tensor and their application to delineat the edges[C].Extended Abstracts of 76th EAGE Conference and Exhibition, 2014,1-3.
[11]
袁园,黄大年,余青露.利用加强水平方向总水平导数对位场全张量数据进行边界识别[J].地球物理学报,2015,58(7):2556-2565.YUAN Yuan,HUANG Danian,YU Qinglu.Using enhanced directional total horizontal derivatives to detect the edges of potential-field full tensor data[J].Chinese Journal of Geophysics,2015,58(7):2556-2565.
[12]
Yuan Y,Huang D,Yu Q,et al.Edge detection of potential field data with improved structure tensor methods[J].Journal of Applied Geophysics,2014,108:35-42.
[13]
MaG,LiuC,HuangD.Theremovalofadditional edges in the edge detection of potential field data[J].Journal of Applied Geophysics,2015,114:168-173.
[14]
Zuo B,Hu X.Edge detection of gravity field using eigenvalue analysis of gravity gradient tensor[J].Journal of Applied Geophysics,2015,114:263-270.
[15]
Yuan Y,Gao J Y,Chen L N.Advantages of horizontal directional theta method to detect the edges of full tensor gravity gradient data[J].Journal of Applied Geophysics,2016,130:53-61.
[16]
Ma G,Huang D,Liu C.Step-edge detection filters for the interpretation of potential field data[J].Pure and Applied Geophysics,2016,173(3):795-803.
[17]
周文纳,杜晓娟,李吉焱.重力数据的平面全张量梯度角度边界识别方法[J].石油地球物理勘探,2013,48(6):1009-1015.ZHOU Wenna,DU Xiaojuan,LI Jiyan.Gravity data edge identification using plane full tensor gradient angle[J].Oil Geophysical Prospecting,2013,48(6):1009-1015.
[18]
邰振华,张凤旭,刘国兴,等.基于差分求导的伪总梯度余弦值法的位场边界检测[J].石油地球物理勘探,2014,49(4):807-814.TAI Zhenhua,ZHANG Fengxu,LIU Guoxing,et al.Potential field edge detection based on difference derivative with cosine of pseudo total gradient[J].Oil Geophysical Prospecting,2014,49(4):807-814.
[19]
Marson I,Klingele E E.Advantages of using the vertical gradient of gravity for 3-D interpretation[J].Geo-physics,1993,58(11):1588-1595.
[20]
Guo L,Meng X,Shi L.3D correlation imaging of the vertical gradient of gravity data[J].Journal of Geophysics and Engineering,2011,8(1):6-12.
[21]
Guo L,Shi L,Meng X.3D correlation imaging of magnetic total field anomaly and its vertical gradient[J].Journal of Geophysics and Engineering,2011,8(2):287-293.
[22]
孟小红,刘国峰,陈召曦,等.基于剩余异常相关成像的重磁物性反演方法[J].地球物理学报,2012,55(1):304-309.MENG Xiaohong,LIU Guofeng,CHEN Zhaoxi,et al.3-D gravity and magnetic inversion for physical pro-perties based on residual anomaly correlation[J].Chinese Journal of Geophysics,2012,55(1):304-309.
[23]
Hou Z,Huang D,Wei X.Fast inversion of probability tomography with gravity gradiometry data based on hybrid parallel programming[J].Journal of Applied Geophysics,2016,124:27-38.
[24]
Hou Z,Huang D.Multi-GPU parallel algorithm design and analysis for improved inversion of probability tomography with gravity gradiometry data[J].Journal of Applied Geophysics,2017,144:18-27.
[25]
马国庆,杜晓娟,李丽丽.利用水平与垂直导数的相关系数进行位场数据的边界识别[J].吉林大学学报(地球科学版),2011,41(增刊1):345-348.MA Guoqing,DU Xiaojuan,LI Lili.Edge detection of potential filed data using correlation coefficients of horizontal and vertical derivatives[J].Journal of Jilin University (Earth Science Edition),2011,41(S1):345-348.
[26]
徐梦龙,杨长保,吴燕冈,等.利用各方向均方差相关系数进行位场边界检测[J].应用地球物理,2015,12(1):23-34.XU Menglong,YANG Changbao,WU Yangang,et al.Edge detection in the potential field using the correlation coefficients of multidirectional standard deviations[J].Applied Geophysics,2015,12(1):23-34.
[27]
Rafael C,Gonzalez,Richard E W.数字图像处理(第三版)[M].北京:电子工业出版社,2017.
[28]
侯振隆.重力全张量梯度数据的并行反演算法研究及应用[D].吉林长春:吉林大学,2016.HOU Zhenlong.Research and Application of the Pa-rallel Inversion Algorithms Based on the Full Tensor Gravity Gradiometry Data[D].Jilin University,Changchun,Jilin,2016.