- •«Основы компьютерных технологий решения геологических задач»
- •Введение
- •Программный пакет surfer
- •Основные возможности пакета Surfer
- •1.2. Подходы к построению карт
- •2. Методы построения сеточных моделей
- •2.1. Метод Крайгинга
- •Метод радиальных базисных функций
- •2.3. Метод обратных расстояний
- •2.4. Метод Шепарда
- •2.5. Метод минимальной кривизны
- •2.6. Метод полиномиальной регрессии
- •2.7. Триангуляция с линейной интерполяцией
- •2.8. Метод ближайшего соседа
- •2.9. Метод естественного соседа
- •2. 10. Сравнительная характеристика основных методов построения сеточной функции
- •3. Инструкция к пакету surfer
- •3.1. Режимы пакета Surfer
- •3.2. Объекты пакета Surfer
- •3.3. Команды пакета Surfer
- •3.4. Пример построения карты
- •4. Задания к лаборторным работам
- •5. Англо-русский словарь основных терминов
- •Подписано к печати Бум. Писч. № 1
- •625000, Г. Тюмень, ул. Володарского, 38
- •625039, Г. Тюмень, ул. Киевская, 52
Метод радиальных базисных функций
В этом методе искомая функция находится как линейная комбинация набора радиальных базисных функций:
, (2.3)
где а -константа, i - индекс точки измерений, i - неизвестные коэффициенты, Ri(x,y) - базисные функции, зависящие от расстояния точки (х, у) до i-ой точки наблюдения.
Существуют несколько типов базисных функций:
Inverse Multiquadric
Multilog
Мультиквадратичная (Multiquadric) , наиболее часто используется;
Natural Cubic Spline
Thin Plate Spline ,
где R2 – фактор сглаживания, чем больше будет параметр, тем более сглаженные будут контура. Разумные значения показателя находятся в интервале от среднего межточечного расстояния выборки до половины этого среднего значения.
2.3. Метод обратных расстояний
Фактически большинство интерполяционных методов используют формулу (2.1) для расчета значений функций в произвольных точках. Различаются принципы, используемые для задания весовых коэффициентов wi и, соответственно, выражения для них расчета. В методе обратных расстояний используется обратно пропорциональная зависимость весовых коэффициентов от некоторой степени расстояния между расчетной точкой и точкой наблюдения (di):
, (2.4)
здесь к - число, обычно принимаемое равным 1, 2 или 3.
Такой подход имеет вполне очевидный смысл - более удаленные точки в меньшей степени определяют значение в расчетной точке и наоборот.
2.4. Метод Шепарда
Метод Шепарда аналогичен методу обратных расстояний. Он также использует обратные расстояния при вычислении весовых коэффициентов. Интерполяция же осуществляется по формуле (2.3) при к=2. Отличие состоит в том, что при построении интерполяционной функции в локальных областях используется метод наименьших квадратов. Это уменьшает вероятность появления на сгенерированной поверхности ложных структур вокруг точек наблюдений.
2.5. Метод минимальной кривизны
Поверхность, построенная с помощью этого метода, аналогична тонкой упругой пластине, проходящей через все экспериментальные точки данных с минимальным числом изгибов. Он генерирует наиболее гладкую поверхность, которая проходит настолько близко к экспериментальным точкам, насколько это возможно, но экспериментальные точки не обязательно принадлежать приближающей поверхности. Окно Parameters позволяет задать критерий сходимости.
Параметр Max Residuals (Максимальная невязка) измеряется в тех же единицах, что и экспериментальные данные. Как показывает опыт, наиболее подходящее значение этого параметра равно, примерно, 10% от точности исходных данных. Если исходные экспериментальные данные измерены с точностью 1.0 единиц измерения, то значение Max Residuals рекомендуется положить равным 0.1. Итерации продолжаются до тех пор, пока максимальная невязка для всей итерации не станет меньше значения параметра Max Residuals.
Параметр Max Iterations (Максимальное число итераций) разумно выбрать в интервале от N до 2N, где N - число узлов генерируемой сети.