CC BY
91
УДК 519.25 (550.831.05)
ОЦЕНКА ФУНКЦИИ АВТОКОРРЕЛЯЦИИ В ВИДЕ СПЛАЙНА ПЕРВОГО ПОРЯДКА
В.Н. Устинова, И.Г. Устинова
Томский политехнический университет E-mail: [email protected]
Предложен метод выделения тренда автокорреляционной функции в виде сплайна первого порядка. На примере модельной сейсмической трассы и функции Бесселя показано, что предложенный метод дает лучшее приближение автокорреляционной функции, чем приближение, полученное с использованием стандартной формулы для оценивания функции автокорреляции.
Ключевые слова:
Автокорреляционная функция, тренд автокорреляционной функции, сплайн первого порядка, сейсмические данные, обработка сейсмической информации.
Key words:
Autocorrelation function, trend of the autocorrelation function, first-order spline, seismic data, processing of seismic data.
В рамках эргодичной и стационарной модели геофизического поля важнейшие сведения о свойствах сейсмического сигнала можно получить по автокорреляционной функции Д[т], энергетическому спектру 5[ю], математическому ожиданию М[х] [1]. Наиболее информативными для оценки свойств составляющих сейсмического сигнала и полей сейсмических параметров являются функция автокорреляции и энергетический спектр [2, 3]. Оценка функции автокорреляции определяется, согласно [4], выражением
R[z] =-
1
-Н
n-|т
|У(у - у)( Уі+т- УІ
(1)
где уьу2,...,уп - текущие значения геофизического по-
1 п
ля В в і-й точке сейсмической трассы; у = — V у,
___ п ,•=—
т- разность аргументов: т=—;, і, у'=1,п Анализ формы автокорреляционной функции (АКФ) нашел применение при решении следующих задач обработки сейсмических данных:
1. Оценка корреляционных свойств сигналов и помех.
2. Расчеты весовых функций и частотных характеристик оптимальных фильтров базируются на знании АКФ сигналов и помех [3].
3. Оценка разрешающей способности сейсмической записи [4].
Часто для оценки функции автокорреляции приходится использовать короткие отрезки сейсмических записей, поэтому вопрос о точности оценивания автокорреляционной функции имеет большое практическое значение. В данной статье выделен тренд функции автокорреляции в виде сплайна первого порядка.
Постановка задачи
Функция автокорреляции Л[т] стационарного случайного процесса является одной из его важнейших характеристик второго порядка [5]. Приведем описание математической модели процесса, которая будет исследована далее.
Пусть у(/) - гауссовский стационарный случайный процесс с математическим ожиданием М[у(0]=0 и функцией автокорреляции Я[т]=Щ—]=М[у(§у(1+т)]. Пусть изучаемый процесс наблюдается на отрезке времени [0; 7], и моменты времени, в которые производятся измерения /;, известны точно. Результатами измерений являются величины у;=у(/;), = 1,И. По результатам измерений {у, 1=1, N требуется построить оценку Дт] функции корреляции Д[т] на интервале значений [0; 7] аргумента т.
Решение задачи
Пусть на отрезке значений аргумента [(к-1)т0;кт0] функция автокорреляции задана выражением кт0 -т „ т - (к - 1)г0
RM = Rk
R
k = 1,2,..., (2)
к-1 к
То Т0
т. е. функция автокорреляции имеет вид сплайна первого порядка. В этом случае решение задачи об аппроксимации функции автокорреляции состоит в нахождении оценок параметров Ло,Ль...,Д„
^ 1 М
параметр Д можно определить как Я о = — ^ у).
N /=1
Заметим, что последнее выражение является несмещенной оценкой Д0. Для нахождения оценки Як рассмотрим выражение
Sk = У У-У/
i, jE Mk
A
^o - (t/ - ti)
+B
(t/ - t,) - (k - ^o
где А, В - неизвестные параметры, которые подлежат определению, а г, у - номера моментов измерений, Мк - множество индексов ( г, у), удовлетворяющих условию Мк={( г, у):(к-1)т0<—<кт0}. Усредняя Зк по реализациям процесса у (/), получим
M[Sk] = У R[tj -1,\
i.j^Mk
A
т - (t/ - ti)
+B (tj - ti) - (k - 1)ro
i=1
или
Л. 10
М[Як] = | Я[т]
(к-1)то
Акт0-т + в т- (к-1К
й?Т.
Подставляя в последнее выражение (2) и вычисляя полученный интеграл, будем иметь
М[Як] = %2А + В)Як-1 + (а + 2В)Я ].
6
Оценку Д параметра ^ найдем из уравнения
Як =Т0[(2А + В)Як-1 + (А + 2В)Я ],
6
из которого получим
X у/у
уеМк
откуда
, кт0 - (0 - К г, - (< ) - (к - 1)Т0
- + В-
= _А[(2 А + В) Як-1 + (А + 2 В) Як ], 6
9 2А + В - 6 1
Як =------------------Як-1 +----------------х
т0 А + 2В
х X уу
/ ./'еМк
А
А + 2В
кт0 - 0/ - ^ К „0/ - */) - (к - 1)Т0
- + В-
ь0 ь0
Заметим, что полученное уравнение является разностным уравнением первого порядка с началь-
9 1 "
ным условием Я0 =—Xу2 [6]. Для того чтобы ре-
N /=1
шение уравнения в конечных разностях было устойчивым, необходимо, согласно [7], чтобы
2 А + В
< 1,
(3)
В
/ А=В
'' А=-В
А=-2В '
Рис. 1. Решение неравенства (3) графическим методом Рассмотрим A и B такие, чтобы 2А +В = 0, на-
А + 2В
пример, А = -—, B=1, тогда
9 2 х-1 3(0 ^) (3к 2)т0 (4)
9 = — X УУ
Т0 /,/еМ,; Т0
причем ^[.1^]=^, то есть полученная оценка является несмещенной.
Имитационное моделирование полученных оценок
Все расчеты и построение графиков выполнены в системе МаШсаё 2000 [8].
В качестве модели сейсмического сигнала возьмем функцию Бесселя первого рода нулевого порядка, реализация которой изображена на рис. 2, а, для которой построим функцию автокорреляции, используя выражение (1) (рис. 2, б) и выражение (4) в виде сплайна первого порядка (рис. 2, в).
Согласно [9] в приложениях при анализе временных последовательностей часто применяются стационарные случайные процессы с показательной функцией автокорреляции
Я[т] = еа\ -да<т<+да, а> 0. (5)
Другим часто встречающимся типом АКФ стационарной случайной функции является показательно-косинусная автокорреляционная функция
Я[т] = е-а соэю0т, -го<т<+<ю, а> 0, (6)
где ю0 - видимый период колебаний, а а - коэффициент затухания. На практике, как правило, выражение (6) используется для описания формы автокорреляционной функции сейсмических колебаний, представленных моделью случайного процесса. Третьей часто встречающейся автокорреляционной функцией стационарного случайного процесса является АКФ, задаваемая выражением
А + 2В
отсюда, используя графический метод, получаем следующее решение неравенства, представленное на рис. 1.
Я[т] = е “Т(соэю0т +у эшю0 |т|),
А а
-да<т<+да, а> 0, у<—.
(7)
Процесс с такой автокорреляционной функцией не имеет угловых точек, поэтому данная АКФ дифференцируема в среднем квадратическом смысле. Автокорреляционные функции (5)-(7) изображены на рис. 2, е. На рис. 2, з изображены все пять рассмотренные автокорреляционные функции (1, 4, 5-7). Из рис. 2 видно, что форму исходной функции (сигнала) отражает лучше автокорреляционная функция, построенная в виде сплайна первого порядка. Аналогичный вывод делаем на основе доверительного интервала для среднего значения функции автокорреляции. Если в качестве доверительной вероятности /взять 0,95, то доверительный интервал для среднего значения АКФ, построенной по формуле (1), есть (-0,016; 0,072), для функции, полученной с использованием выражения (4), - (-0,056; 0,014).
Аналогичные расчеты были проведены и для модельного сейсмического сигнала, реализация которого изображена на рис. 2, г, для которого построена функция автокорреляции в виде сплайна первого порядка, рис. 2, д. На рис. 2, ж изображены три ранее рассмотренные модельные функции автокорреляции (5-7) и одна, построенная по формулам (4).
Рис. 2. Графики моделей сейсмических сигналов, автокорреляционных функций: а) функция Бесселя первого рода нулевого порядка, г) реализация модельного сейсмического сигнала, б), в), д) функция автокорреляции, вычисленная с использованием выражения (1), (4), АКФ модельного сигнала в виде сплайна первого порядка, е) три модельные АКФ (формулы 5-7); p1t обозначена автокорреляционная функция (5), p2t - АКФ (6), иг - АКФ (7), ж) автокорреляционные функции (4-7) p1t/ p2, rt и Rt - автокорреляционная функция, построенная по формуле (4), з) изображены все пять АКФ (1, 4, 5-7) p1, p2, r, R, ft - обозначение автокорреляционной функции (1)
В таблице приведена сумма квадратов отклонений оцененных значений функции автокорреляции от истинных. В верхней строке представлены номера расчетных формул, во второй - искомые значения сумм.
Таблица. Значение суммы квадратов отклонений истинных значений функции автокорреляции от оценок, найденных по формулам (4)-(7)
Номер формулы 4 5 6 7
Значение суммы 0,054 0,227 0,132 0,217
Из приведенной таблицы можно сделать вывод, что наименьшее значение суммы получается при оценивании функции автокорреляции сплайном первого порядка.
СПИСОК ЛИТЕРАТУРЫ
1. Гордин В.Н., Бабаева Т.Н., Михайлов В.О. О статистической параметризации аномальных потенциальных полей // Геофизический журнал. - 1984. - Т. 6. - № 2. - С. 55-63.
2. Устинова В.Н., Устинова И.Г Статистическая параметризация симметричных геофизических объектов // Известия ТПУ. -2003. - №5. - С. 36-41.
3. Устинова В.Н., Устинова И.Г Дискретные иерархические системы в геофизике // Известия ТПУ. - 2012. - Т. 320. - № 1. -С. 91-97.
4. Никитин А.А. Теоретические основы обработки геофизической информации. - М.: Недра, 1986. - 342 с.
Выводы
Выполнен расчет коэффициентов сплайна первого порядка в приложении к оценке функции автокорреляции. Возможности предложенного метода апробированы на модельном сейсмическом сигнале и на модели, применяемой для аппроксимации сейсмических сигналов. В результате проведенного исследования установлено, что автокорреляционная функция, построенная в виде сплайна первого порядка, лучше отражает форму исходного процесса, нежели модельные автокорреляционные функции. Построенные доверительные интервалы для оценки математического ожидания автокорреляционной функции также свидетельствуют о более точном приближении в виде сплайна, чем в опубликованных работах.
5. Идрисов Ф.Ф., Устинова И.Г. Оценка функции корреляции стационарного случайного процесса при случайном числе измерений // Экономика, технология, предпринимательство. -Вып. 1. - Томск: Изд-во ТГПУ, 2000. - С. 75-79.
6. Корн Г, Корн Т. Справочник по математике для научных работников. - М.: Наука, 1978. - 831 с.
7. Гельфонд А.О. Исчисление конечных разностей. - М.: Наука, 1967. - 375 с.
8. Дьяконов В. МаШсаё 2000. - СПб.: Изд-во «Питер», 2001. -592 с.
9. Пугачев В.С., Синицын И.Н. Стохастические дифференциальные системы. Анализ и фильтрация. - М.: Наука, 1990. - 642 с.
Поступила 07.11.2012 г.
