Compressed FWI inversion with non-monotone line search LBFGS

Duan, C., Zhang, F., Han, L., Chang, X. and Yang, F., 2017. Journal of Seismic Exploration, 26: 561-585. Full waveform inversion (FWI) is a high quality seismic imaging method. It is a nonlinear inversion problem which usually needs the monotone line search method to be solved. However, the speed of convergence for such a simple search technique is relatively slow. In this paper, we combine the non-monotone line search technique with the LBFGS method and apply them to the frequency-domain FWI. We test this new method on a two-dimensional Marmousi model. The results show that the method is robust. Comparing with the monotone line search method, the new method could improve the convergence rate of FWI. We also test the new method with a two-dimensional conventional streamer data set and the results show some improvements compared with the conventional FWI method.
- Abdi, H. and Williams, L., 2010. Principal component analysis: Wiley
- Interdisciplinary Reviews. Computat. Statist., 2: 433-459.
- Benhadjali, H., Operto, S. and Virieux, J., 2011. An efficient frequency-domain full
- waveform inversion method using simultaneous encoded sources, Geophysics,
- 76(4): 109-124. DOE:10.1190/1.3581357.
- Brossier, R., 2011. Two-dimensional frequency-domain visco-elastic full waveform
- inversion: Parallel algorithms, optimization and performance. Comput. Geosci.,
- 37: 444-455. DOL:10.1016/j.cageo.2010.09.013.
- Broyden, C.G., 1965. A class of methods for solving nonlinear simultaneous
- equations. Mathemat. Computat., 19(92): 577-593.
- DOI:10.1090/S0025-5718-1965-0198670-6.
- Broyden, C.G., 1967. Quasi-Newton methods and their application to function
- minimization. Mathemat. Computat., 21(99): 368-381. DOI:10.2307/2003239.
- Byrd, R.H. and Nocedal, J., 1989. A tool for the analysis of quasi-Newton methods
- with application to unconstrained minimization. SIAM J. Numer. Analys., 26:
- 727-739. DOI:10.1137/0726042.
- Cheng, W.Y. and Li, D.H., 2010. Spectral scaling BFGS method. J. Optimizat.
- Theory Applicat., 146: 305-319. DOI:10.1007/s10957-010-9652-y.
- Grippo, L., Lampariello, F. and Lucidi, S., 1986. A nonmonotone line search
- technique for Newton’s Method. SIAM J. Numer. Analys., 23: 707-716.
- DOI:10.1137/0723046.
- Grippo, L., Lampariello, F. and Lucidi, S., 1989. A truncated Newton method with
- nonmonotone line search for unconstrained optimization. J. Optimizat. Theory
- Applicat., 60: 401-419. DOI:10.1007/BF00940345.
- Han, J. and Liu, G., 1997. Global convergence analysis of a new nonmonotone BFGS
- algorithm on convex objective functions. Computat. Optimizat. Applicat., 7:
- 277-289. DOI:10.1023/A:1008656711925.
- Huang, Z.H., Hu, S.L. and Han, J.Y., 2009. Convergence of a smoothing algorithm
- for symmetric cone complementarity problems with a non-monotone line search.
- Sci. China, 52: 833-848. DOI:10.1007/s11425-008-0170-4.
- Liu, C.C., Han, M., Han, L.G., Huang, F. and Deng, W., 2012. Application of
- principal component analysis for frequency-domain full waveform inversion.
- Expanded Abstr., 82nd Ann. Internat. SEG Mtg., Las Vegas: 1-5.
- DOI:10.1190/segam2012-0909.1.
- Liu, D.C. and Nocedal, J., 1989. On the Limited Memory BFGS Method for large
- scale optimization. Mathemat. Program., 45: 503-528.
- DOI:10.1007/BF01589116.
- Ma, Y.and Hale, D., 2012. Quasi Newton full-waveform inversion with a projected
- Hessian matrix. Geophysics, 77(5): 207-216. DOI:10.1190/geo02011-0519.1.
- Moore, B., 1981. Principal component analysis in linear systems: Controllability,
- observability, and model reduction. IEEE Transact. Automat. Contr., 26:
- 17-32.
- Nocedal, J., 1980. Updating quasi-Newton matrices with limited storage. Mathemat.
- Comput., 35: 773-782. DOI: 10.2307/2006193.
- Panier, E.R. and Tits, A.L., 1991. Avoiding the Maratos effect by means of a
- non-monotone linear search I , general constrained problems. SIAM J. Numer.
- Analys., 28: 1183-1195.
- Pratt, R.G., 1999. Seismic waveform inversion in the frequency domain, Part 1:
- Theory and verification in a physical scale model. Geophysics, 64: 888-901.
- DOI:10.1190/1.1444597.
- Raydan, M., 1997. The Barzilai and Borwein gradient method for the large scale
- unconstrained minimization problem. SIAM J. Optimizat., 7: 26-33.
- DOI:10.1137/S1052623494266365.
- Schiemenz, A. and Igel, H., 2013. Accelerated 3-D full-waveform inversion using
- simultaneously encoded sources in the time domain: application to Valhall
- ocean-bottom cable data. Geophys. J. Internat., 195: 1970-1988.
- DOI:10.1093/gji/ggt362.
- Shi, Z.J. and Shen, J., 2006. Convergence of non-monotone line search method. J.
- Computat. Appl. Mathemat., 193: 397-412. DOI:10.1016/j.cam. 2005. 06.003.
- Sun, W.Y., Han, J.Y., Sun, J., 2002. Global convergence of non-monotone descent
- method for unconstrained optimization problems. J. Computat. Appl. Mathemat.,
- 146: 89-98. DOI:10.1016/S0377-0427(02)00420-X.
- Tarantola, A., 1984. Inversion of seismic reflection data in the acoustic approximation.
- Geophysics, 49: 1259-1266. DOI:10.1190/1.1441754.
- Wang, Y.H. and Rao, Y., 2009. Reflection seismic waveform tomography. J. Geophys.
- Res., 114(B3): 63-73. DOI:10.1029/2008JB005916
- Zhang, H.C. and Hager, W.W., 2004. A Non-monotone Line search technique and its
- application to unconstrained optimization. SIAM J. Optimizat., 14: 1043-1056.
- DOI:10.1137/S1052623403428208.
- Zhao, H. and Gao, Z.Y., 2005. Equilibrium algorithms with non-monotone line search
- technique for solving the traffic assignment problems. J. Systems Sci.
- Complex., 18: 543-555.
- Zhu, D.T., 2003. Non-monotonic back-tracking trust region interior point algorithm
- for linear constrained optimization. J. Computat. Appl. Mathemat., 155: 285-305.
- DOI:10.1016/S0377-0427(02)00870-1.