Размещено на http://www.allbest.ru/

1

Размещено на http://www.allbest.ru/

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

кандидата физико-математических наук

Моделирование и прогноз возникновения паводковых ситуаций в руслах горно-равнинных рек

05.13.18 - Математическое моделирование, численные методы и комплексы программ

Вандина Наталья Валерьевна

Краснодар - 2012

Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Ставропольский государственный университет»

Научный руководитель: доктор физико-математических наук, профессор СЕМЕНЧИН Евгений Андреевич

Официальные оппоненты:

ЛЕБЕДЕВ Константин Андреевич, доктор физико-математических наук, профессор кафедры прикладной математики ФГБОУ ВПО «Кубанский государственный университет»

УГОЛЬНИЦКИЙ Геннадий Анатольевич доктор физико-математических наук, профессор, заведующий кафедрой прикладной математики и программирования ФГАОУ ВПО «Южный федеральный университет»

Ведущая организация:

ФГБОУ ВПО «Северо-Кавказский государственный технический университет»

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы исследования. В настоящее время остро стоит проблема поиска методов и моделей прогноза возникновения паводковых ситуаций в руслах горно-равнинных рек, наносящих значительный экономических ущерб регионам России. Например, на реке Кубань среднемноголетний ежегодный ущерб от наводнений составляет 2,1 млрд. рублей, а в последние годы наблюдается тенденция к его увеличению.

Проблемам математического моделирования в гидрологии посвящены многочисленные исследования как в нашей стране (Кучмент Л. С., Гельфан А. Н., Корень В. И., Карасев И. Ф., Бураков Д. А., Демидов В. Н., Мотовилов Ю. Г., Найденов В. И. и др.), так и за рубежом (Кюнж Ж. А., Холли Ф. М., Прейсман А., Эббот М. Б., Тсенг, Фрэд Д. Л. и др.).

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

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

Объект и предмет исследования. Объектом исследования является процесс движения потока воды в русле реки. Предметом исследования - процесс возникновения паводковых ситуаций в руслах горно-равнинных рек.

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

Для достижения поставленной цели необходимо решить следующие задачи:

1. Разработать алгоритмы решения (аналитическими и численными методами) моделей движения потока воды, основанные на преобразовании исходной классической математической модели неустановившегося движения потока воды в русле реки. Реализовать предложенные алгоритмы в программных продуктах для ЭВМ.

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

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

4. Применить предложенные модели и методики для расчета и прогноза возникновения паводковых ситуаций на различных участках русла реки Кубань.

Научная новизна диссертации заключается в следующем:

1. Разработан алгоритм численного решения методом расщепления системы дифференциальных уравнений, описывающей нестационарное движение потока воды в русле горно-равнинной реки при характерных для данных рек (больших) величинах уклона дна русла. Данный алгоритм реализован в программном продукте «Calculation of flow characteristics (COFC)», зарегистрированном в Федеральной службе по интеллектуальной собственности, патентам и товарным знакам.

2. Построено аналитическое решение уравнения, описывающего нестационарное движение воды на участке русла горно-равнинной реки при малых глубинах потока.

3. Построено аналитическое решение уравнения, моделирующего движение воды в гидрометрическом створе русла реки произвольной глубины при заданной в этом уравнении зависимости глубины потока от времени.

4. Предложена новая математическая модель прогноза возникновения паводковой ситуации на заданном участке русла реки при больших уклонах дна русла.

5. Предложена новая математическая модель прогноза возникновения паводковой ситуации в гидрометрическом створе русла реки.

6. Предложена новая стохастическая модель прогноза возникновения паводковой ситуации в створе русла реки, учитывающая помехи в результатах измерений характеристик потока. Данная модель основана на построенной в работе оптимальной в среднеквадратическом смысле оценке расхода воды.

7. Разработана методика расчета средней скорости потока воды в русле при неравномерном установившемся движении.

8. Предложена уточненная методика расчета расхода потока на участке русла горно-равнинной реки, основанная на использовании уравнения руслового водного баланса.

9. С помощью геоинформационных систем и программ статистической обработки данных усовершенствована методика исследования морфометрических характеристик русла реки.

10. Предложенные в работе методики и модели применены для расчета характеристик потока воды и прогноза возникновения паводка на участках русла реки Кубань.

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

Практическая значимость. Разработанные методы и модели могут быть использованы для прогноза возникновения паводковых ситуаций в районах, прилежащих к руслу реки, позволяя предотвратить затопление прибрежных построек и сельхозугодий, обезопасить проведение работ по обустройству гидротехнических объектов, находящихся в русле реки. Модели прогноза могут быть применены для расчета экономического ущерба, наносимого паводками; при планировании проведения административно-хозяйственных мер, направленных на предотвращение паводка путем чистки русел, спрямления извилистых потоков, отсыпки дамб и т.д., а также на ликвидацию его последствий.

Результаты диссертационного исследования используются ООО «Вперед» (г. Армавир), а также в учебном процессе ФГБОУ ВПО «Армавирская государственная педагогическая академия», что подтверждено соответствующими актами о внедрении.

Основные положения, выносимые на защиту:

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

Данная методика реализована в программном продукте «COFC», зарегистрированном в Федеральной службе по интеллектуальной собственности, патентам и товарным знакам, и является основой для моделирования прогноза возникновения паводковых ситуаций.

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

Построенные аналитическими методами решения данной упрощенной системы позволяют представить в явном виде зависимости характеристик потока от его длины и времени.

3. Стохастическая модель прогноза возникновения паводковой ситуации в окрестности заданного створа русла горно-равнинной реки.

Построенная модель позволяет увеличить, по сравнению с детерминированными моделями, достоверность прогноза, учитывая наличие помех в результатах измерений характеристик потока.

Апробация работы. Основные результаты диссертационной работы докладывались и обсуждались на следующих научных конференциях: VII объединенной научной конференция студентов и аспирантов факультета компьютерных технологий и прикладной математики КубГУ «Прикладная математика XXI» (г. Краснодар, 2007 г.); V Всероссийской научной конференции молодых ученых и студентов «Современное состояние и приоритеты развития фундаментальных наук в регионах» (г. Анапа, 2008 г.); VI Всероссийской научной конференции молодых ученых и студентов «Современное состояние и приоритеты развития фундаментальных наук в регионах» (г. Анапа, 2009 г); X Всероссийском симпозиуме по прикладной и промышленной математике (Санкт-Петербург, 2009 г. (весенняя сессия); Сочи-Дагомыс, 2009 г. (осенняя сессия)); XI Всероссийском симпозиуме по прикладной и промышленной математике, региональном макросимпозиуме «Насущные задачи прикладной математики на Кубани» (Сочи-Дагомыс, 2010 г. (осенняя сессия)); V Международной конференции «Экологические системы, приборы и чистые технологии» (Москва, 2011 г.).

Область исследования. Содержание диссертационного исследования соответствует пунктам 1, 4 и 7 паспорта специальности 05.13.18 - Математическое моделирование, численные методы и комплексы программ (физико-математические науки).

Публикации. По результатам диссертационных исследований опубликованы 13 печатных работ, в том числе: 7 работ опубликовано в изданиях из перечня ведущих рецензируемых научных журналов и изданий, в которых должны быть опубликованы основные научные результаты диссертаций на соискание ученой степени доктора и кандидата наук. Программный продукт «COFC», используемый для расчета и прогноза глубины и расхода воды на участке русла горно-равнинной реки, зарегистрирован в Федеральной службе по интеллектуальной собственности, патентам и товарным знакам, что подтверждено свидетельством о регистрации программы для ЭВМ.

Структура и объем работы. Диссертационная работа состоит из введения, четырех глав, заключения, списка используемой литературы, содержащего 104 наименования и 3 приложений. Работа изложена на 112 страницах машинописного текста (не включая приложений), содержит 30 рисунков.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе «Обзор математических методов и моделей, используемых для описания речного стока» приведен обзор математических моделей установившегося и неустановившегося движения воды в руслах рек. Рассмотрены основные численные методы решения системы дифференциальных уравнений, описывающей неустановившееся движение потока жидкости. Изложены методы линейной оптимальной фильтрации случайных процессов, используемые в последующих главах для решения поставленных в диссертационной работе задач.

Во второй главе «Анализ системы уравнений Сен-Венана, описывающей неустановившийся поток жидкости» проведен анализ указанной системы аналитическими и численными методами.

Система уравнений Сен-Венана имеет вид:

, (1)

, (2)

при заданных начальных

, , (3)

и граничных

, , (4)

, , (5)

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

При исследовании неустановившегося движения жидкости в достаточно широких водотоках с большим уклоном дна русла инерционными членами , в динамическом уравнении (2) можно пренебречь. Для участка русла горно-равнинной реки с достаточно неизменной по длине потока шириной (), система уравнений сводится к дифференциальному уравнению в частных производных второго порядка:

,(6)

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

Установлено, что при выполнении условий , или , , уравнение преобразуется к уравнению с постоянными коэффициентами (,):

. (7)

Указанные условия выполнимы при малых глубинах: м. Аналитическое решение данного уравнения (при указанных ограничениях) построено методом Фурье:

,

,

, ,

, .

В п. 2.2 построено решение начально-граничной задачи, представляющей собой объединение уравнений конвекции-диффузии (6) и уравнения неразрывности (1) при условиях (3) - (5) методом расщепления ее по физическим процессам.

Разобьем интервал на частичные интервалы , точками , , . Учитывая, что

модель паводковый горный река

,

где - коэффициент шероховатости, построено решение на интервале задачи

, (8)

, (9)

, при , ;(10)

, при , ;(11)

, при ; (12)

, при ; (13)

описывающей процесс конвекции, и задачи

, (14)

, (15)

, при , ;(16)

, при , ;(17)

, при ; (18)

, при , (19)

описывающей процесс диффузии.

Учитывая аддитивность процессов конвекции и диффузии на достаточно малом интервале , приближенное решение общей задачи (1) - (5) имеет вид: , .

Решение каждой из рассматриваемых задач (8) - (13) и (14) - (19) построено конечно-разностным методом. В области введена сетка , , .

Производные в (8) - (13) аппроксимированы разностными отношениями:

; ; .

Будем учитывать изменение ширины потока по длине реки, то есть .

Полученная из (8) - (13) система конечно-разностных уравнений записана в матричном виде:

,(20)

, , , (21)

, , , .

Решение (20), (21) на каждом из интервалов найдено методом матричной прогонки:

, , (22)

где - матрица размерности , - двухкомпонентный вектор.

Прогоночные коэффициенты , вычисляются по формулам:

, , (23)

.

Прямая прогонка выполняется по рекуррентным формулам (23) с начальными условиями , . Затем, с помощью краевого условия , выполняется обратная прогонка и находится , , по формуле (22).

Известно, что метод матричной прогонки устойчив по отношению к случайной ошибке, если выполнено условие ,

которое для рассматриваемой задачи имеет вид: .

Аналогичным образом построено решение задачи (14) - (19). Указано условие устойчивости метода матричной прогонки для этой задачи.

Алгоритм численного решения задачи (1) - (5) методом расщепления реализован в программном продукте «COFC», написанном в среде Maple [14].

Уравнение, полученное из системы уравнений Сен-Венана и описывающее неустановившееся движение жидкости в гидрометрическом створе реки, рассмотрено в п. 2.3:

. (24)

Согласно численным экспериментам, второе слагаемое правой части уравнения (24) пренебрежимо мало по сравнению с остальными слагаемыми и его можно отбросить. Тогда модель движения воды в створе реки предлагается рассматривать в виде:

. (25)

С помощью замены уравнение Риккати (25) сводится к линейному однородному уравнению второго порядка

. (26)

Если , то общее решение уравнения (26) выражается через модифицированные функции Бесселя:

,

где , .

В п. 2.4. изучается задача построения оптимальной в среднеквадратическом смысле оценки расхода воды в створе русла реки.

На первом этапе предполагается, что расход воды описывается с помощью детерминированной модели

, (27)

, (28)

где - ширина рассматриваемого сечения русла, , а результаты наблюдений в каждый момент времени содержат случайные ошибки :

, (29)

где - случайный процесс типа белого гауссовского шума с нулевым средним и корреляционной функцией:

,

- интенсивность белого шума, при . Предполагаем, что функция непрерывно дифференцируема при .

Оптимальная в среднеквадратическом смысле оценка решения , построенная по результатам наблюдений (29) имеет вид:

,

, (30)

,

, (31)

,

,

,

где - приближенное решение (27) при заданном начальном условии .

Затем предполагается, что расход , , является случайным процессом, удовлетворяющим стохастическому уравнению

(32)

с начальным условием (28). При этом - случайный процесс типа белого гауссова шума с нулевым средним, то есть , корреляционная функция которого имеет вид:

,

,

- интенсивность белого шума , - дельта-функция Дирака, и непрерывна на .

Если

,

,

- приближенное решение (27) при заданном начальном условии , то оптимальная в среднеквадратическом смысле оценка процесса , построенная по результатам наблюдений

на интервале является решением фильтра Калмана-Бьюси:

,(33)

,

где определяется из задачи Коши:

, (34)

.

В п. 2.5 предложен аналитический метод расчета средней скорости потока воды в русле реки при неравномерном установившемся движении.

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

Математическая модель прогноза возникновения паводковой ситуации на участке русла горно-равнинной реки при больших уклонах дна русла представлена в п. 3.1 и имеет вид:

,

,

, ,

, ,

, ,

, , ,

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

Определить момент времени начала затопления территории в районе рассматриваемого створа русла реки можно из алгебраических уравнений

или .

В п. 3.2 описана модель прогноза возникновения паводковой ситуации в створе реки, на основе модели неустановившегося движения потока в гидрометрическом створе. Рассматриваемая модель имеет вид:

,

,

,

где , , - критическое для рассматриваемого створа реки значение скорости потока.

Математическая модель прогноза возникновения паводковой ситуации в нижнем створе реки по результатам измерений характеристик потока в верхнем створе приведена в п. 3.3:

,(35)

,(36)

.

В системе (35) - (36) функциональные зависимости и могут быть найдены по результатам аппроксимации наблюдений скорости и глубины потока воды в верхнем по течению -м створе. Из системы линейных уравнений (35) - (36) можно найти функции и в нижнем -м створе и определить их прогнозные значения в момент времени .

Определение момента времени затопления при известном значении при заданном значении сводится к решению уравнения .

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

В п. 3.4 предложена стохастическая модель прогноза возникновения паводковой ситуации на заданном участке русла горно-равнинной реки, использующая методику построения оптимальной в среднеквадратическом смысле оценки расхода воды из задач (27), (28) и (28), (32):

, (37)

, (38)

,

,

где - определяется из (30) - (31) или (33) - (34).

Решение задачи (37) - (38) имеет вид:

,

.

В п. 3.5 описана методика прогноза возникновения паводковой ситуации на участке русла горно-равнинной реки, основанная на использовании уравнения водного баланса. Русло реки предлагается разбить на расчетные участки () таким образом, чтобы их граничными сечениями являлись гидрологические посты, расположенные на реке и ее притоках.

Величина входного гидрографа, то есть функция на некотором участке русла реки, может быть определена, например, с помощью стохастической модели (37), (38). Для прогноза возникновения паводковой ситуации в интересующем створе русла реки необходимо найти функцию , представляющую собой зависимость расхода потока в нижнем створе участка . Для вычисления значений предлагается использовать балансовое соотношение:

,(39)

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

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

,

где - расстояние между створами, - средняя скорость движения потока на заданном участке.

Аналогичным образом могут быть получены величины , , . Для вычисления средней скорости используется соотношение Шези, которое с учетом эмпирической формулы Маннига для коэффициента может быть представлено в виде:

.

В п. 3.6 описаны методики нивелирования профиля русла реки и определения его морфометрических характеристик. Нивелирование продольного профиля осуществляется с помощью геоинформационных систем, позволяющих провести дистанционное зондирование рельефа местности вдоль русла реки и построить эмпирическую зависимость высоты дна русла реки над уровнем моря от продольной координаты . Использование прикладных математических пакетов, программ статистической обработки экспериментальных данных и графопостроительных программ позволяет построить зависимость уклона дна от длины реки , широко используемую в предложенных моделях прогноза возникновения паводковых ситуаций.

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

В четвертой главе «Численные эксперименты прогноза возникновения паводковой ситуации в русле горно-равнинной реки (на примере реки Кубань)» методы и модели, предложенные в главах 2, 3 диссертационной работы, применены для расчета параметров неустановившегося потока воды в русле реки и прогноза возникновения паводковых ситуаций. Выполнено нивелирование продольного профиля русла реки Кубань от ее истоков у подножья г. Эльбрус до впадения в Краснодарское водохранилище с помощью геоинформационных систем.

В п. 4.1 описаны особенности гидрографической сети и водного режима бассейна реки Кубань, отмечены наиболее крупные наводнения за период с 1960 г. по 2010 г.

Расчет расхода потока в русле реки Кубань, выполненный с помощью аналитического метода решения системы уравнений Сен-Венана, описан в п. 4.2. Был рассмотрен участок русла, расположенный в предгорной зоне реки между гидропостами в хуторе Дегтяревском и в селе Успенском. Средняя величина уклона дна на данном участке изменяется от 2,7·10^-3 до 1,8·10^-3, что позволяет использовать для описания движения паводковой волны уравнение (7). Результаты расчетов представлены на рисунке 1 (период времени равен пяти дням).

Рис. 1. Зависимость расхода воды от времени и координаты на участке русла реки от х. Дегтяревский до с. Успенское

Полученное решение позволяет найти величину расхода потока воды в произвольный момент времени в любом интересующем створе рассматриваемого участка русла реки. Согласно проведенным расчетам, критические значения расхода воды на рассматриваемом участке (в х. Дегтяревский (при км) м³/c, в с. Успенское (при км) м³/c) не достигались. Действительно, в рассматриваемый пятидневный период времени паводковой ситуации на данном участке не наблюдалось.

В п. 4.3 приводятся результаты расчетов параметров неустановившегося потока воды на участке русла реки Кубань между г. Армавиром и ст. Темижбекской, произведенные с помощью алгоритма численного решения задачи (1), (3) - (6), описанного в п. 2.2, с использованием программного продукта «COFC» [14]. На рисунках 2 и 3 представлены полученные зависимости расхода и глубины потока от времени в различных точках рассматриваемого участка. Выполнен расчет расхода и глубины потока во внутренних точках участка при заданных граничных условиях (слева от вертикальной прерывистой линии), а также прогноз значений характеристик потока на границах и внутри рассматриваемого участка (на рисунке - справа от прерывистой линии).

Рис. 2. Гидрографы потока (на исследуемом участке реки)	Рис. 3. Зависимость глубины потока от времени (на исследуемом участке реки)

В п. 4.4 описан сравнительный анализ моделей движения воды в створе реки, проведенный на основе результатов численных экспериментов. Подтверждена возможность представления уравнения движения потока в створе русла реки (24) в виде (25) (на основе сравнения результатов аналитических и численных решений уравнений (24) и (25)).

Прогноз изменения основных характеристик потока воды (глубины и средней скорости) на участке русла реки Кубань в районе г. Армавира на основе численного решения системы уравнений Сен-Венана приведен в п. 4.5. В качестве исходных данных взяты результаты измерения скорости и глубины потока во время паводка, наблюдавшегося в июне 2002 г. Графическая интерпретация полученных результатов представлена на рис. 4 и 5.

Рис. 4. Зависимость глубины воды от времени в -м створе	Рис. 5. Зависимость скорости потока от времени в -м створе

Результаты расчетов показывают превышение критического значения глубины потока, которое действительно наблюдалось июня. Определен момент времени начала затопления прилегающих территорий (ч мин). Найдена пропускная способность русла в данном сечении (м³/c), средняя по сечению скорость движения воды при ( м/c) и площадь поперечного сечения ( м²). Полученные результаты хорошо согласуются с данными наблюдений: м/c, м², м³/c.

В п. 4.6. описаны результаты прогноза возникновения паводковой ситуации с 19 на 22 июня 2010 г на участке русла реки Кубань в районе города Армавира с помощью динамико-стохастических моделей (рисунки 6, 7).

Рис. 6. Результаты расчетов, полученные с помощью модели (27), (28)	Рис. 7. Результаты расчетов, полученные с помощью стохастической модели (28), (32)

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

В п. 4.7 приведены результаты расчета расхода воды на участках реки Кубань с помощью уравнения водного баланса, выполненные для бесприточного участка г. Армавир - ст. Ладожская и участка русла реки с боковым притоком с. Успенское - г. Армавир. На рисунке 8 представлен сравнительный анализ данных расхода воды на участке с. Успенское - г. Армавир, полученных по результатам расчетов по формуле (39) и данных непосредственных измерений расхода воды на том же участке в рассматриваемый период времени. Средняя относительная погрешность результатов проведенных вычислений по формуле (39) составляет , что позволяет использовать эту формулу для прогноза значений расхода воды на участке русла реки с боковым притоком - р. Уруп.

Рис. 8. Сравнительный анализ результатов расчета с данными измерений

В п. 4.8 приведены результаты исследования продольного профиля русла реки Кубань от ее истоков до Краснодарского водохранилища с помощью метода, изложенного в п. 3.6. Построена зависимость уклона дна русла реки Кубань:

,

где - пространственная координата, ориентированная по направлению движения воды в русле реки, - уклон дна русла.

Численные эксперименты показывают, что результаты расчетов уклона хорошо согласуются с экспериментальными данными.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

1. Система уравнений Сен-Венана при отсутствии инерционных членов приведена к системе, состоящей из уравнения конвекции-диффузии с переменными коэффициентами и уравнения неразрывности, что позволило разработать алгоритм численного решения полученной системы методом расщепления. Решение каждой из расщепленных задач выполнено методом матричной прогонки. Разработанный алгоритм реализован в программном продукте «COFC» [14].

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

Совокупность результатов, указанных в п.п.1 - 2, позволяет построить решение поставленной в диссертации задачи 1 (которая приведена также и в автореферате).

3. Предложена модель прогноза возникновения паводковой ситуации на предгорных и горных (при больших уклонах дна русла) участках реки. Данная модель основана на использовании уравнений Сен-Венана при отсутствии инерционных членов.

4. Предложена математическая модель прогноза возникновения паводковой ситуации в гидрометрическом створе реки на основе аналитического решения уравнения, описывающего неустановившееся движение потока воды в данном створе.

5. Построена оптимальная в среднеквадратическом смысле оценка расхода воды в гидрометрическом створе реки. Рассмотрены случаи детерминированной и стохастической модели расхода воды в створе.

6. Построена стохастическая модель краткосрочного прогноза возникновения паводковой ситуации на заданном участке русла горно-равнинной реки.

7. Предложена методика прогноза возникновения паводковой ситуации на участке русла горно-равнинной реки, основанная на использовании уравнения водного баланса.

Результаты, указанные в п.п. 3 - 7, позволяют построить решение задачи 2.

8. Предложена методика продольного и поперечного нивелирования профиля русла реки, определения его морфометрических характеристик, нахождения критических значений характеристик потока воды в русле, что позволяет решить задачу 3 диссертационной работы.

9. Предложенные модели и методы применены для расчета и прогноза расхода воды, глубины и средней скорости потока на различных участках русла реки Кубань. Результаты, указанные в п. 9, являются решением задачи 4.

СПИСОК РАБОТ, ОПУБЛИКОВАННЫХ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Вандина Н. В. Расчет скорости потока в открытом русле при неустановившемся движении / Н. В. Вандина // Обозрение прикладной и промышленной математики. - 2009. - Т. 16. - №5. - С. 820-821.

2. Семенчин Е. А. Метод расчета параметров потока на основе решения системы дифференциальных уравнений, описывающей нестационарное движение воды в русле реки / Е. А. Семенчин, Н. В. Вандина // Экологические системы и приборы. - 2009. - №4. - С. 16-20.

3. Семенчин Е. А. Расход воды в сечении русла горно-равнинной реки / Е. А. Семенчин, Н. В. Вандина // Обозрение прикладной и промышленной математики. - 2010. - Т. 17. - №1. - С. 139-140.

4. Семенчин Е. А. Анализ системы уравнений Сен-Венана аналитическими и численными методами / Е. А. Семенчин, Н. В. Вандина // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, - 2010. - №64(10). - С. 1-14. - Шифр Информрегистра: 0421000012/0258. - Режим доступа: http://ej.kubagro.ru/2010/10/pdf/05.pdf.

5. Семенчин Е. А. Построение решения системы уравнений Сен-Венана методом расщепления / Е. А. Семенчин, Н. В. Вандина // Обозрение прикладной и промышленной математики. - 2010. - Т. 17. - №6. - С. 935-937.

6. Семенчин Е. А. Стохастическая модель прогноза возникновения паводковой ситуации в гидрометрическом створе русла реки / Е. А. Семенчин, Н. В. Вандина // Экологические системы и приборы. - 2011. - №8. - С. 42-46.

7. Вандина Н. В. Методика прогноза возникновения паводковой ситуации на участке русла горно-равнинной реки, основанная на использовании уравнения водного баланса / Н. В. Вандина // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, - 2011. - №74(10). - С. 1-13. - Шифр Информрегистра: 0421100012/0409. - Режим доступа: http://ej.kubagro.ru/2011/10/pdf/47.pdf.

8. Ревенко (Вандина) Н. В. Использование геоинформационных систем и программ математической статистики для расчета средней скорости течения и расхода воды в русле Кубани / Н. В. Ревенко // Вестник Армавирского государственного педагогического университета. - 2007. - №1. - С. 78-83.

9. Ревенко (Вандина) Н. В. Методика расчета средней скорости движения воды в русле реки Кубань с использованием геоинформационных систем и статистических программ / Н. В. Ревенко // Прикладная математика XXI в.: Материалы VII объед. науч. конф. студентов и асп. фак. компьютер. технологий и прикл. математики. Часть II. - Краснодар, КубГУ, 2007. - С. 15-17.

10. Ревенко (Вандина) Н. В. Об одном методе решения системы дифференциальных уравнений, описывающей нестационарное движение воды в русле реки / Н. В. Ревенко // Современное состояние и приоритеты развития фундаментальных наук в регионах: Тр. V Всерос. науч. конф. молодых ученых и студентов. Том 2. - Краснодар, 2008. - С. 136-138.

11. Ревенко (Вандина) Н. В. Метод решения системы дифференциальных уравнений Сен-Венана / Н. В. Ревенко // Вестник Армавирского государственного педагогического университета. - 2009. - №3. - С. 83-87.

12. Вандина Н. В. Математическая модель движения жидкости в гидрометрическом створе реки / Н. В. Вандина // Современное состояние и приоритеты развития фундаментальных наук в регионах: Тр. VI Всерос. науч. конф. молодых ученых и студентов. - Краснодар, 2009. - С. 203-205.

13. Вандина Н. В. Математические модели и методы расчета неустановившегося движения воды в русле реки / Н. В. Вандина // Вестник Армавирского государственного педагогического университета. - 2010. - №4. - С. 101-106.

14. Свидетельство о государственной регистрации программы для ЭВМ. Calculation of flow characteristics (COFC) / Н. В. Вандина, Е. А. Семенчин; заявитель и правообладатель Ставропольский государственный университет. - № 2011610596; заявл. 10.11.10; опубл. 11.01.11.- 1 с.

Размещено на Allbest.ru

...