內(nèi)蒙古新達來草原覆蓋區(qū)鉬銅多金屬成礦預測

熊雙才， 夏慶霖， 姚春亮， 張曉軍， 覃　瑩， 張　坤

(1.中國地質(zhì)大學(武漢)資源學院，湖北武漢430074； 2.新疆地質(zhì)礦產(chǎn)局第一地質(zhì)大隊，新疆昌吉831100； 3.緊缺礦產(chǎn)資源勘查協(xié)同創(chuàng)新中心，湖北武漢430074； 4.山東省第六地質(zhì)礦產(chǎn)勘查院，山東威海264209)

內(nèi)蒙古新達來草原覆蓋區(qū)鉬銅多金屬成礦預測

熊雙才1,2， 夏慶霖1,3， 姚春亮1， 張曉軍1， 覃瑩1， 張坤4

(1.中國地質(zhì)大學(武漢)資源學院，湖北武漢430074； 2.新疆地質(zhì)礦產(chǎn)局第一地質(zhì)大隊，新疆昌吉831100； 3.緊缺礦產(chǎn)資源勘查協(xié)同創(chuàng)新中心，湖北武漢430074； 4.山東省第六地質(zhì)礦產(chǎn)勘查院，山東威海264209)

摘要：內(nèi)蒙古新達來草原覆蓋區(qū)位于西伯利亞板塊東南大陸邊緣晚古生代陸緣增生帶，屬古亞洲成礦域的二連—東烏旗銅鉬鐵晚古生代—中生代成礦帶的一部分，找礦潛力巨大。該區(qū)被牧草和第四系大面積覆蓋，基巖出露較差，找礦工作相對滯后。通過研究該地區(qū)區(qū)域成礦地質(zhì)背景和典型礦床地質(zhì)特征，分析控礦因素與找礦標志，構建綜合信息預測模型，提取地質(zhì)、地球化學、地球物理、遙感等找礦信息，借助GIS技術進行覆蓋層“揭蓋”和對找礦信息進行提取，以GeoDAS 礦產(chǎn)預測軟件為平臺，運用多重分形濾波技術區(qū)分復雜背景與異常，提取礦化弱異常，應用模糊證據(jù)權法進行信息綜合集成，形成預測后驗概率圖，圈定找礦遠景區(qū)并進行驗證，探索內(nèi)蒙古草原覆蓋區(qū)內(nèi)生金屬礦床勘查與預測評價新途徑。

關鍵詞：草原覆蓋區(qū)；GeoDAS軟件；非線性模型；模糊證據(jù)權方法；成礦預測；內(nèi)蒙古新達來

doi:10.3969/j.issn.1674-3636.2015.03.422

中圖分類號：P612；P618.4

文獻標識碼：A

文章編號：1674-3636(2015)03-0422-09

收稿日期：2015-05-18；修回日期：2015-05-26；編輯：蔣艷

基金項目：中國地質(zhì)調(diào)查局項目“覆蓋區(qū)找礦靶區(qū)優(yōu)選與驗證”(12120113089300)

作者簡介：熊雙才(1987—)，男，助理工程師，碩士，主要從事成礦規(guī)律與成礦預測研究，E-mail：scXiong@163.com

0引言

內(nèi)蒙古擁有豐富的礦產(chǎn)資源，但大部分地區(qū)為牧草和新生界覆蓋區(qū)，基巖出露較差，找礦工作相對滯后。因此，選擇新達來草原覆蓋區(qū)開展典型示范研究，運用非線性礦產(chǎn)預測理論和方法深層發(fā)掘已有資料，探索解決因覆蓋層影響所造成的找礦信息“弱緩、復合與迭加、缺失與不完整”等難題的途徑，開展綜合信息礦產(chǎn)預測。

1區(qū)域地質(zhì)背景

新達來草原覆蓋區(qū)位于內(nèi)蒙古自治區(qū)錫林郭勒盟蘇尼特左旗，北與蒙古國接壤。研究區(qū)位于西伯利亞板塊東南大陸邊緣晚古生代陸緣增生帶，屬古亞洲成礦域的額濟納旗—興安嶺成礦區(qū)、二連—東烏旗銅鉬鐵晚古生代—中生代成礦帶(邵積東等，2008)。具有古亞洲洋演化與瀕太平洋造山帶疊加的特點(任紀舜等，1992)，古生代—中生代期間，發(fā)生了明顯的碰撞造山過程、碰撞后伸展作用及其間的構造大轉(zhuǎn)換階段(毛景文等，2005；Xiao et al,2003)。因此，該區(qū)大面積分布古生界和中生界，海西期、印支期和燕山期巖漿巖發(fā)育，各種構造形跡復雜。自2006年以來，在內(nèi)蒙古中北部地區(qū)紅格爾礦田及鄰區(qū)陸續(xù)發(fā)現(xiàn)了一系列鉬礦床(點)，如烏蘭德勒(大型)、準蘇吉花(中型)、達萊敖包(中型)、烏花敖包(大型)、寶格達烏拉(大型)、烏日尼圖(大型)，使得該區(qū)可能成為中國繼東秦嶺鉬礦帶、燕遼鉬礦帶及西拉木倫鉬礦帶之后的又一條重要鉬礦帶(劉翼飛等，2012)。

1.1　地層

地層區(qū)劃上古生界屬興安地層區(qū)東烏—呼瑪?shù)貙臃謪^(qū)，中新生界屬大興安嶺—燕山地層分區(qū)博克圖—二連浩特地層小區(qū)。研究區(qū)出露的地層有古生界奧陶系中統(tǒng)巴彥呼舒組(O2b)、泥盆系中下統(tǒng)泥鰍河組(D1-2n)、上石炭—下二疊統(tǒng)寶力高廟組(C2-P1bl)、上二疊統(tǒng)(P3)；中、新生界上侏羅統(tǒng)瑪尼吐組(J3mn)、上侏羅統(tǒng)白音高老組(J3b)、第四系更新統(tǒng)(Qp)、全新統(tǒng)(Qh)。各時代地層的展布方向為NE向，特別是古生代地層，由南到北具有由老變新的變化規(guī)律，反映了地層的分布受區(qū)域構造作用的控制，也說明晚古生代陸緣增生由北而南的遷移規(guī)律(圖1)。

圖1　新達來區(qū)域地質(zhì)圖(據(jù)陶繼雄等，2009a修改)Ⅰ-華北板塊北緣斷裂；Ⅱ-西拉木倫河斷裂；Ⅲ-二連—賀根山斷裂；Ⅳ-烏努爾—鄂倫春斷裂；Ⅴ-德爾布干斷裂；Ⅵ-額爾古納斷裂；Ⅶ-大興安嶺主脊斷裂1-洪沖積砂礫；2-半固結礫石層；3-伊爾丁組；4-白音高老組；5-寶力高廟組；6-泥鰍河組；7-巴彥呼舒組；8-二疊紀堿長花崗巖；9-二疊紀二長花崗巖；10-二疊紀黑云母二長花崗巖；11-二疊紀黑云母花崗巖；12-石炭紀黑云母花崗巖；13-石炭紀黑云母二長花崗巖；14-背斜及編號；15-向斜及編號；16-逆斷層及編號；17-平移斷層及編號；18-實測性質(zhì)不明斷層及編號；19-推測斷層及編號；20-不整合線；21-侵入界線；22-韌性剪切帶及編號；23-礦點；24-研究區(qū)位置Fig.1　Regional geological map of the Xindalai area (modified from Tao et al, 2009a)

1.2　構造

區(qū)域上主體構造線方向為NE向(包括NEE向)并控制著地層和巖體的展布方向，NW向次之，再次為近EW向、近SE與NNE向構造。斷裂構造長度多為數(shù)十千米，小者數(shù)千米—幾百米，主要形成于海西期與燕山期。區(qū)內(nèi)主要控礦構造為NW向斷裂構造。

1.3　巖漿巖

研究區(qū)巖漿活動強烈，自海西期至燕山晚期，均有巖漿侵入，脈巖類從基性到酸性均有發(fā)育。燕山晚期花崗巖類與成礦密切相關，這些巖體呈巖株狀產(chǎn)出，面積零點幾—幾平方千米，出露于大型侵入體內(nèi)部及地層邊部，巖體上部圍巖礦化石英脈發(fā)育。

2礦床模型與找礦概念模型

2.1　礦床模型

通過對烏蘭德勒銅鉬礦(陶繼雄等，2009b)、準蘇吉花銅鉬礦(劉翼飛等，2012)、達來敖包鉬礦(趙志軍等，2014)及烏花敖包鉬礦(孔維瓊等，2012)的礦床產(chǎn)出特征，總結出新達來地區(qū)礦床產(chǎn)出模型(圖2)。區(qū)內(nèi)礦床主要以鉬為主，以燕山晚期花崗巖為中心的巖漿熱液成礦系統(tǒng)，成礦時代除準蘇吉花外主要在燕山晚期，成礦巖體為細粒二長花崗巖，圍巖主要為寶力高廟組一段變質(zhì)粉砂巖。礦化類型與圍巖密切相關，在成礦巖體內(nèi)以浸染狀礦化為主，圍巖中以細脈狀礦化為主。圍巖蝕變主要為為鉀化(黑云母化)、黃鐵絹英巖化、硅化、云英巖化，具有一定的分帶現(xiàn)象。其中黃鐵絹英巖化與成礦密切相關。

圖2　新達來地區(qū)礦床模型圖1-第四系；2-上石炭—下二疊統(tǒng)寶力高廟組一段粉砂質(zhì)板巖；3-侏羅紀(二長)花崗巖；4-二疊紀花崗閃長巖；5-二疊紀石英閃長巖；6-石炭紀花崗巖；7-石英脈；8-花崗巖脈；9-浸染狀礦體；10-細脈狀礦體；11-巖性界線；12-巖相界線；13-鉀化硅化帶；14-絹云母化帶①-烏蘭德勒式鉬銅礦；②-準蘇吉花式鉬銅礦；③-達來敖包式鉬礦；④-烏花敖包式鉬礦Fig.2　Deposit model of the Xindalai area

2.2　找礦概念模型

在礦床模型的基礎上，通過對研究區(qū)地質(zhì)、地球化學、地球物理找礦信息的詳細分析，總結區(qū)域成礦背景和地質(zhì)、地球物理、地球化學找礦標志，將研究區(qū)鉬多金屬礦產(chǎn)成礦預測的找礦概念模型列表歸納如表1。

表1　新達來地區(qū)鉬多金屬礦找礦概念模型表

3找礦信息提取

搜集了新達來地區(qū)1∶5萬地質(zhì)礦產(chǎn)、地球化學、地球物理等數(shù)據(jù)資料，并進行處理、轉(zhuǎn)換，建立相應基礎空間數(shù)據(jù)庫。根據(jù)找礦概念模型進行找礦信息提取。

3.1　構造信息提取

研究區(qū)構造控礦明顯，NE向大型構造為主要導巖導礦構造；次級的NW向、NE向及近EW向斷裂為主要的控礦構造，尤其是2組或多組斷裂交匯部位易成礦。受覆蓋影響，研究區(qū)的實測構造數(shù)量較少，分布不均。通過實地調(diào)查和航片解譯，已查明多條斷裂。

在GeoDAS軟件中對斷裂構造進行緩沖區(qū)分析緩沖，選擇多環(huán)方式，環(huán)間距300 m，共15環(huán)，最佳緩沖距離由軟件中證據(jù)權方法確定。當緩沖距離為2 km時，有50%的礦(床)點落入斷裂緩沖區(qū)；當緩沖距離為4 km時，有82%的礦(床)點落入斷裂緩沖區(qū)內(nèi)。

3.2　巖體信息提取

3.2.1巖體解譯研究區(qū)地形比較平坦，牧草及第四系覆蓋嚴重，少量第三系覆蓋，出露的基巖有限。而巖漿巖對成礦起著至關重要的作用，選取適當?shù)姆椒ㄍ茢嘟庾g新生界覆蓋層下的中酸性巖體是成礦預測的關鍵。研究區(qū)侵入巖體主要為花崗巖、閃長巖等中酸性巖體，其地球化學特征、地球物理特征與圍巖地層有著明顯的差別。與中小比例尺不同的是中小比例尺有重力和高磁數(shù)據(jù)，以及巖體的主要成分的元素化探數(shù)據(jù)，而1∶5萬只有高磁和一些金屬成礦元素。因此，選取能夠代表中酸性巖體主成分的地球化學元素(即在巖體區(qū)呈正異?；蜇摦惓?以及高磁數(shù)據(jù)，運用主成分分析方法，知識驅(qū)動加數(shù)據(jù)驅(qū)動相結合，選取代表中酸性巖體的主成分，并結合地質(zhì)事實，全面綜合地解譯隱伏巖體。

通過各元素異常圖與中酸性巖體和地層吻合的對比可知，元素Co、Ni、Bi、As、Sb、Hg與已知出露的地層吻合較好，而Pb與中酸性巖體(尤其是石炭紀的中酸性巖體)吻合較好，因此選取7個微量元素組合進行主成分分析，解譯中酸性巖體。為了避免單獨地球化學方法的局限性，且高磁異常與巖體吻合也較好，因此又加入高磁異常圖進行主成分分析，選擇總累計方差比率為70%，得到PCA1、PCA2、PCA3、PCA4共4個主成分，其中PCA1方差貢獻為37.9%，PCA2方差貢獻為17.3%。

PCA1方差貢獻最高，是最重要的一個主成分，高磁、Pb和Bi在PCA1中的貢獻率為正，其余5種元素在PCA1中的貢獻率亦均為正，且絕對值較大。在主成分得分圖上，其負異常區(qū)域與出露的中酸性巖體吻合很好，代表了研究區(qū)的中酸性巖體。而其正異常與寶力高廟組一段(C2-P1bl1)即碎屑巖段吻合很好，代表了研究區(qū)寶力高廟組一段地層。PCA2表現(xiàn)為正的Pb、Bi、Ni、Co及高磁，負的Hg、As、Sb。在研究區(qū)西北部，元素Pb與石炭紀中酸性巖體對應很好，故PCA2正異常主要代表石炭紀中酸性巖體，負異常與二疊紀中酸性巖體對應較好。在東南部正好相反，負異常對應石炭紀中酸性巖體。結合PCA1和PCA2解譯出巖體和地層，實現(xiàn)對研究區(qū)覆蓋層的初步“揭蓋”(圖3)。

圖3　新達來地區(qū)覆蓋層揭蓋前后對比圖Fig.3　Comparison of geological maps between original and after the removal of covers in the Xindalai area(a) Original geologica map; (b) Geological map after the removal of covers

3.2.2巖體信息提取研究區(qū)成礦巖體主要為燕山期花崗巖，因此選擇侏羅紀巖體做緩沖分析，選擇多環(huán)方式，環(huán)間距300 m，共15環(huán)。當緩沖距離為2 km時，僅有25%的礦(床)點落入斷裂緩沖區(qū)。但為了突出燕山期巖體的權重貢獻，仍然將其作為證據(jù)權圖層。

3.2.3巖體與地層接觸帶信息提取研究區(qū)成礦巖體主要為燕山期花崗巖，但也有很多礦床以海西期巖體與晚古生代地層為圍巖，因此選擇燕山期及以前的巖體與晚古生代地層在GeoDAS軟件中做交界分析，得到巖體與地層的接觸帶shape文件，并對其進行緩沖分析，選擇多環(huán)方式，環(huán)間距300 m，共15環(huán)。當緩沖距離為1 km時，有45%的礦(床)點落入斷裂緩沖區(qū)；當緩沖距離為2 km時，有70%的礦(床)點落入斷裂緩沖區(qū)內(nèi)；當緩沖距離為3 km時，有88%的礦(床)點落入斷裂緩沖區(qū)內(nèi)。

3.3　巖脈信息提取

諸多礦點與石英脈、石英斑巖脈等巖脈有關，因此在GeoDAS軟件中巖脈進行緩沖分析，選擇多環(huán)方式，環(huán)間距200 m，共15環(huán)。當緩沖距離為1 km時，有52%的礦(床)點落入斷裂緩沖區(qū)；當緩沖距離為2 km時，有70%的礦(床)點落入斷裂緩沖區(qū)內(nèi)；當緩沖距離為3 km時，有81%的礦(床)點落入斷裂緩沖區(qū)內(nèi)。

3.4　地球化學數(shù)據(jù)處理與弱緩異常提取

為討論覆蓋層物質(zhì)組成(各元素含量)與基巖區(qū)及其各單元物質(zhì)組成可能的關系，對各單元元素進行標準化處理(各地質(zhì)單元元素平均含量與全區(qū)元素平均含量的比值，即相對于全區(qū)的富集系數(shù)或襯度系數(shù))，然后作相對富集程度模式圖(圖4)。對比分析覆蓋區(qū)與基巖出露區(qū)可知：Mo、Pb、Au、Zn、Hg等元素更易富集在覆蓋層中，而Bi元素更易富集于巖漿巖出露區(qū)。

圖4　新達來地區(qū)各地質(zhì)單元元素相對富集程度模式圖Fig.4　Models showing relative enrichment of element in each geological unit of the Xindalai area

對研究區(qū)內(nèi)14種元素作主成分分析。選擇相關系數(shù)矩陣計算模型，總累計方差貢獻率>70%，得到PCA1、PCA2、PCA3、PCA4、PCA5共5個主成分。

其中，PCA1主成分的方差貢獻最高，貢獻率達36.6%，是最重要的一個主成分，14種元素在PCA1中的貢獻率均為正，主要有Zn、As、Cu、Sb、W、Ni、Mo、Sn、Co、Ag、Hg等，包括了低溫元素、中溫元素和高溫元素的組合，代表了研究區(qū)最主要的成礦富集事件，該主成分對尋找W、Mo、Cu、Zn、Pb等多金屬共生礦床有直接指示作用。而Au、Pb、Bi在PCA1中的貢獻較小，表明此次成礦事件未造成Au、Pb、Bi等元素的大量富集。

PCA2主成分的方差貢獻較大，貢獻率為15.8%，Bi、Pb、Sn、Mo、W、Ag元素貢獻率均為正，且絕對值較大；Sb、Hg、Co元素元素貢獻率均為負，且絕對值較大，代表研究區(qū)可能存在1期以Pb、Mo、W為主的中高溫成礦事件。

從PCA1和PCA2主成分得分圖(圖5a、b)中可以看出，異常基本分布在基巖出露區(qū)，而在第四系覆蓋區(qū)基本無異?；騼H顯示較弱異常，異常信息受到了明顯的抑制與屏蔽，現(xiàn)已發(fā)現(xiàn)的礦床(點)無一例外地全部分布于基巖當中，而在覆蓋層下找礦潛力巨大，所以合理地識別和加強覆蓋區(qū)的弱緩異常對研究區(qū)找礦尤為重要。S-A法能夠分解出復合疊加的背景和異常(Cheng et al, 2000；Zuo et al, 2013)，通過對PCA1和PCA2進行S-A分解后得到二者的異常分布圖(圖5c、d)，使原來在基巖區(qū)大面積分布的異常被分解成多個點狀異常，尤其在無異常顯示的覆蓋區(qū)也能增強顯示出弱異常，從而可以有效地提取覆蓋區(qū)域的弱、小、隱蔽異常。分解后的異常圖中礦床(點)多分布在異常高值上，尤其是產(chǎn)于巖體中的礦床(點)與異常高值吻合很好。

3.5　預測量構建與優(yōu)化

根據(jù)前面對找礦模型的構建以及對各個找礦信息的提取，采用數(shù)據(jù)驅(qū)動與知識驅(qū)動相結合的方式，確定了以下用于研究區(qū)成礦預測的有利預測變量(表2)：斷裂構造緩沖分析、石英脈、石英斑巖脈緩沖分析、巖體與地層接觸帶緩沖分析、地球化學PCA1主成分S-A分解異常、地球化學PCA2主成分S-A分解異常、侏羅紀巖體緩沖分析。雖然侏羅紀緩沖分析是當緩沖距離為2 km時，僅有25%的礦床(點)落入斷裂緩沖區(qū)，但為了突出燕山期巖體的權重貢獻，仍然將其作為證據(jù)權圖層。

表2　新達來地區(qū)預測變量表

4基于GeoDAS綜合信息成礦預測

選擇研究區(qū)同為訓練區(qū)及預測區(qū)，26個已知的金屬礦床(點)作為訓練點。在GeoDAS中設定預測面積單元為1 km2時，面積單位數(shù)是25.67(基本上為1個礦床(點)分布于1個單元格)，總單元數(shù)是2 896，先驗概率為0.008 9。在設置訓練參數(shù)后，對各個證據(jù)圖層進行模糊證據(jù)權重計算，然后進行后驗概率的計算，得到后驗概率圖。成秋明等(2009)將奇異性理論運用在個舊錫銅礦預測時指出后驗概率與面積服從分形分布。因此，選擇后驗概率與分布面積的雙對數(shù)曲線拐點，將后驗概率大致分為A、B、C、D共4級，并依據(jù)系統(tǒng)優(yōu)化原則、綜合評判原則和地質(zhì)基礎原則(夏慶霖等，2009)，在研究區(qū)圈定A級遠景區(qū)5個，B級遠景區(qū)7個，C級遠景區(qū)10個(圖6)。

圖5　新達來地區(qū)化探元素主成分得分圖Fig.5　Plots showing scoring of geochemical exploration element main component in the Xindalai area

圖6　新達來地區(qū)后驗概率及找礦遠景區(qū)分布圖Fig.6　Map showing posterior probability and prospective areas in the Xinlaida area

A-1遠景區(qū)。該遠景區(qū)位于研究區(qū)中北部烏蘭德勒一帶，A級后驗概率異常面積大且分級明顯，形態(tài)整體呈圓形。區(qū)內(nèi)現(xiàn)已發(fā)現(xiàn)礦床(點)5處，其中包括烏蘭德勒大型鉬銅礦床。區(qū)內(nèi)大面積出露石炭紀花崗巖，二疊紀花崗巖呈橢圓狀分布于石炭紀巖體中，石英脈、石英斑巖脈、花崗巖脈等巖脈較發(fā)育。

A-2遠景區(qū)。該遠景區(qū)位于研究區(qū)中西部準沖格熱格一帶，包含2處幾乎連續(xù)的A級后驗概率異常圍繞二疊紀巖體向南凸出，異常面積大且分級明顯。準蘇吉花中型鉬銅礦床位于其中，區(qū)內(nèi)主要侵入巖為二疊紀花崗閃長巖，侵入寶力高廟組一段地層中；斷裂構造較為發(fā)育，主體構造與區(qū)域構造線一致，呈NE向展布，控巖控礦構造；一系列NW向次級斷裂構造為控礦構造，被后期石英脈所充填，并形成總體走向NW方向雁行式排列的含礦石英脈群。

A-3遠景區(qū)。該遠景區(qū)位于研究區(qū)西北部達來敖包一帶，包含2處A級后驗概率異常，其中一處呈橢圓狀NNE向展布，分布面積大，且達來敖包鉬礦床位于異常中心；另一處鉛鋅礦點位于該異常邊部。區(qū)內(nèi)侵入巖為石炭紀黑云母花崗巖、二疊紀堿長花崗巖以及少量侏羅紀石英斑巖，均侵入于寶力高廟組一段地層中。區(qū)內(nèi)NE向斷裂構造發(fā)育，為主要控礦構造，大多數(shù)礦化與之有密切關系。NW向石英脈、花崗巖脈等巖脈較發(fā)育。

A-4遠景區(qū)。該遠景區(qū)位于研究區(qū)西部烏花敖包一帶，包含1處較大的和幾處較小的A級后驗概率異常，較大的異常面積大且分級明顯。烏花敖包鉬銅礦床位于較大的A級后驗概率異常中心。燕山期石英斑巖體呈NE向脈群狀產(chǎn)出，區(qū)內(nèi)礦化體與細?；◢忛W長巖、石英斑巖有著密切關系。各類巖石中裂隙構造十分發(fā)育，裂隙構造的發(fā)育程度與礦化強度的強弱有著正相關關系。

A-5遠景區(qū)。該遠景區(qū)位于研究區(qū)北部沙拉楚魯圖一帶，A級后驗概率異常面積大且分級明顯，呈NE向展布，區(qū)內(nèi)已發(fā)現(xiàn)1個鉛鋅礦點。區(qū)內(nèi)大面積出露石炭紀花崗巖，東北部出露少量地層為上石炭—下二疊統(tǒng)寶力高廟組一段陸相正常碎屑-砂泥質(zhì)沉積建造，巖性主要為灰綠色變質(zhì)中細粒長石石英砂巖、粉砂質(zhì)板巖。遠景區(qū)構造發(fā)育，總體方向近NE向，擁有NW向、近SN向、近EW向多組斷裂，NE向、NS向為控巖構造，NW向為控礦構造；石英脈、石英斑巖脈、花崗巖脈等巖脈較發(fā)育。

5結論

(1) 通過非線性理論和方法可以實現(xiàn)牧草與新生界覆蓋層的“揭蓋”以及地球化學弱異常的提取。

(2) 以GeoDAS 礦產(chǎn)預測軟件為平臺，建立預測模型，運用模糊證據(jù)權方法進行信息綜合集成，在新達來草原覆蓋區(qū)共圈定A級找礦遠景區(qū)5個、B級遠景區(qū)7個、C級遠景區(qū)10個。

(3) 在A-2遠景區(qū)選取1處成礦條件較好的選區(qū)D10進行大比例尺驗證工作，布置1∶1 萬地質(zhì)草測、土壤地球化學測量、地面高精度磁測。在此基礎上，選擇1處較好的異常進行1∶2 000激電中梯及瞬變電磁剖面測量，發(fā)現(xiàn)該區(qū)中部存低阻高激化體，在L420與L430線地下300 m左右存在多處規(guī)模小的低阻體，推測為隱伏礦化引起。巖石地球化學測量也表明該區(qū)花崗閃長巖屬于地殼伸展作用產(chǎn)生的造山期后高分異的過鋁質(zhì)高鉀鈣堿系列A型花崗巖類侵入巖，巖體中的稀土與烏蘭德勒輝鉬礦稀土元素分配曲線形態(tài)十分相似，因此推斷具有一定的成礦潛力。

(4) 草原覆蓋區(qū)成礦預測的技術流程可概括為：明確找礦目標，從區(qū)域地質(zhì)背景分析和典型礦床研究入手，建立礦床模型和找礦概念模型，通過對覆蓋層“揭蓋”、推斷隱伏成礦要素和提取找礦信息，構建和優(yōu)化預測變量，劃分預測單元，進而計算成礦有利度或后驗概率，圈定找礦遠景區(qū)，最后進行大比例尺查證工作。

參考文獻：

成秋明，趙鵬大，張生元，等.2009.奇異性理論在個舊錫銅礦產(chǎn)資源預測中的應用：綜合信息集成與靶區(qū)圈定[J].地球科學:中國地質(zhì)大學學報，34(2)：243-252.

孔維瓊，劉翠，鄧晉福，等.2012.內(nèi)蒙古二連浩特地區(qū)烏花敖包鉬礦區(qū)火成巖特征和LA-ICP-MS鋯石年代學研究[J].地學前緣，19(5)：123-135.

劉翼飛，聶鳳軍，江思宏，等.2012.內(nèi)蒙古蘇尼特左旗準蘇吉花鉬礦床成巖成礦年代學及其地質(zhì)意義[J].礦床地質(zhì)，31(1)：119-128.

毛景文，謝桂青，李曉峰，等.2005.大陸動力學演化與成礦研究：歷史與現(xiàn)狀：兼論華南地區(qū)在地質(zhì)歷史演化期間大陸增生與成礦作用[J].礦床地質(zhì)，24(3)：193-205.

任紀舜，陳廷愚，牛寶貴.1992.中國東部及鄰區(qū)大陸巖石圈的構造演化與成礦[M].北京：科學出版社.

邵積東，王守光，趙文濤，等.2008.內(nèi)蒙古中部地區(qū)重要成礦區(qū)帶成礦地質(zhì)特征及找礦潛力分析[J].西部資源，(1)：59-61.

陶繼雄，鐘仁，王弢，等.2009a.內(nèi)蒙古達來廟地區(qū)戰(zhàn)略性礦產(chǎn)遠景調(diào)查報告[R].內(nèi)蒙古呼和浩特：內(nèi)蒙古自治區(qū)地質(zhì)調(diào)查院.

陶繼雄，王弢，陳鄭輝，等.2009b.內(nèi)蒙古蘇尼特左旗烏蘭德勒鉬銅多金屬礦床輝鉬礦錸-鋨同位素定年及其地質(zhì)特征[J].巖礦測試，28(3)：249-253.

夏慶霖，成秋明，左仁廣，等.2009. 基于GIS礦產(chǎn)勘查靶區(qū)優(yōu)選技術[J].地球科學:中國地質(zhì)大學學報，34(2)：287-293.

趙志軍，劉立華，魏雅玲，等.2014.內(nèi)蒙古達來敖包鉬多金屬礦地質(zhì)特征及找礦潛力淺析[J].西部資源，(2)：102-105.

CHENG QIUMING， XU YAGUANG, GRUNSKY E. 2000. Integrated spatial and spectrum method for geochemical anomaly separation[J]. Natural Resources Research, 9(1): 43-51.

XIAO WENJIAO, WINDLEY B F, HAO JIE, et al. 2003. Accretion leading to collision and the Permian Solonker suture, Inner Mongolia, China: Termination of the central Asian orogenic belt[J]. Tectonic, 22(6): 1069-1088.

ZUO RENGUANG, XIA QINLIN, ZHANG DAOJUN. 2013. A comparison study of the C-A and S-A models with singularity analysis to identify geochemical anomalies in covered areas[J]. Applied Geochemistry, 33:165-172.

Metallogenic prediction of Mo-Cu polymetalic deposits in the Xindalai grassland-covered area, Inner Mongolia

XIONG Shuang-cai1,2, XIA Qing-lin1,3, YAO Chun-liang1, ZHANG Xiao-jun1, QIN Ying1, ZHANG Kun4

(1.The Faculty of Earth Resources, China University of Geosciences (Wuhan), Wuhan 430074, Hubei, China; 2. No.1 Geological Party, Xinjiang Bureau of Geology and Mineral Resources, Changji 831100, Xinjiang, China; 3.Collaborative Innovation Center for Exploration of Strategic Mineral Resources, Wuhan 430074, Hubei, China; 4. 6th Geology & Mineral Resources Survey Institute, Shandong, Weihai 264209, Shandong, China)

Abstract：The Xindalai grassland-covered area is located in the late paleozoic continental-margin accretionary belt at the southwestern margin of the Siberia Plate. It belongs to the late Paleozoic-Mesozoic Erlian-Dongwuqi Cu-Mo-Fe metallogenic belt, part of the Palaeo-Asian Metallogenetic Domain, with a great prospecting potential. This area is largely covered by pasture and the Quaternary sediments, and the prospecting work has failed to achieve breakthrough because of the poorly outcropped bedrock. Therefore, this study discussed the regional metallogenic background and geological features of typical deposits to analyze ore-controlling factors and prospecting indicators, and established a comprehensive prospecting model to extract useful ore-related geological, geochemical, geophysical and remote sensing information. In addition, we used the GIS technology to uncover the cover layer and extract favorable metallogenic information. In the GeoDAS mineral prediction software platform, we utilized the multi-fractal filtering technology to distinguish between background and anomalies and identify weak mineralization information. Comprehensive ore-related information was integrated by the fuzzy weights of evidence method formation to form a posterior probability prediction map, and prospecting targets were delineated and verified. This study explores a new way of prospecting and assessing the endogenic metallic deposits in grassland-covered areas of Inner Mongolia.

Keywords：grassland-covered areas; GeoDAS software; nonlinear models; fuzzy weights of evidence method; metallogenic prediction; Xindalai in Inner Mongolia