Transient simulation method, device and equipment for liquid fuel pipeline and medium

Abstract

The invention relates to the technical field of liquid fuel pipeline simulation, and discloses a transient simulation method, device and equipment for a liquid fuel pipeline and a medium, which can model a liquid fuel pipeline system to obtain a liquid fuel pipeline model. Spatial discretization is performed on each connecting pipeline to determine at least one in-segment point in each connecting pipeline. Constructing a first operation equation set corresponding to each boundary point according to the station yard type corresponding to each boundary point; and according to the operation characteristics of the points in each pipe section, constructing a second operation equation set corresponding to the points in each pipe section. And based on each first operation equation set and each second operation equation set, liquid fuel pipeline transient simulation is carried out. The transient simulation precision of the liquid fuel pipeline can be improved, the application range of the liquid fuel pipeline can be widened, and reliable support is provided for safe operation of the pipeline under multi-working-condition switching.

Description

Technical Field

[0001] This invention relates to the field of liquid fuel pipeline simulation technology, and in particular to a transient simulation method, apparatus, equipment and medium for liquid fuel pipelines. Background Technology

[0002] Liquid fuel pipelines, as core infrastructure for long-distance and efficient energy transportation, are characterized by long distances, multiple pumping stations, multiple distribution outlets, and multiple operating conditions, and are widely used in cross-regional energy allocation and other scenarios.

[0003] Long-distance liquid pipelines frequently encounter transient conditions during operation, such as pump shutdowns, valve opening and closing, flow fluctuations, and changes in medium properties. These conditions can easily lead to problems like water hammer and sudden pressure spikes and drops, potentially causing pipeline rupture, equipment damage, or even safety accidents. Simulation studies of multi-condition switching in liquid fuel pipelines are a key technical means to optimize pipeline design and formulate safety control strategies.

[0004] However, the simulation methods of related technologies lack accuracy and have limited applicability. The simulation results deviate significantly from actual operating conditions, failing to provide reliable support for the safe operation of pipelines under multiple operating conditions. Summary of the Invention

[0005] This invention provides a transient simulation method, apparatus, equipment, and medium for liquid fuel pipelines, which addresses the shortcomings of related technologies, such as insufficient accuracy, limited applicability, large deviation between simulation results and actual operating conditions, and inability to provide reliable support for the safe operation of pipelines under multiple operating conditions. This invention improves the accuracy and applicability of transient simulation for liquid fuel pipelines, providing reliable support for the safe operation of pipelines under multiple operating conditions.

[0006] In a first aspect, the present invention provides a transient simulation method for a liquid fuel pipeline, comprising: A liquid fuel pipeline system is modeled to obtain a liquid fuel pipeline model; wherein, the liquid fuel pipeline model includes multiple boundary points and connecting pipelines between adjacent boundary points, each boundary point uniquely corresponds to the first station, intermediate station or terminal station in the liquid fuel pipeline system, and each connecting pipeline uniquely corresponds to a pipe segment between adjacent stations; Each of the connecting pipelines is spatially discretized to determine at least one point within a segment of each connecting pipeline; Based on the station type corresponding to each boundary point, a first set of operating equations is constructed for each boundary point; and based on the operating characteristics of each point within the pipe segment, a second set of operating equations is constructed for each point within the pipe segment. Transient simulations of liquid fuel pipelines are performed based on each of the first set of operating equations and each of the second set of operating equations.

[0007] Optionally, the transient simulation of the liquid fuel pipeline based on each of the first set of operating equations and each of the second set of operating equations includes: Based on the operating status of each station and each pipe segment, the operating parameters in each of the first and second operating equation sets are initialized to obtain the initialization parameter sets for each boundary point and each point within the pipe segment; Based on each of the first set of operating equations, each of the second set of operating equations, and the equipment status changes in each of the stations, update the initialization parameter set of each of the boundary points and each of the points within the pipe segment to obtain a new parameter set for each of the boundary points and each of the points within the pipe segment. Based on the new parameter set for each of the first set of operating equations, each of the second set of operating equations, each of the boundary points, and each point within the pipe segment, determine and output the flow rate and head for each of the boundary points and each point within the pipe segment; Based on each of the first set of operating equations, each of the second set of operating equations, the equipment status changes in each of the stations, the flow rate and head of each of the boundary points and each of the points within the pipe segment, update the new parameter set for each of the boundary points and each of the points within the pipe segment until the set iteration termination condition is met.

[0008] Optionally, before initializing the operating parameters in each of the first and second sets of operating equations based on the operating status of each of the stations and each of the pipe segments, the method further includes: In response to the simulation run command, the simulation time step timing starts from 0; When the simulation time step is 0, the step of initializing the operating parameters in each of the first set of operating equations and each of the second set of operating equations according to the operating status of each of the stations and each of the pipe segments is executed; After updating the initial parameter sets for each boundary point and each point within the pipe segment to obtain new parameter sets for each boundary point and each point within the pipe segment, the method further includes: Increment the current simulation time step by 1 to obtain the new simulation time step; Determine whether the new simulation time step has reached the set maximum time step; If not, then it is determined that the iteration termination condition has not been met, and the step of determining and outputting the flow rate and head of each boundary point and each point in the pipe segment based on the new parameter set of each first set of running equations, each second set of running equations, each boundary point and each point in the pipe segment is executed; If so, then the iteration termination condition is satisfied, and the step of determining and outputting the flow rate and head of each boundary point and each point in the pipe segment based on the new parameter set of each first set of running equations, each second set of running equations, each boundary point and each point in the pipe segment is prohibited.

[0009] Optionally, when the station type corresponding to the boundary point is an intermediate station, the first set of operating equations corresponding to the boundary point is: ; in, and The first The inbound and outbound traffic of each intermediate station For the first The injected or distributed flow of an intermediate station; and The first The inbound and outbound pressure heads of each intermediate station; and For the first Characteristic curve parameters of the pump unit at each intermediate station; and The first The sum of valve resistance coefficients before and after the injection or distribution point of each intermediate station; , For the first The point within the second-to-last pipe segment upstream of the intermediate station Characteristic line coefficients; , For the first The second pipe section downstream of the intermediate station Characteristic line coefficients.

[0010] Optionally, when the station type corresponding to the boundary point is the first station, the first set of operating equations corresponding to the boundary point is either the first station constant flow control equation set or the first station constant pressure control equation set, wherein: The control equations for the first station's flow regulation are as follows: ; in, For outbound traffic, For the set injection flow, For exiting the station, , The second inner point of the downstream section of the first station coefficients of characteristics For the pressure head entering the station, This is the sum of the valve resistance coefficients after the initial injection point. , The characteristic curve parameters of the pump unit at the first station; The set of constant pressure control equations for the first station is as follows: ; This is the set value.

[0011] Optionally, when the station type corresponding to the boundary point is a terminal station, the first set of operating equations corresponding to the boundary point is a terminal station constant flow control equation set or a terminal station constant pressure control equation set. The set of control equations for the terminal station constant flow is as follows: ; in, For inbound traffic, For the set distribution flow, For the pressure head entering the station, , The second to last inner point of the upstream section of the terminal station coefficients of characteristics This is the sum of the valve resistance coefficients before the final distribution point. For exiting the station; The set of control equations for the terminal station constant pressure is as follows: ; in, This is the set value.

[0012] Optionally, the second set of operating equations corresponding to the points within the pipe segment is: ; or ; ; ; in, The first segment of the pipe Flow rate at points within a pipe section; The first segment of the pipe Pressure head at points within a pipe section; and The first segment of the pipe Each internal point Characteristic line coefficients; and For the first Each internal point Characteristic line coefficients; , These are the pipe segments at the 1st step of the previous calculation. The, the Flow rate at points within a pipe section; and These are the pipe segments at the 1st step of the previous calculation. The, the The pressure head at each internal point; and For the pipe segment in the previous calculation step The, the Corrected wave velocity at points within each pipe section; It is the acceleration due to gravity; and This refers to the cross-sectional area of ​​the pipe. , For the pipe segment in the previous calculation step The, the The coefficient of friction at each internal point; The wave velocity of the pressure wave. For time step.

[0013] In a second aspect, the present invention provides a transient simulation device for a liquid fuel pipeline, comprising: A modeling unit is used to model a liquid fuel pipeline system to obtain a liquid fuel pipeline model; wherein, the liquid fuel pipeline model includes multiple boundary points and connecting pipelines between adjacent boundary points, each boundary point uniquely corresponds to the first station, intermediate station or terminal station in the liquid fuel pipeline system, and each connecting pipeline uniquely corresponds to a pipe segment between adjacent stations; A spatial discretization unit is used to spatially discretize each of the connecting pipelines to determine at least one point within a segment of each of the connecting pipelines. The first construction unit is used to construct a first set of operating equations corresponding to each boundary point according to the station type corresponding to each boundary point; The second construction unit is used to construct a second set of operating equations corresponding to each point in the pipe segment based on the operating characteristics of each point in the pipe segment. The simulation unit is used to perform transient simulations of liquid fuel pipelines based on each of the first set of operating equations and each of the second set of operating equations.

[0014] Thirdly, the present invention provides a computer device, comprising: a memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the computer instructions to perform the liquid fuel pipeline transient simulation method described in the first aspect or any corresponding embodiment thereof.

[0015] Fourthly, the present invention provides a computer-readable storage medium storing computer instructions for causing a computer to execute the liquid fuel pipeline transient simulation method of the first aspect or any corresponding embodiment described above.

[0016] The present invention provides a method, apparatus, equipment, and medium for transient simulation of liquid fuel pipelines. This method can model a liquid fuel pipeline system to obtain a liquid fuel pipeline model. The model includes multiple boundary points and connecting pipelines between adjacent boundary points. Each boundary point uniquely corresponds to the first, intermediate, or final station in the liquid fuel pipeline system, and each connecting pipeline uniquely corresponds to a pipe segment between adjacent stations. Each connecting pipeline is spatially discretized to determine at least one point within a pipe segment. A first set of operating equations is constructed for each boundary point based on its corresponding station type; and a second set of operating equations is constructed for each point within a pipe segment based on its operating characteristics. Transient simulation of the liquid fuel pipeline is then performed based on each first and second set of operating equations. This invention improves the accuracy and applicability of transient simulation of liquid fuel pipelines, providing reliable support for the safe operation of pipelines under multiple operating conditions. Attached Figure Description

[0017] To more clearly illustrate the technical solutions in this invention or related technologies, the accompanying drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, the accompanying drawings described below are some embodiments of this invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0018] Figure 1 A flowchart of a transient simulation method for a liquid fuel pipeline provided in an embodiment of the present invention; Figure 2 This is a schematic diagram of boundary adjustment solution for an intermediate station provided in an embodiment of the present invention; Figure 3 This is a schematic diagram of the boundary adjustment solution for the first station provided in an embodiment of the present invention; Figure 4 This is a schematic diagram of the boundary adjustment solution for a terminal station provided in an embodiment of the present invention; Figure 5 A topology of an oil pipeline system provided in an embodiment of the present invention; Figure 6 Basic data of pipelines and nodes for an oil pipeline system provided in this embodiment of the invention; Figure 7 This invention provides an embodiment of the pressure and flow rate errors of an oil pipeline system under various simulation conditions. Figure 8 This is a schematic diagram illustrating the upstream and downstream changes over time after a pump is spun out, as provided in an embodiment of the present invention. Figure 9 This is a schematic diagram of the structure of a transient simulation device for a liquid fuel pipeline provided in an embodiment of the present invention; Figure 10 This is a schematic diagram of the structure of a computer device provided in an embodiment of the present invention. Detailed Implementation

[0019] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0020] In related technologies, the increased flexibility requirements of downstream users in fuel transportation have led to a significant increase in the frequency of pipeline operating condition switching. This results in diverse, dynamic, and strongly coupled boundary conditions in the transient simulation of liquid fuel pipelines. On the one hand, boundary types include dynamic boundaries related to pump station start-up / shutdown / frequency conversion, constraint boundaries related to flow regulation at distribution / injection points, and parameter boundaries related to changes in pipe material and medium characteristics. On the other hand, during multi-condition switching, the boundaries of pumps, valves, injection points, and distribution points need to be freely combined according to actual transportation requirements, further increasing the complexity of transient simulation.

[0021] Transient simulation methods in related technologies are mostly developed for single or fixed boundary conditions. They generally simplify the system boundary conditions and derive solution methods applicable only to fixed boundary combinations. They cannot adapt to the free combination scenarios of boundaries such as pumps, valves, injection points, and distribution points under multiple operating conditions. This results in insufficient simulation accuracy, limited applicability, and large deviations between simulation results and actual operating conditions, failing to provide reliable support for the safe operation of pipelines under multiple operating conditions.

[0022] Specifically, in handling boundary conditions, related technologies generally construct and solve equations for a single pump station, valve, distribution point, or injection point. For example, for the boundary of a pump station with multiple pumps, a preset pump station characteristic equation is directly substituted, and the pump start-up scheme is set to a fixed mode, remaining unchanged throughout the simulation cycle. For valve boundaries in the middle or at the end of the pipeline, a fixed valve characteristic curve is used, while ignoring flow changes at intermediate injection and distribution points. In actual engineering, the operating conditions of a pumping station are diverse, and the combination of pumps, valves, distribution points, and injection points dynamically changes according to the operating conditions. Related technologies are only applicable to single or fixed operating conditions. When faced with multi-pump switching, dynamic valve adjustment, and fluctuations in injection / distribution flow rates, their applicability is limited, making it difficult to meet the simulation needs of complex dynamic operating conditions in engineering practice.

[0023] Specifically, related technologies also establish a transient simulation model for refined oil pipelines that couples flow mechanisms with operational data, and simulations are performed based on this model. However, the training of this model is highly dependent on data; its performance relies heavily on a large amount of data covering different transient operating conditions. This data needs to include actual operational data or high-precision simulation data. If the data coverage of operating conditions is incomplete, the sample size is insufficient, or the data contains noise interference, it will directly affect the effectiveness of model training and the accuracy of subsequent simulation predictions, increasing the cost and difficulty of data collection and preprocessing. Secondly, facing the diverse boundary conditions of pipeline systems, the model needs to be trained separately for each specific operating condition, forming multiple sets of models corresponding to different operating conditions. This not only significantly increases the workload of model training but also makes the actual application process more complex. Although the trained model is fast in the transient simulation stage and can meet the real-time requirements of engineering, the model training process itself is extremely time-consuming, especially when training for multiple operating conditions. It is necessary to iterate and optimize the model parameters repeatedly to meet the optimization target of the coupling loss function. In addition, the complexity of automatic differential calculation in the process of incorporating physical prior knowledge leads to a significant increase in the overall training cycle. It is difficult to quickly respond to the simulation requirements of pipeline operating conditions or new scenarios, which restricts the flexible application and efficient deployment of the model in engineering practice.

[0024] The following is combined Figures 1-8 This invention describes a transient simulation method for liquid fuel pipelines.

[0025] like Figure 1 As shown in the figure, this embodiment proposes a first transient simulation method for liquid fuel pipelines, which may include the following steps: S101. Model the liquid fuel pipeline system to obtain a liquid fuel pipeline model; wherein, the liquid fuel pipeline model includes multiple boundary points and connecting pipelines between adjacent boundary points, each boundary point uniquely corresponds to the first station, intermediate station or terminal station in the liquid fuel pipeline system, and each connecting pipeline uniquely corresponds to a pipeline segment between adjacent stations.

[0026] Specifically, a liquid fuel pipeline system may include multiple stations, namely the initial station, the final station, and at least one intermediate station. The aforementioned boundary points in the liquid fuel pipeline model correspond one-to-one with the multiple stations in the liquid fuel pipeline system.

[0027] In the liquid fuel pipeline system, any two adjacent stations are connected by a pipe segment. Any two adjacent boundary points in the liquid fuel pipeline model are connecting pipelines, and each connecting pipeline uniquely corresponds to a pipe segment in the liquid fuel pipeline system.

[0028] S102. Spatial discretize each connecting pipeline to determine at least one point within a segment of each connecting pipeline.

[0029] Specifically, in the construction of the simulation model in this embodiment, spatial nodes can be divided into two categories according to the structural characteristics of the pipeline system: one category is the boundary points located at both ends of the pipeline, such as the first station, the last station, and intermediate stations; the other category is the pipeline internal points located between adjacent stations. Since boundary points usually contain equipment components such as pumps and valves, or involve operational processes such as injection and distribution, corresponding calculation methods need to be established for these two types of nodes during the model development stage.

[0030] S103. Construct the first set of operating equations for each boundary point based on the station type corresponding to each boundary point.

[0031] Specifically, regarding the boundary conditions of intermediate stations, in order to adapt to different types of intermediate station sites, this embodiment can first summarize the general form of intermediate stations based on engineering practice.

[0032] See Figure 2 The diagram illustrates the boundary regulation solution for an intermediate station. The intermediate station's core functional components are injection points, distribution points, pump sets, and valves, following a clear functional division and layout rule. Functionally, a station typically has one injection point or one distribution point, and may be equipped with multiple pump sets, inlet shut-off valves / regulating valves, outlet shut-off valves / regulating valves, and other equipment. In terms of layout and function, pump sets are usually located after the injection / distribution point, their core function being to provide the energy required for downstream pipeline transport of the fluid after injection / distribution regulation; regulating valves and shut-off valves are used to control flow rate and pressure.

[0033] This embodiment can characterize the switching of operating conditions at intermediate stations through flexible component combinations and dynamic parameter adjustment. By combining similar equipment before (or after) the distribution point and calculating combined characteristic parameters, for intermediate stations that do not require additional pressurization, the pump unit can be set to a closed state, and its characteristic curve parameters can be adjusted to 0, thus excluding it from fluid dynamics calculations. For intermediate stations requiring fine-grained control, precise control of transient pipeline flow can be achieved by adjusting the pump unit's operating frequency and the opening of the inlet and outlet regulating valves. This design, based on a general form, with flexible component combinations and dynamically adjustable parameters, reduces the limitations of related technologies that require one-to-one modeling of specific intermediate stations. This reduces the workload of repetitive model construction and provides a unified and scalable topological framework for the subsequent derivation of generalized solution formulas.

[0034] use , , Representing the first The inbound flow, outbound flow, and injection (distribution) flow of each intermediate station. Representing the Liquid fuel is being injected into the stations. This indicates that distribution is in progress. Representing the The operating flow rate of a station is 0 or the station does not have injection or distribution functions. For long-distance pipelines, valves and pumps are usually connected in series, and the characteristic parameters of the equipment group are the coefficients of individual equipment. , These represent the sum of the valve resistance coefficients before and after the injection / distribution point, respectively. and These represent the characteristic curve parameters of the pump unit. The inlet of the intermediate station boundary point is connected to the end of the upstream pipe section, and the inflow rate is... With pressure head It can be started from the second to last inner point of the upstream pipe section. The equation of the characteristic line is updated. coefficients of characteristic line equations and Based on the previous step. Similarly, the outlet of the intermediate station boundary point connects to the starting point of the downstream pipeline segment, and the outflow is... With pressure head It can be obtained from the second inner point of the downstream pipe section. The characteristic line equation is updated. Boundary points can be constructed using the pump unit characteristic curves and valve characteristic curves. The inlet pressure head and exit pressure head The relationship.

[0035] Optionally, when the station type corresponding to the boundary point is an intermediate station, the first set of operating equations corresponding to the boundary point is: ------Formula (1); in, and The first The inbound and outbound traffic of each intermediate station For the first The injected or distributed flow of an intermediate station; and The first The inbound and outbound pressure heads of each intermediate station; and For the first Characteristic curve parameters of the pump unit at each intermediate station; and The first The sum of valve resistance coefficients before and after the injection or distribution point of each intermediate station; , For the first The point within the second-to-last pipe segment upstream of the intermediate station Characteristic line coefficients; , For the first The second pipe section downstream of the intermediate station Characteristic line coefficients.

[0036] Specifically, in this embodiment, the solution can be obtained by solving a system of simultaneous equations. Solve the general formula (2), and then substitute it into the system of equations (1) to obtain the results. , and .in, The general formula (2) for solving this problem is: .

[0037] For details, see Figure 3 The diagram shown illustrates the boundary adjustment solution for the first station. Regarding the boundary conditions for the first station, compared to intermediate stations, the first station lacks an upstream pipe section, therefore it lacks... Characteristic lines, only exist Characteristic line. The valve is located after the first-station oil pump. The first-station has two control modes: constant pressure control and constant flow control. For constant flow control, the outflow rate... Known equal to the set injection flow rate Then, using Solving the characteristic line equation for the head of the first station exit Then, based on the valve and pump characteristic curves, the inlet (pump inlet) pressure head of the first station is calculated in reverse. .

[0038] Optionally, when the station type corresponding to the boundary point is the first station, the first set of operating equations corresponding to the boundary point is either the first station constant flow control equation set or the first station constant pressure control equation set, where: The first station's constant flow control equations are: ------Formula (3); in, For outbound traffic, For the set injection flow, For exiting the station, , The second inner point of the downstream section of the first station coefficients of characteristics For the pressure head entering the station, This is the sum of the valve resistance coefficients after the initial injection point. , The characteristic curve parameters of the pump unit at the first station; The first station's constant pressure control equation set is: ------Equation (4); This is the set value.

[0039] It should be noted that the solution equations for constant pressure control are similar to those for constant flow control, the difference being the pump inlet head under constant pressure control. Known and set value ,but It is unknown and needs to be solved according to equation (4). The solution is obtained by solving a system of simultaneous equations. Solve the general formula (5), and then substitute it into formula (4) to obtain .in, The general formula (5) for solving this problem is: ------Equation (5).

[0040] Specifically, such as Figure 4 The schematic diagram shown illustrates the boundary adjustment solution for the terminal station. Compared to intermediate stations, the terminal station lacks pumps and downstream pipelines, thus lacking... Characteristic lines, only exist Characteristic line. The valve is located before the distribution point. The terminal station also has two control modes: constant pressure control and constant flow control. For constant flow control, the inlet flow rate... The known and equal distribution flow rate Then, using Solving the characteristic line equation for the head of the terminal station's entry point Then, based on the valve and pump characteristic curves, the outlet (after the valve) pressure head of the terminal station is calculated. .

[0041] Optionally, when the station type corresponding to the boundary point is the terminal station, the first set of operating equations corresponding to the boundary point is the terminal station constant flow control equation set or the terminal station constant pressure control equation set. The control equations for the terminal station's constant flow are: ------Formula (6); in, For inbound traffic, For the set distribution flow, For the pressure head entering the station, , The second to last inner point of the upstream section of the terminal station coefficients of characteristics This is the sum of the valve resistance coefficients before the final distribution point. For exiting the station; The control equations for the terminal station's constant pressure are as follows: ------Equation (7); in, This is the set value.

[0042] The solution equations for constant pressure control are similar to those for constant flow control, the difference being the pump inlet head under constant pressure control. Known and set value ,but It is unknown and needs to be solved according to the system of equations (7). The solution is obtained by solving the system of equations. Solve the general formula (8), and then substitute it into the system of equations (7) to obtain .in, The general formula (8) for solving this problem is: ------Equation (8).

[0043] S104. Based on the operating characteristics of each point within a pipe segment, construct the second set of operating equations corresponding to each point within a pipe segment.

[0044] Specifically, in this embodiment, a second set of operating equations can be constructed for points within pipe segments in different connecting pipelines.

[0045] Optionally, the second set of operating equations corresponding to the points within the pipe segment is: ------Equation (9); or ------Formula (10); ------Formula (11); ------Equation (12); in, For the first section Flow rate at points within a pipe section; For the first section Pressure head at points within a pipe section; and For the first section Each internal point Characteristic line coefficients; and For the first Each internal point Characteristic line coefficients; , These are the pipe segments at the 1st step of the previous calculation. The, the Flow rate at points within a pipe section; and These are the pipe segments at the 1st step of the previous calculation. The, the The pressure head at each internal point; and For the pipe segment in the previous calculation step The, the Corrected wave velocity at points within each pipe section; It is the acceleration due to gravity; and This refers to the cross-sectional area of ​​the pipe. , For the pipe segment in the previous calculation step The, the The coefficient of friction at each internal point; The wave velocity of the pressure wave. For time step.

[0046] S105. Based on each first set of operating equations and each second set of operating equations, perform transient simulation of the liquid fuel pipeline.

[0047] Specifically, this embodiment can perform transient simulation of liquid fuel pipelines based on the first set of operating equations corresponding to each boundary point and the second set of operating equations corresponding to each point within the pipe segment.

[0048] Optionally, step S105 above includes: Based on the operating status of each station and each pipe segment, initialize the operating parameters in each first set of operating equations and each second set of operating equations to obtain the initial parameter set for each boundary point and each point within the pipe segment; Based on each first set of operating equations, each second set of operating equations, and the equipment status changes in each station, update the initial parameter set of each boundary point and each point within the pipe segment to obtain the new parameter set of each boundary point and each point within the pipe segment. Based on the new parameter sets for each first set of operating equations, each second set of operating equations, each boundary point, and each point within a pipe segment, determine and output the flow rate and head for each boundary point and each point within a pipe segment; Based on each first set of operating equations, each second set of operating equations, equipment status changes in each station, flow rate and head at each boundary point and point within each pipe segment, update the new parameter set for each boundary point and point within each pipe segment until the set iteration termination condition is met.

[0049] Specifically, this embodiment initializes the operating parameters in each first operating equation set and each second operating equation set based on the operating status of the first station, the last station, each intermediate station and each pipe segment in the liquid fuel pipeline system, so as to obtain the initialization parameter set for each boundary point and each point within the pipe segment.

[0050] It is understandable that the initial parameter group can include the initial parameter values ​​of multiple parameters, including the head, flow rate, friction coefficient, etc. of all points within the pipe section, and the head, flow rate, valve resistance coefficient, pump unit coefficient, etc. of boundary points.

[0051] Specifically, in this embodiment, after updating the new parameter sets for each boundary point and each point within a pipe segment, the latest parameter sets for each boundary point and each point within a pipe segment can be obtained. Then, based on each first operating equation set, each second operating equation set, and the latest parameter sets for each boundary point and each point within a pipe segment, the flow rate and pressure head for each boundary point and each point within a pipe segment are determined and output. Then, based on each first operating equation set, each second operating equation set, the equipment status changes in each station, and the flow rate and pressure head, the latest parameter sets for each boundary point and each point within a pipe segment are updated until the set iteration termination condition is met.

[0052] Specifically, in this embodiment, for any point within a pipe segment, the time step t is calculated according to formulas (11)-(12). and Characteristic line coefficients.

[0053] Specifically, determine whether there is operation at each station. If pumps are operating at a station, update the pump unit coefficients at the corresponding boundary points according to the operating conditions. and If valves are being operated at the station, the valve resistance coefficient at the corresponding boundary point will be updated according to the operating conditions. , If no operation is performed, the parameters will not be updated.

[0054] Specifically, in this embodiment, for all boundary points, the station type corresponding to the boundary point is determined. If it is the first station, the flow rate and head are updated according to formulas (3)-(5); if it is an intermediate station, the flow rate and head are updated according to formulas (1)-(2); if it is the last station, the flow rate and head are updated according to formulas (6)-(8). For all points within the pipe segment, the flow rate and head are first updated according to formulas (9)-(10), and then the flow regime is re-determined and the friction coefficient is updated.

[0055] This embodiment of the transient simulation process for liquid fuel pipelines enables self-adaptive transient simulation of liquid fuel pipelines under multiple boundary conditions, solving the challenge of dynamically changing boundary conditions under various operating conditions in transient simulation of liquid fuel pipelines. This embodiment can derive universal solution formulas for the boundary conditions of the topology configuration of components such as injection points, distribution points, pumps, and valves in the station, automatically matching and calculating the formula parameters without requiring model reconstruction, thus providing support for transient simulation of liquid fuel pipelines under complex multi-operating conditions.

[0056] This embodiment derives a universal solution formula for the boundary conditions of topology configurations. This formula can accurately simulate the complex transient processes of liquid fuel pipelines under long-distance, multi-pump, multi-distribution, and multi-operating-condition conditions. It can simultaneously handle various boundary conditions such as pump start-up and shutdown, flow regulation at distribution points and injection points, and allows for free combination during multi-condition switching. Through this adaptive solution formula for various boundary conditions, this embodiment ensures a high degree of matching between simulation results and actual operating conditions, significantly improving simulation accuracy and applicability, and providing reliable support for the safe operation of pipelines under multiple operating conditions.

[0057] The transient simulation method for liquid fuel pipelines proposed in this embodiment can model the liquid fuel pipeline system to obtain a liquid fuel pipeline model. This model includes multiple boundary points and connecting pipelines between adjacent boundary points. Each boundary point uniquely corresponds to the first, intermediate, or final station in the liquid fuel pipeline system, and each connecting pipeline uniquely corresponds to a pipe segment between adjacent stations. Each connecting pipeline is spatially discretized to determine at least one point within a pipe segment. A first set of operating equations is constructed for each boundary point based on its corresponding station type; and a second set of operating equations is constructed for each point within a pipe segment based on its operating characteristics. Transient simulation of the liquid fuel pipeline is then performed based on each first and second set of operating equations. This embodiment can improve the accuracy and applicability of transient simulation of liquid fuel pipelines, providing reliable support for the safe operation of pipelines under multiple operating conditions.

[0058] based on Figure 1This embodiment proposes a second transient simulation method for liquid fuel pipelines. Before initializing the operating parameters in each first and second set of operating equations based on the operating status of each station and each pipeline segment, this method further includes: In response to the simulation run command, the simulation time step timing starts from 0; When the simulation time step is 0, the above steps are then executed to initialize the operating parameters in each first set of operating equations and each second set of operating equations based on the operating status of each station and each pipe section.

[0059] Optionally, after updating the initial parameter sets for each boundary point and each point within a pipe segment to obtain new parameter sets for each boundary point and each point within a pipe segment, the method further includes: Increment the current simulation time step by 1 to obtain the new simulation time step; Determine whether the new simulation time step has reached the set maximum time step; If not, then it is determined that the iteration termination condition has not been met, and the above steps are performed to determine and output the flow rate and head of each boundary point and each point in the pipe segment based on the new parameter set of each first set of operating equations, each second set of operating equations, each boundary point and each point in the pipe segment; If so, then the iteration termination condition is met, and the above steps of determining and outputting the flow rate and head of each boundary point and each point in the pipe segment based on the new parameter set of each first set of operating equations, each second set of operating equations, each boundary point, and each point in the pipe segment are prohibited.

[0060] A case study analysis is conducted using an oil pipeline system. The system consists of 13 pipe segments, detailed in Table 1. The topology of the oil pipeline system can be found in [reference needed]. Figure 5 N1 is the injection point, and N2 and N3 are the distribution points. Basic data on their pipelines and nodes can be found in [reference needed]. Figure 6 .

[0061] This embodiment tested the system under seven different operating conditions, covering multiple control modes and typical operating scenarios such as distribution, shutdown, intermediate station pump slinging, and abnormal closure of intermediate shut-off valves. The results were compared with relevant software. Figure 7 As shown, the pressure and flow errors under various operating conditions are displayed in detail. The results show that the system performs well overall under various operating conditions.

[0062] like Figure 8The simulation results of this embodiment, using relevant software, show the upstream and downstream changes over time after pump slinging. The average absolute error of pressure is controlled within 0.05 MPa, and the average absolute error of flow rate is also within 2 m³ / h. Taking the fifth operating condition as an example, the solution results of upstream P3 and downstream N2 of disturbance point P4 are selected to analyze the transient process. Pump P4 is instantaneously shut off at the 60th minute, the flow rate decreases instantaneously, and a pressure boosting wave is transmitted upstream, while a pressure depressurization wave is transmitted downstream. The outlet pressure of P3 instantaneously increases from 8.16 MPa to 9.20 MPa, while the inlet pressure of N2 instantaneously decreases from 10.40 MPa to 8.47 MPa. As the water hammer wave propagates upstream and downstream, the system gradually returns to a stable state. The outlet pressure of P3 rises to 8.82 MPa, the inlet pressure of N2 decreases to 7.77 MPa, and the overall flow rate decreases from 892 m³ / h before pump slinging to 811 m³ / h. The trend of change is in perfect agreement with the SPS results, with the mean absolute error of pressure being 0.020 MPa and the mean absolute error of flow rate being 1.99 m³ / h.

[0063] Simulation results show that pipelines frequently switch between multiple operating conditions during daily operation. These changes directly trigger water hammer waves within the system, leading to drastic fluctuations in pressure and flow along the pipeline. This embodiment provides accurate transient simulation analysis, helping maintenance personnel to fully understand the operational status of the entire pipeline and providing crucial technical support for developing scientific control strategies and ensuring the long-term safe and stable operation of the pipeline.

[0064] like Figure 9 As shown, this embodiment proposes a transient simulation device for liquid fuel pipelines, including: Modeling unit 101 is used to model the liquid fuel pipeline system to obtain a liquid fuel pipeline model; wherein, the liquid fuel pipeline model includes multiple boundary points and connecting pipelines between adjacent boundary points, each boundary point uniquely corresponds to the first station, intermediate station or terminal station in the liquid fuel pipeline system, and each connecting pipeline uniquely corresponds to a pipe segment between adjacent stations; Spatial discretization unit 102 is used to spatially discretize each connecting pipeline to determine at least one point within a pipe segment of each connecting pipeline. The first construction unit 103 is used to construct the first set of running equations corresponding to each boundary point according to the station type corresponding to each boundary point; The second building unit 104 is used to build a second set of operating equations corresponding to each point in each pipe segment based on the operating characteristics of each point in the pipe segment. Simulation unit 105 is used to perform transient simulation of liquid fuel pipelines based on each first set of operating equations and each second set of operating equations.

[0065] It should be noted that the processing procedures and beneficial effects of the modeling unit 101, spatial discretization unit 102, first building unit 103, second building unit 104, and simulation unit 105 can be found in the following references: Figure 1 Steps S101 to S105 are not described in detail here.

[0066] Optionally, simulation unit 105 is also used for: Based on the operating status of each station and each pipe segment, initialize the operating parameters in each first set of operating equations and each second set of operating equations to obtain the initial parameter set for each boundary point and each point within the pipe segment; Based on each first set of operating equations, each second set of operating equations, and the equipment status changes in each station, update the initial parameter set of each boundary point and each point within the pipe segment to obtain the new parameter set of each boundary point and each point within the pipe segment. Based on the new parameter sets for each first set of operating equations, each second set of operating equations, each boundary point, and each point within a pipe segment, determine and output the flow rate and head for each boundary point and each point within a pipe segment; Based on each first set of operating equations, each second set of operating equations, equipment status changes in each station, flow rate and head at each boundary point and point within each pipe segment, update the new parameter set for each boundary point and point within each pipe segment until the set iteration termination condition is met.

[0067] Optionally, the above-mentioned device further includes: a first execution unit and a second execution unit; The first execution unit is configured to: before initializing the operating parameters in each first set of operating equations and each second set of operating equations according to the operating status of each station and each pipe segment, in response to the simulation operation command, start the simulation time step timing from 0; when the simulation time step is 0, execute the step of initializing the operating parameters in each first set of operating equations and each second set of operating equations according to the operating status of each station and each pipe segment; The second execution unit is used for: After updating the initial parameter set for each boundary point and each point within a pipe segment to obtain the new parameter set for each boundary point and each point within a pipe segment, increment the current simulation time step by 1 to obtain the new simulation time step; Determine whether the new simulation time step has reached the set maximum time step; If not, then determine that the iteration termination condition has not been met, and execute the steps of determining and outputting the flow rate and head for each boundary point and each point within the pipe segment based on the new parameter set for each first set of running equations, each second set of running equations, each boundary point, and each point within the pipe segment; If so, then the iteration termination condition is satisfied, and the step of determining and outputting the flow rate and head of each boundary point and each point in the pipe segment based on the new parameter set of each first set of operating equations, each second set of operating equations, each boundary point, and each point in the pipe segment is prohibited.

[0068] Optionally, when the station type corresponding to the boundary point is an intermediate station, the first set of operating equations corresponding to the boundary point is: ; in, and The first The inbound and outbound traffic of each intermediate station For the first The injected or distributed flow of an intermediate station; and The first The inbound and outbound pressure heads of each intermediate station; and For the first Characteristic curve parameters of the pump unit at each intermediate station; and The first The sum of valve resistance coefficients before and after the injection or distribution point of each intermediate station; , For the first The point within the second-to-last pipe segment upstream of the intermediate station Characteristic line coefficients; , For the first The second pipe section downstream of the intermediate station Characteristic line coefficients.

[0069] Optionally, when the station type corresponding to the boundary point is the first station, the first set of operating equations corresponding to the boundary point is either the first station constant flow control equation set or the first station constant pressure control equation set, where: The first station's constant flow control equations are: ; in, For outbound traffic, For the set injection flow, For exiting the station, , The second inner point of the downstream section of the first station coefficients of characteristics For the pressure head entering the station, This is the sum of the valve resistance coefficients after the initial injection point. , The characteristic curve parameters of the pump unit at the first station; The first station's constant pressure control equation set is: ; This is the set value.

[0070] Optionally, when the station type corresponding to the boundary point is the terminal station, the first set of operating equations corresponding to the boundary point is the terminal station constant flow control equation set or the terminal station constant pressure control equation set. The control equations for the terminal station's constant flow are: ; in, For inbound traffic, For the set distribution flow, For the pressure head entering the station, , The second to last inner point of the upstream section of the terminal station coefficients of characteristics This is the sum of the valve resistance coefficients before the final distribution point. For exiting the station; The control equations for the terminal station's constant pressure are as follows: ; in, This is the set value.

[0071] Optionally, the second set of operating equations corresponding to the points within the pipe segment is: ; or ; ; ; in, For the first section Flow rate at points within a pipe section; For the first section Pressure head at points within a pipe section; and For the first section Each internal point Characteristic line coefficients; and For the first Each internal point Characteristic line coefficients; , These are the pipe segments at the 1st step of the previous calculation. The, the Flow rate at points within a pipe section; and These are the pipe segments at the 1st step of the previous calculation. The, the The pressure head at each internal point; and For the pipe segment in the previous calculation step The, the Corrected wave velocity at points within each pipe section; It is the acceleration due to gravity; and This refers to the cross-sectional area of ​​the pipe. , For the pipe segment in the previous calculation step The, the The coefficient of friction at each internal point; The wave velocity of the pressure wave. For time step.

[0072] The transient simulation device for liquid fuel pipelines proposed in this embodiment can model a liquid fuel pipeline system to obtain a liquid fuel pipeline model. The liquid fuel pipeline model includes multiple boundary points and connecting pipelines between adjacent boundary points. Each boundary point uniquely corresponds to the first, intermediate, or final station in the liquid fuel pipeline system, and each connecting pipeline uniquely corresponds to a pipe segment between adjacent stations. Each connecting pipeline is spatially discretized to determine at least one point within a pipe segment of each connecting pipeline. A first set of operating equations is constructed for each boundary point based on its corresponding station type; and a second set of operating equations is constructed for each point within a pipe segment based on its operating characteristics. Transient simulation of the liquid fuel pipeline is performed based on each first and second set of operating equations. This embodiment can improve the accuracy and applicability of transient simulation of liquid fuel pipelines, providing reliable support for the safe operation of pipelines under multiple operating conditions.

[0073] In this embodiment, the transient simulation device for liquid fuel pipelines is presented in the form of functional units. Here, a unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and memory that execute one or more software or fixed programs, and / or other devices that can provide the above functions.

[0074] This invention also provides a computer device having the above-described features. Figure 9 The transient simulation device for liquid fuel pipelines shown is illustrated.

[0075] Please see Figure 10The present invention provides a schematic diagram of the structure of a computer device according to an optional embodiment. The computer device includes one or more processors 10, a memory 20, and interfaces for connecting the various components, including high-speed interfaces and low-speed interfaces. The various components are interconnected via different buses and can be mounted on a common motherboard or otherwise installed as needed. The processors can process instructions executed within the computer device, including instructions stored in or on memory to display graphical information of a GUI on an external input / output device (such as a display device coupled to the interface). In some optional embodiments, multiple processors and / or multiple buses can be used with multiple memories, if desired. Similarly, multiple computer devices can be connected, each providing some of the necessary operations (e.g., as a server array, a group of blade servers, or a multiprocessor system). Figure 10 Take a processor 10 as an example.

[0076] Processor 10 may be a central processing unit, a network processor, or a combination thereof. Processor 10 may further include a hardware chip. The hardware chip may be an application-specific integrated circuit (ASIC), a programmable logic device (PLD), or a combination thereof. The programmable logic device may be a complex programmable logic device (CAMP), a field-programmable gate array (FPGA), a general-purpose array logic (GDA), or any combination thereof.

[0077] The memory 20 stores instructions executable by at least one processor 10 to cause at least one processor 10 to perform the method shown in the above embodiments.

[0078] The memory 20 may include a program storage area and a data storage area. The program storage area may store the operating system and applications required for at least one function. The data storage area may store data created based on the use of the computer device. Furthermore, the memory 20 may include high-speed random access memory and may also include non-transitory memory, such as at least one disk storage device, flash memory device, or other non-transitory solid-state storage device. In some alternative embodiments, the memory 20 may optionally include memory remotely located relative to the processor 10, which can be connected to the computer device via a network. Examples of such networks include, but are not limited to, the Internet, intranets, local area networks, mobile communication networks, and combinations thereof.

[0079] Memory 20 may include volatile memory, such as random access memory. Memory may also include non-volatile memory, such as flash memory, hard disk, or solid-state drive. Memory 20 may also include combinations of the above types of memory.

[0080] The computer device also includes a communication interface 30 for communicating with other devices or communication networks.

[0081] This invention also provides a computer-readable storage medium. The methods described above according to embodiments of the invention can be implemented in hardware or firmware, or implemented as computer code that can be recorded on a storage medium, or implemented as computer code downloaded via a network and originally stored on a remote storage medium or a non-transitory machine-readable storage medium and then stored on a local storage medium. Thus, the methods described herein can be processed by software stored on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. The storage medium can be a magnetic disk, optical disk, read-only memory, random access memory, flash memory, hard disk, or solid-state drive, etc.; further, the storage medium can also include combinations of the above types of memory. It is understood that computers, processors, microprocessor controllers, or programmable hardware include storage components capable of storing or receiving software or computer code, which, when accessed and executed by the computer, processor, or hardware, implements the methods shown in the above embodiments.

[0082] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for transient simulation of a liquid fuel pipeline, characterized by, The method comprises the following steps: modeling a liquid fuel pipeline system to obtain a liquid fuel pipeline model, wherein the liquid fuel pipeline model comprises a plurality of boundary points and connecting pipelines between adjacent boundary points, each boundary point uniquely corresponds to a first station, an intermediate station or a terminal station in the liquid fuel pipeline system, and each connecting pipeline uniquely corresponds to a pipe section between adjacent stations; spatially discretizing each connecting pipeline to determine points in at least one pipe section in each connecting pipeline; constructing a first operation equation set corresponding to each boundary point according to the type of the station corresponding to each boundary point, and constructing a second operation equation set corresponding to each point in the pipe section according to the operating characteristics of each point in the pipe section; based on each first operation equation set and each second operation equation set, performing liquid fuel pipeline transient simulation.

2. The method of claim 1, wherein, The liquid fuel pipeline transient simulation based on each first operation equation set and each second operation equation set comprises: initializing the operating parameters in each first operation equation set and each second operation equation set according to the operating state of each station and each pipe section to obtain an initialization parameter set of each boundary point and each point in the pipe section; updating the initialization parameter set of each boundary point and each point in the pipe section based on the change of the device state in each first operation equation set, each second operation equation set and each station to obtain a new parameter set of each boundary point and each point in the pipe section; determining and outputting the flow and pressure head of each boundary point and each point in the pipe section according to the new parameter set of each first operation equation set, each second operation equation set, each boundary point and each point in the pipe section; updating the new parameter set of each boundary point and each point in the pipe section based on each first operation equation set, each second operation equation set, the change of the device state in each station, the flow and pressure head of each boundary point and each point in the pipe section until a set iteration end condition is met.

3. The method of claim 2, wherein, Before the step of initializing the operating parameters in each first operation equation set and each second operation equation set according to the operating state of each station and each pipe section, the method further comprises: starting simulation time step counting from 0 in response to a simulation operation instruction; when the simulation time step is 0, performing the step of initializing the operating parameters in each first operation equation set and each second operation equation set according to the operating state of each station and each pipe section; after the step of updating the initialization parameter set of each boundary point and each point in the pipe section to obtain a new parameter set of each boundary point and each point in the pipe section, the method further comprises: adding 1 to the current simulation time step to obtain a new simulation time step; determining whether the new simulation time step reaches a set maximum time step; and if the new simulation time step does not reach the set maximum time step, repeating the steps of updating the initialization parameter set of each boundary point and each point in the pipe section to obtain a new parameter set of each boundary point and each point in the pipe section, determining whether the new simulation time step reaches the set maximum time step, and adding 1 to the current simulation time step to obtain a new simulation time step until the new simulation time step reaches the set maximum time step. If not, it is determined that the iteration end condition is not met, and the step of determining and outputting the flow rate and pressure head of each boundary point and each point in the pipe segment according to the new parameter set of each first operating equation group, each second operating equation group, each boundary point and each point in the pipe segment is executed. If yes, it is determined that the iteration end condition is met, and the step of determining and outputting the flow rate and pressure head of each boundary point and each point in the pipe segment according to the new parameter set of each first operating equation group, each second operating equation group, each boundary point and each point in the pipe segment is prohibited.

4. The method of claim 3, wherein, When the station type corresponding to the boundary point is an intermediate station, the first operating equation group corresponding to the boundary point is: ； wherein, and are the inlet flow and outlet flow of the intermediate station, respectively, is the injection or distribution flow of the intermediate station; and are the inlet head and outlet head of the intermediate station, respectively; and are the pump set characteristic curve parameters of the intermediate station; and are the sum of the valve resistance coefficients before and after the injection or distribution point of the intermediate station, respectively; , is the characteristic line coefficient of the point in the penultimate pipe section of the upstream pipe section of the intermediate station; , is the characteristic line coefficient of the point in the second pipe section of the downstream pipe section of the intermediate station.

5. The method of claim 3, wherein, When the station type corresponding to the boundary point is a first station, the first operating equation group corresponding to the boundary point is a first station flow control equation group or a first station pressure control equation group, wherein: The first station flow control equation group is: ； wherein, is the outflow, is the set injection flow, is the outflow head, , is the first inner point of the downstream pipe section of the first station characteristic line coefficient, is the inflow head, is the sum of the valve resistance coefficients after the injection point of the first station, , is the pump unit characteristic curve parameter of the first station; The first station pressure control equation group is: ； is set to a value.

6. The method of claim 3, wherein, When the station type corresponding to the boundary point is a last station, the first operating equation group corresponding to the boundary point is a last station flow control equation group or a last station pressure control equation group; The last station flow control equation group is: ； wherein, is the inlet flow rate, is the set distribution flow rate, is the inlet pressure head, , is the last-but-one internal point of the upstream pipe section of the terminal station, characteristic line coefficient, is the sum of the valve resistance coefficients before the distribution point of the terminal station, is the outlet pressure head; The last station pressure control equation group is: ； wherein is a set value.

7. The method of claim 3, wherein, The second operating equation group corresponding to the point in the pipe segment is: ； or ; ； ； in, The first segment of the pipe Flow rate at points within a pipe section; The first segment of the pipe Pressure head at points within a pipe section; and The first segment of the pipe Each internal point Characteristic line coefficients; and For the first Each internal point Characteristic line coefficients; , These are the pipe segments at the 1st step of the previous calculation. The, the Flow rate at points within a pipe section; and These are the pipe segments at the 1st step of the previous calculation. The, the The pressure head at each internal point; and For the pipe segment in the previous calculation step The, the Corrected wave velocity at points within each pipe section; It is the acceleration due to gravity; and This refers to the cross-sectional area of ​​the pipe. , For the pipe segment in the previous calculation step The, the The coefficient of friction at each internal point; The wave velocity of the pressure wave. For time step.

8. A device for simulating transients in a liquid fuel pipeline, characterized in that Comprising: A modeling unit configured to model a liquid fuel pipeline system to obtain a liquid fuel pipeline model, wherein the liquid fuel pipeline model comprises a plurality of boundary points and connecting pipelines between adjacent boundary points, each boundary point uniquely corresponding to a first station, an intermediate station or a last station in the liquid fuel pipeline system, and each connecting pipeline uniquely corresponding to a pipe segment between adjacent stations; A spatial discretization unit configured to spatially discretize each connecting pipeline to determine at least one point in each connecting pipeline; A first construction unit configured to construct a first operating equation group corresponding to each boundary point according to the station type corresponding to each boundary point, respectively; A second construction unit configured to construct a second operating equation group corresponding to each point in the pipe segment according to the operating characteristics of each point in the pipe segment; A simulation unit configured to perform liquid fuel pipeline transient simulation based on each first operating equation group and each second operating equation group.

9. A computer device, comprising: Comprising: A memory and a processor, which are communicatively connected to each other, and the memory stores computer instructions, and the processor executes the computer instructions to perform the liquid fuel pipeline transient simulation method of any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer readable storage medium stores computer instructions for causing a computer to execute the liquid fuel pipeline transient simulation method of any one of claims 1 to 7.

Images