Visualization design of 1D CSAMT forward modeling and research on the induced polarization effect
-
摘要: 针对可控源音频大地电磁法(CSAMT)开发了一套CSAMT一维数据正演可视化软件。其界面友好,操作简单;使用MATLAB进行核心算法的编写,正演计算过程引入线性滤波系数对汉克尔积分进行求解,并对不同线性滤波系数方案的精度进行了比较;基于Java语言进行操作界面设计,提供多参数设置界面,自动化程度高;通过软件自带的绘图功能,将数据处理结果直观地提供给用户。在计算过程中引入Cole-Cole模型参数,以复电阻率代替不考虑地电体极化效应的直流电阻率,对CSAMT场源一维层状激化介质模型进行了正演模拟,讨论了极化层分别位于不同深度时一维层状介质的极化效应。研究发现极化层埋藏越浅,极化效应对于CSAMT的影响频带越宽。这对于认识含极化效应的CSAMT方法的电磁响应特征,提高数据处理便利度与工作效率方面都具有重要的作用。Abstract: This study developed a piece of visualization software of the 1D data forward modeling applicable to the controlled source audio-frequency magnetotellurics (CSAMT). This software has the advantages of friendly interfaces and simple operations. In detail, MATLAB was used to develop the core algorithms, and linear filter coefficients were introduced to solve Hankel integral in the process of forward modeling. Meanwhile, the accuracy of schemes using different linear filter coefficients was compared. The operation interfaces were designed using Java, and multi-parameter setting interfaces with a high automation degree are available. Moreover, the data processing results can be intuitively provided to users using the drawing function of the software. In the calculation process, Cole-Cole model parameters were introduced, and complex resistivity was used to replace the DC resistivity without considering the polarization effect of geoelectric bodies. Furthermore, this study carried out the forward simulation of CSAMT field source on a one-dimensional layered model of induced polarization (IP) media and discussed the IP effect of one-dimensional layered media with polarization layers at different burial depths. It was found that the frequency band of CSAMT affected by the IP effect widened with a decrease in the burial depth of polarization layers. This result is crucial to understanding the electromagnetic response characteristics of the CSAMT method with IP effect and improving the convenience and efficiency of data processing.
-
Key words:
- 1D forward modeling /
- Hankel transform /
- visualization /
- induced polarization /
- Java
-
-
[1] 龚强, 胡祥云, 孟永良. 基于快速汉克尔变换算法的CSAMT一维正演[J]. 煤田地质与勘探, 2008, 36(1):71-73.
[2] Gong Q, Hu X Y, Meng Y L. Linear dimension forward solution of CSAMT by Hankel transformation method[J]. Coal Geology & Exploration, 2008, 36(1):71-73.
[3] 汤井田, 何继善. 可控源音频大地电磁法及其应用[M]. 长沙: 中南大学出版社, 2005.
[4] Tang J T, He J S. Theory and application of CSAMT method[M]. Changsha: Central South University, 2005.
[5] 王佳宝. CSAMT二维反演程序设计与实现[D]. 北京: 中国地质大学(北京), 2016.
[6] Wang J B. Design and implementation of CSAMT two-dimensional inversion program[D]. Beijing: China University of Geosciences, 2016.
[7] 张健. 可控源音频大地电磁法一维正反演研究[D]. 成都: 成都理工大学, 2011.
[8] Zhang J. Research on one-dimensional forward and inversion of CSAMT method[D]. Chengdu: Chengdu University of Technology, 2011.
[9] 王珺璐, 刘明文, 李荡, 等. 基于Cole-Cole复电阻率模型的线源可控源有限元数值模拟[J]. 物探化探计算技术, 2015, 37(1):32-39.
[10] Wang J L, Liu M W, Li D, et al. Finite element numerical simulation on line controlled source based on Cole-Cole model[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2015, 37(1):32-39.
[11] 唐睿, 唐新功, 于鹏. CSAMT正演模拟中激电效应的初步研究[C]// 第十届中国国际地球电磁学术讨论会论文集, 2011.
[12] Tang R, Tang X G, Yu P. The preliminary research on IP effect in CSAMT[C]// Proceedings of The 10th China International Geo-Electromagnetic Workshop, 2011.
[13] 岳安平, 底青云, 王妙月, 等. 油气藏MT激电效应一维正演研究[J]. 石油地球物理勘探, 2009, 44(3):364-370.
[14] Yue A P, Di Q Y, Wang M Y, et al. 1D forward simulation of MT IP effect of reservoir[J]. Oil Geophysical Prospecting, 2009, 44(3):364-370.
[15] 徐玉聪. 含激电效应的CSAMT一维正反演研究[D]. 成都: 成都理工大学, 2016.
[16] Xu Y C. Research on IP effect of one-dimensional forward and inversion of CSAMT[D]. Chengdu: Chengdu University of Technology, 2016.
[17] 侯海涛. Cole-Cole模型等效模拟瞬变电磁响应[D]. 长沙: 中南大学, 2010.
[18] Hou H T. Transient electromagnetic simulation response of Cole-Cole model equivalent[D]. Changsha: Central South University, 2010.
[19] 耿祥义, 张跃平. Java面向对象程序设计[M]. 北京: 清华大学出版社, 2010.
[20] Geng X Y, Zhang Y P. Java object-oriented programming[M]. Beijing: Tsinghua University, 2010.
[21] 余传涛. 含激电效应的CSAMT二维正演与联合反演应用研究[D]. 太原: 太原理工大学, 2012.
[22] Yu C T. 2D Forward modeling and joint inversion of the CSAMT signal incorporation IP effect[D]. Taiyuan: Taiyuan University of Technology, 2012.
[23] 蔡盛. 快速汉克尔变换及其在电法正演计算中的应用研究[D]. 长沙: 中南大学, 2013.
[24] Cai S. The fast hankel transformation and its applications in electrical exploration forward calculations[D]. Changsha: Central South University, 2013.
[25] 蔡盛. 快速汉克尔变换及其在正演计算中的应用[J]. 地球物理学进展, 2014, 29(3):1384-1390.
[26] Cai S. Fast Hankel transformation and its application in forward computing[J]. Progress in Geophysics, 2014, 29(3):1384-1390.
[27] 胡俊, 聂在平. 索末菲尔德积分新方法——快速汉克尔变换[J]. 电子学报, 1998(3):126-128.
[28] Hu J, Nie Z P. A new method to calculate Sommerfeld integral:Fast hankel transform[J]. Acta Electronica Sinica, 1998(3):126-128.
[29] 底青云, 王若. 可控源音频大地电磁数据正反演及方法应用[M]. 北京: 科学出版社, 2008.
[30] Di Q Y, Wang R. Forward and inversion of CSAMT data and its application[M]. Beijing: Science Press, 2008.
-
计量
- 文章访问数: 518
- PDF下载数: 33
- 施引文献: 0
下载: