影响b值计算误差的Monte Carlo实验研究

本文假设地震的发生与时间的关系遵循稳态泊松模型,假设地震震级M服从参数为1/β的指数分布,其中β=b/loge,在此基础上生成地震目录.用最小二乘法估算b值,分别从样本数量、震级间隔、震级误差三方面分析b值误差的大小.统计分析结果表明,b值的标准偏差随着样本数量的增多而减小;随着震级间隔的增大而减小;震级误差使b值的标准误差变大.本文可以为利用真实地震目录求b值时选择适当的方法提供参考.     

内蒙古东部扎兰屯地区的b值时空变化特征

内蒙古东部的扎兰屯地区位于大兴安岭地震带北段,区域地质构造复杂,历史及现代中强地震较为活跃。本文以该地区地震活动作为研究对象,在区域地震序列完整性分析基础上,用最小二乘法进行b值时间扫描计算,用最大似然法进行b值空间扫描计算。时间扫描中的b值为每个扫描窗口内研究区的平均b值,因此其变化幅度不大,b值变化范围基本维系在0.78～1.13,误差范围为0.04～0.065。空间扫描结果显示,扎兰屯地区b值空间分布范围基本维系在0.4～1.8,误差范围为0.045～0.085;低b值区域为阿伦河断裂中北段,b值大致分布在0.5～0.7,表明该区域地壳介质处于高水平应力积累状态,可能是未来发生中强以上地震的危险区域。     

2016年门源6.4级地震序列b值变化特征及周期谱分析

本文对门源区域1970年以来的地震目录及2016年1月21日门源6.4级地震的地震序列、余震分布、不同时段完整震级、b值、地震周期特征等方面进行分析探讨.分析认为该区域的5级以上中强震多发生于b值的下降时段,中等地震存在3～5年的周期特征,6级以上地震的周期约为20年;本次6.4级的地震序列表明,余震密集发生于震后的10余天内,这个时段余震频次达到了5个月以来地震总数的85%以上.主震发生后,随着余震的增多,a值、b值都有上升的迹象,这表明区域应力状态得到了一定缓解,当b值趋于稳定时,可能预示着地震序列将进入结束状态.运用该区域主震前的地震目录分析计算发现,本次6.4级地震就发生在a值高值区及b值低值区,显示空间b值对本次强震的发生地点具有一定的指示意义.     

余震b值的先验概率分布

从世界上M≥7.5浅源地震伴随的81个余震序列,以及日本及其周围大地震的67个余震序列,确定了古登堡—里克特定律中b值的先验分布。我们假设先验统计总体中的b值服从伽马分布Γ(φ,ζ),其中φ和ζ分别代表形状和标度参数。对每个序列所得到的b值都与统计总体中的随机采样一致,但有统计误差。在本文中我们用最大似然法估计总体分布中的φ和平均值,b_(avg)=φζ。对世界上M≥5.0的余震序列估计的φ和b_(avg)分别是55和1.13,而对D<2.55的余震序列的估计值分别为28和1.22,其中D为主震和余震的震级差。对于日本的余震序列我们得到φ=28,b_(avg)=0.97。蒙特卡罗法的数值检验表明,所估计的b_(avg)相当接近于输入值,而φ却往往比输入值大。那些较小的值,如全球序列的φ=26～31,日本序列φ=20～25,或许比由最大似然法得到的估计值更合适些。假设b值的先验分布是Γ(φ,b_(avg)/φ),对一组N个地震,由最大似然法得到的估计值b=(N+φ-1)/(N/b_U+φ/b_(avg)),其中b_U是没用宇津所提出的先验分布而用最大似然法得到的估计值。所估计的值b比b_U更稳定可靠。特别是对于只有很少地震数据时,这里提出的方法很有价值。     

b值统计中的影响因素及危险性分析中b值的选取a

b值是从实际地震资料的统计中得到的.它与实际资料的完整性、可靠性、统计样本量的大小、取样的时空范围、样本的起始震级和取样间隔等因素都有关系. 本文使用累积频度-震级关系,lgN=a-bM,以起始震级为4 3/4和5;震级间隔为1/4,1/2和1级对不同尺度的时空范围进行统计,从中寻找各因素对b值统计结果的影响,为如何合理统计b值取得了一些实验依据. 此外,本文还就历史地震资料和现代仪器记录资料的统一使用问题,lgN-M线性关系两端掉头对b值统计的影响问题也作了探讨,从而提出在地震危险性分析中b值选取的原则及合理选取b值的方法.最后讨论了b值统计中的误差范围及其对危险性分析所产生的影响.      

地震目录的完整性对b值计算的影响

本文通过对华北地区和西南地区完整地震目录按照均匀抽取和非均匀抽取2种原则进行抽取地震,得到了地震目录遗漏率和b值变化量间的经验关系式.结果表明,在进行非均匀抽取地震时,得到的地震目录遗漏率和b值的减小量间呈正比关系;在进行均匀抽取地震时,利用最小二乘法和极大似然法计算得到的b值变化较复杂,b值随地震目录遗漏率的增加,或为增加,或为减少.同时,在进行非均匀抽取地震时,地震目录遗漏率对极大似然法计算的b值影响大于最小二乘法计算的b值影响.     

基于数据驱动的时间序列b值计算新方法(TbDD):以2021年云南漾濞Ms6.4地震序列为例

时间序列的b值在天然地震和工业开采诱发地震的危险性分析中具有重要的应用潜力,但长期以来受到计算规则设置的人为主观性、计算结果的可靠性和时序对突变识别精度不高等问题影响,制约了不同结果的可比较性和共识性科学认识的提炼.本文借鉴基于数据驱动(data-driven)的地震活动参数计算思路,采用连续函数形式的OK1993模型、时间轴随机段落划分、贝叶斯信息准则模型选择等技术环节,构建了基于数据驱动的时间序列b值计算新方法TbDD.利用合成地震目录的理论测试,并分别与固定地震数目的窗长和步长、固定地震数目的步长和累积窗长等传统的固定窗口法进行了比较研究.结果表明,TbDD方法可较好地还原合成地震目录的b.值输入参数,在计算规则设置的客观性和对b值突变过程的准确识别上具有明显优势.此外,我们还对新近发生的2021年5月21日云南漾濞Ms6.4地震序列进行了实际案例应用.结果显示,此次序列的b值在Ms6.4主震前为0.7左右、震前20 h出现了约0.1幅度的下降,表明在序列发生前震区的差应力水平较高.而b值在Ms6.4主震发生后起伏明显、逐渐增加至0.8左右,这一现象可能与震区在主震后早期较为剧烈的应力调整有关.进一步针对随机模型的数量以及时间轴的随机段落划分设置对TbDD方法b值计算结果的影响程度进行了测试,发现b值受随机模型数量影响较小、具备较强的稳定性,时间轴的随机段落划分设置可影响b值时序微观起伏变化的识别.本文发展的TbDD方法在对时间序列b值计算的准确性、余震趋势跟踪的高精度要求,以及工业开采诱发地震风险管控等领域有较好的应用潜力,所获得的2021年云南漾濞Ms6.4地震序列的b值计算结果也对理解此次地震序列的孕育过程有参考价值.     

b值误差范围的研究

从历史地震的角度,对震级误差所引起的b值变化进行了讨论。通过推导,给出了可能产生的误差范围,从而对导致b值变化的原因进行了探讨。     

用B值横截距预报强余震震级的方法探讨

本文在G-R(logN=α—M)关系的基础上,探讨用b值横截距法预报地震序列强余震的震级。选用我国1966—1986年21年内有完整资料的58个Ms≥5.0地震序列,用b值横截距法进行预报检验。这些地震序列的位置分布较广,遍及全国各大区,而且要求每个地震序列满足震级M≥Mo(起始震级)的地震个数不少于50个。统计结果表明,对预报序列中强余震的吻合率为93％以上。另外,我们还用b值横截距法对序列的晚期强余震预测及在主震后尽可能短的时间内预测强余震进行了探索,表明它是可能的。     

汶川地震序列b值的分析研究

采用Aki等（1965）和Shi等（1982）计算b值及标准差的方法,分析了2008年5月12日～2008年8月31日四川地震台网所记录到的汶川地震序列ML≥2．0地震的b值,在空间、时间以及随震源深度的变化特征。结果显示：汶川地震序列b值在震源区上分布是非均匀的,但基本在[0．5,0．9]之间变化,震源区的中段b值高于两端;同时发现,几次M≥6．0余震都发生在b值较低的区域内;震源区b值具有随震源深度增加而减小的趋势,在20km左右处,减小的趋势最为明显;随着时间的推移,b值波动变化幅度逐渐减小,周期变大。     

兰州地震研究所发展的三段历史回顾

对兰州地震研究所的三段历史作回顾,分别为国家地震局、中国地震局兰州地震研究所的体制沿革变迁,中国科学院地球物理研究所;中国科学院兰州地球物理研究所观象台早期创建;兰州地震研究所硕士点早期阶段。     

山西大同数字遥测地震台网记录的震相特征

从地球及地壳的分层结构出发，使用较成熟的理论，从地震波传播的基本原理入手，结合多年的震相观察分析工作实践，选取了山西大同数字遥测地震台网近年来的一些典型震相进行了分析，得出本台网记录到震中距小于100km的地方震震相有Pg，P11，Sg，S11震中距在100km～500km的近震震相有Pn，Pg，P11，Sn，Sg,S11对比分析了在本台网内及周边地区发生的爆破与天然地震在震相上的差异性，并对其特征作了较详细的描述。     

淋巴结结核的影像学表现特点分析

目的:分析淋巴结结核的CT和MRI表现特点,并与病理学对照,探讨两种检查方法的应用价值。方法:回顾性分析2010年1月至2014年12月间34例经手术、穿刺病理或临床证实为淋巴结结核患者的CT或MRI影像学资料,分别从生长方式（融合/非融合）、边缘（光滑/欠光滑）、境界（清晰/欠清晰）、密度及信号、增强后强化方式（均匀/不均匀/环形规则/环形不规则）这5个方面对所受累淋巴结进行分析。结果:34例患者中共累及223枚淋巴结,其中颈部淋巴结67.7%,纵隔内淋巴结27.7%,肺门淋巴1.8%,肺内淋巴结0.5%,腹腔淋巴结0.5%,腹膜后淋巴结1.8%。所有受累淋巴结中,8.5%呈融合性生长,83.9%形态规则,79.4%边缘光滑。CT上发现的171枚病变淋巴结以等密度较多（53.8%）、低密度次之（43.3%）,MRI上发现的66枚病变淋巴结中98.5%为等T1长T2信号。所有CT及MRI观察到的病变淋巴结中,除57%的病灶密度或信号均匀外,仅伴钙化的为12.6%,仅伴部分坏死的26.9%,另外有3.6%同时伴钙化及坏死;增强扫描的204个病变淋巴结中,1%钙化显著、未见明显强化,38.7%均匀强化,伴坏死者内见不强化的低密度区,11.3%呈环形规则强化（均位于颈部）,8.3%环形不规则强化（纵隔内82.4%）,5.9%可见病变淋巴结内有分隔样强化（均位于颈部）。结论:CT与MRI对于淋巴结结核具有重要诊断价值,两者均能清晰显示病变部位、形态以及边缘情况,MRI对于病变边缘及周围情况显示较CT更明确,尤其是增强后,对于包膜及邻近组织受累情况的显示也更加敏感。      

Correction of twilight measurements of influence of night sky

Studia Geophysica et Geodaetica -     

山西夏县中心地震台温泉水化学特征

采集夏县中心地震台温泉水及周边水点样品进行水化学组分特征分析。利用矩形图、Na-K-Mg三角图等方法,初步分析温泉水的水质类型、水-岩平衡状态、热储温度以及循环深度等,并结合氢氧同位素组成特征,初步分析温泉水补给来源。研究表明,夏县中心地震台温泉水化类型属于Na-Cl·SO₄型,水-岩反应属于部分成熟水,热储温度为148.8℃,循环深度为3.91 km,补给源主要为大气降水,温泉为断裂型温泉。以上结果可为该台流体异常分析提供基础研究资料,为夏县地区水化研究奠定基础。     

山东数字地震台网子台背景噪声位移值的分布特征

利用山东数字地震台网的实时记录波形资料,采用在频率域积分方法,计算出各台站背景噪声的位移值。然后进行统计分析,得出了山东地区背景噪声位移值的空间分布和时间变化特征：台站背景噪声幅值与不同地区地质环境、岩性密切相关,位于第四系黄土覆盖层较厚的鲁西聊城、菏泽地区,鲁西北德州地区和鲁北东营、滨州地区各子台的背景噪声位移值较大,鲁中、鲁南山区和胶东半岛地区的基岩台站背景噪声位移值较小;同时,白天噪声值大于夜间噪声值,但深井观测台站的背景噪声值相对稳定。此项研究为绘制台网监测能力图提供基础数据,并为计算各台站的场地响应及震级修正值等测震学研究提供基础资料,促进台网的改建和发展。