[go: up one dir, main page]

CN108986208B - Reconstruction method for caving form of coal mine goaf - Google Patents

Reconstruction method for caving form of coal mine goaf Download PDF

Info

Publication number
CN108986208B
CN108986208B CN201810758623.3A CN201810758623A CN108986208B CN 108986208 B CN108986208 B CN 108986208B CN 201810758623 A CN201810758623 A CN 201810758623A CN 108986208 B CN108986208 B CN 108986208B
Authority
CN
China
Prior art keywords
goaf
coal mine
model
coal
lattice
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
Application number
CN201810758623.3A
Other languages
Chinese (zh)
Other versions
CN108986208A (en
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.)
Liaoning Technical University
Original Assignee
Liaoning Technical University
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 Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN201810758623.3A priority Critical patent/CN108986208B/en
Publication of CN108986208A publication Critical patent/CN108986208A/en
Application granted granted Critical
Publication of CN108986208B publication Critical patent/CN108986208B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/12Computing arrangements based on biological models using genetic models
    • G06N3/126Evolutionary algorithms, e.g. genetic algorithms or genetic programming

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Geometry (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Physiology (AREA)
  • Computer Graphics (AREA)
  • Computer Hardware Design (AREA)
  • Genetics & Genomics (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种煤矿采空区冒落形态的重构方法,涉及矿井安全技术领域。包括:建立煤矿采空区初始形态的物理相似模型,并在模型中布置传感器监测参数数据;随机生成煤矿采空区初始形态并形成煤矿采空区三维网格点阵;构建用于遗传算法的煤矿采空区染色体编码;采用格子Boltzmann法计算随机生成的煤矿采空区模拟模型的参数数据;判断计算结果的绝对误差与监测值的比值是否小于阈值;完成煤矿采空区冒落形态的重构。本发明提供的一种煤矿采空区冒落形态的重构方法,解决了以往建立在REV概念的基础上存在的缝洞组合、多尺度、大块度特征对新生采空区的影响,对研究煤矿采空区自然发火防治、瓦斯抽采及管理,保障煤矿安全生产具有重要的理论意义和实用价值。

Figure 201810758623

The invention provides a method for reconstructing the caving form of a goaf in a coal mine, and relates to the technical field of mine safety. Including: establishing a physical similar model of the initial shape of the goaf in the coal mine, and arranging sensor monitoring parameter data in the model; randomly generating the initial shape of the goaf in the coal mine and forming a three-dimensional grid lattice of the goaf in the coal mine; constructing a genetic algorithm Coal mine goaf chromosome coding; use the lattice Boltzmann method to calculate the parameter data of the randomly generated coal mine goaf simulation model; judge whether the ratio of the absolute error of the calculation result to the monitoring value is less than the threshold; structure. The invention provides a method for reconstructing the caving form of goafs in coal mines, which solves the influence of fracture-cave combination, multi-scale, and large-scale characteristics on the newly-born goafs based on the concept of REV in the past. It has important theoretical significance and practical value to study the prevention and control of natural fire in the goaf of coal mine, gas drainage and management, and to ensure the safe production of coal mine.

Figure 201810758623

Description

一种煤矿采空区冒落形态的重构方法A reconstruction method of caving form in coal mine goaf

技术领域technical field

本发明涉及矿井安全技术领域,具体涉及一种煤矿采空区冒落形态的重构方法。The invention relates to the technical field of mine safety, in particular to a method for reconstructing the caving form of a goaf in a coal mine.

背景技术Background technique

煤层开采后,采煤工作面后方老顶垮落步距尺度范围内直接顶、老顶在矿山压力作用下垮落,采空区呈现动态连续性变化,产生大量的缝隙和大空洞等不同尺度的空隙空间,形成了岩块间缝隙、浮煤碎块间的缝隙、大空洞和遗煤体孔隙交叠嵌套分布的多尺度复杂介质。而采煤工作面后方老顶垮落步距尺度范围内的采空区是自然发火高发区,也是影响采场瓦斯浓度分布的敏感区域。采空区内部流体真实流动情况是研究采空区自然发火及瓦斯分布的前提和基础,因此掌握复杂的采空区几何结构是极其必要的。After the coal seam is mined, the old roof at the back of the coal mining face collapses within the range of the step distance, and the direct roof and the old roof collapse under the action of the mine pressure, and the goaf shows dynamic continuous changes, resulting in a large number of gaps and large cavities of different scales. The interstitial spaces of the rocks form a multi-scale complex medium in which the gaps between rock blocks, the gaps between floating coal fragments, large cavities and the pores of coal remains are overlapped and nested. The goaf within the range of the caving step behind the coal mining face is a high-incidence area of spontaneous combustion, and it is also a sensitive area that affects the gas concentration distribution in the stope. The real flow of fluid inside the goaf is the premise and basis for the study of spontaneous combustion and gas distribution in the goaf, so it is extremely necessary to master the complex geometric structure of the goaf.

以往处理采空区流场假设采空区为连续介质,建立在REV概念的基础上,视为多孔介质范畴。目前描述多孔介质结构的模型主要有:连续介质模型,毛管束模型,网络模型,空间周期性,分形模型,过程法模型、统计重构模型,以及采用CT、扫描电镜以及切片组合等方法的真实结构模型等等,现有的采空区三维重构技术主要采用自动激光扫描系统获取点云数据,从而探测金属矿采空区,煤矿中重构技术主要应用于覆岩裂隙发育研究,揭示上覆岩层力学特性与渗流规律,汪北方、梁冰等假定岩体块体一致条件下选用Menger海绵分形模型建立采空区垮落岩体空隙结构模型,推导采空区垮落岩体空隙率和残余碎胀系数公式,计算了采空区垮落岩体空隙储水量。重构技术多用于岩石等多孔介质微观孔隙的三维重构上,而煤矿采空区随采随落,内部空间无法采用现有先进技术探测,新生采空区大块度和大空洞特点,也无法采用先进技术扫描或切片,因此对多尺度、大块度特征的新生采空区的三维重构还处于空白。In the past, the goaf flow field was treated assuming that the goaf was a continuous medium, which was based on the concept of REV and regarded as the category of porous media. At present, the models describing the structure of porous media mainly include: continuum model, capillary bundle model, network model, spatial periodicity, fractal model, process method model, statistical reconstruction model, and real models using CT, scanning electron microscopy, and slice combination methods. Structural models, etc. The existing 3D reconstruction technology of goafs mainly uses automatic laser scanning The mechanical properties and seepage laws of the overlying strata, Wang Beibei, Liang Bing, etc. used the Menger sponge fractal model to establish the void structure model of the collapsed rock mass in the goaf under the assumption that the rock mass blocks were consistent, and deduced the void ratio and The residual dilation coefficient formula calculates the water storage volume of the collapsed rock mass in the goaf. Reconstruction technology is mostly used for three-dimensional reconstruction of microscopic pores in porous media such as rocks. However, goafs in coal mines fall as they are mined, and the internal space cannot be detected by existing advanced technologies. It is impossible to use advanced technology to scan or slice, so the 3D reconstruction of the new goaf with multi-scale and large-scale features is still blank.

受煤岩物理力学特性、节理、采煤方法、采高等各种因素的影响,煤岩碎块形状各式各样,形成的空洞空间更是千变万化,即使同一个回采工作面其新生采空区的空洞形态也是动态随机变化的。可以说,无论是从时间上、还是从空间上都找不出两个完全相同的新生采空区空洞形态。此外,存在以下问题:(1)岩块的真实空洞结构复杂,并且空洞和缝隙跨尺度、无序分布,难以通过实验手段获得完备的采空区空洞真实结构信息,从理论上也难以准确地定量描述,缺乏相应的数学、物理描述。(2)采用简化的采空区空洞分布或者采用几个简单的统计特征函数来确定空洞结构分布,不能真实表达采空区空洞结构,重构采空区空洞与真实采空区之间必然存在差距。(3)由于试验条件和设备的限制,目前还找不到合理的试验方法来了解采空区真实的内部空洞结构。Affected by various factors such as coal and rock physical and mechanical properties, joints, coal mining methods, and mining heights, the shapes of coal and rock fragments are various, and the void space formed is ever-changing. The shape of the void is also dynamically and randomly changed. It can be said that there are no two completely identical hollow forms in new mined-out areas, either in terms of time or space. In addition, there are the following problems: (1) The real cavity structure of the rock block is complex, and the cavities and cracks are distributed across scales and disorderly. Quantitative description, lack of corresponding mathematical and physical description. (2) Using the simplified goaf cavity distribution or using several simple statistical characteristic functions to determine the cavity structure distribution cannot truly express the cavity structure of the goaf, and there must be a gap between the reconstructed goaf cavity and the real goaf. gap. (3) Due to the limitations of test conditions and equipment, there is no reasonable test method to understand the real internal cavity structure of the goaf.

发明内容Contents of the invention

针对现有技术存在的问题,本发明提供一种煤矿采空区冒落形态的重构方法,利用遗传算法来“穷举猜测”新生采空区冒落形态,采用格子Boltzmann法对猜测出来的新生采空区的空洞流进行解算,解算结果与实际采空区空洞流的气体压力、组份气体浓度分布进行比较,从而得到新生采空区内部结构冒落形态。Aiming at the problems existing in the prior art, the present invention provides a method for reconstructing the caving form of goafs in coal mines. Genetic algorithms are used to "exhaustively guess" the caving forms of newly mined-out areas, and the lattice Boltzmann method is used to estimate the caving forms of goafs. The cavity flow in the new mined-out area is calculated, and the calculation results are compared with the gas pressure and component gas concentration distribution of the actual goaf flow, so as to obtain the caving form of the internal structure of the new mined-out area.

为了实现上述目的,一种煤矿采空区冒落形态的重构方法,包括以下步骤:In order to achieve the above purpose, a method for reconstructing the caving form of a goaf in a coal mine includes the following steps:

步骤1:通过在空间内摆放固体方块的方式建立煤矿采空区初始形态的物理相似模型,并在模型中布置传感器监测参数数据,包括压力、速度和组分浓度;Step 1: Establish a physical similarity model of the initial shape of the coal mine goaf by placing solid blocks in the space, and arrange sensors to monitor parameter data in the model, including pressure, velocity and component concentration;

步骤2:随机生成煤矿采空区的初始形态,将其进行三维网格划分,形成煤矿采空区三维网格点阵,其中,将煤矿采空区三维网格点阵中被煤岩体填充的网格设置为1,未被填充的气相介质设置为0;Step 2: Randomly generate the initial shape of the goaf in the coal mine, divide it into a three-dimensional grid to form a three-dimensional grid lattice of the goaf in the coal mine, and fill the three-dimensional grid lattice of the goaf in the coal mine with coal and rock mass The grid is set to 1, and the unfilled gas phase medium is set to 0;

步骤3:构建用于遗传算法的新生煤矿采空区染色体编码,根据遗传算法随机生成0-1序列的煤矿采空区模拟模型;Step 3: Construct the goaf chromosomal code for the newborn coal mine using the genetic algorithm, and randomly generate the 0-1 sequence of the goaf simulation model of the coal mine according to the genetic algorithm;

步骤4:采用格子Boltzmann法计算随机生成的煤矿采空区模拟模型内部的压力、速度和组分浓度;Step 4: Use the lattice Boltzmann method to calculate the pressure, velocity and component concentration inside the randomly generated coal mine goaf simulation model;

步骤5:判断格子Boltzmann法计算结果的绝对误差与物理相似模型中传感器所监测的参数数据的比值是否小于相对误差阈值,若是,则输出煤矿采空区三维重构模型,即完成煤矿采空区冒落形态的重构,若否,则返回步骤3。Step 5: Determine whether the ratio of the absolute error of the calculation result of the lattice Boltzmann method to the parameter data monitored by the sensor in the physical similarity model is less than the relative error threshold, and if so, output the 3D reconstruction model of the coal mine goaf, that is, complete the coal mine goaf The reconstruction of the falling form, if not, return to step 3.

进一步地,所述步骤5中格子Boltzmann法计算结果的绝对误差的计算公式如下:Further, the calculation formula of the absolute error of the calculation result of the lattice Boltzmann method in the step 5 is as follows:

Δ=|C-L|;Δ=|C-L|;

其中,Δ为格子Boltzmann法计算结果的绝对误差,C为格子Boltzmann法计算值,L为物理相似模型监测点监测参数数据。Among them, Δ is the absolute error of the calculation result of the lattice Boltzmann method, C is the calculation value of the lattice Boltzmann method, and L is the monitoring parameter data of the monitoring points of the physical similarity model.

进一步地,所述步骤5中相对误差阈值根据实际情况具体设定。Further, the relative error threshold in step 5 is specifically set according to the actual situation.

本发明的有益效果:Beneficial effects of the present invention:

本发明提出一种煤矿采空区冒落形态的重构方法,克服了传统技术在进行复杂结构的流动模拟时的局限性,解决了以往建立在REV概念的基础上存在的缝洞组合、多尺度、大块度特征对新生采空区的影响,本发明根据新生煤矿采空区冒落形态重构的采空区气相介质运动过程、弥散过程和传质传热过程,获得能够与现场实际情况相一致的气体流动场、组份气体浓度分布场和温度分布场状态参数,对研究煤矿采空区自然发火防治、瓦斯抽采及管理,保障煤矿安全生产具有重要的理论意义和实用价值。The present invention proposes a method for reconstructing the caving form of gobs in coal mines, which overcomes the limitations of traditional techniques in the flow simulation of complex structures, and solves the problems of fracture-cave combinations, multiple The impact of scale and large block characteristics on the new goaf, the present invention reconstructs the movement process, dispersion process and mass transfer and heat transfer process of the gas phase medium in the goaf according to the caving form of the goaf in the new coal mine, and obtains a process that can be compared with the actual situation on site. The state parameters of gas flow field, component gas concentration distribution field and temperature distribution field with consistent conditions have important theoretical significance and practical value for the study of spontaneous combustion prevention, gas drainage and management in goafs of coal mines, and for ensuring safe production in coal mines.

附图说明Description of drawings

图1为本发明实施例中煤矿采空区冒落形态的重构方法流程图;Fig. 1 is the flow chart of the reconstruction method of the caving form of the coal mine goaf in the embodiment of the present invention;

图2为本发明实施例中构建的煤矿采空区初始形态物理相似模型示意图;Fig. 2 is the schematic diagram of the initial form physical similarity model of the coal mine goaf constructed in the embodiment of the present invention;

图3为本发明实施例中物理相似模型中布置传感器的位置示意图;Fig. 3 is a schematic diagram of the positions of the sensors arranged in the physical similarity model in the embodiment of the present invention;

图4为本发明实施例中最终重构的煤矿采空区三维模型示意图。Fig. 4 is a schematic diagram of the finally reconstructed three-dimensional model of the coal mine goaf in the embodiment of the present invention.

具体实施方式Detailed ways

为了使本发明的目的、技术方案及优势更加清晰,下面结合附图和具体实施例对本发明做进一步详细说明。此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。In order to make the purpose, technical solution and advantages of the present invention clearer, the present invention will be further described in detail below in conjunction with the accompanying drawings and specific embodiments. The specific embodiments described here are only used to explain the present invention, not to limit the present invention.

一种煤矿采空区冒落形态的重构方法,流程如图1所示,具体方法如下所述:A method for reconstructing the caving form of gobs in coal mines, the process flow is shown in Figure 1, and the specific methods are as follows:

步骤1:通过在空间内摆放固体方块的方式建立煤矿采空区初始形态的物理相似模型,并在模型中布置传感器监测参数数据,包括压力、速度和组分浓度。Step 1: Establish a physical similarity model of the initial form of the coal mine goaf by placing solid blocks in the space, and arrange sensors in the model to monitor parameter data, including pressure, velocity and component concentration.

本实施例中,建立的煤矿采空区初始形态物理相似模型如图2所示,布置传感器P的位置如图3所示,其中,P1-P12为布置的12个传感器。In this embodiment, the established physical similarity model of the initial shape of the coal mine goaf is shown in Figure 2, and the position of the sensor P is shown in Figure 3, where P1-P12 are 12 sensors arranged.

步骤2:随机生成煤矿采空区的初始形态,将其进行三维网格划分,形成煤矿采空区三维网格点阵,其中,将煤矿采空区三维网格点阵中被煤岩体填充的网格设置为1,未被填充的气相介质设置为0。Step 2: Randomly generate the initial shape of the goaf in the coal mine, divide it into a three-dimensional grid to form a three-dimensional grid lattice of the goaf in the coal mine, and fill the three-dimensional grid lattice of the goaf in the coal mine with coal and rock mass The grid is set to 1, and the unfilled gas phase medium is set to 0.

如图2所示,本实施例中,网格为被煤岩体填充区域,空白格为未被填充区域。As shown in FIG. 2 , in this embodiment, the grid is the area filled with coal and rock mass, and the blank grid is the unfilled area.

步骤3:构建用于遗传算法的新生煤矿采空区染色体编码,根据遗传算法随机生成0-1序列的煤矿采空区模拟模型。Step 3: Construct the goaf chromosomal code for the newborn coal mine using the genetic algorithm, and randomly generate a 0-1 sequence simulation model of the goaf in the coal mine according to the genetic algorithm.

本实施例中,构建的用于遗传算法的新生煤矿采空区染色体编码为:

Figure BDA0001727348440000031
In this example, the newly constructed goaf chromosome code for the genetic algorithm is:
Figure BDA0001727348440000031

步骤4:采用格子Boltzmann法计算随机生成的煤矿采空区模拟模型内部的压力、速度和组分浓度。Step 4: Use the lattice Boltzmann method to calculate the pressure, velocity and component concentration inside the randomly generated coal mine goaf simulation model.

所述格子Boltzmann法简称LBM,是一种基于介观模拟尺度的计算流体力学方法,该方法相比于其他传统CFD计算方法,具有介于微观分子动力学模型和宏观连续模型的介观模型特点,因此具备流体相互作用描述简单、复杂边界易于设置、易于并行计算、程序易于实施等优势,是一种描述流体运动的有效手段。The Lattice Boltzmann method, referred to as LBM, is a computational fluid dynamics method based on the mesoscopic simulation scale. Compared with other traditional CFD calculation methods, this method has the characteristics of a mesoscopic model between the microscopic molecular dynamics model and the macroscopic continuous model. , so it has the advantages of simple description of fluid interaction, easy setting of complex boundaries, easy parallel computing, and easy implementation of programs, and is an effective means to describe fluid motion.

步骤5:判断格子Boltzmann法计算结果的绝对误差与物理相似模型中传感器所监测的参数数据的比值是否小于相对误差阈值,若是,则输出煤矿采空区三维重构模型,即完成煤矿采空区冒落形态的重构,若否,则返回步骤3。Step 5: Determine whether the ratio of the absolute error of the calculation result of the lattice Boltzmann method to the parameter data monitored by the sensor in the physical similarity model is less than the relative error threshold, and if so, output the 3D reconstruction model of the coal mine goaf, that is, complete the coal mine goaf The reconstruction of the falling form, if not, return to step 3.

所述格子Boltzmann法计算结果的绝对误差的计算公式如公式(1)所示:The calculation formula of the absolute error of the calculation result of the lattice Boltzmann method is as shown in formula (1):

Δ=|C-L|   (1)Δ=|C-L| (1)

其中,Δ为格子Boltzmann法计算结果的绝对误差,C为格子Boltzmann法计算值,L为物理相似模型监测点监测参数数据。Among them, Δ is the absolute error of the calculation result of the lattice Boltzmann method, C is the calculation value of the lattice Boltzmann method, and L is the monitoring parameter data of the monitoring points of the physical similarity model.

本实施例中,相对误差阈值设定为0.05。In this embodiment, the relative error threshold is set to 0.05.

本实施例中,监测位置对应的监测值和采用格子Boltzmann法计算得到的数值结果如表1所示,最终生成的煤矿采空区染色体编码为:

Figure BDA0001727348440000041
In this embodiment, the monitoring values corresponding to the monitoring positions and the numerical results calculated by using the lattice Boltzmann method are shown in Table 1, and the chromosome code of the finally generated coal mine goaf is:
Figure BDA0001727348440000041

表1监测位置对应的监测值和采用格子Boltzmann法计算得到的数值结果Table 1 The monitoring values corresponding to the monitoring positions and the numerical results calculated by the lattice Boltzmann method

Figure BDA0001727348440000042
Figure BDA0001727348440000042

通过公式(1)计算得到格子Boltzmann法计算结果的绝对误差与物理相似模型中传感器所监测的参数数据的比值小于相对误差阈值,因此输出最终重构的煤矿采空区三维模型,如图4所示,完成煤矿采空区冒落形态的重构。Calculated by formula (1), the ratio of the absolute error of the calculation result of the lattice Boltzmann method to the parameter data monitored by the sensor in the physical similarity model is less than the relative error threshold, so the final reconstructed 3D model of the coal mine goaf is output, as shown in Figure 4 It shows that the reconstruction of the caving form of the coal mine goaf is completed.

最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解;其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;因而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand; it still can Modifications are made to the technical solutions described in the foregoing embodiments, or equivalent replacements are made to some or all of the technical features; therefore, these modifications or replacements do not make the essence of the corresponding technical solutions depart from the scope defined by the claims of the present invention.

Claims (3)

1.一种煤矿采空区冒落形态的重构方法,其特征在于,包括以下步骤:1. A reconstruction method of the caving form of a coal mine goaf, characterized in that, comprising the following steps: 步骤1:通过在空间内摆放固体方块的方式建立煤矿采空区初始形态的物理相似模型,并在模型中布置传感器监测参数数据,包括压力、速度和组分浓度;Step 1: Establish a physical similarity model of the initial shape of the coal mine goaf by placing solid blocks in the space, and arrange sensors to monitor parameter data in the model, including pressure, velocity and component concentration; 步骤2:随机生成煤矿采空区的初始形态,将其进行三维网格划分,形成煤矿采空区三维网格点阵,其中,将煤矿采空区三维网格点阵中被煤岩体填充的网格设置为1,未被填充的气相介质设置为0,构建用于遗传算法的新生煤矿采空区染色体编码;Step 2: Randomly generate the initial shape of the goaf in the coal mine, divide it into a three-dimensional grid to form a three-dimensional grid lattice of the goaf in the coal mine, and fill the three-dimensional grid lattice of the goaf in the coal mine with coal and rock mass The grid of the grid is set to 1, and the unfilled gas phase medium is set to 0, and the chromosome coding of the mined-out area of the nascent coal mine is constructed for the genetic algorithm; 步骤3:构建用于遗传算法的新生煤矿采空区染色体编码,根据遗传算法随机生成0-1序列的煤矿采空区模拟模型;Step 3: Construct the goaf chromosomal code for the newborn coal mine using the genetic algorithm, and randomly generate the 0-1 sequence of the goaf simulation model of the coal mine according to the genetic algorithm; 步骤4:采用格子Boltzmann法计算所述步骤3中随机生成的煤矿采空区模拟模型内部气体的压力、速度和组分浓度;Step 4: Using the lattice Boltzmann method to calculate the pressure, velocity and component concentration of the gas inside the coal mine goaf simulation model randomly generated in the step 3; 步骤5:判断所述步骤4中格子Boltzmann法计算结果的绝对误差与所述步骤1中物理相似模型中传感器所监测的参数数据的比值是否小于相对误差阈值,若是,则输出煤矿采空区三维重构模型,即完成煤矿采空区冒落形态的重构,若否,则返回步骤3。Step 5: Determine whether the ratio of the absolute error of the calculation result of the lattice Boltzmann method in the step 4 to the parameter data monitored by the sensor in the physical similarity model in the step 1 is less than the relative error threshold, and if so, output the three-dimensional coal mine goaf Reconstruct the model, that is, complete the reconstruction of the caving form of the coal mine goaf, if not, return to step 3. 2.根据权利要求1所述的煤矿采空区冒落形态的重构方法,其特征在于,所述步骤5中格子Boltzmann法计算结果的绝对误差的计算公式如下:2. the reconstruction method of coal mine goaf caving form according to claim 1, is characterized in that, the calculation formula of the absolute error of lattice Boltzmann method calculation result in the described step 5 is as follows: Δ=|C-L|;Δ=|C-L|; 其中,Δ为格子Boltzmann法计算结果的绝对误差,C为格子Boltzmann法计算值,L为物理相似模型监测点监测参数数据。Among them, Δ is the absolute error of the calculation result of the lattice Boltzmann method, C is the calculation value of the lattice Boltzmann method, and L is the monitoring parameter data of the monitoring points of the physical similarity model. 3.根据权利要求1所述的煤矿采空区冒落形态的重构方法,其特征在于,所述步骤5中相对误差阈值根据实际情况具体设定。3. The method for reconstructing the caving form of coal mine goaf according to claim 1, characterized in that the relative error threshold in step 5 is specifically set according to actual conditions.
CN201810758623.3A 2018-07-11 2018-07-11 Reconstruction method for caving form of coal mine goaf Active CN108986208B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810758623.3A CN108986208B (en) 2018-07-11 2018-07-11 Reconstruction method for caving form of coal mine goaf

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810758623.3A CN108986208B (en) 2018-07-11 2018-07-11 Reconstruction method for caving form of coal mine goaf

Publications (2)

Publication Number Publication Date
CN108986208A CN108986208A (en) 2018-12-11
CN108986208B true CN108986208B (en) 2023-04-07

Family

ID=64536955

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810758623.3A Active CN108986208B (en) 2018-07-11 2018-07-11 Reconstruction method for caving form of coal mine goaf

Country Status (1)

Country Link
CN (1) CN108986208B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111125968B (en) * 2019-12-24 2023-04-14 辽宁工程技术大学 A LDA-based regression method for solid-phase structure and flow field in new mined-out area
CN113902863B (en) * 2021-10-12 2023-08-29 华北科技学院 Method for representing equivalent pore network topology structure of coal
CN114266160B (en) * 2021-12-23 2024-06-04 华中科技大学 Round coal pile internal temperature simulation method and system
CN115935120B (en) * 2023-01-08 2023-06-13 中煤科工西安研究院(集团)有限公司 Goaf Scale Measurement and Calculation Method Based on Miaozi Remote Sensing Data

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB682981A (en) * 1950-04-03 1952-11-19 Markham & Company Ltd An improved method for forming packs between the floor and the roof of a mine
JP2004151790A (en) * 2002-10-28 2004-05-27 Mie Tlo Co Ltd Nurse wopk schedule preparation support software or the like
CN102053945A (en) * 2009-11-09 2011-05-11 中国科学院过程工程研究所 Concurrent computational system for multi-scale discrete simulation
CN103091675A (en) * 2013-01-11 2013-05-08 中南大学 Mining lot exploiting and monitoring method based on interferometric synthetic aperature radar (InSAR) technology
CN106682282A (en) * 2016-12-08 2017-05-17 浙江大学 Wind farm multi-type wind driven generator arrangement optimization method based on genetic algorithm nested in particle swarm optimization
CN107092740A (en) * 2017-04-12 2017-08-25 乌鲁木齐城市轨道集团有限公司 A kind of old goaf residual deformation method for predicting of severe inclined thick coal seam
CN107818425A (en) * 2017-11-23 2018-03-20 辽宁工程技术大学 A kind of dangerous evaluation method of gas explosion at longwall face
CN107942380A (en) * 2017-10-23 2018-04-20 太原理工大学 A kind of numerical simulation method for considering caving zone goaf
CN108128856A (en) * 2018-02-05 2018-06-08 中国矿业大学 Utilize the system of principle of electrophoresis control goaf filling weight metal transport
CN111125968A (en) * 2019-12-24 2020-05-08 辽宁工程技术大学 A solid phase structure and flow field regression method for newly mined goaf based on LDA
CN113705130A (en) * 2021-09-27 2021-11-26 辽宁工程技术大学 Method for determining flow similarity of hybrid ventilation model experiment of excavation roadway
CN114722742A (en) * 2022-05-16 2022-07-08 中国矿业大学(北京) Method for measuring gas dispersion coefficient of porous medium in goaf

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU771240B2 (en) * 1998-09-14 2004-03-18 Aluminium Extrusion And Distribution Pty Limited Lock for double hung sash window
US9892225B2 (en) * 2016-04-01 2018-02-13 International Business Machines Corporation Method for optimizing the design of micro-fluidic devices

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB682981A (en) * 1950-04-03 1952-11-19 Markham & Company Ltd An improved method for forming packs between the floor and the roof of a mine
JP2004151790A (en) * 2002-10-28 2004-05-27 Mie Tlo Co Ltd Nurse wopk schedule preparation support software or the like
CN102053945A (en) * 2009-11-09 2011-05-11 中国科学院过程工程研究所 Concurrent computational system for multi-scale discrete simulation
CN103091675A (en) * 2013-01-11 2013-05-08 中南大学 Mining lot exploiting and monitoring method based on interferometric synthetic aperature radar (InSAR) technology
CN106682282A (en) * 2016-12-08 2017-05-17 浙江大学 Wind farm multi-type wind driven generator arrangement optimization method based on genetic algorithm nested in particle swarm optimization
CN107092740A (en) * 2017-04-12 2017-08-25 乌鲁木齐城市轨道集团有限公司 A kind of old goaf residual deformation method for predicting of severe inclined thick coal seam
CN107942380A (en) * 2017-10-23 2018-04-20 太原理工大学 A kind of numerical simulation method for considering caving zone goaf
CN107818425A (en) * 2017-11-23 2018-03-20 辽宁工程技术大学 A kind of dangerous evaluation method of gas explosion at longwall face
CN108128856A (en) * 2018-02-05 2018-06-08 中国矿业大学 Utilize the system of principle of electrophoresis control goaf filling weight metal transport
CN111125968A (en) * 2019-12-24 2020-05-08 辽宁工程技术大学 A solid phase structure and flow field regression method for newly mined goaf based on LDA
CN113705130A (en) * 2021-09-27 2021-11-26 辽宁工程技术大学 Method for determining flow similarity of hybrid ventilation model experiment of excavation roadway
CN114722742A (en) * 2022-05-16 2022-07-08 中国矿业大学(北京) Method for measuring gas dispersion coefficient of porous medium in goaf

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
An inverse analysis of a transient 2-D conduction–radiation problem using the lattice Boltzmann method and the finite volume method coupled with the genetic algorithm;RanjanDas et al.;《Journal of Quantitative Spectroscopy and Radiative Transfer》;20080731(第11期);2060-2077 *
Hybrid genetic algorithm for the optimization of mine ventilation network;ZHAO DAN et al.;《JOURNAL OF COAL SCIENCE & ENGINEERING(CHINA)》;20091231;第15卷(第4期);389-393 *
INVERSE AEROACOUSTIC DESIGN OF AXIAL FANS USING GENETIC OPTIMIZATION AND THE LATTICE-BOLTZMANN METHOD;Michael Stadler et al.;《Proceedings of ASME Turbo Expo 2013: Turbine Technical Conference and Exposition GT2013》;20131114;1-11 *
地表沉陷预计模型及应用研究;赵迎春 等;《能源技术与管理》;20170430(第4期);127-129 *
基于反演参数的采空区坡沿隆起治理研究;姚培峰等;《武汉理工大学学报》;20170228(第02期);39-44 *
基于格子Boltzmann的采空区顶板冒落仿真模型;薛兵 等;《纺织高校基础科学学报》;20140331(第3期);392-397 *
用遗传算法优化boltzmann机;陈洁 等;《计算机工程与应用》;20010731(第7期);46-48 *

Also Published As

Publication number Publication date
CN108986208A (en) 2018-12-11

Similar Documents

Publication Publication Date Title
CN108986208B (en) Reconstruction method for caving form of coal mine goaf
RU2743121C1 (en) Method of numerical modeling of filtration/damage/stress communication during water injection into coal-bearing mass during regionalization
Zhang et al. Failure mechanism and numerical simulation of zonal disintegration around a deep tunnel under high stress
CN105550441B (en) A kind of engineering rock mass rupture deterioration method for numerical simulation based on continuous media
Hao et al. Bulking factor of the strata overlying the gob and a three-dimensional numerical simulation of the air leakage flow field
CN109241588A (en) A kind of analogy method of the monolete extension based on quasi-continuous geomechanics model
Wilmarth et al. Power density in geothermal fields
CN106649963B (en) Volume fracturing complexity seam net average crack length and equivalent fissure item number determine method
CN114510854B (en) A method for evaluating the accuracy of numerical simulation results of acid fracturing for finding holes by following fractures
CN101476458A (en) Oil pool development simulation system, oil pool model body and its data processing method
CN106226813B (en) Pressure-break net reconstructing method and device based on microseism
NO20110986A1 (en) System and method for predicting fluid flow properties in fractured reservoirs below the surface
NO20110593A1 (en) Multiphase flow in a wellbore and associated hydraulic fracturing
CN104895550A (en) A method for establishing and solving numerical well test models for tight gas fractured horizontal wells
CN103279991A (en) Method for improving oil reservoir exploitation effect by utilizing discrete fracture endpoint deformation numerical simulation
CN105787220A (en) Coal bed high-pressure water injection fracturing-flow seeping value simulation method
CN105512366A (en) Natural-fracture-containing tight-reservoir volume-fracturing tree-shape random fracture net description method
Grant et al. Recovery factor for EGS
NO20131641A1 (en) MULTI-SEGMENT CRACKS
CN116842784A (en) Prediction method for the minimum safe distance between double-line shield tunnels and caves in water-rich karst areas
Aydın et al. Numerical reservoir simulation of Alaşehir geothermal field
CN110259443A (en) A kind of coal seam borehole wall stability prediction method based on 3DEC discrete element
CN110080745B (en) Method and device for predicting productivity of vertical well through separate-layer fracturing
Liu et al. The unsaturated hydromechanical coupling model of rock slope considering rainfall infiltration using DDA
RU2739287C1 (en) Method for analysis and design of stimulation of an underground formation

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