摘要
复电阻率法二维正演采用二次场算法,结合有限单元法进行离散化,然后引入ColeCole模型,完成数值模拟。通过理论模型试算验证了二维正演程序的正确性。当模型剖分较细时,多发射源、多频率的正演计算效率有待提高。复电阻率法二维数值模拟是按发射源或频率计算,且各发射源之间、各频率之间计算是相互独立且互不影响的,因此通过引入MPI并行计算来提高其计算效率,从而减少其正演所需的时间。采用串行、并行复电阻率法二维正演程序来计算同一理论模型,结果验证了该并行程序的正确性。对并行算法的效率进行分析表明,该并行算法是可靠的、高效的。
The complex resistivity method is an artificial source exploration method which is based on the electrical difference of rocks.Two-dimensional forward modeling employs secondary field algorithm,by using the finite element method to discretize the differential equation,and then introduce Cole-Cole model,the numerical simulation is completed.The calculating synthetic model tests,indicate that the accuracy of the forward algorithm is reliable and robust.When the model is divided into smaller parts,the forward modeling efficiency of multi-transmitter and multi-frequency are still to be improved.Twodimensional numerical simulations are calculated at sources or frequencies,the calculations between the sources or between frequencies are independent on each other,so the MPI parallel computing is introduced to improve the calculation efficiency,to reduce the time needed for the forward modeling.In this paper,after calculating the same synthetic model by serial and parallel codes,the results show that the parallel algorithm is robust.By analyzing the efficiency of parallel algorithm,the parallel algorithm is reliable and efficient.
引文
[1]杨进.环境与工程地球物理[M].北京:地质出版社,2011.YANG J.Environmental and engineering geophysics[M].Beijing:Geological Publishing House,2011.(In Chinese)
[2]杨振威,严加永,陈向斌.频谱激电法在安徽沙溪斑岩铜矿中的应用[J].地球物理学进展,2013,28(4):2014-2023.YANG Z W,YAN J Y,CHEN X B.The application of spectral induced polarization in Shaxi porphyry copper in Anhui province[J].Progress in Geophysics,2013,28(4):2014-2023.(In Chinese)
[3]曹蔚杰,单明良,高勇,等.复电阻率法(CR)在铜钼矿勘查中的应用效果[J].工程地球物理学报,2014,11(2):203-207.CAO W J,SHAN M L,GAO Y,et al.The application of complex resistivity method(CR)to coppermolybdenum Mine exploration[J].Chinese Journal of Engineering Geophysics,2014,11(2):203-207.(In Chinese)
[4]徐自生,张丽,唐伟,等.复电阻率法(CR)在内蒙古那吉河铅锌矿探测中的应用[J].矿产勘查,2015(2):165-170.XU Z S,ZHANG L,TANG W,et al.Application of complex resistivity method(CR)in prospecting leadzinc deposit in Inner Mongolia[J].Mineral Exploration,2015(2):165-170.(In Chinese)
[5]REVIl A,KARAOULIS M,JOHNSON T,et al.Review:some low-frequency electrical methods for subsurface characterization and monitor in hydrogeology[J].Hydrogeology Journal,2012,20(4):617-658.
[6]许传建,徐自生,杨志成,等.复电阻率(CR)法探测油气藏的应用效果[J].石油地球物理勘探,2004,A:31-35.XU C J,XU Z S,YANG Z C,et al.The application of complex resistivity method(CR)to oil gas exploration[J].Oil Geophysical Prospecting,2004,A:31-35.(In Chinese)
[7]FLORES OROZCO A,KEMNA A,OBERDRSTER C,et al.Delineation of subsurface hydrocarbon contamination at a former hydrogenation plant using spectral induced polarization imaging[J].Journal of Contaminant Hydrology,2012,136-137:131-144.
[8]BREEDE K,KEMNA A,ESSER O,et al.Spectral induced polarization measurements on variably saturated sand-clay mixtures[J].Near surface geophysics,2012,10(6):479-489.
[9]张武生,薛巍,李建江,等.MPI并行程序设计实例教程[M].北京:清华大学出版社,2009.ZHANG W S,XUE W,LI J J,et al.Parallel programming with MPI[M].Beijing:Tsinghua University Press,2009.(In Chinese)
[10]都志辉,李三立,陈渝,等.高性能计算之并行编程技术——MPI并行程序设计[M].北京:清华大学出版社,2001.DU Z H,LI S L,CHEN Y,et al.High performance parallel computing technology——MPI parallel programming[M].Beijing:Tsinghua University Press,2001.(In Chinese)
[11]NABIGHIAN M.N.Electromagnetic Methods in applied geophysics[M].Theory.Society of Exploration Geophysicists,1988.
[12]徐世浙.地球物理中的有限单元法[M].北京:科学出版社,1994.XU S Z.Finite element method in geophysics[M].Beijing:Science Press,1994.(In Chinese)
[13]张斌.可控源音频大地电磁二维正反演研究[D].北京:中国地质大学,2012.ZHANG B.Research of 2D CSAMT forward and inversion[D].Beijing:China University of Geosciences(Beijing),2012.(In Chinese)
[14]张志勇.复电阻率法二维有限元正演研究[D].北京:中国地质大学,2013.ZHANG Z Y.2DFEM forward modeling research of complex resistivity method[D].Beijing:China University of Geosciences(Beijing),2013.(In Chinese)
[15]王寒冰.极低频电磁法频率域二维正反演研究[D].北京:中国地质大学,2014.WANG H B.Research of 2Dfrequency-domain ELF forward and inversion[D].Beijing:China University of Geosciences,2014.(In Chinese)
[16]桂兵.可控源音频大地电磁法张量数据二维反演研究[D].北京:中国地质大学,2015.GUI B.Research of 2DCSAMT tensor data inversion[D].Beijing:China University of Geosciences,2015.(In Chinese)
[17]PELTON W H,WARD S H,HALLOF P G,et al.Mineral discrimination and removal of inductive coupling with multifrequency IP[J].Geophysics,1978,43(3):588-609.
[18]林绍忠.对称逐步超松弛预处理共轭梯度法的改进迭代格式[J].数值计算与计算机应用,1997(12):266-270.LIN S Z.Improved iterative format of symmetric successive over relaxation-preconditioned conjugated gradient method[J].Journal of Numerical Methods and Computer Applications,1997(12):266-270.(In Chinese)
[19]谭捍东.大地电磁法三维正反演问题研究[D].北京:中国地质大学,2000.TAN H D.A study on 3d MT forward and inversion[D].Beijing:China University of Geosciences,2000.(In Chinese)
[20]胡祥云,李焱,杨文采,等.大地电磁三维数据空间反演并行算法研究[J].地球物理学报,2012,55(12):3969-3978.HU X Y,LI Y,YANG W C,et al.Three-dimensional magnetotelluric parallel inversion using data space method[J].Chinese Journal of Geophysics,2012,55(12):3969-3978.(In Chinese)
[21]林昌洪.三维大地电磁共轭梯度反演算法研究[D].北京:中国地质大学,2004.LIN C H.A study on 3d MT non-linear conjugate gradient inversion[D].Beijing:China University of Geosciences,2004.(In Chinese)