Universidad Mayor de San Simã - N
Universidad Mayor de San Simã - N
Universidad Mayor de San Simã - N
PRÁCTICA N°3
INTEGRANTES:
Cossio Quiñones Angelli Prima
Cruz Guisada Riosmar
Cuqui Inca Thalia
Duran Cadima Diana Michelle
Felix Cari Gary Edgar
Fernandez Rocha Nohelia
Gallego Villca Ivonne Lucia
Guzman Lutina Jhoseline Alejandra
Huanca Alcazar Mariela Lizbeth
Illanes Martinez Jhoelma
Illanes Vargas Estefany
DOCENTE:
MATERIA:
GRUPO:
2
Cochabamba – Bolivia
1
ÍNDICE GENERAL
1. RESUMEN 7
2. INTRODUCCIÓN 7
3. OBJETIVOS 8
-OBJETIVO GENERAL 8
- OBJETIVOS ESPECÍFICOS 8
4. MARCO TEÓRICO 8
4.1. MÉTODOS DE CÁLCULOS DE LAS PÉRDIDAS POR FRICCIÓN: 9
4.1.1. ECUACIÓN GENERAL DE ENERGÍA. 9
4.1.2. TEOREMA DE BERNOULLI. 10
4.1.3. PÉRDIDAS DE ENERGÍA EN UN SISTEMA A PRESIÓN DE UN FLUJO
INCOMPRESIBLE. 10
4.1.4. FLUJO LAMINAR. 11
4.1.5. FLUJO TURBULENTO. 11
4.1.6. NÚMERO DE REYNOLDS. 12
4.1.7. PÉRDIDAS DE CARGA POR FRICCIÓN EN TUBERÍAS 12
4.1.8. ECUACIÓN DE DARCY-WEISBACH. (MATAIX, 1986) 13
4.1.9. COEFICIENTE DE FRICCIÓN (Ƒ) O FACTOR DE PÉRDIDAS POR FRICCIÓN
13
4.1.10. FACTOR DE PÉRDIDAS MENORES 13
4.1.11. EL VALOR DE LOS COEFICIENTES DE PÉRDIDAS MENORES 14
4.1.12. RUGOSIDAD RELATIVA 14
4.1.12.1. RUGOSIDAD ABSOLUTA EN FUNCIÓN DE LOS AÑOS DE SERVICIO. 15
4.1.13. ECUACIÓN DE HAZEN-WILLIAMS 15
4.1.14. PÉRDIDAS SEGÚN MANNING 15
5. PARTE EXPERIMENTAL 16
5.1. DISEÑO Y MONTAJE DEL EQUIPO HIDRÁULICO PARA LA
EXPERIMENTACIÓN DE PÉRDIDAS POR FRICCIÓN, PERDIDAS LOCALIZADAS 16
5.1.1. OBJETIVO ESPECÍFICO 16
5.1.2. MARCO TEÓRICO 16
5.1.2.1. TIPOS DE FLUJOS EN TUBERÍAS 16
5.1.2.2. ECUACIONES GENERALES DE LA MECÁNICA DE FLUIDOS 16
5.1.2.3. PÉRDIDAS DE ENERGÍA POR FRICCIÓN. 17
5.1.2.3.1. FLUJO TURBULENTO EN CONDUCTOS DE PRESIÓN CIRCULAR. 18
5.1.2.4.PERDIDA DE CARGAS EN ACCESORIOS. 18
5.1.2.4.1. ACCESORIO DE TUBERÍAS. 18
5.1.2.4.1.1. BRIDAS. 19
5.1.2.4.1.2. CODOS. 19
5.1.2.4.1.3. VÁLVULAS. 19
5.1.2.4.1.4. VÁLVULA DE BOLA. 20
5.1.3. METODOLOGÍA DEL DISEÑO Y MONTAJE DEL EQUIPO HIDRÁULICO 20
5.1.3.1. ESTUDIOS PREVIOS. 20
2
5.1.3.2. DISEÑO DEL BANCO DE PRUEBAS. 20
5.1.3.2.1. DESCRIPCIÓN DEL SISTEMA. 21
5.1.3.2.2. CONSTRUCCIÓN DEL EQUIPO 23
5.1.4. INSTRUCCIONES 23
5.1.4.1. PROCEDIMIENTOS DE EMPLEO. 23
5.1.4.1.1. MANTENIMIENTO PREVENTIVO. 25
5.1.6. CÁLCULOS Y RESULTADOS DEL DISEÑO: 30
CALCULO DE BOMBA: Pteórico=HDT**Qt 30
5.1.7. CONCLUSIONES: 35
5.2.PÉRDIDAS DE ENERGÍA EN CRUCES DE TUBERÍAS 35
5.2.1. OBJETIVO ESPECÍFICO 35
5.2.2. MARCO TEÓRICO ESPECÍFICO 35
5.2.2.1. RESISTENCIA AL FLUJO EN TUBERÍAS: PÉRDIDAS LOCALES 35
5.2.2.2. PRIMER MÉTODO: FÓRMULA GENERAL DE LAS PÉRDIDAS
SECUNDARÍAS 36
5.2.2.2.1. PÉRDIDA POR ENTRADA 37
5.2.2.2.2. PÉRDIDA POR REJILLA 37
5.2.2.2.3. PÉRDIDA POR AMPLIACIÓN 38
5.2.2.2.4. PÉRDIDA POR REDUCCIÓN 39
5.2.2.2.5. PÉRDIDA POR CAMBIO DE DIRECCIÓN 40
5.2.2.2.6. PÉRDIDA EN VÁLVULAS 41
5.2.2.3. SEGUNDO MÉTODO: LONGITUD DE TUBERÍA EQUIVALENTE 42
5.2.2.4. PÉRDIDAS DE ENERGÍA EN CRUCES DE TUBERÍAS 42
5.2.2.5. VALORES DEL COEFICIENTE K 43
5.2.3. MÉTODO Y MATERIALES 45
5.2.3.1. MATERIALES 45
5.2.3.2. METODOLOGÍA 45
5.2.3.2.1 Diseño de dispositivo experimental 45
5.2.3.2.2. Diseño de las mediciones 48
5.2.4. ANÁLISIS DE RESULTADOS 50
5.2.4.1. NOMENCLATURA 50
5.2.4.2. ESTIMACIÓN DE K 50
5.2.4.2.1. FLUJO PERPENDICULAR 51
5.2.4.2.2. EXPRESIONES PARA ESTIMAR K3 51
5.2.4.2.3. EXPRESIONES PARA ESTIMAR K4 52
5.2.4.2.4. RELACIÓN DE LOS COEFICIENTES KI CON OTROS PARÁMETROS
HIDRÁULICOS 53
5.2.4.2.5. LÍNEAS DE TENDENCIA 57
5.2.5. CONCLUSIONES 58
5.3. PÉRDIDAS POR FRICCIÓN EN LA CONDUCCIÓN DE AGUA CON CONTENIDO
DE SÓLIDOS 60
5.3.1 OBJETIVOS ESPECÍFICOS: 60
5.3.2 MARCO TEÓRICO ESPECÍFICO 60
5.3.2.1 Aplicación de propiedades reológicas al cálculo de las pérdidas de carga 61
3
5.3.2.2 La presencia de las partículas sólidas como causante de variación en el
comportamiento hidráulico 65
5.3.3 MÉTODO 66
5.3.3.1 Aproximación ingenieril y verificación con datos de campo tomados en el
acueducto Chapala-Guadalajara 66
5.3.3.2 Descripción del Acueducto Chapala – Guadalajara 68
5.3.3.3 Mediciones y resultados 69
5.3.4 ANÁLISIS DE RESULTADOS 70
5.3.5 CONCLUSIONES 71
5.4. PÉRDIDAS DE ENERGÍA DEBIDAS A LA FRICCIÓN 73
5.4.1. OBJETIVO ESPECÍFICO 73
5.4.2. MARCO TEÓRICO ESPECÍFICO 73
5.4.2.1. LA ECUACIÓN DE DARCY-WEISBACH 73
5.4.2.2 PÉRDIDAS DE FRICCIÓN EN FLUJO LAMINAR: 74
5.4.2.3 PÉRDIDAS DE ENERGÍA EN FLUJO TURBULENTO 75
5.4.2.3.1 DIAGRAMA DE MOODY 76
5.4.2.4 PÉRDIDAS DE ENERGÍA EN LA ZONA CRÍTICA 77
5.4.2.5.FÓRMULA DE HAZEN-WILLIAMS 78
5.4.2.6 PROGRAMACIÓN EN CASIO BASIC 79
5.4.2.7 ESTRUCTURA SECUENCIAL PARA PROGRAMACIÓN EN BASIC 81
5.4.2.8 DESARROLLO DE UN PROGRAMA PARA EL CALCULO DE LAS
PERDIDAS POR FRICCION 81
5.4.2.8.1 CONFIGURACIÓN DEL SISTEMA DE ADQUISICIÓN DE DATOS: 82
5.4.2.8.2. EJEMPLO DE CÁLCULO DE LA CAÍDA DE PRESIÓN EN UN
INTERCAMBIADOR DE TUBO Y CORAZA 83
5.4.2.8.3. EJEMPLO DEL CÁLCULO DEL CAUDAL EN TUBERÍAS EN SERIE 84
5.4.3. MÉTODO EXPERIMENTAL 86
5.4.3.1.AUTOMATIZACIÓN DEL PROCESO DE CÁLCULO DE LA PÉRDIDA DE
ENERGÍA POR FRICCIÓN. 86
5.4.3.2 OBTENCION DE LA PERDIDA DE ENERGIA DEBIDAS A LA FRICCION EN
TUBERIAS 87
5.4.4. ANÁLISIS DE RESULTADOS 99
5.4.4.1 CODIFICACION DEL PROGRAMA EN LENGUAJE BASIC PARA CASSIO
FX-880P 99
5.4.4.2. TOMA DE DATOS DEL SENSOR MEDIANTE FIELD POINTS Y TRATADOS
CON LABVIEW ANALIZADOS EN TIEMPO REAL 99
5.4.5. CONCLUSIONES 101
5.5. DESARROLLO DE UN ALGORITMO Y PROGRAMA EN MATLAB PARA
SISTEMATIZAR Y AUTOMATIZAR UN MODELO DIGITAL HIDRODINÁMICO EN
ESTADO ESTABLE DE UN OLEODUCTO PARA TRANSPORTE DE CRUDO LIVIANO
101
5.5.1. OBJETIVO ESPECÍFICO 101
5.5.2. MARCO TEÓRICO ESPECÍFICO 101
5.5.2.1 HIDRODINÁMICA DE OLEODUCTOS: 101
5.5.2.2 CARACTERÍSTICAS DE LA TUBERÍA: 101
5.5.2.3 PROPIEDADES FÍSICAS DEL FLUIDO 104
4
5.5.2.4 PARÁMETROS A CONSIDERAR EN EL DISEÑO DE OLEODUCTOS: 106
5.5.2.5 DETERMINACIÓN DE LA CAÍDA DE PRESIÓN EN OLEODUCTOS 106
5.5.2.6 PÉRDIDA DE PRESIÓN POR FRICCIÓN USANDO LA ECUACIÓN DE
DARCY 107
5.5.2.7 MÁXIMA PRESIÓN DE OPERACIÓN ADMISIBLE (MAOP) 107
5.5.3. MÉTODO Y MATERIALES 108
5.5.3.1 DATOS INICIALES CONOCIDOS Todos estos parámetros son datos iniciales
conocidos para el desarrollo de este algoritmo. 108
5.5.3.2 CONDICIONES DE DISEÑO 109
5.5.3.4 RUTA DEL OLEODUCTO 109
5.5.3.5 PERFIL TOPOGRÁFICO DE LA RUTA 109
5.5.3.6 PARÁMETROS DE DISEÑO 110
5.5.3.7 ALGORITMO DE DISEÑO 111
5.5.3.8 CÁLCULO DE LA PRESIÓN MÁXIMA ADMISIBLE DE OPERACIÓN 112
5.5.3.9 HIDRODINÁMICA DE OLEODUCTOS: 112
5.5.3.9.1 CÁLCULO DE LA VELOCIDAD DE FLUJO 112
5.5.3.9.2 CÁLCULO DEL NÚMERO DE REYNOLDS ( RE ) 113
5.5.3.9.3 DETERMINACIÓN DEL FACTOR DE FRICCIÓN (f) 113
CÁLCULO DE PRESIÓN DE ENTRADA REQUERIDA EN EL SISTEMA (PRESIÓN
MÁXIMA NECESARIA EN LA SALIDA DE LA ESTACIÓN “A”) (PA) 114
5.5.3.9.5 CONSIDERACIÓN DE ALTERNATIVAS 116
5.5.2.10 PRIMERA ETAPA DE CÁLCULO 116
5.5.3.11 SEGUNDA ETAPA DE CÁLCULO 120
5.5.4. ANÁLISIS DE RESULTADOS 123
5.5.4.1 PRIMERA ETAPA A-B 123
5.5.4.2 SEGUNDA ETAPA C-B 123
5.5.5. CONCLUSIONES 124
5.6. OPTIMIZACIÓN DEL DESEMPEÑO ENERGÉTICO PARA EL TRANSPORTE DE
HIDROCARBUROS EN EL OLEODUCTO DE LOS LLANOS ORIENTALES MEDIANTE
LA UTILIZACIÓN DEL AGENTE REDUCTOR DE FRICCIÓN EXTREME POWER 125
5.6.1. OBJETIVOS ESPECÍFICOS 125
5.6.2. MARCO TEÓRICO 125
5.6.2.1. Oleoducto 125
5.6.2.2. Reducción de la fricción o reducción de arrastre 125
5.6.2.3. Agentes reductores de arrastre (Drag Reducing Agent DRA) 126
5.6.2.4. Agente reductor de arrastre Extreme Power 126
5.6.2.5. Viscosidad 127
5.6.2.6. Pérdidas por fricción 127
5.6.2.7. Turbulencia 128
5.6.3. METODOLOGÍA Y MATERIALES 128
5.6.3.1. Metodología 128
5.6.3.2. Caracterización y evaluación del crudo con el agente reductor de fricción a
las diferentes concentraciones propuestas 129
5.6.3.3. Caracterización reológica 129
5.6.3.4. Pruebas de flujo 129
5
5.6.3.5. Análisis del efecto del DRA a la temperatura de transporte 130
5.6.3.6. Análisis de fluidez de los sistemas y efectos de los DRA 130
5.6.3.7. Prueba de fluidez con Copa Ford (Vidrio) 131
5.6.3.8. Prueba de fluidez con Copa Ford (Acero) 131
5.6.3.9. Planteamiento y evaluación de las pérdidas de fricción mediante la
simulación analítica de las pérdidas por fricción a condiciones de laboratorio en
ASPEN HYSYS 132
5.6.4. ANÁLISIS DE RESULTADOS 139
5.6.4.1. Análisis experimental 139
5.6.4.2. Caracterización del crudo mediante pruebas de flujo 139
5.6.4.3. Análisis del agente reductor a la temperatura de transporte variando la
concentración 142
5.6.4.4. Análisis de prueba de fluidez en Copa Ford (vidrio) 144
5.6.4.5. Análisis de prueba de fluidez en Copa Ford (acero) 144
5.6.4.6. Simulación 146
5.6.5. CONCLUSIONES 148
6. BIBLIOGRAFÍA 149
6
PÉRDIDAS POR FRICCIÓN
1. RESUMEN
2. INTRODUCCIÓN
Para la ingeniería el estudio de los fluidos siempre ha sido una pieza importante
para comprender todos los fenómenos que afectan el desempeño de una máquina y
el cual está estrechamente relacionado con las pérdidas que dentro de ella se
puedan presentar. Como muchas de las máquinas con las que un ingeniero se
encuentra están conformadas por un sistema de tuberías es vital comprender cuales
son los factores que van a ocasionar pérdidas en un proceso.
El estudio de las pérdidas por fricción es la pérdida de energía mecánica que se
disipa en forma de calor por las paredes de los conductos. También es la energía
mínima que se debe suministrar al sistema de flujo para que el fluido se mueva con
una velocidad o flujo constante.
Es por ello que en este trabajo nos centramos en las pérdidas ocasionadas por la
fricción entre el fluido y la pared del tubo al igual que por los accesorios que
conforman una red de tuberías, teniendo en cuenta la naturaleza del fluido. De esta
manera mejorar nuestras habilidades como ingenieros y ser los más certeros a la
hora de seleccionar un buen material para una tubería y los respectivos accesorios
que la van a conformar, teniendo en cuenta que dichas selecciones deben satisfacer
los requerimientos técnicos de un proyecto y que permitan lograr el objetivo final que
se desee.
7
3. OBJETIVOS
-OBJETIVO GENERAL
● Estudiar y analizar las distintas pérdidas por fricción que se dan o producen
en un sistema de transporte de fluidos por tuberías a través de cálculos y
simulaciones mediante la recopilación de seis artículos seleccionados.
- OBJETIVOS ESPECÍFICOS
4. MARCO TEÓRICO
A medida que un fluido fluye por un conducto, tubo o algún otro dispositivo, ocurren
pérdidas de energía debido a la fricción que hay entre el líquido y la pared de la
tubería o los accesorios; tales pérdidas de energías traen como resultado una
disminución de la presión entre dos puntos del sistema de flujo.
8
En estructuras largas, las pérdidas por fricción son muy importantes, por lo que ha
sido objeto de investigaciones teórico-prácticas para llegar a soluciones
satisfactorias de fácil aplicación.
Las propiedades de los fluidos más relevantes para el estudio de pérdidas de carga
en sistemas de tuberías son: Densidad, Peso específico, Densidad relativa,
Viscosidad cinemática,Viscosidad dinámica, Caudal y Presión.
E Total = (v2/2g) + P γ + z
9
Esta ecuación indica la energía total.
10
El planteamiento de la ecuación de energía, se tiene que:
Es aquel en el cual las partículas del fluido siguen trayectorias paralelas, formando
junto de ellas capas o láminas. La velocidad de estas partículas es mayor cuando
están más alejadas de las paredes del conducto, o sea que la velocidad de dichas
partículas está en función de la distancia de las paredes del conducto.
Es aquel en el cual las partículas del fluido no siguen trayectorias paralelas, es decir
que se mueven en forma desordenada en todas las direcciones. Es imposible
conocer la trayectoria de una partícula individualmente.
11
4.1.6. NÚMERO DE REYNOLDS.
Puede considerarse como la relación de las fuerzas dinámicas de la masa del fluido
respecto a los esfuerzos de deformación ocasionados por la viscosidad.
Re = v ∗ D ∗ ρ/ μ = v ∗ D/θ
Tipo de Flujo Re
La pérdida de energía por fricción es la debida al rozamiento del fluido con las
paredes de la tubería o del conducto. Esta pérdida que continúa en la dirección del
flujo, puede resultar considerable en tramos largos, y por el contrario, puede ser
despreciable en tramos cortos.
12
4.1.8. ECUACIÓN DE DARCY-WEISBACH. (MATAIX, 1986)
Además de las pérdidas por fricción en tuberías, existe otra clase conocidas como
pérdidas menores, que ocurren en las redes hidráulicas como consecuencia de
accesorios presentes en las tuberías. En algunos casos sendas pérdidas llegan a
ser mucho más significativas que las producidas por fricción en la tubería.
13
4.1.11. EL VALOR DE LOS COEFICIENTES DE PÉRDIDAS MENORES
a. Para flujo laminar (0 ≤ Re ≤ 2000) en todas las tuberías y para cualquier fluido,
Hagen y Poiseville determinaron la siguiente ecuación:
f =64/Re
14
4.1.12.1. RUGOSIDAD ABSOLUTA EN FUNCIÓN DE LOS AÑOS DE SERVICIO.
εt = ε0 + at
hf = J ∗ Le
Aunque su empleo es muy generalizado, se debe tener en cuenta que existen varias
restricciones para su utilización y que de no considerarse llevarán a
sobredimensionar la tubería.
15
5. PARTE EXPERIMENTAL
16
diferencial, dependiendo del problema. A este conjunto de ecuaciones dadas en
su forma diferencial también se le denomina ecuaciones de Navier-Stokes.
A medida que un fluido fluye por un conducto, tubo o algún otro dispositivo,
ocurren pérdidas de energía debido a la fricción; tales energías traen como
resultado una disminución de la presión entre dos puntos del sistema de flujo.
(5.3.1)
(5.3.2)
(5.3.3)
Las dos expresiones anteriores muestran que la descarga varía directamente con la
diferencia de las presiones el diámetro del conducto y la aceleración de la gravedad,
e inversamente con la longitud de la tubería y la viscosidad cinemática.
17
5.1.2.3.1. FLUJO TURBULENTO EN CONDUCTOS DE PRESIÓN CIRCULAR.
(5.3.7)
(5.3.8)
Además de las pérdidas de energía por fricción, hay otras pérdidas menores
asociadas con los problemas en tuberías. Se considera que tales pérdidas ocurren
localmente en el disturbio del flujo. La presencia de llaves de paso,
ensanchamientos, codos.
• Bridas
18
• Codos
• Cuellos o acoples
• Válvulas
• Tornillos o niples
5.1.2.4.1.1. BRIDAS.
5.1.2.4.1.2. CODOS.
Son accesorios de forma curva que se utilizan para cambiar la dirección del flujo.
5.1.2.4.1.3. VÁLVULAS.
Es un accesorio que se utiliza para regular y controlar el fluido de una tubería. Las
válvulas pueden ser de varios tipos según sea el diseño del cuerpo y el movimiento
del obturador.
19
5.1.2.4.1.4. VÁLVULA DE BOLA.
El cuerpo de la válvula tiene una cavidad interna esférica que alberga un obturador
en forma de bola o esfera. La bola tiene un corte adecuado (usualmente en V) que
fija la curva característica de la válvula, y gira transversalmente.
Para la investigación del proyecto se pactaron una serie de visitas, las cuales se
analizaron el tipo de montaje, se contaron con asesorías y guías sobre el montaje y
el proceso de funcionamiento del mismo, adicionalmente se contó con asesoría
institucional.
Este montaje se encuentra anclado a una base de acero, la cual fue diseñada
teniendo en cuenta la posibilidad de trasladarlo y que no solo cumpliera su objetivo
principal de sostener el montaje sino que facilitará un lugar a los practicantes de
tomar los apuntes de la experiencia en desarrollo.
20
5.1.3.2.1. DESCRIPCIÓN DEL SISTEMA.
21
Tramo 1: Este es un tramo no experimentable en el cual se impulsa, regula y se
mide el caudal del flujo. Cuenta con un tanque de almacenamiento redondo con
unas dimensiones de altura 0.95 mts y diámetro 0.55 mts, estas dimensiones
garantizan que no se vea afectada la altura o carga de succión, también cuenta con
su respectivo desagüe y su retorno. Un equipo de bombeo que consta de una
bomba centrifuga marca BARNES de ½ HP, con un caudal máximo de 170 Lts/min y
una velocidad angular ω= 3450; dicha bomba es impulsada por un motor eléctrico
marca Siemens directamente acoplado a la bomba, monofásico rebobina a trifásico
de ½ HP. Una válvula de control de flujo, tipo cortina, para controlar el volumen de
flujo en el sistema.
Tramo 4: En este se estudia los accesorios para cambio de direcciones en las redes
hidráulicas, como los semi codos marca Pavco. Cuenta también para
experimentación una válvula Bola de hierro de diámetro ¾ de pulgada marca Italy,
este tramo cuenta con su respectiva válvula bola de cierre.
Tramo 5: En este se estudia las pérdidas sufridas por los cambios bruscos de
diámetro por intermedio tubería de ¾ de pulgada en tubería de PVC a 1 ¼ de
pulgada en tubería de PVC marca Pavco y los sufridos por los cambios de dirección
por intermedio de un codo de ¾ de pulgada, este tramo cuenta con su respectiva
válvula bola de cierre de tramo.
Tramo 6: Último tramo del sistema y es en donde se podrá estudiar las pérdidas
generadas en elementos como válvulas globo y compuerta de entrada de ¾ de
pulgada, en cobre. Este tramo cuenta con su respectiva válvula bola de cierre de
tramo y medidores de presión. En este tramo también se encuentra la tubería de
retorno al tanque de almacenamiento.
22
5.1.3.2.2. CONSTRUCCIÓN DEL EQUIPO
● Base de acero: esto fue realizado por un experto del área con el cual se hizo
un replanteo mínimo para dar más estabilidad de la estructura.
● Sistema de tuberías: con apoyo técnico se realizó la construcción y
ensamblaje de las tuberías.
● Anclaje a tablero: una vez terminado esto se realiza en anclaje con
abrazaderas y tapones de caucho para ajustarlo mejor, se instalaron las
válvulas de purga en los extremos superiores del banco. Se realiza
adecuación de tanque de reserva de agua.
● Acople la bomba y verificación del equipo: se realiza el acople a la bomba e
instalación eléctrica. Se verifica que no haya fugas ni mal funcionamiento del
equipo. Se realiza ensayo sobre el banco.
5.1.4. INSTRUCCIONES
23
● Se conecta la Bomba centrífuga de ½ HP. Se arranca desde el suich parte
posterior (botón Verde – on).
● Se abren todas las válvulas en el tablero con la finalidad de purgar la tubería
con el objeto de eliminar las partículas de aire y gas existentes en la tubería.
● Se abren las válvulas de purga ubicadas en parte superior del tablero y se
deja fluir el agua para terminar de purgar el sistema.
● Una vez purgada la tubería se procede a cerrar las válvulas de purga que se
encuentran en la parte superior del tablero.
● Purgado el sistema, se procede a seleccionar el tramo sobre el cual se
realizará el ensayo. Para ello se cierran las válvulas bolas 2,3,4,5.
● Para estudiar el tramo 2, se abre válvula bola 1, se mantiene cerradas
válvulas bola 2, 3, 4, 5. Se realiza lectura de caudal y de manómetros en el
tramo.
● Para estudiar el tramo 3 (PVC), se abre válvula bola 2, se mantiene cerradas
válvula bola 3, 4,5. Se cierra válvula bola 1. Se realiza lectura de caudal y de
manómetros en el tramo.
● Para estudiar el tramo 4(Semi codo), se abre válvula bola 3, se mantiene
cerradas válvula bola 1, 4,5. Se cierra válvula bola 2. Se realiza lectura de
caudal y de manómetros en el tramo.
● Para estudiar el tramo 5 (expansión), se abre válvula bola 4, se mantiene
cerradas válvulas bola 1, 2, 5. Se cierra válvula bola 3. Se realiza lectura de
caudal y de manómetros en el tramo.
● Para estudiar el tramo 6 (válvulas), se abre válvula bola 5, se mantiene
cerradas válvula bola 1, 2,3. Se cierra válvula bola 4. Se realiza lectura de
caudal y de presión en los manómetros del tramo. Las válvulas 3, 4 deben
estar abiertas. Se recomienda solo experimental en este tramo con rango de
caudales de 0.28 – 0.1 Lts/s.
● Con la válvula 1 se variará el caudal, para así tomar varias lecturas.
● Se debe mantener abierta por lo menos una de las válvulas bola para evitar el
colapso del sistema.
● Al terminar los procedimientos se recomienda tener la precaución de abrir
todas las válvulas bolas, para así evitar pérdida de líquido en los manómetros
de U.
24
Figura 5.3.3 Banco de pruebas Uniminuto – Bogotá.
Se debe cambiar el agua del depósito de almacenamiento cada vez que se utilice el
equipo, quitando el tapón en la base del tanque y descargando directamente hacia
un sifón y volviendo a llenar el tanque por medio de una manguera con agua
completamente limpia. El tanque debe llenarse como mínimo 3 ⁄ 4 partes del mismo.
Se debe revisar periódicamente la instalación para detectar y corregir cualquier
escape que exista en los sitios de conexión, válvulas y los manómetros. Verificar el
ruido del motor, este debe permanecer similar constantemente, de lo contrario se
debe desmontar el motor y revisar los rodamientos. Se debe verificar durante cada
práctica los sellos de los tapones en el tanque de reserva, este sello debe ser
completamente hermético y sin fugas.
25
TABLA DE DATOS PRÁCTICA - PÉRDIDAS POR FRICCIÓN CORPORACIÓN UNIVERSITARIA MINUTO DE DIOS - SEDE GIRARDOT FACULTAD DE INGENIERÍA - PROGRAMA
DE INGENIERÍA CIVIL
TABLA DE DATOS PRÁCTICA - PÉRDIDAS POR FRICCIÓN CORPORACIÓN UNIVERSITARIA MINUTO DE DIOS - SEDE GIRARDOT FACULTAD DE INGENIERÍA - PROGRAMA
DE INGENIERÍA CIVIL
26
TABLA DE DATOS PRÁCTICA - PÉRDIDAS POR FRICCIÓN CORPORACIÓN UNIVERSITARIA MINUTO DE DIOS - SEDE GIRARDOT FACULTAD DE INGENIERÍA - PROGRAMA
DE INGENIERÍA CIVIL
TABLA DE DATOS PRÁCTICA - PÉRDIDAS POR FRICCIÓN CORPORACIÓN UNIVERSITARIA MINUTO DE DIOS - SEDE GIRARDOT FACULTAD DE INGENIERÍA - PROGRAMA
DE INGENIERÍA CIVIL
27
TABLA DE DATOS PRÁCTICA - PÉRDIDAS POR FRICCIÓN CORPORACIÓN UNIVERSITARIA MINUTO DE DIOS - SEDE GIRARDOT FACULTAD DE INGENIERÍA - PROGRAMA
DE INGENIERÍA CIVIL
TABLA DE DATOS PRÁCTICA - PÉRDIDAS POR FRICCIÓN CORPORACIÓN UNIVERSITARIA MINUTO DE DIOS - SEDE GIRARDOT FACULTAD DE INGENIERÍA - PROGRAMA
DE INGENIERÍA CIVIL
28
TABLA DE DATOS PRÁCTICA - PÉRDIDAS POR FRICCIÓN CORPORACIÓN UNIVERSITARIA MINUTO DE DIOS - SEDE GIRARDOT FACULTAD DE INGENIERÍA - PROGRAMA
DE INGENIERÍA CIVIL
TABLA DE DATOS PRÁCTICA - PÉRDIDAS POR FRICCIÓN CORPORACIÓN UNIVERSITARIA MINUTO DE DIOS - SEDE GIRARDOT FACULTAD DE INGENIERÍA - PROGRAMA
DE INGENIERÍA CIVIL
29
5.1.6. CÁLCULOS Y RESULTADOS DEL DISEÑO:
ρ*𝑉*𝐿
CÁLCULO DE Re:𝑅𝑒 = µ
1 𝑘/𝑑 2,54
CALCULO FACTOR DE FRICCIÓN: =− 2𝑙𝑜𝑔( 3,7
+
𝑓 𝑅𝑒 𝑓
2
𝐿 𝑉
CÁLCULO PÉRDIDAS POR FRICCIÓN:ℎ𝑓 = 𝑓 * 𝐷
* 2𝑔
𝑃1−𝑃2
∆ℎ𝑓 = φ*𝑔
2
𝑉
CÁLCULO DE Hm:∆ℎ𝑚 = 𝐾 2𝑔
CONDICIONES INICIALES
TEMP (°C ) : AMBIENTE
MATERIAL: PVC
CÁLCULOS INICIALES
DIÁMETRO 0.75pg
v=1,14E-0,6m/s
CALCULO N REYNOLDS Y RUGOSIDAD
Rugosida absoluta=1,50E-06
RELATIVa
CONDICIONES INICIALES
30
CÁLCULOS INICIALES
DIÁMETRO 0.5pg
v=1,14E-06 m/s
CALCULO N REYNOLDS Y RUGOSIDAD
Re=4,40E+04
RELATIVa
Rugosidad absoluta=6,00E-05
CONDICIONES INICIALES
CÁLCULOS INICIALES
DIÁMETRO 0.5pg
v=1,14E-06 m/s
CALCULO N REYNOLDS Y RUGOSIDAD
Re=4,40E+04
RELATIVa
Rugosidad absoluta=6,00E-05
CONDICIONES INICIALES
CÁLCULOS INICIALES
31
CAUDAL 0,5 L/s
DIÁMETRO 0.5pg
v=1,14E-06 m/s
CALCULO N REYNOLDS Y RUGOSIDAD
Re=4,40E+04
RELATIVa
Rugosidad absoluta=1,18E-04
CONDICIONES INICIALES
CÁLCULOS INICIALES
DIÁMETRO 1,25 pg
v=1,14E-06 m/s
CALCULO N REYNOLDS Y RUGOSIDAD
Re=1,76E+04
RELATIVA
Rugosidad absoluta=4,72E-05
∆ℎ𝑓 = 0, 54 𝑚𝑐𝑎
CANT. ACCES, 6
32
CONDICIONES INICIALES
CÁLCULOS INICIALES
DIÁMETRO 0,75 pg
CONDICIONES INICIALES
CÁLCULOS INICIALES
DIÁMETRO 0,75 pg
CÁLCULO DE Hm K=10,00
ℎ𝑚 =1,57 mca
TIPO DE ACCESORIO:CODO
DIÁMETRO (PULG.) : 0.75
CANT. DE ACCESORIO : 7
CÁLCULOS INICIALES
DIÁMETRO 0,75 pg
ℎ𝑚 =0,22 mca
33
TIPO DE ACCESORIO:SEMICODOS
DIÁMETRO (PULG.) : 0.75
CANT. DE ACCESORIO : 2,00
CÁLCULOS INICIALES
DIÁMETRO 0,75 pg
ℎ𝑚 =0,06 mca
SUMATORIA DE PÉRDIDAS
H accesorio succión 0m
HDT 11,52 m
34
5.1.7. CONCLUSIONES:
Las pérdidas se producen, básicamente, por las perturbaciones que los elementos
mencionados inducen en el flujo a presión, normalmente desarrollado, y que
ocasionan la aparición de turbulencias, torbellinos, vórtices y desprendimientos,
haciendo que parte de la energía del fluido se disipe en forma de calor.
35
Estas pérdidas, a pesar de llamarse secundarías, pueden ser más importantes que
las primarias (pérdidas por fricción), si la conducción es relativamente corta.
Generalmente se considera que la tubería es larga cuando es mayor a 1000
diámetros.
36
5.2.2.2.1. PÉRDIDA POR ENTRADA
A la entrada de las tuberías se produce una pérdida por el efecto de contracción que
sufre la vena líquida y la formación de zonas de separación; En este caso el
coeficiente K varía en mayor o menor medida como brusco sea el ingreso del fluido
a la tubería. Entradas con cantos redondeados producirán pérdidas inferiores a las
ocasionadas por entradas con cantos filosos o angulosos.
Con objeto de impedir la entrada de cuerpos sólidos a las tuberías, suelen utilizarse
estructuras de rejillas formadas por un sistema de barras o soleras verticales,
regularmente espaciadas, que se apoyan sobre miembros estructurales; dichas
rejillas obstaculizan el flujo y producen una pérdida de energía.
Se puede tener una aproximación media del coeficiente K con la fórmula de
Creager:
37
2
𝐾 = 1. 45 − 0. 45 ( ) ( )
𝐴𝑛
𝐴𝑏
−
𝐴𝑛
𝐴𝑏
(5.2.2)
Donde área neta de paro entre rejillas; área bruta de la estructura de rejillas.
38
El coeficiente K depende de la brusquedad de la ampliación y para encontrarlo se
usa la ecuación conocida como fórmula de Borda-Carnot, deducida a partir de las
tres ecuaciones fundamentales de la hidráulica:
2 2
ℎ𝐿 = ( 𝐴2
𝐴1 )
− 1
𝑉2
2𝑔
(5.2.3)
θ4a5 7 10 15 20 25 30 35 40 45 60
39
75 80
k 0.060 0.005 0.16 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.32
0.32 0.34
Tabla 5.2.1 Coeficiente de pérdida por reducción gradual del ángulo θ, según
Kisieliev
𝑔𝐷
𝑡𝑎𝑛 θ = 𝑉
(5.2.4)
40
Si el cambio de dirección es gradual con un curva circular de radio medio R y
rugosidad absoluta ε, para obtener el coeficiente de pérdida K se usa la gráfica de
Hoffman de la Figura 5.2.5 que además toma en cuenta la fricción en la curva
donde:
θ°
𝐾 = 𝐶𝑐 90 °
(5.2.5)
En el caso de válvulas, el valor del coeficiente K depende del tipo de válvula que se
tenga, del diseño particular de cada tipo y del grado de apertura de la misma, por
esta razón el coeficiente de pérdida debe ser proporcionado por los fabricantes. A
falta de estos datos se pueden utilizar los valores medidos que a continuación se
indican.
41
Tabla 5.2.2 Coeficientes de pérdida para válvulas de compuerta de diámetro
D= 50 mm
La unión cruz se puede dividir en unión igual de 4 vías y unión cruz reductora. La
unión igual de 4 vías se utiliza comúnmente para conectar cuatro tubos en un solo
lugar. Su diámetro de tubería principal es igual al diámetro de la tubería de
bifurcación. Los extremos del tubo adaptador en la unión igual de 4 vías son del
mismo tamaño. Con respecto a la unión cruz reductora, los extremos son de
diferentes tamaños, por lo que el tubo adaptador en el tubo de bifurcación es menor
que en la tubería principal.
42
densidad o el teflón. En la Figura 5.2.7 pueden observarse cruces de diferentes
materiales.
Las cruces son ampliamente utilizadas en muchas industrias, tales como la industria
del petróleo, la generación de electricidad, gas natural, productos químicos,
construcción naval, fabricación de papel y metalurgia. También son empleadas
actualmente en gran variedad de redes de distribución de agua a presión, como por
ejemplo sistemas de rociadores contra incendios, sistemas de agua potable y
sistemas de riego.
Las redes contra incendio son muy comunes en edificios de oficinas, almacenes
departamentales, fábricas, teatros y en general, construcciones con grandes
superficies. Estas tienen la particularidad de estar compuestas por tuberías con
diámetros pequeños y generalmente son de acero galvanizado.
Los coeficientes de pérdida local K están bien definidos para la mayoría de los
accesorios de uso común en tuberías. Sin embargo, para los cruces de tuberías es
poca la información sobre la pérdida de energía que se encuentra disponible en la
literatura y en los manuales especializados.
43
Un aspecto hidráulicamente problemático de analizar una cruz es el número de
posibilidades para la dirección del flujo y distribución. Hay cuatro escenarios de flujo
básicos que pueden ocurrir en una cruz:
44
% en cálculos cuando se utiliza un coeficiente de pérdida de energía constante para
una Tee en lugar de los coeficientes correctos para cada tramo.
Con la realización de este estudio, la pérdida de información energética correcta
para cada tramo de una cruz estará disponible para ayudar en los cálculos de
diseño y análisis de redes de tuberías de cruces (Sharp, 2009).
5.2.3.1. MATERIALES
5.2.3.2. METODOLOGÍA
45
volumen), el agua sube hasta dos tinacos de 400 litros de capacidad cada uno,
ubicados a una altura de 3.40m sobre el nivel del piso del laboratorio. Una vez que
los tinacos estuvieran llenos, el agua desciende y circula por las tuberías de tal
forma que al salir de ellas, pueda ser recirculada al cárcamo de bombeo. De esta
forma, no se requerirían fuentes continuas de agua, y se descartaron posibles
desperdicios del líquido.
4). Las válvulas en entradas y salidas, eran un requisito obvio. Se usaron con el fin
de permitir la variación de los gastos y poder obtener múltiples puntos en las
gráficas resultantes del análisis. En la Figura 5.2.11 se puede observar un esquema
del dispositivo experimental.
46
Figura 5.2.9 Sensores de presión conectados en la tubería.
Figura 5.2.10 Modulador para el registro de datos
47
5). Los medidores de flujo utilizados, se ubican en las entradas y salidas de los
cruces. Estos funcionan mediante una turbina alineada con el flujo, y muestran los
gastos instantáneos en una pantalla digital ubicada en su parte superior, con una
incertidumbre de ±1.0 l/h.
6). Para garantizar una medición sin perturbaciones, los medidores se ubicaron lejos
de cualquier tipo de accesorio, válvula o cambio de dirección, a una distancia
aproximada de diez veces el diámetro de la tubería conectada a los mismos. Los
cuatro medidores empleados tenían conexiones en sus extremos de 1”, por lo que
se conectan a tubos de 1” de diámetro nominal (29.80 ± 0.01 mm reales) y de 30.0
cm de longitud libre (ver Figura 5.2.11).
1). Tipo de flujo a estudiar: Se definió que se trabajaría con el flujo perpendicular.
Esto es alimentación del cruce por dos tramos adyacentes.
48
2). Diámetros a considerar: Para aprovechar al máximo la carga suministrada por
los tanques, se eligieron diámetros pequeños, empezando por el de D = 13 mm y D
= 19 mm.
5). Número de lecturas de presión por prueba. Con el propósito de evitar cambios
importantes en los niveles de agua de los tinacos, que llevaran a incumplir la
hipótesis de régimen permanente. Se estimó que era apropiado realizar pruebas de
4 minutos de duración, para obtener un total de 240 datos por prueba.
49
Tabla 5.2.3 Combinación de gastos en los tramos para flujo perpendicular
5.2.4.1. NOMENCLATURA
5.2.4.2. ESTIMACIÓN DE K
Las expresiones para estimar los valores de los coeficientes de pérdidas en cruces
(K), a partir de los datos medidos en el laboratorio, pueden obtenerse mediante un
balance de energía entre dos puntos cualesquiera del cruce, siempre que estén en
tramos diferentes de éste.
50
5.2.4.2.1. FLUJO PERPENDICULAR
En este caso existen dos flujos de entrada y dos de salida. Se asumió que las
pérdidas locales se concentran en los tramos de salida 3 y 4, por lo que existirán
dos coeficientes de pérdidas menores: K3 y K4. Los coeficientes pueden deducirse
haciendo un balance de energía entre los puntos de entrada 1 y 2 y los puntos de
salida 3 y 4.
Haciendo balance de energía entre los puntos de los tramos 1 y 3 (ver Figura
5.2.12), se tiene que:
Donde:
E1 representa el nivel de energía en el punto i;
hfi son las pérdidas por fricción en el tramo comprendido entre el punto i y el centro
del cruce;
hL3 las pérdidas menores correspondientes al tramo que une al punto 3 con el
centro del cruce.
El subíndice 3←1, indica que se trata del coeficiente de pérdidas del tramo 3,
calculado a partir del balance con el tramo 1. Esta forma de marcar los coeficientes
de pérdidas, es muy común para accesorios que implican algún tipo de unión o
separación de varios flujos, como es el caso de tees, yees, bifurcaciones y cruces,
entre otros.
51
Para el caso de cruces alimentados por dos tramos, esto podría llevar a confusiones
a la hora de decidir cuál sería el tramo de partida. Algo más adecuado consiste en
calcular los coeficientes de pérdidas menores desde todos los tramos de entrada,
promediarlos y reportarlos como K del tramo de salida estudiado. En este sentido
está orientado el presente estudio.
De donde:
(𝐸2− 𝐸3)− ℎ𝑓2− ℎ𝑓3
𝐾3←2 = 2 (5.2.8)
𝑉3/ 2𝑔
𝐾3←1+ 𝐾3←2
𝐾3 = 2
(5.2.9)
52
De tal forma que se obtiene:
𝐾4←1+ 𝐾4←2
𝐾4 = 2
(5.2.12)
53
Figura 5.2.14 Gráfica de Re4 vs K4 para el cruce de D = 13 mm
54
Figura 5.2.16 Gráfica de Re4 vs K4 para el cruce de D = 19 mm
Aun cuando la relación individual de los Rei con los Ki, fuera independiente de las
posibles combinaciones de entrada, se encontró que al hacer una relación entre un
Re de salida, con uno de entrada, y graficarla contra los coeficientes de pérdidas, la
nube de puntos se acomodaba mejor que en el primer caso, para los dos diámetros
estudiados (13 mm y 19 mm). En las Figuras 5.2.17 a 5.2.20, se reportan estas
gráficas, en las que se ha trazado una línea indicando la tendencia aproximada de
los puntos. En ellas se relaciona K3 contra la relación Re3/Re1, y K4 contra
Re4/Re2. Nótese que en ambos cocientes, se relaciona el Re de cada entrada con
el de la salida perpendicular a ella.
55
Figura 5.2.17 Gráfica de Re3/Re1 vs K3 para el cruce de D = 13 mm
56
Figura 5.2.20 Gráfica de Re4/Re2 vs K4 para el cruce de D = 19 mm
𝑎
𝐾𝐴𝑗𝑢𝑠𝑡 = 𝑐 + 𝑑 (5.2.13)
( 𝑅𝑒𝑠𝑎𝑙
𝑅𝑒𝑒𝑛𝑡𝑟 )
+𝑏
57
en las Figuras 5.2.17 a 5.2.20, correspondientes a los diámetros de 13 mm y 19
mm. Para cada caso se presenta también el valor del coeficiente de determinación.
Tabla 5.2.4 Ecuaciones de ajuste para los coeficientes de pérdidas K para
13mm y 19 mm
Los valores obtenidos de R2 son muy cercanos a la unidad, esto reafirma que existe
una tendencia clara de los puntos, y que los valores de los coeficientes de pérdidas
locales Ki, son bien explicados por las relaciones (Resalida/Reentrada).
5.2.5. CONCLUSIONES
58
● Cada cruce estaba compuesto por cuatro tramos perpendiculares de igual
diámetro. En ningún momento se emplearon cruces con diámetros
combinados.
● El dispositivo experimental tenía la posibilidad de variar los gastos de entrada
y de salida de los cruces, de tal forma que los flujos generados estuvieron
dentro de un intervalo del número de Reynolds comprendido entre 4,000 y
40,000.
● Las pérdidas por fricción en la vecindad del cruce no fueron despreciadas, y
se estimaron mediante la ecuación de Darcy-Weisbach.
Adicionalmente a estas aclaraciones, es bueno mencionar que los coeficientes de
pérdidas locales, K, calculados en la forma indicada, y reportados en forma de
ecuaciones, partieron de la suposición de un coeficiente K por cada tramo de salida
de los cruces. Entonces, para el caso de alimentación doble o flujo perpendicular se
estimaron dos coeficientes, K3 y K4.
El enfoque adoptado (una K por cada salida) y que puede observarse en las gráficas
de dispersión presentadas. Pese a que las primeras gráficas, que relacionan K con
su respectivo Re, parecen indicar una tendencia, que disminuye el valor del
coeficiente a medida que aumenta Re,. Si bien, los puntos de la gráfica indican
dicha tendencia, presentan una dispersión que no puede desestimarse; además, el
máximo número de Reynolds considerado, no fue superior a 40,000, por lo que
predecir el comportamiento más allá de este número sería un poco especulativo.
59
entre salida y entrada, para cada diámetro y para cada uno de los coeficientes K3 y
K4.
● Calcular las pérdidas por fricción como si se tratase de agua limpia para
después evaluar estas pérdidas en conducciones que transportan FCCSS
mediante el factor de amplificación.
● Proponer una metodología para calcular las pérdidas por fricción en
conducciones que transportan un fluido con cierto contenido de sólidos,
analizando tanto el método del tipo reológico como otros alternos como ser
tablas y gráficas para posterior elaboración del procedimiento ingenieril.
● Desarrollar mediante la metodología propuesta, expresiones numéricas de
fácil evaluación y comprobarlas mediante datos de campo.
El agua, el aceite y la mayoría de los fluidos son newtonianos, lo cual quiere decir
que en condiciones de flujo laminar la pérdida de carga es proporcional a la
velocidad y a la viscosidad cinemática (ν). Conforme la velocidad aumenta y
sobrepasa cierto valor crítico, el flujo pasa a ser turbulento. Los FCCSS tienen un
comportamiento similar al del agua clara, sin embargo, en algunas ocasiones se
comportan como si fuesen no newtonianos, es decir que, para éstos, el esfuerzo
cortante viscoso no varía conforme la ley de Newton prevista. Para los fluidos no
newtonianos existen diferentes procedimientos para la estimación del coeficiente de
pérdida por fricción. En este trabajo se analiza tanto el método basado en
60
consideraciones de tipo reológico, como otro alterno, presentado en Sanks. (diseño
de estación de bombeo) utilizado para elaborar el método denominado Ingenieril.
● el esfuerzo de cedencia, Sy
● el coeficiente de rigidez, η.
61
En las Figuras 5.3.1 y 5.3.2 se muestran los valores típicos de ambos coeficientes
los cuales presentan amplios intervalos de variación (las líneas discontinuas indican
los límites superiores e inferiores), introduciendo errores asociados a la
interpretación misma de estas gráficas.
62
Figura 5.3.2. Coeficiente de rigidez respecto al porcentaje de sólidos totales
(1)
(2)
Donde:
Vale la pena aclarar que las Figuras 5.4.1 y 5.4.2, así como (1) y (2), solamente son
válidas para flujos de agua con concentraciones de sólidos totales entre 0.1 y 10%,
dado que provienen de datos experimentales obtenidos únicamente para este
rango.
(3)
Donde:
63
● D es el diámetro de la tubería en m.
● η es el coeficiente de rigidez en kg/ (m.s).
(4)
(5)
Donde:
64
interior de la tubería, se compara la pérdida de carga obtenida con (5) y la pérdida
de carga para agua, y se toma la mayor de ellas. Así se determina si el fenómeno
predominante en la disipación de energía es un efecto de superficie o de cuerpo.
(6)
(7)
Donde:
65
diferencias de presión constantes aquí apareció una notable anomalía, para flujo
turbulento: cuanto más polímero se agregaba al solvente - dentro de ciertos límites y
con concentraciones no excesivamente altas- más aumentaba el caudal. Lo más
asombroso fue comprobar que existía una región de concentraciones en la cual,
bajo condiciones turbulentas, la solución polimérica parecía experimentar en su flujo
menos resistencia que la que la del solvente puro.
5.3.3 MÉTODO
66
Tabla 5.3.1. Factores de ampliación de pérdida, obtenidos mediante
el enfoque ingenieril
(8-a)
Donde k es el magnificador de las pérdidas en FCCSS y es una constante
definida como:
(8-b)
(9-a)
67
a=5.30879620
b=-6.8728801
c=-0.385149
El
Acueducto Chapala-Guadalajara (ACH-G) fue diseñado para entregar un caudal de
7.5 m3 /s a la ciudad de Guadalajara, siendo su fuente de abastecimiento el lago de
Chapala. El acueducto cuenta con una planta de bombeo (PB) donde se tienen
68
instaladas seis bombas de eje vertical, de las cuales cinco operan continuamente y
una permanece como de respaldo, cada una de éstas con capacidad de bombeo de
1.5 m3 /s para una carga dinámica de 138 m.
Tal como se puede observar en la Tabla 5.3.2, las mediciones de campo se han
realizado de manera periódica durante la presente década: Los aforos se realizaron
mediante la utilización de tubos de Pitot localizados en dos secciones, la primera en
la descarga de los equipos de bombeo y la otra ubicada a unos 700 m aguas arriba
del tanque de entrega. De esta manera se tenía certeza de los caudales transitados
por los tramos a presión y a gravedad. Adicionalmente se verificaban los caudales
mediante las curvas de las bombas, previa medición con manómetros, de la presión
de succión y de descarga de los equipos de bombeo en operación.
69
obstáculos en el interior de la conducción, hecho que se verifica numéricamente a
través de los altos coeficientes de correlación obtenidos, luego de ajustes numéricos
realizados a las mediciones. En la tabla también se muestra los resultados
obtenidos, en ésta, k representa la relación entre los coeficientes globales de
pérdidas de energía (𝒇𝑮) medidos, respecto del coeficiente teórico del acueducto
(correspondiente al medido en el año 1990, fecha en que comienza a operar la
conducción).
70
simultáneamente, antes de depositarse, produce alteraciones en la
concentración de sólidos en la interfaz pared-fluido.
Figura 5.3.5: material que migra y se deposita en las paredes interiores del ACH-G,
México
Figura 5.3.6. Detalle del material sólido (orgánico e inorgánico) que altera la interfaz
líquido- sólido en la conducción del ACH-G, México
71
5.3.5 CONCLUSIONES
72
5.4. PÉRDIDAS DE ENERGÍA DEBIDAS A LA FRICCIÓN
A medida que un fluido fluye por un conducto, tubo o algún otro dispositivo, ocurren
pérdidas de energía debido a la fricción interna del fluido. La ecuación general de la
energía indica que tales pérdidas de energía traen como resultado una disminución
de la presión entre dos puntos del sistema de flujo.
2 2
𝑝1 𝑣1 𝑝2 𝑣2
γ
+ 𝑧1 + α1 2
+ ℎ𝐴 − ℎ𝑙 = γ
+ 𝑧2 + α2 2
(5. 4. 1)
2
𝑓·𝐿·𝑉
ℎ= 2𝑔𝐷
(5. 4. 2)
Donde :
● D : Diámetro de la tubería
● g : Gravedad
73
● L : Longitud De la tubería
● v : Velocidad de fluido
● f : Coeficiente de fricción
Ésta ecuación nos permite calcular las pérdidas de energía tanto en flujos laminares
y turbulentos
32µ𝐿𝑉
ℎ𝑓 = 2 (5. 4. 3)
γ𝐷
2
𝐿 𝑉 32µ𝐿𝑉
𝑓. 𝐷 2𝑔
= 2 (5. 4. 4)
γ𝐷
64µ 64
𝑓= ρ𝑉𝐷
⇒𝑓 = 𝑅𝑒
64
𝑓= 𝑅𝑒
(5. 4. 5)
74
5.4.2.3 PÉRDIDAS DE ENERGÍA EN FLUJO TURBULENTO
0.316
𝑓= 1/4 (5. 4. 6)
𝑅𝑒
75
Si el conducto es rugoso y el flujo es turbulento (2000 < Re < 100000), se puede
aplicar la ecuación de Colebrook-White
1
𝑓
=− 2𝑙𝑜𝑔 ( 𝑘/𝐷
3.7
+
2.51
𝑅𝑒 𝑓 ) (5. 4. 7)
1
𝑓
= 2𝑙𝑜𝑔 ( ) + 1. 74
𝐷
2𝑘
(5. 4. 8)
representa mediante una familia de curvas, una para cada valor del parámetro
{\displaystyle k/D}, donde k es el valor de la rugosidad absoluta, es decir la longitud
(habitualmente en milímetros) de la rugosidad directamente medible en la tubería.
76
Figura 5.1.1 diagrama de moody
En donde:
𝑥1=7FA−FB
𝑋2=0.128−17 FA+2.5FB
𝑋3 =−0.128+13FA−2FB
R=NR/2000
77
FA= (𝑌3) −2
ϵ 5.7
Y= 3.7𝐷
+ 0.9
𝑁𝑅
Y=−2Log ( ϵ
3.7𝐷
+
5.7
0.9
4000 )
En las cuales NR es el número de Reynolds, ε rugosidad de la tubería, D diámetro
5.4.2.5.FÓRMULA DE HAZEN-WILLIAMS
Esta fórmula, establecida tanto para tuberías como para canales abiertos, ha sido
muy usada para tuberías. La selección de los exponentes se hizo con la idea de
conseguir una variación mínima del coeficiente para todos los conductos del mismo
grado de rugosidad.
0.63 0.63
𝑉 = 1. 1318𝐶𝑅ℎ 𝑆𝐹 (5. 4. 8)
Donde:
78
Si sf es igual a hf/L de la ecuación ( 5. 4. 9) se obtiene hf
1.852
10.67𝐿𝑄
ℎ𝑓 = 1.852 4.87. (5. 4. 8)
𝐶 𝐷
79
lenguaje de programación que incorporaba era el BASIC. Más que potente y más
que suficiente para la mayoría de estudiantes, científicos e ingenieros. En el manual
de la casio explicaba lo siguiente: “BASIC es mucho más fácil de usar que otros
lenguajes de programación tales como FORTRAN, lo que lo hace adecuado aún
para los novatos”.
Con Casio BASIC nos referimos específicamente al lenguaje de programación que
nació con la FX-4000P, y se incorporó posteriormente a las calculadoras gráficas de
Casio, cuya primera máquina es la FX-7000G. Sigue vivo en la actualidad, ya que
se ha mantenido durante todos estos años en las diferentes versiones que han ido
saliendo al mercado en todas las calculadoras gráficas.
Hay que decir que una de las principales características de este lenguaje es su
excesiva austeridad en comandos disponibles. Nace con un corto número de
comandos, y este problema no se soluciona en las diferentes versiones de
máquinas, cada vez más avanzadas que van surgiendo a lo largo de los años. Tan
limitado es el número de comandos, que hace inviable la programación de textos,
caracteres o sprites animados, la adecuado presentación de varios datos en la
pantalla, o incluso el poder presentar datos de manera secuencial sin que se
detenga el programa. A modo simple de ejemplo, durante muchos años, un
programa que contará números del 1 al 100 automáticamente en pantalla no era
posible.
Esto no quiere decir que fuera un lenguaje inútil, porque sí permitía, por ejemplo, la
ejecución de complejos cálculos matemáticos que permitían resolver problemas
científicos, de ingeniería o de finanzas. Pero a los ojos de los programadores puros,
quedaba bastante limitado. En la tabla de abajo se puede observar el mismo
programa, el que usamos para ver el rendimiento de las calculadoras, en tres
etapas. El código de la izquierda corresponde a la calculadora FX-4000P, cuando se
inicia el Casio BASIC y tiene la extraña característica de tener que introducir todo el
programa en una sóla línea. El código central puede pertenecer a cualquier
calculadora gráfica desde el año 1985, hasta la aparición de la 2ª generación de
este lenguaje 11 años más tarde. Finalmente, el código de la derecha, es válido
para las calculadoras que ya implementan el Casio Basic de 2ª generación en 1996
con la FX-9750G. Puede observarse que el código es básicamente el mismo con
80
mínimas modificaciones. Esto es debido a que no se emplean las características
gráficas ni de lectura del teclado disponibles a partir de 1996.
81
Desarrollo de la aplicación.- Para desarrollar una aplicación o una subrutina primero
se debe tener un claro conocimiento de que valores se van a utilizar, cuales van a
ser las entradas y cuales las salidas, para así definir cómo se van a entrar y sacar
estos valores. De esta manera se realiza un programa que consta básicamente de
las siguientes partes fundamentales:
● Configuración del sistema de adquisición de datos.
● Acople del sensor: En condiciones de estado estacionario sin la presencia de
fugas y/o aire, se ajustan los voltajes correspondientes a los flujos de entrada
hacia el Field point.
● Análisis de datos.
● Implementación de la aplicación
82
muy esquemáticos y de gran elegancia. Es un programa de mucho poder
donde se cuentan con librerías especializadas para manejos de FieldPoint,
Redes, Comunicaciones, Análisis Estadístico, Comunicación con Bases de
Datos (útil para una automatización de una empresa a nivel total). Como se
programa creando subrutinas en módulos de bloques, se pueden usar otros
bloques creados anteriormente como aplicaciones por otras personas.
● Descripción del programa:
El programa para el cálculo de pérdida de carga en tubería del que está
construido el banco de pruebas se realizó en LabVIEW versión 7.0 Express.
El programa contiene los siguientes subprogramas:
INICIO.VI Programa de Control
EXPERIMENTOS.VI Ventana de Diálogo
PRINCIPAL.VI Programa para escoger Tubería y Presión
PANTALLA1.VI Programa de Control y datos del sensor 161
PERDIDA.VI Cálculo de Pérdidas y coeficientes
REPORTE1.VI Visualización de datos del experimento y curvas asociadas.
GRABAR.VI Ventana de Diálogo
REPORTE AUXILIAR.VI Subprograma para guardar datos. (No visible por el
usuario)
El FieldPoint, se programa en base a la plataforma “Measurement &
Automation” (MAX), que es el software indicado para trabajar con esta clase
de elementos. Dicha programación es importante para el buen desempeño
del sensor.
83
la presión de alimentación y la pérdida de carga en el espacio ínter tubular
(coraza).
(b) Repetir el cálculo anterior mediante un programa en LabVIEW
SOLUCIÓN:
(b) La Fig. obtenida mediante la programación en LabVIEW, muestra los
siguientes resultados: La presión de alimentación es igual a 719,83 kPa y la
pérdida de carga en el intercambiador es de 10095,5 m. La diferencia
se debe a que en el cálculo mediante labVIEW se ha despreciado las
pérdidas menores.
84
caudal?. (b3) ¿Si se modifica la presión inicial a 200 kPa manteniendo las
demás condiciones, cuál será el nuevo caudal?
SOLUCIÓN:
Los cálculos anteriores han sido programados en LabVIEW y el panel frontal
muestra el resultado final.
Como se observa el flujo volumétrico calculado es de 0,086 m3/s. (b) El panel frontal
de LabVIEW se presenta a continuación:
(b1) Modificación de la longitud de las tuberías:
85
(b3) Modificación de la presión inicial:
86
En resumen esta es la metodología que se recomienda para el cálculo de la pérdida
de energía debido a la fricción.
Figura 5.1.3. Diagrama de Flujo para el cálculo de pérdidas por fricción utilizando la ecuación
de Darcy -Weisbach
● Equipo Utilizado:
1. Banco de pruebas para Pérdidas Longitudinales en tubería
2. Bomba de ½ HP
3. Un termómetro
4. Recipiente recolector de agua para purga de tubería
87
5. Elementos de Automatización (Computador, FieldPoint, Sensor)
● Fórmulas Utilizadas:
1. Área
2.Caudal
3. Gradiente hidráulico
4.Numero de Reynolds
5. Velocidad
6. Cálculo de Pérdidas en Longitud (HG en 3⁄4’’ y HG en ½ ‘’) y Ecuación de
Darcy Weisbach
7. Cálculo del Coeficiente C (HG in ¾’’ , HG en ½’’ y PVC en ½’’)
8. Cálculo de Pérdidas en Longitud (PVC en ½ ‘’)
8.1 Veronesse-Datei
● Procedimiento
1. Revisar con un multímetro el valor de la fuente de alimentación del sensor
(12V) sin que se encuentre conectado el sensor.
2. Cebar la bomba de ½ HP.
3. Cerrar válvulas V3 , V6 y V7
4. Encender la bomba con válvulas V1 , V2 , V4 , V5 , V9 , V10 y V11 abiertas
subiendo la presión a 15 psi con válvula V8 durante un minuto para sacar el
aire en las tres tuberías.
5. Seleccionar tubería HG de ¾’’ cerrando llaves V4 , V5 , V9 y V10.
6. Abrir válvula V8 para bajar la presión a 8 psi mostrado en el manómetro
M1. 7. Conectar manguera plástica transparente a los acoples rápidos A1 de
los manómetros M2.
8. Elevar presión a 10 psi con V8.
9. Abrir válvula V3 para purgar y limpiar evitando la existencia de aire en la
manguera. Q L D hf m agua * * 0.00092 ( ) 1.8 4.8 = 158
10. Taponar la salida de manguera con los dedos y cerrar válvula V3. (Sólo
primera práctica).
11. Conectar la manguera a las dos entradas del sensor (al mismo tiempo).
12. Acoplar el conector del sensor a E1.
13. Abrir llave V3 al mismo tiempo para evitar que la diferencia de presión
sea mayor a 5 psi.
14. Encender la fuente del sensor.
88
15. Realizar la toma de datos en el programa variando la presión mostrada en
M1 desde (10 a 20) psi mediante V8 y comparar los datos del programa con
los manómetros M2.
16. Para cada práctica medir la temperatura del agua.
17. Una vez realizada todas las prácticas deseadas bajar la presión de ésta a
10 psi utilizando V8.
18. Cerrar válvula V3 al mismo tiempo
19. Apagar la fuente del sensor
20. Seleccionar tubería HG de ½’’ abriendo V4 , V9 y V10 cerrando llaves V2
, V5, V6 , V7 , V8 y V11
21. Con Válvula V9 bajar a mínima presión 12 psi mostrado en el manómetro
M1 (Repetir paso 7)
22. En caso de observar la presencia de aire desconectar la manguera del
sensor (al mismo tiempo) y repetir los pasos 9, 10 y 11 pero con válvula V6
23. Acoplar el sensor al conector E2 159
24. Abrir válvula V6 al mismo tiempo evitando presiones mayores a 5 psi
25. Repetir paso 14 y realizar la toma de datos en el programa variando la
presión mostrada en M1 desde (12 a 18) psi mediante V9 y comparar los
datos del programa con los manómetros M2.,
26. Repetir paso 16
27. Una vez realizada todas las prácticas deseadas bajar la presión de ésta a
12 psi utilizando V9.
28. Cerrar la válvula V6 al mismo tiempo
29. Apagar la fuente del sensor.
30. Seleccionar tubería PVC de ½’’ abriendo V5 cerrando llaves V1 , V4 y V7
31. Con Válvula V9 bajar a mínima presión 12 psi mostrado en el manómetro
M1 (Repetir paso 7)
32. En caso de observar la presencia de aire desconectar la manguera del
sensor (al mismo tiempo) y repetir los pasos 9, 10 y 11 pero con válvula V7
33. Acoplar el conector del sensor a E3
34. Abrir válvula V7 al mismo tiempo para evitar que exista una diferencia de
presión mayor a 5 psi
89
35. Repetir paso 14 y realizar la toma de datos en el programa variando la
presión mostrada en M1 desde (12 a 18) psi mediante V9 y comparar los
datos del programa con los manómetros M2.,
36. Repetir paso 16 160
37. Una vez realizada todas las prácticas deseadas bajar la presión de ésta a
12 psi utilizando V9.
38. Cerrar la válvula V7 al mismo tiempo
39. Apagar la fuente del sensor
40. Finalizada la práctica en las tres tuberías abrir válvulas V1 , V2 , V4 , V8 y
V11
41. Apagar la bomba
42. Desconectar los acoples rápidos A1 y guardar el sensor.
Descripción del programa.
El programa para el cálculo de pérdida de carga en tubería del que está
construido el banco de pruebas se realizó en LabVIEW versión 7.0 Express.
● Utilización del programa:
-Cargar el programa desde el escritorio. El nombre del programa es Tesis 1.4
90
- Luego de un par de segundos aparecerá una pantalla de diálogo que nos
indica la cantidad de experimentos permitidos.
91
- Al escoger cualquiera de las opciones:
● Hg ¾
● Hg ½
● P.V.C. ½,
- Aparecerá un mensaje que indica que la presión a trabajar en el banco de
pruebas debería ser mínimo de 10 P.S.I. y máximo de 25 P.S.I., esto se debe
a especificaciones de uso del sensor.
- Al indicar la presión y la tubería aparecerá un botón “Continuar”con el cual
se puede desplazar a la siguiente pantalla.
92
- En la siguiente pantalla, se debe presionar el botón “Iniciar”, para empezar
con la toma de datos. Se debe tener en cuenta el canal del FieldPoint en el
que se va a trabajar del combo “Canal del FieldPoint”, este valor, por
sugerencia y omisión se encuentra en el Channel 1.
93
- Se debe ingresar por teclado el número de muestras que tomará el sensor y
el tiempo de espera de una muestra hasta la siguiente, el tiempo
recomendado y mínimo es de 2 segundos. Al presionar el botón “OK”, se
iniciará la toma de datos, y luego del tiempo de espera se cargarán los arrays
con el voltaje que envía el sensor, al lado derecho se cargará un array con la
transformación a presión en metros de agua.
94
- Al presionar los botones “Calcular” el programa despliega todos los datos
realizados y los cálculos pertinentes.
- En esta ventana también se encuentran curvas con las que se han realizado
algunos de los cálculos, tales como, Presión vs Caudal y,
95
Viscosidad Cinemática en función de la Temperatura.
96
- El botón “Salir”, que como su nombre lo indica nos permite salir de la
aplicación.
- El botón “Nuevo”, con el cual todos los datos se ponen en 0 y se reinicia la
práctica desde el paso 3, finalmente;
- El botón “Grabar”, el mismo que envía una ventana de diálogo en la cual
debemos indicar el nombre con el que se grabará la práctica realizada.
97
- Si se elige la opción Salir del Sistema, saldrá al sistema operativo. El
reporte final se verá de la siguiente manera:
98
5.4.4. ANÁLISIS DE RESULTADOS
Figura 5.1.4 Listado del programa en BASIC para el cálculo de pérdidas por fricción
● Tubería HG ¾:
Temperatura para todos los casos 17ºC
99
● Tubería HG ½
Temperatura para todos los casos 17ºC
● Tubería P.V.C. ½
Temperatura para todos los casos 17ºC
100
5.4.5. CONCLUSIONES
101
● Diámetro interno de la tubería (Di)
● Longitud de la tubería (L) y
● Rugosidad relativa de la superficie interna de la pared de la tubería
102
Fig. 5.5.2 diagrama de Moody para determinar el número de Reynolds
Rugosidad relativa = ∊ / Di
Donde:
103
5.5.2.3 PROPIEDADES FÍSICAS DEL FLUIDO
Junto con las características de la tubería, las propiedades físicas del fluido
transportado a través de la tubería afectan el diseño del oleoducto. Hay seis
propiedades del líquido que deben ser consideradas:
Viscosidad
µ
V= ρ
donde:
● v= viscosidad cinetica
● µ = viscosidad absoluta o dinámica
● ρ = densidad
𝑚
ρ= 𝑣
donde:
● ρ = 𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑
● m= masa
● v = volumen
104
Grado api: la gravedad o grado API determinará la densidad de un
hidrocarburo. Entre grados de API se obtiene un hidrocarburo de más baja
densidad.
crudo API
liviano 32 – 42 ⁰ API
pesado 10 – 22 ⁰ API
Punto de fluidez se define como la temperatura más baja a la cual un fluido sigue
siendo líquido vertible, es decir sigue comportándose como fluido. Aunque hay
crudos que pueden ser bombeados a temperaturas por debajo de su fluidez, se
requiere más energía para hacer esto.
105
5.5.2.4 PARÁMETROS A CONSIDERAR EN EL DISEÑO DE OLEODUCTOS:
Estas pérdidas pueden ser mayores o menores, dentro de las pérdidas mayores se
tiene las provocadas por el cambio de elevación, también conocido como “cabeza
estática” y las pérdidas debido a la fricción existente entre el fluido en movimiento y
la pared interior del ducto.
106
5.5.2.6 PÉRDIDA DE PRESIÓN POR FRICCIÓN USANDO LA ECUACIÓN DE
DARCY
MAOP = ( 2 𝑥 𝑆𝑦 𝑥 𝑡
𝐷𝑒 )x F x J x T
donde:
107
5.5.3. MÉTODO Y MATERIALES
El ducto se llevará de manera aérea a lo largo de todo el trayecto entre los puntos
“A” y “B”.
Donde:
108
● Viscosidad cinemática del crudo = v
● Densidad 15⁰C = ρ
Nota 2: Ducto aéreo a lo largo del recorrido (no existen tramos enterrados)
Por lo general existen algunos datos conocidos, que son con los que se realiza las
primeras consideraciones de diseño, normalmente se conoce el sitio desde donde
(explotación y producción) y el sitio hasta donde se desea llevar el crudo
(exportación, almacenamiento, refinación, etc.).
Para determinar la ruta, se deben realizar estudios de algunas clases como son:
● Estudios geográficos.
● Estudios topográficos.
● Estudios geológicos.
● Estudios de impacto ambiental.
● Estudios y análisis de riesgos de toda clase.
● Métodos y sistemas de montaje e instalación de ductos.
109
seleccionada, ya que es necesario para el análisis integral de las presiones
involucradas, por lo que para este problema se puede considerar como dato
conocido.
Existen varios datos que se conocen antes de realizar el diseño del ducto por lo
general, estos parámetros son:
110
5.5.3.7 ALGORITMO DE DISEÑO
Basados en estos datos conocidos iniciales, primero se selecciona por criterio del
diseñador un tamaño de tubería con especificaciones y características tales como
grado, clase, etc.
111
5.5.3.8 CÁLCULO DE LA PRESIÓN MÁXIMA ADMISIBLE DE OPERACIÓN
Con el dato de Caudal (Q) que es constante y depende de las condiciones del
problema, se procede a calcular la velocidad de flujo, que también depende del
diámetro interior del ducto o tubería.
Q = Vel * A
112
A = π * Di2/4
● Donde:
● Vel = velocidad
● Q = caudal de flujo
● Di = diámetro interior
Clase de Flujo:
● Flujo Laminar
● Turbulento
113
f = 64 / Re
f = 0,028
Esta ecuación combina la ecuación de Bernoulli con los factores de pérdida por
fricción, el trabajo realizado por el fluido y la energía agregada al fluido.
Donde:
PA = PB + GPe + GPf
114
Donde:
PB = Presión mínima necesaria en la entrada estación “B” para llevar el crudo hacia
los tanques de almacenamiento y llenarlos, considerando que la estación B es una
estación de Almacenamiento y Refinación
de Bombeo “A”
PA = PB + GPe + GPf
115
En este caso MAOP es menor que PA, por lo que NO cumple la condición primordial
que es que MAOP tiene que ser mayor a PA.
La primera etapa de cálculo tiene como objetivo determinar los parámetros previos
para poder calcular el perfil de presiones.
116
De los valores que se obtienen de la diferencia de presiones entre MAOP y PA, se
debe identificar todos los positivos.En color lila se pueden apreciar los valores
positivos producto de la diferencia entre MAOP y PA.
Este grupo de valores positivos, son de los diferentes espesores existentes para
diámetro de 24 pulgadas.
Para este proyecto se analizarán con diámetros de: 26”, 24”, 22”, 20”, 18” 16”.
o Alternativa 1
117
Fig. 5.6 Resultados en Matlab con tubería de diámetro 16 plg
o Alternativa 2
118
o Alternativa 3
o Alternativa 4
119
Fig. 5.9 Resultados en Matlab con tubería de diámetro 24 pulg
De las cuatro primeras alternativas el número uno es la que más se aproxima a una
solución posible analizando se puede observar que es necesario usar espesores
mayores para lograr que el perfil de presiones sea el adecuado
120
Tabla 5.5.6 datos disponibles de la segunda etapa del cálculo
SECUENCIA DE CÁLCULO
121
Fig. 5.5.10 Resultados en Matlab con tubería de diámetro 14 plg
Diseño hidráulico de un oleoducto que transporta crudo medio desde un punto "C"
(estación de bombeo “C”) hasta un punto "B" (estación de almacenamiento y
refinación “B”) en condiciones de estado estable.
122
fig 5.5.12 Gráfico perfil de presiones requerido para el segundo tramo
Se puede tratar de buscar una alternativa que cumpla el bombeo directo desde el
punto “A” hasta el punto “B”, pero la presión requerida va a ser demasiado alta, por
lo que sería difícil conseguir equipos de bombeo de esa potencia, por tanto se debe
considerar realizar el cálculo completo en dos tramos.
Para el análisis de este tramo “C” - “B” se procede de la misma manera que el
utilizado para el recorrido total es decir el tramo “A” - “B”. Y finalmente seleccionar
una alternativa que complete el recorrido.
123
De acuerdo a los gráficos obtenidos con las alternativas seleccionadas, se analizan
las presiones obtenidas y las necesarias para el objetivo planteado llegando a la
conclusión que es necesario seleccionar las alternativas 3 o 4 para el primer tramo
de diseño, es decir el tramo entre la estación "A" y la nueva estación requerida a la
que se le llamará estación de bombeo “C”.
Para el primer tramo (“A” - “C”) = 70 Km, se puede utilizar ya sea la alternativa 3 o la
alternativa 4.
Al punto “C” se denominará estación de bombeo “C” y se encuentra entre “A” y “B”.
5.5.5. CONCLUSIONES
124
El diseño hidrodinámico o hidráulico de un oleoducto es una parte del diseño
Integral de un ducto, ya que se debe tomar en cuenta muchos temas y aspectos
como seguridad, economía, facilidad de instalación y montaje, operatividad,
mantenimiento, etc. para diseñar completamente un oleoducto.
5.6.2.1. Oleoducto
125
La reducción de arrastre se refiere a la disminución en la resistencia para que un
fluido pueda fluir; esto se lleva a cabo por medio de la adición de ciertos aditivos, los
cuales buscan interferir con la formación de la turbulencia en la tubería y así
disminuir las pérdidas de energía que se generan normalmente cuando se
transporta un fluido.
Los agentes reductores de arrastre son aditivos que se agregan a la tubería para
disminuir las pérdidas de presión que se presentan dentro de ella. Al agregar este
tipo de aditivos, se disminuye la fricción generada dentro de la tubería y de esta
forma la energía que se requiere para bombear el crudo es menor, reduciendo así
los costos de transporte del crudo.
126
El polímero de elección (poliacrilamida) muestra una excelente afinidad por
hidrocarburos de petróleo pesado y es capaz de alcanzar niveles elevados de
reducción de la resistencia aerodinámica en diversos escenarios operativos.
5.6.2.5. Viscosidad
127
5.6.2.6. Pérdidas por fricción
Las únicas pérdidas que se pueden presentar en una tubería son las pérdidas de
presión, ya que la energía potencial solo depende de la posición y la energía
cinética debe permanecer constante si la tubería mantiene sus propiedades. La
pérdida de presión también se puede determinar como la pérdida por altura por
fricción con la ecuación de Darcy Weisbach.
2
𝐿 𝑣
ℎ𝑓 = 𝑓 𝑑 2𝑔
5.6.2.7. Turbulencia
5.6.3.1. Metodología
128
diferentes escenarios propuestos, lo cual ayudará a obtener una solución óptima
para el problema cumplirá el objetivo del proyecto.
Para la evaluación del crudo en el Oleoducto, se utilizó crudo pesado extraído del
campo Rubiales, ubicado en cercanías al municipio de Puerto Gaitán en el
departamento del Meta, el cual tiene una densidad de 12.3° API (987,14 Kg/m3).
Para la evaluación del proceso de reducción de fricción se empleó cuatro diferentes
concentraciones elegidas por la compañía DELRIO S.A.S.
129
Posteriormente, se evaluó el efecto del DRA sobre el crudo mediante estas mismas
pruebas de flujo. Para esto, se hizo un barrido de ida de 0,01 s-1 hasta 10 s-1 y de
vuelta de 10 s-1 hasta 0,01 s-1 en un Reómetro Híbrido DHR-1 de TA instrumentos
utilizando una geometría de platos paralelos de 20 mm de diámetro, a diferentes
condiciones de temperatura en un rango entre 20 y 60°C para evaluar su posible
efecto sobre la propiedad reológica viscosidad.
130
5.6.3.7. Prueba de fluidez con Copa Ford (Vidrio)
Nota: La figura muestra el montaje realizado para la prueba en copa Ford con
tubería de vidrio.
131
Figura 5.6.3.8.1. Prueba de fluidez con Copa Ford (Acero)
Nota. La figura muestra el montaje realizado para la prueba en copa Ford con
tubería de acero.
Planteamiento de la simulación
Para diseñar una simulación analítica se debe elegir un modelo analítico el cual se
adecue más a las características del proyecto, para este caso se diseñará un
escenario donde se simule el comportamiento de las pérdidas de fricción a través
132
del transporte del crudo en el Oleoducto. Para este escenario, se asumió el modelo
Peng – Robinson debido a que se ajustan a los lineamientos del proyecto.
Adicionalmente para facilidad en el uso del software ASPEN HYSYS, se eligieron
los datos del Campo Bachaquero, localizado en Venezuela los cuales se encuentran
precargados en el simulador y presentan características muy similares al crudo
transportado actualmente en el Oleoducto de los Llanos.
133
Figura 5.6.3.9.2. Esquema inicial para la simulación en el Software Aspen HYSYS.
Una de las características del software Aspen HYSYS, es que permite hacer la
carga manual de los datos de un crudo en específico o elegir de una base datos
previamente cargada diferentes tipos de crudos con una amplia variedad de rasgos
y características. A través del software Aspen HYSYS, se realiza la búsqueda y
selección de un crudo que posea características similares a las que actualmente se
tienen en el Oleoducto de los Llanos, debido a que no se cuenta con información
suficiente para realizar la carga de los datos en el simulador.
Importe de datos
134
Figura 5.6.3.9.3. Importación de datos.
Selección de datos
135
Posteriormente en la Figura 5.6.3.9.4, se presentan las principales propiedades del
crudo muestra en la corriente de entrada. Es importante en este paso hacer una
verificación y si es necesario ajustar la viscosidad a los valores obtenidos de las
pruebas de flujo con DRA realizadas anteriormente. Para el caso presentado la
figura 22, se toma el valor de viscosidad 1550 cp @ 800 ppm (concentración actual)
obtenido de la prueba de flujo con DRA a 180°F.
136
Figura 5.6.3.9.5. Propiedades de la corriente de entrada
Accesorios
Accesorios disponibles
137
Figura 5.6.3.9.6. Variedad de accesorios disponibles y sus correspondientes
características.
Después de elegir los accesorios que se necesitan añadir, se ajustan los datos a las
mismas condiciones presentes en el Oleoducto de material, diámetro y longitud
como se muestra en la Figura, para una mayor precisión en los resultados de la
simulación. Cabe resaltar que en algunos accesorios el simulador no tenía valores
similares en algunos parámetros como rugosidad o el material no era el mismo
utilizado en el Oleoducto, por ende, se ingresaron los valores manualmente para
mantener la simulación fiel al escenario actual.
Por último, se procede a correr la simulación para la obtención de las pérdidas por
fricción en el oleoducto y se repite el proceso para cada concentración propuesta .
Simulación finalizada
138
Simulación a las diferentes concentraciones propuestas
139
Figura 5.6.4.2.1. Prueba de flujo No. 1 @ 20 °C – 60 °C.
140
En las Figuras 5.6.4.2.3. y 5.6.4.2.4. @ 20° C se presentan propiedades tixotrópicas.
La tixotropía es la propiedad que tienen algunos fluidos, a mostrar un cambio de su
viscosidad en el tiempo, a medida que se incrementa el esfuerzo de cizalla más se
reduce la viscosidad. Asimismo, tarda un tiempo finito en regresar a tener una
viscosidad en equilibrio, el cual se muestra en las gráficas cuando es sometida a la
ida y vuelta en la prueba de flujo.
141
presencia del esfuerzo prolongado al que se es sometido en la prueba puede
homogeneizar el petróleo. Se demuestra a su vez que la viscosidad es afectada
directamente por la temperatura, en una proporción inversa, por ende, se observa
una disminución uniforme a medida que las temperaturas varían de 20°C a 60°C.
142
Figura 5.6.4.3.2. Prueba de flujo No. 6 @ 180 °F
Para la figura 5.6.4.3.3, el valor se redujo a 1,32 Pa.s @ 2300 ppm con respecto a la
concentración inicial, esta reducción de alrededor del l5% en comparación a la línea
de tendencia de solo crudo (180ºF) se debe a que ya se cuenta con unos valores de
concentración bastante altos que permiten ver una reducción significativa y rentable
para el propósito del proyecto.
143
Figura 5.6.4.3.4. Prueba de flujo No. 8 @ 180 °F
Nota. Flujo másico de prueba de fluidez para Crudo 12º API @ 60°C
144
utilizadas con incrementos bastante considerables en el flujo, estos incrementos de
flujo son del doble o incluso el triple en el caso de la muestra a 2300 ppm, y esto se
debe a que la prueba es realizada en el mismo material utilizado en el oleoducto, lo
cual potencia una de las características principales del agente reductor.
Nota. Flujo másico de Prueba de fluidez para Crudo 12º API @ 60°C
Los resultados de la prueba de fluidez con Copa Ford, están afectados por el
material y dimensiones de los embudos y sus tubos empleados. Por ejemplo, el
diámetro del tubo del embudo de vidrio y metálico difieren significativamente (9 mm
o 0,354 in; y ½ in o 12,7 mm, respectivamente), la longitud de los tubos de vidrio y
metálico (90 mm o 3,54 in; y 500 mm o 19,6 in, respectivamente), además de
irregularidades como la rugosidad relativa, tendrán efecto directo en el desarrollo y
resultados del transporte del crudo pesado.
145
El último escenario propuesto @ 2300 ppm es la concentración con mayor grado de
efectividad en el incremento de flujo en el crudo pesado, con un aumento de hasta 3
veces si se compara con el valor del crudo sin agente reductor.
5.6.4.6. Simulación
146
Figura 5.6.4.6.2 Datos simulación Aspen HYSYS
Valores de pérdidas por fricción obtenidos (kPa) del simulador Aspen HYSYS
147
Figura 5.6.4.6.4 Comportamiento del número de Reynolds a cada concentración
propuesta.
5.6.5. CONCLUSIONES
148
No. 3 obtuvo un valor de 570847,32 kPa y por último el escenario No. 4 521587,01
kPa. Es notorio que el escenario más atractivo es el No. 4 a 2300 ppm.
6. BIBLIOGRAFÍA
APA 7th - Cuellar Muñoz, J. A. (2021) Optimización del desempeño energético para
el transporte de hidrocarburos en el oleoducto de los Llanos Orientales mediante la
utilización del agente reductor de fricción extreme power. [Trabajo de grado,
Fundación Universidad de América] Repositorio Institucional Lumieres.
https://hdl.handle.net/20.500.11839/8277
149
150