Porous Media Equivalents for Networks of Discontinuos Fractures

1. Porous Media Equivalents for Networks of Discontinuos Fractures Equivalentes ao Meio Poroso para Redes de Fratures Descontínuas J. C. S. LONG, J.S. REMER, C.R. WILSON, P. A. WITHERSPOON 1982, Water Resources Research, 18(3), 645-658

2. 1. INTRODUÇÃO

3. Questão? Uma das mais importantes questões levantadas quando se considera o escoamento através das fraturas é: A rede de fraturas comporta-se como meio poroso?

4. Meio poroso? Aqui???

5. Pode-se modelar o sistema por um tensor de permeabilidade equivalente e proceder com a análise do movimento do fluido sob condições iniciais e de contorno conhecidas?

6. Snow [1965] • Fez estudos matemáticos de extensivos sistemas de fraturas. • Desenvolveu uma expressão matemática para o tensor de permeabilidade de uma fratura isolada de orientação e abertura arbitrárias em relação a um sistema de coordenadas. • Ele mostrou que o tensor de permeabilidade de uma rede de fraturas é formado por se adicionar os respectivos componentes do tensor de permeabilidade para cada fratura individual.

7. No campo constata-se que as fraturas tem comprimento finito. • Assim, cada fratura pode contribuir para a permeabilidade da rocha a medida que ela intercepta outras fraturas.

8. Uma fratura isolada, que não intercepta nenhuma outra fratura evidentemente não contribui para a permeabilidade da rocha. • Assim, o escoamento em dada fratura não é independente do escoamento em outras fraturas. • A permeabilidade da rocha então, não éa soma das permeabilidades de cada fratura.

9. 2. PERMEBAILIDADE HOMOGÊNEA ANISOTRÓRICA

10. Considerações • A lei de Darcy foi originalmente postulada para escoamento unidimensional. • A permeabilidade foi apresentada como um simples escalar.

11. Escoamento 3D | Anisotrópico • Ferrandon [1948], Collins [1961] propuseram: • Onde kij é simétrico e a matriz pode ser transformada em diagonal pela rotação do sistema de coordenadas.

12. Em rochas anisotrópicas, o gradiente hidráulico (ch/cl)e o vetor velocidade não são paralelos.Eles o são em rochas isotrópicas! C. W. Fetter – Applied Hydrogeology, página 155

13. Escoamento 3D | Anisotrópico • Assim Knão é um escalar em rochas como em rochas isotrópicas. • Em vez disso, há um tensor de condutividade (Kij) • Os valores do tensor dependem das direções x,y,z, os quais transformam os componentes do gradiente hidráulico nos de velocidade.

15. Escoamento 3D | Anisotrópico • Experimentos comprovam esta teoria. • Contudo, Collins [1961] destaca que não há garantia de que todo material poroso tenha um tensor de permeabilidade ortogonal.

16. 3. ESTATÍSTICA DA GEOMETRIA DAS FRATURAS

17. Sob um dado número de condições de contorno, o comportamento hidráulico de uma rocha fraturada com uma matriz impermeável é determinado inteiramente pela geometria do sistema de fraturas. • Para o escopo deste estudo a a descrição da geometria foi simplificada

18. Descrição • Cada fratura é individualmente descrita em termos de: • Abertura efetiva (ou hidráulica) • Orientação • Localização • Tamanho

19. Comportamento Hidráulico • É função da abertura efetiva. • A determinação da permeabilidade da fratura requer a determinação de sua abertura hidráulica. • Infelizmente, é muito difícil realizar testes hidráulicos em fraturas isoladas no campo. • Gale [1975] tentou fazer estes experimentos mas não foi bem sucedido.

20. 4. MÉTODO DE ANÁLISE NUMÉRICA

21. Um código numérico foi desenvolvido para gerar uma amostra do sistema de fraturas em duas dimensões usando as propriedades geométricas. • Por elementos finitos calcula-se Qg, o componente da vazão na direção do gradiente.

22. Usando a Lei de Darcy: • A condutividade hidráulica na direção do gradiente da amostra gerada é calculada por: • Onde: A é a área total perpendicular ao escoamento.

23. 5. GERAÇÃO DA MALHA

24. SET1 e SET2 INDEPENDENTES!

25. Distribuição Poisson Distribuição Normal Distribuição Log-Normal Distribuição Log-Normal

26. 6. MEDIÇÃO DA CONDUTIVIDADE HIDRÁULICA

27. A condutividade em um meio homogêneo pode ser definida na direção do escoamento ou na direção do gradiente. • Queremos medir a condutividade em um meio heterogêneo de rocha fraturada. • O gradiente pode ser aproximadamente linear.

28. A direção do escoamento é controlada pela direção das fraturas. • Um método que obtenha a condutividade na direção do gradiente precisa ser usado.

29. Kxx=Kg

30. Kxx=Kg

31. 7. ROTAÇÃO DA REGIÃO DO ESCOAMENTO

32. A condutividade pode ser medida em qualquer região escolhida.

33. Meio Homogêneo Anisotrópico PODE SER APROXIMADO COMO MEIO POROSO!

34. Meio Heterogêneo Anisotrópico NÃO PODE SER APROX. COMO MEIO POROSO!

35. 8. PROGRAMA“FRACTURE FLOW”

36. Elementos Finitos • Desenvolvido por Wilson [1970]. • As Fraturas são representadas por elementos lineares. • Escoamento obedecendo a leicúbica. • Matriz da rocha impermeável. • Estado permanente.

37. 9. VALIDAÇÃO DOMÉTODO NUMÉRICO

38. Exemplo 75 cm 110 cm 75 cm 110 cm

39. FLOW REGION 75 cm 75 cm

40. RESULTADOS 3 2 4 1

41. REDE 3 1,019x10-4 Q [cm³/s] 3 ZERO! 1,019x10-4 2 4 8,98x10-11 1 • Q2≠Q4 • Q=0 • ≠ Meio Poroso! ZERO! 5,42x10-10

42. Permeabilidade Direcional 75 cm 75 cm

43. Não há um tensor de Condutividade simétrico!

44. 10. EFEITO DA DENSIDADE DE FRATURA

45. 625 f.m-2 1111 f.m-2 1600 f.m-2 25 cm 40 cm 30 cm 25 cm 30 cm 40 cm

46. Número de interseções cresce de A para C Isopotenciais: Linear de A para C Aproxima-se de uma elipse de A para C

47. 11. EFEITOS DA ABERTURA E DA ORIENTAÇÃO

48. ABERTURA = LOG-NORMAL ORIENTAÇÃO = CONSTANTE ABERTURA = CONSTANTE ORIENTAÇÃO = NORMAL

50. 12. EFEITO DE ESCALA