Deep coalbed methane (CBM) reservoirs hold substantial resource potential and play a crucial role in China’s unconventional natural gas development. However, the vertical propagation behavior of hydraulic fractures in deep CBM formations remains inadequately understood, posing challenges for optimizing fracturing parameters to control fracture height growth and enhance fracture development within the coal seam. To address this, this study establishes numerical simulation models to investigate hydraulic fracture propagation in directional wells, incorporating three typical lithological combinations representative of deep CBM reservoirs. Through these models, the influence mechanisms of bedding density, stress ratio, rock friction coefficient, and fracturing parameters on vertical fracture propagation and post-fracture productivity were systematically analyzed. The results reveal that the fracture propagation characteristics vary significantly with lithological combinations. Initially, hydraulic fractures penetrated adjacent formations near the wellbore while simultaneously generating branched fractures, leading to the formation of a complex fracture network. As propagation continues, branch fractures exhibited reduced width compared to the primary fracture. Well-developed bedding planes in the roof or floor, combined with lower stress ratios and friction coefficients, effectively constrained vertical fracture growth. Furthermore, optimizing fracturing fluid volume, reducing injection rate, and lowering proppant concentration promoted fracture development within the coal seam, thereby enhancing post-fracture well productivity. These findings provide a theoretical foundation for the optimization of hydraulic fracturing strategies in deep CBM reservoirs, contributing to more effective reservoir stimulation and resource recovery.

Hydraulic fracturing is one of the most important enhanced oil recovery technologies currently used to develop unconventional oil and gas reservoirs. During hydraulic fracture initiation, fluid seeps into the reservoir rocks surrounding the wellbore, inducing rock deformation and changes in the stress field. Analyzing the fluid–solid coupling mechanism around the wellbore is crucial to the construction design of fracturing technologies such as pulse fracturing and supercritical carbon dioxide fracturing. In this study, a new transient fluid–solid coupling model, capable of simulating the pore pressure field and effective stress field around the wellbore, was established based on the Biot consolidation theory combined with the finite difference method. The numerical results are in excellent agreement with the analytical solutions, indicating the reliability of the model and the stability of the computational approach. Using this model, the influence of seepage parameters and reservoir properties on the fluid–solid coupling around the open-hole wellbore was investigated. The simulation results demonstrate that, during wellbore pressurization, significant changes occur in the pore pressure field and effective stress field near the wellbore. The fluid–solid coupling effect around the wellbore returns to its initial state when the distance exceeds four times the radius away from the wellbore. As the fluid viscosity and wellbore pressurization rate decrease, the pore pressure field and effective circumferential stress (ECS) field around the wellbore become stronger. Adjusting the fluid viscosity and wellbore pressurization rate can control the effect of seepage forces on the rock skeleton during wellbore fluid injection. For the same injection conditions, rocks with q higher Young’s modulus and Poisson’s ratio exhibit stronger pore pressure fields and ECS fields near the wellbore. This model furnishes a dependable numerical framework for examining the fluid–solid coupling mechanism surrounding the open-hole wellbore in the initiation phase of hydraulic fractures.