| 用改进的光滑约束最小二乘正交分解法实现电阻率三维反演 |
| |
| 引用本文： | 宛新林,席道瑛,高尔根,沈兆武.用改进的光滑约束最小二乘正交分解法实现电阻率三维反演[J].地球物理学报,2005,48(2):439-444. |
| |
| 作者姓名： | 宛新林 席道瑛 高尔根 沈兆武 |
| |
| 作者单位： | 中国科学技术大学地球与空间科学学院,合肥,230026;中国科学技术大学工程科学学院,合肥,230026;安徽建筑工业学院土木工程系,合肥,230022;中国科学技术大学地球与空间科学学院,合肥,230026;中国科学技术大学工程科学学院,合肥,230026 |
| |
| 基金项目： | 国家自然科学基金项目 (4 0 1740 5 0 )资助 |
| |
| 摘 要： | 对三维电阻率反演问题进行了深入研究，提供了一种利用地表观测数据实现三维反演的实用算法.该方法应用有限差分求正演解，并通过对粗糙度矩阵元素进行适当改进，使之适用于各种情况下粗糙度矩阵的求取，进而建立在模型的总粗糙度极小条件下的反演方程.对反演方程采用收敛速度快且稳定的最小二乘正交分解（LSQR）法进行迭代求解，在迭代求解过程中只需利用偏导数矩阵和其转置矩阵乘以一个向量的结果，回避了直接求偏导数矩阵的繁琐计算，节省了内存，加快了反演的计算速度.不同的计算实例表明上述方法是求解大规模三维电阻率反演问题的有效方法.
|
| 关 键 词： | 三维反演 偏导数矩阵 光滑约束 LSQR算法 |
| 文章编号： | 0001-5733(2005)02-0439-06 |
| 收稿时间： | 2003-09-12 |
| 修稿时间： | 2004-09-28 |
3-D resistivity inversion by the least-squares QR factorization method under improved smoothness constraint condition |
| |
| WAN Xin-Lin,XI Dao-ying,GAO Er-Gen,SHEN Zhao-Wu.3-D resistivity inversion by the least-squares QR factorization method under improved smoothness constraint condition[J].Chinese Journal of Geophysics,2005,48(2):439-444. |
| |
| Authors: | WAN Xin-Lin XI Dao-ying GAO Er-Gen SHEN Zhao-Wu |
| |
| Institution: | 1.College of Earth and Space Science, University of Science and Technology of China, Hefei 230026,China 2 College of Engineering Science, University of Science and Technology of China, Hefei 230026,China 3 Department of Civil Engineering, Anhui Institute |
| |
| Abstract: | In this paper, we have deeply studied 3-D resistivity inversion, providing a practical algorithm of \{3-D\} inversion. Using the finite difference method to solve the 3-D forward problem, we have properly improved the elements of roughness matrix in order to form the roughness matrix in different cases, setting up the inversion equation under the condition that the total roughness of the model is minimum. It is fast and stable that using the least-squares QR factorization(LSQR) algorithm to solve inversion equations. The LSQR algorithm only requires the result of the derivative matrix and its transpose multiplying vectors, therefore we avoid direct and complicated computations of the derivative matrix. The above approach reduces the need for computer memory, and speeds up the inversion calculation. Two different calculation examples show that the above approach is efficacious for solving large scale 3-D resistivity inversion. |
| |
| Keywords: | 3-D inversion Partial derivative matrix Smoothness constraint LSQR algorithm |
| 本文献已被 CNKI 维普 万方数据 等数据库收录！ |
| 点击此处可从《地球物理学报》浏览原始摘要信息 |
| 点击此处可从《地球物理学报》下载免费的PDF全文 |
|
