Boundary element methods such as the displacement discontinuity method (DDM) are promising for describing the geo-mechanics of growing fractures as they reduce the dimensionality of the problem by one. Still, this advantage is overshadowed by the computational power required to model stage scale simulations of hydraulic fractures primarily due to the nonlinearity arising from fluid flow inside the fracture. Typical assumptions to improve the computational efficiency lead to the inability of simulators to describe proppant transport which is critical to treatment design. In this paper, we present a novel algorithm to increase the speed of convergence of coupled geomechanics (DDM) and fluid flow problem. The algorithm employs an analytical, local linear approximation of the Reynolds lubrication equation to solve the coupled non-linear system of equations. An initial solution is estimated through the linear system, and it is then used in the iterative solution approach for the non-linear system as an initial guess. This improvement has allowed us to use the algorithm to develop a fully three-dimensional hydraulic fracturing simulator which can run in an acceptable runtime using commonly available computational power.

This content is only available via PDF.
You can access this article if you purchase or spend a download.