基于机器学习的长宁地震三维断层面几何特征建模

龚丽文, 张怀, 陈石, David A. Yuen, 陈丽娟, Brennan Brunsvik, 尹光耀

龚丽文,张怀,陈石,David A. Yuen,陈丽娟,Brennan Brunsvik,尹光耀. 2023. 基于机器学习的长宁地震三维断层面几何特征建模. 地震学报,45(6):1040−1054. DOI: 10.11939/jass.20220079
引用本文: 龚丽文,张怀,陈石,David A. Yuen,陈丽娟,Brennan Brunsvik,尹光耀. 2023. 基于机器学习的长宁地震三维断层面几何特征建模. 地震学报,45(6):1040−1054. DOI: 10.11939/jass.20220079
Gong L W,Zhang H,Chen S,David A. Yuen,Chen L J,Brennan B,Yin G Y. 2023. Geometry features modeling of three-dimensional fault plane of Changning earthquake based on machine learning. Acta Seismologica Sinica45(6):1040−1054. DOI: 10.11939/jass.20220079
Citation: Gong L W,Zhang H,Chen S,David A. Yuen,Chen L J,Brennan B,Yin G Y. 2023. Geometry features modeling of three-dimensional fault plane of Changning earthquake based on machine learning. Acta Seismologica Sinica45(6):1040−1054. DOI: 10.11939/jass.20220079

基于机器学习的长宁地震三维断层面几何特征建模

基金项目: 国家杰出青年科学基金项目(41725017)、中国地震局地球物理研究所基本科研业务费专项(DQJB19A0121,DQJB21R30)和国家重点研发计划课题项目(2020YFA0713401)共同资助
详细信息
    作者简介:

    龚丽文,在读博士研究生,工程师,主要从事计算地球动力学与地震预报研究,e-mail: gongliwen21@mails.ucas.cn

    通讯作者:

    陈石,博士,研究员,主要从事时变重力位场数据处理与反演解释研究,e-mail:chenshi@cea-igp.ac.cn

  • 中图分类号: P315. 63

Geometry features modeling of three-dimensional fault plane of Changning earthquake based on machine learning

  • 摘要:

    结合长宁地区大量的地震精定位数据和其它研究成果,利用监督分类和聚类分析等机器学习算法,基于地震簇的形态特征和地震震源机制解,编写了一套自动化提取三维破裂面形态特征的程序,获取了长宁地区地震破裂面的精细结构,可为相关研究提供可参考的发震构造模型。结果显示通过聚类分析最终获取了四个地震簇,结合对应的震源机制解节面信息,最终拟合出四条破裂面,其中:长宁背斜上的破裂面沿狮子滩背斜下部的高速体呈NW−SE方向展布,破裂面平直,倾角较陡,倾向SE;建武向斜内部的三条破裂面,主要分布在向斜两翼,规模较小,走向分别为NW,NNE和NNW,从外部包围了建武向斜核部的高速体,破裂面的展布方向与该地区三个主要震源机制解节面的产状一致,其中新城镇附近的NNW向破裂面切割深度较深,约为20 km,且倾向ENE,倾角约为70°。此外,结合地质构造背景和速度结构等反演结果推断,地震破裂面主要存在于先期形成的构造薄弱带或断裂带,例如背斜的核部和向斜的翼部因节理面贯通所形成的薄弱带以及高速体周围的软弱带,在构造应力的加载和工业开采下更容易微破裂成核,形成典型的发震构造。

    Abstract:

    In recent years, the seismicity of the Changning area in the Sichuan Province has increased significantly. Seismogenic models and seismogenic structures on the background of structural loading coupled with human activities have gradually become the focus of research in the field. Using abundant and accurate hypocenters in the Changning area, we established a program for automatically extracting morphological fault features by using machine learning algorithms including supervised classification and clustering. The method provides a reliable, detailed model of seismogenic faults for relative researches. As a result, four earthquake clusters were identified by clustering analysis, and four fracture planes were fit based on the distribution of hypocenters. The fracture plane on the Changning anticline spreads in NW-SE direction along a high velocity body beneath the Shizitan anticline. The fracture plane is straight with steep dip angle, and inclines SE. The three fracture planes in the inner part of the Jianwu syncline are mainly distributed in small scale on the limbs of the syncline with strike of NW, NNE, and NNW, respectively. They are also distributed in the periphery of the high-velocity body at the core of the Jianwu syncline, and their spreading directions are consistent with the strike of nodal planes of three main focal mechanism solutions in this area. Among these fracture planes, the Xincheng fracture plane extends deep to about 20 km and dips ENE with dip angle 70°. Based also on the geological tectonic settings and velocity structure, the fracture planes mainly exist in weak tectonic zones, such as the nucleus of the anticline and the limbs of the syncline. In particular, the fragile zone around the high-velocity body is more likely to rupture and nucleate under the loading of tectonic stress and industrial mining, forming new seismogenic structures.

  • 地震观测台网记录到的事件中包含天然地震事件和非天然地震事件(如爆破、塌陷、矿震、核试验等),二者均是地震学的研究对象。天然地震是地壳应力积累和释放的结果,携带了地球深部构造、介质和应力的信息。天然地震的活动性是地震危险性判定的重要依据,但很多情况下天然地震不容易辨识,特别是天然地震与非天然地震事件的区分,一直是监测工作中的重要问题(刘莎等,2012)。天然地震与爆破的波形具有相似性,而在震源性质(曾融生等,2000)、频谱特征等方面有一定的差异,监测实践中,事件类型的判别主要依赖监测工作人员的经验(张帆等,2016),因而受到主观性的影响,无统一标准,也会出现误判的情形。随着地震台网规模的增加,地震台网监测能力的增强,地震记录数据量也大幅增加,人工处理方式在处理效率上面临挑战,因此产生了自动化处理和快速识别的需求。

    通过不同事件类型在时域、频域等方面的差异(杨选辉等,2005刘莎等,2012王婷婷等,2013),可以设计基于特征和判别规则的自动分类方法,这些方法的判据有较明确的物理意义,且规则相对简单。随着计算机技术的发展,机器学习算法也被应用在地震事件分类中,以数据驱动的方式训练模型参数,获得事件分类模型,如支持向量机、神经网络等(边银菊,2002黄汉明等,2010毕明霞等,2011张博等,2014范晓易等,2019蔡杏辉等,2020)。传统机器学习应用的主要困难在于特征工程中有大量需要人工完成的工作,且特征提取过程是基于专家知识的降维过程,可能丢失原始数据蕴含的有价值信息。

    深度学习是机器学习和人工智能研究的最新发展方向之一,在图像和语音识别、自然语言处理等方面已取得较大进步,也被推广应用到各个领域并取得成效。深度学习通过与人脑相似的自下而上的层级结构,将输入端的数据逐层次提取(余凯等,2013),建立了底层数据到高层语义结构的映射。相较传统机器学习,深度学习具有以下主要特点:在输入数据方面,深度学习对原始数据人工加工较少,保留了大部分原始信息;在模型结构方面,主要使用深层神经网络,具有较复杂的模型结构和较多的模型参数;在流程方面,深度学习往往采用端到端的解决方案,而传统机器学习将任务分割为几个模块,分别训练。深度学习在地震事件分类问题中获得了较好的效果(陈润航等,2018隗永刚等,2019)。

    内蒙古地区有丰富的人工爆破事件记录,本文拟基于内蒙古地区天然地震和爆破的波形资料构建试验数据集,通过短时傅里叶变换将事件波形转换为时频域的振幅谱图,作为卷积神经网络的输入,使用计算机视觉领域的卷积神经网络作为分类器实现地震事件分类,并对结果进行测试,以评估该算法的有效性和稳定性。

    本文采用的原始数据为2016年3月至2020年5月内蒙古地震台网记录的事件波形,由于其中天然地震事件较爆破事件多,为了使数据集平衡,仅选取ML2.7以上的天然地震事件。最终采用的数据包含417个爆破事件的900条台站记录和519个地震事件的900条记录,共计1 800条单台事件记录,选取垂直向记录作为原始数据。图1给出了事件空间分布图,天然地震事件在空间上分布较均匀,爆破的空间分布相对集中。图2为爆破事件的震级频次统计,爆破事件震级分布在ML1.5—3.5之间;天然地震的震级分布在ML2.7—4.5之间。

    图  1  地震和爆破事件的空间分布图
    Figure  1.  Spatial distribution of natural earthquakes and explosions
    图  2  天然地震和爆破的震级频次分布
    Figure  2.  Magnitude-frequency distribution of natural earthquakes and explosions

    本文所使用的原始数据为单台垂直向单通道地震记录,事件波形存储在seed格式的文件中,使用MSDP软件解压缩为sac格式,通过obspy读取后在python环境下分析。obspy是一个针对地震领域开发的python库,支持常用的地震数据文件格式的读写,可以实现文件转换功能,集成了许多地震学专有工具包,并提供了简单易用的通用接口(Beyreuther et al,2010)。

    选用数据的震中距在200 km以内,对每条记录截取100 s记录,截取起始点为初至震相前1 s,包含事件完整信息。原始数据的采样率为100 Hz。在数据预处理中,对信号进行高通滤波,保留1 Hz以上的信号并过滤掉低频噪声,之后进行去倾斜和归一化,最后通过短时傅里叶变换,将其转换为时频域的振幅谱。为了使数据分布均匀,再进行取对数操作和归一化操作,取对数操作在保持数据相对关系的同时,可缩小数据尺度,使数据的变化相对平稳,同时也削弱了模型的共线性和方差性等问题。转换后最终的数据为50×100的矩阵,对应频域1—50 Hz,时域0—100 s的对数振幅谱。图3为数据预处理的过程。

    图  3  数据预处理过程
    Figure  3.  The process of data preprocessing

    短时傅里叶变换(short-time Fourier transform)是一种常用的时频分析方法,其优点在于对信号具有一定程度的自适应性,可以还原输入信号在时频域的能量分布,对信号的低频和高频两端都有较为一致的局部描述,得到的非平稳信号的时频分布较为平滑(郝国成等,2019)。STFT是关于频率f和事件t的二元函数,即

    $$ {\rm{STF}}{{\rm{T}}_x}{\text{(}}{t{\text{,}}f} {\text{)}} {\text{=}} \mathop \int \nolimits x{\text{(}}\tau {\text{)}}{g^*}{\text{(}}{\tau {\text{-}} t} {\text{)}}{{\rm{e}}^{ - j2{\rm{\pi }}\tau f}}{\rm{d}}\tau {\text{,}} $$ (1)

    式中:x$ \tau $)为信号;$\tau $表示窗口$ {g}^{*}$中心坐标对应的时刻;t表示待分析时刻;$ {g}^{*} $为窗函数,本文选用汉明(Hamming)窗,以减少频谱泄露问题。窗长决定频谱图的时间分辨率和频率分辨率,窗长越长,截取的信号越长,变换后频率分辨率越高,时间分辨率越差。为了克服数据的关联性,相邻的窗部分重叠,本文选择的窗长为2 s,重叠1 s。

    STFT在时域将信号截取为多段,对每一段作傅里叶变换,求出每一段的频率特性,组合后得到时间-频率域的二维图像信号(时频域的振幅谱),该图像反映了信号能量在时频域的分布特征。

    本文分析了天然地震与爆破在频域和时域能量分布的差异。图4为地震与爆破的原始波形和时频谱的对比,由图可见:爆破的尾波有较多低频信号,含有面波成分(图4a);天然地震的频率较丰富,而爆破的能量集中在低频成分(图4b)。

    图  4  天然地震与爆破的波形(a)及时频域对数振幅谱(b)对比
    Figure  4.  Comparison of waveform (a) and time-frequency log amplitude spectrum (b) between earthquake and explosion

    天然地震和人工爆破可以看作随机信号,频率分布服从一定的统计规律。功率谱密度作为识别判据可以放大地震和爆炸的频域差异(隗永刚等,2019)。本文分别选取天然地震和爆破各300条记录,对归一化后的信号使用周期图法获得功率谱曲线,周期图法将随机序列观测数据视为能量有限的序列,直接计算离散傅里叶变换,取其幅值的平方并除以观测点数,作为真实功率谱的估计(王凤瑛,张丽丽,2006)。图5为天然地震与爆破的功率谱曲线(经平滑)对比,爆破在低频部分能量比例较高,随着频率增加逐渐衰减,地震在低频部分能量比例比爆破低,逐渐升高,在5 Hz左右达到峰值,随后衰减,天然地震能量随频率升高而衰减的速度比爆破缓慢。

    图  5  天然地震与爆破的功率谱对比
    Figure  5.  Power spectrum comparison between natural earthquakes and explosions

    由于天然地震震源深度较爆破更深,天然地震P波频率成分更复杂,爆破的震源较浅,一般在地表附近,其产生的振动的高频成分被松散地层波吸收,因此记录到的爆破P波的周期比较大(张帆,2006)。人工爆破的震源为膨胀源,其破裂机制较简单,频率成分也比较单一,能量快速释放,P波极为发育,当距震中距50 km以上时,会产生瑞雷面波(赵永等,1995)。

    人工神经网络受物神经网络概念的启发,由神经元层组成,这些神经元层通过权重值相互连接。每个神经元从上一层接收输入值,计算这些值的加权总和,加上一个偏置项,经过一个非线性的激活函数的作用,并将该函数值输出到下一层,即

    $$ {{\boldsymbol{x}}}^{l}{\text{=}}{\boldsymbol{f}}{\text{(}}{{\boldsymbol{w}}}^{l}{{\boldsymbol{x}}}^{l-1}{\text{+}}{{\boldsymbol{b}}}^{l}{\text{)}} {\text{,}} $$ (2)

    式中,xl-1为前一层输出特征图输入,wl为权重系数,bl为偏置.

    该过程模拟了生物神经元的行为。训练人工神经网络的过程包括调整权重和偏置项,以使输入信号与所需输出值一致。基于对输出误差进行量化的损失函数对权重和偏差的更新称为反向传播的过程。在此过程中,使用梯度下降算法来迭代更新权重和偏差以使量化损失函数最小化,每个完整数据传递和参数更新过程称为一个迭代周期。

    卷积神经网络与标准神经网络的主要区别在于卷积层的使用。卷积神经网络的前半部分一般为堆叠的卷积层,其间插入池化层构成特征提取器,后半部分将全连接层作为分类器,构成一个端到端的网络模型。卷积层中的神经元从前一层的较小的连续区域接收输入,乘以在整个输入数据上共享权重的滤波器(卷积核),一组共享的滤波器对输入图像进行卷积操作后得到一组特征图(张号奎等,2018)。卷积层输出可以表示为

    $$ \begin{array}{*{20}{c}} {{\boldsymbol{x}}_j^l {\text{=}} {\boldsymbol{f}}\left({ \displaystyle\sum \nolimits_{i \in {M_j}} {\boldsymbol{x}}_j^{l {\text{-}} 1} \times {\boldsymbol{w}}_{ij}^l {\text{+}} {\boldsymbol{b}}_j^l} \right)} \end{array} {\text{,}} $$ (3)

    式中,${\boldsymbol{x}}_j^l$表示第l层卷积层第j个卷积核对应的特征向量,Mj表示当前神经元的接受域,${\boldsymbol{w}}_{ij}^l $表示第l层第j个卷积核的第i个加权系数,${\boldsymbol{b}}_j^l $表示第l层第j个卷积核对应的偏置系数,f是一个非线性激活函数。

    相比全连接层,卷积层具有局部连接、权值共享和输入/输出数据结构化等特点,减少了参数数量和计算复杂度。池化层的操作类似于卷积层,只是不计算加权和,而是对每个区域的输出进行简单的处理,如取最大值(最大池化)或平均值(平均池化),这样减小了特征空间尺度,也保持特征对旋转的不变性,提高了网络的计算效率。

    卷积神经网络的训练过程分为数据前向传播和误差反向传播两个阶段。在前向传播阶段,数据由低层向高层传播。在反向传播阶段,误差从高层向低层进行传播,并根据误差对各层进行权值更新。反向传播算法主要基于梯度下降方法,网络参数首先被初始化为随机值,然后通过梯度下降法向训练误差减小的方向调整,优化卷积核参数k、下采样层网络权重β、全连接层网络权重w和各层的偏置参数b等。

    本文主要利用地震与爆破在频域和时域的差异,将地震事件分类问题转化为图像识别问题,将事件信号经过短时傅里叶变换得到的时频域对数振幅谱图作为输入,使用卷积神经网络实现天然地震事件和爆破事件的分类。

    经过多次试验,对卷积层数量、卷积核数量和全连接层单元数等超参数进行调整。对比不同模型的效果,我们最终选用一个20层的卷积神经网络,包括3个卷积层和2个全连接层,图6为神经网络的结构图示,表1给出了神经网络各层的参数。每个卷积层之后加入一个批量正则化层(batch normalization,缩写为BN)和一个池化层,卷积层使用修正线性单元(rectified linear unit,缩写为RELU)(式(4))。堆叠的卷积层用于特征提取,输入数据为50×100矩阵,经过第一个卷积层的16个卷积核作卷积,经原始输入转换为16个50×100的特征图(50×100×16),卷积层的输出经正则化后作为池化层的输入,池化层对16个特征图降采样,将相邻的四个元素合并为一个元素,经池化层处理后输出为16个25×50的特征图(25×50×16)。第一个池化层的输出作为第二个卷积层的输入,经过32个3×3×16的卷积核处理后转换为32个25×50的特征图(25×50×32),经过第二个池化层降采样处理后输出为12×25×32。第二个池化层的输出作为第三个卷积层的输入,经过卷积层和池化层后输出为6×12×64,合并为长度4 608的特征向量,这个向量提取了时频图的特征,作为全连接层的输入。BN层对中间层的输入(激活函数前)作标准化处理,使得输出服从均值为0,方差为1的正态分布,减少网络训练过程中的中间协变量迁移,使网络训练过程对参数初始值不再敏感(刘建伟等,2020),从而改善反向传播时训练速度变慢和低层神经网络的梯度消失等问题。两个全连接层实现了分类:第一个全连接层包含128个单元,之后设置一个批量正则化层,也采用RELU激活函数;第二个全连接层实现分类功能,采用Softmax激活函数 (柔性最大值传输函数),Softmax函数将第二个全连接层的输出映射成为(0,1)之间的值,两个输出节点之和为1,每个节点的输出可理解为输入信号为该类型的概率,概率最大结点作为分类结果。第一个全连接层前添加了一个丢弃层。丢弃层通过忽略一定比例的特征检测器(将一定比例的神经元的输出置零)减少过拟合现象(Krizhevsky et al,2017)。

    图  6  卷积神经网络结构
    Figure  6.  CNN network structure
    表  1  模型的超参数设置
    Table  1.  Super parameter setting of the model
    序号类型参数输出
    1 输入层 50×100×1
    2 卷积层 16个 3×3卷积核 50×100×16
    3 批量正则化层 50×100×16
    4 RELU激活层 50×100×16
    5 最大池化层 2×2 25×50×16
    6 卷积层 32个3×3×16卷积核 25×50×32
    7 批量正则化层 25×50×32
    8 RELU激活层 25×50×32
    9 最大池化层 2×2 12×25×32
    10 卷积层 64个3×3×32卷积核 12×25×64
    11 批量正则化层 12×25×64
    12 RELU激活层 12×25×64
    13 最大池化层 2×2 6×12×64
    14 Dropout层 比例50% 6×12×64
    15 全连接层 单元数128 128
    16 批量正则化层 128
    17 Relu激活层 128
    18 全连接层 单元数2 2
    19 Softmax激活层 2
    20 输出层
    下载: 导出CSV 
    | 显示表格
    $$ {\rm{RELU}}\left(x \right) {\text{=}} \left\{ {\begin{array}{*{20}{c}} {x{\text{,}}}&{x {\text{>}} \lambda }{\text{,}}\\ {0{\text{,}}}&{x {\text{≤}} \lambda }{\text{.}} \end{array}} \right.$$ (4)

    本文实验的硬件环境:英特尔i7-9850H处理器,64 G内存,英伟达RTX3000显卡,显存6 G。操作系统为Windows10,使用基于python的anaconda科学计算环境,模型的搭建和训练在tensowflow2.0开源框架进行。从1 800个样本中随机选取1 500个样本作为训练集,300个样本作为测试集。训练集是用于模型拟合的数据样本,测试集用来评估模型的泛化能力。

    本文选用交叉熵作为损失函数,交叉熵刻画的是实际输出(概率)与期望输出(概率)的距离,也就是交叉熵的值越小,两个概率分布就越接近,交叉熵H的算法为

    $$ H{\text{(}}{P,Q} {\text{)}} {\text{=}} - \sum \nolimits_x [{\text{(}}{P{\text{(}}x {\text{)}}{\rm{lg}}Q{\text{(}}x {\text{)}} {\text{+}} {\text{(}}{1 {\text{-}} P{\text{(}}x {\text{)}}} {\text{)}}{{\rm{lg}}{\text{(}}{1 {\text{-}} Q{\text{(}}x {\text{)}}} {\text{)}}} } ]{\text{,}} $$ (5)

    式中,P为期望输出,Q为实际输出。

    采用的训练算法为Adam算法(Kingma,Ba,2015),Adam优化算法是随机梯度下降算法(stochastic gradient descent)的扩展,可以更加有效地更新网络权重,使用动量和自适应学习率来加快收敛速度。

    为了评估学习率对训练效果的影响,分别将学习率依次设置为10−3,10−4,···,10−7进行训练,设置迭代次数为880次,图7a为不同学习率的损失曲线(经平滑),图7b为不同学习率的准确率曲线,准确率随着学习率的降低而下降,学习率较高时训练曲线收敛较快,最终准确率也较高,因此我们选取10−3作为初始学习率,并设置分段降低学习率。

    图  7  学习率测试
    (a) 损失函数曲线; (b) 准确率曲线; (c) 准确率随学习率变化
    Figure  7.  Learning rate test
    (a) Loss curves ;(b) Accuracy curves ;(c) Accuracy rate changes with learning rate

    交叉验证方法将训练数据集本身细分成不同的验证数据集去训练模型,验证集可以用于超参数的优化和对模型效果的初步测试。为了评估算法的性能和稳定性,用5折交叉验证检验,将1 500个样本均分为5组,每组300个样本,每轮选取一组作为验证集,将留下的4组作为训练集(图8)。5轮训练的验证集准确率可达95%—99.44%,验证集平均准确率为97.33%。

    图  8  5折交叉验证结果
    Figure  8.  Five fold cross validation results

    将全部训练集样本用于模型参数的训练,使用测试集的300个样本评估网络的泛化能力,最终测试的准确率达到98.02%,模型有较好的泛化能力和准确率。

    选取内蒙古地区2016年3月至2020年5月记录的天然地震和爆破事件中天然地震和爆破的单台垂直向记录各900条,经过预处理后转换成时频域的对数振幅谱分布图作为实验数据集,共包含1 800条样本。使用一个含有3个卷积层的卷积神经网络作为分类器,将实验数据集划分为训练集和测试集进行训练和检验。在训练集上使用5折交叉验证,结果显示,验证集的平均准确率达到97.3%。使用全部训练集数据对网络训练后,测试集的准确率达到98.02%。

    相比于传统机器学习,本文采用的模型在特征提取方面仅进行了简单处理和时频分析,保留了信号在时域和频域分布特征较完整的信息,使用卷积层自动提取特征,避免了使用先验的降维规则,以数据驱动的方式得到分类模型。实验结果表明本文使用的基于短时傅里叶变换和卷积神经网络的分类算法有效提取了地震信号的能量分布的特征,在地震事件分类任务上有较好的效果,与已有研究对比,本文模型达到了较高的准确率。

    在地震类型识别研究中,有多种数据处理方式和分类模型可供选择,在进一步研究中可以对比不同模型的效果。非天然事件包括爆破、塌陷和滑坡等多种类型,在现有数据条件下,本文仅研究了天然地震和爆破的分类问题,在积累了充足的多种类型数据样本后,可以开展更加细致的分类研究。

    评审专家的建议对稿件的提升帮助很大,作者在此表示感谢。

  • 图  1   研究区域的构造地质背景和地震分布(改自易桂喜等,2019Jiang et al,2020

    Figure  1.   The tectonic geological background and earthquake distribution in the studied area (modified from Yi et al,2019Jiang et al,2020

    图  2   长宁地区地震精定位结果分布特征

    (a) 构造和地震的三维分布特征;(b) EW向地震分布剖面图;(c) NS向地震分布剖面图

    Figure  2.   Distribution characteristics of precision earthquake location results in Changning area

    (a) Three-dimensional distribution characteristics of structures and earthquakes;(b) Earthquake distribution profile along EW direction;(c) Earthquake distribution profile along NS direction

    图  3   数据处理流程图

    Figure  3.   Data processing flow chart

    图  4   地震聚类分析及能量粒子簇提取

    (a) 利用轮廓系数评价不同聚类参数的聚类效果,以此获取最优聚类参数;(b) 利用密度聚类分析获取地震簇;(c) 利用K邻近分类对能量粒子进行分类;(d) 能量粒子簇分布特征

    Figure  4.   Earthquake clustering analysis and energy particle cluster extraction

    (a) The Silhouette coefficient used to evaluate the clustering effect of different clustering parameters so as to obtain the optimal clustering parameters;(b) Acquisition of earthquake clusters by DBSCAN;(c) Classification of energy particles using KNN;(d) Distribution characteristics of energy particle clusters

    图  5   长宁地震震源机制解分类结果及对应的震源节面玫瑰图

    Figure  5.   Classification results of focal mechanism solutions of Changning earthquake and the corresponding nodal rose diagrams

    图  6   长宁地区地震破裂面的空间展布特征

    Figure  6.   Spatial distribution characteristics of fracture planes of earthquakes in Changning area

    图  7   S波(a)和P波(b)速度结构和发震构造的分布特征

    Figure  7.   Distribution characteristics of S-wave (a) and P-wave (b) velocity structures and seismogenic structures

    表  1   长宁地区研究成果数据信息汇总

    Table  1   Data information from the research results of Changning area

    数据类型 数据描述 数据来源
    构造背景 基于区域地质与地震资料,结合地表调查结果,获取的节理面、构造线和应力场演化数据 常祖峰等(2 020
    褶皱构造 利用页岩气勘探钻井和反射地震资料进行构造分析,恢复了长宁背斜形成过程和构造
    地质背景
    He et al2 019
    发震断层 使用高分辨率地震反射剖面结合地质、地震和大地测量数据来揭示发震断层的三维分布 Lu et al2 021
    地震重定位 使用TomoDDMC方法联合反演后,重新定位2万1 711次地震,东西向、南北向和垂直方向的误差中值分别为0.201,0.232和0.633 km Zuo et al2 020
    微地震目录 应用最新发展的迁移学习震相识别技术、震相自动关联和定位技术,获取长宁微震目录,水平定位平均误差为(1.45±0.028) km 赵明等(2 021
    历史地震震源机制 对全国地震进行了矩张量反演,获得2 008—2 019年M≥3.0地震的震源机制解 郭祥云等(2 022
    地震序列震源机制 使用CAP (cut and paste)波形反演方法计算长宁地震16次MS≥3.6地震的震源机制解 易桂喜等(2 019
    综合震源机制解 借助于根据同一地震的多个震源机制解确定其中心解的方法,给出与所有震源机制解差别最小的中心震源机制解 刘敬光等(2 019
    速度结构 利用双差地震层析成像方法,获得了长宁—兴文地区高分辨率三维地壳vPvSvP/vS模型及地震位置 Zuo et al2 020
    下载: 导出CSV
  • 常祖峰,张艳凤,王光明,张世民,毛泽斌,刘昌伟. 2020. 2019年四川长宁 MS6.0地震的地质构造成因:区域性构造节理贯通、破裂结果[J]. 地球学报,41(4):469–480. doi: 10.3975/cagsb.2020.030202

    Chang Z F,Zhang Y F,Wang G M,Zhang S M,Mao Z B,Liu C W. 2020. The geological genesis of the 2019 Changning MS6.0 earthquake in Sichuan:Connecting and rupturing of regional structural joints[J]. Acta Geoscientica Sinica,41(4):469–480 (in Chinese).

    宫悦,王宇玺,梁明剑,龙锋,赵敏. 2020. 2019年四川长宁6.0级地震序列时空演化特征及其地震构造环境研究[J]. 地震,40(4):90–102. doi: 10.12196/j.issn.1000-3274.2020.04.007

    Gong Y,Wang Y X,Liang M J,Long F,Zhao M. 2020. Study on the spatio-temporal evolution characteristics and seismic structure environment of the 2019 M6.0 Changning Sichuan earthquake sequence[J]. Earthquake,40(4):90–102 (in Chinese).

    郭祥云,蒋长胜,韩立波,尹海权,赵志远. 2022. 中国大陆及邻区震源机制数据集(2009—2021年)[DB/OL]. [2022-05-01].https://data.earthquake.cn/datashare/report.shtml?PAGEID=datasourcelist&dt=ff8080827e4d6cf5017f1f5b440d0019.

    Guo X Y,Jiang C S,Han L B,Yin H Q,Zhao Z Y. 2022. Focal mechanism data set in Chinese mainland and its adjacent area (2009−2021)[DB/OL]. [2022-05-01].https://data.earthquake.cn/datashare/report.shtml?PAGEID=datasourcelist&dt=ff8080827e4d6cf5017f1f5b440d0019 (in Chinese).

    郭志,高星,路珍. 2020. 2019年6月17日四川长宁地震重定位及震源机制研究[J]. 地震学报,42(3):245–255. doi: 10.11939/jass.20190132

    Guo Z,Gao X,Lu Z. 2020. Relocation and focal mechanism inversion for the Changning,Sichuan,earthquake on 17 June 2019[J]. Acta Seismologica Sinica,42(3):245–255 (in Chinese).

    何登发,鲁人齐,黄涵宇,王晓山,姜华,张伟康. 2019. 长宁页岩气开发区地震的构造地质背景[J]. 石油勘探与开发,46(5):993–1006. doi: 10.11698/PED.2019.05.19

    He D F,Lu R Q,Huang H Y,Wang X S,Jiang H,Zhang W K. 2019. Tectonic and geological background of the earthquake hazards in Changning shale gas development zone,Sichuan Basin,SW China[J]. Petroleum Exploration and Development,46(5):1051–1064.

    胡晓辉,盛书中,万永革,卜玉菲,李振月. 2020. 2019年6月17日四川长宁地震序列震源机制与震源区震后构造应力场研究[J]. 地球物理学进展,35(5):1675–1681. doi: 10.6038/pg2020DD0378

    Hu X H,Sheng S Z,Wan Y G,Bu Y F,Li Z Y. 2020. Study on focal mechanism and post-seismic tectonic stress field of the Changning,Sichuan,earthquake sequence on June 17th 2019[J]. Progress in Geophysics,35(5):1675–1681 (in Chinese).

    胡幸平,崔效锋,张广伟,王甘娇,Zang A,史丙新,姜大伟. 2021. 长宁地区复杂地震活动的力学成因分析[J]. 地球物理学报,64(1):1–17. doi: 10.6038/cjg2021O0232

    Hu X P,Cui X F,Zhang G W,Wang G J,Zang A,Shi B X,Jiang D W. 2021. Analysis on the mechanical causes of the complex seismicity in Changning area,China[J]. Chinese Journal of Geophysics,64(1):1–17 (in Chinese).

    李大虎,詹艳,丁志峰,高家乙,吴萍萍,孟令媛,孙翔宇,张旭. 2021. 四川长宁 MS6.0地震震区上地壳速度结构特征与孕震环境[J]. 地球物理学报,64(1):18–35. doi: 10.6038/cjg2021O0241

    Li D H,Zhan Y,Ding Z F,Gao J Y,Wu P P,Meng L Y,Sun X Y,Zhang X. 2021. Upper crustal velocity and seismogenic environment of the Changning MS6.0 earthquake region in Sichuan,China[J]. Chinese Journal of Geophysics,64(1):18–35 (in Chinese).

    梁姗姗,徐志国,盛书中,张广伟,赵博,邹立晔. 2020. 2019年四川长宁6.0级地震主震及中强余震( MS≥4.0)的震源机制及其应力场[J]. 地震地质,42(3):547–561. doi: 10.3969/j.issn.0253-4967.2020.03.001

    Liang S S,Xu Z G,Sheng S Z,Zhang G W,Zhao B,Zou L Y. 2020. Focal mechanism solutions and stress field of the 2019 Changning,Sichuan mainshock and its moderate-strong aftershocks ( MS≥4.0)[J]. Seismology and Geology,42(3):547–561 (in Chinese).

    刘敬光,万永革,黄志斌,李振月,胡晓辉,李泽潇. 2019. 2019年6月17日四川长宁6.0级地震中心震源机制解及震源区构造应力场研究[J]. 震灾防御技术,14(3):677–685. doi: 10.11899/zzfy20190319

    Liu J G,Wan Y G,Huang Z B,Li Z Y,Hu X H,Li Z X. 2019. Study on central focal mechanism and its surrounding tectonic stress field of the Changning M6.0 earthquake in Sichuan[J]. Technology for Earthquake Disaster Prevention,14(3):677–685 (in Chinese).

    孙权,裴顺平,苏金蓉,刘雁冰,薛晓添,李佳蔚,李磊,左洪. 2021. 2019年6月17日四川长宁 MS6.0地震震源区三维速度结构[J]. 地球物理学报,64(1):36–53. doi: 10.6038/cjg2021O0246

    Sun Q,Pei S P,Su J R,Liu Y B,Xue X T,Li J W,Li L,Zuo H. 2021. Three-dimensional seismic velocity structure across the 17 June 2019 Changning MS6.0 earthquake,Sichuan,China[J]. Chinese Journal of Geophysics,64(1):36–53 (in Chinese).

    万永革. 2020. 震源机制与应力体系关系模拟研究[J]. 地球物理学报,63(6):2281–2296. doi: 10.6038/cjg2020M0472

    Wan Y G. 2020. Simulation on relationship between stress regimes and focal mechanisms of earthquakes[J]. Chinese Journal of Geophysics,63(6):2281–2296 (in Chinese).

    徐志国,梁姗姗,盛书中,张广伟,邹立晔,周元泽. 2020. 2019年四川长宁 MS6.0地震序列重定位和震源特征分析[J]. 地震学报,42(4):377–391. doi: 10.11939/jass.20190170

    Xu Z G,Liang S S,Sheng S Z,Zhang G W,Zou L Y,Zhou Y Z. 2020. Relocation and source characteristics of the 2019 Changning MS6.0 earthquake sequence[J]. Acta Seismologica Sinica,42(4):377–391 (in Chinese).

    易桂喜,龙锋,梁明剑,赵敏,王思维,宫悦,乔慧珍,苏金蓉. 2019. 2019年6月17日四川长宁 MS6.0地震序列震源机制解与发震构造分析[J]. 地球物理学报,62(9):3432–3447. doi: 10.6038/cjg2019N0297

    Yi G X,Long F,Liang M J,Zhao M,Wang S W,Gong Y,Qiao H Z,Su J R. 2019. Focal mechanism solutions and seismogenic structure of the 17 June 2019 MS6.0 Sichuan Changning earthquake sequence[J]. Chinese Journal of Geophysics,62(9):3432–3447 (in Chinese).

    赵明,唐淋,陈石,苏金蓉,张淼. 2021. 基于深度学习到时拾取自动构建长宁地震前震目录[J]. 地球物理学报,64(1):54–66. doi: 10.6038/cjg2021O0271

    Zhao M,Tang L,Chen S,Su J R,Zhang M. 2021. Machine learning based automatic foreshock catalog building for the 2019 MS6.0 Changning,Sichuan earthquake[J]. Chinese Journal of Geophysics,64(1):54–66 (in Chinesc).

    周志华. 2016. 机器学习[M]. 北京:清华大学出版社:114-115.

    Zhou Z H. 2016. Machine Learning[M]. Beijing:Tsinghua University Press:114-115 (in Chinesc).

    Al-Zoubi M B,Rawi M A. 2008. An efficient approach for computing Silhouette coefficients[J]. J Comput Sci,4(3):252–255. doi: 10.3844/jcssp.2008.252.255

    Bergen K J,Johnson P A,de Hoop M V,Beroza G C. 2019. Machine learning for data-driven discovery in solid Earth geoscience[J]. Science,363(6433):1299. doi: 10.1126/science.aau0323

    Brunsvik B,Morra G,Cambiotti G,Chiaraluce L,Di Stefano R,De Gori P,Yuen D A. 2021. Three-dimensional Paganica fault morphology obtained from hypocenter clustering (L’Aquila 2009 seismic sequence,Central Italy)[J]. Tectonophysics,804:228756. doi: 10.1016/j.tecto.2021.228756

    Gong L W,Zhang H,Chen S,Chen L J. 2022. Three-dimensional modeling of the Xichang crust in Sichuan,China by machine learning[J]. Appl Sci,12(6):2955. doi: 10.3390/app12062955

    Jia K,Zhou S Y,Zhuang J C,Jiang C S,Guo Y C,Gao Z H,Gao S S,Ogata Y,Song X D. 2020. Nonstationary background seismicity rate and evolution of stress changes in the Changning salt mining and shale-gas hydraulic fracturing region,Sichuan Basin,China[J]. Seismol Res Lett,91(4):2170–2181. doi: 10.1785/0220200092

    Jiang C S,Han L B,Long F,Lai G J,Yin F L,Bi J M,Si Z Y. 2021. Spatiotemporal heterogeneity of b values revealed by a data-driven approach for the 17 June 2019 MS6.0 Changning earthquake sequence,Sichuan,China[J]. Nat Hazards Earth Syst Sci,21(7):2233–2244. doi: 10.5194/nhess-21-2233-2021

    Jiang D W,Zhang S M,Ding R. 2020. Surface deformation and tectonic background of the 2019 MS6.0 Changning earthquake,Sichuan basin,SW China[J]. J Asian Earth Sci,200:104493. doi: 10.1016/j.jseaes.2020.104493

    Kaven J O,Pollard D D. 2013. Geometry of crustal faults:Identification from seismicity and implications for slip and stress transfer models[J]. J Geophys Res: Solid Earth,118(9):5058–5070. doi: 10.1002/jgrb.50356

    Keller J M,Gray M R,Givens J A. 1985. A fuzzy K-nearest neighbor algorithm[J]. IEEE Trans Syst Man Cybern,SMC-15(4):580–585. doi: 10.1109/TSMC.1985.6313426

    Kong Q K,Trugman D T,Ross Z E,Bianco M J,Meade B J,Gerstoft P. 2019. Machine learning in seismology:Turning data into insights[J]. Seismol Res Lett,90(1):3–14. doi: 10.1785/0220180259

    Kuang W H,Yuan C C,Zhang J. 2021. Real-time determination of earthquake focal mechanism via deep learning[J]. Nat Commun,12(1):1432. doi: 10.1038/s41467-021-21670-x

    Li W,Ni S D,Zang C,Chu R S. 2020. Rupture directivity of the 2019 MW5.8 Changning,Sichuan,China,earthquake and implication for induced seismicity[J]. Bull Seismol Soc Am,110(5):2138–2153. doi: 10.1785/0120200013

    Liu J Q,Zahradník J. 2020. The 2019 MW5.7 Changning earthquake,Sichuan Basin,China:A shallow doublet with different faulting styles[J]. Geophys Res Lett,47(4):e2019GL085408. doi: 10.1029/2019GL085408

    Lu R Q,Xu X W,He D F,John S,Liu B,Wang F Y,Tan X B,Li Y Q. 2017. Seismotectonics of the 2013 Lushan MW6.7 earthquake:Inversion tectonics in the eastern margin of the Tibetan Plateau[J]. Geophys Res Lett,44(16):8236–8243. doi: 10.1002/2017GL074296

    Lu R Q,He D F,Liu J Z,Tao W,Huang H Y,Xu F,Liu G S. 2021. Seismogenic faults of the Changning earthquake sequence constrained by high-resolution seismic profiles in the southwestern Sichuan Basin,China[J]. Seismol Res Lett,92(6):3757–3766. doi: 10.1785/0220200302

    Ouillon G,Ducorbier C,Sornette D. 2008. Automatic reconstruction of fault networks from seismicity catalogs:Three-dimensional optimal anisotropic dynamic clustering[J]. J Geophys Res: Solid Earth,113(B1):B01306.

    Reading A M,Cracknell M J,Bombardieri D J,Chalke T. 2015. Combining machine learning and geophysical inversion for applied geophysics[J]. ASEG Extended Abstracts,(1):1–5.

    Ross Z E,Trugman D T,Hauksson E,Shearer P M. 2019. Searching for hidden earthquakes in southern California[J]. Science,364(6442):767–771. doi: 10.1126/science.aaw6888

    Schubert E,Sander J,Ester M,Kriegel H P,Xu X W. 2017. DBSCAN revisited,revisited:Why and how you should (still) use DBSCAN[J]. ACM Trans Database Syst,42(3):19.

    Sun X L,Yang P T,Zhang Z W. 2017. A study of earthquakes induced by water injection in the Changning salt mine area,SW China[J]. J Asian Earth Sci,136:102–109. doi: 10.1016/j.jseaes.2017.01.030

    Trugman D T,Shearer P M. 2017. GrowClust:A hierarchical clustering algorithm for relative earthquake relocation,with application to the Spanish Springs and Sheldon,Nevada,earthquake sequences[J]. Seismol Res Lett,88(2A):379–391. doi: 10.1785/0220160188

    Wang S,Jiang G Y,Weingarten M,Niu Y F. 2020. InSAR evidence indicates a link between fluid injection for salt mining and the 2019 Changning (China) earthquake sequence[J]. Geophys Res Lett,47(16):e2020GL087603. doi: 10.1029/2020GL087603

    Wang S,Jiang G Y,Lei X L,Barbour A J,Tan X B,Xu C J,Xu X W. 2022. Three MW≥4.7 earthquakes within the Changning (China) shale gas field ruptured shallow faults intersecting with hydraulic fracturing wells[J]. J Geophys Res: Solid Earth,127(2):e2021JB022946. doi: 10.1029/2021JB022946

    Yin J X,Li Z F,Denolle M A. 2021. Source time function clustering reveals patterns in earthquake dynamics[J]. Seismol Res Lett,92(4):2343–2353. doi: 10.1785/0220200403

    Zhang B,Lei J S,Zhang G W. 2020. Seismic evidence for influences of deep fluids on the 2019 Changning MS6.0 earthquake,Sichuan basin,SW China[J]. J Asian Earth Sci,200:104492. doi: 10.1016/j.jseaes.2020.104492

    Zhang Z W,Liang C T,Long F,Zhao M,Wang D. 2022. Spatiotemporal variations of focal mechanism solutions and stress field of the 2019 Changning MS6.0 earthquake sequence[J]. Front Earth Sci,9:797907. doi: 10.3389/feart.2021.797907

    Zuo K Z,Zhao C P,Zhang H J. 2020. 3D crustal structure and seismicity characteristics of Changning-Xingwen area in the southwestern Sichuan Basin,China[J]. Bull Seismol Soc Ama,110(5):2154–2167. doi: 10.1785/0120200085

图(7)  /  表(1)
计量
  • 文章访问数:  182
  • HTML全文浏览量:  60
  • PDF下载量:  58
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-28
  • 修回日期:  2022-08-17
  • 网络出版日期:  2023-12-24
  • 刊出日期:  2023-12-24

目录

/

返回文章
返回