According to the analysis of Jipin Sun , from China Mine University, in terms of mine accidentcategory, the death rate and national mine accident rate: roof 52.7%, and 36.8%, gas 11.3% and 29.7%, transportation: 16.9% and 11.3%, water disaster 3.1% and 8.1%, electrical-mechanic 4.2% and 2.5%, machine boom 2.7% and 1.9%, fire disaster 0.4% and 1.9% from 2004 to 2013. The water and fire accidents caused high death rate, others reversal. We think the high proportion is due to the unpredictable terrain. Compared to gas, water, fire accidents, transportation, machine boom accidents are easily avoidable if we operate properly. In theory, the distribution of gas and water is fractal. In other words, no matter how short of an interval of drilling the gas releasing hole, there is a possibility that we can miss an origin of gas fire. In order to avoid unnecessary loss, it is vital to predict gas peaking time and the location of the gas leakage.
2. Design Theory and Process
2.1. Design Theory
1) Acquisition parameters:
u Positive pressure: converted from negative pressure of drainage
u Temperature: air temperature of drainage pipeline
u Gas level: instantaneous concentration of gas in drainage pipeline
u Working flow: instantaneous flow rate* flow coefficient
2) The explanation of the formulae:
a) Pearson Product-moment Correlation Coefficient: PPCC, describes the linear correlation of two random variables. Karl Pearson put this forward in the 1880s . It is widely accepted in many aspects of science.
The definition of the coefficient:
where the subscribe i runs from 1 to n, n is the number of terms.
b) Hypo-(or Huang) Product-moment Correlation Coefficient: HPCC, measures the nonlinear correlation of two variables. It was put forward by Prof. Steed Huang in 2015 .
The definition of Correlation Coefficient
Huang Product-moment Correlation Coefficient is used to measure the nonlinear correlation of three variables. It can be generalized from the form of two variables to three and can be used in many aspects of industry.
c) Three Factor  Correlation Coefficient (TFCC) definition:
d) Multiple Factor  Correlation Coefficient: If we want to analyze the correlation of four variables, density, temperature, pressure, flow velocity, of gas, then we should define the Multiple Factors Correlation Coefficient (MFCC):
where the subscribe i runs from 1 to n, n is the number of terms, the subscribe j runs from 2 to N, N is the number of variables. It is similar to the Fractional Variance  defined by the Gamma function.
3) Formula comparison
a) The difference and advantage of Huang Product?moment Correlation Coefficient compared with Pearson Product-moment Correlation Coefficient: by using the 1.5 times instead of square, the new definition can reflect the fluctuating asymmetry, fitting the nature of things more accurately and can be used for multiple factors in more aspects of the fields.
b) The difference, disadvantage and advantage by using field between Hypo Product-moment Correlation Coefficient: The new definition, not square, but with 1.5 parties, can reflect the asymmetry of fluctuation and can be extended to the occasion of multiple targets; the disadvantage is that the calculation is relatively complex. If the variable is strictly consistent with the Gauss distribution and the fluctuation is statistically symmetric or only involves two target variables, it is not necessary to use the new definition.
4) The Brain Storm Optimization algorithm
Brain Storm Optimization (BSO) algorithm  and then correlation criterion of gas outburst prediction orientation curve, its principle diagram is as follows in Figure 1.
2.2. Maintaining the Integrity of the Specifications
1) Calculation methods and steps
According to the gas concentration, we treat each pumping station as a set of data units, the ground station data as a high concentration data group, and the
Figure 1. Example of brain storming optimization algorithm.
mobile station data as a low concentration of data.
Then take the entire mining area of low or high concentrations average as a benchmark of respected concentration average for optimization, searching for the difference between the low and high concentration with optimized correlation coefficient, thus inversion derivation and deduction of geological characteristics of coal and gas outburst extrusioncurve.
Assume that the mine space is closed, or said each ventilation are all under observation, it has the integrity of the data, or said outburst only through the monitoring of air discharge, no other unseen conduits, for us to meet the law of conservation for the temperature, pressure, velocity and concentration.
2) Pattern analysis
We found that the gas outburst is cyclical, at about every two to three climaxes a day, estimation that is related to the earth core has the eccentric phase rotary extrusion towards the crust of our planet.
The next task is to establish the monthly benchmark data. According to the data on that day and of input variables: gas concentration. We use brainstorm to optimize and look for the cycle time of today to predict the trends for tomorrow.
The topology of the ground station and the mobile station is needed to calculate the correlation coefficient of the location, and using the position, to predict the periodicity of the space pattern and predict the exact danger area.
3. Design Procedures and Results
3.1. Program Instructions
a(day).xls is the original data contains day 15 to day 31.
a15-a31 compressed package is the original data (standby)
BsoMain.m is the main program and brain storm subroutine
CfdMg.m is a complex number of variance criterion subroutine
chuli.m is the three correlation coefficient calculation subroutine
・ function z = CfdMg(x)
・ n = length(x)-1;
・ z = 0;
・ global x24
・ x24 = round(x*5);
・ gasHig= xlsread('date.xls','a1:x1');
・ gasLow= xlsread('date.xls','a2:x2);
・ THig= xlsread('date.xls','a3:x3')';
・ TLow= xlsread('date.xls','a4:x4')';
・ PHig= xlsread('date.xls','a5:x5')';
・ PLow= xlsread('date.xls','a5:x5')';
・ gasHig= xlsread('date.xls','a6:x6')';
・ THigH=THig *1.8+32*ones(24,1);
・ TLowH=TLow *1.8+32*ones(24,1);
More of the related code is available for download at the MATLAB server site.
3.2. Calculation Results
We found that the gas concentration changes periodically. There are about three peaks every day due to the estimation of the earth’s core and the eccentric rotary phase extrusion relative to the earth crust. The first picture of Figure 2 shows the calculation results of the distance between each round of the brainstorming algorithm and the final optimization goal. The blue (high) curve in the second picture of Figure 3 is our predictive value for the data. The red (low) curve is the
Figure 2. Example of brain storming optimization fitness.
Figure 3. Example of predicted gas emission cycle.
trend curve after the average value is taken every four hours.
This experiment was achieved by performing simulation though the Matlab software . It is more convenient to use the MATLAB software to calculate the fractional variance but in principle, any other scientific computing software can achieve a similar function. The purpose of this study is through the prediction of mine gas using the MATLAB software, we can collect the temperature, pressure, ground station and mobile station of gas flow in the mines for the analysis of correlation, correlation coefficient, and using the brains-storm algorithm to optimize the prediction of the gas emission peaking and the future trend.
Some difficulties are: the algorithm chuli.m changes for ternary correlation coefficient algorithm from binary case, and in CfdMg.m function when modify initial value, call the function parameters, using a formula to turn in between Fahrenheit and Celsius, or the negative pressure is getting to positive pressure. We finally observed in graphics, estimated daily high gas concentration duration of peak hours. We recommend 8 hours a shift so that miners will be safe and can rest during the high gas emission period and return to work on a low gas time period so that the explosion risk is reduced to the lowest level.
Thanks to Professor Yuhui Shi from SUST for providing the brain storming algorithm plus subroutines, thanks to Mr. Zhao and Mr. Ding from China Coal provide the coal mine data sets, and thanks to Ms. Sun Chenxu and 98 Suqian students to assist us in completing the original data calculation and verification for the 25 working days records, that gives us the confidence for seeing the same pattern, finally, thanks to Kenneth Huang from UCSD for editing the English writing.
 Zou, Q., Hu, Y.F. and Huang, J.S. (2015) Definition of Complex Hurst and Fractional Analysis for Stock Market Fluctuation , 2nd International Conference on Industrial Engineering, Management Science and Applications.
 Braun, T., Rich, M. and Kramer, M.F. (2013) Correlation of Three Variables Describing Nasal Patency (HD, MCA, NOSE Score) in Healthy Subjects. Brazilian Journal of Otorhinolaryngology, 79, 354-358. https://doi.org/10.5935/1808-8694.20130062