非均匀介质中地震波传播的数值模拟  

NUMERICAL SIMULATION OF SEISMIC WAVE PROPAGATION IN HETEROGENEOUS MEDIA

作　　者：刘红艳[1] 李小凡[1] 张美根[1] 

机构地区：[1]中国科学院地质与地球物理研究所岩石圈演化国家重点实验室,北京100029

出　　处：《物探化探计算技术》2008年第3期173-177,共5页Computing Techniques For Geophysical and Geochemical Exploration

基　　金：国家自然科学基金重点项目资助(40437018)

摘　　要：地震波场数值模拟方法多种多样,各种方法都有各自的特点。这里推出一种全新地震波场模拟方法——基于Forsyte广义正交多项式的迭积微分算子法,该方法将计算数学中的For-syte多项式,应用到地震波传播的数值模拟中,它同时具有广义正交多项式迭积微分算子的高精度和有限差分短算子算法的高速度。通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息。复杂非均匀介质模型数值模拟结果说明了该方法的可行性。将该方法的计算结果与傅氏变换伪谱法、错格高阶有限差分法相比较,结果说明,该方法在波场模拟方面具有较好的发展潜力,并具有自身独到的优越性。The investigation of seismic numerical simulation is an important issue in seismology. It is well known that each numerical simulation method existed today has its own merits and drawbacks. To improve the accuracy and the efficiency of seismic wave simulation, this paper develops a novel modeling approach referred to as convolutional forsyte polynomial differentiator method by using optimal convolutional operators for spatial differentiation in wave equation. This method has both the high precision of the generalized orthogonal polynomial convolutional differentiator method and the high efficiency of the finite difference short operator method. The local and global information can be considered at the same time by optimizing the coefficients of the operator and adjusting the operator length. The numerical experiments demonstrate that the algorithm can bring reliable results with high precision and can be extended to seismic wave simulation in anisotropic media. By comparing the numerical results with pseudo-spectral method and high order staggered grid finite difference method, it is shown that this method has great potential in seismic numerical simulation and has its own advantages.

关 键 词：佛尔塞多项式 褶积微分算子 地震波数值模拟 

分 类 号：P631.4[天文地球—地质矿产勘探]