APLICACIÓN DE CADENAS DE MARKOV EN TIEMPO DISCRETO PARA LA
SOLUCION DE FALLAS EN VALVULAS DE CONTROL EN OIL COLOMBIA S.A
EMPRESA DE EXTRACCIÓN PETROLERA.
Por:
JULIAN VILLEGAS CASTRILLÓN
Código: 1621980646
JHAIR JOSÉ MUENTES MARTINÉZ
Código: 16110211568
FRANCISCO JOSÉ GOMEZ SERRANO
Código: 1411025403
JULIETH FLOR MONCADA
Código: 1811982052
CINDY JOHANNA RUEDA CUJAVANTE
Código: 1621020644
Tutor:
CARLOS ANDRES QUIROZ CASTELLANOS
MARZO 2020
POLITÉCNICO GRANCOLOMBIANO
PROGRAMACIÓN ESTOCÁSTICA
Grupo [002] / 2020-1
DESARROLLO DE MODELO ESTOCÁSTICO
PRESENTACIÓN............................................................................................................................................3
ACERCA DE LA COMPAÑÍA..........................................................................................................................4
INFORME DE MODELO.................................................................................................................................5
MODELO ESTOCASTICO...............................................................................................................................7
2
PRESENTACIÓN
Este proyecto grupal, consiste en desarrollar un modelo probabilístico para el caso de
OIL Colombia S.A empresa petrolera colombiana. Se espera que el estudiante
construya una CMTD con el objetivo de modelar el sistema y calcularlas principales
medidas de desempeño. Por otra parte, también se busca que el estudiante redacte las
principales conclusiones y recomendaciones para hacer que el sistema sea más
eficiente y confiable.
Una de las aplicaciones con mayor potencial de los modelos probabilísticos, y en
particular de las cadenas de Markov, se presenta en el campo de la fiabilidad
(Reliability). En las compañías petroleras, el comportamiento de la maquinaría de
explotación es fundamental para todos los procesos de extracción y transporte de
crudo. Un derrame de crudo en alguno de los pozos se considera una catástrofe y se
debe evitar a toda costa.
La Empresa Colombiana de Petróleos está preocupada por el funcionamiento del
sistema de protección hidráulico de las torres de extracción. Dicho sistema debe
funcionar de tal manera que evite las fugas de crudo y garantice el flujo correcto en los
momentos de extracción. El sistema está compuesto por un sistema de tuberías de
escape y válvulas de control. Asuma que, dada la complejidad del sistema, el
funcionamiento de cada válvula no se puede monitorear constantemente, sino que se
revisa en ciertos intervalos de tiempo. Si en la revisión se encuentra que la válvula falla,
entonces se cambiará por una nueva (asuma que el tiempo del cambio es
despreciable).
El director de mantenimiento operacional le ha ofrecido a su compañía la oportunidad
de desarrollar un modelo probabilístico que represente, de la mejor manera posible, el
funcionamiento del sistema de protección hidráulico. Para comprobar el funcionamiento
de su modelo, se ha construido una versión simplificada de dicho sistema, compuesta
por tan sólo 10 válvulas de seguridad como se muestra en la figura 1. Cada válvula
tiene una tasa de falla diferente y en caso de falla permitiría fluir el petróleo libremente y
si éste no es interrumpido por alguna otra válvula se presentará un derrame de crudo.
Para establecer los parámetros de su modelo, se le ha proporcionado a su compañía un
archivo anexo ([Link]) con toda la información disponible de fallas y revisiones para
cada tipo de válvula.
Considere que todos los tiempos en el archivo anexo son independientes entre sí.
Además, si se presenta un derrame, se hace una revisión total, en donde todas las
válvulas serán revisadas y el sistema se dejará como nuevo. El tiempo que tarda este
proceso es independiente de todas las demás variables y su comportamiento histórico
se puede encontrar en el archivo anexo ([Link]). Finalmente, considere que cada día
3
que la torre de extracción funciona se generan, en promedio, ingresos netos por
US$40000 y cada día que la torre no funcione se tienen, en promedio, pérdidas netas
por US$ 60000.
ACERCA DE LA COMPAÑÍA
¿Quiénes somos?
OIL COLOMBIA S.A es uno de los productores de petróleo de Colombia, con un
papel dominante y representativo en la industria petrolera del país y de
Latinoamérica. Se establece bajo el título de Sociedad anónima y cuenta con un
gran potencial productivo, siempre comprometida con las buenas prácticas
ambientales, en pro de la conservación del medio ambiente.
Nuestra visión
Para el 2025 OIL COLOMBIA S.A será la empresa líder en la exploración,
explotación, refinación, transporte, almacenamiento, distribución y comercialización
de hidrocarburos, sus derivados y productos de petróleo a nivel nacional, y una de
las competitivas de Latinoamérica.
Nuestros valores
Procuramos ser la personificación de la pasión, la integridad, la innovación y el espíritu
emprendedor, esforzándonos continuamente por llevar a la práctica nuestra visión. Estos
valores nos guían en todo lo que hacemos, y son la roca en la que descansa nuestra
organización.
Objetivo
De acuerdo con los Estatutos Sociales, el objeto social de OIL COLOMBIA S.A.” es
el desarrollo, en Colombia o en el exterior, de actividades comerciales o industriales
correspondientes o relacionadas con la exploración, explotación, refinación,
transporte, almacenamiento, distribución y comercialización de hidrocarburos, sus
derivados y productos.
4
INFORME DE MODELO
El modelo probabilístico busca diseñar un sistema de protección hidráulico, para una de
las torres de producción petrolera de la compañía.
Este análisis inicial incluye la interpretación de los datos suministrados en el Excel
sobre las fallas de las 10 válvulas, las revisiones y revisiones totales, a través de una
prueba de bondad, por medio de la cual pudimos determinar la tasa de fallas por
válvula, tasa de revisiones y la tasa de revisión total, también la mejor distribución de
los datos (exponencial) la cual se halló con el Software Crystal Ball lo que facilitó
moldear el sistema de protección, las variables y el espacio.
PRUEBA DE BONDAD
La prueba de bondad se hizo utilizando la fórmula para hallar tasa de fallas la cual se
puede presentar en términos de porcentaje o con número de fallas observadas en un
tiempo de operación.
Previamente a esto habíamos convertido las fechas a tiempo absoluto (medido
en días)TASA DE FALLAS
VALVULA TOTAL DIAS # FALLAS TASA DE FALLA TIPO DE DISTRIBUCIÓN
1 12107,9383 99 0,008176454 EXPONENCIAL
2 11902,66486 99 0,008317465 EXPONENCIAL
3 18033,38201 99 0,005489819 EXPONENCIAL
4 9069,828618 99 0,010915311 EXPONENCIAL
5 9018,696681 99 0,010977196 EXPONENCIAL
6 9060,614923 99 0,010926411 EXPONENCIAL
7 18218,78212 99 0,005433953 EXPONENCIAL
8 7299,171882 99 0,013563182 EXPONENCIAL
9 6057,553512 99 0,016343232 EXPONENCIAL
10 7247,452565 99 0,013659972 EXPONENCIAL
λ Tasade Falla
TF=
Total Dias
5
TASA DE REVISIONES
VALVUL TOTAL DIAS # REVISIONES TASA DE REVISION TIPO DE DISTRIBUCIÓN
A
1 18088,26451 99 0,0054732 EXPONENCIAL
2 9106,703656 99 0,0108711 EXPONENCIAL
3 9123,964853 99 0,0108505 EXPONENCIAL
4 9091,530332 99 0,0108893 EXPONENCIAL
5 9103,220096 99 0,0108753 EXPONENCIAL
6 9021,004393 99 0,0109744 EXPONENCIAL
7 9106,703656 99 0,0108711 EXPONENCIAL
8 17911,54501 99 0,0055272 EXPONENCIAL
9 18137,46833 99 0,0054583 EXPONENCIAL
10 18015,78354 99 0,0054952 EXPONENCIAL
λ Tasa de Revisiones
TR=
Total Dias
TASA DE REVISION TOTAL
TOTAL DIAS TASA DE
REVISIONES REVISIONES TOTAL
99 699,92077 0,1414446
6
MODELO ESTOCASTICO
Una cadena de Markov es un proceso estocástico en tiempo discreto { X n : n≥ 0 } con
espacio discreto, tal que, para cualquier entero n ≥ 0, y para cualquier estado S0 , … S n S n+1
, de este modo se cumple:
Pr { S n+1|S 0 , … , S n }=Pr { S n+1|Sn }
Siendo n+1 tiempo futuro, n tiempo presente y 0,1 , … . ,n−1 como el pasado. Entonces
la condición anterior establece que la distribución de probabilidad del estado del
proceso al tiempo futuro n+1 depende únicamente del estado del proceso al tiempo n, y
no depende de los estados en los tiempos pasados.
Variables
A raíz de esto se determina que las variables de nuestro sistema son:
B(t) = El estado del sistema total
Espacio de estados (S):
S ( B(t ))
¿ {H =cuando No se necesita Revisión Total , F=cuando Si necesita RevisiónTotal }
Condiciones de transición:
Estado Inicial Estado Final Condición de transición
H H Ninguna válvula falla.
Falla alguna de las válvulas: 2, 3, 4,
5, 6, 8.
Se repara alguna de las válvulas: 2,
3, 4, 5, 6,8.
H F Fallan alguna de las válvulas 1, 7, 9,
10.
F H Se reparan todas las válvulas.
F F No hay transición posible.
