强震复发间隔变异系数的一种计算方法

确定强震复发间隔变异系数在强震发生概率的计算中具有重要意义. 考虑到由小样本量地震序列计算得到变异系数的估计值与实际值可能存在一定的偏差, 本文提出了一种计算强震复发间隔变异系数的新方法. 该方法首先利用最大似然估计法计算出大量小样本地震序列的变异系数估计值; 然后先对这些变异系数的估计值进行标准化处理, 再对这些标准化的变异系数估计值进行统计分析; 最后在广泛搜集的39个地震序列的基础上, 利用该方法计算得到一个通用的变异系数值(α=0.34), 并讨论了该值的物理意义及其应用.      

垂直摆倾斜仪地震效应与倾斜场模拟

在整理垂直摆地震响应的基础上,采用USGS提供的有限断层模型和速度结构模型,运用德国地学研究中心（GFZ）王荣江研究员开发的（PSGRN/PSCMP）计算程序,分析日本M_W 9.0地震和中俄交界M_S 6.6地震倾斜场变化,并与垂直摆记录数据进行对比。通过整理现有地质资料和断层数据,模拟依舒断裂带北段鹤岗周边发生M_S 6.0地震时,该地区同震及震后倾斜场的变化,以了解不同震源性质的地震产生的地倾斜变化及对鹤岗地区产生的影响。     

2022年1月2日宁蒗MS5.5地震序列重定位与发震构造分析

2022年1月2日云南省丽江市宁蒗县发生M_S5.5地震。采用地震编目系统提供的此次地震的震相到时数据,使用双差定位方法对此次宁蒗M_S5.5地震序列进行重新定位,获得了694次地震的高精度相对位置。重定位后的地震空间分布显示:此次地震序列呈NNE至近NS向分布,与震源机制解的节面Ⅱ走向（191°）一致,主震位于地震序列南段;地震序列主体活动区长约11 km,宽约6 km,余震主要分布在4—11 km的深度范围内;地震序列在深度剖面上呈现出两组倾向不同的活动分支,其中东侧分支与震源机制解节面Ⅱ的倾角（81°）一致。此外,本次地震还可能触发了邻区的局部断裂活动。综合分析认为,2022年宁蒗M_S5.5地震的发震构造应该是NNE至近NS向兼具正断层分量的左旋走滑断层,倾向为WNW,倾角约为81°,其活动性质与震源区已知的活动断层均不一致。尽管本次宁蒗M_S5.5地震序列发生在2012年宁蒗—盐源M_S5.7地震序列的北侧,但是两次地震序列的发震断层并不相同。库仑应力反演结果显示,2012年宁蒗—盐源M_S5.7地震对本次宁蒗M_S5.5地震的发生具有促进作用。      

利用GRACE重力和重力梯度变化估计2012年苏门答腊地震断层参数

2012年4月11日苏门答腊北部附近海域发生M_W=8.6地震,国际上采用不同数据获得了该地震的断层滑动模型,但断层走滑性质存在较大差异.卫星重力GRACE（Gravity Recovery and Climate Experiment）观测覆盖震中区域,可以提供很好的断层参数估计约束.本文采用GRACE月重力场模型数据提取了此次地震同震重力和重力梯度变化,表明了北向分量的优越性.并与断层模型理论计算结果进行了比较,分析了T_xx（北-北向重力梯度）分量对断层参数的敏感性,以及基于粒子群算法以及Okubo位错模型反演了该地震断层参数.结果表明GRACE观测到的同震重力梯度变化在空间形态分布上与断层模型模拟结果符合较好,但振幅差异较大.重力和重力梯度的北向分量可以很好地压制条带误差,其中g_N（北向重力变化）和T_xx的误差最小,其次是T_xz（北-上向重力梯度）,误差最大和对条带误差抑制效果最差的是T_xy（北-东向重力梯度）.T_xx对断层的深度不敏感,对其余断层参数较敏感.基于GRACE反演得到的断层的走向角为113.63°,倾角为89.99°,滑移角为175.26°,平均滑移量为28.18 m,相应的矩震为8.71级,且此次地震的走滑性质为右旋走滑.     

新疆于田地区2008年以来3个地震 序列的参数早期特征

考察新疆于田地区2008年以来发生的2008年3月21日M_S7.3、 2012年8月12日M_S6.3和2014年2月12日M_S7.3等3个地震序列的参数早期特征, 利用“时空传染型余震序列”(ETAS)模型进行了参数估计. 统一选用截止震级M_c=M_L3.0, 并使用最大似然法分别估算了这3个序列不同持续时间的α、 p和b等参数. 研究结果表明, 3个地震序列的参数差异明显, 其中: 2014年M_S7.3地震序列α值与2012年M_S6.3地震序列接近, 但高于2008年M_S7.3地震序列; 2014年M_S7.3地震序列p值较低, 表明当前序列衰减缓慢; 2012年M_S6.3地震序列b值最大, 2014年M_S7.3地震序列次之, 2008年M_S7.3地震序列最小. 此外, 2008年M_S7.3地震序列参数在震后20天内变化较为显著, 而2014年M_S7.3地震序列的α值和p值在震后3.5天内则相对稳定.      

华北地区地震序列参数的分布特征

为系统地考察华北地区地震序列参数的分布特征,以期构建适合区域地震活动特征的短期概率预测模型和评估地震危险性,本文利用当前较为前沿的ETAS模型和R-J模型,采用连续滑动、多时段拟合的方式,对华北地区1970年以来的16个地震序列进行了参数拟合,并对参数的整体情况、参数之间的关系、参数与大地热流之间的关系等特征进行了分析研究.研究结果表明,主震发生后的早期阶段,地震序列参数变化相对较为剧烈,误差也较大;两种地震预测模型的序列参数呈现出一种优势分布特征,主要参数的平均值分别为α_ETAS=1.7404±0.3420,p_ETAS=0.9769±0.1396,a_OML=-1.6638±0.5284,b_OML=0.8312±0.1658,p_OML=0.9053±0.1527,这与国际上其他区域的研究具有很强的一致性;大多数情况下,ETAS模型参数p_ETAS高于R-J模型参数p_OML,平均偏高0.0716;序列参数α_ETAS与大地热流值整体上呈现一种负相关关系,激发次级余震的能力与大地热流具有一定的相关性.地震序列参数的分布特征以及与大地热流之间的关系对地震序列类型的判断,及基于地震序列参数构建定量化的短期概率预测模型,均有重要的应用价值.     

汶川地震破裂带断层岩纵波速度与孔隙度关系的实验研究

分布于地震破裂带上的断层岩具有高孔隙度的特征.该特点造成了其弹性波速度与结晶岩石和沉积岩存在明显的差异.确定断层岩的弹性波速度与孔隙度和矿物组成的关系对于利用地震资料探测深部断层和测井资料的解释至关重要.在10~600 MPa条件下,本文对地震断层岩的纵波波速（V_p）和总孔隙度（φ_t）进行了测量,并深入分析了V_p与孔隙度的关系.结果表明在10~600 MPa的压力范围内,V_p（p）随着压力的增高呈现对数增加,其增长率随着压力的上升而逐渐减小,遵从∂V_p（p）/∂p=a_v/p的变化规律.断层岩中的孔隙度随着压力的增高呈对数减小.与传统的认识不同,实验发现在压力高达600 MPa以上,大多数断层岩中仍然可以残留可观的孔隙量.分析显示V_p与总孔隙度及总粘土含量呈负线性相关.该发现有助于认识深部流体的活动通道特征,有助于理解断层带中存在大量粘土矿物、断层带内的物质可被大量带出、围陷波的形成等地质和地球物理现象.     

2013年芦山MS7.0地震序列参数的早期 特征: 传染型余震序列模型计算结果

为考察2013年4月20日芦山M_S7.0地震震后序列参数的早期特征, 利用“传染型余震序列”(ETAS)模型和最大似然法进行了参数估计. 设定截止震级M_c=M_L2.0, 拟合时段为震后0.31—24.12天, 计算获得α=1.89, p=1.22, 同时利用最大似然法估计获得b=0.72. 与中国大陆地区其它中强震的余震序列参数的比较表明, 芦山M_S7.0地震序列参数表现为触发次级余震的能力较弱和序列衰减速率较快的特征, 反映出余震区相对较高的应力水平. 为检测结果的稳定性, 设定不同的截止震级M_c以及不同的拟合截止时间, 分别进行参数拟合和参数标准差估计. 结果表明, M_c的选取对α值影响明显, 对p值影响则较小. 此外, 震后10天内获得的参数拟合结果随时间变化较为明显, 而其后各参数变化总体较为平稳.      

2017年漾濞MS4.8和MS5.1地震序列的微震检测及重定位

基于云南地震台网记录到的2017年漾濞M_S4.8和M_S5.1地震序列的地震目录和波形数据,采用结合波形互相关的双差定位方法得到了M_S5.1主震前7天至主震后30天共258次地震事件更加准确的震源位置,然后利用匹配定位方法获得了更加丰富的地震目录,并对其余震活动的时空演化特征进行了初步分析。双差重定位获得了漾濞地震序列中193次地震的精确位置,水平方向和垂直方向的定位误差均值分别为1.3 km和1.9 km,重定位及震源机制解反演结果显示维西—乔后断裂中南段应为此次地震的发震断层。使用匹配定位方法检测得到了1 625次中小地震事件,事件数量是云南地震台网目录中的6倍,被检测事件主要为主震后早期余震及M_L≤0微震事件。对比2016年云龙M_S5.0地震序列的研究结果后初步分析认为,“震群（多震）型”和“主余型”地震序列的一个主要区别可能是“震群（多震）型”地震序列的余震数量相对更少、余震活动衰减得更快。      

2021年6月10日云南双柏MS 5.1和6月12日云南盈江MS 5.0地震总结

2021年6月10日和6月12日,云南地区先后发生双柏M_S 5.1和盈江M_S 5.0地震,距5月21日云南漾濞M_S 6.4地震发生不足1个月,文中系统总结2次地震前出现的地震活动和地球物理观测异常,其中：地震活动方面：震前存在中等调制地震集中、地震条带等中短期异常;地球物理观测方面：2次地震发生前,震中附近存在形变、电磁和流体等多项异常,其中双柏M_S 5.1地震以形变学科趋势异常为主,盈江M_S 5.0地震以流体异常为主。对2次地震序列跟踪分析,发现双柏M_S 5.1地震余震较为丰富,震后半年内共记录M_L 3.0以上地震7次,最大震级为M_L 4.8,计算得到地震序列活动参数为：h值=1.23,b值=0.85;盈江M_S 5.0地震发生后半年内仅记录M_L 3.0以上余震4次,最大余震为M_L 3.3地震,计算得到地震序列活动参数为：h值=2.45、b值=1....     

論用PKP波測定极远地震的震級

一、引言目前我国测定极远地震（震中距离Δ>110°）震级的方法有两种: 1.面波法计算震级的公式 M_LH=log（A/T）-log（A/T）^*,（1.1） 式中A是水平分向地动位移（μ）,T是相应的周期（秒）,取面波中A/T为最大者,-log（A/T）^*是索洛维约夫-谢巴林1957年发表的量规函数。     

安徽"霍山窗"中小地震的 地震波衰减特征研究

利用安徽合肥台、 六安台和豹子崖台记录到的2001—2010年安徽“霍山窗”发生的M_L≥2.0地震的波形资料, 计算了其振幅衰减系数K值和卓越频率衰减系数a值, 分析研究了K值和a值的时空变化特征, 并讨论了其与华东地区中强地震的对应关系及与Q_c值的相关性. 结果表明, 3个台站均记录到了明显的变化特征: ① “霍山窗”M_L4.0地震、 江西九江M_S5.7地震及河南太康M_L5.0地震前6个月左右时间内, K值和a值均呈现持续下降的异常低值状态, 平均下降幅度K值为20%—35%, a值为45%—51%, 并于上升的过程中发震, 震后趋于正常; ② K值的平面分布特征显示, K值异常地震主要位于断裂带交汇区附近, 并沿断裂带方向分布; ③ K值、 a值与Q_c值基本呈负相关. 上述特征为进一步了解该区域的地壳介质状态和研究中强震前后Q_c的变化特征提供了一定的依据. 该研究结果可能对中长期地震预报和震后趋势判断有一定的指示意义.      

克拉通岩石圈对流减薄的数值模拟

采用二维有限元数值模拟的方法研究了岩石圈的对流减薄过程,特别是克拉通岩石圈的对流减薄过程.模型的主要参数包括增厚岩石圈的宽度x、增厚倍数γ、以及与岩石圈组分变化导致的黏性和密度变化密切相关的黏性比(ηc)和浮力数(B)或等效密度变化(Δρtc).数值计算结果显示,地幔对流将逐渐减薄增厚的岩石圈部分,(1)当B=0和ηc=1时,即对一般地幔岩石圈,增厚岩石圈对流减薄的时间可表示为0.0073γ0.70 x0.26.将数值结果应用于地球,意味着增厚到300km的岩石圈,如宽度为300km,对流移除增厚部分回到初始平衡厚度120km大约需要225 Ma;如宽度为1500km,移除增厚部分大约需要342 Ma.(2)当B和ηc较小,克拉通岩石圈对流减薄过程与一般加厚岩石圈的对流减薄过程类似,但减薄时间受克拉通组分浮力和黏性比的影响而显著增长,克拉通岩石圈对流减薄的时间可表示为0.0057ηc0.52Δρ-0.21tcγ0.78ηc-0.36 x0.04.因而,对300km厚的克拉通岩石圈,如克拉通岩石圈的密度比周围地幔的密度低0.4%(即B=0.1),宽度1500km,若克拉通岩石圈黏性因组分影响比普通地幔岩石圈大10倍,其被对流减薄到120km大约需要1.18Ga.(3)当B和ηc增大到一定量时(如B≥0.2且ηc10),克拉通岩石圈被移除的过程将发生变化,由于组分浮力的影响,对流主要不是将克拉通岩石圈带到软流圈地幔中,而主要是将较厚的岩石圈物质向两边推送.在此情况下,克拉通岩石圈能长时间(3Ga)保持稳定.     

亚热带富铁土的磁学性质及其磁性矿物学

通过５５个亚热带富铁土的磁化率（χ）、频率磁化率(χ_ｆｄ)、非滞后剩磁（ＡＲＭ）和饱 和等温剩磁（ＳＩＲＭ）等磁性参数测定,结合氧化铁化学形态分析和矿物鉴定,初步明确了亚热 带富铁土的磁学特征及其磁性矿物．磁测数据表明富铁土中存在强磁性的矿物,其磁化率χ 与土壤游离氧化铁（Ｆｅ_ｄ）含量呈极显著指数正相关（Ｒ^２＝０．５９７１）,频率磁化率χ_ｆｄ与土壤游离氧 化铁含量呈极显著直线正相关（Ｒ^２＝０．４２８９）．富铁土的χ_ｆｄ和非滞后磁化率χ_ＡＲＭ。值表明土壤 中的磁性矿物以超顺磁性（ＳＰ）和稳定单畴（ＳＳＤ）颗粒为主,富铁上的χ和χ_ＡＲＭ呈极显著直线 正相关（Ｒ^＝ ０．９４２９）,证明富铁土的磁性是由风化成土过程产生的 ＳＰ和 ＳＳＤ磁性颗粒贡献 的．矿物磁测结合Ｘ－衍射证明富铁土中的氧化铁矿物由赤铁矿、磁赤铁矿和针铁矿组成。     

三维物性分布的位场计算

本文在连续层状模型的基础上,将Parker位场公式推广到更具一般性的两种三维物性分布模型中。此时岩石的密度ρ和磁化强度J不仅在水平方向上可以任意变化,而且在深度上可依指数或线性函数形式变化。新的位场公式仍具有严密而简明的特点。结合圣安德烈斯断层带的重力反演,介绍了它的应用、模型分解和双界面同时反演的技术。     

用近震尾波估算昆明及其周围地区的Q值和地震矩

本文应用推广了的Aki尾波理论模式,忽略震源频谱含量的影响,考虑到仪器系统和地球Q-滤波器对尾波频散的影响,计算了f_p-t^*和k-t^*独立的两套尾波响应理论量板。在假定Q与频率有Q（f）=Q₀（f/f₀）ⁿ关系的基础上,用云南区域台网短周期地震仪器记录的小震尾波资料,估算了昆明周围地区频率约1.0Hz的短周期波的地壳Q₀值。结果表明,该区的Q₀值在150-220之间,平均为191;频率相关值η不是常数,多数为0.2-0.3。这些结果与用尾波估算川滇交界地区的结果接近,表明云南地区的地壳Q₀值与美国西部地区类似,是属于低Q₀值地区。 此外,用尾波与地震矩建立的关系,提出用尾波估算地震矩的一种可能方法,讨论了该区震级与地震矩的经验关系。     

各向异性水平岩层对自然电场电位曲线的影响

1) 引言现在自然电流法观察结果的解释,均假设矿体是处在均匀各向同性的介质中。实际上,围岩的不均匀和各向异性对电位曲线起着严重的崎变,因之,对结果的解释也应该作相应的校正。在这里,我们只研究最简单的均匀各向异性水平岩层对自然电场电位曲线的影响,井提出了用普通的方法从自然电位曲线上求矿体埋深时,应作的各向异性校正。     

Recurrence of Large Earthquakes: Bayesian Inference from Catalogs in the Presence of Magnitude Uncertainties

We present a Bayesian method that allows continuous updating the aperiodicity of the recurrence time distribution of large earthquakes based on a catalog with magnitudes above a completeness threshold. The approach uses a recently proposed renewal model for seismicity and allows the inclusion of magnitude uncertainties in a straightforward manner. Errors accounting for grouped magnitudes and random errors are studied and discussed. The results indicate that a stable and realistic value of the aperiodicity can be predicted in an early state of seismicity evolution, even though only a small number of large earthquakes has occurred to date. Furthermore, we demonstrate that magnitude uncertainties can drastically influence the results and can therefore not be neglected. We show how to correct for the bias caused by magnitude errors. For the region of Parkfield we find that the aperiodicity, or the coefficient of variation, is clearly higher than in studies which are solely based on the large earthquakes.