Министерство образования и наук
и Российской Федерации
РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ
ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ (РГГМУ)
На правах рукописи
УДК 556.167 + 556.165
Громова Марина Николаевна
МЕТОДИКА ОЦЕНКИ СТАТИСТИЧЕСКИХ ХАРАКТЕРИСТИК
ГРУППИРОВОК МИНИМАЛЬНОГО СТОКА
Специальность 25.00.27 – гидрология суши, водные ресурсы, гидрохимияАВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата технических наук Санкт-Петербург
Работа выполнена в Российском государственном гидрометеорологическом университете
Научный руководитель: Заслуженный деятель науки РФ доктор технических наук профессор В.В. Коваленко
Научный консультант: кандидат технических наук Н.В. Викторова
Официальные оппоненты: Заслуженный эколог РФ доктор географических наук профессор В.А. Шелутко кандидат технических наук доцент В.Л. Трушевский
Ведущая организация: Государственный гидрологический институт (ГГИ)
Защита диссертация состоится «19» июня 2008 г., в 1530 часов на заседании специализированного совета (Д.212.197.02) Российского государственного гидрометеорологического университета по адресу: г. СанктПетербург, Малоохтинский проспект,
С диссертацией можно ознакомиться в библиотеке Российского государственного гидрометеорологического университета
Автореферат разослан «19» мая 2008 г.
Заслуженный работник высшей школы РФ,
Ученый секретарь специализированного совета, кандидат географических наук Воробьев В.Н.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы. В настоящее время уделяется серьезное внимание возможным тенденциям изменения гидрологического режима под воздействием климатических вариаций и хозяйственной деятельности. Разрабатываются общие и региональные методики долгосрочных прогнозов речного стока. Среди прогнозируемых величин особое место занимают характеристики минимального стока и его многолетних колебаний. Они широко используются в национальном хозяйстве, особенно при водохозяйственном проектировании. Важным является их изучение и оценка в современных условиях возрастающего дефицита водных ресурсов и ухудшения экологического состояния водных систем и отдельных водных объектов в районах интенсивного водопотребления. Особую значимость приобретает постановка задачи о выбросах случайного процесса минимального стока за определенный уровень заданной обеспеченности. С точки зрения практической гидрологии решение этой задачи отвечает на вопрос о группировках лет с водностью ниже некоторого заданного значения (их продолжительности и частоты появления).
Это актуально, например, для гидроэнергетики и водоснабжения. Имея карты подобных характеристик в условиях нового климата, мы фактически прогнозируем статистические характеристики засух.
Цели и задачи исследования. Цель данного исследования – разработка методики оценки характеристик группировок минимального зимнего и летнеосеннего стока России при стационарном гидрометеорологическом режиме и при возможных изменениях климата.
Для достижения поставленной цели решены следующие задачи:
- собраны и статистически обработаны данные о минимальном 30-суточном зимнем и летне-осеннем стоке рек России;
- численно оценены характеристики группировок (частота и длительность), соответствующие стационарному гидрологическому режиму и получен комплект карт с применением ГИС-технологий;
- выполнена параметризация стохастической модели формирования минимального стока для стационарной гидрологической ситуации на территории России;
- выполнена сценарная оценка вероятностных характеристик минимального стока с использованием стохастической модели, основанной на уравнении Фоккера–Планка–Колмогорова (ФПК), при возможном изменении климата и построены соответствующие карты;
- выполнена оценка характеристик группировок минимального стока на различные временные интервалы с учетом изменения климатической ситуации и получен комплект карт возможных их изменений в будущем и карт отклонений значений этих характеристик от таковых, соответствующих стационарному гидрометеорологическому режиму;
- оценены возможные пути использования полученных результатов.
Методика исследований и исходный материал. Решение перечисленных задач проводилось численными методами на основе статистической обработки исходной гидрологической информации. Для расчета характеристик группировок использовался аппарат теории выбросов с помощью программ, написанных на языке Pascal. Решение задачи долгосрочных оценок характеристик минимального стока осуществлялось при помощи математического моделирования на основе уравнения ФПК.
Для построения карт использовались программы Surfer и Idrisi.
Научная обоснованность и достоверность положений и выводов подтверждается статистической оценкой промежуточных результатов, а также использованием теории марковских случайных процессов, приводящих к уравнению Пирсона, описывающему стационарные процессы формирования стока.
Научная новизна и практическая значимость. В ходе проведенного исследования были получены следующие основные результаты:
1. Выполнена параметризация стохастической модели формирования минимального стока на территории России с плотностью 30 км2 на расчетную точку.
2. Произведена оценка характеристик группировок минимального стока, соответствующих стационарному гидрометеорологическому режиму. Построены карты этих характеристик.
3. Выполнена сценарная оценка вероятностных характеристик минимального 30-суточного стока с учетом климатического сценария, при котором происходит повышение температуры воздуха в будущем и почти повсеместное увеличение количества осадков. Построены карты этих характеристик, позволяющие оценивать чувствительность различных отраслей национального хозяйства к возможным изменениям статистических характеристик минимального стока.
4. Выявлены и построены зависимости коэффициента корреляции между смежными рядами минимального стока и такими характеристиками, как модуль минимального стока и коэффициент вариации. Первые использовались для определения коэффициентов r(1) в новой климатической ситуации.
5. Выполнена оценка возможных значений характеристик (частоты и продолжительности) группировок минимального стока с учетом возможных изменений климата. Построены карты этих характеристик.
6. Рекомендованы пути возможного практического применения полученных результатов.
Практическое применение результатов диссертации определяется возможностью использования полученных карт статистических характеристик минимального стока при оценке чувствительности производственных функций различных отраслей экономики к вариациям гидрологической ситуации вследствие изменения климата.
Работа выполнялась в рамках госбюджетной темы кафедры гидрофизики и гидропрогнозов РГГМУ «Моделирование и прогнозирование гидрологических процессов в условиях антропогенного изменения речного стока», а также двух тем («Фрактальная диагностика речного стока для целей частично инфинитного прогнозирования гидроэкологических катастроф» и «Частично инфинитный механизм в моделировании и прогнозировании гидрологических катастроф»), выполняемых по заданию федерального агенства по образованию на проведение научных исследований (№№ гос. рег. 01.2004 09546, 01.2006 03264). Ее результаты внедрены в учебный процесс по специальности «Гидрология» – 07.32.00.
На защиту диссертации выносятся:
- результаты оценок характеристик группировок минимального 30-суточного летне-осеннего и зимнего стока, соответствующих стационарному гидрометеорологическому режиму;
- результаты сценарных оценок вероятностных характеристик минимального 30-суточного стока для территории Сибири и Дальнего Востока по одному из возможных климатических сценариев;
- закартированные результаты сценарных оценок частоты и продолжительности выбросов случайного процесса минимального стока для территории России;
- методика долгосрочной сценарной оценки характеристик группировок минимального стока при возможном изменении климата (для равновесных климатических сценариев).
Апробация работы. Основные положения диссертации докладывались на научных семинарах кафедры гидрофизики и гидропрогнозов, итоговых сессиях ученого Совета РГГМУ в 2005, 2006 г.г., на конференции молодых ученых «Изменение климата и окружающая среда» в 2005 г., на VI международном молодежном экологическом форуме стран Балтийского региона «Экобалтика» в 2006 г., на политехническом симпозиуме «Молодые ученые – промышленности Северо-Западного региона» в 2006 г., на дне молодого ученого в Политехническом университете в 2007 г.
По теме диссертации опубликовано 9 работ.
Структура и объем работы. Диссертация состоит из введения, трех глав, заключения, списка литературы, содержащего 81 источника, 4 приложения. Работа изложена на 165 страницах текста, включая 32 рисунка и таблиц.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении излагается актуальность темы, а также основные цели и задачи диссертации.
В первой главе рассматривается влияние основных физикогеографических факторов на формирование и распределение по территории России минимального стока. Оцениваются возможные изменения основных факторов, формирующих минимальный сток. Приводятся существующие методы оценки гидрологических последствий изменений климата. Описываются существующие прототипы данной методики. Ставятся задачи исследования.
Во второй главе описываются исходные данные и математический аппарат, который лежит в основе разрабатываемой методики.
Экспериментальными данными для расчета и прогноза характеристик группировок служили ряды наблюдений за минимальным 30-суточным зимним и летне-осенним стоком на территории России, приведенные в изданиях Государственного водного Кадастра (Многолетние данные о режиме и ресурсах поверхностных вод суши). При этом использовалась информация по створам на реках с водосборной площадью в пределах от 1200 до 50000 км2 в зависимости от территориального расположения и рельефа местности. Рассматривались реки, сток которых не искажен влиянием антропогенной деятельности. Для расчетов было отобрано 365 рядов наблюдений за минимальным 30-дневным зимним стоком, и 450 рядов – за минимальным 30-дневным летне-осенним. Средняя продолжительность рядов наблюдений составила лет. После статистической обработки выяснилось, что из всех исходных гидрологических рядов репрезентативными оказалось для зимнего стока ряда, для летне-осеннего – 421.
Основной моделью в диссертации служила стохастическая модель формирования стока Фоккера–Планка–Колмогорова (ФПК). Исследование возможности прогнозирования по этой модели характеристик минимального стока было выполнено Н. В. Викторовой в 2002 г. Однако эти исследования остались логически не полностью завершенными, так как, во-первых, рассматривалось только несколько регионов Европейской территории России, а, во-вторых, были опущены вопросы, связанные с прогнозом характеристик выбросов, что является основой данной диссертационной работы.
Уравнение ФПК имеет вид:
где p (Q, t) – плотность вероятности расхода воды;
А(Q, t) и В(Q, t) – соответственно коэффициенты сноса и диффузии, где G c – интенсивность внутреннего шума (~ ) ;
GN – интенсивность внешнего климатического шума N ;
G~N – взаимная интенсивность шумов;
c и N – соответственно статистические нормы бассейновых и климатических шумов.
В установившемся режиме (а климатические сценарии предполагают равновесную ситуацию), уравнение ФПК переходит в уравнение Пирсона где а, b0, b1 и b2 – параметры, которые характеризуют физико-статистические свойства водосбора и рассчитываются по формулам:
Процедура прогнозирования состоит из следующих этапов:
1) параметризация модели, то есть определение численных значений коэффициентов, соответствующих текущей климатической ситуации и уровню хозяйственной деятельности;
2) задание внешнего воздействия (приземной температуры воздуха и нормы годовых осадков из заданного климатического сценария);
3) задание трендов хозяйственной деятельности, демографии и т. д. с помощью коэффициентов, входящих в модель;
4) численное интегрирование уравнения ФПК (прогноз характеристик стока на требуемую заблаговременность).
Обоснованность применимости этой модели объясняется тем, что она подходит для всех видов многолетнего стока, в частности для минимального, так как ее решением в равновесном (квазистационарном) режиме является семейство кривых Пирсона, которыми (как правило) исчерпываются запросы практической гидрологии.
При оценке характеристик группировок минимального стока использовался аппарат теории случайных выбросов. Ряд минимального стока рассматривался как случайная последовательность, так как для него характерны иррегулярные колебания и дискретное задание данных.
условию стационарности. Это означает, что любая m-мерная плотность распределения вероятности px xi1, xi2,..., xim зависит только от взаимного расположения моментов времени i1, i2,..., im и не меняется при сдвиге всех отсчетов на одни и те же интервалы времени :
Считается, что в i-тый момент времени имеет место выброс случайной последовательности {x k } сверху вниз (отрицательный выброс), если В данном случае за значение а принят расход воды в реке, ниже которого не могут быть полностью удовлетворены промышленно-бытовые или какие-либо другие требования (минимальный 30-суточный расход 80 %-ной обеспеченности). В этом случае число выбросов сверху вниз и их продолжительность показывают, сколько раз за интересующий нас промежуток времени будут наблюдаться перебои с водоснабжением и какая будет их длительность. Если данные показатели выше допустимых, то встает вопрос о дополнительных источниках водоснабжения.
В данной работе рассчитывались основные характеристики группировок, такие как средняя их длительность и частота появления.
Частота группировок рассчитывается по формуле где P{X 0 a, X 1 > a} находится по зависимости где рx (x1, x2 ) двумерная плотность распределения смежных членов последовательности {X k }= 0.
Средняя длительность отрицательной группировки или В работе В.А. Румянцева и И.В. Бовыкина было предложено осуществлять предварительную нормализацию исходного ряда {xi }. Нормализованные значения Yi удобно получать с помощью преобразования, которое нашло широкое применение при использовании метода Монте-Карло:
где F(xi ) – одномерная функция распределения исходного ряда;
Ф- 1 – функция, обратная функции распределения случайной величины, распределенной по нормальному закону с параметрами ноль и Преобразованный таким образом ряд значений {Yi } находится во взаимно однозначном соответствии с исходным рядом {Xi } и подчиняется одномерному закону распределения N (0, 1). Если взять величину b, удовлетворяющую условию то каждой группировке значений минимального 30-дневного стока {Xi } относительно уровня а будет однозначно соответствовать аналогичная группировка значений {Yi } относительно уровня b.
Расчет характеристик группировок минимального стока подразумевает выполнение следующих этапов:
1) определение коэффициента вариации исходного гидрологического ряда;
2) расчет коэффициента автокорреляции между смежными членами ряда минимального стока;
3) осуществление нормализации исходного асимметричного ряда и определение коэффициента корреляции между смежными членами нормальной нормированной последовательности;
4) определение характеристик группировок по соответствующим таблицам в зависимости от коэффициента автокорреляции нормализованного ряда и уровня, ниже которого и рассчитываются выбросы минимального стока.
Для оценки изменения характеристик группировок необходимо иметь информацию об изменении значений коэффициентов вариации и автокорреляции r(1) в новых климатических условиях.
Для частоты появления группировок и средней продолжительности отрицательных группировок Еl- входными параметрами для таблиц служат уровень b и коэффициент корреляции между смежными членами нормированной последовательности {yi}. Как таковой методики для прогноза коэффициента корреляции между смежными членами ряда минимального стока найдено не было. В этом случае целесообразно применять графические зависимости и аналитические выражения, связывающие величину коэффициента корреляции между смежными членами ряда и такими характеристиками стока, как среднемноголетний модуль М, коэффициент вариации и т. д. Последние прогнозируются при помощи стохастической модели формирования стока ФПК. Используя найденное значение коэффициента корреляции нормированного ряда, соответствующего новой климатической ситуации, при заданном уровне водности b, определяется частота появлений и средняя длительность группировок минимального 30-дневного стока в будущем при изменении водности.
В третьей главе рассчитываются основные характеристики группировок минимального стока, соответствующие стационарной гидрометеорологической ситуации и дается сценарная оценка статистических характеристик и характеристик группировок на различные временные интервалы с учетом наиболее вероятного климатического сценария. Для сценарной оценки выбросов необходима информация о возможном изменении коэффициента корреляции между смежными членами ряда в будущем. В диссертации был сделан новый шаг в прогнозировании коэффициента автокорреляции с использованием связи, существующей между коэффициентом автокорреляции, модулем минимального стока и коэффициентом вариации. Подобная связь была установлена по имеющимся рядам наблюдений. Полученное поле точек оказалось достаточно хаотичным, однако после разбиения его на определенные градации, удалось получить статистически значимые зависимости (рисунок 1).
Рисунок 1 – Зависимость r(1) = f(M) для минимального 30-суточного зимнего (а) и летне-осеннего (б) стока на территории России Полученные результаты представлены в виде карт вероятностных характеристик минимального стока и карт отклонений характеристик группировок минимального стока в новой климатической ситуации от тех, которые соответствуют стационарному гидрометеорологическому режиму (рисунки 2, 3).
Рисунок 2 – Изолинии возможных отклонений частоты группировок минимального 30-суточного зимнего стока от существующих к 2100 году, % Рисунок 3 – Изолинии возможных отклонений длительности группировок минимального 30-суточного летне-осеннего стока от существующих На основании полученных результатов, можно сделать следующие выводы. Возможное потепление на Земле в значительной степени скажется на вероятностных характеристиках минимального стока России. К 2100 году их изменение составит 20 %. Что касается характеристик выбросов, то влияние изменения климата на них при использованном сценарии скажется меньше. С практической точки зрения, результаты исследований могут быть применены при составлении водохозяйственных балансов на перспективу, при усовершенствовании методик экологического использования речных бассейнов.
В заключении сформулированы основные результаты исследований:
1. Впервые произведена оценка характеристик группировок минимального 30-суточного зимнего и летне-осеннего стока, соответствующих стационарному гидрометеорологическому режиму и построены карты этих характеристик для большей части территории России;
2. Получена сценарная оценка вероятностных характеристик минимального 30-суточного стока при одном из возможных климатических сценариев для территории Сибири и Дальнего Востока, а также построены карты этих характеристик;
3. Разработана методика прогнозирования коэффициента автокорреляции на основе выявленной связи этого коэффициента с модулем минимального стока и коэффициентом вариации;
4. Выполнена сценарная оценка частоты и продолжительности группировок минимального стока для территории России и построены соответствующие карты.
В приложении представлены результаты расчетов, а также справка об эффективности результатов исследований.
По теме диссертации опубликованы следующие работы:
1. Оценка влияния возможных изменений минимального стока на экономическое развитие северо-запада ЕТР // Материалы итоговой сессии ученого совета 27 – 28 января 2004 года. Информационные материалы. Часть I.
Секции метеорологии, гидрологии, экономических и социальногуманитарных наук. – СПб.: изд. РГГМУ, 2004. – С. 66 – 67 (в соавторстве с Н.В. Викторовой);
2. Влияние возможных изменений минимального стока на экономическое развитие Северо-Запада ЕТР // Проблемы современной гидрологии.
Сборник научных трудов РГГМУ. – СПб: изд. РГГМУ, 2004. – С. 45 – 53 (в соавторстве с Н.В. Викторовой);
3. Прогнозы стока на территории Северо-Запада ЕТР // Материалы итоговой сессии ученого совета 25 – 26 января 2005 года. Информационные материалы. – СПб.: изд. РГГМУ, 2005. – С. 8 – 9. (в соавторстве с С.М. Губановым) 4. Оценка характеристик выбросов минимального стока ЕТР.
// Сборник трудов международной школы-конференции молодых ученых, аспирантов и студентов «Изменение климата и окружающая среда». – СПб.:
Изд. «Гранд», 2005. – С. 177 –183.
5. Оценка и прогноз характеристик выбросов минимального стока рек ЕТР // Тезисы докладов итоговой сессии ученого совета 25 – 26 января года. – СПб.: изд. РГГМУ, 2006 – С. 12 – 13.
6. Application of stochastic model for the forecast of the minimal runoff of the rivers of basin of the Baltic sea // Сборник трудов VI международного молодежного экологического форума стран Балтийского региона «Экобалтика».
– СПб: Изд. СПбГПУ, 2006. – С. 139 (в соавторстве с Н.В. Викторовой).
7. Прогноз характеристик минимального стока для целей управления водными ресурсами // Материалы конференций Политехнического симпозиума «Молодые ученые – промышленности Северо-Западного региона». – СПб.: Изд. СПбГПУ, 2006. – С. 13 (в соавторстве с Н.В. Викторовой).
8. Влияние нормы модуля и коэффициента вариации минимального Электронный журнал "Исследовано в России", 230, стр. 2199-2205, 2006 г.
http://zhurnal.ape.relarn.ru/articles/2006/230.pdf (в соавторстве с Н.В. Викторовой).
В рецензируемых журналах из списка ВАК 9. Прогноз характеристик минимального стока для целей управления водными ресурсами // Научно-технические ведомости СПбГПУ. – СПб.: Изд.
СПбГПУ, 2007. – № 2 (50). – С. 284 – 287 (в соавторстве с Н.В. Викторовой).