揭示北爱尔兰加里东期与莫恩花岗岩潜藏生热潜力：被抑制机载放射性生热信号的地化重校准

《Geothermics》：Unveiling the hidden heat potentiality of the caledonian and mourne granites, Northern Ireland: A geochemical recalibration of suppressed airborne radiogenic heat signals

　　

研究人员揭示了北爱尔兰加里东期（Caledonian）与莫恩（Mourne）花岗岩中，由于泥炭覆盖与冰碛沉积导致机载伽马能谱（AGRS）放射性生热（RHP）被系统性低估的现象。研究人员建立了创新幂律校正模型 cf = 5.25 × RHPAGRS?0.828，以准确反映实际基岩生热。研究人员基于AGRS数据结合更新后的RHP评价对莫恩与加里东期花岗岩进行了重新分类。

研究背景与问题提出

区域地热评估中，放射性生热（Radiogenic Heat Production, RHP）是热梯度、热流及岩石圈热模拟的核心参数，源于长寿命同位素铀（U）、钍（Th）、钾（K）衰变。机载伽马能谱（Airborne Gamma-Ray Spectrometry, AGRS）仅采样地表30–40 cm内γ辐射，在北爱尔兰广泛泥炭与冰碛（till）覆盖区，饱和泥炭作为低密度屏蔽强吸收γ射线，冰碛则通过结构混合与空间平滑改变信号，致使AGRS反演的eU、eTh、K%偏低（如航空K?O较土壤低约20%，eU不足土壤三分之一），从而系统性低估基岩RHP。既往利用AGRS（如Tellus 2005–2006测线）与地球化学样品估算北爱尔兰RHP存在显著分歧，且缺乏针对覆盖层衰减的定量校正统计。为此，研究人员以加里东期Newry火成杂岩（Rathfriland、Newry、Cloghoge岩体）与古近纪Mourne山复合体（G1–G5相）为对象，开展双模型地化重校准研究，论文发表于《Geothermics》。

主要关键技术方法

研究人员整合Tellus AGRS数据（2005–2006年，测线距200 m，海拔56–244 m，RADCOR处理获K%、eTh、eU网格化）与85件Mourne花岗岩历史地球化学样品（XRF、ICP-MS、慢γ计数等，精度±10% 1σ），同时引用加里东期少量样品作验证。采用Rybach（1988）与Carmichael（1989）双RHP公式作敏感性交叉验证，前者含密度ρ，后者取地壳均值2700 kg/m3；通过样品对计算校正因子CF = A_AGRS/A_Geo，经对数线性回归推导幂律CF = 5.25 × RHP_AGRS^?0.828（文中另处系数为5.28，属同模型表述），并将该因子迁移应用于加里东期花岗岩AGRS数据。

研究结果

4.1 加里东期花岗岩RHP：AGRS双模型平均RHP 0.94 μW m^?3（σ±0.27），峰值3.48（Carmichael）与3.57（Rybach, ρ=2680 kg/m3），属极低生热，低于大陆地壳均值1.0–1.5 μW m^?3。eTh、eU相关R2≈0.8，K仅0.40，显示近表K迁移性高。双模型偏差2.5–3.0%，低幅值下常数选择引入可观相对不确定。

4.2 莫恩花岗岩RHP：AGRS平均2.4 μW m^?3（σ±0.80），峰值7.58/7.48；双模型与eTh、eU相关R2=0.93，K伪相关R2=0.70源于泥炭厚度/饱和度协同衰减。Rybach-Carmichael parity R2=1，偏差<2%，高U环境模型选择影响微弱。

4.3 地化验证与模型敏感性（Mourne）：85件样品双模型平均6.64 μW m^?3（σ±1.98），峰值12.91/13.3。U贡献47.1%、Th 47.8%、K 5.9%；U统计控制最强（R2≥0.91），K无相关（R2≥0.026）证其为环境假象。Rybach于ρ=2630 kg/m3更贴近Band +1。

4.4 重校准与CF推导：AGRS/地化中位比2.75（压制约63%）。逐样CF分位：低1.59、中位2.52、高3.14。幂律CF=5.28×A_AGRS^?0.828（R2=0.63，对数回归R2=0.51, p<0.001）。中位AGRS 2.4代入CF得2.54≈2.52。校正后AGRS平均6.3 μW m^?3（σ±2.11），贴近地化6.64。

4.5 加里东期验证：套用Mourne CF分级（≥0.9→1.59；0.5–0.9→2.52；<0.5→3.14）得预测1.94 μW m^?3（σ±0.34），单用中位CF得2.27 μW m^?3（σ±0.68），与Noller等（2015）1.5–2.0及Lee（1986）1.9吻合。

讨论部分总结

研究人员指出AGRS对Mourne原平均2.40经CF 2.52校正为6.3 μW m^?3，跨过5.0经济阈值，由中等升为高生热花岗岩（HHP），比肩Cornubian岩基与Carnsore花岗岩；负指数?0.828量化饱和覆层非线性衰减。加里东期校正后1.9–2.27 μW m^?3佐证区域系统性低估。K通道R2=0.026而U=0.92，表明K波动属泥炭/冰碛环境噪声而非岩性。一维稳态热导模型（T_z=T?+(q_s/k)z?(A/2k)z2）取k=2.7 W/m·K、T?=10℃、Mourne延深6 km，RHP 6.3预测5 km处78.9℃、6 km处85.7℃；加里东期RHP 2.27仅49.2℃@5 km，证实Mourne为北爱深部地热靶区。局限含一维平滑、覆层均厚假设、岩性专适性；工作流可移植至全球泥炭/冰碛区。

结论部分翻译

本研究通过Tellus AGRS与Mourne基岩地化样品经Rybach公式调和，克服饱和泥炭与冰碛对AGRS影响，成功重校准北爱加里东期与Mourne花岗岩RHP。双模型揭示AGRS系统性欠估，导出校正因子将Mourne RHP重定级为6.3 μW m^?3，越地热经济阈值，比肩英爱最高产花岗岩。校正模型于加里东期（2.27 μW m^?3）成功验证，彰显该方法可作覆盖层区隐藏高生热靶区识别之可迁移框架，为北爱向可再生深部地热转型提供新地质依据，使区域潜力由边缘可能转为工业级取热可行目标。