Estado del arte de la simulación de la función de transferencia matriz - fractura en modelos de acuíferos con doble porosidad

Abstract

La función de transferencia del flujo entre matriz porosa y red de fracturas en acuíferos con doble porosidad ha sido presentada hace medio siglo (Barenblatt, 1960; Warren y Root, 1963) como el atributo más destacable a incorporar en este tipo complejo de medio subterráneo saturado. Entre las diversas propuestas que le sucedieron a la expresión analítica original, se reconocen dos grandes ramas: una formulación determinística (Zimmerman et al., 1993; Gerke y Van Genuchten, 1993) y, más reciente, una estocástica (Zhang, 2002). Todas ellas hacen depender la transferencia de flujo en forma proporcional a la diferencia de presiones en la región de contacto matriz – fisura. Se pretende así poner en evidencia la discontinuidad del medio saturado con doble porosidad y comprobar cómo influye esta componente, por ejemplo, en la interpretación de ensayos de bombeo en este tipo de acuíferos. Se avanzó primero con la revisión crítica de las distintas expresiones recopiladas y luego con la compilación de casos sencillos en los cuales se ha conseguido una integración del sistema de ecuaciones del flujo. Se realiza para cada caso el cómputo de una conductividad hidráulica equivalente y se discute la representatividad de este coeficiente. Finalmente, se enumeran conclusiones y posibles caminos de futuras investigaciones.

Matrix –fracture transfer function was introduced since fifty years ago (Barenblatt, 1960; Warren and Root, 1963) as the main characteristic to be considered in this particular saturated groundwater model. Afterwards, two branches were developed following that definition: a deterministic way (Zimmerman et al., 1993; Gerke and Van Genuchten, 1993) and a stochastic methodology (Zhang, 2002). Both of them treat the transfer flux as a function of the hydraulic pressure gradient across the contact region. Further objective was to verify the influence of this function in well test analysis. Case studies were compiled in order to establish a comparison among them with a global dual hydraulic conductivity coefficient. Finally, conclusions and future lines of investigations in this area are included.