CN114692524B - Wind tunnel group high-pressure air resource dynamic topological structure modeling method, system and air consumption calculation - Google Patents
Wind tunnel group high-pressure air resource dynamic topological structure modeling method, system and air consumption calculation Download PDFInfo
- Publication number
- CN114692524B CN114692524B CN202210335439.4A CN202210335439A CN114692524B CN 114692524 B CN114692524 B CN 114692524B CN 202210335439 A CN202210335439 A CN 202210335439A CN 114692524 B CN114692524 B CN 114692524B
- Authority
- CN
- China
- Prior art keywords
- wind tunnel
- pressure air
- node
- test
- group
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000012360 testing method Methods 0.000 claims abstract description 62
- 238000005192 partition Methods 0.000 claims description 40
- 238000009826 distribution Methods 0.000 claims description 37
- 239000011159 matrix material Substances 0.000 claims description 29
- 238000010586 diagram Methods 0.000 claims description 14
- 238000004891 communication Methods 0.000 claims description 12
- 238000012544 monitoring process Methods 0.000 claims description 5
- 239000007789 gas Substances 0.000 description 11
- 230000005540 biological transmission Effects 0.000 description 4
- 238000003860 storage Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000000498 cooling water Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 238000010992 reflux Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention is suitable for the field of wind tunnel group test operation scheduling, and particularly relates to a wind tunnel group high-pressure air resource dynamic topological structure modeling method, a system and air consumption calculation. The invention updates and records the topological structure of the tank group and the pipeline in real time, acquires and guarantees the element set of a specific wind tunnel and the sharing relation between the wind tunnel and other wind tunnels, and improves accurate data support for formulating test strategies and calculating the wind tunnel resource consumption; avoiding the inconvenience and inefficiency of manual control and coordination.
Description
Technical Field
The invention relates to the field of wind tunnel group test operation scheduling, in particular to a wind tunnel group high-pressure air resource dynamic topological structure modeling method and system and air consumption calculation.
Background
The large wind tunnel equipment facility is a national major strategic resource, the wind tunnel group generally comprises tens or even tens of production type wind tunnels, and each wind tunnel test needs various power resources such as pure water, electric power, high-pressure air, medium-pressure air, vacuum, nitrogen, hydrogen, oxygen and the like. To meet the test requirements of the wind tunnel group, the configuration of the power system is often very complex.
Taking a high-pressure air system as an example, the high-pressure air system is divided into 32MPa and 22MPa high-pressure air systems, and main equipment comprises a reciprocating compressor, a dryer, a 32MPa high-pressure storage tank, a 22MPa high-pressure storage tank, cooling water and a gas distribution network system. The distribution network system consists of a plurality of valves, pipelines, accessory equipment and the like which are involved in the production, storage, transportation of high-pressure air from the compressor unit to the high-pressure tank group and the transportation of the high-pressure air to the wind tunnel. Because the wind tunnel groups share high-pressure air resources in real time, different wind tunnel test requirements are ensured by switching and combining different tank groups, valves and pipelines.
The high-pressure air source can be divided into different tank groups (subareas) by considering the construction cost and the high-pressure air resource requirements of different dynamic changes of the wind tunnel test, wherein the different tank groups are formed by communicating a plurality of tank bodies through valves and pipelines; the air is conveyed from the tank group to different wind tunnels, wherein each air conveying main pipe is divided into a plurality of branch pipelines, and each branch pipeline can generally guarantee a plurality of wind tunnels; and a communication valve is also generally arranged between the gas transmission main pipes, so that the combined supply of high-pressure air resources of different wind tunnels and different tests is realized. The physical structure reflecting the distribution relation, which consists of a high-pressure compressor set, a valve, a tank group, a gas transmission main pipe, a branch pipeline and a wind tunnel under guarantee implementation, is called a dynamic topology structure of a high-pressure air resource distribution network, and the change of the working state of each node device in the topology means the change of the power resource supply configuration (pressure level, standard square capacity and the like), as shown in fig. 1 and 2.
The high-frequency and multi-type pneumatic tests carried out by the wind tunnel group cause the complexity of high-pressure air resource supply guarantee, and the situations that a plurality of wind tunnels share the same tank group in a time sharing manner, a plurality of wind tunnels share different tank groups in a parallel manner, a plurality of tank groups are used in a large wind tunnel simultaneously and the like exist, so that great challenges are provided for safe operation of the wind tunnel test, timely guarantee of power resources, reasonable allocation, accurate consumption metering and the like.
On the one hand, if a plurality of wind tunnels of the wind tunnel group are started simultaneously, and the same tank group and the same pipeline are used, resonance phenomenon is easy to occur, the safety of the whole high-pressure system and wind tunnel equipment is threatened, the quality of a wind tunnel flow field is reduced, and the quality of test data is influenced. In order to avoid the occurrence of high-pressure air resource conflict, the current wind tunnel test operation mostly adopts a manual coordination mode, communication is carried out through telephone, and the scheduling mode has the problems of incomplete monitoring of the topology structure information of the high-pressure air resource distribution network, insufficient intelligent degree of power resource scheduling, low coordination efficiency and the like, and brings certain potential safety hazard.
On the other hand, in a wind tunnel group sharing high-pressure air resources, two schemes of wind tunnel pipeline additionally provided with a flowmeter or tank group differential pressure metering before and after the test are generally adopted to realize consumption metering of a single wind tunnel and a single test. However, because the tank group pressure changes rapidly during the wind tunnel group test, the wind tunnel pipeline is not accurately monitored by adding the flowmeter, the cost is high, and the flow field of the wind tunnel test can be influenced; by adopting the tank group differential pressure metering scheme, the topological structures of the tank groups and pipelines need to be mastered in real time so as to acquire the volumes (standard square capacity) of the gas supply tank groups and pipelines of the wind tunnel under test to calculate the high-pressure air consumption of the wind tunnel test.
Therefore, the management of the high-pressure air resource distribution network topology structure aiming at the complex and dynamic wind tunnel group is developed, and the method has very important significance for avoiding power resource conflict and realizing wind tunnel test consumption metering.
However, no public report is known for the topology structure management of the high-pressure air resource distribution network of the complex and dynamic wind tunnel group at present.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a wind tunnel group high-pressure air resource distribution network dynamic topology structure modeling method based on a graph database and an adjacent matrix, which establishes a model for automatically identifying, storing and updating the dynamic high-pressure air resource distribution network topology structure and is convenient for the subsequent calculation of air resource consumption.
The wind tunnel group high-pressure air resource dynamic topological structure modeling method is characterized by comprising the following steps of:
s10, acquiring the field operation state of a power system
The power system comprises test dynamics of each wind tunnel of the wind tunnel group, tank group dynamics, compressor unit dynamics and valve states;
s20, acquiring a high-pressure air distribution network adjacency matrix A n×n
Taking a high-pressure air resource distribution network as a directed graph, taking a compressor unit, a pipeline airflow converging point, a tank group and a wind tunnel in the distribution network as nodes, and taking a valve in the distribution network as an edge to obtain an adjacent matrix:
wherein e i For the ith node, e j I epsilon [1, n for j-th node],j∈[1,n]；
S30, calculating a reachable matrix M
S40, decomposing the reachable matrix to obtain the communication branches of the gas distribution network diagram
The distribution network diagram is divided into a plurality of independent, direct or indirect influence free communication branches.
Further, step S40 includes the steps of:
s41, establishing a reachable set R (e i )
Wherein, 
for any one node e i V is the set of nodes;
s42, establishing a look-ahead set Q (e) i )
S43, establishing a bottom node set B of the system
B＝{e i |e i E V and R (e) i )∩Q(e i )＝Q(e i )}
S44, carrying out region division on the distribution network diagram
For any node B, B' in the underlying node set B, if
R(b)∩R(b′)＝φ
The nodes b, b' belong to different areas, so that the system V can be divided into H partitions, which are recorded as
π(V)＝{V 1 ,V 2 ,…,V h ,…，V H }
Wherein each partition is a communication branch.
Further, after step S44, the method further includes the steps of:
s45, establishing a top-level unit set T in each partition
In the subarea
In the process,
T＝{e i |e i ∈V h and R (e) i )∩Q(e i )＝R(e i )}
S46, searching a wind tunnel w k Partition V being the top-level element h The method is characterized by comprising the steps of guaranteeing the collection of a compressor unit, a tank group, a valve and a gas pipeline of the wind tunnel test, wherein the combination of elements and connection relations of the elements in the collection is the dynamic topological structure of the wind tunnel.
A system for executing the wind tunnel group high-pressure air resource dynamic topological structure modeling method is characterized by comprising hardware, an interface and software,
the hardware comprises a detector for monitoring the test dynamic state of the wind tunnel, the tank group dynamic state, the compressor unit dynamic state and the valve state; the interface adopts an OPCserver to upload the information acquired by the hardware to a central server, and the information is processed and visualized by software.
Further, each wind tunnel comprises a local controller, the local controller sends a request signal for requesting test to a central server, the central server makes a test plan according to a scheduling rule, sends a signal for allowing test or not allowing test to the local controller, and the local controller controls the opening or closing of a valve in the wind tunnel according to the feedback signal.
Further, the scheduling rule is that only one wind tunnel is guaranteed to use high-pressure air resources in each partition at the same time.
The invention also provides a method for calculating the consumption of the high-pressure air resources of the wind tunnel group, which is characterized in that a certain wind tunnel w is ensured based on the method for modeling the dynamic topological structure of the high-pressure air resources of the wind tunnel group k Is a dynamic topology of:
calculating and ensuring the wind tunnel w k Sum of volumes Vol of all elements of (C) w ，
This wind tunnel w k High-pressure air resource consumption P of (2) k The method comprises the following steps:
P k tank group pressure at the beginning of the wind tunnel test-tank group pressure at the end of the wind tunnel test-i x Vol = w ×10。
The wind tunnel group high-pressure air resource dynamic topological structure modeling method, the system and the air consumption calculation have at least the following beneficial effects compared with the prior art:
(1) According to the modeling method disclosed by the invention, pipelines and tank groups which are used in a specific wind tunnel can be automatically identified and updated in real time;
(2) The invention can update and record the topological structure of tank groups and pipelines in real time, acquire and ensure the element set of a specific wind tunnel and the sharing relation between the wind tunnel and other wind tunnels, and improve accurate data support for formulating test strategies and calculating the wind tunnel resource consumption; the problems of inconvenience and low efficiency of manual control and coordination are avoided;
(3) By using the modeling method provided by the invention, the high-pressure air resource consumption of a specific wind tunnel can be accurately calculated.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following description will briefly explain the embodiments of the present invention or the drawings used in the description of the prior art, and it is obvious that the drawings described below are only some embodiments of the present invention, and other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic diagram of a dynamic topology of a high-pressure air resource;
FIG. 2 is a schematic diagram of a dynamic topology of a high-pressure air resource;
FIG. 3 is a flow chart of a method for modeling a dynamic topology of a high-pressure air resource of a wind tunnel group according to an embodiment of the invention;
fig. 4 is an example of a method for modeling a dynamic topology structure of a high-pressure air resource of a wind tunnel group according to an embodiment of the present invention.
Detailed Description
The following description provides many different embodiments, or examples, for implementing different features of the invention. The elements and arrangements described in the following specific examples are presented for purposes of brevity and are provided only as examples and are not intended to limit the invention.
Example 1
A wind tunnel group high-pressure air resource dynamic topological structure modeling method, as shown in figure 3, comprises the following steps:
s10, acquiring the field operation state of a power system
The power system comprises test dynamics of each wind tunnel of the wind tunnel group, tank group dynamics, compressor unit dynamics and valve states;
s20, acquiring a high-pressure air distribution network adjacency matrix A n×n
Taking a high-pressure air resource distribution network as a directed graph, taking a compressor unit, a pipeline airflow converging point, a tank group and a wind tunnel in the distribution network as nodes, and taking a valve in the distribution network as an edge to obtain an adjacent matrix:
wherein e i For the ith node, e j I epsilon [1, n for j-th node],j∈[1,n]The method comprises the steps of carrying out a first treatment on the surface of the Element a of adjacency matrix a ij =1 indicates the connection relationship between the node i and the node j, 0 is non-contiguous, and 1 is contiguous. The kth row and the kth column have the same physical meaning and represent the node k, so that all nodes adjacent to the node k need to be found to mainly find whether elements with 1 exist in the kth row or the kth column.
In the invention, the valve is used as an edge because the pipeline determines the connection relation of a gas distribution network, and valve equipment controls switching of tank groups, units and pipelines, the opening and closing of the valve can cause the change of a physical connection model and the change of the adjacent relation of valve nodes, and the state of the valve is sampled in real time, so that the corresponding adjacent state can be obtained, and an adjacent matrix is obtained;
further, to describe the adjacency relationship between nodes of the gas distribution network, the present embodiment creates a graph database containing three relationship tables:
1) Establishing a valve-reflux node table: the table records the actual position relationship between the valve and the manifold node, and the table structure is shown in table 1 (the initial node and the final node in the table can be interchanged).
TABLE 1 valve-confluence node Meter
| Valve number | Start node | Terminating point |
| Val1 | GDx | GDy |
| Val2 | GDp | GDz |
| …… | …… | …… |
2) Establishing a node marking table: the node numbers are uniformly encoded with consecutive integers, corresponding to the subscripts of the adjacency matrix, and the table structure is shown in table 2.
Table 2 node representation table
3) Establishing a valve state table: and (3) operating the monitoring system through the wind tunnel and the high-pressure air system, acquiring the states of all valves in all high-pressure air distribution networks, and recording by using a table. The table records the open/close states of all valves at a certain time, and the table structure is shown in table 3. For a closed valve, the state value is 0; for an open valve, the state value is 1.
TABLE 3 valve State Table
| Time | Valve number | Status value |
| 2021-12-27-12-02-01 | Val1 | 0 |
| 2021-12-27-12-02-01 | |
1 |
| …… | …… | …… |
In the graphic database, the condition inquiry is carried out among the three tables of the valve-confluence node table, the node identification table and the valve state table to generate an adjacency matrix table, and the table structure is shown in table 4.
TABLE 4 adjacency matrix table
| Status value | Initial node identification number | Terminal identification number |
| 0 | 1 | 2 |
| 1 | 3 | 4 |
| …… | …… | …… |
S30, calculating an reachable matrix M to obtain the sum of direct and indirect relations of each node
The reachable matrix refers to the degree of accessibility between nodes of the graph after passing through a certain length of path, and is represented by M, and is also an n×n square matrix:
the ith row and jth column elements M of the reachable matrix M ij =1, meaning that node i has a direct (pipelined direct connection, i points to j) or indirect (through at least one intermediate node) connection relationship with node j, otherwise equal to 0. For example, if i points to g and g points to j, i and j have a 1-step indirect relationship; in the same way, there are 2 steps, 3 steps, and..indirect relations, whether there is a direct or an indirect relation, the corresponding element m ij =1, otherwise 0.
S40, decomposing the reachable matrix to obtain the communication branches of the gas distribution network diagram
The distribution network diagram is divided into a plurality of independent, direct or indirect influence free communication branches.
The method comprises the following steps:
s41, establishing a reachable set R (e i )
e i Having direct or indirect influence on a certain node, thenThese nodes form e i Can be reached by (a):
wherein, 
for any one node e i V is the set of nodes;
s42, establishing a look-ahead set Q (e) i )
Any one node pair e i With direct or indirect influence, then these nodes constitute e i Is a antecedent set of (a):
s43, establishing a bottom layer node set B
The set of nodes that any node points to only it and not to it is satisfied that there is no node that it points to, the underlying node set B that constitutes system V:
B＝{e i |e i e V and R (e) i )∩Q(e i )＝Q(e i )}
S44, carrying out region division on the distribution network diagram
For any node B, B' in the underlying node set B, if
R(b)∩R(b′)＝φ
The nodes b, b' cannot point to the same node and belong to different areas, so that the system V can be divided into H partitions, which are recorded as
π(V)＝{V 1 ,V 2 ,…,V h ,…，V H }
Wherein each partition is a communication branch.
Further, performing step S45 and step S46 may obtain a guarantee of the dynamic topology of a specific wind tunnel:
s45, establishing a top-level unit set T in each partition h
In the subarea
In the method, a set of nodes of which no node is only affected by the node but not affected by the node is satisfied to form a top level unit set T h ：
T h ＝{e i |e i ∈V h And R (e) i )∩Q(e i )＝R(e i )}
S46, searching a wind tunnel w k Partition V being the top-level element h The collection of the compressor unit, the tank group, the valve and the gas transmission pipeline of the wind tunnel test is ensured, and the dynamic topological structure of the wind tunnel is particularly ensured.
Example 2
FIG. 2 is a schematic diagram of the physical connection of a high pressure air resource, the system having three compressor banks: compressor unit 01, compressor unit 02 and compressor unit 03, four tank groups: tank group 01, tank group 02, tank group 03 and tank group 04, three main pipes: main pipe 1, main pipe 2 and main pipe 3, and four wind tunnels: wind tunnel 01, wind tunnel 02, wind tunnel 03 and wind tunnel 04, it should be noted that the real high-pressure air resource system is more complex, and fig. 2 is only an example of simplification, wherein valves are simplified, and a valve (i.e. on a side line) is arranged between each two nodes. During the test, corresponding valves are opened according to the test requirements, and the wind tunnel group high-pressure air resource dynamic topological structure modeling method is constructed by the following method:
s00, acquiring a system connection relation diagram
S01, establishing a node identification table
The node numbers are set according to the arrangement order of the compressor group, the tank group, the main pipe and the wind tunnel, namely, the compressor group 01, the compressor group 02, the compressor group 03, the tank group 01, the tank group 02, the wind tunnel 03 and the wind tunnel 04 are respectively set as numbers 1,2, 3, 4,5, the numbers of the following tables are shown in the following table:
s02, establishing a valve-node table:
| valve number | Start node | Final node | Valve number | Start node | |
| Val1 | |||||
| 1 | 4 | |
6 | 8 | |
| |
1 | 5 | |
6 | 9 |
| |
1 | 6 | |
6 | 10 |
| |
1 | 7 | |
7 | 8 |
| |
2 | 4 | |
7 | 9 |
| |
2 | 5 | |
7 | 10 |
| |
2 | 6 | |
8 | 11 |
| |
2 | 7 | |
8 | 12 |
| |
3 | 4 | |
8 | 13 |
| |
3 | 5 | |
8 | 14 |
| |
3 | 6 | |
9 | 11 |
| |
3 | 7 | |
9 | 12 |
| |
4 | 8 | |
9 | 13 |
| |
4 | 9 | |
9 | 14 |
| |
4 | 10 | |
10 | 11 |
| Val16 | 5 | 8 | Val34 | 10 | 12 |
| Val17 | 5 | 9 | Val35 | 10 | 13 |
| Val18 | 5 | 10 | |
10 | 14 |
The physical connection relation of the wind tunnel group high-pressure air resource system is obtained, the wind tunnel group high-pressure air resource distribution network is regarded as an undirected graph G= < V, E >, wherein V is a node set, E is a set of unordered binary groups formed by elements in V, and the starting nodes correspond to the valves:
V＝{e 1 ,e 2 ,e 3 ,e 4 ,e 5 ,e 6 ,e 7 ,e 8 ,e 9 ,e 10 ,e 11 ,e 12 ,e 13 ,e 14 e, where e 1 Representing the first node, i.e. the compressor group 01, e 2 Representing the second node, compressor string 02, and so on.
E＝{(e 1 ,e 4 ),(e 1 ,e 5 ),(e 1 ,e 6 ),(e 1 ,e 7 ),(e 2 ,e 4 ),(e 2 ,e 5 ),(e 2 ,e 6 ),(e 2 ,e 7 ),(e 3 ,e 4 ),(e 3 ,e 5 ),(e 3 ,e 6 ),(e 3 ,e 7 ),(e 4 ,e 8 ),(e 4 ,e 9 ),(e 4 ,e 10 ),(e 5 ,e 8 ),(e 5 ,e 9 ),(e 5 ,e 10 ),(e 6 ,e 8 ),(e 6 ,e 9 ),(e 6 ,e 10 ),(e 7 ,e 8 ),(e 7 ,e 9 ),(e 7 ,e 10 ),(e 8 ,e 11 ),(e 8 ,e 12 ),(e 8 ,e 13 ),(e 8 ,e 14 ),(e 9 ,e 11 ),(e 9 ,e 12 ),(e 9 ,e 13 ),(e 9 ,e 14 ),(e 10 ,e 11 ),(e 10 ,e 12 ),(e 10 ,e 13 ),(e 10 ,e 14 ) And (e), where 1 ,e 4 ) Representing a communication relationship between the first node and the fourth node, (e) 1 ,e 5 ) Representing the connectivity between the first node to the fifth node, and so on.
S10, acquiring the field operation state of a power system
The open and close states of the valve are collected in the system, and the valve is opened to be 1, and the valve is closed to be 0 as shown in the following table.
S20, acquiring a high-pressure air distribution network adjacency matrix A n×n
It will be appreciated that valve 1 corresponds to the connection path from node 1 to node 4, the valve 1 being in an open state, such that node 1 is contiguous with node 4, and in the contiguous matrix, the elements of the fourth column of the first row are 1; the valve 2 corresponds to a connection path from the node 1 to the node 5, and the valve 2 is in an off state, so that the node 1 is not adjacent to the node 5, and in the adjacent matrix, the element of the fifth column of the first row is 0; valve 3 corresponds to the connection path from node 1 to node 6, valve 3 is in the closed state, so that node 1 is not contiguous with node 6, the elements of the sixth column of the first row are 0 in the adjacency matrix, and so on, to obtain its adjacency matrix a 14×14 ，
S30, calculating a reachable matrix M from the adjacent matrix A:
s40, decomposing the reachable matrix:
s41, establishing a reachable set R (e i )
R(e 1 )＝{e 1 ,e 4 ,e 8 ,e 11 },
R(e 2 )＝{e 2 ,e 5 ,e 8 ,e 11 },
R(e 3 )＝{e 3 ,e 6 ,e 7 ,e 9 ,e 10 ,e 13 ,e 14 },
R(e 4 )＝{e 4 ,e 8 ,e 11 },
R(e 5 )＝{e 5 ,e 8 ,e 11 },
R(e 6 )＝{e 6 ,e 9 ,e 13 },
R(e 7 )＝{e 7 ,e 10 ,e 14 },
R(e 8 )＝{e 8 ,e 11 },
R(e 9 )＝{e 9 ,e 13 },
R(e 10 )＝{e 10 ,e 14 },
R(e 11 )＝{e 11 },
R(e 12 )＝{e 12 },
R(e 13 )＝{e 13 },
R(e 14 )＝{e 14 },
S42, establishing a look-ahead set Q (e) i )
Q(e 1 )＝{e 1 },
Q(e 2 )＝{e 2 },
Q(e 3 )＝{e 3 },
Q(e 4 )＝{e 1 ,e 4 },
Q(e 5 )＝{e 2 ,e 5 },
Q(e 6 )＝{e 3 ,e 6 },
Q(e 7 )＝{e 3 ,e 7 },
Q(e 8 )＝{e 1 ,e 2 ,e 4 ,e 5 ,e 8 },
Q(e 9 )＝{e 3 ,e 6 ,e 9 },
Q(e 10 )＝{e 3 ,e 7 ,e 10 },
Q(e 11 )＝{e 1 ,e 2 ,e 4 ,e 5 ,e 8 ,e 11 },
Q(e 12 )＝{e 12 },
Q(e 13 )＝{e 3 ,e 6 ,e 9 ,e 13 },
Q(e 14 )＝{e 3 ,e 7 ,e 10 ,e 14 },
S43, establishing a bottom layer node set B
R(e 1 )∩Q(e 1 )＝Q(e 1 )＝{e 1 },
R(e 2 )∩Q(e 2 )＝Q(e 2 )＝{e 2 },
R(e 3 )∩Q(e 3 )＝Q(e 3 )＝{e 3 },
The intersection of the remaining node's antecedent set and the reachable set is not equal to its antecedent set, therefore node set b= { e 1 ,e 2 ,e 3 The first node, the second node and the third node are bottom nodes of the system, and correspond to the compressor group 01, the compressor group 02 and the compressor group 03 respectively;
s44, carrying out region division on the distribution network diagram
R(e 1 )∩R(e 2 )≠φ，R(e 1 )∩R(e 3 )＝φ，R(e 2 )∩R(e 3 )＝φ，
Thus, the system is divided into two partitions,
π(V)＝{V 1 ,V 2 },
V 1 ＝{e 1 ,e 2 ,e 4 ,e 5 ,e 8 ,e 11 }，
V 2 ＝{e 3 ,e 6 ,e 7 ,e 9 ,e 10 ,e 13 ,e 14 }，
then node 1,2,4,5,8,11 is a connected branch and node 3,6,7,9,10,13,14 is a connected branch, the result of which is shown in fig. 3.
Further, in order to obtain a dynamic topology structure for guaranteeing a specific wind tunnel, the following steps are executed:
s45, establishing a top-level unit set T in each partition
For partition V 1 ，R(e 11 )∩Q(e 11 )＝R(e 11 )＝{e 11 Intersection of the antecedent set and the reachable set of the remaining nodes is not equal to its reachable set, thus, in partition V 1 In the scheme, the 11 th node, namely wind tunnel 01, is a top layer unit of the partition;
for partition V 2 ，R(e 13 )∩Q(e 13 )＝R(e 13 )＝{e 13 }，R(e 14 )∩Q(e 14 )＝R(e 14 )＝{e 14 Intersection of the antecedent set and the reachable set of the remaining nodes is not equal to its reachable set, thus, in partition V 2 In the scheme, the 13 th node, namely wind tunnel 03 and the 14 th node, namely wind tunnel 04, are top-level units of the partition;
s46, searching a wind tunnel w k Partition V being the top-level element h The collection of the compressor unit, the tank group, the valve and the gas transmission pipeline of the wind tunnel test is ensured, and the dynamic topological structure of the wind tunnel is particularly ensured.
For wind tunnel 01, it is partition V 1 The top layer element of the wind tunnel is ensured to be e 1 ,e 2 ,e 4 ,e 5 ,e 8 I.e. compressor 01, compressor 02, tank farm 01, tank farm 02, collection of main pipes 1;
for wind tunnel 03 and wind tunnel 04, the element is the top layer element of partition V2, and the elements of the two wind tunnels are ensured to be e 3 ,e 6 ,e 7 ,e 9 ,e 10 I.e. compressor 03, tank farm 04, main pipe 2 and main pipe 3.
Based on the method, the total volume of high-pressure air resources of a certain wind tunnel can be calculated, and specifically, the sum of the volumes of all elements of the wind tunnel is guaranteed:
for wind tunnel 01, the total volume of high-pressure air resources of wind tunnel 01 is ensured to be e of all elements of the wind tunnel 3 ,e 6 ,e 7 ,e 9 ,e 10 The sum of the volumes of (1), namely:
Vol 01 ＝V e1 +V e2 +V e4 +V e5 +V e8
for wind tunnel 03, as with wind tunnel 04,
Vol 03 ＝Vol 04 ＝V e3 +V e6 +V e7 +V e9 +V e10 。
it should be noted that, in this embodiment, the open state of the valve marked in the topology structure only represents that the pipeline is connected, and the air resource can flow through the pipeline to connect the elements on two sides, but it does not represent that the wind tunnels in the partition where the air resource is located are all consuming the air resource, and whether the wind tunnels in the partition where the air resource is located consume the air resource depends on whether the valve in the wind tunnel is opened or not. When the wind tunnel is to be tested, the valve in the wind tunnel is opened, and at the moment, the air resource in the partition where the valve is positioned can enter the wind tunnel to consume the air resource.
For example partition V 2 When the wind tunnel 03 is tested, the pipeline communicated with the wind tunnel 04 is filled with air, and the pipeline also belongs to the air resource consumption of the wind tunnel 03 for testing, so that the volume of elements on a path communicated with the wind tunnel 04 is required to be added when the total volume of high-pressure air resources of the wind tunnel 03 is calculated.
In the actual scheduling control process, only one wind tunnel can use power resources in one partition at the same time, namely, for V 2 Wind tunnels 03 and 04 cannot be tested simultaneously.
Example 3
The embodiment provides a system for executing the modeling method as in embodiment 1 or embodiment 2, which comprises hardware, an interface and software, wherein the hardware comprises a detector for monitoring test dynamics of a wind tunnel, tank group dynamics, compressor group dynamics and valve states; the interface adopts an OPCserver to upload the information acquired by the hardware to a central server, and the information is processed and visualized by software.
The software processing method is processed according to the method described in embodiment 1 or embodiment 2, and the implementation state of the software processing method can be displayed on a central server in a picture form, so that a worker can know the state of each element in the wind tunnel group conveniently.
The central server can also play a role of a dispatching center, and specifically, a local controller is arranged in each wind tunnel, the local controller sends a request signal for requesting test to the central server, the central server makes a test plan according to a dispatching rule, sends a signal for allowing test or not allowing test to the local controller, and the local controller also controls the opening or closing of a valve in the wind tunnel according to the feedback signal. When the central server allows the wind tunnel laboratory, the local controller controls the opening of the valve in the wind tunnel, the wind tunnel is communicated with each element in the partition, the high-pressure air resource starts to be supplied to the wind tunnel, and of course, when the test is finished, the local controller sends a test finishing signal to the central server, and after the central server receives the signal, the central server can send a test permitting signal to the rest wind tunnels in the partition requesting the test.
The scheduling rule is that in each partition, only one wind tunnel is guaranteed to use high-pressure air resources at the same time.
For example, in partition V 2 In the method, a wind tunnel 03 and a wind tunnel 04 both send a test request signal to a central server, when the wind tunnel 03 sends a request before the wind tunnel 04, the central server sends a test permission signal to the wind tunnel 03, a test disallowing signal to the wind tunnel 04, a local controller of the wind tunnel 03 opens a valve in the wind tunnel to start a test, when the test ends, a test ending signal is sent to the central server, after the central server receives the signal, the central server sends a test permission signal to the wind tunnel 04, and after receiving the signal, the local controller of the wind tunnel 04 controls the invention in the wind tunnel to start the test, and so on.
Example 4
A method for calculating the consumption of high-pressure air resources in a wind tunnel group, which is obtained by the method described in embodiment 1 or embodiment 2, guarantees a certain wind tunnel w k Firstly, calculating the sum of volumes of all elements of the wind tunnel, and then calculating the high-pressure air resource consumption P of the wind tunnel according to the tank group pressure difference in the partition before and after the test k ：
P k ＝ΔP×Vol k ×10。
Δp= -tank group pressure in the present partition at the start of the test-tank group pressure in the present partition at the end of the test
Therefore, the wind tunnel w for the current test can be calculated according to the pressure difference of the tank group in the current partition acquired at the beginning and ending time of the test k High pressure air resource consumption of (a).
The tank group pressure in this zone is the pressure of all tank groups in this zone, for example, for wind tunnel 03, the tank group pressure is the total pressure of tank groups 03 and 04.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather is intended to cover all modifications, equivalents, and alternatives falling within the spirit and principles of the invention.
Claims (7)
1. The wind tunnel group high-pressure air resource dynamic topological structure modeling method is characterized by comprising the following steps of:
s10, acquiring the field operation state of a power system
The power system comprises test dynamics of each wind tunnel of the wind tunnel group, tank group dynamics, compressor unit dynamics and valve states;
s20, acquiring a high-pressure air distribution network adjacency matrix A n×n
Taking a high-pressure air resource distribution network as a directed graph, taking a compressor unit, a pipeline airflow converging point, a tank group and a wind tunnel in the distribution network as nodes, and taking a valve in the distribution network as an edge to obtain an adjacent matrix:
wherein e i For the ith node, e j I epsilon [1, n for j-th node],j∈[1,n]The method comprises the steps of carrying out a first treatment on the surface of the n is the total node number in the high-pressure air resource distribution network;
s30, calculating a reachable matrix M
S40, decomposing the reachable matrix to obtain the communication branches of the gas distribution network diagram
The distribution network diagram is divided into a plurality of independent, direct or indirect influence free communication branches.
2. The method for modeling a dynamic topological structure of a high-pressure air resource of a wind tunnel group according to claim 1, wherein the step S40 comprises the steps of:
s41, establishing a reachable set R (e i )
Wherein, 
for any one node e i V is the set of nodes;
s42, establishing a look-ahead set Q (e) i )
S43, establishing a bottom node set B of the system
B＝{e i |e i E V and R (e) i )∩Q(e i )＝Q(e i )}
S44, carrying out region division on the distribution network diagram
For any node B, B' in the underlying node set B, if
R(b)∩R(b′)＝φ
The nodes b, b' belong to different areas, so that the system V can be divided into H partitions, which are recorded as
π(V)＝{V 1 ,V 2 ,…,V h ,…，V H }
Wherein each partition is a communication branch.
3. The method for modeling a dynamic topology of a wind tunnel group high-pressure air resource according to claim 2, further comprising, after step S44, the steps of:
s45, establishing a top-level unit set T in each partition
In the subarea
In the process,
T＝{e i |e i ∈V h and R (e) i )∩Q(e i )＝R(e i )}
S46, searching a wind tunnel w k Partition V being the top-level element h The method is characterized by comprising the steps of guaranteeing the collection of a compressor unit, a tank group, a valve and a gas pipeline of the wind tunnel test, wherein the combination of elements and connection relations of the elements in the collection is the dynamic topological structure of the wind tunnel.
4. A system for performing the method for modeling the dynamic topological structure of the high-pressure air resources of the wind tunnel group according to any one of claims 1 to 3, which is characterized by comprising hardware, interfaces and software,
the hardware comprises a detector for monitoring the test dynamic state of the wind tunnel, the tank group dynamic state, the compressor unit dynamic state and the valve state; the interface adopts an OPCserver to upload the information acquired by the hardware to a central server, and the information is processed and visualized by software.
5. The system of claim 4, wherein each wind tunnel includes a local controller, the local controller sends a request signal for "request test" to a central server, the central server formulates a test plan according to the scheduling rules, sends a signal for "allow test" or "not allow test" to the local controller, and the local controller further controls the opening or closing of valves in the wind tunnels according to the feedback signals.
6. The system of claim 5, wherein the scheduling rules are such that, within each partition, only one wind tunnel is guaranteed to use high pressure air resources at the same time.
7. The method for calculating the consumption of the high-pressure air resources of the wind tunnel group is characterized in that a certain wind tunnel w is guaranteed based on the method for modeling the dynamic topological structure of the high-pressure air resources of the wind tunnel group according to claim 3 k Is a dynamic topology of:
calculating and ensuring the wind tunnel w k Sum of volumes Vol of all elements of (C) k ，
This wind tunnel w k High-pressure air resource consumption P of (2) k The method comprises the following steps:
P k ＝ΔP×Vol k ×10，
Δp= -tank group pressure in the present partition at the start of the test-tank group pressure in the present partition at the end of the test.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210335439.4A CN114692524B (en) | 2022-03-31 | 2022-03-31 | Wind tunnel group high-pressure air resource dynamic topological structure modeling method, system and air consumption calculation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210335439.4A CN114692524B (en) | 2022-03-31 | 2022-03-31 | Wind tunnel group high-pressure air resource dynamic topological structure modeling method, system and air consumption calculation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114692524A CN114692524A (en) | 2022-07-01 |
| CN114692524B true CN114692524B (en) | 2023-05-05 |
Family
ID=82140492
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210335439.4A Active CN114692524B (en) | 2022-03-31 | 2022-03-31 | Wind tunnel group high-pressure air resource dynamic topological structure modeling method, system and air consumption calculation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114692524B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116108603B (en) * | 2023-04-12 | 2023-06-27 | 中国空气动力研究与发展中心计算空气动力研究所 | Method for constructing wind tunnel air supply valve unit level information physical system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111738627A (en) * | 2020-08-07 | 2020-10-02 | 中国空气动力研究与发展中心低速空气动力研究所 | Wind tunnel test scheduling method and system based on deep reinforcement learning |
| CN112069463A (en) * | 2020-08-28 | 2020-12-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Method for calculating high-pressure air resource consumption of wind tunnel group |
| CN112069446A (en) * | 2020-08-28 | 2020-12-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Wind tunnel group high-pressure air resource scheduling method and system |
| CN112069462A (en) * | 2020-08-28 | 2020-12-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Large icing wind tunnel high-pressure air resource consumption calculation method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8600572B2 (en) * | 2010-05-27 | 2013-12-03 | International Business Machines Corporation | Smarter-grid: method to forecast electric energy production and utilization subject to uncertain environmental variables |
-
2022
- 2022-03-31 CN CN202210335439.4A patent/CN114692524B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111738627A (en) * | 2020-08-07 | 2020-10-02 | 中国空气动力研究与发展中心低速空气动力研究所 | Wind tunnel test scheduling method and system based on deep reinforcement learning |
| CN112069463A (en) * | 2020-08-28 | 2020-12-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Method for calculating high-pressure air resource consumption of wind tunnel group |
| CN112069446A (en) * | 2020-08-28 | 2020-12-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Wind tunnel group high-pressure air resource scheduling method and system |
| CN112069462A (en) * | 2020-08-28 | 2020-12-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Large icing wind tunnel high-pressure air resource consumption calculation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114692524A (en) | 2022-07-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20230143654A1 (en) | Method for smart gas pipeline network inspection and internet of things system thereof | |
| CN114692524B (en) | Wind tunnel group high-pressure air resource dynamic topological structure modeling method, system and air consumption calculation | |
| CN106372747B (en) | Random forest-based reasonable line loss rate estimation method for transformer area | |
| CN108846089B (en) | Spatial data matching method for urban water supply professional pipelines and comprehensive pipelines | |
| Salomons et al. | Battle of water networks DMAs: Multistage design approach | |
| CN113568954B (en) | Parameter optimization method and system for preprocessing stage of network flow prediction data | |
| CN112069463B (en) | Method for calculating high-pressure air resource consumption of wind tunnel group | |
| CN114198644B (en) | Water supply network leakage detection control method based on DMA (direct memory access) monitoring related flow data | |
| Liu et al. | Multiphase DMA design methodology based on graph theory and many-objective optimization | |
| CN109766337A (en) | Storage method, electronic equipment, storage medium and the system of tree structure data | |
| CN109685684A (en) | A kind of low-voltage network topological structure method of calibration of knowledge based map | |
| CN104965846B (en) | Visual human's method for building up in MapReduce platform | |
| CN110287987B (en) | Method for discovering hierarchical organization structure of business system with hierarchical network structure | |
| Jaillet | Analysis of probabilistic combinatorial optimization problems in Euclidean spaces | |
| CN111368060B (en) | Self-learning method, device and system for conversation robot, electronic equipment and medium | |
| Mu et al. | Multistage iterative fully automatic partitioning in water distribution systems | |
| CN112990776A (en) | Distribution network equipment health degree evaluation method | |
| Li et al. | Optimizing combination of aircraft maintenance tasks by adaptive genetic algorithm based on cluster search | |
| CN113269449B (en) | Node level dispersion method for gas source gas supply flow scheduling of natural gas system | |
| CN110855480B (en) | Network fault fixed factor analysis method and device | |
| CN109872066B (en) | System complexity measurement method and device for sewage treatment plant | |
| CN111369052B (en) | Simplified road network KSP optimization algorithm | |
| CN108959638B (en) | Spatial data matching method for urban gas professional pipeline and comprehensive pipeline | |
| CN102215120B (en) | Data processing method and system for network element | |
| CN106209460B (en) | Power failure system restoration path optimization method based on network flow theory |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |