WO2016074510A1 - 浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法 - Google Patents

浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法 Download PDF

Info

Publication number
WO2016074510A1
WO2016074510A1 PCT/CN2015/086602 CN2015086602W WO2016074510A1 WO 2016074510 A1 WO2016074510 A1 WO 2016074510A1 CN 2015086602 W CN2015086602 W CN 2015086602W WO 2016074510 A1 WO2016074510 A1 WO 2016074510A1
Authority
WO
WIPO (PCT)
Prior art keywords
model
air leakage
fracture
coal
rock
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.)
Ceased
Application number
PCT/CN2015/086602
Other languages
English (en)
French (fr)
Inventor
秦波涛
王奇奇
申宏敏
马立强
鲁义
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
China University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology CUMT
China University of Mining and Technology Beijing CUMTB
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT, China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology CUMT
Priority to AU2015345707A priority Critical patent/AU2015345707B2/en
Publication of WO2016074510A1 publication Critical patent/WO2016074510A1/zh
Priority to ZA2016/06183A priority patent/ZA201606183B/en
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]

Definitions

  • the invention relates to a method for judging the distribution of surface through cracks and the characteristics of air leakage in shallow buried coal seam mining, and belongs to the experimental research on fractured rock mass in the field of underground geotechnical engineering and the method for judging the characteristics of surface fissure leakage.
  • the methods for determining the surface fissure distribution and air leakage characteristics of shallow buried coal seams are mainly tracer gas method and numerical simulation method.
  • the tracer gas method is to release the SF 6 tracer gas at the air leakage source, collect the gas sample at the air leakage, and qualitatively determine the air leakage channel by analyzing the concentration of the gas sample. Due to the limitations of the site environment and measurement methods, the crack distribution and air leakage characteristics of the entire goaf can only be reflected by monitoring some points, and the judgment results are susceptible to factors such as the release amount.
  • the equivalent continuous medium model is generally used to make the fissure and the surrounding rock mass equivalent to a continuous medium with a certain permeability tensor, which is solved by the porous medium theory.
  • it ignores the influence of the longitudinal cracks of the shallow buried coal seams on the surface leakage, and the simulation results often have large deviations from the actual conditions when dealing with such large-scale fractures.
  • the present invention introduces the development of fractures into a numerical model by combining similar material simulation experiments with numerical simulations, thereby overcoming the difficulty in detecting fractures in the rock and over-exaggerating the numerical model. Simplification and other problems; at the same time, the correction of the model can effectively improve the accuracy and reliability of the model, and provide a useful reference for the determination of surface through-fracture distribution and air leakage characteristics after mining in shallow buried mining areas.
  • a method for judging the distribution of surface through cracks and the characteristics of air leakage in shallow buried coal seam mining comprising the following steps:
  • the experimental rock formation model is laid in order according to the stratigraphic relationship and inclination angle of the original rock, and the model is placed, and the resistance strain gauge is arranged in the adjacent rock layer;
  • the resistance strain gauge is disposed between two adjacent rock layers, and the resistance strain gauges in the same horizontal detection plane are arranged in a mesh shape to collect detection data; for example, designing the same level
  • the strain gauges in the detection plane are arranged in a grid pattern of a rectangular array, and the horizontal distance between two adjacent strain gauges on the same lateral or longitudinal straight line is generally designed to be 30 cm.
  • the step (3) it is determined whether the difference between the model strength and the original rock strength is within a threshold range, and the specific method is: before the laying of the model, the mechanical property test is performed to determine that the simulated material reaches a difference between the mechanical properties of the original rock and the threshold value.
  • the water content w 0 in the range; after the model is laid and allowed to stand for a period of time, the water content w of the model material is measured, and when w w 0 , the difference between the model strength and the original rock strength is considered to be within the threshold range.
  • the method for processing the photographed fissure development photograph into a vector graph is: using computer graphics processing technology, including image filtering, sharpening enhancement, image segmentation, noise filtering, and detection refinement.
  • image filtering sharpening enhancement
  • image segmentation image segmentation
  • noise filtering and detection refinement.
  • vectorized fracture data is generated, and the vectorized fracture data is used as a vector graphic.
  • the material properties include a fluid density ⁇ , a hydrodynamic viscosity ⁇ , a coal rock permeability k around the fracture, and a coal rock porosity ⁇ ;
  • the boundary condition is specifically set as: an upper crack inlet
  • the pressure p 0 is set to atmospheric pressure
  • the lower fracture outlet pressure is set to the goaf side pressure
  • the left and right boundaries are set to no flow boundary.
  • the method for solving the permeability k and the porosity ⁇ of the coal rock surrounding the crack is: taking four displacement monitoring points adjacent to each other on the model plane to form a quadrilateral ABCD, the coal seam Mining, when the overburden collapses, the area of the quadrilateral ABCD changes from S to S':
  • the geometric model is described as follows:
  • ⁇ represents fluid density
  • u represents fluid velocity
  • ⁇ represents hydrodynamic viscosity
  • p represents unit fluid pressure difference
  • F unit fluid volume force
  • ⁇ denotes hydrodynamic viscosity
  • k is the permeability of the coal rock mass
  • q fluid flow rate is the unit fluid pressure difference
  • Z is the height change amount
  • the method for judging the distribution of surface through cracks and the characteristics of air leakage in shallow buried coal seams provided by the present invention has the following advantages:
  • Figure 1 is a flow chart of the method of the present invention
  • Figure 2 is a schematic diagram of calculation of the coefficient of expansion of the goaf
  • Figure 3 is a distribution diagram of the resistance strain gauge and the displacement monitoring point
  • Figure 4 is a fracture distribution map after simulated coal seam mining
  • Figure 5 is a vector graphic
  • Figure 6 is a distribution diagram of the fracture air leakage velocity.
  • FIG. 1 is a flow chart showing an implementation method for determining the surface through-fracture distribution and air leakage characteristics of a shallow buried coal seam, and the present invention will be further described below with reference to examples.
  • a coal mine in Shendong mining area is a shallow buried mine.
  • the ground fissures are developed after coal seam mining, and the air leakage in the goaf is serious, causing spontaneous combustion of coal.
  • the air leakage channel sealing scheme is given by the method of the present invention. The specific steps are as follows:
  • the experimental rock formation model is laid down from the bottom to the top according to the stratigraphic relationship and inclination angle of the original rock, and the strain gauge is placed in the adjacent rock layer;
  • the arrangement of the strain gauges is as follows: the strain gauges are arranged in a mesh shape in the same detection plane, and the spacing between adjacent strain gauges in the same detection plane is 30 cm, and the arrangement of the strain gauges is as shown in FIG.
  • Each excavation uses a static strain measurement processor to record relevant data through a computer. After the entire model is excavated, when the data recorded by the computer does not change, the model reaches stress balance, and the professional camera is used after the model is excavated. A photograph of the development of the fracture of the model is shown in Figure 4.
  • the method for solving the permeability k and porosity ⁇ of the coal and rock mass around the fissure is: taking four displacement monitoring points adjacent to each other on the model plane to form a quadrilateral ABCD, coal seam mining, when the overburden layer collapses, The area of the quadrilateral ABCD changes from S to S':
  • the original porosity and permeability of the coal and rock mass are brought into the above formula to obtain the porosity and permeability of the coal and rock mass around the fracture.
  • ⁇ represents the fluid density
  • u represents the fluid velocity
  • ⁇ represents the hydrodynamic viscosity
  • p represents the unit fluid pressure difference
  • F unit fluid volume force
  • coal rock mass around the fracture zone is treated as a porous medium, which belongs to seepage, and is described by Darcy's law:
  • ⁇ denotes hydrodynamic viscosity
  • k is the permeability of the coal rock mass
  • q fluid flow rate is the unit fluid pressure difference
  • Z is the height change amount
  • the corresponding point is taken from the numerical model.
  • the simulation data is in good agreement with the actual data, and the model does not need to be corrected.
  • the simulation results can reflect the crack distribution and air leakage characteristics of shallow coal seams, and are used to guide the sealing of the air leakage passages in the goaf, thus preventing the coal spontaneous combustion in the mine.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Examining Or Testing Airtightness (AREA)
  • Excavating Of Shafts Or Tunnels (AREA)
  • Lining And Supports For Tunnels (AREA)

Abstract

一种浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，包括如下步骤：根据矿井实际地层资料确定相似实验材料配比；根据几何相似和动力相似铺设相似材料模型；材料强度与原岩强度相似后进行煤层开挖，模型开挖稳定后拍摄裂隙发育照片；使用图像处理软件对照片进行灰度化和矢量化处理；将矢量化的裂隙图像导入到数值模拟模拟软件进行计算；将得到的结果与实测数据进行对比分析，通过不断修正，最终获得比较精确的数值模型。该方法可揭示浅埋煤层条件下覆岩裂隙分布及漏风特征，为漏风通道的封堵等现场情况提供一种有益的参考。

Description

浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法 技术领域
本发明涉及一种浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，属于对地下岩土工程领域的裂隙岩体的试验研究及地表裂隙漏风特征的判定方法研究。
背景技术
我国煤炭开采战略西移，西部矿区多为浅埋藏煤层，浅埋藏煤层开采地表漏风严重，易引发采空区煤炭自燃。煤炭自燃不仅会影响矿井的正常生产，还可能引发重特大的火灾或瓦斯爆炸事故。由于西部煤层埋藏浅，在采动应力的作用下形成了大量贯通地表的裂隙，这些裂隙构成了煤自燃漏风供氧的主要通道。在利用黄泥(粉煤灰)浆、砂浆、三相泡沫等防灭火材料对采空区自燃危险区域及漏风点进行防灭火及降温封堵时，因传统的技术手段很难检测到裂隙的分布位置，导致防灭火材料不能及时有效地输送到漏风点附近完成对漏风通道的封堵。所以，研究浅埋藏煤层开采裂隙的分布和漏风特征对封堵漏风通道，防治煤炭自燃有着重要意义。
目前，判定浅埋藏煤层地表裂隙分布和漏风特征的方法主要有示踪气体法和数值模拟法。示踪气体法是在漏风源释放SF6示踪气体，在漏风汇处采集气样，通过对气样浓度的分析来定性判定漏风通道。由于受现场环境和测量方法的限制，只能通过对一些点的监测来反映整个采空区的裂隙分布和漏风特征，且判定结果易受释放量等因素的影响。在对浅埋煤层裂隙漏风进行数值模拟时，一般使用等效连续介质模型，将裂隙及周围岩体等效为具有一定渗透张量的连续介质，利用多孔介质理论进行求解。但它忽略了浅埋煤层贯通地表的纵向裂隙对漏风的影响，在处理这种大规模裂隙时模拟结果与实际情况往往存在较大偏差。
发明内容
发明目的：针对现有数值模拟方法的不足，本发明通过将相似材料模拟实验与数值模拟相结合的方法，将裂隙的发育状况引入数值模型中，克服了岩体内裂隙检测困难、数值模型过度简化等问题；同时，对模型的修正可以有效的提高模型的精度和可靠性，为浅埋藏矿区开采后地表贯通裂隙分布和漏风特征判定提供一种有益的参考。
技术方案：为实现上述目的，本发明采用的技术方案为：
一种浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，包括如下步骤：
(1)确定模型与原岩的比例，根据矿区煤层埋藏的岩层岩性、厚度和物理力学参数，计算模型中模拟各层岩层时不同材料的配比及用量；
(2)根据得到的材料的配比及用量，按照原岩的岩层层位关系及倾角按顺序铺设实验岩层形成模型并静置，在相邻岩层中布置电阻应变片；
(3)当模型强度与原岩强度差在阈值范围内时，模拟现场实际对原型的开采条件，准备对模型中的煤层进行开挖；
(4)根据对原型进行实际开挖时的推进速度和每次开挖的长度，设置对模型进行开挖的推进速度和每次开挖的长度，并且每次开挖结束后，放置40～80min再继续开挖；
(5)在对模型开挖的进行过程中，记录电阻应变片的检测数据，当各个电阻应变片的数据均不再变动或变动幅度均在阈值范围内时，模型达到应力平衡，模型开挖完成后使用相机拍摄模型应力平衡后裂隙发育的照片；
(6)将拍摄所得的裂隙发育照片处理为矢量图形；
(7)将矢量图形导入COMSOL数值模拟软件并设为初始几何模型，调整几何模型大小，设定几何模型材料属性、边界条件；
(8)对设定好的几何模型进行网格剖分后求解计算，获得裂隙漏风风速和压力分布；
(9)将获得的裂隙漏风风速及压力分布与针对原型进行现场实测的各点的漏风数据进行对比分析，通过不断地调整几何模型的设计参数，从而获得与现场实测相吻合的裂隙漏风风速和压力分布规律，为封堵漏风通道提供参考。
具体的，所述步骤(2)中，电阻应变片布置在相邻两个岩层之间，在同一个水平检测平面内的电阻应变片呈网状布置，以采集检测数据；比如设计同一个水平检测平面内的电阻应变片呈矩形阵列的网格状分布，一般设计同一横向或纵向直线上相邻两电阻应变片之间的水平距离为30cm。
具体的，所述步骤(3)中，判断模型强度与原岩强度差是否在阈值范围内，具体方法为：在铺设模型前，通过力学性能实验确定模拟材料达到与原岩力学性能差在阈值范围内时的含水量w0；在模型铺设完成并静置一段时间后，测量模型材料的含水量w，当w＝w0时，即可认为模型强度与原岩强度差在阈值范围内。
更为具体的，所述步骤(3)中，确定材料含水量的方法是称重法，具体为：取一定量的材料作为试样，使用0.1g精度的天平称取试样的重量，记作试样的湿重m，在105℃的烘箱内将试样烘烤至恒重，再次使用0.1g精度的天平称取试样的重量，记作试样的湿重ms，计算含水量w＝ms/m。
具体的，所述步骤(6)中，将拍摄所得的裂隙发育照片处理为矢量图形的方法为：利用计算机图形处理技术，通过包括图像滤波、锐化增强、图像分割、噪音过滤和检测细化在内的处理后，生成矢量化的裂隙数据，将矢量化的裂隙数据作为矢量图形。
具体的，所述步骤(7)中，材料属性包括流体密度ρ、流体动力粘度μ、裂隙周围煤岩体渗透率k和煤岩体孔隙率ε；边界条件的设定具体为：上部裂隙入口压力p0设为大气压力，下部裂隙出口压力设为采空区侧压力，左右边界设为无流动边界。
更为具体的，所述步骤(7)中，裂隙周围煤岩体渗透率k和孔隙率ε的求解方法为：在模型平面上取上下相邻的四个位移监测点构成一个四边形ABCD，煤层开采，当上覆岩层发生塌陷后，四边形ABCD的面积由S变为S'：
计算煤岩体碎胀系数为：Kp＝S'/S；
根据煤岩体碎胀系数计算孔隙率为：
Figure PCTCN2015086602-appb-000001
煤岩体渗透率k和孔隙率ε满足：
Figure PCTCN2015086602-appb-000002
其中，d为破碎煤岩体粒径，C为系数，一般取C＝172.8。
具体的，所述步骤(7)中，几何模型按照如下方式进行描述：
1)裂隙区域内部流体自由流动，采用Navier-Stokes方程描述：
Figure PCTCN2015086602-appb-000003
其中，ρ表示流体密度，u表示流体速度，μ表示流体动力粘度，p表示单位流体压力差，F单位流体体积力；
2)裂隙区域周围煤岩体处理为多孔介质，属于渗流，采用Darcy定律描述：
Figure PCTCN2015086602-appb-000004
其中，μ表示流体动力粘度，k为煤岩体的渗透率，q流体流量，p为单位流体压力差，Z为高度改变量。
有益效果：本发明提供的浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，具有如下优点：
1、提出了通过测定含水量来确定相似材料强度的方法，具有简单、方便等优点；
2、本发明通过相似材料模拟实验可以得到模拟煤层开采后模型上覆岩层裂隙分布特征图，可以直观观察煤层开采后不同区域裂隙的分布情况，可以用于覆岩裂隙漏风特征的数值模拟分析；
3、通过将相似材料模拟实验与数值模拟相结合的方法，可以降低将裂隙处理为等效 连续介质或裂隙网络模型时产生的误差，使得数值模拟结果与实际情况更加吻合。
附图说明
图1为本发明方法流程图；
图2为采空区碎胀系数计算示意图；
图3为电阻应变片及位移监测点分布图；
图4为模拟煤层开采后裂隙分布图；
图5为矢量图形；
图6为裂隙漏风速度分布图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1所示为一种浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法实施流程图，下面结合实例对本发明做出进一步的说明。
神东矿区某煤矿属浅埋矿井，煤层开采后地面裂隙发育，采空区漏风严重，引起煤炭自燃，利用本发明方法给出漏风通道封堵方案，具体步骤如下：
(1)确定模型与原岩的比例为1:100，模型长为2.5m，宽为0.5m，根据矿方提供的岩层岩性、厚度和物理力学参数，计算模型中模拟各层岩层时不同材料的配比及用量，如表1所示：
表1 模型相似材料模拟配比(1:100)
Figure PCTCN2015086602-appb-000005
(2)根据得到的材料的配比及用量，按照原岩的岩层层位关系及倾角按顺序从下至上分层铺设实验岩层形成模型并静置，在相邻岩层中布置电阻应变片；电阻应变片的布置方法为：在同一检测平面内电阻应变片呈网状布置，同一检测平面内相邻两电阻应变片之间的间距为30cm，电阻应变片的布置情况如图3所示。
(3)模型铺设完成后静置10天，然后从模型上部边缘取一小块试样，用称重法测量含水量，与和原岩强度相似时同层试样的含水量进行对比发现二者含水量差别小于5％，在差别范围内，认为模型材料强度与原岩强度相似，可以进行开挖。
(4)由矿井的实际推进速度计算每次模拟开挖的长度为10cm，开挖后放置1小时，再次进行开挖。
(5)每次开挖均使用静态应变测量处理仪通过计算机记录相关数据，整个模型开挖完成后，当计算机记录的数据不再变动时，模型达到应力平衡，模型开挖完成后用专业相机拍摄模型裂隙发育的照片，如图4所示。
(6)利用图形处理软件将拍摄的裂隙发育照片处理为矢量图形，如图5所示。
(7)将步骤(6)得到的矢量图形导入COMSOL数值模拟软件并设为初始几何图形。调整几何模型的大小，设定流体密度ρ＝1.29kg/m3、流体动力粘度μ＝17.9×10-6Pa·s、煤岩体渗透率k和煤岩体孔隙率ε；边界条件的设定具体为：上部裂隙入口压力p0＝1atm，下部裂隙出口压力设为采空区侧压力101.12kpa，左右边界设为无流动边界。
(8)裂隙周围煤岩体渗透率k和孔隙率ε的求解方法为：在模型平面上取上下相邻的四个位移监测点构成一个四边形ABCD，煤层开采，当上覆岩层发生塌陷后，四边形ABCD的面积由S变为S'：
计算煤岩体碎胀系数为：Kp＝S'/S；
根据煤岩体碎胀系数计算孔隙率为：
Figure PCTCN2015086602-appb-000006
煤岩体渗透率k和孔隙率ε满足：
Figure PCTCN2015086602-appb-000007
其中，d为破碎煤岩体粒径，C为系数，一般取C＝172.8。
根据矿方提供的资料将煤岩体原始孔隙率和渗透率带入上式求得裂隙周围煤岩体孔隙率和渗透率。
(9)由于裂隙区域内流体自由流动，采用Navier-Stokes(纳维-斯托克斯)方程描述：
Figure PCTCN2015086602-appb-000008
其中，ρ表示流体密度，u表示流体速度，μ表示流体动力粘度，p表示单位流体压力差，F单位流体体积力。
将裂隙区域周围煤岩体处理为多孔介质，属于渗流，采用Darcy定律描述：
Figure PCTCN2015086602-appb-000009
其中，μ表示流体动力粘度，k为煤岩体的渗透率，q流体流量，p为单位流体压力差，Z为高度改变量。
(10)对模型进行网格剖分，求解系统方程组，获得裂隙漏风风速和压力分布图，如图6所示。
(11)根据现场检测点的位置，从数值模型中取出与之对应的点。通过对比发现现场监测点的风速、风压与数值模型中对应点的风速、风压，发现二者差值小于20％，模拟数据与实际数据吻合良好，不需要对模型进行修正。模拟结果可以反映浅埋煤层裂隙分布和漏风特征，用于指导对采空区漏风通道的封堵，从而防止矿井的煤自燃。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (8)

  1. 一种浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，其特征在于：包括如下步骤：
    (1)确定模型与原岩的比例，根据矿区煤层埋藏的岩层岩性、厚度和物理力学参数，计算模型中模拟各层岩层时不同材料的配比及用量；
    (2)根据得到的材料的配比及用量，按照原岩的岩层层位关系及倾角按顺序铺设实验岩层形成模型并静置，在相邻岩层中布置电阻应变片；
    (3)当模型强度与原岩强度差在阈值范围内时，模拟现场实际对原型的开采条件，准备对模型中的煤层进行开挖；
    (4)根据对原型进行实际开挖时的推进速度和每次开挖的长度，设置对模型进行开挖的推进速度和每次开挖的长度，并且每次开挖结束后，放置40～80min再继续开挖；
    (5)在对模型开挖的进行过程中，记录电阻应变片的检测数据，当各个电阻应变片的数据均不再变动或变动幅度均在阈值范围内时，模型达到应力平衡，模型开挖完成后使用相机拍摄模型应力平衡后裂隙发育的照片；
    (6)将拍摄所得的裂隙发育照片处理为矢量图形；
    (7)将矢量图形导入COMSOL数值模拟软件并设为初始几何模型，调整几何模型大小，设定几何模型材料属性、边界条件；
    (8)对设定好的几何模型进行网格剖分后求解计算，获得裂隙漏风风速和压力分布；
    (9)将获得的裂隙漏风风速及压力分布与针对原型进行现场实测的各点的漏风数据进行对比分析，通过不断地调整几何模型的设计参数，从而获得与现场实测相吻合的裂隙漏风风速和压力分布规律，为封堵漏风通道提供参考。
  2. 根据权利要求1所述的浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，其特征在于：所述步骤(2)中，电阻应变片布置在相邻两个岩层之间，在同一个水平检测平面内的电阻应变片呈网状布置，相邻两电阻应变片之间的水平距离为30cm，以采集检测数据。
  3. 根据权利要求1所述的浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，其特征在于：所述步骤(3)中，判断模型强度与原岩强度差是否在阈值范围内，具体方法为：在铺设模型前，通过力学性能实验确定模拟材料达到与原岩力学性能差在阈值范围内时的含水量w0；在模型铺设完成并静置一段时间后，测量模型材料的含水量w，当w＝w0 时，即可认为模型强度与原岩强度差在阈值范围内。
  4. 根据权利要求3所述的浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，其特征在于：所述步骤(3)中，确定材料含水量的方法是称重法，具体为：取一定量的材料作为试样，使用0.1g精度的天平称取试样的重量，记作试样的湿重m，在105℃的烘箱内将试样烘烤至恒重，再次使用0.1g精度的天平称取试样的重量，记作试样的湿重ms，计算含水量w＝ms/m。
  5. 根据权利要求1所述的浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，其特征在于：所述步骤(6)中，将拍摄所得的裂隙发育照片处理为矢量图形的方法为：利用计算机图形处理技术，通过包括图像滤波、锐化增强、图像分割、噪音过滤和检测细化在内的处理后，生成矢量化的裂隙数据，将矢量化的裂隙数据作为矢量图形。
  6. 根据权利要求1所述的浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，其特征在于：所述步骤(7)中，材料属性包括流体密度ρ、流体动力粘度μ、裂隙周围煤岩体渗透率k和煤岩体孔隙率ε；边界条件的设定具体为：上部裂隙入口压力p0设为大气压力，下部裂隙出口压力设为采空区侧压力，左右边界设为无流动边界。
  7. 根据权利要求6所述的浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，其特征在于：所述步骤(7)中，裂隙周围煤岩体渗透率k和孔隙率ε的求解方法为：在模型平面上取上下相邻的四个位移监测点构成一个四边形ABCD，煤层开采，当上覆岩层发生塌陷后，四边形ABCD的面积由S变为S'：
    计算煤岩体碎胀系数为：Kp＝S'/S；
    根据煤岩体碎胀系数计算孔隙率为：
    Figure PCTCN2015086602-appb-100001
    煤岩体渗透率k和孔隙率ε满足：
    Figure PCTCN2015086602-appb-100002
    其中，d为破碎煤岩体粒径，C为系数,一般取C＝172.8。
  8. 根据权利要求7所述的浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法，其特征在于：所述步骤(7)中，几何模型按照如下方式进行描述：
    1)裂隙区域内部流体自由流动，采用Navier-Stokes方程描述：
    Figure PCTCN2015086602-appb-100003
    其中，ρ表示流体密度，u表示流体速度，μ表示流体动力粘度，p表示单位流体压力差，F单位流体体积力；
    2)裂隙区域周围煤岩体处理为多孔介质，属于渗流，采用Darcy定律描述：
    Figure PCTCN2015086602-appb-100004
    其中，μ表示流体动力粘度，k为煤岩体的渗透率，q流体流量，p为单位流体压力差，Z为高度改变量。
PCT/CN2015/086602 2014-11-11 2015-08-11 浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法 Ceased WO2016074510A1 (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
AU2015345707A AU2015345707B2 (en) 2014-11-11 2015-08-11 Method for determining earth surface interpenetrated crack distribution and air leakage characteristics in shallow burial coal mining
ZA2016/06183A ZA201606183B (en) 2014-11-11 2016-09-06 Method for determining earth surface interpenetrated crack distribution and air leakage characteristics in shallow burial coal mining

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201410633316.4A CN104462654B (zh) 2014-11-11 2014-11-11 浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法
CN201410633316.4 2014-11-11

Publications (1)

Publication Number Publication Date
WO2016074510A1 true WO2016074510A1 (zh) 2016-05-19

Family

ID=52908685

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2015/086602 Ceased WO2016074510A1 (zh) 2014-11-11 2015-08-11 浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法

Country Status (4)

Country Link
CN (1) CN104462654B (zh)
AU (1) AU2015345707B2 (zh)
WO (1) WO2016074510A1 (zh)
ZA (1) ZA201606183B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106869895A (zh) * 2017-04-07 2017-06-20 中国石油大学(华东) 一种室内煤岩水力压裂水泥环破坏的物模制作方法
CN108932365A (zh) * 2017-05-23 2018-12-04 通用电气公司 涡轮机润滑油分析器系统、计算机程序产品和相关方法
CN111027787A (zh) * 2020-01-06 2020-04-17 中国矿业大学(北京) 一种基于断层滑移失稳诱发煤与瓦斯突出的实时预测方法
CN112730730A (zh) * 2020-12-23 2021-04-30 长安大学 一种高速公路下伏燃烧煤层的模型试验装置及试验方法
CN113030430A (zh) * 2021-03-17 2021-06-25 中煤科工集团重庆研究院有限公司 煤矿采动稳定区储层内的游离煤层气浓度测算方法
CN113567061A (zh) * 2021-07-23 2021-10-29 安徽理工大学 一种近距离煤层开采采空区群漏风检测方法
CN114397421A (zh) * 2021-11-29 2022-04-26 中煤科工开采研究院有限公司 采动煤层底板破坏深度监测方法
CN115824853A (zh) * 2022-11-18 2023-03-21 国能神东煤炭集团有限责任公司 缩短地下水库与煤炭开采空间安全距离的方法
CN117072202A (zh) * 2023-08-22 2023-11-17 中国地质大学（北京） 一种黄土丘陵采煤塌陷区防渗堵漏处理方法

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462654B (zh) * 2014-11-11 2017-06-13 中国矿业大学 浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法
CN105093352A (zh) * 2015-08-11 2015-11-25 武汉迈格睿地质环境科技有限公司 一种野外现场测量统计岩体体裂隙率的方法
CN105787220B (zh) * 2016-04-22 2018-10-12 山东科技大学 一种煤层高压注水致裂-渗流数值模拟方法
CN106339546B (zh) * 2016-08-26 2019-03-19 西安科技大学 一种极近距离煤层采空区下工作面支架载荷确定方法及装置
CN106593523B (zh) * 2017-01-23 2018-08-17 中国矿业大学 一种浅埋近距煤层开采过程中采空区漏风控制方法
CN107103143A (zh) * 2017-05-04 2017-08-29 中国矿业大学 采动影响条件下工作面上覆岩体渗透率变化的预测方法
CN108959691B (zh) * 2018-04-28 2022-03-08 湘潭大学 高瓦斯易自燃采空区耦合灾害之动态隔离参数的确定方法
CN108763650B (zh) * 2018-04-28 2022-03-11 湘潭大学 一种覆岩采动裂隙网络模型构建方法
CN109211666B (zh) * 2018-08-31 2019-12-03 山东科技大学 基于ct扫描的预测应力加载条件下煤体渗透率的方法
CN110570401B (zh) * 2019-08-19 2020-04-21 成都理工大学 一种基于图像识别的斜坡碎裂岩体信息获取方法
CN110566276B (zh) * 2019-09-05 2020-09-25 西安科技大学 一种利用采空区漏风判别煤自燃危险区域的方法
CN111504563B (zh) * 2020-04-28 2021-11-09 安徽理工大学 一种抽出式通风浅埋煤层矿井的地表漏风检测方法
CN112378738A (zh) * 2020-11-13 2021-02-19 重庆大学 煤矿采动稳定区煤层气储层裂隙表征方法
CN112364513A (zh) * 2020-11-13 2021-02-12 重庆大学 煤矿采动稳定区煤层气储层范围的界定方法
CN112364519B (zh) * 2020-11-19 2023-08-25 山西工程技术学院 一种用于抽采上隅角瓦斯的大直径钻孔参数确定方法
CN112685926B (zh) * 2020-12-02 2023-12-26 西安科技大学 一种浅埋煤层防治漏风工程的施工周期判别方法
CN114813822B (zh) * 2021-01-18 2025-03-21 神华神东煤炭集团有限责任公司 一种地裂缝漏风监测模拟实验方法及实验装置
CN113006869B (zh) * 2021-03-12 2023-07-25 神华神东煤炭集团有限责任公司 一种采动地裂缝预测方法
CN114419983B (zh) * 2022-01-26 2022-10-25 中南大学 倾斜煤层开采扰动岩层空隙率的随机分布特征的模拟方法
CN116147711B (zh) * 2023-04-17 2023-07-11 山西潞安环保能源开发股份有限公司 一种煤矿上覆岩层裂隙时空演化规律测试装置及方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102279131A (zh) * 2011-07-18 2011-12-14 中国石油大学(北京) 一种煤层水力压裂模拟实验方法
CN202814926U (zh) * 2012-08-31 2013-03-20 中国矿业大学(北京) 一种充填开采相似模拟实验微观裂隙测试系统
CN103674658A (zh) * 2013-12-24 2014-03-26 山东科技大学 一种随机裂隙试验模型的制备方法
CN103698171A (zh) * 2013-11-12 2014-04-02 山东科技大学 底板岩层固流耦合材料的裂隙演化试验方法
CN104462654A (zh) * 2014-11-11 2015-03-25 中国矿业大学 浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101070759B (zh) * 2007-05-04 2010-05-19 中国矿业大学 一种薄基岩浅埋煤层保水开采适用条件分类方法
CN102777201B (zh) * 2012-07-26 2014-06-18 山东科技大学 基于正压通风系统的火区下近距离煤层开采通风方法
CN103806943B (zh) * 2013-12-13 2016-06-22 中国神华能源股份有限公司 用于监测煤层采空区漏风的监测装置及监测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102279131A (zh) * 2011-07-18 2011-12-14 中国石油大学(北京) 一种煤层水力压裂模拟实验方法
CN202814926U (zh) * 2012-08-31 2013-03-20 中国矿业大学(北京) 一种充填开采相似模拟实验微观裂隙测试系统
CN103698171A (zh) * 2013-11-12 2014-04-02 山东科技大学 底板岩层固流耦合材料的裂隙演化试验方法
CN103674658A (zh) * 2013-12-24 2014-03-26 山东科技大学 一种随机裂隙试验模型的制备方法
CN104462654A (zh) * 2014-11-11 2015-03-25 中国矿业大学 浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
XU, HUIJUN ET AL.: "Numerical simulation research on gob air leakage of shallow-buried thin bedrock thick coal seam with fully-mechanized top coal caving technology.", JOURNAL OF CHINA COAL SOCIETY, vol. 36, no. 3, 31 March 2011 (2011-03-31), pages 435 - 441 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106869895A (zh) * 2017-04-07 2017-06-20 中国石油大学(华东) 一种室内煤岩水力压裂水泥环破坏的物模制作方法
CN108932365A (zh) * 2017-05-23 2018-12-04 通用电气公司 涡轮机润滑油分析器系统、计算机程序产品和相关方法
CN108932365B (zh) * 2017-05-23 2024-05-14 通用电气公司 涡轮机润滑油分析器系统、计算机程序产品和相关方法
CN111027787A (zh) * 2020-01-06 2020-04-17 中国矿业大学(北京) 一种基于断层滑移失稳诱发煤与瓦斯突出的实时预测方法
CN111027787B (zh) * 2020-01-06 2023-04-07 中国矿业大学(北京) 一种基于断层滑移失稳诱发煤与瓦斯突出的实时预测方法
CN112730730A (zh) * 2020-12-23 2021-04-30 长安大学 一种高速公路下伏燃烧煤层的模型试验装置及试验方法
CN112730730B (zh) * 2020-12-23 2023-02-03 长安大学 一种高速公路下伏燃烧煤层的模型试验装置及试验方法
CN113030430A (zh) * 2021-03-17 2021-06-25 中煤科工集团重庆研究院有限公司 煤矿采动稳定区储层内的游离煤层气浓度测算方法
CN113030430B (zh) * 2021-03-17 2023-02-17 中煤科工集团重庆研究院有限公司 煤矿采动稳定区储层内的游离煤层气浓度测算方法
CN113567061B (zh) * 2021-07-23 2024-01-12 安徽理工大学 一种近距离煤层开采采空区群漏风检测方法
CN113567061A (zh) * 2021-07-23 2021-10-29 安徽理工大学 一种近距离煤层开采采空区群漏风检测方法
CN114397421A (zh) * 2021-11-29 2022-04-26 中煤科工开采研究院有限公司 采动煤层底板破坏深度监测方法
CN115824853A (zh) * 2022-11-18 2023-03-21 国能神东煤炭集团有限责任公司 缩短地下水库与煤炭开采空间安全距离的方法
CN117072202A (zh) * 2023-08-22 2023-11-17 中国地质大学（北京） 一种黄土丘陵采煤塌陷区防渗堵漏处理方法

Also Published As

Publication number Publication date
CN104462654B (zh) 2017-06-13
ZA201606183B (en) 2017-09-27
CN104462654A (zh) 2015-03-25
AU2015345707B2 (en) 2017-10-19
AU2015345707A1 (en) 2016-12-22

Similar Documents

Publication Publication Date Title
CN104462654B (zh) 浅埋藏煤层开采地表贯通裂隙分布和漏风特征判定方法
CN105239611B (zh) 确定基坑开挖面下方止水帷幕渗漏对周边环境影响的方法
CN110005407B (zh) 盐穴储气库稳定性评价方法
CN103775071B (zh) 采动煤岩体裂隙演化的测量方法
CN109142192B (zh) 可视化异形固井二界面胶结质量测试系统
CN106437844B (zh) 一种超前预报隧道涌水位置的方法
CN102253179B (zh) 煤矿采场固-液-气三相介质耦合作用模拟实验装置
CN103821554B (zh) 基于无煤柱开采y型通风采空区的钻孔布置方法
CN105388265A (zh) 采空区瓦斯富集及运移的三维模拟方法
CN106919770A (zh) 一种基于数值模拟的损伤变量确定方法
Han et al. Response properties of geometries of coal penetrating fracture on seepage behavior
CN104266913A (zh) 一种矿井工作面底板采动破坏模拟试验装置
CN107288607A (zh) 一种利用Gd中子示踪产额成像评价近井压裂裂缝的方法
Qin et al. Optimization of abandoned gob methane drainage through well placement selection
Yu et al. Experimental study on the goaf flow field of the “U+ I” type ventilation system for a comprehensive mechanized mining face
CN202230077U (zh) 煤矿采场固-液-气三相介质耦合作用模拟实验装置
CN106503357A (zh) 基于穿层钻孔瓦斯涌出特征测定有效抽采半径的方法
Li et al. Comprehensive detection technique for coal seam roof water flowing fractured zone height
CN105388091A (zh) 一种矿用泡沫流体材料防灭火特性测试装置及方法
CN105386474A (zh) 确定基坑开挖面上方止水帷幕渗漏对周边环境影响的方法
CN110749533B (zh) 一种基于等效隔水层厚度的保水采煤判别方法
Zhao et al. Centrifuge modeling of loess slope failure induced by rising water level utilizing intact sample
CN102562052B (zh) 密井网浅层套损危害体识别方法
Liu et al. Internal instability evolution mechanism of tunnel face in sand-clay composite strata: Transparent soil model tests and DEM simulations
CN108303158A (zh) 一种浅埋煤层工作面采空区地表漏风量的计算方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15859217

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2015345707

Country of ref document: AU

Date of ref document: 20150811

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 15859217

Country of ref document: EP

Kind code of ref document: A1