Propagation regimes of buoyancy-driven hydraulic fractures with solidification

Abstract

This study investigates the propagation of a semi-infinite buoyancy-driven hydraulic fracture in situations when the fluid is able to solidify along the crack walls. Such problems occur when hot magma ascends from a chamber due to buoyancy forces and solidifies by interacting with colder rock. In the model, the solidification rate is calculated assuming a one-dimensional heat transfer problem, in which case it becomes mathematically equivalent to Carter’s leak-off model, which is commonly used to describe the fluid leak-off from a hydraulic fracture into a porous rock formation. In order to construct a mathematical model for a buoyancy-driven hydraulic fracture with solidification, the aforementioned thermal problem is combined with (i) linear plane-strain elasticity to ensure equilibrium of the rock surrounding the fracture, (ii) linear elastic fracture mechanics to determine the fracture propagation, (iii) lubrication theory to capture the viscous fluid flow inside the crack and to account for the effect of buoyancy, and (iv) volume balance of the magma. To address the problem, the governing equations are first rewritten in terms of one integral equation with a non-singular kernel, which significantly simplifies the analysis and the procedure for obtaining a numerical solution. The latter solution is shown to obey a multiscale behaviour near the fracture tip that is fully resolved by the numerical scheme. In order to understand the structure of the solution and to quantify the regimes of propagation (and the associated transitions), a thorough analysis of the problem has been performed. Finally, the developments are applied to investigate the non-steady propagation of a buoyancy-driven fracture that is fed by a constant flux.