基于CT图像处理的冻结岩石细观结构及损伤力学特性研究

摘要

低温下冻结岩石的损伤是影响寒区岩石工程、地下低温贮存库以及冻结法施工的矿井建设等工程安全稳定的重要因素。论文以寒区岩石工程和人工冻结岩石工程为背景,针对低温状态下的冻结岩石损伤问题,以岩石的细观结构为切入点,以损伤力学、细观力学、体视学、冰力学作为理论依据,以CT扫描实验为研究基础,采用CT图像处理技术、损伤理论分析和细观数值计算相结合的方法,对低温环境下冻结岩石的细观结构及损伤力学特性进行研究,主要内容和成果如下:
　　(1)进行了陕西红砂岩和灰砂岩两种岩石在常温、-2℃、-5℃、-10℃、-20℃、-30℃的CT扫描实验,获得了不同低温下岩石CT扫描图像。冻结岩石CT图像很好地反映了冻结过程中岩石的细观结构、各相组分分布(水、冰、岩石矿物颗粒)及损伤变化情况。低温作用下岩石内部产生新的孔隙和裂隙,同时伴有裂隙、孔隙贯通现象。岩石内部初始损伤随着温度的降低扩展增大,当温度降低至-20℃后,岩石内部细观结构不再变化。
　　(2)将数字图像处理理论应用于冻结岩石CT图像处理中,实现了CT图像的伪彩色增强;获得了不同冻结温度下岩石的CT数直方图;将最大类间方差法和遗传算法相结合,对冻结岩石CT图像进行双阈值分割,实现了冻结岩石CT图像的三值化分割技术;运用Canny算子对冻结岩石CT图像进行边缘检测,获得了冻结岩石细观结构的二值图像。冻结岩石CT图像的伪彩色增强能够提高CT图像的分辨率,减少视觉上的判断误差,根据伪彩色增强图像颜色的变化、CT数直方图可定量分析冻结过程中岩石细观结构、水冰含量以及损伤分布随温度的变化情况。冻结岩石CT图像的三值化分割技术,将冻结岩石内部的水、冰、岩石三相介质区分开,明确地给出了水、冰、岩石空间位置及含量的数字表述。
　　(3)将体视学原理引入到冻结岩石细观结构特性研究中,结合冻结岩石CT图像,分析了冻结岩石细观结构的构成及其相互关系,给出了冻结岩石细观结构参数计算公式。对冻结过程中岩石内部裂隙、孔隙的长度、周长、面积、宽度、圆形度等细观结构参数及水冰含量进行定量计算。细观结构参数能够显示冻结过程中损伤大小及形态的变化情况。冻结岩石内部水和冰含量定量计算结果表明:岩石内部水、冰含量随温度的变化而改变,两者处于动态平衡状态中。
　　引入未冻水含量作为内变量,从与内变量功共轭的相变潜热出发,利用连续介质热力学理论推导出冻结岩石中未冻水、冰含量与冻结温度关系的理论公式。所给出的公式能描述冻结过程中岩石内部未冻水和冰含量的变化情况,同时可预测给定含水量的冻结岩石的纯冰点。
　　(4)根据冻结岩石的细观结构组构及特性,从岩石细观力学机理出发,用混合律方法将含孔隙的冻结岩石视为各向同性介质,将含裂隙的冻结岩石视为横观各向同性介质,推导出不同冻结温度下的等效弹性模量计算公式。
　　针对冻结作用对岩石力学性能的影响,提出“冻结负损伤”的概念,用冻结损伤、荷载损伤描述冻结和荷载两种不同作用下岩石损伤过程,拓展损伤变量内涵。基于连续介质热力学方法,以岩石的初始损伤状态为基准状态,应用推广后的应变等价原理及损伤力学理论,建立荷载作用下冻结岩石的宏—细观损伤本构模型。对模型进行了验证,理论曲线与实验曲线比较接近,所建立的模型能够描述荷载作用下冻结岩石损伤演化规律。
　　(5)提出了冻结岩石数字图像数值分析方法(DIP-FEM),建立不同低温环境、不同荷载作用下岩石破坏过程的有限元数值计算格式。通过CT图像的三维重建技术,获得冻结岩石真实细观结构数字化模型,对其进行矢量化转化,并导入至有限元软件ANSYS中,进行冻结过程中温度场及冰膨胀力分布规律的数值模拟试验研究,分析含裂隙或孔隙岩石在低温环境下的损伤特性,该方法弥补了低温冻结岩石物理实验的不足,为实现冻结岩石的破坏过程研究从细观尺度向宏观尺度过渡提供新的研究途径。

目录

声明

中文摘要

英文摘要

目录

1 绪论

1.1 研究背景及意义

1.2 国内外研究现状及发展趋势

1.3 主要研究内容和技术路线

2 冻结岩石细观结构的CT扫描实验

2.1 CT扫描原理及冻结岩石CT图像

2.2 低温饱和冻结岩石CT扫描实验

2.3本章小结

3 冻结岩石CT图像处理技术

3.1 冻结岩石CT图像的伪彩色增强技术

3.2 冻结岩石的CT数直方图技术

3.3 冻结岩石CT图像的三值化分割技术

3.4 冻结岩石CT图像的边缘检测技术

3.5 本章小结

4 冻结过程中岩石细观结构参数及水冰含量分析

4.1 体视学原理

4.2 基于体视学原理的冻结过程岩石细观结构参数定量计算

4.3 冻结过程岩石内水冰含量理论表达及验证

4.4 本章小结

5 荷载作用下冻结岩石宏—细观损伤力学特性研究

5.1 理论基础

5.2 荷载作用下冻结岩石宏—细观损伤本构模型建立

5.3 模型验证

5.4 本章小结

6 基于数字图像数值分析方法(DIP-FEM)的冻结岩石损伤特性

6.1 冻结岩石CT图像的三维重建

6.2 冻结岩石数字图像数值分析方法（DIP-FEM）的实现过程

6.3 基于DIP-FEM方法的冻结岩石损伤特性分析

6.4本章小结

7 结论及展望

7.1 结论

7.2 主要创新点

7.3 展望

致谢

参考文献

附 录