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

Abstract

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

CLC Number:

P315.7

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.

References

[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.

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

PDF (PC)

Knowledge

Abstract

Cite this article

share this article

References

Related Articles 15

Recommended Articles

Metrics

Comments

[1]	WANG Hui, CAO Jian-ling, YAO Qi, SHI Yu-tao, LIU Yue. Two-dimensional Numerical Simulations of Effects of Static Stress Disturbance on Seismic Rupture Propagation Near Stepovers of Striking-slip Fault [J]. EARTHQUAKE, 2024, 44(2): 86-103.
[2]	LI Yuan, LIU Xi-kang, LIU Xia, DU Xue-song, WAN Yong-kui. Numerical Simulation of the Cross-fault Baseline Anomalies of Trend Break and Annual Cycle Disappearance at the Tangshan Station [J]. EARTHQUAKE, 2019, 39(1): 39-47.
[3]	ZHANG Guo-ling,JIA Li-feng,QIAO Zi-yun,GUO Xue-zeng,NIU Shu-yu. Quantitative Analysis of Georesistivity Disturbances at Xingji Station, Hebei Province [J]. EARTHQUAKE, 2017, 37(4): 102-111.
[4]	SHI Yu-tao, LIU Lan-bo, GAO Yuan. Constructing Anisotropic Model from Directional Heterogeneity Distribution [J]. EARTHQUAKE, 2015, 35(2): 1-10.
[5]	LIU Cui, SHI Yao-lin, QIAO Yan-chao, DENG Jin-fu, LI Ning, DUAN Pei-xin. Numeric Simulation of Early Jurassic Magma-heating Supply in the Yanshan Area [J]. EARTHQUAKE, 2013, 33(4): 257-268.
[6]	JING Hui-min, ZHANG Huai, WU Zhong-liang, XUN Yang, WANG Ji, SHI Yao-lin. Tsunami Constraints on Finite Fault Models: the March 11, 2011 Tohoku Earthquake [J]. EARTHQUAKE, 2013, 33(4): 207-213.
[7]	ZHANG Si-qi, ZHANG Huai, SHI Yao-lin. A Study on Anisotropy of Solid-tidal Strain [J]. EARTHQUAKE, 2013, 33(4): 190-195.
[8]	ZHENG Liang, ZHANG Huai, SUN Yu-jun, CHENG Hui-hong, ZHANG Bei, SHI Yao-lin. Comparison of 2D and 3D Finite Element Simulations in the Research of Reservoir Induced Earthquakes [J]. EARTHQUAKE, 2013, 33(4): 162-171.
[9]	YAN Zhen-zhen, ZHANG Huai, FAN Xiang-tao, DU Xiao-ping, SHI Yao-lin. Numerical Simulation of River Evolution Processes under Strike-slip Faulting [J]. EARTHQUAKE, 2013, 33(4): 105-114.
[10]	CAO Jian-ling, ZHANG Jing, WANG Hui, FANG Yin. Numrical Simulation of Fault Deformation in the Capital Region of China [J]. EARTHQUAKE, 2013, 33(3): 116-123.
[11]	YANG Xing-yue, CHEN Lian-wang, YANG Li-ming, LI Yu-jiang. Effect on Seismic Situation of Gansu Region from Several Strong Earthquakes around the Bayan Har Block [J]. EARTHQUAKE, 2013, 33(3): 77-89.
[12]	GU Shen-yi, LIU Bo-wei, ZHANG Hui, LIU Yang, XIE Xiao-jing, YE Xiang-ding. Numerical Simulation of Coseismic Response and Mechanism of Borehole Water Temperature to Far-field Strong Earthquakes [J]. EARTHQUAKE, 2013, 33(1): 29-39.
[13]	ZHANG Hui, JIAO Ming-ruo, LIU Xia. Numerical Simulations of the Influencing Mechanism of the Pacific Plate Subduction to NE China on Deep and Shallow Earthquakes [J]. EARTHQUAKE, 2012, 32(2): 135-144.
[14]	ZOU Zhen-xuan, FU Jian-fu, ZHU Yuan-qing, LI Jun. Tsunami Simulation for the Coast Area of East China [J]. EARTHQUAKE, 2011, 31(4): 118-124.
[15]	CAO Jian-ling, ZHANG Jing, WANG Hui, FANG Yin. Fault Activities and Their Geodynamic Cause Derived from Cross-fault Observation in the Capital Region of China [J]. EARTHQUAKE, 2011, 31(4): 77-85.