A radial hydraulic fracture driven by a Herschel–Bulkley fluid

• Published in: Journal of non-Newtonian fluid mechanics Vol. 295; p. 104620
• Main Authors: Kanin, E.A., Dontsov, E.V., Garagash, D.I., Osiptsov, A.A.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.09.2021; Elsevier BV
• Subjects: Algorithms; Apexes; Approximate solution; Asymptotes; Chebyshev approximation; Constraining; Crack propagation; Fluid injection; Fracture mechanics; Herschel–Bulkley fluid; Hydraulic fracture; Hydraulic fracturing; Interpolation; Linear elastic fracture mechanics; Numerical methods; Parameters; Power law; Power-law rheology; Propagation; Quadratures; Radial (penny-shaped) model; Rheological properties; Rheology; Shear stress; Stress propagation; Yield strength; Yield stress; Hydraulic fracture; Power-law rheology; Herschel–Bulkley fluid; Radial (penny-shaped) model; Numerical methods; Approximate solution; Yield stress
• ISSN: 0377-0257; 1873-2631
• DOI: 10.1016/j.jnnfm.2021.104620

We analyse the influence of fluid yield stress on the propagation of a radial (penny-shaped) hydraulic fracture in a permeable reservoir. In particular, the Herschel–Bulkley rheological model is adopted that includes yield stress and non-linearity of the shear stress. The rock is assumed to be linear elastic, and the fracture is driven by the point source fluid injection with a constant volumetric rate. The fracture propagation condition follows the theory of linear elastic fracture mechanics, and Carter’s leak-off law is selected to govern the fluid exchange process between the fracture and formation. The numerical solution of the problem is found using the algorithm based on Gauss–Chebyshev quadrature and Barycentric Lagrange interpolation techniques. We also construct an approximate solution with the help of the global fluid balance equation and the near-tip region asymptote. The latter approximation is computationally efficient, and we estimate its accuracy by comparing the primary crack characteristics such as opening, pressure, and radius with those provided by the full numerical solution. We present examples corresponding to typical field cases and demonstrate that the addition of yield stress can lead to a shorter radius and wider opening compared to the corresponding case with simpler power-law fluid rheology. Further, we quantify the limiting propagation regimes (or vertex solutions) characterised by the dominance of a particular physical phenomenon. Relative to the power-law results, there are two new vertices that are associated with the domination of yield stress: storage-yield-stress and leak-off-yield-stress. To understand the influence of various problem parameters, we utilise the constructed approximate solution to investigate the dimensionless parametric space of the problem, in which the applicability domains of the limiting solutions are quantified. This enables one to quickly determine whether yield stress provides a strong influence for given problem parameters. •Analysis of the impact of fluid yield stress on radial fracture propagation.•Comparison with the power-law rheological model.•Numerical and rapid semi-analytical approximate solutions for the problem.•Construction of the dimensionless parameter space for the problem.•Quantitative identification of parameters for which yield stress is important.