
книги из ГПНТБ / Соколов, О. А. Видимость под водой
.pdfПолная облученность находится путем интегрирования выраже ния (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): |
л |
|
|
2к |
|
|
|
Е^{г) = Е(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) |
4г |
|
где Еп — разность облученностей двух сторон площадки 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 |
||
|
|||
4я |
где p = cos0. |
|
|
равной |
|
||
/ р 5 (Ѳ, |
ф ) Л о ’ |
|
|
4я |
|
|
Анализируя экспериментальные данные Тайлера по измерению тела распределения яркости в водной толще, Карелин установил, что ход изменения коэффициента углового распределения излуче ния в неоднородной среде сильно зависит от изменения концентра ции взвеси с глубиной. В то же время для показателей вертикаль ного ослабления эта зависимость незначительна [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