Строительство сети лесовозных дорог, обеспечивающей доступ к лесным ресурсам, наряду с транспортировкой древесины являются одними из наиболее затратных статей расходов для большинства лесозаготовительных предприятий [8], что обусловливает необходимость разработки научно обоснованных и пригодных для применения на практике методов проектирования лесотранспортной сети и оптимизации транспортно-технологических процессов лесозаготовок. Оптимальной может считаться такая лесотранспортная сеть, которая обеспечивает эффективный доступ к лесным ресурсам и при этом требует минимальных затрат на её строительство, дальнейшее обслуживание и транспортировку древесины [4]. Обеспечение экономической эффективности лесозаготовительного предприятия и сохранение экологической безопасности на осваиваемых территориях требует разработки универсального метода поиска оптимальной лесотранспортной сети. При проектировании сети лесовозных дорог необходимо принимать во внимание специфику лесозаготовок, характеризующуюся неравномерностью запасов лесных ресурсов; разрозненностью лесных участков; однонаправленностью лесных грузопотоков; специфичностью грузов; ярко выраженной сезонностью; увеличивающимся расстоянием вывозки и т.д. [2]. Точность определения очередности строительства лесовозных дорог с учетом динамики лесоводственно-таксационных характеристик насаждений и условий их произрастания, выбор наилучшего расположения искусственных сооружений оказывают существенное влияние на рентабельность лесозаготовительного предприятия, а значит, должны быть учтены в процессе проектирования и положены в основу комплексной математической модели для решения задачи по поиску оптимальной схемы сети лесовозных дорог. [1]
С целью оптимизации процессов транспортного освоения лесосырьевых баз предложена методика поиска оптимального пространственного расположения лесовозных дорог на основе методов математического моделирования, в частности, инструментария теории графов (метода остовного покрывающего дерева).
Процесс проектирования лесотранспортной сети включает в себя постановку задачи и задание целевой функции на минимизацию затрат на строительство лесовозных дорог, после чего проводится оценка текущей ситуации с характеристикой местности (особенности рельефа; гидрологические и грунтовые условия; природные и искусственные преграды; наличие особо охраняемых природных территорий и мест обитания животных; другие факторы, оказывающие влияние на расположение дорожной сети и стоимость её строительства). [3]. На основе материалов лесоустройства, таксационных и картографических данных рассчитывается ликвидный запас древесины и определяется годовой объем вывозки леса, а также территориальное расположение транспортных путей, погрузочных пунктов, верхних и нижних складов.
С учетом имеющихся систем и комплектов машин рассчитывается стоимость строительства 1 км дороги для определенных условий местности; вводятся соответствующие коэффициенты, корректирующие стоимость участков дороги в зависимости от различных факторов; выбираются ключевые участки, к которым необходимо проложить транспортные пути. При помощи математической модели, основанной на инструментарии теории графов (ключевые участки являются вершинами связанного неориентированного графа, а возможные пути между ними – нагруженными ребрами, веса которых соответствуют стоимости строительства дорог), производится поиск оптимальной транспортной сети. На последнем этапе рассчитывается экономический эффект от реализации проекта и даются рекомендации к очередности освоения лесосек с учетом динамики лесного фонда.
Использование цифровой модели местности позволяет учесть лесоводственно-таксационные характеристики насаждений и особенности осваиваемой территории: участки с разными типами рельефа и местности, на которых стоимость строительства и содержания участка дороги существенно отличается от базовой; участки-линии (реки, ручьи, овраги), пересечение которых требует дополнительных расходов на строительство искусственных сооружений и введение соответствующих коэффициентов; ключевые участки: места концентрации древесины, пункты примыкания лесовозных дорог, места расположения дорожно-строительных материалов и т.д. [7]. Распределение территорий на несколько категорий обусловлено необходимостью учета затрат на строительство 1 км дороги, которые могут быть выражены в виде базового значения и корректирующих коэффициентов. Все данные сводятся к единой размерности – денежному выражению.
Места концентрации древесины выбираются посредством расчета имеющегося объема лесных ресурсов на потенциальных ключевых участках (территориях, имеющих некоторый запас древесины), а после сравнения их между собой выбираются наиболее предпочтительные с точки зрения экономической эффективности – максимально возможной прибыли от реализации древесины с участка за вычетом расходов на лесозаготовку:
(1)
где SiC – потенциальный доход, который может быть получен от реализации всего запаса древесины, находящегося на i-м участке (i-й вершине графа) за вычетом расходов на заготовку древесины, руб.; – объем i-й породы дерева на i-м участке (i-й вершине графа), м3; Qi – объем заготавливаемой древесины на i-м участке независимо от породы, 1 м3; Pn – цена реализации 1 м3 n-й породы дерева, руб.; C – расходы (в зависимости от комплекта лесосечных машин и рельефа) на заготовку 1 м3 древесины независимо от породы, руб.
Критерием выбора ключевых участков может явиться минимальный порог доходности, устанавливаемый предприятием, а после проведения соответствующих расчетов выбираются все участки, соответствующие выбранному критерию, и обозначаются как вершины связанного неориентированного графа, которые должны быть соединены оптимальной дорожной сетью – ребрами, суммарный вес которых будет минимальным. Критерием оптимальности является минимальная стоимость строительства лесотранспортной сети, соединяющей обозначенные участки (вершины). Общая стоимость строительства сети лесовозных дорог определяется как сумма весов ребер, образующих минимальное покрывающее дерево, связывающее все участки в единую сеть. Вес ребра соответствует стоимости строительства дороги на конкретном участке. В общем виде формула определения веса ребра имеет следующее выражение:
(2)
где wij – вес ребра, соединяющего i-ую вершину с j-й, руб.; Lij,общ – общая длина дороги, соединяющего i-ую вершину с j-й, км; Cстр – стоимость строительства 1 км дороги, руб.
Общая длина дороги может быть выражена следующей формулой:
(3)
где Lбаз – длина участка, строительство дороги на котором имеет базовую стоимость, км; Lсущ – длина участка, имеющего дорогу, пригодную для транспортировки древесины и не требующая строительства и ремонта, км; – длина участка, строительство дороги на котором отличается от базовой стоимости, км;
– длина участка, на котором расположена дорога, требующая восстановления (ремонта), км;
– повышающий коэффициент в зависимости от i-го типа местности на n-м участке дороги;
– понижающий коэффициент в зависимости от i-го состояния дороги, требующей ремонта, на i-м участке дороги.
В ходе строительства дороги возможно осуществление заготовки древесины с близлежащих участков, что может быть принято во внимание при проектировании дорожной сети, т.к. прибыль от реализации заготовленной таким образом древесины напрямую влияет на возможность нахождения оптимального решения [5]. Тогда формула расчета веса ребра будет иметь следующий вид:
(4)
где – вес ребра с учетом доходов от «промежуточной» заготовки древесины на участке дороги ij, руб; Vij,пр – прибыль от реализации древесины, заготовленной на «промежуточном» участке от i-й вершины до j-й.
В свою очередь расчеты Vij,пр включают в себя суммарную прибыль от реализации древесины, заготовленной на «промежуточном» участке, за вычетом стоимости строительства трелевочных волоков, погрузочных пунктов, расходов на лесозаготовку:
(5)
где – объем древесины, заготовленной на базе построенной дороги от вершины i до вершины j; Cij,общ – общие расходы, которые рассчитываются по следующей формуле:
(6)
где lвол – общая длина трелевочных волоков, км; cвол – стоимость строительства 1 км волока, руб; – количество погрузочных площадок, шт.; cпп – стоимость строительства погрузочных площадок, руб.; Qij – общее количество заготовленной древесины, на промежуточном участке от i-й вершины до j-й, м3; Cзаг – стоимость заготовки 1 м3 древесины, руб.
В итоге формула (4) примет следующий вид:
(7)
Ключевые участки могут быть выбраны по методу экспертных оценок и ранжированы на классы от 1 до 5 в соответствии с имеющимися запасами древесины, удаленностью, рельефом, строением почвы и т.д. После нахождения ключевых участков их условные обозначения переносятся на компьютерную платформу, пригодную для расчета оптимальных путей на основе теории графов. Кругами (вершинами графа) обозначаются участки леса, наиболее привлекательные для освоения с точки зрения экономической эффективности (цифрами обозначены номера участков – например, цифра «1» обозначает исходную точку – нижний склад/предприятие/железнодорожный тупик), а соединяющие эти вершины ребра являются возможными транспортными путями. Затем проводится расчет стоимости строительства каждого пути: определяется базовая стоимость строительства 1 км дороги (в млн руб.) и с учетом особенностей местности и климата вводятся дополнительные коэффициенты, корректирующие стоимость строительства всех участков дороги, после чего ребрам присваиваются соответствующие веса. Вес каждого ребра соответствует расходам на строительство дороги на определенном участке. После того, как выбраны вершины и соединены между собой возможными ребрами, запускается алгоритм построения минимального дерева в соответствии с формулой (7).
Минимальное остовное дерево (оптимальная транспортная сеть)
На рисунке показана схема построенного минимального остовного дерева (выбранная сеть выделена жирными линиями), суммарный вес которого соответствует стоимости строительства всей дорожной сети, соединяющий ключевые участки концентрации древесины. Окончательная стоимость строительства в дальнейшем может быть уточнена во время проведения полевых исследований.
Помимо экономического эффекта от лесозаготовок необходимо учитывать социальные и экологические составляющие [6], которые в общем виде можно представить формулой:
E = U + G + C – I, (8)
где E – общая эффективность проекта (экономическая, социальная, экологическая); U – выигрыш для пользователей транспортной инфраструктуры; G – выигрыш для экологии; C – прибыль от реализации заготовленной древесины; I – инвестиции в лесозаготовку. Общая экономическая эффективность капитальных вложений выражается формулой:
(9)
где D – сокращение текущих издержек реализации проекта; T – выигрыш транспортных предприятий; N – снижение расходов в народном хозяйстве; P – прирост валового объема различной продукции в отраслях народного хозяйства; V – выигрыш в социальной сфере; C – улучшение экологии; I – капитальные инвестиции.
Заключение
Последовательность освоения лесного массива может производиться двумя способами: посредством постепенного продвижения лесозаготовок вглубь лесного массива со строительством магистрали небольшой протяженности, обеспечивающей вывозку древесины в краткосрочный период или через строительство магистрали на всю длину (так называемая «глубинная схема»). Второй вариант является наиболее предпочтительным ввиду того, что объемы работ по дорожному строительству концентрируются в одном месте и обеспечивают эффективное применение машин, а единовременные общие капиталовложения позволяют сократить расходы на множество подготовительных работ, которые имели бы место в «пошаговом» освоении лесосырьевых баз. Кроме того, появляется возможность проведения комплексных работ по охране и защите леса на начальном этапе, однако к недостатку второго метода можно отнести большие капиталовложения на первом этапе. В дальнейшем необходимо разработать алгоритм определения очередности освоения ключевых участков в зависимости от сроков, необходимых на строительство дорог, заготовку, вывозку и реализацию древесины, а также наличия у предприятия свободных денежных средств для инвестирования в долгосрочные проекты.
Рецензенты:
Иванов В.А., д.т.н., профессор, академик РАЕ, заведующий кафедрой «Лесные машины и оборудование», ФГБОУ ВПО «Братский государственный университет», г. Братск;
Огар П.М., д.т.н., профессор, проректор по научной деятельности, ФГБОУ ВПО «Братский государственный университет», г. Братск.
Работа поступила в редакцию 14.02.2013