陸春輝，徐 翀，陳永春

(1.煤礦生態(tài)環(huán)境保護(hù)國(guó)家工程實(shí)驗(yàn)室，安徽淮南232033;2.淮南礦業(yè)(集團(tuán))有限責(zé)任公司，安徽淮南232001)

突水水源判別一直是礦山防治水研究熱點(diǎn)內(nèi)容之一，研究方法多種多樣，其中依據(jù)不同含水層的水化學(xué)成分特征，最直接反應(yīng)地下水運(yùn)動(dòng)特征，利用不同的數(shù)學(xué)方法，將含水層水源數(shù)據(jù)進(jìn)行編程處理，從而得到判別效果［1］－［3］，諸如模糊評(píng)判［4］－［6］、灰色關(guān)聯(lián)法［7］、數(shù)量化理論［8］、Bayes 逐步判別［9］、距離判別分析［10］、BP神經(jīng)網(wǎng)絡(luò)分析［11］等。通常情況礦井突水含水層為若干個(gè)，且水質(zhì)特征界限不是十分明顯，難以通過單一因素進(jìn)行準(zhǔn)確判別，而模糊評(píng)判專門針對(duì)此類問題，且原理簡(jiǎn)單、易于實(shí)現(xiàn)，從而在礦山前期突水水源判別中得到了廣泛的應(yīng)用。本文利用模糊評(píng)判方法，結(jié)合淮南某礦山的突水水源為研究實(shí)例，其結(jié)果為類似條件提供借鑒與參考。

1 模糊評(píng)判方法的原理

模糊評(píng)判方法是一種應(yīng)用模糊數(shù)學(xué)原理對(duì)多種因素所影響的事物或現(xiàn)象進(jìn)行綜合評(píng)價(jià)的方法。它實(shí)際上是一種多目標(biāo)決策的模糊數(shù)學(xué)研究方法，其過程是通過模糊變換來完成的，原理為:

式中A=(a1，a2，Λ，ap)為模糊權(quán)向量，它是因素集 X={x1，x2，Λ，xp}上的權(quán)重分配;R為因素與評(píng)價(jià)集的模糊關(guān)系矩陣。若評(píng)價(jià)集為 u={G1，G2，Λ，Gk}，則 R=(rij)p×k表示 xi對(duì) Gj的隸屬度;B=(b1，b2，Λ，bk)為評(píng)判指標(biāo)，表示綜合考慮所有因素影響時(shí)評(píng)價(jià)對(duì)象對(duì)評(píng)判集的隸屬度?！?”表示A與R的一種合成方式，這里取

1)建立評(píng)判集合。設(shè)有 k個(gè)總體 G1，G2，Λ，Gk，則模糊評(píng)判集可記為 u={G1，G2，Λ，Gk}。

2)建立主要影響因素集合。取能夠比較客觀地反映總體特征的p個(gè)指標(biāo)，構(gòu)成主要影響因素集合X=(x1，x2，Λ，xp)。

3)建立評(píng)判集各因子的標(biāo)準(zhǔn)值。從每個(gè)總體Gi中抽取ni個(gè)樣品，每個(gè)樣品測(cè)p個(gè)指標(biāo)，如第i個(gè)總體的第j個(gè)樣品。用第i個(gè)總體Gi的樣本均值生成該總體的評(píng)判因子標(biāo)準(zhǔn)值:

4)隸屬度計(jì)算及模糊評(píng)價(jià)矩陣構(gòu)建。隸屬度一般用rij表示，可通過隸屬函數(shù)計(jì)算求得

式中:rij為第i個(gè)因子對(duì)第類評(píng)價(jià)的隸屬度;xi為第i個(gè)因子的實(shí)測(cè)值為第i個(gè)因子所對(duì)應(yīng)的第j類的評(píng)價(jià)標(biāo)準(zhǔn)值。由rij可構(gòu)成模糊評(píng)價(jià)矩陣

5)評(píng)判因子權(quán)重的計(jì)算。綜合評(píng)判結(jié)果的正確程度，因子權(quán)重的確定是一個(gè)很重要的因素。本次研究采用超標(biāo)加權(quán)法計(jì)算各因子權(quán)重，即變權(quán)法，如(4)式表示:

式中，A(m，n)為第m個(gè)總體、第n個(gè)評(píng)判因子的權(quán)重;x(m，n)為總體第m第n個(gè)因子的實(shí)測(cè)值;xn為評(píng)判集第n個(gè)因子的平均值;p為評(píng)判因子個(gè)數(shù)。

6)模糊評(píng)判。將模糊權(quán)向量A與模糊評(píng)價(jià)矩陣R合成，得到B=A*R=(b1，b2，Λ bk)，若 bj=max{b1，b2，Λ bk}，則未知樣品判歸為第j個(gè)總體。

2 應(yīng)用效果

已有研究認(rèn)為，因素權(quán)重和隸屬函數(shù)的確定對(duì)模糊評(píng)判結(jié)果影響明顯，并對(duì)權(quán)重計(jì)算方法進(jìn)行了一些探討［4］－［6］，而相關(guān)學(xué)者在研究其它判別方法時(shí)，發(fā)現(xiàn)樣本數(shù)據(jù)量及判別因子選取同樣會(huì)對(duì)判別結(jié)果具有影響［8］－［11］。

2.1 突水水源水質(zhì)分析

該礦直接威脅突水水源為:煤層為頂、底板砂巖水、灰?guī)r水及采空區(qū)水。本次采集并分析了15個(gè)已知水樣數(shù)據(jù)，經(jīng)過分析建立標(biāo)準(zhǔn)評(píng)判集，如表1所示，其中砂巖裂隙水、灰?guī)r水和老空水對(duì)應(yīng)的模型輸出結(jié)果分別為I、II、III。利用模糊評(píng)判法對(duì)8個(gè)未知水樣數(shù)據(jù)進(jìn)行歸類判別，如表2。

表1 某礦采集的水源樣本數(shù)據(jù) mg/L

表2 待判水源數(shù)據(jù) mg/L

1)評(píng)價(jià)集及因素集建立。依該突水水源情況及判別要求，建立評(píng)價(jià)集為:

如前所述，主要影響因素集為能較客觀地反映總體特征指標(biāo)集合，根據(jù)以往研究結(jié)果，首先選取 Na++K+，Ca2+，Mg2+，Cl－，SO4

2－，及 HCO3－等構(gòu)建影響因素集。

2)評(píng)判集各因子的標(biāo)準(zhǔn)值。根據(jù)上述表1訓(xùn)練結(jié)果，評(píng)價(jià)集因子標(biāo)準(zhǔn)值參考如表3。

表3 評(píng)價(jià)集因子標(biāo)準(zhǔn)值 mg/L

3)模糊評(píng)價(jià)矩陣的構(gòu)建。在表3的評(píng)判集標(biāo)準(zhǔn)值基礎(chǔ)上，結(jié)合式(2)、(3)式即可建立各待判樣本的隸屬函數(shù)及模糊評(píng)價(jià)矩陣。

4)評(píng)判因子權(quán)重。按照公式(4)，各待判樣本的評(píng)判因子權(quán)重如表4。

表4 評(píng)價(jià)集因子權(quán)重

5)模糊評(píng)判。將模糊權(quán)向量與模糊評(píng)價(jià)矩陣相乘，即可得到各待判因子的評(píng)判集B，如樣本1的評(píng)判集，如下:

根據(jù)最大隸屬度原則，樣本1屬于II類，即灰?guī)r水，其余樣本判別結(jié)果如表5所示。

表5 待判樣本模糊評(píng)判結(jié)果

由表4可以看出，樣本4、7出現(xiàn)了誤判，其余與實(shí)際相符，準(zhǔn)確率80%。

2.2 水源樣本數(shù)量對(duì)評(píng)判結(jié)果影響

此處將水源樣本數(shù)量增加至36個(gè)，重新進(jìn)行評(píng)價(jià)矩陣構(gòu)建及待判樣本分析，判別結(jié)果如表6所示。

表6 增加樣本數(shù)據(jù)后的模糊評(píng)判結(jié)果

由表6可以看出，增加樣本數(shù)據(jù)量后，評(píng)判結(jié)果并沒有明顯改善，反而樣本2、4、7均出現(xiàn)了誤判，由此可見模糊評(píng)判準(zhǔn)確率與樣本數(shù)據(jù)多少并無直接關(guān)聯(lián)，只需選取一定量的典型樣本。

2.3 評(píng)判因子選取對(duì)評(píng)判結(jié)果影響

按前述采用 Na++K+，Ca2+，Mg2+，Cl－，SO2－，4作為判別因子，判別準(zhǔn)確率相對(duì)較低，現(xiàn)增加總硬度、堿度、PH值和礦化度等為新判別因子，以增加樣本數(shù)據(jù)之間的差異性，進(jìn)而改善模糊評(píng)判準(zhǔn)確性，表7為增加判別因子后的15個(gè)樣本數(shù)據(jù)。

增加新判別因子后，建立新的評(píng)價(jià)矩陣進(jìn)行待判樣本評(píng)判，其判別結(jié)果如表8所示。

表8 增加判別因子后的模糊評(píng)判結(jié)果

由此看出，表明增加新的評(píng)判因子數(shù)量后，判別結(jié)果與實(shí)際完全相符。

表7 增加新判別因子

3 結(jié)語(yǔ)

模糊評(píng)判分析方法是該礦突水水源判別的有效方法之一，其樣本數(shù)據(jù)量的多少與模糊評(píng)判分析準(zhǔn)確性沒有直接關(guān)系，而判別結(jié)果準(zhǔn)確性高低與選取的樣本數(shù)據(jù)是否具有代表性有關(guān)，增加適宜的評(píng)判因子，如總硬度、堿度、PH值、礦化度等對(duì)提高水源判別準(zhǔn)確率至關(guān)重要。

［1］李燕，徐志敏，劉勇.礦井突水水源判別方法概述［J］.煤炭技術(shù)，2010，29(11):87 －89

［2］王廣才，王秀輝，李競(jìng)生，等.平頂山礦區(qū)礦井突(涌)水水源判別模式［J］.煤田地質(zhì)與勘探，1998，26(3):47 －50.

［3］王連國(guó)，宋揚(yáng).礦井突水危險(xiǎn)性評(píng)價(jià)模型［J］.工程地質(zhì)學(xué)報(bào)，2000，8(2):160－163.

［4］楊永國(guó)，李賓亭，尚克勤，等.鶴壁礦務(wù)局礦井突水等級(jí)模糊綜合評(píng)判及預(yù)測(cè)［J］.中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，1998，27(2):204 －208.

［5］賁旭東，郭英海，解奕偉.模糊綜合評(píng)判在礦井突水水源的應(yīng)用及探討［J］.礦業(yè)安全與環(huán)保，2006，6:57 －59.

［6］李明山，禹云雷，路風(fēng)光.姚橋煤礦礦井突水水源模糊綜合評(píng)判模型［J］.勘察科學(xué)技術(shù)，2001，2:16－19.

［7］關(guān)英斌，鄭建，豐成.灰色關(guān)聯(lián)分析在牛西礦水源判別中的應(yīng)用［J］.河北工程大學(xué)學(xué)報(bào)，2011，28(1):81－84.

［8］張?jiān)S良，張子戌，彭蘇萍.數(shù)量化理論在礦井突(涌)水水源判別中的應(yīng)用［J］.中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，2003，32(3):251 －254.

［9］張春雷，錢家忠，趙衛(wèi)東.Bayes方法在礦井突水水源判別中的應(yīng)用［J］.煤田地質(zhì)與勘探，2010，38(4):34－37.

［10］王心義，徐濤，黃丹.距離判別法在相似礦區(qū)突水水源識(shí)別中的應(yīng)用［J］.煤炭學(xué)報(bào)，2011，36(8):1354－1358.

［11］錢家忠，呂純，趙衛(wèi)東.Elman與BP神經(jīng)網(wǎng)絡(luò)在礦井水源判別中的應(yīng)用［J］.系統(tǒng)工程理論與實(shí)踐，2010，30(1):145－149.