黑云母温度计能否用于估计花岗质侵入岩的结晶温度.pdf

Doi:10.20086/skw.2024.0103Jan.,2024ETMINERALOGICAACTAPETROLOGICA2024年1月石Vol.43,No.1:3746志心杂矿学岩物第1期第43 卷黑云母温度计能否用于估计花岗质侵入岩的结晶温度?徐鸿雪，汪洋（中国地质大学（北京）地球科学与资源学院，北京100083)摘要：为了探讨黑云母温度计是否适用于估计花岗质侵人岩的结晶温度，汇编了国内外典型花岗质侵人体的岩石化学和矿物化学数据，利用Henry等（2 0 0 5）的黑云母Ti饱和温度计、Li和Zhang（2 0 2 2）的黑云母全组分温度计估算了岩浆温度，并与Shao等（2 0 2 0)的锆饱和温度计估算结果进行了对比。结果表明，黑云母Ti饱和温度计估算的温度比锆饱和温度计和黑云母全组分温度计低50 2 0 0，因此认为黑云母Ti饱和温度计并不适用于估算花岗质侵人岩的结晶温度，而黑云母全组分温度计可以复现S型和I型花岗岩的锆饱和温度计结果。关键词：花岗质侵人体；结晶温度；黑云母Ti饱和温度计；黑云母温度计；锆饱和温度计中图分类号：P588.12*1文献标识码：A文章编号：10 0 0-6 52 4(2 0 2 4)0 1-0 0 37-10Can biotite geothermometer be used to estimate the crystallizationtemperature of granitic intrusive rocks?XU Hong-xue and WANG Yang(School of Earth Sciences and Resources,China University of Geosciences(Beijing),Beijing 100083,China)Abstract:In order to investigate whether the biotite thermometer is suitable for estimating crystallization tempera-ture of granite intrusive rocks,a compilation of the petrochemistry and mineral chemistry data of typical graniticintrusive rocks at domestic and abroad was made.The magma temperature was estimated using Henry et al.(2005)biotite Ti saturation geothermometer and the Li and Zhang(2022)version of the biotite geothermometer,and theresults were compared with those estimated by the zircon saturation geothermometer in Shao et al.(2020).Theresults show that the temperatures estimated by the biotite Ti saturation geothermometer are 50200 l o w e r t h a nthose estimated by the zircon saturation geothermometer and the 2022 version of the biotite geothermometer.Thebiotite Ti saturation thermometer is not suitable for estimating the crystallization temperature of granitic intrusiverocks,while the Li and Zhang(2022)version of the biotite thermometer can reproduce the results of the zircon sat-uration thermometer for S-type and I-type granites.Key words:granitic intrusive rock;crystallization temperature;biotite Ti saturation geothermometer;biotitegeothermometer;zirconium saturation geothermometerFund support:Project of China Geological Survey(D1912)收稿日期：2 0 2 3-0 7-0 5；接受日期：2 0 2 3-11-0 9；编辑：郝艳丽基金项目：中国地质调查局地质调查项目（D1912）作者简介：徐鸿雪（1999-），女，硕士研究生，地质工程专业，E-mail：x u h o n g x u e 16 3 16 3.c o m；通讯作者：汪洋（196 9-），男，主要从事岩石学和地热学研究，E-mail：a l l e n _t h a l a s s a s i n a.c o m。石38心杂矿岩岩第43 卷学物岩浆结晶温度是理解火成岩成因的重要物理参量，环境和岩浆自身温度的变化控制了源岩部分熔融的发生、母岩浆的分离结晶以及岩体侵位后的固结过程（Miller et al.，2 0 0 3；A n d e r s o n e t a l.，2 0 0 8;Putirka，2 0 0 8）。对于特定成分的岩浆，温度和压力的变化对其平衡矿物组合起着决定性的作用，进而控制了岩浆岩的主量元素组成，影响岩浆岩的微量元素含量（Janousek andMoyen，2 0 2 0）。对于广义的花岗质岩浆而言，其温度范围大致在6 50 950 之间。对于花岗质侵人岩，估算其结晶温度的主要方法包括两大类地质温度计。一类是（微量）元素饱和温度计，主要有锆饱和温度计、磷饱和温度计、钛饱和温度计、轻稀土元素饱和温度计等（WatsonandHarrison,1983;Harrison and Watson,1984;Bea etal.,1992;Pichavant et al.，1992;M o n t e l,1993;Wolf and London,1994;Jung and Pfander,2007;Stepanov et al.,2012;Boehnke et al.,2013;Shao etal.，2 0 19，2 0 2 0），近年来发展起来的锆石Ti温度计等(Janousek andMoyen，2 0 2 0)也属于此类。另一类是矿物温度计，包括钛铁矿-磁铁矿温度计、二长石温度计、辉石温度计、角闪石温度计、黑云母温度计等（Fuhrman and Lindsley，198 8；Be n i s e k e t a l.,2004;Henry et al.,2005;Ghiorso and Evans,2008;Putirka,2008,2016;Ridolfi et al.,2010;Wang etal.，2 0 2 1b；Li a n d Zh a n g，2 0 2 2）。由于黑云母在花岗质岩浆岩中广泛存在，所以黑云母Ti温度计（H e n r y e t a l.，2 0 0 5）被不少学者用于估算花岗质侵人岩的结晶温度（Cesare etal.，2 0 0 8；Ba y a t i e t a l.，2017;Moshefi et al.,2018;Azadbakht et al.,2020;Azer et al.,2020;Baidya et al.,2021;Zhao et al.,2023；王栋等，2 0 2 3；邹兴志等，2 0 2 3）。然而，Henry等（2 0 0 5的黑云母Ti温度计是基于泥质变质岩的研究确立的，能否应用于花岗岩尚存疑义，该温度计估算的花岗质岩浆的温度是否准确亦缺乏实质性检验，为此本文应用经验性的方法检验了黑云母温度计估算花岗质侵入岩形成温度的可靠性。1估算花岗质侵人岩成岩温度的常用方法概述锆元素在天然硅酸盐熔体中以微量元素的形式存在,其物理化学行为受稀溶液定律（Henry定律）制约。Watson和Harrison（198 3）基于实验岩石学研究证明，锆在熔体中的饱和度主要取决于熔体的成分和温度。相较于其他元素饱和温度计，锆饱和温度计的优势在于：长英质岩浆体系中锆石-熔体之间Zr分配系数的温度敏感程度高，而且石是花岗质岩浆的早期结晶相，化学稳定性高，结晶后不易被改造，因此研究者对岩浆体系的Zr饱和体系进行了大量实验标定，积累了丰富的实验数据（Dietrich，1968;Watson，197 9;D i c k i n s o n a n d H e s s,198 2;Harrison and Watson,1983;Watson and Harrison,1983;Ellison and Hess,1986;Keppler,1993;Bakeret al.,2002;Hanchar and Watson,2003;Boehnke etal.,2013;Zhang and Xu,2016;Gervasoni et al.,2016；Sh a o e t a l.，2 0 19，2 0 2 0）。同时，锆石是花岗质岩浆中普遍存在的副矿物相,可以在中性-酸性成分的硅饱和熔体中结晶出来。对于含有残留锆石的花岗质侵人岩，锆饱和温度反映其岩浆温度的上限值；对于不含残留锆石的花岗质侵人岩（如A型花岗岩），锆饱和温度可以代表其岩浆温度的下限值。所以，锆饱和温度计是目前用于估算花岗质侵入岩成岩温度的最常用方法，得到了广泛应用（Milleretal.,2003;Collins et al.,2016,2020)。相对于锆饱和温度计，其他的地质温度计在花岗岩成岩温度估算方面存在一些局限。对磷、钛饱和温度计的实验标定研究较少（Harrison andWat-son,1984;Pichavant et al.,1992;Montel,1993),应用实例也比较少。轻稀土元素饱和温度计主要适用于含独居石的S型花岗岩类的温度估计（Montel，1993；St e p a n o v e t a l.，2 0 12），对于不含独居石的I型和A型花岗岩，不能采用轻稀土元素饱和温度计估算其成岩温度。由于不少花岗岩不含辉石、角闪石，加之辉石温度计的精度不高（Wieseretal.，2023），辉石、角闪石温度计在对花岗质岩浆岩中的应用实例有限（Anderson etal.，2 0 0 8）。例如，典型的S型花岗岩不含普通角闪石，无法采用涉及角闪石的地质温度计估算S型花岗岩的温度。同时，针对普通角闪石成分得到的温度计不适用于碱性角闪石，所以，常见的角闪石温度计无法用于含碱性角闪石的A型花岗岩。对于估算花岗质岩浆成岩温度而言，钛铁矿-磁铁矿温度计受固相线温度之下再平衡因素的影响（Fuhrman and Lindsley，198 8；Be n i s e k e tal.，2 0 0 4；G h i o r s o a n d Ev a n s，2 0 0 8），在岩浆固结之后Fe、T i 离子仍然在钛铁矿-磁铁矿之间发生扩散交换，导致计算的温度值明显低于岩浆的固相线温度。二长石温度计存在同样的缺陷。所以，对估计39徐鸿雪等公母温度计能香否用岩的结晶温度？一计花区顶信第1期花岗质侵入岩的成岩温度而言，钛铁矿-磁铁矿温度计和二长石温度计应用得比较少。目前，锆饱和温度计已发表多个版本，主要有Watson和Harrison（198 3）（简称WH83）、Bo e h n k e 等（2 0 13）（简称B13）、G e r v a s o n i 等（2 0 16）（简称G16）、Sh a o 等（2 0 2 0）（简称S20）等。S20在先前模型基础上改进使其适用范围扩展到所有过铝质到过碱质成分的镁铁质到长英质熔体。Shao等（2 0 2 0）的研究表明，在常见的岩浆温度范围内，不管是过铝质还是过碱质熔体，S20模型的估算结果与WH83模型的偏差在5%不确定性范围内，而B13、G 16 模型对相同成分熔体的温度估计值较WH83、S2 0 模型系统偏低(图1)。综合上述，本文采用Shao等（2 0 2 0）锆饱和温度计模型作为计算花岗质侵入岩成岩温度的参照值，以此评估黑云母温度计对花岗质侵入岩的适用性。10001150abO+WH839501100口B13OG169001050过碱性正长岩5%S20_Model1(Na+K)/AIl,(Na+K+Ca)/AI15%850H100010%过铝质流纹岩￥(Na+K)/A1,(Na+K+Ca)/Al1,(Na+K+Ca)/AI1公850名100010%SI过铝质流纹岩800(Na+K)/Al,(Na+K+Ca)/Ak195020%7509007008505%65080010%-60075055020%87005001650500550600650700750800850 900950100065070075080085090095010001050110011501/t/图1锆饱和温度计（S20）与其他版本锆饱和温度计（WH83、B13、G 16)的估算结果对比数据引自Shao等（2 0 2 0）Fig.1 The results of zirconium saturation thermometer(S20)compared with those of other versions of zirconiumsaturation thermometer(WH83,B13,G16)（d a t a f r o m Sh a o e t a l.,2 0 2 0)aS20模型1与版本WH83、B13.G 16 对过铝质流纹岩计算结果对比；bS20模型1与版本WH83、B13、G 16 对过碱性正长岩计算结果对比；cS20模型2 与版本WH83、B13、G 16 对过铝质流纹岩计算结果对比；dS20模型2 与版本WH83、B13、G 16 对过碱性正长岩计算结果对比acomparison of S20 model 1 with other versions(WH83,B13,G16)on the calculation results of peraluminous rhyolite;bcomparison of S20model 1 with other versions(WH83,B13,G16)on the calculation results of peralkaline syenite;ccomparison of S20 model 2 with other versions(WH83,B13,G16)on the calculation results of peraluminous rhyolite;dcomparison of S20 model 2 with other versions(WH83,B13,G16)onthe calculation results of peralkaline syenite石40心杂矿岩学物第43卷本文的目的并非全面评估涉及花岗岩的地质温度计，对其他地质温度计的适用性等问题有待另文详述。2黑云母温度计的经验性检验2.1黑云母温度计概述黑云母（包括狭义的黑云母和金云母）是地壳和地慢中广泛存在的一种含挥发分硅酸盐矿物，是大多数岩浆岩中的主要镁铁质造岩矿物之一，可在较大压力和温度范围内保持稳定。由于黑云母晶体结构比较复杂，能够大量地赋存多种阳离子和阴离子，其化学成分可以指示其形成的温度、氧逸度、水活度等重要参数，被广泛用于反演相关的岩浆性质和岩浆过程(Saha et al.,2021;Samadi et al.,2021;Gionet al.,2022;Li and Zhang,2022)。Henry等（2 0 0 5）根据过铝质变泥质岩体系黑云母Ti含量与其形成温度的相关性，提出：t=ln(Ti)+2.359 4+1.728 3(Xmg)/b 0.33(1)式中，t为温度（），Ti为按2 2 个0 原子为单位计算阳离子数后的原子数，XMg=Mg*/（M g*+Fe 2+Fe*），b=4.6 48 2 10 ,公式适用范围为：Xmg=0.2751.000,Ti4+=0.040.60 a.p.f.u.,t=480 800。黑云母Ti温度计的精度估计，在较低温度范围内为2 4，在较高温度下提高到12。该温度计被国内外许多学者应用于估计花岗质侵人体的结晶温度（李胜荣等，2 0 0 6；Cesareetal.，2 0 0 8；郭耀宇等，2 0 15；Moshefi et al.，2 0 18；A z a d b a k h t e t a l.,2020；王栋等，2 0 2 3；解世雄等，2 0 2 3；Zhaoetal.，2023)。Henry等（2 0 0 5）指出，黑云母Ti温度计适用条件是含石墨、至少一种（黑云母之外的）过铝矿物、含钛铁矿或金红石的过铝质变质泥岩体系。虽然该温度计可能适用于过铝质成分的花岗岩，但对于其他成分的体系（例如准铝质花岗岩）而言，其实际误差可能十分显著（Henry，私人通讯，见Azadbakhtetal.,2020,p.13)。由于黑云母的化学结构复杂性，基于恒定的元素分配系数K。估计黑云母和熔体平衡温度的传统方法不可行，例如，平衡体系中黑云母和熔体之间的Mg/Fe分配系数在黑云母-熔体体系中通常不是一个恒定值。因此,Li和Zhang（2 0 2 2）在高温高压实验大数据的基础上，运用机器学习方法对岩浆成因黑云母矿物主要氧化物组分进行了回归训练，得到黑云母全组分与其形成温度、压力值的统计预测模型，提出了基于黑云母全组分的岩浆体系温压计(https:/lixiaoyan.shinyapps.io/Biotite_thermobarom-eter/）。该文作者认为，黑云母全组分温度计可以较好地复现已有的实验标定温度，而Henry等（2 0 0 5）黑云母Ti温度计则不能2.2数据来源本文汇编了来自全球不同区域的代表性花岗质侵人岩全岩和黑云母成分数据，根据原始文献描述，岩石和矿物样品新鲜，黑云母属于岩浆成因。根据岩相学和岩石地球化学特征，将这些样品分别归为S型、I型和A型花岗岩类（表1）。以Shao等表1岩体名称及来源文献Table 1Names and reference of pluton类型地区岩体/岩基数据来源英国BodminSimons et al.,2016英国CarnmenellisSimonsetal.,2016英国Carn MarthSimons et al.,2016英国Lands EndSimons et al.,2016英国CliggaSimons et al.,2016英国DartmoorSimons et al.,2016英国St.AustellSimons et al.,2016意大利Sattelspitze(Monte Sella)Thony et al.,2009Mount Qomolangma-S型喜马拉雅山Visona et al.,2012Masang Kang阿根廷TuaniDahlquist et al.,2005巴西NanuqueFerreiraetal.,2019阿根廷MazanFerreira et al.,2019阿根廷Senor de la PenaFerreira et al.,2019阿根廷CapillitasFerreiraetal.,2019智利Los TilosFerreira etal.,2019德国MitterteichSiebel,1995中国峡江县金滩Tao et al.,2018瑞士central AarBucher and Seelig,2018中国阿尔泰造山带XiboduCui et al.,2022加拿大NewBrunswickLakeGeorgeYang et al.,2002巴西Serra da GargantaNascimento etal.,2018中国秦岭丹凤蟒岭杨阳等，2 0 2 2I型中国薄竹山Wang et al.,202la中国秦岭沙河湾刘春花等，2 0 13b中国秦岭东江口、作水刘春花等，2 0 13a中国铜山岭Huang et al.,2018韩国Eopyeong钾玄质侵人岩Im et al.,2021巴西Campina doVeado rapakiviGodoy et al.,2021瑞士BristenBucher and Seelig,2018中国阿尔泰造山带WestKouanCui et al.,2022中国陂头Jiang et al.,2022A型中国黄梅尖Zhang et al.,2022尼日利亚Ririwai碱性杂岩体Kinnaird etal.,1985巴西MandiraSiachoque et al.,2021巴西Serrada GraciosaGualda and Vlach,2007喀麦隆Sabongari碱性杂岩体Njonfang etal.,201341徐鸿雪等母温度计能否用于估计花岗质侵入头岩的结晶温度？第1期(2020)的锆饱和温度计计算结果为参照值，同时采用Henry等（2 0 0 5）黑云母Ti饱和温度计、Li和Zhang(2022)黑云母全组分温度计估计对应样品的黑云母结晶温度，将其结果与参照值相比较,判断计算的黑云母结晶温度能否代表岩浆的成岩温度。2.3结果与分析(1)S型花岗岩图2 a为黑云母Ti饱和温度计（Henryetal.，2005）与锆饱和温度计（Shaoetal.，2 0 2 0）估算的S型花岗岩成岩温度对比；图2 b为黑云母全组分温度计（LiandZhang，2 0 2 2）与锆饱和温度计估算的S型花岗岩成岩温度对比。结果显示，黑云母Ti饱和温度计估算结果较S20锆饱和温度计系统偏低50 250，甚至低50 0；黑云母全组分温度计的估算结果与S20锆饱和温度计之间的偏差大多在50 以内。因此，即便是对过铝质的S型花岗岩而言，黑云母Ti饱和温度计（Henryetal.，2 0 0 5）估算的温度值也是明显偏低的，因此该温度计并不适用于估算S型花岗岩的成岩温度。(2)I型花岗岩从图3可知，大多数I型花岗岩样本的黑云母Ti饱和温度较S20锆饱和温度偏低50 2 0 0，有相当数量的样品温度估算值低于6 50（压力大于100MPa、饱和水条件下花岗岩的固相线温度）（图3a）。黑云母全组分温度计的估算结果与S20锆饱和温度计之间的偏差绝大多数在50 以内（图3b）。显然，黑云母Ti饱和温度计不适用于估算I型花岗岩的成岩温度。(3)A 型花岗岩图4a为黑云母Ti饱和温度计与S20锆饱和温度计估算的A型花岗岩成岩温度对比。从图中可以看出，黑云母Ti饱和温度计与S20锆饱和温度计估算1000950a1b290039004*5&6800+78508車9107004118000./SO-1口12￥131460015750H165007004006503006003004005006007008009001000600650700750800850900950T-ZT20/T-zr20/图2锆饱和温度计（Shaoetal.，2 0 2 0）与黑云母Ti饱和温度计（Henryetal.，2 0 0 5）（a）、黑云母全组分温度计（Li a n d Zh a n g，2 0 2 2）（b）计算S型花岗岩温度结果对比图Fig.2 Comparison diagram of S-type granite temperature results calculated by zirconium saturation geothermometer(Shao et al.,2020)and biotite Ti saturation geothermometer(Henry et al.,2005)(a),biotite geothermometer(Li and Zhang,2022)(b)1峡江县金滩（中国）；2 St.Austell（英国）；3Dartmoor（英国）；4Cligga（英国）；5LandsEnd（英国）；6 CarnMarth（英国）；7Carnmenellis（英国）；8 Bodmin（英国）；9Sattelspitze（M o n t e Se l l a）（意大利）；10 Mount Qomolangma-Masang Kang（喜马拉雅山）；11Serior de la Peia(阿根廷）；12 Mazan(阿根廷）；13Capillitas（阿根廷）；14Tuani（阿根廷）；15Mitterteich(德国）；16 Nanuque(巴西）1Jintan,Xiajiang County(China);2St.Austell(UK);3Dartmoor（U K);4Cl i g g a (U K);5La n d s En d （U K);6 Ca r n M a r t h(UK);7Carnmenellis(UK);8Bodmin(UK);9Sattelspitze(Monte Sella)(Italy);10Mount Qomolangma-Masang Kang(Himalayas);11Senor de la Pena（A r g e n t i n a);12 M a z a n （A r g e n t i n a);13Ca p i l l i t a s （A r g e n t i n a);14-T u a n i （A r g e n t i n a);15M i t t e r t e i c h （G e r ma n y);16Nanuque(Brazil)石42心杂矿石岩学第43 卷物900950ab+铜山岭（中国）V薄竹山（中国）850秦岭沙河湾（中国）秦岭丹凤蟒岭（中国）900秦岭东江口、作水（中国）800阿尔泰造山带Xibodu（中国）NewBrunswickLakeGeorge(加拿大)Eopyeong钾玄质侵入岩（韩国）850750XCentralAar(瑞士)米Serra daGarganta(巴西)0./SO-10./-1700800-650750600700550500650500550600650700750800850900650700750800850900950T-zr20/T-z20/图3锆饱和温度计（Shaoetal.，2 0 2 0）与黑云母Ti饱和温度计（Henryetal.，2 0 0 5）（a）、黑云母全组分温度计（Li a n d Zh a n g，2 0 2 2）（b)计算I型花岗岩温度结果对比图Fig.3 Comparison diagram of I-type granite temperature results calculated by zirconium saturation geothermometer(Shao et al.,2020)and biotite Ti saturation geothermometer(Henry et al.,2005)(a),biotite geothermometer(Li and Zhang,2022)(b)11001100ab陂头（中国）1000X黄梅尖（中国）+阿尔泰造山带West-Kouan（中国）Mandira(巴西)1000900Serrada Graciosa(巴西)CampinadoVeadorapakivi(巴西)Sabongari碱性杂岩体（喀寿隆）800Ririwaicomplex(尼日利亚)Bristen(瑞t:)900700O./SO-1/亿-18006005007004003006003004005006007008009001000110060070080090010001100T-zT20/T-zr20/图4锆饱和温度计（Shaoetal.，2 0 2 0）与黑云母Ti饱和温度计（Henryetal.，2 0 0 5）（图a）、黑云母全组分温度计（Li a n d Zh a n g，2 0 2 2）（图b）计算A型花岗岩温度结果对比图Fig.4 Comparison diagram of A-type granite temperature results calculated by zirconium saturation geothermometer(Shao et al.,2020)and biotite Ti saturation geothermometer(Henry et al.,2005)（a),b i o t i t e g e o t e r mo me t e r (Li a n d Zh a n g,2 0 2 2)（b)值之间具粗略的线性相关性，但黑云母Ti饱和温度明显低于锆饱和温度（约10 0 30 0）。黑云母全组分温度计与S20锆饱和温度计估算的A型花岗岩成岩温度之间则不存在明显的相关性（图4b）。锆饱和温度计估算的A型花岗岩成岩温度在6 7 0 110 0 之间变化，黑云母全组分温度计得到的温度值介于7 2 5 8 0 0 之间。A型花岗质岩浆碱质含量高、对Zr的溶解能力43徐鸿雪等母温度计能否用的结晶温度？第1期强(Gervasoni et al.，2 0 16;Sh a o e t a l.，2 0 2 0),因此A型花岗岩往往显示高的锆饱和温度，在某些情况下（如无继承锆石）A型花岗质岩浆中Zr没有饱和，故而锆饱和温度计估算的成岩温度值偏低。对于准铝质和部分过碱性A型花岗岩，黑云母可视为近固相线矿物，其结晶温度低于液相线温度、高于固相线温度（Scailletetal.，2 0 16）。所以黑云母全组分温度计给出的是A型花岗岩中黑云母的结晶温度，约束了相应岩浆体系固相线温度的上限值。3结论（1）无论是S型、I型抑或A型花岗岩类，Henry等（2 0 0 5）黑云母Ti饱和温度计得到的成岩温度估算值较锆饱和温度计和黑云母全组分温度计的估算结果系统偏低50 2 0 0 以上，这表明对于花岗岩体系而言，黑云母Ti饱和温度计的估算结果地质意义不明确。因此Henry等（2 0 0 5）黑云母Ti饱和温度计不适用于估算花岗质侵入岩的结晶温度。（2）Li 和Zhang（2 0 2 2）发表的黑云母全组分温度计可以较好地复现S型和I型花岗质类的锆饱和温度计估算结果。碱质含量高的A型花岗质岩浆对Zr的溶解能力很强，导致A型花岗岩具有变化幅度较大的高Zr含量,其Zr饱和温度变化范围较大。对于A型花岗岩类，黑云母全组分温度计得到的温度值是黑云母在富碱质的花岗质岩浆体系中的结晶温度，在大多数情况下介于富碱质花岗质岩浆体系的液相线温度与固相线温度之间。（3）Li 和Zhang（2 0 2 2）黑云母全组分温压计可作为一种独立于锆饱和温度计的方法用于估算S型和I型花岗质侵人岩的成岩温度，也可以用于约束A型花岗岩的固相线温度。ReferencesAnderson J L,Barth A P,Wooden J L,et al.2008.Thermometers andthermobarometers in granitic systemsJ.Reviews in Mineralogy andGeochemistry,69:121142.Azadbakht Z,Lentz D R,McFarlane C R M,et al.2020.Using magmat-ic biotite chemistry to differentiate barren and mineralized Silurian-Devonian granitoids of New Brunswick,Canada J.Contributions toMineralogy and Petrology,175:69.Azer M K,Surour A A,Madani A A,et al.2020.Mineralogical and geo-chemical constraints on the postcollisional mafic magmatism in the A-rabian-Nubian Shield:An example from the El-Bakriya Area,centralEastern Desert,Egypt J.The Journal of Geology,130(3):209230.Baidya A S,Pal D C and Upadhyay D.2021.Biotite chemistry and min-eral association as an indicator of redox conditions in the iron oxideCu-Au(IOCG)system:Constraints from the Khetri Copper Belt.western IndiaJ.Ore Geology Reviews,139(5):104544.Baker D R,Conte A,Freda C,et al.2002.The effect of halogens on Zrdiffusion and zircon dissolution in hydrous metaluminous graniticmeltsJ.Contributions to Mineralogy and Petrology,142:666678.Bayati M,Esmaeily D,Maghdour-Mashhour R,et al.2017.Geochemis-try and petrogenesis of Kolah-Ghazi granitoids of Iran:Insights intothe Jurassic Sanandaj-Sirjan magmatic arc J.Geochemistry,77(2):281302.Bea F,Fershtater G B and Corretge L G.1992.The geochemistry ofphosphorus in granite rocks and the effects of aluminium J.Lithos,29:43 56.Benisek A,Kroll H and Cemic L.2004.New developments in two-feld-spar thermometry J.American Mineralogist,89(10):1 4961504.Boehnke P,Watson E B,Trail D,et al.2013.Zircon saturation re-revis-itedJ.Chemical Geology,351:324334.Bucher K and Seelig U.2018.Bristen granite:A highly differentiated,fluorite-bearing A-type granite from the Aar massif,Central Alps,SwitzerlandJ.Swiss Journal of Geosciences,111:317340.Cesare B,Satish-Kumar M,Cruciani G,et al.2008.Mi