Estudio Numérico de la Influencia de Fracturas en los Patrones de Flujo en Taludes

César Augusto Torres Paitan
Departamento de Engenharia Civil PUC Rio, Rio de Janeiro, Brasil,
cesar.torres0310@gmail.com

Eurípedes do Amaral Vargas Junior
Prof. Associado Departamento de Engenharia Civil PUC Rio, Rio de Janeiro, Brasil,
vargas@puc-rio.br

Raquel Quadros Velloso
Prof. Adjunto Departamento de Engenharia Civil UFOP, Ouro Preto,
rqvelloso@gmail.com

RESUMO:

Frequentemente, encostas na região sudeste do Brasil possuem camadas/depósitos (em geral de pouca espessura) de solos residuais ou coluvionares em contatos bruscos com maciços rochosos fraturados. As rochas que compõem estes maciços (em geral gnaisses ou granito-gnaisses) se encontram muitas vezes pouco alteradas, mas contendo famílias de fraturas, em geral com aberturas significativas o que lhes confere altas permeabilidades. Estas famílias de fraturas são muitas vezes sub-verticais, provavelmente de origem tectônica e associadas a sistemas das chamadas juntas de alívio de atitude aproximadamente paralela ao relevo. As análises de fluxo e processos de infiltração em tais encostas não levam em conta, de maneira geral, a existência destes sistemas de fraturas e é então desconhecida sua influência nos padrões de fluxo associados. O presente trabalho apresenta uma implementação numérica de fluxo em encostas onde a presença de fraturas é considerada. Detalhes da implementação numérica são apresentados além de resultados de casos hipotéticos analisados através dos quais é possível estabelecer a influência de fraturas subverticais no padrão de fluxo em encostas e em particular no campo de poropressões nas encostas.

PALAVRAS-CHAVE: Fluxo em encostas, análise numérica, fluxo em fraturas.

1. INTRODUÇÃO

As encostas na região sudeste do Brasil frequentemente apresentam problemas de instabilidade devidos comumente às precipitações pluviométricas de alta intensidade, elevando o nível d'água e assim um aumento nas poropressões na superfície potencial de ruptura. Estes taludes típicos destas regiões apresentam uma topografia acidentada, com formações de solos coluvionares e/ou residuais, e apresentam-se bruscamente em contato com rochas fraturadas. Geralmente, as análises de estabilidade de encostas não levam em conta o sistema de faturamento do embasamento rochoso, desconsiderando o comportamento hidráulico destas estruturas. Neste trabalho, propõe-se apresentar um sistema, composto por quatro programas computacionais, capaz de simular o comportamento do fluxo para estas formações (Paitan, 2013).

Como exemplos de aplicações são apresentadas análises de comportamento de fluxo de duas geometrias de taludes formados por solos residuais superpostos a um maciço rochoso que contém fraturas verticais e juntas de alívio paralelas ao relevo, típicas formações em taludes tropicais do Brasil.

2 FLUXO EM MEIOS POROSOS FRATURADOS

2.1 Modelo Conceitual

Existem modelos conceituais que podem ser empregados neste tipo de análise. O modelo frequentemente usado é o modelo de contínuo equivalente, que é aplicável a meios porosos e meios fraturados que apresentem alto grau de fissuramento. Outro modelo é o modelo de fraturas discretas, que tem como virtude levar em conta a localização espacial das fraturas, e é utilizado para redes de fraturas inseridas em matriz rochosa com permeabilidade intrínseca muito baixa. O modelo de dupla porosidade consiste de dois subdomínios, com propriedades hidráulicas diferentes, que interatuam e são acoplados por um termo de transferência que descreve a troca de fluido entre as duas regiões. Para o presente trabalho, será usado um modelo misto: um modelo contínuo de fluxo que leva em conta a geometria e distribuição das fraturas que estão inseridas no meio poroso.

2.2 Modelo Matemático

Os modelos matemáticos para fluxos em meios porosos fraturados seguem basicamente a equação de Richards, que se baseia na lei de Darcy e na equação de continuidade:

A equação (1) é uma forma modificada da equação de Richards (Cooley, 1983 e Huyakorn et al., 1984) e k _ij é o tensor de condutividade hidráulica, k _rw = k _rw ( S_w ) é apermeabilidade relativa do meio com relação ao grau de saturação de água S_w , y é a carga de pressão, z é a carga de elevação, q_s é teor de umidade volumétrica de saturação, e Q pode ser uma fonte ou um sumidouro. A equação (1) é altamente não linear e pode ser resolvida em termosy e S_w , sendo S_w = S_w (

2.3 Método Numérico de Solução

Na solução da equação (1) será aplicado o método dos elementos finitos (MEF) para discretizar o domínio, e o método de diferenças finitas para discretizar o tempo. Para maiores detalhes sobre a formulação do método recomenda-se revisar Neuman (1975), Zienkiewicz (1977) e Pinder & Gray (1977).

3 METODOLOGIA DE MODELAGEM

3.1 Geração das Fraturas

A geração das fraturas é feita com o programa FracGen 3D (Three-dimensional Fracture Generator) (Telles, 2006). O FracGen permite gerar fraturas tridimensionais nos modos determinístico e probabilístico. Na Figura 1 pode ser vista uma região com fraturas geradas de forma probabilística.

Figura 1. Famílias de fratura geradas no FracGen 3D de forma probabilistica.

3.2 Discretização do Domínio

No MEF, o domínio é dividido em pequenos subdomínios que se denominam elementos finitos. Neste trabalho, as geometrias planas serão discretizadas com elementos triangulares, e o espaço tridimensional será discretizado com elementos tetraédricos. Esta discretização gera uma malha compatível, na qual os nós dos elementos planos são coincidentes com os nós dos elementos de volume, como mostrado na Figura 2. Desta forma, os nós, apesar de pertencerem a elementos diferentes, compartilham as mesmas coordenadas e numeração, garantindo uma continuidade da carga hidráulica na interface entre os dois meios (Telles, 2006). Juanes & Samper (2000) indicam que este tipo de modelo tem sido implementado desde os anos setenta, e que foi adaptado por Kiraly (1987) e posteriormente generalizado por Perrochet (1995). Trabalhos relacionados a este modelo podem ser encontrados em Kiraly (1987), Telles (2006) e Alvarenga (2008).

Figura 2. Compatibilidade dos nós dos elementos representativos do meio poroso e dos elementos representativos das fraturas, Telles (2006).

3.2.1 Geração da Malha – ICEM CFD v. 14

Para a geração da malha é usado o programa ICEM CFD v.14 e foi necessário implementar um código de programação para a criação de um script que permita a geração automática de fraturas a partir das geometrias obtidas com o FracGen (Paitan, 2013) (Figura 3). Após a geração das fraturas no ICEM há a necessidade de se fazer uma edição destas.

Figura 3. Geração das fraturas no ICEM a partir do Script. Paitan (2013)

Nas Figuras 4 e 5, podemos ver a geometria inicial e a geometria depois da edição, respetivamente. A partir da malha criada, deve-se verificar e eliminar os erros. Uma vez feita as correções se verificará a qualidade da malha. O programa tem um verificador de qualidade e um suavizador de elementos (Figura 6).

Figura 4. Geometria inicial das fraturas.

Figura 5. Geometria após edição das fraturas.

Figura 6. Verificação da qualidade da malha gerada no

Após a verificação, a correção da malha de elementos finitos e a posterior otimização da largura de banda o passo seguinte é a exportação de dois arquivos com as informações geométricas da malha e as condições de contorno do problema em particular.

3.3 Solução Numérica e Visualização de Resultados

O programa de análise usado é o FTPF-3D Simulating Flow and Solute Transport in Porous and Fractured Media – Threedimensional (Telles, 2006). Este programa resolve as equações de fluxo, transporte de soluto e trajetória de partículas em meios porosos e fraturados, usando como código base o programa SWMS3D Code for Simulating Water Flow and Solute Transport in Three-Dimensional Variably Saturated Media (Šimuneket al., 1995). Para a visualização dos resultados usou-se o programa Pos3D - A Generic ThreeDimensional

Postprocessor (1991-2005) desenvolvido pelo TeCGraf-PUC/RIO – CENPES/PETROBRAS Convention.

Na etapa anterior obtiveram-se dois arquivos, e a partir destes gera-se, através deum programa desenvolvido (Paitan, 2013) um arquivo neutro (neutral file). O programa deanálise FTPF-3D usa como arquivo de entrada o arquivo neutro e como resultado da análise se obtém um arquivo de saída com os resultados de cargas de pressão, cargas totais e vetores de velocidades. Estes dados podem ser visualizadas graficamente usando o programa Pos3D.

4. EXEMPLOS APLICATIVOS

4.1 Analise de Fluxo Aplicado a um Talude com uma Fratura Vertical.

Neste exemplo é analisado um talude que compreende duas regiões: um solo residual e um meio poroso fraturado. O meio fraturado contém uma fratura vertical isolada. A análise de fluxo foi feita para condições de regime permanente de duas formas: uma primeira análise sem considerar a fratura e a segunda considerando a fratura. Nas Figuras 7 e 8 pode-se observar as geometria da fratura e do meio poroso. As características geométricas da fratura estão apresentadas na Tabela 1 e os parâmetros hidráulicos dos meios físicos estão descritos na Tabela 2.

Figura 7. Geometria contendo a fratura isolada e a malha de elementos finitos triangulares.

Figura8. A malha de elementos finitos tetraédricos dos meios e da fratura isolada com elementos finitos triangulares.

Tabela 1: Características geométricas da fratura.

Características geométricas	Fratura isolada
Numero de fraturas	1
Metodologia de geração	Determinística
Orientação:direçãode mergulho e mergulho	90º e 90º
Tamanho e forma da fratura	Tamanho 60 x 25 m. Polígono de 4 lados

Tabela 2: Parâmetros hidráulicos dos meios físicos.

Meios físicos	k(m/s)	S( m ⁻¹)
Meio poroso fraturado	5.79x10^-5	0.02
Fratura
	5.79x10^-4	4.5x10^-6

As condições de contorno em termo de carga de pressão podem ser verificadas na

Figura 9. Os resultados desta análise estão elaborados em termos de cargas totais, cargasde pressão e campo de velocidades, para ambos os domínios, isto é, para o talude com fratura e sem fratura. Quando se analisam estes modelos como um meio contínuo equivalente, isto é, sem modelar a fratura, não há como verificar a influência da fratura nas cargas de pressão que poderiam ocasionar, eventualmente, alguns problemas de instabilidade do talude.

Figura 9. Condições de contorno.

Pode-se observar nas figuras 10 a 17 o resultado diferenciado quando é levada em conta a fratura isolada. A fratura vertical isolada oferece um caminho preferencial ao fluxo e este caminho preferencial resulta na variação do campo de velocidades e sugere uma mudança nas poropressões (aumento da poropressão no pé do talude (Figura11)), indicando um possível problema de instabilidade de talude que não seria observada fazendo uma analise comumente feita sem considerar a fratura.

Figura 10. Distribuição das cargas de pressão no talude sem fratura.

Figura 11. Distribuição das cargas de pressão no talude com fratura.

Figura 12. Distribuição das cargas totais no talude sem fratura.

Figura 13. Distribuição das cargas totais no talude com fratura.

Figura 14. Campo de velocidades no talude sem fratura

Figura 15. Campo de velocidades no talude com fratura

Figura 16. Campo de velocidades na fratura e em no solo residual dos pontos A e B.

Figura 17. Campo de velocidades na fratura e no meio poroso nos pontos A e B, respectivamente.

4.2 Análise de Fluxo Aplicado a um Talude com uma Fratura Vertical e uma Junta de Alívio.

Neste exemplo, apresenta-se um talude de formação geológica típica do Rio de Janeiro. Este consiste em duas regiões: uma camada de solo residual e um meio poroso fraturado, um sobrepondo o outro. O meio poroso fraturado contém 1 fratura isolada vertical e 1 junta de alívio que é paralela à superfície de contato. A análise foi feita em condições de regime permanente. As Figuras 18 e 19 permitem observar a geometria das fraturas e dos meios porosos, respectivamente.

Figura 18. Geometria do exemplo contendo a fratura vertical isolada, a junta do alivio e a malha de elementos finitos.

Figura 19. A malha de elementos finitos dos meios porosos, a fratura e junta de alivio.

As características geométricas das fraturas e junta de alívio estão resumidas na Tabela 3.

Tabela 3: Características geométricas da fratura.

Os parâmetros hidráulicos dos meios físicos estão detalhados na Tabela 4.

Tabela 4: Parâmetros hidráulicos dos meios físicos.

Meio poroso	5.8x10^-5	0.02
Meio poroso fraturado	5.8x10^-8	0.02
Fratura e junta de alivio	1.6x10^-3	4.5x10^-6

As condições de contorno deste caso de estudo podem ser verificadas na Figura 20.

Figura 20. Condições de contorno. Os resultados da análise de fluxo permanente estão expressos em termos de cargas de pressão, cargas totais e campo de velocidades, nas Figuras 21 a 25.

Dos resultados deste exemplo pode-se observar que, eventualmente, estes padrões de fluxo poderiam ocasionar problemas de instabilidade da camada de solo residual. O fluxo neste exemplo hipotético de encosta se mostra seguindo o caminho que apresenta menor resistência que é o que oferecem as fraturas, e é fácil verificar que poderiam ocorrer problemas de piping associados a um enfraquecimento do solo e sua ruptura posterior na zona de contato entre o solo residual e o maciço rochoso, que é o plano de fraqueza desta formação típica.

Figura 21. Distribuição das cargas de pressões no

Figura 22. Distribuição das cargas totais no talude.

Figura 23. Distribuição das cargas totais da fratura e junta de alivio.

Figura 24. Campo de velocidades do talude.

Figura 25. Campo de velocidades da fratura e junta de alivio.

4. CONCLUSÕES

O presente trabalho apresentou a criação de um sistema formado por quatro programas que em conjunto são capazes de modelar fluxo em meios fraturados e meios porosos fraturados, em condições de fluxo permanente ou transiente, em condições saturadas e não saturadas. Este sistema pode ser aplicado também para analisar padrões de fluxo em outras estruturas geológicas e para outras condições como escavações, solo de fundação para barragens, etc.

Nas análises de encostas típicas do Rio de Janeiro foi considerado um regime permanente de águas subterrâneas, com a condição de que a camada de solo residual tenha uma região não saturada. Dos exemplos hipotéticos analisados, pode-se observar a notável influência das fraturas no comportamento hidráulico em virtude da variação dos valores da poropressão e das velocidades de fluxo no meio. Normalmente, nas análises de fluxo não se leva em conta estas estruturas devido a sua complexidade para inserir sua geometria, características e propriedades dentro do meio contínuo. A abordagem de um modelamento que considera estas descontinuidades permite verificar o estado no qual se encontram estas típicas formações, e qualquer outro talude, e conjuntamente com metodologias adequadas de instrumentação poderiam levar a um melhor monitoramento dos problemas de deslizamentos e erosão de taludes naturais ou artificias.

REFERÊNCIAS

Alvarenga, J. E. M. (2008)Modelagem Numérica do Transporte de Vírus em Aqüiferos Fraturados-Porosos. Dissertação de Mestrado. Departamentode Engenharia Civil, Pontifícia Universidade Católica do Rio de Janeiro, Rio de Janeiro, 257pp.

Cooley, R. L. (1984) Some New Procedures for Numerical Solution of Variable Saturated Flow Problems. WaterResour. Res., 19(5), p. 1271 - 1285.

Huyakorn, P. S.; Thomas, S. D.; Thompson, B. M. (1984) Techniques for Making Finite Elements Competitive in Modeling Flow in Variably Saturated Porous Media.WaterResour.Res., 20(8),p.1099 - 1115.

Juanes, R.; Samper, J. (2000). Una Formulación General y Eficiente de Las Fracturas en el MEF: I. Aspectos Teóricos. Revista Internacional de Métodos Numéricos para Cálculo y Diseño emIngenieria, volume 16, 4, p.471 - 491.

Kiraly, L. (1987). Large Scale 3D Groundwater Flow Modelling in Highly Heterogeneous Geologic Medium, Groundwater Flow and Quality Modelling, E Custodio et al. (Eds), Dordrecht,Holland, D. Reidel Publishing Co. p. 761 – 775.

Mualem, Y. A. New Model to Predict The Hydraulic Conductivity of Unsaturated Porous Media. Water Resour. Res., 12: 513–522.

Neuman, S. P. (1975).Galerkin Approach to Saturated-Unsaturated Flow in Porous Media, Chapter 10 in Finite Elements in Fluids, Vol. I Viscous Flow andHydrodynamics, Edited by R. H. Gallagher, J. T. Oden, C. Taylor, and O.C. Zienkiewicz., p. 201-217, Jhon Wiley and Sons, London.

Paitan, C. A.T. (2013). Modelagem Numérica de Fluxo em Meios Fraturados e Meios Porosos Fraturados.

Dissertação de Mestrado. Departamento de Engenharia Civil, Pontifícia Universidade Católica do Rio de Janeiro, Rio de Janeiro, 107p.

Perazzolo L. (2003). Estudo Geotécnico de dois Taludes da Formação Serra Geral, RS. Dissertaçãode Mestrado. Departamento de Engenharia Civil, Universidade Federal do Rio Grande do Sul, Porto Alegre, 150p.

Perrochet, P. (1995). Finite Hyperelements: A 4D Geometrical Framework Using Covariant Bases and Metric Tensors, Communications in Numer, Meth. Engng., Vol. 11, N 6, p. 525-534.

Pinder, G. F.; & Gray, W. G. (1977).Finite Element Simulation in Surface and Subsurface Hydrology.Academic Press, New York.

Šimůnek, J.; Huang, K.; Van Genuchten, M. T.

(1995).The SWMS_3D Code for Simulation Water Flow and Solute Transport in Three-Dimensional Variably-Saturated Media.U.S. Salinity LaboratoryAgricultural Research Service, n°139, version 1.0.

Telles, I. A. (2006).Desenvolvimento de um Sistema Integrado para Modelagem de Fluxo e Transporte em Meios Porosos e Fraturados. Tese deDootorado. Departamento de Engenharia Civil, Pontifícia Universidade Católica do Rio de Janeiro, Rio de Janeiro, 164p.

Therrien, R.; Sudicky, E.A. (1996).Three-dimensional Analysis of Variably-Saturated Flow and Solute Transport in Discretely-fractured Porous Media.Journal of Contaminant Hydrology.23, p. 1 - 44.

Van Genuchten, TH. M. (1980). A Closed-form Equation for Predicting the Hydraulic Conductivity of Unsatured Soils, Soil Science Society of America Journal, volume 44, no 5.

Zienkiewicz, O. C. (1977).The Finite Element Method.3rd ed., McGraw-Hill, London, UnitedKingdom.

Articulo