产生有限差分合成地震图的最佳震源时间函数

震源时间函数在产生有限差分合成地震图中起着十分重要的作用。鉴于此, 我们考虑了5种不同的时间函数的振幅谱。结果表明, 当时间函数的高频信号比主信号携带更多的能量时, 伴随着主信号成分将出现我们所不希望的拖尾现象。本文说明, 在所考虑的几种用于产生有限差分合成地震图的震源时间函数中, 形如f(t)=te-2的时间函数是最佳的。由具有一个覆盖层的半无限空间弹性地球模型, 来产生有限差分合成地震图证实了这个结论。     

新丰江地区部分地震的震源时间函数

广东省数字遥测地震台网建成以来积累了较丰富的数字化地震记录，作者利用这些数据开展震源研究，得到新丰江地区的部分地震震源时间函数的初步结果。得到结果的3组地震它们各自较大地震的震级分别为ML1．8、2．3和3．0级，对应的震源时间函数上升时间分别为0．07s、0．09s和0．11s，函数形态较为简单。     

震源时间函数与震源破裂过程

简要介绍了现代数字地震学的两个重要概念。一个是震源时间函数，另一个是地震震源的时空破裂过程。首先，从地震断层位移表示定量出发，介绍震源时间函数和震源破裂过程的基本概念。然后，分别介绍两种从远场地震记录中提取震源时间函数和获取有限断层面上时空破裂过程图像的反演方法。     

由震源时间函数估算地震辐射能量

用实际地震观测中得到的震源时间函数讨论了地震辐射能量的测定问题。指出从辐射能量测定的角度看，实际地震的震源时间函数可以简化为三角形函数，同时一些复杂的地震可用两个或更多的三角形函数来表示。用三角形函数作为实际震源时间函数的简化，估算了一些地震的辐射能量。将NEIC的辐射能量测定结果与用震源时间函数得到的辐射能量的结果进行比较，结果表明两种估计具有比较好的一致性。     

测量震源辐射的谱震级和谱地震图

为测量地震强度所做的努力已有100多年了。1935年,里克特用伍德-安德森地震仪记录的南加州地震建立了地方震震级标度M_L,所用地震波的周期范围为0.1～2s。1945年,古登堡(1945a)根据远震记录周期为(20±3)s的水平向面波建立了面波震级标度M_s;古登堡(1945b)还分别用平均周期为5s和10s的P波和S波引进了体波震级标度m_b。所有这些标度都是基于某些特定的周     

用经验格林函数方法从红山地震台数字化资料提取隆尧地震的震源时间函数

2002年4月22日在邢台老震区隆尧县境内发生了一次ML4．8地震，这次地震邢台地区震感强烈，随后又发生15次小的余震。我们从这些余震中任意选择了5次余震分别作经验格林函数，提取了该地震的震源时间函数，目的是从这个地震来反映邢台震区的地质构造和震源情况。结果表明：无论以哪次余震的记录作为经验格林函数，结果都相当一致。发生于2002年4月22日的ML4．8地震只是老震区断层的单次破裂，此次事件的持续时间较短，只有0．7s左右，上升时间约为0．5s。为了验证所提取的震源时间函数的稳定性和可靠性，作者用5个小余震的UD，EW，NS方向分别作经验格林函数处理，结果相当一致；同时用全波、P波和SH波进行提取震源时间函数的对比，结果表明，用全波和SH波效果较好，所得结果非常一致，而用P波提取的震源时间函数效果不很理想。即同一台站的不同震相求得的震源时间函数的幅度存在差异。     

远震体波波形反演获得的震源时间函数

地中海东部大地震成核问题的一个有趣的现象，是它们的震源时间函数的急始特性，使用同一反演方法确定了发生在世界其他一些地区破坏了地表的浅源大地震的震源时间函数都呈现出类似的特性，这些地震的断层滑动类型包括：正断层、走滑断层和冲断层型。这些急始型震源时间函数可能反映了破坏障碍体的困难性，即在破裂成核的动力断裂过程完成后，断层的局部地区仍未破坏，也可能反映了通常发生在大地震破裂期间的某些特殊的物理过程（如蠕变不稳定性）。发生在类似地质构造环境中，且震源机制也类似的地震，反映出共同的模式。     

高铁地震震源子波时间函数及验证

基于纵横波产生机理并结合高铁列车及铁轨结构,本文给出了高铁列车产生地震纵波及横波震源（点源及等效源）子波的时间函数;通过实际高铁地震数据及合成高铁地震数据分析,验证了高铁地震震源子波时间函数的合理性;在高铁列车匀速、加速及减速运行情况下,利用高铁地震点源子波合成记录分析了高铁地震震源子波的频谱特征,为实际高铁地震数据处理解释奠定了良好的基础.

     

用经验格林函数确定中小地震的震源时间函数

采用正则化方法,并以小地震近似地代替通常由理论计算出的格林函数,结合近场加速度资料反演了1985年4月18日云南禄劝 Ms=6.1地震4个余震的震源时间函数.结果表明,较大余震除了震源持续时间较长外,其破裂过程也较为复杂,具有明显的阶段性.而较小余震的震源时间函数则近似为一脉冲函数.根据布龙震源模式,我们估计了这4个余震地震断层的破裂速度和断层面上的平均质点运动速度,结果表明,地震断层面的破裂速度小于地震波传播的横波速度;断层面上的平均质点运动速度为几至十几 cm/s,基本保持在同一水平.      

利用视震源时间函数反演尼泊尔M_S8.1地震破裂过程

以2015年4月26日MS7.1余震为经验格林函数事件,利用全国和全球的宽频带记录提取了2015年4月25日尼泊尔MS8.1地震的P波视震源时间函数和Rayleigh波视震源时间函数,并通过联合反演这些视震源时间函数获得了这次地震的时空破裂过程图像.无论是P波视震源时间函数还是Rayleigh波视震源时间函数都呈现出很强的方位依赖性,表明震源断层具有相当的尺度且破裂朝东南方向扩展.时空破裂过程图像清楚地证实了这一特征,并更清晰地显示,破裂几乎是纯粹的单侧破裂,从破裂起始点开始,沿断层面向东南方向扩展~100km,同时沿断层面向深部扩展~80km,形成~125°的破裂优势方向和~5.8m的最大位错.地震的破裂时间历史相对简单,呈非间断性扩展,持续时间约50s.     

芦山地震有限断层混合震源模型模拟

建立了基于有限断层的芦山地震的混合震源模型,突出破裂面上滑动量分布中凹凸体的主导作用,同时借助k平方模型表达对震源复杂性认识尚不够深入所带来的随机性.混合震源模型中全局震源参数均值与局部震源参数均值的确定,主要依据统计的半经验标定公式,进一步借助截断正态分布将区域地震、地质构造、地震活动性提供的参数限值结合起来,生成最终的震源参数值.基于上述思路,本文共建立了30组芦山地震混合震源模型,进而以芦山地震8个近场观测台站为试算点,借助地震动反应谱残差评价的方法,从30组模型中选定了表达“平均”特征的混合震源模型.     

用余震作为经验格林函数从GDSN长周期波形资料中提取1999年集集地震的震源时间函数

1999年9月20日17ｈ47minUTC（北京时间9月21日01ｈ47min）中国台湾省南投县集集镇发生Mw7.6地震，随后又发生许多中强余震．我们从这些余震中选择了Mw6.1和Mw6.2两次余震分别作为经验格林函数，并从包括IRIS，GEOSCOPE和CDSN在内的全球数字地震台网的长周期波形资料中提取了主震的震源时间函数．对于发生在9月22日的Mw6.1余震，有26个台站的78条记录的97对震相清晰可用；对于发生在9月25日的Mw6.2余震，有24个台站的72条记录的81对震相清晰可用．对于每个台站，从Ｐ和Ｓ两种震相中提取震源时间函数，分别称之为Ｐ波震源时间函数和Ｓ波震源时间函数．这样，利用两次余震从长周期波形资料中共提取了178个震源时间函数．从这些震源时间函数可以看到，除在特殊方位的几个震源时间函数由于主震和余震的震源机制的差别引起较大的畸变外，大多数台站的震源时间函数都具有相似性；同时，其形状随方位发生一定程度的变化．这两种特征都是对所提取的震源时间函数的稳定性和可靠性的客观反映．通过分析确认，这次地震由两次子事件构成，总持续时间约26ｓ．第二次事件比第一次事件晚大约7ｓ，第一次事件的矩释放率比第二次事件大15％．      

终止相的理论地震图-兼论唐山地震震源破裂过程

本文首先较为系统地回顾了终止相研究的历史、进展,对震源破裂过程产生终止震相的物理机制进行了讨论,并详细讨论了主事件分析法以及合成包含终止相的理论图的方法。在此基础上应用主事件分析法处理了唐山地震的实际资料,给出唐山地震的震源破裂过程的大致时空图象,同时计算了包括终止事件的合成图。结果表明,考虑破裂过程的终止事件可以较好地拟合P波观测图,从而进一步确认唐山地震震源破裂过程包括两次终止事件的可能性。此外本文还对唐山地震的震源环境以及调整单元的让位情况进行了讨论。     

陕西及周邻地区几个地震丛序列的震源谱振幅相关研究

选取自2008年汶川地震以来山西霍州、陕西宁强、宁夏银川及甘肃古浪地区附近发生4个中小地震丛集,采用谱振幅相关分析法计算了各个地震丛序列的谱振幅相关系数随时间的变化,分析了各序列地区发生强震的危险性。结果表明,山西霍州序列谱振幅相关系数较高,表明霍州地区附近应力水平较高,后续发生强震的危险性较大。其余3个地区附近地震丛序列谱振幅相关系数较低,表明这3个地区应力水平不高,后续发生强震的危险性较小。     

考虑震源破裂过程的青海玉树地震震动图研究

为了在缺少实时台站数据的情况下,研究震后地表地震动分布,本文尝试通过对震源做出某些约束,以提高震动图的精度.由于青海玉树地震的发震断层近似直立(倾向83°),震源破裂过程在时间和空间上呈明显的分区特征,分别考虑了地表破裂的线源发生模型和主震震级分解方法,用考虑了场地效应的震动图快速生成方法生成震动图.将两种方法生成的震动图与调查烈度比较.结果表明,基于震源破裂过程对主震震级进行分解,综合考虑各子事件生成震动图的方法,对震后的应急决策和灾情的快速评估是可行的.     

用功率谱时间函数识别地震前兆异常

多年的前净观测和地震预报实践表明，与前兆观测资料的噪声水平相比，中强地震的前兆信息多属弱信息，信息识别较为困难，用功率谱分析方法，将前兆观测时间函数分解为不同周期的谱密度随时间变化的多个时间函数，通过震例分析，探讨在强干扰背景下识别中强地震前兆信息的可能性。从识别中短期地震前兆这一基本目的出发，采用“数据窗”滑动采样的计算方法，满足不同周期对样本长度的要求，实现了频率域的计算结果在时间域具有一定密     

岩石的双剪摩擦及震源时间函数的实验研究

本文介绍了岩石的双剪摩擦滑动实验。在实验宁测量了不同正应力条件下的粘滑应力降和位移值,求出了静、动摩擦系数;用瞬态波形存储器记录了粘滑过程中位移随时间的变化曲线,得出了震源时间函数,它的形状与布隆提出的震源时间函数极为相似,两者皆以指数形式下降,它的数值与野外实测结果大致相近。为验证实验数据的准确性,本文还对粘滑的理想模型进行了分析和计界。理论、实验和野外观测的结果大体一致。     

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

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个余震的震源时间函数中几乎没有显示．从不同台站或同一台站的不同震相求得的震源时间函数的幅度存在差异,但震源时间函数在其持续时间上的积分却相当稳定,亦即从任何一个台站测得的标量地震矩相当一致．\r\n 以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)一致,表明了作为提取震源时间函数的经验格林函数方法副产品的相对标量地震矩是地震大小的一种很稳定的度量.     

利用改进的粒子群算法反演视震源时间函数

视震源时间函数的提取是研究震源参数的重要途径. 本文提出了利用改进的粒子群(PSO)算法反演视震源时间函数的方法, 以水平线方法得到的结果作为PSO算法的初值, 并对PSO算法的惯性因子和学习因子进行改进, 提高计算效率. 采用改进的PSO算法对模拟数据进行了反演计算, 并与映射Landweber反褶积(PLD)方法和遗传算法(GA)进行了对比分析. 结果表明, 相对于PLD方法, 改进的PSO算法反演结果与真实结果误差更小; 相对于遗传算法, 改进的PSO算法计算效率提高了5倍以上. 最后, 利用改进的算法对2005年10月8日巴基斯坦克什米尔M_W7.6地震的P波视震源时间函数进行了提取, 结果表明此次地震P波视震源时间函数在25 s之内, 震源沿西北向破裂. 该结果与张勇等的结果一致.