Hydraulic fracturing operation as an effective enhancing coalbed gas production method is widely used in ultra-low permeability coal seam. However, complex geo-stresses and high heterogeneity between natural cleats structure lead to difficulty predicting hydraulic fracture patterns. Fracture evolution behavior for fracturing operation in coal seams requires a better understanding. In this study, a 2D model of hydraulic fracture propagation was built based on the cohesive zone model of finite element method. The effect of orthogonal cleat system, in-situ stress, dig angle and construction parameters on fracture geometries were main investigated. The main conclusions were as follows: (1) According to the interaction types between hydraulic fracture and cleat system, ladder-shaped fracture and H-shaped fracture geometry was summarized. The difference between them was whether there were continuous and small pressure fluctuation stages. (2) When the horizontal stress difference coefficient was lower and larger than 3/12, fracture geometry was prone to present Η shape and ladder shape respectively. Besides, the dimensionless fracture length and the dimensionless fracture extension aspect ratio of fracture were increasing with larger horizontal stress difference coefficient. (3) The favorable condition for fracture extension was that the dig angle was 45°. Hydraulic fracture tended to propagate along face cleats direction. (4) Larger fracture fluid displacement was beneficial to form more balanced hydraulic fracture geometry and promote large extension scale. As fracture fluid viscosity increased, the fracture geometries transformed from ladder shape to H shape.@en