《地质流体稳定同位素示踪》由会员分享,可在线阅读,更多相关《地质流体稳定同位素示踪(134页珍藏版)》请在金锄头文库上搜索。
1、地质流体的稳定同位素示踪第一篇 总 论查阅资料重要periodical 1. GCA(Geochimica et Cosmochimica Acta) 2. Chemical geology 3. EPSL (Earth and Planetary Science Letters) .Applied Geochemistry 5. Economic Geology 6. Journal of Geochemical Exploration 7. Lithos 8. Metamorphic Geology 9. Science 10. Nature 序一、同位素及其分类一、同位素及其分类1 1、
2、同位素、同位素(isotope)质子数相同,中子数不同的元素。质子数相同,中子数不同的元素。2 2、放射性同位素、放射性同位素(radioactive isotope) 放射性同位素放射性同位素: :凡能自发地放出粒子而衰变为另一种同凡能自发地放出粒子而衰变为另一种同位素者称为放射性同位素。位素者称为放射性同位素。衰变类型衰变类型(1) (负电子)衰变 一个中子转变成一个质子和一个电子 40 19K4020Ca+ +Q(2)正电子衰变 质子转变成中子 189F188O+ +Q-(3) 电子俘获衰变 40 19K + 4018Ar+Q(4) 衰变 23892U23490Th+42He+Q(5)
3、核裂变 23892U不同产物的核,典型的Zr(40)和Cs(55)伴随其他粒子和大量能量3 3、稳定同位素、稳定同位素(stable isotope)及放射性成及放射性成因同位素因同位素(radiogenic isotope)稳定同位素:无可测放射性的同位素。稳定同位素:无可测放射性的同位素。 其中一部分是放射性同位素衰变的最终稳定产物,其中一部分是放射性同位素衰变的最终稳定产物,称之为称之为放射成因同位素放射成因同位素。 另一部分是另一部分是天然的稳定同位素天然的稳定同位素,即自核合成以来就,即自核合成以来就保持稳定的同位素。保持稳定的同位素。238U234Th234Pa234U230Th2
4、26Ra222Rn218Po214Pb214Bi214Tl (214Po)210Pb210Bi210Po206Pb 二、二、同位素在地质学上的应用同位素在地质学上的应用1 1、放射性成因同位素放射性成因同位素 地质事件的精确年龄;地质事件的精确年龄;地质过程的天然示踪剂;地质过程的天然示踪剂;地地球球主主要要储储库库P Pb b同同位位素素组组成成图图1 1 海洋玄武岩同位素组成变化范海洋玄武岩同位素组成变化范围围地质演化的速率和途径地质演化的速率和途径-利用不同封闭温利用不同封闭温度度 K-Ar:角闪石角闪石520 ;白云母;白云母350;黑云;黑云母母300; 斜长石斜长石250;钾长石;
5、钾长石160 Rb-Sr:斜长石:斜长石600;白云母;白云母550;黑;黑云母云母350;2 2、稳定同位素、稳定同位素 地球化学示踪地球化学示踪(流体和物质流体和物质) 地质温度计地质温度计三、 地质流体的同位素示踪意义 本课程重点应用范围:1、成矿流体的来源示踪与矿床成因 和找矿预测 2、俯冲带流体的来源示踪与壳/幔 演化 3、地表流体的来源示踪与环境地球 化学第一节 稳定同位素基本知识稳定同位素基本知识第一章、稳定同位素基础第一章、稳定同位素基础一、同位素丰度一、同位素丰度绝对丰度:指某一同位素在所有各种稳定同位绝对丰度:指某一同位素在所有各种稳定同位素总量中的相对份额。通常取与素总量
6、中的相对份额。通常取与1H(1012) 或或 28i(106)的比值。的比值。相对丰度:指同一元素各同位素的相对含量。相对丰度:指同一元素各同位素的相对含量。1H(99.985%),2D(0.015%)16O(99.762%),17O(0.038%),18O(0.200%)二、二、 同位素效应同位素效应 由同位素质量差所引起的物理和化学性质上的差异。第二节、 同位素分析结果的表达和标准正是由于同位素效应,造成了地质过程中,两相或多相共存时,不同正是由于同位素效应,造成了地质过程中,两相或多相共存时,不同相中同位素比值差异相中同位素比值差异使同位素示踪成为可能。使同位素示踪成为可能。一、实测一、
7、实测值值 A =1000 (RSa RSt)/RSt (单位:位:) SMOW (RSt):18O/16O: 2005.210-6;Rock( RSa):18O/16O:202010-6则该岩石则该岩石18O = 1000 ( 2020 10-6 2005.210-6 ) 2005.210-6 =+9.4二、 同位素标准要求:组成均一,性质稳定;数量相当大,以便长期使用;化学置备和要求:组成均一,性质稳定;数量相当大,以便长期使用;化学置备和同位素测量的手续简便;大致为天然同位素比值变化范围的中值,以便同位素测量的手续简便;大致为天然同位素比值变化范围的中值,以便用于绝大多数样品的测定;可作为
8、世界范围的零点。用于绝大多数样品的测定;可作为世界范围的零点。H和O同位素标准:SMOW(standard mean ocean water) O同位素有时用PDB18OSMOW=1.0309118OPDB+30.91C同位素标准:Pee Dee Belemite (PDB)南卡罗林纳州S同位素标准:亚历桑那州Canyon Diablo Troilite铁陨石的陨 硫铁(CDT)N同位素标准:空气同位素之比 缩写符号 标准样 2D/1H SMOW 标准平均大洋水 18O/16O SMOW 标准平均大洋水18O/16OPDB 美国南卡罗林纳州白晋系皮狄组的美洲似箭石 13C/12C PDB美国南
9、卡罗林纳州白晋系皮狄组的美洲似箭石34S/32S CDT 美国亚历桑那州卡扬迪阿布洛铁陨石中的陨硫铁 三、自然界一般物质的 同位素值范围第三节、第三节、 同位素分馏同位素分馏(isotope fractionation) 指在一系统中,某元素的同位素以不同的比值分配到两指在一系统中,某元素的同位素以不同的比值分配到两种物质或两相中的现象,是同位素效应的表现。种物质或两相中的现象,是同位素效应的表现。热力学平衡分馏热力学平衡分馏 热力学平衡分馏可以有很多过程热力学平衡分馏可以有很多过程,但都达到平衡。但都达到平衡。例如:化学反应、扩散交换例如:化学反应、扩散交换动力学非平衡分馏动力学非平衡分馏A
10、.交换时未达到平衡交换时未达到平衡后期平衡被破坏后期平衡被破坏非质量相关分馏非质量相关分馏同位素质量差一般与分馏效应成正比。同位素质量差一般与分馏效应成正比。碳质球粒陨石的白色包体碳质球粒陨石的白色包体不满足不满足质量相关分馏质量相关分馏第四节、分馏系数第四节、分馏系数由由定义得定义得:A=(RA/Rst1) 1000,即RA= (10-3A +1) Rst, 同理RB= (10-3B+1) Rst代入定义式定义式A-B= RA/ RB= (10-3A +1)/ (10-3B+1)(1)或或=(A +103)/( B +103).(2)公式(1)两边取自然对数 A-B= (10-3A +1)
11、(10-3B +1),右边泰勒展开得 A-B 10-3A 10-3B 10-3(A B )或 103 A-B A B .(3)或 : 103 A-B, ( 富集系数富集系数= A B )由(2)得或 103 (A-B 1) (4)由上得:由上得: 103 A-B A B 103 ( A-B 1) 第二章、第二章、 稳定同位素分析方法稳定同位素分析方法(略)(略)两个步骤两个步骤:(1)样品制备样品制备(2)质谱测定质谱测定1. 样品制备样品制备具体要求具体要求1.1 氢同位素样品氢同位素样品 (1) 水,铀或锌还原法水,铀或锌还原法 (2) 含羟基矿物含羟基矿物 (3) 气液包裹体气液包裹体1
12、.2 氧同位素样品(1)水(2)碳酸盐磷酸法(3)硅酸盐和氧化物1.3 碳同位素样品(1)有机化合物(2)碳酸盐1.4 硫同位素样品(1)硫化物直接氧化法,SF6法(2) 自然硫和硫酸盐三酸还原法碳还原法热分解法第三章、第三章、 同位素热力学平衡分馏同位素热力学平衡分馏同位素交换反应:同位素交换反应:分馏系数理论计算:分馏系数理论计算:第一节、同位素分馏系数获取第一节、同位素分馏系数获取(校准)方法校准)方法一、理论计算一、理论计算二、实验测定二、实验测定、双向交换、双向交换()()、外推法、外推法(公式,(公式,?)b= (i - f) / (i - e) b (i - e) = (i- f
13、 ),两边乘两边乘b得:得:i - e = (i- f )b,i = e+ (i- f )b或书上:或书上: i = e (f i )b三、经验估计三、经验估计条件:条件:、假定平衡、假定平衡、已知温度(一系列温度,线性拟合)、已知温度(一系列温度,线性拟合)第二节、影响同位素平衡分馏的第二节、影响同位素平衡分馏的因素因素一、温度一、温度103ln=A106/T2+B103/T+C二、压力二、压力三、组分三、组分 物质的同位素分馏行为最终决定于键性(氧化态、离子电荷、原子量、同位素原子及与之结合的元素的电子排布等),所有这些都与振动频率及其变化率有关。一般地,重同位素倾向于富集在一般地,重同位
14、素倾向于富集在键性强键性强的化合物中。的化合物中。 H2O(l)+HDS(g)HDO(l)+H2S(g) H-O键比键比H-S键强,所以键强,所以D富集在水分子富集在水分子中。中。电位高、原子量低电位高、原子量低的阳离子优先结合的阳离子优先结合18O。例如对于例如对于Si4+O键到键到Al3+O到到Fe2+O键,阳离子的电位逐渐降低,原子量逐键,阳离子的电位逐渐降低,原子量逐渐升高,结果是渐升高,结果是18O富集程度逐渐降低。富集程度逐渐降低。因此,在氧化物和硅酸盐矿物中,石英因此,在氧化物和硅酸盐矿物中,石英(SiO2)最富最富18O,而磁铁矿,而磁铁矿(Fe3O4)较贫较贫18O。当铝硅酸
15、盐格架中。当铝硅酸盐格架中Al取代取代Si时,伴时,伴随有随有18O亏损,结果相对于石英,碱性长亏损,结果相对于石英,碱性长石石(K,Na,Rb)AlSi3O8较贫较贫18O,钙长石,钙长石Ca2Al2Si2O8更贫更贫18O。四、物质结构在同一物质的三种不同物态之间,重同在同一物质的三种不同物态之间,重同位素倾向于富集在较位素倾向于富集在较紧密紧密或或有序度高有序度高的的物态中。物态中。冰冰 水水 水蒸汽水蒸汽文石文石 方解石方解石金刚石金刚石 石墨石墨第四章、动力学分馏第四章、动力学分馏第一节、同位素交换机理第一节、同位素交换机理一、单向化学反应 C+16O18OC16O18O C+16O
16、2C16O2二、扩散现象三、溶解再结晶第二节、地质过程中的动力学同第二节、地质过程中的动力学同位素分馏位素分馏瑞利分馏在自然界存在一种特殊的体系,在一定在自然界存在一种特殊的体系,在一定的物理化学条件下发生物相分离。分离的物理化学条件下发生物相分离。分离前不同物相之间保持着热力学平衡并处前不同物相之间保持着热力学平衡并处于封闭状态,但分离后一相物质不断离于封闭状态,但分离后一相物质不断离开体系,不再与另一相保持平衡。这种开体系,不再与另一相保持平衡。这种在开放体系中进行的过程称之为在开放体系中进行的过程称之为瑞利过瑞利过程程,在瑞利过程中发生的同位素分馏称,在瑞利过程中发生的同位素分馏称之为之
17、为瑞利分馏瑞利分馏。在物相分离前,相A和相B处于同位素平衡状态,分离后相A不断离开体系,B为残留体系,f为残留分数,有:RB/R0=f (A-B1).(1)换算成常用的表示由公式:A-BR/R(1+ 10-3A)/(1+ 10-3 B) 得:R/R (1+ 10-3)/(1+ 10-3 )代入(),两边取自然对数得:(1+10-3B)/(1+ 10-30) (A-B1)f 右边展开: (1+10-3B)/(1+ 10-30) 10-3B 10-3010-3B 10-30 (A-B1)f整理得:整理得:B B 0 0 + + 10103 3 (A-B1)f残留分数为f时残留相的的计算公式分离相的
18、(平均值)计算公式:由质量平衡公式得:B fA( f)A ( B f) ( f)B的计算公式代入得:A ( R0f (A-B1 f) ( f)( f ) ( f)A ( f ) ( f)换算成常用的表示:利用 R/R (1+10-3 A)/(1+ 10-3)(1+10-3 A)/(1+ 10-3) ( f ) ( f)左边 10-3 A 10-3A A 103 ( f ) ( f)另一种形式(书上):由(1+10-3 A)/(1+ 10-3)= ( f ) ( f)整理得:1+10-3 A= (1+ 10-3)( f ) ( f)A A= (1000+ = (1000+ )()( f ) (
19、f) 10001000例子书上思考:地质应用思考:地质应用第五章、第五章、 同位素地质测温同位素地质测温第一节、第一节、方法原理方法原理对于两种矿物或一种矿物与水间,分馏满足方程:103ln=A106/T2+B103/T+C,(,(同位素同位素分馏方程)分馏方程)一、单个样品一、单个样品测定两种矿物M1和M2的1和2,计得1-2= 12103ln 1-2= 12 1-2= 12 A106/T2+B103/T+C.(1)方程(1)中、为常数,分馏方程表中查得,将实测的1和2代入解出,单位为, 273.15得二、二、 的加和性与未知分馏方程的获得的加和性与未知分馏方程的获得 三、多组样品或多组矿物
20、对(下节三、多组样品或多组矿物对(下节 同位素平衡检查)同位素平衡检查)第二节、 同位素平衡检查一、共生顺序判别法(一、共生顺序判别法()例如:氧同位素有:例如:氧同位素有:石英石英方解石方解石 碱性长石碱性长石 蓝晶石蓝晶石多硅多硅白云母白云母 钙长石钙长石 白云母白云母 .硫同位素有:硫同位素有:硫酸盐硫酸盐辉钼矿辉钼矿黄铁矿黄铁矿闪锌矿闪锌矿黄铜矿黄铜矿斑铜矿斑铜矿方铅矿方铅矿辉银矿辉银矿二、等温线法二、等温线法 (以以P52分馏方程举例)分馏方程举例)主要基于1-2= A 1-2 106/T2+C 1-2 (忽略B103/T项)选择一种矿物(一般为富重同位素矿物)同位素矿物)为参考矿物
21、,以- - -为纵坐标,以A - -为横坐标矿物有:-= A -106/T2(斜率为106/T2 )矿物有:-= A -106/T2(斜率为106/T2 )矿物有:-= A -106/T2(斜率为106/T2 )假如参考矿物、矿物、间在温度时同位素交同位素交换平平衡衡,则在图上落在斜率为106/T2 的同一条直线上。其实很多分馏方程为1-2= A 1-2 106/T2形式(忽略B103/T+C项,见、表)此时,以以-为纵坐标,以A -为横坐标对于多种于多种矿物物对可根据可根据拟合的直合的直线斜率由斜率由106/T2 计算算出出三、三、-图解法解法这是最常用的一种方法。对于两种矿物Y和X(Y更富
22、重同位素),当两者处于热力学平衡分馏时,有:yx103ln y- x ,即:yx 103ln y- x 在以y和和x为纵、横坐、横坐标的的图解上,解上,样品将落在斜率品将落在斜率为的一条直的一条直线上,上, y轴上的截距上的截距为103ln y- x ,同一地质体不同点的样品、i的矿物Y和X的测试值将将落落在斜率在斜率为、截距、截距为103ln y- x的直的直线上。上。根据根据矿物Y和X的分馏方程: 103ln y- x 可算出温度(截距可算出温度(截距为103ln y- x可直接从图上读出)这就是多个就是多个样品的温度品的温度计算方法算方法T=50T=550T=600T=+T=300T=4
23、00T=600四、图解垂线判别法选择一种物相(例如,热液矿床可选水)为参考物相,以103ln- 为纵坐标、以106T2为横坐标作图,当矿物和处于平衡分馏时,即相同( 106T2相同),具相同横坐标,此时矿物和的投影点连线为垂直线。 第三节、 同位素地质温度计同位素地质温度计一、氧同位素地质温度计一、氧同位素地质温度计、外部测温法只测定一种矿物同位素,另一相(通常为液相)采用假定值),根据已知矿物水间的分馏方程,计算温度。最常用于古温度测定。都建立了经验公式、内部测温法、矿物水、矿物矿物前提:()矿物对间达到分馏平衡;()矿物对间足够大;足够大;()矿物对间103ln随温度变化明显;()有已知的
24、分馏方程;()无后期同位素退化、单矿物测温法主要利用大气降水方程,测定粘土矿物 18 O、D。设某粘土矿物与水间氢、氧同位素间有分馏方程:103ln10/T2+BH.(1)103lnOO10/T2+BO .(2)大气降水: D 18O (3)103ln= DM D,即D= DM 103ln DM 10/T2+BH .()同理 18O = 18 O 103ln 18 O O10/T2+BO .()将()、()代入(),整理得:(O ) 10/T2 18 DM (BH BO) .()方程()为单矿物测温方程二、硫同位素地质温度计二、硫同位素地质温度计硫化物间,同氧同位素内部测温法。三、三、 碳同位
25、素地质温度计碳同位素地质温度计方解石石墨方解石石墨练习、单样品矿物对温度计算、单样品多矿物温度图解计算(等温线法)、多样品矿物对温度图解计算( 图解)第二篇同位素各论第一章 氢和氧同位素地球化学1H(99.985%),2D(0.015%)16O(99.762%),17O(0.038%),18O(0.200%)第一节、天然水一、物态转变过程中的氢氧同位素分馏蒸发冰冻二、海洋水的氢氧同位素组成短周期变化 1、蒸发(盐度变化) 2、冻结(实际上主要因引起盐度变化) 3、混合(其他洋水体加入)长周期变化 1、水岩相互作用(大洋沉积物18O很高、 D低,矿物沉积导致洋水18O降低, D升高)2、 冰川融
26、化二、大气降水二、大气降水纬度效应,即随着纬度升高纬度效应,即随着纬度升高(即即年平均气年平均气温降低温降低), D和和 18O值下降;值下降;大陆效应,即从海岸向大陆内部,大陆效应,即从海岸向大陆内部, D和和 18O值下降;值下降;季节效应,即夏季温度较高,大气降水季节效应,即夏季温度较高,大气降水相对较相对较“重重”,富集,富集18O和和D;冬季反之;冬季反之;高度效应,即随高度增加,高度效应,即随高度增加, D和和 18O值值下降。下降。降雨降雨Rayleigh分馏分馏大气降水线大气降水线中国现代: D=7.9 18O+8.2岩浆侵入浅部地岩浆侵入浅部地壳壳加热围岩和水加热围岩和水导致
27、水岩相互导致水岩相互作用;作用;中性、中性、“氯化物氯化物”地热水地热水H同位素同位素组成与当地大气降组成与当地大气降水类似,但水类似,但 18O值值升高;升高;酸性富硫的地热酸性富硫的地热水水H和和O同位素组同位素组成均不同于成均不同于当地大当地大气降水气降水.三、地热水三、地热水四、其他水初生水(来自地幔未与水圈相遇): D=-6020, 18O =61岩浆水:岩浆水: D=-80 -50, 18O = 5 7变质水:变质水: D=-40 -100, 18O = 5 25第二节、水岩交换作用第二节、水岩交换作用封闭体系封闭体系 Wi- Wf移项:W Cw( Wi- Wf)= R Cr( r
28、f- ri) Wf用用 Rf 代入,变形得:代入,变形得: 18ORf= 18ORi+W/R( 18Owi+ )/(1+W/R) Rf用用 Wf + 代入,变形得:代入,变形得: 18Owf= 18ORi+W/R( 18Owi) /(1+W/R)开放体系开放体系氧氧同同位位素素和和矿矿床床Skaergaard 侵入体侵入体At high temperatures(temperatures of the interior of the Earth or magmatic temperatures), oxygen isotope ratios are minimally affected by
29、chemical processes and can be used as tracers much as radiogenic isotope ratios are.These generalizations lead to an axiom: igneous rocks whose oxygen isotopic compositions show significant variations from the primordial value (+6) must either have been affected by low temperature processes, or must
30、 contain a component that was at one time at the surface of the Earth (Taylor and Sheppard, 1986).火成岩火成岩绝大多数火成岩的绝大多数火成岩的 18O变化范围为变化范围为515, D范围为范围为-40-100。一般来。一般来说,说, 18O值随值随SiO2含量增加而增加。含量增加而增加。年年青青新新鲜鲜玄玄武武岩岩氧氧同同位位素素组组成成MORB:5.56.5OIB:57.5Oxygen isotope ratios in olivines and clinopyroxenesfrom mantl
31、e peridotite xenoliths. Data from Mattey et al. (1994).岛弧火山岩蛇蛇绿绿岩岩洋中脊玄武岩海水相互作用洋中脊玄武岩海水相互作用高温氧同位素交换高温氧同位素交换低温氧同位素交换低温氧同位素交换花岗岩高18O花岗岩( 10)1、源区 2、围岩混染 3、低温水/岩反应低18O花岗岩(400C时,含S物种主要是H2S和SO2,此时, 34S热液热液= 34SH2S+ 34SH2SSO2R/(1+R)R为为SO2/ H2S摩尔比关键要确定R假如成矿热液来自岩浆系统,此时, 34S热液热液与与 34S熔体熔体的关系是:的关系是: 34S热液热液= 34
32、S熔体熔体-34SH2SHS- + 34SH2SSO2 R/(1+R)由此可见,也不能简单用由此可见,也不能简单用 34S热液热液与岩体与岩体 34S类类比,进而判别比,进而判别成矿物质和热液来自岩体岩体图图8.3,但当磁铁矿,但当磁铁矿-磁黄铁矿磁黄铁矿-石英石英-橄榄石共橄榄石共生时,此时生时,此时 34S热液热液 34S熔体熔体 34SH2S二、中低温热液中硫同位素当T+15,地幔:地幔:+5.6. 岛弧岩浆:沉积岛弧岩浆:沉积物物+地幔,地幔, 18O高于地幔高于地幔. 但也可能是地壳混染的结但也可能是地壳混染的结果果.87Sr/86Sr结合结合O和放射成因同位素,可以区分这两种过程和
33、放射成因同位素,可以区分这两种过程. 事事实上,硅酸盐岩石,包括沉积物,岩浆,和橄榄岩实上,硅酸盐岩石,包括沉积物,岩浆,和橄榄岩O含量类似含量类似(close to 50%), 但放射成因同位素但放射成因同位素Sr, Nd, Hf, Pb, etc.变化较大变化较大. 这些元素这些元素(except Os)为为不相容元素,岩浆中的含量高于地幔不相容元素,岩浆中的含量高于地幔. 许多岛弧岩许多岛弧岩浆浆Sr含量比大陆地壳高,但地幔含量比大陆地壳高,但地幔Sr含量比沉积物低含量比沉积物低. 意味着沉积物意味着沉积物+地幔混合线不同于岩浆地幔混合线不同于岩浆+地壳地壳.稳定同位素判别源区混合稳定同
34、位素判别源区混合vs.地壳混染地壳混染沉积物沉积物+地幔混合,地幔混合,r = SrM/OM/SrS/OS1 18O vs. 87Sr/86Sr.设未分馏体系重、轻同位素量分别为设未分馏体系重、轻同位素量分别为n2和和 n, A、分别为生成和残留相,、分别为生成和残留相,n2、n1分别为分别为t时刻残留相的重、轻同位素量,有时刻残留相的重、轻同位素量,有d n2、 d n转移到相转移到相,此时此时 d n2 d n,( n2 d n2)( n1 d n)A-B /d n2 d n/( n2 d n2)( n1 d n)n2 d n2 n2 , n1 d n n1 ,上式写为:,上式写为:A-B
35、 d n2 d n/ n2 n1,即,即A-B (d n n1) d n2 n2积分得(积分得( n n1) A-B n2 n2,两边同除,两边同除n n1得:得:( n n1) ( n2 n1 ) ( n2 n ) ,因为重、轻同位素来说,轻同位素含量因为重、轻同位素来说,轻同位素含量重同位素含量重同位素含量 例如:1H(99.985%),2D(0.015%); 16O(99.762%),17O(0.038%),18O(0.200%)n n1 ( n n)( n1 n)f,即,即f ( n2 n1 ) ( n2 n ) ( n2 n1 )为为t时刻残留相同位素比值时刻残留相同位素比值RB, ( n2 n ) 为体系初始同位为体系初始同位素比值素比值R0RB/R0=f (A-B1)