A Fast Approximation Algorithm for Housner's Spectral Intensity
-
摘要: Housner谱烈度及修正谱烈度作为基于加速度记录时程直接得到的地震动强度参数,与建筑结构破坏及地震宏观烈度存在较高的相关性,是可靠的地震仪器烈度物理参数指标。然而,相对于地面加速度峰值、地面速度峰值等地震动峰值参数,三分量加速度记录对应的谱烈度计算过程较为复杂,耗时相对较长,影响了利用谱烈度确定地震仪器烈度的时效性。基于对强震动加速度记录的统计分析,本文提出了谱烈度的快速近似算法,仅计算4个方向上的谱烈度值,采用其中3点作圆即可获得水平面内谱烈度迹线的近似最大值,使计算速度提高了45倍,且保持了谱烈度作为地震仪器烈度物理指标的精度。利用在汶川MS 8.0地震等386次MS > 3.0地震中获取的2701组强震动加速度记录,经可靠性检验,结果表明所提出的Housner谱烈度快速近似算法的计算误差在±4.5%以内,可以同时满足地震仪器烈度速报的可靠性和时效性需求。Abstract: As a parameter of ground motion strength derived directly from acceleration records, Housner's spectral intensity and its modified version are highly correlated to building destructure and seismic intensity, and has been adopted as a reliable metric of instrumental seismic intensity. But the algorithm for Housner spectum intensity of a three-component acceleration record was relatively complicated than that of peak ground motion parameter such as peak ground acceleration and peak ground velocity, and this defect of time consuming severely influences the timeliness in instrumental seismic intensity determination. Base on the analyses and ststistics of strong-motion acceleration records, a fast calculating algorithm for spectral intensity is proposed, which obtained the approximate value of the maximum spectral intensity value in the horizontal plane by circle fitting of three points among the spectral intensity values in 4 certain directions. The calculation speed of this algorithm is 45 times higher than that of former algorithms, and the precision accuracy of spectral intensity as a metric of instrumental seismic intensity does not decrease. The reliability testing using 2701 sets of acceleration records in 386 MS > 3.0 earthquakes including Wenchuan MS 8.0 eartuquake indicated that the relative error of this fast calculating algorithm is less than 4.5%, which can satisfied the needs of reliability and timeliness for the rapid report of seismic intensity.
-
引言
滇西北香格里拉—丽江—大理地区位于青藏高原东南缘,红河断裂带的西北端(Allen等,1984;虢顺民等,1996),是我国青藏高原与云贵高原分界区,靠近东喜马拉雅构造节,区内发育一系列大型走滑断裂(Allen等,1991;Tapponnier等,1986;Wang等,1997;1998)。区内现今强震频发,地震地质灾害问题突出(国家地震局地质研究所云南省地震局编,1990),其中1515年永胜M7¾级地震为区内历史记录的最强地震事件(国家地震局震害防御司,1995)。但是,本区地形以山高峡深为主要特征,如虎跳峡,盆地小而少,第四纪地层相对不发育。由于隆升快,因此剥蚀速率大。虽然地震众多,但是,地震遗迹保留少,很少见到有关地震地表破裂(冉勇康等,1990)。近年来我们在进行中甸-大具断裂1:50000活动断裂地质填图时,在大具盆地T3阶地上发现一条长约600m的典型地震地表破裂带。地震地表破裂带是发震断层在地震时留在地表的最直接证据,是了解和研究发震断层运动方式、强度等的最好材料。本文将详细介绍大具盆地新发现地表破裂的基本特征,分析讨论有关形成时代、发震断层归属等,供同行参考。
1. 大具地震地表破裂带位置
云南大具盆地位于云南省丽江市玉龙纳西族自治县北部大具乡与迪庆藏族自治州香格里拉县三坝民族乡江边村之间,是长江上游金沙江中游的一个长宽各约为6km的不规则小型盆地。金沙江从西面玉龙雪山下虎跳峡出口流入大具盆地,然后呈弧形转而从盆地北面流出。盆地上游紧接闻名世界的虎跳峡。大具盆地的主体面为金沙江的第三级阶地面,而第一、二级低阶地只发育在金沙江主河道两侧的局部地点,保留不全。
我们新发现的地表破裂带位于大具盆地的西南侧,近下虎跳峡出口处,金沙江的左岸T3阶地面上,如图 1所示。此阶地面非常平坦。由于大具盆地处于干热河谷环境中,周边山地地层以灰岩为主,因此,各级地貌面顶部都形成一套钙化的砂砾石层,与地貌面时代无关,厚度0.5—1.0m,质地坚硬,此阶地面时代为距今33.8~29.3ka(Kong等,2009)。
图 1 云南大具盆地及其邻区主要活动断裂与强震震中分布图(部分内容据徐锡伟等,2003)Figure 1. Major active faults and epicentral distribution of strong earthquakas in Daju basin and its adjacent areas(Part of the content auovding to Xu Xiwei et al, 2003)F1大具断裂;F2玉龙雪山东麓断裂;F3中甸-剑川断裂;F4丽江-小金河断裂;F5鹤庆-洱源断裂;F6程海断裂;F7德钦-中甸断裂;F8长松平-文明断裂;F9金沙江断裂;F10月亮坪断裂2. 地震地表破裂特征
新发现的地震地表破裂带位于大坝子村西北300m处(见图 2),总体走向340°,长约600m。利用旋翼无人机对地震地表破裂带进行了测绘,获得高精度的正射影像(见图 3a)和DEM(见图 3b),结合野外实际观察认为该破裂带由4条地表破裂组成,它们呈左阶羽列状分布(见图 3c),宽120m左右,每条破裂均较为连续,在挤压岩桥区发育多个北西向的挤压垄脊和断层陡坎,破裂带规模总体由南东向北西逐渐变小,至阶地前缘附近尖灭。在破裂带东侧还发育一级更老的陡坎,高约3m,长约230m,与新生地震地表破裂带亦呈左阶羽列关系。紧邻地表破裂带南东侧见一条废弃的古冲沟,位于T3阶地后缘与冲洪积扇的边界处,其左岸形似弯曲的残脊,被右旋断错约12m(见图 3c),形成断层眉脊,垂直位错约3m,此断错地貌可能是断层多次活动的结果,而非地震地表破裂带的同震位移的表示。
该破裂带地表破裂样式主要由左阶雁列状的张裂缝、挤压鼓包、凹陷坑、挤压垄脊和断层陡坎等构成,基本特征及探槽揭露信息如下:
(1)张裂缝及凹陷坑
张裂缝是该地震地表破裂带中较为发育的破裂样式,主要是张剪或压剪作用的结果,与挤压鼓包相间排列,走向15°—40°,长2—10m,一般呈张开状(见图 4a、4b、4c、4d、4f),宽0.1—0.7m,局部伴有北西盘数厘米至数十厘米正断层状下降现象(见图 4c)。张裂缝在地表破裂带中呈左阶雁列状分布,与破裂带走向之间的夹角在35°—60°。地表破裂呈左阶羽列状,在破裂带端部或走向发生变化的位置也可见到由张裂作用形成的凹陷坑(见图 4e),直径约2m,深1m,呈不规则椭圆形。
(2)挤压鼓包、挤压垄脊及断层陡坎
挤压鼓包是走滑型地震地表破裂带上典型且常见的地震破裂样式之一,是新发现的地震地表破裂带中最为发育的破裂单元,主要是剪切、张剪或压剪作用的结果。由于破裂带发育在大具盆地的一套“钙板层”内,野外调查发现沿破裂带的“钙板层”在受到走滑错动形成的局部挤压应力作用下与下部沉积层脱离,形成形态各异的对称鼓包和不对称鼓包,导致“钙板层”挤压翘起(见图 5a、5b、5g)。根据鼓包形态特征,对称鼓包主要为尖顶帐篷状鼓包(见图 5c),在整个破裂带中隆起最高,高度约1.5m,鼓包翼面倾斜45°,对称轴走向295°左右,与破裂带走向之间的夹角约45°,该形态鼓包仅见一处发育。不对称鼓包可分为对冲鼓包和单斜鼓包,其中,对冲鼓包是指两翼倾角不一、翼面不规则的鼓包(见图 5a),多连续发育,翼面倾斜角度≤40°,总体走向300—310°,与破裂带走向之间的夹角为30—40°,鼓包高度小于1m;单斜鼓包主要表现为“钙板层”挤压翘起,仅有单个翼面(见图 5d、5e),多呈独立个体相间发育,翼面倾斜角度≤30°,总体走向300°—310°,鼓包高度0.2—0.8m。由于地震地表破裂带呈左阶羽列状分布,在阶区内可见挤压垄脊(见图 3c、图 5f)和断层陡坎发育。其中,在阶区南侧见2条挤压垄脊,高约0.4m,宽约4m,走向分别为280°和320°;在阶区北侧见两条断层陡坎,近平行,较连续,高0.3—0.5m,走向325°,属挤压性质陡坎。此外,沿单条地表破裂局部见断层陡坎发育(见图 3c、图 5e),长40—60m,高0.2—0.8m。
该地表破裂带由挤压鼓包、张裂缝、挤压垄脊、凹陷坑等破裂单元组合而成,以各种挤压鼓包组合、张剪或压剪破裂为主,这些组合关系及形态(见图 6)可判断沿破裂主线两侧地块发生过明显的右旋走滑运动,也可以排除这些地表破裂不是可能的滑坡造成的。
(3)探槽揭露
为进一步验证和研究新发现的地震地表破裂带,在地表破裂带的最高端开挖了一个探槽(见图 2)。有关探槽揭示的详细信息将另文介绍。在探槽中,断层有过多次活动,最近一次活动切穿了所有地层,到达地表,与我们地表破裂一致,即破裂面之上没有新的堆积物,运动性质以右旋水平走滑为主,兼有少量的垂直分量(尹功明等,2018)。
3. 分析与讨论
3.1 地表破裂时代分析
由于地表破裂和探槽的最新活动都切穿了所有地层,到达地表;地表破裂又发生在大具盆地约33.8~29.3ka的T3阶地上。因此,其时代很难限定。但是,我们可以作以下分析:1、大具盆地位于玉龙雪山—哈巴雪山东侧、下虎跳峡出口处,为典型的干热河谷区。独特的地形地貌和气候使大具盆地内,常有大风天气,因此,易扬沙尘等。如前面介绍,新发现的地表破裂有张裂缝、凹陷坑、鼓包空隙等,而目前这些空隙中基本未见新的沉积物,局部见杂草和落叶,地表破裂的形态、鼓包内夹杂的原生阶地物质基本没有造受破坏或人为改造(见图 4、图 5)。由于是干热地区,周边玉龙雪山—哈巴雪山的主体岩性为灰岩,在地表容易形成钙板和固结。但是,目前,新发表的地表破裂面上新鲜,没有固结。因此,从破裂带的外形和相关堆积物分析,这个破裂带形成时间应该非常年轻。推测地表破裂可能为距今1000a以来的一次破裂事件。
3.2 地表破裂的归属
地表破裂应归于某一次地震事件的产物,而地震事件可能是近代地震,可能是历史地震,也可能是古地震。自有地震记录以来,大具盆地附近共发生3次6级以上地震,分别为:1933年中甸附近61/4级地震、1966年中甸告湾6.4级地震和1996年丽江7.0级地震(见图 1),其中中甸附近61/4级地震和中甸告湾6.4级地震震中位于大具盆地西北方向,地震造成的破坏向南东方向应该减弱,但大具盆地发现的地表破裂程度是南东强,向北西方向减弱,加之6.5级以下地震一般不会产生地表破裂,所以地表破裂应该不是这两次地震造成。
通过对丽江7.0级地震震源破裂特征(王绍晋等,1997)和震源过程体波反演(徐扬等,1998)的研究,认为该地震主震存在两次破裂过程(见图 7),第一次破裂由起始破裂点(微观震中)沿大具断裂,从东南向北西方向破裂,终止于1966年中甸告湾6.4级地震的震源区,形成该段断裂的贯通;第二次破裂由起始破裂点沿玉龙雪山东麓断裂由北向南破裂,终止于该断裂与丽江—小金河断裂的交接部位,这是丽江地震的主破裂面。从丽江地震的震源破裂过程看,新发现的地震地表破裂带表现出的南东强、北西弱的特征与第一次破裂的过程和方向较一致,似乎可以将地表破裂归属为丽江7.0级地震,但仍存在三个显著的问题:①结合震源机制解、地质构造、等震线图、余震分布等资料,在近南北向压应力场作用下(苏生瑞等,2004),丽江地震发震断层性质为正倾滑为主,兼具左旋剪切错动(王绍晋等,1997;韩竹军等,2004),这与地表破裂表现出右旋走滑为主的特征相矛盾;②丽江地震地表破裂现象(见图 7)主要有地震裂缝呈右阶左旋排列,田坎和公路等发生左旋和垂直位错,最大水平和垂直错距分别为30cm和25cm,雁列裂缝阶区内发育地震鼓包,崩塌与滑坡,地表沉陷和跳石等(张建国等,1997;周光全等,1997;皇甫岗,1997),而新发现的地表破裂由左阶雁列状的张裂缝、挤压鼓包、凹陷坑、挤压垄脊和断层陡坎等构成,最大垂直位错为0.8m,显示右旋走滑特征,两者在破裂类型上有较大差异,且造成的位移规模也有较大差距;③如果地表破裂是丽江地震时形成的,那么丽江地震烈度图(见图 7)(韩新民等,1997)中极震区范围的形态应从大具盆地向北西方向凸起,而不是在大具盆地附近呈较规则的弧形。综上所述,将新发现的地表破裂归属于丽江地震的证据还不充分,因此,该地表破裂可能不是丽江地震而是一次独立地震造成的,其时间可能发生在距今1000a以来。
Figure 7. Simplified map of surface rupture, intensity andrupture direction of Lijiang earthquake(according to Wang Shaojin et al, 1997;Zhou Guangquan et al, 1997;Han Xinmin et al, 1997)3.3 震级讨论
地表破裂和大具断裂的性质是右旋走滑为主,兼有正断性质,局部表现为逆冲特征,震级估计应以水平位错量为准,由于目前发现的地表破裂中很难判定水平同震位移,故无法采用相关公式推算震级。但是,我们可以简单分析:①小于6.5级地震,在地表很难产生地表破裂,我们本次发现的地表破破裂长约600米。因此,地震至少大于6.5级地震。②我们测量获得的由于右旋走滑而产生的隆起鼓包最高达1.5米。新发现的破裂组合特征,完全类似于2001年昆仑山口8.1级左旋走滑大地震形成的地表破裂,只是规模有所差异,拉分和鼓包刚好相反。昆仑山口形成的最大隆起鼓包为2m(徐锡伟等,2002;吴中海等,2004)。因此,也可以推测本破裂的地震应该小于8.1级,很可能为7级以上的地震。③本次新发现的地表破裂带中,也有少量的地震陡坎,显示的最大垂直同震位移为0.8m,向北西方向逐渐减小,根据前人(国家地震局地质研究所等,1990)总结的滇西北地区地震位移量与震级的经验关系M=5.1172+0.95931lgD(M为震级,D为一次地震的位错量),可近似求得产生该地表破裂的地震震级为7.0级左右。当然,新发现的地表破裂中主要的运动特征是右旋走滑,因此,根据垂直位移估算的震级,应该小于实际震级。根据上述分析,产生本次新发现地表破裂的地震至少达7级。
3.4 发震断层分析
如图 2所示,本次新发现的地表破裂带位于下虎跳峡出口处。而这个部位恰好为玉龙雪山东麓断层和中甸-大具断裂的交汇处。
玉龙雪山东麓断层北起自大具盆地,向南经玉湖、白沙西,最终与丽江-小金河断裂交汇,全长约52km,总体走向近南北。对该断层的研究目前有一些不同认识:该断层最初是在丽江幅1:20万地质图及区域地质调报告中提出的,认为该断层北起大具盆地,向南沿大沟村,经雪花村西、白水林场、里那比格,终止于丽江玉湖北,根据基岩面上的陡倾擦痕认定为正断层;玉龙雪山东麓断层不仅是川滇菱形块体西南边界断层,也是我国三江挤压区的东部边界,其运动性质是以右旋走滑为主,兼具逆冲分量(尹功明等,2017); 1996年丽江地震发生后,王运生等(2000)针对此断层开展了地球化学探测等工作,认为次断层是向山体内部倾的逆冲断层;也有学者认为玉龙雪山东麓断层为正断层(俞维贤等,2004;张西娟等,2006;石旭华等,2008;吴中海等,2008;Kong等,2009)。
中甸-大具断裂北起自中甸南,向南东经小中甸东、马家村、哈巴、大具、大东,后止于金沙江边,总体走向310—320°,倾向NE或SW,倾角较陡,一般在60°以上,运动性质以右旋走滑为主,兼有正断分量活动,局部显示逆冲性质。该断裂目前是川滇菱形块体的西侧边界断层(沈军等,2001),也是一条大陆活动块体的边界断裂(徐锡伟等,2003)。现今GPS数据表明,川滇菱形块体相对于欧亚大陆有高达10mm/a向南东方向的运移速率(Gan等,2007),在向南东、南东方向平移的同时,川滇块体内部的次级块体还发生了顺时针旋转(徐锡伟等,2003)。因此,中甸-大具断裂作为边界断裂具有很强的活动性。
中甸-大具断裂是川滇菱形块体西边界的重要组成部分,断裂附近发生过多次中强地震,如1933年中甸附近61/4级地震、1961年中甸6.0级地震、1966年中甸告湾6.4级地震。沈军等(2001)根据大具盆地三级阶地面上发育的小冲沟右旋位错,平均右旋错动速率为4~6mm/a;根据中甸盆地北缘、小中甸以东古河道位错量获得平均右旋错动速率为4~6mm/a;尹功明等(2018)根据大具盆地T3级阶地位错,初步获得的右旋错动速率为5.03~8.37mm/a,垂直活动速率为2.37~3.94mm/a。表明中甸-大具断裂的右旋平均位错速率非常高,是一条强烈活动的次级块体的边界断层。
从本次新发现的地震地表破裂所在位置发育的两条断层的运动性质和强度分析,新地表破裂应该归属于中甸-大具断裂。
4. 结论
(1)新发现的地震地表破裂带位于云南大具盆地下虎跳峡出口处、金沙江右岸T3级阶地上,破裂样式主要包括左阶雁列状的张裂缝、挤压鼓包、凹陷坑、挤压垄脊和断层陡坎等,破裂较为连续,长约600m,宽120m左右,表现为右旋走滑特征,兼有倾滑分量。
(2)根据地表破裂长度、隆起鼓包高度和以右旋走滑为主兼有的倾滑分量,震级应该在7.0级左右甚至更大。新发现的地震地表破裂带可能不是丽江7.0地震产生的,而是一次独立的地震事件,其发生时代不会久远,可能不会超过距今1000a。
(3)云南大具新发现的地震地表破裂带的发震断层应该为中甸-大具断裂。
致谢: 中国地震局地质研究所徐锡伟研究员、汪一鹏研究员、冉勇康等亲临野外现场指导,中国地震灾害防御中心李峰研究员、高战武研究员在撰写中提出了许多宝贵意见,在此一并表示衷心感谢。 -
表 1 算法可靠性验证所用强震动记录的数量统计
Table 1. Acceleration records used in reliability veritification of fast calculating algorithm for Housner spectum intensity
震级 PGA/cm·s-2 PGA<3 3≤PGA<10 10≤PGA<30 30≤PGA<100 100≤PGA<300 300≤PGA<500 PGA≥500 MS=8.0 64 135 132 64 44 6 5 MS=7.0 20 39 31 21 7 5 3 MS=6.5 26 31 11 3 2 1 6.0≤MS<6.5 21 84 76 29 7 5.5≤MS<6.0 29 103 65 40 8 1 5.0≤MS<5.5 36 201 89 31 6 4.5≤MS<5.0 74 270 119 27 3 4.0≤MS<4.5 96 296 145 30 7 3.0≤MS<4.0 18 84 43 13 -
崔建文, 刘琼仙, 段建新等, 2014.2014年云南鲁甸6.5级地震强震动观测记录及初步分析.地震研究, 37(4):542-548. doi: 10.3969/j.issn.1000-0666.2014.04.009 金星, 张红才, 李军等, 2013.地震仪器烈度标准初步研究.地球物理学进展, 28(5):2336-2351. http://d.old.wanfangdata.com.cn/Periodical/dqwlxjz201305016 赖敏, 朱建钢, 朱永莉, 2014.芦山地震数字强震动观测记录及其初步分析.四川地震, (2):6-11. http://d.old.wanfangdata.com.cn/Periodical/scdz201402002 李敏, 2010.地震动加速度反应谱与地震烈度的关系研究.哈尔滨:中国地震局工程力学研究所, 44-52. 李世成, 崔建文, 张潜等, 2011.2009年姚安6.0级地震强震动应急观测的数据处理及其震害启示.灾害学, 26(4):82-88. doi: 10.3969/j.issn.1000-811X.2011.04.016 李世成, 林国良, 崔建文等, 2013.宁蒗-盐源5.7级地震数字强震动台网记录及初步分析.地震研究, 36(4):502-507. doi: 10.3969/j.issn.1000-0666.2013.04.015 李山有, 金星, 陈先等, 2002.地震动强度与地震烈度速报研究.地震工程与工程振动, 22(6):1-7. doi: 10.3969/j.issn.1000-1301.2002.06.001 林淋, 孙景江, 2011.不同地震动参数与地震烈度相关性对比研究.地震工程与工程振动, 31(1):6-10. http://d.old.wanfangdata.com.cn/Periodical/dzgcygczd201101002 马强, 李水龙, 李山有等, 2014.不同地震动参数与地震烈度的相关性分析.地震工程与工程振动, 34(4):83-92. http://d.old.wanfangdata.com.cn/Conference/8421967 王玉石, 周正华, 王伟, 2008.基于假设检验的地震动强度(烈度)速报方法.地震工程与工程振动, 28(5):49-54. http://d.old.wanfangdata.com.cn/Periodical/dzgcygczd200805007 王玉石, 周正华, 兰日清, 2010.利用修正谱烈度确定我国西部地区仪器烈度的建议方法.应用基础与工程科学学报, 18(S1):119-129. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201001943620 王玉石, 李小军, 梅泽洪等, 2013.几种仪器烈度算法在汶川地震与芦山地震中的可靠性比较.地震学报, 35(5):759-770. doi: 10.3969/j.issn.0253-3782.2013.05.014 朱永莉, 黎大虎, 龙承厚等, 2015.不同仪器烈度算法在四川地区历次地震中的比较应用.地震学报, 37(2):355-346. doi: 10.3969/j.issn.1000-0844.2015.02.0355 Atkinson G. M., Kaka S. I., 2007. Relationships between felt intensity and instrumental ground motion in the Central United States and California. Bulletin of the Seismological Society of America, 97 (2):497-510. doi: 10.1785/0120060154 Boatwright J., Thywissen K., Seekins L. C., 2001. Correlation of ground motion and intensity for the 17 January 1994 Northridge, California, earthquake. Bulletin of the Seismological Society of America, 91 (4):739-752. doi: 10.1785/0119990049 Bradley B. A., 2009. Site-specific and spatially distributed ground-motion prediction of acceleration spectrum intensity. Bulletin of the Seismological Society of America, 100 (2):792-801. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=4868394826820fe0e80a13b869a474f7 Bradley B. A., Cubrinovski M., Dhakal R. P., et al., 2009a. Intensity measures for the seismic response of pile foundations. Soil Dynamics and Earthquake Engineering, 29 (6):1046-1058. doi: 10.1016/j.soildyn.2008.12.002 Bradley B. A., Cubrinovski M., MacRae G. A., et al., 2009b. Ground-motion prediction equation for SI based on spectral acceleration equations. Bulletin of the Seismological Society of America, 99 (1):277-285. doi: 10.1785/0120080044 Bradley B. A., Dhakal R. P., MacRae G. A., et al., 2010a. Prediction of spatially distributed seismic demands in specific structures:ground motion and structural response. Earthquake Engineering and Structural Dynamics, 39(5):501-520. doi: 10.1002/eqe.954/full Bradley B. A., Dhakal R. P., MacRae G. A., et al., 2010b. Prediction of spatially distributed seismic demands in specific structures:structural response to loss estimation. Earthquake Engineering and Structural Dynamics, 39 (6):591-613. doi: 10.1002/eqe.955/full Bradley B. A., 2011. Empirical correlation of PGA, spectral accelerations and spectrum intensities from active shallow crustal earthquakes. Earthquake Engineering and Structural Dynamics, 40 (15):1707-1721. doi: 10.1002/eqe.v40.15 Chernov Y. K., Sokolov V. Y., 1999. Correlation of seismic intensity with Fourier acceleration spectra. Physics and Chemistry of the Earth, Part A:Solid Earth and Geodesy, 24 (6):523-528. doi: 10.1016/S1464-1895(99)00065-4 Fahjan Y. M., Alcik H., Sari A., 2011. Applications of cumulative absolute velocity to urban earthquake early warning systems. Journal of Seismology, 15 (2):355-373. doi: 10.1007/s10950-011-9229-8 Housner G. W., 1952. Spectrum intensities of strong-motion earthquakes. In: Proceedings of the Symposium on Earthquake and Blast Effects on Structures. Los Angeles, California: Earthquake Engineering Research Institute, 20-36. Iwata T., 1991. Earthquake sensing device by SI (spectrum intensity) value for system shut-down use. Journal of High Pressure Institute of Japan, 29 (3):137-144. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=J-STAGE_531922 Karim K. R., Yamazaki F., 2001. Estimation of the JMA seismic intensity from liquefied records and its correlation with strong motion parameters. In: Proceedings of the 26th Japan Earthquake Engineering Symposium. Sapporo, Japan: J-STAGE, 253-256. Karim K. R., Yamazaki F., 2002. Correlation of JMA instrumental seismic intensity with strong motion parameters. Earthquake Engineering and Structural Dynamics, 31 (5):1191-1212. doi: 10.1002/(ISSN)1096-9845 Katayama T., Sato N., Ohbo N., et al., 1986. Ground shaking severity detector by use of spectrum intensity (SI). In: Proceedings of the 7th Japan Earthquake Engineering Symposium. Tokyo, Japan: J-STAGE, 7: 373-378. Koganemaru K., Shimizu Y., Nakayama W., et al., 2000. Development of a new SI sensor. In: Proceedings of the 12th World Conference of Earthquake Engeering. Auckland, New Zealand: WCEE. Linkimer L., 2008. Relationship between peak ground acceleration and modified Mercalli intensity in Costa Rica. Revista Geológica de América Central, 38:81-94. Margottini C., Molin D., Serva L., 1992. Intensity versus ground motion:a new approach using Italian data. Engineering Geology, 33 (1):45-58. doi: 10.1016/0013-7952(92)90034-V Martínez-Rueda J. E., 1998. Scaling procedure for natural accelerograms based on a system of spectrum intensity scales. Earthquake Spectra, 14 (1):135-152. doi: 10.1193/1.1585992 Nau J. M., Hall W. J., 1984. Scaling methods for earthquake response spectra. Journal of Structural Engineering, 110 (7):1533-1548. doi: 10.1061/(ASCE)0733-9445(1984)110:7(1533) Shabestari K. T., Yamazaki F., 2001. A proposal of instrumental seismic intensity scale compatible with MMI evaluated from three-component acceleration records. Earthquake Spectra, 17 (4):711-723. doi: 10.1193/1.1425814 Shimizu Y., Koganemaru K., Yamazaki F., et al., 2000. Seismic motion observed in Taipei basin by new SI sensors and its implication to seismic zoning. In: Proceedings of the 6th International Conference on Seismic Zonation. Palm Springs, California: Earthquake Engineering Research Institute, 497-502. Sokolov V. Y., Chernov Y. K., 1998. On the correlation of seismic intensity with Fourier amplitude spectra. Earthquake Spectra, 14 (4):679-694. doi: 10.1193/1.1586022 Sokolov V. Y., 2002. Seismic intensity and fourier acceleration spectra:revised relationship. Earthquake Spectra, 18 (1):161-187. doi: 10.1193/1.1469037 Takubo K., Yanada T., Furukawa H., et al., 1999. Innovative spectral intensity transducer using three-axial micromachining accelerometer for earthquake crisis management. In: Proceedings of the 7th IEEE International Conference on Emerging Technologies and Factory Automation. Barcelona, Spain: IEEE, 1: 379-384. https://www.researchgate.net/publication/3831381_Innovative_spectral_intensity_transducer_using_three-axial_micromachining_accelerometer_for_earthquake_crisis_management Trifunac M. D., Brady A. G., 1975. A Study on the duration of strong earthquake ground motion. Bulletin of the Seismological Society of America, 65 (3):581-626. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=954aacfd06f585487542ce832a2bf632 Trifunac M. D., Westermo B., 1977. A note on the correlation of frequency-dependent duration of strong earthquake ground motion with the modified Mercalli intensity and the geologic conditions at the recording stations. Bulletin of the Seismological Society of America, 67 (3):917-927. Ueong Y. S., 2009. A study on feasibility of SI for identification of earthquake damages in Taiwan. Soil Dynamics and Earthquake Engineering, 29 (1):185-193. doi: 10.1016/j.soildyn.2007.12.005 Wald D. J., Quitoriano V., Heaton T. H., et al., 1999. Relationships between peak ground acceleration, peak ground velocity, and modified Mercalli intensity in California. Earthquake Spectra, 15 (3):557-564. doi: 10.1193/1.1586058 Wang Q., Fu J. H., Wang Z. Y., et al., 2012. A seismic intensity estimation method based on the fuzzy-norm theory. Soil Dynamics and Earthquake Engineering, 40:109-117. doi: 10.1016/j.soildyn.2012.03.010 Wang Y. S., Li X. J., Zhou Z. H., 2013. A new instrumental measure of epicentral shaking intensity in Western China. Bulletin of Earthquake Engineering, 11 (4):913-924. doi: 10.1007/s10518-013-9428-5 Wu Y. M., Teng T. L., Shin T. C., et al., 2003. Relationship between peak ground acceleration, peak ground velocity, and intensity in Taiwan. Bulletin of the Seismological Society of America, 93 (1):386-396. doi: 10.1785/0120020097 -