|
|
Analytical algorithms of Jacobian matrix in one-dimensional Occam inversion of central loop transient electromagnetism based on all-time apparent ressivivity |
Song-Wei GUO1, Xiao-Pan LIU2, Kai ZHENG3, Lei ZHANG4 |
1. School of Earth Sciences and Resources, China University of Geosciences (Beijing) , Beijing 100083, China 2. Institute of Geology and Mineral Resources Exploration of Zhejiang Province, Hangzhou 310000, China 3. College of Geophysics and Petroleum Resources, Yangtze University, Wuhan 430100, China 4. Inner Mongolia Land Resources Exploration and Development Co., Ltd., Hohhot 010020, China |
|
|
Abstract In the current engineering application, one-dimensional inversion result of TEM sounding is still an important basis for geological interpretation. In order to improve the computational efficiency and accuracy of one-dimensional inversion, it is necessary to analytically derive Jacobian matrix of one-dimensional forward. In this paper, the analytical formulas of Jacobian matrix in one-dimensional forward of central loop transient electromagnetic based on all-time resisitivity are derived in detail. It is proved that the formulas are correct and worth adopting, as shown by testing one-dimensional Occam inversion of theoretical model and measured data.
|
Received: 08 July 2019
Published: 24 June 2020
|
|
|
|
|
|
Comparison of H-type layered Model inversion results
|
| 解析算法 | 差分算法 | 拟合差 | 5.18 | 169.5 | 对数拉格朗日乘子 | -0.33 | -2.03 | 粗糙度 | 5956.51 | 18101.32 | 调用正演次数 | 56 | 93 | 计算时间/s | 28.34 | 44.45 |
|
Inversion result parameters of H-type layered model
|
|
Comparison of K-type layered Model inversion results
|
| 解析算法 | 差分算法 | 拟合差 | 2.23 | 96.54 | 对数拉格朗日乘子 | -0.02 | -5.75 | 粗糙度 | 195.40 | 3229.25 | 调用正演次数 | 28 | 91 | 计算时间/s | 16.60 | 48.69 |
|
Inversion result parameters of K-type layered model
|
|
Geological map of work area 1—middle low mountain; 2—piedmont inclined plain; 3—Quaternary Holocene alluvial proluvial; 4—Jurassic volcanic rock; 5—Archaean Sanggan group gneiss; 6—Archaean porphyry like granite; 7—fracture fracture Fracture Zone; 8—test TEM sounding line; 9—geological boundary; 10—geomorphic boundary; 11—existing drilling;12—verification of geothermal drilling (DR2); 13—TEM emission line frame position; 14—TEM receiving coil position
|
|
Borehole histogram of work area
|
|
Inversion result of site 820 where is projectionposition of borehole
|
|
Comparison of TEM inversion sections of line 400 a—pseudo section; b—smoke ring inversion section; c—Occam inversion section
|
[1] |
牛之琏. 时间域电磁法原理[M]. 长沙: 中南大学出版社, 2007.
|
[1] |
Niu Z L. Theory of time domain electromagnetic method[M]. Changsha: Press of Central South University, 2007.
|
[2] |
蒋邦远. 实用近区磁源瞬变电磁法勘探[M]. 北京: 地质出版社, 1998.
|
[2] |
Jiang B Y. Practical area near magnetic source transient electromagnetic method prospecting[M]. Beijing: Geology Publishing House, 1998.
|
[3] |
米萨克·N·纳比吉安. 勘查地球物理-电磁法[M].赵经祥译. 北京: 地质出版社, 1992: 195-209.
|
[3] |
Misac N N. Electromagnetic method in applied geophysics[M]. Beijing: Geological Publishing House, 1992: 195-209.
|
[4] |
Guptasarma D, Singh B. New digital linear filters for Hankel J0 and J1 transforms[J]. Geophysical Prospecting, 1997,45:745-762.
|
[5] |
王华军. 正余弦变换的数值滤波算法[J]. 工程地球物理学报, 2004,1(4):329-335.
|
[5] |
Wang H J. Digital filter algorithm of the sine and cosine transform[J]. Chinese Journal of Engineering Geophysics, 2004,1(4):329-335.
|
[6] |
杨云见, 何展翔, 王绪本. 直流电测深法与中心回线瞬变电磁法联合反演[J]. 物探与化探, 2008,32(4):442-444.
|
[6] |
Yang Y J, He Z X, Wang X B. A discussion of the inversion of combined vertical electrical sounding and central-loop TEM data[J]. Geophysical and Geochemical Exploration, 2008,32(4):442-444.
|
[7] |
Constable S C, Parker R L, Catherine G. Constable Occam’ sinversion: apractical algorithm for generating smooth models from electromagnetic sounding data[J]. Geophysics, 1987,52(3):289-300.
|
[8] |
李斌, 赵虎, 钟邱平, 等. 公路深埋隧道音频大地电磁法勘察中的反演方法简介及应用效果分析[J]. 西南公路, 2018(2):6-12.
|
[8] |
Li B, Zhao H, Zhong Q P, et al. Study on inversion and application of AMT in deeply-buried tunnel exploration[J]. Outhwest Highway, 2018(2):6-12.
|
[9] |
冯思臣, 王绪本, 阮帅. 一维大地电磁测深几种反演算法的比较研究[J]. 石油地球物理勘探, 2004,39(5):594-599.
|
[9] |
Feng S C, Wang X B, Ruan S. Comparison among several inversion algorithm of 1D MT[J]. Oil geophysical prospecting, 2004,39(5):594-599.
|
[10] |
郭睿, 武欣, 肖攀, 等. 直升机航空瞬变电磁Occam反演及软件实现[J]. 电子测量技术, 2017,40(9):169-173.
|
[10] |
Guo R, Wu X, Xiao P, et al. Helicopter transient electromagnetic occam inversion and software development[J]. Electronic measurement technology, 2017,40(9):169-173.
|
[11] |
陈小斌, 赵国泽, 汤吉, 等. 大地电磁自适应正则化反演算法[J]. 地球物理学报, 2005,48(4):937-946.
|
[11] |
Chen X B, Zhao G Z, Tang J, et al. An adaptive regularized inversion algorithm for magnetotelluric data[J]. Chinese Journal of Geophysics, 2005,48(4):937-946.
|
[12] |
毛立峰, 王绪本, 陈斌. 直升机航空瞬变电磁自适应正则化一维反演方法研究[J]. 地球物理学进展, 2011,26(1):300-305.
|
[12] |
Mao L F, Wang X B, Chen B. Study on an adaptive regularized 1D inversion method of helicopter TEM data[J]. Progress in Geophysics, 2011,26(1):300-305.
|
[13] |
李刚, 潘和平, 王智, 等. 回线源瞬变电磁法一维反演算法[J]. 煤田地质与勘探, 2017,45(5):161-166,172.
|
[13] |
Li G, Pan H P, Wang Z, et al. One-dimensional inversion for loop source transient electromagnetic method[J]. Coal Geology & Exploration, 2017,45(5):161-166,172.
|
[14] |
白登海, Maxwell A Meju, 卢健, 等. 时间域瞬变电磁法中心回线方式全程视电阻率的数值计算[J]. 地球物理学报, 2003,46(5):698-704.
|
[14] |
Bai D H, Maxwell A Meju, Lu J, et al. Numerical calculation of all-time apparent resistivity for the central loop transient electromagnetic method[J]. Chinese journal of geophysics, 2003,46(5):698-704.
|
[15] |
郭嵩巍, 王绪本. 瞬变电磁全区视电阻率数值计算方法研究[J]. 物探化探计算技术, 2010,32(5):500-507,455.
|
[15] |
Guo S W, Wang X B. Study on numerical calculation methods of transient electromagnetic all-time apparent resistivity[J]. Computing Technique for Geophysical and Geochemical Exploration, 2010,32(5):500-507,455.
|
[16] |
内蒙古地质矿产局. 内蒙古自治区岩石地层[M]. 武汉: 中国地质大学出版社, 1996.
|
[16] |
Inner Mongolia Bureau of Geology and Mineral Resources. Lithostratigraphic of Inner Mongolia. Wuhan: China University of Geosciences Press, 1996.
|
[17] |
任阿龙. 凉城县岱海地区地热资源预可行性勘查项目取得突出成果[J]. 西部资源, 2014(3):32.
|
[17] |
Ren A L. Outstanding achievements in the pre Feasibility Exploration Project of geothermal resources in Daihai area, Liangcheng County[J]. Westen Resource, 2014(3):32.
|
|
|
|