核磁共振与热损伤

岩芯试样取自甘肃省北山地区，平均海拔1600~1700m，年降雨量78.9mm，年蒸发量近3130.9mm，属于干旱戈壁岩漠地区。试样选用宏观均质性良好的花岗岩岩块，加工制作成φ50mm x 100mm的标准样。

图1 花岗岩岩芯试样

通过对芯样进行XRD衍射分析，其矿物主要成分为60.59%长石，34.09%石英，5.32%黑云母。

图2 花岗岩XRD衍射分析

以4℃/min的升温速率加热至预设温度（100℃、200℃、300℃、400℃、500℃、600℃、700℃），保持恒温2h，之后恒定降温速率对样品冷却，每组5个试样，以减少实验误差。

采用纽迈分析的MacroMR12-150H-I大口径核磁共振成像分析仪，对热损伤后的岩芯样品进行了核磁共振测试。

图3 MacroMR12-150H-I 低场核磁共振成像分析仪

图4 不同处理温度下花岗岩样品的弛豫时间T2

从不同温度下的T2谱图中可以看出，T2图谱表现为双峰特性，其中较低的峰值范围主要在10ms以内，大于10ms以上的孔径范围信号幅值偏大，即孔隙主要集中在该范围内。低于400℃以下的图谱双峰形态和走势一致，谱曲线基本重合，表明加热温度低于400℃时，温度对花岗岩内部的微观结构影响极小。

当温度从500℃逐渐升级至700℃时，右侧峰曲线所对应的T2谱范围扩大，左峰和右峰信号幅值显著增强，谱范围持续扩大，表明持续的热损伤加剧导致了岩芯内部微裂纹、穿晶裂纹的出现，使得裂纹的长度和数量不断增加，同时随着温度的升高，较小的裂纹也会逐渐演变为大的裂纹。

图5 T2 谱面积与温度的关系

利用T2谱面积来评价温度作用下岩石损伤程度，从上图可以看出，T2谱面积在低于400℃时，随温度的升高，谱面积并未发生明显改变，当温度高于500℃时，谱面积急剧增大，整体谱面积与温度呈现幂函数形式的增加。

图6 孔隙率与温度的关系

从上图可以看出，核磁共振孔隙度与温度之间的变化规律，与T2谱面积变化规律相似，在高于500℃时，内部结构变化显著，均呈现为幂律函数形式增加。

当像素颜色点颜色越深，孔隙、裂隙所充满的液体也就越多，表明样品内部孔隙、裂隙愈发育，对应像素值愈大，从而可直观反映岩样内部孔、裂隙分布情况。

图7 不同温度处理样品水平截面的质子密度加权像

上图为不同温度下样品的水平截面质子密度加权图像，右侧图例彩条显示内部孔隙含1H质子的相对强度范围。当温度低于500℃时，质子密度分布均匀，从像素分布看，没有明显的质子密度簇；随着温度继续升高，质子密度高的区域逐渐融合成大的连通区域，微裂纹相互作用并聚结成裂纹网。

核磁共振成像像素值大小及像素点深浅时反应损伤程度的介质，图像由不同颜色的像素点组合而成，利用像素值的差异进行岩芯成像的统计分布，通过图像的统计分析是研究岩石受热过程中内部结构变化规律的有效手段。

图8 不同花岗岩试样NMRI像素值概率密度函数分布

上图为不同温度花岗岩试样像素值的概率密度函数分布图，所有对应温度下的成像分布遵循对数正态分布。随温度升高，曲线向右大幅度移动，较大像素值增加，岩石内部损伤加剧，概率密度峰值所对应的像素值增大并向右转移，像素分布迁移规律与T2谱图分布变化特征一致，能定量表示花岗岩在高温作用下内部结构损伤过程。

（1）核磁共振T2谱图在低于400℃以下时没有明显变化，当温度高于500℃时，T2图谱幅值显著增加且向右大幅移动，T2谱面积与孔隙率在升温过程中呈现幂率关系。

（2）通过核磁共振成像可以看到，温度低于500℃时质子密度分布均匀且没有发现明显的质子密度簇，说明岩石内部结构稳定；当温度持续升高，芯样内部晶体裂纹和边界裂纹逐渐产生，融合形成大的连通区域。

（3）花岗岩岩样在不同温度条件下的核磁成像像素概率密度函数都服从对数正态分布，当温度超过500℃时，概率密度函数整体向右移转。