一种基于工业CT的煤岩体多尺度力学特性分析方法技术

一种基于工业CT的煤岩体多尺度力学特性分析方法，是一种能够与CT扫描特征相匹配的煤岩多尺度力学特性分析方法，包括顺序切割、CT精确扫描、基质夹杂重塑、分别三轴实验获力学参数、顺序组合进行数值压缩的一整套方法实现煤岩体多尺度力学特性分析，突破了传统的基于概率统计的随机模拟方法的局限性。该方法首先在开采矿井选取块体切割成1m×1m×1m的立方体，然后再次切割成1000块0.1m×0.1m×0.1m的小煤块并顺序编号和CT扫描，并再次顺序重新组合获得大立方体网格模型；将切后煤体磨成煤粉并进行重塑，分别进行三轴实验获得基质与杂质的力学参数，带入网格模型计算不同体积煤块的相关物理力学参数，从而获得表征单元体积及对应参数。

【技术实现步骤摘要】
所属
本专利技术属于岩石力学领域，尤其是一种煤岩体的多尺度分析方法，是对一般的煤岩体多尺度分析方法的一种补充，尤其适用于多孔介质的多尺度分析方法。
技术介绍
进行煤岩体多尺度的分析，了解其内部结构的组成是基础。现有的煤岩体多尺度分析方法中，为了得到煤岩体内部结构组成，往往首先构建其内部结构面网络。即首先通过现场节理地质调查，针对节理的倾向、倾角、间距、迹长进行概率统计分析及偏差校正；然后，根据校正后的数据，采用随机模拟理论构建节理网络；最后，将节理网络数据与节理偏差校正数据进行对比，验证两者的统计相似性。可以看出，上述内部结构组成的获取主要是基于统计，其优点在于可以以局部代替整体，但是其缺点也在于此，即其获得的内部结构仅仅具有统计学意义，对于其内部真实情况往往难以验证，而由此得到的结果也难以判定正确与否。而煤岩体往往具有显著的非均质性，即使地理相近，也可能受到受到应力、构造等因素的影响，具有显著的变异特征，这时以统计作为基础的传统多尺度分析方法对这种变异分析往往无能为力，这也是现有岩体力学多尺度研究的困境。煤岩体中赋存的大小不一、分布不均的孔隙和夹杂，极大影响着煤体的物理性质，从而给分析的结果带来很大误差。因此，有必要寻找一种以能精确分析煤岩体内部结构特征为基础的多尺度分析方法,该方法要求能够再现煤体内部的真实结构，而以无损化定量化精细化为特点的CT扫描技术可以很好的解决这一问题。因此本专利技术以工业CT扫描煤样获得其细观特征为基础，构建一种与CT扫描特征相匹配的新型煤岩体多尺度力学特性的分析方法。基于以上思想，设计了该种基于工业CT的煤岩体多尺度力学特性的分析方法。
技术实现思路
该专利技术首先在开采矿井选取块体切割成1m×1m×1m的立方体，然后再次切割成1000块0.1m×0.1m×0.1m的小煤块并顺序编号和CT扫描，并再次顺序重新组合获得大立方体网格模型；将切后煤体磨成煤粉并进行重塑，分别进行三轴实验获得基质与杂质的力学参数，带入网格模型计算不同体积煤块的相关物理力学参数，从而获得表征单元体积及对应参数。为了实现上述目的，本专利技术采用如下技术方案：一种基于工业CT的煤岩体多尺度分析方法包括以下步骤：S101：在矿井煤炭开采工作面选择煤体块体，块体何种形状均可，但是要求该块体每条棱边的最小长度均大于1.5m。S102：将S101中所述煤体块体采用矿山锯(型号:KSG-3600)切割成1m×1m×1m的立方体。S103：将S102中所述立方体切割成1000块0.1m×0.1m×0.1m的小煤块，并以每个煤块切割前所在位置为依据，依次给每一个煤块进行编号。S104：将S102步骤中切后剩余的煤体磨成煤粉，分别找到其中的基质部分和夹杂部分。S105：对煤粉中的基质和杂质分别进行重塑，随后分别进行三轴压缩试验，获得煤块中基质与夹杂的力学参数。S106：利用加载式工业CT(型号：ACTIS300-320/225)对S103步骤中所述小煤块按照专利“一种基于CT扫描的受载煤岩损伤本构方程构建方法”所述方法依次进行扫描、网格重建，以获得每个煤块内部基质、孔隙以及夹杂的分布，并导入数值模拟软件中建立煤样网格模型，从而可以得到1000个小煤块的网格模型。S107：以S103步骤中每个煤块切割前所在位置及编号为依据，得到不同体积大小的正方体煤块网格模型。S108：对S106步骤中不同体积大小的正方体煤块网格模型中的基质与夹杂按照S104步骤中力学参数进行赋值，带入数值模拟软件中进行计算，并记录下不同体积煤块的相关物理力学参数。S109：煤体物理参数随煤块体积增大而减小，到达REV体积时物理参数不再改变。因此，找出物理参数随体积不再改变的点即为对应的REV点，而该点所对应的物理参数即为煤体表征单元体积的物理参数。其中，步骤S103所述的大煤块切割及编号方法具体包括以下步骤：步骤1、在1m×1m×1m的正方体大煤块的每一面上分别等距画出九条横向细线和九条纵向细线，即将每一面都划分为100个相同的小方格。步骤2、以正方体煤块的右下端作为原点，相交于原点的三条棱边分别分为X轴、Y轴、Z轴，先沿X轴对划分出的煤块进行排序，从X轴原点到煤块X轴终点依次排序1,2…10。再沿Y轴对划分出的煤块进行排序，并且从Y轴原点到煤块Y轴终点排序1,2…10。再沿Z轴对划分出的煤块进行排序，并且从Z轴原点到煤块Z轴终点排序1,2…10。步骤3、按照上述排序对1000块小煤块进行编号，分别编号为XiYjZk，其中i为该煤块在X轴对应的顺序，j为该煤块在Y轴对应的顺序，k为该煤块在Z轴对应的顺序。步骤4、沿着所画细线将每一个小煤块切割下来，即1000块0.1m×0.1m×0.1m的小煤块，其对应编号为XiYjZk。其中，步骤S104所述的煤粉中基质和夹杂部分的获取方法具体包括以下步骤：步骤1、取一块1dm3左右(近似即可，不需精确)的步骤S102中切下剩余的煤块，采用电子秤称量其质量m。再将其放入一个可以将煤块完全放下的圆柱形量杯中，加入大于1dm3的水，读出初始时水面位置对应的体积V1，将煤块缓缓放入，待水面平静后，再次读出此时水面位置对应的体积V2，可以得到煤块的体积V3为(V2-V1)。煤块密度则为ρ＝m/V3。步骤2、将S102中切下剩余的煤块的用锤式破碎机(型号:PCΦ400×300)破碎成粒度为25mm以下的煤块，再用高压型微粉磨(型号：YGM9531)将破碎后的煤块磨成粒度为90μm以下的煤粉。步骤3、对步骤2中得到的煤粉进行重液离心分离。比重液选用吉林冶金研究所开发的JY型无毒变温型重液，根据步骤1中得到的煤块密度ρ，在(ρ-0.5)g/cm3～(ρ+0.5)g/cm3之间等间距配置十组比重液，每组容量为1000mL。再向每组中加入100cm3混合均匀后的煤粉，观察每组中各级密度物的分离量，找出分离程度最大的一组。步骤4、再配制100L与步骤3中分离程度最大组所对应的密度相同的比重液，加入0.01cm3的煤粉，待其充分分离后，收集得到的各级密度物，包括煤粉基质和夹杂两部分(其中浮在比重液上部的是煤粉基质部分，沉在底部的是夹杂部分)，并分别贴上标签。即完成煤粉中基质和夹杂部分的获取。其中，步骤S105所述的煤粉中基质与夹杂部分分别进行重塑与力学参数的获取方法包括以下步骤:步骤1、将160cm3在S105步骤中获得的煤粉基质放入一个直径×高为50×100mm的圆柱形模具容器中，再加入20cm3的混凝土、20mL并搅拌均匀，使其装满整个容器。步骤2、对步骤1重复15次，完成15块圆柱形煤体基质的重塑。步骤3、放置20天进行干燥，20天后打开模具，获得成型后的的煤粉基质圆柱块。步骤4、对步骤3中的煤粉基质圆柱块进行5MPa、10MPa、20MPa下的围压实验，每个围压进行5块，获得15块煤粉基质圆柱块的强度、弹性模量，利用三个围压下的强度值构建摩尔库伦圆，获得煤粉基质的粘聚力、内摩擦角。步骤5、利用步骤4获得的15块煤粉基质圆柱块的弹性模量平均值作为煤粉基质的弹性模量。步骤6、对煤粉夹杂部分也按照步骤1、步骤2、步骤3、步骤4、步骤5进行重塑、三轴压缩，得到15块0.1m×0.1m×0.1m的煤粉夹杂圆柱块，同样进行5MPa、10本文档来自技高网...

【技术保护点】
一种基于工业CT的煤岩体多尺度分析方法包括以下步骤：S101：在矿井煤炭开采工作面选择煤体块体，块体何种形状均可，但是要求该块体每条棱边的最小长度均大于1.5m。S102：将S101中所述煤体块体采用矿山锯(型号:KSG‑3600)切割成1m×1m×1m的立方体。S103：将S102中所述立方体切割成1000块0.1m×0.1m×0.1m的小煤块，并以每个煤块切割前所在位置为依据，依次给每一个煤块进行编号。S104：将S102步骤中切后剩余的煤体磨成煤粉，分别找到其中的基质部分和夹杂部分。S105：对煤粉中的基质和杂质分别进行重塑，随后分别进行三轴压缩试验，获得煤块中基质与夹杂的力学参数。S106：利用加载式工业CT(型号：ACTIS300‑320/225)对S103步骤中所述小煤块按照专利“一种基于CT扫描的受载煤岩损伤本构方程构建方法”所述方法依次进行扫描、网格重建，以获得每个煤块内部基质、孔隙以及夹杂的分布，并导入数值模拟软件中建立煤样网格模型，从而可以得到1000个小煤块的网格模型。S107：以S103步骤中每个煤块切割前所在位置及编号为依据，得到不同体积大小的正方体煤块网格模型。S108：对S106步骤中不同体积大小的正方体煤块网格模型中的基质与夹杂按照S104步骤中力学参数进行赋值，带入数值模拟软件中进行计算，并记录下不同体积煤块的相关物理力学参数。S109：煤体物理参数随煤块体积增大而减小，到达REV体积时物理参数不再改变。因此，找出物理参数随体积不再改变的点即为对应的REV点，而该点所对应的物理参数即为煤体表征单元体积的物理参数。其中，步骤S103所述的大煤块切割及编号方法具体包括以下步骤：步骤1、在1m×1m×1m的正方体大煤块的每一面上分别等距画出九条横向细线和九条纵向细线，即将每一面都划分为100个相同的小方格。步骤2、以正方体煤块的右下端作为原点，相交于原点的三条棱边分别分为X轴、Y轴、Z轴，先沿X轴对划分出的煤块进行排序，从X轴原点到煤块X轴终点依次排序1,2…10。再沿Y轴对划分出的煤块进行排序，并且从Y轴原点到煤块Y轴终点排序1,2…10。再沿Z轴对划分出的煤块进行排序，并且从Z轴原点到煤块Z轴终点排序1,2…10。步骤3、按照上述排序对1000块小煤块进行编号，分别编号为XiYjZk，其中i为该煤块在X轴对应的顺序，j为该煤块在Y轴对应的顺序，k为该煤块在Z轴对应的顺序。步骤4、沿着所画细线将每一个小煤块切割下来，即1000块0.1m×0.1m×0.1m的小煤块，其对应编号为XiYjZk。其中，步骤S104所述的煤粉中基质和夹杂部分的获取方法具体包括以下步骤：步骤1、取一块1dm3左右(近似即可，不需精确)的步骤S102中切下剩余的煤块，采用电子秤称量其质量m。再将其放入一个可以将煤块完全放下的圆柱形量杯中，加入大于1dm3的水，读出初始时水面位置对应的体积V1，将煤块缓缓放入，待水面平静后，再次读出此时水面位置对应的体积V2，可以得到煤块的体积V3为(V2‑V1)。煤块密度则为ρ＝m/V3。步骤2、将S102中切下剩余的煤块的用锤式破碎机(型号:PCΦ400×300)破碎成粒度为25mm以下的煤块，再用高压型微粉磨(型号：YGM9531)将破碎后的煤块磨成粒度为90μm以下的煤粉。步骤3、对步骤2中得到的煤粉进行重液离心分离。比重液选用吉林冶金研究所开发的JY型无毒变温型重液，根据步骤1中得到的煤块密度ρ，在(ρ‑0.5)g/cm3～(ρ+0.5)g/cm3之间等间距配置十组比重液，每组容量为1000mL。再向每组中加入100cm3混合均匀后的煤粉，观察每组中各级密度物的分离量，找出分离程度最大的一组。步骤4、再配制100L与步骤3中分离程度最大组所对应的密度相同的比重液，加入0.01cm3的煤粉，待其充分分离后，收集得到的各级密度物，包括煤粉基质和夹杂两部分(其中浮在比重液上部的是煤粉基质部分，沉在底部的是夹杂部分)，并分别贴上标签。即完成煤粉中基质和夹杂部分的获取。其中，步骤S105所述的煤粉中基质与夹杂部分分别进行重塑与力学参数的获取方法包括以下步骤:步骤1、将160cm3在S105步骤中获得的煤粉基质放入一个直径×高为50×100mm的圆柱形模具容器中，再加入20cm3的混凝土、20mL并搅拌均匀，使其装满整个容器。步骤2、对步骤1重复15次，完成15块圆柱形煤体基质的重塑。步骤3、放置20天进行干燥，20天后打开模具，获得成型后的的煤粉基质圆柱块。步骤4、对步骤3中的煤粉基质圆柱块进行5MPa、10MPa、20MPa下的围压实验，每个围压进行5块，获得15块煤粉基质圆柱块的强度、弹性模量，利用三个围压下的强度值构建摩尔库伦圆，获得煤粉基质的粘聚力、内摩擦角。步骤5、利用步骤4获得的1...

【技术特征摘要】
1.一种基于工业CT的煤岩体多尺度分析方法包括以下步骤：S101：在矿井煤炭开采工作面选择煤体块体，块体何种形状均可，但是要求该块体每条棱边的最小长度均大于1.5m。S102：将S101中所述煤体块体采用矿山锯(型号:KSG-3600)切割成1m×1m×1m的立方体。S103：将S102中所述立方体切割成1000块0.1m×0.1m×0.1m的小煤块，并以每个煤块切割前所在位置为依据，依次给每一个煤块进行编号。S104：将S102步骤中切后剩余的煤体磨成煤粉，分别找到其中的基质部分和夹杂部分。S105：对煤粉中的基质和杂质分别进行重塑，随后分别进行三轴压缩试验，获得煤块中基质与夹杂的力学参数。S106：利用加载式工业CT(型号：ACTIS300-320/225)对S103步骤中所述小煤块按照专利“一种基于CT扫描的受载煤岩损伤本构方程构建方法”所述方法依次进行扫描、网格重建，以获得每个煤块内部基质、孔隙以及夹杂的分布，并导入数值模拟软件中建立煤样网格模型，从而可以得到1000个小煤块的网格模型。S107：以S103步骤中每个煤块切割前所在位置及编号为依据，得到不同体积大小的正方体煤块网格模型。S108：对S106步骤中不同体积大小的正方体煤块网格模型中的基质与夹杂按照S104步骤中力学参数进行赋值，带入数值模拟软件中进行计算，并记录下不同体积煤块的相关物理力学参数。S109：煤体物理参数随煤块体积增大而减小，到达REV体积时物理参数不再改变。因此，找出物理参数随体积不再改变的点即为对应的REV点，而该点所对应的物理参数即为煤体表征单元体积的物理参数。其中，步骤S103所述的大煤块切割及编号方法具体包括以下步骤：步骤1、在1m×1m×1m的正方体大煤块的每一面上分别等距画出九条横向细线和九条纵向细线，即将每一面都划分为100个相同的小方格。步骤2、以正方体煤块的右下端作为原点，相交于原点的三条棱边分别分为X轴、Y轴、Z轴，先沿X轴对划分出的煤块进行排序，从X轴原点到煤块X轴终点依次排序1,2…10。再沿Y轴对划分出的煤块进行排序，并且从Y轴原点到煤块Y轴终点排序1,2…10。再沿Z轴对划分出的煤块进行排序，并且从Z轴原点到煤块Z轴终点排序1,2…10。步骤3、按照上述排序对1000块小煤块进行编号，分别编号为XiYjZk，其中i为该煤块在X轴对应的顺序，j为该煤块在Y轴对应的顺序，k为该煤块在Z轴对应的顺序。步骤4、沿着所画细线将每一个小煤块切割下来，即1000块0.1m×0.1m×0.1m的小煤块，其对应编号为XiYjZk。其中，步骤S104所述的煤粉中基质和夹杂部分的获取方法具体包括以下步骤：步骤1、取一块1dm3左右(近似即可，不需精确)的步骤S102中切下剩余的煤块，采用电子秤称量其质量m。再将其放入一个可以将煤块完全放下的圆柱形量杯中，加入大于1dm3的水，读出初始时水面位置对应的体积V1，将煤块缓缓放入，待水面平静后，再次读出此时水面位置对应的体积V2，可以得到煤块的体积V3为(V2-V1)。煤块密度则为ρ＝m/V3。步骤2、将S102中切下剩余的煤块的用锤式破碎机(型号:PCΦ400×300)破碎成粒度为25mm以下的煤块，再用高压型微粉磨(型号：YGM9531)将破碎后的煤块磨成粒度为90μm以下的煤粉。步骤3、对步骤2中得到的煤粉进行重液离心分离。比重液选用吉林冶金研究所开发的JY型无毒变温型重液，根据步骤1中得到的煤块密度ρ，在(ρ-0.5)g/cm3～(ρ+0.5)g/cm3之间等间距配置十组比重液，每组容量为1000mL。再向每组中加入100cm3混合均匀后的煤粉，观察每组中各级密度物的分离量，找...

【专利技术属性】
技术研发人员：郝宪杰，袁亮，贾仕刚，卢志国，赵毅鑫，张通，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：北京;11