CC BY
0к.т.н. Павлов В. И., Кусайко Н. П., Кулакова С. И.
(ДонГТИ, г. Алчевск, ЛНР, pavlow2005@rambler.ru)
АНАЛИЗ ИЗМЕНЕНИЯ СОСТАВЛЯЮЩИХ ВОДНОГО БАЛАНСА БАССЕЙНА ВОДОСБОРА ИСАКОВСКОГО ВОДОХРАНИЛИЩА
Проведено сравнительное исследование динамики изменения среднемноголетних годовых притоков в Исаковское водохранилище, нормы атмосферных осадков и температурной нормы. Установлено, что вследствие растущего повышения температуры воздуха среднее значение естественного природного стока с периодом усреднения в 20 лет потеряло свойство стационарности; в последнее десятилетие наметилась устойчивая тенденция уменьшения доли стока в водном балансе бассейна водосбора за счет увеличения доли испаряемой влаги.
Ключевые слова: водохранилище, бассейн водосбора, притоки, подземный сток, поверхностный сток, климатическая норма, среднемноголетнее значение, стационарность среднего, водный баланс, доля стока, доля испарения, изменение климата.
УДК 556.047
Проблема и ее связь с научными и практическими задачами. В последнее время в Луганской Народной Республике возникла острая проблема водоснабжения промышленных и гражданских объектов, вызванная недопустимо большим обмелением используемых водохранилищ. Во второй половине прошедшей летней межени уровень воды в Исаковском водохранилище был ниже уровня гарантированной водоотдачи, и поэтому водозабор на нужды Алчевского металлургического комбината был ограничен. На фоне переживаемого периода малой водности рек ситуация усугубилась переходом к мокрой консервации ряда угольных шахт, в результате чего прекратилась подпитка притока водохранилища шахтной водой. Возникла реальная проблема водопользования Исаковским водохранилищем. Для принятия инженерных решений, обеспечивающих эксплуатационную наполняемость водохранилища, необходимо в первую очередь изучить гидрологические процессы, протекающие в бассейне его водосбора. Важность этих исследований заключается в следующем.
Если критическое снижение притока в водохранилище вызвано прекращением
его подпитки водоотливами шахт «Романовская» и им. Артема, а наблюдаемое снижение естественного стока является возможным отклонением в диапазоне принятой обеспеченности, то очевидным решением будет возобновление подпитки притока шахтной водой.
Если же наблюдаемая низкая величина естественного стока выходит за диапазон значений малой водности и является следствием глобального изменения климатических условий, то возобновление откачки шахтной воды будет недостаточным для решения возникшей проблемы.
Величина естественного стока может уменьшаться за счет возникновения на территории бассейна водосбора зон повышенной инфильтрации, вызванной подтоплением выработанных пространств при прекращении работы шахтных водоотливов.
Анализ исследований и публикаций. Величина естественного стока является основой для проектирования гидротехнических сооружений и рассчитывается исходя из модуля речного стока, стационарность которого является основополагающей научной концепцией гидрологии. На основе этой концепции в прошлом веке были созданы гидрологические карты для всей тер-
ритории бывшего Советского Союза. В работе К. П. Воскресенского [1] отмечено, что гидрологические сооружения на малых реках (пруды, водохранилища), как правило, проектируются исходя из данных стока ближайших средних рек. Такой подход приводит к завышенным оценкам. Для их уточнения необходимы непосредственные гидрологические исследования бассейна водосбора проектируемого сооружения. Первый вопрос, возникающий при этом, — минимально достаточная длительность наблюдений, так как накопление статистических данных весьма трудоёмко. Исходя из большого количества влияющих факторов, величина стока является величиной случайной, и поэтому увеличение длительности периода наблюдений более 10-15 лет не приводит к значимому уточнению среднегодового значения [2].
Водный баланс определяет соотношение между поступлением воды на территорию бассейна водосбора и её расходованием [3]. Из уравнения водного баланса речного водосбора следует, что речной сток является результатом активного во-дообменного цикла и определяется прежде всего атмосферными осадками [1]:
У = X - z ± п, (1)
где У — сток с бассейна водосбора;
X — атмосферные осадки;
г — испарение;
п — запасы влаги на водосборе.
Научными исследованиями установлено, что вследствие потепления скорость обменного водного цикла (атмосфера — осадки — стоки — река — море и снова атмосфера) в целом по планете постоянно из года в год увеличивается [4]. Повышение температуры приводит к увеличению испарения с земной и водной поверхности мирового океана, к возрастанию количества влаги в атмосфере. В водном балансе испарение, как функция температуры воздуха, является отрицательным фактором, поэтому при увеличении осадков сток может уменьшаться [5].
Из (1) следует, что годовой сток может иметь, кроме ежегодного объёма осадков, дополнительные источники пополнения, что существенно усложняет намеченный анализ. Такими источниками могут быть шахтная вода, колодцы, скважины, которые используют воду из замедленного во-дообменного цикла, а также сбросы очистных сооружений населенных пунктов и предприятий.
Научной целью работы является исследование изменения основных составляющих водного баланса бассейна водосбора: поступающего объёма влаги с атмосферными осадками и притока в водохранилище, как следствие естественного стока. В конкретной климатической зоне известна климатическая норма атмосферных осадков: величина постоянная и не зависящая от времени, а наблюдаемые из года в год колебания осадков связываются с естественными случайными факторами. Наличие нормы осадков обуславливает аналогичную привязку к климатической зоне гидрологической нормы — коэффициента стока или модуля речного стока, которая также не должна зависеть от времени. Вместе с тем, согласно [6] основные гидрологические характеристики могут изменяться под влиянием климатических факторов.
Таким образом, в работе решаются такие задачи:
- анализ колебания атмосферных осадков в бассейне водосбора Исаковского водохранилища за весь период его эксплуатации;
- анализ колебаний притоков в водохранилище за этот же период;
- исследование взаимосвязи основных составляющих водного баланса водосбора: притока в водохранилище, как результат (функция Y), и влияния атмосферных осадков, как аргумент Х функции Y.
Изложение материала и его результаты. Исаковское водохранилище было сдано в эксплуатацию в 1953 г., почти семь десятилетий назад. Намного раньше, уже с 1913 г. в бассейне р. Белая работали шахты и попол-
няли её шахтной водой, поэтому приток, заложенный в обоснование объёма водохранилища, был больше естественного речного стока. Следствием ведения горных работ под территорией водосборного бассейна являются зоны с повышенной техногенной инфильтрацией. Объёмы вынимаемых угольных пластов не полностью заполняются обрушенными породами. Образовавшаяся пус-тотность вследствие постепенного оседания подработанных породных слоев распределяется по всему массиву вплоть до земной поверхности. Во время работы шахты эта пус-тотность служит инфильтрационной сетью, создающей приток воды в горные выработки. После прекращения работы шахтных водоотливов по мере заполнения горных выработок и породной пустотности уровень воды повышается до уровня водоносных горизонтов активного водообмена, за счет чего увеличивается и речной сток [7]. На последних стадиях заполнения пустотности из-за влаго-насыщения пород снижается их прочность и возникают масштабные оседания земной поверхности с выдавливанием накопленных объёмов воды в речные стоки. В нижних течениях рек могут происходить даже подтопления отдельных территорий. Водоотливы шахт могут пополнять приток. Шахтные поля, находящиеся за пределами бассейна водосбора, могут увеличивать приток в горные выработки шахт, подпитывающих водосборные водотоки. Создаваемые шахтами ин-фильтрационные зоны являются причиной потерь речного стока. Это влияние существовало весь предшествующий период эксплуатации водохранилища, постепенно убывая по мере закрытия шахт (табл. 1).
Естественный речной сток формируется активным водообменным циклом (рис. 1), включает поверхностный сток (Сп ) и приповерхностный сток (Гл ) подземных водоносных горизонтов. Инфильтрационные процессы перетекания осадков в более глубокие горизонты создают в водном балансе (1) запасы влаги в водосборе ы, которые участвуют в замедленном водообменном
цикле. Для бассейна водосбора их следует принимать со знаком «минус» и трактовать как потери. Особенно значительные потери ы возникают в недрах, подработанных шахтами, где накапливаются застойные, мертвые объёмы воды с высоким содержанием растворенных веществ.
Естественный водообмен можно представить уравнением водного баланса
= Сп + Л_2 + И + Ггл, (2)
где Сп, Г1_2 — годовой поверхностный и подземный объёмы стока в водохранилище;
И — годовой объём испаряющейся воды с бассейна водосбора;
Ггл — годовой объём инфильтрации в глубокие водоносные горизонты;
Жг — годовой объём атмосферных осадков, приходящийся на площадь водосборного бассейна:
_3 з
Жг = h ■ F -10 , млн м ,
где h — годовые осадки, мм;
F — площадь бассейна водосбора Исаковского водохранилища, F = 450 км .
Таблица 1 Шахты, изменявшие приток в Исаковское водохранилище
№ Название шахты Период эксплуатации
начало конец
1 Украина (Романовская) 1932 2012(2020)
2 5-бис (Перевальская) 1913 2017
3 2-бис
4 Селезневская 1971
5 9-бис 1913 1970
6 им. Артема 1913 2020
7 Никанор
8 Комиссаровская 1955 2003
9 им. Косиора 2001
10 Фащевская 1940
11 153 (возле п. Запорожье) 1913 1997
12 106
13 Штеровская 1960 1999
14 Белореченская 1957
Региональная экология и природопользование
направление водообмена:
--активного
----— замедленного
Рисунок 1 Схема активного и замедленного водообмена в бассейне водосбора
Так как норму атмосферных осадков следует считать величиной относительно постоянной и её исходная статистическая база довольно точно и систематически пополняется, то для региона водосбора важной гидрологической характеристикой будет доля речного стока в атмосферных осадках, принимаемых за 100 %. Выражение (2) примет вид
Ср И Г 100 = 100 + —100 + 100.
Ж,
Ж,
Ж,
(3)
Ср
где —-100 = 5р,% — доля речного стока.
Ж р
Приток с бассейна водосбора содержит речной сток Ср и сбросы в основном шахтных вод Вш (рис. 1) и устанавливается при расчете годового водохозяйственного баланса водохранилища.
На начальном этапе анализа были составлены три временных ряда: в каждый текущий год эксплуатации водохранилища рассчитан среднемноголетний годовой объём притока Псрг путем усреднения за прошедшие 20 лет годовых объёмов притока. Аналогично были составлены ряды нормы атмосферных осадков Жср.г и доли притока в атмосферных осадках ёп.сря. Текущее значение ряда Зп.ср.г определялось по формуле:
Региональная экология и природопользование
П,
ср.г
псрг W
'' ср.г
100,%.
(4)
Период усреднения также был принят равным 20-ти годам. Данные по осадкам предоставлены ближайшим метеорологическим постом на территории Луганского государственного аграрного университета.
Изменение показателей Wсp,г и Псрг во времени, полученное методом скользящего среднего [8] (рис. 2), отражает сопряженную динамику развития гидрологических и атмосферных процессов. Тренды этих процессов за весь период эксплуатации водохранилища, отмеченные сплошной и пунктирной линиями, свидетельствуют о противоположных тенденциях — объёмы атмосферных осадков увеличиваются, а притоки в водохранилище уменьшаются.
Такая взаимосвязь показателей Wсp.г и
Пср.г противоречит естественным атмо-сферно-гидрологическим процессам — чем больше объёмы осадков, тем больше должны быть подземные и поверхностные
стоки в поверхностные водные объекты. Если рассмотреть взаимное изменение показателей более детально, то можно выделить периоды, различающиеся по виду их общей динамики.
Весь временной интервал, представленный на рисунках 2 и 3, был разделен на шесть периодов. Периоды № 2 и № 4 характеризуются резким снижением притоков при увеличении осадков.
Причины такого падения притоков очевидны и были обусловлены негативными процессами, происходившими в угольной промышленности Донбасса. Объёмы выкачиваемой шахтной воды резко уменьшались вследствие сокращения количества работающих шахт, соответственно уменьшалась подпитками ими Исаковского водохранилища. Эти негативные процессы в бассейне рассматриваемого водосбора и в целом по отрасли в период № 4 имели общие тенденции, что однозначно наблюдается на рисунке 4. Данные по количеству действующих газовых шахт взяты из работы [9].
Рисунок 2 Тренды среднемноголетних годовых объёмов притоков Псрг в Исаковское водохранилище и осадков Жср.г на водосборную площадь
Региональная экология и природопользование
^ 22,0 -
I 21,0
4 19,0
а
5 18,0 H
17,0
16,0
N 1
V4 н JMsti
№1 Ys2 №3 №4 JfeS
с«*
w
J >
1 s fi У h L J JDnDOL f- Л< 4
1
10,0
' и
g Cl
о
w§ в-
9,0
8*8 i-t §
8.6
8,4 S
S
8.2
1970 1975 1980 1985 1990 1995 2000 2005 2010 2015 2020
~~ - доля притока: -□- - температурная норма
- гидрологический период:----- линия тренда
Рисунок 3 Изменение доли притока в атмосферных осадках и среднегодой температурной нормы
по гидрологическим периодам
В периодах № 1, 3, 5 и 6 (рис. 2 и 3) визуально наблюдается естественная взаимосвязь между осадками Жср.г, притоками Псрг и величиной Зср.г. Статистические характеристики этих показателей по выделенным периодам приведены в таблице 2. Средние, минимальные и максимальные значения осадков увеличивались с периода № 1 по № 5 (табл. 2). В последний период по всем этим характеристикам произошло снижение. Среднеквадратическое отклонение от среднего (СКО) было максимальным в период № 1, а в остальные периоды происходил рост неравномерности осадков.
Объёмы притоков по периодам постоянно снижались, что в первую очередь обуславливалось закрытием угольных шахт. Следует отметить, что в период № 1 статистические характеристики притоков были наибольшими, а характеристики осадков были наименьшими. При этом происходило развитие угольной промышленности. Период № 1 фактически охватывает среднемноголетние значения притоков с двадцатилетним диапазоном ус-
реднения от 1953 г. до 1977 г., когда большинство шахт вели работы на глубине водоносных горизонтов в зоне активного водообмена и шахтные водоотливы увеличивали скорость водообмена. В периодах № 3, 5 и 6 неравномерность (СКО) притоков и осадков росла согласованно (табл. 2). Доля притока в осадках уменьшалась по мере закрытия угольных шахт, но в последние периоды стабилизировалась.
Корреляционно-регрессионный анализ выполнялся с помощью программного приложения Statistica. Исходя из графического представления данных на рисунках 2 и 3 для зависимостей Пср г от Жср.г и ёср,г от текущего времени (Т, лет) проверена линейная форма связи:
У = а ■ Х + Ь. (5)
Следует сразу отметить, что наличие положительной зависимости Псрг от Жср.г отражает изменение притока за счет речного стока, так как доля сточных вод в притоке зависит больше от антропогенных факторов.
Региональная экология и природопользование
м
290
270
250
230
210
190
170
150
о-е-о-е-е^
А Ч
Я V \ * ч \
/ \ \ J»
^ V*
№4
1990
1995
2000
2005
2010
Годы
действующие i притоки
49
47
45 Щ
43
41
39
37
35
2015
Рисунок 4 Взаимосвязь общей динамики сокращения действующих газовых шахт с притоками Пср.г в Исаковское водохранилище (иллюстрация причины резкого снижения притоков)
В таблицах 3 и 4 приведены основные результаты анализа:
- значимость по критерию Стьюдента (при уровне значимости 0,05) коэффициентов корреляции R и детерминации R , коэффициентов регрессии (наличие значимости выделена жирным шрифтом); значимость установленной модели функциональной зависимости по критерию Фишера при уровне значимости 0,05;
- a и Ь — величины коэффициентов регрессии;
- £ — среднеквадратическая ошибка вычисления функции;
- Э — коэффициент эластичности, %.
Установлены следующие зависимости
(табл. 3 и 4):
для периода № 1
5п.сР.г =-0,14 • Т + 24,6; (6)
для периода № 3
Псрг = 0,33 -ЖсРг - 28,66, (7)
5П£р,г = 0,12 • Т +19,34; (8)
для периода № 5
ПсРг = 0,40 • Жсрг - 47,80; (9)
для периода № 6
Пср.г = 0,45 • Wcp z - 63,05, (10)
5п.ср.г =-0,32 • Т +19,01.
(11)
Влияние осадков на приток в водохранилище по периодам № 3, 5, 6 адекватно описывается соответственно линейными зависимостями (7, 9, 10), отражающими положительную связь, увеличение осадков приводит к увеличению притоков. Очевидно, в эти периоды были относительно стабильными сбросы шахтных вод, в результате чего проявились естественные закономерности активного водообмена. Тогда был отмечен рост доли притока в осадках (период № 3 на рисунке 3), что подтверждается зависимостями (7) и (8).
Последний период отличается от всех остальных. Тенденция роста осадков резко меняется на падение. Из уравнений (7), (9) и (10) следует, что наибольшее положительное влияние осадков на среднегодовой приток наблюдается в период № 6. Согласно коэффициенту эластичности, снижение осадков на 1 % приводит к наибольшему уменьшению притока, на 2,58 %. Вместе с тем, как следует из уравнения (11), со временем уменьшается процент притока в водном балансе.
Региональная экология и природопользование
Статистические характеристики притоков и осадков по периодам
Таблица 2
Показатель Период Количество Среднее Мин. Макс. СКО
№ годы значений
1 1972-1977 6 203,67 195,00 214,00 6,25
Осадки, 3 1984-1997 14 215,64 211,00 222,00 3,93
млн. м3/год 5 2002-2011 10 234,06 226,53 239,76 4,34
6 2011-2019 9 226,56 220,00 238,00 5,34
1 1972-1977 6 49,17 48,00 51,00 1,17
Притоки, 3 1984-1997 14 43,57 41,72 46,82 1,65
млн. м3/год 5 2002-2011 10 40,99 38,21 43,72 1,84
6 2011-2019 9 39,50 35,90 43,70 2,78
Доля притока в осадках, % 1 1972-1977 8 23,96 23,30 24,70 0,44
3 1984-1997 13 20,15 19,50 21,10 0,52
5 2002-2011 10 17,51 16,90 18,40 0,51
6 2011-2019 9 17,41 16,13 18,40 0,90
Корреляционно-регрессионный анализ зависимости Жср.г от Пс
1ср.г
Таблица 3
Период R R2 а Ь 5 Э Форма связи
№ годы У = аХ+ Ь
1 1972-1977 0,80 0,64 0,15 -18,58 0,78 - не значима
3 1984-1997 0,80 0,64 0,33 -28,66 1,03 1,66 значима
5 2002-2011 0,90 0,80 0,40 -47,80 0,87 2,17 значима
6 2011-2019 0,87 0,76 0,45 -63,05 1,46 2,58 значима
Таблица 4
Корреляционно-регрессионный анализ зависимости ёлср.г от времени
Период R R2 а Ь 5 Форма связи
№ годы У = аХ + Ь
1 1972-1977 -0,78 0,62 -0,14 24,6 0,29 значима
3 1984-1997 0,88 0,77 0,12 19,34 0,26 значима
5 2002-2011 0,63 0,40 0,11 16,92 0,42 не значима
6 2011-2019 -0,97 0,94 -0,32 19,01 0,22 значима
Причина таких изменений, скорее всего, связана с наступившими климатическими изменениями в нашем регионе, характеристику которых очень убедительно добавляет аналогичный временной ряд изменения температурной нормы воздуха (рис. 3). С третьего периода по пятый включительно среднегодовая температурная норма воздуха совпадает с ростом осадков и положительным влиянием на притоки. Дальнейший рост температуры в последний, шестой период, привел, как подтверждают климатологи, к резким негативным по-
следствиям — в водном балансе бассейна водосбора происходит существенное снижение поверхностного и подземного речного стока за счет увеличения доли осадков, приходящейся на испарение.
Изменение климатических норм во времени (среднемноголетних годовых значений атмосферных осадков, годовых температурных норм), среднемноголетних годовых значений в периоды № 3, 5 и 6 (рис. 2 и 3), наличие значимых связей (7, 9, 10) свидетельствует об отсутствии стационарности средних значений в используемом
двадцатилетнем периоде усреднения. Такие последствия изменения климата подтверждаются исследованиями [6]. Для практического использования результатов гидрологических исследований стандартом [6] требуется уточнение периода стационарности (квазистационарности).
Выводы. В последнее десятилетие наличие зависимостей между среднемноголет-ними годовыми значениями притоков и годовой нормой атмосферных осадков, между долей притока в водном балансе и временем свидетельствует об изменении климатических условий в бассейне водосбора Исаковского водохранилища, нестационарности величины средних годовых притоков.
Проявляющиеся противоположные тенденции увеличения температурной нормы и снижения нормы годовых атмосферных осадков являются причиной снижения поверхностного и подземного стока.
Библиографический список
Дополнительным негативным проявлением климатических изменений является снижение доли притока в водном балансе бассейна водосбора за счет увеличения доли испарения.
Из полученных результатов следует, что для оценки модуля речного стока в бассейне реки Белая необходимо уточнение:
- длительности квазистационарности (неизменчивости) модуля речного стока, обеспечивающего более надёжное прогнозирование наполняемости водохранилища;
- объёмов водопотребления и сбросов в бассейне водосбора;
- зонирования повышенной глубинной инфильтрации атмосферных осадков в пределах шахтных полей на территории бассейна;
- влияния изменения климата на водный режим водотоков путем их посезонного мониторинга.
1. Воскресенский, К. П. Гидрологические расчеты при проектировании сооружений на малых реках и временных водотоках. Методические основы и практика [Текст] / К. П. Воскресенский. — Ленинград : Гидрометеорологическое издательство, 1956. — 467 с.
2. Воскресенский, К. П. Норма и изменчивость годового стока рек Советского Союза [Текст] : монография / К. П. Воскресенский. — Ленинград : Гидрометеорологическое издательство, 1962. — 540 с.
3. Многолетняя динамика показателя «баланс влаги» в Луганской области [Текст] // Е. Д. Долгих, Соколов И. Д. и др. // Зб1рн. наукових праць ЛНАУ. — Луганськ : Елтон-2, 2006. — № 66 (89). — С. 135-142.
4. Антропогенные изменения климата в ХХ1 веке в северной Евразии [Текст] /В. П. Мелешко, В. М. Катцов, В. А Говоркова и др. //Метеорология и гидрология. — 2004. — № 7. — С. 5-26.
5. Изменение климата Луганщины и его прогнозирование. Основания для оптимизма [Текст] : монография / И. Д. Соколов, М. В. Орешкин, О. М. Медведь и др. — Луганск : ФЛП, 2017. — 200 с.
6. СТО ГТИ 52.08.41-2017. Основные гидрологические характеристики при нестационарности временных рядов, обусловленной влиянием климатических факторов [Текст]. — СПб. : ФГБУ «ГГИ», 2010. — 42 с. — (Стандарт организации Росгидромета).
7. Ворхлик, И. Г. Технология закрытия (ликвидации) угольных шахт [Текст] : учеб. пособ. для вузов /И. Г. Ворхлик, В. И. Стрельников, И. Ф. Ярембаш. — Донецк : ДонНТУ, 2004. — 238 с.
8. Кендэл, М. Временные ряды [Текст] / пер. с англ. и предисл. Ю. П. Лукашина. — М. : Финансы и статистика, 1981. — 199 с.
9. Кочерга, В. М. Обгрунтування параметр1в способу дегазацИ' вигмковог дыьниц з використанням свердловин у неконтрольованих виробках [Текст] : дис. .. канд. техн. наук / В. М. Кочерга. — Днтро, 2018. — 254 с.
© Павлов В. И.
© Кусайко Н. П.
© Кулакова С. И.
Региональная экология и природопользование
Рекомендована к печати к.э.н., зав. каф. высшей математики ДонГТИМельничук Д. А., помощником министра природных ресурсов и экологической безопасности ЛНР
Крамаренко А. А.
Статья поступила в редакцию 22.05.2021.
PhD in Technical Sciences Pavlov V. I., Kusayko N. P., Kulakova S. I. (DonSTI, Alchevsk, LPR, pavlow2005@rambler. ru)
ANALYSIS OF CHANGES IN WATER BALANCE COMPONENTS OF THE ISAKOVO RESERVOIR CATCHMENT
A comparative study of the change dynamics in long-term annual inflows into the Isakovo storage lake, the norm of atmospheric precipitation and the temperature norm has been carried out. It was found that due to the growing air temperature increase, the average value of the natural inflow within averaging period of 20 years has lost its stability. For the last decade, there has been a steady tendency towards a decrease in inflow share in the water balance of the catchment due to an increase in the share of evaporated moisture.
Key words: reservoir (storage lake), catchment basin, tributaries, underground runoff, surface runoff, climatic norm, long-term annual value, stability of the average, water balance, runoff share, evaporation share, climate change.
