CN114152980B - 一种快速自动化产出震源机制解的方法与装置 - Google Patents

一种快速自动化产出震源机制解的方法与装置 Download PDF

Info

Publication number
CN114152980B
CN114152980B CN202111451140.7A CN202111451140A CN114152980B CN 114152980 B CN114152980 B CN 114152980B CN 202111451140 A CN202111451140 A CN 202111451140A CN 114152980 B CN114152980 B CN 114152980B
Authority
CN
China
Prior art keywords
seismic
wave
initial motion
facies
polarity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111451140.7A
Other languages
English (en)
Other versions
CN114152980A (zh
Inventor
赵明
肖卓伟
陈石
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
INSTITUTE OF GEOPHYSICS CHINA EARTHQUAKE ADMINISTRATION
Original Assignee
INSTITUTE OF GEOPHYSICS CHINA EARTHQUAKE ADMINISTRATION
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by INSTITUTE OF GEOPHYSICS CHINA EARTHQUAKE ADMINISTRATION filed Critical INSTITUTE OF GEOPHYSICS CHINA EARTHQUAKE ADMINISTRATION
Priority to CN202111451140.7A priority Critical patent/CN114152980B/zh
Publication of CN114152980A publication Critical patent/CN114152980A/zh
Application granted granted Critical
Publication of CN114152980B publication Critical patent/CN114152980B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282—Application of seismic models, synthetic seismograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616—Data from specific type of measurement
    • G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/64—Geostructures, e.g. in 3D data cubes
    • G01V2210/642—Faults

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种快速自动化产出震源机制解的方法与装置，其中该方法为：利用震相识别模型对实时地震监测数据进行震相拾取得到震相类别；利用初动极性识别模型对实时地震监测数据的初动极性进行识别得到P波的初动极性；对震相类别进行震相关联得到震相关联结果；根据震相关联结果得到震中位置；根据震相关联结果、震中位置和P波的初动极性对震相机制进行反演得到震源机制解。本发明通过利用震相识别模型和初动极性识别模型高效准确地获取P、S到时以及P波初动极性等信息，并基于此得到震相关联结果和震中位置，可以大幅降低反演震源机制的工作量，提升了震源机制解的反演速度。

Description

一种快速自动化产出震源机制解的方法与装置
技术领域
本发明涉及地震信息处理技术领域，特别是涉及一种快速自动化产出震源机制解的方法与装置。
背景技术
随着数字化网络项目的推进，我国已建立了1200多个地震台站，实现了全国数据实时传输和共享。地震发生以后，第一时间对实时监测数据进行处理，确定地震震源参数信息后，对于再交相关部门日常地震监测、应急决策、告警信息发布等相关操作具有十分重要的意义。地震震源参数除了大家耳熟能详的发震时刻、震源位置、震级以外，还包括震源机制解。震源机制解主要由断层走向、倾角、滑动角等节面解参数确定，准确快速的震源机制解对海啸预警、震源深度的确定、矩震级计算、断层方位确定以及监测断层活动都有极为重要的作用。目前，由于人工智能等技术在地震事件和震相自动检测、震相关联、地震定位等方面取得的进展，发震时刻、震源位置、震级等震源参数已经在很大程度上可以实现自动化处理，但在震源机制解的自动产出方面，进展要远远落后。这是因为人工智能方法，例如其中有代表性的深度学习技术实际上只是建立输入与输出之间的映射关系，一般不包含物理机制，而震源机制的求解是一个受物理规律约束的反演问题，其解(输出)往往也不是唯一的。因此，目前地震震源机制解的计算仍然主要通过人工完成，这在很大程度上加重了地震监测部门的日常工作负担，还会严重影响大地震发生时和发生后一段时间的应急处理、信息发布的速度。
发明内容
为了克服现有技术的不足，本发明的目的是提供一种应用人工智能技术自动获取全波形、P波初动极性、S/P振幅比等信息，并自动反演震源机制的方法与装置，以解决震源机制解反演效率低，降低解的非唯一性的问题。
为实现上述目的，本发明提供了如下方案：
一种快速自动化产出震源机制解的方法，包括：
获取实时地震监测数据；
利用震相识别模型自动对所述实时地震监测数据进行震相拾取得到震相类别；所述震相类别包括P波震相和S波震相；
利用初动极性识别模型自动对所述实时地震监测数据的初动极性进行识别得到P波的初动极性；
对所述震相类别进行震相关联得到震相关联结果；
根据所述震相关联结果得到震中位置；
根据所述震相关联结果、震中位置和P波的初动极性对震相机制进行反演得到震源机制解；所述震源机制解包括断层走向、倾角、滑动角。
优选的，所述震相识别模型是以三分量地震波形为输入，以震相类别为输出进行训练得到的；其中，震相识别模型在训练过程中的损失函数为：
Figure BDA0003386031500000021
式中，Yi′为二值化编码的标签，i＝1，2，3分别表示噪音、P波和S波三个类别，n为波形采样点数，Yi为U型神经网络模型最后一层的softmax函数计算得到的概率值，zi为U型神经网络模型最后一层的输出张量，Y′ij为第j个样本属于类i的真实概率，Yij为第j个样本属于类i的预测概率值。
优选的，所述初动极性识别模型是以P波到时前后2s的地震垂直分量波形数据为输入，以P波初动的极性为输出进行训练得到的。
优选的，所述根据所述震相关联结果、震中位置和P波的初动极性对震相机制进行反演得到震源机制解，包括：
利用全波形震源机制反演方法对震相机制进行反演得到震源机制解；其中，反演的目标函数为：
Figure BDA0003386031500000031
其中，
Figure BDA0003386031500000032
表示观测地震图，根据所述震相关联结果进行自动截取，
Figure BDA0003386031500000033
表示理论地震图，根据所述震中位置可计算得到，xcorr表示互相关，c表示地震图三分量，s为台站序号；
Figure BDA0003386031500000034
表示观测的P波初动极性，
Figure BDA0003386031500000035
表示理论的P波初动极性，且当
Figure BDA0003386031500000036
和
Figure BDA0003386031500000037
相同时，g值为0，否则为1；
Figure BDA0003386031500000038
表示观测地震图的S/P振幅比，
Figure BDA0003386031500000039
表示理论地震图的S/P振幅比，threshold为设定的阈值；s0表示观测波形和理论波形进行互相关之后获得的偏移量，α1表示第一权重系数，α2表示第二权重系数，α3表示第三权重系数，α4第四权重系数。
本发明还提供了一种快速自动化产出震源机制解的装置，包括：
地震监测数据获取模块，用于获取实时地震监测数据；
震相拾取模块，用于利用震相识别模型对所述实时地震监测数据进行震相拾取得到震相类别；所述震相类别包括P波震相和S波震相；
初动极性识别模块，用于利用初动极性识别模型对所述实时地震监测数据的初动极性进行识别得到P波的初动极性；
震相关联模块，用于对所述震相类别进行震相关联得到震相关联结果；
震中位置确定模块，用于根据所述震相关联结果得到震中位置；
震源机制解计算模块，用于根据所述震相关联结果、震中位置和P波的初动极性对震相机制进行反演得到震源机制解；所述震源机制解包括断层走向、倾角和滑动角。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明提供了一种快速自动化产出震源机制解的方法与装置，其中该方法为：利用震相识别模型对实时地震监测数据进行震相拾取得到震相类别；震相类别包括P波震相和S波震相；利用初动极性识别模型对实时地震监测数据的初动极性进行识别得到P波的初动极性；对震相类别进行震相关联得到震相关联结果；根据震相关联结果得到震中位置；根据震相关联结果、震中位置和P波的初动极性对震相机制进行反演得到震源机制解。本发明通过利用震相识别模型和初动极性识别模型高效准确地获取P、S到时以及P波初动极性等信息，并基于此得到震相关联结果和震中位置，可以大幅降低反演震源机制的工作量，提升了震源机制解的反演速度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明提供的实施例中的一种快速自动化产出震源机制解的方法流程图；
图2为本发明提供的实施例中的U型神经网络算法原理图；
图3为本发明提供的实施例中的smart motion神经网络算法原理图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在本文中提及“实施例”意味着，结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是，本文所描述的实施例可以与其它实施例相结合。
本申请的说明书和权利要求书及所述附图中的术语“第一”、“第二”、“第三”和“第四”等是用于区别不同对象，而不是用于描述特定顺序。此外，术语“包括”和“具有”以及它们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤、过程、方法等没有限定于已列出的步骤，而是可选地还包括没有列出的步骤，或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤元。
本发明的目的是提供一种快速自动化产出震源机制解的方法与装置以解决震源机制解反演效率低的问题。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
请参阅图1，一种快速自动化产出震源机制解的方法，包括：
步骤1：获取实时地震监测数据；
进一步的，本发明采用地震数据采集单元来采集实时地震监测数据。地震数据采集单元由一个Liss实时流、一个宽频带地震数据采集器和一个加速度传感器组成，其中，Liss实时流是架设在野外的台站通过4G网络实时传输并汇聚到数据中心并打包成seed格式的数据流，宽频带地震数据采集器采集的是宽频带地震仪实时记录的数据，加速度传感器记录的是实时三分量加速度数据。
步骤2：利用震相识别模型自动对所述实时地震监测数据进行震相拾取得到震相类别；所述震相类别包括P波震相和S波震相；
需要说明的是，在利用震相识别模型对所述实时地震监测数据进行震相拾取之前需要对所述实时地震监测数据进行预处理得到MiniSeed数据格式文件。进一步的，本发明可对获取的实时地震监测数据包按照时间序列生成对应的MiniSeed数据格式文件；
在本发明中，所述震相识别模型是以三分量地震波形为输入，以震相类别为输出进行训练得到的；其中，震相识别模型在训练过程中的损失函数为：
Figure BDA0003386031500000061
式中，Yi′为二值化编码的标签，i＝1,2,3分别表示噪音、P波和S波三个类别，n为波形采样点数，Yi为U型神经网络模型最后一层的softmax函数计算得到的概率值，zi为U型神经网络模型最后一层的输出张量，Y′ij为第j个样本属于类i的真实概率，Yij为第j个样本属于类i的预测概率值。
请参阅图2，在实际应用中，本发明的震相识别模型是基于U型神经网络构建的。U型网络的基本组件可分为增采样层和降采样层，其中降采样层由两个一维卷积层、一个池化层组成，中间会随机加入一些dropout层来防止过拟合；增采样层由一个转置卷积层、一个裁切层、一个卷积层组成，同样也可根据情况适当加入dropout层。网络输入为截取后的地震波形，其中，降采样层执行卷积和池化操作，用于提取震相的抽象特征，以解决震相定位问题，增采样层用于逐步恢复震相的细节特征，以解决震相分类的问题。最后通过激活函数计算P、S或噪音概率值，并与预设阈值比对即可确定该采样点的类别。
步骤3：利用初动极性识别模型自动对所述实时地震监测数据的初动极性进行识别得到P波的初动极性；具体的，所述初动极性识别模型是以P波到时前后2s的地震垂直分量波形数据为输入，以P波初动的极性为输出进行训练得到的。
请参阅图3，进一步的，本发明的初动极性识别模型是通过smart motion神经网络构建的。smart motion神经网络的网络输入为P波到时前后2s的地震垂直分量波形数据，其基本组件为5个依次相连的block，前两个block各包含2个一维卷积层，后3个block各包含3个一维卷积层，每个block的最后一个卷积层后面加上一个side layer层之后引出“侧输出”，一共得到5个尺度不同的“侧输出”(o1-o5)，本发明将这些侧输出聚合成最终输出(Fuse)，最后通过sigmoid激活函数计算分类概率值，并与预设阈值比较确定P波初动的极性为“UP”、“Down”或是待定。
步骤4：对所述震相类别进行震相关联得到震相关联结果；
在实际应用中，本发明可利用REAL进行自动震相关联，REAL是一种三维网格搜索算法，其功能是将按时间顺序排列的多台P波、S波到时与特定事件相关联，同时通过震相到时计数和与理论到时的走时残差来粗略定位地震位置。REAL算法首先需要将待研究区域划分为三维均匀网格，并以P波最早到达的台站为中心开始向四周逐网格搜索，当关联到的P波、S波到时个数达到预设阈值就认为是找到一个事件，同时到时计数最多的网格所在位置作为该事件的初步定位结果；如果有多个网格拥有着共同最大数目关联震相，则其中走时残差最小的网格为最优位置。
其具体的实现步骤为：
1.利用待研究区域的一维速度结构模型，根据研究区域的大小和最大定位深度，计算待研究区域内的参考走时表。
2.设置合理的REAL参数，包括搜索范围、网格尺寸大小，关联P波、S波个数阈值等。
3.以观测台站经纬度列表，震相识别模型拾取的各台站P、S波到时列表，以及参考走时表作为输入，运行REAL程序即可得到震相关联初步定位结果。
步骤5：根据所述震相关联结果得到震中位置；
本发明可给定台站位置，利用上述REAL的震相关联初步定位结果，以及一维速度参考模型，使用VELEST算法反演得到更加准确的事件位置(震中位置)，包括经纬度和深度。
步骤6：根据所述震相关联结果、震中位置和P波的初动极性对震相机制进行反演得到震源机制解；所述震源机制解包括断层走向、倾角、滑动角。
进一步的，所述步骤6具体包括：
利用全波形震源机制反演方法对震相机制进行反演得到震源机制解；其中，反演的目标函数为：
Figure BDA0003386031500000081
Figure BDA0003386031500000082
其中，
Figure BDA0003386031500000083
表示观测地震图，可根据所述震相关联结果进行自动截取获得，
Figure BDA0003386031500000084
表示理论地震图，可根据VELEST定位结果计算得到，即根据所述震中位置可计算得到，xcorr表示互相关，c表示地震图三分量(东西、南北、垂直)，s为台站序号；
Figure BDA0003386031500000085
表示观测的P波初动极性，
Figure BDA0003386031500000086
表示理论的P波初动极性，且当
Figure BDA0003386031500000087
和
Figure BDA0003386031500000088
相同时，g值为0，否则为1，
Figure BDA0003386031500000089
可由smart motion网络预测得到，
Figure BDA00033860315000000810
可通过对理论P波到时附近的一小段窗的数据进行求和得到求和值，并根据求和值的正负判断极性；
Figure BDA00033860315000000811
表示观测地震图的S/P振幅比，可根据预测的P、S波到时所计算的平均振幅得到，
Figure BDA00033860315000000812
表示理论地震图的S/P振幅比，threshold为设定的阈值；s0表示观测波形和理论波形进行互相关之后获得的偏移量，α1表示第一权重系数，α2表示第二权重系数，α3表示第三权重系数，α4第四权重系数。通过求解目标函数F的最小值即可得到最佳断层参数(走向、倾角、滑动角)。
最后本发明还需要对结果进行质量评级；进一步的，质量评级可根据反演波形匹配相关度(可计算得到一个相关系数)、反演残差(一个值)、结果参数分布、台站数目及分布情况进行评价，并对以上每一项结合实际应用人为设定阈值，将质量评级分为A(非常好，三个以上条件达到阈值)，B(好，两个以上)，C(一般，一个以上)，D(差，均低于阈值)四个质量等级，当出现C和D评级时，需要重新选取训练样本对震相识别模型和初动极性识别模块进行训练，以提高后续反演精度。
本发明公开了一种快速自动化产出震源机制解的方法，包括：获取实时地震监测数据。对所述实时地震监测数据进行预处理得到MiniSeed数据格式文件。利用U型神经网络模型对所述MiniSeed数据格式文件进行震相拾取得到震相类别；利用smart motion神经网络模型对所述P波的初动极性进行识别得到P波的初动极性；利用REAL对所述P、S波震相进行自动震相关联与初步事件定位；利用震相关联结果、台站经纬度、给定的一维速度参考模型等信息，反演得到地震的震中位置；利用全波形震源机制反演方法自动反演震源机制解。本发明通过引入了U型神经网络、smart motion神经网络、REAL自动震相关联等先进技术，可以在30s左右就完成反演震源机制的所有数据准备，再自动调用全波形匹配反演算法，总共在2～10分钟就可反演出震源机制解。
本发明还提供了一种快速自动化产出震源机制解的装置，包括：
地震监测数据获取模块，用于获取实时地震监测数据；
震相拾取模块，用于利用震相识别模型对所述实时地震监测数据进行震相拾取得到震相类别；所述震相类别包括P波震相和S波震相；
初动极性识别模块，用于利用初动极性识别模型对所述实时地震监测数据的初动极性进行识别得到P波的初动极性；
震相关联模块，用于对所述震相类别进行震相关联得到震相关联结果；
震中位置确定模块，用于根据所述震相关联结果得到震中位置；
震源机制解计算模块，用于根据所述震相关联结果、震中位置和P波的初动极性对震相机制进行反演得到震源机制解；所述震源机制解包括断层走向、倾角和滑动角。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明通过利用震相识别模型和初动极性识别模型高效准确地获取P、S到时以及P波初动极性等信息，并基于此得到震相关联结果和震中位置，可以大幅降低反演震源机制的工作量，提升了震源机制解的反演速度。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的方法而言，由于其与实施例公开的装置相对应，所以描述的比较简单，相关之处参见装置部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (4)

1.一种快速自动化产出震源机制解的方法，其特征在于，包括：
获取实时地震监测数据；
利用震相识别模型自动对所述实时地震监测数据进行震相拾取得到震相类别；所述震相类别包括P波震相和S波震相；所述震相识别模型是以三分量地震波形为输入，以震相类别为输出进行训练得到的；其中，震相识别模型在训练过程中的损失函数为：
Figure FDA0003583162500000011
式中，Y′i为二值化编码的标签，i＝1,2,3分别表示噪音、P波和S波三个类别，n为波形采样点数，Yi为U型神经网络模型最后一层的softmax函数计算得到的概率值，zi为U型神经网络模型最后一层的输出张量，Y′ij为第j个样本属于类i的真实概率，Yij为第j个样本属于类i的预测概率值；
利用初动极性识别模型自动对所述实时地震监测数据的初动极性进行识别得到P波的初动极性；
对所述震相类别进行震相关联得到震相关联结果；
根据所述震相关联结果得到震中位置；
根据所述震相关联结果、震中位置和P波的初动极性对震相机制进行反演得到震源机制解；所述震源机制解包括断层走向、倾角、滑动角。
2.根据权利要求1所述的一种快速自动化产出震源机制解的方法，其特征在于，所述初动极性识别模型是以P波到时前后2s的地震垂直分量波形数据为输入，以P波初动的极性为输出进行训练得到的。
3.根据权利要求1所述的一种快速自动化产出震源机制解的方法，其特征在于，所述根据所述震相关联结果、震中位置和P波的初动极性对震相机制进行反演得到震源机制解，包括：
利用全波形震源机制反演方法对震相机制进行反演得到震源机制解；其中，反演的目标函数为：
Figure FDA0003583162500000021
其中，
Figure FDA0003583162500000022
表示观测地震图，
Figure FDA0003583162500000023
表示理论地震图，xcorr表示互相关，c表示地震图三分量，s为台站序号；
Figure FDA0003583162500000024
表示观测的P波初动极性，
Figure FDA0003583162500000025
表示理论的P波初动极性，且当
Figure FDA0003583162500000026
和
Figure FDA0003583162500000027
相同时，g值为0，否则为1；
Figure FDA0003583162500000028
表示观测地震图的S/P振幅比，
Figure FDA0003583162500000029
表示理论地震图的S/P振幅比，threshold为设定的阈值；s0表示观测波形和理论波形进行互相关之后获得的偏移量，α1表示第一权重系数，α2表示第二权重系数，α3表示第三权重系数，α4第四权重系数。
4.一种快速自动化产出震源机制解的装置，其特征在于，包括：
地震监测数据获取模块，用于获取实时地震监测数据；
震相拾取模块，用于利用震相识别模型对所述实时地震监测数据进行震相拾取得到震相类别；所述震相类别包括P波震相和S波震相；所述震相识别模型是以三分量地震波形为输入，以震相类别为输出进行训练得到的；其中，震相识别模型在训练过程中的损失函数为：
Figure FDA00035831625000000210
式中，Y′i为二值化编码的标签，i＝1,2,3分别表示噪音、P波和S波三个类别，n为波形采样点数，Yi为U型神经网络模型最后一层的softmax函数计算得到的概率值，zi为U型神经网络模型最后一层的输出张量，Y′ij为第j个样本属于类i的真实概率，Yij为第j个样本属于类i的预测概率值；
初动极性识别模块，用于利用初动极性识别模型对所述实时地震监测数据的初动极性进行识别得到P波的初动极性；
震相关联模块，用于对所述震相类别进行震相关联得到震相关联结果；
震中位置确定模块，用于根据所述震相关联结果得到震中位置；
震源机制解计算模块，用于根据所述震相关联结果、震中位置和P波的初动极性对震相机制进行反演得到震源机制解；所述震源机制解包括断层走向、倾角和滑动角。
CN202111451140.7A 2021-12-01 2021-12-01 一种快速自动化产出震源机制解的方法与装置 Active CN114152980B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111451140.7A CN114152980B (zh) 2021-12-01 2021-12-01 一种快速自动化产出震源机制解的方法与装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111451140.7A CN114152980B (zh) 2021-12-01 2021-12-01 一种快速自动化产出震源机制解的方法与装置

Publications (2)

Publication Number Publication Date
CN114152980A CN114152980A (zh) 2022-03-08
CN114152980B true CN114152980B (zh) 2022-06-10

Family

ID=80455204

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111451140.7A Active CN114152980B (zh) 2021-12-01 2021-12-01 一种快速自动化产出震源机制解的方法与装置

Country Status (1)

Country Link
CN (1) CN114152980B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114994754B (zh) * 2022-05-18 2023-03-24 成都理工大学 基于直达波和深度震相初动极性的震源机制联合反演方法
CN115373029B (zh) * 2022-10-25 2023-01-31 中国科学院地质与地球物理研究所 基于深度学习的实时微地震震源机制计算方法及系统
CN115826042B (zh) * 2022-11-22 2023-05-12 中国地震局地球物理研究所 一种边云端结合的分布式地震数据处理方法与装置
CN117970488B (zh) * 2024-04-02 2024-05-31 山东省地震局 一种三维复杂地学模型中地震精确定位方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107918157B (zh) * 2016-10-08 2019-07-23 中国石油化工股份有限公司 三分量p波初动震源机制反演方法及装置
CN107843624B (zh) * 2017-09-14 2020-06-02 中国肉类食品综合研究中心 用于检测食品中菌落总数的检测装置及其检测方法
CN107843921A (zh) * 2017-11-03 2018-03-27 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种微地震震源机制反演方法
CN111723329B (zh) * 2020-06-19 2023-04-07 南京大学 一种基于全卷积神经网络的震相特征识别波形反演方法

Also Published As

Publication number Publication date
CN114152980A (zh) 2022-03-08

Similar Documents

Publication Publication Date Title
CN114152980B (zh) 一种快速自动化产出震源机制解的方法与装置
CN108805269B (zh) 一种基于lstm循环神经网络拾取震相到时的方法
CN114063153B (zh) 一种自动反演震源机制解的方法与装置
CN115826042B (zh) 一种边云端结合的分布式地震数据处理方法与装置
CN114509811B (zh) 一种基于深度学习的单台站后方位角估计方法与装置
CN112799128B (zh) 一种地震信号检测和震相提取的方法
CN111983676A (zh) 一种基于深度学习的地震监测方法及装置
CN115032693B (zh) 一种强震的前震自动识别方法及装置
CN110501742A (zh) 一种采用Boosting集成学习算法对地震事件进行区分的方法
CN106874950A (zh) 一种暂态电能质量录波数据的识别分类方法
CN113985479A (zh) 一种基于边缘计算的地震信息处理方法与装置
CN112884134A (zh) 面向地震震相识别的基于时域的卷积神经网络模型及应用
CN111898644A (zh) 一种无故障样本下航天液体发动机健康状态智能识别方法
CN113920255B (zh) 基于点云数据的高效测绘系统
CN116594057B (zh) 一种基于深度学习和边缘计算的地震预警方法与装置
CN114330120A (zh) 一种基于深度神经网络预测24小时pm2.5浓度的方法
Zhang et al. Learnable maximum amplitude structure for earthquake event classification
CN117454100A (zh) 一种面向工程结构的缺失数据重构方法
CN117593601A (zh) 一种基于深度学习的水尺验潮方法
CN115952410A (zh) 一种基于深度学习的滑坡灾害检测系统
CN111624586A (zh) 基于特征提取和神经网络的水声目标测距方法
CN112380198B (zh) 一种基于深度学习的地震接收函数自动挑选方法
CN113806413B (zh) 一种轨迹筛选分类方法、装置、计算机设备和存储介质
CN108051853B (zh) 一种基于单台站初至p波的震中距快速估算方法
CN109350072A (zh) 一种基于人工神经网络的步频探测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant