Рассеяние поля звукового пучка тонким слоем однородных включений Текст научной статьи по специальности «Физика»

КУСТИКА

шашг

Электронный журнал «Техническая акустика» http://www .ejta.org

2014, 7

А. М. Гаврилов, О. А. Савицкий, Д. В. Тимошенко

Инженерно-технологическая академия Южного федерального университета, 347900 г. Таганрог, пер. Некрасовский, 44

gavr_am@mail.ru, osav66@mail.ru, dmitrytim@yandex.ru

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

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

ВВЕДЕНИЕ

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

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

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

Решение перечисленных задач требует разработки модели рассеяния акустических волн слоем малоразмерных объектов. Предложенная в настоящей статье модель основывается на результатах, полученных в работах [1-3]. В частности, в [2] показано, что морские донные осадки с одной стороны, в силу своей водонасыщенности, имеют практически те же акустические характеристики, что и вода; с другой стороны,

Рассеяние поля звукового пучка тонким слоем однородных включений

Получена 28.04.2014, опубликована 29.05.2014

осадочный слой обладает достаточной вязкостью, поэтому присутствующие в нем неоднородности, в частности газовые пузырьки, являются практически неподвижными. Приведенные факты позволяют в качестве одного из предположений при постановке задачи в простейшем случае представить газонасыщенный осадочный слой в виде равномерно распределенного на поверхности дна слоя однородных рассеивателей. Также в [2, 3] показано, что часть звуковых волн, проходящая через слой рассеивателей и достигающая твердой поверхности дна, поглощается и не оказывает существенного влияния на картину рассеянного поля, что повышает достоверность описания характеристик рассеивателей через параметры рассеянного поля. Для описания характера взаимодействия между рассеивателями в предлагаемой модели используется метод самосогласованного поля, особенности применения которого к средам, содержащим неоднородности, изложены в работе [4].

1. ПОСТАНОВКА ЗАДАЧИ

Пусть на плоской поверхности z = 0 произвольным образом распределены одинаковые рассеиватели. На плоскость z = 0 под углом 60 к нормали падает звуковая волна в форме осесимметричного продифрагировавшего гауссова звукового пучка, как показано на рис. 1. Необходимо определить амплитуду и фазу рассеянного поля в произвольной точке верхнего полупространства z > 0, их связь с параметрами падающей волны и характеристиками рассеивателей.

2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ

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

Т (r, R, t) = Т0г (R) • exp[t - kr)], (1)

где ш - круговая частота; t - время; k - волновой вектор; r(x, y, z) - радиус-вектор

точки наблюдения M ; Y0¿ (R) - множитель, учитывающий распределение амплитуды

вдоль фронта волны; R - радиус-вектор отстояния точки M относительно оси пучка OL.

Координатные оси Ox и Oz декартовой системы развернем так, чтобы ось пучка OL оказалась в плоскости xOz . Тогда волновой вектор k будет иметь компоненты:

k = k (- k0 sin0„;0; - k0cos90), (2)

где k0 = k . Обозначим через r0 радиус-вектор точки K, которая является проекцией точки M на ось пучка, тогда r02 + R2 = r2. Выражая длину вектора r0 через скалярное произведение вектора r(x,y, z) и единичного вектора ( - kjk0 ), получим

r0(x,z) = -(k/k0; r)= xsin90 + zcos90. (3)

Рис. 1. Геометрия задачи

Используя соотношения (2) и (3), получим выражение для R2 (квадрат модуля радиус-вектора отстояния):

2 2 2/2 2 2 \{ \2 / \2 2 R = r -r0 = [x + y + z )-(xsin90 + zcos90) = (x• cos90 -z• sin90) + y . (4)

В гауссовом пучке амплитудный множитель в выражении падающей волны (1) имеет вид:

Т (R) = Aexp(- R 2/R0¡). (5)

Здесь параметр R0 - ширина звукового пучка по уровню 1/e. С учетом (4), получим:

Т(x, y, z) = AexP

(x • cos 00 - z • sin 00 )2 +

y

R2

(6)

В плоскости расположения рассеивателей ( z = 0) выражение для потенциала примет вид:

Т (x, y, z = 0) = exp

f 2 2n , 2 Л x • cos 90 + y

R2

exp(( x • sin90).

(7)

В дальнейшем, не уменьшая общности, примем амплитуду падающей волны в выражении (7) равной единице, а временной множитель опустим.

Воспользуемся соотношением (7) для вывода выражения потенциала, создаваемого рассеивателями в верхнем полупространстве, рис. 2. Точку наблюдения М1

рассеянного поля зададим радиус-вектором г1, а положение элементарного участка рассеивающей поверхности определим радиус-вектором р = р(х', у ',0 ).

Для описания потенциала, созданного элементарным участком ёх ' ёу' рассеивающей поверхности, воспользуемся соотношением:

ф(х' у') к н

ёщ = . . • ехр(-/£0 /1 - \\) • п( х' , у') 'ёу ' ;

И

(8)

Ф(х ', у ') = У, ■ (х ', у ',0),

где п(х' , у' )- концентрация рассеивателей в плоскости 2 = 0; У, - коэффициент рассеяния одиночного рассеивателя. Функция двух переменных п( х' , у' ) в настоящей работе не конкретизируется и считается заданной. Тогда выражение (8) можно представить в виде:

ёу = ■

У,

• ехр

(х') ^ + (у') + ,к„ х'^ - «о| К0

п(х' , у' ) ёх'ёу '.

Рис. 2.

(9)

К расчету рассеянного поля

С учетом (9) рассеянное неоднородностями поле в точке наблюдения М1 определяем интегрированием по плоскости х'Оу ':

у=у Я

п( х ', у ')

И

• ехр

(х' )2<Х»'в° + (у )2 + До х^пво - ,*оИ-\) Ко

ёх 'ёу '

(10)

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

Выражение (10) имеет место, когда взаимодействием рассеивателей можно пренебречь. Определим условия, при которых выполнимо данное предположение. С этой целью воспользуемся методом самосогласованного поля [4]. Представим рассеянное поле у в виде суперпозиции поля плоской падающей волны и полей, обусловленных каждым из N рассеивателей, расположенных на рассеивающей поверхности с учетом их взаимодействия:

у(/1) = ехр(-1кгх) + ^ А

ехР(-'к0 И - г„ |)

(11)

п=1

где коэффициенты Ап , учитывающие взаимное влияние рассеивающих

зо

И1 - Гп

неоднородностей, подлежат определению. Если в выражении (11) перенести в левую часть слагаемое, соответствующее к-му рассеивателю, то получим поле, создаваемое всеми остальными рассеивателями совместно с падающей волной в точке наблюдения М1:

I-) л ехР(-ко и - гк |) ( кг) , ехр(-/10 \гх - г„^ ¥(гх) - лк--¡^п-Ь = ехр(-/ки1) +1 Ап--^^-L

п=1 пфк

(12)

Переместим точку наблюдения М1 в точку расположения к -го рассеивателя,

устремив в выражении (12) величину г1 к значению тк. Тогда в правой части (12)

получим поле, создаваемое всеми рассеивателями и плоской волной в точке, где расположен рассеиватель с номером к , исключая сам к -й рассеиватель.

Если теперь выражение для поля, создаваемого всеми рассеивателями кроме к -го, в точке, где расположен к -й рассеиватель, умножить на коэффициент рассеяния то получим амплитуду поля рассеянного к -м рассеивателем:

V

( г- и ехР((\'к- иI)

ехР (кик ) +1 Ап-^-:--

п=1

пфк

И - Гп

■■Ак; к = 1...М.

(13)

Соотношения (13) представляют собой систему уравнений для определения коэффициентов Ап. Введем обозначение:

ехР (-ко \Ик - Ип|)

И - и

к п

ё (ик , Ип ) ёкп , ёкп ё пк

тогда систему для определения коэффициентов рассеяния можно записать в виде:

N / - \

А-I VsAnёkn = V ехр(~1кгк); к = 1...^

п=1 п фк

(14)

Выражение (14) представляет собой систему линейных алгебраических уравнений: ВА=Н

с симметричной матрицей коэффициентов:

1 -^ё12 -^ё13 .. -КчёШ

^?ё21 1 -^ё23 .. - ^ё2 N

В

_ ^ёт VsёN2

и столбцом правых частей:

1

Н

У, • ехр(-/1 • гх) у • ехР(-'к • И2)

У • ехР(-.£ • 4)_

Рассмотрим два очевидных предельных случая - больших и малых расстояний между отдельными рассеивающими включениями. В пределе больших расстояний между неоднородностями | Ип - Ит| ^ з все недиагональные компоненты матрицы В

при фиксированных У3 стремятся к нулевым значениям и В ^ Е, где Е - единичная матрица, а для матрицы Н имеет место предельное выражение:

Н ^ Н =

У, • ехр(-/1 • у • ехр(-.к • И2)

у • ехр(-^ • ^)_ Поэтому выражения для всех коэффициентов Ак примут вид:

Ак = У, • ехр(-/1 • тк ). (15)

Наличие фазового множителя в (15) связано с тем фактом, что центр рассеивающего элемента поверхности находится не в начале координат.

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

При рассмотрении второго предельного случая, когда расстояния между рассеивателями малы и можно положить |И - гк\ ^ 0, необходимо одновременно учитывать асимптотическое поведение зависимости (г0) при малых значениях и0, где И0 - физический размер неоднородности. В частности, если при этом выполняется условие:

^0 У (г0(И,И) = 0,

и ^0

(16)

то матрица коэффициентов В также стремится к единичной матрице (В ^ Е ), для Н имеет место предельное выражение:

Н ^ Н 2

У, • ехр(-/1 • гх) У, • ехр(-/£ • И)

У, • ехр(-/£ • И) и все искомые Ак примут вид: Ак = У, • ехр(-/£ • гх) .

Из выражения (16) следует, что для отражателей, рассеивающая способность которых с уменьшением радиуса ослабевает не медленнее чем И0(1+а), где а > 0,

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

Учитывая полученные ограничения на применение (10), выполним анализ характеристик рассеяния в дальнем поле рассеивающей области. Запишем модуль разности радиус-векторов в декартовых координатах:

Щ-р\=ч1(х -х')2+(у - у')2+2 2,

тогда выражение (10) примет вид:

(17)

V-

и-

( Г (х'2ОС82 в, + у'2) к Г--)—-—)

п (х , у )• ехр -—-^ - кх 81ив0 - (х - х ) + (у - у ) + ;

д/(х - х' )2 +(у - у )2 + ;

(18)

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

И >> ртах будут выполнены неравенства [5, 6]: (хх ' + УУ ') = (?1, ^ .Рты ^ 1

Г1

Г1

Г1

Г1

и асимптотическая формула:

= и,

х ')2 + (у - у ')2 + 22 =4х2 - 2хх ' + (х ')2 + у2 - 2уу ' + (у ')2 + 22 =

2 хх ' + 2 уу ' Р2 ^

1

х2 + у2 + 22 _ 2хх' + 2уу ' . (х')2 +(у')2 =

И

И

И

1 --

■+:г " И

1-

хх ' + уу ' ^

'1

с учетом, которой потенциал рассеянной волны в дальней области поля равен:

у- V, • \\п(х ', у ' )

ехр

(х ') С082 в0 + (у ')2

—0

- ¡к 0 х ' б1п в0 - ¡к 0—

^ - хх ' + уу 'Л

V

И

1 У

( хх ' + уу 'Л 1 ~

ч

йх' йу '.

Пренебрегая членом (хх' + уу')/в амплитудном множителе подынтегрального выражения, после очевидных упрощений получаем:

2

И

V exp(- ik0 r) 7f ( ' , ) V = V —1-- • JJn(x , У )

x exp

t . ____r

(x')2cos2 0о+(y' )2 - ik o x ' sin 0 0 + ik о xx:+yy '

R2

r

dx' dy'

(19)

oo

Выражение (19) дает значения потенциала в дальней области поля, созданного рассеивающей поверхностью, где для каждой точки выполняется условие р «1.

Для удобства анализа рассеянного поля в выражении (19) перейдем к сферической системе координат:

( J \ o 2п ( 2 2 г 2 п , 2 ■ 2 t

( лч т. exp(- ik о Г) ее , ,ч I р COS ф • COS в о + р Sin ф

у(г1,ф,в) = Vs--^J J п(р,ф )• expj---+

м о о I ло (20)

+ ikоp[cos ф' • sin во + sin в • cos(ф - ф ')] }• р dрdф '.

Рассмотрим потенциал рассеянного поля в плоскости падения звукового пучка (плоскости, проходящей через ось пучка и координатную ось Oz). Плоскость падения определяется двумя значениями угла ф : ф = о и ф = к. В результате выражение (2о) принимает вид:

( J \ o 2п Г 2 2 г 2 Л , 2 ■ 2 г

exp(- ik(ír1) f f ( ' ) | р cos ф • cos во + р sin ф +

у(,ф,в) = Vs • -Цп(р,ф ' )• expj

r1 о о I Ro (21)

+ ikQр cos ф '• (sin во ± sin в)}• рdрdф',

где знаки «+» и «-» в фазовом множителе соответствуют случаям ф = о и ф = к.

В соотношении (21) множитель VS является комплексной величиной и содержит фазу коэффициента рассеяния, величина exp (- ik о r1 )/ r1 представляет собой пространственный набег фазы, экспонента в подынтегральном выражении описывает третье фазовое слагаемое - интерференционная фазовая добавка. При проектировании системы распознавания на основе акустической жесткости рассеивателей, в анализе формально нуждаются три перечисленные характеристики рассеянного сигнала. Несмотря на это, два из трех указанных факторов в данной задаче не представляют интереса. Пространственный набег фазы не будет нести полезной информации, поскольку представляет собой шум, поэтому необходимо его подавление. Меры борьбы с пространственным набегом фазы изложены в монографиях [7, 8]. Интерференционная фаза может быть скорректирована за счет изменения расстояния между источником сигнала и исследуемой областью.

ЗАКЛЮЧЕНИЕ

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

единичного рассеивателя на амплитудные и фазовые характеристики рассеянного сигнала.

Из выражения (21) видно, что единственным фазовым множителем является величина:

Фаза в показателе экспоненты обращается в нуль, если:

1). - при любом из двух значений азимутальных углов (ф - 0 или ф - ж) имеем в - 0 (обратное рассеяние по нормали к рассеивающей поверхности);

2). - при ф - ж имеем в - в0 (направление зеркального отражения).

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

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

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

ЛИТЕРАТУРА

1. Андреева И. Б. Рассеяние звука в океанических звукорассеивающих слоях // В кн.: Акустика океана / под ред. Бреховских Л. М. - М.: Наука, 1974. - 300 с.

2. Бреховских Л. М., Лысанов Ю. П. Теоретические основы акустики океана. - Л.: Гидрометеоиздат, 1982. - 264 с.

3. Клей К., Медвин Г. Акустическая океанография. - М.: Мир, 1980. - 580 с.

4. Чабан И. А. Метод самосогласованного поля в применении к расчету эффективных параметров микронеоднородных сред. // Акуст. журн., 1964, Т. 10, №3, с.351 - 358.

5. Басс Ф. Г., Фукс И. М. Рассеяние волн на статистически неровной поверхности. -М: Наука, 1972. - 275 с.

6. Лепендин Л. Ф. Акустика. - М.: Высшая школа, 1978. - 448 с.

7. Гаврилов А. М. Фазозависимые процессы нелинейной акустики: модулированные волны. - Таганрог: ТТИ ЮФУ, 2010. - 352 с.

8. Гаврилов А. М., Савицкий О. А. Фазозависимые взаимодействия акустических волн. - Таганрог: ТТИ ЮФУ, 2010. - 363 с.