Steam simulation method and device based on flow direction pre-solution, equipment and medium

Abstract

The application discloses a steam simulation method and device based on flow direction pre-solution, equipment and medium, relates to the field of steam simulation, and comprises the following steps: setting steam flow for the end node of the current steam pipe network, adjusting the steam source of the pipe network according to the steam flow increase value of the pipe network, performing hydraulic calculation pre-solution on a preset first mathematical model according to the obtained adjusted virtual steam pipe network, correcting the pipe flow direction of the pipe network through the obtained hydraulic pre-solution result, and obtaining a to-be-mapped virtual steam pipe network; mapping the hydraulic pre-solution result to a preset fine grid, determining the leakage of the mapped pipe network, and correcting the leakage based on the leakage; solving a preset second mathematical model through the obtained corrected virtual steam pipe network, obtaining the operation parameters of each node in the pipe network; correcting the condensation of the pipe network by using the condensation distribution of each node determined based on the operation parameters, and performing steam simulation through the obtained target virtual steam pipe network. Therefore, the accuracy of steam simulation can be improved.

Description

Technical Field

[0001] This invention relates to the field of steam simulation, and in particular to a steam simulation method, apparatus, equipment, and medium based on pre-solution of flow direction. Background Technology

[0002] As industrial systems increasingly demand higher levels of energy efficiency management and intelligent operation, refined modeling and simulation of steam pipeline networks have become a key support for ensuring the efficient operation of thermal systems. To achieve system design optimization, energy consumption diagnosis, and operational scheme evaluation, engineering widely relies on steady-state simulation tools to predict pressure, flow rate, temperature, and energy distribution within the pipeline network.

[0003] In the simulation modeling of pipeline networks, it is necessary to solve both the hydraulic balance equations and the heat transfer equations simultaneously. At some nodes, a nodal enthalpy mixing model is typically used, where steam with different enthalpies upstream is uniformly mixed and then transferred downstream. Therefore, the thermal calculation relies on the fluid flow direction to determine the enthalpy transfer path. However, in some branches of the pipeline network or under low-load conditions, the pipeline flow rate may approach zero. During the numerical solution process, due to rounding errors or iteration fluctuations, small but uncertain flow rates can easily occur, leading to inaccurate flow direction judgments. Since the enthalpy transfer direction is strongly correlated with the flow direction, such misjudgments can cause rank deficiencies in the thermal calculation equations, leading to complete computational failure.

[0004] On the other hand, as steam travels through the pipeline, it continuously dissipates heat into the environment, producing condensate. This condensate is periodically discharged from the system through steam traps, constituting a mass and energy loss during actual operation. Existing simulation modeling methods typically estimate condensation based on ideal insulation conditions and standard heat dissipation models, neglecting the influence of actual on-site factors. This results in actual condensate volumes generally exceeding theoretical calculations and being difficult to accurately characterize using fixed parameters. Traditional methods lack an effective compensation mechanism for this deviation, causing the overall network thermal calculation results to deviate from actual operating conditions. This is especially true in long-distance transport or multi-stage pressure reduction systems, where the error amplifies at each stage, weakening the engineering guidance value of the simulation. Summary of the Invention

[0005] In view of this, the purpose of this invention is to provide a steam simulation method, apparatus, equipment, and medium based on flow direction pre-solution, which can effectively improve the accuracy of hydrophobicity modeling of steam pipeline networks through synergistic optimization of flow direction determination and condensation correction. The specific solution is as follows:

[0006] In a first aspect, this application discloses a steam simulation method based on flow direction pre-solution, including:

[0007] Set the steam output of the terminal node in the current steam network to obtain the corresponding initial virtual steam network, and calculate the increase in steam output of the initial virtual steam network to allocate the increase in steam output to the steam source of the initial virtual steam network to obtain the adjusted virtual steam network.

[0008] The adjusted virtual steam network is used to perform hydraulic calculations on the preset first mathematical model, and the pipe flow direction of the adjusted virtual steam network is corrected based on the obtained hydraulic calculation results to obtain the virtual steam network to be mapped.

[0009] The hydraulic pre-solution results of the virtual steam network to be mapped are mapped to a preset fine mesh to obtain the virtual steam network to be corrected, and the leakage of the virtual steam network to be corrected is determined so as to correct the virtual steam network to be corrected based on the leakage to obtain the corrected virtual steam network.

[0010] The preset second mathematical model is solved based on the corrected virtual steam network to obtain the operating parameters of each node in the corrected virtual steam network. The condensation distribution of each node in the corrected virtual steam network is determined based on the operating parameters to correct the condensation of the corrected virtual steam network, so as to obtain the target virtual steam network. Steam simulation is then performed based on the target virtual steam network.

[0011] Optionally, the step of setting the steam output of the end nodes in the current steam network to obtain a corresponding initial virtual steam network, and calculating the increase in the steam output of the initial virtual steam network to allocate the increase in the steam output to the steam source of the initial virtual steam network, in order to obtain an adjusted virtual steam network, includes:

[0012] Identify the terminal nodes of the current steam network and set the steam output of the terminal nodes to a preset steam output to obtain an initial virtual steam network;

[0013] The increase in steam output of the initial virtual steam network is statistically analyzed, and the increase in steam output is allocated to the steam source nodes of the initial virtual steam network based on a preset allocation ratio to obtain the adjusted virtual steam network.

[0014] Optionally, the step of performing hydraulic calculations to pre-solve the preset first mathematical model using the adjusted virtual steam network, and correcting the pipe flow direction of the adjusted virtual steam network using the obtained hydraulic pre-solution results to obtain the virtual steam network to be mapped, includes:

[0015] The hydraulic calculations of the preset first mathematical model are pre-solved using the adjusted virtual steam network to obtain the hydraulic pre-solution results corresponding to the adjusted virtual steam network; the hydraulic pre-solution results include the mass distribution, pressure distribution, and flow distribution corresponding to the virtual steam network;

[0016] Based on the mass distribution, pressure distribution, and flow distribution, the current pipeline flow direction in the adjusted virtual steam network is determined, and the upstream and downstream relationships of the adjusted virtual steam network are adjusted based on the pipeline flow direction to obtain the virtual steam network to be mapped.

[0017] Optionally, mapping the hydraulic pre-solution results of the virtual steam network to be mapped to a preset fine mesh to obtain the virtual steam network to be corrected includes:

[0018] The pre-solved results of the virtual steam pipeline network to be mapped are mapped to a preset fine grid using a preset numerical interpolation algorithm to obtain the virtual steam pipeline network to be corrected; the preset fine grid is a grid with a higher resolution than the current grid of the virtual steam pipeline network to be corrected.

[0019] Optionally, determining the leakage amount of the virtual steam network to be corrected, and correcting the virtual steam network to be corrected based on the leakage amount to obtain a corrected virtual steam network, includes:

[0020] Calculate the current total input flow and total output flow of the virtual steam network to be corrected;

[0021] If the total input flow rate is not less than the total output flow rate, then the virtual steam network to be corrected is taken as the corrected virtual steam network;

[0022] If the total input flow is less than the total output flow, the leakage is determined based on the total input flow and the total output flow, and the total output flow is compressed to be equal to the total input flow based on the leakage to obtain the corrected virtual steam network.

[0023] Optionally, the step of solving the preset second mathematical model based on the corrected virtual steam network to obtain the operating parameters of each node in the corrected virtual steam network includes:

[0024] The preset second mathematical model is spatially discretized based on the preset numerical interpolation algorithm to obtain the discretized mathematical model. The discretized mathematical model is then iteratively solved using the corrected virtual steam network to obtain the operating parameters of each node in the corrected virtual steam network. The preset second mathematical model is a model that introduces an equivalent leakage coefficient into the preset first mathematical model. The equivalent leakage coefficient is a coefficient that characterizes the intensity of condensation mass loss.

[0025] Optionally, the step of using the condensation distribution of each node in the corrected virtual steam network determined based on the operating parameters to correct the condensation of the corrected virtual steam network, so as to obtain the target virtual steam network, and performing steam simulation based on the target virtual steam network includes:

[0026] The condensation distribution of each node in the corrected virtual steam network is determined based on the operating parameters.

[0027] If the condensation distribution is negative, the condensation of each node is reduced based on the steam consumption ratio of each node until the total input flow of the corrected virtual steam network is equal to the total output flow, so as to obtain the target virtual steam network, and steam simulation is performed based on the target virtual steam network.

[0028] Secondly, this application discloses a steam simulation device based on flow direction pre-solution, comprising:

[0029] The steam pipeline network adjustment module is used to set the steam output of the end nodes in the current steam pipeline network to obtain the corresponding initial virtual steam pipeline network, and to count the increase in the steam output of the initial virtual steam pipeline network, so as to allocate the increase in the steam output to the steam source of the initial virtual steam pipeline network to obtain the adjusted virtual steam pipeline network.

[0030] The pipeline flow direction correction module is used to perform hydraulic calculation pre-solution on the preset first mathematical model through the adjusted virtual steam pipeline network, and correct the pipeline flow direction of the adjusted virtual steam pipeline network through the obtained hydraulic pre-solution results to obtain the virtual steam pipeline network to be mapped.

[0031] The steam pipeline network correction module is used to map the hydraulic pre-solution results of the virtual steam pipeline network to be mapped to a preset fine mesh to obtain the virtual steam pipeline network to be corrected, and to determine the leakage of the virtual steam pipeline network to be corrected, so as to correct the virtual steam pipeline network to be corrected based on the leakage, and obtain the corrected virtual steam pipeline network.

[0032] The condensation correction module is used to solve a preset second mathematical model based on the corrected virtual steam network to obtain the operating parameters of each node in the corrected virtual steam network, and to correct the condensation of the corrected virtual steam network using the condensation distribution of each node in the corrected virtual steam network determined based on the operating parameters, so as to obtain a target virtual steam network, and to perform steam simulation based on the target virtual steam network.

[0033] Thirdly, this application discloses an electronic device, including:

[0034] Memory, used to store computer programs;

[0035] A processor is used to execute the computer program to implement the steam simulation method based on flow direction pre-solution as described above.

[0036] Fourthly, this application discloses a computer-readable storage medium for storing a computer program, wherein the computer program, when executed by a processor, implements the steam simulation method based on flow direction pre-solution as described above.

[0037] In this application, an initial virtual steam network is obtained by setting the steam output of the terminal nodes in the current steam network, and the increase in steam output of the initial virtual steam network is statistically analyzed. This increase is then allocated to the steam sources of the initial virtual steam network to obtain an adjusted virtual steam network. The adjusted virtual steam network is used to perform hydraulic calculations on a preset first mathematical model, and the obtained hydraulic calculation results are used to correct the pipe flow direction of the adjusted virtual steam network to obtain the virtual steam network to be mapped. The hydraulic calculation results of the virtual steam network to be mapped are then mapped to a pre-defined virtual steam network. A fine mesh is set to obtain a virtual steam network to be corrected, and the leakage of the virtual steam network to be corrected is determined. The virtual steam network to be corrected is then corrected based on the leakage to obtain a corrected virtual steam network. A preset second mathematical model is solved based on the corrected virtual steam network to obtain the operating parameters of each node in the corrected virtual steam network. The condensation distribution of each node in the corrected virtual steam network, determined based on the operating parameters, is used to correct the condensation of the corrected virtual steam network to obtain a target virtual steam network. Steam simulation is then performed based on the target virtual steam network.

[0038] Therefore, the method of this application allows setting the steam output of the terminal nodes in the current steam network, which is equivalent to modifying the current steam network and obtaining an initial virtual steam network. Then, the increase in steam output of the initial virtual steam network needs to be statistically analyzed and allocated to the steam source of the initial virtual steam network, which is equivalent to modifying the network again to obtain an adjusted virtual steam network. Furthermore, the adjusted virtual steam network is used to perform hydraulic calculations on a preset first mathematical model, and the obtained pre-solution results are used to correct the flow direction of the adjusted virtual steam network to obtain the virtual steam network to be mapped. The network needs to be mapped from the coarse grid of the current pipeline network to a finer grid with higher resolution to obtain the virtual steam pipeline network to be corrected. Then, the leakage of the virtual steam pipeline network to be corrected is used to obtain the corrected virtual steam pipeline network. Finally, the second mathematical model is solved based on the corrected virtual steam pipeline network to obtain the operating parameters of each node in the corrected virtual steam pipeline network. The condensation distribution of each node in the corrected virtual steam pipeline network is determined based on the operating parameters to correct the condensation of the corrected virtual steam pipeline network to obtain the target virtual steam pipeline network. Then, the obtained target virtual steam pipeline network is used for steam simulation. In this way, on the one hand, a preprocessing mechanism can be provided to accurately obtain the flow direction of each pipe in the pipeline network before formal solution, so as to eliminate the uncertainty of flow direction caused by numerical error, thereby improving the stability and robustness of subsequent calculations; on the other hand, by correcting the condensation error, the implicitly increased condensate in actual operation can be reasonably distributed to each section of the pipeline network, making up for the gap between the ideal model and the actual working conditions, and improving the accuracy of thermal simulation results; furthermore, through the synergistic optimization of flow direction judgment and condensation correction, it can be ensured that the simulation model can still maintain good convergence performance and engineering credibility under complex multi-source, ring topology and dynamic working conditions. Attached Figure Description

[0039] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0040] Figure 1 This is a flowchart of a steam simulation method based on flow direction pre-solution disclosed in this application;

[0041] Figure 2 This is a schematic diagram of a steam simulation device based on flow direction pre-solution disclosed in this application;

[0042] Figure 3This is a structural diagram of an electronic device disclosed in this application. Detailed Implementation

[0043] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0044] Currently, in the simulation modeling of pipeline networks, the flow rate in some branches or under low-load conditions may approach zero. During numerical solutions, rounding errors or iteration fluctuations can easily lead to small but uncertain flow rates, resulting in inaccurate flow direction judgments and consequently, calculation errors. Furthermore, steam continuously dissipates heat to the environment during its transmission within the pipeline, producing condensate. Existing simulation modeling methods typically estimate condensation based on ideal insulation conditions and standard heat dissipation models, neglecting the influence of actual on-site factors. This results in actual condensate volumes generally exceeding theoretical calculations and being difficult to accurately characterize using fixed parameters.

[0045] To overcome the above problems, this application discloses a steam simulation method, apparatus, equipment and medium based on flow direction pre-solution, which can effectively improve the accuracy of hydrophobic modeling of steam pipelines through the synergistic optimization of flow direction judgment and condensation correction.

[0046] See Figure 1 As shown, this embodiment of the invention discloses a steam simulation method based on flow direction pre-solution, including:

[0047] Step S11: Set the steam output of the end node in the current steam network to obtain the corresponding initial virtual steam network, and calculate the increase in steam output of the initial virtual steam network to allocate the increase in steam output to the steam source of the initial virtual steam network to obtain the adjusted virtual steam network.

[0048] In this embodiment, the steam output of the terminal nodes in the steam pipeline network needs to be set first. Then, the increase in the steam output of the pipeline network is calculated, and the steam pipeline network is adjusted based on the increase in steam output. Specifically, the terminal nodes of the current steam pipeline network need to be identified first, and the steam output of the terminal nodes is set to a preset steam output to obtain the initial virtual steam pipeline network. It should be noted that the terminal nodes include steam sources without flow, user nodes, pipeline plugs, closed valves in the pipeline network, and other equipment. These nodes are set to user type during the pre-solution process, and users with existing flow rates greater than 0.01 kg / s still use the original collected values ​​for calculation. Other user points are set to a preset steam output of 0.01 kg / s, which is much larger than the simulation flow tolerance. This is significantly higher than the computer's truncation error, but lower than typical user gas consumption. Further, it is necessary to calculate the increase in steam output of the initial virtual steam network and allocate this increase to the steam source nodes of the initial virtual steam network based on a preset allocation ratio to obtain the adjusted virtual steam network. The preset allocation ratio can be set according to actual needs. The purpose of allocating the increase in steam output to the steam source nodes of the initial virtual steam network is to avoid hydraulic solution failure due to changes in the total inlet and outlet pattern.

[0049] Step S12: Perform hydraulic calculations on the preset first mathematical model using the adjusted virtual steam network, and correct the pipe flow direction of the adjusted virtual steam network using the obtained hydraulic pre-solution results to obtain the virtual steam network to be mapped.

[0050] In this embodiment, a hydraulic calculation pre-solution is performed on the preset first mathematical model using the adjusted virtual steam pipe network. Specifically, the adjusted virtual steam pipe network is used to perform a hydraulic calculation pre-solution on the preset first mathematical model to obtain the hydraulic pre-solution results corresponding to the adjusted virtual steam pipe network. The hydraulic pre-solution results include the mass distribution, pressure distribution, and flow distribution corresponding to the virtual steam pipe network. It should be noted that the finite volume method (FVM), finite difference method (FDM), or other numerical methods suitable for discretizing partial differential equations can be used to discretize the mass and momentum conservation equations within the pipe network, and the mass distribution, pressure distribution, and flow distribution of each pipe segment can be solved using an iterative algorithm.

[0051] It should be noted that the first mathematical model is assumed to be as follows:

[0052] ;

[0053] ;

[0054] ;

[0055] Where m is the mass flow rate, in kg / s. The distance along the pipe is expressed in meters (m). The distance along the pipe is expressed in meters (m). The cross-sectional area of ​​the pipe is expressed in units of... , The condensation mass is expressed in kg / s. This is the steam pressure, in Pa. For the pipeline friction coefficient, Latent heat of condensation, unit: J / kg The heat dissipation of the pipe is expressed in W / s, and g is the acceleration due to gravity, expressed in m / s². 2 The value is taken as 9.8, θ is the pipe inclination angle, and h is the gas enthalpy, in J / (kg). K) and z represent the elevation difference of the corresponding micro-element pipe, in meters (m); v represents the flow velocity, in meters per second (m / s); and d represents the pipe diameter, in meters (m).

[0056] Furthermore, it should be noted that if an equivalent leakage coefficient is not introduced during the main simulation stage, the standard hydraulic equations in the preset first mathematical model are used for pre-solution. If an equivalent leakage coefficient is introduced during the main simulation stage, the pre-solution needs to be performed according to the preset second mathematical model. However, it should be noted that since the pipeline flow direction is not yet determined in the initial stage of pre-solution, directly introducing nonlinear leakage terms related to mass flow rate may lead to ill-conditioned equations, iterative divergence, or convergence to non-physical solutions, thereby weakening the stability of flow direction determination. Therefore, to ensure the robustness of the pre-solution process, the standard hydraulic equations in the preset first mathematical model are used for pre-solution before the pipeline flow direction of the steam network is determined.

[0057] After obtaining the pre-solution results, it is necessary to determine the current pipe flow direction in the adjusted virtual steam network based on the mass distribution, pressure distribution, and flow rate distribution in the pre-solution results. Specifically, the flow direction of each pipe in the steam network needs to be determined based on the sign of the flow rate distribution, and the upstream and downstream relationships of the adjusted virtual steam network are adjusted based on the pipe flow direction to reconstruct the upstream and downstream relationships in the entire steam network, thereby obtaining the virtual steam network to be mapped. In this way, by setting virtual low-flow boundaries and performing hydraulic pre-solution before the formal solution, it is ensured that all pipe segments in the network have a clear flow direction, effectively avoiding the problem of misjudgment of flow direction caused by numerical oscillations under low-flow or zero-flow conditions.

[0058] Step S13: Map the hydraulic pre-solution results of the virtual steam network to be mapped to a preset fine mesh to obtain the virtual steam network to be corrected, and determine the leakage of the virtual steam network to be corrected, so as to correct the virtual steam network to be corrected based on the leakage, and obtain the corrected virtual steam network.

[0059] In this embodiment, the hydraulic pre-solution results need to be mapped to a preset fine mesh to improve modeling accuracy. Specifically, a preset numerical interpolation algorithm is used to map the pre-solution results of the virtual steam pipe network to be mapped to a preset fine mesh to obtain the virtual steam pipe network to be corrected. The preset fine mesh is a mesh with a higher resolution than the current mesh of the virtual steam pipe network to be corrected. It should be noted that the pre-solution results obtained after pre-solution constitute a set of physical field solutions that are highly close to the actual operating state. To improve the computational efficiency of the main simulation stage, the pre-solution results can be used as the initial field input for the hydraulic-thermal coupled simulation. Furthermore, since a relatively coarse pipe discrete mesh (coarse mesh) is used in the pre-solution stage, a finer mesh (fine mesh) can be used in the subsequent simulation stage to improve the spatial resolution of thermal and flow parameters. Furthermore, since there may be some network topology differences between coarse and fine grids, such as different numbers of pipe segments or inconsistent node positions, numerical interpolation methods, such as linear interpolation, piecewise constant mapping, or spline interpolation, are used to map the pressure and mass flow rate data on the coarse grid to the nodes of the fine grid, generating a continuous and physically reasonable initial field distribution. This not only enables the effective reuse of pre-solved results but also decouples the grid requirements for flow direction determination from high-precision simulation, further improving the flexibility and computational economy of the simulation system while ensuring the accuracy of flow direction determination. Moreover, flow direction pre-solving can be completed quickly on the coarse grid, while the main simulation can be performed at high resolution on the fine grid, with initial field transfer achieved through interpolation. This multi-level grid strategy effectively reduces overall computational overhead while ensuring computational accuracy.

[0060] Step S14: Solve the preset second mathematical model based on the corrected virtual steam network to obtain the operating parameters of each node in the corrected virtual steam network, and use the condensation distribution of each node in the corrected virtual steam network determined based on the operating parameters to correct the condensation of the corrected virtual steam network to obtain the target virtual steam network, and perform steam simulation based on the target virtual steam network.

[0061] In this embodiment, due to potential underestimation of condensation volume caused by non-explicit modeling factors in the steam pipeline network simulation model, such as insulation aging, additional heat dissipation from condensate branch pipes, and environmental fluctuations, it is necessary to correct the condensation volume of the steam pipeline network using an equivalent leakage calculation and allocation method based on system-level mass balance. Specifically, it is necessary to statistically analyze the current total input flow rate and total output flow rate of the virtual steam pipeline network to be corrected. It should be noted that in actual steam pipeline network operation, condensate is continuously generated and discharged from the system due to factors such as pipeline heat dissipation and equipment condensation, typically resulting in a total input mass flow rate greater than the total output mass flow rate, indicating a reasonable mass loss. Therefore, it is necessary to perform consistency verification on the collected boundary flow data to eliminate mass imbalance problems caused by instrument errors or data anomalies. Thus, it is first necessary to determine the total input flow rate and total output flow rate of the steam pipeline network.

[0062] Furthermore, this can be divided into two specific scenarios. In the first scenario, if the total input flow rate is not less than the total output flow rate, the virtual steam network to be corrected is used as the corrected virtual steam network. In the second scenario, if the total input flow rate is less than the total output flow rate, the leakage rate is determined based on both the total input and output flow rates, and the total output flow rate is compressed to equal the total input flow rate based on the leakage rate to obtain the corrected virtual steam network. It should be noted that when the total input is not less than the total output, this state conforms to physical reality, and the data can be considered reasonable; therefore, the original value is retained for subsequent equivalent leakage calculations. When the total input is less than the total output, a negative mass loss may occur, violating the law of mass conservation, indicating that the user-side measurement value may have an excessively high error or data distortion. However, considering that the steam source end, such as the boiler outlet and main steam header, is usually equipped with high-precision flow meters and undergoes more stringent management and maintenance, its measurement data has higher reliability. Therefore, using the total input flow rate at the steam source as the benchmark value, the total output on the user side is corrected to keep the relative load ratio between users unchanged, compress the total output to be equal to the total input, eliminate the mass surplus, and thus obtain the corrected virtual steam network, avoiding negative values ​​or non-physical calculations in subsequent condensation calculations.

[0063] Next, it is necessary to solve the preset second mathematical model based on the corrected virtual steam pipeline network to obtain the operating parameters of each node in the corrected virtual steam pipeline network. Specifically, it is necessary to spatially discretize the preset second mathematical model based on a preset numerical interpolation algorithm to obtain a discretized mathematical model, and then iteratively solve the discretized mathematical model through the corrected virtual steam pipeline network to obtain the operating parameters of each node in the corrected virtual steam pipeline network. It should be noted that the preset numerical interpolation algorithm in this embodiment adopts the FDM algorithm, which spatially discretizes the pipeline control equations using the FDM algorithm, adopts a backward difference format, and employs a hydraulic and thermal decoupling iterative strategy to improve computational stability and efficiency.

[0064] It should be further explained that the preset second mathematical model is a model that introduces an equivalent leakage coefficient into the preset first mathematical model, and the equivalent leakage coefficient α is a coefficient characterizing the intensity of condensation mass loss. Specifically, to achieve quantitative modeling of the equivalent leakage, an equivalent leakage coefficient α is introduced to characterize the intensity of condensation mass loss under a unit distribution benchmark. Considering that the additional heat dissipation caused by the condensate branch pipes and their connecting structures in the steam pipeline network is one of the main sources of equivalent leakage, and that in actual engineering, condensate valves are usually installed at certain intervals along the pipeline and are relatively evenly distributed, this embodiment preferentially uses the pipeline length as the distribution benchmark for the equivalent leakage. In specific scenarios, if it is known that the main heat loss of the system originates from the degradation of the insulation layer performance, such as aging or moisture, the outer surface area of ​​the pipe can be used as the distribution benchmark to more accurately reflect the distribution characteristics of radial heat dissipation. The formula of the preset second mathematical model is expressed as follows:

[0065] ;

[0066] ;

[0067] ;

[0068] Where α is the equivalent leakage coefficient, and other parameters can be found in the parameter explanation in step S12.

[0069] Finally, after completing the pre-solution of flow direction and the allocation of equivalent leakage, the main simulation calculation stage begins. This stage requires determining the condensation distribution of each node in the corrected virtual steam network based on the operating parameters, in order to correct the condensation of the virtual steam network based on this distribution. It should be noted that since the condensation is a calculated value within the model, it may cause a new imbalance between the total output mass of the entire network (i.e., the sum of user steam consumption and calculated condensation) and the input mass, thus conflicting with the initially set boundary conditions, manifesting as a mass conservation residual. This residual will be incorporated into the calculation of the equivalent leakage coefficient to eliminate it. In certain cases, if the total condensation calculated in the simulation is greater than the actual observed total inflow / outflow difference, the obtained equivalent total leakage will be negative. However, considering that the condensation in the steam network usually accounts for only a small proportion of the total transport flow (generally less than 1%), the absolute value of this negative equivalent leakage coefficient is extremely small. Its physical meaning is equivalent to uniformly injecting a small amount of steam mass into each pipe section to compensate for the small deviation between the model and the actual measurement. Since the injection volume is much smaller than the normal operating flow rate, its impact on the overall pressure, temperature, and velocity distribution is negligible and will not significantly change the engineering reliability of the simulation results. However, to achieve a more refined global mass balance, the steam network can be corrected using the following method: The condensation distribution of each node in the corrected virtual steam network is determined based on the operating parameters. If the condensation distribution is negative, the condensation of each node is reduced based on the steam consumption ratio of each node until the total input flow rate of the corrected virtual steam network equals the total output flow rate, thus obtaining the target virtual steam network. Steam simulation is then performed using this target virtual steam network. In other words, when the equivalent leakage coefficient is negative, the user-side load can be adjusted in reverse, i.e., the calculated condensation value is reduced proportionally according to the steam consumption ratio of each user until the total inflow and total outflow (including condensation) are completely matched. This method avoids the introduction of negative leakage terms, ensuring that all mass changes conform to the physical outflow logic. In this way, to address the discrepancy in condensation between the ideal model and actual operation caused by factors such as insulation aging and heat dissipation from condensate branch pipes, the equivalent leakage rate is used as the basis for correction based on the system-level mass difference, dynamically compensating for additional condensation losses not modeled. This method can achieve a reasonable restoration of the overall energy balance without relying on local high-precision measurements, thereby improving the accuracy of steam simulation.

[0070] In this embodiment, the steam output can be set for the terminal nodes of the current steam network, which is equivalent to modifying the current steam network to obtain an initial virtual steam network. Then, the increase in steam output of the initial virtual steam network needs to be calculated and allocated to the steam source of the initial virtual steam network, which is equivalent to modifying the network again to obtain an adjusted virtual steam network. Furthermore, the adjusted virtual steam network is used to perform hydraulic calculations on a preset first mathematical model, and the pipeline flow direction of the adjusted virtual steam network is corrected based on the obtained pre-solution results to obtain the virtual steam network to be mapped. The pre-solution results need to be mapped from the coarse grid of the current pipeline network to a finer grid with higher resolution to obtain the virtual steam pipeline network to be corrected. Then, the leakage of the virtual steam pipeline network to be corrected is used to obtain the corrected virtual steam pipeline network. Finally, the preset second mathematical model is solved based on the corrected virtual steam pipeline network to obtain the operating parameters of each node in the corrected virtual steam pipeline network. The condensation distribution of each node in the corrected virtual steam pipeline network, determined based on the operating parameters, is used to correct the condensation of the corrected virtual steam pipeline network to obtain the target virtual steam pipeline network. Finally, the obtained target virtual steam pipeline network is used for steam simulation. This approach offers several advantages. First, it provides a preprocessing mechanism to accurately determine the flow direction of each pipe in the network before formal solution calculation, eliminating flow direction uncertainty caused by numerical errors and improving the stability and robustness of subsequent calculations. Second, by correcting for condensation errors, the implicitly increased condensate volume during actual operation can be reasonably distributed to different sections of the network, bridging the gap between the ideal model and real-world conditions and improving the accuracy of thermal simulation results. Third, through the synergistic optimization of flow direction determination and condensation correction, the simulation model can maintain good convergence performance and engineering credibility even under complex multi-source, ring topology, and dynamic operating conditions. Finally, the method described in this application does not rely on manually preset flow directions or prior operating data; only conventional boundary flow information is required to complete the entire modeling process. For situations such as incomplete data or instrument errors, the user-side load is automatically corrected through a total inflow and outflow verification mechanism, ensuring quality conservation closure. Even in the extreme case where the equivalent leakage coefficient is negative, its impact is minimal and can be further optimized through optional strategies to ensure robust algorithm operation.

[0071] See Figure 2 As shown, this embodiment of the invention discloses a steam simulation device based on flow direction pre-solution, comprising:

[0072] The steam pipeline network adjustment module 11 is used to set the steam output of the end node of the current steam pipeline network to obtain the corresponding initial virtual steam pipeline network, and to count the increase in the steam output of the initial virtual steam pipeline network to allocate the increase in the steam output to the steam source of the initial virtual steam pipeline network to obtain the adjusted virtual steam pipeline network.

[0073] The pipeline flow direction correction module 12 is used to perform hydraulic calculation pre-solution on the preset first mathematical model through the adjusted virtual steam pipeline network, and correct the pipeline flow direction of the adjusted virtual steam pipeline network through the obtained hydraulic pre-solution results to obtain the virtual steam pipeline network to be mapped.

[0074] The steam pipeline network correction module 13 is used to map the hydraulic pre-solution result of the virtual steam pipeline network to be mapped to a preset fine mesh to obtain the virtual steam pipeline network to be corrected, and to determine the leakage of the virtual steam pipeline network to be corrected, so as to correct the virtual steam pipeline network to be corrected based on the leakage, and obtain the corrected virtual steam pipeline network.

[0075] The condensation correction module 14 is used to solve a preset second mathematical model based on the corrected virtual steam network to obtain the operating parameters of each node in the corrected virtual steam network, and to correct the condensation of the corrected virtual steam network using the condensation distribution of each node in the corrected virtual steam network determined based on the operating parameters, so as to obtain a target virtual steam network, and to perform steam simulation based on the target virtual steam network.

[0076] In this embodiment, the steam output can be set for the terminal nodes of the current steam network, which is equivalent to modifying the current steam network to obtain an initial virtual steam network. Then, the increase in steam output of the initial virtual steam network needs to be calculated and allocated to the steam source of the initial virtual steam network, which is equivalent to modifying the network again to obtain an adjusted virtual steam network. Furthermore, the adjusted virtual steam network is used to perform hydraulic calculations on a preset first mathematical model, and the pipeline flow direction of the adjusted virtual steam network is corrected based on the obtained pre-solution results to obtain the virtual steam network to be mapped. The pre-solution results need to be mapped from the coarse grid of the current pipeline network to a finer grid with higher resolution to obtain the virtual steam pipeline network to be corrected. Then, the leakage of the virtual steam pipeline network to be corrected is used to obtain the corrected virtual steam pipeline network. Finally, the preset second mathematical model is solved based on the corrected virtual steam pipeline network to obtain the operating parameters of each node in the corrected virtual steam pipeline network. The condensation distribution of each node in the corrected virtual steam pipeline network, determined based on the operating parameters, is used to correct the condensation of the corrected virtual steam pipeline network to obtain the target virtual steam pipeline network. Finally, the obtained target virtual steam pipeline network is used for steam simulation. In this way, on the one hand, a preprocessing mechanism can be provided to accurately obtain the flow direction of each pipe in the pipeline network before formal solution, so as to eliminate the uncertainty of flow direction caused by numerical error, thereby improving the stability and robustness of subsequent calculations; on the other hand, by correcting the condensation error, the implicitly increased condensate in actual operation can be reasonably distributed to each section of the pipeline network, making up for the gap between the ideal model and the actual working conditions, and improving the accuracy of thermal simulation results; furthermore, through the synergistic optimization of flow direction judgment and condensation correction, it can be ensured that the simulation model can still maintain good convergence performance and engineering credibility under complex multi-source, ring topology and dynamic working conditions.

[0077] In some embodiments, the steam pipeline adjustment module 11 may specifically include:

[0078] The steam output setting unit is used to identify the terminal node of the current steam pipeline network and set the steam output of the terminal node to a preset steam output to obtain an initial virtual steam pipeline network.

[0079] The steam source adjustment unit is used to calculate the increase in steam output of the initial virtual steam network and allocate the increase in steam output to the steam source nodes of the initial virtual steam network based on a preset allocation ratio to obtain the adjusted virtual steam network.

[0080] In some embodiments, the pipeline flow direction correction module 12 may specifically include:

[0081] The pre-solution unit is used to perform hydraulic calculations on the preset first mathematical model through the adjusted virtual steam pipe network to obtain the hydraulic pre-solution results corresponding to the adjusted virtual steam pipe network; the hydraulic pre-solution results include the mass distribution, pressure distribution and flow distribution corresponding to the virtual steam pipe network;

[0082] The pipeline flow direction adjustment unit is used to determine the current pipeline flow direction in the adjusted virtual steam network based on the mass distribution, the pressure distribution, and the flow rate distribution, and to adjust the upstream and downstream relationships of the adjusted virtual steam network based on the pipeline flow direction to obtain the virtual steam network to be mapped.

[0083] In some embodiments, the steam pipeline correction module 13 may specifically include:

[0084] The mapping unit is used to map the pre-solved result of the virtual steam network to be mapped to a preset fine grid through a preset numerical interpolation algorithm to obtain the virtual steam network to be corrected; the preset fine grid is a grid with a higher grid resolution than the current grid of the virtual steam network to be corrected.

[0085] In some embodiments, the steam pipeline correction module 13 may specifically include:

[0086] The flow statistics unit is used to count the current total input flow and total output flow of the virtual steam network to be corrected.

[0087] The first pipeline network determination unit is used to determine the virtual steam pipeline network to be corrected as the corrected virtual steam pipeline network if the total input flow rate is not less than the total output flow rate.

[0088] The second pipeline determination unit is used to determine the leakage amount based on the total input flow and the total output flow if the total input flow is less than the total output flow, and to compress the total output flow to be equal to the total input flow based on the leakage amount, so as to obtain a corrected virtual steam pipeline.

[0089] In some embodiments, the condensation correction module 14 may specifically include:

[0090] The operating parameter determination unit is used to spatially discretize the preset second mathematical model based on the preset numerical interpolation algorithm to obtain the discretized mathematical model, and to iteratively solve the discretized mathematical model through the corrected virtual steam network to obtain the operating parameters of each node in the corrected virtual steam network; the preset second mathematical model is a model that introduces an equivalent leakage coefficient into the preset first mathematical model; the equivalent leakage coefficient is a coefficient characterizing the intensity of condensation mass loss.

[0091] In some embodiments, the condensation correction module 14 may specifically include:

[0092] A condensation distribution determination unit is used to determine the condensation distribution of each node in the corrected virtual steam network based on the operating parameters.

[0093] The condensation correction unit is used to reduce the condensation of each node based on the steam consumption ratio of each node if the condensation distribution is negative, until the total input flow rate of the corrected virtual steam network is equal to the total output flow rate, so as to obtain the target virtual steam network, and perform steam simulation based on the target virtual steam network.

[0094] Furthermore, embodiments of this application also disclose an electronic device, Figure 3 This is a structural diagram of an electronic device 20 according to an exemplary embodiment. The content of the diagram should not be construed as limiting the scope of this application.

[0095] Figure 3 This is a schematic diagram of the structure of an electronic device 20 provided in an embodiment of this application. Specifically, the electronic device 20 may include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 stores a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the steam simulation method based on flow direction pre-solution disclosed in any of the foregoing embodiments. Alternatively, the electronic device 20 in this embodiment may specifically be an electronic computer.

[0096] In this embodiment, the power supply 23 is used to provide operating voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and external devices, and the communication protocol it follows can be any communication protocol applicable to the technical solution of this application, and is not specifically limited here; the input / output interface 25 is used to acquire external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs, and is not specifically limited here.

[0097] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk or optical disk, etc. The resources stored thereon can include operating system 221, computer program 222, etc., and the storage method can be temporary storage or permanent storage.

[0098] The operating system 221 is used to manage and control the various hardware devices on the electronic device 20 and the computer program 222, which may be Windows Server, Netware, Unix, Linux, etc. In addition to including a computer program capable of performing the steam simulation method based on flow direction pre-solution disclosed in any of the foregoing embodiments, the computer program 222 may further include computer programs capable of performing other specific tasks.

[0099] Furthermore, this application also discloses a computer-readable storage medium for storing a computer program; wherein, when the computer program is executed by a processor, it implements the aforementioned steam simulation method based on flow direction pre-solution. Specific steps of this method can be found in the corresponding content disclosed in the foregoing embodiments, and will not be repeated here.

[0100] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the apparatus disclosed in the embodiments, since it corresponds to the method disclosed in the embodiments, the description is relatively simple; relevant parts can be referred to in the method section.

[0101] Those skilled in the art will further recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0102] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein can be implemented directly by hardware, a software module executed by a processor, or a combination of both. The software module can be located in random access memory (RAM), main memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art.

[0103] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0104] The technical solutions provided in this application have been described in detail above. Specific examples have been used to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the methods and core ideas of this application. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of this application. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A steam simulation method based on pre-solved flow direction, characterized in that, include: Set the steam output of the terminal node in the current steam network to obtain the corresponding initial virtual steam network, and calculate the increase in steam output of the initial virtual steam network to allocate the increase in steam output to the steam source of the initial virtual steam network to obtain the adjusted virtual steam network. The adjusted virtual steam network is used to perform hydraulic calculations on the preset first mathematical model, and the pipe flow direction of the adjusted virtual steam network is corrected based on the obtained hydraulic calculation results to obtain the virtual steam network to be mapped. The hydraulic pre-solution results of the virtual steam network to be mapped are mapped to a preset fine mesh to obtain the virtual steam network to be corrected, and the leakage of the virtual steam network to be corrected is determined so as to correct the virtual steam network to be corrected based on the leakage to obtain the corrected virtual steam network. The preset second mathematical model is solved based on the corrected virtual steam network to obtain the operating parameters of each node in the corrected virtual steam network. The condensation distribution of each node in the corrected virtual steam network is determined based on the operating parameters to correct the condensation of the corrected virtual steam network, so as to obtain the target virtual steam network. Steam simulation is then performed based on the target virtual steam network.

2. The steam simulation method based on flow direction pre-solution according to claim 1, characterized in that, The process of setting the steam output of the terminal nodes in the current steam network to obtain a corresponding initial virtual steam network, and calculating the increase in steam output of the initial virtual steam network to allocate the increase in steam output to the steam sources of the initial virtual steam network, thereby obtaining an adjusted virtual steam network, includes: Identify the terminal nodes of the current steam network and set the steam output of the terminal nodes to a preset steam output to obtain an initial virtual steam network; The increase in steam output of the initial virtual steam network is statistically analyzed, and the increase in steam output is allocated to the steam source nodes of the initial virtual steam network based on a preset allocation ratio to obtain the adjusted virtual steam network.

3. The steam simulation method based on flow direction pre-solution according to claim 1, characterized in that, The process of performing hydraulic calculations on a preset first mathematical model using the adjusted virtual steam network, and then correcting the pipe flow direction of the adjusted virtual steam network using the obtained hydraulic calculation results to obtain the virtual steam network to be mapped, includes: The hydraulic calculations of the preset first mathematical model are pre-solved using the adjusted virtual steam network to obtain the hydraulic pre-solution results corresponding to the adjusted virtual steam network; the hydraulic pre-solution results include the mass distribution, pressure distribution, and flow distribution corresponding to the virtual steam network; Based on the mass distribution, pressure distribution, and flow distribution, the current pipeline flow direction in the adjusted virtual steam network is determined, and the upstream and downstream relationships of the adjusted virtual steam network are adjusted based on the pipeline flow direction to obtain the virtual steam network to be mapped.

4. The steam simulation method based on flow direction pre-solution according to claim 1, characterized in that, The step of mapping the hydraulic pre-solution results of the virtual steam network to be mapped to a preset fine mesh to obtain the virtual steam network to be corrected includes: The pre-solved results of the virtual steam pipeline network to be mapped are mapped to a preset fine grid using a preset numerical interpolation algorithm to obtain the virtual steam pipeline network to be corrected; the preset fine grid is a grid with a higher resolution than the current grid of the virtual steam pipeline network to be corrected.

5. The steam simulation method based on flow direction pre-solution according to claim 1, characterized in that, The step of determining the leakage amount of the virtual steam network to be corrected, and correcting the virtual steam network to be corrected based on the leakage amount to obtain a corrected virtual steam network, includes: Calculate the current total input flow and total output flow of the virtual steam network to be corrected; If the total input flow rate is not less than the total output flow rate, then the virtual steam network to be corrected is taken as the corrected virtual steam network; If the total input flow is less than the total output flow, the leakage is determined based on the total input flow and the total output flow, and the total output flow is compressed to be equal to the total input flow based on the leakage to obtain the corrected virtual steam network.

6. The steam simulation method based on flow direction pre-solution according to claim 4, characterized in that, The step of solving the preset second mathematical model based on the corrected virtual steam pipeline network to obtain the operating parameters of each node in the corrected virtual steam pipeline network includes: The preset second mathematical model is spatially discretized based on the preset numerical interpolation algorithm to obtain the discretized mathematical model. The discretized mathematical model is then iteratively solved using the corrected virtual steam network to obtain the operating parameters of each node in the corrected virtual steam network. The preset second mathematical model is a model that introduces an equivalent leakage coefficient into the preset first mathematical model. The equivalent leakage coefficient is a coefficient that characterizes the intensity of condensation mass loss.

7. The steam simulation method based on flow direction pre-solution according to any one of claims 1 to 6, characterized in that, The step of correcting the condensation of the corrected virtual steam network by using the condensation distribution of each node in the corrected virtual steam network determined based on the operating parameters to obtain the target virtual steam network, and performing steam simulation based on the target virtual steam network, includes: The condensation distribution of each node in the corrected virtual steam network is determined based on the operating parameters. If the condensation distribution is negative, the condensation of each node is reduced based on the steam consumption ratio of each node until the total input flow of the corrected virtual steam network is equal to the total output flow, so as to obtain the target virtual steam network, and steam simulation is performed based on the target virtual steam network.

8. A steam simulation device based on pre-solved flow direction, characterized in that, include: The steam pipeline network adjustment module is used to set the steam output of the end nodes in the current steam pipeline network to obtain the corresponding initial virtual steam pipeline network, and to count the increase in the steam output of the initial virtual steam pipeline network, so as to allocate the increase in the steam output to the steam source of the initial virtual steam pipeline network to obtain the adjusted virtual steam pipeline network. The pipeline flow direction correction module is used to perform hydraulic calculation pre-solution on the preset first mathematical model through the adjusted virtual steam pipeline network, and correct the pipeline flow direction of the adjusted virtual steam pipeline network through the obtained hydraulic pre-solution results to obtain the virtual steam pipeline network to be mapped. The steam pipeline network correction module is used to map the hydraulic pre-solution results of the virtual steam pipeline network to be mapped to a preset fine mesh to obtain the virtual steam pipeline network to be corrected, and to determine the leakage of the virtual steam pipeline network to be corrected, so as to correct the virtual steam pipeline network to be corrected based on the leakage, and obtain the corrected virtual steam pipeline network. The condensation correction module is used to solve a preset second mathematical model based on the corrected virtual steam network to obtain the operating parameters of each node in the corrected virtual steam network, and to correct the condensation of the corrected virtual steam network using the condensation distribution of each node in the corrected virtual steam network determined based on the operating parameters, so as to obtain a target virtual steam network, and to perform steam simulation based on the target virtual steam network.

9. An electronic device, characterized in that, include: Memory, used to store computer programs; A processor for executing the computer program to implement the steam simulation method based on flow direction pre-solution as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, Used to store a computer program, wherein the computer program, when executed by a processor, implements the steam simulation method based on flow direction pre-solution as described in any one of claims 1 to 7.

Images