綜合找礦方法在河臺金礦找礦預測中的應用

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

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

綜合找礦方法在河臺金礦找礦預測中的應用

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

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

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

河臺金礦;數值模擬;伽馬能譜;地電化學;體積應變;找礦預測

0 引 言

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

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

1 礦床地質特征

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

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

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

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

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

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

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

3 各種異常帶在本區的指示意義

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

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

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

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

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

3.2 地電化學異常帶

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

3.3 體積應變異常帶

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

本文首次運用模擬產生的體積應變異常帶對本區進行找礦預測研究。鑒于本次數值模擬所遵從的地質依據(朱江建等，2011a)，認為模擬結果在理論上是可靠的。考慮到本區各巖性的實測參數不能完全代表糜棱巖化帶形成時的參數及糜棱巖形成過程中相關的參數是變化的，而本次模擬的整個過程中相關參數保持不變等因素，認為體積應變異常帶及其附近有利于發展為本區的含礦糜棱巖化帶(朱江建等，2011a)。

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

4 找礦預測

4.1 找礦預測可行性分析

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

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

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

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

4.2 找礦預測可靠性分析

由K、地電化學、體積應變異常帶的對比研究，同時結合本區的地質實際，在本區圈定了3個可能的含礦異常帶。

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

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

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

綜上所述，圈定的3個可能的含礦異常帶都有地質依據，K、地電化學、體積應變異常帶與這3個異常帶都有較好的對應關系，建議對其進行驗證。

5 結 論

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

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

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

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

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

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

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

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

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

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

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

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

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

林舸，趙重斌，張晏華，王岳軍，單業華.2005.地質構造變形數值模擬研究的原理、方法及相關進展.地球科學進展，20(5):549－555.

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

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

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

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

羅先熔，王衛民，張佩華.1997.隱伏金礦地電化學異常形成機制及異常形態特征.有色金屬礦產與勘查，6(6):364－367.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

朱江建，曾喬松，易金，龔朝陽，李新福，黃棟林，王斯亮，夏永健，龔貴倫，陳廣浩，林舸.2011b.地電化學測量在河臺金礦找礦預測中的應用.地質與勘探，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

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

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

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