田九玲

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

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

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

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

1.1 地 層

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

1.2 構(gòu) 造

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

1.3 巖漿巖

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

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

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

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

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

（2）成礦預測。通過對空間數(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ù)進行處理形成可供預測模型利用的數(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ù)量為：石炭紀多為閃長巖，礦點1個；二疊紀多為花崗巖、閃長巖等，礦點3個；三疊紀多為花崗閃長巖等，礦點4個；侏羅紀多為正長斑巖等，礦點14個；白堊紀多為花崗巖，礦點3個。

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

證據(jù)權(quán)模型是在證據(jù)權(quán)重數(shù)理統(tǒng)計方法的基礎上進行構(gòu)建，該方法是由Agterberg（地質(zhì)學家）提出［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）當證據(jù)因子存在時，

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

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

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

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

3 赤峰—呼倫貝爾成礦帶找礦預測

3.1 找礦預測

通過將空間數(shù)據(jù)庫中的信息進行處理，提取為地質(zhì)信息、地球化學信息、地球物理信息以及其他信息，并將這些信息代入證據(jù)權(quán)預測模型中，得到該區(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ū)域的成礦潛力進行預測。

3.2 找礦靶區(qū)圈定

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

Ⅰ級預測區(qū)內(nèi)礦產(chǎn)資源豐富、成礦量較大，鉛鋅等多金屬的成礦潛力非常巨大。根據(jù)MRAS概率計算結(jié)果，在礦帶上發(fā)現(xiàn)了兩處Ⅰ級預測區(qū)，分別為Ⅰ-1號靶區(qū)和Ⅰ-2號靶區(qū)。Ⅰ-1號靶區(qū)地理坐標為東經(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ù)預測，區(qū)內(nèi)鉛鋅成礦潛力巨大。目前該區(qū)域內(nèi)發(fā)現(xiàn)的有碧流臺等多金屬礦區(qū)［19］。Ⅰ-2號靶區(qū)地理坐標為東經(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ū)。

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

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

4 結(jié) 語

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