全球地震柱的地震层析成像证据^*

 

陈立军, 胡奉湘, 陈晓逢

 

自序：本文于2012年4月投稿国内某刊物，未被录用，故发到博客，以求一乐。

 

摘要：本文依据地震地热说原理，搜集国内外地震层析成像研究的成果，揭示了地震柱的物理属性，得出全球24个地震柱中已有资料的20个地震柱的深部地幔圆锥体均具有P波高速异常体特征，相反，有火山和浅源地震活动而并无中、深源地震活动的地区，比如冰岛和夏威夷，其地幔深部则为P波低速异常体。为解释这个重要的物理现象和构造现象，本文引入空化理论并提出了3个假设。地震柱概念是地震地热说赖以生存的物质基础和构造基础，是与当今流行的地质学派板块构造、地质力学和地幔柱理论的最大区别。深入解读地震柱的物理属性及其构造学意义，或许能为构造地质学和地球动力学带来革命性的大讨论。

关键词：地震地热说，地震柱，地震层析，P波高速异常体，空化理论

 

[中图分类号]P315

 

Imaging Proof of Seismic Tomography in Seismic Cylinders on the Earth

 

Chen Lijun, Hu Fengxiang, Chen Xiaofeng

Earthquake Administration of Hunan Province, 410004 Changsha, China

 

Abstract: According to the principle of the Seismo-geothermics in this paper It was collected some of Chinese and foreign study results of seismic tomography, revealed the physical properties of seismic cylinder and obtained that the deep mantle cone of the existing data of 20 seismic cylinder among 24 seismic cylinder on the Earth is with characteristics of fast velocity anomaly body of P wave, on the contrary, in the regions of volcanoes and shallow seismic activity but deep seismic activity, such as Iceland and Hawaii, the deep mantle cone is lower velocity anomaly body of P-wave. For the interpretation of this important physical phenomena and structural phenomena, the theory of cavitation was introduced in this paper and put forward three hypotheses. It was said that the conception of seismic cylinder was the material basis and tectonic basis for survival of Seismo-geothermics, and was maximum difference with today's popular geological theory such as plate tectonics, Geomechanics and the mantle plume theory. Understanding the physical attributes and its tectonic significance of seismic cylinder, it will bring perhaps revolutionary element for the tectonics and dynamics of the Earth.

Key words: Seismo-geothermics, seismic cylinder, seismic tomography, fast body of P wave, cavitation

CLC Numbers:P315

* 中国地震局2012年度老专家科研基金课题（编号201229）和湖南省地震局2012年度科研课题（编号201201）资助项目

 

1 前言

 

作者在长期学习和研究大陆漂移、板块构造、地质力学、地幔柱理论的基础上，创立了地震地热说^[1-2]。地震地热说认为，全球有两大构造系统，即由24个地震柱组成的M型热机带和由大洋中脊组成的W型冷机带（图1）。热系统主宰全球几乎100%7级以上的壳内强震和近90%活火山的灾变，冷系统则统领全球壳内应力场的调整和缓慢的大陆漂移运动。

 

 

 

图1 全球地震柱分布图

 

地震柱概念区别于传统的地幔柱概念（mantle plume）。地幔柱理论诞生于1971年，Morgan 提出全球约有20个“热点”，并认为热点源区物质对流上升，其结构类似于“直达深部地幔的管子”，植根于地幔深部，且位置相对固定从而解释地表“热点”位置相对固定的特点^[3]。该理论因为无法验证而遭到诸多质疑^[4-5]。

本文的地震柱概念（seismic cylinder），是指由中、深源地震活动、壳内地震和活火山所构成的地幔深部或上地幔内中、深源地震震源密集、顶点向下的圆锥形柱状地质体。因此，地震柱的划分是以地震和火山活动的共同特征来确定的,是地幔深部能量与地表能量交换的稳定通道，具有可视性、唯一性、独立性和稳定性。地震柱形状类似于树，具有树根、树干和树冠，姑且称之为树型构造。同一个地震柱内，可以是单树型、双树型或多树型构造。还有树杈型，即同一个树根和树干上带有多个树杈的。全球共计定义了21个Ⅰ级地震柱，3个Ⅱ级地震柱。Ⅰ级和Ⅱ级地震柱的分支又演变出26个Ⅲ级地震柱。地震柱的编号与命名结果见表1。

研究结果表明，地震柱内部的中、深源地震以热能形式自下而上逐层驱动，到达壳内并推动地壳构造活动，引发壳内强震与火山喷发，最后通过近地表层晚期随机地震活动将剩余能量消减完毕，构成完整的地震链过程。如此往复，周而复始，24个地震柱活动错落交替，构成全球地震与火山活动复杂的过程与格局^[2]。

地震柱的活动机制又区别于“俯冲”和“碰撞”的概念，是一种热能的活动而不是冷物质的机械运动。因此，本文将引用国内外科学家们所做的地震层析成像技术的成果，就地震柱的物理属性进行初步的讨论。

 

表1  全球地震柱命名（修订版）

 

 

2 地震柱的地震层析成像证据

 

据江国明等（2009）介绍^[6]，自从Aki和Lee（1976）的开创性工作以来^[7]，天然地震层析成像方法已成为人们研究地球内部速度结构的一种非常有效的工具。根据所用的地震数据不同，地震层析成像法可分为体波成像和面波成像；根据震源是否在研究区域内，体波成像又可分为近震层析成像和远震层析成像。由于近震射线在地壳内的交叉较好，所以近震成像法适合研究地壳及上地幔顶部的速度结构^[8]，而远震射线在区域地震台网下的上地幔内交叉较好，因此远震层析成像法在研究上地幔速度异常分布、地幔柱形态以及俯冲带速度结构等方面均具有直接、快捷的优势^[9]。

本文收集国内外专家的工作成果200余篇，将其中与地震柱直接有关的结果罗列如下。

 

2.1 南美洲地震柱的地震层析成像证据

 

据J. Trampert（2005）等在《全球地震层析成像定量解释的建立》一文的结果^[10],可以给出南美洲地震柱的物理结构如图2所示。图中a为该地震柱内部中、深源地震震源用Origin软件绘制的三维立体分布图像，包括南智利地震柱和北智利地震柱，b为地震柱内P波高速异常体的大致图像，是由Trampert的平面图像换算为三维坐标后用Origin软件绘制的立体图像（下同）。由图可见，a、b两图的空间图像是大致相似的，说明本地震柱应该是一个P波高速异常体。

 

 

a)                                     b)

图2 南美洲两地震柱及其地震层析结果立体图像示意图

（a.地震目录ANSS，1963-2011；b.层析结果立体图像据J.Trampert等改编）

 

2.2 危地马拉地震柱的地震层析成像证据

 

据S. Husen（2003）等在《由本地地震层析和岩石学模型确立的哥斯达黎加俯冲带结构和岩浆过程》一文的结果^[11],可以给出危地马拉地震柱的物理结构如图3所示。由图可见，本地震柱应该是一个P波高速异常体。图3b地震层析结果的研究区范围偏窄，未能覆盖整个地震柱，因而造成左右两图的位置偏离。

 

 

a)                                     b)

图3 危地马拉地震柱及其地震层析结果立体图像示意图

（a.地震目录ANSS，1963-2011；b.层析结果立体图像据S.Husen等改编）

 

2.3 鄂霍次克海地震柱的地震层析成像证据

 

据Guoming Jiang等（2009）在《太平洋板块边缘勘察加的地震层析》^[12]和江国明等（2009）在《远震层析成像中的地壳校正研究及应用》^[6]的结果,可以给出鄂霍次克海地震柱的物理结构如图4所示。图4a中不包含白令海地震。由图可见，本地震柱也是一个P波高速异常体。但是，图中的a和b存在较大的偏差，可能出于如下原因。

 

 

a)                                     b)

图4 鄂霍茨克地震柱及其地震层析结果立体图像示意图

（a.地震目录ANSS，1963-2011；b.层析结果立体图像据Guoming Jiang等和江国明等改编）

 

（1）地震层析的数学模型问题，或者地震层析的空间范围两者不完全匹配。

（2）或许，本地震柱的中、深源地震主要活动在高速异常体旁的低速层内。

目前来看，第一种可能性要大一些。远震层析中上地幔内的交叉尺度较大，肯定存在一定的系统误差，因而要确定高速异常体的准确位置是有难度的。

 

2.4 汤加地震柱的地震层析成像证据

 

据J.P. Brun等（2009）《俯冲带动力学》、Dapeng Zhao等（2000）《地震层析成像对地球科学的影响》和Serge Lallemand等（2009）《俯冲带动力学》的结果^[13-15],可以给出汤加地震柱的物理结构如图5所示。图5a包括东汤加地震柱和西汤加地震柱。由图可见，本地震柱同样是P波高速异常体。图5b中向上有两个分枝，分别对应东、西两个地震柱。

 

 

a)                                     b)

图5 汤加地震柱及其地震层析结果立体图像示意图

（a.地震目录ANSS，1963-2011；b.层析结果立体图像据Dapeng Zhao等

和Serge Lallemand等改编）

 

2.5 地中海地震柱的地震层析成像证据

 

 

a)                                     b)

图6 地中海地震柱及其地震层析结果立体图像示意图

（a.地震目录ANSS，1963-2011；b.层析结果立体图像据M. Nafi Toksöz等和W. Lowrie改编）

 

据J.P. Brun等（2009）《俯冲带动力学》^[13]、M. Nafi Toksöz等（2010）《阿拉伯-欧洲碰撞带及邻区的地震层析》^[16]和W. Lowrie（2007）《地球物理学基础》^[17]的结果,可以给出地中海地震柱的物理结构如图6所示。由图可见，本地震柱亦是P波高速异常体，而且具有3个分支，即意大利分支、土耳其分支和伊朗分支。3个分支由西向东渐弱，波速异常与地震活动水平相一致。图6a中不含兴都库什地震。

 

3 非地震柱热点的地震层析成像证据

 

在所检索到的200多篇文献中，未发现不存在地震柱的地方出现明显的P波高速异常现象。相反，在冰岛和夏威夷两地，火山活动十分活跃，但是没有中、深源地震活动，不能定义为地震柱，其波速异常情形也与地震柱的情形迥异。

据G.R.Foulger等（2001）《地震层析成像显示冰岛是上地幔狭窄通道内的上涌》^[18]的结果,可以给出冰岛的热柱结构如图7a所示，据Dapeng Zhao（2009）《多道地震层析与地幔动力学》^[19]的结果,可以给出夏威夷的热柱结构如图7b所示。由图可见，这两个热点都是P波低速异常体。

 

 

a)                                     b)

图7 冰岛和夏威夷的地震层析结果立体图像示意图

（a.冰岛，据G.R.Foulger等改编；b.夏威夷，据Dapeng Zhao改编）

 

4 关于地震层析结果的若干假定

 

上述已经考察过的7个地震柱均为P波高速异常体，其中包括南美洲（含01号南智利和02号北智利）、03号危地马拉、06号鄂霍次克海、汤加（含14号西汤加和15号东汤加）和19号地中海地震柱。其实也还有关于04号海地^[10]、05号白令海^[12]、07号日本^[19]、08号中国珲春^[19-20]、马里亚纳海沟（含09号北马里亚纳和10号南马里亚纳）^[21]、11号台湾及琉球^[19、21]、16号印尼（http://davidpratt.info/subduct.htm）、18号兴都库什^[16]、21号南桑威奇^[10、22]、F3号北美洲（Dr. Stephan Husen，http://www.seg2.ethz.ch/boschil/coursetomo08/）等地震柱的零星资料，同样也具有P波高速异常体的特征。相反，在没有中深源地震而有火山活动的地方，比如冰岛和夏威夷，地震层析结果则表明为P波低速异常体。

是什么原因造成如此泾渭分明的差异？这里可能涉及到深源地震的成因。1970年代以前就有关于深源地震机制的相变说^[23]。但是，此说是基于俯冲带冷物质亚稳态相变的，依然是将地表地震的断裂错动机制拓展到了几百公里深度以下。

然而，深源地震全部集中在地震柱之内。地震柱是一种P波的高速异常体，表现为热能上涌的模式^[2]，不能沿用壳内地震的断层破裂机制，必须另辟蹊径。因此，本文引入了C. E. Brennen（1995）的空化理论^[24]：在一定的热力学条件下，地震柱内部可能产生空化现象，由相变—成核—生长—溃灭而完成深源地震的孕育与破裂机制。

沿此思路，并参照李四光（1973）的地质力学理论^[25-26]，本文提出了关于波速异常的如下若干假设。

地震层析中P波速度的定义为（http://en.wikipedia.org/wiki/P-wave）

式中M为P波的弹性模系数，M = K + 4μ / 3，ρ为波动介质密度。在假设存在上涌侵入体而且假设M不变的情况下，Vp只与密度ρ有关，于是可以提出关于P波速度异常的如下假定：

假定1：在地球自转角速度变更的适当相位，部分地幔深部重物质沿着某种习惯性通道向地表方向甩出，形成上地幔的某种热柱，构成地震、火山、大陆漂移的主要能源来源；在地球自转角速度变更的适当相位，地幔深部上涌的重物质由于向心力和重力作用回落；周而复始，构成热柱的活动韵律。

假定2：部分热柱满足空化的热力学条件，由相变—成核产生空泡或空泡流，构成地震柱。地震柱内空泡逐渐生长并上升到一定深度溃灭，形成中、深源地震活动，释放热能，并由下而上构成逐层驱动之势。由于空泡或空泡流的存在，使得地震柱内物质密度低于柱外环境物质，从而形成P波的高速异常体。

假定3：部分热柱不满足空化条件，没有中、深源地震活动，但由于地幔上涌物质密度略高于环境物质，从而形成P波的低速异常体。

假定2类似热机工作原理，适用于全球24个地震柱。引发6级、7级中、深源地震所需空泡尺寸很大，因此地震越活跃的地震柱内P波高速异常的反差越大。假定3适用于冰岛、夏威夷以及大洋中脊带内有火山和浅源地震而无中、深源地震活动的地区。

 

5 讨论和结论

 

本文根据地震地热说的原理，收集国内外大量地震层析成像成果，研究了全球24个地震柱的物理特性。研究结果表明，地震柱属于P波的高速异常体。高速异常体的形态与地震柱内中、深源地震活动的的立体分布基本吻合，因而与地震柱内的本尼奥夫剖面也基本一致。于是，本尼奥夫剖面和P波的高速异常体构成了地震柱的两个重要的物理属性。

在地震地热说的原理中我们知道，本尼奥夫剖面只能存在于地震柱的内部，横向无法扩展，相邻地震柱的本尼奥夫剖面无法相互沟通^[2]。尽管地震层析并不像本尼奥夫剖面那样确定，而且在岩石圈内由于岩体密度差异较大而使得P波速度的扰动较大，但是在一定深度以下，上地幔和/或地幔深部，P波速度的扰动将趋于稳定，根据第3章和第4章的比较研究或许可以推论，深部的P波高速异常体也只能存在于地震柱的内部，横向无法扩展，相邻地震柱的高速异常体也无法相互沟通。

因此，人们将图1上24个地震柱所处的M型热机带定义为“俯冲带”，可能是一种思维方法的错误。尽管图上看起来24个地震柱比较密集分布，但是这些地震柱只是一些深度达到几百公里的圆锥体，除了太平洋西岸部分地震柱间距较小，大多数地震柱间距都在1000km以上，其P波高速异常体和本尼奥夫剖面又都受到空间尺度的限制，根本无法横向勾连，无法成带。人们发现了本尼奥夫剖面，也发现了P波高速异常体，在深度上的图像很漂亮，于是在宽度上加以盲目扩展，构成大量的板块俯冲带模型^[27]，显然是在无意中夸大或误解了地震柱的这两个物理属性。

加之，地震柱内部的活动机制也并非冷物质的俯冲，但这不在本文的讨论之列。但是，无论怎么说，试图用“俯冲”和“碰撞”的概念来万能地解释M型热机带上占全球90%以上的所有壳内强震和火山喷发的灾变，是容易使得人们思想僵化的。人们只要听说日本9级地震是碰撞和俯冲造成的，于是就不求甚解了。这是人类的一种悲剧。

地震柱是地震地热说赖以生存的构造基础和物质基础，是与当今流行的地质学派^[28-31]板块构造、地质力学和地幔柱理论的最大区别。深入解读地震柱的物理属性及其构造学意义，或许能为构造地质学和地球动力学带来革命性的大讨论。

 

本文引用的地震资料取自http://www.ncedc.org/ 网页的ANSS地震目录，火山资料取自http://www.volcano.si.edu/网页的GVP火山月报，谨此致谢。

 

2012-3-16 初稿

 

参考文献（Reference）

 

[1] 陈立军.中国地震震源深度与强震活动状态研究[J].地震地质.2000,22(4)： 360-370

[2] 陈立军. 地震地热说原理与应用[J].内陆地震, 2012, 26(2)

[3] G.F. Davies. 地幔柱存在的依据[J].科学通报,2005.9,50(17):1801-1813

[4] G.R. Foulger.地幔柱：为什么现在怀疑?科学通报, 2005.9,50(17):1814-1819

[5] 牛耀龄. 关于地幔柱大辩论[J]. 科学通报, 2005.9, 50(17):1797-1800

[6] 江国明，赵大鹏，张贵宾．远震层析成像中的地壳校正研究及应用[J]．地球物理学报，2009，52(6)：1508～1514

[7] Aki K，Lee W H K．Determination of three dimensional velocity anomalies under a seismic array using first Parrival times fromlocal earthquakes，J，A homogeneous initial mode1．J．Geophys．Res．1976，81：4381~4399

[8] 田有，赵大鹏，孙若昧，等．1992年美国加州兰德斯地震一地壳结构不均匀性对地震发生的影响[J]．地球物理学报，2007，50(5)：1488-1496

[9] 赵文津，薛光琦，吴珍汉等．西藏高原上地幔的精细结构与构造—地震层析成像给出的启示[J]．地球物理学报，2004，47(3)：449-455

[10] J. Trampert, Robert D. van der Hilst. Towards a Quantitative Interpretation of Global Seismic Tomography. INTERPRETATION OF TOMOGRAPHY. Earth’s Deep Mantle: Structure, Composition, and Evolution Geophysical Monograph Series 160. Copyright 2005 by the American Geophysical Union 10.1029/160GM05. P47-62

[11] S. Husen, R. Quintero, E. Kissling et al. Subduction-zone structure and magmatic processes beneath Costa Rica constained by local earthquake tomography and petrological modelling. Geophys J. Int.(2003)155,11-32

[12] Guoming Jiang, Dapeng Zhao, Guibin Zhang. Seismic tomography of the Pacific slab edge under Kamchatka. Tectonophysics, 465(2009), 190–203

[13] J.P. Brun, O. Oncken, H. Weissert,et al. Subduction Zone Geodynamics.Roma：Springer-Verlag Berlin Heidelberg，2009，p.50-55

[14] Dapeng Zhao, J. R. Kayal. Impact of seismic tomography on earth sciences. SPECIAL SECTION: SEISMOLOGY 2000. CURRENT SCIENCE, VOL. 79, NO. 9, 10 NOVEMBER 2000. 1208-1214

[15] Serge Lallemand • Francesca Funiciello. Subduction Zone Geodynamics，Springer-Verlag Berlin Heidelberg, 2009

[16] M. Nafi Toksöz, Robert D. Van der Hilst, Margaret H. Benoit, et al. SEISMIC TOMOGRAPHY OF THE ARABIAN-EURASIAN COLLISION ZONE AND SURROUNDING AREAS. 29th Monitoring Reseach Review: Ground-Based Nuclear Explosion Monitoring Technologies. Sponsored by Air Force Research Laboratory, (2010) 292-301

[17] W. Lowrie. Fundamentals of Geophysics, Second edition. Cambridge University Press. Published in the United States ofAmerica by Cambridge University Press,New York, 2007.

[18] G.R. Foulger, M.J.Pritchard, B.R.Julian, et al. Seismic tomography shows that upwelling beneath Iceland is confined to the upper mantle. Geophys.J.Int. (2001) 146, 504–530

[19] Dapeng Zhao. Multiscale seismic tomography and mantle dynamics[J]. Gondwana Research 15 (2009), 297–323

[20] 赵大鹏，雷建设，唐荣余. 中国东北长白山火山的起源：地震层析成像证据[J].科学通报，2004，49(14):1439-1446

[21] 刘建华，刘福田，胥颐，等. 西北太平洋俯冲带的形貌特征及动力学状态探讨[J]. 海洋地质地球物理，2012（待刊）

[22] Daria Zandomeneghi. Passive and Active Seismic Tomography of Volcanic Islands Sao Miguel (Portugal) and Deception (Antarctica). Dissertation for the degree of Doctor of Earth Sciences

[23] 刘雷、杜建国、易丽.亚稳态橄榄石相变与深源地震研究进展[J].地震，2007，27(3):41-49

[24] C. E. Brennen. Cavitation and Bubble Dynamics[M]. Oxford University Press, 1995

[25] 李四光.地质力学概论[M].1973.北京：科学出版社

[26] 李四光.地质力学方法[M].1976.北京：科学出版社

[27] 金性春.板块构造学基础[M].1984.上海：上海科学技术出版社

[28] 赵文津.大陆漂移，板块构造，地质力学.地球学报，2009，30(6):717-731

[29] 胡明、廖太平.构造地质学[M].石油工业出版社，2007

[30] Geoffrey F．Davies. 地幔柱存在的依据. 科学通报·前沿[J]，2005，50(17)：1801-1813

[31] 童航寿. 地幔柱构造研究概述[J]. 铀矿地质. 2009，25（4）：193-201