An implicit level set algorithm for hydraulic fracturing with a stress-layer asymptote

The capability to simulate a hydraulic fracturing process is an essential tool that can be used to optimize treatment design and increase the efficiency of field operations. In most practical cases, hydraulic fractures propagate in a multi-layered rock formation. As a result, there is a need to incorporate the effect of such heterogeneities in fracturing models to achieve an accurate prediction. To capture the layered structure of rocks, a hydraulic fracture simulator typically requires a fine mesh, which leads to a drastic reduction in computational performance. An alternative is to use more sophisticated models that are capable of providing reasonably accurate predictions even on a relatively coarse mesh. In the case of fracture growth modeling, the pivotal component of the simulation is a fracture front tracking algorithm that accounts for the layered structure of the formation. Consequently, this paper aims to extend the established Implicit Level Set Algorithm (ILSA) to account for the effect of multiple stress layers within the tip asymptote. The enhanced front tracking algorithm involves the stress-corrected asymptote that incorporates the influence of stress layers within the near-tip region. To further increase the validity region of the stress-corrected asymptote, the stress relaxation factor is introduced, and its accuracy is examined. The numerical algorithm is validated against the reference semi-analytical solutions as well as experimental observations. In addition, we investigate the sensitivity of the fracture geometry to mesh size to demonstrate that the front tracking algorithm based on the stress-corrected asymptote retains its accuracy on a coarse mesh. •A numerical algorithm for a planar hydraulic fracture is presented.•Tip asymptotic solution is used to determine the evolution of fracture boundary.•The near-tip logic includes the effects of stress layers.•A level set algorithm is used to track the moving fracture front.•The developments improve the overall accuracy of the algorithm.