地球物理学进展 ›› 2019, Vol. 34 ›› Issue (2): 568-572.doi: 10.6038/pg2019CC0332

• 应用地球物理学Ⅰ(油气及金属矿产地球物理勘探) • 上一篇    下一篇

聚类分析算法在大地电磁三维解释中的应用

黄颖,王雪秋(),周子坤,邹宗霖,夏广沛,翁爱华   

  1. 吉林大学地球探测科学与技术学院,长春 130026
  • 收稿日期:2018-07-28 修回日期:2019-01-22 出版日期:2019-04-20 发布日期:2019-05-13
  • 通讯作者: 王雪秋 E-mail:wangxueqiu@jlu.edu.cn
  • 作者简介:黄颖,男,1994年生,硕士研究生,从事大地电磁反演与解释研究.(E-mail:huangying16@mails.jlu.edu.cn)
  • 基金资助:
    国家重大科研仪器专项资助(2011YQ05006010)

Application of clustering analysis algorithm in three-dimensional magnetotelluricinterpretation

HUANG Ying,WANG Xue-qiu(),ZHOU Zi-kun,ZOU Zong-lin,XIA Guang-pei,WENG Ai-hua   

  1. College of Geo-Exploration Sciences and Technology, Jilin University, Changchun 130026, China
  • Received:2018-07-28 Revised:2019-01-22 Online:2019-04-20 Published:2019-05-13
  • Contact: Xue-qiu WANG E-mail:wangxueqiu@jlu.edu.cn

摘要:

为了提高对大地电磁三维反演结果的分析和解释能力,完成构造界面识别、异常构造刻画等地质和地球物理解释,本文提出使用非监督的聚类方法分析大地电磁反演结果.根据三维反演模型电阻率值的分布和各种聚类分析方法的特点,选择使用K均值聚类方法对反演模型电阻率值进行聚类分析.在K均值聚类分析过程中,本文采用了RS指数指导选择聚类数目, Kaufman法进行中心初始化.通过对东北地区大地电磁数据三维反演结果使用K均值聚类分析方法,得到了东北地区电性岩石圈的厚度估计,结果表明东北地区岩石圈底部聚类电阻率大约为339 Ω·m,其中松辽盆地岩石圈最薄,约为60 km;大兴安岭地区最厚,约为150 km;佳木斯地体厚度约为100 km;而长白山地区岩石圈厚度不易确定,可能受新生代构造活动影响,电阻率明显减小.聚类方法能够有效地帮助对大地电磁三维反演结果中的地质构造进行识别和归类.

关键词: 大地电磁测深, 聚类分析, 中国东北地区, K均值聚类, 岩石圈厚度

Abstract:

In order to improve the ability of analyzing and interpreting the results of magnetotelluric three-dimensional inversion and complete geological and geophysical interpretation such as structural interface recognition and anomalous structures characterization, this paper proposes an unsupervised clustering method to analyze the results of magnetotelluric 3D inversion. According to the distribution of resistivity values in magnetotelluric 3D inversion model and the characteristics of various clustering analysis methods, the K-means clustering method is selected to analyze the resistivity values in inversion results. In the process of the K-means clustering method, this paper uses RS index to guide the selection of clustering number and Kaufman method to initialize the clustering centers. By using the K-means clustering methods to analyze the 3D inversion results of magnetotelluric data in Northeast China, the thickness estimation of electrical lithosphere in Northeast China is obtained. The result shows that the cluster resistivity at the bottom of the lithosphere in Northeast China is about 339 Ω·m. The lithosphere of Songliao Basin is the thinnest, about 60 km; the lithosphere of Great Xing’an Range is the thickest, about 150 km; Jiamusi Block is about 100 km; and the lithosphere thickness of Changbai Mountain Range is not easy to determine, which may be influenced by Cenozoic tectonic activities, and its resistivity is obviously reduced. The result shows that the clustering methods can effectively help to identify and classify the geotectonic structures in the results of magnetotelluric 3D inversion.

Key words: Magnetotelluric, Cluster analysis, Northeast China, K-means, Lithospheric thickness

中图分类号: