资源描述：
《一阶弹性波方程交错网格高阶差分解法稳定性研究》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第43卷第6期地球物理学报Vol.43,No.62000年11月CHINESEJOURNALOFGEOPHYSICSNov.,2000[文章编号]0001-5733(2000)06-0856-09[中图分类号]P631一阶弹性波方程交错网格高阶差分解法稳定性研究董良国马在田曹景忠(同济大学教育部海洋地质重点实验室,上海200092)[摘要]稳定性问题是数值求解波动方程的基本问题.文中对三维横向各向同性介质中一阶弹性波方程交错网格高阶差分解法的稳定性进行了分析,给出了不同精度差分方程统一的稳定性条件,证明了三维TI介质中一阶弹性
2、波方程交错网格高阶差分解法的稳定性由弹性波在X、Y、Z三个方向上的Courant数共同决定.最后通过几种精度差分方程的稳定性条件,说明了这种一阶弹性波方程高阶差分解法具有高精度、高效率的特点.[关键词]横向各向同性,稳定性条件,高阶差分法,交错网格.1引言在弹性波正演模拟中,除了使用二阶方程外,常常采用一阶速度-应力弹性波方程,其主要优点是勿需对弹性常数进行空间微分.在具体的有限差分解法上,除了规则网格外,[1][2,3]一种较为先进的交错网格最早由Madariaga,R.提出,Virieux在模拟各向同性介质22中的SH波和
3、P-SV波时也使用了这种差分网格,其差分精度为O($t+$x),在不增加计算工作量和存储空间的前提下,和常规差分网格相比局部精度提高了4倍,且收敛速度[4][5]24也较快.Levander又将这种差分网格的精度提高到O($t+$x).[6]为了进一步提高差分精度,减小网格频散,本文作者将速度(应力)对时间的奇数阶高阶导数转化为应力(速度)对空间的导数,从而在不增加所需内存量的前提下,将交错网格和高阶差分法有机结合,成功地运用到求解二维TI介质一阶速度-应力弹性波方程中,并在实际模型的弹性波模拟中取得了较好的模拟效果.这种求解
4、一阶弹性波方程的交错网格高阶差分解法,与常规的有限差分法以及交错网格低阶差分法相比,网格频散显著减小,精度明显提高.在保证一定的差分精度的基础上,可以取较大时间步长和空间网格间距,提高了计算效率.由于差分计算过程中数值参数选择不合理,可能产生无物理意义的按指数增大的数值[7]计算结果,造成模拟结果网格频散严重,影响对问题的分析,严重时会造成溢出而使计算[收稿日期]1999-11-09收到,2000-05-24收到修定稿.[基金项目]教育部重点科技基金项目和海洋863-820青年基金项目.[作者简介]董良国,男,1966年生.1
5、990年毕业于同济大学勘探地球物理专业,现为该校在职博士研究生,主要从事地震波的传播以及地震数据处理等方面的教学和科研工作.6期董良国等:一阶弹性波方程交错网格高阶差分解法稳定性研究857无法进行.因此,对一种数值解法,需要知道计算稳定的离散参数区域,即分析解法的稳定性.2一阶弹性波方程的高阶差分形式2.1三维一阶弹性波方程在不考虑体力的情况下,由速度-应力表示的三维一阶弹性波方程可为5U(t)=QU(t),(1)5tT这里,U=(vx,vy,vz,Sxx,Syy,Szz,Sxy,Sxz,Syz),v为速度,Sij为应力,Q是
6、含有介质弹性常数以及空间微分算子Dx、Dy和Dz的9@9阶方阵.2.2一阶弹性波方程交错网格高阶差分形式$t$t将Ut+和Ut-通过Taylor公式展开后相减,可得2M阶时间差分精度22的差分近似M2m-12m-1$t$t($t/2)5U(t)2MUt+-Ut-=2E2m-1+O($t),(2)22m=1(2m-1)!5t其中,$t为时间步长.当M=1时,(2)式即为传统的二阶精度时间差分近似.2m-15为了提高时间差分精度,同时避免直接计算(2)式中的2m-1U(t)时所涉及过多时5t间层而需求过多的内存量,利用方程(1)可
7、以完全准确地将速度对时间的任意奇数阶高阶导数转为应力对空间的导数,将应力对时间的任意奇数阶高阶导数转为速度对空间的导数.这样,在计算一个时间层上速度(应力)时,只需要前一时间层的速度(应力)场,以及中间节点上的应力(速度)场,不需要过多的时间层,从而节省内存.由(1)式可知2m-15U(t)2m-12m-1=QU(t),(3)5t5f为了提高空间差分精度,采用交错网格法得到的2N阶空间精度的近似式为5xN5f1(N)$x$x2N=ECnf[x+(2n-1)]-f[x-(2n-1)]+O($x),(4)5x$xn=122(N)其
8、中,$x为差分网格间距,Cn为不同阶的交错网格差分系数(该系数的确定方法可参见文献[6]).将(3)式中Q精确的空间微分算子Dx、Dy和Dz用(4)式的2N阶空间精度的近似微分算子代替后,变成矩阵R,这样(3)式就变为如下2N阶空间差分精度的差分形式2m-15U(t)2m-1
