新疆阿希金礦控礦因素三維定量分析與提取

林玉婷， 毛先成， 劉曉霞， 任　佳， 劉占坤， 潘　敏， 羅　丹， 鄧　浩,朱佳瑋， 張苗苗， 夏　芳， 肖　飛， 魏清峰， 韓建民， 樊紅喜， 饒文天， 黃悠悠

(1.有色金屬成礦預測與地質(zhì)環(huán)境監(jiān)測教育部重點實驗室，湖南長沙410083； 2.中南大學地球科學與信息物理學院，湖南長沙410083； 3.新疆有色金屬工業(yè)集團，新疆烏魯木齊830000)

?

新疆阿希金礦控礦因素三維定量分析與提取

林玉婷1,2， 毛先成1,2， 劉曉霞1,2， 任佳1,2， 劉占坤1,2， 潘敏1,2， 羅丹1,2， 鄧浩1,2,朱佳瑋1,2， 張苗苗1,2， 夏芳3， 肖飛3， 魏清峰3， 韓建民3， 樊紅喜3， 饒文天3， 黃悠悠3

(1.有色金屬成礦預測與地質(zhì)環(huán)境監(jiān)測教育部重點實驗室，湖南長沙410083； 2.中南大學地球科學與信息物理學院，湖南長沙410083； 3.新疆有色金屬工業(yè)集團，新疆烏魯木齊830000)

新疆阿希金礦自20世紀90年代投產(chǎn)至今，地表礦、淺部礦已開采殆盡。采用三維空間分析方法提取控礦因素并進行成礦預測是發(fā)現(xiàn)深部隱伏礦體的一條有效途徑。通過阿希金地質(zhì)體三維建模，采用距離場分析、地質(zhì)界面坡度分析以及地質(zhì)界面趨勢-起伏分析方法，提取了控礦主斷裂面距離場、主斷裂面坡度、主斷裂面趨勢-起伏、管道相距離場、管道相趨勢-起伏5個控礦因素，并通過散點圖對礦化指標與控礦因素的關聯(lián)關系進行分析，得到各個因素的礦化富集空間，提取出成礦有利部位。分析結(jié)果可為阿希金礦的深部隱伏礦體三維預測提供成礦信息支持。

控礦因素；距離場分析；坡度分析；趨勢-起伏分析；三維定量分析；深部隱伏礦體;阿希金礦；新疆西天山

0　引　言

新疆阿希金礦床是20世紀80年代末在新疆西天山地區(qū)發(fā)現(xiàn)的一個大型低硫型淺成低溫熱液金礦床。自1993年開始投入生產(chǎn)，經(jīng)過多年開采，地表礦、淺部礦已開采殆盡，找礦主體已經(jīng)向隱伏礦、深部礦逐漸轉(zhuǎn)變，找礦難度不斷加大。為減少深部探礦的盲目性和風險性，采用合適的方法分析控礦因素、進行成礦定量預測顯得尤為緊迫和重要。翟偉等(2007)分析了阿希金礦流體地化特征及成礦機制；魏佳林(2012)根據(jù)阿希金礦的黃鐵礦標型特征及礦體的剝蝕程度等進行了礦體空間定位規(guī)律及找礦方向研究；王慶峰(2013)根據(jù)阿希金礦古火山機構(gòu)和斷裂構(gòu)造的特征及控礦作用進行了找礦前景評價及成礦預測。這些對阿希金礦進行的成礦預測研究是基于地質(zhì)構(gòu)造或地球化學特征的，尚缺乏基于三維地質(zhì)模型的分析。三維成礦預測方法(趙鵬大等，1992；毛先成等，1998a，1998b，2010；毛先成,2006;陳建平等，2007；肖克炎等，2012；Fallara et al.，2006)近年來已發(fā)展成為融合三維地質(zhì)建模、多元地學信息集成的三維定位定量預測方法。通過三維建模能夠形象直觀地反映各個地質(zhì)體的形態(tài)特征以及空間位置關系，通過對塊體模型進行賦值，利用地質(zhì)統(tǒng)計技術對空間數(shù)據(jù)進行計算，可以實現(xiàn)更科學地定位、定量預測的目的(丁建華等，2009)。

基于阿希金礦的三維地質(zhì)模型進行了控礦因素的定量分析與提取，具體思路為建立地質(zhì)體、地質(zhì)界面的TIN模型，采用三維空間分析方法，對礦床的控礦地質(zhì)因素進行定量分析，提取控礦因素，并分析控礦因素與礦化指標的關系，提取成礦有利部位，為深部隱伏礦體三維預測提供成礦信息支持。

1　礦區(qū)地質(zhì)概況

阿希金礦位于中天山—伊犁地塊北部博羅科努古生代島弧南緣的吐拉蘇盆地。吐拉蘇盆地受控于科古琴南坡斷裂和伊犁盆地北緣斷裂。盆地基底具有雙層結(jié)構(gòu)特征，主要由不同變質(zhì)變形程度的構(gòu)造層組成。下部構(gòu)造層為一套淺海相穩(wěn)定大陸邊緣型的碳酸鹽巖和碎屑巖建造，上部構(gòu)造層由未變質(zhì)的中性、基性火山巖和海相沉積巖建造組成(周圣華等，2008；張耀選等，2015)。在元古界和下古生界基礎上，吐拉蘇盆地在石炭紀發(fā)生了大規(guī)模火山活動(王強等,2006；朱永峰等,2006a,2006b；龍靈利等,2008；Zhu et al.,2005；Xia et al., 2008)，是西北天山重要的斑巖-淺成低溫熱液型金礦礦集區(qū)。盆地中產(chǎn)有阿希、東塔吾爾別克等低硫型和京希—伊爾曼得等高硫型淺成低溫熱液型金礦床及阿庇因迪、西塔吾爾別克等斑巖型金礦床(肖龍,2002；董連慧等,2005；周圣華等,2008；翟偉等,2010；魏佳林,2012；徐伯駿等,2014；張耀選等,2015)。

阿希金礦賦礦圍巖主要為下石炭統(tǒng)大哈拉軍山組管道相英安質(zhì)角礫熔巖，局部地段為英安巖、安山巖和霏細斑巖。礦床總體受火山機構(gòu)及其相關斷裂控制，礦區(qū)內(nèi)發(fā)育系列環(huán)形和放射狀斷裂，礦體嚴格受管道邊緣環(huán)形火山斷裂F2的控制。由于F2斷裂與火山管道相產(chǎn)出狀態(tài)密切相關，控制了礦體的產(chǎn)狀，該斷層可能是火山氣液流體運移的通道。礦區(qū)主礦帶總長度>1 280 m，總體呈向西凸出的弧形帶狀分布，走向近南北，傾角60°～80°。礦石類型被分為2類：石英脈型礦石和蝕變巖型礦石。前者主要以大脈形式產(chǎn)出，其次為細脈、網(wǎng)脈帶形式出現(xiàn)；后者產(chǎn)出形式較為單一，主要是熱液蝕變作用疊加作用于英安質(zhì)角礫熔巖或英安巖。礦區(qū)圍巖蝕變主要為火山期后成礦熱液有關的近礦熱液蝕變作用，表現(xiàn)為硅化、黃鐵絹英巖化、綠泥石化、碳酸鹽化、黏土化等。

圖1　新疆阿希金礦區(qū)地質(zhì)簡圖(據(jù)新疆維吾爾自治區(qū)地質(zhì)礦產(chǎn)勘查開發(fā)局第一地質(zhì)大隊，2011)1-阿恰勒河組第三層凝灰質(zhì)中細粒砂巖、凝灰質(zhì)粗粒砂巖；2-阿恰勒河組第二層凝灰質(zhì)粗粒砂巖；3-阿恰勒河組第一層底礫巖；4-大哈拉軍山組第四層火山集塊巖、安山質(zhì)英安巖、英安巖；5-大哈拉軍山組第三層火山質(zhì)細礫、杏仁狀角閃安山巖、安山質(zhì)英安巖；6-大哈拉軍山組第二層含火山彈火山角礫巖、角閃安山巖；7-大哈拉軍山組第一層安山質(zhì)凝灰?guī)r、凝灰質(zhì)火山角礫巖、安山巖；8-斜長斑巖脈；9-金礦體；10-鉛鋅礦體；11-地層界線；12-斷層及編號Fig.1　Simplified geological map of the Axi gold deposit in Xinjiang(after the First Geological Brigade, Xinjiang Bureau of Geology and Mineral Resources Exploration and Development, 2011)

2　控礦因素三維定量分析方法

控礦因素三維定量分析是以地質(zhì)體或地質(zhì)界面三維模型為基礎進行的控礦地質(zhì)條件定量分析。通過控礦因素剖析，把握礦床成礦機制、時空上的產(chǎn)出及分布特征，在此基礎上總結(jié)礦床成礦規(guī)律，進而利用成礦規(guī)律指導預測找礦工作。控礦因素三維定量分析方法主要分為空間距離分析方法和三維形態(tài)分析方法。空間距離分析方法主要研究礦床(點)的分布與地質(zhì)對象存在的空間距離關聯(lián)關系，三維形態(tài)分析方法主要研究礦化空間分布和地質(zhì)體形態(tài)的關系。

本次研究采用空間距離分析方法進行距離場分析，采用三維形態(tài)分析方法進行地質(zhì)界面的坡度和趨勢-起伏的分析。

2.1距離場分析

距離場是指空間某點因幾何接近程度不同而受到的場源作用的影響程度，其值為點與場源的距離。成礦距離場分析是借助于空間分析中的距離模型(De Smith et al.,2007)，描述地質(zhì)對象對礦床(點)的空間定位與分布控制信息。本次研究將地質(zhì)體間的距離約定為最小距離，距離場的數(shù)學模型為歐氏距離場。

以下是空間點分別到點狀地質(zhì)體、線狀地質(zhì)體、面狀地質(zhì)體的距離場數(shù)學模型。

2.1.1點狀地質(zhì)體的成礦距離場

DPG=[(xp-xg)2+(yp-yg)2+(zp-zg)2]1/2

(1)

式(1)中，DPG代表空間某點P(xp,yp,zp)到點狀地質(zhì)對象G(xg,yg,zg)的成礦距離。

2.1.2線狀地質(zhì)體的成礦距離場線狀地質(zhì)對象可看成由n條線段L1，L2，…，Ln組成，首先計算出點到線段的最近距離，分3種情況，采用向量法判斷點到線段垂足的位置(圖2)。

圖2　點到線段的距離示意圖Fig.2　Diagram showing the distance of point to segment

這樣就將點到線段的距離計算轉(zhuǎn)化成點到點的距離計算，用公式(1)可求得。

則空間點P(xp,yp,zp)到線狀地質(zhì)對象G(L1，L2，…，Ln)的最小距離便為空間點P(xp,yp,zp)到構(gòu)成該線狀地質(zhì)對象G(xg,yg,zg)的所有線段的歐氏距離中的最小值，即有：

DPG=min{DPL1,DPL2,DPL3,…,DPLn}

(2)

式(2)中，DPG代表空間某點P(xp,yp,zp)到線狀地質(zhì)對象G(L1，L2，…，Ln)的成礦距離。

2.1.3面狀地質(zhì)體的成礦距離場面狀地質(zhì)對象可看成由n個三角網(wǎng)T1，T2，…，Tn組成。首先要計算出空間點P(xp,yp,zp)到每個三角網(wǎng)的最小距離DPT。

理論上,點到面的最小距離是點到面的垂足點間的距離,但TIN模型各個三角形面并不都在1個平面上,所以需要判斷點到三角形所在平面的垂足是否在三角形內(nèi)。

(1) 垂足在三角形內(nèi)，點到垂足間的距離即為點到面的最小距離，用公式(1)可求得。

(2) 垂足不在三角形內(nèi)，則不能用點到垂足的距離來表示點到面的最小距離，點到三角形面的最小距離必將落在三角形邊上，于是需要求點到3邊的最小距離，用公式(2)可求得。

則空間點到整個面狀地質(zhì)對象的G(T1，T2，…，Tn)最小距離就是該點到所有三角形的距離中的最小值，即有：

DPG=min{DPT1,DPT2,DPT3,…,DPTn}

(3)

式(3)中，DPG代表空間某點P(xp，yp，zp)到面狀地質(zhì)對象G(T1，T2，…，Tn)的成礦距離。

2.2地質(zhì)界面坡度分析

地質(zhì)界面上某點的坡度指的是地質(zhì)界面在該點的陡緩程度，即tanα的值為高程增量h與水平增量l的比值(圖3a)。本次研究中的地質(zhì)界面用TIN模型表示，TIN模型上面某點的坡度，用該三角面片與水平面的夾角大小來表示，該夾角大小在數(shù)值上等于三角形的法向量與水平面法向量的夾角(圖3b)，Δl的法向量vl與Z軸的夾角∠slope的值為Δl的坡度值。

TIN模型的任意△ABC的平面方程：

z=ax+by+c

(4)

其坡度計算公式為：

(5)

圖3　地質(zhì)界面坡度示意圖Fig.3　Diagram showing the slope of geological interface

2.3地質(zhì)界面趨勢-起伏分析

地質(zhì)界面趨勢-起伏分析指的是地質(zhì)界面的趨勢形態(tài)分析以及起伏變化分析。

采用趨勢面分析法進行地質(zhì)界面的趨勢-起伏分析。趨勢面分析就是在空間中已知點的控制下，擬合一個連續(xù)的數(shù)學曲面，并以此研究地質(zhì)變量在區(qū)域上和局部范圍內(nèi)變化規(guī)律的一種統(tǒng)計方法。趨勢面的提取方法有克里金法、多項式擬合法、距離平方反比法等。

采用距離平方反比法提取趨勢面，步驟如下。

(1) 設定一個固定的搜索半徑r，只有位于搜索范圍內(nèi)的三角形頂點參與計算。

(6)

如圖4所示，頂點O的高程擬合值就用以其為圓心，半徑r內(nèi)的a、b、c、d、e、f、g、h點的高程值來估計。

圖4　TIN模型某點高程擬合值計算示意圖Fig.4　Diagram showing calculation of fitted elevation of a point on the TIN model

(3) 用所有點計算出的擬合值構(gòu)建趨勢面的TIN模型。

提取出趨勢面后，構(gòu)建趨勢面分析的形態(tài)結(jié)構(gòu)數(shù)學表達式：

(7)

顯然，

(8)

如圖5所示，εi可以用來描述地質(zhì)界面的起伏情況，若εi>0，表明斷裂面在此處為隆起狀態(tài)；εi<0，表明斷裂面在此處為凹陷狀態(tài)。εi的絕對值越大，表明斷裂面凹陷或隆起的程度越大；εi的絕對值越小，表明斷裂面凹陷或隆起的程度越小。

圖5　地質(zhì)界面起伏示意圖Fig.5　Diagram showing undulation of geological interface

3　阿希金礦控礦因素提取

阿希金礦的成礦主要受控于以下2種礦區(qū)構(gòu)造。

(1) 斷裂構(gòu)造。礦體嚴格受斷裂帶控制，賦存于沿火山管道內(nèi)側(cè)形成的F2弧形張性斷裂中，此斷裂又經(jīng)歷了多期的活動，南段礦體明顯受影響，局部又形成了膨大收縮金礦化相對富集的地段,而且其構(gòu)造變形強度與礦化蝕變也具有一定的相關性。

(2) 火山管道相。阿希礦床金礦體主要賦存于近南北向裂隙式噴發(fā)的火山管道相內(nèi)，具有明顯的巖性控礦特征。火山管道相是成礦的巖性條件，火山熱液又是成礦物質(zhì)的主要載體。火山活動為金礦形成提供了物質(zhì)基礎和熱動力成礦條件。

按照控礦因素分析方法，以阿希金礦區(qū)的礦體三維模型為基礎(圖6)，提取主斷裂面距離場因素(dF)、主斷裂面坡度因素(gF)、主斷裂面趨勢-起伏因素(wrF)、管道相距離場因素(dGV)、管道相趨勢-起伏因素(wrGV)作為阿希礦區(qū)的控礦因素。

圖6　阿希金礦區(qū)礦體三維模型圖Fig.6　Three-dimensional model of ore bodies in the Axi gold deposit

3.1主斷裂面距離場因素(dF)

采用距離場分析方法，先分別建立礦體的TIN模型和主斷裂面(F2)的TIN模型，并將礦體TIN模型轉(zhuǎn)化為柵格模型(即礦化立體單元)，再通過公式(1)、(2)計算每個立體單元質(zhì)心點到每個三角面片的最近距離，最后通過公式(3)取點到所有三角面片最近距離的最小值作為每個立體單元質(zhì)心點到主斷裂面(F2)的最近距離(圖7)。空間的幾何接近程度預示著空間關聯(lián)性的大小，斷裂面的距離場因素能反映距離主斷裂面不同距離的礦化強度分布情況。

圖7　主斷裂面距離場因素可視化圖Fig.7　Visualization of distance field factor of the main faults

圖8　主斷裂面坡度因素可視化圖Fig.8　Visualization of slope factors of the main faults

3.2主斷裂面坡度因素(gF)

采用三維形態(tài)分析方法中的坡度分析，先建立主斷裂面(F2)的TIN模型(圖8)，計算出礦化立體單元到主裂面上距離最近的三角面片，通過公式(4)、(5)求出該三角面片頂點的坡度。斷裂面的坡度能夠局部體現(xiàn)斷裂面的產(chǎn)狀變化以及產(chǎn)狀變化為成礦提供的微觀物理化學環(huán)境，斷裂面坡度的變化影響著礦液的運移和沉淀。該因素能夠模擬反映礦體側(cè)伏的規(guī)律。

3.3主斷裂面趨勢-起伏因素(wrF)

采用三維形態(tài)分析中的趨勢-起伏分析，先建立主斷裂面(F2)的TIN模型(圖9)，再通過公式(6)提取趨勢面，最后通過公式(8)求出剩余值，并用該數(shù)值來表示起伏情況。斷裂面一般呈舒緩波狀，斷裂面的形態(tài)尤其是隆起和凹陷形態(tài)往往控制著礦體的產(chǎn)出及分布,斷裂面形態(tài)變化的地段有利于礦液的富集、礦化的疊加，有利于礦體的形成,該因素主要揭示斷裂面起伏對礦化的影響。

圖9　主斷裂面趨勢-起伏因素可視化圖Fig.9　Visualization of trend-undulation factor of the main faults

3.4管道相距離場因素(dGV)

采用距離場分析方法，先分別建立礦體的TIN模型和火山管道相的TIN模型，并將礦體TIN模型轉(zhuǎn)化為柵格模型(即礦化立體單元)，再通過公式(1)、(2)計算每個立體單元質(zhì)心點到每個三角面片的最近距離，最后通過公式(3)計算每個立體單元質(zhì)心點到管道相的最近距離。通過該因素研究火山管道相對成礦的影響程度，模擬距離管道相界面不同距離的礦化強度分布情況(圖10)。

圖10　管道相距離場因素可視化圖Fig.10　Visualization of distance field factors of volcanic conduit facies

3.5管道相趨勢-起伏因素(wrGV)

采用三維形態(tài)分析中的趨勢-起伏分析，先建立火山管道相的TIN模型，再用式(6)提取趨勢面，最后式(8)求出剩余值，并用該數(shù)值表示起伏情況。火山管道相形態(tài)的外凸內(nèi)凹變化能反映出不同程度的地質(zhì)作用，對礦體的形成有一定的影響,通過該因素可以解釋管道相界面的隆起和凹陷形態(tài)對周圍地質(zhì)地質(zhì)空間的控礦作用影響(圖11)。

圖11　管道相趨勢-起伏因素可視化圖Fig.11　Visualization of trend-undulation factor of the volcanic conduit facies

4　討　論

采用空間相關分析法來分析控礦因素變量與礦化指標之間的關系。依據(jù)已獲得的單元控礦地質(zhì)因素變量和單元礦化指標，生成反映單元控礦因素變量與單元礦化指標關系的散點圖,通過散點的疏密程度和變化趨勢分析關聯(lián)關系。

4.1主斷裂面控礦因素

主斷裂面控礦因素與礦化指標的散點圖如圖12所示。

從dF-Au散點圖(圖12a)可知，單元金平均Au品位與因素dF之間存在著一定的關聯(lián)關系，表現(xiàn)為主要礦化富集空間距斷裂面直線距離為0～150 m，大體趨勢為單元到主裂面的距離越近，其礦化指標Au品位越高，且在0～30 m的空間范圍內(nèi)出現(xiàn)最高峰，距主裂面的直線距離大于150 m的單元為無礦或貧礦空間。

從gF-Au散點圖(圖12b)可知，單元金平均Au品位與因素gF之間存在著一定的關聯(lián)關系，表現(xiàn)為主要礦化富集空間附近的斷裂面坡度為55°～90°，2個礦化指標峰值分別出現(xiàn)在70°和80°附近，隨著斷裂面坡度的逐漸減小，礦化指標的值逐漸降低，在小于55°的范圍內(nèi)，礦化指標的值明顯偏低，這些單元為無礦或貧礦空間。

從wrF-Au散點圖(圖12c)可知，單元金平均Au品位與因素wrF之間存在著一定的關聯(lián)關系，表現(xiàn)為主要礦化富集空間附近斷裂面的起伏程度較小，其范圍為-20～20 m，礦化指標出現(xiàn)了2個峰值：一個是在起伏形態(tài)較為平坦的范圍內(nèi)，起伏程度約為0 m；另一個是在稍微外凸的范圍內(nèi)，起伏程度約為外凸15 m處。隨著斷裂面內(nèi)凹以及外凸程度的逐漸增大(>15 m)，礦化指標逐漸降低。在斷裂面起伏程度<20 m及>20 m處，礦化指標明顯偏低，即這些單元為無礦或貧礦空間。

圖12　主斷裂面控礦因素與礦化指標關系散點圖Fig.12　Scatter diagrams showing the relationship between ore-controlling factors and mineralization index of the main faults(a) scatter diagram of dF-Au; (b) scatter diagram of gF-Au; (c) scatter diagram of wrF-Au

綜合分析后，發(fā)現(xiàn)礦體在斷裂構(gòu)造中就位有以下規(guī)律：在距離斷層0～30 m范圍，坡度較大且起伏程度較小的外凸部位，礦體就位最有利。匯流擴容空間是熱液成因礦床礦體就位的重要條件。其實，匯流擴容空間的形成是地質(zhì)因素耦合作用形成的，其不僅與區(qū)域應力場、局部應力場相關，而且熱力場、流體的循環(huán)也發(fā)揮重要作用(劉亮明等，2008)。阿希金礦礦石以石英脈型礦石為主，其成礦作用機制被認為是沸騰作用(翟偉等，2010)。礦區(qū)主礦體均賦存于斷層上盤，斷層面與下部礦體或火山巖截然的界面，被認為是早期巖石的破裂面和斷層后期發(fā)生滑動疊加形成的，與主礦體中常見的巖石角礫相符，可認為斷層破碎帶屬匯流擴容空間。研究表明，阿希金礦經(jīng)歷了約40 Ma的成礦時間(魏佳林，2012)，與成礦作用的短期性不同，因而其成礦過程是構(gòu)造活動引發(fā)熱液上升、沸騰、充填的循環(huán)過程(馬潤則等，1993)。

循環(huán)的大氣降水是成礦流體的主要來源，長期的熱液作用，形成了以礦體為中心的蝕變分帶，而F2斷層存在多期活動，故礦體的就位部位在距離斷層面相近的地段。F2斷層整體呈“S”形展布，在礦區(qū)南、北段的擴張空間控制了礦體的規(guī)模和產(chǎn)出，原因是其局部應力場的不均一性影響了成礦流體的遷移和匯聚，因此礦質(zhì)沉淀部位為斷層外凸部位。而斷層坡度較大的部位，可能反映了原始巖漿上侵通道的位置，是次級斷裂的相連地段，因而對成礦有利度更高。

4.2管道相控礦因素

圖13為管道相控礦因素與礦化指標的散點圖。

圖13　管道相控礦因素與礦化指標關系散點圖Fig.13　Scatter diagrams showing the relationship between ore-controlling factors and mineralization index of volcanic conduit facies(a) scatter diagram of dGV-Au; (b) scatter diagram of wrGV-Au

從dGV-Au散點圖(圖13a)可知，單元金平均Au品位與因素dGV之間存在著一定的關聯(lián)關系，表現(xiàn)為：主要礦化富集空間距管道相直線距離為0～150 m，礦化指標在50 m左右出現(xiàn)最高峰，隨著距管道相直線距離的增大和減小，礦化指標逐漸降低，在大于150 m的單元內(nèi)，礦化指標為無礦或貧礦空間。

從wrGV-Au散點圖(圖13b)可知，單元金平均Au品位與因素wrGv之間存在著一定的關聯(lián)關系，表現(xiàn)為：主要礦化富集空間附近斷裂面的起伏程度較小，起伏程度在-20～20 m范圍，礦化指標峰值出現(xiàn)在斷裂面整體起伏形態(tài)較為平坦的空間范圍內(nèi)；起伏程度約為0 m處，隨著管道相內(nèi)凹及外凸程度的增大，礦化指標值逐漸降低，且外凸部位的礦化富集程度稍強于內(nèi)凹部位；在管道相起伏程度<20 m及>20 m處，礦化指標明顯偏低，即這些單元為無礦或貧礦空間。

綜合分析后，發(fā)現(xiàn)礦體在管道相巖石中距離火山管道相底板約50 m的外凸部位是金成礦的理想部位。阿希金礦管道相巖石由英安質(zhì)角礫熔、含集塊角礫熔巖組成，局部可觀察到流動構(gòu)造，是金礦體的直接圍巖。火山管道相因斷裂和巖漿的頻繁活動，局部應力場和熱力場不穩(wěn)定，成礦作用常發(fā)生在這種物理化學條件急劇變化的部位。F2斷層產(chǎn)于火山管道相西側(cè)邊緣位置，火山管道相在靠近斷層的位置應力場變化明顯，巖石破碎程度最高，F(xiàn)2可能是成礦熱液的通道，故在距離火山管道相底板不遠的位置礦化品位最高。圖13中Au品位與距離顯示了良好的對稱分布，暗示在距管道相底板50 m附近的位置可能是礦區(qū)擴容匯流空間的中心部位。眾所周知，金在熱液中常以絡合物的形式遷移,成礦熱液在循環(huán)過程中，在凹凸相連之連續(xù)起伏部位構(gòu)造應力場變化最強烈，成礦流體的運移速率或運移方向時常變化，導致礦質(zhì)沉淀，因而礦體就位于該地段。

5　結(jié)　論

(1) 采用三維空間分析方法提取了阿希金礦主斷裂面距離場因素、主斷裂面坡度因素、主斷裂面趨勢-起伏因素、管道相距離場因素、管道相趨勢-起伏因素。

(2) 控礦因素和礦化指標關聯(lián)分析發(fā)現(xiàn)，礦化指標和控礦因素具有顯著的相關性，成礦有利部位在距離斷層0～30 m之間、坡度較大、起伏程度較小的外凸部位以及距離火山管道相底板約50 m的外凸部位。

陳建平，呂鵬，吳文，等，2007. 基于三維可視化技術的隱伏礦體預測[J]. 地學前緣，14(5): 54-62.

董連慧，沙德銘,2005. 西天山地區(qū)晚古生代淺成低溫熱液金礦床[M]. 北京:地質(zhì)出版社.

丁建華，肖克炎，婁德波,等，2009. 大比例尺三維礦產(chǎn)預測[J]. 地質(zhì)與勘探，45(6): 729-734.

劉亮明，疏志明，趙崇斌，等，2008. 矽卡巖礦床的匯流擴容空間控礦機制及其對深部找礦的意義: 以銅陵—安慶地區(qū)為例[J]. 巖石學報，24(8): 1848-1856.

龍靈利，高俊，錢青,等，2008. 西天山伊犁地區(qū)石炭紀火山巖地球化學特征及構(gòu)造環(huán)境[J].巖石學報，24(4): 699-710.

毛先成，陳國珖，1988a. 香花嶺錫礦床數(shù)學模型及立體定量預測初探[J]. 地質(zhì)與勘探，24(10): 25-31.

毛先成，陳國珖，1988b. 香花嶺錫礦田隱伏礦床的立體定量預測[J]. 桂林工學院學報，8(1): 15-22.

馬潤則，王潤民，1993. 新疆阿希金礦近礦蝕變巖的研究[J]. 礦物巖石,13(4):57-67.

毛先成，2006. 三維數(shù)字礦床與隱伏礦體立體定量預測研究[D]. 長沙：中南大學.

毛先成，鄒艷紅，陳進，等，2010. 危機礦山深部、邊部隱伏礦體的三維可視化預測：以安徽銅陵鳳凰山礦田為例[J]. 地質(zhì)通報，29(2/3): 401-413.

王強，趙振華，許繼峰，等，2006. 天山北部石炭紀埃達克巖-高鎂安山巖-富Nb島弧玄武質(zhì)巖:對中亞造山帶顯生宙地殼增生與銅金成礦的意義[J]. 巖石學報，22(1): 11-30.

魏佳林，2012. 新疆阿希金礦床礦體空間定位規(guī)律及找礦方向研究[D]. 武漢：中國地質(zhì)大學(武漢).

王慶峰，2013. 新疆阿希金礦深部成礦預測研究[J]. 新疆有色金屬,36(2):33-36.

肖龍，2002. 新疆伊犁京希—伊爾曼得金礦區(qū)角礫巖特征及成因[J]. 礦物巖石，22(2): 9-12.

肖克炎,李楠,孫莉,等，2012. 基于三維信息技術大比例尺三維立體礦產(chǎn)預測方法及途徑[J]. 地質(zhì)學刊，36(3): 229-236.

徐伯駿，曹新志，魏佳林，等，2014. 新疆伊犁阿希—塔吾爾別克—阿庇因迪成礦區(qū)金-鉛鋅成礦系列和成礦模型研究[J]. 地質(zhì)找礦論叢，29(4): 495-505.

趙鵬大,李紫金,胡光道，1992. 重點成礦區(qū)三維立體礦床統(tǒng)計預測：以安徽月山地區(qū)為例[M]. 武漢:中國地質(zhì)大學出版社.

朱永峰，周晶，郭璇，2006a. 西天山石炭紀火山巖巖石學及Sr-Nd同位素地球化學研究[J]. 巖石學報，22(5): 1341-1350.

朱永峰，周晶，宋彪，等，2006b. 新疆“大哈拉軍山組”火山巖的形成時代問題及其解體方案[J]. 中國地質(zhì)，33(3): 487-497.

翟偉，孫曉明，賀小平，等，2007. 新疆阿希低硫型金礦床流體地球化學特征與成礦機制[J]. 地質(zhì)學報，81(5): 659-669.

周圣華，胡慶雯，田培仁，2008. 論新疆伊犁淺成低溫熱型金礦與斑巖型銅(鉬金)礦的相隨相伴規(guī)律[J]. 礦產(chǎn)與地質(zhì)，22(5): 400-404.翟偉，孫曉明，蘇麗薇,等，2010. 新疆阿希金礦:古生代的低硫型淺成低溫熱液金礦床[J]. 地學前緣，17(2):266-285.張耀選,王樂民,王立社,等,2015. 新疆伊犁吐拉蘇盆地西北部金礦床成因分析[J]. 吉林大學學報(地球科學版),45(增刊1): 646-647.

FALLARA F, LEGAULT M, RABEAU O, 2006. 3D integrated geological modeling in the Abitibi Subprovince (Quebec, Canada): techniques and applications[J]. Exploration and Mining Geology, 15(1/2): 27-43.

DE SMITH M J, GOODCHILD M F, LONGLEY P A, 2007. Geospatial Analysis: A comprehensive Guide to Principles, Techniques and Software Tools[M]. Leicester, GBR: Troubador Publishing Ltd.

XIA L Q, XIA Z C, XU X Y, et al., 2008. Relative contributions of crust and mantle to the generation of the Tianshan Carboniferous rift-related basic lavas, northwestern China[J].Journal of Asian Earth Sciences, 31(4/6): 357-378.

ZHU Y F, ZHANG L F, GU L B, et al., 2005. The zircon SHRIMP chronology and trace element geochemistry of the Carboniferous volcanic rocks in western Tianshan Mountains [J]. Chinese Science Bulletin, 50(19): 2201-2212.

Three-dimensional quantitative analysis and extraction of ore-controlling factors of Axi gold deposit in Xinjiang

LIN Yuting1,2, MAO Xiancheng1,2, LIU Xiaoxia1,2, REN Jia1,2, LIU Zhankun1,2, PAN Min1,2, LUO Dan1,2, DENG Hao1,2, ZHU Jiawei1,2, ZHANG Miaomiao1,2, XIA Fang3, XIAO Fei3, WEI Qingfeng3, HAN Jianmin3, FAN Hongxi3, RAO Wentian3, HUANG Youyou3

(1. Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring, Ministry of Education, Changsha 410083, Hunan, China; 2. School of Geosciences and Info-Physics, Central South University, Changsha 410083, Hunan, China； 3. Xinjiang Nonferrous Metals Industry Group, Urumqi 830000, Xinjiang, China)

Axi gold deposit in Xinjiang began its production in the 1990s. Till now, its surface and shallow ores have been almost exhausted. It is an effective way of discovering deep-seated concealed ore bodies to extract ore-controlling factors by 3D spatial analysis method and then to make metallogenic prognosis. This work used 3D modelling and adopted distance field analysis, slope analysis of geological interface and trend-undulation analysis of geological interface to extract five ore-controlling factors, including distance field factor of main faults, slope factor of main faults, trend-undulation factor of main faults, distance field factor of volcanic conduit facies and trend-undulation factor of volcanic conduit facies. In addition, we analyzed the relationship between mineralization index and ore-controlling factors by scatter diagrams, got the mineralization enrichment space of each factor and extracted the favorable mineralization location. This conclusion can provide information support for 3D metallogenic prediction of deep concealed ore bodies in the Axi gold deposit.Key words： ore-controlling factors; distance field analysis; slope analysis; trend-undulation analysis; 3D quantitative analysis; deep-seated concealed ore bodies; Axi gold deposit; western Tianshan Mountains in Xinjiang

10.3969/j.issn.1674-3636.2016.03.477

2016-06-21；

2016-07-02；編輯：陸李萍

國家自然科學基金項目(41472301),中南大學“創(chuàng)新驅(qū)動計劃”項目(2015CX008),新疆維吾爾自治區(qū)高層次人才工程項目

林玉婷(1993—),女,碩士研究生,研究方向為GIS、地學三維建模與隱伏礦體三維預測,E-mail: linyuting93@126.com

P628； P618.51

A

1674-3636(2016)03-0477-09