НА УЧНЫЕ СТА ТЬИ
Гидрометеорология и экология №3 2009
УДК 556.16.013
РАСЧЕТ ЭМПИРИЧЕСКОЙ ОБЕСПЕЧЕННОСТИ МАКСИМАЛЬНЫХ РАСХОДОВ ПАВОДКОВ
Доктор геогр. наук Г.Н. Трофимов
Рассмотрена возможность расчета эмпирической обеспеченности максимальных в году расходов воды с учетом разности соседних членов ранжированного ряда. Показано, что при таком подходе оценка эмпирической обеспеченности является наименее смещенной по сравнению с оценками по формулам, традиционно используемым в гидрологии. Высказано мнение о необходимости включения в расчеты кривыгх обеспеченности максимумов селевого характера.
Расчеты максимальных расходов являются одной из наиболее сложных гидрологических задач при проектировании и строительстве разного рода сооружений на реках. Занижение максимальных расходов приводит к разрушению сооружений, затоплению освоенных прибрежных территорий и, нередко, к человеческим жертвам. В свою очередь, завышение этих максимумов повышает стоимость сооружений и снижает их экономическую эффективность. Необходимо отметить, что проектировщиков интересуют максимумы редкой повторяемости (малой обеспеченности).
В гидрологии традиционно сложились три направления вычисления расчетных гидрологических характеристик:
- расчеты при наличии достаточно длительного ряда наблюдений;
- расчеты при недостаточности данных гидрометрических наблюдений и
- расчеты при отсутствии этих данных.
Не касаясь здесь двух последних направлений расчетов, остановимся на расчетах максимальных расходов воды малой обеспеченности для рек, имеющих достаточно длительные ряды наблюдений. В этом случае для расчетов максимумов заданной малой обеспеченности есть два пути: а) с помощью того или иного теоретического распределения и б) экстраполяцией эмпирической обеспеченности в область редкой повторяемости.
Как известно, статистические расчеты, в частности, расчеты случайной величины той или иной вероятности превышения, основаны, прежде всего, на требованиях не связности и однородности рядов наблюде-
131
ний. В случаях, если в этих рядах обнаруживаются аномальные (экстремальные) значения переменной, то они классифицируются как «выбросы», ряды рекомендуется «выравнивать», а экстремумы исключать из расчетов. Выявление «выбросов» производится на основании критических значений тех или иных критериев согласия.
Отметим предварительно, что проектировщиков в целях безопасности сооружений интересуют именно такие аномальные расходы воды и поэтому в гидрологических расчетах воспользоваться такими рекомендациями (исключение выбросов) не представляется возможным. Однако, остановимся на оценках однородности рядов наблюдений за максимальными расходами воды. Гидрологи, учитывая сказанное выше, рекомендуют проводить расчеты максимумов с учетом их генезиса (снегового или дождевого). Сразу отметим, что на горных реках, помимо упомянутых причин, максимумы могут формироваться при прорыве плотин высокогорных озер, при опорожнении временных водосборников в теле горных ледников и в ледниковых моренах, при разрушении завалов русел рек, образованных при сходе оползней, и т.п. Другими словами, генезис максимальных расходов воды на горных реках может быть самым разным и, что достаточно важно, не всегда можно с уверенностью установить причину их формирования.
Статистические критерии для оценки «выбросов» применяют, преследуя одну из следующих целей:
- выровнять ряды наблюдений перед анализом (отбрасывание выбросов);
- убедиться, что аномальные значения присутствуют, что указывает на необходимость пересмотра процедуры вычислений и получения новых статистических данных;
- выделить наблюдения, которые могут представлять особый интерес именно из-за их экстремальности.
Теория и практические методы отбрасывания выбросов разработаны слабо, что подтверждается следующим высказыванием Е. Гамбела [1]: «Отбрасывание аномальных значений на чисто статистической основе было и остается весьма опасной процедурой. Само их присутствие может являться доказательством того, что исследуемая совокупность в действительности отличается от предполагаемой».
Статистики экстремальных отклонений, содержат, как правило, разность между экстремальным и последующими значениями, или между экстремальным и выборочным средним значением, а также среднее квадратиче-
132
ское отклонение. К первым можно отнести критерий Ирвина [10], построенный на исследовании двух соседних членов ранжированного ряда хт и хт+1
' Х..И
1 = ^
s
(1)
здесь ах - среднеквадратическое отклонение.
Г.Ф. Лакин приводит критерий оценки сомнительных вариантов (потенциальных выбросов) для максимумов в виде [5]:
t
_ Х1 Х2
(2)
где х1 - наибольшее значение варианты, х2 - следующее за ним значение в ранжированном ряде, хп-1 - второе значение переменной в возрастающем ряде, п - число членов ряда. Предположение, что х1 принадлежит данному ряду отбрасывается при 4 > а, где а- выбранный уровень значимости.
А. Анскомб для аналогичной оценки аномальных значений предложил рассматривать разности УМ = ХМ - Х, где ХМ - максимальное значение, а Х- среднее значение переменной. Исходная гипотеза отбрасывается, если по абсолютному значению | УМ\ > е ах [3]. В свою очередь, с рассчитывается по формуле:
п ■ c
■(у—1)
n ■ c
2
у—1
(3)
где t - распределение Стьюдента, v - число степеней свободы.
В данной работе остановимся на расчете обеспеченности максимальных расходов дождевых паводков, формирующихся на малых низкогорных реках Узбекистана. Этот выбор обусловлен, во-первых, однородностью генезиса паводков на этих реках и, во-вторых, многочисленностью таких водотоков. Так по данным [12] в Узбекистане насчитывается 17777 рек, причем основная их часть это малые реки длиной менее 10 км. Исследуя дождевые максимумы, мы использовали данные по 42 рекам, площадь водосбора которых не превышает 100 км2, а средняя высота бассейнов не более 2,0 км. Для всех этих рек нами проведены статистические расчеты, как их основных статистических характеристик максимумов, так и их эмпирических обеспеченностей. Для иллюстрации полученных результатов приведем данные по малой низкогорной реке - Шаугазсай (левый приток р. Ахан-гаран). Площадь водосбора реки до гидропоста (урочище Караташ) - 65,8
133
2 Хп-1
t
a
2
V
V
км2, длина водотока от его истока до гидропоста 15 км, общая протяженность реки 22 км, средняя высота водосбора 1,66 км.
Нами для максимальных расходов воды р. Шаугазсай - пост Кара-таш проведена проверка принадлежности наибольшего за многолетие расхода воды всей совокупности наблюдений по трем перечисленным критериям (табл. 1). По этому посту имеется ряд наблюдений с 1951 по 2004 гг., т.е. длиной несколько более 50 лет. Максимальный, измеренный на посту, расход воды составил 172 м3/с (27 июля 1964 г.), наименьший годовой максимум - 0,55 м3/с, который наблюдался дважды в 1962 и 1977 гг. Средний из максимальных расходов - 7,46 м3/с, а второй член ряда в возрастающем порядке равен 0,57 м3/с (1986 г.). Коэффициент вариации этого ряда равен 3,20, а коэффициент асимметрии 6,61.
Как видно в табл. 1 все эмпирические значения критериев превосходят их теоретические значения при 5 % уровне значимости. Это говорит о том, что в расчетах и параметрах распределения и самой кривой обеспеченности использовать значение максимального расхода воды 1964 года нельзя и ряд наблюдений неоднороден.
Таблица 1
Проверка принадлежности наибольшего расхода воды р. Шаугазсай общей совокупности с помощью критериев согласия
йщах, м /с Год Критерий принадлежности члена ряда совокупности
по Ирвину по Лакину по Анскомбу
1 ^0,95 г г0,95 ссх
172,00 33,50 0,57 7,46 1964 1969 1986 Среднее 5,87 1,10 0,81 0,19 164,5 160,3
Таким образом, складывается противоречивая ситуация - на гидропосту зарегистрирован резко отличный от общей совокупности годовой максимум. Однако, каких-либо исключительных гидрометеорологических (да и других природных) событий в бассейне реки не происходило. Можно лишь предположить, что максимальный расход 27 июля 1964 года сформировался при выпадении довольно интенсивного дождя, хотя на посту Караташ, расположенному на высоте 1,1 км и на близлежащих постах на эту дату осадков не зафиксировано. Добавим, что в это время года снежный покров в бассейне уже стаял.
Как отмечено выше, для вычисления максимумов малой обеспеченности рекомендуется использовать то или иное теоретическое распределение. В гидрологии чаще всего используются распределения К. Пирсона,
134
С.Н. Крицкого и М.Ф. Менкеля, реже распределение - Е. Гамбела (двойной показательный закон, или закон первого типа) и еще реже распределение -Фишера-Типпета I типа (для максимумов) [4, 7-9]. Если учесть статистические характеристики ряда максимумов р. Шаугазсай, то в таблицах Крицкого-Менкеля при С/С = 2 значения модульного коэффициента для С = 3,20 нет, а распределение Пирсона ограничено величиной С = 6,4. Использование распределения Гамбела для максимумов по сути дела сводятся к спрямлению кривой вероятности не превышения и экстраполяции ее в область высоких значений. Как видно на рисунке, спрямления рядов максимальных расходов воды р. Шаугазсай по распределению Гамбела не произошло. Максимумы, рассчитанные по распределению Фишера-Типпета, в районе обеспеченностей выше 50 % получились отрицательными.
Й 0,8
Й 0,4
& 0,2 т
50
100
150 <3т.х, м /с 200
Рис. Расчет кривой обеспеченности максимальных расходов воды р. Шаугазсай (по Гамбелу).
В условиях, когда теоретические распределения к расчету обеспеченности максимальных расходов воды не применимы, для расчета переменных малой обеспеченности можно использовать их эмпирические значения.
Для описания эмпирической функции распределения и экстраполяции данных расчетов в область редкой повторяемости в гидрологии используется ряд формул. Чаще всего это формулы:
т - 0,5
- А. Хазена
- Н.Н. Чегодаева
- С.Н. Крицкого-М.Ф. Менкеля
- Е.Г. Блохинова
Р
*
Р
*
Р
*
Р*
п
т - 0,3 п + 0,4 т
п +1 т - 0,4
п + 0,2
(4)
(5)
(6) (7)
1,0
0,6
0,0
0
135
1
( т 1 ^
- Д. Коудена р =-=--= + - , (8)
л/п +1 ул/п 2 ^
и др. Здесьр - эмпирические вероятности превышения в долях единицы.
Отметим, что формула Крицкого-Менкеля полностью соответствует ранее предложенной формуле Е. Гамбела, а формулу Д. Коудена зачастую называют формулой «тгп/тах». Все приведенные выше формулы не учитывают ни абсолютные изменения исходных переменных, ни их статистические характеристики. Основными аргументами в них являются - порядковый номер переменной т и общая длина ряда п. Детальный анализ таких формул проделан Ю.Б. Виноградовым [1] и М.А. Мамедовым [6]. Последний, отмечая, что допущение, принятое при выводе вышеприведенных формул о том, что вероятность случайной величины принимается независящей от значения самой случайной величины недопустимо в случаях, когда отклонения ее близки, или превышают саму измеряемую величину. Для таких рядов максимальных расходов он предложил формулу обеспеченности максимумов в виде:
р = т к г , (9)
п +1 + Кт
где Кт - модульный коэффициент для т-го члена ранжированного ряда. Показатель степени г рекомендован в зависимости от величины коэффициента вариации ряда в пределах от 2 до 5. Нужно отметить, что формула, предложенная Мамедовым, пожалуй, единственная на сегодняшний день, когда в ее структуру входят статистические характеристики конкретного ряда. Добавим, что, по его мнению, эта формула дает более стабильные величины обеспеченности при рядах разной продолжительности и, самое главное, «позволяет «отскакивающие» величины включить в общую статистическую совокупность, ... что делает экстраполяцию эмпирических кривых более надежной» [6].
К недостаткам формулы (9), видимо, нужно отнести то, что важнейший ее параметр г задается одинаковым при широком изменении коэффициентов вариации рядов наблюдений, к примеру, при изменении Су более чем в 2 раза (0,3.0,7) г принимается равным 4.
Касаясь формул (4). (8), нужно сказать, что все оценки вероятности превышения, вычисленные по ним, являются смещенными в сравнении с вероятностью, вычисляемой по формуле р* = т . Второй, более суп
щественный недостаток этих формул состоит в том, что для любого ре-
136
ального соотношения соседних членов ряда (в том числе крайних его членов) мы получаем вероятности превышения, обусловленные только номером переменной и общей длиной ряда.
Выше, в качестве критерия принадлежности экстремальных членов ряда общей совокупности, был приведен критерий Ирвина. Структура этого критерия, кроме его «прямого назначения», позволяет давать оценку «близости» соседних членов ряда. Это его свойство использовано нами при расчетах эмпирической обеспеченности максимальных расходов воды по формуле:
* тп -Л2
Р =П-. (10)
п(п + Л2)
Прежде чем перейти к анализу полученных оценок вероятности превышения максимумов р. Шаугазсай заметим, что оценка, вычисленная по (10) будет наименее смещенной по сравнению с оценками, полученными по формулам (4)...(9) и при Qm - Qm+l ® 0 мы получаем практически несмещенную оценку вероятности превышения.
Результаты вычисления эмпирических обеспеченностей максимумов р. Шаугазсай по вышеприведенным формулам для верхней части кривой (первые 10 членов ранжированного ряда) приведены в табл. 2.
Как видим (табл. 2), формула (8) дает наиболее завышенные величины эмпирической обеспеченности по сравнению со всеми остальными. Наименьшее значение вероятности превышения для первого члена ряда дает формула М. Мамедова, почти в два раза большая оценка обеспеченности получена по формуле (10). Для второго и последующих членов ряда по формуле (10) получены, как и предполагалось, практически не смещенные оценки эмпирической вероятности превышения.
Как известно, помимо несмещенности эмпирических оценок статистических характеристик необходимо определить их эффективность. Эффективной является статистическая оценка, имеющая, при заданном объеме выборки, наименьшую дисперсию. Как показали расчеты, наименьшей дисперсией обладает оценка эмпирической обеспеченности, полученная по формуле М. Мамедова, но, как видно в табл. 2 эта оценка весьма существенна смещена относительно несмещенной оценки, вычисленной по
формуле р* = — . Несколько большая дисперсия оценок получена при
п
расчете обеспеченности по формуле (10), однако увеличение дисперсии этой оценки не на столько существенно, чтобы предпочесть более смещенную оценку, полученную по формуле (9).
137
Таблица 2
Эмпирические обеспеченности максимальных расходов воды р. Шаугазсай
Эмпирические обеспеченности максимальных расходов воды
Qmдx,
м /с
т п
т т — 0,3 т — 0,5 т — 0,4
п +1 п + 0,4 п п + 0,2
л[п +11л/п
т
: + 0,5
т
п +1 + К*
тп — 1
п\ п +12
172 33,5 24,0 19,3 17,2 8,79 7,58 6,23 5,29 5,09
0,0189 0,0377 0,0566 0,0755 0,0943 0,1132 0,1321 0,1509 0,1698 0,1887
0,0185 0,0370 0,0556 0,0741 0,0926 0,1111 0,1296 0,1481 0,1667 0,1852
0,0131 0,0318 0,0506 0,0693 0,0880 0,1067 0,1255 0,1442 0,1629 0,1816
0,0094 0,0283 0,0472 0,0660 0,0849 0,1038 0,1226 0,1415 0,1604 0,1792
0,0113 0,0301 0,0489 0,0677 0,0865 0,1053 0,1241 0,1429 0,1617 0,1805
0,0771 0,0937 0,1104 0,1270 0,1436 0,1602 0,1768 0,1935 0,2101 0,2267
0,0017 0,0270 0,0466 0,0659 0,0893 0,1083 0,1272 0,1463 0,1651 0,1836
0,0042 0,0376 0,0565 0,0755 0,0941 0,1132 0,1321 0,1509 0,1698 0,1887
1
Выше было сказано, что с учетом основных статистических характеристик рядов наблюдений за максимумами использование теоретических распределений затруднительно, а порой и невозможно. Тем не менее, для сравнения 1 %-х максимумов мы использовали распределение Пирсона, экстраполируя значения нормированных отклонений от среднего до величины С = 6,6.
Также, учитывая тот факт, что при оценке эмпирических обеспе-ченностей по формулам (4)... (8) мы оперируем только числом членов ряда и порядковым номером переменной. Проведено «разреживание» рядов, путем исключения второго, третьего и т.д. членов ряда в его хронологической последовательности и проанализированы, дополнительно к исходному ряду, состоящему из 54 членов, три выборки, состоящие из 26, 18 и 10 членов ряда. Добавим, что при случайном исключении точек наибольший член ряда включался в каждую выборку (табл. 3). Операция разреживания рядов проведена с целью анализа «устойчивости» (или изменяемости) расчетных квантилей при изменении числа членов ряда п.
Экстраполяция эмпирической кривой обеспеченности в область ее малых значений производилась по экспоненте с использованием первых трех-четырех членов убывающего ряда по формуле
0 = Ьп(д) — Ьп( р*)
^р ъ ' У )
здесь р - эмпирическая обеспеченность, Qp - расход воды заданной обеспеченности, д и Ъ - параметры связи.
Таблица 3
138
Максимальные расходы воды 1 % обеспеченности, вычисленные по
различным формулам
Число точек Максимальный расход воды 1 % обеспеченности, м3/с
по Пирсону т п т т - 0,3 п+0,4 т - 0,5 п т - 0,4 1 ( т , 0 51 т тп -1
п +1 п + 0,2 4п+1 ) п+1+кт п^ п + Я2 ^
53 120 222 220 180 152 166 740 91,2 120
26 161 310 306 251 212 232 764 121 138
18 187 355 348 287 244 266 765 142 143
10 240 430 418 348 297 323 770 174 120
При сопоставлении максимумов, рассчитанных по эмпирическим формулам, с их значением, полученным по распределению Пирсона видно, что величина расходов вычисленных по формуле «тгп/тах» сильно завышена. Вполне понятно завышение максимумов, рассчитанных по формулам (4)...(9), т.к. с уменьшением числа членов ряда неизбежно возрастает обеспеченность в области ее малых значений. Иная картина изменения максимумов 1 % обеспеченности, вычисленных по формуле (10). Добавим, что в среднем наименьшее отклонение эмпирических максимумов от их значения по Пирсону получено по формуле (10) - 22,0 %, несколько большее отклонение (25,1 %) получено по формуле (9). Отметим также, что, по всем формулам, за исключением (10), прослеживается увеличение расчетных величин с уменьшением числа членов ряда.
Помимо максимумов, зарегистрированных на гидропостах, в гидрологических расчетах считается необходимым учитывать, так называемые «исторические максимумы» [2]. Очевидно, что включение этих максимумов в статистические ряды с точки зрения математической статистики неоправдано. Однако потребности практики и, в первую очередь, обеспечение безопасности объектов, находящихся в зоне возможного затопления паводковыми водами, определяют необходимость включения в расчеты именно этих «сверх максимальных» расходов воды. Для горных рек - это обычно селевые расходы.
При прохождении паводков селевого характера оценка их максимумов проводится, как правило, на тех участках реки, где имелись разрушения, причиненные потоком, либо происходила существенная деформация русла (изменение планового рисунка русла, обрушение берегов и т.п.). Другими словами, местоположение так называемых «селевых поперечников» часто не совпадает с расположением гидрологических постов. По этой причине для больших рек включение максимумов такого рода в об-
139
щую совокупность не всегда обосновано. Для малых рек, видимо, можно включать данные Каталога селей и других справочников по максимумам в общий ряд, так как в силу небольшой протяженности водотоков, изменение местоположения створа измерений этого максимума, существенного влияния на его размер не окажет.
На основании обследования следов катастрофического паводка, прошедшего 27 июля 1964 г. на р. Шаугазсай, максимальный расход был оценен в 274 м3/с для створа, расположенного несколько ниже гидропоста [11]. Замена максимума, измеренного на гидропосту, этой величиной, естественно, изменила статистические характеристики ряда - среднее значение расхода - 9,38 м3/с, коэффициенты вариации и асимметрии - 4,00 и 7,00 соответственно. Расчеты эмпирических обес-печенностей по всем формулам для нового ряда приведены в табл. 4.
Таблица 4
Эмпирические обеспеченности максимальных расходов воды р. Шаугазсай (с включением селевого максимума 1964 г.)
Эмпирические обеспеченности максимальных расходов воды
Qmax, м3/с
т п
т т - 0,3 т - 0,5 т - 0,4
п+1 п + 0,4 п п + 0,2
л[п +11л/п
т
: + 0,5
п +1 + К*
тп -1
п\ п +12
274 33,5 24,0 19,3 17,2 8,79 7,58 6,23 5,29 5,09
0,0189 0,0377 0,0566 0,0755 0,0943 0,1132 0,1321 0,1509 0,1698 0,1887
0,0185 0,0370 0,0556 0,0741 0,0926 0,1111 0,1296 0,1481 0,1667 0,1852
0,0131 0,0318 0,0506 0,0693 0,0880 0,1067 0,1255 0,1442 0,1629 0,1816
0,0094 0,0283 0,0472 0,0660 0,0849 0,1038 0,1226 0,1415 0,1604 0,1792
0,0113 0,0301 0,0489 0,0677 0,0865 0,1053 0,1241 0,1429 0,1617 0,1805
0,0771 0,0937 0,1104 0,127 0,1436 0,1602 0,1768 0,1935 0,2101 0,2267
0,0011 0,0300 0,0495 0,0687 0,0872 0,1093 0,1281 0,1469 0,1657 0,1842
0,0024 0,0377 0,0566 0,0755 0,0942 0,1132 0,1321 0,1509 0,1698 0,1887
1
т
При анализе табл. 3, 4 обращает на себя внимание тот факт, что, для расхода 172 м3/с эмпирическая обеспеченность равна 0,0042, а обеспеченность расхода 274 м3/с должна быть существенно меньше и по формуле (10) она равна 0,0024, т.е. практически в 2 раза меньше. Очевидно, что формулы (4)... (8) такой корректировки в оценке эмпирической обеспеченности выполнить не могут.
Расчеты максимума 1 % обеспеченности р. Шаугазсай, с учетом селевого максимума 1964 г. показали, что наименьшее в среднем отклонение расходов 1 % обеспеченности по сравнению с кривой Пирсона получено по формуле (9) - 28,7 %. Однако, по формуле (10) аналогичное отклонение составляет 32,3 % (табл. 5).
140
Таблица 5
Максимальные расходы воды 1 % обеспеченности, вычисленные по различным формулам (с включением селевого максимума 1964 г.)
и Максимальный расход воды 1 % обеспеченности, м3/с
о
н о о о & m m m - 0,3 m - 0,5 m - 0,4 1 f m , 051 m mn -1
ч о к s С n n +1 n+0,4 n n + 0,2 ■Jn +1 V*Jn J n + 1 + KZm n\ n + 12 1
о с V J
53 186 354 351 351 238 261 938 127 159
26 255 496 488 399 337 369 1230 178 188
18 300 568 557 459 389 425 1230 212 199
10 381 687 668 556 475 516 1230 271 174
Кратко сформулируем полученные выводы.
1. В случаях, когда нет возможности использовать то или иное теоретическое распределение, экстраполяцию кривых распределения в зону редкой повторяемости можно проводить по значениям эмпирической обеспеченности. Для этого достаточно использовать первые три-четыре члена убывающего ряда.
2. В соответствии с требованиями, предъявляемыми к статистическим оценкам рядов, оценка эмпирической вероятности, вычисляемая по формуле (10) имеет следующие преимущества:
- она является наименее смещенной по сравнению с остальными смещенными оценками;
- она обладает малой изменчивостью, что говорит об ее эффективности;
- эмпирическая обеспеченность, вычисленная по этой формуле, в наименьшей степени зависит от числа членов ряда.
3. В практических расчетах для малых рек целесообразно учитывать так называемые селевые максимумы, причем использование формулы (10) обеспечивает наименьшее увеличение расчетных величин (примерно в 1,3...1,4 раза), что, в какой-то степени, предотвращает не оправдываемое завышение максимумов и мало обоснованное удорожание проектируемых сооружений.
СПИСОК ЛИТЕРАТУРЫ
1. Виноградов Ю.Б. Математическое моделирование процессов формирования стока. - Л.: Гидрометеоиздат, 1988. - 312 с.
2. Горошков И.Ф. Гидрологические расчеты. - Л.: Гидрометеоиздат, 1979. - 430 с.
3. Длин А.М. Математическая статистика в технике. - М.: Наука, 1958. - 465 с.
141
4. Крицкий С.Н., Менкель М.Ф. Гидрологические основы управления речным стоком. - М.: Наука, 1981. - 255 с.
5. Лакин Г.Ф. Биометрия. - М.: Высшая школа, 1990. - 351с.
6. Мамедов М.А. Анализ условий формирования и методы расчета максимальных расходов горных рек (на примере рек Закавказья и Дагестана): Автореф. дис. ... докт. геогр. наук / Главный гидрологический институт.- Л., 1984. - 39 с.
7. Международное руководство по методам расчета основных гидрологических характеристик. - Л.: Гидрометеоиздат, 1984. - 247 с.
8. Пособие по определению расчетных гидрологических характеристик.
- Л.: Гидрометеоиздат, 1984. - 447 с.
9. Смирнов Н.В., Дунин-Барковский И.В. Курс теории вероятностей и математической статистики (для технических приложений). - М.: Наука, 1965. - 510 с.
10. Хальд А. Математическая статистика с техническими приложениями.
- М.: ИЛ, 1956. - 664 с.
11. Чуб В.Е. Изменение климата и его влияние на гидрометеорологические процессы, агроклиматические и водные ресурсы Республики Узбекистан. - Ташкент: Изд-во Узгидромета, 2007. - 132 с.
12. Чуб В.Е., Трофимов Г.Н., Меркушкин А.С. Селевые потоки Узбекистана. - Ташкент: Изд-во Узгидромета, 2007. - 109 с.
Национальный университет Узбекистана им. Мирзо Улугбека, г. Ташкент
ТАСЦЫНДАРДЬЩ ЕЦ ЖОГАРГЫ СУ ЭТ1МШЩ ЭМПИРИКАЛЬЩ ЦАМТАМАСЫЗДЫГЫН ЕСЕПТЕУ
Геогр. гылымд. докторы Г.Н. Трофимов
Жыл штде ец жогаргы су етгмтц эмпирикалъщ цамтамасъгздъггът барлъщ бащълау щатарындагы айырмашътыцты ескере отырып, есептеу мумшндШ щарастырътды. Мундай ыцгаймен эмпирикалъщ щаматамасъздыщты багалау гидрологияда дэстурл1 щолданътатын формулалармен есептегенге щараганда, аз цозгалмалы келедг. Сонъшен щатар щамтамасъздъщ щисъггътыц есебте су етгмтц селдж сипаттагы ец жогаргы мэндерт енг1зу цажеттыт жайлы тшр айтълады.
142