常犙衰减介质分数阶波动方程优化有限差分模拟-地震学报.PDF

第 卷 第 期 地 震 学 报 ， ３９ ３ Ｖｏｌ．３９ Ｎｏ．３ 　 　 　 　 　 　 年 月 （ ） ＡＣＴＡＳＥＩＳＭＯＬＯＧＩＣＡＳＩＮＩＣＡ ， ２０１７ ５ ３４３３５５ Ｍａ ２０１７ 　 　 　 ｙ 孙成禹，乔志浩，伍敦仕，滕腾 常 衰减介质分数阶波动方程优化有限差分模拟 地震学报， （）： ．２０１７． 犙 ． ３９３ ３４３３５５． ： ／ ｄｏｉ１０．１１９３９ ａｓｓ．２０１７．０３．００４． ｊ ， ， ， ＳｕｎＣＹ ＱｉａｏＺＨ ＷｕＤＳ Ｔｅｎ Ｔ．２０１７．Ｍｏｄｅｌｉｎ ｏｆｗａｖｅｅｕａｔｉｏｎｗｉｔｈｆｒａｃｔｉｏｎａｌｄｅｒｉｖａｔｉｖｅｕｓｉｎ ｏｔｉｍａｌｆｉｎｉｔｅ ｇ ｇ ｑ ｇ ｐ ｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄｉｎｃｏｎｓｔａｎｔ ， （）： ： ／ 犙 ａｔｔｅｎｕａｔｉｏｎｍｅｄｉａ．犃犮狋犪犛犲犻狊犿狅犾狅犻犮犪犛犻狀犻犮犪 ３９ ３ ３４３３５５．ｄｏｉ１０．１１９３９ ａｓｓ． 犵 ｊ ２０１７．０３．００４． 常 衰减介质分数阶波动方程 犙 优化有限差分模拟 １） １）， １） ２） 孙成禹 乔志浩 伍敦仕 滕 腾 　 　 　 　 １）中国山东青岛２６６５８０中国石油大学（华东）地球科学与技术学院 ２）中国天津 ３００４５０中海油能源发展股份有限公司工程技术分公司 摘要 本文基于 常 模型理论，推导了常 衰减介质中黏声波和黏弹性波的速 Ｋａｒｔａｎｓｓｏｎ 犙 犙 　 ｊ 度 应力方程，并采用基于二项式窗函数的优化交错网格有限差分方法进行了数值模拟，同时  引入不分裂的复频移卷积完全匹配层（ ）吸收边界条件，以消除边界反射 使用基于自 ＣＰＭＬ ． 适应时间步长记忆方法的中心差分近似时间分数阶导数，与常用的短时记忆方法相比，提高 了波动方程的离散化精度和计算效率．通过对比均匀模型下声波的数值解与解析解，验证了 算法的精确性，并进一步分析了不同