一种四级四阶龙格—库塔时间高阶伪谱法声波方程模拟  

A four-stage fourth order Runge-Kutta time high-order pseudospectral method for acoustic equation simulation

作　　者：杨怀英[1] 唐小平[1] 刘宽厚[1] 

机构地区：[1]中国地质调查局西安地质调查中心,陕西西安710054

出　　处：《物探与化探》2014年第4期774-780,共7页Geophysical and Geochemical Exploration

摘　　要：龙格—库塔法是常用于求解常微分方程(ODE)的一项技术,该技术具有精度高、稳定性强等特点。笔者以一种新的四级四阶龙格—库塔法为基础,先将其四级计算公式合并为新的两级计算迭代公式,从而达到节约计算内存的目的;再以此为基础推导出声波方程的时间高阶离散形式,并与伪谱法技术相结合,研究四级四阶龙格—库塔时间高阶伪谱法声波高精度、高清晰度的波场模拟技术,进而研究该方法的稳定性与频散特性;最后,分别选取均匀介质、层状介质和透镜体模型进行波场模拟试验。模拟结果表明,该方法具有稳定性强、能有效去除频散,波场清晰度高,能适应更宽的模拟参数范围等特点,是一种高效率的、具有较大应用潜力的波场模拟方法。The Runge-Kutta method is a common technology for solving the Ordinary Differential Equations（ODE） and is characterized by high precision,strong stability and some other advantages. In this paper,based on a new four-stage fourth order Runge-Kutta method,the authors first combined the four-stage calculation formula with a new two-stage iteration formula,thus achieving the purpose of saving computational memory. And then,the time high-order discrete form of the acoustic wave equation was derived and,in combination with the pseudospectral method,some researches on the high-accuracy and high-definition acoustic wave field simulation technology of the four-stage fourth order Runge-Kutta time high-order pseudospectral method were carried out,with an investigation of the stability and dispersion of the method. Finally,homogeneous media,layered media and lens model were selected for wave field simulation test.The simulation results show that the four-stage fourth order Runge-Kutta time high-order pseudospectral method has strong stability and high wave field definition and can effectively remove the dispersion and adapt itself to large simulation parameter range,thus being a high efficient wave field simulation method with great application potential.

关 键 词：四级四阶龙格—库塔法 两级公式 伪谱法 声波方程 频散特性 

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