Метод главных компонент при оценке региональной составляющей поля

Содержание метода главных компонент.

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

Математической моделью метода служит, как и для корреляционно-регрессионного анализа, система случайных величин . При этом каждая случайная величина обычно содержит n – наблюдений, т.е. исходный массив геоданных представлен матрицей X размерностью :

В качестве числа N может быть число профилей съемки, и тогда с помощью главных компонент решается задача оценки регионального тренда (региональной составляющей). Если число N представлено совокупностью различных методов и атрибутов, то путем метода главных компонент решается задача комплексного анализа по разделению исследуемой территории на классы.

Значение n определяет число точек наблюдений по отдельным профилям съемки. Суть метода главных компонент состоит в переходе от системы случайных величин к новой системе случайных величин , ориентируясь на поведение дисперсий . При этом главная компонента определяется как линейная комбинация исходных случайных величин ,

; (3.1),

причем первая главная компонента :

(3.2)

обладает максимальной дисперсией среди всех возможных линейных комбинаций вида (3.1).

Величины являются коэффициентами перехода от одной системы случайных величин к другой системе случайных величин . Дисперсии линейных комбинаций располагаются в убывающем порядке, т.е. .

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

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

Далее для матрицы В или матрицы R находятся их собственные значения и соответствующие этим собственным значениям собственные векторы, которыми являются коэффициенты перехода .

Обычно ограничиваются вычислением первых двух-трех главных компонент, поскольку в этих компонентах сосредоточена основная энергия исходных данных. Так, первая главная компонента содержит 70-90% всей энергии (дисперсии) поля или дисперсий совокупности атрибутов при комплексном их анализе.

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

Ниже рассматриваются алгоритмы метода главных компонент при решении задач комплексного анализа геоданных и оценки региональной составляющей одного поля.

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

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

Алгоритм оценки региональной составляющей на основе метода главных компонент сводится к реализации следующих процедур:

1.Вычисление средних значений поля по каждому профилю , где n – число точек наблюдений, i =1,……. N, N – число профилей, и ковариаций данных различных пар профилей , i,j =1,…. N. Поскольку система x₁,…. x_N представлена измерениями одного и того же поля, нет необходимости рассчитывать коэффициенты корреляции, которые используются при комплексном анализе данных, измеренных в разных физических единицах.

2.Составление ковариационной матрицы исходных данных по их коэффициентам ковариаций

Матрица В симметрична относительно главной диагонали, т.е. , а по диагонали расположены дисперсии значений поля каждого профиля.

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

4.Вычисление значений собственного вектора матрицы , соответствующего максимальному собственному значению из системы линейных уравнений или

Значения собственного вектора определяются с учетом нормировки . Физический смысл этой нормировки состоит в том, чтобы преобразованные данные, т.е. значения региональной составляющей, не отличались бы по масштабу от исходных значений поля, а физический смысл значений заключается в определении весовых коэффициентов для каждого профиля.

5.Нахождение значений первой главной компоненты т.е.

Физический смысл значений первой главной компоненты состоит в том, что они определяют весовые коэффициенты для каждого пикета исходных данных, аналогично тому, как значения определяют весовые коэффициенты для каждого профиля съемки.

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

=

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

Эффективность метода главных компонент иллюстрируется на рис.3.1(стр.95), на котором приведены исходное поле силы тяжести (а), оценка региональной составляющей (б) и локальная составляющая (в). Пунктиром на рисунке показана область рудного объекта.

Рис.3.1.Исходное поле силы тяжести (а), оценка региональной составляющей (б) и локальной составляющей (в) с использованием метода главных компонент.

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

Следует отметить эффективность применения метода главных компонент при решении задач интерполяции. При этом задача интерполяции физического поля, представленного в виде функции двух переменных x и y, сводится к интерполяции функций, зависящих от одного аргумента. Однако чаще всего метод главных компонент используется при решении задач комплексного анализа данных (см.раздел 9.3).

ГЛАВА IV. КОРРЕЛЯЦИОННЫЕ ФУНКЦИИ ГЕОФИЗИЧЕСКИХ ПОЛЕЙ.

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

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

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

Применение модели случайного процесса для описания геолого-геофизических данных существенно расширяет круг решаемых при обработке и интерпретации задач.

Дата добавления: 2015-06-26