震源机制解综述.doc

震源机制解综述1、引言地震学是一门以观测资料为基础的研究地震的成因及其规律已成为地震预报的一种重要手段，它的发展奠定了地震预报的物理基础。地震震源和地震波传播介质的各种参数在强震前的变化早就被当作地震预测的地震学前兆指标，随着地震预测的深入研究，以及我国“十五”台站数字化改造的完成，我们在进一步研究地震时空强分布特征的同时,加强对地震波的运动学和动力学特征的研究,从中提取震源，我们意识到加强对地震波的运动学和动力学的研究，从中提取震源信息，对增强地震预测的物理基础，提高地震预测的水平是十分必要的。地震是地球内部物质运动的结果，这种运动反映在地壳上，使得地壳产生破裂，促成了断层的生成、发育和活动。地震前后的地形变测量和地震波的观测研究等结果确认，天然构造地震是地下岩层的突然错动引起的。发生错动的岩层可称为地震断层。断层活动诱发了地震，地震发生又促成了断层的生成与发育，因此地震与断层有密切联系。地壳中的断层密如织网。实际地震断层的几何形状可能很复杂，但对多数地震，特别是小地震，作为初级近似，总体上可将地震看成是沿一个平面断层发生的突然错动引起的。2、前人对震源机制解的研究历程地震震源处地球介质的运动方式。通常所说的震源机制是狭义的，即专指研究构造地震的机制而言。构造地震的机制是震源处介质的破裂和错动。震源机制研究的内容包括，确定地震断层面的方位和岩体的错动方向，研究震源处岩体的破裂和运动特征，以及这些特征和震源所辐射的地震波之间的关系。对地震震源的研究开始于20世纪初叶。1910年提出的弹性回跳理论，首次明确表述了地震断层成因的概念。在地震学的早期研究中，人们就已注意到P波到达时地面的初始振动有时是向上的,有时是向下的。20世纪的1020年代，许多地震学者在日本和欧洲的部分地区几乎同时发现，同一次地震在不同地点的台站记录,所得的P波初动方向具有四象限分布。日本的中野广最早提出了震源的单力偶力系，第一次把断层的弹性回跳理论和P波初动的四象限分布联系起来。此后，本多弘吉又提出双力偶力系，事实证明它比单力偶力系更接近实际。美国的拜尔利(P.Byerly)发展了最初的震源机制求解法,1938年第一次利用P波初动求出完整的地震断层面解。3、断层及断层面参数3.1、断层参数及分类地震断层通常用断层的走向S、倾角和滑动角三个参数来描述（图2.1）。按目前国际上常用的描述方法，这些参数的定义是：走向S：断层面与水平面交线的方向，但此交线有两个方向，为唯一确定起见，按以下原则确定其中之一为断层的走向：人沿走向看去，断层上盘在右。走向用从正北顺时针量至走向方向的角度S来表示，0S360。 倾角：断层面与水平面的夹角。090。滑动角：在断层面上量度，从走向方向逆时针量至滑动方向的角度为正，顺时针量至滑动方向的角度为负。滑动方向指断层上盘相对于下盘的运动方向。-180180。（仰角：力轴与水平面的夹角（小于90度）方位角：力轴在水平面上的投影线与北方向之间的夹角倾向：节面的上表面的法线在水平面上的投影线与北方向之间的夹角，顺时针量取。）走向S和倾角是断层的几何参数，二者规定了断层的产状；滑动角是断层的运动参数，由这一参数的具体数值，即可描述断层的各种运动类型（图2.2）。有人用断层的倾向代替走向，倾向指下盘断层面向上的法线之水平投影的方向，倾向恒等于走向加90。在地震学中，通常已较少用倾向描述地震断层。按断层节面滑动角判定图2.2滑动角取不同数值所描述的断层类型3.2、断层面上的错动断层滑动是时间和空间的函数，断层面上的错动主要是平行于断层面的剪切位错，描述这些错动的参数如下：（1）地震矩：将地震看成断层面上的突然位错，则形成地震力矩，定义0M如下（2）地震能量：TE （3）应力降：（4）破裂速度：v3.3、震源模型震源机制解（又称地震机理）是指震源区地震发生时的力学过程。鉴于地震机理的研究尚处于探索阶段，目前还属于推断性认识，一般采用各种震源模型进行解析，一种是点源模型，另一种是非点源模型。前者根据点源作用力的不同，又进一步划分为单力偶震源模型和双力偶震源模型；非点源模型也划分为有限移动震源模型和位错震源模型两种。以上震源模型，在分析求解后，提供两组力学参数，一组为断层面走向、倾向和倾角；另一组为最大主应力轴、最小主应力轴和中等主应力轴的方位和产状。根据我国境内150次地震震源机制解，P波初动符号资料确定结果表明，大多数主压应力轴（P轴）和主张应力轴都近水平（T轴），中等应力轴近于直立。地震学的震源理论证明，在均匀弹性介质中，若在一个小的平面断层上发生一个突然的纯剪切错动，则会产生地震波辐射，这样的剪切错动震源产生的远场地震波与在震源处突然有一个双力偶的作用产生的地震波相同。即剪切元位错震源与双力偶点源在产生远场地震波的意义上是等价的。因此，当可将震源近似看成点源时，双力偶点源模型就成为描述发生了剪切错动震源的常用模型。双力偶由一对大小相等、方向相反的力偶组成，（图2.3）是一种合力和合力矩都等于零的集中力系。这样的力系作用于刚体时，不会产生任何运动效果，但在弹性体内部作用，则会使震源区介质产生突然的变形，从而向外辐射地震波。4、震源机制解的测定利用双力偶点源模型，根据地震波观测（或地震前后的地形变测量资料等）求震源模型参数的结果，通常称为震源机制解答，有人称作地震的断层面解。所根据的观测资料可以是P波的初动方向、S波位移的偏振方向、直达P波和直达S波振幅的比值大小，以及P波和S波的波形资料等。双力偶点源模型的独立模型参数只有3个，例如可以是断层面（P波两个节面中的一个）的走向、倾角和滑动角（图2.1），也可以是为确定“震源坐标架”x-y-z（图2.3）相对于“地平坐标架”（例如可分别选为北、东、下三个方向）的空间方位所需要的三zyx个角度值。求解的基本方法是先假定震源模型参数，计算出在给定地球地震波速度结构时，该震源模型在各观测台站所产生的地震波特征，然后与各个台站的实际观测地震波资料进行对比，二者拟合最好的模型参数就作为震源机制的解答。求解过程可以运用反演的数学方法来实现，即选定一种使各个台站的计算结果与观测结果互相拟合的准则，然后使随模型参数变化而变化的准则函数(或目标函数)最优化，即使其最小或最大，而解出模型参数来。4.1P波初动方向法求震源机制解答最简单的方法是根据P波初动方向的观测资料来求解。地表垂直向地震仪记录的初至P波的振动方向，有的向上，即压缩波，记为正号，有的向下，即拉伸波，记为负号。由于介质速度结构的影响，从震源发出的P波一般不是沿直线到达每个台站Si(i=1,2,)的，如图3.1所示意表达的。求震源机制解时，需根据已知的速度结构推算出到达每个台站Si的P波从震源处是沿什么方位Si发出的，即需要将台站Si的记录标在震源球面的相应位置Si上去。震源球面是包围震源的一个球面S，要求球面内的射线不再发生任何弯曲。若将每个台站Si所记录到的P波初动方向都标在震源球面上的相应位置Si上去后，人们发现，对于天然构造地震，只要记录足够多，并且Si在球面上的分布范围足够广，则可以找到过球面中心的两个互相垂直的平面，将震源球面上的正、负号分成4个象限，这两个平面就是上述的双力偶震源的两个节平面。找到两个节平面的空间位置后，震源坐标架的x、y、z轴和P、T轴的空间方位也就知道了。上述求解过程可以通过计算机来实现。两个节面中有一个是断层面，但仅根据P波初动方向记录无法确定哪一个是断层面，还必须根据其他资料，例如现场地质考察资料、余震的空间分布、地震波辐射辐射花样的不对称性（图3.2）和地震波辐射的多普勒效应等，来从两个节面中分析判定实际的断层面，这种判定一般只对大地震才能实现。4.2P波和波振幅比方法利用波初动方向记录反演震源机制解虽简单易行，但也有缺陷。为能将地震波节平面的空间位置约束住，最好要有紧靠节面位置的初动符号观测数据，然而，由于波辐射花样的固有特征，愈靠近节面，波愈弱，初动方向愈不易辨认。再有，由于地表台站布局的限制，观测数据点在震源球面上的覆盖范围经常难以令人满意，当在震源球面上某位置有零散的初动方向数据时，它只能告诉你该点应处在波的正象限或负象限，但不能获得节面与此观测点间的角距离有多大的信息（震源球面上某点观测到的波的振幅大小或波形资料，则含有节面离该点角距离有多远的信息），因而零散分布的初动方向数据对解答起不了多大的约束作用。此外，经常有矛盾的波初动方向读数。上述这些因素常引起解答的不确定性。为解决这些问题，对于较大的远地震，目前多利用在不同方位且具有不同震中距的台站上获得的地震波形记录来反演震源机制解。对地方性的中小地震，有人提出利用从震源向上射出的直达波（Pg）和直达波(Sg)引起的地动位移振幅比求解震源机制的方法5-7。由(2.1)(2.3)式可见，利用直达波（）和SV波（）或SH()的振幅ruuu比，实际是利用它们的辐射花样的比值来求解震源机制参数。振幅比的辐射花样随空间方位的变化比单种波的辐射花样要强烈得多。在力偶平面内（图3.2）振幅比的辐射花样呈瓣分布。从此意义来说，只要有正确的直达波的观测振幅比，且观测值归算到震源球面上后的位置是正确的话，振幅比观测对震源机制参数有较强的约束能力。（但是，振幅比的测量误差和射线在震源球面上出头点位置的误差也容易引起地震节面解的大误差。）此外，对于近震，可近似认为仪器对直达波和波的频率响应是相同的，求振幅比时可以消去仪器影响。但在实用中，也存在一些困难。主要困难是直达波的识别和结构影响的校正。求震源机制解的理论模型分析的皆是直达波和直达S波的振幅比，对于近震，确定初至波振幅常常比较困难，特别要使用垂直向的SV波初动振幅时5,6，测量更困难些。观测量是地动位移，为求震源机制解需要将地动位移的振幅比校正为入射波位移的振幅比。例如，Kisslinger等(1981)的程序5，用的是地方地震的波和SV波垂直向地动位移的振幅比(USV/UP)Z，由于地震较近，用的是上行直达波引起的地动位移，回避了地壳内结构对波的影响，只考虑了自由表面的作用。若令垂直向地动位移振幅比(USV/UP)Z与入射波垂直向位移振幅比(uSV/uP)Z的比值为，图(2.7) 给出了随入射波入射角的变化曲线。观测的振幅比需除以后才能得到入射波的振幅比。由图(2.7)可见，当入射波的入射角为30多度时，值变化剧烈，很难从观测的振幅比求出稳定的入射波振幅比。因为大约35多的角度是SV波入射的临界角(大于此角度后，入射SV在自由面上反射的波出射角变为90，质点呈椭圆振动，振幅随深度衰减)，反射波有相移，由入射SV、反射SV和反射P波三者合成的地表位移振动变得比较复杂。梁尚鸿等(1984)6的程序用的观测量也是垂直向直达SV和P波地动位移振幅比，他们用计算层状介质理论地震图的办法来考虑介质结构的影响，其中也包含自由表面的影响。需要选择符合实际的层状结构才能有效消除结构的影响，否则模型结构的误差会影响机制参数的正确测定。当只用振幅比的大小（振幅比的绝对值）测定震源机制解时5,6，只能求出两个地震节面的空间位置，而不能确定可能断层面的运动特性（旋性），或不能确定波的压缩和膨胀的象限，也即不能确定、轴的具体方位。为确定后者，还必须至少要知道一个波初动方向的可靠读数。如果程序中考虑了观测波和波的初动方向对振幅比取正负值的影响了，则可以求出完整的震源机制参数来。图4.3自由面对垂直向地动位移振幅比的影响因子随直达波入射角的变化。代表地表位移，u代表入射波位移。有人发展只用垂直向SV和P的振幅比测定震源机制的方法5,6，最初是因为单分向的垂直向记录最易得到。如果能获得三分向记录，使用直达SH和P波的振幅比求解震源机制7将会减少结构和自由面对结果的影响，因为结构界面对SH波的影响简单，而SH入射到自由面时，不产生其它转换波，地动位移振幅就简单地是入射波位移振幅的倍。如果使用地动速度的数字记录，而理论分析是用地动位移的振幅比，最好要将速度记录积分成地动位移再求地动位移振幅比，因为一般情况下，地动速度振幅比不一定等于地动位移振幅比。4.3其他方法震源辐射地震波的许多特征都可以用来测定震源机制，例如，利用直达P波或S波振幅的大小8，或在频率域分析直达P波或S波振幅谱的低频幅值9。利用数字地震记录，这些量可以观测得比用笔绘记录更准确些。对于中小地震，利用近距离台站时，对观测振幅或振幅谱一般只需作相对简单的自由面、几何扩散和衰减的校正，而不别作复杂的理论地震图的计算。4.4震源机制参数的表达震源机制解答通常是通过给出震源参数和用图示来表达。作为一例，表3.1给出2001年11月14日中国昆仑山MS8.1地震的震源机制解的参数表达（哈佛大学测定）。力轴的方位（Az）从正北顺时针量度，倾角（PL）为力轴从水平面向下倾的角度。表中的12个参数中，只有3个是独立的，根据其中的任意3个，例如节面或节面的3个参数，或是P轴的2个参数加T轴的一个参数等，可以计算出其他参数来。5个力轴的空间位置用方位角和仰角两个参数来表示：仰角：力轴与水平面的夹角（小于90度）方位角：力轴在水平面上的投影线与北方向之间的夹角将震源球面上表示的数据和震源机制解结果用图示法表达出来时，需要借助于某种将球面上的点与平面上的点一一对应起来的投影法。常用的有乌尔夫网投影和施密特网投影法。二方法的投影平面通常皆是震源球面过球心的水平大圆面（也可以是其它过球心的某个大圆面）（图3.4a），投影点为与大圆面相对的极点。一般只用半个球面的投影，震源机制解通常用下半震源球面的投影，这时投影点为震源球面的上极点z。若某射线穿过下半震源球面的点到达台站，则其在投影图上的投影点为。如果某射线从震源向上发出，穿过震源球面上的点到达台站St（图3.4a），可将射线经震源O向反方向延伸至震源球面上R的对庶点R1，用R1在投影平面上的投影R1来代R的投影。从双力偶点源地震波辐射的空间对称性考虑，这样做是可以的。图4.4震源球面及其投影网。球面上的Q点的空间方位是由过该点矢径的方位角Az和离源角ih确定的，Q在投影网上的投影点Q的位置也是通过这两个角度来确定的（图2.8b），Az从网的正北标记顺时针量（0360），与ih对应的线段OQ的长度d根据投影规则定。对乌尔夫投影，如图2.8c所示，将大圆的半径长度看成为1，dtan(ih/2)3；对施密特投影，该长度应为d=Sin(ih/2)4。乌尔夫投影是一种等角投影，球面上曲2线的交角投影到平面上后保持不变，球面上的圆投影到平面上后仍是一个圆。该投影关系简单，但在表示球