LONG-TERM FORECASTING METHODOLOGY FOR WATER AND HYDROCHEMICAL RUNOFF OF MOUNTAIN RIVERS Yury B. Kirsta1, Alexander V. Puzanov2 1 Institute for Water and Environmental Problems, Siberian Branch, Russian Academy of Sciences, Barnaul, Russia 2 Polzunov Altai State Technical University, Barnaul, Russia Abstract We propose the method of system-analytical modeling in combination with GIS to construct prognostic simulation models of complex natural systems. The structure and parameters of models are determined via the inverse problem solution for system of equations by optimization technique. Using mid-size and small rivers of Altai-Sayan mountain country as a case study, we have developed the complex of universal high quality models for spatial generalization of air temperature and precipitation, and calculating the seasonal and long-term dynamics of water runoff and hydrochemical runoff of rivers. Keywords: system-analytical modeling, river runoff, hydrochemical runoff, mountain rivers, Altai, Sayan МЕТОДОЛОГИЯ ДОЛГОСРОЧНОГО ПРОГНОЗА ВОДНОГО И ГИДРОХИМИЧЕСКОГО СТОКОВ ГОРНЫХ РЕК Кирста Ю.Б.(1)(2), Пузанов А.В.(1) 1 Институт водных и экологических проблем СО РАН, г. Барнаул 2 Алтайский государственный политехнический университет, г. Барнаул Предлагается совмещенный с ГИС метод системно-аналитического моделирования для создания прогнозных имитационных моделей сложноорганизованных природных систем. Структура и параметры моделей определяются путем решения обратной задачи для систем уравнений оптимизационными мето- дами. На примере 34 средних и малых рек Алтае-Саянской горной страны разработан комплекс уни- версальных моделей высокого качества: пространственного обобщения температур воздуха и осадков, расчета сезонной и межгодовой динамики водного и гидрохимического стоков рек. Ключевые слова: системно-аналитическое моделирование, речной сток, гидрохимический сток, горные реки, Алтай, Саяны Введение. Прогноз водного и гидрохимического стоков рек является важнейшей состав- ляющей мониторинга водных ресурсов. Для прогноза стоков разработаны различные методы, в том числе учитывающие физико-географические, гидрографические и гидрохимические осо- бенности речных бассейнов. Серьезную проблему здесь представляют горные территории, имеющие сложную орографию и пеструю пространственно-временную структуру климатиче- ских полей. Современные динамические модели водного и/или гидрохимического стоков рек как правило объединяются в одно целое с ГИС-средствами. ГИС используются для учета про- странственно распределенных входных данных, визуализации результатов расчетов и других целей. Обычно модели основаны на дифференциальных уравнениях гидродинамики и мате- матической физики, дополняемых имитационными зависимостями [1, 2]. Дифференциальные уравнения требуют подробной пространственно распределенной информации об осадках, тем- пературах воздуха, испаряемости, подземных водоносных горизонтах, морфометрии, уклонах, свойствах водонасыщенных почв и пород, показателях сопротивлению течения на склонах и в русловой сети, и т.д., и т.п. В условиях гор подобная информация практически отсутствует, что значительно снижает прогнозную ценность таких моделей. Метод системно-аналитического моделирования. Нами предлагается метод си- стемно-аналитического моделирования (САМ), совмещенный с ГИС и обеспечивающий по- строение эффективных прогнозных моделей водного и гидрохимического стоков горных рек [3–6]. Основой методологии САМ является определение структуры и параметров имита- ционных моделей сложноорганизованных природных систем путем решения обратной мате- матической задачи оптимизационными методами. Модели представляют собой системы ба- лансовых уравнений, решаемых в программной среде MATLAB, и дают проверяемое адекват- ное описании основных физических, гидрологических и гидрохимических процессов, форми- рующих указанные стоки. Эти уравнения обеспечивают значительно меньшую дисперсию расхождений между прогнозными и фактическими значениями водного и гидрохимического стоков по сравнению с прогнозом по их среднемноголетним значениям. В нашем случае САМ используется для поиска и количественной характеристики функ- циональных связей водного и гидрохимического стоков с метеорологическими факторами, морфометрией и ландшафтной структурой речных бассейнов. Имитационные модели водного и гидрохимического стоков строятся из алгебраических уравнений, поскольку, как указано выше, применение дифференциальных уравнений для описания гидрологических и гидрохи- мических процессов в условиях гор крайне затруднительно. Определение характера функцио- нальных связей между процессами и факторами производится через гидрологически и гидро- химически обоснованные подбор и корректировку системы уравнений для минимизации рас- хождений (квадратичной невязки) между рассчитываемыми и наблюдаемыми динамическими характеристиками. Использование в качестве левых частей уравнений наблюдаемых значений 117 водного или гидрохимического стока позволяет через решение обратной задачи найти пара- метры уравнений и оценить квадратичную невязку для тестируемого варианта модели. Мо- дель считается построенной, когда найдена такая система уравнений, которая обеспечивает наименьшее значение невязки. В среде программирования MATLAB нами создан комплекс программ, обеспечивающий выбраковку недостоверных экспериментальных данных и иден- тификацию более 100 параметров проверяемых моделей через оптимизационное решение си- стем из 1500 алгебраических уравнений. В САМ также отсутствует первоначальная фиксация формы искомых зависимостей от факторов среды в ходе решения обратной задачи. Отметим, что в случае дифференциальных уравнений эта форма жестко фиксирована самим выбором тех или иных уравнений. С этой целью нами применяется пробная функция H, задаваемая выражением: 𝑌1 + 𝑍1(𝑋 − 𝑋1), если 𝑋 < 𝑋1 (𝑋 − 𝑋1) + 𝑌1, если {𝑋1 ≤ 𝑋 < 𝑋2 , 𝑌2−𝑌1 𝐻(𝑋1, 𝑋2, 𝑌1, 𝑌2, 𝑍1, 𝑍2, 𝑋) = { (1) 𝑋2−𝑋1 𝑋1 ≠ 𝑋2 𝑌2 + 𝑍2(𝑋 − 𝑋2), если 𝑋 ≥ 𝑋2 где X1, X2, Y1, Y2, Z1, Z2 – параметры, определяемые в процессе решения обратной задачи; X – какая-либо переменная модели. Функция H является непрерывной кусочно-линейной функ- цией из трех произвольных линейных фрагментов и позволяет аппроксимировать широкий спектр различных зависимостей между переменными и факторами среды путем изменения значений своих параметров (рис. 1). (X2,Y2) Y= H (X1,Y1) X Рис. 1. Непрерывная кусочно-линейная функция Y= H(X1,X2,Y1,Y2,Z1,Z2,X) из трех линейных фрагментов с произвольно меняемыми параметрами (см. уравнение (1)). Модели требуют проверки их адекватности. Нами предложен удобный универсальный критерий, позволяющий оценивать адекватность любых расчетных методов и моделей путем сравнения рассчитанных и наблюдаемых значений характеристики [7]: 𝐴 = 𝑆разн ⁄√2𝑆набл , (2) где A – критерий адекватности; Sразн – стандартное (среднеквадратичное) отклонение невязки (разности сравниваемых расчетного и наблюдаемого рядов) модели, Sнабл – стандартное откло- нение для наблюдаемого ряда, 1/√2 – множитель. Согласно (2), критерий A представляет собой погрешность модели, нормированную на стандартное отклонение данных наблюдений. Интервал значений A=0÷0.71 характеризует раз- личную степень адекватности/идентичности расчетных и наблюдаемых значений переменной с их наилучшим совпадением при A~0. Критерий A подобен известным показателям качества моделей: RSR (RMSE-Standard deviation Ratio [8, 9]) и NSE (Nash-Sutcliffe model Efficiency co- efficient – коэффициент Нэша–Сатклиффа [9]), с которыми связан зависимостями RSR=𝐴√2 и NSE = 1–RSR2 = 1–2A2. В САМ рассчитывается также чувствительность моделей к вариациям входных факто- ров, выраженная как вклад конкретного фактора в дисперсию наблюдаемых значений водного или гидрохимического стока [10]. Чувствительность является важным показателем, который необходим для эффективного управления водными ресурсам. По критерию A и чувствитель- ности осуществляется полный компонентный анализ дисперсии невязки моделей (дисперсии 118 разностей между рассчитанными и наблюдаемыми значениями выходных переменных). Такой анализ позволяет адекватно рассчитать критерии RSR и NSE и оценить качество прогноза. Исходные материалы. Исследования выполнены на примере Алтае-Саянской горной страны, являющейся часть мирового водораздела между гумидной областью Северного Ледо- витого океана и аридной бессточной областью Центральной Азии. Для всех рек страны наибольший водный сток (до 80-90% годового) отвечает теплому периоду года. Режим вод- ного стока зависит, в основном, от таяния снега весной и количества осадков летом и осенью. За счет таяния снега обычно формируется более 50% годового стока. САМ водного и гидрохимического стоков выполнялось для 34 бассейнов средних и ма- лых рек указанной страны. Для учета ландшафтной структуры бассейнов выполнена типиза- ция ландшафтных выделов, отражающая условия формирования водного и гидрохимического речных стоков, в том числе высотно-поясную и структурно-ярусную неоднородность терри- тории. Нами выделено 12 типологических групп геосистем и отдельно 13-я для аквальных ландшафтов, имеющих незначительную площадь [11]. Для характеристики гидрохимического стока были выбраны его следующие компоненты: три минеральные формы азота NO− − 2 , NO3 , + 3− NH4 , фосфаты PO4 , общее растворенное железо, минерализация вод, взвешенные вещества. Годовая и межгодовая динамика водного и гидрохимического стоков рассчитывалась для 1951–2003 гг. по 4 гидрологическим периодам/сезонам: первый (зимняя межень, XII-III ме- сяцы), второй (весенне-летнее половодье, IV-VI), третий (летняя межень, VII-VIII), четвертый (осенняя межень с возможными паводками при сильных дождях, IX-XI). Разработка моделей осуществлялась по данным Управления по гидрометеорологии и мо- ниторингу окружающей среды о водном (~5000 значений) и гидрохимическом (~10000) сто- ках, осадках и температурах воздуха (11 реперных метеостанций). Температуры воздуха и осадки характеризовались нами через доли/проценты от их среднемноголетних месячных зна- чений [7]. Помимо указанных данных использовалась картографическая информация о рель- ефе и ландшафтной структуре бассейнов, в том числе площадь и средняя высота каждого ланд- шафта в каждом бассейне, площади пашни, уклоны бассейнов и высоты их замыкающих ство- ров, длины речных русел (от истока до створа). Весь картографический материал обрабаты- вался в программной среде ArcGIS 9.2 с модулем 3D Analyst. Результаты моделирования. Нами разработан комплекс универсальных балансовых моделей сезонной и многолетней динамики водного и гидрохимического стоков для средних и малых рек страны. Комплекс включает климатическую модель, модель водного стока и семь моделей гидрохимического стока (NO− − + 3− 2 , NO3 , NH4 , PO4 , ионов, общего растворенного же- леза, взвешенного вещества). Полученные для всех моделей стоков критерии RSR < 0.60 и NSE > 0.65 означают их хорошее и очень хорошее качество (см., например, [10]). Климатическая модель отражает единообразие внутригодовой и многолетней динамики месячных осадков и среднемесячных температур воздуха по всей территории Алтае-Саянской горной страны при условии нормировки осадков на их среднемноголетнее значение за июль «in situ», а температур – на их среднемноголетнее значение за январь (для X-IV месяцев) и июль (для V-IX) «in situ». Такие пространственно обобщенные нормированные характери- стики уже не зависели от координат или высоты расположения характеризуемых участков, были одинаковы для всех речных бассейнов и, поэтому, применялись в дальнейших расчетах как адекватные площадные характеристики метеорологических полей указанной страны. Модель водного стока представлена балансовым уравнением, учитывающим осадки, температуры воздуха, площадь и высоту ландшафтов: 𝑄 𝑖 = ∑ {𝑎𝑘 𝑆𝑘𝑖 𝑃1 𝐻(𝑐1 , 𝑐1 , 1,1, 𝑐2 , 𝑐3 , 𝑇1 ) 𝐻(𝑐4 , 𝑐4 , 1,1, 𝑐5 , 𝑐6 , ℎ𝑘𝑖 )} + 𝑘 + ∑𝑘{𝑏𝑘 𝑆𝑘𝑖 𝑃2 𝐻(𝑐7 , 𝑐7 , 1,1, 𝑐8 , 𝑐9 , 𝑇2 ) 𝐻(𝑐4 , 𝑐4 , 1,1, 𝑐5 , 𝑐6 , ℎ𝑘𝑖 )} + 𝑐10 (3) 𝑖 где 𝑄 – нормированный на среднемноголетнее значение сезонного стока для замыкающего створа бассейна i, i =1,2,…, 34; первое и второе слагаемые уравнения отвечают вкладам пред- шествующего и текущего гидрологического сезона соответственно, а при расчете 1-го сезона (зимняя межень) – отвечают 3- и 4-му сезонам предшествующего года; ak, bk – параметры, 119 характеризующие вклад k-го ландшафта, k =1–13; 𝑆𝑘𝑖 , ℎ𝑘𝑖 – относительная площадь и высота k-го ландшафта в бассейне i; c1–c9 – параметры, отражающие влияние температур 𝑇1 , 𝑇2 и вы- соты ℎ𝑘𝑖 ; 𝑃1 , 𝑃2 – нормированные месячные осадки в среднем за сезон; 𝑇1 , 𝑇2 – отклонения нормированных температур от 1 в среднем за сезон; c10 – параметр для постоянного пополне- ния (c10>0) или потерь (c10<0) стока в грунтовые воды и воды зон трещиноватых пород (под- земный сток); H –кусочно-линейная функция (1). Все параметры уравнения (3) определяются с помощью оптимизационных методов пакета программ MATLAB по имеющимся ~1300 наблюдаемым значениям стока за соответствующий гидрологический сезон. Модели гидрохимического стока представлены двумя уравнениями, учитывающими рассчитанные по (3) водные стоки, осадки, поперечный уклон бассейна, площадь пашни: для 1-го гидрологического сезона (зимняя межень) – 𝐻𝑅 𝑖 = ∑𝑘{𝑎𝑘 𝑄𝑘𝑖 𝐻(𝑐1 , 𝑐1 , 1,1, 𝑐2 , 𝑐3 , 𝑃) 𝐻(𝑐4 , 𝑐4 , 1,1, 𝑐5 , 𝑐6 , 𝐾 𝑖 )} +b 𝑞 𝑖 + d 𝑆 𝑖 𝑄 𝑖 , (4а) для остальных сезонов – 𝐻𝑅 𝑖 = ∑𝑘{𝑎𝑘 𝑄𝑘𝑖 𝐻(𝑐1 , 𝑐1 , 1,1, 𝑐2 , 𝑐3 , 𝑃) 𝐻(𝑐4 , 𝑐4 , 1,1, 𝑐5 , 𝑐6 , 𝐾 𝑖 )} +b 𝑞 𝑖 + d √𝑆 𝑖 𝑄 𝑖 , (4б) где 𝐻𝑅 𝑖 – сток соответствующего компонента гидрохимического стока; 𝑃 – среднесезонные нормированные месячные осадки (за предыдущие IX-XI мес. для 1-го сезона или за IV-VI, VII- VIII, IX-XI мес. для 2, 3, 4-го сезонов); 𝑄𝑘𝑖 – расчетный водный сток с k-го ландшафта бассейна i, k =1–13, i =1–34; 𝑎𝑘 – параметры, соответствующие постоянной сезонной концентрации вещества в расчетном водном стоке k-го ландшафта; 𝐾 𝑖 – средний поперечный уклон бассейна i, рассчитываемый как средний тангенс угла наклона склонов относительно горизонтали [6]; c1–c6 – параметры, отражающие влияние на сток вещества осадков P и уклона K i ; b – пара- метр, характеризующий сезонную концентрацию соответствующего вещества в расчетном приходящем (или уходящем) сезонном подземном водном стоке 𝑞 𝑖 бассейна i; 𝑆 𝑖 – относи- тельная площадь пахотных земель в бассейне i; d – параметр, характеризующий зависимость концентрации вещества от площади 𝑆 𝑖 в расчетном водном стоке 𝑄 𝑖 бассейна i. Каждое исследуемое вещество в конкретном сезоне характеризуется своими значениями параметров. Все параметры найдены путем решения обратной задачи для систем (4а) или (4б) из 1200–1500 уравнений, характеризующих сезонные стоки исследуемого вещества за отдель- ные годы с каждого речного бассейна. В итоге уравнения (4а) и (4б) позволяют рассчитывать стоки вещества с отдельных ландшафтов и пашни, его общий сток с речного бассейна и кон- центрацию этого вещества в водном стоке реки. Заключение. По моделям водного и гидрохимического стоков легко выполняется их прогноз на следующий гидрологический сезон при подстановке в (3), (4а), (4б) среднемного- летних значений метеорологических факторов этого сезона вместо фактических. Тем самым стоки предсказываются за 3–4 месяца вперед, причем с вдвое уменьшенной дисперсией не- вязки по сравнению с тривиальным прогнозом по их среднемноголетнему значению. Помимо экологического и водно-ресурсного мониторинга прогнозы важны для управле- ния рисками. Например, прогноз водного стока необходим для регуляции попусков воды из водохранилищ горных ГЭС в период весенне-летнего половодья, когда в водохранилища по- ступает основной объем годового стока. Объединенная с прогнозом оценка чувствительности водного и гидрохимического стоков к изменениям факторов среды (ландшафтной структуре, площадям пашни и др. [12]) может использоваться для практической регуляции качества реч- ных вод через соответствующее управление землепользованием. Разработанный комплекс моделей позволяет количественно характеризовать и прогно- зировать сезонную и многолетнюю динамику водного режима и режима качества речных вод для произвольно выбранного речного бассейна Алтае-Саянской горной страны. Для иденти- фикации моделей требуются цифровые картографические данные о ландшафтной структуре бассейна, данные о температурах воздуха и осадках, а также 1–2 года наблюдений за характе- ризуемым стоком. Модели могут быть применены к любой горной территории после ее ланд- шафтной типизации с новой оценкой параметров через решение обратной задачи. 120 ЛИТЕРАТУРА [1] Виноградов Ю.Б., Виноградова Т.А. Математическое моделирование в гидрологии: Учебное по- собие для студентов учреждений высш. проф. образования. М.: Издательский центр «Академия», 2010. 304 с. [2] Румынин В.Г. Геомиграционные модели в гидрологии. СПб: Наука, 2011. 1158 с. [3] Kirsta Yu.B. System-analytical modelling – Part I: General principles and theoretically best accuracies of ecological models. Soil-moisture exchange in agroecosystems // Ecol. Modelling. 2006. Vol. 191. P. 315–330. [4] Кирста Ю.Б., Кирста Б.Ю. Информационно-физический закон построения эволюционных си- стем. Системно-аналитическое моделирование экосистем. Барнаул: Изд-во Алт. гос. ун-та, 2014. 283 с. [5] Кирста Ю.Б., Пузанов А.В., Ловцкая О.В., Лубенец Л.Ф., Кузняк Я.Э., Пахотнова А.Ю. Имита- ционная математическая модель стока средних и малых рек для горных территорий // Известия Самарского научного центра РАН. 2012. Т.14, №1(9). С. 2334-2342. [6] Кирста Ю.Б., Пузанов А.В. Моделирование гидрохимического стока горных рек: 1. Модель стока фосфатов // Вестник Алтайского государственного аграрного университета. 2016. №6 (140). С. 78-84. [7] Кирста Ю.Б. Пространственное обобщение климатических характеристик для горных террито- рий // Мир науки, культуры, образования. 2011. № 3 (28). С. 330-337. [8] Moriasi D.N., Arnold J.G., Van Liew M.W., Bingner R.L., Harmel R.D., Veith T.L. Model evaluation guidelines for systematic quantification of accuracy in watershed simulation // Transactions of the ASABE. 2007. Vol. 50(3). P. 885-900. [9] Koch·M., Cherie N. SWAT-modeling of the impact of future climate change on the hydrology and the water resources in the upper blue Nile river basin, Ethiopia // Proceedings of the 6th International Conference on Water Resources and Environment Research, ICWRER 2013. Koblenz, Germany, June 3-7, 2013. P. 428- 523. [10] Кирста Ю.Б. Моделирование гидрохимического стока горных рек: 2. Оценка качества модели стока фосфатов // Вестник Алтайского государственного аграрного университета. 2016. №6 (140). С. 85-91. [11] Кирста Ю.Б., Лубенец Л.Ф., Черных Д.В. Типизация ландшафтов для оценки речного стока в Алтае-Саянской горной стране // Устойчивое развитие горных территорий. 2011. №2(8). С. 51- 56. [12] Kirsta Yu.B., Puzanov A.V. System-analytical modeling of water quality for mountain river runoff // Вод- ный форум БРИКС: Международная научно-практическая конференция, Москва, 29-30 сент. 2016 г. http://bricswater.org/data/2016/10/20/1107835601/Секция%204_%20Целый.pdf (дата обраще- ния 9.06.2017). 121