WO2020042512A1

WO2020042512A1 - 基于ct扫描的预测应力加载条件下煤体渗透率的方法

Info

Publication number: WO2020042512A1
Application number: PCT/CN2018/125959
Authority: WO
Inventors: 王刚; 江成浩; 程卫民; 孙路路; 刘义鑫; 范酒源
Original assignee: Shandong University of Science and Technology
Current assignee: Shandong University of Science and Technology
Priority date: 2018-08-31
Filing date: 2018-12-30
Publication date: 2020-03-05
Anticipated expiration: 2021-02-28
Also published as: CN109211666A; CN109211666B

Abstract

一种基于CT扫描的预测应力加载条件下煤体渗透率的方法，解决了变形煤体渗透率计算的技术问题，该方法步骤包括：A.制作煤样试件，进行单轴压缩实验，同时对煤样试件进行CT扫描；B.采用数字地形模型法进行阀值分割，得到煤样试件的阀值，将CT扫描图像导入建模软件建立三维数值模型；C.将三维数值模型导入有限元分析软件中并设置渗流条件；D.对三维数值模型施加多个压力梯度值▽P，设置初始流动速度v 0，调整▽P和v 0参数，得到多个渗流模型；E.将渗流模型导入有限元软件中模拟计算，得到计算结果；F.将计算结果导入数据处理软件中，沿渗流方向等距提取渗流速度，得到渗流速度和压力梯度的关系，计算渗透率。

Description

基于 CT扫描的预测应力加载条件下煤体渗透率的方法技术领域

[0001] 本发明涉及煤层注水渗流技术领域，尤其是一种基于 CT扫描的预测并计算应力加载条件下煤体渗透率的方法。

背景技术

[0002] 煤矿生产过程中，需要通过煤层注水实现对工作面粉尘污染、煤与瓦斯突出、冲击地压、自然发火、煤体软化等矿井实际生产问题进行综合治理，而注水后的渗流是自然界最常见的物理现象。在一定压差下，岩石允许流体通过的能力叫渗透率，渗透率是表征土或岩石本身传导液体能力的参数，其大小与孔隙度、液体渗透方向上空隙的几何形状、颗粒大小以及排列方向等因素有关，渗透率用来表示渗透性的大小。

[0003] 传统的渗透率研究方法包括：（1）通过岩石试验系统和真三轴试验系统等对稳定渗流状态下的煤体应力 -渗流耦合作用进行研究, 在轴压和侧压稳定的情况下，通过测取流量、流动长度、压力梯度等参数，并结合达西渗流公式也可以实现对煤体渗透率的计算。但是由于煤体结构的非均质不透明等特性，所以该方法不能对煤体内部微细观孔裂隙结构的渗流过程进行解释和描述，而复杂的外界环境会影响流量等参数的确定，从而降低渗透率计算结果的准确性。（2）对于非稳态渗流过程中的渗透率变化规律, 常采用脉冲衰减法进行研究, 其基本工作原理为通过测试岩样一维非稳态渗流过程中孔隙压力随时间的衰减数据，在此基础上，结合数学模型对煤样脉冲渗透率进行测取。但是该方法在测试高渗透煤样时误差很大，这限制了此方法的测试范围。

[0004] 随着 CT扫描技术在岩石领域的大量应用，其具有的高精度、无损伤等特点，可以帮助人们捕获各种储层岩石的孔裂隙空间几何形状，实现对于孔隙率的测定以及孔裂隙结构特征、矿物形态的分析；同时，借助载荷加载装置还可以实现对于孔裂隙结构演化规律的研究。在气体渗流方面， CT扫描技术主要侧重于气体吸附和解吸作用研究；而在注水渗流方面，利用该技术则可以实现渗流通道的构建以及揭示注水渗流过程中微细观结构的损伤破坏机理。但是受到监测手段的制约，目前的 CT扫描实验仍然无法实现对渗流通道中流速、压力梯度变化情况进行详细的测量统计，同时复杂的实验测试也增大了研究成本。因此，为进一步的提高煤体渗透率计算的准确性，实现对煤体内部渗流分布情况的描述，有效预测不同应力条件下渗透率的变化情况，需要对现有的技术有待于更进一步的改善和发展。

发明概述

技术问题

问题的解决方案

技术解决方案

[0005] 为解决实现直观描述变形煤样内部渗流分布情况，并有效预测不同应力条件下煤体渗透率变化的技术问题，本发明提供了一种基于 CT扫描的预测应力加载条件下煤体渗透率的方法，具体技术方案如下：

[0006] 基于 CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，包括以下步骤：

[0007] A.制作煤样试件，进行单轴压缩实验，同时对煤样试件进行 CT扫描；

[0008] B.采用数字地形模型法进行阀值分割，并得到煤样试件的阀值，将 CT扫描图像导入 Simpleware软件中建立三维数值模型；

[0009] C.将三维数值模型导入 HyperMesh软件中并设置渗流条件；

[0010] D.对三维数值模型施加多个压力梯度值 VP，设置初始流动速度 vO, 调整 VP和 vO参数分别模拟，得到渗流模型；

[0011] E.将渗流模型导入 Ls-dyna软件中模拟计算，得到计算结果；

[0012] F.将计算结果导入 HyperView软件中，沿渗流方向等距提取渗流速度，得到渗流速度和压力梯度的关系，计算渗透率。

[0013] 优选的是，煤样试件为圆柱体，所述单轴压缩实验使用原位拉伸、压缩及温度控制实验装置；所述单轴压缩试验中分别进行加载速度不同的多组煤样试件的单轴压缩，同一煤样试件的加载速度恒定。

[0014] 进一步优选的是，步骤 A中，在 CT扫描前首先根据 X射线稳定性、煤样试件尺寸、煤样试件 X射线衰减分数和曝光时间来确定扫描电压、扫描功率和视场大小 ; 所述 CT扫描是以恒定的速度旋转扫描，检测器捕捉由 X射线源发出穿过煤样试件的 X射线，以电信号的形式储存 CT扫描图像。

[0015] 还优选的是，步骤 B中确定煤样试件的阀值包括：将 CT扫描图像转化为数字地形模型，表示孔裂隙体积和煤体总体积，孔裂隙体积与煤体总体积的比值作为孔隙率；建立孔隙率和图像灰度值的函数关系，计算孔隙率函数所有极小值中的最大值，作为 CT扫描图像的阀值。

[0016] 还进一步优选的是，步骤 B中在 Simpleware软件建立三维数值模型后进行网格划分。

[0017] 还优选的是，步骤 C中的渗流条件包括水场模拟、空气场模拟和煤体模拟；所述水场模拟选用 MAT-NULL本构模型，设置流体沿一个方向流动，在流动方向上的两个侧面施加垂直于流动方向的两个方向的约束，在另外四个侧面施加固定约束；通过布置线性多项式状态方程为水头处流体的运动提供压力。

[0018] 还进一步优选的是，空气场模拟中选用 MAT- VACUUM空隙材料模型，并选择 ALE算法作为运算方法。

[0019] 还进一步优选的是，煤体模拟选用拉格朗日算法作为运算方法。

[0020] 还优选的是，压力梯度值 VP的选择范围为 0~300Pa/mm; 所述初始流动速度 v0 的选择范围为 0.01~0.09mm/s。

[0021] 还进一步优选的是，骤 D中渗流模型施加的压力梯度值 VP分别为 OPa/mm、 10P a/mm、 24Pa/mm和 213Pa/mm; 初始流动速度 vO设置为 0.02mm/s。

发明的有益效果

有益效果

[0022] 本发明的有益效果包括：

[0023] ( 1) 本方法通过室内试验和 CT扫描获得煤体的结构特征，并采用数字地形模型法进行阀值分割，得到煤样阀值并建立三维数值模型，从而降低了试验和测量成本，提高了计算效率，保证了数值模拟的准确性。

[0024] (2) 通过设置渗流条件得到渗流模型，在煤体模拟中选用了拉格朗日算法能够精确地描述结构边界运动；空气场模拟中选用 MAT- VACUUM空隙材料模型，并选择 ALE算法作为运算方法，能够克服因流体穿过煤体骨架结构而引起的自身网格的严重畸变。

[0025] (3) 利用数值模拟的方法对渗流进行模拟，根据模拟结果能够实现对煤体内部渗流分布情况的直观描述，并根据变化过对对渗流变化进行有效的预测。

[0026] (4) 本方法中利用达西定律拟合压力梯度和渗流速度的关系，得到形变煤体结构的渗透率，预测结果准确性高，测试方法的适用范围广。

对附图的简要说明

附图说明

[0027] 图 1是煤样试件形状结构示意图；

[0028] 图 2是 CT扫描结构示意图；

[0029] 图 3是原位拉伸、压缩及温度控制实验装置结构示意图；

[0030] 图 4是煤样试件微观结构示意图；

[0031] 图 5是渗流模型和渗流条件示意图；

[0032] 图 6是裂隙渗流速度曲线图；

[0033] 图 7是孔隙渗流速度曲线图；

[0034] 图中： 1-煤样试件； 2-X射线源； 3 -载物台； 4-CCD探测器； 5-X射线； 6 -刚性压头； 7 -上夹持装置； 8 -上固定螺栓； 9 -下加持装置； 10-下固定螺栓； 11 -煤样试件微观结构； 12 -水场； 13 -空气场。

发明实施例

本发明的实施方式

[0035] 结合图 1至图 7所示，本发明提供的基于 CT扫描的预测应力加载条件下煤体渗透率的方法具体实施方式如下：

[0036] 实施例 1

[0037] 基于 CT扫描的预测应力加载条件下煤体渗透率的方法，包括以下步骤：

[0038] A.制作煤样试件，进行单轴压缩实验，同时对煤样试件进行 CT扫描；

[0039] 其中，煤样试件为圆柱体，所述单轴压缩实验使用原位拉伸、压缩及温度控制实验装置；所述单轴压缩试验中分别进行加载速度不同的多组煤样试件的单轴压缩，同一煤样试件的加载速度恒定。 [0040] 在本步骤中， CT扫描前首先根据 X射线稳定性、煤样试件尺寸、煤样试件 X射线衰减分数和曝光时间来确定扫描电压、扫描功率和视场大小；所述 CT扫描是以恒定的速度旋转扫描，检测器捕捉由 X射线源发出穿过煤样试件的 X射线，以电信号的形式储存 CT扫描图像。

[0041] B.采用 DTM数字地形模型法进行阀值分割，并得到煤样试件的阀值，将 CT扫描图像导入 Simpleware软件建立三维数值模型；在 Simpleware软件建立三维数值模型后进行网格划分。

[0042] 确定煤样试件的阀值包括：将 CT扫描图像转化为数字地形模型，表示孔裂隙体积和煤体总体积，孔裂隙体积与煤体总体积的比值作为孔隙率；建立孔隙率和图像灰度值的函数关系，计算孔隙率函数所有极小值中的最大值，作为 era 描图像的阀值。

[0043] 其中 Simpleware软件是基于核心图形处理平台 ScanIP的图像处理模块，可选择的网格生成 FE Module有限元模块，以及 CAD集成 +CAD模块和 NURBS Module 曲面建模模块，在处理与整合图像、 CAD与仿真技术领域有显著效果。

[0044] C.将三维数值模型导入 HyperMesh软件中并设置渗流条件；其中 Hypermesh软件具有强大的有限元网格划分前处理功能。

[0045] 渗流条件包括水场模拟、空气场模拟和煤体模拟；所述水场模拟选用 MAT-NU LL本构模型，设置流体沿一个方向流动，在流动方向上的两个侧面施加垂直于流动方向的两个方向的约束，在另外四个侧面施加固定约束；通过布置线性多项式状态方程为水头处流体的运动提供压力。空气场模拟中选用 MAT-VACUU M空隙材料模型，并选择 ALE算法作为运算方法。煤体模拟选用拉格朗日算法作为运算方法。

[0046] D.对三维数值模型施加多个压力梯度值 VP，设置初始流动速度 vO, 调整 VP和 vO参数分别模拟，得到渗流模型；

[0047] 压力梯度值 VP的选择范围为 0~300Pa/mm; 所述初始流动速度 v0的选择范围为 0.01~0.09mm/s。在试验中，优选的渗流模型施加的压力梯度值 VP分别为 0Pa/m m、 10Pa/mm、 24Pa/mm和 213Pa/mm; 初始流动速度 vO设置为 0.02mm/s。

[0048] E.将渗流模型导入 Ls-dyna软件中模拟计算，得到计算结果；其中渗流模型为“k ”文件格式，计算后得到“d3plot”文件。 Ls-dyna软件是功能齐全的几何非线性、材料非线性和接触非线性软件，它以 Lagrange算法为主，兼有 ALE和 Euler算法；以结构分析为主，兼有热分析、流体-结构耦合功能；以非线性动力分析为主，兼有静力分析功能的非线性有限元程序。

[0049] F.将计算结果导入 HyperView软件中，沿渗流方向等距提取渗流速度，得到渗流速度和压力梯度的关系，计算渗透率。将“d3plot”文件导入 HyperView软件中，进行后处理，沿模型等距提取渗流模型内部的多个稳定渗流速度。

[0050] 实施例 2

[0051] 本发明提供的基于 CT扫描的预测应力加载条件下煤体渗透率的方法实施步骤包括：使用 nanoVoxel-3502E型扫描仪和原位拉伸、压缩控制实验装置，对多组煤样试件进行不同应力加载条件下的 CT扫描，获取煤样试件压缩过程中的结构 C T扫描图像。将获取的 CT扫描图像通过 Simpleware软件建立三维数值模型，从而将形变煤体结构的 CT扫描图像转化为具有可视化特征的数值模型，以方便对煤体内部渗流分布情况的直观描述。通过 HyperMesh软件进行煤体结构的渗流条件设置，主要包括空气场模拟和水场模拟。通过 Ls-dyna软件进行模拟计算，将计算结果导入到 HyperView软件中提取渗流速度，通过拟合压力梯度和渗流速度记录数据，得到形变煤体结构的渗透率。

[0052] 具体实施步骤包括：

[0053] 步骤 A.首先将煤样试件制作成截面直径 9 mm的煤柱，通过原位拉伸、压缩及温度控制实验装置和 CT扫描试验装置，对煤样试件在进行单轴压缩过程进行 CT 扫描。其中 CT扫描试验装置包括 X射线源 2、载物台 3、 CCD探测器 4和 X射线 5 ; 原位拉伸、压缩及温度控制实验装置包括刚性压头 6、上夹持装置 7、上固定螺栓 8、下加持装置 9和下固定螺栓 10。

[0054] 具体操作过程为：将煤样试件安放在载物台 3上，通过下固定螺栓 10调节下加持装置 9 , 完成对煤样试件的调平；通过上加持装置 7和上固定螺栓 8对煤样试件的顶端进行固定。根据 X射线稳定性、样品尺寸、样品的 X射线衰减分数和曝光时间等因素的影响，对 X射线源 2和 CCD探测器 4进行调节，选择适当的扫描电压、扫描功率和视场大小等实验条件。利用刚性压头 6对煤样试件分别进行单轴压缩加载，加载速度分别选取 Omm/s、 0.001mm/s、 0.002mm/s和 0.003 mm/s; 同时以恒定的扫描速度对载物台 3进行旋转， X射线 5穿过煤柱 1后由 CCD探测器 4进行捕捉，以电信号的形式进行储存并得到 CT扫描图像。

[0055] 步骤 B.使用 DTM法进行阈值分割，确定各煤样的阈值情况。其原理是通过将 C T的灰度图像转化成数字地形模型，基于此对孔裂隙的体积 VE、煤体的总体积 V T进行表征，进而构建孔隙率 cp(x)与灰度值 x之间的变化函数，最终通过求取该函数的极小值中的最大值确定作为 CT扫描图像的阈值。其中孔隙率 cp_(x)与灰度值 _x 的函数关系式为：

[0056] 式中： n为各像素点的灰度值， r_;e[ r_min r_mJ； H( r_;)^I[ r_min r_mJ范围内的灰度直方图，其数值大小等于图像中灰度值为 n的像素的个数与图像像素个数总和的比值。

[0057] 对不同加载速度下获得的 4组 CT扫描图像进行适当裁剪，并将其分别导入到 Si mpleware软件建立三维数值模型，选取适当网格类型对该模型进行网格划分，构建出 4组煤体微细观结构，如图 4所示，煤样试件微观结构示意图，是原始煤体的微细观结构，此时单轴加载速度为 0 mm/s。

[0058] 步骤 C.将构建的三维数值模型煤体微观结构导入 HyperMesh软件中并设置渗流条件。模拟非稳态渗流过程，其中渗流条件模拟包括水场模拟、空气场模拟和煤体模拟。其中水场模拟设置流体并提供动能，空气场模拟和煤体结构耦合接触，空气场是流固耦合的主要场所。

[0059] 在水场模拟中，规定流体只沿 >>方向进行流动，而其余侧面不允许流体流出。

在设置时，对煤体流动方向上的两个侧面施加 ^ z方向上的约束，其它 4个侧面施加 ^ y、 z方向的约束。煤体模拟选用拉格朗日算法作为运算方法，在渗流模拟设置中 Lagrange算法能够精确地描述结构边界运动。由于流体的流动会引起水头和空气场网格的严重变形，所以选择能够克服因流体穿过煤体骨架结构而引起的自身网格严重畸变的 ALE算法作为水头和空气场的运算方法，其控制流体问题的方程是 Navier-Stokes方程的 ALE描述：

应力张量，；?是压力， /是二阶恒等张量；对于牛顿流体，粘性应力张量由

给出，其中 ^是动态粘度； /是质量体力； e是密度。

[0061] 在材料属性设置方面，选择能够使自身结构不发生变形的刚体属性对煤骨架结构进行定义。水场模拟采用适合于流体属性的“MAT-NULL”本构模型；空气场模拟作为流固耦合发生的主要场所，在属性定义时需要同时兼顾空气物理性质以及煤的多孔结构特点，而“MAT-VACUUM”是一种结合 ALE算法的孔隙材料模型，它能够满足模拟的设置要求，因此选其作为空气场的本构模型；通过布置线性多项式状态方程为水头处流体的运动提供压力，其中线性多项式状态方程为：

[0062] 式中： P为流体压力； e _v为内能与初始体积之比；为比体积（^= ^ -1，其中为原始流体密度， P为当前流体密度）；＜^~（：₆为系数；方程参数按照一般压力方程进行取值： C₀= P₀=1.01xl0 ⁵Pa, C₂= C₃= C₄= C₅= C₆=0, （^=2.25x10 ⁹

Pa为水的体积模量。

[0063] 步骤 D.对 4组煤体微细观结构分别进行渗流模拟设置，并施加 0.02 mm/s的初始流动速度 v_Q1，在此基础上设置 v /v v /v V P₃的压力梯度，其中 V AS IO Pa/mm, V _P₂S24 Pa/mm， V _?₃为213 Pa/mm。

[0064] 步骤 E.将建立好的渗流模拟以“K”文件的格式导入到 Ls-dyna软件中进行模拟计算，得到计算结果。

[0065] 步骤 F.

将计算结果导入 HyperView软件中，具体是将生成的“d3plot”文件导入到 HyperVi ew软件中进行后处理。沿模型等距提取内部流体的稳定渗流速度，通过拟合压力梯度与渗流速度的关系，获取形变煤体结构的渗透率。

[0066] 实施例 3

[0067] 以新疆金塔大黄山煤矿 2zwl l工作面气煤煤样为例，对本发明的方法应用及原理做进一步的说明。

[0068] 步骤 A.首先将煤样试件制作成截面直径 9 mm的煤柱，通过原位拉伸、压缩及温度控制实验装置和 CT扫描试验装置，对煤样试件在进行单轴压缩过程进行 CT 扫描。其中 CT扫描试验装置包括 X射线源 2、载物台 3、 CCD探测器 4和 X射线 5 ; 原位拉伸、压缩及温度控制实验装置包括刚性压头 6、上夹持装置 7、上固定螺栓 8、下加持装置 9和下固定螺栓 10。

[0069] 具体操作过程为：将煤样试件安放在载物台 3上，通过下固定螺栓 10调节下加持装置 9 , 完成对煤样试件的调平；通过上加持装置 7和上固定螺栓 8对煤样试件的顶端进行固定。根据 X射线稳定性、样品尺寸、样品的 X射线衰减分数和曝光时间等因素的影响，对 X射线源 2和 CCD探测器 4进行调节，选择扫描电压为 60 kV, 功率为 5W，视场大小为 9.5x9.5 mm ²的实验条件。利用刚性压头 6对煤样试件分别进行单轴压缩加载，加载速度分别选取 Omm/s、 0.001mm/s、 0.002mm/s和 0.003

mm/s；同时以恒定的扫描速度对载物台 3进行旋转， X射线 5穿过煤柱 1后由 CCD 探测器 4进行捕捉，以电信号的形式进行储存并得到 CT扫描图像，经过 4.3 h完成扫描。

[0070] 步骤 B.使用 DTM法进行阈值分割，确定各煤样的阈值为 162Pixel。将获得的 4组 CT图片进行裁剪，得到直径为 0.72 mm的 CT扫描图像，将裁剪后的 CT扫描图像分别导入到 Simpleware软件建立三维数值模型，选取适当网格类型对该模型进行网格划分，构建出 4组煤体微细观结构，如图 4所示，煤样试件微观结构示意图，是原始煤体的微细观结构，此时单轴加载速度为 0 mm/s。 [0071] 步骤 C.将构建的三维数值模型煤体微观结构导入 HyperMesh软件中并设置渗流条件。模拟非稳态渗流过程，其中渗流条件模拟包括水场模拟、空气场模拟和煤体模拟。其中水场模拟设置流体并提供动能，空气场模拟和煤体结构耦合接触，空气场是流固耦合的主要场所。

[0072] 在水场模拟中，规定流体只沿 y方向进行流动，而其余侧面不允许流体流出。

在设置时，对煤体流动方向上的两个侧面施加 ^ z方向上的约束，其它 4个侧面施加 ^ y、 z方向的约束。煤体模拟选用拉格朗日算法作为运算方法，在渗流模拟设置中 Lagrange算法能够精确地描述结构边界运动。由于流体的流动会引起水头和空气场网格的严重变形，所以选择能够克服因流体穿过煤体骨架结构而引起的自身网格严重畸变的 ALE算法作为水头和空气场的运算方法，其控制流体问题的方程是 Na vier- S toke s方程的 ALE描述。

[0073] 煤体模拟选择能够使自身结构不发生变形的刚体属性对煤骨架结构进行定义。

水场模拟采用适合于流体属性的“MAT-NULL”本构模型。空气场模拟选择“MAT -VACUUM”模型。通过布置线性多项式状态方程为水头处流体的运动提供压力

[0074] 步骤 D.对 4组煤体微细观结构进行渗流模拟设置，并施加 0.02 mm/s的初始流动速度 v_Q1，在此基础上设置

, V _P₂为 24 Pa/mm, V _?₃为213 Pa/mm。

[0075] 步骤 E.将建立好的渗流模拟以“K”文件的格式导入到 Ls-dyna软件中进行模拟计算，得到计算结果。

[0076] 步骤 F.通过 Ls-dyna软件进行模拟得到计算结果，将计算结果导入 HyperView软件中，也就是将生成的“d3plot”文件导入到 HyperView软件中进行后处理。利用 H yperView软件沿模型等距提取 3个截面的内部孔隙、裂隙结构中流体速度随时间的变化情况，得到裂隙渗流速度曲线图和孔隙渗流速度曲线图，如图 6和图 7所示。并对 4种渗流模型在压力梯度作用下的稳定渗流速度进行统计，结果如表 1 所示。

[0077] 表 1 压力梯度作用下 4组渗流模型的稳定渗流速度将上述压力梯度、渗流速度数据导入到 Origin软件中，分别采用内置的“线性函数”、 “非线性函数”进行拟合求解，只有“线性函数”的相关系数达到 0.99。由此可以看出，在低压力梯度 10~213 Pa/mm的条件下，随着压力梯度的增加，渗流速度呈线性增大，其关系符合达西定律的描述。通过达西定律对“稳定渗流值-压力梯度”关系进行拟合： y=上. yp

M-

[0078] 式中： v为模型的稳定速度值； VP为压力梯度； k为渗透率; ^为动力粘滞系数。

[0079] 通过计算，分别获得原始模型、 0.001 mm/s、 0.002 mm/s、 0.003 mm/s形变模型结构的渗透率为 37.4 mD、 36.8 mD、 37.1 mD和 37.2 mD。

[0080] 当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。

Claims

权利要求书

[权利要求 1] 基于 CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，包括以下步骤：

A.制作煤样试件，进行单轴压缩实验，同时对煤样试件进行 CT扫描

B.采用数字地形模型法进行阀值分割，并得到煤样试件的阀值，将 C T扫描图像导入 Simpleware软件中建立三维数值模型；

C.将三维数值模型导入 HyperMesh软件中并设置渗流条件；

D.对三维数值模型施加多个压力梯度值▽ P，设置初始流动速度 v_Q，调整▽ P和 ^参数分别模拟，得到渗流模型；

E.将渗流模型导入 Ls-dyna软件中模拟计算，得到计算结果；

F.将计算结果导入 HyperView软件中，沿渗流方向等距提取渗流速度，得到渗流速度和压力梯度的关系，计算渗透率。

[权利要求 2] 根据权利要求 1所述的基于 CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，所述煤样试件为圆柱体，所述单轴压缩实验使用原位拉伸、压缩及温度控制实验装置；所述单轴压缩试验中分别进行加载速度不同的多组煤样试件的单轴压缩，同一煤样试件的加载速度恒定。

[权利要求 3] 根据权利要求 2所述的基于 CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，所述步骤 A中，在 CT扫描前首先根据 X射线稳定性、煤样试件尺寸、煤样试件 X射线衰减分数和曝光时间来确定扫描电压、扫描功率和视场大小；所述 CT扫描是以恒定的速度旋转扫描，检测器捕捉由 X射线源发出穿过煤样试件的 X射线，以电信号的形式储存 CT扫描图像。

[权利要求 4] 根据权利要求 1所述的基于 CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，所述步骤 B中确定煤样试件的阀值包括：将 CT 扫描图像转化为数字地形模型，表示孔裂隙体积和煤体总体积，孔裂隙体积与煤体总体积的比值作为孔隙率；建立孔隙率和图像灰度值的函数关系，计算孔隙率函数所有极小值中的最大值，作为 CT扫描图像的阀值。

[权利要求 5] 根据权利要求 4所述的基于 CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，所述步骤 B中在 Simpleware软件建立三维数值模型后进行网格划分。

[权利要求 6] 根据权利要求 1所述的基于 CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，所述步骤 C中的渗流条件包括水场模拟、空气场模拟和煤体模拟；所述水场模拟选用 MAT-NULL本构模型，设置流体沿一个方向流动，在流动方向上的两个侧面施加垂直于流动方向的两个方向的约束，在另外四个侧面施加固定约束；通过布置线性多项式状态方程为水头处流体的运动提供压力。

[权利要求 7] 根据权利要求 6所述的基于 CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，所述空气场模拟中选用 MAT- VACUUM空隙材料模型，并选择 ALE算法作为运算方法。

[权利要求 8] 根据权利要求 6所述的基于 CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，所述煤体模拟选用拉格朗日算法作为运算方法

[权利要求 9] 根据权利要求 1所述的基于 CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，所述压力梯度值 V P

的选择范围为 0~300Pa/mm; 所述初始流动速度

的选择范围为 0.01~0.09mm/s。

[权利要求 10] 根据权利要求 9所述的基于 CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，所述骤 D中渗流模型施加的压力梯度值 V P分别为 OPa/mm、 10Pa/mm、 24Pa/mm和 213Pa/mm; 初始流动速度 v₀设置为 0.02mm/s。