综述：多孔介质中的磁共振弛豫

《Progress in Nuclear Magnetic Resonance Spectroscopy》：Magnetic resonance relaxation in porous media

　　

研究人员识别了经典磁共振（MR）弛豫解释的局限性；概述了理解多孔介质弛豫的重要进展；通过玻璃微毛细管模型系统验证了Brownstein-Tarr（BT）理论；从单一磁共振测量中估算孔隙尺寸和表面弛豫率；利用BT理论解释了温度依赖的弛豫行为。

1. 引言
多孔介质广泛存在于自然界和工程应用中，其孔隙尺寸跨度从纳米级到微米级，传统表征方法（如压汞法、气体吸附法、扫描电子显微镜和X射线计算机断层扫描）存在破坏性、需干燥样品或分辨率受限等缺点。磁共振（MR）弛豫技术能够非破坏性地、原位地探测饱和流体多孔介质的孔隙结构和动态过程。Bloch方程定义了纵向弛豫时间（T₁）和横向弛豫时间（T₂），而扩散与磁场梯度的耦合会导致MR信号衰减。Brownstein-Tarr（BT）理论通过求解Bloch-Torrey方程，将扩散、表面弛豫率与孔隙几何联系起来，但传统解释多依赖快速交换假设，该假设在很多实际情况中并不成立。

2. 历史概述
2.1 Bloch方程与自旋回波
1946年Bloch和Purcell独立发现磁共振现象，Bloch提出了描述核磁化强度动力学的唯象方程。1950年Hahn引入自旋回波技术，用于消除磁场不均匀性和扩散引起的信号衰减，随后Carr和Purcell改进为多次回波序列，Meiboom和Gill进一步优化得到CPMG序列，成为低场MR弛豫测量中最常用的方法。

2.2 Bloch-Torrey方程
Torrey在Bloch方程中加入了扩散项，得到Bloch-Torrey方程，描述了磁化强度在进动、弛豫和扩散共同作用下的演化。该方程成为脉冲场梯度（PFG）测量扩散、扩散加权成像等技术的理论基础。

2.3 MR表面弛豫机制
Bloembergen-Purcell-Pound（BPP）理论解释了体相弛豫与分子运动的关系。早期实验表明，流体与固体表面的相互作用显著增强弛豫，Korringa-Seevers-Torrey（KST）模型将表面弛豫分为顺磁和非顺磁贡献，顺磁弛豫通常占主导。

2.4 MR弛豫与孔隙结构
Senturia和Robinson通过实验表明，孔隙尺寸控制弛豫时间，较小孔隙对应较短弛豫时间。Brownstein和Tarr进一步推导出多指数弛豫表达式，建立了弛豫时间与孔隙尺寸、表面弛豫率和扩散系数的定量关系。

3. Brownstein-Tarr理论
3.1 Bloch-Torrey方程的解
Brownstein和Tarr忽略了梯度项，采用傅里叶边界条件，将Bloch-Torrey方程转化为标准扩散方程，并将磁化强度表示为多个本征模的叠加，每个本征模对应一个弛豫时间常数。

3.2 自旋弛豫-扩散区间与弛豫时间
BT理论定义了三个弛豫-扩散区间：快速交换区间、中间区间和慢扩散区间。BT数（N_BT）是决定区间的无量纲参数，其定义为ρa/D。在快速交换区间，弛豫时间满足1/T_i = 1/T_i,B + ρ_iS/V。在中间区间，扩散项不可忽略，弛豫时间表达式为1/T_i,n = 1/T_i,B + D·ξ²_i,n/a²。

3.3 弛豫-扩散区间的比较
快速交换区间中，分子扩散足够快，磁化强度在孔隙内保持空间均匀分布；中间区间中，靠近表面的自旋弛豫更快，导致磁化强度空间非均匀，表现为多指数衰减。

3.4 基态和非基态弛豫模式
基态（ground）模式通常主导信号，非基态（nonground）模式在BT数较大时贡献显著。不同孔隙几何（球形、圆柱形、平面形）对应的特征方程不同，非基态模式的相对强度随BT数增大而增加。非基态模式的检测需要高信噪比和适当的正则化参数。

3.5 BT理论的实验支持
利用具有理想圆柱孔几何的玻璃微毛细管（孔径100 μm）进行CPMG测量，T₂分布中观测到基态和第一非基态两个峰，定量计算得到的表面弛豫率约为10 μm/s，非基态弛豫时间与理论预测吻合，验证了中间区间的存在。

4. 应用与最新进展
4.1 快速交换区间中的应用
4.1.1 孔隙尺寸及其分布估算
在快速交换假设下，T₂分布可转换为孔隙尺寸分布，公式为1/T₂ = ρ₂S/V。该方法在岩石物理中广泛使用，但要求样品处于快速交换区间且表面弛豫率已知。

4.1.2 表面弛豫率估算
传统表面弛豫率估算也依赖快速交换假设，需通过独立手段（如压汞法、微CT）获得孔径，再结合T₂计算。对于不在快速交换区间的样品，该假设会显著低估表面弛豫率。

4.1.3 渗透率估算
基于T₂几何均值（T_2,gm）的Schlumberger-Doll Research（SDR）模型给出了渗透率与孔隙度和T₂的经验关系，隐式假设表面弛豫率恒定。

4.2 中间区间中的应用
4.2.1 非基态BT本征模的检测
通过二维T₁-T₂相关测量，在玻璃珠填充柱和天然岩石中检测到非基态本征模，正则化参数的选择对分辨这些弱信号至关重要。

4.2.2 孔隙尺寸和表面弛豫率测量
发展了几种基于中间区间的方法：正演建模拟合CPMG衰减曲线，可同时获得孔径和表面弛豫率；基于内部场扩散衰减（DDIF）的高阶本征模方法；通过温度变化改变扩散系数，利用非线性拟合直接从T₂数据提取孔径和表面弛豫率。

4.2.3 孔隙尺寸分布测量
对于存在孔径分布的样品，通过计算并扣除大孔隙的非基态贡献，重建基态T_2,0分布，进而获得更准确的孔隙尺寸分布。该方法在玻璃珠和砂岩中得到了验证。

4.3 温度依赖的MR弛豫行为
顺磁样品（如钠钙玻璃珠、Berea砂岩）的T₁和T₂随温度升高而降低；抗磁样品（如硼硅玻璃珠、Richemont石灰岩）则相反。该行为可以用BT理论解释：顺磁样品处于中间区间，温度升高增强扩散，导致弛豫加快；抗磁样品处于快速交换区间，表面弛豫率随温度升高而降低，导致弛豫减慢。

4.4 弛豫-扩散区间之间的转换
表面弛豫率的改变（如清洗玻璃珠或水泥水化）可使弛豫行为在区间之间转换。温度变化改变流体扩散系数，同样可引起区间转换，导致某些样品的弛豫时间对温度近乎不敏感。

5. 结论
Brownstein-Tarr理论为解释多孔介质磁共振弛豫提供了统一框架，将分子扩散、表面弛豫率与孔隙几何联系起来。快速交换假设曾长期主导应用，但近年研究揭示了中间区间的重要性。发展更精确的孔隙表征方法需要进一步扩展BT理论，纳入各向异性扩散、内部磁场梯度和多相流体分布等因素，结合先进实验与数值模拟以提升磁共振作为定量探测工具的能力。