书书书
第
５５
卷 第
４
期
２０１２
年
４
月
地
球
物
理
学
报
ＣＨＩＮＥＳＥ
ＪＯＵＲＮＡＬ
ＯＦ
ＧＥＯＰＨＹＳＩＣＳ
Ｖｏｌ．５５
，
Ｎｏ．４
Ａ
ｐ
ｒ．
，
２０１２
严红勇，刘洋
．
黏弹
ＴＴＩ
介质中旋转交错网格高阶有限差分数值模拟
．
地球物理学报，
２０１２
，
５５
（
４
）：
１３５４１３６５
，
ｄｏｉ
：
１０．６０３８
／
ｊ
．
ｉｓｓｎ．０００１５７３３．２０１２．０４．０３１．
Ｙａｎ Ｈ Ｙ
，
ＬｉｕＹ．Ｒｏｔａｔｅｄｓｔａ
ｇｇ
ｅｒｅｄ
ｇ
ｒｉｄｈｉ
ｇ
ｈｏｒｄｅｒｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅｎｕｍｅｒｉｃａｌｍｏｄｅｌｉｎ
ｇ
ｆｏｒｗａｖｅ
ｐ
ｒｏ
ｐ
ａ
ｇ
ａｔｉｏｎｉｎｖｉｓｃｏｅｌａｓｔｉｃ
ＴＴＩｍｅｄｉａ．
犆犺犻狀犲狊犲犑
．
犌犲狅
狆
犺
狔
狊
．
（
ｉｎＣｈｉｎｅｓｅ
），
２０１２
，
５５
（
４
）：
１３５４１３６５
，
ｄｏｉ
：
１０．６０３８
／
ｊ
．ｉｓｓｎ．０００１５７３３．２０１２．０４．０３１．
黏弹
犜犜犐
介质中旋转交错网格高阶
有限差分数值模拟
严红勇
１
，
２
，刘
洋
１
，
２
１
中国石油大学（北京）油气资源与探测国家重点实验室，北京
１０２２４９
２
中国石油大学（北京）
ＣＮＰＣ
物探重点实验室，北京
１０２２４９
摘
要
以
Ｃａｒｃｉｏｎｅ
黏弹各向异性理论为基础，给出了适用 于 黏 弹 性 具 有 任 意 倾 斜 对 称 轴 横 向 各 向 同 性 介 质（黏
弹
ＴＴＩ
介质）的二维三分量一阶速度
应力方程，采用旋转交 错 网 格 任 意 偶 数 阶 精 度 有 限 差 分 格 式 求 解 该 方 程，并
推导出了二维黏弹
ＴＴＩ
介质完全匹配层（
ＰＭＬ
）吸收边界条件公式和 相 应的 旋 转交 错 网格 任 意偶 数 阶精 度 有限 差
分格式，实现了该类介质的地震波场数值模拟
．
数值模拟结果表明：该方法模拟精 度 高，边 界 吸收 效 果好，可以 得到
高精度的波场快照和合成记录；并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征
．
关键词
黏弹性，各向异性，旋转交错网格，有限差分，完全匹配层
ｄｏｉ
：
１０．６０３８
／
ｊ
．ｉｓｓｎ．０００１５７３３．２０１２．０４．０３１
中图分类号
Ｐ６３１
收稿日期
２０１１０３１７
，
２０１２０３１６
收修定稿
基金项 目
国家自然科学基金项 目 （
４１０７４１００
）、国 家 高 技 术 研 究 发 展 计 划 （
８６３
）项 目 （
２００７ＡＡ０６Ｚ２１８
）和 教 育 部 新 世 纪 优 秀 人 才 支 持 计 划
（
ＮＣＥＴ１００８１２
）联合资助
．
作者简 介
严红勇，男，
１９８２
年生，中国石油大学（北京）博士研究生，主要从事 地 震波 传 播与 成 像、地 震 波吸 收 衰 减 及 地 震 资 料 高 分 辨 率 处 理
等方面的研究工作
．Ｅｍａｉｌ
：
ｙ
ａｎｈｏｎ
ｇｙ
ｏｎ
ｇ
＠
１６３．ｃｏｍ
通信作 者
刘洋，
Ｅｍａｉｌ
：
ｗｌｉｕ
ｙ
ａｎ
ｇ
＠
ｖｉ
ｐ
．ｓｉｎａ．ｃｏｍ
犚狅狋犪狋犲犱狊狋犪
犵犵
犲狉犲犱
犵
狉犻犱犺犻
犵
犺狅狉犱犲狉犳犻狀犻狋犲犱犻犳犳犲狉犲狀犮犲狀狌犿犲狉犻犮犪犾犿狅犱犲犾犻狀
犵
犳狅狉狑犪狏犲
狆
狉狅
狆
犪
犵
犪狋犻狅狀犻狀狏犻狊犮狅犲犾犪狊狋犻犮犜犜犐犿犲犱犻犪
ＹＡＮ Ｈｏｎ
ｇ
Ｙｏｎ
ｇ
１
，
２
，
ＬＩＵ Ｙａｎ
ｇ
１
，
２
１
犛狋犪狋犲犓犲
狔
犔犪犫狅狉犪狋狅狉
狔
狅
犳
犘犲狋狉狅犾犲狌犿 犚犲狊狅狌狉犮犲犪狀犱犘狉狅狊
狆
犲犮狋犻狀
犵
，
犆犺犻狀犪犝狀犻狏犲狉狊犻狋
狔
狅
犳
犘犲狋狉狅犾犲狌犿
，
犅犲犻
犼
犻狀
犵
１０２２４９
，
犆犺犻狀犪
２
犆犖犘犆犓犲
狔
犔犪犫狅狉犪狋狅狉
狔
狅
犳
犌犲狅
狆
犺
狔
狊犻犮犪犾犈狓
狆
犾狅狉犪狋犻狅狀
，
犆犺犻狀犪犝狀犻狏犲狉狊犻狋
狔
狅
犳
犘犲狋狉狅犾犲狌犿
，
犅犲犻
犼
犻狀
犵
１０２２４９
，
犆犺犻狀犪
犃犫狊狋狉犪犮狋
ＯｎｔｈｅｂａｓｉｓｏｆＣａｒｃｉｏｎｅ′ｓｔｈｅｏｒｉｅｓｏｆｖｉｓｃｏｅｌａｓｔｉｃｉｔ
ｙ
ａｎｄａｎｉｓｏｔｒｏ
ｐｙ
，
ｔｗｏｄｉｍｅｎｓｉｏｎａｌ
，
ｔｈｒｅｅｃｏｍ
ｐ
ｏｎｅｎｔ
，
ｆｉｒｓｔｏｒｄｅｒｖｅｌｏｃｉｔ
ｙ
ｓｔｒｅｓｓ ｗａｖｅｅ
ｑ
ｕａｔｉｏｎｓ ｏｆｖｉｓｃｏｅｌａｓｔｉｃｔｉｌｔｅｄｔｒａｎｓｖｅｒｓｅｌ
ｙ
ｉｓｏｔｒｏ
ｐ
ｉｃ
（
ｖｉｓｃｏｅｌａｓｔｉｃＴＴＩ
）
ｍｅｄｉａａｒｅ
ｐ
ｒｅｓｅｎｔｅｄａｎｄａｒｏｔａｔｅｄｓｔａ
ｇｇ
ｅｒｅｄ
ｇ
ｒｉｄａｎ
ｙ
ｏｒｄｅｒｆｉｎｉｔｅ
ｄｉｆｆｅｒｅｎｃｅｓｃｈｅｍｅｉｓｕｓｅｄｔｏｎｕｍｅｒｉｃａｌｌ
ｙ
ｓｏｌｖｅｔｈｅｅ
ｑ
ｕａｔｉｏｎｓ．Ｅ
ｑ
ｕａｔｉｏｎｓｏｆｔｈｅ
ｐ
ｅｒｆｅｃｔｌ
ｙ
ｍａｔｃｈｅｄ
ｌａ
ｙ
ｅｒ
（
ＰＭＬ
）
ａｒｅｄｅｒｉｖｅｄｆｏｒｔｈｅ ｗａｖｅｅ
ｑ
ｕａｔｉｏｎｓｉｎ ｖｉｓｃｏｅｌａｓｔｉｃ ＴＴＩ ｍｅｄｉａａｎｄｔｈｅｒｏｔａｔｅｄ
ｓｔａ
ｇｇ
ｅｒｅｄ
ｇ
ｒｉｄａｎ
ｙ
ｏｒｄｅｒｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅｓｃｈｅｍｅｉｓａｌｓｏｕｓｅｄｔｏｓｏｌｖｅｔｈｅｓｅｅ
ｑ
ｕａｔｉｏｎｓ．Ｒｅｓｕｌｔｓｏｆ
ｎｕｍｅｒｉｃａｌ ｍｏｄｅｌｉｎ
ｇ
ｉｎｄｉｃａｔｅｔｈａｔｔｈｅ ｍｏｄｅｌｉｎ
ｇ ｐ
ｒｅｃｉｓｉｏｎｉｓｈｉ
ｇ
ｈａｎｄｔｈｅａｂｓｏｒｂｉｎ
ｇ
ｂｏｕｎｄａｒ
ｙ
ｃｏｎｄｉｔｉｏｎｉｓ
ｇ
ｏｏｄｉｎｔｈｅｖｉｓｃｏｅｌａｓｔｉｃＴＴＩｍｅｄｉａ
，
ａｎｄｈｉ
ｇ
ｈ
ｐ
ｒｅｃｉｓｉｏｎｓｎａ
ｐ
ｓｈｏｔｓｏｆｗａｖｅｆｉｅｌｄａｎｄ
ｓ
ｙ
ｎｔｈｅｔｉｃｓｅｉｓｍｏ
ｇ
ｒａｍｓｃａｎｂｅｏｂｔａｉｎｅｄ
，
ａｎｄｔｈｅ
ｙ
ｃａｎｒｅｆｌｅｃｔｔｈｅｃｈａｒａｃｔｅｒｉｓｔｉｃｏｆｖｉｓｃｏｅｌａｓｔｉｃｉｔ
ｙ
ａｎｄａｎｉｓｏｔｒｏ
ｐｙ
．
- 1
- 2
前往页