Two-dimensional inversion of spectral induced polarization data using MPI parallel algorithm in data space  

复电阻率法二维数据空间反演并行算法研究（英文）

作　　者：张志勇[1,2] 谭捍东[1,2] 王堃鹏[1,2] 林昌洪[1,2] 张斌[3] 谢茂笔 

机构地区：[1]中国地质大学(北京)地球物理与信息技术学院,北京100083 [2]中国地质大学地下信息探测技术与仪器教育部重点实验室,北京100083 [3]有色金属矿产地质调查中心,北京100012

出　　处：《Applied Geophysics》2016年第1期13-24,217,共13页应用地球物理（英文版）

基　　金：jointly sponsored by the National Natural Science Foundation of China(Grant No.41374078);the Geological Survey Projects of the Ministry of Land and Resources of China(Grant Nos.12120113086100 and 12120113101300);Beijing Higher Education Young Elite Teacher Project

摘　　要：Traditional two-dimensional（2D） complex resistivity forward modeling is based on Poisson＇s equation but spectral induced polarization（SIP） data are the coproducts of the induced polarization（IP） and the electromagnetic induction（EMI） effects.This is especially true under high frequencies,where the EMI effect can exceed the IP effect.2D inversion that only considers the IP effect reduces the reliability of the inversion data.In this paper,we derive differential equations using Maxwell＇s equations.With the introduction of the Cole-Cole model,we use the finite-element method to conduct2 D SIP forward modeling that considers the EMI and IP effects simultaneously.The data-space Occam method,in which different constraints to the model smoothness and parametric boundaries are introduced,is then used to simultaneously obtain the four parameters of the Cole-Cole model using multi-array electric field data.This approach not only improves the stability of the inversion but also significantly reduces the solution ambiguity.To improve the computational efficiency,message passing interface programming was used to accelerate the 2D SIP forward modeling and inversion.Synthetic datasets were tested using both serial and parallel algorithms,and the tests suggest that the proposed parallel algorithm is robust and efficient.传统的复电阻率法二维正演是基于泊松方程完成的,而实际的频谱激电数据包括激电效应和电磁效应双重作用的影响,特别是频率较高时电磁效应甚至大于激电效应,因此只考虑激电效应的做法会降低反演结果的可靠性。本文的研究从麦克斯韦方程出发,引入Cole-Cole模型,同时考虑了激电效应和电磁效应,采用有限单元法实现了复电阻率法二维数值模拟。反演采用二维数据空间OCCAM方法,通过施加不同的模型光滑约束和参数界限约束,利用多排列电场数据同时反演Cole-Cole模型四个参数,既提高了反演的稳定性又降低了反演的多解性。为了提高计算效率,我们选用MPI实现了频谱激电二维正反演并行算法。通过理论模型的试算,对比串行程序和并行程序的计算结果和效率,验证了并行算法的稳定性、可靠性和高效性。

关 键 词：Spectral induced polarization 2D inversion data-space method Cole-Cole model MPI parallel computation 

分 类 号：P631.322[天文地球—地质矿产勘探]