非均匀介质散射问题的体积分方程数值解法  

VOLUME INTEGRAL EQUATION METHOD FOR SCATTERING PROBLEMS OF INHOMOGENEOUS MEDIA

作　　者：符力耘[1] 杨慧珠[1] 牟永光[2] 

机构地区：[1]清华大学工程力学系 [2]石油大学

出　　处：《力学学报》1998年第4期488-494,共7页Chinese Journal of Theoretical and Applied Mechanics

摘　　要：将非均匀介质视为某一均匀背景介质中的扰动，可建立用均匀背景介质格林函数作基本解的体积分方程．给出了配置法求解体积分方程的数值方法，首先解得扰动域内各点以速度扰动为权的波场函数，然后回代计算得到观测面上各接收点的散射波场．与边界元法和Ｂｏｒｎ近似法计算结果比较表明该方法具有很高的精度，可得到穿过非均匀介质传播的全波场，包括所有的透射波、反射波、衍射波、首波及多重散射波．最后将体积分方程法用于非均匀介质散射问题的研究和地震模型计算．It is an effective technique to treat inhomogeneous media as perturbations on the back ground of homogeneous media and to treat the perturbations as equivalent source. Two forms of corresponding integral equations are: 1) boundary integral equation with both boundary integral items for boundary displacement and volume integral items for displacement in wavefield; 2) volume integral equation including integral only for inhomogeneous bounded region. Much work has been done in this field. Interpreting inhomogeneous media as a perturbation of the reference medium one can derive a volume integral equation based on the Green's function of the homogeneous reference media. In the present paper, the collocation method is employed to solve the volume integral equation for the velocity weighted wavefield within the scattering region, which includes transmission wave, reflection wave, diffraction wave, head wave and multiscattered wave. Using Fredholm theorem, we can see that the integral equation for the velocity weighted wavefield has unique solution. Then the observed wavefield can be computed through iterated substitutions of the integral equation. Two kinds of techniques for solving the matrix form of the equation are Gaussian elimination and iterative algorithm. The former is more rigorous, but only suitable for smaller scale problems, while the latter is faster and suitable for larger scale problems. In our numerical examples, the volume integral equation method is applied to seismic modeling, with a comparison with those given by the BE method and the approach using the Born approximation, which shows that the volume integral equation is a numerical method with high accuracy for inhomogeneous media.

关 键 词：非均匀介质 散射问题 体积分方程法 地震模拟 

分 类 号：O347[理学—固体力学] P315.3[理学—力学]