Abstract:Finite element partial differential numerical solution is a high accu-racy method for solving wave equation; however, this method perhaps can not be widely used in seismic migration because of its too large quantity of operation. High accuracy migration result can be got fastly so long as finite elements are used in x direction and difference operations are made in y, t directions, that is to say, y, t are considered as parameters, and x the variable,variationaperatidn being conducted in x direction.The research range in x direction is divided into 2N equaI intervals, the even number points being on partition marks, and odd number points on interpolating paints.There are five nonzero elements in even columns and three nonzero elements in odd columns so that two positive definite symmetric matrixes A, B may be solved to accomplish the solution of finite element wave equation in x direction.Taking difference operation on (t, t) plane can bring the equation of quindiagonal positive definite symmetric matrix, The unknown migration wave field can be dedued by conducting LR decomposition of the matrix, This migration method has two advantages:(I)the higher accuracy due to the use of curve interpolation;(2) the Iower rounding error because of the use of quindiagonal symmetric matrix.Unfortunately, this method needs more operation,its computer time is I,6 multiples of that taken by Claerbout's finite difference method.
收稿日期: 1985-07-01
引用本文:
孙彻, 何柏荣, 杨成新, 刘企英. 波动方程有限元—差分法数值解[J]. 石油地球物理勘探, 1986, 21(3): 259-267.
Sun Che. Finite element difference numerical solution of wave equation. OGP, 1986, 21(3): 259-267.