改进牛顿–拉夫逊电阻层析成像图像重建算法

为了满足电阻层析成像图像重建算法应用于两相流领域的精度与实时性要求，提出一种改进牛顿–拉夫逊算法。针对牛顿–拉夫逊算法局部收敛的特点，将线性反投影算法成像结果作为新算法初始电阻率分布估计值，并利用基于改进粒子群算法的平衡法对海森矩阵进行预处理，降低条件数。同时在新算法迭代过程中，当振荡次数超过设定值时，采取自动更新灵敏度矩阵的策略。仿真与实验结果表明，与线性反投影算法、兰德韦伯算法及简化牛顿–拉夫逊算法相比，新算法降低了图像重建误差，有效提高了图像重建质量；而与标准牛顿–拉夫逊算法相比，新算法在满足成像精度前提下，避免了每次迭代均需更新灵敏度矩阵的缺点，提高了算法收敛速度。