Described software system for forest resources modeling on the case study of Irkutsk region on the basis of a model for the calculation of age classes over wide areas, taking into account the economic development of the territory. Various scenarios for the development of forests have been calculated, differing in the volume of natural and anthropogenic factors. A quantitative assessment of the impact of changes in anthropogenic factors on forest dynamics was made. Рассмотрен программный комплекс моделирования динамики лесных ресурсов Иркутской области на основе модели, предназначенной для расчетов по классам возраста на обширных площадях с учетом хозяйственного освоения территории. Рассчитаны различные сценарии развития леса, отличающиеся объемом влияния природных и антропогенных факторов. Проведена количественная оценка влияния изменения величины антропогенного воздействия на динамику лесов.
В работе проведено моделирование зависимости динамики лесных ресурсов от различных видов антропогенного воздействия. На основе результатных данных определена количественная оценка влияния изменения величины антропогенного воздействия на динамику лесов. Полученная оценка предназначена для формирования региональной стратегии использования лесных ресурсов.
Модель. Модель «Динамика древостоев» составлена на основе работ А.К. Черкашина [7]. Здесь динамика описывается системой дифференциальных уравнений. На моделируемой территории выделяются площади земель разных типов: нелесная, не покрытая лесом, молодняки, средневозрастные, приспевающие, спелые и перестойные. Динамика каждого участка описывается формулами вида (1): = − 0 ( ) + ( ); 0 = 0 ( ) − 01 0( ) + 0( ) + 0( ) − 0( ); (1) = −1 −1( ) − +1 ( ) − ( ) − ( ) − ( ), где aij – коэффициенты перехода из одной категории земель или группы возраста в следующую; SN – нелесная площадь; S0 – не покрытая лесом площадь; Si – лесные площади разных классов возраста; unon i – ежегодное увеличение нелесной площади за счет других категорий земель; uncov i – увеличение не покрытой лесом площади; ucut i – площадь рубок.
Ежегодное уменьшение объемов лесных ресурсов происходит из-за воздействия природных неблагоприятных факторов и антропогенного использования. Для лесных ресурсов Иркутской области существует ряд разрешенных видов использования: заготовка древесины; ведение сельского хозяйства; рекреационная деятельность; разработка месторождений полезных ископаемых, работы по геологическому изучению недр; строительство и эксплуатация линейных объектов; строительство и эксплуатация водохранилищ и иных гидротехнических сооружений; переработка древесины и иных лесных ресурсов.
Увеличение нелесной площади в процессе эксплуатации лесов представляется формулой (2):
= ∆ + ∆ + ∆ + ∆ + ∆ + ∆ , (2) где kN – площадь поселков на одного человека, остальные коэффициенты характеризуют увеличение: ∆N – населения лесничества, ∆S – сельскохозяйственных площадей, ∆R – рекреационных зон, ∆G – площади месторождений, ∆Bl – строительства линейных объектов, ∆Bv – обслуживания гидротехнических сооружений.
Переход других категорий земель в нелесные производится случайным образом, в зависимости от потребностей производства. Исходя из этого, распределение unon по остальным категориям земель полагаем пропорциональным текущей площади каждого участка (3):
= ∗ ∑ (3) Ежегодно леса подвергаются воздействию набора природных неблагоприятных факторов. На ослабление и гибель лесных насаждений Иркутской области влияют: пожары; болезни леса; повреждение насекомыми; неблагоприятные погодные условия; другие антропогенные и непатогенные факторы.
Сильнее всего на усыхание и гибель лесных ресурсов влияют пожары (43,7%), болезни леса (25,2%) и повреждения насекомыми (16,2%), поэтому именно эти факторы целесообразно учитывать в расчетах (4):
= + + , (4) где Sg – площадь лесов, погибших от пожара; Snas – площадь лесов, поврежденных насекомыми; Sb – площадь лесов, погибших от болезней.
Распределение uncov по категориям земель считается аналогично (3).
Программный комплекс. Программный комплекс реализован на языке программирования Java. Расчетный блок загружает исходные данные для моделирования из файлов формата Microsoft Excel. В них представлены сведения по лесничествам: распределение площади по категориям земель и классам возраста, объемы рубок, запасов, лесопосадок, количество населения, значения параметров антропогенного использования и природных факторов. При расчетах показатель количества населения увеличивается с каждым годом на известную величину ежегодного прироста численности населения. Значения параметров антропогенного использования и природных факторов для каждого лесничества считаются постоянными на всем периоде моделирования.
После запуска программы пользователь задает начальные значения – период и шаг расчетов. Численное решение системы дифференциальных уравнений производится методом Рунге-Кутта четвертого порядка. Комбинация параметров – объемы рубок, лесопосадок и воздействия неблагоприятных условий – задается в процентном отношении к имеющимся на текущий момент и образует сценарий управления ресурсом.
Рис. 1. Интерфейс программы
Рис. 2. Графики динамики лесных ресурсов по 5 сценариям Мы видим, что при текущих уровнях природного и антропогенного воздействия в долгосрочном периоде площади спелых и перестойных лесов уверенно растут, как и нелесные земли. Леса младших категорий возраста плавно убывают первые 100 лет, затем остаются на одном уровне. Во втором сценарии наблюдаем резкий рост спелых и перестойных лесов при небольшом снижении запасов лесов младших возрастных категорий и не покрытых земель. В третьем сценарии отмечен рост лесных площадей всех категорий при уменьшении не покрытых лесом.
В четвертом сценарии освоение лесосеки приведет к заметному снижению уровня спелых и перестойных лесов на протяжении первых 50 лет при повышении площадей младших возрастных категорий. В последнем сценарии увеличении объема рубок компенсируется одновременным увеличением лесопосадок и уменьшением неблагоприятных факторов: запасы спелых и перестойных в начале периода снижаются менее заметно, затем начинают возрастать.
Результаты моделирования выведены на карту Иркутской области, в качестве показателя взята суммарная площадь лесных ресурсов по всем категориям возраста. На рисунке 3 цвета 1-2 соответствуют отрицательным значениям относительно начального – убывание, цвет 3 – незначительные колебания, цвета 4-5 – положительные значения, рост площадей. Здесь видно, что антропогенному влиянию больше всего подвержены лесничества, расположенные в центральной части области. В больших лесничествах на севере и востоке области из-за низкой транспортной доступности лесные площади сохраняются на одном уровне, независимо от сценария.
Рис. 3. Карты с результатами расчетов по сценариям Рассмотрим сводные данные по изменению общей площади леса при расчетах по этим сценариям. Результаты сравнения приведены в таблице. При расчете по текущим уровням к концу периода моделирования лесные площади спелых и перестойных деревьев увеличатся на 13,1%, леса младших возрастных категорий уменьшатся на величину от 6,7% до 10,7%. При уменьшении воздействия неблагоприятных факторов в 2 раза площади спелых и перестойных возрастут на 19,8%. При увеличении лесопосадок в 2 раза площадь последней возрастной категории увеличится на 13,9%.
Увеличив рубки в 2,5 раза, мы получим снижение площадей спелых и перестойных на 17,9%. В последнем сценарии в этой категории наблюдается снижение на 13,1%. -9,9 13775 13,1 Природные/2
Посадки *2 Вырубка *2,5 Вырубка *2,5; посадки *2; природные/2 10827 11985 14776 14203 -17,4 -8,5 12,7 8,4 11733 12482 15387 15313 -14,8 -9,4 11,7 11,2 5734 5963 7323 7423 -8,9 -5,3 16,3 17,8 29560 28098 20245 21428 19,8 13,9 -17,9 -13,1 Результаты моделирования показали, что увеличение объемов рубок для освоения расчетной лесосеки приведет к заметному сокращению лесов. Если же увеличение рубок дополнить увеличением посадок деревьев и усилением борьбы с природными факторами, особенно лесными пожарами и болезнями леса, то компенсировать сокращение лесов получится гораздо раньше.
Заключение. Разработанный программный комплекс позволяет моделировать динамику лесных ресурсов с учетом воздействия комплекса природных и антропогенных факторов. Он также помогает принимать управленческие решения в лесохозяйственной отрасли, показывая направление изменения площадей конкретной категории земель в зависимости от заданного сценария.
Работа выполнена при финансовой поддержке Интеграционной программы ИНЦ СО РАН. Результаты получены при использовании ЦКП «Интегрированная информационно-вычислительная сеть Иркутского научно-образовательного комплекса» (http://net.icc.ru). [9] Vladimirov I.N., Chudnenko A.K. Forecasting of the spatio-temporal dynamics of the forest resources of the Irkutsk region using GIS technologies. Sun. Earth, Water, and Energy.
Transactions of the Siberian Division of APVN, issue 2. Novosibirsk, Nauka, 2005, 61–68. [10] Cherkashin Е.А., Badmatsyrenova S. B., Vladimirov I.N., Popova A.K., Davydov A. An optimal control module of sustainable natural resources consumption control synthesis for decision support systems // 37th International Convention on Information and Communication Technology, Electronics and Microelectronics (MIPRO), DOI: 10.1109/MIPRO.2014.6859733, 2014. pp. 1100-1105.