Assessment and Prediction of Complex Industrial Steam Network Operation by Combined Thermo-Hydrodynamic Modeling

: Steam network operation stability and reliability is vital for any industrial branch. A combined steam network model comprising a balance and a coupled thermo-hydrodynamic model, including seasonal variations impact and system speciﬁcities, is presented. A balance model can readily be used by a reﬁnery’s operators. The thermo-hydrodynamic model identiﬁes system bottlenecks and cold spots and evaluates proposed operation and investment measures including heat loss reduction. A three-pressure levels reﬁnery steam network served for model testing and validation. Balance model results reveal signiﬁcant misbalance in steam production and consumption, reaching 30.5% in the low-pressure steam system, and heat balance di ﬀ erences in the range of 9.2% to 29.5% on individual pressure levels, attributable both to ﬂow measurement accuracy issues and to heat losses. The thermo-hydrodynamic model results di ﬀ er from the measured steam parameters by less than 5% (temperature) and by less than 4% (pressure), respectively, with the estimated operational insulation heat conductivity exceeding 0.08 W / m / K. Its comparison with that of 0.03 W / m / K for dry insulation material yields the need for pipelines re-insulation and a partial revamp of the steam network. The model is su ﬃ ciently general for any type of industry, pursuing the goal of cleaner and energy-e ﬃ cient steam transport and consumption.


Introduction
Water steam belongs to the most extensively used energetic media. Already in 2006 the share of energy consumption in the form of water steam represented up to 40% of all energetic media used in industry. In the refining industry, up to 28% of total energy consumption belongs to water steam production [1]. A system with a central steam source (combined heat and power (CHP) unit) with decentralized secondary ones is preferred, with steam being transported to its consumers through complicated pipe networks with total length reaching up several tens to hundreds of kilometers within one plant [2]. With such transport distances, it is difficult to keep the steam network in a stable state and to secure desired steam quality to all consumers under all conditions. Fluctuations of steam parameters may affect production units in the refinery. Stable and effective operation of steam networks also leads to a reduction of steam losses to the environment and reduced operational costs for the whole enterprise [3].
Several papers have been published regarding the design of different utility systems and their optimization. Nishio et al. [4] used a simple linear programming (LP) algorithm for turbine allocation.

1.
Despite the documented extensive research in steam networks modeling, most of presented models deal with steam network optimization just as a part of a bigger optimization scheme [31,33]. 2.
Detailed system specificities and topology are seldom considered which might lead to system operation analysis oversimplification [17,18]. 3.
Seasonal ambient air temperature variations impact on system operation is rarely considered which might prevent the performed study to reveal the season-related system bottlenecks [21]. 4.
Most models analyze the steam system under one certain or a few selected loads (transported steam amounts and transport directions) which does not allow for a rigorous assessment of problematic operation situations and troublesome spots in the steam network [34,36]. 5.
Though system operation change/improvement projects (re-insulation, new pipelines commissioning, impact of new production units on system operation) are commonly considered and tested, the testing procedure frequently does not consider all possible network operational states, or, at least, all those states that were recorded in the past [14,19]. This in turn limits the testing ability to identify system bottlenecks possibly resulting from improvement project implementation.
The aim of the presented paper was to focus directly on the issues discussed above, namely closely on the steam system itself, incorporating its specificities and the impact of seasonal variations, thus providing an alternative and more realistic point of view on the steam system's operation limits. More specifically, the main aim can be split into the following: • Creation of a mathematical model that would be able to detect and to predict troublesome spots formation within the steam networks with increased possibility of cold spots or bottlenecks occurrence. • Its testing on a real complex steam system of a refinery; its use to reveal particular operation problems and limits in the existing system. • The use of the developed model to predict the impact of planned investments in the refinery on the operation stability of the steam system, considering seasonal system operation changes.

•
Proposal and simulation of the system changes for its operation improvement.
Firstly, mass and heat balances were performed for each individual steam pipeline followed by a simple balance model calculating steam flow velocities in each pipe segment to assess the possibility of "cold spot" or "bottleneck" formation. Secondly, a combined thermo-hydrodynamic model based on hydrodynamic and heat-transfer equations was created that allows comparing calculated values of steam temperature and pressure to the real field data. This model, after validation, was used for assessment of proposed changes within the steam network.

Problem Description
The SLOVNAFT company is a complex oil refinery designed and built in the 1960s. Gradual development of equipment and technology as well as improved heat recovery meant that the operational parameters of steam networks are significantly different than they were designed for originally. The risk of cold spots occurrence in the steam mains increases especially due to the steadily rising share of steam produced in local steam sources. Non-optimal operation of steam networks leads to significant economic losses. Other facts contribute to this state as well: • The "safe" layout of the refinery-all operational units in the refinery are located relatively far apart. Thus, there is a need to transport steam over long distances, leading to increased heat losses to the surroundings. • Several steam-drives at the SLOVNAFT refinery have been switched to electro-drives, except important devices, that are still being driven by steam turbines as a safety measure. Switching between electro-and steam-drives is important not just from the economical point of view but allows the steam balance to be influenced. With too many devices switched to electro-drives (which is economically effective), steam supply from a CHP unit would drop significantly and result in steam condensation in pipelines. Therefore, the economics is not the only decision-making factor. Operational safety is important as well. • Seasonal impact on operational state of steam networks. Summer months in particular are problematic due to excessive steam supply from secondary energy sources (SEC, all steam supplying units except CHP). The aforementioned switching between electro-and steam-drives reduces the possibilities of its utilization by production units. Steam excess is disposed of, by its condensation in air condensers, or, in worst situations, by its exhaust to the atmosphere.
Operation of steam networks is regulated by an internal directive, whose most important goals are: • Keeping balance between steam production and its consumption.

•
Maximal utilization of steam produced in CHP unit and secondary heat sources.

•
Maximization of waste heat utilization in the form of low-pressure steam (LPS) produced at production units.

•
Preserving steam quality parameters in desired intervals.

•
Elimination of cold spots in steam networks. • Utilization of electro-drives.

•
Minimization of steam let-downs in CHP unit and other production units.
The occurrence of cold spots has the potential to exhibit damaging effect on several steam utilizing devices such as steam turbines. Condensate droplets are carried by flowing steam into turbines and due to their high velocity, turbine blades could get seriously damaged. Therefore, increased attention is paid to particular pipeline sections, deemed as potential cold spots, mainly based on historical data and previous experiences. Steam networks are managed the way the temperature is always higher than a certain value at these locations.
This method is very simple and effective, but has a major disadvantage: This way, new possible cold spots or bottlenecks cannot be detected in advance. This is a motivation for development of a mathematical model of steam networks able to overcome it.
The applied study method includes the following: 1. Creation of mathematical models of steam networks and their application in the refinery on all pressure levels based on available operation data, visual controls and plant layouts.

2.
Mass and heat balances setup in order to determine current operational state of steam networks in the refinery.

3.
Confronting the simulation results with measured values of operational data to identify incorrect local measurement tools (steam pressure, temperature and mass flow).

4.
Assessment of troublesome spots formation risk in terms of low or very high steam flow velocity within individual sections of steam networks.

5.
Reviewing the impact of proposed projects within steam networks on its stability by a mathematical model expressed as steam temperature differences at individual steam-delivery points. 6.
Suitable measures proposal aimed at an improvement in the operational state of steam networks and the associated operating costs reduction.

Operational Data Analysis
Three steam networks of different pressure level (low-pressure steam (LS); middle-pressure steam (MS); high-pressure steam (HS)) are operated in the refinery. The total length of all steam pipes reaches up to 60 km and there are 27, 20 and 18 active production units respectively connected to each steam pressure level. The average share of steam amount supplied by local steam sources to the total steam supply reaches up to 60% in the case of the HS network. Steam flows and steam amounts being consumed or supplied in standard operational modes by individual production units as well as maximal steam flows through individual pipeline sections needed to be determined using past operational data. It was necessary to choose data corresponding to standard operational mode and sort/exclude data corresponding to planned or unplanned technical shutdowns and process data measurement failure.
Analyzed operational data included steam mass flows, temperatures and pressures. Based on values of steam temperature and pressure, the quality of steam (i.e., superheated/saturated/wet steam) can be determined. Moreover, by combination of steam pressure and the Antoine equation, the boiling point and superheating can be determined.
Finally, potentially incorrect values were identified and corrected or removed by detailed analysis of an extensive set of field data. Errors in the set of values can in our experience be caused by measurement failure, signal transmission error, improper placement of the measurement tool, a systematically different method of measuring, or exceeding the measuring scale.
To identify potentially incorrect values of raw process data, both simple logic and authors' industrial experience obtained in previous practice were exploited. Negative process values (of mass flows, temperature and pressure) were treated as follows: The decision for a correction method was always made according to a timespan of such data failure. In case of a timespan shorter than 24 h, negative values have been replaced by the weighted arithmetic mean of the values observed during 24 h before the measurement failure. In a case that such failure persisted for more than 24 h, steam import or export were estimated using unit feed and specific steam consumption or production obtained from the unit s past operation. The same approach was used in case of data failure accompanied by an error messages "No good data for calculation", that were probably the most often occurring issues we had to deal with. Rarely, in the case of steam flows, we observed unit of measurement switch from t/h to kg/h and vice versa. Such values were manually transformed to the original unit of measurement, i.e., they were divided or multiplied by 1000.

Balance Model
The balance model was designed for real time monitoring purposes. In fact, the balance model is as simple as possible, thus shortening the calculation time and reducing hardware requirements and enabling its use by steam network operators.
The input data for the balance model comprise the following: • Steam flows at individual production units (import, export or both, since some production units may act like steam importer and exporter at the same steam pressure level). • Steam temperature and pressure at individual production units.

•
Steam pipeline characteristics (lengths and diameters of individual pipe sections).
The input data were represented by the time-weighted hourly average values of measured field data measured during the whole year 2017. Thus, 8760 data sets containing 106, 68 and 81 measured values corresponding to high-pressure, middle-pressure and low-pressure steam (HPS, MPS and LPS), respectively, were processed. Prior to simulation, all data were analyzed and filtered if necessary (Section 3.1) to obtain reliable results.
The simulation schemes were prepared according to pipeline dimensions and existing technical drawings of the plant to preserve its layout as much as possible. A simplified scheme of HPS network is shown in Figure 1. Based on the output data the overall technical state of steam networks as well as the economic impact of this state was evaluated. It was possible to assess the existence of cold spots as well as bottlenecks based on steam flow velocities w calculated in individual pipe sections by Equation (1), where represents the steam mass flow, v the specific volume of steam and Din is the internal diameter of that pipe section. (1)

Combined Thermo-Hydrodynamic Model
This model was created in a commercial chemical modeling software environment, able to perform simulations of single pipes as well as complicated pipe networks. Other important features are: • Simulation of single-or multiphase flow pressure drop with possible selection of calculation method.

•
Heat-transfer calculation including outside weather conditions or physical properties of pipes.

•
Calculation of profiles of different physical properties in pipe sections.
It is possible to perform simulations from single pipe segments to complex pipe networks. Each pipe segment can be divided into multiple sections with different properties (diameter, material, roughness, …). Number and types of fittings can be specified as well for pressure drop calculation.
This model is built to calculate temperature and pressure of outlet streams, while their mass flows are defined by the user. The inlet streams must be specified by combination of at least three input parameters: inlet pressure, temperature and molar/mass flow.

Pressure Drop Calculation
The used software works with many fluid flow pressure drop correlations. There are many empirical as well as mechanistic correlations, developed due to the calculation of pressure drop in wellbores or crude oil transport. Under common operation conditions, only superheated steam should be present in steam pipelines, but based on experience, saturated steam or even two-phase flow occurs on spots with lower steam flows. Therefore, we chose Beggs and Brill (1979) [43], which is the most widely used and most reliable fluid flow pressure drop correlation in industry suitable for two-phase flow simulations.
All pipe segments require definition of their properties, which are necessary for pressure drop calculation: Specific enthalpy of steam for the purpose of heat balance was calculated by the IAPWS-IF97 method [42] that is considered to be a standard for steam enthalpy calculation published in 1997 by the International Association for the Property of Water and Steam. The specific steam volume was determined by method IAPWS-IF97 as well and it was used for steam flow velocity calculation.
The output data of the balance model are the following: • Steam mass balance difference (i.e., the difference between steam mass supplied and consumed within the steam network). • Steam heat balance difference (i.e., the difference between the heat supplied and consumed within the steam network in the form of steam).

•
Steam flow velocities in individual pipe segments calculated according to steam mass flows in these pipe segments.
Based on the output data the overall technical state of steam networks as well as the economic impact of this state was evaluated. It was possible to assess the existence of cold spots as well as bottlenecks based on steam flow velocities w calculated in individual pipe sections by Equation (1), where . m represents the steam mass flow, v the specific volume of steam and D in is the internal diameter of that pipe section.

Combined Thermo-Hydrodynamic Model
This model was created in a commercial chemical modeling software environment, able to perform simulations of single pipes as well as complicated pipe networks. Other important features are: • Simulation of single-or multiphase flow pressure drop with possible selection of calculation method. It is possible to perform simulations from single pipe segments to complex pipe networks. Each pipe segment can be divided into multiple sections with different properties (diameter, material, roughness, . . . ). Number and types of fittings can be specified as well for pressure drop calculation.
This model is built to calculate temperature and pressure of outlet streams, while their mass flows are defined by the user. The inlet streams must be specified by combination of at least three input parameters: inlet pressure, temperature and molar/mass flow. The used software works with many fluid flow pressure drop correlations. There are many empirical as well as mechanistic correlations, developed due to the calculation of pressure drop in wellbores or crude oil transport. Under common operation conditions, only superheated steam should be present in steam pipelines, but based on experience, saturated steam or even two-phase flow occurs on spots with lower steam flows. Therefore, we chose Beggs and Brill (1979) [43], which is the most widely used and most reliable fluid flow pressure drop correlation in industry suitable for two-phase flow simulations.
All pipe segments require definition of their properties, which are necessary for pressure drop calculation: Inner diameter and outer diameter or tube wall thickness; • Construction material (and its roughness); • Elevation (increase or decrease of geodetic height).
For a single-phase fluid flow a simple pressure drop calculation was used. In such a case, pressure drop is calculated using Bernoulli's Equation (2).
where z is the geodetic height, g is gravitational acceleration, w is the fluid flow velocity, P is absolute pressure, ρ is the fluid density, ε dis is the energy dissipated due to fluid flow and (−ε pump ) is energy added by a pump. Γ is dimensionless parameter with values equal to 0.5 or 1 when the flow regime is laminar or turbulent, respectively. Index 1 refers to inlet and index 2 corresponds to outlet from given pipe segment. Geodetic heights difference is neglected here. The inlet, w 1 and outlet, w 2 steam flow velocities are equal when the pipe diameter and cross-section area remain constant throughout the whole segment. This yields from the continuity Equation (3).
Based on mentioned assumptions, for a single pipe segment without any pump, Equation (2) can be re-written as follows (4): The dissipation of mechanical energy can be easily calculated as the sum of energy dissipation due to friction, ε dis,f and local energy dissipation, ε dis,L (5).
where L is the length of pipe segment, D is inner pipe diameter, w is the average fluid flow velocity in given segment, λ is the fluid flow friction factor and ξ is the coefficient of local dissipation of mechanic energy. The value of friction factor, λ, is a function of the Reynolds number, Re and the relative roughness of pipe wall, n. The Equations (6)-(8), used for calculation of the fluid flow friction factor are listed in Table 1 together with conditions of their use. All these Equations are applicable for fluid flow in pipes of circular cross-section. Alternatively, a novel approach presented in [44], unifying the friction factor calculation for all flow regimes as a result of heuristic and evolutionary algorithms merging, could be used.
Round (1980) [47] As mentioned earlier, the presence of two-phase flow is common in such sections, where steam flows are relatively low. For such situations, Beggs and Brill correlation was used.
The flow regime is determined based on the Froude number, Equation (9) and following dimensionless variables, Equations (10)-(13) [43]: in Table 2. Where N Fr is the Froude number, v m is the flow velocity of two-phase mixture, D is pipe inner diameter and λ L is dimensionless parameter given by the ratio of the flow velocity of liquid phase to the flow velocity of two-phase mixture, Equation (14).
Values of λ L and the Froude number are determining the flow regime. The segregated flow is determined if these conditions, Equations (15) and (16), are met: For transition flow regime, condition represented by Equation (17), must be met: For intermittent flow regime conditions represented by Equations (18) or (19), must be met: Finally, the distributed flow regime is considered if one of conditions, Equations (20) or (21) is met: Liquid holdup Liquid holdup is expressed as Equation (22) and is calculated for all flow regimes, except transition.
where H L(0) , defined by Equations (23) and (24), is the liquid holdup which would exist at the same conditions in a horizontal pipe.
The values of parameters a, b and c are specific for each flow regime and are listed in Table 3. The holdup correcting factor Ψ takes the pipe inclination effect into account and is given by Equation (25): where ϕ is the actual angle of the pipe from the horizontal. The parameter C in Equation (25) is calculated by Equation (26) and its value must be higher than or equal to zero.
The values of parameters d , e, f and j are specific for each flow regime in Table 4. The transition flow regime is a region between the segregated and intermittent flow regimes. The liquid holdup for transition flow regime, H L(transition) , is as a superposition of their liquid holdup factors expressed by Equation (27): where the parameters A and B are determined by Equations (28) and (29): Pressure drop due to friction Frictional pressure drop can be calculated by Equation (30), where f tp is the friction factor calculated by Equation (31): f tp = f n e V (31) f n is the friction factor for smooth pipe and is determined from the smooth pipe curve of the Moody diagram depending on Reynolds number, Re. For laminar flow (Re ≤ 2000), the Hagen-Poiseuille equation is used, Equation (6). The friction factor for turbulent flow is described well by the Blasius equation, (7). The parameter V in Equation (31) can be calculated by Equation (32) if the condition (33) is met. For other values, Equation (34) must be used.

Heat-Transfer Intensity Calculation
For heat loss calculation, there are three heat-transfer sub-models available in the software used.

Zero heat transfer; •
Simple heat transfer-this sub-model assumes, that the heat transfer is directly proportional to the temperature difference between the fluid and the environment, ∆T. The resulting heat loss, is expressed by Equation (35), with U referring to the overall heat-transfer coefficient: • Detailed heat transfer-a complex heat-transfer sub-model that allows specifying more detailed set of heat-transfer parameters to model the heat transfer between the fluid, the vessel and the environment due to conduction and convection.
Simulations were performed with detailed heat-transfer sub-model to achieve high accuracy of results. This sub-model considers all partial phenomena:

•
Heat transfer through natural (or forced) convection between the vessel fluid and the pipe wall; • Heat transfer through conduction through the wall and insulation layers; • Heat transfer through natural or forced convection from the insulation surface to the environment.
The total area, A tot , in the detailed and simple sub-model generally consists of the cylindrical area and the head areas of used device. Steam network was assumed to consist only of pipes, thus the area A refers only to their surface.
In addition to thermal insulation, conduction parameters for the pipe wall are also required: • Pipe wall thickness; • Specific heat capacity; • Mass density of construction material; • Thermal conductivity The value of overall heat-transfer coefficient, U, is continually updated in each time step. It is calculated with auto-selected correlations based on the system conditions. Three convection heat-transfer coefficients for heat transfer were considered: • From the vessel outer wall to ambient air; • From vessel contents to inner wall; • From vapor to liquid phase inside the vessel For heat transfer from outer wall to outside air, the correlation constants are automatically selected based on the Grashof, Gr f and Prandtl numbers, Pr f , defined by Equations (36) and (37), respectively.
where c p is the fluid (air) specific heat capacity, v is the fluid kinematic viscosity, κ is its thermal conductivity, g is the gravitational acceleration, β is the thermal expansion coefficient, ρ is the mass density and X is the characteristic vessel dimension. This depends on the vessel orientation, being its axial length in case of a vertical one, or the vessel diameter in case of a horizontal one. The product of these dimensionless numbers determines the type of convection regime. If Pr f ·Gr f < 10 9 , the flow regime is laminar and the heat-transfer coefficient, α, is defined by Equation (38).
where m = 1.4 and C is parameter depending on vessel orientation. Its value equals C = 1.42 Wm −2 K −1 for vertical vessels and C = 1.32 Wm −2 K −1 for horizontal ones. If the following hold true Pr f ·Gr f > 10 9 , the flow regime is turbulent and the heat-transfer coefficient is expressed as defined by Equation (39): where m = 1.4 and C = 1.31 Wm −2 K −1 for vertical vessel and C = 1.24 Wm −2 K −1 for horizontal ones. The overall heat-transfer coefficient, U, comprises thermal properties of all layers involved into heat transfer and is calculated based on partial heat-transfer coefficients, α. In general, the value of U related to inner wall surface is expressed by Equation (40). This Equation is applicable for number of layers equal to i = 1,2, . . . ,n.
where κ i is the thermal conductivity of i-th pipe wall layer and D LMi is its logarithmic mean diameter. Subscripts inn and out refer to inside and outside pipe wall, respectively.

Model Validation (Analysis of Actual Operation)
One of the issues solved was the mass balance difference. The result of the balance model was used as follows: Measured steam flows at each outlet stream were adjusted in the way, the mass balance difference was distributed through all outlet streams proportionally, i.e., based on the ratio of the steam flow to total steam supply. This modification is expressed as follows Equation (41) where the subscript f,i relates to the balanced mass flow of the i-th outlet stream and the subscript meas,i means the measured mass flow of the i-th outlet stream. The calculation steps are performed by an algorithm depicted in Figure 2.
used as follows: Measured steam flows at each outlet stream were adjusted in the way, the mass balance difference was distributed through all outlet streams proportionally, i.e., based on the ratio of the steam flow to total steam supply. This modification is expressed as follows Equation (41) and ensures an equal distribution of steam balance difference (eventual losses): where the subscript f,i relates to the balanced mass flow of the i-th outlet stream and the subscript meas,i means the measured mass flow of the i-th outlet stream. The calculation steps are performed by an algorithm depicted in Figure 2.
In step 0, the initial steam parameters are imported from the real measured data. Thermodynamic properties are determined. For simplification, steam parameters are assumed to be constant throughout individual pipe sections in this step.
Step 1 involves mass balance, which is crucial part for next calculations. Subsequently, steam flow velocity in each pipe section can be determined as well in step 2 (Equation (1)). Pressure drop in each section is calculated using combination of Equations (6) to (34).
Step 3 involves calculation of heat losses in each pipe section, Equations (35) to (40). This step also includes recalculation of overall heat-transfer coefficient, Uinn, in each iteration. From known values of heat losses, values of the steam temperature decrease can be obtained. This is important for recalculation of steam parameters in step 4. The objective function of the model is the sum of differences of pressure drops in two following iterations squared. The value of parameter k was determined based on previous experiences to obtain results with sufficient accuracy in acceptable calculation time.
Outputs from the thermo-hydrodynamic model are the following: • Temperatures of all inner and outlet streams; In step 0, the initial steam parameters are imported from the real measured data. Thermodynamic properties are determined. For simplification, steam parameters are assumed to be constant throughout individual pipe sections in this step.
Step 1 involves mass balance, which is crucial part for next calculations. Subsequently, steam flow velocity in each pipe section can be determined as well in step 2 (Equation (1)). Pressure drop in each section is calculated using combination of Equations (6)- (34).
Step 3 involves calculation of heat losses in each pipe section, Equations (35)- (40). This step also includes recalculation of overall heat-transfer coefficient, U inn , in each iteration. From known values of heat losses, values of the steam temperature decrease can be obtained. This is important for recalculation of steam parameters in step 4.
The objective function of the model is the sum of differences of pressure drops in two following iterations squared. The value of parameter k was determined based on previous experiences to obtain results with sufficient accuracy in acceptable calculation time.
Outputs from the thermo-hydrodynamic model are the following: Validation/verification of this model was done by comparison of calculated and measured values of temperature and pressure at the outlet streams. The average relative errors for i-th stream are calculated by Equation (42): where a = T,P; indices c and meas refer to calculated and measured values, respectively, N is yearly working time fund, i.e., count of data sets. Values of relative errors can also help us to identify possibly incorrect measurement tools. Calculated temperature and pressure drops are directly related to steam flow velocities. Elevated pressure drop indicates high steam flow and on the other hand, elevated temperature drop indicates low steam flow. Calculated steam flow velocities can be compared to those obtained by the balance model and assessed by temperature and pressure losses as well.

The Impact of Proposed Changes on Steam Network Operation
After model validation (analysis of actual operation), it was used for prediction of transported steam parameters after certain technological changes within the high-pressure steam network. By such simulation it is possible to assess impact of these changes and to identify potential problems. The following changes were considered: A.
Direct connection of units C6 and S2 by a new steam pipeline; B.
Commissioning of new production unit C10; C.
Commissioning new hydrogen production plant S5 thus replacing the old unit S1.
Project A is related only to change of topology, while projects B and C will affect the mass balance of steam networks as well by adding new steam suppliers/consumers. Building a connection between C6 and S2 can change steam flow direction and thus affect steam temperature in adjacent units.
The unit C10 should process part of sour gases from refinery and thus decrease the load of the existing unit sulphur-recovery unit S2 located in its proximity ( Figure 3). As a result, the HPS export from S2 will decrease by more than 11 t/h.

•
Pressures of all inner and outlet streams; • Temperature drops in each pipe segment; • Pressure drops in each pipe segment; • Steam flow velocities in each pipe segment.
Validation/verification of this model was done by comparison of calculated and measured values of temperature and pressure at the outlet streams. The average relative errors for i-th stream are calculated by Equation (42): where a = T,P; indices c and meas refer to calculated and measured values, respectively, N is yearly working time fund, i.e., count of data sets. Values of relative errors can also help us to identify possibly incorrect measurement tools.
Calculated temperature and pressure drops are directly related to steam flow velocities. Elevated pressure drop indicates high steam flow and on the other hand, elevated temperature drop indicates low steam flow. Calculated steam flow velocities can be compared to those obtained by the balance model and assessed by temperature and pressure losses as well.

The Impact of Proposed Changes on Steam Network Operation
After model validation (analysis of actual operation), it was used for prediction of transported steam parameters after certain technological changes within the high-pressure steam network. By such simulation it is possible to assess impact of these changes and to identify potential problems. The following changes were considered: A. Direct connection of units C6 and S2 by a new steam pipeline; B. Commissioning of new production unit C10; C. Commissioning new hydrogen production plant S5 thus replacing the old unit S1.
Project A is related only to change of topology, while projects B and C will affect the mass balance of steam networks as well by adding new steam suppliers/consumers. Building a connection between C6 and S2 can change steam flow direction and thus affect steam temperature in adjacent units.
The unit C10 should process part of sour gases from refinery and thus decrease the load of the existing unit sulphur-recovery unit S2 located in its proximity ( Figure 3). As a result, the HPS export from S2 will decrease by more than 11 t/h. New hydrogen production plant S5 is only in the design stage and should substitute the obsolete hydrogen production plant S1 with connection to the HPS network in its proximity (Figure 3). The HPS export should increase from current 11 t/h (from S1) up to 35 t/h (from S5) on average. However, S5 will be able to produce up to 54 t/h steam at its full load, which can be a critical amount due to the transport capacity of the steam pipeline in this section.
Each of these projects will have a significant impact on steam network's operation. At unchanged loads of other production units, it will be necessary to modify the steam supply from the CHP unit to maintain the steam mass balance. The effect of proposed changes on each steam network is summarized in Table 5. Simulation scheme of HPS network was modified according to information relating to this project as can be seen in Figure 3. Table 5. The effect of proposed changes A, B and C on the steam amount supplied to steam networks. LPS = low-pressure steam, MPD = middle-pressure steam, HPS = high-pressure steam.

Steam Network Steam Export Change, t/h
Data analysis proved that the operation of HPS network varies a lot with season as well as with different loads of other production units. Evidence of this fact is provided in Figure 4 showing a total HPS amount supplied from the CHP unit throughout the year 2017 varying from 0 to 100 t/h. Average ambient temperatures were selected according to local meteorological measurements and they are listed in Table 6. A representative set of values was selected, corresponding to situations, when ambient temperature was equal to the selected one. Table 6 also contains values of overall HPS supply and HPS supply from the CHP unit at that time. New hydrogen production plant S5 is only in the design stage and should substitute the obsolete hydrogen production plant S1 with connection to the HPS network in its proximity (Figure 3). The HPS export should increase from current 11 t/h (from S1) up to 35 t/h (from S5) on average. However, S5 will be able to produce up to 54 t/h steam at its full load, which can be a critical amount due to the transport capacity of the steam pipeline in this section.
Each of these projects will have a significant impact on steam network's operation. At unchanged loads of other production units, it will be necessary to modify the steam supply from the CHP unit to maintain the steam mass balance. The effect of proposed changes on each steam network is summarized in Table 5. Simulation scheme of HPS network was modified according to information relating to this project as can be seen in Figure 3. Table 5. The effect of proposed changes A, B and C on the steam amount supplied to steam networks. LPS = low-pressure steam, MPD = middle-pressure steam, HPS = high-pressure steam.

Steam Network Steam export Change, t/h
Data analysis proved that the operation of HPS network varies a lot with season as well as with different loads of other production units. Evidence of this fact is provided in Figure 4 showing a total HPS amount supplied from the CHP unit throughout the year 2017 varying from 0 to 100 t/h. Average ambient temperatures were selected according to local meteorological measurements and they are listed in Table 6. A representative set of values was selected, corresponding to situations, when ambient temperature was equal to the selected one. Table 6 also contains values of overall HPS supply and HPS supply from the CHP unit at that time.    Simulations were performed with two different values of the thermal conductivity of thermal insulation. The value of κi = 0.038 Wm −1 K −1 refers to a new insulation in proper technical state. On the other hand, the value of κi = 0.080 Wm −1 K −1 refers to real state of thermal insulation and it was determined based on results of analysis of current operation of steam network and visual inspection in SLOVNAFT, a.s.

Data Analysis
The HPS network in the SLOVNAFT refinery is the shortest one with total length of approx. 7 km and there are only 6 steam supplying and 12 steam consuming units connected to it. The share of the amount of HPS supplied from the CHP unit was only 39.52% during year 2017. This means, that this steam network is the least-dependent on the steam supply from the CHP unit. By contrast, the MPS network is the most dependent on the steam supply from the CHP unit with the share of steam amount almost 84% of total one. There is a significant effect of season on steam consumption as well as supply. Figure 5 shows a time-course of LPS supply and the variability of transported steam amounts. Less significant seasonal trends can be observed in the case of the MPS and HPS networks as well. Especially the LPS steam network suffers from exceeding amount of steam supply during summer due to high air temperatures and lower heat utilization. This is common issue in chemical industry, especially in refining and petrochemical industry. In such situations the waste heat-recovery methods should be utilized as much as possible/needed [48]. Moreover, in most critical situations, the exceeding steam is led to the air condenser to decrease its supply into steam network and keep the network balanced. Simulations were performed with two different values of the thermal conductivity of thermal insulation. The value of κi = 0.038 Wm −1 K −1 refers to a new insulation in proper technical state. On the other hand, the value of κi = 0.080 Wm −1 K −1 refers to real state of thermal insulation and it was determined based on results of analysis of current operation of steam network and visual inspection in SLOVNAFT, a.s.

Data Analysis
The HPS network in the SLOVNAFT refinery is the shortest one with total length of approx. 7 km and there are only 6 steam supplying and 12 steam consuming units connected to it. The share of the amount of HPS supplied from the CHP unit was only 39.52% during year 2017. This means, that this steam network is the least-dependent on the steam supply from the CHP unit. By contrast, the MPS network is the most dependent on the steam supply from the CHP unit with the share of steam amount almost 84% of total one. There is a significant effect of season on steam consumption as well as supply. Figure 5 shows a time-course of LPS supply and the variability of transported steam amounts. Less significant seasonal trends can be observed in the case of the MPS and HPS networks as well. Especially the LPS steam network suffers from exceeding amount of steam supply during summer due to high air temperatures and lower heat utilization. This is common issue in chemical industry, especially in refining and petrochemical industry. In such situations the waste heatrecovery methods should be utilized as much as possible/needed [48]. Moreover, in most critical situations, the exceeding steam is led to the air condenser to decrease its supply into steam network and keep the network balanced.

Balance Model
The results of the balance model characterize the overall technical state of steam networks in the refinery and the efficiency of their operation. The results of mass and heat balances of each steam pressure level network are listed in Table 7 as the average balance difference during the whole year 2017. The balance difference means the difference between the amount of supplied and consumed mass and heat within the steam network, respectively. Table 7. Results of mass and heat balances of individual steam pressure level networks.

Balance Model
The results of the balance model characterize the overall technical state of steam networks in the refinery and the efficiency of their operation. The results of mass and heat balances of each steam pressure level network are listed in Table 7 as the average balance difference during the whole year 2017. The balance difference means the difference between the amount of supplied and consumed mass and heat within the steam network, respectively. Graphical interpretation of mass and heat balance is depicted in Figures 6-8.    Graphical interpretation of mass and heat balance is depicted in Figures 6-8.   The average mass balance difference in 2017 in the LPS was 38.13 t/h, i.e., 30.5% of total steam amount supplied from the CHP unit and 55.8% of total steam amount supplied into this network, respectively. This phenomenon could be partially explained due to LPS consumption as steam tracing as well as increased heat losses during the winter months. This consumption is not controlled, and it is very difficult to estimate its amount. There is a significant drop on 15.5.2017, where the steam supply from the CHP dropped by approx. 34 t/h during one single day. However, the total steam measured consumption remained constant. From this reason, the LPS consumption as steam tracing was estimated as high as 34 t/h during the winter months.
Negative LPS mass balance in summer months is seen in Figure 6: when there is no other option to manipulate the LPS steam balance, the LPS is led to air cooled steam condensers. This consumption is not measured which is the explanation of the observed fact.
The mass and heat balance results of MPS network have similar trends as in case of LPS network. Increased balance difference during winter can be explained by increased heat losses to the surroundings and uncontrolled steam tracing. However, the difference between summer and winter is not as sensible as the difference at LPS network. The average mass balance difference is 21.52 t/h. It means the share of 27.70% and 21.16% on the steam supply from the CHP unit and the total steam supply, respectively. The trend of heat balance copies the trend of mass balance. The measurement of temperatures and pressures can be therefore considered as reliable.
The results of HPS mass and steam balance are the most satisfying. The average mass balance difference is only 6.5% of total HPS amount supplied into this network and 16.3% of total HPS amount supplied from the CHP unit, respectively. There are a few drops and peaks, which can be explained by unit shut downs or measurement failures. The biggest drop on 10.12.2017 is highlighted in Figure  8. Input data analysis showed, that during the period from 8.12 to 13.12.2017 all steam flow meters on the CHP unit battery limit showed constant value. The assumption that this is not a failure of the measuring tool, but rather a failure of signal transmission was later confirmed by consultation with energy managers.
Dividing average heat and mass balance difference, the value of theoretical enthalpy of loss of steam, ht, Equation (43) is obtained. Values of theoretical enthalpy calculated for each steam network compared with typical, h, minimal, hmin and maximal, hmax, values of enthalpy of steam at given pressure level defined by internal regulation are listed in Table 8. The average mass balance difference in 2017 in the LPS was 38.13 t/h, i.e., 30.5% of total steam amount supplied from the CHP unit and 55.8% of total steam amount supplied into this network, respectively. This phenomenon could be partially explained due to LPS consumption as steam tracing as well as increased heat losses during the winter months. This consumption is not controlled, and it is very difficult to estimate its amount. There is a significant drop on 15.5.2017, where the steam supply from the CHP dropped by approx. 34 t/h during one single day. However, the total steam measured consumption remained constant. From this reason, the LPS consumption as steam tracing was estimated as high as 34 t/h during the winter months.
Negative LPS mass balance in summer months is seen in Figure 6: when there is no other option to manipulate the LPS steam balance, the LPS is led to air cooled steam condensers. This consumption is not measured which is the explanation of the observed fact.
The mass and heat balance results of MPS network have similar trends as in case of LPS network. Increased balance difference during winter can be explained by increased heat losses to the surroundings and uncontrolled steam tracing. However, the difference between summer and winter is not as sensible as the difference at LPS network. The average mass balance difference is 21.52 t/h. It means the share of 27.70% and 21.16% on the steam supply from the CHP unit and the total steam supply, respectively. The trend of heat balance copies the trend of mass balance. The measurement of temperatures and pressures can be therefore considered as reliable.
The results of HPS mass and steam balance are the most satisfying. The average mass balance difference is only 6.5% of total HPS amount supplied into this network and 16.3% of total HPS amount supplied from the CHP unit, respectively. There are a few drops and peaks, which can be explained by unit shut downs or measurement failures. The biggest drop on 10.12.2017 is highlighted in Figure 8. Input data analysis showed, that during the period from 8.12 to 13.12.2017 all steam flow meters on the CHP unit battery limit showed constant value. The assumption that this is not a failure of the measuring tool, but rather a failure of signal transmission was later confirmed by consultation with energy managers.
Dividing average heat and mass balance difference, the value of theoretical enthalpy of loss of steam, h t , Equation (43) is obtained. Values of theoretical enthalpy calculated for each steam network compared with typical, h, minimal, h min and maximal, h max , values of enthalpy of steam at given pressure level defined by internal regulation are listed in Table 8. Values of theoretical enthalpy of LPS and MPS are inside or close to the interval of typical values defined by internal regulation. Moreover, there is a similar time-course of both mass and heat balances. This shows that temperature and pressure measurements in the LPS network are sufficiently accurate. The value of the theoretical steam enthalpy of HPS is inadequate, since it differs a lot from the typical value of HPS enthalpy. Assuming the enthalpy of 3.1 GJ/t and given mass balance difference, the overall heat losses due to material losses should reach 18.8 GJ/h. The rest of 7.8 GJ/h should be explained by heat loss to surrounding or incorrect measurement of temperature, pressure or both.
The balance model was designed to perform calculations of steam flow velocity in each pipe sections as well, thus aiming to detect the risk of troublesome spots formation within the steam networks. The upper steam flow velocity limit value of 30 m/s was set based on real experience, depicted in Figure 9, as well as in accordance with expert literature [49]. At the moment of highest steam pressure at unit S2 the excess steam was discharged through the pressure relief valve into the environment and the calculated exported steam flow velocity reached value 29.9 m/s as depicted in Figure 9. The following pressure drop resulted from switching a steam-drive into operation, thus reducing the steam amount supplied from the S2 unit into the HPS network.  Values of theoretical enthalpy of LPS and MPS are inside or close to the interval of typical values defined by internal regulation. Moreover, there is a similar time-course of both mass and heat balances. This shows that temperature and pressure measurements in the LPS network are sufficiently accurate. The value of the theoretical steam enthalpy of HPS is inadequate, since it differs a lot from the typical value of HPS enthalpy. Assuming the enthalpy of 3.1 GJ/t and given mass balance difference, the overall heat losses due to material losses should reach 18.8 GJ/h. The rest of 7.8 GJ/h should be explained by heat loss to surrounding or incorrect measurement of temperature, pressure or both.
The balance model was designed to perform calculations of steam flow velocity in each pipe sections as well, thus aiming to detect the risk of troublesome spots formation within the steam networks. The upper steam flow velocity limit value of 30 m/s was set based on real experience, depicted in Figure 9, as well as in accordance with expert literature [49]. At the moment of highest steam pressure at unit S2 the excess steam was discharged through the pressure relief valve into the environment and the calculated exported steam flow velocity reached value 29.9 m/s as depicted in Figure 9. The following pressure drop resulted from switching a steam-drive into operation, thus reducing the steam amount supplied from the S2 unit into the HPS network.

3.
The CHP unit export into MPS pipeline D-the same reason as indicated in point 2. Moreover, this pipeline supplies the MPS steam to an external consumer and must be able to transport significant amounts of steam in the heating period.

Thermo-Hydrodynamic Model
Calculated values of steam temperature and pressure were compared to the measured ones and differences as well as relative errors were evaluated. Steam flow velocities were used to determine potential bottlenecks or cold spots. Such identification is possible based on heat loss values as well. High specific heat loss indicates low steam flow rate and vice versa. There is a similar meaning for specific pressure drop. A high value (over 0.3 kPa/m) indicates high steam flow velocity and a possible bottleneck. By contrast, a low pressure drop indicates low steam flow velocity and a possible cold spot.
Results of simulations proved positive as well as negative temperature and pressure differences. The results can be considered reliable with a value of relative error under 5% and values over 10% were excluded and the measurement was marked as "incorrect". Values of average, minimal and maximal relative errors of temperature and pressure for each network are listed in Table 9. Please note that all values are absolute values, since the difference can be either positive or negative. According to the results in Table 7, the LPS network is in worst technical state, since the corresponding mass balance difference is the largest here. The impact of this finding on the results shown in Table 9 is the following: • A large balance difference distributed to all outlet streams could have caused significant relative errors of the thermo-hydrodynamic model.

•
The complexity of the LPS network and many interconnections between each pipeline may also have caused large errors. • Finally, the nominal values of temperature as well as pressure are lowest in the LPS network which means that the same absolute difference yields in higher relative error in LPS network than in the other ones.

The Impact of Proposed Changes on Steam Network Operation
The purpose of these simulations was to determine, how the proposed changes within steam network would affect the steam network operation and to detect potential operational troubles. The connection of units C6 and S2 (Project A) should solve the problem of comparatively cold HPS delivered to the C6 unit. The temperature at its battery limit is typically within the range 295 to 310 • C. HPS with temperature approx. 300 • C could already be a risk to safety and reliable operation of steam turbines located at this unit based on previous experiences. This is probably the most important parameter observed in these simulations. At present, both mentioned units are located at the ends of main steam pipelines. While unit S2 supplies a significant amount of steam, unit C6 is a less important steam consumer. Moreover, the distance from S2 to C6 is quite long, approx. 350 m. The combination of these two facts leads to low steam flow velocity in this section and accompanied increased heat Results stated below refer to winter conditions, as a worst possible scenario for steam network operation. At the period corresponding to the data taken for simulation in winter conditions, the average HPS temperature at C6 unit was reaching 301.9 • C. Table 10 contains results of all simulations. As expected, these projects would affect mostly units adjacent to the place of realization of proposed projects (C6, C5, C4, C7, C8 and C9). Since the change of steam temperature in other units was negligible, they are not listed in Table 10. The results show, that in case of realization of project A only, the steam temperature at unit C6 will drop by almost 15 • C (assuming unsuitable technical state of thermal insulation). Since the temperature is already low, this decrease would take the steam even closer to its condensation point. Following the results shown in Table 10, one can conclude: • Only the connection of units C6 and S2 itself will not improve the steam temperature at C6 unit. Continuous operation of this pipeline will not be effective, otherwise a significant temperature drop will be observed leading to unacceptable operation conditions. • Also commissioning the new unit C10 will not improve the steam quality at C6 unit (case A+B). Although it will lower the load of unit S2, the supplied steam from S2 will flow in the direction of unit C6 anyway.

•
The operation of new unit S5 would improve the steam temperature at unit C6 only slightly (case A+C). The reason is the high steam pressure of steam supplied from unit S2 that will make the steam from unit S5 to flow in the opposite direction (from unit S5 to S1 and C5).
Significant steam temperature increase at units C7 to C9 in case A and A+B can be explained by decreased steam supply from unit S2 to pipeline F. This will be compensated by increased steam supply from the CHP unit of high quality with temperature over 360 • C. The calculated pressure drop on the section between units S2 and C6 is 5 kPa, which is comparable to pressure drops in other sections (over the same distance). One can conclude, that, supposing constant steam pressure at other units, the steam flow in this section will be in direction from unit S2 to unit C6 and not in the opposite one. According to these statements, we recommend stopping project A.
Project C will cause a slight temperature decrease at units C7 to C9. Although new unit S5 would supply large amount of steam of relatively high temperature, the temperature will decrease due to heat loss through transport over a long distance from units S5 to units C7, C8 and C9.
Another remarkable observation is that, in summer conditions, the steam supply from the CHP unit would be too low and some action would have to be made. At present, there is a return pipeline "C" that is used to allow the steam flow back to the CHP unit to be consumed there. Another, but less economic, way is to use the steam let-downs or steam drives.
Comparison of presented calculation results with real system operation data cannot be made as none of the examined projects is in proper operation yet. Project C has been postponed by at least three years; project B, though commissioned, did not reach a stable long-term operation yet and project A is on hold.

Improvement of Steam Network Operation
There is a significant variability of steam consumption during the year and the steam network is designed to handle different amounts of steam under any conditions. The balance model has proven significant imbalance within the steam network. Certain part of this amount is lost as a part of condensate being disposed into the ground. Few spots with leaking steam into the atmosphere were observed during visual inspections as well; however, these losses are not as high as the result corresponding to the balance model.
Despite steam and condensate losses, the most important factor is the inaccuracy of measurement. A simple algorithm, presented in Figure 10 was developed and used to recommend following general measures, concerning routine and preventive maintenance of steam pipelines, along with the best practice in steam system management [49]: • Regular inspection of steam traps' operation; • Regular visual inspection of steam pipelines and isolation; • Improve the steam condensate recovery system and prevent from its disposal into the ground; • Regular inspection and calibration of measurement tools. There is a significant variability of steam consumption during the year and the steam network is designed to handle different amounts of steam under any conditions. The balance model has proven significant imbalance within the steam network. Certain part of this amount is lost as a part of condensate being disposed into the ground. Few spots with leaking steam into the atmosphere were observed during visual inspections as well; however, these losses are not as high as the result corresponding to the balance model.
Despite steam and condensate losses, the most important factor is the inaccuracy of measurement. A simple algorithm, presented in Figure 10 was developed and used to recommend following general measures, concerning routine and preventive maintenance of steam pipelines, along with the best practice in steam system management [49]: • Regular inspection of steam traps' operation; • Regular visual inspection of steam pipelines and isolation; • Improve the steam condensate recovery system and prevent from its disposal into the ground; • Regular inspection and calibration of measurement tools. The inspection of steam traps is done on a regular basis over a period of 2 years. In the last 2 years, more than 13,000 steam traps have been inspected. More than 1000 were out of operation at that time. Of the remainder, 91.2% were working correctly, 7.7% were leaking and only 1.1% were clogged. Steam traps that are not working correctly are replaced by new ones immediately after inspection.

Conclusions
The work deals with mathematical modeling as a tool for steam network optimization. The authors developed a simple model based on nodes and pipes and, therefore, it is applicable to any steam network. Simultaneous calculation of matrix Equations allows shortening the calculation time. Model was tested on steam network in SLOVNAFT oil refinery with real process data from the year 2017. Required data consisted mainly of steam flow, temperature and pressure measurements on individual production units, which were used as input data for both models used. The most challenging issue was the lack of certain data, production unit shutdowns and other measurement or data transfer-related issues.
The balance model performs the calculation of steam flow velocity in each pipe, in addition to mass and heat balance. This allowed identification of troublesome spots within the steam network in terms of low (cold spots) or high steam flow velocity (bottlenecks). Results were partly used as input data for the thermo-hydrodynamic model. This model uses exact equations and correlations for pressure drop, heat loss or steam properties calculation. A probabilistic model employing a similar The inspection of steam traps is done on a regular basis over a period of 2 years. In the last 2 years, more than 13,000 steam traps have been inspected. More than 1000 were out of operation at that time. Of the remainder, 91.2% were working correctly, 7.7% were leaking and only 1.1% were clogged. Steam traps that are not working correctly are replaced by new ones immediately after inspection.

Conclusions
The work deals with mathematical modeling as a tool for steam network optimization. The authors developed a simple model based on nodes and pipes and, therefore, it is applicable to any steam network. Simultaneous calculation of matrix Equations allows shortening the calculation time. Model was tested on steam network in SLOVNAFT oil refinery with real process data from the year 2017. Required data consisted mainly of steam flow, temperature and pressure measurements on individual production units, which were used as input data for both models used. The most challenging issue was the lack of certain data, production unit shutdowns and other measurement or data transfer-related issues.
The balance model performs the calculation of steam flow velocity in each pipe, in addition to mass and heat balance. This allowed identification of troublesome spots within the steam network in terms of low (cold spots) or high steam flow velocity (bottlenecks). Results were partly used as input data for the thermo-hydrodynamic model. This model uses exact equations and correlations for pressure drop, heat loss or steam properties calculation. A probabilistic model employing a similar approach as the balance model was recently used to assess and predict the transport capacities, bottlenecks and system operation change after inclusion of new infrastructure in the real gas transmission network in Europe [50].
Comparing the results of simulations to the real measured data, one can conclude that the real state of thermal insulation is bad, since the thermo-hydrodynamic model initially estimated higher temperatures than the measured ones. Subsequently, the simulations were performed again with a higher value of thermal conductivity corresponding to obsolete and damaged thermal insulation of pipes. The results obtained were closer to reality.
A thermo-hydrodynamic model was used to assess the operation of steam networks after implementation of the proposed changes. Simulations of multiple variations were performed with results leading to one conclusion. The connection of units S2 and C6 by a common pipeline does not lead to the desired result. Therefore, it is recommended to abort this action if possible. This model can be used in the future for prediction of future operation states as well. A comparison of model results and real operational data cannot be performed, however, as none of the analysed projects impacting the high-pressure steam network management is in proper operation yet.
The proposed thermo-hydrodynamic model is very complex and sufficiently accurate but requires corresponding skills and knowledge among its users. Moreover, its results are only as reliable as the measurements themselves. During the cooperation with the refinery, a few controlled operation experiments were performed to determine the steam flow in pipe interconnections and to assess the accuracy of individual measurement tools as well. It is recommended to continue with these experiments in the future as well and improve the overall state of steam networks.

Acknowledgments:
The authors would like to express many thanks to all SLOVNAFT, a.s., employees who contributed to the final scope and form of this paper.

Conflicts of Interest:
The authors declare no conflict of interest.