1996年2月3日云南丽江地震震源过程的体波反演

利用IRIS的远震资料,经体波反演得到1996年2月3日云南丽江地震震源过程．震源由两个走向不同的正断层型剪切双力偶组成．这两个子事件在空间上相距15 km,破裂起始时间差为7 s．总地震矩为3.811018 Nm (MW=6.3)．由余震资料估计得到的破裂总面积S约为 720 km2,平均位错为0.18 m．综合考虑反演结果和震源周围的构造环境,第1次破裂可能是沿NNW向的中甸——永胜断裂带上1966年中甸6.4级地震破裂方向上的扩展;之后,第2次破裂开始沿似乎更易破裂的NE向雪山东麓断裂进行．      

具有分形结构强度分布的二维地震断层模拟

提出一种二维地震断层破裂模式，这是在一维弹簧-滑块模式的基础上建立的．利用这种断层的二维向量模型，研究了地震断层的动力学平面-应变的破裂问题．形象、直观地模拟了具有均匀或不均匀破裂强度分布的二维地震断层的破裂成核、扩展、传播直至停止的动力学全过程．阐明了在均匀预应力条件下，断层将获得足够的动量以克服破裂传播路径上的高强度障碍体．破裂锋面亦可绕过孤立的障碍体继续向前传播．整个模拟过程表明:断层破裂过程的停止条件对模拟地震破裂的全过程起着重要作用；我们还研究了断层面上破裂强度分布不均匀条件下断层的动态破裂过程，模拟了具有分形结构强度分布的二维地震断层产生的地震序列．它具有类似自然界地震现象的某些特征．这些特征主要取决于断层面上破裂强度的分布及初始应力降的大小．建立在实验和观测基础上的断层破裂过程模拟研究，为解释地震活动性统计规律的某些特征提供了物理基础.      

山西大同震群序列的震源力学特征和有限元数值模拟

统一研究了1989－1993年期间发生的大同震群活动特点及主要地震的震源机制和破裂过程．为了模似大同震群的震源力学特征，从断裂扩展的角度提出了共轭形断裂系破裂的力学模式，并给出了位移场、应力场和应变能分布的有限元计算结果．此外，还着重讨论了三次大同地震高潮期前记录的地热前兆现象．      

地震活跃期和平静期的模型研究

利用有限元模型模拟一条走滑断层上强震重复发生的过程．在模型中引入断层阀机制以考虑流体孔隙压力变化的作用．用5个参数来表征断层的流变、孔隙压变化及破裂准则等基本特性．通过系统地改变这些参数对大量模型进行计算，初步探讨了地震活动的平静期——活跃期现象与模型参数和加载方式之间的关系．在常速率边界位移载荷作用下，处于中等应力和孔隙压状态下的模型表现出有规律的平静——活跃现象，其周期长度为地震原地复发周期的1/3～1．当模型处于高应力或高孔隙压状态时，平静——活跃现象变得无规律或不明显．在周期性变化的边界力作用下，当力变化幅度达到一定程度时，模型地震活动性表现出与外力变化周期相同的平静——活跃现象．      

用遗传算法反演1994年台湾海峡地震的震源过程及其相关研究

将遗传算法应用于地震矩张量解的反演问题．采用数字化台网P波波形资料及台湾台网P波初动方向资料,研究了1994年9月16日台湾海峡7.3级地震及其邻近地区(北纬21～26,东经115～120)ML5.8共5个地震的矩张量解及震源参数．在5个地震中,1994年9月16日台湾海峡7.3级地震是我国东南沿海地区自1918年南澳7.3级地震以来最大的一次地震．结果表明:这次7.3级地震的矩张量解以双力偶成分为主,是断层面接近NW走向的正断层．断层面的走向与大震前地震的条带分布走向及余震分布显示的破裂图象较相像．震源机制的张应力轴接近水平,近NE走向;压应力轴几乎垂直,近NWW走向．似乎表明,这次地震是受菲律宾板块向欧亚板块北西向挤压的力源控制．但从P轴接近竖直而T轴接近水平看,这是一次强烈拉张性的正断层．其它4个强震震源机制解与7.3级大震相互差别较大．这些地震震源机制解的复杂性表明了海峡地震序列发展过程是比较复杂的．      

滇西地区强震前兆源过程的力学模型与数值摸拟

云南西部地区的强震和大地震主要发生在一些板内大型走滑断裂带上，有一个10～20 km的深度优势层位；地震的破裂基本上是以走滑型破裂为特征．联系云南岩石层结构与动力学背景，我们认为，本区的地震可能有这样一种成核过程: 由于岩石层结构的层次性、非均匀性，在整个岩石层板块构造运动背景场中，可能会在岩石层的中深部形成局部剪应变集中区或滑动区，尔后这个剪应变集中区或滑动区沿着断层带边界向地表滑移扩展．当滑动由深部向浅部扩展时，滑动峰将遇到断层面上的最大抗剪强度区而被阻碍闭锁，形成一个地震活动空区．随着构造荷载的逐渐增加，滑动最终要向前扩展导致整个闭锁区的失稳破裂，产生大地震．给出了描述这个地震过程的一个简单的近似积分方程；数值结果表明，该模型具有一个向不稳定非线性加速发展的阶段，这个阶段对地震前兆的形成可能有重要意义．      

微震分析水压致裂的破裂过程

水压致裂技术广泛用于增加石油和天然气产量和开发地热资源．用水压致裂法在地下深部形成破裂带系以储存固液态废物是工业技术与环境保护科学相结合的最新进展．了解水压破裂的力学特性和破裂带的几何形状对于资源开发和贮存废物都十分重要．水压破裂常伴有大量极微震发生．分析伴生微震的时空强分布可得出水压破裂带的几何参数及破裂过程的运动学和动力学参数．利用波形相关分析和时空搜索定位法，我们高精度地确定了157个水压破裂伴生微震的震源位置．微震的空间分布图清晰地刻画出水压破裂带的空间尺度和取向．由微震的时空分布的变化推断出破裂带的扩展方向和速率．应用经验格林函数法分析孪生地震对得到较大微震的震源时间函数和震源参数，如地震矩、破裂半径和应力降．震源时间函数的方向性变化表明破裂向西北方向传播，与破裂带扩展方向一致．微震应力降的较大变化反映出水压破裂带上力学性质的显著不均匀性．      

地震前兆场物理模式与前兆时空分布机制研究(二) ——强震孕育时应力、应变场的演化与地震活动、地震前兆的关系x,auto,auto,415px);}style> hf=httprztl.com >business

分别研究构造块体内有多个震源体和单个震源体孕育时，应力集中过程、时空分布规律及其与地震活动、地震前兆的关系．结果表明:① 当区域应力场增强到相当水平时，破裂开始在某个或某些块体内发生，区域应力场将随着时间而变化，强震的成组性是这个动态应力场演变的结果．在此过程中，多个高应力集中区的存在与发展，有可能导致震前异常与地震关系的多样性；② 孕震块体的应力变化过程一般显示出非线性，在经历了长期弹性变形后，从进入非弹性变形阶段到主破裂前，可能经历不止一次的非弹性变形与断层软化过程；与此相应，孕震块体中的应力、应变场要表现出复杂的时空演化图象；受应力、应变控制的地震前兆场必然表现出多方面的复杂性.      

1975年辽宁海城MS7.3地震的前震序列分析

本文讨论了地震成核问题．加速是成核相的属性，也是地震失稳的一个必要条件．若以方程d／dt=C／(tf-－t)n表征成核过程的加速特性，则可将其概括为地震释放速率比例于失稳破裂剩余时间的负幂．根据这一原理，对1975年辽宁海城MS7.3地震的前震序列进行回顾性分析，表明主震的失稳破裂时间tf和震级可成功地被估算，但要求所测相关系数r2达极大；且主震时间的估算误差一般小于所用最后一个数据点与主震之间剩余时间的一半或更少．      

含硬包体试样微破裂图象演变过程的实验研究

含硬包体试样在双轴加压下的微破裂图象为:加压初期，微破裂先在试样周围活动，中部形成空区；随1增加，微破裂空区逐步变小；1再增加，微破裂在原空区内活跃，而原活跃区的微破裂则减少；临近主破裂时，微破裂逐步形成条带；条带上的微破裂来回跳跃；微破裂空间分维Ds从小到大，临近主破裂前转为降低．这些特征在以往破裂实验中很少提到，但与强震前地震活动图象的长期活动特征有一定的相似性．文中以澜沧-耿马地震为例进行了具体的对比分析．      

岩石破裂孕育过程中的几种Q值变化

采用脉冲传播谱振幅比法、能量法和尾波法分别测定片麻岩破裂孕育过程中的Ｑ值变化．结果表明:谱振幅比法所得Ｑ值呈现相对平稳值——升高至甚高值(极大值)——回降——破裂的变化形式；能量法所得Ｑ值变化形态与上述相反，即平稳值——低值——回升——破裂；而尾波法所得Ｑ值，主频Ｑ值呈现相对高值——下降至低值(极小)——回升——破裂，与能量法Ｑ值变化大体一致．而对于高频(大于主频)Ｑ值则呈现平稳值——升高至极大值——回降——破裂．同时尾波振幅谱还表明:应力达70％破裂应力后，首先是高频Ｓ波迅速减少(QC)变大)，临近主破裂时主频Ｓ波则大幅度减少(QC回升)，这可能是两种尾波QC变化形态不同的原因所致．而P，S波(初至波)振幅谱也表明，随着裂隙出现，Ｓ波频带变窄，高频成分迅速减少．即S波谱变化不同．这使采用不同方法所得的QS的变化特征不同．      

用经验格林函数方法研究澜沧-耿马MS=7.6地震的破裂过程x,auto,auto,415px);} h

运用经验格林函数反褶积方法,用中国数字地震台网(CDSN)6个台记录的MS=6.3余震的地震图作为经验格林函数,研究了1988年11月6日中国云南澜沧-耿马MS=7.6地震的破裂过程．反褶积结果表明,MS=7.6地震是一次比较简单的事件．分析了破裂的方向性,得出MS=7.6地震的破裂开始于微观震中,然后沿西北、东南两个方向近似对称地扩展,破裂总长度为70 km,持续时间19 s,破裂速度约为2 km/s．      

震源孕育模式的数值模拟研究

在充分考虑地质体材料力学性质非均匀性特点的基础上，运用岩石破裂失稳过程有限元分析新方法(RFPA)，以不存在明显孕震构造特征的包体模式为例，模拟研究了挤压边界条件下的震源孕育模式，包括断层的形成过程和地震活动特征，演示出孕震体由变形、破裂到宏观贯通的地震全过程以及它们的时空强基本规律．计算所得到的应力、变形和微震空间分布，形象地给出了变形局部化、震源迁移、膨胀或隆起，弹性回跳和X型共轭变形带等直观图象．      

青海省共和7.0级地震前兆的数值模拟

在共和地区简化的介质模型上，应用固液两相介质的本构关系，模拟计算了1990年共和7.0级地震的前兆特征．结果表明，模拟得到的介质状态分布、扩展变化与共和地震前可靠异常的时空分布规律大体相似．由此得知，共和地震的孕育经历了大面积弹性变形——早期的非弹性膨胀——西宁周围地区的应变软化与弹性恢复——震源附近高密度地体的非弹性膨胀——发震等过程．可以看出，它既不是一个简单的裂隙发育贯通过程，也不是孕震区非弹性体积膨胀、进水过程，而是介质多种状态此起彼伏的复杂进程，其具体程式受震源及周围地区的地质构造环境所决定．由此得出，不同地震孕育的地质构造环境不可能相同，其前兆显示的时空分布规律也必然不同．      

用数字化宽频带波形资料反演共和地震的震源参数

1990年4月26日，在青海省共和县境内发生了一次MS＝6.9的地震．随后于1990年5月7日、1994年1月3日和1994年2月16日在同一地区相继发生了震级分别为MS＝5.5，6.0和5.7的余震．用频率域矩张量反演方法得到了上述地震的矩张量．结果表明，这4次地震的震源机制总体上比较接近，是以逆冲为主的、走向NWW、倾向SSW的断层，其标量地震矩M0分别为:MS＝6.9地震，9.41018 Nm；MS＝5.5地震，8.01016 Nm；MS＝6.0地震，4.91017 Nm；MS＝5.7地震，2.91017 Nm．反演结果还表明，这几次地震的震源过程明显不同，MS＝6.9地震比较复杂，它主要是由两次规模相当的事件构成的，第一次事件从0 s开始至12 s结束，地震矩为4.71018 Nm；第二次事件从31 s开始至41 s结束，地震矩为2.51018 Nm；另外，在上述两次主要事件之间(12～31 s)，还有一次较小的事件发生，地震矩为2.11018 Nm．与MS＝6.9主震相比，3次余震的震源过程比较简单，其震源时间函数呈简单脉冲形状，表明这几次余震都是主要由一次连续的破裂构成．MS＝5.5地震的上升时间为4 s，持续时间11 s；MS＝6.0地震的上升时间为6 s，持续时间16 s；MS＝5.7地震的上升时间约为6 s，持续时间13 s．      

1989年6月9日石棉5.2级地震后横波分裂的时间变化c

1989年6月9日石棉5.2级地震后，在震中附近架设了两个三分向数字地震仪，从1989年至1992年进行了4年的观测．资料分析表明，无论是余震，还是周围地区的微弱地震，在横波窗内的地震记录上都存在横波分裂．横波分裂随时间改变．在余震活动时期，快横波偏振平均为N18W，与5.2级地震主破裂面走向一致；1990年以后，平均为N46W，转向与区域最大水平主应力方向一致．1990年以后两个分裂横波到时差为1989年观测的一半左右．从两个台的资料分析中得到的结果相同，这意味着观测到的横波分裂的时间变化可能与石棉5.2级地震有关．      

1995年10月24日云南武定MS＝6.5地震的震源机制解

本文用地震矩张量反演方法反演了1995年10月24日云南武定MS＝6.5地震的震源机制，得到:两主轴方向近水平，主压力轴位于北偏西30方向，方位角105(西偏北15)的共轭面可能与实际的地震破裂面一致．      

根据环境应力场寻找中国大陆地区的潜在震源区

地壳介质是一种里面存在着许多裂纹的非完整介质．断裂力学的研究结果表明，这种介质的强度远低于完整介质的强度．笔者认为，地震的发生是由于地壳中裂纹在合适的条件下由稳态扩展逐渐发展到失稳扩展的结果，地震破裂即为断裂力学所阐明的低应力快速破裂现象． 快速破裂的出现必需要有局部应力和应变的高度集中．根据Knott(1973)对裂纹尖端进行的弹塑性应力分析表明，裂纹尖端的应力状态对裂纹的扩展起了极其重要的作用．由于应力集中，在裂纹尖端必然形成一个塑性区．塑性区中的位移已不是弹性位移，而是塑性位移．=a20/2y，这里，0是作用剪应力，它相当于地震发生时的初始应力或构造应力，a是裂纹半长，是刚性模量，y是介质剪切屈服强度．而地震位错恰恰是发生在形成塑性区的裂纹末端．于是我们提出了一个关键性的假定，设位错D(1，t)为文中的公式(5)，地震最大位移为Dmax=L20/4g，式中，L是断层长度．若刚性模量取上地壳的数值，=33 GPa，介质剪切屈服强度y取实验室给出的平均值y=30 MPa，那么，根据Dmax和L的观测值，利用此公式就可以估计大震的构造剪应力值．计算结果表明，全球的大震构造剪应力值大多数在5——20 MPa 之间，有地区差异，平均为10 MPa． 笔者在进一步研究了震源谱的性质后，又导出了体波震级mb对0和地震矩M0的依赖关系(即文中的公式(11))，此公式提供了直接利用地震观测资料大量计算构造剪应力值的可能性．我们认为，构造剪应力场是控制地震发生的主要因素，出现高应力值地震的地区是容易发生强震（MS6.0）的地区，即地震危险区．我们对1987年以来全国所有mb3.8的地震计算了它们的剪应力值，据此，划出了我国大陆地区强震危险区并给出了震级范围． 1992年4月——1994年1月31日，中国大陆地区共发生9次MS6.0的地震，其中8次均落在我们事先划定的危险区内，只有1次在危险区外．实践证明，这种新方法作为中期强震预报，是一种有效的预报手段，有着良好的物理基础和令人鼓舞的发展前景，值得深入探讨和研究．      

1996年2月3日云南丽江地震双重破裂的初步估计及其相关问题研究x,auto,auto,415px);}style> hf=httprztl.com >business week launch

利用同时反演两个点源震源机制的原理,估计了1996年2月3日云南丽江主震的震源机制,得到丽江主震由两次主破裂组成,两次破裂的时差约12 s,相距约26 km．第1次破裂形成近南北向的张性倾滑-走滑断裂,第2次破裂的机制需进一步研究．介绍了同时反演两个点源机制的方法及直接分析数字化波形记录的结果,讨论了一些相关问题．认为1995年10月的武定地震和本次丽江地震是同一力学过程在不同时空点上的表现,其发生顺序与动力传递的方向有关．      

用经验格林函数方法从长周期数字波形资料中提取共和地震的震源时间函数

1990年4月26日，在青海省共和县发生了一次震级MS=6.9地震．这次地震之后，于1990年5月7日、1994年1月3日和1994年2月16日在原震区又发生了MS分别为5.5，6.0和5.7的余震．我们分别以这几次余震的中国数字地震台网的长周期记录作为经验格林函数对该地震的长周期记录进行反褶积，提取了该地震的远场震源时间函数．无论以哪次余震的记录作为经验格林函数，反褶积的结果都相当一致．结果表明，发生于1990年4月26日的MS=6.9地震至少由两次规模相当的事件构成，两次事件的发震时刻相差约30 s．第一次事件的持续时间较短，约12 s，震源时间函数的上升时间约5 s；第二次事件的持续时间较长，约17 s，上升时间约8 s．分析分别从P波和SH波提取的震源时间函数，我们注意到，从SH波中提取的震源时间函数比从P波中提取的震源时间函数复杂，表明这次地震的震源过程除了上述两次规模相当的事件外，还存在规模较小的事件．这一结果与我们由矩张量反演得到的结果以及用宽频带波形资料借助于经验格林函数方法所得的结果一致．我们以较小余震对较大余震做反褶积得到较大余震的震源时间函数．结果表明，发生在1994年1月3日的MS=6.0和1994年2月16日的MS=5.7地震的震源过程相当简单，其震源时间函数为一简单脉冲．对于1990年4月26日MS=6.9地震，我们从不同台站的资料中得到的震源时间函数具有方位依赖性．然而，这种方位依赖性在3个余震的震源时间函数中几乎没有显示．从不同台站或同一台站的不同震相求得的震源时间函数的幅度存在差异，但震源时间函数在其持续时间上的积分却相当稳定，亦即从任何一个台站测得的标量地震矩相当一致． 以3个余震的记录作为经验格林函数得到的1990年4月26日地震的相对标量地震矩分别为22(以1994年1月3日MS=6.0地震的记录为经验格林函数)、26(以1994年2月16日MS=5.7地震的记录为经验格林函数)和66(以1990年5月7日MS=5.5地震的记录为经验格林函数)．据此推算，MS=6.0地震与MS=5.7地震的相对标量地震矩为1.18；MS=6.0地震与MS=5.5地震的相对标量地震矩为3.00；MS=5.7地震与MS=5.5地震的相对标量地震矩为2.54．这些结果与我们从资料中直接得到的相对标量地震矩(分别为1.15，3.43和3.05)一致，表明了作为提取震源时间函数的经验格林函数方法副产品的相对标量地震矩是地震大小的一种很稳定的度量.