完全匹配层人工边界条件中的衰减因子分析
|
摘要
为了得到地下地质体的地震波场响应特征,需要进行三维地震波场数值模拟,但其计算量很大。利用有限差分方法进行数值模拟时,减少计算量的途径有两种:一是差分格式设计;另一种是边界条件的选择。从完全匹配层(PML)吸收边界条件的原理出发,总结了目前PML人工边界条件中的常用衰减函数,说明了衰减函数设计的重要性。数值模拟结果表明,余弦型吸收衰减因子的波场模拟效果要优于其它吸收衰减因子。
3D seismic numerical simulation can obtain the more realistic wavefield response of subsurface geologic bodies,but the computation amount is huge.During numerical simulation by using finite difference method,there are two ways to reduce computation amount, one is difference form design,the other is the selection of boundary condition.Starting from the principles of Perfect Matched Layer(PML) absorption boundary condition,the conventional attenuation functions of the PML artificial boundary conditions were summarized,and the importance of attenuation design was elaborated. Numerical simulation results indicate that the wavefield simulation result of cosine absorption is better than other absorption attenuation factors.
3D seismic numerical simulation can obtain the more realistic wavefield response of subsurface geologic bodies,but the computation amount is huge.During numerical simulation by using finite difference method,there are two ways to reduce computation amount, one is difference form design,the other is the selection of boundary condition.Starting from the principles of Perfect Matched Layer(PML) absorption boundary condition,the conventional attenuation functions of the PML artificial boundary conditions were summarized,and the importance of attenuation design was elaborated. Numerical simulation results indicate that the wavefield simulation result of cosine absorption is better than other absorption attenuation factors.
引文
[1]贺振华.反射地震资料偏移处理与反演方法[M].重庆:重庆大学出版社,1989:1-2
[2]Cerjan C,Kosloff D.A none reflection boundary condition for discrete acoustic and elastic wave equation ??[J].Geophysics,1985,50(4):705-708
[3]Clayton R,Engquist B.Absorbing boundary conditions for acoustic and elastic wave equations[J]. Bull Seismic Society America,1977,67(6 ): 1529-1540
[4]Berenger J P.Three-dimensional perfectly matched layer for the absorption of electromagnetic waves [J].Journal of Computational Physics,1996,127(2): 363-379
[5]王守东.声波方程完全匹配层吸收边界条件[J].石油地球物理勘探,2003,38(1):31-34
[6]Collino F,Tsogka C.Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J].Geophysics, 2001,66(1):294-307
[7]Groby J,Tsogka C.A time domain method for modeling viscoacoustic wave propagation[J].Proceedings of the 6~(th) International Conference on Mathematical and Numerical Aspects of Wave Propagation,2003, 631-636
[8]李宾.横向各向同性介质有限差分法波场模拟方法研究[D].北京:中国石油大学,2009
[9]陈可洋.完全匹配层吸收边界条件研究[J].石油物探,2010,49(5):473-477
[2]Cerjan C,Kosloff D.A none reflection boundary condition for discrete acoustic and elastic wave equation ??[J].Geophysics,1985,50(4):705-708
[3]Clayton R,Engquist B.Absorbing boundary conditions for acoustic and elastic wave equations[J]. Bull Seismic Society America,1977,67(6 ): 1529-1540
[4]Berenger J P.Three-dimensional perfectly matched layer for the absorption of electromagnetic waves [J].Journal of Computational Physics,1996,127(2): 363-379
[5]王守东.声波方程完全匹配层吸收边界条件[J].石油地球物理勘探,2003,38(1):31-34
[6]Collino F,Tsogka C.Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J].Geophysics, 2001,66(1):294-307
[7]Groby J,Tsogka C.A time domain method for modeling viscoacoustic wave propagation[J].Proceedings of the 6~(th) International Conference on Mathematical and Numerical Aspects of Wave Propagation,2003, 631-636
[8]李宾.横向各向同性介质有限差分法波场模拟方法研究[D].北京:中国石油大学,2009
[9]陈可洋.完全匹配层吸收边界条件研究[J].石油物探,2010,49(5):473-477
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心