本文使用的软件是内蒙古地震局刘芳老师等编写的评估软件。根据软件说明书,前提条件和准备工作如下:①要求在研究时段内各台站地震震级的定义、台网对地震的触发条件保持不变;②在震级M(本文中震级M是指地方震震级ML)和震中距L的二维分布图布上标注每个台站记录到的地震情况,并作为计算台站检测概率的原始数据;③基于区域台网震级的定义,获得震级M和震中距L的经验换算关系。基于区域地震震级定义和实际产出的震相观测报告,计算每个台站对全部地震在时空上的检测能力[11]。计算原理如下:
地方震震级测定随震中距衰减关系一般表示为:
M=c1logA-c2logL*+c3(1)
式(1)中A为仪器记录的振幅;L*为震中距,c1 、c2和c3为常数。对于某个台站记录到两个地震的振幅均为A的情况下,测得的与两个时间震中距L1和L2相对应的的震级分别为:
M1=c1logA-c2logL1+c3(2)
M2=c1logA-c2logL2+c3(3)
在震级-距离二维图上表示的台站监测概率需要构建M与L的经验转换关系。将公式(2)与(3)相减可得:
ΔM*=│c2logL1-c2logL2│ (4)式(4)中震级差ΔM*仅与震中距L有关。根据规定[12],我国地方震震级(ML)的测定采用仿真短周期位移记录(DD-1或伍德-安德森)的S波或Lg波最大振幅来测定:
ML=lg ( A)+ R( L)(5)
式(5)中R(L)为台站的量规函数。由此,公式(4)中的ΔM*实际上仅与相应震中距L的量规函数R (L)有关。在计算震级-距离二维图上位置(M,L)对应的台站检测概率时,需选定计算所用数据。Schorlemmer和Woessner[11]定义了如下选取原则:
对台站周围发生的某次震级为M′、震中距为L′的地震事件,计算与位置(M,L)对应的震级差ΔM=M′-M ,以及利用公式(4)和(5)计算因不同距离引起的震级差ΔM*=R(L′)-R(L)。数据遴选的度量条件采用如下形式[11,13]:
当符合上述条件的地震事件数Nt≥10时,统计被台站检测到的地震事件数N+和未被检测的地震事件数N-,并计算台站在(M,L)处的检测概率:
“合成检测概率”PE(M,x,t)定义为在位置x、时刻t和震级M的地震能被台网中4个以上台站记录到的概率。综合各台结果,形成时刻t、震级M、位置x时的“合成检测概率”PE(M,x, t )和完整性震级M (P x ,t ) 。其相关计算原理如下:
式(8)中, P0E 、 P1E 、 P2E和P3E分别表示有0个、1个、2个和3个台站记录到的概率。台站没有检测能力的概率为
使用组合式的方法,得到能够被j个台站记录到的概率为
依据式(8)~(10),得到“合成检测概率”PE (M,x,t)为
合成得到基于概率的完整性震级,各震级档下完整性震级为:
其中:取Q=0.0001(误差标准)