CN114692524A - Wind tunnel group high-pressure air resource dynamic topological structure modeling method and system and air consumption calculation - Google Patents

Wind tunnel group high-pressure air resource dynamic topological structure modeling method and system and air consumption calculation Download PDF

Info

Publication number
CN114692524A
CN114692524A CN202210335439.4A CN202210335439A CN114692524A CN 114692524 A CN114692524 A CN 114692524A CN 202210335439 A CN202210335439 A CN 202210335439A CN 114692524 A CN114692524 A CN 114692524A
Authority
CN
China
Prior art keywords
wind tunnel
pressure air
node
test
distribution network
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.)
Granted
Application number
CN202210335439.4A
Other languages
Chinese (zh)
Other versions
CN114692524B (en
Inventor
罗昌俊
王小飞
马永一
何福
付渲理
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202210335439.4A priority Critical patent/CN114692524B/en
Publication of CN114692524A publication Critical patent/CN114692524A/en
Application granted granted Critical
Publication of CN114692524B publication Critical patent/CN114692524B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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 and system and air consumption calculation. The topological structures of the tank groups and the pipelines are updated and recorded in real time, the element set for ensuring a specific wind tunnel and the sharing relation between the wind tunnel and other wind tunnels are obtained, and accurate data support is improved for formulating a test strategy and calculating wind tunnel resource consumption; the problems of inconvenience and low efficiency of manual control and coordination are avoided.

Description

Wind tunnel group high-pressure air resource dynamic topological structure modeling method and system and air consumption calculation
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
Large wind tunnel facilities are major strategic resources in China, wind tunnel groups usually comprise dozens or even dozens of production type wind tunnels, and each wind tunnel test needs multiple power resources such as pure water, electric power, high-pressure air, medium-pressure air, vacuum, nitrogen, hydrogen, oxygen and the like. In order 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 an air distribution network system. The air distribution network system is composed of a plurality of valves, pipelines and accessory equipment which are related to the transportation of high-pressure air from a compressor unit to a high-pressure tank group and to a wind tunnel. Because the wind tunnel groups share the high-pressure air resources in real time, different wind tunnel test requirements are guaranteed by switching and combining different tank groups, valves and pipelines.
In consideration of construction cost and high-pressure air resource requirements of different dynamic changes of a wind tunnel test, a high-pressure air source can be divided into different tank groups (partitions), and the different tank groups are formed by communicating a plurality of tank bodies through valves and pipelines; the system comprises a plurality of air transmission main pipes, a plurality of branch pipelines and a plurality of wind tunnels, wherein the air transmission main pipes are conveyed from a tank group to different wind tunnels; and a communicating valve is generally arranged between the air 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 gas supply and distribution relationship, which is composed of a high-pressure compressor set, valves, a tank group, a gas transmission main pipe, branch pipelines and a wind tunnel implementing guarantee, is called a dynamic topological 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 power resource supply configuration (pressure grade, standard square capacity and the like), as shown in fig. 1 and fig. 2.
The high-frequency and multi-type pneumatic tests carried out by the wind tunnel groups cause the complexity of high-pressure air resource supply guarantee, the conditions that a plurality of wind tunnels share the same tank group in a time-sharing manner, a plurality of wind tunnels use different tank groups in parallel, a large wind tunnel uses a plurality of tank groups simultaneously and the like exist, and great challenges are provided for the safe operation of the wind tunnel tests, the timely guarantee of power resources, the reasonable allocation and the accurate measurement of consumption.
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, the resonance phenomenon is easy to occur, the safety of the whole high-voltage 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 the high-pressure air resource conflict, the current wind tunnel test operation mostly adopts a manual coordination mode and is communicated through a 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, to realize consumption metering of a single wind tunnel and a single test, two schemes of additionally installing a flow meter on a wind tunnel pipeline or metering differential pressure of a tank group before and after the test can be adopted. However, because the pressure of the tank group is changed rapidly during the wind tunnel group test, the scheme of installing a flow meter on the wind tunnel pipeline is inaccurate in monitoring, the cost is high, and the wind tunnel test flow field can be influenced; by adopting the tank group differential pressure method metering scheme, the topological structures of the tank group and the pipelines need to be mastered in real time so as to obtain the volumes (standard square volumes) of the air supply tank group and the pipelines of the wind tunnel under the test to calculate the high-pressure air consumption of the wind tunnel test.
Therefore, the management of the high-pressure air resource distribution network topological 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 aiming at the topological structure management of the high-pressure air resource distribution network of the complex and dynamic wind tunnel group is found at present.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a modeling method of a dynamic topological structure of a wind tunnel group high-pressure air resource distribution network based on a graphic database and an adjacency matrix, establishes a model for automatically identifying, storing and updating the dynamic topological structure of the high-pressure air resource distribution network, and is convenient for the subsequent calculation of air resource consumption.
A wind tunnel group high-pressure air resource dynamic topological structure modeling method is characterized by comprising the following steps:
s10, acquiring the field running state of the power system
The power system comprises the test dynamics of each wind tunnel of the wind tunnel group, the tank group dynamics, the compressor unit dynamics and the valve state;
s20, obtaining a high-pressure air distribution network adjacent matrix An×n
Taking a high-pressure air resource distribution network as a directed graph, taking a compressor set, a pipeline airflow confluence 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 adjacency matrix:
Figure BDA0003576610860000031
Figure BDA0003576610860000032
wherein ,eiIs the ith node, ejFor the j-th node, i ∈ [1, n ]],j∈[1,n]；
S30, calculating a reachable matrix M
Figure BDA0003576610860000033
Figure BDA0003576610860000034
S40, decomposing the reachable matrix to obtain a communication branch of the gas distribution network diagram
The distribution network diagram is divided into a plurality of independent communication branches without direct or indirect influence.
Further, step S40 includes the steps of:
s41, establishing reachable set R (e) of nodesi)
Figure BDA0003576610860000041
wherein ,
Figure BDA0003576610860000042
for any one node eiV is a set of nodes;
s42, establishing a first-class set Q (e) of nodesi)
Figure BDA0003576610860000043
S43, establishing a bottom node set B of the system
B＝{ei|eiE.g. V and R (e)i)∩Q(ei)＝Q(ei)}
S44, carrying out region division on the gas distribution network diagram
For any node B, B' in the bottom node set B, if
R(b)∩R(b′)＝φ
The nodes b, b' belong to different zones, so that the system V can be divided into H partitions, denoted as
π(V)＝{V1,V2,…,Vh,…，VH}
Wherein each partition is a connected branch.
Further, after step S44, the method further includes the following steps:
s45, establishing a top-level unit set T in each partition
In a partition
Figure BDA0003576610860000044
In (1),
T＝{ei|ei∈Vhand R (e)i)∩Q(ei)＝R(ei)}
S46, finding a certain wind tunnel wkPartition V being a top level elementhAnd obtaining a set of a compressor unit, a tank group, a valve and a gas transmission pipeline which guarantee the wind tunnel test, wherein the combination of each element in the set and the connection relation of each element is the dynamic topological structure which guarantees the wind tunnel.
The system for implementing 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 detectors for monitoring the test dynamics, tank group dynamics, compressor unit dynamics and valve states of the wind tunnel; the interface adopts an OPCserver, uploads the information acquired by hardware to a central server, and is processed and visualized by software.
Furthermore, each wind tunnel comprises a local controller, the local controller sends a request signal of 'request test' to the central server, the central server makes a test plan according to a scheduling rule, and sends a 'test permission' or 'test non-permission' signal to the local controller, and the local controller controls the opening or closing of the valve in the wind tunnel according to the feedback signal.
Further, the scheduling rule is that in each partition, only one wind tunnel is ensured to use high-pressure air resources at the same time.
The invention also provides a wind tunnel group high-pressure air resource consumption calculation method which is characterized by ensuring a certain wind tunnel w based on the wind tunnel group high-pressure air resource dynamic topological structure modeling methodkDynamic topology of (2):
calculation ensures the wind tunnel wkVol of the sum of the volumes of all elementsw，
This wind tunnel wkHigh pressure air resource consumption PkComprises the following steps:
Pktank group pressure at the beginning of wind tunnel test-at the end of wind tunnel testPressure of tank group | Volw×10。
Compared with the prior art, the wind tunnel group high-pressure air resource dynamic topological structure modeling method, the system and the air consumption calculation have the following beneficial effects:
(1) according to the modeling method, pipelines and tank groups in use of a specific wind tunnel can be automatically identified, and can be updated in real time;
(2) the method can update and record the topological structures of the tank groups and the pipelines in real time, obtain an element set for ensuring a specific wind tunnel and the sharing relation between the wind tunnel and other wind tunnels, and improve accurate data support for formulating a test strategy and calculating wind tunnel resource consumption; the problems of inconvenience and low efficiency of manual control and coordination are avoided;
(3) by utilizing the modeling method, the high-pressure air resource consumption of the 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 drawings needed to be used in the embodiments of the present invention or in the description of the prior art will be briefly described below, and it is obvious that the drawings described below are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a schematic diagram I 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 wind tunnel group high-pressure air resource dynamic topological structure modeling method according to an embodiment of the present invention;
fig. 4 is an example of a wind tunnel group high-pressure air resource dynamic topological structure modeling method 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 particular examples set forth below are illustrative only and are not intended to be limiting.
Example 1
A wind tunnel group high-pressure air resource dynamic topological structure modeling method is shown in FIG. 3 and comprises the following steps:
s10, acquiring the field running state of the power system
The power system comprises the test dynamics of each wind tunnel of the wind tunnel group, the tank group dynamics, the compressor unit dynamics and the valve state;
s20, obtaining a high-pressure air distribution network adjacent matrix An×n
Taking a high-pressure air resource distribution network as a directed graph, taking a compressor set, a pipeline airflow confluence 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 adjacency matrix:
Figure BDA0003576610860000061
Figure BDA0003576610860000071
wherein ,eiIs the ith node, ejFor the j-th node, i ∈ [1, n ]],j∈[1,n](ii) a Adjacent to the element a of the matrix A ij1 represents the connection relationship between the node i and the node j, 0 is not adjacent, and 1 is adjacent. And the kth row and the kth column have the same physical meaning and all represent the node k, so that the aim of finding all nodes adjacent to the node k is to mainly find whether an element with 1 exists in the kth row or the kth column.
In the invention, the valve is used as the edge because the pipeline determines the connection relation of the distribution network, the valve equipment controls the switching of the tank group, the unit and the pipeline, the switching of the valve can cause the change of a physical connection model and the adjacent relation of valve nodes, the state of the valve is sampled in real time, and the corresponding adjacent state can be obtained, thereby obtaining an adjacent matrix;
further, in order to describe the adjacency relationship between nodes of the distribution network, the present embodiment establishes a graph database including three relationship tables:
1) establishing a valve-backflow node table: the table records the actual position relationship between the valve and the pipeline confluence node, and the structure of the table is shown in table 1 (the start node and the end node in the table can be interchanged).
TABLE 1 valve-confluence node table
Valve number Starting node Terminal 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 structure of the table is shown in table 2.
Table 2 node representation table
Figure BDA0003576610860000072
Figure BDA0003576610860000081
3) Establishing a valve state table: the states of all valves in all high-pressure distribution networks are obtained through the wind tunnel and the high-pressure air system operation monitoring system and are recorded by a meter. The table records the open and close states of all valves at a certain time, and the structure of the table 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 Val2 1
…… …… ……
In the graph database, the condition query 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 structures of the tables are shown in table 4.
TABLE 4 adjacency matrix table
Status value Identification number of starting node Identification number of termination point
0 1 2
1 3 4
…… …… ……
S30, calculating the reachable matrix M to obtain the sum of direct and indirect relations of each node
The reachable matrix refers to a matrix form for describing the reachable degree of each node of the graph after a path with a certain length is passed, and is also an n × n square matrix represented by M:
Figure BDA0003576610860000082
Figure BDA0003576610860000091
up to the ith row and jth column element M of matrix M ij1 means that node i and node j have a direct (with the pipeline directly connected, i points to j) or indirect (at least through one intermediate node) connection, otherwise it is equal to 0. For example, if i points to g and g points to j, i and j have 1-step indirect relationship; in the same way, there are 2 steps, 3 steps, andijotherwise, it is 0.
S40, decomposing the reachable matrix to obtain a communication branch of the gas distribution network diagram
The distribution network diagram is divided into a plurality of independent communication branches without direct or indirect influence.
The method specifically comprises the following steps:
s41, establishing reachable set R (e) of nodesi)
eiHaving a direct or indirect influence on a node, these nodes form eiThe reachable set of (c):
Figure BDA0003576610860000092
wherein ,
Figure BDA0003576610860000093
for any one node eiV is a set of nodes;
s42, establishing a first-class set Q (e) of nodesi)
Any one node pair eiWith direct or indirect influence, these nodes form eiThe first set of (2):
Figure BDA0003576610860000094
s43, establishing a bottom node set B
The condition that there is no set of nodes to which any node points but not to which it points, constitutes the set of bottom nodes B of the system V:
B＝{ei|eiis e.g. V andR(ei)∩Q(ei)＝Q(ei)}
s44, carrying out region division on the gas distribution network diagram
For any node B, B' in the bottom node set B, if
R(b)∩R(b′)＝φ
Then nodes b, b' cannot point to the same node but belong to different regions, so system V can be divided into H partitions, denoted as
π(V)＝{V1,V2,…,Vh,…，VH}
Wherein each partition is a connected branch.
Further, the step S45 and the step S46 are executed to obtain a dynamic topology structure that ensures a specific wind tunnel:
s45, establishing a top-level unit set T in each partitionh
In a partition
Figure BDA0003576610860000101
In the method, the condition that a set of nodes which are only influenced by the node and do not influence the node does not exist, and a top-level unit set T is formedh：
Th＝{ei|ei∈VhAnd R (e)i)∩Q(ei)＝R(ei)}
S46, finding a certain wind tunnel wkPartition V being a top level elementhThe compressor unit, the tank group, the valve and the gas transmission pipeline which guarantee the wind tunnel test are integrated, and the dynamic topological structure which guarantees the wind tunnel is obtained.
Example 2
Fig. 2 is a schematic diagram of a physical connection configuration of a high pressure air resource, the system having three compressor banks: compressor unit 01, compressor unit 02 and compressor unit 03, four jar groups: tank group 01, tank group 02, tank group 03 and tank group 04, the three person in charge: main tube 1, main tube 2 and main tube 3, and four wind tunnels: it should be noted that the actual high-pressure air resource system is more complex, and fig. 2 is only an example of simplifying the system, where valves are simplified, and each node has a valve (i.e., on a side line). During testing, corresponding valves are opened according to test requirements, and a wind tunnel group high-pressure air resource dynamic topological structure modeling method is constructed through 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 sequence of the compressor unit, the tank group, the main pipe and the wind tunnel, namely, the compressor unit 01, the compressor unit 02, the compressor unit 03, the tank group 01, the tank group 02, the jet, the wind tunnel 03 and the wind tunnel 04 are set as numbers 1,2, 3, 4,5, the jet, 13 and 14 respectively, as shown in the following table:
Figure BDA0003576610860000102
Figure BDA0003576610860000111
s02, establishing a valve-node table:
valve number Starting node Terminal node Valve number Starting node Terminal node
Val1
1 4 Val19 6 8
Val2 1 5 Val20 6 9
Val3 1 6 Val21 6 10
Val4 1 7 Val22 7 8
Val5 2 4 Val23 7 9
Val6 2 5 Val24 7 10
Val7 2 6 Val25 8 11
Val8 2 7 Val26 8 12
Val9 3 4 Val27 8 13
Val10 3 5 Val28 8 14
Val11 3 6 Val29 9 11
Val12 3 7 Val30 9 12
Val13 4 8 Val31 9 13
Val14 4 9 Val32 9 14
Val15 4 10 Val33 10 11
Val16 5 8 Val34 10 12
Val17 5 9 Val35 10 13
Val18 5 10 Val36 10 14
Therefore, the physical connection relation of the wind tunnel group high-pressure air resource system is obtained, and 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 disordered binary groups formed by elements in V, and corresponds to the starting nodes of all valves:
V＝{e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14}, wherein ,e1Representing the first node, i.e. compressor group 01, e2Representing the second node, compressor bank 02, and so on.
E＝{(e1,e4),(e1,e5),(e1,e6),(e1,e7),(e2,e4),(e2,e5),(e2,e6),(e2,e7),(e3,e4),(e3,e5),(e3,e6),(e3,e7),(e4,e8),(e4,e9),(e4,e10),(e5,e8),(e5,e9),(e5,e10),(e6,e8),(e6,e9),(e6,e10),(e7,e8),(e7,e9),(e7,e10),(e8,e11),(e8,e12),(e8,e13),(e8,e14),(e9,e11),(e9,e12),(e9,e13),(e9,e14),(e10,e11),(e10,e12),(e10,e13),(e10,e14) Wherein (e)1,e4) Showing the connection relationship from the first node to the fourth node, (e)1,e5) Representing the connectivity from the first node to the fifth node, and so on.
S10, acquiring the field running state of the power system
The open and close states of the valve are collected in the system, as shown in the following table, wherein the open state of the valve is 1, and the close state of the valve is 0.
Figure BDA0003576610860000121
Figure BDA0003576610860000131
S20, obtaining a high-pressure air distribution network adjacent matrix An×n
It will be understood that valve 1 corresponds to the connection path from node 1 to node 4, and that valve 1 is in the open state, so that node 1 is adjacent to node 4, and in the adjacent matrix, the element in the first row and the fourth column is 1; valve 2 corresponds to the connection path from node 1 to node 5, and valve 2 is in the off state, so node 1 is not adjacent to node 5, and in the adjacent matrix, the element in the fifth row and column is 0; the valve 3 corresponds to the connection path from the node 1 to the node 6, and the valve 3 is in an off state, so that the node 1 is not adjacent to the node 6, and in the adjacent matrix, the element in the first row and the sixth column is 014×14，
Figure BDA0003576610860000132
S30, calculating a reachable matrix M by the adjacency matrix A:
Figure BDA0003576610860000141
s40, decomposing the reachable matrix:
s41, establishing an accessible set R (e) of nodesi)
R(e1)＝{e1,e4,e8,e11},
R(e2)＝{e2,e5,e8,e11},
R(e3)＝{e3,e6,e7,e9,e10,e13,e14},
R(e4)＝{e4,e8,e11},
R(e5)＝{e5,e8,e11},
R(e6)＝{e6,e9,e13},
R(e7)＝{e7,e10,e14},
R(e8)＝{e8,e11},
R(e9)＝{e9,e13},
R(e10)＝{e10,e14},
R(e11)＝{e11},
R(e12)＝{e12},
R(e13)＝{e13},
R(e14)＝{e14},
S42, establishing a first-class set Q (e) of nodesi)
Q(e1)＝{e1},
Q(e2)＝{e2},
Q(e3)＝{e3},
Q(e4)＝{e1,e4},
Q(e5)＝{e2,e5},
Q(e6)＝{e3,e6},
Q(e7)＝{e3,e7},
Q(e8)＝{e1,e2,e4,e5,e8},
Q(e9)＝{e3,e6,e9},
Q(e10)＝{e3,e7,e10},
Q(e11)＝{e1,e2,e4,e5,e8,e11},
Q(e12)＝{e12},
Q(e13)＝{e3,e6,e9,e13},
Q(e14)＝{e3,e7,e10,e14},
S43, establishing a bottom node set B
R(e1)∩Q(e1)＝Q(e1)＝{e1},
R(e2)∩Q(e2)＝Q(e2)＝{e2},
R(e3)∩Q(e3)＝Q(e3)＝{e3},
The intersection of the antecedent and reachable sets of the remaining nodes is not equal to its antecedent set, and thus node set B ═ e1,e2,e3That is, the first node, the second node and the third node are bottom nodes of the system, and respectively correspond to the compressor group 01, the compressor group 02 and the compressor group 03;
s44, carrying out region division on the gas distribution network diagram
R(e1)∩R(e2)≠φ，R(e1)∩R(e3)＝φ，R(e2)∩R(e3)＝φ，
Thus, the system is divided into two partitions,
π(V)＝{V1,V2},
V1＝{e1,e2,e4,e5,e8,e11}，
V2＝{e3,e6,e7,e9,e10,e13,e14}，
then, nodes 1,2,4,5,8,11 are a connected branch, and nodes 3,6,7,9,10,13,14 are a connected branch, and the result is shown in fig. 3.
Further, in order to obtain a dynamic topological 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 V1，R(e11)∩Q(e11)＝R(e11)＝{e11The intersection of the antecedent and reachable sets of the remaining nodes is not equal to its reachable set, so in partition V1In 11, the wind tunnel 01 as the node is the top unit of the subarea;
for partition V2，R(e13)∩Q(e13)＝R(e13)＝{e13}，R(e14)∩Q(e14)＝R(e14)＝{e14}, the intersection of the antecedent and reachable sets of the remaining nodes does not equal its reachable set, so in partition V2In the partition, the 13 th node, namely the wind tunnel 03, and the 14 th node, namely the wind tunnel 04, are top-level units of the partition;
s46, finding a certain wind tunnel wkPartition V being a top level elementhThe compressor unit, the tank group, the valve and the gas transmission pipeline which guarantee the wind tunnel test are integrated, and the dynamic topological structure which guarantees the wind tunnel is obtained.
For wind tunnel 01, it is the sector V1The element ensuring the wind tunnel is e1,e2,e4,e5,e8I.e. compressors01, compressor 02, tank group 01, tank group 02, set of main pipes 1;
for the wind tunnel 03 and the wind tunnel 04, the top element of the partition V2 is ensured to be e3,e6,e7,e9,e10Namely, compressor 03, tank group 04, main pipe 2, and main pipe 3.
Based on this, the total volume of the high-pressure air resource of a certain wind tunnel can be calculated, specifically, the sum of the volumes of all elements of the wind tunnel is guaranteed:
for the wind tunnel 01, the total volume of the high-pressure air resource of the wind tunnel 01 is ensured to be e ensuring all elements of the wind tunnel3,e6,e7,e9,e10The sum of the volumes of (a) and (b), i.e.:
Vol01＝Ve1+Ve2+Ve4+Ve5+Ve8
as for the wind tunnel 03, like the wind tunnel 04,
Vol03＝Vol04＝Ve3+Ve6+Ve7+Ve9+Ve10。
it should be noted that, in this embodiment, the state of the valve marked in the topological structure is open, which only represents that the pipeline is communicated, and the air resource can flow through the pipeline to communicate the elements on the two sides, but does not represent that the wind tunnels in the partition where the pipeline is located all consume the air resource, and whether the wind tunnels in the partition where the pipeline is located consume the air resource is determined by whether the valve inside the wind tunnel is open. When the wind tunnel needs to be tested, a valve in the wind tunnel is opened, and at the moment, the air resource in the partition where the wind tunnel is located can enter the wind tunnel to consume the air resource.
E.g. partition V2In the case of the wind tunnel 03 test, the two wind tunnels 03 and 04 share a communication branch, and although the wind tunnel 04 does not consume air resources, the pipeline communicating with the wind tunnel 04 is filled with air, which also belongs to the air resource consumption of the wind tunnel 03 to perform the test, so that the volume of an element on a path communicating with the wind tunnel 04 needs to be added when calculating the total volume of the high-pressure air resources of the wind tunnel 03.
In the actual scheduling control process, it is required to ensure that only one wind tunnel can use power resources in one subarea, namely V2The wind tunnels 03 and 04 cannot be tested simultaneously in this sector.
Example 3
The embodiment provides a system for executing the modeling method as in embodiment 1 or embodiment 2, which includes hardware, an interface, and software, wherein the hardware includes detectors for monitoring test dynamics, tank group dynamics, compressor group dynamics, and valve states of a wind tunnel; the interface adopts an OPCserver, uploads the information acquired by hardware to a central server, and 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 the central server in the form of a picture, so that workers 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, specifically, a local controller is arranged in each wind tunnel, the local controller sends a request signal of 'request test' to the central server, the central server formulates a test plan according to a dispatching rule, sends a 'test permission' or 'test non-permission' signal 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. The method comprises the steps that a central server allows a wind tunnel laboratory, a local controller controls a valve in a wind tunnel to be opened, the wind tunnel is communicated with elements in a subarea, high-pressure air resources begin to be supplied to the wind tunnel, certainly, when a test is finished, the local controller sends a test finishing signal to the central server, and after receiving the signal, the central server can send a test allowing signal to the rest wind tunnels requesting the test in the subarea.
The dispatching rule is that in each subarea, only one wind tunnel is ensured to use high-pressure air resources.
For example, in partition V2In the wind tunnel 03 and the wind tunnel 04 both send 'request test' signals to the central server, and when the wind tunnels03 sends a request before the wind tunnel 04, the central server sends a test permission signal to the wind tunnel 03 and sends a test non-permission signal to the wind tunnel 04, the local controller of the wind tunnel 03 opens a valve in the wind tunnel to start the test, when the test is finished, the central server sends a test finish signal, after receiving the signal, the central server sends a test permission signal to the wind tunnel 04, after receiving the signal, the local controller of the wind tunnel 04 controls the invention in the wind tunnel to open to start the test, and so on.
Example 4
A method for calculating the resource consumption of high-pressure air in a wind tunnel group, which guarantees a certain wind tunnel w obtained according to the method described in embodiment 1 or embodiment 2kThe dynamic topological structure of (1) firstly calculates the volume sum of all elements ensuring the wind tunnel, and then calculates the high-pressure air resource consumption P of the wind tunnel in the test according to the tank group pressure difference in the partition before and after the testk：
Pk＝ΔP×Volk×10。
Pressure of tank group in the partition at the start of the test-pressure of tank group in the partition at the end of the test
Therefore, the wind tunnel w of the test can be calculated according to the pressure difference of the tank group in the partition acquired by the time of starting and ending the testkThe high-pressure air resource consumption.
It should be noted that the tank group pressure in the present division refers to all the tank group pressures in the present division, and for example, in the case of the wind tunnel 03, the tank group pressure is the total pressure of the tank group 03 and the tank group 04.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.

Claims (7)

1. A wind tunnel group high-pressure air resource dynamic topological structure modeling method is characterized by comprising the following steps:
s10, acquiring the field running state of the power system
The power system comprises the test dynamics of each wind tunnel of the wind tunnel group, the tank group dynamics, the compressor unit dynamics and the valve state;
s20, obtaining a high-pressure air distribution network adjacent matrix An×n
Taking a high-pressure air resource distribution network as a directed graph, taking a compressor set, a pipeline airflow confluence 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 adjacency matrix:
Figure FDA0003576610850000011
Figure FDA0003576610850000012
wherein ,eiIs the ith node, ejFor the j-th node, i ∈ [1, n ]],j∈[1,n](ii) a n is the total node number in the high-pressure air resource distribution network;
s30, calculating a reachable matrix M
Figure FDA0003576610850000013
Figure FDA0003576610850000014
S40, decomposing the reachable matrix to obtain a communication branch of the gas distribution network diagram
The distribution network diagram is divided into a plurality of independent communication branches without direct or indirect influence.
2. The modeling method for the dynamic topological structure of the wind tunnel swarm high-pressure air resource according to claim 1, wherein the step S40 comprises the following steps:
s41. buildingReachable set of vertical nodes R (e)i)
Figure FDA0003576610850000021
R(ei)＝{ej|ej∈V,mij＝1}
wherein ,
Figure FDA0003576610850000022
for any one node eiV is a set of nodes;
s42, establishing a first-class set Q (e) of nodesi)
Figure FDA0003576610850000023
Q(ei)＝{ej|ej∈V,mji＝1}
S43, establishing a bottom node set B of the system
B＝{ei|eiE.g. V and R (e)i)∩Q(ei)＝Q(ei)}
S44, carrying out region division on the gas distribution network diagram
For any node B, B' in the bottom node set B, if
R(b)∩R(b′)＝φ
The nodes b, b' belong to different zones, so that the system V can be divided into H partitions, denoted as
π(V)＝{V1,V2,…,Vh,…，VH}
Wherein each partition is a connected branch.
3. The wind tunnel crowd high pressure air resource dynamic topological structure modeling method according to claim 2, further comprising the following steps after step S44:
s45, establishing a top-level unit set T in each partition
In a partition
Figure FDA0003576610850000024
In (1),
T＝{ei|ei∈Vhand R (e)i)∩Q(ei)＝R(ei)}
S46, finding a certain wind tunnel wkPartition V being a top level elementhAnd obtaining a set of a compressor unit, a tank group, a valve and a gas transmission pipeline which guarantee the wind tunnel test, wherein the combination of each element in the set and the connection relation of each element is the dynamic topological structure which guarantees the wind tunnel.
4. A system for executing the wind tunnel swarm high-pressure air resource dynamic topological structure modeling method according to any one of claims 1-3, characterized by comprising hardware, interfaces and software,
the hardware comprises detectors for monitoring the test dynamics, tank group dynamics, compressor unit dynamics and valve states of the wind tunnel; the interface adopts an OPCserver, uploads the information acquired by hardware to a central server, and is processed and visualized by software.
5. The system of claim 4, wherein each wind tunnel comprises a local controller, the local controller sends a request signal of "request test" to the central server, the central server makes a test plan according to the scheduling rule, and sends a "test permission" or "test non-permission" signal to the local controller, and the local controller further controls the opening or closing of the valve in the wind tunnel according to the feedback signal.
6. The system of claim 5, wherein the scheduling rules ensure that, within each zone, only one wind tunnel uses the high pressure air resource at a time.
7. A wind tunnel group high-pressure air resource consumption calculation method is characterized in that a certain guarantee is obtained based on the wind tunnel group high-pressure air resource dynamic topological structure modeling method of claim 3Wind tunnel wkDynamic topology of (2):
calculation ensures the wind tunnel wkVolume sum Vol of all elements of (2)k，
This wind tunnel wkHigh pressure air resource consumption PkComprises the following steps:
Pk＝ΔP×Volk×10。
Δ P ═ tank group pressure in the current zone at the start of the test-tank group pressure in the current zone at the end of the test.
CN202210335439.4A 2022-03-31 2022-03-31 Wind tunnel group high-pressure air resource dynamic topological structure modeling method, system and air consumption calculation Active CN114692524B (en)

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 true CN114692524A (en) 2022-07-01
CN114692524B 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)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116108603A (en) * 2023-04-12 2023-05-12 中国空气动力研究与发展中心计算空气动力研究所 Method for constructing wind tunnel air supply valve unit level information physical system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110307109A1 (en) * 2010-05-27 2011-12-15 Sri-Jayantha Sri M Smarter-Grid: Method to Forecast Electric Energy Production and Utilization Subject to Uncertain Environmental Variables
CN111738627A (en) * 2020-08-07 2020-10-02 中国空气动力研究与发展中心低速空气动力研究所 Wind tunnel test scheduling method and system based on deep reinforcement learning
CN112069446A (en) * 2020-08-28 2020-12-11 中国空气动力研究与发展中心计算空气动力研究所 Wind tunnel group high-pressure air resource scheduling method and system
CN112069463A (en) * 2020-08-28 2020-12-11 中国空气动力研究与发展中心计算空气动力研究所 Method for calculating high-pressure air resource consumption of wind tunnel group
CN112069462A (en) * 2020-08-28 2020-12-11 中国空气动力研究与发展中心计算空气动力研究所 Large icing wind tunnel high-pressure air resource consumption calculation method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110307109A1 (en) * 2010-05-27 2011-12-15 Sri-Jayantha Sri M Smarter-Grid: Method to Forecast Electric Energy Production and Utilization Subject to Uncertain Environmental Variables
CN111738627A (en) * 2020-08-07 2020-10-02 中国空气动力研究与发展中心低速空气动力研究所 Wind tunnel test scheduling method and system based on deep reinforcement learning
CN112069446A (en) * 2020-08-28 2020-12-11 中国空气动力研究与发展中心计算空气动力研究所 Wind tunnel group high-pressure air resource scheduling method and system
CN112069463A (en) * 2020-08-28 2020-12-11 中国空气动力研究与发展中心计算空气动力研究所 Method for calculating high-pressure air resource consumption of wind tunnel group
CN112069462A (en) * 2020-08-28 2020-12-11 中国空气动力研究与发展中心计算空气动力研究所 Large icing wind tunnel high-pressure air resource consumption calculation method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116108603A (en) * 2023-04-12 2023-05-12 中国空气动力研究与发展中心计算空气动力研究所 Method for constructing wind tunnel air supply valve unit level information physical system
CN116108603B (en) * 2023-04-12 2023-06-27 中国空气动力研究与发展中心计算空气动力研究所 Method for constructing wind tunnel air supply valve unit level information physical system

Also Published As

Publication number Publication date
CN114692524B (en) 2023-05-05

Similar Documents

Publication Publication Date Title
WO2022088890A1 (en) Checking method and system for medium-and-long-term maintenance plan in power grid, and device and storage medium
Qing et al. Decentralized unscented Kalman filter based on a consensus algorithm for multi-area dynamic state estimation in power systems
CN110134833A (en) Diagram data modeling and method towards power network topology management
CN106202335B (en) A kind of traffic big data cleaning method based on cloud computing framework
CN108020750A (en) A kind of power distribution network power failure range analysis method and analytical equipment
CN107545085B (en) Automatic generation method for process level network connection diagram of intelligent substation
CN108228724A (en) Power grid GIS topology analyzing method and storage medium based on chart database
CN109873501A (en) A kind of low-voltage network topology automatic identification method
CN114692524A (en) Wind tunnel group high-pressure air resource dynamic topological structure modeling method and system and air consumption calculation
CN112685570B (en) Multi-label graph-based power grid network frame knowledge graph construction method
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
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
CN107729939A (en) A kind of CIM extended method and device towards newly-increased power network resources
Jaillet Analysis of probabilistic combinatorial optimization problems in Euclidean spaces
CN111709102B (en) Water supply network partitioning method based on hierarchical clustering
CN104901328A (en) Multi-terminal flexible DC control mode automatic identification method based on complex control network
CN106209460B (en) Power failure system restoration path optimization method based on network flow theory
CN110855480B (en) Network fault fixed factor analysis method and device
CN111369052B (en) Simplified road network KSP optimization algorithm
CN109872066B (en) System complexity measurement method and device for sewage treatment plant
Changjun et al. A Dynamic Topology Modeling Method in Wind Tunnel Group Based on Adjacency Matrix
Li et al. Advanced split BIRCH algorithm in reconfigurable network
CN109495305A (en) A kind of polymorphic flow network reliability estimation method of the more commodity of mostly distribution and device

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