Одним из основных геостатистических методов исследования окружающей среды является метод кригинга, согласно которому оцениваемое значение регионализованной переменной величины ψр в точке p определяется как взвешенное среднее известных наблюдений в соседних точках по формуле:
где Wi — вес i-го значения регионализованной переменной величины ψi по отношению к оцениваемой точке p из k соседних точек.
Метод кригинга предусматривает решение системы уравнений [1]:
где γ(ξij) — значение полувариограммы для расстояния ξij между точками i и j; γ(ξip) — значение полувариограммы для расстояния ξip между известной точкой i и оцениваемой точкой p, λ— множитель Лагранжа.
Однако применение системы уравнений (2) нередко приводит к тому, что некоторые веса Wi оказываются либо меньше нуля, либо больше 1. В частности, такие результаты приведены Дэвисом [1]. Очевидно, это существенно искажает физический и математический смысл весов Wi и приводит к неправильной оценке величины ψp.
Подобные результаты могут быть вызваны либо неточным подбором модели полувариограммы, либо наличием среди точек наблюдений величины ψi точек, не оказывающих практического влияния на оцениваемое значение ψр. Такие точки должны иметь веса, равные нулю.
Дисперсия оценки методом кригинга может быть выражена формулой Матерона [2]:
где - дисперсия ПЭС относительно ее среднего значения; - ковариация среднего значения и значения ; - ковариация значений и .
Система уравнений (2) выводится из условия минимума дисперсии в формуле (3), условия для весов Wi :
а также заменой ковариаций и соответственно на разности и .
Однако, если подставить в формулы (3) и (4) вместо Wi () и выводить систему уравнений кригинга относительно , а затем сделать обратную подстановку (Wi вместо ), то получается система уравнений кригинга
позволяющая учесть ограничения.
Так как решение системы уравнений (5) может быть затруднено из-за ее нелинейности, то для ее компьютерного решения предлагается следующий алгоритм:
Решить систему линейных уравнений (2).
Если полученное решение не отвечает ограничениям , то исключить из рассмотрения точку наблюдения с наименьшим отрицательным весом (или присвоить ему значение ноль), иначе - закончить вычисления.
Перейти после указанных изменений к пункту 1.
Полученное при реализации данного алгоритма (или использования формулы (5)) решение позволяет учесть ограничения и точнее оценить значение регионализованной переменной величины ψр.
Литература
- Дэвис Дж. С. Статистический анализ данных в геологии. Пер. с англ. В 2 кн. Пер. В.А. Голубевой; Под ред. Д.А. Родионова. Кн. 1. — М.: Недра, 1990. — 319 с; Кн. 2. — М.: Недра, 1990. — 427 с.
- Матерон Ж. Основы прикладной геостатисти- ки. — М.: Мир, 1968. — 408с.
Библиографическая ссылка
Малов А.А., Максимов И.И. Уточнение метода кригинга для исследования геоповерхностей // Успехи современного естествознания. – 2003. – № 12. – С. 97-97;URL: https://natural-sciences.ru/ru/article/view?id=15593 (дата обращения: 17.08.2024).