Geostatistics in the Presence of Extreme Values with Network Models of Regionalization

• Author / Creator

  Harding, Benjamin E

• Geostatistical models are often generated with widely spaced data configurations. Data collection costs prohibit exhaustive sampling and necessitate statistical inference from limited samples. Spatial prediction with sparse data in the presence of extreme values is an enduring challenge in the mining industry. Extreme values may have significant local influence, leading to overstated resources and the risk of production shortfalls. Practitioners are presented with difficult decisions for restricting extreme value influence and characterizing their spatial continuity. Inputs to numerical geologic models consist of observed data, a representative histogram and spatial controls on mineralization. Each of these components presents challenges in the presence of extreme values. Extreme values
  are rare events, making inferences about their probability of occurrence difficult. The influence of extreme values is often controlled in practice through grade capping, which could significantly impact the final resource. Extreme values’ spatial continuity often differs from the barren or lower grade background. Traditional estimation and simulation methodologies are limited in adapting to extreme values and asymmetric spatial continuity features. These challenges motivate the development of a framework for the simulation of continuous variables with explicit consideration of high-order extreme value features. The proposed network model of regionalization (NMR) framework constructs a continuous regionalized variable as a non-linear mixture of latent Gaussian factors and does not require capping or modification of extreme grade values. The network parameters are inferred
  via optimization, considering two- and multi-point connectivity features at grade thresholds. This permits the reproduction of high-order connectivity features and asymmetric spatial continuity of high and low grades that cannot be captured by a single Gaussian random function (RF) model. The latent Gaussian factors are imputed such that they exactly reproduce the observed data values when mixed. The applicability of the proposed methodology is demonstrated on a mineral deposit where the project operators note non-Gaussian, extreme value features in drillhole data. In this deposit the NMR approach shows a 7% improvement in expected metal relative to traditional approaches using a hold-out data set for validation.

• Subjects / Keywords

  • Geostatistics
  • Spatial statistics
  • Extreme values
  • Outliers

• Graduation date

  Fall 2024

• Type of Item

  Thesis

• Degree

  Doctor of Philosophy

• DOI

  https://doi.org/10.7939/r3-211b-yt62

• License

  This thesis is made available by the University of Alberta Library with permission of the copyright owner solely for non-commercial purposes. This thesis, or any portion thereof, may not otherwise be copied or reproduced without the written consent of the copyright owner, except to the extent permitted by Canadian copyright law.