Automated identification of fissure trace in mining roadway via deep learning

The fractured rocks distributed near the structural surface of coal roadway heading face can easily cause water or air conduits,ultimately resulting in accidents that threaten production safety,such as mine flooding,gas leakage,or supporting failure.Therefore,successful detection of fissures is crucial to mine safety.Despite the rising popularity of computer vision in accurate fissures detection,it fails to satisfy the demand of engineering practice.To address this problem,this paper first establishes a 1000-image database of fissures of coal roadway heading face based on data collection,cleaning,and annotation.Then,a framework for fissure detection and segmentation is constructed with deep convolutional neural network(DCNN)named DeepLabv3+ serving as the overall architecture,and a lightweight MobileNetV2,instead of the original Xception,as the main feature extraction network.The database is then employed to train and test the neural network model.Finally,the robustness and adaptability of the model under the common jamming environment in coal mines are evaluated.According to the results,the deep learning algorithm,which performs favorably in identifying various fissures in the coal roadway heading face,is immune to interference such as low illumination,wire mesh,multi-scale edge,cutting mark,and concentrated light beams.Notably,the performance of such method is on par with,or better than,humans in performing individual image segmentation.The approach performs higher segmentation accuracy and calculation speed than the traditional image identification algorithm,thus realizing rapid identification of fissures in coal mines in batch.This study can provide a reference for image semantic segmentation of fissure traces under similar conditions.

基于注汞法的液冷花岗岩微观特性实验研究

在开采高温储层的过程中,微观特征变化重要影响深部热储层人工裂隙网络的建立。因此,研究该过程中微观特征的变化对深部地热能的安全高效开采具有重要意义。采用液体冷却(包括酸冷却和水冷却)处理高温花岗岩,模拟向高温储层注酸情况,以获取相应条件下的花岗岩试样。随后,使用扫描电子显微镜-能量散射谱(SEM-EDS)和高压注汞法对样品进行测试。SEM-EDS结果显示,水冷却几乎不会影响元素的组成,酸冷却导致氧、硅和金属元素减少。在两种情况下,都观察到了孔隙和裂缝的形成。此外,与水冷却相比,酸冷却的样品表面更加不平整、疏松和粗糙。注汞法测试结果显示,随着处理温度的升高,液冷花岗岩样品的孔隙率、孔隙体积和比表面积有增加,但分形维数呈相反的趋势。研究结果表明,与水冷却相比,酸冷却条件对花岗岩的腐蚀作用使各项参数的变化更加显著,在酸冷却过程中产生更多的孔隙,提供更好的运输空间,使孔隙结构的复杂性降低。在高温储层的地热能有效提取方面,酸冷却具有更大应用潜力。

高压电脉冲刺激对于热处理花岗岩影响的实验研究

高压电脉冲技术作为一种环保、可控、高效的破碎坚硬岩石的方法,其潜在应用场景值得进一步深入发掘。本文研究了高压电脉冲诱导干热岩产生预损伤和劣化方法的效果。对热处理后花岗岩薄片试样施加了不同峰值电压的高压电脉冲刺激,探讨了高压电脉冲对于热处理花岗岩的影响。此外,还研究了初始温度对高压电脉冲刺激花岗岩的影响。结果表明:当外加峰值电压超过96 kV时,热处理花岗岩试样发生贯穿性的破裂,峰值电压越高,贯穿裂纹网络越复杂。当电压从96 kV增加到144kV时,花岗岩试样的微裂纹密度从8.63 mm/mm^(2)增加到13.26 mm/mm^(2)。随着热处理温度的升高,花岗岩发生电击穿的难度逐渐降低。在120 kV电脉冲放电刺激3次后,所有热处理温度高于400℃的花岗岩试样均能发生贯穿破裂。高压电脉冲刺激造成花岗岩产生的最大损伤发生在300~400℃的温度区间内。通过扫描电子显微镜观察到,随着试样处理温度的升高,内部孔隙和裂缝的发育情况逐渐加剧,孔隙大小、裂缝宽度和深度也随之增加。这项研究有望为建造干热岩人工储层提供新的见解。

Performances of fissured red sandstone after thermal treatment with constant-amplitude and low-cycle impacts

In the engineering practices,it is increasingly common to encounter fractured rocks perturbed by temperatures and frequent dynamic loads.In this paper,the dynamic behaviors and fracture characteristics of red sandstone considering temperatures(25℃,200℃,400℃,600℃,and 800℃)and fissure angles(0°,30°,60°,and 90°)were evaluated under constant-amplitude and low-cycle(CALC)impacts actuated by a modified split Hopkinson pressure bar(SHPB)system.Subsequently,fracture morphology and second-order statistics within the grey-level co-occurrence matrix(GLCM)were examined using scanning electron microscopy(SEM).Meanwhile,the deep analysis and discussion of the mechanical response were conducted through the synchronous thermal analyzer(STA)test,numerical simulations,one-dimensional stress wave theory,and material structure.The multiple regression models between response variables and interactive effects of independent variables were established using the response surface method(RSM).The results demonstrate the fatigue strength and life diminish as temperatures rise and increase with increasing fissure angles,while the strain rate exhibits an inverse behavior.Furthermore,the peak stress intensification and strain rate softening observed during CALC impact exhibit greater prominence at increased fissure angles.The failure is dominated by tensile damage with concise evolution paths and intergranular cracks as well as the compressor-crushed zone which may affect the failure mode after 400℃.The second-order statistics of GLCM in SEM images exhibit a considerable dependence on the temperatures.Also,thermal damage dominated by thermal properties controls the material structure and wave impedance and eventually affects the incident wave intensity.The tensile wave reflected from the fissure surface is the inherent mechanism responsible for the angle effect exhibited by the fatigue strength and life.Ultimately,the peak stress intensification and strain rate softening during impact are determined by both the material structure and compaction governed by thermal damage and tensile wave.

Mechanical responses of anchoring structure under triaxial cyclic loading

Dynamic load on anchoring structures(AS)within deep roadways can result in cumulative damage and failure.This study develops an experimental device designed to test AS under triaxial loads.The device enables the investigation of the mechanical response,failure mode,instability assessment criteria,and anchorage effect of AS subjected to combined cyclic dynamic-static triaxial stress paths.The results show that the peak bearing strength is positively correlated with the anchoring matrix strength,anchorage length,and edgewise compressive strength.The bearing capacity decreases significantly when the anchorage direction is severely inclined.The free face failure modes are typically transverse cracking,concave fracturing,V-shaped slipping and detachment,and spallation detachment.Besides,when the anchoring matrix strength and the anchorage length decrease while the edgewise compressive strength,loading rate,and anchorage inclination angle increase,the failure intensity rises.Instability is determined by a negative tangent modulus of the displacement-strength curve or the continued deformation increase against the general downward trend.Under cyclic loads,the driving force that breaks the rock mass along the normal vector and the rigidity of the AS are the two factors that determine roadway stability.Finally,a control measure for surrounding rock stability is proposed to reduce the internal driving force via a pressure relief method and improve the rigidity of the AS by full-length anchorage and grouting modification.

Implications for fault reactivation and seismicity induced by hydraulic fracturing

Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.

Kinetic energy and its applications in mining engineering

Reduction of energy consumption in comminution is of significant importance in mining industry. To reduce such energy consumption the energy efficiency in an individual operation such as blasting must be increased. By using both new investigations and previous experimental results, this paper demonstrates that (1) kinetic energy carried by moving fragments in rock fracture is notable and it increases with an increasing loading rate;(2) this kinetic energy can be well used in secondary fragmentation in crushing and blasting. Accordingly, part of the muck pile from previous blast should be left in front of new(bench) face in either open pit or underground blasting. If so, when new blast occurs, the fragments from the new blast will collide with the muck pile left from the previous blast, and the kinetic energy carried by the moving fragments will be partly used in their secondary fragmentation.

Stress distribution and surrounding rock control of mining near to the overlying coal pillar in the working face

The occurrence of overlying coal pillar(OCP)exerts a strong effect on the stress and strain distribution of the surrounding rock in the stope.In this paper,the stress distribution characteristics are analyzed via the numerical calculation with the account of OCP presence or absence.In addition,this study revealed the joint effect of side pressure relief area of the goaf and stress concentration in OCP on the final stress distribution.Furthermore,the rules of abutment stress distribution affected by three influencing factors,namely horizontal-vertical distances between OCP and working face and buried depth of OCP,are analyzed.The functional model linking the peak stress of surrounding rock with the above influencing factors is developed.The field application of the above results proved that the rib spalling and deformation of a 2.95 m-high and 5.66 m-wide roadway could be efficiently controlled by rationally adjusting working states of the support,and adopting the hydraulic prop coordinated with the p type metal beam and anchor cable to strengthen the surrounding rock of working face and roadway,respectively.The proposed measures are considered appropriate to satisfy the safe operation requirements.

Green coal mining technique integrating mining-dressing-gas draining-backfilling-mining

Aiming to address the following major engineering issues faced by the Pingdingshan No. 12 mine:(1) difficulty in implementing auxiliary lifting because of its depth(i.e., beyond 1000 m);(2) highly gassy main coal seam with low permeability;(3) unstable overlying coal seam without suitable conditions for implementing conventional mining techniques for protective coal seam; and(4) predominant reliance on ‘‘under three" coal resources to ensure production output. This study proposes an integrated, closed-cycle mining-dressing-gas draining-backfilling-mining(MDGBM) technique. The proposed approach involves the mining of protective coal seam, underground dressing of coal and gangue(UDCG), pressure relief and gas drainage before extraction, and backfilling and mining of the protected coal seam. A system for draining gas and mining the protective seam in the rock stratum is designed and implemented based on the geological conditions. This system helps in realizing pressure relief and gas drainage from the protective seam before extraction. Accordingly, another system, which is connected to the existing production system, is established for the UDCG based on the dense medium-shallow trough process. The mixed mining workface is designed to accommodate both solid backfill and conventional fully mechanized coal mining, thereby facilitating coal mining, USCG, and backfilling. The results show that: The mixed mining workface length for the Ji15-31010 protected seam was 220 m with coal production capacity 1.2 million tons per year, while the backfill capacity of gangue was 0.5 million tons per year. The gas pressure decreased from 1.78 to 0.35 MPa, and the total amount of safely mined coal was 1.34 million tons. The process of simultaneously exploiting coal and draining gas was found to be safe, efficient, and green.This process also yielded significant economic benefits.

Fault-Induced Coal Burst Mechanism under Mining-Induced Static and Dynamic Stresses

Fault is a common geological structure that has been revealed in the process of underground coal excavation and mining.The nature of its discontinuous structure controls the deformation,damage,and mechanics of the coal or rock mass.The interaction between this discontinuous structure and mining activities is a key factor that dominates fault reactivation and the coal burst it can induce.This paper first summarizes investigations into the relationships between coal mining layouts and fault occurrences,along with relevant conceptual models for fault reactivation.Subsequently,it proposes mechanisms of fault reactivation and its induced coal burst based on the superposition of static and dynamic stresses,which include two kinds of fault reactivations from:mining-induced quasi-static stress(FRMSS)-dominated and seismic-based dynamic stress(FRSDS)-dominated.These two kinds of fault reactivations are then validated by the results of experimental investigations,numerical modeling,and in situ microseismic monitoring.On this basis,monitoring methods and prevention strategies for fault-induced coal burst are discussed and recommended.The results show that fault-induced coal burst is triggered by the superposition of high static stress in the fault pillar and dynamic stress from fault reactivation.High static stress comes from the interaction of the fault and the roof structure,and dynamic stress can be ascribed to FRMSS and FRSDS.The results in this paper could be of great significance in guiding the monitoring and prevention of fault-induced coal bursts.

Dynamic evolution characteristics of mining-induced fractures in overlying strata detected by radon

<正>For environment protection in mining areas in northwest China,we developed a CTSRM(comprehensive test system by radon measurement) to measure radon radioactivity and detect dynamic evolution characteristics of mining-induced fractures in overlying strata.It was used to simulate the relationship between the dynamic evolution characteristics and radon concentrations of 33201~# coalface at Bulianta coal mine in Inner Mongolia,and feasibility of the method was validate.

Quantitative prediction of mining subsidence and its impact on the environment

This study is focused on the prediction of mining subsidence and its impact on the environment in the Hongqi mining area. The study was carried out by means of a probability integral model based, in first instance based on field surveys and the analysis of data collected from this area. Isolines of mining subsidence were then drawn and the impact caused by mining subsidence on the environment was analyzed quantitatively by spatial analysis with Geographic Information System (GIS). The results indicate that the subsidence area of the first working-mine can be as large as 2.54 km 2 , the maximum subsidence is 3440 mm which will cause 1524 houses to be relocated. The entire subsidence area of the mine can reach 8.09 km 2 , with a maximum subsidence of 3590 mm. Under these circumstances the value of the loss of ecosystem services will reach 5.371 million Yuan and the cost of relocating buildings will increase to 6.858 million Yuan.

Stress and deformation analysis of gob-side pre-backfill driving procedure of longwall mining:a case study

At present,non-pillar entry protection in longwall mining is mainly achieved through either the gob-side entry retaining(GER)procedure or the gob-side entry driving(GED)procedure.The GER procedure leads to difficulties in maintaining the roadway in mining both the previous and current panels.A narrow coal pillar about 5-7 m must be left in the GED procedure;therefore,it causes permanent loss of some coal.The gob-side pre-backfill driving(GPD)procedure effectively removes the wasting of coal resources that exists in the GED procedure and finds an alternative way to handle the roadway maintenance problem that exists in the GER procedure.The FLAC^(3D) software was used to numerically investigate the stress and deformation distributions and failure of the rock mass surrounding the previous and current panel roadways during each stage of the GPD procedure which requires"twice excavation and mining".The results show that the stress distribution is slightly asymmetric around the previous panel roadway after the"primary excavation".The stronger and stiffer backfill compared to the coal turned out to be the main bearing body of the previous panel roadway during the"primary mining".The highest vertical stresses of 32.6 and 23.1 MPa,compared to the in-situ stress of 10.5 MPa,appeared in the backfill wall and coal seam,respectively.After the"primary mining",the peak vertical stress under the coal seam at the floor level was slightly higher(18.1 MPa)than that under the backfill(17.8 MPa).After the"secondary excavation",the peak vertical stress under the coal seam at the floor level was slightly lower(18.7 MPa)than that under the backfill(19.8 MPa);the maximum floor heave and maximum roof sag of the current panel roadway were 252.9 and 322.1 mm,respectively.During the"secondary mining",the stress distribution in the rock mass surrounding the current panel roadway was mainly affected by the superposition of the front abutment pressure from the current panel and the side abutment pressure from the previous panel.The floor heave of the current panel roadway reached a maximum of 321.8 mm at 5 m ahead of the working face;the roof sag increased to 828.4 mm at the working face.The peak abutment pressure appeared alternately in the backfill and the coal seam during the whole procedure of"twice excavation and mining"of the GPD procedure.The backfill provided strong bearing capacity during all stages of the GPD procedure and exhibited reliable support for the roadway.The results provide scientific insight for engineering practice of the GPD procedure.

Design of comprehensive test system for detecting overlying strata mining-induced fractures on surface with radon gas

Based on radon gas properties and its existing projects applications, we firstly attempted to apply geo- physical and chemical properties of radon gas in the field of mining engineering, and imported radioac- tive measurement method to detect the development process of the overlying strata mining-induced fractures and their contained water quality in underground coal mining, which not only innovates a more simple-fast-reliable detection method, but also further expands the applications of radon gas detection technology in mining field. A 3D simulation design of comprehensive testing system for detecting strata mining-induced fractures on surface with radon gas (CTSR) was carried out by using a large-scale 3D solid model design software Pro/Engineer (Pro/E), which overcame three main disadvantages of ''static design thought, 2D planar design and heavy workload for remodification design'' on exiting design for mining engineering test systems. Meanwhile, based on the simulation design results of Pro/E software, the sta- bility of the jack-screw pressure bar for the key component in CTSR was checked with a material mechan- ics theory, which provided a reliable basis for materials selection during the latter machining process.

Quantitative assessment of the spatio-temporal correlations of seismic events induced by longwall coal mining

Rock failure process as a natural response to mining activities is associated with seismic events, which can pose a potential hazard to mine operators, equipment and infrastructures. Mining-induced seismicity has been found to be internally correlated in both time and space domains as a result of rock fracturing during progressive mining activities. Understanding the spatio-temporal(ST) correlation of mininginduced seismic events is an essential step to use seismic data for further analysis, such as rockburst prediction and caving assessment. However, there are no established methods to perform this critical task. Input parameters used for the prediction of seismic hazards, such as the time window of past data and effective prediction distance, are determined based on site-specific experience without statistical or physical reasons to support. Therefore, the accuracy of current seismic prediction methods is largely constrained, which can only be addressed by quantitively assessing the ST correlations of mininginduced seismicity. In this research, the ST correlation of seismic event energy collected from a study mine is quantitatively analysed using various statistical methods, including autocorrelation function(ACF), semivariogram and Moran’s I analysis. In addition, based on the integrated ST correlation assessment, seismic events are further classified into seven clusters, so as to assess the correlations within individual clusters. The correlation of seismic events is found to be quantitatively assessable, and their correlations may vary throughout the mineral extraction process.

Digital image processing-based automatic detection algorithm of cross joint trace and its application in mining roadway excavation practice

This paper proposes a digital image processing-based detection algorithm for cross joint traces of coal roadway heading face.Initially,the acquired images were preprocessed,i.e.,adaptive correction was conducted for non-uniform illumination images based on the 2D gamma function.The edge detection algorithm was then applied to extract the edges of the structural plane,followed by the filtration of the non-structural plane noises.Moreover,the Hough transform algorithm was applied to extract the linear edges;finally,the edges were locally connected in accordance with the angle and distance criteria.The experimental results show that this algorithm can be used to reduce the noise caused by non-uniform illumination and avoid the mutual interference of multi-scale edges,so as to effectively extract the traces of the cross joint.Furthermore,Q-system and rock mass rating(RMR),were applied to conduct a quantitative evaluation on the stand-up time of unsupported roof in the four test images.The Q-system quality scores are 26.7,43.3,3.1,and 6.7,and the RMR quality scores are 56.84,58.73,48.42,and 51.42,respectively.The stand-up time of unsupported roofs with a span of 4.6 m are 30,36,7.7 and 14 d,respectively.

Loose Gangues Backfill Body’s Acoustic Emissions Rules During Compaction Test:Based on Solid Backfill Mining

In fully mechanized solid backfilling mining(FMSBM),the loose gangues backfill body(LGBB)that filled into the goaf becomes the main body of bearing the overburden load.The deformation resistance of LGBB is critical for controlling overburden movement and surface subsidence.During the process of load bearing,LGBB will experience grain crushing,which has a significant effect on its deformation resistance.Gangues block will be accompanied with obvious acoustic emissions(AE)features in process of slipping,flipping and damaging.Under confined compression test,monitoring the AE parameters of LGBB can reveal the impact mechanism of grain crushing on LGBB deformation.The study is of great significance for obtaining an in-depth understanding of the mechanical properties of LGBB,and providing guidance to the engineering practice of FMSBM.In order to study the rules of acoustic emissions(AE)of graded Loose gangues backfill body(LGBB)in confined compression test,this article introduces the AE systems to conventional confined compression test to monitor AE signals resulted from the friction and fragmentation among LGBB.The test results show that in the process of LGBB compaction,AE parameters are highly correlated with the strain-stress curve.AE events of balanced-sized graded gangues are more inactive than other two graded samples in different compression stages,AE events of large-particle-dominated graded gangues are most active.In the spatial distribution,AE events are the most active on the edges and the middle part of test samples and the phenomenon of grain crushing is the most obvious in these positions.

Deformation and Failure Analysis of River Levee Induced by Coal Mining and Its Influence Factor

To study the influence of coal mining on the stability of river levees,a mechanical model of mining-induced river levee deformation was established.This was based on the mining-induced deformation characteristics of river levees and the application of a typical surface subsidence function.Meanwhile,a failure criterion was proposed for river levees.Using some examples,the deformation of,and stress distribution through,river levees under the influence of mining were obtained:the maximum tensile stress on the bottom of the river levee was less than the tensile strength,under which circumstance the river levee remained undamaged.Meanwhile,this research analyzed the influence of three factors including the maximum surface subsidence wmax,half-length of surface subsidence basin L,and foundation coefficient k on the stability of river levees.Results showed that reducing the mining height of the working face and the foundation co-efficient,and increasing the strike length of the working face could reduce the influence of mining on river levees.These results provided a theoretical basis for predicting the mining-induced deformation and failure of river levees.

Numerical simulation of coal wall cutting and lump coal formation in a fully mechanized mining face

It is difficult to accurately calculate the lump coal rate in a fully mechanized mining face.Therefore,a numerical simulation of the coal wall cutting process,which revealed the crack expansion,development,evolution in the coal body and the corresponding lump coal formation mechanism,was performed in PFC2D.Moreover,a correlation was established between the cutting force and lump coal formation,and a statistical analysis method was proposed to determine the lump coal rate.The following conclusions are drawn from the results:(1)Based on a soft ball model,a coal wall cutting model is established.By setting the roller parameters based on linear bonding and simulating the roller cutting process of the coal body,the coal wall cutting process is effectively simulated,and accurate lump coal rate statistics are provided.(2)Under the cutting stress,the coal body in the working face underwent three stages—microfracture generation,fracture expansion,and fracture penetration—to form lump coal,in which the fracture direction is orthogonal to the cutting pressure chain.Within a certain range from the roller,as the cutting depth of the roller increased,the number of new fractures in the coal body first increases and then stabilizes.(3)Under the cutting stress,the fractured coal body is locally compressed,thereby forming a compact core.The formation and destruction of the compact core causes fluctuations in the cutting force.The fluctuation amplitude is positively related to the coal mass.(4)Because the simulation does not consider secondary damage in the coal,the simulated lump coal rate is larger than the actual lump coal rate in the working face;this deviation is mainly concentrated in large lump coal with a diameter greater than 300 mm.

Experimental study of coal flow characteristics under mining disturbance in China

With annually increased coal mining depth,gas extraction becomes more and more problematic.The gas extraction efect depends on coal seam permeability,which,in turn,is afected by many factors,including loading and unloading stresses and strains in the coal seam.Stresses induce internal cracks,resulting in cleats and gas emission channels,the coal seam permeability permanently changes accordingly.To clarify the stress-induced efects on coal seam permeability,this survey summarized the available approaches used to link the stress path and seepage law in the coal body seepage law,which can be classifed into two design methods:single load variation and combined feld mining method.The characterization methods used to observe the surface of coal samples and three-dimensional reconstruction include electron microscopy,CT scanning,and Nuclear Magnetic Resonance(NMR).According to the stress paths designed by the above two approaches,the seepage laws and similarities of three kinds of coal samples with the fractured structure were summarized in this paper.The following directions are recommended to study the seepage law of coal bodies with three kinds of fractured structures under stress.Firstly,the stress path of the experimental coal body should be designed by the combined feld mining method.The stressed environment of a deep coal seam is complicated,and the axial and confning pressures change simultaneously.Therefore,one cannot fully refect the real situation on-site by studying permeability evolution alone.Secondly,during the coal seam mining,the stressed state changes from time to time,and the development of coal seam fractures is afected by mining.When studying the stress efect on seepage of coal samples,the fractured structure of coal samples should be considered.Finally,the available structural characterization methods of coal samples can be combined with the 3D printing technology,which would produce artifcial samples with the fractured structure characteristics of natural coal.