FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
"DISEÑO DE UN SISTEMA DE PROTECCIÓN HIDRÁULICO PARA UN POZO
PETROLERO"
PROGRAMACION ESTOCASTICA
PRESENTADO POR:
Yelkin Jose Jimenez Romero
Pedro Andrés Porras Ramírez
Carlos Andres Chaparro Corredor
Geraldin Tamayo
Monica Tatiana Villamil Quiroga
DOCENTE:
Martha Suarez rincón
INSTITUCION UNIVERSITARIA POLITECNICO GRANCOLOMBIANO
FACULTAD DE INGENIERIA INDUSTRIAL
AÑO 2025
1 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
TABLA DE CONTENIDO
INTRODUCCIÓN......................................................................................................................3
DEFINICION DE VARIABLES A MODELAR........................................................................4
PRUEBA DE BONDAD...........................................................................................................10
BIBLIOGRAFÍA.......................................................................................................................12
2 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
INTRODUCCIÓN
La búsqueda de la eficiencia operativa en la industria de los hidrocarburos lidia con
todo lo inherente a la incertidumbre de presentar fallas de los equipos, así como de los
sistemas, uno de los sistemas más críticos en este tipo de industria es el sistema de
protección hidráulico que hace parte de la seguridad para prevenir la contaminación por
derrames de crudo. La integración de los modelos estocásticos en la toma de decisiones ha
sido un acierto cuando se trata de controlar adecuadamente los sistemas para aprovechar al
máximo los recursos y mejorar la confiabilidad del funcionamiento de las operaciones. De
modo que, las Cadenas de Markov en Tiempo Continuo son un procedimiento confiable
para representar el funcionamiento de los sistemas.
El objetivo de este trabajo es el de elaborar un modelo probabilístico basado en las
Cadenas de Markov que represente la operación de un sistema de protección hidráulico que
tiene 10 válvulas de seguridad; a partir de las características de los datos de fallos con sus
correspondientes revisiones se espera generar una representación del sistema para
diferentes condiciones de operación. Igualmente, se generarán diferentes estrategias de
mejoras para ofrecer una mayor seguridad al sistema y lograr un ahorro en los costes de los
tiempos de inactividad; y finalmente un análisis económico en el que se calculará la
expectativa en términos de los ingresos diarios a partir de la operación del sistema.
3 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
DEFINICION DE VARIABLES A MODELAR
Utilizando en las cadenas de Markov donde se estudian las fallas y revisiones de las
válvulas que componen el sistema de un pozo petrolero. La finalidad es presentar la
estructura de este modelo, poniendo el acento en las variables de estado que se aplican a
esta especificidad de CM en un caso de estudio de 10 válvulas que componen un sistema de
protección hidráulico.
Variable de estado para el sistema representa el estado global del sistema hidráulico
en un momento o intervalo de tiempo t dentro del pozo petrolero, esta variable de estado
para cualquier válvula que componga el sistema
Sea Ck(t) el estado de la válvula k en el instante o periodo t. Así, el estado global
del sistema se puede representar como:
r (t)= {C1 (t), C2 (t), C3 (t), C4 (t),…, Ck (t)}
Espacios de estado donde se muestra un modelo de combinaciones parciales e
independientes aplicables a cada válvula, basado en la información previamente analizada:
S (Ck(t)) = espacio de estado para cada válvula: {m: la válvula k opera
correctamente; n: la válvula k no opera correctamente}
El espacio de estado del sistema en su conjunto puede representarse como:
S(r(t)) = {C1(t)(n,m),C2(t) (n,m),C3(t)(n,m),C4(t)(n,m),…,Ck(t)(n,m) }
4 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
Así k solamente puede variar hasta 10 por la composición del sistema de protección
hidráulico y los parámetros del estudio.
Este estudio se compone del análisis de un sistema hidráulico con 10 válvulas. Para
este desarrollo se han realizado pruebas estadísticas de bondad de ajuste en el
comportamiento de las válvulas, llegando a determinar que un 60% de las fallas de las
válvulas cobran un comportamiento aleatorio, en tanto el 40% restante de las fallas no
cobran un carácter aleatorio como resultado del desgaste. Las válvulas que pertenecen al
régimen de revisión son más susceptibles a las variaciones, de hecho, sólo 6 válvulas se
ajustan a la distribución exponencial. Finalmente, la revisión del sistema global da lugar a
que siendo exponencial, los intervalos extremos sugieren que la gestión reactiva es muy
irregular. Los hallazgos del estudio sugieren la implementación de estrategias de
mantenimiento diferenciadas que van desde mantenimiento predictivo para las válvulas no
exponenciales y protocolos de revisión estandarizados en el caso de revisión total que
buscan optimizar la confiabilidad del sistema.
Tenemos los datos de la válvula numero 1 como ejemplo tomamos los xn cantidades
de días que hubo entre el número de fallas por cada año como nuestras variables.
5 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
6 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
Teniendo, así como ejemplo las fallas de las válvulas teniendo en cuenta la varianza
el promedio y la desviación estándar
Vá Varian Desv.E Pro
lvula za stándar medio
15,277 123.
1 .68 123.59 44
16,356 123.
2 .44 127.89 88
10,546 123.
3 .58 102.70 80
11,260 123.
4 .03 106.11 01
11,025 120.
5 .67 105.00 34
16,269 123.
6 .72 127.55 44
7 12,332 111.05 123.
7 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
.92 63
15,838 113.
8 .44 125.85 03
11,987 120.
9 .24 109.49 94
9,059. 115.
10 83 95.18 92
Se evaluaron 34 años de las 10 válvulas que presentan fallas, el parámetro estadístico es el
resultado del promedio de la variable, varianza la desviación estándar y el nivel de
significancia que será de 0,05%.
DISTRIBUCION DE LA PROBABLILIDAD EXPONENCIAL DE LA VALVULA
8 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
De lo expuesto se deduce que el tiempo medio entre fallas para la mayoría de las
válvulas se sitúa levemente bajo los 123 días, aunque con un comportamiento relativamente
estable, ya que las válvulas 5, 8, 9 y 10 tienen un promedio de 113 a 121 días, lo que podría
indicar que la frecuencia de las fallas es más alta. Con respecto a la variabilidad, las
válvulas 2, 6 y 8 destacan por la gran dispersión, ya que tienen una desviación estándar
superior a los 125 días, lo que indica una irregularidad importante en los tiempos de falla.
Mientras que, por el contrario, la válvula 10 tiene la menor varianza (9.059,83 días²), lo que
sugiere un patrón de falla más repetido y predecible. Todos estos resultados hacen posible
detectar que hay válvulas con comportamientos atípicos que podrían ser objeto de una
estrategia de mantenimiento bien diferenciada donde se utilizaría la estrategia de una menor
frecuencia de falla o la estrategia de una gran variabilidad.
9 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
Los tiempos entre revisión se asimilan a los tiempos de fallas, dado que la mayoría
de las válvulas indican un tiempo promedio entre revisiones próximas a los 123 días. Las
válvulas 5,8,9 y 10 tienen tiempos promedio de revisión ajustados sobre un mejor periodo
en unos promedios bajos próximos a los valores de 112 – 119 días. Con respecto a la
variabilidad, las válvulas 5 y 6 son las válvulas que presentan mayor ruido en sus ciclos de
revisión, ya que ofrecen desviaciones típicas de >129 días, un dato que nos indica que
existe una programación asimétrica en la revisión de estas válvulas. Al contrario, la válvula
8 es la que tiene un mejor comportamiento, pues su varianza es la menor en todos los casos
(8,436.85 días²) y sus revisiones son las que se ejecutan de una forma más estandarizada las
una norma general. Una gran parte de los tiempos de revisión indican que es necesario
estipular protocolos para las válvulas que presentan mayor ruido, mientras que en otros
caso, como el de la válvula 8, es posible mantener la lógica de tiempos de revisión que se
habían venido incorporando.
PRUEBA DE BONDAD
La prueba de bondad se toma la variable es discreta, se toma el número de fallas por
de los 34 años con la frecuencia observada se calcula la distribución poisson, La frecuencia
10 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
esperada se encuentra al multiplicar el número de años que observamos por distribución
(Oi−Ei) 2
poisson, en chi cuadrado se calcula con la siguiente formula x =∑
2
Ei
CONCLUSIONES
En conclusión, para la válvula las 10 válvulas se acepta la hipótesis porque los
tiempos entre fallas siguen una distribución. porque el resultado de chi cuadrado estimado
es menor al valor critico de chi cuadrado. El análisis de distribución de fallas indica que
hay seis válvulas (1, 2, 3, 5, 7, 10) cuya distribución del tiempo entre fallas es exponencial,
típico de fallas aleatorias con una tasa constante, y que hay cuatro válvulas (4, 6, 8, 9) cuya
distribución no es exponencial, lo que sugiere el posible efecto del desgaste acumulativo o
del mantenimiento inadecuado. En cuanto a los mantenimientos, la heterogeneidad
observada es mayor: hay seis válvulas (1, 2, 4, 8, 9, 10) cuya distribución es exponencial, lo
11 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
que indica que el mantenimiento es aleatorio, y hay cuatro válvulas (3, 5, 6, 7) cuya
distribución es no aleatoria, probablemente porque siguen un programa de mantenimiento
preventivo.
En particular, la válvula 8, pese a su elevada tasa de fallas, presenta una revisión
exponencial, lo que indica que el mantenimiento es mayoritariamente reactivo.
Globalmente, la distribución de las revisiones totales es exponencial, por lo que sí hay
mantenimiento aleatorio en las revisiones sistémicas. Sin embargo, los intervalos de
revisión notablemente dispersos (con valores extremos de 5.91 a 740.99 días) indican una
necesidad de mejorar la planificación del mantenimiento con el fin de mitigar el riesgo
operacional
12 de 12
FACULTAD DE INGENIERÍA
INDUSTRIAL
Institución Universitaria Politécnico
PROYECTO GRUPAL
Grancolombiano
BIBLIOGRAFÍA
Bolívar Atuesta, S. (s. f.). Introducción a los procesos Estocásticos. En Variable Aleatoria.
Politécnico Grancolombiano.
Bolívar Atuesta, S. (s. f.). Introducción a los procesos Estocásticos. En Introducción A los
Procesos Estocásticos. Politécnico Grancolombiano.
Bolívar Atuesta, S. (s. f.). Cadenas de Markov en tiempo discreto. En Análisis transiente y
estable de las CMTD. Politécnico Grancolombiano.
Bolívar Atuesta, S. (s. f.). Cadenas de Markov en tiempo discreto. En Cadenas de Markov
en tiempo discreto. Politécnico Grancolombiano.
13 de 12
