Este trabalho teve como objetivo desenvolver modelos de crescimento e produção em classes de diâmetro para plantios de Eucalyptus grandis submetidos a desbastes. Os dados utilizados na modelagem são provenientes de plantios de Eucalyptus grandis localizados no município de Telêmaco Borba, estado do Paraná, disponibilizados pela empresa Klabin S.A. e são oriundos de parcelas permanentes, parcelas temporárias e inventário pré-corte, com idades variando de 2,5 a 26,5 anos. A função de densidade de probabilidade Weibull foi escolhida para gerar a distribuição diamétrica. O processo de modelagem de atributos do povoamento foi dividido em etapas, considerando os estágios de desenvolvimento do povoamento. Dessa maneira, dois sistemas de prognose foram produzidos: o primeiro foi utilizado para a predição de atributos na idade do primeiro desbaste e o segundo para a projeção dos atributos na idade do segundo desbaste ou corte final. Os atributos modelados foram aqueles requisitados para recuperação dos parâmetros da função Weibull nas idades de intervenção pelo Método dos Momentos. Duas alternativas de simulação de desbastes foram testadas. Para avaliação das prognoses foram utilizados o teste de aderência de Kolmogorov-Smirnov e o teste de identidade L&O proposto por Leite e Oliveira (2002). Os dois sistemas de prognose apresentaram, de um modo geral, boas estimativas de ajuste e precisão. A equação de desbaste utilizada por Scolforo (1990) foi escolhida para simulação dos dois desbastes, no entanto, foi necessário a utilização de um fator de correção para melhor controle do número de árvores removido, isto propiciou melhores resultados quando comparados às simulações sem esse controle. A metodologia empregada possibilitou simular o crescimento e a produção em plantios de Eucalyptus grandis ao longo da rotação gerando resultados adequados, principalmente para o corte final onde são colhidas as melhores árvores com maiores valores agregados.
This study aimed to develop growth and yield models by diameter classes for Eucalyptus grandis plantation with thinnings. The data used for modeling came from Eucalyptus grandis plantation made available by the company Klabin S.A., located in Telêmaco Borba, Parana State. This data came from permanent plots, temporary plots and pre-harvesting inventory, with ages ranging from 2.5 to 26.5 years. The Weibull probability density function was chosen to generate the diameter distribution. The stand attributes modeling process was divided into stages, considering the development stages of stands. Thus, two prognosis systems were produced: the first was used for the prediction of the attributes in the age of thinning and the second for the projection of the attributes in the age of the second thinning or final cut. The modeled attributes were those required for recovery the Weibull function parameters at the ages of intervention by the Moments Method. Two alternative methods of thinning simulation were tested. To evaluate the prognoses the Kolmogorov-Smirnov adherence test and the L&O identity test proposed by Leite and Oliveira (2002) were used. The two prognosis systems presented, in general, good estimates of adjustment and accuracy. The thinning equation used by Scolforo (1990) was chosen for the simulation of the two thinnings, however, it was necessary to use a correction factor to improve control of the number of trees removed, generated adequately results compared to simulations without this control. The methodology used allowed to simulate the growth and yield in Eucalyptus grandis stands along the rotation, generating adequate results, especially for the final cut where the best trees are harvested with higher aggregate values.