CN102323998B - Method for checking wind-friction resistance value of mine roadway - Google Patents

Abstract

The invention belongs to the field of mine ventilation in mining engineering, and in particular relates to a method for checking a wind-friction resistance value of a mine roadway. According to the technical scheme, the method comprises the following steps: firstly testing the wind-friction resistance of the mine roadway, and then inputting the wind-friction resistance into an intelligent mine ventilation management system; obtaining a critical path in a mine ventilation network and a high-trust branch in the critical path by means of matrix operation; with the wind-friction resistance value of the branch as a standard, checking the wind-friction resistance value of other roadways so as to obtain the wind-friction resistance value of the roadways which is most approximate to a true value in the whole mine ventilation network; and carrying out confirmatory industrial test in the system so as to ensure the accuracy of the technical scheme. By adopting the method, inversion and simulation of the mine ventilation system on a computer screen are realized on the basis of a PC (personal computer), and the maximum error between the obtained wind-friction resistance value and the field test value of a mine is less than 4%, thus providing scientific basis for ventilation management and technological transformation of the mine and powerful guarantee for safety production of the coal mine.

Description

A kind of method of checking wind-friction resistance value of mine roadway

Technical field

The invention belongs to mining engineering mine ventilation field, be specifically related to a kind of method of checking wind-friction resistance value of mine roadway.

Background technology

Coal is the basic energy resource of China, the safety in production in colliery, matters vital to national well-being and the people's livelihood, along with the attention of country to safe coal production, " first take out afterwards adopt, monitoring, with wind fixed output quota " 12 word policies guide under the safe thing ore deposit of coal significantly contained, but along with the increasing of mining rate, the complexity of mining conditions, security situation still can not look down upon.Mine ventilation is a ring important in " one through three prevention ", and ventilating problem solves bad meeting and directly causes a series of safety problems such as gas or spontaneous combustion.At present for the mine ventilation management computer meanses that adopt more, broken away from the management mode of traditional-handwork and empirical formula, in management level and computational accuracy, there is leaping of matter, yet how accurately to obtain fast mine laneway antitripic wind resistance, become lifting computational accuracy, reach the unified bottleneck problem of mine digitizing and actual conditions.

Distinguished and admirable while doing along Cheng Liudong in roadway, because the formed resistance of friction between the friction between fluid layer and fluid and roadway wall is called specific friction resistance, mainly with the mode of mensuration of Mine Ventilation Resistance (MT/T 440-2008), obtain greatly at present the antitripic wind resistance of mine laneway, but in view of in actual test process, exist operate miss, reading error, unstable airflow and air door the factors such as unlatching restriction and cannot obtain its accurate closure result.In order to realize the accurate closure of data, current employing be that path method or loop method carry out parallel error processing to test result.In actual application, the key of path method is choosing of maximum resistance route, by leading limit on critical path, getting optimum value realizes, but the maximum resistance route resistance of only usining mates as verification condition with whole network drag overall, cannot solve the rationality problem of the wind path weights of respectively connecting in crucial path; And the core of loop method is to guarantee in wind net that any one mesh all meets resistance balance, i.e. ∑ h _i=0, so just cannot reflect the impact on whole wind network system before and after regulating of a certain mesh; In addition, judge that certain wind path function and meaning in wind net is also that network resolves the fundamental purpose with ventilating management, be of value to and analyze wind web frame, and loop method is difficult to obtain this information.

Simple path method or loop method, cannot realize the accurate closure of tunnel antitripic wind resistance, usually appear in mine ventilation system Optimizing Reconstruction and occur that local error is large, resolve scheme and practical operation situation is not inconsistent, or even the problem such as local distinguished and admirable back, for instructing the safety in production of mine to hide some dangers for.

Summary of the invention

The problem existing for prior art, the present invention proposes a kind of method of precise calibration mine laneway antitripic wind resistance, solve the related error of means of testing, instrument existing at present and adopt merely path method, loop method cannot realize roadway antitripic wind resistance accurately closed, and then the deficiency of restriction mine numerical simulation accuracy, realize mine ventilation digital dynamic emulation integrated.

Realizing the technical scheme of the inventive method carries out according to the following steps:

First according to < < mensuration of Mine Ventilation Resistance method > > MT/T 440-2008, measure specific mine ventilation resistance h and the air quantity Q in each airway of ventilating system, and each airway network of actual mine is converted in mine ventilation network data structure input computer, according to ventilation resistance law formula (1), try to achieve Rj; According to ventilation resistance law formula (1), try to achieve each tunnel specific friction resistance R _j;

$h_j=R_j\&CenterDot;{\mathrm{<figure-callout\; id="2"\; label="Q\; j"\; filenames="110927142932.png"\; state="\{\{state\}\}">Q</figure-callout>}}_j^{},(j=\mathrm{1,2},\mathrm{\&Lambda;},n)---\left(1\right)$

Wherein: h _irepresent the ventilation resistance of j bar limit or j branch, i.e. the ventilation resistance in the j tunnel in mine ventilation system, the Pa of unit;

R _jrepresent the windage of j bar limit or j branch, i.e. the specific friction resistance parameter in the j tunnel in mine ventilation system, the Ns of unit ²/ m ⁸;

Q _jrepresent the air quantity of j bar limit or j branch, i.e. the air quantity in the j tunnel in mine ventilation system, the m of unit ³/ s;

By R _jin typing mine ventilation intelligent management system, complete the initial assignment of digitization system data;

Then, adopt depth-first search method traversal mine ventilation network, obtain mine ventilation network independent path matrix P, enumerate and contrast each passage resistance H _iwith overall resistance of mine airflow H _trelative error ε _i;

P＝(p _ij) _s×n，i＝1，2，Λ，s；j＝1，2，Λ，n (2)

$p_{\mathrm{ij}}=\left\{\begin{array}{cc}1,& (e_j\&Element;P_i)\\ 0,& (e_j\&NotElement;P_i)\end{array}\right.---\left(3\right)$

$H_i=\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}{p}_{\mathrm{ij}}\&CenterDot;h_j,i=\mathrm{1,2},\mathrm{\&Lambda;},s---\left(4\right)$

ε _i＝|H _i-H _t|/H _t×100％ (5)

Wherein: s represents the number that forms independent path in mine ventilation system;

N represents tunnel number in mine ventilation system;

P _ij1 value is got in the j tunnel of representative in mine ventilation system on i path time, and 0 value is got in the j tunnel in mine ventilation system not on i path time;

H _irepresent the ventilation resistance on i bar path, the Pa of unit;

H _trepresent overall resistance of mine airflow in mine ventilation system, the Pa of unit;

ε _irepresent the relative error of path ventilation resistance and mine ventilation system drag overall, the % of unit;

Select ε _i≤ 8% path, as crucial way set, is selected the mensuration ventilation resistance h in each tunnel in crucial way set _jfor basis, set W _jbe j bar limit or the j bar branch Coefficients of Approximation in its loop, place, according to circuit pressure equilibrium law formula (6), have:

$W_j=\underset{j=1}{\overset{l}{\mathrm{\&Sigma;}}}{c}_{\mathrm{ij}}{h}_j,(i=\mathrm{1,2},\mathrm{\&Lambda;},n)---\left(6\right)$

Wherein: C is independent circuit matrix;

C _ijfor the element in circuit matrix, 1 value is got in the j tunnel in mine ventilation system on i loop and when wind direction is consistent with loop direction,-1 value is got in the j tunnel in mine ventilation system on i loop and when wind direction and loop opposite direction, and 0 value is got in the j tunnel in mine ventilation system not on i loop time;

By W _jthe best as each tunnel in crucial path is approached set of values, in value set, establishes between k data screening area, and each interval territory is determined by formula (8), W _jat territory X _tfrequency by formula (9), obtained.

$X_t=\left\{\begin{array}{c}{X}_0=W_1+\frac{1}{k}(W_{\max }-W_{\min })\\ X_t=X_{t-1}+\frac{1}{k}(W_{m\mathrm{ax}}-W_{\min })\end{array}\right.,(t=\mathrm{1,2},\mathrm{\&Lambda;},k)---\left(8\right)$

$b_j=\left\{\begin{array}{cc}a=a+1,& W_j\&Element;(X_{t-1},X_t)\\ a=a& W_j\&NotElement;(X_{t-1},X_t)\end{array}\right.,(t=\mathrm{1,2},\mathrm{\&Lambda;},k),(j=\mathrm{1,2},\mathrm{\&Lambda;},l)---\left(9\right)$

Wherein: W _maxrepresent the numerical value that approaches most in each tunnel in crucial path;

W _minrepresent each tunnel in crucial path away from numerical value

X _trepresentative data interval function;

L represents independent loop number, in mine ventilation system, forms the number of independent path;

K represents the data interval number of numerical filter;

B _jrepresent W _jat X _tinterval frequency function;

A representative statistics W _jat territory X _tthe frequency occurring;

According to frequency function b _j, by high-frequency, determine the principle in crucial tunnel in way set, obtain the specific friction resistance R of the highest trust test tunnel i _two-ithereby, determined and highly in the crucial path of mine ventilation network and crucial path trusted branch, as test benchmark, by R _two-iin substitution formula (1), obtain the ventilation resistance h of the tunnel i that this method records _two-i;

With h _two-ifor test benchmark, by its substitution formula (1), calculate the specific friction resistance R in other tunnel;

Each tunnel to mine ventilation system repeats said determination and calculating, for each tunnel, obtain altogether i antitripic wind resistance, choose the windage value that wherein frequency of occurrences is the highest, determine that it is the antitripic wind resistance that each tunnel of whole mine ventilation network approaches actual value most;

By the antitripic wind resistance in each tunnel after checking again typing mine ventilation intelligent management system database, generate the mine ventilation intelligent management system after checking; In mine ventilation intelligent management system after check, mine ventilation system is carried out to confirmatory commerical test, open air door, change and adjust window ara, the specific friction resistance actual value obtaining with the test of mine actual industrial compares, verify precision and accuracy that its check is processed, under the prerequisite tallying with the actual situation in commerical test checking, just can use this system to carry out design and the management of mine ventilation, when the situation that occurs not to be inconsistent with actual conditions, return to above-mentioned steps and re-start check.

Compared with prior art, feature of the present invention and beneficial effect are:

The inventive method is relied on PC computing machine, that has realized Reform of Ventilation System in Mine Shaft is calculating inverting and simulation on screen, means of testing and the related error of instrument of current existence have been solved, and it is accurately closed to adopt merely path method, loop method cannot realize roadway characterisitic parameter, and then the not enough problem of restriction mine numerical simulation accuracy, realize mine ventilation digital dynamic emulation integrated, realized the prediction that air current state and main ventilating fan duty after the ventilating system adjustment such as wind, exploiting field adjustment are adjusted in production distribution planning, system; Come comparative analysis ventilating system to change the effect of front and back, to save construction cost, the engineering time of mine, reduce the destruction that construction nonreversibility is brought, this point, pungent ,Zi Sheng colliery, ground coal ore deposit, ,Hui slope, colliery, ,Cao village, colliery etc. of putting in Huozhou Coal Electricity Group group company is all applied, sunykatuib analysis result and on-the-spot test value maximum relative error are less than 4%, and obtain on-the-spot consistent favorable comment, for mine ventilation management and technological transformation provide scientific basis, for the safety in production of coal mine provides powerful guarantee.

Figure of description

Fig. 1 is the technical scheme Implementation Roadmap of the embodiment of the present invention;

Fig. 2 is embodiment of the present invention Chi Yu work area right flank district ventilation network chart;

Fig. 3 is the ventilation network map after right flank exploiting field, embodiment of the present invention Chi Yu work area is adjusted;

Fig. 4 is each tunnel air quantity error contrast after the embodiment of the present invention is checked;

Fig. 5 is mine ventilation intelligent management system network Diffluence Algorithm flow chart of the present invention.

Embodiment

Embodiment

In view of the complicacy of mine ventilation system network structure, the embodiment of the present invention only be take Shanxi Huozhou Coal Electricity Group group, and pungent to put right flank exploiting field, Chi Yu work area, colliery be example, with reference to the work Implementation Roadmap of Fig. 1, carries out.

The implication of the parameter representative that in the present embodiment, list occurs is as follows:

R _very-i: the true antitripic wind resistance of i bar limit or i branch heading, the i.e. true antitripic wind resistance in the i bar tunnel in mine ventilation system, Ns ²/ m ⁸;

R _survey-i: the testing friction windage value of i bar limit or i branch heading, i.e. the testing friction windage value in the i bar tunnel in mine ventilation system, Ns ²/ m ⁸;

R _logical-i: the check antitripic wind resistance of i bar limit or i branch heading while regulating with path method, the check antitripic wind resistance in the i bar tunnel while regulating with path method in mine ventilation system, Ns ²/ m ⁸;

R _two-i: leading to, the check antitripic wind resistance of going back to binary channel method i bar limit or i branch heading while regulating, i.e. the check antitripic wind resistance in the i bar tunnel in mine ventilation system when leading to, returning double circuit method adjusting, Ns ²/ m ⁸;

Q _logical-i: the air quantity of i bar limit or i branch heading while regulating with path method, the air quantity in the i bar tunnel while regulating with path method in mine ventilation system, m ³/ s;

Q _two-i: leading to, the air quantity that goes back to double circuit method i bar limit or i branch heading while regulating, i.e. the air quantity in the i bar tunnel in mine ventilation system when leading to, returning double circuit method adjusting, m ³/ s;

First the actual ventilation network of the mine laneway in embodiment is converted in mine ventilation network data structure input computer, generates ventilation network as shown in Figure 2; According to < < mensuration of Mine Ventilation Resistance method > > MT/T 440-2008, measure the ventilation resistance h of each airway branch of ventilation network of mine _survey-i;

According to ventilation resistance law formula (1), try to achieve the specific friction resistance R of each tunnel branch _survey-i;

By the R obtaining _survey-iin typing mine ventilation intelligent management system, complete the initial assignment of digitization system data;

Then, adopt depth-first search method traversal mine ventilation network, obtain mine ventilation network independent path matrix P, as formula (2)～(5), enumerate and contrast each passage resistance H _iwith overall resistance of mine airflow H _trelative error ε _i;

Select ε _i≤ 8% path, as crucial way set, is selected the mensuration ventilation resistance h in each tunnel in crucial way set _survey-ifor basis, according to formula (6)～(9), determined and highly in the crucial path of mine ventilation network and crucial path trusted branch, as test benchmark, by R _two-iin substitution formula (1), obtain the ventilation resistance h of the tunnel i that this method records _two-i;

With h _two-ifor test benchmark, by its substitution formula (1), calculate the specific friction resistance R in other nine tunnels;

Each tunnel to mine ventilation system repeats said determination and calculating, for each tunnel, obtain altogether 10 antitripic wind resistances, choose the windage value that wherein frequency of occurrences is the highest, determine that it is the antitripic wind resistance that each tunnel of whole mine ventilation network approaches actual value most, by the antitripic wind resistance in each tunnel after checking again typing mine ventilation intelligent management system database, generate the mine ventilation intelligent management system after checking;

In mine ventilation intelligent management system after check, mine ventilation system is carried out to confirmatory commerical test, at former wind net wind path e ₈add ventilation equipment to control its air quantity, connect e simultaneously ₁₁, form new ventilation network structure as shown in Figure 2, adopt technical solution of the present invention and conventional via method to calculate respectively check, the scheme calculation result obtaining gathers, and in Table 2, each wind path air quantity error is shown in Fig. 3, and Network Ventilation drag overall is in Table 3.

Table 1 wind path characterisitic parameter is checked contrast

Table 2 is checked rear each wind path calculation result contrast

Table 3 Network Ventilation drag overall contrast

From Fig. 3 and table 3, can find out: the technical scheme based on leading to, return double circuit method of the present invention is adjusted to low responsive wind path e by test cumulative errors ₈above, do not affect the calculation accuracy of whole system and each wind path, in contrast to path method check result, although the two all can meet static balancing for whole network, for dynamic production, logical, time double circuit method has more directive significance.

Claims (1)

1. a method for checking wind-friction resistance value of mine roadway, is characterized in that carrying out according to following steps:

First measure ventilation resistance h and the air quantity Q in each airway of mine, and each airway network of actual mine is converted in mine ventilation network data structure input computer, according to ventilation resistance law formula (1), try to achieve R _j;

$h_j=R_j\&CenterDot;Q_j^2,(j=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,n)---\left(1\right)$

Wherein: h _jrepresent the ventilation resistance in the j tunnel in mine ventilation system, the Pa of unit;

R _jrepresent the specific friction resistance in the j tunnel in mine ventilation system, the Ns of unit ²/ m ⁸;

Q _jrepresent the air quantity in the j tunnel in mine ventilation system, the m of unit ³/ s;

By R _jin typing mine ventilation intelligent management system, complete the initial assignment of digitization system data;

Then, adopt depth-first search method traversal mine ventilation network, obtain mine ventilation network independent path matrix P, enumerate and contrast each passage resistance H _iwith overall resistance of mine airflow H _trelative error ε _i;

P=(p _ij) _s×n，i=1，2，…，s；j＝1，2，…，n (2)

$p_{\mathrm{ij}}=\left\{\begin{array}{c}1,(e_j\&Element;P_i)\\ 0,(e_j\&NotElement;P_i)\end{array}\right.---\left(3\right)$

$H_i=\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}{p}_{\mathrm{ij}}\&CenterDot;h_j,i=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,s---\left(4\right)$

ε _i=|H _i-H _t|/H _t×100% （5）

Wherein: s represents the number that forms independent path in mine ventilation system;

N represents tunnel number in mine ventilation system;

P _ij1 value is got in the j tunnel of representative in mine ventilation system on i path time, and 0 value is got in the j tunnel in mine ventilation system not on i path time;

H _irepresent the ventilation resistance on i bar path, the Pa of unit;

H _trepresent overall resistance of mine airflow in mine ventilation system, the Pa of unit;

ε _irepresent the relative error of path ventilation resistance and mine ventilation system drag overall, the % of unit;

Select ε _i≤ 8% path, as crucial way set, is selected the mensuration ventilation resistance h in each tunnel in crucial way set _jfor basis, set W _jbe j bar limit or the j bar branch Coefficients of Approximation in its loop, place, according to circuit pressure equilibrium law formula (6), have:

$W_j=\underset{j=1}{\overset{l}{\mathrm{\&Sigma;}}}{c}_{\mathrm{ij}}{h}_j,(i=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,n)---\left(6\right)$

Wherein: C is independent circuit matrix;

C _ijfor the element in circuit matrix, 1 value is got in the j tunnel in mine ventilation system on i loop and when wind direction is consistent with loop direction,-1 value is got in the j tunnel in mine ventilation system on i loop and when wind direction and loop opposite direction, and 0 value is got in the j tunnel in mine ventilation system not on i loop time;

Select the Coefficients of Approximation W in each tunnel in crucial path _j, make it form optimal approximation coefficient collection, in value set, to establish between k data screening area, each interval territory is determined by formula (8), W _jat territory X _tfrequency by formula (9), obtained,

$X_t=\left\{\begin{array}{c}{X}_0=W_1+\frac{1}{k}(W_{\max }-W_{\min })\\ X_t=X_{t-1}+\frac{1}{k}(W_{\max }-W_{\min })\end{array}\right.,(t=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,k)---\left(8\right)$

$b_j=\left\{\begin{array}{cc}a=a+1,& W_j\&Element;(X_{t-1},X_t)\\ a=a& W_j\&NotElement;(X_{t-1},X_t)\end{array}\right.,(t=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,k),(j=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;,l)---\left(9\right)$

Wherein: W _maxrepresent the numerical value that approaches most in each tunnel in crucial path;

W _minrepresent each tunnel in crucial path away from numerical value;

X _trepresentative data interval function;

L represents independent loop number, in mine ventilation system, forms the number of independent path;

K represents the data interval number of numerical filter;

B _jrepresent W _jat X _tinterval frequency function;

A representative statistics W _jat territory X _tthe frequency occurring;

According to frequency function b _j, by high-frequency, determine the principle in crucial tunnel in way set, obtain the highest trust test tunnel specific friction resistance R _two-ithereby, determined and highly in the crucial path of mine ventilation network and crucial path trusted branch, as test benchmark, by R _two-iin substitution formula (1), obtain the Air Resistance in Roadway Ventilation h that this method records _two-i;

With h _two-ifor test benchmark, by its substitution formula (1), calculate the specific friction resistance R in other tunnel;

Each tunnel to mine ventilation system repeats said determination and calculating, for each tunnel, obtain altogether i antitripic wind resistance, choose the windage value that wherein frequency of occurrences is the highest, determine that it is the antitripic wind resistance that each tunnel of whole mine ventilation network approaches actual value most;

By the antitripic wind resistance in each tunnel after checking again typing mine ventilation intelligent management system database, generate the mine ventilation intelligent management system after checking; In mine ventilation intelligent management system after check, mine ventilation system is carried out to confirmatory commerical test, open air door, change and adjust window ara, the specific friction resistance actual value obtaining with the test of mine actual industrial compares, verify precision and accuracy that its check is processed, under the prerequisite tallying with the actual situation in commerical test checking, just can use this system to carry out design and the management of mine ventilation, when the situation that occurs not to be inconsistent with actual conditions, return to above-mentioned steps and re-start check.

Images