CN107462744A - 一种基于电镜扫描的孔隙各向异性估算方法 - Google Patents
一种基于电镜扫描的孔隙各向异性估算方法 Download PDFInfo
- Publication number
- CN107462744A CN107462744A CN201710747525.5A CN201710747525A CN107462744A CN 107462744 A CN107462744 A CN 107462744A CN 201710747525 A CN201710747525 A CN 201710747525A CN 107462744 A CN107462744 A CN 107462744A
- Authority
- CN
- China
- Prior art keywords
- pieces
- debris
- coal sample
- stratification
- disk
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 44
- 238000011156 evaluation Methods 0.000 title claims abstract description 8
- 239000003245 coal Substances 0.000 claims abstract description 84
- 238000012360 testing method Methods 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000004458 analytical method Methods 0.000 claims abstract description 5
- 238000013517 stratification Methods 0.000 claims description 40
- 238000002474 experimental method Methods 0.000 claims description 9
- 239000004575 stone Substances 0.000 claims description 8
- 229910000831 Steel Inorganic materials 0.000 claims description 4
- 238000002389 environmental scanning electron microscopy Methods 0.000 claims description 4
- 238000004826 seaming Methods 0.000 claims description 4
- 239000010959 steel Substances 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 2
- 238000001035 drying Methods 0.000 claims description 2
- 230000005611 electricity Effects 0.000 claims description 2
- 238000000386 microscopy Methods 0.000 claims description 2
- 238000012986 modification Methods 0.000 claims description 2
- 230000004048 modification Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000003672 processing method Methods 0.000 claims description 2
- 239000013589 supplement Substances 0.000 claims description 2
- 238000002604 ultrasonography Methods 0.000 claims description 2
- 238000004140 cleaning Methods 0.000 claims 1
- 230000002035 prolonged effect Effects 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 7
- 239000011435 rock Substances 0.000 abstract description 3
- 238000005065 mining Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 abstract description 2
- 238000010998 test method Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 3
- 230000002265 prevention Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01Q—SCANNING-PROBE TECHNIQUES OR APPARATUS; APPLICATIONS OF SCANNING-PROBE TECHNIQUES, e.g. SCANNING PROBE MICROSCOPY [SPM]
- G01Q30/00—Auxiliary means serving to assist or improve the scanning probe techniques or apparatus, e.g. display or data processing devices
- G01Q30/04—Display or data processing devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
Landscapes
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Dispersion Chemistry (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
孔隙各向异性是研究煤力学性质和煤矿地质内容的重要内容。如何准确且快速地测定煤样的孔隙各向异性，就成为一个亟需解决的问题。而当前对煤体孔隙各向异性的测试方法领域基本空白。因此提出并设计了一种基于电镜扫描的孔隙各向异性估算方法。本专利通过巴西圆盘劈裂实验获得实验所需碎屑片，进而进行电镜扫描、MATLAB处理分析，最终得到相关数据。本专利是一个整体且系统的创新性方法，在相关领域中处于领先地位，对煤岩体动力学性质的研究具有重要意义，在未来工程中的应用具有广阔的前景和重要意义。
Description
所属技术领域
本发明属于煤岩体力学领域，提出了一种基于电镜扫描的孔隙各向异性估算方法，可以快速且较准确地估算出煤体孔隙各向异性的相关数据，填补相关领域的空白，对相关科研具有重要的辅助作用。
背景技术
孔隙各向异性是研究煤样力学性质的重要内容，也是研究煤矿地质内容的重要参考依据。通过研究煤体孔隙各向异性进而为抽采瓦斯或预防瓦斯突出提供重要依据；同时通过研究孔隙各向异性可进一步研究煤体微观结构的特点，为煤矿地质开采提供参考依据，进而达到预防开采过程中地质灾害的发生的目的。如何准确且快速地测定煤样的孔隙各向异性，就成为一个重要的研究内容。而当前对煤体孔隙各向异性的测试方法领域基本空白。
因此，有必要寻找一种可以测量现场煤体孔隙各向异性并且应用于现场的通用方法，目的是为现场测量煤体孔隙各向异性提供一整套的方法，这对煤体孔隙各向异性的测量及相关工程的应用与实践具有重要参考价值。
基于以上思想，提出并设计了一种基于电镜扫描的孔隙各向异性估算方法。该方法是一个整体且系统的创新性方法，在相关领域中处于领先地位，在未来实际工程中的应用具有广阔的前景和重要意义。
发明目的
本发明目的是提供一种可以在不同条件下快速、准确测量现场煤体孔隙各向异性的通用方法。并使得该方法所需器材要求较低，实验所需测量时间较短，实验步骤简单，最终能快速、准确得到煤体孔隙各向异性相关数据，以便对现场煤体地质构造微观结构进行研究，为抽采瓦斯或预防瓦斯突出及预防煤矿地质灾害的发生提供重要参考依据，该专利对煤岩体力学性质的研究具有重要意义。
附图说明
图1为本发明圆盘试件上各参考线参数示意图。图2为A4白纸围栏结构示意图。
实施方案
一种基于电镜扫描的孔隙各向异性估算方法包括以下步骤：
S101：选取所需估算孔隙各向异性区域的大小为50×50×50cm(±5cm)的煤样1块，严格参照国际岩石力学学会(ISRM)试验规程，将所选煤样加工成 7个直径50mm、长径比1:1、层理角度0°的巴西圆盘试件。
S102：找出各个圆盘的圆心，标记出各个圆盘的过圆心的层理参考线，将层理参考线绕圆心分别旋转0°、15°、30°、45°、60°、75°、90°得到新的直径，按照旋转后得到的直径在圆盘侧面画两条母线，母线与圆周交点确定为加载位置。
S103：对S102步骤所得圆盘试件进行巴西圆盘劈裂实验，实验完成后收集煤样抗拉破坏碎屑片，每个试件获得碎屑片不少于5块。
S104：对S103步骤所得的煤样抗拉破坏碎屑片参照《JY/T 010-1996分析型扫描电子显微镜方法通则》进行超声波及干燥处理，使用扫描电镜进行扫描，每块碎屑片获得断面照片不少于16块。
S105：采用MATLAB分析S104步骤所得的每个抗拉破坏碎屑片的每张断面照片中孔隙的几何参数，包括孔隙数目、孔隙面积、孔隙等效直径、孔隙偏心率。
S106：根据S105步骤所得的数据获得孔隙的各向异性曲线，包括孔隙数目的各向异性曲线、孔隙面积的各向异性曲线、孔隙等效直径比的各向异性曲线、孔隙偏心率比的各向异性曲线，四条曲线即是该种煤样的孔隙各向异性曲线。
S101步骤中所述选定所需估算孔隙各向异性区域的典型煤样方法具体为：选取需估算孔隙各向异性煤层的区域作为煤样采集点，通过弱爆破方法炸落一定量煤体，然后人工挑选尺寸为50×50×50cm(±5cm)的煤样1块。
S101步骤中所述层理角度为0°是使得圆盘端面垂直于层理面，其具体加工方法为：将挑选好的煤样放在工作台上，用地质罗盘测量其层理倾角和倾向ψ1，将煤样先旋转再旋转-ψ1(±5°)，采用锯石机对煤样此时的下表面进行切割，采用磨石机对切割的下表面进行磨平，再次使用地质罗盘测量其层理倾角与倾向ψ2，再将煤样先旋转再旋转-ψ2(±2°)，再次采用磨石机、锯石机对煤样下表面进行切割、磨平，重复该步骤，直至地质罗盘测得的煤样层理角度与ψ为0°(±1°)。
步骤S102所述找出各个圆盘的圆心的方法具体包括以下步骤：将加工好的巴西圆盘试件圆形端面C1朝上放置在水平工作台上，用量程为15cm的直尺量出圆盘端面C1上长度最长的弦，用线条粗细为1mm的红色记号笔沿着直尺在圆盘端面C1上画一条长度不超过3cm的线段a1，并且所画线段a1的两端点不与圆盘的圆周相交。再用直尺在圆盘端面C1上量出不与线段a1重合的长度最长的弦，用红色记号笔在圆盘端面C1上画一条长度不超过3cm、两端点不与圆盘圆周相交的线段a2，线段a2与线段a1的交于一点O，点O即为圆盘端面C1的圆心。按此方法依次在7个巴西圆盘试件的端面C1上作出圆心O。
步骤S102所述标记层理线的方法具体包括以下步骤：用红色记号笔在圆盘端面C1上沿着直尺画一条过圆盘圆心O、与圆盘圆周相交的线段b，线段b应与圆盘端面C1上的层理线重合(±2°)。作出的线段b即为圆盘的层理参考线。按此方法依次在7个巴西圆盘试件的端面C1上作出层理参考线b。
步骤S102所述确定加载位置的方法具体包括以下步骤：取一量程为 0°～180°的半圆量角器，将量角器的中心对准巴西圆盘试件端面C1上的圆心 O，量角器的0°刻度线与层理参考线b重合，用红色记号笔在量角器刻度为θ处标记一个点S，再用直尺将圆心O与点S连接并延长，作线段d，与圆盘端面C1的圆周相交于M、N两点。在7个巴西圆盘试件的端面C1上分别作出线段d以及M点与N点，角度θ分别取0°、15°、30°、45°、60°、75°、90°，并按照线段b与线段d之间的夹角分别为0°、15°、30°、45°、60°、75°、90°，依次将巴西圆盘试件命名为T1、T2、T3、T4、T5、T6、T7。用直尺和红色记号笔分别在巴西圆盘试件T1、T2、T3、T4、T5、T6、T7的侧面作出分别过点M与点 N的母线f1与f2，作为巴西圆盘劈裂试验中试验机钢制压条夹取的位置，如图 1。
步骤S103所述的获得碎屑片方法具体包括以下步骤：
步骤1：先用小刷子对RMT-150C型万能试验机的加载平台与劈裂夹具进行清理，目测没有无关碎屑后，将劈裂夹具内的上、下两个钢制压条分别紧贴于步骤S101中的巴西圆盘试件T1侧面的、步骤S102中确定出的加载位置f1与 f2，用V型卡具固定好试样，把装有巴西圆盘试件的劈裂夹具放入试验机的上、下承压板之间。安装好试件后，紧贴圆柱状加载平台的外缘粘贴至少两张 A4白纸，以防止加载过程中产生的碎屑片掉落到加载平台外部无法收集。白纸做成的围栏高度可通过裁剪进行调节，保证其高度不超过试样在试验机中最高点的高度，如图2。
步骤2：进行加载试验。试验完成后，摇动控制试验机上压头升降的旋钮使上压头升起，小心将加载平台的白纸围栏撕下，松开夹持圆盘试件的劈裂夹具，小心地将已破坏的试件及残渣全部倾倒在工作台铺好的A4白纸上，用小刷子把掉落在加载平台的碎屑清理收集到另一张A4白纸上。整个收集过程中，应避免将加载过程产生的碎屑片掉落到其他地方。
步骤3：从步骤3收集到的抗拉破坏碎屑片中挑选大小为1cm×1cm左右的薄片，薄片厚度不超过1cm，用镊子将选中的碎屑片轻轻夹取收集到一张A4 白纸上，要求每个试件获得5块碎屑片，分别编号为Qi1，Qi2，…，Qij，…， Qi5,，其中i为碎屑片对应巴西圆盘试件的编号数，取1，2，3，4，5，6，7；j 为碎屑片的编号数，取1，2，3，4，5。在盛放碎屑片的白纸上用黑色圆珠笔标明每个碎屑片的对应编号，从而获得层理线偏转角度分别为0°、15°、30°、45°、60°、75°、90°煤样的抗拉破坏碎屑片。
步骤S104所述的获得16张电镜扫描图像的方法具体包括以下步骤：采用 Quanta250型扫描电镜对煤样抗拉破坏碎屑片进行观测，采用的放大倍数为 500倍，此倍率下的分辨率为1024×943像素(1个像素的距离为0.291μm)。在碎屑片Qij的断面左上角距上边缘2mm、距左边缘2mm处获取第一张图片，之后横向每隔6872个像素点扫描一次，纵向每隔6872个像素点扫描一次，扫描位置在碎屑片断面上形成一个4×4、行间距为6872个像素点、列间距为6872 个像素点的矩形阵列。即每个煤样抗拉破坏碎屑片至少采集16张图片，分别编号为Pij1，Pij2，…，Pijk，…，Pij16，其中k为图像的编号数，取1～16。
步骤S105所述的MATLAB分析方法具体包括以下步骤：
步骤1：启动MATLAB软件，将步骤S103得到的煤样抗拉破坏碎屑片断面图像Pijk放入工作目录进行操作，输入减除背景灰度的指令，消除碎屑片断面图像Pijk中的不均匀亮度；运用MATLAB软件中的直方图均衡功能，对图像 Pijk进行增强；采用二值化阈值分割法对图像Pijk进行分割。即完成对煤样抗拉破坏碎屑片断面图像Pijk的处理工作。
步骤2：调用MATLAB软件中的bwlabel函数标记步骤1得到的二值图像 Pijk中的连通组，所谓连通组即孔隙，即可返回得到图像Pijk中的孔隙数目nijk。调用MATLAB软件中用于度量图像区域属性的函数regionprops，通过计算二值图像Pijk中的属性‘Area’，得到图像各个区域中像素的个数，各个区域的像素个数乘以1像素表示的面积即可求得各孔隙的面积，进而计算得到孔隙平均面积Sijk与孔隙平均等效直径Dijk；通过计算二值图像Pijk中的属性‘Eccentricity’，得到与图像中孔隙具有相同标准二阶中心距的椭圆的偏心率，即可进一步算出图像中孔隙平均偏心率eijk。将所得数据填入表1中。
表1抗拉破坏碎屑片断面图像Pijk中孔隙几何参数表
步骤S106所述参数计算方法具体为：
步骤1：对于所有同一层理线偏转角度，即i值相等的煤样抗拉破坏碎屑片的电镜扫描图像Pijk，对其中孔隙的孔隙数目nijk、孔隙平均面积Sijk、孔隙平均等效直径Dijk、孔隙平均偏心率eijk四种几何参数各自求平均值，作为该层理线偏转角度的煤样Ti中孔隙的对应参数，即得到煤样T1、T2、T3、T4、T5、T6、 T7中孔隙的孔隙数目ni、孔隙面积Si、孔隙等效直径Di、孔隙偏心率ei。具体为，煤样Ti中孔隙数目孔隙面积孔隙等效直径孔隙偏心率将数据填入表2中。
表2不同层理角度的煤样中孔隙几何参数表
步骤2：将表2中第2行层理线偏转角度的数据0°、15°、30°、45°、60°、 75°、90°作为横轴数据，第3行孔隙数目的数据n1、n2、n3、n4、n5、n6、n7作为纵轴数据，绘制得到的曲线即为煤样中孔隙数目的各向异性曲线。
步骤3：将表2中第2行层理线偏转角度的数据0°、15°、30°、45°、60°、 75°、90°作为横轴数据，第4行孔隙面积的数据S1、S2、S3、S4、S5、S6、S7作为纵轴数据，绘制得到的曲线即为煤样中孔隙面积的各向异性曲线。
步骤4：将表2中第2行层理线偏转角度的数据0°、15°、30°、45°、60°、 75°、90°作为横轴数据，第5行孔隙等效直径的数据D1、D2、D3、D4、D5、 D6、D7作为纵轴数据，绘制得到的曲线即为煤样中孔隙等效直径的各向异性曲线。
步骤5：将表2中第2行层理线偏转角度的数据0°、15°、30°、45°、60°、 75°、90°作为横轴数据，第6行孔隙偏心率的数据e1、e2、e3、e4、e5、e6、e7作为纵轴数据，绘制得到的曲线即为煤样中孔隙偏心率的各向异性曲线。
本发明是为了寻找一种可以快速、准确测量现场煤体孔隙各向异性并且应用于现场的通用方法。而当前对煤体孔隙各向异性的测试方法领域基本空白，因此提出一种基于电镜扫描的孔隙各向异性估算方法。本专利具有以下优点：一是通过巴西圆盘劈裂实验获得实验所需碎屑片，且不对煤体孔隙各向异性造成影响，大大减少了误差；二是通过电镜扫描、MATLAB处理分析获得相关数据，操作步骤快速、便捷，实验所需时间较少，为科研工作减少一定负担；三是实验数据的获取路径较短，减少了大量造成系统误差的因素。本专利对煤岩体动力学性质的研究具有重要意义。
本文中所描述的是对一种基于电镜扫描的孔隙各向异性估算的新方法的说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (7)
1.一种基于电镜扫描的孔隙各向异性估算方法，其特征在于包含以下步骤：
步骤S101：选取所需估算孔隙各向异性区域的大小为50×50×50cm(±5cm)的煤样1块，严格参照国际岩石力学学会(ISRM)试验规程，将所选煤样加工成7个直径50mm、长径比1:1、层理角度0°的巴西圆盘试件。
步骤S102：找出各个圆盘的圆心，标记出各个圆盘的过圆心的层理参考线，将层理参考线绕圆心分别旋转0°、15°、30°、45°、60°、75°、90°得到新的直径，按照旋转后得到的直径在圆盘侧面画两条母线，母线与圆周交点确定为加载位置。
步骤S103：对S102步骤所得圆盘试件进行巴西圆盘劈裂实验，实验完成后收集煤样抗拉破坏碎屑片，每个试件获得碎屑片不少于5块。
步骤S104：对S103步骤所得的煤样抗拉破坏碎屑片参照《JY/T 010-1996分析型扫描电子显微镜方法通则》进行超声波及干燥处理，使用扫描电镜进行扫描，每块碎屑片获得断面照片不少于16块。
步骤S105：采用MATLAB分析S104步骤所得的每个抗拉破坏碎屑片的每张断面照片中孔隙的几何参数，包括孔隙数目、孔隙面积、孔隙等效直径、孔隙偏心率。
步骤S106：根据S105步骤所得的数据获得孔隙的各向异性曲线，包括孔隙数目的各向异性曲线、孔隙面积的各向异性曲线、孔隙等效直径比的各向异性曲线、孔隙偏心率比的各向异性曲线，四条曲线即是该种煤样的孔隙各向异性曲线。
2.根据权利要求1所述的步骤S101，其特征在于包含以下方法：
S101步骤中所述选定所需估算孔隙各向异性区域的典型煤样方法具体为：选取需估算孔隙各向异性煤层的区域作为煤样采集点，通过弱爆破方法炸落一定量煤体，然后人工挑选尺寸为50×50×50cm(±5cm)的煤样1块。
S101步骤中所述层理角度为0°是使得圆盘端面垂直于层理面，其具体加工方法为：将挑选好的煤样放在工作台上，用地质罗盘测量其层理倾角和倾向ψ1，将煤样先旋转 再旋转-ψ1(±5°)，采用锯石机对煤样此时的下表面进行切割，采用磨石机对切割的下表面进行磨平，再次使用地质罗盘测量其层理倾角与倾向ψ2，再将煤样先旋转 再旋转-ψ2(±2°)，再次采用磨石机、锯石机对煤样下表面进行切割、磨平，重复该步骤，直至地质罗盘测得的煤样层理角度与ψ为0°(±1°)。
3.根据权利要求1所述的步骤S102，其特征在于包含以下步骤：
步骤S102所述找出各个圆盘的圆心的方法具体包括以下步骤：将加工好的巴西圆盘试件圆形端面C1朝上放置在水平工作台上，用量程为15cm的直尺量出圆盘端面C1上长度最长的弦，用线条粗细为1mm的红色记号笔沿着直尺在圆盘端面C1上画一条长度不超过3cm的线段a1，并且所画线段a1的两端点不与圆盘的圆周相交。再用直尺在圆盘端面C1上量出不与线段a1重合的长度最长的弦，用红色记号笔在圆盘端面C1上画一条长度不超过3cm、两端点不与圆盘圆周相交的线段a2，线段a2与线段a1的交于一点O，点O即为圆盘端面C1的圆心。按此方法依次在7个巴西圆盘试件的端面C1上作出圆心O。
步骤S102所述标记层理线的方法具体包括以下步骤：用红色记号笔在圆盘端面C1上沿着直尺画一条过圆盘圆心O、与圆盘圆周相交的线段b，线段b应与圆盘端面C1上的层理线重合(±2°)。作出的线段b即为圆盘的层理参考线。按此方法依次在7个巴西圆盘试件的端面C1上作出层理参考线b。
步骤S102所述确定加载位置的方法具体包括以下步骤：取一量程为0°～180°的半圆量角器，将量角器的中心对准巴西圆盘试件端面C1上的圆心O，量角器的0°刻度线与层理参考线b重合，用红色记号笔在量角器刻度为θ处标记一个点S，再用直尺将圆心O与点S连接并延长，作线段d，与圆盘端面C1的圆周相交于M、N两点。在7个巴西圆盘试件的端面C1上分别作出线段d以及M点与N点，角度θ分别取0°、15°、30°、45°、60°、75°、90°，并按照线段b与线段d之间的夹角分别为0°、15°、30°、45°、60°、75°、90°依次将巴西圆盘试件命名为T1、T2、T3、T4、T5、T6、T7。用直尺和红色记号笔分别在巴西圆盘试件T1、T2、T3、T4、T5、T6、T7的侧面作出分别过点M与点N的母线f1与f2，作为巴西圆盘劈裂试验中试验机钢制压条夹取的位置，如图1。
4.根据权利要求1所述的步骤S103，其特征在于包含以下步骤：
步骤1：先用小刷子对RMT-150C型万能试验机的加载平台与劈裂夹具进行清理，目测没有无关碎屑后，将劈裂夹具内的上、下两个钢制压条分别紧贴于步骤S101中的巴西圆盘试件T1侧面的、步骤S102中确定出的加载位置f1与f2，用V型卡具固定好试样，把装有巴西圆盘试件的劈裂夹具放入试验机的上、下承压板之间。安装好试件后，紧贴圆柱状加载平台的外缘粘贴至少两张A4白纸，以防止加载过程中产生的碎屑片掉落到加载平台外部无法收集。白纸做成的围栏高度可通过裁剪进行调节，保证其高度不超过试样在试验机中最高点的高度，如图2。
步骤2：进行加载试验。试验完成后，摇动控制试验机上压头升降的旋钮使上压头升起，小心将加载平台的白纸围栏撕下，松开夹持圆盘试件的劈裂夹具，小心地将已破坏的试件及残渣全部倾倒在工作台铺好的A4白纸上，用小刷子把掉落在加载平台的碎屑清理收集到另一张A4白纸上。整个收集过程中，应避免将加载过程产生的碎屑片掉落到其他地方。
步骤3：从步骤3收集到的抗拉破坏碎屑片中挑选大小为1cm×1cm左右的薄片，薄片厚度不超过1cm，用镊子将选中的碎屑片轻轻夹取收集到一张A4白纸上，要求每个试件获得5块碎屑片，分别编号为Qi1，Qi2，…，Qij，…，Qi5,，其中i为碎屑片对应巴西圆盘试件的编号数，取1，2，3，4，5，6，7；j为碎屑片的编号数，取1，2，3，4，5。在盛放碎屑片的白纸上用黑色圆珠笔标明每个碎屑片的对应编号，从而获得层理线偏转角度分别为0°、15°、30°、45°、60°、75°、90°煤样的抗拉破坏碎屑片。
5.根据权利要求1所述的步骤S104，其特征在于包含以下步骤：
步骤S104所述的获得16张电镜扫描图像的方法具体包括以下步骤：采用Quanta 250型扫描电镜对煤样抗拉破坏碎屑片进行观测，采用的放大倍数为500倍，此倍率下的分辨率为1024×943像素(1个像素的距离为0.291μm)。在碎屑片Qij的断面左上角距上边缘2mm、距左边缘2mm处获取第一张图片，之后横向每隔6872个像素点扫描一次，纵向每隔6872个像素点扫描一次，扫描位置在碎屑片断面上形成一个4×4、行间距为6872个像素点、列间距为6872个像素点的矩形阵列。即每个煤样抗拉破坏碎屑片至少采集16张图片，分别编号为Pij1，Pij2，…，Pijk，…，Pij16，其中k为图像的编号数，取1～16。
6.根据权利要求1所述的步骤S105，其特征在于包含以下步骤：
步骤1：启动MATLAB软件，将步骤S103得到的煤样抗拉破坏碎屑片断面图像Pijk放入工作目录进行操作，输入减除背景灰度的指令，消除碎屑片断面图像Pijk中的不均匀亮度；运用MATLAB软件中的直方图均衡功能，对图像Pijk进行增强；采用二值化阈值分割法对图像Pijk进行分割。即完成对煤样抗拉破坏碎屑片断面图像Pijk的处理工作。
步骤2：调用MATLAB软件中的bwlabel函数标记步骤1得到的二值图像Pijk中的连通组，所谓连通组即孔隙，即可返回得到图像Pijk中的孔隙数目nijk。调用MATLAB软件中用于度量图像区域属性的函数regionprops，通过计算二值图像Pijk中的属性‘Area’，得到图像各个区域中像素的个数，各个区域的像素个数乘以1像素表示的面积即可求得各孔隙的面积，进而计算得到孔隙平均面积Sijk与孔隙平均等效直径Dijk；通过计算二值图像Pijk中的属性‘Eccentricity’，得到与图像中孔隙具有相同标准二阶中心距的椭圆的偏心率，即可进一步算出图像中孔隙平均偏心率eijk。将所得数据填入表1中。
表1抗拉破坏碎屑片断面图像Pijk中孔隙几何参数表
7.根据权利要求1所述的步骤S106，其特征在于包含以下步骤：
步骤1：对于所有同一层理线偏转角度，即i值相等的煤样抗拉破坏碎屑片的电镜扫描图像Pijk，对其中孔隙的孔隙数目nijk、孔隙平均面积Sijk、孔隙平均等效直径Dijk、孔隙平均偏心率eijk四种几何参数各自求平均值，作为该层理线偏转角度的煤样Ti中孔隙的对应参数，即得到煤样T1、T2、T3、T4、T5、T6、T7中孔隙的孔隙数目ni、孔隙面积Si、孔隙等效直径Di、孔隙偏心率ei。具体为，煤样Ti中孔隙数目孔隙面积 孔隙等效直径孔隙偏心率将数据填入表2中。
表2不同层理角度的煤样中孔隙几何参数表
步骤2：将表2中第2行层理线偏转角度的数据0°、15°、30°、45°、60°、75°、90°作为横轴数据，第3行孔隙数目的数据n1、n2、n3、n4、n5、n6、n7作为纵轴数据，绘制得到的曲线即为煤样中孔隙数目的各向异性曲线。
步骤3：将表2中第2行层理线偏转角度的数据0°、15°、30°、45°、60°、75°、90°作为横轴数据，第4行孔隙面积的数据S1、S2、S3、S4、S5、S6、S7作为纵轴数据，绘制得到的曲线即为煤样中孔隙面积的各向异性曲线。
步骤4：将表2中第2行层理线偏转角度的数据0°、15°、30°、45°、60°、75°、90°作为横轴数据，第5行孔隙等效直径的数据D1、D2、D3、D4、D5、D6、D7作为纵轴数据，绘制得到的曲线即为煤样中孔隙等效直径的各向异性曲线。
步骤5：将表2中第2行层理线偏转角度的数据0°、15°、30°、45°、60°、75°、90°作为横轴数据，第6行孔隙偏心率的数据e1、e2、e3、e4、e5、e6、e7作为纵轴数据，绘制得到的曲线即为煤样中孔隙偏心率的各向异性曲线。
本文中所描述的仅仅是对一种基于电镜扫描的孔隙各向异性估算方法的举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710747525.5A CN107462744A (zh) | 2017-08-28 | 2017-08-28 | 一种基于电镜扫描的孔隙各向异性估算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710747525.5A CN107462744A (zh) | 2017-08-28 | 2017-08-28 | 一种基于电镜扫描的孔隙各向异性估算方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107462744A true CN107462744A (zh) | 2017-12-12 |
Family
ID=60550555
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710747525.5A Pending CN107462744A (zh) | 2017-08-28 | 2017-08-28 | 一种基于电镜扫描的孔隙各向异性估算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107462744A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108645876A (zh) * | 2018-03-07 | 2018-10-12 | 中山大学 | 一种内外连续跨尺度量测方法及装置 |
| CN112098298A (zh) * | 2020-09-15 | 2020-12-18 | 安徽明天氢能科技股份有限公司 | 一种测试燃料电池碳纸各向异性渗透率的装置及方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102183410A (zh) * | 2011-01-27 | 2011-09-14 | 中国科学院武汉岩土力学研究所 | 一种测量拉伸条件下岩石弹性参数的巴西劈裂方法 |
| CN102841220A (zh) * | 2012-07-16 | 2012-12-26 | 天津城市建设学院 | 一种基于孔隙率的黏土扫描电镜照片图像分割方法 |
| CN104698020A (zh) * | 2013-12-06 | 2015-06-10 | 中国石油天然气股份有限公司 | 一种松散岩心微观孔隙结构特征参数的采集处理方法 |
| CN105334149A (zh) * | 2015-11-24 | 2016-02-17 | 中国石油大学(北京) | 一种致密储层微观孔隙结构评价与储层分类方法 |
| CN105973706A (zh) * | 2016-06-07 | 2016-09-28 | 中国矿业大学(北京) | 一种基于工业ct的煤岩体多尺度力学特性分析方法 |
-
2017
- 2017-08-28 CN CN201710747525.5A patent/CN107462744A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102183410A (zh) * | 2011-01-27 | 2011-09-14 | 中国科学院武汉岩土力学研究所 | 一种测量拉伸条件下岩石弹性参数的巴西劈裂方法 |
| CN102841220A (zh) * | 2012-07-16 | 2012-12-26 | 天津城市建设学院 | 一种基于孔隙率的黏土扫描电镜照片图像分割方法 |
| CN104698020A (zh) * | 2013-12-06 | 2015-06-10 | 中国石油天然气股份有限公司 | 一种松散岩心微观孔隙结构特征参数的采集处理方法 |
| CN105334149A (zh) * | 2015-11-24 | 2016-02-17 | 中国石油大学(北京) | 一种致密储层微观孔隙结构评价与储层分类方法 |
| CN105973706A (zh) * | 2016-06-07 | 2016-09-28 | 中国矿业大学(北京) | 一种基于工业ct的煤岩体多尺度力学特性分析方法 |
Non-Patent Citations (2)
| Title |
|---|
| BAO-HUA YANG等: "3D characterization and analysis of pore structure of packed ore particle beds based on computed tomography images", 《TRANS. NONFERROUS MET. SOC. CHINA》 * |
| 侯 鹏等: "黑色页岩巴西劈裂破坏的层理效应研究及能量分析", 《岩土工程学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108645876A (zh) * | 2018-03-07 | 2018-10-12 | 中山大学 | 一种内外连续跨尺度量测方法及装置 |
| CN112098298A (zh) * | 2020-09-15 | 2020-12-18 | 安徽明天氢能科技股份有限公司 | 一种测试燃料电池碳纸各向异性渗透率的装置及方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Mora et al. | Sphericity, shape factor, and convexity measurement of coarse aggregate for concrete using digital image processing | |
| CN105352873B (zh) | 页岩孔隙结构的表征方法 | |
| Liu et al. | Multi-scale fractal analysis of pores in shale rocks | |
| US8938045B2 (en) | Method of determining reservoir properties and quality with multiple energy X-ray imaging | |
| EP1813934A1 (en) | Concrete structure crack inspection device and crack inspection method | |
| Muszynski et al. | Particle shape estimates of uniform sands: visual and automated methods comparison | |
| CN1693852A (zh) | 一种仪表指针自动检测识别方法及自动读数方法 | |
| CN110415167A (zh) | 一种基于数字图像技术的粗糙面裂隙生成方法及试验系统 | |
| CN206095779U (zh) | 拉伸试样断后参数测量用夹具 | |
| CN103593668B (zh) | 一种金属板材冲压接头力学性能试验裂纹自动识别方法 | |
| CN102607977A (zh) | 基于数字图像处理的磨损原位测量装置及方法 | |
| Herwegh | A new technique to automatically quantify microstructures of fine grained carbonate mylonites: two-step etching combined with SEM imaging and image analysis | |
| CN210465317U (zh) | 一种适用于ebsd测试的样品座 | |
| CN107462744A (zh) | 一种基于电镜扫描的孔隙各向异性估算方法 | |
| CN101571381A (zh) | 提取材料断面轮廓曲线并实现断面轮廓分形表征的方法 | |
| Al Ibrahim et al. | An automated petrographic image analysis system: Capillary pressure curves using confocal microscopy | |
| CN104268830A (zh) | 基于数字图像确定不均匀岩土材料渗透系数的方法 | |
| Rivera et al. | Automated detection and measurement of cracks in reinforced concrete components | |
| Yong et al. | A rapid field measurement method for the determination of Joint Roughness Coefficient of large rock joint surfaces | |
| CN111624146A (zh) | 一种定量表征致密储层全尺度孔径分布特征的方法 | |
| CN103630093A (zh) | 用于混凝土表面粗糙度检测的图像分析方法 | |
| CN108335310B (zh) | 一种便携式粒形粒度检测方法及系统 | |
| CN102183794A (zh) | 一种井中三分量磁测系统 | |
| CN112991287A (zh) | 一种基于全卷积神经网络的自动压痕测量方法 | |
| CN111028094B (zh) | 一种基于页岩扫描电镜图像提取的等效toc计算方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171212 |
|
| RJ01 | Rejection of invention patent application after publication |