CN109211666A - 基于ct扫描的预测应力加载条件下煤体渗透率的方法 - Google Patents
基于ct扫描的预测应力加载条件下煤体渗透率的方法 Download PDFInfo
- Publication number
- CN109211666A CN109211666A CN201811009279.4A CN201811009279A CN109211666A CN 109211666 A CN109211666 A CN 109211666A CN 201811009279 A CN201811009279 A CN 201811009279A CN 109211666 A CN109211666 A CN 109211666A
- Authority
- CN
- China
- Prior art keywords
- scan
- coal
- coal body
- test specimen
- model
- 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.)
- Granted
Links
- 239000003245 coal Substances 0.000 title claims abstract description 135
- 238000002591 computed tomography Methods 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000035699 permeability Effects 0.000 title claims abstract description 38
- 238000011068 loading method Methods 0.000 title claims abstract description 30
- 238000012360 testing method Methods 0.000 claims abstract description 57
- 238000007906 compression Methods 0.000 claims abstract description 24
- 230000006835 compression Effects 0.000 claims abstract description 23
- 238000002474 experimental method Methods 0.000 claims abstract description 15
- 238000005325 percolation Methods 0.000 claims abstract description 12
- 230000011218 segmentation Effects 0.000 claims abstract description 7
- 239000000284 extract Substances 0.000 claims abstract description 6
- 238000004088 simulation Methods 0.000 claims description 48
- 239000012530 fluid Substances 0.000 claims description 28
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 20
- 239000000463 material Substances 0.000 claims description 6
- 230000000007 visual effect Effects 0.000 claims description 5
- 239000011800 void material Substances 0.000 claims description 4
- 238000012669 compression test Methods 0.000 claims description 3
- 230000008676 import Effects 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 239000011148 porous material Substances 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 239000011435 rock Substances 0.000 description 6
- 230000008878 coupling Effects 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 5
- 238000005859 coupling reaction Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 241000094396 Bolitoglossa carri Species 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011049 filling Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000012916 structural analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- 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/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/0641—Indicating or recording means; Sensing means using optical, X-ray, ultraviolet, infrared or similar detectors
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Fluid Mechanics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种基于CT扫描的预测应力加载条件下煤体渗透率的方法，解决了变形煤体渗透率计算的技术问题，该方法步骤包括：A.制作煤样试件，进行单轴压缩实验，同时对煤样试件进行CT扫描；B.采用数字地形模型法进行阀值分割，得到煤样试件的阀值，将CT扫描图像导入建模软件建立三维数值模型；C.将三维数值模型导入有限元分析软件中并设置渗流条件；D.对三维数值模型施加多个压力梯度值▽P，设置初始流动速度v0，调整▽P和v0参数，得到多个渗流模型；E.将渗流模型导入有限元软件中模拟计算，得到计算结果；F.将计算结果导入数据处理软件中，沿渗流方向等距提取渗流速度，得到渗流速度和压力梯度的关系，计算渗透率。
Description
技术领域
本发明涉及煤层注水渗流技术领域，尤其是一种基于CT扫描的预测并计算应力加载条件下煤体渗透率的方法。
背景技术
煤矿生产过程中，需要通过煤层注水实现对工作面粉尘污染、煤与瓦斯突出、冲击地压、自然发火、煤体软化等矿井实际生产问题进行综合治理，而注水后的渗流是自然界最常见的物理现象。在一定压差下，岩石允许流体通过的能力叫渗透率，渗透率是表征土或岩石本身传导液体能力的参数，其大小与孔隙度、液体渗透方向上空隙的几何形状、颗粒大小以及排列方向等因素有关，渗透率用来表示渗透性的大小。
传统的渗透率研究方法包括：(1)通过岩石试验系统和真三轴试验系统等对稳定渗流状态下的煤体应力-渗流耦合作用进行研究，在轴压和侧压稳定的情况下，通过测取流量、流动长度、压力梯度等参数，并结合达西渗流公式也可以实现对煤体渗透率的计算。但是由于煤体结构的非均质不透明等特性，所以该方法不能对煤体内部微细观孔裂隙结构的渗流过程进行解释和描述，而复杂的外界环境会影响流量等参数的确定，从而降低渗透率计算结果的准确性。(2)对于非稳态渗流过程中的渗透率变化规律，常采用脉冲衰减法进行研究，其基本工作原理为通过测试岩样一维非稳态渗流过程中孔隙压力随时间的衰减数据，在此基础上，结合数学模型对煤样脉冲渗透率进行测取。但是该方法在测试高渗透煤样时误差很大，这限制了此方法的测试范围。
随着CT扫描技术在岩石领域的大量应用，其具有的高精度、无损伤等特点，可以帮助人们捕获各种储层岩石的孔裂隙空间几何形状，实现对于孔隙率的测定以及孔裂隙结构特征、矿物形态的分析；同时，借助载荷加载装置还可以实现对于孔裂隙结构演化规律的研究。在气体渗流方面，CT扫描技术主要侧重于气体吸附和解吸作用研究；而在注水渗流方面，利用该技术则可以实现渗流通道的构建以及揭示注水渗流过程中微细观结构的损伤破坏机理。但是受到监测手段的制约，目前的CT扫描实验仍然无法实现对渗流通道中流速、压力梯度变化情况进行详细的测量统计，同时复杂的实验测试也增大了研究成本。因此，为进一步的提高煤体渗透率计算的准确性，实现对煤体内部渗流分布情况的描述，有效预测不同应力条件下渗透率的变化情况，需要对现有的技术有待于更进一步的改善和发展。
发明内容
为解决实现直观描述变形煤样内部渗流分布情况，并有效预测不同应力条件下煤体渗透率变化的技术问题，本发明提供了一种基于CT扫描的预测应力加载条件下煤体渗透率的方法，具体技术方案如下：
基于CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，包括以下步骤：
A.制作煤样试件，进行单轴压缩实验，同时对煤样试件进行CT扫描；
B.采用数字地形模型法进行阀值分割，并得到煤样试件的阀值，将CT扫描图像导入Simpleware软件中建立三维数值模型；
C.将三维数值模型导入HyperMesh软件中并设置渗流条件；
D.对三维数值模型施加多个压力梯度值▽P，设置初始流动速度v0，调整▽P和v0参数分别模拟，得到渗流模型；
E.将渗流模型导入Ls-dyna软件中模拟计算，得到计算结果；
F.将计算结果导入HyperView软件中，沿渗流方向等距提取渗流速度，得到渗流速度和压力梯度的关系，计算渗透率。
优选的是，煤样试件为圆柱体，所述单轴压缩实验使用原位拉伸、压缩及温度控制实验装置；所述单轴压缩试验中分别进行加载速度不同的多组煤样试件的单轴压缩，同一煤样试件的加载速度恒定。
进一步优选的是，步骤A中，在CT扫描前首先根据X射线稳定性、煤样试件尺寸、煤样试件X射线衰减分数和曝光时间来确定扫描电压、扫描功率和视场大小；所述CT扫描是以恒定的速度旋转扫描，检测器捕捉由X射线源发出穿过煤样试件的X射线，以电信号的形式储存CT扫描图像。
还优选的是，步骤B中确定煤样试件的阀值包括：将CT扫描图像转化为数字地形模型，表示孔裂隙体积和煤体总体积，孔裂隙体积与煤体总体积的比值作为孔隙率；建立孔隙率和图像灰度值的函数关系，计算孔隙率函数所有极小值中的最大值，作为CT扫描图像的阀值。
还进一步优选的是，步骤B中在Simpleware软件建立三维数值模型后进行网格划分。
还优选的是，步骤C中的渗流条件包括水场模拟、空气场模拟和煤体模拟；所述水场模拟选用MAT-NULL本构模型，设置流体沿一个方向流动，在流动方向上的两个侧面施加垂直于流动方向的两个方向的约束，在另外四个侧面施加固定约束；通过布置线性多项式状态方程为水头处流体的运动提供压力。
还进一步优选的是，空气场模拟中选用MAT-VACUUM空隙材料模型，并选择ALE算法作为运算方法。
还进一步优选的是，煤体模拟选用拉格朗日算法作为运算方法。
还优选的是，压力梯度值▽P的选择范围为0～300Pa/mm；所述初始流动速度v0的选择范围为0.01～0.09mm/s。
还进一步优选的是，骤D中渗流模型施加的压力梯度值▽P分别为0Pa/mm、10Pa/mm、24Pa/mm和213Pa/mm；初始流动速度v0设置为0.02mm/s。
本发明的有益效果包括：
(1)本方法通过室内试验和CT扫描获得煤体的结构特征，并采用数字地形模型法进行阀值分割，得到煤样阀值并建立三维数值模型，从而降低了试验和测量成本，提高了计算效率，保证了数值模拟的准确性。
(2)通过设置渗流条件得到渗流模型，在煤体模拟中选用了拉格朗日算法能够精确地描述结构边界运动；空气场模拟中选用MAT-VACUUM空隙材料模型，并选择ALE算法作为运算方法，能够克服因流体穿过煤体骨架结构而引起的自身网格的严重畸变。
(3)利用数值模拟的方法对渗流进行模拟，根据模拟结果能够实现对煤体内部渗流分布情况的直观描述，并根据变化过对对渗流变化进行有效的预测。
(4)本方法中利用达西定律拟合压力梯度和渗流速度的关系，得到形变煤体结构的渗透率，预测结果准确性高，测试方法的适用范围广。
附图说明
图1是煤样试件形状结构示意图；
图2是CT扫描结构示意图；
图3是原位拉伸、压缩及温度控制实验装置结构示意图；
图4是煤样试件微观结构示意图；
图5是渗流模型和渗流条件示意图；
图6是裂隙渗流速度曲线图；
图7是孔隙渗流速度曲线图；
图中：1-煤样试件；2-X射线源；3-载物台；4-CCD探测器；5-X射线；6-刚性压头；7-上夹持装置；8-上固定螺栓；9-下加持装置；10-下固定螺栓；11-煤样试件微观结构；12-水场；13-空气场。
具体实施方式
结合图1至图7所示，本发明提供的基于CT扫描的预测应力加载条件下煤体渗透率的方法具体实施方式如下：
实施例1
基于CT扫描的预测应力加载条件下煤体渗透率的方法，包括以下步骤：
A.制作煤样试件，进行单轴压缩实验，同时对煤样试件进行CT扫描；
其中，煤样试件为圆柱体，所述单轴压缩实验使用原位拉伸、压缩及温度控制实验装置；所述单轴压缩试验中分别进行加载速度不同的多组煤样试件的单轴压缩，同一煤样试件的加载速度恒定。
在本步骤中，CT扫描前首先根据X射线稳定性、煤样试件尺寸、煤样试件X射线衰减分数和曝光时间来确定扫描电压、扫描功率和视场大小；所述CT扫描是以恒定的速度旋转扫描，检测器捕捉由X射线源发出穿过煤样试件的X射线，以电信号的形式储存CT扫描图像。
B.采用DTM数字地形模型法进行阀值分割，并得到煤样试件的阀值，将CT扫描图像导入Simpleware软件建立三维数值模型；在Simpleware软件建立三维数值模型后进行网格划分。
确定煤样试件的阀值包括：将CT扫描图像转化为数字地形模型，表示孔裂隙体积和煤体总体积，孔裂隙体积与煤体总体积的比值作为孔隙率；建立孔隙率和图像灰度值的函数关系，计算孔隙率函数所有极小值中的最大值，作为CT扫描图像的阀值。
其中Simpleware软件是基于核心图形处理平台ScanIP的图像处理模块，可选择的网格生成FE Module有限元模块，以及CAD集成+CAD模块和NURBS Module曲面建模模块，在处理与整合图像、CAD与仿真技术领域有显著效果。
C.将三维数值模型导入HyperMesh软件中并设置渗流条件；其中Hypermesh软件具有强大的有限元网格划分前处理功能。
渗流条件包括水场模拟、空气场模拟和煤体模拟；所述水场模拟选用MAT-NULL本构模型，设置流体沿一个方向流动，在流动方向上的两个侧面施加垂直于流动方向的两个方向的约束，在另外四个侧面施加固定约束；通过布置线性多项式状态方程为水头处流体的运动提供压力。空气场模拟中选用MAT-VACUUM空隙材料模型，并选择ALE算法作为运算方法。煤体模拟选用拉格朗日算法作为运算方法。
D.对三维数值模型施加多个压力梯度值▽P，设置初始流动速度v0，调整▽P和v0参数分别模拟，得到渗流模型；
压力梯度值▽P的选择范围为0～300Pa/mm；所述初始流动速度v0的选择范围为0.01～0.09mm/s。在试验中，优选的渗流模型施加的压力梯度值▽P分别为0Pa/mm、10Pa/mm、24Pa/mm和213Pa/mm；初始流动速度v0设置为0.02mm/s。
E.将渗流模型导入Ls-dyna软件中模拟计算，得到计算结果；其中渗流模型为“k”文件格式，计算后得到“d3plot”文件。Ls-dyna软件是功能齐全的几何非线性、材料非线性和接触非线性软件，它以Lagrange算法为主，兼有ALE和Euler算法；以结构分析为主，兼有热分析、流体-结构耦合功能；以非线性动力分析为主，兼有静力分析功能的非线性有限元程序。
F.将计算结果导入HyperView软件中，沿渗流方向等距提取渗流速度，得到渗流速度和压力梯度的关系，计算渗透率。将“d3plot”文件导入HyperView软件中，进行后处理，沿模型等距提取渗流模型内部的多个稳定渗流速度。
实施例2
本发明提供的基于CT扫描的预测应力加载条件下煤体渗透率的方法实施步骤包括：使用nanoVoxel-3502E型扫描仪和原位拉伸、压缩控制实验装置，对多组煤样试件进行不同应力加载条件下的CT扫描，获取煤样试件压缩过程中的结构CT扫描图像。将获取的CT扫描图像通过Simpleware软件建立三维数值模型，从而将形变煤体结构的CT扫描图像转化为具有可视化特征的数值模型，以方便对煤体内部渗流分布情况的直观描述。通过HyperMesh软件进行煤体结构的渗流条件设置，主要包括空气场模拟和水场模拟。通过Ls-dyna软件进行模拟计算，将计算结果导入到HyperView软件中提取渗流速度，通过拟合压力梯度和渗流速度记录数据，得到形变煤体结构的渗透率。
具体实施步骤包括：
步骤A.首先将煤样试件制作成截面直径9mm的煤柱，通过原位拉伸、压缩及温度控制实验装置和CT扫描试验装置，对煤样试件在进行单轴压缩过程进行CT扫描。其中CT扫描试验装置包括X射线源2、载物台3、CCD探测器4和X射线5；原位拉伸、压缩及温度控制实验装置包括刚性压头6、上夹持装置7、上固定螺栓8、下加持装置9和下固定螺栓10。
具体操作过程为：将煤样试件安放在载物台3上，通过下固定螺栓10调节下加持装置9，完成对煤样试件的调平；通过上加持装置7和上固定螺栓8对煤样试件的顶端进行固定。根据X射线稳定性、样品尺寸、样品的X射线衰减分数和曝光时间等因素的影响，对X射线源2和CCD探测器4进行调节，选择适当的扫描电压、扫描功率和视场大小等实验条件。利用刚性压头6对煤样试件分别进行单轴压缩加载，加载速度分别选取0mm/s、0.001mm/s、0.002mm/s和0.003mm/s；同时以恒定的扫描速度对载物台3进行旋转，X射线5穿过煤柱1后由CCD探测器4进行捕捉，以电信号的形式进行储存并得到CT扫描图像。
步骤B.使用DTM法进行阈值分割，确定各煤样的阈值情况。其原理是通过将CT的灰度图像转化成数字地形模型，基于此对孔裂隙的体积VE、煤体的总体积VT进行表征，进而构建孔隙率φ(x)与灰度值x之间的变化函数，最终通过求取该函数的极小值中的最大值确定作为CT扫描图像的阈值。其中孔隙率φ(x)与灰度值x的函数关系式为：
式中：ri为各像素点的灰度值，ri∈[rmin，rmax]；H(ri)为[rmin，rmax]范围内的灰度直方图，其数值大小等于图像中灰度值为ri的像素的个数与图像像素个数总和的比值。
对不同加载速度下获得的4组CT扫描图像进行适当裁剪，并将其分别导入到Simpleware软件建立三维数值模型，选取适当网格类型对该模型进行网格划分，构建出4组煤体微细观结构，如图4所示，煤样试件微观结构示意图，是原始煤体的微细观结构，此时单轴加载速度为0mm/s。
步骤C.将构建的三维数值模型煤体微观结构导入HyperMesh软件中并设置渗流条件。模拟非稳态渗流过程，其中渗流条件模拟包括水场模拟、空气场模拟和煤体模拟。其中水场模拟设置流体并提供动能，空气场模拟和煤体结构耦合接触，空气场是流固耦合的主要场所。
在水场模拟中，规定流体只沿y方向进行流动，而其余侧面不允许流体流出。在设置时，对煤体流动方向上的两个侧面施加x、z方向上的约束，其它4个侧面施加x、y、z方向的约束。煤体模拟选用拉格朗日算法作为运算方法，在渗流模拟设置中Lagrange算法能够精确地描述结构边界运动。由于流体的流动会引起水头和空气场网格的严重变形，所以选择能够克服因流体穿过煤体骨架结构而引起的自身网格严重畸变的ALE算法作为水头和空气场的运算方法，其控制流体问题的方程是Navier-Stokes方程的ALE描述：
式中：μ是速度矢量；σ＝(-pI+τ(μ))应力张量，p是压力，I是二阶恒等张量；对于牛顿流体，粘性应力张量由τ(u)＝μ(▽u+(▽u)T)给出，其中μ是动态粘度；f是质量体力；ρ是密度。
在材料属性设置方面，选择能够使自身结构不发生变形的刚体属性对煤骨架结构进行定义。水场模拟采用适合于流体属性的“MAT-NULL”本构模型；空气场模拟作为流固耦合发生的主要场所，在属性定义时需要同时兼顾空气物理性质以及煤的多孔结构特点，而“MAT-VACUUM”是一种结合ALE算法的孔隙材料模型，它能够满足模拟的设置要求，因此选其作为空气场的本构模型；通过布置线性多项式状态方程为水头处流体的运动提供压力，其中线性多项式状态方程为：
P＝(C0+C1μ+C2μ2+C3μ3)+(C4+C5μ+C6μ2)ev
式中：P为流体压力；ev为内能与初始体积之比；μ为比体积(μ＝ρ/ρ0-1，其中ρ0为原始流体密度，ρ为当前流体密度)；C0～C6为系数；方程参数按照一般压力方程进行取值：C0＝P0＝1.01×105Pa，C2＝C3＝C4＝C5＝C6＝0，C1＝2.25×109Pa为水的体积模量。
步骤D.对4组煤体微细观结构分别进行渗流模拟设置，并施加0.02mm/s的初始流动速度v01，在此基础上设置▽P1、▽P2、▽P3的压力梯度，其中▽P1为10Pa/mm，▽P2为24Pa/mm，▽P3为213Pa/mm。
步骤E.将建立好的渗流模拟以“K”文件的格式导入到Ls-dyna软件中进行模拟计算，得到计算结果。
步骤F.将计算结果导入HyperView软件中，具体是将生成的“d3plot”文件导入到HyperView软件中进行后处理。沿模型等距提取内部流体的稳定渗流速度，通过拟合压力梯度与渗流速度的关系，获取形变煤体结构的渗透率。
实施例3
以新疆金塔大黄山煤矿2zw11工作面气煤煤样为例，对本发明的方法应用及原理做进一步的说明。
步骤A.首先将煤样试件制作成截面直径9mm的煤柱，通过原位拉伸、压缩及温度控制实验装置和CT扫描试验装置，对煤样试件在进行单轴压缩过程进行CT扫描。其中CT扫描试验装置包括X射线源2、载物台3、CCD探测器4和X射线5；原位拉伸、压缩及温度控制实验装置包括刚性压头6、上夹持装置7、上固定螺栓8、下加持装置9和下固定螺栓10。
具体操作过程为：将煤样试件安放在载物台3上，通过下固定螺栓10调节下加持装置9，完成对煤样试件的调平；通过上加持装置7和上固定螺栓8对煤样试件的顶端进行固定。根据X射线稳定性、样品尺寸、样品的X射线衰减分数和曝光时间等因素的影响，对X射线源2和CCD探测器4进行调节，选择扫描电压为60kV，功率为5W，视场大小为9.5×9.5mm2的实验条件。利用刚性压头6对煤样试件分别进行单轴压缩加载，加载速度分别选取0mm/s、0.001mm/s、0.002mm/s和0.003mm/s；同时以恒定的扫描速度对载物台3进行旋转，X射线5穿过煤柱1后由CCD探测器4进行捕捉，以电信号的形式进行储存并得到CT扫描图像，经过4.3h完成扫描。
步骤B.使用DTM法进行阈值分割，确定各煤样的阈值为162Pixel。将获得的4组CT图片进行裁剪，得到直径为0.72mm的CT扫描图像，将裁剪后的CT扫描图像分别导入到Simpleware软件建立三维数值模型，选取适当网格类型对该模型进行网格划分，构建出4组煤体微细观结构，如图4所示，煤样试件微观结构示意图，是原始煤体的微细观结构，此时单轴加载速度为0mm/s。
步骤C.将构建的三维数值模型煤体微观结构导入HyperMesh软件中并设置渗流条件。模拟非稳态渗流过程，其中渗流条件模拟包括水场模拟、空气场模拟和煤体模拟。其中水场模拟设置流体并提供动能，空气场模拟和煤体结构耦合接触，空气场是流固耦合的主要场所。
在水场模拟中，规定流体只沿y方向进行流动，而其余侧面不允许流体流出。在设置时，对煤体流动方向上的两个侧面施加x、z方向上的约束，其它4个侧面施加x、y、z方向的约束。煤体模拟选用拉格朗日算法作为运算方法，在渗流模拟设置中Lagrange算法能够精确地描述结构边界运动。由于流体的流动会引起水头和空气场网格的严重变形，所以选择能够克服因流体穿过煤体骨架结构而引起的自身网格严重畸变的ALE算法作为水头和空气场的运算方法，其控制流体问题的方程是Navier-Stokes方程的ALE描述。
煤体模拟选择能够使自身结构不发生变形的刚体属性对煤骨架结构进行定义。水场模拟采用适合于流体属性的“MAT-NULL”本构模型。空气场模拟选择“MAT-VACUUM”模型。通过布置线性多项式状态方程为水头处流体的运动提供压力。
步骤D.对4组煤体微细观结构进行渗流模拟设置，并施加0.02mm/s的初始流动速度v01，在此基础上设置▽P1、▽P2、▽P3的压力梯度，其中▽P1为10Pa/mm，▽P2为24Pa/mm，▽P3为213Pa/mm。
步骤E.将建立好的渗流模拟以“K”文件的格式导入到Ls-dyna软件中进行模拟计算，得到计算结果。
步骤F.通过Ls-dyna软件进行模拟得到计算结果，将计算结果导入HyperView软件中，也就是将生成的“d3plot”文件导入到HyperView软件中进行后处理。利用HyperView软件沿模型等距提取3个截面的内部孔隙、裂隙结构中流体速度随时间的变化情况，得到裂隙渗流速度曲线图和孔隙渗流速度曲线图，如图6和图7所示。并对4种渗流模型在压力梯度作用下的稳定渗流速度进行统计，结果如表1所示。
表1 压力梯度作用下4组渗流模型的稳定渗流速度
将上述压力梯度、渗流速度数据导入到Origin软件中，分别采用内置的“线性函数”、“非线性函数”进行拟合求解，只有“线性函数”的相关系数达到0.99。由此可以看出，在低压力梯度10～213Pa/mm的条件下，随着压力梯度的增加，渗流速度呈线性增大，其关系符合达西定律的描述。通过达西定律对“稳定渗流值-压力梯度”关系进行拟合：
式中：v为模型的稳定速度值；▽P为压力梯度；k为渗透率；μ为动力粘滞系数。
通过计算，分别获得原始模型、0.001mm/s、0.002mm/s、0.003mm/s形变模型结构的渗透率为37.4mD、36.8mD、37.1mD和37.2mD。
当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。
Claims (10)
1.基于CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，包括以下步骤：
A.制作煤样试件，进行单轴压缩实验，同时对煤样试件进行CT扫描；
B.采用数字地形模型法进行阀值分割，并得到煤样试件的阀值，将CT扫描图像导入Simpleware软件中建立三维数值模型；
C.将三维数值模型导入HyperMesh软件中并设置渗流条件；
D.对三维数值模型施加多个压力梯度值▽P，设置初始流动速度v0，调整▽P和v0参数分别模拟，得到渗流模型；
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扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，所述压力梯度值▽P的选择范围为0～300Pa/mm；所述初始流动速度v0的选择范围为0.01～0.09mm/s。
10.根据权利要求9所述的基于CT扫描的预测应力加载条件下煤体渗透率的方法，其特征在于，所述骤D中渗流模型施加的压力梯度值▽P分别为0Pa/mm、10Pa/mm、24Pa/mm和213Pa/mm；初始流动速度v0设置为0.02mm/s。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811009279.4A CN109211666B (zh) | 2018-08-31 | 2018-08-31 | 基于ct扫描的预测应力加载条件下煤体渗透率的方法 |
| PCT/CN2018/125959 WO2020042512A1 (zh) | 2018-08-31 | 2018-12-30 | 基于ct扫描的预测应力加载条件下煤体渗透率的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811009279.4A CN109211666B (zh) | 2018-08-31 | 2018-08-31 | 基于ct扫描的预测应力加载条件下煤体渗透率的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109211666A true CN109211666A (zh) | 2019-01-15 |
| CN109211666B CN109211666B (zh) | 2019-12-03 |
Family
ID=64985765
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811009279.4A Active CN109211666B (zh) | 2018-08-31 | 2018-08-31 | 基于ct扫描的预测应力加载条件下煤体渗透率的方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN109211666B (zh) |
| WO (1) | WO2020042512A1 (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110126058A (zh) * | 2019-05-20 | 2019-08-16 | 重庆大学 | 一种基于ct可视化和3d打印的岩样制备方法 |
| CN110146525A (zh) * | 2019-04-24 | 2019-08-20 | 山东科技大学 | 一种基于分形理论及ct扫描的煤体孔渗参数预测方法 |
| CN110162900A (zh) * | 2019-05-28 | 2019-08-23 | 辽宁工程技术大学 | 一种模拟煤层注水渗流过程的方法 |
| CN110702570A (zh) * | 2019-09-24 | 2020-01-17 | 山东科技大学 | 一种实现煤体孔裂隙动态渗流过程可视化的方法 |
| CN113281232A (zh) * | 2020-11-25 | 2021-08-20 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种砂浆试件抗渗压强的计算方法 |
| CN113863400A (zh) * | 2021-12-02 | 2021-12-31 | 中建安装集团有限公司 | 一种水泥土复合桩变形检测系统及检测方法 |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112326524B (zh) * | 2020-10-22 | 2022-07-19 | 中国石油大学(华东) | 一种基于ct扫描图像的岩石孔渗测量方法 |
| CN113495048B (zh) * | 2021-09-08 | 2021-11-16 | 南通华兴石油仪器有限公司 | 一种高压渗流试验装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104535426A (zh) * | 2014-12-04 | 2015-04-22 | 中国科学院武汉岩土力学研究所 | Ct实时扫描的三轴应力、渗流、化学耦合流变试验系统 |
| CN105510204A (zh) * | 2015-12-03 | 2016-04-20 | 大连理工大学 | 一种基于ct图像的渗透率预测方法 |
| CN107036911A (zh) * | 2017-05-17 | 2017-08-11 | 绍兴文理学院 | 一种ct实时三维扫描的渗流、剪切耦合岩石三轴试验系统 |
| CN107462936A (zh) * | 2017-08-28 | 2017-12-12 | 中国石油大学(北京) | 利用压力监测资料反演低渗透储层非达西渗流规律的方法 |
| WO2018084506A2 (ko) * | 2016-11-01 | 2018-05-11 | 한국지질자원연구원 | 가도리니움 화합물을 반응시킨 암석 시료의 전자현미경 이미지를 활용한 공극률 측정방법 |
| CN108362623A (zh) * | 2018-02-09 | 2018-08-03 | 河海大学 | 一种基于μCT扫描的微观岩石耦合渗透试验装置 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9134457B2 (en) * | 2009-04-08 | 2015-09-15 | Schlumberger Technology Corporation | Multiscale digital rock modeling for reservoir simulation |
| CN104462654B (zh) * | 2014-11-11 | 2017-06-13 | 中国矿业大学 | 浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法 |
| CN105717016B (zh) * | 2014-12-05 | 2019-02-15 | 中国石油天然气股份有限公司 | 一种煤样渗透率的计算方法 |
| CN106202695B (zh) * | 2016-07-07 | 2018-04-20 | 清能艾科（深圳）能源技术有限公司 | 一种采用数字岩心模拟计算岩心渗透率的方法 |
| CN106960070B (zh) * | 2016-12-28 | 2020-02-21 | 山东科技大学 | 一种基于有限元-离散元ct重构煤体的渗流模拟方法 |
-
2018
- 2018-08-31 CN CN201811009279.4A patent/CN109211666B/zh active Active
- 2018-12-30 WO PCT/CN2018/125959 patent/WO2020042512A1/zh active Application Filing
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104535426A (zh) * | 2014-12-04 | 2015-04-22 | 中国科学院武汉岩土力学研究所 | Ct实时扫描的三轴应力、渗流、化学耦合流变试验系统 |
| CN105510204A (zh) * | 2015-12-03 | 2016-04-20 | 大连理工大学 | 一种基于ct图像的渗透率预测方法 |
| WO2018084506A2 (ko) * | 2016-11-01 | 2018-05-11 | 한국지질자원연구원 | 가도리니움 화합물을 반응시킨 암석 시료의 전자현미경 이미지를 활용한 공극률 측정방법 |
| CN107036911A (zh) * | 2017-05-17 | 2017-08-11 | 绍兴文理学院 | 一种ct实时三维扫描的渗流、剪切耦合岩石三轴试验系统 |
| CN107462936A (zh) * | 2017-08-28 | 2017-12-12 | 中国石油大学(北京) | 利用压力监测资料反演低渗透储层非达西渗流规律的方法 |
| CN108362623A (zh) * | 2018-02-09 | 2018-08-03 | 河海大学 | 一种基于μCT扫描的微观岩石耦合渗透试验装置 |
Non-Patent Citations (1)
| Title |
|---|
| 王刚等: ""基于CT三维重建的煤层气非达西渗流数值模拟"", 《煤炭学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110146525A (zh) * | 2019-04-24 | 2019-08-20 | 山东科技大学 | 一种基于分形理论及ct扫描的煤体孔渗参数预测方法 |
| CN110126058A (zh) * | 2019-05-20 | 2019-08-16 | 重庆大学 | 一种基于ct可视化和3d打印的岩样制备方法 |
| CN110162900A (zh) * | 2019-05-28 | 2019-08-23 | 辽宁工程技术大学 | 一种模拟煤层注水渗流过程的方法 |
| CN110702570A (zh) * | 2019-09-24 | 2020-01-17 | 山东科技大学 | 一种实现煤体孔裂隙动态渗流过程可视化的方法 |
| CN113281232A (zh) * | 2020-11-25 | 2021-08-20 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种砂浆试件抗渗压强的计算方法 |
| CN113281232B (zh) * | 2020-11-25 | 2022-03-01 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种砂浆试件抗渗压强的计算方法 |
| CN113863400A (zh) * | 2021-12-02 | 2021-12-31 | 中建安装集团有限公司 | 一种水泥土复合桩变形检测系统及检测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2020042512A1 (zh) | 2020-03-05 |
| CN109211666B (zh) | 2019-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109211666B (zh) | 基于ct扫描的预测应力加载条件下煤体渗透率的方法 | |
| Ishutov et al. | 3D printing sandstone porosity models | |
| Ju et al. | Visualization and transparentization of the structure and stress field of aggregated geomaterials through 3D printing and photoelastic techniques | |
| Xiong et al. | Nonlinear flow behavior through rough-walled rock fractures: the effect of contact area | |
| CN108763711B (zh) | 一种基于岩心扫描图像分块数值模拟的渗透率预测方法 | |
| Ju et al. | Visualization of the complex structure and stress field inside rock by means of 3D printing technology | |
| CN207379851U (zh) | 基于透明土模型的循环荷载加载试验装置 | |
| Wang et al. | Meso-damage cracking characteristics analysis for rock and soil aggregate with CT test | |
| Ju et al. | Lattice-Boltzmann simulation of microscale CH 4 flow in porous rock subject to force-induced deformation | |
| Tao et al. | Quantitative characterization of irregular microfracture network and its effect on the permeability of porous media | |
| CN107808049B (zh) | 基于多孔介质三维微观结构模型的dnapl迁移数值模拟方法 | |
| CN109632429A (zh) | 一种土石混合体双轴压缩试验的试样制备方法 | |
| Zhang et al. | High-order streamline simulation and macro-scale visualization experimental studies on waterflooding under given pressure boundaries | |
| Zhou et al. | Investigation of microscopic pore structure and permeability prediction in sand-conglomerate reservoirs | |
| CN105651964B (zh) | 一种确定裂隙岩体表征单元体积的方法 | |
| Butlanska | Cone penetration test in a virtual calibration chamber | |
| CN109696541A (zh) | 模拟含空洞地层隧道开挖沉降的三维模型隧道试验装置 | |
| CN105866376A (zh) | 实时量测开挖放坡滑坡坡体位移应力的试验方法 | |
| CN105717021B (zh) | 一种定量确定裂隙岩体渗透各向异性的方法 | |
| CN115730486A (zh) | 一种含裂隙煤体的有效导热系数测算方法 | |
| YANG et al. | Three-dimensional simulation of pore scale fluid flow in granular ore media with realistic geometry | |
| Wang et al. | Experimental study on the mechanical behavior and failure characteristics of rock analogs with filled internal fractures: A new method by sand powder 3D printing | |
| CN111323356B (zh) | 一种基于数字化图像处理技术测量岩石裂隙开度与流速矢量的方法 | |
| Yuan et al. | Experimental characterizations of contact movement in two-dimensional rod assembly subjected to direct shearing | |
| Li et al. | Meso-level and macro-level mechanical properties of slip zone soil with varying coarse grain contents |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |