Os investimentos e o planejamento econômico na mineração requerem o entendimento, quantificação e avaliação de riscos na determinação de teores e tonelagens de minério. O planejamento de lavra e as operações em usinas de beneficiamento também requerem a estimativa de teores e tonelagens. A geoestatística proporciona as ferramentas necessárias para realizar essas estimativas utilizando adaptações das técnicas clássicas de regressão. O minério de ferro, assim como o de manganês, são exemplos de casos nos quais é necessário determinar múltiplas variáveis para a caracterização do minério, de modo tal, que os teores nos modelos de blocos dos depósitos satisfaçam os balanços de massa entre as frações granulométricas e a estequiometria para as espécies químicas dos dados originais. Esses sistemas são altamente complexos, apresentando múltiplas variáveis correlacionadas e sendo, portanto, apropriados para serem abordados com a utilização de técnicas geoestatísticas multivariadas. Apesar de a cokrigagem apresentar um estimador que desde o ponto de vista teórico é não tendencioso e que minimiza a variância do erro, nos casos de depósitos complexos, com um número elevado de variáveis relacionadas por várias somas constantes, apresenta uma série aspectos problemáticos: (i) estimativas fora do intervalo original das amostras e/ou negativas que precisam ser pós-processadas; (ii) a não satisfação, por parte dos teores estimados, dos balanços de massa e estequiométricos (somas constantes), sendo necessária a distribuição do erro ou carregá-lo em uma variável que seja determinada a partir dessas relações; (iii) dificuldades na modelagem da corregionalização que faz necessários softwares que consigam dar satisfação às condições de definição positiva impostas pelo Modelo Linear de Corregionalização e (iv) modelos variográficos que não se ajustam adequadamente aos variogramas experimentais diretos e cruzados. Essa tese aborda metodologias da geoestatística multivariada, alternativas à cokrigagem, no caso de um depósito de ferro, com um número elevado de variáveis presentes em diversas faixas granulométricas, correlacionadas espacialmente e que satisfazem várias relações de soma constante simultaneamente. Inicialmente, é aplicada a decomposição em Fatores de Autocorrelação Mínimos/Máximos (MAF), uma metodologia semelhante à decomposição em Componentes Principais (PCA), que propõe descorrelacionar as variáveis até um vetor de separação pequeno, geralmente coincidente com o espaçamento amostral ou com o alcance da primeira estrutura do variograma, evitando assim a modelagem da corregionalização, determinando cada fator de forma independente e posteriormente retro-transformando-o ao espaço original. Nessa tese os fatores são estimados por krigagem ordinária como uma metodologia de estimativa simples e rápida, porém aproximada, já que a transformação MAF inclui uma etapa de normalização das informações originais, o que constitui uma transformação não linear. A metodologia fornece resultados adequados, não sendo necessária a utilização de softwares específicos para a modelagem da corregionalização já que o problema se remete à modelagem de variogramas de variáveis independentes, e, ao incluir uma etapa de normalização e retro-transformação ao espaço original, as estimativas permanecem dentro do intervalo original das amostras, não havendo a necessidade de pós-processá-las. Compara-se com a metodologia clássica da krigagem das componentes principais (PCA), que descorrelaciona as variáveis só para o vetor de separação nulo, com resultados superiores para os fatores MAF. Porém, não dá uma solução ao problema do fechamento dos balanços, devendo-se também, distribuir o erro ou carregá-lo numa variável que é determinada por diferença, para satisfazê-los. Em segunda instância, implementa-se a cokrigagem de razões-logarítmicas aditivas (alr), sendo uma metodologia desenvolvida para dados composicionais (que apresentam uma soma constante). Essa metodologia fornece resultados superiores aos obtidos por cokrigagem direta das variáveis originais, com estimativas dentro do intervalo original das amostras e uma satisfação dos balanços considerados para a totalidade dos valores determinados. A consideração mais importante dessa metodologia é que as variáveis que constituem uma composição regionalizada, como é o caso das variáveis consideradas nos depósitos de ferro, têm um espaço amostral que está restringido ao simplex no qual a soma constante é satisfeita. Porém, não oferece uma solução para o problema da modelagem da variabilidade espacial conjunta (modelo linear de corregionalização). No entanto, diminui numa unidade a ordem de magnitude dos sistemas de cokrigagem sendo, portanto, mais simples de modelar. / Determination of grades and tonnages, directly affect the comprehension, quantification and risks evaluation in the investment and economic planning for mining projects. Adapting the classic regression techniques, Geostatistics provide the needed methodologies for these determinations. Iron and manganese ores, are examples of cases in which multiple variables need to be determined for ore characterization, simultaneously satisfying the original mass balances and stoichiometry among granulometric fractions and chemical species, respectively. These highly complex systems, with multiple correlated variables, are appropriate for the utilization of multivariate geostatistics. Although the cokriging provides an unbiased estimator that minimizes the error variance, in the case of complex deposits, with a high number of correlated variables, related by multiple constrained sums, it leads to some problematic aspects: (i) estimates outside the original data interval of values and/or negative values, that need to be pos-processed; (ii) estimates that do not satisfy the mass balances and stoichiometry, with the error having to be distributed among variables or assigned to a single variable determined by difference; (iii) impositions of the Linear Model of Corregionalization in presence of multiple variables, need to be achieved by the utilization of specific softwares and (iv) variogram models that do not adhere to experimental direct and cross-variograms. In this thesis, methodologies of multivariate geostatistics, other than cokriging, are evaluated in the frame of an iron ore deposit, with multiple correlated variables, present in various granulometric fractions, simultaneously satisfying diverse closed sum constraints. In first place, decomposition in Minimum/Maximum Autocorrelation Factors (MAF), a methodology similar to Principal Components decomposition (PCA) is performed. This methodology decorrelates variables up to a small separation vector, generally coincident with sampling spacing or with the range of the first structure of the variogram, thus, allowing estimating each factor individually, avoiding modeling the corregionalization. In this thesis, the MAF factors are estimated through ordinary kriging as it is simple and easily implemented, although it provides approximated estimates because of the nscore transformation of the original data that is embedded in the MAF decomposition and back-transformation, which is a non-linear transformation. It provides adequate results, without specific softwares needed for modeling the corregionalization, because the problem is simplified to modeling and estimating individual variables. It is compared with the classical methodology of kriging the Principal Components (PCA) which decorrelates the variables only for a separation vector equal to zero, with superior results for MAF decomposition. The nscore transformation and back-transformation, leads to estimates restricted to the original samples interval, eliminating the pos-processing step. Nevertheless, MAF decomposition does not provide a solution to the closure of mass balances and stoichiometry and the error has to be distributed among variables or assigned to a single variable obtained by difference, as in the case of cokriging of the original data. Secondly, additive log-ratio (alr) cokriging is implemented, being a methodology developed for compositional data (with closed constant sum). This methodology provides better results when compared to the ones obtained by cokriging of the original data, with all estimates within the original data values interval and satisfying the considered balances. The most important consideration of this methodology is that the variables that conform a regionalized composition, as it is the case of variables considered in iron ore, have a sample space that is restricted to the simplex in which the constant sum condition is satisfied. However, it does not provide a solution for modeling the spatial joint correlation (Linear Model of Corregionalization). But, it leads to cokriging systems that are one unit smaller than the original ones and consequently easier to model.

o processo de avaliação de um depósito mineral é compreendido de duas grandes etapas: a coleta e análise de amostras que fornecerão informações importantes que servirão de base para as estimativas e a estimativa propriamente dita, onde serão aplicados métodos matemáticos aos dados obtidos pela amostragem. Ambas as etapas apresentam diversos problemas e empecilhos. Em relação à amostragem, o método usualmente aplicado é a perfuração rotativa com recuperação de testemunhos. Este método possui um custo bastante elevado, podendo um furo de poucos metros custar centenas de dólares e levar dias para ser feito. Outro problema associado com este método é a probabilidade de não se ter uma elevada recuperação do material, isto é, perda de material durante a execução da perfuração. Já na etapa de análise do material recuperado, uma das dificuldades mais freqüentes é a impossibilidade de se avaliar a flutuabilidade do valor estimado, ou seja, geralmente se estima um valor de atributo para a amostra e se considera esse valor único e real para ela, porém, sabe-se que, com análise em laboratório, é impossível se estimar o valor exato do atributo de uma amostra e, portanto, deve-se realizar uma série de análises e efetuar o calculo da média e desvio padrão do valor do atributo da amostra, porém, muitas vezes não se tem material físico (amostras) disponível para um grande número de análises. Quanto à avaliação do depósito, o principal problema está relacionado á impossibilidade de se considerar essa flutuabilidade do valor do atributo das amostras no momento da estimativa. A grande maioria dos métodos de estimativa usada para avaliação de depósitos carece de mecanismos que permitam a incorporação da incerteza do valor do atributo. Portanto, os objetivos desse trabalho foram: (i) encontrar uma maneira de aprimorar a coleta de informações de um depósito mineral e em contrapartida facilitar a estimativa da flutuabilidade do valor dos atributos das amostras e (ii) encontrar um método de avaliação de depósitos que permita a incorporação dessa incerteza do valor das amostras. Para solucionar os problemas apresentados pela amostragem, esta dissertação apresenta a perfilagem geofísica como uma ferramenta versátil e de fácil manuseio em campo, que permite a obtenção de um grande número de informações de caráter indireto e por meio de ferramentas matemáticas ainda permite o cálculo da flutuação do valor do atributo. Em relação à avaliação de depósitos, será utilizada a krigagem dos indicadores, método de estimativa não linear que permite a incorporação da incerteza dos valores dos atributos de uma maneira rápida no momento do cálculo das áreas de interesse. Descobriu-se que a aplicação da perfilagem geofísica foi bastante eficiente e prática contanto que a equipe responsável tenha o conhecimento das ferramentas utilizadas e como é a interação entre a ferramenta e o meio geológico. A eficácia do método geoestatístico foi comprovada através da comparação do resultado da krigagem dos indicadores que incorporou dados secundários provenientes da perfilagem geofísica com o resultado de uma estimativa usualmente feita em avaliação de depósitos baseada somente em dados primários (sondagem com recuperação de testemunhos). / The process of evaluation of a mineral deposit is constituted of two great stages: the collection and analysis of samples that will supply important information, which it will be the base for the evaluation and the estimating of the deposit where they will be applied mathematical methods on the data obtained from the sampling. 80th stages present several problems and difficulties. In relation to the sampling, the method usually applied it is the rotary drilling with ore recovering. This method is expensive, a hole of few meters will probably cost hundreds of dollars and take days to be done. Another problem associated with this method is the probability of not having a high recovery of the material, that is, materialloss during the execution of the drilling. About the analysis of the recovered material, one of the most frequent difficulties is the impossibility of evaluating the floatability of the estimated value, in other words, is usually estimated an attribute value for the sample and it is considered unique and real value for it, however, it is known that, with laboratory analysis, it is impossible to obtain the exact value of the attribute of a sample and, therefore, it should be done a series of analyses to calculate the average and standard deviation of the value of the attribute of the sample, however, most of times a great quantity of physical material (samples) are not available for a great number of analyses. For the evaluation of the deposit, the main problem is related to the impossibility of considering that floatability of the value of the attribute of the samples in the moment of the estimate. The great majority of the estimate methods used for evaluation of deposits lacks of mechanisms that allow the incorporation of the uncertainty of the value ofthe attribute. Therefore, the objectives of this research project are: to find a way to optimize the collection of information of a mineral deposit and in consequence to facilitate the estimate of the floatability of the value of the attributes of the samples and to find a method of evaluation of deposits that allow the incorporation of that uncertainty of the value of the samples. To solve the problems presented by the sampling, this dissertation presents the borehole logging as a versatile tool and of easy to handle in the field, which allows to obtain a great number of information with indirect features and it allows the calculation of the floatability of the value of the attribute through mathematical tools. In relation to the evaIuation of the deposits it will be use the indicator kriging. A non-linear method that allows the incorporation of the uncertainty of the values of the attributes in a fast way in the moment of the calculation of the areas of interest. It was discovered that the application of the borehole logging with indirect methods was quite efficient and practice as long as the responsible team has the knowledge of the use of the tools and about the interaction between the logging results and the geologicaI formations. The effectiveness of the geostatistic method was proven through the comparison of the result of the indicator kriging executed with secondary obtained with borehole logging with the result of an estimate done usuaIly in evaluation of deposits based only in primary data (survey ore recovery).

CAJAZEIRA, J. P. Caracterização e variabilidade espacial de Atributos Físicos em um Argissolo Amarelo no estado do Ceará. 2007. 80 f. Dissertação (Mestrado em Agronomia/Solos e Nutrição de Plantas) - Centro de Ciências Agrárias, Universidade Federal do Ceará, Fortaleza, 2007. / Submitted by Francisco Lacerda (lacerda@ufc.br) on 2014-09-11T12:58:24Z No. of bitstreams: 1 2007_dis_jpcajazeira.pdf: 7553068 bytes, checksum: 2ef3fd3c2b8712d8c5c627a777afcf8c (MD5) / Approved for entry into archive by José Jairo Viana de Sousa(jairo@ufc.br) on 2014-09-11T23:21:03Z (GMT) No. of bitstreams: 1 2007_dis_jpcajazeira.pdf: 7553068 bytes, checksum: 2ef3fd3c2b8712d8c5c627a777afcf8c (MD5) / Made available in DSpace on 2014-09-11T23:21:03Z (GMT). No. of bitstreams: 1 2007_dis_jpcajazeira.pdf: 7553068 bytes, checksum: 2ef3fd3c2b8712d8c5c627a777afcf8c (MD5) Previous issue date: 2007 / The objective of this research was characterize and determine to spatial variability of physical attributes of a Alfisols, using methods geostatistical e and making of surface maps for the method of Kriging, in two soil depths. They were made soil collections in regular intervals of 2.5 m, in an area of 400 m2, being 81 points in each one of the layers (0 – 0.2 m and 0.2 – 0.4 m), to add 162 samples that were submitted to the physical determinations and the obtained data analyzed by the descriptive statistics to verify the normality of the frequency distribution and later on submitted to the geostatistical. The model that better it was adjusted to the curve of the semivariogram, in most of the variables, it was the spherical. However, for the variable density of the particles and aggregate with diameter 2.00-1.00 mm, in the depth of 0 – 0.2 m and flocculation degree, bulk density and soil-water content with tension 1.5 MPa, in the depth of 0.2 – 0.4 m the most appropriate model that better adjusted it was the exponential. The nugget ratio for the flocculation degree, bulk density and density of the particles was of 100%, characterizing “random variance”, confirming the independence of the data. The largest ranges were observed in the depths of 0.2 – 0.4 m for lime percentage (51.00 m) and percentage of sand (35.06 m). The values of the bulk density were considered relatively high being of 1.69 and 1.65 Mg.m-3 for 0 – 0.2 and 0.2 – 0.4 m, respectively. There was a larger humidity retention in the inferior layers. / O objetivo do trabalho foi caracterizar e determinar a variabilidade espacial de atributos físicos de um Argissolo Amarelo, utilizando métodos geoestatísticos e confecção de mapas de superfície pelo método da Krigagem, em duas camadas de solo. Foram realizadas coletas de solo em intervalos regulares de 2,5 m, numa área de 400 m2, sendo 81 pontos em cada uma das camadas (0 – 0,2 m e 0,2 – 0,4 m), totalizando 162 amostras que foram submetidas às determinações físicas e os dados obtidos analisados pela estatística descritiva para verificar a normalidade da distribuição de freqüência e posteriormente submetidos a geoestatística. O modelo que melhor se ajustou à curva do semivariograma, na maioria das variáveis, foi o esférico. No entanto, para a variável densidade das partículas e agregados com diâmetro 2,00-1,00 mm, na profundidade de 0 – 0,2 m e grau de floculação, densidade do solo e umidade obtida com tensão 1,5 MPa, na profundidade de 0,2 – 0,4 m o modelo que melhor se ajustou foi o exponencial. O efeito pepita relativo para o grau de floculação, densidade do solo e densidade das partículas foi de 100%, caracterizando o efeito pepita puro, confirmando a aleatoriedade dos dados. Os maiores alcances foram observados nas profundidades de 0,2 – 0,4 m para percentagem de silte (51,00 m) e percentagem de areia (35,06 m). Os valores da densidade do solo foram considerados relativamente altos sendo de 1,69 e 1,65 Mg.m-3 para 0 – 0,2 e 0,2 – 0,4 m, respectivamente. Houve uma maior retenção de umidade nas camadas inferiores.

Há, na mineração, no que se refere à estimativa de teores, inúmeras situações em que as variáveis de interesse exibem correlação entre si. Neste caso, devem ser consideradas as correlações no processo a fim de se obter um resultado consistente com a estatística de correlação entre as variáveis presentes no depósito. A maneira clássica de abordar problemas desta natureza é por meio do uso da cokrigagem, que exige a modelagem das covariâncias diretas e cruzadas das variáveis. Esta metodologia é muito trabalhosa quando envolve mais de duas variáveis. Alternativas aparecem na forma de métodos de fatoração, em que as variáveis originais correlacionadas são transformadas em fatores descorrelacionados. Nesta linha, se tem o método de Mínimos-Máximos Fatores de Autocorrelação (Minimum-Maximum Autocorrelation Factors, MAF), no qual, por meio de uma transformação linear nos atributos correlacionados, leva-os a um espaço em que estejam descorrelacionados, possibilitando que a estimativa seja realizada individualmente pela krigagem, desta forma, evitando o uso do modelo linear de corregionalização. Nesta dissertação, é analisada a aplicabilidade do método de Mínimos-Máximos Fatores de Autocorrelação à estimativa de teores em depósitos multivariados, bem como analisadas suas vantagens e limitações. A técnica alternativa é comparada com o método de coestimativa clássico, cokrigagem, com as análises em diferenças de estimativas obtidas por meio de cada metodologia. Os resultados obtidos mostraram, nos estudos de caso realizados, que a solução MAF proveu resultados de maneira mais rápida e simples sem perda de precisão em relação à cokrigagem. / There are in mining, as regards the grade estimation, numerous situations in which the variables of interest exhibit correlation. In this case, the correlation must be considered in the process to obtain a result consistent with the statistical correlation between the variables present in the deposit. The classic way of addressing such problems is through the use of cokriging, which requires the modeling of direct and cross-covariances of the variables. This methodology is very troublesome when it involves more than two variables. Alternative appear as factoring methods, where the original correlated variables are transformed into uncorrelated factors. In this line, you have the method Minimum - Maximum Autocorrelation Factors (Minimum - Maximum autocorrelation Factors, MAF) where through a linear transformation on attributes correlated, leads them to a space where they are uncorrelated, allowing that the estimation be performed individually through kriging, thus avoiding the use of the linear model of coregionalization. In this dissertation is analyzed the applicability of the method of Minimum - Maximum Autocorrelation Factors for grade estimation in multivariate deposits, as well as its advantages and limitations. The alternative technique is compared to the classic coestimation method, cokriging, by analyzing differences in the estimates obtained by each method. The results showed, in the case studies that the MAF solution provided results more quickly and easily without loss of accuracy regarding to the cokriging.

Os modelos geológicos 3D associado as suas compartimentações geomecânicas e hidrogeológicas não são atualmente usados na maioria das minas, particularmente na indústria mineira brasileira. Isso ocorre muitas vezes devido ao investimento em construção da base de dados de dados geomecânicos se inicia em etapas mais madura dos projetos e por consequência o seu volume não é suficiente para ser espacializada. A ausência de informações sobre a qualidade dos maciços rochosos afeta os projetos de mineração, especialmente as mais sensíveis à incerteza associada aos parâmetros geotécnicos. Esta falta de modelos geomecânicos deve-se principalmente à escassa disponibilidade de dados geotécnicos e à consequente dificuldade de interpola-los e discretizá-los em vários locais de depósito. Para superar esta escassez de dados, este estudo propõe a construção de modelo geomecânico 3D associado ao modelo geológico via modelagem indireta de classes de consistência de rochas, usando dados derivados da descrição tátil-visual em testemunhos de sondagem. Estes dados foram agrupados em sete categorias de consistência de rochas e modelados espacialmente através de krigagem de indicadores categóricos. Cada uma destas categorias está associada a uma gama de valores possíveis de resistência à compressão uniaxial. Combinando cada probabilidade de categoria estimada e seu intervalo de possíveis valores de resistência, é possível atribuir a cada bloco uma resistência à compressão uniaxial média. Os valores numéricos podem ser estimados pela média dos histogramas categóricos (classe probabilidade x centro do intervalo de classe) para preencher o modelo de bloco com valores de resistência à compressão uniaxial. A dissertação também apresenta o esquema de interpolação de anisotropia dinâmica que facilita o processo de modelagem da continuidade das camadas em mineralizações dobradas. Este método basicamente modela a direção do mergulho de cada ponto mapeado ou interpolado em mapas de linhas de forma de cada tipo de rocha, permitindo que o variograma e o elipsóide de busca girem adequadamente em cada posição dentro do domínio. O modelo final foi verificado dentro de regiões onde foram realizados ensaios laboratoriais de compressão uniaxiais. A metodologia se mostrou eficiente e foi ilustrada em um grande depósito de minério de ferro no Brasil, apresentando propostas de aplicações, como mapeamento de zonas de instáveis baseadas na resistência e na relação entre descontinuidade e direção dos taludes. / Geomechanical 3D models are not currently used in most mines particularly at the Brazilian mining industry. The absence of rock mass quality information affects mining design specially the ones more sensitive to the uncertainty associated with geotechnical parameters. This lack of geomechanical models is mainly due the scarce availability of geotechnical data and the consequent difficulty to interpolate and discretize them at various deposit locations. To overcome this data shortage, this study proposes building 3D geomechanical models via indirect modelling of rock consistency classes using data derived from tactile-visual description at borehole cores. These data were grouped into seven rock consistency categories and spatially modelled via categorical indicator kriging. Each of these categories are associated with a range of uniaxial compressive resistance possible values. Combining each estimated category probability and its range of possible resistance values allow to assign to each block an average uniaxial compressive strength. Numerical values can be estimated by averaging these categorical histograms (class probability x center of the class interval) to fill the block model with compressive resistance values. The dissertation also presents the dynamic anisotropy interpolation scheme which facilitates the modelling process in folded mineralizations. This method basically models the dip direction of each rock type allowing the variogram and search ellipsoid rotate adequately at each position within the domain. The final model was checked within regions of the model where compressive uniaxial lab tests were conducted. The methodology proved to be efficient and was illustrated at a major iron ore deposit in Brazil, presenting applications proposals, such as mapping unstable zones based on resistance and the relationship between discontinuity and slope direction.

Orientador: Gener Tadeu Pereira / Resumo: O solo submetido as operações agrícolas intensa e repetitivas RESUMO - O fósforo é considerado um dos nutrientes que limitam o desenvolvimento da cultura de cana-de-açúcar em solos brasileiros. Além disso, apresenta grande variabilidade espacial devido aos atributos que influenciam sua adsorção ao solo. O objetivo do trabalho foi mapear o fósforo adsorvido (Pads)com o auxílio da susceptibilidade magnética (SM) e definir a proporção amostral para a caracterização da dependência espacial do Pads. Foram coletadas 371 amostras na profundidade de 0-0,25 m com separação entre as amostras de 145 m a 174 m, abrangendo uma área total de aproximadamente 870 ha. Em cada ponto amostral determinou-se a suscetibilidade magnética (SM), o fósforo adsorvido (Pads). Os dados foram analisados a partir da estatística descritiva e da geoestatística. Na utilização da técnica da cokrigagem de Pads com o auxílio da SM sendo que uma a cada 2.6, assim foram estabelecidas as seguintes proporções amostrais (3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 ha) e foram obtidos os mapas de cokrigagem. A acurácia das estimativas foi realizada utilizando a validação externa e a Raiz Quadrada do Erro Médio (RMSE). Os valores de alcance para as variáveis estudadas variaram de 1861,78 a 2081,68m. A partir desse resultado conclui-se que o Pads apresenta dependência espacial com a SM, indicando seu uso no mapeamento e na identificação de áreas com diferentes potenciais de adsorção de fósforo, e interpolação com co-kri... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: Phosphorus is considered one of the nutrients that limit the development of sugarcane cultivation in Brazilian soils. In addition, it is an element that presents great spatial variability due to the attributes that influence it is adsorption to the soil. The objective of this study was to map the adsorbed phosphorus with the aid of magnetic susceptibility and to define the sample proportion for the characterization of the spatial dependence of the Pads. A total of 371 samples were collected at a depth of 0-0.25 m with a separation between samples of 145 m and 174 m, covering a total area of approximately 870 ha. At each sampling point the magnetic susceptibility (SM), the adsorbed phosphorus (Pads) was determined. Data were analyzed from descriptive statistics and geostatistics. In the use of the cokriging technique of Pads with the aid of SM, one of each being 2.6, the following sample proportions were established (3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 ha), and the cokriging maps were obtained. The accuracy of the estimates was performed using the external validation and the Square Root of the Mean Error (RMSE). The range values for the studied variables ranged from 1861.78 to 2081.68m. From this result, it is concluded that the Pads presents spatial dependence with SM, indicating their use in mapping and identifying areas with different phosphorus adsorption potentials, and interpolation with co-kriging using the SM works up to the sample proportion 6 points per hecta... (Complete abstract click electronic access below) / Mestre

Na maioria das minas, até poucos anos, só havia preocupação em estimar-se o teor de uma ou no máximo duas espécies químicas de interesse, por exemplo, teor de cobre, ouro ou ferro. Hoje, existem padrões de qualidade em produtos e concentrados que exigem um rigoroso controle de outros elementos além do metal/mineral minério. É comum ser necessário a estimativa de múltiplos elementos, possivelmente correlacionados e em algumas vezes com uma combinação de teores que deva fechar uma soma constante (por exemplo, 100%). Dados que somam uma constante são conhecidos como dados composicionais (CODA), carregando consigo informação relativa e não absoluta. Esta situação (soma de uma constante) condiciona a que as estimativas também devam fechar uma constante. As metodologias clássicas, como é o casso da krigagem e a cokrigagem ordinária, podem ser apropriadas para obter a melhor estimativa local do ponto de vista teórico. Mas, cada variável é estimada separadamente (no caso OK) com seus respectivos parâmetros de medidas de continuidade espacial, o que leva a obter pesos diferentes para cada atributo e em contraposto gera problemas para o fechamento final das estimativas. Novas técnicas estatísticas podem ser adaptadas ao estudo de CODA que permitam garantir, além de resultados coerentes, o fechamento após estimativas. Esta dissertação utilizou as transformações de razões logarítmicas isométricas (ilr) que permite projetar o espaço amostral dos dados originais (espaço euclidiano) ao espaço simplex de D-1 partes da composição (SD). Os resultados mostraram-se satisfatórios, já que além de obter bons resultados, foi garantido o fechamento após estimativas. Esses resultados foram comparados com os resultados obtidos na krigagem tradicional utilizando os mesmos critérios de interpolação. Ambas as metodologias produziram resultados semelhantes em termos de qualidade de estimativa, mostrando-se ser uma metodologia alternativa para estimativa em depósitos multivariados. Uma desvantagem nesta transformação é o uso de transformações de tipo não lineares combinada com interpoladores lineares (OK), já que quando fazemos o processo de retro-transformação (ilr) estamos expostos a um viés na média estimada. Para solucionar esse problema, foi utilizada a simulação das transformações (ilr) que não gera viés, já que o dado simulado é obtido a partir de uma tiragem randômica de função de probabilidade local. Comparou-se então, os resultados do e-type das simulações (ilr) com os resultados da krigagem (ilr) para observar a influência do viés da média aritmética localmente. Os resultados mostraram que, embora esse viés exista teoricamente, para este estudo de caso não foi representativo, ou seja, a média local entre o e-type das simulações (ilr) e a krigagem (ilr) não é maior que 5%. / In most of the mines, until a few years ago, there was only concern in estimating the content of one or at most two chemical species of interest, such as percentage of copper, gold or iron. Today, there are standards of quality in products and concentrates, which require careful control of other elements besides metal / ore minerals. It is common to estimate multiple elements, possibly correlated and sometimes with a combination of contents which must close a constant sum (e.g. 100%). Data that add up to a constant are known as compositional data (CODA), carrying information relative and not absolute. This (sum of a constant) determines that the estimates should also close a constant. The classical methods, such as kriging and ordinary cokriging, may be appropriate for the best local estimate. However, each variable is estimated separately (in case OK) with their specific parameters of spatial continuity, which leads to obtain different weights for each attribute and fails to obtain estimates that satisfy the constant sum constraint. New statistical techniques can be adapted to the study of CODA that guarantee consistent results after the closing estimates. This dissertation uses the isometric transformations of logarithmic ratios (ilr) that allows to transform the sample space of the original data (Euclidean space) into the simplex space with D-1 parts of the composition (SD). The results were satisfactory providing closed sums after estimates. These results were compared with the results obtained using the traditional kriging interpolator. Both methods produced similar results in terms of quality of the estimation, proving to be an alternative methodology to estimate multivariate deposits. A disadvantage in this transformation is the use of nonlinear transformations combined with linear interpolation (OK), since when doing this process of retro-transformation (ilr) can lead to a bias in the estimated average. To solve this problem, it was used simulation of the transformation (ilr) and this solution do not generates bias, since the simulated data is obtained from randomly drawing of a local probability function. The average of all simulations wase compared with kriging to observe the impact of the arithmetic bias. The results showed that, although this bias exists theoretically, for this case study it was not significant, i.e., the difference from the local average of the e-type simulations (ilr) and kriging (ilr) is not greater than 5%.

o processo de avaliação de um depósito mineral é compreendido de duas grandes etapas: a coleta e análise de amostras que fornecerão informações importantes que servirão de base para as estimativas e a estimativa propriamente dita, onde serão aplicados métodos matemáticos aos dados obtidos pela amostragem. Ambas as etapas apresentam diversos problemas e empecilhos. Em relação à amostragem, o método usualmente aplicado é a perfuração rotativa com recuperação de testemunhos. Este método possui um custo bastante elevado, podendo um furo de poucos metros custar centenas de dólares e levar dias para ser feito. Outro problema associado com este método é a probabilidade de não se ter uma elevada recuperação do material, isto é, perda de material durante a execução da perfuração. Já na etapa de análise do material recuperado, uma das dificuldades mais freqüentes é a impossibilidade de se avaliar a flutuabilidade do valor estimado, ou seja, geralmente se estima um valor de atributo para a amostra e se considera esse valor único e real para ela, porém, sabe-se que, com análise em laboratório, é impossível se estimar o valor exato do atributo de uma amostra e, portanto, deve-se realizar uma série de análises e efetuar o calculo da média e desvio padrão do valor do atributo da amostra, porém, muitas vezes não se tem material físico (amostras) disponível para um grande número de análises. Quanto à avaliação do depósito, o principal problema está relacionado á impossibilidade de se considerar essa flutuabilidade do valor do atributo das amostras no momento da estimativa. A grande maioria dos métodos de estimativa usada para avaliação de depósitos carece de mecanismos que permitam a incorporação da incerteza do valor do atributo. Portanto, os objetivos desse trabalho foram: (i) encontrar uma maneira de aprimorar a coleta de informações de um depósito mineral e em contrapartida facilitar a estimativa da flutuabilidade do valor dos atributos das amostras e (ii) encontrar um método de avaliação de depósitos que permita a incorporação dessa incerteza do valor das amostras. Para solucionar os problemas apresentados pela amostragem, esta dissertação apresenta a perfilagem geofísica como uma ferramenta versátil e de fácil manuseio em campo, que permite a obtenção de um grande número de informações de caráter indireto e por meio de ferramentas matemáticas ainda permite o cálculo da flutuação do valor do atributo. Em relação à avaliação de depósitos, será utilizada a krigagem dos indicadores, método de estimativa não linear que permite a incorporação da incerteza dos valores dos atributos de uma maneira rápida no momento do cálculo das áreas de interesse. Descobriu-se que a aplicação da perfilagem geofísica foi bastante eficiente e prática contanto que a equipe responsável tenha o conhecimento das ferramentas utilizadas e como é a interação entre a ferramenta e o meio geológico. A eficácia do método geoestatístico foi comprovada através da comparação do resultado da krigagem dos indicadores que incorporou dados secundários provenientes da perfilagem geofísica com o resultado de uma estimativa usualmente feita em avaliação de depósitos baseada somente em dados primários (sondagem com recuperação de testemunhos). / The process of evaluation of a mineral deposit is constituted of two great stages: the collection and analysis of samples that will supply important information, which it will be the base for the evaluation and the estimating of the deposit where they will be applied mathematical methods on the data obtained from the sampling. 80th stages present several problems and difficulties. In relation to the sampling, the method usually applied it is the rotary drilling with ore recovering. This method is expensive, a hole of few meters will probably cost hundreds of dollars and take days to be done. Another problem associated with this method is the probability of not having a high recovery of the material, that is, materialloss during the execution of the drilling. About the analysis of the recovered material, one of the most frequent difficulties is the impossibility of evaluating the floatability of the estimated value, in other words, is usually estimated an attribute value for the sample and it is considered unique and real value for it, however, it is known that, with laboratory analysis, it is impossible to obtain the exact value of the attribute of a sample and, therefore, it should be done a series of analyses to calculate the average and standard deviation of the value of the attribute of the sample, however, most of times a great quantity of physical material (samples) are not available for a great number of analyses. For the evaluation of the deposit, the main problem is related to the impossibility of considering that floatability of the value of the attribute of the samples in the moment of the estimate. The great majority of the estimate methods used for evaluation of deposits lacks of mechanisms that allow the incorporation of the uncertainty of the value ofthe attribute. Therefore, the objectives of this research project are: to find a way to optimize the collection of information of a mineral deposit and in consequence to facilitate the estimate of the floatability of the value of the attributes of the samples and to find a method of evaluation of deposits that allow the incorporation of that uncertainty of the value of the samples. To solve the problems presented by the sampling, this dissertation presents the borehole logging as a versatile tool and of easy to handle in the field, which allows to obtain a great number of information with indirect features and it allows the calculation of the floatability of the value of the attribute through mathematical tools. In relation to the evaIuation of the deposits it will be use the indicator kriging. A non-linear method that allows the incorporation of the uncertainty of the values of the attributes in a fast way in the moment of the calculation of the areas of interest. It was discovered that the application of the borehole logging with indirect methods was quite efficient and practice as long as the responsible team has the knowledge of the use of the tools and about the interaction between the logging results and the geologicaI formations. The effectiveness of the geostatistic method was proven through the comparison of the result of the indicator kriging executed with secondary obtained with borehole logging with the result of an estimate done usuaIly in evaluation of deposits based only in primary data (survey ore recovery).

Orientador: Célia Regina Lopes Zimback / Banca: Paulo Milton Barbosa Landim / Banca: Alessandra Fagioli da Silva / Banca: Bruna Soares Xavier de Barros / Banca: Diego Augusto de Campos Moraes / Resumo: O mapa de solos é uma ferramenta essencial para o planejamento de uso da terra e estudos que envolvem aspectos ambientais relativos a esse importante recurso natural. Técnicas quantitativas e ferramentas de geoprocessamento têm sido aliadas à interpretação dos processos pedogenéticos para possibilitar a elaboração de mapas mais precisos, obtidos por processo mais rápido e menos oneroso. Dentre os modelos aplicados, os denominados modelos híbridos empregam variáveis auxiliares preditoras e autocorrelação espacial, para viabilizar a predição de atributos de solo em locais não amostrados. A iniciativa para mapeamento digital do solo em escala mundial - GlobalSoilMap.net atua no sentido de disponibilizar representações globais de atributos de solo, elaboradas por meio da aplicação de modelo híbrido em dados legados de solos, realizando a prática do Mapeamento Digital de Solos (MDS). Com base nesse princípio, esse trabalho baseou-se na hipótese de que a aplicação da técnica híbrida regressão-krigagem, utilizando dados legados de levantamento de solo e covariáveis de relevo e sensoriamento remoto proveem mapa de atributos de solo representativos de uma área da Cuesta de Botucatu. O modelo foi aplicado localmente, a duas profundidades, para representação contínua do Índice de Avermelhamento (IAV), saturação de bases (V%), teor de areia, teor de argila, CTC e pH dos solos da Fazenda Experimental Edgárdia, para a qual são disponíveis dados de levantamento de solo. As covariáveis ... / Abstract: The soil map is an essential tool for land use planning and studies related to environmental aspects of this important natural resource. Quantitative techniques and geoprocessing tools are currently combined with the interpretation of pedogenic processes to enable the development of more accurate maps obtained by faster and less costly process. Among the models applied to it, the hybrid models employ predictive auxiliary variables and spatial autocorrelation, to enable the prediction of soil attributes in unsampled locations. The digital soil mapping worldwide project - GlobalSoilMap.net acts in order to provide global representations of soil attributes developed through the application of hybrid model in legacy soil data, performing the practice of Digital Soil Mapping (MDS). This work was based on the assumption that the application of the hybrid technique of regression-kriging (RK), using legacy data of soil survey and covariates of relief and remote sensing provide representative map of soil attributes of an area in Cuesta of Botucatu. The goal was to apply locally, in two depths, prediction models and continuous representation of Soil Redness Index (IAV), base saturation index (V%), sand content and clay content, cation-exchange capacity (CTC) and pH of the soils in Edgardia Experimental Farm, for which are available soil survey data. The predictor covariates were derived from an Digital Elevation Model (MDE) and an orbital image. They were all standardized at spatial resolution of 10 m, the methods were selected by checking significant linear correlation between attributes and covariates and spatial autocorrelation of attributes or residues of multiple linear regressions (RLM). The data were separated into training and validation subsets. The correlation coefficients (r) between soil attributes and covariates were significant and ranged from -0.40 to 0.51. The predictors more ... / Doutor

Banca: Ricardo Antônio Ferreira Rodrigues / Banca: Salatiér Buzetti / Banca: Paulo Leonel Libardi / Banca: Carlos Roberto Espindola / Banca: Edson Eiji Matsura / Resumo: Por se tratar de um cultivo perene, a pupunheira se torna bastante exigente quanto aos atributos físico-químicos do solo, refletindo de forma mais imediata ao efeito causado pelo manejo que é adotado. O objetivo do presente trabalho foi identificar por meio de análises lineares, geoestatísticas e multivariadas os atributos do solo que melhor expliquem a variabilidade linear e espacial da cultura da pupunha em Registro, SP. Em 2014, foi analisado os componentes produtivos da cultura da pupunha em função de alguns atributos físicos e químicos de um Cambissolo eutrófico. Foi demarcada a malha geoestatística, com 54 pontos amostrais, numa área de 1 ha. Quanto ao ponto de vista linear, diâmetro do caule e altura da planta foi explicada pelo modelo potencial direto e exponencial, o que foi altamente significativa, devido à resistência mecânica à penetração. A classificação em grupos foi feita por três métodos: método de agrupamentos hierárquico, método não-hierárquico k-means e análise de componentes principais. Assim, pode se concluir que a análise multivariada por meio da análise de componentes principais permitiu identificar dois grupos que explicam 67,8 % da variabilidade total dos dados, onde estes grupos são constituídos pelos atributos físicos densidade do solo, porosidade total, umidade gravimétrica e resistência mecânica à penetração que se mostraram com maior poder de explicação da variação do diâmetro do caule da pupunha. Já para os atributos químicos a classificação em grupos foi feita por três métodos: método de agrupamentos hierárquico, método não-hierárquico k-means e análise de componentes principais. Pode se concluir que a análise multivariada por meio da análise de componentes principais permitiu identificar dois grupos que explicam 56,62 % da variabilidade total dos dados, onde estes grupos são constituídos pelos atributos químicos matéria orgânica, potássio, magnésio.. / Abstract: Because it is a perennial crop, peach palm becomes quite demanding as physico-chemical soil properties, reflecting more immediately the effect caused by the management that is adopted. The objective of this study was to identify through linear analysis, geostatistics and multivariate soil attributes that best explain the linear and spatial variability of the culture of peach palm in Registro, SP. In 2014, was analyzed productive components of the culture of peach palm due to some physical and chemical attributes of a eutrophic Cambisol. It was demarcated the geostatistical grid with 54 sampling points in an area of 1 ha. As the linear point of view, stem diameter, and plant height was explained by the potential direct and exponential model, which was highly significant due to the mechanical resistance. Classification into groups was made by three methods: method of hierarchical clusters, non-hierarchical k-means method and principal component analysis. Thus, it can be concluded that the multivariate analysis using principal component analysis allowed us to identify two groups that explain 67.8% of the total variability of the data, where these groups are made up of the physical attributes of the soil density, total porosity, humidity and gravimetric penetration resistance that demonstrated greater explanatory power of the variation of peach palm stem diameter. As for the chemical attributes the classification in groups was made by three methods: method of hierarchical clusters, non-hierarchical k-means method and principal component analysis. It can be concluded that the multivariate analysis using principal component analysis allowed us to identify two groups that explain 56.62% of the total variability of the data, where these groups are made up of the chemical attributes organic matter, potassium, magnesium, aluminum, exchange capacity cation and aluminum saturation stood out as main components that demonstrated greater explanatory power...