Journal Information
Title: Enfoque UTE
Copyright: 2024, Enfoque UTE
Abbreviated Title: Enfoque UTE
Volume: 15
Issue: 3
ISSN (electronic): 13906542
Copyright statement: License (openaccess,
https://creativecommons.org/licenses/by/3.0/ec/):
Article Information
Date received: 27 febrero 2024
Date revised: 05 abril 2024
Date accepted: 20 mayo 2024
Publication date: July 2024
Publisher: Universidad UTE (Quito, Ecuador)
Pages: 3540
DOI: https://doi.org/10.29019/enfoqueute.1034
http://ingenieria.ute.edu.ec/enfoqueute/
CFD simulation of the surface pumping of heavy crude oil in eastern Ecuador
José CabreraEscobar1, Guillermo Machado Sotomayor2, Diego Cajamarca Carrazco3, María Magdalena Paredes Godoy4, Raúl CabreraEscobar5
Abstract
The research addresses the study of heat exchange between heavy crude oil and the environment during surface pumping, specifically of a crude oil with an API gravity of 17.5, under the particular atmospheric conditions of Eastern Ecuador. The main objective of the study is to evaluate the temperature loss in a 50meter segment of SCH80 pipe, with a diameter of 4 inches, used for the transportation of heavy crude. This aims to understand how heavy crude loses temperature and to determine the convective coefficient, knowing that the heat loss from the fluid to the environment occurs mainly by convection. This is to determine what the temperature losses will be in longer pipe sections. For this purpose, a methodology using computational fluid dynamics (CFD) simulation was employed, a key tool for predicting the thermal behavior of crude in interaction with the environment. It was determined that the convective coefficient is 52 W/m2.K, and there is a temperature loss of 3.2 K in the 50meter section. With this data, future research could evaluate potential heating technologies that facilitate the transport of heavy crude oil. This approach would allow for exploring innovative solutions to improve efficiency and effectiveness in managing heavy crude, facing one of the main challenges in its transport: managing its high viscosity
Keywords
Heavy crude oil, CFD, simulation.
Resumen
La investigación aborda el estudio del intercambio de calor del crudo pesado con el ambiente en el bombeo superficial, específicamente de un crudo con gravedad API de 17.5, bajo las condiciones atmosféricas particulares del Oriente ecuatoriano. El principal objetivo del estudio es evaluar la pérdida de temperatura en un segmento de 50 metros de tubería SCH80, con 4 pulgadas de diámetro, utilizado para el transporte de crudo pesado. Con esto se busca entender cómo el crudo pesado pierde temperatura y determinar el coeficiente convectivo sabiendo que el calor que pierde el fluido con el ambiente se da principalmente por convección. Esto permitirá determinar cuáles van a ser las pérdidas de temperatura en tramos de tubería mayores. Para ello se empleó una metodología que utiliza la simulación por dinámica de fluidos computacional (CFD), una herramienta clave para predecir el comportamiento térmico del crudo en interacción con el ambiente. Se determinó que el coeficiente convectivo es de 52 W/m2.K, y que existe una pérdida de temperatura de 3.2 K en el tramo de 50 metros. Con esta base de datos se podría en investigaciones futuras evaluar potenciales tecnologías de calefacción que faciliten el transporte del crudo pesado. Este enfoque permitiría explorar soluciones innovadoras para mejorar la eficiencia y efectividad en el manejo del crudo pesado, enfrentando uno de los principales desafíos en su transporte: la gestión de su alta viscosidad.
Palabras Clave
Crudo pesado, CFD, simulación.
I. INTRODUCTION
Heavy crude is characterized by its high density and viscosity. This type of oil is located at the lower end of the classification spectrum according to the API scale, indicating that it is denser and more viscous compared to other types of crude oil [1]. The exploitation and transportation of heavy crude present significant challenges for the oil industry, especially in areas with unique environmental conditions, such as the Ecuadorian east. Due to its high viscosity, heavy crude requires specialized technologies for efficient handling. One of the main challenges is to prevent temperature loss during transport, as this can significantly increase its viscosity and hinder its flow [2]. The viscosity of crude oil, which is closely related to its temperature, has a significant impact on both the efficiency and cost of pumping through pipelines. Higher viscosity requires more energy to move the crude oil, thus increasing operational costs [3].
Cazorla et al. [4] they investigated the impact on production and recovery factor in ultraheavy crude reservoirs through the application of electromagnetic waves in horizontal wells. The main objective was to study how electromagnetic heating can contribute to increasing the production of viscous hydrocarbons, specifically in the Orinoco Oil Belt in Venezuela. A numerical simulation was conducted to analyze the behavior of the reservoir and the effectiveness of this technology in reducing crude viscosity, which could lead to significant improvements in the extraction of these resources.
Tolentino et al. [5] they implemented a novel technology for heating oil fields using nuclear reaction capsules. A steadystate thermal simulation of an oil field similar to Cerro Negro was conducted using COMSOL Multiphysics software with the Finite Element Method. A nuclear reaction capsule was designed with specific dimensions to heat the porous medium, and 11 capsules were distributed throughout the reservoir to enhance the extraction of heavy and extraheavy crude oil, thereby demonstrating the viability and effectiveness of this innovative technology.
Ramírez et al. [3] they carried out the prediction of temperature and viscosity profiles in heavy oil production wells by implementing an induction heater at the bottom of the well. The study focuses on the impact of electric heating on the viscosity of heavy oil, with the goal of improving crude mobility and well efficiency. A CFD model was developed to simulate the thermal and viscosity behavior of heavy oil in the well, demonstrating a significant reduction in the dynamic viscosity of the crude with electric heating. The results validate the effectiveness of induction heating in reducing the viscosity of heavy oil, which could have significant implications for oil production in heavy crude wells.
Pinilla et al. [6] they conducted a Computational Fluid Dynamics study on water production in mature heavy oil fields with horizontal wells. Numerical experiments were performed based on a complete set of 3D NavierStokes equations to simulate multiphase flows in porous media in horizontal heavy oil production wells. A sensitivity analysis was conducted on different rock and fluid properties, such as permeability and oil viscosity, assuming homogeneous porous media. The numerical model was verified with field data using two approaches: one assuming homogeneous rock properties and another using heterogeneous rock properties measured from well logs. The validation results showed that the numerical model could reproduce the observed field production trend with a deviation of 20 % compared to actual data. Additionally, critical hydraulic data in the well, such as pressure and velocity profiles, were obtained, which could be used to enhance production efficiency.
Zhang et al. [7] they conducted a threedimensional numerical simulation to study the heat transfer behavior and flow of waxy crude oil in an inclined pipeline. The main objective was to analyze how the temperature and velocity fields of waxy crude oil vary in pipelines with different inclination angles, providing a theoretical guide for the safe operation of pipelines. Factors such as threedimensional heat transfer, heat transfer in the surrounding soil, and changes in physical parameters with temperature were considered to develop a more accurate and realistic mathematical model.
Kumar et al. [8] they conducted a study using computational fluid dynamics analysis to investigate the flow behavior in the annular core of heat exchangers composed of nanoparticles, specifically employing a heat exchanger made up of nanoparticles. The study focused on the flow of crude oil and water in the annular pipe system, aiming to analyze the heat transfer efficiency, pressure drop, and overall performance of the system. Additionally, the flow characteristics and variations between the crude oil and water were explored, along with the molecular behaviors induced by the central channels. The research also aimed to enhance the understanding of heat generation and fluid flow within the annular core, in order to predict variations in annular flow and improve the control of fluid density within the annular core.
While research has focused on the analysis through CFD simulation of crude oil heating methods, as well as velocity, pressure, and temperature profiles, the heat loss of heavy crude during surface transport has received less attention. This critical aspect of heavy crude handling still requires detailed studies that could significantly improve efficiency in the transport and management of such hydrocarbons. In this context, the current study focuses on the simulation of surface pumping of heavy crude with an API gravity of 17.5, under the specific atmospheric conditions of eastern Ecuador. Using CFD, this research aims to analyze the temperature loss over a 50meter pipeline section, in addition to determining the convective coefficient at which the fluid loses heat, knowing that the main mechanism by which heavy crude loses temperature is through convection with the environment.
This introduction to the problem and the methodological approach adopted lays the groundwork for a broader discussion on potential solutions and the practical implications of the findings, ultimately aiming to move towards more efficient and economically viable management of heavy crude.
II. METHODOLOGY
CFD simulation is widely used in the industry due to its multiple advantages, which include lower costs compared to experimental techniques and the ability to simulate complex physical phenomena. However, CFD simulation must be experimentally validated to confirm that the model configurations and parameters are appropriate for the case study [9], CFD has been widely used in the study of processes where there is mass and heat transfer in the air [10],[11]. In Figure 1, the methodology used for the CFD simulation can be observed.
Fig. 1. Methodology used in the CFD simulation [12].
The model used in the CFD simulation has the same dimensions as the actual pipeline. In Figure 2, you can see the model used in the simulation.
Fig. 2. Model used in the simulation.
For meshing, Tetrahedron elements are used as can be seen in Figure 3. The Orthogonal Quality metric is used to measure the quality of the mesh [13],[14], in the simulation, 501 800 elements were used; this is the optimal number of elements because, beyond this number, the analyzed variable, which is temperature, stabilizes, as can be observed in Figure 4.
The properties of the pipeline and heavy crude can be seen in Table 1 and Table 2, respectively.
Table I. Thermal properties of the pipe
Material properties 
Pipeline SCH80 
Density (kg/m3) 
7850 
Thermal Conductivity (W/(m.k) 
53 
Specific Heat (J/(Kg.k) 
0.49 


Table II. Thermal properties of heavy crude
Fluid properties 
Heavy crude 
Density (Kg/m3) 
908.92 
Specific heat (J/Kg.K) 
1880 
Thermal conductivity (W/m.K) 
0.12 
Viscosity (Kg/m.s) 
0.0798 


The fluid properties are determined from the curves obtained in the laboratory, which are shown in Figure 5, Figure 6, and Figure 7. This is done using an average temperature of 347.04 K.
Fig. 5. Kinematic viscosity of heavy crude oil.
Fig. 6. Dynamic viscosity of heavy crude oil.
Fig. 7. Density of heavy crude oil.
In computational fluid dynamics simulations, phenomena are modeled using mathematical formulas, which are usually partial differential equations. To analyze the surface pumping of heavy crude, the NavierStokes equations are used, which describe the movement of fluids. When considering the fluid as incompressible, two fundamental equations are applied: the continuity equation (1), which ensures mass conservation, and the momentum equation (2), which describes the conservation of momentum [15], [16].
(1)
(2)
where is the density of the fluid, is the velocity field, p is the pressure, µ is the dynamic viscosity of the fluid and are external forces.
To model the heat transfer between the fluid flowing inside the pipeline and the external environment, the fluid is considered incompressible and not involved in chemical reactions. In this context, the heat transfer equation (3) is used, which allows for the calculation of the fluid’s thermal variation in response to environmental conditions and the characteristics of the pipeline.
(3)
where is the specific heat at constant pressure, T is the temperature, k is the thermal conductivity, represents heat generation by viscous dissipation and ST is an additional source term that may include effects such as radiation heat transfer or internal or external heat sources.
The algorithm used for the numerical simulation of the equations is SIMPLE [17], for the spatial discretization of pressure, the method used is PRESTO [18], the flow of heavy crude oil through the pipeline is considered laminar since the Reynolds number is less than 2000. The equation used is equation (4).
(4)
where V is the velocity and D is the diameter of the pipe.
Given the exposure of the pipe to the environment, it is imperative to determine the convective heat transfer coefficient. For this purpose, the Nusselt number is used, a dimensionless magnitude that quantifies heat transport by convection. The correlation between the Nusselt number and the convective coefficient is established through equation (5):
(5)
where Nu is the Nusselt number, h is the convection heat transfer coefficient, L is the length of the pipe and k is the thermal conductivity of the pipe.
The Nusselt number for a flow like this case, which is developed both hydraulically and thermally, is 4.36 [19][20].
The boundary conditions used are as follows:
Inlet: the fluid enters at a speed of 0.73043 m/s, due to the analyzed well producing 2225 barrels and the internal diameter of the pipeline being 0.0848 m. The entering temperature is 348.15 K, a data point taken in the field.
Outlet: for simplicity in the simulation and considering that the objective is to calculate the temperature of the fluid at the outlet, it is assumed to be exposed to atmospheric pressure.
Outside of the pipeline: the exterior of the pipeline is exposed to an ambient temperature of 292.15 K.
To obtain temperature data in the field, bimetallic thermometers model E from Ashcroft were used, as shown in Figure 8.
Fig. 8. Ashcroft Bimetallic Thermometer.
III. RESULTS
For the validation of the simulation, the experimental outlet temperature of 345.93 K and the CFD simulation result of 344.95 K were used. With these values, a Mean Absolute Percentage Error (MAPE) of 0.28 % was obtained. Additionally, to validate the convective coefficient value obtained from the simulation of 52 W/m².K, it was compared with the analytically derived value of 52.9 W/m².K, resulting in a MAPE of 1.7 %. For a simulation to be considered a good approximation, the MAPE should not exceed 2 % [21], [22].
In Figure 9, the residuals are presented, which are essential for evaluating the convergence of numerical simulations. The residuals represent the difference between the calculated values and the expected values of the flow variables. For this specific case, an absolute convergence criterion of 1e06 is used, ensuring adequate accuracy in the simulation results.
Fig. 9. Simulation Residual Curves.
In Figure 10, the average temperature curve at the pipeline outlet can be observed. The temperature at which the simulation converges is 344.95 K, which is reached at iteration 293.
Fig. 10. Temperature at the Pipeline Outlet.
In Figure 11, the temperature profile at the pipeline inlet can be observed, representing the highest temperature recorded in this simulation. Specifically, the crude oil temperature at the inlet of the analyzed section is 348.15 K. It can be seen that the closer the heavy crude is to the surface of the pipeline, the lower the temperature.
Fig. 11. Temperature Profile at the Pipeline Inlet.
In Figure 12, the temperature profile at the pipeline outlet can be observed. This area has the lowest temperature recorded in this simulation. Specifically, the crude oil temperature at the outlet of the analyzed section is 344.95 K. Similar to the start of the analyzed section, it can be seen that the closer the heavy crude is to the surface of the pipeline, the lower the temperature.
Fig. 12. Temperature Profile at the Pipeline Outlet.
IV. DISCUSSION
From the experimental and simulation results, it can be observed that there is a temperature reduction from 348.15 K at the starting point of the analysis to 344.95 K at the endpoint of the 50 m pipeline section. This temperature decrease is due to the ambient air temperature being 292.15 K. The primary reason for the heat loss is convection with the environment. The convective coefficient determined through the simulation is 52 W/m²·K. This coefficient can be used to evaluate the temperature loss over long pipeline sections, allowing for the determination of the heat that needs to be added to the fluid. Temperature loss in the pipeline during transport significantly affects the viscosity of heavy crude, which in turn influences the efficiency and costs associated with pumping it over long distances.
The simulation identifies that the variables most influencing this heat loss are the ambient temperature, the crude temperature, the convective heat transfer coefficient, and the crude velocity. Other parameters that also influence but cannot be modified include the area exposed to the environment, and the thermal properties of the pipeline and the fluid. Ambient temperature affects heat loss because, according to Newton’s law of cooling, the rate of heat loss from the pipeline is proportional to the temperature difference between the pipeline and the environment. Crude temperature is important in heat loss because a higher temperature implies a greater potential for heat loss to a cooler environment. The convective heat transfer coefficient is crucial in heat loss, as the loss occurs primarily through convection. The convective transfer coefficient is directly proportional to the heat loss. The crude flow velocity is critical since a higher velocity implies that the crude remains in contact with the environment for less time, resulting in lower heat losses. Conversely, a lower velocity means that the crude is exposed to the environment for a longer time, increasing the heat losses.
V. CONCLUSIONS AND RECOMMENDATIONS
Through the application of computational fluid dynamics, it was possible to quantify the temperature variation along a 50meter segment of heavy crude oil transportation pipeline, determining that there is a 3.2 K reduction in the analyzed section. Additionally, the convective coefficient was found to be 52 W/m²·K. This analysis provided a detailed understanding of the temperature profiles and thermal dynamics involved, which is essential for optimizing transportation processes and ensuring the operational integrity of the system. The results obtained offer valuable insights for the design and improvement of thermal maintenance strategies, significantly contributing to the efficiency and safety in the management of heavy crude transportation.
The results obtained serve as a foundation for future studies that introduce heating technologies for heavy crude to maintain its viscosity within optimal ranges for efficient transportation.
ACKNOWLEDGEMENT
The authors would like to thank the Baker Hughes Company for all the facilities they provided for this research.
1 José CabreraEscobar. Email: omar.cabrera@unach.edu.ec, Orcid: https://orcid.org/0000000201975163, Filiation: Universidad Nacional de Chimborazo.
2 Guillermo Machado Sotomayor. Email: gmachado@unach.edu.ec, Orcid: https://orcid.org/000000015226468X, Filiation: Universidad Nacional de Chimborazo.
3 Diego Cajamarca Carrazco. Email: diego.cajamarca@espoch.edu.ec, Orcid: https://orcid.org/0000000311578900, Filiation: Escuela Superior Politécnica de Chimborazo.
4 María Magdalena Paredes Godoy. Email: maparedes@unach.edu.ec, Orcid: https://orcid.org/0000000282110400, Filiation: Universidad Nacional de Chimborazo.
5 Raúl CabreraEscobar. Email: rvce0001@red.ujaen.es, Orcid: https://orcid.org/0009000585897578, Filiation: Universidad de Jaén.
References
[1] J. G. Speight, “Chapter 1occurrence and formation of crude oil and natural gas,” Subsea and Deepwater Oil and Gas Science and Technology, pp. 1–43, 2015.
[2] C. Zou, “Heavy oil and bitumen,” Unconventional petroleum geology, pp. 345–370, 2017.
[3] J. Ramírez, A. Zambrano and N. Ratkovich, “Prediction of Temperature and Viscosity Profiles in HeavyOil Producer Wells Implementing a Downhole Induction Heater,” Processes, vol. 11, no. 2, 2023, https://doi.org/10.3390/pr11020631
[4] Cazorla and J. J., “Evaluación de la producción y el factor de recobro en yacimientos de crudo extra pesado a través de la aplicación de ondas electromagnéticas en pozos horizontales,” Sep. 2013, http://saber.ucv.ve/handle/10872/4261
[5] S. L. Tolentino, S. Caraballo, Á. Duarte and J. Mendoza, “Calentamiento de yacimientos petrolíferos mediante cápsulas de reacción nuclear,” Universidad, Ciencia y Tecnología, vol. 20, no. 79, pp. 58–68, 2016, Online]. Available https://bit.ly/3WT8qgS
[6] A. Pinilla, M. Asuaje, C. Pantoja, L. Ramirez, J. Gomez and N. Ratkovich, “CFD study of the water production in mature heavy oil fields with horizontal wells,” PLoS One, vol. 16, no. 10, pp. 1–24, Apr. 2021, https://doi.org/10.1371/journal.pone.0258870
[7] L. Zhang, C. Du, H. Wang and J. Zhao, “Threedimensional numerical simulation of heat transfer and flow of waxy crude oil in inclined pipe,” Case Studies in Thermal Engineering, vol. 37, p. 102237, 2022, https://doi.org/10.1016/j.csite.2022.102237
[8] R. Kumar R et al., “Investigation of nano composite heat exchanger annular pipeline flow using CFD analysis for crude oil and water characteristics,” Case Studies in Thermal Engineering, vol. 49, p. 103297, 2023, https://doi.org/10.1016/j.csite.2023.103297
[9] Q. Tu, Z. Ma and H. Wang, “Investigation of wet particle drying process in a fluidized bed dryer by CFD simulation and experimental measurement,” Chemical Engineering Journal, vol. 452, p. 139200, 2023, https://doi.org/10.1016/j.cej.2022.139200.
[10] E. Getahun, M. A. Delele, N. Gabbiye, S. W. Fanta, P. Demissie and M. Vanierschot, “Importance of integrated CFD and product quality modeling of solar dryers for fruits and vegetables: A review,” Solar Energy, vol. 220, pp. 88–110, 2021, https://doi.org/10.1016/j.solener.2021.03.049
[11] P. D. Tegenaw, M. G. Gebrehiwot and M. Vanierschot, “On the comparison between computational fluid dynamics (CFD) and lumped capacitance modeling for the simulation of transient heat transfer in solar dryers,” Solar Energy, vol. 184, pp. 417–425, 2019, https://doi.org/10.1016/j.solener.2019.04.024.
[12] José CabreraEscobar, David Vera Francisco Jurado, Manolo CórdovaSuárez, Gonzalo SantillánValdiviezo, Antonio RodríguezOrta and R. CabreraEscobar, “Optimization of olive pomace dehydration process through the integration of computational fluid dynamics and deep learning,” Energy Sources, Part A: Recovery, Utilization, and Environmental Effects, vol. 46, no. 1, pp. 4756–4776, 2024, https://doi.org/10.1080/15567036.2024.2331563
[13] ANSYS Inc, Ansys Fluent Tutorial Guide, 2023 R1. Canonsburg, 2023.
[14] ANSYS Inc, “ANSYS Meshing User’s Guide,” 2024. [Online]. Available: http://www.ansys.com
[15] K. Garg, S. Singh, M. Rokade and S. Singh, “Experimental and computational fluid dynamic (CFD) simulation of leak shapes and sizes for gas pipeline,” J Loss Prev Process Ind, vol. 84, p. 105112, 2023, https://doi.org/10.1016/j.jlp.2023.105112
[16] Y. Xiao et al., “Bionic pipeline transport characteristics with transverse protuberances in slurry shield circulation system based on CFDDEM,” Powder Technol, vol. 432, p. 119133, 2024, https://doi.org/10.1016/j.powtec.2023.119133.
[17] A. Benhamza, A. Boubekri, A. Atia, T. Hadibi and M. Arıcı, “Drying uniformity analysis of an indirect solar dryer based on computational fluid dynamics and image processing,” Sustainable Energy Technologies and Assessments, vol. 47, p. 101466, 2021, https://doi.org/10.1016/j.seta.2021.101466.
[18] ANSYS Inc, Ansys Fluent Theory Guide, 2023 R1. Canonsburg, 2023.
[19] T. L. Bergman, A. S. Lavine, F. P. Incropera and D. P. DeWitt, Introduction to Heat Transfer. Wiley, 2011. [Online]. Available: https://books.google.com.ec/books?id=YBaNaLurTD4C
[20] A. J. Ghajar, D. Yunus. and A. Cengel, Heat and Mass Transfer: Fundamentals and Applications. McGrawHill Education, 2014. [Online]. Available: https://books.google.com.ec/books?id=B89MnwEACAAJ
[21] F. Hussain, M. Jaskulski, M. Piatkowski and E. Tsotsas, “CFD simulation of agglomeration and coalescence in spray dryer,” Chem Eng Sci, vol. 247, p. 117064, 2022, https://doi.org/10.1016/j.ces.2021.117064
[22] J. CabreraEscobar, D. Vera, F. Jurado and R. CabreraEscobar, “CFD investigation of the behavior of a solar dryer for the dehydration of olive pomace,” Energy Sources, Part A: Recovery, Utilization, and Environmental Effects, vol. 46, no. 1, pp. 902–917, Dec. 2024, https://doi.org/10.1080/15567036.2023.2292242