Проблемы создания нормативных документов в области расчета характеристик максимального стока рек Украины и возможные пути их решения Текст научной статьи по специальности «Строительство и архитектура»

ГИДРОЛОГИЯ

УДК 556.51

Е.Д. Гопченко, В.А. Овчарук, Е.И. Тодорова

ПРОБЛЕМЫ СОЗДАНИЯ НОРМАТИВНЫХ ДОКУМЕНТОВ В ОБЛАСТИ РАСЧЕТА ХАРАКТЕРИСТИК МАКСИМАЛЬНОГО СТОКА РЕК УКРАИНЫ И ВОЗМОЖНЫЕ ПУТИ

ИХ РЕШЕНИЯ

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

Ключевые слова: максимальный сток, паводки, половодье, модель изохрон.

E.D. Gopchenko, V.A. Ovcharuk, E.L. Todorova

PROBLEMS OF THE REGULATORY DOCUMENTS IN CALCULATING CHARACTERISTICS MAXIMUM RUNOFF RIVERS OF UKRAINE AND POSSIBLE SOLUTIONS

Basic calculation charts, in-use in the modern normative documents of Ukraine for the calculation of maximal flow of the rivers, are considered. The variants of improvement of calculation charts are offered on the basis of model of river-bed izokhron. Scientifically-methodical approaches are also offered for the account of possible changes of climate at determination of calculation descriptions of maximal flow of the rivers on statement model which enables to enter «climatic amendments» directly to maximal snow supplies and precipitations in the period of spring tide and floods

K e y w o r d s: maximal flow, floods, high water, model of izokhron.

Введение

В последние десятилетия XX столетия в Советском Союзе были разработаны и получили юридическое оформление две редакции нормативных рекомендаций по расчету характеристик максимального стока дождевых паводков и весенних половодий - сначала СН 435-72 [13], а затем СНиП 2.01.14-83[12]. Различия в их научно-методических подходах связаны с объемом исходной информации и пакетов различного рода вспомогательных карт тех или иных расчетных параметров.

В зависимости от происхождений максимальных расходов воды расчетные методики предусматривают и различную структурную базу [12; 13]. В частности, для весеннего половодья равнинных рек использована известная формула редукционного типа

k Y

--jkbT (1)

где qm - расчетный модуль стока; k0 - коэффициент дружности весеннего половодья; Ym - слой стока за половодье; F - площадь водосборов; b - эмпирический параметр, который учитывает замедление редукции отношения qm/Ym в области малых водосборов; fi - коэффициент, служащий для учета несовпадения статистических параметров рядов максимальных расходов и слоев стока; П1 - показатель степени редукции модуля qm, 8, £1 и £2 - поправочные коэффициенты, которые учитывают озерность, залесенность и заболоченность водосборов.

© Гопченко Е.Д., Овчарук В.А., Тодорова Е.И., 2016

Из недостатков, которые отражает формула 1, следует отметить главный - включение в знаменатель параметра Ь, который нарушает физическую сущность исходного уравнения. Согласно [4], числитель (1) представляет4 собой модуль склонового притока ^, так как

Ч'а = Ym = ^ , (2)

П То

где qа - максимальный модуль склонового притока; - - коэффициент временной

п

неравномерности притока воды со склонов в русловую сеть; Т0 - продолжительность притока воды со склонов в русловую сеть; ко - коэффициент склоновой трансформации (дружности) паводков и половодий.

С учетом формулы (2) и допущении, что ц = 1.0, 8 = 1.0, ¿>1 = 1.0, ¿>2 = 1-0:

аа

ь = (3)

Физическим условиям (qа|q'а = 1.0 ПРИ ¥=0) удовлетворяет только случай, когда Ь=1.0 км2. Следовательно, расчетная формула (1) должна использоваться в редакции

кЛ, „хх

'^2

qа ^А. (4)

(¥+1)

Применительно к дождевым паводкам в редакциях СН 435-72 и СНиП 2.01.14-83 в отдельные категории выделены малые водосборы (¥ < 200км2) и большие водосборы (¥)200км2).

Для водосборов с площадями ¥ < 200км2 нормативные документы рекомендуют формулу, которую многие авторы относят к структурам предельной интенсивности. В общем виде ее можно представить следующим образом [13; 14]:

qа = 16.67^(т)н сП, (5)

где qа - расчетный модуль поводочного стока, м3/скм2; = —— - редукционная кривая средней

Не

интенсивности осадков за расчетное время т(1/мин); Н с - суточный максимум осадков, мм; ^ -сборный коэффициент стока; 16.67 - коэффициент размерности.

Расчетная продолжительность т определяется по сумме времени руслового (/р ) и склонового (¿^ ) добегания, т.е.

т = 1.2?;1+1ск. (6)

С нашей точки зрения, неправомерным в структуре (5) является использование кривой редукции ^ (г) в качестве трансформационной функции распластывания паводочных волн вообще и на малых реках, в частности.

При площадях водосборов ¥)200км2:

( 200 ^ П2

qа = q2001 — 1 , (7)

где q2oo - максимальный модуль стока, приведенный к условной площади водосборов ¥ = 200кж2 ; П2 - показатель степени редукции максимального модуля стока дождевого паводка.

С методической точки зрения, вызывает возражение картирование модуля q2oo в пределах обширных территорий, поскольку при этом не учитываются различия в слоях стока и продолжительности притока, что хорошо видно, например, из более общей редукционной структуры (4).

Постановка и изученность вопроса

В свое время Г.А.Алексеев [1] и АН.Бефани [2], независимо друг от друга, процесс формирования стока на малых водосборах описали дифференциальным уравнением

Гидрология

дО дт ,

—+— = q,,

дх dt

(8)

где Q - расход воды; со - площадь живого сечения; qt - модуль склонового притока воды в русловую сеть.

Интегрирование (8) и последующее обобщение результатов позволило обосновать расчетное выражение для qm , а именно:

= г

(9)

где Ут - слой паводочного или половодного стока; tp - время руслового добегаиия; (р -коэффициент полноты склонового притока, принимающего участие в формировании максимальной ординаты руслового гидрографа и, согласно [4], для одномодальных паводков:

<Р = ~

n +1 tp

1 -

1 1t

n +11 To

(10)

n +1

где - - коэффициент временной неравномерности склонового притока, Т0 - продолжительность

п

притока воды со склонов в русловую сеть.

Если принять, что паводкоформирующие осадки Нт совпадают с их суточным максимумом Нс, то вместо (9) можно записать

qm = 16.67 ^V

(11)

Сравнивая (5) и (11), легко прийти к выводу, что

у{*)=<Р/ = ¿4 ), (12)

/ р

где tp - время руслового добегания, мин; ф[{р) - трансформационная функция, обусловленная распластыванием паводочных волн под влиянием времени руслового добегания, причем: а) при tp|T0 = 0;

~ \ п +11

ч>Чр )=

n T

б) при 0 < tp/T0 < 1.0;

в) при tp/T0 > 1.0;

г) при t >> To ;

~(,p )=

n + 1 1

1 -

1 11

n + 1 [ T

rffp )=r

(13)

(14)

(15)

<p(tp ) = 0. (16)

Изложенное позволяет формулу предельной интенсивности (5) представить в несколько иной редакции, а именно:

qm =<p(tp )н „ч. (17)

Уравнение (17) является универсальным, т.е. может применяться при расчете не только максимального стока паводков, но и весенних половодий, причем во всем диапазоне водосборных площадей.

В случае весеннего половодья

Гидрология

да = )SаV,

(18)

где 8а - максимальные снегозапасы перед началом половодья.

Чтобы расширить область применения дифференциального уравнения (8) относительно больших водосборов, А.Н. Бефани [2] на основе модели русловых изохрон представил его в редакции

дQ дю

-н--= ,

дх д,

(19)

где а - коэффициент густоты гидрографической сети; - модуль склонового притока; В, ширина водосбора по изохронам руслового добегания.

В результате интегрирования (19) получено обобщенное уравнение вида [2]

Y

qа = — V к, 'р

где кг - гидрографический коэффициент, учитывающий форму водосбора, причем: а) при 'р/Т0 = 0;

к = 1.0;

б) при 0</р/Т0 < 1.0;

в) при 1р/Т0 > 1.0;

г) при / >> Т0 ;

кг ,Л Ъср

0 /I 0 к г = В,л |в ср Т\4,л\;

к =

а +1

а +1

(20)

(21) (22)

(23)

(24)

коэффициент,

(25)

где В Ср - средняя ширина водосбора по изохронам руслового добегания, а

зависящий от формы водосбора и равный отношению

а +1 = Ва

а Вср '

где В а - максимальная ширина водосбора по изохронам руслового добегания. Опыт использования формулы (20) для расчета максимального стока паводков и половодий показал, что ее структурой не учитываются в полной мере эффекты руслопойменного и берегового регулирования. Поэтому А.Н. Бефани и Н.Ф. Бефани [3] предложили дополнить (19) параметрами, учитывающими руслопойменное регулирование. С этой целью базовое дифференциальное уравнение было представлено в редакции:

дю 1 дф ^

V-+--= ад, В,,

дх е д,

(26)

где е = -

ю „ю 11 + +

I ю ю

коэффициент руслопойменного и берегового регулирования

(осредненный во времени и по длине реки); ^ - площадь живого сечения потока по пойме; Фа -площадь живого сечения в аллювии; 8 - свободная пористость аллювия.

Интегрирование (26), выполненное по схеме (20), позволило несколько уточнить структуру формулы максимального стока, так как

а

Гидрология

qm =—

(27)

Практическое применение (27) показало, что е не является постоянной величиной, а зависит от размера водосборов. Это же замечание в полной мере относится и к коэффициенту густоты гидрографической сети а .

Теоретическое обоснование предлагаемой методики

На кафедре гидрологии суши ОГЭКУ предложен несколько иной метод построения нормативной базы в области максимального стока, основанный также на модели русловых изохрон. Он заключается в том, что исходя из редукционных гидрографов притока воды со склонов в русловую сеть ирусловых изохронфункции q^ и ft представлены уравнениями:

qt = q.

1 -

(28)

(29)

где То - продолжительность склонового притока, / - межизохронные площади.

Из-за отсутствия сведений о процессах руслопойменного и берегового регулирования будем полагать только, что оно является функцией времени, т.е. е = f ). Опираясь на теорию русловых изохрон:

а) при tp/Tо < 1.0;

б) при tp/Tо > 1.0;

^ =| q' ./а ,

(30)

Qm = | qtffit. (31)

0

Межизохронную площадку можно записать в виде

/ = V* Bt, (32)

где V* - скорость руслового добегания паводочных или половодных волн; Б1 - ширина водосбора

по изохронам руслового добегания. С учетом (28), (29) и (32): а) при tp/Tо <1.0;

^т

-- q'm£tp-— Vg tp т +1

т +1

(п + 1)(т + п +1)1 То

б) при tp/Tо > 1.0;

т +1 п +1

т + 1 п + 1 т т(т + п н

(33)

(34)

п{т + п +1)^ t

где ^т - приближенное решение интегралов (30) и (31), связанное с осреднением функции £с по tp и Т0, соответственно.

п

t

т

t

t

1

т

Т

о

Гидрология

Для корректировки результатов к (33) и (34) вводятся переходные коэффициенты к£ = ^о) . С их учетом

да

■д а

(35)

где в = к в и в _ = к ВТ - коэффициенты русло-пойменного и берегового регулирования;

£ 1р £ о

(' / ]

I - трансформационная функция распластывания паводочных волн под влиянием времени

руслового добегания: а) при 1р/Т0 = 0;

) при 0</р/Т, < 1.0;

в) при 1р/Т0 > 1.0;

Г) при ' „ ;

= 1.0

т +1

V

(п + 1)(от + п +1)

п + 1 ,

а +1 п +1 а а{т + п +1)

(т ^а

г, ^

= 0.

(36)

(37)

(38)

(39)

Входящий в (35) модуль склонового притока д а рассчитывается с использованием уравнения (2). а

Обсуждение результатов

Теоретические и методические преимущества расчетной формулы (35), по сравнению с другими подобными схемами, состоят в следующем:

1. Предложенный вариант опирается на природный процесс трансформации осадков в русловой сток в составе звена «склоновый приток - русловой сток».

2. Она применима во всем диапазоне водосборных площадей - от отдельных склонов до крупных речных систем.

3. Формула является универсальной и с точки зрения ее использования для различных генетических типов паводков и половодий. Для реализации предложенной расчетной формулы максимального стока паводков и половодий на кафедре гидрологии суши ОГЭКУ разработан программный комплекс численного установления в структуре (35) неизмеряемых величин Т , п +1

К операторному виду «склоновый сток - русловой сток» можно привести и структуры предельной интенсивности для малых водосборов (17) и (18). Для этого достаточно числитель и знаменатель в их

п +1 1 „

правой части домножить на ^ =--. В результате для малых водосборов получим обобщенную

п Т

формулу

Т ^

о У

т

0 У

п

^ т^/т0 )т0 ^, (40)

где qт - максимальный модуль склонового притока, равный: а) для дождевых паводков

) для весенних половодии

п +1 1 тт

:--Нл ; (41)

п Т

п +1 1 „

: -— ЗтЧ- (42)

п Т0

Произведение (p{tp|Tо )Т0 представляет собой безразмерную, т.е. нормированную по Т) трансформационную функцию (p{tp ).

Предложенная научно-методическая база для нормирования расчетных характеристик паводков и половодий реализована на примере рек Днепра, Днестра, Дона, Крыма, Карпат, Алжира и др.[5-7; 9; 11; 15].

Учитывая происходящие в последнее время региональные и глобальные изменения климата, при расчетах стоковых характеристик возникает необходимость введения соответствующих поправок. В частности, в [10] предлагается, опираясь на действующий в России нормативный документ СП 33101-2003 [14], рассчитывать соответствующие характеристики стока, а затем вводить к ним так называемую «климатическую поправку». В качестве проектного предлагается брать большее из полученных значений, т.е. с поправкой и без нее.

В целом методика, предложенная учеными РГГМУ, по нашему мнению, может быть использована на практике, но учет «климатической поправки» через начальные моменты является несколько опосредствованной и сложной процедурой.

Операторная модель дает возможность вводить «климатические поправки» непосредственно к максимальным снегозапасам и стокоформирующим осадкам в период весеннего половодья и паводков, а также к коэффициентам стока. Подробно эта проблема рассмотрена в статье авторов [8], а в общем можно отметить, что все прогнозные модели и сценарии ведущих мировых центров, которые представлены на сайте Центра распространения данных (ВВС) Межправительственной группы по климатическим изменениям (1РСС, http://ipcc-data.org), дают возможность получить прогнозные средние значения температуры воздуха и осадков. Авторами проведено исследование взаимосвязей между прогнозируемыми и расчетными величинами - максимальными снегозапасами и осадками во время половодья. Для территории равнинной Украины по данным 103 метеостанций по снегозапасам и 315 метеостанций по осадкам получены соответствующие зависимости, которые имеют достаточно высокие коэффициенты корреляции.Что касается максимального стока дождевых паводков, то здесь ситуация несколько иная. Максимальные суточные осадки 1 % вероятности превышения имеют незначительную тенденцию к увеличению с ростом среднего количества осадков и практически не зависят от средней многолетней температуры воздуха.

Ч

Чт

чт

Выводы

Действовавшие в пределах всей территории бывшего СССР до 2003 г. рекомендации по расчету характеристик максимального стока весенних половодий и дождевых паводков опирались на нормативные документы СН 435-72, а затем СНиП 2.01.14-83. В их основу были положены главным образом формулы редукционного типа (кроме малых водосборов при расчете максимальных расходов воды дождевых паводков, для которых использован упрощенный и весьма схематичный вариант формулы предельной интенсивности). Нормативные документы сопровождались прикладными пакетами карт и вспомогательных таблиц.

Авторы статьи видят свою задачу в разработке предложений по созданию в Украине более теоретически обоснованных расчетных схем. При этом очевидно следующее:

1. Нецелесообразно принимать для нормирования характеристик паводков и половодий различную структуру расчетных формул.

2. Вызывает возражение разделение расчетных схем для паводков на водосборах с площадями Р>200 км2 и Р<200 км2. С другой стороны, формула предельной интенсивности и параметры, которые

Гидрология

в нее входят (в редакции рассмотренных нормативных документов), не имеют соответствующего теоретического и методического обоснования.

3. Авторами обосновывается структура расчетной формулы максимального стока, в которой учитываются характеристики склонового притока (слой стока за паводок или половодье, продолжительность притока и временная неравномерность притока воды со склонов в русловую сеть).

4. Предлагаемая структура дает возможность учитывать «климатические поправки» непосредственно, например, через максимальные снегозапасы и осадки в период половодья.

Статья представляется для обсуждения в кругу специалистов с целью дальнейшего развития научно-методической базы действующих и вновь создаваемых нормативных документов в области максимального стока. В частности, такой опыт уже есть - при создании нормативного документа по расчету характеристик стока рек Молдовы для расчета максимального стока рекомендуется операторная модель [16].

Библиографический список

1. Алексеев Г.А. Анализ и расчет максимального стока // Тр.НИУГМС. 1946. Сер. IV. Вып. 29. С. 33-91.

2. БефаниА.Н. Основы теории ливневого стока // Тр. ОГМИ. 1955. Вып. XIV. Ч. II. 302 с.

3. Бефани А.Н., Бефани Н.Ф.,Гопченко Е.Д. Региональные модели формирования паводочного стока на территории СССР. Обнинск, 1981. Вып. 2. 60 с.

4. Гопченко Е.Д. О редукции максимальных модулей дождевого стока по площади // Метеорология и гидрология. 1975. №2. С. 66-71.

5. Гопченко С.Д., Ладжель Махмуд, Овчарук ВА. Максимальный сток паводков на уэдах северной части Алжира // BicHHK Одеського державного еклолопчного ушверситету. 2007. Вип. 4. С. 287-296.

6. Гопченко Е.Д., Овчарук ВА. Формирование максимального стока весеннего половодья в условиях юга Украины. Одесса: ТЕС, 2002. 110 с.

7. Гопченко Е.Д, Романчук М.Е. Нормирование характеристик максимального стока весеннего половодья на реках Причерноморской низменности. Киев, 2005. 148 с.

8. Гопченко С.Д., Овчарук В.А., Семенова 1.Г. Науково-методичш шдходи до врахування глобальних змш кл1мату при розрахунках максимального стоку р1чок // Вюник Одеського державного еклолопчного ушверситету. 2012. Вип. 14. С. 141-150.

9. Гопченко С.Д., Овчарук В.А., Шаюрзанова Ж.Р. Розрахунки та довгостроков1 прогнози характеристик максимального стоку весняного водошлля в басейш р. Прип'ять. Одесса: Еколог1я, 2011. 335 с.

10. Коваленко В.В. Гидрологическое обеспечение надежности строительных проектов при изменении климата. СПб.: РГГМУ, 2009. 100 с.

11. Овчарук В.А. Особенности применения формулы операторного типа для расчета максимального стока рек в различных физико-географических условиях // Метеоролог1я, кл1матолопя та пдролопя. Одесса: Еколопя, 2008. Вип. 50. Ч. 2. С. 93-101.

12. Пособие по определению расчетных гидрологических характеристик. Л.: Гидрометеоиздат, 1984. 447 с.

13. Руководство по определению расчетных гидрологических характеристик. Л.: Гидрометеоиздат, 1973. 111 с.

14. СП 33-101-2003. Определение основных расчетных гидрологических характеристик. М.: Госстрой России, 2003. 74 с.

15. Eugene Gophenko, Valeriya Ovcharuk. Theoretical ground of normative base for calculation of the characteristics of the maximum runoff and its practical realization // Transboundary Floods: Reducing Risks Through Managment. NATO Sciences Series. IV Earth and Environmental Sciences. Vol. 72. Springer, 2006. P. 91-99.

16. Determinarea caracteristicilor hydrologice pentru Conditiile republicii Moldova. CPD.01.05-2012. 178 p.

Поступила в редакцию 24.12.2016

Сведения об авторах About the authors

Гопченко Евгений Дмитриевич

доктор географических наук, профессор, заведующий кафедрой гидрологии суши Одесского государственного экологического университета;

Украина, 65016, Одесса, ул. Львовская 15; e-mail: dean-gid@odeku.edu.ua

Овчарук Валерия Анатольевна

кандидат географических наук, доцент кафедры гидрологии суши Одесского государственного экологического университета, Украина, 65016, Одесса, ул. Львовская 15; e-mail: dean-gid@odeku.edu.ua

Тодорова Елена Ивановна

аспирант кафедры гидрологии суши Одесского государственного экологического университета, Украина, 65016, Одесса, ул. Львовская 15; e-mail: lenochka.todorova@mail.ru

Eugeniy D. Gopchenko

Doctor of Geographical Sciences, professor, Head of the Department of Land Hydrology, Odessa State Environmental University; 15, Lvovskaya Str., Odessa, 65016, Ukraine; e-mail: dean-gid@odeku.ua

Valeriya A. Ovcharuk

Candidate of Geographical Sciences, associate professor of hydrology, director of hydrometeoro logical Institute, Odessa State Environmental University; 15, Lvovskaya Str., Odessa, 65016, Ukraine; e-mail: dean-gid@odeku.ua

Olena Todorova

Postgraduate student of the Department of Land Hydrology, Odessa State Environmental University; 15, Lvovskaya Str., Odessa, 65016, Ukraine; e-mail: dean-gid@odeku.ua

Просьба ссылаться на эту статью в русскоязычных источниках следующим образом:

Гопченко Е.Д., Овчарук В.А., Тодорова Е.И. Проблемы создания нормативных документов в области расчета характеристик максимального стока рек Украины и возможные пути их решения // Географический вестник. 2016. №1(36). С. 49-57. Please cite this article in English as:

Gopchenko E.D., Ovcharuk V.A., Todorova E.L. Problems of the regulatory documents in calculating characteristics maximum runoff rivers of Ukraine and possible solutions // Geographicheskiy Vestnik. 2016.

№1(36). P. 49-57.