After hydraulic fracturing, some treatments intended for production enhancement fail to yield predetermined effects. The main reason is the insufficient research about the fracture propagation mechanism. There is compelling evidence that I-shaped fracture, two horizontal fractures at the junction of coalbed and cover/bottom layer, and one vertical fracture in the coalbed have formed in part of the coalbed after hydraulic fracturing. Therefore, this paper aims at I-shaped fracture propagation simulation. A novel propagation model is derived on the basis of a three-dimensional (3D) model, and the coupling conditions of vertical fracture and horizontal fractures are established based on the flow rate distribution and the bottom-hole pressure equality, respectively. Moreover, an associated PDA (pressure decline analysis of post-fracturing) model is established. Both models complement with each other and work together to guide fracturing treatment. Finally, a field case is studied to show that the proposed models can effectively investigate and simulate fracture initiation/propagation and pressure decline.