【摘要】 传统的复电阻率法二维正演是基于泊松方程完成的,而实际的频谱激电数据包括激电效应和电磁效应双重作用的影响,特别是频率较高时电磁效应甚至大于激电效应,因此只考虑激电效应的做法会降低反演结果的可靠性。本文的研究从麦克斯韦方程出发,引入Cole-Cole模型,同时考虑了激电效应和电磁效应,采用有限单元法实现了复电阻率法二维数值模拟。反演采用二维数据空间OCCAM方法,通过施加不同的模型光滑约束和参数界限约束,利用多排列电场数据同时反演Cole-Cole模型四个参数,既提高了反演的稳定性又降低了反演的多解性。为了提高计算效率,我们选用MPI实现了频谱激电二维正反演并行算法。通过理论模型的试算,对比串行程序和并行程序的计算结果和效率,验证了并行算法的稳定性、可靠性和高效性。更多还原

【Abstract】 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.更多还原