中低溫成礦流體中氯配合物和硫氫配合物在銅金屬遷移中的行為

周冠軒, 王英濱*, 黃 靖, 張德會

中低溫成礦流體中氯配合物和硫氫配合物在銅金屬遷移中的行為

周冠軒1, 王英濱1*, 黃 靖1, 張德會2

(1. 中國地質(zhì)大學(xué)(北京) 數(shù)理學(xué)院, 北京 100083; 2. 中國地質(zhì)大學(xué)(北京) 地球科學(xué)與資源學(xué)院, 北京 100083)

熱液成礦系統(tǒng)是世界上最主要的金屬礦物來源, 其中的S和Cl是流體中最為重要的兩種配體, 對地殼內(nèi)金屬在流體中的遷移起到重要作用。在中低溫(25～300 ℃)的溫度范圍下, 根據(jù)最新的可用熱力學(xué)性質(zhì), 利用地球化學(xué)模擬軟件GWB (Geochemist’s Workbench), 對Cu-Fe-NaCl-H2O-S體系進(jìn)行了地球化學(xué)模擬, 比較了氯配合物和硫氫配合物在Cu金屬遷移當(dāng)中的作用。模擬流體的各組分濃度均以天然熱液的流體范圍相當(dāng), 該模型解釋了熱液中Cu對于S和Cl的配合模式。在較低的溫度下, 硫氫配合物是體系中主要的Cu遷移物種, 但這種優(yōu)勢也與所處的pH值和氧化還原環(huán)境關(guān)系密切。而較高溫度下則主要為氯配合物形態(tài)進(jìn)行遷移。作為硫氫配合物或氯配合物的輸送則決定了金屬富集的模式, 該模型可用于參考評估自然界中的熱液流體成礦過程。

成礦流體; 銅; 熱力學(xué); 地球化學(xué)模擬

0 引 言

熱液成礦系統(tǒng)是世界最主要的金屬來源, 成礦元素通過對源區(qū)的萃取、遷移和沉淀在熱液系統(tǒng)中富集[1–4]。1977年首次發(fā)現(xiàn)了海底排出的熱液[5]。從那時起, 大量的關(guān)于熱液流體的研究開始進(jìn)行, 而現(xiàn)在水熱溶液在各類金屬礦床的形成中起到的重要作用已經(jīng)得到了廣泛的認(rèn)可。可以通過流體包裹體所顯示的數(shù)據(jù), 來描述在相對簡單的H2O-NaCl-(CO2-H2S-CH4)體系中具有不同成礦溫度和地質(zhì)特征的沉積物的流體[6], 而在流體中的Cl和S元素則是被認(rèn)為金屬在其中進(jìn)行遷移最重要的配體。

對成礦流體研究清楚的表明了它是含有簡單Na離子和復(fù)雜分子種類的多組分電解質(zhì)溶液, 元素在成礦流體中形成配合物的類型可以通過軟硬酸堿理論進(jìn)行解釋, 這是影響配合物形式的內(nèi)在原因[7]。作為配合物配體的堿的硬度, F?> OH?> Cl?> S2?> HS?, 具有體積較大的并且電荷較低特點的堿, 可以被相對較軟的酸如親硫元素(Cu、Ag和Au等)和過渡金屬元素(Fe、Co、Ni、Pt和Pd等)極化; 而如親石元素(K、Na、Ca和Mg等)及含氧酸元素(CO2和SO3等)的相對較硬的酸, 可以極化體積較小并具較高電荷特點的堿。成礦流體是多組分電解質(zhì)溶液, 其中NaCl通常是主要的鹽組分。成礦流體中大量電解質(zhì)鹽與成礦流體系統(tǒng)中的許多過渡金屬元素的配合在遷移/沉積中起到重要作用[8], 而近年來的實驗研究也對Cl元素[9–11]和S元素[12–13]在金屬(Cu、Pb和Zn等)中的遷移、礦化的重要性有了共識。流體中的鹵素元素對金屬元素的溶解、活化和遷移具有重要的作用, 而其中的鹽度、溫度、酸堿度和氧化還原的特性都會對其中金屬離子的活化、遷移和富集等產(chǎn)生重要的影響。

使用數(shù)值模擬技術(shù)模擬自然環(huán)境中發(fā)生的地球化學(xué)過程已成為一種標(biāo)準(zhǔn)技術(shù), 化學(xué)熱力學(xué)一直應(yīng)用于地球化學(xué)的分析和建模, 程序中提供的數(shù)據(jù)庫包含許多水性物種, 其熱力學(xué)參數(shù)通過擬合經(jīng)驗數(shù)據(jù)或通過預(yù)測的方法得出。

本文擬主要根據(jù)現(xiàn)有的化學(xué)熱力學(xué)性質(zhì), 應(yīng)用地球化學(xué)模擬軟件GWB (Geochemist’s Workbench), 對Cu在不同水熱流體體系的中低溫?zé)嵋簵l件下具有不同物理化學(xué)條件的流體進(jìn)行模擬。設(shè)定在中低溫的溫度范圍(25～300 ℃)內(nèi), 探討Cu在不同水熱條件下主要存在形態(tài)以及水熱溶液當(dāng)中Cu配合物的濃度變化, 以期為預(yù)測熱液成礦過程提供依據(jù)。

1 方法與模擬過程

應(yīng)用GWB軟件模擬成礦流體在熱力學(xué)平衡下的化學(xué)反應(yīng)。GWB是最為重要的地球化學(xué)模擬軟件, 由美國Aqueous Solutions LLC公司開發(fā), 廣泛應(yīng)用在環(huán)境地球化學(xué)、油氣地球化學(xué)以及礦床地球化學(xué)等領(lǐng)域。運用計算機(jī)模擬眾多的復(fù)雜地球化學(xué)反應(yīng)過程, 通過量化模型來了解沉積成巖和熱液變化、開發(fā)探明礦床、確定污染物遷移規(guī)律、預(yù)測礦床和油藏物質(zhì)演化。GWB模擬軟件有三個版本, 即GWB Essentials、GWB Standard和GWB Professional, 本研究使用的是GWB Professional 11.0 版本。使用軟件模擬更加適用于難以通過實驗分析得到結(jié)果的領(lǐng)域, 但模擬需要建立在完善的熱力學(xué)數(shù)據(jù)的基礎(chǔ)上, 需要較為準(zhǔn)確的熱力學(xué)數(shù)據(jù)才能得到更加準(zhǔn)確的結(jié)果; 同時, GWB的通用熱力學(xué)參數(shù)在0～300 ℃和有限的壓力條件, 有一定的局限性。

Cl離子(Cl?)和硫氫根離子(HS?)是Cu在水熱流體當(dāng)中最重要的配體, 本研究的主要目的是依據(jù)已有的熱力學(xué)數(shù)據(jù), 模擬Cu在中低溫范圍的成礦流體中, 主要物種形態(tài)以及Cu配合物濃度受不同條件的影響變化。

Cu在天然流體中以兩種氧化態(tài)存在, 預(yù)測Cu(Ⅰ)配合物在流體中的遷移過程中占據(jù)主導(dǎo)地位[14–15]。為了探討中低溫成礦流體中氯配合物和硫氫配合物在銅金屬遷移中的行為, 所使用的的熱力學(xué)數(shù)據(jù), 在勞倫斯利弗莫爾國家實驗室(LLNL)默認(rèn)數(shù)據(jù)庫的基礎(chǔ)上, 選擇最新的熱力學(xué)數(shù)據(jù): Cu(Ⅰ)的氯配合物[16]、Cu(Ⅱ)的氯配合物[17]、以及Cu(Ⅰ)的硫氫配合物[18–19]將其添加至GWB模擬軟件當(dāng)中。選擇成礦流體的特征數(shù)據(jù)[20–21], 使用GWB軟件進(jìn)行模擬, 結(jié)果見表1。

2 結(jié)果與討論

2.1 中低溫鹽水熱液的主要形態(tài)

熱液成礦流體是多組分的電解質(zhì)溶液, 氯化鈉是熱液成礦流體中普遍存在且最主要的鹽組分, 許多元素的穩(wěn)定氯配合物通過流體遷移到整個地殼, 并隨著遷移過程中的理化性質(zhì)變化而沉淀, 對這些元素的遷移并形成礦床起到重要作用, 因此本文首先對最簡單的NaCl-H2O水熱條件下的流體進(jìn)行模擬, 構(gòu)建Eh-pH圖以證明在低鹽度的水熱條件下主要的Cu形態(tài)和礦物質(zhì)隨溫度的變化(圖1)。其中含水物質(zhì)的活度為:Cu=1×10?5mol/L;Na=Cl=0.6 mol/L。

從圖1中可以看出, 在氧化條件下, Cu均以Cu(Ⅱ)與Cl?和OH?形成配合物; 低溫條件下, pH > 6, 開始產(chǎn)生黑銅礦(CuO)沉淀, 隨著溫度的升高, 生成沉淀的pH值由6升至11; [CuCl2]?則是主要在還原環(huán)境中存在, 在100 ℃以上占據(jù)優(yōu)勢, 隨著溫度的升高, 優(yōu)勢區(qū)域所占pH值范圍變大; pH > 8時有赤銅礦(Cu2O)的沉淀生成, 隨著溫度的升高, 生成沉淀的pH值由8升至接近12?？梢缘贸? 溫度的升高更加促進(jìn)Cu在成礦流體中的遷移, 而pH值的升高則更加利于Cu在流體中的沉淀, 氧逸度則是影響Cu在成礦流體中的價態(tài), 從而影響Cu在成礦流體中的存在形式。因為Cu配合物的穩(wěn)定常數(shù)的溫度依賴性導(dǎo)致在冷卻的熱液中, 氯配合物物的穩(wěn)定性急劇下降, 會發(fā)生氯配合銅的極其有效的沉積, 特別是如果水巖相互作用, 導(dǎo)致pH同時增加。

表1 相關(guān)銅配合物的logk

圖1 低鹽度下主要Cu形態(tài)的Eh-pH圖

從圖2可以看出, 在室溫狀態(tài)的還原條件下, 低鹽度環(huán)境中以[CuCl2]?為主要的存在形式, 而[CuCl3]2?則在高鹽度環(huán)境中占據(jù)優(yōu)勢。也有計算表明, 主要存在的配合物的pH值范圍由0～9擴(kuò)大到0～11, 隨著pH的再升高, 會有赤銅礦沉淀生成; 在氧化環(huán)境中, [CuCl]+主要在pH=0～7范圍內(nèi)存在, pH > 7開始沉淀黑銅礦; 在300 ℃和低鹽度的還原態(tài)下, 仍是[CuCl2]?在體系中占優(yōu)勢, 而在氧化態(tài)中以[CuOH]+為體系的主要存在物種; 在高鹽度環(huán)境中, 在pH=0～2的高氧逸度環(huán)境中, CuCl2為體系內(nèi)的主要存在形態(tài); pH=2～10, [CuOH]+在體系內(nèi)占據(jù)優(yōu)勢。pH > 10開始沉淀出黑銅礦。而其他條件下, 均以[CuCl2]?為該環(huán)境的主要Cu形態(tài)。

圖2 不同鹽度環(huán)境下主要Cu形態(tài)的Eh-pH圖

2.2 中低溫含硫熱液的主要物種

2.3 Cu在流體中存在形式隨氯離子和硫化氫濃度的變化

圖3 含硫熱液中主要Cu形態(tài)的Eh-pH圖

Fig.3 Eh-pH diagram of the main Cu form in sulfur-containing hydrothermal fluid

2.4 影響Cu配合物濃度變化的因素

2.4.1 溫度對Cu配合物濃度的影響

圖4 Cu在一定范圍的Cl?和H2S(aq)活度下的主要存在形態(tài)

從圖5中可以看出, [Cu(HS)2]?在低溫下占據(jù)明顯的優(yōu)勢, 為體系當(dāng)中的主要成分。隨著溫度的升高, [CuCl2]?和[CuCl3]2?的濃度有著明顯的上升, 在250 ℃后濃度高于[Cu(HS)2]?。溫度從25 ℃變化到300 ℃, [CuCl2]?的濃度變化達(dá)到8個數(shù)量級, [CuCl3]2?的濃度變化約為6個數(shù)量級, 說明Cu的氯配合物受溫度影響較大。Cu在低溫的成礦流體中主要以硫氫配合物的形式進(jìn)行遷移, 隨著溫度的升高, 更多以氯配合物的形式進(jìn)行遷移。

圖5 Cu配合物濃度隨溫度的變化

2.4.2 氧逸度對Cu配合物濃度的影響

2.4.3 pH值對Cu配合物濃度的影響

圖6 Cu配合物濃度隨氧逸度的變化

圖7 Cu配合物濃度隨pH值的變化

從圖7中可以看出, 隨著pH值的升高, 體系內(nèi)整體Cu配合物濃度均有所升高, 推測一定的酸濃度利于Cu的礦物在體系內(nèi)的溶解, 與配離子形成配合物; 在pH > 4時, [CuCl2]?和[CuCl3]2?濃度開始出現(xiàn)下降, [Cu(HS)2]?濃度基本保持不變, 隨著pH值的增大, 當(dāng)pH > 5時, [Cu(HS)2]?的濃度開始逐漸升高并占據(jù)優(yōu)勢, 與pH=2的環(huán)境下[Cu(HS)2]?的濃度相差8個數(shù)量級。可以看出酸性pH的環(huán)境可以促進(jìn)Cu礦物在成礦流體當(dāng)中的溶解, 而pH的升高則有利于Cu形成[Cu(HS)2]?配合物。

2.4.4 氯離子濃度對Cu配合物濃度的影響

圖8 Cu配合物濃度隨氯離子濃度的變化

從圖8中可以看到, 氯離子濃度較低時, [Cu(HS)2]?為體系內(nèi)的主要配合物, 隨著氯離子濃度從0.6 mol/L升至約1.5 mol/L, Cu的硫氫配合物的濃度緩慢下降, 氯配合物的濃度逐步上升; 在氯離子濃度超過2 mol/L時, [CuCl2]?成為水熱溶液中的優(yōu)勢物種; 氯離子濃度的變化對氯配合物的影響則更大, [CuCl2]?增長了1個數(shù)量級, 而[CuCl3]2?濃度升高了2個數(shù)量級, 硫氫配合物的濃度變化較小。

2.4.5 H2S活度對Cu配合物濃度的影響

從圖9中可以看出H2S活度變化對體系氯配合物幾乎無影響, 而硫氫配合物的濃度則隨著H2S活度的增加得到增長, 在H2S(aq)活度約為0.02 mol/L時[Cu(HS)2]?的濃度高于[CuCl2]?。[Cu(HS)2]?的濃度變化約為1個數(shù)量級。而將圖8和圖9綜合解讀, 可以得到氯離子是在成礦流體當(dāng)中主要控制Cu遷移的配體, 而HS?配體雖然不能達(dá)到氯離子的影響, 但它足以在異常銅礦化的形成當(dāng)中發(fā)揮重要的作用。

圖9 Cu配合物濃度隨H2S活度的變化

3 與天然流體中金屬濃度變化的比較

為了檢查模擬的準(zhǔn)確性, 將模型流體和天然流體的Cu濃度與溫度和鹽度的變化趨勢進(jìn)行了比較; 天然流體數(shù)據(jù)使用文獻(xiàn)中LA-ICP-MS獲得的單個流體包裹體中的金屬濃度[26–30]; 這些天然流體的Cu濃度與溫度、鹽度關(guān)系的散點圖分別繪制在圖10和圖11中。

從圖10看出, 從同一地區(qū)獲得的流體包裹體中Cu的濃度與溫度的變化呈現(xiàn)正相關(guān)關(guān)系的趨勢, 這種現(xiàn)象也與本文中模擬的結(jié)果一致, 但也會出現(xiàn)一些不太符合趨勢的點, 這與流體當(dāng)中存在的鹽度影響有關(guān); 但不同地區(qū)獲得的流體包裹體數(shù)據(jù)隨溫度變化關(guān)系并不明顯, 這可能與流體當(dāng)中存在的其他元素對Cu濃度的影響有關(guān)。

同樣, 從圖11中可以得到, 同一地區(qū)中獲得的流體包裹體數(shù)據(jù)顯示, 隨著鹽度的升高, 整體Cu的濃度變化也出現(xiàn)了升高的趨勢, 這與模擬結(jié)果一致; 但不同地區(qū)中Cu濃度之間比較, 并不會出現(xiàn)非常明顯濃度趨勢, 這同樣也與礦區(qū)類型、一起流體當(dāng)中存在的其他元素對Cu濃度造成的影響有關(guān)。

4 Cu與兩種配體的配合形式

氯離子和硫氫根離子被認(rèn)為是控制成礦流體中銅遷移最重要的配體。銅在自然界以兩種氧化態(tài)的形式存在, 即Cu(Ⅰ)和Cu(Ⅱ)。Cu(Ⅰ)是地殼流體當(dāng)中最為常見的Cu氧化態(tài), Cu(Ⅰ)的100型電子構(gòu)型使其具有獨特的性質(zhì), Cu(Ⅰ)可以形成線性、三角形或四面體結(jié)構(gòu)配合物, 由于它們之間能級相差很低,使得它們可以在配合物的混合物中共存。在室溫下, Cu(Ⅰ)離子不穩(wěn)定, 但可以形成含水配合物穩(wěn)定Cu(Ⅰ)的氧化態(tài), 含水Cu(Ⅰ)具有線性的結(jié)構(gòu)[H2O-Cu-OH2]+ [31]。Cu(Ⅰ)可以與配體形成穩(wěn)定的配合物, 這些配合物根據(jù)配體類型和濃度的不同顯示出不同的幾何形狀。在低鹽度溶液中, Cu(Ⅰ)保留原有的線性結(jié)構(gòu)(如, [H2O-Cu-Cl](aq), [Cl-Cu-Cl]?)[32]; 而在高鹽度溶液中, 具有三角平面結(jié)構(gòu)的[CuCl3]2?變得穩(wěn)定[33]。隨著溫度的升高, 線性[CuCl2]?即使在高鹽度下也占據(jù)主導(dǎo)地位[34]。這是由于溫度和鹽度是影響金屬絡(luò)合物配位變化的兩個關(guān)鍵因素, 溫度的升高則改變了水和氫鍵的網(wǎng)絡(luò)結(jié)構(gòu), 從而改變了金屬離子的溶劑化性質(zhì)和絡(luò)合物的幾何形狀。此外Cu對Cl具有很高的親和力, 這可以通過軟硬酸解理論解釋。Cu(Ⅱ)需要在近地表的高度氧化條件環(huán)境中才能產(chǎn)生, 而Cu(Ⅱ)的水合離子是可溶的, 與氯離子的絡(luò)合會發(fā)生在很高鹽度的條件下。Cu(Ⅱ)的配位化學(xué)以強(qiáng)的楊-特勒效應(yīng)為主, 導(dǎo)致配合物以多幾何形式排列。

圖10 天然流體中Cu濃度與溫度關(guān)系的散點圖

圖11 天然流體中Cu濃度與鹽度關(guān)系的散點圖

Cu(Ⅰ)被歸類為路易斯酸, 可以與HS?形成強(qiáng)鍵。與氯配合物類似, Cu(Ⅰ)硫氫配合物在整個的模擬條件下具有線性的幾何形狀[35]。但可以發(fā)現(xiàn)硫氫化銅配合物的穩(wěn)定性要比氯化銅配合物的穩(wěn)定性高出幾個數(shù)量級, 但兩者的溫度依賴性有顯著的不同。在冷卻的熱液成礦流體中, 氯配合物的穩(wěn)定性急劇下降, 會發(fā)生Cu的有效沉積。相反, 硫氫配合物的穩(wěn)定性隨著溫度的降低而增加, 并且隨著pH的升高, 可以增加Cu的硫氫配合物的在體系內(nèi)的濃度。

5 結(jié) 論

本研究基于檢索到的文獻(xiàn)實驗數(shù)據(jù)和可用的熱力學(xué)數(shù)據(jù), 使用地球化學(xué)模擬軟件GWB模擬在地殼上金屬Cu在熱液成礦流體中遷移行為。模擬結(jié)果表明了成礦流體模型在所處環(huán)境下攜帶金屬遷移的能力。對于氯配合物來說則是在酸性的高溫條件下主要以[CuCl2]?進(jìn)行遷移, 鹽度對氯配合物的濃度產(chǎn)生影響, 在高的鹽度、較低的溫度環(huán)境下則主要以[CuCl3]2?進(jìn)行遷移。硫氫配合物在較低溫度的還原及堿性環(huán)境下發(fā)揮關(guān)鍵作用, 使Cu主要以[Cu(HS)2]?形式在成礦流體當(dāng)中遷移, 其中鹽度顯著影響了硫氫配合物的優(yōu)勢濃度范圍; 而氧逸度的升高, 即是流體越接近地表, 則主要以CuSO4和[CuCl2]?的方式進(jìn)行遷移。模擬的Cu流體模型可與天然地質(zhì)流體模型比較, 對同一礦區(qū)呈現(xiàn)比較一致的變化趨勢, 但對于不礦區(qū)獲得的數(shù)據(jù), 與模擬結(jié)果產(chǎn)生一定差異, 這也與不同礦區(qū)中的元素組成有關(guān)。該模型流體可以用于預(yù)測熱液成礦的過程中Cu在流體當(dāng)中的遷移方式。

[1] 胡瑞忠, 溫漢捷, 蘇文超, 彭建堂, 畢獻(xiàn)武, 陳佑緯. 礦床地球化學(xué)近十年若干研究進(jìn)展[J]. 礦物巖石地球化學(xué)通報, 2014, 33(2): 127–144.

Hu Rui-zhong, Wen Han-jie, Su Wen-chao, Peng Jian-tang, Bi Xian-wu, Chen You-wei. Some advances in ore deposit geochemistry in last decade[J]. Bull Mineral Petrol Geochem, 2014, 33(2): 127–144 (in Chinese with English abstract).

[2] 陳毓川, 裴榮富, 王登紅, 王平安. 論礦床的自然分類——四論礦床的成礦系列問題[J]. 礦床地質(zhì), 2015, 34(6): 1092– 1106.

Chen Yu-chuan, Pei Rong-fu, Wang Deng-hong, Wang Ping-An. Natural classification of mineral deposits: Discussion on minerogenetic series of mineral deposits (Ⅳ)[J]. Mineral Deposits, 2015, 34(6): 1092–1106 (in Chinese with English abstract).

[3] 翟裕生. 試論礦床成因的基本模型[J]. 地學(xué)前緣, 2014, 21(1): 1–8.

Zhai Yu-sheng. A preliminary discussion on fundamental model of metallogenic mechanism[J]. Earth Sci Front, 2014, 21(1): 1–8 (in Chinese with English abstract).

[4] 鄧軍, 高幫飛, 王慶飛, 楊立強(qiáng). 成礦流體系統(tǒng)的形成與演化[J]. 地質(zhì)科技情報, 2005, 24(1): 49–54.

Deng Jun, Gao Bang-fei, Wang Qing-fei, Yang Li-qiang. Formation and evolution of ore-forming fluid system[J]. Geol Sci Tech Inf, 2005, 24(1): 49–54 (in Chinese with English abstract).

[5] Weiss R F, Lonsdale P, Lupton J E, Bainbridge A E, Craig H. Hydrothermal plumes in the Galapagos Rift[J]. Nature, 1977, 267(5612): 600–603.

[6] 張文淮, 張志堅, 伍剛. 成礦流體及成礦機(jī)制[J]. 地學(xué)前緣, 1996, 13(4): 86–93.

Zhang Wen-huai, Zhang Zhi-jian, Wu Gang. Ore-forming fluid and mineralization mechanism[J]. Earth Sci Front, 1996, 13(4): 86–93.

[7] 張德會. 關(guān)于成礦流體地球化學(xué)研究的幾個問題[J]. 地質(zhì)地球化學(xué), 1997, 25(3): 49–57.

Zhang De-hui. Some new advances in ore-forming fluid geochemistry[J]. Geol Geochem, 1997, 25(3): 49–57.

[8] 劉家軍, 柳振江, 楊艷, 石龍, 趙百勝, 毛光劍, 王建平. 鹽類在金屬、非金屬成礦過程中的作用[J]. 地質(zhì)找礦論叢, 2007, 22(3): 161–171.

Liu Jia-jun, Liu Zhen-jiang, Yang Yan, Shi Long, Zhao Bai-sheng, Mao Guang-jian, Wang Jian-ping. The salt roles in the mineralizlng processes of metal and non-metal depostis[J]. Contrib Geol Mineral Resour Res, 2007, 22(3): 161–171.

[9] Lukanin O A, Ryzhenko B N, Kurovskaya N A. Zn and Pb solubility and speciation in aqueous chloride fluids at T-P parameters corresponding to granitoid magma degassing and crystallization[J]. Geochem Int, 2013, 51(10): 802–830.

[10] Borg S J, Liu W. An XAS study of zinc speciation in aqueous acetate solutions at 25–200℃[J]. Nucl Instrum Method Phys Res A, 2010, 619(1): 276–279.

[11] Mei Y, Liu W, Sherman D M, Brugger J. Metal complexation and ion hydration in low density hydrothermal fluids: Ab initio molecular dynamics simulation of Cu(I) and Au(I) in chloride solutions (25–1000℃, 1–5000bar)[J]. Geochimica et Cosmochimica Acta, 2014, 131: 196–212.

[12] Mei Y, Sherman D M, Liu W, Brugger J. Ab initio molecular dynamics simulation and free energy exploration of copper(I) complexation by chloride and bisulfide in hydrothermal fluids[J]. Geochimica et Cosmochimica Acta, 2013, 102: 45–64.

[13] Tagirov B R, Seward T M. Hydrosulfide/sulfide complexes of zinc to 250℃ and the thermodynamic properties of sphalerite[J]. Chemical Geology, 2010, 269(3): 301–311.

[14] Liu W, Mcphail D C. Thermodynamic properties of copper chloride complexes and copper transport in magmatic-hydrothermal solutions[J]. Chemical Geology, 2005, 221(1): 21–39.

[15] Nagaseki H, Hayashi K. Experimental study of the behavior of copper and zinc in a boiling hydrothermal system[J]. Geology, 2008, 36(1): 27–30.

[16] Brugger J, Etschmann B, Liu W, Testemale D, Hazemann J L, Emerich H, van Beek W, Proux O. An XAS study of the structure and thermodynamics of Cu(I) chloride complexes in brines up to high temperature (400℃, 600bar)[J]. Geochimica et Cosmochimica Acta, 2007, 71(20): 4920–4941.

[17] Sverjensky D A, Shock E L, Helgeson H C. Prediction of the thermodynamic properties of aqueous metal complexes to 1000℃ and 5 kb[J]. Geochimica Et Cosmochimica Acta. 1997, 61(7): 1359–1412.

[18] Mountain B W, Seward T M. Hydrosulfide/sulfide complexes of copper(I) : Experimental confirmation of the stoichiometry and stability of Cu(HS)2? to elevated temperatures[J]. Geochimica Et Cosmochimica Acta, 2003, 67(16): 3005–3014.

[19] Mountain B W, Seward T M. The hydrosulphide/sulphide complexes of copper(I): experimental determination of stoichiometry and stability at 22℃ and reassessment of high temperature data[J]. Geochimica et Cosmochimica Acta. 1999, 63(1): 11–29.

[20] Zhong R C, Brugger J, Chen Y J, Li W B. Contrasting regimes of Cu, Zn and Pb transport in ore-forming hydrothermal fluids[J]. Chemical Geology, 2015, 395: 154–164.

[21] Barnes S, Liu W. Pt and Pd mobility in hydrothermal fluids: evidence from komatiites and from thermodynamic modelling[J]. Ore Geology Reviews, 2012, 44(1): 49–58.

[22] 宋國學(xué), 秦克章, 李光明, 王樂. 中硫型淺成低溫?zé)嵋航鸲嘟饘俚V床基本特征、研究進(jìn)展與展望[J]. 巖石學(xué)報, 2018, 34(3): 748–762.

SONG Guo-xue, QIN Ke-zhang, LI Guang-ming, WANG Le, Basic characteristics and research progresses of intermediate sulfidation type epithermal gold poly-metallic deposits, and prospects[J]. Acta Petrologica Sinica.2018, 34(3): 748–762.

[23] 王小林. 壓力和流體組分對硫酸鹽——水體系高溫液—液相分離的影響: 作用機(jī)制和地質(zhì)意義[C]. 南昌: 2017.

WANG Xiao-lin. Effect of Pressure and Fluid Composition on Sulfate-Water System High Temperature Liquid-Liquid Phase Separation: Mechanism of Action and Geological Significance[C]. Nanchang, Jiangxi, China: 2017.

[24] 萬良?xì)J, 喻軍敏, 楊永勝, 鐘源. 淺成低溫?zé)嵋嘿F金屬礦床成礦條件淺析[J]. 西部探礦工程, 2016, 28(9): 133–135.

WAN Liang-qin, Yu Jun-min, YANG Yong-sheng, ZHONG Yuan. Analysis on metallogenic conditions of epithermal depositsin epithermal deposits[J]. West-China Exploration Engineering, 2016, 28(09): 133–135.

[25] Seo J H, Guillong M, Heinrich C A. The role of sulfur in the formation of magmatic–hydrothermal copper–gold deposits[J]. Earth and Planetary Science Letters, 2009, 282(1): 323–328.

[26] Beuchat S, Moritz R, Pettke T. Fluid evolution in the W–Cu–Zn–Pb San Cristobal vein, Peru: fluid inclusion and stableisotope evidence[J]. Chemical Geology, 2004, 210(1): 201–224.

[27] Rusk B G, Reed M H, Dilles J H, Klemm L M, Heinrich C A. Compositions of magmatic hydrothermal fluids determined by LA-ICP-MS of fluid inclusions from the porphyry copper– molybdenum deposit at Butte, MT[J]. Chemical Geology, 2004, 210(1): 173–199.

[28] Kotzeva B G, Guillong M, Stefanova E, Piperov N B. LA-ICP-MS analysis of single fluid inclusions in a quartz crystal (Madan ore district, Bulgaria)[J]. Journal of Geochemical Exploration, 2011, 108(3): 163-175.

[29] Appold M S, Wenz Z J. Composition of Ore Fluid Inclusions from the Viburnum Trend, Southeast Missouri District, United States: Implications for Transport and Precipitation Mechanisms[J]. Economic Geology, 2011, 106(1): 55–78.

[30] Catchpole H, Kouzmanov K, Fontboté L, Guillong M, Heinrich C A. Fluid evolution in zoned Cordilleran polymetallic veins — Insights from microthermometry and LA-ICP-MS of fluid inclusions[J]. Chemical Geology, 2011, 281(3): 293–304.

[31] Blumberger J, Bernasconi L, Tavernelli I, Vuilleumier R, Sprik M. Electronic Structure and Solvation of Copper and Silver Ions: A Theoretical Picture of a Model Aqueous Redox Reaction[J]. Journal of the American Chemical Society, 2004, 126(12): 3928–3938.

[32] Liu W, Brugger J, Etschmann B, Testemale D, Hazemann J. The solubility of nantokite (CuCl(s)) and Cu speciation in low-density fluids near the critical isochore: An in-situ XAS study[J]. Geochimica et Cosmochimica Acta, 2008, 72(16): 4094–4106.

[33] Sherman D M. Complexation of Cu+ in Hydrothermal NaCl Brines: Ab initio molecular dynamics and energetics[J]. Geochimica et Cosmochimica Acta, 2007, 71(3): 714–722.

[34] Fulton J L, Hoffmann M M, Darab J G. An X-ray absorption fine structure study of copper(I) chloride coordination structure in water up to 325℃[J]. Chemical Physics Letters, 2000, 330(3): 300–308.

[35] Etschmann B E, Liu W, Testemale D, Müller H, Rae N A, Proux O, Hazemann J L, Brugger J. An in situ XAS study of copper(I) transport as hydrosulfide complexes in hydrothermal solutions (25–592 ℃, 180–600 bar): Speciation and solubility in vapor and liquid phases[J]. Geochim Cosmochim Acta, 2010, 74(16): 4723–4739.

Behavior of chlorine complexes and sulfur hydrogen complexes transport of copper in medium and low temperature in ore-forming fluids

ZHOU Guan-xuan1, WANG Ying-bin1*, HUANG Jing1and ZHANG De-hui2

1. School of Science, China University of Geosciences (Beijing), Beijing 100083, China; 2. School of Earth Sciences and Resources, China University of Geosciences (Beijing), Beijing 100083, China

Hydrothermal mineralization system is the world’s most important source of metal minerals. Sulfur and chlorine are the two most decisive ligands in the fluid, which play an important role in the migration of metals in the earth’s crust. Geochemical simulation of Cu-Fe-NaCl-H2O-S system bases on the latest available thermodynamic data using geochemical software The Geochemist’s Workbench at medium and low temperature (25 to 300 ℃) comparing the chlorine complexes and sulfur-hydrogen complexes in Cu metal migration. The concentration of each component of the simulated fluid is equivalent to the range of the fluid of the natural hydrothermal fluid. This model explains the mode of coordination of Cu for sulfur and chlorine in the hydrothermal fluid. At lower temperatures, the sulfur-hydrogen complex is the main Cu migration species in the system, but this advantage is also closely related to the pH and redox environment. At higher temperatures, migration occurs primarily in the form of chlorine complexes. The transport of sulfur-hydrogen complexes or chlorine complexes determines the mode of metal enrichment, which can be used as reference for the assessment of hydrothermal fluid mineralization processes in nature.

ore-forming fluid; copper; thermodynamics; geochemical model

P595

A

0379-1726(2021)04-0354-11

10.19700/j.0379-1726.2021.04.003

2019-09-02;

2019-10-28;

2019-11-11

國土資源部公益性行業(yè)科研專項(201411024-4)

周冠軒(1995–), 男, 碩士研究生, 化學(xué)工程專業(yè)。E-mail: 562425446@qq.com

WANG Ying-bin, E-mail: wyb@cugb.edu.cn; Tel: +86-13641309381