Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Соколов, О. А. Видимость под водой

.pdf
Скачиваний:
41
Добавлен:
20.10.2023
Размер:
8.3 Mб
Скачать

Полная облученность находится путем интегрирования выраже­ ния (1.39) в пределах телесного угла 2я:

£(Ѳп,

<рп,

z )=

J В(Ѳ',

cp', z) cos Ѳ'сіш.

(1.40)

 

 

 

2т:

 

 

При оп ~ л (схему

отсчета

углов

см. на рис. 1.8) выражение

(1.40) преобразуется

в

выражение

облученности горизонтальной

Рис. 1.8. Схема отсчета зенитных (Ѳ) и азиму­ тальных (ср) углов.

плоскости, создаваемой нисходящим излучением, т. е. облученно­ сти сверху

2-

т>/2

 

Е^{г)= Е(ъ, z )= j

j В (0, cp, z)cosBdBd<?.

(1.41)

9=0 0=0

 

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

* /2

E^(z)= E (k, z)= 2 k J В (В, z) cos Ѳsin ѲdB.

(1-42)

о

 

Облученность горизонтальной плоскости, создаваемая восходя­ щим излучением, или облученность снизу (Ѳп =0), выразится анало­

гично (1.41) и (1.42):

л

 

 

 

 

Е^{г) = Е(0, z ) = j“

j

В {В, ср, z)cosBdBd<?,

(1.43)

9=0 0 =я/2

 

 

іс

 

 

Ej(z) = E(0, z)= 2n

I

В (В, z) cos Ѳsin ВdB.

(1.44)

 

it/2

 

 

27

2) Пространственная облученность является суммой нормальных облученностей элементарных площадок, имеющих общий центр в исследуемой точке пространства и освещаемых перпендикулярно падающими световыми лучами:

Во (z) = J В (Q, cp, z)dw.

(1.45)

Для осесимметричного светового поля можно записать

Eü(z )= 2k I В(0, z)sin0 dO.

(1-46)

о

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

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

 

2-

іс/2

 

 

El(z) =

j

j Я (0, «5», z)dOd<?,

(1.47)

 

9=0 0 =0

 

 

для глубинного режима

 

-/2

 

 

 

 

 

 

ЕІ (z)=2icj В {В,

z)da.

(1.48)

 

 

о

 

 

Сумма нормальных

облученностей

от излучения,

приходящего

в точку из нижнего полупространства, называется полупространст­ венной облученностью снизу.

2“

-

 

 

Eq(z)= J

j

В (В, <р, z)dBdy,

(1-49)

<р=0 0=-/2

 

для глубинного режима

 

 

 

El (z ) = 2 k

j В (О, z) da.

(1.50)

 

 

*/2

 

3) Средняя сферическая облученность численно равна

средней

плотности светового потока на поверхности сферы исчезающе ма­ лого диаметра с центром в исследуемой точке пространства:

Е4к(г)= 0,25 [ В (0, <р, z) da,

(1.51)

4п

ИЛИ

E ^(z)= 0,25E0(z).

28

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

Для глубинного режима средняя сферическая облу­ ченность равна

(г)=0,5* J В (0, г) sin 0 dO.

о(1.52)

Величину средней сфери­ ческой облученности легко

измерить, применяя К пририс I g Перенос энергии излучения через емнику энергии излучения произвольно ориентированную площадку,

матированные сферические насадки.

4) Объемная плотность энергии излучения равна пределу отно­ шения количества этой энергии к объему, в котором она заключена, при стремлении объема к нулю:

J B d (o = ^ -E 0= -^- Д4п,

(1.53)

4г.

 

где с — скорость света в среде.

Доказательство справедливости соотношения (1.53) читатель может найти в светотехнической литературе [1,2, 50].

5) Вектор переноса энергии излучения (световой вектор) опреде­ ляет значение и направление переноса потока излучения в иссле­ дуемой точке поля через единицу площади, перпендикулярной на­ правлению переноса.

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

Рассмотрим сначала произвольный объем в непоглощающей и нерассеивающей среде, заполненной движущимся в разных направ­ лениях излучением. Внутри этого объема выберем произвольно ори­ ентированную площадку dS (рис. 1.9) с нормалью п в точке Р. Пусть через площадку dS в направлении, определяемом зенитным и азимутальным углами 0 и ср, внутри телесного угла dm проходит поток излучения і 2Ф. Тогда яркость излучения в направлении (0, ср)

29

можно выразить

В (0, с?)

(РФ

(1.54)

dS cos du>

где £ n —угол между направлениями нормали к площадке dS и на­

правлением

(Ѳ, cp); dS cos t,'n— площадь проекции площадки dS на

плоскость, перпендикулярную направлению (0, ф).

 

Из (1.54)

следует

 

 

гі2Ф =5(0, tp) dS cos Cn du.

(1.55)

Поток излучения d2Ф заключен в пределах элементарного телесного угла rfco. Для получения полного потока, проходящего через пло­ щадку dS, необходимо выражение (1.55) проинтегрировать по всем направлениям сферы:

d<è= dS I В (Ѳ, cp)cosCn du.

(1.56)

4тс

Рассмотрим (1.56) внимательно.

При С„< cos£^ >0, подынтегральное выражение имеет по­ ложительный знак.

При cos Сп <0, подынтегральное выражение имеет от­

рицательный знак.

Следовательно, d<$>— есть разность потоков излучения, прохо­ дящих через площадку dS с двух ее сторон.

Выражение (1.56) можно переписать в виде

4 г = Вп = j Я (Ѳ, <р) cos Cn rfw,

(1.57)

 

где Еп — разность облученностей двух сторон площадки dS.

Равенство (1.57) можно интерпретировать как проекцию неко­ торого вектора на нормаль п к площадке dS. Тогда сам вектор мо­ жет быть записан так:

Е = J В (6, ср) d<a.

(1.58)

4Іт

 

Проекции вектора на оси прямоугольной системы координат мо­ жно представить в виде

ЕХ = ^В{Ъ,

tp) cos С* бЦ

 

Ву= j В (Ѳ,

cp) cos Су da>,

 

Дг= |Д (Ѳ ,

cp)cosCz du>,

(1.59)

где t,x, t,y и tz — углы, которые составляет вектор с осями координат.

30

Вектор Е и является вектором переноса энергии излучения или,

как его называл Гершун, световым вектором [50].

 

Так как направление

(Ѳ, <р) было выбрано произвольно, можно

утверждать, что вектор

переноса энергии излучения

существует

в любой точке поля, он не зависит от выбора системы

координат,

а его проекция на какое-либо направление численно равна разности

облученностей двух сторон

элементарной

площадки, помещенной

в исследуемую точку поля

перпендикулярно

этому направлению.

Если площадка расположена перпендикулярно

направлению по­

тока излучения, то разность облученностей

двух ее сторон имеет

максимальную величину, количественно

равную

модулю вектора

потока излучения:

 

 

 

 

 

 

 

 

(1.60)

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

6) Дивергенция (расхождение) вектора переноса энергии излу чения определяет объемную плотность поглощения или излучения светового потока в исследуемой точке поля.

Рассмотрим произвольный объем dV в поглощающей и рассеи­ вающей среде. Если вектор переноса энергии излучения проинте­ грировать по замкнутой поверхности, ограничивающей объем dV, то получим некоторую величину dO, равную

(1.61)

s V

где dSv — векторный элемент поверхности S, ограничивающей ис­ следуемый объем dV.

Очевидно, что величина dO является разностью потоков излу­ чения входящего внутрь объема dV через поверхность dS и выхо­ дящего из него. Так как водная среда сама не излучает энергии в оптическом диапазоне волн, то разность между входящим и выхо­ дящим потоками имеет положительный знак, т. е. гіФ>0. Эта раз­ ность представляет собой поток, поглощенный в объеме dV. Если величину dO отнести к величине объема dV, то получим

(1.62)

гіФ где ——— объемная плотность поглощения потока излучения.

dV

31

Из закона сохранения энергии следует

 

f/ф

у.Е0= —х J В (0, tp)afo),

(1.63)

ІѴ

4s

 

 

 

где Е0— пространственная облученность.

 

Следовательно,

d iv E = —vlE0.

(1.64)

 

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

дЕх

дЕу

дЕг

(1.65)

div Е = дх

~ду~

dz

Для плоскопараллельной среды, каковой практически является толща моря,

дЕх _ дЕу

дх ду

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

 

div Е (г)=- дЕдр{

=

- у.Е0,

(1.66)

или

 

 

 

 

div Е (z) = —у. j

5(0,

<р, z)ckо,

(1.67)

 

4s

 

 

 

где

E(z) — вектор переноса энергии излучения на глубине г.

 

 

Для глубинного режима дивергенция вектора равна

 

'

div E (z )= —2w. I В(ф, z) sin 0 dB.

(1.68)

о

 

 

 

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

£, (2)

 

*<г> = - г щ г -

(1-69)

8) Коэффициент яркости моря определяется отношением ярко­ сти рассеянного излучения, идущего из толщи моря непосредст­ венно под его поверхностью в каком-либо направлении, к яркости

32

идеально белой изотропно рассеивающей поверхности, при освеще­ нии последней естественным излучением:

р(0, ср)=

В (0, <р,

0)

(0, о, 0)

(1.70)

 

Вп

~

Еп

 

где В (Ѳ, ер, 0) — яркость излучения, распространяющегося под водой

в направлении

(0, ф) при глубине z->- 0; В0— яркость идеально бе­

лой изотропно

рассеивающей

 

поверхности; Еп— облученность по­

верхности моря.

 

 

 

 

 

Зависимость коэффициента яркости моря от длины волны излу- ,

чения характеризует истинный цвет моря.

 

нисходящего или во

9)

Коэффициент углового распределения

сходящего излучения выражает отношение

полупространственной

облученности сверху или снизу к

облученности горизонтальной

плоскости сверху или снизу:

 

 

 

 

 

 

 

 

 

 

 

 

(L71)

 

 

gi(z)

 

Ео (*)

 

(1.72)

 

 

 

Ef (*)

 

 

 

 

 

 

 

Рассмотрим выражение величин

(г)

и

g^ (z ), записанные

в следующем виде:

 

 

 

 

 

 

 

* / 2

В (Ѳ, z) sin Ѳd&

 

 

 

 

J

 

 

 

 

ё ^ ) = - Ж ------------------------

(1.73)

 

 

J

В (0, z) sin 0 cos 0 d0

 

 

 

о

 

 

 

 

 

 

 

 

J

В (0,

z) sin 0 db

 

 

 

 

(*)

 

--------------------- •

(1-74)

 

 

J

В (0, z) sin 0 cos 0 d0

 

 

 

r./2

 

 

 

 

 

Если все излучение падает на слой водной среды параллельным

пучком перпендикулярно слою,

т. е. в направлении 0= 0 или Ѳ= я,

то cos 0= 1 и

(z) = 1.

 

 

 

 

 

Если же излучение изотропно, т. е. его яркость равна для всех направлений полусферы, величину В (0, z) можно вынести за знак интеграла, как постоянную. Нетрудно убедиться, что в этом случае gi (г) =2.

Таким образом, величина g (z) характеризует неравномерность

И

в

пределах полусферы.

распределения яркости по направлениям

Коэффициент углового распределения

g

в литературе часто

It

называют величиной, обратной среднему косинусу углов наклона

3 Заказ № 601

33

световых пучков

нисходящего (восходящего) излучения

1

м-U

S В (0,

cp)cfco

 

где p = cos0.

 

равной

 

/ р 5 (Ѳ,

ф ) Л о ’

 

 

 

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

Для гидрооптических расчетов исключительно важное значение имеет определение первичных гндрооптических характеристик по

параметрам светового

поля.

Обоснование расчета

этим

методом

показателя поглощения к было дано Гершуном [52].

 

 

 

Проекция вектора переноса энергии излучения на вертикальную

ось г равна разности облученностей

сверху и снизу

на глубине z:

 

Ez {z)= E i { z ) - E Ji{z).

 

 

 

(1.75)

Тогда дивергенция

этого вектора,

 

согласно

(1.66)

равна

div Е (z) = —£■

(г) -

Е, (г)] =

(z)

dEy (z)

(1.76)

dz

 

dz

 

 

 

 

 

 

Ослабление облученности сверху и снизу при

изменении глу­

бины на величину dz можно выразить:

 

 

 

 

 

 

d E ^ z)= -E ^ (z)a .^ z)d z,

 

 

 

(1.77)

dE^ (z )= —Ef (z) af (2) dz,

 

 

 

(1.78)

где а 1 (2) и at (2) — показатели вертикального

ослабления

сверху

и снизу.

 

 

 

 

 

 

 

 

Из чего следует:

 

 

 

 

 

 

 

 

, V

_ J __

 

dEi <*)

 

 

 

(1.79)

М г)= _

(*)

 

dz

 

 

 

 

 

 

 

 

°Ч

£+(2)

 

dE) (2)

 

 

 

(1.80)

 

dz

 

 

 

 

 

 

 

 

Для глубинного режима

«+(^ )= at (z)=a,

тогда, учитывая (1.79) и (1.80), дивергенцию светового вектора можно выразить следующим образом:

div Е ( г ) = —а (2) — (2)] = —хЕ0(2).

(1.81)

34

Следовательно, показатель поглощения равен

a [ E , { z ) - E , { z ) \

(1.82)

Ео (*)

или, вводя сферическую облученность, можно записать:

4* [£; ( z ) - £ +(z)]

(1.83)

* (*>

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

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

Согласно (1.76) и (1.69),

 

div E (z)= £^ (z )

dE, (z)

 

 

dE. (z)

 

 

 

 

( z ) dz

^ '

' £ | ( г ) dz

 

 

 

 

 

 

 

 

 

 

 

= -

£* (2) [“, (2) -7 ? (г) af (2)],

 

 

(1.84)

где R (z)

Bi (z)

коэффициент диффузного отражения моря.

Ei (z)

 

можно преобразовать

 

 

 

 

 

 

Величину at (z)

 

 

 

 

(г)

 

 

1

dE^ (z)

d \R ( г ) £

; ( г ) I

£+ (г)

 

 

dz

 

at ( 2 ) = - Ef ( г )

dz

R (z)

 

(^ )

dz

R

( г )

£ * ( г )

 

 

7? И

rfg; (г)

7?

 

 

 

 

 

 

 

 

 

dz

 

 

 

% (2),

 

 

(1.85)

 

 

R ( * ) £ *

( г )

Л (г)

 

 

 

 

 

 

 

 

 

 

где /? — производная 7? по глубине (г).

 

в

(1.84),

получим

 

Подставляя полученное значение at (z)

 

 

 

 

 

 

 

Ь /

 

 

 

 

divE(z) = - £

| (z)at (z)

! - £ ( * ) -

 

R{z)

= x (z )£ 0(z),

(1.86)

 

 

(*)

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

к (г)

 

E^ (2Г)

 

 

 

 

x (z)= a; (z) 1 -7 ? (z)

 

 

 

(1.87)

 

“j. (?)

 

£o (2)

 

 

 

 

 

 

 

 

 

 

 

Формула Пелевина пригодна для определения показателя погло­

щения при любых условиях

освещения

поверхности

моря и при

наличии неоднородности в его толще.

 

 

 

 

 

 

 

Связь между показателем вертикального ослабления в глубин­

ном режиме и показателем

поглощения

можно получить также

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

(z)

и(z). Согласно (1.81)

a [£* (z) - Ef (z)] = x [ £ o (z)+£<! (г)] .

( 1.88)

3*

35

Разделив обе части равенства (1.88) на Е\ (z) и решая относи­ тельно а, получим

S* (?) + £+ (2) % (2)

X.

(1.89)

l - R ( 2 )

 

 

При установившемся световом режиме в однородной среде ве­ личины g\, gj и R не зависят от глубины. Так как величина R мала (для районов открытого моря R ^0,02), то можно записать

 

 

(1.90)

Заметим, что в формуле (1.87) множитель

( г )

1

Ео (г)

, где

 

'((г)

 

 

/

gl (z)— разновидность коэффициента углового распределения ни­ сходящего потока излучения, по абсолютному значению близ­ кая к g .

1.4. Спектральный состав излучения в море

При определении характеристик светового поля в водной среде в предыдущем разделе была сделана оговорка, что все ха­ рактеристики рассматриваются для монохроматического излуче­ ния. Но в глубину моря проникает излучение, имеющее сложный спектральный состав, зависящий как от спектрального состава па­ дающего на поверхность моря излучения, так и от избирательных свойств поглощения и рассеяния, присущих среде. При этом состав излучения, проникающего в море, непрерывно изменяется с глуби­ ной [23, 48, 53—63].

Поглощающие и рассеивающие свойства в различных районах Мирового океана различны. Спектральный состав падающего из­ лучения также может меняться в зависимости от атмосферных ус­ ловий и времени дня. Поэтому точно предсказать состав излучения в каком-либо районе моря на какой-либо глубине трудно. Это воз­ можно лишь при условии точного знания оптических характеристик среды и спектрального распределения энергии в падающем излуче­ нии в данный момент времени.

Таблица 1.2

X нм

1

2

3

X нм

1

 

О

3

400

0 ,4 4 7

0 ,6 8

1 ,8 0

560

1 ,0 0 0

1 ,0 0

1 ,0 0

420

0 ,6 7 2

1 ,0 0

1 ,7 3

580

1 ,0 0 0

0 ,9 5

0 ,8 8

440

0 ,7 8 5

1 ,0 3

1 ,7 2

600

1 ,0 1 3

0 ,9 0

0 ,7 7

4 6 0

0 ,9 2 0

1 ,1 7

1 ,7 5

620

1 ,0 0 3

0 ,8 7

0 ,7 0

480

0 ,9 9 7

1 ,1 9

1 ,6 1

640

1 ,0 0 0

0 ,9 2

0 ,6 2

500

1 ,0 2 3

1 ,1 5

1 ,3 8

660

1 ,00

1

0 ,8 8

0 ,5 7

520

1 ,0 0 2

1 ,0 6

1 ,2 2

680

0 ,9 9 5

0 ,8 8

0 ,5 1

540

0 ,9 9 0

1 ,0 3

1 ,1 2

700

0 ,9 6 5

0 ,9 2

0 ,4 6

36

Соседние файлы в папке книги из ГПНТБ