Appendix D – Fractured Rock Appendix
1.0 Introduction
The behavior of LNAPL in fractured bedrock is not necessarily intuitive and is not as easily described using the principles and techniques adopted for more homogeneous porous media. Differences have been observed and theorized for how LNAPL behaves in fractured bedrock formations. While the same physical principles apply for multiphase flow in fractured aquifers as in porous aquifers, unique characteristics of finite and restricted fluid flow paths can lead to unexpected results in fractured settings. Notable differences from unconsolidated materials include:
- Secondary porosity effects (bedrock fracture size and orientation) provide preferential flow pathways for LNAPL migration following a release.
- Fluid flow can be limited to fracture planes, thereby practically constricting aquifer bulk volume. As a result, LNAPL gradients can be greater than in porous media, resulting in increased lateral or vertical migration.
- Initial LNAPL head pressures focused on finite planar features can exceed those exerted on a more voluminous porous media aquifer and result in LNAPL propagation deep into the saturated zone via vertical and/or sub-vertical fractures (Adamski 2005); (Hardisty 2010).
- Fluctuations in groundwater elevation may cause “pumping” of LNAPL that result in greater lateral and vertical spread of LNAPL for longer time frames than for unconsolidated materials (Hardisty 2010).
- Similar to confined or perched LNAPL, the gauged thickness of LNAPL in a monitoring well may not accurately represent the mobile LNAPL interval in the formation.
- LNAPL penetration into the primary rock matrix may preferentially occur along fractures. Following removal of LNAPL from a fracture, there is the potential for LNAPL within the adjacent primary matrix to re-enter the fracture over time.
Approaches to characterization and specific investigation methods for fractured bedrock sites are provided in the ITRC Fractured Bedrock guidance (ITRC 2017a) and other documents. The purpose of this Appendix is to provide information that is focused on the understanding, approaches and methods that uniquely apply to LNAPL releases in fractured bedrock conditions.
2.0 Fractured Rock Concepts and Technical Background
Fractured bedrock can represent a wide variety of rock types with different primary matrix and secondary fracture porosity and permeability influencing the hydrogeologic behavior of groundwater and contaminants in bedrock. To simplify the discussion of how LNAPL behaves in various rock types, LNAPL behavior is discussed for four categories of fractured rock (Nelson 2001):
- Type 1 – minimal matrix permeability and porosity results in fluid flow and storage in the fractures (Single porosity).
- Type 2 – low matrix permeability but high porosity results in fluid flow in fractures, but storage in both the rock matrix and fractures (Dual porosity).
- Type 3 – high matrix permeability and porosity results in fluid flow and storage in both the rock matrix and fractures, but permeability is enhanced by the fractures (Dual permeability).
- Type 4 – high matrix permeability and porosity results in fluid flow and storage in both the rock matrix and fractures, with the fractures not significantly increasing flow (Single porosity or equivalent porous media analogue).
These simplified schematics (Figure FR-1) are intended to present conceptualized system behaviors, but the practitioner is cautioned to develop site-specific LNAPL Conceptual Site Models (LCSMs) as the range of conditions and behaviors is infinitely variable due to the unique setting of each site.
2.1 Fluid Entry Pressure
For LNAPL to migrate within fractured rock, the fluid entry pressure of a fracture must be overcome. This concept is equivalent to pore entry pressure for porous media, or the Type 4 category described above. Fracture entry pressure is a function of the difference in pressure between LNAPL and the other fluid(s) present, as well as the wetting angle, interfacial tension between the fluids, and the aperture of the fracture (i.e., the Laplace Equation) (Corey 1990).
Where:
Pe = Entry pressure
σ = Interfacial tension between LNAPL and water
b = Fracture aperture
= Contact angle with the wetting fluid
In the following sections, the concept of LNAPL head is directly related to this entry pressure, whether it is where LNAPL first enters a fractured rock system, or as it migrates through subsurface fractures of variable size (Kueper and McWhorter 1991).
2.2 Secondary Porosity Effects
LNAPL migration following a release in a fractured rock setting is driven primarily by fracture size (aperture) and orientation. However, where significant primary porosity and permeability of the bedrock matrix exists, fluid behavior similar to that which occurs in porous media can occur, resulting in what is commonly termed a dual porosity system. The range of behaviors is simplified into the Type 1 through Type 4 rock type hierarchy, with Type 1 units demonstrating the greatest departure from porous media fluid behavior and Type 4 units most closely following porous media behavior. The effect of variability in primary matrix permeability and porosity for the rock types is illustrated in Figure FR-1.
Figure FR-1. Rock type and characteristics.
LNAPL migration in Type 1 bedrock (e.g., crystalline bedrock), where primary matrix porosity is negligible, will allow for minimal LNAPL entry into the primary matrix. LNAPL will migrate preferentially along larger more transmissive fractures and can be hydraulically up-gradient from the release location. Increased porosity and permeability will result in increased entry of LNAPL into the rock matrix. For Type 2 and 3 bedrocks (e.g., limestone, sandstone), LNAPL will preferentially flow along fractures, but will increasingly penetrate into the primary matrix. Type 4 bedrock (e.g., highly weathered) can behave more similarly to unconsolidated formations. The effect of variability in primary matrix permeability and porosity on LNAPL distribution is illustrated in Figure FR-2.
Figure FR-2. Rock type and LNAPL distribution.
2.3 Influence of Fracture Plane Angle on LNAPL Migration
The primary rock properties controlling LNAPL penetration and distribution in rock types 1 through 3 are fracture spacing and orientation. Where the volume of LNAPL released to the rock exceeds the volume of porosity near the spill, these properties will control LNAPL penetration and distribution. More steeply dipping fracture planes will allow for greater vertical penetration and lower fracture density (or larger spacing between fractures or lower “secondary porosity” relative to the bulk aquifer volume) will also promote deeper LNAPL penetration by focusing LNAPL head on those more discrete and finite preferential flow pathways. Fracture continuity and pervasiveness, which can be a function of weathering or depth, figure into both flow path direction and the continuum from Type 1 to Type 4 behavior.
In cases where shallow dip angle fracture or bedding planes are present above the water table and vertical fracturing is absent, lateral LNAPL migration will dominate the post-release plume distribution. The migration of the LNAPL body will predominantly occur in the direction of a fracture dip and follow fracture topography. Higher fracture dip angles, lower primary porosity, and open fractures can result in conditions that are not typically observed in unconsolidated formations (i.e., deep LNAPL penetration).
Figure FR-3. LNAPL lateral migration.
In cases where near-vertical fracture attitudes are present, LNAPL can be driven deep into the saturated zone due to the limited pathways for vertical migration and resultant buildup of LNAPL head. LNAPL may penetrate below the apparent groundwater elevation where the angle and aperture of a fracture, combined with the pressure head of the LNAPL (i.e., the height of the LNAPL source in the fracture above groundwater), can result in groundwater within the fracture being displaced to a significant depth (Hardisty et al. 2003). LNAPL has been observed at depths exceeding 30 feet (10m) below the water table in fractured rock settings.
Figure FR-4. LNAPL vertical migration.
2.4 Influence of Groundwater Fluctuations on LNAPL Distribution
Even after reduction or elimination of an LNAPL driving head (i.e., post-release and stabilization), groundwater elevation fluctuations in fractured bedrock have the potential to cause both vertical and lateral redistribution of LNAPL (Hardisty et al. 2003). These effects are most critical in Type 1 and Type 2 bedrock where fracture orientation and frequency are within the zone of groundwater fluctuation.
As groundwater elevations fall, LNAPL will tend to migrate vertically downwards through sub-vertical transmissive fractures via gravity drainage. As groundwater elevations rise, LNAPL can migrate along any fracture that is no longer water saturated—including those with shallow or sub-horizontal dips. The result is that the fracture network could effectively act as a “pump” driving LNAPL in unexpected directions, including apparently up gradient, and may result in LNAPL being “trapped” potentially at significant depths below the static water level.
Figure FR-5. Changes in LNAPL distribution with groundwater fluctuation (modified from (Hardisty et al. 2003).
The difficulty in predicting these movements, combined with the inherent difficulty of constructing monitor wells that intersect the main fractures controlling LNAPL storage and movement, mean that establishing LNAPL body stability may be challenging. This requires a robust LCSM, longer-term monitoring data, and use of multiple lines of evidence.
As with porous media frameworks, changing water levels or saturations over time can also lead to further isolation of LNAPL bodies by alternately allowing for drainage from LNAPL-connected pores (or pathways) or trapping of LNAPL in pore or aperture openings too small for non-wetting fluid penetration absent a driving head—resulting in isolated LNAPL ganglia. These fluctuations can also be more pronounced in fractured rock aquifers with karst conditions being the most extreme (CL:AIRE 2014).
2.5 Gauged LNAPL Thickness in Fractured Bedrock
The gauged thickness of LNAPL in groundwater monitoring wells and open boreholes in fractured bedrock may result in exaggerated LNAPL thicknesses as compared to the mobile interval of LNAPL. These effects are most likely to occur in Type 1 and Type 2 bedrock and can result in confined or perched conditions resulting in exaggerated gauged thicknesses relative to what may be observed in unconsolidated formations (Figure FR-6).
There are several tools that may help to clarify these conditions for the LCSM. After well installation, use of diagnostic gauge plots and evaluation of constant LNAPL discharge during LNAPL baildown testing may help with understanding of the location of transmissive LNAPL fractures and the presence of perched or confined fractures. However, unlike unconsolidated formations, evaluation of diagnostic gauge plot data and LNAPL baildown testing can be complicated by 1) the presence of multiple transmissive fractures, each with different hydrogeologic conditions; and 2) the potential for both lateral and vertical redistribution of LNAPL with groundwater elevation changed as discussed above. If these conditions are expected prior to drilling, some techniques are available prior to well installation, such as NAPL Flexible Liner Underground Technologies (FLUTe™) systems, sampling with packers during drilling, etc.
Figure FR-6. LNAPL vertical distribution and apparent monitoring well thickness.
These conditions illustrate why it is critical to develop the LCSM concurrent with the investigation process, to avoid connecting fractures that may not otherwise be connected. A well connecting two, previously unconnected near horizontal fractures, could cause LNAPL to migrate downward to lower fractured zones.
2.6 Re-Entry of LNAPL into Fractures from Primary Matrix
Entry of LNAPL into the rock matrix can occur where the primary matrix porosity and permeability is sufficiently high and is most likely to occur in Type 2, 3, and 4 bedrocks. In situations where the secondary porosity and permeability is higher than the primary porosity and permeability, mobile LNAPL will preferentially flow along the fractures, but may penetrate into the primary matrix over time as shown on the left side of the Figure FR-7. LNAPL recovery efforts will more readily remove the mobile LNAPL within the higher transmissivity fractures; however, mobile LNAPL may also remain in the primary matrix. After LNAPL recovery is stopped, re-entry of mobile LNAPL from the primary matrix back into the fractures may occur.
Figure FR-7. LNAPL re-entry into fracture from primary matrix.
Conceptually, this poses similar challenges as “back diffusion” of dissolved phase constituents in fractured rock environments. In addition, LNAPL recovery rates and transmissivities may decline during active LNAPL recovery as mobile LNAPL is removed from the secondary porosity, but may rebound if re-entry of LNAPL into the fractures occurs.
3.0 LNAPL Conceptual Site Model
Development of an LNAPL Conceptual Site Model (LCSM) is described in detail in Section 4 of the main body of this document. The level of detail required to develop an LCSM is in part a function of the complexity of environmental conditions. Although multiphase fluid relationships and flow in fractured rock settings are generally governed by the same physical principals as in unconsolidated media, more extreme ranges in quantification of physical properties are found in fractured rock, resulting in increased variability and complexity. Examples include fracture aperture (analogous to porosity) and bulk hydraulic conductivity, where ranges over multiple order of magnitude are common. In addition, the spacing, pervasiveness, amount of “interconnectedness” and attitudes of fractures (idealized as two-dimensional planes) imparts another layer of complexity to fractured rock fluid flow. The result is that an increased level of detail is typically required to develop an LCSM in fractured bedrock.
The ITRC Fractured Bedrock guidance (ITRC 2017a) and DNAPL guidance (ITRC 2011) provide a detailed discussion of the methods to develop a Conceptual Site Model in fractured bedrock. For LNAPL releases in fractured rock, practitioners are encouraged to review both ITRC guidance documents. Other tools for assessing fracture systems have been evaluated by the USGS in their Fractured Rock Geophysical Toolbox Method Selection Tool (FRGT-MST).
Introduction of an LNAPL with sufficient driving head into fractured bedrock often results in: an LNAPL plume distribution with anisotropic dimensions aligned with dip direction in the unsaturated zone, and aligned with bedding strike in the saturated zone; as well as vertical penetration to underlying beds proportional to bed thickness and/or cross-bed fracturing. In addition, LNAPL penetration into primary matrix pore space would be expected in these units and would increase with time and increased primary pore space of the various interbedded units.
Development of an LCSM should include conducting in-depth literature reviews and employing a multiple lines of evidence approach to developing and continually improving LCSMs to match empirical findings. Conventional investigation techniques (i.e., vertical borings and monitoring wells) are typically ill-suited to develop LCSMs in fractured bedrock. However, numerous alternative tools are available and can be used to develop LCSMs in fractured conditions. As real-world fractured aquifer systems exhibit behaviors along a continuum of the simplified Type 1 to Type 4 scenarios, a hybrid LCSM is often the best fit for managing site data and upon which future management decisions can be based.
Consistent with Section 4 of this guidance on the development of the LCSM, the following questions are intended to guide LCSM development related to LNAPL in fractured rock. If all or part of a question is immaterial to matrix lithology, then the reader is directed to Section 4. Only those aspects of the questions that are specific to fractured rock are further discussed.
1. Is current and future land use known? ▼Read more
2. Does the potential for preferential pathways exist? ▼Read more
3. How does stratigraphy relate to affecting impacts and potential migration? ▼Read more
4. Is the extent of the source areas known? ▼Read more
5. Are dissolved or vapor issues expected based on LNAPL composition? ▼Read more
6. Are dissolved or vapor plumes characterized? ▼Read more
7. Do soil or groundwater concentrations exceed criteria? ▼Read more
8. Are exposure pathways complete or incomplete? ▼Read more
9. Is the LNAPL body stable? ▼Read more
10. Is mobile LNAPL hydrogeologic condition known? ▼Read more
11. What “tools” are available to develop an LCSM in fractured Bedrock? ▼Read more
4.0 Remedial Technologies
Remediation in fractured bedrock is subject to the same constraints as porous media, with the added challenge that explicitly mapping the transport pathways (fractures) and mass distribution is more difficult, or in some cases, impossible. This is one of the primary reasons that pump-and-treat remediation is commonly selected in complex fractured rock systems.
A key limitation in fractured bedrock remediation is access to the target compounds and distribution of reagents, when injection-based in-situ remedies are applied. The second matrix in Appendix A summarizes data requirements and data collection methods to evaluate remedial alternatives at fractured rock sites.
An important factor regarding the efficacy and cost for treatment, particularly in source areas, is the degree to which mass has diffused into the matrix and the potential significance of “rebound” that may occur due to back diffusion. When COCs are contained primarily in the fractures, remediation can have significant benefits in a relatively short period of time and pump and treat can be applied as a remedy, not just a containment measure.
In cases where a significant mass of recalcitrant compounds has diffused into the matrix, access to the source material can be inefficient, commonly resulting in rebound after in-situ treatment, or extended duration of pump and treat.
The reverse-diffusion effect on remedial timeframe may not be as significant for contaminants that may degrade in the matrix (e.g., aromatic hydrocarbons) if conditions are conducive for decay in the matrix or upon reverse-diffusion back into the fractures. Mature plumes of recalcitrant compounds, however, often require the most aggressive measures of remediation, such as thermal treatment or chemical oxidation to significantly reduce mass. Table FR-2 summarizes the applicability of conventional remedial technologies and provides a relative ranking as to their cost, reliability and period of performance.
Table FR-2. Summary information for LNAPL remediation technologies in fractured bedrock
View Table FR-2 in Adobe PDF format.