跳到正文
返回基本对称性与中微子
至真探索 · 基本对称性与中微子

地球中微子与反应堆中微子分离分析

2026-07-02

你面前有两种反中微子。它们来自完全不同的源头——一种来自你脚下的地球深处，另一种来自五十公里外核电站的堆芯。它们通过完全相同的物理过程被探测到：反β衰变，ν̄ₑ + p → e⁺ + n。你给探测器装上两万吨液体闪烁体，看着数以千计的微弱闪光。

问题是：你不知道每一个闪光来自哪一种中微子。

分离它们——这个问题的难度被极度低估了。

为什么分离如此之难

地球中微子全部来自铀-238和钍-232的衰变链。它们的能谱在约1.8到3.3 MeV的区间内分布——恰好是核反应堆反中微子谱的低能尾巴所在的位置。在液态闪烁体探测器中，反β衰变产生一个正电子和一个中子。正电子的能量等于入射反中微子能量减去1.8 MeV（反应阈值），所以"探测到的信号"是一段连续谱，而在这个谱的低能端——大约1到4 MeV的"可见能量"窗口——地球中微子和反应堆中微子完全混在一起。

如果你手上有足够多的数据，你可以做一个能谱拟合——利用两个组分之间细微的谱形差异来分开它们。地球中微子来自铀-238和钍-232的β衰变，它们的能量分布由成熟的核衰变数据库定义——谱形的不确定性大约在5%。反应堆反中微子来自四种裂变核素——铀-235、钚-239、铀-238和钚-241——它们的谱形有一个顽固的"5 MeV隆起"问题，系统不确定性在2-5%之间。

但真正让事情变得困难的是动态信号。反应堆的功率随电网负载波动——每天、每周、每个季节都在变化。你今天测到的反应堆中微子通量取决于阳江和台山核电站昨天的运行状态。而地球中微子——好消息是——在过去几十亿年里几乎没有变化。坏消息是：这意味着你没有办法通过时间调制来区分它们。

你只有一个静态的谱形拟合。两个组分的形状有微弱的差异。你有统计误差——事件数是有限的。你有系统误差——反应堆谱形模型的能力有限。你有地质学上的"先验"——对地球地壳中铀和钍含量的不完全了解，这个先验自身带有约20%的不确定性。你把这些全扔进一个分bin的似然函数里，同时拟合振荡参数、反应堆归一化系数、地球中微子铀和钍的归一化系数、以及几十个本底参数。

这就是反β衰变能谱拟合——人类做过的最困难的分离问题之一。

三代探测器：从发现到精密

KamLAND——日本神冈的千吨级液闪探测器——在2005年发表了第一篇地球中微子观测论文。这是地学的革命性时刻：人类第一次用中微子而不是地震波，直接观察了地球内部的放射性。在那之前，地质学家只能从地表岩石样本和热流测量来推断地球内部有多少铀和钍。KamLAND用中微子直接测量了它。

到2020年代，KamLAND和Borexino（意大利Gran Sasso，300吨液闪）各自积累超过十年的数据。两者独立测量出的铀+钍产生的放射性热功率分别是20⁺⁶₋₇ TW（KamLAND）和21⁺⁷₋₈ TW（Borexino）。加上钾-40（其反中微子能量低于反β衰变阈值，无法被直接探测，但从地球化学约束推断贡献约4 TW），地球的总放射性热功率约为24 TW。地球表面流出的总热量是46 ± 3 TW。差值——约23 TW——来自原始吸积热（地球形成时的残余热）和内核结晶潜热。

这个数字的分量很重。如果放射性占比太高（>30 TW），地球的冷却就是缓慢的放射性余烬，板块构造在地球历史早期就已经停止。如果占比太低（<15 TW），地球就是一个仍然在快速冷却的年轻引擎——板块构造可以维持更长的时间。24 ± 8 TW正处于中间：既有足够的放射性来维持对流行星发动机，也有足够的原始热来确保地球还在冷却。

但KamLAND和Borexino能告诉我们的只有这么多。它们的地球中微子通量的总精度约15-18%。而且它们的信号中大约一半来自近场岩石圈——探测器周围几百公里范围内的地壳——这个局部地质的铀钍含量偏差是地表岩石样本无法精确约束的。

2025年，SNO+——加拿大安大略省萨德伯里SNOLAB地下2公里的液闪探测器——发表了西半球首次地球中微子测量。SNO+建在加拿大地盾上——这是一块厚度约40公里的前寒武纪大陆地壳，比KamLAND和Borexino的地壳更厚、更古老。这意味着更大比例的信号来自大陆岩石圈而非地幔——对于约束北美大陆地壳的铀钍预算，这是一个独特的地质窗口。

SNO+团队的完整分析使用2022年5月到2025年7月的685天数据。他们的反β衰变谱拟合将地球中微子通量同时与Δm²₂₁、sin²θ₁₂、反应堆归一化系数和本底一起作为自由参数——这是第一次由单个探测器完成的同时多参数地球中微子提取。结果：总通量49⁺¹³₋₁₂ TNU（TNU = 地球中微子单位 = 每年每10³²个自由质子产生1个事件），铀贡献约38个事件，钍约11个事件，U/Th比值约3.3——与球粒陨石结构的预测值约3.9一致。该结果在女王大学主办的Neutrino Geosciences 2025大会上首次公布，随后在Physical Review Letters（135卷，121801，2025年9月）上正式发表，精细版发表在2026年4月的arXiv:2604.05746上。

JUNO的分离策略：六个层面的同时约束

JUNO（江门地下中微子观测站）在全球地球中微子测量版图上是一个质的跳跃。它每年探测约400个地球中微子——超过KamLAND和Borexino十多年数据的总和。它的物理程序被设计为同样的反β衰变谱中同时做四件事：测量中微子质量顺序、精确化sin²θ₁₂和Δm²₂₁、提取地球中微子通量、以及测量地球内部的铀/钍比值。

JUNO能同时做这四件事，得益于一种仪器：TAO——台山反中微子观测站。

TAO是JUNO的卫星探测器——2.8吨钆掺杂液闪，距离台山核电站一个4.6 GW热功率的反应堆芯仅44米。在-50°C的低温下运行，硅光电倍增管的暗噪声被压制了约三个数量级。其能量分辨率优于2%（1 MeV处），光学覆盖率约95%，每年将记录约70万个反β衰变事件。这个距离使它几乎看不到中微子振荡——所以它测量的是反应堆的"原始"反中微子能谱，不受振荡的调制。JUNO——52.5公里远的母探测器——看到的是这个谱经过质量顺序敏感的振荡调制后的版本。蒋涛、张金龙等在2026年JHEP（arXiv:2512.00172）上的分析展示了一个精细的分bin似然拟合框架：TAO压制反应堆谱形不确定性到约1%的水平，从而释放出JUNO的低能谱区的自由度——那些自由度在2-4 MeV的窗口里，正好是地球中微子活着的区域。

当TAO把反应堆谱形的不确定度从5%压到1%时，JUNO的能谱拟合就不再被"这个隆起到底是反应堆谱本身的形状还是地球中微子贡献"所困扰。地球中微子信号从光谱噪声中浮现出来。

这是JUNO分离策略的核心：用一个近场探测器约束入射谱，使远场探测器中的残余谱形变化能被归因于振荡物理和地球中微子——二者在能量空间中占据不同的位置，一旦反应堆谱形的不确定性被压制到1%，它们就不再"侵蚀"对方的本征空间。

分离还需要五个额外的约束层。第二层是JULOC-I——JUNO的第一个三维综合地壳模型。韩如冰等人在2025-2026年开发了这个模型，结合了地震数据、重力数据、岩石样本地球化学测量和热流数据。关键发现令人吃惊：华南地区的铀和钍浓度高于全球平均水平。这意味着JUNO下方的岩石圈产生的信号比全球均值模型预测的更强。如果你不修正这个局部增强，你会错误地把它归因于地幔信号——从而高估地球内部的放射性。JULOC-I是第一道防线：用局部地质信息减去岩石圈的已知贡献。

第三层是背景约束。JUNO的反β衰变事件选择——在2024-2025年达到成熟阶段——使用了可解释机器学习方法。Agostini、Clementi等人在Physics Letters B（2025年，来自arXiv:2406.12901）上发表了一个全连接神经网络判别器的结果：与传统的截断式选择相比，探测器边缘区域——在这里本底率最高，传统方法丢弃了大量信号——的信号效率显著提高。探测器边缘恰好对地球中微子特别重要：因为地球中微子没有方向性，它们均匀地填充整个探测器体积。提高探测器全体积的均匀选择效率，等于在不用建更大型探测器的情况下增加了有效曝光。

第四层是Th/U分离——这可能是整个领域最棘手的问题。铀和钍的地球中微子谱在能量空间中强烈反相关。如果你试图把它们作为两个独立的自由参数拟合，它们会互相推拉——一个被拉高时另一个被压低。JUNO的灵敏度（Li et al., Chinese Physics C, 2026年2月）表明：固定Th/U比值为球粒陨石值3.9时，十年后地球中微子信号的总精度约8%。放开Th/U比值自由拟合时，铀的精度下降到约30%，钍下降到约35%，总和约15%。这意味着JUNO在没有外部先验的情况下，无法轻易地单独确定铀和钍——但它可以使二者的总约束达到世界最精。

第五层是地幔信号的提取。在大约六年的数据积累后，如果岩石圈通量被约束在约15%以内，JUNO有望以3倍标准偏差看到地幔成分。这是KamLAND和Borexino从未实现的——它们的地幔信号检测仍然是不确定的，取决于岩石圈通量模型。JUNO的统计量和TAO的反应堆谱约束，使地幔信号提取首次从"可能"变为"可规划"。

最新结果与现存张力

在2026年6月22日Neutrino 2026大会（加州大学尔湾分校）上，Cristobal Ignacio Morales Reveco代表JUNO合作组展示了首个地球中微子测量结果的初步比较——与来自Li et al.（CPC 2026）和JULOC-I地壳模型的地质模型预测值进行比对。虽然正式的地球中微子通量数值还在分析中，取数效率超过90%的207天数据已让合作组对"一年内产出世界最精地球中微子通量结果"表达高度信心。

今年早些时候一篇值得注意的arXiv论文（Han、Qian和Chen, 2026）题为"Worldwide Reactor Neutrino Propagation to Underground Labs: Matter Effects and Flux Predictions"。作者指出，随着地球中微子测量迈向亚百分比精度，另一个此前被忽略的误差源必须被正视：地球的物质效应——即中微子在穿过数千公里地球内部的传播过程中经历的MSW-型共振转换——会改变从远处反应堆到达JUNO的能谱形状。当反应堆谱的精度被TAO推到约1%时，物质效应造成的谱形变形——几个千分之一的量级——就不能再安全地忽略了。精密不仅是技术的代价——它是一种物理的精炼：你必须把每一个此前被忽略的小效应找出来、量化它、然后把它放进全局似然函数里，否则它就会伪装成地球中微子信号。

在JUNO走到这一步之前，KamLAND和Borexino已经在同一个问题上形成了一组有趣的对峙。2010年代早期——当时两者都用了当时各自的岩石圈模型——它们反推的地幔铀钍含量存在一些不一致。这件事在2022年被Sammon和McDonough（EPSL, 593, 117684）定量化了：Borexino的近场岩石圈模型低估了Gran Sasso周围钾长石富集火山岩中的铀含量。把修正的岩石圈模型用于Borexino后，KamLAND和Borexino的地球中微子通量收敛到一致的放射性热功率约20 TW（U+Th）。这一修正化解了看似最严重的地球中微子张力——并且确证了分离方法本身在充分理解局部地质的前提下是自洽的。

Wipperfurth、Šrámek和McDonough在2020年发表的JGR综述（125卷，e2019JB018433）量化了这件事的底线：地壳模型的选择可以将总通量预期值整体偏移高达约15%，而地球化学的不确定性相比地球物理模型的不确定性占据了总误差预算的大约70%。这意味着地球中微子的当前瓶颈不在探测器的统计量——而在我们对脚下的岩石圈中到底有多少铀和钍的了解程度。

未来探测器和地幔的最终回答

三个正在发展中的方向值得注视。

一是海底中微子探测器（OBD）。渡边宽人等（东北大学/JAMSTEC）在2024-2025年持续论证在约4公里深的海底放置地球中微子探测器的可行性。夏威夷——远离大陆地壳、下地幔信号流几乎不被岩石圈污染——特别有吸引力。在高地幔放射性的模型中，仅一年观测就可以以3.7σ看到纯地幔信号（McDonough, 2025年1月夏威夷大学研讨会）。

二是JUNO和SNO+将利用它们的地理分散性来联合约束岩石圈与地幔。JUNO主要看华南造山带的厚地壳。SNO+看加拿大地盾的老地壳。KamLAND看日本岛弧的薄地壳。Borexino看意大利造山带的混杂地壳。四个点测出的通量如果与全球岩石圈模型产生系统性差异——比如全部四个点都偏高——问题很可能出在地幔假设。如果只有一个点（比如JUNO下面的华南）偏高——问题大概率是局部地质。

三是钾-40的直接探测。目前没有任何探测器能看到钾-40的反中微子——它的最高能量约1.3 MeV，低于反β衰变的1.8 MeV阈值。但下一代低阈值探测器——基于二次闪烁光或者纯粹溶解钆的高效捕获——可能在2030年代突破这个障碍。如果成功，你将首次直接看到地球内部钾-40的分布——这是地幔对流行为信息中最关键的一块，因为钾的挥发性意味着它主要富集在地壳中，而不是地幔中。确认或否定这个地化预言将深刻影响我们对地球化学分异的理解。

地球中微子是我们拥有的唯一一种可以直接看到地球内部深处化学成分的实验工具。它们的总强度告诉你铀和钍的绝对数量——这直接决定了地球放射性热引擎的马力。它们与反应堆中微子的分离——通过能谱、TAO的高精度谱约束、ML选择器、三维地壳模型、全局似然拟合——是你从地壳和地幔的混合信号中提炼出地幔单独贡献的唯一方法。目前的总放射性热功率收敛在约24 TW——这很不错。但要单独说服地质学家"地幔本身的铀/钍比是球粒陨石值"还是"偏离了球粒陨石值"，需要JUNO级别的统计量——而这恰恰是它在未来几年内即将达到的。

---

参考文献

  1. Li YF, Cao J, Wang YF, Zhan L et al. (JUNO Collaboration). "Prospects for geoneutrino detection with JUNO." *Chinese Physics C* (February 2026). [IHEP, CAS — 10年精度~8%, 3σ地幔检测, TAO约束]
  2. JUNO Collaboration. "Measurement of reactor neutrino oscillation with the first JUNO data." *Nature* 654, 343–348 (2026). DOI: 10.1038/s41586-026-10538-z. [封面文章, 59.1天数据]
  3. 王贻芳 (JUNO Collaboration Spokesperson). "JUNO: First Results on Neutrino Mass Ordering and Geoneutrinos." Neutrino 2026, UC Irvine, June 22, 2026. [207天数据, 正常顺序2.3σ, 取数效率>90%]
  4. SNO+ Collaboration. "Reactor antineutrino oscillations and geoneutrinos in SNO+." arXiv:2604.05746 (2026). [Queen's University / SNOLAB — 685天数据, 49⁺¹³₋₁₂ TNU, U/Th~3.3]
  5. SNO+ Collaboration. "Measurement of reactor antineutrino oscillation at SNO+." *Physical Review Letters* 135, 121801 (2025). arXiv:2505.04469. [首次西半球地球中微子测量]
  6. KamLAND Collaboration. "Partial radiogenic heat model for Earth revealed by geoneutrino measurements." *Nature Geoscience* 4, 647–651 (2011). Updated in multiple analyses through 2022. [Tohoku University / Berkeley]
  7. Borexino Collaboration. "Comprehensive geoneutrino analysis with Borexino." *Physical Review D* 101, 012009 (2020). [INFN / Gran Sasso — U+Th ≈ 21 TW]
  8. Han RB, Qian X, Chen SM et al. "Worldwide reactor neutrino propagation to underground labs: matter effects and flux predictions." arXiv (2026). [物质效应对地球中微子分离的精密挑战]
  9. Han RB et al. (JUNO Collaboration). "JULOC-I: A 3-D integrated refined crustal model for JUNO geoneutrino analysis." (2025-2026). [华南铀/钍高于全球均值 — JUNO近场岩石圈模型]
  10. TAO Collaboration. "The Taishan Antineutrino Observatory." *Journal of High Energy Physics* 02, 109 (2026). [IHEP, CAS — 2.8吨, -50°C SiPM, <2%@1MeV, 2000 IBD/天, 反应堆谱约束~1%]
  11. Agostini F, Clementi A et al. "Interpretable machine learning approach for electron antineutrino selection in a large liquid scintillator detector." *Physics Letters B* 149, 138699 (2025). arXiv:2406.12901. [JUNO ML选择器 — 探测器边缘信号效率提升]
  12. Sammon LG, McDonough WF. "Quantifying Earth's radiogenic heat budget." *Earth and Planetary Science Letters* 593, 117684 (2022). [修正Borexino NFL模型 — 化解KamLAND-Borexino张力]
  13. Wipperfurth SA, Šrámek O, McDonough WF. "Reference models for lithospheric geoneutrino signal." *Journal of Geophysical Research: Solid Earth* 125, e2019JB018433 (2020). [地壳模型不确定性~20%, 化学70%+物理30%]
  14. Huang Y, Chubakov V, Mantovani F, Rudnick RL, McDonough WF. "A reference Earth model for the heat-producing elements and associated geoneutrino flux." *Geochemistry, Geophysics, Geosystems* 14, 2003–2029 (2013). [标准地球模型 — Th/U=3.9, K/U=14000]
  15. Watanabe H et al. (Tohoku University / JAMSTEC). "Ocean-bottom geoneutrino detector: Hawaii as a pure mantle signal site." Seminar, University of Hawaii, January 2025. [3.7σ/1年地幔检测]