Many geological deposits contain nonlinear anisotropic features such as veins, channels, folds or local changes in orientation; numerical property modeling must account for these features to be reliable and predictive. This work incorporates locally varying anisotropy into inverse distance estimation, kriging and sequential Gaussian simulation. The methodology is applicable to a range of fields including (1) mining-mineral grade modeling (2) petroleum-porosity, permeability, saturation and facies modeling (3) environmental-contaminate concentration modeling. An exhaustive vector field defines the direction and magnitude of anisotropy and must be specified prior to modeling. Techniques explored for obtaining this field include: manual; moment of inertia of local covariance maps; direct estimation and; automatic feature interpolation.
The methodology for integrating locally varying anisotropy into numerical modeling is based on modifying the distance/covariance between locations in space. Normally, the straight line path determines distance but in the presence of nonlinear features the appropriate path between locations traces along the features. These paths are calculated with the Dijkstra algorithm and may be nonlinear in the presence of locally varying anisotropy. Nonlinear paths do not ensure positive definiteness of the required system of equations when used with kriging or sequential Gaussian simulation. Classical multidimensional scaling is applied to ensure positive definiteness but is found to be computationally infeasible for large models, thus, landmark points are used for efficiency with acceptable losses in precision. The methodology is demonstrated on two data sets (1) net thickness of the McMurray formation in northern Alberta and (2) gold grade in a porphyry copper deposit. Integrating LVA into numerical modeling increases local accuracy and improves leave-one-out cross validation analysis results in both case studies. / Mining Engineering
Identifer | oai:union.ndltd.org:LACETR/oai:collectionscanada.gc.ca:AEU.10048/1107 |
Date | 06 1900 |
Creators | Boisvert, Jeff |
Contributors | Deutsch, Clayton (Civil and Environmental Engineering), Dimitrakopolous, Roussos (McGill University), Nouri, Alireza (Civil and Environmental Engineering), Schuurmans, Dale (Computer Science), Askari-Nasab, Hooman (Civil and Environmental Engineering), Apel, Derek (Civil and Environmental Engineering) |
Source Sets | Library and Archives Canada ETDs Repository / Centre d'archives des thèses électroniques de Bibliothèque et Archives Canada |
Language | English |
Detected Language | English |
Type | Thesis |
Format | 30703276 bytes, application/pdf |
Relation | Boisvert J and Deutsch C, 2010. Programs for Kriging and Sequential Gaussian Simulation with Locally Varying Anisotropy Using Non-Euclidean Distances. Computers and Geosciences. In press, accepted March 2010., Boisvert J, Manchuk J and Deutsch C, (2009). Kriging in the presence of locally varying anisotropy using non-Euclidean distances. Mathematical Geosciences 41(5): 585-601., Boisvert J and Deutsch C, 2008. Shortest anisotropic path to reproduce complex geological features. In Ortiz J and Emery X (eds.) Geostats 2008: Proceedings of the eighth international geostatistics congress, Gecamin. 1041-1046. |
Page generated in 0.0025 seconds