One-dimensional Occam's inversion for transient electromagnetic data excited by a loop source
-
摘要: 基于Dipole 1D一维正演程序和Occam算法开发了回线源瞬变电磁法一维反演程序,并以四层介质模型验证了其正确性。对于倾斜界面模型,待反演数据由时域矢量有限单元法三维正演计算,通过对比实际模型和等效模型的反演结果,说明了一维反演可直接应用于此类情形。最后,开展了那仁宝力格煤田玄武岩深部探测案例研究。在反演电阻率断面图中,玄武岩整体形态呈“锅底状,与多个钻孔资料吻合良好,进一步验证了该反演程序的正确性。Abstract: A 1D inversion code is developed for the loop-source transient electromagnetic method (TEM) based on the open-source code Dipole1D and Occam's algorithm. This code is tested by a four-layer stratified model. Then, the model with a tilted earth-air interface is considered, for which the synthetic data are calculated by 3D finite-element method. The inversion results show that 1D inversion can be directly used for the scenario with tilted interface. Finally, this 1D inversion code is used for a field case, in which TEM is employed to delineate the 3D distributed domains of a basalt which intruded into shale and sandstone. The inversion result shows that the thickness of the basalt coincides with the drilling data, and the shape of the basalt like a pot bottom.
-
Key words:
- loop source /
- TEM /
- Occam's method /
- 1D inversion
-
-
[1] Flores C, Peralta-Ortega S A. Induced polarization with in-loop transient electromagnetic soundings: A case study of mineral discrimination at El Arco porphyry copper, Mexico[J]. Journal of Applied Geophysics, 2009, 68(3):423-436.
[2] Xue G Q, Qin K Z, Li X, et al. Discovery of a large-scale porphyry molybdenum deposit in Tibet through a modified TEM exploration method[J]. Journal of Environmental and Engineering Geophysics, 2012, 17(1):19-25.
[3] Yang D, Oldenburg D W. Three-dimensional inversion of airborne time-domain electromagnetic data with applications to a porphyry deposit[J]. Geophysics, 2012, 77(2):B23-B34.
[4] Scholl C, Helwig S L, Tezkan B, et al. 1-D multimodel joint inversion of TEM-data over multidimensional structures[J]. Geophysical Journal International, 2009, 176(1):81-94.
[5] Kirkegaard C, Auken E. A parallel, scalable and memory efficient inversion code for very large-scale airborne electromagnetics surveys[J]. Geophysical Prospecting, 2015, 63(2):495-507.
[6] 孙怀凤, 张诺亚, 柳尚斌, 等. 基于L1范数的瞬变电磁非线性反演[J]. 地球物理学报, 2019, 62(12):4860-4873.
[7] Sun H F, Zhang N Y, Liu S B, et al. L1-norm based nonlinear inversion of transient electromagnetic data[J]. Chinese J. Geophys., 2019, 62(12):4860-4873.
[8] Liu Y, Yin C, Qiu C, et al. 3-D inversion of transient EM data with topography using unstructured tetrahedral grids[J]. Geophysical Journal International, 2019, 217(1):301-318.
[9] Auken E, Christiansen A V, Kirkegaard C, et al. An overview of a highly versatile forward and stable inverse algorithm for airborne, ground-based and borehole electromagnetic and electric data[J]. Exploration Geophysics, 2015, 46(3):223-235.
[10] Raiche A P, Jupp D L B, Ruter H, et al. The joint use of coincident loop transient electromagnetic and Schlumberger sounding to resolve layered structures[J]. Geophysics, 1985, 50(10):1618-1627.
[11] 黄皓平, 王维中. 时间域航空电磁数据的反演[J]. 地球物理学报, 1990, 33(1):87-97.
[12] Huang H P, Wang W Z. Inversion of time-domain airborne electromagnetic data[J]. Chinese J. Geophys., 1990, 33(1):87-97.
[13] Farquharson C G, Oldenburg D W. Inversion of time-domain electromagnetic data for a horizontally layered earth[J]. Geophysical Journal International, 1993, 114(3):433-442.
[14] Effersø F, Auken E, Sørensen K I. Inversion of band-limited TEM responses[J]. Geophysical Prospecting, 1999, 47(4):551-564.
[15] Auken E, Christiansen A V. Layered and laterally constrained 2D inversion of resistivity data[J]. Geophysics, 2004, 69(3):752-761.
[16] Auken E, Christiansen A V, Jacobsen L H, et al. A resolution study of buried valleys using laterally constrained inversion of TEM data[J]. Journal of Applied Geophysics, 2008, 65(1):10-20.
[17] Viezzoli V, Christiansen A V, Auken E, et al. Quasi-3D modeling of airborne TEM data by spatially constrained inversion[J]. Geophysics, 2008, 73(3):F105-F113.
[18] Constable S C, Parker R L, Constable C G. Occam’s inversion: A practical algorithm for generating smooth models from electromagnetic sounding data[J]. Geophysics, 1987, 52(3):289-300.
[19] Sudha , Tezkan B, Siemon B. Appraisal of a new 1D weighted joint inversion of ground based and helicopter-borne electromagnetic data[J]. Geophysical Prospecting, 2014, 62(3):597-614.
[20] 毛立峰, 王绪本, 陈斌. 直升机航空瞬变电磁自适应正则化一维反演方法研究[J]. 地球物理学进展, 2011, 26(1):300-305.
[21] 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.
[22] 齐彦福, 殷长春, 王若, 等. 多通道瞬变电磁m序列全时正演模拟与反演[J]. 地球物理学报, 2015, 58(7):2566-2577.
[23] Qi Y F, Yin C C, Wang R, et al. Multi-transient EM full-time forward modeling and inversion of m-sequences[J]. Chinese J. Geophys., 2015, 58(7):2566-2577.
[24] 李海, 薛国强, 钟华森, 等. 多道瞬变电磁法共中心点道集数据联合反演[J]. 地球物理学报, 2016, 59(12):4439-4447.
[25] Li H, Xue G Q, Zhong H S, et al. Joint inversion of CMP gathers of multi-channel transient electromagnetic data[J]. Chinese J. Geophys., 2016, 59(12):4439-4447.
[26] Li M, Cheng J, Wang P, et al. Transient electromagnetic inversion based on the PSO-DLS combination algorithm[J]. Exploration Geophysics, 2019, 50(5):472-480.
[27] Key K. 1D inversion of multicomponent, multifrequency marine CSEM data: Methodology and synthetic studies for resolving thin resistive layers[J]. Geophysics, 2009, 74(2):F9-F20.
[28] Guptasarma D, Singh B. New digital linear filters for Hankel J0 and J1 transforms[J]. Geophysical Prospecting, 1997, 45(5):745-762.
[29] Anderson W L. Fourier cosine and sine transforms using lagged convolutions in double-precision (subprograms DLAGF0/DLAGF1) [R]. U.S. Geological Survey, 1983.
[30] Li J, Lu X, Farquharson C G, et al. A finite-element time-domain forward solver for electromagnetic methods with complex-shaped loop sources[J]. Geophysics, 2018, 83(3):E117-E132.
-
计量
- 文章访问数: 658
- PDF下载数: 22
- 施引文献: 0
下载: