陈先洁, 王绪本, 李德伟, 谢卓良, 乃国茹

(成都理工大学 地球物理学院，成都 610000)

0 引言

在大地电磁测深法中往往存在由于浅地表电性结构不均匀，使得在不均匀体的周围电流密度发生畸变，从而使电场分量发生变化，引起局部的畸变异常[1-4]。为消除此种局部畸变异常带来的影响，Groom和Bailey假设局部三维异常体覆盖在二维区域异常之上，将实际阻抗张量分解为畸变阻抗张量和区域阻抗张量，以此来恢复未受畸变影响的区域阻抗张量(GB分解法)。它的本质是求解由区域阻抗张量和畸变方程组成的非线性超定方程组[6-8]。传统的GB分解法将非线性问题线性化求解，导致求解结果不稳定，易陷入局部极值中，很难得到较为准确的分解结果。李洋等[5]提出来一种混合优化算法，利用模拟退火法得到的全局最优解作为局部非线性最小二乘法的初始值，提高了求解稳定性；尹曜田等[6]研究了基于遗传算法的阻抗张量分解方法，采用遗传算法求解非线性方程组；蔡军涛等[7]提出共轭阻抗变换的概念，不受电场畸变影响，无需地下电性维性信息的假设；谢成良等[8]研究了基于相位张量约束下的阻抗张量分解方法，利用相位张量不受电场畸变影响的优点，分析得到合理的初始模型，提高了计算稳定性和计算效率，但仍然存在对初始值依赖性严重等问题。

笔者在求解GB分解法中的非线性超定方程组时,采用了粒子群优化算法，但传统的粒子群算法极易陷入局部极值当中，无法满足求解要求。……