岩石CT试验数据处理的方法的研究.docVIP

岩石CT试验数据处理的方法的研究 　　摘 要：CT技术在材料科学实验中的应用使得材料的缺陷和内部结构得以检测，同时CT试验成为岩石材料细观破裂过程的热点研究课题。岩石内部矿物组成、致密性等不均一，各部位密度不同，在CT灰度图像中则表现为各部位灰度值不一。一般针对CT实验数据的处理仅仅停留在对不同扫描图幅中灰度频率直方图的研究。文章通过对扫描图像进行裂纹面积、分维计算、三维重建等多种手段的分析，丰富了CT实验数据处理的方法，为多角度分析岩石应力应变过程提供了思路。 　　关键词：CT试验；裂纹面积；分维计算；三维重建 　　中图分类号：TU45 文献标识码：A 文章编号：1671-2064（2017）07-0093-03 　　1 引言 　　20世纪80年代以来，CT技术在材料科学实验中的应用使得材料的缺陷和内部结构得以检测，X射线CT试验成为岩石材料细观破裂过程的热点研究课题。岩石试样内部矿物组成、致密性等不均一，各部位密度不同，在CT灰度图像中则表现为各部位灰度值不一。同时，加载过程中裂纹的萌生、扩展、合并及贯通等都直接影响着试样内部不同部位密度的大小。而由CT扫描原理可知，被扫描材料对X射线的吸收系数与材料密度成正比，故正比于吸收系数的灰度图像可以表征材料内部的密度变化[1]。因此，越来越多的研究人员通过分析CT灰度图像中的影像特征来探索压缩过程中泥岩裂纹的扩展与演化律。 　　王梦蔚等[2]针对混凝土试样开展CT扫描实验并结合数字图像处理技术进行了宏细观建模研究；李磊[3]基于页岩CT图像的数字图像处理技术，结合流固耦合理论和损伤力学，研究了页岩的损伤演化规律；王彦琪等[4]基于图像检索技术对单轴压缩破坏过程中的CT图像进行了优化处理，得到了较好的效果；刘京红等[5]在对裂纹扩展过程的分析中加入分形维数指标，给出了岩石材料从裂纹萌生、扩展到贯通的细观破损的解释；彭瑞东等[6]从能量耗散及释放角度探讨了层状盐岩破坏的细观机制。虽然国内外学者借助CT实验对不同岩石裂纹的演化过程及特征进行了一系列研究，但对于CT数据的处理分析手段较为单一，并不能充分利用到CT扫描实验的技术特点。 　　文章以武汉岩土力学研究所科技人员提供的压缩荷载下泥岩CT实验数据为例进行分析，通过引入裂纹面积、分维计算、三维重建等多种手段，丰富了CT实验数据处理的方法，为多角度分析岩石应力应变过程提供了思路。 　　2 CT实验数据处理传统方法 　　本次试验采用中国矿业大学（北京）煤炭资源与安全开采国家重点实验室与美国BIR公司合作研发的ACTIS300-320/225CT/DR高分辨率工业CT实时成像系统。该系统可以快速扫描得到微米量级分辨率的像素大小为1024×1024的16位灰度切片图像。试验过程中针对试样自上而下间隔50μm依次扫描，共扫描1000层，其典型灰度切片图像如图1所示。 　　传统数据处理的方法多通过对不同应力状态下灰度切片图像进行灰度值的数理统计，得到其灰度频率直方图，再进一步进行对比分析[7]，不同应力水平下灰度频率直方图2所示。 　　3 CT实验数据处理新方法 　　3.1 裂纹面积处理法 　　图像的阈值分割及二值化可以定量的提取细观层面岩石的内部结构。采用单阈值分割法，即 　　（1） 　　即将小于某一阈值T的灰度级降为0，对应像素点为黑色，表征孔隙和裂纹；而大于该阈值的灰度升为65535，对应像素点为白色，表征岩样较完整，没有孔隙、裂纹或裂隙。二值化处理后的图像只有两个灰度级，因此称为二值图像。 　　计算发现：当阈值为26500时，基于初始扫描的CT二值图像求算的岩石孔隙率和压汞试验的结果相近，因此，将灰度阈值26500作为不同荷载水平下灰度图像二值化过程的基准阈值。裂纹面积-荷载水平关系示意图3所示。通过统计二值图像中黑色像素点的个数，进而求算初始孔隙与裂纹所占面积的百分比，并与初始扫描图层（未加载状态）作差，可追踪裂纹面积的变化规律如下： 　　各扫描层裂纹面积求算方法：第次扫描层裂纹面积 　　（2） 　　第层二值图中黑色像素点数；未加载时相应扫描层二值图中黑色像素点数；扫描层圆形统计区域内像素点总数；扫描层圆形统计区域的面积。 　　3.2 灰度图分维运算 　　谢和平[8]研究表明，岩石裂纹结构具有分形特征，且表现出一定的自相似性。几何上，分形维数D表征了岩石裂纹形态的复杂程度，D越大，裂纹图形越弯折，分叉越多，且呈现混乱无序状态；D越小，裂纹图形越光滑，越趋于规则有序排列。即分形维数D能定量表征裂纹的复杂程度与不规则性[9]。对二值化后的灰度图进行截图处理（圆形图幅区域内截取内接四边形区域），仅保留核心区域。截图处理后二值如图4所示。 　　同时采用计盒维数对截图处理后的二值图进行分维计算。原理如