欢迎访问《高校地质学报》官方网站,今天是
分享到:

J4 ›› 2010, Vol. 16 ›› Issue (2): 213-217.

• 矿床岩石地球化学及年代学专栏 • 上一篇    下一篇

磷灰石中氦扩散参数的确定方法及求解过程

张 彦, 陈 文, 刘新宇   

  1. 中国地质科学院 地质研究所氩-氩年代学实验室,北京 100037
  • 收稿日期:2009-06-23 修回日期:2010-03-15 出版日期:2010-06-20 发布日期:2010-06-20
  • 作者简介:张彦,女,1969年生,硕士,副研究员,从事同位素地质年代学研究,E-mail: zhangyan@cags.net.cn
  • 基金资助:

    国家自然科学基金(编号40773043);基本科研业务费项目(J0711;JB0708);国土资源大调查项目(1212010761401)

Detailed Procedure for Determining the Helium Diffusion
Parameter in Apatite

ZHANG  Yan, CHEN  Wen, LIU Xin-yu   

  1. Ar-Ar Geochronology Laboratory, Institute of Geology, Chinese Academy of Geological Sciences, Beijing 100037, China
  • Received:2009-06-23 Revised:2010-03-15 Online:2010-06-20 Published:2010-06-20

摘要:

基于Arrhenius 关系图解, 详细介绍了磷灰石中氦扩散参数及其封闭温度的确定方法及求解过程,并以Wolf et al(1996)
的阶段升温数据为例进行了演示性和验证性计算。主要求解步骤如下:(1)对磷灰石做阶段升温试验,得到每一阶段的氦
累积丢失分数;(2)根据球形颗粒的扩散模型公式计算出每一阶段的lnD/a2;(3)以lnD/a2作为纵坐标,1/ T作为横坐标作
图(T 为阶段升温试验每一阶段的温度),如果在某一温度范围内得到的是一条直线,说明在这一温度范围内磷灰石中氦的
扩散遵从Arrhenius关系式,那么在Arrhenius关系图解上,扩散参数E 和ln D0/a2可通过直线的斜率和截距得到;(4)把E 和
lnD0/ a2代入封闭温度的计算公式,通过迭代计算就可得到磷灰石的封闭温度。求解过程中需要注意两个关键问题:(1)由
于1/ T 为小数,对于磷灰石的氦扩散数据来说,需乘以104 将其化为整数,与纵坐标的数值为同一数量级;(2)计算封闭温
度时,需把冷却速率的单位从“℃/Ma”转化为“ K/s”,转化时需注意1 ℃ /Ma = 1 K/Ma。

关键词: 磷灰石, 氦扩散参数, 封闭温度

Abstract:

In this paper we present the detailed procedure to calculate the helium diffusion parameter and its closure temperature in apatite. We use the step-heating data of Wolf (1996)to show this procedure and its correction. This method is based on the Arrhenius relationship which describes the diffusivity variation with temperature. This procedure is as follows: (1) Do step-heating experiment to get the cumulative gas release fraction. (2) Calculate the lnD/a2 of every step using spherical diffusion equation. (3) Make Arrhenius plot, taking the lnD/a2 as y-axis and the 1/T (T is the temperature of every step) as x-axis. In this plot, if the points defined by lnD/a2 values and the 1/T values within some temperature range show linearity, it indicates that within this temperature range the apatite diffusion obeys the Arrhenius relationship and thus we can get the diffusion parameter E from the slope and the LnD0/a 2 from the intercept in the Arrhenius plot. (4) Insert the value of E and LnD0/a 2 into the closure temperature equation to calculate the closure temperature using iterative method. In addition, when we use the above method to calculate the helium diffusion parameter, we should pay attention to the following aspects: (1) Because 1/T is decimal fraction, for apatite diffusion data, 1/T should be multiplied by 104, making it the same order with ln D/ a2. (2) When calculating the closure temperature, the unit of the cooling rate ℃/Ma should be converted to K/s. For the cooling rate, 1 ℃/Ma equals to 1 K/Ma.

Key words: apatite, diffusion parameter, closure temperature

中图分类号: