收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

求解波动方程的龙格—库塔型方法及其地震波传播模拟

陈山  
【摘要】:地震波传播正演模拟一直是地球物理与石油勘探领域的一个研究热点。本文首先提出了求解波动方程的一种新的数值方法——Runge-Kutta(龙格-库塔)方法。该方法的基本思想是:将二阶波动方程改写为一阶偏微分方程组,在空间方向上使用网格点上的位移、粒子速度及它们的梯度值的组合来近似逼近空间高阶偏导数,使一阶偏微分方程组转化为一个半离散的常微分方程组,然后对这个半离散的常微分方程组采用三阶或四阶Runge-Kutta方法进行时间推进计算,从而获得了求解波动方程的Runge-Kutta方法。本文还针对Runge-Kutta方法存在的不足进行改进,提出了一种加权Runge-Kutta方法,并将这两种方法统称为“Runge-Kutta型方法”。 本文从理论分析和数值计算两个方面对Runge-Kutta型方法进行了系统的研究,主要包括以下研究内容:分别给出了求解波动方程的一维、二维和三维Runge-Kutta方法的具体实现步骤,并使用Fourier分析方法推导了一维、二维和三维Runge-Kutta方法的稳定性条件;对一维、二维和三维Runge-Kutta方法分别进行了误差和频散分析,并与经典的Lax-Wendroff修正方法和交错网格方法进行了数值精度和频散误差的比较;对不同介质中二维和三维的声波和弹性波的传播进行了数值模拟,并给出了模拟结果;给出了加权Runge-Kutta方法的实现步骤以及取不同加权参数时其稳定性条件;并用加权Runge-Kutta方法进行了不同介质中地震波传播的数值模拟,给出了模拟结果,同时分析了计算效率和存储效率。 理论分析和数值实验结果表明,Runge-Kutta方法具有稳定性好,数值精度高,频散误差小的优点,并且在粗网格条件下能有效压制数值频散;加权Runge-Kutta方法可以进一步节省存储空间,更加有效地压制数值频散,从而可以通过使用更大的空间和时间步长以获得更快的计算速度和更小的存储量需求。因此,Runge-Kutta方法及其加权方法将在地球物理与石油勘探领域具有巨大的应用潜力。


知网文化
【相似文献】
中国重要会议论文全文数据库 前10条
1 冯定华;李桦;丁国昊;曹文斌;;TVD Runge-Kutta时间格式在WENO5格式中的应用研究[A];第十四届全国激波与激波管学术会议论文集(下册)[C];2010年
2 蒋呈阅;颜世荣;刘劲松;王可嘉;;基于太赫兹参量过程的非线性耦合波方程组的数值求解[A];第十三届全国红外加热暨红外医学发展研讨会论文及论文摘要集[C];2011年
3 李军;陈予恕;;轴系不对中弯扭耦合共振研究[A];第十三届全国非线性振动暨第十届全国非线性动力学和运动稳定性学术会议摘要集[C];2011年
4 陈艳平;姚明辉;张伟;;高速旋转叶片的非线性动力学研究[A];第十届全国振动理论及应用学术会议论文集(2011)上册[C];2011年
5 冯峰;王强;;两介质流界面-激波相互作用RKDG方法应用分析[A];北京力学会第17届学术年会论文集[C];2011年
6 王永虎;向小军;齐永强;;起飞空中段性能计算改进算法[A];第十三届中国科协年会第22分会场-中国通用航空发展研讨会论文集[C];2011年
7 胡倩;曾大治;谢小刚;;卫星轨道积分算法性能分析[A];第三届中国卫星导航学术年会电子文集——S03精密定轨与精密定位[C];2012年
8 林伟;韩旭;刘世杰;林志勇;周进;;倾斜横向热射流起爆爆震波的二维模拟[A];第十四届全国激波与激波管学术会议论文集(上册)[C];2010年
9 李健;宁建国;;基于基元化学反应模型的H_2/O_2/N_2爆轰数值模拟[A];第十二届现代数学和力学会议论文集[C];2010年
10 王革;谢昌坦;;应用五阶WENO格式和MGFM计算多介质可压缩流[A];第十四届全国激波与激波管学术会议论文集(下册)[C];2010年
中国博士学位论文全文数据库 前10条
1 陈山;求解波动方程的龙格—库塔型方法及其地震波传播模拟[D];清华大学;2010年
2 陈浩;离散与分布型延迟系统的块边值方法[D];华中科技大学;2012年
3 覃婷婷;几类随机与延迟动力学系统的单步离散方法[D];华中科技大学;2010年
4 李计勇;二阶振荡微分方程数值方法研究[D];南京大学;2012年
5 赵国忠;变分迭代方法及Lagrange坐标系下的间断有限元方法研究[D];中国工程物理研究院;2011年
6 王小捷;随机微分方程数值算法研究[D];中南大学;2012年
7 朱为国;矩形薄板非线性热、磁弹性振动与混沌研究[D];燕山大学;2011年
8 李东方;几类常及抛物型微分方程的数值算法研究[D];华中科技大学;2011年
9 刘红良;几类非线性延迟微分代数方程的数值分析[D];湘潭大学;2010年
10 朱华君;基于小波的辛与多辛算法[D];国防科学技术大学;2011年
中国硕士学位论文全文数据库 前10条
1 金永虎;三角指数拟合Runge-Kutta方法[D];华中科技大学;2011年
2 秦军;Runge-Kutta法在求解微分方程模型中的应用[D];安徽大学;2010年
3 赵张益;一维浅水方程的Runge-Kutta间断有限元数值模拟与应用[D];天津大学;2010年
4 康明;求解刚性振荡问题的单调隐式Runge-Kutta方法[D];湘潭大学;2011年
5 巩星田;刚性延迟微分方程Runge-Kutta法处理延迟量的两类不同插值方案比较[D];湘潭大学;2011年
6 龙腾;指数Rosenbrock方法的数值分析[D];哈尔滨工业大学;2011年
7 李海峰;延迟积分微分方程多步Runge-Kutta法及并行实现[D];中南大学;2012年
8 李敏;比例方程指数Runge-Kutta方法的稳定性分析[D];哈尔滨工业大学;2011年
9 张宾;关于延迟微分方程二级Lobatto ⅢC Runge-Kutta法的若干注记[D];湘潭大学;2011年
10 陈丙振;求解振荡常微分方程的辛指数拟合Runge-Kutta(-Nystr(?)m)型方法[D];南京农业大学;2009年
中国重要报纸全文数据库 前1条
1 本报记者 杨曙光;孟山都小麦水稻育种研究国际奖学金项目启动[N];农民日报;2010年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978