CN-121765886-B - High-precision network resolving method based on mass-volume flow double coupling

CN121765886BCN 121765886 BCN121765886 BCN 121765886BCN-121765886-B

Abstract

The invention discloses a high-precision network resolving method based on mass-volume flow double coupling, which relates to the technical field of mine ventilation and adopts Kruskal algorithm to select minimum windage tree Mass air quantity of all branches By Yu Shuxian's quality and air volume Representing, establishing a wind pressure balance equation set, calculating the chord mass and wind quantity of the residual tree Calculating the corrected value of the chord quality air quantity of each residual tree, and the quality air quantity of all branches Iterative correction is carried out, and when convergence condition is reached, high-precision quality air quantity of all branches is obtained And then calculate the high-precision volume air quantity of all branches . According to the invention, the mass flow is used as a calculation reference, and the calculation of the mass flow is not interfered by the change of the air density, so that the airflow dynamics in the mine ventilation network can be more accurately simulated, the calculation error of the volume flow is obviously reduced, and the calculation precision and efficiency are improved.

Inventors

DENG JUN
YAN ZHENGUO
ZHANG LONGCHENG
WANG ZHENPING
WANG YANPING
CHANG LIANG
CAO YUQI
QIN ZHIXIN
YANG BO

Assignees

西安科技大学
江苏智矿科技有限公司

Dates

Publication Date: 20260508
Application Date: 20260302

Claims (9)

1. The high-precision network resolving method based on mass-volume flow double coupling is characterized by comprising the following steps of: S1, inputting a ventilation network topology structure diagram of a target mine, wherein the ventilation network topology structure diagram comprises a plurality of branches Input each branch Calculates each node Relative humidity at Density of air Calculating the volume and the air quantity of the air flow in each branch Mass air quantity M, branch ventilation resistance H and friction air resistance of branch inner air flow ; S2, generating a minimum wind resistance tree by using a Kruskal algorithm ; S3, the mass air quantity of all branches is calculated The air quantity is obtained by the chord quality of the tree Representing, establishing a wind pressure balance equation set, and solving the mass air quantity of the chord of the residual tree Comprising the following steps: s301, in a ventilation network topology structure diagram, the wind resistance tree belongs to the minimum wind resistance tree Branches in (a) Called branch, the topology structure diagram of the ventilation network does not belong to the minimum windage tree Branches of (2) The method is called the chord of the residual tree, each chord of the residual tree corresponds to an independent loop of the ventilation network, the chord air quantity of the residual tree is extracted, the chord air quantity of the residual tree is the air quantity of the corresponding independent loop, and the chord air quantity of the residual tree comprises the chord mass air quantity of the residual tree Volume and air quantity of chord of the tree Calculating the chord quality and air quantity of the residual tree Volume and air quantity of chord of the tree : Mass air quantity with string of remaining tree = Wind current density with remaining tree strings Volume and air quantity of the chord of the residual tree = , Then there is In the above-mentioned method, the step of, , Represent the first The mass air quantity of each remaining tree chord, Represent the first The volume and the air quantity of the chord of each residual tree, Represent the first Wind current density of the remaining tree strings; S302, adopting mass air quantity of the remaining tree strings Establishing a wind pressure balance equation, which comprises the following steps: The wind pressure balance equation for the independent loop is: In the above-mentioned method, the step of, Represent the first Branch and the first Topology of individual loops; Represent the first Wind pressure of each branch; ; Mass air quantity of all branches All written as the chord quality and air quantity of the residual tree Establishing a wind pressure balance equation set, and solving to obtain the mass air quantity of the chord of the residual tree : Represented by Mass air quantity of each remaining tree string A wind pressure balance equation set is formed, Represent the first The wind pressure balance equation is solved to obtain Mass air quantity of each remaining tree string , Representing the first of all branches The friction wind resistance of the branches is high, Representing the first of all branches The wind current density of the individual branches is such that, Representing the first of all branches The wind pressure of the fan contained in each branch, Representing the first of all branches The natural wind pressure contained in each branch is that, Representing the first of all branches The mass air quantity of each branch; s4, calculating a wind pressure error of a wind pressure balance equation, performing Taylor expansion on the wind pressure balance equation, adopting a method of simplifying a Jacobi matrix by adopting a Seidell principle, and obtaining a correction value of the chord quality and the wind quantity of each residual tree based on the wind pressure error; S5, based on correction value of chord quality air quantity of each residual tree, quality air quantity of all branches Performing iterative correction, and obtaining high-precision quality air quantity of all branches when convergence conditions are reached after t iterations ; S6, adopting high-precision quality air quantity of all branches Calculating to obtain high-precision volume air quantity of all branches 。
2. The high-precision network resolving method based on mass-volume flow double coupling according to claim 1, wherein the minimum windage tree is selected in S2 by using Kruskal algorithm The method of (1) comprises the following steps: s201, initializing stage, namely, for all branches According to friction wind resistance Performing ascending sort to form ordered branch set = , Represent the first The number of branches is chosen such that, Represent the first Branches of the Creating an empty set Initializing a full node set Full node set Storing all nodes Initializing and checking the structure to make each node Independent constituent element connected components and checking And The array satisfies And is also provided with , Representing a full node set In presence of nodes , Representing a full node set With a unique node therein ; S202, in the core processing stage, algorithm traversal Each branch of (3) For the current branch Acquiring the initial node through the union operation And end node , Representing current branches Is the starting node of (1) The root of the connected component where it is located, , Representing current branches End node of (a) Root of connected component , Representing a function of finding a node when it is detected In the time-course of which the first and second contact surfaces, Belonging to different connected components, corresponding branches Adding in Aggregation, corresponding initial node Through union of collections Operation merging Updating and searching the parent and rank arrays of the set when the corresponding connected components are detected When skipping the branch Preventing loop formation, adding new branches each time successfully Post-detection termination condition , Representing a set of nodes The number of nodes in the (c) tree, Representation of The number of branches in the set; S203, a result output stage, namely recording the moment when the algorithm meets the termination condition Is the minimum windage tree And output.
3. The high-precision network calculation method based on mass-volume flow double coupling according to claim 1, wherein the method for calculating the wind pressure error of the wind pressure balance equation, performing taylor expansion on the wind pressure balance equation and simplifying the Jacobi matrix by adopting the seidel principle comprises the following specific steps: wind pressure balance equation set After the iteration, the mass air quantity of the remaining tree strings is obtained : In the above-mentioned method, the step of, Represent the first The first iteration The mass air quantity of the strip remaining tree strings; Will be the first The wind pressure balance equation in iteration is developed by a Taylor formula: In the above-mentioned method, the step of, Represent the first The wind pressure balance equation is set up in the form of a wind pressure balance equation, Represent the first +1 Iteration time The mass air quantity of the string of the tree, Represent the first The first iteration The correction value of the chord quality and the air quantity of the strip remaining tree, , Representing the number of nodes to be connected, The number of branches is indicated and, ; The second-order trace of the inner secondary higher order is ignored The sub-linear approximation is: Writing the above as a matrix form: The simplified calculation is as follows: In the above-mentioned method, the step of, Represent the first The wind pressure balance equation pair 1 The chord quality and the air quantity of the residual tree The first derivative Jacobi coefficient matrix has strict principal diagonal advantage, and the same-order trace is omitted, namely, each item on a non-diagonal line in the coefficient matrix is omitted: The above can be written as: Namely: In the above-mentioned method, the step of, Represent the first Wind pressure error of each wind pressure balance equation.
4. The high-precision network resolving method based on mass-volume flow double coupling according to claim 3, wherein the method for obtaining the correction value of the mass air quantity of each remaining tree chord based on the air pressure error in the step S4 is as follows: Will be Described as Is derived from the function of (2) Is defined by the following specific expression: In the above-mentioned method, the step of, Represent the first The mass flow rate of the individual branches, Represent the first Branches and loops Is used for the topology relation of (a), Represent the first The mass air quantity of the string of the tree, Representing the number of nodes to be connected, The number of branches is indicated and, , ; Combining the wind pressure balance equation to obtain the following components: In the above-mentioned method, the step of, Represent the first The wind pressure balance equation is set up in the form of a wind pressure balance equation, Represent the first Branch and the first The topology of the individual loops is such that, Represent the first The friction wind resistance of the branches is high, Represent the first Branches and loops Is used for the topology relation of (a), Represent the first The mass air quantity of the string of the tree, Represent the first The wind current density of the individual branches is such that, Represent the first The wind pressure of the fan contained in each branch, Represent the first The natural wind pressure contained in each branch is that, Representing branches Is a friction windage of (2); will be of the above formula For a pair of And (5) conducting derivation to obtain: In the above-mentioned method, the step of, Represent the first Branches of the Representing a constant, therefore When (when) For a pair of Deriving only when When the derivative is not 0, thus In the same way as the above, the following steps, The method can obtain: Expressed as: In the above-mentioned method, the step of, Represent the first The first iteration The mass air quantity of each remaining tree chord, Represent the first Branches and loops Is a topological relation of (1); Mass air quantity correction value of remaining tree chord The calculation formula of (2) is as follows: 。
5. the high-precision network resolving method based on mass-volume flow double coupling as claimed in claim 1, wherein the mass air quantity for all branches in S5 The iterative correction method comprises the following steps: First, the Second time The mass air quantity of the branches is : In the above-mentioned method, the step of, Representing the number of nodes to be connected, The number of branches is indicated and, , Represent the first The first iteration The mass air quantity of the branch strip, Represent the first Branches and loops Is used for the topology relation of (a), Represent the first Post-iteration first And the chord quality and air quantity correction value of the strip remaining tree.
6. The high-precision network solution method based on mass-volume flow dual coupling according to claim 5, wherein the convergence condition in S5 is: continuously iterating the quality air quantity correction value of each independent loop After the t-th iteration, the maximum independent loop quality air quantity correction value Less than the set precision Namely there is Solving to obtain high-precision quality air quantity of all branches = , 。
7. The method for high-precision network resolution based on dual mass-volume flow coupling as claimed in claim 6, wherein, Is 0.00001-0.01 kg/min.
8. The high-precision network resolving method based on mass-volume flow double coupling as claimed in claim 7, wherein in S6, high-precision mass air quantity of all branches is adopted Calculating high-precision volume air quantity of all branches The method of (1) is as follows: In the above-mentioned method, the step of, , Indicating the t-th iteration after reaching the convergence condition The high-precision volume air quantity of each branch, Represent the first After the iteration reaches the convergence condition The mass air quantity of the branch strip, Represent the first Wind flow density of the individual branches.
9. The high-precision network resolution method based on mass-volume flow dual coupling of claim 1, wherein in S1 at each node The position is provided with measuring points to obtain each branch Is a data of (a) a data of (b).

Description

High-precision network resolving method based on mass-volume flow double coupling Technical Field The invention relates to the technical field of mine ventilation, in particular to a high-precision network resolving method based on mass-volume flow double coupling. Background Coal resources are basic stones for energy safety in China, and along with the increasing frequency and normalization of deep mine exploitation activities, the improvement of the intelligent mine ventilation technology level has become an unprecedented task. The intelligent transformation of the ventilation system is not only a core link of constructing an intelligent mine, but also a key point of guaranteeing the safe production of the mine, and the scientificity and rationality of the design directly influence the safety and the production efficiency of the mine. By introducing an advanced intelligent technology, the dynamic monitoring and the optimized management of the ventilation system can be realized, so that the efficient and sustainable development of mine production is promoted on the premise of ensuring the safety, and the powerful guarantee is provided for national energy safety. Zhong Deyun and the like are used for searching the minimum independent closed loop by adopting a depth-first search (DFS) method aiming at the problems of disordered fan management, circulating wind phenomenon and the like in the traditional mine ventilation design, and a multi-fan multi-stage machine station ventilation network resolving method is established; lei Yajun and the like aiming at the problems of inaccurate measurement and general change of mine ventilation parameters, adopting a multi-situation self-adaptive ventilation network resolving method, realizing high reduction degree of the operation state of a ventilation network by improving a node pressure equation and a clustering noise reduction algorithm, qin Tao and the like aiming at the problems of incapability of intelligent adjustment, regional pollution, unsmooth ventilation and the like of a traditional mine ventilation system, dynamically adjusting the parameters of the ventilation system based on a three-dimensional network model and a mathematical model, combining sensor data, improving the intelligent adjustment and optimization capability of air flow in a mine, tan Xiaolin and the like aiming at the problems of low convergence speed, incapability of convergence and the like in mine ventilation network resolving, adopting a method of firstly detecting whether a self-loop exists in the ventilation network resolving, simplifying the ventilation network and adopting a loop air quantity resolving method if the self-loop exists, effectively improving the calculation efficiency and the accuracy by adopting a node air pressure resolving method if the self-loop exists, xie Biao and the like aiming at the problems of hysteresis and high hysteresis of the traditional mine ventilation network resolving method, adopting a tree-shaped simplified processing and graph theory method, effectively improving the resolving speed and accuracy, realizing the data processing method of the effect of the mine intelligent ventilation system for the mine system, and the traditional ventilation network resolving method of Liu Hui and the like, and the problem of high-dimensional self-adaption ventilation network can not be realized by combining the traditional three-dimensional simulation and the three-dimensional self-adaption ventilation network resolving mode and the problem, the calculation efficiency and the accuracy are improved. Han Xiaohui and the like are used for solving the problems of low resolving speed and poor resolving precision of a mine ventilation network, and the method such as a mesh method, a weight assignment method, a double-channel method and the like is adopted to realize rapid convergence and improve resolving speed and resolving precision. In the existing mine ventilation network air volume resolving technology, the conventional resolving method based on volume flow is widely adopted to ignore the influence of air density change, so that obvious defects exist, the conventional resolving method assumes that the air density is kept constant, but the environmental conditions of all areas are obviously changed along with the increase of the mine depth. In the air volume iterative calculation process, the influence of air density change on ventilation flow is not fully considered by the traditional method, so that larger deviation exists between calculated volume flow and actual flow in mine environments with uneven air density distribution. Such deviations not only prevent high-precision air distribution, but may also adversely affect the design and operation of the ventilation system, thereby affecting the safe production of the mine. Disclosure of Invention The invention aims to solve the problems, and provides a high-precision network resolving method based on mass-volume fl