[发明专利]地震纵波吸收边界条件处理方法及装置

权利要求书

1.一种地震纵波吸收边界条件处理方法，其特征在于，所述方法包括：

分别确定地震纵波传播的三维速度模型在三个维度的吸收边界条件，三维速度模型的三个维度分别表示GPU内存的最慢速方向、第二快速方向和最快速方向；

初始化三维速度模型中每一点在三维地震波数值模拟的初始时刻以及初始时刻的前一时刻的波场值；

针对三维数值模拟的初始时刻至终止时刻之间的每个时刻，执行下述方法计算吸收边界条件内每个点的波场值：

根据当前时刻和当前时刻的前一时刻的波场值，利用声波动方程的数值算法计算当前时刻的下一时刻的第一波场值；

根据第一波场值，利用GPU的寄存器阻塞算法计算当前时刻的下一时刻的单程波波场值，作为第二波场值；利用第一波场值和第二波场值计算GPU内存的最慢速方向和第二快速方向的吸收边界条件内每个点的波场值；

根据第一波场值，利用GPU的共享内存阻塞算法计算当前时刻的下一时刻的单程波波场值，作为第三波场值；利用第一波场值和第三波场值计算GPU内存的最快速方向的吸收边界条件内每个点的波场值。

2.根据权利要求1所述的方法，其特征在于，在初始化三维速度模型中的每一点在三维地震波数值模拟的初始时刻及初始时刻的前一时刻的波场值之后，所述方法包括：

将初始时刻作为当前时刻，计算当前时刻的下一时刻三维速度模型的吸收边界条件内每个点的波场值；

在初始时刻的下一时刻三维速度模型的吸收边界条件内每个点的波场值计算完毕后，将初始时刻的下一时刻作为当前时刻，重新计算当前时刻的下一时刻三维速度模型的吸收边界条件内每个点的波场值。

3.根据权利要求1或2所述的方法，其特征在于，所述利用声波动方程的数值算法计算当前时刻的下一时刻的第一波场值，包括：

根据计算当前时刻的下一时刻的第一波场值

其中，＝为赋值号，y、x和z分别表示当前计算点在空间坐标系中的位置，y表示GPU内存的最慢速方向，x表示GPU内存的第二快速方向，z表示GPU内存的最快速方向，y＝0,1,2,···,ny-1，x＝0,1,2,···,nx-1，z＝0,1,2,···,nz-1，ny、nx和nz分别用于表示三维速度模型在三个维度上的大小；n表示当前时刻，n-1表示当前时刻的上一时刻，n+1表示当前时刻的下一时刻；Δs表示两个相邻计算点在空间上的间隔；Δt表示两个相邻计算点在时间上的间隔；c表示地震纵波在介质中的传播速度；f表示震源函数。

4.根据权利要求3所述的方法，其特征在于，所述利用GPU的寄存器阻塞算法计算当前时刻的下一时刻的单程波波场值，作为第二波场值，包括：

根据当前时刻与当前时刻的前一时刻吸收边界条件中每个点的第一波场值，定义寄存器的初始值；

根据寄存器的初始值，确定当前时刻的下一时刻的单程波波场值。

5.根据权利要求4所述的方法，其特征在于，定义的寄存器为和针对GPU内存的最慢速方向，所述寄存器的初始值为和针对GPU内存的第二快速方向，所述寄存器的初始值为和其中，y₁∈y，y₁＝w-1,w-2,···,0；x₁∈x，x₁＝w-1,w-2,···,0，w用于表示吸收边界条件的长度。

6.根据权利要求5所述的方法，其特征在于，针对GPU内存的最慢速方向，所述根据寄存器的初始值，确定当前时刻的下一时刻的单程波波场值，包括：

根据计算当前时刻的下一时刻的单程波波场值

利用第一波场值和第二波场值计算GPU内存的最慢速方向的吸收边界条件内每个点的波场值，包括：

根据计算GPU内存的最慢速方向的吸收边界条件内点(y₁,x,z)的波场值；

其中，T₁₁＝(2-r)(1-r)/2，T₁₂＝r(2-r)，T₁₃＝r(r-1)/2，r＝cΔt/Δs。