
- •Москва 2008 Аннотация
- •Оглавление.
- •Глава IX. Обработка комплексных геофизических наблюдений. 102
- •Введение
- •Глава I. Оценка статистических и градиентных характеристик геофизических полей.
- •1.1.Статистические характеристики геополей.
- •1.2.Градиентные характеристики геополей.
- •Глава II. Корреляционно-регрессионный анализ и его применение.
- •2.1. Корреляция и регрессия.
- •2.2. Линейная регрессия и ее применение.
- •2.3. Нелинейная регрессия и ее применение.
- •2.3. Множественная регрессия и ее применение.
- •Глава III. Метод главных компонент и его применение
- •3.1. Содержание метода главных компонент.
- •3.2. Метод главных компонент при оценке региональной составляющей поля.
- •Глава IV. Корреляционные функции геофизических полей.
- •4.1. Автокорреляционная функция и ее применение.
- •4.2. Взаимно корреляционная функция и ее применение.
- •4.3. Двумерные корреляционные функции.
- •Глава V. Спектральный анализ геофизических сигналов.
- •5.1. Спектры дискретно заданного сигнала.
- •5.2. Спектры непрерывных сигналов.
- •5.3. Спектры стационарного случайного процесса.
- •5.4. Свойства преобразований Фурье.
- •5.5. Применение спектрального анализа при обработке геофизических данных.
- •5.6. Вейвлет-анализ геофизических полей.
- •Глава VI. Линейная фильтрация геофизических полей.
- •6.1. Понятие о линейном фильтре. Оператор свертки.
- •6.2. Двумерные линейные фильтры.
- •6.3.Адаптивные фильтры.
- •6.3.1.Одномерная адаптивная фильтрация.
- •6.3.2.Двухмерная адаптивная фильтрация.
- •6.3.3.Полиномиальная адаптивная фильтрация.
- •Глава VII. Оптимальные линейные фильтры.
- •7.1. Фильтр Колмогорова-Винера и его применение.
- •7.1.1. Выделение аномалий антиклинального типа в гравиразведке.
- •7.1.2. Построение обратного фильтра.
- •7.1.3. Построение фильтра прогнозирования.
- •7.2. Согласованный фильтр.
- •7.3. Энергетический фильтр.
- •Глава VIII. Теория статистических решений в задачах обнаружения слабых сигналов.
- •8.1. Основные понятия теории статистических решений.
- •8.2. Критерии принятия статистических решений.
- •8.3.Коэффициент правдоподобия и вероятности ошибок I и II рода.
- •8.4. Надежность обнаружения сигнала и ее применение.
- •8.4. Способ обратных вероятностей.
- •8.5. Способ межпрофильной корреляции.
- •8.6. Способ самонастраивающейся фильтрации.
- •8.7. Многомерные аналоги метода обратных вероятностей и самонастраивающейся фильтрации.
- •8.7.1.Многомерный способ обратных вероятностей.
- •8.7.2.Многомерный способ самонастраивающейся фильтрации.
- •Глава IX. Обработка комплексных геофизических наблюдений.
- •9.1.Методы классификации многомерных наблюдений.
- •9.1.1.Метод динамических сгущений (k-средних).
- •9.1.2.Методом общего расстояния.
- •9.1.3.Метод разделения многомерных, нормальных смесей (По Петрову а.В.).
- •9.2.Методы распознавания образов в задачах интерпретации геофизических наблюдений.
- •9.3. Метод главных компонент при расчете комплексного параметра.
- •Заключение.
- •Литература.
Глава VIII. Теория статистических решений в задачах обнаружения слабых сигналов.
Развитие теории статистических решений для обработки данных обусловлено следующими причинами. Критерии оптимальной фильтрации и соответствующие им алгоритмы обработки исходных данных обеспечивают выделение сигналов на фоне помех, однако принятие решения о наличии или отсутствии полезного сигнала после фильтрации осуществляется интерпретатором визуально. Это обстоятельство делает последний этап обработки весьма субъективным, особенно для случая обнаружения слабых сигналов. Под слабым сигналом понимается сигнал, соизмеримый по интенсивности с уровнем помех и ниже этого уровня. Визуальное обнаружение таких сигналов практически исключено. В то же время проблема обнаружения слабых сигналов в разведочной геофизике приобретает все большее значение в связи с поисками объектов, залегающих на больших глубинах, слабоконтрастных по физическим свойствам объектов, малоразмерных залежей углеводородов, а также объектов, аномальные эффекты от которых осложнены интенсивными помехами самой разной природы.
Получив в результате, например, согласованной фильтрации несколько положительных экстремумов, решение о наличии сигнала принимается лишь по максимальному из них. Относительно более слабых экстремумов решение о наличии сигнала принять трудно, поскольку в оптимальных фильтрах отсутствуют пороги принятия решения о наличии сигнала.
Вторая причина, по которой теория статистических решений получила развитие, связана с отсутствием общей теоремы, утверждающей, что максимальное отношение сигнал/помеха на выходе фильтра приводит к максимальному извлечению полезной информации. Под максимальным извлечением полезной информации в настоящее время принято понимать получение апостериорных вероятностей, в частности, для задачи обнаружения сигнала - апостериорных вероятностей наличия или отсутствия сигнала.
В теории статистических решений для решения задачи обнаружения используются те или иные критерии принятия статистических гипотез о наличии (отсутствии) сигнала, которые позволяют установить порог принятия решения о наличии сигнала, а также рассчитать апостериорные вероятности наличия сигнала. Оптимальная согласованная фильтрация при этом входит в качестве составной части алгоритмов теории статистических решений в задачах обнаружения слабых сигналов.
8.1. Основные понятия теории статистических решений.
Статистическое решение связано с понятием статистической гипотезы. При решении задачи обнаружения сигнала имеют место две статистические гипотезы:
- гипотеза
о наличии сигнала,
для которой наблюденные данные
представляют сумму сигнала и помех
.
Эту гипотезу называют ненулевой и
обозначают как Н1,
т.е.
;
- гипотеза
об отсутствии сигнала,
для которой наблюденные данные
представлены только помехой, т.е.
.
Это - нулевая гипотеза
Для каждой гипотезы по значениям
можно построить статистический аналог
плотности распределения - гистограмму.
Такие плотности распределений значений
называются функциями
правдоподобия
и обозначаются как
,
где
.
Функции правдоподобия для задачи
обнаружения сигнала приведены на
рис.8.1.
Рис.8.1. Графики функций правдоподобия для гипотез Н0 и Н1. |






Принятие статистического решения о наличии сигнала не может быть безошибочным из-за осложнения сигнала помехами. При этом возникают ошибки двоякого рода.
Ошибка первого рода заключается в том, что принимается решение о наличии сигнала, т.е. о гипотезе Н1, когда на самом деле сигнал отсутствует, т.е. выполняется гипотеза Н0. Эта ошибка называется ошибкой обнаружения ложного сигнала. Она появляется тогда, когда наложение помех создает ситуацию, похожую на полезный сигнал.
Ошибка второго рода заключается в том, что принимается решение об отсутствии сигнала, т.е. о гипотезе Н0, когда на самом деле сигнал имеется, т.е. выполняется гипотеза Н1. Эта ошибка называется ошибкой пропуска сигнала. Она появляется тогда, когда помехи настолько искажают сигнал, что создается ситуация, отвечающая только помехе.
С ошибками первого и второго рода связаны соответствующие им вероятности.
Вероятность ошибки I рода или вероятность ошибки обнаружения ложного сигнала равна
(8.1)
Вероятность ошибки II рода или вероятность пропуска сигнала равна
(8.2)
Вероятности
и
соответствуют площадям, заштрихованным
на рис. 8.1.
Если для гипотез Н0 и Н1 ввести их априорные вероятности р0 и р1, то общая безусловная вероятность ошибки, связанная как с обнаружением ложного сигнала, так и с пропуском сигнала, будет равна
(8.3)
Величина
называется вероятностью
правильного обнаружения сигнала
при его наличии или надежного
обнаружения сигнала.
Соответственно вероятность правильного
необнаружения сигнала при его отсутствии
равна
.
При выборе порога
h
потери, обусловленные ошибками I
и II
рода, могут быть разными с точки зрения
экономики. Если мы обнаруживаем ложный
сигнал, то должны провести затраты на
его проверку. В случае пропуска сигнала
можно пропустить тот или иной важный
объект поисков. Поэтому вводят цены для
ошибок I
и II
рода -
и
.
Произведение
называют риском
(или потерей),
соответствующим гипотезе Н1.
Средний
риск или
функция потерь при принятии решения
определяется выражением:
(8.4),
в котором величина среднего риска зависит от значения порога h.