Бесплатный автореферат и диссертация по наукам о земле на тему

Разработка методов определения структурно-гидрографических характеристик по данным ЦМР для гидрологического моделирования

ВАК РФ 25.00.27, Гидрология суши, водные ресурсы, гидрохимия

Автореферат диссертации по теме "Разработка методов определения структурно-гидрографических характеристик по данным ЦМР для гидрологического моделирования"

На правах рукописи

Бугаец Андрей Николаевич

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

Специальность 25.00.27 «Гидрология суши, водные ресурсы, гидрохимия»

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

- 6 ОНТ 2011

Санкт-Петербург 2011

4855609

Работа выполнена в Тихоокеанском институте географии ДВО РАН

Научный руководитель: доктор географических наук

Гарцман Борис Ильич (ТИГ ДВО РАН)

Официальные оппоненты: доктор технических наук

Болгов Михаил Васильевич (ИВП РАН)

кандидат технических наук

Викторова Наталья Владимировна (РГГМУ)

Ведущая организация: Институт географии СО РАН (г. Иркутск)

Защита состоится «ч-7^» Р^^^у^ 2011 г. в «/^

_» часов на

заседании диссертационного совета Д 212.197.02 при Российском государственном гидрометеорологическом университете по адресу: 195196, Санкт-Петербург, Малоохтинский пр., д.98

С диссертацией можно ознакомиться в библиотеке Российского государственного гидрометеорологического университета.

Автореферат разослан « ^ » 20И г.

Ученый секретарь диссертационного совета кандидат географических наук

Воробьев В.Н.

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

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

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

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

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

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

2. Разработка и реализация базы данных структурных, гидрографических и геоморфологических характеристик.

3. Апробация методики анализа ЦМР на примере крупных речных систем в бассейне Амура.

4. Усовершенствование модели геоморфологического мгновенного единичного гидрографа (КЛУ-ОШН) за счет учета параметров, полученных по индивидуальным структурным элементам речного

бассейна.

5. Адаптация методики определения эффективных осадков при расчете гидрографов стока с помощью модели KW-GUIH.

6. Апробация усовершенствованной модели с использованием данных бывшей Приморской водобалансовой станции.

Методической основой структурно-гидрографического и морфометрического анализа являются алгоритмы, реализованные на базе ГИС технологий (ERDAS IMAGINE 9.1, ESRI ARCMAP 9.3) и оригинальных программных компонент разработанных в среде DELPHI 7. В качестве базовой СУБД применялась MS Access. Методической основой, используемой при расчете гидрографов дождевого стока, является модель KW-GUIH, разработанная К.Т. Ли с соавторами, прошедшая широкую апробацию за рубежом. Испытание модели в России впервые выполнено автором диссертации. В работе использованы известные методы математического моделирования, а также общепринятые методы статистического и географического анализа.

Объектами исследования являются разнопорядковые речные системы, расположенные в пределах бассейна Амура и прилегающем бассейне северозападного побережья Японского моря. Среднемасштабная ЦМР на всю исследуемую территорию создана на основе данных SRTM (Shuttle Radar Topographic Mission). Для тестовых испытаний различных версий модели KW-GIUH были использованы опубликованные данные бывшей Приморской воднобалансовой станции (ПВБС).

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

Усовершенствована модель KW-GUIH путем учета пространственной неоднородности, как параметров водосборов, так и осадков. В рамках решения частных задач впервые:

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

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

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

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

Архитектура географической базы данных морфометрических, гидравлических и структурно-гидрографических характеристик речных бассейнов, получаемых на основе объективного анализа цифровых моделей рельефа, с применением стандартных процедур, предлагаемых современными ГИС, и оригинального программного комплекса. Результаты опытного структурно-гидрографического анализа речных систем бассейна Амура и принципиальная схема их типизации на основе распределения суммарного индекса напряжения сдвига (shear stress) от расстояния до замыкающего створа

Усовершенствованная модель геоморфологического мгновенного единичного гидрографа (KW-GUIH).

Результаты тестирования различных версий модели KW-GIUH с применением данных приборов с непрерывной записью ПВБС, пространственно-распределенных данных об осадках и гидравлических характеристик подстилающей поверхности.

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

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

Апробация работы. Основные положения и выводы диссертации докладывались автором на следующих конференциях: Всероссийском гидрологическом съезде (СПб, 2004), региональной конференции «Научные основы экологического мониторинга водохранилищ» (Хабаровск, 2005),

региональных молодежных конференциях «Географические и геоэкологические исследования на Дальнем Востоке» (Владивосток, 2006, 2008), XIII совещании географов Сибири и Дальнего Востока (Иркутск,

2007), Генеральной ассамблее европейского геофизического союза (Вена,

2008), региональной конференции «Современные геофизические и географические исследования на Дальнем Востоке России» (Владивосток, 2010), всероссийской конференции с международным участием «Фундаментальные проблемы воды и водных ресурсов» (Барнаул, 2010), международной конференции Pan Ocean Remote Sensing Conference (Килунг, 2010). XIV Совещание географов Сибири и Дальнего Востока (Владивосток, 2011).

Результаты работы применялись при выполнении тематики плана НИОКР 1.1.2 Росгидромета («Разработать, испытать и внедрить новые технологии и методы гидрологических прогнозов для рек и водохранилищ России, включая прогнозы наводнений и других опасных гидрологических явлений. Обеспечить выполнение международных обязательств Росгидромета в области гидрологических прогнозов», ДВНИГМИ, 2007). Также они использованы в рамках проекта «Модернизация и техническое перевооружение учреждений и организаций Росгидромета. Совершенствование системы гидрологических наблюдений, в т.ч. в паводкоопасных районах. Ч.З. Модернизация гидрологической сети в бассейне р.Уссури», при подготовке технического задания на разработку региональной гидрологической модели (РГМ), отдельные блоки которой проходят тестовые испытания в отделе гидрологических прогнозов Приморского УГМС.

Основные положения диссертации изложены в 10 публикациях, в том числе в 3 статьях, опубликованных в журналах списка ВАК.

Диссертация состоит из введения, 5 глав, заключения, и приложений, объемом 215 страниц, включая 55 рисунков, 18 таблиц, списка используемой литературы, включающего 220 наименований и 4-х приложений.

ОСНОВНОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

Первая глава диссертационной работы посвящена изложению современного состояния детерминированного моделирования речного стока и анализа структуры речных систем. Дан обзор работ М.А. Великанова, Г.А. Алексеева, А.Н. Бефани, Г.П. Калинина, В.Д. Комарова, Е.Г. Попова, Ю.Б. Виноградова, JI.K. Шермана, Р.К. Линдслея, М.А. Колера, Д.И. Нэша, М. Киркби и др., которые легли в основу современных представлений о стокообразующих процессах и способствовали внедрению моделирования в практику. Описана иерархия и обзор существующих математических моделей данного типа, применяемых на практике.

Во второй части представлен краткий обзор наиболее значимых публикаций и результатов в области исследования структурно-функциональных отношений в речных системах. Выполнен обзор работ P.E. Хортона, А.Н. Штралера, А.Е. Шайдеггера, P.E. Шрива, В.П. Философова, JIM. Корытного, H.A. Ржаницына, Е.А. Черных, И.Н. Гарцмана, Б.А. Казанского, П. Реджиани, В. Гупты, И. Родригеса-Итурбе и А. Ринальдо и др.

Во второй главе рассмотрены методические аспекты и программные средства структурно-гидрографического анализа. Дано описание методики анализа ЦМР с применением стандартных функций, входящих в общедоступные программные пакеты. Приведена концепция «элементарных структур» положенная в основу базы данных и алгоритмов расчета бассейновых характеристик и моделирования стока.

Впервые метод построения речной сети и определения структурных и геоморфологических параметров по ЦМР был предложен в работе (O'Callaghan & Mark, 1984). С тех пор он широко используется на практике и реализован в большинстве ГИС в виде стандартных функций. Ключевым моментом этого метода является определение пороговой (минимальной) площади формирования водотока. Данный вопрос рассматривается как крайне важный и связан с решением известной проблемы «где начинается река?», ему целиком посвящен второй раздел главы в котором исследуется зависимость результатов оцифровки от разрешения ЦМР и ландшафтных характеристик бассейна, в особенности его рельефа.

Анализ выполнен по территории бывшей ПВБС и включает бассейны р. Комаровки и её правого притока - р. Раковки (1500 км2). В качестве ключевых участков были выделены малые водосборы приблизительно одинакового размера, которые в своих характеристиках достаточно полно отражают разнообразие природных условий водосбора р. Комаровка (рис. 1).

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

Для основных рабочих масштабов карт (1:25000-1:100000) диапазон оптимальной пороговой площади составляет от 0.2-0.5 км2 для первой группы водосборов, 0.5-1.0 км2 для второй, 1.0 км2 для третьей и 1.0-2.0 км2 для четвертой. Наиболее устойчивые результаты сравнения получены для водосборов предгорий (группы 2 и 3), где количество нерусловых элементов гидрографической сети минимально. В среднем величина пороговой площади может быть оценена 0.5 - 1 км2, что отвечает средней площади водосбора 1-го порядка 1 - 2 км2.

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

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

№ п/п Название водотока Я ч ев м 3 2 о * е С Ср. высота, м Ср. уклон водосбора, %0 Группа

1 Комаровская Падь 24.7 338 19.4 1

2 Учхозный Ключ 31.0 366 18.9 1

3 Комаровка 45.7 341 17.9 1

4 Каменушка 31.9 270 17.7 1

5 Раковка 25.2 396 14.7 2

6 Падь Таскаева 25 314 14.7 2

7 Каменка 30.2 351 14.5 2

8 Малая Раковка 30.9 373 14.4 2

9 Ивнячка 28.6 246 13.4 2

10 Лихачевка 40.9 276 11.7 2

11 Барсуковка 36.5 309 10.6 2

12 Глуховка 29.9 91 6.4 3

13 Васильевка 31.1 91 6.0 3

14 Бакарасьевка__ 42.6 94 4.4

15 Падь Дорожная 24.0 90 2.5 4

16 Михайловка 42.5 82 2.0 4

17 Репьевка 32.2 82 1.9 4

а Горные (до 700 м) □ Горные (300 - 500 м) Г~1 Предгорные (до 300м) ЕЗ Равнинные

Рисунок 1. Характеристики экспериментальных водосборов и карта-схема расположения групп в пределах бассейна р. Комаровка.

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

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

виде древовидной сети «элементарных структур» - водотоков, каждый из которых рассматривается вместе со своим частным водосбором.

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

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

Поле осадков Гидрографическая сеть Частные водосборы ЦМР

Рисунок 2. Пространственное объединение данных.

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

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

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

3) Звенья гидрографической сети объединяются в древовидную структуру топологическими отношениями, где каждый элементарный водосбор имеет «указатель» на подобную структуру выше и ниже по течению.

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

Представленная таким образом информация органично размещена в базе данных с точной привязкой к географической основе. БД представляет собой набор сводных таблиц, описывающих различные параметры строения исследуемых речных систем (табл. 1) Записи, содержащиеся в различных таблицах, связаны друг с другом и пространственными данными (тематическими слоями карт) полем «ГО».

Таблица 1. БД структурных, гидрографических и геоморфологических характеристик: имена таблиц и их краткое описание.

Имя таблицы Описание таблицы

ORDERS порядок в системе Штралера, Шрива и Шайдеггера

SumStrOrd количество рек в системе Штралера

Distance дистанция от нижнего узла до истока и до замыкающего створа

Entropy энтропия по формуле Шеннона

StrOrdEntropy энтропия по Шеннону в точках перехода между порядками Штралера

HrtStrGmph общие геоморфологические показатели, вычисленные для замыкающий створов рек в соответствии с системой Штралера

MainRiver Главная река (секции, составляющие главное русло)

SumRvrLng суммарная длина рек

SumWsdArea суммарная площадь водосбора

MainRvrDstArea изменение площади водосбора вдоль главной реки

WidthFunction график ширины бассейна

RvrAvrWgtSlp средний взвешенный уклон рек

WsdAvrWgtHgt средняя взвешенная высота водосбора

WsdAvrWgtSlp средний взвешенный уклон водосбора

TMtxStrOrd матрица перехода рек между порядками

TMtxRvrLng матрица перехода длин рек между порядками

TMtxRvrSlp матрица перехода уклонов рек между порядками

TMtxWsdArea матрица перехода площадей водосборов между порядками

TMtxWsdSlp матрица перехода уклонов водосборов между порядками

В качестве базовой СУБД была выбрана Microsoft Access. Ядро базы данных Microsoft Jet имеет встроенную поддержку интерфейсов OLE DB и расширенный синтаксис языка SQL. Эту систему можно рассматривать и как среду разработки для создания ориентированных на пользователя приложений таких же мощных, как и приложения, написанные на языках программирования. Полная интеграция с другими приложениями MS Office и

встроенный язык Visual Basic for Application позволяет выполнять нестандартные задачи. В итоге получен универсальный продукт, расширение которого может быть ориентировано на практически любую конкретную задачу. Благодаря соответствию внутреннему формату ESRI объекты базы данных могут органично управляться внутри данной ГИС.

Третья глава целиком посвящена практическому применению разработанной методики анализа ЦМР на примере крупных речных систем (от 20 до 70 тыс. км2, субрегиональный масштаб) в пределах бассейна Среднего и Нижнего Амура, а также бассейна Японского моря (табл. 2). Для их выделения назначались, по возможности, «естественные» замыкающие створы - места впадения в более крупную реку или слияния с крупным (равноценным) притоком.

Таблица 2. Основные сведения о водосборах.

Площадь Ср. высота Ср. уклон Длина Уклон

Бассейн водосбора, водосбора, водосбора, главной главной

тыс. км2 м %о реки, км реки, %о

Амгунь 54.8 524 17.5 812 1.08

Бурея 70.5 692 14.4 783 2.46

Горин 22.2 315 10.7 336 1.06

Гилюй 22.4 740 13.7 514 4.08

Зея 59.1 667 11.9 489 2.47

РФ Селемджа 68.7 554 11.6 652 2.46

Тумнин 21.8 610 19.5 352 1.12

Тунгуска 30.1 475 15.5 481 1.08

Хор 24.4 656 21.9 423 3.35

Бикин 21.4 525 15.1 493 3.59

Уссури 27.0 485 17.2 421 1.92

Б. Уссурка 29.5 562 19.8 445 3.78

Сунгач 24.7 193 5.38 305 1.04

Сургари 52.7 577 11.8 692 2.78

Нэньцзян 66.4 507 8.03 658 1.12

Хуланьхэ 36.4 267 3.49 401 0.602

Хулинхэ 38.2 356 8.38 613 5.79

i 1 Лалинхэ 19.6 325 6.87 351 1.26

к Муданьцзян 37.3 579 14.7 682 1.6

Ui Нуоминхэ 25.6 637 12.9 424 1.29

Танванхэ 20.7 474 13,1 424 1.87

Таоэрхэ 32.7 583 9.75 544 0.723

Уюр 41.5 223 1.21 678 0.397

Ялухэ 19.2 560 13.0 347 0.924

Вся исследуемая площадь покрыта однородной ЦМР и все расчеты выполнены по единой автоматизированной методике в одном масштабе. Среди выбранных речных систем имеются водосборы, содержащие в себе крупные озера и водохранилища (Сунгач, Зея, Горин, Амгунь), а так же с обширной равнинной частью, отягощенной антропогенной нагрузкой (Хуланьхэ, Хулинхэ, Лалинхэ, Уюр). В равнинных бассейнах речная сеть

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

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

2 3 4 5 6 Порядок реки

3 4 5 6 Порядок реки

Рисунок 3. Графики связи структурных морфометрических характеристик с порядком в системе Хортона-Штралера.

Из анализа данных зависимостей следует, что закон Хортона для длин водотоков нарушается в бассейнах, содержащих большие водные объекты, в то же время остальные законы продолжают работать. Средняя площадь водосборов низких порядков остается практически постоянной для всех исследуемых рек. Следовательно, увеличение длины водотока первого порядка в озерных бассейнах не порождает пропорционального увеличения площади его водосбора, но структура речной сети при этом, безусловно, претерпевает некоторые изменения. Это особенно заметено на графике зависимости суммарной энтропии от общей протяженности гидрографической сети (Сунгач, Зея, Горин на рис. 4). Данный параметр

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

10000 20000 30000 40000 50000 Протяженность гидрографической сети, км

Рисунок 4. График связи суммарной энтропии и общей протяженности гидрографической сети исследуемых бассейнов.

На долю малых водотоков повсюду приходится значительная часть водосборной площади. Суммарная площадь бассейнов младших порядков вплоть до максимального порядка £2-2, остается практически постоянной, занимая около 60% территории, и затем резко возрастает. Подобным свойством обладает модель бассейна Пеано (Реапо) - одна из первых и наиболее простая модель идеализированной речной сети водосбора (рис. 5).

Рисунок 5. Распределение водосборной площади в порядковой системе Штралера для исследуемых бассейнов и бассейна Пеано 8-го порядка.

В разделе 3.2 предложен новый инструмент анализа и классификации речных бассейнов на основе т.н. «энергетической функции» водосбора.

Определение энергии рельефа в данном случае тесно связанно с понятием эродирующей силы т.н. напряжения сдвига (shear stress). Согласно

Хортону, сопротивляемость эрозии данного участка поверхности количественно может быть выражена в единицах силы, требуемой для появления эрозии, т.е. в кг/м2. Работа, расходуемая в единицу времени на преодоление трения на участке склона 1x1 м при установившемся течении, равна энергии объема воды, протекающей через единицу площади в единицу времени, т.е. произведению веса, уклона и скорости. Таким образом, полная диссипация энергии Е флювиальной системы данной конфигурации s может быть рассчитана по формуле (Rodrigues-Iturbe, Rinaldo, 1997):

£(■*) = S^ÖA" (1)

где р - плотность воды, g - ускорение свободного падения, Q. и bzi -соответственно расход и уклон в произвольной точке i. Поскольку плотность воды величина постоянная, а расход воды пропорционален размеру дренируемой площади, то выражение (1) допустимо записать в виде:

E(s)=ZA,&zl, (2)

где А, и Az, - соответственно площадь и уклон водосбора в точке г.

Применив указанное преобразование к функции распределения площади водосбора вдоль гидрографической сети W(x), получим характеристику Eis) в контексте указанной модели, т.е. распределение напряжения сдвига на расстоянии л от замыкающего створа, рассчитанного вдоль русловой сети -«энергетическую функцию» водосбора Е(х).

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

Связь характера Е{х) с локальными геологическими формами речных долин легко устанавливаются в бассейнах большинства исследуемых рек. Интересным примером с указанной точки зрения является водосбор р. Амгунь (рис. 6, в). Это бассейн неправильной формы, здесь происходили события изменения границ и перестройки гидрографической сети связанные с движениями земной коры. В верхнем течении основное русло огибает литосферный разлом, меняя приблизительно в 100 км от истока основное направление течения с ЮЗ на обратное СВ, сохраняя параллельное направление еще на расстоянии не менее 100 км. В средней и нижней части бассейна участки явного несоответствия речной сети рельефу так же говорят о геолого-структурной предопределенности ее характера. Функция Е(х) особым образом меняет исходное распределение площадей w(x),

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

- w{x) - s(x)- Е(х)

Рисунок 6. График распределения площади w(x), уклонов S(x) и напряжения сдвигаЕ{х) рек: а) Уссури; б) Хулинхэ (КНР); в) Амгунь.

240 I | 60

О 100 200 300 400 500 Расстояние от замыкающего створа, км

Расстояние от замыкающего створа, км

Расстояние от замыкающего створа, км

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

Базовый тип А1 обладает весьма интересным свойством, заключающимся в линейной зависимости средневзвешенного уклона водосбора от расстояния до замыкающего створа (рис. 8). Распределение длин секций русловой сети этого типа хорошо описывается логнормальным законом распределения, часто применяющимся для описания природных систем с ассиметричным распределением параметров. Для описания такой системы подойдет модель бассейна Пеано, в которой уклон водосбора будет линейно возрастать в соответствии с удалением от замыкающего створа, а длины секций распределены логнормально (рис. 9). Остальные типы, очевидно, могут быть получены аналогично, т.е. посредством моделей описывающих распределение площадей IV(х) нуклонов 5(х).

100 200 300 400 Расстояние, км

ь»

II 111 1 и 0-02 illSllih..

I

■■■III

0 1 2 3 4 5 6 7 Длина секций гидрографической сети, км С.] Эмпирическое —Логнормальное

Рисунок 8. Река Уссури, распределение а) средневзвешенного уклона

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

0.02 0.04 0.06 0.08 0.1 0.12 Расстояние до замыкающего створа

— Е(х) — S(x) — W(x)

Рисунок 9. Моделирование типа А1 с помощью бассейна Пеано.

Представляется обоснованным вывод, что функция Е(х) является принципиально новым инструментом анализа морфологии речных бассейнов и их типизации. Кроме этого, функция Е[х) может быть полезным инструментом для разработки или улучшения структуры моделей типа «осадки-сток», основанных на разбиении бассейна на гидрологически однородные части.

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

Группа А. Горные водосборы. Амплитуда максимальных значений 50 - 25

Подгруппа А1. Хорошо выработанные горные бассейны, которые рассматриваются как наиболее простой (базовый) тип. Формируется при взаимодействии одного горного сооружения с одной межгорной впадиной, т.е. в условиях простого тектонического перекоса и умеренной скорости вертикальных движений. Характеризуются правильной формой бассейна и наиболее простой формой функции Е, которая имеет один ясно выраженный смещенный вправо пик, образованный от слияния основных притоков. Речная сеть хорошо согласованна внутренне, а также с главными элементами рельефа.

О 100 200 300 400 500

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

О 100 200 300 400 500

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

Подгруппа А4. Функция Е имеет разбросанные пики и не убывает практически на всем интервале. Главный пик смещен далеко вправо. Бассейны имеют небольшие размеры и неправильную, сильно вытянутую форму, не имеют крупных притоков. Эта группа объединяет бассейны, отставшие от соседних в процессе конкурентного развития и не сформировавши^ _

0 100 200 300 400 500

Подгруппа А5. Эту группу в основном составляют самые крупные водосборы, имеющие преимущественно горный характер. Из-за своих размеров эта группа включает в том или ином качестве черты предыдущих. При относительно высоких значениях, функция Е имеет «размазанный» пик, смещенный вправо. Это обусловлено последовательным впадением крупных притоков

0 200 400 600 800

Группа В. Предгорные водосборы. Значения максимумов располагаются в интервале 15-25,

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

50 40 30 ;

20 '

0 100 200 300 400

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

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

0 100 200 300 400

Группа С. Равнинные водосборы. Максимумов значение Е(х)-2.5.

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

0 150 300 450 600 750

л с к

а, и Ч о С

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

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

В структуре модели КЛ¥-ОШН такая возможность заложена, она заключается в переходе от типов путей добегания воды к их индивидуальному учету. В этом случае, осредненные по порядкам

Гилюй Зея

Амгунь

Бурея

9 Нэньцзян

Нуоминхэ

Та'незно

(уланьхэ

[А4| Бикмн

(улинхэ

.Большая ' Уссурка ,

¡Даоерхэ

Папинхэ

1удань ЦЗЯН ;

■Г и^Л^Ч

■ А1

■ А2

■ АЗ

■ А4

□ А5

□ В1

■ В2

■ вз

■ С1

Рисунок 7. Карта-схема типизации по характеристикам функции е{х).

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

Таким образом, модель КЛУ-ОЮН может быть органично преобразована в пространственно-распределённую модель стока. Для определения эффективных осадков вместо используемой в модели КАУ-вШН иифильтрационной теории Хортона, применяется широко опробованная в пределах юга Дальневосточного региона, процедура в виде части алгоритма модели паводочпого цикла малого речного бассейна (Гарцман, Бугаец, 2003).

Перевод модели на новую структурную основу потребовал внести изменения в ее основные формулы. В новой редакции время добегания по склону Тк (с которого начинается путь IV) определяется по формуле:

" ' г (3)

.'ГлЛ

где /г, - коэффициент шероховатости склона; - средний уклон склонов водосбора; т - показатель степени 5/3 из уравнения Маннинга; ¡е -эффективные осадки. Параметр Л„ обозначает среднюю длину склона элементарного водосбора и определяется, согласно принятому в модели V-образному представлению водосбора, по формуле:

1„ = а/2Лс, (4)

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

=АР11Л /2Ы^С, , (5)

где А - площадь водосбора в замыкающем створе; Ы,, Ъи - соответственно количество водотоков и средняя протяженность водотока ¡-го порядка, РаА -вероятность попадания осадков на частную водосборную площадь ¡-го порядка определяемая по формулам:

Кл = Л, = —;.! 1 У.! , (б)

оА' А{ ' ) " Л',. ^^

где А/ - средняя площадь водосбора ¡-го порядка. Р - вероятность впадения водотока ¡-го порядка в водоток порядка /; - количество водотоков ¡-го порядка, впадающих в водотоки у-го порядка.

Время добегания для каждой секции гидрографической сети определяется расходом в верхнем створе (9, гидравлическими характеристиками русла и боковым притоком </„ с элементарного водосбора

примыкающего к этой секции:

8 к

Яо

Т„ =

где ,УС - средний уклон русла; В - ширина водотока; Ьс - длина секции; « -шероховатость русла. Параметр кса обозначает уровень воды, обусловленный расходом <2 в верхнем створе, сформированным частью водосбора дренируемой данной секцией и гидравлическими характеристиками русла:

В 5,

О

(8)

Изменение уровня воды в нижнем створе участка обусловлено боковым притоком воды ц0 с элементарной водосборной площади а:

<?„ = '>• (9)

Исходные формулы, используемые в прототипе для расчета времени добегания по руслам, выглядят следующим образом:

т. =

2Цо

21епсЬо, Ьс, П. .V.

й„„ =

ХМ .V

(10)

где ис - коэффициент шероховатости русла; /.<„ - средняя длина склона /-го порядка. Остальные параметры означают то же что и в предыдущих формулах но осредненные для i -х порядков в системе Штралера.

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

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

Т. = Г.. + А Г„. + • • • + + ГвГг, (11)

где /,..., П - порядки Хортона - Штралера (О - максимальный порядок) образующие путь ТГ1 - суммарное время добегания по секциям гидрографической сети принадлежащим пути со и имеющим й порядок в системе Хортона - Штралера рассчитанные по формуле:

(12)

где Г - время добегания по секции принадлежащей пути со и имеющим / -й порядок в системе Штралера, рассчитанное по формуле (7), у - номер секции, п - количество таких секций.

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

Пятая глава целиком посвящена тестированию модели и анализу результатов. В качестве объектов были выбраны 6 водосборов в бассейне р. Комаровка (табл.4, рис.11). В процессе модельных экспериментов тестировались как исходная, так и усовершенствованная версии модели KW-GIUH. Источником для структурной основы модели послужила ЦРМ 10x10 м, созданная по оцифрованным горизонталям и высотным отметкам топографических карт м-ба 1:50000.

Источником информации об осадках и измеренных расходов воды послужили данные плювиографов и гидрологических постов ПВБС, из которых выбирались отдельные паводковые события на основе анализа ливневых осадков и прохождения тайфунов. В результате выбраны два паводка 1968 и 1974 года, как наиболее полно освещенные наблюдениями за осадками и стоком по всему бассейну р. Комаровки. Паводок 16-21 августа 1968 связан с выходом на центральную часть Японского моря тайфуна «Полли» (Polly). В результате над территорией Приморья прошли сильные дожди, суммы осадков составили до 200 мм и более. Паводок 2-5 октября 1974 был обусловлен необычно активной циклонической деятельностью. Суточные осадки на отдельных станциях превышали 100 мм.

При расчетах в модели применялся часовой шаг по времени. Интерполяция полей осадков производилась методом обратных взвешенных расстояний (IDW). Поскольку алгоритм основной версии модели предусматривает использование только средних по бассейну осадков, для моделирования стока были использованы осадки, рассчитанные как средневзвешенное значение. Усовершенствованная версия испытывалась как со средними по бассейну, так и с пространственно-распределенными осадками.

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

пространственное распределение данного параметра в виде прямой зависимости от уклона элементарного водосбора:

и, (13)

В таблице 5 представлены данные измерений и использованные для калибровки параметры модели (всего 5) включая: параметры шероховатости русел и склонов, объемы измеренного и расчетного паводка V мм, величины измеренного ()тах и расчетного ()рас, максимума м3/с, а так же положение расчетного максимума относительно измеренного. Отрицательные и положительные значения сдвига соответствуют более раннему или более позднему прохождению расчетного максимума относительно наблюденного.

. Дубки

Сун Ят-сеи

. Дальстрой г Дубининский Михайловский „ Дорожный

Стоковый *. Опытный Пролетарский. Раковский

Сахарный Завод .п 0'В0Й Уссуриискии • # Мостовой Баневуровский• 0 Глуховский

» Дубовый Лечебный • « Кривой

Садовый • Медвежий Лесничий » Центральны« -Доковский " Кондратенковский

Сахарныи\-с& , ваковский

Нижний 'Дальний ° Верхний

•Комаровский Егерьскии

Рисунок 11. Расположение водосборов и осадкомеров по территории ПВБС.

Таблица 4. Основные сведения о водосборах.

Река - пункт Площадь водосбора, км2 Расстояние от истока, км Ср. уклон реки, %0 Ср. высота водосбора, м Ср. уклон водосбора, %о

Комаровка-Центральный 157 25 7.4 300 260

Комаровка-С адовый 395 37 5.7 260 206

Комаровка-Сах. Завод 616 61 3.4 210 105

Глуховка-Мостовой 31.1 5.9 7.2 100 62.7

Раковка-Раковский 198 35 8.3 272 137

Раковка-Опытный 755 69 3.5 130 60.8

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

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

Таблица 5. Характеристики паводков и параметры моделей.

иин Усовершенствованная К\У-СШН Усовершенствованная К\У-ОШН

п^Пц-сот! | пс,п„=((8) К\у-

паводок 16-21.08.1968 вшн паводок 02-05.10.1974

Ср. | Расп. | Ср. | Расп. Ср. 1 Расп. | Ср. 1 Расп.

р.Комаровка-п.Центральный

У=98.4 мм, В=150 м,0„„=165 м'/с, 17.08 15:00 У=63.8 мм, В=100 м, 0„„=58.7 м3/с, 03.10 20:00

Пс 0.1 0.1 0.1 0.1 Бс 0.1 Бс 0.1 0.1 0.1 0.1 Бс 0.1 Эс

по 12.0 12.0 12.0 0.5 Бо 0.5 Бо 11.0 11.0 11.0 0.5 Эо 0.5 Эо

Ррас, м3/с 163 162 162 170 170 60.5 59.0 58.7 60.4 60.1

Сдвиг, час -1 -1 -1 -1 -1 0 +3 +3 +3 +3

V, мм 93.7 93.2 94.0 93.2 94.0 57 56.7 57.1 56.6 57.1

р.Комаровка - п.Садовый

У=86.9 мм, В=400 м, 0„„=297 м'/с, 17.08 15:00 У=49.1 мм, В=300 м, 0„»=120 м3/с, 04.10 00:00

пс 0.1 0.1 0.1 0.5 5с 0.5 Бс 0.05 0.05 0.05 0.05 Бс 0.05 Бс

па 20 20 20 0.5 Бо 0.5 Бо 16 16 16 1.0 1.0

Срас, м'/с 306 291 288 297 291 124.9 120.9 120.2 119.5 118.5

Сдвиг, час +2 +3 +4 +6 +6 + 1 +1 + 1 0 0

У,мм 93.1 92.8 93.4 92.4 92.9 52.2 51.8 51.4 51.6 51.3

р.Комаровка - п.Сахарный завод

У=66 мм, В=600 м, С?„„=412 м /с, 17.08 23:00 У=52.6 мм, В=400 м, 0„„=187 м'/с, 04.10 21:00

пс 0.1 п 0.1 0.1 0.1 Бс 0.1 Бс 0.05 0.05 0.05 0.06 Бс 0.06 Бс

"о 10 10 10 1.0 Бо 1.0 Бо 8.0 8.0 8.0 0.6 Бо 0.6 Эо

Орас, м'/с 402 411 385 402 390 192.9 190.8 183.8 193.3 194.1

Сдвиг, час -2 -2 -1 -5 -6 -19 -19 -19 -22 -22

V, мм 79 79 79 79 79 47.0 47.0 47.3 46.2 47.1

р.Раковка - п.Раковский

У=70.8 мм, В=150 м, 0„„=189 м3/с, 17.08 19:00 У=46.8 мм, В=100 м, 0„„=54.2 м3/с, 04.10 01:00

пс 0.06 0.06 0.06 ^ 0.1 Бс 0.1 Бс 0.06 0.06 0.06 0.1 Бс 0.1 8с

По 6 2.4 2.4 0.4 Бо 0.4 Бо 10.0 6.0 6.0 0.7 Эо 0.7 &о

Орас\ м'/с 192 187 181 190 183 55.6 54.5 54.1 55.3 55.2

Сдвиг, час -4 -2 -4 -4 -4 -1 +2 +2 -1 -2

V, мм 78.7 79.9 79.5 80.0 79.4 42.8 44.1 43.3 43.9 43.1

р.Раковка - п. Опытный

У=20.7 мм, В=600 м, 0„„=128 м'/с, 18.08 19:00 У=8.2 мм, В=400 м, Ош„=39.6 м3/с, 06.10 02.00

пс 0.06 0.06 0.06 0.1 Бс 0.1 Бс 0.06 0.06 0.06 0.1 Бс 0.1 Бс

и„ 7 7 3 1.3 во 1.3 во 16 16 16 2.5 Бо 2.5 Бо

Орас, м '/с 133 130 132 131 131 36.8 39.9 38.5 39.8 41.2

Сдвиг, час -19 -15 -16 -18 -19 -38 -31 -37 -34 -42

V, ми 27.5 27.2 27.7 28.2 27.6 10.9 11.7 12.0 11.7 12.2

р.Глуховка - п.Мостовой

У=46.3 мм, В=50 м, С?„„=47.8 м3/с, 17.08 13:00 У=16.2 мм, В=30 м, 0™=7.40 м3/с, 03.10 21:00

пс 0.06 0.06 0.06 0.01 Бс 0.01 вс 0.05 0.05 0.05 0.1 Эс 0.1 Бс

п0 0.14 0.14 0.14 0.03 Бо 0.03 Бо 5.0 5.0 5.0 0.6 Эо 0.6 Бо

Срас, м'/с 50.6 47.4 47.4 47.5 47.6 6.95 7.24 7.25 7.81 7.78

Сдвиг, час -1 -I -2 -2 -2 -4 -5 -5 -4 -4

У,мм 54.4 55.4 56.5 55.4 56.4 29.9 30.7 31.3 30.7 31.3

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

расчетных максимумов расходов и объемов паводка (таблица 6) рассчитанные по формулам:

0-0 V -V

<7е = - ^ • 100; оу = ишу рас ■ 100 (14)

где Qm и Qrc - соответственно измеренный и рассчитанный максимумы

расхода, а уиы и v - объемы наблюденного и расчетного гидрографов.

Анализ таблицы 6 показывает, что для обеих версий модели значения с и о-,, имеют близкие значения. Средние значения не превышают

6 %. Результаты усовершенствованной версии, по сравнению со своим прототипом, дают меньшие значения ffQm притом что, рассчитывались по большему количеству проведенных экспериментов. Средние значения аг для малых горных и средних водосборов остаются в пределах 10%. Наибольшие отклонения а> отмечены для рек Глуховки и Раковки - п. Опытный.

Таблица 6. Относительные отклонения етд и аг

Река-Пункт KW-GIUH Усов. KW-GIUH KW-G1UH Усов. KW-GIUH Ср. сг

пс,п0—const nc,na=const Пс,По=№ KW-GIUH Усов. KW-GIUH

паводок 16-21.08.1968 паводок 02-05.10.1974

Ср. Расп. Ср. Расп. Ср. Расп. Ср. Расп.

Комаровка -Центральный 1.21 1.82 1.82 3.03 3.03 3.07 0.51 0.00 2.90 2.39 2.14 1.94

(Ту 4.78 5.28 4.47 5.28 4.47 10.7 11.1 10.5 11.3 10.5 7.72 7.87

Комаровка -Садовый 3.03 2.02 3.03 0.00 2.02 4.08 0.75 0.17 0.42 1.25 3.56 1.21

<тг 7.13 6.79 7.48 6.33 6.90 6.31 5.50 4.48 5.09 4.48 6.72 5.91

Комаровка -Сах. Завод 2.43 0.24 6.55 2.43 5.34 3.16 2.03 1.71 3.37 3.80 2.79 3.18

аг 19.7 19.7 19.7 19.7 19.7 10.6 10.6 10.1 12.2 10.5 15.2 15.3

Раковка -Раковский 1.59 1.06 4.23 0.53 3.17 2.58 0.55 0.18 2.03 1.85 2.09 1.70

(Ту 11.2 12.9 12.3 13.0 12.1 8.50 5.80 7.50 6.20 7.90 9.85 9.70

Раковка -Опытный 3.91 1.56 3.13 2.34 2.34 7.07 0.76 2.78 0.51 4.04 5.49 2.18

(Ту 32.9 31.4 33.8 36.2 33.2 32.9 42.7 46.3 42.7 48.8 32.9 39.4

Глуховка-Мостовой 5.86 0.84 0.84 0.63 0.42 6.08 2.16 2.03 5.54 5.14 5.97 2.20

(Ту 17.5 19.7 22.0 19.7 21.8 84.6 89.5 93.2 89.5 93.2 51 56.1

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

При моделировании паводка 1968 г. для практически равнинного водосбора р. Глуховки (рис. 12, г) параметр па получился близким к табличным значениям коэффициентов шероховатости (0.14) назначаемым для покрытых лесом и густым ивняком пойм (Срибный, 1960). Таким образом, можно высказать предположение о том, что действительно имел место сток по «хортоновскому» типу, когда вода, по крайней мере в нижней части склонов, текла слоем.

Для больших и сложных по строению бассейнов (рис. 12, д, е), требуется

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

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

Время, часы Время, часы

а осадки; о измеренные расходы;-К^-СТиН;-усовершенствованная КАУ-вШН.

Рисунок 11. Наблюденный и расчетные гидрографы паводка а), г) - 1621.08.1968 г. и 6), в) -2-5.10.1974; 0 часов на графике соответствует 02:00 16.08.68 и 08:00 02.10.74.

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

В диссертационной работе на основе средств математического моделирования, разработанного комплекса программного обеспечения алгоритмов анализа ЦМР и гидрологической модели, получены следующие научные результаты:

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

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

Усовершенствована модель поверхностного стекания KW-GIUH путем учета пространственной неоднородности параметров водосборов. Разработано программное обеспечение, реализующее алгоритм усовершенствованной версии модели KW-GIUH.

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

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

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

Основные положения диссертации изложены в следующих работах:

1. Модель паводочного цикла речного бассейна в расчетах и прогнозах наводнений // Тезисы докладов VI всероссийского гидрологического съезда. Санкт-Петербург.2004. С 77-78. (в соавторстве с Б.И. Гарцманом и М.А. Макагоновой)

2. Применение модели паводочного цикла малого речного бассейна в расчетах максимального стока // Гидрометеорология и экология Дальнего Востока: Темат. вып. ДВНИГМИ. Владивосток: Дальнаука, 2003. Вып. 4, С. 76-93. (в соавторстве с Б.И. Гарцманом)

3. Прогнозы речного стока в теплый период года на основе модели паводочного цикла. Материалы всероссийской научно-практической конференции «Научные основы экологического мониторинга водохранилищ». Хабаровск, 2005. С. 58-63 (в соавторстве с Б.И. Гарцманом и М.А. Макагоновой)

4. Модель геоморфологического мгновенного единичного гидрографа, основанная на уравнении кинематической волны: теория и применение на Тайване и в России. Материалы XIII научного совещания географов Сибири и Дальнего Востока. Т.1. Иркутск: Изд-во Института географии СО РАН, 2007. С.213-215. (в соавторстве с К.Т.Ли, Н.К.Чен и Б.И.Гарцманом)

5. Цифровые модели рельефа и методы анализа структуры гидрографических систем. Географические и геоэкологические исследования на Дальнем Востоке: сб. науч. Трудов молодых ученых. Вып. 4. Владивосток: Дальнаука, 2008. С. 52-60.

6. Applications of kinematic-wave based géomorphologie IUH model for rain floods prédiction in Taiwan and Russia. Geophysical Research Abstracts, Vol. 10, EGU2008-A-05833, 2008 EGU General Assembly 2008. Vienna (K.T. Lee, N.C. Chen, B.I. Gartsman, A.N. Bugaets).

7. Применение модели геоморфологического мгновенного единичного гидрографа основанной на уравнении кинематической волны для расчета дождевых паводков рек юга Приморья. Материалы Третьей всероссийской научной конференции с международным участием 24-28 августа 2010, Барнаул, С. 42-45.

в том числе в рецензируемых изданиях:

8. Анализ структуры речных систем и перспективы моделирования гидрологических процессов // География и природные ресурсы, 2008. № 2. С. 20-28. (в соавторстве с Гарцманом Б.И., Тегай Н.Д., Краснопеевым С.М.)

9. Краткосрочный прогноз притока воды в водохранилище Бурейской ГЭС // Гидротехническое строительство, 2009. №1. С. 139-145. (в соавторстве с Гарцманом Б.И., Макагоновой М.А., Губаревой Т.С.)

10. Современная версия модели единичного гидрографа и ее применение в Тайване и России // География и природные ресурсы. 2009. - №1. - С. 144151. (в соавторстве с Ли К.Т., Чен Н.К., Гарцманом Б.И.)

Андрей Николаевич БУГАЕЦ

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

Автореферат диссертации на соискание ученой степени кандидата технических наук

Подписано к печати 05.09.2011 г. Печать офсетная. Бумага офсетная. Формат 60x90/16. Усл. п. л. 1,63. Уч.-изд. л. 1,24. Тираж 150 экз. Заказ 97

Отпечатано в типографии ФГУП Издательство «Дальнаука» ДВО РАН 690041, г. Владивосток, ул. Радио. 7

Содержание диссертации, кандидата технических наук, Бугаец, Андрей Николаевич

ВВЕДЕНИЕ.

ГЛАВА 1. МОДЕЛИРОВАНИЕ РЕЧНОГО СТОКА И АНАЛИЗ СТРУКТУРЫ РЕЧНЫХ СИСТЕМ.

1.1. Современные тенденции в разработке бассейновых моделей стока.

1.2. Методические основы анализа структуры речных систем.

ГЛАВА 2. МЕТОДИЧЕСКИЕ АСПЕКТЫ И ПРОГРАММНЫЕ СРЕДСТВА СТРУКТУРНО-ГИДРОГРАФИЧЕСКОГО АНАЛИЗА РЕЧНЫХ СИСТЕМ.

2.1. Анализ речных сетей с использованием ЦМР и ГИС-процедур.

2.2. Пример анализа гидрографической сети бассейна р. Комаровка.

2.3. Программные средства структурно-гидрографического анализа.

ГЛАВА 3. ОПЫТ СТРУКТУРНО-ГИДРОГРАФИЧЕСКОГО АНАЛИЗА И ТИПИЗАЦИИ'РЕЧНЫХ СИСТЕМ БАССЕЙНА АМУРА.

3.1. Структурные закономерности крупных речных систем.

3.2. «Энергетическая функция» как средство типизации речных бассейнов.

ГЛАВА 4'. ВЕРСИИ МОДЕЛИ ГЕОМОРФОЛОГИЧЕСКОГО МГНОВЕННОГО ЕДИНИЧНОГО ГИДРОГРАФА, ОСНОВАННОГО НА УРАВНЕНИИ КИНЕМАТИЧЕСКОЙ ВОЛНЫ (К^-ОЮН).

4.1. Теоретические основы модели К\\МЗШН.ВО

4.2. Усовершенствованная версия модели KW-GIUH.

4.3. Расчет слоя стокоформирования в усовершенствованной версии модели

ГЛАВА 5. ОПЫТ МОДЕЛИРОВАНИЯ ПАВОДКОВ С УЧЕТОМ ПРОСТРАНСТВЕННОЙ НЕОДНОРОДНОСТИ ОСАДКОВ И ПАРАМЕТРОВ

5.1. Объекты и данные моделирования.

5.2. Испытания модели и анализ результатов.

Введение Диссертация по наукам о земле, на тему "Разработка методов определения структурно-гидрографических характеристик по данным ЦМР для гидрологического моделирования"

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

Обобщая, можно сказать,, что весь основной аппарат математического моделирования, разработанный преимущественно в 60-80 годы 20-го века, оказывается малоэффективным в задачах, связанных с экстраполяцией за пределы имеющихся данных наблюдений. Главным препятствием являются трудности корректной параметризации моделей и тривиальность результатов моделирования. Дальнейшее развитие гидрологических моделей, адекватно отражающих сильную нелинейность экстремальной динамики речных систем, особенности их адаптации в условиях нестационарного климата и сложность пространственного строения, требует взаимоувязанного расширения теоретической и информационной базы. Необходимо привлечение нового массива данных, достаточно доступных и надежных, включающих информацию о неизученных бассейнах и редких гидрологических событиях, допускающих органичное развитие теории моделирования - данных о структуре речных систем.

Термин «структурная гидрография» предложен и подробно обоснован Л.М.Корытным (Корьггный, Безруков, 1990). Он определяет его как научное направление, в основе которого лежат структурные закономерности гидрографической системы. Исследования структуры речных систем, основанное на пионерной работе Хортона (1948), широко развивалось в России в 60-80 годах прошлого века, однако во второй половине 90-х они были практически прекращены. В'то же время, за рубежом это направление продолжало развиваться. Фундаментальный труд Родригеса-Итурбе и Ринальдо (Кос^иез-ИшЬе, ШпаМо,

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

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

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

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

2. Разработка и реализация базы, данных структурных, гидрографических и геоморфологических характеристик.

3. Апробация методики анализа ЦМР на примере крупных речных систем в бассейне Амура.

4. Усовершенствование модели геоморфологического мгновенного единичного гидрографа (К"\\^-01ЛН) за счет учета параметров, полученных по индивидуальным структурным элементам речного бассейна.

5. Адаптация^ методики определения эффективных осадков при расчете гидрографов стока с помощью модели К\У-01ЛН.

6. Апробация усовершенствованной модели с использованием данных бывшей

Приморской водобалансовой станции.

Методической основой структурно-гидрографического и морфометрического анализа являются алгоритмы, реализованные на базе ГИС технологий (ERDAS IMAGINE 9.1, ESRI ARCMAP 9.3) и оригинальных программных компонент разработанных в среде DELPHI 7. В качестве базовой СУБД применялась MS Access. Методической основой, используемой при расчете гидрографов дождевого стока, является модель KW-GUIH, разработанная К.Т. Ли с соавторами, прошедшая широкую апробацию за рубежом. Испытание модели в России впервые выполнено автором диссертации. В работе использованы известные методы математического моделирования, а также общепринятые методы статистического и географического анализа.

Объектами исследования являются разнопорядковые речные системы, расположенные в пределах бассейна Амура и прилегающем бассейне северозападного побережья Японского моря. Среднемасштабная ЦМР на всю исследуемую территорию создана на основе данных SRTM (Shuttle Radar Topographic Mission). Для тестовых испытаний различных версий модели KW-GIUH были использованы опубликованные данные бывшей Приморской воднобалансовой станции (ПВБС).

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

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

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

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

Результаты опытного структурно-гидрографического анализа речных систем бассейна Амура и принципиальная схема их типизации на основе распределения суммарного индекса напряжения сдвига (shear stress) от расстояния до замыкающего створа

Усовершенствованная модель геоморфологического мгновенного единичного гидрографа (KW-GUIH).

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

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

Апробация работы. Основные положения и выводы диссертации докладывались автором на следующих конференциях: Всероссийском гидрологическом съезде (СПб, 2004), региональной конференции «Научные основы экологического мониторинга водохранилищ» (Хабаровск, 2005), региональных молодежных конференциях «Географические и геоэкологические исследования на Дальнем Востоке» (Владивосток, 2006, 2008), XIII совещании географов Сибири и Дальнего Востока (Иркутск, 2007), Генеральной ассамблее европейского геофизического союза (Вена, 2008), региональной конференции «Современные геофизические и географические исследования на Дальнем Востоке России» (Владивосток, 2010), всероссийской конференции с международным участием «Фундаментальные проблемы воды и водных ресурсов» (Барнаул, 2010), международной конференции Pan Ocean Remote Sensing Conference (Килунг, 2010). XIV Совещание географов Сибири и Дальнего Востока (Владивосток, 2011).

Результаты работы применялись при выполнении тематики плана НИОКР 1.1.2 Росгидромета («Разработать, испытать и внедрить новые технологии и методы гидрологических прогнозов для рек и водохранилищ России, включая прогнозы наводнений и других опасных гидрологических явлений. Обеспечить выполнение международных обязательств Росгидромета в области гидрологических прогнозов», ДВНИГМИ, 2007). Также они использованы в рамках проекта «Модернизация и техническое перевооружение учреждений и организаций Росгидромета. Совершенствование системы гидрологических наблюдений, в т.ч. в паводкоопасных районах. Ч.З. Модернизация гидрологической сети в бассейне р.Уссури», при подготовке технического задания на разработку региональной гидрологической модели (РГМ), отдельные блоки которой проходят тестовые испытания в отделе гидрологических прогнозов Приморского УГМС.

Диссертация состоит из введения, 5 глав, заключения, списка используемой литературы, включающего 220 наименований и приложений. Работа содержит 215 страниц, основной текст изложен на 151 страницах и включает 55 рисунков и 18 таблиц. Четыре приложения табличного и графического характера помещены, на 64 страницах. По теме диссертации опубликовано 10 работ, из них 3 статьи в журналах, рекомендованных ВАК.

Заключение Диссертация по теме "Гидрология суши, водные ресурсы, гидрохимия", Бугаец, Андрей Николаевич

высоты водосборов

Средние уклоны водосборов в Уссурийском и Михайловском районах Приморского края, включает бассейны р.Комаровки и её правого притока - р.Раковки. Почти со всех сторон указанная территория хорошо очерчена естественными границами и характеризуется большим разнообразием типичных для южного Приморья ландшафтов среднегорья, низкогорья, мелкогорья и Западно-Приморской равнины, и различной степенью антропогенного воздействия на ландшафты. Современный рельеф рассматриваемой территории относится к Суйфуно-Ханкайской депрессии.

1М*в'0"Е 140°0'0"Е

Рис. 2.4. Расположение ПВБС.

Климат территории является типичным для юга Дальневосточного региона, отличается муссонным характером и крайне неустойчивым режимом увлажнения во внутригодовом и многолетнем аспектах. В теплый период года выпадает 85-90% осадков. Для территории характерны тропические циклоны-тайфуны. С ними связанны сильные дожди, выпадающие здесь чаще всего во второй половине лета и начале осени, вызывающие экстремальные дождевые паводки (Краткий физико-географический очерк 1956).

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