Геоинформатика и геостатистика: обзор

Применение геостатистических методов для определения аномальной зоны месторождения медной руды в районе рудника Польковице (регион Любин-Серошовице, юго-западная часть Польши)

Барбара Намысловска-Вильчинска

В данной статье представлены результаты исследований пространственного изменения основных геологических параметров стратегических минеральных ресурсов, а именно месторождений медной руды рудника Польковице, расположенных в регионе Любин-Серошовице (юго-западная часть Польши). Исследования были основаны на данных, полученных путем опробования месторождения руды Cu с использованием бороздковых образцов, распределенных в основном равномерно (с интервалом 15-20 м) по площади P-1 рудника Польковице. Исследования касались содержания Cu, мощности и количества (накопления) (извлекаемого) месторождения, сосредоточенного в песчаниках Вайтлигендес, медных сланцах Цехштейна и известково-доломитовых формациях. Для определения аномальных зон месторождения руды Cu использовались выбранные геостатистические методы, а именно обычный кригинг, связанный индикаторный кригинг и моделирование условных поворотных полос.

На первом этапе пространственных исследований использовался ординарный (блочный) кригинг для оценки оценочных средних значений Z* (с минимальной дисперсией оценки σk 2) параметров месторождений в пределах блока P-1. Изотропные вариограммы были рассчитаны и смоделированы с использованием геостатистических моделей, после чего был применен ординарный кригинг. Далее были рассчитаны индикаторные полувариограммы для следующих пороговых значений изучаемых параметров: содержание Cu>1,7%, 2,2%, 2,7%; мощность>2,5
м, 3,0 м, 3,5 м и количество>20[%mt/m2],>25 [%mt/m2],>30 [%mt/m2]. Значения вероятностей P в центрах блоков элементарной сетки, охватывающей добычной блок P-1, были оценены с помощью связанного индикаторного кригинга с учетом параметров моделей индикаторных полувариограмм. В результате проведенных расчетов показателей получены картины распределения вероятностей значений P на растровых картах превышения предполагаемых пороговых значений.

На последующем этапе анализа была применена методика моделирования условных поворотных полос. Моделирование поворотных полос было выполнено с учетом расчетов теоретических моделей, подобранных к гауссовым полувариограммам параметров месторождения. Разработаны статистические карты распределений моделируемых значений Zs для количества (среднее моделируемых значений Zs, стандартное отклонение реализаций σs, наибольшая реализация Zs, наименьшая реализация Zs) и растровые карты вероятностей P превышения предполагаемых пороговых значений (выше пороговых значений количества > 20[%mt/m2], >25[%mt/m2], > 30[%mt/m2] и ниже пороговых значений < 20[%mt/m2], <25[%mt/m2], < 30[%mt/m2].
Сравнение полученных результатов относительно используемых геостатистических методов позволяет сделать лучший выбор, указав более полезный и эффективный метод оценки параметров минеральных ресурсов и определения аномальной зоны месторождения медной руды.