地球化學(xué)定量預(yù)測法在武夷成礦帶銅資源潛力評價(jià)中的應(yīng)用＊

湛 龍，黃正清，龔 鵬

（1 中國地質(zhì)調(diào)查局南京地質(zhì)調(diào)查中心，南京 210016）（2 中國地質(zhì)大學(xué)信息工程學(xué)院，武漢 430074）

武夷成礦帶位于浙閩贛粵四省邊界的仙霞嶺－武夷山脈地區(qū)，北起浙江衢縣，南至廣東梅縣，南北長約600km，東西寬約160km。區(qū)內(nèi)已全部完成了1：20萬水系沉積物測量和1∶20 萬區(qū)域地質(zhì)及礦產(chǎn)調(diào)查，部分地區(qū)已完成1∶5萬區(qū)域礦產(chǎn)地質(zhì)調(diào)查。區(qū)內(nèi)已發(fā)現(xiàn)大型礦床十余處，中小型礦床數(shù)百處。為了促進(jìn)武夷山成礦帶地質(zhì)找礦工作、闡明區(qū)內(nèi)銅礦產(chǎn)資源潛力，本文采用以地球化學(xué)定量預(yù)測理論為基礎(chǔ)的技術(shù)方法，對銅礦資源的找礦潛力進(jìn)行預(yù)測和探討。

1 區(qū)域地質(zhì)特征

武夷成礦帶大地構(gòu)造位置上隸屬于元古代華夏陸塊活動(dòng)大陸邊緣的武夷－云開島弧系，包括該島弧系北段的火山弧和弧后盆地［1］。東以溫州－南靖－陸豐超巖石圈斷裂與南海－印支陸塊的中生代巖漿活動(dòng)帶為界，北以江紹斷裂與揚(yáng)子陸塊的龍門山－九嶺元古代島弧帶為鄰，西邊是華夏陸塊的加里東褶皺帶，向南與元古代云開島弧帶相接。

圖1 武夷成礦帶大地構(gòu)造位置圖（據(jù)華東地區(qū)大地構(gòu)造圖修編）Fig．1 Geotectonic location of Wuyi metallogenic belt

區(qū)內(nèi)地層自上太古界至第四系發(fā)育齊全。根據(jù)地層巖性、巖相、成巖環(huán)境差異及建造演化史、變質(zhì)程度和構(gòu)造變形特征，將該區(qū)地層分為三個(gè)大的斷代巖系：前泥盆紀(jì)基底巖系、泥盆紀(jì)－中三疊世以海相沉積為主的蓋層巖系和中新生代陸相碎屑及火山巖系，各斷代地層均以明顯的角度不整合為界。

區(qū)內(nèi)構(gòu)造－巖漿巖帶基本呈NE－NNE 向展布，且?guī)r漿活動(dòng)有從NW－SE由老到新的分布特征?？傮w上，除零星出露的元古代花崗巖表征前寒武系結(jié)晶基底的存在外，古生代花崗巖和早中生代花崗巖呈面式分布于區(qū)內(nèi)，而晚侏羅世火山巖漿作用拉開了華南晚中生代大規(guī)模巖漿作用的序幕，形成大面積火山巖和侵入巖，與區(qū)內(nèi)很多成礦活動(dòng)相關(guān)，形成了一系列斑巖型、火山熱液型、隱爆角礫巖型Cu、Pb、Zn、Au、Ag、W、Sn礦床，是我國東部最重要的構(gòu)造－巖漿成礦帶之一。

本區(qū)內(nèi)已知的銅鉛鋅金銀鎢錫鉬等礦種的礦床（點(diǎn)）共一千余個(gè)［1］，如紫金山銅金礦、永平銅礦、冷水坑銀鉛鋅礦、鐵砂街銅礦、行洛坑鎢礦、南平西坑鈮鉭礦、梅縣嵩溪銀銻礦、玉水銅多金屬礦等。整體而言，北部以銅礦為主，也有錫、稀有和金礦出現(xiàn)；南部以鎢錫和稀土、稀有為主。銅礦和鉛鋅礦與同熔型花崗巖有關(guān)，鎢、錫礦則與重熔型花崗巖有關(guān)。

2 定量預(yù)測思路與方法

目前地球化學(xué)塊體法、體積品味法、豐度法、面金屬量法是勘查地球化學(xué)界中應(yīng)用最廣泛的四種定量預(yù)測方法［2］。不同的預(yù)測方法適應(yīng)于不同的預(yù)測尺度。例如，地球化學(xué)塊體法［3］是對大型或巨型礦床而提出的，多適應(yīng)于成礦省或成礦域的資源量預(yù)測；體積品味法更適用于礦田或礦床的資源量預(yù)測。自全國礦產(chǎn)資源潛力評價(jià)項(xiàng)目開展以來，以中國地質(zhì)大學(xué)（武漢）馬振東教授為首的科研團(tuán)隊(duì)經(jīng)過反復(fù)試點(diǎn)研究，在豐度法和面金屬量法的基礎(chǔ)上，完成了礦產(chǎn)資源地球化學(xué)模型建立與定量預(yù)測方法技術(shù)研究??。本文應(yīng)用該技術(shù)方法，對武夷成礦帶銅礦進(jìn)行地球化學(xué)資源量定量預(yù)測研究，取得了良好的應(yīng)用效果。

地球化學(xué)定量預(yù)測包括兩個(gè)核心內(nèi)容，即典型礦床地球化學(xué)模型的建立和定量預(yù)測研究。地球化學(xué)模型實(shí)質(zhì)上由成礦地質(zhì)特征和地球化學(xué)特征兩個(gè)方面組成；而定量預(yù)測研究是在建模的基礎(chǔ)上，根據(jù)預(yù)測區(qū)的礦床類型，選擇相應(yīng)的地球化學(xué)定量預(yù)測計(jì)算方法進(jìn)行資源量估算。具體工作流程包括：①搜集資料：搜集工作區(qū)內(nèi)各尺度勘查地球化學(xué)數(shù)據(jù)、地質(zhì)礦產(chǎn)圖、主要銅礦床勘查儲(chǔ)量報(bào)告等；②建立預(yù)測模型：建立成礦物源、時(shí)空分布規(guī)律、成礦模式、主控礦因素、區(qū)域地球化學(xué)異常特征與礦區(qū)地球化學(xué)特征；③相關(guān)數(shù)據(jù)處理及地球化學(xué)圖件的編制，包括相似度圖、襯值異常圖、剝蝕程度圖等；④圈定預(yù)測靶區(qū)：依據(jù)圈定要素（相似度、剝蝕程度），根據(jù)預(yù)測區(qū)評價(jià)準(zhǔn)則，將遴選出的靶區(qū)進(jìn)行評價(jià)（分成A、B、C類異常）；⑤預(yù)測資源量：依據(jù)礦床類型、預(yù)測尺度，選擇地球化學(xué)預(yù)測方法進(jìn)行儲(chǔ)量計(jì)算。

2．1 資料來源

本次用于武夷成礦帶地球化學(xué)定量預(yù)測的基礎(chǔ)資料為成礦帶內(nèi)涉及江西、浙江、福建、廣東4省范圍內(nèi)相關(guān)1∶20萬水系沉積物測量數(shù)據(jù)，面積約16萬km2，樣品件數(shù)為40120件，分析元素為39種，獲得分析數(shù)據(jù)約為1．56×106個(gè)。依據(jù)中國地質(zhì)調(diào)查局《關(guān)于印發(fā)（區(qū)域化探資料質(zhì)量評估要求）的函》與《1∶20萬區(qū)域化探資料質(zhì)量評估方法》，各省數(shù)據(jù)均通過了專家的評審，數(shù)據(jù)質(zhì)量可靠。

2．2 典型礦床模型的建立

地球化學(xué)資源量預(yù)測的基本思想是相似類比，典型礦床的地質(zhì)、地球化學(xué)模型的建立是進(jìn)行成礦預(yù)測的基礎(chǔ)。武夷成礦帶銅多金屬成礦以巖漿成因?yàn)橹?，與海相火山活動(dòng)、陸相火山活動(dòng)、巖漿侵入活動(dòng)（如斑巖－潛火山活動(dòng)）有關(guān)，尤其與燕山晚期中酸性斑巖及其熱液活動(dòng)關(guān)系密切。本次定量預(yù)測選取7個(gè)銅礦典型礦床，幾乎涵蓋了武夷成礦帶所有的銅礦成礦類型，具體特征見表1。因每個(gè)地球化學(xué)模型的建模思路與方法基本一致，在此以紫金山銅金礦為例，說明相似元素組合、剝蝕程度指標(biāo)的挑選過程，其地質(zhì)特征等參見文獻(xiàn)［1］，本文主要介紹典型礦床模型的地球化學(xué)特征。

紫金山銅金礦區(qū)域地球化學(xué)特征：1∶20 萬化探異常元素組合Cu－Mo－Au－Pb－Zn－W－Ag－Bi－Sn－Cd－As－Sb，異常面積152km2，區(qū)域上該異常呈島狀展布，多元素套合程度高，濃集中心明顯，形成非常醒目的地球化學(xué)異常（圖2），異常與紫金山復(fù)式巖體范圍基本一致，受巖體和北東向構(gòu)造控制。

1∶5萬區(qū)域地球化學(xué)異常（其化探異常特征見表2）與1∶20萬區(qū)域化探異常有繼承性關(guān)系，同樣形成大規(guī)模分布的Cu－Pb－Mo－Au－Ag－Zn－W－As組合異常，其主要成礦元素異常規(guī)模大，元素組合復(fù)雜，元素套合程度高。

圖2 上杭紫金山銅金礦區(qū)1∶20萬區(qū)域化探異常剖析圖Fig．2 1：200000regional geochemical anomaly map of Zijinshan Copper－gold orefield

表1 武夷成礦帶銅典型礦床特征Table 1 Characteristics of typical copper deposits in Wuyi metallogenic belt

區(qū)內(nèi)局部開展過50km21∶1萬原生暈地球化學(xué)測量工作，根據(jù)地表原生暈測量元素定量分析成果，主要異常分布于紫金山復(fù)式巖體及其接觸帶上（圖3）。其原生暈地球化學(xué)測量工作可以看出，其原生暈分帶自上而下為：Au、Ag、Bi、Hg、Mo（As）（前緣暈）→Cu、As、Sb、Pb、Bi、Au、Ag（礦中暈）→Cu、W、Sn、Mo（礦尾暈）。結(jié)合其1∶20萬區(qū)域地球化學(xué)特征、1∶5萬地球化學(xué)特征、礦區(qū)1：1萬地球化學(xué)特征及原生暈地球化學(xué)特征，挑選Cu－Pb－Mo－Bi－Au－Ag作為其相似元素組合，并選擇一組礦尾暈與前緣暈元素組合的累加之比（W＋Sn）／（Ag＋Hg）來評價(jià)該類型礦床的剝蝕程度。

2．3 相關(guān)地球化學(xué)圖件的編制

（1）編制相似度圖

標(biāo)準(zhǔn)樣本的指定：相似度是指地球化學(xué)相似系數(shù)或相近系數(shù)。從元素組合的角度，通過構(gòu)建“標(biāo)準(zhǔn)樣本”元素組合，利用距離公式定量判斷未知成礦元素組合與已知區(qū)成礦元素組合之間的相似程度或相近程度的參數(shù)，它是判別未知礦化與已知礦化關(guān)系的重要參數(shù)。

表2 上杭紫金山1∶5萬化探異常特征［4］Table 2 Characteristics of 1∶50000geochemical anomalies of Zijinshan

圖3 上杭紫金山銅金礦區(qū)1∶1萬化探異常剖析圖［4］Fig．3 1∶10000geochemical anomaly map of Zijinshan copper－gold orefield

制作“標(biāo)準(zhǔn)樣本”對選取的元素組合進(jìn)行賦值，為了取值比較合理，每個(gè)典型礦床采用其各元素異常內(nèi)帶的平均含量進(jìn)行賦值，這樣可以克服極個(gè)別特高值點(diǎn)對相似度計(jì)算的影響。本次定量預(yù)測工作選取的7個(gè)典型礦床的標(biāo)準(zhǔn)樣本參數(shù)見表3。利用制作的典型礦床標(biāo)準(zhǔn)樣本，計(jì)算每個(gè)原始1∶20萬水系沉積物數(shù)據(jù)的相似度值：先計(jì)算樣點(diǎn)與典型礦床的相似距離（D（Si）），再通過公式轉(zhuǎn)為相似系數(shù)（R）。

圖4 紫金山銅金礦床地球化學(xué)異常模型?Fig．4 Geochemical anomaly model of Zijinshan coppergold orefield

表3 武夷成礦帶典型礦床標(biāo)準(zhǔn)樣本參數(shù)Table 3 Standard sample parameter table of typical deposits in Wuyi metallogenic belt

其中，s代表已知礦床（標(biāo)準(zhǔn)樣本），i為需要判別的未知樣點(diǎn)，p 為選取的典型礦床的變量個(gè)數(shù)，Xsk表示標(biāo)準(zhǔn)樣本在第k個(gè)變量（如Cu元素）上的取值，Xik表示第i個(gè)樣本在第k個(gè)變量上的取值，Max（D）為樣點(diǎn)間距離的最大值。每個(gè)樣點(diǎn)計(jì)算出相似度值后，采用7 級(jí)累頻分級(jí)的方法，分級(jí)頻數(shù)為：25%－50%－75%－90%－95%－98%，各典型礦床分級(jí)標(biāo)準(zhǔn)見表4。對應(yīng)于相似度的累頻分級(jí)，可以認(rèn)為98%的含義是2%的樣點(diǎn)與已知礦床相似，屬于成礦信息，相似度的累頻分級(jí)對應(yīng)一定的成礦率。

表4 武夷成礦帶典型礦床相似度圖累頻分級(jí)Table 4 Tired frequency grading of typical deposits in Wuyi metallogenic belt

（2）編制襯值異常圖

襯值處理的基本方法為，以計(jì)算樣點(diǎn)為中心，按給定背景范圍搜索數(shù)據(jù)并計(jì)算背景平均值，計(jì)算點(diǎn)值與所處背景平均值之比值，其中在某一尺寸對高分位值有一個(gè)“拐點(diǎn)”效應(yīng)，即隨著尺寸的繼續(xù)加大，對高異常的突出有限；當(dāng)若窗口開的過大，異常形態(tài)則過于集中，無法區(qū)分不同礦床各自的異常，襯值圖效果過猶不及。經(jīng)試驗(yàn)后，本次圈定襯值異常最終選擇30km 作為背景尺寸。經(jīng)過襯值處理過的數(shù)據(jù)值，再采用累頻分級(jí)的方式成圖，具體分級(jí)的標(biāo)準(zhǔn)與相似度圖一致。

（3）編制剝蝕程度圖

按照原生暈分帶理論，由于元素的活動(dòng)性和沉淀溫度的差異，往往會(huì)在礦體周圍形成低溫遠(yuǎn)礦元素組合（礦頭暈）、中溫近礦元素組合（礦中暈）、高溫礦尾元素組合（礦尾暈）。在不考慮多期疊加的情況下，哪種暈對應(yīng)元素組合出現(xiàn)的強(qiáng)度越高、規(guī)模越大，就說明礦床剝蝕到了哪個(gè)程度。應(yīng)用水系沉積物的元素組合特征可以判別礦床的剝蝕程度，不同的礦床類型具備不同的分帶組合，以礦帶內(nèi)典型礦床的巖石原生暈測量為基礎(chǔ)，結(jié)合礦床的蝕變和礦化分帶，選擇與實(shí)際情況吻合最佳的元素組合（實(shí)驗(yàn)對比），如紫金山陸相火山巖型銅金礦床選擇的元素組合為：（W＋Sn）／（Ag＋Hg），浦城管查熱液脈狀型銅礦床選擇的元素組合為：（Sn＋As）／（Pb＋Zn），會(huì)昌紅山隱爆角礫巖－斑巖型銅礦床和平和鐘騰斑巖型銅礦床選擇的元素組合為（W＋Sn＋Mo）／（As＋Sb＋Hg）；對于梅州玉水海底熱水噴流－沉積疊加改造型銅礦床、弋陽鐵砂街海底熱水噴流疊加改造型銅礦床、鉛山永平層控疊加改造型銅礦床，因巖石原生暈資料缺失，無法確定判別礦頭暈與礦尾暈的元素組合，暫不考慮剝蝕程度。比值計(jì)算完后，按7級(jí)累頻分級(jí)的方法進(jìn)行成圖，分級(jí)頻數(shù)為50%－85%－90%－95．5%－98%－99%。

2．4 預(yù)測指標(biāo)選取與分級(jí)準(zhǔn)則

2．4．1 分級(jí)指標(biāo)

在典型礦床地質(zhì)、地球化學(xué)特征研究及地球化學(xué)找礦模型建立的基礎(chǔ)上，對圈定預(yù)測區(qū)所需的地質(zhì)－地球化學(xué)要素進(jìn)行提取，概括為以下7個(gè)指標(biāo)：

條件1：元素組合與典型銅多金屬礦床的相似度值大于累頻分級(jí)≥98%的對應(yīng)值（如紫金山的對應(yīng)值為0．59）；

條件2：成礦地質(zhì)條件有利（巖體、巖體與圍巖接觸帶、斷層等）；

條件3：已發(fā)現(xiàn)礦點(diǎn)（礦化點(diǎn)）；

條件4：相似典型礦床所對應(yīng)的特征元素組合平均襯值較大（累頻分級(jí)≥95%，如紫金山銅礦對應(yīng)的特征元素組合為Cu－Pb－Au－Ag－Mo－Bi）；

條件5：Cu 襯值異常具備中帶；

圖5 A、B級(jí)預(yù)測區(qū)圈定過程示意圖Fig．5 Schematic diagrams of delineation of A and B grade prediction zones

條件6：相似典型礦床所對應(yīng)的特征元素組合中至少三個(gè)元素襯值不低于累頻分級(jí)95%對應(yīng)值；

條件7：相似典型礦床所對應(yīng)的特征元素組合中至少四個(gè)元素襯值不低于累頻分級(jí)95%對應(yīng)值。

2．4．2 分級(jí)準(zhǔn)則

由以上7個(gè)地質(zhì)地球化學(xué)指標(biāo)，根據(jù)有無相似度的準(zhǔn)則劃分為A、B、C 三級(jí)，具體分級(jí)標(biāo)準(zhǔn)如下：

A級(jí)：滿足條件1－2－3－4－5－6；

B級(jí)：滿足條件1－2－4－5－6（缺少3）；

C級(jí)：滿足條件2－4－5－7（缺少1 和3，必須滿足7）。

其中，可信度級(jí)別為：A＞B＞C。

2．5 各級(jí)預(yù)測區(qū)的圈定

首先將相似度圖、各相關(guān)元素襯值異常圖、礦產(chǎn)圖進(jìn)行疊合，將與典型礦床相似度高，且具備3種以上特征元素襯值異常、異常區(qū)有已知銅礦（化）點(diǎn)的區(qū)域先圈為A 類預(yù)測區(qū)（如圖5中A－YC－22，相似度高、有Cu－Zn－Ag－Cd 襯值異常，且有已知銅礦點(diǎn)）；將僅有相似度高和襯值異常、而無已知礦（化）點(diǎn)的區(qū)域圈為B類預(yù)測區(qū)（如圖5中B－YC－28僅具備相似度高和Cu－Pb－Zn－Cd 異常，無已知銅礦化點(diǎn)）；將相似度較低，但具備4種以上襯值異常的圈為C類預(yù)測區(qū)，再結(jié)合地質(zhì)礦產(chǎn)圖，驗(yàn)證該區(qū)域是否具備典型礦床所需的必要控礦要素。利用此方法，在武夷成礦帶共圈定A 級(jí)預(yù)測靶區(qū)23處，B 級(jí)預(yù)測靶區(qū)31處，C級(jí)預(yù)測靶區(qū)37處。

3 銅資源量預(yù)測結(jié)果

選用豐度法和面金屬量法兩種方法對A、B 級(jí)預(yù)測靶區(qū)進(jìn)行資源量估算，C 級(jí)靶區(qū)由于相似度指標(biāo)達(dá)不到要求，缺乏最佳相似礦床數(shù)據(jù)，不能計(jì)算資源量，只能作為有利的成礦遠(yuǎn)景區(qū)。上述兩種方法均考慮剝蝕系數(shù)，類比法計(jì)算公式為：預(yù)測區(qū)資源量（Vd）＝［預(yù)測區(qū)異常規(guī)模Q×（1－預(yù)測區(qū)剝蝕系數(shù)F）×相似系數(shù)R］／［系數(shù)K× （1－典型礦床剝蝕系數(shù)F）］；面金屬量法計(jì)算公式為：預(yù)測區(qū)資源量（Vs）＝［預(yù)測區(qū)面金屬量P×（1－預(yù)測區(qū)剝蝕系數(shù)F）×相似系數(shù)R］／［系數(shù)K×（1－典型礦床剝蝕系數(shù)F）］；最終取上述兩類加權(quán)平均的方法（類比法加權(quán)系數(shù)為0．6，面金屬量法加權(quán)系數(shù)為0．4來確定資源量：預(yù)測資源量（V）＝0．6×類比法（Vd）＋0．4×面金屬量法（Vs）。

據(jù)全國礦產(chǎn)資源潛力評價(jià)項(xiàng)目最新礦產(chǎn)地?cái)?shù)據(jù)資料顯示，武夷成礦帶已探明銅資源量約573．83萬t（主要包括上杭紫金山196．72萬t、管查2．4萬t、紅山49．95萬t、鐵砂街5．61萬t、永平146．68萬t、玉水11萬t、鐘騰2．24萬t、蘿卜嶺8．13萬t），經(jīng)本次預(yù)測，新增A 級(jí)預(yù)測區(qū)預(yù)測資源量452．3萬t，B級(jí)預(yù)測區(qū)預(yù)測資源量113．4萬t。根據(jù)礦床式分類，在預(yù)測的資源量中，以紅山式斑巖型銅礦227萬t（預(yù)測區(qū)7個(gè)），永平層控疊加改造型銅礦137．7萬t（預(yù)測區(qū)9 個(gè)）和紫金山式的陸相火山巖型銅礦136．1萬t（預(yù)測區(qū)5個(gè)）為主。預(yù)測結(jié)果顯示武夷山成礦帶具有較強(qiáng)的銅礦找礦潛力，為武夷山成礦帶礦床勘察工作部署提供了有利的地球化學(xué)依據(jù)。

4 預(yù)測成果的檢驗(yàn)

在本次圈定的A－YC－13預(yù)測區(qū)內(nèi)，其樣本與紅山銅礦標(biāo)準(zhǔn)樣本顯示出極高的相似度，相似度值達(dá)0．84，銅異常三級(jí)濃度分帶明顯，又具備Cu－Pb－Zn－Ag等多元素套合襯值異常和燕山期中酸性侵入體的存在，因此推測在該區(qū)域具備尋找火山熱液型－斑巖型銅礦的前景。南京地質(zhì)調(diào)查中心在該區(qū)域開展的浙閩松政地區(qū)礦產(chǎn)遠(yuǎn)景調(diào)查工作，已經(jīng)相繼在該區(qū)域發(fā)現(xiàn)了多條隱爆角礫巖脈和青盤巖化、黃鐵絹英巖化等蝕變現(xiàn)象，圈出了3條細(xì)脈浸染狀銅礦化帶。目前已施工三個(gè)鉆孔，發(fā)現(xiàn)多層銅礦體。

5 結(jié)論

（1）在前人研究的基礎(chǔ)上，地球化學(xué)定量預(yù)測首次把剝蝕系數(shù)和相似系數(shù)引入至地球化學(xué)資源量估算公式中，使地球化學(xué)預(yù)測方法更為合理，其研究體系包括基礎(chǔ)地質(zhì)－成礦機(jī)制－理論地球化學(xué)－勘查地球化學(xué)－現(xiàn)代計(jì)算機(jī)技術(shù)多個(gè)方面，能夠整合所有資料為礦床勘察工作部署提供有利的地球化學(xué)依據(jù)。

（2）資源量的地球化學(xué)估算中，典型礦床模型的建立是根本，其資料的完善程度，資源量的可靠程度直接影響定量預(yù)測的結(jié)果。

（3）本次研究在武夷成礦帶圈出了23處A 級(jí)預(yù)測區(qū)、31處B級(jí)預(yù)測區(qū)，A、B 級(jí)預(yù)測區(qū)的銅預(yù)測資源量分別為452．3萬t、113．4萬t。顯示出該成礦帶具有較強(qiáng)的銅礦找礦潛力，為武夷成礦帶的找礦部署提供了有利的地球化學(xué)依據(jù)。

注釋

?馬振東．礦產(chǎn)資源地球化學(xué)模型建立與定量預(yù)測方法技術(shù)一（內(nèi)部培訓(xùn)材料）．2010．

?馬振東．礦產(chǎn)資源地球化學(xué)模型建立與定量預(yù)測方法技術(shù)二（內(nèi)部培訓(xùn)材料）．2010．

?姚敬金．中國主要大型有色、貴金屬礦床綜合信息找礦模型（內(nèi)部資料）．2002．

［1］陳世忠，黃正清，朱筱婷，等．武夷山成礦帶銅多金屬礦勘查選區(qū)評價(jià)項(xiàng)目報(bào)告［R］．南京：南京地質(zhì)調(diào)查中心，2011：97－167．

［2］龔鵬，李娟，胡小梅，等．區(qū)域地球化學(xué)定量預(yù)測方法技術(shù)在礦產(chǎn)資源潛力評價(jià)中的應(yīng)用［J］．地質(zhì)論評，2004，58（6）：1101－1107．

［3］謝學(xué)錦，劉大文，向運(yùn)川，等．地球化學(xué)塊體概念和方法學(xué)的發(fā)展［J］．中國地質(zhì)，2002，29（3）：225－233．

［4］夏春金，黃茜，朱律運(yùn)．福建省礦產(chǎn)資源潛力評價(jià)化探資料應(yīng)用成果報(bào)告［R］．福建：福建省地質(zhì)調(diào)查研究院，2013：57－162．