基于自动检测方法的福建地区断裂带地震活动性研究.pdf

刘芳，孙冬军，周一剑，朱艾斓，魏薇，朴健.2023.基于自动检测方法的福建地区断裂带地震活动性研究.地震学报，45（3）：538549.doi：10.11939/jass.20220188.LiuF，SunDJ，ZhouYJ，ZhuAL，WeiW，PiaoJ.2023.SeismicitycharacteristicsoffaultzonesinFujianareabasedonautomaticseismicdetectionmethod.Acta Seismologica Sinica，45（3）：538549.doi：10.11939/jass.20220188.基于自动检测方法的福建地区断裂带地震活动性研究*刘芳1，2），孙冬军1，2）周一剑3）朱艾斓1，2）魏薇1，2）朴健4）1）中国上海200062上海市地震局2）中国上海200062上海佘山地球物理国家野外科学观测研究站3）美国加利福尼亚州滨河92521加州大学滨河分校地球科学学院4）中国北京100087北京大学地球与空间科学学院摘要利用福建地震台网 88 个地震台站于 2015 年记录的 10 个月的连续波形数据，基于台网约束下的长短窗特征触发式震相检测的模板事件识别方法和匹配滤波自动检测方法，识别出模板事件 919 个，匹配滤波后重新识别和定位了共计 2243 个地震事件，其中 ML1.0 以上地震1991 次，主要集中在 ML1.01.5 之间，这极大扩充了该震级范围内的福建地震台网目录。重新定位的地震震源深度分布在 020km 之间，多数发生在 10km 内的上地壳。沿长乐诏安断裂带和永安晋江断裂带，地震分布集中且表现出明显的分段活动性。在长乐诏安断裂带上检测到较福建台网目录更多的地震，且地震活动性表现为西南段较强、东北段较弱，沿断裂带的震源深度分布由南向北逐渐变浅，与莫霍面深度的变化相一致；在与之共轭的永安晋江断裂带上则北西段地震活动更为集中，而东南段地震活动稍显分散。重新识别与精定位后的地震与福建地区地壳水平形变南强北弱、东强西弱的特点相一致。关键词 自动地震检测方法长乐诏安断裂带永安晋江断裂带地震活动性doi：10.11939/jass.20220188中图分类号：P315.5文献标识码：ASeismicity characteristics of fault zones in Fujian areabased on automatic seismic detection methodLiuFang1，2），SunDongjun1，2）ZhouYijian3）ZhuAilan1，2）WeiWei1，2）PiaoJian4）1）Shanghai Earthquake Agency，Shanghai 200062，China2）Shanghai Sheshan National Geophysical Observatory，Shanghai 200062，China3）Department of Earth Sciences，University of California，Riverside，California 92521，USA4）School of Earth and Space Sciences，Peking University，Beijing 100871，China*基金项目上海佘山地球物理国家野外科学观测研究站海洋地震观测研究室课题（2021SSY01）资助.收稿日期2022-10-07 收到初稿，2022-12-25 决定采用修改稿.作者简介刘芳，硕士，工程师，主要从事测震、噪音和面波成像方面的研究，e-mail：第45卷第3期地震学报Vol.45，No.32023年5月（538549）ACTASEISMOLOGICASINICAMay，2023Abstract：DrivenbythetectonicloadingofTaiwanandTaiwanStrait，thereareseverallargeQuaternaryfracturezonesintheFujianarea，withstronghistoricalandcurrentseismicity.WebuildacomprehensiveseismiccatalogwithanewautomaticworkflowcalledPALM（phasepicking，association，location，matchedfilter），whichcombinesSTA/LTAandmatchedfilter.PALMisappliedto10monthscontinuousdatain2015of88stationsfromtheFujianEarth-quakeAgency.Werecognized919templateeventsandfinallyrelocated2243eventsaftermatchingfilter，amongwhich1991eventswithmagnitudegreaterthanML1.0，anditgreatlyexpandedthecatalogfromFujianEarthquakeAgencyintherangeofML1.01.5.Thedepthofrelocatedearthquakeeventsaredistributedbetween0and20km，withmostofthemoccurredwithin10km.TheeventsareconcentratedalongChangle-ZhaoanfaultzoneandYongan-Jinjiangfaultzoneandshowedobvioussegmentalactivity.WedetectedmoreearthquakesthanthecatalogfromFujianEarthquakeAgency，andtheseismicitywasstrongerinthesouthwestsectionandweakerinthenortheastsection.Thedepthdistributionalongthefaultzonebecameshallowerfromsouthtonorth，whichwasconsistentwiththechangeofMohodepth.Besides，wecanseeahighlyclusteredfeatureofseismicityatnorthwestsectionandslightlydispersedat southeast segment in Yongan-Jinjiang fault zone，which forms a conjugate shape withChangle-Zhaoanfaultzone.TherepickedandrelocatedeventsinFujianareawasconsistentwiththecharacteristicsofhorizontalcrustaldeformation，whichwasstrongerinsouthandeastwhileweakerinnorthandwest.Key words：automatic seismic detection method；Changle-Zhaoan fault zone；Yongan-Jinjiangfaultzone；seismicity 引言1/2受太平洋板块、菲律宾海板块和欧亚板块的俯冲和碰撞作用，我国台湾及台湾海峡地区的地震活动较为强烈。受台湾及台湾海峡强烈的构造运动影响，福建地区构造活动较为复杂。该地区自西而东存在三条 NE 走向的地震构造带，分别为邵武河源、政和海丰和长乐诏安地震构造带（徐锡伟等，2016；郭晓然等，2019），与之共轭的三条断裂带自南向北分别为上杭云霄、永安晋江和闽江断裂带（李祖宁等，2014）（图 1）。福建地区地质构造较为复杂，历史上曾发生过多次 M7.0 以上地震，例如 1600 年广东南澳附近的 M7.0 地震、1604 年 12 月 29 日福建泉州外海 M7地震、1918 年 2 月 13 日南澳附近 M7.3 地震（国家地震局震害防御司，1995；中国地震局震害防御司，1999）。M5.0 以上破坏性地震主要分布在福建西南地区。福建地区及台湾海峡现今中小地震也较为频繁，仅 2015 年福建地震台网记录到的编目地震就达 974 条（除单台记录外），震级主要集中在 ML0.02.9，地震主要分布在各断裂带上。因此，研究福建地区断裂带的地震活动性，预测未来中强地震可能发生的区域，是防震减灾工作的重要内容。目前针对福建省和台湾海峡地质构造和地震活动性的研究十分广泛，包括：利用天然地震 Pg 和 P 震相走时联合反演地壳上地幔顶部结构（李细兵等，2019a）；震源机制解的应力分析（钟继茂，2009）；气枪激发震源的一维 P 波速度结构（李细兵等，2019b）和三维 P 波层析成像（金震等，2018）等。虽然研究者通过人工激发地震的方式获取研究区域的精细速度结构，但精确的三维速度结构模型在目前的地震定位中应用不够广泛；而人工震相识别受人为因素影响存在系统偏差，且受背景噪声干扰，从而导致地震定位精度不高、地震目录完备性欠3期刘芳等：基于自动检测方法的福建地区断裂带地震活动性研究539缺，不能满足现今断裂活动性及地震活动性研究的需要。因此，高精度地震事件自动检测技术成为目前主流的研究方向。近年来，机器学习和匹配滤波技术在获取较完整地震目录方面取得了很好的成效。随着机器学习的突破性发展，基于数据驱动的机器学习微震检测识别技术极大地扩充了地震目录数量，在国内外地震活动性、断裂精细结构及活动习性、诱发地震等研究中得到了广泛应用（Zhanget al，2019；Zhuet al，2019；Liuet al，2020）。但机器学习需要从大量训练样本中提取统计特征，除对地震样本数有一定要求外，还对计算机性能有较高要求。匹配滤波方法（matchedfiltertechnique）（Turin，1960；Gibbons，Ringdal，2006）是通过波形互相关寻找与已知模板事件相似的地震事件，可在信噪比相对较低的条件下获取更为完备的目录（Shellyet al，2007），但其检测能力仍然受限于模板事件的分布和数量，且运行效率较低。为了规避匹配滤波方法在地震自动检测中的局限性，在计算机配置有限的条件下，我们采用基于台网约束下的长短窗特征触发式震相检测的模板事件识别方法和匹配滤波自动检江西广东浙江三明莆田厦门漳州116 E11812028N 262422F4F5F3F6F7F1F2晚更新世断裂2.0ML3.06.0MS7.0全新世断裂历史强震2015年中小地震0ML2.05.0MS6.0前第四纪断裂福建台网台站分布ML4.0早中更新世断裂3.0ML4.0MS7.0台湾海峡台湾福州图1福建地震台网台站、历史中强地震、2015 年中小地震分布和区域断裂带分布图F1：滨海断裂带；F2：长乐诏安断裂带；F3：政和海丰断裂带；F4：邵武河源断裂带；F5：闽江断裂带；F6：永安晋江断裂带；F7：上杭云霄断裂带，下同Fig.1Thedistributionofseismicstations，historymoderateandstrongearthquakes，mid-smallearthquakesin2015andfaultzonesinFujianareaF1：Binhaifaultzone；F2：Changle-Zhaoanfaultzone；F3：Zhenghe-Haifengfaultzone；F4：Shaowu-Heyuanfaultzone；F5：Minjiangfaultzone；F6：Yongan-Jinjiangfaultzone；F7：Shanghang-Yunxiaofaultzone，thesamebelow540地震学报45卷测方法（phasepicking，association，location，matchedfilter，缩写为 PALM）（Zhouet al，2019）对2015 年福建地震台网连续波形进行扫描和地震精定位，以期获得断裂带的地震活动性。1 方法和数据 1.1 方法PALM 方法是一种基于匹配滤波技术的自动目录扩展方法，利用 STA/LTA（Allen，1978）拾取震相，在台网的约束下，通过震相关联算法去除单台误识别，获得完整的模板事件，再由基于匹配滤波技术的目录扩展方法，利用高性价比 GPU 计算地震模板事件与原始连续波形的互相关函数，提高地震识别率，尤其是在福建地区这类背景噪音相对较高的区域，可有效扩展地震目录数目，为后续地震构造背景研究提供可靠数据。PALM 方法在构建 2021 年漾濞 MS6.4 高精度地震序列时（Zhouet al，2021）与 Liu 等（2022）用类似方法获取的地震条目在同一数量级上，且地震分布具有较好的一致性，该方法在 2019 年 7 月 4 日南加州里奇克莱斯特（Ridgecrest）地震序列（Zhouet al，2022）的构建上亦取得了良好效果。PALM 方法由两部分构成：PAL 模板拾取和 MSMS 匹配滤波。PAL 模板拾取部分分为三个步骤：震相拾取（phasepicking）、台网约束下的震相关联（association）及地震事件的定位（location）。首先震相拾取基于 STA/LTA 特征的触发式震相检测和寻峰式震相拾取算法（Allen，1978，1982），利用STA/LTA 长短窗能量比检测 Z 分量的 P 波峰值，再通过偏振滤波器在水平分量上将 P 波能量进行压制后利用 STA/LTA 正确识别出 S 波峰值，该方法可以准确地给出 P 波和 S 波震相到时。其次，震相关联算法基于台站记录在时间和空间上的关联，地震事件在各台站的记录需满足：在时间域，由各台站的 P 波和 S 波到时估算出的发震时刻与实际发震时刻相近；在空间域，台站的震中距与 P 波和 S 波到时差估算出的震中距相近。时间域的关联通过搜索发震时刻丛集实现，而空间域的关联则通过网格搜索走时残差最小的震中位置实现。最后，利用 HypoInverse（Klein，2002）定位法进行地震事件定位，得到 PAL 模板事件目录。MSMS 匹配滤波部分是一种可在 GPU 上进行计算的自动匹配滤波技术方法。该方法共分为四步：计算事件模板与待检测连续波形的互相关函数；根据模板时间走时将各台站的互相关函数对齐；用蒙版的方式拓展互相关函数峰值的范围，由于模板事件与被检测事件的震中位置差异、模板事件震相拾取的误差以及介质波速变化等因素，仅通过走时偏移无法完全对齐互相关峰值，造成后续叠加时相消，而通过应用峰值蒙版可以减少此类偏差，从而提高模板的检测能力；将各台用蒙版拓展过的互相关函数进行线性叠加（Zhouet al，2019），多台叠加可有效压制噪声的互相关值，从而降低误识别率。基于互相关的震相拾取可以准确地确定相对走时：ttcorrtorgttemp，（1）式中：tcorr为互相关峰值位置确定的震相到时；torg为被检测事件的发震时刻；ttemp为模板事件在该台站的走时。根据模板事件与被检测事件的相对走时，即可使用双差重定位算法（HypoDD）确定事件间的相对位置（Waldhauser，Ellsworth，2000）。基于 PALM 方法（Zhou，2022），我们可以构建从原始波形到获取精定位地震事件的整体架构。在 HypoInverse 绝对定位和 HypoDD 相对定位中，我们均使用一维福建 2015 模型（朱元清等，2017）作为初始速度模型。图 2 分别为以 2015 年 1 月 27 日和 2015 年 4 月 19 日为起始时间的不同台站垂向原始波形记录，这两个波形事件在未滤波前初始 P 波震相较为清晰但后续 S 波震相肉眼较难识别，因3期刘芳等：基于自动检测方法的福建地区断裂带地震活动性研究541此均未被福建台网拾取为地震事件，但利用 PALM 方法可识别为模板事件，并参与后续匹配滤波算法，可获得更多地震事件。1.2 数据本文收集了福建台网 88 个台站 2015 年 14 月和 712 月共 10 个月的连续波形数据。台站分布如图 1 所示，中部和南部台站密度相对较大，北部台站密度较为稀疏。将收集到的连续波形数据进行预处理，并在 240Hz范围内进行带通滤波，得到可进行 PAL 计算的波形数据。由于本试验使用的 GPU 具有显存为 8G 的 M4000 显卡，在 MSMS 自动检测过程中，需要将福建地震台网一整天所有台站的连续波形读入内存中，因此，需要将 88 个台站的连续波形数据降采至每秒 25 个数据才能读入显卡内存中，以便快速进行 GPU 运算。另外，为了减少运算时间和降低噪声、脉冲标定等局部影响以确保地震事件识别的准确性，设定至少有 6 个台站同时接收到信号才被当作一个地震事件。PALM 方法可高效地从连续波形中自动检测出地震事件，配合定位程序 HypoInverse 和 HypoDD 获得高精度震源位置，从而给出构造带的地震活动特性。2 计算结果在进行 STA/LTA 特征触发式震相检测和寻峰式震相拾取算法时，将信噪比阈值设置为12，将长短窗平均能量比阈值设置为 8，S 波窗长设置为 20s，以避免福建地震台网外部及台湾地区的地震被计算。由于噪声、脉冲等非天然事件的随机性，当非天然地震事件被当作模板进行后续的 MSMS 自动检测并与原始连续波形进行互相关时，互相关系数会远小于设置的互相关阈值 0.4，不会检测出新的地震事件，即被淘汰。为了使模板事件波形质量更高，更好地服务于后续模板匹配滤波检测，我们规定保留由最少 6 个台站参与计算的模板事件进行后续计算，在进行 PAL 计算后进行目录整理，最终获得地震事件模板 919 个。图 3 所示蓝色波形为 PAL 方法检测到的 3 个事件垂直分量波形示例，即 3 个模板事件，图 3a 蓝色波形为PAL 方法检测出的 2015 年 1 月 27 日 03:26:59 发生在福建南平市附近的 ML1.5 地震模板，为MQDQ0102030FZRXMQXJMQXJFZLQFQDZPNTKPHSGYXHTDSXPNJJSPHJFZPYF（a）（b）t/s0102030t/s幅值/103 counts幅值/103 counts图2起始时间 2015-01-27T03:26:56（a）和 2015-04-19T04:02:47（b）下的不同台站垂直分量原始波形记录Fig.2Verticalcomponentoriginalwaveformsindifferentstationsstartedin2015-01-27T03:26:55（a）and2015-04-19T04:02:47（b）542地震学报45卷图 2a 原始波形进行预处理、滤波、降采后的波形，可从该图中不同台站垂向分量上清晰地识别出 P 波及 S 波震相，该模板事件可参与后续 MSMS 检测；图 3b 蓝色波形为 PAL 方法检测出的 2015 年 4 月 19 日 04:02:49 发生在福建漳州的 ML1.9 地震模板，为图 2b 原始波形进行预处理、滤波、降采后的波形，该地震可看到清晰的 P 波和 S 波震相，该模板事件可参与后续MQDQPPSPSPSPSPSPPPPPPPPSSSSPPPSSSSSSSSFZRXMQXJGTSKFZLQFQDZPHSGDSXPPHSGZAQDPHJFYXHTDSXPNJJSPHJFZPYFPNTK（a）（b）（c）0102030t/s0102030t/s0102030t/s图3PAL 模板事件垂向分量（蓝色）和检测到的地震事件的垂向分量（红色）（a）模板事件为 2015 年 1 月 27 日 ML1.5 地震，检测到的事件为 2015 年 1 月 12 日 ML1.4 地震；（b）模板事件为2015 年 4 月 3 日 ML1.9 地震，检测到的事件为 2015 年 4 月 3 日 ML1.4 地震；（c）模板事件为2015 年 12 月 6 日误识别事件，S 震相识别错误Fig.3ExamplesforPALtemplateeventsinverticalcomponent（bluelines）anddetectedeventsinverticalcomponent（redlines）（a）Thetemplateeventison27January2015ML1.5anddetectedeventison12January2015ML1.4；（b）Thetemplateeventison3April2015ML1.9andthedetectedeventison3April2015ML1.4；（c）Errortemplateeventon6December2015withinaccurateSphase3期刘芳等：基于自动检测方法的福建地区断裂带地震活动性研究543MSMS 检测；图 3c 为 2015 年 12 月 6 日 04:43:23 台湾花莲海域 ML4.0 地震，从图中能看到清晰的初始 P 波震相和被误识别的 S 波震相，由于 STA/LTA 在 S 震相平均能量比较低，故不参与后续 MSMS 检测。利用上述模板事件对福建台网 10 个月连续波形进行 MSMS 检测，为了降低运算量以及噪音对地震事件识别的影响，我们同样规定必须至少有 6 个台站同时检测到地震事件才参与互相关计算。这个规定虽然减少了脉冲和噪音对地震识别的干扰，但同时也降低了对小震级事件的识别能力，使得最小完备震级 Mc升高。我们设置互相关阈值为 0.4，即模板事件和连续波形进行 MSMS 检测时，只有互相关值超过 0.4 才认为检测出新的地震事件。图 3a 和3b 中红色波形为利用模板检测到的新地震事件，分别为 2015 年 1-1 月 12 日 03:41:39ML1.4 地震，震中为（26.20N，119.04E）和 2015 年 4 月 3 日 03:48:16ML1.4 地震，震中为（24.27N，118.07E），与相应模板事件的相关系数分别为 0.6 和 0.5，而图 3c 未检测到新的地震事件。模板匹配自动检测算法共获得 7 万 1678 个事件，但由于同一个地震事件会被不同的模板检测到而被重复计数，需要对重复地震进行筛选剔除。在进行 HypoDD 精定位的过10vP864203020100116 E118120116 E11812022242628N 22242628N 江西广东莆田厦门三明F7F6F5F4F3F2F1漳州福州浙江深度/km3020100深度/km3020100深度/km1.0ML1.9ML2.00ML1.0vP/(kms1)福建台网目录地震事件PALM 目录地震事件（a）（b）（c）（d）台湾海峡台湾图4基于 PALM 地震目录的地震事件分布与一维福建 2015 速度模型（a）地震事件水平分布图；（b）震源深度分布沿纬度的投影；（c）震源深度分布沿经度的投影；（d）一维福建 2015 模型Fig.4EventsdistributionofPALMcatalogand1-DFujian2015Pwavevelocitymodel（a）HorizontaldistributionoftheeventsfromPALMcatalog；（b）Verticaldistributionoftheeventsalonglatitude；（c）Verticaldistributionoftheeventsalonglongitude；（d）1-DFujian2015Pwavevelocitymodel544地震学报45卷程中，选取一维福建 2015 模型（图 4d）作为一维速度模型，地震与台站的最大距离设置为200km，采用共轭梯度法进行计算，共进行 3 组迭代，每组迭代 8 次，最终定位残差为 63ms，说明本次计算选取的参数较为合理，与此同时将发震时刻在 4s 以内的地震事件认为是同一地震事件，以剔除重复地震事件。HypoDD 精定位完成后，地震条目缩减到 2243 条，地震分布如图 4 所示，主要集中分布在福建省中部和南部，尤其是沿长乐诏安断裂带和永安晋江断裂带分布明显。3 结果分析 3.1 PALM 目录分析由于计算过程中，S 波窗长选择、互相关阈值选择及至少 6 个台站记录到地震事件等限制因素，PALM 目录对台网外部较远的海域即台湾海峡区域的滨海断裂带和广东南澳附近海域（23N23.5N，117E117.5E）的地震识别较少，不能如福建编目目录一样识别到广东南澳海域和台湾地区等地的地震。在政和海丰断裂带与永安晋江断裂带的交界处（25.5N，117.7E）附近的地震事件，PALM 方法也无法识别。另外上杭云霄断裂带东南段及福建广东交界位置识别出的地震较少。PALM 方法在南平市和三明市附近（26.5N0.5，117.5E0.5）识别到福建地震台网目录未识别到的地震 235 次，原始连续波形如图 2a 所示，模板波形如图 3a 所示；在长乐诏安断裂带上靠近漳州市厦门市泉州市莆田市位置识别出了更多地震事件，尤其是漳州厦门附近断裂带上识别到 421 次，原始波形如图 2b 所示，模板波形如图 3b 所示。仙游序列（25.6N，118.7E）附近 在 2015 年处于平静期，但仍有小震级地震发生，PALM 目录检测到了仙游序列大部分 ML1.0 以上地震。在福州盆地及周边区域（26.0N，119.1E）附近，PALM 目录识别出了更多地震，可为今后研究福州盆地构造提供依据。此外，PALM 目录还识别出了福建台网网缘的浙江文成 12 次地震。图 4 除了展示 PALM 目录水平分布外，还展示了 PALM 地震目录和台网编目目录中的地震的震源深度分别在不同经度和纬度上的投影（图 4b 和 4c）。福建台网的地震定位主要利用人机交互系统中嵌入的单纯型和 hypoSAT 方法，这两种方法对震源深度的限制都依赖于含有深度信息的震相，如 Pn 震相（首波）和 PmP 震相（莫霍面反射波），当台站没有接收到这类震相时对震源深度的约束较差。福建地震台网目录的震源深度主要集中在 525km，小于5km 的震源较少，且深度只保留整数位。主要集中在广东南澳海域的地震在沿经度和纬度投影的震源分布图中显示为有一定间隔（间隔为 1km）的短横线。PALM 目录利用 HypoDD 方法进行地震精定位后，可获得精度更高的震源位置，尤其对震源深度有更好的约束。PALM 目录震源深度主要集中在 020km，且大部分地震深度小于 10km，表明地震主要发生在上地壳。在 2243 条 PALM 地震目录中，ML1.0 以上地震为 1991 条，ML1.1 地震数目最多，为346 条，整体呈现高斯分布特征，主要集中在 ML1.01.5（图 5a）。PALM 目录中的震级由各个台站 S 波振幅的中位数获得，由于滤波的原因，其整体震级较台网编目地震震级偏小，存在系统偏差。图 5 同时给出了 PALM 目录和福建台网编目目录不同震级的数目及累计地震数目的统计和对比。由于 PALM 目录的地震事件需至少 6 个台站记录，而福建编目目录只需要不少于 4 个台站记录，因此利用 zmap 软件（Wiemer，2001）中的最大曲率法（Wiemer，Wyss，2000）计算福建编目目录最小完备震级 Mc为 ML0.8，而 PALM 目录最小完备震级 Mc仅为3期刘芳等：基于自动检测方法的福建地区断裂带地震活动性研究545ML1.1。福建编目有更丰富的 ML3.0 以上目录，累计震级-频度曲线斜率更为平缓，而 PALM目录对 ML1.01.5 目录进行了大量补充，而小于 ML1.0 小震和大于 ML2.0 地震数目较少，使得累计震级-频度值增长较快，其斜率较大，平均 b 值增大至 1.5。b 值平均值较大说明PALM 目录不完整，因此结合福建编目目录中小于 ML1.0 和大于 ML2.0 的地震研究福建地区地震活动特性更可靠。3.2 福建地区断裂带地震活动性长乐诏安断裂带是福建地区最活跃的断裂带，呈东北走向，内部发育一系列次级断裂带，以左旋压性活动为主，走滑速率为（0.71.3）mm/a，福建地块呈现出顺时针运动状态（占惠等，2015）。在 PALM 自动检测获得的地震目录中，有 755 次地震发生在长乐诏安断裂带上，占地震总数的 1/3 左右；而福建台网编目目录中长乐诏安断裂带上的地震仅59 次。地震分布以永安晋江断裂带为界，南西段地震条目 571 条，北东段地震条目184 条，呈现南西段较强、北东段相对较弱的地震活动特性（图 4），与福建沿海地壳水平形变南强北弱、东强西弱（陈光等，2010）的特征相一致。为了研究震源深度分布特征，我们沿着断裂走向提取了一条震源深度分布剖面 AB（图 6a），并将剖面附近 20km 范围内的地震进行等深度投影。从图 6b 可以看出，地震震源深度主要集中在 010km 范围内，从 A 点到 B 点震源深度呈逐渐变浅的趋势，与福建沿海莫霍面深度由西南向东北方向变浅的特征（金震，2018）一致。永安晋江断裂带的走向与福建地区主应力方向一致，整体构造活动较长乐诏安断裂带弱，但历史上依然发生过中强地震（国家地震局震害防御司，1995；中国地震局震害防御司，1999）。永安晋江断裂带的左旋走滑速率和倾滑速率分别为（0.61.8）mm/a 和（0.40.5）mm/a（占惠等，2015），其地震动力作用与菲律宾海板块对台湾岛弧作用密切相关，动力作用波前峰大致位于泉州湾，断裂带东南段以拉张走滑运动为主（林松建，黄宗林，2013）。根据 PALM 目录，其整体发生地震较长乐诏安断裂带少，主要分布在断裂带两端，地震数目为 180 条。图 6a 中剖面 CD 为沿着该断裂带从北西至南东提取的地震分布剖面，同样将剖面附近 20km 范围内的地震等深度投影至剖面上（图 6c）。从图 6c 中可以看出，剖面频次MLMc累计震级频度分布非累计震级频度分布00.5（a）（b）1.01.5ML00.51.01.52.02.53.03.54.02.02.53.0103102101100频次103102101100Mc累计震级频度分布非累计震级频度分布图5PALM 地震目录（a）和福建地震台网编目目录（b）震级-频度分布及最小完备震级 Mc对比Fig.5Frequency-magnitudedistributionofPALMcatalog（a）andFujianEarthquakeAgencycatalog（b）andtheircompletemagnitudesMc546地震学报45卷CD 上的地震分布以政和海丰断裂带为界，可分为两段：北西段地震分布更为集中，深度主要集中在 10km 范围内；南东段地震分布较为稀疏，深度主要在 20km 范围内，中部存在地震平静，与其动力环境具有一致性。4 讨论与结论采用 PALM 自动检测方法对福建台网 2015 年 10 个月的连续波形数据进行自动识别与重新定位，共获得 2243 条地震目录。震级主要集中在 ML1.01.5 之间，震源深度主要集中在020km 范围内，且大部分地震发生在上地壳内，深度小于 10km。PALM 目录对福建台网ML1.01.5 地震事件进行了极大的扩充，在台网的西北边缘南平、三明一带识别出台网目录丢失的 235 个地震事件，在福建沿海长乐诏安断裂带多识别出 696 个地震事件，为地震活动性分析提供了较福建台网更为完备、定位精度更高的地震目录。但信噪比阈值和至少 6 个台站接收波形等参数设置仍在一定程度上限制了震级下限及网外海域地震的检测能力。重新识别与精定位后的地震活动沿长乐诏安断裂带与永安晋江断裂带分布集中，且表现出明显的分段活动性：沿长乐诏安断裂带震源深度分布由南向北逐渐变浅，与莫霍面深度的变化相一致；在永安晋江断裂带北西段地震活动更为集中。福建地区地震活动性与水平形变南强北弱、东强西弱的特点相一致。PALM 方法可以在福建地区取得较为可靠的地震目录，为后续依靠地震活动性展开的研究提供了重要依据。审稿专家为本文提出了宝贵的修改意见和建议，福建省地震局提供了 2015 年台网连续波形数据，文中图件使用 GMT（Wessel，Smith，1998）制作生成，作者在此一并表示感谢。116 E11812022242628N 江西广东F7ABDCF6F5F4F3F2F1莆田台湾厦门三明漳州福州浙江（a）0501001502002502015105距离/km深度/kmAB0（b）050100150200距离/km15105深度/km0CD（c）台湾海峡图6剖面 AB 和 CD 位置（a）及沿两剖面的震源深度分布图（b，c）Fig.6LocationofprofileABandprofileCD（a）andeventsdistributionwithdepthalongthesectionsAB（b）andCD（c）3期刘芳等：基于自动检测方法的福建地区断裂带地震活动性研究547参 考 文 献陈光，李祖宁，王紫燕，钟继茂，陈超贤.2010.福建地壳运动与地震关系初探 J.华北地震科学，28（1）：2125.ChenG，LiZN，WangZY，ZhongJM，ChenCX.2010.StudyontherelationshipbetweencrustalmovementandseismicityinFujianareaJ.North China Earthquake Sciences，28（1）：2125（inChinese）.国家地震局震害防御司.1995.中国历史强震目录（公元前 23 世纪公元 1911 年）M.北京：地震出版社：141144，453455.DepartmentofEarthquakeDisasterPrevention，StateSeismologyBureau.1995.TheCatalogueofChineseHistoricalStrongEarthquakes（The23rdCentury BC1911 AD）M.Beijing：SeismologicalPress：141144，453455（inChinese）.郭晓然，谢志招，闫培，李普春，王笋.2019.福建台湾海峡深部地壳结构探测与滨海断裂带性质研究 J.华南地震，39（2）：3442.GuoXR，XieZZ，YanP，LiPC，WangS.2019.TheexplorationonthedeepcrustalstructureintheFujian-WesternTaiwanStraitandstudyonthelittoralfaultzoneJ.South China Journal of Seismology，39（2）：3442（inChinese）.金震.2018.福建及台湾海峡中南部三维地壳结构研究 D.哈尔滨：中国地震局工程力学研究所：8384，8789.JinZ.2018.StudyonThree-DimensionalCrustalStructureintheCentralandSouthernFujianStraitandTaiwanStraitD.Harbin：InstituteofEngineeringMechanics，ChinaEarthquakeAdministration：8384，8789（inChinese）.金震，李山有，蔡辉腾，李培，李海艳，徐嘉隽.2018.利用气枪地震资料对福建及台湾海峡南部地壳三维 P 波速度结构研究 J.地球物理学报，61（7）：27762787.JinZ，LiSY，CaiHT，LiP，LiHY，XuJJ.2018.3DP-wav