3-Entrega Proyecto
3-Entrega Proyecto
PRESENTADO POR:
Subgrupo N° 27
Tutor:
Mayo 2020
Programación Estocástica
TABLA DE CONTENIDO
LISTADO DE TABLAS
Se pone a consideración en este proyecto la situación por la que está pasando las compañías
petroleras, la cual está dada por las pérdidas económicas considerables y daños al medio
ambiente que están afectando a este sector industrial debido a los derrames de crudo que
están presentando, todo esto debido a la falta de un sistema de protección hidráulico el cual
posee unos tubos que extraen el crudo y son controlados por 10 válvulas a las cuales se les ha
requerido un plan de mantenimiento exhaustivo.
La compañía ha decidido contratar profesionales para realizar un análisis estadístico y
probabilístico en todo su sistema de extracción del crudo para evitar estas pérdidas y daño al
medio ambiente, el cual iniciaremos analizando el histórico de fallas y el plan de verificación
que posee la compañía por cada válvula para evidenciar cual es la probabilidad de que
vuelvan a fallar las válvulas y cuáles son las más propensas a fallas.
Para la segunda entrega de este proyecto se plantea la siguiente situación:
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 4 válvulas (1, 8, 9, 10) de seguridad como se
muestra en la figura.
OBJETIVOS
OBJETIVO PRIMERA ENTREGA:
Modelar sistemas reales, a través de la caracterización probabilística de estos, como Cadenas
de Markov.
Acorde con los datos históricos entregados se trabajará con el registro histórico de fallas, con
estos datos se pueden dar respuestas a las preguntas planteadas. El modelo probabilístico
elegido es una cadena de Markov de tiempo discreto; en donde la probabilidad del evento
futuro está condicionada únicamente por el evento inmediatamente anterior. La variable de
estado es:
X n ,d =Estado de la válvula n en el día d
n=( 1,2,3,4,5,6,7,8,9,10 )
d ≥0
Teniendo dos posibles estados, bueno (B) o malo (M).
De esta manera, el espacio muestral para el sistema resulta ser una combinación de los
posibles estados (que son 2) de las válvulas (10 en total), siendo 210 posibles estados (1024
puntos muéstrales). Sin embargo, el espacio muestral se puede restringir a las combinaciones
que resultan en un derrame, ya que lo solicitado por el cliente tiene que ver con el sistema
cuando sucede un derrame. Limitando el espacio muestral a las combinaciones resultantes en
derrames se puede responder a las preguntas planteadas. Dichas combinaciones son
S={( M , B , B , B , B , B , B , B , B , B ) , ( M , B , B , B , B , B , B , M , M , B ) ,… ,
, … , ( M , M , M , B , M , M , M , M , M , M ) , (M , M , M , M , M , M , M , M , M , M )
Considerando todos los escenarios donde puede ocurrir derrame se tiene un espacio muestral
con 109 posibles escenarios. Este espacio muestral puede reducirse aún más cuando se
consideren las probabilidades de que alguno de esos eventos suceda.
Definiendo Xn,d de la forma hecha y limitando el espacio muestral a los puntos de interés
podremos predecir si el sistema sufrirá un derrame de petróleo en el día n, como también se
podrán estimar los días que pasarán antes de que ocurra un derrame.
1 Paso:
El primer paso que se realiza es analizar los datos de entrada para así hallar un modelo
probabilístico que nos ayude a llegar al potencial de falla, lo realizamos de la siguiente
manera; restamos la medición actual de cada casilla con la inmediatamente anterior para
hallar P luego pasamos la información que está en formato de fecha a minutos luego
dividimos por la suma total de la revisión.
Programación Estocástica
2 Paso:
Con los datos obtenidos elaboramos una matriz estocástica que nos ayude a definir que
válvulas presentan el mayor tiempo de revisión es decir las válvulas que han presentado más
fallas y con esto la causa de las pérdidas:
Programación Estocástica
3 Paso:
Realizamos las pruebas de bondad y ajuste para analizar los datos y encontrar la zona con
mayores dificultades:
Acorde con los datos entregados (valores lambda) se trabajará con el fin de dar respuesta a las
preguntas planteadas. El modelo probabilístico elegido es una cadena de Markov de tiempo
discreto; en donde la probabilidad del evento futuro está condicionada únicamente por el
evento inmediatamente anterior.
1. Con los datos brindados en este documento genere el modelo específico para el sistema
con las 4 válvulas (1, 8, 9, 10). Debe definir la variable de estado Xn, el espacio de estados S
y la matriz de probabilidades de transición de un paso (Guíese por el ejemplo resuelto).
Estas son las variables de estado de las válvulas las cuales pueden estar de 2 maneras que
son; B (Buena) y M (Mala), en el día n.
Donde,
La variable A es la válvula 1
La variable B es la válvula 8
La variable C es la válvula 9
La variable D es la válvula 10
S=BBBB,BBBM,BBMB,BBMM,BMBB,BMBM,BMMB,BMMM,MBBB,MBBM,MBM
B,MBMM,MMBB,MMBM,MMMB,MMMM
En el cual encontramos 16 posibles eventos o combinaciones por las que pueden estar las
válvulas
Programación Estocástica
Para hallar la matriz de probabilidades de transición primero se hallan las probabilidades (P)
de cada evento con la fórmula: FT (t) = P (T ≤ t) = 1 − e ^−λt, para t ≥ 0.
Luego se multiplican las probabilidades que ocurra el evento por el que pasa cada válvula
desde el estado inicial (columna de estados vertical) al estado final (fila de estados
horizontal).
Programación Estocástica
2.1.3 El número de días, en promedio, que se habrán utilizado para el proceso de revisión
total justo antes de cada auditoria.
Para hallar este número de días en promedio que se habrían utilizado para revisiones
tendríamos que encontrar primero el estado estable que para nuestro ejercicio seria en el día
256, esto lo hallamos elevando la matriz de transición de probabilidades de un paso a la 256,
se sabe que esta en estado estable porque todas sus probabilidades de una columna son
iguales, así las cosas decimos que para la primera auditoria obtenemos un numero promedio
de días para el proceso de revisión de:
2.1.4 El número de días, en promedio, que transcurrirán hasta el próximo derrame de crudo.
Para la solución de este punto aplicamos la metodología de tiempos de primera pasada, en el
cual obtenemos que el número de días en promedio que transcurrirían para el próximo
derrame seria de:
Para mayor detalle por favor revisar el archivo anexo Excel. (VER ANEXO)
3. Un análisis económico de largo plazo, que muestre el valor esperado de los ingresos
diarios para una torre de extracción que cuente con el sistema de protección hidráulico.
Programación Estocástica
(Asuma que la única causa para que la torre no funcione es la falla del sistema de protección
hidráulico).
Sabiendo que:
Lo primero que se realiza para la solución de este punto es encontrar las revisiones totales
antes de llegar a estado estable es decir en los 256 días y encontramos que se realizan:
Revisiones Totales:
Sabemos que la revisión total dura un promedio de 7 días de acuerdo a los datos
históricos que se envían en el documento Excel.
Las pérdidas que se generan por NO tener una torre funcionando OK son de US$
60.000 por día.
Las ganancias que se generan por tener la torre funcionando correctamente son de
US$ 40.000 por día.
Tenemos 101.28 días en que la torre NO funciona correctamente por derrame de
crudo y mantenimiento del sistema de válvulas.
Tenemos 154.72 días en que la torre funciona correctamente.
Luego de obtener estos datos procedemos a hallar solución a esta pregunta de los costos
generados y de perdidas:
Programación Estocástica
CONCLUSIONES
Encontramos que el recorrido donde están instaladas las válvulas 1, 8, 9, 10 tiene una gran
oportunidad de mejora debido a que tienen promedios de revisiones altos y son las válvulas
que están consecutivas según el plano entregado por la compañía
El diseño del sistema no es el adecuado ya que en este recorrido crítico se está dejando
mucha capacidad de la extracción del crudo en la válvula 9 ya que esta recibe el paso de
crudo de las válvulas 1 y 8 y cuando se satura sobre carga la válvula 10.
ANEXOS
ANEXO.xlsx