基于GPU并行加速的叠前逆时偏移方法.pdfVIP

东 北 石 油 大 学 学 报 第 36卷 第 4期 2012年 8月 JOURNALOFNORTHEASTPETROLEUM UNIVERSITY Vo1．36 No．4 Aug． 2012 基于 GPU并行加速的叠前逆时偏移方法 石 颖 ，陆加敏 ，柯 璇 ，田东升 ，王 菲。 (1．东北石油大学 地球科学学院 ，黑龙江 大庆 163318； 2．大庆 油田有 限责任公 司 勘探开发研究院 ，黑龙江 大庆 163712； 3．中国石油大学 (北京)地球科 学学 院，北 京 102200) 摘 要 ：为 了提高复杂地下介质 的成像精度和偏移算法的计算效率 ，提 出可高效对地下复杂构造进行准确成像 的 GPU加速叠前逆时偏移方法．该方法采用双程声波方程进行波场延拓，突破倾角限制，借助于高阶有限差分方法实现叠 前逆时偏移成像 ；利用 GPU(GraphicProcessingUnit)并行加速技术对波场延拓和成像进行计算 ，相 比于传统算法 ，其计 算效率有较大提高，可 以解决叠前逆时偏移算法计算量过大 问题 ；在获取波场信息过程中，也采用随机边界条件 ，实施 以 计算换存储策略，解决逆时偏移计算中的海量存储问题．模型测试结果表 明，该方法能够高效和高精度地对地下复杂地 质体成像． 关 键 词：逆时偏移；GPU；加速 ；高阶有限差分；随机边界条件；复杂构造 中图分类号 ：TE132．1 文献标识码 ：A 文章编号 ：2095—4107(2012)04—0l11—05 0 引言 随着勘探 目标的 日趋复杂化 ，叠前时间偏移 无法完全满足地震勘探对成像技术 的需求．在叠前深度 偏移方法中，叠前逆时偏移作为一种高精度地震成像方法 ，可以对回转波 、棱柱波和多次波等成像 ，突破倾 角的限制 ，能够对复杂地质构造进行精确成像．逆时偏移应用于地震勘探领域始于 20世纪 8O年代 。]， 受到当时硬件计算和存储条件约束影响，并未得到充分发展．近年来 ，随着低成本 的并行计算设备和有效 存储硬件的出现 ，逆时偏移在地震成像领域获得新的发展，逐渐成为主要的成像复杂构造 的叠前深度偏移 方法．然而 ，计算成本高、存储量大以及低频噪音 问题是 困扰逆时偏移方法发展 的瓶颈．刘红伟Ⅲ和李 博I8等通过有 限差分算法及 GPU 实现逆时偏移 ，提高逆时偏移算法 的计算效率．RobertGCE 提出随 机边界的思想 ，使得边界反射波以随机噪音的形式成像 ，进而可以通过波场反向传播计算模拟各个时刻的 波场 ，该方法在很大程度上节约了存储量．ZhangY_】和刘 红伟 。利用拉普拉斯算子滤波方法消除低频 成像噪音． 笔者利用 GPU／CPU协同并行加速实现逆时偏移算法中的波场延拓和相关成像计算 ，与传统方法相 比，能够大幅度提高逆时偏移的计算效率 ，缩短计算时间，节省地震偏移成像 的处理周期 ，可以解决叠前逆 时偏移计算量 巨大的问题 ；利用随机边界解决存储 问题 ，利用拉普拉斯算子滤波方法去除低频噪音 ，使得 叠前逆时偏移技术在工业领域得 以广泛应用成为可能． 1 逆时偏移原理 叠前逆时偏移算法的实现步骤包括震源波场沿时间方 向的正向延拓、检波点波场沿时间方 向的反向 延拓和正确应用成像条件．在反向延拓中，需要计算从最大时刻到零时刻的所有时间点上地下介质的波 场．逆时偏移的成像条件通常包括激发时刻成像条件、互相关成像条件及振幅比值成像条件．较为常用的 互相关成像条件需要使用在同一时刻的震源波场和检波点波场．由于前者是正传波场 ，后者是反传波场 ； 若想同时得到相同时刻的2个波场 ，则必须存储其中一个波场的整个传播过程的波场信息，即每一时刻的 收稿 Et期 ：2012—05—18；编辑 ：任志平 基金项 目：国家 “863”高技术研究发展计划项 目(2O12AAO612o2)；国家 自然科学基金青年基金项 目41004057)；黑龙江省教 育厅科学技术研究项 目；黑龙江省博士后科学基金项 目(LBH—Z11272) 作者简介 ：石 颖(1976一)，女 ，博士 ，副教授 ，主要从事地震资料处理方面的研究． 东 北 石 油 大 学