The SASIS model: validation of furrow irrigation simulation
Rev. bras. eng. agríc. ambient. vol.12 no.1 Campina Grande Jan./Feb. 2008
http://dx.doi.org/10.1590/S1415-43662008000100002
O modelo SASIS: validação da simulação da irrigação por sulco
Roberto V. PordeusI; Carlos A. V. de AzevedoII; José Dantas NetoII; Vera L. A. de LimaII; Márcia R. de Q. A. AzevedoIII; José de A. de MatosI
IDCA/UFERSA, BR 110, km 47, Costa e Silva, CEP 59625-900, Mossoró, RN. Fone: (84) 3315-1739. Email: rpordeus@ufersa.edu.br
IIUAEAg/UFCG, Av. Aprígio Veloso 882, Bodocongó, CEP 58109-970, Campina Grande, PB. Fone: (83) 3310-1056; (83) 3310-1373; (83) 3310-1349. Email: cazevedo@deag.ufcg.edu.br; zedantas@deag.ufcg.edu.br; antunes@deag.ufcg.edu.br.
IIIUEPB, CEP 58109-753, Campina Grande, PB. Fone: (83) 3333-3860. Email: marciarqaa@ibest.com.br
Nesta pesquisa se objetivou validar a simulação da irrigação por sulco, realizada pelo modelo SASIS através do procedimento matemático de ondas cinemáticas. Utilizam-se, na validação do modelo SASIS, 7 conjuntos de dados que representam diferentes condições de campo, em que o comprimento de sulco variou de 67 a 403 m, a declividade de 0,0016 a 0,0173 m m-1, a vazão de 1,1 a 2,0 L s-1 e o tipo de textura do solo de franco-arenoso a franco silte-argiloso. Para maior consistência da validação, as simulações realizadas por este modelo foram comparadas com as do modelo SIRMOD, além das comparações com as medições de campo. Para a maioria dos dados de campo estudados, as simulações da fase de avanço pelo modelo SASIS mostraram discrepâncias no tempo de avanço, no final da área, inferiores àquelas dadas pelo modelo SIRMOD, em relação às medições de campo. Novamente, para grande parte das condições de campo analisadas as discrepâncias entre os valores simulados pelo modelo SASIS e os medidos em campo do tempo de avanço no final da área, foram pequenas, a ponto de não comprometerem o prognóstico de parâmetros de dimensionamento, avaliação e manejo dos sistemas de irrigação por sulco.
Palavras-chave: irrigação superficial, modelo de ondas cinemáticas, manejo de água
The aim of this research was to validate the simulation of furrow irrigation performed by the SASIS model through the mathematical procedure of kinematic waves. In the validation of the SASIS model, 7 data sets representing different field conditions were used, in which the furrow length varied from 67 to 403 m, the field slope from 0.0016 to 0.0173 m m-1, the inflow from 1.1 to 2.0 L s-1 and the soil texture type from sandy loam to silty clay loam. In order to give more consistency to the validation, the simulations accomplished by this model were compared with those from the SIRMOD model, and also compared with data from the field measurements. For most field data studied, the simulations of the advance phase by SASIS model presented discrepancies in the advance time at the end of the area, being inferior to the ones given by SIRMOD model, in relation to field measurements. Again, for the greater part of field conditions analyzed, the discrepancies between the values simulated by SASIS model and those measured in field of the advance time at the end of the area were small, to the point of not compromising the prognostic of design, evaluation and management parameters of furrow irrigation systems.
Key words: surface irrigation, kinematic wave model, water management
A simulação matemática da irrigação por superfície tem processos complexos na hidráulica do escoamento superficial, os quais têm sido simulados por modelos computacionais com amplo grau de complexidade e precisão: (Strelkoff & Katopodes, 1977; Elliott et al., 1982; Walker & Humpherys, 1983; Strelkoff & Souza, 1984; Rayej & Wallender, 1985; Walker & Skogerboe, 1987; Azevedo, 1992). Tais modelos simulam as fases de avanço, armazenamento, depleção e recessão da água ao longo da superfície do solo e o volume de água infiltrado, escoado e percolado.
O modelo hidrodinâmico resolve completamente as equações de Saint-Venant (equação da continuidade e do movimento), isto é, sem assumir aproximações simplificativas; o modelo de zero inércia baseia-se também nessas equações negligenciando, porém, na equação do movimento, os termos inerciais; já o modelo de ondas cinemáticas se baseia apenas na equação da continuidade e na equação de Manning. Esses três modelos são baseados na solução numérica das equações diferenciais da conservação de massa e de energia e, segundo vários estudos, predizem mais adequadamente o processo da irrigação por superfície que o modelo do balanço de volume de água, Raghuwanshi & Wallender (1996).
O modelo hidrodinâmico é o mais preciso dos modelos, razão por que é considerado padrão, porém os modelos zero inércia e ondas cinemáticas, apesar de serem mais simples apresentam, em termos gerais, precisão satisfatória sendo, portanto, boas alternativas na simulação da irrigação por superfície (Strelkoff & Katopodes, 1977; Rayej & Wallender, 1985; Azevedo, 1992; Valiantzas, 1999).
Os diferentes modelos de simulação da irrigação por superfície foram desenvolvidos para simular um evento de irrigação isolado, assumindo-se que não existe variabilidade espacial nos parâmetros de campo (infiltração, rugosidade, declividade e seção transversal de fluxo); na prática, tem-se verificado a validade desta hipótese, tendo em vista que as simulações se têm aproximado bastante das medições de campo das fases, porém as variabilidades temporais nesses parâmetros são sempre levadas em conta, uma vez que, para a avaliação de qualquer evento de irrigação ao longo da estação de cultivo, novas medições dos parâmetros de campo são realizadas.
Muitos aspectos de uma irrigação por sulco devem ser considerados na qualidade de sua simulação, dentre os quais se citam o modelo matemático empregado e o grau de precisão dos dados de campo, como mais importantes. A qualidade da simulação reflete na precisão do prognóstico de parâmetros extremamente relevantes no dimensionamento, manejo e desempenho da irrigação por sulco, tais como: volumes infiltrados, escoados e percolados, e parâmetros de desempenho da irrigação, que são eficiência de aplicação d'água, perdas d'água por percolação e escoamento, eficiência de armazenamento e uniformidade de distribuição.
Objetivou-se, com esta pesquisa, validar a simulação da irrigação por sulco com fluxo contínuo realizada pelo modelo matemático computacional SASIS - Software Aplicado à Simulação da Irrigação por Superfície.
Equações do modelo de ondas cinemáticas
No modelo matemático de ondas cinemáticas, utilizado nesta pesquisa, se assume que não existe variação da altura de fluxo com a distância, isto é, ¶y/¶x = 0, negligenciando-se completamente a equação do movimento sem a qual nada pode ser dito com respeito à dinâmica da forma do perfil superficial de fluxo, ficando a equação da continuidade indeterminada no termo ¶A/¶t; para solucionar este problema, ressalta-se haver uma relação única que descreve vazão em função da área de fluxo, substitui-se, então, a equação do movimento pela equação de Manning. O escoamento, estudado apenas cinematicamente, assemelha-se à propagação de uma onda cinemática que entra em colisão, razão por que este modelo é designado ondas cinemáticas. A designação de modelo de escoamento uniforme, também usada, se deve, como é óbvio, à fundamentação da equação de Manning, que descreve o fluxo uniforme.
Assim sendo, as equações que constituem o modelo de ondas cinemáticas, são:
Equação da continuidade:
Equação de Manning:
A - área da seção transversal de fluxo, m2
t - tempo de ocorrência, min
x - distância de avanço da água, m
t - tempo de oportunidade de infiltração, min
Z - volume infiltrado acumulado por unidade de comprimento de sulco, m3 m-1
Q - vazão, m3 s-1
n - coeficiente de rugosidade de Manning, m-1/3 s
So - declividade do terreno, m m-1
R - raio hidráulico, m
A equação de Manning foi utilizada, nesta análise, para gerar relação única entre vazão e seção hidráulica. Elliott et al. (1982) propuseram relações empíricas para a seção hidráulica, dadas por:
y - altura de fluxo, m
A - área da seção transversal de escoamento, m2
s1, s2, r1 e r2 - parâmetros empíricos que dependem da forma do sulco
Os parâmetros geométricos s1, s2, r1, e r2, são considerados constantes para um sulco prismático, ou seja, em um evento de irrigação não variam ao longo do comprimento do sulco, válidos apenas para um evento de irrigação específico, haja vista que seus valores podem mudar entre uma irrigação a outra.
Esta hipótese garante que as funções potenciais descrevem adequadamente relações entre altura de fluxo, área, largura da superfície livre d'água, raio hidráulico etc; enfim da equação de Manning e se atendendo à Eq. 4, chega-se a:
em que: Q é a descarga de fluxo em m3 s-1; n é o coeficiente de Manning em m-1/3 s e So é a declividade média do sulco em m m-1, de onde se tira Q cuja derivada, em conjunto com a equação da infiltração, se substitui na Eq. 1, seguindo-se uma cadeia de transformações que, afinal, compõem o desenvolvimento do modelo (Walker & Humpherys, 1983).
Com base nas Eqs. 2 e 5, este tipo de modelo não é aplicável a sulcos quando a declividade é muito pequena, ou seja, tende a zero; na realidade, sua precisão diminuirá quando So se aproximar de zero. Strelkoff & Katopodes (1977) acharam que ele simula melhor as condições de escoamento quanto maior for a declividade longitudinal.
Usando-se a relação da Eq. 3, a equação de Manning se tornará:
sendo n o coeficiente de rugosidade de Manning em m-1/3 s, So a declividade do sulco em m m-1 e r1, e r2 constantes empíricas ajustadas às medições de campo atuais da geometria do sulco; a e m são constantes empíricas.
Solução numérica do modelo
A solução numérica das equações do modelo de ondas cinemáticas se baseia no conceito do controle de volume constituído de células individuais deformáveis, descrevendo o perfil superficial de fluxo em um evento de irrigação. A deformação mais acentuada do perfil superficial de fluxo ocorre nas regiões de contorno do sistema. A solução dessas equações leva em consideração que a maior deformação das células ocorre na seção de entrada do fluxo ou, ainda, na seção de saída. Quando se baseia na maior deformação na seção de entrada, as células desenvolvem uma velocidade avante definindo, assim, um sistema Lagrangeano, ou seja, com células em movimento, ocorrendo sua expansão tanto ao longo da distância como do tempo, enquanto se a solução das equações é baseada na deformação na seção de saída, as células se tornam estacionárias, estabelecendo um sistema Euleriano no qual a expansão das células ocorre apenas ao longo do tempo.
Para solução numérica espacial das equações do modelo de ondas cinemáticas utilizou-se, nesta pesquisa, o procedimento de integração Euleriana com aproximação de primeira ordem, apresentado por Walker & Humpherys (1983) e Wallender (1986), o qual resulta em duas equações algébricas mais estáveis e mais fáceis de serem resolvidas em microcomputadores. Conceitualmente, a aproximação considera o perfil superficial e subsuperficial de água ao longo da porção da área umedecida durante etapas seqüenciais de cálculo. A Figura 1 ilustra os perfis superficial e subsuperficial de fluxo nos tempos ti-1 e ti, identificando as células que os compõem. Durante cada etapa de cálculo o fluxo d'água avança a uma distância incremental, dx; por exemplo, durante o primeiro intervalo de tempo (primeira etapa de cálculo), se estende a uma distância dx1; no segundo intervalo de tempo, a uma distância dx2, e assim sucessivamente, podendo-se generalizar, para a distância da frente de avanço, xi , no tempo ti , da seguinte maneira:
em que dxk é o k-ésimo incremento de espaço, definido pelo avanço durante o intervalo de tempo, quando i = k, donde k é número de incrementos de tempo.
Uma célula típica é ilustrada na Figura 2, exibindo os perfis nas etapas de cálculo ti-1 e ti . As anotações J, M, L e R são introduzidas em cada célula para identificar as variáveis que descrevem as condições de fluxo com relação ao tempo e ao espaço. Desta maneira, as variáveis subscritas por J ou M se referem às condições de fluxo no tempo ti-1 e nas fronteiras à esquerda (montante) e à direita (jusante) da célula, respectivamente. Similarmente, L e R são subscritos nas fronteiras à esquerda e à direita da célula no tempo ti . Combinando-se as células de todos os incrementos de tempo ter-se-á uma grade de cálculo no plano (x,t), na qual as trajetórias de avanço e recessão podem ser traçadas (Figura 3). Observa-se, neste plano, que durante a fase de avanço as células são retangulares, exceto as da frente de avanço, que são triangulares por não existir fluxo na fronteira à jusante dessas células, nos tempos ti e ti-1 . Durante as fases de armazenamento e depleção todas as células são retangulares; já durante a recessão, a célula da extremidade à montante é triangular, pelas mesmas razões do comportamento do fluxo na frente de avanço mas as demais células são retangulares. A largura de cada célula é determinada pela distância da frente de avanço durante cada etapa de cálculo dt, que se torna, então, uma das incógnitas do problema, uma vez que se considerou dt constante; a isto se chama uma solução espacial das equações do modelo de ondas cinemáticas, porque dx é uma incógnita, enquanto o valor de dt é definido pelo usuário do modelo.
Para se obter a forma integrada da equação da continuidade (Eq. 1), primeiro se integra, analiticamente, com relação à distância e depois numericamente, quanto ao tempo. Tem-se, então:
em que a primeira integração resulta em:
e a segunda em:
Uma vez que durante a fase de avanço o fluxo vai diminuindo ao longo da distância, devido ao efeito da taxa de infiltração d'água no solo, caindo para zero na frente de avanço, a não linearidade no perfil, tanto superficial como subsuperficial, se torna bastante acentuada fazendo com que o fator ½ não seja apropriado para ponderar as condições de fluxo; desta forma, uma ponderação adequada deverá estar entre ½ a 1. Criam-se, então, os fatores q e f para ponderarem a não linearidade no perfil, respectivamente, com relação ao tempo e ao espaço. Esses fatores são extremamente importantes apenas durante a fase de avanço, porque depois desta fase ocorre tendência de se estabelecer condições estáveis de fluxo. A ponderação deve ser maior no tempo t + dt visto que durante a fase de avanço a não linearidade aumenta com o tempo, sendo máxima quando a água atinge o final da área (condições de fluxo nulas); além disso, também deve ser maior na seção de entrada uma vez que a altura de fluxo varia mais rapidamente nesta seção, e é nela que a vazão está sendo aplicada.
Substituindo-se, na Eq. 12, ½ por q, (1 - q), f, (1 - f) e se escrevendo as variáveis de fluxo em termos das anotações de uma célula computacional (Figura 2), tem-se:
sendo q e f coeficientes de ponderação temporal e espacial, respectivamente, cujos valores variam entre ½ e 1; em geral, tem-se tomado valores iguais a 0,65 e 0,51 para q e f, respectivamente. Walker (2001) em seu software SIRMOD III, usa 0,60 para ambos, isto é, o mesmo valor adotado nesta pesquisa.
A solução numérica é obtida resolvendo-se a Eq. 13 para cada célula na malha computacional (Figura 3) começando-se horizontalmente da esquerda para a direita, em cada etapa de cálculo. As únicas incógnitas em cada célula, são QR e AR; entretanto, uma vez que Q está sendo calculado explicitamente pela equação de Manning (Eq. 6), não se o considera uma incógnita; o mesmo caso é aplicado a Z (infiltração) o qual está sendo calculado pela equação de Kostiakov-Lewis; o que se faz é substituir as equações de Q e Z na Eq. 13 ficando esta, então, apenas com uma incógnita (AR).
Substituindo-se a Eq. 6 na Eq. 13 e se dividindo por qa/dt, chegar-se-á à seguinte equação não linear em AR:
Objetivando-se simplificar a Eq. 14, isolam-se as constantes e variáveis com valores conhecidos a partir da etapa de cálculo anterior nos coeficientes C1 e C2. Tem-se, então:
chegando-se à equação
A Eq. 17 é utilizada para células interiores e para a primeira célula, depois da primeira etapa de cálculo. Uma vez que a Eq. 17 é solucionada implicitamente (pelo método de Newton-Raphson) para cada célula separadamente, não se tem, então, uma matriz. Esta equação é usada de forma implícita para determinar AR e, depois, se determina explicitamente QR, pela equação de Manning (Eq. 6)
As Eqs. 15 a 17 geram soluções para todas as células, exceto para aquelas na frente de avanço (célula da extremidade do avanço), em que as condições de contorno resultam em AR, AJ, AM, ZR, ZJ e ZM iguais a zero. Para essas células, a formulação matemática das Eqs. 15 a 17, se reduz a:
Como AL é conhecida da solução da célula anterior e ZL é função única do tempo de oportunidade de infiltração, a única incógnita nas células da frente de avanço é, então, a distância incremental (dx); desta forma, a Eq. 18 se reduz a:
Quando o suprimento de água na entrada da área é cortado, assume-se que a área no contorno esquerdo, AL, cai imediatamente a zero. A partir do corte da água se inicia a fase de recessão e os cálculos são feitos de forma similar aos da fase de avanço e terminam quando a área molhada é reduzida a um valor igual ou menor que 5% do valor da área na entrada do sulco.
Os sistemas de irrigação por superfície apresentam duas condições de contorno de jusante: 1) uma célula triangular durante a fase de avanço e 2) uma célula retangular para a fronteira do campo. No caso (2), várias configurações são possíveis mas se tem usualmente, uma condição de QN = 0 (sulco fechado) ou um fluxo de saída uniforme; entretanto, em algumas condições a descarga pode passar sobre ou através de uma estrutura hidráulica com característica diferente caso em que é necessário se estabelecer uma avaliação da descarga da estrutura para descrever o escoamento na condição de contorno.
O contorno à montante é definido pelas condições de fluxo na cabeceira da área irrigada, na qual A = 0 e Q = 0 para t = 0 e t > tap (tempo de aplicação da água) e A = Ao para 0 < t < tap. As condições de fluxo no contorno à jusante são descritas pelas Eqs. 6 e 19, respectivamente, para as fases de avanço e armazenamento, com ocorrência de escoamento superficial. Quando o fluxo atinge o final do campo, as variáveis M e R não serão zero para todos os intervalos de tempo, para cujas condições, como pode ser generalizado para a drenagem livre, dique, ou condições controladas, a condição de contorno pode ser descrita pela equação de Manning (Eq. 6).
Dentre as fases da irrigação superficial simuladas o avanço é a mais importante, em virtude de contribuir, de forma mais decisiva, para as variações no tempo de oportunidade de infiltração principalmente na irrigação por sulco, o que resulta na desuniformidade da lâmina infiltrada; assim sendo, o modelo SASIS faz a simulação apenas da fase de avanço, considerando o tempo de depleção constante ao longo da área e igual ao tempo de corte do fluxo d'água no sulco, hipótese que leva em conta o fato de que na irrigação por sulco a duração da fase de recessão é muito curta, em razão não apenas da declividade da superfície do solo mas, também, do volume de água armazenado no sulco, no momento de corte da vazão, que é muito pequeno comparado ao de uma bacia ou de uma faixa, não contribuindo significativamente para o volume de água infiltrado.
Os dados de campo utilizados na validação do modelo SASIS corresponderam a quatro conjuntos de dados (PISG1, PISG2, PISG3 e PISG4) coletados nesta pesquisa, referentes às avaliações de campo de eventos de irrigação por sulco no Perímetro Irrigado de São Gonçalo, no município de Sousa, PB, publicados por Azevedo et al. (2001) e Pordeus et al. (2003); dois conjuntos de dados (AMALGACQ, propriedade privada, e GUFCQ, fazenda da Utah State University em Logan, USA) publicados por Azevedo (1992), empregados na demonstração do modelo SIRTOM, e um conjunto de dados (KWF-Kimberly Wheel Furrow) publicados por Walker & Skogerboe (1987). Esses conjuntos de dados (Tabela 1) representam diferentes condições de campo, em que o comprimento de sulco variou de 67 a 403 m, a declividade de 0,0016 a 0,0173 m m-1, a vazão de 1,1 a 2,0 L s-1 e o tipo de textura do solo de franco arenoso a franco silte-argiloso, e dados de avanço medido em campo (Tabela 2). Para maior consistência da validação do modelo SASIS, as simulações realizadas por este modelo foram comparadas com as do modelo SIRMOD (Walker, 1989), além das comparações com as medições de campo.
As curvas de avanço correspondentes às medições de campo e simulações pelos modelos SIRMOD e SASIS para a irrigação por sulco com fluxo contínuo, são apresentadas na Figura 4. Observa-se semelhança quando se comparam as curvas simuladas por ambos os modelos com as medições de campo. Na maioria das simulações, o tempo de avanço é superestimado pelos dois modelos. Para os dados de campo PISG1 observa-se, na Tabela 3 que, no final da área, o tempo de avanço é superestimado em 8 min, pelo modelo SASIS, e em 35,2 min pelo SIRMOD correspondendo, portanto, a uma variação de 20 e 110%, respectivamente; verifica-se, na Figura 4A, que o avanço simulado pelo SASIS ao longo da área se aproxima mais dos valores medidos em campo que o SIRMOD; já a curva de avanço simulada pelo modelo SIRMOD se acerca mais dos valores medidos em campo, apenas nos primeiros 50% do comprimento do sulco; no restante, a curva de avanço se torna menos linear em relação à simulada pelo SASIS e aos valores medidos em campo, à medida que se avizinha do final da área; os valores simulados pelo SASIS se aproximam dos valores medidos em campo do início da área até o seu final, o que não ocorre com o SIRMOD, o qual simulou 27,2 min a mais no final da área que o SASIS; para este caso, verifica-se melhor simulação do tempo de avanço pelo modelo SASIS.
Observa-se, no exemplo PISG2 e na Figura 4B, que mais uma vez os valores simulados pelo modelo SASIS se aproximam mais dos valores medidos em campo que o SIRMOD; o tempo de avanço no final da área é subestimado pelo SASIS em 5,8% e superestimado pelo SIRMOD em 64,6%; os tempos de avanço medidos em campo e simulados por esses modelos foram, respectivamente, de 69,65 e 113,6 min (Tabela 3); a curva de avanço simulada pelo SIRMOD apresenta o mesmo comportamento do exemplo anterior, isto é, nos primeiros 50% do comprimento do sulco os valores simulados se aproximam mais dos valores medidos com valores mais próximos no início da área, enquanto no restante da área a diferença é crescente sempre que se aproxima do final do campo.
Para os dados de campo PISG3 (Figura 4C), quando se comparam os valores de avanço medidos com os simulados vê-se que o SIRMOD se aproxima mais dos dados de campo que o SASIS; o SIRMOD subestimou o tempo de avanço no final da área, em 2,5%, enquanto o SASIS o superestima em 25%; para este exemplo, os tempos de avanço medido e simulado pelo SIRMOD e SASIS no final da área, foram de 24,23,4 e 30 min (Tabela 3) respectivamente e, embora a diferença percentual entre os tempos de avanço medido e simulado pelo SASIS tenha sido de 25%, a diferença absoluta foi de apenas 6 min e a diferença entre o tempo simulado pelos dois modelos foi de apenas 6,6 min, o que mostra consistência na simulação realizada pelo modelo SASIS.
Observa-se, na Figura 4D, para os dados de campo PISG4, que o avanço simulado pelo SIRMOD se acerca mais dos valores de avanço medidos que o SASIS, que superestimou o avanço praticamente em todo o comprimento do sulco, com valores maiores em torno do centro da área, o que ocorreu também com o SIRMOD, porém com menores diferenças; nesta mesma figura ainda se constata que, exatamente nesta região, ocorre mudança no comportamento da curva de avanço medido em campo mas, no final da área, o tempo de avanço simulado pelo SASIS se assemelha mais ao de campo; para este exemplo, os tempos de avanço medido no final da área e simulado pelos modelos SASIS e SIRMOD foram, respectivamente, de 66,65 e 56,6 min (Tabela 3); conclui-se, daí, que o SASIS subestimou o avanço em apenas 1 min, enquanto o SIRMOD em 9,4 min, correspondendo a variações percentuais de 1,5 e 14,2%, respectivamente, em relação aos dados de campo. Embora o SIRMOD tenha apresentado variação maior no final da área, o avanço simulado ao longo da área se assemelhou mais aos valores medidos, resultando em um prognóstico do perfil de infiltração mais coerente com o que aconteceu em campo.
No exemplo KWF o tempo de avanço medido no final da área foi de 200 min e os simulados pelos modelos SASIS e SIRMOD foram, respectivamente, 130 e 190 min constatando-se, então, que o modelo SASIS subestimou o avanço em 35% e o SIRMOD em 5% (Tabela 3); verifica-se, na Figura 4E, que nos primeiros 50% do comprimento da área os avanços simulados pelos dois modelos se assemelham aos medidos em campo, subestimando-os no restante da área, com maiores diferenças sempre que se aproxima do final da área.
Quanto aos dados de campo AMALGACQ nota-se, na Figura 4F, que o avanço simulado pelo modelo SASIS é similar, ao longo de toda a área, ao medido em campo, o que não ocorre com o avanço simulado pelo modelo SIRMOD; o tempo de avanço medido no final da área e o simulado pelos modelos SASIS e SIRMOD foi, respectivamente, 232, 265 e 67,3 min (Tabela 3), verificando-se que o modelo SASIS superestimou o avanço em 33 min, enquanto o SIRMOD o subestimou em 164,7 min, correspondendo a variações percentuais de 14,2 e 71%, respectivamente; desta forma, o modelo SASIS apresentou maior precisão na simulação.
Observa-se, para os dados GUFCQ, que na Figura 4G, as curvas de avanço obtidas pelos modelos SASIS e SIRMOD se assemelham, no geral, à curva de avanço medida em campo; entretanto, o tempo de avanço medido no final da área e simulado pelos modelos SASIS e SIRMOD, foi 28, 45 e 57,3 min (Tabela 3), respectivamente, apresentando superestimativas de 17 (SASIS) e 29,3 min (SIRMOD), que se referem a variações percentuais de 60,7 e 104,6%. O modelo SIRMOD o superestimou em 43,9% a mais que o SASIS.
Na maioria dos exemplos, os tempos de avanço simulados se assemelham mais aos valores medidos em campo nos primeiros 50% do comprimento do sulco porque, provavelmente, a partir desse trecho o avanço se vai tornando mais lento, devido ao efeito da infiltração d'água na velocidade do fluxo que diminuída ao longo do sulco, aumenta o grau de não linearidade na curva de avanço. É bem mais fácil simular uma reta que uma curva. No restante da área a diferença entre os valores simulados e os medidos cresce sempre que se aproxima do final da área e, no geral, os valores simulados pelos dois modelos não apresentam grandes diferenças com relação aos medidos em campo.
De acordo com Levien (2003), o avanço é muito mais rápido no início que no final do sulco, em virtude da infiltração contribuir para redução na velocidade de avanço da água ao longo do sulco; segundo o autor, a uniformidade de irrigação depende da uniformidade do tempo de oportunidade de infiltração da água nos diferentes pontos da parcela irrigada; por sua vez, esta uniformidade é maior quanto mais rápido for o tempo de avanço.
Constatou-se, ante o exposto, que em cinco exemplos (PISG1, PISG2, PISG4, AMALGACQ e GUFCQ), o tempo de avanço no final da área simulado pelo SASIS se aproximou mais dos valores medidos em campo contra dois exemplos (PISG3 e KWF) simulados pelo SIRMOD. Entre os exemplos estudados a maior diferença percentual entre o tempo de avanço medido e o simulado no final da área (Tabela 3) ocorreu para o exemplo GUFCQ que, pelo SASIS, foi de 60,7% e, pelo SIRMOD, de 104,6%.
A diferença entre os dados observados e os simulados pelos modelos (SASIS e SIRMOD) depende decisivamente do grau de precisão e representatividade dos parâmetros de entrada requeridos (seção transversal, declividade, rugosidade e parâmetros da equação de infiltração d'água no solo), conseguido quando da realização dos ensaios de campo.
Os resultados apresentados na Tabela 3 e na Figura 4 constatam a viabilidade da aplicação do modelo desenvolvido (SASIS) nesta pesquisa para simulação e racionalização da irrigação por sulco com fluxo contínuo, uma vez que as simulações se assemelham bastante aos valores medidos em campo e, ainda, o referido modelo apresenta grande facilidade de operação, possuindo integração total ao Windows em todas as suas versões, permitindo ao usuário não apenas a entrada de dados de maneira clara e prática mas, também, mudanças em alguns parâmetros do modelo, ajustando-se à realidade de campo em questão. Existem alguns módulos no software SASIS que utilizam as simulações na determinação de parâmetros de dimensionamento, avaliação e manejo da irrigação por sulco, a exemplo do módulo que investiga a vazão ótima, isto é, aquela que minimiza as perdas de água por percolação e escoamento.
1. Para a maioria dos dados de campo estudados as simulações da fase de avanço pelo modelo SASIS apresentaram discrepâncias no tempo de avanço no final da área, inferiores às identificadas pelo modelo SIRMOD, em relação às medições de campo.
2. Novamente, para grande parte das condições de campo analisadas as discrepâncias entre os valores simulados pelo modelo SASIS e medidos em campo no tempo de avanço no final da área foram pequenas, a ponto de não comprometerem o prognóstico de parâmetros de dimensionamento, avaliação e manejo dos sistemas de irrigação por sulco.
Os autores agradecem ao CNPq, pelo apoio financeiro ao trabalho.
Azevedo, C. A. V. de. Real-time solution of the inverse furrow advance problem. Agricultural and irrigation engineering department, Logon: Utah State University, 263p. 1992. Ph.D. Dissertation [ Links ]
Azevedo, C. A. V. de; Pordeus, R. V.; Lima, V. L. A. de; Dantas Neto, J. Behavior of the optimal discharge in furrow discharge in furrow irrigation with continuous flow. ASAE Annual International Meeting, Sacramento, Calinornia, paper n. 012165, 2001, p.1-13. [ Links ]
Elliott, R. L.; Walker, W. R.; Skogerboe, G. V. Zero inertia modeling of furrow irrigation advance. Journal of the Irrigation and Drainage Division, ASCE, v.108(IR3), p-179-195, 1982. [ Links ]
Levien, S. L. A. Irrigação por superfície: Sistemas de irrigação por superfície. Mossoró: ESAM, 2003. 20p. [ Links ]
Pordeus, R. V.; Azevedo, C. A. V. de; Lima, V. L. A. de; Dantas Neto, J. Field evaluation of water infiltration profile in opened and blocked furrow irrigation with continuous flow. ASAE Annual International Meeting, Las Vegas, paper n. 032244, 2003, p.1-12. [ Links ]
Raghuwanshi, N.; Wallender, W. W. Modeling of seasonal furrow irrigation. Journal of the Irrigation and Drainage Engineering, ASCE v.122, n.4, p.235-242. 1996. [ Links ]
Rayej, M.; Wallender, W. W. Furrow Irrigation Simulation Time Reduction. Journal of the Irrigation and Drainage Division, ASCE, v.111, n.2, p.134-146, 1985. [ Links ]
Strelkoff, T.; Katopodes, N. D. End Depth Under Zero-Inertia Conditions. Journal of the Hydraulics Division. ASCE, 103(HY1), p.699-611, 1977. [ Links ]
Strelkoff, T.; Souza, F. Modeling effect of depth on furrow infiltration. Journal of Irrigation and Drainage Engineering. ASCE. v.110, n.4, p.375-387, 1984. [ Links ]
Valiantzas, J. D. Explicit time of advance formula for furrow desing. Journal of Irrigationand Drainage Engineering. Davis, v.125, n.1, p.19-25, 1999. [ Links ]
Walker, W. R. Sirmod, Surface irrigation simulation software. Logan: Utah State University. 1989. [ Links ]
Walker, W. R. Software for simulation, design, and evaluation of surface irrigation - User's guide and technical documentation. Logan: Utah State University, Biological and Irrigation Engineering Department, 2001. 63p. [ Links ]
Walker, W. R.; Humpherys, A. S. Kinematic-wave furrow irrigation model. Journal of the Irrigation and Drainage Division, ASCE, New York, 109(IR4), p.377-392, 1983. [ Links ]
Walker, W. R.; Skogerboe, G. V. Surface irrigation: Theory and practice. 1.ed. New Jersey: Prentice-Hall, Inc., Englewood Cliffs, 1987, 470p. [ Links ]
Wallender, W. W. Furrow model with spatially varying infiltration. Transactions of the ASAE, St. Joseph, v.29, n.4, p.1012-1016, 1986. [ Links ]