Programación Estocástica

DISEÑO DE UN SISTEMA DE PROTECCIÓN HIDRÁULICO PARA UN POZO

PETROLERO

PRESENTADO POR:

Anderson Fabián Murcia Fandiño

Héctor Antonio Rincón Herrera
Iván Darío Hurtado Lara
Judith Daniela Sierra Barragán
Millerlay Moreno Calderón

Proyecto de Programación Estocástica

Subgrupo N° 27

Tutor:

Carlos Andrés Quiroz Castellanos

Institución Universitaria Politécnico Gran Colombiano

Facultad De Ingeniería Y Ciencias Básicas
Campus Virtual

Mayo 2020
 Programación Estocástica

TABLA DE CONTENIDO

1 DESCRIPCIÓN DEL PROBLEMA ................................................................... 2

2 OBJETIVOS ..................................................................................................... 3
3 SOLUCION Y CORRECIONES PRIMERA ENTREGA………………………..5
4 SOLUCION SEGUNDA ENTREGA……………………………………………. 9
5 SOLUCION TERCERA ENTREGA ……………………………………….. 13
6 RECOMENDACIONES Y CONCLUSIONES............................................. …15

LISTADO DE TABLAS

Tabla 1. Análisis de datos de entrada ............................................................................. 3

Tabla 2. Matriz estocástica ............................................................................................ 6
Tabla 3. Pruebas de bondad y ajuste .............................................................................. 7
Tabla 4.Probabilidad de Evento………………………………………………………...10

Tabla 5. Matriz de Probabilidades………………………………………………………11

Tabla 6. Auditoria 1……………………………………………………………………..11
Tabla 7.Auditoria 2……………………………………………………………………..12
 Programación Estocástica

DESCRIPCIÓN DEL PROBLEMA

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.

Se resolverán diferentes situaciones con las variaciones de estas válvulas.

 Programación Estocástica

OBJETIVOS
OBJETIVO PRIMERA ENTREGA:
Modelar sistemas reales, a través de la caracterización probabilística de estos, como Cadenas
de Markov.

OBJETIVO SEGUNDA ENTREGA:

• Calcular correctamente las medidas de desempeño relacionadas con el sistema modelado
• Analizar y proponer alternativas de solución, con base en el diagnóstico realizado y los
resultados del modelo probabilístico.
 Programación Estocástica

SOLUCION PRIMERA ENTREGA

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:
𝑿𝒏,𝒅 = 𝐸𝑠𝑡𝑎𝑑𝑜 𝑑𝑒 𝑙𝑎 𝑣á𝑙𝑣𝑢𝑙𝑎 𝒏 𝑒𝑛 𝑒𝑙 𝑑í𝑎 𝒅

𝑛 = (1,2,3,4,5,6,7,8,9,10)
𝑑≥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
𝑆 = {(𝑀, 𝐵, 𝐵, 𝐵, 𝐵, 𝐵, 𝐵, 𝐵, 𝐵, 𝐵), (𝑀, 𝐵, 𝐵, 𝐵, 𝐵, 𝐵, 𝐵, 𝑀, 𝑀, 𝐵), … ,
, … , (𝑀, 𝑀, 𝑀, 𝐵, 𝑀, 𝑀, 𝑀, 𝑀, 𝑀, 𝑀), (𝑀, 𝑀, 𝑀, 𝑀, 𝑀, 𝑀, 𝑀, 𝑀, 𝑀, 𝑀)
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

Tabla 1. Análisis de datos de entrada

Fuente: Los autores

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

Tabla 2. Matriz estocástica

Fuente: Los autores

3 Paso:
Realizamos las pruebas de bondad y ajuste para analizar los datos y encontrar la zona con
mayores dificultades:

Tabla 3. Pruebas de bondad y ajuste

Programación Estocástica
 Programación Estocástica

SOLUCION SEGUNDA ENTREGA

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).

- La variable de estado para las 4 válvulas (1,8,9,10) es:

O en términos de las válvulas de la siguiente manera:

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

- Los espacios de estado para este ejercicio serían los siguientes:

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

Matriz de probabilidades de transición:

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.

Dándonos como resultado:

Tabla 4. Probabilidades del evento

Fuente: Los autores

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

Tabla 5. Matriz de probabilidad

Fuente: Los autores

2. Adicionalmente, el director de mantenimiento operacional instalará el sistema de protección

hidráulico en un nuevo pozo que acaba de presentar un derrame de crudo (estado MMM) (por
fortuna, sin consecuencias mayores para la compañía) y realizará dos auditorías para
comprobar su funcionamiento. La primera en tres meses y la segunda en 6 meses (Asuma que
un mes tiene exactamente 30 días). Por lo tanto, quiere que su modelo calcule:
2.1.1 La probabilidad de que en la primera auditoria el sistema esté funcionando como nuevo.

La probabilidad de que en la primera auditoria el sistema está funcionando correctamente es

de un 19% para esto se eleva la matriz P^1 (matriz de probabilidades de transición hallada
anteriormente) hasta llegar al número de días en que se realiza la primera auditoria (90 días) y
buscar la casilla donde las 4 válvulas están MMMM como estado inicial al estado final BBBB.
Donde lo que simulamos es que el sistema pase de estar Malo a estar Bueno.
Tabla 6. Resultado auditoria 1

Fuente: Los autores

 Programación Estocástica

2.1.2 La probabilidad de que en la segunda auditoria el sistema esté funcionando como nuevo.

La probabilidad de que en la segunda auditoria el sistema está funcionando correctamente es

de un 17% para esto se eleva la matriz P^1 (matriz de probabilidades de transición hallada
anteriormente) hasta llegar al número de días en que se realiza la segunda auditoria (180 días)
y buscar la casilla donde las 4 válvulas están MMMM como estado inicial al estado final
BBBB. Donde lo que simulamos es que el sistema pase de estar Malo a estar Bueno.
Tabla 7. Resultado auditoria 2

Fuente: Los autores

 Programación Estocástica

SOLUCION TERCERA ENTREGA

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. (Asuma que la
única causa para que la torre no funcione es la falla del sistema de protección hidráulico).
Sabiendo que:
 Programación Estocástica

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

Ilustración 1. Sistema crítico

Fuente: Los autores

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.

La distribución de las válvulas 2, 3, 4, 5, 6 y 7 se adecua más a un sistema óptimo, aunque al

final del ciclo se tiene una única válvula encargada de dirigir el crudo a la superficie.

ANEXOS

ANEXO.xlsx