Regular Article
Influence of accounting for thermodynamic processes on the processes of mixture formation during sequential pumping of petroleum products
National University of Oil and Gas « Gubkin University », 65 Leninsky Prospekt, Moscow 119991, Russia
^{*} Corresponding authors: imvanchugov@yandex.ru, shestakov.r@gubkin.ru
Received:
9
November
2020
Accepted:
15
December
2020
In this paper, we study the processes of mixing when using the technology of transportation of light oil productssequential pumping by direct contact. Modeling of mixing processes is carried out taking into account the influence of parametric and thermodynamic factors. For numerical modeling, a software package was developed that allows not only modeling and calculating the parameters of the oil pipeline operation in real time with subsequent graphical visualization, but also comparing them with real data processed by operators at production facilities.
© I.M. Vanchugov & R.A. Shestakov, published by IFP Energies nouvelles, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Currently, one of the most costeffective and efficient types of transportation of petroleum products is pipeline transport [1–3]. The main oil pipelines that are already in operation are being transferred to the transportation of oil products [4]. However, despite its many advantages, it also has disadvantages, especially with regard to the sequential pumping of various grades of petroleum products through one pipeline, namely, the formation of a mixture of neighboring grades of petroleum products during sequential pumping [5–13]. In modern practice, the calculation methods and algorithms used in the design and operation of pipelines through which oil and petroleum products are pumped are based on the works of major pipeline transport scientists V.G. Shukhov, L.S. Leibenzon, M.V. Lurie, and others, as well as on works published in relatively recent times [5, 14–20]. Research is also being carried out to reduce the volume of a mixture of differentgrade oil products that occurs during sequential pumping [21], expand the range of oil products pumped by direct contact [6, 22], and improve the accuracy of monitoring the boundaries of the mixing zone using instrument methods [23]. Also, the issue of managing, determining the volume of shipments and their distribution during pipeline transportation has been widely studied in the framework of multicriteria management and logistics tasks [24–30].
This fact demonstrates that the subject is relevant, but even more important is research into the influence of nonisothermal transfer in the processes of mixture formation and the parameters of the pumping mode for pipelines, working on a method of direct contact, but the pumping oil temperature effect is studied very well [5, 14–17, 31–33].
Thus, to assess the influence of temperature accounting, it is necessary to build a mathematical model of pipeline transportation of petroleum products by direct contact on the section of the main oil product pipeline between pumping stations. To do this, it is necessary to take into account not only the parameters of oil products, the pipeline and the pumping mode [34], but also the temperature distribution over the oil product pipeline section and the initial temperature of the pumped products. This is due to the fact that the hydraulic calculation is influenced by the parameters and results of thermal calculation, so this work eliminates simplifications and errors that are not taken into account in the isothermal calculation of the operation of the oil product pipeline and the processes of mixing as a result of its operation.
2 Mathematical model and algorithm
The primary task in modeling nonisothermal pumping is to simulate the movement of each oil product throughout the entire section of the oil product pipeline separately (assuming that this oil product completely fills the pipeline), to conduct intermediate calculations and obtain data that will be used in hydraulic calculation and in calculating the volume of the mixture.
The main essence of the mathematical model of nonisothermal sequential pumping of oil and petroleum products by direct contact method is that this model consists of several simpler models, as well as refinement by taking into account additional pumping parameters. The connection of petroleum products into a single mathematical model is carried out when modeling the operation of the entire pipeline.
Let’s assume that oil products are pumped through a section of the main oil product pipeline, and n ≥ 2, and the entire section of the oil product pipeline is filled with transported oil products at the initial time.
The simulation is performed in several iterations. First of all, the flow of each oil product is modeled separately, assuming that the entire oil product pipeline is filled with jm oil products, j = 1 … n, during which the thermodynamic parameters of pumping are determined – the heat transfer coefficient from the flow to the external environment, the primary temperature distribution along the axis of the oil product pipeline, as well as the parameters of the jth oil product. During the next iteration, a more accurate distribution of the temperature of the jth oil product along the axis of the oil product pipeline will be determined using the refined thermodynamic parameters of pumping.
To determine the state of the entire system, it is necessary to place batches of products in the oil product pipeline, assuming that at the initial time the first pumped product has reached the end of the linear section, and taking into account the mass of each oil product, as well as its density, obtained by refined thermodynamic calculation. The contact boundaries of petroleum products will be the middle of the mixture formed by them during sequential pumping.
Based on the main parameters obtained during the refined thermodynamic calculation, the volume of the mixture, its boundaries with each of the oil products, as well as the distribution of concentrations, density and viscosity of the mixture are determined. Thus, simple mathematical models of nonisothermal pumping of each oil product are combined into a single complex model, but this is not all.
To achieve a better connection between the model and industry, the final model was integrated with a model of operation of main and backup pumps with the possibility of implementing various connection schemes and operating modes.
In order to perform hydraulic and thermodynamic calculations based on the presented model, it is necessary to know the parameters of the pipeline, the transported products, and the temperature distribution of the environment surrounding the oil product pipeline.

Parameters of the linear part of the pipeline:
x – coordinate, m;
z(x) – height, m;
D – outer diameter, mm;
δ – wall thickness, mm;
d = (D – 2 · δ)/1000 – outer diameter, m;
Δ – the absolute roughness of the wall, mm;
T_{surr. area} (x) – the temperature of the surrounding area (soil), K;
λ_{wl.pl.} (x) – coefficient of thermal conductivity of the pipeline wall material, W/(m K);
λ_{surr. area} (x) – soil thermal conductivity coefficient, W/(m K);
h_{0} (x) – the depth of the pipeline, m;
T_{wl. pl.} (x) – the wall temperature of the pipeline, K.

Parameters of pumped oil products:
ρ_{0} – density of petroleum product under normal conditions, kg/m^{3} ;
ν_{0} – kinematic viscosity of petroleum products under normal conditions, cSt;
T_{n.c.} – temperature of normal conditions, K;
m – weight of oil product, t;
– initial specific heat of the oil product, W/(m K);
– coefficient of thermal conductivity of petroleum product under normal conditions, W/(m K);
α_{0} – average temperature correction of density when the oil product temperature changes by 1 K;
u_{0} – coefficient of steepness of the viscogram;
T_{0} – initial temperature of the oil product, K.
The calculation algorithm is described below, during which the main technological parameters are determined using the proposed mathematical model, taking into account the above clarifications.
First of all, the external heat transfer coefficient from the pipe wall to the environment is determined by the formula (1):
After that, the distribution of parameters of each transported oil product along the pipeline length is calculated. The Isobaric heat capacity is determined in accordance with (2), the density is (3), and the kinematic viscosity coefficient is (4):
where
j – oil product number in the batch,
T_{wl.pl.} (x)_{j} – temperature of the pipeline wall, in crosssection with x coordinate, if x = 0, then .
Since the physical parameters of the transported products are not constant and change along the length of the pipeline, and the ambient temperature can change along the length of the route, the parameters of the process of TV transmission from the transported products to the environment will be represented as functions for each transported product in the form of a distribution along the length of the pipeline. The main thermal parameter of the pipeline is the transfer coefficient, which is determined according to (9), but its distribution along the length of the pipeline for each product to calculate the distribution coefficient of heat conductivity of oil according to (5), Prandtl number (6), as the velocity of the oil – (7), the internal heat transfer coefficient from the fluid to the wall of the pipe (8):
where ϵ_{k} (x) is correction factor that depends on the ratio of the length of the pipe section to its internal diameter in turbulent conditions ϵ_{k} (x) = 1:
Next, the temperature distribution of the pumped product along the length of the pipeline is determined, taking into account heat exchange with the environment and energy dissipation to overcome friction forces with the release of heat. To do this, it is necessary to perform a calculation algorithm according to formulas (10)–(15) with a step Δx. This algorithm is described in [5, 14], but it is used to determine the average pipeline parameters of oil, which is not suitable for carrying out refined calculations of sequential pumping of petroleum products.
Shukhov coefficient of petroleum products:
Reynolds number of petroleum product:
The coefficient of relative roughness of pipes:
The coefficient of hydraulic resistance, in a section with a coordinate x according to [5, 15]. Hydraulic slope ratio:
where g = 9.80665 м/c^{2} – acceleration of free fall.
Coefficient b for determining the influence of viscous friction forces in the cross section with the coordinate:
The temperature distribution of an oil product is a set of temperature values:
Next step of calculation:
Calculations are performed for each of the points with coordinates x ∈ [0; L] and step Δx = L/1000, where L – is length of the oil product pipeline section.
After the determination of temperature distribution of products along the length of the plot is the updated recalculation of the distribution of values of the physical properties of products (2)–(4). That allows to increase the accuracy of determining the volume of each batch of petroleum product during its movement through the pipeline, the accuracy of calculating the parameters of the transportation mode of petroleum products and the volume of the mixture of the neighboring oil, which is determined by (17) [5, 17]:
where is internal volume of the part of the pipeline passed by the party.
So, after constructing a mathematical model, it is necessary to numerically evaluate the effect of changes in the temperature of petroleum products along the length of the pipeline on the volume of the resulting mixture and the parameters of the pumping mode.
3 Results and discussion
As an example, let’s consider three variants of a pipeline with a variety ∅630 × 7 mm, ∅530 × 7 mm and ∅426 × 7 mm, the other parameters are assumed to be the same for all three varieties: ∆ = 0.01 mm, λ_{wl.pl.} = 44 W/(m K), λ_{soil} = 1.05 W/(m K), T_{wl.pl.} (summer) = 288.15 K, T_{wl.pl.} (winter) = 263.15 K. The route profile is shown in Table 1, the ambient temperature distribution along the pipeline route is shown in Table 2, the parameters of transported products are shown in Table 3.
Profile of the linear part of the pipeline.
Ambient temperature along the pipeline route at different times of the year.
For numerical simulation, a software package was developed in accordance with the mathematical algorithm described above, which shows the numerical values of the pumping parameters in a crosssection (Figs. 2 and 3). The input data window is shown in Figure 1.
Fig. 1 Form for entering characteristics of pumped liquids and parameters of the linear part of the pipeline. 
Fig. 2 Form for output of calculation of oil product parameters along the pipeline axis. 
Fig. 3 Form for output of calculation of oil product parameters along the pipeline axis. 
As the results of numerical modeling of sequential pumping of petroleum products, the volumes of a mixture of contacting petroleum products – automobile gasoline brands were determined AI76, AI80, AI92 and AI95, as well as the pressure and head drop required to ensure the specified flow rate.
To conduct a numerical study, a simulation was carried out taking into account different pumping conditions by time of yearsummer and winter, which will allow us to clarify the volume of mixtures obtained during pumping, as well as the pressure and head losses for pumping in a given mode. As a result of various factors, it may be necessary to perform reverse pumping, so to study the direction of pumping within the framework of changing climatic parameters, numerical modeling of sequential pumping for selected pipeline grades from South to North and from North to South during the summer and winter periods of operation was carried out. However, in order to compare and determine the degree of influence of the inhomogeneity of the temperature distribution on the pipeline route, isothermal calculations were performed during the summer and winter periods of operation.
The numerical simulation results to pipeline export product mix ∅426 × 7 mm with a volume flow of 1044 m^{3}/h are presented in Table 4, a graphical rsepresentation of the dependence of the volume of the mixture from pump – in Figure 4, and values of the required differential pressure and head – Figure 5.
Fig. 4 Volume of the mixture formed under different pumping conditions. 
Fig. 5 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 
The results of the studies under Q = 1044 m^{3}/h.
Numerical results of simulation of pipeline pumping by assortment ∅530 × 7 mm with the expense 1800 m^{3}/h, is shown in Table 5, a graphical representation of the dependence of the volume of the mixture on the pumping conditions is shown in Figure 6, and the values of the required pressure drop and head are shown in Figure 7.
Fig. 6 Volume of the mixture formed under different pumping conditions. 
Fig. 7 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions 
Research result where Q = 1800 m^{3}/h.
Numerical results of modeling pumping through a pipeline with a range of ∅630 × 7 mm with a flow rate of 2880 m^{3}/h are shown in Table 6, a graphical representation of the dependence of the mixture volume on pumping conditions is shown in Figure 8, and the values of the required pressure drop and head are shown in Figure 9.
Fig. 8 Volume of the mixture formed under different pumping conditions. 
Fig. 9 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 
Research result where Q = 2880 m^{3}/h.
Based on the results of calculating the sequential pumping of four brands of automobile gasoline, it is worth noting the general and generally obvious trend for all grades that energy losses will increase with a decrease in temperature, which will lead to the need for more precise adjustment of pumping equipment at oil pumping stations.
As a result of the updated nonisothermal calculation for summer conditions, it should be noted that the required head for pumping in the specified modes is less than for the isothermal calculation for all grades, which may lead to an incorrect decision on regulating the pumping equipment and reducing the efficiency of the entire system. This is primarily due to the friction forces and temperature distribution along the length of the oil product pipeline. In winter, the opposite trend is observed, despite the forces of viscous friction, the head loss for pumping has become greater than for nonisothermal. This fact is due to the fact that for isothermal pumping, the average temperatures of the pumped product were taken in each of the periods. Failure to take into account thermalhydraulic factors will lead to incorrect adjustment of the pumping equipment and the inability to provide the specified pipeline capacity, but increasing the accuracy of the model and, as a result, calculations will allow you to determine more accurately the energy loss at a given pumping mode and make the required adjustment.
One of the main results of modeling is the dependence of the volume of the resulting mixture on the viscosity of the contacting products, which is primarily due to the uneven displacement of the first product by the next product, and secondarily to the turbulent diffusion in the contact zone of two neighboring products.
As is known, the viscosity of petroleum products strongly depends on the temperature of the oil product under consideration, so it is important to take into account its effect on the mixture formation as accurately as possible. In the course of numerical modeling of mixture formation, the volume of the mixture was the smallest in summer when pumping from South to North, and the largest when pumping from North to South in winter. The calculated volume of the mixture during isothermal pumping for summer conditions turned out to be more than the specified one, and for winter time – less.
So, usually at the final point of the main oil product pipeline, the layout of the mixture of petroleum products is carried out by net, i.e. adding the mixture to the oil products that participated in the [5, 6, 10]. If the mixture is formed more than can be decomposed into individual oil products, then it must be separated again, but at the refinery [27, 29]. For this reason, the inaccuracy of modeling the mixture formation in winter will lead to the need to send part of the mixture for recycling, which will lead to additional problems and reduce the efficiency of the entire enterprise.
In the summer, the specified volume of the mixture has become smaller, and the head loss for pumping has also decreased, which will allow for more flexible regulation of pumping through the pipeline and, if necessary, slightly reduce the minimum possible volume of the batch of oil products and increase the pumping speed.
The simulation of petroleum products similar in their physical properties – automobile gasoline was performed above, but it is necessary to assess the impact of new factors taken into account when pumping several gasoline and Diesel Fuel (DF), for example, AI95, AI92, AI80 and Diesel fuel. With sequential pumping of gasoline and diesel fuel, diesel fuel is located in a cycle between batches of the worstquality gasoline [5, 6, 17].
We will conduct a simulation of sequential pumping along the section of the oil product pipeline presented earlier, with the same parameters and flow rate, but we will replace AI76 with DF, and in summer we will consider Summer Diesel Fuel (SDF), and in winter – Winter Diesel Fuel (WDF). Their main differences that are of interest to us are the different values of kinematic viscosity and density in SDF and WDF [36] (Tab. 3), as well as their application in different temperature conditions.
Numerical results of modeling pumping through a pipeline with a range of ∅426 × 7 mm with a flow rate of 1044 m^{3}/h are shown in Table 7, a graphical representation of the dependence of the mixture volume on pumping conditions is shown in Figure 10, and the values of the required pressure drop and head are shown in Figure 11.
Fig. 10 Volume of the mixture formed under different pumping conditions. 
Fig. 11 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 
Research result where and pumping of gasoline and diesel fuel.
From the results obtained, it can be seen that during sequential pumping of various petroleum products by direct contact, the properties of which are quite different from each other, for example, automobile gasoline and diesel fuel, the volume of the mixture between gasoline (AI92 + AI80) for summer operating conditions was 68.73 m^{3} and 69.53 m^{3}, depending on the direction of pumping, i.e. did not change. However, the volume of the mixture formed in the contact zone of automobile gasoline and diesel fuel (AI80 + Diesel fuel) was significantly larger and amounted to 88.11 m^{3} and 91.34 m^{3} for summer operating conditions, depending on the pumping direction. This is due to the fact that diesel fuel is a more viscous petroleum product compared to automobile gasoline.
Consider the worst pumping conditions – in winter and the direction of pumping from North to South. Under these conditions, the volume of the mixture will increase in the contact zone of all petroleum products and for AI92+AI80 will be 84.27 m^{3}, and the increase will be 22.6% of the most favorable conditions. This fact is associated with an increase in viscosity due to a decrease in the temperature of the pumped product and the environment. For AI80+winter diesel fuel, the volume of the mixture will be 146.83 m^{3}, which corresponds to an increase of 74.24% when pumping in the summer period of operation in the direction from South to North. This is due to the fact that the viscosity of diesel fuel has a more pronounced dependence on temperature.
In addition, due to the participation in the sequential pumping of more viscous and dense oil products, the losses for pumping have sharply increased. this fact is most pronounced for winter pumping conditions and leads to the fact that the existing pumping equipment suitable for pumping in the summer time simply will not be able to pump over a given section. There are two ways out of this situation:
Necessary to install additional pumping equipment that will only be used in winter, which entails additional capital investments, but the increase in pressure at the outlet of the pumping station may not be able to withstand the pipeline itself, so it will be necessary to lay a safety margin by increasing the wall thickness or using other steel grades that are stronger and more expensive, which will also increase the overall cost of the project. This option may be relevant for the planned or projected oil product pipeline.
Necessary to reduce the volume flow rate on the section of the oil product pipeline. This option is relevant for already commissioned and operated oil product pipelines. In this case, there is no need to increase the outlet head and design the pipeline with a safety margin, but there are also disadvantages.
Consider the option of reducing the volume flow, the results of the modeling pipeline, range ∅426 × 7 mm presented in Table 8, a graphical representation of the dependence of the volume of the mixture from pump – in Figure 12, and the values of the required differential pressure and head in Figure 13. The reduced flow rate in the considered variant was 900 m^{3}/h, which led to a reduction in losses from 695.62 m to 532.75 m.
Fig. 12 Volume of the mixture formed under different pumping conditions. 
Fig. 13 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 
Research result where and pumping of gasoline and diesel fuel.
According to the results obtained, when the volume flow rate of pumping decreases, the volume of the mixture of transported products will increase. So for the contact zone AI92+AI80 under winter conditions, pumping will be 88.77 m^{3} (against 84.27 m^{3} before the reduction of pumping), which corresponds to an increase of 5.34%. For AI80+winter diesel fuel – 156.2 m^{3} (against 146.83 m^{3} before the reduction of pumping), which corresponds to 6.38%.
The formation of the mixture as a whole is a negative consequence of the sequential pumping of various petroleum products, not to mention its increase, which will require recalculation of the minimum possible batch volume for its layout at the final point, provided that there is a quality reserve for the transported fuels.
In the following, numerical results of modeling the operation of a section of the main oil product pipeline for sequential pumping for grades ∅530 × 7 mm and ∅630 × 7 mm at the initial flow rate are presented in Tables 9 and 11, respectively, and a graphical representation in Figures 14, 15 and 18, 19, respectively. Also these modeling options are characterized by an increase in losses when the AI76 is replaced by diesel fuel, in summer – by summer, in winter – by winter.
Fig. 14 Volume of the mixture formed under different pumping conditions. 
Fig. 15 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 
Fig. 16 Volume of the mixture formed under different pumping conditions. 
Fig. 17 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 
Fig. 18 Volume of the mixture formed under different pumping conditions. 
Fig. 19 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 
Research result where and pumping of gasoline and diesel fuel.
Research result where and pumping of gasoline and diesel fuel.
Research result where and pumping of gasoline and diesel fuel.
Therefore, for assortments ∅530 × 7 mm and ∅630 × 7 mm also been carried out by lowering the volumetric flow of 1800 m^{3}/h to 1620 m^{3}/h, and 2880 m^{3}/h to 2520 m^{3}/h, the results of which are presented not only for winter but for the summer mode of operation in Tables 10 and 12, respectively, and in Figures 16, 17 and 20, 21 respectively. As a result, for the other two grades, a decrease in the volume flow rate of pumping also led to an increase in the volume of the mixture formed during pumping. Changes in the volume of the mixture at a reduced mode are shown in Table 13.
Fig. 20 Volume of the mixture formed under different pumping conditions. 
Fig. 21 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 
Research result where and pumping of gasoline and diesel fuel.
Change in the volume of the mixture when the pumping mode is lowered.
4 Conclusion
Despite the fact that sequential pumping of petroleum products of various brands by direct contact is one of the most costeffective modes of transport for large volumes of transportation, this method has a number of disadvantages, so taking into account any additional factors is important, which showed:
the volume of the mixture when pumping in the NorthSouth direction is greater than in the SouthNorth direction;
when pumping in winter, the volume of the mixture increases;
in the implementation of pumping in the winter there is an increase of head losses and pressure;
when pumping heavier fuels, the properties of which are more dependent on temperature, there are large fluctuations throughout the year in the values of the mixture volumes, pressure losses and head losses on the section of the oil product pipeline.
Taking into account the change in the temperature distribution of the soil surrounding the pipeline along its route makes it possible to more accurately calculate the volume of the resulting mixture of oil products. This allows you to accurately estimate the required volume of the batch required for the implementation for the complete layout of the mixture for the transported oil products in the presence of their quality stock [10].
Based on the strong fluctuations during the year in the values of head losses and pressure losses, it was proved that it is necessary to regulate pumping equipment or pumping flow rate when transporting heavier fuels with sequential pumping for the implementation of the pumping process itself.
The proposed algorithm makes it possible to establish the possibility of tracking the process of transportation of oil products during sequential and nonisothermal pumping by direct contact in real time, taking into account parametric factors.
References
 Deineko S.V., Alikhashkin A.S., Shestakov R.A., Ulanov V.V. (2018) Main technological equipment and processes for oil and petroleum products transportation, National University of Oil and Gas «Gubkin University», Moscow, Russia. [Google Scholar]
 Kitaev S.V. (2018) Technical diagnostics of oil and oil products pipeline transport facilities, Ufa State Petroleum Technological University, Ufa, Russia. [Google Scholar]
 Vasilkovsky V.V. (1991) Development of energysaving technology for sequential pumping through branched oil product pipelines, PhD Thesis, National University of Oil and Gas «Gubkin University», Moscow. [Google Scholar]
 Sunagatullin R, Timofeev F, Kuznetsov A, Oludina Y (2019) Relevant issues on quality evaluation of petroleum pipeline preparation for oil product transportation, Oil & Gas Science and Technology – Rev. IFP Energies nouvelles 74, 44. https://doi.org/10.2516/ogst/2018098. [EDP Sciences] [Google Scholar]
 Lurie M.V., Mastobayev B.N., RevelMuroz P.A., Sosenko A.E. (2019) Design and operation of oil pipelines, National University of Oil and Gas «Gubkin University», Moscow, Russia. [Google Scholar]
 Korolenok A.M., Lurie M.V., Timofeev F.V. (2012) Expansion of the range of light oil products transported through pipelines using the sequential method, Science and Technology of Pipeline Transportation of Oil and Petroleum Products 8, 4, 40–43. [Google Scholar]
 Zorya E.I., Korolenok A.M., Losenkova O.V., Yu N.K. (2018) Fundamentals of resource conservation in the turnover of hydrocarbons, Max Press, Moscow, Russia. [Google Scholar]
 Magomedov I.R. (2019) Problems of terminology in customs control of oil and petroleum products, Economics 180, 11, 65–70. [Google Scholar]
 Stigneev I.V. (2017) Analysis of Russian oil and oil products exports under sanctions, Customs and foreign economic activity of companies 3, 2, 1–8. [Google Scholar]
 Lurie M.V., Timofeev F.V., Sereda S.V. (2017) Mix layout for sequential pumping of petroleum products, Science and Technology of Pipeline Transportation of Oil and Petroleum Products 7, 2, 42–47. [Google Scholar]
 Yablonsky V.S., Yufin V.A., Budarov I.P. (1959) Sequential pumping of petroleum products and oils through main pipelines, Gostoptekhizdat, Moscow, Russian. [Google Scholar]
 Lurie M.V. (1979) Optimization of sequential pumping of petroleum products, Nedra, Moscow, Russia. [Google Scholar]
 Nechval M.V., Novoselov V.F., Tugunov P.I. (1976) Sequential pumping of oil and oil products through main pipelines, Nedra, Moscow, Russia. [Google Scholar]
 Polyakov V.A., Shestakov R.A. (2018) Systematic approach to pipeline design, National University of Oil and Gas «Gubkin University», Moscow, Russia. [Google Scholar]
 Didkovskaya A.S., Lurie M.V. (2016) Iterative algorithm for hydraulic calculation of nonisothermal oil pumping, Science and Technology of Pipeline Transportation of Oil and Petroleum Product 2, 50–55. [Google Scholar]
 Lurie M.V. (2014) Algorithm for calculating quasisteadystate modes of nonisothermal oil pumping, Science and Technology of Pipeline Transportation of Oil and Petroleum Product 2, 128–131. [Google Scholar]
 Lurie M.V. (2012) Mathematical modeling of pipeline transportation of oil, oil products and gas, National University of Oil and Gas «Gubkin University», Moscow, Russia. [Google Scholar]
 Gong J., Wang Q., Wang W., Guo Y. (2010) The calculation method of mixing volume in a products pipeline, in: 8th International Pipeline Conference, Calgary, Canada, Vol. 3, pp. 393–398. [Google Scholar]
 Dombalov D.U., Melnikov D.I., Drozdov D.A. (2016) Determination of the effective diffusion coefficient and volume of the mixture during sequential pumping of fuel of various groups and grades by direct contact on collapsible pipelines, Science and Technology of Pipeline Transportation of Oil and Petroleum Product 6, 1, 30–38. [Google Scholar]
 Shestakov R.A. (2020) Research of distribution of oil flow in the pipeline with looping, J. Phys. Conf. Series 1679, 052035. https://doi.org/10.1088/17426596/1679/5/052035. [Google Scholar]
 Golunov N.N. (2018) Parameters of sequential pumping of petroleum products using small antiturbulence additives to reduce the volume of the resulting mixture, Territory of Neftegaz 13, 5, 68–72. [Google Scholar]
 Mayorov N.S., Istomova M.A. (2020) Sequential pipeline pumping of petroleum products, Business Magazine Neftegaz. RU 99.5, 3.5, 45–50. [Google Scholar]
 Sereda S.V., Lyapin A.Y., Dubovoy E.S., Timofeev F.V., Vishnevskaya Yu.A. (2020) Analysis of the possibility of using the method of infrared spectroscopy to determine the boundaries of the zone of mixing of petroleum products, Science and Technology of Pipeline Transportation of Oil and Petroleum Product 10, 1, 64–69. [Google Scholar]
 Li Y., Chen M., Qiang X. (2020) Simultaneous scheduling of multiproduct pipeline distribution and depot inventory management for petroleum refineries, Chem. Eng. Sci. 220, 115618. [Google Scholar]
 Dimas D., Murata V.V., Neiro S.M.S., Relvas S., BarbosaPóvoa A.P. (2018) Multiproduct pipeline scheduling integrating for inbound and outbound inventory management, Comput. Chem. Eng. 115, 377–396. [Google Scholar]
 Zhang H., Liang Y., Liao Q., Shen Y., Yan X. (2018) A selflearning approach for optimal detailed scheduling of multiproduct pipeline, J. Comput. Appl. Math. 327, 41–63. [Google Scholar]
 Kirschstein T. (2018) Planning of multiproduct pipelines by economic lot scheduling models, Eur. J. Oper. Res. 264, 327–339. [Google Scholar]
 Zhou X., Zhang H., Xin S., Yan Y., Long Y., Yuan M., Liang Y. (2020) Future scenario of China’s downstream oil supply chain: Low carbonoriented optimization for the design of planned multiproduct pipelines, J. Clean. Prod. 244, 118866. [Google Scholar]
 Chen H., Zuo L., Changchun W., Li Q. (2019) An MILP formulation for optimizing detailed schedules of a multiproduct pipeline network, Transp. Res. Part E 123, 142–164. [Google Scholar]
 Mostafaei H., Castro P.M., GhaﬀariHadigheh A. (2015) A novel monolithic MILP framework for lotsizing and scheduling of multiproduct treelike pipeline networks, Indus. Eng. Chem. Res. 54, 37, 9202–9221. [Google Scholar]
 Zhapbasbayev U.K., Makhmotov E.S., Bekibayev T.T., Ramazanova G.I., Rziyev S.A. (2015) Calculation of the optimal pumping temperature for oil transportation, Science and Technology of Pipeline Transportation of Oil and Petroleum Product 5, 4, 61–65. [Google Scholar]
 Didkovskaya A.S., Lurie M.V. (2015) Universal algorithm for numerical calculations of stationary modes of oil pipelines operation, Science and Technology of Pipeline Transportation of Oil and Petroleum Product 5, 4, 86–91. [Google Scholar]
 Zhapbasbaev U.K., Makhmotov E.S., Ramazanova G.I., Bekibaev T.T. (2018) Calculations of energysaving modes of sequential pumping of oil mixtures on the section of the main oil pipeline, Science and Technology of Pipeline Transportation of Oil and Petroleum Product 8, 3, 30–38. [Google Scholar]
 Vanchugov I.M. (2019) Computer simulation of mixing processes during sequential pumping of oil and petroleum products, in: 73rd International Youth Scientific Conference «Oil and Gas2019», National University of Oil and Gas «Gubkin University», Moscow, Russia, pp. 106–107. [Google Scholar]
 GOST 32513–2013 Motor Fuels. Unleaded petrol. Technical conditions. [Google Scholar]
 GOST 305–2013 Diesel Fuel. Technical conditions. [Google Scholar]
All Tables
All Figures
Fig. 1 Form for entering characteristics of pumped liquids and parameters of the linear part of the pipeline. 

In the text 
Fig. 2 Form for output of calculation of oil product parameters along the pipeline axis. 

In the text 
Fig. 3 Form for output of calculation of oil product parameters along the pipeline axis. 

In the text 
Fig. 4 Volume of the mixture formed under different pumping conditions. 

In the text 
Fig. 5 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 

In the text 
Fig. 6 Volume of the mixture formed under different pumping conditions. 

In the text 
Fig. 7 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions 

In the text 
Fig. 8 Volume of the mixture formed under different pumping conditions. 

In the text 
Fig. 9 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 

In the text 
Fig. 10 Volume of the mixture formed under different pumping conditions. 

In the text 
Fig. 11 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 

In the text 
Fig. 12 Volume of the mixture formed under different pumping conditions. 

In the text 
Fig. 13 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 

In the text 
Fig. 14 Volume of the mixture formed under different pumping conditions. 

In the text 
Fig. 15 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 

In the text 
Fig. 16 Volume of the mixture formed under different pumping conditions. 

In the text 
Fig. 17 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 

In the text 
Fig. 18 Volume of the mixture formed under different pumping conditions. 

In the text 
Fig. 19 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 

In the text 
Fig. 20 Volume of the mixture formed under different pumping conditions. 

In the text 
Fig. 21 Changing the pressure and head parameters at the ends of the linear section under different pumping conditions. 

In the text 