在高放废物地质处置、地热岩石的开发与利用、隧道、洞库、涵洞等经受火灾后修复等工程中，都涉及到高温岩石问题。该类工程除了承受常规静荷载外，也可能面临地震、爆炸、武器冲击以及车辆荷载等多种动力荷载作用^[1]。这需要我们加深对高温前后岩石静力学与动力学性能的了解，为此国内外学者进行了大量的研究。在静力学方面，徐小丽等^[2]、Chen等^[3]、Liu等^[4]、Dwivedi等^[5]研究了不同温度热处理后，花岗岩单轴压缩下物理力学性质。蔡燕燕等^[6]研究了高温后花岗岩在卸荷路径下的力学性质。在动力学方面，刘石等^[7]、李明等^[8]、陈腾飞等^[9]、尹土兵等^[10-11]分别对高温后大理岩和砂岩进行了冲击动力实验，获得了岩石高温后动力学性能参数。刘石等^[12]、支乐鹏等^[13]对经温度处理后的花岗岩进行了波动特性与动力特性的研究。卢志堂等^[14]探讨了应变率对不同温度下花岗岩动力学参数的影响。

虽然以上实验研究考虑了温度对岩石的影响，但主要是对高温岩石自然冷却后再进行实验。然而在实际工程中往往会出现高温岩石遇水冷却的情形，如对于遭受火灾的隧道、洞库等地下工程，采取喷水灭火；高温岩石施工过程中利用钻井液或泥浆对岩石冷却；高地温地质条件下先借助水冷却再开挖修建隧道等等。水的急剧冷却将导致围岩温度迅速降低，势必产生热破裂，这显然不利于围岩的稳定性。该问题在近几年逐渐受到重视，黄真萍等^[15-16]对高温大理岩、石灰岩遇水冷却后的力学、声学以及波动性能进行了研究；郤保平等^[17]发现600 ℃以内花岗岩遇水冷却后，其单轴抗压强度、抗拉强度以及弹性模量均随着温度增大而减小；朱振南等^[18]研究了500 ℃内高温花岗岩遇水冷却后的物理力学特性；梁铭等^[19]通过实验得到花岗岩经自然冷却与遇水冷却后的温度和孔隙率变化曲线。以上研究主要是从静力学角度讨论水冷却对高温岩石的影响，然而经水冷却后的高温岩石在后期仍面临着地震、爆炸、施工振动等动力影响，所以从岩体工程的安全性、可靠性、耐久性等角度出发，必须了解高温岩石经水冷却后的动力学性能。但目前只有翟越等^[20]研究了高温混凝土遇水冷却后的动力学性能，对岩石材料的实验研究尚未见报道。鉴于此，本文中选择洞库、隧道工程中常见的花岗岩，对其加热至预定温度后，采取水冷却，随后进行波速测试、核磁共振扫描、分离式霍普金森压杆冲击实验和扫描电镜观测，深入研究水冷却对高温岩石的细观结构损伤及动态力学性能影响，并与自然冷却结果对比，得到一些有益成果，为高温后岩体工程的研究提供借鉴与参考。

1.   岩样制备及实验方案

1.1   岩样制备

岩样取自湖南岳阳完整花岗岩大样，岩石状态较好，致密无裂纹。根据国际岩石实验试样尺寸的标准，按照2∶1比例，切割加工成直径50 mm、高25 mm的圆柱体，并对其两端进行仔细打磨，保证其表面平整度在0.02 mm以内。加热前试样平均密度2.61 g/cm³，纵波波速在3 859～4 074 m/s之间，试样均匀性与一致性较好。

1.2   岩样处理与实验方案

为了研究高温花岗岩经水冷却后的动态力学性能，设定5组温度条件：25 ℃（常温）、200 ℃、400 ℃、600 ℃和800 ℃。作为对比，加热后岩石分别采取浸水冷却与自然冷却两种方式，每组温度每种冷却方式各3块试样（常温组3块），共计27块。

采用ZBL-U510非金属超声检测仪对高温前后试样进行波速测试，获得波速变化值。根据编号将每组6块试样同时放入型号为SX2-2.5-12N的箱式电阻炉内加热到预定温度并恒温2 h后，取出3块试样放入预先注满水的容积为20 L的铁皮水桶中进行浸水冷却。由于600 ℃与800 ℃时，岩石温度较高，为保证冷却效果，将高温岩石放入水桶中后，随即再往桶中注入水，将升温后的热水置换出来，以保证水冷却效果，余下3块留在电阻炉中自然冷却。对冷却后的岩石试样采用MesoMR23−060H−I型核磁共振成像分析仪进行核磁孔隙度测试，随后进行动力冲击实验，最后对冲击破碎后的试样进行扫描电镜观测，分析温度、冷却方式与动力冲击耦合作用下花岗岩损伤特征。由于本次实验不考虑应变率效应，为保证冲击速度一致，系统气压统一设为0.4 MPa，每次实验前用直尺将子弹送至同一位置，最终获得每组冲击速率平均值为6.25～6.57 m/s。为获得可靠的应力应变关系，根据前人研究，采用波形整形技术，在入射杆端部黏贴紫铜波形整形器，并在试样与杆件接触部位涂抹凡士林以减小摩擦。图1为实验获的得典型SHPB实验应力波形。

图  1  高温花岗岩典型SHPB应力波形

Figure  1.  Typical SHPB stress waveforms for high-temperature granite specimens

下载: 全尺寸图片 幻灯片

从图1中可以看出入射波缓慢上升，应力波在传播过程中无明显震荡，表明较好满足应力波一维传播假设。

2.   超声波及孔隙度结果分析

在测试前，先用游标卡尺测量出各个试样的高度，并在其表面涂抹凡士林作为耦合剂，随后在主机界面输入高度尺寸，即可直接读取波速结果，每个试样读取6个点。

2.1   超声波形图分析

根据超声波形图实验结果发现25～400 ℃热损伤较小，超声波形状基本相似，仅幅值大小不同，800 ℃时高温引起损伤已经很大，冷却方式影响不明显，故限于篇幅，此处只列出部分超声波形图。

从图2中可以看出，从室温到600 ℃自然冷却，随着温度升高，超声波幅降低，超声波形比较规则，波峰间距均衡，表明此时花岗岩内部虽然产生新的损伤，但岩样仍能维持一定均匀性。600 ℃花岗岩经水冷却后超声波形不规则，随着温度升高到800 ℃，超声波形更加凌乱，表明岩石内部萌生更多孔隙，试样均匀性变差，声波传播到裂隙部位将进行多次反射、折射等，使超声波形变得凌乱。

图  2  高温花岗岩不同温度后纵波波形图

Figure  2.  Longitudinal waveforms of granite sample exposed to different temperatures

下载: 全尺寸图片 幻灯片

2.2   波速变化分析

超声波在穿过不同介质时，波速会发生变化，若内部存在裂缝、空洞等缺陷则超声波速会降低。图3所示为岩石高温前后超声波速对比情况，图中数值为三个试样均值。

图  3  高温前后纵波波速均值对比

Figure  3.  Comparison of average longitudinal wave velocity before and after exposure to high temperature

下载: 全尺寸图片 幻灯片

从图3中可以看出，随着温度升高，水冷却处理后的高温花岗岩超声波速分别下降到初始超声波速平均值的88%、81%、68%和29%，整体表现为非线性下降趋势。从图中还可以发现高温花岗岩经水冷却后的超声波速均比自然冷却后低，在600 ℃时超声波速差达到最大值(344 m/s)，而在800 ℃时超声波速差最小，仅为36 m/s。分析原因是水冷却使花岗岩表面温度急剧下降，试样内部将形成较大温度梯度，从而产生二次温度应力引起更大热损伤，加剧高温花岗岩的劣化。但在800 ℃时，由于花岗岩内部矿物在高温作用下已经充分膨胀，产生较多微孔隙或原生孔隙进一步扩大、贯通，使得力学性能大幅降低。水冷却虽能进一步引起岩石力学性能劣化，但空间不大，故二者差异不明显。

2.3   孔隙度

由于岩石结构的骨架不会产生核磁信号，故需要先对岩样进行饱和处理，让水浸入岩石内部，利用水的原子核被极化产生的横向磁化矢量衰减时间（即${T_2}$弛豫时间）与饱和岩样孔隙体积存在的对应关系，通过对${T_2}$谱总面积进行标准刻度转换，从而计算出岩石的孔隙度和孔径分布。所以，核磁共振检测得到的孔隙度实质上是岩样内部的孔隙与裂隙之和，孔隙度的变化也可反应岩样内裂隙的变化情况，二者具有一致性^[21-23]，故接下来均从孔隙角度进行讨论。表1为对应的${T_2}$谱面积，图4所示为4个温度等级下，高温花岗岩冷却后孔径分布图，从中可以看出：

表  1  核磁共振T₂谱面积

Table  1.  T₂ spectrum area of NMR

冷却方式	T2 谱面积
	25 ℃	200 ℃	400 ℃	600 ℃	800 ℃
自然冷却	1 852	1 910	2 086	3 072	4 029
水冷却	1 852	1 981	2 254	3 544	4 887

下载: 导出CSV 

| 显示表格

图  4  核磁共振孔径分布及断面扫描图

Figure  4.  Pore size distribution curves and scanning images of cross-sectional area

下载: 全尺寸图片 幻灯片

（1）原始岩样的${T_2}$谱曲线仅有一谱峰，对应孔径大小为0.01 μm左右，表明原始岩样中的孔隙以小孔径孔隙为主，大孔径孔隙为辅。高温花岗岩经水冷却处理后，${T_2}$谱面积随着温度增大而增大，${T_2}$谱曲线出现两个谱峰，第一个谱峰对应孔径仍然为0.01 μm左右，第二个谱峰对应孔径则随着温度增加而变大，在800 ℃时最大孔径达到16 μm。表明温度引起热损伤，使得试样孔径扩张，孔径大小与孔隙数量增加。

（2）小孔径孔隙数量随着温度增大呈现先减少后增多的变化规律。在200 ℃时冷却后的试样内小孔径孔隙数量大小关系为：自然冷却＜水冷却＜原始岩样，在其他温度条件下则为：原始岩样＜自然冷却＜水冷却。这是因为在200 ℃时温度对岩石的主要影响是使其所含矿物体积发生膨胀，从而使试样内原生微孔隙闭合，导致小孔隙减少。而水冷却加速了高温岩石冷却，减少了矿物膨胀时间，所以此时小孔径孔隙数量大于自然冷却。当温度升高到400 ℃后，矿物间不均匀热膨胀趋于显著，由此引起的温度应力将会超过岩石强度极限，因而萌生出比原始岩样更多的小孔径孔隙。

（3）水冷却后试样大孔径尺寸和数量随着加热温度的增大而增大，且均大于同等温度条件下自然冷却的结果，两者之间差异随着温度增大而增大，但800 ℃时水冷却与自然冷却后试样的大孔径基本相同。这是因为随着温度升高，试样内微孔隙不断发展、贯通形成大孔隙。水的急剧冷却产生显著的温度梯度，因而在岩样内产生二次温度应力，诱发出更多孔隙。这种急剧冷却效果与温度直接相关，温度越高效果越明显。当温度达到800 ℃时，由于加热到预定温度后保温2 h，故孔隙得到充分发展，水的急剧冷却主要起到萌生孔隙的作用，所以在裂隙数量上明显多于自然冷却，而孔径大小基本不变。

3.   SHPB实验结果及分析

3.1   应力应变特征

图5所示为高温花岗岩采用水冷却后的动态应力应变曲线，作为对比同步给出自然冷却后的实验结果。从中可以看出，采用水冷却处理后，随着温度的增大，花岗岩动态峰值应力减小，动态峰值应变增加。与自然冷却结果对比可发现，在各个温度条件下，水冷却后花岗岩的动态峰值应力均小于自然冷却，应力应变曲线峰后应力跌落迅速，表现出一定脆性破坏特征。

图  5  花岗岩动态应力应变曲线图

Figure  5.  Dynamic stress-strain curves of granite

下载: 全尺寸图片 幻灯片

从图5中还可以发现，自然冷却后花岗岩应力应变曲线在低于200 ℃时压密段不明显，随着温度继续升高，压密段越加明显。这主要是因为当温度不高时，不同矿物质热膨胀系数引起矿物间不均匀热膨胀变形较小，此时温度主要引起花岗岩内部空气以及部分矿物质膨胀，使得一些初始微裂缝闭合，提高试样的密实性。而随着温度继续升高，矿物间不均匀膨胀越加显著，由此产生的热应力超过花岗岩抗拉强度时则会诱发新裂缝，导致内部孔隙增多，且矿物内结晶水的逃逸也使得矿物的晶格结构发生破坏而造成内部缺陷增大并联通，所以表现出明显的压密段。采用水冷却的应力应变曲线不管在何种温度条件下，均表现较强的压密段。分析其原因是水冷却使高温岩石表面温度急剧下降，产生显著温度梯度，相比较自然冷却而言则会引起额外的热应力，加剧岩石的损伤，产生更多的孔隙（核磁扫描结果可证明这一点）。所以水冷却花岗岩最终表现出来为相对低的动态峰值应力和显著的压密段。

3.2   动态峰值应力与动态峰值应变

图6所示为高温花岗岩经水冷却后动态峰值应力${\sigma _{\rm{VW}}}$ 随温度变化，以及与自然冷却结果 ${\sigma _{\rm{VN}}}$ 对比情况。通过对数据进行拟合，发现两种方式冷却后花岗岩的动态峰值应力均随着温度增大而呈线性下降趋势，且水冷却后动态峰值应力下降速度相对更为剧烈。

图  6  花岗岩动态峰值应力

Figure  6.  Dynamic peak stress of granite

下载: 全尺寸图片 幻灯片

从图6中还可以发现，与自然冷却结果相比，各个温度条件下高温花岗岩经水冷却后的动态峰值应力均低于自然冷却的结果，从200 ℃到800 ℃，强度分别降低29、20、42、10 MPa，降低幅度达到了13.4%、10.2%、24.1%、11.9%。由此可见，水冷却处理对高温花岗岩的峰值应力具有明显的劣化作用，会进一步降低花岗岩的强度。这也可以解释部分岩质工程在火灾过程中未倒塌破坏，但经射水灭火降温后发生倒塌的现象。

图7所示为花岗岩动态峰值应变与温度关系。从中可以看到经水冷却后的高温花岗岩动态峰值应变随着温度升高而增大。与自然冷却结果对比可以发现在200 ℃时水冷却结果大于自然冷却，随后则小于自然冷却。这是因为水冷却对高温岩石有两方面影响：一是扩大温度差，产生温度应力，加速孔隙萌生与发展；二是具有冷淬作用提高试样硬度。前者导致花岗岩致密程度降低，在冲击荷载作用下，动态峰值应变显著增大，而后者则会降低花岗岩塑性变形能力，引起动态峰值应变减小，二者综合作用使得水冷却后花岗岩峰值应变表现出与自然冷却不同的规律。冷淬作用跟温度有密切关系，在200 ℃时冷淬作用不显著，此时水冷却的影响主要体现在引发温度应力，加速孔隙萌生，起到劣化岩石的作用，所以动态峰值应变大于自然冷却。而从400 ℃开始，水的冷淬作用愈加显著，在一定程度上提高了岩样的硬度，降低了塑性变形能力，故在冲击荷载作用下表现出较自然冷却小的动态峰值应变。

图  7  花岗岩动态峰值应变

Figure  7.  Dynamic peak strain of granite

下载: 全尺寸图片 幻灯片

3.3   弹性模量

由于高温后岩石动态应力应变曲线均存在一定长度压密段，所以为保证结果具有对比性，本文中采用弹性模量E_C^[7]，即选取应力应变曲线上升段应力为40%和60%峰值应力的两点连线的斜率，结果如图8所示。

图  8  弹性模量随温度变化曲线

Figure  8.  Variation of dynamic elastic modulus EC with temperature

下载: 全尺寸图片 幻灯片

从图中发现水冷却后花岗岩弹性模量E_CW先增大后减小，整体上呈非线性下降趋势，拟合结果如下：

同时还可以看到从400 ℃开始，采用水冷却处理的花岗岩的弹性模量大于自然冷却后的结果E_CN，其原因与前文所述相同，是因为水的冷淬作用提高了花岗岩硬度。

4.   冲击破坏形态及损伤机理

4.1   冲击破坏形态

图9所示为高温花岗岩经两种方式冷却、冲击后的破坏形态对比图（编号第1个数字表示冷却方式（1：自然冷却、2：水冷却），第2个数字表示温度等级、第3个数字表示块号）。从图中可以看出花岗岩冲击破坏程度随温度的增大而剧烈，由室温时的块状逐步变成800 ℃时的粉末状。在200～400 ℃，试样破坏均表现出轴向劈裂现象，随着温度的升高，试样周边破坏加剧，水冷却方式对高温花岗岩破坏形态影响不明显。但在600 ℃时，水冷却后花岗岩出现了轴向破裂面，四周破裂以块为主，而自然冷却轴向破裂面不明显，外围破坏出现较多粉末。800 ℃时自然冷却的花岗岩粉碎更为彻底，水冷却则存在相当数量的大颗粒状碎块，通过测量发现自然冷却最大颗粒尺寸（长×宽×深）为：2.1 cm×1.2 cm×0.8 cm；水冷却最大颗粒尺寸：2.5 cm×1.4 cm×1.2 cm。

图  9  花岗岩冲击破坏形态

Figure  9.  Failure forms of granite under different temperatures and different cooling methods

下载: 全尺寸图片 幻灯片

4.2   水冷却热损伤机理

花岗岩内部含有石英（8.1×10⁻⁶ K⁻¹）、斜长石（3.5×10⁻⁶ K⁻¹）、方解石（13.1×10⁻⁶ K⁻¹）等多种矿物，热膨胀系数存在显著差异，当温度达到一定阈值时，矿物间将发生不均匀热膨胀。石英在573 ℃会发生从$\alpha$相到$\beta$相位移型相变，在此过程中石英体积将增大5%^[24]，进而导致试样内裂纹增多。矿物不均匀膨胀与石英相变综合作用引起微观层面上晶粒面滑移，产生微裂缝。同时岩石内部原有孔隙进一步发育、贯通导致动态力学性能下降。

此外，将高温岩石试样视作整体，采用不同冷却方式，则会得到不同温度场与温度梯度，由此引起的不协调变形将产生二次温度应力。为从理论上探讨水冷却对高温岩石损伤机理，将岩样内部温度应力的力学模型分为热膨胀应力与二次温度应力两部分。把岩石试样划分成若干个有限区域，并假设每个区域内仅有两种矿物，则各区域温度应力计算如下：

其中矿物间不均匀变形所引起的热膨胀应力^[17]：

水冷却引起的二次温度应力：

式中：${\alpha _{ij}}$、${\alpha _1}$、${\alpha _2}$为$ij$区域内岩石、矿物1和2的热膨胀系数；${E_{ij}}$、${E_1}$、${E_2}$为$ij$区域内岩石、矿物1和2的弹性模量；${\nu _{ij}}$为$ij$区域内材料的泊松比；$\Delta {T_{1ij}}$为$ij$区域从室温加热到某温度的温差；$\Delta {T_{2ij}}$为$ij$区域与相邻区域的的温度梯度。

由式（3）～（5）可知，矿物间不均匀膨胀变形所引起的温度应力将使试样内部萌生孔隙，而水冷却引起的二次温度应力则进一步导致已有孔隙尖端应力集中，有利于裂纹的扩展，导致岩石内部的破坏更加严重。同时水的入侵又削弱矿物颗粒间的联系，引起力学性能的进一步降低，所以水冷却后高温花岗岩无论在孔隙数量与孔隙孔径尺寸上均比自然冷却的大。核磁共振扫描结果即可验证这一结论。图10所示为高温花岗岩经两种方式冷却后核磁共振扫描孔隙率分布伪彩图，红色表示孔隙。由于低于400 ℃的孔隙率值较小，从图上看差异不明显，故限于篇幅仅给出600 ℃与800 ℃结果。从中可以直观看出水冷却后花岗岩孔隙数量明显多于自然冷却，这表明水冷却会加速高温岩石力学性能的劣化。

图  10  花岗岩样品成像伪彩图

Figure  10.  pseudo-color image of granite sample

下载: 全尺寸图片 幻灯片

由式（5）可知，二次温度应力大小主要受温度梯度的影响。因岩石导热性能较差，当高温岩石遇水冷却，水的温度并不能快速传递到试样内部，所以整个冷却过程为非稳态热传导。为直观对比冷却方式的影响，采用ABAQUS软件，模拟高温花岗岩温度分布情况。由于高温物体冷却实质上是与外界换热过程，受热换系数的控制，水与岩石接触时热换系数从几千到上万不等^[25]（本文取两个值作为对比），而自然环境则一般为几十。故花岗岩相关参数取为：导热系数3.0 W/(m·K)，试样温度600 ℃，比热635 J/(kg·℃)，热换系数${\alpha _{\rm{H}}}$：水冷却1 000 W/(m²·K)和100 W/(m²·K)，自然冷却25 W/(m²·K)。根据路径提取温度结果见图11。

图  11  温度模拟结果

Figure  11.  Simulation results for temperature

下载: 全尺寸图片 幻灯片

从图11中可以看出，高温岩石经水冷却后，试样表面温度急剧降低，尤其在表面与侧面交点A处，温度降低最快，故此处产生的温度梯度最大，从而将产生最大二次温度应力，且为拉应力^[25]。由于岩石的抗拉强度较低，所以水冷却高温花岗岩，当温度差引起的拉应力足够大时，花岗岩表面将产生裂缝。在实验过程中发现，从600 ℃开始，高温花岗岩经水冷却后，率先在边缘出现肉眼可见裂缝，随着温度升高裂缝逐步向中心发展，且边缘出现岩样脱落现象，如图12所示。

图  12  水冷却花岗岩表面裂缝

Figure  12.  Surface cracks in granite specimen after water cooling

下载: 全尺寸图片 幻灯片

4.3   冲击破坏微观裂缝形态

对冲击破碎后的回收试样进行扫描电镜分析，探讨温度、冷却方式以及冲击荷载耦合作用下对裂缝发展的影响，限于篇幅仅展示部分图片。

从图13中可以看出，花岗岩经200 ℃处理后，裂纹小且短。随着温度升高到400 ℃，裂纹沿着一个主方向延长，四周伴随从裂纹并向着主裂纹汇集。在600 ℃和800 ℃时，通过对比可以发现自然冷却后的断面呈蜂窝状，凌乱无规则，周边伴随小碎块，水冷却后的断面则相对整齐，主要是因为温度应力作用使得岩石内部已经萌生大量裂纹，故在冲击荷载作用下，裂纹迅速扩展贯通形成裂缝。而水的冷淬作用在一定程度上提高了表层花岗岩硬度，在冲击荷载作用下，具有一定脆性破碎特征，从而获得相对规则的破裂断面。

图  13  冲击破坏后的SEM形貌

Figure  13.  SEM image after impact damage

下载: 全尺寸图片 幻灯片

5.   结　论

（1）高温花岗岩，不论采用自然冷却还是水冷却，与常温时对比，其波速均呈现非线性下降趋势。采用水冷却花岗岩，在各个温度等级其纵波波速均低于自然冷却。最大纵波波速差异在600 ℃时出现，大小为344 m/s，最小纵波波速差异则发生在800 ℃时，高温岩石经水冷却后纵波波速仅比自然冷却后纵波波速低36 m/s。

（2）核磁共振扫描结果显示，经温度作用后，花岗岩岩样内${T_2}$谱曲线均出现两个峰值。小孔径孔隙在200 ℃时最小，小于原始岩样。水冷却处理后的高温花岗岩试样大孔隙的孔径大小与数量均随着温度增大而增大，且大孔隙的孔径大小与数量均比自然冷却后岩样大。

（3）随着温度升高，水冷却后花岗岩动态峰值应力减小，峰值应变增大，弹性模量先增大后减小。水的急剧降温诱发更多的孔隙，使得水冷却高温花岗岩应力应变曲线上压密段比自然冷却的要显著。水的冷淬作用降低了花岗岩强度、提高了其硬度，因而水冷却后花岗岩的峰值应力与应变小于自然冷却，而弹性模量大于自然冷却。

（4）温度大小和冷却方式对高温花岗岩冲击破碎形态以及断面形态有一定影响，自然冷却后花岗岩冲击破碎以粉末和小颗粒碎块为主，而水冷却则含有一定数量的大颗粒碎块。从600 ℃开始自然冷却后花岗岩冲击断面呈不规则蜂窝状，而水冷却后冲击断面则相对平整。