The truncated Newton (TRN) method has been demonstrated as an important optimization method for multiparameter full waveform inversion (FWI). In this study, we propose an implementation of this strategy for timedomain FWI based on the 1st-order VTI viscoacoustic wave system. Compared with previous frequency-domain FWI studies with 2nd-order wave equation using TRN, the time-domain counterpart allowing for anisotropy and attenation is appealing due to widely available well-developed preprocessing tools, straightforward data windowing and the simultaneous inversion within broad frequency bands. To overcome the computational challenge for efficiently building the gradient and Hessian-vector product in the presence of seismic attenuation, we use wavefield reconstruction technique with checkpointing-assisted reverse-forward simulation (CARFS) algorithm. The TRN method equipped with the pseudo Hessian preconditioner and unity-based normalization, allows fast convergence when updating multiple parameters. The 2D Valhall synthetic study demonstrates that preconditioned TRN (PTRN) significantly improves the multiparameter reconstruction under (Vp,rho, 1/Q) parameterization, while mitigating the crosstalk issue from parameters of different classes. Moving to field data application using PTRN is an attractive direction in the future.


Article metrics loading...

Loading full text...

Full text loading...


  1. Fichtner, A. and Trampert, J.
    [2011] Hessian kernels of seismic data functionals based upon adjoint techniques. Geophysical Journal International, 185(2), 775–798.
    [Google Scholar]
  2. Gardner, G.H.F., Gardner, L.W. and Gregory, A.R.
    [1974] Formation velocity and density—The diagnostic basics for strati-graphic traps. Geophysics, 39, 770–780.
    [Google Scholar]
  3. Innanen, K.A.
    [2014] Seismic AVO and the inverse Hessian in precritical reflection full waveform inversion. Geophysical Journal International, 199(2), 717–734.
    [Google Scholar]
  4. Korta, N., Fichtner, A. and Sallarcs, V.
    [2013] Block-diagonal approximate hessian for preconditioning in full waveform inversion. In: Expanged Abstracts, 75th EAGE Conference & Exhibition incorporating SPE EUROPEC 2013.
    [Google Scholar]
  5. Métivier, L. and Brossier, R.
    [2016] The SEISCOPE Optimization Toolbox: A large-scale nonlinear optimization library based on reverse communication. Geophysics, 81(2), F11–F25.
    [Google Scholar]
  6. Métivier, L., Brossier, R., Operto, S. and Virieux, J.
    [2015] Acoustic multi-parameter FWI for the reconstruction of P-wave velocity, density and attenuation: preconditioned truncated newton approach. In: Expanded Abstracts, 85th Annual Meeting. SEG, 1198–1203.
    [Google Scholar]
  7. Métivier, L., Brossier, R., Virieux, J. and Operto, S.
    [2013] Full Waveform Inversion and the truncated Newton method. SIAM Journal On Scientific Computing, 35(2), B401–B437.
    [Google Scholar]
  8. Nocedal, J. and Wright, S.J.
    [2006] Numerical Optimization. Springer, 2nd edn.
    [Google Scholar]
  9. Pratt, R.G., Shin, C. and Hicks, G.J.
    [1998] Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion. Geophysical Journal International, 133, 341–362.
    [Google Scholar]
  10. Wang, Y., Dong, L., Liu, Y. and Yang, J.
    [2016] 2D frequency-domain elastic full-waveform inversion using the block-diagonal pseudo-Hessian approximation. Geophysics, 81(5), R247–R259.
    [Google Scholar]
  11. Yang, P., Brossier, R., Métivier, L. and Virieux, J.
    [2016a] A review on the systematic formulation of 3D multiparameter full waveform inversion in viscoelastic medium. Geophysical Journal International, 207(1), 129–149.
    [Google Scholar]
  12. [2016b] Wavefield reconstruction in attenuating media: A Checkpointing-assisted reverse-forward simulation method. Geophysics, 81(6), R349–R362.
    [Google Scholar]

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error