山东乳山震群遗漏地震检测及震源区活动特征分析
2021-07-21王鹏郑建常
王鹏 郑建常
1)山东省地震局，济南 250102 2)东营市地震监测中心，山东东营 257000
0 引言
2013年以来华北地区震群活动频发，其中山东乳山震群尤为突出。2013年10月1日乳山地区地震丛集活动持续至2019年，共记录到小震1.3万余次，其中ML≥3.0地震31次，最大震级为2015年5月22日ML5.0，期间伴随着多次地震活动增强—减弱的过程。
乳山震群所在的胶东半岛地区，位于华北地块东缘与苏鲁造山带山东段的交汇位置，在大地构造划分上属于华北克拉通与扬子克拉通的结合部位(晁洪太等，2001；周瑶琪等，2015)。受苏鲁造山带与古太平洋板块俯冲的双重影响，该区域岩浆活动十分频繁，自太古代到新生代均有发现岩浆活动，其中以燕山晚期岩浆活动最为剧烈，同期还发生了大规模的岩浆侵入作用。在火山活动区域，高热流区会导致较小的低粘度累积弹性应变，相对于主余型序列而言，往往会产生更多的震群(Ben-Zion et al，2006；Yang et al，2009)。根据历史地震记录，胶东半岛地区历史上大震较少，却经常发生震群活动(李冬梅等，2011)，但这些小震序列(震群)持续时间通常很短，一般持续几天，最多不会超过1个月。而乳山震群持续6年，且频度和强度均较大，这在大陆东部地区甚为罕见。对乳山震群震源区活动特征的深入研究可为区域地震危险性分析和地震趋势判断提供有效的依据。
火山或地热地区的震群发震机制通常可以用流体运移模型来解释，在地壳中运移的流体包括水和挥发性气体(CO2)等，其通常被认为是由中地壳到上地幔这一深度的岩浆源提供的(Takahata et al，2003；Ohba et al，2011)。在前人对乳山震群发震机制的研究中，也猜测流体可能是震群的触发因素(Zheng et al，2017；郑建常等，2019)，但始终缺乏直接的观测证据，而在震群发震构造(曲均浩等，2014、2015、2019；郑建常等，2015a、2015b；李铂等，2017；张斌等，2017)、震源参数(苗庆杰等，2016；王鹏等，2016；郑建常等，2016；刘方斌等，2018)等问题的讨论上也存在不同见解。究其原因为，乳山地区断裂错综复杂，次级断裂(近SN向)、隐伏断裂(NE、NW向)较为发育，断裂格局总体呈现羽状或雁列式形态(刘善宝，2005)，震中区并无活动断层，距离最近的乳山断裂15km。
统计学方法可以从复杂的现象中探究本质规律，自Ogata(1988)在“大森-宇津”公式(Omori，1894；Utsu，1961)的基础上提出了“传染型余震序列”(Epidemic-type Aftershock Sequence，ETAS)模型后，该模型已受到广泛的关注，其在研究地震活动变化规律(Ogata，1992、2001；Zhuang，2000；蒋海昆等，2007；蒋长胜等，2013a、2013b、2013c、2015、2017)、探测应力变化(雷兴林等，2008；Lei et al，2017)、识别与流体作用有关的诱发地震(Hainzl et al，2005；Lei et al，2008；龙锋等，2010；蒋海昆等，2011、2012)等方面得到了长足发展。
统计地震学方法很大程度上依赖于地震目录，地震目录的不完备势必会影响到基于地震目录的相关统计学结果的可靠性(李智超等，2014)。地震目录的完备性主要依赖于研究区域地震台网的监测能力，而乳山震群周围台站分布格局不够理想，地震序列短时间的丛集发生通常会造成地震的遗漏。通过遗漏地震检测来使地震目录更加完整(Schaff，2008)，可以提高分析结论的可靠性。用于微震检测的模板匹配滤波技术(Gibbons et al，2006)提出后，该技术发展迅速并得到了广泛的应用，在检测地脉动(Shelly et al，2007)、研究余震序列变化(Peng et al，2009；侯金欣等，2017；王鹏等，2016)、确定发震构造(Yang et al，2009；谭毅培等，2016)、分析震源区应力变化(Meng et al，2013)、识别重复地震(马腾飞，2015)、监视核爆(Zhang et al，2015)等诸多方面发挥着重要的作用。
本文利用基于GPU并行计算的模板匹配方法检测乳山震群的遗漏地震，基于完备地震目录对该震群ETAS模型的统计参数进行分析，并初步推测乳山震群的发震机制。
1 理论方法
本文所用模板匹配方法的处理流程参照Peng等(2009)，并与侯金欣等(2017)的处理步骤类似，具体方法和处理流程可见Wang等(2019)，在此不再赘述，仅对ETAS模型做简单介绍。
ETAS模型条件强度函数λ(t)(Ogata，1988；蒋海昆等，2012)为
(1)
式中，右侧第二项为第i次地震对地震发生率的贡献，i遍历地震序列中的所有地震，Mz为参考震级，其必须大于等于最小完整性震级Mc；p代表地震序列衰减的快慢，p大则衰减快，反之则衰减慢；α表示触发次级余震的能力，α越大表示越强的高阶余震的激发能力；μ(t)最初的物理含义为剔除地震丛集后单位时间内的理论地震数目，也表示由于外力作用触发微震活动的强弱，与地震间的自激发过程无关(Hainzl等，2005)；c为数值极小的正值，其保证了右侧第二项的分母不为0，与主震后短时间内数据的完备程度(蒋海昆等，2007)和震源区的应力状态有关(Narteau et al，2009)；K为与主震震级及b、p、α、c、μ参数有关的地震活动水平，其中b为G-R关系中的比例系数(宋金等，2009)。
对于同一组观测数据，数学上可采用多个统计模型对其进行拟合，并可得到多组符合条件的参数。Ogata(1988)利用最大似然法估计ETAS模型参数，并利用最小化模型的信息量准则(Akaike information criterion，AIC)进行最佳模型参数的判别(Akaike，1974)。AIC准则既考虑模型对观测数据拟合的优劣，也会惩罚为了提高拟合程度而无限制地增加模型参数的行为。ETAS模型的对数似然值为
(2)
AIC=-lgL+2k
(3)
其中，tb和te为用于AIC计算的起止时间；NAIC为计算时段内的地震数目，为了保证计算结果的可靠性，计算时段内的地震目录必须保证基本完备；k为待拟合参数个数，本文中k=5，取AIC最小时对应的模型参数作为最佳模型。
由式(1)可知，ETAS模型所表征的地震活动率包括两部分，其中右侧第二项是由于“余震激发余震”所导致的地震丛集，右侧第一项是与余震激发无关、因外力作用而产生的地震活动。因而，通过对ETAS模型两部分地震活动进行分离，将外力触发地震比例Rb用式(4)表示，而自激发地震所占比例Rt=1-Rb，可以根据各自所占的比重来讨论外因触发与地震活动的关系(蒋海昆等，2011)
(4)
2 数据处理
乳山震群周围100km内共有9个固定地震台站，均位于震群的北部、西部和东部，而南侧无固定台站。乳山流动台阵于2014年5月架设完成，共计18套宽频带地震仪，采样率为100Hz(图1(b))，平均台间距约为0.5km，距离震群最大距离约10km。乳山台阵对震群的包络性较好，在震群南侧架设的宫家岛台(GJDT)弥补了震群南侧海岛无固定台站的缺陷。乳山台阵自2014年5月投入使用，乳山震群持续时间较长，震群最大地震发生在2015年5月，地震活动主要集中在2015年6月，后续地震活动开始逐渐衰减。因此，选择2014年5月至2015年6月期间固定地震台站和流动台阵的连续波形进行遗漏地震检测。
图 1 乳山震群和地震台站空间分布(a)中黑色三角形代表固定台站，蓝色三角形代表乳山台阵，灰色圆圈为台网定位结果，红色圆圈为重定位结果，红色方框为研究区域，红色五角星为乳山震群；(b)为图(a)中黑色方框区域放大
首先，对连续波形数据去均值、去线性趋势，并采用零相移4阶Butterworth带通滤波器对2～10Hz进行带通滤波。选取台网目录(中国地震台网中心)中ML≥2.0的地震作为模板地震，根据台网目录中发震时刻截取乳山台阵中距离乳山震群最近且分布较好的RSH、GJDT、DFS、NHNC、WHUT和XJZT等6个台站的三分量记录作为模板地震波形，再依据震相报告读取模板地震的P和S波到时，计算并选取信噪比大于5的109个地震作为模板地震，最大震级为ML4.5。然后，按照模板匹配方法的数据处理流程(Wang et al，2019)，利用模板地震与连续波形互相关计算后做叠加得到互相关系数波形，为了减少误检测事件，本文选取12倍的绝对离差中位数(Median Absolute Deviation，MAD)作为判定新检测地震的阈值，并设定5s时间范围内仅检测1次地震事件，选取互相关系数最大的遗漏事件，最后确定检测事件的震级和发震时刻。
图 2为检测到的一次ML2.7遗漏地震的示例。其中，图2(a)为互相关系数和MAD的倍数，可见在1s的时间内互相关系数出现多次高值，且均高于12倍MAD的阈值，但根据假定只取互相关系数最大的点，即图中红色圆圈所对应的时间；图2(b)表示了遗漏地震所在的波形序列的位置，可见在遗漏地震之前连续发生了2次较大地震，分别是2014年9月16日14时42分59秒的ML4.1和14时43分33秒的ML3.9地震，而该遗漏地震几乎淹没在前面较大地震的尾波里，幅值是前面较大地震的1/10，不易分辨，这也是地震漏检测的原因；图2(c)、2(d)为RSH台和DFS台记录的模板地震在连续波形上的扫描结果，可以看到模板地震与检测出的遗漏事件相关性较好，相关系数为0.61。同时，模板地震的S波的幅值与遗漏事件相当，因此根据S波振幅比计算的遗漏地震的震级也为ML2.7。依据检测事件与所对应的模板事件具有相同走时的假设，确定检测地震事件的发震时刻为2014年9月16日14时44分03秒。
图 2 模板匹配方法检测遗漏地震示例(以ML2.7地震为例)(a)为经过偏移、叠加和取平均值后的互相关系数和MAD的倍数图，其中相关系数大小用竖线表示，MAD的倍数用橙色点表示；(b)中红色点线框为遗漏地震的位置，即图(a)所示的时间范围；(c)、(d)为DFS台和RSH台记录的模板地震和遗漏地震事件的波形对比，其中红色为2014年11月19日3时24分30秒的ML2.7模板地震，蓝色为该时段内的连续波形，检测地震事件为2014年9月16日14时44分3秒的ML2.7地震
3 结果分析
3.1 遗漏地震检测结果
2014年5月5日至2015年6月30日期间，共检测到7500个地震事件，其中109个为模板事件的自检测，7391个为检测地震事件，检测地震事件也包括了台网已记录但未用作模板的地震事件。山东地震台网记录到地震1937个，检测地震事件是台网目录的近4倍，由模板匹配方法获得的地震目录简称为检测目录。根据检测目录和台网目录的频度对比(图3)可见，与台网目录相比，检测目录中的遗漏地震主要集中在ML≤1.0的震级范围，1.0
图 3 检测目录和台网目录频度和累积频度对比
利用最大曲率法对比分析了台网目录和检测目录的完备震级，该方法选取震级-频率曲线(Frequency Magnitude Curve)中斜率最大值所对应的震级作为Mc(Wiemer et al，2000)。在实际应用中，Mc震级通常对应非累积频度-震级分布(Non Cumulative Frequency Magnitude Distribution)中具有最大地震频度的震级(黄亦磊等，2016)。
由图3 可以看出，在非累积频度-震级分布中，检测目录和台网目录的最大值对应的完备震级分别为0.1级和0.9级，由于最大曲率方法得到的完备震级通常偏小，利用最大似然法进行拟合，发现检测目录和台网目录分别在震级0.2和1.0级以上b值的拟合较为稳定，最小完整性震级分别取0.2级和1.0级，计算的b值分别为 0.95±0.01和 0.89±0.03。基于模板匹配方法的遗漏地震检测将研究时段内乳山震群的完备震级从1.0级降至0.2级，弥补了台网目录中1.0级以下地震的缺失，使地震目录更加完整。
3.2 利用ETAS模型分析地震活动特征
乳山震群自2013年10月1日开始，期间经历了多个增强—减弱的阶段，根据日频度和累积频度随时间的演化特征(图4)，将整个震群发展过程分为5个活动时段：第一阶段为2013年10月1日至2014年1月6日，2013年10月1日12时7分55秒震群初始地震ML3.8地震发生后，地震活动短时间内集中发生，日频次达到高值，并在接下来的2个月内逐渐减弱至2014年1月6日；第二阶段开始于2014年1月7日，乳山震群发生了ML4.7地震，随后地震活动再一次增强，日活动频度逐渐减弱后又出现了短时间的增强，至2014年4月3日又趋于减弱；第三阶段则开始于2014年4月4日的ML4.6地震，与前一阶段类似，地震活动在短时间的增强后迅速减弱；第四阶段为震群最大地震2015年5月22日ML4.6地震发生后，日频度达到最高值，随后地震活动经过半年多时间的调整，期间也伴随着增强—减弱的过程，至2015年底基本趋于稳定；第五阶段为2016年以后，地震活动较弱，总体较为稳定，期间仅有个别地震活动小幅度增强的现象。
图 4 乳山震群日频度和累积频度随时间的演化左侧纵坐标表示日频度，用黑色点线表示；右侧纵坐标表示累积频度，用红色线表示；蓝色点线将整个发展过程分为5个阶段，黑色点线又将每个阶段分为若干个子阶段
图 5 乳山震群震级、频度-时间图(a)为检测目录(黑点)和台网目录(灰点)的震级随时间变化图，其中红色点线框为本文利用模板匹配方法检测的时间段； (b)为图(a)红色框的放大，并加入累积频度的统计，其中蓝点为检测目录，黑点为台网目录，红色线为累积频度曲线
根据震群发展过程，每个阶段又分为若干个小的阶段(图4)。将检测得到的完备震级加入到台网目录中，利用ETAS模型对这些不同的阶段分别进行反演，得到各模型参数的变化特征，从而得出外力作用和自激发作用对震群发展的影响。
不同的截止震级对ETAS模型参数的反演影响较大(蒋海昆等，2012；蒋长胜等，2013c)，要保证每个阶段的ETAS模型参数反演中截止震级必须大于完备震级。因此首先利用震级-序号法考察完备震级在震群发展过程中的变化(图6)，震级-序号法是按地震发生的次序排列，通过分析不同震级的地震概率密度分布来定性分析完备震级Mc的方法，其中地震概率密度最大的位置所对应的震级即为完备震级Mc。之所以使用地震发生次序，而不是地震时间，主要是为了避免余震或震群中“丛集事件”的影响(蒋长胜等，2013a)。这里以500个地震作为窗长，以0.1级为间隔，计算窗内各个震级地震的概率密度，可以看到在震群发生初期完备震级在0.9级左右，在检测目录加入后使得完备震级明显降低，在检测的2014年5月至2015年6月期间，完备震级基本维持在0.1级左右，在2015年6月之后完备震级又上升到0.9级左右。为了分析乳山震群整个发展过程中ETAS模型参数的趋势变化，并能更好地与其他阶段的参数变化进行对比，因此统一将截止震级取为0.9级，按照图4 中的阶段划分，反演了各个阶段的ETAS模型参数，见表1，其中第一阶段(Ⅰ)至第四阶段(Ⅳ)以及每个阶段的子阶段(Ⅲ1、Ⅲ2，Ⅳ1、Ⅳ2)的拟合图见图7。
图 6 乳山震群震级-序号图色标表示计算窗内地震在某个震级段的概率密度分布
表 1乳山震群不同阶段的ETAS模型参数
图 7 乳山震群各阶段ETAS模型拟合结果
对于第一阶段，背景地震或外力作用的比例为35.5%，说明有35.5%的地震是由外力作用触发的，剩余64.5%的地震由余震的自激发产生。α值为0.947，低于典型的构造地震(2.0)，震群型地震的α值通常小于1(Ogata，2001)。日本地区震群的α值一般介于0.35～0.85之间，而一般的主余序列则在1.2～3.1之间(Ogata，1992)。α值表征余震的激发能力，这一阶段的α值接近1，表明有较强的余震的激发能力。p值为余震衰减率，这里p值(0.99)与典型构造地震(1.0左右)相当，表明序列属于正常衰减。同时，在这一阶段模型理论值与观测值拟合较好，表明参数估计较为准确，外力作用在震群发生初期占一定比重，但主要还是受余震的自激发作用影响。
第二阶段的地震活动伴随着较大地震(2014年1月7日ML4.7地震)短时间内集中发生，随后逐渐衰减，期间发生过1次日频度突然增强的现象。这一阶段外力作用所占比例有所下降，从第一阶段的35.5%降至24.3%，表明外力作用较弱，p值无变化，仍为0.99，表明序列仍处于正常衰减，但α值较高(1.258)，且高于第一阶段，表明余震的激发能力增强。
第三阶段与第二阶段类似，伴随着较大地震(2014年4月4日ML4.6地震)的发生，地震活动出现丛集，随后趋于衰减，由于地震活动速率变化较大，将该阶段分为2个子阶段(Ⅲ1、Ⅲ2，图4)，其中Ⅲ1阶段外力作用所占比例有所上升(39.9%)，表明外力作用有增强的趋势，这种短期丛集的地震受外力作用触发的较多。α值明显下降，由1.258降至0.362，表明余震的激发能力明显减弱，与外力作用增强一致。同时p值也明显增大，变为1.5，表明序列的衰减作用增强。在这一阶段，模型理论值与观测值拟合较差，实际地震发生率(图7(c)蓝色实线)明显高于拟合的理论曲线(图7(c)灰色点线)，表明大森定律型自触发只有很小的促进作用，外力作用引起的地震活动的比例可能被低估(Lei et al，2008)。在Ⅲ2阶段，外力作用所占比例继续增大(45.9)，说明将近一半的地震是由于外力作用触发的，同时α值明显升高(1.444)，表明这一活动时期地震活动的自激发作用也相应增强，p值也再次回到1.01，序列呈现天然地震的正常衰减速率。
第四阶段地震活动速率明显增加，也伴随着震群最大地震(2015年5月22日ML5.0地震)的发生，日频次达到最高值，随后快速下降。将这一活动时段也分为2个子阶段(Ⅳ1、Ⅳ2，图4)，其中Ⅳ1阶段外力作用较上一阶段有所降低(35.2%)，约有 1/3 的地震是由外力因素触发的，而α值依然保持较高值(1.002)，表明这一阶段更多的是受最大地震(ML5.0)的影响，激发了短期内的余震发生，p值为1.06，表明仍在正常衰减。Ⅳ2阶段在衰减(p值)变化不大的情况下，α值略微下降(0.969)，说明受最大地震(ML5.0)的影响减弱，同时受外力作用明显增强(43.6%)，与Ⅲ1阶段类似，该阶段的实际地震发生率(图7(f)蓝色实线)也明显高于拟合的理论曲线(图7(f)灰色点线)，表明外力作用也可能被低估(Lei et al，2008)。外力作用在这一阶段可能起主要作用，而地震的自激发作用次之。
第五阶段随着地震活动的减弱，外力作用比例减弱，α值和p值也有所下降，说明这一阶段外力作用和自激发作用均在减弱，同时衰减也相应减弱，与这一时期持续较长时间的弱地震活动水平有关。
我们还考察了乳山震群的ETAS模型参数Rb、α和p随时间的演化规律(图8)，可以看到除个别点外，Rb、α的变化趋势较为一致，总体呈现先上升后下降的趋势，而p值除在个别点外，其他时段变化较小，基本维持在1.0左右，表明乳山震群一直保持着天然地震正常的平均衰减速率。Rb值在整个过程中基本维持在20%～50%之间，表明外力作用在震群的发展过程起到一定的作用，其在某些阶段甚至接近50%，但总体来看仍是震群前期地震的自激发起主导作用。
图 8 ETAS模型参数随时间的变化
4 讨论
4.1 不同截止震级对ETAS模型反演结果的影响
通过对2014年5月至2015年6月期间的遗漏地震进行检测，使得完备震级大幅减小，在完备震级随时间的演化中也能明显看到这种特征。2014年5月至2015年6月这一时间段对应Ⅲ2、Ⅳ1阶段(图4)，因此，测试这2个阶段在截止震级分别取0.1级和0.9级时ETAS模型的拟合结果，比较不同参数的结果(图9)。从图9 可以看到，外力触发地震的比例Rb差别较大，截止震级取0.1级时2个阶段的Rb值明显低于截止震级取0.9级的情况。由于任何一个地震均能以一定概率触发后续地震，低于Mc的大量微震事件将可能显著改变背景地震概率(Rb值)(蒋长胜等，2013b)。在III2阶段，随着截止震级的减小，α值和μ值明显减小，表明序列触发微震的能力和外力作用触发微震活动的能力均有所减弱。在Ⅳ1阶段，随着截止震级的减小，α值和μ值则有所增加，表明自激发作用和外力作用均有所增强。这种现象可能与2个阶段的地震活动特征不同有关，III2阶段地震活动趋势较为平稳，而Ⅳ1阶段由于震群最大地震的发生，地震活动频度迅速增加，日频度达到最大值(图5)，震级越大，激发高阶余震的能力越强，对后续地震的影响也就越大。尽管α值和μ值均有所增加，受震群最大地震的影响，外力触发地震的比例却很低(3.7%)，地震活动中仍是自激发作用的贡献更大。另一方面由于震群中最大地震发生后，短期内ETAS模型参数随时间变化较为明显且波动较大，ETAS模型参数与地震序列持续时间也有一定关系(蒋长胜等，2013c)，2个阶段的持续时间相差较大，可能也是一方面因素。此外，震源区的破裂特征、断层的愈合速度、应力加载、周边构造场的调整等诸多因素均会使地震序列的衰减特征和激发余震能力出现明显差异。
由于ETAS模型假定每一次地震均可触发其后发生的任何震级的地震，因此，截止震级在确保地震目录完整性的同时，也可能切断了截止震级以下地震与以上地震的触发关系，从而造成“链接缺失”现象(Wang et al，2010)。通过遗漏地震检测，使得完备震级减小，而截止震级的减小尽可能地减少了这种现象的出现，因此地震目录的完整性对ETAS模型参数反演至关重要。
经过遗漏地震检测，在保证地震目录完整性的基础上，降低了截止震级，计算得到的ETAS模型参数也更加可靠。从理论曲线和实际地震活动的拟合来看，截止震级取0.1级时拟合较好，同时AIC值也更低(图9)。但由于乳山台阵于2014年5月才架设运行，而乳山震群2013年10月开始集中活动，因此基于台阵波形数据的遗漏地震检测的时间段仅是震群发展过程中的一部分(图5)。为了考察乳山震群整个发展过程中ETAS模型参数的趋势变化(图8)，同时能更好地与其他阶段进行对比(图7)，对于整个震群发展过程取截止震级0.9级以上的完备地震目录进行反演计算。将截止震级设为高于完备震级(0.9级)的某个震级，利用ETAS模型进行反演，尽管各参数的数值不完全相同，但反映出各参数随时间的变化规律是一致的。
图 9 不同截止震级下ETAS模型拟合结果对比(a)、(b)截止震级为0.1级时2个阶段的拟合结果；(c)、(d)截止震级为0.9级时2个阶段的拟合结果
4.2 发震机制讨论
本文利用外力作用引起的地震与自激发地震的比例系数Rb来考察外力因素对震群活动的影响。除个别阶段外，Rb和α的总体变化趋势较为一致，呈现出先升高后降低的过程，而p值除个别阶段外，总体变化较小，维持在天然地震正常的衰减速率。Rb值在整个过程中基本维持在20%～50%之间，表明外力作用在震群的发生、发展过程中均占不小的比例，尤其在震群初期或较大地震发生之前，地震可能是受外力因素触发或间接引起的应力调整引发的，但从整个过程的比例来看，余震的自激发作用所占比重更大，发挥着更大的作用。震群发展后期(V阶段)外力作用和自激发作用均在减弱，地震活动呈现较长时间的弱活动水平。
前人在对乳山震群空间分布、震中迁移速率、ETAS模型、高频衰减系数等方面的研究中推测乳山震群可能是由流体入侵触发的(郑建常等，2015b、2016；Zheng et al，2017)，但目前仍没有确切的观测证据。通过分析ETAS模型参数特征和其随时间演化的规律，也可以看到外力因素在乳山震群的发生、发展过程中起着一定的作用，而流体入侵或运移也可能是外力因素之一，今后还需结合其他依据做详细的论证。
5 结论
本文利用基于GPU并行运算的模板匹配滤波技术，对于乳山震群2014年5月5日至2015年6月30日期间的连续地震波形进行检测，识别出7500个地震事件，其中109个为模板事件的自检测，7391个为检测地震事件。检测事件是台网目录的4倍，利用较为准确的b值拟合方法得到的震群完备震级从1.0级降至0.2级。
乳山震群2013年10月开始集中活动，而乳山台阵自2014年5月5日才开始运行，因此，将检测目录补充到整个乳山震群的地震目录中，并划分整个震群发展过程为5个活动时段，借助ETAS模型分析了乳山震群的地震活动性以及在各个阶段的变化特征。结果表明，乳山震群外力触发的地震所占比例较大，结合前人研究，初步推测流体作用可能是外力因素之一，且不同时期外力对微震活动的触发强度不同。
