利用CUDA的PML边界三维地震波传播模拟方法

摘要

本发明提供了一种利用CUDA的PML边界三维地震波传播模拟方法，方法包括：读入初始参数和雷克子波；将三维地震波的第一声波波动方程中的应力分解成相互垂直三个方向的应力分量，并加入PML边界条件形成第二声波波动方程；运用有限差分方法将第二声波波动方程离散化形成有限差分声波波动方程；给第一三维速度模型加边形成一第二三维速度模型；根据有限差分声波波动方程在处理器CPU中计算第二三维速度模型的棱、角、面的衰减系数；根据所述衰减系数在图像处理器GPU中计算所述第二三维速度模型网格点的应力；根据所述应力输出最后时刻的三维地震波模拟记录。本发明可以获得较高的加速比，缩短模拟时间。

主权项

一种利用计算统一设备架构CUDA的最佳匹配层PML边界三维地震波传播模拟方法，其特征在于，所述利用CUDA的PML边界三维地震波传播模拟方法包括：读入初始参数和雷克子波；将三维地震波的第一声波波动方程中的应力分解成相互垂直三个方向的应力分量，并加入PML边界条件形成第二声波波动方程；运用有限差分方法将所述第二声波波动方程离散化形成有限差分声波波动方程；给所述第一三维速度模型加边形成一第二三维速度模型；根据所述有限差分声波波动方程在中央处理器CPU中计算所述第二三维速度模型的棱、角、面的衰减系数；根据所述衰减系数在图像处理器GPU中计算所述第二三维速度模型网格点的应力；以及根据所述应力输出最后时刻的三维地震波模拟记录。