有限元法弹性波反演

ELASTIC WAVE INVERSION BY FINITE ELEMENT METHOD

大连理工大学,大连116023

Dalian University of Technology, Dalian 116023

本文在弹性波传播的有限元解法的基础上,通过求残差最小和引入实际地质构造的约束条件,将弹性地震波反演转化二次规划问题;应用Lcmke算法求解二次规划问题,并在新的初始背景下重复上述过程,即可得到弹性波反演问题的解.本文的方法能灵活应用于各种复杂地质构造,对有大量约束条件的反问题计算效率较高.

With the finite element method in solving clastic wave propagation problems, an elastic waveform inversion problem is transformed into a quadratic programming one by minimizing residuals and introducing constrains on stratigraphic structures, The quadratic programming is then iterativcly solved using the Lemke algorithm with a new background value for each iteration. In this way, the solution of the clastic waveform inversion is finally obtained. The method presented here can be appliad to various complex geological and stratigraphic structures in a flexible manner, and it is efficient for inverse problems with a lot of constrains.