一种基于透射边界的混合吸收边界方法技术

本发明专利技术属于地震勘探技术领域，具体涉及一种基于透射边界的混合吸收边界方法，1、设置震源、检波点、地震记录等信息；2、对内部计算区域(区域一)、粘声波场与透射波场的过渡区域(区域二)计算，获得波场值；3、对内部计算区域和透射边界计算，获得透射边界波场值；4、构建双程和透射函数的权重函数并计算最终波场结果。本方法相比于传统吸收边界，本方法仅需较少的吸收层即可获得良好的吸收效果；在边界处引入了自适应权重函数，可同时压制内边界和外边界处的边界反射，进一步提高了吸收效果；在边界处仅需使用当前时刻及上一时刻的波场，即可获得下一时刻吸收区域的波场值，有利于多线程同时计算，降低了时间复杂度。降低了时间复杂度。降低了时间复杂度。

【技术实现步骤摘要】
一种基于透射边界的混合吸收边界方法


[0001]本专利技术属于地震勘探
，具体涉及一种用于时间域粘滞介质数值模拟的混合吸收边界方法。

技术介绍

[0002]地震波数值模拟是利用计算机模拟地震波在地下介质中的传播情况，其有助于正确认识地下介质形态及性质。该技术在天然地震预测和矿产资源普查中受到广泛的应用，尤其是地震资料处理中最为至关重要的环节之一。由于地下介质是非弹性介质，考虑地层吸收衰减作用的粘滞波场模拟技术更符合真实情况。其中，求解时间域分数阶拉氏方程的波场模拟技术(Zhu 2014)，是目前被应用最广的技术之一。然而，受计算能力限制，计算机无法模拟地下介质的全部区域。因此，一定存在截断边界，在数值模拟中必须压制这种截断噪声，否则在地震资料处理中会产生错误的成像结果。因此，在粘滞介质的数值模拟中必须有效解决这一问题。
[0003]目前，对地震波数值模拟中边界反射的压制可分为三大类：一是通过引入一定数量的衰减层达到吸收边界反射的目的，例如指数衰减边界、sponge吸收边界条件、完全匹配层(PML)边界条件等。其中PML边界条件在理本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于透射边界的混合吸收边界方法，其特征在于：所述混合吸收边界方法包括以下步骤：步骤一、设置震源、检波点和模型大小地震记录信息；步骤二、对内部计算区域、粘声波场与透射波场的过渡区域计算粘声波场，获得双程波波场值；内部计算区域和过渡带双程波动方程的计算，用二阶常Q粘声方程模拟地震波在粘滞介质中传播，其表示为：式中，p1和分别表示双程波粘滞波场和拉普拉斯算子，c0为定义在参考频率ω0处的参考相速度；Q为品质因子；步骤三、对过渡区和最外层边界计算透射波方程，获得透射方程波场值；过渡区和最外层边界波场的计算，利用多次透射公式可以预测地震波在边界处的传播情况，t+Δt时刻的波场值可表示为：式中，(x，y，z)表示空间位，p2表示透射波场，t表示当前时刻，Δt为时间步长，N表示透射阶数，F
k
，E
k
分别为递归系数矩阵和波场矩阵：递归系数矩阵和波场矩阵：式中：t1＝(2
‑
X)(2
‑
X)/2t2＝X(2
‑
X)t3＝X(X
‑
1...

【专利技术属性】
技术研发人员：王宁，石颖，李松龄，王维红，宋利伟，王银凤，
申请(专利权)人：东北石油大学，
类型：发明
国别省市：