Regularized inversion of 3D gravity data:a new GPU parallelized method based on CUDA
LI Wu-Yang1,2, ZHANG Jian1,2, LIN Wei3
1. College of Earth Science, University of Chinese Academy of Sciences, Beijing 100049;
2. Key Laboratory of Computational Geodynamics, Chinese Academy of Sciences, Beijing 100049;
3. Department of Geology and Geophysics, University of Utah, Salt Lake City 84112
We introduce a new 3D parallelgravity inversion method based on CUDA GPU in PGI Fortran,which could be used in PCs and Laptops with NVDIA graphic cards to accelerate iteration speed for Re-Weight Regularized Conjugate Gradient method up to hundreds times.Storage and threads optimization was made for visual OS, leading a more than 0.1 billion cell's number for inversion in PCs.The result of model study shows that this method is an efficient and believable parallel computing algorithm.
李午阳, 张健, 林巍. 基于CUDA的GPU并行优化重力三维反演[J]. 物探与化探, 2016, 40(1): 179-184.
LI Wu-Yang, ZHANG Jian, LIN Wei. Regularized inversion of 3D gravity data:a new GPU parallelized method based on CUDA. Geophysical and Geochemical Exploration, 2016, 40(1): 179-184.
[1] Lü Q,Qi G,Yan J.3D geologic model of Shizishan ore field constrained by gravity and magnetic interactive modeling:A case history[J].Geophysics,2012,78(1):B25-B35.[2] Wang G,Zhang S,Yan C,et al.Mineral potential targeting and resource assessment based on 3D geological modeling in Luanchuan region,China[J].Computers & Geosciences,2011,37(12):1976-1988.[3] 姚长利,郝天珧.重磁遗传算法三维反演中高速计算及有效存储方法技术[J].地球物理学报,2003,46(2):252-258.[4] 陈召曦,孟小红,刘国峰,等.基于GPU的任意三维复杂形体重磁异常快速计算[J].物探与化探,2012,36(1):117-121.[5] Chen Z,Meng X,Guo L,et al.GICUDA:A parallel program for 3D correlation imaging of large scale gravity and gravity gradiometry data on graphics processing units with CUDA[J].Computers & Geosciences,2012,46:119-128.[6] 林巍,张健,李家彪.南沙中业群礁地区中生代残留盆地GPU重力异常三维反演[J].热带海洋学报,2013,32(4):36-42.[7] Cuma M,Wilson G A,Zhdanov M S.Large-scale 3D inversion of potential field data[J].Geophysical Prospecting,2012,60(6):1186-1199.[8] Cuma M,Zhdanov M S.Massively parallel regularized 3D inversion of potential fields on CPUs and GPUs[J].Computers & Geosciences,2014,62:80-87.[9] Zhdanov M S.New advances in regularized inversion of gravity and electromagnetic data[J].Geophysical Prospecting,2009,57(4):463-478.[10] Portniaguine O,Zhdanov M S.Focusing geophysical inversion images[J].Geophysics,1999,64(3):874-887.[11] Zhdanov M S.Geophysical inverse theory and regularization problems[M].Elsevier,2002.[12] Cook S.CUDA programming:a developer's guide to parallel computing with GPUs[M].Newnes,2012.[13] Hoshino T,Maruyama N,Matsuoka S,et al.CUDA vs openACC:Performance case studies with kernel benchmarks and a memory-bound CFD application[C]//Cluster,Cloud and Grid Computing (CCGrid),201313th IEEE/ACM International Symposium on.IEEE,2013:136-143.[14] Zhdanov M S,Smith R B,Gribenko A,et al. Three-dimensional inversion of large-scale EarthScopemagnetotelluric data based on the integral equation method:Geoelectrical imaging of the Yellowstone conductive mantle plume[J].Geophysical Research Letters,2011,38(8):1-7.[15] Cox L H,Wilson G A,Zhdanov M S.3D inversion of airborne electromagnetic data using a moving footprint[J].Exploration Geophysics,2010,41(4):250-259.