MT time series simulation based on constructed system function
Zhang Baoqiang1, Pei Jianxin1,2, Wang Qi3
1. College of Marine Geosciences, Ocean University of China, Qingdao, Shandong 266100, China; 2. Key Laboratory of Submarine Geosciences and Prospecting Techniques, Ministry of Education, Qingdao, Shandong 266100, China; 3. Guangzhou Marine Geological Survey, Guangzhou, Guangdong 510760, China
Abstract:For magnetotelluric (MT) method, an effective geophysical tool, data impedance estimation and noise elimination are key steps, which have great influence on the impedance characteristics.So it is necessary to test impedance estimation and subsequently signal processing in the time domain with simulated electromagnetic time series.Based on constructed system function, the magnetic field in the frequency domain is obtained by the discrete Fourier transform (DFT) of random sequence, and the electric field in the frequency domain is calculated with model impedance and simulated magne-tic field.Finally, the symmetry processing of complex conjugate is performed on the unilateral spectrum of electromagnetic field and the impedance sequence in the frequency domain, and the simulated time series of electromagnetic field is obtained with inverse fast Fourier transform (IFFT).Test results show that the proposed method can be used to simulate the time series of one-dimensional and two-dimensional horizontal heterogeneous model, and compared with measured data, the electromagnetic field data obtained with this method can better reflect MT frequency spectrum characteristics.
张宝强, 裴建新, 王启. 基于构造系统函数的大地电磁时间序列模拟方法[J]. 石油地球物理勘探, 2018, 53(2): 410-417.
Zhang Baoqiang, Pei Jianxin, Wang Qi. MT time series simulation based on constructed system function. Oil Geophysical Prospecting, 2018, 53(2): 410-417.
张全胜,王家映.大地电磁测深资料的去噪方法.石油地球物理勘探,2004,39(增刊1):17-23.Zhang Quansheng,Wang Jiaying.An method of noise elimination for magnetotelluric sounding data.OGP,2004,39(S1):17-23.
[2]
汤井田,张弛,肖晓等.大地电磁阻抗估计方法对比.中国有色金属学报,2013,23(9):2351-2358.Tang Jingtian,Zhang Chi,Xiao Xiao et al.Comparison of methods for magnetotelluric impedance estimation.The Chinese Journal of Nonferrous Metals,2013,23(9):2351-2358.
[3]
蔡剑华,王先春,胡惟文.基于经验模态分解与小波阈值的MT信号去噪方法.石油地球物理勘探,2013,48(2):303-307.Cai Jianhua,Wang Xianchun,Hu Weiwen.De-noising of MT signal based on empirical mode decomposition and wavelet threshold method.OGP,2013,48(2):303-307.
[4]
张翔,严良俊,苏朱刘等.山区MT资料处理新技术及效果.石油地球物理勘探,2007,42(4):454-456.Zhang Xiang,Yan Liangjun,Su Zhuliu et al.New technology and effects of mountainous MT data processing.OGP,2007,42(4):454-456.
[5]
Goubau W M,Gamble T D,Clarke J.Magnetotelluric data analysis:removal of bias.Geophysics,1978,43(6):1157-1166.
[6]
Larsen J C,Mackie R L,Manzella A et al.Robust smooth magnetotelluric transfer functions.Geophysical Journal International,2007,124(3):801-819.
[7]
Yee E.The reconstruction of the magnetotelluric impedance tensor:An adaptive parametric time-domain approach.Geophysics,1988,53(3):1080-1087.
[8]
Filipo W A S,Hohmann G W.Computer simulation of low-frequency electromagnetic data acquisition.Geophysics,1982,48(9):1219-1232.
[9]
Lamarque G.Improvement of MT data processing using stationary and coherence tests.Geophysical Prospecting,1999,47(6):819-840.
[10]
Ernst T.A comparison of two methods of the transfer function calculation using the least-square criterion in time and frequency domain.Publications of Institute of Geophysics Polish Academy of Sciences,1981,143:13-24.
[11]
Kunetz G.Processing and interpretation of magnetotelluric soundings.Geophysics,1972,37(6):1005-1021.
[12]
Mcmechan G A,Barrodale I.Processing electromagnetic data in the time domain.Geophysical Journal International,1985,81(1):277-293.
[13]
Spagnolini U.Time-domain estimation of MT impe-dance tensor.Geophysics,1994,59(5):712-721.
[14]
Chen J.Using empirical mode decomposition to process magnetotelluric data.Christian-Albrechts-Universität,2013.
[15]
Loddo M,Schiavone D,Siniscalchi A.Generation of synthetic wide-band electromagnetic time series.Geophysics,2009,45(2):289-301.
[16]
Varentsov I M,Sokolova E Y.Generation of synthetic magnetotelluric data.Izvestiya Physics of the SolidEarth,1995,30(6):554-562.
[17]
严家斌.大地电磁信号处理理论及方法研究[学位论文].湖南长沙:中南大学,2003.Yan Jiabin.The Study on Theory and Method of Magnetotelluric Signal Processing[D].Central South University,Changsha,Hu'nan,2003.
[18]
蔡剑华,李晋.基于频率域小波去噪的大地电磁信号工频干扰处理.地质与勘探,2015,51(2):353-359.Cai Jianhua,Li Jin.Suppression of power line interfe-rence on MT signals based on the frequency domain wavelet method.Geology and Exploration,2015,51(2):353-359.
[19]
陈知富,邓居智,陈辉等.基于小波变换的大地电磁信号去噪研究.工程地球物理学报,2012,9(6):732-737.Chen Zhifu,Deng Juzhi,Chen Hui et al.De-noising research of magnetotelluric signal based on wavelet transform.Chinese Journal of Engineering Geophy-sics,2012,9(6):732-737.
[20]
凌振宝,王沛元,万云霞等.强人文干扰环境的电磁数据小波去噪方法研究.地球物理学报,2016,59(9):3436-3447.Ling Zhenbao,Wang Peiyuan,Wan Yunxia et al.A combined wavelet transform algorithm used for de-noising magnetotellurics data in the strong human noise.Chinese Journal of Geophysics,2016,59(9):3436-3447.
[21]
Bracewell R,Kahn P B.The Fourier Transform and Its Applications.McGraw-Hill,2000.
[22]
Rao K R,Kim D N,Hwang J J.Fast Fourier Transform-Algorithms and Applications.Springer Netherlands,2010.
[23]
程乾生.数字信号处理.北京:北京大学出版社,2010.Chen Qiansheng.Digital Signal Processing.Peking University Press,Beijing,2010.
[24]
Labson V F,Becker A,Morrison H F et al.Geophysical exploration with audiofrequency natural magnetic fields.Geophysics,1985,50(4):656-664.
[25]
Serson P H.Instrumentation for induction studies on land.Physics of the Earth & Planetary Interiors,1973,7(3):313-322.
[26]
Filloux J H.Techniques and instrumentation for study of natural electromagnetic induction at sea.Physics of the Earth & Planetary Interiors,1973,7(3):323-338.
[27]
Oppenheim A V,Schafer R W.Digital Signal Processing.Prentice-Hall,New Jersey,1975.