A Tightly Coupled GIS and Spatiotemporal Modeling for Methane Emission Simulation in the Underground Coal Mine System

: Mine safety is of primary concern in the underground coal mining system. At present, there is a lack of an e ﬃ cient platform to manage the numerical simulation procedure and inherent spatiotemporal data for coal mine disasters. This necessitates the coupling of spatiotemporal model with geographic information system (GIS) in practical application. Here, a novel spatiotemporal model tightly coupled with GIS is presented to improve the model-data integration. Such tight coupling is achieved by developing a lattice Boltzmann method (LBM) based turbulent model with an underlying shared FluentEntity model within the LongRuanGIS platform. The case study and comparison with the traditional computational ﬂuid dynamics (CFD) method demonstrated that the platform is capable and e ﬀ ective in providing functionalities for lattice domain decomposition, simulation, visualization and analyses, as well as improving the computational e ﬃ ciency. The proposed approach and platform, promising for the disaster prevention, o ﬀ er a template for future GIS-Model integration and also applicable for other underground coal mine disasters.


Introduction
Mine safety issue has long been a paramount concern in the underground coal mining business, since mine accidents can lead to serious injuries for personnel, substantial economic losses, and delayed production. Numerical modelling offers a useful tool for prediction and control of various disasters in mine workings. The principles of computational fluid dynamics are widely applied to the underground coal mine systems and numerous computational fluid dynamics (CFD) methods are utilized to simulate various ventilation-related safety and health issues [1][2][3][4]. It is remarkably noted that the procedure of CFD simulation interrelated with geometric modeling and numerical calculation, is inherently geospatial. The spatiotemporally simulated data requires effective management and analysis, which is essential to reduce the risk associated with mine accidents. However, the current simulations of coal mine disasters mostly depend on the independent third-party numerical simulation software, and there is no professional simulation platform developed for specific problems. For instance, numerical platforms such as Fluent, COSFLOW, FLAC3D, and AutoReaGas are utilized by many researchers in simulating gas emission in mining face, the hole wall, and other ventilation related activities [5][6][7][8]. There is a lack of an integrative and effective platform to store and manage the massive amount of spatiotemporal simulated data, as well as the data visualization and analysis, which makes the whole analysis procedure difficult to direct the on-site production. According to researches proposed by Goodchild, one of the key challenges in the application of physics based models is large amount of location-based information by spatiotemporal query and analysis to guide the coal mine operation.
Longruan GIS (LrGIS) developed by our team, is professionally designed and ideally suited to assist mining professionals in meeting the complex challenges of running the mine operation. With tools to compile, process, display, analyze, and archive massive volumes of data, the LrGIS platform is increasingly applied in the business of mining industry in China, which allows us to develop the mathematical model and manage volumes of data through its secondary development mode [38]. Thus, this study chooses LrGIS as the supporting GIS software package, and it takes the methane gas emission as the case study, which is one of the most crucially important mine safety issues. It is worth noted that the numerical model can be extended and improved to simulate other mining activities, such as ventilation solution, mine fire, and methane explosion.
We performed the model coupling and numerical simulation from the following ways. First, the U-type ventilation system-15116 in Sijiazhuang coal mine is specified as the computational domain, and a novel turbulent LBM based velocity-concentration coupling model is proposed to investigate the spatiotemporal characteristics of the airflow behavior and the methane gas emission in a coal mine working face. Second, the numerical model is integrated with LrGIS platform via FluentEntity model, and the graphical representations of simulated field data as well as various map layers are demonstrated. From this, the spatial characteristics of methane distribution can be intuitively visualized, and detailed location-based information for the practical mining production can be provided by spatiotemporal query and analysis. Third, the numerical simulation results are verified by field measurements in Sijiazhuang coal mine and compared with traditional CFD methods, demonstrating the accuracy and reliability of the proposed model. Thus, the developed system can assist decision makers to simulate, understand, and forecast ventilation-related activities and its dynamic change inside the laneway, which is expected to be an efficient tool for improving and maintaining the health standards in the underground coal mining industry.

Case Study
A fully mechanized U-type working face (15116 laneway) in Sijiazhuang coal mine was selected as the physical prototype to implement the GIS-coupled spatiotemporal modeling, which is located in Shanxi province, China. This mine working face is selected for access to precisely and continuously recorded mine data. The layout of the working face laneway is shown in Figure 1. It has a working face zone, wind intake laneway, and wind outtake laneway with the rectangular cross-section of 4.2 m wide and 4 m high. The working face is 220 m long, while the wind intake and outtake laneways are all 600 m long.
According to the safety operation regulations in Sijiazhuang coal mine, the airflow quantity in the 15116 working face laneway is 53.3 m 3 /s, which is calculated based on the requirement of the mine gas and carbon dioxide emission quantity, the temperature and the number of maximum miners in the working face area, and it has been checked by the maximum and minimum wind speed. The estimated methane emission flow rate is about 0.2 m 3 /s. During the procedure of coal mining in the U-type laneway system, the working face is the main source of the methane gas emission. A large amount of the methane gas generated from the working face disperses with the pressing air from the wind inflow side, finally, part of it is removed through the wind inflow laneway, while the rest stays in the laneway space. During the delivering process, movement of methane coupled with the airflow exhibits a species transportation flow. Here, we have three model assumptions consistent with the published research [39]: (i) the airflow in the working face zone belongs to the turbulent flow; (ii) the airflow and methane gas are continuous medium gas; (iii) heat and mass transfer are ignored during the process of air flow.

Model-GIS Coupling Strategy
In this study, the numerical model and LrGIS platform are linked to achieve the model-GIS integration. Specifically, the proposed numerical model is the velocity-concentration coupled 2D LBM, which is integrated with LrGIS platform by developing the FluntEntity model, so that both GIS platform and model functions can access the geospatial simulation results. The GIS platform, the numerical model, and the integration strategies are described in detail in the following sections.

LrGIS Platform
LrGIS (Longruan GIS 3.2) is a professional coal mine GIS platform (Beijing LongRuan Technologies Inc.) with its own intellectual property right, which is widely applied in the business of mining industry in China and well supported by professional developers. It is well designed for coal mine system with respect to the geospatial features of coal mine laneway, ventilation network, and etc., and professionally focuses on large amount of mine spatial data, such as surveying, geology, hydrology reserves, transportation, ventilation network, mine design schedule. Besides, it fulfills the integrated management of massive spatial data and information sharing among different platforms, and provides flexible secondary development mode, such as API, DLL, ActiveX, and application template [40]. Therefore, LrGIS is suitable to perform GIS-numerical model integration in this study.
Despite all these advantages, these functions rely on inherent and advanced GIS capabilities, and there is a lack of defined conceptual model associated with spatiotemporal modeling to simulate various kinds of disasters to satisfy some actual needs, such as the methane gas emission simulation, the gas explosion prediction in emergency response, ventilation network design based on the presimulation, and mine fire simulation to guide miners to escape from the disaster and etc.. The integration of the spatiotemporal modeling with GIS can provide great support for practical production needs. These inbuilt GIS functions are insufficient in terms of spatiotemporal modeling, because the data require specific treatments to be used in spatiotemporal modeling. In this study, explicit mesh creation is not needed, which is usually troublesome in traditional CFD method and inconsistencies of input geometries always lead to labor intensive and time-consuming process for mesh generation [41,42]. The integration mode proposed in this study avoids this problem and improves the computational efficiency.

Mathematical Models
The essential regularity of methane dispersion and its coupling with air flow is a key issue for the methane gas emission problem, and its flowing regulation belongs to viscous Newtonian fluid, which is governed by the Navier-Stokes (N-S) equations. There are many numerical approaches to solve the N-S equations, such as the traditional CFD methods and the LBM. Here the LBM based numerical method is adopted. In phase space, the Boltzmann equation is discretized in threefold,

Model-GIS Coupling Strategy
In this study, the numerical model and LrGIS platform are linked to achieve the model-GIS integration. Specifically, the proposed numerical model is the velocity-concentration coupled 2D LBM, which is integrated with LrGIS platform by developing the FluntEntity model, so that both GIS platform and model functions can access the geospatial simulation results. The GIS platform, the numerical model, and the integration strategies are described in detail in the following sections.

LrGIS Platform
LrGIS (Longruan GIS 3.2) is a professional coal mine GIS platform (Beijing LongRuan Technologies Inc.) with its own intellectual property right, which is widely applied in the business of mining industry in China and well supported by professional developers. It is well designed for coal mine system with respect to the geospatial features of coal mine laneway, ventilation network, and etc., and professionally focuses on large amount of mine spatial data, such as surveying, geology, hydrology reserves, transportation, ventilation network, mine design schedule. Besides, it fulfills the integrated management of massive spatial data and information sharing among different platforms, and provides flexible secondary development mode, such as API, DLL, ActiveX, and application template [40]. Therefore, LrGIS is suitable to perform GIS-numerical model integration in this study.
Despite all these advantages, these functions rely on inherent and advanced GIS capabilities, and there is a lack of defined conceptual model associated with spatiotemporal modeling to simulate various kinds of disasters to satisfy some actual needs, such as the methane gas emission simulation, the gas explosion prediction in emergency response, ventilation network design based on the pre-simulation, and mine fire simulation to guide miners to escape from the disaster and etc. The integration of the spatiotemporal modeling with GIS can provide great support for practical production needs. These inbuilt GIS functions are insufficient in terms of spatiotemporal modeling, because the data require specific treatments to be used in spatiotemporal modeling. In this study, explicit mesh creation is not needed, which is usually troublesome in traditional CFD method and inconsistencies of input geometries always lead to labor intensive and time-consuming process for mesh generation [41,42]. The integration mode proposed in this study avoids this problem and improves the computational efficiency.

Mathematical Models
The essential regularity of methane dispersion and its coupling with air flow is a key issue for the methane gas emission problem, and its flowing regulation belongs to viscous Newtonian fluid, which is governed by the Navier-Stokes (N-S) equations. There are many numerical approaches to solve the N-S equations, such as the traditional CFD methods and the LBM. Here the LBM based numerical method is adopted. In phase space, the Boltzmann equation is discretized in threefold, including space, time, and velocities. The movement and distributions of the fluid are described as particle distribution functions residing at the sites of a regular lattice of points, which encompasses the entire computational field. The particle distribution functions stand for the probability of particle presence with a specific velocity at each lattice or grid site. The macroscopic variables of the fluid, such as velocity or density, can be derived from these distribution functions.
The LBM consists of two fundamental steps, namely, the streaming-step and the collision-step. For every time step, distribution functions are distributed and streamed from each site to the neighboring sites in the streaming-step, and then the collision-step is conducted with the distribution functions relaxed towards a local equilibrium on the basis of the new macroscopic variables [32]. Here, the coupling flow of the methane gas and airflow belongs to the species transportation. The coupling model of the velocity and concentration based on 2D LBM was deduced, and the simulation procedure for velocity field and concentration field is presented. The Bossinesq approximation method is utilized to couple these two fields.

Lattice Boltzmann Method for Velocity Field
The LBM is commonly labelled as DdQq, where d stands for the space dimension and q is the number of microscopic lattice velocity directions. The possible nodes for 2D lattices are D2Q5, D2Q9. In this study, the D2Q9 model is employed to carry out the airflow velocity field, because it maintains good isotropy of the lattice. The airflow particles distribution on a lattice is shown in Figure 2a. including space, time, and velocities. The movement and distributions of the fluid are described as particle distribution functions residing at the sites of a regular lattice of points, which encompasses the entire computational field. The particle distribution functions stand for the probability of particle presence with a specific velocity at each lattice or grid site. The macroscopic variables of the fluid, such as velocity or density, can be derived from these distribution functions.
The LBM consists of two fundamental steps, namely, the streaming-step and the collision-step. For every time step, distribution functions are distributed and streamed from each site to the neighboring sites in the streaming-step, and then the collision-step is conducted with the distribution functions relaxed towards a local equilibrium on the basis of the new macroscopic variables [32]. Here, the coupling flow of the methane gas and airflow belongs to the species transportation. The coupling model of the velocity and concentration based on 2D LBM was deduced, and the simulation procedure for velocity field and concentration field is presented. The Bossinesq approximation method is utilized to couple these two fields.

Lattice Boltzmann Method for Velocity Field
The LBM is commonly labelled as DdQq, where d stands for the space dimension and q is the number of microscopic lattice velocity directions. The possible nodes for 2D lattices are D2Q5, D2Q9. In this study, the D2Q9 model is employed to carry out the airflow velocity field, because it maintains good isotropy of the lattice. The airflow particles distribution on a lattice is shown in Figure 2a. The fictitious fluid particles are subject to the lattice Boltzmann equation which describes the evolution of distribution function i f on Cartesian grids with a selective velocity i e , c is the movement velocity, ( ) i f Ω denotes the discrete collision operator, and t Δ is the time step. The complex collision operator is approximated by using the standard Bhatnagar-Gross-Krook (BGK) scheme [43], which states that the distribution functions is close to a local equilibrium distribution function eq i f and relaxes toward this equilibrium with a characteristic time τ. The evolution of the distribution functions using the BGK collision is described by the following equation: where τ is the relaxation factor (viscosity of fluid is related to the relaxation time), eq i f depends on the density ρ and velocity u of the gas, and is typically defined as The fictitious fluid particles are subject to the lattice Boltzmann equation which describes the evolution of distribution function f i on Cartesian grids with a selective velocity e i , c is the movement velocity, Ω i ( f ) denotes the discrete collision operator, and ∆t is the time step. The complex collision operator is approximated by using the standard Bhatnagar-Gross-Krook (BGK) scheme [43], which states that the distribution functions is close to a local equilibrium distribution function f eq i and relaxes toward this equilibrium with a characteristic time τ. The evolution of the distribution functions using the BGK collision is described by the following equation: where τ is the relaxation factor (viscosity of fluid is related to the relaxation time), f eq i depends on the density ρ and velocity u of the gas, and is typically defined as where ω i is the model-dependent weight coefficient c s = √ RT (R is the gas constant) is the lattice sound speed. For isothermal flows, c s is set to be c/ √ 3 with c = ∆x/∆t, where ∆x is the lattice spacing (c = 1 in this paper). Through the Chapman-Enskog expansion, the macroscopic fluid density ρ, and velocity u, can be derived as the zeroth and first order moments of f i respectively, The fluid pressure is defined directly as p = c 2 s ρ, the velocity set and the corresponding weight coefficients are defined as The above model can deduce the Navier-Stokes equation, which describes airflow regularity in the laneway.

Concentration Field through a Coupled Model
A coupled model is used to simulate the methane concentration field, where a D2Q9 lattice with a BGK collision operator is also used to solve the velocity and density, while the methane gas concentration is solved on a smaller and separate D2Q5 lattice. The distribution of concentration particles on a lattice is shown in Figure 2b.
The evolution equation for the methane concentration field is described by the evolution of distribution function C i on the Cartesian grids with selective velocity e i , where c denotes the velocity, ∆t is the time step, and Ψ i (C) is the discrete collision operator, which is given by, τ c = 2∆tD/∆x 2 + 0.5 (11) in which R is the source item, C eq i is the equilibrium distribution function of the methane concentration, and depends on the velocity u, temperature T of the gas, which is typically defined as τ c is the relaxation factor, D is the methane gas diffusion coefficient, c s and ∆x is the same as mentioned above, the direction set is defined as C(x, t) is the methane concentration in lattice x at time t, which is calculated by This model can also deduce the macroscopic equations for gas concentration diffusion.
The Boussinesq approximation method is employed to implement the coupling mode of the airflow velocity and the methane concentration field; the evolution equation after coupling is where g i , k i and C 0 are the gravitational acceleration, the volume expansion coefficient of concentration and reference concentration constant respectively.

The Turbulence Model
Airflow behavior in underground coal mine laneway is generally turbulent, hence our flow solver based on LBM should take the effect of turbulence into account. Therefore, a wide range of scales of fluid motion existed in the coupling flow of the air and methane concentration is necessary to resolve. The simulation in all scales in a turbulence flow requires a very fine lattice and long computation time. Instead, a Smagorinsky sub-lattice model is utilized to simulate the effects of the unresolved sub-lattice motion on the resolved motion, which is similar to the large eddy simulation (LES) in traditional CFD methods [44]. As for the LBM evolution equation, the effect of the sub-lattice can be addressed by local relaxation time [45][46][47]. Then the modified relaxation time is adopted in the relaxation process, and every node of lattice would relax at different rates.

Implementation of the LBM Integrated with LrGIS Platform
Numerical model and GIS integration can be achieved using one of the three potential coupling strategies outlined in Table 1 [28]. Here, a tight coupling approach is used to integrate LBM in LrGIS platform, which requires a straightforward connection from the geometric laneway data to LBM based numerical model. An applicable numerical coupling model was developed, which was specifically designed to adapt geometric laneway data to be integrated with LrGIS platform. The implementation procedure focuses on improving the geometry of the input features that configure the model geometry, the boundary conditions, and other aspects of spatiotemporal modeling, which avoids explicit mesh generation procedure that is necessary in traditional CFD method. It will help modelers in the development of future projects. The coupled platform merges different modules in a single powerful system, which avoids inconsistency and data loss originating from redundancy and heterogeneity of method base.
Data exchange is automatic between GIS platform and model.
Programming and data management is significantly complex, and the code base is not easy to be changed due to embedded large source code structure.
Steerable simulation according to the changes of parameters or processes.
Examples [48] [49] [50] Figure 3 shows the entire architecture of numerical model integrated with LrGIS platform. The platform consists of three layers: Data source, GIS kernel, and desktop application. Different kinds of data source can be quickly stored, retrieved, indexed, and searched by SDE interface layer. GeoLattice can be constructed on LrGIS platform to prepare for the numerical simulation, like other geometric models such as GeoPoint, GeoLine. In terms of the data management, LrGIS contains many entity models, such as GeoEntity, RasterEntity. Similarly, FluentEntity is specifically designed to manage the field data obtained by fluid simulation. The physical attributes consist of velocity, pressure, temperature, gas concentration, and display attributes comprise geometrical shape, color rendering, profile. The simulated result can be visualized on LrGIS GUI and execute the spatial query function. Specifically, two main parts have been solved: (i) GeoLattice and FluentEntity setup, which constructs the laneway model and configures the geometry and attributes of each lattice to make them valid for numerical computation; (ii) Database organization, which controls and traces the output data that are obtained from the numerical simulation. The data is saved as time series grid table files, which can be called by FluentEntity.  The lattice attributes and the corresponding boundary conditions are described in details: different kinds of lattice nodes on the geometry are associated with different material numbers in their own coordinates. Figure 4 presents a 2D top view of the laneway, with material number one, two, three, four and five corresponding to fluid region, inflow, outflow, methane mass flow, and The lattice attributes and the corresponding boundary conditions are described in details: different kinds of lattice nodes on the geometry are associated with different material numbers in their own coordinates. Figure 4 presents a 2D top view of the laneway, with material number one, two, three, four and five corresponding to fluid region, inflow, outflow, methane mass flow, and bounce back boundary, respectively. Then different material numbers are defined with the corresponding dynamic conditions, which are described, as follows: (i) fluid region: all species can move in this area, namely, the mesoscopic streaming and collision step of particles happened; (ii) inflow: the velocity inflow boundary condition is adopted, which means that the airflow is through this boundary with specified velocity; (iii) outflow: the outlet of the laneway is set to be under the constant pressure boundary condition; (iv) methane mass flow: the methane gas is released evenly from working face area with a total flow rate of 0.2 m 3 /s; (v) bounced back boundary: the no-slip boundary is prescribed at all laneway walls.  The detailed algorithm of LBM based numerical model integrated with LongRuanGIS platform is implemented as follows (Table 2), where the computational complexity of this numerical model is the product of lattice amount and number of iterations.

Algorithm. The LBM based turbulent velocity-concentration coupling model integrated with LrGIS platform
Step 1: Geometric model construction with specific geospatial coordinates.
Step 2: Lattice attributes configuration for two different lattice (D2Q9 and D2Q5) Step  The detailed algorithm of LBM based numerical model integrated with LongRuanGIS platform is implemented as follows (Table 2), where the computational complexity of this numerical model is the product of lattice amount and number of iterations. Table 2. Algorithm steps of the numerical model.

Algorithm. The LBM Based Turbulent Velocity-Concentration Coupling Model Integrated with LrGIS Platform
Step 1: Geometric model construction with specific geospatial coordinates.
Step 2: Lattice attributes configuration for two different lattice (D2Q9 and D2Q5) Step 3: Main loop starts Parameters initialization: f i , f eq i , C i , C eq i , ∆x, ∆t, maxIter For (it = 0; it < maxIter; ++it) If (it = 0) Set static boundary conditions (fluid region, wall). Else if (it != 0) Set dynamic boundary conditions(inflow, outflow and methane emission rate). Do collision step Step 4: Save the result data as time series grid files Step 5: FluentEntity call time series grid files and display field data on LrGIS GUI integrated with various coal mine map objects

Spatiotemporal Characteristics of Airflow and Methane Distribution
The established model was numerically simulated and all terms of the evolution equations, including streaming step and collision step were solved. The coupling mode is solved using Boussinesq approximation method. The detailed parameters are presented in Table 3. The simulation required around 1.6 h on workstations with four core processors and 8 GB RAM. The velocity field of airflow is one of the main factors which directly influence the methane distribution and mitigation. As presented in Figures 5 and 6, the predicted air velocity and methane concentration profiles were investigated at different time points with a fixed methane emission rate of 0.2 m 3 /s. Here, we analysed the overall spatiotemporal characteristics of velocity behaviour and methane distribution inside the laneway. The pressing air is flowing from the intake side of the laneway at high velocity, and there exists an obvious jet flow, which goes along the laneway and flows into the working face. When the airflow arrives at the working face, it is obstructed by the face and the airflow changes its flow direction and flows along the working face area, as can be seen in Figure 5a-c about the airflow distribution at different times. Correspondingly, the methane gas gushes from the coal seam and disperses into the working face with the impact of airflow after about 5 min. (Figure 6a), when the airflow is distributed fully in the working face area. Then the methane gas gradually diluted and decreased to a low concentration level with the effect of airflow in the wind outtake side of the laneway and approximately becomes stable at 10 min. (Figure 6b-d).

The Spatiotemporal Analysis of Methane Concentration Based on LrGIS
Based on the proposed integration method, the CFD simulation function is developed on LrGIS platform, which supports the numerical simulation and can display the simulated data on the ventilation system map on the basis of coordinate matching, and also facilitate the visualization of additional relevant data. The effective user access, the location-based information, and means of viewing and analyzing for the simulated data are also the main focuses. Thus, the numerical model is integrated with LrGIS platform, of which the query and analysis functions can also be operated

The Spatiotemporal Analysis of Methane Concentration Based on LrGIS
Based on the proposed integration method, the CFD simulation function is developed on LrGIS platform, which supports the numerical simulation and can display the simulated data on the ventilation system map on the basis of coordinate matching, and also facilitate the visualization of additional relevant data. The effective user access, the location-based information, and means of viewing and analyzing for the simulated data are also the main focuses. Thus, the numerical model is integrated with LrGIS platform, of which the query and analysis functions can also be operated

The Spatiotemporal Analysis of Methane Concentration Based on LrGIS
Based on the proposed integration method, the CFD simulation function is developed on LrGIS platform, which supports the numerical simulation and can display the simulated data on the ventilation system map on the basis of coordinate matching, and also facilitate the visualization of additional relevant data. The effective user access, the location-based information, and means of viewing and analyzing for the simulated data are also the main focuses. Thus, the numerical model is integrated with LrGIS platform, of which the query and analysis functions can also be operated based on the numerical simulation results. Figure 7 shows the map of ventilation system of 5# Sijiazhuang coal mine, the study area of 15116 working face in this case is located in the bottom right corner on the map, with the latitude and longitude ranging from 105,800 to 106,020 m and from 67,125 to 67,725 m. The roam and move functions allow users to observe more detailed features on the coal mine map.  To demonstrate the applicability of the CFD simulation function in LrGIS platform, the spatiotemporal analysis for our case study is performed. First, based on the specified numerical parameters, the geometry is prepared and lattice is generated, which are then used to perform the numerical simulating procedure. It is noted that the properties of lattices include the lattice ID, material numbers, and the coordinates. The output data of the simulation is saved as a series of time ordered text files, which contains the X, Y coordinates, the airflow velocity, and methane concentration of each lattice. The tool of methane emission analysis in the CFD simulation model allows for us to demonstrate the spatiotemporal distribution of the airflow behavior and methane concentration. By specifying the time point, the corresponding spatiotemporal distribution of the methane concentration can be visualized on the ventilation system map, together with various spatial coal mine layers, including geographical locations of fans, methane sensors and first aid stations, contour line, measured faultage, river and railway, and also the fictitious graticule, as presented in Figure 8. These information provide managers and miners the access to the large quantities of geographic field data, as well as its query and analysis for the location-based information integrated with simulated data.
Basically, the managers or decision makers can use the program to perform several queries according to the specifications, such as, (i) basic mine mapping operations integrated with simulated field data, (ii) finding operation, and (iii) location-based information operation. For example, the maximum methane concentration with its coordinates in an area of interest will be demonstrated immediately by the function of polygonal query, which also shows whether the predicted concentration exceeds the threshold value, as can be seen in Figure 9. To demonstrate the applicability of the CFD simulation function in LrGIS platform, the spatiotemporal analysis for our case study is performed. First, based on the specified numerical parameters, the geometry is prepared and lattice is generated, which are then used to perform the numerical simulating procedure. It is noted that the properties of lattices include the lattice ID, material numbers, and the coordinates. The output data of the simulation is saved as a series of time ordered text files, which contains the X, Y coordinates, the airflow velocity, and methane concentration of each lattice. The tool of methane emission analysis in the CFD simulation model allows for us to demonstrate the spatiotemporal distribution of the airflow behavior and methane concentration. By specifying the time point, the corresponding spatiotemporal distribution of the methane concentration can be visualized on the ventilation system map, together with various spatial coal mine layers, including geographical locations of fans, methane sensors and first aid stations, contour line, measured faultage, river and railway, and also the fictitious graticule, as presented in Figure 8. These information provide managers and miners the access to the large quantities of geographic field data, as well as its query and analysis for the location-based information integrated with simulated data.
Basically, the managers or decision makers can use the program to perform several queries according to the specifications, such as, (i) basic mine mapping operations integrated with simulated field data, (ii) finding operation, and (iii) location-based information operation. For example, the maximum methane concentration with its coordinates in an area of interest will be demonstrated immediately by the function of polygonal query, which also shows whether the predicted concentration exceeds the threshold value, as can be seen in Figure 9.  Besides, if a specific point is of interest in an actual field, for example, the position is where miners always work, then the point can be clicked or the coordinate can be input on the interface to get the methane concentration curve changed over time, which will be displayed on the interface and provides intuitive information for the dynamic regularity of methane concentration. The time series data can also be validated with monitored data based on the database management system. Here, one of the methane sensors located 10 m away from the working face in the outtake side of the  Besides, if a specific point is of interest in an actual field, for example, the position is where miners always work, then the point can be clicked or the coordinate can be input on the interface to get the methane concentration curve changed over time, which will be displayed on the interface and provides intuitive information for the dynamic regularity of methane concentration. The time series data can also be validated with monitored data based on the database management system. Here, one of the methane sensors located 10 m away from the working face in the outtake side of the laneway is selected to perform the data validation with the monitored data, as presented in Figure   Figure 9. Methane concentration visualization and polygonal query analysis.
Besides, if a specific point is of interest in an actual field, for example, the position is where miners always work, then the point can be clicked or the coordinate can be input on the interface to get the methane concentration curve changed over time, which will be displayed on the interface and provides intuitive information for the dynamic regularity of methane concentration. The time series data can also be validated with monitored data based on the database management system. Here, one of the methane sensors located 10 m away from the working face in the outtake side of the laneway is selected to perform the data validation with the monitored data, as presented in Figure 10. The numerical simulation shows that, after the algorithm converges in 10 min, and the airflow distributes fully and stably within the laneway. The simulated methane concentration maintains at a stable level and agrees well with the monitored data after 10 min. Thus, some detailed information can be provided for decision makers, such as how long the methane distribution will be stable with the effect of the ventilation system, the peak value of methane concentration, and its specific position, which can guide miners about when to avoid the methane accident if the value is beyond the maximum allowable value. Besides, the methane distribution can also guide the production intensity, and reasonable arrangement of underground coal mine equipment. distributes fully and stably within the laneway. The simulated methane concentration maintains at a stable level and agrees well with the monitored data after 10 min. Thus, some detailed information can be provided for decision makers, such as how long the methane distribution will be stable with the effect of the ventilation system, the peak value of methane concentration, and its specific position, which can guide miners about when to avoid the methane accident if the value is beyond the maximum allowable value. Besides, the methane distribution can also guide the production intensity, and reasonable arrangement of underground coal mine equipment.

The Comprehensive Comparison of Traditional CFD Method and LBM
According to the kinetic theory [40], both traditional CFD method and LBM can satisfy the macroscopic Navier-Stokes equation in their continuous forms. However, Navier-Stokes equations are the target for LBM while the starting point for the traditional CFD methods, and these two methods also differ, not only in the spatial discretization forms, but also in the temporal discretization forms. Although a number of studies have been conducted to comprehensively compare their difference in essence, there is still a lack of comparison based on the actual simulation example. In this case study, comprehensive comparison and analysis between the proposed turbulent LBM based velocity-concentration coupling model and traditional CFD method are conducted from the following three aspects: accuracy, simulation efficiency, simplicity, and advantages of the integration with LrGIS platform.
Here, the traditional CFD method was implemented based on Computational Fluid Dynamics platform (Fluent) by employing the Euler-Lagrange model, all terms of the governing equations and constitutive relations were solved with the semi-implicit pressure-linked equation (SIMPLE) algorithm and second-order upwind scheme. The turbulence model and boundary conditions were solved using finite volume solver. The dynamic boundary parameters are the same with the LBM.

The Comprehensive Comparison of Traditional CFD Method and LBM
According to the kinetic theory [40], both traditional CFD method and LBM can satisfy the macroscopic Navier-Stokes equation in their continuous forms. However, Navier-Stokes equations are the target for LBM while the starting point for the traditional CFD methods, and these two methods also differ, not only in the spatial discretization forms, but also in the temporal discretization forms. Although a number of studies have been conducted to comprehensively compare their difference in essence, there is still a lack of comparison based on the actual simulation example. In this case study, comprehensive comparison and analysis between the proposed turbulent LBM based velocity-concentration coupling model and traditional CFD method are conducted from the following three aspects: accuracy, simulation efficiency, simplicity, and advantages of the integration with LrGIS platform.
Here, the traditional CFD method was implemented based on Computational Fluid Dynamics platform (Fluent) by employing the Euler-Lagrange model, all terms of the governing equations and constitutive relations were solved with the semi-implicit pressure-linked equation (SIMPLE) algorithm and second-order upwind scheme. The turbulence model and boundary conditions were solved using finite volume solver. The dynamic boundary parameters are the same with the LBM. Here, a line at 3 m high in the middle of the outtake side of laneway was chosen to investigate the airflow velocity and methane concentration at two different time points. As can be seen in Figures 11 and 12, it shows that both methods agree well with each other in the overall tendency of airflow velocity and methane concentration.   The traditional CFD simulation requires about 20 min on computers with four core processors and 8 GB RAM, while the model adopted in this paper uses around 15 min. It is noted that Fluent belongs to the third-party professional fluid dynamics simulation software, its independent geometric modeling and numerical simulation process makes it significantly difficult to directly guide the on-site problems, and the required data have to be preset and saved manually in the early stage of the simulation, which shows the poor practicability when applied in the on-site guidance. However, in this paper, the property of the unnecessary explicit mesh generation and the inherent simplicity of the algorithm make LBM more efficient and easy to integrate with the LrGIS platform. Specifically, the integrated system essentially implemented four objectives: (i) enabling managers to simulate case study directly on the GIS system and visualize the simulated results intuitively; (ii) storing the simulated spatiotemporal field data in the database system, which can be used for post analysis; (iii) demonstrating the simulated field data with various coal mine map objects, and providing miners the spatiotemporal analyze and query operation tools to obtain the useful locationbased information; (iv) validating the results based on the coal mine monitoring system of Sijiazhuang coal mine. Therefore, by utilizing the synergies between GISs and spatiotemporal modeling, the integrated system is effective in providing decision supports for safety production, which is expected to improve the health standards in the underground coal mines.
The comparison case study demonstrates that the model and integration mode carried out in   The traditional CFD simulation requires about 20 min on computers with four core processors and 8 GB RAM, while the model adopted in this paper uses around 15 min. It is noted that Fluent belongs to the third-party professional fluid dynamics simulation software, its independent geometric modeling and numerical simulation process makes it significantly difficult to directly guide the on-site problems, and the required data have to be preset and saved manually in the early stage of the simulation, which shows the poor practicability when applied in the on-site guidance. However, in this paper, the property of the unnecessary explicit mesh generation and the inherent simplicity of the algorithm make LBM more efficient and easy to integrate with the LrGIS platform. Specifically, the integrated system essentially implemented four objectives: (i) enabling managers to simulate case study directly on the GIS system and visualize the simulated results intuitively; (ii) storing the simulated spatiotemporal field data in the database system, which can be used for post analysis; (iii) demonstrating the simulated field data with various coal mine map objects, and providing miners the spatiotemporal analyze and query operation tools to obtain the useful locationbased information; (iv) validating the results based on the coal mine monitoring system of Sijiazhuang coal mine. Therefore, by utilizing the synergies between GISs and spatiotemporal modeling, the integrated system is effective in providing decision supports for safety production, which is expected to improve the health standards in the underground coal mines.
The comparison case study demonstrates that the model and integration mode carried out in this study is more efficient in improving the computational efficiency and simplifying the integration mode of numerical simulation model with the GIS platform. It can also be concluded that the The traditional CFD simulation requires about 20 min on computers with four core processors and 8 GB RAM, while the model adopted in this paper uses around 15 min. It is noted that Fluent belongs to the third-party professional fluid dynamics simulation software, its independent geometric modeling and numerical simulation process makes it significantly difficult to directly guide the on-site problems, and the required data have to be preset and saved manually in the early stage of the simulation, which shows the poor practicability when applied in the on-site guidance. However, in this paper, the property of the unnecessary explicit mesh generation and the inherent simplicity of the algorithm make LBM more efficient and easy to integrate with the LrGIS platform. Specifically, the integrated system essentially implemented four objectives: (i) enabling managers to simulate case study directly on the GIS system and visualize the simulated results intuitively; (ii) storing the simulated spatiotemporal field data in the database system, which can be used for post analysis; (iii) demonstrating the simulated field data with various coal mine map objects, and providing miners the spatiotemporal analyze and query operation tools to obtain the useful location-based information; (iv) validating the results based on the coal mine monitoring system of Sijiazhuang coal mine. Therefore, by utilizing the synergies between GISs and spatiotemporal modeling, the integrated system is effective in providing decision supports for safety production, which is expected to improve the health standards in the underground coal mines.
The comparison case study demonstrates that the model and integration mode carried out in this study is more efficient in improving the computational efficiency and simplifying the integration mode of numerical simulation model with the GIS platform. It can also be concluded that the integrated system is more practical in on-site application compared with numerical simulation software.

Conclusions
Due to the lack of an effective platform for the practical application of numerical simulation methods in the underground coal mine systems, the concept of tight coupling of GIS and spatiotemporal modeling is introduced in this area for the first time.
Specifically, a turbulent LBM based velocity-concentration coupling model is proposed and tightly integrated with LrGIS platform by developing the FluentEntity model, which allows mine managers to simulate ventilation-related case study directly on the LrGIS system, visualize the simulated field data on the ventilation system map and other coal mine layers intuitively, and operate the spatiotemporal analyze and query to get location-based information, and therefore to guide the practical operation. Besides, the simulated result is validated with monitored data of Sijiazhuang underground coal mine, which shows the accuracy and reliability of the proposed model. In addition, by comparing with the traditional CFD method, the developed coupling scheme is demonstrated to be more efficient and easy to be integrated with GIS system.
In conclusion, this study provides a possibility to integrate the numerical simulation model with the GIS platform effectively, and shows that GIS could not only serve to carry out general constructive safety, but also generate case scenarios, quantify spatiotemporal processes, and therefore improve the safety standards for coal mines by GIS-coupled spatiotemporal modeling, analysis, and query tools. Furthermore, the system is more flexible and potentially applicable to other case study, like ventilation network solution, including both local network and the whole mine ventilation system.