Совершенствование методов расчета капельного орошения плодовых культур в условиях Египта

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

ФГОУ ВПО

Московский государственный университет природообустройства

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

УДК 631.6

Совершенствование методов расчета капельного орошения плодовых культур в условиях Египта

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

АВТОРЕФЕРАТ

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

махмуд мохамед али абдель азим

• Москва 2010
• Работа выполнена на кафедре мелиорации, рекультивации и охраны земель ФГОУ ВПО «Московского государственного университета природообустройства».
• Научный руководитель: доктор технических наук, профессор
• Голованов Александр Иванович
• Федеральное государственное образовательное учреждение высшего профессионального образования «Московский государственный университет природообустройства» (ФГОУ ВПО МГУП)
• Официальные оппоненты:
• доктор сельскохозяйственных наук, профессор
• Шуравилин Анатолий Васильевич
• (ФГОУ ВПО “Российский университет дружбы народов”)
• кандидат технических наук, доцент
• Рыбкин Владимир Николаевич
• (ФГОУ ВПО “ Московский государственный университет природообустройства ”)
• Ведущая организация: Российский государственный аграрный университет - МСХА им. К.А. Тимирязева
• Защита состоится «____»_________ 2010 года в 15 часов на заседании диссертационного совета Д 220.045.01 ФГОУ ВПО «Московский государственный университет природообустройства» по адресу: 127550, г. Москва, ул. Прянишникова, 19 (1 учебный корпус, аудитория. 201).
• Тел./факс: 8(495) 976-10-46.
• С диссертацией можно ознакомиться в библиотеке ФГОУ ВПО «Московский государственный университет природообустройства».
• Автореферат диссертации размещен «04» марта 2010 года на официальном сайте ФГОУ ВПО МГУП по адресу: http://www.msuee.ru/html/19_1.html
• Автореферат разослан «___»________________2010 г.
• Учёный секретарь Диссертационного совет Сурикова Т.И.
• ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

эвапотранспирация капельный орошение почва

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

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

Задачи работы:

· оценить применимость различных способов оценки потенциальной эвапотранспирации (суммарного испарения);

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

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

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

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

Научная новизна работы:

· выявлены наиболее подходящие к условиям Египта способы оценки потенциального суммарного испарения в зависимости от местоположения метеопунктов относительно побережья Средиземного моря; для ряда метеостанций применима формула Н.И. Иванова (приморская зона), для пустынной зоны - метод Блейни и Криддла и стандартный метод ФАО РАД;

· на основании разработок кафедры мелиорации и рекультивации земель МГУП предложен и верифицирован вариант модели двумерного влагопереноса при специфическом (траншейном) капельном увлажнении;

· исследовано влияние предполивной влажности на размер оросительной нормы и продуктивность плодовых деревьев; оптимальным уровнем влажности является значение 0,62…0,65 ППВ; во вневегетационный период для поддержания жизнедеятельности растений редкие поливы следует назначать при снижении предполивной влажности до 0,5 ППВ;

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

На защиту выносятся:

· результаты моделирования влагопереноса при капельном орошении в условиях Египта; размеры оросительных норм при траншейном капельном орошении;

· результаты исследования рациональной предполивной влажности в период вегетации и во вневегетационный период для плодовых культур;

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

Практическое значение и реализация работы:

· результаты исследований позволяют принимать научно-обоснованные решения при оценке рациональности и продуктивности орошения плодовых культур в условиях Египта;

· получать разумные урожаи плодов и фруктов при экономном расходовании поливной воды;

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

Достоверность результатов основана на: детальном учете природных условий орошаемых земель Египта; использовании современных разработок в области влагопереноса при орошении, использовании практического опыта орошаемого садоводства, статистической обработке экспериментов.

Апробация результатов исследований и публикации

Основные методические положения и полученные результаты докладывались и обсуждались на научно-технических в г. Москве (МГУП-2008, 2009 гг.); на Международной конференции молодых учёных-мелиораторов 24 сентября 2009г. Коломне, (ВНИИ «Радуга»); the First International Conference on: Economists and Management of Water in Arab World and Africa 18-19 November, 2009 Assiut. Egypt; на международной научной конференции "Водное хозяйство - состояние и перспективы развития", 15-16 апреля 2010 года в Национальном университете водного хозяйства и природопользования, г. Ровно (Украина).

Публикации

По результатам исследований опубликовано 6 печатных работ, в том числе 3 работы в изданиях, рекомендованных ВАК РФ и 1 работа в трудах First International Conference on: Economists and Management of Water in Arab World and Africa.

Структура и объём диссертационной работы. Диссертационная работа состоит из введения, четырех глав, выводов и приложений, списка литературных источников из 103 наименований. Общий объём диссертации составляет 144 страниц компьютерного текста, 32 рисунка и 23 таблицы.

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

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

В основе исследований по ирригации находятся работы Аверьянова С.Ф., Айдарова И.П., Бездниной С.Я., Голованова А.И., Губера К.В.; Добрачева Ю.П., Иванова Н.Н., Кирейчевой Л.В., Ковды В.А., Костякова А.Н., Константинова А.Р., Маслова Б.С., Ольгаренко Г. В. Плюснина И.И., Рекса Л.М., Селянинова Г.Т., Шабанова В.В., Харченко С.И.; вопросы капельного орошения изложены в работах Алексашенко А.А., Алпатьева С.М., Болцату И.Н., Григорова М.С., Кузнецова Е.Н., Мекки М., Нестеровой Г.С. и Зона И.С., Кружилина И.П., Ромашенко М.И., Шумакова Б.Б., Brant A., Breasler E., Feddes R.A. и других ученых в области мелиорации сельскохозяйственных земель, растениеводства, климатологии, гидрогеологии, земледелия и почвоведения.

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

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

Во второй главе приводится расчет и моделирование капельного орошения. Основной задачей мелиорации А.Н. Костяков считал управление круговоротом воды и зольных питательных элементов в целях прогрессивного повышения плодородия почв (Костяков, 1960). Совместное управление геологическим и биологическим круговоротом воды и зольных питательных веществ может быть достигнуто при условии комплексной взаимосвязи мелиоративно-гидротехнических и агротехнических методов, составляющих единую систему мелиоративных мероприятий, отвечающих природным и хозяйственным условиям данного района. А.И. Голованов развил эти представления и применительно к землям сельскохозяйственного назначения сформулировал: «цель мелиорации заключается в расширенном воспроизводстве плодородия почвы, получении оптимального урожая определенных сельскохозяйственных культур при экономном расходовании всех ресурсов, недопущении или компенсации ущерба природным системам и другим землепользователям» (Основы природообустройства, 2001). Однако цели мелиорации земель могут быть достигнуты только при выполнении набора требований к управляемым факторам почвообразования, роста растений и воздействия на окружающую среду, которые должна обеспечивать система мелиоративных мероприятий. Этот набор требований А.И. Голованов и И.П. Айдаров назвали мелиоративным режимом (Айдаров, Голованов, 1986).

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

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

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

Табл. 1. Исходные данные для моделирования

Культура	Площадь дерева, м2	Площадь кроны, м2	Kб / Kзат	Оптимальная влажность, доли ППВ	Предполивная влажность, доли ППВ	Период вегетации
Персик	24	18	0,61/ 0,4	0,82	0,6 … 0,7	1.01-31.05
Яблоко	24	18	0,63 / 0,4	0,82	0,6 … 0,7	1.02-30.06
Цитрус	21	15,7	0,58 / 0,4	0,82	0,6 … 0,7	1.01-30.07
Виноград	10,5	7,8	0,53 / 0,4	0,82	0,6 … 0,7	1.01-30.06
Манго	15	11.3	0,62/0,4	0,82	0,6 … 0,7	1.02-15.06

Обозначения в таблице: Kб - биологический коэффициент Kзат-коэффициент затененности почвы растительным покровом.

Модель, использованная для расчетов, представляет собой уравнение двумерного потока влаги, реализованное в виде конечно-разностной численной схемы. Исследуемая толща разбивалась на элементарные слои переменной толщины, от 0,1 м вблизи поверхности до 1 м вблизи водоупора, . Для учета разной площади, обслуживаемой одной капельницей, рассматриваемый пласт шириной (нормально к плоскости чертежа) B разбивался вертикальными плоскостями для образования столбцов и расчетных блоков. Ширина этих блоков (по длине катены) bi принималась различной в зависимости от ее длины, при этом b0 = bNy = 0.

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

где

- напор, м, на расчетный момент времени ; при отсчете напоров от поверхности земли

;

- напор, м, эквивалентный каркасно-капиллярному давлению в зоне неполного насыщения () и эквивалентный гидростатическому давлению в зоне полного насыщения;

- коэффициент влагоемкости, м3в/м4:

;

где , - объемная влажность почвы, м3в/м3 (м3в - кубический метр почвенной влаги). При полном влагонасыщении Сw = 0. Связь между каркасно-капиллярным потенциалом и влажностью почвы принята в виде:

;

где m - пористость, м3/м3; - максимальная гигроскопичность, м3/м3; - максимальная высота капиллярного поднятия, м; м и n - коэффициенты, зависящие от механического состава и структуры почвы, для суглинистых почв принято , показатель степени n = 1.

Коэффициент влагоемкости при n = 1 равен:

.

- вертикальное сопротивление потоку влаги между центрами и блоков, сут;

;

- горизонтальное сопротивление потоку влаги между центрами и блоков, сут;

;

K(щ) - коэффициент влагопроводности м3в/м2/сут, зависящий от коэффициента фильтрации Kф; пористости почвы m; объемной влажности почвы ,; для определения коэффициента влагопроводности в зависимости от влажности почвы пользуются формулой С.Ф. Аверьянова, предложенной им в 1947 году:

,

где щ* - максимальная молекулярная влагоемкость или влажность разрыва капилляров, для более широкого диапазона влажности эту зависимость можно заменить на

.

Расходование влаги на испарение принято зависящим от погодных условий и от влажности почвы, оно разделялось на испарение с поверхности почвы, которое учитывалось как граничное условие, и на транспирацию, последняя распределялась по корнеобитаемому слою пропорционально влажности почвы и плотности корней и входила в уравнение в виде интенсивности влагоотбора корнями растений из единичного объема почвы ei,j, м3в/м3/сут. С этой целью для каждой декады теплого периода по известным средней температуре воздуха T, оС и относительной влажности воздуха a, % подсчитывалось потенциальное (при оптимальной влагообеспеченности) суммарное испарение (эвапотранспирация) Epot по формуле Н.Н. Иванова:

, мм/сут

где Kб - биологический коэффициент, учитывающий особенности конкретной культуры.

Потенциальная эвапотранспирация разделялась на потенциальное испарение с поверхности почвы и потенциальную транспирацию пропорционально затененности почвы растительным покровом fр, которая изменялась по декадам: = (1 - fр) и = fр. Эти потенциальные величины испарения редуцировались на каждом временном шаге:

; ; ;

при влажности поверхностного 2…5 см слоя почвы п >0,8p ' = 1; эти зависимости согласуются, например, с исследованиями А.И. Будаговского. Фактическая транспирация редуцируется в зависимости от неоптимальности средней влажности корнеобитаемого слоя почвы:

, где , ;

w - коэффициент, учитывающий уменьшение транспирации при отклонении влажности почвы от оптимальной, вид этой зависимости соответствует исследованиям А.Р. Константинова (1968); k - средняя влажность корнеобитаемого слоя почвы, переменная во времени; - то же, оптимальная в данную декаду; ВЗ - влажность завядания.

Скорректированная величина транспирации распределялась по глубине каждого столбца в заданном корнеобитаемом слое пропорционально влажности почвы и массы корней в виде интенсивности влагоотбора корнями растений из единичного объема почвы ei,j, м3в/м3/сут.

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

Uj = | ; ; ; . . . | при этом j= 0,1,2,3...Nx.

С помощью этого вектора система уравнений запишется в матричном виде:

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

Левое граничное условие, т.е. отсутствие потока в центре рассматриваемого пласта учитывается особыми правилами вычисления элементов этой матрицы D11 и D12; D11 = -D12; . Аналогичное правое граничное условие учитывается при вычислении последних элементов этой матрицы ; . В случае, если в каком-то блоке i,j имеется источник или сток (дрена или канал), они учитываются при вычислении соответствующих элементов матрицы .

Вектор Fj объединяет все свободные члены.

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

Верхнее граничное условие учитывают, особым образом вычисляя первую матрицу прогоночных коэффициентов PP0 и первый вектор-столбец QQ0 при обратной прогонке. Так, если через верхнюю границу (х=0) нет потока влаги, то , тогда элементы диагонали матрицы PP0 равны , а остальные - нулевые. Все элементы вектора-столбца QQ0 равны нулю. При физическом испарении через поверхность почвы поток влаги равен:

.

Поэтому диагональные элементы матрицы PP0 равны , а остальные - нулевые. Элементы вектора-столбца QQ0 равны: . Если поверхность почвы увлажняется поливом, то в приведенных выражениях величина физического испарения заменяется интенсивностью водоподачи (м/сут), взятой с обратным знаком. Если поливается часть поверхности, то эту замену осуществляют для соответствующих номеров столбцов i.

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

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

Для обустройства модели использованы некоторые данные полевых исследований, которые проводились кафедрой мелиорации и рекультивации земель МГУП на ОПУ системы капельного орошения «Таврия» в базовом хозяйстве УКРНИИГиМ Красногвардейского района Крымской области. В табл.2 дана сводка значений водно-физических свойств для условий Египта, использованных в расчетах.

Табл.2. Водно-физические свойства для условий Египта

Станция	Бехера
Гор	Н	p	щм	Kф	Hk	Почвы по мех. составу
A	1	0,48	0,1	0,25	3	Суглинок
B1	1,5	0,45	0,09	0,1	2,8
B2	3	0,44	0,08	0,1	2,7
C	глубже	0,45	0,09	0,1	2,8
Станция	Гиза
A	1	0,38	0,08	1,5	1	Песок
B1	1,5	0,35	0,07	1,5	0,9
B2	3	0,34	0,06	1,5	0,8
C	глубже	0,35	0,07	1,5	0,9
Станция	Исмаила
A	1	0,45	0,09	0,8	1,5	Супесь
B1	1,5	0,42	0,08	0,5	1,3
B2	3	0,41	0,07	0,5	1,2
C	глубже	0,42	0,08	0,5	1,3
Станция	Фаем
A	1	0,45	0,09	0,8	1,5	Супесь
B1	1,5	0,42	0,08	0,5	1,3
B2	3	0,41	0,07	0,5	1,2
C	глубже	0,42	0,08	0,5	1,3
Станция	Вади
A	1	0,48	0,1	0,25	3	Суглинок
B1	1,5	0,45	0,09	0,1	2,8
B2	3	0,44	0,08	0,1	2,7
C	глубже	0,45	0,09	0,1	2,8

Обозначения в таблице: Н - толщина горизонта, м; p - пористость; щм - максимальная гигроскопичность; Hk - высота капиллярной каймы, м; Kф - коэффициент фильтрации, м/сут.

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

Суммарным водопотреблением (эвапотранспирацией) называют количество воды, расходуемое на транспирацию растениями и испарение почвой за расчетный период. Обозначим эту величину Е, а ее потенциальное значение или испаряемость, т.е. при неограниченном количестве почвенной влаги E0 ( в зарубежной литературе она часто обозначается как ETо); её измеряют в мм или в м3/га за период вегетации, месяц, декаду, неделю, сутки.

Е = Т + И;

где Т - транспирация воды растениями; И - испарение воды из почвы.

Величины Т и И зависят от метеорологических условий, вида и фазы развития растений, площади листьев, состояния почвы - разрыхленности, затененности листьями и очень сильно от увлажненности почвы. Соотношение между Т и И изменяется по мере развития растений и изменения условий испарения. За период вегетации для большинства сельскохозяйственных культур на испарение с поверхности почвы затрачивается около 30…40 % воды, на транспирацию растений - примерно 60…70 % от суммарного водопотребления. Послеполивное рыхление почвы (культивация) позволяет существенно сократить непроизводительные потери воды на испарение.

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

E = еE0 , е = 2,72•и/иopt ·Exp(-и/иopt);

где и - относительные доступные влагозапасы в корнеобитаемом слое, и = (щ - ВЗ)/(p -ВЗ); иopt - то же оптимальные для конкретной культуры; p - пористость почвы в расчетном слое; ВЗ - влажность завядания, примерно равная 1,3…1.5 максимальной гигроскопичности почвы.

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

Была выполнена оценка применимости некоторых методов оценки потенциального суммарного испарения в условиях Египта. Цель этого исследования состояла в том, чтобы оценить точность восьми методов, сравнивая их со стандартным методом ФАО 56 - ПМ для двенадцати метеостанций Египта при различных климатических условиях. Рассмотрены метеостанции Матрох, Бехера, Фаюм, Вади, Харга, Асван, Кафр-саад, Ага, Исмаила, Гиза, Миня, Манфалот. Были использованы все имеющиеся метеоданные за 10 лет.

Методы были сгруппированы в 3 группы:

А) стандартный метод ФАО 56- Пенман-Монтейс (ПМ);ФАО Блэйни-Кридла (БК), Иванова (ИВ), Тюрка, Дженсен-Хэйса (ДХ) и Харгривса (ХГ);

Б) радиационные методы ФАО Рад, Пристли Тэйлор (ПТ);

В) метод с использованием водных испарителей 3(И-Пан).

Значения испарения, полученные различными методами были сравнены со стандартным методом 56 ФАО-ПМ следующим образом:

1) Была вычислена относительная ошибка PE делением значений испаряемости различными методами (метод) на значения стандартного метода ФАО 56 - ПМ (станд):

.

2) Определена стандартная ошибка оценки SEE:

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

Для зоны 1 (западная пустыня), включающей метеостанции Матрох, Бохера, Фаюм, Вади, Харга и Асван, наиболее подходящим для оценки потенциального суммарного водопотребления оказался метод Блэйни-Кридла (БК).

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

Для зоны 3 (восточная пустыня), включающей метеостанции Гиза, Миня и Манфалот, наиболее подходящим для оценки потенциального суммарного водопотребления оказался метод ФАО Рад

Была выполнена оценка способа расчета оросительных норм при капельном орошении сада в условиях Египта.

Целью данного исследования было сравнение фактических и расчетных оросительных норм, т.е. проверка жизнеспособности модели. Сравнение оросительных норм приведено в принятых в Египте размерностях [м3/феддан] (1феддан=4200м2). Надо также иметь в виду, что оросительная норма относится только к увлажняемой площади (см. табл. 1).

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

В третьей главе приводятся результаты исследования влияния предполивной влажности на размер оросительной нормы сада в условиях Египта.

Для исследований использовалась математическая модель, описанная выше. Для всех 5 культур и всех 5 метеостанций были рассмотрены различные значения предполивной влажности в пределах 0,62…0,75 ППВ в период вегетации и 0,4..0,55 ППВ во вневегетационный период, в том числе и обеспечивающие 85 и 75 % максимальной относительной урожайности, для чего предварительно были построены графики зависимости относительной урожайности от предполивной влажности. В табл. 4 приведены значения предполивной влажности а период вегетации и вне его, принятые при моделировании.

Относительная урожайность определялась по модели продуктивности растений, разработанной проф. В.В. Шабановым.

Затем были подсчитаны размеры оросительных норм при тех же значениях предполивной влажности. Объединяя эти зависимости, получили значения оросительных норм при разных уровнях урожайности: 85 и 75 % максимальной. При этом использовалась известная закономерность: заметное снижение оросительных норм при несильном отклонении влажности от оптимальной. Эта закономерность в качестве примера иллюстрируется рис. 3.

Табл. 3. Сравнение расчетных оросительных норм и относительных урожайносте (по модели) при оптимальной влажности с фактической

Станция	Культура	Оросительная норма по результатам моделирования, м3/феддан	Фактическая оросительная норма, м3/феддан	Урожайность по модели *	Фактическая урожайность **
Бехера	персик	5198	5600	0,93	0,73
	яблоко	5651	6000	0,94	0,75
	цитрус	5765	6200	0,93	0,76
	виноград	5036	5530	0,93	0,79
	манго	5451	5860	0,95	0,77
Гиза	персик	7299	6900	0,94	0,73
	яблоко	7645	7130	0,95	0,75
	цитрус	8261	7940	0,94	0,76
	виноград	7273	6970	0,94	0,79
	манго	7280	7000	0,96	0,77
фаем	персик	5944	6000	0,94	0,73
	яблоко	6587	6400	0,95	0,75
	цитрус	7027	6840	0,93	0,76
	виноград	6187	6300	0,92	0,79
	манго	6304	6450	0,94	0,77
вади	персик	5198	5500	0,93	0,73
	яблоко	5651	5980	0,94	0,75
	цитрус	5765	6190	0,93	0,76
	виноград	5038	5400	0,93	0,79
	манго	5451	5700	0,95	0,77
Исмаила	персик	5169	5600	0,94	0,73
	яблоко	5840	6000	0,95	0,75
	цитрус	6264	6200	0,94	0,76
	виноград	5532	5530	0,94	0,79
	манго	5498	5860	0,95	0,77

Примечания : * при оптимальной влажности.

** при фактической влажности.



Pис.1. Сравнение расчетной (по модели) оросительной нормы с фактической



Pис. 2. Сравнение расчетной (по модели) относительной урожайности с фактической (среднемноголетние данные)

Табл. 4 Исследованные варианты предполивной влажности

Вариант	Предполивная влажность в период вегетации	Предполивная влажность вне периода вегетации
1	0,70	0,56
2	0,65	0,52
3	0,62	0,49
4	0,60	0,48
5	0,65	0,49
6	0,62	0,46
7	0,62	0,40
8	0,62	0,65
9	0,60	0,65
10	0,60	0,62
11	0,60	0,60
12	0,58	0,58
Оросительная норма		Относительная урожайность
Предполивная влажность	Предполивная влажность

Рис. 3. Зависимость оросительных норм (слева) и относительной урожайности (справа) от предполивной влажности для Бехера яблоко.

Оросительная норма, м3/ф
Варианты предполивной влажности

Рис. 4 . Влияние предполивной влажности на размер оросительной нормы

Относительная урожайность
Варианты предполивной влажности

Рис. 5. Влияние предполивной влажности на относительный урожай

Исследования показали, что оптимальным уровнем влажности является значение 0,62…0,65 ППВ, которое обеспечивает получение около 85 % урожайности и сокращение оросительной нормы примерно на 40%, что и имеет место на практике; во вневегетационный период для поддержания жизнедеятельности растений редкие поливы следует назначать при снижении предполивной влажности до 0,5 ППВ.

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

где kщ- коэффициент влагопроводности, м/сут.; - коэффициент анизотропии грунта, т.е. отношение kфгор/kфверт; kф - коэффициент фильтрации увлажняемого слоя почвы, обычно определяемый в полевых условиях методом налива, то есть при вертикальной фильтрации, в связи с этим необходимо учитывать возможную анизотропию почвы, то есть неодинаковость значений коэффициента фильтрации в вертикальном и горизонтальном направлениях, что может быть существенным для почв, образовавшихся на лессовидных и аллювиальных породах, поэтому в формулу вводится коэффициент анизотропии , для лессовидных грунтов он может быть равен 0,3...0,6, а для слоистых аллювиальных он больше единицы; - расчетная объемная влажность увлажняемого слоя, зависящая от предполивной влажности щпп и поливной нормы m, мм

где ППВ - предельная полевая влагоемкость в долях от пористости почвы p; - максимальная гигроскопичность; - увлажняемый слой почвы в м.

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

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

(1)

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

,

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

, м-1

- высота капиллярного поднятия, м. (см табл. 5)

Расчетные значения водно-физических свойств почв для каждой метеостанции приведены в табл. 5.

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

Станция				, м	, м/сут	, м		, м/сут	, м-1	, сут
Фаем	0,45	0,09	0,65	1,5	0,80	0,50	0,21	0,005	0,143	7,20
Гиза	0,38	0,08	0,65	1	1,20	0,50	0,18	0,004	0,095	6,10
Бехера	0,48	0,10	0,65	3	0,25	0,50	0,22	0,001	0,039	13,7
Исмаила	0,45	0,09	0,65	1,5	0,80	0,50	0,21	0,005	0,143	7,20
Вади	0,48	0,10	0,65	3	0,25	0,50	0,22	0,001	0,039	13,7

Табл. 6. Зависимость коэффициента Щ в формуле (1) от поливной нормы m и водно-физических свойств почв применительно к метеостанциям

m, мм	Фаем	Гиза	Бехера	Исмаила	Вади
30	0,75	0,37	0,73	0,75	0,74
40	0,70	0,33	0,66	0,70	0,66
50	0,56	0,29	0,60	0,56	0,58
60	0,51	0,28	0,55	0,51	0,53

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

Табл. 7. Средняя продолжительность цикла полив-межполив в зависимости от размера поливной нормы

m, мм	Время межу поливами, сут
30	3,5
40	4,7
50	5,8
60	6,0

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

Рис. 6.а- динамика влажности в корнеобитаемой зоне в период вегетации и вне ее. и б- боковой отток нарастающим итогом

На рис. 8 в качестве примера показаны результаты проверки формулы (1) результатами моделирования на примере метеостанции Фаем (культура яблоко) для характерных по погодным условиям лет. Слева приведены кривые мгновенного (мм/сут.) и суммарного оттока (мм) за средний цикл полив-межполив и суммарный отток, рассчитанный по формуле (1) для разных поливных норм. Этот рисунок показывает хорошее совпадение рассчитанных значений оттока с данными моделирования. Справа показаны корреляционные зависимости между динамикой оттока по формуле (1) и по данным моделирования внутри периода полив-межполив. Угловые коэффициенты уравнений связи близки к единице, а коэффициент корреляции R2 практически равен единице, что также подтверждает правильность формулы (1).

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

Метеостанция	Культура	Отток за вегетацию	Отток за вне-вегетацию	Отток за год	Отток по формуле (1)
Бехера	яблоко	275	-75	200	260
	цитрус	344	-67	277	295
	манго	274	-75	199	236
	персики	252	-78	174	206
	виноград	304	-72	232	242
Фаем	яблоко	454	-510	-56	447
	цитрус	410	-451	-41	413
	манго	461	-520	-59	434
	персики	447	-493	-46	454
	виноград	459	-516	-57	475
Гиза	яблоко	290	-424	-134	305
	цитрус	267	-419	-152	268
	манго	291	-418	-127	297
	персики	295	-416	-121	305
	виноград	286	-419	-133	288
Исмаила	яблоко	420	-475	-55	447
	цитрус	359	-437	-78	379
	манго	432	-502	-70	434
	персики	416	-508	-92	396
	виноград	419	-494	-75	434
Вади	яблоко	274	-73	201	263
	цитрус	344	-66	278	299
	манго	275	-76	199	239
	персики	252	-78	174	197
	виноград	304	-72	232	245

Рис. 7. Рассчитанные отток по модели и по формуле (1) в мм для все культура и метеостанция.

Отток, мм
Относительное время
Отток по формуле (1)
Суммарный отток

Рис. 8. а) Сравнение мгновенного и суммарного оттока по данным моделирования и по формуле (1); б) Корреляционные зависимости между значениями оттока по данным моделирования и по формуле (1)

Общие выводы и предложения

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

2. Обобщены данные по погодным условиям 12 метеостанций, выполнен сравнительный статистический анализ результатов по 9 распространенным в мире способам расчетов потенциального суммарного испарения (эвапотранспирации), для условий Египта наиболее подходящими являются формула Н.И. Иванова (для приморской зоны), стандартный метод ФАО РАД и метод Блэйни-Кридла (для пустынной зоны).

3. Обобщены многолетние данные практического применения капельного орошения садовых культур (персика, яблока, винограда, цитруса, манго) в климатических условиях метеостанций Бехера, Гиза, Исмаила, Фаем, Вади; анализировались средние за 10 лет значения относительной урожайности и размеры оросительной нормы.

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

5. Исследовано влияние предполивной влажности на размер оросительной нормы и продуктивность плодовых деревьев; оптимальным уровнем влажности является значение 0,62…0,65 ППВ, которое обеспечивает получение около 85 % урожайности и сокращение оросительной нормы примерно на 40%, что и имеет место на практике; во вневегетационный период для поддержания жизнедеятельности растений редкие поливы следует назначать при снижении предполивной влажности до 0,5 ППВ.

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

Публикации по теме диссертации

1. Абдель-азим М.М. Сравнение некоторых методов оценки суммарного испарения в условиях Египта [Текст]/Абдель-азим М.М.// Материалы международной научно-практической конференции «Роль мелиорации и водного хозяйства в реализации национальных проектов»: М.: ФГОУ ВПО МГУП, 2008.-Ч.I.- C.7-11. ISBN 978-5-892 31-243-1

2. Абдель-азим М.М. Оценка способа расчета оросительных норм при капельном орошении сада в условиях Египта. [Текст]/ Голованов А.И., Абдель-азим М.М.// Мелиорация и водное хозяйство: научно- практический журнал. - М., 2009, двухмесячник,- ISSN 0235-2524. № 3 С.24-25

3. Абдель-азим М.М. Исследования влияния предполивной влажности на размер оросительной нормы сада в условиях Египта. [Текст]/Голованов А.И., Абдель-азим М.М.// Мелиорация и водное хозяйство: научно- практический журнал. - М., 2009, двухмесячник. - ISSN 0235-2524. № 6 С.48-49

4. Абдель-азим М.М. Evaluation the Irrigation Rate in Drip Irrigation and Identify Standards Irrigation in Terms of Egypt. [Текст] /Голованов А.И., Абдель-азим М.М.// The First International Conference «Economists and Management of Water in Arab World and Africa» 2009 Assiut University, Egypt. C. 223-229

5. Абдель-азим М.М. Влияние капельного орошения садовых культур на прилегающие земли в условиях Египта. [Текст]/ Голованов А.И., Абдель-азим М.М.// Мелиорация и водное хозяйство: научно- практический журнал. - М., 2010, двухмесячник. - ISSN 0235-2524. № 2 (в печати)

6. Абдель-азим М.М. Creation New Model of Drip Irrigation for Egyptian Conditions. [Текст]/ Голованов А.И., Абдель-азим М.М.// Материалы международной научно-практической конференции «Водное хозяйство - состояние и перспективы развития» г. Ровно, Национальный университет водного хозяй...