逆时深度偏移成像方法及其在CPU%2fGPU异构平台上的实现.pdfVIP

岩性油气藏 2010年1勇 CEG会议专刊 UTHOLOGIC RESERVOIRS ·地震处理· 逆时深度偏移成像方法及其在CPU／GPU异构 平台上的实现 赵 磊，王华忠，刘守伟 (同济大学海洋与地球科学学院) 摘要：逆时深度偏移的优点可以用来对全波场逆时外推，使得来自高角度反射界面的反射渡，甚至来自 超过90。反射界面的反射波精确成像。逆时深度偏移涉及到2个重要问题：精确且高效的逆时波场外推算 法及反射波的成像条件。文中提出了修改的激励时间成像条件：在反射界面处。反射波出发时等于入射渡 到达时。认为用反射界面上某成像点与法线方向对称的波矢量所定义的波场进行相关成像是完善的成像 条件。边界条件的处理借鉴Robert提出的随机边界，这可以提高波场外推算法的并行性。关于存储，认为 牺牲一次正向波场外推，先把波场正传到最大时间，然后与观测波场同时逆时外推并成像。在GPU／CPU计 算机平台上是比较合适的方案，这样就彻底抛弃了正向外推波场的存储。数值试验验证了分析的正确性。 关键词：逆时偏移：高阶有限差分；GPU／CPU平台：随机边界条件 中图分类号：P631．4 文献标识码：A 0 引言 逆时偏移成像条件方面，Claerboutcz]提出的上、 下行波相关成像条件是被研究者普遍采用的成像 目前常用的叠前深度偏移方法主要有2种：单 条件。但是，逆时偏移成像是用双向波方程逆时间 程波波动方程偏移和克希霍夫积分偏移。单程波波 外推波场与沿正时间方向外推的震源波场在空间 动方程偏移基于双向波方程的单向波分解。这样的 任何一点进行零延迟的互相关。从传播路径看，这 分解只有在常速情况下才是精确成立的。利用差分 2个波场在从炮点到反射点．再到检波点的整个路 方法求解单向波方程，需要对单向波方程进行傍轴 径上是重合在一起的，就是说它们在空间上任何一 近似，不能对陡倾角反射界面准确成像。为此．对陡 点都正相关。这就是逆时偏移产生很强的低频噪音 倾角成像时通常采用克希霍夫积分偏移。但克希霍 夫积分偏移不方便处理单成像点多到达时的成像 完全适用于逆时偏移，因为震源波场和检波点波场 问题。另外，克希霍夫积分法只能描述波在光滑介 中既包含上行波又包含下行波。单向波偏移时．震 质中的传播，而且焦散现象难以克服。因此，单程波 波动方程偏移和克希霍夫积分偏移都无法对复杂 源波场中仅有下行波，而检波点波场中仅有上行波， 介质特别是盐下构造实现很好的成像，于是业界把 不会出现从炮点到反射点。再到检波点的整个路径 目光重新转向逆时偏移。 上2个波场重合在一起的现象。因此，Claerbout 提出的上、下行波相关成像条件在单向波偏移时不 逆时偏移最早由Whitemoretl]提出。相关研究 工作主要集中在2个方面：①波动方程的逆时外推 算法；②逆时偏移成像条件及其计算方法。 结果高通滤波可以移除假象：Yoon等[4]提出通过对 波动方程