Abstract：Full-waveform inversion (FWI) obtains high-resolution parameters of subsurface models by minimizing the errors between observed data and calculated data. Regularization techniques are o-ften used to overcome the ill-posedness of FWI. However, complex geological models are likely to have sharp boundaries and smooth characteristics at the same time, in which case a single regularization method usually fails to deliver satisfactory inversion results. Given the non-differentiability of the objective function involving hybrid regularization, this paper proposes a modified Orthant-Wise Limited-memory Quasi-Newton method for solving the corresponding regularization objective function. Numerical simulation experiments are conducted on the complex-constructed modified Marmousi model and BG Compass model, and the proposed method is compared with FWI without regularization and the proximal Limited-memory Quasi-Newton method. The results show that the proposed algorithm has obvious advantages in computational efficiency and quantitative analysis.
VIRIEUX J，OPERTO S. An overview of full-waveform inversion in exploration geophysics[J]. Geophysics，2009，74(6)：WCC1-WCC26.
[2]
LAILLY P. The seismic inverse problem as a sequence of before stack migrations[C]. Theory and Applications，Society for Industrial and Applied Mathematics，1983，206-220.
[3]
TARANTOLA A. Inversion of seismic reflection data in the acoustic approximation[J]. Geophysics，1984，49(8)：1259-1266.
[4]
PRATT R G. Seismic waveform inversion in the frequency domain：Part 1，theory and verification in a physical scale model[J]. Geophysics，1999，64(3)：888-901.
[5]
张子良，李振春，张凯，等. 地质模型约束的全波形速度建模反演及在复杂断块区的应用[J]. 石油地球物理勘探，2020，55(3)：599-606.ZHANG Ziliang，LI Zhenchun，ZHANG Kai，et al. Research of geological model-constrained FWI and application in complex fault-block zones[J]. Oil Geophysical Prospecting，2020，55(3)：599-606.
[6]
杜泽源，吴国忱，王玉梅. 基于测井约束的地震全波形反演方法[J]. 石油地球物理勘探，2017，52(6)：1184-1192.DU Zeyuan，WU Guochen，WANG Yumei. Full waveform inversion based on well logging data constraint[J]. Oil Geophysical Prospecting，2017，52(6)：1184-1192.
[7]
王豆豆，王守东，邹少峰，等. 基于混合快速共轭梯度法的有限差分对比源反演[J]. 石油地球物理勘探，2020，55(2)：351-359.WANG Doudou，WANG Shoudong，ZOU Shaofeng，et al. Finite difference comparison source inversion based on hybrid fast conjugate gradient method[J]. Oil Geophysical Prospecting，2020，55(2)：351-359.
[8]
TIKHONOV A N，GONCHARSKY A V，STEPANOV V V，et al. Numerical Methods for the Solution of Ill-Posed Problems[M]. Kluwer Acade-mic Publishers，Boston，1995.
[9]
RUDIN L I，OSHER S，FATEMI E. Nonlinear total variation based noise removal algorithms[J]. Physica D：Nonlinear Phenomena，1992，60(1-4)：259-268.
[10]
ESSER E，GUASCH L，VAN LEEUWEN T，et al. Total variation regularization strategies in full-waveform inversion[J]. SIAM Journal on Imaging Sciences，2018，11(1)：376-406.
[11]
AGHAMIRY H S，GHOLAMI A，OPERTO S. ADMM-based multi-parameter wavefield reconstruction inversion in VTI acoustic media with TV regula-rization[J]. Geophysical Journal International，2019，219(2)：1316-1333.
[12]
DURAND S，FROMENT J. Reconstruction of wavelet coefficients using total variation minimization[J]. SIAM Journal on Scientific Computing，2003，24(5)：1754-1767.
[13]
BENNING M，BURGER M. Modern regularization methods for inverse problems[J]. Acta Numerica，2018，27：1-111.
[14]
LIN Y Z，HUANG L J. Acoustic- and elastic-waveform inversion using a modified total-variation regularization scheme[J]. Geophysical Journal International，2014，200(1)：489-502.
[15]
PENG Y，LIAO W Y，HUANG J P，et al. Total vriation regularization for seismic waveform inversion using an adaptive primal dual hybrid gradient method[J]. Inverse Problems，2018，34(4)：045006.
[16]
AGHAMIRY H S，GHOLAMI A，OPERTO S. Compound regularization of full-waveform inversion for imaging piecewise media[J]. IEEE Transactions on Geoscience & Remote Sensing，2019，58(2)：1192-1204.
[17]
苗永康. 基于L-BFGS算法的时间域全波形反演[J]. 石油地球物理勘探，2015，50(3)：469-474.MIAO Yongkang. Full waveform inversion in time domain based on limited-memory BFGS algorithm[J]. Oil Geophysical Prospecting，2015，50(3)：469-474.
[18]
SCHEINBERG K，TANG X C. Practical inexact proximal quasi-Newton method with global complexity analysis[J]. Mathematical Programming，2016，160(1)：495-529.
[19]
GONG P H，YE J P. A modified Orthant-Wise limited memory Quasi-Newton method with convergence analysis[C]. Proceedings of the 32nd International Conference on International Conference on Machine Learning，2015，276-284.
[20]
DAI M X，CHEN J B，CAO J. l1-Regularized full-waveform inversion with prior model information based on orthant-wise limited memory quasi-Newton method[J]. Journal of Applied Geophysics，2017，142：49-57.
[21]
ANDREW G，GAO J F. Scalable training of L1-regularized log-linear models[C]. Proceedings of the 24th International Conference on Machine Learning，2007，33-40.
