Estimativa da radiação global incidente em superfícies inclinadas por modelos isotrópicos e índice de claridade

Estimation of the incident global radiation on tilted surfaces using isotropic models and clearness index

AGROMETEOROLOGIA

Estimativa da radiação global incidente em superfícies inclinadas por modelos isotrópicos e índice de claridade

Estimation of the incident global radiation on tilted surfaces using isotropic models and clearness index

Adilson Pacheco de Souza^I; João Francisco Escobedo^I; Alexandre Dal Pai^II; Eduardo Nardini Gomes^I

^IUniversidade Estadual Paulista, Faculdade de Ciências Agronômicas, Departamento de Recursos Naturais, Fazenda Experimental Lageado, Rua José Barbosa de Barros, nº 1.780, CEP 18610-307 Botucatu, SP. E-mail: [e-mail], [e-mail], [e-mail]

^IIFaculdade de Tecnologia de Botucatu, Departamento de Informática, Avenida José Ítalo Bacchi, s/nº, Jardim Aeroporto, CEP 18606-855 Botucatu, SP. E-mail: [e-mail]

RESUMO

O objetivo deste trabalho foi avaliar o desempenho de modelos isotrópicos de estimativa do total de radiação incidente em superfícies inclinadas e propor estimativas com base nas correlações entre os índices de claridade horizontais e inclinados, em diferentes condições de cobertura de céu, em Botucatu, SP. Foram avaliadas superfícies com inclinação de 12,85º, 22,85º e 32,85º, pelos modelos isotrópicos propostos por Liu & Jordan, Revfeim, Jimenez & Castro, Koronakis, a teoria Circunsolar, e a correlação entre os índices de claridade horizontais e inclinados, para diferentes condições de cobertura de céu. O banco de dados de radiação global utilizado corresponde ao período de 1998 a 2007, com intervalos de 4/1998 a 8/2001 para a inclinação de 22,85º, de 9/2001 a 2/2003 para 12,85º e de 1/2004 a 12/2007 para 32,85º. O desempenho dos modelos foi avaliado pelos indicadores estatísticos erro absoluto médio, raiz quadrada do quadrado médio do erro e índice "d" de Wilmott. Os modelos de Liu & Jordan, Koronakis e de Revfeim apresentaram os melhores desempenhos em dias nublados, em todas as inclinações. As coberturas de céu parcialmente difuso e parcialmente aberto, nos maiores ângulos de inclinação, apresentaram as maiores dispersões entre valores estimados e medidos, independentemente do modelo. As equações estatísticas apresentaram bons resultados em aplicações com agrupamentos de dados mensais.

Termos para indexação: ângulo de inclinação, indicativos estatísticos, isotropia, radiação solar, transmissividade atmosférica.

ABSTRACT

The objective of this work was to evaluate the performance of isotropic models estimative of the global radiation on tilted surfaces and to propose estimations based on correlation between the clearness index for horizontal and tilted surfaces, for different sky conditions, in Botucatu, SP, Brazil. The isotropic model proposed by Liu & Jordan, Revfeim, Jimenez & Castro, Koronakis, the Circunsolar theory and the correlation between the clearness index for horizontal and tilted surfaces, for different sky conditions, evaluated tilted surfaces at 12.85º, 22.85º, and 32.85º. Global irradiation database used corresponds to the period from 1998 to 2007, with intervals from 4/1998 to 8/2001 for the surface tilted at 22.85º, from 9/2001 to 2/2003 for the surface tilted at 12.85º and from 1/2004 to 12/2007 for the surface tilted at 32.85º. The performance of the models was analyzed using the statistical indicators mean bias error, root mean square error and d-Wilmott. The models of Liu & Jordan, Koronakis and Revfeim showed the best performance on cloudy days considering all the tilted surfaces. The sky coverage partially diffuse and partially clear with the largest tilted angles showed the highest scattering between estimated and measured values, irrespective of the model. Statistical equations showed good results in monthly average data applications.

Index terms: angle of inclination, indicative statistical, isotropy, solar irradiation, atmospheric transmission.

Introdução

A radiação solar vem se destacando nos últimos anos em pesquisas que objetivaram explorar o rendimento potencial das culturas (Lindquist et al., 2005; Amarante et al., 2007; Kunz et al., 2007; Souza et al., 2009).

Atualmente, no Brasil, encontram-se em pleno desenvolvimento tecnológico as fontes renováveis de energia, como conversão fototérmica e fotovoltaica, biomassa e biodisel, entre outras. Por conseguinte, existe uma crescente demanda por informações sobre a disponibilidade da radiação solar e a sua sazonalidade.

Os monitoramentos e medições da radiação solar ainda são limitados, e restritos principalmente às universidades, com número reduzido de estações de superfície que registram a radiação solar, em que normalmente são realizadas medidas de apenas uma das componentes da radiação solar (radiação global), o que inviabiliza o estudo da sazonalidade das radiações, mesmo para uma única região (Soares et al., 2004; Souza et al., 2005; Tiba et al., 2005; Codato et al., 2008; Borges et al., 2010).

Estudos das intensidades e variações sazonais da radiação incidente em superfícies inclinadas por meio de séries temporais necessitam da realização de medidas por um período mínimo de cinco anos, ou então de estimativas com elevada qualidade e confiabilidade, visto que, se a radiação solar disponível em certa localidade é menor que a estimada, o desempenho dos sistemas será deficiente e as metas econômicas não serão alcançadas (Scolar et al., 2003).

Por causa da falta de medidas em superfícies inclinadas, vários modelos para sua estimativa foram desenvolvidos. Muitos desses modelos utilizam medidas das radiações global, difusa e direta em superfícies horizontais. Normalmente, nessas estimativas são empregados dois tipos de modelos: os estatísticos, ou de decomposição, e os paramétricos. Os modelos paramétricos necessitam de informações mais detalhadas das condições atmosféricas locais, conforme a componente da radiação a ser estimada, ou então a utilização de percentuais fixos de ocorrência de cada componente da radiação global (Liu & Jordan, 1960; Notton et al., 2006; Turco & Rizzatti, 2006; Jiang, 2009; El-Sebaii et al., 2010). Os modelos estatísticos são representados por equações matemáticas ajustadas por meio de regressão linear ou polinomial, em correlações com as frações radiométricas.

A maioria dos modelos propostos para estimativa da radiação solar foram elaborados levando-se em consideração o comportamento isotrópico da radiação (Liu & Jordan, 1960; Revfeim, 1978; Jimenez & Castro, 1982; Koronakis, 1986), entre outros, com enfoques principalmente na estimativa da componente difusa (Diez-Mediavilla et al., 2005; Notton et al., 2006; Loutzenhiser et al., 2007; Noorian et al., 2008), em função do seu comportamento anisotrópico em algumas condições de nebulosidade. No Brasil, os estudos mais recentes sobre o comportamento da radiação global em superfícies inclinadas e suas estimativas foram apresentados por Scolar et al. (2003) e Turco & Rizzatti (2006).

O objetivo deste trabalho foi avaliar o desempenho de modelos isotrópicos de estimativa do total de radiação incidente em superfícies inclinadas a 12,85º, 22,85º e 32,85º, e propor estimativas com base nas correlações entre os índices de claridade horizontais e inclinados, em diferentes condições de cobertura de céu, em Botucatu, SP.

Material e Métodos

Os dados de radiância solar global, direta e difusa, referentes ao período de abril de 1998 a dezembro de 2007 foram coletados na Estação de Radiometria Solar da Universidade Estadual Paulista (Unesp), em Botucatu, SP (22º51'S, 48º26"W, altitude 786 m). De acordo com a classificação climática de Köppen, o clima de Botucatu é Cwa, caracterizado como temperado úmido, com inverno seco e verão quente (Centro de Pesquisas Meteorológicas e Climáticas Aplicadas à Agricultura, 2009).

As leituras foram feitas a cada segundo e as médias armazenadas a cada 5 minutos, em um sistema de aquisição de dados Microlloger, modelo CR23X (Campbell Scientific, Inc., Utah, USA), na frequência de 1 Hz. Na transferência dos dados, foi empregado o módulo de memória SM192, com interface SC532 (Campbell Scientific, Inc., Utah, USA) para microcomputadores, operado pelo software PC 208W (Campbell Scientific, Inc., Utah, USA).

As radiâncias globais no plano horizontal 〈I_GH〉 foram medidas com auxílio de piranômetros Eppley-PSP (The Eppley Laboratory, Rhode Island, USA), com fator de calibração de 7,45 µV W^-1 m^-2 e linearidade de ±0,5% (0 a 2.800 W m^-2). Para a radiância global inclinada 〈I_Gb〉, os piranômetros foram posicionados paralelamente sobre planos inclinados de 12,45º 〈|L| - 10º〉, 22,45º 〈|L|〉 e 32,45º 〈|L| + 10º〉, com a face voltada para o norte. As radiâncias diretas foram obtidas por um pireliômetro Eppley-NIP, acoplado a um rastreador solar Eppley, modelo ST3 (The Eppley Laboratory, Rhode Island, USA), com fator de calibração de 7,59 µV W^-1 m^-2 e linearidade de ±0,5% (0 a 1.400 W m^-2). A projeção da radiância direta na horizontal foi dada pelo produto entre a radiância direta e o cosseno do ângulo zenital. A radiância difusa na horizontal foi obtida pela diferença entre a radiação global e direta.

Os valores das radiações diária global horizontal e inclinada , direta horizontal e difusa horizontal foram obtidos pela integração das respectivas curvas de radiância, entre o nascer e o pôr do sol (Iqbal, 1983).

Na inclinação de 12,45º, as medidas foram realizadas entre 9/2001 e 2/2003; na inclinação de 22,85º, entre 4/1998 e 8/2001, e na inclinação de 32,45º, entre 1/2004 e 12/2007. Esses dados foram submetidos a uma análise de consistência, e os valores discrepantes oriundos de erros de leitura dos sensores ou mal funcionamento do sistema de aquisição de dados foram removidos das séries.

O desempenho dos seguintes modelos para estimativa da radiação global que incide em superfície inclinada com ângulo (β) a partir de medidas na horizontal, foi analisado: Liu & Jordan (1960), Revfeim (1978), Koronakis (1986), Jimenez & Castro (1982), e Circunsolar (Iqbal, 1983).

A radiação diária global inclinada segundo o modelo Circunsolar foi calculada pela equação (Iqbal, 1983):

em que: R_B é um fator geométrico decorrente da mudança no ângulo de incidência dos raios solares sobre a superfície inclinada, quando comparada com a superfície horizontal. Segundo Iqbal (1983), na partição diária, esse fator pode ser obtido diretamente por,

R_B = [(π/180) ω'_s sin δ sin (φ - β) + cos δ cos (φ - β) sin ω'_s]/[ (π/180) ω_s sin δ sin φ + cos δ cos φ sin ω_s],

em que: ω_s é o ângulo horário solar diário horizontal, calculado pela equação:

e ω'_s é o ângulo horário solar diário inclinado, calculado pela equação:

ω'_s = min{cos^-1 (-tan φ tan δ) cos^-1 [-tan δ tan ( φ - β )]},

em que: φ é a latitude local; β é o ângulo de inclinação; δ é a declinação solar. Esta última é calculada pela equação:

,

em que DJ é o número do dia ao longo do ano (dia juliano). A radiação diária global inclinada pelo modelo de Liu & Jordan (1960) foi calculada pela equação:

em que ρ é o albedo da superfície.

A estimativa das radiações diária global inclinada pelo modelo de Koronakis (1986) foi calculada pela equação:

Essa estimativa se assemelha à de Liu & Jordan (1960), todavia assume que a componente isotrópica da radiação difusa que atinge o plano inclinado equivale a 66,7% do total de radiação difusa da superfície horizontal.

O modelo de Jimenez & Castro (1982), que estima a radiação diária global inclinada pela equação:

também emprega a base teórica de Liu & Jordan (1960), porém com a pressuposição de que a radiação difusa em um plano horizontal corresponde a 20% da radiação global.

O modelo desenvolvido por Revfeim (1978) emprega a fração radiométrica difusa e as limitações do plano de horizonte pela inclinação de uma superfície:

em que: K_D é a razão entre a radiação difusa e a global na horizontal; f_β é a proporção do plano de horizonte bloqueada pela superfície inclinada, quando comparada com a superfície horizontal, estimada por f_β = (1 - β/180). Essa fração K_D também pode ser estimada apenas com base em medidas da radiação global na horizontal (Scolar et al., 2003).

A radiação incidente no topo da atmosfera para uma superfície horizontal pode ser estimada pela seguinte equação:

em que H_SC é conhecida como a constante solar, com valor igual a 4.921 kJ m^-2 h^-1, E₀ é o fator de correção da excentricidade da órbita da Terra. Segundo Iqbal (1983), o fator de excentricidade pode ser calculado por:

Além desses modelos de estimativa, foram propostas neste trabalho correlações entre as frações da radiação global inclinada . Os índices são obtidos pela razão entre a radiação global e a radiação no topo da atmosfera para superfícies inclinadas:

Nas correlações entre os valores estimados pelos modelos de estimativa da radiação global inclinada supracitados e os valores medidos de , foram usados dados de janeiro a dezembro de 2002, para a inclinação de 12,85º, de janeiro de 1999 a dezembro de 2000, para 22,85º, e de janeiro de 2004 a dezembro de 2005 para 32,85º. Foram empregados, nas análises de desempenho das estimativas, dados dos seguintes períodos: de abril de 1998 a dezembro de 1998 e de janeiro a agosto de 2001, para a inclinação 22,85º; de setembro a dezembro de 2001 e de janeiro a fevereiro de 2003, para a inclinação 12,85º; e de janeiro a dezembro de 2006, para a inclinação 32,85º.

Na avaliação do desempenho das equações de estimativa diária dos seis modelos e das correlações propostas, foram empregados os indicadores estatísticos erro absoluto médio ("Mean bias error" MBE), raiz quadrada do quadrado médio do erro ("Root mean square error" RMSE), recomendados por Iqbal (1983), e o índice de ajustamento "d" de Willmott, para diferentes condições de coberturas de céu (Diez-Mediavilla et al., 2005; Notton et al., 2006; Loutzenhiser et al., 2007; Escobedo et al., 2009; Evseev & Kudish, 2009; El-Sebaii et al., 2010). Esses indicadores foram calculados pelas seguintes equações:

em que: Pi representa os valores estimados, Oi os valores medidos, N é o número de observações, |P'i| é o valor absoluto da diferença Pi - , e |O'i| representa o valor absoluto da diferença Oi - .

A classificação da cobertura de céu em função do índice de claridade empregada neste trabalho seguiu a proposta de Escobedo et al. (2009), dada pela seguinte classificação de cobertura de céu para a partição diária: <0,35, a radiação direta é praticamente nula, sendo a radiação global igual à radiação difusa (céu nublado); 0,35<

>0,55, a radiação global é composta principalmente pela fração difusa e em menor escala pela radiação direta (céu parcialmente difuso); 0,55<

<0,65, é o inverso do segundo intervalo, visto que a radiação difusa diminui gradativamente (céu parcialmente aberto); >0,65, a radiação global possui a maior parcela composta pela radiação direta e mínima difusa, portanto, a cobertura do céu é considerada aberta. Para valores de acima de 0,78, não ocorre uma sequência lógica de evolução das três radiações em função da transmissividade atmosférica e, por essa razão, foram eliminados das análises e também dos modelos de estimativas gerados.

Resultados e Discussão

Em relação à frequência de distribuição dos valores de em diferentes condições de cobertura de céu, observou-se que, no município de Botucatu, SP, foi verificado no máximo 70 dias de céu nublado – valores de inferiores a 0,35 – em 2004, o que corresponde a 19,18% dos dias ao longo do ano (Figura 1). Em torno de 65,48, 58,08, 61,37, 52,52 e 55,07% dos dias dos anos de 1999, 2000, 2002, 2004 e 2005, respectivamente, apresentaram índices de claridade acima de 0,55, indicando que na maioria dos dias a atmosfera local transmite acima de 55% da radiação incidente no topo da atmosfera para o plano horizontal. Escobedo et al. (2009) verificaram, nesse mesmo município, que, entre janeiro de 2001 e dezembro de 2004, ocorreram 217 dias nublados, 215 dias com céu parcialmente difuso, 271 dias com céu parcialmente aberto e 490 dias com céu aberto.

Os modelos isotrópicos Liu & Jordan (1960), Revfeim (1978), Jimenez & Castro (1982), Koronakis (1986) apresentaram ótimas correlações entre os valores medidos e estimados de radiação global incidente, nas três inclinações estudadas (Figura 2). Os coeficientes de determinação nas correlações desses modelos foram em torno de 0,98, com tendência de diminuição dos valores de R² com o incremento do ângulo de inclinação, como pode ser evidenciado para o modelo Jimenez & Castro (1982), que apresentou valores de R² de 0,99; 0,98 e 0,97, para as inclinações de 12,85º, 22,85º e 32,85º, respectivamente.

A estimativa de com base na teoria circunsolar apresentou correlações semelhantes aos modelos isotrópicos, para pequenas inclinações (Figura 2). Todavia, com o aumento do ângulo de inclinação, a dispersão gráfica dos dados aumentou significativamente, possibilitando menores valores do coeficiente de determinação.

Nos menores ângulos de inclinação (12,85º), não ocorreu variação no comportamento dos modelos isotrópicos nas diferentes coberturas de céu (Figura 2). Porém, em ângulos maiores (32,85º), as maiores dispersões foram originadas em dias com céu parcialmente aberto (0,55 <>0,65) e aberto ( > 0,65). Para estimativas com o modelo Circunsolar, as maiores dispersões ocorreram em dias com cobertura de céu parcialmente nublado, tanto difuso quanto aberto.

Em dias com céu aberto, a disponibilidade de densidade de fluxo de radiação solar é maior, normalmente superior a 15 MJ m^-2 dia^-1, decorrente do aumento das componentes direta e refletida, em contrapartida da diminuição da radiação difusa. As maiores dispersões para céu aberto podem ser explicadas pelo aumento da componente refletida, decorrentes do termo ou seja, é grandemente dependente dos maiores níveis de . Essas dispersões ainda são decorrentes da forma com que a radiação difusa é originada nesses dias (Iqbal, 1983; Gueymard, 2009). Scolar et al. (2003) verificaram que modelos isotrópicos de estimativa de apresentam correções da componente difusa com bons resultados, em dias com céu nublado, pois desconsideram o comportamento anisotrópico de dias com céu aberto.

Contudo, o modelo Circunsolar apresentou menores dispersões em dias de céu aberto do que em dias de céu nublado ou parcialmente difuso, em consequência de a maior percentagem da radiação difusa em dias de céu aberto ser proveniente do disco solar, ou seja, a componente circunsolar (Iqbal, 1983; Scolar et al., 2003).

Em relação aos desempenhos dos cinco modelos isotrópicos de estimativa, nos agrupamentos em diferentes condições de cobertura de céu e com todos os dados, observou-se que, independentemente das classes de , todos os modelos apresentaram tendência de subestimar os valores de para a inclinação de 12,85º (Figura 3).

Para os agrupamentos com todos os dados, em todas as inclinações estudadas, o modelo Revfeim (1978) subestimou os valores de . Ainda nesse agrupamento de dados, as maiores superestimativas foram obtidas pelo modelo Circunsolar em 32,85º (0,347 MJ m^-2 dia^-1), que corresponderam em termos percentuais de desvios entre os valores médios mensais medidos e estimados em -5,15 e 1,93%, respectivamente. Nas diferentes condições de cobertura de céu, o modelo Revfeim (1978) superestimou os valores de apenas em céu aberto, para a inclinação de 32,85º (0,329 MJ m^-2 dia^-1 ou 1,37%). Todos os modelos apresentaram tendências de superestimativas e subestimativas similares, em todas as inclinações e coberturas de céu, que variaram entre -0,621 MJ m^-2 dia^-1 ou -3,78%, no modelo Revfeim (1978), em 〈|L|+10º〉 para céu parcialmente difuso, e 0,510 MJ m^-2 dia^-1 ou 6,61%, no modelo Jimenez & Castro (1982), em 〈|L|+10º〉 para céu nublado.

As maiores dispersões entre valores estimados e medidos de foram verificadas na aplicação do modelo Circunsolar, sendo crescente com o aumento do índice de claridade e equivalentes a 0,762, 1,606 e 2,012 MJ m^-2 dia^-1 para as inclinações de 12,85, 22,85 e 32,85º, em dias com céu parcialmente difuso, respectivamente. Nos cinco modelos isotrópicos estudados, os índices de ajustamento apresentaram diminuição com o aumento da inclinação e diminuição da nebulosidade. Com exceção do modelo Circunsolar, os ajustamentos permaneceram superiores a 90% em todas as condições de cobertura de céu. Essas variações nos valores de MBE, RMSE e "d" são verificadas na maioria dos trabalhos na partição diária da radiação global (Diez-Mediavilla et al., 2005; Notton et al., 2006; Loutzenhiser et al., 2007; Escobedo et al., 2009; Evseev & Kudish, 2009).

Quanto às estimativas de , com base nas correlações entre e , observou-se que as medidas de nos ângulos de 12,85 e 32,85º foram maiores do que as verificadas no plano horizontal, em função dos coeficientes angulares serem superiores a 1 (Figura 4). Na inclinação de 22,85º, ocorreu uma inversão, com indicação de que a radiação global na horizontal é maior do que na inclinada. Nota-se uma linearidade entre os valores de e , visto que para pequenos valores de são encontrados pequenos valores de .

No desempenho das estimativas da radiação global, em superfícies inclinadas com base na radiação global na horizontal, foram encontrados coeficientes de determinação (R²) acima de 0,88, que em conjunto com os demais índices estatísticos permitiram boas estimativas. Essas correlações apresentaram tendências a superestimar os valores de , com MBE de 1,98, 5,36 e 6,93%, e índices de ajustamento em torno de 0,99. Todavia, quando os modelos são validados em função das frações estimadas, não são possíveis análises em termos de energia. As três equações ajustadas apresentaram tendência de superestimar os valores de , com valores de MBE crescentes com o aumento do ângulo de inclinação. Em 12,85, 22,85 e 32,85º , o MBE foi de 0,394 MJ m^-2 (1,93%), 1,025 MJ m^-2 (5,30%) e 1,485 MJ m^-2 (7,86%), respectivamente. As dispersões entre valores estimados e medidos de também seguiram essa tendência de aumento com a intensidade de inclinação, e variaram de 1,186 MJ m^-2 (5,81%) a 3,733 MJ m^-2 (19,76%). Consequentemente, os índices de ajustamento foram reduzidos em função das inclinações, sendo de 0,9901, 0,95227 e 0,87918, para 12,85, 22,85 e 32,85º, respectivamente. As maiores dispersões verificadas com o aumento da inclinação foram decorrentes do aumento dos níveis da radiação refletida em superfícies inclinadas com maiores ângulos e também do ganho de radiação direta nas maiores inclinações, no período de inverno.

Os valores médios mensais medidos e estimados de correspondem aos mesmos períodos de medidas na avaliação do desempenho das estimativas com correlações entre e . Observou-se que, para a estimativa na inclinação de 12,85º nos meses de janeiro, fevereiro, novembro e dezembro, todos os modelos subestimaram os valores de , sendo os valores medidos superiores em até 24,14%, em novembro, para o modelo Liu & Jordan (1960) (Figura 5). Ainda para essa inclinação, o modelo Jimenez & Castro (1982) apresentou melhores resultados, em estimativas no mês janeiro (-2,66%), fevereiro (-0,94%) e novembro (-0,68%), enquanto o modelo Revfeim (1978) foi melhor em outubro (+0,67%) e dezembro (-0,98%), e o modelo Circunsolar foi superior em setembro (+0,02 MJ m^-2 dia^-1).

Para a inclinação de 22,85º, os modelos Liu & Jordan (1960) e Revfeim (1978) subestimaram, enquanto o modelo Revfeim (1978) superestimou, em todos os meses, os valores de (Figura 5). Os modelos gerados pelas correlações entre os índices de claridade horizontal e inclinados e Circunsolar apresentaram tendências de superestimar os valores de no verão e subestimar nos meses de inverno. As maiores subestimativas foram de 27,49%, ou 5,42 MJ m^-2 dia^-1, em fevereiro, para o modelo Liu & Jordan (1960), e as maiores superestimavas foram propiciadas pelo método das correlações entre osíndices de claridade horizontal e inclinados no mês de junho (+9,33%, ou 1,67 MJ m^-2 dia^-1). Por conseguinte, para estimativas na inclinação de 32,85º, os modelos Koronakis (1986) e Liu & Jordan (1960) apresentaram comportamentos semelhantes, superestimando os valores de nos meses de novembro e dezembro em, no máximo, 0,91 MJ m^-2 dia^-1. Os piores resultados foram verificados pelos modelos das correlações entre os índices de claridade horizontal e inclinados, nos meses de inverno, com superestimativas de até 23,09% (4,48 MJ m^-2 dia^-1) no mês de julho.

De Miguel et al. (1995) verificaram que os modelos Koronakis (1986) e Liu & Jordan (1960) apresentaram resultados similares na estimativa de em Burgos, na Espanha, região localizada na latitude 42º20'25", para superfícies inclinadas até 40º, com face voltada para o sul. Entre os modelos Koronakis (1986) e Jimenez & Castro (1982), esses mesmos autores verificaram similaridades nas estimativas apenas nos meses de verão. Tian et al. (2001) concluíram que o modelo Revfeim (1978) estimou satisfatoriamente os valores de , considerando as variações de cobertura de céu para as cidades de Kataia e Invercargill (latitudes 35º e 46º S), na Nova Zelândia.

Kamali et al. (2006) estudaram oito modelos de estimativa de para superfícies inclinadas a 45º e 40º, com faces voltadas para o sul e leste, respectivamente, em Karaj, Irã (latitude 35º55'N) e verificaram que o modelo Koronakis (1986) apresentou os melhores desempenhos na predição da componente difusa para a face leste.

Para melhorar o desempenho das equações geradas pelas correlações entre e , foram propostas equações geradas com agrupamentos mensais (Tabela 1), dadas por regressão linear simples, como . Neste caso, existem dois períodos distintos que acompanham as evoluções das radiações global e extraterrestre. No verão, as superfícies horizontais apresentam níveis incidentes maiores que os verificados na inclinada, com coeficientes angulares maiores que 1,0 e, entre os equinócios, essa situação se inverte. As piores correlações foram verificadas no mês de abril, dadas pelo início da entrada das frentes frias na região, o que propicia uma grande variação no número de dias com céu nublado e céu aberto. Nos meses de verão, dezembro e janeiro, as equações para as inclinações de 22,85 e 32,85º apresentaram tendências de subestimar, dadas pelos valores de MBE negativos (Tabela 2).

Na inclinação de 32,85º, foram verificados os maiores desvios das médias nas estimativas de , com equações mensais de inverno apresentando variação entre 1,30 e 1,90 MJ m^-2 entre os meses de maio e agosto. Acima de 85% dos modelos mensais superestimam os valores de . Para menores períodos no agrupamento de dados, foram verificados melhores ajustamentos em todas as inclinações. Todavia, nos meses em que ocorre maior número de dias com céu parcialmente nublado, os modelos mensais apresentaram os piores desempenhos, para os três indicativos analisados.

Conclusões

1. Quanto maior a transmissividade atmosférica do dia, maior a dispersão entre valores medidos e estimados de radiação global inclinada por modelos isotrópicos.

2. Menores ângulos de inclinação propiciam melhores desempenhos de modelos isotrópicos na estimativa da radiação global inclinada.

3. Os modelos de Liu & Jordan, Jimenez & Castro e Koronakis podem ser aplicados independentemente da condição de cobertura do céu.

4. Em agrupamentos mensais, as equações de estimativa com base nas correlações entre os índices de claridade horizontal e inclinados possibilitam bons resultados.

Recebido em 1 de julho de 2009 e aprovado em 12 de março de 2010

• AMARANTE, C.V.T. do; STEFFENS, C.A.; MOTA, C.S.; SANTOS, H.P. dos. Radiação, fotossíntese, rendimento e qualidade de frutos em macieiras 'Royal Gala' cobertas com telas antigranizo. Pesquisa Agropecuária Brasileira, v.42, p.925-931, 2007.
• BORGES, V.P.; OLIVEIRA, A.S. de; COELHO FILHO, M.A.; SILVA, T.S.M. da; PAMPONET, B.M. Avaliação de modelos de estimativa da radiação solar incidente em Cruz das Almas, Bahia. Revista Brasileira de Engenharia Agrícola e Ambiental, v.14, p.74-80, 2010.
• CENTRO DE PESQUISAS METEOROLÓGICAS E CLIMÁTICAS APLICADAS A AGRICULTURA. Clima dos municípios paulistas Disponível em: <http://www.cpa.unicamp.br>. Acesso em: 30 jun. 2009.
• CODATO, G.; OLIVEIRA, A.P.; SOARES, J.; ESCOBEDO, J.F.; GOMES, E.N.; PAI, A.D. Global and diffuse solar irradiances in urban and rural areas in southeast Brazil. Theoretical and Applied Climatology, v.93, p.57-73, 2008.
• DE MIGUEL, A.; BILBAO, J.; DIEZ-MEDIAVILLA, M. Solar radiation incident on tilted surfaces in Burgos, Spain: isotropic models. Energy Conversion and Management, v.36, p.945-951, 1995.
• DIEZ-MEDIAVILLA, M.; DE MIGUEL, A.; BILBAO, J. Measurement and comparison of diffuse solar irradiance models on inclined surfaces in Valladolid (Spain). Energy Conversion and Management, v.46, p.2075-2092, 2005.
• EL-SEBAII, A.A.; AL-HAZMI, F.S.; AL-GHAMDI, A.A.; YAGHMOUR, S.J. Global, direct and diffuse solar radiation on horizontal and tilted surfaces in Jeddah, Saudi Arabia. Applied Energy, v.87, p.568-576, 2010.
• ESCOBEDO, J.F.; GOMES, E.N.; OLIVEIRA, A.P.; SORAES, J. Modeling hourly and daily fractions of UV, PAR and NIR to global solar radiation under various sky conditions at Botucatu, Brazil. Applied Energy, v.86, p.299-309, 2009.
• EVSEEV, E.G.; KUDISH, A.I. The assessment of different models to predict the global solar radiation on a surface tilted to the south. Solar Energy, v.83, p.377-388, 2009.
• GUEYMARD, C.A. Direct and indirect uncertainties in the prediction of tilted irradiance for solar engineering applications. Solar Energy, v.83, p.432-444, 2009.
• IQBAL, M. An introduction to solar radiation Toronto: Academic Press, 1983. 389p.
• JIANG, Y. Estimation of monthly mean daily diffuse radiation in China. Applied Energy, v.86, p.1458-1464, 2009.
• JIMENEZ, J.I.; CASTRO, Y. Solar radiation on sloping surfaces with different orientations in Granada, Spain. Solar Energy, v.28, p.257-262, 1982.
• KAMALI, G.A.; MORADI, I.; KHALILI, A. Estimating solar radiation on tilted surfaces with various orientations: a study case in Karaj (Iran). Theoretical and Applied Climatology, v.84, p.235-241, 2006.
• KORONAKIS, P.S. On the choice of the angle of tilt for south facing solar collectors in the Athens basin area. Solar Energy, v.36, p.217-225, 1986.
• KUNZ, J.H.; BERGONCI, J.I.; BERGAMASCHI, H.; DALMAGO, G.A.; HECKLER, B.M.M.; COMIRAN, F. Uso da radiação solar pelo milho sob diferentes preparos do solo, espaçamento e disponibilidade hídrica. Pesquisa Agropecuária Brasileira, v.42, p.1511-1520, 2007.
• LINDQUIST, J.L.; ARKEBAUER, T.J.; WALTERS, D.T.; CASSMAN, K.G.; DOBERMANN, A. Maize radiation use efficiency under optimal growth conditions. Agronomy Journal, v.97, p.72-78, 2005.
• LIU, B.Y.H.; JORDAN, R.C. The interrelationship and characteristic distribution of direct, diffuse and total solar radiation. Solar Energy, v.4, p.1-19, 1960.
• LOUTZENHISER, P.G.; MANZ, H.; FELSMANN, C.; STRACHAN, P.A.; FRANK, T.; MAXWELL, G.M. Empirical validation of models to compute solar irradiance on inclined surfaces for building energy simulation. Solar Energy, v.81, p.254-267, 2007.
• NOORIAN, A.M.; MORADI, I.; KAMALI, G.A. Evaluation of 12 models to estimate hourly diffuse irradiation on inclined surfaces. Renewable Energy, v.33, p.1406-1412, 2008.
• NOTTON, G.; CRISTOFARI, C.; POGGI, P. Performance evaluation of various hourly and daily slope irradiation models using Mediterranean experimental data of Ajaccio. Energy, Conversion and Management, v.47, p.147-173, 2006.
• REVFEIM, K.J.A. A simple procedure for estimating global daily radiation on any surface. Journal of Applied Meteorology, v.17, p.1126-1131, 1978.
• SCOLAR, J.; MARTINS, D.; ESCOBEDO, J.F. Estimativa da irradiação total sobre uma superfície inclinada a partir da irradiação global na horizontal. Revista Brasileira de Geofísica, v.21, p.249-258, 2003.
• SOARES, J.; OLIVEIRA, A.P.; BOZNAR, M.Z.; MLAKAR, P.; ESCOBEDO, J.F.; MACHADO, A.J. Modeling hourly diffuse solar-radiation in the city of São Paulo using a neural-network technique. Applied Energy, v.79, p.201-204, 2004.
• SOUZA, J.L. de; NICÁCIO, R.M.; MOURA, M.A.L. Global solar radiation measurements in Maceio, Brazil. Renewable Energy, v.30, 1203-1220, 2005.
• SOUZA, P.J. de O.P. de; RIBEIRO, A.; ROCHA, E.J.P. da; FARIAS, J.R.B.; LOUREIRO, R.S.; BISPO, C.C.; SAMPAIO, L. Solar radiation use efficiency by soybean under field conditions in the Amazon region. Pesquisa Agropecuária Brasileira, v.44, p.1211-1218, 2009.
• TIAN, Y.Q.; DAVIES-COLLEY, R.J.; GONG, P.; THORROLD, B.W. Estimating solar radiation on slopes of arbitrary aspect. Agricultural and Forest Meteorology, v.109, p.67-74, 2001.
• TIBA, C.; AGUIAR, R. de; FRAIDENRAICH, N. Analysis of a new relationship between monthly global irradiation and sunshine hours from a database of Brazil. Renewable Energy, v.30, p.957-966, 2005.
• TURCO, J.E.P.; RIZZATTI, G.S. Avaliação de modelo matemático para estimar a radiação solar incidente sobre superfícies com diferentes exposições e declividades. Engenharia Agrícola, v.26, p.257-267, 2006.

Datas de Publicação

• Publicação nesta coleção
  12 Jul 2010
• Data do Fascículo
  Abr 2010

Histórico

• Aceito
  12 Mar 2010
• Recebido
  01 Jul 2009