CN105804739A - Prediction method and prediction system for gas emission quantity - Google Patents

Abstract

The invention discloses a prediction method for predicting gas emission quantity of a horizontal layered mining working face of an inclined coal seam, which comprises the following steps: s001: investigating and calculating working face gas prediction basic data; s002: calculating the gas emission quantity q of the mining layer₁(ii) a S003: calculating the gas emission quantity q of the adjacent coal bed adjacent to the mined coal bed in the stope face₂(ii) a S004: calculating the stoperPressure relief gas emission q at the lower part of the working face₃(ii) a S005: calculating the gas emission quantity q of the mining working face_MiningWherein q is_Mining＝q₁+q₂+q₃. The invention also discloses a prediction system for predicting the gas emission quantity of the horizontal layered mining working face of the inclined coal seam. According to the prediction method and the prediction system for the gas emission quantity, provided by the invention, a prediction model for the gas emission quantity of the horizontal layered mining working face of the steep ultra-thick coal seam is established, and the defects of the prior art are overcome.

Description

The Forecasting Methodology of gas emission and prognoses system

Technical field

The present invention relates to Predicting Gas technical field in coal mining, particularly relate to Forecasting Methodology and the prognoses system of a kind of gas emission.

Background technology

Forecast of Gas Emission is to carry out the key factor of mine, exploiting field, horizontal ventilation design, gas administration and gas pumping design.Forecast of Gas Emission is as a kind of important technological means, and it predicts whether accurately and order of accuarcy all drastically influence safety in production and the economic benefit of mine, thus Forecast of Gas Emission has very important meaning for mine safety production.

Current face gas prediction of gas emission method is a lot, mainly has following three kinds according to forecast model and target prediction method: the first is the mine statistic law being primarily adapted for use in prediction geological conditions simply, based on mathematical statistics；The gas bearing capacity method that the second is is fundamental forecasting parameter with coal seam gas-bearing capacity；The third is widely used in point source predicted method of New Built Mine.

Wherein, mine statistic law needs sufficiently large sample space, i.e. the abundant actual measurement gas data of this mine or contiguous mine accumulation, and the most important thing is to produce the similarity degree of district and Target area, and its accuracy and range of application are limited.

When adopting gas bearing capacity meter algorithm predicts gas emission, it is possible to illustrate the separate sources of gas, but the distribution wanting to make coal seam gas-bearing capacity measured value, content point and density reach certain accuracy and have certain difficulty.In implementing at the scene, Application comparison difficulty, error is big, is difficult to ensure its precision of prediction of card.

Point source predicted method Emission Law according to gas sources different in Mine Production, the size of each source gas emission of gushing out and each gas source, is predicted each face gas outburst amount of mine.

For severe inclined thick coal seam horizontally grouped top-coal drawing mining method work surface, current industry still lacks the standard of the Forecast of Gas Emission method under this condition, if utilizing existing standard to be predicted there is certain defect.

Summary of the invention

It is an object of the invention to overcome defect of the prior art, it is provided that a kind of Forecasting Methodology and prognoses system carrying out Forecast of Gas Emission for inclined seam horizontally grouped top-coal drawing mining method work surface.

Technical solution of the present invention provides a kind of Forecasting Methodology carrying out Forecast of Gas Emission for inclined seam horizontally grouped top-coal drawing mining method work surface, comprises the steps:

S001: investigate evaluation work face Predicting Gas master data, including:

Country rock Gas COEFFICIENT K₁, work surface lose coal Gas COEFFICIENT K₂, primitive coalbed gas bearing capacity X₀, remaining coal seam gas-bearing capacity X_c, working seam the thickness m of i-th adjacent coal seam adjacent with described working seam for thickness m_i, described working seam working thickness m₀, adjacent with described working seam i-th adjacent coal seam gas bearing capacity X_i, adjacent with described working seam i-th adjacent coal seam residual gas content X_ic, adjacent with described working seam i-th adjacent coal seam is by the drawing-off gas rate ζ of mining influence_i, exploitation layering height M, described working seam tilt angle alpha, exploitation layering bottom coal body mining influence degree of depth h, working seam gas bearing capacity gradient X_t, working seam drawing-off gas rate gradient λ_t, exploitation layering the bottom maximum drawing-off gas rate λ of coal body_max；

S002: calculate the gas emission q of described exploitation layering₁；

S003: calculate the gas emission q of adjacent coal seam adjacent with described working seam in stope₂；

S004: calculate the bottom gas pressure relief outburst amount q of described stope₃；

S005: calculate the gas emission q of exploitation work surface_Adopt, wherein, q_Adopt=q₁+q₂+q₃。

Further, the gas emission q of described exploitation layering₁Calculated by equation below:

$q_1=K_1\×K_2\×\frac{m}{m_0}\×(X_0-X_C).$

Further, the gas emission q of described adjacent coal seam₂Calculated by equation below:

$q_2=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\frac{m_i}{m_0}\×{\mathrm{\ζ}}_i(X_i-X_{ic}).$

Further, bottom gas pressure relief outburst amount q₃Calculated by equation below:

$q_3=\frac{1}{M}{\cos }^{-1}\mathrm{\α}\[(X_0-X_c){\mathrm{\λ}}_{max}h+\frac{X_t{\mathrm{\λ}}_{max}+{\mathrm{\λ}}_t(X_0-X_C)}{2}{h}^2+\frac{X_t{\mathrm{\λ}}_t}{3}{h}^3\].$

Technical solution of the present invention also provides for a kind of prognoses system carrying out Forecast of Gas Emission for inclined seam horizontally grouped top-coal drawing mining method work surface, including:

Data collection computing module, for collecting and evaluation work face Predicting Gas master data, described data include:

Country rock Gas COEFFICIENT K₁, work surface lose coal Gas COEFFICIENT K₂, primitive coalbed gas bearing capacity X₀, remaining coal seam gas-bearing capacity X_c, working seam the thickness m of i-th adjacent coal seam adjacent with described working seam for thickness m_i, described working seam working thickness m₀, adjacent with described working seam i-th adjacent coal seam gas bearing capacity X_i, adjacent with described working seam i-th adjacent coal seam residual gas content X_ic, adjacent with described working seam i-th adjacent coal seam is by the drawing-off gas rate ζ of mining influence_i, exploitation layering height M, described working seam tilt angle alpha, exploitation layering bottom coal body mining influence degree of depth h, working seam gas bearing capacity gradient X_t, working seam drawing-off gas rate gradient λ_t, exploitation layering the bottom maximum drawing-off gas rate λ of coal body_max；

First computing module, for the data obtained according to described data collection computing module, calculates the gas emission q of exploitation layering₁；

Second computing module, for the data obtained according to described data collection computing module, calculates the gas emission q of adjacent coal seam adjacent with described working seam in stope₂；

3rd computing module, for the data obtained according to described data collection computing module, calculates the bottom gas pressure relief outburst amount q of described stope₃；

4th computing module, for the result of calculation according to described first computing module, described second computing module and described 3rd computing module, calculates the gas emission q of exploitation work surface_Adopt, wherein, q_Adopt=q₁+q₂+q₃。

Further, described first computing module adopts equation below to calculate q₁:

$q_1=K_1\×K_2\×\frac{m}{m_0}\×(X_0-X_C).$

Further, described second computing module adopts equation below to calculate q₂:

$q_2=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\frac{m_i}{m_0}\×{\mathrm{\ζ}}_i(X_i-X_{ic}).$

Further, described 3rd computing module adopts equation below to calculate q₃:

$q_3=\frac{1}{M}{\cos }^{-1}\mathrm{\α}\[(X_0-X_c){\mathrm{\λ}}_{max}h+\frac{X_t{\mathrm{\λ}}_{max}+{\mathrm{\λ}}_t(X_0-X_C)}{2}{h}^2+\frac{X_t{\mathrm{\λ}}_t}{3}{h}^3\].$

Adopt technique scheme, have the advantages that

The Forecasting Methodology of gas emission provided by the invention and prognoses system, it is applied in stope working surface of coal mines Forecast of Gas Emission, it is particularly applicable in severe inclined thick coal seam horizontal slice coal mining stope Forecast of Gas Emission, it establishes applicable severe inclined thick coal seam horizontally grouped top-coal drawing mining method face gas prediction of gas emission model, compensate for the deficiency that existing Emit Quantity Prediction Methods In Coal Mines is applied in severe inclined thick coal seam horizontally grouped top-coal drawing mining method face gas prediction of gas emission.

Accompanying drawing explanation

The step schematic diagram of the Forecasting Methodology of the gas emission that Fig. 1 provides for one embodiment of the invention；

The schematic diagram of the prognoses system of the gas emission that Fig. 2 provides for one embodiment of the invention.

Detailed description of the invention

The specific embodiment of the present invention is further illustrated below in conjunction with accompanying drawing.

Forecasting Methodology provided by the invention and prognoses system, carry out Predicting Gas mainly for the anxious super high seam that inclines.Suddenly incline the thickness of super high seam at more than 8m, and inclination alpha is between 45 °-90 °.

As it is shown in figure 1, a kind of Forecasting Methodology carrying out Forecast of Gas Emission for inclined seam horizontally grouped top-coal drawing mining method work surface that one embodiment of the invention provides, comprise the steps:

S001: investigate evaluation work face Predicting Gas master data, including:

Country rock Gas COEFFICIENT K₁, work surface lose coal Gas COEFFICIENT K₂, primitive coalbed gas bearing capacity X₀, remaining coal seam gas-bearing capacity X_c, working seam the thickness m of i-th adjacent coal seam adjacent with working seam for thickness m_i, working seam working thickness m₀, adjacent with working seam i-th adjacent coal seam gas bearing capacity X_i, adjacent with working seam i-th adjacent coal seam residual gas content X_ic, adjacent with working seam i-th adjacent coal seam is by the drawing-off gas rate ζ of mining influence_i, exploitation layering height M, working seam tilt angle alpha, exploitation layering bottom coal body mining influence degree of depth h, exploitation coal body gas bearing capacity gradient X_t, exploitation coal body drawing-off gas rate gradient λ_t, exploitation layering the bottom maximum drawing-off gas rate λ of coal body_max；

S002: calculate the gas emission q of exploitation layering₁；

S003: calculate the gas emission q of adjacent coal seam adjacent with working seam in stope₂；

S004: calculate the bottom gas pressure relief outburst amount q of stope₃；

S005: calculate the gas emission q of exploitation work surface_Adopt, wherein, q_Adopt=q₁+q₂+q₃。

Wherein, country rock Gas COEFFICIENT K₁Refer to, exploitation work surface roof and floor institute rich in the gas such as gas, Gas after being destroyed by seam mining, this part gush out gas with can desorbing gas proportional, this proportionality coefficient employing country rock Gas coefficient represents.

Work surface loses coal Gas COEFFICIENT K₂Referring to, work surface is at recovery process, owing to can not reclaim coal body completely, loses in coal and can produce a part of Gas, this part is gushed out gas with can desorbing gas ratio as losing coal Gas coefficient, generally adopt rate of extraction inverse to represent.

Primitive coalbed gas bearing capacity X₀Refer to, the gas bearing capacity that the coal seam of mining level is original.

Remaining coal seam gas-bearing capacity X_cRefer to, mine and transported the gas bearing capacity of coal seam remaining after mine.

Drawing-off gas rate ζ_iRefer to, owing to adjacent coal seam distance working seam has different distances, therefore it is different that different adjacent coal seams pours into the gas emission of exploitation work surface, using different adjacent coal seams to the gas emission in exploitation space and adjacent layer self can coal seams gas discharging quantity ratio as the next layer gas emission index.

Exploitation layering bottom coal body mining influence degree of depth h refers to, after work surface seam mining, and the depth capacity that lower coal body stress changes relative to the stress of primary rock.

Working seam gas bearing capacity gradient X_tReferring to, due to the increase of coal seam buried depth, gas bearing capacity is gradually increased, and gas bearing capacity gradient refers to the gas bearing capacity rate of change with buried depth direction, it is common that gas bearing capacity and pipe laying are carried out linear fit, and its fitting coefficient is gas bearing capacity gradient.

Working seam drawing-off gas rate gradient λ_tRefer to, due to severe inclined thick coal seam horizontally grouped top-coal drawing mining method, bottom to be exploited coal seam layering with exploit exploitation layering be connected, therefore the drawing-off gas rate of underlying seam layering has certain change along with the underlying seam degree of depth, and with the rate of change of distance mined bed distance, drawing-off gas rate is defined as drawing-off gas rate gradient.

The exploitation layering bottom maximum drawing-off gas rate λ of coal body_maxRefer to, owing to exploitation layering bottom coal body is to exploiting space Gas, coal body gas remainder quantity change over time in bottom is caused to be continually changing, maximum drawing-off gas rate be bottom coal body gas emission maximum with can desorbing mash gas content ratio (have under the conservation of mass: raw coal gas bearing capacity=remnants can desorbing mash gas content+gas emission+residual gas amount, usual gas emission maximum in theory can be equal to raw coal gas bearing capacity deduct residual gas amount).

The thickness m of working seam refers to, the coal seam integral thickness exploited.

The height M of exploitation layering, the height being namely layered, it is possible to be called height of lift.Wherein, working seam carries out multiple-zone production, and every layer of exploitation or layering coal seam to be exploited are referred to as exploitation layering.

The gas emission q of exploitation layering₁Refer to, the gas emission in exploitation layering coal seam, it is possible to be referred to as the gas emission in layering or segmentation coal seam.

The gas emission q of adjacent coal seam₂Refer to, the gas emission of the adjacent coal seam adjoining with working seam.

Bottom gas pressure relief outburst amount q₃Referring to, severe inclined thick coal seam is exploited, and work surface underlying seam is subject to mining influence, can produce release, cause that gas can pour into digging space, and to work surface this part gas amount of gushing out, work surface underlying seam is defined as bottom gas pressure relief outburst amount.

Thus, the Forecasting Methodology of gas emission provided by the invention, it is applied in stope working surface of coal mines Forecast of Gas Emission, it is particularly applicable in severe inclined thick coal seam horizontal slice coal mining stope Forecast of Gas Emission, it establishes applicable severe inclined thick coal seam horizontally grouped top-coal drawing mining method face gas prediction of gas emission model, compensate for the deficiency that existing Emit Quantity Prediction Methods In Coal Mines is applied in severe inclined thick coal seam horizontally grouped top-coal drawing mining method face gas prediction of gas emission.

It is preferred that the gas emission q of exploitation layering₁Calculated by equation below:

$q_1=K_1\×K_2\×\frac{m}{m_0}\×(X_0-X_C).$

It is preferred that the gas emission q of adjacent coal seam₂Calculated by equation below:

Wherein n refers to the number of plies of adjacent coal seam, and it is the natural number be more than or equal to 1, and 1≤i≤n.

It is preferred that bottom gas pressure relief outburst amount q₃Calculated by equation below:

$q_3=\frac{1}{M}{\cos }^{-1}\mathrm{\α}\[(X_0-X_c){\mathrm{\λ}}_{max}h+\frac{X_t{\mathrm{\λ}}_{max}+{\mathrm{\λ}}_t(X_0-X_C)}{2}{h}^2+\frac{X_t{\mathrm{\λ}}_t}{3}{h}^3\].$

As in figure 2 it is shown, a kind of prognoses system carrying out Forecast of Gas Emission for inclined seam horizontally grouped top-coal drawing mining method work surface that one embodiment of the invention provides, including:

Data collection computing module 1, for collecting and evaluation work face Predicting Gas master data, data include:

Country rock Gas COEFFICIENT K₁, work surface lose coal Gas COEFFICIENT K₂, primitive coalbed gas bearing capacity X₀, remaining coal seam gas-bearing capacity X_c, working seam the thickness m of i-th adjacent coal seam adjacent with working seam for thickness m_i, working seam working thickness m₀, adjacent with working seam i-th adjacent coal seam gas bearing capacity X_i, adjacent with working seam i-th adjacent coal seam residual gas content X_ic, adjacent with working seam i-th adjacent coal seam is by the drawing-off gas rate ζ of mining influence_i, exploitation layering height M, working seam tilt angle alpha, exploitation layering bottom coal body mining influence degree of depth h, working seam gas bearing capacity gradient X_t, working seam drawing-off gas rate gradient λ_t, exploitation layering the bottom maximum drawing-off gas rate λ of coal body_max。

First computing module 2, for the data obtained according to data collection computing module 1, calculates the gas emission q of exploitation layering₁；

Second computing module 3, for the data obtained according to data collection computing module 1, calculates the gas emission q of adjacent coal seam adjacent with working seam in stope₂；

3rd computing module 4, for the data obtained according to data collection computing module 1, calculates the bottom gas pressure relief outburst amount q of stope₃；

4th computing module 5, for the result of calculation according to first computing module the 2, second computing module 3 and the 3rd computing module 4, calculates the gas emission q of exploitation work surface_Adopt, wherein q_Adopt=q₁+q₂+q₃。

The prognoses system of the gas emission that invention provides, it is applied in stope working surface of coal mines Forecast of Gas Emission, it is particularly applicable in severe inclined thick coal seam horizontal slice coal mining stope Forecast of Gas Emission, it establishes applicable severe inclined thick coal seam horizontally grouped top-coal drawing mining method face gas prediction of gas emission model, compensate for the deficiency that existing Emit Quantity Prediction Methods In Coal Mines is applied in severe inclined thick coal seam horizontally grouped top-coal drawing mining method face gas prediction of gas emission.

It is preferred that the first computing module 2 adopts equation below to calculate q₁:

$q_1=K_1\×K_2\×\frac{m}{m_0}\×(X_0-X_C).$

It is preferred that the second computing module 3 adopts equation below to calculate q₂:

Wherein n refers to the number of plies of adjacent coal seam, 1≤i≤n.

It is preferred that the 3rd computing module 4 adopts equation below to calculate q₃:

$q_3=\frac{1}{M}{\cos }^{-1}\mathrm{\α}\[(X_0-X_c){\mathrm{\λ}}_{max}h+\frac{X_t{\mathrm{\λ}}_{max}+{\mathrm{\λ}}_t(X_0-X_C)}{2}{h}^2+\frac{X_t{\mathrm{\λ}}_t}{3}{h}^3\].$

To sum up, the Forecasting Methodology of gas emission provided by the invention and prognoses system, it is applied in stope working surface of coal mines Forecast of Gas Emission, it is particularly applicable in severe inclined thick coal seam horizontal slice coal mining stope Forecast of Gas Emission, it establishes applicable severe inclined thick coal seam horizontally grouped top-coal drawing mining method face gas prediction of gas emission model, compensate for the deficiency that existing Emit Quantity Prediction Methods In Coal Mines is applied in severe inclined thick coal seam horizontally grouped top-coal drawing mining method face gas prediction of gas emission.

As required, it is possible to above-mentioned each technical scheme is combined, to reach best-of-breed technology effect.

Above is only principles of the invention and preferred embodiment.It should be pointed out that, for the person of ordinary skill of the art, on the basis of the principle of the invention, it is also possible to make other modification some, also should be regarded as protection scope of the present invention.

Claims (8)

1. the Forecasting Methodology carrying out Forecast of Gas Emission for inclined seam horizontally grouped top-coal drawing mining method work surface, it is characterised in that comprise the steps:

S001: investigate evaluation work face Predicting Gas master data, including:

Country rock Gas COEFFICIENT K₁, work surface lose coal Gas COEFFICIENT K₂, primitive coalbed gas bearing capacity X₀, remaining coal seam gas-bearing capacity X_c, working seam the thickness m of i-th adjacent coal seam adjacent with described working seam for thickness m_i, described working seam working thickness m₀, adjacent with described working seam i-th adjacent coal seam gas bearing capacity X_i, adjacent with described working seam i-th adjacent coal seam residual gas content X_ic, adjacent with described working seam i-th adjacent coal seam is by the drawing-off gas rate ζ of mining influence_i, exploitation layering height M, described working seam tilt angle alpha, exploitation layering bottom coal body mining influence degree of depth h, working seam gas bearing capacity gradient X_t, working seam drawing-off gas rate gradient λ_t, exploitation layering the bottom maximum drawing-off gas rate λ of coal body_max；

S002: calculate the gas emission q of described exploitation layering₁；

S003: calculate the gas emission q of adjacent coal seam adjacent with described working seam in stope₂；

S004: calculate the bottom gas pressure relief outburst amount q of described stope₃；

S005: calculate the gas emission q of exploitation work surface_Adopt, wherein, q_Adopt=q₁+q₂+q₃。

2. Forecasting Methodology according to claim 1, it is characterised in that the gas emission q of described exploitation layering₁Calculated by equation below:

$q_1=K_1\×K_2\×\frac{m}{m_0}\×(X_0-X_C).$

3. Forecasting Methodology according to claim 1, it is characterised in that the gas emission q of described adjacent coal seam₂Calculated by equation below:

$q_2=\underset{i=1}{\overset{n}{\Σ}}\frac{m_i}{m_0}\×{\mathrm{\ζ}}_i(X_i-X_{ic}).$

4. Forecasting Methodology according to claim 1, it is characterised in that bottom gas pressure relief outburst amount q₃Calculated by equation below:

$q_3=\frac{1}{M}{\cos }^{-1}\mathrm{\α}\[(X_0-X_c){\mathrm{\λ}}_{max}h+\frac{X_t{\mathrm{\λ}}_{max}+{\mathrm{\λ}}_t(X_0-X_C)}{2}{h}^2+\frac{X_t{\mathrm{\λ}}_t}{3}{h}^3\].$

5. the prognoses system carrying out Forecast of Gas Emission for inclined seam horizontally grouped top-coal drawing mining method work surface, it is characterised in that including:

Data collection computing module, for collecting and evaluation work face Predicting Gas master data, described data include:

Country rock Gas COEFFICIENT K₁, work surface lose coal Gas COEFFICIENT K₂, primitive coalbed gas bearing capacity X₀, remaining coal seam gas-bearing capacity X_c, working seam the thickness m of i-th adjacent coal seam adjacent with described working seam for thickness m_i, described working seam working thickness m₀, adjacent with described working seam i-th adjacent coal seam gas bearing capacity X_i, adjacent with described working seam i-th adjacent coal seam residual gas content X_ic, adjacent with described working seam i-th adjacent coal seam is by the drawing-off gas rate ζ of mining influence_i, exploitation layering height M, described working seam tilt angle alpha, exploitation layering bottom coal body mining influence degree of depth h, exploitation coal body gas bearing capacity gradient X_t, exploitation coal body drawing-off gas rate gradient λ_t, exploitation layering the bottom maximum drawing-off gas rate λ of coal body_max；

First computing module, for the data obtained according to described data collection computing module, calculates the gas emission q of exploitation layering₁；

Second computing module, for the data obtained according to described data collection computing module, calculates the gas emission q of adjacent coal seam adjacent with described working seam in stope₂；

3rd computing module, for the data obtained according to described data collection computing module, calculates the bottom gas pressure relief outburst amount q of described stope₃；

4th computing module, for the result of calculation according to described first computing module, described second computing module and described 3rd computing module, calculates the gas emission q of exploitation work surface_Adopt, wherein, q_Adopt=q₁+q₂+q₃。

6. prognoses system according to claim 5, it is characterised in that described first computing module adopts equation below to calculate q₁:

$q_1=K_1\×K_2\×\frac{m}{m_0}\×(X_0-X_C).$

7. prognoses system according to claim 5, it is characterised in that described second computing module adopts equation below to calculate q₂:

$q_2=\underset{i=1}{\overset{n}{\Σ}}\frac{m_i}{m_0}\×{\mathrm{\ζ}}_i(X_i-X_{ic}).$

8. prognoses system according to claim 5, it is characterised in that described 3rd computing module adopts equation below to calculate q₃:

$q_3=\frac{1}{M}{\cos }^{-1}\mathrm{\α}\[(X_0-X_c){\mathrm{\λ}}_{max}h+\frac{X_t{\mathrm{\λ}}_{max}+{\mathrm{\λ}}_t(X_0-X_C)}{2}{h}^2+\frac{X_t{\mathrm{\λ}}_t}{3}{h}^3\].$