Un modelo fractal para estimar la conductividad hídráulica no saturada de rocas fracturadas
Palavras-chave:
Conductividad hidráulica, rocas fracturadas, flujo no saturadoResumo
Durante las últimas dos décadas el estudio y la modelación del flujo de agua en rocas fracturadas no saturadas ha recibido considerable atención por parte de investigadores de distintas áreas como la geología, geofísica e ingeniería. Una de las razones principales para su estudio radica en la búsqueda de formaciones geológicas de baja permeabilidad para la construcción de repositorios de residuos nucleares. Para modelar el flujo de agua cuando se utiliza la aproximación del continuo es necesario contar con relaciones constitutivas del medio poroso considerado. Las relaciones constitutivas son curvas de saturación (S) y conductividad hidráulica (K) en función de la altura de presión (h). En general, la determinación en laboratorio de la curva de saturación no suele presentar mayores inconvenientes. En cambio, la determinación experimental de K(h) resulta difícil y costosa, por lo que se suele recurrir a modelos teóricos para su estimación. Entre los modelos teóricos más utilizados se destacan los de Burdine (Burdine, 1953) y Mualem (Mualem, 1976), los cuales permiten predecir la conductividad hidráulica K a partir de la curva de saturación S. Estos modelos predictivos han sido desarrollados para medios porosos de tipo sedimentario, no existiendo en la literatura modelos específicos para rocas fracturadas. En este trabajo se presenta un modelo fractal simple para predecir
la conductividad hidráulica no saturada de rocas fracturadas. La deducción del modelo se basa en la hipótesis de que la red de fracturas puede ser descrita mediante un objeto fractal clásico denominado carpeta de Sierpinski. La expresión propuesta de K es cerrada y depende únicamente de tres parámetros independientes: la dimensión fractal y las aperturas máxima y mínima de las fracturas. Una de las características que presenta el modelo propuesto, y que permite validarlo en forma teórica, es que la curva de conductividad hidráulica relativa es equivalente a la que se obtiene con el modelo de Burdine cuando se asume un factor de tortuosidad unitario. Finalmente, se realiza un estudio comparativo entre el modelo analítico propuesto y las relaciones constitutivas obtenidas mediante simulación numérica por Liu y Bodvarsson (2001). La comparación muestra que el modelo propuesto puede predecir dichas relaciones constitutivas en un amplio rango de valores de saturación de agua, y constituye por lo tanto una validación adicional. El modelo de K(h) propuesto es único en su tipo ya que ha sido derivado íntegramente a partir de leyes físicas clásicas y propiedades geométricas de la red de fracturas.
Downloads
Referências
Assouline S., 2001. A model for soil relative hydraulic conductivity based on the water retention characteristic curve. Water Resour. Res. 37 (2), 265–271.
Barton, C.A., Zoback, M.D., 1992. Self-similar distribution and properties of macroscopic fractures at depth in crystalline rock in the Cajon pass scientific drill hole. J. Geophys. Res. 97, 5181–5200.
Bear, J., 1988. Dynamics of Fluids in Porous Media, Elsevier, New York.
Berkowitz, B., Hadad, A., 1997. Fractal and multifractal measure of natural and synthetic fracture networks. J. Geophys. Res. 102 (B6), 205–218.
Buckingham, E., 1907. Studies on the movement of soil moisture. Washington D.C., U.S. Department of Agriculture. Bureau of Soil Bulletin 38, 29-61.
Burdine, N.T., 1953. Relative permeability from pore size distribution data. Trans. Am. Inst. Min. Metall. Pet. Eng. 198, 71 78.
Guarracino, L., 2006. A fractal constitutive model for unsaturated flow in fractured hard rocks. J. Hydrol. 324, 154–162.
Guarracino, L. Quintana, F., 2008. A constitutive model for flow in unsaturated fractured rocks. Hydrological Processes. In press.
Liu, H. H., Bodvarsson, G.S., 2001. Constitutive relations for unsaturated flow in fracture network. J. Hydrol. 252, 116–125.
Liu, H. H., Doughty, C., Bodvarsson, G.S., 1998. An active fracture model for unsaturated flow and transport in fractured rocks. Water Resour. Res. 12 (3), 513-522.
Mualem, Y., 1976. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 12 (3), 513–522.
Mandelbrot B. B., 1983. The Fractal Geometry of Nature. W. H. Freeman, Ney York, 488 pp.
Obuko, P.G., Aki, K., 1987. Fractal geometry in the San Andreas fault system. J. Geophys. Res. 92, 345–355.
Tuller M., Or D., 2002. Unsaturated hydraulic conductivity of structured porous media: a review of liquid configuration-based models. Vadose Zone Journal 1:14-37.
Turcotte D.L., 1986. Fractals and fragmentations. J. Geophys. Res.91, 1921–1926.
Turcotte D.L., 1997. Fractals and Chaos in Geology and Geophysics. Cambridge University Press, New York, 398 pp.
Tyler, S.W., Wheatcraft, S.W., 1990. Fractal process in soil water retention. Water Resour. Res. 26 (5), 1047–1054.
Downloads
Publicado
Como Citar
Edição
Seção
Licença
Copyright (c) 2010 Leonardo B. Monachesi, Luis Guarracino
Este trabalho está licenciado sob uma licença Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License.
A partir de 2022 (Vol. 43 número 2) os artigos serão publicados na revista sob uma licença Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International (CC BY-NC-SA 4.0)
De acordo com esses termos, o material pode ser compartilhado (copiado e redistribuído em qualquer meio ou formato) e adaptado (remixado, transformado e criado a partir do material outra obra), desde que a) a autoria e a fonte original de sua publicação (revista e URL da obra), b) não é utilizado para fins comerciais ec) os mesmos termos de licença são mantidos.
Antes desta data os artigos foram publicados na revista sob uma licença Creative Commons Attribution (CC BY)
Em ambos os casos, a aceitação dos originais pela revista implica a cessão não exclusiva dos direitos patrimoniais dos autores em favor do editor, que permite a reutilização, após edição (postprint), sob a licença que corresponde de acordo com a edição. .
Tal atribuição significa, por um lado, que após a sua publicação (postprint) na Revista GEOACTA da Associação de Geofísicos e Geodesistas, os autores podem publicar o seu trabalho em qualquer idioma, suporte e formato (nestes casos, solicita-se que fique registrado que o material foi publicado originalmente nesta revista); por outro, a autorização dos autores para que o trabalho seja colhido por SEDICI, o repositório institucional da Universidade Nacional de La Plata, e divulgado nas bases de dados que a equipe editorial considere adequadas para aumentar a visibilidade da publicação e seus autores.
Da mesma forma, a revista incentiva os autores para que após sua publicação na Revista da Associação de Geofísicos e Geodesistas, depositem suas produções em outros repositórios institucionais e temáticos, sob o princípio de que oferecer à sociedade produção científica e bolsa acadêmica irrestrita contribui para uma maior intercâmbio de conhecimento global.