粘弹介质中地震波交错网格法数值模拟.pdfVIP

粘弹介质中地震波交错网格法数值模拟 邵志刚1傅容珊’刘财2 (1中国科学技术大学地空学院，合肥，230026；2吉林大学地探学院，长春，130026) 1． 引言 虽然弹性介质是传统的模拟地震波传遍的介质模型，但实际上没有一种介质是理想的弹性体。而与以弹性波理 论为基础的理论结果相比，实际记录中的地震波存在明显的衰减。导致衰减的因素很多，至尽还没有一种极力可以 描述所有环境条件下的衰减。但有一点可以确定的是，介质的粘滞性是其中的重要因素。因此，在正演，反演以及 成像过程中不应该忽略介质粘性对波传播的影响。事实上，粘弹介质比弹性介质更接近实际的地球介质。粘弹性介 质的本构方程有微分形式和积分形式两种表达形式。微分型是根据理想模型导出的，模型能只管地模拟更适合地震 波场的正演算模拟计算。本人将采取一个折衷的方案，即用数值的方法求解以模型参数表达的积分本构方程为基础 的波动方程。 2．数值模型 个弹簧并联所表示的线性模型)，下面是一维情况下对应的波动方程 0：土塑 、(1)、LlJ V=一一 a)c p 0：彳堡+厂． (2) 苏 ；：B堡一Cr (3) ax 度，o一应力，r-ig忆变量，E“)为体积模量，E(t)为剪切模量，①f，r)=E、Ef，O)、E(O)、①f，0J 分别对应它们在t=0时的值，p为密度，k．、k，、77为模型弹性和黏性参数。 出于对计算速度和求解精度的考虑，本文使用了时间上二阶精度近似，空间上四阶精度近似的交错网格有限 差分法的模拟方案。其稳定条件是C…，譬号，C。。，为最大相速度，At、gx分别为时间步长和空间步长。 3． 结论 本文采用的记忆变量与介质参数有直接的关系，更便于从力学的角度来分析地震波的传播。本文首先从理论 上分析认为粘弹波频率域内的品质因子的最小值极其对应的频率都与介质参数有着密切的关系；在高频和低频段， 粘弹波的相速度和品质因子都趋于弹性波。本文用交错网格有限差分法模拟了=维粘性很强的粘弹介质中的地震 波的传播，模拟结果表明了粘弹波与弹性波在振幅与传播上的差异；模拟结果还表明粘弹横波要比粘弹纵波有更 大的能量衰减；所的的模拟结果与理论分析相一致。 参考文、献 wave in Virieux，J．，P-SV finite-difference propagationheterogeneousmedia：velocity-stress method。Geophysics，51，889-901，1986 finite-difference Blanch,J．O．eta1．，Viscoelastic modeling，Geophysics，59，9，1444-1456，1994 -A-本研究由中科院创新项目(KzCX3．SW．131)和国家自然科学基金赞助 ·522· 粘弹介质中地震波交错网格法数值模拟 作者： 邵志刚， 傅容珊， 刘财 作者单位： 邵志刚,傅容珊(中国科学技术大学地空学院,合肥,230026)， 刘财(吉林大学地探学院,长春 ,130026) 本文读者也读过(6条) 1. 王东.林伟军.王秀明.张海澜 有限差分数值模拟声场中自由边界的处理[会议论文]-2002 2. 刘浩 自由边界问题的数值算法研究[学位论文]2006 3. 李景叶.陈小宏.芮振华.LI Jing-ye.CHEN Xiao-hong.RUI