A complex ventilation network solving method and a solving device

Abstract

The application discloses a kind of complex ventilation network solving method and solving device, it is related to mine ventilation technical field;In the method of the application, node wind pressure equation set and its auxiliary equation, fan wind pressure characteristic function etc. are unified as a generalized matrix equation, the composition of the equation, independent variable and dependent variable vector, jacobian matrix etc. are defined, so as to realize direct solution using Newton method, and the amount of programming is small;Since the number of equation and the number of independent variable are not limited by the method, as long as the number of the two is related, therefore, by adding fan wind pressure characteristic function, direct solution of wind network can be realized when fan operating point is not fixed;By adding auxiliary equation reflecting the parameter relationship of multiple inlet and outlet of wind network, direct solution of multiple inlet and outlet wind network can be realized.

Description

Technical Field

[0001] This invention relates to the field of mine ventilation technology, and in particular to a method and apparatus for solving complex ventilation networks. Background Technology

[0002] With the increasing emphasis placed on mine safety and green mining by the state, mine ventilation solutions are becoming increasingly refined and real-time. With the development of algorithms and the improvement of computing power, not only can the calculation of air volume in traditional mine shafts be done in real time, but the simulation of air leakage in areas such as goafs and waste rock piles can also be solved using the concept of equivalent ventilation network calculation. Unlike traditional ventilation networks, the equivalent ventilation networks abstracted from goafs and waste rock piles have the characteristics of multiple sources and sinks. Traditional ventilation network calculation methods based on loop analysis are no longer suitable for solving complex ventilation networks in the following situations due to the complexity of virtual loop assumptions, the lack of intuitiveness in obtaining independent loop matrices, and the difficulty in external coupling: 1. Non-closed ventilation networks; 2. Ventilation networks with fans whose operating points are not fixed (local ventilators and main ventilators); 3. Ventilation networks with multiple inlets and outlets. Summary of the Invention

[0003] In view of the above problems, the purpose of this invention is to provide a method and apparatus for solving complex ventilation networks.

[0004] The first aspect of this application provides a method for solving complex ventilation networks, the method comprising:

[0005] Establish an equivalent ventilation network diagram of the complex wind network to be solved, establish a basic correlation matrix based on the equivalent ventilation network diagram, determine the branch wind resistance column vector based on the basic correlation matrix and the wind resistance values ​​of each branch, and establish a solution function for the nodal wind pressure equation system containing unknowns P and Q based on the basic correlation matrix and the branch wind resistance column vector.

[0006] Establish a correction vector for the solution function of the nodal wind pressure equations, and establish an auxiliary solution function based on the introduced unknowns other than P and Q.

[0007] The solution function for the nodal wind pressure equations is modified and combined with the auxiliary solution function to form the solution function for the equations to be solved.

[0008] The independent and dependent variables of the solution function of the system of equations to be solved are preprocessed. Specifically, the Jacobian matrix of the solution function of the system of equations is assembled based on the derivative of the dependent variable with respect to the independent variable.

[0009] Given the initial values ​​of the branch wind volume column vector Q, the nodal wind pressure column vector P, and other auxiliary unknowns, the wind network is solved using Newton's method for the solution function of the problem under study. Under the given termination conditions, the approximate solution of the equation is obtained as the approximate value of the actual wind volume of the wind network.

[0010] The beneficial effects of this invention are as follows: This application provides a method for solving complex ventilation networks, unifying the nodal wind pressure equations and their auxiliary equations, fan wind pressure characteristic functions, etc., into a generalized matrix equation. The composition, independent and dependent variable vectors, Jacobian matrix, etc., of this equation are defined, thus enabling direct solution using Newton's method with minimal programming complexity. Since this method does not limit the number of equations or independent variables, as long as they are related, by adding the fan wind pressure characteristic function and solving using this method, direct solution of the ventilation network when the fan operating point is not fixed can be achieved. Adding auxiliary equations reflecting the parameter relationships of multiple inlets and outlets of the ventilation network enables direct solution of multi-inlet and outlet ventilation networks. Furthermore, taking a non-closed ventilation network, represented by a local mine, as an example, this invention, based on the actual situation, makes isolated branches or isolated nodes equivalent at the boundary of the actual ventilation network. The air volume balance equation is corrected by adding an external air volume compensation vector, and the branch resistance equation is corrected by adding an external wind pressure compensation vector. This yields a corrected nodal wind pressure equation set to describe the non-closed ventilation network. Attached Figure Description

[0011] Figure 1 This is a flowchart of a complex ventilation network calculation method provided in an embodiment of this application.

[0012] Figure 2 This is a schematic diagram of a local wind network that conforms to scenario one in the embodiments of this application.

[0013] Figure 3 yes Figure 2 Equivalent ventilation network diagram of a local ventilation network.

[0014] Figure 4 This is a schematic diagram of a local wind network that conforms to scenario two in the embodiments of this application.

[0015] Figure 5 yes Figure 4 Equivalent ventilation network diagram of a local ventilation network.

[0016] Figure 6 This is a schematic diagram of a local wind network that conforms to scenario three in the embodiments of this application.

[0017] Figure 7 yes Figure 6 Equivalent ventilation network diagram of a local ventilation network.

[0018] Figure 8 This is a schematic diagram of a local wind network that conforms to scenario four in the embodiments of this application.

[0019] Figure 9 yes Figure 8 Equivalent ventilation network diagram of a local ventilation network.

[0020] Figure 10 This is a schematic diagram of a local wind network that conforms to scenario five in the embodiments of this application.

[0021] Figure 11 yes Figure 10 Equivalent ventilation network diagram of a local ventilation network. Detailed Implementation

[0022] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0023] The technical problem to be solved by this application and related concepts will be explained below:

[0024] With the increasing emphasis placed on mine safety and green mining by the state, mine ventilation solutions are becoming increasingly refined and real-time. With the development of algorithms and the improvement of computing power, not only can the calculation of air volume in traditional mine shafts be done in real time, but the simulation of air leakage in areas such as goafs and waste rock piles can also be solved using the concept of equivalent ventilation network calculation. Unlike traditional ventilation networks, the equivalent ventilation networks abstracted from goafs and waste rock piles have the characteristics of multiple sources and sinks. Traditional ventilation network calculation methods based on loop analysis are no longer suitable for the complex ventilation network calculation described in this application due to the complexity of virtual loop assumptions, the lack of intuitiveness in obtaining independent loop matrices, and the difficulty in external coupling. The following section will specifically introduce the specific types and characteristics of the complex ventilation networks described in this application.

[0025] In existing technologies, wind network calculations are based on the following three fundamental laws:

[0026] (1) Law of Conservation of Mass: The total airflow into a node must be equal to the total airflow out of that node. The matrix equation of this law is as follows:

[0027] (a)

[0028] In the formula, B is the basic correlation matrix of the wind network, and Q is the branch wind volume column vector.

[0029] (2) Law of Conservation of Energy: The total voltage drop of any closed loop must be equal to zero. The matrix equation of the law is as follows:

[0030] C( (b)

[0031] In the formula, P is the nodal wind pressure column vector, and R is the branch wind resistance column vector. Let C be the column vector of fan pressure, and C be the independent loop matrix. This represents the fan's airflow. The subscript `diag` indicates that the vector has been diagonalized into a matrix, and the symbol `||` indicates taking the absolute value.

[0032] (3) Branch Resistance Law: The pressure difference between the two ends of a branch is equal to the algebraic sum of the resistance and the driving force on the branch. The matrix equation of the law is as follows:

[0033] (c)

[0034] In the formula, the superscript T indicates the transpose of the matrix.

[0035] The fundamental incidence matrix B and the independent loop matrix C have the following relationship:

[0036] (d)

[0037] Substituting equation d into equation c, it can be proven that equations b and c are equivalent. Therefore, solving equations a and b simultaneously, or equations a and c simultaneously, can be used to solve the ventilation network.

[0038] By combining equations a and d, we can obtain:

[0039] (e)

[0040] In the formula This refers to the amount of airflow passing through the remaining branches of the windbreak.

[0041] Substituting equation e into equation b, we obtain the wind network solution equation as shown in equation f. This method is called the loop analysis method, which is equivalent to solving equations a and b simultaneously.

[0042] (f)

[0043] The above method has the following drawbacks:

[0044] 1. The unknowns in the solution equation of this method are: , These two unknowns are essentially wind volume, but they cannot be unified. Therefore, this equation is difficult to solve directly using Newton's method or similar methods. It is often solved indirectly using a three-step method, the idea of ​​which is as follows:

[0045] 1.1 First, assume that the fan is at a fixed air pressure operating point, that is, take the fan ventilation pressure. =c, calculate the total inflow (or outflow) air volume in the mine. and Further based on the mine ventilation resistance equation ,make Calculate the mine air resistance .

[0046] 1.2, Order , Simultaneous equations for mine ventilation resistance And fan ventilation pressure function Find the value at this time. . ) is a wind turbine pressure characteristic curve function fitted using an algorithm.

[0047] 1.3. Based on the calculated operating point of the fan According to the proportional relationship and Solve for the column vector of air volume of the wind network branches. Further, the branch airflow column vector is obtained according to equation e. .

[0048] 2. This method requires finding the cotree first. The currently recognized best method for finding the cotree is the depth (or breadth) first search algorithm, but the programming implementation of this algorithm is relatively complex.

[0049] 3. This method requires the establishment of an independent loop matrix C. Whether it is obtained by various methods based on the cotree or directly by transforming the basic incidence matrix, the implementation is quite complicated.

[0050] 4. Since real-world wind networks are not completely closed, it is necessary to assume virtual loops to connect the various entrances and exits between the wind network and the outside world, thereby achieving the theoretical closure of the wind network. For equivalent wind networks represented by goaf areas, the number of airflow entrances and exits is very large, making artificial assumptions extremely complicated.

[0051] 5. This method is difficult to use for local ventilation fan problems.

[0052] Furthermore, by combining equations a and c, we obtain a system of equations in the form of equation h. The method of solving the wind network using this system of equations is called the nodal pressure method.

[0053] (h1)

[0054] (h2)

[0055] In this system of equations, Q, P, The unknowns have been verified, and the system of equations has a solution. This model has a long history, with equation h1 called the branch resistance equation and equation h2 called the nodal airflow balance equation, together known as the nodal air pressure equation system.

[0056] Compared with loop analysis, this method has the following advantages: 1. It does not require finding the remainder tree or the independent loop matrix, saving preparation time for equation calculation parameters. 2. This system of equations does not require the wind network to be closed.

[0057] However, this method also has the following drawbacks:

[0058] 1. This system of equations has two matrix equations and three unknowns: P, Q, ... Therefore, similar to the loop analysis method, its disadvantage is that it cannot be directly solved using Newton's method. The conventional solution method is the elimination method: all sub-equations are listed, and the nodal wind pressure or branch air volume is eliminated by substitution, resulting in a nonlinear equation system for solution. This method cannot fully utilize the convenience of matrix operations, and it requires a large amount of programming and is prone to errors.

[0059] 2. The number of unknowns is greater in the loop analysis method than in the loop analysis method, resulting in greater computational resource consumption during the solution process. Currently, after comparative analysis, the conclusion is that if the parameter preparation time of the loop analysis method is also taken into account, the difference in solution time between the two methods is minimal, and the loop analysis method is superior.

[0060] 3. Although this system of equations does not require the wind network to be closed, from the essence of each sub-equation, there cannot be isolated nodes (i.e., nodes connected to only a single branch) or isolated branches (i.e., branches connected to only a single node).

[0061] Based on the foregoing introduction of related technologies, the purpose of this application is to overcome the shortcomings of the nodal wind pressure method, and thereby utilize the advantages of this method to solve complex ventilation networks with the following characteristics:

[0062] 1. The ventilation network is not closed.

[0063] 2. There are fans (local ventilators and main ventilators) in the ventilation network whose operating points are not fixed.

[0064] 3. A ventilation network with multiple entrances and exits.

[0065] In view of this, this application provides a method for solving complex ventilation networks. It unifies the nodal wind pressure equations and their auxiliary equations, fan wind pressure characteristic functions, etc., into a generalized matrix equation. The composition, independent and dependent variable vectors, and Jacobian matrix of this equation are defined, thus enabling direct solution using Newton's method with minimal programming complexity. Since this method does not limit the number of equations or independent variables, as long as the total number is equal, by adding the fan wind pressure characteristic function and solving using this method, direct solution of the ventilation network when the fan operating point is not fixed can be achieved. Adding auxiliary equations reflecting the parameter relationships of multiple inlets and outlets of the ventilation network enables direct solution of multi-inlet and outlet ventilation networks. Furthermore, taking a non-closed ventilation network, represented by a local mine, as an example, this invention, based on the actual situation, creates equivalent isolated branches or isolated nodes at the boundary of the actual ventilation network. The air volume balance equation is corrected by adding an external air volume compensation vector, and the branch resistance equation is corrected by adding an external wind pressure compensation vector. This yields a corrected nodal wind pressure equation set to describe the non-closed ventilation network.

[0066] It should be noted that the execution subject of the solution method provided in this application embodiment can be a solution device. This solution device can be implemented as part or all of a computer device through software, hardware, or a combination of both. The computer device can be a server or a terminal. The server can refer to a single server or a server cluster composed of multiple servers. The terminal can be a desktop computer, tablet computer, laptop computer, or smartphone, etc. This application embodiment does not limit the specific type of terminal. In the following method embodiments, the execution subject is always described using a computer device as an example.

[0067] Please refer to Figure 1 It shows a flowchart of a complex ventilation network calculation method provided in an embodiment of this application, such as... Figure 1 As shown, the solution method for this complex ventilation network includes the following steps:

[0068] Step S1: Establish an equivalent ventilation network diagram of the complex wind network to be solved. Establish a basic correlation matrix based on the equivalent ventilation network diagram. Determine the branch wind resistance column vector based on the basic correlation matrix and the wind resistance values ​​of each branch. Establish a solution function for the nodal wind pressure equation system containing unknowns P and Q based on the basic correlation matrix and the branch wind resistance column vector.

[0069] Step S101: Establish an equivalent ventilation network diagram of the complex wind network to be solved. Specifically, based on the inlet and outlet conditions connected to the outside world in the complex wind network to be solved, and the wind pressure characteristics of the fans in the wind network, perform equivalent processing on the complex wind network to be solved. Then, assume the airflow direction and make all equivalent branches into directed branches. Finally, connect the branches and nodes to construct the equivalent ventilation network diagram.

[0070] The method in this application is mainly aimed at the following five types of complex ventilation networks. The method for establishing the equivalent ventilation network diagram of the complex ventilation network to be solved in the five types of scenarios is described in detail below.

[0071] Scenario 1: The inflow or outflow of air to the outside of the complex ventilation network to be calculated is a fixed value. The method for constructing an equivalent ventilation network diagram is as follows. For ease of understanding, examples will be provided below.

[0072] like Figure 2 The diagram shown is a schematic of a local wind network that conforms to scenario one; Figure 3 What is shown is Figure 2 Equivalent ventilation network diagram of a local ventilation network. The equivalent treatment process is described in detail below.

[0073] First, such as Figure 2As shown, hollow arrows represent the direction of the working face advance, and solid arrows represent the direction of airflow. Point a is the inlet of this local air network, and point i is the outlet of this local air network. The air volume at these two points is a fixed value, both being [missing information]. .

[0074] Figure 2 The method for constructing an equivalent ventilation network diagram of the local airflow network shown is as follows:

[0075] A single wind path (such as bh) or multiple wind paths in series (such as bc-ce-de) in the wind network is equivalent to a branch, and the intersection of two or more wind paths (such as point e) is equivalent to a node.

[0076] Specifically, the entrances and exits connecting the wind network to the outside world are treated as isolated nodes (i.e., nodes connected to only a single branch), and used as the objects for subsequent airflow compensation. In this example, Figure 2 Points a and i in the equation are equivalent to: Figure 3 Nodes in and The specific equivalence relationships are shown in Table 1 below.

[0077] Table 1

[0078] Wind Road Branches Intersection or location node ab b bcde h bh e eh g efg a hg i gi

[0079] Next, assuming the airflow direction is as shown in the diagram, all branches are made directed branches; finally, the branches and nodes are connected to construct an equivalent ventilation network diagram as shown. Figure 3 As shown.

[0080] Scenario 2: The pressure difference between the inlets and outlets of the complex ventilation network to be solved and the external environment is constant. The method for constructing the equivalent ventilation network diagram is as follows. For ease of understanding, examples will be provided below.

[0081] like Figure 4 The diagram shown is a schematic of a local wind network that conforms to scenario two. Figure 5 What is shown is Figure 4 Equivalent ventilation network diagram of a local ventilation network. The equivalent treatment process is described in detail below.

[0082] like Figure 4 As shown, hollow arrows represent the direction of the working face advance, and solid arrows represent the direction of airflow. Point a is the inlet of this local air network, and point i is the outlet of this local air network. The pressure difference between these two points is a fixed value. The air pressure at point a is used as a reference point and its value is set to 0; the air pressure at point i is less than 0 and its value is -. Airflow moves from areas of high wind pressure to areas of low wind pressure.

[0083] Figure 4 The method for constructing an equivalent ventilation network diagram of the local airflow network shown is as follows:

[0084] A single wind path (such as bh) or multiple wind paths connected in series (such as bc-ce-de) in the wind network is equivalent to a branch, and the intersection of two or more wind paths (such as point e) is equivalent to a node.

[0085] Specifically, the channels through which the wind network's airflow enters and exits the outside world are treated as isolated branches (i.e., branches connected to only a single node), and these branches are used as the objects of subsequent wind pressure compensation. For example... Figure 4 The wind paths ab and gi in the diagram are equivalent to Figure 5 Branches in and At this point, points a and i have no corresponding nodes in the equivalent wind network and are therefore discarded. The specific equivalent relationships are shown in Table 2 below.

[0086] Table 2

[0087] Wind Road Branches Intersection or location node ab b bcde h bh e eh g efg hg gi

[0088] Next, assuming the airflow direction, all branches are made directed. Finally, the branches and nodes are connected to construct an equivalent ventilation network diagram as follows: Figure 5 As shown.

[0089] Scenario 3: A main ventilation fan is located at the connection point between the complex ventilation network to be solved and the outside world, and the operating point of the fan is yet to be determined. The method for constructing an equivalent ventilation network diagram in this case is as follows. For ease of understanding, examples will be provided below.

[0090] like Figure 6 The diagram shown is a schematic of a local wind network that conforms to scenario three. Figure 7 What is shown is Figure 6 Equivalent ventilation network diagram of a local ventilation network. The equivalent treatment process is described in detail below.

[0091] like Figure 6 As shown, hollow arrows represent the direction of the working face advance, and solid arrows represent the airflow direction. Points a and i are the inlet and outlet of this local ventilation network. The air pressure at point a is taken as a reference point and set to 0. A main ventilation fan is located at point i; the operating point of this fan is undetermined, and the air pressure at this point caused by the fan is less than 0, a value that is unknown and set as -. The airflow moves from a position of high pressure to a position of low pressure. The air volume passing through the fan is an unknown quantity, denoted as [missing information]. .

[0092] Figure 6 The method for constructing an equivalent ventilation network diagram of the local airflow network shown is as follows:

[0093] A single branch (such as bh) and multiple connected wind paths (such as bc-ce-de) in the wind network are equivalent to a single branch, and the intersection of two or more wind paths (such as point e) is equivalent to a node.

[0094] Specifically, the main ventilation fan is represented as a branch with zero air resistance and nodes at both ends. The node directly connected to the outside is represented as an isolated node (i.e., a node connected only to a single branch) and is used as the object of subsequent airflow compensation; the airflow on this branch is... The wind pressure at the isolated node is - This wind pressure is also the ventilation power of the equivalent branch. and The relationship is given by the wind pressure characteristic curve. The other inlet / outlet connecting the wind network to the outside world is also equivalent to an isolated node. Figure 6 Points a and i in the equation are equivalent to: Figure 7 Nodes in and The wind turbine is equivalent to a node. and And the branches connecting the two nodes The specific equivalence relationships are shown in Table 3 below.

[0095] Table 3

[0096] Wind Road Branches Intersection, location or fan Nodes or combinations ab b bcde h bh e eh g efg a hg Fan gi

[0097] Next, assuming the airflow direction, all branches are made directed. Finally, the branches and nodes are connected to construct an equivalent ventilation network diagram as follows: Figure 7 As shown.

[0098] Scenario 4: The complex ventilation network to be solved has local ventilators on its branches. The method for constructing the equivalent ventilation network diagram is as follows. For ease of understanding, an example will be given below.

[0099] like Figure 8 The diagram shown is a schematic of a local wind network that conforms to scenario four. Figure 9 What is shown is Figure 8 Equivalent ventilation network diagram of a local ventilation network. The equivalent treatment process is described in detail below.

[0100] like Figure 8 As shown, hollow arrows represent the direction of the working face advance, and solid arrows represent the airflow direction. Points a and i are the inlet and outlet of this local air network, and the airflow at these two points is a fixed value, both being [missing information]. A local ventilation fan is installed on the air duct fg. The operating point of this fan is undetermined, and the ventilation pressure generated by the fan is an unknown quantity, denoted as - The air volume is an unknown quantity, set as follows: .

[0101] Figure 8 The method for constructing an equivalent ventilation network diagram of the local airflow network shown is as follows:

[0102] A single wind path (such as bh) and multiple wind paths connected in series (such as bc-ce-de) in the wind network are equivalent to a branch, and the intersection of two or more wind paths (such as point e) is equivalent to a node.

[0103] Specifically, the entrances and exits connecting the wind network to the outside world are treated as isolated nodes (i.e., nodes connected only to a single branch) and used as the objects of wind volume compensation; such as Figure 8 Points a and i in the equation are equivalent to: Figure 9 Nodes in and A local ventilation fan can be represented as a branch with zero air resistance and nodes at both ends, such as... Figure 8 The local fan in the middle is equivalent to Figure 9 Nodes in and And the branches connecting the two nodes The airflow on this branch is Ventilation pressure is - , and The relationship is given by the wind pressure characteristic curve. The specific equivalent relationship is shown in Table 4 below.

[0104] Table 4

[0105] Wind Road Branches Intersection, location or fan Nodes or combinations ab b bcde h bh e eh Fan ef-fan g Fan-g i hg a gi

[0106] Next, assuming the airflow direction, all branches are transformed into directed branches. Finally, the branches and nodes are connected to construct an equivalent ventilation network diagram, as shown below. Figure 9 As shown.

[0107] Scenario 5: The complex ventilation network to be solved has more than two entrances / exits to the outside. The method for constructing the equivalent ventilation network diagram is as follows. For ease of understanding, examples will be provided below.

[0108] like Figure 10 The diagram shows a local wind network that meets scenario five. Figure 11 What is shown is Figure 10 Equivalent ventilation network diagram of a local ventilation network. The equivalent treatment process is described in detail below.

[0109] like Figure 10 As shown, hollow arrows represent the direction of face advancement, and solid arrows represent the airflow direction. Points a and i are the inlet and outlet of this local ventilation network. Additionally, calculations indicate that airflow leaks into point e in the direction of the goaf; assuming the leakage channel is airflow path ef. Therefore, the total airflow entering the ventilation network from points a and j is equal to the airflow exiting from point i. This total airflow is a fixed value, its magnitude being... Assume the leakage air volume from point j is... Then the air volume of a entering the ventilation network is .

[0110] Figure 10 The method for constructing an equivalent ventilation network diagram of the local airflow network shown is as follows:

[0111] A single wind path (such as bh) and multiple wind paths connected in series (such as bc-ce-de) in the wind network are equivalent to a branch, and the intersection of two or more wind paths (such as point e) is equivalent to a node.

[0112] Specifically, the entrances and exits connecting the wind network to the outside world are treated as isolated nodes (i.e., nodes that are only connected to a single branch) and are used as the objects of subsequent wind volume compensation. Figure 10 Points a, i, and j in the equation are equivalent to: Figure 11 Nodes in , and The specific equivalence relationships are shown in Table 5 below.

[0113] Table 5

[0114] Wind Road Branches Intersection or location node ab b bcde h bh e eh g efg a hg i gi j je

[0115] Next, assuming the airflow direction, all branches are made directed. Finally, the branches and nodes are connected to construct an equivalent ventilation network diagram as follows: Figure 11 As shown.

[0116] Step S102: Establish the basic correlation matrix B based on the equivalent ventilation network diagram. This includes the following steps:

[0117] (1) Based on the number of nodes M and branches N in the equivalent ventilation network diagram, construct an M-row N-column all-zero matrix. A Record all nodes in the node number set. Node numbering in the equivalent ventilation network diagram v and A line number r and the branch numbers in the ventilation network diagram. e and A column number c Establish the corresponding relationship. If the number is... The node is located at number At the tail of the branch, according to the correspondence, the matrix is... A The corresponding OK The number in the column position is rewritten as 1. If the number is The node is located at number The head of the branch, according to the correspondence, will be the matrix. A The corresponding OK The number in the column position is rewritten as -1.

[0118] (2) Obtain the rank of the matrix A , if = M, assign to A . If < M, according to the difference between and M, delete the corresponding number of rows from , and delete the node numbers corresponding to these rows from the node number set. Let the obtained new matrix be A , and its number of rows is , , the rank of is , if = = = , assign the new matrix to B . At this time, B is the basic incidence matrix, and its number of rows is denoted as , and the number of columns remains N.

[0119] (3) According to the basic incidence matrix B , establish the corresponding relationship between the node number in the node number set and the row number v of B . Inherit the corresponding relationship between the branch number r and the column number A of e as the corresponding relationship of the matrix c B .

[0120] R。 Step S103, according to the basic incidence matrix and the wind resistance values of each branch, determine the branch wind resistance column vector B Specifically, according to the number of columns N of the basic incidence matrix e , establish a column vector with N rows and 1 column. According to the corresponding relationship between the branch number c and the column number c , fill in the wind resistance value of branch e at the R -th row of this column vector to obtain the branch wind resistance column vector .

[0121] B Step S104, according to the basic incidence matrix R , the branch wind resistance column vector P , establish a node wind pressure equation group solving function containing unknowns Q and :

[0122] (1a)

[0123] (1b)

[0124] in, Q Let N be the number of rows and 1 column of the branch air volume column vector. P for The nodal wind pressure column vector in row 1 and column 1. (Index) diag This indicates that a vector is diagonalized into a matrix, and the superscript is used. T The symbol || represents the transpose of the matrix, and || denotes taking the absolute value. Equation 1a is called the solution function for the branch resistance equation, and Equation 1b is called the solution function for the nodal airflow balance equation.

[0125] Step S2: Establish the correction vector of the solution function for the nodal wind pressure equation system, and establish the auxiliary solution function based on the introduced unknowns other than P and Q.

[0126] Step S2, more specifically, involves establishing a correction vector for the solution function of the nodal wind pressure equations. Specifically, this means the wind volume compensation column vector for the five scenarios listed in step S1, including scenarios one, three, four, and five. Wind pressure compensation column vector in scenario two The main ventilation fan ventilation pressure compensation column vector in scenario three The column vector for local ventilation fan pressure compensation in scenario four. Based on the introduced unknowns other than P and Q, auxiliary solution functions are established. Specifically, auxiliary solution functions with the same name are established in cases three and five. The following section will further explain step S2 in detail with reference to a specific example from step S1.

[0127] Scenario 1, Figure 3 In the equivalent ventilation network diagram shown, the airflow into (or out of) the outside is a fixed value. Based on this parameter setting, airflow compensation is performed on the corresponding isolated nodes in the ventilation network diagram, and an airflow compensation column vector is established.

[0128] In this case, the methods for establishing the air volume compensation column vector include: according to the basic correlation matrix. B Node number v with line number r The corresponding relationship, look up the number as The row number corresponding to the isolated node ,Establish A column vector with 1 row and 1 column all zeros; if the isolated node is a wind network inlet, then the first zero in this column vector... Rewrite the line as If the isolated node is the outlet of the wind network, then the first node in the column vector... Rewrite the line as Establish an air volume compensation column vector .

[0129] if Figure 3 The basic correlation matrix of the wind network is obtained by deleting nodes. The corresponding row obtained only requires... Perform airflow compensation; since this node is the inlet of the air network, its airflow compensation value should be [value missing]. Therefore, the airflow compensation column vector of this ventilation network is:

[0130]

[0131] Scenario 2, Figure 5 In the ventilation network diagram shown, the pressure difference between the inlets and outlets connecting the ventilation network to the outside is constant, with a value of [value missing]. Assuming the inlet air pressure is 0 and the outlet air pressure is... Based on this parameter setting, wind pressure compensation is performed on the corresponding isolated branches in the ventilation network diagram, and a wind pressure compensation column vector is established.

[0132] In this case, the method for establishing the wind pressure compensation column vector includes: according to the correspondence between branch number e and column number c in the basic correlation matrix B, searching for the number... The column number corresponding to the isolated branch Create an N x 1 column vector consisting entirely of zeros. The inlet air pressure remains constant, while the outlet air pressure is the opposite. Then, the Nth column vector is... The row is rewritten as the wind pressure value at the entrance / exit where the branch connects to the outside world, and a wind pressure compensation column vector is established. .

[0133] Figure 5 Middle branch As the inlet branch of the wind network, its wind pressure compensation value is 0. As the outlet branch of the wind network, its wind pressure compensation value is Therefore, the wind pressure compensation column vector of this ventilation network is:

[0134]

[0135] Scenario 3, Figure 7 The ventilation network shown has a main ventilation fan installed at its connection point with the outside environment, and the fan's operating point is yet to be determined. The air volume passing through the fan is unknown. Corresponding to this air volume, the ventilation pressure on the fan pressure characteristic curve is: Based on this parameter setting, establish the air volume compensation column vector, the ventilation pressure column vector, and the auxiliary solution function.

[0136] In this case, the methods for establishing the air volume compensation column vector include: according to the basic correlation matrix. B Node number v with line number r The corresponding relationship, look up the number as The row number corresponding to the isolated node ,Establish A column vector with 1 row and 1 column all zeros; if the isolated node is a wind network inlet, then the first zero in this column vector... Rewrite the line as If the isolated node is the outlet of the wind network, then the first node in the column vector... Rewrite the line as Establish an air volume compensation column vector .

[0137] further, Figure 7 The isolated nodes in the data are As the wind network entry point, an isolated node For the wind network outlet, the air volume compensation values ​​at these two locations are respectively and If the basic association matrix is ​​obtained by deleting nodes... If the value is obtained from the row in question, then the wind volume compensation vector for that wind network is:

[0138] =

[0139] In this case, the method for establishing the ventilation pressure column vector includes: fitting the wind pressure characteristic curve of the main ventilation fan to a vector based on the ventilation pressure of the main ventilation fan. and The expression generally represents a polynomial of degree no more than four. Based on the correspondence between branch number n and column number c, find the polynomial with the following number: The column number corresponding to the equivalent branch of the main fan Create an N x 1 column vector containing only zeros. Then, define the Nth column vector containing zeros. Rewrite the line as Establish the main ventilation fan ventilation pressure column vector .

[0140] In the formula It is a symbolic function.

[0141] further, Figure 7 In and If the relationship is fitted as: Then the ventilation pressure column vector for:

[0142] = =

[0143] Because this situation introduces unknowns. Therefore, it is necessary to add an auxiliary solution function based on the relationship between the number of unknowns and the number of equations. Specifically, the method for establishing the auxiliary solution function is as follows: assume that the isolated node corresponding to the wind turbine is... Another entry and exit node of the wind network is The auxiliary solution function, as shown in Equation 2, is used to indicate the wind pressure relationship between the two inlets and outlets:

[0144] (2)

[0145] Scenario 4, such as Figure 9 The air network shown has a fixed inflow and outflow volume. In a ventilation network, a local ventilation fan is installed in a certain ventilation path (assumed to be a dedicated roadway). The ventilation pressure generated by the local ventilation fan is... The air volume through the fan is Based on the above parameter settings, establish the air volume compensation column vector and the local ventilation fan ventilation pressure column vector.

[0146] In this case, the methods for establishing the air volume compensation column vector include: according to the basic correlation matrix. B Node number v with line number r The correspondence was found, and the number was... The row number corresponding to the isolated node .Establish A column vector with 1 row and 1 column all zeros. If the isolated node is a wind network inlet, then the first zero in this column vector is... Rewrite the line as If the isolated node is the outlet of the wind network, then the first node in the column vector... Rewrite the line as Establish an air volume compensation column vector .

[0147] Specifically, if Figure 9 The basic correlation matrix of the wind network is obtained by deleting nodes. The corresponding row is obtained, so only the row needs to be checked. Perform airflow compensation; since this node is the outlet of the air network, its airflow compensation value should be [value missing]. Therefore, the airflow compensation column vector of this ventilation network is:

[0148]

[0149] In this case, the method for establishing the ventilation pressure column vector includes: fitting the wind pressure characteristic curve of the local fan to a vector based on the ventilation pressure of the main fan. and The expression generally represents a polynomial of degree no more than four. Based on the correspondence between branch number n and column number c, find the polynomial with the following number: The column number corresponding to the equivalent branch of the main fan Create an N x 1 column vector containing only zeros. Then, define the Nth column vector containing zeros. Rewrite the line as Establish main ventilation pressure column vector .

[0150] Specifically, Figure 9 In and If the relationship is fitted as: Then the ventilation pressure column vector for:

[0151] = =

[0152] Scenario 5 Figure 11 The wind network shown has two or more air inlets and outlets, and the total airflow into (or out of) the outside is a fixed value. Based on this parameter setting, airflow compensation is performed on the corresponding isolated nodes in the ventilation network diagram, an airflow compensation column vector is established, and an auxiliary solution function is created.

[0153] In this case, the methods for establishing the air volume compensation column vector include: according to the basic correlation matrix. B Node number v with line number r The correspondence was found, and the number was... The row number corresponding to the isolated node .Establish A column vector with 1 row and 1 column all zeros, if the isolated node corresponds to the wind network inlet, then the first zero in the column vector... Rewrite the row with the airflow value of that inlet, taking the negative value. If the isolated node corresponds to the outlet of the air network, then the row in the column vector... Rewrite the row with the air volume value at that entrance, taking positive numbers, and establish an air volume compensation column vector. .

[0154] Specifically, Figure 11 isolated nodes in As the entry node, As the outlet node of the wind network, its total inlet and outlet air volume is .assumed The air volume is ,but The air volume is ,if Figure 11 The basic correlation matrix of the wind network shown is obtained by deleting nodes. The corresponding row is obtained, therefore it is necessary to... , , Perform airflow compensation, with the following compensation values: , , Therefore, the airflow compensation column vector of this ventilation network is:

[0155]

[0156] In this case, the method for solving the auxiliary solution function is to add an auxiliary solution function of the form Equation 3 to indicate the wind pressure relationship between the two air inlet nodes.

[0157] (3)

[0158] Step S3: Modify the solution function of the nodal wind pressure equations and combine it with the auxiliary solution function to form the solution function of the equations to be solved. It should be noted that auxiliary solution functions are established in cases three and five; in other cases, it is not necessary to combine them with auxiliary solution functions to form the solution function of the equations to be solved. The following will provide a detailed explanation of each case with specific examples.

[0159] Scenario 1: Modify the solution function of the nodal wind pressure equations, including: for wind networks with isolated nodes, modify the solution function of the nodal air volume balance equation in the solution function of the nodal wind pressure equations. This modification is: adding an external air volume compensation vector. The corrected solution function for the nodal wind pressure equations is shown in Equation 4:

[0160]

[0161] (4)

[0162] Scenario 2: For wind networks with branches connected to only a single node, the solution function for the branch drag equation in the nodal wind pressure equation system is modified. This modification is: adding an external wind pressure compensation vector. .

[0163] The corrected solution function for the nodal wind pressure equations is shown in Equation 5:

[0164]

[0165] (5)

[0166] Scenario 3: For a ventilation network with only a main ventilation fan and undetermined main ventilation fan operating conditions, the solution function for the branch resistance in the nodal wind pressure equations is modified. This modification is: The external main ventilation pressure column vector... The solution function for the nodal airflow balance equation is modified, and this modification is: external airflow compensation vector. Together with the auxiliary solution function of the form 2 introduced, they constitute the solution function for the problem under study, as shown in Equation 6:

[0167]

[0168]

[0169] (6)

[0170] Scenario 4: For cases with only local ventilators and undetermined operating points, where the ventilation power of the air network is a fixed air volume, the branch resistance equation in the nodal wind pressure equation system is modified. This modification is as follows: External main ventilation pressure column vector... The solution function for the nodal airflow balance equation is modified, and this modification is: external airflow compensation vector. The corrected solution function for the nodal wind pressure equations is shown in Equation 7:

[0171]

[0172] (7)

[0173] Scenario 5: For a multi-inlet ventilation network with three wellheads, assuming the ventilation dynamic is a fixed air pressure, the solution function for the nodal air volume balance equation in the nodal air pressure equation system is modified. This modification is: an external air volume compensation vector. Together with the auxiliary solution function of the form shown in Equation 3, this constitutes the solution function for the problem under study; the modified solution function for the problem under study is shown in Equation 8:

[0174]

[0175]

[0176] (8)

[0177] Step S4: Preprocess the independent and dependent variables of the solution function of the system of equations to be solved. Specifically, based on the derivative of the dependent variable with respect to the independent variable, assemble the Jacobian matrix of the solution function of the system of equations. The following will provide a detailed explanation with specific examples.

[0178] Scenario 1, such as Figure 3 The ventilation network diagram shown places the wind pressure column vectors. P Air volume column vector Q The columns of independent variables X are concatenated to form the solution function for the problem under study. They take the form:

[0179] (9)

[0180] The dependent variables of the solution function for the problem under study are combined into a column vector Y of dependent variables (functions). It takes the form:

[0181] (10)

[0182] According to the definition of a Jacobian matrix, the partial derivatives of each component of the dependent variable (function) column vector Y with respect to each component of the independent variable column vector X are taken, and these partial derivative expressions are combined to form a Jacobian matrix. It takes the form:

[0183] (11)

[0184] Scenario 2 Figure 5 The ventilation network diagram shown places the wind pressure column vectors. P Air volume column vector Q The columns of independent variables X are concatenated to form the solution function for the problem under study. They take the form:

[0185] (12)

[0186] The dependent variables of the solution function for the problem under study are combined into a column vector Y of dependent variables (functions). It takes the form:

[0187] (13)

[0188] According to the definition of a Jacobian matrix, the partial derivatives of each component of the dependent variable (function) column vector Y with respect to each component of the independent variable column vector X are taken, and these partial derivative expressions are combined to form a Jacobian matrix. It takes the form:

[0189] (14)

[0190] Scenario 3, such as Figure 7 The ventilation network diagram shown places the wind pressure column vectors. P Air volume column vector Q and auxiliary unknowns The columns of independent variables X are concatenated to form the solution function for the problem under study. They take the form:

[0191] (15)

[0192] The dependent variables of the solution function for the problem under study are combined into a column vector Y of dependent variables (functions). It takes the form:

[0193] (16)

[0194] definition ,but Given an N x 1 column vector, this column vector corresponds to Location is All other positions are 0.

[0195] definition ,but For one The column vector in row 1 and column 1, which corresponds to The positions are -1 and 1 respectively, and the rest are 0.

[0196] definition ,but Given a row vector with 1 row and N columns, this row vector, excluding the corresponding... The positions are 1 and -1 respectively, and the rest are 0.

[0197] According to the definition of a Jacobian matrix, the partial derivatives of each component of the dependent variable (function) column vector Y with respect to each component of the independent variable column vector X are taken, and these partial derivative expressions are combined to form a Jacobian matrix. It takes the form:

[0198] (17)

[0199] Scenario 4, such as Figure 9 The ventilation network diagram shown places the wind pressure column vectors. P Air volume column vector Q The columns of independent variables X are concatenated to form the solution function for the problem under study. They take the form:

[0200] (18)

[0201] The dependent variables of the solution function for the problem under study are combined into a column vector Y of dependent variables (functions). It takes the form:

[0202] (19)

[0203] Defined, then Given a square matrix with N rows and N columns, except for the branches on the diagonal... The corresponding position is Except for the position , all other positions are 0. = = +...).

[0204] According to the definition of a Jacobian matrix, the partial derivatives of each component of the dependent variable (function) column vector Y with respect to each component of the independent variable column vector X are taken, and these partial derivative expressions are combined to form a Jacobian matrix. It takes the form:

[0205] (20)

[0206] Scenario 5, such as Figure 11The ventilation network diagram shown has wind pressure column vectors. P Air volume column vector Q and auxiliary unknowns The columns of independent variables X are concatenated to form the solution function for the problem under study. They take the form:

[0207] (twenty one)

[0208] The dependent variables of the solution function for the problem under study are combined into a column vector Y of dependent variables (functions). It takes the form:

[0209] (twenty two)

[0210] definition ,but For one The column vector in row 1 and column 1, which corresponds to The positions are -1 and 1 respectively, and the rest are 0.

[0211] definition ,but Given a row vector with 1 row and N columns, this row vector, excluding the corresponding... The positions are 1 and -1 respectively, and the rest are 0.

[0212] According to the definition of a Jacobian matrix, the partial derivatives of each component of the dependent variable (function) column vector Y with respect to each component of the independent variable column vector X are taken, and these partial derivative expressions are combined to form a Jacobian matrix. It takes the form:

[0213] (twenty three)

[0214] Step S5: Given the initial values ​​of the branch wind volume column vector Q, the nodal wind pressure column vector P, and other auxiliary unknowns, use Newton's method to solve the wind network problem using the solution function of the problem under study. Under the given termination conditions, obtain the approximate solution of the equation as the approximate value of the actual wind volume of the wind network.

[0215] Specifically, taking the solution function represented by Equation 8 as an example, step S5 includes the following steps:

[0216] 5.1 Choosing an arbitrary initial value , , according to Equation 21, combine the upper and lower parts to form a column vector .

[0217] 5.2, will , and The input parameters are fed into the solution function of the problem under study, which calculates the results according to Formula 8. , and Combine these three vectors according to Equation 22. As the output of the function that solves the system of equations.

[0218] 5.3, will and Assemble the Jacobian matrix using the form shown in Equation 23. .

[0219] 5.4. According to Newton's iterative method formula Calculate ,judge Does the value satisfy the given iteration exit condition? If it does, proceed to step 5.5. If it does not, Assigned Further Split into , Proceed to step 5.2.

[0220] 5.5, will Split into , This is the solution to the problem we are studying, which is the function we are trying to solve.

[0221] Based on the same inventive concept, this application also provides a solution apparatus for implementing the solution method described above. The solution provided by this apparatus is similar to the solution described in the above method; therefore, the specific limitations in one or more solution apparatus embodiments provided below can be found in the limitations of the solution method described above, and will not be repeated here.

[0222] In one embodiment, this application provides a solution apparatus, including: a setup module, an auxiliary module, a construction module, an assembly module, and a solution module.

[0223] The establishment module is used to: establish an equivalent ventilation network diagram of the complex wind network to be solved; establish a basic correlation matrix based on the equivalent ventilation network diagram; determine the branch wind resistance column vector based on the basic correlation matrix and the wind resistance values ​​of each branch; and establish a solution function for the nodal wind pressure equation system containing unknowns P and Q based on the basic correlation matrix and the branch wind resistance column vector.

[0224] The auxiliary module is used to: establish a correction vector for the solution function of the nodal wind pressure equations, and establish an auxiliary solution function based on the introduced unknowns other than P and Q.

[0225] The module is used to modify the solution function of the nodal wind pressure equation system and combine it with the auxiliary solution function to form the solution function of the equation system to be solved.

[0226] The assembly module is used to: preprocess the independent and dependent variables of the solution function of the system of equations to be solved. Specifically, it assembles the Jacobian matrix of the solution function of the system of equations based on the derivative of the dependent variable with respect to the independent variable.

[0227] The solution module is used to: given the initial values ​​of the branch wind volume column vector Q, the node wind pressure column vector P, and other auxiliary unknowns, use Newton's method to solve the wind network problem using the solution function of the problem under study, and obtain the approximate solution of the equation as the approximate value of the actual wind volume of the wind network under the given termination conditions.

[0228] Each module in the aforementioned calculation device can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in the processor of a computer device in hardware form or independent of it, or stored in the memory of a computer device in software form, so that the processor can call and execute the operations corresponding to each module.

[0229] In one embodiment, this application provides a computer device, which may be a server. The computer device includes a processor, memory, and a network interface connected via a system bus. The processor provides computing and control capabilities. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The network interface is used to communicate with external terminals via a network connection. When the computer program is executed by the processor, it implements a computational method.

[0230] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments described above. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM). The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0231] The following application will provide a complete explanation of specific examples of the above five scenarios.

[0232] Example 1

[0233] like Figure 2 As shown, this embodiment represents a local ventilation network in a mine. This network has two entrances / exits connected to the outside, and its internal geometry and physical parameters are known. The total air volume at the entrances / exits is fixed. Using the ventilation network calculation method described in this invention, the air volume information of each branch within the ventilation network can be calculated through the following steps.

[0234] Step S1: Establish an equivalent ventilation network diagram of the complex wind network to be solved. Establish a basic correlation matrix based on the equivalent ventilation network diagram. Determine the branch wind resistance column vector based on the basic correlation matrix and the wind resistance values ​​of each branch. Establish a solution function for the nodal wind pressure equation system containing unknowns P and Q based on the basic correlation matrix and the branch wind resistance column vector.

[0235] Step S101: Establish the equivalent ventilation network diagram of the complex air network to be solved; the specific process has been described above, and the established equivalent ventilation network diagram is as follows. Figure 3 As shown.

[0236] Step S102: Establish the basic correlation matrix B based on the equivalent ventilation network diagram. The specific process has been described above, including deleting nodes. The corresponding rows form the basic association matrix as follows:

[0237]

[0238] Step S103: Determine the branch drag column vector based on the basic correlation matrix and the drag values ​​of each branch. R The specific process has been described above. The following branch drag column vector is established, with units of [unit missing]. .

[0239]

[0240] The superscript T represents the transpose of the matrix.

[0241] Step S104, based on the basic association matrix B Branch drag column vector R Establish a system containing unknowns P and Q The solution function for the nodal wind pressure equations is described above.

[0242] Step S2: Establish the correction vector of the solution function for the nodal wind pressure equations, and establish an auxiliary solution function based on the introduced unknowns other than P and Q. This embodiment only involves the airflow compensation vector, the establishment process of which has been described above. The airflow compensation vector can be written in the following form, with units of [unit missing]. .

[0243]

[0244] Step S3: Modify the solution function of the nodal wind pressure equations and combine it with the auxiliary solution function to form the solution function of the equations to be solved. The specific solution function has been described above.

[0245] Step S4: Preprocess the independent and dependent variables of the solution function of the system of equations to be solved. Specifically, assemble the Jacobian matrix of the solution function of the system of equations based on the derivative of the dependent variable with respect to the independent variable. The specific formula has been described above.

[0246] Step S5: Given the initial values ​​of the branch wind volume column vector Q, the nodal wind pressure column vector P, and other auxiliary unknowns, use Newton's method to solve the wind network problem using the solution function of the problem under study. Under the given termination conditions, obtain the approximate solution of the equation as the approximate value of the actual wind volume of the wind network.

[0247] initial value Take a random column vector of 5 rows and 1 column. Take a random column vector with 7 rows and 1 column.

[0248] The final solution result in this embodiment is rounded to two decimal places. .

[0249]

[0250] Example 2

[0251] like Figure 4 As shown, this embodiment is a local ventilation network of a mine. Its internal structure and physical parameters are known. The inlet and outlet air pressures of the ventilation network are fixed by actual measurement. Using the ventilation network calculation method described in this invention, the air volume and other information of each branch inside the ventilation network can be calculated through the following steps.

[0252] Step S1: Establish an equivalent ventilation network diagram of the complex wind network to be solved. Establish a basic correlation matrix based on the equivalent ventilation network diagram. Determine the branch wind resistance column vector based on the basic correlation matrix and the wind resistance values ​​of each branch. Establish a solution function for the nodal wind pressure equation system containing unknowns P and Q based on the basic correlation matrix and the branch wind resistance column vector.

[0253] Step S101: Establish the equivalent ventilation network diagram of the complex air network to be solved; the specific process has been described above, and the established equivalent ventilation network diagram is as follows. Figure 5 As shown.

[0254] Step S102: Establish the basic correlation matrix B based on the equivalent ventilation network diagram. The specific process has been described above. The basic correlation matrix for this ventilation network diagram does not require deleting rows corresponding to nodes. The established basic correlation matrix is ​​as follows:

[0255]

[0256] Step S103: Determine the branch drag column vector based on the basic correlation matrix and the drag values ​​of each branch. R The specific process has been described above. The following branch drag column vector is established, with units of [unit missing]. .

[0257]

[0258] Step S104, based on the basic association matrixB Branch drag column vector R Establish a system containing unknowns P and Q The solution function for the nodal wind pressure equations is described above.

[0259] Step S2: Establish the correction vector for the solution function of the nodal wind pressure equations, and establish an auxiliary solution function based on the introduced unknowns other than P and Q. This embodiment only involves the wind pressure compensation vector, the establishment process of which has been described above. The air volume compensation vector can be written in the following form, with units of... .

[0260]

[0261] Step S3: Modify the solution function of the nodal wind pressure equations and combine it with the auxiliary solution function to form the solution function of the equations to be solved. The specific solution function has been described above.

[0262] Step S4: Preprocess the independent and dependent variables of the solution function of the system of equations to be solved. Specifically, assemble the Jacobian matrix of the solution function of the system of equations based on the derivative of the dependent variable with respect to the independent variable. The specific formula has been described above.

[0263] Step S5: Given the initial values ​​of the branch wind volume column vector Q, the nodal wind pressure column vector P, and other auxiliary unknowns, use Newton's method to solve the wind network problem using the solution function of the problem under study. Under the given termination conditions, obtain the approximate solution of the equation as the approximate value of the actual wind volume of the wind network.

[0264] initial value Take a random column vector of 4 rows and 1 column. Take a random column vector with 7 rows and 1 column.

[0265] The final solution result in this embodiment is rounded to two decimal places. .

[0266]

[0267] Example 3

[0268] like Figure 6 As shown, this embodiment represents a local ventilation network in a mine. This network has two inlets / outlets connected to the outside environment. Its internal structure and physical parameters are known. A main ventilation fan is installed on the branch connecting the network to the outside. The operating point of this fan is unknown and is given as a wind pressure characteristic curve. The ventilation network can be used to solve for parameters such as air volume through the following steps using the method described in this invention.

[0269] Step S1: Establish the equivalent ventilation network diagram of the complex wind network to be solved. Based on the equivalent ventilation network diagram, establish the basic correlation matrix. Based on the basic correlation matrix and the wind resistance value of each branch, determine the branch wind resistance column vector. Based on the basic correlation matrix and the branch wind resistance column vector, establish the solution function of the nodal wind pressure equation system containing unknowns P and Q.

[0270] Step S101: Establish the equivalent ventilation network diagram of the complex air network to be solved; the specific process has been described above, and the established equivalent ventilation network diagram is as follows. Figure 7 As shown.

[0271] Step S102: Establish the basic correlation matrix B based on the equivalent ventilation network diagram. The specific process has been described above, including deleting nodes. The corresponding rows form the basic association matrix as follows:

[0272]

[0273] Step S103: Determine the branch drag column vector based on the basic correlation matrix and the drag values ​​of each branch. R The specific process has been described above. The following branch drag column vector is established, with units of [unit missing]. .

[0274]

[0275] Step S104, based on the basic association matrix B Branch drag column vector R Establish a system containing unknowns P and Q The solution function for the nodal wind pressure equations is described above.

[0276] Step S2: Establish the correction vector of the solution function for the nodal wind pressure equations, and establish an auxiliary solution function based on the introduced unknowns other than P and Q. Unknowns are introduced in this embodiment. This represents the airflow at the inlet and outlet of the fan connected to the outside environment. An airflow compensation vector is established based on this unknown quantity. As mentioned earlier.

[0277] The wind pressure characteristic curve of the wind turbine is fitted to a polynomial function. In this embodiment, only the quadratic term was fitted, and the obtained parameter values ​​are: a=3100, b=-0.12, c=-0.035.

[0278] Write the fan pressure vector In actual solution, it is written as:

[0279]

[0280] In the formula It is a symbolic function.

[0281] Because there is an additional independent variable This leads to an inconsistency between the number of independent variables and the number of equations. In this invention, an auxiliary equation is introduced to represent the wind pressure relationship between the air outlets and inlets of the air supply network. In this embodiment, it is assumed that the wind pressure at the two inlets and outlets is equal. The specific form of the auxiliary solution function has been described above.

[0282] Step S3: Modify the solution function of the nodal wind pressure equations and combine it with the auxiliary solution function to form the solution function of the equations to be solved. The specific solution function has been described above.

[0283] Step S4: Preprocess the independent and dependent variables of the solution function of the system of equations to be solved. Specifically, assemble the Jacobian matrix of the solution function of the system of equations based on the derivative of the dependent variable with respect to the independent variable. The specific formula has been described above.

[0284] in

[0285] Step S5: Given the initial values ​​of the branch wind volume column vector Q, the nodal wind pressure column vector P, and other auxiliary unknowns, use Newton's method to solve the wind network problem using the solution function of the problem under study. Under the given termination conditions, obtain the approximate solution of the equation as the approximate value of the actual wind volume of the wind network.

[0286] initial value Take a random column vector of 6 rows and 1 column. Take a random column vector of 8 rows and 1 column. Generate a single random number.

[0287] The final solution result in this embodiment is rounded to two decimal places.

[0288]

[0289] =254.20

[0290] Example 4

[0291] like Figure 8 As shown, this embodiment represents a local ventilation network in a mine. The internal structure and physical parameters of this ventilation network are known, and the airflow at the entrances and exits connecting the network to the outside is fixed. A local ventilator, occupying a dedicated roadway, is installed for auxiliary ventilation. The local ventilator is represented by a wind pressure characteristic curve. According to the present invention, parameters such as airflow within this ventilation network can be determined through the following steps.

[0292] Step S1: Establish the equivalent ventilation network diagram of the complex wind network to be solved. Based on the equivalent ventilation network diagram, establish the basic correlation matrix. Based on the basic correlation matrix and the wind resistance value of each branch, determine the branch wind resistance column vector. Based on the basic correlation matrix and the branch wind resistance column vector, establish the solution function of the nodal wind pressure equation system containing unknowns P and Q.

[0293] Step S101: Establish the equivalent ventilation network diagram of the complex air network to be solved; the specific process has been described above, and the established equivalent ventilation network diagram is as follows. Figure 9 As shown.

[0294] Step S102: Establish the basic correlation matrix B based on the equivalent ventilation network diagram. The specific process has been described above, including deleting nodes. The corresponding rows form the basic association matrix as follows:

[0295]

[0296] Step S103: Determine the branch drag column vector based on the basic correlation matrix and the drag values ​​of each branch. R The specific process has been described above. The following branch drag column vector is established, with units of [unit missing]. .

[0297]

[0298] Step S104, based on the basic association matrix B Branch drag column vector R Establish a system containing unknowns P and Q The solution function for the nodal wind pressure equations is described above.

[0299] Step S2: Establish the correction vector of the solution function for the nodal wind pressure equation system, and establish the auxiliary solution function based on the introduced unknowns other than P and Q.

[0300] In this embodiment, the total air volume entering and exiting the air network is assumed to be 100, and the air volume compensation column vector can be written as follows:

[0301]

[0302] The wind pressure characteristic curve of the wind turbine is fitted to a polynomial function. In this embodiment, only the quadratic term was fitted, and the obtained parameter values ​​are: a=3100, b=-0.12, c=-0.035.

[0303] Write the fan pressure vector In actual solution, it is written as:

[0304]

[0305] Step S3: Modify the solution function of the nodal wind pressure equations and combine it with the auxiliary solution function to form the solution function of the equations to be solved. The specific solution function has been described above.

[0306] Step S4: Preprocess the independent and dependent variables of the solution function of the system of equations to be solved. Specifically, assemble the Jacobian matrix of the solution function of the system of equations based on the derivative of the dependent variable with respect to the independent variable. The specific formula has been described above.

[0307] in

[0308]

[0309] Step S5: Given the initial values ​​of the branch wind volume column vector Q, the nodal wind pressure column vector P, and other auxiliary unknowns, use Newton's method to solve the wind network problem using the solution function of the problem under study. Under the given termination conditions, obtain the approximate solution of the equation as the approximate value of the actual wind volume of the wind network.

[0310] initial value Take a random column vector of 7 rows and 1 column. Take a random column vector with 9 rows and 1 column.

[0311] The final solution result in this embodiment is rounded to one decimal place. .

[0312]

[0313] Example 5

[0314] like Figure 10 As shown, this embodiment represents a local ventilation network in a mine. Its internal structure and physical parameters are known. The network has three inlet and outlet ports connected to the outside (one inlet is a hypothetical leakage inlet), and the outlet airflow is fixed. The network does not have a fan. Its airflow and other parameters can be obtained through the steps described in this invention.

[0315] Step S1: Establish the equivalent ventilation network diagram of the complex wind network to be solved. Based on the equivalent ventilation network diagram, establish the basic correlation matrix. Based on the basic correlation matrix and the wind resistance value of each branch, determine the branch wind resistance column vector. Based on the basic correlation matrix and the branch wind resistance column vector, establish the solution function of the nodal wind pressure equation system containing unknowns P and Q.

[0316] Step S101: Establish the equivalent ventilation network diagram of the complex air network to be solved; the specific process has been described above, and the established equivalent ventilation network diagram is as follows. Figure 11 As shown.

[0317] Step S102: Establish the basic correlation matrix B based on the equivalent ventilation network diagram. The specific process has been described above, including deleting nodes. The corresponding rows form the basic association matrix as follows:

[0318]

[0319] Step S103: Determine the branch drag column vector based on the basic correlation matrix and the drag values ​​of each branch. R The specific process has been described above. The following branch drag column vector is established, with units of [unit missing]. .

[0320]

[0321] Step S104, based on the basic association matrix B Branch drag column vector R Establish a system containing unknowns P and Q The solution function for the nodal wind pressure equations is described above.

[0322] Step S2: Establish the correction vector of the solution function for the nodal wind pressure equation system, and establish the auxiliary solution function based on the introduced unknowns other than P and Q.

[0323] Introducing an unknown quantity Assume that the unknown represents a node. The compensation air volume, assuming the outlet If the compensation air volume is -100, then the node The compensation air volume is 100- Based on this, write the compensation air volume column vector for this wind network:

[0324]

[0325] Step S3: Modify the solution function of the nodal wind pressure equations and combine it with the auxiliary solution function to form the solution function of the equations to be solved. The specific solution function has been described above.

[0326] Step S4: Preprocess the independent and dependent variables of the solution function of the system of equations to be solved. Specifically, assemble the Jacobian matrix of the solution function of the system of equations based on the derivative of the dependent variable with respect to the independent variable. The specific formula has been described above.

[0327] Step S5: Given the initial values ​​of the branch wind volume column vector Q, the nodal wind pressure column vector P, and other auxiliary unknowns, use Newton's method to solve the wind network problem using the solution function of the problem under study. Under the given termination conditions, obtain the approximate solution of the equation as the approximate value of the actual wind volume of the wind network.

[0328] initial value Take a random column vector of 6 rows and 1 column. Take a random column vector of 8 rows and 1 column. Choose a random number.

[0329] The final solution result in this embodiment is rounded to two decimal places. .

[0330]

[0331] =45.15

[0332] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0333] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A method for solving a complex ventilation network, characterized in that, The method includes: Establish an equivalent ventilation network diagram of the complex wind network to be solved, establish a basic correlation matrix based on the equivalent ventilation network diagram, determine the branch wind resistance column vector based on the basic correlation matrix and the wind resistance values ​​of each branch, and establish a solution function for the nodal wind pressure equation system containing unknowns P and Q based on the basic correlation matrix and the branch wind resistance column vector. Establish a correction vector for the solution function of the nodal wind pressure equations, and establish an auxiliary solution function based on the introduced unknowns other than P and Q; The solution function of the nodal wind pressure equation system is modified and combined with the auxiliary solution function to form the solution function of the equation system to be solved. The independent and dependent variables of the solution function of the system of equations to be solved are preprocessed. Specifically, the Jacobian matrix of the solution function of the system of equations is assembled based on the derivative of the dependent variable with respect to the independent variable. Given the initial values ​​of the branch wind volume column vector Q, the nodal wind pressure column vector P, and other auxiliary unknowns, the wind network is solved using Newton's method for the solution function of the problem under study. Under the given termination conditions, the approximate solution of the equation is obtained as the approximate value of the actual wind volume of the wind network. The correction vector for establishing the solution function of the nodal wind pressure equation system includes: The air volume compensation column vector in the case one, case three, case four and case five The air pressure compensation column vector in the case two The main fan ventilation pressure compensation column vector in the case three The local fan ventilation pressure compensation column vector in the case four The auxiliary solving function is established according to the introduced unknown quantity except P and Q, specifically, the auxiliary solving functions with the same name are established in the case three and case five ​ The following is stated: Scenario 1: The inflow or outflow of air to the outside of the complex wind network to be calculated is a fixed value; Scenario 2 is: the pressure difference between the inlet and outlet of the complex wind network to be solved and the outside world is constant; Scenario 3 is as follows: The main ventilation fan is installed at the connection point between the complex ventilation network to be solved and the outside world, and the operating point of the fan is yet to be determined; Scenario 4 is: the complex air network to be solved has local ventilation fans on its branches; Scenario 5 is when the complex wind network to be solved has more than two entrances and exits to the outside world; Case 3 establishes an auxiliary solution function: for a wind network with only a main ventilator and an undetermined main ventilator operating condition, the solution function for the branch resistance equation in the nodal wind pressure equation set is modified; Case 5 establishes an auxiliary solution function: For a multi-inlet ventilation network with three wellheads, assuming the ventilation dynamic is a fixed wind pressure, the solution function for the nodal air volume balance equation in the nodal wind pressure equation set is modified; Other auxiliary unknown quantities are the air volume passing through the fan and the air intake of one of the three intake air shafts of the two-in-one-out three-intake-air-shaft system .

2. The method of claim 1, wherein, The process of establishing an equivalent ventilation network diagram for a complex wind network to be solved includes: equating a single wind path or multiple wind paths in series in the wind network to a branch, and equating the intersection of two or more wind paths to a node; then, assuming the direction of airflow, making all branches into directed branches; and finally, connecting the branches and nodes to construct an equivalent ventilation network diagram.

3. The method of claim 1, wherein, The establishment of the basic correlation matrix based on the equivalent ventilation network diagram includes: According to the node number M and branch number N in the equivalent ventilation network diagram, a full zero matrix A of M rows and N columns is established, and all nodes are recorded in a node number set ; the node number v in the equivalent ventilation network diagram and the row number r of A, and the branch number e in the equivalent ventilation network diagram and the column number c of A establish a corresponding relationship; if the node numbered is located at the tail of the branch numbered , according to the corresponding relationship, the number in the corresponding row column position of the matrix A is rewritten as 1; if the node numbered is located at the head of the branch numbered , according to the corresponding relationship, the number in the corresponding row column position of the matrix A is rewritten as -1; rank of matrix A is obtained , if =M, A is assigned to B; if <M, A is deleted corresponding number of rows according to the difference between and M, and the node numbers corresponding to the rows are deleted from the node number set ; let the obtained new matrix be , the number of rows of which is , , the rank of which is , if = = , the new matrix is assigned to B; at this time, B is the basic incidence matrix, the number of rows of which is recorded as , and the number of columns remains N; According to the basic incidence matrix B, a node number set is established The correspondence between the node number v and the row number r of B; the correspondence between the branch number e of the inheritance A and the column number c is the correspondence of the matrix B.

4. The method of claim 3, wherein, The step of determining the branch drag column vector based on the basic correlation matrix and the drag values ​​of each branch includes: establishing an N-row, 1-column column vector based on the number of columns N of the basic correlation matrix B; filling the drag value of branch e in row c of the column vector according to the correspondence between branch number e and column number c; and obtaining the branch drag column vector R.

5. The method of claim 4, wherein, The solution function for the nodal wind pressure equation system containing unknowns P and Q, based on the basic correlation matrix B and the branch wind resistance column vector R, is as follows: （1a） （1b） where Q is a N-row 1-column branch air volume column vector, P is a N-row 1-column node air pressure column vector; subscript diag represents performing a diagonalization operation on a vector to become a matrix, superscript T represents the transpose of a matrix, and symbol || represents taking an absolute value; formula 1a is called a branch resistance equation solving function, and formula 1b is called a node air volume balance equation solving function.

6. A solving device characterized by comprising: include: The establishment module is configured to establish an equivalent ventilation network diagram of a complex ventilation network to be solved, establish a basic incidence matrix according to the equivalent ventilation network diagram, determine a branch resistance column vector according to the basic incidence matrix and the wind resistance values of the branches, and establish a node wind pressure equation set solving function containing unknowns P and Q according to the basic incidence matrix and the branch resistance column vector; The auxiliary module is configured to establish a correction vector of the node wind pressure equation set solving function, and establish an auxiliary solving function according to the introduced unknowns other than P and Q; The construction module is configured to perform correction processing on the node wind pressure equation set solving function, and combine the auxiliary solving function to form an equation set solving function to be solved; The assembly module is configured to perform preprocessing on independent variables and dependent variables of the equation set solving function to be solved, and specifically, assemble a Jacobian matrix of the equation set solving function according to a derivative function formula of the dependent variables with respect to the independent variables; The solving module is configured to give initial values of a branch wind volume column vector Q, a node wind pressure column vector P and other auxiliary unknowns, solve the solving function of the studied problem by using a Newton method, and obtain an approximate solution of the equation as an approximate value of the actual wind volume of the ventilation network under a given end condition. The correction vector of the node wind pressure equation set solving function includes: The air volume compensation column vector in the case one, case three, case four and case five The air pressure compensation column vector in the case two The main fan ventilation pressure compensation column vector in the case three The local fan ventilation pressure compensation column vector in the case four The auxiliary solving function is established according to the introduced unknown quantity except P and Q, and specifically, the auxiliary solving functions with the same name are established in the case three and case five ​ The correction vector of the node wind pressure equation set solving function includes: Case one is that the inflow or outflow of the complex ventilation network to the outside world is a fixed value; Case two is that the pressure difference between the entrance and exit of the complex ventilation network connected with the outside world is constant; Case three is that a main ventilation fan is arranged at the connection point of the complex ventilation network with the outside world, and the working point of the fan is to be determined; Case four is that a local ventilation fan is arranged on a branch in the complex ventilation network; Case five is that the complex ventilation network has more than two entrances to the outside world; The auxiliary solving function is established for case three: for the ventilation network in which only a main ventilation fan exists and the working point of the main ventilation fan is undetermined, the branch resistance equation solving function in the node wind pressure equation set solving function is corrected; The auxiliary solving function is established for case five: for the multi-entrance ventilation network, in the case of three wellheads, the ventilation power is assumed to be a fixed wind pressure, and the node wind volume balance equation solving function in the node wind pressure equation set solving function is corrected. Other auxiliary unknown quantities are the air volume passing through the fan and the air intake of one of the three intake air shafts of the two-in-one-out three-intake-air-shaft system . 7.A computer device, comprising a memory and a processor, wherein the memory stores a computer program, and the computer device is configured to perform the method according to any one of claims 1-6 when the computer program is executed by the processor. The computer program is executed by the processor to realize the steps of the method of any one of claims 1 to 6.

8. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to realize the steps of the method of any one of claims 1 to 6.

Images