CC BY
32
УДК 519.688:630*245
UDK 519.688:630*245
ОПТИМИЗАЦИЯ ПАРАМЕТРОВ АВТОМАТИЗИРОВАННОЙ СИСТЕМЫ УПРАВЛЕНИЯ ВЫРАЩИВАНИЕМ СОСНОВЫХ НАСАЖДЕНИЙ
Мурзинов Юрий Валерьевич
Воронежская государственная лесотехническая академия, Воронеж, Россия
В статье рассмотрена возможность повышения эффективности автоматизированной системы управления выращиванием сосновых насаждений при помощи оптимизации режимов и параметров рубок ухода за лесом. Представлены регрессионные уравнения для определения параметров роста деревьев. Выделен критерий оптимизации рубок ухода в сосновых насаждениях. Показано увеличение выхода полезной древесины и снижения возраста технической спелости деревьев при использовании автоматизированной системы управления
Ключевые слова: АВТОМАТИЗИРОВАННАЯ СИСТЕМА УПРАВЛЕНИЯ, РЕГРЕССИОННЫЕ МОДЕЛИ, РУБКИ УХОДА, СОСНОВЫЕ НАСАЖДЕНИЯ, ВЫРАЩИВАНИЕ ЛЕСА
OPTIMIZATION OF PARAMETERS OF AUTOMATED CONTROL SYSTEM OF CULTIVATION OF PINE PLANTINGS
Murzinov Jury Valerjevich
Voronezh State Academy of Forestry and Technologies, Voronezh, Russia
In the article, the opportunity of increase of efficiency of automated control system of cultivation of pine plantings by means of optimization of modes and parameters of cleaning cutting of forestry is considered. The regression equations for definition of parameters of growth of trees are received. The criterion of optimization of cleaning cutting in pine plantings is allocated. The increase of the forest yield and reduction in age of technical maturity of trees is shown, with a help of automated control system
Keywords: AUTOMATED CONTROL SYSTEM, REGRESSION MODELS, CLEANING CUTTING, PINE PLANTINGS, CULTIVATION OF A WOOD
Повышение выхода деловой древесины и улучшение качества состава лесов в настоящее время является одной из основных проблем лесного хозяйства РФ.
Рубки ухода за лесом были и остаются самым трудоёмким, сложным и потому проблематичным лесохозяйственным мероприятием. Уход за лесом является распределённым во времени процессом на протяжении всего периода лесовыращивания и на каждом возрастном этапе имеет свои особенности, которые заключаются, прежде всего, в возрастной оптимизации площади питания каждого дерева [1]. Для эффективного планирования проведения рубок ухода за лесом необходимо иметь достаточно точные математические модели, наиболее полно описывающие процессы выращивания насаждений.
Повышение эффективности выращивания сосновых насаждений возможно за счёт оптимизации режимов и параметров рубок ухода,
которая основывается на анализе текущего состояния насаждения описываемого таблицами хода роста.
Таблицы находят применение в научных исследованиях при проектировании лесохозяйственных мероприятий, прогнозировании общей продуктивности лесов и установлении эффективности деятельности лесохозяйственных предприятий. Статистическая информация о ходе роста деревьев, заключенная в таблицах хода роста, может использоваться для построения эмпирических зависимостей, связывающих параметры роста деревьев и таксационные показатели насаждений. Обработка статистических данных методом наименьших квадратов позволяет построить регрессионные модели, наиболее точно описывающие процессы роста и развития насаждений.
В ФБГОУ ВПО «Воронежская государственная лесотехническая академия» на кафедре «Автоматизация производственных процессов» разработаны двухпараметрические регрессионные модели хода роста сосновых насаждений, а так же модели, описывающие реакции на проведение внутренних рубок ухода. Двухпараметрические регрессионные модели наиболее полно и точно отражают динамику развития насаждений.
Таблица 1
Коэффициенты регрессионной модели средней высоты в насаждениях
сосны обыкновенной
Бонитет Ш0 Ш1 Ш2 Ш3 Ш4 Ш5 Шб Ш7 Ш8
1б -2.908^10"2 3.957 4.837^ 10-2 -5.917 0.353 -0.299 5.363^10"4 -1.342^10"4 3.569^10"7
1а -1.38Ы0"2 2.149 3.115^ 10-2 -3.134 0.2 -0.17 -2.86^10"4 9.055^ 10-5 1.885^10"7
I -2.044^10"2 1.743 3.588^ 10-2 -2.771 0.347 -0.244 2.718^10"5 8.036^10"5 -7.563^10_/
II -1.376^ 10-2 1.738 2.002^10"2 -3.546 0.246 -0.152 -4.002^10"4 8.80Ы0"5 4.60Ы0"7
III -3. 174^ 10-2 1.828 1.838^ 10-2 -4.987 0.279 -9.713^10"2 -2.548^10"4 1.06Ы0"4 4.582^10"7
Для построения автоматизированной системы управления выращиванием одновозрастных сосновых насаждений получены следующие регрессионные модели:
Двухпараметрическая модель средней высоты дерева
H (Т, О) = m0 + шхТ + m2T 2 + m3 D + m4 D 2 +
+ m5TD + m6TD 2 + m7T 2 О + Ш8Т 2 О 2
(1)
где Н - средняя высота дерева, м; Т - возраст насаждения, лет; О - средний диаметр деревьев, м; ш0...ш8 - коэффициенты регрессионной модели, значения которых показаны в табл.1.
Таблица 2
Коэффициенты регрессионной модели количества деревьев в насаждениях
сосны обыкновенной
Бонитет к0 к: к2 к3 к4 К5 к6
1б 13347.38 -872.23 26.88 -0.46 4.448^10"3 -0.227^10"4 0.478^10"7
1а 14050.86 -915.23 28.44 -0.49 4.792^10"3 -2.469^10"5 5.208^ 10-8
I 14416.66 -895.17 26.81 -0.45 4.249^10"3 -2.134^10"5 4.406^10"8
II 13254.96 -553.28 9.46 -7.226^10"2 1.993^10"4 3.724^10"8 0
III 13235.77 -461.2 6.62 -4.308^10"2 1.044^10"4 0 0
Модель количества деревьев в насаждении
N = к0Т6 + к1Т5 + к2Т 4 + к3Т3 + к4Т2 + к5Т + к6 (2)
где N - количество деревьев в насаждении в возрасте Т, шт; к0...кб -коэффициенты регрессионной модели (табл.2).
Двухпараметрическая модель общего текущего прироста 7 оЬ (Т, Н) = х0 + х1Т + х2Т2 + х3 Н + х4Н 2 + х5ТН + х6ТН 2 + х7Т2 Н + х8Т2 Н 2 (3) где 7°ъ - общий текущий прирост насаждения, м3; х0...х8 - коэффициенты регрессионной модели (табл.3).
Двухпараметрическая модель прироста растущей части
7Роа (Т, Н) = у + УТ + у2Т2 + У3 Н + у4 Н 2 +
+ У5ТН + УбТН 2 + У7Т2 Н + У8Т2Н 2 где 7расГ - прирост растущей части насаждения, м3; у0...у8 - коэффициенты регрессионной модели (табл. 4).
Таблица 3
Коэффициенты регрессионной модели общего текущего прироста в насаждениях сосны обыкновенной
Бонитет Х0 Х1 Х2 Х3 Х4 Х5 Х6 Х7 Х8
1б 2.246^10-2 -2.479 1.336^10-2 6.402 -0.199 5.592^10-2 6.36^10-4 -1.262^10-4 -7.224^10-6
1а 1.242^ 10-2 -1.696 3.578^ 10-3 5.361 -0.164 4.151^10-2 4.792^10-4 -6.644^10-5 -2.335^10-6
I 3.949^10-3 -1.054 -1.414^ 10-2 4.277 -5.946^10-2 2.058^10-2 -4.435^10-5 3.214^10-4 4.565^10-6
II 0.139^ 10-2 2.04 -8.06^10-3 -9.034 -1.297 1.045^ 10-2 -1.717^ 10-4 -1.86Ы0"4 8.658^10-6
III 1.381^10-3 -0.958 -3.418^ 10-3 4.569 -0.127 3.419^ 10-2 -7.669^10-5 2.335^10-/ 5.934^ 10-7
Таблица 4
Коэффициенты регрессионной модели прироста растущей части в насаждениях сосны обыкновенной
Бонитет У0 У1 У2 У3 У4 У5 У6 У7 У8
Iб 6.408^10-2 -1.831 2.004^10-2 4.912 -0.227 6.38^10-2 6.911^10-4 -2.273^10- 4 -1.135^ 105
Iа 2.382^10-2 -1.169 1.219^10-2 3.929 -0.229 6.63^10-2 5.218^10-4 -3.019^ 104 -6.386^10- 6
I 5.772^10-2 -1.364 1.756^10-2 4.662 -0.263 4.756^10-2 1.198^10-3 -4.206^10- 4 -8.971^10- 6
II 3.986^10-2 -1.273 1.488^10-2 4.039 -0.234 7.006^10-2 9.526^10-4 -3.744^10- 4 -7.16^ 10-6
III -4.452^103 -1.06 1.826^10-2 5.318 1.222^10-2 -0.13 3.922^10-3 -1.494^ 104 -1.355^ 105
Полученные регрессионные модели (1), (2), (3) и (4) позволяют получить оптимальные показатели режимов и параметров рубок ухода. Для реализации оптимизации выращивания хвойных насаждений было получено уравнение, определяющее запас разреживаемого насаждения в возрасте Т
Т Т п
гт - Гто + I гГ + I (^ -г(Ко) - ±ьу} , (5)
*-Т0 * -Т0 ]-1
где Ут - запас разреживаемого насаждения в возрасте Т, м3; ¥То - запас
насаждения до проведения промежуточных рубок в некотором начальном периоде Т0, м ; К0=0,1...0,2% - среднегодовой процент отпада после проведения промежуточных рубок; Т1 - время проведения первой рубки, лет; п - количество рубок ухода, проведенных за период от Т0 до Т; Д V} -часть запаса, удаляемая при р-й промежуточной рубке и равная объёму
3
промежуточного пользования, м
= Р] • V,
где Р] - процент изреживания насаждения по запасу.
В настоящее время существуют методики оптимизации выращивания насаждений путём повышения эффективности проведения рубок ухода. Анализ восстановления полноты насаждений, проведённый С.Н. Сенновым [2] показывает, что прирост насаждения, пройденного правильными рубками ухода, исходя из природы снижения полноты, можно считать неизменным. При этом отпад оставленных в покое деревьев заметно снижается, и в идеальном случае практически равен нулю, что дает насаждению восстановиться до исходной полноты. С.Н. Сеннов утверждает, что при обоснованном проведении промежуточных рубок величина отпада составляет примерно 0.1% в год от величины общего запаса. Поэтому каждую последующую рубку назначают при восстановлении полноты насаждения до 1.
Ч. Система автоматизированного управления ускоренным выращиванием одновозрастных сосновых древостоев
Файл Оптимизация Вычисления Справка
Начальный возраст древостоя, Т лет
Запас насаждения, V мЗ/га
Бонитет, В
Целевой диаметр, Ои
Время первой рубки, Т1 лет
Таблица!. Режимы промежуточных и сплошной рубок
№ рубки Время ПРОЕ Объем уда; Удаляемая Выход пиле Диаметр д< Высота дер
1 25 41.68 21 0 12.1 12.6
2 34 63.03 21 0 16 16.7
3 45 85.37 21 0 20.5 20.8
4 58 84.55 13 0 25.6 24.7
5 74 581.67 0 31.5 28
Итого 866.33 0
Таблица2. Объем/Стоимость заготавливаемых от рубок сортиментов
№ рубки і 2 3| 4 5| По каждому виду сортимента !
Крупный пиловочник 1 то сорта 0/0 0/0 0/ 0 3.75/ 8.26 110.3/ 242.66 114 05/ 250.82
Крупный пиловочник 2-го сорта 0/0 0/0 о о о о 321.78/585.28 321.78/585.28
Крупный пиловочник 3-го сорта 0/0 0/0 о о о о о о о о
Средний пиловочник 1 то сорта 0/0 0/0 о о о о 0/0 0/0
Средний пиловочник 2-го сорта 0/0 0/0 о о о о о о о о
Средний пиловочник 3-го сорта 0/0 0/0 о о о о о о о о
Баланс 1-го сорта 0/0 0/0 о о о о о о о о
Баланс 2-го сорта 0/0 0/0 0/ 0 22.36/ 22.36 431.86/431.86 454 22/ 454.22
Баланс 3-го сорта 0/0 0/0 о о о о о о о о
Рудничная стойка 0/0 0/0 о о о о о о о о
Жерди 0/0 0/0 о о о о о о о о
Итого 0/0 0/0 0/0 26.12/30.62 863.84/1268.82 880.06/1300.44
4.74ХІІ
III
Рис. 1. Интерфейс системы автоматизации поиска оптимальных режимов и параметров рубок ухода
Процесс лесовыращивания развёрнут во времени и временные интервалы между рубками ухода являются главным фактором, влияющим на объём получения деловой древесины. Время проведения рубки в значительной степени зависит от полноты насаждения и запаса древесины.
На основе двухпараметрических регрессионных моделей хода роста (1)-(4) и уравнения (5) была разработана компьютерная программа для ЭВМ [3]. Интерфейс работы программы, показан на рис.2. В процессе работы программы варьируются проценты изреживания запаса насаждения Р] и возраст проведения первой рубки Т1, и, следовательно, возрасты проведения рубок ухода Тх. Оптимизация режимов и параметров рубок ухода будет проводиться по общему объёму заготавливаемой древесины
(7).
V(Т, Р] ) ® тах,
Т е Т . (7)
Р е р
Разработанная система автоматизированного управления выращиванием одновозрастных сосновых насаждений осуществляет имитационное моделирование хода роста сосновых насаждений (1б, 1а, I, II и III бонитетов) позволяет отобразить динамику их основных таксационных параметров. Преимуществами применения компьютерной оптимизации планирования рубок ухода является увеличение пользования древесиной, снижение возраста получения технически спелой древесины и улучшение качества древесных стволов.
Для сосновых насаждений Учебно-опытного лесхоза ФБГОУ ВПО «Воронежская государственная лесотехническая академия» с помощью программы «Система автоматизации инженерных расчётов режимов и параметров рубок ухода за лесом» были проведены оптимизационные расчёты. На рис.2 показана эффективность использования машинных
методов расчёта оптимальных режимов и параметров рубок ухода за лесом.
V, м3 1000
800
600
400
200
0
1а 1 2 3
□ Режим-оптимальный ПРежим-обычный Бонитет
Рис.2. Объём выхода деловой древесины при среднем возрасте главной рубки 82 года
Для наглядной оценки эффективности системы автоматизации инженерных расчётов режимов и параметров рубок ухода за лесом используем целевую функцию Q, характеризующую отношение увеличения объема получаемой от проведения рубок ухода древесины к сроку выращивания.
W = D-, м3/год, (8)
где А— - увеличение объёма запасаемой древесины после проведения рубок ухода по сравнению с контрольным объёмом запасаемой древесины.
Целевая функция (8) показывает, насколько увеличится пользование древесины с 1 га площади за 1 год лесовыращивания по сравнению с контрольным насаждением. Наиболее выгодным будут являться те режимы проведения рубок ухода, при которых W ® max .
По расчётам, проведенным для учебно-опытного лесхоза ФБГОУ ВПО «Воронежская государственная лесотехническая академия»,
пользование древесиной сосны увеличится в среднем на 3,2 м3/га в год (на 16,7 %), что показано на рис.3.
П, м3/год
Рис.3. Среднегодовое увеличение пользования древесины с 1га
Моделирование процесса лесовыращивания позволило определить оптимальные параметры автоматизированной системы управления этим процессом. Внедрение этого метода расчёта параметров позволит для действующих лесхозов получить дополнительный объём деловой древесины, снизить затраты на обслуживание техники и получить реальную экономию средств.
Литература
1. Малышев В.В., Петровский В.С. Разработка и проверка теории метода моделирования динамики распределенных технологических систем многолетнего формирования плантационных древостоев [Текст]. - Вестник Воронежского государственного технического университета. Том 2, №12, 2006. - С. 156-159
2. Сеннов С.Н. Уход за лесом (экологические основы). М.: Лесн. пром-ть, 1984. 128 с.
3. Патент 2010610420 РФ. Программа для расчёта режимов оптимального ускоренного выращивания сосновых древостоев: программа для ЭВМ. М.: Роспатент, 2010.
