黃金波 王軼多
(浙江省建工集團有限責任公司, 浙江 杭州 310012)
隨著城市化進程的不斷推進,建筑業(yè)也不斷發(fā)展,并不斷嘗試在采空區(qū)上方進行建設[1]。但采空區(qū)存在地表變形,會威脅建筑物的安全。因此,為了采空區(qū)建筑的安全使用,對采空區(qū)建筑地基進行穩(wěn)定性評價顯得尤為主要[2]。
長期以來,許多專家學者對采空區(qū)建筑地基的穩(wěn)定性進行了研究。其中,劉洋等[3]利用層次分析法和模糊評價理論對采空區(qū)建筑地基穩(wěn)定性進行綜合評價,所得評價結(jié)果與實際礦山情況相吻合;趙超等[4]將層次分析法與可拓評價理論法相結(jié)合,建立了建筑地基穩(wěn)定性評價的物元模型,準確地判別出穩(wěn)定性等級,其穩(wěn)定性評價結(jié)果對于采空區(qū)的治理提供了科學、有效的依據(jù)。上述研究雖然取得了一定的成果,但對采空區(qū)建筑地基進行穩(wěn)定性評價是一個復雜的過程[5],僅采用單一賦權(quán)方法具有一定的主客觀性,降低評估結(jié)果的可靠性[6]。
為此,本文先結(jié)合已有研究構(gòu)建了建筑地基穩(wěn)定性評價指標體系,并基于層次分析法、改進的標準間重要性相關性(criteria importance though intercrieria correlation,CRITIC)法對評價指標進行主客觀賦權(quán),最后再采用灰云模型理論,建立了一套建筑地基穩(wěn)定性評價模型,避免了專家主觀決策對評價結(jié)果的影響,從而獲得了科學、準確的建筑地基穩(wěn)定性評價結(jié)果。
由于影響采空區(qū)上方建筑地基穩(wěn)定性的指標眾多,且部分指標存在一定的冗余性,為了對建筑地基合理地進行穩(wěn)定性評價,通過綜合考慮采空區(qū)的實際情況以及參考相關文獻,選取出涵蓋停采時間、深厚比等7個指標構(gòu)建了穩(wěn)定性評價指標體系,如圖1所示[7-8]。將采空區(qū)建筑地基穩(wěn)定性劃分為Ⅰ、Ⅱ、Ⅲ和Ⅳ級,分別代表穩(wěn)定、基本穩(wěn)定、欠穩(wěn)定和不穩(wěn)定,各評價指標分級標準如表1所示[9]。
圖1 采空區(qū)建筑地基穩(wěn)定性評價指標體系
表1 評價指標分級表
在本文中,以層次分析法對評價指標進行主觀賦權(quán),依據(jù)1~9標度表[10]對評價指標進行兩兩比較,確定評價指標的相對重要性,從而構(gòu)建判斷矩陣F,如式(1)所示[11]。
(1)
式中,flt為指標l與指標t的比值;n為指標個數(shù)。
根據(jù)矩陣F計算最大特征根λmax相對應的特征向量,并對向量進行歸一化處理,得到各評價指標的權(quán)重P,實現(xiàn)對評價指標的主觀賦權(quán)。為保證所得權(quán)重的合理性,必須對矩陣進行一致性檢驗。一致性比例CR可定義為
式中,CI為一致性指標;RI為平均隨機一致性指標。
RI隨著評價指標數(shù)量n的變化而變化,其具體取值如表2所示。當一致性比例CR<0.1時,認為判斷矩陣F具有良好的一致性,若CR>0.1則需對判斷矩陣F進行適當調(diào)整,直至其滿足一致性。
表2 平均隨機一致性指標RI
CRITIC不僅考慮了指標之間的信息量,而且顧及了指標之間的相關性。但是,不同指標間的量綱、數(shù)量級往往不同,CRITIC法用標準差衡量指標的變異性存在一定的不足[12]。因此,采用變異系數(shù)對CRITIC法進行改進,更好地實現(xiàn)指標賦權(quán),其主要步驟如下[13]:
(1)以原始評價指標值構(gòu)建矩陣X,如式(4)所示。
(4)
式中,xij為第i個評價對象第j個評價指標值;m為評價對象數(shù)目;n為評價指標數(shù)目。
(2)對矩陣X中的元素進行標準化處理,得到標準化矩陣X*=(xij*),標準化計算如式(5)所示。
(5)
(3)計算指標體系中各指標的變異系數(shù)vj,如式(6)所示。
(6)
(4)通過SPSS軟件,可計算出矩陣X*中各指標的Pearson相關系數(shù),并構(gòu)造出相關系數(shù)矩陣R=(rkj)n×n(k=1,2,…,n;j=1,2,…,n),再根據(jù)式(7)計算各指標的獨立性系數(shù)ηj。
(7)
(5)結(jié)合所求的變異系數(shù)和獨立性系數(shù),由式(8)計算各指標的重要性系數(shù)Hj。
(8)
最終可通過式(9)計算第j個評價指標的權(quán)重Qj。
(9)
基于博弈論的綜合賦權(quán)法,將由層次分析法確定的主觀權(quán)重與由改進的CRITIC確定的客觀權(quán)重相結(jié)合,從而得到評價指標的綜合權(quán)重[14]。具體算法如下:
對于多指標評價系統(tǒng),假設采用L種方法得到指標的權(quán)重,即權(quán)重向量W為
(10)
以此得到一個權(quán)重集W,L種權(quán)重向量的線性組合如式(11)表示。
(11)
式中,αk為線性組合系數(shù),αk>0;W為所有的權(quán)重向量集。
為了選擇出最合適的權(quán)重向量W*,需對線性組合系數(shù)進行優(yōu)化,使得W與每個wk的離差極小化,并導出其目標函數(shù)。
(12)
按照矩陣微分性質(zhì),與上式等價的最優(yōu)一階導數(shù)條件為
(13)
對式(14)進行求解可得(α1,α2,…,αL),然后通過式(15)對其進行歸一化處理,即可得到綜合權(quán)重向量W*。
灰色模型是在云模型理論的基礎上引進灰色理論,解決了模糊概念的定量化處理,從而實現(xiàn)定性與定量之間的轉(zhuǎn)換[15]?;以颇P筒捎弥笜朔旨壍纳舷逻吔?Tx,Px)、峰值Rx、熵En和超熵He等數(shù)字特征進行轉(zhuǎn)換。不同穩(wěn)定性等級所對應的灰云模型數(shù)字特征的計算方法如表3所示,其中k為常數(shù),其可以調(diào)節(jié)云模型云化程度,通常取值為7[16]。
表3 灰云模型數(shù)字特征
根據(jù)表3可確定出灰云模型的數(shù)字特征,采用MATLAB軟件可生成評價指標的灰云模型,再結(jié)合博弈論所得的組合權(quán)重,基于最大隸屬度原則便可得到各評價對象所隸屬于各穩(wěn)定性等級的確定度。
本文以徐州夾河煤礦與龐莊煤礦的井田為例[17],運用所構(gòu)建的采空區(qū)建筑地基穩(wěn)定性評價體系以及模型,對該研究區(qū)域進行穩(wěn)定性評價。
首先利用層次分析法,對7個評價指標進行兩兩比較,所構(gòu)造的判斷矩陣為
(16)
層次分析法確定的權(quán)重如表4所示。根據(jù)判斷矩陣F計算得到的特征值λmax=7.249 4和一致性比例CR=0.030 6<0.1,判斷矩陣滿足一致性,所得的權(quán)重有效。
再根據(jù)改進的CRITIC方法計算評價指標權(quán)重,結(jié)果如表4所示。根據(jù)博弈論的組合原則,得到各評價指標的組合權(quán)重,如表4所示。
表4 評價指標權(quán)重
根據(jù)灰云模型理論,編寫MATLAB程序求解各建筑地基隸屬于各穩(wěn)定等級的確定度,并基于最大隸屬度原則確定出各建筑地基的穩(wěn)定性等級,結(jié)果如表5所示。
表5 建筑地基穩(wěn)定性等級判定
基于組合賦權(quán)法和灰云模型理論構(gòu)建了建筑地基穩(wěn)定性評價模型,并將其應用于具體工程實例中,其主要結(jié)論如下:
(1)采用層次分析法和改進的CRITIC法分別確定評價指標的主客觀權(quán)重,能有效減小單一賦權(quán)造成的誤差,消除指標之間的差異性,提高評判結(jié)果的可靠性。
(2)采用灰云模型理論構(gòu)建地基穩(wěn)定性評價模型,實現(xiàn)定性指標與定量數(shù)據(jù)的轉(zhuǎn)換,從而得到各地基的穩(wěn)定性等級,綜合反映了地基穩(wěn)定性評判過程中的灰性和隨機性,為采空區(qū)建筑地基穩(wěn)定性的評價提供了一種新方法。