Введение

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

Исследования в области математического моделирования продукционного процесса растений интенсивно развивается, о чём свидетельствует множество работ в этой области. Будучи эффективным средством интеграции имеющихся знаний о процессе формирования урожая сельскохозяйственных культур, изучения динамики процесса в различных условиях внешней среды, модели продуктивности находят все более широкое применение для решения задач прогноза и управления [2,5,9,10].

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

Продуктивность севооборотов определяется продуктивностью культур, входящих в их структуру, которые по-разному реагируют на количество осадков и температуру в отдельные периоды вегетации. Исследованиями многих учёных доказано, что за ротацию севооборота продуктивность посевов достоверно зависит от природно-климатических условий и структуры посевных площадей [1,4,7].

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

Материалы и методы исследования

Исследования проводились в  длительных стационарных опытах по изучению чередования сельскохозяйственных культур в севооборотах прошедших более пяти ротаций. Стационар расположен в лесостепной зоне Западной Сибири на поля ФГБНУ «СибНИИСХ». Размещение делянок в опытах рендомизированное, размер делянок 0,275 га  (110×25 м) и 0,138 га (110×12,5 м), 4-х кратная повторность. Сорта полевых культур районированы для Западной Сибири: пшеница яровая – «Омская 36», ячмень яровой – «Омский 95», соя – «Золотистая», рапс – «Юбилейный», овес – «Иртыш 21»,  кукуруза – «Омка 135». Почва опытного участка − чернозем слабо выщелоченный среднегумусовый среднемощный тяжелосуглинистый. Содержание гумуса в пахотном горизонте 6,5 – 7 %. Система агротехнических мероприятий строилась с учетом рекомендаций СибНИИСХ для зоны лесостепи Западной Сибири [6]. Методы регрессионного, корреляционного анализа и математического моделирования [3]. Оценка экономических и агротехнических мероприятий и севооборотов проводилась согласно методическим рекомендациям сибирского отделения ВАСХНИЛ [ ].

 

Рис.1. Стационар в лесостепной зоне Западной Сибири на поля ФГБНУ «СибНИИСХ»

Результаты и их обсуждение

В целом введение в севооборот фуражных культур и бобовых (сои) с долей 15 – 25 % позволит увеличить выход КПЕ кормопротеиновых единиц на 0,3 – 0,7 тонны с 1 га севооборотной площади при сохранении высокого уровня урожайности зерна в сравнении с яровой пшеницей (таблица 1).

Таблица 1-  Продуктивность полевых севооборотов в зависимости от длины ротации и насыщенности полевыми культурами в 2010-2016гг.

В тоннах  

Севообороты

Урожайность

зерновых

Выход на 1 га пашни

зерна

кормовых единиц

переваримого протеина

КПЕ

Зернопаровые севообороты (Контроль)

Пар – пшеница – пшеница - овес

2,17

2,01

2,43

0,19

3,21

Занятый пар – пшеница –  пшеница – овес (без химизац.)

1,65

1,31

2,54

0,12

3,91

Занятый пар – пшеница –  пшеница  - овес (с химизац.)

2,32

1,95

3,87

0,31

6,78

Плодосменные севообороты.

Соя – пшеница – ячмень – овес

1,99

2,06

2,55

0,23

4,00

Занятый пар – пшеница – ячмень – овес (без химизац.)

1,87

1,4

2,87

0,28

5,42

Занятый пар – пшеница – ячмень - овес (с химизац.)

2,36

2,11

4,18

0,38

9,75

Бессменные посевы

Бессменная пшеница

1,10

1,10

1,92

0,15

2,24

Бессменный ячмень

1,37

1,37

3,36

0,25

5,49

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

Для повышения эффективного использования пашни при развитии животноводства возможно увеличение ротации севооборотов и введение в структуру площадей культур увеличивающих сбор кормовых единиц и переваримого протеина (соя и зернофуражные культуры). Введение этих культур в севообороты увеличивает продуктивность одного гектара пашни на 20 – 25 % [8].

В целом введение в севооборот зернобобовых и зернофуражных культур приводит к увеличению выхода зерна на 0,3 – 0,5 т/га, а КПЕ на 0,5 – 0,8 тонны с гектара севооборотной площади в сравнении с яровой пшеницей. Применение средств интенсификации способствовало увеличению продуктивности как зернопаровых, так и плодосменных севооборотов. В плодосменных севооборотах прибавка от их применения составила 0,15 – 0,22 т/га зерна, 2,5 – 3,8 т/га периваримого протеина.

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

В качестве переменных определены площади отдельных культур в севооборотах (хik, - площадь i-ой культуры в k-ом севообороте или звене), площади севооборотов (хk ), общая площадь пашни, объемы выхода продукции (Рj - объем j-ro вида продукции), издержек производства, чистого дохода (Yj - объем j-го стоимостного показателя; Qj - объем j-ro вида ресурса), энергозатрат, приращения энергии и другие.

Критерием оптимальности служит максимум производства продукции j-ro вида:

где N i - множество возделываемых культур; N 2 - множество севооборотов и их звеньев; b j i k - выход продукции j-ro вида с единицы площади i-ой культуры k-го севооборота или звена.

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


где M 4 - множество производственных ресурсов; с ik - доля i-ой сельскохозяйственной культуры в k-ом севообороте; d ik - выход валовой продукции с единицы площади i-ой культуры в k-м севообороте; a jik - затраты j-ro ресурса на единицу площади i-ой культуры в k-ом севообороте или звене; М5 - множество видов продукции; М6 - множество производственных затрат; М7 - множество валовой продукции; В j - гарантированный объем производства продукции j-ro вида.

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

Зернопаровые севообороты и плодосменное возделывание зерновых позволяют получить от 2,11 до 1,10 т/га товарного зерна (табл.1, 2).

Таблица 2 -  Структура севооборотов для получения товарного зерна на (N30P30  и солома) в среднем за 2010 – 2016 гг.

 

Севообороты

Бессменные посевы

Выход продовольственного зерна, т/га

Рентабельность по зерну, %

Энергетический коэффициент

Производительность агроэкосистемы, МДж/дн/ГДж

Баланс гумуса, т/га

Зернопаровые севообороты

Пар – пшеница – пшеница - овес

2,01

87,6

1,84

0,086

-0,78

Пар – пшеница –  пшеница  

1,95

94,9

1,72

0,098

- 0,21

Плодосменные севообороты

Соя – пшеница – ячмень – овес

2,06

91,0

2,01

0,123

+ 0,28

Рапс–пшеница – ячмень – соя - овес

2,11

94,9

2,14

0,127

+ 0,15

Бессменный посев

Бессменная пшеница

1,10

-12,6

1,34

0,040

- 0,18

При этом доля зерновых культур составляет 67-100% (таблица 3). На производство требуется затратить материально-денежных средств от 1094 до 1636,5 тыс.руб/га, совокупной энергии – 23,0-44,6 ГДж/га и эксплуатационной – 5,8-10,4 ГДж/га. Рентабельность производства по кормовым единицам с учетом побочной продукции варьирует от 17,0 до 43,9%, по зерну от 87,6 до 94,9 %. Только бессменное возделывание пшеницы оказывается нерентабельным. Коэффициент энергетической эффективности составляет 1,34-2,14 показатель производительности агроэкосистемы – 0,040 – 0,127 МДж/дн./ГДж. В зернопаровых севооборотах формирование урожая происходит за счёт естественного плодородия почвы, о чем свидетельствует отрицательный баланс гумуса, который достигает 0,18 – 0,78 т/га. При плодосменном чередовании урожай формировался на фоне воспроизводства органического вещества почвы, за счет увеличения растительных остатков от применения удобрений, баланс гумуса был положительным – 0,15 – 0,28 т/га.

Таблица 3 – Оценка структуры посевных площадей для получения продовольственного зерна на фоне N30P30  + солома (по результатам решения модели)

 

Показатели

Структура посевных площадей

Пар25%- зерновые 75%

Пар33% - зерновые 67%

Плодосмен (соя) 100%

Плодосмен (рапс) 100%

Бессменный посев 100%

Валовая энергия, Гдж/га

зерно

31,8

29,4

31,9

36,7

18,4

продукции

84,9

74,3

94,3

93,2

49,5

Выход кормовых единиц, т к.ед./га

продукции

2,43

2,26

2,55

3,01

1,92

Содержание протеина, т прот. ед./га

продукции

0,19

0,17

0,23

0,27

0,15

Затраты на производство, тыс.руб./га

1424,0

1279,4

1589,3

1636,5

1094,8

Чистый доход по к.ед., тыс.руб./га

588,7

578,3

696,2

709,7

-186,2

Рентабельность по к.ед., %

41,3

42,5

43,9

43,4

17,0

Затраты на зерно тыс.руб./га

1299,0

1154,6

1589,3

1636,5

1094,8

Чистый доход от зерна, тыс.руб./га

1138,5

1095,2

1446,3

1553,0

-137,9

Энергозатраты, Гдж/га

эксплуатационные

8,9

8,5

9,5

10,4

5,8

овеществлённые

27,4

26,2

28,6

34,2

17,2

совокупные

36,3

34,8

38,1

44,6

23,0

Приращение валовой энергии, Гдж/га

без учёта

86,8

78,8

96,2

105,3

58,5

с учётом гумуса

66,9

60,0

75,6

87,2

54,6

Энергопотенциал, Гдж/га

11228,6

9982,4

12563,8

13475

7493,2

Удобрения

N, ц.д.в./га

0,26

0,20

0,30

0,3

0,2

P, ц.д.в./га

0,26

0,20

0,30

0,3

0,2

Пестициды, кг д.в./га

1,0

0,67

2,00

2,00

0,67

Введение условия, не допускающего отрицательный баланс гумуса, исключает возможность использования на фоне N30P30  + солома зернопаровых севооборотов, включающих пар и одну-три пшеницы, так как они сопровождаются значительными потерями гумуса. Воспроизводство плодородия почв при максимальном выходе 2,11—2,06 т/га продовольственного зерна возможно при следующей структуре посевных площадей: Соя и рапс на масло семена 25 – 40 %, зерновые – 60 – 75 %. При этом рентабельность производства товарного зерна составила более 90 %.

Заключение

Для расширенного воспроизводства плодородия черноземных почв лесостепи Западной Сибири, по решению модели следует увеличить в структуре посевных площадей зернобобовые (соя, горох) и масличные культуры (рапс, подсолнечник) до 40%, при этом доля зерновых снижается до 65,5—78,0 %, а выход товарного зерна – до 2,06—2,11 т/га.

Товаропроизводитель может выбирать различные варианты структуры посевных площадей, например, снижение доли чистого пара до 14,3 % и пшеницы до 43,2 %, увеличение доли озимой ржи и технических культур (рапс, соя) до 21,3 %, позволяет получать 2,15 т/га продовольственного зерна. При использовании наряду с чистым паром занятого, сурепицу на сидерат, а также органическое удобрение солому зерновых культур,  получено 2,01 т/га товарного зерна, при условии что, под чистый пар отведено 16,8 %, под занятый – 11 %, под яровые зерновые – 72,2 % площади пашни. Увеличивая в севооборотах площадь с использованием в качестве органического удобрения соломы и сурепицы на сидерат, можно добиться воспроизводства плодородия почвы с некоторым снижением выхода товарного зерна до 2 т/га.

 

Библиографическая ссылка

 
Чибис В. В., Моделирование продуктивности севооборотов в условиях лесостепи Западной Сибири // «Живые и биокосные системы». – 2017. – № 20; URL: http://www.jbks.ru/archive/issue-20/article-5