刘小凤,梅秀苹,冯建刚,程建武,张 辉
(甘肃省地震局,甘肃 兰州 730000)
区域地震活动状态判定在地震趋势预测中具有重要意义。而以往研究大多基于定性的经验总结,地震活动性异常的确认与描述基本处于定性的状态;其次是不同区域小震定位精度存在较大差异;第三是不同地区历史地震目录完整性差异较大。这些导致我们对区域地震活动状态判定具有一定的局限性。
本文根据区域构造背景、应力场状态及动力学环境等,将青藏高原北部地区划分出4个与构造及地震活动相关的统计单元,在对各构造单元地震监控能力分时段评估的基础上,从不同震级范围的地震频次、应变释放、强度变化、G-R关系b值、地震时间间隔等基本参数入手,定量判定青藏高原北部地区不同构造区带平均状态下地震活动的强弱,在更精细的层次上对不同构造单元地震活跃和平静状态进行定量、半定量评估。
青藏高原北部地区地质构造复杂,地震活动频度高、强度大。根据区域活动构造定量数据、GPS和跨断层形变观测、震源机制解等研究结果[1-4],将青藏高原北部地区划分为祁连山地震带、甘东南地区、柴达木-共和地区和库玛地震带4个相对独立的统计单元(图1)。
图1 青藏高原北部地区构造单元划分图Fig.1 Tectonic zoning of northern region of Qinghai-Xizang plateau.
祁连山地震带四周被深大走滑活动断裂包围,在NNE向区域应力场作用下地块内部发生挤压逆冲,形成一系列推覆体构造,区内发生过1920年海原8.5级地震、1927年古浪8级地震。震源机制解显示该地区中强以上地震断层错动类型以左旋走滑为主兼有一定的逆冲分量,区域主压应力方向自西向东由NNE、NE向到NEE和近EW向。
甘东南地区是青藏高原东北缘一系列NWNWW向构造和青藏块体东部近SN向构造的过渡地区,地壳活动十分活跃,区内构造复杂,地震活动强烈,曾发生1654年天水8级地震、1879年武都8级地震。震源机制解显示,该地区中强以上地震断层错动类型以左旋走滑为主兼有一定的逆冲分量,区域主压应力P轴优势方向主要为NW和近SN向,整体构造应力场在较长时间内基本上保持稳定状态。自晚更新世以来一直处于NE-NEE向区域构造应力场作用下,导致块体内部向东南方向的挤出。
柴达木—共和地块内部构造和介质性质比较简单,周边受区域深大活动断裂围限,内部发育NNW和NW向次级断裂。在NE向区域应力场作用下NW-NWW向边界断裂发生左旋走滑活动,NNW向次级断裂发生右旋走滑活动,导致共和-兰州地区发生顺时针旋转和整个块体向南东方向的挤出。块体内发生的最大地震为1990年共和7.0级。震源机制解显示柴达木-共和地块中强以上地震断层错动类型以左旋走滑为主兼有一定的逆冲分量,区域主压应力方向主要为NE向,晚更新世以来该地区的运动状态为向右的旋转运动。
库玛活动断裂带是青藏块体内部长度超过千里的一条重要的左旋走滑活动断裂带,1900年以来库玛地震带及附近地区共发生7级以上地震11次,最大地震为2001年昆仑山口西8.1级地震。震源机制解显示库玛地震带中强以上地震断层错动类型以左旋走滑为主,区域主压应力方向主要为NENEE向。
地震目录的完整性是指某一震级以上的地震都能被观测或记录到,如果这一震级以上地震出现遗漏,即称为不完整地震目录[5]。已有的研究结果表明,分析历史地震目录完整性最常见的方法是利用古登堡-里克特的震级-频度关系[6],计算出最小完整性震级,在时间和震级尺度基本不变的条件下,利用震级-频度关系对地震目录进行不完整性定量估计,确定目录的遗失率等。本文采用古登堡-里克特的震级-频度关系(G-R关系)对青藏高原北部地区历史地震目录的完整性进行了定量分析,通过地震频度和b值拟合,大致确定了不同下限震级的完整性目录起始时间。
2.1.1 频度曲线
本研究中,分别采用预防性护理和常规护理对两组继发性癫痫患者进行干预,并比较两组护理效果,结果显示,住院期间不良事件发生率方面,研究组(2.2%)显著低于对照组(15.6%),护理总有效率方面,研究组(97.8%)显著高于对照组(86.7%),说明对继发性癫痫患者采用预防性护理干预,不仅能降低其住院期间不良事件发生率,还能提高其整体预后效果。
利用各个构造单元地震目录,删除5级以上强余震,分别以M5.0和M6.0做为震级下限,20年窗长,2年滑动步长做频度曲线,得到地震频度随时间的变化特征。根据西北地区各构造单元不同震级下限频度曲线,祁连山地震带及甘东南地区从1570年开始历史地震记载较多,1920年之后地震目录相对完整;柴达木-共和地块和库玛地震带从1900年以后才有地震记载,柴达木-共和地块1935年之后地震目录相对完整。表明频度曲线存在一个相对明显的变化过程,显示出地震目录由不完整趋向完整,处在一个稳定的变化范围(图2)。
图2 青藏高原北部各构造单元不同震级下限频度图Fig.2 Different completeness magnitude of tectonic units in the northern region of Qinghai-Xizang plateau.
2.1.2 b值拟合
b值拟合主要反映出了低震级端地震目录的完整性及起始时间。利用各个统计单元5级以上地震目录,删除5级以上强余震,以5.0级地震做为震级下限,得到b值拟合图。从中可以看出,1900年以来青藏高原北部地区5级以上地震记录基本完整,四个区域b值差异不大,库玛地震带b值略低,这反映了该区域的应力水平相对较高,库玛地震带处于青藏块体内部,构造活动较强烈(图3)。
根据震级-频度关系[6-7]对青藏高原北部地区历史地震目录完整性进行了分析,通过地震频度、b值拟合的估计,初步确定了不同下限震级的完整目录起始时间:祁连山地震带6级以上地震自1569年以后是相对完整的,5级以上地震目录自1920年以后是完整的;甘东南地区6级以上地震自1575年以后是相对完整的,5级以上地震目录自1920年以后是完整的;柴达木-共和地块6级以上地震自1935年以后是相对完整的,5级以上地震目录自1950年以后是完整的;库玛地震带6级以上地震自1935年以后是相对完整的,5级以上地震目录自1950年以后是完整的。
仪器记录地震目录是地震活动分析研究中的关键资料,随着地震监测台网的不断改造升级,地震监测能力有了很大的提高,完整地震目录的震级下限也在不断提高。对4个构造单元基于古登堡-里克特的震级-频度关系对每个区域不同时间段小震目录进行完整性分析。采用1970年以来全国小震目录,利用K-K法进行余震删除,再删除震群,最后用最小二乘法拟合不同区域、不同时期的b值,根据b值的拟合情况分析完整目录的震级下限。
图3 青藏高原北部1900年以来5级以上地震b值拟合图Fig.3 Fitting of b value of M≥5.0earthquakes from 1900in the northern region of Qinghai-Xizang plateau.
青藏高原北部地区不同构造单元分时段b值拟合结果显示,完整地震目录的震级下限随着台网的改造升级而不断降低,4个统计单元中完整目录震级下限有一定的差异,祁连山地震带和甘东南地区震级下限要低于柴达木-共和和库玛地震带,各区域完整地震目录的震级下限分别为:祁连山地震带ML2.6、甘东南地区 ML2.3、库玛地震带 ML3.1、柴达木-共和地块 ML3.2[8]。
1900年以来5级以上地震的研究结果显示,青藏高原北部地区大体可以划分为6个集中丛集活跃时段,每个时段持续时间不等,地震频度和强度也差异较大[9]。
(1)祁连山地震带
1900年以来祁连山地震带5级以上地震表现出4个集中丛集活动的阶段(图4(a)),第一阶段为1920-1934年,连续发生了1920年海原8.5级、1927年古浪8.0级和1932年昌马7.6级地震,间隔时间分别为7年和5年,并依次由东向西迁移;第二阶段为1951-1962年,最大地震为1954年山丹7.3级,发生在祁连山地震带北部的龙首山断裂带上;第三阶段为1980-1993年,以中强地震活动为主,包含了3次6级地震,最大地震为1986年门源6.4级;第四阶段始于2000年,发生了5次地震,具有自东向西再折返向东迁移的特征,到目前为止最大地震为2003年民乐6.1级。每个活跃阶段持续10年左右,一般至少发生5.5级以上地震5次,6.0级以上地震3次,每个活跃阶段之间为15年左右的平静期。
图4 祁连山地震带和甘东南地区5级以上地震M-t图Fig.4 M-t diagrams of M≥5.0earthquakes in Qilian mountain seismic belt and southeastern Gansu province.
(2)甘东南地区
甘东南-甘青川交界地区5级以上地震具有平静和活跃交替活动特征,大体可划分为6个活跃阶段(图4(b))。第一阶段为1573-1654年,发生4次6级地震和1次8级地震,最大地震为1654年天水8.0级;第二阶段为1718-1765年,发生2次6级地震和1次7级地震,最大地震为1718年通渭7.5级;第三阶段为1819-1850年,发生1次6级地震,最大地震为1837年岷县6.0级;第四阶段为1879-1936年,发生4次6级地震和1次8级地震,最大地震为1879年武都8.0级;第五阶段为1957-2006年,发生3次6级地震和2次7级地震,最大地震为1976年松潘7.2级。每个活跃阶段持续时间一般为31-81年,平均约50年,每个活跃阶段至少发生4次5.5级以上地震或1次6级以上地震。活跃阶段之间为21-64年的平静时间,平静时段内一般无5级以上地震发生。
(3)柴达木-共和地块
图5 1900年以来柴达木-共和地块和库玛地震带5级以上地震M-t图Fig.5 M-t diagrams of M≥5.0earthquakes from 1900 of Qaidam-Gonghe mass and Kuma seismic belt.
1900年以来柴达木-共和地块5级以上地震大体可以划分出5个阶段,每个阶段持续时间不等(图5(a))。第一阶段为1930-1938年,持续时间8年,最大地震为1938年天峻6.0级;第二阶段为1949-1964年,持续时间15年,最大地震为1962年霍布逊湖6.8级;第三阶段为1972-1980年,持续时间8年,最大地震为1977年茫崖6.4级;第四阶段为1985-1995年,持续时间10年,最大地震为1990年共和7.0级;第五阶段为2000-2008年,持续时间8年,最大地震为2000年兴海6.6级、2003年德令哈6.6级地震。由此可见,该地区地震活动强度不大,但频次较高,最大地震为1990年共和7.0级。活跃期一般持续10年左右,至少发生2次6级地震。平静期持续时间不等,分别为11年、12年、15年和5年,平静期无5级以上地震发生。
(4)库玛地震带
1900年以前库玛地震带历史地震资料遗漏,未记录到5级以上地震,1900年以后特别是1950年以前地震记录也不完整。根据现有的地震资料,大体可以划分出4个活跃阶段,第一阶段为1931-1940年,持续时间为9年,最大地震为1937年托索湖7.5级;第二阶段为1947-1952年,持续时间为5年,最大地震为1947年达日7.7级;第三阶段为1961-1971年,持续时间为10年,最大地震为1963年阿兰湖7.0级;第四阶段为1995-2006年,持续时间为12年,最大地震为2001年昆仑山口西8.1级(图5(b))。在第四活跃阶段中若不考虑8.1级地震的余震活动,则持续时间为7年。活跃期一般持续10年左右,发生5级以上地震5-9次,至少发生1次7级以上地震。平静期持续时间不等,分别为7年、9年和23年。
根据有记载以来5级以上地震目录,结合前述地震目录完整性分析结果,确定了各个构造单元统计地震的起始时间。采用 Mapseis2.9.3软件中的计算方法,分别获得了4个分区地震活跃与平静期的地震活动状态参数,其中m值计算采用蒋海昆提供的DA模式。DA模式的含义是,大震前临界区范围内地震活动呈加速态势,同一时期在孕震区较小范围内则出现地震活动的减弱。震前孕震区范围地震活动的减弱被归因于由于震前滑移所导致的应力松弛过程[10-11]。
不同构造单元5级以上地震活跃与平静期地震活动状态参数统计结果显示,祁连山地震带和甘东南地区地震活跃期持续时间明显大于平静期,柴达木-共和地块和库玛地震带地震活跃期与平静期持续时间差异不大。活跃期一般发生5级地震5~27次,6级地震1~6次,7级以上地震1~2次,5级以上地震年频次为0.5~1.8,6级以上地震年频次为0.03~0.7,折合震级为6.1~8.5,平均震级为5.5~7.7,b值为0.4~0.8,m 值为1.2~1.5。平静期发生5级以上地震1~3次或不发生5级以上地震,均无6级以上地震发生(表1)。
根据1970年以来ML2.0以上地震目录,采用Mapseis2.9.3软件中的计算方法,分别获得了4个区的地震活动状态参数。需要说明的是,分别使用不删除余震和用K-K法删除余震两套地震目录分别计算各单元地震地震活跃、平静时段小震活动参数特征,其结果显示删除余震后活跃期与平静期各参数的差异性减小,而不删除余震的差异性较大,因此以下定量指标判据均为未删除余震计算的结果。
不同构造单元小震活跃与平静时段地震活动状态参数统计结果显示(表2),祁连山地震带和柴达木-共和地块活跃时段平均持续时间较长,而甘东南与库玛地震带平静时段平均持续时间较长。活跃时段ML≥3.0、ML≥4.0和ML≥5.0地震平均年频次和月频次明显大于平静时段,活跃时段平均年应变释放折合震级一般大于5级,平均月应变释放折合震级一般大于4级,b值小于0.8,m值也较低。平静时段平均年应变释放折合震级小于5级,平均月应变释放折合震级小于4级,b值大于0.8,m值一般较高。
表1 不同构造单元5级以上地震活跃与平静期地震活动状态参数统计表
表2 1970-2008年不同构造单元小震活跃与平静时段地震活动状态参数统计表
(1)根据区域构造地质背景、形变观测及震源机制解揭示出的应力场状态及区域动力学环境等,结合微震精定位后呈现出的时空强和震源深度分布特征,将青藏高原北部地区划分为四个独立的统计构造单元,即祁连山构造带、甘东南地区、柴达木-共和地块和库玛构造带。
(2)通过地震频度和b值拟合,确定了各个地区不同下限震级的完整性目录起始时间,其中祁连山地震带6级以上地震自1569年以后是相对完整的,5级以上地震目录自1920年以后是完整的;甘东南地区6级以上地震自1575年以后是相对完整的,5级以上地震目录自1920年以后是完整的;柴达木-共和地块6级以上地震自1935年以后是相对完整的,5级以上地震目录自1950年以后是完整的;库玛地震带6级以上地震自1935年以后是相对完整的,5级以上地震目录自1950年以后是完整的。仪器记录地震精度分析结果显示,1970年以来区域完整地震目录的震级下限分别为,祁连山地震带ML2.6,甘东南地区 ML2.3,库玛地震带 ML3.1,柴达木-共和地块ML3.2。
(3)各构造单元中强以上地震存在活跃与平静交替活动的特征。祁连山地震带5级以上地震表现出4个集中丛集活动的阶段,每个活动阶段持续10年左右,每个活跃阶段内一般至少发生5.5级以上地震5次,6.0级以上地震3次,每个活跃阶段之间为15年左右的平静期;甘东南地区可划分为6个活跃阶段,每个活跃阶段持续时间一般为31-81年,平均约50年,每个活跃阶段至少发生4次5.5级以上地震或1次6级以上地震;活跃阶段之间为21~64年的平静时间,平静时段内一般无5级以上地震发生;柴达木-共和地块可以划分出5个阶段,每个阶段持续时间不等,活跃期一般持续10年左右,至少发生2次6级地震。平静期持续时间不等,分别为11年、12年、15年和5年,平静期无5级以上地震发生;库玛地震带可以划分出4个活跃阶段,活跃期一般持续10年左右,发生5级以上地震5~9次,至少发生1次7级以上地震。平静期持续时间不等,分别为7年、9年和23年。
(4)定量计算了祁连山地震带、甘东南地区、柴达木-共和地块及库玛地震带中强以上地震活动特征参数及平静期地震活动特征参数。不同构造单元5级以上地震活跃与平静期地震活动状态参数统计结果显示,地震活跃期一般发生5级地震5~27次,6级地震1~6次,7级以上地震1~2次;5级以上地震年频次为0.5~1.8,6级以上地震年频次为0.03~0.7;折合震级为6.1~8.5,平均震级为5.5~7.7;b值为0.4~0.8,m 值为1.2~1.5。平静期发生5级以上地震1~3次或不发生5级以上地震,均无6级以上地震发生。获得了各个区域地震活跃与平静状态下的小震活动状态参数,祁连山地震带和柴达木-共和地块活跃时段平均持续时间较长,而甘东南与库玛地震带平静时段平均持续时间较长;活跃时段ML≥3.0、ML≥4.0和 ML≥5.0地震平均年频次和月频次明显大于平静时段;活跃时段平均年应变释放折合震级一般大于5级,平均月应变释放折合震级一般大于4级;b值小于0.8,m值也较低。平静时段平均年应变释放折合震级小于5级,平均月应变释放折合震级小于4级,b值大于0.8,m值一般较高。
[1]袁道阳,张培震,刘百篪,等.青藏高原东北缘晚第四纪活动构造的几何图像与构造转换[J].地质学报,2004,78(2):270-278.
[2]刘百篪,曹娟娟,袁道阳,等.青藏高原北部活动地块内部的活断层定量资料[J].地震地质,2008,30(1):161-175.
[3]张希,张晓亮,王双绪,等.青藏块体东北缘近期地壳水平运动与应变积累[J].大地测量与地球动力学,2008,28(4):12-16.
[4]刘小凤,刘百篪,杨立明.祁连山中东段三维构造物理模型及其在地震预报中的应用[J].西北地震学报,2000,22(2):-.
[5]王海涛,李莹甄,屠弘伟.新疆历史地震目录完整性分析[J].内陆地震,2006,20(1):10-17.
[6]Gutenberg B,Richter C F.Earthquake magnitude,intensity,energy and acceleration[J].Bull.Seism.Soc.Am.,1956,46:105-143.
[7]Richter C F.Elementary seismology[M].San Francisco:Freeman,1958.
[8]冯建刚,姚家骏,代炜.青藏高原北部地区小震目录完整性分析[J].高原地震,2010,22(2):10-14.
[9]杨立明,肖丽珠,张小美,等.青藏高原北部地震活动形势及其与中国大陆地区地震活动的关系研究[J].地震,2001,23(3):243-249.
[10]蒋海昆,吴琼,董祥,等.不同温压条件下声发射应变能释放特征——加速模型参数物理含义的初步讨论[J].地球物理学报,2009,52(8):2064-2073.
[11]蒋海昆,董祥,苗青壮,等.中国大陆7级强震前中小地震应变释放特征[J].地震,2009,29(3):1-10.