富含有机质页岩中控制黏土束缚水核磁共振弛豫的界面分子机制及其对核磁共振测井的启示

《Molecules》：Interfacial Molecular Mechanisms Governing the NMR Relaxation of Clay-Bound Water in Organic-Rich Shales with Implications for NMR Logging

　　

富含有机质页岩含有化学非均质性矿物—有机界面，会产生强且空间可变的质子表面弛豫，使黏土束缚水识别以及核磁共振（NMR）弛豫时间向孔径转换变得困难。传统解释常把表面弛豫率视为常数，而保水界面化学、分子受限和顺磁中心各自贡献尚未被充分分离。研究人员以永川区块五峰组—龙马溪组页岩为对象，采用矿物学与孔隙结构表征、傅里叶变换红外光谱（FTIR）、X射线光电子能谱（XPS）、阳离子交换容量（CEC）、ζ电位、电子顺磁共振（EPR）、受控水化、一维和二维低场时域NMR，以及分子动力学（MD）模拟开展研究。在本研究流体与采集条件下，强界面关联水操作上主要界定为T2 < 1.6 ms、T1 < 85 ms，有效横向表面弛豫率（ρ2,eff）为3.6—9.1 μm/s。保水更密切关联CEC与表面–OH/O–C环境，而弛豫效率更多由EPR可检测顺磁中心和界面水受限分子运动控制。对Na-蒙脱石、绿泥石、伊利石及有机质覆盖伊利石的模拟显示，水密度层化、吸附强度、氢键持续性、驻留、平动扩散、转动重取向和质子—质子偶极关联存在系统差异。将EPR所得顺磁中心密度与表面积加权分子受限指数结合的理化模型解释了86%的实测弛豫率变化，调整R2为0.83，留一交叉验证RMSE与MAE分别为0.71和0.57 μm/s。岩心尺度上，变弛豫率解释把特征孔径平均绝对百分比误差从21.9%降至4.5%，黏土束缚水分数RMSE从3.4降至0.4个百分点。转入NMR测井后降低了岩性相关偏差。结果表明页岩表面弛豫率是磁性与分子耦合控制产生的突现界面属性，为化学非均质页岩NMR分析提供机理基础。

论文发表在《Molecules》。富含有机质页岩的低场时域核磁共振（NMR）可非破坏性评价含氢组分与孔隙流体分布，但短T₂信号可能同时来自黏土束缚水、毛细管约束水、可动水、烃类、干酪根和固体沥青。传统方法常把横向表面弛豫率（ρ₂）当作常数，导致相同弛豫时间在不同地层不一定对应相同孔径或水态。已有研究证明Fe承载顺磁矿物、黏土含量、内部磁场梯度都会影响弛豫，但 bulk 黏土或黄铁矿含量不能单独区分可接近顺磁位点、表面电荷、含氧官能团、分子限域和内部磁场梯度的贡献，因此需要从实验和分子层面直接分离这些控制因素。

研究人员以四川盆地永川区块五峰组—下龙马溪组12块岩心（YC-01至YC-12，深度约3450—3550 m）为样本队列，开展矿物、孔隙、界面化学、顺磁中心和NMR弛豫的综合研究，目的是用受控水化而非通用截断值识别强界面关联水，并分离保水化学、顺磁中心与分子受限对NMR弛豫的作用，建立理化知情的变表面弛豫率模型并推广到测井解释。

主要关键技术方法包括：全岩与黏土X射线衍射（XRD）用于矿物组成；低温N₂/CO₂吸附用于孔隙结构与比表面；衰减全反射傅里叶变换红外光谱（ATR-FTIR）与X射线光电子能谱（XPS）用于表面官能团和氧环境；阳离子交换容量（CEC）与ζ电位用于电荷相关保水性质；X波段电子顺磁共振（EPR）用于顺磁中心密度；12 MHz低场时域NMR一维T₂与二维T₁–T₂用于水态与弛豫；经典分子动力学（MD）模拟Na-蒙脱石、绿泥石、伊利石、干酪根覆盖伊利石界面；再用EPR顺磁中心密度与MD分子受限指数建立化学知情表面弛豫率模型并做留一交叉验证和测井盲验。

2.1 矿物学、孔隙结构与界面化学特征

2.1.1 矿物与有机质组成：XRD显示样品以石英和黏土为主，伊利石占黏土分数56—66%，TOC为2.44—5.42 wt%，R_o,eq为2.58—2.81%，成熟度集中。研究人员据此认为样品成分差异主要来自黏土类型、有机质覆盖和Fe相，而非热演化。

2.1.2 孔隙结构与比表面特征：N₂等温线为IV型伴H3—H4迟滞，介孔主峰3.4—4.6 nm，CO₂吸附解析0.55—0.75 nm微孔。孔体积与总黏土含量正相关，BET比表面与TOC中等相关，说明有机微孔贡献比表面、黏土孔系统贡献孔体积，但气体可接近表面不等于盐水可接近表面。

2.1.3 表面官能团、电荷性质与顺磁中心：ATR-FTIR与O 1s XPS显示表面–OH/O–C组分占24.6—35.6%；CEC为11.7—16.1 cmol(+)/kg，ζ电位为?22.1至?28.9 mV；EPR给出g=2.003有机自由基、g≈2.0 Fe相关信号和g≈4.3高自旋Fe³⁺，有效顺磁中心密度为2.7×10¹⁸—5.8×10¹⁸ spins/g。研究人员得出结论：CEC、含氧环境与ζ电位主要控制保水，EPR中心控制磁性弛豫增强。

2.2 黏土束缚水的NMR弛豫行为：受控水化下，33%相对湿度为单短T₂组分，97%变双峰，真空饱和后出现更长组分。强界面关联水操作边界为T₂=0.2—1.6 ms、T₁=12—85 ms。短T₂几何均值与EPR顺磁中心密度强负相关（r=?0.86），97%吸水量与CEC强正相关（r=0.84）。结论是保水量由带电羟基表面决定，弛豫效率由EPR中心决定。

2.3 界面水的分子结构与动力学

2.3.1 界面水分布与吸附构型：MD显示第一水密度峰位于距表面0.24—0.31 nm，水化层厚度Na-蒙脱石0.94 nm > 绿泥石0.82 nm > 伊利石0.76 nm > 干酪根覆盖伊利石0.58 nm；吸附能依次为?46.8、?40.6、?34.7、?24.9 kJ/mol。干酪根覆盖掩盖亲水界面，不单纯增加吸附面。

2.3.2 氢键网络与分子驻留动力学：水—表面氢键密度Na-蒙脱石8.3 nm^?2最高，干酪根覆盖伊利石降至3.9 nm^?2；氢键寿命与首层驻留时间同序。研究人员说明短T₂水是动态交换集合，不是固定单分子层。

2.3.3 平动—转动动力学与NMR弛豫分子起源：界面平动扩散和二阶转动关联时间均被抑制，质子—质子偶极关联积分相对体相水分别增至3.8、3.2、2.6、1.7倍。MD只给分子受限对偶极弛豫的本征贡献，顺磁增强由EPR独立约束。

2.4 有效表面弛豫率的化学控制：ρ_2,eff为3.6—9.1 μm/s，与EPR顺磁中心密度正相关（r=0.87），与表面积加权分子受限指数I_MD正相关；BET面积和总黏土量单独不显著。两预测变量标准化模型解释86%变异，调整R²=0.83，留一验证RMSE=0.71 μm/s、MAE=0.57 μm/s。结论是ρ₂不是页岩或单类黏土常数，而是条件依赖的界面参数。

2.5 化学知情的NMR测井解释与验证

2.5.1 岩心尺度与独立参考比对：固定ρ₂=6.2 μm/s时特征孔径平均绝对百分比误差21.9%，矿物校正10.9%，化学知情模型4.5%；黏土束缚水分数RMSE由3.4降至0.4个百分点。原因是矿物丰度不能区分保水环境与顺磁环境。

2.5.2 现场NMR测井应用与不确定性：YC-A井建深度相关ρ_2,chem，YC-B井3610—3670 m盲验显示化学知情预测贴近1:1线并在95%置信区间内。主要不确定来自测井元素Fe在绿泥石、黄铁矿等相间的分配及短T₂检测限转换。

讨论部分总结：研究人员强调表面弛豫率应视为由分子受限、EPR可检测磁环境和界面水—孔心交换共同决定的突现界面属性；固定ρ₂会系统性误判化学非均质页岩的孔径与流体划分。化学知情框架把实验室EPR和MD描述符转为测井代理变量，在不改总振幅前提下降低岩性偏差，但测井代理不是井下直接测得的顺磁中心密度或分子动力学量。

研究结论部分翻译：有效横向表面弛豫率（ρ_2,eff）在3.6至9.1 μm/s之间，强界面关联水操作上主要位于T₂ < 1.6 ms与T₁ < 85 ms。保水更密切关联阳离子交换容量和表面–OH/O–C环境，而弛豫效率更多由电子顺磁共振可检测顺磁中心与界面水受限分子运动控制。Na-蒙脱石、绿泥石、伊利石和干酪根覆盖伊利石的模拟揭示水密度层化、吸附强度、氢键持续性、分子驻留、平动扩散、转动重取向和质子—质子偶极关联存在系统差异。将电子顺磁共振派生顺磁中心密度与表面积加权分子受限指数结合的理化模型解释86%实测弛豫率变异，调整R²为0.83，留一交叉验证RMSE与MAE分别为0.71和0.57 μm/s。岩心尺度上，变弛豫率解释把特征孔径平均绝对百分比误差从21.9%降到4.5%，黏土束缚水分数RMSE从3.4降到0.4个百分点；转入核磁共振测井进一步降低孔径转换与黏土束缚水划分中的岩性相关偏差。这些结果把页岩表面弛豫率定义为由耦合磁性与分子控制产生的突现界面属性，并为化学非均质页岩材料的核磁共振分析提供机理基础。