[发明专利]一种地震台网定位能力评估方法

权利要求书

1.一种地震台网定位能力评估方法，其特征在于，包括以下步骤：

步骤1、获取地震台网中每个地震台站对震源位于x处震级为M的地震P波信号检测概率，记为P，且P的取值介于0-1之间；定义任意y处的地震台站定位参与概率p_i等于震源位于x处震级为M的地震被y处地震台站检测到P波信号的概率P，即p_i＝P(y,x|M)，其中，i为1,2,3,……,N，N为地震台网中地震台站的总数；

步骤2、随机生成一个数值介于0-1之间，长度为N的随机数组ζ；其中，数组元素ζ_i表示第i个地震台站的当前定位概率，当p_i＞ζ_i时，第i个地震台站参与定位，依此确定参与当前地震定位的地震台站；

步骤3、判断参与当前地震定位的地震台站个数是否小于3；若否，则执行步骤4；若是，则表示此时地震台网对该地震事件没有定位能力，执行步骤6，并将当前定位覆盖椭圆面积设为无限大值；

步骤4、根据确定的参与当前地震定位的地震台站，计算震源的水平位置定位误差协方差矩阵的两个特征值λ₁和λ₂，要求λ₁≤λ₂；其中，矩阵元素B₂₂、B₂₃、B₃₂、B₃₃由定位时的走时偏导矩阵确定；

步骤5、根据下式计算当前定位覆盖椭圆的长半轴a和短半轴b：

其中，σ_T为走时模型误差的方差；σ_t为信号到时测量误差的方差；表示自由度为2的χ²分布的上α分位数；

步骤6、计算当前定位覆盖椭圆面积S＝πab，并保存结果；

步骤7、重复步骤2-6共K次，得到K个定位覆盖椭圆面积；

步骤8、将K个定位覆盖椭圆面积从小到大排序，取其中第q个定位覆盖椭圆面积为所述地震台网对震源位于x处震级为M的地震的定位能力；其中，q＝int(K×(1-α))，1-α为定位误差分布置信度；int()表示对括号内的计算结果向无穷小方向取整。

2.根据权利要求1所述的地震台网定位能力评估方法，其特征在于：步骤7中所述K不小于800。

3.根据权利要求1或2所述的地震台网定位能力评估方法，其特征在于：步骤8中所述1-α＝0.9，即定位覆盖置信度为90％。