THE DESIGN AND APPLICATION OF TOPEAK: A THREE-DIMENSIONAL MAGNETOTELLURIC INVERSION CLOUD COMPUTING SYSTEM

doi:10.3969/j.issn.0253-4967.2022.03.015

Abstract

Abstract:

Magnetotelluric(MT)three-dimensional inversion has the advantages of simple data preprocessing, the model is close to actual situation, and the inversion result is more reliable and stable. It is one of the most advanced research topics and would take the place of the dominant two-dimensional inversion definitely. With the improvement of computing capability of computers and the breakthrough in inversion methods, great progress was made in MT three-dimensional inversion in recent years, from the theoretical research and test of this method at the beginning to the current application to practical data interpretation. For the great computation amount of MT three-dimensional inversion, current MT three-dimensional inversion algorithm programs are all implemented in parallel way and it is recommended to do three-dimensional inversion calculations on supercomputing system to make better use of computing resources and improve the inversion efficiency.

Different from the MT three-dimensional inversion algorithm programs which have basically realized the utility function, the practical application of MT three-dimensional inversion is still in an early stage. Users should be familiar with the use of multiple software and fulfill the function manually with the help of the software as follows: generating the files required for the inversion program, connecting to the supercomputer to upload data, inputting the command to perform the inversion, etc. The process of manually connecting and operating calculations is the most primitive cloud computing. All processes need to be done manually, which would cause not only heavy workload and the complicated operation, but also the problems for the long-term effective preservation and management of complex inversion data.

To conquer this, we develop independently a three-dimensional magnetotelluric inversion cloud computing system, toPeak, using Delphi language. This paper introduces some main features of toPeak. To begin with, system design and analysis are carried out in combination with the current situation and system structure and functions are realized. The main idea is to realize a set of cloud computing system platform based on server-client(C/S), on the basis of perfect inversion data management, integrate the most advanced three-dimensional magnetotelluric inversion algorithm program in the cloud, and connect through the Internet to realize all the system functions of three-dimensional magnetotelluric inversion. Then, the different parts of toPeak are introduced separately, including design structures and designs. The server is deployed on the supercomputer system(supercomputing)to receive the data for inversion tasks, configure and manage the storage of the inversion result data. Combined with the Internet connection, the server and the Internet together constitute a computing cloud. The client is deployed on the users’ windows operating system, including Windows visual data integration processing software and Internet operation middleware. The client is designed on the basis of object-oriented programming ideas, with data as the core, using data engineering objects to encapsulate and store all MT data, process and interpret the results, realize data processing inversion and other operations around this data project, and display the process and results of these processing and inversion in graphics using visualization technology. Internet operation middleware connects the client and server based on the SSH protocol to realize data processing and inversion, transmission and command sending and receiving. Furthermore, the whole work flow of inversion using toPeak and parts of procedure of it are shown. At last, some inversion results from toPeak are displayed. toPeak has realized the full functions require for implementing three-dimensional inversion and can grid, process and select, inverse and explain the data. It is a good tool for the practical use of three-dimensional inversion.

Key words: magnetotelluric, object oriented, three-dimensional inversion, visualization, cloud computing, supercomputer

摘要：

三维反演是大地电磁领域研究的热点问题, 已有的三维反演算法程序使用繁琐、操作不便, 导致三维MT反演的实用化进度缓慢。为此, 我们采用Delphi语言自主研发了大地电磁三维反演云计算系统toPeak。文中对三维MT反演云计算系统进行了需求分析和系统设计, 确定了系统结构和功能模块; 阐述了toPeak各部分的设计结构和主要功能; 介绍了云反演的流程, 并展示了使用toPeak进行三维云反演的案例。 toPeak已实现较大规模的大地电磁三维反演所需的全部功能, 可对较大规模的大地电磁三维实测资料进行网格化、流程化资料处理、反演与解释, 为大地电磁三维反演实用化提供了一种有力工具。

关键词: 大地电磁, 面向对象, 三维反演, 可视化, 云计算, 超算

CLC Number:

P315.721

LIU Zhong-yin, CHEN Xiao-bin, CAI Jun-tao, CUI Teng-fa, ZHAO Guo-ze, TANG Ji, OUYANG Biao. THE DESIGN AND APPLICATION OF TOPEAK: A THREE-DIMENSIONAL MAGNETOTELLURIC INVERSION CLOUD COMPUTING SYSTEM[J]. SEISMOLOGY AND GEOLOGY, 2022, 44(3): 802-820.

刘钟尹, 陈小斌, 蔡军涛, 崔腾发, 赵国泽, 汤吉, 欧阳飚. 大地电磁三维反演云计算系统toPeak的设计与实现[J]. 地震地质, 2022, 44(3): 802-820.

Figures/Tables 14

Fig. 1 Structure of the three-dimensional magnetotelluric inversion cloud computing system.

Fig. 2 Server data storage structure.

Fig. 3 The encapsulation and polymorphism of the middleware.

Fig. 4 Main design idea of client.

Fig. 5 Structure of client.

Fig. 6 Client’s data project structure.

Fig. 7 Function of cloud configuration and application.

Fig. 8 Function of initial model building.

Fig. 9 Workflow of toPeak.

Fig. 10 Software’s functions.

Fig. 11 Software’s function.

Fig. 12 Three-dimensional electrical structure distribution map of crust of an earthquake(CUI Teng-fa et al., 2020).

Fig. 13 Three-dimensional inversion result of MT profile of an oilfield.

Fig. 14 Three-dimensional inversion result of AMT data of a railway tunnel.

References 28

[1]	蔡军涛, 陈小斌. 2010a. 大地电磁资料精细处理和二维反演解释技术研究(二): 反演数据极化模式选择[J]. 地球物理学报, 53(11): 2703—2714.
	CAI Jun-tao, CHEN Xiao-bin. 2010a. Refined techniques for data processing and two-dimensional inversion in magnetotelluric(Ⅱ): Which data polarization mode should be used in 2D inversion[J]. Chinese Journal of Geophysics, 53(11): 2703—2714 (in Chinese).
[2]	蔡军涛, 陈小斌, 董泽义, 等. 2020. 1976年唐山 M_S7.8 地震震源区三维深部电性结构和发震机制研究 [G]//2020年中国地球科学联合学术年会论文集(十五):电磁地球物理学研究应用及其新进展.
	CAI Jun-tao, CHEN Xiao-bin, DONG Ze-yi, et al. 2020. Study on 3-D deep electrical structure and seismogenic mechanism of Tangshan M_S7.8 earthquake in 1976[G]//Proceedings of Annual Meeting of Chinese Geoscience Union(15): Application and New Progress of Electromagnetic Geophysics (in Chinese).
[3]	蔡军涛, 陈小斌, 赵国泽. 2010b. 大地电磁资料精细处理和二维反演解释技术研究(一): 阻抗张量分解与构造维性分析[J]. 地球物理学报, 53(10): 2516—2526.
	CAI Jun-tao, CHEN Xiao-bin, ZHAO Guo-ze. 2010b. Refined techniques for data processing and two-dimensional inversion in magnetotelluric(Ⅰ): Tensor decomposition and dimensionality analysis[J]. Chinese Journal of Geophysics, 53(10): 2516—2526 (in Chinese).
[4]	陈小斌, 蔡军涛, 王立凤, 等. 2014. 大地电磁资料精细处理和二维反演解释技术研究(四): 阻抗张量分解的多测点-多频点统计成像分析[J]. 地球物理学报, 57(6): 1946—1957.
	CHEN Xiao-bin, CAI Jun-tao, WANG Li-feng, et al. 2014. Refined techniques for magnetotelluric data processing and two-dimensional inversion(4): Statistical image method based on multi-site, multi-frequency tensor decomposition[J]. Chinese Journal of Geophysics, 57(6): 1946—1957 (in Chinese).
[5]	陈小斌, 崔腾发, 蔡军涛, 等. 2017a. 川滇东部三维电性结构及动力学 [G]//2017年中国地球科学联合学术年会论文集(六):陆陆碰撞带深部结构和动力学意义.
	CHEN Xiao-bin, CUI Teng-fa, CAI Jun-tao, et al. 2017a. Zoning characteristics and dynamic interpretation of lithospheric electrical structure in Ordos block [G]// Proceedings of the Annual Meeting of Chinese Geoscience Union, 2017, (6): The Deep Structure and Dynamic Significance of the Continental Collision Zone (in Chinese).
[6]	陈小斌, 郭春玲. 2017b. 大地电磁资料精细处理和二维反演解释技术研究(五): 利用阻抗张量成像识别大地线性构造[J]. 地球物理学报, 60(2): 766—777.
	CHEN Xiao-bin, GUO Chun-ling. 2017b. Refined techniques for data processing and two-dimensional inversion in magnetotelluric(Ⅴ): Detecting the linear structures of the Earth by impedance tensor imaging[J]. Chinese Journal of Geophysics, 60(2): 766—777 (in Chinese).
[7]	陈小斌, 叶涛, 蔡军涛, 等. 2019. 大地电磁资料精细处理和二维反演解释技术研究(七): 云南盈江—龙陵地震区深部电性结构及孕震环境[J]. 地球物理学报, 62(4): 1377—1393.
	CHEN Xiao-bin, YE Tao, CAI Jun-tao, et al. 2019. Refined techniques for data processing and two-dimensional inversion in magnetotelluric(Ⅶ): Electrical structure and seismogenic environment of Yingjiang-Longling seismic area[J]. Chinese Journal of Geophysics, 62(4): 1377—1393 (in Chinese).
[8]	陈小斌, 赵国泽, 詹艳. 2004. MT资料处理与解释的Windows可视化集成系统[J]. 石油地球物理勘探, 39(S1): 11—16, 169.
	CHEN Xiao-bin, ZHAO Guo-ze, ZHAN Yan. 2004. A windows visual integrated system for MT data processing and interpretation[J]. Oil Geophysical Prospecting, 39(S1): 11—16, 169 (in Chinese).
[9]	崔腾发. 2019. 郯庐—大别构造区地壳三维电性结构及其动力学意义[D]. 北京: 中国地震局地质研究所.
	CUI Teng-fa. 2019. Three-dimensional electrical structure of crust beneath the Tanlu-Dabie tectonic region and its dynamical significance[D]. Institute of Geology, China Earthquake Administration, Beijing (in Chinese).
[10]	崔腾发, 陈小斌, 詹艳, 等. 2020. 安徽霍山地震区深部电性结构和发震构造特征[J]. 地球物理学报, 63(1): 256—269.
	CUI Teng-fa, CHEN Xiao-bin, ZHAN Yan, et al. 2020. Characteristics of deep electrical structure and seismogenic structure beneath Anhui Huoshan earthquake area[J]. Chinese Journal of Geophysics, 63(1): 256—269 (in Chinese).
[11]	邓琰. 2019. 玉树地震区三维电性结构及孕震环境研究[D]. 北京: 中国地震局地质研究所.
	DENG Yan. 2019. The study of three-dimensional electrical resistivity structure and the seismogenesis in Yushu earthquake area[D]. Institute of Geology, China Earthquake Administration, Beijing (in Chinese).
[12]	段新华, 王宏勇, 丁汨. 2009. 集群作业管理系统的关键技术分析及比较[J]. 计算机技术与发展, 19(2): 87—90.
	DUAN Xin-hua, WANG Hong-yong, DING Mi. 2009. A core technology analysis and comparison for cluster job management system[J]. Computer Technology and Development, 19(2): 87—90 (in Chinese).
[13]	郭春玲, 陈小斌. 2018. 大地电磁资料精细处理和二维反演解释技术研究(六): 交错模型的大地电磁二维反演[J]. 地球物理学报, 61(6): 2548—2559.
	GUO Chun-ling, CHEN Xiao-bin. 2018. Refined processing and two-dimensional inversion of magnetotelluric(MT)data(Ⅵ): Two-dimensional magnetotelluric inversion based on the staggered model[J]. Chinese Journal of Geophysics, 61(6): 2548—2559 (in Chinese).
[14]	叶涛, 陈小斌, 严良俊. 2013. 大地电磁资料精细处理和二维反演解释技术研究(三): 构建二维反演初始模型的印模法[J]. 地球物理学报, 56(10): 3596—3606.
	YE Tao, CHEN Xiao-bin, YAN Liang-jun. 2013. Refined techniques for data processing and two-dimensional inversion in magnetotelluric(Ⅲ): Using the impressing method to construct starting model of 2D magnetotelluric inversion[J]. Chinese Journal of Geophysics, 56(10): 3596—3606 (in Chinese).
[15]	赵凌强, 詹艳, 王庆良, 等. 2020. 祁连山东端冷龙岭隆起及邻区深部电性结构与孕震构造背景[J]. 地球物理学报, 63(3): 1014—1025.
	ZHAO Ling-qiang, ZHAN Yan, WANG Qing-liang, et al. 2020. The deep electrical structure and seismogenic background of Lenglongling uplift and its adjacent areas in the eastern end of Qilian Mountains[J]. Chinese Journal of Geophysics, 63(3): 1014—1025 (in Chinese).
[16]	Compagnie Générale de Géophysique. 2013. GeoTools[EB/OL].(2019-01-28)[2020-06-22]. https://www.cgg.com/en/What-We-Do/Multi-Physics/Imaging/Software/Geotools-Magnetotellurics.
[17]	Gamma E, Helm R, Johnson R, et al. 1994. Design Patterns: Elements of Reusable Object-Oriented Software[M]. Reading, Massachusetts, Addison-Wesley Professional.
[18]	Fu H, Liao J, Yang J, et al. 2016. The Sunway TaihuLight supercomputer: System and applications[J]. Science China Information Sciences, 59(7): 1—16.
[19]	Henderson R L. 1995. Job scheduling under the Portable Batch System [G]. Workshop on Job Scheduling Strategies for Parallel Processing. Springer Berlin Heidelberg, Berlin: 279—294.
[20]	Kelbert A, Egbert G D, Schultz A. 2008. Non-linear conjugate gradient inversion for global EM induction: Resolution studies[J]. Geophysical Journal International, 173(2): 365—381. DOI URL
[21]	Kelbert A, Meqbel N, Egbert G D, et al. 2014. ModEM: A modular system for inversion of electromagnetic geophysical data[J]. Computers & Geosciences, 66: 40—53. DOI URL
[22]	Madden T M, Mackie R L. 1989. Three-dimensional magnetotelluric modelling and inversion[J]. Proceedings of the IEEE, 77(2): 318—333. DOI URL
[23]	Newman G A, Alumbaugh D L. 2000. Three-dimensional magnetotelluric inversion using non-linear conjugate gradients[J]. Geophysical Journal International, 140(2): 410—424. DOI URL
[24]	Schlumberger. 2013. WinGLink geophysical interpretation software[EB/OL].(2019-07-29)[2020-06-22]. https://www.slb.com/services/seismic/seismic-reservoir-characterization/electromagnetics/emsoftware/winglink.aspx.
[25]	Siripunvaraporn W, Egbert G. 2009. WSINV3DMT: Vertical magnetic field transfer function inversion and parallel implementation[J]. Physics of the Earth and Planetary Interiors, 173(3-4): 317—329.
[26]	Siripunvaraporn W, Egbert G, Lenbury Y, et al. 2005. Three-dimensional magnetotelluric inversion: Data-space method[J]. Physics of the Earth and Planetary Interiors, 150(1-3): 3—14. DOI URL
[27]	Turkoglu E. 2009. A magnetotelluric investigation of the Arabia-Eurasia collision in eastern Anatolia[D]. University of Alberta, Edmonton.
[28]	Yoo A B, Jette M A, Grondona M. 2003. SLURM: Simple Linux utility for resource management [G]. Workshop on Job Scheduling Strategies for Parallel Processing. Springer Berlin Heidelberg, Berlin: 44—60.