04 Simulacion Multivariable
04 Simulacion Multivariable
04 Simulacion Multivariable
Resumen
• Simulación multivariable
• Introducción simulación
• Repaso simulación univariable
• Cosimulación secuencial gaussiana
• Aplicaciones
Motivación
• El kriging suaviza: los valores estimados son menos dispersos que los valores verdaderos
• no se puede predecir la ocurrencia de valores extremos
• puede dar la impresión de que todo es mineral (leyes estimadas sobre la ley de corte,
cuando esta última es baja)
• Simulación condicional: reproduce además los valores de los datos en los sitios de
muestreo. Produce realizaciones de la función aleatoria condicionada a los valores de los
datos (noción de distribución condicional o a posteriori).
más realista
en un sitio con dato, todos los escenarios reproducen este dato
lejos de los datos, la simulación se vuelve no condicional
Motivación
No condicional Condicional
Modelo plurigaussiano
Algoritmos de Simulación
• Método matricial (Descomposición LU): poco usada debido a las restricciones de tamaño
(debe resolverse una matriz de N x N donde N podría corresponder a millones). LUSIM
• Bandas Rotantes: simula la variable en líneas (1-D) y combina los valores simulados (3-D).
Requiere condicionar a datos locales.
• Fractales: simulación no condicional rápida que requiere procesar a posteriori para
respetar los datos locales.
• Espectral: Los métodos que usan FFTs pueden ser rápidos… pero no son intuitivos y
tienen limitaciones con respecto al tamaño de la grilla y datos condicionantes.
• Annealing: Basado en convergencia de una función objetivo a cero. Extremadamente
flexible (y potencialmente muy lento), no tiene una base teórica. SASIM
• Simulación Secuencial: ampliamente usada y recomendada para simulación condicional.
SGSIM, SISIM
Principios de simulación
• La idea de la simulación es agregar al mapa de kriging la variabilidad perdida
• Para ello, se agrega al valor estimado por kriging simple, un residuo de media cero y
varianza igual a la varianza de kriging simple
• Esto permite que la relación entre un dato y un valor simulado sea correcta se
reproduce la variabilidad espacial
Principios de simulación
• La covarianza no cambia:
Caso 1 Caso 2
𝜎2 𝜎2
𝜇 = 𝑧∗ 𝑢 𝜇 𝜇
2
𝜎 2 = 𝜎𝑆𝐾
Modelo Gaussiano para simular leyes
Distancia
Modelo Gaussiano para simular leyes
Distancia
Modelo Gaussiano para simular leyes
• Noción de fluctuación estadística: variogramas de los valores simulados vs. variograma
teórico (100 realizaciones)
• Cada realización constituye un escenario posible, por lo que se puede trabajar con ella como
si fuera la realidad. En consecuencia, cada realización entrega una respuesta insesgada al
problema planteado.
• La respuesta del modelo de kriging podría estar sesgada, debido a que el kriging no
reproduce la variabilidad de las leyes.
Aplicaciones
Ejemplo 1: variabilidad en la ley extraída según una secuencia de extracción
predeterminada
0,9000 Sim_01
Sim_02
Sim_03
Sim_04
0,8000 Sim_05
Sim_06
Sim_07
Sim_08
Sim_09
0,7000 Sim_10
Sim_11
Sim_12
Grade
Sim_13
0,6000 Sim_14
Sim_15
Sim_16
Sim_17
Sim_18
0,5000 Sim_19
Sim_20
Sim_21
Sim_22
0,4000 Sim_23
Sim_24
Sim_25
Avg_Sim
Kriging
0,3000
2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016
Year
Título
Cantidad de columnas
Coordenada x
Coordenada y
Coordenada z
Variable 1
Variable 2
Geología
01-dofile
• Script Python: separación de los datos en diferentes archivos según categoría
geológica.
02-nscore
• Gslib: transformación de variables continuas a distribución gaussiana. (mostrar
.par y una validación de scatplt)
• Para ejecución, arrastrar archivo .par al ejecutable.
03-mergefiles
• Script Python: Unión de los datos transformados (nscore) en un solo archivo con
formato para rutinas posteriores.
04-gamv
• Gslib: rutina para calcular variogramas experimentales (directos y
cruzados para multivariable).
05-varfit
• Gslib: generación de modelos semiautomáticos para variogramas
experimentales
06-gamvindicadores
• Gslib: Creación de variograma de indicadores (códigos1-0) para la utilización de
Blocksis.
• Cambio en último parámetro de Gamv:
07-varfitindicadores
• Gslib: Generación de modelo de variograma de indicadores para Blocksis
08-blocksis
• Gslib: Rutina para generar simulaciones de variables categóricas. Permite realizar
suavizamiento y simulación en función de proporciones.
08-blocksis
• Gslib: Rutina para generar simulaciones de variables categóricas. Permite realizar
suavizamiento y simulación en función de proporciones.
09-usgsim
• Gslib: Para cada categoría (en la geología) se simulan las leyes (variable continua).
09-usgsim
• Captura de relaciones multivariables
𝜌 𝑐𝑢, 𝑎𝑢 = +0.7