利用震源扫描方法计算2008年大柴旦6．3级地震震源破裂过程

2008年11月10日在青海大柴旦地区发生Ms6．3级地震。运用甘、青区域数字地震台网观测资料,采用震源扫描方法对此次地震进行亮度函数归一化扫描分析研究,得到研究区域内分层归一化亮度函数的空间成像结果,最大亮度空间分布与青海省野外考察结果基本吻合。震源扫描方法不需要计算理论地震图,利用包括相对振幅和到时在内的波形信息,无需事先假定震源的几何参数,最大限度地利用波形数据信息,得到地震的整个空间分布,提高了获取中强地震震源运动空间分布信息的速度为震后快速趋势判定、应急救援提供技术支持。     

模糊数学方法在识别地震和确定震源参数中的某些应用

本文把模糊聚类分析方法用于确定P波到时不十分清晰的情况下,震源参数计算前的慢度(以及P波到时)、观测点之间的距离、波的周期、F—P时间以及最大振幅等特征量。经过模糊聚类分析之后。就可用最大一组的数据来计算震源参数,在各组数目相同的情况下,则用初至P波到时最早的台所在的一组数据来进行震源参数计算。     

用震源扫描算法（SSA）研究列车源的运动

用震源扫描算法（Source-Scanning Algorithm,SSA）对2004年5月山东兖州试验时地震仪记录到的列车产生的一次振动波形进行研究.结果表明,在3个短暂时间段内,我们可以清楚地看到列车运动所产生的鸟巢形波形的精细结构.通过分段分离的方法,我们得到了3张有关震源分布的图像,图像显示列车正在由NW往SE行走穿过本区.扫描结果表明,列车震源经过适当预处理可以当作许多分散的小震源来分别处理,也显示了SSA方法处理象列车这种无法识别震相的数据波形具有的潜力.     

领域算法对震源特性之描述

对震源特性完整描述需要说明震源位置和震源机制。采用邻域算法经两个步骤就可实现非线形反演，邻域算法是一种参数空间中的直接搜索法，它能优先采集具最小数据拟合误差的区域，仅采用两个控制参数，无需导数计算。第一步是采用一个四维搜索空间的集中搜索方案来进行震源定位，第二步是对P和S波的早期波形反演来精确定位震源深度及提取源机制信息。这第二步就是采用一个矩张量表达式扫描一个八维空间。就震源和源机制反演而言，邻域算法（NA法）仅用几个测站的资料，就给出了快速有效的结果，正如对新疆南部一次事件处理所演示的一样。     

山西地区震源机制一致性参数时空特征分析

利用山西数字地震台网记录到的中小地震波形资料,采用层状介质中点源位错模型的广义透射系数的快速算法和理论地震图拟合直达波最大振幅比来求取小震震源机制解的方法,计算了2001年—2012年山西地区281个中小地震震源机制解,根据优势分布得到山西地区现今平均构造应力场分布,计算震源机制解应力主轴与相应的区域构造应力场主轴之间的震源机制一致性参数,分析一致性参数时空分布与中等地震活动的关系。结果显示,2001年以来山西5次ML≥5.0地震前均有震源机制趋于一致性现象,震中分布在一致性参数低值附近或高低值过渡区域。     

用震源扫描算法(SSA)进行微震的定位

引入了一种新的方法震源扫描算法(Source-Scanning Algorithm, SSA)来对地震定位和地震破裂面进行研究, 这种方法充分利用了数字地震波形资料, 在不用精确拾取到时和计算理论地震图的情况下达到比较理想的定位效果。 首先, 文中阐述了2004年由Kao等人发展的SSA方法的原理; 其次对扫描算法的分辨率进行了理论上的模拟, 用5点试验证明了SSA的很强的分辩能力; 最后将SSA应用于南北地震带微震的定位试验中去, 发现它能够在没有人工干预的情况下较好地定出地震的位置。 理论模拟和实际定位结果都证明, SSA是普通微小地震定位的一种好工具。     

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

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

滇中地区震源机制一致性参数时空分布与强震活动

利用云南数字地震台网记录到的地震直达波Pg、Sg振幅资料,采用层状介质中点源位错模型的广义透射系数的快速算法,计算理论地震图的Pg、Sg最大振幅比值,将其与观测值拟合,反演得到1999年至2009年11月14日滇中地区206个中小地震的震源机制解,计算震源机制解应力主轴与相应的区域构造应力场主轴之间的一致性参数,分析震源机制一致性参数时空分布与强震活动的关系.分析表明,强震发生前3年至数月,强震震源区附近出现多个震源释放应力场与区域应力场一致或接近的中小地震,强震就发生在中小震震源机制一致性参数低值分布区内或其边缘附近.     

用宽频带区域地震图确定震源参数

近期发展的震源反演方法并没有充分利用区域地震图的宽频带性质,原因是用平行层模型计算的格林函数和真实地壳的复杂传播现象不同。这种反演通常去掉了短周期成分,并且只反演记录中选定的部分。在本文中,我们介绍一种估算震源参数的方法,它可以较好地应用整个宽频带记录,而且只采用目前仅有的不完善的格林函数即可。这个方法通过分别拟合格林函数的各个部分来使得记录图与合成图主要的地壳到达波到时差别不大。确定震源参数的原     

利用数字化速度波形和仿真位移波形资料求解震源机制解的比较

对2007年1～7月云南发生的19次ML3～4.6地震的实例分析表明,数字地震仪的速度记录波形与仿真位移波形的形态相似,振幅分布比较一致,最大振幅比值(ASV/AP)也基本相等;由数字地震仪速度记录振幅资料和仿真位移振幅资料测定的震源机制解具有显著的一致性.在采用Pg、Sg波最大振幅比测定震源机制解时,直接量取数字地震仪速度记录相应波形的最大振幅值,获取相应振幅比数据,这个方法是可行的.     

基于地震波反演研究南天山中西段的震源深度

南天山中西段位于塔里木块体与西昆仑之间,地震频发。本文利用基于波形反演的CAP方法对南天山中西段地震震源深度进行了反演,得到的震源深度与用到时定位得到的震源深度进行了对比,发现用到时定位得到的震源深度整体偏浅。同时,南天山中段山前地区震源深度比南天山中段中部地区震源深度浅,南天山由西往东呈现震源深度变浅的趋势。南天山中西段震源深度整体分布在15—23km,有明显的地震活动深度下界。     

2004年东乌珠穆沁旗地震震源参数研究

利用波形反演方法, 研究并确定了2004年3月24日发生在内蒙古东乌珠穆沁旗与西乌珠穆沁旗交界处地震的震源机制及深度。 近震波形对震源机制有很好的约束, 而远震体波对于震源深度分辨率很高, 因此综合利用了近震和远震波形数据。 对于近震宽频带波形记录, 用CAP方法反演震源机制解。 利用远震直达P波、 pP波和sP波的到时能较好的确定震源深度。 综合利用近震宽频带波形和远震体波波形记录, 反演得出2004年03月24日地震为逆冲型, 震级为M_W5.3。 其两个节面分别为: 节面I的走向、 倾角、 滑动角分别为147°、 22°、 87°, 节面Ⅱ的走向、 倾角、 滑动角分别为 330°、 68°、 91°; 震源深度为(12±2) km, 是一个典型的浅源上地壳地震。 与前震序列对比, 节面Ⅰ是可能的发震构造。     

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

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

基于重构激发信号的电磁式可控震源相关检测参考信号方法研究

针对电磁式可控震源地震数据的相关检测,研究发现,在地下结构复杂、基板-大地耦合不佳时,常规方法——基于震源控制信号或基板附近信号作为参考信号检测得到的地震记录中,存在子波到时误差和虚假多次波问题.本文分析了上述问题的理论原因,并提出基于重构激发信号的相关检测参考信号方法(Correlation Detection Reference Signal Based on the Reconstructed Excitation Signal,CDRSBRES).首先,利用直达波与其他地震波到时不一致的特点,从震源基板附近信号中分离、提取直达波.然后,利用直达波重构震源激发信号并作为参考信号对地震数据进行相关检测.最后,应用谱白化技术提高检测结果质量.数值模拟研究表明,重构激发信号与理想激发信号的相关系数为0.9869,达到高度线性相关,CDRSBRES方法检测的地震记录在子波到时和波形特征上均与模型相符.随后,在某金属矿区开展了可控震源对比实验.与液压式可控震源MiniVib T15000检测结果相比,电磁式可控震源PHVS 500的检测结果中:基于震源控制信号的检测结果存在子波到时误差约0.012s,对应垂向精度误差约11.16m;基于基板附近信号的检测结果部分区域出现虚假多次波,信噪比降低;而CDRSBRES方法的检测结果子波到时误差约0.001s,对应垂向精度误差约0.93m,波形特征一致,相同区域无虚假多次波.综上,本方法适用于电磁式可控震源地震数据的高精度检测,尤其对于地下结构复杂区域的高分辨率地震勘探具有重要意义.     

青海地区中强地震震源过程的宽频带资料分析

１９８８年至１９９０年,青海地区发生了４次震级大于６．０的地震,本文利用全球数字地震台网（ＧＤＳＮ）的宽频带波形资料,通过波形模拟,结合地质结造的背影资料对这几个地震的震源破裂和发震构造背景进行了研究。通过台站的视震源时间函数（ａＳＴＦ）和视时间差（ａＴＤ）的分析,特别对震源的复杂性进行了讨论。本文的研究从震源分析的角度进一步支持了青藏高原东北部构造应力的压力轴为接的水平的ＮＥ方向,而随着向青藏高     

震源机制与烈度分布的关系

本文综合了哈斯克尔（N.A.Haskell）的有限长度的运动断层模型和安艺敬一（K.Aki）的所谓ω²模型（他的模型用了一个指数衰减函数作为位错速度的自相关函数）,略去了一些次要的因素和项,考虑了与烈度的联系,建立了一个计算烈度分布的模型。用此模型和震源机制工作给出的震源参数,计算了四个典型地震的理论烈度分布。同实际观测的烈度分布比较,两者基本一致。说明我们建立的模型,不仅可以解释许多实际的烈度分布（而有些是从地质构造角度无法解释的）,而且可以用于计算大面积的烈度分布。     

宁夏及邻区地震震源衰减模型及震源参数研究

选取2003年1月至2011年12月宁夏及邻区发生的201个M_L2.0以上地震事件进行反演, 用经验格林函数计算得到本区地震震源谱衰减模型。 在已知非弹性衰减系数和场地响应的条件下, 计算得到震源谱并研究了近震震级与拐角频率、 地震矩、 震源破裂半径及应力降之间的关系。 详细分析了震源参数误差大小及各种误差影响因素。 对剩余应力降时空演化特征分析研究结果表明: 高剩余应力降集中分布在宁夏北部宁蒙交界地区、 低剩余应力降分布在宁夏中部同心、 中卫一带, 银川盆地中南部的吴灵地区和宁夏南部海原断裂带中东段剩余应力降变化较为明显, 在时间分布上这两个区域剩余应力降在中强地震前均有降低—增强的趋势。     

用P波、S波初动和振幅比计算中小地震震源机制在云南地区的应用

利用云南数字地震台网记录的区域波形资料,用P、S波初动和P、S波的振幅比联合计算了发生在云南省内的91次中小地震的震源机制.通过与矩心矩张量法和P波初动法的结果对比分析,认为用P、S波初动和振幅比求解震源机制的方法是可行可靠的,用这种方法不仅可以得到5级以上地震的震源机制,也可以计算台网分布较好地区的3～4级地震的震源机制.另外,用91个地震的震源机制资料得到了云南地区震源机制类型、主压应力轴的平面分布及其统计图,结果表明云南地区以走滑型地震为主,主压应力轴优势分布方向为NNE向,这些结果与利用哈佛得到的结果是一致的.     

首都圈地区中小地震全波形震源机制反演

目前对首都圈地区中小地震的震源机制研究所用的方法主要有P波初动极性法、 振幅比法及波形反演方法。 前两种方法只利用了波形中很少量的信息, 且对于稀疏台阵效果较差, 而传统的波形反演法也仅利用了波形的部分信息, 且对于震级、 数据资料等有一定的要求。 为克服上述传统方法中的局限并进一步提高震源机制解的准确性, 本文采用一种新的综合利用上述各种信息的全波形匹配方法来反演首都圈地区中小地震的震源机制解, 提高反演解的可靠性。 为了检验所使用的全波形匹配震源机制反演方法的可靠性, 我们首先进行了合成数据测试, 结果表明全波形匹配方法可以得到稳定可靠的震源机制解。 利用新的全波形匹配方法计算了2019年4月北京发生的两次地震震源机制解, 并结合活动断裂进行分析。 将新的全波形震源机制反演方法应用于首都圈地区2015年以来的22个地震(2.2≤M≤4.3), 并对22个地震的震源机制解及应力轴分布进行了分析, 与前人研究结果具有较好的一致性。     

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

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