Steam network condensate distribution prediction method, device, equipment and medium

Abstract

The present invention discloses a method, device, equipment, and medium for predicting condensate distribution in a steam pipeline network, relating to the field of smart city operation technology. The method comprises: establishing a condensate distribution prediction model; the condensate distribution prediction model includes a thermal runaway condensation calculation model and an instantaneous condensation calculation model; inputting gas source parameters, pipeline parameters, and steam state parameters at each node in the pipeline network into the condensate distribution prediction model; calculating the thermal runaway condensation amount using the thermal runaway condensation calculation model, and calculating the instantaneous condensation amount using the instantaneous condensation calculation model; and outputting a condensate distribution prediction result based on the sum of the thermal runaway condensation amount and the instantaneous condensation amount. This method can achieve real-time calculation within seconds, quickly and accurately predicting pipeline condensate distribution, reducing average error, and is low-cost and highly adaptable, providing data support for steam pipeline network operation optimization.

Description

Technical Field

[0001] The present invention relates to the field of smart city operation technology, and in particular to a method, device, equipment and medium for predicting condensate distribution in a steam network. Background Art

[0002] Steam pipeline networks are crucial heat transport systems in industrial production, and their operational efficiency directly impacts energy utilization and production costs. During steam transport, heat loss and pressure fluctuations cause steam to continuously condense, producing condensate. This accumulation of condensate can lead to problems such as water hammer, pipeline corrosion, and equipment damage, seriously impacting the safe operation and thermal efficiency of the pipeline network.

[0003] In existing technical solutions, steam pipe networks generally lack effective condensate distribution prediction methods, relying primarily on empirical judgment or periodic manual discharge. This results in an inability to accurately understand the condensate generation situation in each section of the pipe network. Furthermore, parameter coupling effects are often ignored, leading to large prediction errors. This ultimately results in either excessive condensate discharge, resulting in steam waste, or insufficient condensate discharge, resulting in safety hazards.

[0004] It can be seen that how to accurately predict the distribution of condensate in the steam pipe network is a technical problem that people in this field urgently need to solve. Summary of the Invention

[0005] The purpose of the present invention is to provide a method, device, equipment and medium for predicting the distribution of condensate in a steam pipeline network, which can quickly and accurately predict the distribution of condensate in the pipeline network in real time, with low cost and strong adaptability, and provide data support for the optimization of steam pipeline network operations.

[0006] In order to solve the above technical problems, the present invention provides a method for predicting condensate distribution in a steam network, comprising:

[0007] Establishing a condensate distribution prediction model; the condensate distribution prediction model includes a thermal runaway condensation amount calculation model and an instantaneous condensation amount calculation model;

[0008] Inputting gas source parameters, pipeline parameters and steam state parameters of each node of the pipeline network into the condensate distribution prediction model;

[0009] Calculating the thermal runaway condensation amount using the thermal runaway condensation amount calculation model, and calculating the instantaneous condensation amount using the instantaneous condensation amount calculation model;

[0010] According to the sum of the thermal runaway condensation amount and the instantaneous condensation amount, a condensation water distribution prediction result is output.

[0011] In a first aspect, in the steam network condensate distribution prediction method provided by the present invention, calculating the amount of thermal runaway condensation includes:

[0012] Calculate the heat loss per unit pipe length including insulation parameters;

[0013] The amount of heat runaway condensation is calculated based on the phase change heat, the pipe length, and the heat loss per unit pipe length including the insulation parameters.

[0014] On the other hand, in the steam network condensate distribution prediction method provided by the present invention, calculating the heat loss per unit pipe length including the insulation parameters includes:

[0015] Calculate the heat loss per unit pipe length including insulation parameters based on steam temperature, ambient temperature, thermal conductivity of insulation material, pipe outer diameter, insulation layer outer diameter, and ambient convection coefficient.

[0016] On the other hand, in the steam network condensate distribution prediction method provided by the present invention, calculating the instantaneous condensation amount includes:

[0017] Calculating the two-phase pressure drop and the flow pressure drop, and obtaining the total pressure drop according to the sum of the two-phase pressure drop and the flow pressure drop;

[0018] The instantaneous condensation amount is calculated based on the total pressure drop, the steam dryness change value, the steam quality, the pipeline input pressure and the condensation gain factor.

[0019] On the other hand, in the steam network condensate distribution prediction method provided by the present invention, calculating the two-phase pressure drop includes:

[0020] Calculate the coefficient of friction;

[0021] Calculate the pipeline elevation difference based on the pipeline length and pipeline inclination angle;

[0022] The two-phase pressure drop is calculated based on the friction coefficient, the pipeline elevation difference, the pipeline length, the pipeline outer diameter, the fluid density, and the average flow velocity of the fluid.

[0023] On the other hand, in the steam network condensate distribution prediction method provided by the present invention, calculating the friction coefficient includes:

[0024] Calculate the Re number based on the fluid density, average flow velocity of the fluid, outer diameter of the pipe and dynamic viscosity of the fluid;

[0025] The friction coefficient is calculated by an implicit equation based on the Re number, the outer diameter of the pipe and the roughness of the pipe;

[0026] While calculating the friction coefficient, the Newton iteration method is used to gradually approximate the roots of the implicit equation through local linearity, and the friction coefficient is updated in each iteration.

[0027] On the other hand, in the steam network condensate distribution prediction method provided by the present invention, calculating the friction coefficient includes:

[0028] Calculate the Re number based on the fluid density, average flow velocity of the fluid, outer diameter of the pipe and dynamic viscosity of the fluid;

[0029] Calculating a first intermediate variable based on the Re number and the pipeline roughness; and calculating a second intermediate variable based on the Re number;

[0030] The friction coefficient is calculated using the indicated equation based on the Re number, the first intermediate variable, and the second intermediate variable.

[0031] In order to solve the above technical problems, the present invention further provides a steam network condensate distribution prediction device, comprising:

[0032] A model building module is used to establish a condensed water distribution prediction model; the condensed water distribution prediction model includes a thermal runaway condensation amount calculation model and an instantaneous condensation amount calculation model;

[0033] A data input module, configured to input basic pipe network data including gas source parameters and pipeline parameters, as well as steam state parameters of each node in the pipe network, into the condensate distribution prediction model;

[0034] a condensation amount calculation module, configured to calculate the thermal runaway condensation amount using the thermal runaway condensation amount calculation model, and calculate the instantaneous condensation amount using the instantaneous condensation amount calculation model;

[0035] The result output module is used to output the condensation water distribution prediction result according to the sum of the thermal runaway condensation amount and the instantaneous condensation amount.

[0036] In order to solve the above technical problems, the present invention further provides an electronic device, comprising:

[0037] Memory for storing computer programs;

[0038] A processor is used to implement the steps of the above-mentioned steam network condensate distribution prediction method when executing the computer program.

[0039] In order to solve the above technical problems, the present invention also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the above-mentioned steam network condensate distribution prediction method are implemented.

[0040] It can be seen from the above technical solution that the method for predicting the condensate distribution in a steam pipeline network provided by the present invention includes: establishing a condensate distribution prediction model; the condensate distribution prediction model includes a thermal runaway condensation amount calculation model and an instantaneous condensation amount calculation model; inputting gas source parameters, pipeline parameters and steam state parameters of each node of the pipeline network into the condensate distribution prediction model; calculating the thermal runaway condensation amount using the thermal runaway condensation amount calculation model, and calculating the instantaneous condensation amount using the instantaneous condensation amount calculation model; outputting the condensate distribution prediction result based on the sum of the thermal runaway condensation amount and the instantaneous condensation amount.

[0041] The beneficial effect of the present invention is that the above-mentioned steam pipeline condensate distribution prediction method provided by the present invention takes into account multiple condensation mechanisms, establishes a condensate distribution prediction model between gas source parameters, pipeline parameters, steam state parameters and condensate generation, and uses the thermal runaway condensation amount calculation model to calculate the thermal runaway condensation amount, and uses the instantaneous condensation amount calculation model to calculate the instantaneous condensation amount. According to the sum of the thermal runaway condensation amount and the instantaneous condensation amount, real-time calculation in seconds can be achieved, and the condensate distribution of the pipeline network can be predicted quickly and accurately, the average error is reduced, the cost is low, and the adaptability is strong, providing data support for the optimization of steam pipeline operation.

[0042] In addition, the present invention also provides a corresponding steam pipeline condensate distribution prediction device, electronic device and computer-readable storage medium for the steam pipeline condensate distribution prediction method, which has the same or corresponding technical features as the above-mentioned steam pipeline condensate distribution prediction method and has the same effect as above. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] In order to more clearly illustrate the embodiments of the present invention, the following is a brief introduction to the drawings required for use in the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0044] Figure 1 A flow chart of a method for predicting condensate distribution in a steam network provided by an embodiment of the present invention;

[0045] Figure 2 A schematic diagram of the structure of a steam network condensate distribution prediction device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0046] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making any creative efforts shall fall within the scope of protection of the present invention.

[0047] It should be noted that, in the description of the present invention, the terms "comprises," "includes," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a process, method, article, or apparatus comprising a series of elements includes not only those elements but also other elements not explicitly listed, or elements inherent to such process, method, article, or apparatus. The terms "first," "second," etc., in the present invention are used to distinguish similar objects, and are not used to describe a particular order or precedence.

[0048] In order to enable those skilled in the art to better understand the solutions of the present invention, the present invention is further described in detail below with reference to the accompanying drawings and specific implementation methods.

[0049] In conjunction with the specific application environment architecture or specific hardware architecture on which the execution of the steam network condensate distribution prediction method depends, the specific application environment architecture or specific hardware architecture is described herein.

[0050] An embodiment of the present invention provides a method for predicting condensate distribution in a steam network. The method is described in detail in conjunction with the execution flow of the method for predicting condensate distribution in a steam network. Figure 1 The flow chart of the method for predicting the distribution of condensed water in a steam network provided by the embodiment of the present invention is as follows: Figure 1 As shown, the method includes:

[0051] S101. Establish a condensation water distribution prediction model; the condensation water distribution prediction model includes a thermal runaway condensation amount calculation model and an instantaneous condensation amount calculation model.

[0052] It should be noted that the condensate distribution prediction model established by the present invention is a model for predicting the distribution of condensate in a steam pipe network. In the condensate distribution prediction model, the calculation of thermal runaway condensation and instantaneous condensation are parallel core calculation items. The condensate distribution prediction model can be composed of a thermal runaway condensation calculation model and an instantaneous condensation calculation model. Among them, the thermal runaway condensation refers to the amount of condensate generated due to heat dissipation of the pipe wall; the thermal runaway condensation calculation model is a model for quantifying the amount of condensate generated due to heat dissipation of the pipe wall. The instantaneous condensation refers to the amount of condensate generated during pressure changes (especially during sudden drops); the instantaneous condensation calculation model is a model for quantifying the amount of condensate generated during pressure changes (especially during sudden drops).

[0053] S102. Input gas source parameters, pipeline parameters, and steam state parameters of each node in the pipeline network into a condensate distribution prediction model.

[0054] It should be noted that, in the present invention, gas source parameters refer to a series of quantitative indicators of the physical, chemical, and thermodynamic properties of the gas source in the steam pipeline network. Gas source parameters may include phase change heat (i.e., the collective term for the heat absorbed / released during the phase change of a substance), ambient temperature, ambient convection coefficient, fluid density, average fluid flow rate, fluid dynamic viscosity, etc. Pipeline parameters refer to a series of quantitative indicators that describe the physical characteristics, geometric dimensions, operating status, and transportation capacity of the steam pipeline. Pipeline parameters may include pipeline length, pipeline outer diameter, pipeline inclination angle, insulation layer outer diameter, insulation material thermal conductivity, etc. The steam state parameters of each node in the pipeline network refer to the quantitative indicators of the physical state of steam at each node in the steam pipeline network. The steam state parameters of each node in the pipeline network may include steam temperature, steam quality, steam dryness change value, etc.

[0055] S103. Calculate the thermal runaway condensation amount using a thermal runaway condensation amount calculation model, and calculate the instantaneous condensation amount using an instantaneous condensation amount calculation model.

[0056] In practice, the heat runaway condensation amount calculation model can be used to calculate the heat runaway condensation amount in real time. , the instantaneous condensation amount can be calculated in real time using the model .

[0057] S104: Output the condensate distribution prediction result based on the sum of the thermal runaway condensation amount and the instantaneous condensation amount.

[0058] In practice, the condensate distribution prediction results can be .

[0059] In the above-mentioned steam pipeline condensate distribution prediction method provided by the embodiment of the present invention, multiple condensation mechanisms are taken into account, and a condensate distribution prediction model is established between gas source parameters, pipeline parameters, steam state parameters and condensate generation. The thermal runaway condensation amount is calculated using a thermal runaway condensation amount calculation model, and the instantaneous condensation amount is calculated using an instantaneous condensation amount calculation model. According to the sum of the thermal runaway condensation amount and the instantaneous condensation amount, real-time calculation in seconds can be achieved, and the condensate distribution in the pipeline network can be predicted quickly and accurately, reducing the average error. It has low cost and strong adaptability, providing data support for steam pipeline operation optimization.

[0060] In practical applications, before executing step S103, the pipe network topology can be verified for compliance, including isolated nodes, pipe materials, and connectivity. Specifically, the closure of pipe segment connections (no isolated nodes) can be checked. This closure requires that all pipe segments in the pipe system are properly connected, forming a complete network with no disconnections or disconnections. A node is a connection point between pipe segments, and an isolated node is a node that is not connected to any other pipe segments. Such isolated points cannot exist in steam pipes, as they will affect the overall system functionality. This verification ensures that the pipe system's structure meets design requirements, providing an accurate basis for subsequent fluid flow analysis and pressure calculations. Furthermore, material parameter matching can be performed. For example, a hydrophobic cotton database can be established to capture the variation in its thermal conductivity under different temperature conditions and present it as a curve. This helps understand the changes in its thermal insulation performance under different operating temperatures and provides a basis for optimal selection of hydrophobic cotton. Using computer programs or related algorithms, based on input basic information such as pipe diameter, the recommended insulation thickness can be quickly and accurately determined, improving work efficiency and accuracy.

[0061] After executing step S104, a condensate distribution prediction report can be generated based on the condensate distribution prediction results, including the condensate amount of each pipe section, recommended drainage strategy, etc., to achieve a closed loop from prediction to application.

[0062] Furthermore, in a specific implementation, in the above-mentioned steam network condensate distribution prediction method provided in an embodiment of the present invention, step S103 calculates the thermal runaway condensation amount, which may specifically include: calculating the heat loss per unit pipe length including insulation parameters; calculating the thermal runaway condensation amount based on phase change heat, pipeline length, and heat loss per unit pipe length including insulation parameters.

[0063] In practice, when calculating the amount of heat runaway condensation using the heat runaway condensation calculation model, the heat loss per unit pipe length including the insulation parameters can be calculated first. Then, based on the phase change heat, the pipe length, and the heat loss per unit pipe length including the insulation parameters, the heat runaway condensation can be calculated using the following formula (1):

[0064] ; (1)

[0065] in, Indicates the amount of thermal runaway condensation, which refers to the mass of matter that changes from gas to liquid per hour, and the unit is usually kilograms per hour (kg / h); Indicates the heat loss per unit pipe length including insulation parameters; Indicates the pipe length, that is, the actual length of the pipe through which the fluid flows; Represents phase change heat, which refers to the phase change process of matter ( ) can be obtained by looking up the steam table.

[0066] Furthermore, in the specific implementation, in the above steps, the heat loss per unit pipe length including insulation parameters is calculated, which can specifically include: calculating the heat loss per unit pipe length including insulation parameters based on steam temperature, ambient temperature, thermal conductivity of insulation material, outer diameter of pipe, outer diameter of insulation layer and ambient convection coefficient.

[0067] In practice, based on the heat transfer principle of the cylindrical wall, according to the steam temperature, ambient temperature, thermal conductivity of the insulation material, outer diameter of the pipe, outer diameter of the insulation layer and ambient convection coefficient, the heat loss per unit pipe length including insulation parameters can be calculated using the following formula (2):

[0068] ; (2)

[0069] in, Indicates steam temperature (°C), Indicates the ambient temperature (°C), Indicates the thermal conductivity of hydrophobic cotton (W / m·K), with a reference value of 0.035-0.045; Indicates the outer diameter of the pipe (m); represents the outer diameter of the insulation layer (m), where is the thickness of the insulation layer; Indicates the ambient convection coefficient (W / m²·K), with a value of 5-25 (depending on wind speed).

[0070] It should be noted that the above formula (1) combined with formula (2) can be used as a calculation model for the amount of thermal runaway condensation.

[0071] The present invention can use the energy conservation equation to To verify, the energy conservation direction can be expressed as the following formula (3):

[0072] ; (3)

[0073] in, 、 They represent the specific enthalpy of the fluid at the inlet and outlet of the pipeline, respectively, in joules per kilogram (J / kg). Specific enthalpy is the enthalpy per unit mass of the working fluid, a thermodynamic state parameter that reflects the energy state of the fluid; 、 The flow rates of the fluid at the inlet and outlet of the pipe are respectively in meters per second (m / s), reflecting the flow speed of the fluid at different positions; 、 are the heights of the fluid at the inlet and outlet of the pipe, respectively, in meters, used to determine the vertical position of the fluid; is the potential energy term, reflecting the change in pipe height ( represents the acceleration due to gravity, is the node elevation); is the mass flow rate (kg / s).

[0074] Furthermore, in specific implementation, in the above-mentioned steam network condensate distribution prediction method provided in an embodiment of the present invention, step S103 calculates the instantaneous condensation amount, which can specifically include: calculating the two-phase pressure drop and the flow pressure drop, and obtaining the total pressure drop based on the sum of the two-phase pressure drop and the flow pressure drop; calculating the instantaneous condensation amount based on the total pressure drop, steam dryness change value, steam quality, pipeline input pressure and condensation gain factor.

[0075] In practice, when calculating the instantaneous condensation using the instantaneous condensation calculation model, the two-phase pressure drop can be calculated first. and flow pressure drop , two-phase voltage drop is specific to two-phase flow systems and includes friction pressure drop as a component; flow pressure drop , refers to the total pressure drop in fluid motion (including friction, gravity, acceleration, and form resistance), where represents the density of the fluid, represents the acceleration due to gravity; Indicates the pipeline elevation difference and flow pressure drop Based on the principles of fluid statics, the flow pressure drop is proportional to the fluid density, gravitational acceleration, and the elevation difference between the two ends of the pipe, and is related to the inclination angle. For example, in a vertical pipe transporting liquid, the greater the liquid density and the greater the height difference of the pipe, the greater the flow pressure drop. and flow pressure drop The total pressure drop is obtained by adding The total pressure drop is composed of two parts: the two-phase pressure drop and the flow pressure drop. It refers to the total pressure difference between two points along the flow direction when the fluid flows in the pipeline. It reflects the total pressure loss caused by the fluid overcoming the pipeline resistance and the change in potential energy.

[0076] Then, according to the total pressure drop , steam dryness change, steam quality, pipeline input pressure and condensation gain factor, the following formula (4) can be used to calculate the instantaneous condensation amount:

[0077] ; (4)

[0078] in, Indicates the instantaneous condensation amount, which refers to the mass of the substance that changes from gas to liquid during the instantaneous condensation process; Indicates the steam mass, that is, the mass of steam involved in the instantaneous condensation process, in kilograms (kg); 、 Respectively represent the dryness of steam at the process inlet and outlet; dryness is a parameter that measures the proportion of steam mass in wet steam, with a value range of 0-1, where 0 represents all liquid and 1 represents all gas; is the change in steam dryness; Indicates the inclination gain factor, which is a correction factor based on the Re number (i.e., Reynolds number) or inclination angle and has no fixed unit; Indicates the input pressure in Pascal (Pa) or Megapascal (MPa).

[0079] It should be noted that the instantaneous condensation amount can include the tilt angle correction, indicating that the calculation of the condensation amount takes into account the influence of the tilt angle factor. This formula (4) is established based on the relevant theories or experiences that consider the influence of the tilt angle. Some parameters in the formula (such as etc.) may be determined or modified by considering the effect of the tilt angle on the system (such as the gravity component along the tilt direction affecting fluid flow, energy exchange, etc.), so that the calculated It can more accurately reflect the instantaneous setting amount under conditions with tilt angles.

[0080] It should be noted that when the pipe is tilted at an angle If the angle is less than or equal to the set angle (such as 10°), you can directly use Perform instantaneous condensation; when the pipe is tilted at an angle When the angle is greater than the set angle (such as 10°), the tilting gain factor of the activated tilting section can be calculated using the following formula (5):

[0081] ; (5)

[0082] in, It represents the tilt gain factor of the tilt section. It is based on the basic coefficient (coefficient related to system characteristics) based on the pipe inclination angle Formula (5) shows that the tilt angle The bigger, relatively The more it increases, the more it reflects the enhancement effect of the inclination angle on instantaneous setting. ), the calculation of instantaneous coagulation is corrected. In the calculation of instantaneous coagulation, the Replace the original coefficient , thereby more accurately considering the influence of the inclination angle on the instantaneous setting process, so that the calculation results can better reflect the instantaneous setting situation under actual working conditions, and improve the accuracy of related engineering calculations and analysis.

[0083] Furthermore, in a specific implementation, in the above steps, the two-phase pressure drop is calculated, including: calculating the friction coefficient; calculating the pipeline elevation difference based on the pipeline length and the pipeline inclination angle; calculating the two-phase pressure drop based on the friction coefficient, the pipeline elevation difference, the pipeline length, the pipeline outer diameter, the fluid density and the average flow rate of the fluid.

[0084] In practice, when calculating the two-phase pressure drop, the friction coefficient can be calculated first; then the pipeline elevation difference can be calculated based on the pipeline length and pipeline inclination angle using the following formula (6):

[0085] ; (6)

[0086] in, Indicates the pipeline elevation difference; Indicates the pipeline inclination angle, which refers to the angle between the pipeline centerline and the horizontal direction.

[0087] Then, the two-phase pressure drop can be calculated using the following formula (7) based on the friction coefficient, pipeline elevation difference, pipeline length, pipeline outer diameter, fluid density, and average flow velocity of the fluid:

[0088] ; (7)

[0089] in, Indicates the two-phase voltage drop; It represents the friction coefficient, which is related to factors such as the roughness of the inner wall of the pipe and the flow state of the fluid (laminar flow or turbulent flow). It reflects the magnitude of the frictional resistance between the fluid and the pipe wall. Indicates the density of the fluid, the mass per unit volume of fluid; Indicates the average flow velocity of the fluid, that is, the average speed at which the fluid flows in the pipe; represents the acceleration due to gravity; It indicates the pipeline elevation difference, which refers to the height difference between the two ends of the pipeline in the vertical direction, and is used to measure the position difference of the pipeline in the vertical direction.

[0090] It should be noted that the present invention takes into account that when the pipeline is tilted, the fluid is affected not only by the friction resistance with the pipe wall, but also by the force component of gravity along the pipeline direction. This term corrects the additional effect of gravity on the pressure drop caused by the inclination of the pipe, so that the equation can more accurately calculate the friction pressure drop of the fluid in the inclined pipe.

[0091] When a fluid flows in a pipe, the friction coefficient is an important parameter that measures the frictional resistance between the fluid and the pipe wall. It plays a key role in calculating pressure drop, flow rate, etc. In practice, in the above steps, the friction coefficient can be calculated using two specific implementation methods:

[0092] In a first specific implementation, calculating the friction coefficient may specifically include: calculating the Re number based on the fluid density, the average flow velocity of the fluid, the outer diameter of the pipe and the dynamic viscosity of the fluid; calculating the friction coefficient through an implicit equation based on the Re number, the outer diameter of the pipe and the roughness of the pipe; while calculating the friction coefficient, using the Newton iteration method to gradually approximate the roots of the implicit equation through local linearity, and updating the friction coefficient in each iteration.

[0093] In practice, when calculating the friction coefficient, the Re number can be calculated using the following formula (8) based on the fluid density, the average flow velocity of the fluid, the outer diameter of the pipe and the dynamic viscosity of the fluid:

[0094] ; (8)

[0095] in, The above formula (8) is used to determine the flow state of the fluid (laminar or turbulent), and comprehensively reflects the influence of factors such as fluid velocity, density, viscosity and pipe size on the flow.

[0096] Then, according to the Re number, the outer diameter of the pipe and pipe roughness, the friction coefficient is calculated using the implicit equation (9):

[0097] ; (9)

[0098] in, represents the pipe roughness, in meters. The above formula (9) is an implicit equation, that is, the friction coefficient Since the coefficient appears on both sides of the equation, it cannot be solved directly. Therefore, the present invention utilizes the Newton iteration method to iteratively calculate the friction coefficient. Specifically, this method uses local linear stepwise approximation to the roots of the implicit equation, updating the friction coefficient with each iteration. This method is suitable for high-precision industrial-grade standard calculations. Alternatively, software tools (such as Excel iterative calculations or specialized fluid calculation software) can be used for numerical solutions, enabling accurate analysis and design of pipes and conduits in fluid transport systems.

[0099] In a second specific implementation, calculating the friction coefficient may specifically include: calculating the Re number based on the fluid density, the average flow velocity of the fluid, the outer diameter of the pipe and the dynamic viscosity of the fluid; calculating the first intermediate variable based on the Re number and the roughness of the pipe; and calculating the second intermediate variable based on the Re number; calculating the friction coefficient through the displayed equation based on the Re number, the first intermediate variable and the second intermediate variable.

[0100] In practice, in the process of calculating the friction coefficient, the Re number can be first calculated based on the fluid density, the average flow velocity of the fluid, the outer diameter of the pipe and the dynamic viscosity of the fluid using the above formula (8); then, based on the Re number and the pipe roughness, the first intermediate variable can be calculated using the following formula (10); and based on the Re number, the second intermediate variable can be calculated using the following formula (11); finally, based on the Re number, the first intermediate variable and the second intermediate variable, the friction coefficient can be calculated using the following formula (12) by displaying the equation:

[0101] ; (10)

[0102] ; (11)

[0103] ; (12)

[0104] in, represents the first intermediate variable, The above formula (12) avoids implicit solution. It is an explicit equation. Given the parameters such as Re, the friction coefficient can be directly calculated. , no iterative solution is required, and compared with the implicit equation, the error is less than , capable of rapid calculation, real-time simulation or programming implementation.

[0105] In practical applications, if high accuracy is required, formula (9) can be used to calculate the friction coefficient; if the calculation result needs to be obtained quickly, formula (12) can be used to calculate the friction coefficient.

[0106] It should be noted that the above formulas (5)-(12) can be used as a calculation model for instantaneous coagulation.

[0107] In the present invention, the characteristic comparison between the thermal runaway condensation amount and the instantaneous condensation amount can be seen in Table 1.

[0108] Table 1 Comparison of characteristics between thermal runaway condensation and instantaneous condensation

[0109]

[0110] In actual operating conditions, the formula for calculating thermal runaway condensation can be used. By introducing empirical correction factors for insulation degradation and dynamic effects, the formula for calculating instantaneous condensation can be used. The total condensate is the sum of thermal runaway condensation and instantaneous condensation.

[0111] In the method for predicting condensate distribution in a steam network provided by the present invention, some parameters can be subjected to sensitivity analysis, as shown in Table 2 for details.

[0112] Table 2 Parameter sensitivity analysis

[0113]

[0114] The present invention can conduct an experimental design on the steam network condensate distribution prediction method according to the content of Table 3, and Table 4 is a comparison of the results.

[0115] Table 3 Experimental design

[0116]

[0117] Table 4 Comparison of results

[0118]

[0119] The present invention establishes a mathematical model based on the basic principles of thermodynamics, takes into account multiple condensation mechanisms, and makes the prediction results more accurate. The multi-parameter coupling model reduces the average error from 22.6% to 4.3%.

[0120] The present invention can achieve predictions without the need for additional hardware equipment, utilizing Supervisory Control and Data Acquisition (SCADA) system data. It can also quickly respond to changes in pipe network topology and adjustments to operating parameters, support dynamic predictions, and provide a scientific basis for optimizing steam trap layout and formulating operating strategies. This significantly improves steam system energy efficiency and supports complex operating conditions such as inclined pipes and insulation deterioration.

[0121] Through the description of the above implementation methods, those skilled in the art can clearly understand that the method according to the above embodiment can be implemented by means of software plus the necessary general hardware platform, and of course it can also be implemented by hardware, but in many cases the former is a better implementation method.

[0122] An embodiment of the present application also provides a device for predicting condensate distribution in a steam network. Figure 2 This is a schematic diagram of the structure of the steam network condensate distribution prediction device provided by the embodiment of the present invention. This embodiment is based on the perspective of functional modules, such as Figure 2 As shown, the device includes:

[0123] The model building module 10 is used to establish a condensed water distribution prediction model; the condensed water distribution prediction model includes a thermal runaway condensation amount calculation model and an instantaneous condensation amount calculation model;

[0124] The data input module 11 is used to input the basic data of the pipeline network including gas source parameters and pipeline parameters, as well as the steam state parameters of each node of the pipeline network into the condensate distribution prediction model;

[0125] The condensation amount calculation module 12 is used to calculate the thermal runaway condensation amount using the thermal runaway condensation amount calculation model and calculate the instantaneous condensation amount using the instantaneous condensation amount calculation model;

[0126] The result output module 13 is used to output the condensation water distribution prediction result according to the sum of the thermal runaway condensation amount and the instantaneous condensation amount.

[0127] In the above-mentioned steam network condensate distribution prediction device provided in an embodiment of the present invention, a condensate distribution prediction model can be established between the gas source parameters, pipeline parameters, steam state parameters and condensate generation through the interaction of the above-mentioned four modules. The thermal runaway condensation amount is calculated using the thermal runaway condensation amount calculation model, and the instantaneous condensation amount is calculated using the instantaneous condensation amount calculation model. According to the sum of the thermal runaway condensation amount and the instantaneous condensation amount, real-time calculation in seconds can be achieved, and the condensate distribution in the network can be predicted quickly and accurately, the average error can be reduced, the cost is low, and the adaptability is strong, providing data support for the optimization of steam network operations.

[0128] Since the embodiments of the steam network condensate distribution prediction device correspond to the embodiments of the steam network condensate distribution prediction method, the description of the features of the corresponding embodiments of the steam network condensate distribution prediction device can be found in the description of the corresponding embodiments of the steam network condensate distribution prediction method, and will not be repeated here. The embodiments have the same beneficial effects as the aforementioned steam network condensate distribution prediction method.

[0129] An embodiment of the present application also provides an electronic device, including a memory and a processor, wherein the memory stores a computer program, and the processor is configured to run the computer program to execute the steps in any of the above-mentioned steam network condensate distribution prediction method embodiments.

[0130] An embodiment of the present invention further provides a computer-readable storage medium storing a computer program, wherein the computer program is configured to execute the steps of any of the above-mentioned steam network condensate distribution prediction method embodiments when running.

[0131] In an exemplary embodiment, the computer-readable storage medium may include, but is not limited to, various media that can store computer programs, such as a USB flash drive, a read-only memory (ROM), a random access memory (RAM), a mobile hard disk, a magnetic disk, or an optical disk.

[0132] An embodiment of the present invention further provides a computer program product, which includes a computer program. When the computer program is executed by a processor, the steps in any of the above-mentioned steam network condensate distribution prediction method embodiments are implemented.

[0133] An embodiment of the present invention also provides another computer program product, including a non-volatile computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, it implements the steps of any of the above-mentioned steam network condensate distribution prediction method embodiments.

[0134] Professionals may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the above description has generally described the components and steps of each example according to their functions. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians may use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the present invention.

[0135] The above is a detailed introduction to the steam network condensate distribution prediction method, device, equipment and medium provided by the present invention. Specific examples are used herein to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only intended to help understand the method and core concept of the present invention. It should be pointed out that for ordinary technicians in this technical field, without departing from the principles of the present invention, several improvements and modifications can be made to the present invention, and these improvements and modifications also fall within the scope of protection of the present invention.

Claims

1. A method for predicting condensate distribution in a steam network, characterized in that: include: Establishing a condensate distribution prediction model; the condensate distribution prediction model includes a thermal runaway condensation amount calculation model and an instantaneous condensation amount calculation model; Inputting gas source parameters, pipeline parameters and steam state parameters of each node of the pipeline network into the condensate distribution prediction model; Calculating the amount of thermal runaway condensation using the thermal runaway condensation amount calculation model; The instantaneous condensation amount calculation model is used to calculate the two-phase pressure drop and the flow pressure drop, and the total pressure drop is obtained based on the sum of the two-phase pressure drop and the flow pressure drop. The instantaneous condensation amount is calculated using the following formula based on the total pressure drop, the steam dryness change, the steam quality, the pipeline input pressure, and the condensation gain factor: ； in, Indicates the instantaneous condensation amount, which refers to the mass of the substance that changes from gas to liquid during the instantaneous condensation process; Indicates steam quality; 、 Respectively represent the dryness of steam at the process inlet and outlet; is the change in steam dryness; It represents the inclination gain factor, which is a correction factor based on the inclination angle; Indicates input pressure, represents the total pressure drop; According to the sum of the thermal runaway condensation amount and the instantaneous condensation amount, a condensation water distribution prediction result is output.

2. The method for predicting condensate distribution in a steam network according to claim 1, characterized in that: Calculate the amount of heat runaway condensation, including: Calculate the heat loss per unit pipe length including insulation parameters; The amount of heat runaway condensation is calculated based on the phase change heat, the pipe length, and the heat loss per unit pipe length including the insulation parameters.

3. The method for predicting condensate distribution in a steam network according to claim 2, wherein: Calculates heat loss per unit pipe length with insulation parameters, including: Calculate the heat loss per unit pipe length including insulation parameters based on steam temperature, ambient temperature, thermal conductivity of insulation material, pipe outer diameter, insulation layer outer diameter, and ambient convection coefficient.

4. The method for predicting condensate distribution in a steam network according to claim 1, wherein: Calculates two-phase voltage drop, including: Calculate the coefficient of friction; Calculate the pipeline elevation difference based on the pipeline length and pipeline inclination angle; The two-phase pressure drop is calculated based on the friction coefficient, the pipeline elevation difference, the pipeline length, the pipeline outer diameter, the fluid density, and the average flow velocity of the fluid.

5. The method for predicting condensate distribution in a steam network according to claim 4, characterized in that: Calculates the coefficient of friction, including: Calculate the Re number based on the fluid density, average flow velocity of the fluid, outer diameter of the pipe and dynamic viscosity of the fluid; The friction coefficient is calculated by an implicit equation based on the Re number, the outer diameter of the pipe and the roughness of the pipe; While calculating the friction coefficient, the Newton iteration method is used to gradually approximate the roots of the implicit equation through local linearity, and the friction coefficient is updated in each iteration.

6. The method for predicting condensate distribution in a steam network according to claim 4, characterized in that: Calculates the coefficient of friction, including: Calculate the Re number based on the fluid density, average flow velocity of the fluid, outer diameter of the pipe and dynamic viscosity of the fluid; Calculating a first intermediate variable based on the Re number and the pipeline roughness; and calculating a second intermediate variable based on the Re number; The friction coefficient is calculated using the indicated equation based on the Re number, the first intermediate variable, and the second intermediate variable.

7. A steam network condensate distribution prediction device, characterized in that: include: A model building module is used to establish a condensed water distribution prediction model; the condensed water distribution prediction model includes a thermal runaway condensation amount calculation model and an instantaneous condensation amount calculation model; A data input module, configured to input basic pipe network data including gas source parameters and pipeline parameters, as well as steam state parameters of each node in the pipe network, into the condensate distribution prediction model; The condensation amount calculation module is used to calculate the thermal runaway condensation amount using the thermal runaway condensation amount calculation model; calculate the two-phase pressure drop and the flow pressure drop using the instantaneous condensation amount calculation model, and obtain the total pressure drop based on the sum of the two-phase pressure drop and the flow pressure drop; calculate the instantaneous condensation amount using the following formula based on the total pressure drop, the steam quality change value, the steam quality, the pipeline input pressure, and the condensation gain factor: ； in, Indicates the instantaneous condensation amount, which refers to the mass of the substance that changes from gas to liquid during the instantaneous condensation process; Indicates steam quality; 、 Respectively represent the dryness of steam at the process inlet and outlet; is the change in steam dryness; It represents the inclination gain factor, which is a correction factor based on the inclination angle; Indicates input pressure, represents the total pressure drop; The result output module is used to output the condensation water distribution prediction result according to the sum of the thermal runaway condensation amount and the instantaneous condensation amount.

8. An electronic device, characterized in that: include: memory for storing computer programs; A processor is configured to implement the steps of the method for predicting condensate distribution in a steam network as described in any one of claims 1 to 6 when executing the computer program.

9. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method for predicting condensate distribution in a steam network according to any one of claims 1 to 6 are implemented.

Images