Simulación geoestadística aplicada al modelamiento de yacimientos de petróleo
Autor: Castaño Agudelo, Andres Felipe
Asesor: Alzate Espinosa, Guillermo Arturo (Thesis advisor)
Fecha: 2004
Programa: Bogotá - Medicina -Fisioterapia
Facultad: Facultad de Medicina

Resumen: En ingenieria de petróleos la caracterización y el entendimiento de los parámetros de un yacimiento son las claves para poder realizar estudios que permitan tomar decisiones acerca de la forma más efectiva en la que este se debe explotar. Para cumplir con dicho objetivo es necesario obtener información de las características que definen al yacimiento a través de métodos como los registros eléctricos, las pruebas de presión, etc. Lástimosamente, por cuestiones técnicas y ecónomicas, no es posible tener esa información para todos los puntos del yacimiento. Por eso se hace necesario estudiar técnicas que permitan generar valores de propiedades del yacimiento en lugares donde no se ha realizado ningún tipo de muestreo, a partir de los datos que se han medido. Para el estudio de los datos medidos, que en particular siempre tienen asociados una posición particular en el espacio y son conocidos como variables regionalizadas, se usan actualmente módelos matemáticos dentro de los cuales la geoestadística ha jugado un papel importante en las últimas tres decadas y ha probado su superioridad entre muchos tipos de estimación. En este trabajo se presentan los elementos fundamentales de la geoestadística y en particular el de simulación condicional, la cual usa al kriging que es el estimador por excelencia usado en geoestadística. Después de presentar dichos conceptos, se muestra un ejemplo de aplicación para obtener valores de porosidad y permeabilidad, y se usan los modelos obtenidos para plantear la forma de realizar un análisis económico de un proyecto de inyección de agua. / Abstract: In Petroleum Engineering, the characterization and understanding of reservoir parameters are the keys to develop studies that let to make decisions about the most effective way the reservoir must be produced. In order to reach such goal it is necessary to obtain information about the characteristics that define the reservoir using methods like the open hole electric logs, pressure well testing, etc. Pitifully, by technical and economical issues, it is not possible to obtain such information in all the points that cover the reservoir. That is why it is necessary to study techniques that let generate values of the reservoir properties in places where there is not any measurement, starting from the measure data. For the study of the measure data, that in particular are always associated with a position in the space and are knowing as regional variables, mathematical models are used nowadays, in particular the geoestatistics has played an important role in the last three decades and it has probe its superiority among many kinds of estimation. In this work, fundamental elements of geostatistics are presented. In particular conditional simulation, that uses kriging which is the essential estimator used by geostatistics. After presenting such concepts, an application example is showed where Porosity and Permeability values are obtained, and the final models are used to suggest the way to realize and economical analysis in a Water flooding project.

SIMULACIÓN GEOESTADÍSTICA APLICADA AL MODELAMIENTO DE
YACIMIENTOS DE PETRÓLEO
ANDRÉS FELIPE CASTAÑO AGUDELO
FRANCISCO VERGARA ELORZA
UNIVERSIDAD NACIONAL DE COLOMBIA
SEDE MEDELLÍN
FACULTAD DE MINAS
2004SIMULACIÓN GEOESTADÍSTICA APLICADA AL MODELAMIENTO DE
YACIMIENTOS DE PETRÓLEO
ANDRÉS FELIPE CASTAÑO AGUDELO
FRANCISCO VERGARA ELORZA
Trabajo dirigido de grado para optar al título
de ingeniero de petróleos
Asesor
GUILLERMO ALZATE
Ingeniero de Petróleos
UNIVERSIDAD NACIONAL DE COLOMBIA
SEDE MEDELLÍN
FACULTAD DE MINAS
2004AGRADECIMIENTOS
Queremos agradecer a todas las personas que de una o otra manera hicieron posible la
realización de este trabajo, muy en particular al Profesor Guillermo Arturo Alzate (M.Sc
Ingeniería de Petróleos) por haber servido como asesor del proyecto y más aún por su
incondicional ayuda que sin duda alguna será muy importante para nosotros poder
cumplir con las nuevas metas que nos hemos trazado por conseguir.
También agradecer a nuestros compañeros de la Universidad que nos ayudaron a lo largo
de estos años. En particular aquellos que nos sirvieron de puente de comunicación con la
Universidad cuando no pudimos estar presentes. Les hacemos partícipes de este logro:
Jaime, Pablo, Juan Guillermo, Carlitos, David (Yeti), Isabel, Marcelo, Eduard, y a los
demás. Su apoyo es invaluable.
Al profesor Dr. Richard Hughes de la Universidad de Oklahoma quien de muy buena
manera nos facilitó la información para el ejemplo de aplicación y sirvió de medio para
poder acceder al uso del programa 3DSL para simulación de flujo tipo streamline.
Andrés Felipe Castaño Agudelo
Francisco Vergara ElorzaRESUMEN
En ingenieria de petróleos la caracterización y el entendimiento de los parámetros de un
yacimiento son las claves para poder realizar estudios que permitan tomar decisiones acerca
de la forma más efectiva en la que este se debe explotar. Para cumplir con dicho objetivo es
necesario obtener información de las características que definen al yacimiento a través de
métodos como los registros eléctricos, las pruebas de presión, etc. Lástimosamente, por
cuestiones técnicas y ecónomicas, no es posible tener esa información para todos los puntos
del yacimiento. Por eso se hace necesario estudiar técnicas que permitan generar valores de
propiedades del yacimiento en lugares donde no se ha realizado ningún tipo de muestreo, a
partir de los datos que se han medido. Para el estudio de los datos medidos, que en
particular siempre tienen asociados una posición particular en el espacio y son conocidos
como variables regionalizadas, se usan actualmente módelos matemáticos dentro de los
cuales la geoestadística ha jugado un papel importante en las últimas tres decadas y ha
probado su superioridad entre muchos tipos de estimación.
En este trabajo se presentan los elementos fundamentales de la geoestadística y en
particular el de simulación condicional, la cual usa al kriging que es el estimador por
excelencia usado en geoestadística. Después de presentar dichos conceptos, se muestra un
ejemplo de aplicación para obtener valores de Porosidad y Permeabilidad, y se usan los
modelos obtenidos para plantear la forma de realizar un análisis económico de un proyecto
de inyección de agua.
Palabras clave: geoestadística, semivariograma, anisotropía geoestadística, estimación,
kriging, simulación geoestadística, simulación gaussiana secuencial, simulación stream-
line, distribución de probabilidad.ABSTRACT
In Petroleum Engineering, the characterization and understanding of reservoir parameters
are the keys to develop studies that let to make decisions about the most effective way the
reservoir must be produced. In order to reach such goal it is necessary to obtain information
about the characteristics that define the reservoir using methods like the open hole electric
logs, pressure well testing, etc. Pitifully, by technical and economical issues, it is not
possible to obtain such information in all the points that cover the reservoir. That is why it
is necessary to study techniques that let generate values of the reservoir properties in places
where there is not any measurement, starting from the measure data. For the study of the
measure data, that in particular are always associated with a position in the space and are
knowing as regional variables, mathematical models are used nowadays, in particular the
geoestatistics has played an important role in the last three decades and it has probe its
superiority among many kinds of estimation.
In this work, fundamental elements of geostatistics are presented. In particular conditional
simulation, that uses kriging which is the essential estimator used by geostatistics. After
presenting such concepts, an application example is showed where Porosity and
Permeability values are obtained, and the final models are used to suggest the way to
realize and economical analysis in a Water flooding project.
Key Words: Geostatistics, semivariogram, geostatistical anisotropy, estimation, kriging,
geostatistical simulation, sequential gaussian simulation, streamline simulation, probability
distribution.TABLA DE CONTENIDO
Pág
LISTA DE TABLAS I
LISTA DE FIGURAS II
LISTA DE SÍMBOLOS Y VARIABLES III
RESUMEN IV
ABSTRACT V
1. INTRODUCCIÓN 1
2. FUNDAMENTO TEÓRICO DE LA GEOESTADÍSTICA. 3
2.1. GENERALIDADES 3
2.2. EL SEMIVARIOGRAMA 9
2.2.1. Modelamiento del Semivariograma. 12
2.2.1.1. Modelo Esférico. 14
2.2.1.2. Modelo Exponencial. 15
2.2.1.3. Modelo Gaussiano. 15
2.2.1.4. Modelos Combiandos. 17
2.2.1.5. Modelamiento de Efectos de Oscilación. 18
2.2.1.6. Otros Modelos de Variograma. 20
2.2.2. Admisibilidad Para Modelos de Variograma. 21
2.3. CONCEPTO DE ANISOTROPÍA GEOESTADÍSTICA. 22
2.4. UN EJEMPLO SIMPLE. 25
2.5. MODELOS DE ESTIMACIÓN. 27
2.5.1. Estimación Lineal – Kriging. 29
2.5.1.1. El Modelo de Función Aleatoria y el Concepto de Estimador Insesgado. 29
2.5.1.2. El Modelo de Función Aleatoria y la Varianza del Error. 31
2.5.1.3. El Parámetro de Lagrange. 32
2.5.1.4. Minimización de la Varianza del Error. 332.5.2. Corregionalización. 36
2.5.2.1. Modelamiento de una Corregionalización. 37
2.5.3. Kriging con Parámetro Externo. 37
2.5.4. Kriging Colocado. 38
2.6. SIMULACIÓN GEOESTADÍSTICA. 39
2.6.1. Métodos de Simulación Condicional. 39
2.6.1.1. Métodos Secuenciales. 39
2.6.1.2. Simulación Gaussiana Secuencial. 40
2.7. FUNCIONES DE PÉRDIDA O ERROR. 41
2.7.1 Decisiones Basadas en Funciones de Pérdida. 42
3. DEFINICIÓN DEL PROBLEMA. 45
4. METODOLOGÍA DE TRABAJO. 48
5. RESULTADOS. 51
5.1. MAPAS DE VARIOGRAMA. 51
5.2. VARIOGRAFÍA. 53
5.3. KRIGING. 56
5.4. SIMULACIÓN GEOESTADÍSTICA. 58
5.5. SIMULACIÓN DE FLUJO TIPO “STREAM-LINE”. 66
6. ANÁLISIS DE RESULTADOS. 70
6.1. SUPERFICIES VARIOGRÁFICAS. 70
6.2. VARIOGRAFÍA. 70
6.3. KRIGING Y SIMULACIÓN GEOESTADÍSTICA. 71
6.4. APLICACIONES. 73
6.4.1. Evaluación de un Proyecto de Ïnfill-Drilling”. 74
7. CONCLUSIONES 77
8. RECOMENDACIONES 79
REFERENCIAS.
BIBLIOGRAFÍA.
ANEXOSLISTA DE TABLAS
Pág
TABLA 1. Diferentes estadísticos de Dispersión 6
TABLA 2. Parámetreos de las ecuaciones de los modelos de Variograma 57
TABLA 3. Propiedades de los Fluidos Utilizadas para la Simulación de Flujo 67
TABLA 4. Estadísticos de la Distribución de Probabilidad del Petróleo Final Recuperado 69LISTA DE FIGURAS
Pág
FIGURA 1. Diagramas H de Dispersión. 5
FIGURA 2. Coeficiente de Correlación, Covarianza y Momento de Inercia 7
FIGURA 3. Área de Busqueda de Datos. 10
FIGURA 4. Principales Parámetros en un Semivariograma. 11
FIGURA 5. Modelo Esférico 14
FIGURA 6. Modelo Exponencial 15
FIGURA 7. Comparación del Modelo Esférico y Exponencial 16
FIGURA 8. Modelo Gaussiano 16
FIGURA 9. Modelo Combinado 18
FIGURA 10. Modelo Cosenoidal para Fenómenos de Oscilación 19
FIGURA 11. Modelo Sinusoidal para Fenómenos de Oscilación 20
FIGURA 12. Modelo Cosenoidal para Fenómenos Gradualmente Atenuados 20
FIGURA 13. Vista Isométrica de una Anisotropía 22
FIGURA 14. Reducción de una Anisotropía Geométrica 23
FIGURA 15. Anisotropía Geométrica 25
FIGURA 16. Semivariograma Omnidireccional 26
FIGURA 17. Ajuste Esférico para el Semivariograma Omnidireccional 26
FIGURA 18. Ajuste Exponencial para el Semivariograma Omnidireccional 27
FIGURA 19. Histograma de Petróleo Recuperado 43
FIGURA 20. Función de Frecuencia Acumulada de Petróleo Recuperado 43
FIGURA 21. Función de Pérdida 44
FIGURA 22. Minimización de la Función de Pérdida 44
FIGURA 23. Ubicación de los Pozos en el Campo 45
FIGURA 24. Diagrama de Flujo: Kriging Ordinario 48
FIGURA 25. Diagrama de Flujo: Simulación Geoestadística 50
FIGURA 26. Histograma de Frecuencias de Porosidad 51FIGURA 27. Histograma de Frecuencias de Log-Permeabilidad 52
FIGURA 28. Mapa de Variograma 2-D para Porosidad 53
FIGURA 29. Mapa de Variograma 3-D para Porosidad 54
FIGURA 30. Mapa de Variograma 2-D para Permeabilidad 54
FIGURA 31. Mapa de Variograma 3-D para Permeabilidad 55
FIGURA 32. Modelo Ajustado del Variograma de Porosidad 56
FIGURA 33. Modelo Ajustado del Variograma de Permeabilidad 56
FIGURA 34. Mapa de Porosidad Elaborado con Kriging 58
FIGURA 35. Mapa de Permeabilidad Elaborado con Kriging 59
FIGURA 36. Simulaciones Geoestadísticas de Porosidad y Permeabilidad (19479) 61
FIGURA 37. Simulaciones Geoestadísticas de Porosidad y Permeabilidad (37333) 62
FIGURA 38. Simulaciones Geoestadísticas de Porosidad y Permeabilidad (67495) 63
FIGURA 39. Simulaciones Geoestadísticas de Porosidad y Permeabilidad (99999) 64
FIGURA 40. Simulaciones Geoestadísticas de Porosidad y Permeabilidad (123451) 65
FIGURA 41. Arreglo de Pozos Para Inyección de Agua 66
FIGURA 42. Histograma del Petróleo Final Recuperado a 15 Años 68
FIGURA 43. Distribución de Frecuencia Acumulada del Petróleo Final Recuperado 69
FIGURA 44. Mapa de Localización de Datos de Porosidad 71
FIGURA 45. Mapa de Localización de Datos de Permeabilidad 72
FIGURA 46. Árbol de Decisiones para un Proyecto de “Infill-Drilling” 76LISTA DE SIMBOLOS Y VARIABLES
h: vector de distancia entre un par de datos
ρ: coeficiente de correlación
ρ(h): función de correlación o correlograma
C(h): función de covarianza
γ(h): semivariograma
σ: desviación estándar
m: media de una distribucion
N(h): número de pares separados por una distancia h
θ: ángulo del vector h
C : valor del efecto nugget
0
a: rango del variograma
C: sill del variograma
Z(x): función aleatoria estacionaria
λ: pesos de una combinación lineal
i
Ψ: combinación lineal de pesos λ
i
r: error de estimción
E{V}: valor esperado
R(X ): error en el punto X
o o
µ: parámetro de lagrange
Z: valor medido de una variable
Z*: valor estimado de una variable1. INTRODUCCIÓN
En el campo de las geociencias es muy frecuente encontrar variables que se manifiestan
distribuidas espacialmente. Para el estudio de estas variables son usados actualmente
diversos procedimientos de estimación y simulación. Esto es, a partir de un conjunto de
muestras tomadas en localizaciones del domino en que se manifiesta un fenómeno y
consideradas representativas de su realidad, que por lo general es siempre desconocida,
estos procedimientos permiten su descripción o caracterización con dos fines diferentes: El
primero es proporcionar valores estimados en localizaciones de interés y el segundo,
generar valores que en conjunto presenten iguales características de dispersión y
estadísticas que los datos originales.
La estimación requiere un modelo de como se comporta el fenómeno que se va a estudiar
en los lugares donde este no se ha muestreado: sin un modelo no se podrían realizar
inferencias acerca de los valores desconocidos de la variable regionalizada. En la ingeniería
de yacimientos y tal vez en todo estudio científico el factor incertidumbre nunca dejará de
tener peso en el éxito o fracaso de las decisiones tomadas. Esta incertidumbre existe debido
a la falta de conocimiento de lo que realmente yace en el subsuelo. Para el ingeniero de
yacimientos sería valioso tener un alto grado de certeza en las características y propiedades
de la formación productora a la hora de hacer una estimación de reservas o cuando se
realiza una predicción de la producción, pues de esta forma la viabilidad económica de un
proyecto sería fácilmente determinada y las consecuencias por sobrestimar o subestimar el
potencial del yacimiento serían mínimas. Solo por citar otro ejemplo, considérese los
sofisticados simuladores de flujo disponibles a la fecha. Uno de estos poderosos programas
sería completamente inútil si el modelo de los parámetros del yacimiento que alimenta al
simulador no está fundado en información de buena calidad. Esto quiere decir que el peso
de la incertidumbre parece ser controlable en tanto que la calidad de la información
recopilada sea llevada al máximo.
Particularmente en la geología y en la ingeniería de yacimientos, existe la tendencia a crear
algoritmos que generan una imagen suavizada de la realidad estudiada que solo revela
tendencias a gran escala, sin tener en cuenta las fluctuaciones siempre presentes en la
práctica. Pero el problema a la hora de estudiar un campo petrolífero no se limita solo a
comprender las características geológicas del mismo debido a que el flujo de fluidos en el
yacimiento es también de vital importancia y bien es sabido que la respuesta del flujo es
enormemente afectada por la distribución y el contraste de altos y bajos valores de
permeabilidad, por ejemplo. Es entonces en este punto donde la simulación geoestadística
aparece atractiva para el ingeniero de yacimientos ya que esta genera múltiples "escenarios
geoestadísticos" de una variable que reproducen la misma "realidad estadística" de ella
minimizando al tiempo el mapa de varianzas de los valores estimados. De esta forma la
respuesta al flujo de tales escenarios se podría acercar un poco más a la realidad. Porejemplo, de entre muchos otros beneficios, un ajuste histórico de la información disponible
sería más probable de obtenerse cuando se dispone de los modelos de yacimiento generados
por simulación geoestadística.
En este trabajo se presentan los elementos fundamentales de la teoría de simulación
condicional, el cual usa en el proceso de condicionamiento al kriging que es el
procedimiento para estimación por excelencia que propone la geoestadística. Además se
presenta un ejemplo de aplicación en donde se generan modelos geoestadísticos de las
variables porosidad y permeabilidad de un campo de petróleo, a partir de la información
disponible de dichas variables. La intención del trabajo no es encontrar el verdadero
modelo de las propiedades de yacimiento, sino enfatizar en las aplicaciones que se pueden
lograr después de realizar el estudio geoestadístico del mismo. Es así como al final del
trabajo se encontrará una aplicación en la que se muestra como se pueden encontrar
distribuciones de probabilidad de variables de interés en ingeniería de yacimientos a la hora
de analizar económicamente un proyecto de explotación de hidrocarburos.2. FUNDAMENTO TEÓRICO DE LA GEOESTADISTICA
2.1 GENERALIDADES
El lector de un tema cualquiera en geoestadística debe conocer y estar familiarizado con los
términos y conceptos básicos empleados en este tipo de estudio. Algunos de estos términos
pueden parecer extraños en el contexto de la ingeniería de petróleos pero por consistencia
con la literatura serán mantenidos aquí. En esta parte del capítulo se presentan algunos
conceptos básicos para poder entender la teoría geoestadísitica.
Una variable distribuida en un espacio o área es una variable regionalizada. Ejemplos de
estas variables son el espesor de una formación, el gradiente geotérmico, la producción
acumulada, la porosidad y la permeabilidad. El espacio en el cual se distribuye la variable
puede incluso considerar el tiempo. Por ejemplo, los valores de temperatura a través del
tiempo en un punto dado es una variable regionalizada.
Una variable aleatoria toma valores que siguen una distribución cualquiera de
probabilidad. Considérese como ejemplo un valor de permeabilidad medido en un
horizonte dado de un pozo. Entonces dicho valor de permeabilidad es una realización
particular de una variable aleatoria. Cada valor de permeabilidad observado representa
una variable aleatoria diferente y el conjunto formado por estas variables constituye una
función aleatoria. Para estimar un nuevo valor de esta función en un punto nuevo
cualquiera, se usan los valores ya observados en conjunto con algún tipo de información
extra que se posea del problema.
Dicha función aleatoria tiene dos componentes:
• Un componente regional, que manifiesta cierto grado de autocorrelación espacial y
dependencia entre valores próximos.
• Un componente aleatorio local.
Los geólogos usualmente asumen la presencia del primer componente mientras desprecian
el segundo. Por el contrario, la geoestadística cuantifica la contribución relativa de ambos
componentes.
Uno de los aspectos característicos de los datos utilizados en un estudio geoestadístico es
que estos están siempre localizados espacialmente. Tal característica le añade al conjunto
de datos ciertas propiedades particulares como lo son: localización de valores extremos,
tendencias a través del espacio y un grado de continuidad espacial. Pocas o ninguna de
las técnicas aportadas por la estadística univariada o bivariada pueden capturar las
3propiedades particulares mencionadas. A continuación se presentan los conceptos básicos
que permitirán la introducción de la localización de una variable en el tratamiento de un
grupo de datos recolectados.
La continuidad espacial es de esperarse cuando se trabaja con datos localizados
espacialmente. Aquellos valores situados a escasa distancia tendrán mayor probabilidad de
ser similares en valor que aquellos que están distantes el uno del otro. Cuando se miran
datos en un sistema coordenado, los valores no parecen estar ubicados aleatoriamente, más
bien, los valores bajos tienden a estar cerca de valores igualmente bajos y los valores altos
tienden a estar cercanos a valores altos.
Las herramientas usadas para describir la relación entre dos variables también pueden ser
usadas para describir la relación de una variable con valores de esta misma en puntos
cercanos. Una de tales herramientas son los diagramas H de dispersión.
Un diagrama H de dispersión muestra todos los pares de datos posibles cuyas ubicaciones
están separadas por una cierta distancia en una dirección particular. La coordenada en x
corresponde a un valor ocurrido en un punto particular y la coordenada en y corresponde al
valor ocurrido a una distancia h de dicho punto. Por ejemplo, cuando h = (0,1) significa
que se ha tomado cada punto y se ha apareado con aquellos puntos cuya coordenada Este es
la misma y su coordenada Norte es una unidad mayor. Dicho de otra forma, debido a que h
es un vector, este indicará una distancia y una dirección. Para este ejemplo, h = (0,1), la
distancia indicada es 1 unidad y la dirección es 90 grados.
La forma de la nube de puntos en un diagrama H de dispersión revela la continuidad
espacial que presentan los datos a una distancia y dirección dada por el vector h. Si los
valores ubicados en puntos separados por una distancia h son muy similares, entonces la
nube de puntos se distribuirá cercana a la línea de pendiente 45E que pasa por el origen.
Contrariamente, si los valores se hacen menos similares, la nube de puntos lucirá mucho
más dispersa.
En la Figura 1 se hace evidente el cambio de continuidad en la dirección Norte al realizar
varios diagramas H de dispersión para una variable cualquiera. Los cuatro diagramas, que
van desde una distancia de 1 m hasta 4 m en la dirección Norte, dan evidencia de cómo la
nube de puntos se hace más dispersa a medida que la distancia de los puntos apareados
aumenta. Aunque la línea de 45E pasa a través de la nube de puntos en cada uno de los
diagramas, la nube no se muestra simétrica respecto a dicha línea.
Nótese que en los cuatro diagramas de la Figura 1, algunos puntos se alejan mucho más de
los demás que forman la nube. Por ejemplo, en cada uno de los diagramas existe un punto
ubicado cerca a la línea del eje x. Puntos como este sugieren un comportamiento inusual
haciendo que sea buena práctica chequear como luciría el mismo análisis si puntos como
este son removidos.
4Como en todo análisis estadístico, es necesario un resumen cuantitativo de la información
contenida en los diagramas H. Una de las características esenciales de este tipo de diagrama
es la dispersión de la nube de puntos obtenida. Esta característica se puede cuantificar a
partir del coeficiente de correlación r . A medida que la nube de puntos de hace más
dispersa, el coeficiente de correlación tiende a decrecer. El coeficiente de correlación para
cada uno de los diagramas de la Figura 1 se muestra en la Tabla 1. Como se esperaba, el
coeficiente de correlación decrece al aumentar la dispersión de la nube de puntos, lo cual
significa que es un buen índice para cuantificar la característica de dispersión.
FIGURA 1. Diagramas H de dispersión. (a) h = (0,1), (b) h = (0,2), © h = (0,3), (d) h = (0,4)
Tomado de Srivastava e Isaaks, “Applied Geostatistics”, 1989
La relación entre el coeficiente de correlación de un diagrama H de dispersión y el valor de
h para ese diagrama, se denomina la función de correlación o correlograma. El
coeficiente de correlación depende de h el cual, siendo un vector, tiene una magnitud y una
dirección. Una gráfica de la función de correlación, muestra el coeficiente de correlación
como una función de la magnitud y dirección de h. Usualmente se hacen gráficas separadas
del coeficiente de correlación versus la magnitud de h para varias direcciones.
La Figura 2 parte (a) muestra como el coeficiente de correlación decrece con el aumento de
5la distancia en la dirección norte. Gráficas similares en las demás direcciones nos darían
una idea de cómo el coeficiente de correlación cambia en función de la distancia y la
dirección.
TABLA 1. Diferentes estadísticos de dispersión
Tomado de Srivastava e Isaaks, “Applied Geostatistics”, 1989
Coeficiente de
H Covarianza Momento de Inercia
Correlación
0,1 0.742 448.8 312.8
0,2 0.59 341 479.2
0,3 0.56 323.8 521.4
0,4 0.478 291.5 652.9
Un índice alternativo de la continuidad espacial es la covarianza. La covarianza para cada
uno de los diagramas en la Figura 1 está dada en la Tabla 1. Estos valores muestran como
también la covarianza decrece en forma similar al coeficiente de correlación. La relación
entre la covarianza de un diagrama H de dispersión y el valor de h para ese diagrama es
llamada la función de covarianza C(h). La Figura 2 parte (b) muestra la función de
covarianza en la dirección norte.
Otro índice para cuantificar la dispersión de la nube de puntos es el momento de inercia M
respecto a la línea de 45E que pasa por el origen. Dicho momento de inercia puede ser
calculado de la siguiente ecuación:
1 ∑n ( )2
M = x - y
2n i i (1)
i=1
El momento de inercia calculado con la ecuación (1), es la mitad del promedio del cuadrado
de la diferencia entre la coordenada x y la coordenada y para cada par de puntos en el
diagrama H de dispersión. El factor ½ es consecuencia de que se esta calculando la
distancia perpendicular de los puntos a la línea de 45E .
La relevancia del momento de inercia respecto a la línea de 45E radica en el hecho de estar
apareando puntos de la misma variable. Considérese el caso de aparear cada punto consigo
mismo, es decir, cuando h = (0,0). En tal caso, todos los puntos del diagrama deberían
coincidir con la línea de 45E grados. Entonces el significado de dispersión de los puntos a
medida que h se incrementa es que se tiene otro índice de continuidad a partir del momento
de inercia.
En forma contraria a los demás índices de continuidad ya expuestos, el momento de inercia
6crece a medida que la nube de puntos se hace más dispersa. Lo anterior se puede observar
en la Tabla 1, donde el momento de inercia crece al tiempo que el coeficiente de
correlación y la covarianza decrecen.
La relación entre el momento de inercia de un diagrama H de dispersión y su valor de h es
llamado semivariograma γ(h). El semivariograma, o simplemente el variograma, en la
dirección norte se muestra en la Figura 2 parte (c).
.
Los tres métodos estadísticos propuestos para cuantificar la dispersión de los puntos en la
nube son sensibles a datos erróneos. Siempre es importante la evaluación de un cambio en
la estadística de los datos cuando tales puntos son removidos. En la práctica, tanto la
función de correlación como la función de covarianza y el variograma podrían no describir
satisfactoriamente la continuidad espacial debido a unos pocos datos erróneos. Nuevamente
se recomienda como práctica, chequear la existencia de datos erróneos que puedan alterar el
estudio de la continuidad espacial.
FIGURA 2. Figuras de (a) coeficiente de correlación, (b) covarianza y (c) momento de inercia.
Tomado de Srivastava e Isaaks, “Applied Geostatistics”, 1989
Aunque los diagramas H de dispersión contienen mucha más información que cualquiera
de los tres estadísticos propuestos, es una práctica común evitar el cálculo de tales
diagramas e ir directamente al cálculo de r (h), C(h) o γ(h) para describir la continuidad
espacial. Es conveniente, por lo tanto, tener las fórmulas de estas funciones expresadas en Comentario [Comment1]: fu
términos de los valores de los datos recolectados y no en términos de las coordenadas de ncion
cada valor.
7La función de covarianza, C(h), puede ser calculada de:
C( h)= N1 ( h) ∑ v i (cid:215) v j- m -(cid:215) h m + h (2)
(i,j)hi,j=h
Donde n ,.....,n representan los datos; la sumatoria se aplica sobre los N(h) pares cuyas
1 n
ubicaciones están separadas por una distancia h. m es la media de todos los datos cuyas
-h
ubicaciones están separadas por una distancia -h, mientras que m es la media de todos los
+h
datos cuyas ubicaciones están separadas por una distancia +h.
La función de correlación, r (h), es la función de covarianza estandarizada por las
desviaciones estándar de todos los datos cuyas ubicaciones están separadas por una
distancia -h y +h.
C(h)
r(h)=
s s(cid:215) (3)
- h + h
El semivariograma, γ(h), es la mitad del promedio del cuadrado de la diferencia en
distancia de dos puntos apareados:
1 ∑( )2
g(h) = v - v
2N(h) i j (4)
(i,j)hij=h
Los valores de r (h), C(h) y γ(h) no se ven afectados si los subíndices i y j son
intercambiados en las ecuaciones anteriores. Esto implica que:
r (h) = r (-h) , C(h) = C(-h) y g (h) = g (-h) (5)
Es decir, el valor calculado en una dirección cualquiera de las tres funciones anteriores, será
idéntico al calculado en la dirección exactamente opuesta. Es por esto que cuando se
describe la continuidad espacial se habla en direcciones combinadas, norte - sur por
ejemplo, en lugar de norte solamente.
Ahora bien, aunque una serie de diagramas H provee la más completa descripción espacial,
usualmente estos contienen demasiada información que necesita cierto tipo de resumen. Ya
se han introducido tres funciones que ayudan en este propósito: la función de correlación,
r (h); la función de covarianza, C(h); y el variograma ,γ(h). Todas ellas utilizan
información de los diagramas H de dispersión para describir como la continuidad espacial
de un grupo de datos cambia en función de la distancia y la dirección. Así pues, cualquiera
de ellos es adecuado para el propósito descriptivo, sin embargo, el variograma es la función
8tradicionalmente usada.
2.2 EL SEMIVARIOGRAMA
El primer paso en un análisis geoestadístico es la variografía. Esto implica calcular y
modelar el variograma. El semivariograma es la herramienta básica de la geoestadística en
la medición de la autocorrelación espacial de una variable regionalizada. Como su nombre
implica, un semivariograma es un método de medir varianza. Si bien existen
procedimientos para modelar un semivariograma a través de métodos iterativos de mínimos
cuadrados, en la práctica se recomienda una detallada observación del modelo ajustado al
comportamiento exhibido por los datos. Finalmente, un modelo adecuadamente ajustado
permitirá el cálculo de estimativos que reflejarán la correlación espacial de la variable
tratada.
Una vía sencilla para cuantificar como una variable z cambia en valor entre un sitio x y otro
situado a una distancia de h unidades, es decir, a una distancia x + h, es calcular la
diferencia z(x) - z(x + h). Si la superficie representada por los dos puntos es continua y h es
una distancia corta, se esperaría un valor pequeño para la diferencia en cuestión.
Contrariamente, a medida que h se incrementa, dicha diferencia se hace mayor.
Trasladando esta idea intuitiva a una fórmula, es obvio que se estará interesado en el
comportamiento de:
∑[ z(x)- z(+x h)]2
g( h)= (6)
2n
Al igual que la bien conocida varianza de la estadística básica, el semivariograma es una
suma de cuadrados dividida por el numero n de diferencias en la muestra. A diferencia de
la varianza de la respecto a una media usada en estadística clásica, el semivariograma mide
la diferencia de un dato respecto a otro.
Cuando una variable es medida a intervalos regulares dentro de un trayecto, se puede
calcular el valor del semivariograma para valores enteros de h. Pero bien sabemos que en
la ingeniería de petróleos y la geología, muestras tomadas en intervalos irregularmente
espaciados es lo normal. Así que el primer problema que se tiene cuando se trata de
construir el semivariograma es que para cualquier valor de h escogido se tendrán muy
pocos datos que estén separados exactamente por una distancia h. De hecho, este problema
ocurre igualmente cuando se trata de definir una dirección para el semivariograma. Pero ya
que esta es la realidad con la que se debe tratar, podría pensarse en una solución práctica
como la siguiente: Aceptar cualquier par de datos cuya separación y dirección es más o
menos cercana a h y q , respectivamente. Es decir, se ha introducido un valor de tolerancia
y así cualquier par de puntos que se encuentren en el área sombreada que se muestra en la
Figura 3 se asume que están exactamente a una distancia h y a una dirección q .
9El primer paso en la construcción del variograma en una dirección particular es la selección
del incremento para la distancia h y su respectivo valor de tolerancia. Mientras más
pequeño el incremento, se tendrá la posibilidad de calcular mayor cantidad de puntos para
una mejor definición del variograma. Pero al mismo tiempo esto implicará un menor
número de pares al momento de calcular cada punto. Obviamente, el número de
incrementos está limitado por la extensión del área bajo estudio. Respecto al valor de
tolerancia, es una práctica común el asignarle una distancia constante. Sin embargo, una
alternativa podría ser el uso de una tolerancia pequeña para valores de h pequeños y
conforme el valor de h incrementa, el valor de la tolerancia puede incrementarse también.
En la práctica, un buen modelo de variograma es el resultado de un trabajo iterativo en el
cual se intenta mantener un equilibrio entre la obtención de una curva suave a partir de un
tamaño de muestra grande y el uso de un buen número de puntos que describan el
comportamiento de la curva especialmente en el origen.
FIGURA 3. Ejemplo del área de búsqueda de datos después de ser asignados los respectivos
valores de tolerancia.
Tomado de Srivastava e Isaaks, “Applied Geostatistics”, 1989
Ahora es apropiado introducir la terminología que es usada para describir las características
relevantes de un variograma (Figura 4).
Rango: A medida que la distancia de separación entre los puntos apareados aumenta, el
correspondiente valor del variograma generalmente también se incrementa.
Eventualmente, los incrementos en h no causarán más incremento en el valor del
variograma, de tal forma que este alcanzará una zona de valores constante sin importar los
cambios en h. La distancia a la cual el variograma alcanza esta zona es llamada rango del
10variograma.
Sill: corresponde al valor que alcanza el variograma cuando este llega al rango.
Efecto Nugget: aunque el valor del variograma cuando h = 0 es estrictamente 0, varios
factores, tales como el error al momento de tomar los datos y la variabilidad a pequeña
escala, pueden causar que valores separados por distancias muy cortas sean
extremadamente diferentes. Esto causa una discontinuidad en el origen del variograma.
Este salto vertical desde el origen hasta el valor exhibido por el variograma es llamado el
efecto nugget. La razón entre el efecto nugget y el sill es conocida como el efecto nugget
relativo y es usualmente reportado en porcentaje.
La presencia del efecto nugget implica, en la variable estudiada, abruptos cambios de
continuidad en pequeñas distancias, alta variabilidad a una escala más pequeña que la de
muestreo o baja precisión en la medición. La presencia o ausencia de efecto nugget puede
entonces depender de la escala relativa de muestreo respecto a la escala de variación del
fenómeno estudiado.
FIGURA 4. Principales parámetros en un semivariograma.
Considérese como ejemplo el caso de un yacimiento naturalmente fracturado1. Pozos
adyacentes pueden exhibir comportamientos muy diferentes dependiendo de su proximidad
relativa a las fracturas naturales. El rango de influencia de una fractura sobre el variograma
depende de la probabilidad de que un pozo intercepte a las fracturas y tal probabilidad
11depende de varios factores como lo son el espesor de la formación o el grado de
verticalidad de la fractura. Por ejemplo, la creación de fracturas a través de tecnología de
completamiento incrementará la probabilidad de que un pozo se intercomunique con las
fracturas naturales. De esta forma, el grado de influencia de una fractura natural puede
verse incrementado. Si el yacimiento es muy fracturado en relación al espaciamiento de
pozos, la variabilidad debida al fracturamiento no tendrá una influencia detectable.
Igualmente, si los pozos están tan cerca como sea físicamente posible, la probabilidad de
interceptar una fractura natural aumentará y así el efecto nugget debido a fracturas
disminuirá.
La discusión anterior sugiere que el efecto nugget no es más que un nuevo modelo de
variograma con un rango muy pequeño que podría usarse para modelar el aporte de la suma
de todas las variaciones cuyo rango es más pequeño que las distancias de muestreo. Sin
embargo, este concepto no resulta muy práctico pues es difícil contar con un muestreo
adicional a pequeña escala que permita aplicarlo.
Cuando un semivariograma resulta ser completamente plano o simplemente varia
aleatoriamente alrededor de una línea, se dice que el grupo de datos presenta un efecto
nugget puro, es decir, la muestra es espacialmente independiente.
Finalmente puede decirse que el efecto nugget representa el componente aleatorio de una
variable regionalizada. Típicamente es llamado ruido y es modelado con un único número
C , el cual representa el sill de un variograma cuyo rango es muy pequeño. Si al momento
0
de analizar un determinado grupo de datos se puede estar seguro del motivo de variación a
pequeña escala, la palabra efecto nugget podría ser reemplazada. Por ejemplo, en lugar de
efecto nugget podría llamársele efecto de fractura o efecto de estratigrafía. En general,
sería valioso reducir este tipo de efecto para mejorar la precisión de los estimativos.
Formas en las que uno puede reducir este efecto incluyen mejorar el control sobre los
medios que proporcionan los datos que serán analizados, reexaminar el grupo de datos para
hallar las posibles fuentes de error y aplicar factores de corrección.
2.2.1 Modelamiento del Semivariograma
El análisis de continuidad espacial se inicia con un variograma omnidireccional para el cual
la tolerancia direccional es lo suficientemente grande, tal que la dirección de un vector de
separación h cualquiera pierde importancia2. Así pues, con todas las posibles direcciones
ij
combinadas en un solo variograma, solo la magnitud del vector h es importante. Un
ij
variograma omnidireccional puede ser interpretado como un tipo de promedio de todos los
variogramas direccionales. No es estrictamente un promedio puesto que la localización de
los datos puede causar un opacamiento de aquellas direcciones en las que se tienen pocos
datos o, lo que es igual, una sobre representación de aquellas direcciones donde se presenta
el mayor número de datos.
12El cálculo de un variograma omnidireccional no implica la suposición de que la
continuidad espacial es la misma en todas las direcciones; este tipo de variograma
simplemente sirve como punto de partida para establecer algunos de los parámetros
requeridos en el modelamiento de la continuidad espacial. Por ejemplo, debido a que la
dirección no juega un papel importante en los cálculos de un variograma omnidireccional,
el análisis podría concentrarse en hallar los valores de g( h) y su respectiva tolerancia de tal
forma que se produzca la estructura más clara.
Otra razón por la cual es ventajoso iniciar con la construcción de un variograma
omnidireccional es que en este los datos erróneos son fácilmente detectables. El
variograma omnidireccional contiene más pares apareados que ningún otro variograma
direccional, por lo tanto, es de esperarse que este presente una estructura mas clara e
interpretable. Si el variograma omnidireccional no presenta una estructura interpretable, no
debería esperarse mucho éxito con un variograma para una dirección más particular. Un
análisis de los diagramas H de dispersión o el mapa de localización de los datos puede
revelar las causas de una estructura distorsionada. Una vez identificadas las fuentes de
error, se debe adaptar el modelo de variograma para que tenga en cuenta tales anomalías.
Un proceso de adaptación puede ir desde remover definitivamente algunos datos del grupo
o simplemente algunos pares en alguna dirección particular.
La estimación y posterior mapeo de una variable llevados a cabo en un método
geoestadístico debe ser precedido por el ajuste de un modelo al variograma observado. El
modelo escogido para ajustar el grupo de datos estudiado depende de consideraciones tanto
teóricas como prácticas. La mayoría de los semivariogramas experimentales pueden ser
descritos por unos cuantos modelos, y por esto la idea de un ajuste mediante modelos más
elaborados y complicados puede ser ignorada. Como se describirá en los siguientes
párrafos, un modelo para el semivariograma tendrá un uso crítico en el proceso de
estimación, el cual exige ciertas características e impone restricciones a las propiedades del
variograma. Una de las características fundamentales que debe satisfacer el semivariograma
es el criterio de admisibilidad. Este criterio se describirá después de presentar el grupo de
modelos más comunes.
Las dos características del semivariograma observado que servirán de guía para el ajuste de
un modelo son:
• Presencia o ausencia de sill
• Comportamiento cercano al origen
El sill ya ha sido previamente definido. En cuanto al comportamiento cercano al origen,
este se agrupa en dos tipos: comportamiento lineal y parabólico.
Otra característica importante que se debe observar en el semivariograma, es el valor que
este toma al aproximarse al origen. (Magnitud del vector h igual a 0). Esta característica ya
ha sido definida como el efecto nugget. Este mencionado efecto es importante y debe ser
13tenido en cuenta al momento de modelar el semivariograma.
Los tres modelos de ocurrencia mas frecuente en la literatura son:
• Esférico
• Exponencial
• Gaussiano
2.2.1.1 Modelo Esférico
Cerca al origen, el modelo esférico se comporta linealmente. Esta es una de las principales
características del modelo que puede ser explotada a la hora de ajustar el variograma3. En
este modelo, una tangente trazada al variograma desde el origen intercepta el sill a 2/3 del
rango.
La ecuación para el modelo esférico es la siguiente:
( )  3h  h3 
g h =C -    para h £ a

2a 2a3  
(7)
( )
g h =C para h >a
donde C y a corresponden al sill y rango respectivamente.
Sill = 480
Rango = 9
FIGURA 5. Modelo Esférico
Tomado de Hohn, “Geoestatistics and Petroleum Geology” 2da edición , 1999
14En el ejemplo mostrado en la Figura 5, el rango a es igual a 9 unidades y el sill es igual a
480 unidades.
2.2.1.2. Modelo Exponencial
Al igual que el modelo esférico, el modelo exponencial presenta un comportamiento lineal
en el origen. La diferencia es que en el modelo exponencial el variograma se aproxima de
forma asintótica al sill3. Este modelo esta descrito por la ecuación 8.
( )   3h
g h =C 1- exp-    (8)
a
La Figura 6 muestra un modelo exponencial con un sill de 480 unidades y un rango de 9
unidades. El rango de 9 es un rango práctico ya que el rango real es tres veces este.
Sill = 480
Rango = 9
FIGURA 6. Modelo Exponencial
Tomado de Hohn, “Geoestatistics and Petroleum Geology” 2da edición , 1999
La Figura 7 muestra ambos modelos con el fin de comparar el respectivo comportamiento.
Como se puede apreciar, el modelo exponencial se levanta más rápido en el origen pero se
estabiliza más lentamente, mientras que el modelo esférico muestra un cambio abrupto en
comportamiento cerca al sill. Este tipo de observaciones son las que se deben tener
presentes al momento de seleccionar un modelo.
2.2.1.3 Modelo Gaussiano
Este modelo, como los dos anteriores, posee un sill pero se comporta parabólicamente cerca
al origen3. La ecuación (9) sirve para representar un variograma de este tipo.
15FIGURA 7. Comparación del modelo Esférico y exponencial.
Tomado de Hohn, “Geoestatistics and Petroleum Geology” 2da edición , 1999
  h2 
( )
g h =C1- exp- 
(9)
 a2 
 
El modelo Gaussiano de la Figura 8 ajusta bien el comportamiento parabólico de un
semivariograma experimental que ha sido el resultado de datos tomados en la presencia de
una anomalía regional, la cual afecta en forma notable el comportamiento del
semivariograma cerca del origen. En ingeniería de petróleos, el modelo Gaussiano es
usado para analizar datos sísmicos en los cuales el espaciamiento de muestreo es muy
pequeño y el grado de continuidad es alto.
FIGURA 8. Modelo Gaussiano.
Tomado de Hohn, “Geoestatistics and Petroleum Geology” 2da edición , 1999
162.2.1.4 Modelos Combinados
Pocos semivariogramas experimentales tienen una apariencia similar a la de los modelos ya
descritos. Es común que como mínimo un efecto nugget mas un modelo esférico,
exponencial o gaussiano sea la forma más simple de ajustar el semivariograma
experimental. Un modelo escrito como la suma de dos o más modelos de menor
complejidad es lo que constituye un modelo combinado. Estrictamente hablando, la suma
del modelo más simple y un efecto nugget genera un modelo combinado aunque
usualmente no se le llame así. Por ejemplo:
( ) ( ) ( )
g h =g hg+ h
1 2 (10)
Donde
( h)=C  3h  h3  
g -   parah£ a
1 2a 2a3
(11)
g ( h)=C parah>a
1
y
( )
g h =C paratodoa
2 0 (12)
De hecho, en una forma más típica, un modelo de la forma siguiente es llamado un modelo
combinado:
g( h)=wg ( h)+wg ( h)+...+gw ( h)
1 1 2 2 k k (13)
Donde g (h) tiene en cuenta el efecto nugget y los demás modelos sumados describen las
1
demás características del semivariograma experimental. Los w son factores de peso cuya
k
suma debe ser igual al sill del variograma experimental.
Como un ejemplo de aplicación de un modelo combinado, considérese el caso de un
yacimiento naturalmente fracturado que produce de una formación que consta de varios
estratos geológicos intercalados. Es claro que en un caso como este son varios los factores
de variación de porosidad, permeabilidad o potencial inicial de gas o petróleo.
• A un nivel cercano al pozo, la variabilidad es atribuible al método de completamiento y
a los errores en la medición.
• A una escala de 0.5 Km., por ejemplo, la presencia de fracturas es un factor influyente.
17• A una escala de 10 Km., el espesor de las intrusiones impermeables juega un papel
importante.
Cada uno de estos factores de variación puede ser modelado dada una densidad de pozos
suficiente e incluso podemos incluir un efecto nugget para describir el tipo de variación
cercana al pozo. El modelo hipotético de la Figura 9 incluye tres términos:
• g (h) es el término de variación debida al error de medición.
1
• g (h) modela el comportamiento debido a la intersección de fracturas.
2
• g (h) tiene en cuenta la presencia de estratos impermeables.
3
FIGURA 9. Modelo Combinado
Tomado de Hohn, “Geoestatistics and Petroleum Geology” 2da edición , 1999
2.2.1.5 Modelamiento de Efectos de Oscilación
Algunos fenómenos a diario estudiados presentan una marcada periodicidad, por ejemplo el
tope de una formación en un estrato ondulado constituye un caso típico. Nótese que si el
estrato es ondulado, el valor del tope máximo se repetiría periódicamente a lo largo de la
estructura. Lo mismo pasaría para el valor mínimo. La consecuencia de esto, es que el
variograma, obtenido presentará una oscilación en los valores de varianza4 haciendo que el
gráfico del variograma luzca con un marcado efecto oscilatorio en donde no se alcanza un
grado de estabilización, es decir un valor de sill constante cuando el valor del vector h
aumenta.
Existen dos modelos comúnmente usados para describir el efecto de oscilación:
18( )  h
g h =C1- cos  (14)
a
 
( )
( )  sen ah  (15)
g h =C1- 
 ah 
donde h esta expresado en radianes.
La amplitud del efecto de oscilación es el valor mínimo de covarianza dividido por el valor
del sill:
( )
minC h
amplitud = ( ) (16)
C 0
El modelo de variograma de la ecuación (14), está representado con el esquema que se
muestra en la Figura 10. Nótese que la periodicidad de la variable que representa este
variograma, esta modelada con un efecto cosenoidal incluido en el lado derecho de la
ecuación (14).
FIGURA 10. Modelo cosenoidal de variograma para fenómenos de oscilación.
Tomado de Hohn, “Geoestatistics and Petroleum Geology” 2da edición , 1999
Un efecto de oscilación puede presentar una disminución gradual de la amplitud. El
modelo sinusoidal de la ecuación (15) presenta este tipo de atenuación ya que el valor del
seno (que varia entre -1 y 1) esta dividido por un valor que siempre esta aumentando. El
modelo cosenoidal obviamente no presenta esta característica pero puede agregársele
haciendo uso de un modelo combinado, como el siguiente de la ecuación 17.
19g( h) = exp- h  cos h  +C 3h h3
1-  -  (17)
  10 1  2 2 
En dicho modelo el término exponencial adiciona la característica en cuestión.
Las Figuras 11 y 12 muestran un esquema de las ecuaciones (15) y (17) respectivamente.
FIGURA 11. Modelo sinusoidal de variograma para fenómenos de oscilación
Tomado de Hohn, “Geoestatistics and Petroleum Geology” 2da edición , 1999
FIGURA 12. Modelo cosenoidal para modelos de oscilación gradualmente atenuados.
Tomado de Hohn, “Geoestatistics and Petroleum Geology” 2da edición , 1999
2.2.1.6 Otros Modelos de Variograma
Además de los modelos presentados anteriormente, existen otros que presentan
características menos comunes3.
El modelo más simple y que carece de sill es el modelo lineal:
20( )
g h =ph
(18)
donde p es la pendiente de una línea.
Un modelo más general, está representado por:
g( h)=phr
(19)
donde p es una constante y 0 < r < 2. Modelos con valores de r < 1 presentan
comportamiento semejante al modelo exponencial o esférico en las cercanías del origen.
Los modelos con r > 1 son parabólicos cerca al origen.
2.2.2 Admisibilidad Para Modelos de Variograma
Todos los modelos de semivariograma hasta aquí descritos exhiben una propiedad
importante en la posterior etapa de estimación; todos ellos son positivamente definidos5.
Dada una función aleatoria estacionaria Z(x) con covarianza C(h), sea Ψ una combinación
lineal con pesos l :
i
Y= ∑n l Z( x )
i i (20)
i=1
la cual es una variable aleatoria con varianza dada por:
n n ( )
VAR{Y =} ∑∑ll C x- x
i j i j (21)
i=1 j=1
y que debe ser mayor o igual a cero para que así C(h) sea positivamente definida. Así
pues, un modelo cualquiera de covarianza ajustado a un semivariograma experimental debe
cumplir que VAR{Ψ} > 0. Esta necesidad se debe al hecho de que las futuras estimaciones
serán combinaciones lineales de valores conocidos. Es decir que los sistemas de
ecuaciones que se tendrán que resolver para hallar los estimativos, deberán tener solución
real y además esta solución ha de ser única. El concepto de admisibilidad garantiza un
sistema de ecuaciones con estas características.
Por lo tanto, un modelo cualquiera para γ(h) debe probar admisibilidad. Afortunadamente,
modelos como los descritos cumplen admisibilidad. En cuanto a los modelos combinados
se tiene lo siguiente: “Un modelo resultado de una combinación lineal de modelos
21admisibles es también admisible”.
Este concepto de admisibilidad advierte que debe tenerse cuidado a la hora de ajustar
mediante polinomios arbitrarios un semivariograma experimental, ya primero se debería
probar la admisibilidad de dichos polinomios.
2.3 CONCEPTO DE ANISOTROPÍA GEOESTADÍSTICA.
Muchos fenómenos geológicos presentan una clara anisotropía espacial en cuanto a la
continuidad. Es decir, el semivariograma exhibe mayores cambios en sus parámetros como
lo son el sill y el rango a medida que la dirección cambia. El sill y el rango son ahora una
función de la dirección. El ejemplo en la Figura 13 parte (a) muestra una vista isométrica
de la superficie generada por un variograma en la cual el rango cambia con la dirección
pero el valor del sill permanece constante. Este tipo de anisotropía es conocida como
anisotropía geométrica. En el caso de anisotropía zonal el valor del sill cambia en función
de la distancia mientras el rango del semivariograma permanece constante. El ejemplo de la
Figura 13 parte (b) es la muestra una mezcla de estas dos anisotropías donde tanto el rango
como el valor del sill cambian con la dirección.
FIGURA 13. Vista isométrica de una anisotropía (a)Geométrica, (b) Zonal.
Tomado de Srivastava e Isaaks, “Applied Geostatistics”, 1989
Dado un grupo de variogramas que muestren el rango y/o el sill cambiando con la
dirección, se debe empezar por identificar los ejes de anisotropía. Dichos ejes
corresponden a las direcciones en las cuales se presentan el mayor y menor rango o el
máximo y mínimo sill en caso de anisotropía zonal. Mapas de contorno son
extremadamente útiles para determinar estas direcciones. Otra alternativa, además de los
mapas de contorno, es una gráfica en coordenadas polares de todos los variogramas
experimentales disponibles.
22Por otro lado, información cualitativa como la orientación litológica o planos de
buzamiento es también útil al momento de identificación de los ejes de anisotropía.
Igualmente, un conocimiento de la génesis del fenómeno bajo estudio también puede ser de
utilidad.
Una vez identificados los ejes de anisotropía, el próximo paso es la creación de un modelo
que describa el comportamiento del variograma con los cambios en distancia y dirección.
Por ahora, se trabajará en el sistema coordenado definido por los ejes de anisotropía. Una
vez deducido como construir un modelo en este sistema, se pasará a buscar el método que
habilitará al modelo creado para ser usado en el sistema coordenado original de los datos.
Un método para combinar varios semivariogramas direccionales en uno solo que describa
todas las direcciones, es definiendo una transformación que los reduce a un modelo común
con un rango estandarizado de 1. El procedimiento consiste en transformar la distancia h
de separación de tal forma que el modelo estandarizado arrojará un valor de variograma
que es idéntico a cualquiera de los modelos direccionales para esa distancia h de separación
en particular.
Por ejemplo, dos modelos de semivariograma son mostrados con idéntico valor de sill en la
Figura 14. Uno de ellos tiene un rango de 1 mientras el otro tiene un rango de a. Nótese
que si se evalúa el modelo cuyo rango es 1a una distancia de h/a se obtiene el mismo valor
que si se evaluara el modelo de rango a a una distancia h. De esta forma, se ha reducido el
modelo con un rango de a a un modelo equivalente de rango 1 simplemente mediante el
cambio de la distancia h por una distancia reducida h/a. Esta equivalencia puede ser
expresada como:
g  h =g ( h)
1a a (22)
o lo que es equivalente,
( ) ( )
g h =g ah
1 a (23)
FIGURA 14. Reducción de una anisotropía geométrica
Tomado de Srivastava e Isaaks, “Applied Geostatistics”, 1989.
23Este concepto de un modelo equivalente y distancia reducida puede ser extendido a dos
dimensiones. Si a es el rango en la dirección x y a es el rango en la dirección y, entonces
x y
el modelo para este variograma anisotrópico puede ser expresado como:
g( h)=g ( h ,h ) g= ( h )
x y 1 1 (24)
donde la distancia reducida h esta dada por:
1
h 2 h 2
h 1 =  ax  +   ay   (25)
x y
siendo h y h las componentes del vector h en el eje x y y respectivamente.
x y
Similarmente, un modelo equivalente para un variograma isotrópico en tres dimensiones,
con rangos a , a y a , puede ser expresado como:
x y z
g( h)=g ( h ,h ,h ) g= ( h )
x y z 1 1 (26)
donde la distancia reducida h esta dada por:
1
h 2 h 2 h 2
h 1 =  ax  +   ay   +  az  (27)
x y z
Este método de usar modelos equivalentes y distancias reducidas también funciona con
modelos que no alcanzan un sill, es decir, con modelos lineales. La aplicación que esto
tendría en el modelamiento de una anisotropía geométrica se presenta en el siguiente
ejemplo:
Recuérdese que la anisotropía geométrica de un grupo de datos esta caracterizada por
semivariogramas direccionales que presentan aproximadamente el mismo sill pero difieren
en rango (en el caso de variogramas lineales, estos presentarán un cambio de pendiente con
la dirección). La Figura 15 muestra tres variogramas direccionales a lo largo de los tres
ejes perpendiculares de anisotropía. Cada modelo direccional consta de solo una estructura
(no combinado) y los tres alcanzan idéntico valor de sill. Sin embargo, sus rangos son
diferentes. Un modelo tridimensional equivalente para el caso presentado en la Figura 12
esta dado por:
24( ) ( )
g h =wg h (28)
1 1 1
y la distancia reducida h esta dada por la ecuación 27, donde a , a y a son los rangos de
1 x y z
los modelos direccionales a lo largo de los ejes de anisotropía; h , h y h son las
x y z
componentes de h en las direcciones x, y y z de anisotropía y γ (h ) es el modelo
1 1
equivalente con rango estandarizado de 1.
FIGURA 15. Anisotropía Geométrica
Tomado de Srivastava e Isaaks, “Applied Geostatistics”, 1989
2.4 UN EJEMPLO SIMPLE
Si bien la teoría hasta aquí desarrollada es sencilla, no hay mejor forma de asimilar el
concepto de variografía que a través de un ejemplo práctico. El siguiente análisis
corresponde al potencial inicial de un campo de gas. De este grupo de datos se sabe que la
producción puede provenir de una o más formaciones productoras6. El tamaño y la
estratigrafía de cada zona completada obviamente varía de pozo a pozo. Idealmente, uno
estaría interesado en estudiar el potencial de cada zona productora por separado, pero en la
realidad la producción total de un pozo es la contribución de varias zonas productoras. Por
lo tanto, los volúmenes están mezclados en los datos reportados. Así pues, en un caso
como este, un análisis geoestadístico modelará el comportamiento típico de un pozo y no el
de una formación productora en particular.
La gráfica del semivariograma omnidireccional de la Figura 16 muestra la forma general
de un modelo esférico o tal vez exponencial con un posible efecto nugget. El primer intento
de ajustar un modelo usa un modelo esférico como el presentado en la Figura 17; la
ecuación de este modelo es:
( )  3h h3 
g h =1.8· 106  2(cid:215) 0.9 - 2(cid:215) 0.93 (29)
25FIGURA 16. Semivariograma Omnidireccional
Tomado de Hohn, “Geoestatistics and Petroleum Geology” 2da edición , 1999
En general, los parámetros usados en este primer intento parecen ajustar bien, pero el
modelo crece demasiado rápido para valores pequeños de h. Un rango más largo provee un
mejor ajuste en esta región pero alcanza el sill a una distancia muy corta. Por dichas
características el modelo esférico parece estar equivocado sugiriendo que el uso de un
modelo exponencial podría dar mejores resultados.
FIGURA 17. Ajuste del semivariograma omnidireccional con un modelo esférico.
Tomado de Hohn, “Geoestatistics and Petroleum Geology” 2da edición , 1999
La Figura 18 muestra el ajuste dado mediante un modelo exponencial. Los resultados en
este caso son mucho mejores. El modelo exponencial tiene un sill de 1'800.000 Mpcpd2
(miles de pies cúbicos por día) y un rango de 0.4 Km. Por lo tanto, la ecuación para el
modelo es:
g( h) =1.8· 106 -1 ex-p  h    (30)
 0.4
 
El ajuste es tan bueno que se podría pensar que los datos fueron intencionalmente
calculados para el ejemplo. Sin embargo, datos reales fueron usados y la buena apariencia
26del semivariograma probablemente se debe al gran número de datos disponibles, a la
homogeneidad del potencial en el área de estudio y al buen número de pozos cercanos que
permiten una excelente definición del variograma cerca del origen. El número de pares por
punto calculado del semivariograma esta reportado en la tabla del Anexo A. Es notable que
incluso para pequeñas distancias de separación, el número de pares sea bastante alto.
FIGURA 18. Ajuste final del semivariograma mediante modelo exponencial.
Tomado de Hohn, “Geoestatistics and Petroleum Geology” 2da edición , 1999
Los resultados indican que el modelo de la ecuación 30 representado en la Figura 18, es el
mejor ajuste para los datos del semivariograma omnidireccional de la Figura 16. Es de
destacar que el semivariograma crece rápidamente y alcanza el sill a una distancia de
separación de 2 Km. El significado práctico de esta observación es que se podría descartar
cualquier estimación utilizando datos a más de 2 Km. de distancia del punto de control.
2.5 MODELOS DE ESTIMACIÓN
Cuando en ingeniería de yacimientos se realiza el muestreo de una variable regionalizada,
es evidente que por cuestiones económicas no es posible tener un dato para cada punto
dentro del área que en la que se distribuye la variable. La estimación permite predecir los
valores de este tipo de variables en sitios donde no se tiene una muestra mediante el uso de
los valores de la variable en donde se ha realizado un muestreo.
La mayoría de los métodos de estimación tienen en cuenta una combinación lineal de los
datos disponibles por medio de la asignación de un peso de combinación w. Estos pesos
i
son generalmente estandarizados para que su suma sea igual a 1, sin embargo esto no es
siempre necesario. La forma en que se asignen los pesos a los valores de los datos da origen
a diversas metodologías de estimación.
27n
estimado=v~ =∑ w (cid:215) v (31)
i i
i=1
La estimación requiere un modelo de como se comporta el fenómeno que se va a estudiar
en los lugares donde este no se ha muestreado: sin un modelo no se podrían realizar
inferencias acerca de los valores desconocidos de la variable regionalizada.
La estimación de los valores desconocidos requiere una información adicional o la
elaboración de algunos supuestos. Lo ideal sería conocer los procesos físicos y/o químicos
involucrados en la generación del fenómeno estudiado. En las situaciones donde se conoce
este tipo de información lo apropiado seria utilizar un modelo determinístico de estimación.
Desafortunadamente, en los procesos involucrados en las ciencias de la tierra, no se tiene
un conocimiento tal que se pueda utilizar este tipo de modelos. Para la gran mayoría de los
datos asociados con fenómenos de geociencias, es necesario admitir que existe una
incertidumbre de como se comporta el fenómeno entre las ubicaciones de las muestras.
Debido a que los modelos determinísticos parecen ser inútiles en este tipo de estudios, es
necesario traer a este trabajo una breve descripción de los modelos probabilísticos de
estimación.
La aproximación geoestadística de estimación utiliza los modelos probabilísticos ya que
estos tienen en cuenta la incertidumbre de lo sucedido en la generación del fenómeno que
se quiere estudiar. En un modelo probabilístico, la información disponible es vista como el
resultado de un proceso aleatorio. Por ejemplo, si tratáramos de describir todos los procesos
involucrados en la formación de un yacimiento de petróleos, tendríamos un
desconocimiento tan grande de todo la situación, que su complejidad aparecería como la de
un comportamiento aleatorio; pero esto no significa que los procesos son en verdad
aleatorios, simplemente lo que sucede es que se ignora lo sucedido.
Es importante notar el hecho de que la palabra aleatorio no significa impredecible, y
tampoco significa que considerar un problema de geociencias como un fenómeno aleatorio
quiera decir que el resultado de la estimación pueda generarnos cualquier valor, en
particular uno que carezca de sentido físico.
La mayoría de la teoría de geoestadística es referida en términos de variables aleatorias y
las herramientas más utilizadas son la media y la varianza de una combinación lineal de
variables aleatorias.
El conjunto de valores de salida de las variables aleatorias y sus correspondientes
probabilidades son conocidos como distribución de probabilidad de una variable aleatoria.
Si la distribución de probabilidad es conocida, se podrían calcular muchos parámetros de
interés de la variable aleatoria, como por ejemplo, el valor máximo, el mínimo, la media y
la desviación estándar.
Una distribución de probabilidad no puede ser determinada con solo unos parámetros; la
media y la desviación estándar, aunque de gran ayuda, no son suficientes, ya que dos
28variables aleatorias podrían tener la misma media y la misma desviación, pero aun así sus
distribuciones pueden ser muy diferentes. Sin embargo, si conocemos que una variable
aleatoria se distribuye de una forma Gaussiana la media y la varianza serán suficientes para
definirla completamente.
2.5.1 Estimación Lineal - Kriging
La estimación lineal es simplemente un proceso de interpolación que busca determinar una
aproximación del valor de una variable en un punto donde esta no es conocida, a partir de
los valores de la misma variable en lugares que si poseen un valor medido. El kriging es un
método de estimación lineal que calcula los estimados de una variable regionalizada en un
punto, sobre un área o dentro de un volumen, utilizando un criterio de minimización de una
varianza de estimación. Los valores estimados mediante kriging pueden ser utilizados para
dibujar mapas de contorno cuando estos son calculados en los nodos de una malla regular.
Existen varios algoritmos de kriging que pueden ser usados dependiendo de la información
que se tiene en el momento de comenzar con el análisis geoestadístico. En esta sección se
mostrará el desarrollo de las ecuaciones fundamentales del kriging ordinario.
El método de kriging ordinario esta asociado con la sigla BLUE (best linear unbiased
estimator) que significa mejor estimador lineal imparcial. Es lineal porque los estimados
que resultan son combinaciones lineales de los datos disponibles, es imparcial porque trata
de tener el error igual a cero y es el mejor porque trata de minimizar la varianza de los
errores. Estos dos últimos objetivos resultan algo ambiciosos ya que en un sentido practico
no se conocen los valores reales de los datos estimados, por lo tanto el valor del error como
de la varianza de los errores es desconocido. Lo que hace el kriging ordinario es usar un
modelo de probabilidad en donde se escogen pesos para las muestras cercanas que aseguren
que el error sea igual a cero y que la varianza sea minimizada. El desarrollo del modelo
tendrá en cuenta el uso de una función aleatoria que permitirá expresar el error, su valor
medio y su varianza
2.5.1.1 El Modelo de Función Aleatoria y el Concepto de Estimador Insesgado:
En todo punto donde no se tenga una muestra, se estimará el valor desconocido utilizando
una combinación lineal de los valores de las muestras disponibles,
n
u$=∑ w u(cid:215)
j (32)
j=1
El símbolo ^ es usado para denotar un estimado de un modelo y distinguirlo de los datos
medidos. Los w representan los pesos para el modelo.
j
29Si definimos el error, r, de un valor estimado como la diferencia entre este y el valor real en
la misma posición,
r =u$ - u
i i i (33)
Entonces el error promedio de un número k de estimados será:
m =1 ∑n u$ -u
r k i i (34)
i=1
El problema de la ecuación anterior es que no se conocen los valores reales u ,.......,u . La
1 k
solución probabilística a este problema consiste en considerar a los valores desconocidos
como datos de salida de un proceso aleatorio. Asumiendo que la función aleatoria es
estacionaria, entonces cada una de las variables aleatorias tendrá un valor esperado E{V}, y
la covarianza entre parejas de variables aleatorias separadas por una distancia particular, h,
será, Ĉ (h). La clave de este modelo es tener en cuenta que las muestras son valores de una
V
variable aleatoria, y como los estimados son combinaciones lineales de las muestras,
entonces por lo tanto estas también serán variables aleatorias:
n
V$ =∑
(x ) wV(x )
0 i i (35)
i=1
De forma similar, el error de estimación, definido como la diferencia entre el estimado y el
valor real de la variable aleatoria modelada, será también una variable aleatoria:
R(x ) =V$ (x )- V(x ) (36)
0 0 0
Sustituyendo la ecuación (35) en la (36) se obtiene una expresión para el error que se
comete cuando se estima el valor desconocido en el punto x
0.
n
R(x )=∑ w (cid:215) V(x )- V(x )
0 i i 0 (37)
i=1
La forma de asegurar que el error de estimación en una ubicación particular tiende a cero,
es determinando la forma que tomaría la ecuación (37) cuando el valor esperado del error
es cero.
30n
E{R(x )} =E{∑ wV(x ) - V(x )}
0 i i 0
i=1
n
E{R(x )} =∑ w E{V(x )} - E{V(x )}
0 i i 0
i=1
n
E{R(x )} =E{V}∑ w - E{V}
0 i
i=1 (38)
La condición de imparcialidad se obtendrá cuando el valor esperado del error sea igual a
cero. Igualando la ecuación (38) a cero se concluye que la condición de imparcialidad del
modelo se cumplirá cuando la suma de los pesos de la combinación lineal sea igual a 1.
n
0=E{V}∑ w - E{V}
i
i=1
n
∑ w =1
i
i=1
(39)
2.5.1.2 El Modelo de Función Aleatoria y la Varianza del Error.
El kriging ordinario tiene como finalidad generar estimados de tal forma que se minimice la
varianza del error. Como en la sección anterior, se utilizaran las variables aleatorias para
encontrar un conjunto de ecuaciones que represente la minimización de la varianza del
error.
Ya hemos definido el error del estimado por medio de la ecuación (36). Además es
conocido por medio de la ecuación (38), la forma en que se puede expresar la varianza de la
combinación lineal. Por lo tanto la varianza del error esta dada por:
Var{R(x )} =Cov{V$ (x )V$ (x0)} - 2Cov{V$ (x )V(x )}+ Cov{V(x )V(x )} (40)
0 0 0 0 0 0
La covarianza de un parámetro consigo mismo es igual a la varianza del parámetro, por lo
tanto el primer término de la parte derecha de la ecuación (40) pueden ser expresado de la
siguiente manera:
Cov{V$ (x )V$ (x )} =Var{V$ (x )} =Var{∑n w (cid:215) V }= ∑n ∑n w w C~
0 0 0 i i i j ij (41)
i=1 i=1 j=1
31Asumiendo que todas las variables aleatorias tienen la misma varianza entonces el tercer
término del lado derecho de la ecuación (41) puede ser reemplazado por:
Cov{ V( 0) V( 0)} =Var{ 0)} =s~2
x x V(x (42)
Utilizando la definición de la ecuación (35) se puede escribir el segundo término de la
derecha de la ecuación (41) como:
n
2Cov{V$ =2Cov{(∑
(x )V(x )} wV )V }
0 0 i i o
i=1
n
=2∑ w (cid:215) Cov{VV }
i i 0
i=1
=2∑n C~
w
i i0
i=1
(43)
Combinando los tres términos obtenidos en las ecuaciones (41), (42) y (43) y
reemplazándolos en la ecuación (40) se obtiene la siguiente expresión para la varianza del
error:
s~2 =s ~2 +∑n ∑n w w (cid:215) C~ - 2∑n w C~
R i j ij i i0 (44)
i=1 j=1 i=1
Una vez escogidos los parámetros de la función aleatoria, esto es, la varianza σ2 y las
covarianzas Ĉ , la ecuación (44), se tendrá una expresión para calcular la varianza del error
ij
como una función de los pesos w ,.......,w . La minimización de la varianza del error se
1 n
puede realizar calculando las primeras derivadas parciales e igualándolas a cero. Este
procedimiento generaría un sistema de n ecuaciones con n incógnitas que puede ser resuelto
con un algoritmo para un sistema de ecuaciones lineales. Desafortunadamente este
procedimiento resulta incorrecto ya que la solución obtenida esta limitada por el concepto
de imparcialidad (la suma de los pesos debe ser igual a 1). Para evitar este problema se
tendrá en cuenta la técnica de los multiplicadores de Lagrange.
2.5.1.3 El parámetro de Lagrange.
Si al sistema generado por las n primeras derivadas parciales de la ecuación (44) se le
agrega la ecuación del concepto de imparcialidad, se estaría creando un sistema de n+1
ecuaciones con solo n incógnitas (sistema mal definido y que no tiene una solución única).
Para evitar este problema, se puede introducir otra incógnita dentro de la ecuación (44) por
32medio de una nueva variable denominada como el parámetro de Lagrange (µ). Escribiendo
la ecuación (44) con el nuevo término que incluye el parámetro de Lagrange se tiene:
s~2 =s ~2 +∑n ∑n w w C~ - 2∑n w C~ +m 2 (∑n w- 1)
R i j ij i i0 i (45)
i=1 j=1 i=1 i=1
Hay dos cosas importantes que resaltar en la ecuación (45): la primera es que la
introducción del nuevo termino no afecta la ecuación ya que este es igual a cero, y la
segunda es que se esta teniendo en cuenta el concepto de imparcialidad al introducir la
sumatoria de los pesos y dar por entendido que esta sumatoria es igual a 1.
Ahora bien, si se expresan las n+1 derivadas parciales de la ecuación (45) se obtendrá un
sistema de n+1 ecuaciones con n+1 incógnitas, en donde la nueva incógnita seria el
multiplicador de Lagrange. En particular si se calcula la derivada parcial de σ 2 con
R
respecto a µ y se minimiza (igualar la derivada a cero) se obtendrá la condición de
imparcialidad de la ecuación (39).
2.5.1.4 Minimización de la Varianza del Error.
Para minimizar la varianza del error se calculan las n+1 primeras derivadas de la ecuación
(45). La derivación con respecto a w se muestra en detalle y las demás derivadas con
1
respecto a los demás pesos se calculan de una forma análoga.
Expandiendo la doble sumatoria de la ecuación de la ecuación (45) y eliminando todos los
términos que no incluyen w se obtiene:
1
¶∑n ∑n w w C~  ¶w2C~ +2w ∑n w C~ 
 i=1 j=1 i j ij  1 11 1 j=2 j 1j
=
¶w ¶w
1 1
C~ +2∑n C~
=2w w
1 11 j 1j
j=2
=2∑n C~
w
j 1j
j=1
(46)
El tercer término de la derecha de la ecuación (45) solo contiene un término que involucra a
w . La derivada quedará de la siguiente manera:
1
¶ ∑ in =1w iC~ i0  =¶( w 1C~ 10)
=C~ (47)
¶w ¶w 10
1 1
33De igual manera el último término del lado derecho de la ecuación (45) solo contiene un
termino que involucra w . La derivada estará representada por
1
¶m   ∑ in =1w - 1   =¶m( 1)
i w (48)
=m
¶w ¶w
1 1
Haciendo uso de las ecuaciones (46) a (48) se puede escribir la primera derivada de σ 2 con
R
respecto a w de la siguiente forma:
1
( )
¶s ~2
R =2∑n w C~ - 2C~ + 2m (49)
¶w j 1j 10
1 j=1
Minimizar la función de la varianza del error es igualar las derivadas parciales a cero.
Igualando la ecuación (49) a cero, se obtiene la siguiente expresión:
∑n C~ =C~
w +m
j 1j 10 (50)
j=1
La derivación con respecto a los demás pesos producirá ecuaciones similares a la ecuación
(50). También se vió en la sección anterior que la derivación con respecto a m genera la
condición de imparcialidad. Estas dos situaciones llevan a concluir que el conjunto de pesos
que minimizan la varianza del error son aquellos que cumplen con las siguientes dos
ecuaciones:
∑n C~ =C~
w +m " =i 1,.......,n
j ij i0
j=1
(51)
n
∑ w =1
i
i=1
El sistema de ecuaciones referido comúnmente como “kriging ordinario” puede ser escrito
en notación matricial como:
34C~ C~ C~
L 1 w  
 11 1n   1   10
M O M M M M
 (cid:215)  = 
 C~ L C~ 1 w   C~ 
 n1 nn   n  n0
 1 L 1 0 m   1 
C(cid:215) w =D (52)
Todas las covarianzas en el sistema de kriging de la ecuación (52) se calculan a partir del
semivariograma y de las distancias entre las muestras y los puntos a ser estimados. Cada
entrada en la matriz C es una covarianza de muestra a muestra, es decir de solo valores en
donde la variable corregionalizada ha sido medida. Dicha matriz depende solo de las
ubicaciones de las muestras y no contiene ninguna información acerca de la ubicación
geográfica en la que se encuentra el estimado deseado. Si dos muestras están muy cerca la
una de la otra, esto será registrado como un valor alto de covarianza en la matriz C; lo
contrario sucederá cuando las muestras están muy lejanas. Esta matriz es generalmente
constante, ocasionando un ahorro en el tiempo computacional.
El vector D es un vector de covarianzas entre los puntos de las muestras y el punto a ser
estimado. Cabe recordar que para calcular las covarianzas se requiere tener el dato de la
distancia entre el punto de la muestra y el punto a ser estimado. Es claro que el sistema de
ecuaciones del kriging ordinario dependen de las covarianzas que sean escogidas para
calcular las matrices C y D. Escoger un modelo de covarianza es un requisito necesario del
kriging, lo que hace que los cálculos se tornen un poco tediosos. Sin embargo es
importantísimo resaltar que la inclusión del modelo de covarianzas genera un mecanismo
que permite incluir el patrón de anisotropías de la variable regionalizada que se este
estudiando.
En la práctica, el patrón de continuidad espacial que se escoge para la función aleatoria es
el mismo que se determina del grupo de datos de la muestra. Una vez que se calcule el
variograma, se le ajusta un modelo mediante una función, por medio de la cual se podrán
determinar las covarianzas para las matrices del sistema de kriging. Sin embargo esta
práctica no siempre puede resultar exitosa, debido a que la elaboración del variograma, y
por consiguiente del modelo de continuidad espacial, se puede dificultar por falta de datos o
errores en el muestreo, o por posibles valores erráticos.
Otra forma de encontrar un modelo de continuidad espacial es mediante una apreciación
cualitativa de los datos que se tienen. La experiencia con otros campos similares al que se
este estudiando (datos similares de las propiedades), puede ser una buena guía para
seleccionar el patrón de continuidad espacial, en particular en campos nuevos donde por
falta de información no se puede determinar un buen variograma.
35La selección del modelo del variograma es quizás el paso más importante en todo el
proceso geoestadístico ya que este afectará de una forma directa los cálculos de estimación
que se realicen con el kriging. Algunos parámetros importantes del variograma, que afectan
los resultados del kriging son:
Efecto de forma: Un comportamiento parabólico al origen del variograma es un indicativo
de un fenómeno muy continuo, por lo tanto el mecanismo de estimación en un punto utiliza
mucho más los valores de muestras cercanos. Cuando esto ocurre, suelen calcularse pesos
cercanos a cero o en ocasiones negativos. Este fenómeno se conoce como efecto pantalla.
Se dice que una muestra fue afectada por la pantalla cuando otra muestra se interpone entre
ella y el punto a ser estimado. El grado en que una muestra pierde su influencia sobre un
estimado a causa de la pantalla depende del patrón de continuidad espacial. La ventaja de
obtener pesos negativos o pesos mayores que 1, radica en el hecho de que se pueden
generar estimados que sean menores que el mínimo valor muestreado, o mayores que el
máximo valor muestreado; esto es importante porque es muy poco probable que los valores
extremos de la variable estudiada hayan sido muestreados.
2.5.2 Corregionalización.
Existen varias situaciones en las que seria conveniente sacar provecho de la covarianza
entre dos o más variables regionalizadas:
• La variable de interés es una combinación lineal de otras variables regionalizadas.
• La cantidad de datos de la variable de interés es poca, pero sin embargo esta altamente
correlacionada con otra variable de la cual el numero de datos es alto. Se puede sacar
ventaja de la correlación para mejorar la estimación de la variable que presenta
deficiencia en el número de datos.
• Una variable exhibe baja autocorrelación espacial pero muestra alta correlación con otra
variable de buena autocorrelación espacial. Nuevamente, esta situación es
aprovechable para mejorar la estimación de aquella variable que presenta baja
continuidad espacial.
Al análisis del comportamiento espacial apareado de variables regionalizadas es conocido
como la técnica de cokriging. Dentro de esta técnica, que permite la estimación de una
variable a partir de otra o más, existen varios métodos como cokriging, kriging con
parámetro externo y kriging colocado.
Considérese un caso de estimación del valor z* a partir de dos grupos de datos z , z , ..., z y
1 2 n
y , y , ..., y . Llamaremos a este valor estimado z* para indicar que ha sido calculado a
1 2 m k
través de alguna de las técnicas de corregionalizacion. La tarea que se debe hacer en este
momento es calcular dos grupos de pesos l , l , ..., l y k , k , ..., k que serán utilizados
1 2 n 1 2 m
para el cálculo del nuevo valor estimado:
36n
z* =∑l z +∑k y
ck i i i i (53)
i=1
Cokriging hace uso de la covarianza entre la variable principal o primaria z y la variable
secundaria y. Si el muestreo de la variable primaria es seriamente pobre comparado con el
caso de la variable y, la intención es que la variable secundaria sea lo suficientemente bien
muestreada y que la covarianza entre las dos variables, z y y, lo suficientemente alta para
mejorar la estimación de la variable primaria. Dicho de otra forma, la idea es minimizar la
incertidumbre al momento de la estimación.
Por ultimo, es muy importante señalar que gran parte de la importancia de la estimación a
partir del modelamiento de una corregionalizacion radica en que este es menos susceptible
a efectos de homogenización. Consideremos un caso simple en el que dos pozos están
perforados en el mismo estrato productor. El primer pozo registra saturación de agua de
0.25 y el segundo registra saturación de agua igual a 0.4. Una técnica de estimación
cualquiera, incluso una como kriging, predecirá un valor entre este rango para un nuevo
pozo perforado en medio de los dos anteriores. Ahora bien, supóngase que a causa de
fenómenos de capilaridad la zona en medio presenta alta presión capilar, es decir, alta
saturación de agua. Entonces un método simple de estimación arrojará valores erróneos a
menos que de alguna forma involucremos datos de presión capilar en la técnica de
estimación. La técnica que nos permite hacer esto es la corregionalización.
2.5.2.1 Modelamiento de una Corregionalización7,8
Continuando con el caso más común de una variable primaria y una variable secundaria, el
primer paso en cokriging es el cálculo y modelamiento del semivariograma para cada
variable. Además, el semivariograma cruzado debe ser calculado y modelado. Si más de
una variable secundaria ha de ser incluida, un semivariograma cruzado debe ser calculado
para cada par de variables incluida la variable primaria. Siendo n el número total de
variables primarias y secundarias, el número de semivariogramas cruzados ha examinar
será n(n + 1)/2. Esto puede constituir una de las principales desventajas de la técnica pues
el modelamiento de los semivariogramas es la etapa que demanda más tiempo y cuidado.
2.5.3 Kriging con Parámetro Externo8
Considérese un caso en el cual la variable secundaria es conocida en todo punto del área
estudiada. Por lo tanto, se conoce el valor de tal variable en cada nodo de una malla
arbitraria como también en todo punto donde la variable primaria es conocida. Por
37ejemplo, sea la profundidad de una formación productora la variable primaria; esto debido a
que los datos de profundidad son solo observables directamente en los pozos. La variable
secundaria, debido a que un barrido sísmico puede proveer una cobertura total del área,
podría ser el tiempo de viaje de una onda sísmica hasta el horizonte productor.
Kriging con parámetro externo es un estimador lineal como todos los anteriores que se han
visto hasta aquí:
n
z* =∑l z
kpe i i (54)
i=1
La diferencia respecto al kriging ordinario esta en el sistema de ecuaciones:
∑n C( ) y( ) =C( ,V)
l v ,v +m m+ v v paratodoi =1,n
j i j 1 2 i i (55)
j=1
n
∑l =1
j (56)
j=1
∑n y( ) =y( V)
l v
j j (57)
j=1
Nótese que solo hay un grupo de pesos; y que la variable secundaria no es usada en la
estimación directamente pero si se muestra en una de las restricciones.
El método de kriging con parámetro externo requiere dos grupos de datos. El primero es el
valor de ambas variables, tanto primaria como secundaria, en cada pozo. El segundo grupo
lo constituyen las observaciones de la variable secundaria en cada punto de la malla donde
la variable primaria ha de ser estimada. Debe proveerse también un semivariograma para la
variable primaria. Finalmente, el método no requiere un variograma cruzado o un
variograma de la variable secundaria.
2.5.4 Kriging Colocado9
El cokriging ordinario puede ser utilizado para mejorar las estimaciones de, por ejemplo,
porosidad a partir de propiedades sísmicas. El hecho de que los datos sísmicos poseen un
muestreo más denso comparado con los datos obtenidos de pozo, combinado con el efecto
de pantalla del kriging, conduce a que unos pocos valores de la variable secundaria son
necesarios en las ecuaciones de cokriging. Solo el valor de la variable secundaria más
38cercano al punto de estimación recibe un alto peso. Si la malla usada para la estimación es
igual a aquella donde la variable secundaria ha sido medida, entonces los nodos son
llamados colocados.
Esta observación significa que las ecuaciones de cokriging pueden ser significativamente
simplificadas si se incluyen valores adyacentes de la variable primaria y solo un valor
colocado de la variable secundaria.
2.6 SIMULACIÓN GEOESTADISTICA
En la caracterización de yacimientos la estimación propuesta por la geoestadística ha tenido
un amplio uso a través del procedimiento del kriging. Este método, de acuerdo a las
condiciones con que es obtenido, así como la amplia gama de interpoladores conocidos,
muestran una imagen suave de la realidad que se estudia, sin tener en cuenta las
fluctuaciones que se presentan en la práctica. Una alternativa para esta situación es simular
el comportamiento de las variables analizadas10. Considerando las variables como
funciones aleatorias, se puede obtener una de las posibles realizaciones del fenómeno con
las mismas características de variabilidad y correlación espacial de los datos originales y
por tanto de la realidad, características que son reveladas de la información disponible a
través del cálculo de los semivariogramas experimentales y el ajuste de modelos teóricos.
La simulación no debe ser considerada superior a la estimación, sino como una alternativa
en el grado de descripción que se desee realizar. Los objetivos de estos procedimientos son
diferentes, pero en muchos casos el análisis de sus resultados paralelamente puede ayudar a
enriquecer la caracterización que se desarrolla en el estudio de variables regionalizadas.
2.6.1 Métodos de Simulación Condicional
La siguiente discusión de métodos para simulación considera solo el caso en el cual se
pretende crear una realización del fenómeno estudiado que además de proveer nuevos
valores estimados, reproduzca exactamente los datos proporcionados. A este tipo de
simulación se le conoce como simulación condicional. Usualmente, en la práctica, se está
interesado en este tipo de simulación.
Los métodos de simulación geoestadística se agrupan en aquellos cuyo resultado en el
primer paso es un resultado condicionado y aquellos que primero producen un resultado no
condicionado y en un paso posterior condicionan la simulación a los datos proporcionados.
Al primer grupo pertenecen los llamados métodos secuenciales.
2.6.1.1 Métodos secuenciales
En la simulación secuencial, un valor simulado en cada ubicación u es tomado de una
distribución de probabilidad la cual es calculada a partir de los datos suministrados y de
39valores previamente simulados que se hallan en la vecindad de u. El algoritmo se inicia
con la selección aleatoria de un punto u y se avanza secuencialmente a través de la malla en
que se ha divido el área bajo estudio. El orden de la secuencia en la cual se recorre la malla
es aleatorio.
Para cada punto u de la malla se buscan puntos en una región de vecindad definida por el
modelo de continuidad espacial. Estos puntos pueden incluir tanto datos previamente
simulados como los suministrados inicialmente. Hecho esto, una distribución de
probabilidad es calculada a partir de tales puntos. La forma de cálculo de esta distribución
de probabilidad es la que da origen a diversos tipos de simulación secuencial. Dos de ellos
son simulación gaussiana secuencial y simulación de indicador secuencial.
2.6.1.2 Simulación gaussiana secuencial
Este método utiliza la técnica kriging ordinario aplicado a los datos iniciales después de
trasformarlos a un espacio normal. Este tipo de simulación requiere un solo modelo de
semivariograma pero basado en los datos transformados. Cuando la simulación se ha
realizado, los resultados son transformados nuevamente al espacio original.
Ambos tipos de simulación requieren algunas decisiones respecto a las colas de las
distribuciones. Debido a que el número de nodos es mucho mayor que los datos
originalmente suministrados, estos no son suficientes para deducir un comportamiento
acertado de esta parte de la distribución. En muchas aplicaciones de la ingeniería de
petróleos, las distribuciones son de larga cola pero definida por muy pocos valores. Como
resultado de esto, la distribución de probabilidad es vagamente conocida en su cola por lo
cual debe introducirse cierta información ganada de la experiencia o de las condiciones que
rodean a la consecución de los datos. Por ejemplo, en casos donde el control de pozos es
bueno, se puede realizar interpolación lineal entre cuantiles y así definir el comportamiento
de la distribución de probabilidad en sus extremos. Pero en casos de bajo control de pozo
la extrapolación se hace necesaria. Obviamente los modelos de extrapolación determinarán
si la distribución tiene una cola extendida o corta. Es entonces en este punto donde el
criterio de quien realiza la simulación entra en juego. Por ejemplo, en el caso de estimación
de reservas se podría preferir el uso de extrapoladores que produzcan una cola
conservativamente corta para no incurrir en una exagerada sobrestimación de reservas.
Otro caso similar se presenta cuando valores extremos de permeabilidad generan problemas
durante una simulación de recobro secundario, en esta situación una cola extensa es
apropiada para generar todos los escenarios posibles.
Sin importar el método escogido para calcular la función de probabilidad, la forma, la
media y la varianza de la distribución dependerán del modelo de continuidad adoptado, la
localización de los puntos vecinos y los valores de la variable en estos puntos.
Una vez determinada la función de probabilidad, el valor estimado será tomado de esta
40distribución a través de simulación montecarlo11. El rango de datos que puedan ser
seleccionados será amplio en áreas de bajo control de pozo mientras que el rango será
estrecho en las zonas de buen control. Los resultados de repetidas simulaciones de un área
diferirán más, una de otra, en aquellas zonas donde se presenta bajo control.
2.7 FUNCIONES DE PÉRDIDA O ERROR.
Dados los muchos posibles valores de una propiedad simulada geoestadísticamente, la
pregunta obvia es cual de todas las realizaciones generadas debe ser escogida como la
mejor estimación de la variable en cuestión. Pues bien, muchas personas encogerían un
valor medio de todas las realizaciones; otros preferirían la mediana o incluso la moda de los
datos generados. Con tantos candidatos diferentes para mejor estimado otras técnicas de
elección son preferibles. Por ejemplo, el desarrollo de simulación de flujo en cada
escenario generado con el fin de obtener el mejor ajuste a la historia de producción del
campo. Por otro lado, desde un punto de vista de evaluación económica, el concepto de
función de perdida12 es aplicable al momento de elegir una realización como la mejor. Tal
concepto es explicado a continuación.
Siempre que alguna decisión esta basada en un valor estimado más que en el valor
verdadero, se incurre en sobrestimación o subestimación que finalmente generarán un
castigo o pérdida generalmente cuantificado en dinero. Esta pérdida representa la
diferencia entre la ganancia o pérdida al tomar decisiones usando un valor estimado y no en
el valor verdadero. La pérdida asociada a la estimación de estas propiedades puede ser
vista como el castigo de no poseer información exacta de las variables en cuestión. Es
decir, la perdida es función de un error.
Nótese que una función de pérdida puede presentar asimetría; es decir, la penalización por
sobrestimación no es igual a la penalización por subestimación. Otra característica
importante de las funciones de error es que dependen no de aquello que esta siendo
estimado sino de aquello para lo que serán usadas. Por ejemplo, dos yacimientos evaluados
para recobro mejorado, uno para inyección de agua y otro para inyección de gas, tendrán
diferente función de pérdida sin importar que la variable estimada (volumen poroso
interconectado) sea la misma. Esto es claro pues un error en la cantidad de gas inyectado
tendrá un costo diferente al costo que presentaría un error en la cantidad de inyección de
agua.
El concepto de función de error es importante porque provee un criterio de minimización de
estimadores. El interés se centraría en saber si el valor estimado, tratado como el
verdadero, minimiza o no la pérdida.
Desafortunadamente, poco trabajo se ha hecho en la teoría de funciones de pérdida y
generalmente se aborda el problema como un asunto de minimización de la varianza de
41estimación, es decir, se asume que las funciones de error son cuadráticas incurriendo así en
el siguiente inconveniente: No importa la forma de la distribución de probabilidad de los
estimados, una función cuadrática será siempre minimizada por la media. Sin embargo, en
muchas situaciones prácticas, otro estimativo sería más razonable de usar. Por ejemplo, en
el diseño de una plataforma marina de perforación no se usaría el estimado de la carga
media sino uno de la carga máxima.
2.7.1 Decisiones Basadas en Funciones de Pérdida
Considérese el caso donde se está interesado en el análisis económico del recobro final
generado mediante una inyección de agua, sabiendo que el valor real del factor de recobro
permanecerá desconocido hasta que el campo sea definitivamente cerrado. Sin embargo, la
incertidumbre en el factor de recobro puede ser estimada mediante el siguiente
procedimiento:
• Generar L realizaciones geoestadísticas de las propiedades del yacimiento necesarias
para una simulación de flujo.
• Realizar simulación de flujo en cada una de los escenarios geoestadísticos generados.
• Cada uno de los factores de recobro, resultado de cada simulación de flujo, pueden ser
utilizados para construir el histograma y función de probabilidad acumulada del factor
de recobro.
Con esta información de la incertidumbre de la variable en cuestión más una función de
pérdida que cuantifique el impacto de cometer un error en la estimación del recobro, puede
tomarse una decisión respecto al éxito del proyecto de inyección y más aun, puede
seleccionarse la realización geoestadística más conveniente para futuras simulaciones.
Para ilustrar lo dicho en el párrafo anterior se presenta el siguiente ejemplo: en la Figura 20
se presentan el histograma y la distribución acumulada del petróleo recuperado en cada una
de las cien realizaciones generadas mediante simulación geoestadística. Una función de
pérdida asumida previamente aparece también en la Figura 21. Esta función de pérdida
cuantifica el impacto económico de sobrestimar o subestimar el recobro final en campo
simulado. Finalmente, la Figura 22 se obtiene a partir de la distribución acumulada y la
función de pérdida. Esta Figura representa la perdida asociada a cualquiera de los
estimativos de recobro hechos. Utilizando el concepto de función de pérdida, el mejor
estimativo es aquel que minimiza la pérdida asociada, es decir, la realización que de como
resultado un recobro de 800 MMbbl barriles es la óptima.
42FIGURA 19. Histograma de Petróleo Recuperado.
Tomado de Srivastava, “An application of Geostatistical Methods for Risk Analysis in Reservoir
Management”, SPE paper 20608, 1990.
FIGURA 20 . Función defrecuencia Acumulada de Petróleo Recuperado.
Tomado de Srivastava, “An application of Geostatistical Methods for Risk Analysis in Reservoir
Management”, SPE paper 20608, 1990.
43FIGURA 21. Función de Pérdida.
Tomado de Srivastava, “An application of Geostatistical Methods for Risk Analysis in Reservoir
Management”, SPE paper 20608, 1990.
FIGURA 22. Minimización de la función de pérdida
Tomado de Srivastava, “An application of Geostatistical Methods for Risk Analysis in Reservoir
Management”, SPE paper 20608, 1990.
443. DEFINICIÓN DEL PROBLEMA
Una de las ayudas mas significativas que puede ofrecer la herramienta geoestadística en los
estudios de ingeniería de yacimientos, es sin duda alguna la manera en que se puede
manejar la información que usualmente se trabaja en dicho tipo de estudios (permeabilidad,
porosidad, etc.). Para el caso particular de este trabajo, se cuenta con una cantidad
considerable de información petrofísica que involucra 241 pozos de un mismo campo
petrolífero. En particular, se encuentran disponibles las ubicaciones de los pozos
pertenecientes al campo, los datos de espesor y del tope de la formación productora de
petróleo, la información de porosidad y de permeabilidad, y algunos adicionales obtenidos
de estudios geofísicos (velocidad de propagación de ondas y densidad de la roca). Es
importante resaltar que los datos con los que se cuentan son puntuales, es decir, que para
241 puntos dentro del yacimiento se tiene dicha información.
El campo al que se le quiere encontrar el mapa de porosidad y de permeabilidad tiene un
área de 20,000*40,000 pies2, es decir unos 4,591 acres y los pozos se encuentran
espaciados aproximadamente 1000 pies uno del otro, es decir que se encuentran
distribuidos en áreas de 20 acres aproximadamente. Un mapa de la ubicación de los datos
dentro del campo se presenta a continuación en la Figura 23.
FIGURA 23. Ubicación de los pozos en el campo.
45Es ilógico pensar que las propiedades petrofísicas de un yacimiento son constantes en todos
los puntos que este contiene, y para la muestra basta solo con mirar con algún detalle los
datos a los que se están haciendo mención (Anexo B). Es aquí donde entra en juego la
técnica geoestadística que se ha venido mencionando a lo largo de este texto. Lo que se
pretende es dar un manejo apropiado a la información con la que se cuenta, para generar un
modelo geoestadístico que caracterice dicho yacimiento. A continuación se presentan dos
opciones que permiten la solución del problema planteado, es decir la generación de mapas
geoestadísticos de porosidad y permeabilidad a partir de la información que se tiene del
campo.
1. Encontrar un modelo único de las propiedades petrofísicas del yacimiento a través
de la técnica del kriging ordinario. Sería muy ambicioso pensar que un modelo de
variables aleatorias pudiese dar una respuesta única a un problema tan complejo que
requeriría de estudios más avanzados en los que se incluyeran análisis de covarianza
para analizar como interactúan unas variables con otras; por ejemplo, sería
necesario verificar los modelos obtenidos con un kriging ordinario para una variable
como la porosidad, comparando las respuestas con las obtenidas por el cokriging
entre la misma porosidad y un parámetro geofísico como la impedancia acústica.
Sin embargo, reacuérdese que la finalidad de este proyecto es mostrar una serie de
técnicas geoestadísticas que en una eventualidad puedan ser usadas en estudios de
caracterización de yacimientos, mas no se busca obtener un modelo exacto del
campo al que se hace referencia en párrafos anteriores. Por lo tanto, es prioridad
para este trabajo de grado mostrar los pasos fundamentales para desarrollar mapas
de propiedades petrofísicas obtenidos con kriging ordinario, haciendo uso de
programas geoestadísticos16, que se encuentran disponibles para el uso de quienes
pueden estar interesados en este tipo de estudios.
2. La simulación gaussiana da la posibilidad de obtener diferentes mapas para una
misma variable que de igual forma pueden ser los que representen dicho parámetro
en la formación a la que pertenecen, en otras palabras, por medio de dicha
simulación se podría encontrar una cantidad de modelos equiprobables que
identifiquen la variable que se este estudiando. Sin embargo, una propiedad no
puede estar representada por una cantidad indefinida de mapas, y por el contrario un
único mapa sería suficiente para representar la variable deseada dentro del
yacimiento o formación de interés.
Se plantea entonces la inquietud de definir una finalidad para la cantidad de modelos que se
pueden obtener a través de una simulación geoestadística. Hipotéticamente, las
realizaciones obtenidas podrían ser utilizadas en un simulador de yacimientos para acoplar
un modelo que permita realizar simulaciones de flujo. Pero entonces de igual forma se
tendría la posibilidad de elegir entre una cantidad indefinida de realizaciones para incluir en
los parámetros de entrada del simulador. Una posibilidad de selección de un modelo único
de parámetros petrofísicos, sería la de realizar un ajuste histórico de la producción del
campo con la ayuda de un simulador de flujo, utilizando todas las realizaciones que se
46obtienen de la simulación gaussiana; posteriormente se determinaría que combinación de
realizaciones es la que menos error genera con respecto a los datos históricos de
producción. Suponiendo que las variables petrofísicas son independientes, entonces se
debería realizar una cantidad de realizaciones tal que el trabajo de determinar el conjunto de
mejores realizaciones consumiría demasiado tiempo. Este trabajo podría realizarse, pero
para los fines de este proyecto se ha considerado que es un trabajo tedioso y además no es
de tanta relevancia determinar el verdadero modelo de parámetros petrofísicos que
representan el campo que se esta trabajando.
Por el contrario, se considerará en este trabajo, la manera en que se pueden utilizar todas las
realizaciones que se obtienen, en la búsqueda de información de utilidad para la evaluación
de un proyecto petrolífero determinado. Lo que se pretende es utilizar los conceptos de la
geoestadística para obtener las realizaciones de las propiedades petrofísicas y con estas
realizaciones poder estimar un rango de valores para variables más importantes y que
obligatoriamente se necesitan en la evaluación económica de un proyecto. Para explicar lo
anterior, basta con pensar en el hecho de que la gerencia de un proyecto siempre querrá
saber la cantidad de petróleo que se puede producir en un campo en un determinado tiempo,
pero ellos también querrán saber cuál es la confianza que tiene ese valor. Habitualmente, en
la evaluación de proyectos, se requiere conocer cuál sería la máxima producción esperada,
la mínima producción esperada, y cuál la más probable. El problema entonces será plantear
una metodología que permita determinar dichos rangos para variables necesarias en una
evaluación de proyectos (Petróleo producido, tasas de producción, cortes de agua, tasas de
inyección de agua o de gas si las hay, etc.)
Será prioridad de este trabajo mostrar como se obtienen las realizaciones de los parámetros
petrofísicos por medio de la simulación geoestadística gaussiana. Además, ahora también
será prioritario desarrollar una metodología que permita obtener rangos de los valores de
una propiedad de interés utilizada en la evaluación de proyectos de petróleo.
474. METODOLOGÍA DE TRABAJO
Ya se han definido los problemas que se pretenden atacar en este proyecto. Ahora, se
plantea la metodología de trabajo a seguir para cumplir con los objetivos que se han
planteado. La primera opción de generación del modelo geoestadístico, es quizás la mas
simple de las dos y la que históricamente se ha utilizado para obtener dichos modelos. La
metodología para obtener un modelo de parámetro petrofísico mediante el uso del kriging
ordinario se presenta a continuación.
FIGURA 24. Diagrama de flujo para obtención de modelo geoestadístico usando kriging
ordinario.
La Figura 24 puede ser explicada de la siguiente manera:
El primer paso que se debe seguir es determinar por medio de la variografía cual es el
variograma que modela la correlación espacial que dicha variable puede tener dentro del
área de estudio. Para poder obtener un buen variograma se hará uso de programas
geoestadísticos, que en primer lugar permiten analizar si existen anisotropías
geoestadísticas, y si tales anisotropías existen, entonces obtener los variogramas en la
mínima y en la máxima dirección de anisotropía. Recuérdese que isotropía geoestadística,
48no significa que los valores de la variable son constantes dentro del área de estudio, lo que
significa es que los cambios de dicha variable se pueden apreciar de una manera constante
en cualquiera de las direcciones en las que se estudie el fenómeno; por lo tanto anisotropía
geoestadística significa que los valores de la variable de estudio cambian de diferente
manera a medida que se varia la dirección en la que se analiza el variograma. En caso de
que la variable sea isotrópica estadísticamente hablando, bastará con un variograma para
observar su correlación espacial. El variograma (o variogramas si la variable es
anisotrópica) puede ser obtenido utilizando el programa GSLIB. En dicho paquete también
se puede obtener el modelo respectivo para el variograma, cuyas ecuaciones serán
utilizadas en la matriz del kriging ordinario.
El segundo paso será utilizar las ecuaciones obtenidas en el modelo del variograma, para
alimentar la matriz del kriging ordinario. Dichas ecuaciones pueden ser ingresadas por
medio del GSLIB; el mismo programa realiza tanto la solución de las matrices para cada
nodo dentro de la malla como un gráfico en el que esquemáticamente se puede apreciar la
solución del problema. Debido a que las matrices del kriging son completamente lineales la
solución del sistema no tomará demasiado tiempo computacional y por lo tanto la malla que
se puede utilizar en la obtención del modelo geoestadístico por kriging puede ser tan
refinada como se quiera.
También se definió en el capítulo anterior una opción para modelar geoestadísticamente un
yacimiento mediante el uso de la simulación gaussiana. Se ha dicho que no basta con
obtener múltiples realizaciones de un parámetro petrofísico, ya que carecería de sentido
entregar como resultado muchos modelos geoestadísticos sin decidir cual es el definitivo.
La idea planteada era la generación rangos de valores para las variables que se incluyen en
la evaluación económica de un proyecto. La segunda opción se puede esquematizar con la
ayuda del diagrama de flujo de la Figura 25.
La Figura 25 puede ser explicada de la siguiente manera.
Los resultados obtenidos de la variografía y el kriging serán aplicados para alimentar un
archivo de datos de entrada en el programa GSLIB, con el fin de realizar la simulación
geoestadística gaussiana. En el marco teórico de la simulación gaussiana se ha dicho que se
podrán obtener múltiples realizaciones, es decir múltiples modelos de una misma variable,
dependiendo de un número aleatorio que definirá el recorrido por medio del cual se irán
obteniendo los valores particulares de dicha variable de estudio en cada una de las celdas
del mallado que contiene el área de estudio. Para fines prácticos se decidió generar cien
realizaciones para las variables porosidad y permeabilidad.
Utilizando los cien modelos de porosidad y permeabilidad, se realizarán igual número de
simulaciones de flujo, con las que se pretende modelar un proceso de inyección de agua.
Tales simulaciones se plantean hacer con un simulador de flujo de yacimientos del tipo
stream-line14, el cual permitirá realizar dicho trabajo con un menor tiempo computacional
que el que se gastaría utilizando un simulador de yacimientos convencional. Los
49parámetros de salida del simulador stream-line incluyen el petróleo producido durante el
periodo de inyección de agua en el yacimiento en un tiempo determinado.
FIGURA 25. Diagrama de flujo para obtención de modelo geoestadístico usando
simulación geoestadística.
La cantidad de simulaciones de flujo que se plantean realizar implica que el número de
valores de petróleo producido durante un tiempo determinado sea el mismo. Los cien
valores de petróleo producido serán utilizados para realizar un histograma de frecuencias
que permita visualizar las probabilidades y los rangos en que dicha producción se puede
encontrar. El histograma que se plantea realizar, es una herramienta fundamental en la
evaluación de proyectos de petróleo para la toma de decisiones de inversión en la industria,
ya que este ilustra el valor más esperado como también los escenarios más pesimistas y
optimistas de producción de petróleo. Con esta técnica no se obtendrá un modelo único de
porosidad o permeabilidad, pero si se obtendrán resultados útiles para la toma de decisiones
económicas, que en últimas instancias es lo que define si un proyecto petrolífero se lleva a
cabo o si por el contrario no se hace.
505. RESULTADOS
A continuación se presentarán los resultados más relevantes que se obtuvieron siguiendo
la metodología propuesta para solucionar el problema de obtener un modelo de porosidad
y permeabilidad.
Inicialmente se presenta un análisis exploratorio de los datos. Este consistirá
exclusivamente de una prueba de normalidad que permita aplicar el algoritmo de
simulación gaussiana expuesto e incluso detectar la presencia de mas de una familia en el
grupo de datos pues tanto este como el área de estudio es bastante grande.
El mismo GSLIB tiene la opción de trabajar los datos de la muestra con la estadística
clásica. Esta función se utilizó para encontrar los histogramas de ambas variables. Se
encontró que la porosidad se distribuye de una formal normal y que la permeabilidad se
distribuye log-normal (logaritmo de permeabilidad se distribuye normal). Los
histogramas para porosidad y para log-permeabilidad, se muestran a en las Figuras 26 y
27 respectivamente.
FIGURA 26. Histograma de frecuencias de Porosidad.
51FIGURA 27. Histograma de frecuencias de Log - Permeabilidad.
5.1 MAPAS DE VARIOGRAMA
El primer paso en el análisis geoestadístico es desarrollar mapas de variograma para
determinar si una propiedad presenta anisotropía geoestadística o no. Ya se ha definido,
que existirá anisotropía geoestadística si los variogramas que se elaboran en diferentes
direcciones presentan diferencias sustanciales en sus propiedades. Un mapa de
variogramas muestra en un plano cartesiano como se comportan todos los variogramas en
una cobertura de 3600. Teniendo en cuenta que el variograma que representa la
correlación espacial de la dirección norte es el mismo que representa a la dirección sur, y
así con todos los demás y sus correspondientes en las direcciones opuestas, entonces el
mapa de variogramas será simétrico con cualquier eje que se tome como referencia en el
cuadrante cartesiano. Por lo tanto el mapa de variogramas parece ser la herramienta
indicada para mostrar si tales anisotropías existen o si por el contrario la variación de una
propiedad es la misma en cualquier dirección.
Para realizar los mapas de variograma se usó el software GS+. Los datos de entrada que
se necesitan para este programa son simplemente las coordenadas de las ubicaciones de
los datos y los valores de las propiedades en dichos puntos. Solo esta información basta
para obtener una primera aproximación de un mapa de variograma, sin embargo, se
52puede realizar una mejor aproximación manipulando las magnitudes del vector h de
búsqueda y de la tolerancia que se le quiere dar a dicho vector. Se ha encontrado como
regla general que cuando el valor de la tolerancia se hace igual a la magnitud del vector h,
se pueden suavizar los picos que habitualmente se presentan en los variogramas que se
obtienen con diferentes tolerancias. Sin embargo, para encontrar el valor del vector h
adecuado, es necesario probar diferentes valores hasta que a juicio personal se crea que se
ha encontrado el variograma o mapa de variograma adecuado.
Ya se ha definido que se modelaran en este proyecto las propiedades de porosidad y
permeabilidad. Los mapas de variograma que se obtuvieron mediante el uso del GS+ se
muestran a continuación en las Figuras 28 a 31. Es importante señalar que los valores que
se usaron para los vectores h de búsqueda y para la tolerancia del vector fueron 1522 y
2000 pies para porosidad y permeabilidad respectivamente. Básicamente lo que se
encontró en dichos mapas, es que ninguna de las propiedades presentó anisotropía
geoestadística. Un análisis más detallado de la razón por la cual se consideró que las
propiedades porosidad y permeabilidad son isotrópicas geoestadísticamente se presentará
en el siguiente capítulo. Por ahora solo basta saber que se puede considerar un
variograma omnidireccional para modelar la continuidad espacial tanto de porosidad
como de permeabilidad.
FIGURA 28. Mapa de variograma 2D para porosidad.
53FIGURA 29. Mapa de variograma 3D para porosidad.
FIGURA 30. Mapa de variograma 2D para permeabilidad.
54FIGURA 31. Mapa de variograma 3D para permeabilidad
5.2 VARIOGRAFÍA
Para elaborar un variograma omnidireccional basta con considerar que la tolerancia del
ángulo de búsqueda del variograma es igual a 900. Esto significa, que la búsqueda de
datos para calcular el variograma se hace sin tener en cuenta una dirección especifica;
solo se tiene en cuenta la distancia definida por el vector h y por el valor de tolerancia de
dicho vector.
Utilizando el GSLIB se puede obtener el omnivariograma. Los datos de entrada que se
utilizaron fueron simplemente las coordenadas de la ubicación de los datos, los datos
mismos, definir la tolerancia del ángulo de búsqueda como 900 y utilizar los valores de
magnitudes del vector de búsqueda y de la tolerancia que se obtuvieron con los mapas de
variograma (1522 pies para porosidad y 2000 para permeabilidad). Después de obtener el
variograma, se utilizó el mismo paquete para encontrar el mejor modelo que lo
represente; se dice que es el mejor porque es el que menos error presenta entre los valores
del variograma y los valores del modelo. Los variogramas para porosidad y
permeabilidad, como también los modelos que los representan se pueden apreciar en las
Figuras 32 y 33.
55FIGURA 32. Modelo ajustado del variograma de porosidad.
FIGURA 33. Modelo ajustado del variograma de permeabilidad.
56La ecuación 58 representa la gráfica de un variograma del tipo esférico con nugget. El sill
en este caso no es igual al valor de C, sino que es igual al valor de C más el valor del
nugget (Ngt). La ecuación indica que cuando h es igual a 0 entonces el variograma tendrá
un valor igual al del nugget, y cuando h sea mayor que el valor del rango a entonces el
valor del variograma siempre será igual al del sill. La ecuación también muestra que el
variograma se comporta linealmente para valores pequeños de h. Los valores de los
parámetros de la ecuación 58 para las variables porosidad y permeabilidad se muestran en
la tabla 2.
TABLA 2. Parámetros de las ecuaciones de los modelos de variograma.
Parámetro de la ecuación (58) Variograma Porosidad Variograma Permeabilidad
C (Semi-sill) 0.651 1
a (Rango) 9140 18000
Ngt (Efecto Nugget) 0.389 0.3
Es importante mencionar que los variogramas que se obtuvieron con el GSLIB se
encuentran estandarizados, ya que el valor real del variograma se divide por la varianza
de la muestra de datos, por eso los valores de variograma no son muy grandes y como se
puede ver oscilan entre 0.389 y 1.04 para porosidad y 0.3 y 1.3 para permeabilidad. Lo
anterior se hizo de esta manera ya que los datos de entrada que necesita la opción kriging
del GSLIB requieren tener los parámetros de un variograma estandarizado. El programa
realiza internamente los cálculos de estandarización, por eso no se considera importante
revelar el valor de la varianza en este momento. Más adelante se presentaran unos
estadísticos clásicos (media, varianza, etc.) de las muestras de porosidad y permeabilidad,
que son necesarios conocer para el desarrollo de la simulación geoestadística gaussiana.
5.3 KRIGING
La elaboración del kriging se realiza con la opción de kriging para GSLIB. En los datos
de entrada de esta opción, se requiere entrar los parámetros del variograma que se
obtuvieron en el paso anterior. Esto con el fin de que el programa calcule los pesos
apropiados w para estimar los valores de la variable deseada en los bloques de la malla
i
donde no se tiene un dato.
Recordando que el campo tiene 20,000 pies en la dirección este-oeste y 40,000 pies en la
dirección norte-sur se decidió elaborar una malla con bloques de 100*100 pies2. Es decir
que en total se trabajaron 80,000 bloques. Como consecuencia de que las ecuaciones del
sistema de kriging son lineales, el tiempo computacional que tarda el GSLIB en elaborar
un mapa de kriging es poco significativo, por lo tanto la malla se puede hacer tan fina
como se desee sin tener problema alguno. Sin embargo, si se escogiesen los modelos
geoestadísticos de porosidad y permeabilidad que se obtienen del kriging para alimentar
un simulador de yacimientos, se tendría que evaluar cual seria el tamaño de malla
57adecuada para trabajar. En particular, es bien conocido que el tiempo computacional de
un simulador de flujo de yacimientos es significativo, y que entre más refinada sea la
malla, más tiempo se estará consumiendo, por lo tanto el simulador seria quien defina el
tamaño de la malla que se debería realizar con kriging. Una malla con 80,000 bloques
podría ser demasiado grande para ser ingresada en un simulador de flujo; sin embargo a
manera de solo ilustración se muestra este kriging pues este no se usará en la simulación
de flujo que se plantea realizar con los modelos obtenidos a partir de simulación
geoestadística.
FIGURA 34. Mapa de Porosidad elaborado con kriging.
El resultado del kriging para porosidad y permeabilidad que se obtuvo con el GSLIB para
la malla descrita anteriormente se muestra en las Figuras 34 y 35 respectivamente. Un
análisis de las gráficas de kriging se muestra en el siguiente capitulo de este texto.
58FIGURA 35. Mapa de permeabilidad elaborado con kriging.
5.4 SIMULACIÓN GEOESTADÍSTICA
Para poder elaborar la simulación geoestadística gaussiana es necesario que las variables
que se van a estudiar se distribuyan normalmente y esto ya se demostró en la sección de
análisis exploratorio de los datos.
Ahora que se sabe la forma en que se distribuyen las muestras de las variables porosidad
y permeabilidad y de que se puede realizar la simulación geoestadística gaussiana, se
puede describir la manera de encontrar las diferentes realizaciones. Primero es necesario
encontrar la distribución normal estandarizada (media igual 0 y varianza igual a 1) ya que
esta es la forma en que trabaja el GSLIB. La función “Normal scores” permite encontrar
los valores normalizados de una distribución gaussiana, por lo tanto esta función se
utilizó para normalizar los datos de porosidad y de log-permeabilidad. El archivo de datos
normales es el que contiene los datos de entrada para hacer las diferentes simulaciones
geoestadísticas y encontrar las realizaciones que se requieren para entrar en el simulador
de flujo del tipo stream-line. El archivo de datos de entrada para hacer la simulación
geoestadística es prácticamente el mismo para todas las realizaciones, lo único que difiere
es la selección de un número aleatorio que define el recorrido que hará el simulador para
ir encontrando los estimados en todos los bloques de la malla.
59La malla que se seleccionó para elaborar las simulaciones geoestadísticas es de 350*350
pies2, es decir que comprende 57 bloques en la dirección este-oeste y 114 en la dirección
norte-sur (6,318 bloques en total). Nótese que la malla que se plantea trabajar para la
simulación geoestadística no es tan refinada como lo era la malla que se utilizó en el
kriging ordinario. Esto porque las realizaciones que se obtengan de las simulaciones
serán utilizadas en el simulador de flujo y por lo tanto una malla de 80,000 bloques seria
tediosa de simular por efectos de tiempo computacional, mucho más teniendo en cuenta
que por tener 100 realizaciones se requiere de 100 simulaciones para desarrollar la
metodología que se ha planteado. El archivo de datos de entrada para la simulación
geoestadística incluye los datos de porosidad normalizada o log-permeabilidad
normalizada según sea el caso, junto con la ubicación correspondiente de dichos datos, un
número aleatorio cualquiera que definirá la ruta de búsqueda de los estimados, y los
parámetros del variograma obtenidos anteriormente. Es importante hacer saber que el
variograma de la propiedad también representa a la propiedad normalizada y también al
logaritmo de la propiedad normalizada, ya que este no depende del valor de la propiedad
sino de la forma en que se correlacionen los datos de la muestra dentro del campo, y por
lo tanto el mismo variograma puede estar representando a la propiedad y a una de sus
transformadas. Con el mismo número aleatorio que define la ruta de la estimación se han
obtenido dos realizaciones, una para porosidad y otra permeabilidad. Las realizaciones
que se obtienen con el archivo de datos del que se habla en el párrafo anterior son
realizaciones de la transformada de las variables estudiadas. Por lo tanto es necesario
realizar una transformación inversa para obtener la realización deseada con las unidades
correspondientes. En el caso de la permeabilidad es necesario obtener además de la
transformada inversa el antilogaritmo para llevarlo a las unidades de mD. La
transformación inversa se puede realizar de una forma simple con el comando “Back-
Transform” de GSLIB.
Una pareja de realizaciones de las dos variables se utilizará en la elaboración de un
modelo de yacimiento para simular un proceso de inyección de agua. Este proceso se
llevo a cabo 100 veces. A continuación se presentan 5 parejas de modelos de porosidad y
permeabilidad de las 100 que se obtuvieron en las Figuras 36 a 40.
60FIGURA 36. Simulaciones geoestadísticas de porosidad y permeabilidad.
61FIGURA 37. Simulaciones geoestadísticas de porosidad y permeabilidad.
62FIGURA 38. Simulaciones geoestadisticas de porosidad y permeabilidad.
63FIGURA 39. Simulaciones geoestadísticas de porosidad y permeabilidad.
64FIGURA 40. Simulaciones geoestadísticas de porosidad y permeabilidad.
655.5 SIMULACIÓN DE FLUJO TIPO “STREAM-LINE”
Después de que fueron obtenidas las 100 realizaciones de porosidad y permeabilidad con
el procedimiento de simulación geoestadístico, se procedió a simular un proceso de
inyección de agua para el campo estudiado. El objetivo fundamental de este trabajo se ha
cumplido hasta este momento: estudiar una técnica para generar datos de porosidad y
permeabilidad en donde no se tienen, usando la información de las mismas propiedades
que se encontraban disponibles en otros puntos del mismo campo. Por lo tanto el modelo
de simulación de flujo que se realizó es muy simple y de carácter demostrativo, por lo
tanto está basado en hipótesis que posiblemente no corresponden a las mismas
características que se pueden presentar en dicho campo, exceptuando las propiedades
petrofísicas que se estimaron con el proceso geoestadístico.
FIGURA 41. Arreglo de pozos para inyección de agua.
66Una de las primeras hipótesis que tiene para este ejemplo de simulación es la selección
del arreglo de pozos para la inyección de agua. Debido a la simetría de las líneas de los
pozos en la dirección vertical, se planteo un arreglo de línea directa que alterna inyectores
en una línea horizontal y productores en la siguiente línea. Un esquema del arreglo de
inyección se presenta en la Figura 41, en donde los pozos inyectores están representados
como puntos azules y los productores con puntos rojos.
Para realizar las simulaciones de flujo se utilizo el programa 3DSL (Simulador tipo
stream-line). Dentro de los parámetros que se requieren para el archivo de entrada se
encuentran las propiedades petrofísicas, que como se ha dicho serán tomadas del
resultado de la simulación geoestadística. Otros parámetros necesarios para la simulación
son las propiedades de los fluidos involucrados en el proceso, la presión del yacimiento al
tiempo inicial de simulación, las curvas de permeabilidad relativa, las características de
los pozos inyectores y productores entre otros.
La malla que se utilizó para la simulación stream-line fue de 57 bloques en la dirección
este-oeste y 114 bloques en la dirección norte-sur (modelo 2D). Se consideró un petróleo
tipo “Black oil” inmiscible con el agua; las propiedades usadas para el petróleo y el agua
se pueden ver en la tabla 3.
TABLA 3. Propiedades del petróleo y el agua usados para simulación de flujo.
Propiedad Petróleo Agua
Viscosidad (cp) 2.5 1.5
Densidad (lbm/pie3) 45 63
Saturación Irreducible (%) 0.2 0.22
Saturación Inicial (%) 0.78 0.22
Se consideró que la permeabilidad en la dirección x era igual a la permeabilidad en la
dirección y. Se aplicaron restricciones tanto a la presión de fondo de los pozos
productores como a la presión de inyección que determinaron las tasas de producción y
de inyección de cada uno de los pozos. Para los pozos productores se uso una restricción
de 200 psi de presión de fondo de pozo, esto significa que si la presión de fondo de pozo
bajaba de 200 psi entonces el programa daría la orden para que este se cerrara
automáticamente, de otro lado para pozos productores la restricción se fijó en 2,000 psi,
en este caso la restricción de cierre aplicaría para presiones mayores de este valor. La
presión del campo en el tiempo cero de simulación fue fijada en 3,600 psi. Todos los 241
pozos se consideraron verticales y con un diámetro de 6 pulgadas para la tubería de
producción. El tiempo de inyección simulado en el campo fue de 15 años con un paso de
tiempo de 75 días.
67El resultado obtenido de la simulación de flujo que se consideró para entrar en este
estudio fue el petróleo final recuperado al cabo de los 15 años de producción. Ya que se
realizaron 100 simulaciones de flujo, por lo tanto se obtuvieron 100 datos de petróleo
final recuperado con los cuales se construyo el histograma de frecuencias de la Figura 42
y el histograma de frecuencias acumuladas de la Figura 43. Los parámetros estadísticos
más importantes de los histogramas se presentan en la tabla 4. Los valores de petróleo
final recuperado con los que se elaboraron los histogramas se encuentran en la tabla del
Anexo C. Estos valores son el resultado final de este estudio y se analizaran en el
siguiente capitulo, como también se hará mención a su utilidad en la ingeniera de
petróleos.
FIGURA 42. Histograma de frecuencias del petróleo final recuperado a 15 años.
68FIGURA 43. Histograma de frecuencia acumulada del petróleo final recuperado a 15 años.
TABLA 4. Parámetros estadísticos del histograma de frecuencia del petróleo final recuperado.
Parámetro Valor
Media (MMSTB) 50.5
Desviación estándar (MMSTB) 1.69
Valor mínimo (MMSTB) 47.1
Valor máximo (MMSTB) 54.4
696. ANÁLISIS DE RESULTADOS
6.1 SUPERFICIES VARIOGRÁFICAS
En las Figuras 28 a 31 se presentan las superficies variográficas correspondientes a los
datos de porosidad y permeabilidad. En la sección teórica de este trabajo se mencionó que
la utilidad de estos mapas radica en que en ellos es fácilmente detectable una dirección de
anisotropía. Pues bien, observando ambos mapas de variograma 2D (Figuras 28 y 30) es
notable un contorno circular que crece desde el centro del mapa hacia afuera. Tal contorno
conserva una forma aproximadamente circular hasta distancias de 9000 pies para el mapa
de permeabilidad y 7000 pies para el de porosidad. Esta es la característica típica de una
superficie variográfica generada por datos isotrópicos geoestadísticamente hablando, es
decir, no importa la dirección en que se comparen un par de puntos, el rango del
variograma es siempre el mismo. Sin embargo, si se presta atención al comportamiento de
las superficies variográficas más alla de las distancias anotadas, se podría estar tentado a
sugerir cierto tipo de anisotropía pues la altura del variograma es variable. Pero
considérense un los siguientes aspectos: 1) Recuérdese que los datos disponibles están
regularmente espaciados a una distancia media de 1000 pies así que para el caso más
pesimista (porosidad) se tiene un radio de 9000 pies en el que se puede utilizar 9 pozos
para estimar un nuevo valor. Por otro lado, en un área de alta densidad de pozos como la
estudiada, cuando se estima un nuevo valor se hace innecesario el uso de datos a distancias
lejanas del punto que va a ser estimado.
Así pues, a partir de las superficies variográficas queda determinado que el comportamiento
de ambos grupos de datos, porosidad y permeabilidad, es isotrópico. De esta conclusión
queda también establecido que el modelamiento del variograma se hará basado en un
cálculo omnidireccional del mismo.
6.2 VARIOGRAFÍA
Para los variogramas que se obtuvieron de las variables porosidad y permeabilidad, se
ajustaron modelos esféricos con efecto nugget. Tal vez estos modelos no describen el
comportamiento total del variograma en todas las distancias para las que fue calculado,
pero si lo hacen en forma muy exacta en la zona cercana al origen que es la zona de mayor
importancia al momento de estimar nuevos valores de las variables en cuestión. Además,
aplicando el mismo concepto de que puntos muy distantes son innecesarios a la hora de
estimar, se ha optado por no recurrir a un modelo con comportamiento oscilatorio como los
mostrados en las Figuras 10 a 12.
Combinados los resultados de los variogramas y las superficies variográficas, quedan
70determinadas las siguientes áreas de búsqueda para las variables en cuestión:
Porosidad: área circular definida por un radio de 9140 pies
Permeabilidad: área circular definida por un radio de 18000 pies.
Puntos que estén separados a una distancia mayor de dichos valores no serán tenidos en
cuenta para calcular los estimados con las ecuaciones de kriging o con el proceso de
simulación geoestadística.
6.3 KRIGING Y SIMULACIÓN GEOESTADISTICA
En la Figuras 34 y 35 (Mapas de kriging) es clara la tendencia suavizada de la variable
estimada.
FIGURA 44. Mapa de localización de datos de porosidad.
71Comparando con el mapa de localización de los datos utilizados para generar tanto el mapa
de porosidades como permeabilidades (Figuras 34 y 35), es notable como kriging genera
valores con bajo contraste entre dos puntos ya conocidos y es precisamente este fenómeno
el que conlleva a recurrir a nuevos métodos que generen más contraste entre valores
estimados y medidos pero que aun cumplan con características importantes como el
histograma original del grupo de datos (misma media y misma varianza), la correlación
espacial original (el variograma original) y los datos originalmente recolectados.
FIGURA 45. Mapa de localización de datos de permeabilidad.
Vale la pena anotar que especialmente en variables como la permeabilidad esta
característica de suavidad en el mapa generado no es representativa del fenómeno
estudiado, pues en un yacimiento donde son comunes las diferencias litológicas,
tortuosidades de canales en la roca, porosidad no efectiva, etc. las condiciones son mas
propicias para encontrar mayor variación en la propiedad a pequeña escala. En conclusión,
72una simulación de flujo ha de esperarse mas acertada en un campo de permeabilidades que
presente mayor contraste que en uno suavizado como el generado mediante kriging.
Pasando ahora a los mapas generados mediante simulación geoestadística (Figuras 36 a 40)
se observa como la característica de suavidad ha desaparecido y en cambio se tienen mapas
que presentan altos contrastes de valores. Incluso en bloques vecinos los valores ya no
parecen estar tan altamente relacionados como lo sugerían los mapas generados mediante
kriging. No obstante, si bien se ha mejorado en la representatividad del fenómeno, ahora se
tiene que enfrentar un problema de unicidad. Kriging proporcionaba un modelo de
permeabilidad y porosidad único, que a pesar de ser suavizado cumplía con todas las
restricciones impuestas (histograma, variograma, datos originales), en cambio ahora con la
simulación geoestadística se tienen múltiples realizaciones de porosidad y permeabilidad
que cumplen con las mencionadas restricciones. Así pues, la pregunta obvia que se genera
es cuál de todos estos modelos es que realmente describe la variable analizada.
6.4 APLICACIONES
Para el tratamiento de toda la información generada con la simulación geoestadística se
tiene una de las siguientes posibilidades: elegir el mapa de estimativos que mejor represente
el yacimiento, o bien utilizar todas las realizaciones como fuente de información.
En el primer caso varias son las vías de análisis.
• Simulación de flujo en cada escenario generado y posterior ajuste histórico de la
producción. El escenario que proporcione un mejor ajuste será el indicado a ser elegido
como mejor representación del yacimiento.
• Uso de funciones de pérdida.
• Corregionalización.
• Métodos de estimación que aplican como restricción el ajuste de otra variable conocida
en el área de estudio como atributos sísmicos, por ejemplo (Kriging colocado)
En el segundo caso, la aplicación típica se centra en la generación de distribuciones de
probabilidad para cada una de las propiedades derivadas o aquellas que sean influidas por
las variables estimadas. Por ejemplo las respuestas de flujo de un simulador o cálculos
volumétricos de petróleo original in situ o el volumen de petróleo recuperado mediante un
método de recobro mejorado.
Como fue manifestado al principio de este texto, debido al extenso trabajo que demandaría
hacer un ajuste histórico de datos, está fuera del alcance de este trabajo decidir cual de las
realizaciones generadas es la mejor representación del área estudiada. Sin embargo, para
finalizar y por ser problemas de menor complejidad, a continuación se citarán algunos de
73los ejemplos de aplicación de una distribución de probabilidad obtenida vía simulación
geoestadística, es decir, el histograma de frecuencias de petróleo recuperado que se obtuvo
después de usar las realizaciones geoestadísticas y realizar las simulaciones de flujo.
El histograma de frecuencias que se obtuvo permite visualizar las probabilidades de
ocurrencia de diferentes volúmenes de petróleo recuperado. Basados en estas
probabilidades se pueden generar varios escenarios para realizar el estudio económico del
proyecto, los cuales pueden incluir un escenario pesimista uno optimista y por un último un
caso más probable de producción que corresponde al valor medio del histograma. En este
trabajo se mostró como obtener el histograma de frecuencias de una sola variable de interés
en el análisis económico de la evaluación de un proyecto de inyección de agua para
recuperación de petróleo, de igual manera se pueden obtener las funciones de probabilidad
de diferentes parámetros que deben ser tenidos en cuenta en dicho análisis económico como
por ejemplo la producción de agua acumulada en el tiempo del proyecto. Este mismo
procedimiento se puede realizar para estudiar otro proceso de recuperación diferente al de
inyección de agua para comparar los resultados y decidir cual sería el método de recobro
mas conveniente desde un punto de vista económico, o también analizar diferentes patrones
de inyección diferentes al de línea directa que se utilizó en las simulaciones de flujo de este
trabajo.
6.4.1 Evaluación de un Proyecto de “Infill Drilling”
Para visualizar como entrar a jugar las probabilidades de ocurrencia de eventos en un
análisis económico, se presenta el siguiente ejemplo ilustrativo de un proyecto de “infill
drilling”
En primer lugar debe notarse que todo proyecto antes de ser ejecutado debe ser revisado en
el departamento económico y pasar los análisis necesarios para ser declarado
económicamente beneficioso. Segundo, es claro que cualquier tipo de información tanto
adquirida mediante medición como estimada tiene un grado de incertidumbre inherente.
En un proceso típico de evaluación de proyectos interviene todo tipo de información
estadística disponible con el fin de simular las posibilidades de éxito que podrían tenerse. Y
es en este punto donde la geoestadística cobra importancia permitiendo producir
estimativos y funciones de probabilidad que combinen diferentes tipos de información
(sísmica, petrofísica, geológica y de producción) y al mismo tiempo reduzcan los niveles de
incertidumbre.
Para exponer un ejemplo sencillo, considérese un proyecto de “infill drilling” en el cual la
caracterización de las zonas no barridas en el yacimiento es crítica15. Usualmente un
proyecto económico es analizado siguiendo un esquema de eventos y sus respectivas
probabilidades como el mostrado en la Figura 46. En este esquema están anotados una
serie de eventos y consecuencias cada uno de ellos con una respectiva probabilidad de
74ocurrencia. Al final del mismo se han calculado una serie de probabilidades de éxito para
tres eventos A, B y C diferenciados por la cantidad de petróleo recuperable después de la
ejecución del proyecto de “infill drilling”. Como se puede ver, el éxito de encontrar
reservas recuperables es función de las probabilidades de ocurrencia de cada uno de los
eventos que preceden. Así pues, la probabilidad de éxito aumentará si por alguna vía es
posible aumentar dichas probabilidades. La vía que se propone en este trabajo es la
simulación geoestadística. Pero no se debe pensar que el único objetivo es aumentar
probabilidades, en definitiva se trata de disminuir incertidumbre y aumentar confiabilidad,
de hecho las probabilidades de las que se ha venido hablando pudrían incluso disminuir,
pero lo importante es que se esta teniendo en cuenta la probabilidad correcta en el análisis,
y obviamente es mejor usar un valor de probabilidad bajo de alta confiabilidad que un valor
sobreestimado que conlleve a tomar decisiones equivocadas. En conclusión, el objetivo de
las distribuciones de probabilidad generadas vía geoestadística es, sin importar la decisión
final, que esta sea la más correcta.
75FIGURA 46. Árbol de Decisiones para un Proyecto de “Infill-Drilling”
767. CONCLUSIONES
• En el proceso de simulación geoestadística se obtiene una representación con iguales
características de variabilidad y correlación espacial que la realidad estudiada. Este
proceso se concentra fundamentalmente en la generación de valores a partir de una
función que describe la correlación espacial del fenómeno, obtenida del análisis de la
información experimental disponible, valores que posteriormente son condicionados a
dicha información. La simulación geoestadística, a diferencia de la estimación que
proporciona una representación suave o lisa de la realidad, es considerada como una de
las posibles realizaciones de una función aleatoria, de la cual la realidad se considera
también como una realización particular.
• El éxito de la simulación geoestadística está dado en que permite mostrar una
representación de la realidad, nunca conocida totalmente en la mayor cantidad de
problemas prácticos de la ingeniería de petróleos. La simulación geoestadística no
proporciona el mejor valor a obtener en una localización, pero si muestra lo que podría
ser el fenómeno real en una forma general, pues reproduce las dos principales
características de los datos experimentales que son el histograma y el semivariograma.
• No se debe distinguir entre simulación geoestadística o estimación como dos
procedimientos que se contradicen, aunque sus objetivos son diferentes, pero
dependiendo de los intereses particulares en una aplicación o investigación se deben
usar com dos alternativas que complementan un estudio: la simulación geoestadística
para obtener una representación del fenómeno en forma general y la estimación para
obtener el mejor valor estimado en cualquier localización puntual o de bloque.
• Respecto a la evaluación de proyectos, la toma de decisiones requiere de la definición
de objetivos y un detallado entendimiento del proyecto bajo estudio. Tanto las opciones
disponibles con las que cuenta quien toma la decisión como las formas en las que el
proyecto puede ser mejorado, deben ser totalmente entendidas. Debe definirse un
método para medir la utilidad del proyecto. Tal utilidad es una función de indicadores
críticos cuyos valores dependen de parámetros de control. Una gran variedad de
técnicas de simulación pueden ser usadas para obtener las distribuciones de
probabilidad de la utilidad del proyecto y el riesgo en que incurre puede ser reducido
mediante mejoras de la información recolectada y estimada. La simulación
geoestadística juega un papel crítico en la evaluación del riesgo, sin embargo, las
limitaciones de la simulación geoestadística no deben ser olvidadas: el método es
vulnerable a la calidad de la información utilizada, y algo mucho más importante aún es
que si la distribución de la variable que se pretende estudiar no es normal el
77procedimiento de simulación Gaussiana no podría ser llevado a cabo.
• Los resultados obtenidos de continuidad espacial de las propiedades porosidad y
permeabilidad del ejemplo de aplicación revelan que ambas propiedades presentan
isotropía geoestadística en el área de estudio. Este hecho implicó que solo se necesitase
de un variograma para modelar el comportamiento de las variables, pero esto no quiere
decir que este fenómeno siempre se vaya a presentar. Los mapas de variografía son la
herramienta indicada para determinar si existe o no la isotropía, pero no puede olvidarse
que estos se deben combinar con el buen juicio y criterio de quien está realizando el
estudio para tomar la decisión de la existencia de dicho fenómeno.
• Uno de los problemas que se tenían después de realizar la simulación geoestadística es
que se obtuvieron múltiples realizaciones de una misma propiedad y todos ellas con la
misma probabilidad de ocurrencia. Siempre se dijo que era necesario encontrar una
utilidad o aplicación a las realizaciones obtenidas con la simulación geoestadística. En
este trabajo se demostró que es posible encontrar distribuciones de probabilidad de
propiedades de interés a la hora de evaluar un proyecto petrolífero. En este trabajo se
encontró una distribución de frecuencias para la variable petróleo final recuperado a
quince años. Dicha distribución mostró el comportamiento de una variable normal.
Este tipo de distribuciones permiten tener la posibilidad de estudiar varios escenarios de
producción y de esta forma analizar cual sería el más conveniente para poder tomar la
decisión más acertada económicamente hablando.
788. RECOMENDACIONES
• La principal recomendación derivada de este trabajo consiste en la utilización de más de
una variable en los procesos de estimación y simulación geoestadística. Es decir,
siempre que sea posible se debe realizar un estudio de corregionalización para
garantizar que el modelo final generado reproduce la mayor cantidad de características
experimentalmente medidas como lo son atributos sísmicos, parámetros petrofísicos y
parámetros de producción. Dependiendo del tipo de estudio que se vaya a desarrollar se
entraría a decidir que es lo importante, si obtener un único modelo de yacimiento o si
obtener distribuciones de frecuencia como las obtenidas en este trabajo. Si el caso es el
de resolver el problema de unicidad del modelo, entonces la corregionalización podría
ser de gran utilidad ya que mediante ella se puede probar que modelo de porosidad o
permeabilidad cumple no solo con las condiciones de histograma y correlación espacial,
sino también con la correlación que existe entre ella y otro tipo de variables propias de
la formación.
• Seria conveniente trabajar el ejemplo de aplicación con áreas más pequeñas ya que la
selección de áreas de gran tamaño ocasionan que la simulación geoestadística
enmascare variaciones a pequeña escala. La selección del tamaño de dichas áreas debe
ser tal que permita buscar un variograma para encontrar un modelo de continuidad
espacial que refleje los fenómenos a pequeña escala, pero debe considerarse que dicha
área debe contener datos suficientes que faciliten el cálculo del variograma. Lo ideal es
encontrar un balance entre las dimensiones del área de estudio y el número de datos
contenidos en ella.
• A pesar de que en este trabajo se explicó el porqué de la selección del variograma
esférico para porosidad y permeabilidad, se recomienda modelar estos datos con un
modelo combinado o uno que ajuste los efectos de oscilación que se presentaron en los
variogramas experimentales. El programa utilizado en el ajuste de variogramas fue la
limitante que no permitió modelar el efecto de oscilación presentado por el variograma,
es recomendable programar modelos de variograma sinusoidales y cosenoidales para
aplicarlos al grupo de datos utilizados en este trabajo.
79BIBLIOGRAFIA
• R. MOHAN, Srivastava y EDWARD H., Isaaks. “Applied Geostatistics” Ed. OXFORD
UNIVERSITY PRESS, 1989.
• HOHN, Michael Edward. “Geostatistics and Petroleum Geology” Ed. KLUWER
ACADEMIC PUBLISHERS, 1999. Segunda Edición.
• Kashib, Tarum. “Iterative Integration of Dyamic Data in Reservoir Models” SPE 2003. pp
1 -12
• DEUTSCH, C. V. y JOURNEL, A. G. “GSLIB: Geostatistical Software Library and User’s
Guide”. Ed. OXFORD UNIVERSITY PRESS, New York 1986.
• L. Wang y P. M. Wong. “Modeling Porosity Distribution in the A’nan Oil Field: Use of
Geological Quantification, Neural Networks and Geostatistics” SPE Reservoir Evaluation &
Engineering vol. No. 6, Diciembre 2 de 1999. pp 1 - 6
• Hisham M., Al Qassab. “Cross-Discipline Integration in Reservoir Modeling: The Impact
on Fluid Flow Simulation and Reservoir Managment” SPE 62902, 2000. pp 1 - 12
• Lang, Zhang y Yanis C. Yortsos. “A direct Method for the Identification of the Permeability
Field of an Anisotropic Porous Medium” SPE 62976, 2000.
• Queipo, V. Nestor y Pintos, Salvador. “Surrogate Modeling – Based Optimization for the
Integration of Static and Dynamic Data Into a Reservoir Description” SPE 63065, 2000.
• Landa, J.L. y Horne, R.N. “Reservoir Characterization Constrained to Well-Test Data: A
Field Example” SPE Reservoir Evaluation & Engineering vol. 3 No. 4, Agosto de 2000. pp
1 - 10
• Yoon, Seongesik. “A Multiscale Approach to Production Data Integration Using Stream-
Line Models” SPE Journal, Junio 2001.
• Cayeux, Eric. “Well Plannig Quality Improved Using Cooperation Between Drilling and
Geosciences” SPE 71331, 2001. pp 1 – 9
• Zabalza-Mezghani, I. “Constraining Reservoir Facies Models to Dynamic Data – Impact of
Spatial Distribution Uncertainty on Production Forecast” SPE 71335, 2001. pp 1 – 15
• Landa, Jorge. “Technique to Integrate Production and Static Data in a Self-Consisted Way”.
SPE 71597, 2001. pp 1 – 11”• Pelgrain de Lestang, A. “Geologically Oriented Geostatistics: An Integrated Tool for
Reservoir Studies” SPE 74371, 2002. pp 1 – 10
• Liu, Xiaohuan. “Geological Characterization of Naturally Fractured Reservoir Using
Multiple Point Geostatistics” SPE 75246, 2002. pp 1 – 12
• Arpat, Burc. “Featured Based Geostatistics: An Application To A Submarine Channel
Reservoir” SPE 77426, 2002. pp 1 – 9
• Caers, Jef. “Geostatistical History Matching Under Training-Image Based Geological
Model Constraints” SPE 77429, 2002. pp 1 – 6.
• Spat, Jeff. “Waterflood Optimization Using A combined Geostatistical – 3D Streamline
Simulation Approach: A Field Example” SPE 38355, 1997. pp 1 - 12••••••••••••••••
REFERENCIAS
1 A. Quenes y S. Richardson. “Fractured Reservoir Characterization and Performance
Forecasting Using Geomechanics and Artificial Intelligence” SPE 30572, 1995. pp 429 –
431
2 E. Gringarten. “Methodology for Variogram Interpretation and Modeling for Improved
Reservoir Characterization” SPE 56654, 1999. pp 1 – 13
3 R. MOHAN, Srivastava y EDWARD H., Isaaks. “Applied Geostatistics” pp 370 - 397
4 Journel, A.G. y Froidevaux, R.”Anisotropic Hole-Effect Modeling” Matemáticas para
Geología vol. 14 pp 217 - 247, 1982.
5 Cressie, N. A. C. “Statistics for Spatial Data” Ed. John Wiley & Sons, 1980. pp 900
6 HOHN, Michael Edward. “Geostatistics and Petroleum Geology” Ed. KLUWER
ACADEMIC PUBLISHERS, 1999. Segunda Edición. Anexo A
7 Almeida, A. S. y Journel, A.G. “Joint Simulation of Multiple Variables with a Markov-Type
corregionalization Model”. Matemáticas para Geología vol. 5. pp 565 - 588
8 HOHN, Michael Edward. “Geostatistics and Petroleum Geology” Ed. KLUWER
ACADEMIC PUBLISHERS, 1999. Segunda Edición. pp 190 - 203
9 Xu, W., Tran, T. y Srivastava R. M. “Integrating Seismic Data in Reservoir Modeling: The
Collocated Cokriging Alternative” Proceeding of the Society of Petroleum Engineers
Technical Conference, 1992. pp 833 - 842
10 HOHN, Michael Edward. “Geostatistics and Petroleum Geology” Ed. KLUWER
ACADEMIC PUBLISHERS, 1999. Segunda Edición. pp 180 - 203
11 Isaaks, E. “The Application of Monte Carlo Methods to the Analysis of Spatially Correlated
Data”. Disertación de PhD. Standford University, 1990.
12 R.M. Srivastava. “An Application of Geostatistical Methods for Risk Analisys in Reservoir
Management” SPE 20608, 1990. pp 1 - 10
13 G. James.“Prediction of Oil Production With Confidence Intervals” SPE 66350, 2001.pp 1- 614 S. Yoon y A. Datta-Gupta. “A multiscale approach to Production Data Integration Using
Stream-Lines Models” SPE 56653, 1999. pp 1 -12
15 R. Pathatk y D. O. Ogbe. “Application of Geostatistical and Flow Simulation to Evaluate
Options for Well Placement” SPE 62554, 2000. pp 1 - 11
16 www.statios.comANEXO A
NÚMERO DE PARES POR PUNTO CALCULADO DEL SEMIVARIOGRAMA DE
LA FIGURA 16. (Datos tomados de Hohn, Michael Edward)
Distancia (Km) No. Pares
0.18 64
0.56 1791
1.02 3456
1.51 4863
2.01 6220
2.51 7448
3.01 8685
3.5 9733
4.01 10730
4.5 11510
5 12329
5.5 13044
6 13869
6.5 14481
7 15147
7.5 15523
8 16160
8.5 16928
9 17288
9.5 17648
10 17957ANEXO B
INFORMACIÓN PETROFÍSICA DE LOS 241 POZOS DEL CAMPO ESTUDIADO
Espesor de Tope de Porosidad K
Pozo E (ft) N (ft)
Formacion (ft) Formacion (ft) (%) (mD)
1 11564 5691 80 4275 7.32 3.51
2 10679 13706 93 4232 8.79 5.81
3 6311 36307 89 4317 7.12 0.77
4 9386 24799 105 4262 9.57 6.29
5 10761 24755 95 4278 7.74 2.9
6 10769 23463 95 4264 10.89 12.32
7 9145 23413 107 4252 6.9 0.79
8 8201 23557 93 4258 9.86 8.19
9 6727 23472 93 4200 7.22 1.45
10 3834 23446 101 4120 7.55 1.6
11 2780 23494 94 4104 7.6 1.57
12 3994 22098 101 4112 7.6 1.51
13 5482 22214 91 4153 7.35 1.79
14 6692 22138 86 4191 9.55 8.21
15 9324 22199 98 4230 9.74 4.48
16 10817 22305 101 4253 8.62 2.73
17 10761 20774 87 4264 8.63 2.68
18 9350 20761 94 4225 9.27 3.42
19 8237 20748 93 4207 7.74 3.82
20 5354 20804 87 4138 6.66 2.43
21 4013 20873 98 4110 7.21 2.72
22 4214 19455 92 4110 7.65 2.7
23 5372 19467 93 4136 9.3 9.63
24 10748 19443 106 4251 9.52 6.53
25 10728 18125 102 4246 7.35 2.58
26 9348 18125 96 4225 9.77 7.27
27 8298 18147 91 4211 11.84 25.78
28 5368 18157 93 4132 9.28 7.44
29 4101 18042 97 4126 9.18 4.56
30 3910 16690 96 4128 7.45 2.59
31 5341 16788 99 4128 9.03 8.34
32 6685 16852 90 4154 9.85 13.15
33 8300 16803 91 4197 11.82 19.29
34 9385 16770 90 4221 9.3 6.37
35 10776 16781 103 4255 8.49 2.78Espesor de Tope de Porosidad K
Pozo E (ft) N (ft)
Formacion (ft) Formacion (ft) (%) (mD)
36 10740 15545 97 4245 9.82 11
37 9241 15332 95 4204 11.77 15.94
38 8225 15385 95 4190 11.25 15.15
39 6637 15384 90 4161 10.59 13.36
40 5370 15489 95 4140 10.25 13.55
41 4110 15328 94 4129 6.02 0.96
42 4222 14324 93 4142 8.48 2.44
43 5503 14290 98 4158 9.18 8.58
44 6683 14094 97 4176 8.69 5.8
45 7640 18817 90 4197 10.95 20.52
46 12775 24773 90 4315 8.08 2.39
47 12260 23473 88 4286 9.03 3.34
48 12185 22070 88 4275 9.18 4.45
49 12180 19494 90 4274 8.98 4.72
50 12148 18259 101 4268 7.93 3.69
51 12144 16998 95 4299 9.47 5.72
52 12128 15660 99 4317 11.63 14.36
53 11277 28705 93 4361 10.01 4.62
54 12283 28712 93 4396 8.88 1.44
55 13036 36272 89 4482 6.43 0.31
56 10323 32598 93 4337 7.64 0.88
57 17500 31524 81 4610 8.07 1
58 12185 29848 86 4393 7.2 0.91
59 9497 28845 94 4307 9.45 2.53
60 9443 27493 97 4309 11.33 10.4
61 12051 26104 93 4353 10.06 3.58
62 13436 26010 94 4410 10.25 8.53
63 16926 29997 83 4587 8.39 2.38
64 14643 25975 88 4446 10.46 4.67
65 15963 26086 86 4475 9.41 3.73
66 17138 25993 84 4502 8.97 2.3
67 9265 12786 97 4204 9.63 12.36
68 13046 11699 81 4310 12.45 19.94
69 10670 10061 86 4237 7.43 3.56
70 10557 37687 92 4443 6.95 0.61
71 11555 37376 94 4457 7.9 0.79
72 1325 36860 94 4253 9.56 2.73
73 2694 36849 104 4297 9.93 2.78
74 4032 36875 100 4315 8 1.49
75 10213 36236 90 4379 8.1 0.84
76 11622 36323 93 4431 8.39 1.37
77 1375 35515 90 4216 8.88 4.82
78 2618 35342 96 4205 9.59 3.64
79 3915 35341 90 4258 8.63 2
80 5444 35644 91 4312 8.32 1.42Espesor de Tope de Porosidad K
Pozo E (ft) N (ft)
Formacion (ft) Formacion (ft) (%) (mD)
81 9804 35321 90 4352 8.06 0.97
82 11150 35267 93 4404 9.93 1.59
83 1393 34225 86 4186 8.08 3.47
84 2628 34162 87 4182 6.24 0.84
85 4064 34105 98 4225 8.59 2.15
86 5478 34237 93 4272 7.41 1.78
87 6905 33995 103 4280 5.48 0.72
88 8260 34066 90 4298 7.27 0.46
89 9172 34206 87 4310 6.92 0.34
90 10420 33813 103 4352 10.88 2.66
91 11693 33847 102 4407 10.81 2.89
92 1356 32818 85 4179 8.47 2.59
93 2647 32772 88 4198 6.91 2.19
94 4047 32789 94 4212 6.03 0.58
95 5391 32728 95 4240 5.58 0.36
96 6765 32748 94 4235 8.05 2.31
97 8295 32709 89 4256 7.78 0.67
98 8890 32531 93 4275 8.66 0.91
99 11587 32674 92 4380 7.81 0.77
100 12931 33125 89 4441 8.7 1.46
101 1424 31212 91 4167 8.57 4.37
102 2641 31413 90 4230 8.04 3.66
103 3893 31373 87 4219 7.28 3.09
104 5387 31450 92 4230 7.65 2.83
105 6885 31391 97 4228 7.68 1.79
106 7933 31309 92 4256 8.72 1.62
107 9155 31127 97 4282 7.55 0.83
108 11651 31433 90 4373 8.29 3.59
109 1439 30286 81 4178 7.42 3.03
110 3986 30123 90 4255 7.75 2.29
111 5301 30022 90 4258 7.24 1.96
112 6856 30278 84 4242 6.81 0.6
113 8154 30096 96 4262 6.05 0.58
114 9197 30071 97 4285 9.29 3.58
115 10546 29983 95 4326 10.21 8.25
116 1393 28730 94 4169 8.74 4.19
117 2650 28895 83 4225 8.27 4.14
118 3717 28860 95 4243 7.72 2.34
119 5362 28775 98 4242 6.71 1.34
120 6898 28887 87 4257 7.64 1.13
121 8009 28756 81 4290 8.05 1.45
122 14187 29319 86 4475 10.54 5.25
123 15482 29316 78 4516 8.17 1.88
124 16214 28518 82 4532 8.21 2.53
125 17015 28642 86 4568 7.6 1.92Espesor de Tope de Porosidad K
Pozo E (ft) N (ft)
Formacion (ft) Formacion (ft) (%) (mD)
126 1198 27160 89 4145 8.12 2.73
127 2645 27530 88 4178 8.61 3.4
128 3927 27420 98 4200 7.25 1.23
129 5308 27552 88 4234 6.5 1.58
130 6751 27428 94 4268 6.74 0.73
131 7989 27440 91 4294 9.32 4.9
132 10582 27315 91 4335 9.52 1.99
133 13544 27547 99 4425 9.32 5.42
134 14655 27507 88 4477 8.59 1.37
135 1555 26108 88 4131 8.59 2.98
136 2634 26118 93 4145 7.32 1.79
137 3698 26033 100 4151 7.04 1.23
138 5396 26135 93 4197 7.4 1.85
139 6759 26116 86 4259 8.07 1.49
140 7963 26117 88 4290 9.28 5.31
141 9400 26109 97 4290 10.94 11.08
142 10870 26129 99 4320 9.02 1.9
143 1702 24835 94 4110 8.91 3.45
144 2595 24846 98 4124 8.06 2.92
145 3993 24811 96 4141 7.69 1.81
146 5369 24829 91 4189 7.97 1.38
147 6730 24792 90 4220 8.44 3.25
148 7956 24795 98 4260 7.57 0.84
149 1936 23645 99 4098 8.99 2.75
150 1708 22288 104 4086 8.44 3.33
151 2593 21886 100 4098 8.08 2.13
152 1871 20786 118 4107 8.11 7.57
153 3105 19578 93 4112 8.43 3.39
154 3086 18238 98 4136 6.44 0.9
155 3061 16831 105 4138 7.74 2.82
156 10701 12849 94 4221 10.45 10.76
157 12792 12843 95 4292 11.84 18.07
158 9346 11530 100 4209 10.27 7.6
159 10705 11397 95 4223 10.16 5.32
160 11623 11384 93 4254 12.15 16.97
161 9281 10037 93 4220 9.28 4.03
162 13534 10053 74 4337 10.15 10.08
163 9337 14054 87 4222 12.11 15.03
164 11913 14054 100 4278 12.02 18.48
165 5362 12801 95 4136 8.06 5.62
166 6695 12811 97 4166 10.49 13.13
167 7999 12841 106 4184 6.14 0.58
168 3905 11547 103 4169 6.73 1.34
169 5357 11550 95 4144 8.67 4.13
170 6678 11442 100 4164 10.42 10.54Espesor de Tope de Porosidad K
Pozo E (ft) N (ft)
Formacion (ft) Formacion (ft) (%) (mD)
171 4288 9779 102 4166 7.58 1.05
172 5395 10112 90 4166 9.1 4.06
173 7874 10134 93 4190 7.25 1.03
174 1269 8630 89 4169 6.38 0.73
175 2826 8844 105 4153 5.61 0.51
176 4128 8977 87 4170 7.04 0.88
177 5360 8829 84 4180 6.86 0.85
178 6569 8935 84 4194 9.51 3.13
179 7907 8624 93 4214 6.83 1.63
180 9353 8720 92 4231 8.08 1.89
181 10998 8611 87 4248 8.53 2.52
182 12455 8916 77 4280 6.97 1.24
183 14080 9042 83 4362 6.26 0.88
184 17431 8946 87 4454 7.12 1.12
185 2960 7413 96 4177 5.26 0.47
186 5316 7480 99 4178 7.26 0.91
187 6604 7308 83 4192 6.87 0.97
188 7905 7493 95 4216 6.68 2
189 12809 7543 76 4271 6.74 1.15
190 14028 7628 92 4309 8.32 1.46
191 14971 7466 79 4336 7.7 2.05
192 16328 7362 75 4372 7.33 1.25
193 2844 6197 105 4206 4.81 0.31
194 5323 6081 92 4192 6.3 0.56
195 6718 5962 87 4204 7.1 1.71
196 8238 6067 88 4210 9.1 2.06
197 9093 6225 89 4230 8.29 1.9
198 9975 6612 82 4234 7.07 1.43
199 14754 5984 89 4324 9.18 2.08
200 16150 6138 97 4386 7.08 0.94
201 4018 4797 93 4248 7.42 1
202 5288 4742 91 4230 8.83 1.43
203 6565 4968 90 4209 8.8 2.05
204 7949 4764 90 4204 7.93 1.85
205 9196 4716 89 4250 9.14 2.28
206 10497 4742 87 4262 7.92 1.94
207 12000 4781 80 4291 9.03 2.95
208 13533 4639 79 4302 10.2 5.48
209 16146 4701 75 4340 7.42 1.08
210 5339 3353 99 4252 9.05 5.07
211 6688 3341 93 4240 9.48 4.05
212 7807 3550 91 4231 8.87 1.84
213 8869 3544 89 4238 6.5 0.64
214 10477 3612 84 4252 9.28 2.55
215 12048 3536 80 4260 8.35 2.89Espesor de Tope de Porosidad K
Pozo E (ft) N (ft)
Formacion (ft) Formacion (ft) (%) (mD)
216 13274 3536 82 4292 8 1.43
217 14710 3541 80 4299 7.3 1.11
218 15982 3522 74 4312 7.52 1.5
219 18437 3665 85 4452 6.82 0.73
220 6752 1913 81 4269 8.31 1.24
221 9096 2089 91 4246 7.18 0.79
222 12144 2106 82 4245 8.95 4.01
223 13536 2096 75 4282 9.4 3.85
224 16014 2014 85 4300 8.98 1.86
225 9310 769 103 4237 7.04 2.16
226 13539 807 86 4267 7.31 1.18
227 13329 17265 86 4318 8.58 4.35
228 13329 16061 87 4357 9.43 6.8
229 13209 14475 84 4428 9.68 5.1
230 1894 7377 98 4178 5.47 0.47
231 6800 35201 102 4312 6.02 0.62
232 1670 6207 103 4210 7.01 0.95
233 1545 3624 101 4264 4.9 0.33
234 2828 1906 109 4270 6.74 0.68
235 9142 36445 91 4350 8.15 0.79
236 12404 35086 87 4454 6.77 0.28
237 8736 37564 96 4356 6.39 0.48
238 9579 37746 99 4390 7.64 0.91
239 18799 25955 89 4587 9.92 4.5
240 4899 993 101 4286 6.13 0.53
241 3821 769 93 4266 6.4 0.7ANEXO C
RESULTADOS DE LA SIMULACIÓN DE FLUJO UTILIZADOS EN LA
CONSTRUCCIÓN DEL HISTOGRAMA DE FRECUENCIAS ACUMULADAS DE
PETRÓLEO RECUPERADO
Petróleo
Frecuencia Acumulada
Recuperado
(%)
(STB)
47,093,190 1.02
47,159,740 2.041
47,331,800 3.061
47,417,980 4.082
47,575,960 5.102
47,807,800 6.122
48,165,840 7.143
48,272,390 8.163
48,293,720 9.184
48,326,380 10.204
48,412,830 11.224
48,498,390 12.245
48,560,090 13.265
48,649,190 14.286
48,730,190 15.306
48,759,620 16.327
48,765,200 17.347
48,828,700 18.367
48,910,260 19.388
48,982,510 20.408
49,032,200 21.429
49,144,860 22.449
49,153,250 23.469
49,208,960 24.49
49,215,400 25.51
49,306,800 26.531
49,312,630 27.551
49,331,560 28.571
49,439,680 29.592
49,441,770 30.612
49,481,290 31.633
49,497,780 32.653Petróleo
Frecuencia Acumulada
Recuperado
(%)
(STB)
49,617,510 33.673
49,649,000 34.694
49,669,330 35.714
49,754,180 36.735
49,780,410 37.755
50,012,430 38.776
50,107,760 39.796
50,212,310 40.816
50,230,370 41.837
50,294,420 42.857
50,303,500 43.878
50,303,500 44.898
50,314,090 45.918
50,341,920 46.939
50,344,400 47.959
50,407,120 48.98
50,506,980 50
50,540,090 51.02
50,603,210 52.041
50,633,960 53.061
50,653,490 54.082
50,673,120 55.102
50,700,890 56.122
50,770,890 57.143
50,776,810 58.163
50,892,830 59.184
51,019,270 60.204
51,028,620 61.224
51,107,200 62.245
51,177,370 63.265
51,280,020 64.286
51,330,990 65.306
51,341,090 66.327
51,451,070 67.347
51,664,220 68.367
51,723,780 69.388
51,753,020 70.408
51,794,350 71.429
51,826,990 72.449
51,827,440 73.469
51,846,360 74.49
51,875,150 75.51
51,875,760 76.531
51,930,980 77.551Petróleo
Frecuencia Acumulada
Recuperado
(%)
(STB)
51,987,430 78.571
51,997,480 79.592
52,080,910 80.612
52,124,120 81.633
52,298,620 82.653
52,390,160 83.673
52,430,720 84.694
52,488,750 85.714
52,505,150 86.735
52,548,980 87.755
52,612,470 88.776
52,720,060 89.796
52,816,000 90.816
52,889,470 91.837
53,035,390 92.857
53,084,930 93.878
53,198,400 94.898
53,226,250 95.918
53,311,020 96.939
53,444,270 97.959
53,938,960 98.98
54,376,130 100