Recent examples of hydraulic fracture diagnostic data suggest complex, multi-stranded hydraulic fractures geometry is a common occurrence. This reality is in stark contrast to the industry-standard design models based on the assumption of symmetric, planar, bi-wing geometry. The interaction between pre-existing natural fractures and the advancing hydraulic fracture is a key condition leading to complex fracture patterns. Performing hydraulic fracture design calculations under these less than ideal conditions requires modeling fracture intersections and tracking fluid fronts in the network of reactivated fissures. Whether a hydraulic fracture crosses or is arrested by a pre-existing natural fracture is controlled by shear strength and potential slippage at the fracture intersections, as well as potential debonding of sealed cracks in the near-tip region of a propagating hydraulic fracture. We present a complex hydraulic fracture pattern propagation model based on the Extended Finite Element Method (XFEM) as a design tool that can be used to optimize treatment parameters under complex propagation conditions. Results demonstrate that fracture pattern complexity is strongly controlled by the magnitude of anisotropy of in situ stresses, rock toughness, and natural fracture cement strength as well as the orientation of the natural fractures relative to the hydraulic fracture. Analysis shows that the growing hydraulic fracture may exert enough tensile and shear stresses on cemented natural fractures that they may be debonded, opened and/or sheared in advance of hydraulic fracture tip arrival, while under other conditions, natural fractures will be unaffected by the hydraulic fracture. Detailed aperture distributions at the intersection between fracture segments shows the potential for difficulty in proppant transport under complex fracture propagation conditions.