VR-based digital twin for remote monitoring of mining equipment:Architecture and a case study

Background Traditional methods for monitoring mining equipment rely primarily on visual inspections,which are time-consuming,inefficient,and hazardous.This article introduces a novel approach to monitoring mission-critical systems and services in the mining industry by integrating virtual reality(VR)and digital twin(DT)technologies.VR-based DTs enable remote equipment monitoring,advanced analysis of machine health,enhanced visualization,and improved decision making.Methods This article presents an architecture for VR-based DT development,including the developmental stages,activities,and stakeholders involved.A case study on the condition monitoring of a conveyor belt using real-time synthetic vibration sensor data was conducted using the proposed methodology.The study demonstrated the application of the methodology in remote monitoring and identified the need for further development for implementation in active mining operations.The article also discusses interdisciplinarity,choice of tools,computational resources,time and cost,human involvement,user acceptance,frequency of inspection,multiuser environment,potential risks,and applications beyond the mining industry.Results The findings of this study provide a foundation for future research in the domain of VR-based DTs for remote equipment monitoring and a novel application area for VR in mining.

A systematic review of occupational exposure to respirable coal mine dust(RCMD)in the U.S.mining industry

Cumulative inhalation of respirable coal mine dust(RCMD)can lead to severe lung diseases,including coal worker's pneu-moconiosis(CWP),silicosis,mixed dust pneumoconiosis,dust-related diffuse fibrosis(DDF),and progressive massive fibrosis(PMF).Statistics from the number of reported cases showed a significant decrease in the progression of respiratory diseases in the 1990s.However,an unexpected increase in the number of CWP cases was reported in the late 1990s.To date,there has been no comprehensive systematic review to assess all contributing factors to the resurgence of CWP cases.This study aims to investigate the effects of various mining parameters on the prevalence of CWP in coal mines.A systematic review using the preferred reporting items for systematic reviews and meta-analysis(PRISMA)method was conducted to investigate the health effects of RCMD exposure and identify the factors that may contribute to the recent resurgence of CWP cases.The systematic review yielded a total of 401 papers,which were added to the database.The total number of 148 and 208 papers were excluded from the database in the process of screening and eligibility,respectively.Then,18 papers were considered for data selection and full-text assessment.The review revealed that factors including geographic location,mine size,mining operation type,coal-seam thickness,coal rank,changes in mining practices,technology advancement,and engi-neering dust control practices are contributing to the recent resurgence of CWP among coal workers.However,the evidence for root causes is limited owing to the methodological constraints of the studies;therefore,further detailed studies are needed.

Tomographic reconstruction of the Earth’s magnetosheath from multiple spacecraft:a theoretical study

Following our earlier work on tomographic reconstruction of the magnetosheath soft X-ray emissions with superposed epoch analysis of many images recorded from a single spacecraft we now explore the instantaneous reconstruction of the magnetosheath and magnetopause using a few images recorded simultaneously from a few spacecraft.This work is motivated by the prospect of possibly having two or three soft X-ray imagers in space in the coming years,and that many phenomena which occur at the magnetopause boundary,such as reconnection events and pressure pulse responses,do not lend themselves as well to superposed epoch analysis.If the reconstruction is successful-which we demonstrate in this paper that it can be-this collection of imagers can be used to reconstruct the magnetosheath and magnetopause from a single image from each spacecraft,allowing for high time resolution reconstructions.In this paper we explore the reconstruction using,two,three,and four spacecraft.We show that the location of the subsolar point of the magnetopause can be determined with just two satellites,and that volume emissions of soft X-rays,and the shape of the boundary,can be reconstructed using three or more satellites.

Micro-macro evolution of mechanical behaviors of thermally damaged rock:A state-of-the-art review

The influence of thermal damage on macroscopic and microscopic characteristics of different rocks has received much attention in the field of rock engineering.When the rocks are subjected to thermal treatment,the change of macroscopic characteristics and evolution of micro-structure would be induced,ultimately resulting in different degrees of thermal damage in rocks.To better understand the thermal damage mechanism of different rocks and its effect on the rock performance,this study reviews a large number of test results of rock specimens experiencing heating and cooling treatment in the laboratory.Firstly,the variations of macroscopic behaviors,including physical parameters,mechanical parameters,thermal conductivity and permeability,are examined.The variations of mechanical parameters with thermal treatment variables(i.e.temperature or the number of thermal cycles)are divided into four types.Secondly,several measuring methods for microstructure,such as polarizing microscopy,fluorescent method,scanning electron microscopy(SEM),X-ray computerized tomography(CT),acoustic emission(AE)and ultrasonic technique,are introduced.Furthermore,the effect of thermal damage on the mechanical parameters of rocks in response to different thermal treatments,involving temperature magnitude,cooling method and thermal cycle,are discussed.Finally,the limitations and prospects for the research of rock thermal damage are proposed.

Prediction of rock fragmentation in a fiery seam of an open-pit coal mine in India

Spontaneous combustion of coal increases the temperature in adjoining overburden strata of coal seams and poses a challenge when loading blastholes.This condition,known as hot-hole blasting,is dangerous due to the increased possibility of premature explosions in loaded blastholes.Thus,it is crucial to load the blastholes with an appropriate amount of explosives within a short period to avoid premature detonation caused by high temperatures of blastholes.Additionally,it will help achieve the desired fragment size.This study tried to ascertain the most influencial variables of mean fragment size and their optimum values adopted for blasting in a fiery seam.Data on blast design,rock mass,and fragmentation of 100 blasts in fiery seams of a coal mine were collected and used to develop mean fragmentation prediction models using soft computational techniques.The coefficient of determination(R^(2)),root mean square error(RMSE),mean absolute error(MAE),mean square error(MSE),variance account for(VAF)and coefficient of efficiency in percentage(CE)were calculated to validate the results.It indicates that the random forest algorithm(RFA)outperforms the artificial neural network(ANN),response surface method(RSM),and decision tree(DT).The values of R^(2),RMSE,MAE,MSE,VAF,and CE for RFA are 0.94,0.034,0.027,0.001,93.58,and 93.01,respectively.Multiple parametric sensitivity analyses(MPSAs)of the input variables showed that the Schmidt hammer rebound number and spacing-to-burden ratio are the most influencial variables for the blast fragment size.The analysis was finally used to define the best blast design variables to achieve optimum fragment size from blasting.The optimum factor values for RFA of S/B,ld/B and ls/ld are 1.03,1.85 and 0.7,respectively.

Diagenetic evolution and reservoir quality of the Oligocene sandstones in the Baiyun Sag, Pearl River Mouth Basin, South China Sea

The Oligocene Zhuhai sandstones are significant reservoirs for hydrocarbons in the Baiyun Sag, South China Sea.For effective appraisal, exploration and exploitation of such a deep-water hydrocarbon sandstone, samples of five wells from depths of 850 m to 3 000 m were studied. A series of comprehensive petrographic and geochemical analyses were performed to unravel the diagenetic features and their impact on the reservoir quality.Petrographically, the sandstones are dominated by feldspathic litharenites and lithic arenites with fine to medium grain sizes and moderate to good sorting. The reservoir quality varies greatly with a range of porosity from 0.2% to 36.1% and permeability from 0.016 ×10~(–3) μm~2 to 4 301 ×10~(–3) μm~2, which is attributed to complex diagenetic evolution related to sedimentary facies;these include compaction, cementation of calcite, dolomite, siderite and framboidal pyrite in eogenetic stage;further compaction, feldspar dissolution, precipitation of ferrocalcite and ankerite, quartz cements, formation of kaolinite and its illitization, precipitation of albite and nodular pyrite, as well as hydrocarbon charge in mesogenetic stage. The dissolution of feldspar and illitization of kaolinite provide internal sources for the precipitation of quartz cement, while carbonate cements are derived from external sources related to interbedded mudstones and deep fluid. Compaction is the predominant factor in reducing the total porosity, followed by carbonate cementation that leads to strong heterogeneity. Feldspar dissolution and concomitant quartz and clay cementation barely changes the porosity but significantly reduces the permeability.The high-quality reservoirs can be concluded as medium-grained sandstones lying in the central parts of thick underwater distributary channel sandbodies(>2 m) with a high content of detrital quartz but low cement.

Enhancing XRF sensor-based sorting of porphyritic copper ore using particle swarm optimization-support vector machine(PSO-SVM)algorithm

X-ray fluorescence(XRF)sensor-based ore sorting enables efficient beneficiation of heterogeneous ores,while intraparticle heterogeneity can cause significant grade detection errors,leading to misclassifications and hindering widespread technology adoption.Accurate classification models are crucial to determine if actual grade exceeds the sorting threshold using localized XRF signals.Previous studies mainly used linear regression(LR)algorithms including simple linear regression(SLR),multivariable linear regression(MLR),and multivariable linear regression with interaction(MLRI)but often fell short attaining satisfactory results.This study employed the particle swarm optimization support vector machine(PSO-SVM)algorithm for sorting porphyritic copper ore pebble.Lab-scale results showed PSO-SVM out-performed LR and raw data(RD)models and the significant interaction effects among input features was observed.Despite poor input data quality,PSO-SVM demonstrated exceptional capabilities.Lab-scale sorting achieved 93.0%accuracy,0.24%grade increase,84.94%recovery rate,57.02%discard rate,and a remarkable 39.62 yuan/t net smelter return(NSR)increase compared to no sorting.These improvements were achieved by the PSO-SVM model with optimized input combinations and highest data quality(T=10,T is XRF testing times).The unsuitability of LR methods for XRF sensor-based sorting of investigated sample is illustrated.Input element selection and mineral association analysis elucidate element importance and influence mechanisms.

Mechanical behavior and damage constitutive model of sandstone under hydro-mechanical (H-M) coupling

Underground engineering often passes through water-rich fractured rock masses, which are prone to fracture and instability under the long-term coupling of in-situ stress field and pore water(P-W) pressure, ultimately threatening the stability of underground structures. In order to explore the mechanical properties of rocks under H-M coupling, the corresponding damage constitutive(D-C) model has become the focus of attention. Considering the inadequacy of the current research on rock strength parameters,energy evolution characteristics and D-C model under H-M coupling, the mechanical properties of typical sandstone samples are discussed based on laboratory tests. The results show that the variation of characteristic stresses of sandstone under H-M coupling conforms to the normalized attenuation equation and Mohr-Coulomb(M-C) criterion. The P-W pressure mechanism of sandstone exhibits a dynamic change from softening effect to H-M fracturing effect. The closure stress is mainly provided by cohesive strength, while the initiation stress, damage stress, and peak stress are jointly dominated by cohesive strength and friction strength. In addition, residual stress is attributed to the friction strength formed by the bite of the fracture surface. Subsequently, the energy evolution characteristics of sandstone under H-M coupling were studied, and it was found that P-W pressure weakened the energy storage capacity and energy dissipation capacity of sandstone, and H-M fracturing was an important factor in reducing its energy storage efficiency. Finally, combined with energy dissipation theory and statistical damage theory, two types of D-C models considering P-W pressure are proposed accordingly, and the model parameters can be determined by four methods. The application results indicate that the proposed and modified D-C models have high reliability, and can characterize the mechanical behavior of sandstone under H-M coupling, overcome the inconvenience of existing D-C models due to excessive mechanical parameters,and can be applied to the full-range stress–strain process. The results are conducive to revealing the deformation and damage mechanisms of rocks under H-M coupling, and can provide theoretical guidance for related engineering problems.

Effect of slope angle on fractured rock masses under combined influence of variable rainfall infiltration and excavation unloading

Intense precipitation infiltration and intricate excavation processes are crucial factors that impact the stability and security of towering and steep rock slopes within mining sites.The primary aim of this research was to investigate the progression of cumulative failure within a cracked rock formation,considering the combined effects of precipitation and excavation activities.The study was conducted in the Huangniuqian eastern mining area of the Dexing Copper Mine in Jiangxi Province,China.An engineering geological investigation was conducted,a physical model experiment was performed,numerical calculations and theoretical analysis were conducted using the matrix discrete element method(Mat-DEM),and the deformation characteristics and the effect of the slope angle of a fractured rock mass under different scenarios were examined.The failure and instability mechanisms of the fractured rock mass under three slope angle models were analyzed.The experimental results indicate that as the slope angle increases,the combined effect of rainfall infiltration and excavation unloading is reduced.A novel approach to simulating unsaturated seepage in a rock mass,based on the van Genuchten model(VGM),has been developed.Compared to the vertical displacement observed in a similar physical experiment,the average relative errors associated with the slope angles of 45,50,and 55were 2.094%,1.916%,and 2.328%,respectively.Accordingly,the combined effect of rainfall and excavation was determined using the proposed method.Moreover,the accuracy of the numerical simulation was validated.The findings contribute to the seepage field in a meaningful way,offering insight that can inform and enhance existing methods and theories for research on the underlying mechanism of ultra-high and steep rock slope instability,which can inform the development of more effective risk management strategies.

Efficient desorption and reuse of collector from the flotation concentrate:A case study of scheelite

Flotation is the most common method to obtain concentrate through the selective adsorption of collectors on target minerals to make them hydrophobic and floatable.In the hydrometallurgy of concentrate,collectors adsorbed on concentrate can damage ion-exchange resin and increase the chemical oxygen demand(COD)value of wastewater.In this work,we proposed a new scheme,i.e.,desorbing the collectors from concentrate in ore dressing plant and reusing them in flotation flowsheet.Lead nitrate and benzohydroxamic acid(Pb-BHA)complex is a common collector in scheelite flotation.In this study,different physical(stirring or ultrasonic waves)and chemical(strong acid or alkali environment)methods for facilitating the desorption of Pb-BHA collector from scheelite concentrate were explored.Single-mineral desorption tests showed that under the condition of pulp pH 13 and ultrasonic treatment for 15 min,the highest desorption rates of Pb and BHA from the scheelite concentrate were 90.48%and 63.75%,respectively.Run-of-mine ore flotation tests revealed that the reuse of desorbed Pb and BHA reduced the collector dosage by 30%for BHA and 25%for Pb.The strong alkali environment broke the chemical bonds between Pb and BHA.The cavitation effect of ultrasonic waves effectively reduced the interaction intensity between Pb-BHA collector and scheelite surfaces.This method combining ultrasonic waves and strong alkali environment can effectively desorb the collectors from concentrate and provide“clean”scheelite concentrate for metallurgic plants;the reuse of desorbed collector in flotation flowsheet can reduce reagent cost for ore dressing plants.

Deformation mechanism and roof pre-splitting control technology of gob-side entry in thick hard main roof full-mechanized longwall caving panel

This paper explores the deformation mechanism and control technology of roof pre-splitting for gob-side entries in hard roof full-mechanized longwall caving panel(LTCC).The investigation utilizes a comprehensive approach that integrates field monitoring,theoretical analysis,and numerical simulation.Theoretical analysis has illuminated the influence of the length of the lateral cantilever beam of the main roof(LCBM)above the roadway on the stability of the gob-side entry behind the panel.Numerical simulations have further revealed that the longer LCBM results in heightened vertical stress within the coal pillar,developed cracks around the roadway,and more pronounced damage to the roadway.Moreover,numerical simulations also demonstrate the potential of roof pre-splitting technology in optimizing the fracture position of the hard roof.This technology significantly reduces the length of the LCBM,thereby alleviating stress concentration in the coal pillars and integrated coal rib while minimizing the destruction of the gob-side entry.Therefore,this manuscript first proposes the use of roof pre-splitting technology to control roadway deformation,and automatically retain the entry within a hard roof LTCC panel.Field implementation has demonstrated that the proposed automatically retained entry by roof pre-splitting technology effectively reduces gob-side entry deformation and achieves automatically retained entry.

Change in Grain-Size Composition of Lignite under Cyclic Freezing-Thawing and Wetting-Drying

The paper presents the change in grain-size composition of lignite under cyclic freezing-thawing (FTC) and wetting-drying (WDC). The article shows that in the spring and autumn periods the lignites can be subjected to repeated freezing-thawing and wetting-drying, which determines the possibility of changing their grain-size composition and structure. Experimental studies in laboratory conditions on the influence of cyclic freezing-thawing (FTC) and wetting-drying (WDC) on the quality indicators of lignites have been carried out, their granulometric (fractional) composition has been studied. Freezing-thawing cycle conditions are as follows (FTC): minimum exposure temperature: -20°C;maximum: +5°C;relative humidity: 30%;number of processing cycles: 3. Wetting-drying cycles are as follows (WDC): drying temperatures are +20, +40, +60, +80°C, drying time 90 minutes, the coals are further subjected to rain (soaking) for a period of water saturation to humidity of 30% - 40% and dry again. The number of wetting-drying cycles is 3 times. The tests have revealed the destructive effects of FTC and WDC on the samples of lower metamorphic grade coal, and the cycles of wet-dry lead to the much higher yield of fine sizes (-6+0;-13+0 mm) than the cycles of freeze-thaw. Furthermore, it is found that the increase in the yield of fines depends on the heating temperature: coal disintegration proceeds more intensively at a higher temperature of drying.

Theoretical design and dynamic simulation of new mining paths of tracked miner on deep seafloor

With comprehensive considerations of the operational safety and collection efficiency for the tracked miner collecting the seafloor poly-metallic nodules, two new improved mining paths for the miner on the deep seafloor were proposed. Compared to the conventional mining path, the design principles and superiorities of the two new paths are that the miner turning with relative long radius should avoid large sinkage and high slip, so as to ensure its operational safety, while the space between its straight-line trajectories before and after the turning is optimum, which is designed as the total width of the miner, and collect nodules as more as possible, so as to ensure its collection efficiency. To realize the new mining paths, theoretical designs and quantitative calculations were carried out to determine the exact positions for the speed controls of the miner during its whole operation process. With the new dynamic model of the miner, and through regulations of the speeds of the left and right tracks of the miner on the exact motion positions according to the theoretical calculations, the two new improved mining paths for the miner on the seafloor were successfully simulated, thus the turning radius of the miner in the simulation is about 21.8 m, while the distance between the straight-line trajectories before and after the turning is about 5.2 m. The dynamic simulation results preliminarily prove the feasibility of these two new mining paths, and further can provide important theoretical guidance and useful technical reference for the practical tracked miner operation and control on the seafloor.

Similar material simulation of time series system for induced caving of roof in continuous mining under backfill

With the help of similar material simulation test,time series system for induced caving of roof in continuous mining under complex backfill in ore body No.92 of Tongkeng Tin Mine was studied. According to the similarity theory,a two-dimensional similar simulation test-bed was constructed. The stress and displacement that change along with the advance of mining were acquired and analyzed automatically by data system. The processes of continuous mining of ore-block in 5 intervals and artificial induced caving of roof were simulated. The results of the test show that ore body remained as safety roof in thickness of 15 m guarantees the safe advance of stoping work face. Caving of safety roof puts in practice at the first two mining intervals when the third interval of continuous mining is finished,and one interval as the safety distance should be kept all the time between stopping and caving. While mining in the last interval,pre-slotting should be implemented first of all,and the roof of the last two mining intervals is caved simultaneously. Only this kind of time series system can be an efficient and safe way for induced caving of roof in continuous mining.

Analysis on Shock Wave Speed of Water Hammer of Lifting Pipes for Deep-Sea Mining

Water hammer occurs whenever the fluid velocity in vertical lifting pipe systems for deep-sea mining suddenly changes. In this work, the shock wave was proven to play an important role in changing pressures and periods, and mathematical and numerical modeling technology was presented for simulated transient pressure in the abnormal pump operation. As volume concentrations were taken into account of shock wave speed, the experiment results about the pressure-time history, discharge-time history and period for the lifting pipe system showed that： as its concentrations rose up, the maximum transient pressure went down, so did its discharges; when its volume concentrations increased gradually, the period numbers of pressure decay were getting less and less, and the corresponding shock wave speed decreased. These results have highly coincided with simulation results. The conclusions are important to design lifting transporting system to prevent water hammer in order to avoid potentially devastating consequences, such as damage to components and equipment and risks to personnel.

An Integrated Dynamic Model of Ocean Mining System and Fast Simulation of Its Longitudinal Reciprocating Motion

An integrated dynamic model of China＇s deep ocean mining system is developed and the fast simulation analysis of its longitudinal reciprocating motion operation processes is achieved. The seafloor tracked miner is built as a three-dimensional single-body model with six-degree-of-freedom. The track-terrain interaction is modeled by partitioning the track-terrain interface into a certain number of mesh elements with three mutually perpendicular forces, including the normal force, the longitudinal shear force and the lateral shear force, acting on the center point of each mesh element. The hydrodynamic force of the miner is considered and applied. By considering the operational safety and collection efficiency, two new mining paths for the miner on the seafloor are proposed, which can be simulated with the established single-body dynamic model of the miner. The pipeline subsystem is built as a three-dimensional multi-body discrete element model, which is divided into rigid elements linked by flexible connectors. The flexible connector without mass is represented by six spring-damper elements. The external hydrodynamic forces of the ocean current from the longitudinal and lateral directions are both considered and modeled based on the Morison formula and applied to the mass center of each corresponding discrete rigid element. The mining ship is simplified and represented by a general kinematic point, whose heave motion induced by the ocean waves and the longitudinal and lateral towing motions are considered and applied. By integrating the single-body dynamic model of the miner and the multi-body discrete element dynamic model of the pipeline, and defining the kinematic equations of the mining ship, the integrated dynamic model of the total deep ocean mining system is formed. The longitudinal reciprocating motion operation modes of the total mining system, which combine the active straight-line and turning motions of the miner and the ship, and the passive towed motions of the pipeline, are proposed and simulated with the developed 3D dynamic model. Some critical simulation results are obtained and analyzed, such as the motion trajectories of key subsystems, the velocities of the buoyancy modules and the interaction forces between subsystems, which in a way can provide important theoretical basis and useful technical reference for the practical deep ocean mining system analysis, operation and control.

Design and implementation of model test installation of heave compensation system of deepsea mining

In order to validate the simulation model and develop heave compensation control strategy,heave compensation model tests were performed.The model test installation includes themining ship motion simulator,the heave compensation system,the lifting pipe simulator,the buffer simulator and the water pool.The tests ofmining ship motion simulator show that it is able to perform under the predetermined attitude path smoothly and can meet the requirements of themining ship motions.The heave compensation effect is more than 60% under random wave and the goal is set to be 50%.The model test results indicate that this heave compensation system is effective and feasible.

A probe into“mining technique in the condition of floor failure”for coal seam above longwall goafs

Targeting at the coal seam with useful value discarded above goafs,attempted to explore the feasibility of'mining technique in the condition of floor failure' from theoretical point of view,and predicted.It indicated that mining technique in the condition of floor failure used above Longwall Goafs in Baijiazhuang Mining is totally feasible.At law,the deformation of the floor in the mining technique by means of probability-integral method.And it is discov- ered that deformed basin can emerge in the footwall of No.6 coal seam and its maximum subsidence was possibly 1 633 mm or so and its maximum positive curvature is 61.74/10^(-3). At last,it therefore suggests appropriate ground pressure control measures as strengthening observation of ground pressure and adopting false slope for exploitation and strengthening support for reasonable push and slide based on the adverse ground pressure behaviors possibly occurring in the mining technique.This serves to gather data and lay sturdy founda- tion for further probe into the mining technique,and offers theoretical and technical grounds for concrete implementation of the mining technique.

Underground mining of thick coal seams

This paper reviews underground mining methods for total thickness of a thick coal seam in single lift （TI＇rCSSL）. Review shows the required engineering for extraction of thick seams needs to be fitted with thickness of the seam, behavior of rock-mass and surrounding stress conditions for efficient mining. Variants of TI＇rCSSL are able to extract a maximum 10-12 m thickness only. An improvement in bending moment of the overlying coal band in longwall top coal caving （LTCC） provides better under-winning opportunity for the roof coal band. An acceptable limit of 25 MPa compressive strength of coal for the success of LTCC may be increased under favorable geo-technical conditions. Bord and pillar in India adopted induced caving of roof coal band for single lift depillaring of total thickness （SLDTr） of a compe- tent thick coal seam developed along floor. Case studies are given to arrest the adverse effects of extrac- tion height on pillars.

An analysis of contributing mining factors in coal workers’pneumoconiosis prevalence in the United States coal mines,1986-2018

In the United States,an unexpected and severe increase in coal miners’lung diseases in the late 1990s prompted researchers to investigate the causes of the disease resurgence.This study aims to scrutinize the effects of various mining parameters,including coal rank,mine size,mine operation type,coal seam height,and geographical location on the prevalence of coal worker’s pneumoconiosis(CWP)in surface and underground coal mines.A comprehensive dataset was created using the U.S.Mine Safety and Health Administration(MSHA)Employment and Accident/Injury databases.The information was merged based on the mine ID by utilizing SQL data management software.A total number of 123,589 mine-year observations were included in the statistical analysis.Generalized Estimating Equation(GEE)model was used to conduct a statistical analysis on a total of 29,707,and 32,643 mine-year observations for underground and surface coal mines,respectively.The results of the econometrics approach revealed that coal workers in underground coal mines are at a greater risk of CWP comparing to those of surface coal operations.Furthermore,underground coal mines in the Appalachia and Interior regions are at a higher risk of CWP prevalence than the Western region.Surface coal mines in the Appalachian coal region are more likely to CWP development than miners in the Western region.The analysis also indicated that coal workers working in smaller mines are more vulnerable to CWP than those in large mine sizes.Furthermore,coal workers in thin-seam underground mine operations are more likely to develop CWP.