DETERMINAÇÃO DA CONDUTIVIDADE HIDRÁULICA E DA SORVIDADE DE UM SOLO NÃO-SATURADO UTILIZANDO-SE PERMEÂMETRO A DISCO

DETERMINATION OF UNSATURATED HYDRAULIC CONDUCTIVITY AND SORPTIVITY OF A SOIL USING A DISK PERMEAMETER

Resumos

Dois métodos, utilizando permeâmetros a disco, foram usados para medir a sorvidade e a condutividade hidráulica em três horizontes de um solo de Santa Maria da Boa Vista, Pernambuco, utilizando-se potenciais de fornecimento de água de 0, -2,5, -5 e -10 cm de água. No primeiro método, a condutividade hidráulica foi obtida através de estimativas da sorvidade e do fluxo estacionário, utilizando-se apenas um permeâmetro. No segundo, as medidas da sorvidade e da condutividade hidráulica foram feitas utilizando-se valores dos fluxos estacionários, obtidos a partir de dois permeâmetros de diferentes raios. Esse segundo método apresentou resultados mais consistentes dos tempos -- gravitacional e geométrico --, e do raio característico de poros que o primeiro método. A sorvidade foi tão importante quanto a condutividade hidráulica para caracterizar o processo de infiltração. A redução brusca dos valores dos raios característicos de poros do horizonte A2/B, com mudanças no potencial de fornecimento de água, revelou a heterogeneidade da estrutura deste horizonte, permitindo identificar a influência dos volumes argilosos compactados sobre o processo de infiltração.

permeâmetro; infiltração

DETERMINAÇÃO DA CONDUTIVIDADE HIDRÁULICA E DA SORVIDADE DE UM SOLO NÃO-SATURADO UTILIZANDO-SE PERMEÂMETRO A DISCO¹

EDVANE BORGES², ANTONIO CELSO DANTAS ANTONINO³, ATTILIO DALL'OLIO⁴, PIERRE AUDRY⁵ e CLEMENTE JOSÉ GUSMÃO CARNEIRO⁶

RESUMO - Dois métodos, utilizando permeâmetros a disco, foram usados para medir a sorvidade e a condutividade hidráulica em três horizontes de um solo de Santa Maria da Boa Vista, Pernambuco, utilizando-se potenciais de fornecimento de água de 0, -2,5, -5 e -10 cm de água. No primeiro método, a condutividade hidráulica foi obtida através de estimativas da sorvidade e do fluxo estacionário, utilizando-se apenas um permeâmetro. No segundo, as medidas da sorvidade e da condutividade hidráulica foram feitas utilizando-se valores dos fluxos estacionários, obtidos a partir de dois permeâmetros de diferentes raios. Esse segundo método apresentou resultados mais consistentes dos tempos ¾ gravitacional e geométrico ¾, e do raio característico de poros que o primeiro método. A sorvidade foi tão importante quanto a condutividade hidráulica para caracterizar o processo de infiltração. A redução brusca dos valores dos raios característicos de poros do horizonte A2/B, com mudanças no potencial de fornecimento de água, revelou a heterogeneidade da estrutura deste horizonte, permitindo identificar a influência dos volumes argilosos compactados sobre o processo de infiltração.

Termos para indexação: permeâmetro, infiltração.

DETERMINATION OF UNSATURATED HYDRAULIC CONDUCTIVITY AND SORPTIVITY OF A SOIL USING A DISK PERMEAMETER

ABSTRACT - Sorptivity and hydraulic conductivity were determined by two methods using disc permeameters. Infiltration experiments were carried out on three horizons of a representative soil of the semi-arid region, located in Santa Maria da Boa Vista, State of Pernambuco, Brazil. Water was supplied to this soil by disc permeameters at potentials of 0, -2.5, -5 and -10 cm of water. In the first method, hydraulic conductivity was measured using estimates of sorptivity and steady state fluxes. In the second, sorptivity and hydraulic conductivity were determined using measurements of steady state fluxes with two disc permeameters of different radius. Results of the second method, regarding characteristic times and characteristic pore radius, were more consistent than those of the first one. Sorptivity was as important as hydraulic conductivity to characterize the infiltration process. The abrupt reduction of the characteristic pore radius in the horizon A2/B with changes in the water supply potential indicates the structure heterogeneity of this horizon and allows the identification of the role of the clay compact volumes on the infiltration process.

Index terms: permeameter, infiltration.

INTRODUÇÃO

O processo de infiltração é caracterizado por dois parâmetros: a sorvidade (S₀) e a condutividade hidráulica (K₀). O valor de S₀ é associado à fase inicial da infiltração. Nessa fase, a infiltração é praticamente independente dos efeitos gravitacionais e geométricos. Os valores S₀ e K₀ na camada superficial do solo definem a partição dos aportes externos de água (precipitação, irrigação), afetando assim tanto o armazenamento de água no solo, como o escoamento superficial. É importante considerar que os macroporos são funcionais quando o solo está muito próximo à saturação, sendo capazes de transportes preferenciais de água e substâncias químicas de forma rápida durante o processo de infiltração. Por causa da natureza frágil dos macroporos é necessária a realização de medidas de infiltração in situ usando técnicas que minimizem as perturbações do solo. Do ponto de vista mais prático, pode-se avaliar a importância dos macroporos sobre a infiltração em condições de inundação (Watson & Luxmoore, 1986) através da diferença entre S₀ e K₀ a um potencial matricial y₀, próximo à saturação, e S e K medidos à saturação.

O processo de infiltração pode ser estudado utilizando-se dispositivos chamados permeâmetros a disco. Técnicas baseadas nesses permeâmetros são de grande interesse, pois permitem determinar simultaneamente, S₀ e K₀in situ, de forma simples e rápida, minimizando as perturbações no solo. A análise das relações entre S₀ e K₀ pode ser usada para inferir a importância relativa das forças capilares, da gravidade e da geometria da fonte de água em controlar a densidade de fluxo de água que entra no solo, permitindo caracterizar o seu comportamento hidráulico. A baixa penetração da frente de umedecimento permite caracterizar o solo em pequenas espessuras (camadas de superfície, interfaces entre dois horizontes).

Inicialmente, os permeâmetros a disco foram usados para medir as taxas de infiltração mantendo-se um potencial positivo e constante na superfície do solo (Talsma, 1969). Posteriormente, foram desenvolvidos dispositivos que permitem manter esse potencial constante e negativo (Dirksen, 1975; Clothier & White, 1981), possibilitando a determinação das propriedades físico-hídricas de solos fissurados. Perroux & White (1988) aperfeiçoaram esse dispositivo, a fim de se obter um maior controle do potencial negativo mantido na base do disco.

Dois métodos são usados para determinar S₀ e K₀ a partir de medidas das taxas de infiltração, obtidas com permeâmetros a disco, tanto nos instantes iniciais desse processo, onde se aplica a teoria de Philip (1969), como no regime estacionário, empregando a equação proposta por Wooding (1968). O primeiro método de Perroux & White (1988) é baseado no uso de um único permeâmetro, e requer estimativas da sorvidade e do fluxo estacionário. O fato de esse método estimar a sorvidade, a partir do comportamento inicial da infiltração acumulada na razão da raiz quadrada do tempo, representa uma desvantagem na obtenção deste parâmetro, sobretudo em solos úmidos e altamente permeáveis. O segundo, de Smettem & Clothier (1989), necessita dos valores dos fluxos estacionários, medidos a partir de dois permeâmetros de raios diferentes situados em locais tão próximos quanto possível, visando diminuir o efeito da variabilidade espacial do solo.

O objetivo deste trabalho foi comparar dois métodos de determinação dos valores de sorvidade e de condutividade hidráulica, com o uso de permeâmetros a disco, em três horizontes de um solo do Sertão do Estado de Pernambuco.

MATERIAL E MÉTODOS

Teoria

No início da infiltração, quando a capilaridade é dominante, tem-se que a infiltração acumulada em função do tempo (I(t)) é determinada pela equação (Philip, 1957):

onde S₀ representa a sorvidade (cm/s^1/2).

De acordo com Wooding (1968), o fluxo de água que passa através do permeâmetro a disco de raio r, tende para o valor estacionário:

é a equação transformada de Kirchoff, ou o potencial de fluxo matricial (Gardner, 1958); K_0, a condutividade hidráulica correspondente ao potencial da água y_0;y_n, o potencial da água do solo antes de ser iniciado o experimento, e q_n e q₀ as umidades volumétricas inicial e final.

Assumindo-se a expressão exponencial para D(q), relacionada com a sorvidade S₀, sugerida por Reichardt et al. (1972):

onde b e g são parâmetros que podem ser obtidos através de ajustes por mínimo quadrado, e integrando-se a Equação 3, obtém-se

com

onde b é um parâmetro de forma da difusividade (1/2£b£p/4), que foi considerado igual a 0,55 (White & Sully, 1987).

Durante o processo de infiltração tridimensional, os fatores que podem afetar o fluxo de entrada de água no solo são: a geometria da fonte de água, a capilaridade e a gravidade. Esses fatores são caracterizados pelas escalas de tempo ¾ gravitacional (t_grav) e geométrico (t_geom) ¾ e de comprimento ¾ comprimento capilar (l_c) e raio característico de poros (l_m).

Os tempos gravitacional e geométrico são determinados respectivamente pelas equações (Thony et al., 1991):

O tempo gravitacional pode ser interpretado como o intervalo durante o qual a dinâmica da infiltração é controlada mais pelas forças capilares do que pela ação da gravidade; enquanto o tempo geométrico representa o intervalo durante o qual o caráter unidimensional do processo de infiltração não é alterado. Talsma (1969) observou que a influência da gravidade pode ser considerada desprezível para t £ 0,02t_grav. Em estudos realizados no campo, foram obtidos valores de t_grav variando entre 0,08 a 34 horas (White & Sully, 1987). Segundo Philip (1969), o intervalo de tempo no qual a Equação 1 pode ser aplicada deve ser menor que t_geom.

A escala de comprimento capilar é dada por White & Sully (1987):

Essa escala representa a importância relativa das forças capilares em relação à gravidade, quando a água é transmitida de uma fonte mantida a q₀ através do solo, com umidade inicial q_n.

A partir da teoria da capilaridade, Philip (1987) definiu um raio característico de poros:

onde s é a tensão superficial, r é a massa específica da água e g a aceleração da gravidade. O raio característico de poros define a dimensão média dos poros que participam do processo de infiltração submetido à pressão aplicada y₀. Quanto maior for o raio característico, maior será o efeito da gravidade em comparação ao da capilaridade.

Experimentos

Os experimentos de infiltração foram realizados utilizando dois permeâmetros, P1 e P2, de concepção idêntica, representados esquematicamente na Fig. 1, diferindo apenas com relação às dimensões de suas bases, 25 cm e 8 cm de diâmetro, respectivamente.

O modo de operação foi o mesmo para os dois permeâmetros. O reservatório lateral foi preenchido com água, através do tubo C₁, utilizando-se uma seringa hipodérmica. A pressão imposta y₀ é dada pela diferença entre z₂ e z₁ (Fig. 1), de tal forma que

A altura z₂ é fixa, sendo de 1 cm para P1 e 5 cm para P2.

Para o preenchimento do reservatório de alimentação, a base do permeâmetro a disco foi imergida em uma bacia cheia de água e em seguida foi feita a aspiração do ar pelo tubo acoplado na rolha de borracha B₀ (Fig. 1). Quando o nível da água alcançou o topo do reservatório, a pinça de Mohr (PM) foi fechada. O ar que permaneceu na base do disco foi retirado com movimentos que permitiam que as bolhas escoassem pelo reservatório de alimentação (R_A). A leitura da variação do nível de água no reservatório de alimentação em determinados intervalos de tempo, permitiu a obtenção da infiltração acumulada. As observações foram interrompidas quando a variação do nível de água no reservatório, em função do tempo, tornou-se constante, ou seja, quando o regime estacionário foi alcançado.

Inicialmente, as superfícies escolhidas para a realização dos experimentos foram limpas, para retirar qualquer tipo de vegetação presente, e niveladas, eliminando assim as irregularidades mais proeminentes. Para uma adequada conexão entre a base do permeâmetro e a superfície do solo, foi necessária a colocação de um material de contato. Neste caso, utilizou-se uma fina camada (5 mm) de areia, umedecida a um teor volumétrico de água de 3% a 5%, como sugerido por Thony et al. (1991), e distribuída uniformemente na superfície do solo. Para a determinação da umidade inicial (q_i) foram coletadas amostras de solo a uma distância de aproximadamente 30 cm da borda do permeâmetro. Para a determinação de q₀, amostras do solo úmido foram coletadas, imediatamente após a remoção do permeâmetro e da camada de areia, no próprio local onde estava assentado o permeâmetro.

As medidas para determinação de S₀ e K₀ foram efetuadas em um solo situado no município de Santa Maria da Boa Vista, Sertão de Pernambuco. O solo é classificado como Podzólico "intergrade" planossolo, e apresenta uma sucessão de três horizontes, que possuem, respectivamente, as seguintes características: A, humífero, arenoso e homogêneo; B, de tipo Podzólico, argilo-arenoso, com estrutura poliédrica fina, pouco desenvolvida, homogêneo; e A2/B, profundo, também globalmente argilo-arenoso, mas de estrutura heterogênea, constituída de volumes mais argilosos, compactos e freqüentemente endurecidos (B), inseridos em uma rede arenosa, porosa, lavada e embranquecida, sem consistência (A2). Os percentuais de areia, silte e argila nesses horizontes são apresentados na Tabela 1.

Nas profundidades de 10, 70 e 180 cm, consideradas representativas de cada horizonte, foram realizados experimentos de infiltração com a água a potenciais que variaram de 0 até -10 cm. Para a determinação de S₀ e de K₀, em cada uma dessas profundidades, foram utilizados tanto o método de Perroux & White (1988), como o de Smettem & Clothier (1989). No primeiro método, o valor de S₀ foi estimado pelo ajuste, pelo método dos mínimos quadrados, da Equação 1. No segundo caso, K₀ e S₀ foram determinados através das Equações 12 e 13, resultantes da resolução simultânea de duas equações da forma da Equação 2, a partir das medidas dos fluxos estacionários q₁ e q₂ em cada permeâmetro:

RESULTADOS E DISCUSSÃO

A Tabela 2 apresenta os valores de S₀, K₀, l_m, t_grav e t_geom obtidos nas três profundidades do solo, consideradas representativas dos horizontes A, B e A2/B, pelos métodos de Perroux & White (PW) e Smettem e Clothier (SC), nos diferentes potenciais de infiltração (y₀) aplicados. Em relação ao método SC, quando os valores de y₀ foram mais negativos, observou-se uma diminuição da condutividade hidráulica, da sorvidade e do raio característico de poros em cada uma dessas profundidades. No método PW, observou-se um aumento dos valores dos raios característicos de poros quando os valores de y₀ tornaram-se mais negativos; este resultado é fisicamente incoerente pois os poros de maiores dimensões participam mais do processo de infiltração quanto menos negativo for o potencial y₀. Por outro lado, os resultados do método SC apresentam correta relação entre os raios característicos de poros e o potencial y₀, sugerindo a adoção preferencial deste último método. Além disso, a análise comparativa entre os valores de S₀ e K₀ mostra que no método SC existem relações monótonas crescentes, de K₀ em relação a S₀, bem definidas em todos os horizontes, enquanto no método PW tais relações são menos claras e até em sentido oposto. O horizonte A2/B apresentou um comportamento diferenciado em relação aos horizontes A e B.

O método SC permitiu evidenciar, para cada horizonte, uma boa variação na condutividade hidráulica nas proximidades imediatas da saturação (Tabela 2). Por exemplo, para o horizonte A a condutividade hidráulica à saturação (y₀=0) foi da ordem de nove vezes maior que a condutividade hidráulica para y₀ = -10,0 cm.

Os valores dos raios característicos de poros obtidos pelo método SCvariaram de 0,31 a 0,99 mm nos três horizontes estudados. Esses valores estão praticamente na faixa entre 0,10 e 0,90 mm, que corresponde aos valores obtidos por outros autores (White & Sully, 1987; Sauer et al., 1990; Thony at al., 1991; Cook & Boeren, 1994).

Em virtude dos altos valores de t_grav nas profundidades de 10 e 70 cm, pode-se afirmar que, neste solo, a sorvidade é tão importante quanto a condutividade hidráulica para caracterizar o processo de infiltração. Considerando-se que, para um experimento de infiltração vertical, o regime estacionário quase sempre é obtido após oito a dez vezes o t_grav, seria necessário prosseguir o experimento durante várias horas e mesmo alguns dias para obter a estabilização do fluxo e, portanto, a condutividade hidráulica. Os valores de t_geom, muito inferiores a t_grav, mostram que o caráter tridimensional imposto pela geometria da fonte ao escoamento resulta no estabelecimento de um regime estacionário muito antes que a gravidade desempenhe um papel importante, possibilitando a realização de experimentos mais rápidos que os tradicionais de infiltração vertical. Isto pode ser verificado na Fig. 2, que representa a evolução do fluxo com o tempo para um conjunto de três repetições na profundidade de 70 cm, com y₀ = -5 cm, utilizando-se o disco P1. Para este caso, o tempo gravimétrico foi de 208 min e o geométrico de 37 min. O regime estacionário foi alcançado de forma muito rápida em, aproximadamente, t_geom/9.

A comparação entre os resultados obtidos nas profundidades de 10 e 70 cm, respectivamente nos horizontes A e B, mostram uma diminuição da condutividade menor do que era esperada pelo aumento da fração argilosa. Tal fato pode ter sido causado pela porosidade interagregados, relacionada com a microestrutura do horizonte B, que pode contribuir de maneira importante para o aumento da velocidade de infiltração, como mostram os valores do raio característico de poros.

O comportamento diferenciado do horizonte A2/B, profundidade de 180 cm, pode ser interpretado referindo-se à estrutura já assinalada: 1) nos potenciais menos negativos, de 0 a -2,5 cm, de fornecimento da água, a porosidade interagregados, isto é, a rede arenosa, embranquecida e porosa, parece ser também a responsável por uma infiltração aparentemente mais rápida em relação às características granulométricas globais, já que se obteve um maior raio característico de poro e um menor t_grav nessa profundidade; 2) a partir do potencial de fornecimento de água de -5 cm, pela participação dos volumes argilosos compactados no processo, tem-se uma diminuição importante do raio característico de poros, acarretando aumento da sorvidade e do t_grav, e a predominância das forças capilares no processo de infiltração, contrariamente ao observado em relação aos potenciais de 0 e -5 cm, onde as forças gravitacionais apresentam uma maior influência sobre o processo.

CONCLUSÕES

1. O método de Smettem & Clothier apresenta melhores resultados nas determinações da sorvidade e da condutividade hidráulica que o método de Perroux & White.

2. A sorvidade é tão importante quanto a condutividade hidráulica para caracterizar o processo de infiltração.

3. A redução brusca dos valores dos raios característicos de poros do horizonte A2/B, com mudanças no potencial de fornecimento de água, revela a heterogeneidade da estrutura desse horizonte, permitindo identificar a influência dos volumes argilosos compactados sobre o processo de infiltração.

REFERÊNCIAS

• CLOTHIER, B.E.; WHITE, I. Measurement of sorptivity and water diffusivity in the field. Soil Science Society of America. Journal, Madison, v.45, p.241-245, 1981.
• COOK, F.J.; BROEREN, A. Six methods for determining sorptivity and hydraulic conductivity with disc permeameters. Soil Science, Baltimore, v.157, p.3-11, 1994.
• DIRKSEN, C. Determination of soil water diffusivity by sorptivity measurements. Soil Science Society of America. Proceedings, Madison, v.39, p.22-27, 1975.
• GARDNER, W.R. Some steady-state solutions of unsaturated moisture flow equation with application to evaporation from a water table. Soil Science, Baltimore, v.85, p.228-232, 1958.
• PERROUX, K.M.; WHITE, I. Designs for disc permeameters. Soil Science Society of America. Journal, Madison, v.52, p.1205-1215, 1988.
• PHILIP, J.R. The quasi-linear analysis, the scattering analog, and other aspects of infiltration and seepage. In: FOK, Y.S. (Ed.). Infiltration development and application Honolulu: Water Resources Research Center, 1987. p.1-27.
• PHILIP, J.R. Theory of infiltration. Advances in Hydroscience, San Diego, v.5, p.215-297, 1969.
• PHILIP, J.R. The theory of infiltration: 4. Sorptivity and algebraic infiltration equations. Soil Science, Baltimore, v.84, p.257-265, 1957.
• REICHARDT, K.D.; NIELSEN; D.R.; BIGGAR, J.W. Scaling of horizontal infiltration into homogenous soils. Soil Science Society of America. Proceedings, Madison, v.36, p.241-245, 1972.
• SAUER, T.J.; CLOTHIER, B.E.; DANIEL, T.C. Surface measurements of the hydraulic properties of tilled and untilled soil. Soil and Tillage Research, Amsterdam, v.15, p.359-369, 1990.
• SMETTEM, K.R.J.; CLOTHIER, B.E. Measuring unsaturated sorptivity and hydraulic conductivity using multiple disc permeameters. Journal of Soil Science, Oxford, v.40, p.563-568, 1989.
• TALSMA, T. In situ measurements of sorptivity. Australian Journal of Soil Research, Melbourne, v.7, p.269-276, 1969.
• THONY, J.L.; VACHAUD, G.; CLOTHIER, B.E.; ANGULO-JARAMILLO, R. Field measurement of the hydraulic properties of soil. Soil Technology, Amsterdam, v.4, p.111-123, 1991.
• WATSON, K.W.; LUXMOORE, R.J. Estimating macroporosity in a forest watershed by use of a tension infiltrometer. Soil Science Society of America. Journal, Madison, v.50, p.578-582, 1986.
• WHITE, I.; SULLY, M.J. Macroscopic and microscopic capillary length and time scales from field infiltration. Water Resources Research, v.23, p.1514-1522, Washington, 1987.
• WOODING, R.A. Steady infiltration from a shallow circular pond. Water Resources Research, Washington, v.4, p.1259-1273, 1968.

Datas de Publicação

• Publicação nesta coleção
  21 Nov 2005
• Data do Fascículo
  Nov 1999