Файл: Мясников, Л. Л. Новые методы измерений в подводной акустике и радиотехнике.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 17.10.2024

Просмотров: 74

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Квадрат комплексного волнового числа равен

k2 = k l ^ + i ^ y

(3.12>

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

Z = |0, О, Zz),

причем z — это вертикальная ось; Zz — функция координат г и г. В дальнейшем будем ее обозначать просто через Z.

Решение уравнения (ЗЛО) должно удовлетворять следующему

условию.

Если ввести радиус-вектор /?, имеющий длину R =

= | / r 2 +

z2, то в точке О решение должно обращаться в бесконеч­

ность определенным образом (как HR), а на бесконечности должно обращаться в нуль. Величины, относящиеся к воздушному полу­ пространству (т. е. при z > 0), будем отмечать индексом нуль. Таким образом, k 0 будет волновым числом для воздуха (оно веще­ ственно и равно ю/с). Тогда решение уравнения (3.10) будет иметь вид

2k2

pik0R

при

z >

0;

— ,2

Zo

r + K

 

 

 

(3.13)

2

eikR

 

 

2k\

 

 

 

0

Z'

при

z <

0.

Z ~ k 2 + kl

R

 

 

 

Функции Zo и Z' должны удовлетворять уравнению Гельмгольца для амплитуд. Будем искать выражения для этих функций в виде интеграла Зоммерфельда

Zo = J /о (Я) J0 (Яг) е ~ ^ ~ гdX]

°.

(314

Z' =f { XJ) J o ( X r ) e V v - v

о

где Jo (Яг) — функция Бесселя первого рода нулевого порядка. Данные выражения удовлетворяют дифференциальному уравне­

нию (3.10) и требуемым условиям,

если f 0(X) и f(X) убывают по край­

ней мере как

ИХ2 при

X —>оо.

 

 

Желательно выразить целиком Z0 и Z через интегралы. Это вполне

возможно,

так

как

из

теории

бесселевых

функций следует:

 

 

 

03

Л ---- уо (Яг)

1 Г 2

2

 

J

Г

v х ~feo g dX при г > 0;

R

О

V x 2- k l

 

 

(3.15)

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

4 r e ik* =

Г

,

Я

JolXrt eVv^ *dX при г < 0 .

R

 

J

Y x 2 — k2 ov

 

 

о

55


Если

2 = 0,

то

R = г.

 

 

 

 

В качестве предельных ставятся следующие условия:

1) Zokl =

Zkl

при 2 =

0;

 

 

 

2)

производная функции Z по z должна быть непрерывной.

В

результате

получается

 

 

 

 

 

 

?

X

(lr)e

- V

x 2 - k 2 Z

 

 

Z0= 2k2J

 

1

k° z dx-

(ЗЛ6)

Z = 2kl j- \ j , { X r ) e V v - k * z d K

о

где N = k2 ] / V — ko k20j / V — k2.

Выполняя операции над вектором Герца по формулам (3.11), можно найти компоненты напряженности магнитного поля Н и напряженности электрического поля Е.

Приближенные формулы для интегралов приводят к так на­ зываемой формуле Вейля — Ван-дер-Поля [37].

Вектор Герца для воздушной среды равен

 

 

 

 

 

 

 

 

 

Jk„R

 

 

 

 

 

 

 

 

 

Z 0 =

A _ _ s ,

(3.17)

где 5 — коэффициент

ослабления,

определяемый выражением

 

 

 

 

 

 

 

 

 

 

х + т

 

 

 

 

5 = 2 —4хе~<и+т>2 |

da,

(3.18)

в

котором

 

 

 

 

 

 

 

 

—/со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

------т 711

 

 

 

 

 

 

 

_

R

е

4

 

 

] / '

к2 sin2 О ;

 

 

 

 

 

4

 

 

 

 

У Ч

 

 

sin

ft

(3.19)

 

 

 

 

 

 

 

 

 

3 .

 

 

 

 

 

 

 

 

 

 

Г Т Ш

 

 

 

 

 

 

 

 

 

 

^ —5- cosft.

 

 

 

 

 

 

 

 

 

 

sin й

 

 

 

Обычно

формулу

(3.18)

пишут

в

упрощенном виде,

полагая

 

 

 

 

 

 

 

 

 

 

^2

 

 

•sin 0 = 1 и

пренебрегая

 

величиной ~

под радикалом.

 

 

Обозначив

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_

 

v 2

__

k 0R

/ k 0 \2

 

 

 

 

Р

 

 

 

 

 

 

 

и

имея в виду,

что т ^ О

(так как

c o sft^ O ), получим

 

 

 

 

 

 

2 — 4 У ре -р

Ур

 

 

 

 

5 =

j

ea2da.

(3.20)

— /со

Величина р носит название «числового расстояния».

56


На рис. 3.2 показана зависимость коэффициента ослабления 5 от «числового расстояния» р. Начиная с некоторых достаточна больших значений р коэффициент ослабления выражается в логариф­

мическом масштабе прямой линией.

s' 2

 

S/2J5

«Числовое

расстояние»

 

р

с

учетом

 

проводимости

воды,

когда

можно поло-

 

 

 

жить ' 2 < < i A

r

 

можно

 

 

приравнять

 

 

 

 

 

 

 

 

 

 

 

 

 

вещественному числу

 

 

 

 

 

 

 

 

 

я R

 

 

 

 

 

 

 

 

(3.21)

от -

 

 

Р =---------Т

 

 

 

 

 

 

V

*

+ (

2яе0/

 

 

 

 

0,0010,01 0,1

1,0 10 too

woo р

Если полагать

морскую

 

воду идеаль­

 

Рис. 3.2. Зависимость коэф­

ным проводником, то коэффициент 5 бу­

дет постоянным

и

S -----

2.

 

 

 

 

фициента ослабления от «чис­

Если

к +

т велико

по

сравнению с

лового расстояния»

для воз­

 

духа.

 

единицей,

то на том

участке,

где

подын­

 

 

 

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

cos г1) —

k

l

----sm“

 

____

cos {1

A

l /

i _ J ^

k

V

k2 sin2 d

Здесь — коэффициент отражения Френеля для магнитных ком­ понент.

Случай скользящего падения волны (# близко к л/2) требует применения формулы Вейля — Ван-дер-Поля.

§ 3.2. ИССЛЕДОВАНИЕ РАСПРОСТРАНЕНИЯ СВЕРХДЛИННЫХ РАДИОВОЛН ПОД ВОДОЙ

Рассмотрим теперь распространение электромагнитной волны, излучаемой диполем на границе раздела воздух—вода, в морской воде. В этой области г < 0 , и, для того чтобы с абсолютным возра­ станием z поле не увеличивалось, выбраны те знаки, которые при­ ведены во второй из формул (3.16). Преобразуя, имеем

Z = 2kl

и 0 (\r) е^§**-** zdX__

(3.22)

Вектор Герца для среды, представляющей собой морскую воду, может быть записан в виде

Z =

i k R

(3.23)

S.

5Т


Коэффициент ослабления S также выражается формулой Вейля— Ван-дер-Поля

S = 2 — 4xe-(*+T)‘ и+Jт ea2da,

— * со

но параметры х и т будут другими: з .

------ --------------Г - Я 1

 

___________________

kR

е 4

k

Г .

k2

sin2 ф

'

sin d

' k0 ] /

~~

*2

S1

 

 

3

.

(3.24)

 

 

ти

 

I /

4

 

 

г ‘

 

 

V

2 '

sin -O' COS Ф.

 

Коэффициент ослабления

5

для

морской

воды гораздо меньше

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

множителю eikR). Ослабление

амплитуды в е раз

достигается при

толщине скин-слоя б, определяемой формулой

 

 

_

 

S =

2 у,

(3.25)

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

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

 

 

Zr = Ur,

(3.26)

то

U удовлетворяет уравнению

 

 

 

A U + k2U = 0.

 

Решение этого

уравнения записывается в виде

 

 

 

ikR

(3.27)

 

 

U ^ ^ - S ,

где

R — У г2 +

а2— 2arcos'&, причем а есть радиус

Земли.

 

Когда величина ka велика и велико абсолютное значение диэлек­

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

58


переходит в выражение S = 1 + £, где £ — коэффициент Френеля. Однако для очень длинных волн таким приближением ограничиться нельзя.

В. А. Фоком задача была решена с учетом сферичности Земли. В этом случае функция Герца равна

 

 

 

 

p ika■0'

У-

9)-

(3.28)

 

 

 

 

 

где

 

 

 

 

 

 

 

 

V(x,

 

.

Я

 

elx1w (t — у) dt.

 

 

у, q) = e

 

 

 

 

 

 

 

 

w ' (t ) — qw (t )

 

Здесь интегрирование производится

по

замкнутому

контуру С

в комплексной плоскости t,

охватывающему корни знаменателя

в подынтегральном выражении; х и

у — безразмерные

величины,

через которые г

и О выражаются следующим образом:

 

 

 

 

r= = a{ l + T v ) ’‘

* =

(3-29>

,

 

1

 

 

 

 

( ka\ 3

; Я— параметр, зависящий от кривизны земной

причем А =

(-тг)

 

поверхности в данной точке.

Функция w (/), введенная В. А. Фоком, выражается через функ­

цию Ганкеля первого рода порядка 1/3 по формуле

 

w(t) = e 3

(3.30)

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

Ограничим задачу тем кругом условий, которые имеют отношение к измерениям в области низкочастотных электромагнитных колеба­ ний. Для стандартной звуковой частоты 1000 Гц длина электро­ магнитной волны в воздухе равна 300 км, и, таким образом, можноограничиться рассмотрением поля на расстоянии, малом по сравне­ нию с длиной электромагнитной волны в воздухе или сравнимом с ней, т. е. ограничиться ближним полем. Специальные вопросы дальнего приема к указанной области измерений не относятся. Следовательно, имеются основания пренебречь сферичностью Земли и рассматривать распространение радиоволн над поверхностью моря как над плоскостью. При этом, учитывая скользящее движение волны, нужно исходить из теории Зоммерфельда.

Благодаря значительной электропроводности морской воды и относительно большой величине модуля комплексной диэлектри­

59