Estimativa da irradiação total sobre uma superfície inclinada a partir da irradiação global na horizontal

Estimativa da irradiação total sobre uma superfície inclinada a partir da irradiação global na horizontal

José Scolar^I; Dinival Martins^II; João Francisco Escobedo^III

^IPesquisador aposentado do Instituto de Pesquisas Meteorológicas da UNESP de BAURU - Rua Martim Afonso 11-28 Vila Souto - CEP 17051-040 Bauru, SP, Brasil Fone (14) 32272832 - E-mail: [e-mail]

^IIFaculdade de Ciências Agronômicas – Departamento de Ciências Ambientais da UNESP- Botucatu - SP - Fazenda Experimental Lageado – Cx Postal 237 CEP 18603-970 Botucatu, SP, Brasil. Fone (14) 38117162 – E-mail: [e-mail]

^IIIFaculdade de Ciências Agronômicas – Departamento de Ciências Ambientais da UNESP- Botucatu - SP - Fazenda Experimental Lageado – Cx Postal 237 CEP 18603-970 Botucatu, SP, Brasil. Fone (14) 38117162 - E-mail: [e-mail]

RESUMO

Neste trabalho realizou-se um estudo de modelagem da irradiação total que incide sobre uma superfície inclinada, com face voltada para o Equador e inclinação de 22,85º, igual a da latitude local, usando partição diária de energia em Botucatu – SP. Dos vários modelos de estimativa da irradiação total em uma superfície inclinada existentes na literatura, três deles (LIU; JORDAN, 1963; HAY, 1979; PEREZ et al., 1987) foram utilizados e modificados de duas maneiras diferentes, com o propósito de se obter um modelo que utilize somente a irradiação global na estimativa da irradiação total na superfície inclinada. Na primeira, a irradiação difusa foi estimada em função da irradiação global, através do ajuste empírico utilizando a técnica de regressão linear, na forma K_d=f(K_t) , fração da difusa contida na global (K_d), com o índice de claridade (K_t). Nesse caso, o melhor ajuste é dado por um polinômio de quarto grau. Na segunda modificação, a irradiação direta K_BH, fração da irradiação direta contida na irradiação global, foi ajustada através da técnica de regressão linear, em função do índice de claridade, na forma K_BH=f(K_t ). Nesse caso, o melhor ajuste foi linear, obtido para intervalos específicos de K_t. Os resultados obtidos com essas modificações mostram que os modelos têm melhores desempenhos para a modificação feita com a função K_BH=f(K_t ), e o modelo modificado de Hay apresenta o melhor desempenho na estimativa da irradiação total em Botucatu-SP.

Palavras-chave: Irradiação solar, Superfície inclinada, modelos.

ABSTRACT

A modeling study of the total irradiation reaching a tilted surface facing towards the Equator at a place with the same latitude (-22.85º), using the Daily partitioning energy, is presented. In this study three total irradiation incident on tilted surface estimate models (LIU; JORDAN, 1963; HAY, 1979; PEREZ et al., 1987) were modified in two different manners, to improve irradiation estimative model for tilted surface knowing only the global irradiation. The first modification consisted to apply empirical functions to determine the diffuse irradiation by the global radiation measurements. These empirical functions were tuned using a linear regression technique, the fraction of the diffuse irradiation (K_d) from the global irradiation, was correlated with the clearness index (K_d). In this case the highest correlation found was a fourth order polynomial function. In the second modification the fraction of the direct irradiation K_BH from the global irradiation was correlated assuming a linear correlation at given intervals of K_t. Model results showed that the empirical relationship K_BH=f(K_t ) presented better results and the Hay model presented a better performance compared with the two others models for the Botucatu region.

Keywords: Solar irradiation, Tilted surface, models.

INTRODUÇÃO

A crescente demanda por informações sobre a utilização da radiação solar como fonte suplementar e alternativa às fontes convencionais de energia, requer o conhecimento detalhado das características climáticas da radiação solar incidente sobre superfícies inclinadas em relação à superfície horizontal, para serem usadas em uma grande variedade de aplicações, incluindo projetos de engenharia para coletores solares, projetos de arquitetura, planejamento urbano, estudos agronômicos de insolação sobre vegetação e em estudos micrometeorológico sobre circulação local. Para satisfazer qualquer um desses estudos é necessário conhecer a intensidade da radiação solar que incide sobre uma superfície inclinada e sua variação sazonal por um período mínimo de um ano. Contudo, o alto custo instrumental e a necessidade de pessoal operacional dificultam a implantação e operação de uma rede de estações de radiometria solar. O resultado é um número limitado de estações ao redor do globo terrestre para observações das componentes da radiação solar na superfície horizontal e principalmente a observação da radiação solar em uma superfície inclinada, já que nem todas as estações de radiometria solar fazem essas medidas.

Por causa da falta de observações da radiação solar que incide em uma superfície inclinada, vários modelos teóricos de estimativa foram desenvolvidos, utilizando medidas da radiação global e da componente direta da radiação obtidas na superfície horizontal (LIU; JORDAN, 1963; HAY, 1979; TEMPS; COULSON, 1977; KLUCHER, 1979, PEREZ et al., 1987). Esses modelos utilizam individualmente os valores da radiação direta, difusa e da refletida pela superfície do solo como componentes da radiação total que incide sobre uma superfície inclinada. A única diferença entre os vários modelos existentes está no tratamento da radiação difusa que incide sobre a superfície inclinada, devido ao seu comportamento anisotrópico.

Em muitas localidades a componente direta da radiação solar não é medida, existindo somente observações da radiação global na superfície horizontal. Nessa circunstância, para estimar a radiação total que incide sobre uma superfície inclinada, as componentes difusa e direta da radiação solar devem ser estimadas, e isto pode ser feito usando modelos conceituais que relacionam a fração da radiação difusa contida na radiação global com o índice de claridade, fração entre as irradiações global e do topo da atmosfera, K_d=f(K_t) , ou a fração da direta contida na global com o índice de claridade K_BH=f(K_t ). São exemplos dessas correlações os modelos de Erbs Klein e Duffie (1982); Oliveira e outros (2002); Ricieri (1988), Gueymard (1993); Feuermann e Zemel (1992). Nesse caso, são utilizados para estimar a radiação total na superfície os modelos empíricos de estimativa da radiação difusa ou da direta, desenvolvidos em uma determinada localidade, em conjunto com os valores da radiação solar global medida localmente na superfície horizontal inclinada (KLEIN, 1977; BUGLER, 1977; NOTTON, MUSELLI; LOUCHE, 1996; BEHR, 1997).

O objetivo deste trabalho é obter um modelo para estimar a irradiação total diária que incide em uma superfície inclinada para Botucatu-SP, usando somente a radiação global medida na superfície horizontal, utilizando modelos modificados de estimativa da radiação total incidente na superfície inclinada, com as funções K_d=f(K_t) , e por K_BH=f(K_t ), e verificar o desempenho dos modelos modificados por essas relações empíricas.

CONJUNTO DE DADOS

Os dados das irradiâncias solar direta e global para a superfície horizontal e da total em um plano inclinado igual a latitude local (22º54') e a face voltada para o Equador, foram coletados na Estação de Radiometria Solar da UNESP de Botucatu-SP (latitude 22º51' Sul, longitude 48º26' Oeste, altitude 786m), entre setembro de 1997 e agosto de 2001. As leituras foram feitas a cada segundo e as médias armazenadas a cada cinco minutos. Cerca de 9% dos dados foram desprezados na análise de consistência. Os modelos (LIU; JORDAN, 1963; HAY, 1979; PEREZ et al., 1987) de estimativa da irradiação total na superfície inclinada foram avaliados, e os modelos modificados, desenvolvidos com os dados coletados no período de setembro de 1997 a agosto de 2000 (1034 dias), validados com os dados obtidos no período setembro de 2000 a agosto de 2001 (332 dias).

A irradiância direta na incidência normal foi obtida por meio de um pireliômetro EPPLEY-NIP, acoplado a um rastreador solar EPPLEY modelo ST-3, fixo na direção norte-sul geográfica e ajustado na latitude local. A irradiância global no plano horizontal foi obtida com auxílio de um piranômetro EPPLEY-PSP, posicionado num plano horizontal e, no plano inclinado foi obtida também com auxílio de um piranômetro EPPLEY-PSP, posicionado paralelamente sobre o plano inclinado.

As irradiações diárias (), no topo da atmosfera para a superfície horizontal, no topo da atmosfera para a superfície inclinada, direta projetada na horizontal, global no plano horizontal, difusa na superfície horizontal obtida pela diferença entre as irradiações global e direta na horizontal e total observada na superfície inclinada, foram obtidas através da integração das respectivas curvas de irradiância, integradas do nascer ao pôr-do-sol.

.

MODELOS DE ESTIMATIVA DA IRRADIAÇÃO TOTAL

A irradiação total diária , que incide em uma superfície inclinada, com um ângulo b de inclinação em relação a superfície horizontal, é dada pela soma das irradiações direta , difusa e refletida pelo solo , isto é:

A irradiação direta que atinge a superfície inclinada, segundo Iqbal (1983), é calculada através de correção da irradiação direta projetada na horizontal , dada pela seguinte expressão:

onde, R_B= cos (q_b) / cos (q_z) é o fator de correção devido à mudança no ângulo de incidência dos raios solares sobre a superfície inclinada, q_z é o ângulo zenital do Sol e q_b é o ângulo de incidência dos raios solares sobre um superfície inclinada com a face voltada para o Equador.

Considerando-se que a reflexão pela superfície do solo seja isotrópica, a fração da energia refletida pelo solo que atinge a superfície inclinada é dada por:

onde, r é o albedo da superfície definido como a irradiação refletida por uma superfície sobre a irradiação incidente nessa superfície, b é o ângulo de inclinação e o termo [(1 – cos (b))] é a razão entre a irradiação incidente na superfície inclinada sobre a radiação refletida pelo solo.

Esse tratamento matemático para as irradiações direta e refletida pelo solo, que incidem numa superfície inclinada, é comum à maioria dos modelos, excetuando-se uns poucos, como os de Temps e Coulson (1977) e Gueymard (1987), que dão tratamento anisotrópico para a irradiação refletida pelo solo. A diferença básica entre os vários modelos de estimativa da irradiação total, está na modelagem da irradiação difusa, desde que a sua distribuição espacial e temporal geralmente é desconhecida devido ao seu comportamento anisotrópico, tornando-a uma variável difícil de ser modelada.

Modelos de estimativa da irradiação difusa

Normalmente são utilizadas três subcomponentes para avaliar o comportamento anisotrópico da radiação difusa que incide na superfície inclinada: circunsolar, brilho horizontal e radiação difusa isotrópica. A irradiação difusa que vem da região ao redor do disco solar é chamada de circunsolar, sendo predominantemente resultante do espalhamento devido aos aerossóis existente na atmosfera. O aumento na intensidade da radiação difusa próximo ao horizonte é chamado de brilho horizontal, em virtude de a radiação incidente percorrer uma trajetória maior na atmosfera próximo ao horizonte e também por causa das múltiplas reflexões na atmosfera terrestre. Isotrópica é o restante da radiação difusa uniformemente distribuída na porção restante da hemisfério celeste.

Vários modelos têm sido propostos para estimar a irradiação difusa que atinge uma superfície inclinada; no entanto, nem todos consideram essas três subcomponentes. Desses vários modelos existentes na literatura apenas três deles foram selecionados e utilizados para estimar a irradiação difusa, (LIU; JORDAN, 1963; HAY, 1979; PEREZ et al., 1987).

Liu e Jordan (1963) assumem que a intensidade da irradiação difusa é uniforme em todo o hemisfério celeste. Sendo portanto, válida para condições de céu nublado, condição na qual a irradiação difusa é suposta isotrópica. Seguindo o mesmo procedimento adotado para a irradiação refletida pelo solo, Liu e Jordan (1963) propõem a seguinte equação para o cálculo da irradiação difusa que atinge uma superfície inclinada:

onde o termo [(1 – cos (b))] é a razão entre a irradiação difusa incidente na superfície inclinada sobre a irradiação difusa vinda da cúpula do céu.

Esse modelo apresenta bons resultados para condições de céu nublado, quando a irradiação difusa assume comportamento isotrópico. No entanto, em condições de céu claro, a irradiação difusa tem um comportamento anisotrópico e o modelo subestima os valores calculados. Portanto, para o cálculo da irradiação difusa que atinge uma superfície inclinada, é necessário que os modelos incorporem os efeitos tanto da isotropia quanto da anisotropia da irradiação difusa.

Hay (1979) considera a irradiação difusa composta das subcomponentes circunsolar vinda diretamente do disco solar e a isotropicamente distribuída vinda do resto do hemisfério celeste. Para cada componente é dado um peso específico, através do índice anisotrópico definido como a razão entre a irradiação direta na superfície horizontal e a irradiação extraterrestre: (A_I = / ). Este índice define a percentagem da irradiação difusa a ser tratada como circunsolar, com a percentagem restante considerada isotrópica. Em condições de céu claro, o índice anisotrópico tende a um, e a irradiação difusa circunsolar é mais fortemente ponderada do que a irradiação difusa isotrópica. Em condições de céu nublado, o índice anisotrópico tende a zero e toda a irradiação difusa é tratada como isotrópica. Esse comportamento é consistente com a irradiação difusa observada. A irradiação difusa circunsolar ponderada através do índice anisotrópico é projetada na superfície inclinada da mesma maneira que a componente da irradiação direta, sendo a porcentagem restante tratada como isotrópica. Portanto, a irradiação difusa que atinge o plano inclinado, no modelo proposto por Hay (1979) é dada pela seguinte equação:

Utilizando o mesmo desenvolvimento de Hay (1979), Perez e outros (1987) consideram a irradiação difusa sendo composta pela soma das componentes isotrópica, circunsolar e o brilho horizontal. Na versão mais simples do modelo (PEREZ et al., 1987), válida para superfícies com face inclinada para o Equador, a irradiação difusa circunsolar é considerada originar-se em um ponto no centro do disco solar. Toda energia decorrente do brilho horizontal origina-se na região infinitesimal de (0º) de elevação. Assim, a expressão de Perez e outros (1987) é escrita da seguinte maneira:

Onde, F₁ representa a contribuição da irradiação difusa circunsolar e F₂ é a contribuição da irradiação difusa horizontal. Esses dois parâmetros são considerados independentes um do outro e são determinados empiricamente através de dados observacionais.

MODELOS MODIFICADOS DE ESTIMATIVA DA IRRADIAÇÃO TOTAL

As medidas de irradiação global na superfície horizontal, em valores horários ou diários, são feitas em várias localidades, mas o mesmo não ocorre para as medidas de irradiação total em superfícies inclinadas. No entanto, avaliações de irradiação total em superfícies inclinadas são de fundamental importância nas estimativa da radiação absorvida por superfícies topográficas com inclinação natural, ou superfícies com inclinação forçada (caso de aquecedores solares).

Com a finalidade de estimar a irradiação incidente sobre superfícies inclinadas com a face voltada para o Equador, é testado um modelo matemático simples de estimativa que utiliza apenas a irradiação solar global observada na superfície horizontal em conjunto com modelos conceituais de estimativa das irradiações difusa e direta.

Utilizando as equações (2, 3 e 4), o modelo isotrópico de Liu e Jordan (1963) para a estimativa da irradiação total em uma superfície inclinada é escrito da seguinte maneira:

Como sugerido por Klein (1977), dividindo a equação (7) pela irradiação global medida na horizontal (), a Equação (7) pode ser escrita da seguinte maneira:

onde K_d = , é a razão entre as irradiações difusa e global na superfície horizontal.

A Equação (7) também pode ser escrita em função da fração da irradiação direta contida na irradiação global da seguinte maneira:

onde, K_BH = , é a razão entre as irradiações direta projetada na superfície horizontal e a global na superfície horizontal.

O mesmo procedimento foi adotado para os modelos de Hay (1979) e Perez e outros (1987).

RESULTADOS E DISCUSSÃO

Modelos conceituais de estimativa das irradiações difusa e direta

Os modelos modificados de estimativa da irradiação total em uma superfície inclinada exigem o conhecimento da irradiação global medida na superfície horizontal e da irradiação difusa na horizontal estimada através de modelos conceituais que correlacionam a componente difusa, através da fração difusa (K_d), com a irradiação global através do índice de claridade K_t = / , razão entre as irradiações global e no topo da atmosfera, isto é, K_d = f(K_t), ou o conhecimento da irradiação direta na horizontal estimada através de modelos conceituais que correlacionam a componente direta, através da fração direta (K_BH), com a irradiação global através do índice de claridade K_t, isto é, K_BH=f(K_t ).

Os modelos conceituais de estimativa das irradiações difusa e direta baseados no índice de claridade, são chamados do tipo Liu e Jordan (1960), a partir do momento em que esses autores foram os primeiros a propor tal correlação. Esse tipo de correlação é encontrado em diversos trabalhos para diferentes localidades (KLEIN, 1976; ERBS; KLEIN; DUFFIE., 1982; GUEYMARD, 1993; RICIERI, 1998; OLIVEIRA et al., 2002). No desenvolvimento dos modelos conceituais de estimativa da irradiação difusa e direta, foram utilizados os dados coletados de setembro de 1997 a agosto de 2000, totalizando 1034 valores que passaram no teste de consistência.

A Figura 1 mostra os valores diários de K_d, em função de K_t, e a curva ajustada por regressão linear múltipla. Nota-se que quando K_d tende a um, K_t tende a zero, e quando K_d tende a zero, K_t tende a um. No intervalo intermediário, a correlação entre essas duas variáveis pode assumir tendência linear, exponencial, ou polinomial, utilizando polinômios de diferentes graus. Dessa maneira, alguns autores correlacionam K_d=f(K_t) , em função de intervalos específicos de K_t.

Oliveira e outros (2002), trabalhando com dados coletados na cidade de São Paulo, concluem que o melhor ajuste é obtido com um polinômio de quarto grau, impondo limites máximos e mínimos para K_t.

Ricieri (1998) fez extensiva revisão dos modelos conceituais que correlacionam essas duas variáveis e concluiu que quando é considerada a partição diária, das irradiações global, difusa e no topo da atmosfera, a regressão linear múltipla polinomial de quarto grau fornece o melhor ajuste para K_d em função de K_t, para Botucatu-SP.

A equação de regressão linear múltipla obtida com valores diários para Botucatu-SP, no ajuste da curva polinomial de K_d em função de K_t, de quarto grau é dada pela Equação (10), com o coeficiente de determinação, R²=0,91.

A Figura 2 mostra os valores diários de K_BH, em função de K_t e a curva ajustada por regressão linear. Nota-se que, para pequenos valores de K_tK_BH tende a zero, e quando K_BH tende a um, K_ttambém tende a um. Portanto, acima de um determinado valor de K_t existe linearidade entre as duas variáveis, sugerindo que a correlação K_BH = f(K_t), deva ser feita em função de intervalos específicos de K_t. Nesse caso, K_BH= 0 para valores de K_t < 0,25, condições de céu nublado. Para valores de K_t > 0,25 a correlação entre essas duas variáveis assume tendência linear.

A equação de regressão linear obtida com valores diários para Botucatu-SP, no ajuste da curva polinomial de K_BH em função de K_t, é dada pela Equação (11), com o coeficiente de determinação, R²=0,90.

Avaliação dos modelos de estimativa da irradiação total na superfície inclinada

Nos modelos modificados de estimativa da irradiação total diária através K_d = f(K_t) e por K_BH = f(K_t), o procedimento adotado foi substituir, nos modelos de Liu e Jordan (1963), Hay (1979) e Perez e outros (1987), os valores observados da irradiação global na superfície horizontal R_G em conjunto com os valores da irradiação difusa estimada através da Equação (10), e, similarmente, substituir os valores observados de R_G em conjunto com os valores da irradiação direta estimada através da Equação (11).

Na Tabela 1 estão apresentados os valores diários da Raiz Quadrada do Erro Quadrático Médio (REQM) e do Erro Médio (EM) calculados para a irradiação total estimada na superfície inclinada obtida com os vários modelos de estimativa modificados por K_BH = f(K_t), K_d = f(K_t) e os respectivos valores calculados para os modelos originais.

Na Tabela 1 verifica-se que o modelo original de Liu e Jordan (1963) tem o maior valor de dispersão REQM = 0,611 MJ/m² e o maior erro-médio em módulo EM = -0,38 MJ/m², indicando que esse modelo subestima os valores da irradiação total em uma superfície inclinada. Esse resultado é esperado porque o modelo de Liu e Jordan corrige a irradiação difusa apenas com o fator geométrico 0,5. (1+cos(b)) = 0,96, correção suficiente apenas para condições de céu nublado onde a irradiação difusa vinda do céu apresenta comportamento uniforme, ou seja isotrópico. Contudo, para condições de céu parcialmente nublado ou claro, o modelo de Liu e Jordan subestima a quantidade da irradiação difusa que incide sobre a superfície inclinada, principalmente em condições de céu claro em conseqüência de que a maior percentagem da irradiação difusa é proveniente do disco solar, ou seja a componente circunsolar. Esse resultado, demonstra claramente que o modelo isotrópico é deficiente na estimativa da irradiação total em uma superfície inclinada em condições de céu parcialmente nublado ou claro, principalmente por causa do comportamento anisotrópico da irradiação difusa, tornando-se necessário que a correção seja maior do que a simples correção geométrica. O modelo de Hay (1979) superestima pouco os valores da irradiação total EM=0,01 MJ/m² e também tem o menor valor de dispersão REQM=0,36 MJ/m², demostrando que a introdução do índice anisotrópico (A_I = / ) no modelo, para caracterizar o estado do céu, e subdividir a irradiação difusa nas componentes circunsolar e isotrópica ponderadas por esse índice, produz excelentes resultados na estimativa da irradiação total em uma superfície inclinada. Em condições de céu claro, o índice anisotrópico tende a 1 e a componente circunsolar da irradiação difusa ponderada por esse índice tem um peso maior do que a difusa isotrópica. Em condições de céu nublado, o índice anisotrópico tende a zero e praticamente toda a irradiação difusa é tratada como isotrópica. Na formulação matemática do modelo de Perez e outros (1987), são incluídas as subcomponentes da irradiação difusa, circunsolar, isotrópica e brilho horizontal. Nota-se que de maneira geral, esse modelo superestima os valores da irradiação total na superfície inclinada e apresenta valores de dispersão comparável com o modelo de Hay (1979).

Comparando os resultados, obtidos para os modelos modificados, com os originais, verifica-se que para o modelo de Liu e Jordan (1963) e para as duas diferentes modificações feitas, existe melhoria significativa com EM e REQM, um pouco menores do que os resultados obtidos com o modelo original, embora continue subestimando os valores. Como esse modelo, na forma original, não considera condições que caracterizam o estado do céu, com o uso da irradiação difusa ou da direta estimada em função do índice de claridade, esta, pode ser a causa da melhoria no modelo de Liu e Jordan (1963). De maneira geral, a análise dos erros demostra que eles aumentam ligeiramente para todos os outros modelos, com pouco aumento da superestimativa e dos valores da dispersão, com similaridade para ambos os métodos de modificação quando comparados com os erros obtidos com os modelos originais, sendo que esses erros são menores para a modificação feita com a função K_BH = f(K_t), provavelmente devido à correlação linear entre as duas variáveis, e pode estimar a irradiação total com razoável exatidão.

VALIDAÇÃO DOS MODELOS

Para verificar o comportamento dos modelos no decorrer do ano, os mesmos foram validados para um intervalo mais curto de tempo – mensalmente.

A Figura 3 mostra os valores de EM obtidos mensalmente com os modelos originais de estimativa da irradiação total diária. Nota-se claramente que o modelo de Liu e Jordan (1963) subestima os valores de irradiação em praticamente todos os meses do ano, exceto janeiro e dezembro. A subestimativa, em torno de –0,7 MJ/m², é grande para os meses de outono e inverno, decorrência de que no modelo de Liu e Jordan a irradiação difusa é corrigida apenas pelo fator geométrico – uma correção válida somente para irradiação difusa isotrópica. No entanto, durante o inverno, na região de Botucatu – SP, onde prevalecem condições de céu claro e a irradiação difusa tem caráter anisotrópico, e a correção geométrica é insuficiente, torna-se necessário incluir os efeitos da anisotropia na correção da irradiação difusa. Os modelos de Hay (1979) e Perez e outros (1987) superestimam os valores em todos os meses do ano, com valores máximos em torno de 0,4 MJ/m² durante o inverno, sendo os valores de EM ligeiramente menores para o modelo de Hay (1979).

Os valores de REQM, obtidos mensalmente com os modelos de estimativa da irradiação total diária, estão mostrados na Figura 4. Observa-se que para todos os modelos, os REQM são maiores durante outono e inverno do que durante a primavera e verão, e o modelo de Hay tem os menores valores de REQM em todos os meses, sendo este o melhor modelo para estimar a irradiação total diária na superfície inclinada em Botucatu – SP, por apresentar os menores erros-médios e menores valores de dispersão.

Na Figura 5 estão mostrados os EM resultantes da validação mensal da irradiação total diária na superfície inclinada utilizando os modelos de estimativa de Liu e Jordan (1963), Perez e Seals (1987) e Hay (1979), modificados por K_d = f(K_t).

Na análise dos erros verifica-se que, novamente, o modelo de Liu e Jordan subestima os valores da irradiação, principalmente durante o outono e o inverno, enquanto que os outros modelos superestimam os valores em todos os meses do ano. Quando os erros são comparados com os obtidos nos modelos de estimativa originais (Figura 4), verifica-se que os valores de EM permanecem quase inalterados para o modelo de Liu e Jordan (1963) durante os meses do inverno e outono e aumentam no começo do verão e final da primavera. O modelo de Hay (1979), apresenta pouco aumento dos EM durante a primavera e verão, e atinge um máximo de 0,5 MJ/m² durante o inverno, o mesmo comportamento é observado para o modelo de Perez e outros (1987), onde durante o inverno o modelo apresenta erros da ordem de 0,7 MJ/m².

Na Figura 6 estão mostrados os REQM resultantes da validação mensal da irradiação total diária na superfície inclinada utilizando os modelos de estimativa de Liu e Jordan (1963), Perez e outros (1987). e Hay (1979), modificados por K_d = f(K_t). Na análise, verifica-se que os REQM seguem o mesmo comportamento dos EM, permanecem quase que inalterados para o modelo de Liu e Jordan (1963), aumentam em todos os meses para o modelo de Perez e outros (1987) principalmente durante o outono e inverno atingindo um pico máximo no mês de julho em torno de 1,2 MJ/m², e têm pouco aumento durante a primavera e verão, o modelo de Hay (1979) também apresenta aumento nos valores da dispersão com pico em torno de 0,8 MJ/m² durante o outono e inverno e ficam praticamente inalterados durante a primavera e verão.

Quando os modelos são modificados através de K_d , excetuando-se o modelo de Liu e Jordan (1963), onde os erros permanecem praticamente inalterados, tanto os EM quanto os REQM aumentam para todos os outros modelos, sendo extremamente elevados durante os meses de inverno e outono.

Na Figura 7 estão mostrados os EM resultante da validação mensal da irradiação total diária na superfície inclinada utilizando os modelos de estimativa de Liu e Jordan (1963), Perez e outros (1987), Hay (1979), modificados por K_BH = f(K_t) Na análise dos erros verifica-se que o modelo de Liu e Jordan continua subestimando os valores da irradiação total, principalmente durante o outono e inverno, e superestima em alguns meses da primavera e verão, no entanto são encontradas melhorias significativas, tornando esse modelo compatível com os outros modelos modificados de estimativa da irradiação total diária na superfície inclinada. Os modelos de Perez e outros (1987) e Hay (1979) superestimam a irradiação total em todos os meses do ano, mas nesse caso, os EM são similares e praticamente iguais aos obtidos com os modelos originais, sendo novamente o modelo de Hay (1979) o que apresenta os menores erros.

Na Figura 8 estão mostrados os REQM resultantes da validação mensal da irradiação total diária na superfície inclinada utilizando os modelos de estimativa de Liu e Jordan (1963), Perez e outros (1987), Hay (1979), modificados por K_BH= f(K_t) Na análise dos REQM verifica-se que o modelo de Liu e Jordan apresenta valores de dispersão menores do que os obtidos com o modelo original tornando o modelo modificado compatível com o de Hay tendo praticamente os mesmos valores de dispersão e com menores valores do que o modelo de Perez e outros.

CONCLUSÕES

Os resultados obtidos com os modelos demonstram que o modelo de Liu e Jordan (1963) subestima a irradiação total principalmente durante o outono e inverno, enquanto que os modelos de Perez e outros (1987) e Hay (1979) superestimam em todos os meses do ano.

Para os modelos modificados, a comparação entre os erros com as duas modificações feitas, mostra que os erros têm um aumento considerável para todos os modelos, exceto para o de Liu e Jordan (1963), quando a modificação e feita através de K_d = f(K_t), demonstrando que quando a irradiação difusa é estimada através de regressão linear múltipla com um polinômio de 4º grau, e usada nos modelos modificados, não produz resultados satisfatórios quando se utiliza um conjunto de dados diferente daqueles nos quais ela foi obtida.

Para os modelos modificados através de K_BH = f(K_t), os erros aumentam pouco em comparação com os obtidos com os modelos originais, obtendo-se melhorias significativas para o modelo de Liu e Jordan tornando-o compatível com os outros modelos, esse resultado demonstra que quando a irradiação direta na horizontal é estimada através de regressão linear com um polinômio de 1º grau, e usada nos modelos modificados, produz resultados satisfatórios quando se utiliza um conjunto de dados diferente daqueles em que ela foi obtida, e podem ser usados para estimar a irradiação total em uma superfície inclinada com razoável precisão, sendo o modelo de Hay ligeiramente superior aos outros e o mais indicado para ser utilizado na estimativa da irradiação total em Botucatu-SP.

Agradecimentos

Os autores agradecem a FAPESP pelo apoio financeiro que possibilitou a montagem da Estação de Radiometria Solar da UNESP em Botucatu-SP.

Recebido em 21 out. 2003 / Aceito em 22 jun. 2004

Received oct. 21, 2003 / Accepted jun. 22, 2004

NOTAS SOBRE AUTORES

Dr. José Scolar é bacharel em Física pela Pontifícia Universidade Católica de São Paulo em 1976. Mestrado em meteorologia pelo Instituto Astronômico e Geofísico da Universidade de São Paulo em 1983. Doutorado pela Faculdade de Ciências Agronômicas da UNESP – Campus de Botucatu em 2003. Atualmente é pesquisador aposentado pelo Instituto de Pesquisas Meteorológica (IPMet-UNESP). Suas atividades concentram-se na investigação dos processos de transporte turbulento na camada limite planetária, e na investigação do comportamento da radiação solar recebidas em superfícies inclinadas.

Dr. Dinival Martins é engenheiro florestal pela Escola Superior da Agricultura " Luiz de Queiroz " – USP, Piracicaba-SP em 1977. Mestrado em Energia Nuclear aplicada à Agricultura pelo centro de Energia Nuclear na Agricultura – USP, Piracicaba em 1981. Doutorado em Solos e Nutrição de Plantas – USP, Piracicaba em 1986. É professor da Faculdade de Ciências Agronômicas da UNESP, Campus de Botucatu-SP onde leciona as disciplinas de graduação: Meteorologia e Climatologia; Climatologia Aplicada e auxilia na disciplina de Climatologia Agrícola. No curso de Pós-graduação em Energia na Agricultura, é responsável pelas disciplinas: Termodinâmica Atmosférica; e Aproveitamento de Energia Solar. Desenvolve pesquisas com Fontes Alternativas de Energia Natural.

Dr. João Francisco Escobedo é licenciado em Física pela UNESP em 1974 e doutorado em Física pelo Instituto de Física de São Carlos/USP em 1987. Atualmente é professor adjunto pelo Departamento de Recursos Naturais da UNESP-Botucatu e desenvolve atividades acadêmicas para os cursos de graduação em agronomia, lecionando a disciplina Climatologia Agrícola, e na Pós-graduação é responsável pelas disciplinas Radiação Solar, Técnicas e Medidas de Parâmetros Agrometeorológicos e Processos Térmicos da Radiação Solar para as áreas Energia na Agricultura e Irrigação e Drenagem. Desenvolve pesquisas na área de radiação solar com linhas direcionadas para Instrumentação, Medidas, Modelos e Conversão Energética.

• BEHR, H.D. Solar radiation on tilted south oriented surface validation of transfer-models. Solar Energy, [S.l.], v. 61, n. 6, p. 399-413, 1997.
• BUGLER, J.W. The determination of hourly insolation on an inclined plane using a diffuse irradiance model based on hourly measured global horizontal insolation. Solar Energy, [S.l.], v. 19, p. 477-91, 1977.
• FEUERMANN, D.; ZEMEL, A. Validation of models for global irradiance on inclined planes. Solar Energy, [S.l.], v. 48, n. 1, p. 59-66, 1992.
• ERBS, D. G.; KLEIN, S. A.; DUFFIE, J. A. Estimation of the diffuse radiation fraction for hourly, daily and monthly-average global radiation. Solar Energy, [S.l.], v. 28, p. 293-302, 1982.
• GUEYMARD, C. An anisotropic solar irradiance model for tilted surfaces and its comparison with selected engineering algorithms. Solar Energy v. 38 n. 5 p. 367-386. 1987.
• _____. Mathematically integrable parameterization of clear-sky beam and global irradiances and its use in daily irradiation applications. Solar Energy, [S.l.], v. 50, n. 2, p. 385-97, 1993.
• HAY, J. E. Calculation of monthly mean solar radiation for horizontal and inclined surfaces. Solar Energy, [S.l.], v. 23, p. 301-07, 1979.
• IQBAL, M. An introduction to solar radiation. Canada: Academic Press, 1983.
• KLEIN, S.A. Calculation of monthly average insolation on tilted surfaces. Solar Energy, [S.l.], v. 19, p. 325-29, 1977.
• KLUCHER, T. M. Evaluation of models to predict insolation on tilted surfaces. Solar Energy, [S.l.], v. 23, p. 111-14, 1979.
• LIU, B. Y. H.; JORDAN, R. C. The interrelationship and characteristic distribution of direct, diffuse and total solar radiation. Solar Energy, [S.l.], v. 4, p. 1-19, 1960.
• LIU, B. Y. H.; JORDAN, R. C. The long term average performance flat plate solar energy collectors. Solar Energy, [S.l.], v. 7, p. 53-74, 1963.
• NOTTON, G.; MUSELLI, M.; LOUCHE, A. Two estimation methods for monthly mean hourly total irradiation on tilted surfaces from monthly mean daily horizontal irradiation from solar radiation data of Ajaccio, Corsica. Solar Energy, [S.l.], v. 57, n. 2, p. 141-53, 1996.
• OLIVEIRA, A. P. et al. Correlation models of diffuse solar-radiation applied to the city of São Paulo, Brazil. Applied Energy, [S.l.], v. 71, p. 59-73, 2002.
• PEREZ, R.; SEALS, R.; INEICHEN, P.; STEWART, R.; MENICUCCI, D. A new simplified version of the Perez diffuse irradiance model for tilted surfaces. Solar Energy, [s.l.], v. 39, n. 3, p. 221-31, 1987.
• RICIERI, R. P. Modelos de estimativas e avaliação dos métodos de medida da radiação solar difusa 1998. 89f. Tese (Doutorado em Energia na Agricultura)-Faculdade de Ciências Agronômicas, Universidade Estadual Paulista, Botucatu, 1998.
• TEMPS, R. C.; COULSON, K. L. Solar radiation incident upon slopes of different orientations. Solar Energy, [S.l.], v. 19, p. 179-84, 1977.

Datas de Publicação

• Publicação nesta coleção
  24 Jan 2007
• Data do Fascículo
  Dez 2003

Histórico

• Aceito
  22 Jun 2004
• Recebido
  21 Out 2003