Автореферат (1151760), страница 3
Текст из файла (страница 3)
Во внутригодовом разрезе температура и относительная влажность имеют коэффициенты вариации допустимые для прогнозирования (соответственно 0,08-0,37 и 0,03-0,07), а внутригодовое распределение осадков неустойчиво (коэффициент вариации 0,5-0,6), имеет стохастический характер и не может достоверно прогнозироваться.
Для повышения точности прогнозных расчетов метеорологических параметров автором разработана эмпирическая модель прогнозирования гидрометеорологических факторов с использованием уравнения Фурье, которая учитывает цикличность изменения метеорологических параметров. Результат расчета по модели реализован в виде графика (рисунок 3), позволяющего проследить общую тенденцию изменения коэффициента природного увлажнения в течении длительного ряда метеонаблюдений.
С целью роста эффективности реализации модели предложена структура базы данных, обеспечивающая предварительную обработку метеорологических показателей для формирования эмпирической зависимости, характеризующей их цикличность.
Ку
годы
Коэффициент природного увлажнения
Рисунок 3 - Цикличность коэффициента природного увлажнения по уравнению Фурье
Для вычисления характера изменчивости данных достаточно рассчитать сумму первых членов ряда Фурье:
f(x)=a0/2+a1cos(πx/l)+b1sin(πx/l)+a2cos(2πx/l)+b2sin(2πx/l) (1)
где а0/2 – среднее значение метеорологических параметров, l – половина рассматриваемого периода , годы; x - текущий год;
/2 – начальный год.
Коэффициенты a1, a2, b1, b2 вычислены по формулам:
где интегралы заменяются на конечные суммы, рассчитываемые на отрезке [a, b] (dx = Δx = 1 год).
В качестве интегрального показателя изменчивости климата был выбран коэффициент природного увлажнения Ку, который наиболее объективно характеризует тепло-влагообеспеченность периодов в различные годы.
Ку=2,26/2+0,05cos(πx/l)+0,1sin(πx/l)+0,2cos(2πx/l)+0,03sin(2πx/l), (4)
где l - половина рассматриваемого периода, в годах; x – шаг, через который определяется коэффициент природного увлажнения(Ку).
где: Е – испаряемость (потенциальная эвапотранспирация) за период с температурой воздуха >50C, мм; Wa – активные влагозапасы в метровом слое почвы на начало расчетного периода, мм; P – осадки за период с температурой воздуха >50C, мм
Дальнейшее совершенствование расчетной методики прогнозирования гидрометеорологических условий и расширение ее функциональных возможностей осуществлялось за счет использования «сплайн-функций» для моделирования кривых периодических колебаний гидрометеорологических показателей и определения тепло-, влагообеспеченности вегетационного периода. В этом случае, зависимость, характеризующая изменчивость гидрометеорологических факторов с учетом цикличности колебаний за многолетний период имеет вид:
F(t)=a0+a1t+a2t2+aK1cos(tw1)+aK2 cos (tw2)+aK3 cos (tw3)+bK1 sin (tw1)+ +bK2 sin(tw2)+bK3 sin(tw3), (6)
где а0, а1, а2, аК1, аК2, аК3, bK1, bK2, bK3 - коэффициенты, определяющиеся методом наименьших квадратов; w1,2,3 - частота периода, определяющаяся при анализе всего диапазона Т1,2,3 [1,N] с шагом 0,25, где N – количество расчетных периодов.
Из полученных «сплайн-функций» для периодов Т1,2,3 выбираем тот Ку, среднеквадратичное отклонение от истинных значений которого минимально.
Ку=0,8868+0,013∙t-0,0001∙t2+0,094∙sin∙(t∙1,79)+0,102∙sin∙(t∙0,48)+
+0,187∙sin∙(t∙0,23)+0,07∙cos∙(t∙1,79) -0,003∙cos∙(t∙0,48)+0,087∙cos∙(t∙0,23), (7)
где t – порядковый номер члена ряда.
На рисунке 4 представлен коэффициент природного увлажнения, полученный на основе использования «сплайн-функции».
годы
Ку
Коэффициент природного увлажнения Ку
Рисунок 4 - Цикличность изменчивости коэффициента природного увлажнения, полученная на основе использования метода «сплайн-функции»
Полученное уравнение Ку предлагается использовать при долгосрочных прогнозах режимов орошения и прогнозировании урожайности сельскохозяйственных культур.
Сложность решения задачи оперативного планирования режимов орошения обусловливается большим количеством тесно взаимосвязанных вероятностных агроклиматических процессов, формирующих водный и пищевой режимы почвы. Основной информацией для выбора решений по нормам и срокам поливов служит: текущая информация об агроклиматических характеристиках, влагообмене в зоне аэрации, агрофизических и гидрологических константах почвы, фазах развития растений, динамике влажности почвы, получаемой на основе решения уравнения водного баланса.
В результате проведенных исследований сформирована информационная база данных, которая дает возможность учесть большинство протекающих на поле процессов при разработке комплексных математических моделей для расчета оросительных норм и оперативного планирования поливов (рисунок 5). Разработаны блок-схема (рисунок 6), алгоритмы решения задач по оперативному расчету режимов орошения и компьютерная программа.
Методика оперативного расчета режимов орошения позволяет выполнить расчет динамики запасов влаги в почве на каждом поливном участке; оперативно планировать режимы орошения для нескольких севооборотов с учетом их структуры; оценить возможность влияния на экологическую
Рисунок 5- Информационная база данных для оперативного планирования режимов орошения
Рисунок 6 – Блок-схема расчета эксплуатационных режимов орошения

обстановку величины оросительных норм и режима поливов; осуществить контроль за режимом поливов, прогнозировать оросительные нормы с учетом вероятностного характера гидрометеорологических факторов.
Модель прогнозирования водопотребления и планирования режима орошения позволяет провести расчеты режима поливов, обеспечивающего получение экономически оправданного уровня урожайности при минимизации потерь воды на инфильтрацию и сброс, максимальную замкнутость водного баланса.
Четвертая глава посвящена уточнению параметров (биоклиматических коэффициентов и ветровой функции) модели оперативного управления поливами на основе полевых исследований, изучению динамики водного баланса и пищевого режима посевов моркови, а также установлению закономерностей изменения суммарного испарения и урожайности при различном уровне влагообеспеченности и минерального питания, требующихся для формирования параметров модели оперативного управления поливами и прогнозирования продуктивности овощных культур в условиях ограниченности природных и материальных ресурсов.
Обработка и статистическая оценка полученных за период исследований показателей суммарного испарения E, испаряемости Ev, дефицита водопотребления dEv, осадков P проводились с использованием компьютерной программы «Расчет параметров режимов орошения сельскохозяйственных культур – «ROCK.xls». Определение изменчивости коэффициента природного увлажнения выполнялось по компьютерной программе «Raduga Irrigation» по результатам многолетнего ряда метеоданных за 58 лет для метеостанции Коломна (таблица 2).
Таблица 2 – Характеристика гидрометеорологических показателей по метеостанции « Коломна».
Год наблюдений | E | Ev | dEv | P | Kу | |||||
Значение, мм | Вероятность, % | Значение, мм | Вероятность, % | Значение, мм | Вероятность, % | Значение, мм | Обеспеченность, % | Значение, мм | Обеспеченность, % | |
2005 | 286 | 5 | 209 | 3 | 46,7 | 45 | 340 | 75 | 1,70 | 15 |
2006 | 322 | 20 | 231 | 15 | 32,0 | 25 | 254 | 95 | 1,67 | 20 |
2007 | 411 | 70 | 259 | 45 | 58,8 | 55 | 419 | 40 | 1,19 | 60 |
2008 | 312 | 15 | 228 | 12 | 8,4 | 8 | 620 | 1 | 1,85 | 5 |
Для поддержания расчетной влажности почвы в пределах от 0,8 до 1,0 НВ на опытных участках в различные по метеоусловиям годы потребовались следующие оросительные нормы: в 2005 – 62мм; в 2006 – 89мм; в 2007 – 190мм; в 2008 – 32мм.
Различный уровень минерального питания и водообеспеченности обусловили разницу в урожайности моркови по вариантам опыта (таблица 3). Как показали результаты экспериментальных исследований внесение минеральных удобрений дозой N100P100K185 положительно влияет на урожайность моркови, прирост составляет от 5,19% до 22,69 %. При снижении дозы минеральных удобрений до N80P80K150 прирост урожайности колеблется от 2,60% до 10,84%, а при увеличении дозы минеральных удобрений на 20% прирост урожайности колеблется от 12,99 % до 18,74%.
Максимальный прирост урожайности моркови был отмечен в 2007 году и составил по вариантам опыта 28 % - 58,31%, а в 2008 году отмечен наименьший прирост 1,3 % - 7,56%. Снижение урожайности в 2008 году обусловлено низкими температурами в течении вегетационного периода и малой активностью солнечной радиации.
Таблица 3 - Урожайность моркови по вариантам опыта, т/га.
Годы, процент обеспеченности Ку | Удобрения Орошение | Без удобрений | 0,8 NPK (N80P80K150) | NPK (N100P100K185) | 1,2 NPK (N120P120K222) |
2005 (15%) | Без орошения | 45,4 | 49,2 | 54,3 | 55,7 |
М (22-31мм) | 64,2 | 66,1 | 67,2 | 69,0 | |
0,8М | 63,7 | 66,4 | 66,1 | 67,1 | |
0,6М | 53,9 | 55,3 | 56,1 | 56,7 | |
2006 (20%) | Без орошения | 44,3 | 46,9 | 52,6 | 53,1 |
М (22-31мм) | 63,9 | 66,8 | 55,6 | 69,1 | |
0,8М | 61,8 | 65,3 | 66,8 | 68,2 | |
0,6М | 53,9 | 55,3 | 56,1 | 56,7 | |
2007 (60%) | Без орошения | 41,5 | 43,7 | 45,1 | 45,4 |
М (22-31мм) | 65,7 | 68,4 | 69,1 | 71,5 | |
0,8М | 64,2 | 66,1 | 67,2 | 69,0 | |
0,6М | 53,2 | 54,6 | 56,4 | 60,7 | |
2008 (5%) | Без орошения | 46,3 | 48,2 | 46,1 | 51,2 |
М (22-31мм) | 49,8 | 55,2 | 55,7 | 58,6 | |
0,8М | 47,3 | 52,1 | 51,2 | 55,3 | |
0,6М | 46,9 | 49,9 | 48,2 | 52,4 |
По результатам полевых опытов получена методом наименьших квадратов зависимость урожайности моркови от дозы вносимых удобрений, величины поливных норм и коэффициента природного увлажнения.