印度喀拉拉邦伊都基水库附近小震的精确定位和震源机制:地震成因的启示

印度喀拉拉邦伊都基水库附近小震的精确定位和震源机制：地
震成因的启示
Utpal Saikia;S.S.Rai;M.Subrahmanyam;Satyajit Dutta;Somasish
Bose;Kajaljyoti Borah;Rishikesh Meena;任家琪
【摘要】应用印度南部新建的一个拥有21个台站的临时台阵的地震波形,显著提高了喀拉拉邦伊都基水库附近地区小震的检测阀值和参数.我们给出了这一地区16个地震的精确定位,地方震级为1.5～3.6,震源深度7.2～9.9km.所选6个最好地震的断层面解显示走滑断层错动兼右旋运动.水库加载通常会导致应力的产生和在浅层深度(＜5km)引发地震,但是在伊都基水库区域内却没有该深度的地震.记录到的地震仅限于靠近Karur-Kambam-Painavu-Trichur (KKPT)剪切带西北—东南走向的断层.这些结果表明,伊都基水库区域内的地震均属于天然构造地震,而与水库没有联系.
【期刊名称】《世界地震译丛》
【年(卷),期】2016(047)003
【总页数】11页(P177-187)
【关键词】地震定位;断层面解;水库;构造地质
【作者】Utpal Saikia;S.S.Rai;M.Subrahmanyam;Satyajit Dutta;Somasish Bose;Kajaljyoti Borah;Rishikesh Meena;任家琪
【作者单位】;;;;;;;宁夏地震局
【正文语种】中文
作为印度深部地球成像实验(INDEX)的一部分1，我们于2011年1月至2012年
3月运行了一个由21个宽频带地震仪组成的台阵，台阵位于印度南部的泰米尔纳
德邦和喀拉拉邦，用于绘制地震图像和区域深层结构(图1)。

在这个低应变(10-12～10-13/年)地区2，我们的仪器在横跨喀拉拉邦贝里亚尔河的科钦以东约
120km的伊都基水库(9°50′N；76°58.5′E)附近记录了一些小地震。

伊都基水库附近最近发生的地震活动Rajendran等3已报道过。

然而，由于仪器数据不足，他
们不能提供精确的震中位置、震源深度和震级。

伊都基是亚洲最高的拱坝之一(169m)，自1975年10月开始启用。

在历史上，喀拉拉邦地区发生过几次小到中等强度的地震4，5，包括(M<5)的1988年伊都基地震6，7，1994年Wadakachery地震8，2000年和2001年的Erattupeta/Pala双震9，10。

伊都基区域的低增益地震台网从1971年开始运行，已经记录了大量的地震7。

1977～1983年在伊都基水库附近发生的地震(M≤3.5)，被认为是水库引起的11。

另一种观点认为，这些地震是被该地区早就存在的西北—东南走向的断裂所再次
激活的12，13。

由于多数地震都很小，波形质量很差且波形数量不足，因此地震参数无法精确测定。

在这种情况下，我们给出台网运行期间所发生地震的精确震源位置和震源机制解，并且讨论这些地震与水库附近地表地质特征之间的可能联系。

伊都基水库位于西高止山脉的南部，走向与印度西海岸平行。

图1显示了印度南
部的主要区域结构，地质上称为南麻粒岩(SGT)。

地形被良好映射的北部Moyar
剪切带和珀瓦尼剪切带(MSZ/BSZ)和南部的诺伊尔—高韦里剪切带(NKSZ)及KKPT剪切带所分割，印度南部的诺伊尔—高韦里剪切带南边是泛非地形，包括马杜赖北部地块(MB)和南部喀拉拉邦孔兹岩带(KKB)。

伊杜基区的主要岩石类型为紫苏花岗岩、片麻岩和孔兹岩系。

主要构造特征是西北—东南走向的贝里亚尔河轮
廓和KKPT剪切带14，15(图2)。

这个区域的这两个突出的轮廓有右旋走滑分量(右旋)。

2011年1月至2012年3月，我们在该区域部署了一个拥有21个数字宽带地震
台站的临时台阵(图1)。

台站参数在表1中列出。

这些台站都安装了Guralp
CMG-3T/ESP宽带传感器，配备4GB插拔硬盘和全球定位系统的24位REFTEK(RT130/01)数字记录器。

所有的台站都是连续运行模式，波形记录的速度为50个样本/秒。

为了探测和提取当地的地震，我们利用布设的RTVIEW中的REFTEK应用软件手
动检查该区域的连续波形，要求最少有3个代表性的台站来准备当地的地震时间
序列。

此列表用来为所有台站选取地震时间序列。

最后，我们做出了高信噪比
(S/N)和P波震相清晰的地震波形的事件台站列表。

这些地震用来校正仪器响应。

我们使用地震分析软件SEISAN16进行了P波和S波震相的拾取。

根据清楚的震相，我们指定一个不确定的时间范围，P波到时为0.05～0.50s，S波到时为
0.1～1.5s。

为了避免震相影响，到时都是手动从未经过滤的地震记录中拾取的。

P 波在垂直分量拾取，S波在水平分量。

为了给数据质量一个概念，图3呈现了台阵记录的一个地震的垂直分量。

地震(表2)最少挑选12个P波和10个S波的到时，使用HYPOCENTER17程序
来确定震源位置。

震源测定精度取决于记录台站的分布以及震源和台站之间的速度结构。

用于震源定位的初始地壳速度模型(表3)，是一个最初反演结果的简化版本，来自台网台站的接收函数研究(Ritima Das，私人通讯)。

波速比VP/VS的平均值1.71±0.001 7是通过P波的到时和S-P波到时差从和达曲线中得出的(图4)。

这
些地震的方位角范围为120°～125°。

此16个地震的时间残差均方根、震源深度
和震中详细误差统计(H)、震源深度(Z)和时间残差均方根(s)都在图5中显示。

震中误差、震源深度和时间变化范围分别为为0.8～2.2km，1.5～2.5km和0.10～
0.20s。

这些地震的震源深度很好地限制在7.2和9.9km之间。

这些地震的细节列于表2，震中在图2中显示(插图)。

要评估这些地震的深度可靠性，我们选定了2
个高质量的地震，每个地震至少被有清晰的P和S波震相12个台站记录到。

在地震资料分析中，震源深度和发震时刻作为一个耦合参数，这会导致它们根据到时数据估计不准确。

因此，我们进行了一系列的只有3个自由参数的地震反演(纬度、经度和发震时刻)同时保持每一次深度都固定不变(图6)。

这确保了深度和发震时刻是互相解耦的，可以提高震源深度的全局检索。

我们计算地震位置时，震源深度范围1～25km，其中增量变化为1km。

将时间残差均方根最小时的解，作为最后选定的解。

两次地震的结果显示在图6，其中最小均方根的变化值是根据不同震相数据子集的震源深度(组合的P和S波以及单独的P波)标绘的。

针对这两次地震，我们粗略地估算了震源深度，均方根和深度函数都围绕最小值呈陡峭状态。

最低配达到深度10km和9km。

相应的反演结果，震源深度为9.7km和9.1km。

这种敏感性试验证实计算震源深度是可靠的。

地震的里氏震级是由伍德—安德森地震仪记录，根据波的振幅(A)的对数(log)确定的，利用公式ML=log(A/A0)计算，其中A0是在参考距离的振幅。

相应地，震级的调整包括弥补各种地震仪和地震的震中之间距离的变化。

震级由水平分量的最大振幅值确定。

为了获得地方震级(ML)，对地震图进行数字滤波来模拟伍德—安德森地震仪，用记录的最大水平振幅直接确定ML。

最终的震级大小等于所有单个站点震级的中值。

地方地级(ML)使用公式
18ML=a×log(amp)+b×log(dist)+c×dist+d计算，其中振幅(amp)是伍德—安德森仪记录的地震波形中的最大水平振幅(零峰值；nm)，dist是震源距(km)。

我们使用的校准参数(a，b，c，d)，由Alsaker等19提供，其中a=1.0，b=0.91，c=0.000 87，d=1.67。

为了建立适合当地地质的校准常量，以及校正衰减参数，我们把震级作为距离的函数进行计算(图7)。

观察到，震级随震中距(变化)并没有太大的波动(标准偏差约为±0.108)。

没有任何系统计算震级的距离依赖性表明该校准参数可以适当地描述该地区的衰减。

在表4中，我们比较了由国际地震中心和
印度气象局测定的本研究的4个地震的参数。

我们计算的地震震级接近由国际地
震中心和印度气象局目录确定的值。

作为示例我们在图7中展示了一个地震的结果，图中所示在SEISAN默认参数的帮助下计算的振幅没有显著的偏移。

使用所
述参数，我们计算了16个地震的的震级(表2)。

我们使用程序FPFPIT20计算了6个记录良好的地震的震源机制解，都有明确的P 波初动，台站方位角<125°，这个网格搜索程序把节面和观测到的初动数据的拟合误差达到最小。

对于每个地震，我们使用12～16个P波初动，初动从原始地震记录中手动拾取。

这些地震震级范围为2～3.6。

在图8中，将最大震级事件
(ML=3.6)的P波初动绘制于等积投影的下半球。

震源机制解如图9所示，并在表
5中列出。

断层类型根据Aki和Richards21的约定计算。

因为这些分类是基于倾角，因此有必要利用地质信息来确定这两个节面之间的断层面。

不考虑节面是否是断层面时，所有的震源机制解有相似的断裂特征，震源机制解的节面都显示北西—南东和北东—南西方向表示走滑断层。

这些震源机制解表明，是沿北西—南东
方向的右旋走滑运动，并且这个节面被认为是断层面，地质证据也证实了这一点。

我们得出的结果与Rastogi等7通过复合震源机制解析出的节面一致。

该震源机
制解析显示出右旋运动的平均走滑断层作用，P轴近乎垂直，而T轴几乎是水平的。

伊都基地震起因的科学观点可分为：伊都基构造地震或水库触发地震(RTS)。

我们
下面讨论这两类地震的某些特征，理解地震的本质，并区分这两类地震。

(i)全球没有水库触发例子表明，地震发生在水库开始蓄水20年之后22。

伊都基
附近近期的地震活动发生在1977年活动开始集中爆发后的近35年，因此是否为
水库触发是可疑的。

(ii)2007～2011年这5年里，湖泊水位的周变化和颤动频度这些参数无显著变化3。

此外，地震活动似乎与水库蓄水的速率或者降雨强度无关。

所以我们不能把这个与水库影响和地震活动关联起来。

(iii)湖水水位变化引起的应力变化，可能使断层弱化引起水库触发地震，这一直是一个长期争论的问题。

已有知识表明1.5m/周的水位变化引起在浅震深度(2～
5km)23，240.1bar或更少的应力变化。

这一般适用于大多数水库触发地震，由于最近伊都基地区的地震震源深度达9km，因此这些地震的发生与湖水水位的变化无关。

(iv)就水库触发地震而言，前震活动通常先于主震6，在这种情况下，没有观测到这种情况的报告。

(v)我们计算出的地震震源机制是向走滑兼右旋运动，与Rastogi等71988年报道的伊都基地震一样。

地震震中与伴随KKPT剪切带产生的北西—南东走向的断层密切相关，表明了它们之间存在内在联系。

我们研究了喀拉拉邦伊都基水库西南部附近发生的16个小地震(ML1.5～3.6)的位置和震源机制。

这些地震发生于2011年1月21日至2012年3月，由一个21个台站的临时台阵记录到。

用这些足够数量的地震台站，我们给出了这些事件准确的震源位置和断层面解。

这些地震的震源深度达7.2～9.9km，不太可能是由水库加载所产生的应力触发的，水库加载通常产生的地震在浅的深度(<5km)。

前震活动的缺失和与水位时间序列变化不相关进一步证明了这一点。

近期的地震都局限在北西—南东走向的断层和KKPT剪切带内，在伊都基水库西南方大约6～8km。

利用清晰的P波初动，得到了6个事件的断层面解。

计算的断层面解表明沿着西北—东南走向的走滑兼右旋断层运动；P轴明显沿北—南方向垂直，T轴几乎水平(南方向)。

几个大的轮廓和断层横贯伊都基区域。

贝里亚尔河和KKPT是这个地区的主要轮廓，轮廓主要为剪切带。

该地区的轮廓结构显示了以北西—南东和北东—南西走向的断层为主，震源机制解也支持这个观点。

综合历史和近期的地震记录以及断层面解的结果表明，北西—南东走向的断层和喀拉拉邦中部KKPT剪切带可能是未来发生小到中等地震的潜在区域。

任家琪(1992—)，男，成都理工大学地球物理专业本科，主要从事宁夏地震局监测工作。

E-mail：****************。