大地电磁三维反演云计算系统toPeak的设计与实现

doi:10.3969/j.issn.0253-4967.2022.03.015

地震地质 ›› 2022, Vol. 44 ›› Issue (3): 802-820.DOI: 10.3969/j.issn.0253-4967.2022.03.015

• 极低频地震电磁专题文章 • 上一篇下一篇

大地电磁三维反演云计算系统toPeak的设计与实现

刘钟尹¹⁾(), 陈小斌¹^,²⁾(), 蔡军涛¹⁾, 崔腾发¹⁾, 赵国泽¹⁾, 汤吉¹⁾, 欧阳飚¹⁾

1)中国地震局地质研究所, 地震动力学国家重点实验室, 北京 100029
2)应急管理部国家自然灾害防治研究院, 北京 100085

收稿日期:2021-02-12 修回日期:2021-07-02 出版日期:2022-06-20 发布日期:2022-08-02
通讯作者: 陈小斌
作者简介:刘钟尹, 男, 1990年生, 现为中国地震局地质研究所固体地球物理专业在读博士研究生, 研究方向主要为电磁测深方法应用和软件开发, E-mail: liuzy@ies.ac.cn。
基金资助:
中国地震局地壳应力研究所基本科研业务专项(ZDJ2019-26);国家重点研发计划项目(2017YFC1500204-1)

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

LIU Zhong-yin¹⁾(), CHEN Xiao-bin¹^,²⁾(), CAI Jun-tao¹⁾, CUI Teng-fa¹⁾, ZHAO Guo-ze¹⁾, TANG Ji¹⁾, OUYANG Biao¹⁾

1) State Key Laboratory of Earthquake Dynamics, Institute of Geology, China Earthquake Administration, Beijing 100029, China
2) National Institute of Natural Hazards, Ministry of Emergency Management of People’s Republic of China, Beijing 100085, China

Received:2021-02-12 Revised:2021-07-02 Online:2022-06-20 Published:2022-08-02
Contact: CHEN Xiao-bin

摘要/Abstract

摘要：

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

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

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

中图分类号:

P315.721

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

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.

图/表 14

图 1 大地电磁三维反演云计算系统的结构

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

图 2 服务端的数据储存结构

Fig. 2 Server data storage structure.

图 3 中间件的封装与多态

Fig. 3 The encapsulation and polymorphism of the middleware.

图 4 客户端主要的设计思想

Fig. 4 Main design idea of client.

图 5 客户端的架构

Fig. 5 Structure of client.

图 6 客户端的数据工程结构

Fig. 6 Client’s data project structure.

图 7 云端配置与应用功能

Fig. 7 Function of cloud configuration and application.

图 8 初始模型构建功能

Fig. 8 Function of initial model building.

图 9 toPeak操作流程示意图

Fig. 9 Workflow of toPeak.

图 10 部分功能的示意图 a 测点分布经纬度显示; b 测点分布相对直角坐标显示; c 初始模型构建; d 地形显示; e 反演数据拟合显示; f 反演拟合残差显示

Fig. 10 Software’s functions.

图 11 部分功能的示意图 a 选择云端计算节点和计算核; b 云端任意数据工程指定反演任务的拟合情况; c 当前数据工程指定反演的实时计算信息

Fig. 11 Software’s function.

图 12 某地震地壳的三维电性结构分布图(崔腾发等,2020)

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

图 13 某油田MT剖面的三维反演结果

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

图 14 某铁路隧道AMT数据的三维反演结果

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

参考文献 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.

大地电磁三维反演云计算系统toPeak的设计与实现

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

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 14

参考文献 28

相关文章 15

编辑推荐

Metrics

本文评价

[1]	李东臣, 任俊杰, 张志文, 刘亮. 基于高分辨率无人机影像的地震地表破裂半自动提取方法--以2021年M_S7.4青海玛多地震为例[J]. 地震地质, 2022, 44(6): 1484-1502.
[2]	赵凌强, 胡亚轩, 王庆良, 祝意青, 操聪, 李仲巍, 綦伟, 汶宇龙. 吉林龙岗火山区深部电性结构特征分析[J]. 地震地质, 2022, 44(4): 845-858.
[3]	董泽义, 汤吉, 赵国泽, 陈小斌, 崔腾发, 韩冰, 姜峰, 王立凤. 首都圈极低频电磁台网区地下电性结构探测[J]. 地震地质, 2022, 44(3): 649-668.
[4]	赵凌强, 詹艳, 王庆良, 孙翔宇, 韩静, 操聪, 张松, 蔡妍. 大地电磁数据揭示的1303年洪洞8级地震区精细结构[J]. 地震地质, 2022, 44(3): 686-700.
[5]	韩静, 詹艳, 孙翔宇, 赵国泽, 刘雪华, 包雨鑫, 孙建宝, 彭远黔. 强电磁干扰环境下的大地电磁数据特征及处理[J]. 地震地质, 2022, 44(3): 736-752.
[6]	张赟昀, 王培杰, 陈小斌, 詹艳, 韩冰, 王立凤, 赵国泽. 强干扰环境下的大地电磁时间序列处理过程[J]. 地震地质, 2022, 44(3): 786-801.
[7]	阮帅, 汤吉, 董泽义, 王立凤, 邓琰, 韩冰. 基于三维大地电磁AR-QN反演的长白山天池火山区电性结构[J]. 地震地质, 2020, 42(6): 1282-1300.
[8]	孙翔宇, 詹艳, 赵国泽, 赵凌强, 邓琰, 胡亚轩, 胡久常, 向小娟. 琼东北马鞍岭-雷虎岭火山区深部岩浆系统大地电磁三维探测[J]. 地震地质, 2020, 42(3): 640-653.
[9]	孙翔宇, 詹艳, 赵凌强, 陈小斌, 李陈侠, 孙建宝, 韩静, 崔腾发. 东昆仑断裂带东端和2017年九寨沟7.0级地震区深部电性结构探测[J]. 地震地质, 2020, 42(1): 182-197.
[10]	姜峰, 陈小斌, 董泽义, 崔腾发, 刘钟尹, 王培杰. 利用单剖面大地电磁三维反演识别沙德和玉农希断裂[J]. 地震地质, 2019, 41(6): 1444-1463.
[11]	张继红, 赵国泽, 董泽义, 王立凤, 韩冰, 王庆林, 唐廷梅, 王梅. 郯庐断裂带安丘、莒县电磁台地壳电性结构研究[J]. 地震地质, 2019, 41(5): 1239-1253.
[12]	赵凌强, 詹艳, 王庆良, 孙翔宇, 杨皓, 陈小斌. 1954年甘肃民勤7级地震区深部电性结构特征及地震构造环境研究[J]. 地震地质, 2018, 40(3): 552-565.
[13]	韩江涛, 王天琪, 刘文玉, 刘国兴, 韩松, 刘立家. 阿尔山火山群深部“拱桥式”岩浆系统及其稳定性分析[J]. 地震地质, 2018, 40(3): 590-610.
[14]	周耀明, 朱文斌, 陈正乐, 朱炳玉, 薛峰. 准噶尔盆地克-百断裂带火山岩分布特征的重磁资料解释[J]. 地震地质, 2018, 40(3): 641-655.
[15]	翁爱华, 李建平, 范小平, 李斯睿, 韩江涛, 李大俊, 李亚彬, 赵祥阳, 唐裕. 大地电磁测深揭示的1668年郯城8.5级地震震中地壳精细结构[J]. 地震地质, 2018, 40(2): 396-409.