Tight glutenites are typically composed of heterogeneous sandstone and gravel. Due to low or ultra-low permeability, it is difficult to achieve commercial production in tight glutenites without hydraulic fracturing. Efficient exploitation requires an in-depth understanding of the fracturing behavior of these reservoirs. This paper provides a numerical method that integrates the digital image processing (DIP) technique into a numerical code rock failure process analysis (RFPA). This method could consider the glutenite heterogeneities, including intrarock and interrock heterogeneities, and the practicability is verified through two numerical tests. Two-dimensional (2D) simulations show hydraulic fractures (HFs) can penetrate or deflect to propagate along the gravels, depending on the magnitude of stress anisotropy and gravel strength. Three-dimensional (3D) simulations with the consideration of gravel distribution orientation, gravel size and axial ratio show HFs could propagate past the gravel with no deflection, forming a bypass fracture that is not easy to observe in common laboratory experiments. HFs could also deflect to propagate along the gravels. The impacts of the gravel distribution orientation, gravel size and axial ratio are discussed in detail. The main propagation modes of HFs intersecting the gravels are summarized as: (1) penetrating directly; (2) deflecting to propagate along the gravels to form distorted HFs; (3) propagating to bypass the gravels; (4) a combination of (1) and (2), or (2) and (3).