基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移.PDF

第 卷 第 期 地 球 物 理 学 报 ， ６０ ４ Ｖｏｌ．６０ Ｎｏ．４ 　 　 　 　 　 年 月 ， ２０１７ ４ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ａ ｒ．２０１７ 　 　 　 ｐ 吴玉，符力耘，陈高祥 基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移 地球物理学报， （）： ．２０１７． ． ６０４ ， ： ／ １５２７１５３７ｄｏｉ１０．６０３８ ｊｇ ， ， ＷｕＹ ＦｕＬＹ ＣｈｅｎＧＸ．２０１７．Ｆｏｒｗａｒｄｍｏｄｅｌｉｎ ａｎｄｒｅｖｅｒｓｅｔｉｍｅｍｉｒａｔｉｏｎｏｆｖｉｓｃｏａｃｏｕｓｔｉｃｍｅｄｉａｕｓｉｎ ｄｅｃｏｕｌｅｄ ｇ ｇ ｇ ｐ （ ）， （）： ， ： ／ ｆｒａｃｔｉｏｎａｌＬａｌａｃｉａｎｓ．犆犺犻狀犲狊犲犑．犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ６０４ １５２７１５３７ｄｏｉ１０．６０３８ ｐ 狆狔 ｊｇ 基于分数阶拉普拉斯算子解耦的黏声介质地震 正演模拟与逆时偏移 ， ， １２ １ １２ 吴玉 ，符力耘 ，陈高祥 １中国科学院地质与地球物理研究所，中国科学院油气资源研究重点实验室，北京　１０００２９ ２中国科学院大学，北京 １０００４９ 　 摘要 时间域常 黏声波方程，由于含分数阶时间导数项，数值求解需要大量内存，计算效率低，不利于地震偏移 犙 　 的实施 通过一系列近似，可将该方程简化为介质频散效应和衰减效应解耦的分数阶拉普拉斯算子黏声波方程，数 ． 值求解内存需求少，计算效率高 本文采用交错网格有限差分逼近时间导数，改进的伪谱法计算空间导数， 吸 ． ＰＭＬ 收边界去除边界反射，对该方程进行数值离散和地震正演模拟 开展地震数据的黏声介质逆时偏移，实现波场逆时 ． 延拓过程中同时完成频散校正和衰减补偿 改善深层构造的成像精度，数值结果表明，基于分数阶拉普拉斯算子解 ． 耦的黏声介质地震正演模拟与逆时偏移可大幅度提高地震模拟计算效率，偏移剖面明显优于常规声波偏移剖面， 极大改善深层构造的成像品质． 关键词 时间域常 黏声波方程；分数阶拉普拉斯算子；频散与衰减解耦；黏声介质地震模拟与逆时偏移