基于不等式约束的最小二乘法三维电阻率反演及算法优化.pdfVIP

第 卷 第 期 地 球 物 理 学 报 ， ５５ １ Ｖｏｌ．５５ Ｎｏ．１ 　 　 　 　 　 年 月 ， ２０１２ １ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ｊａｎ．２０１２ 　 　 　 刘斌，李术才，李树忱等 基于不等式约束的最小二乘法三维电阻率反演及其算法优化 地球物理学报， ， （）： ， ． ． ２０１２５５１ ２６０２６８ ： ／ ｄｏｉ１０．６０３８ ．ｉｓｓｎ．０００１５７３３．２０１２．０１．０２５． ｊ ， ， ， ＬｉｕＢ ＬｉＳＣ ＬｉＳＣ ｅｔａｌ．３Ｄｅｌｅｃｔｒｉｃａｌｒｅｓｉｓｔｉｖｉｔｉｎｖｅｒｓｉｏｎｗｉｔｈｌｅａｓｔｓｕａｒｅｓｍｅｔｈｏｄｂａｓｅｄｏｎｉｎｅｕａｌｉｔｃｏｎｓｔｒａｉｎｔａｎｄｉｔｓ ｙ ｑ ｑ ｙ （ ）， ， （）： ， ： ／ ｃｏｍ ｕｔａｔｉｏｎｅｆｆｉｃｉｅｎｃｏｔｉｍｉｚａｔｉｏｎ．犆犺犻狀犲狊犲犑．犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ２０１２ ５５１ ２６０２６８ｄｏｉ１０．６０３８ ．ｉｓｓｎ．０００１５７３３． ｐ ｙ ｐ 狆狔 ｊ ２０１２．０１．０２５． 基于不等式约束的最小二乘法三维电阻率反演 及其算法优化 刘 斌，李术才，李树忱，聂利超，钟世航，李利平，宋 杰，刘征宇 　 　 山东大学岩土与结构工程研究中心，济南 ２５００６１ 　 摘 要 基于光滑约束的最小二乘法是三维电阻率反演的主要方法，但该方法在某些情况下存在着多解性较强的 　 　 问题，且普遍耗时较长，严重制约了三维反演方法的推广与发展 为改善上述问题，将表征模型参数变化范围的不 ． 等式约束作为先验信息引入最小二乘线性反演方法中，有效地改善了反演结果的精度，降低了反演的多解性问题． 为了解决耗时较长的问题，基于预条件共轭梯度（ ）算法和 分解法的特点提出了一套优化三维电阻率 ＰＣＧ Ｃｈｏｌｅｓｋｙ 反演计算效率的计算方案 在该方案中， 分解法被用来求解敏感度矩阵计算中的多个点源场的正演问题， ． Ｃｈｏｌｅｓｋｙ 分解法只需对总体系数矩阵进行一次分解，然后对不同的右端向量进行回代即可 将预条件共轭梯度法 Ｃｈｏｌｅｓｋｙ ． 引入到三维电阻率反演方程的求解中，将雅可比迭代中的对角阵作为预处理矩阵，其具有求逆方便、无需内存空间 的特点，有效地加快了收敛速