close

Вход

Забыли?

вход по аккаунту

?

ПРОБЛЕМА МОДЕЛИРОВАНИЯ

код для вставкиСкачать
ПРОБЛЕМА МОДЕЛИРОВАНИЯ
ПРОДУКТИВНОСТИ
АГРОЭКОСИСТЕМ: ОТ ЭВОЛЮЦИИ
ИДЕЙ К ЭВОЛЮЦИИ АЛГОРИТМОВ
Р.А. Полуэктов, А.Г. Топаж
(Лаборатория математического
моделирования агроэкосистем АФИ)
Фотосинтез (история)
М.Монси и Т.Саэки ˝Ueber den Lichtfactor in den
Pflanzengeselschaften und seine Bedeutung fuer die
Stoffproduktion (1953). Впервые поставлена задача
моделирования процесса фотосинтеза в условиях с.-х.
поля.
Дальнейшие работы : Гаастра (Gaastra, 1959), М.И.
Будыко (Будыко, 1964), Шартье (Chartier, 1970), и Лайск
(Лайск, 1977).
Теория фотосинтеза на уровне листа:
В.Л. Калер (Калер и др., 1976), Х. Тооминг, (Тооминг,
1977), де Вит. (Wit, 1978),
Фотосинтез (реализация)
Р.А.Полуэктов, 1991 г.:
три реакции
(1) возбуждения акцептора, (2) окислительного
карбоксилирования и (3) регенерации акцептора,
проходящую в присутствии дополнительного фермента
X Q Ф X ,
*
k2
P
X CO 2 Ф X
k3
X P B X
B
.
k1
*
Соотношения, определяющие
интенсивность фотосинтеза
Скорость фотосинтеза единичной
поверхностью листа (мг CO2/(м2•с))
где
QФ-поглощенная ФАР,
Cw – концентрация CO2 в хлоропласте
Фmm- макс. скор. фотосинтеза
Rd –интенсивность фотодыхания
rx – сопротивление карбоксилирования
α - наклон световой кривой ф.-с.
Интенсивность газообмена:
где
Ca - концентрация CO2 в атмосфере
rCΞ - суммарное сопротивление
rx
1
1
n C w Q mm
n Ca Cw
rC 1
Rd
Зависимость фотосинтеза листа от
интенсивности ФАР и концентрации CO2
0 .8
2
Б р у тто -ф о то си н те з, м г/(м с)
0 .7
0 .6
0 .5
0 .4
0 .3
0 .2
0 .1
13 .8
12
10 .2
Р и с.2.7.
3
0
35
70
2
4 .8
105
Ф А Р , В т /м
140
175
6 .6
210
24 5
8 .4
2 80
315
3 50
0
К о н ц е н т ра ц и я
C O 2 ,10
-7
г/см
3
Влагоперенос в системе почва – растение –
воздух (история)
Уравнение Ричардса (Ричардс, 1931)
( x, t )
t
w
P ( x, t )
1 fˆ ( x , t ),
k (P)
x x
С функцией стока, начальными и граничными условиями (Нерпин и др.,
1976; Diepen, 1988)
( x , 0 ) w ( x ),
Pr,
w
q x ( 0 , t ) q up t Es ,
если
Pr 0 ,
если
Pr 0 ,
NR
Er i 1
i
( p si PR )
P ( H s , t ) PH
Транспирация и физическое испарение
(реализация)
Потенциальная эвапотранспирация (Penman, 1948; Monteith, 1957)
Моделирование реальных значений транспирации и физического
испарения (Полуэктов и др., 1992)
E p f ( R n , Ta , q a , u , PL )
NR
Er i ( p si PR )
i1
Ep » Er
PL » PR
П ог ло щ е ни е в од ы к ор ня м и (E r ) и тран с пир а ция (E p ), м м /сут
Баланс воды в растении
5
P
4
*
L
Er
3
2
1
Ep
0
-1 ,5 0
-1 ,2 5
-1 ,0 0
-0 ,7 5
-0 ,5 0
П оте н ц и ал в од ы в ли сте (P
Р и с. 6.9 .
L
), М П а
-0 ,2 5
0 ,0 0
Динамика влагозапаса в почве
(Яровой ячмень, Белогорка, 2000 г.)
а
В лаг оз ап ас, см ( 1) ; оса дки , м м ( 2 )
35
1
30
25
20
2
15
10
5
0
10
20
30
40
50
60
70
80
90
10 0
110
Рост растений (история)
Детерминированные ростовые функции – distribution keys
(Росс,1971; Сиротенко, 1975; Penning de Vries et al., 1982)
y
1 ,0
С теб ли
0 ,8
0 ,6
0 ,4
К оло с
Л и сть я
0 ,2
x
0 ,0
0 ,0
0 ,5
1 ,0
1 ,5
2 ,0
Адаптивная ростовая функция (Thornley, 1972)
Рост растений - shoot : root relation
(реализация, Полуэктов, Топаж, 2005)
Распределение ассимилятов soot – root:
W r c rs W s (1 c rs ) 0 c rs i 1
Потребность в азоте:
N s ( k ) W s ( k ) KS c N shoot ,
N r ( k ) W r ( k ) KR c N root ,
Доступный азот
N av V N ( k 1) ( k ) ,
Баланс азота
N s ( k ) N r ( k ) N av ( k )
( k ) S root W r ( k )
Баланс азота в растении
2
3
N itrogen absorp tion b e cro p, kg h a
-1
d
-1
(реализация)
2
1
1
0
0
0 ,2
0,4
0,6
0,8
P art of assim ilates allocated to roots
1
Влияние азотных удобрений на
отношение root:shoot
0,35
О тн о ш ен и е roo t:shoo t
0,3
0,25
0,2
1
0,15
0,1
2
Зависимость отношения сухой
биомассы корней к сухой
надземной биомассе от дозы
азотных удобрений в посеве
ячменя. Представлены средние
арифметические из 7
компьютерных экспериментов
при вводе в модель погодных
условий за 1983-1989 гг.
Δ – в фазе цветения,
О – в фазе полной спелости.
0,05
0
0
20
40
60
Д оза азотны х удобрени й, кг/га
80
1 00
Модель с имитаций влияния доз азотных
удобрений
V N ( i ) V max
Поглощение азота в i- слое
почвы:
N con ( i )
N con ( i ) K N
i NR
VN Суммарное поглощение:
Результат
компьютерного
эксперимента:
Аппроксимирующая
зависимость
V
N
(i )
i 1
N=0
VN~0.2 Vmax
N=45 кг/га
VN~0.5 Vmax
N=90 кг/га
VN~0.7 Vmax
V N V max ( c 0 c1 N c 2 N )
2
Идентификация модели
cor * 0
1. общая продуктивностькоэффициент cor:
Озимая рожь
Яровой ячмень
Картофель
0.9
0.85
1.4
2. Влияние азотных удобренийкоэффициенты
CO, C1, C2
Расчетные и фактические урожаи озимой ржи
(МОС, 1982-2005 гг.)
Годы
вегетации
82
83
84
85
86
87
88
89
01
05
Опыт без удобрений
Расчет
Эксп.
38.8
30.7
21.3
21.7
20.8
20.8
19.8
19.8
12.6
12.7
18.0
18.0
27.0
26.9
25.3
25.2
29.0
33.4
35.5
58.0-
Средняя доза
Расчет
Эксп.
45.8
45.6
33.8
33.4
32.0
32.0
31.5
31.5
18.6
18.6
22.7
22.7
42.2
42.1
38.5
38.3
38.0
38.0
41.2
61.2-
Высокая доза
Расчет
Эксп.
47.8
47.6
38.9
38.8
32.8
32.8
36.5
36.520.5
20.5
26.0
25.9
48.2
48.0
42.5
42.5
44.4
65.0
Расчетные и фактические урожаи озимой
ржи
(МОС, 1982-2005 гг.)
70
Р а сч ет н ы й у рож а й , ц /га
60
50
40
30
20
10
0
0
10
20
30
40
50
Ф а к ти ч еск и й у р ож а й , ц /га
60
70
Расчетные и фактические урожаи ячменя
(МОС, 1983-2004 гг.)
Годы
вегетации
83
84
85
86
87
88
89
00
01
02
04
Опыт без удобрений
Расчет
Эксп.
20.1
24.1
18.8
16.6
16.7
22.6
13.7
17.4
19.7
20.2
14.1
15.8
22.9
16.0
21.0
22.4
22.3
21.0
14.3
15.4
15,7
15.1±1.4
Средняя доза
Расчет
Эксп.
29.4
33.1
27.4
27.5
23.4
27.6
21.7
24.8
27.1
25.7
22.0
21.8
33.6
25.2
28.0
32.4
31.2
27.0
23.6
21.6
21,7
21.4±2.4
Высокая доза
Расчет
Эксп.
37.8
35.0
30.2
30.4
25.6
32.2
24.8
27.7
29.5
27.2
25.1
22.8
36.5
26.6
30.1
36.0
33.7
31.2
28.1
22.6
23,5
22.1±4.4
Расчетные и фактические урожаи ячменя
(МОС, 1983-2004 гг.)
40
Р а сч ет н ы й ур ож а й , ц /га
35
30
25
20
15
10
5
0
0
5
10
15
20
25
30
Ф ак ти чески й ур ож а й, ц/га
35
40
Динамика роста органов растения
Биомасса органов растений, кг/га
(яровой ячмень)
а
7000
3
6000
5000
4000
4
3000
2
2000
1000
1
0
0
20
40
60
80
100
МОС, 2000 г. а – опыт с однократной дозой удобрений,
б – опыт без удобрений. 1, 2, 4 – биомасса листьев,
стеблей и колоса соответственно; 3 –надземная биомасса;
– эксп. данные.
б
Р и с . 7 .7 а
5000
3
4000
3000
4
2000
2
1000
1
0
0
20
40
60
80
Д н и , н ачи н ая с сева
100
Фактические и расчетные сроки
наступления фенофаз
аа
Э ксп ерим ен тальн ы е даты , сут
250
200
150
Яровая пшеница (Саратов, 1972–1987 гг.) (а)
и яровой ячмень (МОС, 1983–1989, 2001–2002 гг. (б).
На осях отложены номера дней по календарю.
100
100
150
200
250
б
260
220
180
140
100
100
140
180
220
Р асчетн ы е даты , сут.
260
Прогноз сроков начала полевых работ
Н ачало
п
о
л
О тн оси тельная влаж ность п очвы
р
е
а
б
в
ы
о
х
т
30
2 7 ,5
25
В Р К - спелость п очвы
2 2 ,5
20
1 7 ,5
15
1 2 ,5
10
0
50
100
150
Н ом ер дн я с 1-го октбря
200
250
Прибавка урожая от различных доз
азотных удобрений
Я ров ой я чм е нь, прибав ка урож ая
П ри б а вк а , ц /га
12
10
8
N30
6
N60
N90
4
2
0
1983 2000 1987 1985 2001 1984 2002 1988 1989
Го д
Фактические и прогнозируемые урожаи
(МОС, озимая рожь)
Ра сч етн ы й и п р ог но зи р у еы й у р о ж ай , ц /га
70
60
50
40
30
20
10
10
20
30
40
50
Ф акт ич ески й урож ай , ц /г а
60
70
Динамический прогноз урожая
Ф а к т и ч е с к и й ур о ж ай
П рогнозируемы й урож ай, ц/га
35
30
25
У б о р ка
20
15
0 5 .0 7
1 5 .0 7
2 5 .0 7
0 5 .0 8
1 ,6
1 ,7
1 0 .0 8 1 3 .0 8
10
5
0
1
1 ,1
1 ,2
1 ,3
1 ,4
1 ,5
Ф и з и о л о ги ч е с к о е в р е м я , б /р
1 ,8
1 ,9
2
Прогноз дозы азотных подкормок
Пр о г н о з д о зы аз о т н ы х п о д ко р мо к п р и о сн о в н о м уд о б р ен и и N45P 90K50
45
40-45
40
35-40
35
30-35
30
25-30
25
20-25
Д оза а зот но й подк орм к и, к г/ га
15-20
20
10-15
15
5-10
10
0-5
5
45
0
40
Же ла е м ый
5
35
10
15
20
30
25
Биом а с с а на м ом е нт
к о лош е ния , ц/ га
30
у рож а й, ц/ га
Диаграмма классов модели
(Rational Rose)
Система имитационного моделирования
AGROTOOL, v.4
Модель - Rational Rose,
объектный Паскаль;
Оперативная база данных (ОБД)
- Excel;
Стационарная база данных
(СБД) - СУБД Access;
Интерфейс пользователя Delphi
Верификация
1. Ленинградская область:
агрометеостанции – Белогорка, Волосово, Волхов, Сосново,
Николаевское, Тихвин;
культуры – яровой ячмень, яровая пшеница, озимая рожь, овес,
картофель, многолетние травы, 5-6 лет вегетации.
2. Саратовская область – яровая пшеница, 12 лет вегетации.
3. Краснодарский край – озимая пшеница, кукуруза, 6 лет вегетации.
4. Алтайский край – люцерна, 5 лет вегетации
5. Шекинский р-н Азербайджана – озимая пшеница, 6 лет вегетации.
6. Мюнхеберг, Германия – озимая пшеница, яровой ячмень, сахарная
свекла, 3 года вегетации.
7. Бадлаухштадт, Германия - озимая пшеница, яровой ячмень,
сахарная свекла, 6 лет вегетации.
8. Советский р-н Калининградской обл. – озимая пшеница,
многолетние травы – 5 лет вегетации.
Сотрудники лаборатории, принимавшие
участие в разработке модели
1.Пых Ю.А.
2.Вол И.А.
3.Финтушал С.М.
4.Заславский Б.Г.
5.Заславская Л.А.
6. Неусыпина Т.А.
7. Нагиев А.Т.
8. Чувашина Н.В.
9.Топаж А.Г.
10.Опарина И.В.
11.Захарова Е.Т.
12.Терлеев .В.В.
13.Бакалнко Б.И.
14.Глядченкова .Н.А.
Благодарю за
внимание
Документ
Категория
Презентации
Просмотров
14
Размер файла
558 Кб
Теги
1/--страниц
Пожаловаться на содержимое документа