基于MRAS的鞍山—本溪地區鐵礦資源潛力分析

姚志宏 孫鵬慧 劉長純 李志鋒 王成龍

(遼寧省地質礦產調查院，遼寧 沈陽 110032)

·地質與測量·

基于MRAS的鞍山—本溪地區鐵礦資源潛力分析

姚志宏 孫鵬慧 劉長純 李志鋒 王成龍

(遼寧省地質礦產調查院，遼寧 沈陽 110032)

簡要介紹了遼寧省鞍山—本溪地區區域地質特征及航磁、重力場特征，在此基礎上建立了鞍本地區鞍山式沉積變質型鐵礦區域預測模型：鐵礦產于鞍山—本溪太古宙綠巖盆地；成礦時代是新太古代；含礦建造是新太古代鞍山群櫻桃園巖組、茨溝巖組；礦物組合以磁鐵礦和赤鐵礦為主；磁鐵礦具強磁性，磁異常為其重要的找礦標志。基于礦產資源評價系統(MRAS)，采用綜合信息法對鞍本地區鞍山式鐵礦進行了預測單元的圈定。根據預測模型，構置了10個預測要素變量，經篩選保留了磁鐵石英巖存在標志、磁鐵石英巖面積、磁異常存在標志、磁異常面積、含礦建造存在標志及鐵礦床點存在標志等6個變量。對各預測單元成礦有利度進行分類，從而實現了預測單元的優選，最終圈定A類預測區11個，B類預測區9個，C類預測區24個。采用地質經濟模型法(非總和式)對區內鞍山式鐵礦資源量進行了定量分析，共確定預測資源量476億t。

鞍山式鐵礦 MRAS 預測單元圈定 資源量定量分析

遼寧省鐵礦資源豐富，以新太古代鞍山式鐵礦為主，其儲量占全國總量的23%，鞍山—本溪地區鞍山式鐵礦占全省總儲量的96%以上。該地區鐵礦資源經過50多a的勘查和開發，目前累計探明儲量154億t。

Mineral Resource Assessment System(簡稱MRAS)是由中國地質科學院礦產資源研究所開發的基于GIS的區域礦產資源評價系統。該系統經過多個大學、研究所及地勘單位的使用，取得了較好的效果，可以作為本次鞍本地區鞍山式鐵礦資源評價的工具。

礦產資源定量預測是在一個較大區域內對某種礦產資源的總量進行估算。人們不僅希望得到總的資源量，而且還希望明確知道某種礦產資源潛力的具體分布形式及找礦靶區。本次資源潛力分析主要是通過對研究區預測要素的分析而建立找礦模型、預測單元，對多元地質變量進行提取、賦值和篩選，通過特征分析法對預測區進行優選，進而實現礦產資源定位預測[1]；采用地質經濟模型法對各預測區進行資源量定量分析，最終實現整個研究區資源量定量預測。

1 研究區概況

鞍山—本溪地區是遼寧境內重要的鐵礦成礦帶，目前已發現鐵礦床(點)近百處，其中已探明的超大型礦床6處，大型礦床12處，中型礦床達十幾處，小型礦床幾十處。

研究區大地構造單元屬于華北地臺、膠遼臺隆、太子河渾江臺陷、遼陽—本溪凹陷及營口—寬甸臺拱、鳳城凸起的西北部[2-4]。該地區地層屬遼東地層分區，出露地層主要有太古界鞍山群，下元古界遼河群，上元古界青白口系和古生界及新生界第四系。鞍山式鐵礦賦存于太古界鞍山群的茨溝組和櫻桃園組中，主要為磁鐵石英巖和角閃磁鐵石英巖型鐵礦(即鞍山式鐵礦)[5-6]。茨溝組地層發育在遼陽、本溪等地。疊置厚度大于2 200 m。原巖以基性熔巖為主，夾含鐵硅質巖，少量陸源碎屑沉積，上部中酸性凝灰巖增多。原巖建造為基性火山巖鐵硅質巖建造，主要分布于歪頭山—梨樹溝和弓長嶺—本溪等地。鞍山群上部櫻桃園組含鐵礦層分布于櫻桃園，眼前山，東、西鞍山等地，貧鐵礦賦存在櫻桃園巖組上部。巖石類型主要有白云石英片巖、綠泥片巖、絹云石英綠泥片巖、千枚巖、二云片巖、黑云變粒巖夾斜長角閃巖及磁鐵石英巖，疊置厚近1 000 m。原巖主要為黏土巖、泥質砂巖，部分英安質疑灰巖、酸性凝灰巖夾厚層含鐵硅質巖。原巖建造為碎屑巖—黏土巖—含鐵硅質巖建造。

鞍本地區為平穩高磁異常區，區內航磁異常具有規模大、強度高、成群成帶、集中分布的特點。異常多為太古界鞍山群磁鐵石英巖引起。在太古宙基底出露的區域中，已發現鐵礦床引起的航磁異常僅是鞍本地區整個航磁異常的一部分，而其他蓋層分布區同樣具有強度高、面積大的航磁異常[7]，由此可見，該區域深部隱伏的鐵礦資源潛力巨大。

研究區重力場的變化總體呈北東向展布，海城—歪頭山地區重力場呈現北東向條帶狀中心遞增式重力高，異常較穩定，幅值為23 mgal。營口—橋頭地區重力場以重力低為主，局部有重力高出現，異常特征形態不明顯，中心幅值為42 mgal。推斷重力低場區為侵入巖體分布，重力高場區由花崗片麻巖及太古宙變質巖引起。從鐵礦分布情況看，大部分鞍山式鐵礦的背景場為重力高場，并多分布在重力梯度帶上。

2 預測區圈定

2.1 建立預測模型

通過對研究區內西鞍山和歪頭山2個代表性礦床的研究，總結鞍本地區鞍山式鐵礦預測模型。①成礦構造背景為古華北板塊的東北部鞍山—本溪太古宙綠巖盆地；②成礦時代為新太古代；③含礦建造為新太古代鞍山群櫻桃園巖組、茨溝巖組的巖石組合，如綠泥片巖、千枚巖、斜長角閃巖、磁鐵石英巖等；④礦物組合以磁鐵礦和赤鐵礦為主；⑤磁鐵礦的強磁性成為這類鐵礦床獨特的物理性質，磁異常為其重要的找礦標志；⑥磁鐵礦石具有較高的密度，鐵礦體與圍巖間具有較大的密度差，因此，重力異常在預測該類型鐵礦床時，可作為一項重要預測要素來考慮。

2.2 圈定預測單元

預測單元的圈定通過建模器實現。首先根據預測模型構置相應的圖層，具體包括：①構置磁異常存在范圍。將磁異常所控制的范圍以磁異常面文件進行表達。②構置磁鐵石英巖的分布范圍。從遼寧省鞍山—本溪太古宙變質建造構造圖中提取磁鐵石英巖的面文件。③構置礦床點緩沖區控制范圍。求取礦床點作緩沖區分析，構置礦床點緩沖區的面文件。④將磁鐵石英巖圖層和礦點緩沖區圖層在MapGIS中疊加(求并集)，構成一個新的礦產圖層。

在MARS2.0中用建模器建模，使磁異常的面文件與新構置礦體文件疊加求并。每個單元網格的面積為1 km2，以滿足1∶5萬精度的預測工作需要。

大千世界，人海茫茫，少年血氣方剛但缺乏經驗，老年人經驗豐富但氣血兩虧，惟有中年人，既有經驗閱歷，穩重沉著，又如日中天，年富力強，是行行業業的中堅力量，左突右沖，擔當起方方面面的重任。而且承上啟下，前后呼應。中年人的狀態如何，決定著國家的命運如何；中年人的前途怎樣，關系著民族的前途怎樣。因而，我們看到成千上萬的中年人在兢兢業業工作，孜孜矻矻奮斗，在創造一個個奇跡，在刷新一個個記錄。

2.3 預測要素變量的構置與選擇

根據前面總結的鞍本地區鞍山式鐵礦各預測要素，選取了與鞍山式鐵礦有關的6個專題圖層作為預測區優選的變量，分別為磁鐵石英巖(面文件)、致礦磁異常(面文件)、太古宙含礦變質建造(面文件)、重力剩余異常平面等值線(線文件)、航磁化極平面等值線(線文件)、鐵礦床點(點文件)。

在MRAS中，對這6個圖層分別提取10個變量作為鞍本地區鞍山式鐵礦的預測要素變量：磁鐵石英巖存在標志、磁鐵石英巖面積、磁異常存在標志、磁異常面積、含礦建造存在標志、重力剩余異常平面等值線的最高值、重力剩余異常平面等值線的均值、航磁化極平面等值線的最高值、航磁化極平面等值線的平均值、鐵礦床點的存在標志。

應用相似系數法對預測變量進行統計分析，重力剩余異常的最高值和均值以及航磁異常的最大值和均值這4個變量的相關性較差，可以剔除。篩選后的變量分別是磁鐵石英巖存在標志、磁鐵石英巖面積、磁異常存在標志、磁異常面積、含礦建造存在標志及鐵礦床點存在標志6個變量。

2.4 預測區優選

預測區優選主要通過對各預測區成礦有利程度的分析來實現，剔除成礦不利的區域，并對優選后的預測區分別用3種顏色來區分表示成礦優劣的A、B、C 3個級別。

鞍山—本溪地區經過優選后的預測區為44個，其中成礦條件十分有利、成礦匹配程度高、資源潛力大的預測區(A類)11個；成礦條件較為有利、成礦匹配程度較高、預測潛力較大的預測區(B類)9個；具有成礦條件、有可能發現資源、可作為探索的其他預測區(C類)24個，見圖1。

圖1 優選后的預測區Fig.1 Prediction regions after optimally selecting

3 資源量定量分析

(1)構造鞍本地區品位-噸位模型。在MRAS系統中雖然已經提供了全國沉積變質型鐵礦的品位-噸位模型，但是考慮到鞍本地區的鐵礦床成礦較好，大型和超大型礦床較多，其總體成礦性能要優于全國的平均水平，因此，使用鞍本地區內已知礦床的數據來構造新的品位-噸位模型。本次共收集到鞍本地區內58個鞍山式鐵礦床的品位及相應的資源量的資料。結合地質勘探報告，對礦區的外圍和2 000 m以淺深度的資源量進行了推測，并將該推測資源量與探明儲量的和作為模型噸位值，見表1。基于表中數據，在MRAS中構造鞍本地區的品位-噸位模型。

表1 鞍山—本溪地區鞍山式鐵礦品位-噸位數據 Tab.1 Grade-Tonnage data table of the Anshan-type iron deposits in the Anshan-Benxi area

(2)礦床數估計。在MRAS中采用非總和式資源量估算方法。由地質專家根據礦產地信息、磁異常等多種地質證據對各預測區中可能存在的礦床數按90%、50%、10%的概率分別進行估計[10]。

(3)蒙特卡羅資源量模擬。通過蒙特卡羅模擬，確定了第1個預測區的預測結果，包括礦點數、品位及礦石量。對于鐵礦而言，其資源量為礦點數與礦石量的乘積。保存第1個預測區的預測結果并對其他預測區進行估算，分別獲得各個子區的預測資源量，見表2。經估算，各個子區的預測資源量合計為226.5億t，模型區的外圍和深部的資源量預測值為249.5億t，因此，鞍本地區預測資源量為476億t。

表2 預測資源量 Tab.2 Prognosis of resources table

4 結 論

(1)通過對鞍山—本溪地區成礦地質條件及重、磁場特征的分析，建立了鞍山式鐵礦預測模型。

(2)基于礦產資源GIS評價系統(MRAS)在遼寧省鞍山—本溪地區共圈定預測區44個，其中A類最小預測區11個，B類最小預測區9個，C類最小預測區24個。

(3)采用地質經濟模型法(非總和式)預測鞍山式沉積變質型鐵礦資源量為476億t。

[1] 肖克炎,婁德波,陰江寧,等.中國鐵礦資源潛力定量分析[J].地質通報,2011,30(5):650-651. Xiao Keyan,Lou Debo,Yin Jiangning,et al.Quantified research on iron potential of China[J].Geological Bulletin of China,2011,30(5):650-651.

[2] 楊占興.遼寧省成礦系列與成礦區帶研究[D].北京:中國地質大學,2006：16-18. Yang Zhanxing.Metallogenic Series and Metallogenic Belts in Liaoning Province[D].Beijing:China University of Geosciences,2006:16-18.

[3] 遼寧省地質礦產局.遼寧省區域地質志[M].北京:地質出版社.1989:707-727. Bureau of Geology and Mineral Resources of Liaoning Province.Regional Geology of Liaoning Province[M].Beijing:Geological Publishing House.1989:707-727.

[4] 孫鵬慧.遼寧省成礦區帶劃分及特征[J].中國科技博覽,2013(31):277-281. Sun Penghui.The division and character of mineralization belt,Liaoning province[J].China Science and Technology Review,2013(31):277-281.

[5] 鄭峻慶,張寶化,蔡一廷,等.遼寧省本溪地區北臺—歪頭山一帶太古宙鞍山群主要構造特征及其對鐵礦的控制[J].地質找礦論叢,1986,1(1):20-24. Zheng Junqing,Zhang Baohua,Cai Yiting,et al.Tectonic characteristics of the archean Anshan group and their effects on iron ore deposits in area of Beitai to Waitoushan,Benxi,Liaoning province[J].Contributions to Geology and Mineral Resources Research,1986,1(1):20-24.

[6] 張 朋,彭明生,歐陽兆灼,等.遼寧鞍本地區鐵礦床地質特征及找礦標志分析[J].地質與資源,2012,12(6):516-521. Zhang Peng,Peng Mingsheng,Ouyang Zhaozhuo,et al.Geological characteristics and ore-searching guides of the iron deposits in Anshan-Benxi area,Liaoning province[J].Geology and Resources,2012,6(12):516-521.

[7] 王福利.淺談鞍山式富鐵礦地質特征及找礦方向[J].吉林地質,2008,27(2):45-47. Wang Fuli.Geologic features and prospecting direction of the Anshan type iron-rich mineral[J].Jilin Geology,2008,27(2):45-47.

[8] 肖克炎,葉天竺,李景朝,等.礦床模型綜合地質信息預測資源量的估算方法[J].地質通報,2010,29(10):1404-1412. Xiao Keyan,Ye Tianzhu,Li Jingchao,et al.Integrated geo-information modeling reserve estimation method[J].Geological Bulletin of China,2010,29(10):1404-1412.

[9] 肖克炎,張曉華,李景朝,等.全國重要礦產總量預測方法[J].地學前緣,2007,14(5):21-24. Xiao Keyan,Zhang Xiaohua,Li Jingchao,et al.Prediction methods of the national important mineral total[J].Earth Science Frontiers,2007,14(5):21-24.

[10] 張 璟,李守義,徐 山,等.基于MRAS證據權重法的遼西錦州—阜新金礦化帶金礦潛力評價[J].中南大學學報:自然科學版,2012,43(9):3565-3566. Zhang Jing,Li Shouyi,Xu Shan,et al.Potential gold evaluation with weighting of evidence based on MRAS in Jinzhou-Fuxin gold metallogenic belt,western Liaoning province,China[J].Journal of Central South University：Natural Science Edition,2012,43(9):3565-3566.

(責任編輯 鄧永前)

Analysis on the Iron Resource Potential in Anshan-Benxi Area Based on MRAS

Yao Zhihong Sun Penghui Liu Changchun Li Zhifeng Wang Chenglong

(LiaoningSurveyAcademyofGeologyandMineralResources,Shenyang110032,China)

The regional geological features，aeromagnetic and gravity field features in Anshan-Benxi area are introduced briefly.The regional prediction model of the Anshan-type iron deposits is established according to the prediction results：Iron ores deposit is located in Archean greenstone Basin of Anshan-Benxi area，and the metallogenic epoch is in Neo-archean.The ores occur in metamorphic rock of Yingtaoyuan formation and Cigou formation.Minerals are mainly magnetite and hematite.With the characteristics of ferromagnetism of magnetite，magnetic anomalies can be as an important prospecting sign.The predicting units for Anshan-type iron ores in Anshan-Benxi Area are delineated by adopting the comprehensive information method based on the mineral resource assessment system (MRAS).According to the predicting model，10 predictive variables are structured，among which six ones are retained by screening，including magnetite quartzite existence symbol and its area，magnetic anomaly existence symbol and its area，ore formation existence symbol and iron ore deposits existence symbol.The metallogenic favorability degree of each prediction unit is classified，so as to realize the optimal selection of prediction units.So，11 A-level predicting areas，9 B-level predicting areas and 24 C-level predicating areas are delineated eventually.The iron ore resources of Anshan type is quantitatively analyzed by using the geological economic method (non-total evaluating)，and the about 4.76×1010t of resources is predicted.

Anshan-type iron deposits，MRAS，Predicting units delineation，Quantitative analysis of resources

2013-11-11

全國礦產資源潛力評價項目(編號：1212011121004)。

姚志宏(1983—)，女，工程師。

P624.6

A

1001-1250(2014)-01-084-04