基于模糊數(shù)學(xué)的露天鐵礦采空區(qū)地球物理探測評價研究

賈三石 付建飛 王恩德 郭 凱 門業(yè)凱

（1.東北大學(xué)秦皇島分校資源與材料學(xué)院，河北秦皇島066004；2.東北大學(xué)資源與土木工程學(xué)院，遼寧沈陽110819）

采空區(qū)作為一種礦產(chǎn)資源開采后的遺留產(chǎn)物，大量存在于金屬和非金屬礦山［1］，特別是在鐵礦開采歷史悠久的鞍本地區(qū)表現(xiàn)得最為典型和突出。該地區(qū)既有礦業(yè)整合開發(fā)前偷采濫挖、采富棄貧和無序開采留下的采空區(qū)，又有日偽時期掠奪式開采留下的采空區(qū)，且它們數(shù)量龐大、規(guī)模不等和賦存空間不明。但從20 世紀(jì)末以來，隨著資源開發(fā)利用的科學(xué)技術(shù)水平不斷提高，主要是采掘運輸能力和選礦工藝水平的大幅提升，礦產(chǎn)資源整合露天開采成為一種趨勢，具有開采量大、礦石回收率高、生產(chǎn)安全、效益高、成本低等優(yōu)點［2］。由于歷史原因和資源開采方式的轉(zhuǎn)變，露天礦山多由散亂小規(guī)模井下開采變?yōu)檎洗笠?guī)模露天開采，導(dǎo)致采場多處于不明采空區(qū)的發(fā)育區(qū)。而采空區(qū)頂板厚度會隨著持續(xù)的開采擾動不斷減小，直至突然塌陷形成深坑，輕則造成采掘設(shè)備墜毀，重則造成生產(chǎn)運輸車輛和人員損失。且隨著露天采礦活動的持續(xù)進(jìn)行，采坑深度不斷加大，形成了高陡邊坡作用下的礦山開采應(yīng)力場不穩(wěn)多變，導(dǎo)致礦區(qū)采空區(qū)安全問題日益突出，不斷誘發(fā)坍塌、變形移位或含水空區(qū)突水等礦山次生地質(zhì)災(zāi)害，成為時刻威脅露天金屬礦山生產(chǎn)安全的主要難題［3-6］。

為解決上述遺留采空區(qū)誘發(fā)的生產(chǎn)安全技術(shù)難題，近年來，科研院所和生產(chǎn)單位采用了多種技術(shù)方法和手段對采空區(qū)進(jìn)行探測，如人工調(diào)查法、地球物理探測、工程鉆探、放射性氡氣測量、三維激光掃描和水下聲吶測量等［7］，尤以地球物理探測方法最受青睞。該類方法具有種類多、成本低、方便快捷高效、無損非接觸的地質(zhì)CT 探測成像等特點，屬于超前的“綠色探測”范疇，適用于對面積性采場區(qū)域內(nèi)的采空區(qū)異常進(jìn)行快速掃面篩查和定位，確保后期高效的鉆探驗證和三維激光或水下聲吶精準(zhǔn)測量工作有效開展。但地球物理探測成果影像顯示異常具有間接性、多解性和圈定指標(biāo)的不確定性，這使得探測成果的解譯和應(yīng)用需要有豐富經(jīng)驗的地質(zhì)與地球物理專家級的人才聯(lián)合參與?，F(xiàn)階段，在露天礦山快速持續(xù)不斷的平臺作業(yè)生產(chǎn)中，每天采集的大量多方法物探影像數(shù)據(jù)難以得到具有豐富經(jīng)驗的專家級人才及時處理，致使采礦作業(yè)平臺采空區(qū)異常得不到及時、準(zhǔn)確圈定，造成鉆探驗證和采空區(qū)安全處理工作嚴(yán)重滯后，在制約礦山正常生產(chǎn)的同時威脅了礦山安全。針對該類應(yīng)用地球科學(xué)問題僅依靠專家經(jīng)驗解決的局限性，隨著計算機(jī)科學(xué)技術(shù)和數(shù)學(xué)地質(zhì)建模技術(shù)的發(fā)展，大量學(xué)者提出了多種基于專家經(jīng)驗和統(tǒng)計分析驅(qū)動的數(shù)據(jù)分析處理判別模型，實現(xiàn)了大部分地球科學(xué)數(shù)據(jù)成果的定性和半定量化分析研究與應(yīng)用普及［8-9］。特別是模糊數(shù)學(xué)理論的出現(xiàn)，使得地球物理異常數(shù)據(jù)的處理和定量解譯成為可能［10］。采空區(qū)地球物理異常作為礦致地球物理異常之上的二次人造異常，更具多樣性、變異性和復(fù)雜性，造成的不確定性具有模糊性，尤其是采空區(qū)地球物理異常的最終解譯成果只能反映異常圈定區(qū)內(nèi)有無采空區(qū)存在，而這類“有、“無”之間就是一種模糊關(guān)系。

為解決采空區(qū)地球物理異常圈定的不確定性和專家經(jīng)驗圈定的滯后性，便于及時開展鉆探驗證和安全處理工作，基于多年的BIF 型鐵礦山采空區(qū)探測成果，特別是采空區(qū)類型和地球物理特征研究成果的總結(jié)，應(yīng)用模糊數(shù)學(xué)理論和模型評價方法，基于GIS 平臺構(gòu)建采空區(qū)空間地球物理異常影像快速處理與定量評價模型，以期實現(xiàn)金屬礦山不明采空區(qū)地球物理異常圈定的模型化和智能化。

1 露天鐵礦采空區(qū)類型和地球物理識別特征

1.1 露天鐵礦采空區(qū)類型

采空區(qū)作為一種地表以下一定深度礦產(chǎn)資源開采后形成的空間，目前主要存在如下類型劃分思路［1］。

（1）基于資源開采時間的采空區(qū)分類。根據(jù)采空區(qū)的形成時間，采空區(qū)可以分為3類：第1類為歷史時期形成的老采空區(qū)，第2類為目前正在開采形成的現(xiàn)采空區(qū)，第3類為未來計劃開采區(qū)域形成的采空區(qū)。

（2）基于開采礦產(chǎn)資源類型的采空區(qū)分類。根據(jù)開采礦產(chǎn)資源類型的不同，采空區(qū)可以分為金屬礦山采空區(qū)和非金屬礦山采空區(qū)兩大類，其中金屬礦山采空區(qū)依據(jù)采礦方法又可分為空場法采空區(qū)、崩落法采空區(qū)和充填法采空區(qū)。

（3）基于采空區(qū)規(guī)模大小分類?？梢院唵畏譃楠毩⒌V山采空區(qū)和整體礦山采空區(qū)兩類。其中，獨立礦山采空區(qū)按體積可以分為0.1～1 萬m3、1～3 萬m3和3～10萬m33類，整體礦山采空區(qū)按體積規(guī)模可以分為5類，分別為＜50萬m3、50～100萬m3、100～500萬m3、500～1 000 萬m3和＞1 000 萬m3。此外，還可以劃分為小規(guī)模的巷道式采空區(qū)和大規(guī)模的采場式采空區(qū)。

（4）基于采空區(qū)埋藏深度分類。根據(jù)采空區(qū)埋藏深度，可以將采空區(qū)分為兩類：第1 類為淺部采空區(qū)，埋藏深度在60 m 以內(nèi)，是礦山開采和工程施工必須及時探明和安全處理的采空區(qū)；第2類為深埋采空區(qū)，埋藏深度大于60 m，需要進(jìn)行安全評價、監(jiān)測與處理的采空區(qū)。

（5）基于采空區(qū)賦存狀態(tài)分類。根據(jù)采空區(qū)的賦存狀態(tài)，既可以分為已知采空區(qū)和未知采空區(qū)，也可以分為已塌陷采空區(qū)和未塌陷采空區(qū)。未知采空區(qū)，即為不明采空區(qū)或隱伏采空區(qū)，是地球物理探測和礦山安全防治的重點對象。

（6）基于采空區(qū)充填介質(zhì)類型的物性分類。不論何種采空區(qū)，一旦形成，在未被人為充填和完全塌陷之前，總保持一定的形態(tài)。對于其空間內(nèi)的充填介質(zhì)，除了常見的圍巖和礦石碎塊外，對其物性影響最大的介質(zhì)即為四處流動的空氣和水。依據(jù)對采空區(qū)物性差異的影響程度，一般可以將不同充填介質(zhì)類型的采空區(qū)劃分為兩大類：第1類為空氣充填型采空區(qū)，其內(nèi)充填介質(zhì)以空氣為主；第2類為水充填型采空區(qū)，其內(nèi)充填介質(zhì)以水為主。該分類方式也是不明采空區(qū)可以進(jìn)行地球物理探測的前提條件和理論基礎(chǔ)。

1.2 采空區(qū)地球物理識別特征

不同類型的采空區(qū)具有不同的地球物理特征，而不同地球物理特征的采空區(qū)來自于不同物性差異對采空區(qū)類型及賦存圍巖類型的可區(qū)分程度，這種區(qū)分程度也是選擇不同地球物理探測方法的理論基礎(chǔ)。基于物性差異的認(rèn)識，本研究將采空區(qū)分為兩大類，即空氣充填型采空區(qū)和水充填型采空區(qū)。

（1）空氣充填型采空區(qū)?？諝獬涮钚筒煽諈^(qū)以一種極低密度絕緣體空氣為充填介質(zhì)。該類采空區(qū)與賦存空間的巖礦石相比，具有極低密度、高電阻率和低電導(dǎo)率的地球物理識別特征，可以采用高密度電阻率法、瞬變電磁法和地震映像法進(jìn)行組合探測。相關(guān)研究表明，空氣充填型采空區(qū)異常在高密度電阻率法和瞬變電磁法數(shù)據(jù)影像上均表現(xiàn)為高阻值異常圈閉空間，且阻值從異常邊界到異常核心逐漸增大，呈正向遞增梯度變化趨勢［11-12］；在對地震映像法探測數(shù)據(jù)進(jìn)行解譯時，可通過抽取最佳偏移距道數(shù)據(jù)來快速定位采空區(qū)異常位置，具體表現(xiàn)為地震波形圖同相軸的突然增減和錯斷，且主要以同相軸的錯斷增多為采空區(qū)異常存在特征［13-14］。

（2）水充填型采空區(qū)。水充填型采空區(qū)以一種來源于礦山水循環(huán)系統(tǒng)、富含多種礦物質(zhì)元素的水為充填介質(zhì)，對外顯示為低密度的良導(dǎo)體。該類采空區(qū)與賦存空間的巖礦石相比，具有低密度、低電阻率和高電導(dǎo)率的地球物理識別特征，也可以選用分辨率高的高密度電阻率法、瞬變電磁法和地震映像法進(jìn)行組合探測。相關(guān)研究表明，水充填型采空區(qū)異常在高密度電阻率法和瞬變電磁法數(shù)據(jù)影像上均表現(xiàn)為低阻值異常圈閉空間，且阻值從異常邊界到異常核心逐漸減小，呈負(fù)向遞減梯度變化趨勢［7，11］；對于地震映像法探測數(shù)據(jù)，水充填型采空區(qū)與空氣充填型采空區(qū)的地球物理異常識別特征一致［13-14］。

2 模糊數(shù)學(xué)理論評價方法

模糊數(shù)學(xué)評價方法作為一種綜合評價方法，主要對受到多個不確定性因素影響的目標(biāo)進(jìn)行整體性評價［15-16］，這就需要首先確定評價對象的影響因素。該類因素對評價結(jié)果的影響程度不同，而需要確定該類因素對上層評價目標(biāo)的權(quán)重，并經(jīng)過模糊變換得到該類因素對評價總目標(biāo)的隸屬度，最終完成對研究目標(biāo)的綜合評價［17］。假設(shè)2 個有限集W=，且W 與K 的模糊關(guān)系可表示為K=WY。其中，Y 是W 到K 的一個模糊變換；W 為影響因素權(quán)重集；K 為模糊綜合評價成果集，表示最終綜合評價結(jié)果；Y 為模糊判別矩陣，描述因素權(quán)重集W和評價成果集K之間的模糊關(guān)系。

本研究考慮到物探數(shù)據(jù)的異常不確定性和間接性，為最大限度減少人為干預(yù)對評價結(jié)果的影響，故采用全程模糊綜合評價模型。首先，依據(jù)前期研究成果選出采空區(qū)異常評價指標(biāo)因素集，通過數(shù)據(jù)重構(gòu)提取每個評價指標(biāo)因素，形成評價指標(biāo)數(shù)據(jù)影像集；然后，對于獲取的每個評價指標(biāo)數(shù)據(jù)進(jìn)行自然間斷分級和模糊賦權(quán)分類，形成評價指標(biāo)概率分布影像集；最后，通過模糊綜合評價模型判別采空區(qū)綜合地球物理異常，形成采空區(qū)異常分布概率影像圖。

3 露天鐵礦采空區(qū)異常評價指標(biāo)

根據(jù)大量的BIF 型鐵礦采空區(qū)地球物理探測研究和鉆探驗證成果經(jīng)驗總結(jié)［7，11-12，18］，優(yōu)化組合出集高密度電阻率法、瞬變電磁法和地震映像法于一體的采空區(qū)異常探測方法，并利用方法中采空區(qū)與周圍巖礦石物性差異最大的一組物性參數(shù)組合來評價采空區(qū)的存在與否，如電法中的電阻率和電磁法中的電導(dǎo)率等。此外，依據(jù)探測數(shù)據(jù)中的隱含數(shù)據(jù)信息和不明采空區(qū)異常的空間賦存規(guī)律，如空氣充填型采空區(qū)的高阻異常圈閉性和地表的高阻異常不可能為采空區(qū)異常等，建立了如圖1所示的采空區(qū)地球物理異常評價指標(biāo)體系。對于具體評價指標(biāo)的數(shù)據(jù)特征和數(shù)據(jù)重構(gòu)及其綜合評價，本研究以空氣充填型采空區(qū)為例進(jìn)行分析。

3.1 高密度電阻率法數(shù)據(jù)特征和數(shù)據(jù)重構(gòu)

利用高密度電阻率法探測采空區(qū)，主要是利用采空區(qū)與周圍巖礦石的電阻率差異，具體表現(xiàn)在電阻率數(shù)據(jù)二維影像上顯示為相對高阻值塊區(qū)和相對低阻值塊區(qū)（圖2（a）和圖2（b））。對于空氣充填型采空區(qū)，高密度電阻率法數(shù)據(jù)特征主要表現(xiàn)為高阻值異常圈閉區(qū)，從異常區(qū)的邊界到核心，阻值逐漸增高，表現(xiàn)為正向梯度變化。考慮到空氣充填型采空區(qū)具有正向梯度變化特征，可以采用離散傅里葉變換垂向求導(dǎo)獲得正向梯度值圈閉區(qū)域［19］（圖2（c）），即評價指標(biāo)體系中的采空區(qū)異常梯度矢量特征G2T。對于電阻率高低的評價指標(biāo)G1R，本研究按照傳統(tǒng)地球科學(xué)異常評價方法，在無遺漏采空區(qū)異常標(biāo)準(zhǔn)指導(dǎo)下，以每幅影像的電阻率平均值為背景值，高于背景值的區(qū)塊為正異?！咦柚等﹂]區(qū)→空氣充填型采空區(qū)異常，低于背景值的區(qū)塊為負(fù)異?！妥柚等﹂]區(qū)→水充填型采空區(qū)異常。此外，探測區(qū)域的地表接地條件、不明采空區(qū)的賦存位置和探測方法數(shù)據(jù)的空間分布等均可對采空區(qū)異常圈定的可靠性產(chǎn)生極大影響，如地表的正向梯度異常和高阻值圈閉難以成為采空區(qū)異常。為此，本研究提出了采空區(qū)異常埋深條件G3H作為評價指標(biāo)，評價電阻率法探測數(shù)據(jù)點的空間分布可信度。

上述3 個評價指標(biāo)G1R、G2T和G3H，由于數(shù)據(jù)格式和數(shù)據(jù)取值范圍存在很大差異，需要進(jìn)行標(biāo)準(zhǔn)化的數(shù)據(jù)重構(gòu)處理，主要是對各個評標(biāo)指標(biāo)涉及的參數(shù)值區(qū)間進(jìn)行分類賦值。為使各個類所涵蓋的數(shù)值區(qū)間差異最大化，并體現(xiàn)出固有的自然分組圈閉特征，本研究采用了自然間斷點10 級分類法，形成了高密度電阻率法的3 個采空區(qū)異常評價指標(biāo)影像分級圖（圖2（d）、圖2（e）和圖2（f））。

3.2 瞬變電磁法數(shù)據(jù)特征和數(shù)據(jù)重構(gòu)

瞬變電磁法通過測量地下介質(zhì)由于一次發(fā)射脈沖電磁場產(chǎn)生的二次感應(yīng)電磁場隨時間變化的衰減特征，來獲取地下不均勻體的導(dǎo)電性能和位置（圖3（a）），其反演結(jié)果經(jīng)過處理后如圖3（b）所示，但由于瞬變電磁法存在關(guān)斷時間效應(yīng)，其淺部為探測盲區(qū)?？諝獬涮钚筒煽諈^(qū)在瞬變電磁法成果圖中也表現(xiàn)為高阻值異常圈閉特征，從圈閉異常的邊緣到核心，電阻率逐漸增高，同樣可以采用離散傅里葉變換垂向求導(dǎo)獲得正向梯度值圈閉區(qū)域（圖3（c）），即評價指標(biāo)體系中的采空區(qū)異常梯度矢量特征S2T。對于電阻率評價指標(biāo)S1R，可采用探測成果電阻率的平均值為背景值，高于背景值的區(qū)塊即可判定為空氣充填型采空區(qū)異常。由于瞬變電磁法存在探測盲區(qū)，且隨著探測信號的衰減造成深部成果的可信度較低，故本研究提出采空區(qū)異常埋深條件S3H來綜合評價探測數(shù)據(jù)點的空間異?？尚哦?。

為使瞬變電磁法的探測成果可與高密度電阻率法的探測成果進(jìn)行模糊綜合評價，需對其3個評價指標(biāo)（S1R、S2T和S3H）進(jìn)行標(biāo)準(zhǔn)化的數(shù)據(jù)重構(gòu)處理。數(shù)據(jù)重構(gòu)處理同樣采用自然間斷點10 級分類法，形成了瞬變電磁法的3 個評價指標(biāo)影像分級圖（圖3（d）、圖3（e）和圖3（f））。

3.3 地震映像法數(shù)據(jù)特征和數(shù)據(jù)重構(gòu)

基于最佳偏移距的地震映像法主要用于快速確定采空區(qū)異常存在與否和空間位置，其埋深和規(guī)模可通過與高密度電阻率法和瞬變電磁法綜合評價獲得。采空區(qū)在地震映像法中主要表現(xiàn)為同相軸的錯斷和增多（圖4（a）），個別顯示為同相軸的錯斷和減少。為此，本研究將空氣充填型采空區(qū)在地震映像法中的數(shù)據(jù)異常特征分為3 個類別進(jìn)行位置區(qū)間賦值重構(gòu)，無采空區(qū)異常位置區(qū)段賦值為1，采空區(qū)異常顯示一般區(qū)段賦值為3，采空區(qū)異常顯示強區(qū)段賦值為5，得到地震映像法的采空區(qū)異常評價指標(biāo)D1W（圖4（b））。

4 基于模糊數(shù)學(xué)模型的露天鐵礦采空區(qū)探測地球物理異常評價

依據(jù)高密度電阻率法、瞬變電磁法和地震映像法的探測成果，提取了7個采空區(qū)地球物理異常評價指標(biāo)，分別為G1R、G2T、G3H、S1R、S2T、S3H和D1W，且當(dāng)它們每個指標(biāo)在一個指定區(qū)段內(nèi)都達(dá)到權(quán)重概率預(yù)定值以上，方可綜合判定為采空區(qū)異常。同時，需在保證無遺漏采空區(qū)的條件下縮小探測數(shù)據(jù)異常圈定范圍，以降低鉆探驗證成本。因此，有必要對7 個不確定性評價指標(biāo)合理提取關(guān)鍵參數(shù)并進(jìn)行賦權(quán)，確保評價模型的可靠性和拓展性，從而可適用于復(fù)雜多變的礦山地球物理探測評價，達(dá)到采空區(qū)地球物理異常評價的定量化和智能化。本研究采用已驗證區(qū)域的空氣充填型采空區(qū)的震—電—磁探測數(shù)據(jù)來開展研究（圖5），同時，為了確保一致性和系統(tǒng)性，文中用到的震—電—磁地球物理探測成果數(shù)據(jù)為同一組數(shù)據(jù)。

4.1 評價指標(biāo)模糊分類

7 個采空區(qū)異常評價指標(biāo)可以分為4 類：第1 類為探測數(shù)據(jù)值指標(biāo)，為G1R和S1R；第2類為探測數(shù)據(jù)提取指標(biāo)，為G2T和S2T；第3 類為有效探測區(qū)間指標(biāo)，為G3H和S3H；第4類為采空區(qū)探測異常位置指標(biāo)，為D1W。對于上述每類評價指標(biāo)的模糊分類，可使已進(jìn)行數(shù)據(jù)重構(gòu)的評價指標(biāo)歸一化和概率化，即將采空區(qū)異常評價指標(biāo)分類數(shù)據(jù)轉(zhuǎn)換到0～1 的概率區(qū)間內(nèi)。該區(qū)間內(nèi)，“0”分配給確定不屬于指定采空區(qū)異常集合的分類值，“1”分配給確定屬于指定采空區(qū)異常集合的分類值，0～1 區(qū)間內(nèi)的概率值分配給可能為采空區(qū)異常的分類值。

4.1.1 第一類采空區(qū)異常評價指標(biāo)模糊分類

G1R代表高密度電阻率法的直接電阻率測量值，對于空氣充填型采空區(qū)，在無遺漏采空區(qū)異常標(biāo)準(zhǔn)指導(dǎo)下，探測的電阻率值若大于背景值或探測影像的電阻率平均值，即可認(rèn)定為采空區(qū)異常，可以采用模糊較大值計算模型進(jìn)行分析：

式中，x為輸入計算值；f1為定義過渡區(qū)形狀和特征的展開參數(shù)，本研究取3；f2為定義的起算中點，用于確定交叉點（起算中點分配概率值為0.5，大于起算中點的值為采空區(qū)異常的概率較高，小于起算中點的值為采空區(qū)異常的概率較低）。

該模型以評價指標(biāo)G1R中探測影像電阻率平均值所代表的分類值（圖2（e））為起算中點f2，其起算概率為0.5，將柵格數(shù)據(jù)G1R（圖2（e））圖像值作為輸入值進(jìn)行模糊計算，得到采空區(qū)異常評價指標(biāo)G1R指示采空區(qū)異常存在的概率分布圖（圖6（a）），其概率值越大，屬于采空區(qū)的可能性便越大。同理，可以采用模糊較大值計算模型對瞬變電磁法采空區(qū)異常評價指標(biāo)S1R（圖3（e））進(jìn)行采空區(qū)異常概率分布計算，其起算中點f2為探測影像電阻率均值所代表的分類值，可以得到采空區(qū)異常評價指標(biāo)S1R指示采空區(qū)異常存在的概率分布圖（圖6（b））。

4.1.2 第2類采空區(qū)異常評價指標(biāo)模糊分類

第2 類采空區(qū)異常評價指標(biāo)是高密度電阻率法的G2T和瞬變電磁法的S2T。G2T為電阻率的空間變化梯度值，正梯度值為高阻值異常圈閉區(qū)顯示，代表空氣充填型采空區(qū)異常。該評價指標(biāo)可以采用模糊高斯模型進(jìn)行計算：

式中，x為輸入計算值；f1為定義過渡區(qū)形狀和特征的展開參數(shù)，本研究取3；f2為定義的起算中點，用于確定交叉點（起算中點分配概率值為1，大于起算中點的值或小于起算中點的值對應(yīng)為采空區(qū)異常的概率均較低）。

模糊高斯模型以指標(biāo)G2T平均值所代表的分類值為起算中點f2，起算概率為1，小于或大于起算中點的概率均較低，該模型可以有效消除地表接觸不良引起的正梯度值超大“采空區(qū)異?！?，也可以規(guī)避將正梯度值過低的采空區(qū)異常邊界過渡區(qū)確定為采空區(qū)異常。本研究將柵格數(shù)據(jù)G2T（圖2（d））圖像值作為輸入計算值進(jìn)行模糊計算，得到了采空區(qū)異常評價指標(biāo)G2T指示采空區(qū)異常存在的概率分布圖（圖6（c））。

對于瞬變電磁法的采空區(qū)異常評價指標(biāo)S2T，由于存在關(guān)斷時間導(dǎo)致的淺部探測盲區(qū)，可有效消除因地表接觸不良造成的“高阻采空區(qū)異常圈閉”，其正梯度值（圖3（d））變化越大，顯示為采空區(qū)異常的概率便越大，也可采用模糊較大值計算模型進(jìn)行計算，其起算中點f2為正梯度值平均值所代表的分類值，計算得到的采空區(qū)異常評價指標(biāo)S2T指示采空區(qū)異常存在的概率分布如圖6（d）所示。

4.1.3 第3類采空區(qū)異常評價指標(biāo)模糊分類

第3 類采空區(qū)異常評價指標(biāo)為有效探測區(qū)間指標(biāo)，為G3H和S3H（圖2（f）和圖3（f））。該類評價指標(biāo)可以有效篩除淺部接觸不良導(dǎo)致的異常數(shù)據(jù)、可信度低的深邊部探測數(shù)據(jù)和關(guān)斷時間導(dǎo)致的探測盲區(qū)數(shù)據(jù)，以降低采空區(qū)異常誤判概率。該類指標(biāo)可采用模糊高斯模型（式（2））進(jìn)行計算，并依據(jù)前期探測和驗證成果，將有效探測數(shù)據(jù)深度區(qū)間的中心深度值所代表的分類值作為起算中點f2，賦權(quán)概率為1，而淺部或深部異常數(shù)據(jù)為采空區(qū)的概率均較低，得到G3H和S3H的采空區(qū)異常概率分布圖（圖6（e）和6（f））。

4.1.4 第4類采空區(qū)異常評價指標(biāo)模糊分類

地震映像法的采空區(qū)異常位置指標(biāo)D1W，可依據(jù)其異常分類特征采用模糊較大值模型（式（1））計算，f1取值為2，其起算中點f2為3 級異常分類中的1 級和2 級采空區(qū)異常的中間值2（圖4（b）），起算概率為0.5，圖像計算得到的采空區(qū)異常評價指標(biāo)D1W的采空區(qū)異常概率分布圖如圖6（g）所示。

4.2 評價模型構(gòu)建

本研究采用數(shù)據(jù)重構(gòu)方法和模糊數(shù)學(xué)計算模型提取了7 個采空區(qū)異常評價指標(biāo)G1R、G2T、G3H、S1R、S2T、S3H和D1W，綜合評價采空區(qū)異常存在與否。鑒于上述每個采空區(qū)異常評價指標(biāo)均來自不確定性判別數(shù)據(jù)集或依據(jù)專家經(jīng)驗總結(jié)得出，且評價指標(biāo)模糊分類中的概率賦值也選用了最低可能性，以免遺漏采空區(qū)異常，具體表現(xiàn)在有效探測數(shù)據(jù)區(qū)間內(nèi)，空區(qū)存在的區(qū)域必有采空區(qū)異常，而有采空區(qū)異常的區(qū)域未必全是空區(qū)。因此，考慮到指標(biāo)體系的不確定性，本研究空區(qū)異常評價采用模糊綜合評價模型（圖6（h））［20］，即模糊疊加模型，具體可為模糊And 函數(shù)評價模型：

式中，arg1,arg2,arg3,…,argn 為不同的評價指標(biāo)取值，本研究n ∈[ ]1,7，F(xiàn)uzzy And 疊加類型將返回單元位置所屬集合的最小值。

模糊And疊加模型會創(chuàng)建一個輸出柵格，為像元位置所屬的每個集合的每個像元值給定了最小模糊值。本研究中，采空區(qū)異常評價模型是7個指標(biāo)條件中的每個條件都已相對于其屬于合適集合的隸屬度被模糊化，模糊And模型會識別多個條件內(nèi)像元屬于一個合適集合的最小可能性。對于采空區(qū)地球物理異常評價，主要依據(jù)無遺漏采空區(qū)異常標(biāo)準(zhǔn)和模糊分類的概率賦值原則，只需要選擇符合所有指標(biāo)中概率值大于或等于0.5 的具有采空區(qū)異常的空間位置集合，疊加生成采空區(qū)異常指標(biāo)綜合評價概率分布圖，圖中概率值大于0.5 的區(qū)塊圈閉即為采空區(qū)異常，需要進(jìn)一步進(jìn)行鉆探驗證。

4.3 評價結(jié)果

基于上述綜合評價結(jié)果，采用ArcGIS 軟件中的空間分析模塊進(jìn)行弓長嶺鐵礦區(qū)多采空區(qū)異常評價指標(biāo)圖像矩陣運算，計算模型為Fuzzy And，結(jié)果見圖7。圖7中測線起點以高密度電阻率法G8測線起點為0 點，與其平行的瞬變電磁法和地震映像法測線則以相對于電法起點疊加距離作為起點計算空間坐標(biāo)。在圖7 中，根據(jù)0.5 及其以上概率值作為異常圈定標(biāo)準(zhǔn)，圈定出了一個采空區(qū)異常區(qū)（即空氣充填型采空區(qū)）。在對異常區(qū)進(jìn)行驗證時，可將采空區(qū)異常區(qū)上部邊界與驗證得出的采空區(qū)頂板埋深進(jìn)行對比分析。圖7 中采空區(qū)異常區(qū)圈定的頂板埋深為35.7 m，經(jīng)過鉆探驗證埋深為33.8 m，兩者相差1.9 m，相比于傳統(tǒng)經(jīng)驗法的圈定誤差（±5.0 m 以內(nèi)），精度提高了3.0 m以上。

通過對空氣充填型采空區(qū)的地球物理探測數(shù)據(jù)進(jìn)行指標(biāo)化提取，并采用本研究提出的綜合評價模型，在ArcGIS 平臺上定量圈定了采空區(qū)異常邊界，通過鉆探成果相比，可以在無遺漏采空區(qū)的前提下，大大提高采空區(qū)圈定效率和準(zhǔn)確度。在上述評價方法指導(dǎo)下，對弓長嶺鐵礦區(qū)水充填型采空區(qū)進(jìn)行了綜合評價，圈定的采空區(qū)異常頂板埋深18.3 m，鉆探驗證結(jié)果為20.1 m，兩者相差在2 m 以內(nèi)?？梢?，本研究選擇的采空區(qū)異常圈定指標(biāo)具有合理性，即對于不確定性的地球物理異常可以采用模糊數(shù)學(xué)方法進(jìn)行定量研究，實現(xiàn)地球物理探測數(shù)據(jù)評價的指標(biāo)化、定量化和智能化。此外，本研究評價方法可以通過ArcGIS 平臺對采集到的大量探測數(shù)據(jù)進(jìn)行流程式的智能化快速綜合評價，并可依據(jù)評價結(jié)果進(jìn)行快速鉆探驗證，確保礦山生產(chǎn)安全。

5 結(jié) 論

（1）劃分了多種采空區(qū)類型，依據(jù)采空區(qū)充填介質(zhì)和物性差異的不同，將采空區(qū)分為空氣充填型采空區(qū)和水充填型采空區(qū)兩大類，并分別建立了相應(yīng)的地球物理識別特征和震—電—磁探測方法組合。

（2）根據(jù)高密度電阻率法、瞬變電磁法和地震映像法的采空區(qū)異常數(shù)據(jù)識別特征，構(gòu)建并提取了7個采空區(qū)異常評價指標(biāo)G1R、G2T、G3H、S1R、S2T、S3H和D1W。

（3）通過對7個采空區(qū)異常指標(biāo)進(jìn)行模糊分類和模糊數(shù)學(xué)模型綜合評價能夠定量圈定出采空區(qū)異常邊界，可使探測精度保持在±2.0 m以內(nèi)，為采空區(qū)異?？焖俣炕椭悄芑u價提供了技術(shù)手段。