Internet Archive Scholar

Fast non-convex low-rank matrix decomposition for separation of potential field data using minimal memory

Dan Zhu, ,Institute of Geophysics & Geomatics, China University of Geosciences (Wuhan), Wuhan, MO 430074, China, Rosemary A. Renaut, Hongwei Li, Tianyou Liu, ,School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ 85287-1804, USA, ,School of Mathematics & Physics, China University of Geosciences (Wuhan), Wuhan, MO 430074, China
2020 Inverse Problems and Imaging  

Preserved Fulltext

fulltext thumbnail
Web Archive Capture PDF (9.0 MB)
https://web.archive.org/web/20210428062307/https://www.aimsciences.org/article/exportPdf?id=242f180e-bbf2-45d6-a5c6-056fca8f11d9

A copy of this work was available on the public web and has been preserved in the Wayback Machine. The capture dates from 2021; you can also visit the original URL. The file type is application/pdf.

A fast non-convex low-rank matrix decomposition method for potential field data separation is presented. The singular value decomposition of the large size trajectory matrix, which is also a block Hankel matrix, is obtained using a fast randomized singular value decomposition algorithm in which fast block Hankel matrix-vector multiplications are implemented with minimal memory storage. This fast block Hankel matrix randomized singular value decomposition algorithm is integrated into the Altproj
more » ... algorithm, which is a standard non-convex method for solving the robust principal component analysis optimization problem. The integration of this improved estimation for the partial singular value decomposition avoids the construction of the trajectory matrix in the robust principal component analysis optimization problem. Hence, gravity and magnetic data matrices of large size can be computed and potential field data separation is achieved with better computational efficiency. The presented algorithm is also robust and, hence, algorithm-dependent parameters are easily determined. The performance of the algorithm, with and without the efficient estimation of the low rank matrix, is contrasted for the separation of synthetic gravity and magnetic data matrices of different sizes. These results demonstrate that the presented algorithm is not only computationally more efficient but it is also more accurate. Moreover, it is possible to solve far larger problems. As an example, for the adopted computational environment, matrices of sizes larger than 205 × 205 generate "out of memory" exceptions without the improvement, whereas a matrix of size 2001 × 2001 can now be calculated in 1062.29s. Finally, the presented algorithm is applied to separate real gravity and magnetic data in the Tongling area, Anhui province, China. Areas which may exhibit mineralizations are inferred based on the separated anomalies.
doi:10.3934/ipi.2020076 fatcat:uozh33maszgytb5olmbvulgkem
Citation

Dan Zhu, ,Institute of Geophysics & Geomatics, China University of Geosciences (Wuhan), Wuhan, MO 430074, China, Rosemary A. Renaut, Hongwei Li, Tianyou Liu, ,School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ 85287-1804, USA, ,School of Mathematics & Physics, China University of Geosciences (Wuhan), Wuhan, MO 430074, China. "Fast non-convex low-rank matrix decomposition for separation of potential field data using minimal memory." Inverse Problems and Imaging 15.1 (2020) 1-25