《分配系数理论》PPT课件.ppt

高等地球化学,高等地球化学,马 东 升,南京大学 地球科学与工程学院 研究生课程,2011年 秋季版,第二章 分配系数理论和地球化学定量模型,参考教材：地球化学，科学出版社，2004，第11章,第三节,主要内容,一、元素的分布和分配二、分配系数理论三、岩浆作用的微量元素地球化学模型 四、主要造岩矿物和副矿物对岩浆微量元素含量的影响五、分配系数的选择,分布：元素在某一种地质单元或天体中（1）整体统计含量（对大单元特称丰度abundance）或（2）及其变化分配：元素在不同地质体、天体或物质相中的含量配置,一、元素的分布和分配,二次电子像,背散射电子像,阴极发光像,Si 面分布,K 面分布,Ca 面分布,矿物电子探针分析,化学元素的太阳系丰度,元素在太阳系中的分布,(%),1,2,3,（约1%）,镓在华南某黑云母花岗岩造岩矿物中的分配,（1）测定岩体中各矿物百分含量(A)（2）测定各矿物中某元素的百分含量（B）（3）求出岩体中某元素在某矿物中的量（AB）（4）求出岩体某元素的平均含量（C）（5）某矿物中某元素的量占岩体该元素总量的百分数（X）X(%)=(AB)/C,有关元素分配的几个基本规律和概念,1、戈尔德施密特三法则 Goldschmidts three rules（岩浆、热液过程）,(1)若两种离子具有相同的半径和电价，则它们进入特定晶格位置的能力也相等。,Zr:193 ppm(UC)0.80(+4,VI)Hf:5.3 ppm(UC)0.79(+4,VI),Nb:12 ppm(UC),0.72(+5,VI)Ta:0.9 ppm(UC),0.72(+5,VI),Zr:193 ppm(UC),0.80(+4,VI);Hf:5.3 ppm(UC),0.79(+4,VI).,Nb:12 ppm(UC),0.72(+5,VI);Ta:0.9 ppm(UC),0.72(+5,VI).,(2)若两种离子有相似的半径和相同的电价，则半径较小者优先从液相进入固相。,Isobaric T-X phase diagram at atmospheric pressure.After Bowen and Shairer(1932),Amer.J.Sci.5th Ser.,24,177-213.From Winter(2001)An Introduction to Igneous and Metamorphic Petrology.Prentice Hall.,铁橄榄石,镁橄榄石,Fe2:0.86(1010 m)Mg2+:0.80(1010 m),完全类质同相的镁铁橄榄石,具有相同电价和相似半径的离子相互置换,钾长石：Rb K；碳酸盐：Sr 2Ca 2；橄榄石：Ni 2Mg 2；铝硅酸盐：Ga 3Al 3+；锆石：Hf 4Zr 4；硅酸盐和石英：Ge 4Si 4；铌铁矿：Ta 5 Nb 5。,微量元素不单独形成矿物，而“隐蔽”在性质相似的高丰度元素的矿物中。,隐蔽（Camouflage）：,(3)若两种离子的半径相似，但电荷不同，则具较高电荷者优先从液相进入固相。相对于较小半径而言，离子半径差不超过15 电荷平衡：K+Si4+Sr2+(Ba2+)+Al3+（钾长石）2LiF MgO（云母）,捕获（capturing）：较低离子势（Z/r）的主元素被较高离 子势的离子置换（优先 早期置换）独居石：Th 4+Ce 3；钛的矿物：Nb 5 Ti 4；镁硅酸盐：Sc 3 Mg 2；长石：Sr 2或Ba 2 K允许（admission）：较高离子势的主元素被较低离子势的离 子置换（“勉强”晚期置换）Mg 2被Li置换，O 2 被F 置换,2、Ringwood的补充法则和负电性,具有较低负电性的离子优先进入固相,负电性 Electronegativity:The power of an atom in a molecule to attract electrons to itself.两种具有相同负电性的原子形成非极性共价键：EN 2。,1.金属元素的负电性低，非金属元素负电性高；2.在周期表中负电性从左到右增高；3.在周期表中同一族元素的负电性从上到下降低。,Pauling scale,元素负电性的周期变化,3、晶体场理论简介,对元素在成岩成矿过程中分布规律的认识和推断必须建立在热力学基础之上。Orgel（1952）应用晶体场理论（Bethe,1929）解释了第一过渡系列金属离子的水合热变化趋势。之后，该理论被广泛应用于过渡金属化学和地球化学之中。Mineralogical Applications of Crystal Field TheoryRoger G.Burns,Cambridge University Press,1970.晶体场理论的矿物学应用，R.G.伯恩斯，科学出版社，1977.,Hans Bethe 1906-2005,过渡元素：具有部分充填的d 或f 壳层的元素.,d电子层2种类型5个轨道,电子云的分布,八面体晶体场分裂参数o,简併,分裂,t2g轨道中的每一个电子使离子稳定了2/5 0eg中每一个电子使离子稳定性降低了3/5 0,Z轴,X或y轴,电子云,晶体场稳定能CFSE（Crystal Field Stabilization Energy):CFSE=(2/5)oN(t2g)+(3/5)oN(eg)N(t2g)和N(eg)分别为相应轨道中的电子数。,晶体场分裂,四面体和立方体配位,t2g轨道能量增高,eg轨道能量降低,八面体位置优先能OSPEOctahedral Site Perference Energy,c,t,o,o:c:t=1:(8/9):(4/9)OSPE:八面体场中的离子比它处于四面体场中时的能量降低，即稳定程度的增加。,CFSEt=(2/5)tN(t2g)(3/5)tN(eg),立方体晶体场分裂参数,四面体晶体场分裂参数,八面体晶体场分裂参数,八面体位置优先能OSPE,优 先 能,元素进入早期结晶矿物的顺序二价离子：Ni(Cr)(Cu)CoFeMnGa,Zn;三价离子：Cr(Mn)CoVTi FeSc,Ga():畸变效应影响,岩浆：四面体八面体 造岩矿物：八面体,二、分配系数理论,分配系数理论本身源于化学热力学。它将描述稀溶液的亨利定律应用于自然界中微量组份在岩浆、溶液和类质同像固溶体中的分配，从而达到热力学计算的简化。因此，特别适用于微量元素的地球化学研究，使许多重要地球化学元素对比值和同位素比值建立在坚实的热力学基础之上。在此基础上已经发展起来了许多我们目前赖以进行地球化学研究的理论和方法，如地质温度计、压力计、氧逸度计，同位素分配理论以及将要介绍的岩浆作用微量元素地球化学模型。,1、亨利定律和分配系数,拉乌尔定律和理想溶液 1887年法国化学家拉乌尔在研究溶液的蒸汽压与溶液组成关系时发现，某些溶液在恒温下的蒸汽压和它的组成之间呈线性关系，即如果以摩尔分数Xi表示溶液的组成，符合以下线性关系：Pi=Pi0 Xi 式中Pi为组份i在Xi浓度时的蒸汽分压，Pi0为纯组份i在相同温度下的饱和蒸汽压。该实验规律称为拉乌尔定律。性质十分相近的组份混合组成的溶液往往在整个浓度范围内都符合这一规律。这类溶液称为理想溶液。不符合拉乌尔定律的溶液则称为非理想溶液。,Francois Marie Raoult 1830一1901,理想溶液 在理想溶液中，具有相同粒子体积和晶格键力的组份混合不造成体系分子能态和体积的任何变化。在这种情况下，体系在混合过程既不吸热也不放热，因此，在理想状态下混合组份的活度（ai）等于它们的浓度：ai=Xi,当非理想组份越来越被稀释时，组份在溶剂中越来越分散，以至于最终其周围的离子或分子氛围变成相同状态。因此，在高度稀释时，组份的活度与其浓度之间呈正比例变化，在右图中用虚线表示，即有：ai=hi Xi式中hi是比例常数,又称亨利常数。,William Henry 1775-1836,亨利定律 Henrys law（1803）,亨利定律Henrys law,在组份i的稀溶液中（即摩尔分数Xi 0），物质i的蒸汽压Pi与稀溶液的Xi呈线性关系：Ki=Pi/Xi式中Ki为比例常数。该定律又称为“稀溶液定律”。,拉乌尔定律和亨利定律可用来解释许多地球物质的溶液行为。固溶体的许多常量组份混合行为符合拉乌尔定律，如橄榄石中的镁橄榄石和铁橄榄石组份；而亨利定律则常用来描述微量元素，如橄榄石中的镍。亨利定律常用于描述微量元素分配系数它是一种与亨利常数有关的可测定的量。,亨利定律的地球化学表述（简单）分配系数,在一定温度和压力条件下，当两个共存相A、B相平衡时，以相同形式均匀赋存于其中的微量组分i在这两相中的浓度比值为一常数：ki A/B=CiA/CiB(1)Ci A和Ci B分别为元素i在A、B相中的浓度，ki A/B为分配系数（Distribution Coefficient or Partition Coefficient），也称简单分配系数（Simple partition coefficient）、常规分配系数(Conventional distribution coefficient)或Nernst分配系数。,Walther HermannNernst,18641941,分配系数的热力学意义,当A、B两相平衡时，元素i在其间的化学势相等。即有：iA=iB。理想气体的化学势 i=i*(T)RTlnai，有 i*ARTlnaiA=i*BRTlnaiB(2)lnaiAlnaiB=(i*Ai*B)/RT去对数后，得：aiA/aiB=exp(i*B i*A)/RT=KiA/B(3)如果i浓度很低,则公式(1)ki A/B=Ci A/Ci B 与(3)式相等，浓度等于活度。,低浓度理想气体,复合分配系数Compound partition coefficient,简单分配系数易随体系的基本成份而变化。为了减小这一影响，应用复合分配系数（Compound partition coefficient）或标准化分配系数(Normalized distribution coefficient)可部分达到这一目的。为了与分配系数相区别，以下用D表示复合分配系数，其定义为：,下标cr：携带元素,即可被微量元素置换的常量元素；tr：置换常量元素的微量元素；下标s和l分别表示固、液相；ktrs/l和kcrs/l分别为有关微量元素和常量元素在固液相中的简单分配系数。,(4),复合分配系数的热力学意义,Ds/ltr-cr实质上是理想溶液中tr和cr之间的平衡置换常数例如，在KAlSi3O8Rb=RbAlSi3O8K的反应中，Rb（tr）置换钾长石中的K（cr）。该反应的平衡常数为产物的活度积除以反应物的活度积。如果将AlSi3O8简写为A，该反应的平衡置换常数为：K=a strAa lcr/a scrA a ltr(5)根据溶度积常数的定义，微量元素tr的纯晶体溶度积为：Ksp(trA)=(atrl/aAl)/atrs=atrlaAl（下标sp表示钾长石）。因此，对难溶电解质而言，平衡常数为：K=(acrl/acrAs)/(atrl/atrAs)=K sp(crA)/K sp(trA)(6)即cr和tr纯晶体的溶度积之比。,总体分配系数,bulk distribution coefficient,当固体由多相组成时，固液相共存平衡体系中元素i的分配系数称为总体分配系数（bulk distribution coefficient）式中Di为元素i在某单一矿物中的分配系数；上标a,b,n分别代表矿物a,b,n；上标l 表示与矿物平衡共存的液相；Xa，Xb，Xn分别为矿物a,b,n在岩石中的比例分数。,分配系数：kiA/B=CiA/CiB(1)i*ARTlnaiA=i*BRTlnaiB(2)R：气体常数1.987卡度摩尔aiA/aiB=exp(i*B i*A)/RT=KiA/B(3)如果i的浓度很低，（1）（3）,小结：,计算总体分配系数,已知与花岗质岩浆平衡的结晶岩石矿物组成为：石英(Q)20，钾长石(Or)10，斜长石(Pl)50，角闪石(Ho)20。各矿物的Rb,Sr,Ba分配系数分别如下表。求结晶岩石的总体分配系数。,结晶岩石的Rb的总体分配系数等于0.09。Sr的总体分配系数等于?Ba的总体分配系数等于?哪个是相容元素，哪个是不相容元素？,样品分离不纯对分配系数测定的影响,当岩浆结晶时，晶出的矿物与熔体之间的分离有时并不很完全，矿物之间甚至矿物中往往包捕获有熔体。这类似于样品人工分离不完全时的情况。设测定含基质的斑晶得到的为元素i的表观分配系数kis/l，它和实际分配系数kis/l的关系为：kis/l=kis/l(1-a)+a式中a 为被捕获熔体的量l与有关矿物固体的量s之比，即l/s。上式表明，当固相中的粒间熔体含量增高或样品分离不完全时，测得的表观分配系数总是趋于1。例如，当某一分配系数为2的矿物中含25%的熔体时，其表观分配系数降至1.75。,三、岩浆作用的微量元素 地球化学模型,（一）岩浆结晶作用,1.分离结晶作用 Fractional Crystallization 瑞利分馏定律的推导已知 k=Cs/Cl,令m为相的质量，x为摩尔数，则有：,k=,设体系质量为mo，元素总量即总摩尔数为xo,可得：,(1),积分：,定义 ml/mo=F(熔融度，即液相所占的质量百分数)：,将x转换为浓度C:C=x/m，F=ml/mo，去对数，两边同 乘以1/F：,，,，,(2);,(3);,(4)*,分离结晶作用公式,2.平衡结晶作用 Equilibrium Crytallization,k=Cs/Cl，Cl=Cs/k=,Co=xo/mo，Cl=xl/ml;,Cl=,平衡结晶作用公式,，Cl=Co/k(1-F)+F(5),或：Cl/Co=1/k(1 F)+F,Cs=kCo/k(1 F)+F,（二）部分熔融作用,实比过程（modal melting）：岩石按原有矿物的比例进行熔融的过程。在整个熔融进程中，固相矿物的比例不变，从而总体分配系数也不改变。这是一种为了计算方便而采取的简化处理，或因难以获得岩浆结晶历史中实际矿物组成变化的信息而采取的一种近似方法。非实比过程：指岩石的不按照源岩矿物组成的比例进行熔融的过程，即在熔融的不同阶段，固相矿物的组成是变化的。总体分配系数也相应改变。这显然更符合实际的情况。但是需要用其它手段确定熔融过程中矿物组成的变化。,（二）部分熔融作用,1、分离熔融作用 Fractional Melting 瑞利熔融过程的推导,概念：只有刚熔融的，无限小部分的熔体dml中的微量元素与已结晶出来的固相相平衡，即无限小的连续分离：,整理后积分,去对数后得：,ms/mo=1F（凝固度）,转化为浓度C（两边同乘mo/ms,即1/(1 F)）得到:,，,Cs=Cl k，,(6),分离熔融作用公式,实比过程（modal melting）:,非实比过程：,2、累积熔融作用 Incremented Fractional Melting,连续分离熔体在体系中的聚集和混合,（7）,解方程得：,相应的非实比方程为：,3、平衡（批次）熔融作用 Equilibrium or Batch Melting,与平衡结晶作用互为逆过程，因此可用同一个方程描述。通过对分配系数公式的整理和与上述同样方法定义熔融度，不难导出微量元素在完全平衡条件下结晶或熔融时分配的公式：Cl/Co=1/F+k(1F)在非实比过程的情况下，公式变为Cl/Co=1/k+F(1D)其中D代表当熔融度为F 时被熔融掉的矿物组合的总体分配系数。,共同特征：D(k)Co;D 1时，Cl Co D(k)=1 时，Cl=Co;D0时，Cl/Co=1/F 分离结晶作用:D(k)1时，F减小，Cl 迅速降低。当D(k)1时？,这两种过程的差别主要在于分配系数在分离结晶中呈指数形式，而平衡部分熔融时则在分母中。因此，结晶分离所产生的元素分异程度远比部分熔融大得多。这一差别尤其在强烈不相容元素上表现得特别显著。,在10、50和92.5的平衡部分熔融情况下，D=10的强相容元素在熔体中的含量分别降至其源岩的11、18%和60%；然而在相同熔融度的分离结晶作用中，残余熔体中D=10的强相容元素的浓度则仅为其母岩浆的109、2 和50。,Cl/Co=F D-1 分离结晶作用Cl/Co=1/F+D(1F)平衡部分熔融,相容元素,两类不相容元素：1、小半径高电荷的高场强元素（HFS，High Field Strength)REE，Th，U，Zr，Hf，Nb，Ta;2、低场强大离子亲石元素（LIL，LargeIon Lithophile)K，Rb，Cs，Ba，Pb2，Sr，Eu2,不相容元素,Cl/Co=1/F+k(1F),F1，C l/C o1,（三）饱和含水岩浆的热液分异,初 始岩 浆,热液H2O,硅酸盐熔体,结晶矿物,含水岩浆,K1=CH2O/Cmelt,K2=Csolid/Cmelt,在含水岩浆的结晶作用中，如某微量元素i基本上不进入固相，即其kis/l 1时，可近似地只考虑 i 仅在硅酸盐熔体和热液之间分配：k=C/(C)（下标 i 省略，下同）其中C为元素在热液中的浓度，（C）为元素在熔体相中的浓度。热液分异也有两种极端情况：1、围岩是可渗透的：与岩浆不平衡的热液分离作用2、围岩是不渗透的：与岩浆平衡的热液分异作用,1、可渗透围岩的不平衡热液分离作用,元素在m质量熔体中的量为m(C),其变化量为d(m(C)；它在热液相中相应的变化量为Cwd(C)，即有：d(m(C)Cwd(C)m：硅酸盐熔体的质量；（C）：岩浆中元素的浓度；C：热液中元素的浓度；w：岩浆中水的含量分数。根据k=C/(C)C k(C)和瑞利分馏条件可得：Ck（Co）F（kw-1）形式上与分离结晶作用的公式相似。,d(m(C)Cwd(C)C k(C)，,整理后得：,（ml)=(mo)F d(ml)=modF,两边积分,去对数,2、不渗透岩浆房的热液分异（平衡过程）,在平衡的热液分异过程中，初始岩浆中元素的减少量等于热液中的量：（Como)-(C l ml)=w(mo-ml)C根据Ck（C）的条件，推导可得：Ck(Co)/F+kw(1-F),Ck(Co)/F+kw(1-F),虚线：,实线：,平 衡,非 平 衡,元素在含水10岩浆热液分异中的富集度C(Co),k=C/(C),在k10时,因为kw=1,元素在热液相中的浓度不变；当k10时，元素在热液中的浓度随固结度增高而降低。在非平衡条件下浓度降低很快；在平衡条件下趋于kw=1。在k80%)出现差异。,Ck(Co)/F+kw(1-F),（四）岩浆结晶和岩石熔融的 平衡分析,实际自然过程能否符合建立数学模型时所做的平衡假设?1、结晶作用 一般规模的岩浆房分异结晶过程延续时间大致为0.51106年。在11001200C的温度下，元素在晶体中的典型扩散系数d为10-1110-13cm2/s。根据粒子在固体中的扩散距离S与扩散系数d和扩散时间t的近似关系：S(dt)1/2，因此，估算在岩浆房结晶时间内元素在固体矿物中的扩散距离为110 cm。因此，一旦矿物从熔体中结晶与熔体隔离，元素在固体矿物中就很难均匀分布而与残余岩浆达到平衡。,2、岩 浆 对 流？,流体能否在地球物理场中产生对流取，决于浮力是否能克服因流体粘度而产生的运动阻力。这两种力的比率称为瑞利数（Rayleigh Number）：,Ra是一个无量纲值。公式中的符号分别为：a热膨胀系数，DT流体垂直高度上的温度变化，d对流体系的高度，g重力加速度，r流体密度，K热传导率，h粘度。当浮力足以克服粘性阻力，而刚刚可能使流体产生对流时的瑞利数称为临界瑞利数Rac。必须在RaRac的条件下才能发生对流。,因此，如果知道有关物质的性质和相关环境的物理、几何参数，就可以预测是否可能形成对流。例如，若取地幔的相关参数分别为：a=310-5 K-1，g=10 m/s，r=4.66 Mg/m3，K=10-6 m2s-1，=1021 poise，d=2.9106m则可得Ra106。取决于不同的边界条件，估算的地幔临界瑞利数约为1,0005,000，大致为103数量级。因此，岩浆在地幔中显然能产生对流,根据火山喷出的熔岩体积和对深成岩体大小的直接测定，一个岩浆房的体积可达若干平方公里。在这种条件下，即使岩浆的粘度达到108泊，也可能发生对流。在熔体中（1200C）元素的扩散系数大约为10-7cm2/s*，所以岩浆熔体的成份基本上是均匀的。由此可见，岩浆结晶过程一般应遵循瑞利分馏定律，即符合分离结晶作用方程。根据上述单晶尺度的扩散距离，粘度较小的基性岩浆熔体可能与结晶矿物达到平衡。实验岩石学已证实了这一点，并且实际观察也表明，基性岩中矿物晶体的分带现象远比中酸性岩中更少见。,*元素在晶体中的典型扩散系数d为10-1110-13cm2/s,3、部分熔融：在部分熔融过程中，因为元素在熔体中的扩散较快，并且在岩浆房中可能存在对流，所以熔体往往是均质的。对于一个粒度为0.11 cm的固体矿物而言，如果元素在其中的扩散系数为1012cm2/s，40040,000年后固液相之间即可达到平衡。板块的平均扩张速度大致为2cm/年，所以地幔熔融物质上升及其在岩浆房中的停留时间也应有同样的速度。如果岩浆形成的深度为100200公里，不难算出它需要大约1107年才能到达地表。在这段时间内固体矿物和熔体应早已达到平衡。,4、花 岗 岩 类？,花岗质岩石的部分熔融发生在相对低温的地壳中，并且粘度更大。因此，花岗岩质与玄武质岩浆相比更难达到平衡。目前虽然对地壳深熔过程的动力学了解不多，但对变质岩的研究发现，许多高级变质作用就已达到岩相平衡。所以元素在深熔作用期间也应当达到平衡。因此，岩石部分熔融过程大多服从平衡熔融作用的公式。综上所述，尽管岩浆作用包括许多复杂的过程，但其中分离结晶和平衡部分熔融是两种在自然界中最常发生的过程。,片麻岩角闪石和黑云母中的V2O3含量关系,样品取自400平方英里的片麻岩中。这两种矿物中的Mn（低含量时）和Ti的分布也呈直线关系，表明岩体已达到平衡。,如果微量元素在岩体各部分样品的矿物对中有几乎一致的分配系数，则表明该岩体是在平衡条件下形成的,（五）分配系数和微量元素 地球化学模型的应用,1、部分熔融与分离结晶作用的甄别,在火成岩及其有关矿产成因研究中经常遇到的问题是：所研究的岩石或岩套是部分熔融形成的还是经分离结晶作用产生的？它们与源岩或母岩浆的关系如何？利用分离结晶和平衡部分熔融作用模型公式的性质有可能加以鉴别。,（1）利用亲岩浆元素和超亲岩浆元素,亲岩浆元素 M（magmatophile elememt)总体分配系数明显小于1（0.40.1），如 HREE，Zr，Hf。超亲岩浆元素 H（hypermagmatophile element)总体分配系数近于0（明显小于0.20.5）,如 Ta,Th,La,Ce（LREE）。,在分离结晶作用中(Cl=Co,lF(k-1)：CMl CMo,l/F,CHl CHo,l/F CHl/CMl CHo,l/CMo,l 常数 因此，在CHl/CMl对CHl作图的图解中，样品投影点呈水平分布。在平衡熔融作用中(Cl/Co,s=1/F+k(1F)：kH(1F)0,CHl CHo,s/F；kM F 0,CMl CMo,s/(kM+F)所以有：CHl/CMl=(CHo,s/F)/CMo,s/(kM+F)=(CHo,s/F)(kM+F)/CMo,s=(CHl kM+CHo,s)/CMo,s=CHl(kM/CMo,s)+CHo,s/CMo,s 因此，在CHl/CMl对CHl作图的图解中，样品投影点沿斜率为kM/CMo,s的直线分布。,邻近冰岛的雷克雅奈斯洋中脊样品呈批次部分熔融的趋势；而在冰岛则可能在浅部发生分离结晶作用,from Treuil&Joron,1975,在CHl/CMl对CHl图解中，分离结晶作用样品的投影点呈水平分布。,在CHl/CMl对CHl图解中，平衡部分熔融的样品投影点沿斜率为kM/CMo,s的直线分布。,在熔融度为00.5的范围内：分配系数明显大于1的相容元素，在平衡部分熔融作用中的含量变化比较稳定，变化小于强不相容元素。在熔融度为10.5的范围内：分配系数明显大于1的相容元素，在分离结晶作用中的含量变化明显大于强不相容元素（D0）。因此，当源岩先部分熔融，然后再经分离结晶固结成岩后，在所形成的一系列岩石中，这一对微量元素的含量变化将反映出这一演化过程。,（2）利用分配行为相反的微量元素,在熔融度为00.5的范围，分配系数明显大于1的相容元素在平衡部分熔融作用中的含量变化比较稳定，变化小于强不相容元素；在熔融度为10.5的范围内，分配系数明显大于1的相容元素在分离结晶作用中的含量变化明显大于强不相容元素（D0）。,实 例 考 察,设一种岩石含：Ce（铈）10 ppm(D=0，强不相容元素)Sr（锶）278 ppm(D=4，相容元素)开始发生40的平衡部分熔融；然后所形成的岩浆发生分离结晶作用，直至完全结晶；考察不相容元素和相容元素含量的变化。,40的平衡部分熔融过程,C l/C o=1/F+D(1F),分离结晶作用,Cl/Co=F D-1,当某一源岩先部分熔融后，再通过分离结晶固结成岩所形成的一系列岩石的过程，将会在这一对微量元素的含量变化上反映出来。相容元素（D=4）和不相容元素（D=0）在两套成因岩石中有不同的分布。图中百分数为熔融度F。由于Ce的分配系数为零，它在部分熔融和分离结晶作用中的浓度按Cl/Co1/F的趋势变化。,（据Hanson，1978修改）,F,F,四、主要造岩矿物和副矿物对 岩浆微量元素含量的影响,（一）花岗质岩浆,钾长石：钾长石有高的Sr、Eu分配系数和低REE（Eu除外）分配系数（DSr3.87，DEu=1.13，DREE=0.00n)，与斜长石相似。但是它的Ba分配系数(DBa6.12)比斜长石高得多。钾长石会造成熔体有很高的Sr/Ba比值，并且使熔体中的Rb/Sr增加，但增加的速率不像斜长石那样快。钾长石的分离使熔体的K/Rb比值明显降低。,斜长石：锶是相容元素（D4.4），铷是不相容元素（D=0.04），斜长石的分离使熔体中Sr贫化，Rb增高。锶和铷的分配系数在斜长石中相差两个数量级，所以斜长石的不断分离使熔体Rb/Sr比值急剧增大。斜长石的高Eu正异常（D=2.15）造成熔体的Eu负异常。斜长石Eu异常强度随氧逸度和温度降低而增高。斜长石的K分配系数（D=0.1)明显高于Rb，但因其绝对值都很小，所以斜长石分离对熔体K/Rb比值影响不大。除了Sr、Rb和Eu之外，斜长石对其它元素在熔体中的变化都影响都不大。如果斜长石中有较高的K呈固溶体存在，Rb和Ba的分配系数会明显增高（钾长石：DRb0.66，DBa6.12，Philpotts et al.,1970；歪长石：DRb0.45，DBa5.04，Sun et al.，1976）。,黑云母：它具有比较平缓的稀土分配系数变化，因此，对熔体的REE浓度影响不大。黑云母与钾长石是花岗岩类中主要的两种富钾矿物。它们的存在能对熔体中的钾含量起缓冲作用，而使熔体保持稳定的钾含量。黑云母的分离将使熔体的Rb/Sr比值降低，Sr/Ba比值增加。DRb=3.26，DSr=0.12，DBa=6.36,角闪石：对熔体Rb、Sr、Ba和K含量影响不大。K和Rb在角闪石中的分配系数都很小（DK=0.081，DRb=0.014），所以它的分离虽然能使熔体的K/Rb比值有所增加，但幅度不大。角闪石的结晶主要对熔体的稀土分配有影响。在中、酸性岩浆中，REE对角闪石是相容元素，并且在GdEr之间有最高的分配系数。因此，角闪石的分离使熔体内的中稀土降低。重稀土也有较小幅度降低，并产生铕正异常。当岩石含较高的角闪石时（2030），对岩石的稀土元素含量有重要影响。,辉石：对熔体微量元素的影响类似于角闪石。但因为其中稀土元素的分配系数相对较低，尤其是紫苏辉石，所以对REE的影响比角闪石更小。紫苏辉石的分离能使熔体略富集轻稀土和造成铕正异常；单斜辉石将造成熔体的重稀土亏损，较小程度的重稀土亏损，和弱的铕正异常。,石榴子石：具有极低的Rb、Sr、Ba、K和轻稀土分配系数，以及高的重稀土分配系数。当它以显著的量与熔体分离时，造成熔体明显的重稀土亏损，并产生铕正异常。REE在石榴子石中的分配系数变化极大。因此，当石榴子石作为与熔体平衡的残留矿物存在于岩浆源区时，火成岩中的重稀土相对于轻稀土会产生明显的亏损。在玄武质岩浆中，因为Lu在石榴子石中的分配系数是La的三个数量级（DLu=39.755，DLa=0.390），所以重稀土相对于轻稀土会产生极端的亏损。,其它副矿物：即使它们的含量较低，但因为在中酸性岩浆中都有较高的REE分配系数，如锆石、磷灰石、榍石等，因此会对稀土元素分配有较大影响。锆石产生重稀土亏损，磷灰石主要形成中稀土亏损，两者都产生铕正异常；榍石也形成中稀土亏损；而独居石和褐帘石则造成轻稀土亏损,（二）玄武质岩浆,在玄武质岩浆中，除了2、3价过渡元素和石榴子石中的重稀土之外，大多数微量元素具有较低的矿物熔体分配系数。因为单斜辉石、石榴子石和斜长石等对不相容元素具有最高的分配系数，所以它们是控制基性岩熔融和结晶期间不相容元素分配的主要矿物。,角闪石：因其在低压和高温（1100）下不稳定，故在玄武岩中较少出现。橄榄石：是地幔中最常见矿物。它的不相容元素分配系数很低，对不相容元素只产生较小分馏；但它是控制相容过渡元素分异的主要矿物。尖晶石：岩石中含量通常很低，它对微量元素浓度的影响不大。石榴子石：因为Lu在石榴子石中的分配系数是La的三个数量级（DLu=39.755，DLa=0.390），所以石榴子石的分离将使重稀土相对于轻稀土产生极端亏损。,五、分配系数的选择,尽量选择那些在与被研究体系的性质，尤其是化学成份最接近条件下测定的分配系数。经验表明，无论是直接测定天然样品或经实验得到的分配系数往往偏高。因此一般可选择较低的值。为了避免污染，最理想的是使用原地微量分析技术。如用离子探针获得的数据显然优于全岩或全矿物分析。尽量选用不掺有与研究对象无关的其它杂质的实验数据。,4.在研究多元素时（如一套LIL、HFS或稀土元素），不同元素分配系数之间精确的相对关系比其绝对值更重要。因此应尽量选择那些在相同实验条件下一次测得的，并同时测定较多元素分配系数的数据。5.研究基性岩浆时，应首选在12001400 条件下用天然玄武质熔体进行实验测定获得的分配系数。6.尽量排除那些分配系数强烈受成份或温度影响的元素，如Ni。,Thank you for your attention!,一个花岗岩体主要由石英、钾长石和斜长石组成。这三种矿物在岩石中的比例大致各占1/3。已测得该岩体的Rb、Sr、Ba含量分别为30ppm、500ppm和600ppm。设该岩体先发生部分熔融,待其全部熔为岩浆后再发生分离结晶作用。在这些过程中岩石的矿物比例保持不变,即为实比过程。已知各元素的分配系数如下表所列。试计算熔融度分别为10%、20%、40%和80%时这两种过程所产生的熔体中Rb、Sr和Ba的含量,并以表格的形式列出。,铷、锶和钡在矿物中的分配系数,练习题 1.,对一种含橄榄单斜辉石岩包体的玄武岩成因进行研究。测定从出玄武岩和橄榄单斜辉石岩包体分离出来的单矿物(单斜辉石和橄榄石)中若干微量元素与钾的含量分别如下:,(上表中CeN和YbN分别为铈和镱含量球粒陨石标准化后的数据，其余元素含量单位为ppm)。计算中涉及到的玄武岩和包体矿物中元素的分配系数如下：,已知橄榄单斜辉石岩包体含80%橄榄石和20%辉石。求:(1)上述数据是否表明 该包体是来自玄武 岩源区的捕虏体?(2)在对问题(1)求解基 础上估算这种包体 形成时的熔融度F(假定源区稀土元素 丰度是球粒陨石的 12倍)。(3)估算源区的钾,铷,锶 含量。,练习题 2.,