Dimensionamento de terraços de infiltração pelo método do balanço volumétrico

Infiltration terrace design by volumetric balance method

ENGENHARIA DE IRRIGAÇÃO E DRENAGEM

Dimensionamento de terraços de infiltração pelo método do balanço volumétrico

Infiltration terrace design by volumetric balance method

Jarbas H. de Miranda^I; Sergio N. Duarte^II; Késia O. da Silva^III; Quirijn de J. van Lier^IV; Nilson A. Villa Nova^V

^IDepartamento de Ciências Exatas - ESALQ/USP. Av. Pádua Dias 11, CP 09, CEP 13418-900, Piracicaba,SP, Fone:(19) 3429-4283 Ramal: 210, E-mail: [e-mail]

^IIDepartamento de Engenharia Rural, ESALQ/USP, E-mail: [e-mail]

^IIIInstituição Moura Lacerda, E-mail: [e-mail]

^IVDepartamento de Ciências Exatas, ESALQ/USP, E-mail: [e-mail]

^VDepartamento de Ciências Exatas, ESALQ/USP, E-mail: [e-mail]

RESUMO

As práticas de conservação de solo permitem o controle de perdas de solo e água em áreas agriculturáveis, objetivando assim a maximização do lucro sem provocar a redução da capacidade produtiva. Ante essa importância, há necessidade de se adotar procedimentos técnicos para o dimensionamento de estruturas conservacionistas. Assim, objetivou-se, com o presente trabalho, o desenvolvimento de um modelo computacional em linguagem Visual Basic 6.0, para dimensionamento de terraços em nível. O modelo foi aplicado para 3 localidades do Estado de São Paulo (Araçatuba, Piracicaba e Natividade da Serra). Pelos resultados obtidos, os valores da largura da seção parabólica de canais de terraço em nível de 0,5 m de profundidade para as 3 localidades, com período de retorno de 5 anos, foram de 3,30, 3,34 e 3,80 m, respectivamente.

Palavras-chave: conservação de solo e agua, modelagem

ABSTRACT

Soil conservation practices allow soil and water loss control in agricultural areas, aiming the maximization of profit without provoking reduction of yield potential. There is need to use technical procedures to design soil conservation structures. Therefore, the aim of the present work was the development of a computational model, in Visual Basic 6.0 language, for design of contour terrace. The model was applied for three municipalities of São Paulo State (Araçatuba, Piracicaba and Natividade da Serra). From the results obtained, it may be concluded that the width of a parabolic section contour terrace channel for these localities, with a return period of 5 years, was 3.30, 3.34 and 3.80 m, respectively, with 0,5 m of depth.

Key words: soil and water conservation, modeling

INTRODUÇÃO

As práticas conservacionistas permitem o controle de perdas de solo e água em áreas agriculturáveis, objetivando a maximização do lucro sem provocar redução da capacidade produtiva. A erosão consiste no processo de desprendimento e arraste das partículas do solo causado pela ação da água e do vento, constituindo-se na principal causa da degradação de terras agrícolas (Pruski & Grieseler, 1996). As práticas de conservação do solo podem ser divididas em vegetativas, edáficas e mecânicas. As de caráter vegetativo são aquelas em que se utiliza a vegetação para proteger o solo contra a ação direta das chuvas e, conseqüentemente, minimizar o processo erosivo, enquanto as práticas edáficas são as que, com modificações no sistema de cultivo, além do controle da erosão, mantêm ou melhoram a fertilidade do solo. De acordo com o que se observa na prática, nem sempre as práticas edáficas e vegetativas são suficientes para o controle da erosão, principalmente em regiões em que ocorrem chuvas de grande intensidade. Neste caso, a adoção de procedimentos complementares torna-se necessária para reduzir a velocidade do escoamento superficial e, conseqüentemente, a capacidade de transporte de sedimentos, o que pode ser alcançado mediante a implantação de barreiras mecânicas como, terraços, canais escoadouros, bacias de captação de águas pluviais e barragens, entre outros.

O terraceamento em terras agrícolas é uma das práticas mais difundidas entre os agricultores brasileiros para o controle da erosão hídrica, consistindo basicamente na construção de terraços, os quais são estruturas compostas por um dique e um canal, dispostos no sentido transversal à declividade do terreno, formando obstáculos físicos com a intenção de reduzir a velocidade do escoamento superficial e ordenar o movimento da água sobre a superfície do terreno. Segundo Lombardi Neto et al. (1994), a taxa de escoamento superficial é o parâmetro mais importante para dimensionamento de terraços, porém ainda são freqüentes os erros na aplicação desta prática, principalmente no seu dimensionamento correto (Sparovek & Silva, 1997).

O custo de construção e manutenção de um sistema de terraceamento é relativamente alto; portanto, antes da adoção dessa tecnologia deve-se fazer um estudo criterioso sobre as condições locais, como clima, solo, sistema de cultivo, culturas a serem implantadas, relevo do terreno e equipamento disponível, para que se tenha segurança e eficiência no controle da erosão. O rompimento de um terraço pode levar à destruição dos demais que estiverem a jusante, com grandes prejuízos para a área cultivada (Cruciani, 1987).

O critério corrente para o dimensionamento de terraços em gradiente é a adoção de uma seção capaz de escoar a vazão de escoamento superficial causada por uma chuva de duração igual ao tempo de concentração da área de contribuição; já para os terraços de infiltração ou em nível, dimensiona-se a seção capaz de armazenar o volume de escoamento superficial causado por uma chuva mais longa. Cruciani (1989) sugere que se adote uma chuva de 24 h com um período de retorno de 5 ou 10 anos. Entretanto, este critério gera seções excessivamente grandes, sobretudo quando o escoamento superficial é calculado pelo método do Número da Curva.

Com o desenvolvimento da ferramenta computacional vêm surgindo modelos capazes de considerar um número maior de fatores intervenientes permitindo realizar o dimensionamento de forma mais precisa. Dentre os modelos utilizados para este fim pode-se citar o modelo matemático "Terraço for Windows", que é um programa computacional desenvolvido por Pruski & Griebeler (1996), aplicado para o dimensionamento de terraços em gradiente ou em nível utilizando-se, para tal, dados específicos da região e da cultura em questão, como precipitação pluvial máxima esperada para períodos de retorno escolhidos pelo projetista, tipo de solo, taxa de infiltração básica de água no solo, declividade do terreno, manejo do solo, de culturas e de resíduos culturais e altura da crista do terraço, que pode ser construído em função das condições topográficas do terreno e do equipamento disponível para a sua construção.

Entretanto, este programa depende da adoção de determinado período de retorno para a chuva de projeto, obtida por meio da relação intensidade-duração-freqüência da localidade em questão. Como se sabe, muitos estados brasileiros ainda não dispõem dessas relações espacializadas; além disso, o escoamento superficial nem sempre possui o mesmo período de retorno da chuva que o gerou. Assim, considerando-se esses aspectos objetivou-se, com o presente trabalho, o desenvolvimento de um modelo computacional, que pode ser aplicado para avaliação do desempenho de terraços em nível, levando-se em conta os processos de escoamento superficial, infiltração de água no solo e armazenamento de água no terraço.

DESENVOLVIMENTO DO MODELO

O modelo denominado 'Terradim' foi desenvolvido no Departamento de Ciências Exatas da Escola Superior de Agricultura "Luiz de Queiroz" - ESALQ/USP, em linguagem de programação Visual Basic 6.0 (Fig. 1), O programa opera realizando um balanço volumétrico diário do volume de água armazenado no canal do terraço, durante um período equivalente a uma série de vários anos de dados de precipitação e consiste basicamente de duas etapas: 1. determinação do escoamento superficial diário para uma série longa de dados pluviométricos, mediante a aplicação da metodologia do Número da Curva (Curve Number) e 2. determinação do número de rompimentos da capacidade total de armazenamento de água no canal do terraço.

1. Determinação do escoamento superficial

O primeiro procedimento a ser tomado é determinar o valor diário do escoamento superficial, mediante a aplicação do método do Número da Curva (CN).

O modelo permite ao usuário marcar os meses em que será realizado o preparo do solo e, em função disso, interpola a tabela do SCS (USA, 1972) para determinar os valores de CN classes I e III, mediante a Tabela 1.

A retenção potencial é função do Número da Curva (CN), cujo valor é tabelado e depende da cobertura do solo, da condição hidrológica, tipo de solo e umidade antecedente do solo. Quanto maior o valor de CN maior também é o escoamento superficial e menor a retenção potencial. Por outro lado, quando se diminui o valor deste parâmetro o escoamento decresce e a retenção aumenta (Eq. 1).

(1)

em que:

R - retenção potencial máxima, mm

CN_I,II,III - número da curva, adimensional

Após o cálculo da retenção potencial máxima, o escoamento superficial é obtido pela Eq. 2.

(2)

em que:

E_(i) - escoamento superficial diário, mm

ppt_(i) - precipitação total diária, mm

2. Estimativa do período de retorno para a seção adotada

Uma vez obtidos os valores diários de escoamento superficial, o modelo determina o número de rompimentos do terraço durante a série histórica de dados de chuva, ou seja, quantifica o número de vezes que o volume máximo de capacidade de contenção de água pelo terraço foi ultrapassado. Com base neste valor é possível se obter o período de retorno pela relação entre o tamanho da série de dados e o número de rompimentos (Eq. 3).

(3)

em que:

TR - período de retorno, anos

TS - tamanho da série de dados, anos

NR - número de rompimentos

Para a determinação do número de rompimentos do terraço é conveniente a obtenção de alguns parâmetros, os quais serão abordados a seguir:

2.1. Parâmetros de entrada: primeiro o modelo é alimentado pelos seguintes dados de entrada - área de contribuição do terraço (Ac, ha), largura superior e altura do terraço com seção parabólica (B e H, m), valor da velocidade de infiltração básica (VIB, mm h^-1), comprimento do terraço (L, m) e Número da Curva (CN).

2.2. Determinação da seção do terraço, volume máximo de acumulação e volume diariamente escoado: após o modelo ter sido alimentado pelos parâmetros de entrada, o procedimento na seqüência é a obtenção da seção do terraço, obtida pela Eq. 4:

(4)

em que:

S - seção do terraço, m²

B - base superior do terraço, m

H - altura da lâmina de água no terraço, m

Uma vez obtida a seção do terraço, o volume máximo de acumulação de água pelo canal é obtido pelo produto entre a seção e o comprimento do terraço (Eq. 5).

V_máximo = SL

(5)

em que:

V_máximo - volume máximo de capacidade do terraço, m³

L - comprimento do terraço, m

De posse do valor do escoamento é possível se determinar o volume diariamente escoado, o qual é obtido pelo produto entre o escoamento superficial diário e a respectiva área de contribuição (Eq. 6).

VDE_(i) = E_(i)Ac

(6)

em que:

VDE_(i) - volume diariamente escoado, no dia (i), m³

E_(i) - escoamento superficial, no dia (i), m

Ac - área de contribuição do terraço, m²

2.3. Determinação dos parâmetros geométricos da área molhada do canal do terraço: A seção molhada, base, altura e perímetro molhado foram determinados paras o dias anterior (i-1) e posterior (i+1) pelas seguintes equações:

(7)

em que:

s_(i-1) - seção molhada do dia anterior, m²

VAT_(i-1) - volume de água no terraço no dia anterior, m³

(8)

em que:

s_(i+1) - seção molhada do dia posterior, m²

VDE_(i) - volume diariamente escoado, m³

(9)

em que:

b - base superior molhada do terraço, para os dias anterior e posterior, respectivamente, m

A altura de lâmina de água no terraço é obtida pela Eq. 10:

(10)

em que:

h -altura da lâmina de água no terraço, para os dias anterior e posterior, m

A determinação do perímetro molhado é obtida pela Eq. 11:

(11)

em que:

P_{(i - 1)} - perímetro molhado do dia anterior, m

De posse desses valores determina-se o valor do perímetro molhado médio pela Eq. 12:

(12)

em que:

P_(i) - perímetro molhado médio, m

Com base no valor do perímetro molhado médio, determina-se o valor do volume infiltrado (VI) (Eq. 13):

VI_(i) = VIBP_(i)L

(13)

Conhecendo-se o volume infiltrado pode-se determinar o volume de água no terraço para o dia corrente por meio de um balanço de volumes (Eq. 14).

VAT_(i) = VAT_(i-1) + VDE_(i) - VI_(i)

(14)

em que:

VAT_(i) - volume de água no terraço, m³

Uma vez obtido o valor do volume de água contido pelo terraço, compara-se com o volume máximo de capacidade do terraço e, assim, consegue-se obter quantas vezes ocorreu o rompimento do terraço. De posse do número de vezes em que o terraço foi rompido e no número de anos da série, é possível se estimar o período de retorno para a seção adotada pelo projetista.

As etapas do modelo podem ser conferidas no fluxograma apresentado na Figura 2.

APLICAÇÃO DO MODELO

O presente modelo, "Terradim", foi aplicado para 3 regiões com comportamentos distintos de chuva no Estado de São Paulo. As localidades selecionadas seguiram o critério de ter baixa, média e alta intensidade de chuvas e foram, respectivamente: Araçatuba (Lat.: 21º12'S, Long.: 50º25'W Altitude: 390 m), Piracicaba (Lat.: 22º42'S, Long.: 47º38' W e Altitude: 546 m) e Natividade da Serra (Lat.: 23º23' S, Long.: 45º27'W e Altitude: 750 m), cujos valores médios anuais de chuva da série de dados de 30 anos (1970 a 1999) foram, respectivamente, 1.281,48; 1.344,07 e 1.570,23 mm ano^-1. Os parâmetros de entrada utilizados nas simulações foram, área de contribuição (Ac), em ha (1 ha), características geométricas do terraço com seção parabólica (B variando de 2 a 5 m, com intervalo de 0,1 m e H = 0,5 m), valor da velocidade de infiltração básica (VIB = 50 mm h^-1), comprimento do terraço (L, 300 m) e o Número da Curva (CN = 75).

Para se calcular a duração da chuva que levaria à obtenção das mesmas seções obtidas pelo modelo "Terradim", calculou-se o volume do escoamento superficial retido no terraço como o produto da seção pelo comprimento do terraço; obteve-se a lâmina de escoamento dividindo-se o volume pela área de contribuição do terraço; encontrou-se a altura da chuva pelo pelo método do Número da Curva e, finalmente, de forma iterativa, sua duração pelas equações de intensidade-duração-freqüência para as localidades de Araçatuba (Eq. 15), Piracicaba (Eq. 16) e Natividade da Serra (Eq. 17), obtidas no programa "Pluvio 1.3" (Pruski & Silva, 2000).

(15)

(16)

(17)

em que:

I - Intensidade média da chuva, mm h^-1;

T - período de retorno, anos; e

t - duração da chuva, min.

A Figura 3 apresenta o número de rompimentos do terraço e valores de período de retorno em função dos valores de B (m), para um valor de H fixado em 0,5 m, para as 3 localidades em questão.

Verifica-se, como esperado, que o número de rompimentos diminui à medida em que se adotam valores de B maiores. Para Araçatuba, por exemplo, o número de rompimentos variou de 30 a 2, para valores de B variando de 2 a 5 m. A ocorrência de trechos horizontais nos gráficos apresentados para determinada faixa de valores de B, significa que os terraços com essas dimensões apresentam o mesmo número de rompimentos dentro da série testada porém, caso se dispusesse de uma série com maior número de dados, poder-se-ia obter maior refinamento das curvas de número de rompimentos ou período de retorno em função de B.

Fixando-se um período de retorno em 5 anos obtém-se, para as localidades de Araçatuba, Piracicaba e Natividade da Serra, valores de B de 3,30, 3,34 e 3,80 m, respectivamente, cujos resultados estão coerentes, tendo em vista que as localidades apresentam pluviosidade crescente. As seções correspondentes a esses terraços são, respectivamente, 1,10, 1,11 e 1,27 m², e estão dentro da ordem de grandeza sugerida por Lombardi Neto (1994) que, para terraços em nível, varia de 0,60 a 1,50 m².

Realizando-se o cálculo inverso, ou seja, verificando-se que chuva levaria à adoção dessas seções, obtiveram-se as alturas de 88,81, 89,23 e 95,89 mm, para as localidades de Araçatuba, Piracicaba e Natividade da Serra, respectivamente. Para um período de retorno 5 anos, segundo as relações intensidade-duração-freqüência, as durações correspondentes são de 2,79, 2,81 e 3,33 h, respectivamente. Verifica-se, portanto, que são valores muito aquém da duração de 24 h e que sua adoção como critério levaria a superestimativas sensíveis da seção (8,19, 9,63 e 10,90 m², respectivamente).

CONCLUSÕES

1. A metodologia de avaliação de desempenho de estruturas hidráulicas por meio da análise de uma longa série de anos, pode também ser aplicada ao dimensionamento de terraços.

2. O modelo forneceu, como resultado, seções de ordem de grandeza coerentes com os dados da literatura e apresentou sensibilidade segundo o regime pluviométrico de cada localidade.

3. O critério da adoção da chuva de 24 h para o dimensionamento, se aplicado com o método do Número da Curva, geraria seções superdimensionadas.

4. A metodologia proposta permite o dimensionamento de terraços em nível em localidades que não dispõem ainda de relações intensidade-duração-freqüência de chuvas intensas.

LITERATURA CITADA

Protocolo 87 - 9/5/2003 - Aprovado em 1/3/2004

• Cruciani, D.E. A drenagem na agricultura. 4. ed. São Paulo: Nobel, 1987. 337p.
• Cruciani, D.E. Dimensionamento de sistemas de drenagem superficial e terraços com base nas características hidrológicas locais. In: Simpósio sobre Terraceamento Agrícola, Campinas : Fundação Cargill, 1989. 266p.
• Lombardi Neto, F.; Bellinazzi Júnior, R.; Lepsch, I.F.; Oliveira, J.B. de; Bertolini, D.; Galeti, P.A.; Drugowich, M.I. Terraceamento agrícola. Campinas: Secretaria da Agricultura e do Abastecimento do Estado de São Paulo, Coordenadoria da Assistência Técnica Integral, 1994. 39p. Boletim Técnico, 206
• Pruski, F. F. Desenvolvimento de metodologia para o dimensionamento de canais de terraços. Viçosa: Universidade Federal de Viçosa, 1993. 97p. Tese Doutorado
• Pruski, F.F.; Griebeler, N.P. Comparação entre métodos para a determinação do volume de escoamento superficial. In: Congresso Brasileiro de Engenharia Agrícola, 25 e Congresso Latino-Americano de Engenharia Agrícola, 2, 1996, Bauru. Anais... Bauru:Sociedade Brasileira de Engenharia Agrícola, 1996. p.422.
• Pruski, F.F.; Silva, D.D. Pluvio 1.3: Chuvas intensas para o Brasil. http://www.ufv.br/dea/gprh/pluvio 2000.
• Sparovek, G.; Silva, A. C. Dimensionamento hidrológico de terraços de drenagem e canais escoadouros. Revista da Universidade de Alfenas, Alfenas, v.3, p.137- 143, 1997.
• USA. Soil Conservation Service. National engineering handbook: IV. Hydrology. Washington: Department of Agriculture. 1972, n.p.

Datas de Publicação

• Publicação nesta coleção
  26 Jan 2005
• Data do Fascículo
  Dez 2004

Histórico

• Aceito
  01 Mar 2004
• Recebido
  09 Maio 2003