地球物理学报  2018, Vol. 61 Issue (1): 183-198   PDF    
红河断裂带库仑应力演化及未来地震危险性估计
尹凤玲1, 蒋长胜1, 韩立波1, 张怀2, 张贝1,2     
1. 中国地震局地球物理研究所, 北京 100081;
2. 中国科学院计算地球动力学重点实验室, 中国科学院大学, 北京 100049
摘要:现今地震活动性显示红河断裂带强震主要发生在断裂带的北段,中段长期存在地震空区.为了更好地理解红河断裂带不同段落地震活动性差异和评估未来潜在地震危险性,本研究基于分层半无限空间黏弹性地球模型计算了自1833年嵩明M8.0地震以来,红河断裂带上及其周边共25次强震由于同震应力阶变、震后黏滞松弛和震间构造应力加载的综合作用,红河断裂带上库仑应力变化的演化过程.结果表明,在近180年红河断裂带南北两段得到加载,中段始终位于应力影区,这可能部分解释红河断裂带南北段地震相对活跃、中段长期存在的地震平静现象.假设未来三十余年该地区不再发生大震,考虑震后和震间作用,红河断裂带中段应力影区仍然存在,但范围在缩小;洱源附近、大理至大斗门以北地区、元江以南地区应力增加超过0.1 MPa,可能仍是地震潜在危险区段.
关键词: 红河断裂带      震后黏弹松弛      震间构造加载      库仑应力变化      地震危险性     
Seismic hazard assessment for the Red River fault: insight from Coulomb stress evolution
YIN FengLing1, JIANG ChangSheng1, HAN LiBo1, ZHANG Huai2, ZHANG Bei1,2     
1. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
2. Key Laboratory of Computational Geodynamics of Chinese Academy of Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: Historic strong earthquakes ruptured the northern segment of the Red River fault, while left a seismic gap in the middle segment. To better understand this seismicity pattern and assess the seismic hazard in this region, we investigated the evolution of Coulomb failure stress on the Red River fault since 1833 taking coseismic, postseismic and interseismic deformation into account based on a stratified viscoelastic earth model. In total, 25 large earthquakes were analyzed. The results revealed that both the northern and the southern segments of the Red River fault are loaded, while the stress changes of the central segment remained negative in recent 180 years. This may contribute to the high level of seismicity near the northern and southern segments of the Red River fault and the low level of seismicity in the central segment. The combined prediction of the postseismic and interseismic stress changes on the fault over the next 30 years suggests that this seismic gap of the central segment will be shortened gradually. The segments that surrounding Eryuan, between Dali and north of Dadoumen and south of Yuanjiang may suffer high seismic risk due to the increased stress level.
Key words: Red River fault    Postseismic viscoelastic relaxation    Interseismic tectonic loading    Coulomb stress change    Seismicity    
0 引言

红河断裂带位于青藏高原东南缘,从NW到SE跨越整个云南地区并向南延伸至越南北部,是一条巨型的右旋走滑断裂,全长超过900 km,在中国云南境内长约600 km(Allen et al., 1984; 张建国等, 2009a).前人根据断裂的几何结构和地震活动性等差异,将中国境内的红河断裂带从NW到SE依次分成了北(洱源—苴力断裂段)、中(苴力—大斗门断裂段)、南(春元—河口断裂段)3段(虢顺民等, 2001; 张建国等, 2009a; 李西等, 2016).

红河断裂带地处川滇菱形块体的西南边界,是扬子地块与印支地块间的边界断裂(Tapponnier et al., 1990),两侧强烈的板块运动、周围极其复杂的构造背景和频繁的地震活动,使红河断裂带成为板块运动和地震学研究的重要窗口和热点区域.长达数百年的地震记载历史表明,红河断裂带北段地震频发,5级以上地震有几百次,有记载以来的强震与大地震集中发生在弥渡及其以北的北段,共发生过8次MS≥6.5的地震(张建国等, 2009a).断裂带中南段,历史上未记载或记录到强震和大地震,有现代地震仪器记录以来,存在长期缺乏强震的地震空区,中段即使小震也鲜有发生(图 1).但野外地质地貌考察和古地震研究发现,红河断裂带中段具有明确的更新世和全新世以来的强震活动证据(Allen et al., 1984)及M≥7古地震遗迹(Replumaz et al., 2001),说明该断裂段具备发生M≥7大地震的能力.根据地震空区理论,位于已发生大地震震源区之间的、尚未发生大地震的地震空区,是未来大地震最有可能发生的地段(Sykes, 1971).现今红河断裂带中段表现出来的地震相对平静现象,可能正处在大震孕育期,只是它的复发周期比较长(Allen et al., 1984).

图 1 1970年以来红河断裂带M≥2.0的地震震中分布 Fig. 1 Seismicity (M≥2.0) in the study region since 1970

分析断层上库仑应力演化是研究断层上应力状态和判定区域强震趋势的重要手段.地震学家们广泛采用断层上库仑应力的变化来分析断层变得更危险还是更安全.在过去的20多年时间里,研究断层间相互作用的库仑应力变化的概念发展起来并得到广泛应用(King et al., 1994; Harris and Simpson, 1998; Marsan, 2003; Ma et al., 2005; Toda et al., 2008),岩浆等流体的侵入也会引起附近断层上应力变化,进而加快或延迟未来地震的发生(Thatcher and Savage, 1982; Hill et al., 2002).近年来,根据断层上库仑应力演化过程对破裂带进行区域地震活动性分析, 国内外开展了广泛研究,并取得了大量成果(Deng and Sykes, 1997; Tabrez, 2008; Mouyen et al., 2010; Shan et al., 2011, 2013; Verdecchia and Carena, 2015; Shao et al., 2016).

地震的发生会引起区域应力场调整,包括弹性的同震应力阶变和由于地球介质的震后黏滞松弛造成的震后应力场变化.地震造成的同震应力阶变对短距离地震触发和可能的余震分布估计等有重要意义(Reasenberg and Simpson, 1992; Toda et al., 1998; 石耀霖, 2001).同震应力场的计算是基于地球模型为完全弹性体的假定,但实际上,岩石圈下地壳和上地幔并非完全弹性,而是黏弹性的.由于热的下地壳和上地幔的流变性质,在震后的数年至数百年时间里(达到稳态所用时间取决于地层的黏滞系数)会将下地壳和上地幔中的同震应力场变化向上传递到上地壳的孕震层中,传递上来的应力进一步作用于周围断层(Pollitz, 1992; Parsons and Dreger, 2000; Chery et al., 2001; Marsan and Bean, 2003; Lorenzo-Martín et al., 2006; Freed et al., 2007),对周围断层加载还是卸载这与发震断层和接收断层的几何特征、滑动性质等有关(Ali et al., 2008).大量研究认为震后黏滞松弛造成的应力转移在区域应力演化方面起着非常重要的作用(Deng et al., 1999; Freed and Lin, 2001; Pollitz et al., 2003; Lorenzo-Martin et al., 2006; Smith and Sandwell, 2006; Freed et al., 2007; Ali et al., 2008).震间构造应力的加载效应长期作用于断层,对断层上应力演化的贡献也不可忽略(Scholz, 1998; Nalbant et al., 2002),甚至对地震的发生起了主导作用,例如Shao等(2016)在计算鲜水河断裂带上强震发生前的累计库仑应力变化对强震的触发作用,发现对于该断裂带上的大多数地震来说震间构造应力加载起了主要作用.因此,在计算数百年尺度的应力场变化时同震应力阶变、震后黏滞松弛效应和震间构造应力加载作用都不容忽视.

红河断裂带邻区的安宁河—则木河—小江断裂带、楚雄—建水断裂、丽江—小金河断裂、无量山断裂、澜沧江断裂等诸多活动断裂上地震频繁且强烈,发生过多个灾难性的大地震(图 2).周围强震的发生对红河断裂带有何影响?该断裂带上的应力是如何演化的?这些信息对红河断裂带上未来地震的灾害评估有什么帮助?本文基于黏弹性地球介质模型,考虑同震应力阶变、震后黏滞松弛和震间构造加载三方面的综合效应,计算1833年以来红河断裂带上库仑应力变化随时间的演化,并分析其与现今地震活动的相应关系,最后对未来三十年红河断裂带的地震危险性趋势进行判断.

图 2 1833年以来红河断裂带(M≥6.0)及其周围(M≥7.0)的强震震源机制解 图中灰色震源球为红河断裂带上的地震震源机制解,红色震源球为周围断裂带上发生的强震震源机制解.以红色字体标出的时间和地点为死亡人数在300人以上的破坏性大地震(陈尚平,2008).地震目录资料来自《中国历史强震目录》(国家地震局震害防御司, 1995),中国近代地震目录(中国地震局震害防御司, 1999)和由中国地震台网中心提供的1991年1月1日至今的《中国大震速报目录》.绿色曲线为中国境内的红河断裂带,黑色曲线为云南地区主要活动断裂带,资料来自邓起东等(2002). Fig. 2 Focal mechanism of strong earthquakes on the Red River fault(M≥6.0) and surrounding faults (M≥7.0) since 1833 The gray beach balls denote the focal mechanism of events on the Red River fault, and the red beach balls denote the focal mechanism of events on surrounding faults. The earthquakes which resulted in a death toll more than 300 (Chen, 2008) are marked in red with origin time and loaction. The Data refers to "The catalogue of Chinese historical strong earthquakes" (Department of Earthquake Disaster Prevention, State Seismological Bureau, 1995), "The catalogue of Chinese modern earthquakes" (Department of Earthquake Disaster Prevention, China Earthquake Administration, 1999) and "Fast Reports on Chinese Large Earthquakes" supplied by China Earthquake Networks Center. The green lines indicate the Red River fault zone within China, and the black lines indicate the main active faults within Yunnan (Deng et al., 2002).
1 研究方法和数据 1.1 研究方法

本研究主要基于库仑破裂应力的概念(公式(1)),计算断层上库仑应力的变化量(公式(2)),

(1)

其中,CFS为库仑应力,τ为断层面上的剪应力,σn为断层面上的正应力,拉伸为正,压缩为负.μ为岩石摩擦系数,p为流体孔隙压力(Jaeger and Cook, 1979; Schloz, 2002).

静态库仑应力的变化量的表达式为

(2)

其中,ΔCFS为库仑应力变化量,Δτ为断层面上的剪应力变化量,Δσn为断层面上的正应力变化量,方向同上述的正应力方向的定义.计算得到的ΔCFS值为负时,表示断层变得更稳定,反之则断层趋于失稳.μ′为有效摩擦系数,不同研究中取值有所不同,但不会改变库仑应力变化的整体分布形态(孙玉军等, 2013; 万永革等, 2009),本研究取强震引起的库仑应力变化研究中的常用值0.4(Stein et al., 1992; King et al., 1994).

本文采用汪荣江研究员开发的PSGRN/PSCMP计算程序(Wang et al., 2006),该程序基于分层、重力作用下的黏弹性地球模型可用来研究强震引起的同震及震后形变场和应力场变化.震间构造应力加载作用参考前人所用的方法(Tabrez et al., 2008; 徐晶等, 2013; Verdecchia and Carena, 2015; Shao et al., 2016),基于负位错理论(Savage and Burford, 1973; Savage, 1983; Matsu′ura et al., 1986)进行计算,负位错量用断层位移滑动亏损速率和时间间隔的乘积表示,将负位错量加载到相应断层段,用与计算同震相同的方法来计算震间构造加载作用造成的应力场变化.

1.2 地层介质模型参数

地球介质模型采用分层半无限空间黏弹性模型,分为弹性的上地壳、黏弹性的下地壳和上地幔.本研究参考王夫运等(2014)在滇中南地区布设的跨越红河断裂带的主动源高分辨宽角反射/折射探测剖面结果,给出了红河断裂带地区的壳幔分层结构和P波速度(图 3).S波速度由经验公式VS=VP/3得到.地层介质密度参考朱思林等(1994)按照Nafe-drake密度-波速经验转换公式(Talwani et al., 1959)结合滇西人工地震测深资料给出的介质密度和P波速度之间的关系式给出.剪切模量G的表达式(Aki and Richards, 2002)为

图 3 分层黏弹性地球模型 其中黑色实线为VP随深度变化曲线,黑色虚线为剪切模量G随深度变化曲线.ρ表示地层介质密度, 单位为kg·m-3.η1η2分别为地层的瞬时和稳态黏滞系数,单位为Pa·s. Fig. 3 Viscoelastic stratified earth model Stratified earth model comprised of upper crust, viscoelastic lower crust and viscoelastic mantle. The solid line indicates the velocity of P wave, while the dashed line indicates the share modulus. ρ is the rock density, and its unit is kg·m-3. η1η2 are transient viscosity and steady-state viscosity, respectively,and their unit is Pa·s.

(3)

本模型中弹性体厚度为22 km,黏弹性的下地壳采用Burgers体模型(邵志刚等, 2007Ryder et al., 2011),根据石耀霖和曹建玲(2008)计算的云南南部下地壳的稳态黏滞系数为1.0×1020 Pa·s,瞬时与稳态黏滞系数一般相差1~1.5个数量级(Ryder et al., 2011),本研究下地壳瞬时黏滞系数取1.0×1019 Pa·s.地幔黏滞系数比下地壳大很多,为1.0×1021 Pa·s(石耀霖和曹建玲, 2008).川滇地区的大多数地震发生在5~25 km深度,这部分属于脆性的中上地壳,有条件形成地震活动带的孕震区(王椿镛等, 2002),本研究选取地下10 km作为计算断层面上库仑应力变化的深度.

1.3 红河断裂带及周围的强震目录和同震位错参数

红河断裂带北区段有较全地震记载以来,大致可分为3个活动时段,分别为1489—1518年,1566—1663年和1839—1943年(张建国等, 2009a),中区段从1900年至1964年云南有较完全的破坏性地震记载以来仅有一次发生于1929年2月9日的5.25级地震.而1833年9月6日发生在小江断裂西支的M8.0嵩明大地震,距红河断裂带也就百余公里,震级之大对红河断裂带必定会产生一定影响,因此,本文主要计算包括北区段第三个活动时段在内自1833年以来红河断裂带上的库仑应力变化.自1833年至今,红河断裂带上共发生M≥6.0强震5次,附近断裂带(50 km范围内)M≥6.0地震和(300 km范围内)M≥7.0地震共20次,因此参与库仑应力计算的地震事件共25个(图 2).

计算强震造成的同震和震后影响,须已知每一地震的同震位错量、滑动角和断层几何特征(长度、宽度、走向、倾角)等信息.地震目录给出了每一地震的发震时刻、位置和震级.其它参数主要参考地质考察的文献资料结果,大震的地质考察结果较为详尽,例如1833年M8.0嵩明地震(陈睿和李坪, 1988; 俞维贤等, 1996).对于没有地质考察资料的事件,断层破裂长度、宽度和平均位错量则根据全球地震定比定律(Wells and Coppersmith, 1994)确定,断层面的走向、倾角、深度,则参考美国哈佛大学Global CMT(http://www.globalcmt.org/CMTsearch.html)或近年来附近小震震源机制解的结果.最终给出每一地震事件的长度、宽度、走向、倾角、沿走向和倾向的滑动量(表 1).

表 1 红河断裂带及其周围的强震同震位错模型参数 Table 1 Coseismic slip parameters of strong earthquakes along the Red River fault and surrounding areas
1.4 断层滑动亏损速率

红河断裂带由NW向SE方向伸展,逐渐由NNW变为NWW方向,根据断层的几何展布特征,在计算震间构造应力加载时,将断层分为以下4段考虑,同时,各个子断层也是用来作为计算库仑应力变化的接收断层,计算深度为10 km,具体参数如表 2所示.断层倾角较陡,倾角70°以上的占一半以上(王绍晋等, 2010),计算中采用断层平均倾角70°.根据刘耀辉等(2015)利用1999—2013年GPS反演结果,红河断裂带北段16~20 km以下逐步转为蠕滑,中段和南段在12~20 km以下逐步转为蠕滑,所以取断层最大闭锁深度为20 km,根据倾角换算成断层闭锁宽度最大为21 km,与图 3中给出的弹性体厚度22 km基本吻合.红河断裂带在6~12 km深度不同分段的滑动亏损速率分别为北段2.1 mm·a-1,中段1.5 mm·a-1,南段1.8 mm·a-1.断层滑动亏损速率和时间间隔的乘积即为计算震间所用的负位错量,只是符号相反.

表 2 负位错反演参数(王绍晋等, 2010; 刘耀辉等, 2015) Table 2 Parameters of "virtual negtive displacements" along the Red River Fault(Wang, et al., 2010; Liu, et al., 2015)
2 库仑应力演化 2.1 同震库仑应力变化

本文计算了由于红河断裂带及其周围强震的发生对红河断裂带造成的累积同震库仑应力变化随时间的演化,如图 4所示.由图中可以看出,断层上同震应力场的整体分布主要受1833年嵩明M8.0地震控制,该地震也是研究区内近200年来发生的最大的地震.1833年嵩明地震造成的同震库仑应力变化,对南北两段有不同程度的加载,这两个区段也是近200年以来红河断裂带上地震活动相对活跃的地区;同时形成了断裂带中区段的应力影区,这一应力影区范围与中区段长期存在的地震空区基本吻合,从一定程度上可以解释现今中区段的地震空区现象.

图 4 自1833年以来,由于强震的同震应力阶变在红河断裂带上造成的累积库仑应力变化 各子图分别为(a)1839a洱源地震、(b)1839b洱源地震、(c)1901邓川地震、(d)1925a大理地震、(e)1925b弥渡地震、(f)1950年、(g)2000年、(h)2015年之前,所有已发生强震在红河断裂带上造成的累积同震库仑应力变化.接收断层参数如表 2所示,计算深度为10 km,断层有效摩擦系数取0.4.红色五角星表示1833年嵩明地震震中位置;黑色五角星表示当前时刻地震的震中位置. Fig. 4 Cumulative coseismic Coulomb stress changes on the Red River fault from all preceding events since 1833 Cumulative coseismic Coulomb stress changes resolved on the Red River fault immediately before (a) 1839a Eryuan earthquake, (b) 1839b Eryuan earthquake, (c) 1901 Dengchuan earthquake, (d) 1925a Dali earthquake, (e) 1925b Midu earthquake, (f) the year of 1950, (g) the year of 2000, (h) the year of 2015 from all preceding events since 1833. The parameters of the receiver fault are listed in Table 2. Stresses are resolved at a detpth of 10 km with a friction of 0.4. The red pentagram indicates the epicenter of 1833 Songming earthquake, and the black pentagrams indicate the epicenter of the earthquake marked in each subplot.

1833年嵩明M8.0地震的发生,对红河断裂带中段(约23.5°—25°N)卸载明显,最大同震库仑应力卸载量为0.07 MPa;对北段影响不明显,最大同震库仑应力加载为0.005 MPa;对南段最大加载量达0.038 MPa(图 4a).1839年2月23日洱源M7.0地震,与上次地震(1839a)同在洱源且震源机制解相同,相隔仅16天相继破裂,2月16日的M7.0洱源地震断层破裂已经释放了震源区附近的大部分应力,因此1839b地震震源位于应力影区(图 4b).1901年以前发生的地震对南段有进一步加载,在红河州附近超过0.1 MPa(图 4c).1925a大理地震前,在前面所有地震即1833年嵩明地震、1839a洱源地震、1839b洱源地震和1901年邓川地震的同震累积作用下,震源位置附近同震库仑应力累积加载量达到0.01 MPa(达到地震触发阈值0.01 MPa)(图 4d).1925b弥渡地震以后,红河断裂带上没有再发生强震,1950年(图 4f)、2000(图 4g)和2015(图 4h)的库仑应力变化主要是在周围断裂带上强震作用下的结果,但周围断裂带上离红河断裂带较近的地震在M6.0~6.5之间,M7.0地震发生在红河断裂带西南200余公里,所以红河断裂带上的同震库仑应力变化量级并没有显著变化,应力状态并未受到周围断裂带上地震的显著影响.

2.2 震后库仑应力变化

图 5(a—h)显示了由于震后黏弹松弛作用在红河断裂带上造成的累积库仑应力变化随时间的演化过程.图 5图 4比较可以看出,在震后黏滞松弛和同震完全弹性两种截然不同的物理机制作用下,1833年嵩明M8.0地震在震后应力调整中仍起到了主导作用,造成了断裂带南北段库仑应力的增加区和断裂带中段库仑应力变化的低值区.

图 5 自1833年以来,强震的震后黏滞松弛作用在红河断裂带上造成的累积库仑应力变化 计算的震后累积库仑应力変化,其它同图 4. Fig. 5 Cumulative postseismic Coulomb stress changes. See Fig. 4 for more details

图 6显示了强震的同震应力阶变和震后黏滞松弛效应共同作用下在红河断裂带上造成的累积库仑应力变化.正如上面讨论的,在同震效应基础上,震后黏滞松弛作用进一步增加了库仑应力变化的幅度.例如,1925a大理地震前,与累积同震库仑应力变化(约0.01 MPa)(图 4d)相比,震源位置附近同震和震后库仑应力共同累积加载量达到0.02 MPa(图 6d).1901年地震前,南段0.1 MPa量级的库仑应力变化增加区域从元江延伸到红河(图 6c6f).截止至2015年,中段的应力影区,与同震库仑应力最大卸载值(-0.16 MPa)(图 4h)相比,幅度增加,最大卸载值达-0.27 MPa(图 6h).

图 6 自1833年以来,红河断裂带上强震造成的库仑应力变化(同震+震后) 计算(同震+震后)效应,其它同图 4 Fig. 6 Cumulative combined coseismic and postseismic Coulomb stress changes, see Fig. 4 for more details
2.3 震间库仑应力变化

基于负位错理论,利用分段断层负位错参数(表 2)计算了震间构造长期持续加载下,红河断裂带上库仑应力年变化速率(图 7).计算结果表明,约600 km长的红河断裂带均得到加载,25°N以北的断层段应力加载速率最大(>1.5 kPa·a-1),苴力至大斗门地区(24.2°N—25°N)以南加载速率为1.0~1.5 kPa·a-1,个旧以南和中段向西南凸出的段落加载速率最小(<1 kPa·a-1).与地震造成的应力场变化相比,震间构造加载年变化速率虽然量级较小,但它的作用是长期且持续的,在地震周期的应力调整中不可轻易忽视.

图 7 红河断裂带上震间构造加载造成的库仑应力年变化速率 分段断层参数如表 2所示.计算库仑应力变化的深度为10 km,断层有效摩擦系数取0.4. Fig. 7 The rate of Coulomb stress changes of the Red River fault due to interseismic tectonic loading The parameters of the receiver fault are listed in Table 2. Stresses are resolved at a depth of 10 km with a friction of 0.4.

假设震间构造加载在历史时间尺度上是均匀的,则库仑应力变化速率与时间间隔的乘积即可得到相应时间段内的库仑应力变化量.本文计算了震间库仑应力变化随时间的演化过程(图 8).图 8(a—h)展示了近180年震间构造应力加载是如何随时间演化并作用于断层的.地震发生前震间构造运动对不同断层段起到了不同程度的加载作用,构造应力持续不断地加载,经过近180年的加载作用,到2015年在北段达到了0.1~1 MPa的量级范围.事实上,震间构造运动在上地壳中积累的应力最终大部分会以地震的形式释放掉,如果将图 8中断层上库仑应力变化量的正负符号反过来,则各子图可以理解为,假设断层发生完全破裂释放应力的过程.

图 8 自1833年以来,红河断裂带上震间构造加载引起的累积库仑应力变化随时间的演化过程 Fig. 8 Cumulative interseismic Coulomb stress changes, see Fig. 4 for more details
2.4 红河断裂带上库仑应力变化随时间的演化过程(同震+震后+震间)

上述三节分别计算了弹性同震、黏弹性震后和震间构造加载造成的红河断裂带上库仑应力变化情况,本节综合同震、震后和震间三方面效应,给出了自1833年以来红河断裂带上累积库仑应力变化随时间的演化过程(图 9).红河断裂带上累积库仑应力变化显示出两端加载、中段卸载的整体布局,这与云南地区小地震重定位结果所显示的红河断裂带上小震多分布在断裂带的南北两端,而在中段小震相对稀少(张广伟, 2016)的分布特点也是一致的.与强震造成的累积库仑应力变化水平(图 6)相比,震间构造加载使红河断裂带南北区段库仑应力增加区域进一步增强,中区段的应力影区范围自1925年以后有所减小.截至2015年,除去北区段的大理附近区域,对断裂带南北两段加载明显,达到0.1~1.0 MPa的水平,断裂带中段应力影区范围为23.7°N—24.4°N,仍处于应力水平较低的稳定状态(<-0.14 MPa),这与现今观测到的中段小震空区基本吻合.在不考虑其它条件下,综合中段应力影区的库仑应力卸载量和震间累积速率,影区南端最开始消失,需要约4年时间,影区中部消失最慢,需要约213年,平均消失时间需约96年.

图 9 自1833年以来,红河断裂带上库仑应力变化演化过程(同震+震后+震间) Fig. 9 Cumulative combined coseismic, postseismic and interseismic Coulomb stress changes, see Fig. 4 for more details

2015年以来在红河断裂带上的个旧附近共发生5次M≥3.0地震,即2016年5月4日在个旧连续发生M4.6、M3.8、M4.8、M3.6地震,2016年6月8日又发生M3.0地震,本次地震是5月4日地震的余震,5次地震均位于红河断裂带的累积库仑应力变化高值区(图 10).假设2015年以后红河断裂带及周围未发生强震,则在震后黏滞松弛和震间构造持续加载作用下,截至2050年,红河断裂带上的累积库仑应力变化分布,如图 10所示.结果表明,自2015年经过35年的应力积累,断裂带中段应力影区依然存在,洱源附近、大理至大斗门以北地区、元江以南地区应力增加超过0.1 MPa,发震危险性增加.

图 10 截至2050年红河断裂带上累积(同震+震后+震间)库仑应力变化分布 绿色空心圆为发生在个旧的共5个M≥3.0地震的震中位置. Fig. 10 Cumulative Coulomb stress changes in 2050 due to coseismic and postsesimic and interseismic deformation Green emtpty circles indicate the epicenters of the earthquakes (M≥3.0) in Gejiu
3 讨论

由于下地壳和地幔采用黏弹流变模型,地层的黏滞系数对震后的应力场演化可能会有一定影响.在前面计算中主要参考了石耀霖和曹建玲(2008)的结果,但由于研究区震后形变观测数据缺少,很难对下地壳和地幔的黏滞性进行很好的约束,为了保证计算结果的可靠性,本文选取了另一组黏滞性参数来进行对照和检验.选取值比前面相应地低一个数量级,即下地壳的稳态黏滞系数为1.0×1019 Pa·s,瞬时黏滞系数为1.0×1018 Pa·s,地幔黏滞系数为1.0×1020 Pa·s,计算了黏滞性系数对断层上震后累计库仑应力变化的影响,结果如图 11所示.

图 11 自1833年以来,强震的震后黏滞松弛作用在红河断裂带上造成的累积库仑应力变化 下地壳的稳态黏滞系数为1.0×1019 Pa·s,瞬时黏滞系数为1.0×1018 Pa·s,地幔黏滞系数为1.0×1020 Pa·s,其他同图 5. Fig. 11 Same as Fig. 5, but for different viscosity The stable viscous coefficient of the lower curst is 1.0×1019 Pa·s, the transient, 1.0×1018 Pa·s. The viscous coefficient of the mantle is 1.0×1020 Pa·s.

由于下地壳和上地幔中的流变性质,震后会将下地壳和上地幔中的同震应力变化向上传递到上地壳的孕震层中,应力传递速度主要与地层的黏滞系数有关,黏滞系数越低,则传递越快(Shan et al., 2013),再次达到稳态所需的时间越短,反之,黏滞系数越高,则传递越慢,再次达到稳态需要更长的时间.图 11图 5相比,可以发现,在震后相同时间间隔内,二者库仑应力变化整体正负分布相同,但在低的黏滞系数作用下,图 11中的应力变化幅度更大,即库仑应力变化正值更大,负值更小.图 12为低黏滞系数作用下,同震、震后和震间综合作用在断裂带上造成的累积库仑应力变化结果,与图 9相比,除中段地震空区范围有所变大(23.7°N—24.7°N)外,黏滞系数差异导致的震后库仑应力变化的差异,并未对断层上累积库仑应力变化的整体分布和量级等造成显著影响.

图 12 自1833年以来,红河断裂带上库仑应力变化演化过程(同震+震后+震间) 下地壳的稳态黏滞系数为1.0×1019 Pa·s,瞬时黏滞系数为1.0×1018 Pa·s,地幔黏滞系数为1.0×1020 Pa·s,其他同图 9. Fig. 12 Same as Fig. 9, but for different viscosity The stable viscous coefficient of the lower curst is 1.0×1019 Pa·s, the transient, 1.0×1018 Pa·s. The viscous coefficient of the mantle is 1.0×1020 Pa·s.

本研究中仅考虑M≥6.0地震的影响,并未考虑数目更多的小地震,这样可能对计算结果带来一定的不确定性.然而小震对计算结果可能造成约10%的影响,并不影响库仑应力变化的趋势分布(Mouyen et al., 2010).

自1925年3月17日弥渡地震前,断裂带的北区段25.5°N—26°N之间大理附近开始落入应力影区,但近年来大理附近的小震活动比较频繁,两者看似矛盾,一个可能的原因是由于附近历史地震可考资料有限(比如1925年3月16日大理地震参数的不确定性)导致附近的计算结果不准确,另外本文未考虑1833年以前断裂带的应力场状态,这些因素也是以后需要探讨的.

4 结论

本文计算了1833年以来红河断裂带上的库仑应力变化随时间的演化,得出主要结论如下:

(1) 在历史强震的同震、震后效应和震间构造加载综合作用下,近180年红河断裂带中区段长期处于应力影区,这与中区段长期缺乏强震、甚至小震也很少发生的地震平静现象相一致.而南、北两段的库仑应力累积了正的库仑应力,与现今地震活动主要发生在南北段的观测结果也是一致的.

(2) 红河断裂带的发育演化,尽管与欧亚板块和印度板块的碰撞、青藏高原的动力演化密切相关(Molnar and Tapponnier, 1975; Leloup et al., 1995; Tapponnier et al., 2001),但断裂带中段与南北两段截然相反的地震活动性差异,与近百年来强震的发生和构造运动对断裂带不同段落造成的不同程度的应力加载或卸载作用是不容忽视的,尤其1833年嵩明M8.0大地震对近180年来研究区库仑应力变化的整体布局起了主导作用.

(3) 假设2015年以后红河断裂带及周围不再发生强震,则在震后黏滞松弛和震间构造持续作用下,截至2050年,断裂带中段应力影区依然存在,但范围在缩小;洱源附近、大理至大斗门以北地区、元江以南地区应力增加超过0.1 MPa,可能仍是地震潜在危险段.综合中段应力影区的库仑应力卸载量和震间累积速率,影区完全消失最长需要约213年.

(4) 计算中存在一些不确定因素和基本假设,如十九世纪没有仪器记载的历史地震的震源机制,断层摩擦系数,同震位错模型采用均一的平均位错等,对计算结果可能会产生一定影响,这是未来研究中需要进一步考虑的.尽管如此,该研究通过计算地震孕育和发生过程中同震、震后和震间三个重要阶段的库仑应力变化,对分析红河断裂带上的现今应力演化,进而科学评估未来地震潜在危险性,提供了参考依据.

致谢

感谢汪荣江老师和胡岩老师对本研究提出的建设性建议.

参考文献
Aki K, Richards P G. 2002. Quantitative Seismology, second edition.University Science Books, California, USA.
Ali S T, Freed A M, Calais E, et al. 2008. Coulomb stress evolution in Northeastern Caribbean over the past 250 years due to coseismic, postseismic and interseismic deformation. Geophysical Journal International, 174(3): 904-918. DOI:10.1111/gji.2008.174.issue-3
Allen C R, Gillespie A R, Yuan H, et al. 1984. Red River and associated faults, Yunnan Province, China:Quaternary geology, slip rates, and seismic hazard. GSA Bulletin, 95(6): 686-700. DOI:10.1130/0016-7606(1984)95<686:RRAAFY>2.0.CO;2
Chen R, Li P. 1988. Slip rates and earthquake recurrence intervals of the western branch of the Xiaojiang fault zone. Seismology and Geology, 10(2): 1-13.
Chen S P. 2008. Anthology of Historical Destructive Earthquakes in Modern China (1900-1950). Beijing: China Ocean Press.
Chery J. 2001. A physical basis for time clustering of large earthquakes. Bulletin of the Seismological Society of America, 91(6): 1685-1693. DOI:10.1785/0120000298
Deng J S, Sykes L R. 1997. Evolution of the stress field in southern California and triggering of moderate-size earthquakes:A 200-year perspective. Journal of Geophysical Research:Solid Earth, 102(B5): 9859-9886. DOI:10.1029/96JB03897
Deng J S, Hudnut K, Gurnis M, et al. 1999. Stress loading from viscous flow in the lower crust and triggering of aftershocks following the 1994 Northridge, California, earthquake. Geophysical Research Letters, 26(21): 3209-3212. DOI:10.1029/1999GL010496
Deng Q D, Zhang P Z, Ran Y K, et al. 2002. Basic characteristics of active tectonics of China. Science in China Series D:Earth Sciences, 46(4): 356-372.
Department of Earthquake Disaster Prevention, State Seismological Bureau. 1995. The Catalogue of Chinese Historical Strong Earthquakes. Beijing: Seismological Press.
Department of Earthquake Disaster Prevention, China Earthquake Administration. 1999. The Catalogue of Chinese Modern Earthquakes (1912~1990, MS ≥ 4. Beijing: China Science and Technology Press.
Freed A M, Lin J. 2001. Delayed triggering of the 1999 Hector Mine earthquake by viscoelastic stress transfer. Nature, 411(6834): 180-183. DOI:10.1038/35075548
Freed A M, Ali S T, Bürgmann R. 2007. Evolution of stress in Southern California for the past 200 years from coseismic, postseismic and interseismic stress changes. Geophysical Journal International, 169(3): 1164-1179. DOI:10.1111/gji.2007.169.issue-3
Guo S M, Ji F J, Xiang H F, et al. 2001. The Honghe Active Fault Zone. Beijing: Ocean Publishing House.
Guo S M, Xu X W, Xiang H F, et al. 2002. Segmentation of earthquake rupture and earthquake prediction along the Longling-Lancang fault zone in the southwestern Yunnan Province. Seismology and Geology, 24(2): 133-144.
Han X M, Mao Y P. 1993. Seismic risk analysis to Shiping-Jianshui fault zone in the next 30 years. Journal of Seismological Research, 16(1): 52-59.
Hao P, Fu Z X, Tian Q J, et al. 2004. The tress triggering of large aftershocks following the 1970 Tonghai great earthquake (M=7.7). Earthquake, 24(2): 38-46.
Harris R A, Simpson R W. 1998. Suppression of large earthquakes by stress shadows:A comparison of Coulomb and rate-and-state failure. Journal of Geophysical Research:Atmospheres, 103(B10): 24439-24451. DOI:10.1029/98JB00793
Hill D P, Pollitz F, Newhall C. 2002. Earthquake-volcano interactions. Physics Today, 55(11): 41-47. DOI:10.1063/1.1535006
Jaeger J C, Cook N G W. 1979. Fundamentals of Rock Mechanics. (3rd ed). London: Chapman and Hall.
Ji S C, W Q, Sun S S, et al. 2008. Continental extrusion and seismicity in China. Acta Geologica Sinica, 82(12): 1644-1667.
King G C P, Stein R S, Lin J. 1994. Static stress changes and the triggering of earthquakes. Bulletin of the Seismological Society of America, 84(3): 935-953.
Leloup P H, Lacassin R, Tapponnier P, et al. 1995. The Ailao Shan-Red River shear zone (Yunnan, China), tertiary transform boundary of Indo-China. Tectonophysics, 251(1-4): 3-10. DOI:10.1016/0040-1951(95)00070-4
Li X, Ran Y K, Chen L Y, et al. 2016. The Holocene seismic evidence on southern segment of the Red River fault zone. Seismology and Geology, 38(3): 596-604.
Li X G, Ran Y K, Peng G, et al. 1986. Active faults and their relationship with seismic activities in Dali of Yunnan Province and its adjacent regions. Seismology and Geology, 8(4): 51-61.
Liu Y H, Li J P, Wang L W. 2015. Inversion study on locking and slip deficit of Red River fault zone. Engineering of Surveying and Mapping, 24(8): 20-22, 26.
Lorenzo-Martín F, Roth F, Wang R J. 2006. Elastic and inelastic triggering of earthquakes in the North Anatolian Fault zone. Tectonophysics, 424(3-4): 271-289. DOI:10.1016/j.tecto.2006.03.046
Ma K F, Chan C H, Stein R S. 2005. Response of seismicity to Coulomb stress triggers and shadows of the 1999 MW=7.6 Chi-Chi, Taiwan, earthquake. Journal of Geophysical Research:Solid Earth, 11(B5): B05S19. DOI:10.1029/2004JB003389
Matsu'Ura M, Jackson D D, Cheng A. 1986. Dislocation model for aseismic crustal deformation at Hollister, California. Journal of Geophysical Research Solid Earth, 91(B12): 12661-12674. DOI:10.1029/JB091iB12p12661
Marsan D, Bean C J. 2003. Seismicity response to stress perturbations, analysed for a world-wide catalogue. Geophysical Journal International, 154(1): 179-195. DOI:10.1046/j.1365-246X.2003.01963.x
Marsan D. 2003. Triggering of seismicity at short timescales following Californian earthquakes. Journal of Geophysical Research:Atmospheres, 108(B5): 2266. DOI:10.1029/2002JB001946
Molnar P, Tapponnier P. 1975. Cenozoic Tectonics of Asia:Effects of a Continental Collision:Features of recent continental tectonics in Asia can be interpreted as results of the India-Eurasia collision. Science, 189(4201): 419-426. DOI:10.1126/science.189.4201.419
Molnar P, Deng Q D. 1984. Faulting associated with large earthquakes and the average rate of deformation in central and eastern Asia. Journal of Geophysical Research:Solid Earth, 89(B7): 6203-6227. DOI:10.1029/JB089iB07p06203
Mouyen M, Cattin R, Masson F. 2010. Seismic cycle stress change in western Taiwan over the last 270 years. Geophysical Research Letters, 37(3): L03306. DOI:10.1029/2009GL042292
Mozaffari P, Xu L S, Wu Z L, et al. 1999. Moment tensor inversion of the November 6, 1988 MS7.6, Lancang-Gengma, China, earthquake using Long-Period body waves data. Acta Seismologica Sinica, 21(4): 344-353.
Nalbant S S, Mccloskey J, Steacy S, et al. 2002. Stress accumulation and increased seismic risk in eastern Turkey. Earth and Planetary Science Letters, 195(3-4): 291-298. DOI:10.1016/S0012-821X(01)00592-1
Paradisopoulou P M, Garlaouni C G, Jin X S, et al. 2007. Application of the stress evolutionary model along the Xiaojiang fault zone in Yunnan Province, Southeast China. Acta Geophysica, 55(4): 577-593.
Parsons T, Dreger D S. 2000. Static stress impact of the 1992 Landers earthquake sequence on nucleation and slip at the site of the 1999 Mw7.1 Hector Mine earthquake, southern California. Geophysical Research Letters, 27(13): 1949-1952. DOI:10.1029/1999GL011272
Pollitz F F. 1992. Postseismic relaxation theory on the spherical earth. Bulletin of the Seismological Society of America, 82(1): 422-453.
Pollitz F, Vergnolle M, Calais E. 2003. Fault interaction and stress triggering of twentieth century earthquakes in Mongolia. Journal of Geophysical Research:Solid Earth, 108(B10): 2503. DOI:10.1029/2002JB002375
Qin J Z, Ye J Q, Qian X D, et al. 2004. Source parameters of the Yao'an earthquakes in 2000. Chinese Journal of Geophysics, 46(5): 633-641.
Qin S L, Zhang J G, Liao L X. 2012. Analysis on focal mechanism and seismogenic structures of Yaoan earthquake with MS6.0, 2009. Inland Earthquake, 2009.
Reasenberg P A, Simpson R W. 1992. Response of regional seismicity to the static stress change produced by the Loma Prieta earthquake. Science, 255(5052): 1687-1690. DOI:10.1126/science.255.5052.1687
Replumaz A, Lacassin R, Tapponnier P, et al. 2001. Large river offsets and Plio-Quaternary dextral slip rate on the Red River fault (Yunnan, China). Journal of Geophysical Research:Solid Earth, 106(B1): 819-836. DOI:10.1029/2000JB900135
Ryder I, Bürgmann R, Pollitz F. 2011. Lower crustal relaxation beneath the Tibetan plateau and Qaidam basin following the 2001 Kokoxili earthquake. Geophysical Journal International, 187(2): 613-630. DOI:10.1111/gji.2011.187.issue-2
Savage J C, Burford R O. 1973. Geodetic determination of relative plate motion in central California. Journal of Geophysical Research, 78(5): 832-845. DOI:10.1029/JB078i005p00832
Savage J C. 1983. A dislocation model of strain accumulation and release at a subduction zone. Journal of Geophysical Research:Solid Earth, 88(B6): 4984-4996. DOI:10.1029/JB088iB06p04984
Scholz C H. 1998. Earthquakes and friction laws. Nature, 391(6662): 37-42. DOI:10.1038/34097
Scholz C H. 2002. The Mechanics of Earthquakes and Faulting. Cambridge: Cambridge University Press.
Shan B, Xiong X, Zheng Y, et al. 2011. The co-seismic Coulomb stress change and expected seismicity rate caused by 14 April 2010 Ms=7.1 Yushu, China, earthquake. Tectonophysics, 510(3-4): 345-353. DOI:10.1016/j.tecto.2011.08.003
Shan B, Xiong X, Wang R J, et al. 2013. Coulomb stress evolution along Xianshuihe-Xiaojiang fault system since 1713 and its interaction with Wenchuan earthquake, May 12, 2008. Earth and Planetary Science Letters, 377-378: 199-210. DOI:10.1016/j.epsl.2013.06.044
Shao Z G, Fu R S, Xue T X, et al. 2007. Simulating postseismic viscoelastic deformation based on Burgers model. Journal of Geodesy and Geodynamics, 27(5): 31-37.
Shao Z G, Xu J, Ma H S, et al. 2016. Coulomb stress evolution over the past 200years and seismic hazard along the xianshuihe fault zone of Sichuan, China. Tectonophysics, 670: 48-65. DOI:10.1016/j.tecto.2015.12.018
Shi Y L. 2001. Stress triggers and stress shadows:How to apply these concepts to earthquake prediction. Earthquake, 21(3): 1-7.
Shi Y L, Cao J L. 2008. Effective viscosity of China continental lithosphere. Earth Science Frontiers, 15(3): 82-95. DOI:10.1016/S1872-5791(08)60064-0
Smith B R, Sandwell D T. 2006. A model of the earthquake cycle along the San Andreas Fault System for the past 1000 years. Journal of Geophysical Research Solid Earth, 111: B01405. DOI:10.1029/2005JB003703
Stein R S, King G C P, Lin J. 1992. Change in failure stress on the southern San Andreas fault system caused by the 1992 magnitude=7.4 Landers earthquake. Science, 258(5086): 1328-1332. DOI:10.1126/science.258.5086.1328
Sun Y J, Dong S W, Fan T Y, et al. 2013. The effect of Tohoku MW9.0 earthquake on the near-field seismic activity from the coseismic and postseismic deformation. Progress in Geophysics, 28(3): 1131-1139.
Sykes L R. 1971. Aftershock zones of great earthquakes, seismicity gaps, and earthquake prediction for Alaska and the Aleutians. Journal of Geophysical Research, 76(32): 8021-8041. DOI:10.1029/JB076i032p08021
Tabrez A S, Freed A M, Calais E, et al. 2008. Coulomb stress evolution in Northeastern Caribbean over the past 250 years due to coseismic, postseismic and interseismic deformation. Geophysical Journal International, 174(3): 904-918. DOI:10.1111/gji.2008.174.issue-3
Talwani M, Worzel J L, Landisman M. 1959. Rapid gravity computations for two-dimensional bodies with application to the Mendocino submarine fracture zone. Journal of Geophysical Research, 64(1): 49-59. DOI:10.1029/JZ064i001p00049
Tapponnier P, Lacassin R, Leloup P H, et al. 1990. The Ailao Shan/Red River metamorphic belt:tertiary left-lateral shear between Indochina and South China. Nature, 343(1): 431-437.
Tapponnier P, Xu Z Q, Roger F, et al. 2001. Oblique stepwise rise and growth of the Tibet plateau. Science, 294(5547): 1671-1677. DOI:10.1126/science.105978
Thatcher W, Savage J C. 1982. Triggering of large earthquakes by magma-chamber inflation, Izu peninsula, Japan. Geology, 10(12): 637-640. DOI:10.1130/0091-7613(1982)10<637:TOLEBM>2.0.CO;2
Toda S, Stein R S, Reasenberg P A, et al. 1998. Stress transferred by the 1995 Mw=6.9 Kobe, Japan, shock-effect on aftershocks and future earthquake probabilities. Journal of Geophysical Research:Solid Earth, 103(B10): 24543-24565. DOI:10.1029/98JB00765
Toda S, Lin J, Meghraoui M, et al. 2008. 12 May 2008 M=7.9 Wenchuan, China, earthquake calculated to increase failure stress and seismicity rate on three major fault systems. Geophysical Research Letters, 35(17): L17305. DOI:10.1029/2008GL034903
Verdecchia A, Carena S. 2015. One hundred and fifty years of Coulomb stress history along the California-Nevada border, USA. Tectonics, 34(2): 213-231. DOI:10.1002/2014TC003746
Wan Y G, Shen Z K, Sheng S Z, et al. 2009. The influence of 2008 Wenchuan earthquake on surrounding faults. Acta Seismologica Sinica, 31(2): 128-139.
Wang C Y, Mooney W D, Wang X L, et al. 2002. Study on 3-D velocity structure of crust and upper mantle in Sichuan-Yunnan region, China. Acta Seismologica Sinica, 24(1): 1-16.
Wang F Y, Pan S Z, Liu L, et al. 2014. Wide angle seismic exploration of Yuxi-Lincang profile——The research of crustal structure of the Red River fault zone and southern Yunnan. Chinese Journal of Geophysics, 57(10): 3247-3258. DOI:10.6038/cjg20141013
Wang R J, Lorenzo-Martín F, Roth F. 2006. PSGRN/PSCMP-a new code for calculating co-and post-seismic deformation, geoid and gravity changes based on the viscoelastic-gravitational dislocation theory. Computers & Geosciences, 32(4): 527-541.
Wang S J, Long X F, Li Z H. 2002. Focal fault and stress field of Yao'an M6.5 earthquake on January 15, 2000. Earthquake Research in China, 18(1): 59-66.
Wang S J, Zhang J G, Yu Q K, et al. 2010. Focal mechanism of strong and medium-small earthquakes and modern tectonic stress-field of the Red-River fault zone. Journal of Seismological Research, 33(2): 200-207.
Wang X L, Liu J, Zhang G M, et al. 2006. Accurate location on the Yao'an earthquake sequence and the Yongsheng earthquake sequence. Earthquake Research in China, 21(3): 386-397.
Wang Y L, Qin Y X, Gu Z G, et al. 1994. The seismic rupture of the 1988 Lancang-Gengma earthquake with double main shock type and the recent tectonic movement in the seismic region. Journal of Seismological Research, 17(1): 46-55.
Wang Y R, Li Y L, Si S P. 2015. The geometrical characteristics of offset streams and the Shiping-Jianshui fault activity. Quaternary Sciences, 35(1): 100-108.
Well D L, Coppersmith K J. 1994. New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement. Bulletin of the Seismological Society of America, 84(4): 974-1002.
Wen X Z, Du F, Long F, et al. 2011. Tectonic dynamics and correlation of major earthquake sequences of the Xiaojiang and Qujiang-Shiping fault systems, Yunnan, China. Science China Earth Sciences, 54(10): 1563-1575. DOI:10.1007/s11430-011-4231-0
Xiang H F, Han Z J, Guo S M, et al. 2004. Processing about quantitative study of large-scale strike-slip movement on Red River fault zone. Advance in Earth Sciences, 19(S1): 56-59.
Xu J, Shao Z G, Ma H S, et al. 2013. Evolution of Coulomb stress and stress interaction among strong earthquakes along the Xianshuihe fault zone. Chinese Journal of Geophysics, 56(4): 1146-1158. DOI:10.6038/cjg20130410
Yan N L, Li H, Shen C Y. 2003. Analyses on gravity variation before and after the Lijiang earthquake based on a finite rectangular dislocation model. Acta Seismologica Sinica, 25(2): 172-181.
Yang X P, Chen L C, Ma W T, et al. 2008. Structural analysis and interpretation of the surface deformation associated with the Ning'er, Yunnan Province, China MS6.4 earthquake of June 3, 2007. Acta Seismologica Sinica, 30(2): 165-175.
Yu W X, Wang Y P, Song F M, et al. 1996. Survey study on surface rupture zone to Yunnan's Songming M=8 earthquake in 1833. Journal of Seismological Research, 19(4): 385-390.
Zhang G W. 2016. Research on earthquake relocation and b-value in the Yunnan area. Earthquake Research in China, 32(1): 54-62.
Zhang J G, Xie Y Q, Jin M P, et al. 2009a. Red River Fault Activity in China & Vietnam. Kunming: Science and Technology Publishing House.
Zhang J G, Liu L F, L X, et al. 2009b. Preliminary study on the seismotectonics of Yaoan & Dayao moderate strong earthquake area. Seismology and Geology, 31(3): 536-543.
Zhang S C, Liu B C. 1978. Seismic geological characteristics of Tonghai earthquake in 1970. Chinese Journal of Geology, 13(4): 323-335.
Zhou G Q, Su Y J, Wang S J. 2003. Study on the cause of the Yongsheng earthquake with M6.0 on October 27, 2001. Earthquake Research in China, 19(2): 166-174.
Zhou R Q, Yu W X, Gu Y S, et al. 1990. A study on rupture zone of the 1988 Gengma earthquake with magnitude 7.2 in Yunnan Province. Seismology and Geology, 12(4): 291-301.
Zhu C N, Zhou R Q. 1980. Basic analysis on the surface rupture of 1979 Pu'er, Yunnan M6.8 earthquake. Northwestern Seismological Journal, 2(2): 64-71, 97.
Zhu S L, Gan J S, Xu J S, et al. 1994. Three dimensional inversion of gravity anomalies in the western Yunnan. Crustal Deformation and Earthquake, 14(1): 1-10.
陈尚平. 2008. 中国近代(1900-1950)破坏性地震史料文图选集. 北京: 海洋出版社.
陈睿, 李玶. 1988. 小江西支断裂的滑动速率与强震重复周期. 地震地质, 10(2): 1–13.
邓起东, 张培震, 冉勇康, 等. 2002. 中国活动构造基本特征. 中国科学:地球科学, 32(12): 1020–1030, 1057.
邓起东, 张培震, 冉勇康, 等. 2003. 中国活动构造与地震活动. 地学前缘, 10(S1): 66–73.
国家地震局震害防御司. 1995. 中国历史强震目录(公元前23世纪~公元1911年). 北京: 地震出版社.
虢顺民, 计凤桔, 向宏发, 等. 2001. 红河活动断裂带. 北京: 海洋出版社.
虢顺民, 徐锡伟, 向宏发, 等. 2002. 龙陵-澜沧新生断裂带地震破裂分段与地震预测研究. 地震地质, 24(2): 133–144.
韩新民, 毛玉平. 1993. 石屏-建水断裂带未来三十年内七级以上大地震危险性分析. 地震研究, 16(1): 52–59.
郝平, 傅征祥, 田勤俭, 等. 2004. 1970年通海7. 7级大地震强余震触发.地震, 24(2): 38–46.
嵇少丞, 王茜, 孙圣思, 等. 2008. 亚洲大陆逃逸构造与现今中国地震活动. 地质学报, 82(12): 1643–1667.
李西, 冉勇康, 陈立春, 等. 2016. 红河断裂带南段全新世地震活动证据. 地震地质, 38(3): 596–604.
李祥根, 冉勇康, 彭贵, 等. 1986. 云南大理及其邻近地区活断层和地震活动特征. 地震地质, 8(4): 51–61.
刘耀辉, 李金平, 王刘伟. 2015. 红河断裂带闭锁程度和滑动亏损分布特征研究. 测绘工程, 24(8): 20–22, 26.
MozaffariP, 许力生, 吴忠良, 等. 1999. 用长周期体波数据反演1988年11月6日澜沧-耿马MS7.6地震的矩张量. 地震学报, 21(4): 344–353.
秦嘉政, 叶建庆, 钱晓东, 等. 2004. 2000年姚安地震的震源参数. 地球物理学报, 46(5): 633–641.
秦双龙, 张建国, 廖丽霞. 2012. 2009年云南姚安6.0级地震震源机制与发震构造的分析研究. 内陆地震, 26(1): 52–61.
邵志刚, 傅容珊, 薛霆虓, 等. 2007. 以Burgers体模型模拟震后粘弹性松弛效应. 大地测量与地球动力学, 27(5): 31–37.
石耀霖. 2001. 关于应力触发和应力影概念在地震预报中应用的一些思考. 地震, 21(3): 1–7.
石耀霖, 曹建玲. 2008. 中国大陆岩石圈等效粘滞系数的计算和讨论. 地学前缘, 15(3): 82–95.
孙玉军, 董树文, 范桃园, 等. 2013. 从同震和震后形变分析日本东北MW9.0级大地震对近场地震活动性的影响. 地球物理学进展, 28(3): 1131–1139. DOI:10.6038/pg20130304
万永革, 沈正康, 盛书中, 等. 2009. 2008年汶川大地震对周围断层的影响. 地震学报, 31(2): 128–139.
王椿镛, MooneyW D, 王溪莉, 等. 2002. 川滇地区地壳上地幔三维速度结构研究. 地震学报, 24(1): 1–16.
王夫运, 潘素珍, 刘兰, 等. 2014. 玉溪-临沧剖面宽角地震探测——红河断裂带及滇南地壳结构研究. 地球物理学报, 57(10): 3247–3258. DOI:10.6038/cjg20141013
王绍晋, 龙晓帆, 李忠华. 2002. 2000年1月15日姚安6.5级地震的震源断层与震源应力场. 中国地震, 18(1): 59–66.
王绍晋, 张建国, 余庆坤, 等. 2010. 红河断裂带的震源机制与现代构造应力场. 地震研究, 33(2): 200–207.
王新岭, 刘杰, 张国民, 等. 2006. 姚安地震序列与永胜地震序列的高精度定位. 中国地震, 21(3): 386–397.
王洋龙, 覃玉玺, 顾正光, 等. 1994. 1988年澜沧-耿马双主震型地震破裂与震区的现代构造运动. 地震研究, 17(1): 46–55.
王怡然, 李有利, 司苏沛. 2015. 石屏-建水断裂水平错断河流的几何特征与断层活动性研究. 第四纪研究, 35(1): 100–108.
闻学泽, 杜方, 龙锋, 等. 2011. 小江和曲江-石屏两断裂系统的构造动力学与强震序列的关联性. 中国科学:地球科学, 41(5): 713–724.
向宏发, 韩竹军, 虢顺民, 等. 2004. 红河断裂带大型右旋走滑运动定量研究的若干问题. 地球科学进展, 19(S1): 56–59.
徐晶, 邵志刚, 马宏生, 等. 2013. 鲜水河断裂带库仑应力演化与强震间关系. 地球物理学报, 56(4): 1146–1158. DOI:10.6038/cjg20130410
燕乃玲, 李辉, 申重阳. 2003. 丽江地震前后重力场变化的有限矩形位错模型分析. 地震学报, 25(2): 172–181.
杨晓平, 陈立春, 马文涛, 等. 2008. 2007年6月3日宁洱6.4级地震地表变形的构造分析和解释. 地震学报, 30(2): 165–175.
俞维贤, 汪一鹏, 宋方敏, 等. 1996. 1833年云南嵩明8级大地震地表破裂带的考查研究. 地震研究, 19(4): 385–390.
张广伟. 2016. 云南地区地震的重新定位及b值研究. 中国地震, 32(1): 54–62.
张建国, 谢英情, 金明培, 等. 2009a. 中越红河断裂活动性研究. 昆明: 云南科技出版社.
张建国, 刘丽芳, 李西, 等. 2009b. 姚安、大姚中强震区地震构造初析. 地震地质, 31(3): 536–543.
张四昌, 刘百篪. 1978. 1970年通海地震的地震地质特征. 地质科学, 13(4): 323–335.
中国地震局震害防御司. 1999. 中国近代地震目录(公元1912~1990 MS ≥ 4. 北京: 中国科学技术出版社.
周光全, 苏有锦, 王绍晋. 2003. 2001年10月27日永胜6.0级地震的成因探讨. 中国地震, 19(2): 166–174.
周瑞琦, 俞维贤, 谷一山, 等. 1990. 云南耿马7.2级地震地表破裂带研究. 地震地质, 12(4): 291–301.
朱成男, 周瑞琦. 1980. 1979年云南普洱6. 8级地震地表破坏的基本分析. 西北地震学报, 2(2): 64–71, 97.
朱思林, 甘家思, 徐菊生, 等. 1994. 滇西试验场区三维重力反演研究. 地壳形变与地震, 14(1): 1–10.