CC BY
45
УДК 550.834
ОПРЕДЕЛЕНИЕ ВРЕМЕННОГО СДВИГА МЕЖДУ СЕЙСМИЧЕСКИМИ ТРАССАМИ НА ОСНОВЕ АНАЛИЗА ИХ ФАЗОЧАСТОТНЫХ ХАРАКТЕРИСТИК
Иванченков Виктор Павлович1,
am@am.tpu.ru
Кочегуров Александр Иванович1,
kaicc@tpu.ru
Нгуен Суан Хунг1,
hungnx@tpu.ru
Злобина Елена Вячеславовна1,
elena_zlobina@list.ru
Лаевский Виктор Евгеньевич (V. Geringer)2,
geringer@dhbw-ravensburg.de
1 Национальный исследовательский Томский политехнический университет, Россия, 634050, г. Томск, пр. Ленина, 30.
2 Государственного объединенного университета Баден-Вюртемберг, Германия, D-88045, Friedrichshafen, Fallenbrunnen, 2. E-mail:
Актуальность работы. Процедуры определения временного сдвига, получившие широкое распространение на практике, не всегда гарантируют правильность определения статических поправок. Возникающие при этом ошибки приводят к искажению кинематических параметров отраженных волн, что отрицательно сказывается на всех дальнейших построениях. Поэтому разработка эффективного метода определения временных сдвигов является актуальной проблемой, с успешным решением которой связываются перспективы повышения достоверности построения структур, подготавливаемых для разведки на нефть и газ. Цель исследования: разработать и исследовать надежность алгоритма определения временных сдвигов между сейсмическими трассами на основе оптимальной и субъоптимальной обработки фазочастотных характеристик отраженных волн, регистрируемых при поиске нефтяных и газовых месторождений.
Методы исследования: статистический синтез и анализ оптимальных решающих правил, цифровая обработка сигналов, методы дискретного преобразования Фурье, математическое моделирование и вычислительный эксперимент. Результаты. Показана возможность извлечения полезной информации из фазочастотных характеристик сейсмических сигналов для получения эффективных оценок временных сдвигов между трассами. На основе теории статистических решений построен оптимальный алгоритм определения временных сдвигов по мгновенным взаимным фазовым спектрам. Для практической реализации алгоритма предложен переход от оптимальной к субоптимальной обработке компонент мгновенного взаимного фазового спектра путем равновесного суммирования. Показано, что равновесная обработка приводит к уменьшению суммарного отношения сигнал/шум для практически используемой полосы частот почти вдвое. Однако даже такие потери позволяют получать достаточно надежные оценки временных сдвигов при наличии интенсивных нерегулярных помех и симметричных распределениях временных сдвигов с большой дисперсией.
Ключевые слова:
Сейсмические волны, статические поправки, функция взаимной корреляции, взаимные фазовые спектры, фазочастотное прослеживание.
Введение
В процессе обработки сейсмических наблюдений учет статических поправок, связанных с рельефом дневной поверхности и неоднородностями верхней части разреза, является одной из важнейших процедур интерпретации данных сейсморазведки. Качество выполнения этой процедуры непосредственно влияет как на надежность выделения и корреляцию отраженных волн, так и на достоверность последующих глубинных построений. Наличие в наблюденных временах статических сдвигов препятствует синфазному суммированию полезных волн, и именно поэтому первой важной проблемой в обработке сейсмических данных явля-
ется устранение этих сдвигов посредством введения процедур коррекции статических поправок [1-6].
Одним из основных этапов при проведении коррекции статических поправок является определение временного сдвига между сейсмическими трассами. Различаясь в деталях, большинство алгоритмов автоматического определения временных сдвигов построены на основе вычисления функций взаимной корреляции (ФВК) между трассами сейсмограммы или трассами сейсмограммы и набором эталонных трасс [7-10]. Вместе с тем такие способы определения временного сдвига, получившие широкое распространение на практике, не всегда гарантируют правильность определения
статических поправок. Возникающие при этом ошибки приводят к искажению кинематических параметров отраженных волн, что отрицательно сказывается на всех дальнейших построениях [7].
Таким образом, разработка новых методов определения временных сдвигов является актуальной проблемой, с успешным решением которой связываются перспективы повышения достоверности построения структур, подготавливаемых для разведки на нефть и газ.
Одним из возможных путей решения данной проблемы является использование новых информативных параметров сейсмических сигналов, а именно их фазовых спектров. В фазу сигнала, точнее в сложный закон изменения его фазового спектра, заложена полная информация о временном положении сигнала, поэтому оптимальная обработка фазочастотной характеристики (ФЧХ) сигнала реализует оптимальный метод определения его временного положения [11-13]. Вопросам извлечения информации из фазочастотных характеристик сейсмических сигналов для оценки временных сдвигов между трассами был посвящен ряд работ [14-16]. Например, в [14] показано, что оценки временного сдвига т можно определить из следующего выражения:
(а) -фI (а)]
т = -
X
а
(1)
Функция правдоподобия в этом случае может быть представлена в виде
Цт) = W (ДфЮ, ДфК),..., Аф(ат) / т), (2)
где Ж(Дф(ок)/т) - многомерная условная плотность вероятности Дф(ак); Дф(а1),Дф(а2),...,Дд)(ая) - отклонения ФЧХ смеси сигнала с помехой от ФЧХ сигнала при различных значениях т.
Одномерная плотность вероятности для Дф(ак) приведена в [20]. Тогда, рассматривая ее как функцию от т и учитывая, что т является неэнергетическим параметром, запишем:
2л
exp
р2(тк Л
W (Дф(ак)/ т) = Р(ак)
) + 42м
^СДф(тк) -ткт) х
ХФ[Р(®£)cos(ДФ(ак) -акт)] х
( р2(®.) ^ х exp |--8ш(Дф(тк) -ткт) ) (3)
где Ф - преобразование Лапласа.
Теперь, выбирая Да так, чтобы Дф(ак) (к=1,т) были статистически независимы, выражение (2) можно переписать в виде
1(т) = ^ {ДфК )/ т},
(4)
П
где фэт.(а) - фазочастотная характеристика эталонной трассы; ф(а) - ФЧХ 1-й трассы сейсмограммы.
Однако, к сожалению, подобные алгоритмы до настоящего времени не нашли широкого применения при обработке и интерпретации реальных сейсмических материалов, что, по-видимому, связано прежде всего с трудностью идентификации реальных ФЧХ сейсмических сигналов [15-17], а вследствие этого - небольшим количеством исследований в области синтеза и анализа фазочастот-ных алгоритмов [18, 19], для обработки сейсмической информации.
Фазочастный алгоритм определения
временных сдвигов
В настоящей работе производится дальнейшее изучение информационных возможностей ФЧХ сейсмических трасс для определения временных сдвигов, создание на этой основе новых алгоритмов и исследования их в ситуациях, где применяемые на практике алгоритмы на основе ФВК не обеспечивают требуемой точности оцениваемых параметров.
Сначала построим процедуру определения временного положения сейсмических сигналов т, когда оптимальной обработке подвергаются лишь их ФЧХ. Для этого примем, что сигнал слабый (р(ак)<1, р(ак) - отношение сигнал/помеха на частоте ак) и регистрируется на фоне аддитивного нормального шума.
где W(Дф(аk)/т) - определяется из (3); т =--
Да
число анализируемых частотных компонент; Да -шаг дискретизации по частоте.
Прологарифмируем (4) и, опуская члены, не зависящие от т, для случая слабого сигнала (р(ак)<1) получим:
т
1П Цт) =Хр2(®к Дф(®к)-®кт). (5)
к =1
Тогда оптимальную оценку топт. можно найти из решения уравнения
д
—1п Цт) дт
= 0.
(6)
Уравнение (6), где 1иЦт) описывается (5), не позволяет в явном виде получить топт.. Однако, если учесть, что т принимает конечное число значений из некоторого интервала то топт=т0, при котором
1п Ь(т0) = max{ln Ь(т)},
теК (2]. (7)
Покажем, что критерий (7) является также оптимальным и для сильных сигналов. При сильном сигнале флуктуации ФЧХ незначительны, поэтому (5) можно представить в виде
1ПЬ(т) =Хр\®к)(Дф(®к) -®кт)2.
(8)
Тогда из решения уравнения (6) с учетом (8) получим:
к=1
а
к =1
Топт. —
£р2(тк )юк Дф(®к )
т
£р2(тк )тк2
(9)
Выражение (9) описывает известную оценку временного положения сильного сигнала [13], когда оптимальной обработке подвергается его ФЧХ. Таким образом, критерий (5) является оптимальным как для случая слабых сигналов, так и для случая сильных сигналов, в то время как критерий (9) оптимальным для слабых сигналов не является. При этом дисперсия оценки Топт. по критерию (9) для сильных сигналов совпадает с предельно достижимой и составляет:
я! — [^(шо)2]-1,
где чк — £р1(тк) -
суммарное отношение сиг-
1п ^(т) — £ С05( Дф(Юк ) -®кТ>-
(11)
V — (чк / ?Ё)2 -1, (12)
где дк - суммарное отношение сигнал/помеха для критерия (5); д'к - суммарное отношение сигнал/помеха для критерия (11).
Тогда, как следует из [18], ^тм для (12) составляет:
т
— £{4к->/&"—1 )2 -1
или
(С)2 — -
(чк )2
£ (VI )2
к—1
Теперь, заменяя в (10) (д'Е) на (дЁ), получим дисперсию оценки т для критерия (11):
т
4£(>/к -4к-л)2
нал/помеха; от - среднеквадратическая ширина спектра сигнала.
Принимая Дф(юк)=ф8Т.(юк)-ф;(юк), нетрудно показать, что выражение (1) является частным случаем критерия (9) и дает оптимальную оценку временного положения сильного сигнала только для случая, когда отношение сигнал/помеха в анализируемой полосе частот является константой.
Для слабых сигналов дисперсия оптимальной оценки, по сравнению со случаем сильного сигнала, возрастает в 4/жраз:
2 4 4 о ——2-- —-2-- (10)
Т». щ1(от)2 (чк) (от)
Для реализации оптимального алгоритма, как видно из (5), необходимо знать распределение отношения сигнал/помеха р(т) в полосе частот П. На практике получить достаточно надежные оценки рр(т) для ®ёП весьма сложно. Поэтому предлагается принять в (5) р(ш)=ео^ для теП, т. е. перейти к равновесной обработке компонент фазового спектра. В этом случае, естественно, можно говорить только о квазиоптимальности получаемых результатов.
Оценим точность получаемых оценок временного положения сигналов Т при переходе к равновесной обработке. Для этого примем в выражении (5) р(ш)=ео^=1 (теП) и получим:
Как следует из (10), дисперсия оценки топреде-ляется суммарным отношением сигнал/помеха д'к и среднеквадратической шириной спектра сигнала от. Так как при переходе от (5) к (11) от не изменяется, найдем чему будет равно дк при использовании критерия (11). Для этого воспользуемся методикой оценки потерь V при равновесной обработке многочастотного сигнала [21]. В нашем случае имеем:
Т лд^(от)
Как и следовало ожидать, потеря оптимальности приводит к снижению точности получаемых оценок, причем о- зависит от числа учитываемых гармоник в спектре сигнала. С ростом т дисперсия оценки о- по отношению к предельно достижимой для учитываемого числа частотных компонент увеличивается, так, при т=3 77тах=0,27, а при т=10 ?7тах=0,754. Однако даже при т=20 точность получаемых оценок достаточно высока. Рассмотрим практическую реализацию критерия с равновесной обработкой (11) для определения временных сдвигов между сейсмическими трассами.
Как следует из (11), для его реализации необходимо вычислить Дф(тк)=ф8т.(тк)-ф;(тк) для (к=1,т) при различных - диапазон возмож-
ных значений временных сдвигов. Эту процедуру можно осуществить либо путем вычисления аргументов преобразования Фурье - преобразований отдельно эталонной трассы и трассы сейсмограммы, либо вычислением взаимного фазового спектра между этими трассами. Второй путь требует дополнительного вычисления ФВК между трассами, однако он позволяет получить более надежные оценки Др(тк) для (к=1,т). Поэтому остановимся на этом пути. Тогда процедура определения временного сдвига, в соответствии с (11), сводится к вычислению мгновенного взаимного фазового спектра (МВФС) в интервале Мгновенный
фазовый спектр является функцией, изменяющейся во времени и по частоте, и определяется как Фэт. I (т,т) — Дф(т)-тт —
Т+Т г — Аг% | |Д
■эт.,/ V)е - - Т -Т| Ж
(13)
где Д8т.,;(0 - ФВК между эталонной трассой и к-й трассой сейсмограммы; Т - интервал расчета взаимного фазового спектра (окно анализа), включающий 3-5 видимых периодов ФВК.
к—1
к —1
к—1
к—1
к—1
Таким образом, как следует из (13), мгновенный взаимный фазовый спектр находится «скользящим» интегрированием: интервал интегрирования имеет постоянную длину Т и перемещается по ФВК с шагом Дt Д - шаг дискретизации трассы) в окрестности ее глобального экстремума; для каждого мгновенного положения окна анализа (конкретного те[^,У) вычисляется взаимный фазовый спектр, при этом базисные функции, независимо от
положения окна, считаются на отрезке
Т Т 2 , 2
В качестве оценки временного сдвига между трассами принимается т, при котором достигается максимум логарифма функции качества [16] для рассчитанного мгновенного взаимного фазового спектра в интервале те[^^2]. Нетрудно увидеть, что определяемая таким образом оценка временного сдвига является точкой центра симметрии окрестности ФВК, включающей глобальный экстремум. Поскольку центральная симметричная оценка носит интегральный характер, то она должна быть более устойчивой, чем положение максимума ФВК [22, 23].
Теперь проведем исследование эффективности данного критерия для определения временных сдвигов в ситуациях, где применяемые на практике алгоритмы, основанные на расчете и анализе функций взаимной корреляции (ФВК), не обеспечивают требуемой точности получаемых оценок временных сдвигов.
Рассмотрим ситуации, при которых возникают существенные погрешности при оценке сдвига по положению максимумов ФВК [7]:
• симметричные распределения временных сдвигов тс дисперсией, превышающей половину видимого периода волны;
• распределения с резкой асимметрией;
• наличие интенсивных нерегулярных помех. Будем считать, что в процессе обработки сейсмических данных производится вычисление ФВК между трассами сейсмограммы и набором эталонных трасс. Тогда взаимная спектральная плотность между трассой сейсмограмм и эталонной трассой без учета нерегулярных помех имеет вид [7]
(а) = уа 2( А») Т»,
где у - среднее число сейсмических сигналов, появляющихся в единицу времени; -2=М{а2) - математическое ожидание квадратов амплитуд сигналов; Д(а) - амплитудный спектр сигнала; Тт(а) -характеристическая функция временного сдвига т (сопряженное преобразование Фурье от плотности вероятностей W (т)).
В этом случае взаимный фазовый спектр запишется как
Дф(а) = arctg
Re{Т»}
(14)
Если дисперсия скоростей полезных волн в интервале вычисления ФВК равна нулю, то Дф(а) = arctg[tgат] = ат,
и средний сдвиг между трассами т - определяется точно из условия максимизации (11):
т = т.
При наличии дисперсии волн по скорости проанализируем последовательно ситуации симметричного распределения временных сдвигов и распределения с резкой асимметрией. Возьмем в качестве примера симметричного распределения равномерное распределение временных сдвигов. В[7] показано, что при параметре разброса 5>Твид/2 определение сдвига по положению максимума ФВК становится ненадежным. Запишем характеристическую функцию для равномерного распределения:
Тт (а) = 8^т е(28а) еа, (15)
и, подставляя (15) в (14), а найденное выражение в (11), окончательно получим:
т
1ПЬ1(т) =Хcos(a т -ат). (16)
к=1
Из (16) видно, что максимум 1и!1( т) независимо от величины 5 достигается при т= т, то есть сдвиг между трассами определяется верно. Следовательно, в данном случае предлагаемый алгоритм устраняет ошибки, возникающие при определении сдвига по положению максимума ФВК.
Результаты экспериментов
Исследование эффективности алгоритма определения сдвигов по МВФС путем максимизации (11) для распределений сдвигов с резкой асимметрией проводилось с помощью статистического моделирования на ЭВМ. При этом в качестве примера распределения с резкой асимметрией было выбрано распределение Рэлея.
На рис. 1. и 2 приведены оценки математического ожидания и дисперсии временного сдвига в зависимости от величины параметра разброса в распределении Рэлея, найденные, соответственно, по положению максимума ФВК и предлагаемому алгоритму. Из рисунков видно, что для обоих алгоритмов получаемые оценки оказываются смещенными относительно среднего временного сдвига между трассами, и тем в большей степени, чем больше разброс распределения Рэлея. Однако смещение оценок и их дисперсии у предлагаемого алгоритма (рис. 1, б и 2, б) несколько ниже, чем в случае использования ФВК (рис. 1, а и 2, а). Следовательно, предлагаемый алгоритм при распределениях сдвигов с резкой асимметрией не устраняет полностью ошибок, свойственных ФВК в этом случае.
Исследование эффективности предлагаемого алгоритма определения сдвигов при интенсивных нерегулярных помехах проводилось также путем
Рис. 2. Оценка дисперсии временных сигналов между трассами: а) ВКФ; б) МВФС Fig. 2. Assessment of dispersion of temporary signals between traces: a) MCF; b) IMPhS
m (mc)
AT
12 10 8 6 4 2 0
тЛт(мс)
95% У И*
S --
✓ — — _ - -
95%
10
9
95% ✓
> S
" ■ »*" """
✓ ✓
Л Ф
95%
2 4 6 7 8 10
8 10
a/a 6/b
Рис. 1. Оценка математического ожидания временных сдвигов между трассами: а) ФВК; б) МВФС
Fig. 1. Assessment of a population mean of temporary shifts between traces: a) mutual correlation functions (MCF); b) instantaneous mutual phase spectra (IMPhS)
О - (MC2)
Дт 4 y
70 60 50 40 30 20 10 0
95% У
У ✓
✓ ✓
✓ ✓ ✓
f У ✓ у ? ^^^ 4 «ч» ** * 95%
У У . ** **
3 4 6
a/a
8 10
a - (mc2)
Дт 4 y
60 50 40 30 20 10 0
95%
✓
у S г
S S У
✓ ✓ ✓ 95%
** **
3 4 6
б/Ь
8 10
статистического моделирования на ЭВМ. На рис. 3 показаны оценки дисперсии временного сдвига, найденного по двум сопоставляемым алгоритмам, когда анализируемые трассы содержали аддитивную помеху. Сравнение результатов показывает, что алгоритм определения сдвигов по МВФС (рис. 3, б) обеспечивает более высокую помехоустойчивость при низких отношениях сигнал/помеха (^<2), чем алгоритм на основе ВКФ (рис. 3, а).
Этот эффект объясняется тем, что, как уже отмечалось выше, центральная симметричная оценка (11) носит интегральный характер и, следовательно, является более устойчивой при нерегуляр-
ных помехах по сравнению с положением глобального максимума ФВК.
Заключение
Таким образом, проведенные исследования показали, что предлагаемый фазочастотный алгоритм позволяет в ряде практически важных случаев, а именно при наличии интенсивных нерегулярных помех и симметричных распределений временных сдвигов с большой дисперсией, существенно повысить точность и достоверность определения временных сдвигов по сравнению с методами, базирующимися на анализе ФВК.
О • (мс )
Дт 4 7
40 35 30 25 20 15 10 5 0
V 95%
■ч ч ч
V V
Ч ч
N. ч» N^ NJ
95% V
•s.
0,5
a/a
а1- (мс )
Дт х '
30 25 20 15 10 5 0
% N ^95%
-w '■ч
* ^
95% ^ ~ — ---
0,5
б/Ъ
Рис. 3. Оценка дисперсии временных сдвигов между трассами при наличии помех: а) ВКФ; б) МВФС Fig. 3. Assessment of dispersion of temporary shifts between traces at noises: a) MCF; b) IMPhS
СПИСОК ЛИТЕРАТУРЫ
1. Козырев В.С., Коротков И.П., Жуков А.П. Способы учета сложно построенной верхней части разреза в Западной Сибири // Технологии сейсморазведки. - 2014. - № 3. - C. 66-78.
2. Kearey Ph., Brooks M., Hill I. An Introduction to Geophysical Exploration. - 3d ed. -US: Wiley-Blackwell, 2002. - 272 p.
3. Taner M.T., Koehler F., Alhilali K.A. Estimation and correction of near-surface time anomalies // Geophysics. - 1974. - V. 39. -P. 441-463.
4. Douma H., Haney M. Surface-wave inversion for nearsurface shear-wave velocity estimation at Coronation field // 81st SEG Annual Meeting. Expanded Abstracts. - San Antonio, 2011. -P.1411-1415.
5. Finite Difference Modeling to Evaluate Acquisition and Processing Parameters on Resolution of Land Seismics / O. Adamovich, K. Epov, I. Korotkov, P.C.H. Veeken // 72nd EAGE Conference and Exhibition incorporating SPE EUROPEC, 2010. - Barcelona, 2010. DOI: 10.3997/2214-4609.201401101.
6. Korotkov I., Veeken P.C.H., Kozyrev V. Interactive Static Corrections to Avoid Mis-interpretation of Seismic Data Due to Near-surface Velocity Changes // 72nd EAGE Conference and Exhibition incorporating SPE EUROPEC, 2010. - Barcelona, 2010. DOI: 10.3997/2214-4609.201401294.
7. Яновский А.К. Использование функций взаимной корреляции сейсмических трасс для увязки сейсмических наблюдений // Дискретная корреляция сейсмических волн. - Новосибирск: Наука, 1971. - C. 50-72.
8. An interactive solution for resolving mid-wavelength statics anomalies / V. Kozyrev, Al. Zhukov, I. Korotkov, Ar. Zhukov // The Leading Edge. - 2004. - V. 23. - № 3. - P. 270-274.
9. Салищев М.В., Сысоев А.П. Уточнение модели коррекции статистических поправок // Технологии сейсморазведки. -2014.- №3. - C. 61-65.
10. Shaun Strong, Steve Hearn. Statics correction methods for 3D converted-wave (PS) seismic reflection // Exploration Geophysics. URL: http://dx.doi.org/10.1071/EG15115 (дата обращения: 05.03.2016).
11. Худяков Г.И. О потенциальной точности определения временного положения флюктуирующих сигналов // Вопросы радио-
электроники. Серия: Общие вопросы радиоэлектроники. -
1984. - Вып. 8. - C. 55-60.
12. Иванченков В.П., Кочегуров А.И. Определение временного положения сейсмических сигналов по оценкам их фазочастот-ных характеристик // Геология и геофизика. - 1988. - № 9. -С. 77-83.
13. Kochegurov A.I., Kochegurova E.A., Kupina N.A. Detection accuracy of the temporary state of complex signals using phase-frequency tracking methods with equilibrium and non-equilibrium processing // Advances in Intelligent Systems and Computing. - 2015. - V. 342. - P. 27-36.
14. Зенов А.А. Вычисление оптимального временного сдвига между сейсмическими трассами //Разведочная геофизика. -1977.- Вып. 76. - C. 18-21.
15. Данквардт Э., Патцер У. Проблемы использования фазовых спектров сейсмических данных // Труды 30-го Международного географического симпозиума. Стендовые доклады. - М.,
1985. - Ч. 2. - C. 132-142.
16. Митрофанов Г.М. Обработка фазовых спектров многоканальных сейсмограмм // Геология и геофизика. - 1986. - № 10. -С. 99-109.
17. Каныгин М.Ю., Петрова Н.В. Согласование фазовых спектров сигналов при формировании обобщённых сейсмических проектов // Технологии сейсморазведки. - 2009. - № 1. -С. 56-59.
18. Методы фазочастотного прослеживания отраженных волн и их применения в задачах обработки сейсмической информации / В.П. Иванченков, А.И. Кочегуров, Н.А. Купина, О.В. Орлов // Технология сейсморазведки. - 2013. - № 3. - С. 5-10.
19. Иванченков В.П., Кочегуров А.И., Орлов О.В. Применение методов фазочастотного анализа сейсмических волн для изучения околоскважинного пространства // Нефтяное хозяйство. -2015. - № 10. - С. 58-63.
20. Гольдин С.В. К теории спектрально-статистического метода обработки сейсмограмм // Геология и геофизика. - 1976. -№ 1.- C. 138-149.
21. Поиск, обнаружение и измерение параметров сигналов в радионавигационных системах / под ред. Ю.М. Казаринова. -М.: Советское Радио. - 1975. - 296 с.
22. Hezo A., Schwatz S. Alternatives cross correlation for time delay estimation // ICA SSP 84; IEEE Int. conf. Acoust., Speech, and Signal Process. - San Diego, Calif., 1984. - V. 1. -P. 15.411-15.414.
23. Miller L.E., Lee J.S. Error analysis of time delay estimation using a finite integration time correlator // IEEE Transactions on Acoustics, Speech, and Signal Processing. - 1981. - V. 29. -P. 490-496.
Поступила 26.04.2016 г.
Информация об авторах
Иванченков В.П., кандидат технических наук, доцент кафедры прикладной математики Института кибернетики Национального исследовательского Томского политехнического университета.
Кочегуров А.И., кандидат технических наук, доцент кафедры прикладной математики Института кибернетики Национального исследовательского Томского политехнического университета.
Нгуен Суан Хунг, аспирант кафедры прикладной математики Института кибернетики Национального исследовательского Томского политехнического университета.
Злобина Е.В., магистрант кафедры прикладной математики Института кибернетики Национального исследовательского Томского политехнического университета.
ЛаевскийВ.Е. (V. Geringer),кандидат технических наук, зав. лаб. «Регенеративные источники энергии и автомобильная электроника» Государственного объединенного университета Баден-Вюртемберг.
UDC 550.834
DEFINITION OF TEMPORARY SHIFT BETWEEN SEISMIC TRACES ON THE BASIS OF THE ANALYSIS OF THEIR PHASE-FREQUENCY CHARACTERISTICS
Viktor P. Ivanchenkov1,
am@am.tpu.ru
Alexander I. Kochegurov1,
kaicc@tpu.ru
Nguen Xuan Hung1,
hungnx@tpu.ru
Elena V. Zlobina1,
elena_zlobina@list.ru
Viktor Geringer2,
geringer@dhbw-ravensburg.de
1 National Research Tomsk Polytechnic University, 30, Lenin Avenue, Tomsk, 634050, Russia.
2 Baden-Wuerttemberg Cooperative State University, Germany, D-88045, Friedrichshafen, Fallenbrunnen, 2.
The relevance of the study. The procedures for determining temporary shift, which are widely spread in practice, do not always guarantee the accuracy in defining static corrections. The arising errors lead to a distortion of kinematic parameters of the reflected waves. This affects adversely all subsequent constructions. Therefore, development of the effective method for determining the temporary shift is the relevant issue. The prospects of increasing the reliability of designing the structures, prepared for oil and gas exploration, are related to the successful solution of the problem.
The aim of the study is to develop and explore the reliability of the algorithm for determining the temporary shifts between seismic traces based on the optimal and suboptimal processing of phase-frequency characteristics of the reflected waves recorded when searching for oil and gas fields.
The methods used in the study: statistical synthesis and analysis of optimal decision rules, digital signal processing, methods of discrete Fourier transform, mathematical modeling and computer experiment.
The results. The authors showed the possibility of extracting useful information from the phase-frequency characteristics of seismic signals for obtaining the effective estimates of temporary shifts between traces. Based on the theory of statistical solutions the optimal algorithm for determining the temporary shifts by the instantaneous mutual phase spectra was built. For practical implementation of the algorithm, the authors proposed a transition from optimal to suboptimal processing of the components of the instantaneous mutual phase spectra by summing the equilibrium. It is shown that the equilibrium processing leads to almost twice reduction of the total signal/noise ratio for practically used strip of frequency. Even such losses allowed obtaining rather reliable estimates of temporary shifts at the intense irregular interference and symmetric distributions of temporary shifts with large dispersion.
Key words:
Seismic waves, static corrections, function of mutual correlation, mutual phase ranges, phase-frequency tracing.
REFERENCES
1. Kozyrev V.S., Korotkov I.P., Zhukov A.P. Method of accounting complicated constructed top part of a section in Western Siberia. Tekhnologii seysmorazvedki, 2014, vol. 3, pp. 66-78. In Rus.
2. Kearey Ph., Brooks M., Hill I. An Introduction to Geophysical Exploration. 3'd ed. US, Wiley-Blackwell, 2002. 272 p.
3. Taner, M.T., Koehler, F. and Alhilali, K.A. Estimation and correction of near-surface time anomalies. Geophysics, 1974, vol. 39, pp. 441-463.
4. Douma H., Haney M. Surface-wave inversion for nearsurface shear-wave velocity estimation at Coronation field. 81st SEG Annual Meeting, Expanded Abstracts. San Antonio, 2011. pp. 1411-1415.
5. Adamovich O., Epov K., Korotkov I., Veeken P.C.H. Finite Difference Modeling to Evaluate Acquisition and Processing Parameters on Resolution of Land Seismics. 72nd EAGE Conference and
Exhibition incorporating SPE EUROPEC. Barcelona, 2010. DOI: 10.3997/2214-4609.201401101.
6. Korotkov I., Veeken P.C.H., Kozyrev V. Interactive Static Corrections to Avoid Mis-interpretation of Seismic Data Due to Near-surface Velocity Changes. 72nd EAGE Conference and Exhibition incorporating SPE EUROPEC, 2010. Barcelona, 2010. DOI: 10.3997/2214-4609.201401294.
7. Yanovskiy A.K. Ispolzovanie funktsiy vzaimnoy korrelyatsii se-ysmicheskikh trass dlya uvyazki seysmicheskikh nablyudeniy [Use of functions of mutual correlation of seismic routes for coordination of seismic supervision]. Diskretnaya korrelyatsiya seysmicheskikh voln [Descrete correlation of seismic waves]. Novosibirsk, Nauka Publ., 1971, pp. 50-72.
8. Kozyrev V., Zhukov Al., Korotkov I., Zhukov Ar. An interactive solution for resolving mid-wavelength statics anomalies. The Leading Edge, 2004, vol. 23, no. 3, pp. 270-274.
9. Salishchev M.V., Sysoev A.P. Specification of correction model of statistical amendments. Tekhnologii seysmorazvedki, 2014, vol. 3, pp. 61-65. In Rus.
10. Shaun Strong, Steve Hearn. Statics correction methods for 3D converted-wave (PS) seismic reflection. Exploration Geophysics. Available at: http://dx.doi.org/10.1071/EG15115 (accepted: 5 March 2016).
11. Khudyakov G.I. O potentsialnoy tochnosti opredeleniya vremen-nogo polozheniya fluktuiruyushchikh signalov [About the potential accuracy of defining provisional regulations of fluctuating signals]. Voprosy radioelektroniki. Seriy: Obshchie voprosy radio-elektroniki, 1984, no. 8, pp. 55-60.
12. Ivanchenkov V.P., Kochegurov A.I. Opredelenie vremennogo po-lozheniya seysmicheskikh signalov po otsenkam ikh fazochastot-nykh kharakteristik [Determination of temporal position of seismic signals to measure their phase response]. Russian Geology and Geophysics, 1988, vol. 9, pp. 77-83.
13. Kochegurov A. I., Kochegurova E. A., Kupina N. A. Detection accuracy of the temporary state of complex signals using phase-frequency tracking methods with equilibrium and non-equilibrium processing. Advances in Intelligent Systems and Computing, 2015, vol. 342, pp. 27-36.
14. Zenov A.A. Vychislenie optimalnogo vremennogo sdviga mezhdu seysmicheskimi trassami [Calculation of optimum temporary shift between seismic routes]. Razvedochnay geofizika, 1977, no. 76, pp. 18-21.
15. Dankvardt E., Pattser U. Problemy ispolzovaniy fazovykh spek-trov seysmicheskikh dannykh [Problems of using phase spectra of seismic data]. Trudy 30 Mezhdunarodnogo geograficheskogo sim-poziuma [Proc. of the 30 International geographical symposium]. Moscow, 1985. No. 2, pp. 132-142.
16. Mitrofanov G.M. Obrabotka fazovykh spektrov mnogokanalnykh seysmogramm [Processing of phase ranges of multichannel seis-
mographic records]. Russian Geology and Geophysics, 1986, vol. 10, pp. 99-109.
17. Kanygin M.Yu., Petrova N.V. Soglasovanie fazovykh spektrov signalov pri formirovanii obobshchennykh seysmicheskikh pro-ektov [Coordination of signal phase ranges when forming generalized seismic projects]. Tekhnologii seysmorazvedki, 2009, vol. 1, pp. 56-59. In Rus.
18. Ivanchenkov V.P., Kochegurov A.I., Kupina N.A., Orlov O.V. Phase/frequency picking of reflection events: method and application to cmp data processing. Tekhnologii seysmorazvedki, 2013, vol. 3, pp. 5-10. In Rus.
19. Ivanchenkov V.P., Kochegurov A.I., Orlov O.V. Application of the phase-frequency analysis methods of seismic waves to study borehole environment. Neftyanoe khozyaystvo, 2015, vol. 10, pp. 58-63. In Rus.
20. Goldin S.V. K teorii spektralno-statisticheskogo metoda obrabot-ki seysmogramm [To the theory of a spectral and statistical method of processing seismographic records]. Russian Geology and Geophysics, 1976, vol. 1, pp. 138-149.
21. Poisk, obnaruzhenie i izmerenie parametrov signalov v radionavi-gatsionnykh sistemakh [Search, detection and measurement of parameters of signals in radio navigational systems]. Ed. by Yu.M. Kazarinov. Moscow, Sov. Radio Publ., 1975. 296 p.
22. Hezo A., Schwatz S. Alternatives cross correlation for time delay estimation. ICA SSP 84; IEEE Int. conf. Acoustics, Speech, and Signal Process. San Diego, Calif., 1984. Vol. 1, pp. 15.411-15.414.
23. Miller L.E., Lee J.S. Error analysis of time delay estimation using a finite integration time correlator. IEEE Transactions on Acoustics, Speech, and Signal Processing, 1981. Vol. 29, pp. 490-496.
Received: 26 April 2016.
Information about the authors
Viktor P. Ivanchenkov, Cand. Sc., associate professor, National Research Tomsk Polytechnic University. Alexander I. Kochegurov, Cand. Sc., associate professor, National Research Tomsk Polytechnic University. Nguen Xuan Hung, postgraduate student, National Research Tomsk Polytechnic University. Elena V. Zlobina, graduate student, National Research Tomsk Polytechnic University. Viktor Geringer, Dr. Sc., Head of the Laboratories, Baden-Wuerttemberg Cooperative State University.
