Generalization of PhaseNet in Shandong and its application to the Changqing M4.1 earthquake sequence

Waveforms of seismic events, extracted from January 2019 to December 2021 were used to construct a test dataset to investigate the generalizability of PhaseNet in the Shandong region. The results show that errors in the picking of seismic phases(P-and S-waves) had a broadly normal distribution, mainly concentrated in the ranges of-0.4–0.3 s and-0.4–0.8 s, respectively. These results were compared with those published in the original PhaseNet article and were found to be approximately 0.2–0.4 s l...     

基于模板匹配定位法的江苏盐城附近微震检测和构造分析

利用基于GPU加速的匹配定位法和双差定位法,对江苏盐城及邻区18个台站记录的2009~2018年共10年的连续地震资料进行分析。首先从台网目录中挑选211个地震事件作为模板事件,使用匹配定位技术对江苏盐城附近连续10年的地震进行检测和识别,共识别出1349个地震事件,约为台网目录地震事件的3倍,最小完备震级由台网目录的M_L1.9降为M_L1.2。然后利用双差定位法对检测到的地震事件进行精定位,精定位的结果揭示：建湖地区的地震密集带与洪泽-沟墩断裂有关,震源深度优势分布为5~20km,断裂两侧震源深度有显著差异,断裂带倾向NW;射阳震群震源深度比建湖震群有所加深,优势分布为10~25km,震源深度由南东向西北逐渐变浅;宝应地区地震丛集分布;东台地区由于模板事件相对较少,扫描定位后,地震事件在陈家堡-小海断裂带附近零星分布。研究结果为研究盐城地区的地震活动性、发震断层的深部构造提供了基础数据支撑。     

地震检测与震相自动拾取研究

针对微震事件易受噪声干扰等特点,本文将STA/LTA方法和基于方差的AIC方法（var-AIC）相结合,在震相到时初步拾取的基础上,使用台站的德洛内（Delaunay）三角剖分及台站间最大走时差约束来减少噪声干扰的影响. 利用到时进行地震定位之后,根据台站预测到时,在设定的时间窗内对地震震相进行更精细的分析. 特别是针对微震事件信噪比低的特点,设计了基于偏振分析的拾取函数,根据窗内STA/LTA方法和var-AIC方法的拾取结果自动选择合适的值作为震相到时. 最后,对西昌流动地震台阵2013年304个单事件波形数据的分析处理和检验结果表明,本文方法较传统方法具有更高的地震事件检测能力和更高的震相拾取精度.      

强噪声环境下基于信噪比的地震P波到时自动提取方法

大数据量、强噪声环境给地震P波到时的自动提取带来很大挑战.针对此问题,本文通过构建特殊的特征函数,建立SNR与STA/LTA的内在联系,提出两种基于SNR的地震P波到时自动提取方法,即基于SNR的STA/LTA方法与基于SNR的综合方法.这两种方法分别是运用SNR概念对传统STA/LTA方法和STA/LTA与AIC综合方法的改进.仿真分析结果表明:对于弱噪声环境(10dB)和一般噪声环境(6dB),本文方法较传统STA/LTA方法对地震P波到时提取的准确度更高;而对于强噪声环境(3dB),本文方法仍能准确提取地震P波到时,而传统STA/LTA方法则出现了较大的误判率(10%)与漏判率(65%).本文方法为STA/LTA赋予了明确的物理意义,使其阈值的选取建立在严密的数学推导之上.另外,本文方法在进行地震P波到时自动提取的同时,兼具数据预处理功能,无需额外的基线校正或高通滤波,因而具有较好的实时性.     

基于深度学习到时拾取自动构建长宁地震前震目录

将深度学习到时拾取、震相关联技术与传统定位方法联系起来,构建一套连续波形自动化处理与地震目录自动构建流程,对于高效充分利用地震资料,提升微震检测能力具有十分重要的意义.我们应用最新发展的迁移学习震相识别技术、震相自动关联技术,对长宁M S6.0地震震中附近21个台站震前半个月(6月1日—6月17日)的连续记录波形进行P、S震相识别、震相自动关联和初步定位,并应用传统绝对定位和相对定位技术得到了长宁地震震前微震活动的绝对和相对定位目录.其中绝对定位目录能在较小的误差范围匹配85%的人工处理目录,其发震时刻平均误差为0.36±0.07 s,震级平均误差为0.15±0.024级,水平定位平均误差为1.45±0.028 km,其识别的1.0级以下微震数目是人工的8倍以上,将长宁地震震前微震目录的检测下限提升至M L-1左右,证明了基于深度学习到时识取和REAL(Rapid Earthquake Association and Location,快速震相关联和定位技术)震相自动关联来构建微震目录具有较好的实用性.我们的自动地震目录揭示了长宁M S6.0主震所发生的区域震前异常频繁的微震活动,以及与区域内盐矿注水井的关联性,更好地描绘了这些微震活动的时空演化特征,其空间活动性分布特征与长宁M S6.0余震序列的分布一致.     

东祁连山北缘断裂带基于深度学习的密集台阵地震事件快速检测与定位研究

为监测东祁连山北缘断裂带附近的地震活动性,布设包含240台短周期地震仪的面状密集台阵,进行约30 d的连续观测。首先使用基于深度学习的多台站地震事件检测算法（CNNDetector）进行地震事件检测,然后使用震相拾取网络（PhaseNet）对地震事件进行P波和S波到时拾取,其次使用震相关联算法（REAL）进行震相关联及初定位,最后使用双差定位（hypoDD）进行地震重定位,最终的精定位地震目录中共有517个地震。在密集台阵观测期间,中国地震台网正式地震目录中共有39个位于台阵内的地震事件,相比而言,密集台阵检测到大量小于0级的地震。因此通过布设密集台阵,可提高活动断裂微地震活动性的监测能力。与历史地震空间分布相比,密集台阵地震精定位分布具有较好的一致性,表现出更明显的线性分布特征。基于地震分布,发现研究区域存在与地表断层迹线走向不同的隐伏活跃断裂。     

近震信号检测参数设置研究

地震信号检测是进行各种地震数据分析和处理的首要任务,STA/LTA方法具有算法简单、便于实时处理等特点,被广泛应用于地震信号检测.结合实际震例数据研究STA/LTA方法进行地震信号检测的各种影响因素,得到该方法进行检测时最合理的参数设置范围.     

基于机器学习算法的银川台阵资料自动化震相分析和定位

近年来快速发展的机器学习算法显著提高了震相拾取的精度和效率.采用卷积神经网络和递归神经网络的震相识别方法对银川台阵2019年6～7月的连续波形数据进行事件检测和P、S震相拾取,并通过快速震相关联和事件定位得到了银川地区较全的地震目录.结果表明,当震相数小于10时,虽然可以检测出较多事件,但分布呈弥散状,与区域地震活动特...     

A high-resolution seismic catalog for the 2021 MS6.4/MW6.1 Yangbi earthquake sequence,Yunnan, China: Application of AI picker and matched filter

We present a high-resolution seismic catalog for the 2021 M_S6.4/M_W6.1 Yangbi sequence. The catalog has a time range of 2021-05-01 to 2021-05-28, and contains ~8,000 well located events. It captures the features of the whole foreshock sequence and the early aftershocks. We designed a detection strategy incorporating both an artificial intelligent (AI) picker and a matched filter algorithm. Here, we adopt a hybrid AI method incorporating convolutional and recurrent neural network (CNN & RNN) for event detection and phase picking respectively (i.e. CERP), a light-weight AI picker that can be trained with small volume of data. CERP is first trained with detections from a STA/LTA and Kurtosis-based method called PAL, and then construct a rather complete template set of ~4,000 events. Finally, the matched filter algorithm MESS augments the initial detections and measures differential travel times with cross-correlation, which finally results in precise relocation. This process gives 9,026 detections, among which 7,943 events can be well relocated. The catalog shows as expected power-law distribution of frequency magnitude and reveals detailed pattern of seismicity evolution. The main features are: (1) the foreshock sequence images simple fault geometry with consistent strike, but also show a variable event depth along strike; (2) the mainshock ruptures the same fault of the foreshock sequence and activate conjugate faults further to the southeast; (3) complex seismicity are developed in the post-seismic period, indicating complex triggering mechanisms. Thus, our catalog provides a reliable basis for further investigations, such as b-value studies, rupture process, and triggering relations.     

用于地震预警的P波震相到时自动拾取

P波震相的自动拾取可用于地震预警中地震事件判别和地震定位,是实现基于地震台网地震预警的首要条件.针对地震预警中P波震相拾取的特点,本文发展了一套基于长短时平均(STA/LTA)和池赤准则(AIC)算法的多步骤P波自动拾取技术,应用Delaunay三角剖分提出了一种非几何相关的干扰信号剔除方法,并应用福建省数字地震台网记录对方法进行了验证,目前方法已经用到了福建省地震预警试验系统中.     

岩石超声检测信号初至的AIC自动拾取方法

岩石超声检测中最重要的一个环节是初至的拾取,然而该项工作往往费时费力,拾取精度受人为因素影响较大。为提高声波速度检测、声发射定位、以及超声层析成像的应用效率和精度,本研究将地震学中应用比较广泛的AIC初至自动提取技术引入到岩石超声检测中,并进行了适当改进。利用改进前后的AIC方法,自动拾取仿真信号和实际信号的初至,并利用长短时窗比方法(STA/LTA)和手动方法拾取了初至,同时分别与设定的实际初至进行对比。根据实验结果,对于信噪比较低的信号AIC方法要优于STA/LTA方法;改进前的AIC方法适用于起跳干脆、幅度变化大的信号,而改进后的AIC方法则适用于起跳较平缓的信号,且拾取到的初至与手动拾取的初至更加接近。      

基于U形卷积神经网络的震相识别与到时拾取方法研究

精确获取震相到时是地震定位和地震走时成像等研究的重要基础.近年来,随着地震台站的不断加密,地震台网监测到的地震数量成倍增长,发展快速、准确、适用性强的震相到时自动拾取算法是地震行业的迫切需求.本文在前人工作基础上,发展了Pg、Sg震相自动识别与到时拾取的U网络算法（Unet_cea）,使用汶川余震和首都圈地震台网记录的89344个不同震级、不同信噪比的样本进行训练和测试.研究表明,U网络能够较好地识别Pg、Sg震相类型和拾取到时,Pg、Sg震相的正确识别率分别为81%和79.1%,与人工标注到时的均方根误差分别为0.41 s和0.54 s.U网络在命中率、均方根误差等性能指标上均明显优于STA/LTA和峰度分析自动拾取方法.研究获得的最优模型可以为区域地震台网的自动处理提供辅助.     

Automatic onset phase picking for portable seismic array observation

Automatic phase picking is a critical procedure for seismic data processing, especially for a huge amount of seismic data recorded by a large-scale portable seismic array. In this study is presented a new method used for automatic accurate onset phase picking based on the proporty of dense seismic array observations. In our method, the Akaike＇s information criterion （AIC） for the single channel observation and the least-squares cross-correlation for the multi-channel observation are combined together. The tests by the seismic array observation data after triggering with the short-term average/long-term average （STA/LTA） technique show that the phase picking error is less than 0.3 s for local events by using the single channel AIC algorithm. In terms of multi-channel least-squares cross-correlation technique, the clear teleseismic P onset can be detected reliably. Even for the teleseismic records with high noise level, our algorithm is also able to effectually avoid manual misdetections.     

基于深度学习卷积神经网络的地震波形自动分类与识别

发展高效、高精度、普适性强的自动波形拾取算法在地震大数据时代背景下显得越来越重要.波形自动拾取算法的主要挑战来自如何适应不同区域的不同类型地震事件的分类与筛选.本文针对地震事件-噪音分类这一问题,使用13839个汶川地震余震事件建立数据集,应用深度学习卷积神经网络(CNN)方法进行训练,并用8900个新的汶川余震事件作为检测数据集,其训练和检测准确率均达到95%以上.在对连续波形的检测中,CNN方法在精度和召回率上优于STA/LTA和Fbpicker传统方法,并能找出大量人工挑选极易遗漏的微震事件.最后,我们应用训练好的最优模型对选自全国台网的441个台站8天的连续波形数据进行了识别、到时挑取及与参考地震目录关联,CNN检出7016段波形,用自动挑选算法拾取到1380对P,S到时,并与540个地震目录事件成功关联,对1级以上事件总体识别准确率为54%,二级以上为80%,证明了CNN模型具有泛化能力,初步展示了CNN在发展兼具效率、精度、普适性算法,实时地震监测等应用上具有巨大潜力.     

一种地震P波和S波初至时间自动拾取的新方法

地震P波、S波初至时间的拾取是地震波分析的一项基础性工作.本文提出了一种新的地震波初至时间自动拾取的方法:首先,把地震波的三分量时程曲线变换为一组空间向的能量变化率时程曲线;然后对能量变化率时程曲线进行STA/LTA(Short Time Average/Long Time Average,短时间的均值/长时间的均值)处理,拾取地震P波和S波的大致初至时间;最后提出采用一种二次方自回归模型对初至附近的能量变化率曲线进行二次方自回归处理,精确拾取出P波和S波的初至时间.本文采用了10组芦山地震的记录数据和150组汶川地震的记录数据对此方法的可靠性进行了检验.以人工拾取结果为参考,此方法具有很高的准确率和稳定性,同时,相比于常用的STA/LTA方法和AIC(Akaike Information Criterion,Akaike信息准则)方法,此方法在计算时间效率方面稍微逊色,但是对S波初至时间的拾取精度和可靠性更高.此方法丰富了地震P波、S波初至时间的自动拾取方法.     

基于深度学习的地震检测模型在区域台网的泛化性研究

将识别地震的深度学习算法PhaseNet应用于四川台网和首都圈台网,对该模型的泛化能力进行了测试和评估.首先利用2010年1月至2018年10月首都圈台网199个地震台站记录的29 328个事件(ML0～ML4)所对应的126761段事件波形,以及2019年4-9月四川及邻省部分台网227个地震台站记录的16595个事...     

Active Source Seismic Identification and Automatic Picking of the P-wave First Arrival Using a Convolutional Neural Network

In seismic data processing, picking of the P-wave first arrivals takes up plenty of time and labor, and its accuracy plays a key role in imaging seismic structures. Based on the convolution neural network (CNN), we propose a new method to pick up the P-wave first arrivals automatically. Emitted from MINI28 vibroseis in the Jingdezhen seismic experiment, the vertical component of seismic waveforms recorded by EPS 32-bit portable seismometers are used for manually picking up the first arrivals (a total of 7242). Based on these arrivals, we establish the training and testing sets, including 25,290 event samples and 710,616 noise samples (length of each sample:2s). After 3,000 steps of training, we obtain a convergent CNN model, which can automatically classify seismic events and noise samples with high accuracy (> 99%). With the trained CNN model, we scan continuous seismic records and take the maximum output (probability of a seismic event) as the P-wave first arrival time. Compared with STA/LTA (short time average/long time average), our method shows higher precision and stronger anti-noise ability, especially with the low SNR seismic data. This CNN method is of great significance for promoting the intellectualization of seismic data processing, improving the resolution of seismic imaging, and promoting the joint inversion of active and passive sources.     

江苏数字地震台网中心记录波形分析

对江苏数字地震台网的波形资料,从震相、频谱、振幅比3个方面进行对比分析,找出天然地震、人工爆破、塌陷的不同之处,能够在地震定位过程中快速识别,并将此方法在江苏地震台网实际工作中应用。     

Comparison of two earthquake early warning location methods

According to earthquake catalog records of Fujian Seismic Network, the T _now method and the four-station continuous location method put forward by Jin Xing are inspected by using P-wave arrival information of the first four stations in each earthquake. It shows that the four-station continuous location method can locate more seismic events than the T _now method. By analyzing the results, it is concluded that the reason for this is that the T _now method makes use of information from stations without being triggered, while some stations failed to be reflected in earthquake catalog because of discontinuous records or unclear records of seismic phases. For seismic events whose location results can be given, there is no obvious difference in location results of the two methods and positioning deviation of most seismic events is also not significant. For earthquakes outside the network, the positioning deviation may amplify as the epicentral distance enlarges, which may relate to the situation that the seismic stations are centered on one side of epicenter and the opening angle between seismic stations used for location and epicenter is small.     

Local Seismic Events in the Area of Poland Based on Data from the PASSEQ 2006–2008 Experiment

PASSEQ 2006–2008 (Passive Seismic Experiment in TESZ; Wilde-Piórko et al. 2008) was the biggest passive seismic experiment carried out so far in the area of Central Europe (Poland, Germany, the Czech Republic and Lithuania). 196 seismic stations (including 49 broadband seismometers) worked simultaneously for over two years. During the experiment, multiple types of data recorders and seismometers were used, making the analysis more complex and time consuming. The dataset was unified and repaired to start the detection of local seismic events. Two different approaches for detection were applied for stations located in Poland. The first one used standard STA/LTA triggers (Carl Johnson’s STA/LTA algorithm) and grid search to classify and locate the events. The result was manually verified. The second approach used Real Time Recurrent Network (RTRN) detection (Wiszniowski et al. 2014). Both methods gave similar results, showing four previously unknown seismic events located in the Gulf of Gdańsk area, situated in the southern Baltic Sea. In this paper we discuss both detection methods with their pros and cons (accuracy, efficiency, manual work required, scalability). We also show details of all detected and previously unknown events in the discussed area.