Generación de patrones de corte óptimos para árboles individuales a partir de productos demandados en plantaciones comerciales

Herrera, J.F., Barrios, A., López, A.M. & Nieto, V. (2015). Generación de patrones de corte óptimos para árboles individuales a partir de productos demandados en plantaciones comerciales. Colombia Forestal, 18(2), 193-206.

GENERACIÓN DE PATRONES DE CORTE ÓPTIMOS PARA ÁRBOLES INDIVIDUALES A PARTIR DE PRODUCTOS DEMANDADOS EN PLANTACIONES COMERCIALES

Generating optimal tree cutting patterns from commercial plantationsa

Jonnier F. Herrera¹, Alonso Barrios², Ana M. López³ & Víctor Nieto⁴

¹ Corporación Nacional de Investigación y Fomento Forestal, Bogotá, Colombia. Jfabian88@gmail.com Autor para correspondencia.
² Facultad de Ingeniería Forestal, Universidad del Tolima, Ibagué, Colombia. abarriost@ut.edu.co
³ Facultad de Ingeniería Forestal, Universidad del Tolima, Ibagué, Colombia. amlopeza@ut.edu.co
⁴ Corporación Nacional de Investigación y Fomento Forestal, Bogotá, Colombia. victornieto@conif.org.co

Recepción: 27 de octubre de 2014 / Aprobación: 13 de abril de 2015

RESUMEN

En el manejo de plantaciones forestales, la decisión sobre el trozado de árboles es importante, ya que está relacionada directamente con la rentabilidad de la operación. Para maximizar los beneficios económicos derivados del negocio forestal, es necesaria una adecuada y eficiente utilización de los fustes, pues existen diversas combinaciones posibles para dimensionar y asignar los diferentes productos forestales, de acuerdo con las exigencias del mercado. El objetivo de este estudio es desarrollar una metodología basada en técnicas heurísticas y de programación dinámica para la optimización del trozado de fustes individuales que permita maximizar el valor económico del árbol. Los modelos desarrollados se evalúan poniéndolos en práctica en el trozado de fustes de Eucalyptus tereticornis, que crece en la costa atlántica colombiana. Se utilizó un algoritmo para la generación de puntos potenciales de corte a lo largo del fuste comercial del árbol, y un algoritmo basado en programación dinámica determinó la combinación óptima de productos que debe ser obtenida en el fuste evaluado, imponiéndose al modelo restricciones debidas a las dimensiones de los productos demandados. Los resultados obtenidos fueron comparados con la simulación de un trozado secuencial, en la cual no se encontró ninguna diferencia entre los dos métodos, debido principalmente a la marcada diferenciación entre productos empleados en el ejemplo práctico. La metodología desarrollada puede ser utilizada para la optimización del trozado de cualquier especie que cuente con funciones de ahusamiento que permitan predecir diámetros y alturas fustales y volúmenes comerciales a lo largo del fuste.

Palabras clave: algoritmo, modelos de ahusamiento, Eucalyptus tereticornis, optimización de trozado, programación dinámica.

ABSTRACT

In the management of forest plantations decisions on tree bucking are important as they directly affect profitability. Proper and efficient use of stems is necessary in order to maximize the economic benefits derived from the forestry sector. There are various possible combinations for sizing and assigning different forest products according to market requirements. This study aims to develop a methodology based on heuristics and dynamic programming for bucking optimization of individual stems to maximize the economic value of the tree. The evaluation of bucking models was conducted through a practical application on stems of Eucalyptus tereticornis growing on the Colombian Atlantic coast. An algorithm for generating potential points along the commercial bole of the tree was used and an algorithm based on dynamic programming determined the optimal combination of products to be obtained in the evaluated stem, by imposing constraints on the model based on the size of the demanded products. A comparison of the results obtained with a simulation of sequential bucking did not find any difference between the two methods, mainly due to the marked differentiation between products used in the practical example. This methodology can be used to optimize the bucking of any species that has taper functions to predict upper-stem diameters and heights and commercial volumes along the stem.

Key words: algorithm, taper models, Eucalyptus tereticornis, bucking optimization, dynamic programming.

INTRODUCCIÓN

Generar técnicas precisas para estimar los volúmenes comerciales y la distribución de productos derivados de las plantaciones forestales es uno de los principales objetivos del manejo forestal, así como buscar la maximización del volumen de productos obtenidos y aumentar los ingresos económicos. Por esta razón, es necesario una adecuada y eficiente utilización de los fustes, pues existen diversas combinaciones posibles para decidir cómo dimensionar y asignar los diferentes productos forestales específicos (Arce et al., 2004).

Dentro de las etapas de transformación (trozado y aserrío), la asignación eficiente presenta el mayor potencial de mejoramiento, ya que las decisiones tomadas en estas operaciones influyen directamente en el nivel de volumen maderable aprovechado (Garzón, 2009). Actualmente, las decisiones de corte para la obtención de productos recaen sobre el operario encargado, quien define la mejor manera de trozar de acuerdo con su experiencia, por lo que se establece de esta forma un sistema de transformación empírico. Estudios demuestran una pérdida de valor comercial del fuste entre el 10% y 35% al usar un trozado normal o empírico frente a un corte óptimo (Cancino, 1993).

La decisión de trozar es compleja, más aun cuando se parte de las múltiples combinaciones que se generan de la relación entre fustes y productos con diferentes dimensiones que demanda el mercado. Generalmente, las decisiones de trozado no son las mejores, aún más cuando se aumenta la variedad de productos; por ello, para el óptimo manejo de los recursos forestales es necesario implementar técnicas de investigación operativa u otras técnicas de planificación que contribuyan a resolver los problemas de perdida de volumen comercial, en este caso en el proceso de trozado. Para Gilmore & Gomory (1961), esta dificultad se enmarca dentro de los denominados problemas de asignación o de corte (Cutting stock problem), mientras que Laroze (1999)lo denomina como problema de optimización de trozado (Bucking optimization problem) en el sector forestal.

Entre los primeros estudios realizados en el tema, Pnevmaticos & Mann (1972) abordan el problema de trozado por medio de programación dinámica (PD), que consiste en asociar una secuencia de etapas representadas por la longitud mínima de una serie de productos demandados y los estados asociados a las etapas, representados por la longitud de cada producto, donde estos son múltiplos de la longitud mínima del producto (Briceño, 2005).

Faaland & Briggs (1984) asocian los procesos de trozado y aserrío mediante la PD para maximizar la utilidad económica tanto de las trozas obtenidas por fuste, así como la de productos escuadrados. Cancino (1993) implementó un modelo optimizador de trozado basado en funciones de ahusamiento como entrada para estructurar la solución mediante un algoritmo de programación dinámica. Esta herramienta también fue implementada por Arce et al. (2002) para la especie Pinus taeda L., junto con un algoritmo heurístico; el primero determinaba la combinación de trozas por producto a obtener, mientras que la heurística definía los cortes potenciales a realizar a lo largo del fuste.

Kivinen (2007) asocia la programación dinámica para determinar los largos de troza óptimos de acuerdo con la demanda y tipos de productos rollizos, y resuelve el problema de la asignación óptima de productos mediante la aplicación de algoritmos genéticos, teniendo en cuenta restricciones de demanda, volumen en el inventario de fustes, precio del transporte de los productos, entre otros.

Además de la utilización de la programación dinámica en la solución de este tipo de problemas, se ha empleado la programación lineal (PL). Smith & Harrell (1961) y Jackson & Smith (1961), a partir de una tabla de rodal, generaron patrones de corte factibles por clase diamétrica, en donde la variable de decisión fue el número de árboles a asignar un tipo de corte, teniendo en cuenta restricciones de volumen demandado y oferta de materia prima en campo. En ese mismo sentido, Marshall et al. (2006) aplicaron programación entera mixta en la planificación de la producción de artículos según las restricciones de mercado e inventario de volumen aprovechable en el patio de apilado de fustes, y obtuvieron como resultado patrones de trozado óptimos de acuerdo con el beneficio económico.

Laroze (1999) maximizó el valor global de una muestra de fustes cosechados, sujeto a restricciones, especificaciones de exportación de productos y calidad de individuos mediante la implementación de una combinación de técnicas entre la programación lineal y la meta heurística denominada búsqueda de tabú (tabu search): mientras el modelo de PL genera un conjunto de patrones de corte de acuerdo con las restricciones planteadas, la meta heurística resuelve el problema de asignar los patrones de corte a cada fuste, tras lo que se obtuvo un nivel de eficiencia de alrededor del 97%.

La literatura también reporta el uso de combinación de las técnicas de PD y PL planteada por Eng & Daellenbach (1985) y Sessions et al. (1989), metodología denominada Dantzig–Wolfe, en donde la PD genera la maximización del beneficio por individuo originando patrones de corte óptimos para cada uno de ellos, mientras que la PL maximiza el valor del rodal al asignar un patrón de corte determinado por individuo.

En Colombia, la estimación de la oferta volumétrica de las plantaciones forestales comerciales se hace a través de modelos de crecimiento y rendimiento. Sin embargo, estas herramientas a menudo generan información a nivel de rodal y/o de árbol entero, lo que no permite la segregación de los volúmenes comerciales a nivel productos. Esto se traduce en incertidumbre sobre la cantidad de productos ofertados por el bosque y lleva a ineficiencias en la planificación de la cosecha. Teniendo en cuenta que este tipo de investigaciones han sido poco abordadas en el país, este trabajo tiene por objetivo el desarrollo de una metodología basada en técnicas heurísticas y de programación dinámica para la optimización del trozado de fustes individuales que permita maximizar el valor económico del árbol. Este modelo de optimización puede ser incorporado en sistemas de simulación de crecimiento y rendimiento y de procesamiento de inventarios forestales. En el estudio se presenta un caso práctico de aplicación de los algoritmos de trozado para árboles de Eucalyptus tereticornis Sm.

MATERIALES Y MÉTODOS

Modelo de optimización de trozado

El modelo empleado en el presente estudio se fundamenta en un enfoque basado en investigación operativa (Hillier & Lieberman, 1997), en donde se establecen las relaciones entre los elementos focales y los predictores con el fin de determinar el sistema en sus entradas y salidas (Figura 1). El modelo de optimización está diseñado para operar a partir de la tabla de rodal obtenida de inventarios forestales o producto de un simulador de crecimiento y rendimiento. Para cada clase diamétrica de la tabla de rodal, los diámetros índices de utilización, alturas fustales y volúmenes comerciales son predichos a través de una función de ahusamiento para la especie específica. El modelo de optimización de trozado también requiere información de los productos demandados como dimensiones (largo y diámetro mínimo) y precios de mercado (Figura 1).

El modelo está integrado por un generador de puntos potenciales de corte a lo largo del fuste el cual utiliza un algoritmo heurístico basado en la formulación de Arce et al. (2004), y un algoritmo de programación dinámica que identifica la combinación óptima de productos a ser obtenidos del fuste (Pnevmaticos & Mann, 1972). El modelo de optimización maximiza económicamente el valor de fustes individuales mediante la generación de patrones de corte óptimos, a partir de las dimensiones de productos, precios y dimensiones (D y altura total) de los fustes evaluados (Figura 1). La metodología empleada no consideró la aplicación de descuentos por defectos en el fuste, bifurcaciones, curvaturas, etc., y, por lo tanto, parte del supuesto de que los árboles a trozar presentan las siguientes características: a) tienen un fuste recto y no presentar bifurcaciones y b) no tienen daños mecánicos o problemas fitosanitarios.

Generación de puntos potenciales de corte

Con el fin de aumentar la eficiencia del proceso de optimización, se utilizó una heurística similar a la propuesta por Arce et al. (2004) para disminuir el número de puntos potenciales de corte. La Figura 2 muestra el diagrama de flujo del algoritmo empleado considerando tres productos. Si se tiene un fuste de longitud comercial (L_c), descontando la altura del tocón y hasta un diámetro índice de utilización igual al menor diámetro mínimo solicitado: d_i ≥ d_min, y considerando m productos de diferente largo l₁, l₂,…, l_m, es decir l_i = {1, 2,…, m}, a través de este algoritmo se obtienen todos los posibles puntos potenciales de corte. El algoritmo busca puntos factibles de corte desde la base hasta la altura comercial del fuste, iniciando por el producto de mayor longitud (l_m), y evaluando las posibles combinaciones con los productos de menor longitud (l_m-1, l_m-2).

Generación de patrones de corte óptimos

De acuerdo con las características de la PD, cada una de las etapas (X_n) tiene asociado un número de estados (s) que contienen información importante (precio de venta, longitud y diámetro mínimo admisible) en la búsqueda del óptimo. En este sentido, los estados están representados por el número de productos (m) a obtener del fuste (Figura 3).

Una vez definidas las etapas, los estados asociados y parámetros del sistema, como diámetro, longitud comercial del fuste y dimensiones mínimas aceptables de productos demandados, se plantea la función de recurrencia (1). Por medio de esta ecuación se espera obtener el patrón óptimo y algunos factibles de trozado para las diferentes marcas de clase que conforman la distribución diamétrica del rodal a cosechar.

Donde n es el índice de la etapa actual, n - 1 es la etapa previa, s es el estado del sistema en la etapa actual, x_n es la variable de decisión en la etapa actual, es el valor total óptimo de la mejor combinación de productos obtenida para la etapa x_n y el estado s, es la función de valor para la decisión x_n cuando inició en el estado s, es el valor de la mejor combinación para la etapa n - 1 para la decisión x_n cuando empezó en el estado s.

La función de recurrencia planteada en [1] solo tiene en cuenta especificaciones de precios y longitud, por tanto, es necesario agregar una restricción a partir de los diámetros mínimos aceptables por producto l_i y calcular el diámetro en el punto x_n mediante un modelo de ahusamiento, de esta forma el producto l_i, se acepta siempre y cuando:

De esta forma, con la ecuación [2] se descartan patrones de corte no factibles dado que por especificaciones de diámetro no pueden ser asignados en esa sección del fuste, siendo el diámetro del producto requerido mayor que el diámetro encontrado en la sección evaluada.

Comparación con un trozado secuencial

Para efectos de evaluación y comparación del modelo de optimización de trozado desarrollado, se obtuvo una simulación de trozado secuencial o por prioridad. Este sistema de trozado simula la decisión que debe tomar un operario al momento de trozar un árbol. En una primera instancia, el operario tratará de obtener del fuste el producto de mayor valor; una vez este sea agotado, continuará buscando el producto que le sigue y así sucesivamente hasta que se agote la longitud comercial del árbol. Por lo tanto, este sistema de trozado solo tiene en cuenta las dimensiones de los productos a obtener, mientras el precio de los productos solo es empleado para fijar el orden de prioridad de cada uno.

Ejemplo práctico

A modo de ilustración se presenta una aplicación práctica de los algoritmos de trozado desarrollados en la cuantificación de productos comerciales a partir de tablas de rodal en labores de inventario de plantaciones de Eucalyptus tereticornis Sm. La información base para este ejemplo práctico ha sido compilada por la Corporación Nacional de Investigación y Fomento Forestal (CONIF) en alianza con la Corporación para el Desarrollo Forestal de la Costa Norte de Colombia. Los datos empleados en el estudio proceden de una plantación de E. tereticornis, a la edad de 8 años, ubicada en la Reforestadora San Sebastián (RESS), en el departamento de Magdalena, Colombia. En el rodal fue instalada una parcela permanente de 500 m², en la cual todos los árboles presentaron condiciones ideales, sin ningún tipo de defecto. A partir de las mediciones de diámetro a la altura del pecho y altura total, se obtuvo una tabla de rodal, la cual incluye variables como: diámetro a la altura del pecho (D), altura total (H) y número de árboles por hectárea (N) (Tabla 1).

El modelo optimizador de trozado se implementó para árboles representativos de cada marca de clase de la tabla de rodal de los datos de entrada (Tabla 1). Debido a que el modelo optimizador de trozado requiere estimaciones de diámetros y alturas fustales, así como volúmenes comerciales, se requiere el empleo de una función de ahusamiento. En este estudio se empleó el modelo fustal polinomial simple desarrollado por Bruce et al. (1968), ajustado por López et al. (2011) para la determinación de volúmenes comerciales de E. tereticornis, desarrollado a partir de perfiles fustales de la especie en las mismas plantaciones objeto de estudio, en la costa atlántica colombiana. Este modelo de ahusamiento permite determinar la posibilidad de asignar un tipo de producto de acuerdo con el diámetro de la sección del fuste en donde se esté evaluando esta posibilidad.

Donde d es el diámetro fustal sin corteza (cm) a una altura fustal h, D es el diámetro a la altura del pecho con corteza (cm), H es la altura total del árbol (m), y X=(H-h)(H-1.3).

Para la estimación de volúmenes comerciales a nivel de productos, así como del desperdicio, se calculó el volumen fustal y el volumen de las trozas aplicando técnicas de integración numérica de la función [3]. La altura comercial se determinó como aquella altura fustal en donde se pueda asignar el producto con el diámetro mínimo requerido de los productos demandados de transformación primaria (Tabla 2), por lo que la longitud corresponde entonces a un diámetro de 6 cm.

RESULTADOS

Generación de puntos potenciales de corte

La generación de una secuencia de etapas o puntos factibles de corte se implementó en primera fase a manera ilustrativa para el árbol medio de la clase diamétrica cinco de la Tabla 1. Considerando un árbol con altura media de 26.4 m, el cual presenta una altura comercial (hasta el índice de utilización 6 cm) igual a 22.9 m (evaluada empleando la función de ahusamiento [3]) y descontando 0.2 m de tocón, quedaría L_c = 22.7 m. De acuerdo con la metodología planteada por Pnevmaticos & Mann (1972) y Sessions et al. (1989), para calcular los puntos potenciales de corte se debe considerar el mínimo común divisor de los productos a asignar; en este sentido, los tres productos a evaluar en la presente investigación son: l₁ = 8 m, l₂ = 4 m y l₃ = 2.5 m, donde el mínimo común divisor entre ellos es de 0.5 m. Para el fuste comercial evaluado L_c = 22.7 m surgen (22.7 m/ 0.5m) ≈ 45 puntos potenciales de corte. Por otro lado, al utilizar el algoritmo generador de puntos potenciales de corte se reduce a 32 puntos, lo que equivale a una disminución del 29% (Figura 4). Esta reducción podría aumentar la eficiencia del algoritmo de programación dinámica en la búsqueda del máximo beneficio, ya que este tendría que evaluar menos puntos de corte factibles a lo largo del fuste o menos etapas (x_n).

Obtención de patrones de corte óptimos

La Figura 4 ilustra un diagrama de red que relaciona los patrones de corte generados a partir de la asociación de los puntos potenciales de corte, donde cada nodo representa las diferentes etapas y los arcos contienen información asociada a cada estado (producto: precios, longitud y diámetros admisibles). Es importante recalcar que muchos patrones de corte no son factibles, ya que son producto de las combinaciones entre longitud de producto a obtener. Una vez evaluada la restricción impuesta por la ecuación (2), se descartan los patrones de corte que no son factibles.

Para seleccionar el patrón de corte óptimo se implementó la función de recurrencia planteada en la ecuación (1), a fin de seleccionar la ruta con mayor beneficio económico; dicha solución se construye a partir de la sumatoria de los óptimos encontrados en cada etapa (Tabla 3). En este caso, a modo de ejemplo, se evaluó el algoritmo de PD para un individuo con altura total de 22.8 m y D = 24.9 cm, el cual corresponde a la clase diamétrica tres de la Tabla 1. A partir de la función de ahusamiento se predijo una altura comercial de 18.9 m (altura hasta un diámetro límite de utilización de 6 cm), y descontando la altura del tocón de 0.2 m, queda un fuste utilizable de 18.7 m.

En la Tabla 3, la última celda subrayada corresponde a la etapa 20 (x_n = 17 m) asociada al estado s = 1. Este valor corresponde al máximo beneficio económico obtenido al trozar el fuste, con un sobrante de 1.5 m (8% longitud fustal). A partir de este punto se construye la trayectoria para establecer los tipos de productos a extraer; en este sentido, los superíndices en la etapa 20 muestran el tipo de productos a ser retirados: 1 unidad del producto l₁ (l = 8 m y d_min = 16 cm), 1 unidad del producto l₂ (l = 4 m y d_min = 12 cm), 2 unidades del producto l₃ (l = 2.5 m y d_min = 6 cm), es decir, 0.2 (tocón) + 8 + 4 + 2.5 + 2.5 = 17.2 m, con un retorno de utilidad de 64 280 COP.

Comparación con un trozado secuencial

Una simulación de un modelo de trozado secuencial para el mismo árbol de la clase diamétrica tres (Tabla 1), teniendo como prioridad la obtención de postes de 8 m, limatones de 4 m y, por último, varas de 2.5 m, encontró el mismo resultado que el modelo optimizador de trozado (Tabla 4). Este resultado se debe a la clara diferenciación entre los productos, de acuerdo a sus largos, y especialmente a sus precios de mercado.

Aplicando la función de recurrencia planteada y los puntos potenciales de corte a los diferentes árboles que representan la marca de clase en la distribución diamétrica del rodal de estudio, se obtuvieron los resultados que se muestran en la Tabla 5. En general, se obtuvieron 25 patrones de corte diferentes, con un nivel de aprovechamiento promedio del 95%. Los niveles de aprovechamiento más bajos de 91 a 95% se encuentran en las dos primeras clases diamétricas, con marcas de clase de 14.9 y 19.9 cm, respectivamente.

DISCUSIÓN

El estudio presenta una metodología para la optimización del trozado de fustes como herramienta básica de apoyo a las decisiones en el manejo y aprovechamiento de plantaciones forestales (Cancino, 1993). En el proceso de optimización, el modelo determina el patrón de trozado óptimo, señala la posición de los cortes en el fuste y los productos a obtener, además, cubica y clasifica las trozas obtenidas de cada producto. Este tipo de modelo optimizador de trozado ha sido implementado en otros países como Brasil (Arce, 2000), Chile (Garzón 2009) y Argentina (Yapura et al., 2014) para la cuantificación y valoración de volúmenes comerciales y el mejoramiento de la planificación de la cosecha forestal. La disponibilidad de un modelo para la optimización del trozado de fustes permite que cada productor pueda analizar para su situación cuál será el modo más conveniente de segregar su producción maderera (Yapura et al., 2014). El algoritmo desarrollado para la generación de puntos potenciales de corte permitió obtener una disminución considerable (29%) en el número de puntos potenciales a lo largo del fuste. Este resultado está de acuerdo con los reportados por Arce et al. (2004) para Pinus sp. donde se presentan disminuciones hasta de 98% en el número de puntos potenciales de corte.

Al maximizar el ingreso por fuste trozado se asume la utilización óptima del rodal y el aumento de los ingresos económicos. Garzón (2009) indica que la segregación de productos tiene un efecto sobre la recuperación de volumen y el valor del rodal. El autor encontró para Pinus radiata que al incrementar de un solo producto a tres no hubo una diferencia notoria en el valor del árbol y el rodal. Sin embargo, al llegar a siete productos se produjo un aumento significativo en el valor bruto del rodal. En el presente estudio no se observaron diferencias entre el empleo de un algoritmo de trozado basado en programación dinámica y un modelo de trozado secuencial, el cual según Corvalán (1992) simula las decisiones que debe tomar un operario al momento del trozado de un fuste. La falta de diferencias entre los dos métodos para los datos del caso práctico se debe presumiblemente por la marcada diferencia en dimensiones y en precios de los productos evaluados y al reducido número de productos analizados, lo que es similar a lo encontrado por Garzón (2009). Corvalán (1992) reporta que la formulación de las prioridades en modelos de trozado secuencial tiene mucha importancia, en las que se encuentran diferencias si se ordenan los productos por precios estrictamente o si, por el contrario, se ordenan solo por longitud de la troza.

En este estudio, los niveles de aprovechamiento a partir de los patrones de corte óptimos oscilaron entre 91 a 99%, siendo las clases diamétricas más bajas las que presentaron los menores rendimientos en volumen. Esto se debe a que los niveles de aprovechamiento presentan una relación directa con el diámetro índice de utilización e inversa con respecto a la conicidad de las trozas (Cancino, 1993). Similarmente, estudios realizados en Argentina para Eucalyptus grandis han demostrado que el rendimiento está relacionado positivamente con el diámetro índice de utilización y la forma de las trozas (Mastrandrea & Alberti, 2009).

Los resultados encontrados en el presente estudio concuerdan con los obtenidos por Arce (2000) en plantaciones de Pinus sp. en Argentina, quien registró rendimientos de alrededor del 92%. Es importante aclarar que estos resultados no incluyen desperdicios derivados de la corteza, tocones y defectos presentes en las características anatómicas del fuste. De acuerdo con Lledó (2004), el porcentaje asociado a estos desperdicios se encuentra en el orden del 23% (fuste 6%, corteza 11%, tocón 6%), entonces los índices reales de aprovechamiento deben oscilar entre el 63 y 70% de acuerdo con los diferentes patrones de corte obtenidos por el modelo desarrollado. Según Juacida (2004), el porcentaje de aprovechamiento en árboles de Eucalyptus nitens de 12 y 21 años osciló entre 57 y 56%, respectivamente.

La metodología descrita en el presente estudio puede ser adaptada a otras especies comerciales en el país, lo cual requiere el desarrollo de funciones de ahusamiento para las especies de interés (Methol, 2003) que permitan predecir diámetros y alturas fustales y volúmenes comerciales (Prodan et al., 1997). A partir de la metodología propuesta se puede generar una estrategia de aplicación mediante la creación de tablas de esquemas de trozado que serían distribuidas a las cuadrillas de operarios para mejorar las operaciones de trozado en la cosecha forestal (Corvalán, 1992). Sin embargo, estas deben ser actualizadas cada vez que se requieran nuevos productos demandados por el mercado.

CONCLUSIONES

La utilización de técnicas de investigación de operaciones en los procesos de cosecha forestal constituye una gran herramienta para la toma de decisiones en la planificación del trozado y su articulación con la oferta del bosque y la demanda del mercado, la cual ofrece una base sólida para el análisis de procesos industriales en el sector forestal colombiano. A pesar de que estas técnicas han sido reportadas en la literatura desde la década de los setenta, en Colombia no han sido muy utilizadas en el proceso de la cuantificación de volúmenes comerciales para la planificación forestal.

En la presente investigación, el modelo de optimización desarrollado permitió determinar el patrón de trozado óptimo para fustes de E. tereticornis, e identificar la posición de los cortes en el fuste y los productos a obtener que maximizan el valor económico del árbol. El algoritmo propuesto para la generación de puntos potenciales de corte permitió obtener una disminución del 29% en el número de puntos potenciales a lo largo del fuste, lo cual se traduce en mayor eficiencia en el proceso de optimización.

En condiciones donde el número de productos demandados no sea considerable y estos sean fácilmente diferenciables de acuerdo con sus dimensiones y precios, un algoritmo de trozado secuencial o de prioridad generaría los mismos resultados que un algoritmo de optimización de trozado, con la ventaja de su facilidad de entendimiento y menor dificultad en la implementación informática.

La metodología de optimización de trozado desarrollada en el presente estudio puede ser fácilmente adaptable a otras especies. Para su implementación se requiere de funciones de ahusamiento para la especie de interés que permitan predecir diámetros, alturas fustales y volúmenes comerciales. En próximos estudios se debe considerar la incorporación de variables de calidad del fuste para permitir mejorar la segregación de rollizos y hacer descuentos por defectos en el fuste, considerando que estas situaciones son comunes en la cosecha forestal.

AGRADECIMIENTOS

Esta investigación se llevó a cabo bajo el programa de investigación “Fortalecimiento del sector forestal productivo de la Costa Atlántica” derivado del convenio CIAT – Ministerio de Agricultura y Desarrollo Rural (MADR) No. 054 de 2008 y el Programa Cierre de Brechas 2013, financiado por Colciencias, ejecutados por la Corporación Nacional de Investigación y Fomento Forestal (CONIF).

REFERENCIAS BIBLIOGRÁFICAS

Arce, J. (2000). Um sistema de análise, simulação e otimização do sortimento florestal em função da demanda por multi productos e dos custos de transporte (Tese Doutorem Ciências Florestais). Curitiba: Universidade Federal do Paraná. 125 p.
Arce, J., Carnieri, C., Sanquetta, C., & Figueiredo, A. (2002). A forest-level bucking optimization system that considers customers demand and transportation costs. Forest Science, 48(3), 492-503.
Arce, J., MacDonagh, P., & Friedl, R. (2004). Generação de padrões ótimos de corte através de algoritmos de traçamento aplicados a fustes individuais. árvore, 28(2), 207-217.
Briceño, M. (2005). Revisión bibliográfica de modelación del trozado (trabajo de pregrado, Ingeniería Forestal). Talca: Universidad de Talca, Facultad de Ciencias Forestales. 59 p.
Bruce, D., Curtis, R. O. & Vancoevering, C. (1968). Development of a system of taper and volume tables for red alder. Forest Science, 14, 339-350.
Cancino, J. (1993). Modelo optimizador de trozado de árboles basado en funciones de ahusamiento y programación dinámica (Tesis de Maestría).Turrialba: Centro Agronómico Tropical de Investigación y Enseñanza (CATIE), Ciencias Agrícolas y Recursos Naturales. 255 p.
Corvalán, P. (1992). Un modelo optimizador de corta y trozado en Eucalyptus globulus. En: Modelos forestales. Taller internacional / Ipinza, R., Peters, R., Loewe, V. 1992., 266-284.
Eng, G. & Daellenbach, H. (1985). Forest outturn optimization by Dantzig-Wolfe decomposition and dynamic programming column generation. Operations Research, 33(2), 459-464.
Faaland, B. & Briggs, D. (1984). Log bucking and lumber manufacturing using dynamic programing. Management Science, 30(2), 245-257.
Garzón, S. (2009). Efecto del uso de distintos esquemas en la simulación de trozado por prioridad de un rodal podado de Pinus radiata (tesis pregrado, Ingeniería Forestal). Valdivia: Universidad Austral de Chile, Facultad de Ciencias Forestales. 50 p.
Gilmore, P. & Gomory, R. (1961). A linear programming approach to the cutting stock problem. Operations Research, 9(6), 849-859.
Hillier, F. & Lieberman, G. (1997). Introducción a la investigación de operaciones. Atlampa: Editorial McGraw Hill. 1214 p.
Jackson, N. & Smith, G. (1961). Linear programming in lumber production. Forest Products Journal, 39, 272-274.
Juacida, R. (2004). Eucalyptus nitens en Chile: procesos industriales de la madera. Valdivia: Instituto Forestal, Informe Técnico N°164. 74 p.
Kivinen, V.P. (2007). Design and testing of stand-specific bucking instructions for use on modern cut-to-lenght harvesters (Ph.D. Thesis). Department of Forest Resource Management, Faculty of Agriculture and Forestry, University of Helsinky. 65 p.
Laroze, A. (1999). A linear programming, tabu search method for solving forest-level bucking optimization problems. Forest Science, 45(1), 108-116.
Lledó, G. (2004). Evaluación del trozado para rodales de pino insigne en canchas de forestal Bio-Bio S.A. Facultad de Ciencias Forestales. Santiago, Universidad de Santiago de Chile. 74 p.
López, A.M., Barrios, A., Nieto, V. & Trincado, G. (2011). Monitoreo y modelamiento del crecimiento para el manejo de plantaciones forestales comerciales. Corporación Nacional de Investigación y Fomento Forestal (CONIF) -Ministerio de Agricultura y Desarrollo Rural. Bogotá D.C., Colombia. 90 p.
Mastrandrea, C., & Alberti, S. (2009). Características de las trozas de Eucalyptus grandis (W. Hill ex Maiden) y su relación con el rendimiento de aserrado. En: XIII Congreso Forestal Mundial 2009. Buenos Aires, Argentina. 3 p.
Marshall, H.D., Murphy, G. & Boston, K. (2006). Three mathematical models for bucking-to-order. Silva Fennica, 40(1), 127-142.
Methol, R. (2003). SAG grandis: sistema de apoyo a la gestión de plantaciones de Eucalyptus grandis. Montevideo: INIA, Serie Técnica N°131. 42 p.
Pnevmaticos, S. & Mann, S. (1972). Dynamic programming in tree bucking. Forest Products Journal, 22(2), 26-30.
Prodan, M., Peters, R., Cox, F., & Real, P. (1997). Mensura forestal. San José: IICA-GTZ. 561 p.
Sessions, J., Olsen, E., & Garland, J. (1989). Tree bucking for optimal stand value with log allocation constrains. Forest Science, 35(1), 271-276.
Smith, G. & Harrell, C. (1961). Linear programming in log production. Forest Products Journal, 37, 8-11.
Yapura, P.F., Fassola, H.E., Crechi, E.H., Keller, A.E., Sañudo, G., Caraballo, H., Gonzalez, C. & Altamirano, R. (2014). Optimización del trozado de fustes de Pinus taeda, “pino híbrido” (Pinus elliottii x Pinus caribaea F2) y Eucalyptus grandis, en las provincias de Misiones y noreste de Corrientes. Misiones: INTA, Informe final, PIA 10107. 42 p.