A method and system for predicting and analyzing wind temperature in multiple wind paths of metal mines

Abstract

The present application relates to the field of mining technology, and provides a method and system for predicting and analyzing the wind temperature of multiple wind paths in metal mines. The method establishes a dimensionless model of the surrounding rock temperature field of the tunnel surrounding rock of the metal mine, determines the dimensionless excess temperature of the wall of the tunnel surrounding rock, calculates the unstable heat exchange coefficient between the tunnel surrounding rock and the wind flow of each branch tunnel, and then determines the heat dissipation of the surrounding rock of the branch tunnel; at the same time, calculates the heat dissipation of different types of heat sources in each branch tunnel of the metal mine, constructs a multi-wind path wind temperature prediction model that includes a single branch tunnel wind temperature prediction and a wind node wind temperature prediction under the condition of multiple heat source coupling, and determines the tunnel wind flow temperature of each branch tunnel of the metal mine and the wind flow temperature at the wind node based on the finite difference method. Thus, through the prediction of wind temperature of multiple wind paths in the entire mine, the wind flow temperature prediction under the complex wind network conditions of the metal mine is realized, and the comprehensive influence of the multi-wind path configuration on the underground temperature is determined to establish the optimal wind path configuration.

Description

Technical Field

[0001] The present application relates to the field of mining technology, and in particular to a method and system for predicting and analyzing wind temperatures in multiple wind paths in a metal mine. Background Art

[0002] With the gradual depletion of shallow mineral resources and the increasing demand for resources for economic and social development, the mining of deep mines has become an inevitable trend.

[0003] The increase in mining depth is accompanied by a sharp increase in the temperature of the bottom layer. In addition, the large-scale use of high-power electromechanical equipment has led to a rapid deterioration of the underground thermal environment in recent years, and the high temperature problem encountered during mining has become increasingly serious. Long-term work in a high-temperature environment may cause serious health problems such as skin diseases, heat stroke, and heat exhaustion, threatening the lives of miners.

[0004] In addition, in a high temperature environment, the normal operation of electromechanical equipment is also challenged, not only the equipment stability is reduced, the failure rate is increased, but also the maximum life span is greatly shortened. Therefore, the high temperature problem has become one of the main challenges restricting the safe mining of deep mineral resources. Summary of the invention

[0005] The purpose of this application is to provide a method and system for predicting and analyzing wind temperature in multiple wind paths in a metal mine, so as to solve or alleviate the problems existing in the above-mentioned prior art.

[0006] In order to achieve the above objectives, this application provides the following technical solutions:

[0007] The present application provides a method for predicting and analyzing the wind temperature of multiple wind paths in a metal mine, which is used to predict the wind temperature of multiple wind paths in a metal mine under multi-heat source coupling. The method comprises: step S101, establishing a dimensionless model of the surrounding rock temperature field of the tunnel surrounding rock of the metal mine, and determining the dimensionless excess temperature of the wall of the tunnel surrounding rock of each branch tunnel of the metal mine based on the finite volume method; step S102, calculating the unstable heat transfer coefficient between the tunnel surrounding rock and the wind flow of each branch tunnel according to the dimensionless excess temperature of the wall, so as to determine the heat dissipation of the surrounding rock of the branch tunnel; and, According to the characteristic parameters of each of the branch tunnels, the heat source heat dissipation of different types of heat sources in each branch tunnel of the metal mine is calculated; in step S103, according to the heat source heat dissipation of different types of heat sources in each branch tunnel of the metal mine and the heat dissipation of surrounding rocks, a multi-wind path wind temperature prediction model including a single branch tunnel wind temperature prediction and a combined wind node wind temperature prediction under the condition of multi-heat source coupling is constructed, and based on the finite difference method, according to the tunnel nodes in the ventilation network diagram of the metal mine, the tunnel airflow temperature of each of the branch tunnels of the metal mine and the airflow temperature at the combined wind node are determined.

[0008] Preferably, the dimensionless model of surrounding rock temperature field is:

[0009]

[0010] In the formula, are respectively the dimensionless abscissa and dimensionless ordinate of the surrounding rock temperature field of the tunnel surrounding rock; is the dimensionless ventilation time of the surrounding rock temperature field of the tunnel surrounding rock; is the dimensionless excess temperature of the wall of the surrounding rock of the tunnel; It is the integrated area of ​​the surrounding rock temperature field of the tunnel surrounding rock.

[0011] Preferably, according to the formula:

[0012]

[0013] Determine the unstable heat transfer coefficient between the tunnel surrounding rock and the wind flow of the branch tunnel ;

[0014] In the formula, is the Biot number in the branch lane, is the thermal conductivity of the surrounding rock of the branch tunnel, is the radius of the lane; It is the dimensionless excess temperature of the wall of the surrounding rock of the branch tunnel.

[0015] Preferably, in step S102, according to the formula:

[0016]

[0017] Calculate the heat dissipation of different types of heat sources in the branch tunnel;

[0018] In the formula, is the heat dissipation of the surrounding rock of the branch roadway, is the unstable heat transfer coefficient between the tunnel surrounding rock and the wind flow of the branch tunnel, is the cross-sectional perimeter of the branch tunnel, is the initial cross-sectional area of ​​the branch tunnel, is the average temperature of the wind flow in the branch lane, is the length of the branch lane;

[0019] To release heat for the electromechanical equipment in the branch lanes, is the operating time coefficient of the electromechanical equipment of the branch lane, For the The rated power of the electromechanical equipment; is the total number of the electromechanical equipment; wherein, is a positive integer;

[0020] is the air self-compression heat in the branch tunnel, is the air flow mass flow rate of the branch lane, is the acceleration due to gravity, is the inclination angle of the branch roadway; is the heat dissipation of personnel in the branch tunnel, is the number of labor individuals in the branch lane, is the average heat dissipation of the working individuals in the branch lane; is the heat dissipation of the transported ore in the branch tunnel, is the average specific heat of the ore transported in the branch tunnel, is the ore transport volume in the branch tunnel, is the average temperature of the ore transported in the branch tunnel, is the average temperature of the wind flow in the branch lane;

[0021] is the heat dissipated by oxidation of the tunnel wall of the branch tunnel, is the oxidation heat dissipation coefficient of the wall of the branch tunnel; Heat is released to the groundwater in the branch tunnel, is the specific heat capacity of water, The temperature in the branch lane is higher than The amount of high-temperature water gushing out, is the water outlet temperature of the water inflow port in the branch tunnel; is the heat dissipation of the local ventilator in the branch lane, is the air flow mass flow in the branch lane, is the working wind pressure of the local ventilator in the branch tunnel, is the working efficiency of the motor of the local ventilator, is the heat transfer coefficient of the local ventilator;

[0022] is the heat dissipation of the air inlet shaft of the branch tunnel, For the The inlet air flow temperature of the branch lane is is the height of the air inlet shaft, is the elevation of the air inlet shaft wellhead, where: is a positive integer; The heat released by underground blasting of the branch tunnel is is the heat conversion coefficient of the explosive during blasting, is the explosive mass corresponding to the branch tunnel, The explosion heat of a unit explosive.

[0023] Preferably, the multi-wind path wind temperature prediction model is:

[0024]

[0025] In the formula, is the heat dissipation of the surrounding rock of the branch roadway, is the inclination angle of the branch roadway; is the heat dissipation of other heat sources per unit length; is the air flow mass flow rate of the branch lane; is the acceleration due to gravity;

[0026] Select the wind flow enthalpy value on the micro-element segment for the branch lane changes, is the distance between the branch tunnel and the starting section of the branch tunnel;

[0027] is the wind flow temperature at the wind node of the metal mine, is the number of branch lanes corresponding to the wind converging node, is the specific heat of dry air at constant pressure, for Latent heat of vaporization of water; For the The airflow mass flow rate of the branch lanes is For the The wind flow temperature of the branch lane is For the The moisture content in the branch tunnel described in Article 1, For the The moisture content in the branch tunnel described in Article 1; wherein, is a positive integer, ;

[0028] To release heat for the electromechanical equipment in the branch lanes, is the air self-compression heat in the branch tunnel, is the heat dissipation of personnel in the branch tunnel, is the heat dissipation of the transported ore in the branch tunnel, is the heat dissipated by oxidation of the tunnel wall of the branch tunnel, Heat is released to the groundwater in the branch tunnel, is the heat dissipation of the local ventilator in the branch lane, is the heat dissipation of the air inlet shaft of the branch tunnel, The heat is released for underground blasting of the branch tunnel.

[0029] Preferably, according to the formula:

[0030]

[0031] Determine the wind flow enthalpy in the branch lane ;

[0032] In the formula, is the specific heat of dry air at constant pressure, are the average airflow temperatures in the branch lanes, for The latent heat of vaporization of water, is the specific heat of water vapor at constant pressure; is the air humidity in the branch lane; is the saturated vapor pressure in the branch tunnel; is the distance between the branch tunnel and the starting section of the branch tunnel Relative humidity of wind flow at the location; are the relative humidity of the wind flow at the initial section and the final section of the branch tunnel, respectively. is the length of the branch lane; is the atmospheric pressure in the branch tunnel, It is the elevation below sea level.

[0033] The embodiment of the present application also provides a metal mine multi-wind path wind temperature prediction and analysis system, which is used to predict the wind temperature of multiple wind paths in a metal mine under multi-heat source coupling, and the system includes:

[0034] A temperature field unit is configured to establish a dimensionless model of the surrounding rock temperature field of the tunnel surrounding rock of the metal mine, and determine the dimensionless excess temperature of the wall of the tunnel surrounding rock of each branch tunnel of the metal mine based on the finite volume method;

[0035] The multi-heat source heat calculation unit is configured to calculate the unstable heat transfer coefficient between the tunnel surrounding rock and the wind flow of each branch tunnel according to the dimensionless excess temperature of the wall surface, so as to determine the heat dissipation of the surrounding rock of the branch tunnel; and calculate the heat source heat dissipation of different types of heat sources in each branch tunnel of the metal mine according to the characteristic parameters of each branch tunnel;

[0036] The wind temperature determination unit is configured to construct a multi-wind path wind temperature prediction model including a single branch tunnel wind temperature prediction and a combined wind node wind temperature prediction under multi-heat source coupling conditions according to the heat source heat dissipation of different types of heat sources in each branch tunnel of the metal mine and the heat dissipation of surrounding rocks, and based on the finite difference method, determine the tunnel wind flow temperature of each branch tunnel of the metal mine and the wind flow temperature at the combined wind node according to the tunnel nodes in the ventilation network diagram of the metal mine.

[0037] Beneficial effects:

[0038] The method for predicting and analyzing the wind temperature in multiple wind paths of a metal mine provided in an embodiment of the present application is used to predict the wind temperature in multiple wind paths of a metal mine under multi-heat source coupling. In this method, a dimensionless model of the surrounding rock temperature field of the tunnel surrounding rock of the metal mine is established, and the dimensionless excess temperature of the wall of the tunnel surrounding rock of each branch tunnel of the metal mine is determined based on the finite volume method to calculate the unstable heat transfer coefficient between the tunnel surrounding rock and the wind flow of each branch tunnel, and then determine the heat dissipation of the surrounding rock of the branch tunnel; at the same time, the heat source heat dissipation and the surrounding rock heat dissipation of different types of heat sources in each branch tunnel of the metal mine are calculated, and then a multi-wind path wind temperature prediction model including a single branch tunnel wind temperature prediction and a combined wind node wind temperature prediction under multi-heat source coupling conditions is constructed, and based on the finite difference method, the tunnel wind flow temperature of each branch tunnel of the metal mine and the wind flow temperature at the combined wind node are determined. Therefore, by predicting the wind temperature of multiple wind paths in the whole mine, the wind flow temperature under the complex wind network conditions of the metal mine can be predicted, and the comprehensive impact of multiple wind path configurations on the underground temperature can be determined to establish the optimal wind path configuration, providing a basis for mine refrigeration and cooling; at the same time, according to the wind temperature prediction results, the mutual influence and constraints between the wind paths are clarified, and the high-temperature wind paths are focused on and managed, which promotes the optimal allocation of resources in the mine, helps to reduce the thermal environment hazards of metal mines, and provides an important theoretical basis for the optimal design of metal mine ventilation systems and safe production. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] The drawings constituting part of the present application are used to provide a further understanding of the present application. The exemplary embodiments and descriptions of the present application are used to explain the present application and do not constitute an improper limitation on the present application. Among them:

[0040] Figure 1 A schematic flow chart of a method for predicting and analyzing wind temperature in multiple wind paths of a metal mine provided according to some embodiments of the present application;

[0041] Figure 2 A schematic diagram of a surrounding rock temperature field calculation area of ​​a mine tunnel provided according to some embodiments of the present application;

[0042] Figure 3 A schematic diagram of grid division of a tunnel surrounding rock temperature field calculation area provided according to some embodiments of the present application;

[0043] Figure 4 A schematic diagram of the division of the internal unit control volume of the tunnel surrounding rock temperature field according to some embodiments of the present application;

[0044] Figure 5 A schematic diagram of the division of boundary unit control volumes of the tunnel surrounding rock temperature field provided according to some embodiments of the present application;

[0045] Figure 6 A schematic diagram of a mine ventilation network provided according to some embodiments of the present application;

[0046] Figure 7 A schematic diagram of node division for single lane wind temperature prediction provided according to some embodiments of the present application;

[0047] Figure 8 A schematic diagram of branch lane wind flows converging at a wind converging node according to some embodiments of the present application;

[0048] Fig. 9 A schematic diagram of a ventilation network of a metal mine provided according to some embodiments of the present application;

[0049] Fig.10 for Fig. 9 A comparison diagram of the predicted and measured airflow temperatures at the starting nodes of each lane in the ventilation network provided in the embodiment;

[0050] Fig.11 for Fig. 9 A comparison diagram of the predicted and measured temperatures of the starting points of each lane in the ventilation network provided in the embodiment;

[0051] Fig.12 for Fig. 9 A comparison diagram of the predicted and measured temperatures of the wind flow at each wind-combining node in the ventilation network provided in the embodiment;

[0052] Fig.13 This is a schematic diagram of the structure of a metal mine multi-wind path wind temperature prediction and analysis system provided according to some embodiments of the present application. DETAILED DESCRIPTION

[0053] The present application will be described in detail below with reference to the accompanying drawings and in conjunction with embodiments. Each example is provided by way of explanation of the present application and does not limit the present application. In fact, it will be clear to those skilled in the art that modifications and variations may be made in the present application without departing from the scope or spirit of the present application. For example, a feature shown or described as a part of an embodiment may be used in another embodiment to produce yet another embodiment. All other embodiments obtained by those of ordinary skill in the art based on the embodiments in the embodiments of the present invention should belong to the scope of protection of the embodiments of the present invention.

[0054] In order to reduce the temperature of the mine, a large amount of money needs to be invested to build a large-scale refrigeration and cooling system, which greatly increases the production cost and makes it difficult for small and medium-sized enterprises to afford it. Therefore, in order to reduce the impact of extreme high temperatures underground and improve the safety of deep mine mining, multi-wind path wind temperature prediction analysis can not only determine the comprehensive impact of multi-wind path configuration on underground temperature and find the optimal wind path configuration plan, but also clarify the mutual influence and restriction relationship between each wind path based on the prediction results, focus on and manage the high-temperature wind path, and effectively promote the optimal allocation of mine resources.

[0055] At present, due to the long ventilation lines, large ventilation resistance, and many electromechanical equipment in metal mines, high temperature heat damage has occurred underground. Especially in summer, the surface temperature is high, which leads to the continuous rise in the temperature of underground workplaces and the degree of heat damage. The determination of the multi-wind path ventilation wind temperature prediction and analysis method provides a basis for the refrigeration and cooling of mines, helps to reduce the thermal environment hazards of metal mines, and has important theoretical and practical significance for the optimization design of metal mine ventilation systems and safe production. It not only provides a scientific basis for the adoption of cooling measures and the design of cooling solutions, but also provides guidance for the precise placement of the apex of underground high temperature heat damage, which can maximize the guarantee of underground safe production.

[0056] However, the complex on-site wind network environment restricts the development of underground multi-wind path wind temperature prediction technology. The existing mine tunnel wind temperature prediction and analysis methods have limitations. Usually, only the wind temperature of a single branch tunnel is predicted, ignoring the influence of factors such as wind flow separation and intersection, resulting in low efficiency in heat damage control and a large amount of resource waste. It is difficult to adapt to the wind temperature prediction needs of the complex wind network environment of mines.

[0057] Based on this, the embodiment of the present application proposes a method for predicting and analyzing the wind temperature of multiple wind paths in metal mines, establishes a multi-wind path wind temperature prediction model under the conditions of multiple heat source coupling, which includes two parts: single branch tunnel wind temperature prediction and combined wind node wind temperature prediction. The airflow temperature of each tunnel is obtained based on the finite difference method, and the wind temperature prediction of multiple wind paths in metal mines under the conditions of multiple heat source coupling is realized, which effectively reduces the wind temperature prediction error in complex wind network environments and improves the accuracy and reliability of wind temperature prediction in complex tunnels of metal mines.

[0058] like Figures 1 to 8 As shown, the metal mine multi-wind path wind temperature prediction and analysis method for predicting the wind temperature of multiple wind paths in a metal mine under multi-heat source coupling includes:

[0059] Step S101: Establish a dimensionless model of the surrounding rock temperature field of the tunnel surrounding rock of the metal mine, and determine the dimensionless excess temperature of the wall of the tunnel surrounding rock of each branch tunnel of the metal mine based on the finite volume method.

[0060] In this application, the integral expression of the dimensionless model of the surrounding rock temperature field of the tunnel surrounding rock of the metal mine is constructed:

[0061]

[0062] In the formula, They are the dimensionless abscissa and dimensionless ordinate of the surrounding rock temperature field of the tunnel surrounding rock respectively; is the dimensionless ventilation time of the surrounding rock temperature field of the tunnel surrounding rock; It is the dimensionless excess temperature of the wall of the tunnel surrounding rock; It is the integral area of ​​the surrounding rock temperature field of the tunnel surrounding rock.

[0063] The integral expression of the dimensionless model of the tunnel surrounding rock temperature field is simplified by Green's formula:

[0064]

[0065] In the formula, It is the perimeter of the integral area of ​​the surrounding rock temperature field of the tunnel surrounding rock.

[0066] When determining the dimensionless excess temperature of the tunnel surrounding rock wall through the integral expression of the dimensionless model of the tunnel surrounding rock temperature field, the calculation area of ​​the surrounding rock temperature field must be determined first. ,in, represents the convection heat transfer boundary between the tunnel wall and the wind flow, It represents the distal isothermal boundary inside the surrounding rock, where the boundary temperature is equal to the original rock temperature.

[0067] After determining the calculation area of ​​the surrounding rock temperature field, the calculation area of ​​the surrounding rock temperature field is meshed with triangular units as the basic units. At the tunnel wall, the heat exchange between the surrounding rock and the wind flow is intense. The geometric increase method is used to increase the unit density around the tunnel wall to effectively improve the calculation accuracy and efficiency.

[0068] After the unit grid is divided, the node heat balance method of the internal grid node unit in the surrounding rock temperature field calculation area is established. , which is associated with the six surrounding triangular elements, so the node The heat balance equation of is the sum of the heat contributions of the six associated units. Triangular unit The center of gravity, point Side The midpoint of the triangle unit For example, a quadrilateral Node In the triangle unit The sub-control area in the quadrilateral The area of ​​the triangular unit Furthermore, the polygon Finally, the node The total control area in all its associated units. Therefore, nodes within the total control area can be established The heat balance equation is:

[0069]

[0070] In the formula, For the Sub-control region to node The contribution to the energy equation, Respectively The corresponding boundary of the sub-control area is The projection length in the direction, For the The area of ​​the sub-control region.

[0071] For triangular elements After analysis, we know that the node 、 node 、 node The temperatures at , the coordinates are , triangular unit The area is , then the triangular unit The temperature at any location inside is:

[0072]

[0073] in,

[0074]

[0075] Node 、 node 、 node The dimensionless excess temperature at .

[0076] Furthermore, based on Fourier's law, the triangular unit Inner edge The heat flux densities in the directions are:

[0077]

[0078] With triangular element For Node Contribution for:

[0079]

[0080] The same triangle unit can be obtained For Node ,node Contribution . The triangular element For Node ,node ,node Contribution In matrix form:

[0081]

[0082] in,

[0083]

[0084] When the node When the node is located at the boundary of the computational region, the heat flux passing through the boundary will also affect the node The boundary is also expressed in matrix form. For Node ,node The contribution of the heat balance equation is:

[0085]

[0086] in,

[0087]

[0088] Then, in the calculation area, all the associated triangle element pairs of nodes are calculated in turn by the element search method. Contributions to the heat balance equation are superimposed and nodes are established The complete form of the heat balance equation at . Similarly, repeat the above operation for all nodes in the calculation area to construct a heat balance linear equation system with all node temperatures as unknown quantities:

[0089]

[0090] In the formula, are the dimensionless temperatures of each node, are the elements of the principal coefficient matrix, is the non-steady-state term of the dimensionless temperature change of each node with time, are the elements of the non-stationary coefficient matrix, is the constant term matrix. Finally, the backward difference format is used to convert The simplified equation is solved by Gaussian elimination method, and the dimensionless excess temperature of the surrounding rock wall is finally obtained. .

[0091] Step S102: Calculate the unstable heat transfer coefficient between the tunnel surrounding rock and the wind flow of each branch tunnel according to the dimensionless excess temperature of the wall to determine the heat dissipation of the surrounding rock of the branch tunnel; and calculate the heat source heat dissipation of different types of heat sources in each branch tunnel of the metal mine according to the characteristic parameters of each branch tunnel.

[0092] After obtaining the dimensionless excess temperature of the surrounding rock wall, according to the formula:

[0093]

[0094] Determine the unstable heat transfer coefficient between the tunnel surrounding rock and the air flow in the branch tunnel In the formula, is the Biot number in the branch lane, is the thermal conductivity of the surrounding rock of the branch tunnel, is the radius of the roadway. According to the formula:

[0095]

[0096] Determine the Biot number in the branch lane In the formula, is the convective heat transfer coefficient of the branch roadway, It is the characteristic length of the branch tunnel, which is equal to 4 times the ratio of the cross-sectional area of ​​the branch tunnel to its perimeter.

[0097] In this application, the analysis of the change of air flow temperature in the mine is realized through parameters such as the air path type, tunnel basic parameters, original rock temperature, and thermodynamic properties of air in the metal mine tunnel. In the internal environment of the mine, the thermodynamic parameters of the air flow will continue to change due to the combined influence of multiple factors such as geological conditions, the design and operation status of the ventilation system, and the mining activities of the working face. Here, by analyzing the thermodynamic parameters related to the wet air in the mine, an accurate and comprehensive basis is provided for the change of the air flow temperature in the mine. Among them, the thermodynamic parameters include: ventilation time, tunnel circumference, tunnel area, tunnel length, thermal conductivity of surrounding rock, convection heat transfer coefficient, original rock temperature of surrounding rock, initial temperature of air flow, initial humidity of air flow, air volume, etc.

[0098] The wind flow path in the mine is not a combination of all possible paths, but is limited by the mine structure and ventilation system design. The ventilation network diagram of the mine reflects the separation and combination relationship of the wind flow in the mine and the mutual connection between the wind flow direction and branches, and describes the connection between the wind flow propagation path and branches in the mine ventilation system. By constructing a ventilation network diagram under the conditions of multi-wind path ventilation in metal mines, the temperature rise caused by various heat sources in the mine during the flow of the mine wind in the tunnel is determined. Specifically, according to the formula:

[0099]

[0100] Determine the heat dissipation of the surrounding rock in the branch tunnel In the formula, is the unstable heat transfer coefficient between the tunnel surrounding rock and the air flow of the branch tunnel, is the cross-sectional perimeter of the branch tunnel, is the original rock temperature of the surrounding rock of the branch tunnel, is the average temperature of the wind flow in the branch lane, is the length of the branch lane.

[0101] At the same time, according to the characteristic parameters (thermodynamic parameters) of each branch tunnel, the heat dissipation of different types of heat sources in each branch tunnel of the metal mine is calculated. Specifically, according to the formula:

[0102]

[0103] Calculate the heat dissipation of different types of heat sources in the branch tunnel. To release heat for the electromechanical equipment in the branch lanes. is the operating time coefficient of the electromechanical equipment in the branch lane, For the Rated power of each electromechanical device; is the total number of electromechanical equipment. is the air self-compression heat in the branch tunnel, is the air flow mass flow rate of the branch lane, is the acceleration due to gravity, is the inclination angle of the branch tunnel. The heat dissipation for personnel in the branch tunnels, is the number of labor individuals in the branch lane, is the average heat dissipation of working individuals in the branch lane.

[0104] is the heat dissipation of the transported ore in the branch tunnel, is the average specific heat of ore transported in the branch tunnel, is the ore transport volume in the branch tunnel, is the average temperature of the ore transported in the branch tunnel, is the average temperature of the wind flow in the branch tunnel. is the heat dissipated by oxidation on the wall of the branch tunnel, is the oxidation heat dissipation coefficient of the wall of the branch tunnel. Release heat to the groundwater in the branch tunnels. is the specific heat capacity of water, The temperature in the branch tunnel is higher than The amount of high-temperature water gushing out, It is the water outlet temperature of the water outlet in the branch tunnel. is the heat dissipation of the local ventilator in the branch tunnel, is the air flow mass flow in the branch lane, is the working wind pressure of the local ventilator in the branch tunnel, is the motor efficiency of the local ventilator, is the heat transfer coefficient of the local ventilator. is the heat dissipation of the air inlet shaft of the branch tunnel, For the The inlet air flow temperature of the branch lanes, is the height of the air inlet shaft, It is the elevation of the air inlet shaft mouth. The heat released by underground blasting of branch tunnels, is the heat conversion coefficient of the explosive during blasting, is the explosive mass corresponding to the branch tunnel, The explosion heat of a unit explosive.

[0105] Step S103: According to the heat source heat dissipation of different types of heat sources in each branch tunnel of the metal mine, a multi-wind path wind temperature prediction model is constructed, which includes a single branch tunnel wind temperature prediction and a combined wind node wind temperature prediction under multi-heat source coupling conditions; and based on the finite difference method, according to the unstable heat transfer coefficient and the tunnel nodes in the ventilation network diagram of the metal mine, the tunnel airflow temperature of each branch tunnel of the metal mine and the airflow temperature at the combined wind node are determined.

[0106] Combining the heat dissipation of different types of heat sources in each branch tunnel of a metal mine, the heat dissipation of other heat sources per unit length of the gas heat source in the branch tunnel except for the heat dissipation of the surrounding rock can be obtained, that is:

[0107]

[0108] Then, according to the unstable heat transfer coefficient between the tunnel surrounding rock and the air flow in the metal mine , heat dissipation of other heat sources per unit length and other thermodynamic parameters in the branch tunnel, arbitrarily select a branch tunnel micro-element segment dx , construct the heat balance equation for predicting the wind temperature in a single branch lane:

[0109]

[0110] In the formula, Select the wind flow enthalpy value on the micro-element segment for the branch roadway changes, is the inclination angle of the branch roadway; is the air flow mass flow rate of the branch lane; is the acceleration due to gravity. According to the formula:

[0111]

[0112] Determine the enthalpy of air flow In the formula, is the specific heat of dry air at constant pressure, are the average airflow temperatures in the branch lanes, for The latent heat of vaporization of water, is the specific heat of water vapor at constant pressure; is the air humidity content in the branch tunnel; is the saturated vapor pressure in the branch tunnel; The distance between the branch tunnel and the starting section of the branch tunnel Relative humidity of wind flow at the location; are the relative humidity of the wind flow at the beginning and end sections of the branch tunnel, is the length of the branch lane; is the atmospheric pressure in the branch tunnel, It is the elevation below sea level.

[0113] Then, the wind temperature prediction equation of a single branch roadway is discretized according to the finite difference method, and the length of the branch roadway is divided into The lane is divided into several nodes, and the distance between adjacent nodes is , and replace the differential quotient with the difference quotient, in any The discrete finite difference prediction equation of the wind temperature in a single supporting tunnel is:

[0114]

[0115] In the formula, For the The node number of the branch lane.

[0116] The airflow in the lane is dispersed to multiple branch lanes through the air distribution node. The air volume changes while the temperature and humidity remain unchanged. The air distribution node only affects the air volume distribution; the wind converging node gathers the airflow, and both the air volume and temperature change. According to the wind temperature prediction results of each branch lane, according to the formula:

[0117]

[0118] Determine the air flow temperature at the wind node in a metal mine In the formula, is the number of branch lanes corresponding to the wind converging node, is the specific heat of dry air at constant pressure, for Latent heat of vaporization of water; For the The airflow mass flow rate of the branch lanes, For the The air flow temperature of the branch lanes, For the The moisture content in the branch tunnel, For the The moisture content in the branch tunnel; Is a positive integer.

[0119] To release heat for the electromechanical equipment in the branch lanes. is the air self-compression heat in the branch tunnel, The heat dissipation for personnel in the branch tunnels, is the heat dissipation of the transported ore in the branch tunnel, is the heat dissipated by oxidation on the wall of the branch tunnel, Release heat to the groundwater in the branch tunnels. is the heat dissipation of the local ventilator in the branch tunnel, is the heat dissipation of the air inlet shaft of the branch tunnel, It releases heat for underground blasting of branch tunnels.

[0120] In a specific example, the wind flow in branch lane 1 and branch lane 2 converges to branch lane 3 through the wind node; the wind flow temperature in branch lane 1 and branch lane 2 can be obtained by the wind temperature prediction equation of a single branch lane. The wind flow in branch lane 1 and branch lane 2 follows the mass conservation before and after mixing, that is:

[0121]

[0122] The wind flow follows the law of energy conservation before and after mixing:

[0123]

[0124] The air flow temperature of branch tunnel 3 (air converging node) is:

[0125]

[0126] in, are the air flow mass flow at nodes ①, ②, and ③ respectively, Wind flow enthalpy values ​​at nodes ①, ②, and ③ respectively.

[0127] In another specific example, Figures 9 to 12 As shown in the figure, there are 487 branch tunnels and 314 nodes in the metal mine, and the ventilation network diagram of the metal mine is drawn. Then, the simulation data and actual measurement data of the starting nodes of each branch tunnel are selected for analysis, and the predicted temperature and actual temperature of the wind flow of each wind node in the ventilation network are compared; the wind flow flows in from the 1# wind shaft and then is diverted to each branch tunnel.

[0128] In this embodiment, the relevant parameters such as the type of mine tunnel wind path, tunnel basic parameters, original rock temperature, and air thermodynamic properties are first determined and predicted, and a dimensionless mathematical model of the tunnel surrounding rock temperature field is established. The unstable heat exchange coefficient between the surrounding rock and the wind flow is solved based on the finite volume method; the mine ventilation network diagram is drawn, and the heat dissipation of various heat sources in each branch tunnel of the mine is calculated respectively. Then, a multi-wind path wind temperature prediction model under multi-heat source coupling is established, and discretization and settlement are performed based on the finite difference method.

[0129] Therefore, by predicting the wind temperature of multiple wind paths in the whole mine, the wind flow temperature under the complex wind network conditions of the metal mine can be predicted, and the comprehensive impact of multiple wind path configurations on the underground temperature can be determined to establish the optimal wind path configuration, providing a basis for mine refrigeration and cooling; at the same time, according to the wind temperature prediction results, the mutual influence and constraints between the wind paths are clarified, and the high-temperature wind paths are focused on and managed, which promotes the optimal allocation of resources in the mine, helps to reduce the thermal environment hazards of metal mines, and provides an important theoretical basis for the optimal design of metal mine ventilation systems and safe production.

[0130] like Fig.13 As shown, the embodiment of the present application also provides a metal mine multi-wind path wind temperature prediction and analysis system, which is used to predict the multi-wind path wind temperature of a metal mine under multi-heat source coupling, and the system includes:

[0131] The temperature field unit 1301 is configured to establish a dimensionless model of the surrounding rock temperature field of the tunnel surrounding rock of the metal mine, and determine the dimensionless excess temperature of the wall of the tunnel surrounding rock of each branch tunnel of the metal mine based on the finite volume method;

[0132] The multi-heat source heat calculation unit 1302 is configured to calculate the unstable heat transfer coefficient between the tunnel surrounding rock and the wind flow of each branch tunnel according to the dimensionless excess temperature of the wall surface, so as to determine the heat dissipation of the surrounding rock of the branch tunnel; and calculate the heat source heat dissipation of different types of heat sources in each branch tunnel of the metal mine according to the characteristic parameters of each branch tunnel;

[0133] The wind temperature determination unit 1303 is configured to construct a multi-wind path wind temperature prediction model including a single branch tunnel wind temperature prediction and a combined wind node wind temperature prediction under multi-heat source coupling conditions according to the heat source heat dissipation of different types of heat sources in each branch tunnel of the metal mine and the heat dissipation of surrounding rocks, and based on the finite difference method, determine the tunnel wind flow temperature of each branch tunnel of the metal mine and the wind flow temperature at the combined wind node according to the tunnel nodes in the ventilation network diagram of the metal mine.

[0134] The metal mine multi-wind path wind temperature prediction and analysis system provided in the embodiment of the present application can implement the steps and processes of the metal mine multi-wind path wind temperature prediction and analysis method of any of the above-mentioned embodiments, and achieve the same technical effects, which will not be repeated here one by one.

[0135] In the description of the present invention, it is to be understood that the terms “center”, “longitudinal”, “lateral”, “length”, “width”, “thickness”, “up”, “down”, “front”, “back”, “left”, “right”, “vertical”, “horizontal”, “top”, “bottom”, “inside”, “outside”, “clockwise”, “counterclockwise”, “axial”, “radial”, “circumferential”, etc., indicating orientations or positional relationships based on the orientations or positional relationships shown in the accompanying drawings, are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the referred device or element must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be understood as limiting the present invention.

[0136] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined as "first" and "second" may explicitly or implicitly include at least one of the features. In the description of the present invention, the meaning of "plurality" is at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0137] In the present invention, unless otherwise clearly specified and limited, the terms "installed", "connected", "connected", "fixed" and the like should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral one; it can be a mechanical connection, an electrical connection, or communication with each other; it can be a direct connection, or an indirect connection through an intermediate medium, it can be the internal connection of two elements or the interaction relationship between two elements, unless otherwise clearly defined. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0138] In the present invention, unless otherwise clearly specified and limited, a first feature being "above" or "below" a second feature may mean that the first and second features are in direct contact, or the first and second features are in indirect contact through an intermediate medium. Moreover, a first feature being "above", "above" or "above" a second feature may mean that the first feature is directly above or obliquely above the second feature, or simply means that the first feature is higher in level than the second feature. A first feature being "below", "below" or "below" a second feature may mean that the first feature is directly below or obliquely below the second feature, or simply means that the first feature is lower in level than the second feature.

[0139] In the present invention, the terms "one embodiment", "some embodiments", "examples", "specific examples", or "some examples" mean that the specific features, structures, materials, or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described may be combined in any one or more embodiments or examples in a suitable manner.

[0140] The above description is only a preferred embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for predicting and analyzing wind temperature in multiple wind paths of a metal mine, characterized in that: The method is used to predict the wind temperature of multiple wind paths in a metal mine under multi-heat source coupling, and comprises: Step S101, establishing a dimensionless model of the surrounding rock temperature field of the tunnel surrounding rock of the metal mine, and determining the dimensionless excess temperature of the surrounding rock wall of each branch tunnel of the metal mine based on the finite volume method; Step S102, calculating the unstable heat exchange coefficient between the tunnel surrounding rock and the wind flow of each branch tunnel according to the dimensionless excess temperature of the wall surface, so as to determine the heat dissipation of the surrounding rock of the branch tunnel; And, according to the characteristic parameters of each branch tunnel, the heat source heat dissipation of different types of heat sources in each branch tunnel of the metal mine is calculated; wherein, according to the formula: Calculate the heat dissipation of different types of heat sources in the branch tunnel; In the formula, is the heat dissipation of the surrounding rock of the branch roadway, is the unstable heat transfer coefficient between the tunnel surrounding rock and the wind flow of the branch tunnel, is the cross-sectional perimeter of the branch tunnel, is the original rock temperature of the surrounding rock of the branch tunnel, is the average temperature of the wind flow in the branch lane, is the length of the branch lane; To release heat for the electromechanical equipment in the branch lanes, is the operating time coefficient of the electromechanical equipment of the branch lane, For the The rated power of the electromechanical equipment; is the total number of the electromechanical equipment; wherein, is a positive integer; is the air self-compression heat in the branch tunnel, is the air flow mass flow rate of the branch lane, is the acceleration due to gravity, is the inclination angle of the branch roadway; is the heat dissipation of personnel in the branch tunnel, is the number of labor individuals in the branch lane, is the average heat dissipation of the working individuals in the branch lane; is the heat dissipation of the transported ore in the branch tunnel, is the average specific heat of the ore transported in the branch tunnel, is the ore transport volume in the branch tunnel, is the average temperature of the ore transported in the branch tunnel, is the average temperature of the wind flow in the branch lane; is the heat dissipated by oxidation of the tunnel wall of the branch tunnel, is the oxidation heat dissipation coefficient of the wall surface of the branch tunnel; Heat is released to the groundwater in the branch tunnel, is the specific heat capacity of water, The temperature in the branch lane is higher than The amount of high-temperature water gushing out, is the water outlet temperature of the water inflow port in the branch tunnel; is the heat dissipation of the local ventilator in the branch lane, is the working wind pressure of the local ventilator in the branch tunnel, is the working efficiency of the motor of the local ventilator, is the heat transfer coefficient of the local ventilator; is the heat dissipation of the air inlet shaft of the branch tunnel, For the The inlet air flow temperature of the branch lane is is the height of the air inlet shaft, is the elevation of the air inlet wellhead; wherein, is a positive integer; The heat released by underground blasting of the branch tunnel is is the heat conversion coefficient of the explosive during blasting, is the explosive mass corresponding to the branch tunnel, is the explosion heat per unit of explosive; Step S103: According to the heat source heat dissipation of different types of heat sources in each branch tunnel of the metal mine and the heat dissipation of surrounding rocks, a multi-wind path wind temperature prediction model including a single branch tunnel wind temperature prediction and a wind node wind temperature prediction under the condition of multi-heat source coupling is constructed, and based on the finite difference method, according to the tunnel nodes in the ventilation network diagram of the metal mine, the tunnel wind flow temperature of each branch tunnel of the metal mine and the wind flow temperature at the wind node are determined; wherein, the multi-wind path wind temperature prediction model is: In the formula, is the inclination angle of the branch roadway; is the heat dissipation of other heat sources per unit length; is the acceleration due to gravity; Select the wind flow enthalpy value on the micro-element segment for the branch lane changes, is the distance between the branch tunnel and the starting section of the branch tunnel; is the wind flow temperature at the wind node of the metal mine, is the number of branch lanes corresponding to the wind converging node, is the specific heat of dry air at constant pressure, for Latent heat of vaporization of water; For the The airflow mass flow rate of the branch lanes is For the The wind flow temperature of the branch lane is For the The moisture content in the branch tunnel described in Article 1, For the The moisture content in the branch tunnel described in Article 1; wherein, is a positive integer, .

2. The method for predicting and analyzing wind temperature in multiple wind paths of a metal mine according to claim 1, characterized in that: The dimensionless model of surrounding rock temperature field is: In the formula, are respectively the dimensionless abscissa and dimensionless ordinate of the surrounding rock temperature field of the tunnel surrounding rock; is the dimensionless ventilation time of the surrounding rock temperature field of the tunnel surrounding rock; is the dimensionless excess temperature of the wall of the surrounding rock of the tunnel; It is the integrated area of ​​the surrounding rock temperature field of the tunnel surrounding rock.

3. The method for predicting and analyzing wind temperature in multiple wind paths of a metal mine according to claim 2, characterized in that: According to the formula: Determine the unstable heat transfer coefficient between the tunnel surrounding rock and the wind flow of the branch tunnel ; In the formula, is the thermal conductivity of the surrounding rock of the branch tunnel, is the radius of the lane; is the dimensionless excess temperature of the wall of the surrounding rock of the branch tunnel; is the Biot number in the branch lane.

4. The method for predicting and analyzing wind temperature in multiple wind paths of a metal mine according to claim 1, characterized in that: According to the formula: Determine the wind flow enthalpy in the branch lane ; In the formula, is the specific heat of dry air at constant pressure, are the average airflow temperatures in the branch lanes, for The latent heat of vaporization of water, is the specific heat of water vapor at constant pressure; is the air humidity in the branch lane; is the saturated vapor pressure in the branch tunnel; is the distance between the branch tunnel and the starting section of the branch tunnel Relative humidity of wind flow at the location; are the relative humidity of the wind flow at the initial section and the final section of the branch tunnel, respectively. is the length of the branch lane; is the atmospheric pressure in the branch tunnel, It is the elevation below sea level.

5. A metal mine multi-wind path wind temperature prediction and analysis system, characterized in that: The system is used for predicting wind temperature in multiple wind paths of a metal mine under multi-heat source coupling, and comprises: A temperature field unit is configured to establish a dimensionless model of the surrounding rock temperature field of the tunnel surrounding rock of the metal mine, and determine the dimensionless excess temperature of the wall of the tunnel surrounding rock of each branch tunnel of the metal mine based on the finite volume method; The multi-heat source heat calculation unit is configured to calculate the unstable heat transfer coefficient between the tunnel surrounding rock and the wind flow of each branch tunnel according to the dimensionless excess temperature of the wall surface, so as to determine the heat dissipation of the surrounding rock of the branch tunnel; and calculate the heat source heat dissipation of different types of heat sources in each branch tunnel of the metal mine according to the characteristic parameters of each branch tunnel; wherein, according to the formula: Calculate the heat dissipation of different types of heat sources in the branch tunnel; In the formula, is the heat dissipation of the surrounding rock of the branch roadway, is the unstable heat transfer coefficient between the tunnel surrounding rock and the wind flow of the branch tunnel, is the cross-sectional perimeter of the branch tunnel, is the original rock temperature of the surrounding rock of the branch tunnel, is the average temperature of the wind flow in the branch lane, is the length of the branch lane; To release heat for the electromechanical equipment in the branch lanes, is the operating time coefficient of the electromechanical equipment of the branch lane, For the The rated power of the electromechanical equipment; is the total number of the electromechanical equipment; wherein, is a positive integer; is the air self-compression heat in the branch tunnel, is the air flow mass flow rate of the branch lane, is the acceleration due to gravity, is the inclination angle of the branch roadway; is the heat dissipation of personnel in the branch tunnel, is the number of labor individuals in the branch lane, is the average heat dissipation of the working individuals in the branch lane; is the heat dissipation of the transported ore in the branch tunnel, is the average specific heat of the ore transported in the branch tunnel, is the ore transport volume in the branch tunnel, is the average temperature of the ore transported in the branch tunnel, is the average temperature of the wind flow in the branch lane; is the heat dissipated by oxidation of the tunnel wall of the branch tunnel, is the oxidation heat dissipation coefficient of the wall surface of the branch tunnel; Heat is released to the groundwater in the branch tunnel, is the specific heat capacity of water, The temperature in the branch lane is higher than The amount of high-temperature water gushing out, is the water outlet temperature of the water inflow port in the branch tunnel; is the heat dissipation of the local ventilator in the branch lane, is the working wind pressure of the local ventilator in the branch tunnel, is the working efficiency of the motor of the local ventilator, is the heat transfer coefficient of the local ventilator; is the heat dissipation of the air inlet shaft of the branch tunnel, For the The inlet air flow temperature of the branch lane is is the height of the air inlet shaft, is the elevation of the air inlet wellhead; wherein, is a positive integer; The heat released by underground blasting of the branch tunnel is is the heat conversion coefficient of the explosive during blasting, is the explosive mass corresponding to the branch tunnel, is the explosion heat per unit of explosive; The wind temperature determination unit is configured to construct a multi-wind path wind temperature prediction model including a single branch roadway wind temperature prediction and a wind node wind temperature prediction under a multi-heat source coupling condition according to the heat source heat dissipation of different types of heat sources in each branch roadway of the metal mine and the heat dissipation of surrounding rocks, and based on the finite difference method, determine the roadway wind flow temperature of each branch roadway of the metal mine and the wind flow temperature at the wind node according to the roadway nodes in the ventilation network diagram of the metal mine; wherein the multi-wind path wind temperature prediction model is: In the formula, is the inclination angle of the branch roadway; is the heat dissipation of other heat sources per unit length; is the acceleration due to gravity; Select the wind flow enthalpy value on the micro-element segment for the branch lane changes, is the distance between the branch tunnel and the starting section of the branch tunnel; is the wind flow temperature at the wind node of the metal mine, is the number of branch lanes corresponding to the wind converging node, is the specific heat of dry air at constant pressure, for Latent heat of vaporization of water; For the The airflow mass flow rate of the branch lanes is For the The wind flow temperature of the branch lane is For the The moisture content in the branch tunnel described in Article 1, For the The moisture content in the branch tunnel described in Article 1; wherein, is a positive integer, .

Images