黃保有 張波** 張進(jìn)江 譚舟 張磊 陳思雨 李曉蓉

1. 造山帶與地殼演化教育部重點實驗室，北京大學(xué)地球與空間科學(xué)學(xué)院，北京 1008712. 巖石圈演化國家重點實驗室，中國科學(xué)院地質(zhì)與地球物理研究所，北京 100029

基于變質(zhì)巖石內(nèi)特定礦物組合估算變質(zhì)溫度和P-T-t軌跡是揭示造山帶熱演化過程的重要途徑(Spear, 1995；魏春景等，2001；Wuetal., 2002；Liuetal., 2004；劉福來等，2006；吳春明等，2007；Weietal., 2009)，然而，在很多變質(zhì)帶中常缺乏可用于估算變質(zhì)溫度的合適礦物組合，碳質(zhì)物拉曼光譜溫度計可作為一個新的方法給予補(bǔ)充(Beyssacetal., 2002a；Rahletal., 2005；Aoyaetal., 2010)。碳質(zhì)物廣泛存在于變質(zhì)沉積巖中，在成巖和變質(zhì)過程中逐漸轉(zhuǎn)變?yōu)槭?，即石墨?Buseck and Beyssac, 2014)，在石墨化過程中，碳質(zhì)物的結(jié)晶度逐漸提高，是巖石變質(zhì)程度(尤其是溫度)的可靠指標(biāo)。

對于碳質(zhì)物結(jié)晶度與溫度關(guān)系研究，早期的學(xué)者通過X射線衍射(XRD)技術(shù)進(jìn)行觀測(Itaya, 1981；Yuietal., 1996；Beyssacetal., 2002a；Rahletal., 2005；Aoyaetal., 2010；Lahfidetal., 2010；Kouketsuetal., 2014)，隨后通過鏡質(zhì)體反射率(French, 1964；Landis, 1971；Kisch, 1980；Itaya, 1981；Tagiri, 1981；Mori and Taguchi, 1988；Georgeetal., 2001)以及利用高分辨率透射電子顯微鏡(Buseck and Huang, 1985；Jehlicka and Rouzaud, 1990；Beyssacetal., 2002b；Nakamura and Akai, 2013)，這些研究結(jié)果表明碳質(zhì)物結(jié)晶度與變質(zhì)溫度之間具有正相關(guān)性。近年來，高空間分辨率激光拉曼光譜定量分析技術(shù)已成為表征碳質(zhì)物結(jié)晶度的有效工具(Beyssacetal., 2002a)。Beyssacetal.(2002a)利用該技術(shù)測量變質(zhì)巖中碳質(zhì)物的結(jié)晶度，并提出碳質(zhì)物拉曼光譜溫度計，由于方法簡單實用，引起廣泛關(guān)注(d’Alessioetal., 2003；Satoetal., 2006；Yamadaetal., 2007；Robertetal., 2010；Tagami, 2012；Cooperetal., 2013； Yoshidaetal., 2016；Laurentetal., 2018)。隨后有研究者基于實驗和天然樣品的分析，對該溫度計進(jìn)行了修訂和擴(kuò)展，提出了不同版本的變質(zhì)溫度計算公式，拓寬了該方法的使用范圍(Rahletal., 2005；Aoyaetal., 2010；Lahfidetal., 2010；Kouketsuetal., 2014)。作為新型溫度計，其影響因素也被關(guān)注，在低溫變質(zhì)環(huán)境(溫度<300℃)，碳質(zhì)物的結(jié)構(gòu)不僅受溫度影響，還受變質(zhì)壓力、變質(zhì)持續(xù)時間、變形、催化物質(zhì)和碳質(zhì)物前體類型的影響(Landis, 1971；Diesseletal., 1978；Bonijolyetal., 1982；Wopenka and Pasteris, 1993；Wadaetal., 1994；Bustinetal., 1995；Suchyetal., 1997；Buseck and Beyssac, 2014)。由于很難從天然巖石樣品中評估每一種因素的影響(Lahfidetal., 2010)，因此目前尚未建立統(tǒng)一的碳質(zhì)物拉曼光譜溫度計計算公式。盡管如此，激光拉曼光譜可以原位、準(zhǔn)確、高效的獲取碳質(zhì)物結(jié)晶度相關(guān)參數(shù)進(jìn)而估算變質(zhì)峰期溫度，解決了傳統(tǒng)方法繁瑣、費時等問題，大大提高了工作效率(Beyssacetal., 2002a)。相比基于特定礦物組合和礦物組成的溫度計(魏春景等，2001；Wuetal., 2002；吳春明等，2007；Weietal., 2009)，拉曼光譜溫度計在以下幾方面極具優(yōu)勢：(1)碳質(zhì)物廣泛存在于變質(zhì)沉積巖中，巖片及巖石薄片均可用于拉曼光譜分析(Beyssacetal., 2002a；Rahletal., 2005；Aoyaetal., 2010)；(2)激光拉曼光譜可對巖石薄片中的碳質(zhì)物進(jìn)行原位無損的探測分析，可同時觀測到碳質(zhì)物與周圍礦物間的結(jié)構(gòu)、構(gòu)造與組合關(guān)系，便于將構(gòu)造與變質(zhì)溫度與事件相關(guān)聯(lián)(Kouketsuetal., 2014)；(3)激光拉曼光譜分析準(zhǔn)確度高(分辨率為～1μm)、效率高(單點光譜獲取時間～10s)，因此激光拉曼可在短時間內(nèi)獲取大量、乃至海量光譜測量數(shù)據(jù)，以消除測量樣品的異質(zhì)性(Aoyaetal., 2010)；(4)變質(zhì)過程中碳質(zhì)物石墨化具有不可逆性，碳質(zhì)物結(jié)構(gòu)記錄峰期變質(zhì)，不會遭受后期退變質(zhì)改造(Buseck and Beyssac, 2014)，因此碳質(zhì)物拉曼光譜變質(zhì)溫度計估算的溫度為峰期變質(zhì)溫度(Beyssacetal., 2002a；Rahletal., 2005；Aoyaetal., 2010)。

1 碳質(zhì)物拉曼光譜溫度計原理與方法

變質(zhì)沉積巖中碳質(zhì)物隨變質(zhì)溫度增加而逐漸石墨化是碳質(zhì)物拉曼溫度計的理論基礎(chǔ)(Beyssacetal., 2002a, 2004；Rantitschetal., 2004)。沉積巖通常含有痕量的有序度低、結(jié)晶程度低的碳質(zhì)物，但隨變質(zhì)程度的增加，碳質(zhì)物逐漸轉(zhuǎn)變?yōu)橛行蚨雀叩慕Y(jié)晶石墨(Grew, 1974；Pasteris and Wopenka, 1991；Wopenka and Pasteris, 1993；Yuietal., 1996)。通過激光拉曼光譜獲取碳質(zhì)物結(jié)晶程度相關(guān)量化參數(shù)，再根據(jù)結(jié)晶度與溫度間的相關(guān)公式(Beyssacetal., 2002a；Rahletal., 2005；Aoyaetal., 2010；Lahfidetal., 2010；Kouketsuetal., 2014)，可估算變質(zhì)巖石變質(zhì)溫度(Beyssacetal., 2002a；Rahletal., 2005)，值得注意的是，由于石墨化過程的不可逆性(碳質(zhì)物傾向于熱力學(xué)穩(wěn)定相(石墨))，即碳質(zhì)物結(jié)構(gòu)對退變質(zhì)不敏感，因此利用該方法計算出的溫度可以被解釋為樣品整個變質(zhì)過程中達(dá)到的變質(zhì)峰期溫度(Landis, 1971；Pasteris and Wopenka, 1991；Wopenka and Pasteris, 1993；Wadaetal., 1994；Beyssacetal., 2002a, b, 2004, 2007；Rahletal., 2005；Negroetal., 2006；Angiboustetal., 2009；Gabaldaetal., 2009；Lahfidetal., 2010)。

1.1 碳質(zhì)物拉曼溫度計原理

拉曼光譜是一種散射光譜，當(dāng)激光發(fā)出的入射光子與作為散射中心的樣品分子相互作用時會產(chǎn)生彈性碰撞與非彈性碰撞，在非彈性碰撞過程中，入射光子與分子之間存在能量交換，導(dǎo)致入射光子方向與頻率發(fā)生改變，這種散射稱之為拉曼散射。拉曼散射依據(jù)樣品分子原始能級狀態(tài)又可進(jìn)一步分為斯托克斯散射與反斯托克斯散射(圖1)。根據(jù)玻爾茲曼定律，在室溫下，分子絕大多數(shù)處于振動能級基態(tài)，所以拉曼光譜記錄的一般都是斯托克斯散射過程(Reed and Simon, 1979；賀福，2005；張延會等，2006)。

圖1 拉曼散射示意圖Fig.1 Sketch diagram of Raman scattering

在拉曼光譜中橫坐標(biāo)拉曼位移是斯托克斯和反斯托克斯散射光頻率與入射光頻率之差(△v)，拉曼位移取決于分子振動能級的變化，它與入射光頻率無關(guān)，不同的化學(xué)鍵或基態(tài)有不同的振動方式，決定了其能級間的能量變化，因此與之對應(yīng)的拉曼位移也具特征(賀福，2005；張延會等，2006)，這是拉曼光譜進(jìn)行分子結(jié)構(gòu)定性分析的理論依據(jù)。拉曼光譜的縱軸代表譜峰強(qiáng)度，譜峰的強(qiáng)度與產(chǎn)生譜峰特定物質(zhì)的濃度有關(guān)(Tuinstra and Koenig, 1970；張延會等，2006)。

石墨為層狀結(jié)構(gòu)，共有六支色散曲線(圖2a)(Malardetal., 2009)，分別為三個光學(xué)支(面內(nèi)縱向光學(xué)支iLO、面內(nèi)橫向光學(xué)支iTO和面外橫向光學(xué)支oTO)和三個聲學(xué)支(面內(nèi)縱向聲學(xué)支iLA、面內(nèi)橫向聲學(xué)支iTA和面外橫向聲學(xué)支 oTA)。面內(nèi)(i)和面外(o)分別為原子的振動方向平行或者垂直于石墨烯平面，縱向(L)和橫向(T)即為原子的振動方向平行或者垂直于A-B碳碳鍵的方向(吳娟霞等，2014)。這些晶格振動產(chǎn)生的色散曲線表示出各種振動模式的頻率與波矢散射的關(guān)系(圖2a)。

碳質(zhì)物典型的拉曼光譜由一級序區(qū)(1100～1800cm-1)和二級序區(qū)(2500～3100cm-1)兩部分組成(Tuinstra and Koenig, 1970；Nemanich and Solin, 1979)(圖2c)。在一階序區(qū)，碳質(zhì)物拉曼光譜特征峰包括D1(D)峰(1350cm-1)、G峰(1580cm-1)和D2(D’)峰(1620cm-1)，結(jié)晶度低的碳質(zhì)物中通常還有D3峰(1520cm-1)和D4峰(1245cm-1)，其中，G峰代表石墨峰，D峰代表缺陷峰。二階序區(qū)對應(yīng)于倍頻和散射組合，主要為S1峰(2700cm-1)和S2峰(2900cm-1)(Nemanich and Solin, 1979；Wopenka and Pasteris, 1993)。

圖2 碳質(zhì)物拉曼光譜產(chǎn)生原理(a)單層石墨烯振動曲線(據(jù)吳娟霞等，2014修改)；(b)石墨晶體空間群振動模式及振動頻率(據(jù)賀福，2005)；(c)碳質(zhì)物典型拉曼光譜(據(jù)吳娟霞等，2014修改)；(d)碳質(zhì)物各個拉曼特征峰產(chǎn)生過程(據(jù)吳娟霞等，2014)Fig.2 Formation principle of carbon material Raman spectroscopy(a) phonon dispersion relation in graphite (modified after Wu et al., 2014)；(b) vibration mode and frequency of graphite group (after He, 2005)；(c) example of Raman spectrum of carbon material (modified after Wu et al., 2014)；(d) generate process of carbon material Raman characteristic peak (after Wu et al., 2014)

圖2d展示了碳質(zhì)物典型拉曼特征峰產(chǎn)生的過程，G峰是與布里淵區(qū)中心雙重簡并的iTO和iLO光學(xué)聲子相互作用的結(jié)果(Tuinstra and Koenig, 1970)。Tuinstra and Koenig(1970)最初估計G峰的位置約為1575cm-1，后人的研究也表明G峰的位置及強(qiáng)度隨入射波長的變化而變化(Wangetal., 1990；Matthewsetal., 1999；Satoetal., 2006)。D1峰最初歸因于碳質(zhì)物小晶體或大晶體邊緣的A1g振動模式(Tuinstra and Koenig, 1970)，但進(jìn)一步的研究發(fā)現(xiàn)D1峰是一個iTO聲子與一個缺陷谷間散射的結(jié)果(Beny-Bassez and Rouzaud, 1985；Malardetal., 2009；吳娟霞等，2014)。D2峰則為谷內(nèi)雙共振過程，兩次散射過程分別為與缺陷的谷內(nèi)散射和與K點附近的iLO聲子的非彈性谷內(nèi)散射過程(吳娟霞等，2014)，常出現(xiàn)在G峰段肩部，在結(jié)晶度低的碳質(zhì)物中，常不能將二者分開，只在1600cm-1附近現(xiàn)一個寬帶(Kouketsuetal.(2014)稱其為Gl帶)。D3峰推測來自無定形碳(Sadezkyetal., 2005)，Beny-Bassez and Rouzaud(1985)將該帶歸因于四面體碳等面外缺陷，認(rèn)為這種缺陷在石墨化過程早期出現(xiàn)。D4峰由sp2-sp3鍵或多烯類結(jié)構(gòu)的C-C和C=C伸縮振動產(chǎn)生。S1峰是與K點附近的iTO光學(xué)聲子發(fā)生兩次谷間非彈性散射產(chǎn)生的，S1峰拉曼位移約為D1峰的兩倍，因此通常表示為2D峰，但S1峰的產(chǎn)生與缺陷無關(guān)，并非D1峰的倍頻信號(吳娟霞等，2014)。在石墨化的最終階段S1峰將會分裂成兩部分，這種分裂被解釋為碳質(zhì)物二維結(jié)構(gòu)朝向三維結(jié)構(gòu)過渡的結(jié)果(Lespadeetal., 1984)。至于S2峰，它的含義還不是很清楚，通常只出現(xiàn)在結(jié)晶度低的碳質(zhì)物中，Tsuetal.(1978)將其歸因于C-H鍵的存在，但是這種假設(shè)仍存在爭議。

如上所述，產(chǎn)生D峰的缺陷種類很多，如結(jié)構(gòu)缺陷、面內(nèi)雜原子等。最初Tuinstra and Koenig(1970)根據(jù)D峰強(qiáng)度與晶體尺寸的關(guān)系，并給出著名的TK(Tuinstra Koenig)公式：

ID/IG=C(vL)/La

(i)

該公式首次建立了碳質(zhì)物石墨化程度與拉曼光譜的關(guān)系。式中ID/IG是D峰和G峰的強(qiáng)度比，La是納米晶體的橫向尺寸(平行于石墨烯平面)?？紤]到D峰可由其他類型的缺陷所激活(如局部雜質(zhì))，新近研究量化了ID對離子轟擊在石墨上引起的缺陷密度的依賴性(Ferreiraetal., 2010；Luccheseetal., 2010)。在這種情況下，TK關(guān)系變成了(ID/IG)∝(1/La2)，其La2是石墨烯平面中缺陷的平均距離。已知隨著石墨化程度(結(jié)晶度)的增加，La逐漸增大，因此可以通過ID/IG來量化表示碳質(zhì)物的石墨化程度(結(jié)晶度)(Wopenka and Pasteris, 1993；Beyssacetal., 2002b)。

變質(zhì)沉積巖中碳質(zhì)物拉曼光譜表現(xiàn)出隨變質(zhì)等級系統(tǒng)變化(圖3)，通常變質(zhì)等級越高，拉曼光譜譜峰值的數(shù)量越少(Pasteris and Wopenka, 1991；Wopenka and Pasteris, 1993；Yuietal., 1996；Beyssacetal., 2002a)，并且實驗和天然樣品觀測表明，碳質(zhì)物質(zhì)的石墨化過程中溫度是主要控制因素，而壓力為次級影響因素(圖4a, b)(Landis, 1971；Wopenka and Pasteris, 1993；Wadaetal., 1994；Beyssacetal., 2002a, b；Rahletal., 2005；Aoyaetal., 2010；Lahfidetal., 2010；Kouketsuetal., 2014)。對天然地質(zhì)樣品觀測顯示，隨變質(zhì)程度的增加，La(石墨結(jié)晶度參數(shù))隨之增加，即石墨化程度增加(圖4c)(Wopenka and Pasteris, 1993)。在此基礎(chǔ)上Beyssacetal.(2002a)根據(jù)拉曼參數(shù)R2(R2=D1/(G+D1+D2)area)與變質(zhì)溫度之間的正相關(guān)性，擬合出碳質(zhì)物拉曼參數(shù)與變質(zhì)溫度之間的關(guān)系方程，建立了碳質(zhì)物溫度公式(如圖4d)。

1.2 碳質(zhì)物拉曼光譜溫度計公式

Beyssacetal.(2002a)通過激光拉曼(514.5nm入射激光)分析了全球10個區(qū)域變質(zhì)地區(qū)的54個樣品，發(fā)現(xiàn)碳質(zhì)物結(jié)晶度與變質(zhì)峰期溫度密切相關(guān)，而與變質(zhì)壓力無關(guān)?；谧冑|(zhì)礦物組合獲得的變質(zhì)溫度和R2參數(shù)之間的線性關(guān)系，構(gòu)建了溫度計算公式(公式1，表1)。

表1變質(zhì)巖中碳質(zhì)物拉曼溫度計公式

Table 1 Various formulas for carbon material Raman thermometers

序號碳質(zhì)物拉曼光譜溫度計算公式應(yīng)用條件激發(fā)波長(nm)可計算溫度范圍(℃)參考文獻(xiàn)1T(℃)=-445(R2)+641區(qū)域變質(zhì)巖514.5300～650Beyssac et al.(2002a)2T(℃)=737.3+320.9R1-1067R2-80.638R12區(qū)域變質(zhì)巖532100～700Rahl et al.(2005)3T(℃)=1250RA1-469.75低級變質(zhì)巖532200～320Lahfid et al.(2010)4T(℃)=222.2RA2-60低級變質(zhì)巖532200～320Lahfid et al.(2010)56T(℃)=-2.15FWHMD1+478T(℃)=-6.78FWHMD2+535區(qū)域變質(zhì)巖R1大于1.0532532200～400150～200Kouketsu et al.(2014)7T(℃)=221R22-637.1R2+672.3接觸變質(zhì)巖532340～655Aoya et al.(2010)8T(℃)=91.4R22-556.3R2+676.3區(qū)域變質(zhì)巖514.5300～650Aoya et al.(2010)

注：R1=D1/Gintensity；R2=D1/(G+D1+D2)area；RA1=(D1+D4)/(D1+D2+D3+D4+G)area；RA2=(D1+D4)/(D1+D2+D3+D4+G)area；FWHMD1和FWHMD2分別是D1峰和D2峰的全峰半高寬

圖3 變質(zhì)巖中碳質(zhì)物拉曼光譜隨溫度增加演化圖數(shù)據(jù)來源: Beyssac et al., 2002b；Lahfid et al., 2010Fig.3 The evolution of Raman spectra of carbon material with increasing metamorphic temperature in metamorphic rockAll the samples come from Beyssac et al., 2002b；Lahfid et al., 2010

該溫度計適用于區(qū)域變質(zhì)巖，基于碳質(zhì)物石墨化過程不可逆性(Landis, 1971；Pasteris and Wopenka, 1991；Wopenka and Pasteris, 1993；Wadaetal., 1994；Buseck and Beyssac, 2014)，Beyssacetal.(2002a)解釋該溫度代表變質(zhì)峰期溫度。

為了將碳質(zhì)物拉曼溫度計外延至更低的溫度，Rahletal.(2005)通過激光拉曼(532nm激光源)對來自華盛頓州奧林匹克山脈的11個經(jīng)歷變質(zhì)的硅質(zhì)碎屑巖樣品進(jìn)行分析，拉曼數(shù)據(jù)結(jié)果顯示，R1(R1=D1/Gintensity)在低溫變質(zhì)過程中(100～350℃)仍然保持變化，從115℃的0.5增加到250℃的2.1，隨后在更高的溫度條件下減小，這些學(xué)者們結(jié)合Beyssacetal.(2002a)的數(shù)據(jù)，提出新的經(jīng)驗性溫度估算方程，將碳質(zhì)物拉曼溫度計的應(yīng)用范圍擴(kuò)展至100～700℃(公式2，表1)。

Rahletal.(2005)對該公式進(jìn)行了方差和殘差的檢驗，得出該公式在95%的置信度的情況下誤差為±50℃。但其也強(qiáng)調(diào)該公式僅僅是對已有數(shù)據(jù)擬合得到的一個經(jīng)驗公式，還需要在實際的地質(zhì)運用中進(jìn)行修正和完善。

圖4 碳質(zhì)物石墨化程度與溫度、壓力之間的關(guān)系(a)碳質(zhì)物石墨化程度與變質(zhì)溫度之間的關(guān)系(據(jù)Wada et al., 1994)；(b)溫度和壓力對碳質(zhì)物石墨化程度的影響(據(jù)Landis, 1971)；(c)變質(zhì)等級與碳質(zhì)物石墨化程度之間的關(guān)系(據(jù)Wopenka and Pasteris, 1993)；(d)碳質(zhì)物拉曼光譜參數(shù)(R2)與變質(zhì)溫度之間的關(guān)系(據(jù)Beyssac et al., 2002a)Fig.4 The relationships between the degrees of graphitization of carbon material and temperatures/pressures(a)relationship between the degrees of graphitization and the metamorphic temperatures(after Wada et al., 1994)；(b)effect of temperatures and pressures on the degrees of graphitization(after Landis, 1971)；(c)the relationship between metamorphic grades and degrees of graphitization of carbon material(after Wopenka and Pasteris, 1993)；(d)the relationship between the Raman parameter(R2)and the peak metamorphic temperature(after Beyssac et al., 2002a)

Lahfidetal.(2010)對格拉魯斯阿爾卑斯山低級變質(zhì)巖(鈣質(zhì)片巖、石英片巖)的研究發(fā)現(xiàn)，低級變質(zhì)巖中的碳質(zhì)物拉曼光譜與高級變質(zhì)巖中碳質(zhì)物拉曼光譜明顯不同且更加復(fù)雜，最顯著的區(qū)別是高等級變質(zhì)巖碳質(zhì)物拉曼光譜中不存在D3、D4峰。低級變質(zhì)巖中碳質(zhì)物拉曼光譜的復(fù)雜性也暗示Beyssacetal.(2002a)基于中級變質(zhì)程度提出的溫度計不能簡單的擴(kuò)展至低溫環(huán)境。進(jìn)而，Lahfidetal.(2010)通過比較各種參數(shù)(如R1、R2)描述碳質(zhì)物拉曼光譜隨變質(zhì)等級增加演變的能力，提出了兩個新的參數(shù)RA1(RA1=(D1+D4)/(D1+D2+D3+D4+G)area)和RA2(RA2=(D1+D4)/(D1+D2+D3+D4+G)area)，用于定量化描述格拉魯斯阿爾卑斯山低級變質(zhì)巖(鈣質(zhì)片巖、石英片巖)在200～320℃范圍內(nèi)的碳質(zhì)物拉曼光譜隨溫度演化(圖5)，通過簡單變換可得到應(yīng)用于200～320℃范圍內(nèi)的溫度計方程式(公式3、4，表1)。該溫度計誤差為±25℃，然而該溫度計僅由一個變質(zhì)區(qū)域的樣品構(gòu)建，正如Lahfidetal.(2010)所闡述的，用于校準(zhǔn)該溫度計的樣品變質(zhì)溫度在一個有限的范圍內(nèi)，因此，無法將該溫度計外延到溫度高于320℃的變質(zhì)環(huán)境。

近年來，Kouketsuetal.(2014)利用激光拉曼(532nm入射激光)對日本中央構(gòu)造帶內(nèi)的19個變質(zhì)沉積巖中碳質(zhì)物進(jìn)行分析，將獲得拉曼光譜參數(shù)(R1、R2等)與樣品已知的變質(zhì)溫度進(jìn)行比較，發(fā)現(xiàn)變質(zhì)溫度與D1峰和D2峰的全峰半高寬(FWHMD1、FWHMD2)呈負(fù)相關(guān)，并利用二者之間的相關(guān)性，構(gòu)建了適用于150～400℃溫度范圍的地質(zhì)溫度計，即半寬高(FWHMD1、FWHMD2)與變質(zhì)溫度的關(guān)系方程(公式5、6，表1)。

Kouketsuetal.(2014)認(rèn)為，在200～400℃的變質(zhì)條件下，公式5比公式6更為準(zhǔn)確，而在150～200℃的變質(zhì)條件下，公式5更準(zhǔn)確，公式5的誤差為±30℃，而公式6的誤差為±50℃(在FWHM柱狀圖不出現(xiàn)寬峰或者雙峰的條件下)。當(dāng)變質(zhì)溫度接近或者超過400℃時使用該公式需要謹(jǐn)慎，Kouketsuetal.(2014)提出了在此情況下使用該溫度計的條件，當(dāng)D1/G(強(qiáng)度)大于1.0時，公式5、6是適用的，當(dāng)D1/G(強(qiáng)度)小于1.0時，使用公式5或6時應(yīng)當(dāng)與前人提出的溫度計進(jìn)行相互檢驗和驗證。值得注意的是，前人提出的溫度計均是以比率作為參數(shù)，這樣做的好處是可以減少對實驗設(shè)備的依賴性，但Kouketsuetal.(2014)通過比較不同實驗室的測量結(jié)果，認(rèn)為FWHM測定不受實驗室變化的影響。此外由于該溫度計使用的溫度范圍為150～400℃，貫穿了碳質(zhì)物從無定形碳向結(jié)晶石墨轉(zhuǎn)化的整個過程，因此在選擇拉曼譜峰擬合方法時尤需嚴(yán)謹(jǐn)以便獲得可靠的溫度，拉曼譜峰分解方法可參考Kouketsuetal.(2014)提出的譜峰擬合程序。

Aoyaetal.(2010)將溫度計適用范圍從區(qū)域變質(zhì)巖擴(kuò)展至接觸變質(zhì)巖。接觸變質(zhì)巖與區(qū)域變質(zhì)巖最大的不同是巖石變形弱且加熱時間長，這種屬性導(dǎo)致接觸變質(zhì)巖中碳質(zhì)物結(jié)晶程度高，表現(xiàn)為更低的R2值(Aoyaetal., 2010)，由于使用532nm激光比使用514.5nm激光能夠獲得相對更大的R2值(Wangetal., 1990；Matthewsetal., 1999；Satoetal., 2006；Beyssac and Lazzeri, 2012)，因此Aoyaetal.(2010)使用532nm激光拉曼對來自日本Daimonji和Kasuga地區(qū)的領(lǐng)家變質(zhì)帶中的11個接觸變質(zhì)巖樣品(原巖為泥質(zhì)巖和砂巖)進(jìn)行分析，并提出適用于接觸變質(zhì)作用過程的溫度估算的公式(公式7，表1)。

表2碳質(zhì)物溫度計的影響因素

Table 2 Affecting factors of the carbonaceous Raman spectroscopy thermometers

碳質(zhì)物拉曼光譜溫度計影響因素誤差來源誤差減小措施參考文獻(xiàn)光譜擬合引入的誤差碳質(zhì)物自身特性引入的誤差拉曼光譜系統(tǒng)使用時的誤差基線校正、譜峰擬合函數(shù)的選擇樣品制備、樣品的非均質(zhì)性及碳質(zhì)物結(jié)構(gòu)各向異性入射激光波長及拉曼光譜獲取時間基線校準(zhǔn)中明確所有假設(shè), 使用洛倫茲函數(shù)擬合譜圖測量薄片內(nèi)部碳質(zhì)物,單個樣品測量數(shù)大于50根據(jù)計算公式合理選擇激光波長,墊塊銅片導(dǎo)熱或在 薄片中測量Lünsdorf et al.(2014)

圖5 格拉魯斯阿爾卑斯山地區(qū)低級變質(zhì)巖內(nèi)碳質(zhì)物拉曼光譜參數(shù)隨變質(zhì)溫度演化圖(據(jù)Lahfid et al., 2010)(a)RA1隨變質(zhì)溫度演化圖；(b)RA2隨變質(zhì)溫度演化圖Fig.5 Evolution of carbon material Raman spectral parameters with metamorphic temperatures in low-grade metamorphic rocks in the Glarus Alps(after Lahfid et al., 2010)(a)evolution of the RA1 with temperature；(b)evolution of the RA2 with temperature

這些作者強(qiáng)調(diào)該溫度計適用范圍是340～655℃的接觸變質(zhì)巖，測試過程使用的入射激光波長為532nm。此外，Aoyaetal.(2010)發(fā)現(xiàn)Beyssacetal.(2002a)所提出的溫度計在溫度小于330℃及大于550℃時估算存在明顯偏差。因此利用Beyssacetal.(2002a)報道的數(shù)據(jù)對區(qū)域變質(zhì)巖的溫度計進(jìn)行了進(jìn)一步修正(入射激光514.5nm)，如公式8(表1)。

修正后的溫度計算公式，適用溫度范圍為300～650℃(Aoyaetal., 2010)，Aoyaetal.(2010)認(rèn)為Beyssacetal.(2002a)提出的溫度計不僅適用于區(qū)域變質(zhì)巖，也適用于接觸變質(zhì)巖，只不過在接觸變質(zhì)巖溫度估算中誤差較大(±50℃)，但該溫度計(公式7，表1)可通過增加測試點數(shù)量來降低溫度估算誤差(主張單個樣品測量的數(shù)量應(yīng)大于50個)，誤差可降至±30℃。

1.3 碳質(zhì)物拉曼光譜溫度計的影響因素

碳質(zhì)物拉曼光譜測試從樣品制備到光譜處理過程中產(chǎn)生誤差的原因可歸結(jié)為三大類(表2)：(1)光譜的擬合引入的誤差；(2)碳質(zhì)物自身特性引入的誤差；(3)拉曼光譜系統(tǒng)使用時的誤差(Lünsdorfetal., 2014)。光譜擬合過程中的誤差來自于不同的基線校正，用于峰匹配的不同數(shù)學(xué)函數(shù)以及不同的拉曼光譜模型。碳質(zhì)物自身特性引入的誤差主要來自碳質(zhì)物結(jié)構(gòu)的各向異性、樣品的制備及樣品的非均質(zhì)性。拉曼光譜系統(tǒng)使用產(chǎn)生的誤差主要與入射激光的波長及測量時間有關(guān)，下文將對這三類影響因素做進(jìn)一步的闡釋。

圖6 碳質(zhì)物拉曼光譜擬合實例(使用Voigt函數(shù)三波段擬合)Fig.6 Example of peak-fitting procedures for Raman spectroscopy of carbon material(three-bands fitting by using Voigt function)

圖7 薄片表面碳質(zhì)物與內(nèi)部碳質(zhì)物拉曼光譜對比(a)青藏高原嘉黎韌性走滑剪切帶內(nèi)的糜棱巖化片巖單偏光顯微照片；(b)碳質(zhì)物反射光及透射光下的顯微照片；(c)表面及內(nèi)部碳質(zhì)物拉曼光譜Fig.7 Comparison of Raman spectra of carbon material obtained at the surface of a thin section(triangle)and at inside within the thin section(circle)(a)photomicrograph under single polarized light for the mylonite schist from Jiali ductile strike-slip shear zone of the Tibet；(b)photomicrographs under reflected and transmitted light for carbon material from Fig.7a；(c)Raman spectroscopy of carbon material obtained at the surface of a thin section and at inside within the thin section

從拉曼光譜中獲取碳質(zhì)物結(jié)晶度定量參數(shù)需要進(jìn)行光譜擬合，光譜擬合過程中用于峰匹配的數(shù)學(xué)函數(shù)以及不同的拉曼光譜模型的選擇受個人影響較大，建立一種可靠且重復(fù)的碳質(zhì)物拉曼光譜擬合的方法是開發(fā)適用于各種變質(zhì)溫度的碳質(zhì)物拉曼溫度計需要克服的最大問題之一(Kouketsuetal., 2014)。目前，光譜擬合普遍使用Peakfiting和Fityk軟件，通常不同的光譜擬合軟件對擬合結(jié)果不會產(chǎn)生影響(Lünsdorfetal., 2014)。光譜擬合第一步是基線校正，這對任何擬合程序的質(zhì)量和可重復(fù)性都至關(guān)重要，對于石墨很簡單，因為基線是線性的，如果結(jié)構(gòu)很混亂，基線會變得很復(fù)雜，通常伴隨著熒光性組分的影響，熒光可能是由相鄰的礦物質(zhì)導(dǎo)致的，基線校準(zhǔn)過程中必須明確規(guī)定做出的所有假設(shè)。關(guān)于光譜擬合，使用最多的函數(shù)是Voigt(洛倫茲和高斯函數(shù)的組合)和洛倫茲函數(shù)。對于高結(jié)晶度的石墨，使用Voigt函數(shù)的三波段擬合(G、D1和D2)將產(chǎn)生優(yōu)異的結(jié)果(圖6；Beyssacetal., 2003, 2004)；隨著碳質(zhì)物結(jié)晶度的降低，在拉曼光譜一階序區(qū)出現(xiàn)D3峰，對于低結(jié)晶度碳質(zhì)物光譜，使用五個Voigt函數(shù)將會增加系統(tǒng)的自由度并導(dǎo)致多種擬合方案。另一種方法是使用純洛倫茲或高斯函數(shù)來降低系統(tǒng)的自由度，以此獲得唯一且可重復(fù)的擬合方案(Sadezkyetal., 2005；Lahfidetal., 2010)。Lünsdorfetal.(2014)認(rèn)為對于低結(jié)晶度碳質(zhì)物拉曼光譜只使用洛倫茲函數(shù)擬合將得到更好的結(jié)果，而隨著碳質(zhì)物結(jié)晶度的增加，擬合函數(shù)的選擇不會影響擬合結(jié)果。由于洛倫茲函數(shù)擬合對低結(jié)晶樣品和高結(jié)晶樣品光譜擬合均表現(xiàn)良好，因此Lünsdorfetal.(2014)建議使用該函數(shù)擬合所有結(jié)晶度的樣品。

碳質(zhì)物自身固有的誤差主要來源于樣品的制備、樣品的非均質(zhì)性及碳質(zhì)物結(jié)構(gòu)的各向異性(Beyssac and Lazzeri, 2012；Lünsdorfetal., 2014)。Pasteris(1989)研究發(fā)現(xiàn)，除碳質(zhì)物固有結(jié)構(gòu)缺陷外，鋸切或拋光可能產(chǎn)生新的結(jié)構(gòu)缺陷，對薄片表面碳質(zhì)物進(jìn)行直接測試將會受新產(chǎn)生的結(jié)構(gòu)缺陷影響。這些缺陷在拉曼光譜中無法區(qū)分，且隨著碳質(zhì)物的結(jié)晶度增加鋸切或拋光產(chǎn)生的影響也增加，這種效應(yīng)對于晶格完整的石墨最為顯著(Beyssacetal., 2003)。Pasteris(1989)認(rèn)為通過測量薄片內(nèi)部的碳質(zhì)物可避免這種誤差，即將激光聚焦在透明/半透明相鄰礦物顆粒下方的碳質(zhì)物來實現(xiàn)。筆者的研究也從“反面”進(jìn)一步證實了這種效應(yīng)，如圖7所示，圖中表層碳質(zhì)物其實并未真正的暴露在薄片上表面，相反，內(nèi)部的碳質(zhì)物已經(jīng)暴露在薄片的下表面，受拋光的影響產(chǎn)生缺陷峰。因此將激光聚焦在透明/半透明相鄰礦物顆粒下方的碳質(zhì)物時，也應(yīng)注意透明顆粒下方的碳質(zhì)物是否暴露在巖石薄片的下表面。

除了薄片制備過程會產(chǎn)生誤差外，碳質(zhì)物取向(結(jié)構(gòu)各向異性)相對于入射激光的方向，對拉曼光譜也有影響，尤其是結(jié)晶程度高的石墨(Katagirietal., 1988；Wangetal., 1989；Compagninietal., 1997)。通常入射激光垂直于C軸比平行C軸產(chǎn)生的G峰更小，D1峰更大，即R2值在垂直于C軸的測量中較大(Aoyaetal., 2010；Beyssac and Lazzeri, 2012)。Aoyaetal.(2010)比較了入射激光平行和垂直于C軸的測量結(jié)果，以檢驗取向的影響，結(jié)果表明對于結(jié)晶度高的碳質(zhì)物，平均R2比率不受取向的影響，但R2非零值，在垂直于C軸和平行于C軸方向上存在很大的差異，入射激光垂直于C軸會比平行于C軸產(chǎn)生更大的R2值，但只要測量的數(shù)量足夠多(N≥25)，即使是結(jié)晶度高的碳質(zhì)物，R2的平均值也不受碳質(zhì)物取向的影響。上述的討論表明，結(jié)晶度高的碳質(zhì)物R2取向效應(yīng)更為顯著，這是由于結(jié)晶度高的石墨具有較強(qiáng)的各向異性(Beyssacetal., 2002a)。

相比碳質(zhì)物的非均質(zhì)性，碳質(zhì)物取向的影響可以忽略(Aoyaetal., 2010)。天然碳質(zhì)物結(jié)構(gòu)和化學(xué)成分在微米尺度上通常是不均勻的(Planeixetal., 1994；Dreyer and Bielawski, 2011)。很多學(xué)者通過透射電鏡(TEM)和X射線吸收近邊光譜(XANES)在納米尺度上對碳質(zhì)物的非均質(zhì)性進(jìn)行了觀測(Buseck and Huang, 1985；Beyssacetal., 2002b；Bernardetal., 2007, 2008；van Zuilenetal., 2012)，結(jié)果表明碳質(zhì)物結(jié)構(gòu)非均質(zhì)性非常強(qiáng)，結(jié)晶度高的石墨和結(jié)晶度低的無定型碳共存于同一塊巖石中。如隕石中Fe-Ni顆粒中包含石墨片，而隕石基質(zhì)中則存在低結(jié)晶度無定型碳(Pianietal., 2012)，以及變質(zhì)沉積巖中由有機(jī)物石墨化產(chǎn)生的碎屑石墨和無定型碳共存(Beyssacetal., 2004)。在這兩種情況下，石墨和無定型碳不是同時形成的，但卻經(jīng)歷相同的石墨化過程。因此，在拉曼研究過程中，對碳質(zhì)物進(jìn)行巖石學(xué)研究以描述其形態(tài)和與礦物的結(jié)構(gòu)共生關(guān)系至關(guān)重要。Aoyaetal.(2010)通過實驗表明這種來自于碳質(zhì)物非均質(zhì)性的誤差可以通過大量的測試數(shù)據(jù)(單個樣品測量數(shù)N>50)來減小。

拉曼系統(tǒng)固有誤差可能來自于拉曼光譜的獲取時間及入射波長的大小。由于黑色的碳質(zhì)物會吸收可見光范圍內(nèi)的光(激光)，不同的拉曼光譜獲取時間即不同的加熱時間，碳質(zhì)物受熱會導(dǎo)致光氧化或局部燃燒，從而改變局部化學(xué)成分和結(jié)構(gòu)(Beyssac and Lazzeri, 2012)。激光誘導(dǎo)的加熱是G峰、D峰左移的原因(Everalletal., 1991；Kagietal., 1994)。Tuinstra and Koenig(1970)測得G峰位置為1575cm-1可能是由于這種原因造成的(現(xiàn)在普遍接受的是1580cm-1)。對于從巖石中提取出的碳質(zhì)物，激光加熱效應(yīng)尤為顯著，解決方法是在提取出的碳質(zhì)物下方墊塊銅片并置于氬氣環(huán)境中(Beyssac and Lazzeri, 2012)。對于薄片樣品，激光加熱引起的誤差可以忽略(Aoyaetal., 2010)。這是由于薄片中碳質(zhì)物不會直接暴露在空氣中，并且誘導(dǎo)產(chǎn)生的熱量也可以通過上方或周圍的礦物傳導(dǎo)散去。

此外，入射激光的能量對碳質(zhì)物拉曼光譜具有顯著的影響，尤其是缺陷峰(D峰)，通常入射激光能量的減小會使D峰向左偏移(Satoetal., 2006)。碳質(zhì)物拉曼光譜常用的入射激光能量為2.33eV(532nm)，地球科學(xué)領(lǐng)域常用的激光能量為2.41eV(514nm)，但由于在低等級(低溫)變質(zhì)巖中，更高的激光能量能夠減少熒光背景并獲得較大的R1、R2比值，因此近年來所提出的溫度計多使用2.33eV(532nm)的入射激光來獲取拉曼圖像(Aoyaetal., 2010；Kouketsuetal., 2014)。

2 碳質(zhì)物拉曼光譜溫度計與傳統(tǒng)溫度計對比

圖8 碳質(zhì)物溫度計與傳統(tǒng)溫度計估算結(jié)果對比Fig.8 Comparison of calculated results from the carbon material thermometer and traditional thermometers

上述不同學(xué)者所提出的溫度計已被應(yīng)用于不同的變質(zhì)巖區(qū)，特別是Beyssacetal.(2002a)提出的基于R2構(gòu)建的溫度計應(yīng)用最為廣泛。為評估碳質(zhì)物拉曼溫度計可靠性，本次研究對已發(fā)表的同一地區(qū)同一地體碳質(zhì)物拉曼溫度計所估算的溫度與傳統(tǒng)溫度計計算的溫度進(jìn)行了對比，包括尼泊爾低喜馬拉雅泥質(zhì)片巖(Beyssacetal., 2004)，尼泊爾DZAKAACHU地區(qū)的藏南拆離系中的糜棱巖化大理巖(Cottleetal., 2011)，印度阿魯納恰爾地區(qū)喜馬拉雅造山帶中的長英質(zhì)片巖、云母片巖、千枚巖(Mathewetal., 2013)，龍門山片巖(Robertetal., 2010)，西阿爾卑斯Oisans地塊的灰?guī)r、頁巖及綠泥石片巖(Bellangeretal., 2015)，新西蘭奧塔哥片巖(Huetal., 2015)，西摩洛哥Jebilet地塊綠片巖、角閃巖相變質(zhì)的泥質(zhì)巖(Delchinietal., 2016)，日本西紀(jì)伊半島泥質(zhì)片巖(Yoshidaetal., 2016)，希臘基克拉澤斯藍(lán)片巖和榴輝巖(Laurentetal., 2018)(表3)。

對比結(jié)果如圖8所示，圖中橫坐標(biāo)為碳質(zhì)物拉曼溫度計(基于Beyssacetal.(2002a)和Rahletal.(2005)溫度計公式)估算的峰期變質(zhì)溫度，縱坐標(biāo)為傳統(tǒng)溫度計(相平衡模擬、礦物組合)獲取的變質(zhì)溫度，虛線為y=x直線。結(jié)果顯示碳質(zhì)物拉曼溫度計與傳統(tǒng)溫度計所獲得的變質(zhì)溫度值均分布于y=x斜直線附近，總體上碳質(zhì)物拉曼溫度計與傳統(tǒng)溫度計估算的溫度較為接近，二者之間的差別也在分析和校準(zhǔn)誤差范圍內(nèi)。拉曼溫度計相對誤差通常小于15℃，遠(yuǎn)小于傳統(tǒng)溫度計(±50℃)。但Robertetal.(2010)和Cottleetal.(2011)在龍門山、尼泊爾DZAKAA CHU地區(qū)的研究表明，碳質(zhì)物拉曼溫度計與傳統(tǒng)溫度計存在較大差異，二者結(jié)果相差約100℃，原因是Cottleetal.(2011)和Robertetal.(2010)用傳統(tǒng)溫度計估算的溫度代表退變質(zhì)溫度，而碳質(zhì)物拉曼溫度計結(jié)果為峰期變質(zhì)溫度，且不受后期退變質(zhì)的影響(Landis, 1971；Pasteris and Wopenka, 1991；Wopenka and Pasteris, 1993；Wadaetal., 1994；Beyssacetal., 2002a, b, 2004, 2007；Rahletal., 2005；Negroetal., 2006；Angiboustetal., 2009；Gabaldaetal., 2009；Lahfidetal., 2010)，因此，二者之間存在較大差異。

3 碳質(zhì)物拉曼光譜溫度計在地學(xué)中的應(yīng)用

解析造山帶形成過程的空間和時間熱演化是造山帶研究的核心。變質(zhì)溫度估算、熱年代學(xué)和地質(zhì)年代學(xué)已成為約束造山帶熱演化的關(guān)鍵信息。傳統(tǒng)溫度計對于特定礦物組合的要求限制了其在造山帶中的應(yīng)用，而每個造山帶均會產(chǎn)生增生楔沉積物質(zhì)，這些沉積單元中富含碳質(zhì)物，為碳質(zhì)物拉曼溫度計的使用提供了基本的物質(zhì)基礎(chǔ)。 碳質(zhì)物拉曼光譜溫度計一經(jīng)提出便得到廣泛使用，如在造山帶熱結(jié)構(gòu)和熱演化史方面(Beyssacetal., 2004；Robertetal., 2010；Cottleetal., 2011；Kellett and Grujic, 2012；Cooperetal., 2013； Mathewetal., 2013；Scharfetal., 2013；Bellangeretal., 2015；Boutouxetal., 2016；Delchinietal., 2016；Yoshidaetal., 2016；Laurentetal., 2018；Mollietal., 2018)、成礦流體的古流體熱異常(Clercetal., 2015；Delchinietal., 2017)、以及強(qiáng)應(yīng)變帶剪切熱過程等(Furuichietal., 2015；Morietal., 2015；Itoetal., 2017)。在實際運用的過程中，其有效性也得到了進(jìn)一步的檢驗(Beyssacetal., 2004；Scharfetal., 2013；Delchinietal., 2016；Laurentetal., 2018)。

表3碳質(zhì)物溫度計與傳統(tǒng)溫度計對比表

Table 3 Calculated results from carbon material geothermometers and conventional geothermometers

地區(qū)碳物質(zhì)拉曼溫度計傳統(tǒng)溫度計溫度(℃)標(biāo)準(zhǔn)誤差(℃)溫度(℃)系統(tǒng)誤差(℃)碳物質(zhì)拉曼溫度計公式傳統(tǒng)溫度計方法參考文獻(xiàn)尼泊爾低喜馬拉雅541.716.1558.820.4Beyssac et al.(2002a)石榴石-云母Fe2+-Mg交換溫度計(最高變質(zhì)溫度)Kohn et al. (2001)、Beyssac et al. (2004)尼泊爾DZAKAA CHU地區(qū)藏南拆離系338.0479.0493.0545.0533.0544.0530.04.06.05.05.06.04.07.0350.0380.0365.0430.0430.0450.0550.050.040.040.020.020.040.030.0Beyssac et al.(2002a)多平衡礦物組合溫壓計,石英、長石微觀結(jié)構(gòu)(退變質(zhì)溫度)Cottle et al. (2011)印度阿魯納恰爾喜馬拉雅造山帶659.0681.8650.2669.3598.2422.6607.5617.0603.2610.2564.2413.47.87.38.510.913.412.33.32.33.84.59.66.7685.0677.0594.0603.0552.0458.0685.0677.0594.0603.0552.0458.055.035.058.070.045.011.055.035.058.070.045.011.0Rahl et al.(2005)Beyssac et al.(2002a)石榴石-云母Fe2+-Mg交換溫度計(最高變質(zhì)溫度)Goswami et al.(2009)、 Mathew et al. (2013)龍門山599.0562.0515.0504.0547.010.010.010.010.010.0490.0480.0470.0417.0428.040.070.056.025.036.0Beyssac et al.(2002a)石榴石-云母Fe2+-Mg交換溫度計、多平衡礦物組合溫壓計(退變質(zhì)溫度)Huang et al.(2003a,b)、de Sigoyer et al.(2008)、 Robert et al. (2010)西阿爾卑斯Oisans地塊306.0345.0335.0308.0306.0345.0335.0308.013.09.05.05.013.09.05.05.0294.0296.0297.0288.0252.0313.0294.0269.017.09.012.013.061.055.045.051.0Beyssac et al.(2002a)石榴石-云母Fe2+-Mg交換溫度計、綠泥石溫度計(最高變質(zhì)溫度)Bellanger et al. (2015)新西蘭奧塔哥239.0304.0226.0283.0433.0339.08.012.08.012.0200.0300.0200.0400.0400.0400.050.050.050.050.050.050.0Beyssac et al.(2002a)礦物組合(最高變質(zhì)溫度)Mortimer. (2000)、Hu et al. (2015)西摩洛哥Jebilet地塊618.0613.0616.014.012.011.0620.0620.0625.012.012.025.0Beyssac et al.(2002a)基于相圖計算的溫壓計(最高變質(zhì)溫度)Delchini et al. (2016)日本西紀(jì)伊半島285.0289.07.013.0275.0270.025.030.0Beyssac et al.(2002a)礦物組合(最高變質(zhì)溫度)Yoshida et al. (2016)、Banno and Sakai (1989)希臘基克拉澤斯536.124.5530.030.0Beyssac et al.(2002a)石榴石-云母Fe2+-Mg交換溫度計、基于相圖計算的溫壓計(最高變質(zhì)溫度)Laurent et al. (2018)

Beyssacetal.(2004)利用該溫度計限定了低喜馬拉雅熱演化，研究表明低喜馬拉雅(LH)經(jīng)歷了一次大規(guī)模的熱變質(zhì)作用，溫度從頂部540℃逐漸降低到最深折返層內(nèi)的330℃。Mathewetal.(2013)利用Beyssacetal.(2002a)和Rahletal.(2005)提出的溫度計估算了西阿魯納恰爾拉邦低喜馬拉雅和西瓦里克熱演化史，研究發(fā)現(xiàn)大喜馬拉雅序列(GHS)在主中央逆沖區(qū)溫度從700～800℃降至650～700℃，并在西瓦里克的中上新世序列中進(jìn)一步降低至200℃，據(jù)此預(yù)測GHS倒置的變質(zhì)梯度演化與渠道流模型緊密相關(guān)。在阿爾卑斯造山帶，Bellangeretal.(2015)使用碳質(zhì)物拉曼光譜溫度計、綠泥石溫度計結(jié)合白云母40Ar/39Ar原位測年對西阿爾卑斯Oisans地塊的演化進(jìn)行了溫度和年齡的限制，提出了西阿爾卑斯外部結(jié)晶塊(ECM)新的縮短模式。其后有學(xué)者結(jié)合鋯石的(U-Th-Sm)/He與碳質(zhì)物拉曼溫度計數(shù)據(jù)，約束了西阿爾卑斯山ECM的熱演化和折返過程(Boutouxetal., 2016)。Mollietal.(2018)將碳質(zhì)物拉曼光譜溫度計應(yīng)用于亞平寧山脈，其搜集的數(shù)據(jù)修正了亞平寧山脈熱結(jié)構(gòu)及阿普亞內(nèi)山核部變形與變質(zhì)結(jié)構(gòu)之間的關(guān)系，解釋認(rèn)為殘存的古熱結(jié)構(gòu)及部分折返與地殼增厚有關(guān)，提出阿普亞內(nèi)雜巖為冷的變質(zhì)核雜巖。Laurentetal.(2018)利用石榴石平衡模型、P-T相圖結(jié)合碳質(zhì)物拉曼光譜測溫數(shù)據(jù)恢復(fù)了希臘基克拉迪藍(lán)片巖(CBU，Greece)峰期變質(zhì)溫度，表明該地區(qū)經(jīng)歷了多階段的折返過程，并認(rèn)為這種折返由板塊回卷驅(qū)動制約。碳質(zhì)物拉曼光譜溫度計還被成功運用在洋殼深俯沖型低溫-超高壓變質(zhì)雜巖帶。通過碳質(zhì)物拉曼光譜溫度計和石英包裹體彈性壓力計的聯(lián)合應(yīng)用，Bayetetal.(2018)對中亞造山帶西端的西南天山低溫-超高壓變質(zhì)雜巖帶進(jìn)行了區(qū)域大尺度的P-T繪圖，發(fā)現(xiàn)前人所厘定的HP/UHP混雜巖的北部整體上經(jīng)歷了較為均一的低溫超高壓榴輝巖相峰期變質(zhì)條件(2.5±0.2GPa，536±11℃)。此外，在對區(qū)域內(nèi)前人工作進(jìn)行了系統(tǒng)總結(jié)的前提下，Tanetal.(2019)基于激光微區(qū)40Ar/39Ar測年法、鋯石U-Pb定年、相平衡P-T計算及碳質(zhì)物拉曼光譜溫度計等手段，開展了區(qū)域大尺度P-T-t-D繪圖，指出前人所厘定的高壓/超高壓混雜帶從北到南應(yīng)劃分為具有不同峰期變質(zhì)P-T-t-D演化歷史的四個次級構(gòu)造單元，分別是：北部超高壓單元(～2.75GPa，480～560℃，320±1Ma)，中部藍(lán)片巖/榴輝巖相基性單元(～2.1GPa，505℃，332±2Ma)，中部低藍(lán)片巖相超基性-基性單元(～1.45GPa，485℃，359±2Ma)及南部低藍(lán)片巖-綠片巖相過渡帶(>0.7～1.0GPa，470～520℃，ca. 280～310Ma)，所揭示出的四期短暫的深部巖石折返事件暗示了俯沖地溫梯度存在時空上的漸進(jìn)變化(～12℃/km到6～7℃/km)，這可能反映了俯沖帶系統(tǒng)隨著時間發(fā)生冷卻。

許多礦床存在于富含碳質(zhì)物的變質(zhì)沉積物中。因此對碳質(zhì)物結(jié)晶度的研究提供了估算巖石變質(zhì)峰值溫度的精確工具。地層溫度通常由地層的埋藏剝蝕、巖漿侵入和流體循環(huán)控制。在熱液礦床中，碳質(zhì)物可以提供熱事件及之后最高溫的準(zhǔn)確信息，因此，碳質(zhì)物可以用于研究古溫度并闡明流體流動引起的熱異常(Clercetal., 2015；Delchinietal., 2017)。為了檢驗碳質(zhì)物拉曼光譜溫度計在熱液礦床環(huán)境中檢測熱異常的能力，以美國加利福尼亞州的弗朗西斯卡雜巖Lucia微地體為例進(jìn)行了研究，表明碳質(zhì)物拉曼光譜溫度計是確定由熱液流動引起的熱異常的可靠工具(Lahfidetal., 2010)。碳質(zhì)物拉曼溫度計在北比利牛斯Arguenos-Moncaup地區(qū)估算的變質(zhì)溫度也與前人研究結(jié)果吻合(Clercetal., 2015)，碳質(zhì)物拉曼光譜溫度數(shù)據(jù)集清楚地表明了地幔折返和熱異常之間的聯(lián)系，并且證實了該地區(qū)變質(zhì)作用與白堊紀(jì)地殼減薄之間的聯(lián)系。此外，碳質(zhì)物拉曼光譜溫度計提供了一個方便快捷獲取巖石古溫度的方法，并且能夠廣泛用于礦產(chǎn)勘查(Delchinietal., 2017)。Delchinietal.(2017)將碳質(zhì)物拉曼光譜溫度計運用至礦產(chǎn)勘探，通過碳質(zhì)物拉曼溫度計解析了Jebilet地塊礦床的熱歷史，表明該地區(qū)經(jīng)歷了構(gòu)造、變質(zhì)和熱液事件，這些因素共同影響著礦床的形成，同時也檢驗了碳質(zhì)物拉曼光譜溫度計對于檢測古熱異常的有效性。

圖9 MTL地區(qū)熱異常成因(據(jù)Mori et al., 2015修改)(a)垂直MTL斷層的碳質(zhì)物拉曼溫度剖面；(b)MTL地區(qū)構(gòu)造演化模型Fig.9 Causes of thermal anomalies in MTL areas(modified after Mori et al., 2015)(a)cross-sections across the MTL fault with carbon material Raman spectrum temperaturesand；(b)tectonic evolution model in MTL region

在大陸地殼發(fā)育長期活動的斷層體系，如走滑斷層，轉(zhuǎn)換斷層和俯沖邊界。對主斷層的剪切應(yīng)力的估算是了解地殼強(qiáng)度和大地震的產(chǎn)生以及與運動相關(guān)剪切生熱程度的重要前提。出露地表的斷裂帶記錄了深部地殼的變形，是觀測剪切生熱和斷裂強(qiáng)度的天然實驗室。由于缺乏合適的礦物以及溫度變化較小，基于傳統(tǒng)巖相學(xué)的溫度計不再適用(d’Alessioetal., 2003；Yamadaetal., 2007；Tagami, 2012)。Itoetal.(2017)通過分析天然和實驗合成的粘土巖中假熔巖的拉曼光譜，認(rèn)為假熔巖中的碳質(zhì)物拉曼光譜特征可以用于評估共震斷層強(qiáng)度。比較碳質(zhì)物質(zhì)結(jié)晶動力學(xué)實驗的結(jié)果與熱模擬結(jié)果表明，碳質(zhì)物拉曼光譜溫度計非常適合估算主斷層帶的剪切生熱(Morietal., 2015)，這些作者對日本中央構(gòu)造帶(MTL)泥質(zhì)片巖分析發(fā)現(xiàn)在垂直于MTL斷層面的150m附近存在約60℃的峰值溫度上升(圖9a)，這種與斷層空間相關(guān)聯(lián)熱異常意味它是由剪切加熱引起的(圖9b)。然而，近年來摩擦實驗研究表明，鏡質(zhì)體反射率(R0)和碳質(zhì)物拉曼光譜溫度計均可以用于檢測斷層摩擦生熱，但常規(guī)鏡質(zhì)體反射率(R0)和碳質(zhì)物拉曼光譜溫度計不能精準(zhǔn)估算地震斷層滑動期間的峰值溫度，需要結(jié)合高位移速率下快速加熱的反應(yīng)動力學(xué)及碳質(zhì)物原始微觀結(jié)構(gòu)和成分相關(guān)研究，建立適用于斷層摩擦加熱的碳質(zhì)物拉曼光譜溫度計(Furuichietal., 2015)。

4 結(jié)論

本文綜述了碳質(zhì)物拉曼溫度計的基本原理和影響因素，對比分析了碳質(zhì)物拉曼溫度計與傳統(tǒng)溫度計估算結(jié)果，并討論了該方法在造山帶熱結(jié)構(gòu)恢復(fù)、熱演化、以及斷層剪切生熱等方面的研究實例。

該方法是基于碳質(zhì)物石墨化程度(結(jié)晶度)與溫度之間的關(guān)系構(gòu)建的溫度計，其計算公式是通過分析大量天然與實驗樣品統(tǒng)計歸納得到的經(jīng)驗性公式，估算的溫度代表峰期變質(zhì)溫度，不受后期退變質(zhì)的影響，可估算的溫度范圍為100～700℃(相對誤差±15℃，絕對誤差±30℃)，與傳統(tǒng)溫度計相比更為便捷。

對于本文總結(jié)歸納的不同版本的碳質(zhì)物拉曼光譜溫度計，在使用的過程中我們需要根據(jù)研究區(qū)的地質(zhì)條件和環(huán)境給予選擇不同的溫度計公式，就目前碳質(zhì)物拉曼光譜溫度計使用情況來看，普遍使用的是Beyssacetal.(2002a)提出的溫度計及Rahletal.(2005)提出的溫度計，但需要特別注意的是，Beyssacetal.(2002a)和Rahletal.(2005)提出的溫度計適用于區(qū)域變質(zhì)巖，當(dāng)樣品來自于接觸變質(zhì)帶時，選擇Aoyaetal.(2010)提出的溫度計估算溫度可能會得到更加準(zhǔn)確的結(jié)果。

碳質(zhì)物拉曼光譜溫度計相比傳統(tǒng)溫度計其最大的特點是記錄最高溫變質(zhì)溫度，不受后期退變質(zhì)影響，因此結(jié)合傳統(tǒng)溫壓計可以更好的描述巖石多期次變質(zhì)歷史。但同時也需注意，影響碳質(zhì)物結(jié)晶度的因素不僅僅只有溫度，在低溫環(huán)境中，壓力、變質(zhì)作用持續(xù)時間、碳質(zhì)物前體類型甚至是變形均會影響碳質(zhì)物結(jié)晶度，但目前為止，尚未有人量化這些因素對碳質(zhì)物拉曼光譜溫度計的影響，這也是碳質(zhì)物拉曼光譜溫度計今后需要進(jìn)一步完善的地方。

致謝本研究成果還得到北京大學(xué)造山帶與地殼演化教育部重點實驗室SEM-EBSD實驗室和TIMA實驗室的資助；審稿專家及俞良軍老師嚴(yán)謹(jǐn)細(xì)致審閱全文，提出了很多建設(shè)性的修改意見與建議，對論文質(zhì)量的提高大有幫助；在此一并表示衷心感謝。