基于CUDA的地震数据层位面自动追踪算法

英文篇名：CUDA-based Automatic Seismic Horizon Tracking Algorithm
中文刊名：计算机技术与发展
英文刊名：Computer Technology and Development
作者：魏强 ; 唐杰
英文作者：WEI Qiang ; TANG Jie(Department of Computer Science and Technology ; Nanjing University ; Nanjing 210093 ; China)
中文关键词：CUDA ; 层位面 ; 地震数据处理 ; 自动追踪算法
英文关键词：CUDA ; horizon ; seismic data processing ; automatic tracking algorithm
出版日期：2012-09-10
机构：南京大学计算机科学与技术系;
年：2012
期：09
出版单位：计算机技术与发展

摘要

层位面追踪是三维地震数据解释的一项关键操作,对于地质勘探和油藏工程有重要的指导意义。传统的层位面自动追踪算法采用的是种子点蔓延法,该方法是一个不断进行迭代计算的过程。由于三维地震数据的数据量很大,传统的基于CPU的层位面自动追踪算法的计算效率不高。文中提出了一种基于CUDA的三维层位面自动追踪算法。该算法对种子点蔓延法进行了改进,实现了层位面自动追踪的GPU并行解决方案。实验表明,与传统的基于CPU的层位面自动追踪算法相比,本算法的计算效率有明显的提升。
Horizon tracking is an important operation in 3-D seismic data processing,and the tracked horizon can be used by the experts to detect the earth's subsurface structure.Processing speed is a major bottleneck in seismic horizon tracking since the data size of 3-D seismic data is usually very huge.In this paper,present a new automatic horizon tracking algorithm,which is based on CUDA technique.The algorithm fully utilizes the many-core structure and the high performance computing ability of GPU,realizing a large scale parallel method for horizon tracking.The results of experiments show that the algorithm has an obvious speedup than the traditional CPU-based horizon tracking algorithm.

引文

[1]陆基孟.地震勘探原理(上册)[M].东营:石油大学出版社,1993.
    [2]Faraklioti M,Petrou M.Horizon picking in 3D seismic datavolumes[J].Machine Vision and Applications,2004,15(4):216-219.
    [3]李雪峰,阎建国,赵州,等.利用相干属性剖面特征进行层位解释[J].物探化计算技术,2011,33(2):134-139.
    [4]Bahorich M,Farmer S.The coherence cube[J].The LeadingEdge,1995,14(10):1053-1058.
    [5]NVIDIA Corporation.NVIDIA CUDA C Programming Guide,Version 4.0[EB/OL].2011.http://developer.nvidia.com/nvidia-gpu-computing-documentation.
    [6]Sanders J,Kandrot E.CUDA by Example:An Introduction toGeneral-purpose GPU Programming[M].北京:清华大学出版社,2010.
    [7]Deschizeaux B,Blanc Jean-Yves.Imaging Earth's SubsurfaceUsing CUDA[C]//GPU Gems 3.[s.l.]:Addison WesleyProfessional,2007.
    [8]Keskes N,Zaccagnino P,Rether D,et al.Automatic extractionof 3-D seismic horizons[C]//Annual Meeting Expanded Ab-stracts.[s.l.]:SEG,1983:557-559.
    [9]dGB Earth Sciences.OpendTect User Documentation Version4.0[EB/OL].2009.http://opendtect.org/rel/doc/User/base/.
    [10]Harris M.Parallel Prefix Sum with CUDA[C]//GPU Gems3.[s.l.]:Addison Wesley Professional,2007.
    [11]Harish P,Narayanan P J.Accelerating large graph algorithmson the GPU using CUDA[C]//High performance computing-HIPC 2007,lecture notes in computer science.[s.l.]:[s.n.],2007:197-208.
    [12]Kirk D B,Hwu W W.Programming Massively Parallel Proces-sors[M].[s.l.]:Morgan Kaufmann,2010:96-103.