Physical processes and governing equations

Highlights
Description of major physical processes in limestone aquifers Governing equations for flow and transport in limestone

Physical processes Flow in fractured limestone can be very complex, because strongly contrasting flow conditions between fractures and matrix can prevail. The flow in the fractures can be very fast, whereas the limestone matrix has usually low flow velocities, which is characterized by relatively low hydraulic conductivities. However, the fate of contaminants in fractured limestone aquifer is determined by both, the fractures and the matrix. Transport of substances in porous media can be subdivided in advective transport (due to the groundwater flow) and diffusive/dispersive transport (related to concentration gradients). The fractures are the primary transport pathways for contaminants, where mainly advective transport dominates. The hydraulic conductivity of the fractures is mainly depending on their aperture (width). For an estimation of the flow velocity in fractures, the cubic law is often used as approximation. With that, the water flux scales with the third power of the fracture aperture. Hence, fractures with a big aperture are far more important as small fractures, since the total water flow and the flow velocity velocity are greater in them. The connectivity of the fractures is another very important parameter for the spreading of a substance in the aquifer. While being transported through the fractures, the contaminant diffuses into the surrounding matrix. Once diffused into the matrix, the contaminant is difficult to get out again, since the matrix provides usually a high porosity and can store substantial amounts of contaminants. Due to the extremely low flow velocities in the limestone matrix, the removal of contaminants from the matrix happens again due to diffusion back to the fractures. Part of the contaminant can be sorbed to the surface of the limestone matrix (and the fracture walls). The sorption behavior is often quantified by sorption coefficients. Furthermore, contaminants can be degraded by microbial activity, if the proper microbes are present and the conditions are favorable (for PCE: anaerobic conditions). Governing equations Flow Groundwater flow in porous media can be described by Darcy's law, giving a relation between hydraulic head gradient and groundwater flow, with the hydraulic conductivity as proportionality factor. The flow field is calculated by solving the mass balance for incompressible fluid flow in combination with Darcy's law. Mass balance equation for incompressible fluid flow: [math] S_s \frac{\partial h}{\partial t} - \nabla \cdot (K \nabla h) = q [/math] Darcy's law (also used in the fractures): [math] q_{darcy} = -k \nabla h [/math] The conductivities in the fractures can be several orders of magnitude higher than the conductivity of the limestone matrix. The fracture conductivity is usually calculated via the fracture aperture using the cubic law, which gives a relation between fracture aperture and fracture conductivity. Cubic law to calculate fracture conductivity: [math] k_f = \frac{b^2 \rho_w g}{12 \mu} [/math] Variable Name [math] S_s [/math] specific storage [math] t [/math] time [math] K [/math] hydraulic conductivity [math] q [/math] sources/sinks [math] q [/math] Darcy flux [math] k_f [/math] fracture conductivity [math] b [/math] fracture aperture [math] \rho_w [/math] density of water [math] g [/math] gravity constant Contaminant transport The advection-dispersion equation is a balance equation for the transport of a substance in a porous medium. The first term describes the storage (here including retardation due to linear sorption), the second and third terms describe advective and dispersive transport of a substance. Degradation can be included as additional source/sink term. Transport equation for a dissolved species: [math](n + \rho_b k_d) \frac{\partial c}{\partial t} + \nabla \cdot n(vc) - \nabla \cdot n(D_m \nabla c) = q_m[/math] Transport in fractures: [math] b \frac{\partial c_f}{\partial t} + \nabla \cdot b(v_f c_f) - \nabla \cdot b(D_f \nabla_T c_f) = q_f[/math] Variables are: Variable Name [math] n [/math] porosity [math] q_m [/math] mass sources and sinks [math] c_f [/math] concentration in the fractures [math] \rho_b [/math] bulk density [math] k_d [/math] sorption coefficient [math] c [/math] concentration [math] v [/math] flow velocity in the matrix [math] v_f [/math] flow velocity in the fractures [math] D_m [/math] dispersion coefficient Return to Content