潘吉成
(神東煤炭集團(tuán)布爾臺煤礦,內(nèi)蒙古 鄂爾多斯 017000)
目前已知有16項因素影響綜放工作面著火的發(fā)生與否,分別是煤田地質(zhì)史與儲量、工作面用電安全保護(hù)措施、煤層上下附近巖石性質(zhì)、工作面周邊煤的開采進(jìn)程、采空區(qū)頂部巖石自然垮落的進(jìn)度、煤的性質(zhì)與發(fā)火期長短、工作面火災(zāi)警報系統(tǒng)、工作面著火安全應(yīng)急管理措施、含矸量低的煤層可開采厚度、煤層的傾斜角度、煤層距離地面所埋高度、煤層周圍巖石的溫度大小、所使用采煤方法的出煤率、工作面通風(fēng)系統(tǒng)、工作面每日的出煤量、工作面的采掘速率[1-3],可分別用X1~X16表示。根據(jù)以上諸多因素可以得知,綜放工作面著火原因的綜合評價是一個繁雜困難且模糊的過程。因此,需要確定出各影響因素對引起綜放工作面著火原因的權(quán)重,如果可以將引起綜放工作面著火的各個因素根據(jù)其所占比例大小進(jìn)行排序,安全管理工作者就可以對綜放工作面著火問題更加有針對性地進(jìn)行防范治理工作。
目前有專家調(diào)查法、相對比較法、層次分析法等方法可以將影響著火的因素進(jìn)行相關(guān)性大小的排序評判[4-5]。專家調(diào)查法是邀請一定數(shù)量經(jīng)驗豐富的專家,針對每個影響因素單獨(dú)發(fā)表意見,并打分評價,最后將分?jǐn)?shù)整理綜合,就可以得到每個影響因素的平均值即代表其與著火原因的相關(guān)性大?。幌鄬Ρ容^法與序列比較法相仿,即對影響因素進(jìn)行兩兩比較,然后把比較后的各因素分值相加以此作為各指標(biāo)的權(quán)重。層次分析法,是近些年誕生的一種比較科學(xué)的評判方法,其克服了前述2種方法的缺點(diǎn),充分考慮了各影響因素之間的內(nèi)在聯(lián)系,不容易受到專家知識水平和經(jīng)驗的限制,可以比較真實地處理這類問題。鑒于此,運(yùn)用層次分析法評價各種因素對引起綜放工作面著火的權(quán)重,以期定量地給出影響綜放工作面著火各因素的權(quán)重因子。
決策問題的應(yīng)用過程一般采用定量分析方法,但是此方法還存在著一些問題有待完善。由于多因素影響,需要決策的問題一般較為復(fù)雜,這樣定量分析法難以構(gòu)造出適合解決特定問題的模型,而且通過經(jīng)驗公式構(gòu)造出的模型難以根據(jù)特定的問題進(jìn)行說明并執(zhí)行程序。此外,在決策問題的過程中決策者通常帶有相當(dāng)一部分的主觀性,從而難以確定最優(yōu)化模型。針對上述定量分析方法所存在的問題,決策者重新思考社會科學(xué)中數(shù)量方法所起到的作用,特別是在決策問題方面,如何將規(guī)范的數(shù)學(xué)分析所具備的精確性與人類大腦產(chǎn)生的靈活的決策思維過程及思維規(guī)律有機(jī)地結(jié)合起來,簡單說就是將定性與定量的方法巧妙連接,由此層次分析法應(yīng)運(yùn)而生。層次分析法首次提出于20世紀(jì)70年代,提出者為美國運(yùn)籌學(xué)家SAATY T L教授,該方法的提出旨在找到一種更為簡便的方式來解決一些較為復(fù)雜、模糊的問題并作出決策處理,特別適用于那些在定量方面不能完全分析的問題。
層次分析法的英文縮寫為AHP,1971年美國國防部研究“應(yīng)急計劃”時,SAATY T L提出該方法并加以應(yīng)用,為AHP法的首次亮相;1977年國際數(shù)學(xué)建模會議,SAATY T L發(fā)表了“無結(jié)構(gòu)決策問題的建模—層次分析法”一文,至此AHP法進(jìn)入問題決策方法領(lǐng)域內(nèi)眾人的視野,得以在決策問題的許多領(lǐng)域得到應(yīng)用,并且自身也在不斷地深化和發(fā)展[6-8]。1982年中美召開圍繞能源、資源、環(huán)境等方面問題的會議,會議上SAATY T L教授的學(xué)生GHOLAMNEZHAD H向中國學(xué)者介紹了AHP這種新型決策方法,較傳統(tǒng)決策方法更為準(zhǔn)確與便利,AHP法借此契機(jī)傳入中國。通過創(chuàng)新思想的學(xué)者許樹柏等人的更深層次研究,并與本國實際問題相結(jié)合,發(fā)表了國內(nèi)AHP法的首篇文章“層次分析法—決策的一種實用方法”。文章發(fā)表之后,AHP法在我國得到迅速發(fā)展,并且在1987年9月我國召開了第1屆AHP法學(xué)術(shù)討論會,次年我國順利邀請到國際該領(lǐng)域知名專家共同參加了第1屆國際AHP學(xué)術(shù)會議,探討存在問題與發(fā)展前景。時至今日,AHP法已基本成熟,但其應(yīng)用和理論方面還在隨著社會經(jīng)濟(jì)科學(xué)的發(fā)展不斷地進(jìn)步與完善。
層次分析法的基本原理是先列舉影響事件發(fā)生的因素,然后根據(jù)各種因素與事件發(fā)生的相關(guān)性大小進(jìn)行排序,得到包含先后順序的一組數(shù)據(jù),并且為決策提供數(shù)據(jù)支持,這種原理叫排序原理。具體來說,該方法首先將需要進(jìn)行決策的問題看作是一個由多種因素共同作用的大系統(tǒng),將系統(tǒng)中的各個因素按照相互之間關(guān)聯(lián)性的大小分成一系列階層并進(jìn)行高低排序,即為優(yōu)劣排序,得到構(gòu)造層次結(jié)構(gòu)圖。之后再請相關(guān)領(lǐng)域?qū)<覍ε判蚪Y(jié)果進(jìn)行比較分析,判斷該排序的現(xiàn)實性與可行性,進(jìn)行輔助決策。層次分析法之所以從眾多方法中脫穎而出,最重要的原因是它可以合理地融合定性與定量分析,做到將數(shù)學(xué)方式的科學(xué)性和人的主觀感知的靈活性結(jié)合起來,并摒棄數(shù)學(xué)方式的死板與主觀感知的不可靠性,較為精確地反映復(fù)雜的社會科學(xué)領(lǐng)域的實際情況。雖然該法的理論基礎(chǔ)非常深厚,但其表現(xiàn)形式簡易明了,易于理解與操作,因此這一方法得到了廣泛應(yīng)用。
一般由4個步驟來完成層次分析法的建模。首先建立進(jìn)階層次結(jié)構(gòu)的評價模型,接著構(gòu)造出各層次中的所有判斷矩陣,然后進(jìn)行層次單排序及一致性檢驗,最后完成層次總排序及一致性檢驗[9-10]。
為了明確引起綜放面著火原因的各種影響因素各自的影響權(quán)重,在前述綜合分析的基礎(chǔ)之上構(gòu)建評價層次結(jié)構(gòu)模型。設(shè)定A為目標(biāo)層、B為準(zhǔn)則層、X為子準(zhǔn)則層,由左至右展開,為了使數(shù)據(jù)更為直觀,操作更為簡便,規(guī)定同一層次的不同因素從上至下用X1~X16代表,得到層次結(jié)構(gòu)的模型,見表1。
表1 綜放面綜合評價火災(zāi)危險性的分層結(jié)構(gòu)模型Table 1 Layered structure model for comprehensive evaluation of fire risk in fully mechanized caving face
經(jīng)過前人長期積累的經(jīng)驗,綜合考慮陽泉礦區(qū)綜放面著火原因的實際情況,將影響綜放面著火的各種因素逐層分解。采用SAATY T L教授所創(chuàng)建的9標(biāo)度法,分別對同一層次中各個影響事件發(fā)生的因素進(jìn)行對比排序,通過排序結(jié)果對因素進(jìn)行評價[11-12],9標(biāo)度法含義見表2。其依據(jù)是通過對心理學(xué)的研究,進(jìn)行大量測試與問卷調(diào)查得出結(jié)果,大多數(shù)人對屬性相同的事物有1~9級的感應(yīng)差別;此外,隨著科技發(fā)展,人們對安全知識認(rèn)識水平越來越高,9標(biāo)度法的使用越來越普遍,越來越多的專家對該系統(tǒng)方法進(jìn)行改進(jìn)使之趨向成熟;目前,9標(biāo)度法能夠完全找出因素中引起人感知差別的方面。
表2 9標(biāo)度法含義Table 2 Meaning of 9-scale method
賦值結(jié)果見表3。
表3 綜放面著火原因通過層次分析法得出的矩陣模型Table 3 The matrix model of the cause of fire in fully mechanized caving face by AHP
可通過MATLAB軟件編程得到矩陣的一般屬性,再對矩陣進(jìn)行一致性檢驗和特征向量歸一化,經(jīng)過軟件計算,得到該矩陣的權(quán)系數(shù)向量。其求解結(jié)果如下
μA=(0.895 7,0.376 2,0.237 0)T;
μB1=(-0.701 2,-0.064 5,-0.598 9,-0.273 8,-0.244 2,-0.072 6,-0.074 5)T;
μB2=(0.557 3,0.732 7,0.127 2,0.119 9,0.298 0,0.182 1)T;
μB3=(-0.879 8,-0.120 2,-0.459 9)T。
通過以上數(shù)據(jù)進(jìn)行各矩陣特征向量μA,μB1,μB2,μB3的歸一化處理,即可得到在單一條件約束下,矩陣中每個元素之間相對排序權(quán)重的大小,即向量ωA,ωB1,ωB2,ωB3,具體結(jié)果如下:ωA=(0.593 6,0.249 3,0.157 1)T;ωB1=(0.345 5,0.031 8,0.295 1,0.134 9,0.120 3,0.035 8,0.036 7)T;ωB2=(0.276 3,0.363 2,0.063 1,0.059 4,0.147 7,0.090 3)T;ωB3=(0.602 6,0.082 3,0.315 0)T。
一致性檢驗所應(yīng)用到的經(jīng)典公式
(1)
式中,CR為一致性比率;CI為一致性指標(biāo);RI為一致性檢驗的平均隨機(jī)指標(biāo),見表4;λmax為判斷矩陣的最大特征值。
表4 一致性檢驗的平均隨機(jī)指標(biāo)表Table 4 Average random index table of consistency test
通過軟件處理,可得到各個矩陣的一致性檢驗結(jié)果,即
CRA=0.046 2;CRB1=0.069 8;CRB2=0.042 4;CRB3=0.001 7。通過觀察可以發(fā)現(xiàn),以上矩陣的CR值都低于0.1。因此可以認(rèn)為各個判斷矩陣所得結(jié)果皆為滿意一致性。
層次分析法歸根結(jié)底是要得到方案層中每個元素相對于目標(biāo)層的相關(guān)性大小的順序排列,這就需要在階層概念上從上方向下方依次進(jìn)行每個層面中所包含各種元素對目標(biāo)的相關(guān)性大小的計算。計算結(jié)果為
ω=(0.205 1,0.018 9,0.175 2,0.080 1,0.071 4,0.021 3,0.021 8,0.068 9,0.090 5,0.015 7,0.014 8,0.036 8,0.022 5,0.094 7,0.012 9,0.049 5)T
需要對層次分析法進(jìn)行一致性檢驗來確定各層中各個元素對目標(biāo)層的向量的相關(guān)性大小所達(dá)到的滿意接受的程度高低進(jìn)行綜合排列順序,其檢驗過程為
CI=(0.069 8,0.042 4,0.001 7)×ωA
=0.052 3
RI=(1.32,1.24,0.58)×ωA=1.183 8
通過檢驗得出層面水平相關(guān)性大小排序達(dá)到滿意的一致性,最終得到各因素對影響綜采放頂煤工作面著火事件發(fā)生的權(quán)重因子,如圖1所示,以下是具體數(shù)值:煤田地質(zhì)史與儲量(0.205 1)、煤層上下附近巖石性質(zhì)(0.018 9)、工作面周邊煤的開采進(jìn)程(0.068 9)、采空區(qū)頂部巖石自然垮落的進(jìn)度(0.090 5)、工作面用電安全保護(hù)措施(0.094 7)、工作面火災(zāi)警報系統(tǒng)(0.012 9)、工作面著火安全應(yīng)急管理措施(0.049 5)、煤的性質(zhì)與發(fā)火期長短(0.175 2)、含矸量低的煤層可開采厚度(0.080 1)、煤層的傾斜角度(0.071 4)、煤層距離地面所埋高度(0.021 3)、煤層周圍巖石的溫度大小(0.021 8)、所使用采煤方法的出煤率(0.015 7)、工作面每日的出煤量(0.014 8)、工作面通風(fēng)系統(tǒng)(0.036 8)、工作面的采掘速率(0.022 5)。
圖1 各因素對煤層自燃造成影響的相關(guān)性大小Fig.1 The correlation of various factors on coal spontaneous combustion
在綜放面著火的各種因素的基礎(chǔ)之上,結(jié)合層次分析法的思想構(gòu)建了綜放面著火的綜合評價模型,子準(zhǔn)則層X中每個元素相對目標(biāo)層A的相關(guān)性大小排列順序可以通過定量的數(shù)學(xué)分析方法得到更加準(zhǔn)確的表達(dá)。通過權(quán)重分析可發(fā)現(xiàn),煤層賦存地質(zhì)構(gòu)造(0.205 1)、煤自然發(fā)火性(0.175 2)、工作面電氣設(shè)備及防爆情況(0.094 7)、采空區(qū)冒落及充填情況(0.090 5)、煤層厚度(0.080 1)這5個因素是引起礦區(qū)綜采放頂煤工作面著火的主要因素。當(dāng)然,對于具體工作面的火災(zāi)案例,上述5個因素或5個因素中的某幾個因素相互耦合,共同促成了工作面火災(zāi)。