CC BY
58
9. Zhuravlev, V.B. Opredelitelj rihb basseyjna Verkhneyj Obi: monografiya / V.B. Zhuravlev, S.L. Lomakin, S.N. Satyukov. - Barnaul, 2010.
10. [Eh/r]. - R/d: http://www.eurolab.ua/encyclopedia/trip/1397
11. Pravdin, I.F. Rukovodstvo po izucheniyu rihb. - M., 1966.
12. Tipovihe metodiki issledovaniya produktivnosti vidov rihb v predelakh ikh arealov. - Viljnyus, 1974. - Ch. 1.
13. Tipovihe metodiki issledovaniya produktivnosti vidov rihb v predelakh ikh arealov. - Viljnyus, 1978. - Ch. 2.
14. Nikoljskiyj, G.V. Ehkologiya rihb. - M., 1963.
15. Metodicheskie ukazaniya po opredeleniyu vozbuditeleyj geljmintozoonozov v presnovodnihkh rihbakh. - 1999. - № 13-4-2/1751.
16. SanPiN 15-6/44. Sanitarnihe pravila po sanitarno-geljmintologicheskoyj ehkspertize rihbih i usloviyam ee obezzarazhivaniya ot lichinok difillobotriid i opistorkhisov. - M., 1990.
Статья поступила в редакцию 16.11.14
УДК 556.531
Kudishin A.V. CONSTRUCTION OF ONE-DIMENSIONAL HYDRAULIC MODEL FOR A WATER FLOW USING CURVILINEAR COORDINATES. The construction of one-dimensional hydraulic model for a water flow based on the integration of the Saint-Venant plane equations recorded in a curvilinear coordinate system is discussed. While compared with a standard model, the model constructed by the author makes it possible to calculate the river's volume and the surface area together with flood plains of significant width more precisely. The conclusion of the paper has a list of advantages of the produced model: 1) the correct calculation of the volume and areas of the surface of the river and its flood plains; 2) in case of a complicated geometry of the mouth of the river with many of the islands, the researcher expects to obtain precise number of the live mouth of the river and the module of the river discharge; 3) a more precise definition of volume-based characteristics of the flood plain and the correction of the figures that show the location of the water waves in floods.
Key words: hydrology, hydraulics, mathematical modeling, curvilinear coordinates.
А.В. Кудишин, канд. физ.-мат. наук, с.н.с. ИВЭП СО РАН, Барнаул, E-mail: bezmater@iwep.ru
СПОСОБ ПОСТРОЕНИЯ ОДНОМЕРНОЙ ГИДРАВЛИЧЕСКОЙ МОДЕЛИ ТЕЧЕНИЯ С ИСПОЛЬЗОВАНИЕМ КРИВОЛИНЕЙНЫХ КООРДИНАТ*
Рассмотрен способ построения одномерной гидравлической модели течения воды, основанный на интегрировании записанных в криволинейной системе координат плановых уравнений Сен-Венана. С целью разработки некоторой методологической основы построения одномерной гидравлической модели течения представлен альтернативный вариант, основанный на сформулированной в двухмерной криволинейной системе координат гидравлической модели течения. Полученная одномерная гидравлическая модель течения по сравнению со стандартной позволяет более точно вычислять объем и площадь зеркала реки вместе с пойменными участками значительной ширины.
Ключевые слова: гидрология, гидравлика, математическое моделирование, криволинейные координаты.
Для решения многих задач речной гидравлики широко и успешно применяются модели течения, основанные на теории «мелкой воды». В рамках этой теории течение воды описывается системой гиперболических уравнений (уравнения Сен-Венана), получаемых интегрированием по глубине уравнений Эйлера при определенных условиях [1]. Применение конечно-разностных методов для численного решения уравнений «мелкой воды» позволяет с достаточной точностью определять гидравлические параметры течений для решений ряда важных практических задач. Это затопление территорий и прохождение паводков в речной сети, задачи качества воды, ледотермика, расчет процессов руслоформирования.
В практическом плане существенную роль играет размерность используемых уравнений по пространству. Более универсальными являются двухмерные (плановые) уравнения Сен-Ве-нана, которые при интегрировании по поперечной (к основному направлению течения) координате преобразуются к упрощенным одномерным уравнениям Сен-Венана, в которых продольная координата ориентирована вдоль русла.
В практике часто встречаются ситуации с существенно двухмерным течением. Это геометрически сложные извилистые участки русел, многорукавность, наличие поймы, затопление урбанистической территории. Плановая модель в таких случаях проявляет себя хорошо при наличии достаточно подробной цифровой модели рельефа (ЦМР), но при этом необходимо отметить высокую требовательность к вычислительным ресурсам доступных многоядерных систем (4-12
дЪ дъ * ий]
--\--= и , сумма по j=1,2
д^ дх.
ядер). Последнее накладывает существенные ограничения на размер области определения задачи как по пространству (около 1000 км2, шаг 10-20 м), так и по временному расчетному периоду (1-2 месяца). Развитие вычислительных систем уменьшает эти ограничения, но пока не в достаточной степени. Если увеличить область опрелеления задачи на порядок (по пространству и по времени), то плановая модель теряет свою привлекательность как излишне громоздкая, и для этих задач используют уравнения Сен-Венана в одномерном приближении, позволяющие более или менее правильно воспроизводить продольную структуру течения [2]. Одномерные нестационарные модели течения характеризуются высокой скоростью выполнения численных алгоритмов и достаточно высокой точностью определения гидравлических характеристик водотоков. Эти качества обусловливают эффективность использования одномерных моделей при расчетах на больших временных промежутках. Основной недостаток одномерных моделей прямо следует из способа их получения - теряется информация по двухмерной структуре течения. В работе делается попытка методологически уменьшить потери этой информации на стадии построения одномерной модели интегрированием двухмерных уравнений Сен-Венана.
Общепринятый подход построения одномерной гидравлической модели течения (уравнения Сен-Венана) основан на использовании концепции «прямоугольных ортогональных координат». Запишем в декартовой системе двухмерные уравнения Сен-Венана [1]. Уравнение неразрывности:
(1)
Уравнения сохранения импульса (для i=1,2):
dh * Udi dh * udi*u dg „ _
d ■ t----- + g*h*^~ + Fdi = 0 , сумма по j=1,2 (2)
dt dxj dx.
= 2 (3)
c
Здесь h - глубина, t- время, ^ — уровень поверности воды, С - коэффициент Шези, - компоненты вектора скорости в исходных декартовых координатах х1=х и х2=у. Предполагая русло слабоизвилистым, ось х=х расположим! примерно вдоль основного направления течения (рис. 1).
X
Рис. 1. Русло реки в прямоугольных координатах
Береговые границы реки в координатах (и,у) соответствуют Ь1(х) и Ь2(х). Отбросим уравнение сохранения импульса в проекции на ась у. Проинтиг рируем по кнординате у от Ь1(ч) до Ь2(х) оставшиеся уравнения (1)-(2) с учетом дифференцирования по параметру. При этом предполагаем поперечный градиент уровня воды равным нулю и используем кинематическое уеловие на береговой границе ЬЬ,х)=Ь1,2 реки [1]:
дь db _ = 0
dt + ud_(xU=b) dx Ud= (x0 =b) = 0
(4)
Получим стандартную запись уравнений Сен-Венана [2]:
да д( ^
— + — = О да дх
(5)
dQ du*Q DC \Q\Qp
— +-— + g * — * CO + g * CO * = 0
¿5а dx dx K2
(6)
Также имеется эквивалинтный выше приведеннооу способ вывода уравнений (4)-(5), осно ванный на рассмотрени) баланса массы и сил в объеме воды ограниченным русловым отсеком. Извилистость русла не учитывается, в общем случае координата х направлена по руслу и не соответствует исходной декартовой системе координат. Концепция хорошо зарекомендовала себя при решении многих практических задач. Однако для извилистого русла сложной геометрии или при наличии поймы такой подход может приводеть к ошибкам при определении триведенного сопротивления русла, вычислении объемов в расчетнеlх ячейках конечно-разностной сетки и опредвлении приведениого размера ячеок по продольной координ ате.
С целью разоаботки некоторой методологической основы! построения однооерной гидравличеткой модели течения рассмотяен альтернативный вариант, основоннеlй на сформулированной в двухмеиной криволинейной системе координат гидравлической редели течения.
Пусть (хй, X.) - декартовы! координаты!, (¿г1,. ) - система ктиволинейнеlх X^ = X^ (б^, ), I = 1,2 - не зависящее от времени преобразование координат и [3-4]:
координат,
d(ХЬ x2) _ _
J =---—- - якобиан преобразования;
d(e\e2)
gjj = , i, j = 1,2 - ковариантный метрический тензор;
k=1
g = *ßjk, ^ j = 1,2 - контрвариантный метрический тензор;
k=1
a = dxk_ ß = de aki - d i > Pik -
de'
dx.
k
2 2
dl2 (dxk )2 g ij * de 'de J - квадратичная дифференциальная форма.
k=1 i=1 j=1
3]:
При переходе к криволинейной системе координат уравнения (1)-(2) в проекции на ковариантный базис записываются [1
dh 1 dj * h * и] 0
--1--*-= 0 , сумма по j=1,2
dt j de
(7)
dh * и1 1 dJ * h * и1 * U} h j dC Гi h k j Fi 0
-+ —*-+ —— + TL*h*u *uJ + F =0
dt J deJ de }
В уравнении (8) «g» без индексов - ускорение свободного падения, i=1,2 и сумма по
(8)
И
= ^ udk * ßik - контрвариантная компонента скорости
* i i ii
(9)
rkj - символ Кристоффеля 2-рода
F =
g * \U\ * и'
(10)
Будем рассматривать ортогональные системы кооодина^ для которых: g12 = 0 , g12 = 0 , g" = Уg
J = Jgn* g22 (не зависят от времени); символы Кристоффеля 2-рода при этом записываются в общем виде:
Г' =
lkj -
2 * gи
Ogti ^ß ^k]
de dek de
, ¡,k,j=1,2.
(ii)
Г1 = ^
;n_ r1 = _ 1 ' 1 22 _
T--
11 de
1
1
2*g„ de
dg22 г1 =Г1 = 1 , _
1 ' 2 * „ - d„2
11
11
11 Г2
p2 -All —
■>22 Г
12 21 An — x о i —
' 22
2 , d 2 ' 22 2 , d 2 ' M2 21 2 , d 1
2 * g22 de 2 * g22 de 2 * g22 de
/12^,1 2 Координатные линии 15 ,5 I (5 - вдоль русла, 5 - поперек русла) в физическом пространстве (х, у) проецируются
в семейство линий S2) соответственно (рис. 2).
Рис. 2. Построение криволинейной координатной сетки
2 1-1
Отбросим уравнение сохранения импульса (8) в проекции на ось 5 (для i=2). Умножим уравнение (7) на а (8) для 1=1
на ( J * ), получим для ортогональной системы координат:
dJ * h dJ * h * u1 „ .
+-:-= 0, J = 1,2
dt dsJ dJ gn * J * h * u1
dt
- +
_ 0^ _ _
+ Я * Ь *л[я22 + * J * н1 = 0
05
В (12)-(13) сумма по ^=1,2.
Преобразуем слагаемые в квадратных скобках в (13):
I- dJ * h * u1 * uJ 1- r „1 , —f —j
/gii *-—-+ Vgii * J*гз *h*u *u
dsJ
+
(12)
(13)
dJ * h * u1 * u} dJg11 * J * h * U * u} dJg11 —r —
1 *-:-= —--:----—:—* J * h * U * U
dS dS dS
g1
(14)
1 0а 0л/ а
Учитывая --;= *-=- и (11), имеем:
2л а 05 05
dVgn
dsJ
= Tij*4g~, Г/, =г;1; j=1,2
1j ~ i J1
Подставляя (14) в (13) и используя (15),(9), преобразуем (13):
d^Jgn * J * h * u1 d-jg11 * J * h * u1 * u3 dt dS
+ И А и 2
* Р 9" 5и * Р N 522 и1Р18- -и 2 9д
+
сумма по j=1,2.
(16)
+ 5 А И ^^52^ ^ТТ + А 3 А р 1 = 0
98
Обозначим 8 = 8, 7] = 82 . Береговые границы реки (линии уреза воды) в координатах (8,]) соответствуют Т]](ё)
и Проинтегрируем по координате 7] от ]/(£) до Т]^(е) уравнения (12) и (16) с учетом дифференцирования по
параметру. При этом предполагаем поперечный градиент уровня воды равным нулю и используем кинематическое условие на береговой границе 1]12(Ъ 8) реки:
% + и\а,%) р9^-и 28,%) = 0 91 98
После интегрирования по«1^» уравнение неразрывности (12) имеет дивергентный вид:
90 9Q п — + — = 0 9t 9а
В 2
|И а ёВ - площадь поперечного сечения, построенного по кривой S2=S2(x,y).
(17)
(18)
В2
0=| И А
В1
]2
]2
В2
В 2
0 = 13 А И А ё] = ^ л/517 А И А ё] = ^ л/577 А И а ёВ = аР | И а ёВ = а.ю (19)
В1
В1
Здесь ёВ = л/522 а ё] - элемент длины по кривой Б2=Б2(х,у), в соответствии с определением квадратичной дифференциальной формой.
В 2
В = |ёВ - ширина зеркала воды в поперечном сечении, построенном по кривой Б2=Б2(х,у).
1]2 _7 72 и * В2
2 = 13 А И А и А ё] = ^д/577 Ад/522 А И А . 1 А ё] = | И А и. А ёВ
] ] \511 В1
Ц - физический расход воды через поперечное сечение, построенное по кривой S2=S2(x,y). Используем стандартные формы осреднений:
]2 _ В 2 В2
|л/511 а 3 а И а и1 а ё] = 511 а И а и* А ёВ = у1 |И А и* А ёВ =ух А 2 ,
]1
]2
В1
В1
В2
В2
I 11 * * *
IV511 А 3 А И А и А и А ё] = I И А и* А и* А ёВ = ~1 I И А и* А ёВ = у2 А ~1 А 0 ,
]1 В1
В1
и, =
= б/
(20)
Л2
В 2
{Т^/ * 3 * Р * ё1 = {Т5п
В1
5 * и * и,
-1 * ёВ
с
(21)
Ы2
| И И и*
Ы1
и1 И
ду 711
да2
- и2 И
Ч 7 22
да1
: тТы = Ч (а)
(22)
С учетом (20)-(22) после интегрирования по «□» уравнение импульса (16) имеет вид:
дг, И О ду2 И ~ И П дС
———— л— -^^ + 7 и0ИА + '
да да да
В 2
Ч(а) + 7 -л[тп
и, И и, \ 1 И ТВ = 0
(23)
В1
Итоговая система уравнений Сен-Венана в одномерном приближении имеет вид (18)-(23). Уравнение (18) является точным для расчета баланса массы. Уравнение (23) довольно грубо моделирует закон сохранения импульса для реальных течений. Для малых чисел Фруда определяющими являются силы давления и трения [5], это 3-й и 4-й члены в уравнении (23).
В этом случае можно положить близкие к 1 коэффициенты = 1. При и* ^ 0 Ч(а) ^ 0 Ч(а) , в результате
уравнение (23) по своей структуре становится похожим на (6). Это позволяет использовать для численного решения (18)-(23) алгоритмы, применяемые для решения (5)-(6).
Семейство реперных для координаты «£» линий S:L (ограниченное количество, см. рис. 2), строится в транзитной зоне по линиям тока (либо приблизительно по линиям тока). Далее по реперным линиям строится ортогональная система координат. С учетом специфики решаемой задачи требования к построению сетки достаточно слабые, основное значение имеют линии
S2=S2(x,y) (по которым проводится интегрирование, при этом значения g22 не нужны). Допустимо построения сетки по областям (по длине русла, границу областей нужно располагать в местах, где русло близко к призматическому), требования к сгущению по « Ы» отсутствуют. Пример построения сетки для реального русла представлен на рисунке 3. В данном случае для построения сетки использовался алгебраический метод.
С
Рис. 3. Вариант построения координатной сетки для участка р. Обь
Выводы
Предложен способ построения одномерной гидравлической модели течения воды, основанный на интегрировании записанных в криволинейной ортогональной системе координат плановых уравнений Сен-Венана. Полученная одномерная гидравлическая модель течения имеет следующие достоинства по сравнению со стандартной:
- корректное вычисление объемов и площади зеркала реки вместе с пойменными участками значительной ширины;
- в случае сложной геометрии русла (много межостровных
проток) можно ожидать более точного определения живого сечения русла и модуля расхода;
- более точное определение емкостных характеристик поймы, соответственно, уточнение расчетов прохождения волн паводков (половодья).
*Работа выполнена при финансовой поддержке в рамках проекта СО РАН VШ.76.1.1 «Исследование процессов формирования стока и разработка информационно-моделирующих систем оперативного прогнозирования опасных гидрологических ситуаций для крупных речных систем Сибири».
Библиографический список
1. Вольцингер, Н.Е. Длинноволновая динамика прибрежной зоны / Н.Е. Вольцингер, К.А. Клеванный, Е.Н Пелиновский. - Л., 1989.
2. Васильев, О.Ф. Математическое моделирование качества воды в системах открытых русел / О.Ф. Васильев, А.Ф. Воеводин // Динамика сплошной среды. Вып. 2. - Новосибирск, 1975.
3. Лисейкин, В.Д. Технология построения разностных сеток / В.Д. Лисейкин, Ю.И. Шокин, И.А. Васева, Ю.В. Лиханова. - Новосибирск, 2009.
4. Корн, Г. Справочник по математике / Г. Корн, Т. Корн. - М., 1977.
5. Кюнж, Ж.А. Численные методы в задачах речной гидравлики / Ж.А. Кюнж, Ф.М. Холли, А. Вервей. - М., 1985.
Bibliography
1. Voljcinger, N.E. Dlinnovolnovaya dinamika pribrezhnoyj zonih / N.E. Voljcinger, K.A. Klevannihyj, E.N Pelinovskiyj. - L., 1989.
2. Vasiljev, O.F. Matematicheskoe modelirovanie kachestva vodih v sistemakh otkrihtihkh rusel / O.F. Vasiljev, A.F. Voevodin // Dinamika sploshnoyj sredih. Vihp. 2. - Novosibirsk, 1975.
3. Liseyjkin, V.D. Tekhnologiya postroeniya raznostnihkh setok / V.D. Liseyjkin, Yu.I. Shokin, I.A. Vaseva, Yu.V. Likhanova. - Novosibirsk, 2009.
4. Korn, G. Spravochnik po matematike / G. Korn, T. Korn. - M., 1977.
5. Kyunzh, Zh.A. Chislennihe metodih v zadachakh rechnoyj gidravliki / Zh.A. Kyunzh, F.M. Kholli, A. Verveyj. - M., 1985.
Статья поступила в редакцию 16.11.14
УДК 597.5 + 591.639 + 591.5
Popov P.A. ICHTHYOCENOSES OF THE YENISEI RIVER MOUTH. The paper studies a topical question of the study of the structure and functioning of ecosystems in river mouths in theory and in connection with the negative influence of anthropogenic factors upon them. At the mouth section of the Yenisei River in the fish habitat a total number of 32 species of freshwater fishes is found. This is 30 per cent of ichthyofauna of Siberia. The situation of the state of ichthyocenoses in the mouth section of the Yenisei River is formed by the fishes of the family of Coregonidae. The commercial stock is not enough for intensive fishing. The number of Acipenser baerii Brandt, Coregonus autumnalis (Pallas) and Stenodus leucichthys nelma (Gueldenstaedt) decreased significantly as a result of the intensive fishing in the past. The further investigation of the ichthyocenoses in this territory under research is important because of the necessity of working out a strategy of protection, preservation and rational usage of the fish stock in the Yenisei. Key words: Subarctic region of Siberia, the Yenisei River, ichthyocenosis, ecological of fishes.
П. А. Попов, д-р биол. наук, проф., в.н.с, Институт водных и экологических проблем СО РАН, проф. Новосибирского гос. ун-та, г. Новосибирск, E-mail: icg-adm@bionet.nsc.ru
ИХТИОЦЕНОЗЫ УСТЬЕВОЙ ОБЛАСТИ ЕНИСЕЯ
В статье затрагиваются актуальные вопросы изучения структуры и функционирования экосистем устья рек в теоретическом отношении, и в связи с отрицательным воздействием на них антропогенного фактора. В состав ихтиоценозов устьевой области Енисея входит в общей сложности 32 пресноводных вида, или 30% ихтиофауны Сибири. Облик ихтиоценозов формируют рыбы семейства сиговых, промысловые запасы которых существенно подорваны в результате интенсивного вылова. Дальнейшее изучение ихтиоценозов этой территории актуально в связи с необходимостью разработки стратегии охраны и рационального использования рыбных запасов Енисея. Ключевые слова: Субарктика Сибири, Енисей, ихтиоценоз, экология рыб.
Устьевые области Великих сибирских рек - Оби, Енисея и Лены - играют важнейшую роль в глобальных процессах взаимодействия вод суши и океана, в круговороте вещества и энергии в гидросфере северного полушария Земли [1]. Благодаря обилию приносимых рекой биогенных веществ, повышенной турбулентности вод, особенностям солевого режима устьевые области рек являются одними из наиболее биопродуктивных районов на Земле [1-2]. Изучение структуры и функционирования экосистем устья рек актуально как в теоретическом отношении, так и в связи с отрицательным воздействием на них антропогенного фактора. В частности, в бассейн Енисея негативное влияние на гидробионтов связано с зарегулированием рек плотинами ГЭС, повышенным, а на ряде участков Енисея и высоким содержанием в воде загрязняющих веществ [2-3], вселением в водохранилища Енисея леща (Abramis brama), распространившегося к настоящему времени вплоть до Полярного круга, сазана (Cyprinus carpió) и амурской формы серебряного карася (Carassius auratus) [4].
Цель настоящей работы - анализ и обобщение сведений о структуре ихтиоценозов устьевой области Енисея, которая относится к эстуарно-дельтовому типу и состоит из устьевого участка реки (869 км) и устьевого взморья - Енисейского залива (около 350 км) и прилегающей части Карского моря. Речная граница устьевой области располагается в районе впадения в Енисей р. Нижней Тунгуски [1; 5].
Условия обитания гидробионтов, включая рыб, в устьевой области Енисея весьма неоднородны, сложны и определяются особенностями уровенного режима, притоком соленых вод в период преобладания нагонных течений, коротким периодом открытой воды и ее низкими температурами. Минерализация речных вод устьевой области колеблется от 120 мг/л (в половодье) до 230 мг/л (зимой). Кислородный режим в реке благоприятный для гидробионтов в течение всех сезонов года [6-7].
В общей сложности в бассейне Енисея обитает 48 видов рыб [8]. В состав ихтиофауны водоемов устьевой области Енисея входит 32 вида пресноводных и 7 видов морских рыб (табл.). В русле Енисея на отрезке от устья Нижней Тунгуски до дельты и губы включительно зарегистрирован 31 вид, из которых 7 видов являются полупроходными и поднимаются с мест нагула -
южной части залива, губы и дельты - на нерест в Енисей. Это сибирский осетр, омуль, муксун, чир, сибирская ряпушка, нельма и азиатская зубатая корюшка, или зубатка. Остальные виды пресноводных рыб, живущие в русловых водах Нижнего Енисея, являются туводными и больших по протяженности миграций в течение жизни не совершают [9-10].
В левых притоках устьевой области Енисея обитает в общей сложности 29 видов рыб, в т.ч. в Турухане - 24, в Большой Хете -21, в Танаме - 18. По мере продвижения на север в левобережье из состава ихтиоценозов или полностью исчезают (пескарь, караси, язь), или становятся весьма малочисленными (елец, плотва) сравнительно теплолюбивые виды рыб семейства карповых. Из других представителей бореального равнинного комплекса в рассмотренных реках левобережья сравнительно многочисленны щука и ерш. Рыбы семейства сиговых в более северных реках левобережья становятся многочисленнее, и их роль в промысле возрастает [11].
В правобережных притоках устьевой области Енисея отмечено в общей сложности 28 видов рыб, в т.ч. в Подкаменной Тунгуске - 22, в Нижней Тунгуске - 20, в Курейке - 24, в Хан-тайке - 23. Облик ихтиофауны в этих притоках формируют холо-долюбивые виды-реофилы: таймень, ленок, хариус и сиг-валек. Однако лишь хариус сравнительно многочислен здесь, запасы тайменя и ленка на большинстве участков этих рек истощены промыслом, а сиг-валек вообще весьма редок [12].
Видовой состав ихтиофауны Енисейского залива определяется, прежде всего, отношением обитающих здесь рыб к солености, а характер их распределения - динамикой взаимодействия соленой и пресной масс воды. Постоянно встречаются на акватории залива 6 видов морских и 6 видов пресноводных полупроходных рыб. Кроме того из морских рыб в заливе изредка вылавливается полярная камбала. Из пресноводных рыб в дельте и на приустьевых участках притоков преимущественно в период весеннего паводка встречаются в небольшом числе туводные формы сига-пыжьяна и чира, а также сибирский хариус, щука и ледовито-морская рогатка. Из дельты в залив в период половодья или во время продолжительных сгонных ветров заходят единично осетр и налим. Рыбы семейства карповых в заливе не обнаружены [9].
