方位各向异性黏弹性介质波场有限元模拟.PDFVIP

第 卷 第 期 年 月 物 理 学 报 $ ##+ $ ， ， ， :;6 ) M; ) $ GNNO= ##+ （ ） !###/+*#E##+E #$ E#!#/# G-HG 1IJKL-G KLML-G !##+ -P89 ) 1PQO ) K;R ) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 方位各向异性黏弹性介质波场有限元模拟! ）） ） 杜启振! 杨慧珠! !）（清华大学工程力学系，北京 !###$% ） ）（石油大学（华东）地球资源与信息学院，东营 ’#(!） （ 年 月 日收到； 年 月 日收到修改稿） ## ! ## ! ( 随着裂缝性储集层的地位越来越重要，黏弹性各向异性介质中地震波传播理论的研究也越来越受到重视 在 ) 弱黏滞性条件下，针对标准线性黏弹性体模型，采用时间增量法，在每一个时间步长确定黏弹性变形及载荷，并在 计算下一个时间步长的黏弹性变形时把它加载进去，建立类似于解弹性问题的黏弹性方位各向异性介质波动方程 的有限元递推公式) 数值模拟给出了黏弹性波在地球介质中的波场传播特征) 关键词：黏弹性，方位各向异性，有限元，时间增量法 ： ， ， !## (%# *!( *!+#, 程的有限元格式) 并以二维波动问题为例，给出了黏 !4 引 言 弹性波在地球介质中的波场传播特征) 随着实际地震勘探的逐步深入以及工作地区地 4 有限元解法递推格式 层的日趋复杂，复杂地层以及各向异性介质中地震 波传播的研究更加受到重视 波动方程反映了波在 ) $ % % 有限元运动方程 介质中的传播特性，能够较全面地反映波场特征，可 ［］ 以通过数值求解方法来解波动方程) 近年来，有限元 考虑黏弹性各向异性介质的本构关系 ，坐标 ［］ ( 数值模拟方法越来越受到国内外研究者的重视) 有 旋转 以后得到以裂缝方位角为参数的应力应变关 限元方法将连续问题转变为离散问题，有限区域内 系方程