伪谱算法在声波数值模拟中的数值频散分析

地震 ›› 2017, Vol. 37 ›› Issue (2): 135-146.

伪谱算法在声波数值模拟中的数值频散分析

武晔,赵晓燕,石砚斌,俞子钊,何欣娟,胡泊,熊仲华

防灾科技学院, 河北燕郊 101601

收稿日期:2016-11-08 出版日期:2017-04-20 发布日期:2019-08-14
通讯作者: 赵晓燕, 副教授。 E-mail: Zxy@cidp.edu.cn。
作者简介:武晔(1975-), 女, 内蒙集宁人, 目前主要从事地震波正演模拟、数据处理方面的研究工作。
基金资助:
防灾科技学院基金项目(JY2015B11、 JPJS2016003)联合资助

Dispersion Analysis Using Pseudo-spectral Algorithm in Acoustic Wave Numerical Simulation

WU Ye,ZHAO Xiao-yan,SHI Yan-bin,YU Zi-zhao,HE Xin-juan, HU Po,XIONG Zhong-hua

Institute of Disaster Prevention Science And Technology, Hebei Sanhe 101601, China

Received:2016-11-08 Online:2017-04-20 Published:2019-08-14

摘要/Abstract

摘要： 本文推导了声波方程频散函数, 分析了伪谱法的空间网格大小和采样周期对数值频散的影响, 通过数值模拟实验得到了最佳空间参数选择方法。结果表明: 伪谱方法稳定数值模拟的最大空间采样间距选取原则是使中波长(奈奎斯特频率的一半)的采样点数为2个; 对于所有维度, 稳定性随空间采样间距的增加而增加, 但不易变化太大, 变化太大时需要适当减小震源子波的主频, 以满足空间合理采样; 空间采样间距的大小设置, 需要考虑满足采样定理和稳定性计算条件, 并且稳定性条件对空间采样间距的要求更加严格; 伪谱法数值模拟的最佳(数值频散最小)空间参数选择为中波长2个采样点, 对应主波长约6~7个采样点。以上研究对于采用伪谱法进行声波方程数值模拟过程中, 如何合理选择模拟参数提供一些参考。

关键词: 数值频散, 伪谱法, 数值模拟, 声波

Abstract: The dispersion function of acoustic wave equation for pseudo-spectral algorithm is derived. And, selection principle of the spatial sampling space for the pseudo spectral method and the influence of the spatial sampling interval and the time sampling period on the numerical dispersion are analyzed. Numerical experiments indicates that the selection principle of the spatial sampling space for stable numerical simulation is that the sampling points of the medium wavelength which frequency is half of Nyquist-frequency is 2; for all dimensions, the stability increases with the increase of the space sampling space, but the increment can not be too large; when setting the value of the space sampling interval, it is needed to consider the sampling theorem of the time domain and the stability condition, and the stability condition is more strict to the grid size; the optimal(smaller numerical dispersion)spatial grid number for the numerical simulation of pseudo spectral method is 2 sampling points for the medium wavelength, corresponding to about 6~7 sampling points for the main wavelength. The above research provides some reference for how to select the reasonable simulation parameters using pseudo spectral method.

Key words: Numerical simulation, Numerical dispersion, Pseudo-spectral method

中图分类号:

P315.7

武晔,赵晓燕,石砚斌,俞子钊,何欣娟,胡泊,熊仲华. 伪谱算法在声波数值模拟中的数值频散分析[J]. 地震, 2017, 37(2): 135-146.

WU Ye,ZHAO Xiao-yan,SHI Yan-bin,YU Zi-zhao,HE Xin-juan, HU Po,XIONG Zhong-hua. Dispersion Analysis Using Pseudo-spectral Algorithm in Acoustic Wave Numerical Simulation[J]. EARTHQUAKE, 2017, 37(2): 135-146.

参考文献

[1] Alford R M, Kelly K R, Boore D M. Accuracy of finite difference modelling of the acoustic wave equation[J]. Geophysics, 1974, 39(6): 834-842.
[2] Dablain M A. The application of high-differencing to the scalar wave equation[J]. Geophysics, 1986, 51(1): 54-66.
[3] Virieux J. P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method(shear waves)[J]. Geophysics, 1986, 51(4): 889-901.
[4] Igel H, Mora P, et al. Anisotropic wave propagation through finite-difference grid[J]. Geophysics, 1995, 60(4): 1203-1216.
[5] 侯安宁, 何樵登. 各向异性介质中弹性波记录的傅里叶变换法正演模拟[J]. 石油地球物理勘探, 1994, 29: 64-71.
[6] 董良国, 李培明. 地震波传播数值模拟中的频散问题[J]. 天然气工业, 2004, 24(6): 53-56.
[7] 李胜军, 孙成禹, 高建虎, 等. 地震波数值模拟中的频散压制方法分析[J]. 石油物探, 2008, 47(5): 444-450.
[8] 孙成禹, 宫同举, 张玉亮,等. 波动方程有限差分法中的频散与假频分析[J]. 石油地球物理勘探,2009, 44(1): 43-48.
[9] Mullen R, Belytschko T. Dispersion analysis of finite element semidiscretizations of the two-dimensional wave equstion[J]. International Journal for Numerical Methods in Engineering, 1982, 18: 11-29.
[10] Abboud N N, Pinsky P M. Finite element dispersion analysis for the three-dimensional second-order scalar wave equation[J]. Int J Numer Meth Engng, 1992, 35(6): 1183-1218.
[11] Liu J B, Sharan S K, Yao L. Wave motion and its dispersive properties in a finite element model with distortional elements[J]. Comput Struct, 1994, 52(2): 205-214.
[12] Christon M A. The influence of the mass matrix on the dispersive nature of the semi-discrete, second-order wave equation[J]. Comput methods Appl Mech Engng, 1999, 173(1): 147-166.
[13] 房营光, 莫海鸿. 有限元网格中波动的频散与稳定性的一种改进方法[J]. 地震工程与工程振动, 2000, 20(1): 21-26.
[14] 吴国忱, 王华忠.波场模拟中的数值频散分析与校正策略[J]. 地球物理学进展, 2005, 20(1): 58-65.
[15] 薛东川, 王尚旭.利用组合质量矩阵压制数值频散[J]. 石油地球物理勘探, 2008, 43(3): 318-320.
[16] 印兴耀, 周建科, 吴国忱, 等. 有限元算法在声波方程数值模拟中的频散分析[J]. 地震学报, 2014, 36(5): 944-955.
[17] De Basabe J D, Sen M K. Grid dispersion and stability criteria of some common finite-element methods for acoustic and elastic wave equations[J]. Geophysics, 2007, 72(6): T81-T95.
[18] Seriani G, Oliveira S P. Dispersion analysis of spectral element methods for elastic wave propagation[J]. Wave Motion, 2008, 45: 729-744.
[19] Fornberg B. On a Fourier method for the integration of hyper-bolic equations[J]. SIAM J. Num. Anal. 1975, 12: 509-528.
[20] Gazdag J. Modeling of the acoustic wave equation with transform methods[J]. Geophysics, 1981, 46: 854.
[21] Dan D Kosloff, Edip Baysal. Forward modeling by a Fourier method[J]. Geophysics, 1982, 40(10): 1402-1412.
[22] Dault C R, et al. A comparison of finite-difference and Fourier method calculations of synthetic seismograms[J]. Bull. Seis. Soc. Am., 1989, 79: 1210-1230.
[23] 何樵登, 张中杰. 横向各向同性介质中地震波及其数值模拟[M]. 长春: 吉林大学出版社, 1996.
[24] 侯安宁, 何樵登. VSP数据逆时偏移[J]. 长春地质学院学报, 1990, 20: 227-234.
[25] 张文生, 何樵登. 二维横向各向同性介质的伪谱法正演模拟[J]. 石油地球物理勘探[J], 1998, 33(3): 310-319.
[26] 傅旦丹, 何樵登. 正交各向异性介质地震弹性波场的伪谱法正演模拟[J]. 石油物探, 2001, 40(3): 8-14.
[27] Sheriff R E. Encyclopedic dictionary of exploration geophysics[M]. Tulsa, SEG, 1973, p.3.
[28] Hamming A R. Numerical methods for scientists and engineers[M]. New York. McGraw-Hill Book Company. Inc. 1962.