Simulation of hydraulic fracture utilizing numerical manifold method

A 2nd order numerical manifold method (NMM) based method is developed to simulate the hydraulic fractures propagating process in rock or concrete. The proposed method uses a weak coupling technique to analyze the fluid phase and solid phase. To study the seepage behavior of the fluid phase, all the fractures in solid are identified by a block cutting algorithm and form a flow network. Then the hydraulic heads at crack ends are solved. To study the deformation and destruction of solid phase, the 2-order NMM and sub-region boundary element method are combined to solve the stress-strain field. Crack growth is controlled by the well-accepted criterion, including the tension criterion or Mohr-Coulomb criterion for the initialization of cracks and the maximum circumferential stress theory for crack propagation. Once the crack growth occurs, the seepage and deformation analysis will be resolved in the next simulation step. Such weak coupling analysis will continue until the structure becomes stable or is destructed. Five examples are used to verify the new method. The results demonstrate that the method can solve the SIFs at crack tip and fluid flow in crack network precisely, and the method is effective in simulating the hydraulic facture problem. Besides, the NMM shows great convenience and is of high accuracy in simulating the crack growth problem.