This paper develops a closed-form approximate solution for a penny-shaped hydraulic fracture whose behaviour is determined by an interplay of three competing physical processes that are associated with fluid viscosity, fracture toughness and fluid leak-off. The primary assumption that permits one to construct the solution is that the fracture behaviour is mainly determined by the three-process multiscale tip asymptotics and the global fluid volume balance. First, the developed approximation is compared with the existing solutions for all limiting regimes of propagation. Then, a solution map, which indicates applicability regions of the limiting solutions, is constructed. It is also shown that the constructed approximation accurately captures the scaling that is associated with the transition from any one limiting solution to another. The developed approximation is tested against a reference numerical solution, showing that accuracy of the fracture width and radius predictions lie within a fraction of a per cent for a wide range of parameters. As a result, the constructed approximation provides a rapid solution for a penny-shaped hydraulic fracture, which can be used for quick fracture design calculations or as a reference solution to evaluate accuracy of various hydraulic fracture simulators.

No abstract

This study revisits the problem of a steadily propagating semi-infinite hydraulic fracture in which the processes of toughness-related energy release, viscous dissipation and leak-off compete on multiple length scales. This problem typically requires the solution of a system of integro-differential equations with a singular kernel, which is complicated by the need to capture extremely disparate length scales. In this study the governing equations are rewritten in the form of one non-singular integral equation. This reformulation enables the use of standard numerical techniques to capture the complete multiscale behaviour accurately and efficiently. This formulation also makes it possible to approximate the problem by a separable ordinary differential equation, whose closed-form solution captures the multiscale behaviour sufficiently accurately to be used in practical applications. We also consider a similar reformulation of the equations governing the propagation of a buoyancy-driven semi-infinite hydraulic fracture. The resulting numerical solution is able to capture the near-tip multiscale behaviour efficiently and agrees well with published solutions calculated in the large-toughness limit.

No abstract

The goal of this study is to analyse the steady flow of a Newtonian fluid mixed with spherical particles in a channel for the purpose of modelling proppant transport with gravitational settling in hydraulic fractures. The developments are based on a continuum constitutive model for a slurry, which is approximated by an empirical formula. It is shown that the problem under consideration features a two-dimensional flow and a boundary layer, which effectively introduces slip at the boundary and allows us to describe a transition from Poiseuille flow to Darcy's law for high proppant concentrations. The expressions for both the outer (i.e. outside the boundary layer) and inner (i.e. within the boundary layer) solutions are obtained in terms of the particle concentration, particle velocity and fluid velocity. Unfortunately, these solutions require the numerical solution of an integral equation, and, as a result, the development of a proppant transport model for hydraulic fracturing based on these results is not practicable. To reduce the complexity of the problem, an approximate solution is introduced. To validate the use of this approximation, the error is estimated for different regimes of flow. The approximate solution is then used to calculate the expressions for the slurry flux and the proppant flux, which are the basis for a model that can be used to account for proppant transport with gravitational settling in a fully coupled hydraulic fracturing simulator.

Please cite this article as: Dontsov, E.V., Peirce, A.P., An enhanced pseudo-3D model for hydraulic fracturing accounting for viscous height growth, non-local elasticity, and lateral toughness, Engineering Fracture Mechanics (2015), doi: http://dx. AbstractThe goal of this paper is to develop a more accurate pseudo-3D model for hydraulic fracturing. The primary weaknesses of the classical pseudo-3D model are: 1) its inability to capture viscous height growth, and 2) its failure to include lateral fracture toughness. These flaws are addressed respectively by: 1) introducing an apparent fracture toughness in the vertical direction, and 2) using an approximate non-local elasticity operator. To evaluate the accuracy and the level of improvement of the model we have developed, the results are compared to the predictions calculated using a recently developed fully planar hydraulic fracturing simulator.

A hybrid discrete-continuum numerical scheme is developed to study the behavior of a hydraulic fracture crossing natural fractures. The fully coupled hybrid scheme utilizes a discrete element model for an inner domain, within which the hydraulic fracture propagates and interacts with natural fractures. The inner domain is embedded in an outer continuum domain that is implemented to extend the length of the hydraulic fracture and to better approximate the boundary effects. The fracture is identified to propagate initially in the viscosity-dominated regime, and the numerical scheme is calibrated by using the theoretical plane strain hydraulic fracture solution. The simulation results for orthogonal crossing indicate three fundamental crossing scenarios, which occur for various stress ratios and friction coefficients of the natural fracture: (i) no crossing, that is, the hydraulic fracture is arrested by the natural fracture and makes a T-shape intersection; (ii) offset crossing, that is, the hydraulic fracture crosses the natural fracture with an offset; and (iii) direct crossing, that is, the hydraulic fracture directly crosses the natural fracture without diversion. Each crossing scenario is associated with a distinct net pressure history. Additionally, the effects of strength contrast and stiffness contrast of rock materials and intersection angle between the hydraulic fracture and the natural fracture are also investigated. The simulations also illustrate that the level of fracturing complexity increases as the number and extent of the natural fractures increase. As a result, we can conclude that complex hydraulic fracture propagation patterns occur because of complicated crossing behavior during the stimulation of naturally fractured reservoirs.

No abstract