PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
PROGRAMACIÓN ESTOCASTICA
DISEÑO DE UN SISTEMA DE PROTECCIÓN HIDRÁULICO PARA UN POZO
PETROLERO
INTEGRANTES:
ANDRES MAURICIO ARENAS PAREDES COD: 100010521
NATHALY TICORA OBANDO COD: 100080914
NATALIA SANCHEZ AGUILAR COD: 1511980689
LUIS MIGUEL GARZON SCARPETTA COD: 1611020414
PRESENTADO A:
EDWARD PARRA FLOREZ
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
INSTITUCIÓN UNIVERSITARIA POLITECNICO GRAN COLOMBIANO
FACULTAD DE INGENIERIA, DISEÑO E INNOVACIÓN
SEPTIEMBRE 2021
DESCRIPCIÓN DEL PROBLEMA
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, además si este no es
interrumpido por alguna otra válvula, se presentará un derrame de crudo.
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
CONSTRUCCIÓN DEL MODELO PROBABILISTICO
Inicialmente debemos encontrar a que tipo de distribución de probabilidad está relacionado el
funcionamiento del sistema, teniendo en cuenta los datos de tiempos de falla, tiempos de revisión
de cada válvula y la revisión general del sistema.
Para esto utilizamos el procedimiento de Pruebas de bondad de ajuste realizando lo siguiente:
1. Definimos la hipótesis nula (Ho) y la hipótesis alternativa (Ha)
Ho = Tiempo de falla válvula n y que sigue una distribución exponencial
Ha = Tiempo de falla válvula n y que no sigue una distribución exponencial
Hallamos el máximo y mínimo de los datos para calcular la amplitud del intervalo, también se
calcula el valor de k que es el número de clases, esto se realizó para cada válvula (en el ejemplo
esta dado para la válvula 1)
X min 0
X máx 1194
N (número de muestras) 100
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
Amplitud Intervalo 156
X max −X min
La amplitud del intervalo esta dado por: ∆ X=
k
K esta dado por la fórmula de Sturges como:k =1+3,32∗log 10 (n)
k (# Clases) 7,64
Procedemos a elaborar los cálculos que nos permitirá encontrar el valor estadístico de prueba,
que ayuda a determinar si el comportamiento de los datos corresponde a una distribución
exponencial.
Empleamos la siguiente tabla e iremos describiendo como se halla cada valor contenido en ella.
Rango Límite Limite Clase Frecuenci Probabilidad Frecuenci Estadístico
Inferior Superio a Exponencial a Chi-
r Observada Esperada Cuadrado
(Ei)
1 0 156 0- 76 0,7226 72 0,1934
156
2 157 313 157 - 15 0,2011 20 1,2973
313
3 314 470 314 - 9 0,0763 8 0,2454
470
4
5
6
7
8
100 1 100 1,7362
Los rangos del 5 al 8 se omiten ya que se ajustan las frecuencias para que su valor sea al menos
de 5. Ahora se procede a realizar la descripción de la obtención de cada uno de los datos de la
tabla:
Rango: De acuerdo con el número de clases que se calcula con la Fórmula de Sturges descrita
anteriormente.
Limites: La amplitud del intervalo define los límites para cada uno de los rangos, este dato nos
permite establecer los límites para cada rango
Grados de libertad: Es el total de clases del rango menos 1 con las que se trabaja.
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
Significancia: Se define un intervalo de confianza de 95%, por lo que el valor es del 5%.
Valor Critico: Se elige de las tablas de chi-Cuadrado, al tener los grados de libertad y el valor
de significancia definido, este valor se muestra en la siguiente imagen:
Con los valores ya definidos se tienen los siguientes datos:
Grados de Libertad 2
Significancia 0,05
Valor Critico 5,991
5
Por otra parte, se calcula el tiempo promedio de falla que es sumar todas las muestras y dividirlas
por 100 que es el total de muestras por válvula.
La tasa de falla que es el inverso del promedio hallado anteriormente. Los resultados son:
Tiempo Promedio de 122 días
falla
Tasa 0,00822 falla/día
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
Frecuencia observada: Es la cantidad de datos dentro de los limites que se tiene para cada uno
de los rangos, para esto se hace un conteo de datos que este entre cada uno
De los datos para la válvula 1 se tienen los siguientes resultados:
Clase Frecuencia Observada
0 - 156 76
157 - 313 15
314 - 470 9
Probabilidad Exponencial: Se elige asumiendo que los datos son exponenciales y con el uso de
la fórmula de DISTR.EXP.N en Excel se hallan cada uno de los valores. Para el segundo valor
se resta el dato anterior para el tercero la unidad menos los datos anteriores obtenidos
Frecuencia Esperada: Es el producto de cada valor de la probabilidad exponencial
multiplicada por el total de muestras
Estadístico Chi-Cuadrado: Ya con los datos anteriores se emplea la ecuación de Chi-Cuadrado
que es la siguiente:
k
( Frec esp−Frecobs )2
Chi−Cuadrado=∑
i=1 Frec obs
El resultado anterior es nuestro estadístico de prueba donde probamos si la distribución es o no
exponencial de acuerdo con el siguiente criterio
Si el Estadístico de prueba es mayor al valor critico elegido la Hipótesis nula es rechazadas.
Para cada uno de los datos hallados para cada válvula, el valor critico siempre fue mayor, por
lo cual la distribución es exponencial.
Ya con todo realizado elaboramos la tabla con la tasa y el tiempo promedio de falla y reparación
de cada válvula y también para el total de reparaciones
Estas son las siguientes
Válvula Distribución Tasa Tiempo
(Fallas/día) Promedio
Falla (días)
1 Exponencial 0,00822 122
2 Exponencial 0,00822 122
3 Exponencial 0,00548 182
4 Exponencial 0,01096 91
5 Exponencial 0,01096 91
6 Exponencial 0,01096 91
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
7 Exponencial 0,00548 182
8 Exponencial 0,01370 73
9 Exponencial 0,01644 61
10 Exponencial 0,01370 73
Tabla de Fallas y distribución de probabilidad
Válvula Distribuci Tasa Tiempo
ón (Revisión/dí Promedio
a) Revisión
(días)
1 Exponencia 0,00548 182
l
2 Exponencia 0,01096 91
l
3 Exponencia 0,01096 91
l
4 Exponencia 0,01096 91
l
5 Exponencia 0,01096 91
l
6 Exponencia 0,01096 91
l
7 Exponencia 0,01096 91
l
8 Exponencia 0,00548 182
l
9 Exponencia 0,00548 182
l
10 Exponencia 0,00548 182
l
Tabla de Revisiones y distribución de probabilidad
Válvula Distribución Tasa Tiempo
(Revisión Promedio
Total/día) Revisión
Total
(días)
1 Exponencial 0,14287 7
Tabla de Reparación total y distribución de probabilidad
PROGRAMACIÓN ESTOCASTICA
Ya con
los
datos
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
definidos como distribución exponencial procedemos a realizar el modelamiento del sistema.
Tomamos como variables de estado el comportamiento del sistema protección de extracción
hidráulico W(t) en función del estado de las válvulas determinando como punto importante, si
hay o no derrame por el daño de algunas válvulas que es el objetivo del modelo.
Estados de transición
Para reducir la cantidad de posibles combinaciones para reparar algunas de las válvulas se
agrupan estas en estados donde implicaría hacer una intervención o no del sistema, sería lo
siguiente:
X = Ninguna válvula requiere reparación y no hay derrames de crudo
Y = Algunas válvulas requieren reparación y no hay derrames de crudo
Z = Válvula fallando requiere reparación y hay derrame de crudo
En cada uno de los estados se involucrarían las siguientes válvulas que implicarían realizar una
revisión total del sistema y parar la extracción
Estado Válvulas Observación
X V1-V10 Todas las válvulas funcionando
Y V2, V3, V4, V5, V6, V8 Algunas válvulas de la 2 a la 6 o la 8 fallando
Z V1, V7, V9, V10 Cualquiera de las válvulas 1, 7, 9 o 10 fallando
Definiendo las condiciones de transición
Estado Estado Condición
Inicial Final
X Y Falla de una o varias válvulas de la 2 a la 6 o la 8
X Z Falla de una o varias válvulas 1, 7, 9 o 10
Y X Reparación de una o varias válvulas de la 2 a la 6 o la 8
Z X Reparación de una o varias válvulas 1, 7, 9 o 10
De acuerdo con la tabla anterior el espacio de estado estaría definido por:
W(t) = {(X,Y), (X,Z), (Y,X), (Z,X)}
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
Segunda entrega:
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 solo 4 válvulas 1,8,9,10 de seguridad como se muestra en la
figura
VARIABL
E Xn=((A(t)B(t)C(t)D(t)))
B si la válvula 1 está funcionando
A(t)=
M si la válvula 1 NO está funcionando
B si la válvula 8 está funcionando
B(t)=
M si la válvula 8 NO está funcionando
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
B si la válvula 9 está funcionando
C(t)=
M si la válvula 9 NO está funcionando
B si la válvula 10 está funcionando
D(t)= M si la válvula 10 NO está
funcionando
POSIBLES
HALLAZGO VALVUL VALVUL VALVUL VALVUL
S A1 A8 A9 A 10
Sin
1 B B B B derrame
Sin
2 B B B M derrame
Sin
3 B B M B derrame
Sin
4 B B M M derrame
Sin
5 B M B B derrame
Sin
6 B M B M derrame
Sin
7 B M M B derrame
Sin
8 B M M M derrame
9 M B B B derrame
10 M B B M derrame
11 M B M B derrame
12 M B M M derrame
13 M M B B derrame
14 M M B M derrame
15 M M M B derrame
16 M M M M derrame
TIEMPO ENTRE TIEMPO DE REVISIÓN
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
FALLAS SE VALVULAS SE
DISTRIBUYE DISTRIBUYE
EXPONENCIAL EXPONENCIAL
VALVUL
A
µ VALVUL
(Fallas/día) A (Válvulas /día)
1 0,008 1 0,005
8 0,014 8 0,005
9 0,016 9 0,005
10 0,014 10 0,005
Teniendo las probabilidades de los eventos individuales que pueden ocurrir en el transcurso de
un día, que es el periodo que se escogió para observar el sistema, calculamos las probabilidades
de la matriz
Tiempo entre fallas se distribuye exponencial P(X<1)
Probabilidad de que la Válvula 1 falle en un tiempo menor a 1 día 0,008
Probabilidad de que la Válvula 8 falle en un tiempo menor a 1 día 0,014
Probabilidad de que la Válvula 9 falle en un tiempo menor a 1 día 0,016
Probabilidad de que la Válvula 10 falle en un tiempo menor a 1 día 0,014
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
Tiempo de revisión de válvula se distribuye exponencial P(X<1
)
Probabilidad de que la Válvula 1 se haya revisado y reparado en un tiempo menor a 1
día 0,005
Probabilidad de que la Válvula 8 se haya revisado y reparado en un tiempo menor a 1
día 0,005
Probabilidad de que la Válvula 9 se haya revisado y reparado en un tiempo menor a 1
día 0,005
Probabilidad de que la Válvula 1 se haya revisado y reparado en un tiempo menor a 1
dia 0,005
Las válvulas son independientes una de las otras por lo que la probabilidad de ocurrencia de los
eventos es independiente, multiplicamos las probabilidades de cada evento individual que este
involucrado hallando así la probabilidad de que ocurra un evento conjunto
PROBABILIDA
EVENTO D
E1: La válvula 1 falla en un tiempo menor a 1 dia 0,008
E2: La válvula 1 falla en un tiempo mayor a 1 dia 0,992
E3: La válvula 8 falla en un tiempo menor a 1 dia 0,014
E4: La válvula 8 falla en un tiempo mayor a 1 dia 0,986
E5: La válvula 9 falla en un tiempo menor a 1 dia 0,016
E6: La válvula 9 falla en un tiempo mayor a 1 dia 0,984
E7: La válvula 10 falla en un tiempo menor a 1 dia 0,014
E8: La válvula 10 falla en un tiempo mayor a 1 dia 0,986
E9: La válvula 1 se revisa repara en un tiempo menor a 1 0,005
E10: La válvula 1 se revisa y repara en un tiempo mayor a 1 0,995
E11: La válvula 8 se revisa y repara en un tiempo menor a 1 0,005
E12: La válvula 8 se revisa y repara en un tiempo mayor a 1 0,995
E13: La válvula 9 se revisa y repara en un tiempo menor a 1 0,005
E14: La válvula 9 se revisa y repara en un tiempo mayor a 1 0,995
E15: La válvula 10 se revisa y repara en un tiempo menor a 1 0,005
E16: La válvula 10 se revisa y repara en un tiempo mayor a 1 0,995
E17: Todo el sistema se revisa y repara en un tiempo menor 0,133
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
E18: La válvula 10 se revisa y repara en un tiempo mayor 0,867
Teniendo en cuenta estas probabilidades realizamos los cálculos, tomamos las probabilidades
con tres cifras decimales, obteniendo la siguiente matriz de probabilidad de un paso P:
(BBBB) significa que cuando se va a observar el sistema de válvulas, las 4 válvulas se
encuentran funcionando. De igual manera, el estado (BMBB) significa que cuando se va a
observar el sistema de las 4 válvulas, la válvula 1,9,10 están funcionando, pero la válvula 8 no
está funcionando por falla. A partir de esta información y según la relación entre las graficas se
puede decir que hay derrame de crudo en los estados que se resaltan en rojo, considerando que
las válvulas son independientes, si se llega a presentar un derrame, se hará una revisión total, en
donde todas las válvulas serán revisadas y el sistema se dejara como nuevo.
PROGRAMACIÓN ESTOCASTICA
PROYECTO GRUPAL Institución Universitaria Politécnico Gran colombiano
CONCLUSION Y RECOMENDACION
en donde se presentan fallas en las válvulas 5 y 6 tiene la mayor probabilidad de aparición dentro del
modelo con una probabilidad de 0.0314El sistema de protección hidráulico presenta una distribución
de válvulas errónea puesto que el subsistema donde se ubican las válvulas 1, 8, 9 y 10 tiene una ruta
directa a la superficie solo con la falla en la válvula 1 lo que hace a las demás inservibles para controlar
un derrame y aumenta la probabilidad que se presente ese evento. En relación al otro subsistema donde
aparecen las válvulas del 2 al 7 se observa que aparecen válvulas que cumplen una misma función
lo cual hace de estas válvulas innecesarias generando un costo adicional al sistema de protección.
El modelamiento de un problema y para el caso particular de una cadena de Márkov permite entender
el funcionamiento de un sistema y estar en capacidad de disminuir la complejidad de su solución
por medio de la eliminación de estados similares en comportamiento.
