Orientador: Carlos Alberto C. Altemani / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-31T21:58:36Z (GMT). No. of bitstreams: 1 Bazani_MarcioAntonio_D.pdf: 8200093 bytes, checksum: 0eca0ed6afa5487d222dd86214f3f802 (MD5) Previous issue date: 2001 / Resumo: Este trabalho considerou o resfriamento convectivo de placas planas localizadas na superfície de incidência de um jato confinado de ar. Inicialmente o problema foi simulado numericamente, considerando a superfície de incidência do jato aquecida e isotérmica, utilizando dois modelos de turbulência de alto Reynolds distintos: k-e RNG (grupos de renormalização). As simulações foram efetuadas para um escoamento bidimensional incompressível de ar. O método dos volumes de controle foi usado para resolver iterativamente as equações de conservação de massa, quantidade de movimento e energia, bem como as equações de energia cinética turbulenta e taxa de dissipação de energia cinética turbulenta. Os resultados obtidos da simulação de escoamento e da transferência de calor foram comparados com resultados numéricos e experimentais obtidos na literatura numa faixa do número de Reynolds na seção de entrada (8000 < Rej < 20000) e da razão de aspecto (1,0 < H/w < 8,0) do jato. Em seguida, uma montagem experimental foi construída, com três placas planas aquecidas montadas na superfície do jato. Os resultados de testes experimentais de laboratório foram comparados com resultados correspondentes de simulação numérica. Sob condições distintas de aquecimento das três placas, o coeficiente adiabático de transferência de calor, descrito por Moffat (1998), foi utilizado e comparado com resultados do coeficiente convectivo baseado na temperatura de entrada do jato. Nos testes experimentais, o número de Reynolds do jato variou na faixa (9000 < Re < 16000) e a razão de aspecto foi mantida fixa (H/w = 6,0). Foi verificado, pelos resultados experimentais e de simulação, que o coeficiente adiabático (had) depende apenas das condições do escoamento, enquanto que o coeficiente baseado na temperatura de entrada do ar (hin) depende tanto do escoamento quanto da potência dissipada nas placas a montante. Os resultados deste estudo tem aplicações no resfriamento de componentes eletrônicos / Abstract: The convective cooling of discrete heated plates flush mounted on the incidence surface of a confined slot jet was investigated in this Thesis. Initially numerical simulations were performed considering a heated isothermal impingement plate, using two high-Reynolds turbulence models: k-e and RNG (Groups of Renormalization). A control volumes method was used to solve a two-dimensional incompressible airflow. The conservation equations (mass, momentum and energy), and those associated to the turbulence mo del (the turbulent kinetic energy equation and that of its dissipation rate) were solved iteratively by the contral volumes method. The numerical simulation results for the flow and heat transfer were compared with numerical and experimental data obtained from the literature to explore the effects of the jet Reynods number (8000 < Re < 20000), and the aspect ratio (1 < H/w < 8). An experimental apparatus was built, with three heated plates flush mounted on the impingement plate. The obtained experimental results were compared with corresponding numerical simulations. Under distinct heating conditions of the three plates, the adiabatic heat transfer coefficient, described by Moffat (1998), was compared with the convective coefficient based on the inlet air temperature. The experimental tests were performed for a range of 9000 < Rej < 16000 and for a single aspect ratio, Hlw = 6. It was verified, by the experimental results and the numerical simulations, that the adiabatic heat transfer coefficient (had) depends on the flow conditions, while that based on the inlet air temperature (hin) depends both on the flow conditions and power dissipated in the upstream plates. The results of this investigation are relevant to applications in electronics cooling / Doutorado / Termica e Fluidos / Doutor em Engenharia Mecânica

Magíster en Minería / Este documento aborda el tema de la validación de las simulaciones geoestadísticas empleadas para generar modelos numéricos en geociencias. Estos modelos son de gran importancia en la minería, toda vez que muchas de las decisiones que en esta industria se toman, parten del conocimiento generado por estas interpretaciones de la realidad. En ese orden de ideas, el documento examina el estado del arte en cuanto a diversas métricas empleadas para validar las simulaciones geoestadísticas, haciendo énfasis en las simulaciones multi-punto. En primer lugar se hace un recuento de diversas propuestas efectuadas por algunos autores quienes han presentado sus propuestas o enfoques relacionados con criterios que se deberían tener en cuenta al momento de validar este tipo de modelos. Esto se clasificó en 4 niveles de validación: un primer nivel de condicionamiento de los datos, un segundo nivel de validación de estadística, un tercer nivel de validación geoestadística y finalmente un nivel de validación de predicción. Se propone acorde a esto una metodología para abordar el tema, consistente en la aplicación a dos casos prácticos de diferentes métricas de validación para determinar su potencial uso en diferentes situaciones. Un primer caso está relacionado con una variable continua, la cual fue simulada a partir de información extraída de un modelo sintético realístico de plegamientos. Las simulaciones se efectuaron con dos tipos de algoritmos. El primero corresponde a un algoritmo paramétrico que emplea variogramas para su inferencia como es el caso de SGSIM (simulación secuencial gaussiana). El segundo es un algoritmo no paramétrico (simulación multi-punto) que usa una imagen de entrenamiento para inferir la función de distribución y simular las diferentes realizaciones (FILTERSIM simulación basada en puntuación por filtros). El segundo caso es una variable categórica de un modelo de formaciones aluviales, del cual fue extraída la información. En este caso se emplearon dos algoritmos no paramétricos. Uno de ellos es SISIM (simulación secuencial de indicadores) que emplea los variogramas de una variable indicador para inferir el modelo. El otro corresponde a SNESIM (simulación ecuación normal simple) que emplea una imagen de entrenamiento para determinar la función de distribución de probabilidad para poder generar el modelo. A partir de estas simulaciones se ponen en práctica diversas métricas para validar los modelos obtenidos, lo cual permite determinar el desempeño de estas en el proceso de validación. Con dicha información y a partir de revisar los supuestos fundamentales relacionados con la Estacionariedad de las variables y de la inferencia de los modelos, se hace una propuesta de pasos a seguir para validar un modelo a partir de simulaciones geoestadísticas.

Magíster en Ciencias, Mención Computación / El modelamiento geológico es uno de los primeros pasos en la estimación de recursos en un yacimiento. Este modelamiento se realiza a partir de una interpolación basada en un número limitado de muestras obtenidas en sondajes y tiene por objetivo reproducir de la manera más exacta lo que ocurre en el subsuelo. Esta representación es de suma importancia debido a que influye directamente sobre la toma de decisiones y, por lo tanto, en la metodología para explotar el yacimiento, la rentabilidad y otros asuntos operacionales. El objetivo de este trabajo es implementar una nueva metodología para el modelamiento de cuerpos geológicos utilizando geoestadística para estimar y simular la distancia de las muestras al contacto entre el cuerpo mineralizado y la roca caja que lo contiene. Esta nueva metodología permite incorporar información interpretativa o blanda, en los algoritmos de estimación mediante kriging ordinario y simulación secuencial gaussiana. El trabajo consistió en modificaciones a GSLIB, ésta es una librería que contiene los algoritmos más usados en geoestadística. Las rutinas modificadas fueron dos: \textit{kt3d}, que corresponde al algoritmo de kriging, y \textit{sgsim}, que corresponde a la simulación secuencial gaussiana. Estos algoritmos se modificaron para incluir la información blanda en dos modalidades. La primera consiste en sondajes artificiales, o pseudosondajes, creados por el usuario y que explicitan lo que el experto infiere con respecto a la continuidad y conectividad del cuerpo y/o de la roca caja. La segunda es usar direcciones variables, las cuales definen las direcciones de continuidad en el algoritmo. Estas direcciones se utilizan de dos maneras: orientando la búsqueda y el variograma según el campo de direcciones o bien, utilizándolas en una nueva definición de distancia no euclidiana que las incorpora. Los resultados obtenidos, para una prueba en un yacimiento real, muestran una mejora de entre 1% y 5% para el acierto total. Si bien, para dicho indicador es una leve mejora, sigue siendo una mejora, mas aún considerando que el yacimiento no posee plegamientos importantes que permitan explotar todo el potencial de estos algoritmos. De los algoritmos de estudiados, cabe destacar el comportamiento de la simulación utilizando anisotropía variable, pues mostró una mejora importante en la precisión y en el calce con el volumen real del cuerpo, aunque acompañada de una baja considerable en la cobertura. Esta característica hace de este método una buena alternativa para evaluar yacimientos de una manera conservadora, cuando la inversión se ve restringida y el riesgo debe ser el menor posible, es decir, no tanto alcance, pero con gran presición. El uso de estos nuevos algoritmos constituyen una nueva alternativa para el modelamiento, donde el usuario puede asistir con su conocimiento experto, la creación del modelo, sobre todo en yacimientos donde exista estructuras complejas que dificulten el proceso normal. La ventaja de éstos sobre otras técnicas que persiguen el mismo objetivo es lo intuitivo de las ideas detrás de las modificaciones y de la información adicional que se genera. Estos algoritmos contituyen parte del backend de una aplicación que permitirá al usuario interactuar con los datos para guiar el modelamiento de cuerpos geológicos.

Magíster en Minería / Ingeniero Civil de Minas / El presente trabajo corresponde al desarrollo de una herramienta que permite simular valores de una variable regionalizada considerando que tales valores tienen una variación sistemática en el espacio. En este contexto, se desarrolla una nueva herramienta de simulación consistente en un algoritmo de simulación Gaussiana secuencial con rechazo considerando una deriva de referencia como input, bajo la hipótesis que esta herramienta permite respetar tal deriva, obteniendo resultados representativos de la base de datos en cuanto a sus estadísticos de orden 1 (histograma) y orden 2 (variograma). La metodología del algoritmo comienza definiendo la secuencia de visitas de nodos a simular de manera aleatoria. Se acepta o rechaza el nodo simulado en base a la deriva de referencia considerando un rechazo determinístico o probabilístico y una tolerancia dinámica. Para cada nodo se considera una vecindad de búsqueda de datos condicionantes para la simulación y una vecindad de búsqueda de datos para el cálculo de una media local simulada. El algoritmo permite ajustar el número aceptable de rechazos, el tamaño de la vecindad de búsqueda de la media local, la tolerancia y el tipo de rechazo. Se presentan dos casos de estudio. El primero consiste en un ejemplo sintético de una coordenada con deriva lineal. En este primer caso se tiene que, a mayor tolerancia o mayor vecindad de búsqueda de la media local, los valores simulados se distribuyen con mayor dispersión en torno a la deriva de referencia. El segundo estudio de caso consiste en una zona de interés del yacimiento Compañía Minera Cerro Colorado donde se realiza el proceso de simulación en seis unidades de estimación considerando diecisiete sensibilizaciones de los parámetros del algoritmo más una simulación basada en Kriging Simple (SK) y otra basada en Kriging de residuos (BT). En el caso de presencia de deriva se obtiene en general mejores resultados con el algoritmo propuesto que con el SK o BT cuando la deriva se ve reflejada de manera clara en el variograma como en la unidad de estimación cuatro. Las estadísticas de validación en términos de desempeño de las simulaciones como estimación (coeficiente de determinación R2, pendiente de la regresión de datos reales versus simulados y error medio) y en términos de cuantificación de la incertidumbre de los datos originales (accuracy plot) mejoran en relación al SK y BT. De esta manera, la herramienta desarrollada ofrece una alternativa flexible que mejora los estadísticos de validación en comparación al enfoque tradicional frente a un escenario de simulación con presencia de deriva clara en el variograma.

Magíster en Minería / Durante los últimos años, el desarrollo y aplicación de algoritmos basados en estadísticas multipuntos ha ido en aumento. La presente tesis considera la incorporación de derivas externas, asociadas a la media local, al algoritmo Filter-Based Simulation (Filtersim). Para lograrlo se programó en Python el algoritmo original. Luego, una modificación interna al algoritmo permite incorporar derivas en la media. Para ello un nuevo ponderador, w0, mide el nivel de influencia que tendrá en la elección del patrón en el proceso de simulación. Dos casos de estudio son llevados a cabo. El primero busca reproducir el comportamiento sinusoidal que presentan los valores altos de una variable continua en la imagen del corte pulido de una roca sedimentaria. De ella se infiere el comportamiento estructural mientras que de otra zona conocida, con características similares, son extraídos los datos condicionantes y las derivas externas. Se simula a distintos grados de información por parte de la deriva y distinto número de datos condicionantes. Se deriva de lo anterior que un óptimo porcentaje de influencia de la deriva en la media oscila entre un 40% y un 60%, logrando caracterizar adecuadamente el comportamiento sinusoidal y la continuidad espacial de estas. De la variación en la cantidad de datos condicionantes se concluye que considerar entre un 0.5% y un 1.0% del dominio informado ayuda a mejorar los resultados, y las derivas obtenidas de ellos tienen un alto grado de información. Además, se reafirma la capacidad de los algoritmos multipuntos de reproducir estructuras complejas en contraste con la incapacidad de los métodos de estimación/simulación convencionales basados en el variograma. El segundo caso de estudio se aplica en el yacimiento de Compañía Minera Cerro Colorado, específicamente en el sector Sur-Oeste entre los bancos 2395 y 2355. Se busca reproducir la ley de cobre total de los bancos 2385 al 2355 extrayendo del banco 2395 la imagen de entrenamiento para la obtención de patrones mediante Filtersim. Las simulaciones se condicionan con los sondajes de exploración y se contrastan con los pozos de tronadura asumiendo su representación de la realidad. Se utiliza un 40% de influencia en la media local. Los resultados muestran una buena aproximación a la realidad al utilizar Filtersim convencional, y al incorporar derivas se disminuye la dispersión entre los resultados obtenidos, mejorando con esto la reproducción de las leyes. No se observa un beneficio claro y considerable a la hora de incorporar derivas en la media. Esto debido, posiblemente, a que el proceso geológico de formación asociado al yacimiento es mas bien masivo careciendo de estructuras geológicas con forma y continuidad bien marcadas. En específico se obtiene una mejor cuantificación de los recursos, a través de Curvas Tonelaje-Ley, al incorporar deriva local que sin ella. Así como también, la capacidad de reproducir el comportamiento de la deriva en los distintos bancos.

Orientador : Douglas Eduardo Zampieri / Tese (doutorado) : Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-18T18:48:56Z (GMT). No. of bitstreams: 1 Barbieri_Nilson_D.pdf: 9919500 bytes, checksum: 26df2a1a469601e74854c805db8b3025 (MD5) Previous issue date: 1993 / Resumo: O comportamento dinâmico de um veículo de quatro rodas é analisado através de um modelo com sete graus de liberdade para duas situações de carga do veículo, carregado e descarregado, e para dois diferentes tipos de via, asfáltica e pavimentada. Varia-se a velocidade do veículo entre 5 e 30 m/ s. Para encontrar as características dos elementos suspensivos usa-se a teoria do controle estocástico linear ótimo para minimizar um índice de desempenho cujas partes integrantes são: aceleração vertical do veículo (conforto), espaço de trabalho das suspensões (limitação estrutural), deslocamentos relativos entre os centros das rodas e a via (segurança) e forças de controle. A influência destas partes no índice de desempenho é relatada através de constantes de peso. Otimiza-se primeiramente um sistema de suspensão passiva com constantes de peso favorecendo a segurança. Estuda-se a influência da velocidade nos coeficientes de rigidez e amortecimento. Em seguida otimizam-se vários sistemas de suspensão com elementos passivos em paralelo com elementos ativos com realimentação de estado e de saída. Dois tipos de constantes de peso favorecendo o conforto e a segurança são usados. Otimizam-se também as constantes de peso das forças de controle e os elementos passivos das suspensões. Verifica-se o comportamento próprio de todos os sistemas através da obtenção dos autovalores. Compara-se o desempenho de todos os sistemas quando o veículo passa por uma via com perfil aleatório / Abstract: The dynamical behaviour of a four wheeled car is analysed by the use of a seven degrees of freedom model for two load conditions of the vehicle,loaded and unloaded, and for two different types of road, asphalt and pavement. The vehicle velocity is varied between 5 and 30m/s. To find the characteristics of the suspension elements stochastic optimallinear control theory is used through the minimization of the performance index whose constituent parts are: vehicle vertical acceleration (comfort), suspension working space (structurallimitation), road holding (safety) and the control force. The infiuence of these parts in the performance index are related through the weighting constants. The passive linear suspension system is firstly optimized with weighting constants favouring safety. The infiuence of the vehicle speed on the linear spring and damper coefficients are studied. Afterwards some systems with passive elements in parallel with active elements with state and output feedback are optimized. Two types of weighting constants favouring comfort and the safety are used. The weighting constants of the control forces and the values of the coefficients of passive suspensions are also optimized. The dynamical eigenbehaviour of the whose system is studied through the derivation of the eigenvalues. The performance of every system is compared when the vehicle is driving on a ramdom road / Doutorado / Doutor em Engenharia Mecânica

Orientador: Armando Zaupa Remacre / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Geociencias / Made available in DSpace on 2018-07-18T18:40:17Z (GMT). No. of bitstreams: 1 Paraizo_PauloLopesBrandao_M.pdf: 5427666 bytes, checksum: 9a5c378fa53b0dd82cff128b79ab3ab6 (MD5) Previous issue date: 1993 / Resumo: Neste trabalho foi utilizada uma metodologia geoestatística, - a simulação condicional-, em duas escalas de heterogeneidades geológicas do reservatório, tendo como objetivo poder medir as incertezas associadas a estas escalas num modelo simplificado de fluxo de fluidos. A primeira escala se refere à distribuição espacial das eletrofácies, que foi simulada com o algoritmo de componentes principais indicadores, e a segunda à distribuição espacial das propriedades petrofísicas, internamente à cada eletrofácies, tendo sido utilizado um algoritmo gaussiano seqüencial. Para cada escala de heterogeneidade foram selecionadas 10 imagens para a simulação de fluxo, onde se verificou, através de um modelo simples, que as incertezas dos parâmetros geológicos dependem não só da sua variabilidade, mas sim da forma como eles interagem com as condições de fluxo as quais o reservatório é submetido. Foi efetuado também um estudo do empilhamento vertical das eletrofácies, onde através de ferramentas simples como curvas de proporção vertical e variogramas, foi possível identificar padrões de comportamento vertical distintos para as várias zonas do reservatório em estudo. Esta abordagem forneceu uma visão abrangente do pacote sedimentar, e se mostrou útil para auxiliar a interpretação de caráter explotatório da área / Abstract: A geostatistica1 method - conditional simu1ation - was used in two sca1es of geo1ogica1 heterogeneity, with the objective of assessing the uncertainty of these sca1es in a simplified fluid flow model. The fust scale is associated with spatial distribuition of eletrofacies, simulated with a indicator principal component algorithm, and the second with the spatial distribuition of petrophysical properties for each e1etrofacies independent1y, simulated with a gaussian sequential algorithm. Ten realizations of each scale were chosen for fluid flow simulation, revealing that the uncertainties are not on1y associated with the variability of the images, but also with the interaction of these with the conditions of flow to which reservoir is submitted. A study of the vertical stacking of e1etrofacies was done with the use of simple too1s like proportion curves and variograms. Itwas possib1e to recognize different patterns of vertical stacking. The approach proved to be useful by providing a general view of the reservo ir, and helpful for the explotatory interpretation / Mestrado / Geologia do Petroleo / Mestre em Geoengenharia de Reservatorios

Memoria para optar al título de Ingeniero Civil de Minas / En la explotación de los yacimientos minerales en la actualidad, a menudo, se carece de una cuantificación de la incertidumbre en las unidades geológicas (UG), en donde principalmente están presentes modelos manuales típicamente determinísticos. Existen métodos para generar un mayor conocimiento de la distribución de las unidades geológicas en el espacio, que involucran una colección de técnicas geoestadísticas, con la finalidad de realizar un apropiado diseño minero con la información obtenida. El objetivo de esta memoria es modelar la incertidumbre geológica en un yacimiento minero mediante la utilización y comparación de siete modelos propuestos. Estos modelos se basan en estimar las proporciones de unidades geológicas con estadísticas geográficamente ponderadas mediante vecindades móviles y k datos más cercanos, para luego aplicar una simulación plurigaussiana, obteniendo de esta forma realizaciones de las unidades geológicas en el espacio. Para validar los modelos de simulación, se utiliza el método "jack knife", ocultando información de sondajes con datos conocidos de UG, para luego simular en dichos puntos del espacio, comparando los valores reales en los sitios de muestreo, contra los valores obtenidos en las realizaciones. Los resultados obtenidos muestran que la simulación plurigaussiana con los modelos propuestos es una buena herramienta para calcular la ocurrencia de las unidades geológicas (UG) en el espacio. La incorporación de modelos no estacionarios en la estimación de proporciones muestra una capacidad de adaptación a los posibles escenarios, ya que si bien, por una parte la simulación plurigaussiana permite un manejo de los contactos entre UG a través de una regla de truncación, los modelos de proporciones no estacionarios permiten un manejo de las transiciones más preciso, logrando contactos más abruptos o suaves en el caso de estadísticas geográficamente ponderadas, siendo éste un modelo más versátil que el de vecindad móvil y k datos más cercanos. Para los modelos de simulación, el conjunto de realizaciones muestran coincidencia visual en comparación a la estimación de probabilidades de UG, a pesar de haber conseguido una baja cantidad de realizaciones debido al costo computacional del método. Los resultados para la validación fueron positivos para los siete modelos de simulación utilizados, aunque los valores estuvieron por debajo de lo esperado, mostrando en muchos casos que para un cierto rango de probabilidad de ocurrencia de UG, el porcentaje de acierto estaba por debajo del rango, pero no muy lejano. Aun así se observan mejores resultados en cuanto a calidad y realismo de la simulación en comparación al modelo estacionario (que considera proporciones de UG constantes en el espacio) destacando por sobre el resto el modelo de proporciones calculadas con estadísticas geográficamente ponderadas. Finalmente, los modelos de simulación mostrados presentan una forma versátil para modelar la probabilidad de ocurrencia de UG, pudiendo extender estos modelos a la simulación conjunta con leyes.

Magíster en Minería / Diversos modelos geoestadísticos permiten simular tipos de roca o, más generalmente, unidades geológicas en los yacimientos, entre los cuales se cuenta el modelo de simulación plurigaussiana y el de simulación secuencial de indicadores. Estos modelos se basan en la definición de reglas de contacto entre unidades geológicas y/o modelos de correlación espacial de sus variables indicadores. El presente trabajo plantea entonces abordar los vínculos existentes entre las relaciones de contacto entre unidades geológicas y los variogramas de sus indicadores, y ver en qué medida estas relaciones de contacto pueden ser reproducidas al fijar los variogramas de indicadores, y recíprocamente, cómo los variogramas de indicadores son reproducidos al fijar relaciones de contacto. Más específicamente, interesa examinar la existencia o no de contactos, así también si el contorno de la unidad geológica en extensión, es continuo o no. Se cuenta con datos reales de cuatro yacimientos tipo pórfido. Para realizar el análisis se trabajó mediante dos etapas, la primera donde se despliegan los datos (litología, alteración y mineralización) del modelo de bloques, que es una interpretación geológica, y los datos de muestras de sondajes. Mediante el análisis de los mapas se deducen las propiedades esperadas de los variogramas de indicadores en el origen. Luego en la segunda etapa se determina las relaciones de contacto a partir de los variogramas experimentales de indicadores del modelo de bloques, donde además se corrobora las propiedades deducidas en la anterior etapa y, de esta forma, la consistencia del modelo interpretado, comparando casos con los variogramas de indicadores calculados a partir de las muestras de sondajes. En los resultados de este trabajo se muestra que existe una relación entre los variogramas de indicadores y los contactos geológicos (tipo de roca, alteración y mineralización), mediante el análisis del comportamiento en el origen de los variogramas directos y cruzados de indicadores se logra validar la consistencia entre un modelo de bloques interpretado y datos experimentales. Además, en la práctica, si solo se tiene datos de sondajes, es sencillo calcular variogramas de indicadores y así deducir propiedades que debería tener el modelo interpretado, lo cual ayuda en la consistencia al momento de construir dicho modelo.

Magíster en Minería / Usualmente los yacimientos son separados en distintos dominios geológicos, en los cuales las propiedades de las variables a simular son espacialmente y estadísticamente homogéneas. Sin embargo, cuando se consideran casos multivariables, la partición del yacimiento normalmente es llevada a cabo de acuerdo a los criterios de la variable cuyo estudio es más relevante, induciendo un criterio que no necesariamente se ajusta a las otras variables. Esto es particularmente problemático cuando las variables involucradas tienen distinta naturaleza, como lo son las leyes de elementos, que son resultados de procesos geológicos, y consumo de ácido, que es el resultado de procesos metalúrgicos dependientes de la naturaleza geológica. El objetivo principal de esta tesis es cuantificar la incertidumbre y variabilidad del consumo de ácido y de leyes de carbonato en un yacimiento mediante modelos geoestadísticos de simulación que utilizan la información de sondajes de exploración. Para la definición de dominios se propone una metodología que define los dominios de acuerdo a cada variable de interés, pudiendo ser iguales, traslapados o disjuntos para las distintas variables. Adicionalmente, se considera la dependencia espacial de los datos que pertenecen a distintos dominios, permitiendo aportar información a otros dominios, lo cual es particularmente significativo en zonas de transición entre dos dominios. El enfoque propuesto es aplicado a un caso de estudio, utilizando la base de datos del depósito Antucoya del norte de Chile, donde se simulan las leyes de carbonato y el consumo de ácido para distintas unidades mineralizadas. Distintos escenarios de la distribución de unidades mineralizadas son generados mediante simulaciones plurigaussianas, para obtener la variabilidad de las unidades en el yacimiento. Posteriormente, se realizan cosimulaciones multigaussianas de la ley de carbonato y del consumo de ácido, para poder cuantificar su incertidumbre. Estos resultados fueron validados mediante validación cruzada, donde no se observaron diferencias significativas, en términos globales, entre los distintos modelos de simulación puestos a prueba. Sin embargo, se observa diferencias locales en las simulaciones para los sitios que se encuentran cercanos al contacto entre los distintos dominios.