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, which is characterized by relatively low hydraulic conductivities, has usually low flow velocities. 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, which provides usually a high porosity and can store substantial amounts of contaminants. As a consequence, the concentration in the fracture decreases. The continuous diffusion from the fractures into the matrix slows the propagation speed of the plume more and more down. Once diffused into the matrix, it is difficult to remove the contaminant again. Due to the extremely low flow velocities in the limestone matrix, the removal of contaminants from the matrix with a pump-and-treat remediation happens mainly due to back-diffusion from the matrix into the fractures, which can take very long. Part of the contaminant can be sorbed to the surface of the limestone matrix (and the fracture walls), which has a retarding effect on plume propagation. The sorption behavior is often quantified by sorption coefficients or retardation 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 The basic equations for flow and transport in fractured porous media are shown below. For more details, be referred to standard literature like e.g. Bear, J. (1972) . 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_w [/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_w [/math] water sources/sinks [math] q_{darcy} [/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, also degradation [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