綜合找礦方法在河臺金礦找礦預(yù)測中的應(yīng)用

朱江建，陳廣浩，曾喬松，龔貴倫，4，易金，龔朝陽，黃棟林，林舸

(1.中國科學(xué)院 廣州地球化學(xué)研究所 礦物學(xué)與成礦學(xué)重點實驗室，廣東 廣州510640;2.中國科學(xué)院研究生院，北京100049;3.中國科學(xué)院南海海洋研究所邊緣海地質(zhì)重點實驗室，廣東廣州510301;4.中國人民武裝警察部隊黃金第十一支隊，湖南寧鄉(xiāng)410600;5.廣東高要河臺金礦，廣東肇慶526127)

綜合找礦方法在河臺金礦找礦預(yù)測中的應(yīng)用

朱江建1，2，陳廣浩3，曾喬松1，龔貴倫1，2，4，易金1，2，龔朝陽5，黃棟林5，林舸1

(1.中國科學(xué)院 廣州地球化學(xué)研究所 礦物學(xué)與成礦學(xué)重點實驗室，廣東 廣州510640;2.中國科學(xué)院研究生院，北京100049;3.中國科學(xué)院南海海洋研究所邊緣海地質(zhì)重點實驗室，廣東廣州510301;4.中國人民武裝警察部隊黃金第十一支隊，湖南寧鄉(xiāng)410600;5.廣東高要河臺金礦，廣東肇慶526127)

河臺金礦是一個典型的與韌性剪切帶有關(guān)的金礦床，也是目前粵西、桂東南已發(fā)現(xiàn)的最大的金礦床。為了給深邊部的找礦提供科學(xué)依據(jù)，本文應(yīng)用伽馬能譜、地電化學(xué)及數(shù)值模擬產(chǎn)生的體積應(yīng)變異常帶對本區(qū)進(jìn)行了找礦預(yù)測實驗研究。伽馬能譜K異常帶指示了本區(qū)糜棱巖化及熱液蝕變過程中引起的鉀化蝕變作用，金礦體與糜棱巖化帶主要分布在K異常帶中;地電化學(xué)異常帶可能代表了本區(qū)不同深度礦體在地表的垂直投影;體積應(yīng)變異常帶及其附近有利于發(fā)展為本區(qū)的含礦糜棱巖化帶。在同一構(gòu)造成礦體系中，根據(jù)相似類比原則，表明以上三種方法在本區(qū)進(jìn)行找礦預(yù)測可行。綜合以上三種方法，圈定了3個可能的含礦帶，建議對其進(jìn)行驗證。

河臺金礦;數(shù)值模擬;伽馬能譜;地電化學(xué);體積應(yīng)變;找礦預(yù)測

0 引 言

粵西河臺金礦是典型的與韌性剪切帶有關(guān)的金礦床(段嘉瑞等，1992a)。含金韌性剪切帶是廣寧－羅定深層次推覆構(gòu)造前緣的一個陡傾帶(張伯友等，1995)，位于華南云開大山北翼、羅定－廣寧和吳川－四會斷裂帶交匯部位。韌性剪切的動力來源于印支板塊與華南板塊三疊紀(jì)碰撞效應(yīng)向華南內(nèi)陸的延伸(圖1左上角)(Zhang and Cai，2009;Cai and Zhang，2009)。鐘增球等(1997)指出它的成礦模式與Bonnemaison和Marcoux(1990)總結(jié)的剪切帶三階段金礦化模式相仿。糜棱巖化過程中的成礦溫度主要集中于300～340℃，成礦時代主要集中于197～168 Ma(朱江建等，2011a)。然而，糜棱巖化可能并未形成達(dá)到工業(yè)品位的金礦體，其后的熱液成礦期才是河臺金礦的主要成礦期(Zhang et al.，2001;劉偉，2004;翟偉等，2006;朱江建等，2011a)。由于糜棱巖化帶，尤其是硅化糜棱巖(王斯亮和李得欽，2000;伍思洪，2005)，在后期脆性變形作用下易產(chǎn)生脆性破裂，所以熱液成礦作用往往疊加于糜棱巖化帶上。

鑒于糜棱巖化帶與礦體的密切關(guān)系及目前的找礦勘查需要，在我們應(yīng)用地電化學(xué)測量(夏永健等，2008;朱江建等，2011b)、EH4測量(王彪等，2009)等單一找礦方法對本區(qū)進(jìn)行找礦研究后，本文通過對河臺云西礦床開展伽馬能譜、地電化學(xué)和數(shù)值模擬的綜合研究，尋找糜棱巖化帶和識別可能的礦化異常信息，進(jìn)而更好地服務(wù)于礦山深邊部找礦勘探。

1 礦床地質(zhì)特征

河臺金礦由北往南依次出露震旦系、奧陶系、志留系。震旦系是本區(qū)主要的含礦地層，以局部混合巖化的石英云母片巖、石英巖為主;奧陶系與志留系以薄層淺變質(zhì)砂巖、粉砂巖及薄層板巖為主。F1斷裂是礦區(qū)主要的導(dǎo)礦構(gòu)造之一，位于礦區(qū)南部，傾向NW，傾角55°～70°。含礦構(gòu)造為糜棱巖化帶及發(fā)育于其中的脆性斷裂(劉偉等，2006)，其中糜棱巖化帶傾向NW，局部反轉(zhuǎn)，傾角60°～85°(李新福和湛培任，2007)。導(dǎo)礦構(gòu)造(F1)與含礦構(gòu)造在剖面上呈“y”字型(圖1右下角)。礦區(qū)西部出露黑云母斜長花崗巖，東部出露巨斑狀黑云母二長花崗巖(圖1)。

圖1 河臺金礦區(qū)地質(zhì)簡圖(據(jù)陳駿和王鶴年，1993;朱江建等，2011a;蔡建新，2009修改)Fig.1 Simplified geological map of the Hetai goldfield，Guangdong province，South China

糜棱巖中發(fā)育S面理和糜棱巖C面理，變形程度與破碎程度由中心向邊界遞減，即以超糜棱巖為中心，向兩側(cè)漸變?yōu)槊永鈳r、初糜棱巖、糜棱巖化巖石。原巖有變質(zhì)巖、偉晶巖、黑云母斜長花崗巖(礦區(qū)西北部)、混合巖等。韌性剪切帶分布在F1斷裂南北兩側(cè)，呈平行條帶狀，成群成帶地夾在相對剛性的巖塊中間，總寬度約3000 m，其中北側(cè)約2500 m，南側(cè)約500 m。糜棱巖化帶平面形態(tài)大部分為狹長帶狀，長一般500～1000 m，寬數(shù)十厘米至數(shù)十米。根據(jù)圍巖條件和形態(tài)規(guī)模，可分為北帶、中帶和南帶。其中，北帶距F1斷裂1500～2000 m，全長約5500 m，由9、19、11、12、13 等糜棱巖化帶組成，圍巖為混合巖化片巖，是礦區(qū)主要的礦化帶;中帶距F1斷裂500～1500 m，全長大于 9000 m，由 15、5、10、40等糜棱巖化帶組成，圍巖為混合巖，形態(tài)復(fù)雜，有“小螞蝗”之稱;南帶距F1斷裂500 m以內(nèi)，全長大于9000 m，由36、38、58等糜棱巖化帶組成，圍巖以混合巖化片巖、千枚狀粉砂質(zhì)頁巖、變質(zhì)石英粉砂巖為主，部分為混合巖(圖1)。與礦化相關(guān)的蝕變主要有硅化與黃鐵礦化。礦石類型主要有浸染狀硅化千糜巖金礦石和具網(wǎng)脈狀金屬硫化物疊加的硅化千糜巖金礦石。伴生金屬礦物主要包括黃銅礦、黃鐵礦、磁黃鐵礦、菱鐵礦、毒砂;非金屬礦物主要有石英、絹云母;副礦物有電氣石、磷釔石、鋯石、榍石、獨居石等(劉偉等，2006)。

2 含礦糜棱巖的特征及其形成條件推斷

目前所發(fā)現(xiàn)的金礦體基本位于糜棱巖化帶中(51號含礦石英脈除外)，但不是所有的糜棱巖化帶都含礦。能形成工業(yè)礦體的糜棱巖須具備下列5個條件:(1)糜棱巖化帶產(chǎn)于礦區(qū)內(nèi)最主要的礦源層震旦系中(王鶴年等，1989;陸建軍，1993;張伯友等，1993);(2)糜棱巖具有一定的規(guī)模，長度一般在幾百米以上;(3)有沿脈大斷裂貫穿整個糜棱巖化帶，深部要有導(dǎo)礦通道;(4)糜棱巖化帶中有次生張性裂隙導(dǎo)礦和容礦;(5)糜棱巖形成較早，糜棱巖化程度較高(李新福和湛培任，2007)。目前所發(fā)現(xiàn)的金礦體主要分布于9、19號糜棱巖化帶(云西礦床)、11號糜棱巖化帶(高村礦床)、12、13號糜棱巖化帶(后逕礦床)中。這些糜棱巖化帶平面上呈左階斜列產(chǎn)出(圖1)。礦脈厚度從幾厘米到幾十米不等，礦體主要沿糜棱巖的C面理(圖2a)及糜棱巖中的脆性裂隙分布(圖2b)。

圖2 糜棱巖化帶與礦體的空間關(guān)系Fig.2 The spatial relationship between the mylonitized zone and orebodies

蔡建新(2009)通過求解石英的捕獲溫度及Si含量，同時結(jié)合多硅白云母P-T條件與Si含量的經(jīng)驗關(guān)系(Massonne and Schreyer，1987)，推斷河臺金礦剪切變形壓力為420～650 MPa。差異壓力可通過石英亞顆粒進(jìn)行估測(Mercier，1977)，多數(shù)樣品在160～210 MPa;應(yīng)變速率可通過石英動態(tài)重結(jié)晶顆粒與亞晶進(jìn)行估測(Twiss，1977)，其范圍為1.13 ×10－11～3.72 ×10－11s－1(段嘉瑞等，1992b)，一般區(qū)域性應(yīng)變速率為 10－13～10－15s－1(Pfiffner and Ramsay，1982)，說明糜棱巖化帶形成時的應(yīng)變速率較快。Stipp and Heilbronner(2002)、Frederico et al.(2010)認(rèn)為石英發(fā)生動態(tài)重結(jié)晶的起始溫度為300℃，Pryer(1993)認(rèn)為長石發(fā)生動態(tài)重結(jié)晶的起始溫度為500℃。河臺金礦石英發(fā)生動態(tài)重結(jié)晶，而長石基本上只是發(fā)生脆性變形，因此限定糜棱巖化帶形成時的溫度為300～500℃。數(shù)值模擬結(jié)果表明，糜棱巖形成于108°－288°方位的擠壓作用(朱江建等，2011a)。

3 各種異常帶在本區(qū)的指示意義

數(shù)值模擬區(qū)的邊界如圖1所示，伽馬能譜測量在研究區(qū)的39、8、56號勘探線開展;地電化學(xué)測量在研究區(qū)的 55、39、11、8、56 號勘探線開展(圖3)。

3.1 伽馬能譜鉀(K)異常帶

圖3 地電化學(xué)與伽馬能譜測量所在勘探線位置圖Fig.3 The exploration lines for gamma spectrometry and geoelectrical measurement

絹云母是含鉀礦物，可引起伽馬能譜測量的K高異常(潘勇飛，1994)。結(jié)合河臺金礦的地質(zhì)實際，提出本區(qū)的K異常有兩種原因。一是糜棱巖化形成的K異常。主要證據(jù)有:(1)河臺金礦絹云母石英千糜巖的流體包裹體中K+濃度較高，為形成大量絹云母提供了物質(zhì)條件(何文武和張文淮，1993);(2)剪切造成的細(xì)?；纱偈沽黧w向糜棱巖化帶中匯聚(Oliver，2001;Musumeci，2002)，指示了糜棱巖化過程中水巖反應(yīng)的進(jìn)行，而模擬實驗表明水巖反應(yīng)過程中長石向絹云母的轉(zhuǎn)化很容易發(fā)生(Zhao et al.，2001);鏡下觀察也表明河臺金礦糜棱巖化過程中長石向絹云母轉(zhuǎn)化(圖4a、b)是普遍現(xiàn)象(岳石，1990)。鏡下統(tǒng)計表明，初糜棱巖中絹云母含量為5%～25%，糜棱巖中為20%～50%，超糜棱巖中為20%～65%，說明絹云母含量與剪應(yīng)變強(qiáng)度為正相關(guān)關(guān)系。結(jié)合含金糜棱巖的礦化強(qiáng)度與剪應(yīng)變強(qiáng)度也為正相關(guān)關(guān)系(段嘉瑞等，1992a)，及大量絹云母在動力變質(zhì)期的生成有利于熱液成礦期硅化作用的進(jìn)行(何文武和張文淮，1993)，表明糜棱巖化形成的K異常帶與糜棱巖化帶中礦體的位置一致。本區(qū)絹云母的另一種成因是熱液期的絹云母化。雖然，絹云母化與硅化、黃鐵礦化(與礦化密切相關(guān))不是一個期次(劉偉，2004)，但糜棱巖化帶，尤其是硅化糜棱巖(王斯亮和李得欽，2000;伍思洪，2005)，在后期(脆性變形時期)應(yīng)力作用下易于產(chǎn)生脆性破裂，表明熱液期的絹云母化也同樣主要發(fā)育于糜棱巖帶及其附近。這較好地解釋了研究區(qū)伽馬能譜K異常帶的寬度大于9號含礦糜棱巖化帶的寬度(圖5)。綜上所述，伽馬能譜的K異常帶，揭示了本區(qū)糜棱巖化或熱液蝕變過程中引起的鉀化蝕變作用，含金糜棱巖帶主要發(fā)育在K異常帶中。

圖4 糜棱巖中的長石絹云母化Fig.4 Micrographs showing sericitization of feldspar in the mylonites

3.2 地電化學(xué)異常帶

地電化學(xué)找礦實踐表明，地電化學(xué)異常帶代表了不同深度礦體在地表的垂直投影(李曉勇，2002)。與常規(guī)化探方法相比，地電化學(xué)找礦方法的主要優(yōu)點是靈敏度高、預(yù)測深度大(羅先熔和楊曉，1989;譚克仁，2000;文美蘭等，2010;曹中煌等，2010);與物探方法相比，地電化學(xué)找礦方法的主要優(yōu)點是避免了地球物理找礦方法的多解性。據(jù)前人研究(羅先熔，1994;羅先熔等，1997)，地電化學(xué)找礦適用于地表較厚土壤層覆蓋的糜棱巖化蝕變巖型的金礦。測量結(jié)果表明，本區(qū)9號含礦糜棱巖化帶上方或傾向側(cè)上方地電化學(xué)高異常明顯(圖5)，表明地電化學(xué)異常帶揭示了河臺金礦的含金糜棱巖帶的存在與產(chǎn)出位置。值得注意的是，地電化學(xué)在河臺金礦的應(yīng)用，要結(jié)合實際情況，考慮可能的非礦致異常與非地電化學(xué)異常區(qū)可能含礦的情況(朱江建等，2011b)。

3.3 體積應(yīng)變異常帶

構(gòu)造－流體成礦數(shù)值模擬實驗，前人已經(jīng)開展了一些有益的探討(Liu and Peng，2003;Zhang et al.，2007;龔發(fā)雄等，2009;Liu et al.，2010;胡才志等，2010;桑廣森等，2010;趙義來和劉亮明，2011;錢建平等，2011)。地質(zhì)構(gòu)造變形數(shù)值模擬是采用數(shù)學(xué)物理方法對相關(guān)地質(zhì)學(xué)問題的科學(xué)描述，并采用理論分析和數(shù)值模擬實驗方法，對所描述的相關(guān)地質(zhì)學(xué)過程進(jìn)行定量化求解(林舸等，2005)。

本文首次運用模擬產(chǎn)生的體積應(yīng)變異常帶對本區(qū)進(jìn)行找礦預(yù)測研究。鑒于本次數(shù)值模擬所遵從的地質(zhì)依據(jù)(朱江建等，2011a)，認(rèn)為模擬結(jié)果在理論上是可靠的?？紤]到本區(qū)各巖性的實測參數(shù)不能完全代表糜棱巖化帶形成時的參數(shù)及糜棱巖形成過程中相關(guān)的參數(shù)是變化的，而本次模擬的整個過程中相關(guān)參數(shù)保持不變等因素，認(rèn)為體積應(yīng)變異常帶及其附近有利于發(fā)展為本區(qū)的含礦糜棱巖化帶(朱江建等，2011a)。

體積應(yīng)變異常帶與糜棱巖化帶、K異常帶、地電化學(xué)異常帶對比如下:Ⅰ號體積應(yīng)變異常帶與9號糜棱巖化帶的位置、總體方位、寬度與長度都較為一致;Ⅱ號體積應(yīng)變異常帶東半段與19號糜棱巖化帶東半段的位置、方位、規(guī)模與長度都較為一致，西半段與19號糜棱巖帶西半段的方位相差較大;Ⅲ號體積應(yīng)變異常帶處未發(fā)現(xiàn)糜棱巖化帶，但其西半段在地電化學(xué)異常帶與伽馬能譜異常帶內(nèi)，東半段也與伽馬能譜異常帶較為吻合，只是稍偏南;Ⅳ號體積應(yīng)變異常帶處未發(fā)現(xiàn)糜棱巖化帶，但其在地電化學(xué)高異常帶內(nèi)且與伽馬能譜K高異常帶穩(wěn)合(圖5)。綜上所述，體積應(yīng)變異常帶對本區(qū)進(jìn)行找礦預(yù)測具有參考意義。

4 找礦預(yù)測

4.1 找礦預(yù)測可行性分析

根據(jù)39、8、56號勘探線的伽馬能譜測量結(jié)果(圖3)，在研究區(qū)圈定了三條主要的K異常帶(圖5);根據(jù) 55、39、11、8、56 號勘探線的地電化學(xué)測量結(jié)果(圖3)，在研究區(qū)圈定了三條地電化學(xué)異常帶(圖5);3條體積應(yīng)變異常帶也標(biāo)示在圖5中。

9號糜棱巖帶是本區(qū)已開采的含金礦糜棱巖帶(已從160 m中段開采至地下－140 m中段)，而19號糜棱巖化帶的含礦性(尤其是深部的)還未確定?？紤]到19號糜棱巖化帶與本區(qū)主要的含礦糜棱巖化帶(如9號糜棱巖化帶)的方位相差較大(圖5)，且與K、地電化學(xué)、體積應(yīng)變異常帶的方位相差也較大(K、地電化學(xué)、體積應(yīng)變異常帶的方位與9號糜棱巖化帶的一致)，說明19號糜棱巖化帶真正的延伸方位還待商榷(19號糜棱巖化帶可能由多條尖滅側(cè)現(xiàn)的糜棱巖化帶組成)。鑒此，本次主要通過與9號糜棱巖化帶的對比研究，探討K、地電化學(xué)、體積應(yīng)變異常帶在本區(qū)找礦預(yù)測的可行性。

9號糜棱巖化帶及其附近如上三種異常帶的空間位置與9號糜棱巖化帶的空間位置較為吻合;9號糜棱巖化帶圍巖的地質(zhì)特征(主要為云母石英片巖)與該帶以外的三種異常帶的地質(zhì)特征較為一致;9號糜棱巖化帶及其附近三種異常帶的特征(方位與規(guī)模)分別與該帶以外同種異常帶的特征較為一致(圖5)。在同一構(gòu)造成礦體系中，根據(jù)相似類比原則，說明了以上三種方法在本區(qū)進(jìn)行找礦預(yù)測的可行性。

圖5 K、地電化學(xué)、體積應(yīng)變異常帶與預(yù)測的含礦帶Fig.5 Anomalous belts of K，geo-electrochemical value，volumetric strain increment and the forecasting ore-bearing zones

4.2 找礦預(yù)測可靠性分析

由K、地電化學(xué)、體積應(yīng)變異常帶的對比研究，同時結(jié)合本區(qū)的地質(zhì)實際，在本區(qū)圈定了3個可能的含礦異常帶。

Ⅰ號異常帶長約855 m，最寬處約85 m。主要圈定依據(jù)是其總體方位(68°)與礦區(qū)主要的含礦糜棱巖化帶一致(如9號糜棱巖化帶的總體方位約為71°)，也與19號糜棱巖化帶0線以東的位置與方位一致;該異常帶39線在－140 m巷道處發(fā)現(xiàn)了糜棱巖及黃鐵礦化，雖然金品位不高(0.01 g/t)，但提供了深部找礦線索。地表出現(xiàn)了較多的小糜棱化帶也預(yù)示深部可能有較大的糜棱巖化帶;該異常帶55、39、8線有地電化學(xué)異常;56線與39線之間有體積應(yīng)變異常帶(圖5)。

Ⅱ號異常帶長約1136 m，最寬處約72 m。主要圈定依據(jù)是其總體方位(68°)與礦區(qū)主要的含礦糜棱巖化帶一致;11號糜棱巖化帶(河臺金礦主要的含礦糜棱巖化帶)沿走向延伸正好位于Ⅱ號異常帶內(nèi)，說明該異常帶可能是11號糜棱巖化帶西部的延伸;25線以東，19號糜棱巖化帶已基本在Ⅱ號異常帶之外，較好地解釋了已劃定的19號糜棱巖化帶從35線開采至25線礦體變貧難以繼續(xù)開采(可能沿Ⅱ號異常帶開采效果會更好);Ⅱ號異常帶39、8、56線都有伽馬能譜異常;55、39、11、8線有地電化學(xué)異常;Ⅲ號體積應(yīng)變異常帶有近一半位于該異常帶中(圖5)。

Ⅲ號異常帶長約1146 m，最寬處約65 m。主要圈定依據(jù)是其總體方位(68°)與礦區(qū)主要的含礦糜棱巖化帶一致;目前11號糜棱巖化帶北部約100 m處發(fā)現(xiàn)了12號含礦糜棱巖化帶(－230中段128線)。而12號糜棱巖化帶向西沿走向延伸則正好位于Ⅲ號伽馬能譜異常帶。所以同樣認(rèn)為，Ⅲ號伽馬能譜異常帶可能為12號糜棱巖化帶向西延伸的結(jié)果;Ⅲ號異常帶39、8、56線有伽馬能譜異常;55、39、11、8、56線有明顯的地電化學(xué)異常;Ⅳ號體積應(yīng)變異常帶全部位于該異常帶中(圖5)。

綜上所述，圈定的3個可能的含礦異常帶都有地質(zhì)依據(jù)，K、地電化學(xué)、體積應(yīng)變異常帶與這3個異常帶都有較好的對應(yīng)關(guān)系，建議對其進(jìn)行驗證。

5 結(jié) 論

(1)伽馬能譜K異常帶指示了本區(qū)糜棱巖化及熱液蝕變過程引起的鉀化蝕變作用，金礦體與糜棱巖化帶主要分布在鉀化蝕變帶中;地電化學(xué)異常帶可能代表了本區(qū)不同深度礦體在地表的垂直投影;體積應(yīng)變異常帶及其附近有利于發(fā)展為本區(qū)的含礦糜棱巖化帶。

(2)9號糜棱巖化帶及其附近如上三種異常帶的空間位置與9號糜棱巖化帶的空間位置較為吻合;9號糜棱巖化帶圍巖的地質(zhì)特征(主要為云母石英片巖)與該帶以外的三種異常帶的地質(zhì)特征較為一致;9號糜棱巖化帶及其附近三種異常帶的特征(方位與規(guī)模)分別與該帶以外同種異常帶的特征較為一致。在同一構(gòu)造成礦體系中，根據(jù)相似類比原則，說明以上三種方法在本區(qū)進(jìn)行找礦預(yù)測可行。結(jié)合以上三種方法，圈定了3個可能的含礦異常帶，建議對其進(jìn)行驗證。

致謝:感謝兩位審稿老師及編輯部對本文提出的修改意見。

蔡建新.2009.華南不規(guī)則大陸邊緣三疊紀(jì)楔合和復(fù)理石盆地充填.北京:中國科學(xué)院研究生院博士學(xué)位論文.

曹中煌，羅先熔，王培培，王發(fā)明.2010.不同覆蓋區(qū)地電化學(xué)提取法尋找銅鎳礦對比研究.地質(zhì)與勘探，46(3):476－482.

陳駿，王鶴年.1993.廣東省河臺含金剪切帶中REE及其它微量元素的含量和分布特征.礦床地質(zhì)，12(3):202－211.

段嘉瑞，何紹勛，周崇智.1992a.剪切帶型金礦——以廣東河臺金礦為例.地質(zhì)與勘探，28(6):6－11.

段嘉瑞，何紹勛，周崇智，徐靄君.1992b.廣東河臺金礦構(gòu)造特征及剪切帶型金礦模式探討.中南礦冶學(xué)院學(xué)報，23(3):245－253.

龔發(fā)雄，單業(yè)華，李自安，劉士林.2009.松遼盆地中部深部斷層影響淺部次級斷層發(fā)育的數(shù)值研究.大地構(gòu)造與成礦學(xué)，33(3):396－401.

何文武，張文淮.1993.廣東河臺金礦成礦物理－化學(xué)條件及找礦方向.礦床地質(zhì)，12(2):120－128.

胡才志，曾佐勛，任建業(yè)，祁鵬，李志勇，胡燁.2010.黃驊坳陷新生代沉降中心遷移:數(shù)值模擬證據(jù).大地構(gòu)造與成礦學(xué)，34(4):573－584.

李曉勇.2002.地電提取測量法尋找隱伏礦方法研究及應(yīng)用.大地構(gòu)造與成礦學(xué)，26(1)，92－96.

李新福，湛培任.2007.河臺金礦田糜棱巖帶含礦性分析.中國礦山工程，36(1):13－16.

林舸，趙重斌，張晏華，王岳軍，單業(yè)華.2005.地質(zhì)構(gòu)造變形數(shù)值模擬研究的原理、方法及相關(guān)進(jìn)展.地球科學(xué)進(jìn)展，20(5):549－555.

劉偉.2004.廣東河臺金礦礦床成因及成礦規(guī)律研究.長沙:中南大學(xué)碩士學(xué)位論文.

劉偉，戴塔根，黃滿湘，賀輝.2006.廣東河臺金礦礦體賦存規(guī)律及找礦前景.黃金，27(3):9－13.

陸建軍.1993.河臺韌性剪切帶型金礦床成礦作用地球化學(xué)研究.南京大學(xué)學(xué)報，29(2):293－302.

羅先熔.1994.再論地電化學(xué)測量法尋找隱伏礦床.桂林冶金地質(zhì)學(xué)院學(xué)報，14(3):295－302.

羅先熔，王衛(wèi)民，張佩華.1997.隱伏金礦地電化學(xué)異常形成機(jī)制及異常形態(tài)特征.有色金屬礦產(chǎn)與勘查，6(6):364－367.

羅先熔，楊曉.1989.地電化學(xué)測量找尋隱伏礦床的研究及找礦預(yù)測.地質(zhì)與勘探，25(12):43－51.

潘勇飛.1994.伽馬能譜法找金礦.物探與化探，18(6):444－456.

錢建平，陳宏毅，吳小雷，王自國，蒙勇.2011.膠東望兒山金礦成礦構(gòu)造分析和成礦預(yù)測.大地構(gòu)造與成礦學(xué)，35(2):221－231.

桑廣森，夏斌，張勝利，蔡周榮，梁正中，萬志峰.2010.松遼盆地徐家圍子三維構(gòu)造應(yīng)力場數(shù)值模擬研究.大地構(gòu)造與成礦學(xué)，34(2):196－203.

譚克仁.2000.金礦地電化學(xué)勘查新技術(shù)、新方法研究進(jìn)展.黃金科學(xué)技術(shù)，8(1):23－31.

王彪，任廣利，王核，林舸，劉建平，朱江建.2009.EH4連續(xù)電導(dǎo)率法在廣東省河臺金礦的應(yīng)用.礦物學(xué)報(增刊):562－563.

王鶴年，張景榮，戴愛華，凌井生，陳礎(chǔ)庭，季明鈞.1989.廣東河臺糜棱巖帶蝕變巖型金礦床的地球化學(xué)研究.礦床地質(zhì)，8(2):61－71.

王斯亮，李得欽.2000.河臺金礦云西礦床富礦包特征及預(yù)測.黃金學(xué)報，2(1):17－21.

文美蘭，羅先熔，熊健，曾南石，侯寶宏，John Keelin.2010.地電化學(xué)法在南澳大利亞尋找隱伏金礦的研究.地質(zhì)與勘探，46(1):153－159.

伍思洪.2005.河臺金礦田富礦包的特征.黃金科學(xué)技術(shù)，13(3):33－35.

夏永健，陳友東，王斯亮，龔朝陽，龔貴倫，張德圣，陳廣浩，林舸.2008.地電化學(xué)測量在河臺金礦找礦中的試驗研究.黃金，29(9):7－10.

岳石.1990.構(gòu)造動力成巖成礦作用的實驗研究.大地構(gòu)造與成礦學(xué)，14(4):325－332.

張伯友，俞鴻年，郭令智，董平.1993.粵西深層次推覆構(gòu)造帶的成巖成礦規(guī)律.地質(zhì)找礦論叢，8(3):1－13.

張伯友，趙振華，夏斌，石滿全，俞鴻年，郭令智.1995.廣寧深層次推覆構(gòu)造前緣陡傾帶變形構(gòu)造研究.地震地質(zhì)，17(4):405－415.

翟偉，李兆麟，孫曉明，黃棟林，梁金龍，苗來成.2006.粵西河臺金礦鋯石SHRIMP年齡及其地質(zhì)意義.地質(zhì)論評，52(5):690－699.

趙義來，劉亮明.2011.復(fù)雜形態(tài)巖體接觸帶成礦耦合動力學(xué)三維數(shù)值模擬:以安慶銅礦為例.大地構(gòu)造與成礦學(xué)，35(1)，128－136.

鐘增球，周漢文，游振東.1997.云開隆起區(qū)剪切帶陣列及其與金礦成礦關(guān)系.地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報，22(1):20－26.

朱江建，陳廣浩，龔貴倫，易金，蔡建新，龔朝陽，黃棟林，曾喬松，林舸.2011a.廣東河臺金礦糜棱巖化過程構(gòu)造－流體成礦研究.地學(xué)前緣，18(5):67－77.

朱江建，曾喬松，易金，龔朝陽，李新福，黃棟林，王斯亮，夏永健，龔貴倫，陳廣浩，林舸.2011b.地電化學(xué)測量在河臺金礦找礦預(yù)測中的應(yīng)用.地質(zhì)與勘探，47(5)，894－902.

Bonnemaison M and Marcoux E.1990.Auriferous mineralization in some shear zones:A three-stage model of metallogenesis.Mineralium Deposita，25(2):96－104.

Cai J X and Zhang K J.2009.A new model for the Indochina and South China collision during the Late Permian to the Middle Triassic.Tectonophysics，467:35－43.

Frederico M F，Ginaldo A C，Rosa M S B and Kazuo F.2010.Quartz recrystallization regimes，c-axis texture transitions and fluid inclusion reequilibration in a prograde greenschist to amphibolite facies mylonite zone(Ribeira Shear Zone，SE Brazil).Tectonophysics，485:193－214.

Liu L M and Peng S L.2003.Prediction of hidden ore bodies by synthesis of geological，geophysical and geochemical information based on dynamic model in Fenghuangshan ore field，Tongling district，China.Journal of Geochemical Exploration，81(1－3):81－98.

Liu L M，Zhao Y L and Zhao C B.2010.Coupled geodynamics in the formation of Cu skarn deposits in the Tongling-Anqing district，China:Computational modeling and implications for exploration.Journal of Geochemical of Exploration，106:146－155.

Massonne H J and Schreyer W.1987.Phengite geobarometry based on the limiting assemblage with K-feldspar，phlogopite，and quartz.Contrib Mineral Petrol，96:212－224.

Mercier R J.1977.Stress in the lithosphere:inferences from steadystate flow of rocks.Pure and Applied Geophysics，115:199－226.

Musumeci G.2002.Sillimanite-bearing shear zones in syntectonic leucogranite:Fluid assisted brittle ductile deformation under amphibolites facies conditions.Journal of Structural Geology，24(9):1491－1505.

Oliver N H S.2001.Linking of regional and local hydrothermal systems in the mid-crust by shearing and faulting.Tectonophysics，335:147－ 161.

Pfiffner O A and Ramsay J G.1982.Constraints on geological strain rates:Arguments from finite strain states of naturally deformed rocks.Journal of Geophysical Research，87:311－321.

Pryer L L.1993.Microstructures in feldspars from a major crustal thrust zone:The Grenville Front，Ontario，Canada.Journal of Structural Geology，15(1):21－36.

Stipp M S H and Heilbronner R.2002.The eastern Tonale fault zone:a"natural laboratory"for crystal plastic deformation of quartz over a temperature range from 250 to 700℃.Journal of Sructural Geology，24:1861－1884.

Twiss R J.1977.Theory and applicability of a recrystallized grain size paleopiezometer.Pure and Applied Geophysics，115:337－340.

Zhang G L，Boulter C A and Liang J C.2001.Brittle origins for disseminated gold mineralization in mylonite:Gaocun gold deposit，Hetai goldfield，Guangdong Province，South China.Economic Geology，96:49－59.

Zhang K J and Cai J X.2009.NE-trending Hepu-Hetai dextral shear zone in southern China:Penetration of the Yunkai promontory of South China into Indochina.Journal of Structural Geology，31:737－748.

Zhang Y H，Lin G，Roberts P and Ord A.2007.Numerical modelling of deformation and fluid flow in the Shuikoushan district，Hunan Province，South China.Ore Geology Reviews，31:261－278.

Zhao C B，Hobbs B E，Mühlhaus H B and Ord A.2001.Finite element modeling of rock alteration and metamorphic process in hydrothermal systems.Communications in Numerical Methods of Engineering，17:833－843.

Implication of Synthetical Methods in Prospecting for the Hetai Goldfield，Western Guangdong

ZHU Jiangjian1，2，CHEN Guanghao3，ZENG Qiaosong1，GONG Guilun1，2，4，YI Jin1，2，GONG Chaoyang5，HUANG Donglin5and LIN Ge1
(1.CAS Key Laboratory of Mineralogy and Metallogeny，Guangzhou Institute of Geochemistry，Chinese Academy of Sciences，Guangzhou510640，Guangdong，China;2.Graduate University of Chinese Academy of Sciences，Beijing100049，China;3.CAS Key Laboratory of Marginal Sea Geology，South China Sea Institute of Oceanology，Chinese Academy of Sciences，Guangzhou510301，Guangdong，China;4.The11thBranch of Chinese Armed Police Gold Force，Ningxiang410600，Hunan，China;5.Hetai Gold Mine，Zhaoqing526127，Guangdong，China)

The Hetai gold deposit is the largest gold mine in Western Guangdong and Southeastern Guangxi，which is related to the ductile shear zone.The anomalous belts of K，geo-electrochemical value and volumetric strain increment were used for prospecting in order to find more gold orebodies in the Hetai goldfield.The anomalous belts of K indicate the K alteration took place during mylonitization and hydrothermal alteration.Auriferous veins and mylonitized zones commonly occur within or nearby the anomalous high value belts of K.The anomalous belts of geo-electrochemical value represent the vertical projection of allopelagic orebodies.The anomalous belts of volumetric strain increment and adjacent areas are possible auriferous mylonitized zone.The anomalous belts of K，geoelectrochemical value and volumetric strain increment correspond well with 9#mylonitized zone.These anomalous belts have been partially or locally proved to be closely related to gold mineralization.These demonstrate that the methods are feasible for prospecting of the deposit.By the contrast researches of the anomalous belt of K，geo-electrochemical value and volumetric strain increment，three possible auriferous anomalous belts are outlined，which might be plausible Au-bearing mineralized mylonitized zones.

Hetai goldfield;mylonitized zone;numerical modeling;Gamma spectrometry;geoelectrical measurement;volumetric strain increment;prospecting

P542;P611.5

A

1001-1552(2012)02-0240-008

2011－08－11;改回日期:2011－10－31

項目資助:國家自然科學(xué)基金項目(41172082)、國家基礎(chǔ)研究計劃(973)項目(2007CB411402)及中國科學(xué)院廣州地球化學(xué)研究所創(chuàng)新項目(A15-030513)聯(lián)合資助。

朱江建(1985－)，男，博士生，構(gòu)造與找礦預(yù)測研究方向。Email:zjj3xsjl@163.com

陳廣浩，男，研究員。Email:chengh@gzb.ac.cn