基于精细积分法的高性能三维地震正演模拟.pdfVIP

219 SPG/SEG 深圳 2011 国际地球物理会议论文集 基于精细积分法的高性能三维地震正演模拟 * 段玉婷 胡天跃 李靖宇 (北京大学地球与空间科学学院,地球物理系) 摘要 精细积分方法在空间域采用有限差分，时间域采用精细积分，提高了波动方程求解的精度，同 时保持了高效的特点，其稳定性条件也较有限差分法宽松。更适合复杂地质体和起伏地表的地震正演 模拟。为提高程序效率，本文引入了基于图形处理器（GPU）的优化技术，实现了三维高性能高精度的 地震正演模拟。 关键词 高性能计算 精细积分法 正演模拟 分法在在时间域得到的解是精确的，因此解决 1 前言 了有限差分法频散严重的问题，从而提高了整 个计算格式的精度。同时空间域采用单点子域 中国西北部的复杂地质情况给油藏预测 可以克服有限差分带宽较大的弱点，达到计算 和地球物理方法的实现带来了很多的困难,包 精度和效率的折衷。该方法数据量小，稳定性 括近地表的低速带,深层的多孔介质和油水混 好，同时又能较好地适应非均匀介质和复杂边 合成藏的情况。二维的任意精细积分法(ADPI) 界条件，对复杂的地质体模型和起伏地表均有 在复杂结构地区的成功应用,已经证明了其对 很好的效果。 计算精确度的提高。王润秋(2004) 将精细积分 法成功应用到地球物理正演模拟中来.唐跟阳 考虑方程组 (2007)利用改进型精细积分法，对二阶声波方    MU CU KU F  程进行求解。但是,由于三维模型的计算产生了 U(0) U0 (1)   巨大的计算量,提高计算效率成为了一个必须 U(0) U0  解决的问题。 M K 其中 是质量矩阵，C 是阻尼矩阵， 是  F 本文中将三维精细积分法正演模拟在 GPU 刚度矩阵。 是等效外力载荷，U0 和U0 都是 的平台上实现,达到了高性能计算的目的,通过 关于矩阵U 的初始条件。