颜丙乾等：基于PCA和MCMC的贝叶斯方法的海下矿山水害源识别分析 ·1415 Ⅳ类从参考文献[26]所提供的从大量的矿井 p统计表)，其中突（涌）水水源为I类水层的有17个， 突水案例中选取40个典型突水水样作为训练样本 Ⅱ类含水层的有4个，Ⅲ类含水层的有15个，W类 (如表1所示为不同水样的主要离子成分质量浓度 含水层的有4个 表1矿井突水水源判别的变量资料 Table 1 Variable data for discrimination of water inrush sources in mine 序号pHCO)MmgL)p(C)mgL-)pSO方)(mgL)p(K+Na)/(mgL-)p(Ca)(mgL-)pMg(mgL-)水样水源结果贝叶斯法判别结果 268.4 20585.5 2713.7 11707.5 716.4 1175.5 I 2 295.3 23158.4 3097.9 13025.5 651.3 1427.6 I 244.0 22644.0 3001.9 13025.5 676.4 1412.4 237.9 24702.3 2977.9 13442.5 876.8 1576.5 5 262.3 21408.6 2929.8 11957.5 651.3 1366.9 6 274.5 23981.9 3122.0 13783 751.5 1549.1 I 244.0 39112.3 4418.8 19692.5 2204.4 3049.7 m Ⅲ 8 268.4 21099.8 2881.8 11957.5 576.2 1351.7 274.5 21614.6 2833.8 12395.5 626.3 1369.9 10 286.7 24393.5 3314.1 13783 776.6 1582.5 11 268.4 27790.3 3362.1 16362.5 891.8 1767.8 12 317.3 32936.6 4130.6 18505 1052.1 2126.3 Ⅲ 13 323.4 36024.3 4514.8 19448 1262.5 2423.9 Ⅲ 14 329.5 38082.9 4682.9 21295 1067.1 2530.2 Ⅲ 15 335.6 33965.7 4156.0 19365 926.9 2098.9 m Ⅲ 16 268.4 46317.2 5283.3 25574 1703.4 3462.8 ) Ⅲ 少 280.6 29848.9 3001.9 15760 1392.8 1889.3 《 Ⅱ 18 207.4 39626.7 4755.0 21170 1853.7 3007.1 Ⅲ 19 280.6 31392.7 3962.5 17490 1112.2 2272.1 Ⅲ 219.6 33348.2 3266.0 16582 2148.3 2196.7 M 21 286.7 33183.7 3746.3 17362 1242.5 2308.5 ⑧ Ⅲ 22 250.1 28654.9 3688.7 15928 1042.1 1905.1 I 23 195.2 26184.4 3073.9 13518 1162.3 1715.6 I 24 250.1 23467.2 2881.8 12912 681.4 1458.0 25 244.0 23467.2 2881.8 12692 825.6 1492.0 250.1 23879.1 3035.5 12892 881.8 1569.8 27 250.1 43641,1 4956.7 23226 1603.2 3134.7 Ⅲ 225.7 35406.8 4169.0 19250 1362.7 2478.6 Ⅲ Ⅲ 29 268.4 34171.7 4265.1 18910 1162.3 2269.6 四 30 323.4 37053.8 4514.8 20530.8 1202.4 2527.2 ⑧ Ⅲ 31 274.5 39112.3 4923.1 21270 1468.9 2694.3 Ⅲ 32 158.6 16468.3 2439.9 9412.5 380.8 1125.1 33 152.5 16303.5 2267.1 9150 392.8 1122.7 I 1 34 317.3 20585.5 2881.8 12847.5 551.1 1184.6 I 35 225.7 19762.0 2473.5 11098.8 551.1 1233.2 I 36 176.9 17456.3 2228.6 10112.5 320.6 1093.5 37 140.3 89.0 76.8 61.2 100.2 14.6 N W 38 256.2 212.3 405.4 200 138.3 42.5 N W 39 209.9 21.3 55.7 322 81.0 4.4 V W 40 219.6 257.0 134.5 133.7 144.3 26.7 NⅣ类[25] . 从参考文献 [26] 所提供的从大量的矿井 突水案例中选取 40 个典型突水水样作为训练样本 （如表 1 所示为不同水样的主要离子成分质量浓度 ρ 统计表），其中突（涌）水水源为Ⅰ类水层的有 17 个， Ⅱ类含水层的有 4 个，Ⅲ类含水层的有 15 个，Ⅳ类 含水层的有 4 个. 表 1 矿井突水水源判别的变量资料[26] Table 1 Variable data for discrimination of water inrush sources in mine[26] 序号 ρ(HCO3 − )/(mg·L−1) ρ(Cl− )/(mg·L−1) ρ(SO4 2−)/(mg·L−1) ρ(K++Na+ )/(mg·L−1) ρ(Ca2+)/(mg·L−1) ρ(Mg2+)/(mg·L−1) 水样水源结果 贝叶斯法判别结果 1 268.4 20585.5 2713.7 11707.5 716.4 1175.5 Ⅰ Ⅰ 2 295.3 23158.4 3097.9 13025.5 651.3 1427.6 Ⅰ Ⅰ 3 244.0 22644.0 3001.9 13025.5 676.4 1412.4 Ⅰ Ⅰ 4 237.9 24702.3 2977.9 13442.5 876.8 1576.5 Ⅰ Ⅰ 5 262.3 21408.6 2929.8 11957.5 651.3 1366.9 Ⅰ Ⅰ 6 274.5 23981.9 3122.0 13783 751.5 1549.1 Ⅱ Ⅱ 7 244.0 39112.3 4418.8 19692.5 2204.4 3049.7 Ⅲ Ⅲ 8 268.4 21099.8 2881.8 11957.5 576.2 1351.7 Ⅰ Ⅰ 9 274.5 21614.6 2833.8 12395.5 626.3 1369.9 Ⅰ Ⅰ 10 286.7 24393.5 3314.1 13783 776.6 1582.5 Ⅱ Ⅱ 11 268.4 27790.3 3362.1 16362.5 891.8 1767.8 Ⅱ Ⅱ 12 317.3 32936.6 4130.6 18505 1052.1 2126.3 Ⅲ Ⅲ 13 323.4 36024.3 4514.8 19448 1262.5 2423.9 Ⅲ Ⅲ 14 329.5 38082.9 4682.9 21295 1067.1 2530.2 Ⅲ Ⅲ 15 335.6 33965.7 4156.0 19365 926.9 2098.9 Ⅲ Ⅲ 16 268.4 46317.2 5283.3 25574 1703.4 3462.8 Ⅲ Ⅲ 17 280.6 29848.9 3001.9 15760 1392.8 1889.3 Ⅱ Ⅱ 18 207.4 39626.7 4755.0 21170 1853.7 3007.1 Ⅲ Ⅲ 19 280.6 31392.7 3962.5 17490 1112.2 2272.1 Ⅲ Ⅲ 20 219.6 33348.2 3266.0 16582 2148.3 2196.7 Ⅲ Ⅲ 21 286.7 33183.7 3746.3 17362 1242.5 2308.5 Ⅲ Ⅲ 22 250.1 28654.9 3688.7 15928 1042.1 1905.1 Ⅰ Ⅰ 23 195.2 26184.4 3073.9 13518 1162.3 1715.6 Ⅰ Ⅰ 24 250.1 23467.2 2881.8 12912 681.4 1458.0 Ⅰ Ⅰ 25 244.0 23467.2 2881.8 12692 825.6 1492.0 Ⅰ Ⅰ 26 250.1 23879.1 3035.5 12892 881.8 1569.8 Ⅰ Ⅰ 27 250.1 43641.1 4956.7 23226 1603.2 3134.7 Ⅲ Ⅲ 28 225.7 35406.8 4169.0 19250 1362.7 2478.6 Ⅲ Ⅲ 29 268.4 34171.7 4265.1 18910 1162.3 2269.6 Ⅲ Ⅲ 30 323.4 37053.8 4514.8 20530.8 1202.4 2527.2 Ⅲ Ⅲ 31 274.5 39112.3 4923.1 21270 1468.9 2694.3 Ⅲ Ⅲ 32 158.6 16468.3 2439.9 9412.5 380.8 1125.1 Ⅰ Ⅰ 33 152.5 16303.5 2267.1 9150 392.8 1122.7 Ⅰ Ⅰ 34 317.3 20585.5 2881.8 12847.5 551.1 1184.6 Ⅰ Ⅰ 35 225.7 19762.0 2473.5 11098.8 551.1 1233.2 Ⅰ Ⅰ 36 176.9 17456.3 2228.6 10112.5 320.6 1093.5 Ⅰ Ⅰ 37 140.3 89.0 76.8 61.2 100.2 14.6 Ⅳ Ⅳ 38 256.2 212.3 405.4 200 138.3 42.5 Ⅳ Ⅳ 39 209.9 21.3 55.7 32.2 81.0 4.4 Ⅳ Ⅳ 40 219.6 257.0 134.5 133.7 144.3 26.7 Ⅳ Ⅳ 颜丙乾等： 基于 PCA 和 MCMC 的贝叶斯方法的海下矿山水害源识别分析 · 1415 ·