Таблица 5.11

Содержание золота в пробах, г/т

Пункт измерений x, м	Содержание золота в пробах f(x), г/т	Пункт измерений x, м	Содержание золота в пробах f(x), г/т
1	5,2	14	4,6
2	3,5	15	4,7
3	4,6	16	6,7
4	5,2	17	6,8
5	6,8	18	4,9
6	6,1	19	3,0
7	7,0	20	3,5
8	7,4	21	3,8
9	7,8	22	5,4
10	7,6	23	4,1
11	6,2	24	4,0
12	5,2	25	5,5
13	2,6	26	5,1

Пример 5.11. В горной выработке определено содержание зо­лота через 1 м (табл.5.11, рис.5.22). Требуется построить вариограмму содержаний и аппроксимировать ее какой-либо моделью.

Эмпирическая вариограмма, вычисленная по формуле (5.4), показана на рис.5.18. Она содержит 13 точек, т.е. рассчитана до 12-го шага (половина длины ряда наблюдений).

Начальная часть ва­риограммы (первые четыре точки) расположена практически на одной линии, поэтому для построения теоретической вариограммы используем линейную модель, как на рис.5.19. Уравнение теоретической вариограммы, рассчитанное по четырем точкам, γ(h) = 0,8h. Радиус автокорреляции R = 2,5 м получен при пересечении линии теоретической вариограммы с линией дисперсии исходных данных D = 2,03.

5.5.2. Влияние на вариограмму геометрической

базы измерений

Многие характеристики пространственной переменной зависят от геометрической базы измерений, т.е. от формы, размеров, а в анизотропных телах и от ориентировки области измерений. К таким характеристикам, в первую очередь, относятся дисперсия и вариограмма.

Если область, в которой производится измерение, настолько мала, что ее размерами можно пренебречь, то измерение рассматривается как точечное. Характеристики, полученные из таких измерений, рассматриваются как измерения на точечной геометрической базе.

Пусть в геологическом поле имеется множество точек, в каждой из которых пространственная переменная имеет значение f(x). Если геологическое поле разделить на области объемом v и в каждом из них найти среднее из N точечных значений, то получим новую пространственную переменную f_v(x):

при , (5.26)

где p(m) – весовая функция, зависящая от взаимного расположения точек в объеме v.

При увеличении размеров геометрической базы происходит усреднение значений пространственной переменной, соответственно, будут меняться некоторые статистические характеристики. При переходе к бесконечному множеству точек в объеме v выражение (5.26) преобразуется в интегральную форму:

при . (5.27)

Нахождение средних значений f_v(x) по точечным значениям f(x) называется регуляризацией (сглаживанием) пространственной переменной в объеме v. Меняя объем v, будем получать различные регуляризованные пространственные переменные, соответственно, будут меняться и их характеристики.

Изменение дисперсии по мере увеличения объема v описывается формулой

, (5.28)

где D – дисперсия пространственной переменной на точечной базе; N – число точек в объеме v.

Вычитаемое в формуле (5.28) представляет собою среднее значение вариограммы, получаемое при всех возможных положениях точек i и j в объеме v. При переходе к бесконечному множеству точек формула (5.28) приобретает интегральный вид:

. (5.29)

Из формул (5.28) и (5.29) следует, что по мере увеличения объема v дисперсия пространственной переменной уменьшается.

Подобные формулы существуют и для вариограммы:

(5.30)

или в интегральной форме

(5.31)

Первое слагаемое в этих формулах характеризует среднюю ковариограмму, которая получается как средняя из вариограмм между всеми точками в объеме v, а вычитаемое – в таком же объеме, смещенном на расстояние h.

Формулы (5.28)-(5.31) выведены в предположении, что объем v мал по сравнению с объемом изучаемого геологического объекта, иначе сказывается граничный эффект, который можно учесть введением специальной весовой функции – геометрической ковариограммы p(h):

, (5.32)

где k(х) = 1, когда точка х находится внутри геологического объекта, k(х) = 0, когда точка х находится за его пределами.

Геометрический смысл ковариограммы заключается в степени перекрытия двух объемов, смещенных относительно друг друга на расстояние h. С учетом функции p(h) формулы (5.29) и (5.31) приобретают следующий вид:

; (5.33)

. (5.34)

Если вариограмма аппроксимирована каким-либо алгебраическим выражением, т.е. задан вид математической модели, то интегралы во многих случаях могут быть вычислены и заменены соответствующими алгебраическими выражениями.

Пример 5.12. Имеются исходные данные по содержанию золота в горной выработке (табл.5.11). Требуется показать изменение дисперсии и вариограммы при увеличении размера проб.

Будем объединять и усреднять исходные данные по два, три и четыре соседних измерения путем сглаживания методом скользящего окна, что аналогично увеличению размеров проб исходных данных v. Количество исходных данных будет уменьшаться. Рассчитаем дисперсию D усредненных данных по формуле (2.2) и вариограмму (h) по формуле (5.4) по каждому варианту усреднения. Результаты расчета изображены на рис.5.23 и 5.24. На рисунках видно, что при увеличении размера проб дисперсия уменьшается, по сравнению с рис.5.18 убывает также и вариограмма.

5.5.3. Понятие о кригинге

На основе геостатистических моделей создан новый метод интерполяции результатов между пунктами измерений, который получил название кригинг (крайгинг) в честь автора метода Д.П.Крига [7]. Как показано Д.Матероном [12], кригинг дает минимальную сумму квадратов отклонений прогнозных значений от фактических значений пространственной переменной. Существует несколько видов кригинга. Вначале рассмотрим простейший – точечный кригинг.

Пусть имеется область v, в которой произведены измерения пространственной переменной по дискретной сети (рис.5.25). Каждое измерение считаем точечным, т.е. пренебрегаем геометрической базой измерений.

Н еобходимо рассчитать прогнозное значение пространственной переменной в точке Z. Предположим, что геологический объект изотропный, тогда вокруг точки Z можно провести окружность с радиусом влияния R. Все измерения внутри окружности влияют на значение пространственной переменной в точке Z. За пределами окружности результаты измерений в прогнозировании не участвуют.

Прогнозное значение находится в результате выполнения ряда последовательных операций:

1. Рассчитывается эмпирическая вариограмма, как на рис.5.18.

2. Осуществляется аппроксимация эмпирической вариограммы теоретической вариограммой (см. рис.5.19-5.21), т.е. выбирается геостатистическая модель.

3. В области с радиусом R рассчитываются все расстояния между пунктами измерений, а также между пунктами измерений и пунктом прогноза.

4. Составляется симметричная система уравнений кригинга следующего вида:

. (5.35)

Здесь ₁₂ – значение теоретической вариограммы между пунктами 1 и 2 (учитывается расстояние между этими пунктами); m – количество пунктов, участвующих в расчете (в данном случае m = 5).

5. Решение системы дает весовые коэффициенты p₁, p₂, p₃, …, p_m.

6. Сумма весовых коэффициентов должна быть равна единице, для чего все коэффициенты делятся на их сумму.

7. С помощью весовых коэффициентов вычисляется прогнозное значение свойства в пункте Z:

, (5.36)

где z_i – значения свойства в пунктах c 1-го по m-й.

Если понадобится вычислить прогнозное значение в следующем пункте, то операции 3-6 придется выполнить снова. Если в пределах радиуса влияния не окажется ни одного пункта измерения, то кригинг нельзя применять. Приходится обращаться к другим методам интерполяции данных, чаще всего используется метод обратных расстояний.

Иногда при составлении системы уравнений (5.35) в матрице коэффициентов оказывается много нулей (расстояния между точками больше радиуса автокорреляции). В этом случае система уравнений является несовместимой, т.е. не имеет решения (или имеет много решений). Тогда весовые коэффициенты лучше всего находить путем деления значений D – γ_1z, D – γ_2z, D – γ_3z, , D – γ_mz на дисперсию D и потом приводить сумму коэффициентов к единице.

У точечного кригинга имеется один недостаток – в окрестностях пунктов измерений касательная к прогнозным значениям ориентирована горизонтально, как в полиномиальной модели (см. рис.5.2). Чтобы устранить этот недостаток, применяется универсальный кригинг – комбинация тренда с кригингом. Вначале вычисляется тренд f(x), потом остаток от тренда и по остатку осуществляется кригинг. Прогнозное значение находится по формуле

(5.37)

В качестве тренда используется аппроксимирующий полином не выше третьего порядка. Использование тренда частично снимает явление анизотропии, что позволяет применять изотропную схему кригинга.

Точечный кригинг можно распространить на объем v. Во всех точках объема с помощью кригинга рассчитывают прогнозное значение свойства, а потом находят среднее из всех значений. Если пункты измерений не точечные, а имеют какой-то объем, то при прогнозировании учитывают все точки в этом объеме. Практически все задачи с объемными геометрическими базами решаются путем вычисления определенных интегралов по этим объемам численными методами.

Пример 5.13. Нужно провести точечный кригинг, используя схему рис.5.25 и дополнив недостающие данные.

Примем изотропную линейную модель, радиус влияния R = 14 м, дисперсию исходных данных D = 0,42. Исходные данные приведены в табл.5.12. Вариограмма для такой модели имеет следующий вид: до 14 м и после 14 м соответственно

γ(h) = D/R = 0,03h;	(5.3)
γ(h) = D = 0,420.