THE BATCH PROCESSING TECHNIQUE FOR MAGNETOTELLURIC DATA WITH SSMT-2000 BASED ON MATLAB
GUO Ze-qiu1,2, YE Gao-feng1,2,3, WANG Hui1,2, WEI Wen-bo1,2,3
1. School of Geophysics and Information Technology, China University of Geosciences, Beijing 100083, China; 2. Geo-detection and Instruments Laboratory, Ministry of Education, Beijing 100083, China; 3. State Key Laboratory of Geological Processes and Mineral Resources, Beijing 100083, China
Abstract:Data-processing plays an important role in the magnetotelluric (MT) sounding work. For several simultaneous sounding sites, in particular, it is required to process field data with magnetic remote reference, which is time-consuming and demands a huge amount of calculation. The professional processing software named SSMT-2000 is commonly used to carry out the work; nevertheless, because of the restriction of the software, the batch processing for data from several sites can't be carried out efficiently. Hence data-processing has become a tedious and long time work for MT workers. In this paper, the authors analyzed the format and the calling pattern of parameter files, probed into the problem as to how the internal executables are called by SSMT-2000 and, on such a basis, developed a smart code based on the MATLAB language. The code has the utilities of running a batch processing for fast Fourier transform, Robust estimation of impedance tensors for one single site and utilization of magnetic remote reference technique in three ways. The authors can apply the code to process data from many MT sounding sites at one time, which makes multiple-site data processing procedures much simpler and more efficient.
郭泽秋, 叶高峰, 王辉, 魏文博. 基于MATLAB的SSMT-2000大地电磁资料批处理技术[J]. 物探与化探, 2013, 37(5): 921-925.
GUO Ze-qiu, YE Gao-feng, WANG Hui, WEI Wen-bo. THE BATCH PROCESSING TECHNIQUE FOR MAGNETOTELLURIC DATA WITH SSMT-2000 BASED ON MATLAB. Geophysical and Geochemical Exploration, 2013, 37(5): 921-925.
[1] Cagniard L.Basic theory of the magnetotelluric method of geophysical prospecting[J].Geophysics,1953,18(3):605-635. [2] Fiona Simpson,Karsten Bahr.Practical Magnetotellurics[M].Cambridge University Press,2005:165-180. [3] 魏文博.我国大地电磁测深新进展及瞻望[J].地球物理学进展,2002,17(2):245-254. [4] 赵国泽,陈小斌,汤吉.中国地球电磁法新进展和发展趋势[J].地球物理学进展,2007,22(4):1171-1180. [5] 仇根根,方慧,朱正君,等.基于MATLAB平台对V5-2000大地电磁测深仪TBL文件的分析与应用[J].物探化探计算技术,2012,34(2):178-181. [6] 柳建新,刘春明,马捷,等.V5-2000大地电磁测深仪文件头数据格式研究[J].物探化探计算技术,2007,29(4):359-362. [7] 陈乐寿,刘任,王天生.大地电磁资料处理与解释[M].北京:石油工业出版社,1989:28-31. [8] Sims W E,Bostick F X,Smith H W.The estimation of magnetotelluric impedance tensor elements from measured data[J].Geophysics,1971,36(5):938-942. [9] Gamble T D,Goubau W M,Clarke J.Magnetotellurics with a remote magnetic reference[J].Geophysics,1979,44(1):53-68. [10] Gamble T D,Goubau W M,Clarke J.Error analysis for remote reference magnetotellurics[J].Geophysics,1979,44(5):959-968. [11] Jones A G,Jodicke H.Magnetotelluric transfer function estimation improvement by a coherence-based rejection technique[J].SEG Expanded Abstracts,1984,3(1):51-55. [12] Egbert G D,Livelybrooks D W.Single station magnetotelluric impedance estimation:Coherence weighting and the regression M-estimate[J].Geophysics,1996,61(4):964-970. [13] Egbert G D,Booker J R.Robust estimation of geomagnetic transfer function[J].Geophysical Journal International,1986,87(1):173-194. [14] Egbert G D.Robust multiple-station magnetotelluric data processing[J].Geophysical Journal International,1997,130(2): 475-496. [15] Larsen J C.Transfer functions: smooth robust estimates by least-squares and remote reference methods[J].Geophysical Journal International,1989,99(3):645-663.