A multi-species liquid pipeline sequential transportation hydraulic simulation method and system based on dynamic interpolation

Abstract

The present application relates to the field of liquid fuel pipeline simulation, and discloses a kind of multi-variety liquid pipeline sequential transport hydraulic simulation method and system based on dynamic interpolation method, can be in the multiple division pipe sections of multi-variety liquid fuel pipeline sequential transport simulation model, select the two adjacent pipe sections for transporting different varieties of liquid fuel. According to the maximum fuel wave speed in the model, each adjacent pipe section is respectively grid divided, and at least one grid node in each adjacent pipe section is determined. The relative position relationship between the oil mixing interface and the grid node in the two adjacent pipe sections is monitored. According to the relative position relationship, the fuel wave speed of the grid node, the maximum fuel wave speed and the dynamic interpolation method, the flow and pressure head of at least one grid node are determined. The present application can accurately depict the multi-variety liquid fuel pipeline sequential transport process, realize the accurate simulation of sequential transport scene, and effectively enrich the multi-variety liquid pipeline sequential transport hydraulic simulation method.

Description

Technical Field

[0001] This invention relates to the field of liquid fuel pipeline simulation, and in particular to a hydraulic simulation method and system for sequential transportation of multiple types of liquid fuels in pipelines based on dynamic interpolation. Background Technology

[0002] With the advancement of the dual-carbon strategy, the diversification of energy consumption structure, and the deepening of the construction of a nationwide oil and gas pipeline network, the terminal market's demand for differentiated supply of multiple types of liquid fuels is becoming increasingly urgent, and a single medium transportation mode is difficult to adapt to the industry's development needs.

[0003] Sequential transport technology for multiple liquid fuels can significantly improve pipeline utilization and transport capacity by alternating different media in specific batches within the same pipeline, greatly reducing the huge investment required for laying new pipelines, while effectively reducing overall energy consumption and carbon emissions. It has become a core technological path to break through the bottleneck of energy storage and transportation.

[0004] However, the hydraulic characteristics of sequential transportation of multiple types of liquid fuels are fundamentally different from those of single-medium transportation, and there are few hydraulic simulation methods for sequential transportation of multiple types of liquid fuels in pipelines. Summary of the Invention

[0005] This invention provides a hydraulic simulation method and system for sequential transport of multi-variety liquids through pipelines based on dynamic interpolation, which addresses the deficiency of limited hydraulic simulation methods for sequential transport of multi-variety liquids in related technologies and effectively enriches the hydraulic simulation methods for sequential transport of multi-variety liquids through pipelines.

[0006] In a first aspect, the present invention provides a hydraulic simulation method for sequential pipeline transport of multiple types of liquids based on dynamic interpolation, comprising:

[0007] In the simulation model of sequential pipeline transportation of multiple types of liquid fuels, two adjacent pipeline segments are selected from the multiple segments to transport different types of liquid fuels. Based on the maximum fuel wave velocity in the model, each adjacent pipe segment is divided into grids, and at least one grid node in each adjacent pipe segment is determined. Monitor the relative positional relationship between the oil mixing interface in the two adjacent pipe sections and the grid nodes; Based on the relative positional relationship, the fuel wave velocity of the grid node, the maximum fuel wave velocity, and dynamic interpolation, determine the flow rate and pressure head of at least one of the grid nodes.

[0008] Optionally, determining the flow rate and pressure head of at least one of the grid nodes based on the relative positional relationship, the fuel wave velocity of the grid node, the maximum fuel wave velocity, and dynamic interpolation includes: When the oil mixing interface is located at the position corresponding to the grid node, the flow rate and pressure head of at least one of the grid nodes are determined based on the fuel wave velocity of the grid node, the maximum fuel wave velocity, and the time interpolation method. When the oil mixing interface is located between two adjacent grid nodes, the flow rate and pressure head of at least one grid node are determined based on the fuel wave velocity of the grid node, the maximum fuel wave velocity, and spatial interpolation.

[0009] Optionally, when the mixing interface is located at the position corresponding to the grid node, determining the flow rate and head of at least one of the grid nodes based on the fuel wave velocity of the grid node, the maximum fuel wave velocity, and the time interpolation method includes: When the oil mixing interface is located at the position corresponding to the grid node, for any one of the grid nodes: If the fuel wave velocity of the grid node at the current time step is equal to the maximum fuel wave velocity, then the initial flow rate and initial pressure head of the grid node in the next time step are determined based on the flow rate and pressure head of the upstream and downstream grid nodes at the current time step, and are directly used as the final flow rate and final pressure head of the grid node in the next time step. If the fuel wave velocity of the grid node at the current time step is less than the maximum fuel wave velocity, then based on the flow rate and pressure head of the upstream and downstream grid nodes at the current time step, the initial flow rate and initial pressure head of the grid node at the next time step are determined. The flow rate and pressure head of the grid node at the current time step, and the initial flow rate and initial pressure head of the grid node at the next time step are then input into the first model created according to the time interpolation method for interpolation to obtain the final flow rate and final pressure head of the grid node at the next time step.

[0010] Optionally, the first model is: ; ; ; or ; ; ; in, and For the upstream grid node At the current time step of Characteristic line coefficients; For the downstream grid node At the current time step of Characteristic line coefficients; and For the upstream grid node At the current time step Flow rate and pressure head; and The upstream grid nodes are respectively At the current time step The inertial water hammer coefficient and friction coefficient, and The downstream grid nodes are respectively At the current time step The inertial water hammer coefficient and friction coefficient; To calculate the distance step; and For the grid node In the next step The final flow rate and final head; For the grid node At the current time step, the wave speed The maximum wave speed; and For the grid node At the current time step Flow rate and pressure head; and For the grid node In the next step The initial flow rate and initial pressure head.

[0011] Optionally, when the oil mixing interface is located between two adjacent grid nodes, determining the flow rate and head of at least one grid node based on the fuel wave velocity of the grid node, the maximum fuel wave velocity, and spatial interpolation includes: When the mixing interface is located between the grid node and the upstream grid node, the fuel average wave velocity of the mixing grid is determined, the flow rate and pressure head of the grid node at the current time step are determined, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are determined. The mixing grid is the grid node. The average fuel wave velocity, the maximum fuel wave velocity, the flow rate and pressure head of the grid node at the current time step, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are input into the second model created according to the spatial interpolation method for interpolation to obtain the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step. Based on the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step, the flow rate and pressure head of the grid node at the next time step are determined.

[0012] Optionally, the second model is: ; ; ; ; ) ; ) ; ; in, For the mixed oil grid at the current time step The average wave velocity of fuel, For the mixed oil interface at the current time step Time and the grid node The distance; The step size is obtained by multiplying the maximum fuel wave velocity by the duration. and For the upstream grid node and the grid node at the current time step fuel wave speed; and For the mixed oil grid at the current time step The average friction coefficient and the average inertial water hammer coefficient, and For the upstream grid node at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the grid node at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the upstream grid node at the current time step The final flow rate and the final pressure head, and For the downstream grid node at the current time step The final flow rate and the final pressure head, and The grid node and the upstream grid node are respectively at the current time step. Initial traffic, and The grid node and the upstream grid node are respectively at the current time step. The initial pressure head, This represents the pressure head difference caused by the density difference at the interface of the oil mixture.

[0013] Optionally, when the oil mixing interface is located between two adjacent grid nodes, determining the flow rate and head of at least one grid node based on the fuel wave velocity of the grid node, the maximum fuel wave velocity, and spatial interpolation includes: When the oil mixing interface is located at the grid node and downstream grid nodes When determining the position between the nodes, the average fuel wave velocity of the mixed oil grid is determined, the flow rate and head of the grid node at the current time step are determined, and the initial flow rate and initial head of the upstream and downstream grid nodes at the current time step are determined, wherein the mixed oil grid is the grid node. With the downstream mesh node The average fuel wave velocity, the maximum fuel wave velocity, the flow rate and pressure head of the grid node at the current time step, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are input into a third model created according to spatial interpolation for interpolation to obtain the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step. Based on the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step, the flow rate and pressure head of the grid node at the next time step are determined.

[0014] Optionally, the third model is: ; ; ; ; ) ; ) ; ; ; ; in, For the mixed oil grid at the current time step The average wave velocity of fuel, For the mixed oil interface at the current time step Time and the downstream grid node The distance; The step size is obtained by multiplying the maximum fuel wave velocity by the duration. and For the downstream mesh node and the mesh node at the current time step fuel wave speed; and For the mixed oil grid at the current time step The average friction coefficient and the average inertial water hammer coefficient, and For the grid node At the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the downstream grid node at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the upstream grid node at the current time step The final flow rate and the final pressure head, and For the downstream grid node at the current time step The final flow rate and the final pressure head, and The grid node and the upstream grid node are respectively at the current time step. Initial traffic, and The grid node and the upstream grid node are respectively at the current time step. The initial pressure head, This represents the pressure head difference caused by the density difference at the interface of the oil mixture.

[0015] Optionally, after determining the flow rate and head of at least one of the grid nodes based on the relative positional relationship, the fuel wave velocity of the grid node, the maximum fuel wave velocity, and dynamic interpolation, the method further includes: Monitor operational changes at the first, intermediate, and final stations, and update the flow and pressure at the boundary points in the model based on the monitored operational changes at the first, intermediate, and final stations.

[0016] Secondly, the present invention provides a hydraulic simulation system for sequential pipeline transport of multiple types of liquids based on dynamic interpolation, comprising: Selecting a unit is used to select two adjacent pipe segments for transporting different types of liquid fuels from multiple divided pipe segments in a simulation model of sequential pipeline transportation of multiple types of liquid fuels. A partitioning unit is used to partition each of the adjacent pipe segments into grids according to the maximum fuel wave velocity in the model, and to determine at least one grid node in each of the adjacent pipe segments. The monitoring unit is used to monitor the relative positional relationship between the oil mixing interface in the two adjacent pipe sections and the grid node; The determining unit is used to determine the flow rate and pressure head of at least one of the grid nodes based on the relative positional relationship, the fuel wave velocity of the grid node, the maximum fuel wave velocity, and dynamic interpolation.

[0017] 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 hydraulic simulation method for sequential transportation of multi-variety liquids in pipelines based on dynamic interpolation as described in the first aspect or any corresponding embodiment.

[0018] Fourthly, the present invention provides a computer-readable storage medium storing computer instructions for causing a computer to execute the hydraulic simulation method for sequential transport of multi-variety liquids in pipelines based on dynamic interpolation, as described in the first aspect above or any corresponding embodiment.

[0019] This invention provides a hydraulic simulation method and system for sequential transport of multiple types of liquid fuels in pipelines based on dynamic interpolation. It allows selection of two adjacent pipe segments from multiple subdivided pipe segments in a simulation model of sequential transport of multiple types of liquid fuels, for transporting different types of liquid fuels. Based on the maximum fuel wave velocity in the model, each adjacent pipe segment is meshed, and at least one mesh node is determined within each adjacent segment. The relative positional relationship between the mixing interface in the two adjacent pipe segments and the mesh nodes is monitored. Based on the relative positional relationship, the fuel wave velocity of the mesh nodes, the maximum fuel wave velocity, and the dynamic interpolation method, the flow rate and head of at least one mesh node are determined. This invention enables hydraulic simulation of sequential transport of multiple types of liquid fuels in pipelines based on dynamic interpolation, accurately depicting the sequential transport process, achieving precise simulation of sequential transport scenarios, and effectively enriching the hydraulic simulation methods for sequential transport of multiple types of liquid fuels in pipelines.

[0020] This invention addresses the challenges of hydraulic simulation in pipelines transporting multiple types of liquid fuels sequentially, where the distribution of oil within the pipeline changes constantly, and the varying water hammer velocities due to differences in density and bulk modulus among different fuels increase the complexity of hydraulic simulation. It develops a dynamic interpolation method in both time and space, and combines it with the method of characteristics to create a hydraulic simulation approach suitable for pipelines transporting multiple types of liquid fuels sequentially. This method adapts to a unified time step and distance step partitioning mode, eliminating the need for complex mesh adjustments or time step switching. While improving the simulation accuracy of parameters such as pressure and flow rate, it also ensures the stability of the calculation process, effectively addressing the complex hydraulic simulation requirements of sequential transport of multiple types of liquid fuels. It provides accurate transient simulation analysis of sequential liquid fuel transport, helping maintenance personnel to fully understand the operational status of the entire pipeline and providing key technical support for formulating scientific control strategies and ensuring the long-term safe and stable operation of the pipeline. Attached Figure Description

[0021] 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.

[0022] Figure 1 A flowchart of a hydraulic simulation method for sequential pipeline transportation of multiple liquids based on dynamic interpolation, provided in an embodiment of the present invention; Figure 2 A schematic diagram of a grid distribution based on a pure oil segment grid time interpolation method is provided for an embodiment of the present invention; Figure 3 A schematic diagram of a grid distribution based on the spatial interpolation method on the left side of the oil mixing interface, provided for an embodiment of the present invention; Figure 4 A schematic diagram of a grid distribution based on the right-side spatial interpolation method of the oil mixing interface is provided for an embodiment of the present invention; Figure 5 This is a schematic diagram of a refined oil pipeline system provided in an embodiment of the present invention; Figure 6 A schematic diagram of pressure and flow rate errors under different operating conditions is provided for an embodiment of the present invention; Figure 7 A flow / pressure comparison chart for the entire line provided in an embodiment of the present invention; Figure 8 This invention provides a time-varying trend diagram of upstream and downstream changes after pumping stops at N3 station; Figure 9 A schematic diagram of a hydraulic simulation system for sequential transport of multiple liquids in pipelines based on dynamic interpolation, 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

[0023] 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.

[0024] The following is combined with Figures 1-8 This invention describes a hydraulic simulation method for sequential pipeline transport of multiple liquids based on dynamic interpolation.

[0025] like Figure 1 As shown, this embodiment proposes a first hydraulic simulation method for sequential pipeline transportation of multiple liquids based on dynamic interpolation. This method may include the following steps: S101. In the simulation model of sequential pipeline transportation of multiple types of liquid fuels, select two adjacent pipeline segments for transporting different types of liquid fuels from multiple pipe segments.

[0026] Specifically, this embodiment can construct a corresponding simulation model for the sequential transportation of multiple types of liquid fuels based on the actual parameters of the pipelines. Then, in this embodiment, the model is divided into several pipe segments according to the boundary elements such as pumps, valves, injection points, and distribution points, thus obtaining multiple segmented pipe sections.

[0027] Specifically, in this embodiment, based on the location distribution of each type of liquid fuel in the model at the current time step, two adjacent pipe segments for transporting two different liquid fuels can be selected from all the divided pipe segments of the model.

[0028] S102. Based on the maximum fuel wave velocity in the model, divide each adjacent pipe segment into grids and determine at least one grid node in each adjacent pipe segment.

[0029] Specifically, this embodiment determines the maximum fuel wave velocity among the wave velocities of various types of liquid fuels in the model. Based on the maximum fuel wave velocity, meshing is performed on two selected adjacent pipe segments to divide each segment into multiple meshes, with the boundaries between adjacent meshes serving as mesh nodes. It should be noted that this embodiment, by meshing the pipe segments based on the maximum fuel wave velocity, ensures adaptation to the flow dynamics corresponding to the fastest wave velocity.

[0030] Optionally, in this embodiment, when dividing the pipe segment into grids based on the maximum fuel wave velocity, if the number of grids is not an integer, the number of grids can be rounded down, and the fuel wave velocity of the pipe segment can be corrected based on the rounded number of grids. This ensures that the grid division both fits the actual calculation requirements and reflects the characteristics of the oil.

[0031] S103. Monitor the relative positional relationship between the oil mixing interface and the grid nodes in two adjacent pipe sections.

[0032] Specifically, in this embodiment, during model operation, the oil mixing interface in the two adjacent pipe segments is monitored, and the relative positional relationship between the oil mixing interface and the nearby grid nodes is also monitored.

[0033] S104. Based on the relative positional relationship, the fuel wave velocity of the grid node, the maximum fuel wave velocity, and the dynamic interpolation method, determine the flow rate and pressure head of at least one grid node.

[0034] Specifically, in this embodiment, the flow rate and pressure head of at least one grid node can be determined based on the monitored relative positional relationship, maximum fuel wave velocity, dynamic interpolation, and the fuel wave velocity of the grid node close to the oil mixing interface.

[0035] Optionally, step S104 includes: S1041. When the oil mixing interface is located at the position corresponding to the grid node, determine the flow rate and pressure head of at least one grid node based on the fuel wave velocity, maximum fuel wave velocity and time interpolation of the grid node. S1042. When the oil mixing interface is located between two adjacent grid nodes, determine the flow rate and pressure head of at least one grid node based on the fuel wave velocity, maximum fuel wave velocity, and spatial interpolation of the grid node.

[0036] Optionally, step S1041 includes: When the mixing interface is located at the position corresponding to the grid node, for any grid node: If the fuel wave velocity of a grid node at the current time step is equal to the maximum fuel wave velocity, then the initial flow rate and initial pressure head of the grid node in the next time step are determined based on the flow rate and pressure head of the upstream and downstream grid nodes at the current time step, and are directly used as the final flow rate and final pressure head of the grid node in the next time step. If the fuel wave velocity of a grid node at the current time step is less than the maximum fuel wave velocity, then the initial flow rate and initial pressure head of the grid node in the next time step are determined based on the flow rate and pressure head of the upstream and downstream grid nodes at the current time step. The flow rate and pressure head of the grid node at the current time step, and the initial flow rate and initial pressure head of the grid node in the next time step are then input into the first model created according to the time interpolation method for interpolation to obtain the final flow rate and final pressure head of the grid node in the next time step.

[0037] Optionally, the first model is: ----------Formula (1); ----------Formula (2); ----------Formula (3); or ----------Formula (4); ----------Formula (5); ----------Formula (6); in, and upstream grid node At the current time step of Characteristic line coefficients; For downstream grid nodes At the current time step of Characteristic line coefficients; and upstream grid node At the current time step Flow rate and pressure head; and Upstream grid nodes At the current time step The inertial water hammer coefficient and friction coefficient, and Downstream grid nodes At the current time step The inertial water hammer coefficient and friction coefficient; To calculate the distance step; and For grid nodes In the next step The final flow rate and final head; For grid nodes At the current time step, the wave speed Maximum wave speed; and For grid nodes At the current time step Flow rate and pressure head; and For grid nodes In the next step The initial flow rate and initial pressure head.

[0038] See Figure 2 This demonstrates the sequential delivery of two fuels with different wave velocities through a pipeline, with the mixing interface currently at a grid node. Above. Node Left side liquid fuel wave speed Greater than the wave speed on the right The pipe mesh is divided based on the wave velocity on the left. As the reference, the step size is set to For the left side filled with wave velocity... fuel pipeline inside point In terms of its upstream and downstream internal point water hammer waves, After reaching the interior point The right side is filled with wave speeds of fuel pipeline inside point In terms of its upstream and downstream internal point water hammer waves, After reaching the interior point , .

[0039] The numerical solution of the water hammer characteristic equation is applicable to points inside a pipe without boundary points, requiring that the medium inside the pipe must be a single fluid. In specific calculations, this embodiment can first calculate the mesh nodes using formulas (1)-(2). and Characteristic line coefficients, grid nodes Traffic With water head From the upstream grid node of the previous time step and downstream grid nodes The calculation is shown in formulas (3)-(4). However, for the right-side grid nodes, the result obtained by solving using the method of characteristics is actually the flow rate at the gray point. and pressure head Therefore, time interpolation is required to convert it into time steps. Point within a time Corresponding traffic and pressure head This ensures that the calculation results of all grid nodes are interpolated to the same time step, guaranteeing the time consistency of the simulation calculation of the entire pipeline system. This embodiment introduces formulas (5)-(6) to adapt to the calculation of points within pure oil pipe sections with different wave velocities. When the grid nodes... Wave speed Equal to the reference wave speed hour, , That is, no insertion is required within a certain time; when hour, Depend on , Linear weighted calculation, Depend on , Linear weighted calculation.

[0040] Optionally, step S1042 above includes: When the mixing interface is located between a grid node and an upstream grid node, determine the fuel-average wave velocity of the mixing grid, the flow rate and head of the grid node at the current time step, and the initial flow rate and initial head of the upstream and downstream grid nodes at the current time step. The mixing grid is defined as the grid node. For the grid between upstream and downstream grid nodes, the average fuel wave velocity, maximum fuel wave velocity, flow rate and pressure head of the grid node at the current time step, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are input into the second model created by spatial interpolation for interpolation to obtain the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step. Based on the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step, the flow rate and pressure head of the grid node at the next time step are determined.

[0041] Alternatively, the second model is: ----------Formula (7); ----------Formula (8); ----------Formula (9); ----------Formula (10; ) ----------Formula (11); ) ----------Formula (12); ----------Formula (13); in, For the mixed oil grid at the current time step The average wave velocity of fuel, For the mixed oil interface at the current time step Time and grid nodes The distance; The step size is obtained by multiplying the maximum fuel wave velocity by the duration. and For upstream grid nodes and grid nodes at the current time step fuel wave speed; and For the mixed oil grid at the current time step The average friction coefficient and the average inertial water hammer coefficient, and For the upstream grid node at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the grid node at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the upstream grid node at the current time step The final flow rate and the final pressure head, and For downstream grid nodes at the current time step The final flow rate and the final pressure head, and These represent the current time steps of the grid node and the upstream grid node, respectively. Initial traffic, and These represent the current time steps of the grid node and the upstream grid node, respectively. The initial pressure head, This represents the pressure head difference caused by the density difference at the oil mixing interface.

[0042] It should be noted that during sequential transport, the oil mixing interface inside the pipe will continuously and dynamically move with the flow of oil, and cannot always completely coincide with the pre-divided discrete grid nodes. Therefore, the interface will frequently fall between two grid nodes.

[0043] See Figure 3 The oil-mixing interface is in the grid ( Between ), assuming the mixed oil interface is at the current time step. Distance to right grid node distance Mixed oil mesh ( fuel average wave speed in ) Average friction coefficient Average inertial water hammer coefficient Calculated according to the oil volume ratio within the grid using equations (7)-(9). At that time, the oil mixing interface is at a grid node. Above, fuel properties between mixed oil grids and grid nodes The fuel properties are consistent. When At this time, the average fuel properties in the grid node are jointly determined by the location of the mixing interface and the fuel properties on both sides of the mixing interface. Grid node The mixed-oil wave velocity on the left side is necessarily less than the reference wave velocity. Grid nodes The water shock wave will be in time step Then reach the grid node .because The wave velocities of the grids on the left and right sides of the point are inconsistent, before interpolation , The time-step difference of the characteristic lines makes it impossible to calculate using the characteristic equation. Therefore, spatial interpolation is needed to convert the mesh nodes. Interpolate to the gray dot, making it The feature lines intersect at the time step. The interpolated distance is determined by the fuel-average wave velocity of the mixed-oil grid. Compared with the reference wave speed The ratio is determined as shown in formulas (10)-(11). In formula (11), This represents the pressure head difference caused by the density difference at the oil-mixing interface. This is because of the mesh nodes. and The pressure head difference originates from the pressure head difference caused by frictional resistance along the path and density difference. Similarly, the right-side grid nodes... The wave velocity of that grid node is then used. Compared with the reference wave speed The ratio determines the spatial interpolation, see formulas (12)-(13). If , , That is, no interpolation is required; if , By respectively and , and Linear weighted calculation.

[0044] Optionally, the second model described above may also include: ----------Formula (14); ----------Formula (15); ----------Formula (16).

[0045] Specifically, formula (14) calculates the mesh nodes. Inserted to the right The characteristic line coefficients, and formula (15) are used to calculate the mesh nodes. Inserted to the left Characteristic line coefficients. Finally, substitute the characteristic line coefficients into formula (3) to calculate the mesh nodes. In time step Traffic The grid nodes are calculated from the right side using formula (16). pressure head .

[0046] Optionally, step S1042 above includes: When the oil mixing interface is located at a grid node and downstream grid nodes When determining the position between the nodes, the fuel-average wave velocity of the mixed oil grid is determined, the flow rate and pressure head of the grid nodes at the current time step are determined, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are determined. The mixed oil grid is the grid node. With downstream grid nodes The average fuel wave velocity, maximum fuel wave velocity, flow rate and pressure head of the grid nodes at the current time step, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are input into the third model created by spatial interpolation for interpolation to obtain the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step. Based on the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step, the flow rate and pressure head of the grid nodes at the next time step are determined.

[0047] Alternatively, the third model is: ----------Formula (17); ----------Formula (18); ----------Formula (19); ----------Formula (20); ) ----------Formula (21); ) ----------Formula (22); ----------Formula (23); ----------Formula (24); ----------Formula (25); in, For the mixed oil grid at the current time step The average wave velocity of fuel, For the mixed oil interface at the current time step Time and downstream grid nodes The distance; The step size is obtained by multiplying the maximum fuel wave velocity by the duration. and For downstream grid nodes and grid nodes at the current time step fuel wave speed; and For the mixed oil grid at the current time step The average friction coefficient and the average inertial water hammer coefficient, and For grid nodes At the current time step The coefficient of friction and the coefficient of inertial water hammer, and For downstream grid nodes at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the upstream grid node at the current time step The final flow rate and the final pressure head, and For downstream grid nodes at the current time step The final flow rate and the final pressure head, and These represent the current time steps of the grid node and the upstream grid node, respectively. Initial traffic, and These represent the current time steps of the grid node and the upstream grid node, respectively. The initial pressure head, This represents the pressure head difference caused by the density difference at the oil mixing interface.

[0048] See Figure 4 This demonstrates that the mixed oil interface is in a grid ( In scenarios involving mixing oil, interpolation is performed on the right-side grid nodes of the mixing oil interface. We still assume the mixing oil interface is at a distance from the right-side grid node... distance Mixed oil mesh ( fuel average wave speed in ) Average friction coefficient Average inertial water hammer coefficient The calculation is still performed using equations (17)-(19). At this point, the mesh nodes... The velocity of the mixed oil wave on the right side must be less than the reference wave velocity. Therefore, it is necessary to set the grid nodes. Interpolate to the gray dot, making it The feature lines intersect at the time step. See formulas (20)-(21). Similarly, the left inner point The wave velocity of that grid node is then used. Compared with the reference wave speed The ratio is determined by formulas (22)-(23). Formula (24) is used to calculate the mesh nodes. Inserted to the right The characteristic line coefficients, and formula (25) are used to calculate the mesh nodes. Inserted to the left Characteristic line coefficients. Finally, substitute the characteristic line coefficients into formula (3) to calculate the mesh nodes. In time step Traffic The grid nodes are calculated from the left side using formula (26). pressure head Formula (26) is: ----------Formula (26).

[0049] It is understood that, in this embodiment, for any two adjacent pipe segments in the multiple divided pipe segments of the simulation model for sequential transportation of different types of liquid fuels, the following approach can be adopted: Figure 1 The steps shown are for determining the flow rate and head of the grid nodes.

[0050] It should be noted that the hydraulic characteristics of sequential transportation of multiple types of liquid fuels differ fundamentally from those of single-medium transportation, posing significant challenges to system design and operational control. First, the density, viscosity, and other physical properties of different media vary significantly, causing dynamic changes in pipeline friction coefficients and pump station characteristic curves during sequential transportation. This results in highly nonlinear characteristics in overall pressure distribution, flow stability, and energy supply and demand relationships. Second, oil mixing due to convection and diffusion inevitably occurs at the interfaces of adjacent media. Simulation methods using this approach, which assume the mixing interface lies on a defined grid, oversimplify the mesh treatment for areas with mixing interfaces, leading to substantial errors in sequential transportation scenarios.

[0051] This embodiment addresses the challenges of sequential pipeline transport of multiple types of liquid fuels, where the distribution of oil within the pipeline constantly changes, and the varying density and bulk modulus of different fuels lead to different water hammer velocities, increasing the complexity of hydraulic simulation. To address this, a dynamic interpolation method in both time and space was developed, combined with the method of characteristics, to create a hydraulic simulation approach suitable for pipelines transporting multiple types of liquid fuels sequentially. This embodiment adapts to a unified time step and distance step partitioning mode, eliminating the need for complex mesh adjustments or time step switching. While improving the simulation accuracy of parameters such as pressure and flow rate, it also ensures the stability of the calculation process, ultimately effectively meeting the complex hydraulic simulation requirements for the sequential transport of multiple types of liquid fuels.

[0052] This embodiment proposes a hydraulic simulation method for sequential transport of multiple types of liquid fuels in pipelines based on dynamic interpolation. From multiple segmented pipe sections in a simulation model of sequential transport of multiple types of liquid fuels, two adjacent pipe sections are selected for transporting different types of liquid fuels. Based on the maximum fuel wave velocity in the model, each adjacent pipe section is meshed, and at least one grid node is determined in each adjacent pipe section. The relative positional relationship between the oil mixing interface in the two adjacent pipe sections and the grid nodes is monitored. Based on the relative positional relationship, the fuel wave velocity of the grid nodes, the maximum fuel wave velocity, and the dynamic interpolation method, the flow rate and head of at least one grid node are determined. This embodiment can perform hydraulic simulation of sequential transport of multiple types of liquid fuels in pipelines based on dynamic interpolation, accurately depicting the sequential transport process of multiple types of liquid fuels, achieving precise simulation of sequential transport scenarios, and effectively enriching the hydraulic simulation methods for sequential transport of multiple types of liquid fuels in pipelines.

[0053] based on Figure 1 In the second hydraulic simulation method for sequential transport of multi-variety liquids in pipelines based on dynamic interpolation proposed in this embodiment, after step S104, the method further includes: Monitor operational changes at the first, intermediate, and last stations, and update the flow and pressure at boundary points in the model based on these changes.

[0054] It should be noted that the first station of a sequential transport pipeline generally adopts constant flow control, with the flow control value being... Then the outlet pressure can be measured at the point inside the first pipeline. The characteristic line is calculated as shown in equation (27). (Terminal station) Constant pressure control is adopted, and the pressure control value is... Connect the last grid node The characteristic lines and the valve characteristic equations of the terminal station yielded formula (28), which further led to the derivation of the terminal station flow rate. The general formula for solving this problem is (29). Intermediate boundary points are generally for flow control; assuming boundary points... Operational traffic is (Positive values ​​represent injection, negative values ​​represent distribution). Formula (30) is obtained by simultaneously solving the flow balance equations at the intermediate boundary point, the characteristic line equations at the upstream and downstream inner points, and the pump characteristic equations of this station. The general formula (31) for solving the inflow rate at the intermediate boundary point is further derived. Formulas (27) to (31) are shown below: ----------Formula (27); ----------Formula (28); ----------Formula (29); ----------Formula (30); in, The calculation formula (31) is as follows: .

[0055] in, for Time Step Flow control values ​​at each boundary point for Time Step Inbound / outbound traffic at each boundary point for Pressure control value at the end of the time step. for Time-step terminal valve characteristic parameters for Time Step Characteristic parameters of boundary point pump units for Time Step The inbound / outbound pressure heads at each boundary point. Other symbols are the same as above.

[0056] Based on the above, this embodiment proposes a third hydraulic simulation method for sequential pipeline transportation of multiple liquids based on dynamic interpolation. This method tracks the movement of the mixed oil interface in real time and achieves automatic spatiotemporal interpolation of the grid nodes by updating the liquid fuel property parameters of the grid nodes. This avoids complex grid adjustments and time-step switching, improving the simulation accuracy of parameters such as pressure and flow rate while ensuring the stability of the calculation process. Specifically, it includes the following steps: Step 1: Pipe Modeling and Meshing A simulation model was constructed based on the actual parameters of the pipeline. The maximum wave velocity of the sequentially transported liquid fuel was used as the benchmark for mesh division, and the wave velocity was corrected in reverse according to the mesh rounding.

[0057] Step 2, Update pipeline point parameters: Based on the location and physical characteristics of the liquid fuel batch, determine the relevant physical property parameters of the liquid fuel at each pipeline point and in the mixing grid, including wave velocity, friction coefficient, inertial water hammer coefficient, etc.

[0058] Step 3, Boundary point parameter update: Determine the first... At any given time, check if there are any operational changes at each station, and update the pressure / flow control values, pump unit / valve characteristic parameters, etc., at each station as needed.

[0059] Step 4, Update the flow rate and pressure at points within the pipeline: For all points within the pure oil pipeline section, perform time interpolation according to formulas (1)-(6) and calculate the first... The flow rate and pressure head at each step; for the inner point on the left side of the mixing interface, spatial interpolation is performed according to formulas (3) and (7)-(16) respectively, and the first step is calculated. The flow rate and pressure head at each step; for the inner point on the right side of the mixing interface, spatial interpolation is performed according to formulas (3) and (17)-(26) respectively, and the first step is calculated. The flow rate and pressure head of the time step.

[0060] Step 5, Boundary Point Flow and Pressure Update: For the first station, calculate the flow and pressure of the boundary point according to formula (27). The flow rate and pressure head at each time step; for the terminal station, calculate the flow rate and pressure head according to formulas (28)-(29) respectively. The flow rate and pressure head at each time step; for intermediate stations, the first step is calculated according to formulas (30)-(31). The flow rate and pressure head of the time step.

[0061] Step 6: Batch Interface Update: Calculate the flow velocity based on the flow rate of the pipe segment in the batch interface, and update the batch interface. The position of the batch interface during the time step is the distance between each grid node on the batch interface.

[0062] Step 7, Time Update: Determine if the maximum number of calculation steps has been reached. If so, terminate the calculation; otherwise, jump to Step 2.

[0063] To better explain the above content, such as Figure 5 As shown, the inventors of this invention selected an actual operating refined oil pipeline as the test object. The pipeline has a total length of 1840.85 km and 10 stations along its route. Among them, station N4 is a distribution station without a booster function. The first station N1 uses a two-main-transfer pump series operation mode, while stations N2-N3 and N5-N9 are booster pump stations, each equipped with a separate main-transfer pump. The pipeline sequentially transports methanol and gasoline. The two types of oil have significant differences in density, viscosity, and elastic modulus. Gasoline has a lower density than methanol, a dynamic viscosity of only 51.8% of methanol, and an elastic modulus that is about 72.8% higher than methanol. These differences are key factors affecting the hydraulic characteristics of the pipeline.

[0064] The core settings for the simulation test are as follows: The control strategy adopts the first-station flow setpoint control (1100 m³ / s). 3 The N4 station uses a combination of pressure setpoint control (7 MPa) and pressure setpoint control at the terminal station ( / h) at 100 m 3 / h flow rate distribution. Initially, the oil batches in the pipeline were distributed as batch 1 (gasoline), batch 2 (methanol), batch 3 (gasoline), and batch 4 (methanol), corresponding to initial oil head locations of 1840.85 km, 1105.43 km, 610.74 km, and 115.24 km from the first station, respectively. The total simulation duration was set to 50 hours. To verify the reliability and applicability of the simulation method, this embodiment designed four differentiated operating conditions, comprehensively covering sequential transportation quasi-steady-state conditions, intermediate distribution station start-up conditions, intermediate distribution station shutdown conditions, and sudden pump stoppage at intermediate pumping stations, among other practical engineering scenarios. The simulation results of the proposed method were compared and analyzed with field operation data.

[0065] See Figure 6 , Figure 7 and Figure 8 The results show that the system performs well under various operating conditions, with the pressure and flow rate errors detailed. Specifically, the average absolute error of pressure is controlled within 0.098 MPa, and the average absolute error of flow rate is within 0.86 m³. 3 Within / h. Taking the intermediate station pump shutdown bypass operation as an example: At time 1 hour, station N3 stops pumping and opens the bypass valve, simulating the pump shutdown bypass operation. The hydraulic parameter solution results of upstream station N2 and downstream station N7 of station N3 are selected to analyze the transient process. The pump unit of station N3 shuts down instantaneously at the 60th minute, and the flow rate decreases instantaneously. The outlet pressure at station N2 instantaneously increased from 12.49 MPa to 14.43 MPa, eventually stabilizing at 15.78 MPa. This was due to the sudden shutdown of pump N3, which reduced flow and transmitted a pressure surge upstream, causing a rapid increase in pressure along the upstream pipeline. Subsequently, as the bypass valve opened, the oil accumulation in the upstream pipeline gradually eased, and the pressure gradually stabilized. In this new steady state, the overall flow decreased, reducing friction along the pipeline and leading to a decrease in pressure at stations upstream of station N3. Similarly, the inlet pressure at station N7 instantaneously decreased from 4.3 MPa to 3.76 MPa, eventually stabilizing at 4.29 MPa. This was due to the sudden shutdown of pump N3 transmitting a pressure reduction surge downstream. As the bypass valve opened, oil flowed into the downstream pipeline through the bypass valve, and the pressure drop along the pipeline decreased as the flow decreased, causing the downstream pipeline pressure to rise and eventually stabilize. The overall flow rate showed an instantaneous decrease followed by a rebound, eventually stabilizing. The steady-state flow rate at station N2 upstream of the distribution station reached 1100 m³ / h. 3 / h, the flow rate at station N7 downstream of the distribution station stabilized at 1000 m³ / h. 3 / h. The trend of change shows good consistency with the field operation data, with the average absolute error of pressure being 0.095 MPa and the average absolute error of flow rate being 0.35 m³ / h. 3 / h.

[0066] The sequential transport of various liquid fuels in pipelines causes dynamic changes in the friction coefficient, pump station characteristic curves, and water hammer propagation velocity, leading to significant nonlinear characteristics in the pressure distribution, flow stability, and energy supply and demand relationship along the entire pipeline. This embodiment provides accurate transient simulation analysis of sequential liquid fuel transport, helping operation and maintenance personnel to fully understand the operational status of the entire pipeline and providing key technical support for formulating scientific control strategies and ensuring the long-term safe and stable operation of the pipeline.

[0067] like Figure 9 As shown, this embodiment proposes a hydraulic simulation system for sequential pipeline transport of multiple liquids based on dynamic interpolation, including: Select unit 101 is used to select two adjacent pipe segments for transporting different types of liquid fuels from multiple divided pipe segments in the simulation model of sequential pipeline transportation of multi-variety liquid fuels. The partitioning unit 102 is used to partition each adjacent pipe segment into a grid according to the maximum fuel wave velocity in the model, and to determine at least one grid node in each adjacent pipe segment. Monitoring unit 103 is used to monitor the relative positional relationship between the oil mixing interface and the grid nodes in two adjacent pipe sections; The determination unit 104 is used to determine the flow rate and pressure head of at least one grid node based on the relative position relationship, the fuel wave velocity of the grid node, the maximum fuel wave velocity, and the dynamic interpolation method.

[0068] It should be noted that the processing procedures for selecting unit 101, dividing unit 102, monitoring unit 103, and determining unit 104, and the beneficial effects thereof, can be referred to respectively. Figure 1 Steps S101 to S104 in the process will not be described again.

[0069] Optionally, the determining unit 104 is also used for: When the mixing interface is located at the position corresponding to the grid node, the flow rate and pressure head of at least one grid node are determined based on the fuel wave velocity, maximum fuel wave velocity and time interpolation of the grid node. When the mixing interface is located between two adjacent grid nodes, the flow rate and pressure head of at least one grid node are determined based on the fuel wave velocity, maximum fuel wave velocity, and spatial interpolation of the grid node.

[0070] Optionally, the determining unit 104 is also used for: When the mixing interface is located at the position corresponding to the grid node, for any grid node: If the fuel wave velocity of a grid node at the current time step is equal to the maximum fuel wave velocity, then the initial flow rate and initial pressure head of the grid node in the next time step are determined based on the flow rate and pressure head of the upstream and downstream grid nodes at the current time step, and are directly used as the final flow rate and final pressure head of the grid node in the next time step. If the fuel wave velocity of a grid node at the current time step is less than the maximum fuel wave velocity, then the initial flow rate and initial pressure head of the grid node in the next time step are determined based on the flow rate and pressure head of the upstream and downstream grid nodes at the current time step. The flow rate and pressure head of the grid node at the current time step, and the initial flow rate and initial pressure head of the grid node in the next time step are then input into the first model created according to the time interpolation method for interpolation to obtain the final flow rate and final pressure head of the grid node in the next time step.

[0071] Optionally, the first model is: ; ; ; or ; ; ; in, and upstream grid node At the current time step of Characteristic line coefficients; For downstream grid nodes At the current time step of Characteristic line coefficients; and upstream grid node At the current time step Flow rate and pressure head; and Upstream grid nodes At the current time step The inertial water hammer coefficient and friction coefficient, and Downstream grid nodes At the current time step The inertial water hammer coefficient and friction coefficient; To calculate the distance step; and For grid nodes In the next step The final flow rate and final head; For grid nodes At the current time step, the wave speed Maximum wave speed; and For grid nodes At the current time step Flow rate and pressure head; and For grid nodes In the next step The initial flow rate and initial pressure head.

[0072] Optionally, the determining unit 104 is also used for: When the mixing interface is located between a grid node and an upstream grid node, determine the fuel-average wave velocity of the mixing grid, the flow rate and head of the grid node at the current time step, and the initial flow rate and initial head of the upstream and downstream grid nodes at the current time step. The mixing grid is defined as the grid node. For the grid between upstream and downstream grid nodes, the average fuel wave velocity, maximum fuel wave velocity, flow rate and pressure head of the grid node at the current time step, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are input into the second model created by spatial interpolation for interpolation to obtain the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step. Based on the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step, the flow rate and pressure head of the grid node at the next time step are determined.

[0073] Alternatively, the second model is: ; ; ; ; ) ; ) ; ; in, For the mixed oil grid at the current time step The average wave velocity of fuel, For the mixed oil interface at the current time step Time and grid nodes The distance; The step size is obtained by multiplying the maximum fuel wave velocity by the duration. and For upstream grid nodes and grid nodes at the current time step fuel wave speed; and For the mixed oil grid at the current time step The average friction coefficient and the average inertial water hammer coefficient, and For the upstream grid node at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the grid node at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the upstream grid node at the current time step The final flow rate and the final pressure head, and For downstream grid nodes at the current time step The final flow rate and the final pressure head, and These represent the current time steps of the grid node and the upstream grid node, respectively. Initial traffic, and These represent the current time steps of the grid node and the upstream grid node, respectively. The initial pressure head, This represents the pressure head difference caused by the density difference at the oil mixing interface.

[0074] Optionally, the determining unit 104 is also used for: When the oil mixing interface is located at a grid node and downstream grid nodes When determining the position between the nodes, the fuel-average wave velocity of the mixed oil grid is determined, the flow rate and pressure head of the grid nodes at the current time step are determined, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are determined. The mixed oil grid is the grid node. With downstream grid nodes The average fuel wave velocity, maximum fuel wave velocity, flow rate and pressure head of the grid nodes at the current time step, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are input into the third model created by spatial interpolation for interpolation to obtain the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step. Based on the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step, the flow rate and pressure head of the grid nodes at the next time step are determined.

[0075] Alternatively, the third model is: ; ; ; ; ) ; ) ; ; ; ; in, For the mixed oil grid at the current time step The average wave velocity of fuel, For the mixed oil interface at the current time step Time and downstream grid nodes The distance; The step size is obtained by multiplying the maximum fuel wave velocity by the duration. and For downstream grid nodes and grid nodes at the current time step fuel wave speed; and For the mixed oil grid at the current time step The average friction coefficient and the average inertial water hammer coefficient, and For grid nodes At the current time step The coefficient of friction and the coefficient of inertial water hammer, and For downstream grid nodes at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the upstream grid node at the current time step The final flow rate and the final pressure head, and For downstream grid nodes at the current time step The final flow rate and the final pressure head, and These represent the current time steps of the grid node and the upstream grid node, respectively. Initial traffic, and These represent the current time steps of the grid node and the upstream grid node, respectively. The initial pressure head, This represents the pressure head difference caused by the density difference at the oil mixing interface.

[0076] Optionally, the above system also includes: The update unit is used to monitor the operational changes of the first, intermediate, and last stations after determining the flow rate and pressure head of at least one grid node based on the relative position relationship, the fuel wave velocity of the grid node, the maximum fuel wave velocity, and the dynamic interpolation method, and to update the flow rate and pressure at the boundary points in the model based on the monitored operational changes of the first, intermediate, and last stations.

[0077] This embodiment proposes a hydraulic simulation system for sequential transport of multi-variety liquid fuels via pipelines based on dynamic interpolation. From multiple segmented pipe sections in a simulation model of sequential transport of multi-variety liquid fuels, two adjacent pipe sections are selected for transporting different types of liquid fuels. Based on the maximum fuel wave velocity in the model, each adjacent pipe section is meshed, and at least one grid node is determined within each adjacent pipe section. The relative positional relationship between the mixing interface in the two adjacent pipe sections and the grid nodes is monitored. Based on the relative positional relationship, the fuel wave velocity of the grid nodes, the maximum fuel wave velocity, and the dynamic interpolation method, the flow rate and head of at least one grid node are determined. This embodiment can perform hydraulic simulation of sequential transport of multi-variety liquid fuels via pipelines based on dynamic interpolation, accurately depicting the sequential transport process of multi-variety liquid fuels, achieving precise simulation of sequential transport scenarios, and effectively enriching the hydraulic simulation methods for sequential transport of multi-variety liquid fuels via pipelines.

[0078] In this embodiment, the hydraulic simulation system for sequential transport of multi-variety liquids in pipelines based on dynamic interpolation is presented in the form of functional units. Here, a unit refers to an ASIC (Application Specific Integrated Circuit), a processor and memory that execute one or more software or fixed programs, and / or other devices that can provide the above functions.

[0079] This invention also provides a computer device having the above-described features. Figure 9 The system shown is a hydraulic simulation system for sequential pipeline transport of multiple liquids based on dynamic interpolation.

[0080] 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.

[0081] 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.

[0082] 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.

[0083] 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.

[0084] 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.

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

[0086] 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.

[0087] 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 hydraulic simulation method for sequential pipeline transport of multiple liquids based on dynamic interpolation, characterized in that, include: In the simulation model of sequential pipeline transportation of multiple types of liquid fuels, two adjacent pipeline segments are selected from the multiple segments to transport different types of liquid fuels. Based on the maximum fuel wave velocity in the model, each adjacent pipe segment is divided into grids, and at least one grid node in each adjacent pipe segment is determined. Monitor the relative positional relationship between the oil mixing interface in the two adjacent pipe sections and the grid nodes; Based on the relative positional relationship, the fuel wave velocity of the grid node, the maximum fuel wave velocity, and dynamic interpolation, determine the flow rate and pressure head of at least one of the grid nodes.

2. The method according to claim 1, characterized in that, The step of determining the flow rate and pressure head of at least one of the grid nodes based on the relative positional relationship, the fuel wave velocity of the grid node, the maximum fuel wave velocity, and dynamic interpolation includes: When the oil mixing interface is located at the position corresponding to the grid node, the flow rate and pressure head of at least one of the grid nodes are determined based on the fuel wave velocity of the grid node, the maximum fuel wave velocity, and the time interpolation method. When the oil mixing interface is located between two adjacent grid nodes, the flow rate and pressure head of at least one grid node are determined based on the fuel wave velocity of the grid node, the maximum fuel wave velocity, and spatial interpolation.

3. The method according to claim 2, characterized in that, When the oil mixing interface is located at the position corresponding to the grid node, determining the flow rate and head of at least one of the grid nodes based on the fuel wave velocity of the grid node, the maximum fuel wave velocity, and the time interpolation method includes: When the oil mixing interface is located at the position corresponding to the grid node, for any one of the grid nodes: If the fuel wave velocity of the grid node at the current time step is equal to the maximum fuel wave velocity, then the initial flow rate and initial pressure head of the grid node in the next time step are determined based on the flow rate and pressure head of the upstream and downstream grid nodes at the current time step, and are directly used as the final flow rate and final pressure head of the grid node in the next time step. If the fuel wave velocity of the grid node at the current time step is less than the maximum fuel wave velocity, then based on the flow rate and pressure head of the upstream and downstream grid nodes at the current time step, the initial flow rate and initial pressure head of the grid node at the next time step are determined. The flow rate and pressure head of the grid node at the current time step, and the initial flow rate and initial pressure head of the grid node at the next time step are then input into the first model created according to the time interpolation method for interpolation to obtain the final flow rate and final pressure head of the grid node at the next time step.

4. The method according to claim 3, characterized in that, The first model is: ； ； ； or ; ； ； in, and For the upstream grid node At the current time step of Characteristic line coefficients; For the downstream grid node At the current time step of Characteristic line coefficients; and For the upstream grid node At the current time step Flow rate and pressure head; and The upstream grid nodes are respectively At the current time step The inertial water hammer coefficient and friction coefficient, and The downstream grid nodes are respectively At the current time step The inertial water hammer coefficient and friction coefficient; To calculate the distance step; and For the grid node In the next step The final flow rate and final head; For the grid node At the current time step, the wave speed The maximum wave speed; and For the grid node At the current time step Flow rate and pressure head; and For the grid node In the next step The initial flow rate and initial pressure head.

5. The method according to claim 2, characterized in that, When the oil mixing interface is located between two adjacent grid nodes, the flow rate and head of at least one grid node are determined based on the fuel wave velocity of the grid node, the maximum fuel wave velocity, and spatial interpolation, including: When the mixing interface is located between the grid node and the upstream grid node, the fuel average wave velocity of the mixing grid is determined, the flow rate and pressure head of the grid node at the current time step are determined, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are determined. The mixing grid is the grid node. The average fuel wave velocity, the maximum fuel wave velocity, the flow rate and pressure head of the grid node at the current time step, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are input into the second model created according to the spatial interpolation method for interpolation to obtain the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step. Based on the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step, the flow rate and pressure head of the grid node at the next time step are determined.

6. The method according to claim 5, characterized in that, The second model is: ； ； ； ； ) ； ) ； ； in, For the mixed oil grid at the current time step The average wave velocity of fuel, For the mixed oil interface at the current time step Time and the grid node The distance; The step size is obtained by multiplying the maximum fuel wave velocity by the duration. and For the upstream grid node and the grid node at the current time step fuel wave speed; and For the mixed oil grid at the current time step The average friction coefficient and the average inertial water hammer coefficient, and For the upstream grid node at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the grid node at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the upstream grid node at the current time step The final flow rate and the final pressure head, and For the downstream grid node at the current time step The final flow rate and the final pressure head, and The grid node and the upstream grid node are respectively at the current time step. Initial traffic, and The grid node and the upstream grid node are respectively at the current time step. The initial pressure head, This represents the pressure head difference caused by the density difference at the interface of the oil mixture.

7. The method according to claim 2, characterized in that, When the oil mixing interface is located between two adjacent grid nodes, the flow rate and head of at least one grid node are determined based on the fuel wave velocity of the grid node, the maximum fuel wave velocity, and spatial interpolation, including: When the oil mixing interface is located at the grid node and downstream grid nodes When determining the position between the nodes, the average fuel wave velocity of the mixed oil grid is determined, the flow rate and head of the grid node at the current time step are determined, and the initial flow rate and initial head of the upstream and downstream grid nodes at the current time step are determined, wherein the mixed oil grid is the grid node. With the downstream mesh node The average fuel wave velocity, the maximum fuel wave velocity, the flow rate and pressure head of the grid node at the current time step, and the initial flow rate and initial pressure head of the upstream and downstream grid nodes at the current time step are input into a third model created according to spatial interpolation for interpolation to obtain the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step. Based on the final flow rate and final pressure head of the upstream and downstream grid nodes at the current time step, the flow rate and pressure head of the grid node at the next time step are determined.

8. The method according to claim 7, characterized in that, The third model is: ； ； ； ； ) ； ) ； ； ； ； in, For the mixed oil grid at the current time step The average wave velocity of fuel, For the mixed oil interface at the current time step Time and the downstream grid node The distance; The step size is obtained by multiplying the maximum fuel wave velocity by the duration. and For the downstream mesh node and the mesh node at the current time step fuel wave speed; and For the mixed oil grid at the current time step The average friction coefficient and the average inertial water hammer coefficient, and For the grid node At the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the downstream grid node at the current time step The coefficient of friction and the coefficient of inertial water hammer, and For the upstream grid node at the current time step The final flow rate and the final pressure head, and For the downstream grid node at the current time step The final flow rate and the final pressure head, and The grid node and the upstream grid node are respectively at the current time step. Initial traffic, and The grid node and the upstream grid node are respectively at the current time step. The initial pressure head, This represents the pressure head difference caused by the density difference at the interface of the oil mixture.

9. The method according to any one of claims 1 to 8, characterized in that, After determining the flow rate and head of at least one of the grid nodes based on the relative positional relationship, the fuel wave velocity of the grid node, the maximum fuel wave velocity, and dynamic interpolation, the method further includes: Monitor operational changes at the first, intermediate, and final stations, and update the flow and pressure at the boundary points in the model based on the monitored operational changes at the first, intermediate, and final stations.

10. A hydraulic simulation system for sequential pipeline transport of multiple liquids based on dynamic interpolation, characterized in that, include: Selecting a unit is used to select two adjacent pipe segments for transporting different types of liquid fuels from multiple divided pipe segments in a simulation model of sequential pipeline transportation of multiple types of liquid fuels. A partitioning unit is used to partition each of the adjacent pipe segments into grids according to the maximum fuel wave velocity in the model, and to determine at least one grid node in each of the adjacent pipe segments. The monitoring unit is used to monitor the relative positional relationship between the oil mixing interface in the two adjacent pipe sections and the grid node; The determining unit is used to determine the flow rate and pressure head of at least one of the grid nodes based on the relative positional relationship, the fuel wave velocity of the grid node, the maximum fuel wave velocity, and dynamic interpolation.

Images