Abstract: In response to the key issues of unclear fracture propagation patterns and production dynamics of deep coalbed methane (CBM) horizontal wells with staged fracturing in the Shenfu block, this study aims to reveal the fracture propagation and gas production behavior under multi-field coupling effects, providing theoretical support for the efficient development of deep CBM. Based on geological data, well logging, fracturing, and production data from the Shenfu block, numerical simulation methods are employed, using the Benxi Formation 8+9# coal as the research object. A gas-water two-phase flow fracturing-production integrated numerical simulation method, considering adsorption/desorption effects in the flow-solid-damage (HMD) model, is proposed. This study investigates the control of fracturing section design parameters (cluster spacing) and geological parameters (stress difference, elastic modulus, and Poisson’s ratio) on fracture propagation and evaluates the effects of reservoir fracturing on CBM development and optimization of the fracturing scheme. The numerical simulation results show that the HMD model is effectively applied to the study of fracture propagation and CBM development effects in deep CBM horizontal wells with staged fracturing. The ranking of geological and engineering parameters in terms of their control over fractures, from most to least influential, is as follows: stress difference > Poisson’s ratio > elastic modulus > cluster spacing. After fracturing, the coal seam permeability and cumulative gas production increased by 167 times and 5.18 times, respectively, compared to before fracturing. Production simulation results indicate that the highest production is achieved when the fracture cluster spacing is between 20 and 22.5 meters. The study reveals the fracture propagation patterns under multi-field coupling effects in deep CBM horizontal wells, providing valuable insights for optimizing deep coalbed fracturing schemes and evaluating development effects.