- •1.1. Геологические объекты и их свойства
- •Характеристика текстуры руды
- •1.2. Понятие о математическом моделировании
- •2.1. Одномерная статистическая модель
- •Содержание никеля в руде х
- •Расчет начальных моментов случайной величины
- •Расчет центральных моментов случайной величины
- •Частота и частость содержания железа в руде
- •Расчет статистических характеристик по сгруппированным данным (по данным гистограммы рис.2.1)
- •2.2. Законы распределения случайных величин
- •Функции нормального распределения
- •Функция распределения Стьюдента Fk (t) в зависимости от числа степеней свободы k
- •Плотность вероятности распределения Стьюдента fk(t) в зависимости от числа степеней свободы k
- •Коэффициенты вероятности t распределения Стьюдента при заданной вероятности и степени свободы k
- •Коэффициенты вероятности распределения 2 при заданных вероятности и числе степеней свободы k
- •Расчет графика плотности вероятности нормального закона
- •Подбор функции преобразования случайной величины х (содержание TiO2 в магнетите)
- •2.3. Геологические приложения
- •Пример выявления аномальных значений
- •Результат вычисления истинных характеристик
- •3.1. Двухмерная статистическая модель
- •Расчет коэффициента корреляции
- •Расчет параболической зависимости, аппроксимирующей изменение мощности рудного тела
- •Сравнение фактической и расчетной (теоретической) мощности
- •Дисперсии отклонений
- •3.2. Геологические приложения двухмерной статистической модели
- •Расчет случайной погрешности химического анализа
- •Расчет систематической погрешности анализов серы
- •Результаты опытов по измерению радиоактивности препарата
- •Содержание железа в руде и плотность руды
- •Расчет частот n
- •4.1. Многомерная статистическая модель
- •Расчет статистических характеристик многофакторной зависимости
- •Сравнение фактических и расчетных содержаний золота
- •Химические анализы проб руды
- •Матрица коэффициентов корреляции
- •Состав магнетита, %
- •Преобразованная матрица коэффициентов корреляции
- •Матрица коэффициентов корреляции
- •Расчет первого собственного числа
- •Расчет матрицы коэффициентов первого фактора
- •Первая остаточная матрица коэффициентов корреляции
- •Факторные нагрузки
- •Нормированный состав магнетита
- •Главные компоненты z
- •Средний состав групп магнетита, %
- •Состав циркона, %
- •Сложные признаки
- •Результаты распознавания объектов
- •Результаты измерения мощности рудного тела
- •Данные по скважинам
- •5.4. Вероятностные модели
- •Результаты измерения мощности рудного тела
- •Значения пространственной переменной
- •Содержание молибдена, %
- •Сглаживание исходных данных
- •Результаты измерений пространственной переменной и некоторые вычисленные данные
- •Результаты расчета всех гармоник
- •Расчет периодической изменчивости
- •Дисперсия различных длин волн
- •5.5. Основы геостатистики
- •Содержание золота в пробах, г/т
- •Данные для кригинга, м
- •Значения вариограммы (h)
- •Реестр разведочных скважин
- •Замеры искривлений скважин
- •Банк данных геологической документации
- •Банк опробования рядовых проб
- •Координаты пунктов искривлений
- •Банк рудных пересечений
- •Банк геологических границ
- •Глава 1. Общие сведения 7
- •Глава 2. Одномерная статистическая модель и ее применение в геологии 23
- •Глава 3. Двухмерная статистическая модель и ее применение в геологии 78
- •Глава 4. Многомерная статистическая модель и ее применение в геологии 115
- •Глава 5. Математическое моделирование пространственных геологических закономерностей 159
- •Глава 6. Основы математического моделирования месторождений 203
5.5. Основы геостатистики
5.5.1. Вариограмма и ее аппроксимации
В основе геостатистической группы математических моделей лежит гипотеза о том, что случайный результат измерений обусловлен случайным расположением сети наблюдений. При перемещении сети наблюдений результаты измерений будут другие, но сохраняется одна характеристика – средний квадрат разности между результатами измерений на расстоянии h. Это возможно при эргодичном характере пространственной переменной. На основе гипотезы введена вариограмма γ(h) – главная характеристика в геостатистике. Она равна полусумме среднего квадрата разности между результатами измерений при шаге h и выражается формулой (5.4).
Следует отметить, что вариограмма тесно связана со случайными функциями. Сумма вариограммы и ковариации (автокорреляционной функции) равна дисперсии исходных данных:
γ(h) + K(h) = D. (5.25)
График вариограммы зависит от характера дискретности пространственной переменной. Для непрерывных пространственных переменных (например, мощности пласта) вариограмма начинается с нулевой отметки и возрастает до дисперсии исходных данных (рис.5.16). Для дискретных пространственных переменных вариограмма начинается с некоторой величины С (рис.5.17), называемой эффектом самородков, потому что он был вначале установлен на месторождениях золота.
Вариограмма имеет радиус влияния R, который идентичен радиусу автокорреляции в случайной функции. За пределами радиуса влияния вариограмма постоянная и равна дисперсии D, т.е. ее влияние отсутствует. Как упоминалось, радиус R является векторной величиной. В изотропных геологических телах по всем направлениям радиус влияния описывает окружность, за пределами которой влияние вариограммы отсутствует. В анизотропных геологических телах значения радиуса R по разным направлениям столь различны, что напоминают восьмерку.
В пункте с координатой R вариограмма может иметь горизонтальную касательную и плавно переходить в линию дисперсии.
Эмпирическая вариограмма, получаемая на основе дискретных измерений, строится по отдельным точкам и имеет вид ломаной линии (рис.5.18). Чтобы использовать ее для дальнейших вычислений, необходимо выполнить аппроксимацию вариограммы какой-либо теоретической кривой. Вид аппроксимирующей функции определяет вид геостатистической модели. Существует много видов аппроксимирующих функций. Наибольшее распространение получили четыре функции и, соответственно, четыре геостатистические модели.
Первая модель линейная. В ней вариограмма аппроксимируется прямой линией (рис.5.19), что нередко близко к действительности. Во второй модели вариограмма аппроксимируется кубическим полиномом, который не имеет общей касательной с линией дисперсии (рис.5.20). Эта модель лучше всего соответствует действительности. Третья модель (сферическая) также аппроксимируется кубическим полиномом, но он имеет горизонтальную касательную при соприкосновении с линией дисперсии (см. рис.5.16 или рис.5.17). Четвертая модель (модель Де-Вийса) характеризуется тем, что шаг вариограммы откладывается в логарифмическом масштабе, а сама вариограмма представляется в виде прямой линии (рис.5.21).
Выбор варианта аппроксимации носит произвольный характер. Наилучшим вариантом можно считать тот, который дает наименьшую дисперсию отклонений эмпирических значений от теоретических.
