摘要
利用六面体单元和三线性插值函数 ,推导出在任意弹性各向异性介质中三维三分量波动方程所满足的有限元方程 ,同时给出也适用于矩形单元和双线性插值的二维三分量各向异性介质中的波动方程模拟的一般性公式 ,最后讨论两个数值模拟的结果。
The finite element equation, which the three dimensional three component wave equation satisfies, is derived if the domain of an arbitrary anisotropic medium is divided into hexahedron elements and a linear interpolation function in each coordinate direction is used. The finite element equation is a general formula because it is also suitable for two dimensional three component wave equations.
引文
1 Crampin S and Lovell JH. A decade of shear-wavespilitting in the earth' s crust:whatdoes itmean?whatcan wemake of itand whatshould we do next?[J]Geophys JInt,1 991 ,1 0 7:3 87~ 4 0 7
2 Faria E L,Stoffa P L.Finite-difference modeling in TIM. [J]Geophysics,1 994 ,59( 2 ) :2 82~ 2 89
3 Kosloff D,etal. Propagation modeling in elastic anisotropic media.[J]SEG59Expanded Abstracts,1 989( SM1 ) :1
4 周 辉 ,何樵登 .利用 Hartley变换模拟各向异性介质地震波场 . [J]石油地球物理勘探 ,1 995,3 0 :593~ 60 2
5 Tarantola A.Inverse problem theory:methods for data fitting and model parameter estimation. [M][s.l.]:Elsevier Sciences Publishing Co,1 987
6 胡德绥 .弹性波动力学 .北京 :地质出版社 ,1 989. 1~ 2
7 Aki K,Richards P. Quantative seismology,1 . [M][s.l.]:W H Freeman and Co,1 980
8 徐世浙著 .地球物理中的有限单元法 . [M]北京 :科学出版社 ,1 994 . 76~ 77