1971年以来土耳其强震时空特征研究

马 兰 曹广超 曹生奎

（中国西宁810008青海师范大学生命与地理科学院，青海省自然地理与环境过程重点实验室）

引言

地震预测是一个公认的世界性难题，是地球科学研究追求的终极目标之一（张晓东等，2008）.20世纪60年代以来，中期与长期地震预测研究取得了一些进展，如板块边界大“地震空区”的确认、“应力影区”、地震活动性图像、图像识别等（陈运泰，2007）.但由于地震具有突发性和瞬时性的特点，地震预测又受到地球内部的“不可入性”、大地震的“非频发性”和地震物理过程复杂性等困难的制约（陈运泰，2007），世界地震物理预报的准确率仅有10%—20%左右（延军平，延娜，2008），在没有弄清地震物理机理的前提下，可以通过探求地震现象、地震发生时间、地震发生地点的规律，利用地震数学分析预测方法进行地震发生趋势的预测.地震预测方法分为3种类型，即地震统计法、地震地质法和地震前兆法（延军平，延娜，2008）.可公度法是自然特性有序性中隐性规律的特殊表现.利用可公度法来解释地震发生的规律，是对其内在机理的外在体现（延军平，延娜，2008）.翁文波院士等首先将可公度性预测方法从天文学引入信息预测学（翁文波，1981，1984；董国胜，2001；白晶等，2010；延军平等，2010，2011；延军平，闫娜，2008；徐道一，2001；徐道一等，2007），并利用可公度信息系对1976年唐山地震、1982年华北干旱和1991年长江流域洪水等自然灾害现象作过成功的预测（龙小霞等，2006）.

上一次大地震引起的地壳变动，会引起下一次的地震从而形成连锁反应（董颂声，2001）.美国地质学家纳非·托克瑟兹认为，1999年8月17日土耳其大地震发生的地方，恰巧就是上一次（1967年）地震引起的断层裂缝终止的地方（谷利源，1999）.因此，受连锁反应的影响土耳其地震的发生趋势在时间和地理位置上都有一定的内在规律性和较强的预测性.基于此，本文采用“可公度法”和“蝴蝶结构法”对土耳其地区近42a来强震发生趋势进行了判断，以此为地震趋势预测方法的研究以及今后灾害应对提供一些参考.

1 数据与方法

1.1 研究区概况

土耳其地跨亚、欧两洲，位于42°15′—35°40′N，25°36′—44°53′E之间.由于土耳其两侧的阿拉伯板块和欧亚板块的缓慢移动碰撞，北安纳托利亚断层一直很活跃（Monastersky，1999）.这就使土耳其成了一个地震多发国家.

1.2 数据来源

本文所用数据主要来源于国际地震中心（http：∥www.isc.ac.uk／）和中国地震台网中心（2011）关于1971—2011年土耳其地区的地震参数，同时也参考了陈化然等（2000）的研究结果.对获得数据的准确性和完整性做了进一步验证，确保了可公度趋势判断的准确性.

1.3 研究方法

可公度性（commensurability）源于天文学，表示自然界事物之间的一种秩序，是自然现象周期性的一种客观外在反映，是一种信息系（翁文波，1984）.这种信息系不仅在天文学中有明显的表现，而且在其它的学科领域中也广泛存在，因此可以把事物之间的这种可公度性作为对自然现象进行预测研究的依据和方法.

可公度的计算方法包括三元可公度法、四元可公度法和五元可公度法，常用的公式主要有3个（翁文波，1981）：

式中，A，B，C，D，E为以前的重要历史数据，N为预测的未来时间.

蝴蝶对称结构图主要用来解决三元、四元、五元可公度计算中，两个年份的信号都比较强时该如何选择的问题（李双双，延军平，2012）.蝴蝶结构的随机性概率为T＝M／N，其中T为灾害事件时间序列中预测年份发生的可能概率，N为总灾害事件次数，M为参与实际预测的灾害次数，即与主周期关联的年份统计（李双双，延军平，2012）.

2 结果与分析

2.1 土耳其地区强震发生特点

表1为土耳其地区1971—2012年所发生的MS≥6.4地震的记录.由表1可知土耳其地区发生MS≥6.4地震的年份分别为1971，1975，1983，1992，1995，1998，1999，2002，2003年和2011年.

表1 土耳其地区MS≥6.4地震参数Table 1 The parameters of MS≥6.4earthquakes in Turkey

2.1.1 可公度性判断

根据表1的数据，对土耳其地区地震的时间序列分别进行三元可公度法、四元可公度法和五元可公度法的判断.

1）按照三元可公度的计算方法，设x1＝1971，x2＝1975，x3＝1983，x4＝1992，x5＝1995，x6＝1998，x7＝1999，x8＝2002，x9＝2003，x10＝2011，根据公式（1）计算的结果见附表.结果表明，2012年2组，2014年3组，2015年4组，2017年1组，2019年2组，2020年2组，2022年2组，2023年2组.其中2015年的信号最强.

2）按照四元可公度的计算方法，根据公式（2）计算的结果列于附表中.其中，x1＋x3－x2－x2＝4的“4”是等间隔（D），把D加2011得出未来地震的可能发生年份，其它各等式的结果（D）也如此.结果表明，2012年5组，2013年7组，2014年6组，2015年3组，2016年10组，2017年5组，2018年3组，2019年9组，2020年2组，2021年3组，2022年2组，2024年5组，2025年4组，2026年1组，2027年3组.其中2016年的信号最强.

3）按照五元可公度的计算方法，根据公式（3）部分计算结果列于附表中.结果表明，2015年72组，2023年69组.与三元可公度法计算的结果相同，2015年的信号最强.

2.1.2 蝴蝶对称结构图趋势判断

根据土耳其地区发生MS≥6.4地震的时间序列和可公度信息绘制的蝴蝶对称结构图，如图1所示.结果显示，1971年与1983年相差12a，1983年与1995年相差12a，1995年与1999年相差4a，1999年与2003年相差4a.按照三元、五元可公度法的计算结果，如果在2015年发生地震，则按照上述规律，2003年与2015年相差12a，各地震年份之间的差值呈现12，12，4，4a依次交替的规律；而如果按照四元可公度预测的2016年计算，2003年与2012年相差13a，不符合12，12，4，4a依次交替的规律.所以按照这一规律，根据时间序列蝴蝶结构图（李双双，延军平，2012）可以判断土耳其地区未来发生强烈地震的年份为2015年.

图1 土耳其地区MS≥6.4强震时间序列蝴蝶结构图Fig.1 The butterfly structure diagram of MS≥6.4earthquake time series in Turkey

从图2和图3中可看出，与2015年相关的年份有6组，分别是4，12，13，16，20，40a；与2016年相关的周期有2组，分别是13和24a，因此可以得到N＝10，M2015＝10，M2016＝8，则2015年、2016年土耳其地区MS≥6.4强震的蝴蝶结构随机性概率分别为100%和80%.所以可以进一步证明2015年土耳其地区发生强震的可能性最大.

图2 土耳其地区MS≥6.4强震的2015年预测时间序列蝴蝶结构图Fig.2 The butterfly structure diagram of MS≥6.4earthquake time series in Turkey in 2015

图3 土耳其地区MS≥6.4强震的2016年预测时间序列蝴蝶结构图Fig.3 The butterfly structure diagram of MS≥6.4earthquake time series in Turkey in 2016

2.2 土耳其强震发生的空间分布特征及震中迁移规律

图4 中的黑色圆点为土耳其1972—2011年地震的震中（经纬度）位置和依次顺序.圆点1表示1971年地震的震中，2表示1975年地震的震中，3表示1983年地震的震中，4表示1993年地震的震中，5表示1995年地震的震中，6表示1998年地震的震中，7表示1999年地震的震中，8表示2002年地震的震中，9表示2003年地震的震中，10表示2011年地震的震中.其中相同年份发生2次地震的，选取其中震级较大的经纬度.纬度以38.8°N为坐标轴，经度以37.4°E为坐标轴，整个地震点的分布以斜线划分为两个区域，第一次至第四次地震发生在坐标轴的东侧，第五次至第八次地震迁移至坐标轴的西侧，而第九次和第十次地震再次迁移到坐标轴的东侧位置，以每4次地震为一个周期.按照这个规律，则可以预测土耳其地区下一次的地震也发生在坐标轴的东侧位置.

图4 土耳其地震地区空间分布对称图Fig.4 Spatially symmetrical distribution of earthquake areas in Turkey

如果土耳其地区强震震中迁移的物理机制是当一次地震发生后或其断层传播到另一个孕震区时，使其发震或是已发生地震引起的应力场变化影响另一个孕震区发生（郭增建，秦保燕，1983）.那么，下一次强震发生的区域将有两种可能：第一种可能在2011年（第十次地震）附近小幅度迁移；第二种可能是向东北方向迁移.

地震震中迁移是早已为人们所注意的重要现象，国内外曾有多次用震中迁移现象预报大地震的实例（郭增建，秦保燕，1983；李双双，延军平，2012；张国民等，2001）.例如图4中，从第一次地震到第四次地震的迁移路线可以看出，下一次地震的震中以前一次地震的震中为坐标呈现北—南—北—南的回旋迁移规律，即第二次地震的震中位于第一次地震震中的下（偏南）方，第三次地震的震中回旋到第二次地震的上（偏北）方.而第四次地震的震中再次位于第三次地震的下（偏南）方，同理第五次到第八次地震的迁移路线，也呈现北—南—北—南的回旋迁移规律，即第十次地震的震中位于第九次地震的下（偏南）方.所以按照这个规律，可以预测土耳其下一次的强震应发生在土耳其第十次地震，即2011年地震震中的北部，坐标轴的东侧位置，即土耳其的东北部地区.

3 讨论与结论

土耳其位于喜马拉雅——地中海地震带上，地震活动频繁，在地震趋势预测研究上具有很强的代表性.1971年以来，土耳其地区MS≥6.4强震具有明显的空间迁移规律以及良好的可公度性.土耳其强震时空特征的分析可以为地震趋势预测方法的研究以及今后应对重大地震灾害，减少经济、社会等诸方面的损失提供一些参考依据.而对于土耳其地区震中迁移规律和可公度性之间的联系可以作为未来进一步研究的切入点，进行进一步的探索.

本文将“可公度法”和“蝴蝶对称结构法”相结合，可以提高地震预测的准确性和未来趋势判断的可靠性.通过选取土耳其作为强震时空对称特征的研究区域，利用“可公度法”对土耳其地区1971年以来MS≥6.4强震做出趋势判断，预测未来可能再次发生强震的时间，并利用“蝴蝶对称结构法”对1971年以来MS≥6.4强震的震中空间分布及迁移规律进行分析，判断下一次强震发生的区域范围，得到主要结论如下：

1）1971年以来土耳其地区MS≥6.4强震多为浅源地震，大多发生在40°N左右的地区，自1998年以来地震发生较为频繁.

2）土耳其地区强震发生时间呈现较好的“可公度”特征，2015年和2016年强震信号比较强，有可能发生强烈地震，其中2015年信号要强于2016年.

3）强震发生的时间呈现一定的规律性，其中纬度以38.8°N为对称轴，经度以37.4°E为坐标轴.震中主要分布在西南、东北方向，呈北西—南东向对称.

4）震中呈现南北回旋迁移的特征，下一次地震的震中以前一次地震的震中为坐标，呈现北—南—北—南的回旋迁移规律.未来地震的震中可能发生在土耳其第十次地震（2011年地震）震中的北部，坐标轴的东侧位置，即土耳其的东北部，38.8°N以北，37.4°E以东的区域.

白晶，延军平，苏坤慧.2010.新疆MS≥7地震时空对称性及未来趋势判断［J］.内陆地震，24（4）：324－－329.

Bai J，Yan J P，Su K H.2010.The time-space symmetry and the tendency judgement of future MS≥7earthquake in Xinjiang［J］.Inland Earthquake，24（4）：324－－329（in Chinese）.

陈运泰.2007.地震预测：进展、困难与前景［J］.地震地磁观测与研究，28（2）：62－－65.

Chen Y T.2007.Earthquake prediction：Progress，difficulties and prospect［J］.Seismological and Geomagnetic Observation and Research，28（2）：62－－65（in Chinese）.

陈化然，刘文兵，张国民.2000.土耳其地震研究与监测情况［J］.国际地震动态，（1）：35－－38.

Chen H R，Liu W B，Zhang G M.2000.The seismological research and observation in Turkey［J］.Recent Developments in World Seismology，（1）：35－－38（in Chinese）.

董国胜.2001.地震预测与防灾减灾［M］.北京：地震出版社：95－－98.

Dong G S.2001.Earthquake Prediction and Disaster Prevention and Reduction［M］.Beijing：Seismological Press：95－－98（in Chinese）.

董颂声.2001.1999年土耳其伊兹米特地震［J］.地震学刊，21（1）：62－－65.

Dong S S.2001.The 1999Izmit，Turkey earthquake［J］.Journal of Seismology，21（1）：62－－65（in Chinese）.

谷利源.1999.美科学家谈土耳其大地震：三个板块挤压是导因［N／OL］.人民日报.［1999-08-22］.http：∥www.people.com.cn／rmrb／199908／22／newfiles／wzb＿19990822001033＿4.html.

Gu L Y.1999.US scientists talk about the Turkey earthquake，the extrusion from three tectonic plates is the main reason［N／OL］.People’s Daily.［1999-08-22］.http：∥www.people.com.cn／rmrb／199908／22／newfiles／wzb＿19990822001033＿4.html（in Chinese）.

郭增建，秦保燕.1983.地震迁移问题综述［J］.国际地震动态，（1）：8－－10.

Guo Z J，Qin B Y.1983.Review of earthquake migration［J］.Recent Developments in World Seismology，（1）：8－－10（in Chinese）.

李双双，延军平.2012.西北太平洋俯冲带北部地区MS≥8强震时空对称特征［J］.地球物理学进展，27（3）：960－－965.

Li S S，Yan J P.2012.Space-time symmetry of MS≥8earthquake in north-western Pacific plate subduction［J］.Progress in Geophys，27（3）：960－－965（in Chinese）.

龙小霞，延军平，孙虎，王祖正.2006.基于可公度方法的川滇地区地震趋势研究［J］.灾害学，21（3）：81－－84.

Long X X，Yan J P，Sun H，Wan Z Z.2006.Study on earthquake tendency in Sichuan－－Yunnan region based on commensurability［J］.Journal of Catastrophology，21（3）：81－－84（in Chinese）.

翁文波.1981.可公度性［J］.地球物理学报，24（2）：151－－154.

Weng W B.1981.Commensurability［J］.Acta Geophysics Sinica，24（2）：151－－154（in Chinese）.

翁文波.1984.预测论基础［M］.北京：石油工业出版社：47－－49.

Weng W B.1984.Prediction Theory Basis［M］.Beijing：Petroleum Industry Press：47－－49（in Chinese）.

徐道一.2001.大地震发生网络性质：兼论有关地震预测的争论［J］.地学前缘，8（2）：211－－216.

Xu D Y.2001.The network features of large earthquake occurrence and some words on the debate of earthquake prediction［J］.Earth Science Frontiers，8（2）：211－－216（in Chinese）.

徐道一，王明太，耿庆国，汪纬林.2007.翁文波院士的信息预测理论体系的创新性及其意义［J］.地球物理学进展，22（4）：1375－－1379.

Xu D Y，Wang M T，Geng Q G，Wang W L.2007.The creativity of informative forecasting theory and it’s significance［J］.Progress in Geophys，22（4）：1375－－1379（in Chinese）.

延军平，闫娜.2008.关于地震预测体系构建可能性及验证［J］.陕西师范大学学报：哲学社会科学版，37（5）：19－－23.

Yan J P，Yan N.2008.Probability and verification of constructing an earthquake predicting system［J］.Journal of Shaanxi Normal University：Philosophy and Social Sciences Edition，37（5）：19－－23（in Chinese）.

延军平，闫军辉，白晶，王娟，王新华.2010.基于可公度方法的陕西及邻近地区重大自然灾害发生趋势探讨［J］.灾害学，25（2）：18－－20.

Yan J P，Yan J H，Bai J，Wang J，Wang X H.2010.Discussion on trends of major natural disasters in Shaanxi Province and its vicinities based on commensurable method［J］.Journal of Catastrophology，25（2）：18－－20（in Chinese）.

延军平，白晶，苏坤慧，王娟，王新华.2011.对称性与部分重大自然灾害趋势研究［J］.地理研究，30（7）：1159－－1167.

Yan J P，Bai J，Su K H，Wang J，Wang X H.2011.Research on symmetry and tendency of several major natural disasters［J］.Geographical Research，30（7）：1159－－1167（in Chinese）.

张国民，傅征祥，桂燮泰.2001.地震预报引论［M］.北京：科学出版社：127－－129.

Zhang G M，Fu Z X，Gui X T.2001.A General Introduction to Earthquake Prediction［M］.Beijing：Science Press：127－－129（in Chinese）.

张晓东，蒋海昆，黎明晓.2008.地震预测与预警探讨［J］.中国地震，24（1）：62－－65.

Zhang X D，Jiang H K，Li M X.2008.Discussion on earthquake forecast and warning［J］.Earthquake Research in China，24（1）：62－－65（in Chinese）.

中国地震台网中心.2011.中国地震台网（CSN）地震目录［EB／OL］.［2011-11-23］.http：∥www.csndmc.ac.cn／newweb／data／csn＿catalog＿p002.jsp？page＝1＆date＿min＝1971／10／01＆date＿max＝2014／10／31＆longitude＿min＝25.00＆longitude＿max＝45.00＆latitude＿min＝35.00＆latitude＿max＝43.00＆depth＿min＝0＆depth＿max＝1000＆magnitude＿min＝3.0＆magnitude＿max＝10.0＆events＿sortedby＝descent＆region＝no＆region＿en＝China＆region＿cn＝%CD%C1%B6%FA%C6%E4＆region＿mode＝region＿cn.

China Earthquake Networks Center.2011.China Seismic Network（CSN）earthquake catalogue［EB／OL］.［2011-11-23］.http：∥www.csndmc.ac.cn／newweb／data／csn＿catalog＿p002.jsp？page＝1＆date＿min＝1971／10／01＆date＿max＝2014／10／31＆longitude＿min＝25.00＆longitude＿max＝45.00＆latitude＿min＝35.00＆latitude＿max＝43.00＆depth＿min＝0＆depth＿max＝1000＆magnitude＿min＝3.0＆magnitude＿max＝10.0＆events＿sortedby＝descent＆region＝no＆region＿en＝China＆region＿cn＝%CD%C1%B6%FA%C6%E4＆region＿mode＝region＿cn（in Chinese）.

Monastersky R.1999.Turkish earthquake：A wobbly domino falls［J］.Science News，156（9）：132.

附表 利用可公度法计算的土耳其地区未来地震的结果Appendix Future earthquakes in Turkey calculated by commensurability method