Эффекты интерференции волновых полей в смежных геофизических средах

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Эффекты интерференции волновых полей в смежных геофизических средах

М.С. Хайретдинов, Г.М. Воскобойникова

Институт вычислительной математики и

математической геофизики СО РАН

Введение

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

Ключевые слова: сейсмические вибраторы, сейсмические волны, акустические волны, гидроакустические волны, интерференция.

При вибрационном зондировании Земли с помощью мощных сейсмических вибраторов наряду с сейсмофизическим полем источники порождают акустическое поле [1]. Возникающие при этом акустические волны при определенных условиях могут распространяться на большие расстояния, приводя к интерференции сейсмических и акустических колебаний. Если же вибросейсмические волны переходят через границу «дно-море», то они могут взаимодействовать с гидроакустическим волнами, порождаемыми промышленными бортовыми источниками вибрационных колебаний. Изучение развивающихся особенностей волновновой интерференции при этом составляет цель настоящей работы.

Интерференция сейсмических и акустических колебаний. Явление возникновения акустических колебаний от сейсмических вибраторов обусловлено тем, что в процессе работы последних на инфразвуковых частотах наблюдаются интенсивные синхронные колебания почвы на значительной площади вблизи вибратора. Такие сейсмические колебания являются вторичным источником, генерирующим акустическое излучение. Применительно к двум граничащим средам - упругому нижнему полупространству и верхнему газообразному - в предположении наличия плоской границы между ними получено соотношение, количественно оценивающее долю энергии акустических волн, излучаемых сейсмическим вибратором [2]. Если принять амплитуду возбуждающей силы источника F, скорости волн в средах и их плотности соответственно V₁, V₂ и ₁, _2, то мощность акустического излучения при V₁= 3 V₂ может быть оценена соотношением

N_a= 3.16 (₂/₁) F^{2 2} /₁V₁³ (1)

Сравнение этой мощности с мощностью упругого излучения продольной волны Р показывает, что мощность первой составляет около 2% от второй, т.е. от уровня упругого излучения.

Распространение акустических волн на расстояния в несколько десятков километров и более возможно благодаря явлению рефракции звуковых волн в атмосфере. Один из механизмов этого явления связан с тем, что при наличии попутного ветра в направлении «источник-приемник» образуется приповерхностный волновой канал, способствующий дальнему распространению звука [3]. При этом оба типа колебаний -сейсмические и акустические - регистрируются с помощью сейсмических датчиков. В качестве примера, подтверждающего сказанное, на рис.1а представлены оба типа волн от вибратора ЦВ-100 [1]: головные сейсмические волны регистрируются на 8-ой сек. и акустические на 145 с. Результаты соответствуют случаю попутного ветра. Направление и сила ветра приведены на рисунке справа. При этом излученные вибратором в среду сигналы представляли собой свип-сигнал в полосе частот 6,25.5-9.57 Гц. Расстояние «вибратор-приемник» составило 50 км. Представленные на рисунке результаты обработки в виде двух коротких импульсов соответствуют взаимокорреляционной свертке двух сигналов - опорного, являющегося копией излученного сигнала, и исходного сигнала, зарегистрированного с помощью сейсмоприемников. Вычисление откликов производится в соответствии с алгоритмом

(2)

Здесь М - число дискретных отсчетов взаимокорреляционной свертки, L - число усреднений по повторяющимся сеансам зондирования, S(t_n) - опорный сигнал с линейной частотной модуляцией (ЛЧМ) вида , параметрами которого являются A(t) - огибающая, f₀ - начальная частота развёртки, - скорость развёртки по частоте, равная где - максимальная частота, Т - длительность развёртки.

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

Рис.1 Взаимодействие сейсмических и акустических волн на расстоянии источник-приемник 50 км.: а)- совместная регистрация сейсмических и акустических волн

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

(3)

(4)

где A_s (t), A_a (t), _s, _a - соответственно огибающие и начальные фазы установившихся сейсмического и акустического сигналов,

Измерение амплитуд установившихся гармонических колебаний на каждой из частот зондирования осуществляется путем оценивания статистики R_N с помощью итерационного квадратурного алгоритма вида [4]

(5)

Здесь , , u_N= S_N +n_N,, где S_N и n_N -текущие значения сигнала и шума. В этом случае среднее значение . В реальной ситуации достаточно, чтобы N Nз, где Nз - множество дискретных отсчетов на участке установления отклика среды. При условии равномерной плотности распределения начальной фазы установившегося колебания в пределах [0, 2]. значение ее определяется в виде .

Пример результата обработки интерференционного сигнала, образованного путем взаимодействия продолжительных во времени сейсмического и акустического сигналов в гармоническом режиме излучения на частоте 6.75 Гц, представлен на рис.1b. Здесь в верхней части рисунка представлен результат накопления по амплитуде слабого гармонического сигнала на фоне микросейсм в соответствии с (4), в нижней части - поведение фазы суммарного сигнала _s + _a.. В случае, когда акустическая волна отсутствует в связи с меняющимися метеоусловиями (рис.1с), на выходе среды имеют место устойчивые оценки параметров сейсмического поля {A, } в установившемся режиме колебаний (рис.1d). Напротив, из рис.1b видно, что поведение тех же параметров колебания, полученного в результате интерференции сейсмического и акустического колебаний, отличается неустойчивостью во времени в силу тех факторов, о которых говорилось выше. В первую очередь это проявляется в поведении оценки фазы (нижняя часть рис.1b) с момента вступления акустической волны, которая приходит с задержкой на время распространения вдоль трассы «источник-приемник», в данном случае равному 145 с. Очевидно, что вариации оценок фазы колебания в рассматриваемом случае могут служить в определенных условиях мерой неустойчивости атмосферы на трассе распространения волн. Из рис.1 следует, что учет явления интерференции сейсмических и акустических колебаний важен в задачах вибросейсмического мониторинга очаговых зон природных катастроф, основанного на измерениях вариаций амплитудно-фазовых характеристик передаточной функции среды

Интерференция сейсмических и гидроакустических волн. Данное явление возникает при взаимодействии двух типов смежных волн - сейсмических и гидроакустических, которые одновременно регистрируются датчиком давления (например, пьезодатчиком в море. Источником сейсмических волн является наземный вибрационный источник, а гидроакустических - бортовые механизмы, например различные вентиляционные установки. При этом последние выступают как источники помех при регистрации полезных сейсмических волн. Такая ситуация возникает при сейсморазведке на углеводороды в шельфовой зоне. Если оба источника излучают гармонические колебания с близкими частотами, то явление интерференции в этом случае может быть представлено в терминах биений двух гармонических колебаний на фоне широкополосного внешнего шума n(t). с параметрами N (0, _n): Как известно], огибающая биений двух компонент равна

(6)

где =₂ -₁ Максимальное значение огибающей достигается при cost=1 и равно А₁ +А_2.В случае равных амплитуд исходных колебаний значение огибающей удваивается и соответственно возрастает уровень полезного сигнала над шумом. Соответственно соотношение “биения”/шум составит Этим обусловлено превышение амплитуды биений над шумом на рис.2.

Приведенная запись интерференционного сигнала соответствует ₁=20Гц - частота излучаемого сейсмического колебания, ₂=21Гц - частота мешающего гидроакустического сигнала. Следующий за картиной биений участок записи соответствует выключению вентиляционной установки и соответственно маскированию внешним шумом полезного сигнала. Таким образом, мы здесь сталкиваемся с явлением усиления на фоне шумов основного сейсмического колебания второстепенным (мешающим) гидроакустическим колебанием близко расположенной частоты.

Рис.2 Явление усиления на фоне шумов основного сейсмического колебания второстепенным акустическим колебанием близко расположенной частоты

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

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

(7)

где - амплитуда и фаза колебания, распространяющегося вдоль i-го луча и меняющиеся случайным образом в пространстве сейсмического поля; n - количество лучей. По отношению к акустическому и гидроакустическому полям эти параметры дополнительно будут меняться случайным образом во времени. Соответственно амплитуда и фаза результирующего установившегося колебания будут принимать случайные значения. Огибающую интерференционного сигнала определим равенством , которая, в общем случае описывается плотностью распределения Накагами (m-распределением), приводимое к известному - распределению [5]:

.(8)

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

Заключение

интерференция волна геофизический

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

Литература

B.M. Glinskii, V.V. Kovalevskii, M.S. Khairetdinov. Relationship of wave fields from powerful vibrators with atmospheric and geodynamic processes. Geology and Geophysics, vol.40, No.3, pp.422-431,1999.

Ю.М. Заславский. К оценке мощности инфразвука, побочно излучаемого в атмосферу при вибрационном просвечивании Земли// Известия АН СССР. Сер. Физика Земли. -1982. - №9. - С. 86-89.

М.А. Исакович. Общая акустика. М.: “Наука”, с.495

М.С. Хайретдинов. Методы распределенной обработки вибросейсмических сигналов. Матер.научно-технич. конф. »Информационные системы и технологии» НГТУ, Новосибирск, 2000, с.607-612.

В.И. Тихонов. Статистическая радиотехника.-М.:Сов.радио.1966.

Размещено на Allbest.ru