FACULTAD DE INGENIERÍA Y CIENCIAS

BÁ SICAS

PROYECTO GRUPAL Institució n Universitaria Politécnico Grancolombiano

PROGRAMACIÓN ESTOCASTICA

TRABAJO COLABORATIVO
Diseño de un sistema de protección hidráulico para un pozo petrolero

TUTOR
CARLOS ANDRES QUIROZ

INTEGRANTES

INSTITUCIÓN UNIVERSITARIA POLITECNICO GRANCOLOMBIANO

INGENERIA INDUSTRIAL
BOGOTA D.C
2020
 FACULTAD DE INGENIERÍA Y CIENCIAS
BÁ SICAS

PROYECTO GRUPAL Institució n Universitaria Politécnico Grancolombiano

ENUNCIADO DETALLADO PARA LA SEGUNDA ENTREGA:

Descripción del Proyecto Una de las aplicaciones con mayor potencial de los modelos
probabilísticos, y en particular de las cadenas de Márkov, 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 4
válvulas (1, 8, 9, 10) de seguridad como se muestra en la figura 1.

RESUELVA

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.

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á
FACULTAD DE INGENIERÍA Y CIENCIAS
BÁ SICAS

PROYECTO GRUPAL Institució n Universitaria Politécnico Grancolombiano

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 auditoría el sistema esté

funcionando como nuevo.
2.1.2 La probabilidad de que en la segunda auditoría el sistema esté
funcionando como nuevo.
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 auditoría.
2.1.4 El número de días, en promedio, que transcurrirán hasta el próximo
derrame de crudo.
 FACULTAD DE INGENIERÍA Y CIENCIAS
BÁ SICAS

PROYECTO GRUPAL Institució n Universitaria Politécnico Grancolombiano

VARIABLE DE ESTADO
A ( t )=[ Estado de la válvula1]

B ( t )=[ Estado de la válvula 8]

X n={(A9]( t ) B ( t ) C ( t ) D (t ) )}
C ( t )=[Estado de la válvula

D ( t ) =[ Estado de la válvula 10]

La variable será estudiada cada día

A ( t )= [ MBsisilala válvula1
válvula 1 está funcionando
no está funcionando ]

B ( t )=[ B si la válvula 8 está funcionando ]

M sila válvula 8 no está funcionando

C ( t )=[ B si laválvula 9 está funcionando ]

M si la válvula 9 no está funcionando

B sila válvula 10 está funcionando

D ( t ) =[
M sila válvula 10 no está funcionando ]

S={( BBBB),(BBBM) ,( BBMB),( BBMM ) ,(BMBB),( BMBM ) ,(BMMB),

( BMMM ) , ( MBBB ) , ( MBBM ) , ( MBMB ) , ( MBMM ) , ( MMBB ) , ( MMBM ) ,

( MMMB),( MMMM )}

Rojo = Presencia de derrame de crudo

Negro = No hay ningún derrame de crudo
 FACULTAD DE INGENIERÍA Y CIENCIAS
BÁ SICAS

PROYECTO GRUPAL Institució n Universitaria Politécnico Grancolombiano

TASAS DE FALLAS, REPARACIÓN Y REVISIÓN TOTAL

Válvulas Tiempo entre fallas Válvulas Tiempo entre fallas

se distribuye se distribuye
exponencial exponencial
1 μ1=0,008 fallas /día 1 μ1=0,008 fallas /día

8 μ8=0,014 fallas /día 8 μ8=0,014 fallas /día

9 μ9 =0,016 fallas /día 9 μ9 =0,016 fallas/día

10 μ10=0,014 fallas/día 10 μ10=0,014 fallas/día

Tiempo de revisión general se distribuye exponencial

Todo el sistema λ=0,143 sistemas /día

Posteriormente, se calculan las probabilidades de los acontecimientos individuales que pueden

ocurrir en el transcurso de un día.
PROBABILIDAD DE LOS EVENTOS INDVIDUALES

Tiempo menor a un día, se distribuye exponencial

Probabilidad de que falle la válvula P ( X <1 )=1−e−0,008 (1)=0,008
1

Probabilidad de que falle la válvula P ( X <1 )=1−e−0,014 (1)=0,014

8
 FACULTAD DE INGENIERÍA Y CIENCIAS
BÁ SICAS

PROYECTO GRUPAL Institució n Universitaria Politécnico Grancolombiano

Probabilidad de que falle la válvula P ( X <1 )=1−e−0,016 (1)=0,016

9

Probabilidad de que falle la válvula P ( X <1 )=1−e−0,014 (1)=0,014

10

Probabilidad de que la válvula 1 se P ( X <1 )=1−e−0,005 (1)=0,005

haya revisado y reparado
Probabilidad de que la válvula 8 se P ( X <1 )=1−e−0,005 (1)=0,005
haya revisado y reparado
Probabilidad de que la válvula 9 se P ( X <1 )=1−e−0,005 (1)=0,005
haya revisado y reparado
Probabilidad de que la válvula 10 se P ( X <1 )=1−e−0,005 (1)=0,005
haya revisado y reparado
Probabilidad de todo el sistema se P ( X <1 )=1−e−0,143 (1)=0,133
haya revisado y reparado

Como las válvulas son independientes las unas de las otras la probabilidad de ocurrencia de los
eventos es independiente. Por lo tanto, la probabilidad de ocurrencia de un evento conjunto es la
multiplicación de las probabilidades de cada evento individual que esté involucrado. Para este
problema, los eventos de interés son 18, como se muestra a continuación:

PROBABILIDAD DE OCURRENCIA

Evento Probabilidad de que ocurra el

evento
E1: la válvula 1 falla en un tiempo menor a 1 día P( E 1)=0,008
E2: la válvula 1 falla en un tiempo mayor a 1 día P( E 2)=1−P(E 1)=0,992
E3: la válvula 8 falla en un tiempo menor a 1 día P( E 3)=0,014
E4: la válvula 8 falla en un tiempo mayor a 1 día P( E 4)=1−P( E3)=0,986
E5: la válvula 9 falla en un tiempo menor a 1 día P( E 5)=0,016
E6: la válvula 9 falla en un tiempo mayor a 1 día P( E 6)=1−P( E 5)=0,984
E7: la válvula 10 falla en un tiempo menor a 1 día P( E 7)=0,014
E8: la válvula 10 falla en un tiempo mayor a 1 día P( E 8)=1−P( E 7)=0,986
 FACULTAD DE INGENIERÍA Y CIENCIAS
BÁ SICAS

PROYECTO GRUPAL Institució n Universitaria Politécnico Grancolombiano

E9: la válvula 1 se revisa y repara en un tiempo menor a 1 día P( E 9)=0,005

E10: la válvula 1 se revisa y repara en un tiempo mayor a 1 P( E 10)=1−P ( E 9)=0,995
día
E11: la válvula 8 se revisa y repara en un tiempo menor a 1 P( E 11)=0,005
día
E12: la válvula 8 se revisa y repara en un tiempo mayor a 1 P( E 12)=1−P(E 11)=0,995
día
E13: la válvula 9 se revisa y repara en un tiempo menor a 1 P( E 13)=0,005
día
E14: la válvula 9 se revisa y repara en un tiempo mayor a 1 P( E 14)=1−P( E 13)=0,995
día
E15: la válvula 10 se revisa y repara en un tiempo menor a 1 P( E 15)=0,005
día
E16: la válvula 10 se revisa y repara en un tiempo mayor a 1 P( E 16)=1−P( E 15)=0,995
día
E17: Todo el sistema se revisa y repara en un tiempo menor a P( E 17)=0,133
1 día,
E18: Todo el sistema se revisa y repara en un tiempo mayor a P( E 18)=1−P ( E 17)=0,867
1 día,

Con la tabla anterior, es posible conformar una matriz de transiciones que modele el sistema y
sea capaz de calcular probabilidades de eventos en un tiempo cercano. A continuación, se
muestra la utilidad de dicha matriz:
El director de operaciones quiere instalar el sistema de protección hidráulico en un nuevo pozo
debido a que acaba de presentar un derrame de crudo (estado MMMM) y se realizarán dos
auditorías en dos periodos de tiempo diferentes (3 y 6 meses).
 FACULTAD DE INGENIERÍA Y CIENCIAS
BÁ SICAS

PROYECTO GRUPAL Institució n Universitaria Politécnico Grancolombiano

(BMMM) (MBBB) (MBBM) (MBMB) (MBMM) (MMBB) (MMBM) (MMMB) (MMMM)

(BBBB) 0,000 0,008 0,000 0,000 0,000 0,000 0,000 0,000 0,000
(BBBM) 0,000 0,000 0,008 0,000 0,000 0,000 0,000 0,000 0,000
(BBMB) 0,000 0,000 0,000 0,008 0,000 0,000 0,000 0,000 0,000
(BMBM) 0,016 0,000 0,000 0,000 0,000 0,000 0,008 0,000 0,000
(BMMB) 0,014 0,000 0,000 0,000 0,000 0,000 0,000 0,008 0,000
(BBMM) 0,014 0,000 0,000 0,000 0,008 0,000 0,000 0,000 0,000
(BMBB) 0,000 0,000 0,000 0,000 0,000 0,008 0,000 0,000 0,000
(BMMM) 0,977 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,008
(MBBB) 0,000 0,867 0,000 0,000 0,000 0,000 0,000 0,000 0,000
(MBBM) 0,000 0,000 0,867 0,000 0,000 0,000 0,000 0,000 0,000
(MBMB) 0,000 0,000 0,000 0,867 0,000 0,000 0,000 0,000 0,000
(MBMM) 0,000 0,000 0,000 0,000 0,867 0,000 0,000 0,000 0,000
(MMBB) 0,000 0,000 0,000 0,000 0,000 0,867 0,000 0,000 0,000
(MMBM) 0,000 0,000 0,000 0,000 0,000 0,000 0,867 0,000 0,000
(MMMB) 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,867 0,000
(MMMM) 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,867

PRIMERA AUDITORIA
En la primera auditoría se quiere saber la probabilidad de que el sistema esté funcionando como
nuevo, en otras palabras, que tan probable es que se pase del estado MMMM a un estado BBBB
en el día 90. Al multiplicar la matriz de transición de un paso 90 veces se puede obtener la matriz
en el día 90:
Pˆ90 = 0,1906
 FACULTAD DE INGENIERÍA Y CIENCIAS
BÁ SICAS

PROYECTO GRUPAL Institució n Universitaria Politécnico Grancolombiano

La probabilidad de que en la primera auditoría el sistema esté funcionando como nuevo es igual
al 19,06%

(BBBB) (BBBM) (BBMB) (BBMM) (BMBB) (BMBM) (BMMB) (BMMM)

(BBBB) 0,186 0,091 0,107 0,104 0,102 0,103 0,103 0,148
(BBBM) 0,173 0,103 0,085 0,124 0,076 0,123 0,077 0,183
(BBMB) 0,177 0,076 0,113 0,104 0,104 0,104 0,105 0,162
(BMBM) 0,167 0,080 0,082 0,143 0,077 0,090 0,084 0,220
(BMMB) 0,173 0,071 0,089 0,087 0,125 0,105 0,111 0,183
(BBMM) 0,167 0,080 0,082 0,090 0,084 0,143 0,078 0,219
(BMBB) 0,173 0,071 0,089 0,105 0,110 0,087 0,126 0,183
(BMMM) 0,165 0,068 0,076 0,101 0,085 0,101 0,086 0,263
(MBBB) 0,191 0,094 0,111 0,102 0,103 0,102 0,104 0,138
(MBBM) 0,191 0,094 0,111 0,102 0,103 0,102 0,104 0,138
(MBMB) 0,191 0,094 0,111 0,102 0,103 0,102 0,104 0,138
(MBMM) 0,191 0,094 0,111 0,102 0,103 0,102 0,104 0,138
(MMBB) 0,191 0,094 0,111 0,102 0,103 0,102 0,104 0,138
(MMBM) 0,191 0,094 0,111 0,102 0,103 0,102 0,104 0,138
(MMMB) 0,191 0,094 0,111 0,102 0,103 0,102 0,104 0,138
(MMMM) 0,191 0,094 0,111 0,102 0,103 0,102 0,104 0,138

SEGUNDA AUDITORIA
De igual manera, en una segunda auditoría se quiere saber la probabilidad de que el sistema esté
funcionando como nuevo, en otras palabras, que tan probable es que se pase del estado MMMM
a un estado BBBB en el sexto mes. Al multiplicar la matriz de transición del paso 90 por ella
misma se puede obtener la matriz en el día 180:
P^180 = 0,1727
Por lo tanto, se determina que la probabilidad de que en la segunda auditoría el sistema esté
funcionando como nuevo es igual al 17,27%
FACULTAD DE INGENIERÍA Y CIENCIAS
BÁ SICAS

PROYECTO GRUPAL Institució n Universitaria Politécnico Grancolombiano

(BBBB) (BBBM) (BBMB) (BBMM) (BMBB) (BMBM) (BMMB) (BMMM)

(BBBB) 0,175 0,081 0,092 0,106 0,096 0,106 0,097 0,191
(BBBM) 0,174 0,081 0,091 0,107 0,094 0,107 0,094 0,196
(BBMB) 0,175 0,080 0,092 0,106 0,097 0,106 0,097 0,193
(BMBM) 0,174 0,080 0,090 0,108 0,094 0,105 0,095 0,199
(BMMB) 0,174 0,080 0,091 0,105 0,097 0,106 0,098 0,194
(BBMM) 0,174 0,080 0,090 0,105 0,094 0,108 0,095 0,199
(BMBB) 0,174 0,080 0,091 0,106 0,097 0,104 0,098 0,194
(BMMM) 0,173 0,079 0,090 0,106 0,094 0,105 0,095 0,201
(MBBB) 0,175 0,081 0,093 0,106 0,097 0,106 0,097 0,190
(MBBM) 0,175 0,081 0,093 0,106 0,097 0,106 0,097 0,190
(MBMB) 0,175 0,081 0,093 0,106 0,097 0,106 0,097 0,190
(MBMM) 0,175 0,081 0,093 0,106 0,097 0,106 0,097 0,190
(MMBB) 0,175 0,081 0,093 0,106 0,097 0,106 0,097 0,190
(MMBM) 0,175 0,081 0,093 0,106 0,097 0,106 0,097 0,190
(MMMB) 0,175 0,081 0,093 0,106 0,097 0,106 0,097 0,190
(MMMM) 0,175 0,081 0,093 0,106 0,097 0,106 0,097 0,190

CONCLUSIONES
