← 返回 资料库 › 核化学与放射化学 › 应用 核化学与放射化学
放射性同位素测年 放射性衰变常数与母体-子体同位素比如何构建地球年代学时钟,涵盖从考古碳-14 测年到等时线及地球年龄测定的广阔领域。
直觉 直觉:封存在岩石中的放射性原子沙漏 沙漏之所以能够精确计时,是因为细沙以恒定速率由上部容器落入下部容器。在自然界中,放射性原子核构成了天然的微观原子沙漏:不稳定的母核依据严格的核量子跃迁概率衰变为放射成因的子核,这一过程完全不受外界温度、地层压力或化学键环境的影响。当岩浆冷却结晶为固体矿物或生物生命活动停止呼吸代谢的瞬间,沙漏被正式倒置启动:母体原子被牢固锁死在晶格节点中,而子体产物的持续累积便忠实记录下地质历史的浩瀚时间跨度。
选择地质年代计以观察其衰变轨迹:5 万年尺度内 14C 剩余份额的指数衰减(标有 5.73、11.5 和 50 千年标记)、直至 45.4 亿年地球年龄的放射成因 40Ar/40K 累积比值,以及追溯数十亿年行星演化历史的 207Pb /206Pb* 比值曲线。 中学阶段 中学阶段:基本年代方程与碳-14 测年法 定义: 基本放射性年代方程与封闭体系假设
对于衰变常数为 λ \lambda λ 的放射性母核 P P P ,其衰变积累的放射成因子核数目为 D ∗ = P ( e λ t − 1 ) D^* = P(e^{\lambda t} - 1) D ∗ = P ( e λ t − 1 ) 。解出结晶年龄 t t t 得到 t = 1 λ ln ( 1 + D ∗ P ) t = \frac{1}{\lambda}\ln\left(1 + \frac{D^*}{P}\right) t = λ 1 ln ( 1 + P D ∗ ) 。该方程所求得的年龄要代表真实的地质年代,前提是该矿物自冷却至其特有的封闭温度 (T c T_c T c )以下以来,晶格系统必须保持严格的封闭体系 (既无母体或子体元素的流失,也无外界原子的侵入)。
放射性碳测年依托于高层大气中的宇宙射线核反应:次级热中子轰击氮原子核产生碳-14:X 7 14 X 2 7 2 14 N + X 0 1 X 2 0 2 1 n → X 6 14 X 2 6 2 14 C + X 1 1 X 2 1 2 1 p \ce{^{14}_7N + ^1_0n -> ^{14}_6C + ^1_1p} X 7 14 X 2 7 2 14 N + X 0 1 X 2 0 2 1 n X 6 14 X 2 6 2 14 C + X 1 1 X 2 1 2 1 p 。新生的碳-14 迅速被氧化为 X 14 X 2 2 14 C O X 2 \ce{^{14}CO2} X 14 X 2 2 14 CO X 2 并通过植物光合作用进入全球食物网,使活体生物内的 X 14 X 2 2 14 C / X 12 X 2 2 12 C \ce{^{14}C/^{12}C} X 14 X 2 2 14 C / X 12 X 2 2 12 C 原子比恒定维持在约 1.2 × 10 − 12 1.2 \times 10^{-12} 1.2 × 1 0 − 12 的稳态水平。生物死亡后新陈代谢立即终止,体内的 X 14 X 2 2 14 C \ce{^{14}C} X 14 X 2 2 14 C 便以 t 1 / 2 = 5730 ± 40 年 t_{1/2} = 5730\pm 40\text{ 年} t 1/2 = 5730 ± 40 年 (剑桥半衰期)通过 β − \beta^- β − 衰变自发恢复为氮-14。按照国际地质年代学公约,未经校正的常规年龄采用利比半衰期(5568 年 5568\text{ 年} 5568 年 )并以“距今放射性碳年龄 BP”(以 1950 年为基准点)表示。鉴于古地磁场波动和太阳耀斑活动,必须应用基于树木年轮与湖泊纹泥建立的 IntCal20 树轮校准曲线转换为真实日历年。现代加速器质谱技术(AMS)能直接在毫克级微量样品上计数碳-14 粒子,使测年极限拓展至约 5 万年。
主要地质年代测定同位素体系与目标矿物 衰变体系 衰变模式 母核半衰期 有效测年范围 典型分析矿物与材料 X 14 X 2 2 14 C → X 14 X 2 2 14 N \ce{^{14}C -> ^{14}N} X 14 X 2 2 14 C X 14 X 2 2 14 N Beta- 负衰变 5730 年 300 至 50,000 年 木炭、木材、古骨胶原、贝壳、泥炭 X 40 X 2 2 40 K → X 40 X 2 2 40 A r \ce{^{40}K -> ^{40}Ar} X 40 X 2 2 40 K X 40 X 2 2 40 Ar 轨道电子俘获(10.7%) 12.48 亿年 10 万年至 45 亿年 透长石、黑云母、白云母、角闪石 X 87 X 2 2 87 R b → X 87 X 2 2 87 S r \ce{^{87}Rb->^{87}Sr} X 87 X 2 2 87 Rb X 87 X 2 2 87 Sr Beta- 负衰变 488 亿年 1000 万年至 46 亿年 云母、钾长石、花岗岩类全岩等时线 X 238 X 2 2 238 U → X 206 X 2 2 206 P b \ce{^{238}U->^{206}Pb} X 238 X 2 2 238 U X 206 X 2 2 206 Pb Alpha 与 Beta 级联衰变链 44.68 亿年 100 万年至 46 亿年 锆石、斜锆石、独居石、榍石 X 235 X 2 2 235 U → X 207 X 2 2 207 P b \ce{^{235}U->^{207}Pb} X 235 X 2 2 235 U X 207 X 2 2 207 Pb Alpha 与 Beta 级联衰变链 7.038 亿年 100 万年至 46 亿年 与 238U 配对构建谐和曲线与铅铅等时线 X 147 X 2 2 147 S m → X 143 X 2 2 143 N d \ce{^{147}Sm->^{143}Nd} X 147 X 2 2 147 Sm X 143 X 2 2 143 Nd Alpha 衰变 1060 亿年 1 亿年至 46 亿年 基性玄武岩、月球斜长岩、球粒陨石
例题: 考古木炭的放射性常规碳-14 年龄计算
某考古发掘遗址出土的一根古雪松木横梁,其实测比放射性活度恰好为工业革命前现代木材活度标准的 12.5 % 12.5\% 12.5% (即 1 / 8 1/8 1/8 )。试使用标准的剑桥物理半衰期 t 1 / 2 = 5730 年 t_{1/2} = 5730\text{ 年} t 1/2 = 5730 年 ,计算该木制文物在进行树木年轮曲线校准前的常规同位素测年年龄。
解答 放射性剩余活度分数刚好为 N / N 0 = 0.125 = ( 1 / 2 ) 3 N/N_0 = 0.125 = (1/2)^3 N / N 0 = 0.125 = ( 1/2 ) 3 ,表明母核经历了整整 3 个半衰期的衰变过程:t = 3 × t 1 / 2 = 3 × 5730 = 17190 年 t = 3 \times t_{1/2} = 3 \times 5730 = 17190\text{ 年} t = 3 × t 1/2 = 3 × 5730 = 17190 年 。代入指数衰变方程严格计算:t = − t 1 / 2 ln 2 ln ( N N 0 ) = − 5730 0.69315 × ( − 2.07944 ) ≈ 17190 年 t = -\frac{t_{1/2}}{\ln 2}\ln\left(\frac{N}{N_0}\right) = -\frac{5730}{0.69315}\times (-2.07944) \approx 17190\text{ 年} t = − l n 2 t 1/2 ln ( N 0 N ) = − 0.69315 5730 × ( − 2.07944 ) ≈ 17190 年 。在考古学实践中,这一未经校正的常规放射性碳年龄需进一步通过 IntCal20 曲线校准,以消除历史大气浓度波动,最终换算为日历公历年份(cal BP)。
大学 大学阶段:K-Ar 分支衰变、40Ar/39Ar 阶梯加热与 Rb-Sr 等时线 t = 1 λ ln ( 1 + λ λ e 40 Ar ∗ 40 K ) with λ = 5.543 × 10 − 10 y − 1 , λ e = 0.580 × 10 − 10 y − 1 t = \frac{1}{\lambda}\ln\left(1 + \frac{\lambda}{\lambda_e}\frac{^{40}\text{Ar}^*}{^{40}\text{K}}\right) \quad \text{with} \quad \lambda = 5.543 \times 10^{-10}\ \text{y}^{-1},\ \lambda_e = 0.580 \times 10^{-10}\ \text{y}^{-1} t = λ 1 ln ( 1 + λ e λ 40 K 40 Ar ∗ ) with λ = 5.543 × 1 0 − 10 y − 1 , λ e = 0.580 × 1 0 − 10 y − 1 钾-40 发生双通道分支衰变:轨道电子俘获衰变为惰性气体 X 40 X 2 2 40 A r \ce{^{40}Ar} X 40 X 2 2 40 Ar (分支比 10.72 % 10.72\% 10.72% ,λ e = 0.580 × 10 − 10 年 − 1 \lambda_e = 0.580 \times 10^{-10}\ \text{年}^{-1} λ e = 0.580 × 1 0 − 10 年 − 1 ),以及 β − \beta^- β − 负衰变为 X 40 X 2 2 40 C a \ce{^{40}Ca} X 40 X 2 2 40 Ca (分支比 89.28 % 89.28\% 89.28% ,λ β = 4.963 × 10 − 10 年 − 1 \lambda_\beta = 4.963 \times 10^{-10}\ \text{年}^{-1} λ β = 4.963 × 1 0 − 10 年 − 1 ),总衰变常数 λ = 5.543 × 10 − 10 年 − 1 \lambda = 5.543 \times 10^{-10}\ \text{年}^{-1} λ = 5.543 × 1 0 − 10 年 − 1 (半衰期 12.48 亿年)。在地下炽热岩浆中,挥发性极强的氩气完全逸出脱气,使矿物初始氩含量清零(X 40 X 2 2 40 A r 0 ≈ 0 \ce{^{40}Ar}_0 \approx 0 X 40 X 2 2 40 Ar 0 ≈ 0 )。一旦岩浆喷发冷却为透长石或黑云母晶体,放射成因的 X 40 X 2 2 40 A r ∗ \ce{^{40}Ar}^* X 40 X 2 2 40 Ar ∗ 便被封闭在矿物晶格中。现代 X 40 X 2 2 40 A r / X 39 X 2 2 39 A r \ce{^{40}Ar/^{39}Ar} X 40 X 2 2 40 Ar / X 39 X 2 2 39 Ar 测年技术先将样品送入反应堆快中子轰击,通过 X 39 X 2 2 39 K ( n , p ) X 39 X 2 2 39 A r \ce{^{39}K}(n,p)\ce{^{39}Ar} X 39 X 2 2 39 K ( n , p ) X 39 X 2 2 39 Ar 反应建立同位素关联,再在惰性气体质谱计中进行激光阶梯升温释放分析,由平坦的年龄坪谱能够准确剔除后期变质加热导致的放射氩流失或过剩氩干扰。
87 Sr 86 Sr = ( 87 Sr 86 Sr ) 0 + 87 Rb 86 Sr ( e λ t − 1 ) \frac{^{87}\text{Sr}}{^{86}\text{Sr}} = \left(\frac{^{87}\text{Sr}}{^{86}\text{Sr}}\right)_0 + \frac{^{87}\text{Rb}}{^{86}\text{Sr}}\left(e^{\lambda t} - 1\right) 86 Sr 87 Sr = ( 86 Sr 87 Sr ) 0 + 86 Sr 87 Rb ( e λ t − 1 ) 当矿物结晶时含有非零的初始子体同位素 D 0 D_0 D 0 时,单一样品的放射性测量在数学上无法求解年龄。等时线法 (Isochron method)通过引入子体元素的非放射成因稳定同位素进行比值归一化,圆满解决了这一难题。例如在铷-锶测年体系(X 87 X 2 2 87 R b → X 87 X 2 2 87 S r \ce{^{87}Rb -> ^{87}Sr} X 87 X 2 2 87 Rb X 87 X 2 2 87 Sr ,t 1 / 2 = 48.8 亿年 t_{1/2} = 48.8\text{ 亿年} t 1/2 = 48.8 亿年 ,λ = 1.42 × 10 − 11 年 − 1 \lambda = 1.42 \times 10^{-11}\ \text{年}^{-1} λ = 1.42 × 1 0 − 11 年 − 1 )中以稳定的 X 86 X 2 2 86 S r \ce{^{86}Sr} X 86 X 2 2 86 Sr 作为分母。同一岩浆房同源结晶的多种矿物,由于结晶化学分异导致元素比 R b / S r \ce{Rb/Sr} Rb / Sr 各不相同,但在固化瞬间共享完全均一的初始锶同位素比 ( X 87 X 2 2 87 S r / X 86 X 2 2 86 S r ) 0 (\ce{^{87}Sr/^{86}Sr})_0 ( X 87 X 2 2 87 Sr / X 86 X 2 2 86 Sr ) 0 。以 X 87 X 2 2 87 R b / X 86 X 2 2 86 S r \ce{^{87}Rb/^{86}Sr} X 87 X 2 2 87 Rb / X 86 X 2 2 86 Sr 为横坐标、X 87 X 2 2 87 S r / X 86 X 2 2 86 S r \ce{^{87}Sr/^{86}Sr} X 87 X 2 2 87 Sr / X 86 X 2 2 86 Sr 为纵坐标作图,数据点连成一条笔直的等时线,其斜率 m = e λ t − 1 ≈ λ t m = e^{\lambda t} - 1 \approx \lambda t m = e λ t − 1 ≈ λ t 直接决定了矿物的结晶年龄 t t t ,完全无需事先假设初始子体丰度。
例题: 火山凝灰岩透长石晶体的钾-氩同位素年代测定
对火山凝灰岩夹层中分选出的透长石单晶进行高精度同位素质谱分析,测得放射成因氩与母体钾的原子数比为 X 40 X 2 2 40 A r ∗ X 40 X 2 2 40 K = 0.0950 \frac{\ce{^{40}Ar}^*}{\ce{^{40}K}} = 0.0950 X 40 X 2 2 40 K X 40 X 2 2 40 Ar ∗ = 0.0950 。已知总衰变常数 λ = 5.543 × 10 − 10 年 − 1 \lambda = 5.543 \times 10^{-10}\ \text{年}^{-1} λ = 5.543 × 1 0 − 10 年 − 1 ,电子俘获分支常数 λ e = 0.580 × 10 − 10 年 − 1 \lambda_e = 0.580 \times 10^{-10}\ \text{年}^{-1} λ e = 0.580 × 1 0 − 10 年 − 1 ,试计算该火山岩层的结晶地质年龄(以百万年 Ma 为单位)。
解答 代入 K-Ar 分支衰变测年公式:t = 1 λ ln ( 1 + λ λ e X 40 X 2 2 40 A r ∗ X 40 X 2 2 40 K ) t = \frac{1}{\lambda}\ln\left(1 + \frac{\lambda}{\lambda_e}\frac{\ce{^{40}Ar}^*}{\ce{^{40}K}}\right) t = λ 1 ln ( 1 + λ e λ X 40 X 2 2 40 K X 40 X 2 2 40 Ar ∗ ) 。首先计算分支比修正系数:λ λ e = 5.543 × 10 − 10 0.580 × 10 − 10 ≈ 9.5569 \frac{\lambda}{\lambda_e} = \frac{5.543 \times 10^{-10}}{0.580 \times 10^{-10}} \approx 9.5569 λ e λ = 0.580 × 1 0 − 10 5.543 × 1 0 − 10 ≈ 9.5569 。代入对数项真数:1 + 9.5569 × 0.0950 = 1 + 0.9079 = 1.9079 1 + 9.5569 \times 0.0950 = 1 + 0.9079 = 1.9079 1 + 9.5569 × 0.0950 = 1 + 0.9079 = 1.9079 。计算自然对数值:ln ( 1.9079 ) ≈ 0.6460 \ln(1.9079) \approx 0.6460 ln ( 1.9079 ) ≈ 0.6460 。除以总衰变常数 λ \lambda λ :t = 0.6460 5.543 × 10 − 10 年 − 1 ≈ 1.1654 × 10 9 年 = 1165.4 Ma t = \frac{0.6460}{5.543 \times 10^{-10}\ \text{年}^{-1}} \approx 1.1654 \times 10^9\text{ 年} = 1165.4\text{ Ma} t = 5.543 × 1 0 − 10 年 − 1 0.6460 ≈ 1.1654 × 1 0 9 年 = 1165.4 Ma (即 11.654 亿年)。
常见错误. 在经历后期变质作用或风化淋滤的岩石中盲目假设封闭体系会导致严重的地质年代判读错误。若岩石被再次加热超过矿物的封闭温度,晶格内的子体原子(如气态的 40Ar 或不相容离子 Pb2+)将通过固相扩散逃逸到晶界之外,导致测得的表观年龄仅代表后期的变质冷却事件而非岩石最初的岩浆结晶年代;相反,岩浆在深部高压流体环境下结晶时极易在流体包裹体中捕获未排尽的“过剩氩”,导致测算出的年龄比实际地质时代虚假偏老数亿年。在放射性碳测年领域,未经树轮校准的“距今放射性碳年龄 BP”绝不能与真实日历年画等号。 进阶 高级阶段:U-Pb 谐和线、Jack Hills 锆石与地球的年龄 锆石(Z r S i O X 4 \ce{ZrSiO4} ZrSiO X 4 )被公认为同位素地质年代学的金标准。在岩浆结晶过程中,锆石四方岛状硅酸盐晶格中的四价锆位点(r = 84 pm r = 84\text{ pm} r = 84 pm )容许电荷与半径相近的 U X 4 + \ce{U^4+} U X 4 + (r = 100 pm r = 100\text{ pm} r = 100 pm )以类质同象方式大量赋存,而以超过 10 5 10^5 1 0 5 的排斥倍数彻底排斥半径较大且为二价的 P b X 2 + \ce{Pb^2+} Pb X 2 + (r = 119 pm r = 119\text{ pm} r = 119 pm ),从而使新形成锆石中的初始普通铅几乎为零(P b 0 ≈ 0 \ce{Pb}_0 \approx 0 Pb 0 ≈ 0 )。锆石内部同时运行着半衰期迥异的两个独立铀衰变时钟:X 238 X 2 2 238 U → X 206 X 2 2 206 P b \ce{^{238}U -> ^{206}Pb} X 238 X 2 2 238 U X 206 X 2 2 206 Pb (t 1 / 2 = 44.68 亿年 t_{1/2} = 44.68\text{ 亿年} t 1/2 = 44.68 亿年 ,λ 238 = 1.55125 × 10 − 10 年 − 1 \lambda_{238} = 1.55125 \times 10^{-10}\ \text{年}^{-1} λ 238 = 1.55125 × 1 0 − 10 年 − 1 )和 X 235 X 2 2 235 U → X 207 X 2 2 207 P b \ce{^{235}U -> ^{207}Pb} X 235 X 2 2 235 U X 207 X 2 2 207 Pb (t 1 / 2 = 7.038 亿年 t_{1/2} = 7.038\text{ 亿年} t 1/2 = 7.038 亿年 ,λ 235 = 9.8485 × 10 − 10 年 − 1 \lambda_{235} = 9.8485 \times 10^{-10}\ \text{年}^{-1} λ 235 = 9.8485 × 1 0 − 10 年 − 1 )。以 X 207 X 2 2 207 P b ∗ / X 235 X 2 2 235 U \ce{^{207}Pb}^*/\ce{^{235}U} X 207 X 2 2 207 Pb ∗ / X 235 X 2 2 235 U 为横轴、X 206 X 2 2 206 P b ∗ / X 238 X 2 2 238 U \ce{^{206}Pb}^*/\ce{^{238}U} X 206 X 2 2 206 Pb ∗ / X 238 X 2 2 238 U 为纵轴绘制的参数轨迹称为韦瑟里尔(Wetherill)谐和曲线 。保持封闭体系的样品数据点严密落在谐和曲线上;若岩石在后期热事件中发生放射铅的扩散流失,测点将构成一条笔直的不谐和线 (Discordia),其与谐和线的上交点代表岩浆初生结晶年龄,下交点则代表变质扰动年代。采自西澳大利亚杰克山(Jack Hills)的碎屑锆石测得了高达 44.04 亿年 44.04\text{ 亿年} 44.04 亿年 的谐和年龄,成为地球冥古宙原始地壳存在液态水与长英质陆壳的直接见证。
对两个铀-铅衰变方程两端求比值,便巧妙消除了对母体元素铀绝对化学浓度的测量依赖:X 207 X 2 2 207 P b ∗ X 206 X 2 2 206 P b ∗ = X 235 X 2 2 235 U X 238 X 2 2 238 U ( e λ 235 t − 1 e λ 238 t − 1 ) = 1 137.88 ( e λ 235 t − 1 e λ 238 t − 1 ) \frac{\ce{^{207}Pb}^*}{\ce{^{206}Pb}^*} = \frac{\ce{^{235}U}}{\ce{^{238}U}}\left(\frac{e^{\lambda_{235} t} - 1}{e^{\lambda_{238} t} - 1}\right) = \frac{1}{137.88}\left(\frac{e^{\lambda_{235} t} - 1}{e^{\lambda_{238} t} - 1}\right) X 206 X 2 2 206 Pb ∗ X 207 X 2 2 207 Pb ∗ = X 238 X 2 2 238 U X 235 X 2 2 235 U ( e λ 238 t − 1 e λ 235 t − 1 ) = 137.88 1 ( e λ 238 t − 1 e λ 235 t − 1 ) ,其中 137.88 137.88 137.88 为现代地球与太阳系的天然 X 238 X 2 2 238 U / X 235 X 2 2 235 U \ce{^{238}U/^{235}U} X 238 X 2 2 238 U / X 235 X 2 2 235 U 同位素比值。1956 年,克莱尔·帕特森(Clair Patterson)测定了不含铀的峡谷暗黑(Canyon Diablo)铁陨石陨硫铁相(F e S \ce{FeS} FeS )中的铅同位素组成,确立了太阳系初生普通铅的原始基准。通过将多种石陨石与现代大洋深海沉积物铅同位素联合投在一条全球共生的铅铅等时线上,帕特森最终测定地球的年龄为 45.5 ± 0.7 亿年 45.5 \pm 0.7\text{ 亿年} 45.5 ± 0.7 亿年 (即 45.4 亿年 45.4\text{ 亿年} 45.4 亿年 )。现代对原始球粒陨石中富钙铝包裹体(CAI)的超高精度铀-铅测年,将太阳星云最早固相物质凝聚时间精确锚定在 45.67 亿年 45.67\text{ 亿年} 45.67 亿年 前。
历史注记
1905 年,欧内斯特·卢瑟福(Ernest Rutherford)首次敏锐地提出利用天然放射性衰变建立定量地质时钟的构想,根据矿物中积累的铀-氦比测出数亿年的古老年龄,彻底打破了开尔文勋爵(Lord Kelvin)依据热传导冷却算出的狭隘年轻地球假说。1913 年,弗雷德里克·索迪(Frederick Soddy)提出“同位素”划时代概念,证明化学性质完全相同的元素可具有不同原子质量和独立核衰变路径。此后,伯特拉姆·博尔特伍德与阿瑟·霍姆斯建立了首个铀-铅地质年代表,威拉德·利比因发明碳-14 测年荣获 1960 年诺贝尔化学奖,而克莱尔·帕特森最终测定了地球 45.5 亿年的真实年龄。
欧内斯特·卢瑟福 弗雷德里克·索迪
参考文献 Atmospheric Helium Three and Radiocarbon from Cosmic Radiation · W. F. Libby, 1946 Precision Measurement of Half-Lives and Specific Activities of 235U and 238U · A. H. Jaffey, K. F. Flynn, L. E. Glendenin, W. C. Bentley, A. M. Essling, 1971 The IntCal20 Northern Hemisphere Radiocarbon Age Calibration Curve (0-55 cal kBP) · P. J. Reimer et al., 2020