Publications
Detailed Information
Improved Hessian-Free optimization for acoustic full waveform inversion : 음향파 파형역산을 위한 개선된 헤시안-프리 최적화 기법
Cited 0 time in
Web of Science
Cited 0 time in Scopus
- Authors
- Advisor
- 민동주
- Major
- 공과대학 에너지시스템공학부
- Issue Date
- 2017-08
- Publisher
- 서울대학교 대학원
- Keywords
- seismic full waveform inversion ; truncated Newton method ; Hessian-free method ; acoustic wave equation
- Description
- 학위논문 (박사)-- 서울대학교 대학원 공과대학 에너지시스템공학부, 2017. 8. 민동주.
- Abstract
- Seismic full waveform inversion (FWI) is a method to reconstruct material properties of subsurface structures by minimizing the objective function based on residuals between modeled and observed seismic data. For seismic inverse problem, various kinds of optimization methods have been introduced. The truncated Newton method, also known as the Hessian-free (HF) optimization method, has been chosen to optimize large-scale inverse problems. The HF does not need to explicitly compute, store and invert the Hessian matrix. Instead of the Hessian matrix itself, the product of Hessian matrix and column vector is used for the linear conjugate-gradient loop during FWI process. To calculate the product of the Hessian matrix and column vector, the second-order adjoint (SOA) method or finite difference approximation (FDA) method has been widely used. The FDA is easy and intuitive to use in the linear conjugate-gradient method compared with SOA. The accuracy of FDA is dependent on not only the approximation interval but also the inversion settings, such as the model parameter, initial model, frequencies, etc.
To overcome dependency of HF optimization on the approximation method and inversion setting, an improved method is proposed for a stable HF optimization method. The derivations of the improved method are based on not the FDA method but the limit of a function, which is independent of epsilon value. In other words, the improved HF method stably and accurately approximates the matrix-vector product of the Hessian matrix and column vector without any selection of epsilon value. In addition, computational cost of the improved HF optimization method is much lower than the conventional HF optimization method because additional construction and factorization of modeling operator are not needed during the linear conjugate-gradient method in the improved HF optimization method. To demonstrate the feasibility of the improve HF method, numerical examples for the Marmousi and acoustic Overthrust models are performed. Numerical examples indicate that the improved HF method shows better computational efficiency and stability than the conventional HF method without any degradation of inversion results.
- Language
- English
- Files in This Item:
Item View & Download Count
Items in S-Space are protected by copyright, with all rights reserved, unless otherwise indicated.