A pipeline gas column transient flow simulation method, device, electronic equipment and medium

Abstract

The present application provides a method, device, electronic equipment and medium for simulating transient flow of gas column in pipeline, which relates to the field of hydraulic numerical simulation calculation. The method includes: constructing elastic water control algorithm, gas column control algorithm and gas-liquid boundary control algorithm; establishing interface tracking algorithm at the gas-liquid boundary according to the transient state of gas-liquid boundary and rigid interface assumption; using moving wave to calculate the target pressure and target flow of any node to be measured of the target pipeline; the present application can combine the characteristics and mechanism of transient flow of impact gas column at the end of the pipeline nozzle, enrich and improve the theory of transient flow of impact gas column, and provide scientific judgment for the engineering design and safe operation of water pipeline system. The present application can also solve the problem of poor simulation effect of existing models and actual engineering, so as to simulate the actual pipeline in real time and accurately, grasp the hydraulic safety problem of the engineering, and timely check the hidden dangers, which is of great significance to the safety and stability of the actual pipeline.

Description

Technical Field

[0001] The present application relates to the field of hydraulic numerical simulation calculation, and in particular, to a method, device, electronic equipment and medium for simulating transient flow of a pipeline gas column. Background Art

[0002] Water pipelines are an indispensable part of hydropower stations, long-distance water diversion projects, urban pipe networks and other projects. At present, the design of water pipelines pays little attention to air columns, and the hazards of air columns are often ignored. When air columns remain in water pipelines, they are destructive to the pipelines and can easily cause engineering accidents.

[0003] The existing transient flow model of the impact air column considering the air nozzle does not comprehensively consider whether the water body is compressible, whether the air column is compressible, and the dynamic capture of the gas-liquid boundary. As a result, the boundary conditions and initial conditions used in the numerical simulation are limited, and the simulation effect on more complex practical problems is poor. Summary of the invention

[0004] The purpose of the embodiments of the present application is to provide a pipeline gas column transient flow simulation method, device, electronic equipment and medium to solve the above-mentioned problems existing in the prior art and accurately obtain the pressure and flow of the node to be measured in the pipeline.

[0005] In a first aspect, a pipeline gas column transient flow simulation method is provided, which may include:

[0006] Divide the target pipeline according to the predetermined interval to obtain N pipeline nodes;

[0007] If there is no air column in any node to be measured of the target pipeline relative to the target pipeline, the pressure and flow rate of the pipeline node adjacent to the node to be measured are calculated according to a predetermined pipe flow elastic water body control algorithm to obtain the target flow rate and target pressure of the node to be measured;

[0008] If there is an air column in any node to be measured of the target pipeline relative to the target pipeline, the target flow and target pressure of the node to be measured are determined according to the predetermined pipe flow elastic water body control algorithm, air column control algorithm and gas-liquid boundary control algorithm.

[0009] In a possible implementation, if there is an air column in any node to be measured of the target pipeline relative to the target pipeline, the method further includes:

[0010] The target pressure of the node to be measured is the same as the absolute pressure of the air column obtained according to the air column control algorithm.

[0011] In a possible implementation, the pipe flow elastic water body control algorithm is:

[0012] HP =H A +MQ A -(M+R|Q A |)Q P

[0013] H P =H B -MQ B +(M+R|Q B |)Q P

[0014] Where M and R are constants, H A is the water pressure at the pipeline node A, H B is the water pressure at pipeline node B, Q A is the water flow rate at pipeline node A, Q B is the water flow at pipeline node B, H p is the target pressure of the node P to be tested, Q P is the target flow of the node P to be tested, pipeline node A and pipeline node B are the pipeline nodes adjacent to the node P to be tested, and | | is the absolute value.

[0015] In a possible implementation, the gas column control algorithm is:

[0016]

[0017] Among them, Q a is the gas flow rate of the gas column at the gas nozzle, H a is the absolute pressure of the gas column in the target pipeline, V a is the volume of the air column in the pipe, k is the air column constant, and t is the time.

[0018] In a possible implementation, the method further includes:

[0019] When H a <1.89H atm When H is , the gas flow rate of the gas column at the gas nozzle is calculated based on the gas nozzle cross-sectional area, the atmospheric pressure outside the target pipeline, and the resistance coefficient of the gas nozzle according to the gas nozzle control algorithm, where H atm is the atmospheric pressure outside the target pipe.

[0020] In a possible implementation, the air nozzle control algorithm is:

[0021]

[0022] Among them, Q a is the gas flow rate of the gas column at the gas nozzle, A0 is the cross-sectional area of ​​the gas nozzle, ρ w is the water density, ρ a is the air density, Hatm is the atmospheric pressure outside the target pipe, Y is the resistance coefficient of airflow through the nozzle, C d is the flow coefficient of the nozzle, H a is the absolute pressure of the gas column in the target pipeline, and g is the acceleration due to gravity.

[0023] In a possible implementation, the gas-liquid boundary control algorithm is:

[0024]

[0025] H wa =H a

[0026] Among them, x w is the target distance between the gas-liquid boundary position and the upstream node of the pipeline at time t+△t, x w0 is the distance between the gas-liquid boundary at time t and the most upstream node of the pipeline, v c is the water velocity at time t, H wa is the water pressure at time t, △t is the time interval, H a is the absolute pressure of the gas column in the target pipeline.

[0027] In a second aspect, a pipeline gas column transient flow simulation device is provided, which may include:

[0028] The device comprises:

[0029] A partitioning unit, used for partitioning the target pipeline according to a predetermined interval to obtain N pipeline nodes;

[0030] A calculation unit, configured to calculate the pressure and flow rate of a pipeline node adjacent to the node to be measured, if no air column exists in any node to be measured relative to the target pipeline, according to a predetermined pipe flow elastic water body control algorithm, to obtain a target flow rate and a target pressure of the node to be measured;

[0031] Furthermore, if there is an air column in any node to be measured of the target pipeline relative to the target pipeline, the target flow and target pressure of the node to be measured are determined according to a predetermined pipe flow elastic water body control algorithm, an air column control algorithm and a gas-liquid boundary control algorithm.

[0032] In a third aspect, an electronic device is provided, the electronic device comprising a processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface and the memory communicate with each other via the communication bus;

[0033] Memory, used to store computer programs;

[0034] The processor is used to implement any method step described in the first aspect when executing the program stored in the memory.

[0035] In a fourth aspect, a computer-readable storage medium is provided, wherein a computer program is stored in the computer-readable storage medium, and when the computer program is executed by a processor, any method step described in the first aspect is implemented.

[0036] The present application provides a method for simulating transient flow of gas column in pipeline, the method comprising: constructing elastic water control algorithm of pipe flow, gas column control algorithm and gas-liquid boundary control algorithm; establishing interface tracking algorithm at the gas-liquid boundary according to transient state of gas-liquid boundary and rigid interface assumption; using moving wave to calculate target pressure and target flow of any node to be measured of target pipeline; the present application can enrich and improve the theory of transient flow of impact gas column by combining the characteristics and mechanism of transient flow of impact gas column at the end of pipeline nozzle, and provide scientific judgment for engineering design and safe operation of water pipeline system. The present application can also solve the problem of poor simulation effect between existing models and actual engineering, so as to simulate actual pipeline in real time and accurately, grasp engineering hydraulic safety problems, and timely check hidden dangers, which is of great significance to the safety and stability of actual pipeline. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the embodiments of the present application will be briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present application and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0038] Figure 1 A system architecture diagram for a method for simulating transient flow of a gas column in a pipeline provided in an embodiment of the present application;

[0039] Figure 2 A schematic diagram of a flow chart of a pipeline gas column transient flow simulation method provided in an embodiment of the present application;

[0040] Figure 3 A schematic diagram of pipeline node division provided in an embodiment of the present application;

[0041] Figure 4 Schematic diagram of the experimental platform provided for the embodiment of the present application;

[0042] Figure 5 A comparison chart of the pressure fluctuation curves of the experiment and numerical simulation at FT-1#, a working condition provided in the embodiment of the present application;

[0043] Figure 6A comparison chart of the pressure fluctuation curves of the experiment and numerical simulation at FT-1# of the working condition 2 provided in the embodiment of the present application;

[0044] Figure 7 A schematic diagram of the structure of a pipeline gas column transient flow simulation device provided in an embodiment of the present application;

[0045] Figure 8 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0046] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments of the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0047] For ease of understanding, the terms involved in the embodiments of the present application are explained below:

[0048] Rigid interface hypothesis: In order to avoid the tediousness and complexity of the elastic water model in tracing the gas-liquid boundary, the following assumption will be introduced, that is, the mechanical inertia of the thin layer of water column adjacent to the gas-liquid boundary is ignored, that is, the thin layer of water column is regarded as a rigid water body, which is called the "rigid interface" hypothesis.

[0049] Transient flow is the transition process of unstable flow state when water flow changes from one stable state to another stable state.

[0050] The pipeline gas column transient flow simulation method provided in the embodiment of the present application can be applied to Figure 1 In the system architecture shown in Figure 1 As shown, the system may include: a server and a terminal. The server may be a physical server, or a server cluster or distributed system composed of multiple physical servers, or a cloud server that provides cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, content delivery networks (Content Delivery Network, CDN), as well as basic cloud computing services such as big data and artificial intelligence platforms. The terminal may be a sensor that can detect data such as water flow rate and head pressure. The terminal and the server may be directly or indirectly connected via wired or wireless communication methods, which is not limited in this application.

[0051] The terminal is used to receive data such as the water flow rate, water head pressure, pipe diameter of the target pipe at the water head position of the target pipe, and air column flow at the air nozzle position at the end of the target pipe, and send the above data to the server.

[0052] The server is used to receive data sent by the terminal and execute a pipeline air column transient flow simulation method provided in the present application.

[0053] Water pipelines are an indispensable part of hydropower stations, long-distance water diversion projects, urban pipe networks and other projects. At present, the design of water pipelines pays little attention to air columns, and the hazards of air columns are often ignored. When air columns remain in water pipelines, they are destructive to the pipelines and can easily cause engineering accidents.

[0054] The existing transient flow model of the impact air column considering the air nozzle does not comprehensively consider whether the water body is compressible, whether the air column is compressible, and the dynamic capture of the gas-liquid boundary. As a result, the boundary conditions and initial conditions used in the numerical simulation are limited, and the simulation effect on more complex practical problems is poor.

[0055] Therefore, the present application provides a method for simulating transient flow of pipeline air columns, specifically studying the water-gas interaction laws of air columns in water pipelines and their influence on the transient of water pipelines, clarifying and improving the research on the transient flow mechanism of air columns, and providing scientific judgment basis for the engineering design and safe operation of water pipeline systems.

[0056] The preferred embodiments of the present application are described below in conjunction with the drawings in the specification. It should be understood that the preferred embodiments described herein are only used to illustrate and explain the present application, and are not used to limit the present application. In addition, the embodiments and features in the embodiments of the present application may be combined with each other if there is no conflict.

[0057] Figure 2 The following is a flow chart of a pipeline gas column transient flow simulation method provided in an embodiment of the present application. Figure 2 As shown, the method may include:

[0058] Step S210: Divide the target pipeline according to a predetermined interval to obtain N pipeline nodes.

[0059] Specifically, due to the calculation characteristics of the moving wave, the target pipeline needs to be divided into segments with equal intervals to determine the position of the node to be measured (computation node) in the target pipeline.

[0060] The target pipeline is divided into N-1 segments with equal spacing according to a predetermined spacing Δx, including the head and tail nodes of the target pipeline, to obtain N pipeline nodes, that is, N computing nodes.

[0061] Step S220: Determine the target flow and target pressure of the node to be tested according to a predetermined integrated control algorithm.

[0062] Among them, the integrated algorithms include pipe flow elastic water body control algorithm, air column control algorithm and gas-liquid boundary control algorithm.

[0063] It should be noted that the node to be tested is any pipeline node.

[0064] Before executing step S220, A, the control equation of the elastic water body of the pipe flow, can be specifically expressed as:

[0065]

[0066]

[0067] Among them, H is the head pressure of the target pipe, V is the water velocity, c is the propagation speed of the pipe flow wave, and g is the gravitational acceleration; combined with Figure 3 The horizontal and vertical axes are shown, is the horizontal distance in the mathematical sense, is the mathematical time, D is the pipe diameter, f is the pipe friction coefficient, To find the sign of the partial derivative.

[0068] Afterwards, the control equations of the elastic water body of the pipe flow (1) and (2) are transformed into ordinary differential forms, and the moving wave equation is:

[0069]

[0070]

[0071] That is to say, the above-mentioned pipe flow elastic water body control equation can be understood as a relationship equation between the parameters related to the water body in the target pipe, and d is the symbol of the differential.

[0072] Then, combined with Figure 3 As shown, the calculation node is integrated along the + and - moving waves to obtain the elastic water control algorithm of the pipe flow, which is as follows:

[0073] H P =H A +MQ A -(M+R|Q A |)Q P (5)

[0074] H P =H B -MQ B +(M+R|Q B |)Q P (6)

[0075] Among them, M and R are constants, H A is the water pressure at the pipeline node A, H B is the water pressure at pipeline node B, Q A is the water flow at pipeline node A, QB is the water flow at pipeline node B, H p is the target pressure of the node P to be tested, Q P is the target flow of the node P to be measured, pipeline node A and pipeline node B are the pipeline nodes adjacent to the node P to be measured, and || is the absolute value.

[0076] Furthermore, the calculation process of constants M and R is M = c / gA; R = f△x / (2gDA 2 ), where g is the acceleration of gravity, A is the cross-sectional area of ​​the target pipeline, c is the velocity of pipe flow wave propagation, f is the pipeline friction coefficient, △x is the spacing between the pipeline nodes of the target pipeline, △x=c△t, and △t is the time interval.

[0077] B. Air column control algorithm, which can be specifically expressed as:

[0078]

[0079] Among them, Q a is the gas flow rate of the gas column at the gas nozzle, H a is the absolute pressure of the gas column in the target pipeline, V a is the volume of the air column in the pipeline, k is the air column constant, generally k = 1.4, and t is the time.

[0080] Furthermore, the gas flow rate of the gas column at the gas nozzle is divided into the following two cases:

[0081] Case 1: When the air flows out at a subsonic speed, that is, H a <1.89H atm The gas flow rate of the gas column at the gas nozzle is calculated based on the gas nozzle control algorithm, the gas nozzle cross-sectional area, the atmospheric pressure outside the target pipeline, and the resistance coefficient of the gas nozzle, where H atm is the atmospheric pressure outside the target pipe.

[0082] The air nozzle control algorithm is:

[0083]

[0084] Among them, Q a is the gas flow rate of the gas column at the gas nozzle, A0 is the cross-sectional area of ​​the gas nozzle, ρ w is the water density, ρ a is the air density, Y is the resistance coefficient of airflow through the air nozzle, C d is the flow coefficient of the nozzle, H a is the absolute pressure of the gas column in the target pipeline.

[0085] Furthermore, the resistance coefficient Y of the airflow through the air nozzle is solved by the following formula:

[0086]

[0087] Where n is the polynomial index of the compressible flow through the nozzle, generally n = 1.4. Although the polynomial index n and k in the previous equation have the same value, their meanings are different.

[0088] Case 2: When the air flows out at the critical flow rate, that is, H a ≥1.89H atm , the air nozzle is blocked, and its air flow rate remains unchanged and no longer changes with the air column pressure in the pipeline. At this time, the resistance coefficient Y of the air flow through the air nozzle is solved by the following formula:

[0089]

[0090] C. Determine the gas-liquid boundary control algorithm, which can be specifically expressed as:

[0091]

[0092] H wa =H a (12)

[0093] Among them, x w is the target distance between the gas-liquid boundary position and the upstream node at time t+△t, x w0 is the distance between the gas-liquid boundary position and the most upstream node at time t, v c is the water velocity at time t, H wa is the water pressure at time t.

[0094] Step S220 specifically includes: if there is no air column in any node to be tested of the target pipeline relative to the target pipeline, indicating that the fluid inside the node to be tested relative to the target pipeline is water, then according to a predetermined pipe flow elastic water body control algorithm, the pressure and flow of the pipeline node adjacent to the node to be tested are calculated to obtain the target flow and target pressure of the node to be tested;

[0095] Specifically, by combining equations (5) and (6), the target flow rate and target pressure of the node to be measured relative to the water body in the target pipeline can be obtained.

[0096]

[0097]

[0098] Among them, A1 is the first characteristic parameter of pipeline node A, B1 is the first characteristic parameter of pipeline node B; A2 is the second characteristic parameter of pipeline node A, and B2 is the second characteristic parameter of pipeline node B. It should be noted that pipeline node A is the upstream node of the node P to be tested, and pipeline node B is the downstream node of the node P to be tested.

[0099] Furthermore, A1=H A +MQ A ; B1=H B -MQ B ; A2=M+R|Q A |;B2=M+R|Q B |.

[0100] If there is an air column in any node to be tested of the target pipeline relative to the target pipeline, it indicates that there is gas in the node to be tested relative to the target pipeline. Then, the target flow and target pressure of the node to be tested are determined according to the predetermined pipe flow elastic water body control algorithm, air column control algorithm and gas-liquid boundary control algorithm.

[0101] Specifically, there is a gas column in the target pipeline relative to the node to be tested. First, the gas-liquid boundary position is determined according to the gas-liquid boundary control algorithm.

[0102] Secondly, for the water body part, the water pressure and water flow are determined according to the pipe flow elastic water body control algorithm.

[0103] Secondly, for the gas column part, the gas pressure and gas flow in the gas column will be affected by the relevant parameters of the gas nozzle at the end of the target pipeline. Therefore, the target pressure and target flow of the gas column are finally determined according to the gas column control algorithm.

[0104] Finally, for the gas-liquid boundary part, according to the transient state of the gas-liquid boundary and the rigid interface assumption, an interface tracking algorithm is established at the gas-liquid boundary to determine the pressure and flow rate of the gas-liquid boundary.

[0105] Furthermore, according to the rigid interface assumption, the pressure of the water body and the air column at the rigid interface is equal, and the flow rate of inflow and outflow is equal. The interface tracking algorithm can be specifically expressed as:

[0106] H P’ =H a (15)

[0107] Q P’ =Q a (16)

[0108] Among them, H a is the gas pressure at time t+△t; Q a is the gas flow rate at time t+△t, H P’ is the pressure at the gas-liquid boundary at time t+△t, Q P’ It is the flow rate at the gas-liquid boundary at time t+△t.

[0109] That is to say, the present application can determine data such as the gas-liquid boundary position at a future time (according to a time interval), the target flow rate and target pressure of any pipeline node in a pipeline containing a gas column, and the like.

[0110] Example 1: The experimental results of the hydraulics laboratory of Hohai University are selected to verify the effectiveness of this application. Figure 4 As shown. The horizontal pipe length is 8.83m, the inner diameter is 0.040m, the average resistance coefficient of the target pipe is 0.085, and the measured pipe flow wave velocity is 400m / s. This application uses two sets of working conditions with different boundary conditions for model verification.

[0111] Condition 1: The upstream of the target pipeline is set as the pressure inlet boundary (pressure tank), the pressure is 0.08MPa, and the initial gas column length in the pipeline is 1.058m (i.e. Figure 4 The 1# ball valve is closed, and the 2# and 3# ball valves are opened. It is considered that the air column from the 1# ball valve to the air nozzle at the end of the pipeline is full of air. The diameter of the air nozzle is 1mm. At the beginning of the experiment, the 1# ball valve is opened instantaneously, the pipeline begins to be filled, and the water impacts the air column in the pipeline, causing a transient change.

[0112] Comparing the calculation results of this application with the experimental data results, we can obtain the following Figure 5 The pressure fluctuation comparison curve at the sensor (pressure sensor) position shown in the figure is Figure 5 It can be seen that the present application can well predict the fluctuation trend, cycle and pressure peak of the pipeline pressure change curve, especially the first pressure peak.

[0113] Condition 2: The upstream of the target pipeline is set as the pressure inlet boundary, the pressure is 0.08 MPa, and the initial gas column length in the pipeline is 6.45 m (i.e. Figure 4 The 3# ball valve is closed, and the 1# and 2# ball valves are opened. It is considered that the air column from the 3# ball valve to the air nozzle at the end of the pipeline is full of air. The diameter of the air nozzle is changed to 2mm. At the beginning of the experiment, the 3# ball valve is opened instantaneously, the pipeline begins to be filled, and the water impacts the air column in the pipeline, causing a transient change.

[0114] Comparing the calculation results of this application with the experimental data results, we can obtain the following Figure 6 A comparison graph of pressure fluctuations at the sensor (pressure sensor) location shown.

[0115] For the simulation of transient flow of impact gas column of gas nozzle at the end of pipeline in actual engineering, the accuracy of system pressure peak simulation is the first consideration, which is related to whether the maximum pressure of pipeline exceeds the standard value of pipeline design and whether it may cause accidents such as pipe burst. Therefore, according to the results of Example 1, it can be seen that the present application can better simulate the transient flow problem of impact gas column of gas nozzle at the end of pipeline, especially the simulation result of pressure peak is more accurate.

[0116] The present application provides a method for simulating transient flow of gas column in pipeline, the method comprising: constructing elastic water control algorithm of pipe flow, gas column control algorithm and gas-liquid boundary control algorithm; establishing interface tracking algorithm at the gas-liquid boundary according to transient state of gas-liquid boundary and rigid interface assumption; using moving wave to calculate target pressure and target flow of any node to be measured of target pipeline; the present application can enrich and improve the theory of transient flow of impact gas column by combining the characteristics and mechanism of transient flow of impact gas column at the end of pipeline nozzle, and provide scientific judgment for engineering design and safe operation of water pipeline system. The present application can also solve the problem of poor simulation effect between existing models and actual engineering, so as to simulate actual pipeline in real time and accurately, grasp engineering hydraulic safety problems, and timely check hidden dangers, which is of great significance to the safety and stability of actual pipeline.

[0117] Corresponding to the above method, the present application embodiment also provides a pipeline gas column transient flow simulation device, such as Figure 7 As shown, the device comprises:

[0118] A division unit 710 is used to divide the target pipeline according to a predetermined interval to obtain N pipeline nodes;

[0119] The calculation unit 720 is used to calculate the pressure and flow rate of the pipeline node adjacent to the node to be measured according to a predetermined pipe flow elastic water body control algorithm if there is no air column in any node to be measured of the target pipeline relative to the target pipeline, so as to obtain the target flow rate and target pressure of the node to be measured;

[0120] Furthermore, if there is an air column in any node to be measured of the target pipeline relative to the target pipeline, the target flow and target pressure of the node to be measured are determined according to a predetermined pipe flow elastic water body control algorithm, an air column control algorithm and a gas-liquid boundary control algorithm.

[0121] The functions of each functional unit of a pipeline air column transient flow simulation device provided in the above embodiment of the present application can be achieved through the above method steps. Therefore, the specific working process and beneficial effects of each unit in a pipeline air column transient flow simulation device provided in the embodiment of the present application will not be repeated here.

[0122] The present application also provides an electronic device, such as Figure 8 As shown, it includes a processor 810 , a communication interface 820 , a memory 830 and a communication bus 840 , wherein the processor 810 , the communication interface 820 , and the memory 830 communicate with each other via the communication bus 840 .

[0123] Memory 830, used for storing computer programs;

[0124] The processor 810 is used to execute the program stored in the memory 830 to implement the following steps:

[0125] Divide the target pipeline according to the predetermined interval to obtain N pipeline nodes;

[0126] If there is no air column in any node to be measured of the target pipeline relative to the target pipeline, the pressure and flow rate of the pipeline node adjacent to the node to be measured are calculated according to a predetermined pipe flow elastic water body control algorithm to obtain the target flow rate and target pressure of the node to be measured;

[0127] If there is an air column in any node to be measured of the target pipeline relative to the target pipeline, the target flow and target pressure of the node to be measured are determined according to a predetermined pipe flow elastic water body control algorithm, an air column control algorithm and a gas-liquid boundary control algorithm.

[0128] The communication bus mentioned above can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. The communication bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, only one thick line is used in the figure, but it does not mean that there is only one bus or one type of bus.

[0129] The communication interface is used for communication between the above electronic device and other devices.

[0130] The memory may include a random access memory (RAM) or a non-volatile memory (NVM), such as at least one disk memory. Optionally, the memory may also be at least one storage device located away from the aforementioned processor.

[0131] The above-mentioned processor can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it can also be a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field programmable gate array (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.

[0132] The implementation methods and beneficial effects of the components of the electronic device in the above embodiments to solve the problems can be seen in Figure 2 The various steps in the illustrated embodiment are implemented, therefore, the specific working process and beneficial effects of the electronic device provided by the embodiment of the present application are not repeated here.

[0133] In another embodiment provided in the present application, a computer-readable storage medium is also provided, in which instructions are stored. When the computer-readable storage medium is run on a computer, the computer executes a pipeline air column transient flow simulation method described in any of the above embodiments.

[0134] In another embodiment provided in the present application, a computer program product including instructions is also provided, which, when executed on a computer, enables the computer to execute a pipeline air column transient flow simulation method described in any one of the above embodiments.

[0135] Those skilled in the art will appreciate that the embodiments in the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt a complete hardware embodiment, a complete software embodiment, or a form of an embodiment combining software and hardware. Moreover, the present application may adopt a form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0136] The embodiments of the present application are described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0137] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0138] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process in the computer or other programmable device. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0139] Although the preferred embodiments of the present application have been described, those skilled in the art may make additional changes and modifications to these embodiments once they are aware of the basic creative concepts. Therefore, the present application embodiments are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the present application embodiments.

[0140] Obviously, those skilled in the art can make various changes and modifications to the embodiments in the present application without departing from the spirit and scope of the embodiments in the present application. Thus, if these modifications and variations of the embodiments in the present application are within the scope of the embodiments in the present application and their equivalents, the embodiments in the present application are also intended to include these modifications and variations.

Claims

1. A method for simulating transient flow of pipeline gas column, characterized in that: The method comprises: Divide the target pipeline according to the predetermined interval to obtain N pipeline nodes; If there is no air column in any node to be measured of the target pipeline relative to the target pipeline, the pressure and flow rate of the pipeline node adjacent to the node to be measured are calculated according to a predetermined pipe flow elastic water body control algorithm to obtain the target flow rate and target pressure of the node to be measured; If there is an air column in any node to be measured of the target pipeline relative to the target pipeline, the target flow and target pressure of the node to be measured are determined according to a predetermined pipe flow elastic water body control algorithm, an air column control algorithm and a gas-liquid boundary control algorithm; The pipe flow elastic water body control algorithm is: H P1 =H A +MQ A -(M+R|Q A |)Q P H P2 =H B -MQ B +(M+R|Q B |)Q P Where M and R are constants, H A is the water pressure at the pipeline node A, H B is the water pressure at pipeline node B, Q A is the water flow at pipeline node A, Q B is the water flow at pipeline node B, H p1 To determine the target pressure of the node P to be tested through the pipeline node A on the upstream side of the node P to be tested, H p2 To determine the target pressure of the node P to be tested through the pipeline node B on the downstream side of the node P to be tested, H p1 =H p2 , Q P is the target flow of the node P to be tested, || is the absolute value; The gas column control algorithm is: Among them, Q a is the gas flow rate of the gas column at the gas nozzle, H a is the absolute pressure of the gas column in the target pipeline, V a is the volume of the air column in the pipe, k is the air column constant, and t is the time; The method further comprises: When H a <1.89H atm When H is , the gas flow rate of the gas column at the gas nozzle is calculated based on the gas nozzle cross-sectional area, the atmospheric pressure outside the target pipeline, and the resistance coefficient of the gas nozzle according to the gas nozzle control algorithm, where H atm is the atmospheric pressure outside the target pipeline; The air nozzle control algorithm is: Among them, Q a is the gas flow rate of the gas column at the gas nozzle, A0 is the cross-sectional area of ​​the gas nozzle, ρ w is the water density, ρ a is the air density, H atm is the atmospheric pressure outside the target pipe, Y is the resistance coefficient of airflow through the nozzle, C d is the flow coefficient of the nozzle, H a is the absolute pressure of the gas column in the target pipeline, and g is the acceleration due to gravity; The gas-liquid boundary control algorithm is: H wa =H a Among them, x w is the target distance between the gas-liquid boundary position and the upstream node of the pipeline at time t+△t, x w0 is the distance between the gas-liquid boundary at time t and the most upstream node of the pipeline, v c is the water velocity at time t, H wa is the water pressure at time t, △t is the time interval, H a is the absolute pressure of the gas column in the target pipeline; Wherein, if there is an air column in any node to be measured of the target pipeline relative to the target pipeline, the method further includes: The target pressure of the node to be measured is the same as the absolute pressure of the air column obtained according to the air column control algorithm.

2. A pipeline gas column transient flow simulation device, characterized in that: The device comprises: A partitioning unit, used for partitioning the target pipeline according to a predetermined interval to obtain N pipeline nodes; A calculation unit, configured to calculate the pressure and flow rate of a pipeline node adjacent to the node to be measured, if no air column exists in any node to be measured relative to the target pipeline, according to a predetermined pipe flow elastic water body control algorithm, to obtain a target flow rate and a target pressure of the node to be measured; And, if there is an air column in any node to be measured of the target pipeline relative to the target pipeline, the target flow and target pressure of the node to be measured are determined according to a predetermined pipe flow elastic water body control algorithm, an air column control algorithm and a gas-liquid boundary control algorithm; The pipe flow elastic water body control algorithm is: H P1 =H A +MQ A -(M+R|Q A |)Q P H P2 =H B -MQ B +(M+R|Q B |)Q P Where M and R are constants, H A is the water pressure at the pipeline node A, H B is the water pressure at pipeline node B, Q A is the water flow at pipeline node A, Q B is the water flow at pipeline node B, H p1 To determine the target pressure of the node P to be tested through the pipeline node A on the upstream side of the node P to be tested, H p2 To determine the target pressure of the node P to be tested through the pipeline node B on the downstream side of the node P to be tested, H p1 =H p2 , Q P is the target flow of the node P to be tested, || is the absolute value; The gas column control algorithm is: Among them, Q a is the gas flow rate of the gas column at the gas nozzle, H a is the absolute pressure of the gas column in the target pipeline, V a is the volume of the air column in the pipe, k is the air column constant, and t is the time; Also includes: When H a <1.89H atm When H is , the gas flow rate of the gas column at the gas nozzle is calculated based on the gas nozzle cross-sectional area, the atmospheric pressure outside the target pipeline, and the resistance coefficient of the gas nozzle according to the gas nozzle control algorithm, where H atm is the atmospheric pressure outside the target pipeline; The air nozzle control algorithm is: Among them, Q a is the gas flow rate of the gas column at the gas nozzle, A0 is the cross-sectional area of ​​the gas nozzle, ρ w is the water density, ρ a is the air density, H atm is the atmospheric pressure outside the target pipe, Y is the resistance coefficient of airflow through the nozzle, C d is the flow coefficient of the nozzle, H a is the absolute pressure of the gas column in the target pipeline, and g is the acceleration due to gravity; The gas-liquid boundary control algorithm is: H wa =H a Among them, x w is the target distance between the gas-liquid boundary position and the upstream node of the pipeline at time t+△t, x w0 is the distance between the gas-liquid boundary at time t and the most upstream node of the pipeline, v c is the water velocity at time t, H wa is the water pressure at time t, △t is the time interval, H a is the absolute pressure of the gas column in the target pipeline; Wherein, if there is an air column in any node to be measured of the target pipeline relative to the target pipeline, the method further includes: The target pressure of the node to be measured is the same as the absolute pressure of the air column obtained according to the air column control algorithm.

3. An electronic device, characterized in that: The electronic device comprises a processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface and the memory communicate with each other via the communication bus; Memory, used to store computer programs; The processor is used to implement the method steps described in claim 1 when executing the program stored in the memory.

4. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method steps described in claim 1 are implemented.

Images