CN112417663B - 煤矿采动稳定区抽采过程中煤层气运移路径表征方法 - Google Patents
煤矿采动稳定区抽采过程中煤层气运移路径表征方法 Download PDFInfo
- Publication number
- CN112417663B CN112417663B CN202011271744.9A CN202011271744A CN112417663B CN 112417663 B CN112417663 B CN 112417663B CN 202011271744 A CN202011271744 A CN 202011271744A CN 112417663 B CN112417663 B CN 112417663B
- Authority
- CN
- China
- Prior art keywords
- gas
- stable area
- mining
- mining stable
- porosity
- 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.)
- Active
Links
- 238000005065 mining Methods 0.000 title claims abstract description 61
- 239000003245 coal Substances 0.000 title claims abstract description 55
- 238000000605 extraction Methods 0.000 title claims abstract description 41
- 238000013508 migration Methods 0.000 title claims abstract description 32
- 230000005012 migration Effects 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000035699 permeability Effects 0.000 claims abstract description 24
- 239000002245 particle Substances 0.000 claims abstract description 20
- 238000013178 mathematical model Methods 0.000 claims abstract description 6
- 238000011160 research Methods 0.000 claims abstract description 6
- 238000002474 experimental method Methods 0.000 claims abstract description 4
- 230000001052 transient effect Effects 0.000 claims abstract description 4
- 239000004575 stone Substances 0.000 claims description 27
- 239000011148 porous material Substances 0.000 claims description 20
- 238000012360 testing method Methods 0.000 claims description 19
- 230000008859 change Effects 0.000 claims description 18
- 239000011435 rock Substances 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000008092 positive effect Effects 0.000 abstract description 4
- 230000000007 visual effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F7/00—Methods or devices for drawing- off gases with or without subsequent use of the gas for any purpose
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Mining & Mineral Resources (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种煤矿采动稳定区抽采过程中煤层气运移路径表征方法,包括步骤:1)通过实验获得渗透系数K和非达西流因子β关于孔隙率n和粒度d的函数关系式;2)根据目标采动稳定区实际情况,融合渗透系数K与非达西流因子β连续函数建立瓦斯抽采渗流数学模型;3)根据目标采动稳定区现场情况,建立数值模型构建的假设条件;4)选用软件对研究目标建立采动稳定区瓦斯抽采物理数值模型,并在PED模块中输入建立的数学模型;5)对模型进行瞬态解算,得到采动稳定区抽采过程中煤层气在不同时间的运移路径。本发明实现了直观的、符合实际的展示采动稳定区抽采过程中煤层气运移路径,对研究采动稳定区抽采过程中煤层气运移规律有积极作用。
Description
技术领域
本发明涉及煤矿瓦斯抽采技术领域,特别涉及一种可直接应用于煤矿采动稳定区瓦斯抽采井地面位置选取的。
背景技术
在研究采动稳定区抽采过程中煤层气运移规律中,如果能直观的展示采动稳定区抽采过程中煤层气运移路径,能让人能更直观的了解煤层气运移规律,将对研究煤层气运移规律的工作产生积极作用,能帮助选取煤矿采动稳定区瓦斯抽采井的位置。但是如何能符合实际的展示采动稳定区抽采过程中煤层气的运移路径,是一个技术难题。
发明内容
有鉴于此,本发明的目的是提供一种煤矿采动稳定区抽采过程中煤层气运移路径表征方法,以解决目前不能直观的、符合实际的展示采动稳定区抽采过程中煤层气运移路径的技术问题。
本发明煤矿采动稳定区抽采过程中煤层气运移路径表征方法,其包括以下步骤:
(1)连续性方程
(2)动量守恒方程
(2)式中,μ为气体的粘性系数,uns为工作面及巷道中气体的速度矢量;ubr是采动稳定区中气体的速度矢量,ρ为气体的质量密度,pns为工作面及巷道中气体的压力,pbr为采动稳定区气体的压力,F为自定义外部体积力,ρ为气体的质量密度,K为采动稳定区的渗透率,n为采动稳定区的孔隙率;
(3)组分守恒方程
(3)式中,cs、ρcs和Ds分别是组分s的体积浓度、质量浓度和该组分的扩散系数;Ss是系统内部单位时间内单位体积通过化学反应产生的该组分的质量,即生产率;μ是气体的粘性系数,ρ是气体的质量密度;
(4)气体状态方程
(4)式中,M是气体的分子量;R是气体的普适常数;P是气体压力,T是气体温度,ρ是气体的质量密度;
3)根据目标采动稳定区现场情况,建立数值模型构建的假设条件,包括:假设不考虑各组分之间的化学反应,假设混合气体为不可压缩气体;
4)选用COMSOL Multiphysics软件对研究目标建立采动稳定区瓦斯抽采物理数值模型,并在PED模块中输入建立的数学模型;
5)设置合适的边界条件与初始条件后,对模型进行瞬态解算;在结果中选择三维绘图中的体箭头,选择速度v的相关数据,并设置合适的显示参数,得到采动稳定区抽采过程中煤层气在不同时间的运移路径;同时,在结果中选择三维绘图中的切面、面上箭头组合,选择速度v的相关数据,并设置合适的显示参数,得到采动稳定区抽采过程中煤层气在不同时间、不同切面的运移路径。
进一步,所述步骤1)通过实验获得渗透系数K关于孔隙率n和粒度d的函数关系式及非达西流因子β关于孔隙率n和粒度d的函数关系式,其包括:
a)采集目标工作面围岩试样,测出完整岩样的单轴抗压强度与密度,以单轴抗压强度设定碎石承压强度;
b)制备不少于三种粒径的与步骤a)实验结果相符的同一岩性碎石样品,并烘干;
c)对不同粒径碎石样品进行承压渗流试验:通过压缩碎石样品使其孔隙率减少,以模拟采动稳定区不同地点冒落岩石由于受到矿压作用而产生的孔隙率变化,并在压缩碎石样品过程还向碎石样品中通入测试渗流情况的气体;通过改变压头位移量及通入气体的压力,得到若干组轴向载荷、压头位移、进气压力、出气压力、进气口流量和出气口流量的实验数据;
d)试验完成后选取若干组典型数据,做孔隙压力梯度与渗流速度的变化曲线,将拟合得到的孔隙压力梯度与渗流速度的变化曲线与Forchheimer型非达西方程表示的孔隙压力梯度与渗流速度关系曲线进行比较,判断其是否符合Forchheimer型非达西方程;所述Forchheimer型非达西方程如下:
式中,J为孔隙压力梯度,单位为Pa/m;ρ为流体密度,单位为kg/m3;g为重力加速度,单位为N/kg;β为非达西流因子,单位为m-1;v为渗流速度,单位为m/s;K为非达西渗透系数,单位为m/s,其中k为渗透率,单位为m-2;u为流体的动力黏性系数,单位为Pa·s;
e)在确定孔隙压力梯度与渗流速度的变化曲线符合Forchheimer型非达西方程后,将拟合得到的孔隙压力梯度与渗流速度的变化曲线按式(1)处理,得到非达西渗透系数K和非达西流因子β;然后进行二元拟合得到K与承压渗流试验中的碎石样品孔隙率n和碎石样品粒度d的函数关系式f(n,d)以及β与n和d的关系式g(n,d),
本发明的有益效果:
本发明煤矿采动稳定区抽采过程中煤层气运移路径表征方法,其实现了通过COMSOL Multiphysics直观的展示采动稳定区抽采过程中煤层气运移路径;并且其通过碎石承压渗流试验,得到渗流速度、孔隙压力梯度、粒径等实验数据,并对数据进行拟合,得到渗透系数K关于孔隙度和粒径的函数关系式,本方法获得的渗透系数K与实际情况更符合,使得煤层气渗流数学模型更贴合实际情况。本发明方法能够更符合实际的表征煤矿采动稳定区抽采过程中煤层气运移路径,对研究采动稳定区抽采过程中煤层气运移规律有积极作用,能帮助选取煤矿采动稳定区瓦斯抽采井的位置。
附图说明
图1为采动稳定区煤层气抽采模型。
图2为采动稳定区抽采过程中煤层气运移路径图。
图3为采动稳定区抽采过程中某平面煤层气运移路径图。
具体实施方式
本实施例煤矿采动稳定区抽采过程中煤层气运移路径表征方法,其包括以下步骤:
(1)连续性方程
(2)动量守恒方程
(2)式中,μ为气体的粘性系数,uns为工作面及巷道中气体的速度矢量;ubr是采动稳定区中气体的速度矢量,ρ为气体的质量密度,pns为工作面及巷道中气体的压力,pbr为采动稳定区气体的压力,F为自定义外部体积力,ρ为气体的质量密度,K为采动稳定区的渗透率,n为采动稳定区的孔隙率。
(3)组分守恒方程
(3)式中,cs、ρcs和Ds分别是组分s的体积浓度、质量浓度和该组分的扩散系数;Ss是系统内部单位时间内单位体积通过化学反应产生的该组分的质量,即生产率;μ是气体的粘性系数,ρ是气体的质量密度。
(4)气体状态方程
(4)式中,M是气体的分子量;R是气体的普适常数;P是气体压力,T是气体温度,ρ是气体的质量密度。
3)根据目标采动稳定区现场情况,建立数值模型构建的假设条件,包括:假设不考虑各组分之间的化学反应,假设混合气体为不可压缩气体。
4)选用COMSOL Multiphysics软件对研究目标建立采动稳定区瓦斯抽采物理数值模型,并在PED模块中输入建立的数学模型。
5)设置合适的边界条件与初始条件后,对模型进行瞬态解算;在结果中选择三维绘图中的体箭头,选择速度v的相关数据,并设置合适的显示参数,得到采动稳定区抽采过程中煤层气在不同时间的运移路径;同时,在结果中选择三维绘图中的切面、面上箭头组合,选择速度v的相关数据,并设置合适的显示参数,得到采动稳定区抽采过程中煤层气在不同时间、不同切面的运移路径。
本实施例中,所述步骤1)通过实验获得渗透系数K关于孔隙率n和粒度d的函数关系式及非达西流因子β关于孔隙率n和粒度d的函数关系式,其包括:
a)采集目标工作面围岩试样,测出完整岩样的单轴抗压强度与密度,以单轴抗压强度设定碎石承压强度。
b)制备不少于三种粒径的与步骤a)实验结果相符的同一岩性碎石样品,并烘干。
c)对不同粒径碎石样品进行承压渗流试验:通过压缩碎石样品使其孔隙率减少,以模拟采动稳定区不同地点冒落岩石由于受到矿压作用而产生的孔隙率变化,并在压缩碎石样品过程还向碎石样品中通入测试渗流情况的气体;通过改变压头位移量及通入气体的压力,得到若干组轴向载荷、压头位移、进气压力、出气压力、进气口流量和出气口流量的实验数据。本实施例中对不同粒径碎石样品进行承压渗流试验是将碎石样品装入申请号为202010014728.5公开的测试破碎岩石渗透性参数的装置中,并在GCTS试验系统平台进行承压渗流试验。
d)试验完成后选取若干组典型数据,做孔隙压力梯度与渗流速度的变化曲线,将拟合得到的孔隙压力梯度与渗流速度的变化曲线与Forchheimer型非达西方程表示的孔隙压力梯度与渗流速度关系曲线进行比较,判断其是否符合Forchheimer型非达西方程;所述Forchheimer型非达西方程如下:
式中,J为孔隙压力梯度,单位为Pa/m;ρ为流体密度,单位为kg/m3;g为重力加速度,单位为N/kg;β为非达西流因子,单位为m-1;v为渗流速度,单位为m/s;K为非达西渗透系数,单位为m/s,其中k为渗透率,单位为m-2;u为流体的动力黏性系数,单位为Pa·s。
e)在确定孔隙压力梯度与渗流速度的变化曲线符合Forchheimer型非达西方程后,将拟合得到的孔隙压力梯度与渗流速度的变化曲线按式(1)处理,得到非达西渗透系数K和非达西流因子β;然后进行二元拟合得到K与承压渗流试验中的碎石样品孔隙率n和碎石样品粒度d的函数关系式f(n,d)以及β与n和d的关系式g(n,d),
本实施例中煤矿采动稳定区抽采过程中煤层气运移路径表征方法,其实现了通过COMSOL Multiphysics直观的展示采动稳定区抽采过程中煤层气运移路径;并且其通过碎石承压渗流试验,得到渗流速度、孔隙压力梯度、粒径等实验数据,并对数据进行拟合,得到渗透系数K关于孔隙度和粒径的函数关系式,本方法获得的渗透系数K与实际情况更符合,使得煤层气渗流数学模型更贴合实际情况。本发明方法能够更符合实际的表征煤矿采动稳定区抽采过程中煤层气运移路径,对研究采动稳定区抽采过程中煤层气运移规律有积极作用。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (1)
1.煤矿采动稳定区抽采过程中煤层气运移路径表征方法,其特征在于:包括以下步骤:
1)通过实验获得渗透系数K关于孔隙率n和粒度d的函数关系式及非达西流因子β关于孔隙率n和粒度d的函数关系式:
式中K为渗透系数;β为非达西流因子;n为碎石样品的孔隙度;d为碎石样品的粒度;a1、a2和a3为通过试验数据拟合渗透系数K关于孔隙率n和粒度d变化曲线的相关拟合常数;b1、b2和b3为通过试验数据拟合非达西流因子β关于孔隙率n和粒度d变化曲线的相关拟合常数;
2)根据目标采动稳定区实际情况,融合渗透系数K与非达西流因子β连续函数建立瓦斯抽采渗流数学模型,包括:
(1)连续性方程
(1)式中n为孔隙率,ρ为气体的质量密度,V为气体在多孔介质内的渗流速度,q为源强度,▽是拉普拉斯算子;
(2)动量守恒方程
(2)式中,μ为气体的粘性系数,uns为工作面及巷道中气体的速度矢量;ubr是采动稳定区中气体的速度矢量,ρ为气体的质量密度,pns为工作面及巷道中气体的压力,pbr为采动稳定区气体的压力,F为自定义外部体积力,K为采动稳定区的渗透率,n为采动稳定区的孔隙率;
(3)组分守恒方程
(3)式中,cs、ρcs和Ds分别是组分s的体积浓度、质量浓度和该组分的扩散系数;Ss是系统内部单位时间内单位体积通过化学反应产生的该组分的质量,即生产率;μ是气体的粘性系数,ρ是气体的质量密度;
(4)气体状态方程
(4)式中,M是气体的分子量;R是气体的普适常数;P是气体压力,T是气体温度,ρ是气体的质量密度;
3)根据目标采动稳定区现场情况,建立数值模型构建的假设条件,包括:假设不考虑各组分之间的化学反应,假设混合气体为不可压缩气体;
4)选用COMSOL Multiphysics软件对研究目标建立采动稳定区瓦斯抽采物理数值模型,并在PED模块中输入建立的数学模型;
5)设置合适的边界条件与初始条件后,对模型进行瞬态解算;在结果中选择三维绘图中的体箭头,选择速度v的相关数据,并设置合适的显示参数,得到采动稳定区抽采过程中煤层气在不同时间的运移路径;同时,在结果中选择三维绘图中的切面、面上箭头组合,选择速度v的相关数据,并设置合适的显示参数,得到采动稳定区抽采过程中煤层气在不同时间、不同切面的运移路径;
所述步骤1)通过实验获得渗透系数K关于孔隙率n和粒度d的函数关系式及非达西流因子β关于孔隙率n和粒度d的函数关系式,其包括:
a)采集目标工作面围岩试样,测出完整岩样的单轴抗压强度与密度,以单轴抗压强度设定碎石承压强度;
b)制备不少于三种粒径的与步骤a)实验结果相符的同一岩性碎石样品,并烘干;
c)对不同粒径碎石样品进行承压渗流试验:通过压缩碎石样品使其孔隙率减少,以模拟采动稳定区不同地点冒落岩石由于受到矿压作用而产生的孔隙率变化,并在压缩碎石样品过程还向碎石样品中通入测试渗流情况的气体;通过改变压头位移量及通入气体的压力,得到若干组轴向载荷、压头位移、进气压力、出气压力、进气口流量和出气口流量的实验数据;
d)试验完成后选取若干组典型数据,做孔隙压力梯度与渗流速度的变化曲线,将拟合得到的孔隙压力梯度与渗流速度的变化曲线与Forchheimer型非达西方程表示的孔隙压力梯度与渗流速度关系曲线进行比较,判断其是否符合Forchheimer型非达西方程;所述Forchheimer型非达西方程如下:
式中,J为孔隙压力梯度,单位为Pa/m;ρ为流体密度,单位为kg/m3;g为重力加速度,单位为N/kg;β为非达西流因子,单位为m-1;v为渗流速度,单位为m/s;K为非达西渗透系数,单位为m/s,其中k为渗透率,单位为m-2;u为流体的动力黏性系数,单位为Pa·s;
e)在确定孔隙压力梯度与渗流速度的变化曲线符合Forchheimer型非达西方程后,将拟合得到的孔隙压力梯度与渗流速度的变化曲线按式(1)处理,得到非达西渗透系数K和非达西流因子β;然后进行二元拟合得到K与承压渗流试验中的碎石样品孔隙率n和碎石样品粒度d的函数关系式f(n,d)以及β与n和d的关系式g(n,d),
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011271744.9A CN112417663B (zh) | 2020-11-13 | 2020-11-13 | 煤矿采动稳定区抽采过程中煤层气运移路径表征方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011271744.9A CN112417663B (zh) | 2020-11-13 | 2020-11-13 | 煤矿采动稳定区抽采过程中煤层气运移路径表征方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112417663A CN112417663A (zh) | 2021-02-26 |
CN112417663B true CN112417663B (zh) | 2022-09-20 |
Family
ID=74830965
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011271744.9A Active CN112417663B (zh) | 2020-11-13 | 2020-11-13 | 煤矿采动稳定区抽采过程中煤层气运移路径表征方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112417663B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101115906A (zh) * | 2005-02-04 | 2008-01-30 | 烷材料股份有限公司 | 制备支撑剂的组合物及方法 |
CN102735600A (zh) * | 2012-07-05 | 2012-10-17 | 重庆大学 | 真三轴状态下煤样渗流测试方法 |
CN102735547A (zh) * | 2012-07-05 | 2012-10-17 | 重庆大学 | 真三轴状态下煤岩水压致裂试验方法 |
CN103134448A (zh) * | 2013-01-29 | 2013-06-05 | 河南理工大学 | 一种钻孔瓦斯抽采有效半径的测试方法 |
JP2016159264A (ja) * | 2015-03-04 | 2016-09-05 | 新日鐵住金株式会社 | 連続式固定床触媒反応装置およびその装置を用いたガス改質方法 |
CN105980516A (zh) * | 2013-08-02 | 2016-09-28 | 埃莫瑞油田矿产公司 | 包括高岭土的支撑剂和抗回流添加剂 |
CN107462936A (zh) * | 2017-08-28 | 2017-12-12 | 中国石油大学(北京) | 利用压力监测资料反演低渗透储层非达西渗流规律的方法 |
CN207636453U (zh) * | 2017-12-21 | 2018-07-20 | 西安科技大学 | 松散煤体内co2渗透系统测定装置 |
CN111189755A (zh) * | 2020-01-07 | 2020-05-22 | 重庆大学 | 煤矿采动稳定区瓦斯有效储存空间的数值模拟获取方法 |
CN111520131A (zh) * | 2020-04-08 | 2020-08-11 | 中国矿业大学 | 一种超远距离原位测定煤层瓦斯压力装置及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8614171B2 (en) * | 2006-01-04 | 2013-12-24 | Halliburton Energy Services, Inc. | Compositions for stimulating liquid-sensitive subterranean formations |
-
2020
- 2020-11-13 CN CN202011271744.9A patent/CN112417663B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101115906A (zh) * | 2005-02-04 | 2008-01-30 | 烷材料股份有限公司 | 制备支撑剂的组合物及方法 |
CN102735600A (zh) * | 2012-07-05 | 2012-10-17 | 重庆大学 | 真三轴状态下煤样渗流测试方法 |
CN102735547A (zh) * | 2012-07-05 | 2012-10-17 | 重庆大学 | 真三轴状态下煤岩水压致裂试验方法 |
CN103134448A (zh) * | 2013-01-29 | 2013-06-05 | 河南理工大学 | 一种钻孔瓦斯抽采有效半径的测试方法 |
CN105980516A (zh) * | 2013-08-02 | 2016-09-28 | 埃莫瑞油田矿产公司 | 包括高岭土的支撑剂和抗回流添加剂 |
JP2016159264A (ja) * | 2015-03-04 | 2016-09-05 | 新日鐵住金株式会社 | 連続式固定床触媒反応装置およびその装置を用いたガス改質方法 |
CN107462936A (zh) * | 2017-08-28 | 2017-12-12 | 中国石油大学(北京) | 利用压力监测资料反演低渗透储层非达西渗流规律的方法 |
CN207636453U (zh) * | 2017-12-21 | 2018-07-20 | 西安科技大学 | 松散煤体内co2渗透系统测定装置 |
CN111189755A (zh) * | 2020-01-07 | 2020-05-22 | 重庆大学 | 煤矿采动稳定区瓦斯有效储存空间的数值模拟获取方法 |
CN111520131A (zh) * | 2020-04-08 | 2020-08-11 | 中国矿业大学 | 一种超远距离原位测定煤层瓦斯压力装置及方法 |
Non-Patent Citations (2)
Title |
---|
基于非达西渗流的采空区自然发火数值模拟;秦跃平 等;《煤炭学报》;20120731;第37卷(第7期);1177-1183 * |
采空区气体三维多场耦合规律研究;车强;《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅰ辑》;20110115(第01期);B021-13 * |
Also Published As
Publication number | Publication date |
---|---|
CN112417663A (zh) | 2021-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108590601B (zh) | 一种注水扩容膨胀施工参数优选的实验方法 | |
Ma et al. | Compaction and seepage properties of crushed limestone particle mixture: an experimental investigation for Ordovician karst collapse pillar groundwater inrush | |
Chen et al. | Generalized lattice Boltzmann model for flow through tight porous media with Klinkenberg's effect | |
Karacan | Prediction of porosity and permeability of caved zone in longwall gobs | |
Swami et al. | Non-Darcy flow in shale nanopores: do we have a final answer? | |
Ma et al. | An experimental investigation of permeability measurement of water flow in crushed rocks | |
CN108956435B (zh) | 一种高温高压储层溶蚀的模拟实验方法和装置 | |
Pine et al. | The development of a new numerical modelling approach for naturally fractured rock masses | |
Ma et al. | Experimental investigation on hydraulic properties of granular sandstone and mudstone mixtures | |
Fu et al. | Extended finite element method simulation and experimental test on failure behavior of defects under uniaxial compression | |
Reyes et al. | Empirical correlation of effective stress dependent shale rock properties | |
Ma et al. | Effect of particle mixture on seepage properties of crushed mudstones | |
Duan et al. | Discrete element modelling of stress-induced instability of directional drilling boreholes in anisotropic rock | |
Wu et al. | Coupled effect of cement hydration and temperature on hydraulic behavior of cemented tailings backfill | |
Herr et al. | Experimental studies of mass transport in porous media with local heterogeneities | |
He et al. | Parametric study and dimensional analysis on prescribed hydraulic fractures in cave mining | |
Marina et al. | Modelling of hydraulic fracturing process by coupled discrete element and fluid dynamic methods | |
Lin et al. | Numerical study on clogging mechanism of slurry infiltration in porous media based on coupled CFD-DEM method | |
Tan et al. | On non-Darcian flow behavior in a rough-walled fracture filled with a porous medium | |
Pang et al. | Experimental investigation of water-gas mixed seepage parameters characteristics in broken coal medium | |
Karev et al. | Experimental studies of the deformation, destruction and filtration in rocks: a review | |
CN112417663B (zh) | 煤矿采动稳定区抽采过程中煤层气运移路径表征方法 | |
Ayawah et al. | A novel approach for determining cutting geometry for TBM using full-scale laboratory linear rock cutting and PFC3D-based numerical simulations | |
Yu et al. | Numerical simulation study on the seepage characteristics of coal seam infusion effected by mining-induced stress | |
Banerjee et al. | Influence of fluid viscosity and flow transition over non-linear filtration through porous media |
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 |