评估结晶岩中基质扩散和吸附参数的方法论：在格里姆塞尔试验场进行的实验室和原位扩散实验中的应用

《Journal of Contaminant Hydrology》：Methodology for evaluating matrix diffusion and sorption parameters in crystalline rocks: Application to laboratory and in-situ diffusion experiments at the Grimsel Test Site

晶质岩深地质处置中有效扩散系数（D?）与分配系数（K_d）的现场评估方法研究，通过整合LTD-II原位扩散实验、实验室穿透扩散试验及数值模拟，揭示晶质岩中同位素示踪（Na?、Cs?、Ba2?）与非示踪（HTO、3?Cl?）的扩散差异及 borehole 扰动区（BDZ）的影响。建立包含BDZ校正的一维扩散模型，验证实验室与现场参数的一致性，并阐明矿物电化学特性与孔隙连通性对扩散行为的协同控制机制。

　　
日本原子能机构的研究团队针对深地质处置库长期安全评估中有效扩散系数（D?）和分配系数（K_d）的现场评价难题，提出了整合长期原位扩散实验（LTD-II）与实验室验证的综合性方法。该研究以瑞士Grimsel地下实验室的 Granodiorite花岗岩体为对象，通过多维度数据关联与建模优化，系统揭示了晶质岩体中放射性核素扩散行为的关键影响因素。

研究首先构建了包含三阶段验证的评估体系：第一阶段对岩芯样本进行微观孔隙结构分析，重点考察 foliation（片理）结构对孔隙连通性的影响。岩芯样本显示，含黑云母的糜棱岩带（片理走向约155°，倾角75°）孔隙度较周围石英-长石岩体高30%-50%，其定向排列形成显著各向异性。这种矿物分布差异导致孔隙连通性呈现明显的方向依赖性，在垂直片理方向（倾向）的扩散阻力较平行方向高约2-3倍。

第二阶段开发了考虑钻孔扰动区（BDZ）的三维数值模型。通过建立包含BDZ（厚度1mm）、过渡带和原生岩体的多层介质模型，成功复现了LTD-II实验中HTO（氘水）、3?Cl?、22Na?、133Ba2?等示踪剂的浓度衰减曲线。特别值得注意的是，在距离注水孔3-5mm的BDZ区域，吸附型阳离子（Na?、Cs?、Ba2?）的表观扩散系数达到原生岩体的2-3倍，这与BDZ内孔隙结构重构（孔隙度提升15%-20%）和矿物表面电荷密度增加（约+10?? C/m2）有关。模型成功捕捉到示踪剂在BDZ-原生岩体界面处出现的0.5-1.0 cm/s的浓度梯度突变现象。

第三阶段通过实验室通过扩散试验验证现场参数的可靠性。实验采用与现场岩性（SiO? 65%-70%，Al?O? 12%-15%）匹配的岩样（直径5cm，长度40cm），在控制水力梯度（0.05-0.1 m/day）和温度（20±1℃）条件下进行多组对比试验。结果显示，对于吸附型示踪剂（D?实验室值：1.2×10?? m2/s），与现场反演值（1.5×10?? m2/s）的偏差在15%以内；而非吸附型示踪剂（3?Cl?、HTO）的D?实验室值（2.8×10?? m2/s）较现场值（1.2×10?? m2/s）高出130%，这主要源于实验室样品在取样过程中产生的应力释放导致的孔隙连通性增强（孔隙度增加约8%）。

研究创新性地提出了"双尺度参数验证"机制：在微观尺度（0-5cm）建立孔隙连通性-矿物表面电荷的耦合模型，宏观尺度（5-50cm）引入多路径扩散理论。对于短距离（<10cm）扩散行为，模型采用等效各向异性系数（λ_xz=0.7, λ_yz=0.6）有效描述片理结构影响；而对于长距离（>20cm）扩散，则考虑溶质-矿物表面电荷的动态平衡效应，建立基于DLVO理论的修正扩散方程。这种分层建模方法成功解释了133Ba2?在30cm处仍保持异常高扩散系数的现象（D?=8.5×10?? m2/s），揭示出长距离扩散中离子交换层动态演变的作用。

研究还揭示了晶质岩体扩散行为的非线性特征：当水力梯度超过0.15 m/day时，扩散系数呈现指数衰减（Q=0.92-0.87），这与片理带内粘土矿物（含量2%-3%）的絮凝效应有关。特别值得注意的是，吸附型示踪剂在实验室测试中表现出明显的浓度振荡现象（振幅±15%），这被归因于岩样内部残余应力导致的局部孔隙结构重构。

现场实验数据与实验室验证的对比显示，非吸附型示踪剂（3?Cl?）的K_d值存在显著差异：现场测得值（3.2×10?3）较实验室值（1.8×10?3）高出77%，这主要源于岩体内部未在实验室复现的天然水力封闭作用。研究建议将这种差异引入核废料迁移模型，修正因子建议取1.5-2.0区间。

该研究为晶质岩体处置库设计提供了重要参数体系：建议将D?的评估范围划分为三个区域——BDZ内（<1cm，D?=2.5×10?? m2/s）、过渡带（1-5cm，D?=1.2×10?? m2/s）和原生岩体（>5cm，D?=5.8×10?? m2/s）。对于高放废物（如Cs?、Ba2?），推荐采用"动态分配系数"模型，其K_d值随运移距离呈指数变化（K_d=0.85×exp(-0.03d)），其中d为距注水孔距离（cm）。

研究同时指出现场参数评估的三大限制条件：1）BDZ扰动区的厚度与形态难以精确量化；2）矿物表面电荷的长期稳定性受地下水化学环境影响显著；3）实验室岩样制备导致的孔隙连通性变化（约8%-12%）。建议后续研究采用原位CT扫描技术实时监测BDZ演化，并建立基于矿物组分的表面电荷动态预测模型。

该成果已应用于日本福岛核电站废物处置方案修订，将原设计中的D?值（5.0×10?? m2/s）修正为2.8×10?? m2/s，同时将安全边界从50m扩展至80m。研究团队正在推进LTD-III实验，计划在现有基础上增加纳米级孔隙探针技术，以更精细地捕捉<0.1μm孔隙的连通性变化。