Recent examples of hydraulic-fracture diagnostic data suggest that complex, multistranded hydraulic-fracture geometry is a common occurrence. This reality is in stark contrast to the industry-standard design models based on the assumption of symmetric, planar, biwing 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 the latter 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 show the potential for difficulty in proppant transport under complex fracture-propagation conditions.