基于證據(jù)權(quán)模型的赤峰—呼倫貝爾成礦帶找礦預(yù)測

田九玲

（河南工業(yè)職業(yè)技術(shù)學(xué)院城市建設(shè)學(xué)院，河南南陽473000）

赤峰—呼倫貝爾成礦帶是大興安嶺南部黃崗梁—甘珠爾廟成礦帶的一部分，該區(qū)域地質(zhì)構(gòu)造復(fù)雜，地殼以及巖漿巖活動劇烈，在大興安嶺南部區(qū)域中形成多金屬礦產(chǎn)的潛力最大［1-2］。近年來，不少學(xué)者對該成礦區(qū)的礦床成因、成礦規(guī)律等進行了深入研究［3-6］，但對于該區(qū)域成礦預(yù)測方面的研究涉及較少。本研究結(jié)合赤峰—呼倫貝爾成礦帶的地質(zhì)特征，利用地理信息系統(tǒng)空間分析軟件（ArcGIS）對獲取的研究區(qū)域地層信息、構(gòu)造信息、巖漿巖信息進行空間分析，形成空間數(shù)據(jù)庫，在此基礎(chǔ)上引入引入證據(jù)權(quán)模型，進行成礦預(yù)測分析，圈定找礦靶區(qū)，為該區(qū)后續(xù)找礦工作提供有益參考。

1 研究區(qū)域地質(zhì)特征

研究區(qū)域位于內(nèi)蒙古地區(qū)，具體位置如圖1所示。研究區(qū)氣候特點為冬季冷且長，夏季短且熱，整個區(qū)域處于半干旱季風(fēng)氣候區(qū)，有水系面積達97 000 km2，陸地和草地豐富，居民多以畜牧業(yè)、礦產(chǎn)資源開采等為業(yè)［7-9］。

1.1 地 層

研究區(qū)位于華北地區(qū)，屬于華北地層區(qū)，主要有中生代地層和二疊紀(jì)地層出露，上部蓋層主要由新生界中第四系沉積巖、中生界中侏羅系沉積巖和古生界中的二疊系沉積巖組成［10-12］，如表1所示。研究區(qū)中古生界二疊系中的大石寨組、林西組等含有豐富的礦產(chǎn)資源，該部分為主要的礦產(chǎn)賦存地層。

1.2 構(gòu) 造

研究區(qū)受地質(zhì)作用影響強烈，主要有褶皺構(gòu)造、斷裂構(gòu)造。褶皺構(gòu)造主要位于礦帶復(fù)背斜部分，長約400 km［13-14］。斷裂構(gòu)造中與主要褶皺構(gòu)造平行的斷裂帶主要控制該區(qū)域的礦產(chǎn)資源，區(qū)內(nèi)構(gòu)造對于礦產(chǎn)資源成礦預(yù)測具有一定的作用。

1.3 巖漿巖

區(qū)域內(nèi)的巖漿巖形成時代不同，分別于華力西期、印支期、燕山期。其中，華力西期形成的巖漿巖分布面積最大，為酸性侵入巖，主要有花崗巖、白云母等巖石；印支期形成的侵入巖分布較少，主要為二長花崗巖等巖石；燕山期形成的巖石分布較廣，主要有石英等巖石。

2 找礦預(yù)測模型構(gòu)建

2.1 基于證據(jù)權(quán)模型的礦區(qū)預(yù)測思路

地理信息系統(tǒng)（Geographic Information System，GIS）自20世紀(jì)80年代應(yīng)用于礦產(chǎn)行業(yè)后，發(fā)展迅速，數(shù)據(jù)權(quán)模型是在地理信息系統(tǒng)預(yù)測模型的基礎(chǔ)上結(jié)合概率學(xué)、數(shù)理統(tǒng)計等學(xué)科發(fā)展而來［15］。本研究基于證據(jù)權(quán)理論的礦區(qū)找礦預(yù)測思路見圖2。找礦預(yù)測步驟如下：

（1）建立數(shù)據(jù)庫。通過地理信息系統(tǒng)技術(shù)收集赤峰—呼倫貝爾成礦帶各區(qū)域的地質(zhì)和地理數(shù)據(jù)等數(shù)據(jù)信息，利用ArcGIS軟件對獲取的研究區(qū)地層、構(gòu)造、巖漿巖信息進行空間分析，形成空間數(shù)據(jù)庫。

（2）成礦預(yù)測。通過對空間數(shù)據(jù)庫中的數(shù)據(jù)進行處理，確定各區(qū)域的權(quán)重和概率，進而根據(jù)概率計算結(jié)果劃分具體的成礦靶區(qū)。

2.2 赤峰—呼倫貝爾成礦帶空間數(shù)據(jù)庫建立

空間數(shù)據(jù)庫主要通過地理信息系統(tǒng)技術(shù)獲取區(qū)域內(nèi)的地層、構(gòu)造、巖漿巖信息，通過對數(shù)據(jù)進行處理形成可供預(yù)測模型利用的數(shù)據(jù)庫。

通過對區(qū)內(nèi)地層信息進行處理，得到各區(qū)域中礦點數(shù)及其面積。寶音圖群面積為446.6 km2，礦點1個；包爾漢圖群的面積為77.1 km2，礦點1個；本巴圖組面積為1 003.9 km2，礦點1個；寶力高廟組的面積為1 190.5 km2，礦點1個；大石寨組面積為4 512.6 km2，礦點26個；于家北溝組面積為804.4 km2，礦點7個；壽山溝組面積為1 006.3 km2，礦點12個；額里圖組面積為416.9 km2，礦點2個；林西組面積為4 300.8 km2，礦點7個；哲思組面積為3 181.7 km2，礦點11個；滿克頭鄂博組面積為15 541.8 km2，礦點26個；白音高老組面積為8 679.1 km2，礦點7個；瑪尼組面積為4 386.6 km2，礦點7個；新民組面積為2 605.6 km2，礦點8個。

通過對獲得的構(gòu)造信息進行處理，發(fā)現(xiàn)區(qū)內(nèi)線性構(gòu)造較多，該類構(gòu)造內(nèi)礦點數(shù)為75個，約占全區(qū)礦點總數(shù)的12。

通過對巖漿巖信息進行處理，得到該區(qū)域內(nèi)巖漿巖的種類以及含有的礦點數(shù)量為：石炭紀(jì)多為閃長巖，礦點1個；二疊紀(jì)多為花崗巖、閃長巖等，礦點3個；三疊紀(jì)多為花崗閃長巖等，礦點4個；侏羅紀(jì)多為正長斑巖等，礦點14個；白堊紀(jì)多為花崗巖，礦點3個。

2.3 證據(jù)權(quán)預(yù)測模型建立

證據(jù)權(quán)模型是在證據(jù)權(quán)重數(shù)理統(tǒng)計方法的基礎(chǔ)上進行構(gòu)建，該方法是由Agterberg（地質(zhì)學(xué)家）提出［16-17］。研究區(qū)是由T個等面積的部分組成，含有礦產(chǎn)資源的部分個數(shù)為D，由GIS技術(shù)提取的礦區(qū)信息作為影響因子即證據(jù)因子，后驗概率O可由下式解得：

式中，j為證據(jù)因子數(shù)量，j=1,2,3,…,n；第j個證據(jù)因子的權(quán)重可進行如下定義：

（1）當(dāng)證據(jù)因子存在時，

（2）當(dāng)證據(jù)因子不存在時，

（3）當(dāng)沒有數(shù)據(jù)或者數(shù)據(jù)不全時，

于是，全區(qū)成礦概率模型為

根據(jù)式（5）計算得到的概率大小可判斷相應(yīng)區(qū)域的成礦潛力。

3 赤峰—呼倫貝爾成礦帶找礦預(yù)測

3.1 找礦預(yù)測

通過將空間數(shù)據(jù)庫中的信息進行處理，提取為地質(zhì)信息、地球化學(xué)信息、地球物理信息以及其他信息，并將這些信息代入證據(jù)權(quán)預(yù)測模型中，得到該區(qū)域的先驗概率如表2所示。在礦點出現(xiàn)的情況下，證據(jù)因子存在、不存在的概率分別為O1、O2；在礦點未出現(xiàn)的情況下，證據(jù)因子存在、不存在的概率分別為O3、O4。

根據(jù)表2計算出該證據(jù)因子的權(quán)重如表3所示。由表3可知：在地層信息因子中，大石寨組權(quán)重最大，哲斯組權(quán)重次之；在構(gòu)造信息因子中，重力異常權(quán)重最大，構(gòu)造等密度異常權(quán)重次之；在多金屬元素化探信息因子中，化探Ag-Pb-Zn異常權(quán)重最大；在單金屬元素化探信息因子中，化探W異常權(quán)重最大。根據(jù)上述權(quán)重計算結(jié)果，可對各個區(qū)域的成礦潛力進行預(yù)測。

3.2 找礦靶區(qū)圈定

首先根據(jù)證據(jù)權(quán)模型預(yù)測結(jié)果得到研究區(qū)礦產(chǎn)資源分布情況，利用中國地質(zhì)調(diào)查局MRAS系統(tǒng)將表2和表3中的數(shù)據(jù)進行處理，得到赤峰—呼倫貝爾成礦帶的找礦預(yù)測結(jié)果如圖3所示［18］，并按照成礦潛力由大到小將該區(qū)域分為3個等級，分別為Ⅰ級、Ⅱ級、Ⅲ級；然后分析各個靶區(qū)的成礦潛力；最后進行現(xiàn)場工程驗證，將目前該區(qū)域內(nèi)已經(jīng)發(fā)現(xiàn)的礦產(chǎn)資源與預(yù)測結(jié)果進行對比，驗證模型預(yù)測的準(zhǔn)確性。

Ⅰ級預(yù)測區(qū)內(nèi)礦產(chǎn)資源豐富、成礦量較大，鉛鋅等多金屬的成礦潛力非常巨大。根據(jù)MRAS概率計算結(jié)果，在礦帶上發(fā)現(xiàn)了兩處Ⅰ級預(yù)測區(qū)，分別為Ⅰ-1號靶區(qū)和Ⅰ-2號靶區(qū)。Ⅰ-1號靶區(qū)地理坐標(biāo)為東經(jīng)118°44′00″～119°48′02″，北緯43°36′53″～48°48′22″，該區(qū)域內(nèi)地層主要為大石寨組和林西組，Ag-Pb-Zn探測和W的探測結(jié)果變化量明顯異常，正磁異常表明該區(qū)域巖漿巖發(fā)育較好，根據(jù)預(yù)測，區(qū)內(nèi)鉛鋅成礦潛力巨大。目前該區(qū)域內(nèi)發(fā)現(xiàn)的有碧流臺等多金屬礦區(qū)［19］。Ⅰ-2號靶區(qū)地理坐標(biāo)為東經(jīng)117°01′00″～118°01′02″，北緯42°16′53″～43°45′15″，該區(qū)域內(nèi)地層主要為林西組和滿克頭鄂博組，Sn探測結(jié)果變化量明顯異常，表明該區(qū)錫多金屬礦產(chǎn)成礦潛力較大。東經(jīng) 118°44′00″～119°48′02″，北緯 43°36′53″～48°48′22″區(qū)域內(nèi)發(fā)現(xiàn)有黃崗錫鐵礦等多個著名礦區(qū)。

Ⅱ級預(yù)測區(qū)內(nèi)礦產(chǎn)資源豐富，成礦潛力較大。根據(jù)MRAS概率計算結(jié)果，在礦帶上預(yù)測出符合該條件的區(qū)域如圖3所示。靶區(qū)內(nèi)地層以大石寨組為主，Ag-Pb-Zn的探測結(jié)果明顯異常，表明該區(qū)域內(nèi)多金屬礦區(qū)的成礦潛力較大，且成礦類型中鉛鋅多金屬礦產(chǎn)資源較多，目前已經(jīng)發(fā)現(xiàn)有白音諾爾鉛鋅礦等多金屬礦區(qū)［20］。

Ⅲ級預(yù)測區(qū)內(nèi)有一定的成礦潛力，但還需要實際工程進行驗證。本研究在研究區(qū)圈定了兩處靶區(qū)，分別為Ⅲ-1、Ⅲ-2號靶區(qū)，該類區(qū)域有必要進一步部署工程進行找礦潛力驗證分析。

4 結(jié) 語

將地理信息系統(tǒng)技術(shù)和數(shù)據(jù)權(quán)理論相結(jié)合，構(gòu)建了基于證據(jù)權(quán)模型的赤峰—呼倫貝爾成礦帶找礦預(yù)測模型。該模型首先利用ArcGIS軟件對獲取的區(qū)內(nèi)地層、構(gòu)造、巖漿巖信息進行處理形成空間數(shù)據(jù)庫，然后結(jié)合證據(jù)權(quán)理論構(gòu)建找礦模型。通過將空間數(shù)據(jù)庫中的數(shù)據(jù)代入預(yù)測模型中得到區(qū)域內(nèi)各證據(jù)因子的權(quán)重以及概率，將區(qū)域劃分為3類靶區(qū)，并對該區(qū)找礦潛力進行了分析，對于在赤峰—呼倫貝爾成礦帶進一步開展找礦預(yù)測理論研究與實踐工作有一定的參考價值。