1974年和1979年溧阳两次地震同震库仑应力变化及其发震构造研究

储飞　王韶稳　张毅　谈昕　袁泉　徐如刚

摘要：

1974年和1979年在苏皖交界的溧阳地区分别发生了两次破坏性地震，造成了大量人员伤亡和财产损失。分别计算这两次地震引起的同震库仑应力变化，并与余震分布情况进行对比分析，认为两次地震的发震断层均为NW走向。计算得出溧阳1979年地震震中位置处于1974年地震引发的静态库仑应力增加的区域，应力变化值约0001 MPa，表明1974年地震的发生对1979年地震的发生具有一定的积极影响，但由于应力变化值较小，未能造成直接的触发。

关键词：溧阳地震；库仑应力变化；震源机制解；苏皖交界

中图分类号：P31572文献标识码：A文章编号：1000-0666（2016）03-0397-07

0引言

1974年4月22日在苏皖交界的溧阳地区发生了MS55地震，震源深度16 km（贺楚儒等，1988），5年后的1979年7月9日几乎在同一位置又发生了MS60地震，震源深度12 km（贺楚儒，张德齐，1990）。5年内原地重复发生两次破坏性地震的情况在中国东部的板内地震中非常少见（叶洪等，1980）。

溧阳两次地震造成了大量的人员伤亡和财产损失。地震发生后专家们先后对其进行了大量研究分析。对1974年MS55地震的发震构造问题，贺楚儒等（1983）、叶洪等（1980）和Chung等（1995）意见较为一致，认为其发震断裂应为NWW向断裂。而对1979年MS60地震的发震断裂问题则一直没有较为统一的意见。叶洪等（1980）应用零矢量方法研究了1979年MS60地震与其余震的震源机制解，分析认为该地震是在NEE向挤压应力作用下，沿着NE向的茅山断裂带东支发生了右旋走滑兼有正断层性质的错动。林邦慧等（1982）、刘万琴和黄家正（1982）分别利用P波频谱和瑞利波方向性函数对溧阳MS60地震进行研究，认为此次地震是NW向断层向SE方向发生的左旋走滑破裂。Chung 等（1995）利用2°～90°之间的长周期P波、SH波及短周期远震PD波的波形数据，确定了此次地震具有逆冲分量的走滑震源机制，但对于此次地震的发震断裂没有给出具体的意见。对于两次地震之间的关系，朱春俊（1992）、胡连英（1997）和叶洪等（1980）认为这不是两次孤立的地震事件，两次地震之间必然存在一定的联系。

本文分别计算了溧阳两次地震引起的同震库仑应力变化，并与余震分布情况进行对比分析，根据分析结果提出1979年溧阳MS60地震发震断裂的意见，判断1974年地震是否对1979年地震具有触发作用。

1库仑应力变化研究概述

库仑破裂应力是指作用在断层面上，克服摩擦阻力使断层滑动的力，而库仑破裂应力变化是指库仑破裂应力的改变量（王连捷等，2008）。地震发生时所释放的断层上所积累的应力并没有完全凭空消失，而是通过传输，转移至周边地区，导致邻近区域内的应力升高或者降低，进而触发或者抑制这些地区的地震活动（单斌等，2010）。近几年来，大地震产生的库仑破裂应力变化对后续的地震活动造成的影响受到广泛关注（Harris，1998；Parsons，Dreger，2000；Toda et al，2011；石耀霖，曹建玲，2010），并已被诸多震例研究所证实（Stein et al，1994；King et al，1994；Toda et al，2008）。许多震例表明，大地震后库仑应力的变化对后续地震活动有促进作用，很小的库仑应力变化就可能对地震的活动性及时空分布产生较大的影响（Stein et al，1997；Deng，Sykes，1997）。

自2008年汶川大地震发生后，地震引发的库仑应力变化研究迅速受到了研究学者们的关注，分别从不同角度对库仑应力计算过程中的若干问题、库仑应力变化与余震分布的对应关系、庫仑应力变化对周边断裂活动的影响等方面进行了深入的探讨。

库仑破裂应力变化的表达式如下（Okada，1992；King et al，1994）：

ΔCFS=Δτ+μ′Δσn（1）

式中，ΔCFS为库仑破裂应力变化，Δτ为剪切应力变化，Δσn为正应力变化，μ′为有效摩擦系数。根据库仑破裂准则，断层面上的剪切应力大于或等于阻碍滑动的摩擦力时，沿断层面就会发生摩擦滑动。诸多震例研究成果表明，地震所引发的库仑应力增加的区域有利于后续地震的发生，而负值区则明显抑制后续地震的发生（Lin，Stein，2004；Toda et al，2011）。

2震区地震地质背景

溧阳地震处于华南块体的扬子地台接近与华北块体接壤的北部边界（Chung et al，1995）。这一地区的典型特点是一系列晚中生代到早新生代和中新生代时期形成的隆起和凹陷相间展布（贺楚儒，张德齐，1990）。溧阳MS60地震即发生在茅山东侧断陷盆地内，该盆地是一个新生代断陷盆地（叶洪等，1980）。隆起和凹陷多被一系列规模较大的断裂分隔开来。

震区与此次地震活动相关的主要断裂有（图1）：

（1）茅山断裂带。茅西断裂位于茅山西侧，总体呈NE30°方向延伸，倾向SE，总长约110 km（江苏省地质矿产局，1984）。茅东断裂为NNE向，切割NE向的逆冲断裂和NNE向推覆构造，长约85 km，断裂走向NE10°～20°，倾向多为SE，主要出露在金坛致和至高淳花山。该断裂具有多期活动特征，属左行兼逆冲和旋转运动特质（胡连英，徐学思，1992）。

（2）金坛—南渡断裂。呈NNE向延伸，长约80 km，1974年和1979年溧阳2次地震的余震多沿该断裂呈带状分布。跨越金坛—南渡断裂的水准观测资料表明，在溧阳1979年MS60地震前后具明显差异（侯康明等，2012）。根据上述资料及地震活动情况，可推断该断裂应属于第四纪晚更新世活动断裂。

（3）溧阳—南渡断裂。该断裂向近EW向延伸长约20 km，与金坛—南渡断裂相交，大部分为第四系所覆盖，属隐伏断裂（侯康明等，2012）。

31979年溧阳MS60地震引起的库仑应力变化及其发震断裂讨论

1979年7月9日溧阳地区发生MS60地震，其后又发生了多次余震。自1979年7月9日主震后至1980年1月31日共记录到余震336次（MS≥0），其中MS≥40的余震有2次，均为走滑型错动（叶洪等，1980）。由图2可以看出，溧阳MS60地震震中位于NE向的金坛—南渡断裂与近EW向的溧阳—南渡相交的位置，余震大多分布在主震震中附近，大致沿EW向溧阳—南渡断裂分布，另有少数余震分布较为分散。整体上来看，几乎所有的余震都处在金坛—南渡断裂的东侧，极少数余震发生在其西侧。

溧阳MS60地震发生后，国内外相关科研机构和学者陆续对该地震进行了深入的研究，并给出了震源机制解（表1）。由表1可以看出，震源机制解显示有两组节面，一组节面走向为NW向，另一组节面走向为NE向，2组节面的滑动方式均以走滑错动为主，兼有少许逆冲分量。NW向节面走向与此区域近EW向的溧阳—南渡断裂较为接近，且溧阳MS60地震发生之后的余震活动主要沿此断裂分布，而NE向节面走向则与NE走向的金坛—南渡断裂较为吻合。

表1中所列2组不同机构提供的震源机制解数据虽然有所不同，但差别极小，因此本文选取第二组震源机制解参数（Chung et al，1995）作为发震断层，利用Coulomb33程序进行同震库仑应力变化计算，根据中国地震台网数据设定震源深度为12 km，断层面参数根据Coulomb33程序所提供的经验公式计算获得，计算深度设定为12 km。

由计算出的应力变化与余震分布对比图样（图3）可以看出：（1）震中附近区域主要为应力减小区域，在距震中约数千米的外围区域，应力变化呈现不同分布，应力增加区域和应力减小区域呈间隔分布；（2）对比计算所得的同震库仑应力变化与余震分布可得，由NE走向的节面II所得的计算结果中只有少数余震处在应力增加的区域，其余则处在应力减小区域和边缘地带；而由NW走向的节面I所得的计算结果中绝大多数余震都分布在应力增加的区域，且余震发生、发展的趋势与应力增加区域的方向趋势和范围非常一致。因此，相比于NE走向的节面II，NW走向的节面I计算所得的应力变化增加的区域与余震的分布更加吻合。

通过将地震造成的静态库仑应力变化与余震分布进行对比，比较结果更支持1979年溧阳地震的发震断层为NW走向，这与林邦慧等（1982）以及刘万琴和黄家正（1982）的研究结果较一致。

41974年溧阳MS55地震引起的库仑应力变化及其对1979年地震的触发作用分析

41库仑应力变化

1974年4月22日，在江苏省溧阳西部地区发生了MS55地震。由于震中区人口集中，河网密布，尽管只是一次中强地震，仍然造成了较大的损失。此后的1979年在距离非常近的同一地区又发生了MS60地震，这种地震重复性发生的现象在中国东部地区是很少见的（Chung et al，1995），因此，时至今日依然受到地震研究学者们的高度关注。图4给出了1974年溧阳MS55地震及其余震震中的空间分布。谢瑞征等（1980）和高祥林等（1993）利用P波初动对此次地震进行了研究，认为本次地震是带有正滑分量的走滑错动，基于对余震和最大烈度等震线分布的研究，认为发震断层面应该为NWW向。Chung 等（1995）通过对地震波波形研究的结果表明本次地震活动为NWW走向断层具有逆冲分量的走滑错动，并给出了震源机制解数据。

笔者利用Coulomb33程序对1974年溧阳地震所引发的同震库仑应力变化进行了计算，根据中国地震台网（CSN）数据设定震源深度为16 km，断层面参数根据Coulomb33程序所提供的经驗公式计算获得，计算深度设定为18 km。

图5可知：（1）在震中附近应力变化相对较为复杂，距震中稍远处应力变化趋势较为明显；（2）震后余震在震中附近分布不多，主要集中分布在震中SE方向，这与计算出的静态库仑应力增加区域的分布比较一致；（3）与NE向节面相比，将NW向节面设为发震断层计算出的结果与余震分布更为吻合。

421974年溧阳地震对1979年溧阳地震的触发作用分析

大量研究通过计算地震引起的库仑应力变化发现地震之间存在相互影响、相互作用的联系（King et al，1994；Freed，2005；Stein ，Lisowski，1983；Toda et al，1998；万永革等，2007；刘桂萍，傅征祥，2000）。地震的发生释放了聚集的应力，但应力并不会完全凭空消失，部分应力通过传输、转移至其他地区，进而影响后续地震的发生（王连捷等，2008）。库仑应力的增强相当于断层额外负荷的加载，促进地震的发生；反之，将导致地震的发生时间推迟。

以Chung等（1995）提供的1974年溧阳地震震源机制解NW向节面数据作为源断层，1979年溧阳地震发震断层面作为接收断层，计算了1974年溧阳地震引起的库仑应力变化对1979年地震的影响，计算深度设为12 km，计算出的应力变化如图6所示。结果显示，1979年溧阳地震震中位置处于应力增加的区域，但应力变化值在0001 MPa左右，小于001 MPa（01 bar），而001 MPa通常被认为是能够有效促进后续大地震发生的阀值（Harris，1998）。因此，笔者认为溧阳1974年地震引起的静态库仑应力变化对1979年地震的发生具有一定的积极影响，但未能造成直接的触发。

5讨论与结论

溧阳1974年和1979年地震是发生于中国东部大陆稳定地区的2次破坏性地震，因其破坏性及短时期内地震原地重复发生而备受关注。本文对溧阳1974年和1979年2次地震序列和震区构造背景进行了简要分析，分别计算了2次地震引起的静态库仑应力变化，并与其余震分布情况进行了对比分析，得到以下结论和讨论：

（1）2次地震的计算结果均显示将NW向震源机制解节面数据设为发震断层时的计算结果与余震分布的吻合度要明显优于NE向节面的计算结果，从余震分布的趋势与应力增加区域的发展方向和范围比较也呈现同样的结果，因此本文认为溧阳两次地震的发震断层均应为NW走向断层。

（2）以1974年地震的发震断层作为源断层，以1979年地震的发震断层作为接收断层，计算了1974年地震引发的静态库仑应力变化对1979年地震的影响，计算结果显示1979年溧阳地震震中位置处于应力增加的区域，但应力变化值相对较小，未能达到应力触发的阀值。因此，笔者认为溧阳1974年地震的发生对1979年地震的发生具有一定的积极影响，但由于应力变化值较小，未能造成直接的触发。

（3）两次地震的发生对周围断裂的活动性产生了一定的影响，尤其是对距离较近的溧阳—南渡断裂、茅山断裂带和金坛—南渡断裂影响较大，后续跨断层水准观测资料（侯康明等，2012）也证实这一点。

（4）地震的发生是一个十分复杂的过程，难以用一种简单的模型和方法来解释，库仑应力的变化只是地震的一个触发因素，能否发生地震还取决于斷层所积累的构造应力状态、大小及其变化情况。虽然如此，库仑应力的变化对研究余震发展趋势和周边区域地震活动仍具有重要意义。

参考文献：

高祥林，胡连英，徐学思等.1993.茅山东侧断裂的断层运动及其与溧阳地震关系的动力学模型[J].地震研究，16（4）：401-409.

贺楚儒，彭纯一，张四维等.1983.溧阳六级地震的磁异常与断层面的破裂特征[J].地震地质，5（2）：52-58.

贺楚儒，张德齐，张绍治.1988.1974年4月22日江苏省溧阳5.5级地震[M]//张肇诚，罗兰格，李海华等.中国震例（1966-1975）.北京：地震出版社.

贺楚儒，张德齐.1990.1979年7月9日江苏溧阳6.0级地震[M]∥张肇诚，罗兰格，李海华等.中国震例（1976-1980）.北京：地震出版社.

侯康明，熊振，李丽梅.2012.对江苏省溧阳2次破坏性地震发震构造的新认识[J].地震地质，34（2）：303-312.

胡连英，徐学思.1992.茅东断裂带运动与溧阳地震成因机理研究新进展[J].国际地震动态，（6）：5-8.

胡连英.1997.江苏溧阳地震震源构造[J].江苏地质，21（3）：169-176.

江苏省地质矿产局.1984.江苏省及上海市区域地质志[M].北京：地质出版社.

林邦慧，魏富胜，刘万琴等.1982.溧阳—介休—五原北西地震带强震的破裂特征[J].地震学报，4（2）：116-126.

刘桂萍，傅征祥.2000.1976年7月28日唐山7.8级地震触发的区域地震活动和静应力场变化[J].地震学报，22（1）：17-26.

刘万琴，黄家正.1982.用瑞雷波广义方向性函数研究5-6级地震的破裂过程[J].地震学报，4（1）：27-34.

单斌，熊熊，郑勇.2010.2010年玉树地震同震库仑应力变化及地震活动性预测研究[C]∥中国地球物理学会年会、中国地震学会第十三次学术大会.

石耀霖，曹建玲.2010.库仑应力计算及应力过程中若干问题的讨论——以汶川地震为例[J].地球物理学报，53（1）：102-110.

万永革，沈正康，曾跃华等.2007.青藏高原东北部的库仑应力积累演化对大地震发生的影响[J].地震学报，29（2）：115-129.

王连捷，周春景，孙东生等.2008.汶川5.12地震引起的库仑应力变化及其对周边地震活动性的影响[J].地质力学学报，14（3）：193-200.

谢瑞征，陈晓明，李端路.1980.江苏溧阳6级地震的发震构造[J].地震学刊，3（4）：27-32.

叶洪，张文郁，于之冰等.1980.1979年溧阳6级地震震源构造的研究[J].地震地质，2（4）：27-38.

朱春俊.1992.溧阳两次地震的构造背景和发震构造[M]∥徐玉明.溧阳地震工作文集.北京：地震出版社.