Two-dimensional reactive transport modeling of the alteration of a fractured limestone by hyperalkaline solutions at Maqarin (Jordan)

Abstract

Two-dimensional reactive transport modeling of the Maqarin Eastern Springs site, a natural analogue for the alteration of a fractured limestone by high-pH Portland cement waters, has been performed using the CrunchFlow code. These 2D calculations included transport by advection–dispersion–diffusion along a single fracture and diffusion in the wall rock. Solute transport was coupled to mineral dissolution and precipitation. A limited sensitivity analysis evaluated the effect of different values of primary mineral surface areas, flow velocity and sulfate concentration of the inflowing high-pH solution.Major secondary minerals include ettringite–thaumasite, C–S–H/C–A–S–H and calcite. C–S–H/C–A–S–H precipitation is controlled by the dissolution of primary silicates. Ettringite precipitation is controlled by diffusion of sulfate and aluminum from the wall rock to the fracture, with aluminum provided by the dissolution of albite. Calcite precipitation is controlled by diffusion of carbonate from the wall rock. Extents of porosity sealing along the fracture and in the fracture-wall rock interface depend on assumptions regarding flow velocity and composition of the high-pH solution. The multiple episodes of fracture sealing and reactivation evidenced in the fracture infills were not included in the simulations. Results can qualitatively reproduce the reported decrease in porosity in the fractures and in the wall rock next to the fractures. Instances of porosity increase next to fractures caused by carbonate dissolution were not reproduced by the calculations.