随着“十五”测震软件JOPENS在全国的推广应用, 各区域地震台网的地震分析处理能力有了明显提高，JOPENS—MSDP地震分析软件提供了多种地震定位方法，实际工作中该如何根据地震台站分布的实际情况选择合适的定位算法，从而使定位结果的精度更高或者更可信，一直困扰着地震研究者。多个地震台网的技术人员曾对同一地震几种算法的不同定位结果进行对比研究，如：陈贵美等（2009）认为，网内地震各种定位方法的差异不大，效果较好；张炳等（2012）分析认为，对于网内地震，Hyp2000程序稍优于单纯型定位方法；对于网缘及网外地震，单纯型定位方法优于Hyp2000程序；管勇等（2017）认为，随着空隙角增大，单纯型定位方法结果较差，其余3种方法差别不大。如果有事件发生的精确时刻和地点，将其与地震台网人机交互结果进行对比，可以估计出发震时刻和地理位置的绝对误差，从而检验测震台网所用定位软件的定位能力。王平等（2006）和李桂华等（2006）分别使用陕西地震台网和云南地震台网资料检验过HYPO81、GEIGER-P、GENETIC和KM-n定位程序；吕作勇等（2014）和黄文辉等（2017）也曾讨论过广东地震台网的定位精度。以上大部分研究是利用地震事件进行定位方法的对比, 或不区分地震定位方法的定位精度准则统计，实际工作中可以使用抽样概率方法，更合理地研究各种地震定位方法的不同特点。

本文利用广东地震台网记录的2个已知起爆时刻和爆破位置的人工爆破事件，使用JOPENS—MSDP提供的Hyposat、Hyp2000、Locsat和单纯型定位算法，进行分析处理2次爆破资料，检验4种定位方法的结果偏差，获取一些经验性认识，供大家参考。

在科技部国家科技支撑计划项目资助下，中国地震局地震预测研究所、中国地震局地球物理勘探中心与广东省地震局合作，在新丰江水库开展主动源深地震探测工作，于2009年12月31日及2010年1月1日，分别在广东省陆丰市陆河镇和清远市英德大华镇进行人工地震爆破激发，2次爆破均采用三孔同时起爆方式，每孔间距10 m，孔深50 m，25 m以下放600 kg的TNT，25 m以上填充砂石。

爆破前在现场使用手持GPS测量爆破点的经度、纬度和高程，爆破时使用起爆器记录爆破时刻，爆破后查阅广东地震台网实时记录资料发现，2次爆破记录均比较理想，爆破波形被距爆破点150 km内的所有地震台记录到，但S波不发育，台站分布见图 1。英德爆破处于网内，对爆破波形有42个台站记录，英德台（YND）距离最近，震中距34 km，普宁台（PUN）距离最远，震中距271 km，可以识别27个Pg震相和22个Pn震相；陆河爆破处于网缘，对爆破波形有37个台站记录，揭西台（JIX）距离最近，震中距31 km，肇庆台（ZHQ）距离最远，震中距312 km，可以识别出31个Pg震相和21个Pn震相。

图 1(Fig. 1)

图 1 爆破位置与接收台站分布 (a)广东英德爆破；(b)广东陆河爆破 Fig.1 Explosion location and the stations distribution

在GPS同步时钟控制下，起爆器引爆电雷管的时间精度可以达0.001 s。地震台站端的GPS授时精度优于±1 μs，时间服务精度钟差不超过1 ms（中国地震局监测预报司，2003），广东地震台网在收集地震信号时采样频率为100 Hz，因此决定了使用JOPENS—MSDP软件分析每个台站记录时，到时误差不小于10 ms。总体来说，起爆时刻精度至少比记录系统时间精度高一个数量级，实测走时t_i可以用地震台站记录到时减去起爆时刻获得。

各模型下的理论走时T_i可依据震中距查询J—B、IASPEI91、AK135和华南走时表获取。震中距（相当于炮点距）计算公式为式（1），式中λ₁、φ₁、λ₂和φ₂分别为2点的纬度和经度。按照震相分类绘制各模型下的理论走时T_i与实测走时t_i走时曲线，见图 2。由图 2可见，各模型下的理论走时T_i与实测走时t_i之间存在明显差，可用式（2）计算RMS值定量表述。

图 2(Fig. 2)

图 2 2次爆破Pg和Pn震相实测走时及理论走时 (a) Pg走时；(b)Pn走时 Fig.2 Pg and Pn phase observational and theoretical for two explosion

计算得到各模型下RMS值为：J—B走时表为0.474 s，IASPEI91为0.393 s，AK135为0.333 s，华南走时表为0.234 s。RMS值越大，说明假设模型与真实模型偏差越大，显然，本区域使用华南走时表测定时模型误差较小，而使用J—B走时表测定时模型误差较大。

Hyposat、Hyp2000与Locsat均采用传统的Geiger法，具体计算方法相同，即观测方程组降维后求解，不化为正规方程组，直接用奇异值分解最小二乘法方程组，实际计算中采用多种数据加权。3种定位方法的主要区别在于，所参与计算的资料权重设置与初始值确定。单纯型定位法利用数学的单纯型算法搜索残差最小位置即为震中，该方法在极值附近收敛比较慢，对初始值较敏感。几种定位方法各有长短，本次测试初始深度统一设置为1 km，均使用JOPENS—MSDP默认的参数设置对2次爆破事件进行定位，英德爆破事件的定位结果见表 1，陆河爆破事件的定位结果见表 2。

表 1(Tab. 1

表 1 JOPENS—MSDP中4种定位程序测定的英德爆破位置（网内）Tab. 1 The location results of Yingde explosion by 4 methods in JOPENS-MSDP(Intra network) 定位方法 模型 t/s φN/(°) λE/(°) 深度/km 平均残差/s 爆破时刻误差/s 水平参考误差/km 垂直参考误差/km 与GPS实测水平误差/km 单纯型 IASPEI91 15.87 24.234 8 113.665 7 8.3 0.166 -0.101 - - 1.462 单纯型 华南 16.39 24.233 7 113.674 3 6.6 0.160 0.419 - - 1.214 Locsat IASPEI91 16.27 24.224 0 113.680 0 0 1.028 0.300 4.400 0 2.199 Locsat 华南 16.43 24.234 0 113.674 0 5.4 0.189 0.459 3.100 3.80 1.192 Hyp2000 IASPEI91 17.35 24.232 0 113.673 7 5.0 1.930 1.379 3.300 17.91 1.405 Hyp2000 华南 16.35 24.238 8 113.677 0 1.3 0.080 0.381 0.685 1.17 0.685 Hyposat AK135 15.67 24.232 0 113.674 0 5.9 0.196 -0.301 0.612 0.69 1.397 Hyposat IASPEI91 15.67 24.232 0 113.675 0 5.9 0.195 -0.301 0.612 0.69 1.374 Hyposat J—B 14.74 24.224 0 113.672 0 0.1 0.450 -1.231 1.256 2.25 2.285 Hyposat 华南 16.41 24.238 0 113.676 0 4.7 0.128 0.439 0.411 0.53 0.732 注：①时刻以00:00:00.000作为零时换成秒表示；②水平参考误差和垂直参考误差是对应定位方法所提供的误差。表 2与此相同。

表 1 JOPENS—MSDP中4种定位程序测定的英德爆破位置（网内） Tab.1 The location results of Yingde explosion by 4 methods in JOPENS-MSDP(Intra network)

表 2(Tab. 2

表 2 JOPENS—MSDP中4种定位程序测定的陆河爆破位置（网缘）Tab. 2 The location results of Luhe explosion by 4 methods in JOPENS-MSDP(margin network) 定位方法 模型 t/s φN/(°) λE/(°) 深度/km 平均残差/s 爆破时刻误差/s 水平参考误差/km 垂直参考误差/km 与GPS实测水平误差/km 单纯型 IASPEI91 15.00 23.148 2 115.614 0 6.1 0.169 -0.755 - - 3.488 单纯型 华南 15.90 23.165 3 115.595 2 7.0 0.258 0.145 - - 1.428 Locsat IASPEI91 16.34 23.152 0 115.613 0 18.5 0.352 0.585 4.600 4.200 3.075 Locsat 华南 15.86 23.156 0 115.603 0 5.9 0.326 0.105 4.100 3.900 2.515 Hyp2000 IASPEI91 17.12 23.160 8 115.632 8 5.0 1.960 1.365 3.950 16.37 2.714 Hyp2000 华南 17.12 23.160 8 115.632 8 5.0 1.960 1.365 3.950 16.37 2.714 Hyposat AK135 14.89 23.161 0 115.586 0 2.8 0.264 -0.865 0.956 0.91 1.863 Hyposat IASPEI91 14.90 23.161 0 115.586 0 2.8 0.264 -0.855 0.967 0.91 1.863 Hyposat J—B 14.04 23.205 0 115.552 0 0.6 0.513 -1.715 1.634 2.30 3.340 Hyposat 华南 16.04 23.175 0 115.587 0 4.3 0.148 0.285 0.534 0.57 0.306

表 2 JOPENS—MSDP中4种定位程序测定的陆河爆破位置（网缘） Tab.2 The location results of Luhe explosion by 4 methods in JOPENS-MSDP(margin network)

从表 1和表 2中的具体定位结果可以看出：无论使用区域走时表还是全球走时表，每种定位方法的定位结果与实际爆破相比，爆破时刻误差绝大多数在±1 s以内；水平误差未超过5 km，根据《地震编目规范》(DB/T 66—2016)，达到1类精度要求；垂直误差较大，最小0 km最大18.5 km；定位方法提供的参考误差，水平向与实测偏差在同一数量级上，垂直向偏差较大不具备太大的参考价值。网内的英德爆破水平误差在0.7—2.3 km，网缘的陆河爆破水平误差0.3—3.5 km。对于同种定位算法，使用华南走时表（模型）水平误差最小，使用J—B走时表误差最大，特别是使用Hyposat方法进行定位时，在区域走时模型下测定水平误差小于1 km，爆破时刻误差在±0.5 s以内。

为了进一步比较4种定位方法在测定不同地震台站布局时结果的稳定性，选取台站分布较好的英德爆破，采用抽取不同台站组合模拟不同台站布局的方法，抽取155万多种台站组合，即从42个可使用台站中分别抽取3个台的11 480种组合、4个台的111 930组合、5个台的85 067种组合（随机）、6个台的52 458种组合（随机）、37个台的85 067种组合（随机）、38个台的111 930种组合、39个台的11 480种组合、40个台的861种组合、41个台的42种组合和42个台的1种组合，6—36个台随机抽取一百多万种组合。选择华南走时模型，分别对每种台站组合使用4种方法进行爆破定位，定位结果及台站布局见图 3，可见：水平误差随着使用台站数增多台站布局变好而减小。

图 3(Fig. 3)

图 3 随机抽取台站组合情况及定位结果偏差 Fig.3 Random station combination and deviation of location results

最大空隙角是描述台站布局好坏的重要参数之一，简单的台站布局分类以最大空隙角划分：小于180°为网内，180°附近为网缘，大于180°为网外。与此类似，按照最大空隙角，以15°为一档对抽取的台站组合进行分类，统计各类别抽取台站组合样本数及各类别中4种方法的定位结果满足指定精度的次数，两者相除得到满足指定精度的抽样概率，具体统计结果见表 3。相同台站布局下，满足指定精度的抽样概率反映对应定位结果的可信程度；不同台站布局下，满足指定精度和置信度的台站布局档数反映定位方法的适用范围，即稳定程度。爆破时刻误差概率、水平误差概率及垂直误差概率随最大空隙角的变化分别见图 4、图 5和图 6。

表 3(Tab. 3

表 3 JOPENS—MSDP中4种定位方法的抽样概率Tab. 3 The sampling probability of four location methods in JOPENS-MSDP 最大空隙角/（°） 定位算法 发震时刻精度类别及概率(%) 水平误差精度类别及概率(%) 垂直误差精度类别及概率(%) 抽取样本数 ≤0.5 s ≤1 s ≤5 s ≤1.0 km ≤5 km ≤15 km ≤30 km ＞30 km ≤1 km ≤5 km 30—45 Hyposat 96.6 100.0 100.0 83.6 100.0 100.0 100.0 0.0 0.0 48.7 349 256 30—45 Hyp2000 92.9 100.0 100.0 83.9 100.0 100.0 100.0 0.0 15.0 99.5 349 256 30—45 Locsat 99.5 100.0 100.0 8.7 100.0 100.0 100.0 0.0 0.0 0.6 349 256 30—45 单纯型 2.5 100.0 100.0 6.0 100.0 100.0 100.0 0.0 0.0 4.2 349 256 45—60 Hyposat 76.7 100.0 100.0 54.2 100.0 100.0 100.0 0.0 0.0 33.9 259 749 45—60 Hyp2000 81.6 99.9 100.0 50.6 100.0 100.0 100.0 0.0 25.4 96.8 259 749 45—60 Locsat 92.3 100.0 100.0 17.0 100.0 100.0 100.0 0.0 0.0 2.5 259 749 45—60 单纯型 5.9 100.0 100.0 19.5 100.0 100.0 100.0 0.0 0.0 5.7 259 749 60—75 Hyposat 87.1 100.0 100.0 66.3 100.0 100.0 100.0 0.0 0.0 65.2 261 096 60—75 Hyp2000 82.4 99.9 100.0 53.5 99.9 100.0 100.0 0.0 25.2 97.1 261 096 60—75 Locsat 95.4 100.0 100.0 17.3 100.0 100.0 100.0 0.0 0.0 12.1 261 096 60—75 单纯型 24.6 100.0 100.0 14.8 100.0 100.0 100.0 0.0 0.0 7.7 261 096 75—90 Hyposat 78.8 100.0 100.0 57.1 100.0 100.0 100.0 0.0 0.0 59.4 161 782 75—90 Hyp2000 75.7 98.5 100.0 41.0 99.1 100.0 100.0 0.0 28.5 96.2 161 782 75—90 Locsat 90.0 100.0 100.0 20.3 100.0 100.0 100.0 0.0 0.0 13.8 161 782 75—90 单纯型 22.5 100.0 100.0 17.2 100.0 100.0 100.0 0.0 0.0 13.3 161 782 90—105 Hyposat 80.0 100.0 100.0 46.7 100.0 100.0 100.0 0.0 0.3 60.1 83 773 90—105 Hyp2000 74.0 96.7 100.0 43.3 96.8 99.8 100.0 0.0 26.5 91.9 83 773 90—105 Locsat 88.2 100.0 100.0 19.4 99.9 100.0 100.0 0.0 0.2 19.3 83 773 90—105 单纯型 33.5 99.9 99.9 18.8 99.9 99.9 99.9 0.1 0.1 19.3 83 773 105—120 Hyposat 74.7 100.0 100.0 46.7 100.0 100.0 100.0 0.0 0.6 61.0 76 651 105—120 Hyp2000 56.6 88.1 99.8 17.0 87.1 98.5 99.8 0.2 32.1 87.4 76 651 105—120 Locsat 82.7 99.9 99.9 20.1 99.8 99.9 99.9 0.1 0.3 21.9 76 651 105—120 单纯型 19.3 99.9 99.9 21.3 99.9 99.9 99.9 0.1 0.1 16.2 76 651 120—135 Hyposat 66.7 99.7 100.0 39.7 99.7 99.9 99.9 0.1 1.5 51.9 55 507 120—135 Hyp2000 49.2 84.6 99.4 17.6 82.6 97.4 99.4 0.6 22.8 77.2 55 507 120—135 Locsat 72.3 99.2 99.3 25.4 99.2 99.2 99.3 0.7 1.0 23.4 55 507 120—135 单纯型 31.3 99.5 99.5 15.7 99.3 99.5 99.5 0.5 0.1 16.5 55 507 135—150 Hyposat 70.6 99.6 100.0 33.1 99.2 99.8 99.9 0.1 1.7 51.9 57 800 135—150 Hyp2000 45.9 79.8 98.8 13.5 74.3 94.4 98.6 1.4 23.9 72.6 57 800 135—150 Locsat 76.2 98.7 98.8 22.4 98.3 98.8 98.8 1.2 1.4 25.2 57 800 135—150 单纯型 43.7 99.3 99.4 12.2 98.2 99.3 99.3 0.7 0.1 16.8 57 800 150—165 Hyposat 68.4 99.5 99.9 26.7 98.8 99.8 99.9 0.1 2.2 55.0 57 171 150—165 Hyp2000 42.2 75.0 98.1 5.8 63.8 90.9 97.5 2.5 24.3 73.4 57 171 150—165 Locsat 71.9 98.1 98.2 17.9 97.4 98.2 98.2 1.8 1.8 30.0 57 171 150—165 单纯型 35.2 99.1 99.4 12.9 95.8 99.4 99.4 0.6 0.1 16.5 57 171 165—180 Hyposat 64.9 99.2 99.9 23.0 97.4 99.7 99.9 0.1 3.5 58.3 41 816 165—180 Hyp2000 44.1 80.4 98.1 7.5 65.1 89.9 96.8 3.2 19.6 69.0 41 816 165—180 Locsat 65.7 96.5 96.8 17.9 94.6 96.5 96.6 3.4 3.1 38.1 41 816 165—180 单纯型 33.6 98.2 98.9 9.9 91.6 98.4 98.6 1.4 0.2 18.2 41 816 180—195 Hyposat 58.3 99.0 99.9 21.3 95.4 99.7 99.8 0.2 2.9 47.5 28 049 180—195 Hyp2000 40.3 72.9 96.7 5.8 55.2 84.4 93.9 6.1 18.3 65.7 28 049 180—195 Locsat 59.4 95.5 95.9 16.6 91.2 95.2 95.4 4.6 2.8 29.1 28 049 180—195 单纯型 32.8 97.9 98.9 8.1 85.6 97.8 98.4 1.6 0.2 19.2 28 049 195—210 Hyposat 58.6 99.2 99.9 17.0 96.1 99.6 99.8 0.2 2.2 51.5 31 898 195—210 Hyp2000 33.7 59.4 97.2 2.7 42.3 79.0 94.1 5.9 24.2 67.1 31 898 195—210 Locsat 62.1 94.8 95.5 13.2 91.3 94.5 94.7 5.3 3.2 30.5 31 898 195—210 单纯型 27.2 90.6 99.1 7.1 81.1 97.8 98.3 1.7 0.1 17.3 31 898 210—225 Hyposat 53.6 96.6 99.8 16.5 91.8 99.2 99.8 0.2 4.1 57.8 25 374 210—225 Hyp2000 33.1 62.1 96.7 4.6 46.7 79.2 92.4 7.6 18.3 65.5 25 374 210—225 Locsat 50.7 90.6 94.1 13.4 85.3 92.1 92.6 7.4 4.9 38.6 25 374 210—225 单纯型 26.3 80.0 98.0 6.7 71.3 93.9 94.9 5.1 0.1 17.9 25 374 225—240 Hyposat 54.9 94.2 99.8 11.0 86.1 99.0 99.7 0.3 5.8 62.4 17 204 225—240 Hyp2000 26.0 47.2 95.1 3.0 31.4 69.1 89.3 10.7 18.1 64.8 17 204 225—240 Locsat 53.1 87.1 93.5 9.1 78.0 90.1 90.9 9.1 6.8 44.0 17 204 225—240 单纯型 24.9 63.7 97.2 3.8 56.0 91.6 92.9 7.1 0.1 16.6 17 204 240—255 Hyposat 45.8 87.7 99.9 9.8 80.2 98.9 99.8 0.2 3.3 58.5 14 883 240—255 Hyp2000 15.7 32.6 89.3 1.0 21.6 58.8 82.1 17.9 18.6 62.5 14 883 240—255 Locsat 42.1 77.3 91.6 7.6 70.8 87.2 88.1 11.9 4.5 39.6 14 883 240—255 单纯型 10.6 31.8 91.3 2.4 32.9 83.3 85.9 14.1 0.2 15.8 14 883 255—270 Hyposat 44.7 81.3 99.7 6.1 70.4 97.1 99.6 0.4 7.9 60.5 10 873 255—270 Hyp2000 18.1 34.7 84.8 1.2 22.2 54.9 76.4 23.6 14.8 62.4 10 873 255—270 Locsat 39.6 70.9 88.4 5.5 59.5 84.4 86.3 13.7 7.6 44.0 10 873 255—270 单纯型 8.6 26.5 88.7 1.4 28.0 80.3 84.3 15.7 0.2 15.9 10 873 270—285 Hyposat 38.3 70.4 99.3 5.4 54.4 95.0 99.2 0.8 9.9 57.1 9 253 270—285 Hyp2000 9.9 19.5 66.1 0.8 11.5 36.7 57.9 42.1 18.3 63.7 9 253 270—285 Locsat 31.2 56.8 83.6 4.4 42.4 77.3 81.7 18.3 8.9 45.9 9 253 270—285 单纯型 4.6 14.0 81.4 0.9 15.9 70.7 79.2 20.8 0.2 15.3 9 253 285—300 Hyposat 25.4 51.4 98.7 2.9 38.0 85.4 98.4 1.6 14.5 59.3 3 837 285—300 Hyp2000 4.8 11.2 50.4 0.2 5.8 24.6 44.5 55.5 13.1 64.4 3 837 285—300 Locsat 17.8 34.7 77.5 2.2 25.8 61.4 76.5 23.5 18.0 54.6 3 837 285—300 单纯型 1.6 5.5 82.7 0.2 5.3 59.2 78.9 21.1 0.1 18.3 3 837 300—315 Hyposat 17.7 38.0 96.2 1.6 30.2 77.0 95.4 4.6 15.9 50.0 3 238 300—315 Hyp2000 2.8 5.2 30.0 0.2 2.6 12.2 25.7 74.3 9.9 58.4 3 238 300—315 Locsat 7.8 16.7 67.4 1.1 12.4 46.4 66.6 33.4 22.4 51.1 3 238 300—315 单纯型 1.2 3.3 78.5 0.1 3.3 43.4 74.1 25.9 0.1 15.2 3 238 315—330 Hyposat 20.7 36.7 88.2 1.9 26.0 69.5 87.6 12.4 15.6 48.0 1 394 315—330 Hyp2000 2.4 3.7 20.9 0.2 2.2 7.7 19.4 80.6 9.7 61.1 1 394 315—330 Locsat 4.1 10.8 63.6 0.3 5.4 38.4 62.0 38.0 33.9 62.2 1 394 315—330 单纯型 1.9 5.6 80.3 0.1 3.4 39.2 73.0 27.0 0.5 23.2 1 394 330—345 Hyposat 16.7 33.5 79.3 0.7 20.1 61.7 79.1 20.9 11.0 51.8 546 330—345 Hyp2000 0.0 0.0 4.9 0.0 0.0 1.5 4.2 95.8 7.1 63.0 546 330—345 Locsat 2.4 5.9 63.4 0.0 2.0 30.0 60.6 39.4 45.6 70.5 546 330—345 单纯型 1.1 6.8 78.9 0.0 3.8 39.4 69.0 31.0 0.2 17.8 546 345—360 Hyposat 28.3 39.1 71.7 0.0 11.6 58.7 68.8 31.2 10.1 61.6 138 345—360 Hyp2000 0.0 0.7 8.7 0.0 0.0 1.4 5.8 94.2 2.2 55.1 138 345—360 Locsat 0.7 5.8 42.0 0.0 0.0 6.5 34.8 65.2 36.2 68.1 138 345—360 单纯型 2.9 5.8 87.7 0.7 2.9 23.2 65.2 34.8 0.0 25.4 138

表 3 JOPENS—MSDP中4种定位方法的抽样概率 Tab.3 The sampling probability of four location methods in JOPENS-MSDP

图 4(Fig. 4)

图 4 爆破时刻误差概率随最大空隙角变化 (a)爆破时刻误不超过1 s的概率；(b)爆破时刻误不超过5 s的概率 Fig.4 The variation diagram of the explosion time error probability with the maximum station gap

图 5(Fig. 5)

图 5 水平误差概率随最大空隙角变化 (a)水平误差不超过1 km的概率变化；(b)水平误差不超过5 km的概率变化； (c)水平误差不超过15 km的概率变化；(d)水平误差不超过30 km的概率变化 Fig.5 The variation diagram of horizontal error probability with maximum station gap

图 6(Fig. 6)

图 6 垂直误差概率随最大空隙角变化 Fig.6 The variation diagram of the vertical error probability with the maximum station gap

由图 4可见：以90%的抽样概率为衡量基准，爆破时刻的测量误差要求不超过1 s时，Hyposat、Hyp2000、Locsat和单纯型定位算法容许的最大空隙角上限分别为240°、105°、225°和210°；爆破时刻测量误差要求不超过5 s时，定位算法容许的最大空隙角上限分别为315°、240°、255°和255°。

由图 5可见：以90%的抽样概率为衡量基准，水平误差要求不超过1 km时，4种方法均不理想，仅有Hyposat在最大空隙角小于90°时和Hyp2000在最大空隙角小于75°时抽样概率超过50%；水平误差要求不超过5 km时，4种定位算法容许的最大空隙角上限分别为225°、105°、210°和180°；水平误差要求不超过15 km时，4种定位算法容许的最大空隙角上限分别为285°、180°、240°和240°；水平误差要求不超过30 km时，4种定位算法容许的最大空隙角上限分别为315°、225°、240°和240°。

由图 6可见：垂直误差要求不超过5 km时，仅Hyp2000在最大空隙角上限不超过105°时置信度高于90%；Hyposat大多稳定在50%—70%；其他2种方法多数在50%以下。值得注意的是，在最大空隙角较小，即台站布较好时Hyposat、Locsat和单纯型定位方法的置信度偏低，结合图 3可知最大空隙角减小，使用台站增多，也就是更多远台资料被使用，可能是因为Hyp2000和Hyposat对远台料资进行权重设置，而Locsat和单纯型并未做距离权重设置，所以出现该反常现象。

综合分析可知：Hyposat在各种情况下均较稳定；Hyp2000在台站布局较好时可信度较高，但台站布局变差其结果快速变得最不可靠，即稳定性较差；Locsat和单纯型水平向表现介于Hyposat和Hyp2000之间，垂直向明显较差。

从上述2个爆破的定位结果可以看出，广东地震台网研发的地震常规分析软件JOPENS—MDSP提供的Hyposat、Hyp2000、Locsat和单纯型常规定位程序给出的定位结果之间差异不大。采用不同模型不同定位算法的定位结果与实测结果相比，最大水平误差为3.5 km，最小水平误差仅0.3 km，网内或网缘均能达到1类精度要求，垂直方向误差较大。通过分析155万多种随机抽取台站组合的定位结果，可以得到以下结论。

（1）Hyposat定位方法最稳定，具体现表为相同台站布局下Hyposat定位结果的水平向置信度最高，同一置信区间内，相同精度要求下Hyposat容许的最大空隙角上限最大；Hyp2000在台站布局好时可信度较高，但随着台站布局变差其结果快速变得最不可靠，即稳定性较差；Locsat和单纯型水平向介于两者之间，垂直向明显较差。

（2）在置信度不低于90%时，4种方法的定位结果达到1类、2类和3类精度要求的最大空隙角上限不相同。其中Hyposat最大空隙角上限分别不超过225°、285°和315°；Hyp2000最大空隙角上限分别不超过105°、180°和225°；Locsat最大空隙角上限分别不超过210°、240°和240°；单纯型最大空隙角上限分别不超过180°、240°和240°。

（3）综合考虑，Hyposat适用范围广, 稳定性好，Hyp2000的稳定性受最大空隙角的影响较大。最大空隙角小于105°时4种方法测定的发震时刻和水平方向的差别可以不考虑，但垂直向Hyp2000置信度较高，其次为Hyposat，单纯型可信度较低。

影响地震定位精度的因素较多，如台站与台站之间的距离(即台网分布密度)；台网监控范围内介质的不均匀性及地壳的分层结构；震相识别及到时量取准确程度等。本次测试主要考虑近震台网（震中距小于300 km）中台站方位分对不同定位算法的影响，其他影响因素还值得进一步分析。