[发明专利]一种三维各向异性弹性波数值模拟方法及系统

说明书

技术领域

本发明涉及一种三维各向异性弹性波数值模拟方法及系统，属于地球物理勘探领域。

背景技术

地震波数值模拟以地震波在地下介质中的传播理论为基础，在勘探地震学与天然地震学中得到了广泛的应用。目前常规的三维声波方程、弹性波各向同性数值模拟已经广泛的应用于数值模拟、成像、以及反演等各个地球物理领域中。但是对于如何有效的实现大尺度采集方式(例如宽方位采集)和复杂各向异性介质(例如水平横向各向同性HTI或者是正交各向异性)的三维弹性波波场数值模拟研究仍然存在很大的挑战，在实际应用中未得到广泛的使用。此外常规的基于CPU的大规模三维声波、弹性波数值模拟，通常需要大量专用集群计算资源。无论从硬件成本还是计算能耗上来说其成本都十分昂贵。因此快速的提高计算性能、显著的降低计算成本对于实现三维弹性波各向异性数值模拟的实际应用具有重要意义。

近十年中，利用图形处理单元(GPU)进行计算密集型应用的加速实现已经得到了突飞猛进式的发展。图形处理单元(GPU)由于其具有高速的内存带宽，相较于CPU至少高出两个数量级的计算处理核心，更适合并行计算的单指令多数据(SIMD)计算模式，以及更低的能耗成本，正广泛的应用到计算科学的相关领域。对于勘探地球物理领域，对于使用图形处理单元(GPU) 的兴趣也在显著增强，越来越多的研究已经将GPU用于加速地震处理中的核心算法，例如地震数值模拟、地震成像、地震高精度反演等。

发明内容

经过研究发现，现有技术及其存在以下问题：

将时间域有限差分方法应用到GPU设备上实现复杂介质波传播算法上的研究很少，尤其对于如何处理存在大量多节点数据交换的大规模三维问题研究就更少(Heinrich et al，2014)。

目前存在一些GPU地震波传播模拟的实现方案，对于三维GPU集群的实现方案我们以(龙桂华等)的方案为例进行描述。该方案利用区域分解技术将单个GPU上不能计算的地质体沿Z轴方向进行粗粒度分解，采用消息传递接口(MPI)交换边界数据，从而运用MPI+CUDA的方式实现了大尺度三维地震波场的数值模拟。但是该方法存在很大的一个问题是，采用GPU集群计算的加速比相较于单个GPU与CPU来说显著下降，造成这一结果的原因是GPU 计算的大部分耗时消耗在GPU集群节点间的GPU到CPU以及CPU到GPU设备之间的数据拷贝上。

本发明所要解决的技术问题是，针对现有技术的不足，提供一种将GPU Direct技术与有限差分数值模拟方法结合得到的基于GPU Direct优化通信的三维各向异性弹性波数值模拟方法及系统，为大规模数值模拟提供有力保证。

本发明解决上述技术问题的技术方案如下：一种三维各向异性弹性波数值模拟方法，具体包括以下步骤：

步骤1：建立介质模型，对介质模型进行网格离散得到多个网格点；

步骤2：计算震源函数，根据震源函数计算每个网格点上的压力值；

步骤3：将三维各向异性弹性波方程转换为传播方程，将每个网格点上的压力值带入传播方程进行计算，得到每一刻的波场值；

步骤4：根据波场值确定每个网格点的计算区域，进行分区并对分区边界数据进行数据交换，完成弹性波数值模拟。

本发明的有益效果是：从应力应变方程出发，实现了利用图形处理单元 (GPU)加速复杂介质下的三维弹性波数值模拟，并针对三维问题引入GPU设备所面对的由于区域分解所造成的多节点间、节点内通信瓶颈问题进行了深入的研究和分析，提出了利用GPU Direct技术加速数据传输的实现方案，避免了大量的CPU到GPU，GPU到CPU的数据拷贝，实现了优化通信瓶颈的问题。通过GPU图形加速设备并引入GPU Direct通信优化策略，可以显著的提升整体计算性能，可以用更低的硬件成本和更少的时间实现的三维各向异性弹性数值模拟，为各种依赖于波动方程正演模拟的算法如逆时偏移，全波形反演的应用提供有力的保证。

在上述技术方案的基础上，本发明还可以做如下改进。

进一步，所述步骤4具体包括以下步骤：

步骤4.1：根据所有波场值确定每个网格点的计算区域，并采用吸收边界的方式确定计算区域边界，在计算区域内模拟地下介质中波的传播；

步骤4.2：对所有计算区域进行分区，对相邻分区的边界数据进行数据交换，得到模拟弹性波数据，完成弹性波数值模拟。

进一步，所述步骤4.2中边界数据的采用GPU-Direct技术进行数据交换。