This study presents development and application of a fully coupled two-phase (methane and water) and poromechanics numerical model for the analysis of geomechanical impact on coalbed methane (CBM) production. The model considers changes in two-phase fluid flow properties, i.e., coal porosity, permeability, water retention, and relative permeability curves through changes in cleat fractures induced by effective stress variations and desorptioninduced shrinkage. The coupled simulator is first verified for poromechamics coupling and simulation parameters of a CBM reservoir model are calibrated by history matching against one year of CBM production field data from Shanxi Province, China. Then, the verified simulator and calibrated CBM reservoir model are used for predicting the impact of geomechanics on production rate for twenty years of continuous CBM production. The simulation results show that desorption-induced shrinkage is the dominant process in increasing permeability in the near wellbore region. Away from the wellbore, desorptioninduced shrinkage is weaker and permeability is reduced by pressure depletion and increased effective stress. A sensitivity analysis shows that for coal with a higher sorption strain, a larger initial Young's modulus and smaller Poisson's ratio promote the enhancement of permeability as well as the production rate. Moreover, the conceptual model of the cleat system, whether dominated by vertical cleats with permeability correlated to horizontal stress or with permeability correlated to mean stress can have a significant impact on the predicted production rate. Overall, the study clearly demonstrates and confirms the critical importance of considering geomechanics for an accurate prediction of CBM production. (1) A fully coupled two-phase flow and poromechanics model for methane recovery (2) Simulation parameters are calibrated by history matching against field data (3) The geomechanics behaviors significantly affect the prediction of CBM production