Информац_технологии_Черных_2000
.pdf- 91 -
N 40000 / A 1 0.0015 A ,
Где A-возраст,лет.
Таблица 5.1
Коэффициент оптимальности числа стволов
Возраст,лет |
20 |
40 |
60 |
80 |
100 |
120 |
140 |
Кn |
0.98 |
0.94 |
0.90 |
0.86 |
0.82 |
0.78 |
0.74 |
Для сосны и ели оптимальное число стволов в зависимости от класса бонитета, возраста и оптимальной относительной полноты после 70 лет можно вычислить по следующим формулам:
N опт(1 бон)= (40 1600) / A K n Р,
N опт(1 бон)= (40 1800) / A K n P ,
N опт(1 бон)= (40 2000) / A K n P,
где A-возраст,лет;
P-относительная полнота после 70лет:для 1 класса бонитета равна 0,95,для 2 класса бонитета-0,85 и 3 класса бонитета-0,75;
Kn-коеффициент оптимального числа стволов в зависимости от возраста(табл.5,1)
5.1.2.3. Стратегия рубок ухода по К.Ф.Тюрмеру
На основе исследований известный лесовод К.Ф. Тюрмер ввел правило: в древостое после рубок ухода расстояние между деревьями должно быть равно 18-20 средним диаметрам (табл. 5.2). При проведении рубок ухода необходимо стремиться к соблюдению соотношения числа стволов и среднего диаметра древостоя согласно данным после проведения расчетов для определения интенсивности из-реживания конкретного древостоя.
Таблица 5.2 Оптимальное число стволов после рубки в зависимости от среднего
диаметра древостоя по К.Ф. Тюрмеру
D,см |
N,шт |
D,см |
N,шт |
D,см |
N,шт |
8 |
4822 |
16 |
1205 |
25 |
494 |
9 |
3955 |
17 |
1067 |
26 |
457 |
10 |
3086 |
18 |
952 |
27 |
423 |
11 |
2550 |
19 |
854 |
28 |
394 |
12 |
2143 |
21 |
700 |
29 |
367 |
13 |
1826 |
22 |
637 |
30 |
343 |
14 |
1574 |
23 |
584 |
31 |
321 |
15 |
1371 |
24 |
535 |
32 |
301 |
5.1.2.4.Стратегия рубок ухода по С.Н. Сеннову
Врезультате многолетних исследований и практического опыта
- 92 -
профессор С.Н.Сеннов предложил нормативы относительной интенсивности изреживания по запасу при рубках ухода и их повторяемости (табл.5.3).
Таблица 5.3 Показатели рубок ухода в сосновых древостоях таежной
лесорастительной зоны
|
Минимальная |
|
|
|
|
сомкнутость |
Интенсивность |
|
|
Возраст,лет |
полога, 0,1ед. |
рубки,в% от |
Повторяемость,лет |
|
|
До |
После |
запаса |
|
|
рубки |
рубки |
|
|
|
Осветления и прочистки |
|
||
10-15 |
1.0 |
0.8- |
10-20 |
5-15 |
|
|
0.7 |
|
|
|
|
Прореживание |
|
|
21-40 |
0.9 |
0.8 |
10-20 |
15 |
|
|
Проходные |
|
|
41 и более |
0.8 |
0.7 |
10-20 |
20-30 |
5.1.3. Результаты работы имитационной модели "Сосна"
Имитационная модель С», на" реализована на алгоритмическом языке TURBO - BASIC. Для запуска программы необходимо в командной строке набрать имя программы - "SOSNA99". После запуска программы пользователь ведет диалог с ПК, отвечает на поставленные вопросы.
Врезультате имитации рубок ухода необходимо получить наибольший запас в 80 лет.
Варианты имитации студент выбирает исходя из рекомендаций С.Н.Сеннова, С.В.Белова, К.Ф.Тюрмера или «Наставлений по рубкам ухода....» [60]. После анализа результатов имитации роста сосновых древостоев пользователь выбирает собственную стратегию интенсивности изреживания древостоя.
Втабл. 5.4-5.7 приведены примеры имитации роста сосновых древостоев, произрастающих на территории Нолькинского лесничества УНПКЛП МарГТУ. Изреживание высокополнотного древостоя сосны проводилось в 55 лет и 65 лет, с выборкой по запасу 4 м3 и 167 м3.
- 93 -
Таблица 5.4
Прогноз роста соснового древостоя
Область-Марий Эл
Лесхоз-УНПКЛПМарГТУ
Лесничество-Нолькинское Квартал-112 Выдел-3 Пробная площадь№9
|
. |
|
|
|
|
|
|
|
|
|
Число стволов,шт |
|
|
|
|
|
|
2 |
|
Возраст,лет |
Высота,м |
Диаметр,см |
Длина кроным, |
Видовое число,0,001 |
Сомкнутость полога, 1.0 |
Полнота,1.0 |
Суммаплощадей сечений,м |
3 |
|
Запас,м |
|||||||||
|
|
|
|
|
|
|
|
|
|
55 |
1000 |
25.0 |
27.0 |
|
490 |
0.933 |
1.32 |
57.2 |
701 |
55 |
995 |
25.0 |
27.0 |
5.0 |
490 |
0.929 |
1.31 |
56.9 |
697 |
65 |
852 |
27.6 |
29.6 |
5.5 |
480 |
0.916 |
1.27 |
58.7 |
778 |
65 |
742 |
27.3 |
28.2 |
5.2 |
483 |
0.745 |
1.01 |
46.4 |
611 |
75 |
639 |
29.9 |
31.8 |
6.1 |
473 |
0.789 |
1.05 |
50.8 |
718 |
|
|
|
|
|
|
|
|
|
|
Таблица 5.5 Распределение числа стволов по ступеням толщины
Возраст,лет |
|
|
|
Ступени толщины,см |
|
|
|
||
|
16 |
20 |
24 |
28 |
32 |
36 |
40 |
44 |
48 |
55 |
54 |
166 |
273 |
257 |
158 |
65 |
18 |
10 |
0 |
55 |
54 |
166 |
731 |
257 |
158 |
61 |
17 |
10 |
0 |
65 |
17 |
73 |
721 |
230 |
193 |
110 |
43 |
13 |
0 |
65 |
17 |
73 |
72 |
230 |
193 |
37 |
15 |
4 |
0 |
75 |
0 |
27 |
88 |
160 |
168 |
116 |
55 |
18 |
0 |
- 94 -
Прогноз роста сосны обыкновенной
Рис 5.1. Распределение числа стволов по ступеням толщины при прогнозе роста с учетом изреживания сосны обыкновенной
Прогноз роста соснового древостоя Область - Марий Эл Лесхоз – УНПКЛПМарГТУ
Лесничество— Нолькинское Квартал — 12 Выдел - 6 Пробная площадь № 11
Таблица 5.6
лет,Возраст |
стволовЧисло,шт |
см,Высота |
см,Диаметр |
кроныДлина,м |
Видовое ,0.001число |
Сомкнутость ,1,0полога |
,1.0Полнота |
площадейСумма м,сечения |
м,Запас |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
99 |
455 |
25.0 |
28.0 |
9.0 |
0.463 |
0.700 |
0.68 |
28.0 |
324 |
99 |
451 |
25.0 |
28.9 |
9.0 |
0.463 |
0.602 |
0.67 |
27.6 |
319 |
109 |
447 |
26.0 |
29.4 |
8.1 |
0.468 |
0.613 |
0.73 |
30.4 |
370 |
120 |
435 |
26.9 |
30.5 |
7.7 |
0.469 |
0.608 |
0.74 |
31.8 |
402 |
|
|
|
|
|
|
|
|
|
|
- 95 -
Таблица 5.7 Распределение числа стволов по ступеням толщины
Возраст, |
|
|
|
Ступени толщины,см |
|
|
|
|||
лет |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12 |
16 |
20 |
24 |
28 |
32 |
36 |
40 |
44 |
48 |
|
|
|
|
|
|
|
|
|
|
|
99 |
11 |
33 |
68 |
96 |
96 |
74 |
45 |
21 |
8 |
3 |
99 |
11 |
33 |
68 |
96 |
96 |
74 |
42 |
20 |
7 |
3 |
109 |
6 |
21 |
52 |
87 |
100 |
85 |
56 |
28 |
11 |
3 |
120 |
0 |
13 |
42 |
77 |
100 |
91 |
63 |
33 |
13 |
4 |
Общая закономерность распределения числа стволов по ступеням толщины с увеличением возраста сохраняется (рис. 5.1). К возрасту 75 лет древостой имеет запас 718 м3 на 1 га и полноту 0.789.
Результаты имитации не согласуются с рекомендациями К.Ф.Тюрмера, Интенсивность изреживания была занижена. Число стволов необходимо уменьшить в два раза (табл. 5.2). Для достижения требуемого результата по прогнозу роста сосны имитацию необходимо повторить.
В таблицах 5.6 и 5.7 приведены результаты прогноза роста с выборкой по запасу в 99 лет 5 м3. Этот пример имитации роста сосны отражает закономерности роста в возрасте спелости.
Следует отметить, что при наличии рядов распределения числа стволов по ступеням толщины можно рассчитать товарную структуру будущего древостоя.
5.2. Автоматизированное построение моделей роста Древостоев
5.2.1. Теоретические сведения
Таблицы хода роста (ТХР) имеют большое значение для лесоустройства. Сведения о ходе и производительности насаждений служат основой многим лесохозяйственным расчетам. Есть несколько видов ТХР, подразделяемых по назначению:
1)ТХР нормальных насаждений;
2)ТХР модальных насаждений;
3)ТХР оптимальных насаждений;
4)ТХР с различной густотой и разными суммами площадей сечений.
ТХР предназначены:
• для характеристики и прогнозирования роста и развития древостоев;
-
• установления спелости леса и обоснования возраста рубки;
-96 -
•проектирования лесохозяйственных мероприятий;
•выявления закономерностей роста и развития древостоев;
•составления региональных таксационно-справочных нормативов и тд. В разработку методики составления таблиц хода роста внесли большой
вклад видные лесоводы страны: Орлов М.М., Варгас де Бедемар, Тюрин А.В., ТретьяковН.В., АнучинН.П, АнтанайтисВ.В., Багинский В.Ф., Бочков И.М., Галицкий В.В., Дворецкий М.Л., ДракинВ.Н., Вуевский Д.И., ЗагреевВ.В , Зейде Б.Б., Кивисте А.К., Кор-сунь Ф.Д., Кофман Г.Б., Кузьмичев В.В., Макаренко А.А., Моисеев B.C., Моисеенко Ф.П., Мошкалев А.Г., Никитин К.Е., СаликовН.Я., Свалов Н.Н. и многие другие. ЗагреевВ.В. (1971,1978) на основании 400 таблиц хода роста выявил возможность систематизации таблиц. Оказалось, что для определения хода роста сосны, например, нужна всего одна таблица с тридцатью типовыми рядами. Такие индексные таблицы (графики) служат для сравнительной оценки и группировки таблиц по степени сходства и различия в характере (типе кривых) хода роста.
Наиболее известные таблицы хода роста оптимальных насаждений составили Е. Assmann и F. Franz (1967). Одно из новшеств в таблицах заключается в том, что в пределах каждого класса бонитета введены три уровня производительности (верхний, средний и низший). ТХР составляются для каждого уровня производительности.
Ход роста описывается множеством функций хода роста, которые были сгруппированы А.К. Кивисте (1988) в семь групп и проанализированы с учетом общих закономерностей роста живых организмов.
К функциям роста предъявляются следующие требования:
1)График должен проходить через начало координат, f(x) =0 прих=0.
2)Функция роста должна быть возрастающей,/'(х) больше или равно 0
при х>0.
3)Предел функции при неограниченном возрастании аргумента должен стремиться к асимптоте: lim fix) = max x —> оо
106
4)График текущего прироста должен исходить из начала координат,/'^ =0 при х=0.
5)Функция роста должна иметь точку nepera6a,/"(xj=0 прих>0. Анализ известных функций роста позволил выбрать функцию
А. Мичерлиха (Mitscherlich). Эта функция удовлетворяет необходимым требованиям для описания общих закономерностей роста живых организмов.
Функция имеет вид:
Т мод Т мах (1 е А С1 )С2
где Т - моделируемый таксационный показатель, диаметр (см), высота (м), сумма площадей сечений (м:), запас (м3);
Т - асимптотическое значение таксационного показателя для данmax
но го естественного ряда развития древостоя;
- 97 -
А -возраст, лет; е - основание натурального логарифма;
C1 - параметр роста;
С2 - параметр формы кривой.
Кривая роста имеет S-образный вид при С2> 1 с точкой перегиба m=lnC2/Cr.
т" Если С2 находится в пределах от 0 до 1, то функция не имеет точки перегиба. Такой вид кривой хорошо описывает значения сумм площадей сечений и диаметров, a S - образный вид кривой описывает значение высот и запасов. Следовательно, для моделирования основных таксационных показателей можно использовать один вид функции, изменяя параметры роста и формы кривой.
Следует отметить, что функция А. Мичерлиха приемлема для описания динамики изменения таксационных показателей отдельного дерева и древостоя.
5.2.2. Алгоритм моделирования хода роста древостоев
Функция роста А.Мичерлиха является нелинейной по параметрам. Для поиска оптимальных параметров С; и С2 использован итерационный метод, который можно описать по этапам:
А. В каждой точке естественного ряда по возрасту для моделируемых таксационных показателей вычисляется Су. при фиксированных значениях С2 (С2 = 1,2,3,4,...,10):
С |
|
C2 ln(1 Ton / Tmax ) |
, |
|
|||
1i |
|
A |
|
|
|
- 98 -
Где Топ( А) -фактическое значение таксационного признака в возрасте
А
I=1,2,3…,n;n-число наблюдений
Б.По наименьшей абсолютной ошибке функции выбирается значение параметра С2,которое принимается за константу.
В.При постоянном значении параметра С2 производится сглаживание параметра С1 для всех точек естественногои ряда по полиному третьей степени,С1=f(A).
Г.Параметр С2 расчиываем по уравнению:
С2 |
|
(Т |
оп |
/ Tmax ) |
, |
|
C , |
ln(1 e) |
|||||
|
|
|
||||
|
|
1 |
|
|
|
Где С’1-расчетные значения параметра С1 по полиному третьей степени. Д.Производится сглаживание параметра С2 по полиному третьей
степени,С2=f(A).
Е.Для моделируемого таксационного показателя решается задача поиска значения асимптоты при минимальной среднеквадратической ошибке
(m) на отрезке [a,b], F(X,C1,C)=Tвыкi m=((Tвычi-Ton)^2)/n-->min,где n-число экспериментальных значений моделируемого таксационнго показателя;Т м' ах-
наибольшее значение таксационного показателя из n экпериментальных
точек; a Tmax' 1.05; b Tmax' 1.30.
Оптимальное значение асимптоты таксационного показателя определяется по методу золотого сечения, который реализуется по следующему алгоритму:
1.Вычисляется коэффициент дробления отрезка [а,b], К=(5^1/2- 1)/2=0.6180339.
2.Вычисляется абсцисса Х1=а+(1-К)'(Ъ-а) и вычисляется F(X1,C1,C). Cj и С2 определяются согласно п. А.
3.Вычисляется абсциссаХ2=а+К'(b-а).
4.Проверяется условие (X2-XJ<E, где Е - заданная погрешность вычисления Хт При выполнении этого условия вычисляется Xm=(X; +XJ/2 yiF(Xm,CyCJ и останавливается счет.
Если условие не выполняется, идѐм к п. 5.
5.EcmFpCj. С;}, С) < F(X2, С,, С), то а=Хг Х=Х2, F(Xy Cp CJ = F(X2,
С,, С2) и идем к п. 3 и п. 4.
6.Если^(Т7, C]t CJ > F(XT Cj, CJ, то b=X, X2=X, f(X? Cr CJ =f(Xr CJ: CJ идѐм к п. 2 и п. 4.
- 99 -
Рис. 5.2. Модельная структура программы MICHXOD
Ж. Моделированию подвергаются четыре таксационных показателя: средний диаметр, средняя высота, сумма площадей сечений и запас.
- 100 -
Модель хода роста древостоя элемента леса принимает вид:
Tij Tmax (1 e A C1 y )C2 y ,
Гдеj=1-индекс высоты; J=2-индекс диаметра; J=3-индекс запаса;
J=4-индекс суммы площадей сечения; I=1,2,3…,n-индекс возраста с шагом 1,5 и 10 лет.
3
С1кij bkj' Akj ,
0
3
C2kij bkj'' Akj
0
Динамика отпада, общей производительности древостоя, прироста по запасу, среднего объема среднего ствола и его видового числа рассчитываются по общепринятым в лесной таксации формулам.
5.2.3. Программа моделирования хода роста древостоя
По приведенному алгоритму разработаны две программы в среде
"TURBO-BASIC'-MICHVVOD.BAS hMICHXOD.BAS.
Программа MICHVVOD предназначена для формирования файла экспериментальных данных на гибкий магнитный диск.
Структуры данных имеют вид:
•имя файла данных;
•название объекта (древесная порода, тип леса и т. д.);
•число экспериментальных точек (количество пробных площадей);
•значения: возраст (лет), высота (м), диаметр (см), запас (м3), сумма площадей сечений (м2).
Программа MICHXOD предназначена для автоматизированного построения модели хода роста древостоя на базе функции Мичерли-ха.
Модельная структура программы приведена на рис 5.2. Программа занимает 72 Кбайт дисковой памяти. Время счета по одному естественному ряду 1... 2 минуты.