CAJ | 학술논문

在地震波场数值模拟中，交错网格有限差分技术得到了广泛的应用，但是在弹性模量变化较大时，通常会因插值而导致模拟误差增大。旋转交错网格可以很好地克服这个缺点，因而适合于各向异性介质正演模拟。但是对于同样大小的网格单元，旋转交错网格需要的步长比常规交错网格要大，这会使梯度和散度算子的误差增大因而更易产生空间数值频散。针对这些问题，本文提出了旋转交错网格与紧致有限差分相结合的方法，并基于模拟退火算法进行全局优化，压制数值频散，拓宽波数范围。数值模拟结果表明，此方法可以有效地压制数值频散，且具有较高的模拟精度。
재지진파장수치모의중，교착망격유한차분기술득도료엄범적응용，단시재탄성모량변화교대시，통상회인삽치이도치모의오차증대。선전교착망격가이흔호지극복저개결점，인이괄합우각향이성개질정연모의。단시대우동양대소적망격단원，선전교착망격수요적보장비상규교착망격요대，저회사제도화산도산자적오차증대인이경역산생공간수치빈산。침대저사문제，본문제출료선전교착망격여긴치유한차분상결합적방법，병기우모의퇴화산법진행전국우화，압제수치빈산，탁관파수범위。수치모의결과표명，차방법가이유효지압제수치빈산，차구유교고적모의정도。
Staggered grid finite difference techniques have been widely used in numerical simulation of seismic wave field,but the interpolation might enlarge the simulation errors when the elastic moduli change dramatically.Rotated staggered grid can overcome this drawback thus be more applicable for modeling the anisotropic media.However,for the grid cells with the same size,rotated staggered grid needs larger step than conventional staggered grid does,conse-quently,the gradient operator and divergence operator generate larger errors which are more likely to cause spatial numerical dispersion.Aiming at the above problems,a combination of rotated staggered grid and compact finite difference method is proposed,simultaneously,global optimization is preformed based on simulated annealing algorithm to suppress numerical dispersion and to broaden the range of wavenumber.Numerical simulation results show that this method can effectively suppress the numerical dispersion and exhibits high simulation accuracy at the same time.