CN112883542B - 考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法 - Google Patents
考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法 Download PDFInfo
- Publication number
- CN112883542B CN112883542B CN202110033385.1A CN202110033385A CN112883542B CN 112883542 B CN112883542 B CN 112883542B CN 202110033385 A CN202110033385 A CN 202110033385A CN 112883542 B CN112883542 B CN 112883542B
- Authority
- CN
- China
- Prior art keywords
- permeability
- loading
- coal
- rock
- fracture
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 164
- 230000035699 permeability Effects 0.000 title claims abstract description 151
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000005065 mining Methods 0.000 title claims abstract description 35
- 239000003245 coal Substances 0.000 claims abstract description 123
- 238000004088 simulation Methods 0.000 claims abstract description 38
- 238000012360 testing method Methods 0.000 claims abstract description 18
- 238000005452 bending Methods 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000000926 separation method Methods 0.000 claims description 10
- 238000005192 partition Methods 0.000 claims description 5
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 4
- 239000010410 layer Substances 0.000 description 11
- 230000006835 compression Effects 0.000 description 7
- 238000007906 compression Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 239000011241 protective layer Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000004451 qualitative analysis Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000004445 quantitative analysis 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
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法，包括以下步骤：a、实测获得卸压开采各岩层力学参数，并进行相似模拟配比；b、采用相似模拟开采获得应力演化特征；c、根据卸压开采覆岩损伤度将覆岩分为破碎煤岩体、裂隙煤岩体和弹性煤岩体；d、根据卸压开采过程中不同区域的应力演化特征获得破碎煤岩体加载阶段，e、采用渗透率测试系统分别进行不同损伤程度煤岩体的应力‑渗透率测试，获得破碎煤岩体加载渗透率公式、裂隙煤岩体第一次加卸载渗透率公式和第二次加载渗透率公式、完整岩体的第一次加卸载渗透率公式；f、构建卸压开采数值模型更新渗透率。
Description
技术领域
本发明涉及煤层开采领域，具体涉及一种考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法。
背景技术
对于高瓦斯煤层卸压开采覆岩渗透率发育情况的研究主要采用现场实测，实验室相似模拟及数值模拟研究。现场实测一般视为最直接可靠的描述覆岩裂隙及渗透率演化规律的手段，国内外均有大量使用。但是在很多情况下由于费用高，安全性差，可行性低等原因，不能全方位，大面积的进行测量，只能进行定点定性的分析。国内外很多学者采用实验室相似模拟或者数值模拟等手段结合现场观测进行覆岩裂隙发育及渗透率分布特征的研究。但数值模拟研究同样存在多种假设，参数的合理选取对模拟结果至关重要，因此，实验室不同损伤程度煤样应力渗透率研究成为数值模拟参数选取的基础。
综上可以看出，对于卸压开采过程中覆岩渗透率演化规律大多仍处于定性分析阶段，缺乏进一步的定量描述。现有的数值模拟方法主要通过围岩破坏情况间接的反映保护层开采对本煤层和被保护层瓦斯渗流的影响，或者根据经验直接设置渗透率的分布情况。
综上所述，现有技术中存在以下问题：深部高瓦斯煤层卸压开采过程中应力路径复杂多变、各区域损伤程度差别大，渗透率无法量化计算。
发明内容
本发明提供一种考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法，以解决深部高瓦斯煤层卸压开采过程中应力路径复杂多变、各区域损伤程度差别大、渗透率无法量化计算的问题。
为此，本发明提出一种考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法，包括以下步骤：
步骤a、实测获得卸压开采各岩层力学参数，并进行相似模拟配比；
步骤b、采用相似模拟开采获得卸压开采过程中各区域的应力演化特征；
步骤c、根据卸压开采覆岩损伤度将覆岩分为破碎煤岩体、裂隙煤岩体和弹性煤岩体，破碎煤岩体位于垮落带，裂隙煤岩体位于裂隙带，弹性煤岩体位于弯曲下沉带；
步骤d、根据卸压开采过程中不同区域的应力演化特征设计渗透率测试的应力路径，获得破碎煤岩体加载阶段，其中裂隙煤岩体为卸载-加载阶段，弹性煤岩体为加载-卸载阶段；
步骤e、采用渗透率测试系统分别进行不同损伤程度煤岩体的应力-渗透率测试，获得破碎煤岩体加载渗透率公式、裂隙煤岩体第一次加卸载渗透率公式和第二次加载渗透率公式、完整岩体的第一次加卸载渗透率公式；
步骤f、构建卸压开采数值模型，在模拟过程中根据获得的不同应力状态、不同损伤程度煤岩体的渗透率公式更新渗透率。
进一步地，步骤f中，在模拟过程中利用FISH语言更新渗透率。
进一步地，步骤f中，在模拟过程中利用FISH语言更新渗透率具体包括：
步骤f1、首先判断卸压开采过程中覆岩所处分区；
步骤f2、判断为弯曲下沉带的块体进一步判断块体的膨胀变形，当膨胀变形大于等于0.3％时，则认为该区域处于离层裂隙带；
步骤f3、损伤程度判断完毕后，进一步判断块体所处的加卸载状态；
步骤f4、加卸载判断完毕后，进一步判断块体所经历的加卸载次数；
步骤f5、判断为垮落带的块体的水平渗透率和垂直渗透率均采用破碎煤岩体渗透率模型计算；
步骤f6、判断为裂隙带的块体的水平渗透率和垂直渗透率均采用裂隙煤岩体渗透率模型计算；
步骤f7、判断为离层裂隙带的块体的水平渗透率采用裂隙煤岩体渗透率模型计算，垂直渗透率采用弹性煤岩体渗透率计算模型；
步骤f8、膨胀变形小于0.3％的弹性块体的水平渗透率和垂直渗透率采用弹性煤岩体渗透率计算模型。
进一步地，步骤f1中，垮落带分区根据采高和直接顶系数在模拟过程中直接设置；裂隙带分区则根据模拟过程中的弹塑性状态判断：其中，弹性区属于弯曲下沉带，塑性区属于裂隙带。
进一步地，步骤f1中，垮落带分区确定采用如下方法：
f11、根据直接顶单轴抗压强度对直接顶分成3类，分别是：Ⅰ类，单轴抗压强度大于40MPa；Ⅱ类，单轴抗压强度在20～40MPa；和Ⅲ类，单轴抗压强度小于20MPa；
f12、Ⅰ类直接顶c1，c2和c3分别为2.1、16和2.5
f13、Ⅱ类直接顶c1，c2和c3分别为4.7、19和2.2
f14、Ⅲ类直接顶c1，c2和c3分别为6.2、32和1.5
f15、根据公式
计算获得垮落带范围高度H1。
进一步地，具体包括以下步骤：
步骤S1、实测获得卸压开采各岩层力学参数，并进行相似模拟配比；
步骤S2、采用相似模拟开采获得卸压开采过程中各区域的应力演化特征；
步骤S3、根据卸压开采覆岩损伤度将覆岩分为破碎煤岩体、裂隙煤岩体和弹性煤岩体，破碎煤岩体位于垮落带，裂隙煤岩体位于裂隙带，弹性煤岩体位于弯曲下沉带；
步骤S4、根据卸压开采过程中不同区域的应力演化特征设计渗透率测试的应力路径，获得破碎煤岩体加载阶段，裂隙煤岩体为卸载-加载阶段，弹性煤岩体是加载-卸载阶段；
步骤S5、采用渗透率测试系统分别进行不同损伤程度煤岩体的应力-渗透率测试，获得破碎煤岩体加载渗透率公式、裂隙煤岩体第一次加卸载渗透率公式和第二次加载渗透率公式、完整岩体的第一次加卸载渗透率公式；
步骤S6、构建卸压开采数值模型，在模拟过程中运行利用FISH语言更新渗透率；
步骤S7、首先判断卸压开采过程中覆岩所处分区，垮落带分区根据公式
根据直接顶单轴抗压强度对直接顶进行分类：Ⅰ类，单轴抗压强度大于40MPa；Ⅱ类，单轴抗压强度在20～40MPa；Ⅲ类，单轴抗压强度小于20MPa；Ⅰ类直接顶c1，c2和c3分别为2.1、16和2.5；Ⅱ类直接顶c1，c2和c3分别为4.7、19和2.2；Ⅲ类直接顶c1，c2和c3分别为6.2、32和1.5；其中，c1是硬岩垮落系数，c2是中硬岩垮落系数，c3是软岩垮落系数；
步骤S8、裂隙带分区则根据模拟过程中的弹塑性状态判断：弹性区属于弯曲下沉带，塑性区属于裂隙带。
步骤S9、判断为弯曲下沉带的块体进一步判断块体的膨胀变形，当膨胀变形大于等于0.3％时，则认为该区域处于离层裂隙带；
步骤S10、损伤程度判断完毕后，进一步判断块体所处的加卸载状态；
步骤S11、加卸载判断完毕后，进一步判断块体所经历的加卸载次数；
步骤S12、判断为垮落带的块体的水平垂直渗透率和垂直渗透率均采用破碎煤岩体渗透率模型计算；
步骤S13、判断为裂隙带的块体的水平垂直渗透率和垂直渗透率均采用裂隙煤岩体渗透率模型计算；
步骤S14、判断为离层裂隙带的块体的水平渗透率采用裂隙煤岩体渗透率模型计算，垂直渗透率采用弹性煤岩体渗透率计算模型；
步骤S15、膨胀变形小于0.3％的弹性块体的水平垂直渗透率和垂直渗透率采用弹性煤岩体渗透率计算模型。
本发明涉及深部高瓦斯煤层卸压开采效果评价方法中的渗透率计算方法，本发明按照覆岩“三带”损伤裂隙煤岩体特征将煤岩体进行分类并进行了循环加卸载应力-渗透率实验得到了不同加卸载阶段不同损伤程度煤样渗透率的拟合公式，从而建立了卸压开采覆岩渗透率更新模型。本发明对于卸压开采过程中覆岩渗透率演化规律实现了定量分析，能够准确高效的实现卸压开采覆岩渗透率的更新，对卸压开采效果评价及瓦斯抽采具有重要意义。
附图说明
图1为本发明的覆岩渗透率更新方法的工作原理流程图；
图2为本发明的覆岩损伤分区示意图；
图3为本发明的破碎岩样的拟合曲线；
图4为本发明的裂隙煤样的拟合曲线；
图5为本发明的完整煤样的拟合曲线。
附图标号说明：
1、采煤工作面；2、保护层；3、垮落带；4、裂隙带；5、弯曲下沉带；6、地表；7、直接顶；8、离层裂隙带；9、被保护层。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现说明本发明。
本发明通过实验室实测获得卸压开采矿井各岩层力学参数，运用相似模拟获得采空区垮落带、裂隙带和弯曲下沉带内煤岩层的应力路径。按照覆岩损伤程度将煤岩体分成完整煤岩体(弯曲下沉带以上)、裂隙煤岩体(裂隙带)和破碎煤岩体(垮落带)。根据实际采动应力路径实测不同损伤岩体的应力-渗透率演化规律，获得加卸载条件下不同损伤程度渗透率模型。基于不同加卸载阶段不同损伤程度渗透率模型利用FISH语言嵌入至FLAC3D数值模拟软件中对其进行二次开发。采用开发的数值模型进行卸压开始数值模拟，并实时判断岩体损伤程度和加卸载状态更新渗透率模拟卸压开采过程中渗透率与瓦斯压力的演化规律，进而评价被保护层的卸压效果及卸压范围。
如图1和图2所示，本发明提供一种考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法，具体如下：
a、实测获得卸压开采各岩层力学参数，包括保护层2和被保护层9，并进行相似模拟配比；
b、采用相似模拟开采获得卸压开采过程中各区域的应力演化特征；
c、根据卸压开采覆岩损伤度将覆岩分为破碎煤岩体、裂隙煤岩体和弹性煤岩体，破碎煤岩体位于垮落带3、裂隙煤岩体位于裂隙带4、弹性煤岩体位于弯曲下沉带5；
d、根据卸压开采过程中不同区域的应力演化特征设计渗透率测试的应力路径，获得破碎煤岩体加载阶段，裂隙煤岩体卸载-加载阶段，弹性煤岩体加载-卸载阶段；
e、采用渗透率测试系统分别进行不同损伤程度煤岩体的应力-渗透率测试，获得破碎煤岩体加载渗透率公式
裂隙煤岩体第一次加载渗透率公式
第一次卸载渗透率公式
和第二次加载渗透率公式
完整岩体的第一次加载渗透率公式
第一次卸载渗透率公式
公式中k1为渗透率，单位md，σ1为有效应力，单位MPa。
上述公式可以按照建立模型和拟合得到，例如具体过程如下：
基于火柴棍模型，建立了同样适用于本文不规则裂隙煤岩体以及孔隙煤岩体的渗透率模型。为了使得拟合出来的公式尽量覆盖整个有效应力范围，将模型中初始应力状态假设为0，该处的渗透率也由拟合曲线得出。3个模型具体拟合公式如下：
式中，kf0为应力状态为0时的渗透率；cf0是原始裂隙的压缩系数，αf是裂隙压缩系数随有效应力的改变比率，σ1为有效应力，单位MPa。
基于火柴棍模型图提出了基于有效应力的渗透率理论模型：
式中，kf为煤层节理裂隙的渗透率；cf为裂隙压缩系数；σ1为有效应力；kf0，σ10为初始状态的渗透率及有效应力，σ10单位为MPa。
裂隙压缩系数cf并不是一个定值，随着有效应力的变化而变化，一般利用数据拟合得出的平均裂隙压缩系数
来代替：
式中，cf0是原始裂隙的压缩系数，αf是裂隙压缩系数随有效应力的改变比率。将式(5-2)带入式(5-1)得到简化公式：
根据公式(1)结合实验室不同应力条件下渗透率的测试数据拟合获得不同加卸载不同损伤程度的渗透率公式，具体如图3、图4和图5所示。
f、构建卸压开采数值模型，在模拟过程中运行利用FISH语言更新渗透率；
g、首先判断卸压开采过程中覆岩所处分区，垮落带3分区根据公式
计算获得(该公式根据以往实测数据统计获得，为经验公式)。根据直接顶7单轴抗压强度对直接顶7进行分类：Ⅰ类(单轴抗压强度大于40MPa)、Ⅱ类(单轴抗压强度在20～40MPa)和Ⅲ类(单轴抗压强度小于20MPa)。Ⅰ类直接顶c1，c2和c3分别为2.1、16和2.5；Ⅱ类直接顶c1，c2和c3分别为4.7、19和2.2；Ⅲ类直接顶c1，c2和c3分别为6.2、32和1.5；其中，c1是硬岩垮落系数，c2是中硬岩垮落系数，c3是软岩垮落系数；
h、裂隙带4分区则根据模拟过程中的弹塑性状态判断：弹性区属于弯曲下沉带5，塑性区属于裂隙带4。
i、判断为弯曲下沉带5的块体进一步判断块体的膨胀变形，当膨胀变形大于等于0.3％时，则认为该区域处于离层裂隙带8；
j、损伤程度判断完毕后，进一步判断块体所处的加卸载状态；
k、加卸载判断完毕后，进一步判断块体所经历的加卸载次数；
l、判断为垮落带3的块体水平和垂直渗透率均采用破碎煤岩体渗透率模型计算；
m、判断为裂隙带4的块体水平和垂直渗透率均采用裂隙煤岩体渗透率模型计算；
n、判断为离层裂隙带8的块体水平渗透率采用裂隙煤岩体渗透率模型计算，垂直渗透率采用弹性煤岩体渗透率计算模型；
o、膨胀变形小于0.3％的弹性块体水平和垂直渗透率采用弹性煤岩体渗透率计算模型。
本发明按照覆岩“三带”损伤裂隙煤岩体特征将煤岩体进行分类并进行了循环加卸载应力-渗透率实验得到了不同加卸载阶段不同损伤程度煤样渗透率的拟合公式，从而建立了卸压开采覆岩渗透率更新模型，能够高效快速精确的更新渗透率。
以上所述仅为本发明示意性的具体实施方式，并非用以限定本发明的范围。为本发明的各组成部分在不冲突的条件下可以相互组合，任何本领域的技术人员，在不脱离本发明的构思和原则的前提下所作出的等同变化与修改，均应属于本发明保护的范围。
Claims (4)
1.一种考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法，其特征在于，所述考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法包括以下步骤：
步骤a、实测获得卸压开采各岩层力学参数，并进行相似模拟配比；
步骤b、采用相似模拟开采获得卸压开采过程中各区域的应力演化特征；
步骤c、根据卸压开采覆岩损伤度将覆岩分为破碎煤岩体、裂隙煤岩体和弹性煤岩体，破碎煤岩体位于垮落带，裂隙煤岩体位于裂隙带，弹性煤岩体位于弯曲下沉带；
步骤d、根据卸压开采过程中不同区域的应力演化特征设计渗透率测试的应力路径，获得破碎煤岩体加载阶段，其中裂隙煤岩体为卸载-加载阶段，弹性煤岩体为加载-卸载阶段；
步骤e、采用渗透率测试系统分别进行不同损伤程度煤岩体的应力-渗透率测试，获得破碎煤岩体加载渗透率公式、裂隙煤岩体第一次加卸载渗透率公式和第二次加载渗透率公式、完整岩体的第一次加卸载渗透率公式；
步骤f、构建卸压开采数值模型，在模拟过程中根据获得的不同应力状态、不同损伤程度煤岩体的渗透率公式更新渗透率；步骤f中，在模拟过程中利用FISH语言更新渗透率；
步骤f中，在模拟过程中利用FISH语言更新渗透率具体包括：
步骤f1、首先判断卸压开采过程中覆岩所处分区；
步骤f2、判断为弯曲下沉带的块体进一步判断块体的膨胀变形，当膨胀变形大于等于0.3％时，则认为该区域处于离层裂隙带；
步骤f3、损伤程度判断完毕后，进一步判断块体所处的加卸载状态；
步骤f4、加卸载判断完毕后，进一步判断块体所经历的加卸载次数；
步骤f5、判断为垮落带的块体的水平渗透率和垂直渗透率均采用破碎煤岩体渗透率模型计算；
步骤f6、判断为裂隙带的块体的水平渗透率和垂直渗透率均采用裂隙煤岩体渗透率模型计算；
步骤f7、判断为离层裂隙带的块体的水平渗透率采用裂隙煤岩体渗透率模型计算，垂直渗透率采用弹性煤岩体渗透率计算模型；
步骤f8、膨胀变形小于0.3％的弹性块体的水平渗透率和垂直渗透率采用弹性煤岩体渗透率计算模型。
2.如权利要求1所述的考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法，其特征在于，步骤f1中，垮落带分区根据采高和直接顶系数在模拟过程中直接设置；裂隙带分区则根据模拟过程中的弹塑性状态判断：其中，弹性区属于弯曲下沉带，塑性区属于裂隙带。
3.如权利要求1所述的考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法，其特征在于，步骤f1中，垮落带分区确定采用如下方法：
f11、根据直接顶单轴抗压强度对直接顶分成3类，分别是：Ⅰ类，单轴抗压强度大于40MPa；Ⅱ类，单轴抗压强度在20～40MPa；和Ⅲ类，单轴抗压强度小于20MPa；
f12、Ⅰ类直接顶c1，c2和c3分别为2.1、16和2.5
f13、Ⅱ类直接顶c1，c2和c3分别为4.7、19和2.2
f14、Ⅲ类直接顶c1，c2和c3分别为6.2、32和1.5
f15、根据公式计算获得垮落带范围高度H1。
4.如权利要求1所述的考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法，其特征在于，具体包括以下步骤：
步骤S1、实测获得卸压开采各岩层力学参数，并进行相似模拟配比；
步骤S2、采用相似模拟开采获得卸压开采过程中各区域的应力演化特征；
步骤S3、根据卸压开采覆岩损伤度将覆岩分为破碎煤岩体、裂隙煤岩体和弹性煤岩体，破碎煤岩体位于垮落带，裂隙煤岩体位于裂隙带，弹性煤岩体位于弯曲下沉带；
步骤S4、根据卸压开采过程中不同区域的应力演化特征设计渗透率测试的应力路径，获得破碎煤岩体加载阶段，裂隙煤岩体为卸载-加载阶段，弹性煤岩体是加载-卸载阶段；
步骤S5、采用渗透率测试系统分别进行不同损伤程度煤岩体的应力-渗透率测试，获得破碎煤岩体加载渗透率公式、裂隙煤岩体第一次加卸载渗透率公式和第二次加载渗透率公式、完整岩体的第一次加卸载渗透率公式；
步骤S6、构建卸压开采数值模型，在模拟过程中运行利用FISH语言更新渗透率；
步骤S7、首先判断卸压开采过程中覆岩所处分区，垮落带分区根据公式
计算获得；
根据直接顶单轴抗压强度对直接顶进行分类：Ⅰ类，单轴抗压强度大于40MPa；Ⅱ类，单轴抗压强度在20～40MPa；Ⅲ类，单轴抗压强度小于20MPa；Ⅰ类直接顶c1，c2和c3分别为2.1、16和2.5；Ⅱ类直接顶c1，c2和c3分别为4.7、19和2.2；Ⅲ类直接顶c1，c2和c3分别为6.2、32和1.5；
步骤S8、裂隙带分区则根据模拟过程中的弹塑性状态判断：弹性区属于弯曲下沉带，塑性区属于裂隙带；
步骤S9、判断为弯曲下沉带的块体进一步判断块体的膨胀变形，当膨胀变形大于等于0.3％时，则认为该区域处于离层裂隙带；
步骤S10、损伤程度判断完毕后，进一步判断块体所处的加卸载状态；
步骤S11、加卸载判断完毕后，进一步判断块体所经历的加卸载次数；
步骤S12、判断为垮落带的块体的水平垂直渗透率和垂直渗透率均采用破碎煤岩体渗透率模型计算；
步骤S13、判断为裂隙带的块体的水平垂直渗透率和垂直渗透率均采用裂隙煤岩体渗透率模型计算；
步骤S14、判断为离层裂隙带的块体的水平渗透率采用裂隙煤岩体渗透率模型计算，垂直渗透率采用弹性煤岩体渗透率计算模型；
步骤S15、膨胀变形小于0.3％的弹性块体的水平垂直渗透率和垂直渗透率采用弹性煤岩体渗透率计算模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110033385.1A CN112883542B (zh) | 2021-01-12 | 2021-01-12 | 考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110033385.1A CN112883542B (zh) | 2021-01-12 | 2021-01-12 | 考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112883542A CN112883542A (zh) | 2021-06-01 |
| CN112883542B true CN112883542B (zh) | 2023-05-05 |
Family
ID=76044588
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110033385.1A Active CN112883542B (zh) | 2021-01-12 | 2021-01-12 | 考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112883542B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114813371B (zh) * | 2022-04-15 | 2024-05-14 | 安徽理工大学 | 一种基于渗透-损伤关系的覆岩损伤分区及高度确定方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108344675A (zh) * | 2018-02-08 | 2018-07-31 | 四川大学 | 模拟保护层开采条件下煤体采动渗流力学规律的试验方法 |
| CN109655392A (zh) * | 2018-12-03 | 2019-04-19 | 中国矿业大学(北京) | 一种破碎煤岩样可视化伺服加载渗流实验测试方法 |
| CN110221036A (zh) * | 2018-03-01 | 2019-09-10 | 中国矿业大学 | 带渗流装置的保水开采“声发射-红外辐射”实验系统 |
| CN110501232A (zh) * | 2019-07-04 | 2019-11-26 | 同济大学 | 基于双轴流变仪实现可视化的真三轴水力耦合试验装置 |
-
2021
- 2021-01-12 CN CN202110033385.1A patent/CN112883542B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108344675A (zh) * | 2018-02-08 | 2018-07-31 | 四川大学 | 模拟保护层开采条件下煤体采动渗流力学规律的试验方法 |
| CN110221036A (zh) * | 2018-03-01 | 2019-09-10 | 中国矿业大学 | 带渗流装置的保水开采“声发射-红外辐射”实验系统 |
| CN109655392A (zh) * | 2018-12-03 | 2019-04-19 | 中国矿业大学(北京) | 一种破碎煤岩样可视化伺服加载渗流实验测试方法 |
| CN110501232A (zh) * | 2019-07-04 | 2019-11-26 | 同济大学 | 基于双轴流变仪实现可视化的真三轴水力耦合试验装置 |
Non-Patent Citations (2)
| Title |
|---|
| 尹光志 等.真三轴加卸载应力路径对原煤力学特性及渗透率影响.《煤炭学报》.2018,第43卷(第1期),131-136. * |
| 张村.高瓦斯煤层群应力—裂隙—渗流耦合作用机理及其对卸压抽采的影响.《中国博士学位论文全文数据库 工程科技Ⅰ辑》.2018,(第1期),B021-39. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112883542A (zh) | 2021-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhou et al. | Damage analysis of rock mass coupling joints, water and microseismicity | |
| CN113820750B (zh) | 基于弹塑性力学定量预测泥岩构造裂缝的方法 | |
| CN103913772B (zh) | 基于储层地质力学参数的微地震事件正演模拟方法 | |
| EP3789941A1 (en) | Method and apparatus for predicting production of oil and gas obtained from shale oil in-situ exploitation | |
| Tulu et al. | Verification of a calibrated longwall model with field measurements | |
| Khanal et al. | Prefeasibility study—Geotechnical studies for introducing longwall top coal caving in Indian mines | |
| CN105589999B (zh) | 应用于地下工程围岩注浆方案的确定方法 | |
| CN107991188B (zh) | 一种基于岩心残余应力水平预测水力裂缝复杂性的方法 | |
| Song et al. | Experimental investigation and numerical simulation of surrounding rock creep for deep mining tunnels | |
| CN115267905B (zh) | 一种复杂构造区钻井工程中裂缝性井漏预测方法 | |
| CN112883542B (zh) | 考虑采动应力路径与围岩损伤程度的覆岩渗透率更新方法 | |
| CN116127752B (zh) | 一种冲击地压的预测方法及系统 | |
| Elahi* et al. | Characterization of fracture length and conductivity from tracer test and production data with Ensemble Kalman filter | |
| CN110705168A (zh) | 构造应力场的模拟方法 | |
| CN109064016B (zh) | 一种低渗煤层水力压裂增透效果评价方法 | |
| CN109268003A (zh) | 一种地质工程一体化框架下的致密油增产数学建模方法 | |
| CN115270533A (zh) | 重复压裂设计方法、装置、存储介质及电子设备 | |
| CN102609780A (zh) | 矿井瓦斯涌出量预测新方法 | |
| CN106291715A (zh) | 一种基于断裂自相似理论的低级序断层发育规律预测方法 | |
| Bakhtiyari et al. | Determination of shear strength parameters of rock mass using back analysis methods and comparison of results with empirical methods | |
| Ou et al. | Three-dimensional discrete network modeling of structural fractures based on the geometric restoration of structure surface: Methodology and its application | |
| CN114692456A (zh) | 一种诱发矿井强动力灾害地应力临界判据的确定方法 | |
| Souley et al. | Damage around a keyed URL excavation: change in permeability induced by microcracks growth | |
| Li et al. | Calculating the potential for coal pillar bumps using a local mine stiffness criterion | |
| Kaifang et al. | Using field measurements to quantify the post-failure behavior of coal |
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 |