Potential field separation based on the low rank and sparse characteristics
ZHU Dan1, LIU Tianyou1, LI Hongwei2
1. Institute of Geophysics and Geomatic, China University of Geosciences(Wuhan), Wuhan, Hubei 430074, China;
2. School of Mathematics and Physic, China University of Geosciences(Wuhan), Wuhan 430074, China
Abstract:The robust principal component analysis (RPCA) extends the single objective function optimization problem for low-dimensional signal rank minimization of principal component analysis (PCA) to the two-objective function optimization problem for low-dimensional signal rank minimization with high-dimensional signal sparse.Conventional potential field separation methods such as matched filtering usually have overfitting or underfitting problems.To solve this kind of problems,this paper analyzes the low rank structures of regional field and sparse structures of local field.The RPCA is used to the potential field separation,also makes the separation more robust.Based on synthetic model tests,the proposed method is an effective and practical spatial-domain method with weight-ing parameters easy to be valued,and it can avoid errors caused by the Fourier transform.The proposed method is used to deal with data in Weiningbeishan area,Ningxia Province.Separated local low gravity anomalies and local high magnetic anomalies match very well with known underground igneous rocks,and six potential igneous rock zones are delineated.
朱丹, 刘天佑, 李宏伟. 利用数据低秩性和稀疏性的位场分离[J]. 石油地球物理勘探, 2019, 54(4): 925-936.
ZHU Dan, LIU Tianyou, LI Hongwei. Potential field separation based on the low rank and sparse characteristics. Oil Geophysical Prospecting, 2019, 54(4): 925-936.
Lin Z C,Chen M M,Wu L Q,et al.The augmented large range multiplier method for exact recovery of corrupted low-rank matrices[C].https://arxiv.org/PS_cache/arxiv/pdf/1009/1009.5055v1.pdf,2010.
[3]
Candes E J,Recht B.Exact matrix completion via convex optimization[J].Foundations of Computational Mathematics,2009,9(6):717-772.
[4]
Candes E J,Tao T.The power of convex relation:Near-optimal matrix completion[J].IEEE Transaction on Information Theory,2010,56(5):2053-2080.
[5]
Candes E J,Li X D,Ma Y,et al.Robust principal component analysis?[J].Journal of the ACM,2011,58(3):1-37.
[6]
史加荣,郑秀云,魏宗田,等.低秩矩阵恢复算法综述[J].计算机应用研究,2013,30(6):1601-1605.SHI Jiarong,ZHENG Xiuyun,WEI Zongtian,et al.Survey on algorithms of low-rank matrix recovery[J].Application Research of Computers,2013,30(6):1601-1605.
[7]
杨宇山,李媛媛,刘天佑.高阶统计量在地震弱信号及"磁亮点"识别中的应用[J].石油地球物理勘探,2005,40(1):103-107.YANG Yushan,LI Yuanyuan,LIU Tianyou.Application of high-order statistics to identify weak seismic signal and "magnetic bright spot"[J].Oil Geophysical Prospecting,2005,40(1):103-107.
[8]
马龙,孟军海,付强,等.独立分量分析在重磁数据处理中的应用[J].石油地球物理勘探,2017,52(6):1344-1353.MA Long,MENG Junhai,FU Qiang,et al.Application of independent component analysis in gravity and magnetic data processing[J].Oil Geophysical Prospecting,2017,52(6):1344-1353.
[9]
Spector A,Grant F S.Statistical model for interpreting aeromagnetic data[J].Geophysics,1970,35(2):293-302.
[10]
程方道,刘东甲,姚汝信.划分重力区域场与局部场的研究[J].物化探计算技术,1987,9(1):1-9.CHENG Fangdao,LIU Dongjia,YAO Ruxin.A study on the identification of regional and local gravity fields[J].Computing Techniques for Geophysical and Geochemical Exploration,1987,9(1):1-9.
[11]
Mickus K L,Aiken C L V,Kennedy W D.Regional-residual gravity anomaly separation using the minimum-curvature technique[J].Geophysics,1991,56(2):279-283.
[12]
王万银,邱之云,刘金兰,等.位场数据处理中的最小曲率扩边和补空方法研究[J].地球物理学进展,2009,24(4):1327-1338.WANG Wanyin,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.
[13]
纪晓琳,王万银,邱之云.最小曲率位场分离方法研究[J].地球物理学报,2015,58(3):1042-1058.JI Xiaolin,WANG Wanyin,QIU Zhiyun.The research to the minimum curvature technique for potential field data separation[J].Chinese Journal of Geophysics,2015,58(3):1042-1058.
[14]
Fedu M,Quarta T.Wavelet analysis for the regional-residual and local separation of potential field anomalies[J].Geophysical Prospecting,1998,46(5):507-525.
[15]
侯遵泽,杨文采.中国重力异常的小波变换与多尺度分析[J].地球物理学报,1997,40(1):85-95.HOU Zunze,YANG Wencai.Wavelet transform and multi-scale analysis on gravity anomalies of China[J].Chinese Journal of Geophysics,1997,40(1):85-95.
[16]
杨文采,施志群,侯遵泽,等.离散小波变换与重力异常多重分解[J].地球物理学报,2001,44(4):534-541.YANG Wencai,SHI Zhiqun,HOU Zunze,et al.Discrete wavelet transform for multiple decomposition of gravity anomalies[J].Chinese Journal of Geophysics,2001,44(4):534-541.
[17]
朱丹,刘天佑,李宏伟.基于奇异谱分析的位场分离方法[J].地球物理学报,2018,61(9):3800-3811.ZHU Dan,LIU Tianyou,LI Hongwei.Separation of potential field based on singular spectrum analysis[J].Chinese Journal of Geophysics,2018,61(9):3800-3811.
[18]
Broomhead D S,King G P.Extracting qualitative dynamics from experimental data[J].Physica D:Nonli-near Phenomena,1986,20(2):217-236.
[19]
刘天佑,崔宁,蔡鑫,等.重磁异常宽度幅值特征滤波[J].地球物理学报,1996,31(2):225-231.LIU Tianyou,CUI Ning,CAI Xin,et al.Characteristic filtering for the width and amplitude of gravimetric-magnetic anomaly[J].Chinese Journal of Geophysics,1996,31(2):225-231.
[20]
刘天佑,吴昭才,詹应林,等.磁异常小波多尺度分解及危机矿山的深部找矿:以大冶铁矿为例[J].地球科学,2007,32(1):135-140.LIU Tianyou,WU Zhaocai,ZHAN Yinglin,et al.Wavelet multi scale decomposition of magnetic ano-maly and its application in searching for deep buried minerals in crisis mines:a case study from Daye iron mines[J].Earth Science,2007,32(1):135-140.
[21]
刘天佑,盛秋红,杨宇山.复小波频谱分析在地震数据处理中的应用[J].石油地球物理勘探,2007,42(增刊1):72-75.LIU Tianyou,SHENG Qiuhong,YANG Yushan.Application of frequency spectrum analysis of complex wavelet to seismic data processing[J].Oil Geophysical Prospecting,2007,42(S1):72-75.
[22]
邓文彬,尚海滨,许闯.小波基选取对重力异常多尺度分解的影响[J].测绘科学,2018,43(4):30-37.DENG Wenbin,SHANG Haibin,XU Chuang.Effect of wavelet bases on multi-scale decomposition for gravity anomaly[J].Science of Surveying and Mapping,2018,43(4):30-37.
[23]
宋新华,尹秉喜,闫红,等.航磁资料在卫宁北山寻找多金属矿中的应用[J].物探与化探,2010,34(3):289-293.SONG Xinhua,YIN Bingxi,YAN Hong,et al.The relationship between the aeromagnetic weak anomalies and the polymetallic ore deposits (ore spots) in the Beishan mountain,Weining[J].Geophysical and Geochemical Exploration,2010,34(3):289-293.
[24]
艾宁.宁夏卫宁北山金场子金矿矿床地质与地球化学研究[D].陕西西安:西北大学,2014.AI Ning.Studies on the Geochemical and Geological Characteristics of Weiningbeishan Jingchangzi Gold Deposit[D].Northwest University,Xi'an,Shaanxi,2014.
[25]
刘志坚.宁夏卫宁北山金、铅、银多金属矿成矿地质特征[D].四川成都:成都理工大学,2013.LIU Zhijian.Metallogenic Characteristics of Weiningbeishan Au-Pb-Ag Polymetallic Deposit in Ningxia Province[D].Chengdu University of Technology,Chengdu,Sichuan,2013.
[26]
张琪,张英堂,李志宁,等.改进的低纬度化极稳定算法[J].石油地球物理勘探,2018,53(3):606-616.ZHANG Qi,ZHANG Yingtang,LI Zhining,et al.An improved stable algorithm of the reduction-to-the-pole at low latitudes[J].Oil Geophysical Prospecting,2018,53(3):606-616.