本文将Tikhonov正则化方法与activeset算法相结合,利用双频电磁波电导率成像原理,求解其反演成像方程.不仅对现有算法进行了改进,也促进了算法的实际应用.本文研究了在双频电磁波电导率成像方程建立后,如何根据其严重病态性质,选择合适的算法求解矩阵成像方程.针对电导率非负的特性,引入正则化参数,将问题转化为一个非负最小二乘问题,并用active—set算法求解.采用改进后的迭代算法对理论模型进行了数值模拟计算,验证了该方法的有效性.应用到实际电导率成像反演,与常规的LSQR、SP-LSQR、Tikhonov正则化等算法进行比较,取得了满意的结果.
Using a combination of Tikhonov regularization and active-set method, this paper introduces an algorithm to numerically resolve electrical conductivity imaging equation using dualfrequency EM data, so that we develop not only the existing algorithms, but also its implementation. The linear-equation system, however, generated from dual-frequency electrical conductivity imaging, is exceedingly ill posed, and this paper just focuses on its resolution. Considering the non-negativity of the electrical conductivity, we introduce a proper regularization parameter to reduce the problem to a non negative least square problem, and then solve it with active-set algorithm. Afterwards some numerical examples are carried out to test the validity of the proposed method. And finally, we implement it to the actual inversion problem of conductivity imaging using dual-frequency EM data, and gain satisfactory results compared to normal methods such as LSQR, SP-LSQR and shear Tikhonov regularization.