|
|
Research on the Two-dimensional Magnetotelluric Regularization Inversion |
YIN Yin-xuan1,ZHANG Lian-wei1 , TONG Xiao-zhong2 |
1 .The Third Railway Survey and Design Institute Group Corporation,Tianjin 300142,China; 2. Central South University,Changsha, Hunan 410083, China |
|
|
Abstract Research purposes : In the study of magneto telluric inversion methods,models and observational data are generally non-linear, limited and noisy observation data makes the inversion results non-unique, approximate geophysical modeling theory and lack of priori information makes it difficult to solve the problem of global convergence and convergence efficiency simultaneously by all inversion methods. A stable solution of the ill-posed inverse problem can be obtained by utilizing the regularization methods in the objective function. Therefore,it is particularly important to find a suitable dimensional regularization inversion method.
Research conclusions:(1)Solving large scale linear equation of inverse problem,the damped Gauss-Newton algorithm was adopted,which can improve local convergence of Gauss-Newton method. Through the synthetic model simulation, the inversion results based on damped Gauss-Newton algorithm truly reflected the ge。-electrical parameters of the model and accurately showed the depth and size of the abnormal body,So the inversion algorithm has been proven to be effective.(2)Compared with the actual observation data, whether the apparent resistivity data predicted or the impedance phase data,are basically consistent. So the inversion algorithm is accuracy for the applying in two- dimensional magneto telluric data interpretation.(3)Using GCV method can effectively determine the optimum regularization parameter, the ill-posed problem is transformed into well-posed problem for solving , keeping the balance between the data objective function and the objective function constrained by model,the model variance and the image resolution can achieve the best compromise. (4) The inversion method is suitable for two-dimensional magneto telluric forward modeling under the condition of approximate level terrain and can be used to analysis two-dimensional magneto telluric data in the field of railways,highway or other engineering survey.
|
Received: 29 November 2013
|
|
|
|
|
[1]王家映.地球物理反演理论[M].武汉:中国地质大学出 版社,1998. [2]
Wang Jiaying. Inverse Theory in Geophysics[M]. Wuhan; China University of Geosciences Press, 1998.
[2] 陈乐寿,刘任,王天生.大地电磁测深资料处理与解释 [M].北京:石油工业出版社,1989.
Chen Leshou, Liu Ren, Wang Tiansheng. The Magneto telluric Sounding Data Processing and Interpretation Methods[M〕.Beijing; Petroleum Industry Press, 1989.
[3] 严良俊,胡文宝.大地电磁测深资料的二次函数逼近非线性反演[J].地球物理学报,2004(5);935-940.
Yan Liangjun, Hu Wenbao. Non-linear Inversion with the Quadratic Function Approaching Method for Magneto telluric Data[J ]Chinese Journal of Geophysics , 2004 ( 5 ) ; 935-940.
[4] 戴世坤,徐世浙.MT二维和三维连续介质快速反演[J]. 石油地球物理勘探,1997(3);305-317.
Dai Shikun, Xu Shizhe. Rapid Inversion of Magnetotelluric Data for 2-D and 3-D Continuous Media[J〕.Oil Geophysical Prospecting , 1997(3): 305-317.
[5] 柳建新,童孝忠,杨晓弘,等实数编码遗传算法在大地电磁测深二维反演中的应用[J].地球物理学进展,2008 (6):1936-1942.
Liu Jiamdn, Tong Xiaozhong, Yang Xiaohong, etc. Application of Real Coded Genetic Algorithm in Two- dimensional Magnetotelluric Inversion[J].Progress in Geophysics , 2008 ( 6 ) :1936-1942.
[6] Rodi W,Mackie R L. Nonlinear Conjugate Gradients Algorithm for 2-D Magneto telluric Inversion[J〕. Geophysics , 2001(1);174-187.
[7] 余年,李坚,雷旭友,等.大地电磁约束反演在隧道勘探中的研究[J].铁道工程学报,2012(7);41 -46.
Yu Nian , Li Jian , Lei Xuyou , etc. Study on Application of Magnetotelluric Constrained Inversion in Tunnel Exporation [ J ].Journal of Railway Engineering Society , 2012 ( 7);41-46.
[8]Dennis J E , Schnabel R B. Numerical Methods for Unconstrained Optimization and Nonlinear Equations [M].Philadelphian; Society for Industrial and Applied Mathematics , 1996.
[9] 柳建新,蒋鹏飞,童孝忠,等.不完全LU分解预处理的 BICGSTAB算法在大地电磁二维正演模拟中的应用LJ}. 中南大学学报(自然科学版),2009 ( 2 ) ; 484 - 491.
Liu Jianxin , Jiang Pengfei , Tong Xiaozhong , etc. Application of BICGSTAB Algorithm with Incomplete LU Decomposition Preconditioning to Two- dimensional Magneto telluric Forward Modeling[J].Journal of Central South University(Science and Technology), 2009(2);484-491.
|
|
|
|