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

Ромко, Алексей Викторович

Бесплатный автореферат и диссертация по сельскому хозяйству на тему
Обоснование водного и питательного режима мелиорируемых земель гумидной зоны с использованием интегрированной модели агрогеосистемы
ВАК РФ 06.01.02, Мелиорация, рекультивация и охрана земель

Автореферат диссертации по теме "Обоснование водного и питательного режима мелиорируемых земель гумидной зоны с использованием интегрированной модели агрогеосистемы"

□03057042

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

РОМКО АЛЕКСЕИ ВИКТОРОВИЧ

ОБОСНОВАНИЕ ВОДНОГО И ПИТАТЕЛЬНОГО РЕЖИМА МЕЛИОРИРУЕМЫХ ЗЕМЕЛЬ ГУМИДНОЙ ЗОНЫ С ИСПОЛЬЗОВАНИЕМ ИНТЕГРИРОВАННОЙ МОДЕЛИ АГРОГЕОСИСТЕМЫ

Специальность 06.01.02 - Мелиорация, рекультивация и охрана земель

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

Москва - 2007

003057042

Работа выполнена в отделе природоохранных технологий Всероссийского научно-исследовательского института гидротехники и мелиорации им. А.Н.Костякова (ГНУ ВНИИГиМ Россельхозакадемии)

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

доктор технических наук, профессор Кирейчева Людмила Владимировна Официальные оппоненты:

доктор технических наук Добрачев Юрий Павлович кандидат технических наук, доцент Сухарев Юрий Иванович

Ведущая организация: РГАУ - МСХА им. К.А. Тимирязева

Защита состоится « мая 2007 года в 13 часов на заседании диссертационного совета Д 220.045.01 в Московском государственном университете природообустройства в ауд. 1/201 по адресу: 127550, Москва, ул. Прянишникова, д. 19, факс (495)976-10-46

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

Автореферат разослан «

<3> апреля 2007 года

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

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

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

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

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

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

Для реализации цели поставлены и решены следующие задачи:

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

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

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

4. Сопряжение модели агроценоза с моделью пространственной

геофильтрации.

5. Создание интегрированной модели агрогеосистемы для мелиоративного объекта.

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

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

Адаптация интегрированной модели агрогеоценоза и проведение сценарных исследований выполнены на объекте Куликово Яхромской поймы Московской обл.

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

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

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

- предложена методология проведения сценарных исследований для обоснования мелиоративного режима почв в гумидной зоне с учетом экологических ограничений;

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

Практическая значимость. Разработанная модель может быть использована при обосновании технических решений при реконструкции или модернизации мелиоративных систем гумидной зоны, а также для обоснования агротехнических и мелиоративных мероприятий при интенсивной сельскохозяйственной деятельности. Рекомендации по режиму орошения и нормам внесения азотных удобрений переданы и будут использованы в ЗАО «Куликово» Дмитровского района Московской обл.

Апробация работы. Результаты исследований и основные положения диссертационной работы были доложены на секциях Ученого Совета ВНИИГиМ (2003, 2004, 2005, 2007 гг.), Костяковских чтениях (2005 г.), на Молодежном Форуме «Агробиотехнологии и экологическое земледелие»

(Владимир, 2005).

Публикации. Основные положения диссертационной работы опубликованы в 5 печатных работах,' в том числе 1 в журнале, рекомендованном ВАК.

Объем и структура диссертации. Диссертационная работа состоит из введения, четырех глав и списка литературы из 117 наименований, в том числе 63 работы на иностранном языке. Она изложена на 145 страницах машинописного текста, иллюстрирована 26 рисунками, содержит 22 таблицы.

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

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

Разработка моделей фильтрации и влагопереноса началась более 50 лет назад. Модели создавались разными авторами и институтами и, в зависимости от поставленных задач, моделировали процессы, не предполагая дальнейшего взаимодействия с другими моделями. Это привело к большому количеству моделей с сильно различавшимся набором входных и выходных параметров. Однако, для практических расчетов зачастую требовались результаты моделирования по нескольким моделям, что вызвало необходимость объединения различных моделей в единые системы, позволяющие моделировать более широкий круг взаимосвязанных природных процессов. Для таких систем приходилось писать конверторы входной/выходной информации, чтобы передавать информацию из одной модели в другую, взаимодействие которых друг с другом изначально не предполагалось. В дальнейшем развитие моделей было ' связано с интеграцией с ГИС, базами данных, средствами анализа и представления данных. Стоимость разработки и поддержки программного обеспечения сильно возросла, что усложнило создание сложных моделей отдельными специализированными институтами. Наметились тенденции к стандартизации используемых в моделировании данных и кооперации деятельности различных профильных институтов для создания систем нового поколения. Многие институты отказались от поддержки собственных систем и присоединились к существующим проектам.

В настоящее время достаточно хорошо разработана теория фильтрации в водонасыщенных грунтах. Моделированию процессов геофильтрации посвящены фундаментальные работы П.Я. Полубариновой-Кочиной, В.В. Веригина, А.И. Силина-Бекчурина, В.А. Мироненко, В.М. Шестакова, Д.Ф. Шульгина, Р. Винстона, К. Восса, К. Киппа и многих других. Теоретические основы влагомассопереиоса в почвогрунтах заложены в основополагающих трудах С.Ф. Аверьянова, Е. Букингема, В.В. Ведерникова, H.H. Веригина, А. Дарси, А.Н. Костякова, Л.А. Ричардса и других ученых. Дальнейшее развитие теории осуществлялось в работах И.П. Айдарова, А.И. Голованова,

JI.M. Рекса, В.М. Шестакова, Д.А. Манукъяна, И.О. Пашковского, B.C. Борисова, Муалема, Ван Генухтена, И. Шимунека, Й. Ван Дама и других ученых. Существующие модели фильтрации основываются на законах Дарси и предпосылке неразрывности потока. Различия между отдельными моделями состоят в основном в методе решения (конечные разности, конечные элементы) и в способах задания граничных условий (van Waveren R.H.,2003).

Модели массопереноса в грунтовых водах условно можно разделить на две категории. Первый класс - "общие" модели, в которых химические процессы задаются упрощенно в виде коэффициентов (дисперсия, сорбция и т.д.). К этому классу относятся такие модели как MT3D, МОС и HST3D, разработанные американским геологическим обществом (USGS). Как и в моделях фильтрации, различия между моделями массопереноса лежат в основном в методе решения и в способах задания граничных условий. Второй класс моделей позволяет описывать химические реакции в явном виде. Такие модели нельзя применять без серьезной химической экспертизы. К ним относятся такие модели как модель процессов равновесной химической сорбции (Ю.В. Шваров, 1992) и модель PHREEQC (Parkhurst D.L. 1999, Simunek J., 2004), которая разработана американским геологическим обществом.

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

• Балансовые - наиболее простые модели, не требующие большого количества входных данных, но и не дающие детального решения (van Waveren R.H., 2003). Этот тип моделей применяется, как правило, для того, чтобы получить общее представление о процессах инфильтрации. Примером таких моделей является модель WATBAL, разработанная в университете Вагенингена.

• Псевдо-стационарные - комбинация балансовых и динамических моделей влагопереноса. Примером таких моделей является модель MUST, разработанная в университете Вагенингена (De Laat P.J.M., 1995), имитирующая вертикальный перенос влаги в зоне аэрации с учётом суммарного водопотребления сельскохозяйственных культур и оттока воды в дренаж.

• Динамические модели, в основе которых лежит уравнение Ричардса (van Waveren R.H., 2003). Такие модели требуют достаточно большого количества входных данных, но дают детальные результаты. Этот тип моделей применяется, как правило, в моделях малых масштабов (уровня поля), чтобы получить детальное представление о процессах инфильтрации. Примером таких моделей являются модели А.И. Голованова, JI.M. Рекса, A.M. Якиревича, модель SWAP и др.

Исследование природных систем показывает, что одномерные и псевдодвухмерные имитационные модели не позволяют учесть все возможные процессы, возникающие на мелиорированных

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

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

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

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

Для расчета влагопереноса в интегрированной модели автором используется обобщенная модель, базирующаяся на уравнении Ричардса (Голованов А.И. 1974, Van Dam, J.C., 1997, Kroes, J.G., 2003):

где С - водоемкость , К - коэффициент влагопроводности

(см/сут), в - содержание влаги в почве (см3/см3), Б - отбор влаги корнями растений (см/сут), Ь - напор почвенной влаги (см).

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

S(h),

(1)

растений.

Запасы влаги" рассчитываются по аналитической формуле Ван Генухтена:

е = еге2+е^~9^, (2)

1+НТ

3 3

где £>„1 - содержание влаги в насыщенной почве (см /см ), вк$ -остаточное содержание влаги в почве (см3/см3), а (см'1), п (-) и т (-) -эмпирические параметры.

Коэффициент влагопроводности рассчитывается по формуле Муалема:

^saf-

-(i-S^-JT

(3)

где Ksat - коэффициент фильтрации в насыщенной почве, X -эмпирический параметр (-), Se- относительное насыщение:

Se = 6~6res . (4)

fl — А

"sat °res

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

Отбор влаги корнями растений рассчитывается по зависимости, успешно апробированной на фактических данных, в которой отбор влаги корнями зависит от толщины корнеобитаемого слоя и распределения корневой системы (Van Dam, J.C., 1997, Kroes, J.G., 2003):

Sp(z)= 0Lroot(z) Tp, (5)

iLroot(z)dz

где DM - толщина корнеобитаемого слоя (см); Lr00t - функция, характеризующая распределение корневой системы; Тр - потенциальная скорость транспирации (см/сут).

Недостаток влаги или чрезмерное увлажнение почвы могут уменьшить отбор влаги корнями растений. Действительный отбор влаги корнями Sa(z) (сут"1) рассчитывается из условия:

Sa(z) = drwSp(z), (6)

где d№ - коэффициент снижения отбора влаги корнями из-за водного стресса (-). Этот коэффициент зависит от напора влаги.

В интегрированной модели для расчёта суммарного водопотребления использовано уравнение Пенмана-Монтейта. По результатам сравнения моделей суммарного водопотребления по данным лизиметров с 11 станций во всем мире в разных климатических условиях, эта зависимость была признана лучшей для всех зон и погодных условий (Jensen М.Е., 1990). За

рубежом эта модель является стандартом расчета потенциальной эвапотранспирации, так как она наиболее полно отражает физические процессы, обусловленные солнечной радиацией, аэродинамикой и транспирацией растений. Зависимость имеет следующий вид (Van Dam, J.C., 1997, Kroes, J.G., 2003):

10~4Av(Rn -G)+8.64106pairCair(esat -eacl)-i-^wETp --7---ч-—>. (7)

Лу+Уа

1 +

rcrop

где - энергия парообразования (Дж/гр), ETP - потенциальная транспирация (см/сут), Av - наклон кривой давления пара - поток

солнечной радиации (Дж м"2 сут'1), G - поток почвенного тепла (Дж м"2 сут'1), pair - плотность воздуха (гр/см3), Cair - теплоёмкость воздуха (Дж гр"2 "С"1), esaI - давление насыщенного пара (кПа), eact - фактическое давление пара (кПа), Гсгор - сопротивление листовой поверхности (сек/м), га;г - аэродинамическое сопротивление (сек/м), yair — психометрическая константа (кПа "С"1).

Сопротивление листовой поверхности зависит от скорости ветра и высоты растений (Kroes, J.G., 2003):

ln|bn-d

zh -d

_ ч -ош у ч zoh у /о\

crop =-~2-' W

kvku

где г™ - высота точки измерения скорости ветра (м), zh - высота точки измерения температуры и давления (м), d - точка отсчета профиля ветра (м), Z6m - коэффициент неточности для импульса (м), Zoh - коэффициент неточности для тепла и давления (м), k„k - константа фон Кармана = 0,41, и -скорость ветра на высоте zm (м/с).

Параметры d, zom и z0h определяются как (Kroes, J.G., 2003): 2

d = ~hcrop, (9)

zom=0,123hcrop) (10)

zoh=0,lzom, (11)

где hcrop - высота растений (м). Для расчета аэродинамического сопротивления используется зависимость (8), при этом высота растений принимается равной 1 мм.

Входной информацией для расчета по формуле Пенмана-Монтейта являются среднесуточная температура воздуха, солнечная радиация, скорость ветра и атмосферное давление.

В качестве модели урожайности принята упрощенная зависимость, в которой развитие растения задается через индекс листовой поверхности, высоту растения и глубину корнеобитаемого слоя как линейные по времени функции фазы развития растений. Расчёт урожайности производится по формуле (Van Dam, J.C., 1997, Kroes, J.G., 2003): -

1а. = П Уш i=l

1-Kyi

1 Ек

. ETmJ,

(12)

где j - фаза развития растения; п - число фаз развития; уа - фактический урожай (т/га); ут - максимальный урожай (т/га); ETaj - фактическое суммарное водопотребление в определенную фазу (мм); ETmj - максимальное суммарное водопотребление в определенную фазу (мм); Ку; - фактор урожайности в определенную фазу (-).

Существуют более сложные и точные модели роста и развития растений, такие как модели Ю.П. Добрачева (Добрачев Ю.П., 1998) и В.В. Шабанова (Шабанов В.В., 2005). Однако эти модели требуют большого количества труднополучаемой исходной информации и иногда могут давать физически невозможные результаты.

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

Перенос азотных соединений в почве описывается зависимостями (Николаенко, 1999, 2002, P. Groenendijk, 1999): для органического азота:

5(6CqN) _ dJs,ON R R . nr.

-^----^--Kd,ON + Kp,ON > v '

для аммонийного азота:

d(6CNH4) , aXe,NH4 . 5Js,NH4

-^ +Pd ^ -----Ku,NH4 ~Kd,NH4 +Kp,NH4> (14J

для нитратного азота:

-^----^--Ku,ON3 -Kd,ON3 +Kp,N03>

где Ci - концентрация i-го вещества (кг/мЗ), Ru,i — отбор азота корнями растений (кг м-3 сут-1), Rdj — разложение i-ro вещества (кг м"3 сут"1), Rp i -источник i-ro вещества (кг м'3 сут"1), Jsi — вертикальная миграция i-ro вещества (кг м"2 сут"1), pj - плотность почвы (кг/м3), Хе,ын4 — объем сорбированного аммония (кг/кг), в - содержание влаги в почве (см3/см3).

Разложение азотных соединений характеризуется кинетикой первого порядка (P. Groenendijk, 1999):

Rd = kjGjCj, (16)

где ki — коэффициент распада первого порядка (сут"1). Отбор азотных соединений корнями растений описывается зависимостью (P. Groenendijk, 1999):

Ru=aStLCi, (17)

Д z

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

минерализация органического азота, которая описывается уравнением (Р. Groenendijk, 1999):

1 nf

Rp,NH4 = — i (I(fN,fh -afN,hu)fhkfnOMfn(t) + (fnL -afN)hu)kseC0M + At t fn=l

(fN,ex " afN(hu )kexEX(t) + fN>hukhuHU(t))dt

(18)

где ffjfn - фракция азота в fn-том органическом материале (долей единицы), а - фактор ассимиляции вещества в гумус (долей единицы), f^u -фракция азота в гумусе (долей единицы), fi, - процент разлагающихся твердых органических соединений, идущий на формирование гумуса (долей единицы), kfn - коэффициент распада fn-той органической фракции (сут"1), ОМ - органические соединения в твердой фазе (кг м"3), fNL - фракция азота в растворенных органических соединениях (долей единицы), k<,x, khu, ks -коэффициенты распада остатков корней, гумуса и растворенных органических соединений соответственно (сут"1), EX, HU - количество остатков корней и гумуса соответственно (кг м'3), fN,cx — процент содержания азота в отмирающих частях растений (долей единиц), Сот - концентрация органических соединений в почвенной влаге (кг/мЗ).

Источником поступления нитратов в почву являются процессы нитрификации, которые описываются уравнением (P. Groenendijk, 1999): 1 1„+д'

Rp,N03=— J kniteCNH4(t)dt, (19)

где кпц — коэффициент нитрификации (сут"1). В анаэробных условиях происходит разложение нитратов (денитрификация).

Сорбция аммония в равновесной фазе рассчитывается по формуле (Р. Groenendijk, 1999):

Xe,NH4 =Ke,NH4cnh4> (20)

где K^njm - линейный коэффициент сорбции (м3 кг'1). Предполагается что объем сорбированного вещества пропорционален его концентрации в жидкой фазе.

Для моделирования геофильтрации и геомиграции используется трёхмерная модель движения подземных вод и переноса загрязнений в зоне насыщения nMtWolf, разработанная Санкт-Петербургским отделением ИГЭ РАН (Румянцев, 1997). Модель предназначена для решения трехмерных задач подземной гидродинамики и переноса загрязнений в стационарной и нестационарной постановке при напорной, напорно-безнапорной или безнапорной фильтрации в любом из слоев модели. Возможно подключение набора гидрологических и гидрогеологических объектов, таких как реки различного порядка, горизонтальные и вертикальные скважины, горизонтальные дрены и некоторые другие, на которых реализуются граничные условия II-IV рода. В ее • основу положено трехмерное нестационарное уравнение геофильтрации вида:

где кХ) ку, к2 - коэффициенты фильтрации по осям х, у, г, параллельным главным осям проницаемости; Н - напор; <3 - расход, приведенный к единице объема и представляющий внешние источники поступления и/или отбора воды; т] - упругоемкость; I - время.

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

где [i] - индекс координатной оси (1[1]=х, 1[2]=у, l[3]=z), D[i] -коэффициент гидродисперсии, С и N - концентрация загрязнителя в хорошо и плохо проницаемых зонах, v - скорость фильтрации, а - коэффициент массообмена.

Для совместного моделирования процессов влаго- массопереноса в зоне аэрации и полного насыщения необходимо сопряжение соответствующих моделей. В модели nMtWolf нет верхнего граничного условия в явном виде. Все граничные условия (I, II и III родов) реализуются j на уровне алгоритма через источники/стоки (II и III род), включая и | инфильтрационное питание. Иными словами, то, что поступает в ячейки модели (в верхний слой) сверху - это расход источника. Дополнительный подъем уровня грунтовых вод от действия этого источника определяется уже на этапе решении задачи. Учитывая это обстоятельство, сопряжение моделей зоны аэрации и зоны насыщения было реализовано автором алгоритмически. Имея заданный начальный УГВ в качестве нижнего граничного условия, модель зоны аэрации рассчитывает поток в зону насыщения. Этот поток задается в качестве источника/стока в верхней ячейке модели зоны насыщения. УГВ, рассчитанный в модели зоны насыщения, используется в качестве нижнего граничного условия модели зоны аэрации на следующем временном шаге.

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

В третьей главе «Создание интегрированной модели агрогеосистемы и описание ее алгоритма» описан алгоритм и верификация как отдельных частей, так и интегрированной модели в целом.

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

(22)

расчета подземных вод находится в модели пЖ№о1£ Все остальные процессы реализованы в модели, разработанной автором и названной таЗБ. Управляющий модуль отвечает за согласованную совместную работу модулей зоны аэрации и модели пространственной геофильтрации и геомиграции пМ£^о1£ Он обеспечивает поток данных между ними, формируя файлы с входными данными для пМ£,Лго1£ на основании результатов расчетов процессов в зоне аэрации, запускает ее и считывает из ее выходных файлов результаты расчета. При передаче данных между моделями, управляющий модуль преобразует их единицы измерения в соответствии с принятыми в моделях.

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

Для обеспечения ввода исходной информации для интегрированной модели была разработана оригинальная ГИС, так как многие ГИС созданы для конкретных моделей и их использование для разработанной интегрированной модели не представляется возможным. Есть ГИС разрабатываемые как универсальные. К таким ГИС можно отнести А^шОпе. Указанная ГИС использовалась автором для изучения моделей MODFLOW и МОСЗБ. При этом обнаружилось, что АгдиэОпе неверно создавал входные файлы к указанным моделям - их содержимое не соответствовало отображаемым в ГИС данным. Поэтому была специально разработана ГИС, адекватно работающая для выбранных моделей. Для удобства задания входных параметров в работе был использован аналоговый подход для задания пространственных и временных данных.

Рисунок 1 - Схема моделируемых процессов и их взаим

обрения

.ТИП

(1 -Ык„Рлом„

рг. материя

"С

1 /0+А „у Л» ,п й>=1

Растворенная орг. материя

ч,ои->ни

->ни = сскехЕХ

10 /"=1

<ЛЙОС

к.всп,

/г

Р,оос->ни

. /и+д/

Гумус

ИТ

^^d,нu - ~k.hu чи

Минерализация

Перенос

-О// _ шг,ОК

дг

■Я

■¿.он

■К

■р.оы

/лет/чивание при несении удобрений

Нитрификация\ денитрификация

р.ЫОЪ

1МН4

./„ — а/и.ии (() + — а/д, А„ )к!0Сом +

дсорбированн ый N144

Сорбция\десорбция

^ е,ИН 4 = ^е.Л'Я4<-пЛ4

Перенос

Л/

/ кп„вСт4{(),

N03

Отбор корнями растений

Аг

-► <-

адг.

а/.

эг

,Л7/4 К.<1.НН4 + Кр,Ж4

Перенос

двСцог __ „ г> , о

■'Ч.ОЛГ3 л<*.ОЛГ3 Л/7.Л

Эг

зеи

Рисунок 2 - Блок-схема алгоритма интегрированной модели arpo reo системы

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

азотных соединений с учетом всех моделируемых процессов; моделирование сложных задач для экстремальных условий (интенсивный дождь на сухую почву, сильная эвапорация с мокрой почвы, флуктуация УГВ вблизи поверхности). Максимальные отклонения при верификации модели влагопереноса (исключая задачи для экстремальных условий) составили не более 4%. Максимальные отклонения при верификации модели миграции и трансформации азотных соединений составили не более 8%. Это подтверждает достаточно высокую степень сходства результатов на тестовых задачах с моделями SWAP и ANIMO.

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

Для верификации разработанного модельного комплекса в целом, а также для апробации интегрированной модели был выбран участок К-14 совхоза «Куликовский» Дмитровского района Московской области (рис. 3). Главными особенностями участка, обусловившими его выбор, являются наличие двустороннего регулирования водного режима (орошение и дренаж). Вся территория поймы осушается закрытым горизонтальным дренажем с глубиной 0,8-1,2 м и междренным расстоянием 12, 14 и 18 м. Орошение осуществляется дождевальными агрегатами ДДА-100В. Пойма используется для выращивания овощей и картофеля, кормовых культур и многолетних трав. Использование интенсивных технологий приводит к существенному выносу биогенных загрязнителей в грунтовые воды и реки, что увеличивает связь биологического круговорота воды и веществ с геологическим и вызывает ухудшение экологического состояния агроландшафта, (загрязнение поверхностных и грунтовых вод биогенными веществами). Двустороннее регулирование водного режима позволяет управлять мелиоративным режимом объекта.

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

На рис. 4 приведен литолого-геологический разрез Яхромской поймы. При схематизации гидрогеологических условий в модели выделено пять слоев. За региональный водоупор приняты доледниковые древнеаллювиальные глины.

пжумзда В6 га орошнне ддд-юов

элкрытав проснге1ж>нв9 сета • 2.41 км (ВТ ЗОЛ. 200) стпфы'Е&н оросительная си[ъ • 3,422 ьм отсритэй ос51ютв1*.мя ешь - 2 илнэлэ. 11843 км /-еакрмтая кодлокггормщронажиая сеть ■ коллектора 110 и ]£0 нм, донны 63мм - 2&37 км

Рисунок 3 - Схема участка К-14 совхоза «Куликовский».

Рисунок 4 — Литолого-геологический разрез Яхромской поймы.

Условные обозначения: а1 - аллювиальные суглинки прируслового вала, р] с1 П^1 -покровные и делювиальные пески и суглинки, К - болотные торфяные почвы, I <3^ -озерные отложения сапропеля, а|] I Ош - озерно-аллювиальные и аллювиальные отложения песка суглинка, ¿1 Он т — моренный суглинок Московского оледенения, Гц1[ £1 <2п Йп-ш - моренный суглинок Днепропетровского оледенения и межледниковые флювиогляциальные пески с г-равием и галькой, а! О! - доледниковые древне аллювиальные пески и глины, К - меловые пески, глины, опоки, .Гз - вер хне горе кие глины, С - каменноугольные известняки.

Для обоснования водного и питательного режимов мелиоративного

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

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

2. Расчёт по интегрированной модели водного и питательного режима для всех вариантов.

3. Анализ вариантов, и выбор сценариев при которых вынос загрязняющих веществ не превышает ПДК

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

С целью анализа влияния орошения на ирригационное питание грунтовых вод и вынос биогенных элементов были проведены сценарные исследования режима орошения по интегрированной модели (Ромко A.B., 2006). Моделировались варианты с различными пределами влажности корнеобитаемого слоя, по которым рассчитывалась величина инфильтрационного питания грунтовых вод, урожайность сельскохозяйственных культур и вынос аммония и нитратов в грунтовые воды, реки и дрены. Сценарные исследования проводились для условий 2005 года, характеризующегося как средне-влажный. В указанный год на участке выращивали картофель, посев которого производился 10 мая, уборка урожая — 25 сентября. Параметры для картофеля были получены с Помощью модели WOFOST с учетом фактического урожая картофеля, составляющего 40 т/га. Учитывая тот факт, что это очень высокая урожайность для рассматриваемых условий, это значение было, принято за 100% урожая. Результаты экспериментов приведены в таблице 1.

Таблица 1 - Результаты сценарных исследований на участке К-14

Показатели Варианты расчета

1 2 3 4

Пределы регулирования влажности корнеобитаемого слоя, в долях от НВ 0,6...0,7 0,7...0,8 0,8...0,9 0,9... 1

Относительная урожайность, ^¡ЛУтах 0,90 0,95 0,99 0,99

Влагообмен между зоной аэрации и грунтовыми водами, мм 4,13 4,33 8,62 11,93

Вынос аммония из зоны аэрации в грунтовые воды, кг/га 0,7 0,7 1,3 1,9

Вынос нитратов из зоны аэрации в грунтовые воды, кг/га 4,7 4,8 9,5 12,6

Вынос аммония в. реки, мг с га 6000 6000 6400 7000

Вынос аммония в дрены, мг с га 37 37 75 75

Вынос нитратов в реки, мг с га 97000 97000 101000 106000

Вынос нитратов в дрены, мг с га 758 758 1137 1516

В вариантах с пределом регулирования влажности корнеобитаемого слоя 0,8-0,9 и 0,7-0,8 НВ, влагообмен зоны аэрации с грунтовыми водами стремительно уменьшается по сравнению с вариантом 0,9-1 НВ. В вариантах 0,6-0,7 и 0,7-0,8 ГШ он практически не меняется. В вариантах 0,8-0,9 и 0,7-0,8 НВ вынос аммония и нитратов в грунтовые воды, реки и дрены в сравнении с вариантом 0,9-1 НВ резко снижается. Начиная с варианта 0,7-0,8 НВ эти показатели практически не изменяются. Таким образом, наилучший сценарий - вариант 2, при регулировании влажности в пределах 0,7-0,8 НВ. При этом вынос в грунтовые воды аммония составляет 0,7 кг/га, а нитратов 4,8 кг/га.

Вынос биогенных веществ в дрены существенно ниже выноса в реки. Это объясняется тем фактом, что расчетный период работы дрен в сценарных исследованиях ограничен (не превышает 15 суток).

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

сутки

Рисунок 5 - Расчетные значения максимальных концентраций аммония выносимых на УГВ по вариантам сценарных исследований на участке К-14. Условные обозначения: 1 -для влажности 0,6-0,7 НВ, 2 - для влажности 0,7-0,8 НВ, 3 - для влажности 0,8-0,9 НВ, 4 -для влажности 0,9-1 НВ.

Максимальная концентрации аммония за вегетационный период не превышает 2,5 мг/л. Максимальная концентрации нитратов - 16,6 мг/л. Таким образом, максимальная концентрации аммония превышает ПДК почти в пять раз. При этом максимальная концентрация за весь период практически не опускается ниже ПДК. Максимальная концентрация нитратов не превышает ПДК. Из графиков видно, что максимальные концентрации

выносимого аммония уменьшаются с уменьшением влажности корнеобитаемого слоя. Однако, при снижении влажности ниже 0,7-0,8 НВ концентрация практически не падает.' Максимальные концентрации выносимых нитратов увеличиваются с уменьшением влажности корнеобитаемого слоя. При увеличении влажности, выше 0,8-0,9 НВ концентрация практически не увеличивается.

— 1 — 2

— 3

— 4

1 12 23 34 45 56 67 78 89 100 111 122 133 144 сутки

Рисунок б - Расчетные значения максимальных концентраций нитратов выносимых на УГВ по вариантам сценарных исследований на участке К-14. Условные обозначения: 1 -для влажности 0,6-0,7 НВ, 2 - для влажности 0,7-0,8 НВ, 3 - для влажности 0,8-0,9 НВ, 4 -для влажности 0,9-1 НВ.

За счет снижения влажности корнеобитаемого слоя происходит некоторое снижение урожайности картофеля - от 10% для 0,6-0,7 НВ до 1% при 0,8-0,9 НВ (табл. 1). Результаты моделирования показали, что в вариантах с меньшей влажностью при одинаковых нормах внесения удобрений, отбор азота корнями растений был выше. Это можно объяснить тем фактом, что при меньшем содержании влаги концентрация азотных соединений в почвенном растворе выше. Таким образом, при одинаковом уровне транспирации количество поглощенного азота будет выше, что может отразиться на качестве продукции.

Динамика влажности корнеобитаемого слоя для расчетного 2005 года приведена на рисунке 7. Большую часть вегетационного периода влажность колебалась около значения НВ. Увеличение до 1,1 НВ в середине сезона обусловлено поливами, снижение до 0.8 НВ в конце сезона обусловлено сухой и жаркой погодой в сентябре. Расчеты показывают, что при таком режиме влажности за вегетационный период влагообмен составил 16,28 мм, вынос аммония ,и нитратов - 2,9 и 15,65 кг/га соответственно. Среднее расчетное значение концентрации аммония, выносимого непосредственно в

и,и |ии

0,0158 0,0156 0,0154 0,0152 0,0150 0,0148 0,0146 0,0144 0,0142 0,0140 0.0138

реки, составило 1,05 мг/л. Это соответствует измеренным значениям - 0,93 мг/л (Стрельбицкая Е.Б., 2003).

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

4,0000 3,5000 3,0000 г 2,5000

О 2.0000 +

О 1,5000 1,0000 0,5000 0,0000

" А "

1.10 1,05 1,00 0.95 0,90 0,85 0,80 0,75

IОс+Ор -Влажность

1 15 29 43 57 71 85 99 113127141 155 сутки

Рисунок 7 - Динамика влажности корпсобитаемого слоя и 2005 году. Условные обозначения: Ос+Ор - осадки и орошение (мм/сут), влажность (долей ИВ).

На основе созданной модели участка Куликовский были проведены прогнозные расчеты выноса биогенных веществ с мелиоративного объекта на 2006 год. В качестве исходных данных были использовапы фактические погодные условия 2006 года, фактические нормы и сроки внесения минеральных удобрений и проведения поливов (табл. 2). 2006 год характеризовался как влажный (обеспеченность по осадкам 83%).

Таблица 2 - Агрономические данные по участку Куликовский за 2006 год

Выращи- Внесение Посев Поливы Уборка

ваемая удобрений урожая

культура {площадь, Виды и нормы (кг/га) сроки нормы сроки Нормы (мм) Сроки величина, т/га сроки

Капуста Азофоска 15.04- 33 тыс. 28.04- 30 5.05 95 октяб

(74) -600 Туран -0,02 (HHj - 42 кг д.в., NOj - 60 30.05 шт/га 30.05 30.05 15.06 10.07 рь 1

кг д.в.)

Морковь (24) . Азофоска -600 Туран -0,02 (Ш4-42 кг д.в., ИОз - 60 кг д.в.) Апре ль-май 1 млн. шт./га 1.055.05 30 5.06 25.06 • 15.07 55 05.0915.09

Картофель (30) Диаммофо ска-600 (ЫН4-30 кг д.в., Ы03-30 кг д.в.) Май 2,8 т/га 12.0514.05 41 10.0915.09

Прогнозные расчеты показали, что при фактическом мелиоративном режиме влажности в 2006 году влагообмен между зоной аэрации и грунтовыми водами составил 25 мм, вынос аммония и нитратов с моделируемого участка - 4,8 и 17,7 кг/га соответственно (при внесении д.в. 42 и 60 кг/га соответственно). Однако, с экологических позиций это не самый благоприятный вариант. Выполненные прогнозные расчеты показали, что при снижении урожайности на 5% представляется возможным уменьшить , величину влагообмена до 6,9 мм за вегетационный период, при этом вынос аммония и нитратов снизится до 1,1 и 7,5 кг/га соответственно. Таким образом, при применении оптимального водного режима можно снизить -вынос аммония на 3,7 кг/га, а нитратов на 10,2 кг/га, соответственно на 77 и 57%. Применение оптимального режима позволило бы снизить максимальную концентрацию выносимого аммония с 2,8 до 2,5,мг/л (рис. 8), а максимальную концентрацию нитратов с 17 до 16 мг/л (рис. 9). Однако, следует отметить, что даже при таком режиме, максимальная концентрация аммония превышает ПДК в пять раз.

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

п 0,0030 2

* 0,0025 о;

сутки

Рисунок 8 - Расчетные значения максимальных концентраций аммония, выносимых на УГВ. Условные обозначения: 1 - для фактической, 2 - для влажности 0,7-0,8 НВ.

Рисунок 9 - Расчетные значения максимальных концентраций нитратов, выносимых на УГВ. Условные обозначения: 1 - для фактической, 2 - для влажности 0,7-0,8 НВ.

Выводы и рекомендации

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

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

3. Обоснован выбор математических моделей для моделирования процессов в зоне аэрации: влагопереноса, переноса и трансформации азотных соединений, водопотребления и урожайности. Для моделирования пространственной геофильтрации и массопереноса, была использована трёхмерная модель пМ1\Уо1Г, разработанная Санкт-Петербургским отделением ИГЭ РАН, которая предназначена для решения трехмерных задач подземной гидродинамики и переноса загрязнений в стационарной и нестационарной постановке при напорной, напорно-безнапорной или безнапорной фильтрации.

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

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

6. Проведена верификация как отдельных компонентов интегрированной модели, так и модели зоны аэрации в целом. Для этого использовались задачи, применяемые для верификации и оценки моделей, на основании концепций которых была построена разработанная модель. Предельные отклонения в расчетах испарения и инфильтрации на УГВ составили не более 4%. Предельные отклонения концентраций аммония и нитратов составили не более 8%.

7. Выполнены сценарные исследования на опытном участке Яхромской поймы и даны рекомендации о выборе водного режима, поддерживающего влажность почвы в диапазоне 0,7-0,8 НВ, что позволяет получать высокие урожаи (картофеля до 40 т/га, капусты до 90 т/га) при достаточно низких ' уровнях загрязнения грунтовых вод. При применении рекомендуемого

водного режима можно снизить вынос в грунтовые воды: аммония на 3,7 кг/га, а нитратов на .10,2 кг/га (соответственно на 77 и 57%) за вегетационный период по сравнению с применяемой технологией.

По материалам диссертации опубликованы следующие работы:

1. Ромко A.B., Применение математических моделей для управления агроценозом на мелиорированных землях, Вопросы мелиорации № 1-2, 2005 г. с. 99-103.

2. Яшин В.М., Ромко A.B., Формирование дренажного стока с мелиорируемых земель Яхромской поймы, Материалы молодежного форума "Агробиотехнологии и экологическое земледелие" (13-16 апреля 2005 г., г. Владимир), Владимир 2005, с. 102-104.

3. Ромко A.B., Создание интегрированной модели агрогеоценоза на мелиорированных землях, Материалы международной конференции "Наукоемкие технологии в мелиорации", М. 2005, с. 385-389.

4. Ромко A.B., Перенос и трансформации азотных соединений на мелиорируемых землях, Вопросы мелиорации № 1-2, М. 2006, с. 92-100.

5. Ромко A.B., Обоснование режима орошения и норм азотных удобрений, Агрохимический вестник №5, 2006, с. 16-17.

Подписано к печати 09.04.2007 г. Формат 60x84 1/32 Тираж 100. Зак. № 33 ГНУ ВНИИГиМ Россельхозакадемии 127550, Москва, Б. Академическая, 44

Содержание диссертации, кандидата технических наук, Ромко, Алексей Викторович

Введение.

Глава 1 Существующие методы моделирования мелиоративного режима почв, модели фильтрации и влагомассопереноса для ненасыщенно-насыщенной зон.

1.1 Модели фильтрации в насыщенной зоне.

1.2 Модели массопереноса в грунтовых водах.

1.3 Модели влаго- массопереноса в зоне аэрации.

1.4 Псевдо-стационарные и имитационные модели влагопереноса и массопереноса в насыщенно-ненасыщенных почвах и грунтах.

1.5 Сложные имитационные модели.

1.6 Интеграция моделей и построение моделирующих комплексов.

1.7 Графические оболочки.

1.7 Базы данных.

1.8 Средства анализа и представления данных.

1.9 Средства для решения обратных задач.

Выводы по 1 главе.

Глава 2 Выбор математических зависимостей для создания интегрированной модели агрогеосистемы в условиях двухстороннего регулирования водного режима гумидной зоны.

2.1 Краткая характеристика гумидной зоны.ч

2.2 Создание интегрированной модели агрогеосистемы для гумидной зоны

2.3 Обоснование выбора математических моделей.

2.3.1 Расчёт движения влаги в зоне аэрации.

2.3.2 Расчёт суммарного водопотребления.

2.3.3 Расчёт урожайности сельскохозяйственных культур.

2.3.4 Перенос азотных соединений.

2.3.5 Фильтрация и массоперенос в подземных водах.

Математическая модель фильтрации.

Модель переноса загрязняющих веществ в грунтовых водах.

2.3.6 Сопряжение моделей зоны аэрации и зоны полного насыщения.

Выводы по 2 главе.

Глава 3 Создание интегрированной модели агрогеосистемы и описание ее алгоритма.

3.1 Описание ГИС интегрированной модели.

3.2 Верификация модели зоны аэрации.

Выводы по 3 главе.

Глава 4 Моделирование водного и питательного режима с использованием интегрированной модели для Яхромской поймы Московской области.

4.1 Выбор объекта исследований.

4.2 Краткая характеристика Яхромской поймы.

4.3 Построение фильтрационной модели мелиоративного объекта на участке Куликовский.

4.4 Построение модели зоны аэрации участка Куликовский.

4.5 Калибровка интегрированной модели участка Куликовский.

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

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

Выводы по 4 главе.

Введение Диссертация по сельскому хозяйству, на тему "Обоснование водного и питательного режима мелиорируемых земель гумидной зоны с использованием интегрированной модели агрогеосистемы"

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

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

Для реализации цели поставлены и решены следующие задачи:

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

4. Сопряжение модели агроценоза с моделью пространственной геофильтрации.

5. Создание интегрированной модели агрогеосистемы для мелиоративного объекта.

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

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

Адаптация интегрированной модели агрогеоценоза выполнена на объекте Куликово Яхромской поймы Московской обл.

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

- выполнено обоснование мелиоративного режима для объекта Куликово Яхромской поймы Московской обл.

Апробация работы. Результаты исследований и основные положения диссертационной работы были доложены на секциях Ученого Совета ВНИИГиМ (2003, 2004, 2005 гг.), Костяковских чтениях (2005 г.), на Молодежном Форуме «Агробиотехнологии и экологическое земледелие» (Владимир, 2005).

Публикации. Основные положения диссертационной работы опубликованы в 5 печатных работах, в том числе 1 в журнале, рекомендованном ВАК.

Объем и структура диссертации. Диссертационная работа состоит из введения, четырех глав и списка литературы из 118 наименований, в том числе 64 работы на иностранном языке. Она изложена на 145 страницах машинописного текста, иллюстрирована 26 рисунками, содержит 22 таблицы.

Заключение Диссертация по теме "Мелиорация, рекультивация и охрана земель", Ромко, Алексей Викторович

Основные выводы

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

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

3. Обоснован выбор математических моделей для моделирования процессов в зоне аэрации: влагопереноса, переноса и трансформации азотных соединений, водопотребления и урожайности. Для моделирования пространственной геофильтрации и массопереноса, была использована трёхмерная модель пШ\¥оШ, разработанная Санкт-Петербургским -отделением ИГЭ РАН, которая предназначена для решения трехмерных задач подземной гидродинамики и переноса загрязнений в стационарной и нестационарной постановке при напорной, напорно-безнапорной или безнапорной фильтрации.

5. Разработана ГИС, обеспечивающая ввод входных данных для интегрированной модели. Для удобства задания входных параметров в работе был выбран вариант с использованием аналогового подхода для задания пространственных и временных данных.

6. Проведена верификация как отдельных компонентов интегрированной модели, так и модели зоны аэрации в целом. Для этого использовались задачи, применяющиеся для верификации и оценки моделей, на основании 4 концепций которых была построена разработанная модель. Предельные отклонения в расчетах испарения и инфильтрации на УГВ составили не более 4%. Предельные отклонения концентраций аммония и нитратов составили не более 8%.

7. Выполнены сценарные исследования на опытном участке Яхромской поймы и даны рекомендации о выборе водного режима, поддерживающего влажность почвы в диапазоне 0,7-0,8 HB, что позволяет получать высокие урожаи (картофеля до 40 т/га, капусты до 90 т/га) при достаточно низких уровнях загрязнения грунтовых вод. При применении рекомендуемого водного режима можно снизить вынос в грунтовые воды: аммония на 3,7 кг/га, а нитратов на 10,2 кг/га (соответственно на 77 и 57%) за вегетационный период.

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

Библиография Диссертация по сельскому хозяйству, кандидата технических наук, Ромко, Алексей Викторович, Москва

1. Аверьянов С.Ф., Борьба с засолением орошаемых земель. М., "Колос", 1978.

2. Айдаров И.П. Очерки по развитию орошения в СССР и России. Изд. МГУП, М. 2006, 269с.

3. Айдаров И.П. перспективы развития комплексных мелиораций в России. Изд. МГУП,М., 2004, 137с.

4. Агроэкология / В.А. Черников, P.M. Алексахин, A.B. Голубев и др.; Под ред. В.А. Черникова, А.И. Черкеса. М.: Колос, 2000. - 536 с.

5. Бахвалов Н.С. Численные методы / Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков. 3-е изд., перераб. и доп. - М.: БИНОМ. Лаборатория знаний, 2003. - 632 е., илл.

6. Бондарик И.Г., Корягин В.А. Моделирование процессов миграции и * трансформации азота в системе почва-растение-атмосфера, современные проблемы мелиораций и пути их решения. Том II (99).- М.: 1999.- 367 с.

7. Борисов B.C. Моделирование водно-солевого режима мелиорируемых земель на основе плоских задач теории массопереноса Дисс. на соис. к.т.н. М. 1983

8. Борисов B.C., Манукьян Д.А. О моделировании водного и солевого режимов осваиваемых территорий при автоматизаци проектирования, ГМС. Сб. научных трудов "Автоматизация подготовки и управления водохозяйственым строительством", М., ВНИИГиМ, 1985, с. 110-119

9. Ю.Борисов B.C., Манукьян Д.А. Совершенствование модели водного режима134на мелиорируемых землях. Труды ВНИИГиМ, № 75 "Теория и практика мелиораций", М., 1989, с. 193-199.

10. П.Веригин H.H., Шестаков В.М. Методы расчета движения грунтовых вод в двухслойной среде. Изд. ВНИИ ВОДГЕО, Госстойиздат, М. 1954.

11. Голованов А.И. Оптимизация режимов орошения черноземов. "Почвоведение", 1993, №6

12. Голованов А.И. Прогноз водно-солевого режима и расчет дренажа на орошаемых землях. Дисс. на соис. д.т.н. М. 1974

13. Голованов А.И., Зимин Ф.М. Природообустройство (курс лекций). Изд. МГУП. М. 2000,149с.

14. Голованов А.И., Кожанов Е.С., Сухарев Ю.И. Ландшафтоведение.-М.: Колос, 2005.-216 с.

15. Голованов А.И. Прогноз водно-солевого режима и расчет дренажа на орошаемых землях. Дисс. на соис. д.т.н. М. 1974

16. Громов A.B. Моделирование процессов миграции соединений азота в дерново-подзолистых почвах при поливах животноводческими стоками. Дисс. на соис. к.т.н. М. 1991

17. Добрачев Ю.П. Теория и технология управления орошением на основе эколого-физиологических моделей. Дисс. на соис. уч. ст. д.т.н. М.: 199819.3айдельман Ф.Р. Эколого-мелиоративное почвоведение гумидных ландшафтов. ВО-Агропромиздат, 1991, с.320.

18. Карпенко Н.П., Манукьян Д. А, Кирейчева JI.B., Методология ансамблевого прогнозирования состояния природно-мелиоративных систем, Мелиорация и водное хозяйство №4,2004, с. 44-46

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

20. Мелиоративная энциклопедия том 2, ФГНУ «Росинформагротех», М.,2004

21. Мироненко В. А., Румынии В. Г. Проблемы гидроэкологии: Монография: В 3 т. Т. 1. Теоретическое изучение и моделирование геомиграционныхпроцессов. М.: Изд-во Моск. Горн, ун-та, 1998 - 610с.: ил.

22. Мироненко В.А., Проблемы численного моделирования переноса загрязнений в подземных водах, Известия Высших учебных заведений, Геология и разведка, 1990 N5.

23. Мироненко В.А., Шестаков В.М. Теория и методы интерпретации опытно-фильтрационных работ. М., "Недра", 1978. 325 с.

24. Николаенко А.Н. Математическая модель миграции и трансформации азотных соединений в системе почва-раствор-растение. Современные проблемы мелиораций и пути их решения. Том 2 (99).-М.:1999. с 134-154.

25. Николаенко А.Н. Математическое описание физико-химических процессов азотно-солевого комплекса почв. Доклады Российской академии сельскохозяйственных наук. 2002 N2. с 38-41.

26. Никольский Б.П. Методы исследования адсорбции почвами из растворов. Современные методы исследования физико-химических свойств почв. Изд. АН СССР. 1948.

27. Пашковский И.С. Методы определения инфильтрационного питания по расчетам влагопереноса в зоне аэрации Изд. МГУ, 1973

28. Плюснин И.И., Голованов А.И. Мелиоративное почвоведение (учебник). "Колос", М., 1983

29. Полубаринова-Кочина П.Я. Теория движения грунтовых вод. М., "Наука", 1977

30. Пчелкин В.В., Обоснование мелиоративного режима на осушаемых пойменных землях. М.: КолосС, 2003.253 с.

31. Рекс JI.M., Системные исследования мелиоративных процессов и систем. -М.: Издательство "Аслан", 1995. с. 192.

32. Ромко A.B., Обоснование режима орошения и норм азотных удобрений, Агрохимический вестник №5,2006, с. 16-17.

33. Ромко A.B., Создание интегрированной модели агрогеоценоза на мелиорированных землях, Материалы международной конференции "Наукоемкие технологии в мелиорации", М. 2005, с. 385-389.

34. Румянцев В.А. Отчет по НИР "Создание объектно-ориентированных средств компьютерного моделирования", ООО "BEXEST", Санкт-Петербург, 1997 г., 210 стр.

35. Румянцев В.А., Другов Д.В. Научно-методическое сопровождение -построения модели фильтрации и переноса загрязнений на объекте <Куликовский>, Информационный отчет по теме: № 07/03, Санкт-Петербург 2003.

36. Сергеев А.И., Математическое моделирование влияния орошаемого земледелия на гидрологический и гидрохимический режим рек, Автореферат дисс. на соискание уч. ст. д.т.н., Ташкент 2001

37. Силин-Бекчурин А. И., Динамика подземных вод, 2 изд., М., 1965

38. Сотнева Н.И. Оценка направленности водно-солевого процесса в почвахсеверного Прикаспия и возможность их выборочного орошения. Дисс. на соис. к.т.н. М. 2004

39. Стрельбицкая Е.Б. Влияние осушительных мелиораций на процессы эвтрофикации малых рек. Дисс. на соис. к.б.н. М. 2003

40. Тооминг Х.Г., Солнечная радиаци и формирование урожая. Ленинград: Гидрометеоиздат, 1977. 200 с.

41. Трифонов В.А. Вынос химических веществ дренажным стоком с осушаемых пойменных почв и его регулирование, Диссертация на соискание уч. ст. ктн, М. 1989.

42. Шабанов В.В., Обоснование технических систем для восстановления природного состояния автотрофного звена экосистемы водосбора, Мелиорация и водное хозяйство №5, М.: 2005, с. 33-39.

43. Шестаков В.М. Динамика подземных вод. М. МГУ, 1979

44. Шестаков В.М. Теоретические основы оценки подпора, водопонижения и дренажа Изд. МГУ, 1965 *

45. Шишов К.Н., Панов Е.П. Изменение агрохимических свойств торфяных почв под влиянием осушения и длительного сельскохозяйственного использования, Тр. ВНИИГиМ, том 51,1972, с. 91-95.

46. Якиревич А.М. Оценка качества воды методом математического моделирования. Мелиорация и водное хозяйство, 1989 N10.

47. Яшин В.М., Загрязнение дренажного стока с мелиорируемых пойм,

48. Материалы международной конференции "Наукоемкие технологии в мелиорации", М. 2005, с. 457-463.

49. Яшин В.М., Ромко А.В., Формирование дренажного стока с мелиорируемых земель Яхромской поймы, Материалы молодежного форума "Агробиотехнологии и экологическое земледелие" (13-16 апреля 2005 г., г. Владимир), Владимир 2005, с. 102-104.

50. Anderman, E.R., Hill, М.С., MODFLOW-2000, The U.S. Geological Survey ^ Modular Ground-Water Model Documentation of effective-porosity parameters in the advective-transport observation (ADV2) package. Denver Colorado 2003,9 p.

51. ArgusONE User's Guide. Argus Open Numerical Environment A GIS

52. Modelling System, Version 4.0, Jericho, NY 1998,320 p.

53. Blind, M.W., R.H. van Waveren, A. Bresser, P. Gijsbers, T. van der Wal, H. Scholten, J. Noort, "Application of models in water management in The Netherlands: past, present and future", Proceedings Scientific and Technical Review EurAqua, Madrid 2000.

54. Bogdanova T.A. Antropological Impact on Groundwater System. GIS based Groundwater Modelling, The Netherlands 8 p.

55. Danuso F., Gani M. e Giovanardi R., 1995. Field water balance: BIdriCo 2. In "Crop Water models in practice" L.S. Pereira, B.J. Van den Broek, P. Kabat and R.G. Allen, pag. 49-73, Wageningen Pers (ISBN 90-74134-26-2).

56. De Laat, P.J.M. 1995. Design and operation of a subsurface irrigation scheme with MUST. In Pereira, L.S., B.J. van den Broek, P. Kabat and R.G. Allen (Editors). Crop-water-simulation models in practice. Wageningen Pers, The Netherlands: 123-140.

57. Dogrul E.C. and Kadir T., Integrated Groundwater-Surfacewater Model 2 (IGSM2), International Groundwater Modelling Center, Volume XXII ISSUE 1, page 2

58. Ellingson C. and Schwartzman P., Integration of a Detailed Groundwater Model into a Regional HSPF Model, International Groundwater Modelling Center, Volume XXII ISSUE 1, page 3

59. Giudici M., Inverse Modelling for Flow and Transport in Porous Media, Lecture given at the College on Soil Physics Trieste, 3-21 March 2003

60. Goode, D.J. Age, Double Porosity, and Simple Reaction Modifications for the MOC3D Ground-Water Transport Model U.S. Geological Survey Water-Resources Investigations Report 99-4041, Lemoyne, Pennsylvania 1999,34 p.

61. Graham D., Integrated Groundwater\Surface Water Modeling with MIKE SHE,1.ternational Groundwater Modelling Center, Volume XXII ISSUE 1, page 5

62. Groenendijk P. and J.G. Kroes, 1999. Modelling the nitrogen and phosphorus leaching to groundwater and surface water with ANIMO 3.5. Wageningen (The Netherlands), Winand Staring Centre. Report 144. 138 pp.; 50 Figs; 13 Tables; 77 Refs.

63. Guo, W., and Langevin, C.D., User's guide to SEA WAT: A computer program for simulation of three-dimensional variable-density ground-water flow: Techniques of Water-Resources Investigations of the U.S. Geological Survey, Book 6, Chapter A7,77 p.

64. Harbaugh, A.W. and McDonald, M.G. Programmer's Documentation for MODFLOW-96, an update to the U.S. Geological Survey Modular Finite-Difference Ground-Water Flow Model U.S. Geological Survey Open-File Report 96-486, Reston, Virginia 1996,220 p.

65. Harbaugh, A.W. and McDonald, M.G. Users's Documentation for MODFLOW-96, an update to the U.S. Geological Survey Modular Finite-Difference Ground-Water Flow Model U.S. Geological Survey Open-File Report 96-485, Reston, Virginia 1996, 56 p.

66. Hill, M.C., Methods and Guidelines for effective model calibration U.S. Geological Survey Water-Resources Investigations Report 98-4005 Denver, Colorado 1998,90 p.

67. Jan W. Hopmans and colleagues Estimation of Multi-Dimensional Root Water and Nutrient Uptake Rates under Time-Varying Soil Moisture Conditions, Advances Agronomy, vol 77, University of California, Davis, 2002

68. Jensen, M.E., R.D. Burman and R.G. Allen, 1990. Evapotranspiration and irrigation water requirements. ASCE mauals and reports on engineering practice 70, ASCE, New York. 332 p.

69. Kipp K.L., Jr. Guide to the Revised Heat and Solute Transport Simulator: HST3D Version 2, U.S. Geological Survey Water-Resources Investigations Report 97-4157, Denver Colorado, 1997,149 p.

70. Kipp, K.L., Konikow, L.F., and Hornberger, G.Z. An Implicit Dispersive Transport Algorithm for the U.S. Geological Survey MOC3D Solute-Transport Model U.S. Geological Survey Water-Resources Investigations Report 984234, Reston, Virginia 1998, 54 p.

71. Kireycheva L. Mathematical models of soil irrigation and salting, Encyclopedia of Life Support Systems EOLSS, vol 6, Chapter 3., 2001

72. Konikow, L.F., Goode, D.J., Hornberger, G.Z. A Three-Dimensional Methof of Characteristics Solute-Transport Model (MOC3D), U.S. Geological Sur-vey Water-Resources Investigations Report 96-4267, Reston, Virginia 1996, 87 p.

73. K00I, J.B., and M.Th. van Genuchten HYDRUS, One-dimensional variably satuated flow and transport model including hysteresis and root water uptake. Research Report 124, U.S. Salinity Laboratory, USDA, ARS, Riverside, CA 1991.

74. Kool, J.B., and M.Th. van Genuchten, 1991. HYDRUS, One-dimensional variably saturated flow and transport model including hysteresis and root water uptake. Research Report 124, U.S. Salinity Laboratory, USDA, ARS, Riverside, CA.

75. Kroes, J.G. and J. Roelsma, 1997. ANIMO 3.5; User's guide for the ANIMO version 3.5 nutrient leaching model. Wageningen, the DLO Winand Staring

76. Centre. Technical Document 46; 98 pp.; 15 tab.; 16 fig.; 27 ref.

77. Kroes, J.G. and J.C. Van Dam (eds), 2003. Reference Manual SWAP version 3.0.3. Wageningen, Alterra, Green World Research. Alterra-report 773. Reference Reference Manual SWAP version 3.0.3.doc. 211 pp. 39 figs.; 6 tables.; 17 appendices.

78. Kronvang В., Larsen S.E., Jensen J.P., Andersen H.E., de Miguel Fernandez I., Gomez S., Catchment Report: Guadiamar, Spain trend analysis, retention and source apportionment, EUROHARP 11-2005,27 p.

79. Kyllmar K., Nitrogen Leaching in Small Aricultural Catchments Modelling and monitoring for assessing state, trends and effects of counter-measures. Doctoral thesis, 2004.

80. Mualem, Y. 1976. A new model for prediciting the hydraulic conductivity ofunsaturated porous media. Water Resour. Res. 12: 513-522.

81. Poeter E.P. and Hill, M.C. Documentation of UCODE, A Computer Code for Universal Inverse Modelling, U.S. Geological Survey Water-Resources Investigations Report 98-4080, Denver, Colorado 1998,116 p.

82. Poeter E.P., Hill M.C., Banta E.R., Mehl S., Christensen S. UCC)DE2005 and six other computer codes for universal sensivity analysis, calibration, and uncertainty evaluation, U.S. Geological Survey, Reston, Virginia: 2005,283 p.

83. Pollock D.W., Users's Guide for MODPATH/MODPATH-PLOT, Version 3: A particle tracking post-procesing, package for MODFLOW, the U.S. Geologi-cal

84. Survey finite-difference ground-water flow model, U.S. Geological Survey Open-File Report 94-464, Reston, Virginia 1994,214 p.

85. Prucha R.H., SOFTWARE REVIEW MIKE SHE, International Groundwater Modelling Center, Volume XXII ISSUE 1, page 5

86. Schoumans F.S., Silgram M., Review and Literature Evaluation of Quantification Tools for the Assessment of Nutient Losses at Catchment Scale. EUROHARP Report no 1-2003, NIVA Report SNO 4739-2003, Oslo, Norway, 120 p.

87. Silgram M., Schoumans O.F. Modelling approaches: Model parametrisation, calibration and performance assessment methods in the EUROHARP project, EUROHARP 8-2004,18 p.

88. Simunek J., van Genuchten M.Th., Jacques D., Schaap M. and Mattson E.D. Recent Development in the Hydrus Software: Overland Flow and Biogeochemical Modules, 2004.

89. Simunek J., van Genuchten M.Th., Sejna M., Toride N., and Leij F.J. The STANMOD Computer Software for Evaluating Solute Transport in Porous Media Using Analytical Solutions of Convection-Dispersion Equation, Riverside, California 1999

90. Solo-Gabriele H. and Nemeth M., MODBRANCH: A means to couple Surface Water Model (BRANCH) with a Groundwater Model (MODFLOW), International Groundwater Modelling Center, Volume XXII ISSUE 1, page 4

91. SWAN, электронная документация, Wageningen, 1998

92. Teixeira J.L., Fernando R.M., Pereira L.S. 1995. RELREG: a model for real time irrigationscheduling. .In: Pereira L.S., van den Broek B, Kabat P, Allen RG (eds), Crop-WaterSimulation Models, in Practice, pp 3-15, 1995,

93. Wageningen Pers, Wageningen

94. Tiktak A., Beusen A.H.W., Boumans L.J.M., Groenendijk P., de Haan B.J., Portielje R., Schottem C.G.J., Wolf J., Toets van STONE versie 2.0, RIVM rapport 718201007/2003, Bilthoven 2003.

95. Tiktak A., F. Van den Berg, J.J.T.I. Boesten, D. Van Kraalingen, M. Leistra and A.M.A. van der Linden., RIVM Report 711401 008, Alterra Report 28. Manual of FOCUS PEARL version 1.1.1; November 2000.

96. Tiktak A., van der Linden A.M.A. and Boesten J.J.T.I. The GeoPEARL model Model description, applications and manual RIVM Report 716601007/2003 Wageningen 2003, 79 p.

97. Van Dam J.C., Field-scale water flow and solute transport. SWAP model concepts, parameter estimation and case studies. Doctoral Thesis Wageningen University, 2000.

98. Van Dijk A., Moene A.F., de Bruin H.A.R., The principles of surface flux physics: theory, practice and description of the ECPACK library, Internal report 2004/1, Wageningen 2004,98 p.

99. Van Waveren R.H., Groot S., Scholten S., van Geer F.C., Wosten J.H.M., Koeze R.D., Noort J.J. Good Modelling Practice Handbook, STOWA Report 99-05.

100. VanderKwaak J. and Loague K., InHM (Integrated Hydrology Model), International Groundwater Modelling Center, Volume XXII ISSUE 1, page 2

101. Voss C.I., Boldt D., and Shapiro A.M., A Graphical-User Interface for the U.S. Geological Survey's SUTRA Code using ArgusONE (for Simulation of Variable-Density Saturated-Unsaturated Ground-Water Flow with Solute or

102. Energy Transport) U.S. Geological Survey Open-File Report 97-421 Reston, Virginia 1997,106 p.

103. Vrugt J.A., van Wijk M.T., Hopmans J.W. and Simunek J. One-, two-, and three-dimensional root water uptake functions for transient modelling, WATER RESOURCES RESEARCH, VOL. 37, NO. 10, PAGES 2457-2470, OCTOBER 2001

104. WASY Aktuell, Volume 10, No. 3 September 2004

105. Winston R.B. A Graphical-User Interface for MODFLOW and MT3D 1.5 Using ArgusONE, DRAFT: beta 0.92, New York, 1998.

106. Winston R.B. Graphical User Interface for MODFLOW version 4. U.S. Geological Survey Open-File Report 00-315, Reston, Virginia 2000,27 p.

107. Winston R.B. Upgrade to MODFLOW-GUI: Addition of MODPATH,, ZONEBDGT, and additional MODFLOW packages to the U.S. Geological Survey MODFLOW-96 Graphical-User Interface. U.S. Geological Survey Open-File Report 99-184 Reston, Virginia 1999, 63 p.

108. Winston R.B., Programs for Simplifying the Analysis of Geographic Information in U.S. Geological Survey, Ground-Water Models, U.S. Geological Survey Open-File Report 01-392, Reston, Virginia 2001,67 p.

109. Wolf J., Hack-ten Broecke M.J.D., Rotter R., Simulation of nitroen leaching in sandy soils in The Netherlands with the ANIMO model and the integrated modelling system STONE. Agriculture, Ecosystems and Environment 105 (2005), pp. 523-540

110. Yakirevich A., Borisov V., Sorek S. A quasi three-dimensional model for flow and transport in unsaturated and saturated zones: 1. Implementation of the quasi two-dimensional case. - Advances in Water Resources, Vol. 2, No. 8, 1998, pp. 679-689.

Информация о работе

Ромко, Алексей Викторович
кандидата технических наук
Москва, 2007
ВАК 06.01.02

Диссертация

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

Автореферат

Похожие работы