论文部分内容阅读
在过去的几十年里,基于单程波波动方程的传统单程波算子广泛应用在理论研究和实际生产中。与双程波算子相比,单程波算子效率更高,需要的计算资源更少。在解决大型三维问题时,单程波算子的这两个优势显得尤其重要。然而由于与生俱来的广角限制,传统单程波算子不能模拟传播角度接近或超过90°的地震波。由于单程波算子只能沿着特定的方向传播(一般是深度方向),所以模拟回转波和反射波的问题超出了单程波算子的能力范围。我们发展了三维超广角单程波传播算子来克服传统单程波算子的这种广角限制,并且在带有大角度结构的复杂模型中使用超广角算子来模拟传播角度为超广角(即大于90°)的回转波,以此解决模拟和成像问题。具体实现中,通过特定的叠加格式,将下行波场和水平传播的波场叠加起来,保留优势方向上的正确波场,使得波场计算更加准确。为了将不同方向的波场进行叠加,设计了依赖于传播角度的权函数,并讨论了权函数的性质。由于三维问题计算的复杂性,在GPU上进行并行计算,以提高效率。为此,设计了一套适用于GPU架构的工作流程。数值算例表明,超广角单程波算子比传统单程波算子具有更高的计算精度,在含有竖直结构的成像中有更好的表现。事实上,超广角单程波传播算子可以基于任何单程波方法,提高地震波模拟和成像的效果。虽然超广角算子的计算效率和内存消耗比单程波算子差,但是优于双程波方法,实际中也是可以接受的。