In the magnetic exploration theory, total-field anomaly ΔT is regarded as the component Tap of the magnetic anomaly vector Ta on the main field (T0) direction and thus constitutes the theoretical basis. However, there is an error in this approximation. Theoretical calculations and experiments have proved that this approximation error will increase rapidly as the Ta increases. When the magnetic anomaly Ta is much smaller than T0, the influence of the error is small and negligible. In the case of a strong magnetic anomaly, the error is large, and the processing interpretation accuracy of the ΔT anomaly is greatly affected. For high-precision magnetic exploration, ΔT must be converted to a magnetic anomaly component Tap for processing and interpretation. In this paper, the method of accurately calculating the magnetic anomaly component using Tap based on the Limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm is proposed. Firstly, the authors derived the forward formula for ΔT from Tap, and then constructed the objective function of Tap inversion by the difference function between ΔT and Tap. L-BFGS algorithm was used to solve the Tap from ΔT. Model experiments show that the Tap calculated by this method is very close to the real value, which can reduce the error by two orders of magnitude. This method also yields good results in the presence of noise and background fields. The method was applied to the processing of ΔT magnetic survey data of the Yangshan iron mine in Fujian Province, and the results of processing and interpretation which are more consistent with the actual results were obtained.
Guan Z N . Research and development of magnetic law exploration in China[J]. Chinese Journal of Geophysics, 1997,40(1):299-307.
[3]
Blakely R J. Potential theory in gravity and magnetic applications[M]. Cambridge: Cambridge University Press, 1995.
[4]
管志宁 . 地磁场与磁力勘探[M]. 北京: 地质出版社, 2005.
[4]
Guan Z N. Geomagnetic field and magnetic exploration [M]. Beijing: Geological Press, 2005.
[5]
刘天佑 . 磁法勘探[M]. 北京: 地质出版社, 2013.
[5]
Liu T Y. Magnetic exploration[M]. Beijing: Geological Press, 2013.
[6]
Hinze P W J, Von Frese R R B, Saad A H . Gravity and magnetic exploration: principles, practices, and applications[M]. Cambridge: Cambridge University Press, 2013.
[7]
塔费耶夫 . 金属矿床的地球物理勘探文选[M]. 北京: 地质出版社. 1954.
[7]
Tafeev. A geophysical exploration of metal deposits[M]. Beijing: Geological Publishing House. 1954.
Yuan X Y, Yao C L, Zheng Y M , et al. Error analysis of ΔT anomaly calculation of ferromagnetic body[J]. Chinese Journal of Geophysics, 2015,58(12):4756-4765.
Yuan X Y . Research on calculation error and high-precision processing conversion analysis of strong magnetic anomaly ΔT [D]. China University of Geosciences ( Beijing), 2016.
[10]
Baranov W A . New method for interpretation of aeromagnetic maps: pseudo-gravimetric anomalies[J]. Geophysics, 1957,22(2):359-383.
doi: 10.1190/1.1438369
[11]
Bhattacharyya B K . Two-dimensional harmonic analysis as a tool for magnetic interpretation[J]. Geophysics, 1965,30(5):829-857.
doi: 10.1190/1.1439658
[12]
袁亚湘, 孙文瑜 . 最优化理论与方法[M]. 北京: 科学出版社, 1997.
[12]
Yuan Y X, Sun W Y. Optimization theory and method[M]. Beijing: Science Press, 1997.
[13]
Sun W, Yuan Y X . Optimization theory and methods[M]. New York: Springer US, 2006: 175-200.
[14]
姚姚 . 地球物理反演[M]. 北京: 中国地质大学出版社, 2002.
[14]
Yao Y. Geophysical inversion[M]. Beijing: China University of Geosciences Press, 2002.
[15]
王家映 . 地球物理反演理论[M]. 北京: 高等教育出版社, 2002.
[15]
Wang J Y , Geophysical inversion theory[M]. Beijing: Higher Education Press, 2002.
[16]
Avdeeva A, Avdeev D . A limited-memory quasi-Newton inversion for 1D magnetotellurics[J]. Geophysics, 2006,71(5):G191-G196.
doi: 10.1190/1.2236381
Zhang S Q, Liu C C, Han L G , et al. Frequency multi-scale full waveform inversion based on L-BFGS algorithm and simultaneous excitation source[J]. Journal of Jilin University: Earth Science Edition, 2013,43(3):1004.
[18]
Broyden C G . The convergence of a class of double-rank minimization algorithms[J]. Journal of the Institute of Mathematics and Its Applications, 1970(6):76-90.
[19]
Fletcher R . A new approach to variable metric algorithms[J]. Computer Journal, 1970,13(3):317-322.
doi: 10.1093/comjnl/13.3.317
[20]
Goldfarb D . A family of variable metric updates derived by variational means[J]. Mathematics of Computation, 1970,24(109):23-26.
doi: 10.1090/S0025-5718-1970-0258249-6
[21]
Shanno D F . Conditioning of quasi-Newton methods for function minimization[J]. Mathematics of Computation, 1970,24(111):647-656.
doi: 10.1090/S0025-5718-1970-0274029-X
[22]
Zheng Y K, Wang Y, Chang X . Wave-equation traveltime inversion: comparison of three numerical optimization methods[J]. Computers & Geosciences, 2013,60(10):88-97.
[23]
Fabien-Ouellet G, Gloaguen E, Giroux B . A stochastic L-BFGS approach for full-waveform inversion[J]. Seg Technical Program Expanded, 2017: 1622-1627.
[24]
Dong C, Liu J N . On the limited memory BFGS method for large scale optimization[J]. Mathematicial Programming, 1989,45(1-3):503-528.
doi: 10.1007/BF01589116
[25]
Nocedal J . Upgrading quasi-newton matrices with limited storage[J]. Mathematics of Computation, 1980,35(151):773-782.
doi: 10.1090/S0025-5718-1980-0572855-7
[26]
Nocedal J, Wright S J . Numerical optimization[M]. New York: Springer, 2006.
[27]
Pan W, Innanen K A, Liao W . Accelerating hessian-free gauss-newton full-waveform inversion via L-BFGS preconditioned conjugate-gradient algorithm[J]. Geophysics, 2017,82(2):R49-R64.
doi: 10.1190/geo2015-0595.1