全部 标题 作者
关键词 摘要

OALib Journal期刊
ISSN: 2333-9721
费用:99美元

查看量下载量

相关文章

更多...

大地电磁三维数据空间反演并行算法研究

DOI: 10.6038/j.issn.0001-5733.2012.12.009, PP. 3969-3978

Keywords: 大地电磁,三维反演,数据空间,并行计算,MPI

Full-Text   Cite this paper   Add to My Lib

Abstract:

目前大地电磁三维反演实际应用的主要问题是计算效率低.在对大地电磁三维数据空间反演算法进行深入分析的基础上,本文提出了基于频点和矩阵划分的大粒度并行反演方案和具体实现步骤,并在曙光TC5000A高性能计算平台上实现了基于MPI的大地电磁三维数据空间反演并行算法.该算法实现了包括三维正演、灵敏度矩阵、叉积矩阵以及模型改正量的并行执行,不仅计算效率高,而且每个节点机上灵敏度矩阵的存储空间只需原来微机上的2/N(N是参加并行计算的节点机个数),大大地减少了内存开销.通过两个理论模型合成的数据对实现的三维数据空间反演并行算法进行试算,对比分析了多个节点机下程序的执行效率.测试结果表明,所实现的三维数据空间反演并行算法是可行的、高效的,与单机相比,不仅可以提高运行速度,缩短计算时间,而且还可以扩大计算规模,极大地推动了大地电磁三维反演的实用化.

References

[1]  林昌洪, 谭捍东, 佟拓. 利用大地电磁三维反演方法获得二维剖面附近三维电阻率结构的可行性. 地球物理学报, 2011, 54(1): 245-256. Lin C H, Tan H D, Tong T. The possibility of obtaining nearby 3D resistivity structure from magnetotelluric 2D profile data using 3D inversion. Chinese J. Geophys. (in Chinese), 2011, 54(1): 245-256.
[2]  Mackie R L, Madden T R. Three-dimensional magnetotelluric inversion using conjugate gradients. Geophys. J. Int., 1993, 115(1): 215-229.
[3]  Spichak V, Popova I. Artificial neural network inversion of magnetotelluric data in terms of three-dimensional earth macroparameters. Geophys. J. Int., 2000, 142(1): 15-26.
[4]  Yuan C Z, Paulson K V. Magnetotelluric inversion using regularized Hopfield neural networks. Geophysical Prospecting, 1997, 45(5): 725-743.
[5]  New man G A, Alumbaugh D L. Three-dimensional massively parallel electromagnetic inversion-Ⅰ. Theory. Geophys. J. Int, 1997, 128(2): 345-354.
[6]  Zyserman F I, Santos J E. Parallel finite element algorithm with domain decomposition for three-dimensional magnetotelluric modelling. J. Appl. Geophys., 2000, 44(4): 337-351.
[7]  陈金窗, 戴光明. 微机网络并行计算及2.5维CSAMT正演的并行实现. 物探化探计算技术, 1997, 19(2): 103-107. Chen J C, Dai G M. Micro-computer networked computing and 2.5-D csamt forward parallel computing. Computing techniques for Geophysical and Geochemical Exploration (in Chinese), 1997, 19(2): 103-107.
[8]  Tan H D, Tong T, Lin C H. The parallel 3D magnetotelluric forward modeling algorithm. Applied Geophysics, 2006, 3(4): 197-202.
[9]  刘羽, 王家映, 孟永良. 基于PC机群的大地电磁Occam反演并行计算研究. 石油物探, 2006, 45(3): 311-315. Liu Y, Wang J Y, Meng Y L. PC cluster based magnetotelluric 2-D Occam''s inversion parallel calculation. Geophysical Prospecting for Petroleum (in Chinese), 2006, 45(3): 311-315.
[10]  李焱、胡祥云、金钢燮等. 基于MPI 的一维大地电磁并行计算研究. 地球物理学进展, 2010, 25(5): 1612-1616. Li Y, Hu X Y, Kim K, et al. Research of 1-D magnetotelluric parallel computation based on MPI. Progress in Geophys. (in Chinese), 2010, 25(5): 1612-1616.
[11]  Boonchaisuk S, Vachiratienchai Ci, Siripunvaraporn W. Two-dimensional direct current (DC) resistivity inversion: Data space Occam''s approach. Physics of the Earth and Planetary Interiors, 2008, 168(3-4): 204-211.
[12]  张林波, 迟学斌, 莫则尧等. 并行计算导论. 北京: 清华大学出版社, 2006. Zhang L B, Chi X B, Mo Z Y, et al. Introduction to Parallel Computing (in Chinese). Beijing: Tsinghua University Press, 2006.
[13]  都志辉. 高性能计算并行编程技术: MPI并行程序设计. 北京: 清华大学出版社, 2001. Du Z H. High-powered computing parallel programming technique-MPI parallel program design (in Chinese). Beijing: Tsinghua University Press, 2001.
[14]  Siripunvaraporn W, Gary E. An efficient data-subspace inversion method for 2-D magnetotelluric data. Geophysics, 2000, 65(3): 791-803.
[15]  魏文博. 我国大地电磁测深新进展及瞻望. 地球物理学进展, 2002, 17(2): 245-254. Wei W B. New advance and prospect of Magnetotelluric Sounding (MT) in China. Progress in Geophysics (in Chinese), 2002, 17(2): 245-254.
[16]  谭捍东, 余钦范, Booker J等. 大地电磁法三维快速松弛反演. 地球物理学报, 2003, 46(6): 850-855. Tan H D, Yu Q F, Booker J, et al. Three-dimensional rapid relaxation inversion for the magnetotelluric method. Chinese J. Geophys. (in Chinese), 2003, 46(6): 850-855.
[17]  Newman G A, Recher S, Tezkan B, et al. 3D inversion of a scalar radio magnetotelluric field data set. Geophysics, 2003, 68(3): 791-802.
[18]  Sasaki Y. Three-dimensional inversion of static-shifted magnetotelluric data. Earth, Planets and Space, 2004, 56(2): 239-248.
[19]  Zhdanov M S, Tolstaya E. Minimum support nonlinear parametrization in the solution of a 3D magnetotelluric inverse problem. Inverse Problems, 2004, 20(3): 937-952.
[20]  Siripunvaraporn W, Uyeshima M, Egbert G. Three-dimensional inversion for Network- Magnetotelluric data. Earth, Planets and Space, 2004, 56(9): 893-902.
[21]  Siripunvaraporn W, Egbert G, Lenbury Y, et al. Three-dimensional magnetotelluric inversion: data-space method. Physics of the Earth and Planetary Interiors, 2005, 150(1-3): 3-14.
[22]  胡祖志, 胡祥云. 大地电磁三维反演方法综述. 地球物理学进展, 2005, 20(1): 214-220. Hu Z Z, Hu X Y. Review of three dimensional magnetotelluric inversion methods. Progress in Geophysics (in Chinese), 2005, 20(1): 214-220.
[23]  胡祖志, 胡祥云, 何展翔. 三维大地电磁数据的二维反演解释. 石油地球物理勘探, 2005, 40(3): 353-359. Hu Z Z, Hu X Y, He Z X. Using 2-D inversion for interpretation of 3-D MT data. Oil Geophysical Prospecting (in Chinese), 2005, 40(3): 353-359.
[24]  胡祖志, 胡祥云, 何展翔. 大地电磁非线性共轭梯度拟三维反演. 地球物理学报, 2006, 49(4): 1226-1234. Hu Z Z, Hu X Y, He Z Z. Pseudo-Three-Dimensional magnetotelluric inversion using nonlinear conjugate gradients. Chinese J. Geophys. (in Chinese), 2006, 49(4): 1226-1234.
[25]  Avdeev D B, Avdeeva A D. A rigorous three-dimensional magnetotelluric inversion. Progress in Electromagnetics Research, 2006, 62: 41-48.
[26]  杨迪琨, 胡祥云. 含噪声数据反演的概率描述. 地球物理学报, 2008, 51(3): 901-907. Yang D K, Hu X Y. Inversion of noisy data by probabilistic methodology. Chinese J. Geophys. (in Chinese), 2008, 51(3): 901-907.
[27]  Han N, Nam M J, Kim H J, et al. Efficient three-dimensional inversion of magnetotelluric data using approximate sensitivities. Geophysical Journal International, 2008, 175(2): 477-485.
[28]  Lin C H, Tan H D, Tong T. Three-dimensional conjugate gradient inversion of magnetotelluric sounding data. Applied Geophysics, 2008, 5(4): 314-321.
[29]  Avdeev D, Avdeeva A. 3D magnetotelluric inversion using a limited-memory quasi-Newton optimization. Geophysics, 2009, 74(3): 45-57.
[30]  Siripunvaraporn W, Egbert G. WSINV3DMT: Vertical magnetic field transfer function inversion and parallel implementation. Physics of the Earth and Planetary Interiors, 2009, 173(3-4): 317-329.
[31]  徐凯军, 李桐林, 刘展. 基于积分方程法的大地电磁三维反演. 物探化探计算技术, 2009, 31(6): 564-567. Xu K J, Li T L, Liu Z. Three dimensional magnetotelluric inversion based on integral equation method. Computing Techniques for Geophysical and Geochemical Exploration (in Chinese), 2009, 31(6): 564-567.
[32]  Lin C H, Tan H D, Tong T. Parallel rapid relaxation inversion of 3D magnetotelluric data. Applied Geophysics, 2009, 6(1): 77-83.
[33]  Maris V, Wannamaker P E. Parallelizing a 3D finite difference MT inversion algorithm on a multicore PC using OpenMP. Computers & Geosciences, 2010, 36(10): 1384-1387.
[34]  Siripunvaraporn W. Three-dimensional magnetotelluric inversion: An introductory guide for developers and users. Surveys in Geophysics, 2012, 33(1): 5-27.
[35]  Newman G A, Alumbaugh D L. Three-dimensional magnetotelluric inversion using non-linear conjugate gradients. Geophys. J. Int., 2000, 140(2): 410-424.
[36]  Zhdanov M S, Fang S, Hursán G. Electromagnetic inversion using quasi-linear approximation. Geophysics, 2000, 65(5): 1501-1513.
[37]  刘羽. 基于机群的二维大地电磁Occam反演的并行计算研究. 武汉: 中国地质大学, 2006. Liu Y. PC-Cluster Based Parallel Computation Research on 2-D Magnetotelluric Occam Inversion (in Chinese). Wuhan: China University of Geosciences, 2006.
[38]  李焱、胡祥云、吴桂桔等. 基于MPI的二维大地电磁正演的并行计算. 地震地质, 2010, 32(3): 392-401. Li Y, Hu X Y, Wu G J, et al. Parallel computation of 2-D magnetotelluric forward modeling based on MPI. Seismology and Geology (in Chinese), 2010, 32(3): 392-401.
[39]  李焱. 基于MPI的大地电磁三维正反演并行算法研究. 武汉: 中国地质大学, 2011. Li Y. Parallel Computation Research of 3-D Magnetotelluric Forward Modeling and Inversion Based on MPI[Master''s thesis](in Chinese). Wuhan: China University of Geosciences, 2011.
[40]  de Groot-Hedlin C, Constable S. Occam''s inversion to generate smooth, two-dimensional models from magnetotelluric data. Geophysics, 1990, 55(12): 1613-1624.

Full-Text

Contact Us

service@oalib.com

QQ:3279437679

WhatsApp +8615387084133