吳春明 劉嘉惠WU ChunMing and LIU JiaHui

中國(guó)科學(xué)院大學(xué)地球與行星科學(xué)學(xué)院，北京 100049

College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China

早期的野外觀察和隨后的一系列巖石學(xué)實(shí)驗(yàn)，證明壓力(P, pressure)、溫度(T, temperature)不僅決定了變質(zhì)巖、巖漿巖中各種“相”(礦物、流體、熔體)出現(xiàn)和消失的先后順序，還直接制約它們化學(xué)成分的變化。例如，花崗巖中共生的斜長(zhǎng)石和堿性長(zhǎng)石的成分變化情況，就明顯受到溫度條件的控制，據(jù)此人們建立了溫度與兩種長(zhǎng)石化學(xué)成分之間的函數(shù)關(guān)系，即二長(zhǎng)石溫度計(jì)(Barth, 1934; 轉(zhuǎn)引自Putirka, 2008)。二長(zhǎng)石溫度計(jì)是人類歷史上最早建立的地質(zhì)溫度計(jì)。另一方面，如果壓力占據(jù)了主導(dǎo)地位即壓力條件是控制礦物組合中各種相化學(xué)成分的主要因素，就可據(jù)此建立地質(zhì)壓力計(jì)。有些地質(zhì)溫度計(jì)和壓力計(jì)，建立在涉及流體的“模式反應(yīng)”(model reaction)基礎(chǔ)上，例如綠泥石-白云母-黑云母壓力計(jì)(Bucher-Nurminen, 1987)。但是，絕大多數(shù)地質(zhì)溫度計(jì)、壓力計(jì)與流體無關(guān)，因此分別稱為礦物溫度計(jì)、礦物壓力計(jì)。

天然變質(zhì)巖、巖漿巖中的大多數(shù)礦物，都是礦物晶格中混有“雜質(zhì)”離子的“固態(tài)溶液”(固溶體，solid solution)。在一定的P-T范圍內(nèi)，溫度或壓力的變化，會(huì)導(dǎo)致礦物組合中各種礦物的化學(xué)成分發(fā)生變化，但既不會(huì)有舊相的消失，也不會(huì)有新相的出現(xiàn)。礦物溫度計(jì)、壓力計(jì)就是利用了此類“緩沖”特性建立的。礦物成分連續(xù)變化的此種模式礦物反應(yīng)，又稱為連續(xù)反應(yīng)(continuous reaction)、滑動(dòng)反應(yīng)(sliding reaction)、緩沖反應(yīng)(buffer reaction)。

例如，變質(zhì)泥質(zhì)巖中常用的石榴子石-黑云母(GB)溫度計(jì)，是建立在石榴子石-黑云母之間Fe2+-Mg2+離子交換的滑動(dòng)反應(yīng)(Ferry and Spear, 1978)基礎(chǔ)上：

(1)

滑動(dòng)反應(yīng)(1)左側(cè)為相對(duì)低溫礦物組合，右側(cè)為相對(duì)高溫礦物組合(本文采用沈其韓(2009)、Whitney and Evans (2010)推薦的礦物縮寫代碼)。對(duì)天然巖石的研究及石榴子石-黑云母之間Fe2+-Mg2+離子交換的相平衡實(shí)驗(yàn)都證實(shí)，給定壓力條件下，隨著溫度的升高，石榴子石中Fe2+/Mg2+比值越來越小，黑云母中Fe2+/Mg2+比值越來越大，實(shí)際上就是石榴子石(Grt)中的鐵鋁榴石(Alm)、鎂鋁榴石(Prp)與黑云母(Bt)中的鐵云母(Ann)、金云母(Phl)之間Fe2+、Mg2+濃度的此消彼長(zhǎng)過程。石榴子石、黑云母化學(xué)成分的變化，受壓力條件的影響很小。通常1kabr的壓力變化，引起石榴子石-黑云母溫度計(jì)的溫度變化僅有3～5℃(Holdaway, 2000)。

天然變質(zhì)泥質(zhì)巖中，石榴子石+黑云母組合穩(wěn)定的溫度范圍很寬，從綠片巖相(Ghent, 1975; McLellan, 1985; St?ubli, 1989; Owonaetal., 2011)、角閃巖相、麻粒巖相(Indares and Martignole, 1990; Bhowmik and Spiering, 2004; Vránaetal., 2005; Millonigetal., 2008; Estrada-Carmonaetal., 2009; Hallett and Spear, 2014; Yinetal., 2014)乃至超高溫麻粒巖相(Mohan and Windley, 1993; Patio Douceetal., 1993; Guoetal., 2012)條件下，它們都可穩(wěn)定共生。這表明，石榴子石-黑云母礦物對(duì)只要隨時(shí)保持Fe2+-Mg2+離子交換的平衡狀態(tài)，即不停地調(diào)整各自的化學(xué)成分，它們就能夠自綠片巖相一直“攜手”共存、滑動(dòng)演化到麻粒巖相乃至超高溫麻粒巖相條件下。因此，石榴子石-黑云母溫度計(jì)適用的溫度范圍很寬廣。實(shí)際上，天然石榴子石中還含有Ca2+、Mn2+離子，黑云母中還含有Al3+、Ti4+離子。為了準(zhǔn)確刻畫它們之間Fe2+-Mg2+離子交換的熱力學(xué)效應(yīng)、建立準(zhǔn)確的石榴子石-黑云母溫度計(jì)，需要準(zhǔn)確描述鐵鋁榴石、鎂鋁榴石、鐵云母、金云母的活度(activity)，即它們參與平衡反應(yīng)的“有效濃度”或曰“活性”。

變質(zhì)泥質(zhì)巖中常用的石榴子石-Al2SiO5礦物-斜長(zhǎng)石-石英(GASP)壓力計(jì)，建立在GASP礦物組合的滑動(dòng)反應(yīng)(Ghent, 1976)基礎(chǔ)上：

(2)

滑動(dòng)反應(yīng)(2)左側(cè)為相對(duì)低壓礦物，右側(cè)為相對(duì)高壓礦物組合。天然變質(zhì)泥質(zhì)巖中，Al2SiO5礦物(紅柱石，And；藍(lán)晶石，Ky；夕線石，Sil)和石英(Qz)為不含“雜質(zhì)”的純物質(zhì)。該滑動(dòng)反應(yīng)表明，在礦物組合達(dá)到熱力學(xué)平衡狀態(tài)下，如果固定了溫度條件，那么隨著壓力的升高，石榴子石中Ca2+離子濃度逐漸升高，斜長(zhǎng)石中Ca2+離子濃度逐漸降低。壓力逐漸降低的情況下，礦物成分演變規(guī)律與之相反。對(duì)于GASP礦物組合，溫度條件對(duì)石榴子石、斜長(zhǎng)石化學(xué)成分的控制程度較弱。通常50℃的溫度變化，引起的GASP的壓力變化為0.8～1.0kbar(Holdaway, 2001)。

對(duì)天然變質(zhì)泥質(zhì)巖的研究表明，GASP礦物組合穩(wěn)定的壓力范圍較大。從低壓(1～2kbar; Briggs and Foster, 1992; Percival and Skulski, 2000; Zeh and Holness, 2003)、中壓(5～8kbar)到高壓(10～17kbar; Lemennicieretal., 1996; Liu and Zhong, 1997; Ding and Zhong, 1999; Ghebreab, 1999; Fosteretal., 2002; Vránaetal., 2005, 2013; Owonaetal., 2011; Yinetal., 2014; 范文壽等, 2018)條件下，該礦物組合都可穩(wěn)定共生，因此GASP壓力計(jì)應(yīng)用的壓力范圍也比較寬廣。石榴子石、斜長(zhǎng)石中并非只含有Ca2+離子，還含有其它種類離子。只有準(zhǔn)確刻畫石榴子石中鈣鋁榴石的活度，以及斜長(zhǎng)石中鈣長(zhǎng)石的活度，才能建立準(zhǔn)確的GASP壓力計(jì)。

變質(zhì)泥質(zhì)巖中，除了少數(shù)溫度計(jì)(Green and Usdansky, 1986; Vidaletal., 2005; Wu and Chen, 2015a, b)與石榴子石無關(guān)外，常用的溫度計(jì)(石榴子石-黑云母溫度計(jì), Holdaway, 2000; 石榴子石-堇青石溫度計(jì), Dwivedietal., 1998; 石榴子石-紫蘇輝石溫度計(jì), Aranovich and Berman, 1997)都與石榴子石直接有關(guān)；除了極少數(shù)壓力計(jì)(Bucher-Nurminen, 1987; Wu, 2020)與石榴子石、斜長(zhǎng)石無關(guān)外，其余壓力計(jì)(Ghent and Stout, 1981; Bohlenetal., 1983; Robinson, 1983; Hodges and Crowley, 1985; Bohlen and Liotta, 1986; Koziol and Newton, 1988; Hoisch, 1990, 1991; Koziol and Bohlen, 1992; Nicholsetal., 1992; Holdaway, 2001; Wuetal., 2004a, b; Wu and Zhao, 2006a, b, 2007a, b; Wu, 2015, 2017, 2018, 2019)都離不開石榴子石或斜長(zhǎng)石。因此，石榴子石、黑云母、斜長(zhǎng)石活度的準(zhǔn)確限定，對(duì)溫度計(jì)和壓力計(jì)的準(zhǔn)確度至關(guān)重要。本文以變質(zhì)泥質(zhì)巖中常用的GB溫度計(jì)和GASP壓力計(jì)為例，探討活度的效應(yīng)。當(dāng)然，限于實(shí)例數(shù)量，本文只能起到拋磚引玉的作用。

1 活度的描述

1.1 固溶體的活度

混合前后，如果溶液(溶劑+溶質(zhì))總的體積不發(fā)生改變、混合過程中又沒有熱效應(yīng)，這就是理想溶液(包括液態(tài)溶液、固態(tài)溶液)。理想固溶體當(dāng)然也要求礦物晶體中不同種類陽離子(離子對(duì))或陰離子(陰離子對(duì))發(fā)生相互替代、混合時(shí)，沒有引起礦物晶體體積的變化，混合過程中也沒有熱效應(yīng)的存在。實(shí)際上，只有稀溶液的熱力學(xué)性質(zhì)才接近理想溶液。一系列研究證明，在允許誤差條件下，采用理想固溶體模型研究地質(zhì)問題，也能取得比較符合地質(zhì)事實(shí)的結(jié)論。但是，很多情況下，由于固溶體中離子(離子對(duì))置換、混合的熱力學(xué)效應(yīng)并非一定是“理想”的，采用理想固溶體模型會(huì)引起較大的誤差，此時(shí)必須采用“非理想”固溶體模型來描述一般固溶體中離子(對(duì))置換、混合的熱力學(xué)效應(yīng)，即引入活度。

對(duì)于非理想溶液，人們采用“活度”替代理想溶液中的濃度，活度因此類比為“有效濃度”。一般將活度表達(dá)為理想項(xiàng)(“理想活度”Xi)與非理想項(xiàng)(“活度系數(shù)”γi)兩項(xiàng)的乘積：ai=Xiγi。這是個(gè)已經(jīng)廣泛接受的習(xí)慣作法，是對(duì)活度的近似數(shù)學(xué)表述?；疃认禂?shù)γi一般描述為礦物成分和P-T條件的多項(xiàng)式。也可以說，活度系數(shù)把固溶體偏離理想溶液性質(zhì)的所有非理想因素(“超額項(xiàng)”)，都“一股腦”地包括進(jìn)去了。

1.2 石榴子石、黑云母、斜長(zhǎng)石固溶體的晶體化學(xué)特征

2 活度對(duì)石榴子石-黑云母溫度計(jì)的作用

2.1 根據(jù)不同的活度模型組合建立石榴子石-黑云母溫度計(jì)

溫度計(jì)中的礦物固溶體模型及其熱力學(xué)參數(shù)，大多是在建立溫度計(jì)的同時(shí)確定的，所以許多固溶體礦物的活度是相互關(guān)聯(lián)的。此種條件下，單獨(dú)討論某種礦物固溶體活度對(duì)溫度計(jì)的影響，往往難以實(shí)現(xiàn)。

這里以建立石榴子石-黑云母溫度計(jì)的逆轉(zhuǎn)反應(yīng)(reversed experiment)實(shí)驗(yàn)數(shù)據(jù)(Ferry and Spear, 1978; Perchuk and Lavrent’eva, 1983; Gessmannetal., 1997)為例。采用不同的石榴子石、黑云母活度模型組合，算得溫度計(jì)中各種熱力學(xué)參數(shù)(表1)，由此可以建立不同版本的石榴子石-黑云母溫度計(jì)(圖1)。表1中，“理想固溶體”的活度等于濃度；“Holdaway (2000)”指的是該文獻(xiàn)所描述的Fe-Mg-AlⅥ-Ti四元非理想黑云母固溶體模型，以及Fe-Mg-Ca-Mn四元非理想石榴子石活度模型。

圖1 采用不同活度模型組合(表1)標(biāo)定的石榴子石-黑云母溫度計(jì)計(jì)算溫度與實(shí)驗(yàn)溫度對(duì)比各小圖圖號(hào)分別對(duì)應(yīng)表1中各行所示的活度模型組合代碼. 圖中代碼FS78、PL83、G97分別代表Ferry and Spear (1978)、Perchuk and Lavrent’eva (1983)、Gessmann et al. (1997)實(shí)驗(yàn)數(shù)據(jù). 短線代表溫度計(jì)重現(xiàn)實(shí)驗(yàn)溫度的誤差范圍(±50℃)Fig.1 Experimental versus calculated temperatures determined by different garnet-biotite geothermometers adopting different activity models listed in Table 1Names of the sub-figures correspond to the activity sets listed in Table 1. Symbols FS78, PL83, G97 represent the experimental temperatures of Ferry and Spear (1978), Perchuk and Lavrent’eva (1983) and Gessmann et al. (1997), respectively. Dashed lines represent ±50℃ deviations

采用不同的活度組合所建立的4種石榴子石-黑云母溫度計(jì)，在回算實(shí)驗(yàn)溫度方面的質(zhì)量有所不同(表1、圖1)。單從計(jì)算溫度與實(shí)驗(yàn)溫度的相關(guān)系數(shù)、標(biāo)準(zhǔn)差而論，活度模型組合自優(yōu)而劣的順序是：(d)>(b)>(a)>(c)。將石榴子石、黑云母同時(shí)作為非理想固溶體(d)得到的溫度計(jì)，重現(xiàn)實(shí)驗(yàn)溫度的效果最佳；將石榴子石作為理想固溶體的計(jì)算效果最差，這意味著石榴子石的確是非理想固溶體。

表1 基于實(shí)驗(yàn)數(shù)據(jù)，采用不同活度模型獲得的石榴子石-黑云母溫度計(jì)的熱力學(xué)參數(shù)Table 1 Thermodynamic parameters of different garnet-biotite thermometers obtained by adopting different activity models of garnet and biotite, based on experimental data

2.2 不同版本石榴子石-黑云母溫度計(jì)的應(yīng)用情況

礦物溫度計(jì)的準(zhǔn)確度，表現(xiàn)在它能準(zhǔn)確反映客觀地質(zhì)事實(shí)，也是其實(shí)用性的必要條件。例如，優(yōu)質(zhì)的礦物溫度計(jì)應(yīng)該能準(zhǔn)確反映遞增變質(zhì)帶、倒轉(zhuǎn)變質(zhì)帶、熱接觸變質(zhì)暈圈中，不同變質(zhì)地帶變質(zhì)溫度的規(guī)律性變化(Wu and Cheng, 2006)。這里就以表1中根據(jù)4種不同活度組合建立的石榴子石-黑云母溫度計(jì)為例，通過對(duì)天然變質(zhì)帶的應(yīng)用，判斷不同版本溫度計(jì)的實(shí)用性。以下P-T計(jì)算中，石榴子石-黑云母溫度計(jì)與GASP壓力計(jì)迭代求解，即同時(shí)計(jì)算P-T條件。

2.2.1 遞增變質(zhì)帶

巴羅型(Barrovian type)遞增變質(zhì)帶屬于中壓變質(zhì)相系，例如美國(guó)北愛達(dá)荷州Snow Peak地區(qū)(圖2a)、中國(guó)四川丹巴地區(qū)(圖2b)、美國(guó)紐約州Dutchess地區(qū)(圖2c)、哥倫比亞Silgará地區(qū)(圖2d)。巴肯型(Buchan type)遞增變質(zhì)帶屬于低壓變質(zhì)相系，例如加拿大Wopmay造山帶(圖2e)、美國(guó)北卡羅來納州藍(lán)嶺地區(qū)(圖2f)。不同版本的石榴子石-黑云母溫度計(jì)都反映了各個(gè)遞增變質(zhì)帶溫度逐漸增加的地質(zhì)事實(shí)(圖2)。

圖2 不同版本石榴子石-黑云母溫度計(jì)應(yīng)用于遞增變質(zhì)帶的情況圖中橫坐標(biāo)自左至右表示變質(zhì)程度逐漸增加的變質(zhì)帶劃分情況，不同顏色線條分別對(duì)應(yīng)表1中各行所示的活度模型組合代碼. (a)美國(guó)北愛達(dá)荷州Snow Peak地區(qū)巴羅型遞增變質(zhì)帶(Lang and Rice, 1985)；(b)中國(guó)四川丹巴巴羅型遞增變質(zhì)帶(Huang et al., 2003)；(c)美國(guó)紐約州Dutchess地區(qū)巴羅型遞增變質(zhì)帶(Whitney et al., 1996)；(d)哥倫比亞Santander地體Silgará地區(qū)巴羅型遞增變質(zhì)帶(Ríos et al., 2003)；(e)加拿大Wopmay造山帶巴肯型遞增變質(zhì)帶(St-Onge and Davis, 2017)；(f)美國(guó)北卡羅來納州藍(lán)嶺地區(qū)巴肯型遞增變質(zhì)帶(Eckert et al., 1989)Fig.2 Application of the different versions of the garnet-biotite geothermometers to prograde metamorphic terranesThe metamorphic grade increases from left to right, as depicted in the abscissa. Different color curves in each sub-figure correspond to the activity sets listed in Table 1. (a) Barrovian-type metamorphic sequence of the Snow Peak area, Northern Idaho, USA (Lang and Rice, 1985); (b) Barrovian-type metamorphic sequence of the Danba area, Sichuan, China (Huang et al., 2003); (c) Barrovian-type metamorphic sequence of the Dutchess County, New York, USA (Whitney et al., 1996); (d) Barrovian-type metamorphic sequence of the Silgará Formation, Santander Massif, Colombian Andes (Ríos et al., 2003); (e) Buchan-type metamorphic sequence of the Wopmay Orogen, Canada (St-Onge and Davis, 2017); and (f) Buchan-type metamorphic sequence of the southern Blue Ridge, North Carolina, USA (Eckert et al., 1989)

但是，不同溫度計(jì)之間仍然存在明顯的差別。除了美國(guó)藍(lán)嶺地區(qū)遞增變質(zhì)帶(圖2f)之外，模型(d)計(jì)算的溫度比其它組合得到的溫度都高，而模型(a)得出的溫度普遍比其它模型都低(圖2)。其中，模型(d)建立的石榴子石-黑云母溫度計(jì)回算實(shí)驗(yàn)溫度最準(zhǔn)確。由此看來，模型(a)建立的溫度計(jì)計(jì)算結(jié)果系統(tǒng)性偏低，即把石榴子石、黑云母同時(shí)作為理想固溶體的作法，是不可取的。

2.2.2 倒轉(zhuǎn)變質(zhì)帶

倒轉(zhuǎn)變質(zhì)帶是指以今天的相對(duì)高程作為參照，出露位置越高的巖石變質(zhì)程度越高的現(xiàn)象。這種現(xiàn)象與地下“層位越低，巖石變質(zhì)程度越高”的現(xiàn)象相悖，因此被稱為倒轉(zhuǎn)變質(zhì)帶(inverted metamorphic sequence/terrane)，其實(shí)質(zhì)一般是正常的遞增變質(zhì)帶被后期的地質(zhì)構(gòu)造改變所致。倒轉(zhuǎn)變質(zhì)帶的實(shí)例有錫金高喜馬拉雅倒轉(zhuǎn)變質(zhì)帶(圖3a)、美國(guó)阿拉斯加州Juneau地區(qū)倒轉(zhuǎn)變質(zhì)帶(圖3b)、尼泊爾高喜馬拉雅倒轉(zhuǎn)變質(zhì)帶(圖3c)?？梢钥闯?，4種版本的石榴子石-黑云母溫度計(jì)都準(zhǔn)確識(shí)別了自高而低的溫度遞降現(xiàn)象(圖3, 自右向左)。

圖3 不同版本石榴子石-黑云母溫度計(jì)應(yīng)用于倒轉(zhuǎn)變質(zhì)帶的情況圖中橫坐標(biāo)自右至左表示變質(zhì)程度逐漸降低的變質(zhì)帶劃分情況，不同顏色線條分別對(duì)應(yīng)表1中各行所示的活度模型組合代碼. (a)錫金高喜馬拉雅倒轉(zhuǎn)變質(zhì)帶(Dasgupta et al., 2004)；(b)美國(guó)阿拉斯加州Juneau地區(qū)倒轉(zhuǎn)變質(zhì)帶(Himmelberg et al., 1991)；(c)尼泊爾高喜馬拉雅倒轉(zhuǎn)變質(zhì)帶(Imayama et al., 2010)Fig.3 Application of the different versions of the garnet-biotite geothermometers to inverted metamorphic terranesThe metamorphic grade decreases from right to left, as depicted in the abscissa. Different color curves in each sub-figure correspond to the activity sets listed in Table 1. (a) inverted metamorphic sequence in the Sikkim Himalayas (Dasgupta et al., 2004); (b) inverted metamorphic isograds in the western metamorphic belt, Juneau, Alaska, USA (Himmelberg et al., 1991); and (c) inverted metamorphic sequence of the Nepal Himalaya (Imayama et al., 2010)

除了尼泊爾高喜馬拉雅倒轉(zhuǎn)變質(zhì)帶(圖3c)之外，模型(a)得出的溫度普遍比其它模型都低(圖3)。這里再次看到把石榴子石、黑云母同時(shí)作為理想固溶體(模型a)的溫度計(jì)，計(jì)算溫度系統(tǒng)性偏低。

2.2.3 接觸變質(zhì)帶

接觸變質(zhì)作用是指侵入巖體烘烤圍巖，在巖體外接觸帶巖石中形成變質(zhì)程度向巖體逐漸遞增的變質(zhì)帶。盡管接觸變質(zhì)帶寬度不大(一般在500m以內(nèi))，但巖石變質(zhì)程度的分帶現(xiàn)象一般比較明顯。這里以法國(guó)東Rouergue地區(qū)接觸變質(zhì)帶(圖4a)、愛爾蘭Ardara接觸變質(zhì)帶(圖4b)、美國(guó)緬因州Farmington區(qū)域熱接觸變質(zhì)帶(圖4c)為例。從圖4可以看出，雖然不同活度模型組合建立的石榴子石-黑云母溫度計(jì)都很好地區(qū)別了變質(zhì)溫度逐漸增加的現(xiàn)象，但模型(a)計(jì)算的溫度仍然系統(tǒng)性偏低(圖4)。

圖4 不同版本的石榴子石-黑云母溫度計(jì)應(yīng)用于接觸變質(zhì)帶的情況圖中橫坐標(biāo)自左至右表示變質(zhì)程度逐漸增加的變質(zhì)帶劃分情況，不同顏色線條分別對(duì)應(yīng)表1中各行所示的活度模型組合代碼. (a)法國(guó)French Massif Central 東Rouergue地區(qū)接觸變質(zhì)帶(Delor et al., 1984)；(b)愛爾蘭Ardara接觸變質(zhì)帶(Homam, 2005)；(c)美國(guó)緬因州Farmington區(qū)域熱接觸變質(zhì)帶(Holdaway et al., 1988)Fig.4 Application of the different versions of the garnet-biotite geothermometers to thermal contact aureolesThe metamorphic grade increases from left to right, as depicted in the abscissa. Different color curves in each sub-figure correspond to the activity sets listed in Table 1. (a) thermal aureole of the east Rouergue area, French Massif Central (Delor et al., 1984); (b) the Ardara aureole, NW Ireland (Homam, 2005); and (c) regional thermal contact aureoles in the Farmington area, west-central Maine, USA (Holdaway et al., 1988)

3 活度在GASP壓力計(jì)中的作用

3.1 GASP壓力計(jì)中的石榴子石與斜長(zhǎng)石活度模型組合

數(shù)十年來，人們進(jìn)行過多次涉及GASP組合的相平衡實(shí)驗(yàn)(Hays, 1966; Hariya and Kennedy, 1968; Hensenetal., 1975; Cresseyetal., 1978; Schmidetal., 1978; Wood, 1978, 1988; Goldsmith, 1980; Gasparik, 1984; Koziol and Newton, 1988, 1989; Koziol, 1990, 1996; Tropperetal., 2005)，其中有少數(shù)實(shí)驗(yàn)采用了含微量鈣鋁榴石成分的石榴子石固溶體。

為了分別考察石榴子石、斜長(zhǎng)石活度對(duì)該壓力計(jì)的效應(yīng)，本文采用由純鈣鋁榴石和純鈣長(zhǎng)石進(jìn)行的逆轉(zhuǎn)反應(yīng)(Hays, 1966; Hariya and Kennedy, 1968; Goldsmith, 1980; Gasparik, 1984; Koziol and Newton, 1988)所建立的GASP壓力計(jì)。這些實(shí)驗(yàn)結(jié)果相互恰合得很好，該平衡反應(yīng)的焓(ΔH)、熵(ΔS)變化量分別為-41.0kJ/mol、-146.0J/K·mol(McKenna and Hodges, 1988)。這兩個(gè)數(shù)據(jù)與根據(jù)最新內(nèi)洽性熱力學(xué)數(shù)據(jù)庫(Holland and Powell, 2011)推出的數(shù)據(jù)(-41.6kJ/mol、-138.070J/K·mol)在誤差范圍內(nèi)一致。因此，本文采用其平均值(-41.3kJ/mol、-142.035J/K·mol)，同時(shí)采用GASP壓力計(jì)模式反應(yīng)的體積變化量為-6.605J/bar。因此，只含端元礦物的GASP(Ky)壓力計(jì)表達(dá)式為(其中平衡常數(shù)K=1)：

-41295.0+142.035T(K)-6.605P(bar)+RTlnK=0

(GASP-Ky)

由于天然石榴子石、斜長(zhǎng)石都是固溶體礦物，GASP (Ky)壓力計(jì)表達(dá)為：

(GASP-Ky)

結(jié)合熱力學(xué)數(shù)據(jù)庫(Holland and Powell, 2011)，得到含有夕線石、紅柱石的GASP壓力計(jì)表達(dá)式分別如下：

(GASP-Sil)

(GASP-And)

人們已建立了不少石榴子石活度模型，其中有三種比較準(zhǔn)確的模型(Berman and Aranovich, 1996; Gangulyetal., 1996; Mukhopadhyayetal., 1997)。這三種石榴子石活度模型都采用了非對(duì)稱參數(shù)來描述石榴子石中陽離子混合的熱力學(xué)性質(zhì)。不同的是，Berman and Aranovich (1996)、Mukhopadhyayetal. (1997)模型描述的是Fe-Mg-Ca三元非理想石榴子石固溶體，但不同活度模型中對(duì)應(yīng)的各種熱力學(xué)參數(shù)卻有較大差別。Gangulyetal. (1996)模型描述的是Fe-Mg-Ca-Mn四元非理想石榴子石固溶體，但其Fe-Mg-Ca次級(jí)三元系活度模型對(duì)應(yīng)的熱力學(xué)參數(shù)同前兩種也有大的差異。常用的三種斜長(zhǎng)石固溶體活度模型(Fuhrman and Lindsley, 1988; Elkins and Grove, 1990; Beniseketal., 2010)描述的是Ca-Na-K三元非理想斜長(zhǎng)石固溶體，但斜長(zhǎng)石不同活度模型中對(duì)應(yīng)的各種熱力學(xué)參數(shù)也各自明顯不同。因此，不同的石榴子石、斜長(zhǎng)石活度模型組合能夠“組裝”出不同版本的GASP壓力計(jì)表達(dá)式。

與此同時(shí)，在建立不同版本的GASP壓力計(jì)時(shí)，Holdaway (2001)采用了Fe-Mg-AlVI-Ti四元非理想黑云母活度模型(H00)，但與不同石榴子石活度模型匹配的黑云母熱力學(xué)參數(shù)略有不同，以減少溫度計(jì)匹配壓力計(jì)時(shí)的誤差。表2列出了13種石榴子石、斜長(zhǎng)石、黑云母活度模型的組合，對(duì)應(yīng)13種版本的GASP壓力計(jì)。

表2 GASP壓力計(jì)中不同的活度模型組合Table 2 Different activity models in different GASP geobarometers

3.2 不同版本GASP壓力計(jì)的應(yīng)用效果

3.2.1 含Al2SiO5礦物的變質(zhì)泥質(zhì)巖

根據(jù)保持了熱力學(xué)平衡狀態(tài)下礦物成分的一些石榴子石+黑云母+斜長(zhǎng)石+石英+Al2SiO5礦物組合，包括含紅柱石(Hoisch, 1991; Zeh and Holness, 2003; 張阿利等, 2004; Pattison and Vogl, 2005)、夕線石(Grew, 1981; Hoisch, 1991; Raase and Schenk, 1994; Percival and Skulski, 2000; Zehetal., 2004; Estrada-Carmonaetal., 2009)和藍(lán)晶石(McClellandetal., 1991; Lemennicieretal., 1996; Fraseretal., 2000; Rütti, 2001; Hacker and Gans, 2005)的變質(zhì)泥質(zhì)巖樣品，通過應(yīng)用判斷GASP壓力計(jì)的準(zhǔn)確度。從圖5a-g可以看出，除了活度組合(b)、(d)外，其余將石榴子石或斜長(zhǎng)石當(dāng)作理想固溶體的GASP壓力計(jì)，都明顯過高估算了壓力，將含紅柱石或夕線石的變質(zhì)泥質(zhì)巖錯(cuò)誤地投進(jìn)了藍(lán)晶石穩(wěn)定域(圖5a, c, e-g)。采用非理想活度模型組合的各種GASP壓力計(jì)，都將變質(zhì)泥質(zhì)巖準(zhǔn)確投入了Al2SiO5礦物各自的穩(wěn)定域(圖5h-m)。

圖5 不同活度模型組合建立的GASP壓力計(jì)應(yīng)用于含Al2SiO5礦物變質(zhì)泥質(zhì)巖的情況小圖編號(hào)(a)-(m)分別對(duì)應(yīng)于表2中的活度模型組合. Al2SiO5礦物平衡線：實(shí)線據(jù)Holdaway and Mukhopadhyay (1993)，虛線據(jù)Pattison (1992)Fig.5 Application of the different GASP geobarometers to Al2SiO5 phase-bearing metapeliteNames of the sub-figures (a)-(m) correspond to the respective activity sets (a)-(m) listed in Table 2. Al2SiO5 phase diagram: the sold lines are from Holdaway and Mukhopadhyay (1993), whereas the dashed line is from Pattison (1992)

3.2.2 愛爾蘭Ardara接觸變質(zhì)帶

范圍不寬的接觸變質(zhì)帶，其變質(zhì)壓力條件應(yīng)該基本一致(Wu and Cheng, 2006)。愛爾蘭Ardara接觸變質(zhì)帶寬度約500m，分為藍(lán)晶石帶、紅柱石帶、夕線石帶(Homam, 2005)。不同版本的石榴子石-黑云母溫度計(jì)都反映了溫度的遞增變化(圖4b、圖6)。除藍(lán)晶石帶以外，六種GASP壓力計(jì)均反映這些巖石的變質(zhì)壓力基本一致(4.03±1.07kbar, 圖6)。藍(lán)晶石帶壓力明顯高于紅柱石帶和夕線石帶(圖6)，說明藍(lán)晶石帶可能來自接觸變質(zhì)帶的更深部位(?)，并且熱源較淺即中-上地殼層次？值得注意的是，這6種版本的GASP壓力計(jì)估算的含紅柱石巖石變質(zhì)壓力均偏高(圖6)。

圖6 不同活度模型組合建立的GASP壓力計(jì)應(yīng)用于愛爾蘭Ardara接觸變質(zhì)帶的情況小圖編號(hào)(a)-(f)分別對(duì)應(yīng)于表2中的活度模型組合(h)-(m). Al2SiO5礦物平衡線：實(shí)線據(jù)Holdaway and Mukhopadhyay (1993)，虛線據(jù)Pattison (1992)Fig.6 Application of the different GASP geobarometers to metapelite in the Ardara aureole, NW Ireland (Homam, 2005)Names of the sub-figures correspond to activity models (h)-(m) listed in Table 2. Al2SiO5 phase diagram: the sold lines are from Holdaway and Mukhopadhyay (1993), whereas the dashed line is from Pattison (1992)

3.2.3 美國(guó)緬因州Farmington區(qū)域熱接觸變質(zhì)帶

該地區(qū)有多個(gè)泥盆紀(jì)-石炭紀(jì)花崗巖體侵入，圍繞這些巖體外接觸帶形成一系列接觸變質(zhì)帶，稱為“區(qū)域熱接觸變質(zhì)帶”(Holdawayetal., 1988)。六種版本的GASP壓力計(jì)得出的變質(zhì)壓力介于4.55±0.59kbar～4.96±0.59kbar之間(圖7)，在誤差范圍內(nèi)基本一致，也都正確反映了這些巖石位于夕線石穩(wěn)定域的地質(zhì)事實(shí)(圖7)。

圖7 不同活度模型組合建立的GASP壓力計(jì)應(yīng)用于區(qū)域熱接觸變質(zhì)帶的情況小圖編號(hào)(a)-(f)分別對(duì)應(yīng)于表2中的活度模型組合(h)-(m). Al2SiO5礦物平衡線：實(shí)線據(jù)Holdaway and Mukhopadhyay (1993)，虛線據(jù)Pattison (1992)Fig.7 Application of the different GASP geobarometers to metapelite in the regional contact aureoles, west-central Maine (Holdaway et al., 1988)Names of the sub-figures (a)-(f) correspond to activity models (h)-(m) listed in Table 2. Al2SiO5 phase diagram: the sold lines are from Holdaway and Mukhopadhyay (1993), whereas the dashed line is from Pattison (1992)

4 討論

盡管礦物溫度計(jì)、壓力計(jì)模式反應(yīng)的焓(ΔH)、熵(ΔS)、體積(ΔV)等變化量隨P-T而變化，但大量研究實(shí)踐表明，純固相平衡反應(yīng)系統(tǒng)的這些參數(shù)隨P-T演變的變化量并不大。非超高溫、非超高壓條件下，將它們當(dāng)作常數(shù)，帶來的誤差可以忽略不計(jì)。因此，礦物溫度計(jì)、壓力計(jì)研究中一般采用礦物“標(biāo)準(zhǔn)狀態(tài)下”的熱力學(xué)參數(shù)，忽略礦物的熱膨脹系數(shù)、壓縮系數(shù)、熱容。我們?cè)谟懻摶疃饶Ｐ?表1)對(duì)于GB溫度計(jì)的影響、活度模型(表2)對(duì)GASP壓力計(jì)的影響時(shí)，也是這么處理的。

溫度升高的情況下，礦物晶格增大。這樣，半徑差別大的離子或者電價(jià)不同的離子，相互替代更為容易。因此，麻粒巖相、超高溫麻粒巖相情況下，礦物固溶體的混合性質(zhì)接近理想溶液。升壓的情況與之相反。對(duì)于一般P-T條件下形成的變質(zhì)巖石，礦物固溶體不再是理想溶液，其活度需要盡可能準(zhǔn)確描述。

4.1 活度誤差被掩蓋的情況

建立GB溫度計(jì)時(shí)，對(duì)于同樣的實(shí)驗(yàn)數(shù)據(jù)，由于采用了不同的活度模型(表1)，得出的溫度計(jì)模式反應(yīng)的焓、熵、體積的變化量之間有很大的差別，也都與熱力學(xué)數(shù)據(jù)庫(Holland and Powell, 2011)的數(shù)據(jù)也有明顯差別(表1)。換句話說，建立溫度計(jì)時(shí)得到的這些熱力學(xué)數(shù)據(jù)，誤差顯然都很大。但是，除活度組合(c)之外(圖1c)，不同版本的GB溫度計(jì)回算實(shí)驗(yàn)溫度的效果近似(圖1a, b, d)，標(biāo)準(zhǔn)差介于±22～31℃(表1)。從“實(shí)用”的角度看，這幾個(gè)版本的GB溫度計(jì)效果相近。

無論是將石榴子石、黑云母同時(shí)作為理想固溶體，還是同時(shí)都作為非理想固溶體，活度模型的巨大差異(表1)，并沒有造成GB溫度計(jì)過大的誤差(圖1)，說明活度誤差可以被模式反應(yīng)的焓、熵、體積變化量的誤差所“吸收”，由此建立起“實(shí)用”的GB溫度計(jì)。實(shí)際上，大多數(shù)礦物溫度計(jì)、壓力計(jì)都是采用此類方式建立的，這也就是為什么同一種礦物在不同溫度計(jì)、壓力計(jì)中的活度差別那么大，但卻沒有造成各個(gè)溫度計(jì)、壓力計(jì)誤差過大的原因。

此外，從圖1可以看出，活度組合(c)的效果相對(duì)最差(圖1c)，計(jì)算溫度與實(shí)驗(yàn)溫度的標(biāo)準(zhǔn)差為±37℃(表1)。該模型組合已經(jīng)充分考慮黑云母的非理想混合性質(zhì)，但卻將石榴子石作為理想固溶體處理，由此造成GB溫度計(jì)誤差過大。在實(shí)際應(yīng)用中，將石榴子石、黑云母同時(shí)作為理想固溶體的GB溫度計(jì)，計(jì)算溫度總是系統(tǒng)性偏低(圖2、圖3、圖4)，說明石榴子石、黑云母的確屬于非理想固溶體。

一般說來，在建立礦物溫度計(jì)或壓力計(jì)過程中，不同礦物的活度模型如果是同時(shí)確定的，它們的活度誤差會(huì)被“掩埋”，各種礦物活度的準(zhǔn)確性是難以分別評(píng)判的。

圖8 礦物溫度計(jì)-壓力計(jì)計(jì)算結(jié)果與流體無關(guān)假想的變質(zhì)巖或巖漿巖中的礦物組合.接近平直的線條代表全部由固相礦物參與的反應(yīng)(其中彩色線條代表礦物溫度計(jì)與壓力計(jì))，曲線代表有流體參與的反應(yīng). (a-c)含流體的礦物組合中的反應(yīng):(a)達(dá)到熱力學(xué)平衡前的各種反應(yīng)；(b)該流體條件下，礦物組合中各種反應(yīng)達(dá)到完全熱力學(xué)平衡，構(gòu)成平衡反應(yīng)的“全集”；(c)僅由平緩斜率(壓力計(jì))、陡傾斜率(溫度計(jì))的固相平衡反應(yīng)線構(gòu)成的平衡反應(yīng)集合的一部分，即平衡反應(yīng)集合的“子集”；(d-f):流體成分改變后礦物組合中的各種反應(yīng): (d)達(dá)到熱力學(xué)平衡前的各種反應(yīng)；(e)新的流體條件下，礦物組合中各種反應(yīng)達(dá)到完全熱力學(xué)平衡，構(gòu)成新的平衡反應(yīng)新的“全集”；(f)僅由溫度計(jì)與壓力計(jì)組成的新的平衡反應(yīng)“子集”. 注意(b)與(e)、(c)與(f)中的溫度與壓力條件不同F(xiàn)ig.8 P-T conditions determined by fluid-absent geothermobarometers are independent of fluidTake hypothetical metamorphic or magmatic phase assemblages as examples. Near straight curves stand for fluid-absent reactions, among which the color curves are geothermobarometers, whereas the curved lines are reactions involving fluid phases. (a-c) Fluid-related reactions of the assemblage: (a) reactions prior to thermodynamic equilibria; (b) at stable composition of the fluid, all the reactions reach thermodynamic equilibrium and constitute the full cluster; (c) the geobarometer (shallow curves) and geothermometer (steep curves) constitute the partial thermodynamic cluster; (d-f) The reactions under new composition of fluid: (d) reactions prior to new thermodynamic equilibrium; (e) new full thermodynamic cluster under new fluid composition; (f) new partial thermodynamic cluster consisted by geothermobarometers. The P-T conditions of (b) and (e) or (c) and (f) are different

4.2 可以逐個(gè)檢驗(yàn)活度作用的情況

在非超高溫、非超高壓的一般變質(zhì)作用條件下，石榴子石、斜長(zhǎng)石均表現(xiàn)出明顯的非理想混合性質(zhì)，其活度的準(zhǔn)確刻畫就顯得非常必要。由于GASP壓力計(jì)中的石榴子石、斜長(zhǎng)石活度模型，不是在標(biāo)定壓力計(jì)的同時(shí)建立的，因此可以通過實(shí)際應(yīng)用檢驗(yàn)其準(zhǔn)確性。

將石榴子石、斜長(zhǎng)石同時(shí)作為理想固溶體(表2)所建立的GASP壓力計(jì)，計(jì)算壓力嚴(yán)重偏高，將含有紅柱石、夕線石的巖石全部錯(cuò)誤地投到了藍(lán)晶石穩(wěn)定域(圖5a)。將石榴子石作為理想固溶體，采用不同的斜長(zhǎng)石非理想活度模型(Fuhrman and Lindsley, 1988; Beniseketal., 2010)所建立的GASP壓力計(jì)，已能將大部分樣品投入合適的Al2SiO5礦物穩(wěn)定域(圖5b, d)，只是采用Elkins and Grove (1990)斜長(zhǎng)石活度模型的GASP壓力計(jì)效果較差(圖5c)。與此不同的是，無論采用什么樣的石榴子石非理想固溶體活度模型，只要將斜長(zhǎng)石作為理想固溶體，那么得到的GASP壓力計(jì)往往會(huì)過高估算壓力(圖5e-g)，說明斜長(zhǎng)石是明顯偏離理想溶液性質(zhì)的非理想固溶體。

如果同時(shí)將石榴子石(Berman and Aranovich, 1996; Gangulyetal., 1996; Mukhopadhyayetal., 1997)、斜長(zhǎng)石(Fuhrman and Lindsley, 1988; Elkins and Grove, 1990)作為非理想固溶體，那么無論采用什么樣的活度模型組合(表2)，所得到的不同版本的GASP壓力計(jì)，均能將含藍(lán)晶石、紅柱石、夕線石的變質(zhì)泥質(zhì)巖，準(zhǔn)確投進(jìn)各自的穩(wěn)定域(圖5h-m)。這些石榴子石活度模型、斜長(zhǎng)石活度模型中，對(duì)應(yīng)的熱力學(xué)參數(shù)分別都有很大區(qū)別；盡管活度模型差別這么大，但將它們“組裝”入GASP壓力計(jì)(表2)之后，不同版本壓力計(jì)的計(jì)算效果卻幾乎一致(圖5)。這說明，這些活度模型都有類似的可靠的準(zhǔn)確度。

4.3 缺少通用型活度模型

不同巖性中的同一種礦物，其化學(xué)成分一般有明顯差異。例如，變質(zhì)泥質(zhì)巖中的石榴子石富FeO而貧MgO，離子混合發(fā)生在六面體結(jié)點(diǎn)上；變質(zhì)基性巖(角閃巖、基性麻粒巖)中的石榴子石既富FeO也富MgO，離子混合發(fā)生在六面體結(jié)點(diǎn)上；變質(zhì)超基性巖中的石榴子石貧FeO而富MgO，離子混合既發(fā)生在六面體結(jié)點(diǎn)上，也發(fā)生在八面體結(jié)點(diǎn)上；超高壓變質(zhì)的榴輝巖中的石榴子石，離子混合既發(fā)生在六面體結(jié)點(diǎn)上，還發(fā)生在四面體結(jié)點(diǎn)上。這些情況導(dǎo)致了適用于變質(zhì)泥質(zhì)巖的石榴子石-黑云母溫度計(jì)、適用于斜長(zhǎng)角閃巖的石榴子石-角閃石溫度計(jì)、適用于基性麻粒巖的石榴子石-單斜輝石溫度計(jì)、適用于變質(zhì)超基性巖的石榴子石-橄欖石溫度計(jì)，其中的石榴子石活度都不同。目前，還沒有廣泛適用的通用型石榴子石活度模型。

從變質(zhì)演化角度看，石榴子石中不僅可以保存早期變質(zhì)階段的礦物組合，還易于保存退變質(zhì)階段的礦物組合，石榴子石本身還可以用于定年，因此石榴子石是變質(zhì)地質(zhì)學(xué)研究中最為看重的研究對(duì)象。從結(jié)晶學(xué)角度看，石榴子石是對(duì)稱性最高的等軸晶系礦物。從化學(xué)角度看，其晶體化學(xué)規(guī)律早已查明。但是，就是這樣一個(gè)應(yīng)該標(biāo)定好活度的礦物，其活度描述還很不理想。此外，角閃石既是廣泛見于角閃巖相和麻粒巖相變質(zhì)巖的礦物，還是晶體化學(xué)結(jié)構(gòu)最復(fù)雜的常見礦物，其活度模型研究更是遠(yuǎn)不理想。其它礦物活度模型的研究，當(dāng)然談不上完善。由于人們開展的巖石學(xué)相平衡實(shí)驗(yàn)還嚴(yán)重不足，導(dǎo)致高質(zhì)量的平衡礦物組合數(shù)據(jù)很是缺乏，這也直接導(dǎo)致通用型活度模型的匱乏。開發(fā)通用型礦物活度模型，應(yīng)該是今后重點(diǎn)開展的研究方向。

4.4 礦物溫度計(jì)、壓力計(jì)計(jì)算結(jié)果并不受流體的影響

既然礦物溫度計(jì)、壓力計(jì)未考慮流體，一個(gè)自然而然的疑問是，它們是否真實(shí)反映了客觀地質(zhì)事實(shí)？計(jì)算結(jié)果是否反映了含有流體的礦物組合的平衡P-T條件？

眾所周知，在變質(zhì)作用、巖漿作用過程中，如果存在流體，流體會(huì)參與巖石中的許多種類的反應(yīng)。此處以圖8為例簡(jiǎn)單說明。圖8表示任意變質(zhì)巖、巖漿巖中的礦物+流體組合演化到熱力學(xué)平衡狀態(tài)的過程。在給定溫度、壓力、流體條件下，巖石達(dá)到熱力學(xué)平衡之前，涉及固相-固相、固相-流體相的各種反應(yīng)，各自獨(dú)立，反應(yīng)線并未“聚合”在一起(圖8a)。巖石達(dá)到完全的熱力學(xué)平衡時(shí)，所有反應(yīng)線交匯于一點(diǎn)，這些平衡反應(yīng)線組成了平衡狀態(tài)的“全集”(圖8b)。從“全集”中抽提出來的不涉及流體相的溫度計(jì)與壓力計(jì)平衡反應(yīng)線(圖8c)，即全部平衡反應(yīng)線的“子集”，肯定與“全集”的平衡壓力與溫度(P1,T1)條件完全一致。當(dāng)流體條件發(fā)生改變時(shí)(圖8d-f)，從未平衡演化到完全熱力學(xué)平衡，規(guī)律仍然和圖8a-c類似，最終達(dá)到新的平衡壓力與溫度(P2,T2)條件。因此可以說，只要礦物組合達(dá)到了完全的熱力學(xué)平衡，礦物溫度計(jì)與壓力計(jì)的計(jì)算結(jié)果，的確能反映客觀地質(zhì)事實(shí)，并不受是否考慮流體的影響。

類似地，如果巖石中還存在熔體的話，那么，采用只考慮固相礦物的礦物溫度計(jì)與壓力計(jì)，也同樣能準(zhǔn)確反映客觀地質(zhì)事實(shí)，仍然不受“撇開”熔體、流體的影響。

需要說明的是，圖8中所示的與流體和熔體無關(guān)的、只由固相礦物構(gòu)成的平衡反應(yīng)中，那些P/T斜率<10bar/℃的平衡反應(yīng)受溫度影響很小，才是理想的礦物壓力計(jì)；P/T斜率>80bar/℃的平衡反應(yīng)受壓力影響很小，才是理想的礦物溫度計(jì)。其余的固相平衡反應(yīng)，既不能作為壓力計(jì)，也不能作為溫度計(jì)。這也就是礦物溫度計(jì)、壓力計(jì)數(shù)量有限的原因。此外，由于微量元素、穩(wěn)定同位素、稀土元素等在固相礦物之間的分配受到壓力的影響較小，因此至少目前還難以采用這些元素建立礦物壓力計(jì)。

5 結(jié)論

(1)將石榴子石、黑云母、斜長(zhǎng)石當(dāng)作理想固溶體，得到的石榴子石-黑云母溫度計(jì)往往會(huì)過低估算溫度條件，得到的GASP壓力計(jì)往往會(huì)過高估算壓力條件，它們都存在較大的系統(tǒng)誤差；

(2)一般變質(zhì)作用條件下，固溶體礦物屬于非理想固態(tài)溶液。如果其活度刻畫不準(zhǔn)確，將直接導(dǎo)致溫度計(jì)、壓力計(jì)的系統(tǒng)誤差；

(3)比較接近真實(shí)“有效濃度”的各種礦物活度模型，盡管其中熱力學(xué)參數(shù)有所不同，但其實(shí)用性通?？梢詽M足；

(4)建立適用于多種巖石類型、溫度與壓力適用范圍大的通用型礦物活度模型，是值得開展的重要研究任務(wù)；

(5)雖然礦物溫度計(jì)、壓力計(jì)不涉及流體、熔體，但對(duì)于達(dá)到并保持了熱力學(xué)平衡狀態(tài)條件下化學(xué)成分的礦物組合而言，其計(jì)算結(jié)果的確能準(zhǔn)確反映客觀地質(zhì)事實(shí)，并不受是否考慮流體、熔體的影響。

致謝陳意教授和向華博士提出了寶貴的修改建議，提高了本文的學(xué)術(shù)水平。作者向他們致以真摯的感謝。

謹(jǐn)以拙文祝賀著名前寒武紀(jì)地質(zhì)學(xué)家、巖石學(xué)家沈其韓院士百歲大壽，敬祝恩師健康長(zhǎng)壽！