К ВОПРОСУ РАСЧЕТА МАКСИМАЛЬНОГО РАСХОДА СЕЛЕВЫХ ПАВОДКОВ ПО "КОСВЕННОМУ МЕТОДУ" (НА ПРИМЕРЕ ВОДОТОКОВ ФЕРГАНСКОЙ ДОЛИНЕ РЕСПУБЛИКИ УЗБЕКИСТАН) Текст научной статьи по специальности «Строительство и архитектура»

UDC 627.51

TO THE QUESTION OF CALCULATION OF THE MAXIMUM FLOW RATE BY THE "INDIRECT METHOD"(ON THE EXAMPLE OF WATERCOURSES IN THE FERGANA VALLEY OF THE REPUBLIC OF UZBEKISTAN)

Abdukakhkhar TULYAGANOV, Candidate of Technical Sciences, Docent

Bekhzodjon MAKHKAMOV*, Leading Specialist of LLC "Road engineering bureau",

Tashkent State University of Transport

1, Temiryulchilar str., Tashkent, 100167, Uzbekistan

Tel. +99893 515-68-25

*E-mail: bekhzodzhon@inbox.ru

Abstract. The article presents the empirical dependences included in the computational formula for determining the maximum discharge of mudflow floods: the coefficient of selenium density and the discharge of the liquid component. It is recommended to calculate the flow rate of the liquid component according to the "limiting intensity" formula based on the "vertical pattern".

Key words: mudflow, maximum discharge, selenium content coefficient; maximum 10-minute intensity, channel slope, flow rate

УУК 627.51

СЕЛ ТОШЦИНЛАРИНИ МАКСИМАЛ САРФИНИ ^ИСОБЛАШДА "БИЛВОСИТА УСУЛИНИ" ЦУЛЛАШГА ДОИР (УЗБЕКИСТОН РЕСПУБЛИКАСИ ФАРГОНА ВОДИЙСИНИНГ КИЧИК ДАРЁЛАРИ МИСОЛИДА)

Абдукаххар ТУЛЯГАНОВ, техника фанлари номзоди, доцент

Бехзоджон МАХКАМОВ*, "Йул Лойиха Бюроси" МЧЖ етакчи мутахассиси

Тошкент давлат транспорт университети

100167. Узбекистан, Тошкент, Темирйулчилар, 1

Тел. +99893 515-68-25

*E-mail: bekhzodzhon@inbox.ru

Аннотация. Маколада сел тошкинларининг максимал чикиндисини аниклаш учун хисоблаш формуласига киритилган эмпирик богликликлар келтирилган: селен зичлиги коеффитсиенти ва суюк компонентнинг чикиши. Суюк компонентнинг оким тезлигини "вертикал накш" асосида "чеклаш интенсивлиги" формуласи буйича хисоблаш тавсия этилади.

Калит сузлар: сел, максимал чикинди сув, селен таркибининг коеффитсиенти; максимал 10 дакикалик интенсивлик, канал мойиллиги, оким тезлиги.

УДК 627.51

К ВОПРОСУ РАСЧЕТА МАКСИМАЛЬНОГО РАСХОДА СЕЛЕВЫХ ПАВОДКОВ ПО «КОСВЕННОМУ МЕТОДУ» (НА ПРИМЕРЕ ВОДОТОКОВ ФЕРГАНСКОЙ ДОЛИНЫ РЕСПУБЛИКИ УЗБЕКИСТАН)

Абдукаххар ТУЛЯГАНОВ, кандидат технических наук, доцент

Бехзоджон МАХКАМОВ*, ведущий специалист ООО «Бюро проектирования дорог»,

Ташкентский Государственный Транспортный университет,

100167, Узбекистан, Ташкент, ул. Темирйулчилар, 1

Тел. +99893 515-68-25

*E-mail: bekhzodzhon@inbox.ru

Аннотация. В статье приведены эмпирические зависимости, входящие в вычислительную формулу определения максимального расхода селевых паводков: коэффициент селеностности и расхода жидкой составляющей. Расход жидкой составляющий рекомендуется вычислять по формуле «предельной интенсивности» основанной на «вертикальной закономерности».

Ключевые слова: селевой поток, максимальный расход, коэффициент селеносности; максимальный 10-минутная интенсивность, уклон русла, скорость потока

1. ВВЕДЕНИЕ

В селеопасных странах, в том числе в Узбекистане ущерб от прохождения селевых потоков в отдельные годы весьма ощутим. Естественно, ущерб от прохождения селей может быть снижен и её вредное воздействие сведено к минимуму [6]. Для правильной научно-обоснованной оценки селеопасности учет основных факторов влияющих на формирование селевых паводков и определение величины селевых максимумов приобретает научно-практический интерес.

2. ЦЕЛЬ И ЗАДАЧИ РАБОТЫ

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

Усовершенствованию подлежит расчет максимального расхода селевых паводков рто «косвенным методом» предложный в 1947 году Д.Л. Соколовским [12], представляющий собой водный расход Отв умноженный на коэффициент селеностности у.

Ото = у Отв . (1)

Как видно из соотношения (1), точность вычисления максимального расхода селевого паводка во многом зависит от величины коэффициента селеностности у.

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

3. МЕТОДИКА ОПРЕДЕЛЕНИЯ РАСЧЕТ КОЭФФИЦИЕНТА СЕЛЕНОСНОСТИ у

Рассмотрим применение соотношения (1) при определении максимального расхода воды на малых предгорных реках и саях Ферганской долины Узбекистана.

В 1978 году С.М. Флешманом [13] на основе статистических обработок большого количества данных и по наивысшим селевым расходам в самых различных природных условиях, получено примерное значение у изменяющаяся от 2-3 до 13-18 и более величин. Как показал анализ и расчет, эти значения коэффициента селеностности у в условиях для малых предгорных рек Узбекистана, оказались завышенными. Это обстоятельство может быть объяснено особенностью как рельефа местности, так и физико-географическими, геологическими факторами. Селевой поток, проходящий на изучаемой территории, относится к несвязанным (несвязанным селевым потоком называют количество воды, которое велико по отношению к количеству грунтовых частиц) и по типу относится к разряду грязевых или грязе-щебенистых потоках [7]. Диаметр наносов не превышает 30-35 мм. При этом селевой паводок сопровождается обилием твердых материалов, величина твердой составляющей из общего объёма селевого потока колеблется от 7-8 до 30%. Прохождение селей отмечается даже при незначительном слое дождя (4 и более мм) с высокой интенсивностью. Характерным также является наличие связи модуля смыва с интенсивностью этого дождя.

Для предварительной оценки значения коэффициента селеностности по исследуемой территории, воспользуемся расчетной зависимостью Тер-Минасьяна (Армения):

У =-Ро-1--(2)

Ро +Е5-Рс (1+е5)

где: ро - плотность селевых выносов;

рс - плотность селевого потока;

5 - пористность грунта;

£ - коэффициент водонасыщенности грунтов;

Расчетная величина уполученная для водотоков Ферганской долины Узбекистана оказалась равным у=1,4. Установленная величина коэффициента селеносности у характеризует в целом весь район Ферганской долины, поскольку параметры, включенные в зависимость (2) были взяты средними по исследуемой территории. Таким образом, для получения у в отдельных селевых водосборах необходимы трудоёмкие полевые материалы: о плотности селевых выносов селевого потока, пористности грунта и коэффициента водонасыщенности грунтов. В связи с этим, была поставлена цель получения достаточно простых по сравнению с другими методами, не требующих больших затрат и дополнительных исследований, в то же время, не уступающих по точности определения принятым в практике проектирования.

4. ПОЛУЧЕННЫЕ РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ

Для решения этой задачи использовались архивные полевые материалы Узгипроводхоза (ныне АО «Узсувлойиха») по селевым паводкам на водотоках предгорных рек и материалы подадырных саев Ферганской долины [7].

Анализ и обработка полевых материалов. Представленных в таблице1 показало характерное для этого района уменьшение значения у от 1,14-1,16 до 1,36 в зависимости от площади водосбора. Такая закономерность была также отмечена в работе [13].

Таблица 1

Сведения о материалах, принятые в анализ предгорных рек

_Ферганской долины_

Количество Количество Величины Площадь Длина русла, Уклон

водотоков зафиксированных расходов, водосборов, км русла

расходов м3/с км2

10 16 14-174 6-169 5,5-40 0,006-0,035

Сравнение полученных величин коэффициента селеносности с материалами Флешмана показало (табл.2), что величина коэффициента селеностности в малых предгорных реках Ферганской долины ниже до 4х кратного значения, чем данные, представленные в монографии. Это может быть объяснено особенностями, отмеченными выше в формировании селевого паводка изучаемого района.

Таблица 2

Сравнение значения коэффициента селеносности у_

Русловая сеть и характер движения селей Площадь бассейна, 2 км2 Чисто эрозионные механизмы формирования селей (склоновая и русловая эрозия) , коэффициент у

Водосборы с одним селевым руслом 10 4-5

по[13] 10-50 3-4

50 2-3

Водосборы с одним селевым руслом в предгорных саях и реках Ферганской долины меньше 10 10 10-50 50 и более 1,30-1,40 1,25 1,20 1,16-1,18

Для выявления общих закономерностей изменения коэффициента селеностности от основных факторов, было предпринято получение связи с морфометрическими характеристиками водосбора, в частности, площадью водосбора и уклона русла [10] (рис.1). При выборе этих характеристик, как определяющие факторы, влияющие на вынос мелкозема, основывались, что площадь водосбора и уклон русла являются интегральными показателями рельефа, характеризующими условия формирования селевых паводков и смыва. Иными словами, считали, что величина площади и уклон русла косвенно определяют длину пути перемещения продуктов смыва и влияют на размеры мелкозема. Наше предположение о связи коэффициента селеностности, как видно из рис.1, оправдано. Однако эта связь оказалась слабой. Коэффициент достоверности связи соответственно для площади водосбора и уклоном русла составил R2=0,52 и R2=0,44. При этом с ростом площади водосбора отмечалось уменьшение значения коэффициента селеносности и рост с увеличением уклона русла.

Й и £ Б

■е- I

<П О) О ^ 5С 5 и

1,3 1,25 1,2 1,15 1,1 1,05 1

100

200

Площадь водосбора, км 2

0

х

01 ^

<П

о

5С

1,3 1,25 1,2 1,15 1,1 1,05 1

0,01 0,02 0,03 Уклон русла

0,04

0

Рис. 1.Зависимость коэффициента селеностности от площади водосбора и уклона русла

Далее анализировались совместное влияние площади водосбора и уклона русла на коэффициент селеностности. Для этой цели воспользовались графическим методомсвязи у=f Эта зависимость (рис.2) оказалась достаточно тесной (коэффициент достоверности равен R2=0,70), где с ростом значения F•J обнаружено уменьшение коэффициента у.

1,3 1,25

1- \ ♦

01 и ЧТ

и I 1- 1,2 1,15 ♦

(П 0 1 01 е; 01 и * ♦

ас #

1,1

0 1 2 3 4 Р • J

Рис.2. Зависимость коэффициента селеностности от площади водосбора и уклона русла (Ц)

Уравнение полученной зависимости имеет вид:

_ 1,22

Т (р^|)0,04" (3)

Для оценки этой расчетной зависимости было сопоставлено фактическое и расчетное значения коэффициента селеностности, не вошедшие в обработку при выводе уравнения (3). В табл.3 представлены результаты сопоставления фактических и расчетных величин по семи водотокам Ферганской долины.Результаты сравнения показали хорошую корреляцию при этом отклонение фактических и расчетных данных, это составило с = ± 3%.

5. РАСЧЕТ ДОЖДЕВОЙ СОСТАВЛЯЮЩЕЙ МАКСИМАЛЬНОГО РАСХОДА СЕЛЕВОГО ПАВОДКА

В литературе известны ряд зависимостей вычисления дождевых максимумов [1,2,4,5,6]. Для расчета Отв при определении максимального расхода дождевого паводка принята предложенная одним из автора формула, основанная на «методе предельной интенсивности» выведенная для условий малых водосборах предгорий Узбекистана [9]:

Отах = 16,67 Я10 к F а ф. (4)

Здесь а10 - интенсивность ливня 10 минутной продолжительности, мм/мин; к— коэффициент перехода от иненсивности 10 минутной продолжительности к интенсивности ливня расчетной продолжительности; F -площадь водосбора;а -коэффициент стока;ф- коэффициент редукции.

Таблица 3

Сравнение фактических и вычисленных по уравнению (4) значений коэффициента селеностности по

некоторым водотокам Ферганской долины

п/н Водоток F, ¥фак. ¥выч. Разность

км2 ¥выч.-¥фак. в %

1 Чартаксай 169 0,022 1,18 1,16 0,02 - 1,69

2 Шишакесай 50 0,013 1,20 1,24 + 0,04 +3,33

3 Сай Бешбуз -№13 5,55 0,053 1,29 1,28 -0,01 -0,78

4 Сай Бешбуз -№13 5,55 0,053 1,30 1,28 +0,02 +1,54

5 Сай Бешбуз -№2 2,50 0,040 1,36 1,34 -0,02 -1,47

6 Сай Бешбуз -№33 1,30 0,088 1,36 1,33 -0,03 -2,26

7 Сай Бешбуз -№36 1,02 0,065 1,30 1,36 -0,06 -4,62

Среднее квадратическое отклонение равно с ~ 3 % П р и м е ч а н и е: F-площадь водозабра; J- уклон русло; ¥фак.и ¥ выч.- соответственно фактический и вычисленный значения коэффициента селеностности; Во второй графе цифрами указаны номера саев

адыра с запада на северо-восток [ 6 ].

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

L = К F п , (5)

где К коэффициент, для предьгоря Узбекистана равен К=2,39; п - показатель степени, п = 0,53. Изложим расчет параметров входящих формулу (4).

А.Максимальный 10 минутный интенсивность дождя ^10).. Максимальный 10-минутный интенсивность дождя определяется по выражению

яю = т! (6)

Здесь т-переходной коэффициент от средней к максимальной интенсивности дождя, принимается равным т = 4.

При отстутствии плювиографов среднюю интенсивность дождя находят по известной формуле:

Д

1 = ^ > (7)

где Д- метеорологическая сила дождя; Т-продолжительность дождя; п - показатель степени. Изменчивые параметры в формуле (7) А, п и Т определяют по эмпирическим формулам полученные для территории предгорья Узбекистана подчинаяющиеся закону "вертикальной зональности" т.е изменению этих величин с ростом высоты местности. Приведем расчетные формулы и ее корреляции (Я): Метеорологическая сила дождя

. 496,7

Д = , Я = 0,66 (8)

ъ0,70 '

где Ъ - высота водосбора, м.

Показатель степени п и продолжительность дождя Т

п = 0,52 А °'15 ,

Я = 0,57

(9)

Тд = Агть

Я = 0,91

(10)

где Аг - единица измерения времени;

Ь - параметр; В среднем Аг =15,1 и Ь= 2,01; г - относительная влажность воздуха.

г = 0, 584 Ъ 099.. Я=0,82 (11)

Б.Пререходной коэффициент 10 минутной интенсивности к интенсивности расчетной продолжительности ( к). Этот параметр тесно связана со скоростью потока (у) и длины основного русла (Ь). Величина коэффицента определяют по формуле:

*0,66

10 V.

к К —

Структура формулы расчета этой скорости принята в виде [8]:

V = кД™ км/мин, где к - коэффициент связанный шероховатости русла к = 0,22; J - уклон русла;

(12)

(13)

т - показатель степени, т = 0,201. В.Коэффициент стока (а). Согласно наших исследование [7,8] коэффициент стока а принимается в зависимости от интенсивности дождя, так для средней интенсивности:

Ьр< 0,5 мм/мин а = 0,5 и ьр > 0,5 мм/мин а = 0,7. Г.Коэффициентредукции (ф).Коэффициент редукции (ф) определяется по известной формуле:

1

ф

— (14)

( Р+С ) пг , ( )

где С - поправочный коэффициент к площади Б, менее резкое убывание максимальных модулей дождевого стока в зоне малых водосборных площадей;

пг - показатель стпени редукции максимальных модулей стока от изменения площади водосбора. Для предгорных рек Узбекистана Пг = 0,35 .

По данным [5] значение С для площадей 5 км2 может быть принять от 1 до 2,30. На основании изложенного метода по формуле (4) рассчитаны нормы модулей максимальных расходов по 21 водотокам расположенных в различных гидрологических районах предгорья Узбекистана. В табл.4 представлены материалы сопоставления между расчетных и фактических нормы модулей максимальных расходов по некоторым рекам Чирчик-Ахангаранского гидрологичесго района Узбекистана.

Таблица 4

п/н Река (сай) Площадь водосбора, км2 Длина, км Уклон, %

1. Каранкульсай 15,6 7,2 120

2. Акташсай 19,3 6,0 130

3. Паркентсай 39,7 14 66

4. Альтинбелсай 19,1 8,2 111

Таблица 4.6

Сравнение рассчитанных и фактических величин модулей максимальных расходов воды по некоторым рекам _Чирчик -Ахангаранского гидрологического района_

Река- пост Б, км2 q т м3/с- км2 q2m м3/с- км2 Разница рассчитенных и фактических величин

м3/с- км2 %

Каранкульсай-Каранкуль 15,6 1,60 1,16 -0,44 -27,5

Акташсай-Акташ 19,3 0,88 0,95 + 0,07 +7,95

Паркентсай-Киргиз 39,7 0,82 0, 69 - 0,13 -15,8

Альтинбельсай- Киргиз 19,1 1,06 1,10 +0,04 +3,77

Примечание: F-площадь водосбора; qm - фактическая норма модуль максимального расхода; q2m-рассчитанные нормы модуля максимальных расходов воды по формуле (4).

Сопоставление между расчетными и фактическими величинами нормы максимального расхода в целом по водотокам Узбекистана показало хорошую сходимость при этом среднеквадратическое отклонение составил с = 27,6%.

Если учесть точность определения фактической величины измеренных максимальных расходов воды с=±18% (по оценке экспертов), то погрешность определения максимального расхода воды по предлагаемой методе будет равна:

а = /а

2ыч. - °фак. = V27,62 - 18,02 = ±18,0%

фак.

Здесь отметим, чтобы рассчитать максимальные дождовые расходы необходимой обеспеченности, следует определить их коэффициенты вариации Сущ и асимметрии С™. Коэффициенты вариации максимальных расходов для рассматриваемых рек территории Узбекистана оказались хорошо связанными с высотами водосборов:

„ 1,89

Су = — . Я=0,71. (15)

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

Кят = Сsm / Сут

или

Сsm = КятСут. (16)

Согласно исследованиям рек предгорья Узбекистана величина К8т изменяется от 1,50 до 3,77 [8].

6. ВЫВОДЫ

5. Усовершенствован метод расчёта максимального расхода селевого паводка «косвенным методом» по предложенной эмпирической зависимости, зависящей от площади водосбора и уклона селевого русла.,

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

7. Выявлена «вертикальная закономерность» интенсивности дождя, присущая для исследуемой территории.

8. Точность определения селевого максимума по предлагаемой методике не привышает с = 20%., что вполне допустимы для практического использования.

ИСПОЛЬЗОВАННАЯ ЛИТЕРАТУРА

1. Алексеев Г.А. Расчет вероятных максимальных расходов воды и обьемов стока снеговых и дождевых паводков.// Труды ГГИ.- 1963. -Вып.39(92).- с. 3-65 [In Russian: Alekseev G.A. Calculation of the probable maximum water discharges and volumes of runoff of snow and rain floods. GGI Proceedings, 1963].

2. Волчек А.А., Шелест Т.А. Моделирование гидрографов дождевых паводков рек Белорусии при отсутствии данных гидрометрических наблюдений // Вестник БГТУ, -2013.-№2. -с3-7 [In Russian: Volchek A.A., Shelest T.A. Modeling of hydrographs of rain floods of rivers of Belarus in the absence of hydrometric observation data. BSTU Bulletin, 2013].

3. Денисов В.М. О расчете максимальных расходов воды дождевых поводков // Тр. САНИИГМИ. 1986. - вып. 119. - с. 23 - 43 [In Russian: Denisov V.M. About calculating the maximum water consumption of rain leashes.

SANIIGMI Proceedings, 1986].

4. Денисов Ю.М. Общая структура формулы расчета максимальных расходов селевых паводков // Труды НИГМИ.-2007.-Вып.4(249).-с.3-14 [In Russian: Denisov Yu.M. General structure of the formula for calculating the maximum discharge of mudflow floods. Proceedings of NIGMI, 2007].

5. Перевозников Б.Ф.Расчет максимального стока при проектировании дорожных сооружений. М., Транспорт, 1975.- 304 с [In Russian: Perevoznikov BF Calculation of the maximum flow in the design of road structures. Moscow: Transport, 1975].

6. Рождественский А.В., Лобанова Т.А.Современная проблема инженерных гидрологических расчетов по обобщению гидрологической информации в России //Метеорология и гидрология.-2011.-№7.-с.82-86 [In Russian: Rozhdestvensky A.V., Lobanova T.A. The modern problem of engineering hydrological calculations for the generalization of hydrological information in Russia. Meteorology and Hydrology, 2011].

7. Салимова Б.Д. Метод расчёта максимальных расходов дождевых вод с малых водосборов.-Т.:- 2011.- 96 с

[In Russian: Salimova B.D. Method for calculating the maximum discharge of rainwater from small catchments. Tashkent, 2011].

8. Тулаганов А.Х.,Тулаганов С.Х. Селлар гидрологияси ва сел о^имларидан автомобиль йулларини химоялаш. -Тошкент: "Ивдисод- Молия". -2014. - 94 бет [In Uzbek: Tulaganov A.X., Tulaganov S.X. Flood hydrology and protection of highways from flood currents. Tashkent: Economy - Finance, 2014].

9. Туляганов А.Х, Салимова Б.Д. Автомобиль йулларидаги сув утказувчи иншоотларни лойихалашда сув ва сел тош^инларининг тавсифларини хисоблаш. -Тошкент, "Iqtisod-Moliya", 2016.-156 бет [In Uzbek: Tulyaganov A.X, Salimova B.D. Calculation of flood and flood characteristics in the design of drainage structures on highways. Tashkent: Economy - Finance, 2016].

10. Туляганов А.Х. Кичик тоголди дарёлари хавзаларида ёмгир тош^инларининг максимал сарфини хисоблаш масаласи/ -Тошкент: Узбекистан география жамияти ахбороти, 56-жилд, 2019. -250 -254 бет [In Uzbek: Tulyaganov A.X. The problem of calculating the maximum consumption of rainwater in the basins of small mountain rivers. Tashkent: Information of the Geographical Society of Uzbekistan , 2019].

11. Туляганов А.Х., Махкамов Б.Р., Мамедова К.З. Значение и расчет коэффициента селеносности в предгорных малых реках Ферганской долины// Universum: Технические науки: электрон.научн.журн. 2020. № 6(75). - с.2-7 [In Russian: Tulyaganov A.Kh., Makhkamov B.R., Mamedova K.Z. The value and calculation of the coefficient of selenium content in the foothill small rivers of the Fergana Valley. Universum: Engineering Sciences: Electronic Science Journal, 2020].

12. Туляганов А.Х., Махкамов Б.Р. К вопросу определения уровня высокой воды по следам паводков на опорах существующих мостов (на примере малых предгорных рек Узбекистана/Штуугеит: тeхнических наук. 2020, №9 (78). С.5-8 [In Russian: Tulyaganov A.Kh., Makhkamov B.R. On the issue of determining the level of high water in the wake of floods on the supports of existing bridges (on the example of small foothill rivers of Uzbekistan.

Universum: technical sciences, 2020].

13. Чеботарев А.И. Общая гидрология. - Л.: Гидрометеоиздат, 1975. - 544 с [In Russian: Chebotarev A.I. General hydrology. Leningrad, Hydrometeoizdat, 1975].

14. Флейшман С.М. Сели.-Изд.2-в, перераб, и доп. - Л.:Гидрометеоиздат, 1978. -310 с [In Russian: Fleishman S.M. They sat down-Publishing house 2-in, revised, and additional. Leningrad, Hydrometeoizdat, 1978].