湖南花垣地區(qū)鉛鋅礦床稀土元素地球化學(xué)特征初步研究

周 云，段其發(fā)，曹 亮，彭三國，甘金木

ZHOU Yun,DUAN Qi-Fa,CAO Liang,PENG San-Guo,GAN Jin-Mu

（中國地質(zhì)調(diào)查局武漢地質(zhì)調(diào)查中心，武漢430205）

(Wuhan Center of China Geological Survey,Wuhan 430205,China)

湖南花垣礦集區(qū)位于揚(yáng)子地塊東南邊緣與雪峰（江南）造山帶的過渡部位（圖1），構(gòu)造上位于NE向湘黔斷裂帶中部，經(jīng)歷了晉寧期碰撞造山、加里東期拉張裂陷－被動(dòng)大陸邊緣盆地－前陸盆地以及印支期－燕山期陸內(nèi)變形等演化過程，形成了多種沉積/層控金屬礦產(chǎn)和非金屬礦產(chǎn)，是湘西－鄂西成礦帶的重要組成部分[1]。

圖1 花垣礦集區(qū)地質(zhì)礦產(chǎn)簡(jiǎn)圖(據(jù)文獻(xiàn)[1]修編)Fig.1 Geological and mineral map of Huayuan ore area(revised according to ref.[1])

花垣礦田是一個(gè)鉛鋅資源儲(chǔ)量達(dá)千萬噸級(jí)的超大型鉛鋅礦田，也是整個(gè)湘西地區(qū)鉛鋅礦床的典型代表。自取得找礦突破以來，前人對(duì)其成礦地質(zhì)背景、礦床地質(zhì)特征、巖相古地理特征、地層古生物、稀土元素地球化學(xué)特征、C、H、O、S、Pb、Sr等穩(wěn)定同位素、成礦流體、成礦地質(zhì)時(shí)代、成礦物質(zhì)來源及成礦模式等方面進(jìn)行了大量較為深入的研究[2-24]，但該區(qū)稀土元素的相關(guān)研究報(bào)道僅限于花垣獅子山礦床的成礦期方解石和碳酸鹽圍巖[25]，未對(duì)其礦石礦物進(jìn)行該方面的研究。礦石礦物稀土元素特征取決于礦物沉淀時(shí)成礦流體稀土元素的組成和沉淀時(shí)的物理化學(xué)條件。因此，根據(jù)礦石礦物的稀土元素地球化學(xué)研究，可以探討成礦物質(zhì)來源或者成礦流體來源[26]。

本文分析了花垣地區(qū)鉛鋅礦床閃鋅礦、方鉛礦、黃鐵礦三種礦石礦物與圍巖灰?guī)r的稀土元素地球化學(xué)特征，探討了賦礦圍巖與成礦的關(guān)系、成礦流體來源及成礦物質(zhì)來源，從而為整個(gè)湘西地區(qū)鉛鋅礦床成因的研究提供參考。

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

花垣礦集區(qū)總體走向北北東，位于花垣-張家界深大斷裂與麻栗場(chǎng)斷裂之間（圖1），長(zhǎng)38 km，寬4～16 km，面積約215 km2[21]。礦集區(qū)內(nèi)主要出露了從青白口系至寒武系的一系列地層。下寒武統(tǒng)清虛洞組為區(qū)內(nèi)主要的容礦地層，在礦區(qū)內(nèi)分布廣泛。礦區(qū)內(nèi)巖漿巖不發(fā)育，未見巖體出露。礦集區(qū)內(nèi)斷裂構(gòu)造和褶皺構(gòu)造較發(fā)育，斷裂主要有花垣-保靖-張家界深大斷裂、兩河-長(zhǎng)樂斷裂和麻栗場(chǎng)斷裂（保靖-銅仁-玉屏深大斷裂北東段）（圖1）。褶皺構(gòu)造主要為寬緩的背（向）斜構(gòu)造。

鉛鋅礦體主要賦存于下寒武統(tǒng)清虛洞組下段灰?guī)r中，以層狀、似層狀礦體為主，次為脈狀。礦石礦物簡(jiǎn)單，以閃鋅礦和方鉛礦為主，次為黃鐵礦，脈石礦物主要為方解石（圖2），次為重晶石，少量的石英和螢石。閃鋅礦主要呈黃棕色、淺黃色、黃綠色，片狀晶形，具半金屬光澤。結(jié)晶顆粒較粗大，具半自形或他形晶粒狀結(jié)構(gòu)，粒徑一般為0.5～2 mm，最小粒徑為0.05 mm，最大粒徑可達(dá)4～6 mm，以粒狀或脈狀集合體形式分布于礦石中，脈狀閃鋅礦多分布于與灰?guī)r接觸的方解石脈體邊緣。方鉛礦呈鉛灰色，顆粒大小為0.01～0.2 mm，具半自形晶～他形晶顆粒結(jié)構(gòu)，在脈石礦物和閃鋅礦的裂隙中主要呈不規(guī)則粒狀及細(xì)脈狀分布，有時(shí)可呈粗大團(tuán)塊狀（圖2）。圍巖蝕變類型主要為方解石化，還常伴隨有弱的黃鐵礦化、螢石化、重晶石化等[1,27]，其中方解石化與礦化關(guān)系密切。

2 樣品采集與測(cè)試方法

本文分析了花垣鉛鋅礦田22件原生礦石中的閃鋅礦、方鉛礦、黃鐵礦硫化物和圍巖灰?guī)r的稀土元素含量，包括7件閃鋅礦單礦物樣品，6件方鉛礦單礦物樣品，2件黃鐵礦單礦物樣品，以及7件賦礦圍巖下寒武統(tǒng)清虛洞組下段灰?guī)r樣品。采集的閃鋅礦、黃鐵礦、方鉛礦樣品分別來自花垣由北而南依次分布的團(tuán)結(jié)、李梅、土地坪、蜂塘、大石溝鉛鋅礦床（圖1），閃鋅礦和方鉛礦、黃鐵礦多沿脈石礦物方解石邊緣分布，或呈團(tuán)塊狀、斑狀、浸染粒狀分布于方解石脈體中或其邊緣（圖2）。

首先將樣品粉碎至40目，在雙目鏡下挑純閃鋅礦、方鉛礦和黃鐵礦的單礦物，確保單礦物純度高于98%，用去離子水清洗挑純的硫化物單礦物，低溫干燥后，用瑪瑙研缽研磨至200目置于容器中備用。稀土元素含量分析由武漢地質(zhì)調(diào)查中心巖礦測(cè)試實(shí)驗(yàn)室完成，采用電感耦合等離子體質(zhì)譜ICP-MS方法測(cè)試，分析精度一般優(yōu)于5%，檢測(cè)下限為（0.n～n）×10-9。

3 稀土元素分析結(jié)果及地球化學(xué)特征

樣品REE元素含量分析結(jié)果及相關(guān)參數(shù)計(jì)算見表1，圖3為相應(yīng)的球粒隕石標(biāo)準(zhǔn)化配分模式圖。由表1和圖3可見，7件閃鋅礦的稀土元素總量 ΣREE變化范圍為 0.16×10-6～0.33×10-6，平均值為0.23×10-6；ΣLREE變化范圍為0.12×10-6～0.28×10-6，ΣHREE 變化范圍為 0.04×10-6～0.05×10-6，輕重稀土元素比值LREE/HREE變化范圍為 2.96 ～ 5.73，相對(duì)富集輕稀土；（La/Yb）N＝0.23～3.98，平均為1.76＞1，富集輕稀土；δEu值為1.65～2.80，平均為2.30，存在中等程度的正Eu異常，δCe值變化范圍為1.16～2.41，平均為1.56，具有弱的正Ce異常。以球粒隕石為標(biāo)準(zhǔn)的稀土元素分布模式圖顯示稀土元素模式曲線為右傾型。

6件方鉛礦的ΣREE變化范圍為0.88×10-6～6.34×10-6，平均值為 2.26×10-6；ΣLREE 變化范圍為 0.82×10-6～6.23×10-6，ΣHREE 變化范圍為0.06×10-6～0.12×10-6，輕重稀土元素比值LREE/HREE變化范圍為13.53～53.77，相對(duì)富集輕稀土；（La/Yb）N＝7.72 ～ 200.52，平均為 48.60＞1，輕重稀土分餾明顯，富集輕稀土；δEu值為0.57～1.09，平均為0.91，存在弱的負(fù)Eu異常，δCe值變化范圍為0.81～3.35，平均為2.16，具有Ce正異常特征。以球粒隕石為標(biāo)準(zhǔn)的稀土元素分布模式圖顯示稀土元素模式曲線為右傾型。

2件黃鐵礦的ΣREE變化范圍為0.46×10-6～0.49×10-6，平均值為 0.47×10-6；ΣLREE 變化范圍為 0.40×10-6～0.43×10-6，ΣHREE 變化范圍為0.0603×10-6～0.0627×10-6，輕重稀土元素比值LREE/HREE變化范圍為6.55～6.82，相對(duì)富集輕稀土；（La/Yb）N＝6.92～7.83，平均為 7.37＞1，輕重稀土分餾明顯，富集輕稀土；δEu值為1.03～1.14，平均為1.08，存在弱的正Eu異常，δCe值變化范圍為1.17～1.25，平均為1.21，具有弱的Ce正異常。以球粒隕石為標(biāo)準(zhǔn)的稀土元素分布模式圖顯示稀土元素模式曲線為右傾型。

表1 湖南花垣地區(qū)鉛鋅礦床單礦物與灰?guī)r稀土元素含量表（ωB/1 0-6）Table 1 REE contents in minerals and limestone of the lead-zinc deposits in Huayuan area, Hunan

表2 湖南花垣地區(qū)鉛鋅礦床單礦物與灰?guī)r稀土元素特征參數(shù)表Table 2 REE statistics parameters of the minerals and limestone in the lead-zinc deposits in Huayuan area, Hunan

圍巖灰?guī)r的 ΣREE 為 3.70×10-6～11.32×10-6，平均值為 5.75×10-6；ΣLREE變化范圍為3.18×10-6～ 8.92×10-6，ΣHREE 變 化 范 圍 為0.51×10-6～ 2.40×10-6，輕重稀土元素比值LREE/HREE變化范圍為3.72～7.95，相對(duì)富集輕稀土；（La/Yb）N＝2.66 ～ 10.43，平均為 7.25＞1，輕重稀土分餾明顯，富集輕稀土；δEu值變化范圍為0.63～1.15，平均為0.86，δCe值變化范圍為0.74～1.00，平均為0.82，存在中等程度的負(fù)Eu異常和負(fù)Ce異常。以球粒隕石為標(biāo)準(zhǔn)的稀土元素分布模式圖顯示稀土元素模式曲線為右傾型。

（La/Sm）N比值反映了輕稀土之間的分餾程度，比值越大，輕稀土越富集[28]。花垣鉛鋅礦床的礦石硫化物（La/Sm）N比值的平均值為4.76＞1.0，賦礦圍巖灰?guī)r的（La/Sm）N比值的平均值為4.41＞1.0，表明輕稀土分餾明顯，礦石硫化物和灰?guī)r均屬于輕稀土富集型。

（Gd/Yb）N比值反映了輕稀土之間的分餾程度，比值越小，重稀土富集程度越高[28]花垣鉛鋅礦床的礦石硫化物（Gd/Yb）N比值的平均值為1.77＞1.0，賦礦圍巖灰?guī)r的（Gd/Yb）N比值的平均值為1.71＞1.0，表明重稀土發(fā)生了一定程度的分餾作用，礦石硫化物和灰?guī)r均屬于重稀土虧損型。

圖3 花垣地區(qū)鉛鋅礦床硫化物與賦礦圍巖灰?guī)r稀土元素球粒隕石標(biāo)準(zhǔn)化配分模式圖Fig.3 Chondrite-normalized REE patterns of sulfides and hostrock in the lead-zinc deposits in Huayuan area

4 討論

4.1 稀土模式成因

總體而言，花垣地區(qū)鉛鋅礦床礦石硫化物表現(xiàn)出較低的總稀土含量特征，閃鋅礦的ΣREE平均值為0.23×10-6，方鉛礦的ΣREE平均值為2.26×10-6，黃鐵礦的ΣREE平均值為0.47×10-6，明顯低于圍巖灰?guī)r的稀土元素總量平均值5.75×10-6。圖4可以看出，礦石硫化物的ΣREE顯著低于圍巖灰?guī)r。整個(gè)礦床中ΣREE含量變化范圍不大，稀土元素的分布較為均勻。熱液礦物由于受礦物晶體結(jié)構(gòu)及物理化學(xué)性質(zhì)等因素的影響，會(huì)導(dǎo)致稀土元素總量的差別主要控制。含Ca的礦物通常具有較高的含量，因?yàn)橄⊥猎乜梢灶愘|(zhì)同象形式替代含鈣質(zhì)礦物[29]。而稀土元素?zé)o法交代閃鋅礦等硫化物中的陽離子進(jìn)入硫化物晶格，只能在其中以包裹體的形式存在[30]。因此，成礦流體中的稀土元素主要分布在方解石等含Ca礦物中，而石英、重晶石、閃鋅礦等礦物則具有較低的稀土元素總量[29]。這可能是導(dǎo)致花垣地區(qū)鉛鋅礦床硫化物礦物中稀土元素含量低的原因。

花垣地區(qū)鉛鋅礦床中閃鋅礦和黃鐵礦的稀土元素球粒隕石標(biāo)準(zhǔn)化配分模式重疊性相對(duì)較好，反映了這兩種硫化物的同源性。方鉛礦與閃鋅礦、黃鐵礦的稀土元素球粒隕石標(biāo)準(zhǔn)化配分模式也較為相似，但是方鉛礦稀土元素總量ΣREE多高于閃鋅礦和黃鐵礦的ΣREE平均值，反映了礦床成礦物質(zhì)可能是多來源的。方鉛礦主要出露于花垣礦集區(qū)的南西端，如蜂塘、大石溝礦床，其形成晚于閃鋅礦和黃鐵礦，多穿插交代閃鋅礦礦物（圖2e，圖2 f），且一般呈脈狀沿構(gòu)造裂隙分布。由于花垣礦集區(qū)成礦流體的運(yùn)動(dòng)方向?yàn)橛杀辈苛飨蚰喜縖16]，因此，方鉛礦應(yīng)是后期多階段和構(gòu)造熱液等多種成礦作用共同疊加的結(jié)果。

圖4 花垣地區(qū)鉛鋅礦床硫化物與賦礦圍巖灰?guī)r稀土元素球粒隕石標(biāo)準(zhǔn)化配分模式疊合圖Fig.4 Chondrite-normalized REE patterns of sulfides and host rock in the lead-zinc deposits in Huayuan area

三種礦石硫化物與近礦圍巖灰?guī)r雖均具有富輕稀土元素組成的特點(diǎn)，稀土配分曲線均為右傾型，但閃鋅礦、方鉛礦、黃鐵礦與圍巖灰?guī)r完全無重疊性，三種硫化物的稀土配分曲線呈鋸齒狀，推斷其成礦物質(zhì)可能部分來源于古老基底地層[26]?；?guī)r的稀土配分曲線則相對(duì)較為平滑，兩者區(qū)別明顯（圖3，圖4）。閃鋅礦存在中等程度的正Eu異常，方鉛礦存在弱的正Eu異常，個(gè)別樣品為弱的負(fù)Eu異常，黃鐵礦具有弱的正Eu異常，圍巖灰?guī)r則為中等程度的負(fù)Eu異常。閃鋅礦、方鉛礦和黃鐵礦三種硫化物均具有弱正Ce異常，圍巖灰?guī)r則存在中等程度的負(fù)Ce異常。兩者的δEu和δCe也具有一定的區(qū)別。硫化物中的稀土元素組成可以近似地代表成礦流體的稀土組成[30]。閃鋅礦、方鉛礦和黃鐵礦的正Eu異常和正Ce異常應(yīng)來源于對(duì)源巖的繼承或淋濾，圍巖清虛洞組灰?guī)r存在中等程度的負(fù)Eu異常和負(fù)Ce異常，因此，硫化物的正Eu異常和正Ce異常不可能完全來源于對(duì)圍巖灰?guī)r的繼承或淋濾。反映了花垣地區(qū)鉛鋅礦成礦與圍巖的無直接親緣關(guān)系，圍巖清虛洞組灰?guī)r為本區(qū)鉛鋅礦床提供成礦物質(zhì)的可能性不大，成礦物質(zhì)來源并非為含礦層位本身。

在硫化物與灰?guī)r的Sm/Nd-LREE/HREE變異圖（圖 5）和Sm/Nd-Tb/La變異圖（圖 6）上，硫化物與灰?guī)r不具明顯的相關(guān)性，沒有明顯的成因聯(lián)系，表明花垣地區(qū)鉛鋅礦床硫化物與灰?guī)r之間不具同源性，成礦物質(zhì)并非來源于含礦圍巖。

Y和Ho兩個(gè)元素具有相同的價(jià)態(tài)和離子半徑，它們8次配位的離子半徑也幾乎完全相同，因此該對(duì)元素在大部分的地球化學(xué)環(huán)境及其地球化學(xué)行為具有緊密的一致性[31]。Bau和Dulski[32]認(rèn)為Y/Ho值在不同類型巖漿巖、沉積巖及球粒隕石中無明顯的變化，球粒隕石中Y/Ho為24～36。花垣地區(qū)鉛鋅礦床灰?guī)r的Y/Ho值范圍為25.86～38.60，與球粒隕石的Y/Ho值較為接近，而該區(qū)礦石硫化物的Y/Ho值范圍為10.80～47.06，主要介于11～23，與賦礦圍巖灰?guī)r的Y/Ho值范圍差異較大，反映了硫化物與灰?guī)r之間的非同源性。另外，在La/Ho-Y/Ho圖解中（圖7），數(shù)據(jù)投點(diǎn)分布較為散亂，未體現(xiàn)出水平分布的特點(diǎn)，從而進(jìn)一步證實(shí)了花垣地區(qū)鉛鋅礦床硫化物與灰?guī)r之間的非同源性。

圖5 花垣地區(qū)鉛鋅礦床硫化物與灰?guī)r的Sm/Nd-LREE/HREE變異圖解Fig.5 Sm/Nd-LREE/HREE graphic of sulfides and hostrock in the lead-zinc deposits in Huayuan area

圖6 花垣地區(qū)鉛鋅礦床硫化物與灰?guī)r的Sm/Nd-Tb/La變異圖解Fig.6 Sm/Nd-Tb/La graphic of sulfides and hostrock in the lead-zinc deposits in Huayuan area

4.2 成礦環(huán)境與流體來源指示

礦石礦物的稀土元素地球化學(xué)特征可以代表成礦流體的稀土元素特征，反映出物質(zhì)來源、成礦流體的來源及演化、形成環(huán)境和介質(zhì)的性質(zhì)等方面的重要信息[29]，可成為解決物質(zhì)來源、成礦條件和成礦流體的有效手段[33]。

在水/巖反應(yīng)過程中，相對(duì)氧化條件下生成的礦物常出現(xiàn)負(fù)Eu異常、高ΣREE含量和低的ΣLREE/ΣHREE比值，相對(duì)還原條件下生成的礦物則出現(xiàn)正Eu異常、低ΣREE含量和高的ΣLREE/ΣHREE比值，這就是稀土元素地球化學(xué)演化的氧化還原模式[34]。鈰異常也取決于沉積環(huán)境的氧化還原條件，在氧逸度較高的條件下，成礦流體體系虧損Ce[29]，花垣礦集區(qū)內(nèi)的礦石礦物主要為閃鋅礦、方鉛礦和黃鐵礦，含礦熱液中富含高活動(dòng)性的S2－，ΣREE含量低，具有高的ΣLREE/ΣHREE比值，礦石硫化物δCe主要介于1.2～3.3，說明成礦流體處于相對(duì)還原環(huán)境，為還原性流體。

圖7 花垣地區(qū)鉛鋅礦床硫化物與灰?guī)r的La/Ho-Y/Ho圖解Fig.7 La/Ho-Y/Ho graphic of sulfides and hostrock in the lead-zinc deposits in Huayuan area

海底熱液及其沉積物球粒隕石標(biāo)準(zhǔn)化曲線具正Eu異常、明顯右傾的特征[35]?；ㄔ貐^(qū)鉛鋅礦床閃鋅礦、方鉛礦、黃鐵礦三種礦石硫化物球粒隕石標(biāo)準(zhǔn)化曲線具正Eu異常、弱負(fù)Ce異常、配分曲線右傾的特點(diǎn)，表明其成礦流體來源為具有海水混合的熱液流體。周云等[16]通過花垣地區(qū)鉛鋅礦床碳?xì)溲跬凰氐难芯刻岢龅V床成礦流體主要來源是大氣降水和經(jīng)過海水滲濾等作用演化形成的地層水，證明了其混合的海水熱液來源。

周云等[17]還研究了花垣鉛鋅礦床的Sr同位素地球化學(xué)特征，也論證了礦床成礦物質(zhì)來源并非為賦礦圍巖清虛洞組灰?guī)r。閃鋅礦礦物的鍶同位素初始比值（87Sr/86Sr）i值為 0.70912～0.70940，高于圍巖下寒武統(tǒng)清虛洞組灰?guī)r的鍶同位素比值平均值0.70904，認(rèn)為該礦床較高的鍶同位素初始比值反映了成礦物質(zhì)來源主要為陸殼，推斷成礦流體應(yīng)是流經(jīng)了下伏地層牛蹄塘組，與其中具有高鍶同位素比值的碎屑巖和泥巖等進(jìn)行了水巖反應(yīng)及同位素交換，從而導(dǎo)致沉淀出來的硫化物具有比賦礦圍巖高的鍶同位素組成。因此，花垣地區(qū)鉛鋅礦床的成礦物質(zhì)應(yīng)主要來源于下伏地層牛蹄塘組的陸源碎屑巖。

受加里東運(yùn)動(dòng)伸展構(gòu)造作用影響，封存于下覆地層中的由地層水、大氣降水和殘余海水混合形成的流體受構(gòu)造及熱動(dòng)力等的驅(qū)動(dòng)，發(fā)生大規(guī)模流動(dòng)，萃取了古老基底地層中的成礦物質(zhì)后，沿深大斷裂帶開始往上運(yùn)移，礦田內(nèi)發(fā)育的深大斷裂為深部熱液的大規(guī)模運(yùn)移和熱量傳輸提供了良好的通道[1]。在熱流體系統(tǒng)的活動(dòng)過程中，流體淋濾萃取底寒武統(tǒng)牛蹄塘組礦源層中的鉛鋅等成礦物質(zhì)，從而演化為成礦流體，形成成礦元素的最終富集。

5 結(jié)論

花垣地區(qū)鉛鋅礦床礦石硫化物的ΣREE含量低，相對(duì)富集輕稀土，具有高的ΣLREE/ΣHREE比值，與圍巖灰?guī)r相比，前者稀土元素總量明顯低于后者，Eu異常和Ce異常均不相同，稀土元素球粒隕石標(biāo)準(zhǔn)化配分模式也各不相同，礦石礦物與圍巖稀土元素組成具有較大的差異性，顯示該區(qū)含礦層并不是成礦物質(zhì)的主要來源，成礦物質(zhì)應(yīng)來源于古老基底和下伏牛蹄塘組的陸源碎屑巖。成礦流體為還原性流體，來源為具有海水混合的熱液流體。

[1]段其發(fā).湘西-鄂西地區(qū)震旦系-寒武系層控鉛鋅礦成礦規(guī)律研究[D].武漢:中國地質(zhì)大學(xué)(武漢)，2014.

[2]夏新階,舒見聞.李梅鋅礦床地質(zhì)特征及其成因 [J].大地構(gòu)造與成礦學(xué)，1995，19(3)：197-204.

[3]羅衛(wèi)，尹展，孔令，戴塔根.花垣李梅鉛鋅礦集區(qū)地質(zhì)特征及礦床成因探討 [J].地質(zhì)調(diào)查與研究，2009，33(3)：194-202.

[4]楊紹祥,勞可通.湘西北鉛鋅礦床的地質(zhì)特征及找礦標(biāo)志[J].地質(zhì)通報(bào)，2007，26(7)：899-908.

[5]毛黨龍.湖南省花垣縣大腦坡鉛鋅礦地質(zhì)特征及成因[J].現(xiàn)代礦業(yè)，2016，562(2)：36-45.

[6]彭能立，劉偉，梁恩云，劉庚寅，袁甫，熊苗，黃樂清.湘西地區(qū)北洛塔和花垣鉛鋅礦田對(duì)比研究[J].地質(zhì)找礦論叢，2016，31(2)：190-198.

[7]趙佳進(jìn)，許明珠，劉建平，唐宇薔，孔華，張強(qiáng)，楊哲.湘西花垣鉛鋅礦田的成礦構(gòu)造 [J].地質(zhì)找礦論叢，2016，31(3)：346-354.

[8]匡文龍，楊紹祥，余沛然，勞可通.湘西北花垣縣下寒武統(tǒng)清虛洞組濁積巖沉積特征及其地質(zhì)意義 [J].地質(zhì)科學(xué)，2008，43(2)：347-358.

[9]梁薇，牟傳龍，周懇懇，葛祥英.湘西花垣排碧寒武系花橋組上段—車夫組沉積環(huán)境的探討 [J].地學(xué)論評(píng)，2012，58（2）：159-175.

[10]薛長(zhǎng)軍，呂古賢，高偉利，楊陽.湘西花垣李梅礦田含礦層清虛洞期巖相古地理分析及成礦預(yù)測(cè) [J].地學(xué)前緣，2017，24（2）：159-175.

[11]張欣平，鄧華龍.花垣縣漁塘鉛鋅礦含礦層時(shí)代及早、中寒武世的三葉蟲[J].湖南地質(zhì)，1984，3(3)：194-202.

[12]朱才伐，李國祥，董熙平.湘西花垣中寒武統(tǒng)具表面裝飾的管狀微體化石[J].北京大學(xué)學(xué)報(bào)(自然科學(xué)版)，2004，40(2)：193-202.

[13]王辰，劉建朝，張海東，楊陽，楊飛，張旭陽.湘西花垣地區(qū)奧陶－志留系沉積巖稀土元素地球化學(xué)特征及地質(zhì)意義[J].礦物巖石地球化學(xué)通報(bào)，2017，36（3）：516-522.

[14]周云.湘西花垣MVT型鉛鋅礦集區(qū)成礦作用研究[D].成都:成都理工大學(xué)，2017.

[15]周云，段其發(fā)，陳毓川，唐菊興，曹亮，彭三國，甘金木.湘西花垣鉛鋅礦田成礦物質(zhì)來源的C、O、H、S、Pb、Sr同位素制約[J].地質(zhì)學(xué)報(bào)，2016，90(10)：2786-2802.

[16]周云，段其發(fā)，陳毓川，唐菊興，曹亮，彭三國，甘金木.湘西花垣地區(qū)鉛鋅礦床碳?xì)溲跬凰靥卣骷捌鋵?duì)成礦流體來源的指示[J].地質(zhì)通報(bào),2017，41(3):823-833.

[17]周云，段其發(fā)，曹亮，于玉帥，甘金木.湘西花垣鉛鋅礦田下寒武統(tǒng)清虛洞組灰?guī)r與熱液礦物的鍶同位素研究[J].地層學(xué)雜志,2017，41(3):335-343.

[18]周云,段其發(fā),唐菊興.湘西地區(qū)鉛鋅礦的大范圍低溫流體成礦作用-流體包裹體研究[J].地質(zhì)與勘探,2014，50(3)：515-532.

[19]劉文均,鄭榮才.花垣鉛鋅礦床成礦流體特征及動(dòng)態(tài)[J].礦床地質(zhì)，2000，19(2)：173-181.

[20]周皓迪，邵擁軍，葉周，魄含濤.湘西花垣大腦坡鉛鋅礦床流體包裹體研究[J].礦物學(xué)報(bào)，2015，35（增刊）：652.

[21]楊紹祥，勞可通.湘西北鉛鋅礦床碳?xì)溲跬凰靥卣骷俺傻V環(huán)境分析[J].礦床地質(zhì)，2007，26(3)：330-340.

[22]蔡應(yīng)雄,楊紅梅,段瑞春.湘西—黔東下寒武統(tǒng)鉛鋅礦床流體包裹體和硫、鉛、碳同位素地球化學(xué)特征[J].現(xiàn)代地質(zhì)，2014，28(1):29-41.

[23]李堃,吳昌雄,湯朝陽.湘西黔東地區(qū)鉛鋅礦床C、O同位素地球化學(xué)特征及其對(duì)成礦過程的指示 [J].中國地質(zhì)，2014，41(5)：1608-1619.

[24]段其發(fā),曹亮,曾健康.湘西花垣礦集區(qū)獅子山鉛鋅礦床閃鋅礦Rb-Sr定年及地質(zhì)意義 [J].地球科學(xué)-中國地質(zhì)大學(xué)學(xué)報(bào)，2014，39(8):977-999.

[25]王彥博，王敏芳，胡太平.湖南花垣獅子山礦床方解石稀土元素和碳氧同位素地球化學(xué)特征研究 [J].礦物學(xué)報(bào)，2015，35（增刊）：728-729.

[26]歐 正，韋重韜，莫紹文，高賀鳳，Santigie Kekuda Sesay，鄒明俊.廣西下巴鉛鋅礦床稀土元素地球化學(xué)特征研究[J].現(xiàn)代地質(zhì)，2013，27(4)：841-846.

[27]周云，段其發(fā)，曹亮，甘金木.湘西花垣鉛鋅礦稀土元素地球化學(xué)特征與指示意義[J].礦物學(xué)報(bào)，2015，35（增刊）：751.

[28]王中剛，于學(xué)元，趙振華.稀土元素地球化學(xué)[M].北京：科學(xué)出版社，1989，92.

[29]劉淑文，石順，李榮西，高云寶，劉玲芳，段立志，陳寶赟，張少妮.揚(yáng)子板塊北緣馬元鉛鋅礦床稀土元素地球化學(xué)研究[J].礦床地質(zhì)，2013，32(5)：979-988.

[30]顧尚義.黔西北鉛鋅礦稀土元素組成特征-兼論黔西北地區(qū)鉛鋅礦成礦與峨眉山玄武巖的關(guān)系 [J].貴州地質(zhì),2006，89(4)：274-277.

[31]鄧明國，林冰霞，梁徐文，李文昌.云南鎮(zhèn)康蘆子園鉛鋅多金屬礦床稀土元素地球化學(xué)特征及意義[J].中國稀土學(xué)報(bào),2014，32(1):117-128.

[32]Bau M,Dulski P.Anthropogenic origin of positive gadolinium anomalies in river waters[J].Earth Planet.Sci.Lett.，1996，143：245.

[33]Hanson G N.Rare earth elements in petrogenetic studies of igneous systems[J].Ann.Rev.Earth＆Planet Scil.，1980，8：371.

[34]周家喜,黃智龍,周國富.黔西北天橋鉛鋅礦床熱液方解石C、O同位素和REE地球化學(xué) [J].大地構(gòu)造與成礦學(xué)，2012，36(1)：93-101.

[35]陳寶赟，李榮西，劉淑文，張少妮，段立志，丁磊.馬元鉛鋅礦成礦流體類型-來自稀土元素的證據(jù)[J].礦床地質(zhì)，2012，31(增刊)：531-532.