论文部分内容阅读
地震层析成像技术是反演地下速度结构的一种重要方法。与医学CT不同,一般地震波射线对目的区覆盖不均匀,所获得的数据不充分。在这种情况下,若初始猜测不准,常规的成像方法,如阻尼最小二乘法、ART法等难以获得理想效果。为此,我们研究了带有约束条件的成像方法。这种方法是通过加入不等式约束条件后再求解最小二乘问题的。文中得出了这种约束条件下的最优化问题的解。处理该问题的方法也适用于其它非线性约束最优化问题。最后的数值结果及解的评价证明了该方法在处理射线数据不足情况下的有效性及优越性,且有一定的抗噪能力。
Seismic tomography is an important method to invert subsurface velocity structure. Unlike medical CT, seismic wave radiation generally does not cover the target area unevenly and the data obtained are not sufficient. In this case, if the initial guess is not accurate, conventional imaging methods, such as damping least squares, ART method is difficult to obtain the desired effect. To this end, we study the imaging method with constraints. This method solves the least-squares problem by adding inequality constraints. In this paper, we get the solution of the optimization problem under such constraints. The method of dealing with this problem is also applicable to other nonlinear constrained optimization problems. Finally, the numerical results and solution evaluations demonstrate the effectiveness and superiority of this method in dealing with the shortage of ray data, and have certain anti-noise ability.