Уточнение метода кригинга для исследования геоповерхностей > Полезные советы
Тысяча полезных мелочей    

Уточнение метода кригинга для исследования геоповерхностей

Уточнение метода кригинга для исследования геоповерхностей

Малов А.А. Максимов И.И. Статья в формате PDF 105 KB

Одним из основных геостатистических методов исследования окружающей среды является метод кригинга, согласно которому оцениваемое значение регионализованной переменной величины ψр в точке p определяется как взвешенное среднее известных наблюдений в соседних точках по формуле:

где Wi — вес i-го значения регионализованной переменной величины ψi по отношению к оцениваемой точке p из k соседних точек.

Метод кригинга предусматривает решение системы уравнений [1]:

 

где γ(ξij) — значение полувариограммы для расстояния ξij между точками i и j; γ(ξip) — значение полувариограммы для расстояния ξip между известной точкой i и оцениваемой точкой p, λ— множитель Лагранжа.

Однако применение системы уравнений (2) нередко приводит к тому, что некоторые веса Wi оказываются либо меньше нуля, либо больше 1. В частности, такие результаты приведены Дэвисом [1]. Очевидно, это существенно искажает физический и математический смысл весов Wi и приводит к неправильной оценке величины ψp.

Подобные результаты могут быть вызваны либо неточным подбором модели полувариограммы, либо наличием среди точек наблюдений величины ψi точек, не оказывающих пpaктического влияния на оцениваемое значение ψр. Такие точки должны иметь веса, равные нулю.

Дисперсия оценки методом кригинга может быть выражена формулой Матерона [2]:

 

где - дисперсия ПЭС относительно ее среднего значения; - ковариация среднего значения и значения ; - ковариация значений и .

Система уравнений (2) выводится из условия минимума дисперсии в формуле (3), условия для весов Wi :

 

а также заменой ковариаций и соответственно на разности и .

Однако, если подставить в формулы (3) и (4) вместо Wi () и выводить систему уравнений кригинга относительно , а затем сделать обратную подстановку (Wi вместо ), то получается система уравнений кригинга

позволяющая учесть ограничения.

Так как решение системы уравнений (5) может быть затруднено из-за ее нелинейности, то для ее компьютерного решения предлагается следующий алгоритм:

Решить систему линейных уравнений (2).

Если полученное решение не отвечает ограничениям , то исключить из рассмотрения точку наблюдения с наименьшим отрицательным весом (или присвоить ему значение ноль), иначе - закончить вычисления.

Перейти после указанных изменений к пункту 1.

Полученное при реализации данного алгоритма (или использования формулы (5)) решение позволяет учесть ограничения и точнее оценить значение регионализованной переменной величины ψр.

Литература

  1. Дэвис Дж. С. Статистический анализ данных в геологии. Пер. с англ. В 2 кн. Пер. В.А. Гoлyбевой; Под ред. Д.А. Родионова. Кн. 1. — М.: Недра, 1990. — 319 с; Кн. 2. — М.: Недра, 1990. — 427 с.
  2. Матерон Ж. Основы прикладной геостатисти- ки. — М.: Мир, 1968. — 408с.


ОПЫТ ПРОГНОЗИРОВАНИЯ ДИНАМИКИ ЦЕН НА ЖИЛЬЕ В ИВАНОВСКОЙ ОБЛАСТИ

ОПЫТ ПРОГНОЗИРОВАНИЯ ДИНАМИКИ ЦЕН НА ЖИЛЬЕ В ИВАНОВСКОЙ ОБЛАСТИ На основе построения тренд-сезонных моделей исследуется динамика цен на первичном и вторичном рынках жилья Ивановской области в период 2000-2007 гг. В статье освещаются основные этапы построения моделей, приводятся количественные оценки их параметров. Особое внимание уделяется присутствию S – образной кривой роста в динамике цен на жилье. В результате использования методики с учетом индексов сезонности получены средние прогнозные значения цен на жилье Ивановской области. ...

27 04 2026 8:53:44

О ЮРИДИЧЕСКИХ ЛИЦАХ

О ЮРИДИЧЕСКИХ ЛИЦАХ Статья в формате PDF 152 KB...

21 04 2026 9:59:48

Сравнительный анализ динамики формирования алкогольной зависимости у крыс линий Вистар и WAG/Rij с генотипом А1/А1 по локусу Taq 1A дофаминового рецептора второго типа (DRD2)

Сравнительный анализ динамики формирования алкогольной зависимости у крыс линий Вистар и WAG/Rij с генотипом А1/А1 по локусу Taq 1A дофаминового рецептора второго типа (DRD2) Изложенные в статье результаты свидетельствуют о тождестве динамики формирования алкогольной зависимости у неэпилептической линии крыс (Вистар) и крыс с абсансной эпилепсией (WAG/Rij). ...

17 04 2026 16:59:35

Правовые аспекты эвтаназии

Правовые аспекты эвтаназии Статья в формате PDF 102 KB...

08 04 2026 18:44:36

ВОЗМОЖНОСТИ МОДЕЛИРОВАНИЯ ЭКВАТОРИАЛЬНОЙ ЭЛЕКТРОСТРУИ И ЕЁ СЛЕДСТВИЙ

ВОЗМОЖНОСТИ МОДЕЛИРОВАНИЯ ЭКВАТОРИАЛЬНОЙ ЭЛЕКТРОСТРУИ И ЕЁ СЛЕДСТВИЙ Для исследования вариаций параметров живых существ, обитающих в биосфере в разных широтных регионах, в частности экваториальных, построена модель экваториального электроджета, основанная на численном решении дифференциальных уравнений второй степени для потенциала, вызванного прострaнcтвенным зарядом. ...

28 03 2026 23:13:30

Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::