交错网格地震波场模拟及频散校正策略
|
摘要
交错网格高阶差分解法是地震波场模拟的一种有效方法。较大的时间延拓步长可以提高效率,但可能引起算法不稳定,使模拟无法进行。因此,对于空间网格大小和时间步长的选择要适中。虽然震源子波的主频对模拟精度的影响不大,但适当提高震源子波的主频,可以提高地震波对于薄层的分辨率。FCT方法能够有效地压制数值频散,将FCT方法与交错网格有限差分法相结合,可以提高波场模拟的精度和运算效率。
The staggered-grid high-order differential method is an effective method in simulating seismic wave propagation.In this thesis,we find that a long time elapse will improve the efficiency of calculating,but may cause instability of the algorithm and make simulating fail.So the grid of space and the step of time should be selected appropriately.Although the dominant frequency of the seismic wavelet doesn't have much impact on the accuracy of simulating,enhancing it properly can improve the resolution for thin bed.Moreover the method flux-corrected transport(FCT) can suppress numerical dispersion effectively.Combining the method of FCT with staggered-grid properly can improve the accuracy of the wave field simulation and processing effectiveness.
The staggered-grid high-order differential method is an effective method in simulating seismic wave propagation.In this thesis,we find that a long time elapse will improve the efficiency of calculating,but may cause instability of the algorithm and make simulating fail.So the grid of space and the step of time should be selected appropriately.Although the dominant frequency of the seismic wavelet doesn't have much impact on the accuracy of simulating,enhancing it properly can improve the resolution for thin bed.Moreover the method flux-corrected transport(FCT) can suppress numerical dispersion effectively.Combining the method of FCT with staggered-grid properly can improve the accuracy of the wave field simulation and processing effectiveness.
引文
[1]VIRIEUX JEAN.P-SV wave propagation in heteroge-neous media:Velocity-stress finite-difference method[J].Geophysics,1986,51(4):889.
[2]STOLL R D.Sediment Acoustics[J].Springer Verlag,1989:154.
[3]GASSMAN FRITZ.Elastic Waves through a Packing ofSpheres[J].Geophysics,1951,16(4):673.
[4]AKI K I,RICHARDS P G.Quantitative seismology[M]:W.H.Freeman and Co;1980.
[5]牛滨华,孙春岩.固体弹性介质与地震波传播[M].北京:地质出版社,2005.
[6]牟永光,裴正林.三维复杂介质地震数值模拟[M].北京:石油工业出版社,2003.
[7]董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):412.
[8]董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法稳定性研究[J].地球物理学报,2000,43(6):856.
[9]董良国.弹性波数值模拟中的吸收边界条件[J].石油地球物理勘探,1999,34(1):45.
[10]陆基孟.地震勘探原理[M].北京:石油大学出版社,1993.
[2]STOLL R D.Sediment Acoustics[J].Springer Verlag,1989:154.
[3]GASSMAN FRITZ.Elastic Waves through a Packing ofSpheres[J].Geophysics,1951,16(4):673.
[4]AKI K I,RICHARDS P G.Quantitative seismology[M]:W.H.Freeman and Co;1980.
[5]牛滨华,孙春岩.固体弹性介质与地震波传播[M].北京:地质出版社,2005.
[6]牟永光,裴正林.三维复杂介质地震数值模拟[M].北京:石油工业出版社,2003.
[7]董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):412.
[8]董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法稳定性研究[J].地球物理学报,2000,43(6):856.
[9]董良国.弹性波数值模拟中的吸收边界条件[J].石油地球物理勘探,1999,34(1):45.
[10]陆基孟.地震勘探原理[M].北京:石油大学出版社,1993.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心