In full waveform inversion （FWI）, Hessian information of the misfit function is of vital importance for accelerating the convergence of the inversion; however, it usually is not feasible to directly calculate the He...In full waveform inversion （FWI）, Hessian information of the misfit function is of vital importance for accelerating the convergence of the inversion; however, it usually is not feasible to directly calculate the Hessian matrix and its inverse. Although the limited memory Broyden-Fletcher-Goldfarb-Shanno （L-BFGS） or Hessian-free inexact Newton （HFN） methods are able to use approximate Hessian information, the information they collect is limited. The two methods can be interlaced because they are able to provide Hessian information for each other; however, the performance of the hybrid iterative method is dependent on the effective switch between the two methods. We have designed a new scheme to realize the dynamic switch between the two methods based on the decrease ratio （DR） of the misfit function （objective function）, and we propose a modified hybrid iterative optimization method. In the new scheme, we compare the DR of the two methods for a given computational cost, and choose the method with a faster DR. Using these steps, the modified method always implements the most efficient method. The results of Marmousi and overthrust model testings indicate that the convergence with our modified method is significantly faster than that in the L-BFGS method with no loss of inversion quality. Moreover, our modified outperforms the enriched method by a little speedup of the convergence. It also exhibits better efficiency than the HFN method.展开更多
The sensitivity problem to mesh distortion and the low accuracy problem of the stress solutions are two inherent difficulties in the finite element method.By applying the fundamental analytical solutions (in global Ca...The sensitivity problem to mesh distortion and the low accuracy problem of the stress solutions are two inherent difficulties in the finite element method.By applying the fundamental analytical solutions (in global Cartesian coordinates) to the Airy stress function of the anisotropic materials,8-and 12-node plane quadrilateral hybrid stress-function (HS-F) elements are successfully developed based on the principle of the minimum complementary energy.Numerical results show that the present new elements exhibit much better and more robust performance in both displacement and stress solutions than those obtained from other models.They can still perform very well even when the element shapes degenerate into a triangle and a concave quadrangle.It is also demonstrated that the proposed construction procedure is an effective way for developing shape-free finite element models which can completely overcome the sensitivity problem to mesh distortion and can produce highly accurate stress solutions.展开更多
基金financially supported by the National Important and Special Project on Science and Technology(2011ZX05005-005-007HZ)the National Natural Science Foundation of China(No.41274116)
文摘In full waveform inversion （FWI）, Hessian information of the misfit function is of vital importance for accelerating the convergence of the inversion; however, it usually is not feasible to directly calculate the Hessian matrix and its inverse. Although the limited memory Broyden-Fletcher-Goldfarb-Shanno （L-BFGS） or Hessian-free inexact Newton （HFN） methods are able to use approximate Hessian information, the information they collect is limited. The two methods can be interlaced because they are able to provide Hessian information for each other; however, the performance of the hybrid iterative method is dependent on the effective switch between the two methods. We have designed a new scheme to realize the dynamic switch between the two methods based on the decrease ratio （DR） of the misfit function （objective function）, and we propose a modified hybrid iterative optimization method. In the new scheme, we compare the DR of the two methods for a given computational cost, and choose the method with a faster DR. Using these steps, the modified method always implements the most efficient method. The results of Marmousi and overthrust model testings indicate that the convergence with our modified method is significantly faster than that in the L-BFGS method with no loss of inversion quality. Moreover, our modified outperforms the enriched method by a little speedup of the convergence. It also exhibits better efficiency than the HFN method.
基金supported by the National Natural Science Foundation of China(Grant No.10872108,10876100)the Program for New Century Excellent Talents in University(Grant No. NCET-07-0477)+1 种基金the National Basic Research Program of China(Grant No. 2010CB832701)ASFC
文摘The sensitivity problem to mesh distortion and the low accuracy problem of the stress solutions are two inherent difficulties in the finite element method.By applying the fundamental analytical solutions (in global Cartesian coordinates) to the Airy stress function of the anisotropic materials,8-and 12-node plane quadrilateral hybrid stress-function (HS-F) elements are successfully developed based on the principle of the minimum complementary energy.Numerical results show that the present new elements exhibit much better and more robust performance in both displacement and stress solutions than those obtained from other models.They can still perform very well even when the element shapes degenerate into a triangle and a concave quadrangle.It is also demonstrated that the proposed construction procedure is an effective way for developing shape-free finite element models which can completely overcome the sensitivity problem to mesh distortion and can produce highly accurate stress solutions.
