Geração de padrões ótimos de corte através de algoritmos de traçamento aplicados a fustes individuais

Optimal bucking pattern generation through cutting algorithms applied to individual stems

Geração de padrões ótimos de corte através de algoritmos de traçamento aplicados a fustes individuais

Optimal bucking pattern generation through cutting algorithms applied to individual stems

Julio Eduardo Arce^I; Patricio MacDonagh^II; Ramón Alejandro Friedl^III

^IDepartamento de Ciências Florestais, Universidade Federal do Paraná, Av. Pref. Lothário Meissner, 3400, 80210-170 Curitiba-PR, Tel.: [phone] - Fax: [phone], <[e-mail]>

^IICátedra de Aprovechamiento Forestal, Fac.de Ciencias Forestales, Universidad Nacional de Misiones, Bertoni 124, Caixa Postal 3380, Eldorado, Misiones, Argentina, Tel.: [phone], <[e-mail]>

^IIICátedra de Dasometría, Fac. de Ciencias Forestales, Universidad Nacional de Misiones, <[e-mail]>

RESUMO

No presente trabalho descreve-se uma metodologia de otimização do corte ou traçamento de árvores para obtenção de multiprodutos em nível de fustes individuais. A geração dos padrões ótimos de corte foi realizada através de dois módulos: 1) um algoritmo heurístico que identifica os cortes potenciais (números úteis), que devem ser avaliados ao longo do fuste; e 2) um algoritmo com base na programação dinâmica (PD), que determina a combinação ótima de toras de cada tipo de comprimento a serem retiradas do fuste, visando maximizar o seu comprimento utilizado. Um exemplo da obtenção dos padrões ótimos de corte para um conjunto de 25 fustes é apresentado. O traçamento ótimo e o traçamento real obtido pelo motosserrista foram comparados. O resíduo de madeira deixado na floresta foi reduzido de 16 para 5% através dos padrões de corte ótimos.

Palavras-chave: Otimização, programação dinâmica, heurísticas e Pinus.

ABSTRACT

This work describes a method of optimizing individual tree stem bucking operations. The individual optimal cutting patterns were generated by means of two modules: 1) a heuristic algorithm which identifies the potential cuts (useful numbers) to be evaluated along the stem, and 2) a Dynamic Programming based-algorithm to determine the optimum combination of products to be obtained to maximize stem profit. An example is provided to illustrate the optimum cutting patterns for 25 individual stems. Optimal and conventional bucking patterns were compared. Residues were reduced from about 16% to 5% by means of the optimal bucking patterns.

Key words: Optimization, dynamic programming, heuristics and Pinus.

1. INTRODUÇÃO

A colheita representa a operação final de um ciclo de produção florestal, na qual são obtidos os produtos mais valiosos, constituindo um dos fatores que determinam a rentabilidade florestal. A maneira como é conduzido o traçamento das árvores obtidas na colheita ou corte raso é fundamental para estabelecer as quantidades de cada um dos diferentes tipos de produtos que podem ser obtidos da floresta. Porém, esta tarefa é executada quase exclusivamente pelo motosserrista, com base na sua intuição, o que pode comprometer a rentabilidade do negócio florestal.

O traçamento das árvores derrubadas nas operações de colheita florestal pode ser considerado como um processo de tomada de decisões de curto prazo. De fato, a posição dos cortes que definirão as toras deve ser decidida no local, após a visualização da árvore abatida, e esta tarefa muitas vezes é realizada sem o suporte de ferramentas de otimização. Quando um único comprimento de tora é retirado da floresta estas decisões tornam-se bem mais simples, porém no caso de multiprodutos florestais a complexidade combinatória necessária para atingir o traçamento otimizado aumenta. A utilização de técnicas de otimização aplicada à obtenção de multiprodutos florestais pode contribuir de forma significativa para aumentar a utilização da matéria-prima, reduzindo os resíduos deixados na floresta ou subutilizados, ao serem aproveitados como energia.

Em algumas indústrias, como as de vidro, metalúrgica e têxtil, o planejamento da produção depara-se freqüentemente com os denominados problemas de corte. Nestas indústrias, a redução dos custos de produção é freqüentemente obtida através da seguinte estratégia: a matéria-prima utilizada é inicialmente produzida em tamanhos grandes padronizados, possivelmente estocada e, somente mais tarde, reduzida a tamanhos menores, para então ser utilizada pela indústria, ou para atender às demandas externas de tamanhos variados, muitas vezes não-padronizados (Arenales & Morabito, 1997). Esta estratégia gera a necessidade de um estágio adicional de produção, denominado corte ou desdobramento da matéria-prima. As perdas decorrentes destas operações de corte são inevitáveis, motivando o seu planejamento para minimizar os efeitos negativos gerados pelo desperdício sobre os custos de produção.

Sob essa mesma ótica, as florestas comerciais podem ser vistas como fábricas produtoras de grandes peças de matéria-prima – árvores – de dimensões conhecidas – distribuições diamétricas e hipsométricas – e com uma grande capacidade de estoque. De fato, por ser a árvore fábrica e produto ao mesmo tempo (Gregersen & Contreras, 1980), a totalidade da produção da floresta vai sendo estocada nas árvores até o momento da colheita, quando elas devem ser derrubadas e cortadas em peças menores – toras –, atendendo aos diversos pedidos realizados pelos clientes.

Alguns sistemas de inventário florestal e prognose do crescimento e da produção existentes no mercado possuem módulos de sortimento que informam a quantidade de cada tipo de produto que o talhão possui. Mas para tal fim, estes sistemas consideram somente a otimização da floresta em pé – otimização em função da oferta –, uma vez que não dispõem de cadastros de pedidos ou demandas mínimas por tipo de produto. Dentre estes sistemas podem ser mencionados o sistema de inventário MARVL (Method for assessment of recoverable volume by log types), desenvolvido na Nova Zelândia a partir do trabalho de Deadman & Goulding (1978), e o SISPINUS (Oliveira, 1995), desenvolvido especificamente pela EMBRAPA-CNPF (Empresa Brasileira de Pesquisa Agropecuária – Centro Nacional de Pesquisas de Florestas) para plantios de Pinus spp. no Sul do Brasil, os quais não são alimentados com informações de demanda.

No trabalho pioneiro de Borges (1981) foi desenvolvido, para Pinus taeda L., no Estado do Paraná, Brasil, um programa de computação para o seccionamento do fuste, visando obter, de cada árvore, o número máximo de peças serradas com as maiores dimensões. Para tal fim, foram construídas funções de afilamento-padrão relativas por classes do fator de forma natural, baseadas na função polinomial de quinto grau. A partir dos mesmos dados utilizados pelo autor, foram utilizadas curvas de afilamento em séries de quocientes naturais transformados. O objetivo foi estimar volumes parciais através da integração, entre os pontos relativos desejados, das equações que descrevem estas séries (Jorge, 1982).

Dentre os trabalhos desenvolvidos no Brasil, foi apresentado em 1998, por Andrade & Leite, um método para quantificar multiprodutos de árvores individuais na unidade estere. Em Assis et al. (2001), encontra-se uma comparação entre modelos polinomiais segmentados e não-segmentados na estimativa de diâmetros e volumes para Pinus taeda. Foram comparados neste trabalho o modelo segmentado de Clark et al., o modelo segmentado de Max e Burkhart, o polinômio de potências fracionárias e inteiras de Hradetzky e o modelo de Goulding e Murray.

No trabalho de Faaland & Briggs (1984), foi utilizada a PD para resolver o problema do corte ótimo das árvores em produtos serrados terminados. Este problema se caracteriza por duas etapas: o desdobramento da árvore em toras e a subseqüente transformação destas toras em produtos como tábuas, sarrafos e demais peças obtidas no desdobro. Estas duas atividades foram consideradas como um sistema único de produção.

Uma abordagem de subalgoritmos hierárquicos de PD foi utilizada por Reinders & Hendriks (1989), para conversão ótima de árvores em peças serradas. Esta conversão foi descrita em três níveis: 1) árvores em toras; 2) toras em peças – fatias –, obtidas por cortes paralelos realizados nas toras; e 3) peças menores, obtidas por cortes paralelos entre si e perpendiculares aos cortes anteriores.

Maness & Adams (1991) utilizaram um modelo estruturado como um problema de PL. O modelo utilizou estratégias de traçamento de árvores em toras e de toras em peças serradas, através de um procedimento de resolução de três estágios. O volume dos produtos elaborados foi controlado através de relações preço-volume, simulando curvas de demanda.

O sistema OPTICORT, desenvolvido pela equipe chilena de pesquisadores liderada pelo professor Weintraub, para um grupo de empresas florestais do Chile (Epstein et al., 1999), faz parte de uma verdadeira malha de sistemas de Pesquisa Operacional aplicados ao planejamento e desenvolvimento do setor florestal. A malha está formada por um sistema de planejamento florestal de longo prazo (PLANEX), um sistema de otimização do sortimento (OPTICORT) e um sistema de programação do transporte florestal (ASICAM). Em particular, o OPTICORT consiste em um módulo de PL combinado com um algoritmo de geração de colunas. Este algoritmo, na verdade, é um gerador de padrões de corte tipo branch-and-bound, acrescentando, a partir da ramificação e avaliação iterativa de uma árvore de alternativas, opções às soluções ótimas encontradas pelo módulo de PL.

O objetivo deste trabalho foi descrever uma metodologia para otimizar o traçamento de fustes de árvores florestais, baseada em um algoritmo heurístico e na programação dinâmica, que se valora através da quantificação das perdas de madeira para serraria ou pasta celulósica, obtidas em campo pelos métodos tradicionais.

2. MATERIAL E MÉTODOS

Para avaliar a eficácia da metodologia proposta, foi rigorosamente quantificada uma operação real de colheita e traçamento, utilizando 25 fustes de Pinus spp., selecionados de forma aleatória. Os dados foram obtidos na província de Misiones, Argentina. Para cada fuste, foi obtido o comprimento útil da porção utilizável para serraria, considerando um diâmetro comercial de 18 cm. No Quadro 1 estão as características dos fustes utilizados.

Após a medição do comprimento útil de todos os fustes selecionados, foi executado o traçamento por um motosserrista idôneo, sem a presença da equipe de medição, segundo o critério-padrão, ou seja, da maneira como habitualmente é feito. É importante frisar que o operário desconhecia a avaliação que estava sendo realizada. Os comprimentos de tora aceitos pela serraria foram l₁ = 12' (pés), l₂ = 14' e l₃ = 16'. Entretanto, na prática, ao converter estes comprimentos ao sistema métrico decimal e incluir a folga para a destopa na serraria, estes três comprimentos passam a ser l₁ = 3,75 m, l₂ = 4,35 m, l₃ = 4,99 m.

Após o traçamento foram medidas, em cada um dos fustes, todas as toras que o motosserrista obteve, caracterizando, desta maneira, a situação ou o cenário real de traçamento.

A situação otimizada de traçamento foi obtida através de padrões de corte preestabelecidos. Para tanto, inicialmente foram definidos alguns volumes e resíduos em nível de fustes individuais, conforme mostrado na Figura 1.

A geração dos padrões de corte ótimos foi efetivada através de dois módulos: 1) um algoritmo heurístico que identifica os cortes potenciais (números úteis), que devem ser avaliados ao longo do fuste, e 2) um algoritmo com base na programação dinâmica (PD) que determina a combinação ótima de produtos a serem retirados do fuste, com o objetivo de maximizar o seu valor (Gilmore & Gomory, 1961, 1963, 1965; Carnieri et al., 1994a, 1994b).

2.1. Algoritmo Heurístico para Geração dos Números Úteis (Algoritmo 1)

O algoritmo heurístico para geração dos números úteis foi adaptado a partir da proposta de Carnieri et al. (1994a), em seu trabalho de otimização do corte bidimensional de painéis de madeira compensada. Na seqüência descreve-se o algoritmo proposto para geração do conjunto útil de pontos onde serão avaliadas as estratégias de corte ao longo do fuste de comprimento L, considerando m produtos de comprimentos diferentes l₁, l₂,..., l_m.

Inicialmente é obtido o menor comprimento b a partir da lista de diferentes produtos; é definido também o primeiro número útil como sendo o início do comprimento útil (Passo 1). Na seqüência, a partir de cada número útil encontrado, devem ser calculados os pontos de corte ao longo do fuste para cada novo produto k extraindo j toras (j = 1, 2, ...), desde que a porção remanescente do fuste seja suficientemente longa para retirar pelo menos uma tora do produto de menor comprimento b (Passo 2). Os números úteis encontrados desta forma – F^(k) – são acrescentados na lista existente – U^{(k + 1)} –, removendo-se quaisquer números repetidos ou equivalentes (Passo 3). Finalmente, o procedimento se repete até considerar o último produto (k = 1), concluindo a geração dos números úteis (Passo 4), os quais são armazenados no vetor U ⁽¹⁾. A Figura 2 apresenta o pseudocódigo deste algoritmo.

2.2. Algoritmo de PD para Geração do Padrão de Corte Ótimo (Algoritmo 2)

A programação dinâmica (PD), também denominada programação em múltiplos estágios, caracteriza-se pela avaliação sucessiva da evolução da otimização em estágios através de uma equação recursiva. Em cada um destes estágios são comparados vários estados ou alternativas, e é guardado somente o valor – label – do melhor destes. Desta maneira, conforme a proposta de Gilmore & Gomory (1961) e Carnieri et al. (1994a), os diferentes produtos que são considerados na otimização de um determinado fuste ou classe diamétrica representam as etapas do problema de PD, e os números úteis gerados através do algoritmo constituem os estados. A seguinte equação recursiva tem por função maximizar o valor do fuste.

em que x = ponto do fuste no qual é realizada a avaliação; F_s(x) = valor da melhor combinação de produtos obtida até o comprimento x, utilizando somente os primeiros s produtos; l_s = comprimento do produto s; P_s = valor da tora do produto s que está sendo considerado; e F_s(x – l_s) = valor da melhor combinação de produtos obtida até o comprimento (x – l_s), utilizando somente os primeiros s produtos.

2.3. Padrões Ótimos de Corte

A partir dos dois algoritmos descritos nas seções anteriores foram gerados os padrões ótimos de corte para os fustes analisados. Para cada fuste foram gerados os números úteis, que indicam os cortes potenciais que poderiam ser realizados, através do algoritmo 1. A partir destes números úteis, foi posteriormente gerado o padrão de corte ótimo, através do algoritmo 2.

2.4. Comparação Entre as Operações de Traçamento Real e Otimizado

Finalmente, foram comparados o traçamento real e o traçamento otimizado, como forma de validar a metodologia proposta. Para realizar esta comparação, foi utilizada a diferença agregada, de modo a quantificar a diferença total ou cumulativa entre as operações de traçamento real e otimizado.

3. RESULTADOS E DISCUSSÃO

3.1. Geração dos Números Úteis

A geração dos números úteis é ilustrada através do seguinte exemplo. Seja um fuste de comprimento útil L = 18,35 m até um diâmetro comercial de 18 cm. Os m = 3 produtos diferentes a serem obtidos devem possuir, respectivamente, os seguintes comprimentos: l₁ = 3,75 m, l₂ = 4,35 m, l₃ = 4,99 m. Ao realizar uma busca exaustiva, considerando, por exemplo, o máximo divisor comum dos três comprimentos de tora, que neste caso é 1 cm, surgem (18,35 m/0,01 m) = 1.835 pontos (0,01; 0,02; 0,03; ...;18,33; 18,34; e 18,35). Através do algoritmo proposto foram obtidos 19 pontos, o que representa uma redução considerável no número de pontos a avaliar. Estes pontos, que representam possíveis cortes ao longo do fuste do exemplo, são 3,75; 4,35; 4,99; 7,50; 8,10; 8,70; 8,74; 9,34; 9,98; 11,25; 11,85; 12,45; 12,49; 13,05; 13,09; 13,69; 13,73; 14,33; e 18,35. O último valor (18,35) não é obtido a partir do pseudocódigo apresentado; seu ingresso é forçado para poder quantificar, posteriormente, as perdas de aproveitamento no traçamento do fuste.

Durante a geração dos números úteis não é estabelecida a ordem em que as diferentes toras poderiam ser obtidas do fuste; este ordenamento é realizado posteriormente, na geração dos padrões de corte.

3.2. Obtenção do Padrão Ótimo de Corte Através do Algoritmo de PD

No Quadro 2 está o procedimento de solução através do algoritmo de PD proposto, para o mesmo fuste de 18,35 m utilizado como exemplo na geração dos números úteis.

Em um problema de PD a solução ótima geralmente deve ser reconstruída a partir da equação recursiva. Este procedimento, denominado de reconstrução da trajetória, é indicado no Quadro 2 através das células hachuradas. O valor ótimo encontrado aparece na célula inferior direita da última linha do Quadro 2, correspondente ao último estado (x = 18,35) no último estágio (s = 3). O valor desta célula (18,08 m) indica que traçando o fuste de forma ótima são utilizados 18,08 m. O resíduo associado a esta solução ótima é de 0,27 m, obtido através da diferença entre o comprimento útil (18,35 m) e o comprimento efetivamente utilizado (18,08 m). A partir deste vértice começa a reconstrução da trajetória. O último produto retirado na tora que gera o valor ótimo aparece sobrescrito na célula. Neste caso, o número 3 indica que a última tora retirada é do produto 3, com um comprimento l₃ de 4,99 m. Portanto, a célula que segue na reconstrução da trajetória é aquela localizada no estado x = 13,09 m, subtraindo o comprimento da tora (l₃ = 4,99 m) do valor ótimo obtido no último estado (18,08 m), ou seja, analisando a porção remanescente do fuste após a retirada de uma tora de 4,99 m. Este comprimento remanescente corresponde ao estado x = 13,09 m.

A célula que segue na reconstrução da trajetória localiza-se no estágio s = 3, estado x = 13,09. Nesta célula aparece novamente o número 3 (hachurado), indicando que a penúltima tora do padrão ótimo de corte também é do produto 3, cujo comprimento é de 4,99 m.

Na seqüência, prosseguindo com a reconstrução da trajetória, é obtida a célula correspondente ao estado x = 8,10 m, subtraindo o comprimento da tora anterior (4,99 m) ao estado anterior (13,09 m - 4,99 m = 8,10 m). Neste caso, a célula hachurada com o número sobrescrito 2, correspondente ao estado x = 8,10 m e ao estágio s = 3, indica que a antepenúltima tora a ser retirada é do produto 2, cujo comprimento é de 4,35 m. Deste modo, mantendo o estado x = 8,10 m, a célula hachurada se desloca para a coluna da esquerda, correspondente ao estágio s = 2.

Finalmente, a partir da célula anterior é subtraído o comprimento da tora anterior (l₂ = 4,35) ao estado anterior x = 8,10 m. O valor obtido de 3,75 m indica o estado no qual deve se localizar a tora seguinte no padrão de corte ótimo. Observa-se no Quadro 2 que a célula hachurada correspondente ao estado x = 3,75 m no estágio s = 2 indica que a tora seguinte é do produto 1 (l₁ = 3,75 m).

Em síntese, a partir da reconstrução da trajetória é obtido o padrão de corte ótimo para o fuste utilizado como exemplo, que consiste em retirar uma tora de 3,75 m, uma tora de 4,35 m e duas toras de 4,99 m, totalizando 3,75 + 4,35 + 2 x 4,99 = 18,08 m.

3.3. Implementação dos Padrões Ótimos de Corte Através de Regras de Traçamento

A implementação prática dos resultados deste trabalho pode ser feita, por exemplo, através da elaboração de regras de traçamento, as quais seriam distribuídas de forma impressa e plastificadas entre os motosserristas. Estas regras deveriam ser atualizadas cada vez que algum novo sortimento com diâmetro e, ou, comprimento diferente for demandado por parte dos compradores das toras. O Quadro 3 contém os padrões de corte ótimos, gerados a partir dos algoritmos 1 e 2, considerando três comprimentos diferentes de tora (l₁ = 3,75 m; l₂ = 4,35 m; l₃ = 4,99 m), para fustes com comprimentos úteis teóricos de 4,00 a 24,00 m, com intervalos de 50 cm. Somente por uma questão de espaço não se apresentam os padrões de corte para fustes com intervalos menores, por exemplo 10, 5 e 1 cm. Porém, para validar a metodologia proposta, os padrões de corte foram gerados para fustes com intervalos de 1 cm no que se refere ao comprimento útil. Estes padrões de corte constam na base de dados do programa utilizado para geração dos padrões ótimos de corte.

No exemplo considerado, foi analisado um único tipo de produto com seu diâmetro comercial (d), o que significa que, analogamente, apenas um único comprimento útil foi considerado em cada fuste. Para o caso de múltiplos d's, simplesmente deveriam ser utilizados os algoritmos 1 e 2 em forma repetitiva, tantas vezes quanto diferentes d devam ser avaliados, começando pelas toras de maior d e finalizando com as de menor d. Ao comprimento útil do fuste correspondente a cada d deve ser acrescentado o comprimento útil residual das toras com d imediatamente anterior. Esta situação não foi considerada no exemplo apresentado.

No exemplo foram tratados os diferentes comprimentos de tora como se todos os produtos tivessem o mesmo preço unitário (R$/m³). No caso de se querer utilizar preços diferenciados por produto, além dos algoritmos descritos, deveriam ser utilizadas funções de afilamento^1/ para o cálculo dos volumes parciais correspondentes. Embora existam modelos de afilamento que permitem estimar diretamente as alturas em diferentes diâmetros comerciais, como é o caso das equações apresentadas em Demaershalk & Kozak (1977), em se tratando dos polinômios de grau elevado deve ser utilizado algum algoritmo específico, como o método da bisseção, para obter as alturas nas quais se localizam os diferentes d no fuste (Arce et al., 1999; Arce, 2000). É importante ressaltar que existem ainda sistemas de processamento de inventário florestal comercializados no Brasil que requerem do usuário o ajuste prévio de duas funções polinomiais de quinto grau independentes: uma para a estimativa de d_i = f (h_i) e outra para a estimativa de h_i = f (d_i). As curvas que definem o perfil de uma mesma árvore a partir destes dois modelos, embora muito semelhantes, dificilmente apresentarão um ajuste perfeito.

3.4. Comparação entre as Operações de Traçamento Real e Otimizado

No Quadro 4 estão o resultado do traçamento realizado na floresta pelo motosserrista (situação real) e o resultado otimizado em nível de tora. Neste quadro é indicado o comprimento útil dos fustes medidos na floresta, a quantidade de toras obtidas em cada fuste de cada um dos três comprimentos utilizados, o comprimento total utilizado obtido através da soma dos comprimentos das toras e os resíduos de cada um dos fustes analisados.

Observa-se (Quadro 4) que em alguns casos os resíduos, obtidos como a diferença entre o comprimento útil até d igual a 18 cm e o comprimento aproveitado pelo motosserrista, apresentam valores negativos, indicando que foi utilizado um comprimento superior ao indicado pelo d. Esta situação implica o fato de, provavelmente, algumas das toras obtidas possuírem um diâmetro em ponta fina inferior ao estipulado para as toras para serraria (18 cm), podendo eventualmente ser rejeitadas pela indústria. Porém, por tratar-se de uma situação real de colheita, estes casos foram mantidos com a finalidade de caracterizar a situação real da melhor forma possível. Obviamente, padrões ótimos de corte nunca geraram estes resíduos negativos, uma vez que os dois algoritmos utilizados efetuam comparações de comprimentos até o comprimento máximo definido como útil.

No Quadro 4 observa-se que os padrões ótimos de corte não conseguem evitar por completo a geração de determinados resíduos nos fustes. Apenas para um único fuste, cujo comprimento útil foi de 11,85 m, foi atingida a redução teórica completa deste resíduo, através do traçamento do mesmo em duas toras de 3,75 m e uma tora de 4,35 m (2 x 3,75 + 4,35 = 11,85). Porém, o resíduo máximo segundo os padrões de corte otimizados foi de 2,41 m, para um fuste com comprimento útil de 7,40 m. Paradoxalmente, neste caso o motosserrista conseguiu obter uma situação mais favorável, retirando duas toras de 3,75 m em vez de uma tora de 4,99 m. O resíduo negativo obtido pelo motosserrista (-0,18 m) nunca poderia ter sido alcançado utilizando os padrões de corte, uma vez que os algoritmos 1 e 2 baseiam-se em critérios precisos de avaliação de possibilidades potenciais de corte. De qualquer forma, a habilidade do motosserrista em contornar situações difíceis, considerando que alguns fustes possuem defeitos como tortuosidade, bifurcações, dentre outros, é fundamental para atingir, conjuntamente com os padrões ótimos de corte, a maior eficiência nas operações de colheita florestal.

No Quadro 5 apresenta-se uma análise comparativa dos resíduos que efetivamente foram deixados na floresta após as tarefas de traçamento e os que poderiam ter sido deixados no caso de utilizar os padrões ótimos de corte.

Dos resíduos deixados na floresta foram considerados apenas os resíduos úteis, ou seja, aquelas toras cujo d atingiu a exigência da indústria (d > 18 cm), mas por questões de comprimento foram abandonadas.

Observa-se no Quadro 5 que os mesmos comprimentos totais foram considerados para o conjunto de dados reais e otimizados. Isto foi possível devido ao fato de utilizar padrões ótimos de corte gerados, considerando intervalos de 1 cm no comprimento útil dos fustes. Neste sentido, o comprimento útil total, ou seja, a soma dos comprimentos úteis de todos os fustes, foi de 296,50 m.

Analisando o fato de que em determinados fustes o resíduo deixado pelo motosserrista foi menor do que o resíduo sugerido pelos padrões ótimos de corte (fustes 4, 6, 14 e 19), verifica-se que estas situações foram possíveis somente porque foi utilizado um comprimento útil maior que o efetivamente medido na floresta antes das operações de traçamento. O que deve ser aproveitado deste fato é a experiência prática do motosserrista, que em determinadas circunstâncias utiliza algum tipo de "tolerância intuitiva" para executar sua tarefa.

Sem considerar os resíduos negativos, em todos os demais casos os resíduos sugeridos pelos padrões ótimos de corte foram menores que os efetivamente obtidos na floresta. Enquanto o conjunto de dados reais apresentou um resíduo total de 48,55 m, com os padrões ótimos de corte este resíduo foi reduzido a 15,07 m, o que equivale a um aumento no comprimento total utilizado de 247,95 a 281,43 m.

A Figura 3 apresenta o número de casos observados por classe de magnitude do resíduo. Observa-se que os resíduos otimizados distribuem-se predominantemente nas primeiras duas classes, com até 5 e 10% do comprimento útil desperdiçado, respectivamente, enquanto os resíduos reais apresentam números de casos observados mais elevados nas classes de 20, 25, 30, 45 e 50% do comprimento útil desperdiçado.

Em termos relativos, o Quadro 5 mostra que o resíduo deixado na floresta, equivalente a 16% do comprimento útil total, foi reduzido para 5% com a utilização dos padrões ótimos de corte. Estes números ilustram o ganho potencial de utilizar este tipo de procedimento próprio da Pesquisa Operacional nas atividades florestais de traçamento de fustes.

A experiência prática do motosserrista, que o possibilitou utilizar certos níveis de tolerância nos comprimentos e nos diâmetros em ponta fina, deveria ser incorporada de alguma maneira neste tipo de soluções matemáticas, visando aprimorar ainda mais a eficiência das operações florestais de colheita.

4. CONCLUSÃO

A utilização de técnicas de Pesquisa Operacional foi eficaz para obtenção de padrões de corte ótimos para fustes de árvores comerciais. O algoritmo de geração de números úteis reduziu o número potencial de cortes que devem ser avaliados ao longo do comprimento do fuste. O algoritmo de programação dinâmica garante a obtenção de padrões ótimos de corte, considerando os fustes existentes e os comprimentos de tora dos produtos demandados pela indústria.

A utilização dos padrões ótimos de corte permitiu redução nos resíduos reais, equivalentes a 16% do comprimento total utilizado, para 5%, considerando os resíduos teóricos.

5. REFERÊNCIAS BIBLIOGRÁFICAS

ARCE, J. E. Um sistema de análise, simulação e otimização do sortimento florestal em função da demanda por multiprodutos e dos custos de transporte. 2000. 129 f. Tese (Doutorado em Ciências Florestais) – Universidade Federal do Paraná, Curitiba, 2000.

ARCE, J. E. et al. Patrones individuales óptimos de corte para la obtención de multiproductos forestales. In: CONGRESO FORESTAL LATINOAMERICANO, 1., 1999, Lima, Peru. Anais... Lima, Perú: 1999. (CD ROM Abstract).

ARENALES, M. N.; MORABITO, R. O problema de corte e empacotamento e aplicações industriais. In: CONGRESSO NACIONAL DE MATEMÁTICA APLICADA E COMPUTACIONAL, 2., OFICINA NACIONAL DE PCE, 2., 1997, Gramado-RS. Anais... Gramado: 1997. 146 p.

ASSIS, A. L. et al. Comparação de modelos polinomiais segmentados e não-segmentados na estimativa de diâmetros e volumes ao longo do fuste de Pinus taeda. Cerne, v. 7, n. 1, p. 20-40, 2001.

BORGES J. F. Seccionamento do fuste de Pinus taeda L. para obtenção do volume de madeira serrada através da função de forma polinomial. 1981. 81 f. Dissertação (Mestrado em Ciências Florestais) – Universidade Federal do Paraná, Curitiba, 1981.

CARNIERI, C.; MENDOZA, G. A.; GAVINHO, L. G. Solution procedures for cutting lumber into furniture parts. European Journal of Operational Research, v. 73, p. 495-501, 1994a.

CARNIERI, C.; MENDOZA, G. A.; GAVINHO, L. G. Optimal cutting of lumber and particleboards into dimension parts: some algorithms and solution procedures. Wood and Fiber Science, v. 26, n. 1, p. 131-141, 1994b.

DEADMAN, M. W.; GOULDING, C. J. A method for assessment of recoverable volume by log types. New Zeland Forest of Science, v. 9, n. 2, p. 225-239, 1978.

DEMAERSHALK, J. P.; KOZAK A. The whole bole system: a conditioned dual-equation system for precise prediction of tree profiles. Canadian Journal Forest of Research, v. 7, p. 488-497, 1977.

EPSTEIN R. et al. Use of OR Systems in the Chilean Forest Industries. Interfaces, v. 29, n. 1, p. 7-29, 1999.

FAALAND, B.; BRIGGS, D. Log bucking and lumber manufacturing using dynamic programming. Management Science, v. 30, n. 2, p. 245-257, 1984.

GILMORE, P. C.; GOMORY R. E. A linear programming approach to the cutting stock problem. Operations Research, v. 9, p. 848-859, 1961.

GILMORE, P. C.; GOMORY R. E. A linear programming approach to the cutting stock problem, Part II. Operations Research, v. 11, p. 863-888, 1963.

GILMORE, P. C.; GOMORY R. E. Multi-stage stock cutting problems of two or more dimensions. Operations Research, v. 13, p. 94-120, 1965.

GREGERSEN, H.; CONTRERAS A. Análisis económico de proyectos forestales. Roma: FAO Montes, 1980.

JORGE, L. A. B. Utilização de curvas de forma em séries relativas transformadas nas estimativas de volumes totais e parciais de árvores de Pinus taeda L. 1982. 91 f. Dissertação (Mestrado em Ciências Florestais) – Universidade Federal do Paraná, Curitiba, 1982.

MANESS T. C.; ADAMS D. A. The combined optimization of log bucking and sawing strategies. Wood and Fiber Science, v. 23, n. 2, p. 296-314, 1991.

OLIVEIRA, E. B. Um sistema computadorizado de prognose do crescimento e produção de Pinus taeda L., com critérios quantitativos para a avaliação técnica e econômica de regimes de manejo. 1995. 134 f. Tese (Doutorado em Ciências Florestais) – Universidade Federal do Paraná, Curitiba, 1995.

REINDERS, M. P.; HENDRIKS, Th. H. B. Lumber production optimization. European Journal of Operational Research, v. 42, p. 243-253, 1989.

Recebido para publicação em 4.4.2002 e aceito para publicação em 30.4.2004.

• ANDRADE, V. C. L.; LEITE, H. G. Um método para quantificar multiprodutos de árvores individuais na unidade estere. Revista Árvore, v. 22, n. 3, p. 299-306, 1998.

1

As funções de afilamento, funções de forma ou funções

taper, estimam diâmetros a diferentes alturas

d_i = f (h_i)

_.

Datas de Publicação

• Publicação nesta coleção
  10 Ago 2004
• Data do Fascículo
  Abr 2004